承压水完整井井群互阻计算中流量削减法的一种求解方法
【专利摘要】承压水完整井井群互阻计算中流量削减法的一种求解方法，属于市政工程技术领域。针对以往流量削减法尚不能正确计算干扰抽水的出水量减少系数、进而不能正确计算干扰出水量的问题，本发明建立了干扰抽水的附加降深与单独抽水的附加降深的关系式，并基于此关系式，建立了实用的干扰抽水时各井出水量减少系数与水位降深关系的方程组，由其可计算干扰抽水时各井的出水量减少系数；为提高计算效率进一步建立了干扰抽水出水量减小系数的矩阵解析式，利用此式可方便的计算干扰抽水时各井的出水量减少系数，进而计算各井的干扰出水量。结合一个实施例论述了实施步骤，并进一步验证了本发明方法的正确性。
【专利说明】承压水完整井井群互阻计算中流量削减法的一种求解方法 -、技术领域
[0001] 本发明属于市政工程技术领域。承压水完整井井群互阻计算中流量削减法常用于 水源地的井群规划，其求解方法及其结果直接影响水源地的井群布置、抽水设备的选择等。 二、【背景技术】
[0002] 水源地常由很多眼管井组成，当井距大于或等于两倍的影响半径，各井抽水不产 生干扰。但占地面积大，井群分散，供电线路和井间连络管投资大，且不便管理。因此，一般 井群系统井距小于两倍的影响半径W，运样各井同时抽水时会产生干扰现象，即井群互阻。 井群互阻计算常见问题是:在干扰抽水情况下，当各井水位降深一定时，求干扰出水量;或 反之。目前，井群互阻计算有稳定流干扰井法、非稳定流干扰井法。而稳定流干扰井法又分 为水位削减法和流量削减法。
[0003] 对于干扰抽水，若各井的水位降深与其独立抽水时的水位降深相同，因存在井群 互阻影"^心WAAn心'''S则各井的出水量必定减小，出水量减少系数为[2^6]
[0004] (1)
[0005] 式中：Cii--干扰抽水时第i号井受其他各井抽水共同影响的出水量减少系数，i = l，2，…，n;Qi、Q/i--第i号井单独抽水、干扰抽水时的出水量，m3/d，i = l，2，…，n。
[0006] 若求得了出水量减少系数Cii，贝。
[0007] Q'I = Qi(I-Qi) (2)
[000引式(2)即为流量削减法计算干扰抽水量Q/1。可见，流量削减法关键在于推求出水 量减少系数曰1。
[0009] W两眼承压水完整井互阻抽水为例，分析出水量减少系数，并推广到n眼井。
[0010] 方法一:先明确单独抽水相应的水位降深，在再此降深下分析干扰抽水的降深构 成，进而由该降深相应的干扰出水量、单独抽水的出水量确定出水量减少系数。分析得 出口-4]:
[0011]
(3)
[0012] 式中：ai-2-一干扰抽水时第1号井受第2号井影响的出水量减少系数;t/1-2-一第 2号干扰抽水流量皆2在第1号井产生的附加降深，m;SS-一干扰抽水时第1号井的水位降 深，或简称降深(下同），是干扰抽水流量分i、Q/2综合作用的结果，m。
[OOK]推广至巧扰抽水时的n眼井:第巧井受其他各井抽水共同影响的出水量减少系数 ai[3](i = l，2,…，n)为：
[0014]
(4)
[0015] 式中：ay--干扰抽水时第i号井受第j号井影响的出水量减少系数，i = l，2,…， n，j = l，2，…，n，i声j;t/d-一干扰抽水时第j号井抽水流量皆J在第i号井产生的附加降 深，m，i = 1，2,…，n，j = 1，2,…，n，i声j ; S' i--干扰抽水时第i号井的水位降深，i = I， 2,…，n，是各井干扰抽水流量分1综合作用的结果，也即各井的设计降深，m。
[0016] 对于干扰抽水的任意两眼井：
[0017]
(5)
[0018] 式中各符号含义同前。
[0019] 式(5)表明，任意刚艮井i，j，干扰抽水第j眼井对第i眼井干扰影响的出水量减小 系数等于第j眼井在第i井产生的附加降深与第i眼井的水位降深之比。
[0020] 目前未能利用式(4)解决已知各井的降深^ 1，计算各井的〇1，进而计算干扰出水量 的问题。因为方程组(4)中各个附加降深油应出水量为Q/J，而皆功待求的，故tVj 也未知，且同时抽水时也无法测出t/1-2,因此无法确定〇1。故目前式（4)未能用于实际计 算口-4]。
[0021] 方法二:先明确干扰抽水相应的水位降深构成，在再此降深下分析干扰出水量、单 独抽水的出水量，进而由该降深相应的干扰抽水与单独抽水的出水量确定出水量减少系 数。对于干扰抽水的两眼承压水完整井，文献[2-引得出：
[0022]
(6)
[0023] 式中：ti-2-一第2号井单独抽水时出水量化在第1号井产生的附加降深，m;Si-一 单独抽水时第1号井出水量化相应的水位降深，m。
[0024] 推广到一般：
[0025]
(7)
[0026] 式中：td-一单独抽水时第j号井出水量如在第i号井产生的附加降深，m，i = l， 2，…，n，j = l，2，…，n，i声j;Si--单独抽水时第i号井出水量化相应的水位降深，m，i = 1， 2,…，n;其他符号含义同前。
[0027] 目前采用式(6)、式(7)确定干扰抽水出水量减小系数〇d。然而，分析表明式(6)、 式(7)是不正确的，混淆了干扰抽水与独立抽水时的出水量、水位降深。现结合附图进行分 析，说明如下。
[0028] 若两井分别WQ/ i，Q/ 2干扰抽水，I号井的水位降深S/等于按Q/1抽水时的有效降 深（即产生出水量Q/1相应的降深)iy与2号井抽水Q/2对1号井产生的附加降深t/1-2之和， 良PS/ly+t/1-2。设1号井的单位出水量(指降深Im时管井在单位时间的出水量)为qi，且 承压水完整井出水量与降深呈直线关系，则由有效降深少Iy,得干扰抽水的出水量为Q^i = qiS' IyD
[0029] 单独抽水时若水位降深等于S/，则出水量化=qiS/ =ql(S/iy+t/i-2)。因此根据式 (1 )，1号井受2号井抽水影响的出水量减少系数为
[0030]
(8)
[0031] 式中各符号含义同前。
[0032] 可见，式(8)与式(3)相同，即采用方法二分析出水量减少系数，所得结论与方法一 必定是相同的。因此，目前采用式(6)、式(7)确定干扰抽水出水量减少系数是不正确的。 H、
【发明内容】
[0033] 为了克服目前流量削减法无法确定干扰抽水的附加降深、不能准确计算出水量减 少系数及干扰出水量的问题，本发明提出承压水完整井井群互阻计算中流量削减法的一种 求解方法，即已知干扰抽水的降深^ 1，准确求解出水量减少系数〇1，进而由Qi计算干扰抽水 的出水量。
[0034] 本发明采用如下技术方案。
[0035] (1)建立干扰抽水附加降深i-j与单独抽水附加降深ti-j的关系式
[0036] 由稳定流承压水完整井公式可得2号井抽水在1号井产生的附加降深的计算 式口'引：
[0037] (9)
[00 測
[0039]式（9)、式（10)中：Q/ 2--干扰抽水时2号井出水量，m3/s ; R--影响半径，m; ri-2一一1、2号井井孔中屯、之间的水平距离，m;K-一承压含水层渗透系数，m/d;M-一承压 含水层厚度，m;Q2-一2号井独立抽水时降深等于干扰抽水降深^ 2的出水量，m^d;其他符 号含义同前。
[0040]利用式(2)，将式(9)中分2由化表示，并根据式（10)，得： 幽]
(…
[0042] 依此类推，可得：
[0043] t'i-j = (l-aj)ti-j (12)
[0044] 其中；
[00 例
(13)
[0046] 式中：Qj-一单独抽水时降深等于干扰抽水降深^北勺出水量，m3/d，j = l，2,…，n; ri-j-第i、j号井井孔中屯、之间的水平距离，m，i = l，2,…，n，j = l，2,…，n，i声j;其他符 号含义同前。
[0047] 式（12)、式（13)分别包含了 nX(n-l)个方程。
[004引因 Qj是降深等于干扰抽水降深非寸的单独抽水的出流量，则依据裘布依公
式口'8]，得： (14)
[0049]
[0050] 式中：邮一一第j号井过滤器的半径，111〇 = 1，2,3，。'，11;其他符号含义同前。
[0051] (2)建立实用的干扰抽水出水量减少系数与降深关系的方程组
[0052] 前已叙及，式(4)所反映的干扰抽水出水量减少系数与降深的关系不能用于实际 计算口-"。
[0053] 将式（12)t/W带入方巧组（4)，并蒋理得：
[0054]
(15)
[0055] 式中各符号含义同前。
[0056] 式（15)即为能用于实际计算的干扰抽水出水量减少系数与降深关系的方程组。该 式对于承压水完整井井群，各井降深、井径相同或不同的情况均是适用的。对于干扰抽水， 当已知各井的降深^ 1，少2,…少。时，只要由式（13)确定了ti-j，利用方程组（15)，贝阿求解 日i。
[0057] (3)建立干扰抽水出水量减小系数的矩阵解析式 [005引为便于求解Qi，将方程组（15)写成矩阵形式： W 别
(16)
[0060]巧
[0061 ]则：
[0062] Q=A-Ir (17)
[0063] 式中:A^i-方阵A的逆阵。
[0064] 式（17)即为干扰抽水出水量减少系数的矩阵解析式，借助Mtalab软件W计算逆阵 A^i,进而易于求解出水量减少系数矩阵曰。
[0065] 至此，当求得了干扰出水时各井的出水量减小系数〇1，则可由式(2)计算干扰出水 量分i。
[0066] 本项发明的有益效果是，针对承压水完整井井群干扰抽水，解决了 W往的流量消 减法无法准确计算出水量减小系数、进而不能正确确定干扰出水量的问题;建立了实用的 干扰抽水时各井出水量减少系数与水位降深关系的方程组，由其可计算各井的出水量减少 系数日1;建立了出水量减少系数日1矩阵解析式，借助Matlab软件计算各井出水量减少系数， 计算效率高。本项发明已利用有关案例进行了验证，计算易于实现、效果很好。 四、
【附图说明】
[0067] 附图是承压水完整井干扰抽水互阻影响示意图。
[0068] 图中：数字1、2为1、2号井的井号;两条虚线分别为1号井抽水流量Q/1、2号井抽水 流量皆2形成的有效水位降落曲线;实线为两井干扰抽水形成的水位降落曲线;r〇i、r〇2分别 为第1、2号井过滤器的半径;ri-2为1、2号井井孔中屯、之间的水平距离，m;图中其他符号含义 同前。 五、
【具体实施方式】
[0069] 案例W :拟在中细砂承压含水层中建造井径400mm的完整井4眼，直线排列，相邻两 眼井的井距为100m。已知承压含水层厚度为M = 20m，渗透系数7m/d，井的影响半径为350m。 若共同抽水时各井设计水位降深为5m，试求共同抽水时各井的出水量。
[0070] (1)采用本发明提出的流量削减法的求解方法
[0071] 1)由式（14)计算各井单独抽水时降深等于干扰抽水降深少j = 5m( j = l，2,3,4)时 的出水量：
[0072]
[0073] 2)根据案例所给资料，由式(13)计算td，并得矩阵A、B为：
[0074]
[0075] 由式（17)，解得日1 =日4 = 0.1835，日2 =日3 = 0.3135。
[0076] 3)由式（2)计算各井的干扰出水量为:Q' 1 = Q' 4 = 481.2m^d，Q' 2 = Q' 3 = 404.6m^ d。
[0077] (2)采用水位削减法
[0078] 采用水位削减法建立方程注：W往文献中水位削减法干扰出水量符号为化， 本文采用前述引入的符号分i，i = l，2,…，n):
[0079
(18)
[0080」联解方程组（18)，得:Q' 1 = Q' 4 = 481. ImVd，Q' 2 = Q' 3 = 404.6mVd (本发明方法Q' 1 =Q' 4=481.2m3/d，水位削减法Q' I = Q' 4=481. ImVd,属舍入影响）。
[0081] 因此本发明方法与水位削减法计算结果相同，进一步验证了本发明方法的正确 性。
[0082] (3)文献[2,3]方法的计算结果
[0083] 按文献[2,3]方法，根据案例所给资料，文献[2]得:Q/I = Q/4=451.7mVd，Q/2 = 9^3 = 378.6m^d，与水位削减法结果不同，进一步说明W往文献出水量减少系数的公式及 解法是不正确的。
[0084] 有关文献
[00化][1化B 50296-2014管井技术规范.北京：中国计划出版社，2014:10-13.
[0086] [2]刘自放，张廉均，邵丕红.水资源与取水工程，北京：中国建筑工业出版社， 2000:91-94.
[0087] [3]徐得潜主编.水资源利用与保护，北京:化学工业出版社，2013:151-153.
[0088] [4]严賊世，范逵初主编.给水工程(第四版），北京：中国建筑工业出版社，1999: 162-166.
[0089] [5]上海市政工程设计研究院主编，给水排水设计手册第3册《城镇给水》（第二 版）.北京：中国建筑工业出版社，2004:125.
[0090] [6]虎胆?吐马尔白.地下水利用（第4版）.北京：中国水利水电出版社，2008:55- 56.
[0091] [7]左建，溫庆博主编.工程地质与水文地质，北京：中国水利水电出版社，2004: 106-107.
[0092] [引张子贤，袁涛编著.水资源与取水工程，北京:化学工业出版社，2016.
[0093] [9]谢中华，李国栋，刘焕进，等.MATLAB从零到进阶[M].北京:北京航空航天大学 出版社，2012.275-276.
【主权项】
1.承压水完整井井群互阻计算中流量削减法的一种求解方法，其特征是利用式（1)计 算干扰抽水时各井的出水量减少系数ai (i = l，2,…，η)，进而计算干扰出水量；(1) 式中：S'i-一干扰抽水时第i号井的水位降深，或简称降深（下同），m，i = l，2,…，η; a,一一干扰抽水时第i号井受其他各井抽水共同影响的出水量减少系数，i = l，2,…，η; ti-j--单独抽水时第j号井出水量Qj在第i号井产生的附加降深，m，i = l，2,…，n，j = l， 2，…，η，i辛j，Qj为独立抽水时第j号井降深等于干扰抽水降深S' j时的出水量，m3/d，j = 1， 2,···，n; 其特征在于，具体操作步骤为： 步骤1，计算独立抽水时第j号井降深等于干扰抽水降深S、时的出水量^式中：K一一承压含水层渗透系数，m/d;M-一承压含水层厚度，m;R-一影响半径，m; r〇J一一第j号井过滤器的半径，111;」=1，2,3，一，11;其他符号含义同前； 步骤2,计算第j号井出水量仏在第i号井产生的附加降深式中：ri-j--第i、j号井井孔中心之间的水平距离，m，i = l，2,…，n，j = l，2,…，n，i辛 j ;其他符号含义同前； 式⑶包含了 nX(n-l)个方程； 步骤3，已知干扰抽水时各井的水位降深S' i，解方程组式（1)计算干扰抽水各井的出水 量减少系数ai，? = 1，2，···，η，其特征是： α=Α-咕 (4)Α-1一一方阵 A的逆阵;矩阵a、A、B中的各符号含义同前； 步骤4,根据单独抽水的出水量Qi及干扰抽水各井的出水量减少系数ai，计算各井的干 扰出水里Q i = Qi(l_ai)，i = l，2r"，n。
【文档编号】G06F19/00GK105956366SQ201610250446
【公开日】2016年9月21日
【申请日】2016年4月15日
【发明人】张子贤
【申请人】江苏建筑职业技术学院