首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 233 毫秒
1.
深入探究孔隙尺度下的流体流动特性和溶质运移规律对石油开采、农田养分管理、地下水污染修复有着重要意义。以人工构建的多孔介质结构和同步辐射X射线显微CT扫描的土壤团聚体(分辨率3.7μm)为研究对象,在空间节点数多达64 000 000的情况下,基于格子Boltzmann模型和GPU并行技术计算得到多孔介质流体运动和溶质运移过程的关键参数,并据此探究多孔介质空间异质性对水力学特性的影响。通过对3组不同结构的多孔介质比较发现,结构复杂程度最高的土壤样品和不规则堆叠的圆球结构的渗透率在100 mD(即10^-13m^2)量级,远低于规则堆叠的圆球结构(>20 000 mD);土壤的迂曲度为1.40~1.60,明显高于规则堆叠的圆球结构。研究结果表明,渗透率大的样品具有较小的迂曲度,这与结构的空间异质性有较强的关系;土壤的渗透率和迂曲度呈现各向异性;在水力梯度一定的前提下,渗透率较大的样品,纵向弥散系数也较大;同时,结构的异质性也会影响溶质的穿透曲线。本研究提出的模拟方法可在土壤结构中进行高效的水流运动和溶质运移模拟,可用于土壤多孔介质在孔隙尺度下的水力学特性研究。  相似文献   

2.
孔隙尺度多孔介质流体流动与溶质运移高性能模拟   总被引:2,自引:0,他引:2       下载免费PDF全文
深入探究孔隙尺度下的流体流动特性和溶质运移规律对石油开采、农田养分管理、地下水污染修复有着重要意义。以人工构建的多孔介质结构和同步辐射X射线显微CT扫描的土壤团聚体(分辨率3.7 μm)为研究对象,在空间节点数多达64 000 000的情况下,基于格子Boltzmann模型和GPU并行技术计算得到多孔介质流体运动和溶质运移过程的关键参数,并据此探究多孔介质空间异质性对水力学特性的影响。通过对3组不同结构的多孔介质比较发现,结构复杂程度最高的土壤样品和不规则堆叠的圆球结构的渗透率在100 mD(即10-13m2)量级,远低于规则堆叠的圆球结构(>20 000 mD);土壤的迂曲度为1.40~1.60,明显高于规则堆叠的圆球结构。研究结果表明,渗透率大的样品具有较小的迂曲度,这与结构的空间异质性有较强的关系;土壤的渗透率和迂曲度呈现各向异性;在水力梯度一定的前提下,渗透率较大的样品,纵向弥散系数也较大;同时,结构的异质性也会影响溶质的穿透曲线。本研究提出的模拟方法可在土壤结构中进行高效的水流运动和溶质运移模拟,可用于土壤多孔介质在孔隙尺度下的水力学特性研究。  相似文献   

3.
多孔介质广泛存在于地下岩土结构、工程材料、生物体中,由于多孔介质微结构的复杂性,多孔介质渗流问题受到广泛关注。总结介绍了多孔介质渗流理论在国内外的研究现状,对多孔介质的描述、渗流问题的研究方法以及多孔介质中溶质运移、热量运移以及数值模型等方面的研究进展作了评述,并就多孔介质渗流理论在环境、农业、生态等领域的应用与发展方向提出看法。  相似文献   

