地学前缘 ›› 2024, Vol. 31 ›› Issue (3): 59-67.DOI: 10.13745/j.esf.sf.2023.9.32

• 沉积作用与沉积环境 • 上一篇    下一篇

不同胶结度的多孔介质中溶质横向弥散的孔隙尺度模拟研究

侯玉松1,*(), 胡晓农1, 吴吉春2,*()   

  1. 1.济南大学 水利与环境学院, 山东 济南 250022
    2.南京大学 地球科学与工程学院, 江苏 南京 210023
  • 收稿日期:2023-08-10 修回日期:2023-09-30 出版日期:2024-05-25 发布日期:2024-05-25
  • 通信作者: *侯玉松(1987—)(第一作者),男,博士,主要从事地下水数值模拟、渗流及溶质运移等研究。E-mail: stu_houys@ujn.edu.cn; 吴吉春(1968—),男,教授,博士生导师,主要从事水文地质学及地下水数值模拟的研究。E-mail: jcwu@nju.edu.cn
  • 基金资助:
    国家自然科学基金项目(42002257);山东省自然科学基金项目(ZR2020QD123);山东省新高校院所创新团队项目(2021GXR070);山东省高校院所创新团队项目(2018GXRC012)

Pore scale simulation study of transverse dispersion of solute in porous media with different cementation degrees

HOU Yusong1,*(), HU Xiaonong1, WU Jichun2,*()   

  1. 1. School of Water Conservancy and Environment, University of Jinan, Jinan 250022, China
    2. School of Earth Sciences and Engineering, Nanjing University, Nanjing 210023, China
  • Received:2023-08-10 Revised:2023-09-30 Online:2024-05-25 Published:2024-05-25

摘要:

深入理解、精确刻画溶质孔隙尺度横向弥散行为是揭示溶质运移现象内在机制、提高溶质运移模拟精度的关键。基于孔隙尺度模拟方法研究了不同胶结度的多孔介质中溶质横向弥散特征。结果表明:当胶结度增大时,溶质横向弥散程度有所增强,但与纵向弥散相比增强幅度较小,且在模拟的胶结度范围内横向弥散未展现显著的尺度效应。研究还发现多孔介质中横向弥散系数DT与纵向弥散系数DL间的差异程度不是固定不变的,而是随着溶质运移距离逐渐增大最终达到稳态的过程,并且DL/DT的渐进值及达到渐进值所需时间均受胶结度的影响。因此,将DT视为比DL小一个数量级的值与溶质实际弥散行为存在较大偏差。同时,通过对比不同介质中流体流速横向、纵向分量概率密度及溶质对流速率概率密度的差异,揭示了以上溶质横向弥散特征变化及其与纵向弥散特征变化不同的内在机制。

关键词: 横向弥散, 多孔介质, 胶结度, 孔隙尺度, 数值模拟

Abstract:

Deeply understanding and accurately characterizing the transverse dispersion behavior of solutes at the pore scale is crucial for unveiling the internal mechanisms of solute transport phenomena and enhancing the precision of solute transport simulations. In this study, utilizing a pore-scale simulation method, we have investigated the transverse dispersion behavior of solutes in porous media with varying degrees of cementation. The findings of this research are as follows: The transverse dispersion of solutes increases with the cementation degree, although the degree of enhancement is lower compared to longitudinal dispersion. The transverse dispersion does not exhibit a significant scale effect within the range of simulated cementation degrees. The study revealed that the difference between the transverse diffusion coefficient (DT) and the longitudinal diffusion coefficient (DL) in all media is not constant but gradually increases, eventually reaching a steady state with solute transport. Both the asymptotic value of DL/DT and the time required to reach this value are influenced by the cementation degree. There is a notable discrepancy between assuming DT to be one order of magnitude smaller than DL and the actual solute dispersion behavior. By comparing the probability density of the lateral and longitudinal components of flow velocity and the probability density of the advective rate of solutes in different media, this paper elucidates the internal mechanisms underlying the variations in solute transverse dispersion characteristics and highlights the differences between transverse and longitudinal dispersion behaviors.

Key words: transverse dispersion, porous media, cementation degree, pore scale, numerical simulation

中图分类号: