《干扰井群计算分析程序》编制.pdf

《干扰井群计算分析程序》编制.pdf

《《干扰井群计算分析程序》编制.pdf》由会员分享，可在线阅读，更多相关《《干扰井群计算分析程序》编制.pdf（44页珍藏版）》请在冰豆网上搜索。

干扰井群计算分析程序干扰井群计算分析程序-张焕智干扰井群计算分析程序干扰井群计算分析程序编制：

张焕智编制：

张焕智1编制依据1编制依据1.1稳定流理论计算公式1.1.1承压完整干扰井计算公式=njiijiiirRQKMS1,lg73.21有+=+=+=nnnnnnnnnnnnnrRQrRQrRQKMSrRQrRQrRQKMSrRQrRQrRQKMSlglglg73.21lglglg73.21lglglg73.2122211122222211121122211111LLLLLLLLLLLLLLLLLLLLLLLL式中Si：

i井干扰水位降深（m）Qi：

i井干扰出水量（m3/d）K：

含水层渗透系数（m/d）M：

含水层厚度（m）Ri：

i井影响半径（m）rij：

i井至j井距离，i=j时为该井半径r0（m）n：

布井数1.1.2潜水完整干扰井计算公式=njiijiiiirRQKSSH1,lg366.11）2（式中，H为从含水层底板算起的含水层厚度。

1.1.3承压非完整干扰井计算公式干扰井群计算分析程序干扰井群计算分析程序-张焕智+=+=+=CmnnnnnnnnnnCmnnnCmerRQrRQrRQSKMrRQerRQrRQSKMrRQrRQerRQSKMlglglg73.2lglglg73.2lglglg73.222211122222211121122211111LLLLLLLLLLLLLLLLLLLLLLLL式中AlMrMllMCmiiiii24lg=AlMMlMrllMeriiiiiiCmii217.04lg）1（434.0lg3.23.1）lg（+=）125.01（）875.01（）125.0（）875.0（lg）（CCCCCfA=MlCi=式中，li为i井过滤器长度，M为含水层厚度，（C）为伽马函数。

1.1.4潜水非完整干扰井计算公式在厚度较大的潜水含水层地区，一般不需揭穿整个含水层，因此单井抽水时可采用下式计算水位降深（非淹没过滤器）：

01lg73.2rRKHQS=式中，H1为至过滤器底部的含水层厚度。

按照水位叠加原理，其干扰井群公式为=njiijiiirRQKHS1,1lg73.211.2非稳定流理论计算公式1.2.1承压完整干扰井计算公式（同步抽水）=njiijiiratQS1,225.2ln41上式可展开如下：

干扰井群计算分析程序干扰井群计算分析程序-张焕智+=+=+=22222222225.2ln25.2ln25.2ln425.2ln25.2ln25.2ln425.2ln25.2ln25.2ln422112222211211221111nnnnnnnnnnratQratQratQSratQratQratQSratQratQratQSLLLLLLLLLLLLLLLLLLLLLLLL式中Si：

i井水位降深（m）Qi：

i井出水量（m3/d）T：

含水层导水系数（m2/d）a：

含水层压力传导系数（m2/d）t：

抽水时间（同步,d）rij：

i井至j井距离，i=j时为该井半径r0（m）n：

布井数依据以上公式，按各井分配出水量Qi，可计算任意井（检验点）产生的水位降深Si；另外，当给定各井水位降深Si，可计算各井出水量Qi，进而求得干扰井群总干扰出水量。

1.2.2潜水完整干扰井计算公式（同步抽水）=njiijiiiratQKSSH1,225.2lg21）2（上式可展开如下：

（）（）（）+=+=+=22222222225.2ln25.2ln25.2ln2225.2ln25.2ln25.2ln2225.2ln25.2ln25.2ln222211222221122112211111nnnnnnnnnnnratQratQratQSSHKratQratQratQSSHKratQratQratQSSHKLLLLLLLLLLLLLLLLLLLLLLLL式中，K为含水层渗透系数（m/d），H1为含水层厚度（m）。

1.3实验干扰井计算（干扰井计算的相关分析法）*3干扰井群计算分析程序干扰井群计算分析程序-张焕智试验资料表明，管井抽水中出水量Q与其相应的水位降深S一般都具有一定的函数关系，如对数：

Q=A+BlgS

（1）幂：

Q=ASB

（2）抛物线：

Q=AS+BS2（3）双曲线：

1/Q=A+B/S（4）等，这种函数关系均可通过单井抽水试验确定下来。

另外，当抽水井以一定流量抽水时，观测孔水位降深与抽水井出水量Q以及距离d的多元关系可综合表述为下式：

t=b0+b1lgQ+b2（1+lgQ）lgd（5）或t=b0+b1lgQ+（b2+b3lgQ）lgd（6）上式亦可通过多孔抽水试验确定下来。

如上，当有n口井同时进行抽水时，抽水井中水位降深不仅受本井出水量影响，而且还将同时受到其它各抽水井出水量的影响，故各井水位降深可表示为:

+=+=+=nnnnnnnSSSSSSSSSSSSLLLLLLLLLLLLLLLLLLLLL21222212112111其中，抽水井降深可以由

（1）式确定，其它井对抽水井所产生的干扰水位降可以由（5）或（6）式确定文中均以（5）式为例叙述，则有+=+=+=）（）,（）,（）,（）（）,（）,（）,（）（2211222112112211nnnnnnnnQfdQfdQfSdQfQfdQfSdQfdQfQfSLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL（7）通过求解（7）式，即可直接计算各干扰井的降深。

