孔隙尺度多孔介质流体流动与溶质运移高性能模拟  

Pore-scale simulations of fluid flow and solute transport in porous media by high-performance Lattice Boltzmann Method

作　　者：周鸿翔 郑延丰 吴劳生 陈铖 曾令藻 ZHOU Hongxiang;ZHENG Yanfeng;WU Laosheng;CHEN Cheng;ZENG Lingzao(College of Metrology and Measurement Engineering,China Jiliang University,Hangzhou 310018,China;Department of Mining and Minerals Engineering,Virginia Tech,Blacksburg,VA 24061,USA;College of Environmental and Resource Sciences,Zhejiang University,Hangzhou 310058,China;College of Civil Engineering and Architecture,Zhejiang University,Hangzhou 310058,China;Department of Environmental Sciences,University of California,Riverside,CA 92521,USA)

机构地区：[1]中国计量大学计量测试工程学院,浙江杭州310018 [2]Department of Mining and Minerals Engineering,Virginia Tech,Blacksburg,VA 24061,USA [3]浙江大学环境与资源学院,浙江杭州310058 [4]浙江大学建筑工程学院,浙江杭州310058 [5]Department of Environmental Sciences,University of California,Riverside,CA 92521,USA

出　　处：《水科学进展》2020年第3期422-432,共11页Advances in Water Science

基　　金：国家重点研发计划资助项目(2018YFC1800503);国家自然科学基金资助项目(41571215)。

摘　　要：深入探究孔隙尺度下的流体流动特性和溶质运移规律对石油开采、农田养分管理、地下水污染修复有着重要意义。以人工构建的多孔介质结构和同步辐射X射线显微CT扫描的土壤团聚体(分辨率3.7μm)为研究对象,在空间节点数多达64 000 000的情况下,基于格子Boltzmann模型和GPU并行技术计算得到多孔介质流体运动和溶质运移过程的关键参数,并据此探究多孔介质空间异质性对水力学特性的影响。通过对3组不同结构的多孔介质比较发现,结构复杂程度最高的土壤样品和不规则堆叠的圆球结构的渗透率在100 mD(即10^-13m^2)量级,远低于规则堆叠的圆球结构(>20 000 mD);土壤的迂曲度为1.40~1.60,明显高于规则堆叠的圆球结构。研究结果表明,渗透率大的样品具有较小的迂曲度,这与结构的空间异质性有较强的关系;土壤的渗透率和迂曲度呈现各向异性;在水力梯度一定的前提下,渗透率较大的样品,纵向弥散系数也较大;同时,结构的异质性也会影响溶质的穿透曲线。本研究提出的模拟方法可在土壤结构中进行高效的水流运动和溶质运移模拟,可用于土壤多孔介质在孔隙尺度下的水力学特性研究。Understanding the mechanism of fluid flow and solute transport at the pore scale is of great importance for oil recovery,crop nutrient management and groundwater pollution restoration. This study employed lattice Boltzmann model combining with GPU parallel technology to investigate the porous media of computer-generated structures and synchrotron-based X-ray micro-CT scans of soil aggregates(resolution 3.7 μm). The key parameters of fluid flow and solute transport in the porous media were obtained,and the influence of spatial heterogeneity of porous media on hydraulic properties was explored by high-performance simulation(spatial nodes up to 64 000 000). By comparing the three groups of porous media with different structures,it was found that the permeabilities of the soil sample with the highest structural complexity and beads irregularly stacked are on the order of 100 mD(i.e. 10^-13m^2),which is much lower than that of the regularly stacked beads(>20 000 mD);The soil sample has a tortuosity of 1.40~1.60,which is significantly higher than that of the regularly stacked beads. Our results show that the porous media with high permeabilities have small degree of tortuosity,indicating that the permeabilities of porous media are related to the spatial heterogeneity of the structure. The permeability and tortuosity of soil aggregate are anisotropic. At given pressure gradient,the longitudinal diffusion coefficient is greater for a sample with higher permeability. The heterogeneity of the pore structure also affects the breakthrough curve. The method established in this work can simulate water flow and solute migration in real soil structure,and can be used to study the hydraulic characteristics of porous media at the pore scale.

关 键 词：孔隙尺度 格子Boltzmann法 流体流动 溶质运移 

分 类 号：TV139.1[水利工程—水力学及河流动力学]