非承压含水层定水头抽水两区井流数值模型研究  

Numerical model of two-region flow by constant-head pumping in an unconfined aquifer

作　　者：史鹏钰 宗一杰 滕开庆 刘健军 肖良 SHI Pengyu;ZONG Yijie;TENG Kaiqing;LIU Jianjun;XIAO Liang(College of Civil Engineering and Architecture,Guangxi University,Nanning 530004,China;Guangxi Key Laboratory of Disaster Prevention and Engineering Safety,Nanning 530004,China;Key Laboratory of Disaster Prevention and Structural Safety of Ministry of Education,Guangxi University,Nanning 530004,China)

机构地区：[1]广西大学土木建筑工程学院,广西南宁530004 [2]广西防灾减灾与工程安全重点实验室,广西南宁530004 [3]广西大学工程防灾与结构安全教育部重点实验室,广西南宁530004

出　　处：《煤田地质与勘探》2023年第10期124-133,共10页Coal Geology & Exploration

基　　金：国家自然科学基金青年基金项目(41807197);广西自然科学基金项目(2018GXNSFAA138042);河北高层次人才资助项目(B2018003016)。

摘　　要：为了揭示在非承压含水层中定水头抽水试验引起的达西−非达西两区流动机理,提出基于有限差分法的地下水定水头抽水井流数值模型。该模型根据抽水的流态特征将含水层分为2个区域:靠近抽水井的有限非达西渗流区域和远离抽水井的半无限达西渗流区域,其中非达西流区域流态的模拟基于Izbash方程实现。通过与COMSOL Multiphysics的有限元数值解进行比较,验证了所提出数值解的可靠性。最后,研究有限非达西流效应对水头和抽水井抽水速率的影响以及井内水头对抽水井抽水速率的影响。研究表明:在抽水试验中非达西区域的影响不可忽略,湍流会分别导致两区流中水头较纯非达西流和纯达西流的水头偏大和偏小,且随抽水时间的增加逐渐变大;通过减小抽水井井内水头或增大非达西系数可提高抽水速率,但该影响会随抽水时间的增加而逐渐减弱;断面流量随径向距离的增大而不断减小,断面流量与径向距离曲线下降速率不断减小,且在转换界面处会出现转折点。该模型为定量研究在非达西流和达西流耦合作用下抽水井附近的井流水头特征提供了一种简洁的方法,并为调查定水头抽水测试期间的抽水速率提供理论依据。In order to reveal the mechanism of two-region Darcian and non-Darcian flow induced by constant-head pumping tests in unconfined aquifers,a numerical model of well flow by constant-head pumping underground based on finite difference solution was proposed.In the model,it is assumed that the aquifer is divided into two regions according to the pumping flow characteristics:the finite non-Darcian flow region near the pumping well and the semi-infinite Darcian flow region far away from the pumping well.Specifically,the flow regime in the non-Darcian flow region was simulated by the Izbash’s equation,and the reliability of the proposed solution was verified by comparison with finite element numerical solutions by COMSOL Multiphysics.Finally,the influence of finite non-Darcian effect on hydraulic head and pumping rate,as well as the influence of hydraulic head in the pumping well on pumping rate,was especially studied.The results show that the influence of non-Darcian flow region in the pumping test cannot be ignored.The turbulent flow makes the hydraulic head of the two-region flow larger than that of the pure non-Darcian flow and smaller than that of pure Darcian flow.Besides,such difference of the hydraulic head is increased with the pumping time.The pumping rate can be increased by reducing the hydraulic head in the pumping well or increasing the non-Darcian coefficient,but the effect is gradually decreased as the pumping continues.The flow rate at the cross-section is decreased with the increase of radial distance,the gradient of the flow rate-radial distance curve is decreased with time,and a turning point will appear at the conversion interface between the two regions.The proposed model provides a simple method for quantitative studies on the characteristics of hydraulic head near the pumping well under the coupling of non-Darcian and Darcian effects,and provides a theoretical basis for determining the pumping rate during constant-head pumping test.

关 键 词：非承压含水层 定水头抽水 两区流动 非达西效应 有限差分解 数值模型 

分 类 号：P345[天文地球—水文科学]