当给定各井降深后，计算各井出水量时，继续展开（7）式有干扰井群计算分析程序干扰井群计算分析程序-张焕智+=+=+=nnSSSBlgQA）lgdlgQ（1blgQbb）lgdlgQ（1blgQbb）lgdlgQ（1blgQbbBlgQA）lgdlgQ（1blgQbb）lgdlgQ（1blgQbb）lgdlgQ（1blgQbbBlgQAn222210n1121102nn2n102211211021nn2n10122221011LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL整理上式+=+=+=n2n221n2201n121n120n2n212n2021212121202n1n211n2021221122011BlgQA）lgQlgdbb（lgdbb）lgQlgdb（blgdbb）lgQlgdbb（lgdbbBlgQA）lgQlgdb（blgdbb）lgQlgdbb（lgdbb）lgQlgdb（blgdbbBlgQALLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLnSSS则+=+=+=+n2n2211n121n220n120n2n212121212n2021202n1n212122111n2012201BlgQ）lgQlgdbb（）lgQlgdb（b）Algdbblgdbb（）lgQlgdbb（BlgQ）lgQlgdb（b）lgdbblgdbb（）lgQlgdbb（）lgQlgdb（bBlgQ）lgdbblgdbbA（LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLnSASS亦即lg）lg（）lg（lg121120jnjiijnjiijiiQdbbdbbQBAS+=（8）或lg）lg（lg）lg（121120jnjiijinjiijiQdbbQBdbbAS+=+（8）上式中i,j=1,2,n。

当取（6）式计算其它井对抽水井的干扰水位降时，（8）或（8）式方程右侧各变量系数（b1+b2lgdij）中的b2改为b3即可。

当用于预测井群出水量时，由于分别求得的Q=f（S）与S=f（Q）相互计算时易出现较大误差，难以取得较好效果。

因此，程序编制中首先预定设计方案下各井出水量，然后依据确定的Q=f（S）及t=f（Q,d）计算预定设计方案下各井可能产生的水位降深；并通过井群干扰降深的反复调整计算，最终确定基本能够满足设计要求的井群布置方案。

2程序功能与应用2程序功能与应用干扰井群计算分析程序干扰井群计算分析程序-张焕智本程序可进行干扰井群抽水试验计算分析。

包括稳定流的承压完整井及潜水非完整井、潜水完整井、承压非完整井的干扰井群计算，非稳定流的承压完整井及潜水完整井干扰井群计算，实验公式法（Q=f（S）及t=f（Q,d））的干扰井群计算；可分别按照预定各井出水量（Q）或降深（S）进行干扰井群计算分析，为开采井群设计提供依据。

本程序使用QBASIC算法语言编制，编译为.EXE可执行文件。

支撑环境主要是：

硬件386以上微机,打印机,鼠标器等；软件DOS3.30以上操作系统,UCDOS3.0以上汉字系统。

本程序可安装在D：

盘目录（文件夹）GRJJS中（其他盘也可以，但需修改批命令中的所在盘符），UCDOS汉字系统与本程序应安装在同一盘中。

批命令GRJQ.BAT如下（程序安装在D：

盘）：

ECHOOFFd:

d:

UCDOSRD16%1d:

UCDOSKNL%2d:

UCDOSPYd:

UCDOSLIMDWBd:

UCDOSRDFNTd:

UCDOSprntd:

UCDOSTXcdgrjjsgrjq.execdd:

ucdosquit本程序有建立数据文件的功能。

计算所用的数据文件也可使用“记事本”或“写字板”，按照数据文件模板（GRJ01.DATGRJ08.DAT）建立并编辑（见图），更加方便。

数据文件存放在D：

GRJJSDAT子目录（文件夹）中。

数据文件模板如下：

干扰井群计算分析程序干扰井群计算分析程序-张焕智grj01.dat：

grj01.dat：

*以上存放承压完整井与潜水井按照预定降深（S）计算数据,其格式（左边无空格）为:

N,M,H,K回车*干扰井数,检验点数（为0）,含水层厚度,渗透系数Z$

（1）,S

（1）,R0

（1）,R

（1）,XY（1,1）,XY（1,2）打回车Z$

（2）,S

（2）,R0

（2）,R

（2）,XY（2,1）,XY（2,2）打回车Z$（N）,S（N）,R0（N）,R（N）,XY（N,1）,XY（N,2）打回车井编号,降深,井半径,影响半径,井坐标X,井坐标Ygrj02.dat：

*以上存放承压非完整井按照预定降深（S）计算数据,其格式（左边无空格）为:

N,M,H,K回车*干扰井数,检验点数（为0）,含水层厚度,渗透系数Z$

（1）,S

（1）,L

（1）,R0

（1）,R

（1）,XY（1,1）,XY（1,2）打回车Z$

（2）,S

（2）,L

（2）,R0

（2）,R

（2）,XY（2,1）,XY（2,2）打回车Z$（N）,S（N）,L（N）,R0（N）,R（N）,XY（N,1）,XY（N,2）打回车井编号,降深,过滤器长度,井半径,影响半径,井坐标X,井坐标Ygrj03.dat：

grj03.dat：

*以上存放承压完整井与潜水井按照预定出水量（Q）计算数据,其格式（左边无空格）为:

N,M,H,K回车*干扰井数,检验点数,含水层厚度,渗透系数Z$

（1）,Q

（1）,R0

（1）,R

（1）,XY（1,1）,XY（1,2）打回车Z$

（2）,Q

（2）,R0

（2）,R

（2）,XY（2,1）,XY（2,2）打回车Z$（N）,Q（N）,R0（N）,R（N）,XY（N,1）,XY（N,2）打回车井编号,降深,井半径,影响半径,井坐标X,井坐标YM$

（1）,MX（1,1）,MX（1,2）打回车M$

（2）,MX（2,1）,MX（2,2）打回车M$（M）,MX（M,1）,MX（M,2）打回车点编号,点坐标X,点坐标