进气值面,井流模型,完整井,Boussinesq方程,解析解,参数反演模型


," /> 进气值面,井流模型,完整井,Boussinesq方程,解析解,参数反演模型


,"/> <p class="MsoNormal"> 考虑毛细饱水带的非承压含水层中仿<span>Theis</span>井流模型及解析解

吉林大学学报(地球科学版) ›› 2025, Vol. 55 ›› Issue (2): 563-574.doi: 10.13278/j.cnki.jjuese.20230151

• 地质工程与环境工程 • 上一篇    下一篇

考虑毛细饱水带的非承压含水层中仿Theis井流模型及解析解

晋浩颖1,2,程大伟1,2,詹红兵3,杨胜科1,2,张晓斐4,张琳1,2   

  1. 1.长安大学水利与环境学院,西安710064

    2.旱区地下水文与生态效应教育部重点实验室(长安大学),西安710064

    3.德克萨斯农工大学地理与地球物理系,美国德克萨斯卡城 TX77843-3115

    4.西安曲江渼陂湖投资建设有限公司,西安710300

  • 出版日期:2025-03-26 发布日期:2025-05-10
  • 基金资助:

    国家重点研发计划项目(2020YFC1808304);国家自然科学基金项目(41602237)


Modified Theis Well Model and Its Analytical Solution in Unconfined Aquifers Considering Capillary Saturated Zone

Jin Haoying1,2,Cheng Dawei 1,2Zhan Hongbing3Yang Shengke1,2Zhang Xiaofei4Zhang Lin1,2   

  1. 1. School of Water and Environment, Chang’an University, Xi’an 710064,China

    2. Key Laboratory of Subsurface Hydrology and Ecological Effects in Arid Region (Chang’an University),Ministry of Education,

    Xi’an 710064,China

    3. Department of Geology and Geophysics, Texas A&M University, College Station, TX 77843-3115, USA

    4. Xi’an Qujiang Meibei Lake Investment and Construction Co., Ltd., Xi’an 710300,China

  • Online:2025-03-26 Published:2025-05-10
  • Supported by:
    Supported by the National Key Research and Development Program of China (2020YFC1808304) and the National Natural Science Foundation of China (41602237)

摘要: 受制于经典饱和渗流理论中自由面定位在潜水面处的局限性,非承压含水层中仿Theis井流模型存在理论基础不坚实、参数物理意义不明确等问题。为克服上述局限性,首先将自由面由潜水面移动到进气值面,以定位在进气值面处的Boussinesq方程为控制方程,修正仿Theis井流模型,并推导模型的解析解;之后,构建参数反演模型,以校准模型参数,并通过分别将修正模型计算的降深和经典模型计算的理论降深与实测降深进行对比,检验修正模型的合理性;然后,利用参数敏感性分析探讨修正模型下进气值面高程的变化特征;最后,讨论进气值面处铅直方向渗流分速度的计算方法和变化特征。结果表明:修正模型计算的降深曲线与实测降深曲线基本吻合;进气值面高程随与井中心的距离、给水度增大而增大,随抽水流量增大而减小。在抽水早期,进气值面高程随渗透系数的增大而减小,在抽水后期则相反。基于进气值面处水量均衡关系(基于线性或非线性方程,采用完整解或近似解)及进气值面处渗流连续性方程(基于线性方程,采用完整解)等5种情况,推导的铅直方向渗流分速度的解析表达计算结果基本吻合,曲线趋势一致,表现为随抽水时间的延长呈非线性增长。

关键词: 进气值面')">

进气值面, 井流模型, 完整井, Boussinesq方程, 解析解,
')">参数反演模型


Abstract:

