地球科学进展 doi: 10.11867/j.issn.1001-8166.2023.010

   

基于负表皮层影响的径向溶质运移模型构建与新求解方法
张开鑫 1,黄璟胜 1,2*,王晨 1,童晨晨 1,王子成 1   
  1. (1. 河海大学水文水资源学院,江苏 南京 210098;2. 河海大学长江保护与绿色发展研究院,江苏 南京 210098)
  • 通讯作者: 黄璟胜(1984-),男,中国台湾新北人,教授,主要从事地下水机理研究与数值模拟. E-mail:cshuang0318@hhu.edu.cn
  • 基金资助:
    国家自然科学基金项目“耦合大—小尺度介质地下水流的滞后理论”(编号:52079042)和“淮河流域遥感土壤水分数据产品验证及暴雨洪水数据同化系统构建”(编号:41830752)资助.

A New Analytical Method for Modeling Radially Divergent Solute Transport in Two-zone Confined Aquifers with Negative Skin Effects#br#

ZHANG Kaixin 1, HUANG Chingsheng 1, 2*, WANG Chen 1,TONG Chenchen 1, WANG Zicheng 1   

  1. (1. College of Hydrology and Water Resources, Hohai University, Nanjing, 210098, China; 2. Yangtze Institution for Conservation and Development, Hohai University, Nanjing 210098, China
  • Contact: HUANG Chingsheng (1984-), male, New Taipei City, Taiwan Province, China, Professor. Research areas include groundwater mechanism and numerical simulation. E-mail: cshuang0318@hhu.edu.cn
  • About author:ZHANG Kaixin (1998-), female, Lianyungang City, Jiangsu Province, Master student. Research areas include groundwater mechanism and numerical simulation. E-mail: zkx_0111@qq.com
  • Supported by:
    Project supported by the National Natural Science Foundation Program of China “The lagging theory coupling groundwater flows in large- and small-scale media” (Grant No. 52079042) and “Verification of remote sensing soil moisture data products and construction of rainstorm and flood data assimilation system in the huaihe river basin” (Grant No. 41830752).
对于模拟示踪试验的溶质运移,传统方法采用细网格离散小尺度表皮层,存在网格数量多和计算时间长等问题。由于表皮层外缘的参数变化剧烈,即使采用细网格,仍然存在显著误差。因此发展基于完整井的新示踪试验溶质运移模型,提出新的瞬态Robin 边界条件反映表皮层吸附/解吸溶质的影响,实现表皮层的无网格。应用拉普拉斯变换导出模型的解析解,通过有限元法建立基于非完整井的数值解。结果显示新瞬态Robin 边界条件精准反映负表皮层的影响。对于求解反演问题,表皮层宽度的估计范围由0.45~0.54 m缩至0.47~0.48 m,纵向弥散度由0.6~10.0 m缩至6.4~7.7 m,参数估计的可靠度显著提升。表皮层的无网格节省97%的数值解计算时间。
Various models have been developed for radially divergent tracer tests in two-zone confined aquifers of the skin and formation zones. However, existing numerical solutions require considerable computing time because of the fine spatial discretization of skins. The abrupt change in parameters near the skin-formation interface produces significant errors while predicting the spatiotemporal concentration near the interface, despite fine spatial discretization. In this study, a new model was developed for conducting radially divergent tracer tests in a partially penetrating well in a two-zone-confined aquifer. The skin was treated as a new transient Robin boundary condition specified at the skin-formation interface to reflect the effect of solute adsorption/release in the
skin and achieve no skin discretization. A finite element solution for the model was developed. The analytical solution of the model modified for full penetration of the well was developed using the Laplace transform. These results suggest that the transient Robin boundary condition leads to accurate concentration predictions affected by negative skin. The analytical solution predicts reliable ranges of 0.47~0.48 m for the skin width w and 6.4~7.7 m for the longitudinal dispersivity α'l, whereas a traditional solution exhibits a range of 0.45 m≤w≤0.54 m and 0.6 m≤α'l≤10 m. The finite element solution required only 3% of the computing time for obtaining a finite element solution based on fine skin discretization. In conclusion, this study provides implications not only
for theoretical advances but also for useful numerical methods.

中图分类号: 

[1] 张开鑫, 黄璟胜, 王晨, 童晨晨, 王子成. 基于负表皮层影响的径向溶质运移模型构建与新求解方法[J]. 地球科学进展, 2023, 38(4): 429-440.
阅读次数
全文


摘要