4.
张建桥  窦智  张学羿 《地球科学》2023,(9):3444-3453
多孔介质污染物运移对于明晰地下水污染很重要,但在多孔介质中粒径变异系数(coefficientofvariation,COV)对内部微观孔隙结构中污染物运移过程影响的研究还存在不足.为此,基于随机算法,提出了一种不同COV且孔隙度一致的多孔介质几何模型构建方法,通过对Navier-Stokes(N-S)方程和对流-扩散方程(advection-diffusion equation,简称ADE)进行耦合求解得到多孔介质地下水流场及污染物浓度场,引入克里斯琴森均匀系数,定量评价流场流速分布的均匀性,基于MIM(mobile and immobile)模型和ADE模型分析耦合求解得到的穿透曲线特征.结果表明:随着粒径COV的增大,流场流速分布的不均性增强,MIM模型中的溶质流动区域占比β、无量纲传质率α*均增大;MIM模型拟合优度高于ADE模型,且随COV增大,ADE模型的拟合全局误差Ei逐步增大.总体上,粒径COV控制了溶质流动区域和非流动区域的大小及其之间的溶质交换强度,造成了多孔介质内部溶质运移的“非费克”行为,使得ADE模型的误差逐步增大,对于较大COV的多孔介...  相似文献   

5.
岩溶区地下水数值模拟研究进展   总被引:2,自引:2,他引:0  
岩溶含水介质的不均一性导致岩溶地下水流动、溶质运移和热量迁移的数学模拟研究成为地下水模拟的难点。本文综述了岩溶区地下水流模拟的几种方法,重点阐述了等效多孔介质法、双重连续介质法和三重介质法的定义、发展过程和适用范围,并回顾了这几种方法的研究成果。从等效多孔介质法到三重介质法,模拟精度不断提高,适用范围也逐渐由大区域实际问题向小区域理论研究过渡。介绍了溶质运移模拟和热迁移模拟的研究方法及实例。溶质运移模拟以对流弥散方程为基础,其中尺度效应是溶质运移模拟的重点研究问题;热量迁移模拟应考虑地下热水密度变化对地下热水运动的影响。溶质运移模拟和热量迁移模拟往往是将迁移模型和已经调试成功的地下水流动模型相耦合,从而达到模拟溶质及热量迁移的目的。由于溶质运移和热量迁移的复杂性,现阶段水流模型多数处于等效多孔介质模型阶段。综合理论及实际应用,指出精确刻画裂隙及管道和注重基础数学算法是岩溶水数值模拟进步的关键。   相似文献   

6.
李娜  任理 《水科学进展》2012,23(6):881-886
近年来,基于连续时间随机游动(Continuous Time Random Walk, CTRW)理论所建立的模拟非均质多孔介质中溶质运移的方法已在大量的数值实验、室内实验、野外实验中得到了广泛的验证,为非均质多孔介质中的溶质运移行为提供了一种有效的模拟方法。简述了提出和发展CTRW的研究背景、基础理论以及与经典的对流-弥散方程等其他模拟方法的关系,综述了该理论在模拟溶质运移中的发展和应用,分析了实际应用中的关键问题,并展望了将其进一步发展应用于模拟反应性溶质运移的前景。  相似文献   

7.
多尺度非均质多孔介质中溶质运移的蒙特卡罗模拟   总被引:4,自引:0,他引:4       下载免费PDF全文
探讨了将蒙特卡罗(Monte Carlo)方法应用于多尺度非均质含水层中溶质运移模拟的方法。所研究的含水层由两种具有不同渗透系数统计特征的多孔介质所组成,每一种多孔介质是非均质的,且其渗透系数场符合平稳假设,而整个模拟区的渗透系数是非平稳的。Monte Carlo方法要求参数是平稳的,因此,分别对两种多孔介质产生若干随机渗透系数场后,用两种方法进行组合,并进行溶质运移的模拟计算。通过对计算结果的分析,综合考虑计算精度、计算时间等因素,得出了处理多尺度非均质多孔介质中溶质运移问题的较好方法。  相似文献   

8.
刘咏  张琪  钱家忠  吴盾  张文永 《地学前缘》2022,29(3):248-255
采用图像分析法确定溶质浓度,以硫酸铜和EDTA二钠作为双分子反应物,在多孔介质模型中开展了不同粒径(1.52.0,2.53.0,3.54.0 mm)和流量(1.0,1.5,2.0 mL/s)下反应性溶质运移实验,探讨了应用不完全混合的对流弥散模型(IM-ADRE)对双分子反应溶质运移的模拟和预测,并进行了参数敏感性分析。结果表明:图像分析法可准确获取多孔介质中显色反应性溶质的浓度,灰度值与浓度的决定系数R2大于0.96;用IM-ADRE模型能够准确预测双分子反应性溶质硫酸铜和EDTA二钠在3种不同多孔介质中的运移过程,误差低于3.71%;实验条件的改变对IM-ADRE模型参数Dmβ0的影响显著,说明模型参数依赖于环境条件,其变化规律需要根据实际环境条件进一步率定,便于IM-ADRE模型的进一步推广应用。  相似文献   

9.
多孔介质的微观结构对可溶性污染物的运移具有重要影响。文章利用随机生长四参数生成法(Quartet Structure Generation Set, QSGS),构造具有各向异性孔隙结构的黏土多孔介质,且其固相颗粒具有随机性的粗糙表面。通过Stokes方程与对流扩散方程(ADE)耦合计算得到微观孔隙结构的非均匀水流场与浓度场。对黏土微观结构的非均匀水流场和浓度场进行分析,研究不同截面位置上水流分布情况,在脉冲注入的条件下讨论黏土微观结构对于可溶性污染物运移的影响。结果表明:具有随机性的黏土微观孔隙结构是导致水流场水流速度分布不均匀的主要因素;黏土中微小空隙或死端孔隙往往容易残留高浓度污染物,造成污染穿透曲线的“托尾”现象,也是污染后期的主要污染残留位置。  相似文献   

10.
李博  黄嘉伦  钟振  邹良超 《岩土力学》2019,40(9):3670-3768
对交叉裂隙渗流传质特性的定量描述是研究整个裂隙网络渗透传质特性的基础。为真实模拟水流及溶质在三维交叉裂隙中的运移过程,首先通过三维轮廓仪获取天然岩石裂隙表面的形貌数据,再应用三维重构技术生成相应的三维交叉裂隙模型,随后求解Navier-Stokes方程,假定溶质运移满足Fick定律,模拟水流和溶质在三维交叉裂隙中的运移过程。通过对比粗糙裂隙模型与平行平板模型的模拟结果发现:粗糙度对流体的分布及流动状态存在显著的影响;不同进、出口工况下的流体流动及溶质运移状态亦表明:裂隙交叉的几何形貌会显著地影响溶质混合行为。这些结果表明,目前被广泛采用的平行平板模型在评估岩体内特别是交叉口的物质运移特性时将导致较大的偏差,在将来的研究中有必要针对裂隙交叉口的几何特征建立修正的模型以提高评估的准确性。  相似文献   

11.
A study of the fracture distribution, hydraulic properties, groundwater levels and the transport of bromide was conducted to characterize vertical transport in the oxidized and reduced zones of a fractured glacial till. Detailed vertical profiles of groundwater levels and solute concentrations were obtained over a 4.5-year period. Vertical migration occurred at several time scales, as a low concentration front was rapidly transported at rates of 100–500 m/year ahead of a slower moving main plume, which advanced at rates of 0.2–0.8 m/year. Concentrations in the leading edge of the plume displayed a high degree of spatial variability over short vertical distances through day 1,000. Late in the test, the influence of matrix diffusion became apparent as concentration patterns developed from being irregular to more uniform distributions. Calculations show that the mass within the low concentration plume front accounts for less than 1% of the total solute mass. Simulation of the breakthrough curves using a simple one-dimensional advection-dispersion model of transport in porous media indicates that vertical transport is dominated by advection. Furthermore, the results indicate that vertical transport of solutes in oxidized and reduced zones of the till can be adequately simulated using an equivalent porous media.  相似文献   

12.
On the basis of site investigation and data collection of a certain electro plating factory, the groundwater flow and solute transport coupled models were established by applying Visual MODFLOW 4.1 software, which was used to conduct a numerical simulation that forecast the transport process of Cr6+ in groundwater. The results show that contamination plume of Cr6+ transports with groundwater flow direction. Without control measures, in 3 650 days, 19 wells for drinking would be contaminated, and the range of transport would be 52 172 m2, the maximum contamination would be 35.8 mg/L  相似文献   

