潜水完整干扰井非稳定流抽水试验计算水文地质参数.doc

潜水完整干扰井非稳定流抽水试验直线图解法(水文地质参数计算) (发表于《内蒙古水利》( 2011 年第 1期) ,荣获“ 2010 — 2011 年度科学发展与构建和谐社会理论实践成果”一等奖, 颁奖单位: 四川西部文献编译研究中心、《环球人文地理》( 理论版)编辑部、中国西部科技杂志社) 内容摘要: 对于潜水完整干扰井非稳定流抽水试验计算水文地质参数, 可应用类似的承压-- 无压水定解问题的太斯解析解公式, 依据“势能叠加原理”求出干扰井解析式, 再进行承压水降深 S 和潜水含水层天然厚度 H 与抽水计算时刻观测井中水柱高 h 的平方差(H 2 -h 2 )近似代换之后,得到相应的潜水干扰井太斯解析解公式。然后,再进行雅可布式简化, 得到便于图解计算的潜水干扰井雅可布公式, 类似于承压水雅可布图解法, 求得水文地质参数 K 的较精确值和μ的可参考值。引言: 在水文地质勘查实践中,为更清楚地了解含水层水动力学特性,能够掌握较切合实际的水文地质参数, 为开采井布局和设计提供较为可靠的出水量及其与水位变化关系的依据性资料, 常布置实施多主孔完整干扰井非稳定流抽水试验。对于潜水含水层, 至近年来的规范规程和教科书, 都未直接给出完整干扰井非稳定流抽水试验计算水文地质参数的具体公式, 在生产实践中, 常需水文地质人员自行推演计算公式, 缺乏表达的统一性和应用的普遍性。本文根据地下水动力学教程提示, 对潜水完整干扰井非稳定流抽水试验计算水文地质参数方法和公式进行论述。一、水文地质模型 1 、含水层在平原区或其它潜水含水层分布较广、连续发育的地区, 含水层厚度分布较稳定, 岩性较均匀,底板倾斜度较小,适用于地下水动力学方法计算抽水试验水文地质参数。 2 、边界条件边界距抽水井距离大于 500m ,可按无限边界处理。图1潜水含水层完整干扰井非稳定流抽水水文地质模型示意图二、单井数学模型及其解析解对于分布连续广泛,厚度相差不大,岩性较均匀的承压—无压含水层中完整井井流,常用如下数学模型描述: 其中: t—抽水时间( T); r—计算水位点与抽水井距离( L); Q—井孔抽水量( L 3T -1)。求解上述定解问题得太斯解析解: 式中: S—任意点观测孔水位降深; Q—单井出水量; T—含水层导水系数; —承压- 无压井函数; —抽水孔到水位观测点距离(L); t—抽水延续时间(T); —含水层释水系数或给水度。三、干扰井基本公式推导及图解计算公式设有 n 个潜水完整井,同时开泵进行非稳定流干扰抽水,各抽水井距观测孔的距离为 r i (任意定值) ,各抽水井出水量为 Q i (任意定值), 抽水延续能满足计算所需要时间 t ,观测孔处水位降深为 S ,根据势能叠加原理,有: 式中: S i—单井抽水量为 Q i时, 任一井单独抽水时,在观测孔处形成的水位将深值。用太斯解析式可表示为: 式中: ,, 其余含义同前。对于潜水,可以对 S 作如下等量代换: 令抽水试验中含水层平均厚度,H 为含水层抽水前厚度,h 为抽水延续 t 时间时观测孔中水柱高。由于,则有: 。用替换 S ,潜水完整干扰井非稳定流抽水定解问题的太斯解析式为: 式中: K—含水层渗透系数。为便于计算,对上式进行雅可布式简化: 当抽水时间足够长到各