Theis井函数计算方法及井模型参数优化计算研究
抽水试验规范方法及计算公式
观测孔的布置应符合下列要求： (1) 为了计算水文地质参数，在抽水孔的一侧宜垂直地下水的流向布置 2~3 个观测孔。 (2) 为了测定含水层不同方向的非均质性或确定抽水影响半径，可以根据含水层的不同情 况，以抽水孔为中心布置 1~4 条观测线；如有两条观测线，一条垂直地下水流向，另一条宜平行 地下水流向。 (3) 群孔干扰抽水试验和试验性开采抽水试验应在抽水孔组中心布置一个观测孔；为查明相 邻已采水源地的影响，应在连接两个开采中心方向布置观测孔。为确定水位下降漏斗形态和补给 (或隔水)边界，应在边界和外围一定范围内布设一定数量的观测孔。 (4) 多孔抽水孔组的第一个观测孔应尽量避开三维流的影响，相邻两观测孔的水位下降值相 差不小于 0.1 m，最远观测孔的下降值不宜小于 0.2 m，各观测孔应在对数数轴上呈均匀分布。 (5) 在半承压水含水层进行抽水试验时，宜在观测孔附近覆盖层(半透水层或弱含水层)中布 置副观测孔。 (6) 在进行试验性开采抽水试验时，应在水位下降漏斗范围内的重要建筑物附近增设工程地 质、环境地质观测点。
(3) 水温、气温宜 2～4 h 观测一次，读数应准确到 0.5 ℃，观测时间应与水位观测时间相对 应。
4.3.4 恢复水位观测要求
停泵后应立即观测恢复水位，观测时间间隔与抽水试验要求基本相同。若连续 3 h 水位不变， 或水位呈单向变化，连续 4 h 内每小时水位变化不超过 1 cm，或者水位升降与自然水位变化相一 致时，即可停止观测。
注意：①要消除区域水位下降值；②在基岩地区要消除固体潮的影响；③傍河抽水要消除河 水位变化对抽水孔水位变化的影响。
多孔抽水试验、群孔干扰抽水试验和试验性开采抽水试验均应编写试验小结，其内容包括： 试验目的、要求、方法、获得的主要成果及其质量评述和结论。
中国石油大学(北京)现代试井分析-第一章 绪论
现 代 试 井 分 析014一、试井概念¾试井是对油、气、水井进行测试和分析的总称。
在不同工作制度下测量井底压力和温度等信号的工艺。
测试内容包括：产量、压力、温度、取样等。
分析（试井解释）：应用渗流力学理论，分析测试数据，反求油层和井的动态参数。
¾试井是一种以渗流力学为基础，以各种测试仪表为手段，通过对油井、气井或水井生产动态的测试来研究和确定油、气、水层和测试井的生产能力、物性参数、生产动态，判断测试井附近的边界情况，以及油、气、水层之间的连通关系的方法。
举例：不稳定试井压力和产量对应关系图二、试井的分类就研究的目的来说⎧⎧⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎩⎪⎧⎧⎨⎨⎪⎪⎩⎪⎪⎪⎪⎨⎪⎪⎧⎪⎪⎨⎪⎪⎩⎩⎩系统试井等时试井产能试井修正等时试井一点法试井试井压力降落试井单井不稳定试井压力恢复试井不稳定试井干扰试井多井不稳定试井脉冲试井按地层类型分类均质油藏试井非均质油藏试井双孔介质油藏试井双渗介质油藏试井复合油藏油藏试井¾按井类别分类，可分为垂直井、水平井、压裂井、定向井和分支井等试井方法。
¾按流动形态分类，可分为线性流、非线性流的试井。
二、试井的分类常规试井分析按分析方法分现代试井分析数值试井分析压降试井分析压恢试井分析变产量叠加试井分析典型图版手动拟合分析典型图版自动拟合分析针对油气藏和油气井研究的严密的测试设计；应用高精度的仪器设备进行现场测试；压力计精度， 分辨率，在井下高温高压条件下连续记录、存储压力数据量测试过程中要求产油气井配合测试进程反复的开关井，准确计量产气量，并处理好产出的气体；以复杂油气藏为背景的渗流力学理论和方法的研究；以解数理方程中的反问题为基础的试井解释方法及软件；结合地质、物探、测井、油藏及工艺措施的油藏动静态描述。
四、试井的作用2014姚约东2014试井研究贯穿于油气田勘探开发全过程2014四、试井的作用试井的作用总结为以下几点：（1）估算测试井的井底污染情况，判断是否需要采取增产措施（如酸化、 压裂），分析增产措施的效果；（2）估算测试井的地层参数、产能；（3）平均地层压力计算、压力分布；（4）判断和预测油气藏类型，均质、非均质油气藏，边底水等；（5）判断和预测油气藏范围，河道油藏，断层距离，透镜体，油（气）层边界，估算控制储量；（6）判断和评价断层的性质，包括密封性等；（7）判断井间连通性；（8）描述井筒周围油藏特性，包括渗透率、孔隙度、厚度、饱和度分布等。
地下水动力学期末考试复习题(井附近地下水流动)
承压水井的Dupuit（裘布依）公式
（1）假设条件（适用条件）
1）水井布置于均质、各向同性、水平分布、等厚的圆形岛屿状承压含水层的中心，岛屿半径为R，岛屿周围自含水层底面起算的水头H0保持不变；——Dupuit模型（圆岛模型）
21、什么条件下会产生承压-无压井流？推到出承压-无压井流公式？
承压水井中大降深抽水时，如果井中水位低于含水层顶板，井附近含水层中水位也将低于含水层顶板而呈现为无压水流，此时就变为承压—潜水井（承压—无压水井）。
承压—潜水井公式：
22、什么是影响半径？
R—影响半径，即从抽水井开始到实际观测（或可忽略）不到水位降深处的径向距离（Thiem的影响半径的定义），m。
2)计算步骤(以降深-时间距离配线法为例)
在双对数坐标纸上作标准曲线W(u)-1/u；
根据实际观测资料，在另一张同模数透明双对数纸上作s－t/r2实际曲线；
将实际曲线叠放于标准曲线上，保持对应坐标轴平行，平移曲线，直至二曲线最大限度地重合；
任取一匹配点(在曲线上或曲线外均可)，记下匹配点对应坐标[W(u)]、[1/u]、[s]、[t/r2]，代入Theis公式计算参数：
2）抽水前含水层水位面水平，水头为H0；
3）抽水过程中地下水运动符合Darcy定律。
数学模型:
数学模型的解——Dupuit公式
采用分离变量法求解，在rw至R区间上进行积分，得到方程的通解，再利用边界条件确定通解中的积分常数，便得上述数学模型的解：
或
公式符号含义：
sw—井中水位降深，m；
Q—抽水井流量，m3/d；
Theis解析模型的数值解及其在MatLab中的实现
望
摘要 ： 了如何应用数值积分方法求解含水层 中地下水非稳定流抽水试验的解析模型 以获 得含水层 的水文 地质 讨论 参数 ， 以 Ｔ ｅ 公 式为例讨论求解原理 。引入 了 Ｍ ｔａ 并 ｈｉ ｓ ａＬｂ的数值积 分工具和 优化工具 ， 轻松 实现 了对 水文地质参 数 的求解 ， 此举既保证 了解 的可靠性和唯一性 ， 也大 大提高 了方法 的可操作性 。最后 以一个具 体 的承压水非稳定 流抽水试 验为例验证 了这个方法 ， 所得 的结果 与应 用配线 法得 到 的结 果基本 一致 ， 而肯定 了方法 的应用 价值 。 从 研究 结果尽管是针对 Ｔ ｅ 模型进行 的 ， ｈｉ ｓ 但求解原理 同样适用于其他各种模型。
１ 研究思路
通过 对含水 层进 行地下 水非 稳定 抽水试 验 以求 解水 文地 质参数 是一 种水文 地质基 础研 究 的通用 方
来进 行 ， 例如矩 形法 、 梯形 法 、 物线 法等 ， 抛 但积 分 区
间无 法 事先 给定 。
②积分区间 ［ ，＋∞］ ｕ 的确定 ， 根据 （ ） 可 ２ 式，
数 值积分 是 通 过 在 有 限 个 采 样 点 上 计 算 ｆ ） （ 的值来 逼 近 ） 区 问 ［ ， ］ 的定 积分 ， 据 离 在 ａｂ内 根 散 区间 内引入 的插值 函数 的不 同 ， 分 为 （ 可 组合 ） 梯 形 公式 、 组合 辛 普 生 公 式 （ 又称 抛 物 线 公 式 ） 布 尔 、 公 式 、 斯一 勒让 德 公式 等 ， 比较 而 言 ， 普 生 公 高 相 辛 式 应用较 为广 泛且 只需对 少量 的函数求 解就 能得 到 准 确 的答 案 。文献 ［ ］ 组 合 辛普 生 公 式 的讨 论 简 ４对
关键 词 ： ｔａ ； Ｍａ ｂ 数值 积分 ；ｈｉ解析模型 ； Ｌ Ｔｅ ｓ 水文地质参数 ； 优化
基于抽水试验的渗透系数计算方法对比与研究
2017年27期科技创新与应用Technology Innovation and Application方法创新基于抽水试验的渗透系数计算方法对比与研究朱保坤，苟联盟，柴永进，张建涛(中设设计集团股份有限公司，江苏南京210014)摘要：以无锡梅梁湖湖区抽水试验为例，对几种常用渗透系数计算方法进行初步的比较和分析。
传统的稳定流公式法计算方法简单，但需对结构进行简化处理，计算结果误差较大；非稳定流配线法计算严谨，准确性高，但通常需借助计算机程序帮助，过程较复杂，对试验数据及 其准确性要求较高。
关键词：抽水试验;稳定流；非稳定流中图分类号：P 641.73 文献标识码：A 文章编号院2095-2945 (2017)27-0071-021概述根据现场抽水试验资料确定含水层渗透系数是查明水文地 质工作中经常采用的有效方法之一 ™。
用抽水试验数据分析水文 ±地质参数的方法已经比较成熟，早在1863年J .Dupuit 就提出了 完整井稳定流的Dupuit 公式，1970年W .C .Walton 提出了承压含 水层完整井有越流、无边界影响时的井流公式，此外还有Theis 公式、Copper-Jacob 法、Boulton 模型及 HantushJacob 法等\ 不 同的计算方法各有优缺点及适用条件，并且存在计算误差，本文 采用稳定流公式法、非稳定流配线法及水位恢复法计算渗透系 数。
2抽水试验简介本次抽水试验项目位于无锡梅梁湖湖区。
试验设置1个抽 水孔，2个观测孔，2个观测孔在抽水孔的一侧，呈直线排列，距离 抽水孔分别为5m 、15m 。
承压含水层厚度约7m ，含水层土性为稍 密状粉土，上下均有厚层隔水层。
桔水孔 观測孔1 观*孔2图1抽水试验剖面简图40 000-*-S Z X 301»5 <m>140030 0CQ警臓8.00 ^10 000\6 00 g 4 W 2008S S 88S S S 888S 8S 888S 8885S 88S 8S S 888S S 88888«爾图2抽水孔的S -t 曲线及Q -t 曲线3渗透系数计算 3.1 稳定流公式法由于场地条件限制，本次试验无法观测抽水孔井外地下水 位，无法消除水跃值对参数计算的影响，故采用两个观测孔观测4000030 000MHASZK302AS-1lll<ESatatt-m -Q (m 3/d >-*-SZK302K3K <m >100090oeo|o e o «~ 20 0005-/050 s0 40 S 100000300200.10l l l l i i i l l l l i i i i ia 888a «8R a a 88S R R R 888sn »图3观测孔1的S -t 曲线及Q -t 曲线数据计算，公式如下：K = 0.366Q &n H ^-S ,) 5r ,式中:K 为渗透系数(m /d ) ;Q 为抽水孔涌水量(m 3/d ) ;m 为含 水层厚度(m )^为1号观测孔水位下降值(m );S 2S 2号观测孔 水位下降值(m );ri为1号观测孔到抽水孔距离(m );r 2为2号观测孔到抽水孔距离(m )3.2Theis 配线法在双对数坐标系下分别绘制W (u )-1/u 标准曲线和S -t /r 2实 测曲线，将实测曲线置于标准曲线上，在保持对应坐标轴彼此平 行的条件下相对平移，直至两条曲线重合为止，任取一匹配点，记 录匹配点的对应坐标值:W (u )，代人式计算参数。
第四章 地下水向完整井的非稳定运动
第四章 地下水向完整井的非稳定运动§1、承压含水层中的完整井流一、定流量抽水时的Theis(泰斯公式)假定条件：（1）含水层均质各向同性、等厚，侧向无限延伸，产状水平； （2）抽水前天然状态下水力坡度为零； （3）完整井定流量抽水，井径无限小； （4）含水层中水流服从Darcy 定律；（5）水头下降引起的地下水从贮存量中的释放是瞬时完成的。
教学模型的建立和求解：抽水形成以井轴为对称轴的下降漏斗，将坐标原点放在含水层底板抽水井的井轴处，井轴为z 轴，建立坐标系。
则以降深表示的微分方程为：教学模型为： 其解为：其中：s —抽水影响范围内，任一点任一时刻的水位降深；Q —抽水井的流量； T —导水系数；t —自抽水开始导计算时刻的时间； r —计算点到抽水井的距离； u —含水层的贮水系数。
利用上述W(u)和u 的关系制定P 93表，W(u)可查表得。
首先由 计算出u 値，然后查表得相应的W(u)，再求r 处的s 値。
tsT urs r rs ∂∂=∂∂+∂∂*221⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧-=∂∂>=∂∂=∞∞<<=∞<<>∂∂=∂∂+∂∂→∞→T r s r t r s t s r r s r t t s T u rs r r s r r π2Qlim 00,0),(00)0,(0,010*22()()⎰∞-===uy dy y e u W Tt u r u u W T Q s 44*2πTtu r u 4*2=二、流量变化时的计算公式流量随时间的变化，可分为阶梯变化和连续渐变，相应的流量过程线为台阶状和连续光滑的曲线。
连续光滑的曲线应概括成阶梯状折线。
概化的原则：矩形面积等于曲线与横坐标所围成的面积。
每一个阶梯视为定流量，用Theis 公式计算降深,然后将各降深叠加起来，得流量变化的总降深値。
如图，连续抽水，概括为4个阶梯，若求t 时刻的水位降深，则可分解为四个亚问题，第一个亚问题以Q 1流量抽水从t 0→t ；第二个亚问题以Q 2－Q 1流量抽水，从t 1→t ；第三个亚问题以Q 3－Q 2流量抽水从t 2→t ；第四个亚问题以Q 4－Q 3流量抽水从t 3→t 。
渤海Q油田提液井生产规律分析
石油地质与工程2021年11月PETROLEUM GEOLOGY AND ENGINEERING 第35卷第6期文章编号：1673–8217（2021）06–0114–05渤海Q油田提液井生产规律分析张锋利1，高智梁1，李权1，欧阳雨薇1，周日1，张利健2（1．中海石油（中国）有限公司秦皇岛32–6作业公司，天津300450；2．中海油能源发展股份有限公司工程技术分公司，天津300450）摘要：为改善渤海Q油田提液生产效果，对油田实施大泵提液井生产动态参数进行分析，总结出油井提液前后含水率变化、后期递减率变化和含水上升速度与含水阶段、能量补充类型、储层厚度和提液幅度有一定的相关性。
据此确定Q油田提液井筛选及优化方向为：优选能量补充类型为底水油藏、储层厚度超10 m、含水率达到95%的井实施提液。
从理论上对提液井生产的合理性进行分析探讨，并对两口井的提液情况进行了验证，为油田下一步提液井筛选、提液时机优化等提供指导。
关键词：大泵提液；因素分析；提液幅度；提液时机中图分类号： TE2 文献标识码：AAnalysis on production law of liquid extraction wells in Bohai Q oilfieldZHANG Fengli1，GAO Zhiliang1，LI Quan1，OUYANG Yuwei1，ZHOU Ri1，ZHANG Lijian2（1. QHD32-6 Operating Company of CNOOC（China) Co.，Ltd.，Tianjin 300450，China;2. Engineering Technology Company of CNOOC Energy Development Co., Ltd.,Tianjin 300450，China)Abstract: In order to improve the production effect of oilfield liquid extraction, according to the production performance parameters of large pumping liquid extraction wells in the oilfield, it is analyzed and summarized that the change of water cut before and after oil well liquid extraction, the change of decline rate in the later stage and the rising rate of water cut are related to the water cut stage, energy supplement type, reservoir thickness and liquid extraction amplitude. Therefore, the selection and optimization direction of liquid extraction wells in Q oilfield is proposed as follows: the preferred energy supplement type is bottom water reservoir, the reservoir thickness exceeds 10 m, and the water cut reaches 95%. To analyze the rationality of the production law of the liquid extraction wells, the theoretical mechanism has been discussed, and then the liquid extraction conditions of two wells are verified, which provides a good guidance for the selection of liquid extraction wells and the optimization of liquid extraction time in the next step.Key words: large pumping liquid extraction; factor analysis; amplitude of liquid lifting; liquid extraction time 渤海Q油田油藏埋深较浅，为疏松砂岩高孔高渗储层，产液能力高，目前油田平均单井日产液量可达600 m3以上，综合含水接近95%，处于特高含水期。
4-1承压含水层中的完整井流
1 s 当 r一 定 时 ， 和e u 起 着 增 、 减 速 度 两 个 向 方的 作 用 ， 所 以 不 是 t t t的 单 调 函 数 ， s t曲 线 不 能 沿 着 同 一 斜 变 率化 ， 存 在 着 拐 点 。
2.拐 点 位 置 ： 说 明 st位 置 存 在 拐 点 ， v 2 s Q 1 u r 2 μ* 2 2e ( 1 ) 0 t t 4πT t 4Tt r 2 μ* 当u 1， 即 1， 出 现 拐 点 即 最 大 降 速 点 的。 4Tt r
6
（4）确 定 系 数 C dH H Q r ，（当u 0时） du r 2T u 0 dH 上 式 变 为 ： 2u 2Ce u，两 边 同 时 求 极 限 ， du dH Q u （据 边 界 条 件 ） 2 Ce 2 C ； 2 u lim lim du 2T u 0 u 0 Q Q 方 程 左 边 ＝ 右 边 ， 即 2C 所以C , 2T 4T dH Q e u 把C代入（3）得 即 du 4T u H0 Q e u 分离变量积分 ： H dH 4T u u du，（当H (u ) H 0时，u ） Q e u Q e u H0 H du，得H H 0 du u u 4T u 4T u 即 为 定 解 问 题 （ 1） 的 解 ， 也 是 著 名 的 泰 斯式 公。 lim 2u
r *
19
（2）水位下降速度的变化规律
u s Q e Q 1 v du e u t u 4T u t 4T t 1.当 t一 定 时 ， r增 大 v减 小 ； r增 大 ， 降 速 越小， r ，v 0；
u r 2 4Tt
aquifertest软件与人工配线法计算矿区水文地质参数
水文地质参数是矿山预测矿坑涌水量的重要基 础数据，参数计算的精度直接影响了矿山涌水量预测 的准确性，影响矿山的防排水设计和防治水方案的选 择。在实际生产中，矿山为了获得比较准确的水文地 质参数往往利用钻孔稳定流和非稳定流抽水试验获 得数据，然后稳定流采用公式法计算，非稳定流采用 直线图解法和配线法计算求解。稳定流试验简单，计 算方便，但是获取参数较少，故对于水文地质条件中 等及以上矿山，一般采用非稳定流试验获取矿山含水 层的水文地质参数。但是非稳定流试验后期数据计 算过程比较繁琐，导致人工计算参数时效率低，且结 果往往因人而异，难以有统一的标准，最终影响矿山 涌水量预测结果的准确性。近年来随着计算机技术 的进步，水文地质参数计算技术也得到了飞速发展， 加拿大滑铁卢公司开发的 ＡｑｕｉｆｅｒＴｅｓｔ水文计算软件 针对非稳定流抽水试验可以采用多种方法求得含水 层水文地质参数，大大提高了数据处理效率和计算精 度，已经逐步成为矿山水文地质参数计算的主要应用 工具 。 ［１３］
４
０．００６
３５
６
０．００６
４０
８
０．００７
４５
１０
０．００８
５０
１２
０．００９５
６０
１４
０．０１９
７０
１６
０．０２６ຫໍສະໝຸດ １００２００．０４３
１３０
２５
０．０７４５
１６０
０．１０５ ０．１４３ ０．１８ ０．２２３ ０．２４８ ０．２９１ ０．３４７ ０．４５８ ０．５４５ ０．６１３
１９０ ２２０ ２５０ ２８０ ３１０ ３４０ ３７０ ４００ ４３０ ４６０
15第十五章 地下水向井的非稳定运动
一、承压含水层中的完整井流
（一）、定流量抽水时的Theis(泰斯公式)
一、承压含水层中的完整井流
（一）、定流量抽水时的Theis(泰斯公式)
2数学模型的建立 在上述假设条件下，抽水后将形成以井轴为对称轴的下降漏斗，如图 后示，建立坐标轴，以井底为中心为原点坐标，井底板向右为r轴向。 假设抽水由t时刻至t+ t时刻时，井中水位变化为H 据水均衡原理（Qr+dr－Qr) · t=V入－V出＝ V贮……..(1) 又据高数知识：f(x+ x)=f(x)+f’(x) · x 所以，Qr+dr=Qr+Qr’ · ……………….. (2) dr
r 2* 4Tt
r 2 * 1 4T t 2
r 2 * 4Tt 1 0
r ti 4T
2*
所以，当r一定时，t＜ti时，断面的水头降速由小逐渐增大；当 t=ti(u=1)时，达到最大；当t＞ti时，下降速度又由大变小。 不同的断面拐点出现的时间ti不同。 拐点的降深：
一、承压含水层中的完整井流
（一）、定流量抽水时的Theis(泰斯公式) 2数学模型的建立
1 r H H 2 r r r T t
2 *
2 H 1 H * H t 0, 0 r r 2 r r T t 0r H ( r ,0) H 0 H H0 t0 H ( , t ) H 0 , r r H Q lim r r0 r 2T
)......... ......( 4)
T KM
由（4）求导得：
H Qr 2T r r r
一、承压含水层中的完整井流
一种求解承压含水层水文地质参数的新配线法
一种求解承压含水层水文地质参数的新配线法周志芳;王萍;李雅冰;吴建;夏文俊;王峻;方忠强【摘要】In order to make full use of the continues pumping test data to estimate the hydraulic conductivity (K), the transmissivity (T) and the storage coefficient (μ*) of different head depression cones in a confined aquifer, a dimensionless analytical solution of drawdown transients with variable pumping rates was derived based on the Theis formula and the principle of superposition.Furthermore, a new type curve method was proposed.The in-situ pumping test was conducted to obtain the drawdown record perturbed by stepwise rates.Then, the drawdown curve was used to match the corresponding type curve.After that, this method selects one match point and records its coordinate values on both the type curve and the real drawdown curve.By substituting these recorded values into related formulas, the hydraulic parameters of the confined aquifer can be determined.This method was applied to a field pumping test conducted at a southern highway in Wuxi, Jiangsu Province, China.The type curve matching process is implemented in stages.Field test results reveal that the hydraulic conductivity of the confined aquifer with the 2-3 silt layer are 5.12×10-4cm/s, 2.54×10-3cm/s and 2.83×10-3cm/s respectively, which are corresponding to three-stage stepwise pumping rates.%为了能充分利用抽水试验规范所规定的连续性抽水试验数据资料,原位确定承压含水层不同水头降落漏斗范围的渗透系数(K)、导水系数(T)和储水系数(μ*),基于Theis井流公式及叠加原理导出一种流量变化时无量纲形式的承压含水层降深随时间变化的解析解,提出一种含有参变量的一组无量纲降深-时间标准曲线及相应的配线方法.将变流量抽水试验得到的降深随时间变化的实测曲线与标准曲线进行配线,选择任一匹配点记录其在标准曲线、实测曲线中的坐标值和配线的参变量值,将这些记录值代入相关公式,确定对应某一流量时承压含水层的水文地质参数.以苏锡常南部高速无锡境内的场地抽水试验为例进行实例验证,运用配线法分阶段完成,最终求得该场地2-3粉土层承压含水层在变流量下的渗透系数分别为5.12×10-4cm/s、2.54×10-3cm/s、2.83×10-3cm/s.【期刊名称】《河海大学学报（自然科学版）》【年(卷),期】2019(047)001【总页数】6页(P7-12)【关键词】连续性抽水试验;承压含水层;Theis井流;叠加原理;参数解析解;配线法;参变量;水文地质参数;苏锡常南部高速公路;无锡境内场地【作者】周志芳;王萍;李雅冰;吴建;夏文俊;王峻;方忠强【作者单位】河海大学地球科学与工程学院,江苏南京 210098;河海大学地球科学与工程学院,江苏南京 210098;河海大学地球科学与工程学院,江苏南京 210098;河海大学地球科学与工程学院,江苏南京 210098;江苏省交通工程建设局,江苏南京 210004;江苏省交通工程建设局,江苏南京 210004;中设设计集团股份有限公司,江苏南京 210014【正文语种】中文【中图分类】P64;TU443水文地质参数对于地下工程抗浮设计、基础工程围护结构设计与施工、工程安全运行评价等有着至关重要的作用[1-2]。
底水油藏中水平井的合理位置优化研究
．
文献标 识码 ： Ａ
关键 词 ：底水 油藏 ； 平井 ； 向位 置 ；多 目标 优化 水 垂
在矿 场实践 中 ， 底水 油藏水 平井 的垂 向位置 是影 响水 平 井 开发效 果 的 重要 因 素 。水 平 井 离油 水 界 面越 远， 其临界 产量 越高 ， 反之 越靠 近油藏 顶部 ， 底水 油藏 渗流动 力越 小而 阻力增 大 ， 而 降低水平 井 的产能 。因 进
产量 随 渗透 率各 向异性 比的增 加而 降低 。因此 ， 油藏 垂 向渗 透率 的提 高将 有 利于 底水 油藏 水平 井 产 能 的提
高。
３ 底 水 油 藏 水 平 井 临 界 产 量
考虑 油藏 各 向异性 的影 响 ， 水 油藏 水平 井 的临界 产量 为 底
㈤
Ｉ ｎ
ｎ ｎ Ｊ ｌ Ｉ ＿
无 因次水平井垂 向位置ｈ 。
／ ：— — ０ ２ ； … Ｋ ５ １ ０； … ２ ５ … ８ ０ ．；
４ 底 水 油 藏 水 平 井 见 水 时 间
考 虑 油藏 各 向异性 的影 响 ， 水 油 藏水 平 井 见 底
水 时 间 的计 算公 式ｌ 为 ６
一
图 ３ 水 平 井 临界 产 量 与 水 平 井 垂 向位 置 关 系 曲线
考虑水 平井 位置 对水平 井产 量 、 临界产 量和 见水 时间 的影 响 ， 合 多 目标最 优 化 理论 ， 结 建 立 了以底水 油藏 水平井 位 置为优 化变 量 ， 以水平井 产量 、 临界 产量 和见水 时 间为 目标 函数
的多 目标 优化模 型 ， 出了水平井 合 理位置 的优 化方 法 。结 果表 明 ， 水油藏 中水 平井 的 提 底
多 目标 优化是 指在 系统 中要考 虑 的 目标 不止一 个 ， 而且 各 目标 之 间又 互相 冲突 、 消彼 长 的情 况 下 ， 此 找 出一个 能 同时满 足所有优 化 目标 的解 。但在 多 目标 优化 问题 中 ， 多个 目标 同时达 到最 优是不 可 能的 ， 只 使 而 能从 所 有可 能的解 集 （ ａｅｏ解集 ） Ｐ ｒｔ 中找 到最 合理 、 靠 的解 ， 各 目标 尽 可 能达 到 最优 。多 目标优 化 数 学 可 使
钱塘江涌潮对承压水位势的波动响应机理探究
1引言在深基坑工程中,造成基坑工程事故的最直接原因就是地下水,特别是承压水,承压水基坑容易产生基坑突涌、坑底凸起等破坏现象。
因此,承压水基坑应采取一定措施防止基坑在开挖深度过深或压强水头过高时发生失稳。
对于承压水失稳的研究,目前有压力平衡理论、土体剪切破坏理论、土体扰曲破坏理论等,关于基坑突涌稳定性的判断分析方法,不管是哪种理论,都用到水头压力。
影响承压水水头压力变化的因素有多种,如大面积抽水、潮汐能等,地表水体也有潮汐现象,例如,海洋潮汐及江河潮汐等。
钱塘江潮就是属于潮汐能的一种。
2工程概况拟建杭州地铁9号线沿线场地属钱塘江河岸及两岸的钱塘江河口相冲海积沉积平原;杭州段钱塘江河床由于近海,为河口地形,属感潮河,受下游潮汐及上游流水冲刷作用。
根据详勘勘探孔揭露的地层结构、岩性特征、埋藏条件及物理力学性质,沿线岩土层按其成因分类如表1所示。
表1地基土层位划分表岩土名称层厚/m 地下水类型1层填土3层砂性土6层淤泥质粉质黏土9层黏性土12-1层粉砂12-4层圆砾20层泥质粉砂岩0.00~2.002.00~19.0019.00~25.0025.00~31.0031.00~38.0038.00~50.0068.00潜水含水层潜水含水层隔水层隔水层承压含水层承压含水层基岩含水层3钱塘江涌潮及波浪特性3.1钱塘江涌潮特性钱塘江涌潮除地理地势因素外,潮汐本身的变化以及钱塘江径流强弱也助长涌潮的生成,钱塘江河口段的下游,潮位、流速过程线均接近简谐曲线。
即当钱塘江潮位上升时为涨潮、涨至最高水位为高潮位,落至最低水位时的水位为低潮位,相邻两高潮位或者低潮位的时间间隔为潮周期。
潮位一日一升一降为全日潮,一日两升一降为半日潮,钱塘江涌潮周期【作者简介】陈军（1983~），男，福建漳州人，高级工程师，从事岩土工程勘察设计与研究。
钱塘江涌潮对承压水位势的波动响应机理探究Study on the Fluctuation Response Mechanism of Qiantang River Tidal Bore tothe Pressure Water Level Potential陈军，范大军，庄兴岳（浙江省地矿勘察院有限公司，杭州310063）CHEN Jun,FAN Da-jun,ZHUANG Xing-yue(Institute of Geology and Mineral Exploration Zhejiang,Hangzhou 310063China)【摘要】根据前人关于压力含水层有越流补给时所建立完整井的公式，得出了压力含水层在钱塘江潮汐荷载作用下水势响应的解析式，将不同潮水位、抽水量带入解析式中，对承压水的变化进行曲线绘制；并与杭州地铁9号线及钱江新城片区的承压水观测资料进行分析对比，得出理论与实测两者之间存在的联系，以及钱塘江潮汐对承压水水位的变化规律。
岩土工程渗流：第5章 地下水井流理论
5.1.3 其他井流情 况与相关的约定
向井中注水的情况实际上和抽水的情况一样，在计算与分 析中可以同样利用抽水的相关理论。
本章除特别提到的条件外，一般都用了以下的假定： (1)含水层中水流服从Darcy定律 (2)在水头下降的瞬间水就释放出来（专指潜水浸润面） (3)含水层均质、各向同性、无限延伸 (4)含水层底部水平，承压含水层等厚 (5)抽水前的地下水面是水平 (6)忽略弱透水层的贮水性
17
计算影响半径R的经验公式和经验值
18
19
2）两个观测孔的井流公式
Q 2 KM H0 hw
ln R rw
sw
H0
hw
Q
2 KM
ln R
rw
边界条件并不一定在r=R和r=rw，可以是任何已知点的 水头值（ri，Hi）或降深值（ri，si）。
承压井流量 一般形式为：
Q
2
KM H1
计算机广泛应用 前的时代产物
14
5.2.2 承压井的齐姆（Thiem）公式
1）引用影响半径
影响半径R的取值，理论上没有严格证明，实际上 是困扰人们的一个问题。 该问题分为两类：
无限大区域中R取值； 有界区域R取值
2）两个观测孔的井流公式
由于井壁影响，应用时有误差，工程中常用两个 观测孔的井流公式。
s(x, y,t) H0 (x, y, 0) H (x, y,t)
井中心处降深最大，离井越远降深越小，整个水头下降 区呈漏斗状，称为降落漏斗。
潜水井的降落漏斗在含水层内部扩展，有自由面。
承压水井的降落漏斗没有自由面，是水头的降低区。
如无其他来源，潜水井抽出的水量相当于降落漏斗的含 水层体积的重力疏干，而承压水井的水来自于因降深漏 斗处的水头降低造成含水层的弹性释水。（非稳定井）
抽水实验定渗透率方法总结
用抽水试验确定渗透系数1.抽水试验资料整理试验期间，对原始资料和表格应及时进行整理。
试验结束后，应进行资料分析、整理，提交抽水试验报告。
单孔抽水试验应提交抽水试验综合成果表，其内容包括：水位和流量过程曲线、水位和流量关系曲线、水位和时间(单对数及双对数)关系曲线、恢复水位与时间关系曲线、抽水成果、水质化验成果、水文地质计算成果、施工技术柱状图、钻孔平面位置图等。
并利用单孔抽水试验资料编绘导水系数分区图。
多孔抽水试验尚应提交抽水试验地下水水位下降漏斗平面图、剖面图。
群孔干扰抽水试验和试验性开采抽水试验还应提交抽水孔和观测孔平面位置图(以水文地质图为底图)、勘察区初始水位等水位线图、水位下降漏斗发展趋势图(编制等水位线图系列)、水位下降漏斗剖面图、水位恢复后的等水位线图、观测孔的S－t、S－lg t曲线[注]、各抽水孔单孔流量和孔组总流量过程曲线等。
注意：（1）要消除区域水位下降值；（2）在基岩地区要消除固体潮的影响；3）傍河抽水要消除河水位变化对抽水孔水位变化的影响。
多孔抽水试验、群孔干扰抽水试验和试验性开采抽水试验均应编写试验小结，其内容包括：试验目的、要求、方法、获得的主要成果及其质量评述和结论。
2.稳定流抽水试验求参方法求参方法可以采用Dupuit公式法和Thiem公式法。
(1)只有抽水孔观测资料时的Dupuit公式承压完整井：潜水完整井：式中K——含水层渗透系数(m/d);Q——抽水井流量(m3/d);s w——抽水井中水位降深(m);M——承压含水层厚度(m);R——影响半径(m);H——潜水含水层厚度(m);h——潜水含水层抽水后的厚度(m);r w——抽水井半径(m)。
(2)当有抽水井和观测孔的观测资料时的Dupuit或Thiem公式式中h w——抽水井中水柱高度(m);h1、h2——与抽水井距离为r1和r2处观测孔（井）中水柱高度(m)，分别等于初始水位H0与井中水位降深s之差，h1= H0–s1；h2= H0–s2。
地下水动力学课程中的基础数学运用
地下水动力学课程中的基础数学运用作者：王超月来源：《教育教学论坛》2020年第45期[摘要]在地下水动力学教学过程中，发现不少学生存在数学基础不够扎实的问题，在学习这门课程时，很难与学过的数学知识联系。
据此，通过将教材中用到的数学知识进行系统的梳理，主要包括基本方法、基本理论、数学模型及其求解三部分，期望对学生或从事该课程教学的教师有所帮助。
[关键词]地下水;水文地质;动力学;数学;基本原理[基金项目]河北地质大学2018年度校级教学发展与改革实践项目“地下水动力学课程教学资源库建设”（2018JF06）[作者简介]王超月（1988—），男，博士，河北地质大学水资源与环境学院讲师，主要从事水文地质学理论与教学研究。
[中图分类号] G640 [文献标识码] A [文章编号] 1674-9324（2020）45-0-04 [收稿日期] 2020-06-12地下水动力学是地下水科学与工程、水文与水资源工程、地质工程等专业本科生的一门基础理论课。
其先修课程有高等数学、普通地质学、水文地质学基础、水力学等。
它是地下水科学或水文地质学工作者必要的知识储备，也是地下水数值模拟的理论基础。
课程主要讲述了地下水运动的基本原理以及计算方法。
该课程的特点是用到的数学知识多，理论性强，并且具有鲜明的实际应用背景。
笔者仔细翻阅了《地下水动力学》[1，2]以及《高等数学》[3]等教材，并结合多年授课经验，梳理了课程中主要的数学知识点。
其中也涉及笔者对一些基本概念的理解，有不当之处，恳请读者给予指正。
一、基本方法与高等数学纯数学相比，作為专业课，多了实际背景。
差别首先体现在变量符号的使用上。
数学里通常用y表示因变量，x表示自变量。
而在地下水动力学课程中，各物理量均有实际的物理意义，x，y，z用来表示空间自变量，t表示时间自变量，H，p，v分别表示地下水水头、压力、渗流速度（来源于相应英文单词首字母），它们是时空变量的函数，如H（x，y，z，t）。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
由于 Theis 井函数为指数积分式, 数学上 无法 直
W (u)=
i -ln(u)+ða• i u i=0 4
5
式 中 ：γ =0.5772156649015329， 是 欧拉 常 数 ；N 理论 上 趋近于无穷㊂ 展开 项 N㊂ 国内 绝 大 多数 研究 者 都 是采用的这种方 似 数 值 计算 方法会 因 为 级 数 项 的取 舍 而导 致实际应 5~11 时的 u-W (u)图， 可以很明显地看出当 u 大于 2 时 大,其后果是模型参数最终求解结果出现较大误差㊂
数 W (u)的方法,其替代公式如下： 的误差㊂
2 • (1-G ) ,h = (1-G ) • (G 2-6 • G +12) ɕ G• (2-G ) 3 • G• (2-G )2 • b
如图 2 示,拟合式(6)当 u 大于 2 时 W (u)出现明显 笔者 基 于 最小 二乘原理 , 得出 如 下 井 函 数 W (u)替
⎧ | | | | ⎨ | | | | ⎩
解几乎完全重合： W (u)=
代公式,具体系数见表 2,效果图见图 2,和高精度数值 b 1+b 2 • ub 3+b 4 • ln(u)+b 5 • exp(b 6 • u)
p
系数如表 1㊂
Srivastava 于 1998 年提出的方法如公式(4)示,其
4.156969E-6 -2.624277E-5 4.156478E-6 4.156969E-6 4.156969E-6 4.112166E-6 0.7894027 100.17575
10
9.835525E-11 -52247.83580 4027.53816
程先云 1， 郑凡东 2， 杨
浩 2， 曹大岭 1,杨
勇2
优缺点及对模型最终参数计算的影响， 提出了新的简易计算方法； 以经典实际案例为对象， 对比研究 数优化计算问题， 并可广泛用于其它水文模型参数求解㊂ 关键词:Theis 井模型；Theis 井函数；参数优化；1stOpt 中图分类号:P345 文献标识码: A 应慎用广为采纳的级数展开井函数计算方式；(2)1stOpt 通用优化计算平台可很好处理 Theis 井模型参 文章编号:1000-0852(2015)03-0008-06 即采用该平台进行模型参数优化计算㊂ 了优化软件平台在 Thies 井模型参数计算中的实际应用效果；计算结果表明：(1)Thies 井模型参数估算
2.5 1.5 0.5 -0.5 -1.5 2 1 0 N=5 N=7
(2)
1 • i=0 4 u • exp(u)
ðb•u ðc•u
i=0 i
i
u>1
(4)
实际应用中, 可以根据不同的精度需要选取级数
序号 a b c
法,级数项 N 一般都取 3~5,如文献[1-4]等㊂ 然而该近 用 中 产 生较大 误差 问题 , 如 图 1 为 级 数 项 N 分 别 取 W (u)值 即处于 震荡 发 散 不 稳 定 状态 ,N 值 越 小 误差 越
-0.57722 3.95850 0.2677
0
表1 Srivastava方法系数表 Table1 Coefficient of Srivastava method 0.99999 1 -0.24991 25.63296 18.05902 2 0.05519 3 -0.00976 1 1 4
21.09965
能通过数值计算方法近似求解㊂ 关于井函数 W (u)， 在数 称 为 Theis 井 函 数或井 函 数， 其 数 值 求解 方式国内外 学者已进行了众多的探讨和研究， 除了传统的误差较 主要可归纳总结如下㊂ 大且 实施繁琐 的配 线 法和 Jacob 直 线图 解 法外， 其 它
收稿日期:2014-04-26 基金项目:北京市西郊地区地下水战略储备关键技术研究与示范(Z141100003614060)；水利部公益性行业科研专项经费项目 山洪灾害监测预警 系统标准化研究 (201201058) 作者简介:程先云 (1965-),男， 四川宜宾人， 工学博士， 高级工程师， 主要从事水信息研究㊂ E-mail:xycheng@
⎞ | | | | | | | | | | | | | ⎠
u>4
﹙
其中：G =exp(-γ),γ=0.5772156649015328606 ；
ɿ 26 1 • q ,q = 20 • h= + hɕ u 1+ q 47 1+u •ɿ u 2.3 Srivastava 方法[6]
31
b=
ɿ
﹙
﹚
﹚
(3)
误差 小 于 1E-10 的 高 精 度 数 值 解 , 以 此 结果 作 为参 表 3㊂ 照 目 标 值 检 验 上 述 不同 井 函 数 计算 方法的 精 度 , 见
u 高精度 级数 展开 W (u) 切比雪夫 Srivastava 拟合 方程 式(6) 式(7) Barry N=5 N=30
表3 不同方法井函数计算精度对比 Table3 Accurate comparsion of various methods for well function 10.9357198 10.9357195 10.9357195 10.9357198 10.9357198 10.9357198 9.46056389 10.9357198 1E-5 4.03792958 4.03792924 4.03792924 4.03783522 4.03792958 4.03787773 4.03792515 4.03793067 0.01 0.21938393 0.21958956 0.21938360 0.21927152 0.21938393 0.21938351 0.21938407 0.221 1
第3期
程先 云等：Theis 井 函数计算方法及井 模型 参数优 化计算研究
⎧ | | | | | | | | | | | ⎨ | | | | | | | | | | | ⎩
9 uɤ1
i i
2.1 级数展开
接得出其积分表达式,为便于数值计算,通常将其以无 穷级数形式展开,如下式： • u W (u)=-γ-ln(u)-ð (-1)n n • n! n=1
-i
2.2 Barry 方法 法,概要如下：
[5]
3,i
-i
Barry 于 2000 年提出了一种估计井函数 W (u)的方 G 1-G exp(-u)ln 1+ u ( h + b • u)2 W (u)= -u G +(1-G ) • exp 1-G
式中：p 1㊁p 2㊁p 3㊁q 1㊁q 2 和 q 3 为常数项系数, 具体值可参 考文献[7]㊂ 2.5 拟合方程替代 巩彦文[8]等提出 了 利用 拟 合 方程 来近 似计算 井 函 W (u)=18.109 • u0.256-31.310 • u0.16+13.482 (6)
2
Theis 井 模 型
表示为：
定流量非稳定流承压水完整井 Thies 井模 型可以 S=
ɕ i
其中：W (u)=
显 而 易 见 ， 虽然国内外 还 有不 少 其 它 替 代 计算 方式 ， 但都缺乏较为全面和系统的对比介绍㊂ 第二点是模型 参数估算， 虽然自行编程亦 可实现 ， 但 由 于涉及 到 优 化和数值算法以及编程技巧， 对一般水文工作者而言 有相当难度， 而随着当今计算 软硬 件 的 快速 发 展 ， 利 用成熟优化建模计算平台不失为一可用之道， 这类平 台林林总总㊁ 性能各异， 有的优化能力强但无内置井 函数计算函 数 ， 如 Lingo㊁GAMS 等； 有的 井函 数为内 应用效果已被国内外各行业广大研究者认可㊂ 本研究 1stOpt 却可很好满足两者的要求， 使用上也非常简单， 置 函 数 ， 但 优 化 能力 却 相对 较 弱 ， 如 Matlab 等 ； 而
p 5 p• • p 8) u>1 1 exp(1+p• 2 u 3+p• 4 u +p• 6 ln(u))+p• 7 exp(u
uɤ1
(7)
10
4.3 3.8 3.3 2.8 2.3 1.8 1.3 0.8 0.3 -0.2
水 文
表2 拟合替代公式 (7) 系数表 Table2 The coefficient of regression function 序号 1 2 3 4
契比雪夫有理式逼近可用于高精度计算复杂函数
-ln(u)+
ðp •u ðq •u
i=0 1,i i=0 n1 1,i
n1
i
i
0<uΒιβλιοθήκη 1-1 0 1 2 u 3 4 5
N=9 N=11 N=10 N=8 N=6
• W (u)= exp(-u)
ðp
i=0 n2
⎛ | | | | | | | | | | | | | ⎝
式 中 ：S 为 实 测 降 深 (m)；Q 为 抽 水 流 量 (m3/d)；T 为 承 压含水层导水系数(m2/d)；W (u)为 Theis 井函数；us 为承
ʏ
exp(-x) dx,u= r2 • us x 4 • T • t
Q • W (u) 4 • π • T
(1)
压含水层贮水系数；t 为抽水时间(d)；r 为观测 孔到抽 水井的距离(m)㊂ T 和 us， 难点之一是公式中的井函数 W (u)无解析解， 只 学上也叫指数积分函数 (ExpInt)， 在水文地质中则常 根据抽 水 试 验 资料 和 模 型公式 (1)确定 模 型参数
0.564514611673438
5 6 7 8
Fig.2
图 2 拟合替代公式井函数计算示意图 The regression function for well function