13.
The combined influence of dip angle and adsorption heterogeneity on solute transport mechanisms in heterogeneous media can be understood by performing simulations of steady-state flow and transient transport in a heterogeneous aquifer with dipping anisotropy. Reactive and non-reactive contaminant transport in various types of heterogeneous aquifer is studied by simulations. The hydraulic conductivity (K) of the heterogeneous aquifer is generated by HYDRO_GEN with a Gaussian correlation spectrum. By considering the heterogeneity of the adsorption distribution coefficient (K d), a perfect negative correlation between lnK and lnK d is obtained by using the spherical grains model. The generated K and K d are used as input to groundwater flow and transport models to investigate the effects of dipping sedimentary heterogeneity on contaminant plume evolution. Simulation results showed that the magnitude of the dip angle strongly controls the plume evolution in the studied anisotropic and heterogeneous aquifer. The retarded average pore-water velocity (v/R) of the adsorption model significantly controls the horizontal spreading of the plume. The bottom plume is intensively retarded in the zones between the dipping lenses of lower hydraulic conductivity and the no-flow bottom boundary. The implications of these findings are very important for the management of contaminated heterogeneous aquifers.  相似文献   

14.
蒋立群  孙蓉琳  梁杏 《地球科学》2021,46(11):4150-4160
为探讨含水层非均质性不同刻画方法对地下水流和溶质运移预测的影响,基于非均质含水层砂箱实验,分别用传统等效均质模型、克立金插值和水力层析刻画含水层渗透系数场,并探讨了先验信息对水力层析结果的影响.将不同方法估算的渗透系数场用以预测地下水流和溶质运移过程,以此判断不同方法估算结果的优劣,分析含水层非均质性对地下水流和溶质运移的影响.结果表明:与克立金插值法相比,水力层析法可以更好地刻画含水层非均质性,较准确地预测地下水流和溶质运移过程;钻孔岩心渗透系数样本值作为先验信息可以提高水力层析法估算结果的精度;传统等效均质模型无法准确预测地下水流和溶质运移过程.含水层非均质性的增强将导致溶质污染羽分布形态和运移路径的空间变异性增强,并且优势通道直接决定溶质的分布及运移路径.   相似文献   

15.
溶质运移模拟是污染风险预测的一种有效手段。本研究为查明研究区内某有机污染场地作为污染源对拟建水库的影响,基于对区内流场的精确测量,建立了准确的地下水流模型和化工厂泄露有机污染物的溶质运移模型,预测分析了研究区拟建水库修建前、后地下水流场及污染物迁移情形的变化。模拟结果表明:水库和周边两种设计水头的截渗沟建成后,改变了自然条件下的流场,形成水库向截渗沟排泄、截渗沟向周边排泄以及水库和周边区域均向截渗沟排泄两种地下水流场。某些污染羽扩散范围虽比自然条件下有所增大,但向水库方向的扩散均未超出300 m,不会对水库水质构成污染风险。  相似文献   

16.
Sedimentological processes often result in complex three-dimensional subsurface heterogeneity of hydrogeological parameter values. Variogram-based stochastic approaches are often not able to describe heterogeneity in such complex geological environments. This work shows how multiple-point geostatistics can be applied in a realistic hydrogeological application to determine the impact of complex geological heterogeneity on groundwater flow and transport. The approach is applied to a real aquifer in Belgium that exhibits a complex sedimentary heterogeneity and anisotropy. A training image is constructed based on geological and hydrogeological field data. Multiple-point statistics are borrowed from this training image to simulate hydrofacies occurrence, while intrafacies permeability variability is simulated using conventional variogram-based geostatistical methods. The simulated hydraulic conductivity realizations are used as input to a groundwater flow and transport model to investigate the effect of small-scale sedimentary heterogeneity on contaminant plume migration. Results show that small-scale sedimentary heterogeneity has a significant effect on contaminant transport in the studied aquifer. The uncertainty on the spatial facies distribution and intrafacies hydraulic conductivity distribution results in a significant uncertainty on the calculated concentration distribution. Comparison with standard variogram-based techniques shows that multiple-point geostatistics allow better reproduction of irregularly shaped low-permeability clay drapes that influence solute transport.  相似文献   