Due to the limitation of setting the upper boundary at the free surface located at the water table, the Theis well model in unconfined aquifers has some problems such as unsound theoretical basis and unreasonable physical meaning of parameters. To overcome these issues, firstly, the well model is revised by moving the free surface from the water table to the air entry plane, and takes the Boussinesq equation located at the air entry plane as the governing equation, and derives the analytical solution of the model; Then, a parameter inversion model is constructed to calibrate the model parameters. By comparing the revised model drawdown and the theoretical drawdown of the classical model with the measured drawdown, it verifies the rationality of the well revised model proposed in this paper. Then, the variation characteristics of the air entry plane elevation under the modified well model are discussed by using parameter sensitivity analysis. Finally, the computational method and variation characteristics of seepage velocity in the vertical direction of the air entry plane are discussed. It is found that the drawdown curve obtained by the revised model coincides with the measured drawdown curve; The air entry plane elevation increases with increasing distance from the center of the well and the specific yield and decreases with increasing pumping flow rate. In the early stage of pumping, the air entry plane elevation decreases with increasing hydraulic conductivity, but the opposite is true in the late stage of pumping. The calculation results of the analytical expressions of the seepage velocity in the vertical direction based on the water balance relation at the air entry plane (based on the linearized or non-linearized equation, using complete solutionor approximate solution) and based on continuity equation of seepage flow at the air entry plane (based on the linearized equation, using complete solution) under five conditions are basically consisent, and the trends of the curve are consistent with each other, showing a nonlinear increase with the prolongation of pumping time.

Key words: air entry plane, well model, fully penetrating well, Boussinesq equation, analytical solution, parameter inversion model

中图分类号: 

  • P641.2
[1] 王锦博, 贺志毅, 王亚楠, 程诗雨, 刘明柱. 某污染场地土壤与地下水协同修复阈值[J]. 吉林大学学报(地球科学版), 2025, 55(4): 1288-1297.
[2] 雷怡, 帅飞翔, 孙红月, 张泽坤, 熊超. 边坡负压排水非稳定流研究[J]. 吉林大学学报(地球科学版), 2022, 52(1): 162-.
[3] 刘国庆, 吴时强, 范子武, 周志芳, 谢忱, 乌景秀, 柳杨. 回灌与回扬物理过程的解析推导及灌压变化规律[J]. 吉林大学学报(地球科学版), 2016, 46(6): 1799-1807.
[4] 汤井田, 辛会翠, 王冉. 点电源下复杂角域地形影响及校正[J]. J4, 2012, 42(1): 254-261.
[5] 叶栋成,慕山, 陶月赞. 地下水补给对河流水质模型的影响[J]. J4, 2008, 38(4): 644-0648.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 崔建军, 刘晓春, 胡娟, 曲玮. 桐柏杂岩中印支期变质岩包体的变质作用[J]. J4, 2009, 39(4): 618 -629 .
[2] 尤敏鑫,刘建民. 同位素地球化学在峨眉山大火成岩省研究中的应用现状与进展[J]. 吉林大学学报(地球科学版), 2014, 44(4): 1231 -1243 .
[3] 杨春梅, 李洪奇,陆大卫,张方礼,高 原,邵英超. 不同驱替方式下岩石电阻率与饱和度的关系[J]. J4, 2005, 35(05): 667 -671 .
[4] 祝洪臣,张炯飞,权 恒. 大兴安岭中生代两期成岩成矿作用的元素、同位素特征及其形成环境[J]. J4, 2005, 35(04): 436 -0442 .
[5] 朱建伟, 赵刚, 刘博, 郭巍, 成俊. 油页岩测井识别技术及应用[J]. J4, 2012, 42(2): 289 -295 .
[6] 陈力,梁海安,张文娟,荣帆. 模糊数学方法在城市工程地质环境区划中的应用--以抚顺市城区为例[J]. J4, 2008, 38(5): 837 -0840 .
[7] 高桂梅,苏 克,王文颖,甘树才,刘招君. 吉林省桦甸油页岩中稀土元素和微量元素的研究[J]. J4, 2006, 36(6): 974 -0979 .
[8] 吴孔运,蒋忠诚,叶 晔. 不同植物群落对灰岩试块溶蚀速率的影响[J]. J4, 2007, 37(5): 967 -0971 .
[9] 周彦章,迟宝明,刘中培. 山东夏甸金矿床充水机理构造控制模式[J]. J4, 2008, 38(2): 255 -0260 .
[10] 张渊,刘连登,孙景贵,陈国华,张洪喜,闫复传,杨开春. 胶东西北部黄埠岭金矿床两期次叠加成矿[J]. J4, 2008, 38(1): 21 -0026 .