17.
A recursion formulation for the transverse spreading of a solute is developed, and under conditions of steady flow in a stratified aquifer, the transport of a linearly sorbing solute undergoing nonequilibrium sorption is studied. The effect of spatial variability in the velocity field and the sorption kinetics are modeled to see the combined effect of the two processes on the spreading of the solute injected at a point in the aquifer. The main result of this work is a transport model based on a discrete formulation that includes local dispersion and leads to nonasymptotic behavior in the spreading of the plume in a direction normal to the mean flow velocity.  相似文献   

18.
Combining groundwater flow models with solute transport models represents a common challenge in groundwater resources assessments and contaminant transport modeling. Groundwater flow models are usually constructed at somewhat larger scales (involving a coarser discretization) to include natural boundary conditions. They are commonly calibrated using observed groundwater levels and flows (if available). The groundwater solute transport models may be constructed at a smaller scale with finer discretization than the flow models in order to accurately delineate the solute source and the modeled target, to capture any heterogeneity that may affect contaminant migration, and to minimize numerical dispersion while still maintaining a reasonable computing time. The solution that is explored here is based on defining a finer grid subdomain within a larger coarser domain. The local-grid refinement (LGR) implemented in the Modular 3D finite-difference ground-water flow model (MODFLOW) code has such a provision to simulate groundwater flow in two nested grids: a higher-resolution sub-grid within a coarse grid. Under the premise that the interface between both models was well defined, a comprehensive sensitivity and uncertainty analysis was performed whereby the effect of a parameter perturbation in a coarser-grid model on transport predictions using a higher-resolution grid was quantified. This approach was tested for a groundwater flow and solute transport analysis in support of a safety evaluation of the future Belgian near-surface radioactive waste disposal facility. Our reference coarse-grid groundwater flow model was coupled with a smaller fine sub-grid model in two different ways. While the reference flow model was calibrated using observed groundwater levels at a scale commensurate with that of the coarse-grid model, the fine sub-grid model was used to run a solute transport simulation quantifying concentrations in a hypothetical well nearby the disposal facility. When LGR coupling was compared to a one-way coupling, LGR was found to provide a smoother flow solution resulting in a more CPU-efficient transport solution. Parameter sensitivities performed with the groundwater flow model resulted in sensitivities at the head observation locations. These sensitivities identified the recharge as the most sensitive parameter, with the hydraulic conductivity of the upper aquifer as the second most sensitive parameter in regard to calculated groundwater heads. Based on one-percent sensitivity maps, the spatial distribution of the observations with the highest sensitivities is slightly different for the upper aquifer hydraulic conductivity than for recharge. Sensitivity analyses were further performed to assess the prediction scaled sensitivities for hypothetical contaminant concentrations using the combined groundwater flow and solute transport models. Including all pertinent parameters into the sensitivity analysis identified the hydraulic conductivity of the upper aquifer as the most sensitive parameter with regard to the prediction of contaminant concentrations.  相似文献   

19.
A finite element model was developed in order to solve for both regional groundwater flow and conservative solute transport in porous medium. The model was applied to a 55-sq-km groundwater basin in the Ruehen region of Germany using a network of 1450 elements and 780 nodes. This model was used in simulating a contaminant plume done through injection. Similarly, the model was applied to a 4750-sq-km portion in the eastern Nile Delta aquifer in Egypt. The model was applied to this portion of the delta using a network of 543 elements and 310 nodes with the main objective of simulating the problem of salt water intrusion.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号