共查询到20条相似文献，搜索用时 171 毫秒
1.
本文针对饱和多孔介质中渗流特点,在分析渗流连续介质力学模型的基础上,采用描述多孔介质中固相颗粒和液相流体耦合的细观力学模型,其中固相颗粒采用离散元的颗粒流理论(PFC)模拟,液相流体通过求解平均Navier-Stokes方程的计算流体动力学(CFD)技术计算,并利用建议的模型对二维渗流问题进行了模拟验证。数值模拟结果表明,采用PFC-CFD耦合细观力学模型能够描述从低雷诺流到高雷诺流范围很大的流体运动,可以用于多孔介质中的渗流问题分析。 相似文献
2.
:断裂多孔介质中的流体流动和组分输运过程对于理解金属成矿过程、油气运移、污染物质的输运和分布十分重要 ,近 2 0年来得到了广泛的研究。文中综述了这一领域的进展 ,涉及断裂多孔介质中的流体流动和组分输运的连续介质模型、多孔介质和断裂多孔介质的渗透率关系 ,特别是重点介绍了渗滤理论及其在断裂多孔介质中的流体流动和组分输运中的应用 相似文献
3.
4.
非常规油气资源的孔隙结构及其连通性非常复杂,其孔隙尺度从毫米到纳米跨越多个量级.多孔介质中气体的输运过程不仅依赖于介质的多尺度微观结构特征,还依赖于气体的相关属性.气体在多尺度多孔介质中的输运过程包括无滑流、滑脱流和过渡流,涉及分子扩散和努森扩散等多种机制,因此很难用唯一的连续介质理论来描述气体的输运特征.大量的数据表明真实多孔介质中的内部孔隙具有分形标度特征,因此采用分形几何表征多尺度多孔介质的孔隙结构,引入孔隙分形维数和迂曲度分形维数定量表征多孔介质的微结构和弯曲流道,建立多尺度多孔介质气体输运过程的细观模型;推导了多尺度多孔介质中气体的有效渗透率和有效扩散系数,并讨论了多尺度多孔介质微结构参数和气体属性对于气体等效输运特性的定量影响.该研究不仅可以丰富渗流理论,且有利于深入理解非常规油气藏的产出机制. 相似文献
5.
泥水盾构泥膜动态形成机制研究 总被引：2，自引：0，他引：2
随着大直径泥水盾构的广泛应用,盾构掘进过程切削成层土的可能性越来越大,常规的静力平衡的方法不能有效地控制开挖面的稳定。从泥膜形成的机制入手,将泥水仓中的泥浆视为由液相和颗粒相组成的两相流,将开挖地层视为多孔介质,利用多相流在多孔介质中的渗流特性提出了泥膜动态形成的过滤模型。此模型将泥膜视为控制体,并认为泥膜形成过程中遵守质量守恒和体积守恒定律,从边界条件上认为单位时间内从单元体内排出水量与土体积的压缩量相等。过滤模型较好地揭示了泥水盾构掘进过程中泥膜的动态形成过程和开挖地层力学性质的动态变化过程,对泥水盾构正面稳定有较为重要的意义。 相似文献
6.
7.
针对复杂流体运动中物质输运方程的数值求解面临地形复杂、数值阻尼过大以及数值振荡等难题,建立了Godunov格式下求解二维水流-输运方程的高精度耦合数学模型,提出了集成输运对流项的HLLC (Harten-Lax-van Leer-Contact)型近似黎曼算子,可同时计算水流通量及输运通量,不仅有效模拟了复杂地形上水流运动,而且解决了输运方程中对流项产生的数值阻尼过大和不稳定振荡等难题。采用水深-水位加权重构技术和Minmod限制器,提高了模型处理复杂混合流态的能力,同时结合Hancock预测-校正方法,使模型具有时空二阶精度。算例结果表明,模型精度高、稳定性好,能有效抑制数值阻尼,适合模拟实际复杂流体运动中物质的输运过程,具有较好的推广应用价值。 相似文献
8.
数学模型是研究相对渗透率与饱和度关系曲线的重要方法。采用自行开发设计的人工平面多孔介质模型,测定了相对渗透率与饱和度的关系曲线。多孔介质选择粒径为0.5~1mm、1~2mm的标准砂,纯净的水为湿润相,用3号苏丹红染色的93#汽油为非湿润相,组成多孔介质油水两相流动系统。采用Van Genuchten and Mualeum(VGM)和Brooks-Corey-Burdine(BCB)两种数学模型计算相对渗透率与饱和度的关系曲线,通过比较两种数学模型计算结果之间和模型计算结果与实测结果的差异以及模型的应用、多相渗流系统自身特征,得出VGM、BCB两种数学模型计算结果符合实际情况,VGM模型应用过程更为简便,但VGM模型具有一定适用条件;在砂性多孔介质中,BCB模型计算相对渗透率与饱和度关系曲线更准确。 相似文献
9.
计算流体地球化学研究的进展 总被引：3，自引：0，他引：3
成矿作用的化学机理可以通过实验和计算机模拟进行研究。随着计算机运算能力的不断增强 ,在地球化学中正在形成一门新兴学科———计算地球化学。其中热质输运模拟、化学质量迁移数值模拟和流体输运化学反应耦合动力学研究取得了显著进展。建立在Darcy定律和守恒方程基础上的多孔介质热质输运模拟通过流函数图、等温线图及速率矢量图等 ,从古水文学和流体地球化学方面高度动态研究成矿作用。根据化学和热力学原理进行的化学质量迁移数值模拟则通过矿物和流体中化学物种的热力学数据 ,预测多组分体系中发生的流体岩石相互作用 ,定量揭示经历了复杂化学反应进程的成矿作用的化学行为。将上述两方面结合的流体输运化学反应耦合动力学 ,可以从时间和空间上模拟真实成矿流体系统复杂的动力学行为 ,是计算流体地球化学的发展方向。 相似文献
10.
重非水相（DNAPLs）是地下水常见的有机污染物,理解其运移机制对于污染物修复具有重要意义。为解释实验中多孔介质润湿性对DNAPL的流动速度及残余饱和度的影响,使用多相流相场法,在构造的二维孔隙中模拟DNAPL液滴的下落过程。结果表明：模拟能准确刻画多相流界面的不稳定性,重现与实验室尺度相似的现象。文章提出的双界面模型相较于以往的微观模型能够更好的解释介质润湿性对DNAPL运移的影响。DNAPL的下落速度受前后界面共同影响,而饱和度主要由后界面决定。同时,润湿性对运移的影响大小受制于多孔介质的非均质性。 相似文献
11.
Gour-Tsyh Yeh Yilin Fang Fan Zhang Jiangtao Sun Yuan Li Ming-Hsu Li Malcolm D. Siegel 《Computational Geosciences》2010,14(1):149-170
Subsurface contamination problems of metals and radionuclides are ubiquitous. Metals and radionuclides may exist in the solute
phase or may be bound to soil particles and interstitial portions of the geologic matrix. Accurate tools to reliably predict
the migration and transformation of these metals and radionuclides in the subsurface environment enhance the ability of environmental
scientists, engineers, and decision makers to analyze their impact and to evaluate the efficacy of alternative remediation
techniques prior to incurring expense in the field. A mechanistic-based numerical model could provide such a tool. This paper
communicates the development and verification of a mechanistically coupled fluid-flow thermal-reactive biogeochemical-transport
model where both fast and slow reactions occur in porous and fractured media. Theoretical bases, numerical implementations,
and numerical experiments using the model are described. A definition of the “rates” of fast/equilibrium reactions is presented
to come up with a consistent set of governing equations. Two example problems are presented. The first one is a reactive transport
problem which elucidates the non-isothermal effects on heterogeneous reactions. It also demonstrates that the rates of fast/equilibrium
reactions are not necessarily greater than that of slow/kinetic reactions in the context of reactive transport. The second
example focuses on a complicated but realistic advective–dispersive–reactive transport problem. This example exemplifies the
need for innovative numerical algorithms to solve problems involving stiff geochemical reactions. It also demonstrates that
rates of all fast/equilibrium reactions are finite and definite. Furthermore, it is noted that a species-versus-time curve
cannot be used to characterize the rate of homogeneous fast/equilibrium reaction in a reactive transport system even if one
and only one such reaction is responsible for the production of this species. 相似文献
12.
Peter C. Lichtner 《Geochimica et cosmochimica acta》1985,49(3):779-800
A time-space continuum model for transport of hydrothermal fluids in porous media is presented which provides for simultaneous, reversible and irreversible chemical reactions involving liquids, gases and minerals. Homogeneous and heterogeneous reactions are incorporated in the model in a similar fashion through source/sink terms added to the continuity equation. The model provides for moving reaction fronts through surfaces of discontinuity across which occur jump discontinuities in the various field variables satisfying generalized Rankine-Hugoniot relations. Reversible reactions including aqueous complexing, oxidation-reduction reactions, mineral precipitation and dissolution reactions and adsorption are explicitly accounted for by imposing chemical equilibrium constraints in the form of mass action equations on the transport equations. This is facilitated by partitioning the reacting species into primary and secondary species corresponding to a particular representation of the stoichiometric reaction matrix referred to as the canonical representation. The transport equations for the primary species combined with homogeneous and heterogeneous equilibria result in a system of coupled, nonlinear algebraic/partial differential equations which completely describe the evolution of the system in time. Spatially separated phase assemblages are accommodated in the model by altering the set of independent variables across surfaces of discontinuity. Constitutive relations for the fluid flux corresponding to primary species are obtained describing transport of both neutral and charged species by advection, dispersion and diffusion. Numerical implementation of the transport equations is considered and both explicit and implicit finite difference algorithms are discussed. Analytical expressions for the change in porosity and permeability with time are obtained for an assemblage of minerals reacting reversibly with a hydrothermal fluid under quasi-steady state conditions. Fluid flow is described by Darcy's law employing a phenomenological expression relating permeability and porosity. Finally an expression for the local retardation factor of solute species is derived for the case of advective transport in a single spatial dimension which accounts for the effects of homogeneous and heterogeneous equilibria including adsorption on the rate of advance of a reaction front. The condition for the formation of shock waves is given. 相似文献
13.
对饱和多孔介质提出了一个含溶混污染物输运（传质）过程的混合元方法,其中污染物输运过程数学模型包含了对流、机械逸散、分子弥散和吸附等机制。固相位移、应变和有效应力,孔隙水压力、压力空间梯度和Darcy速度,污染物浓度、浓度空间梯度和浓度流量在单元内均为独立变量分别插值。基于胡海昌-Washizu三变量广义变分原理,结合可以滤掉虚假振荡的特征线方法,推导出饱和土中水力-力学-传质耦合问题控制方程的单元弱形式,并导出了混合元计算公式。数值模拟证明了所提出的方法可以提供与传统4点积分方案同样精度,同时能够提高计算效率。 相似文献
14.
A non‐equilibrium sorption—advection—diffusion model to simulate miscible pollutant transport in saturated–unsaturated soils is presented. The governing phenomena modelled in the present simulation are: convection, molecular diffusion, mechanical dispersion, sorption, immobile water effect and degradation, including both physical and chemical non‐equilibrium processes. A finite element procedure, based on the characteristic Galerkin method with an implicit algorithm is developed to numerically solve the model equations. The implicit algorithm is formulated by means of a combination of both the precise and the traditional numerical integration procedures. The stability analysis of the algorithm shows that the unconditional stability of the present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results illustrate good performance of the present algorithm in stability and accuracy, and in simulating the effects of all the mentioned phenomena governing the contaminant transport and the concentration distribution. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
15.
This paper presents a numerical model for simulating free surface flow in porous media with spatially varying porosity. The governing equations are based on the mixture theory. The resistance forces between solid and fluid is assumed to be nonlinear. A multiphase SPH approach is presented to solve the governing equations. In the multiphase SPH, water is modeled as a weakly compressible fluid, and solid phase is discretized by fixed solid particles carrying information of porosity. The model is validated by several numerical examples including seepage through specimen, fast flow through rockfill dam and wave interaction with porous structure. Good agreements between numerical results and experimental data are obtained in terms of flow rate and evolution of free surface. Parameter study shows that (1) the nonlinear resistance law provides more accurate results; (2) particle size and porosity have significant influence on the porous flow. 相似文献
16.
The transport of chemically reactive solutes (e.g. surfactants, CO2 or dissolved minerals) is of fundamental importance to a wide range of applications in oil and gas reservoirs such as enhanced
oil recovery and mineral scale formation. In this work, we investigate exponential time integrators, in conjunction with an
upwind weighted finite volume discretisation in space, for the efficient and accurate simulation of advection–dispersion processes
including non-linear chemical reactions in highly heterogeneous 3D oil reservoirs. We model sub-grid fluctuations in transport
velocities and uncertainty in the reaction term by writing the advection–dispersion–reaction equation as a stochastic partial
differential equation with multiplicative noise. The exponential integrators are based on the variation of constants solution
and solve the linear system exactly. While this is at the expense of computing the exponential of the stiff matrix representing
the finite volume discretisation, the use of real Léja point or the Krylov subspace technique to approximate the exponential
makes these methods competitive compared to standard finite difference-based time integrators. For the deterministic system,
we investigate two exponential time integrators, the second-order accurate exponential Euler midpoint (EEM) scheme and exponential
time differencing of order one (ETD1). All our numerical examples demonstrate that our methods can compete in terms of efficiency
and accuracy compared with standard first-order semi-implicit time integrators when solving (stochastic) partial differential
equations that model mixing and chemical reactions in 3D heterogeneous porous media. Our results suggest that exponential
time integrators such as the ETD1 and EEM schemes could be applied to typical 3D reservoir models comprising tens to hundreds
of thousands unknowns. 相似文献
17.
A three-dimensional, three-phase numerical model is presented for simulating the movement of immiscible fluids, including
nonaqueous-phase liquids (NAPLs), through porous media. The model is designed to simulate soil flume experiments and for practical
application to a wide variety of contamination scenarios involving light or dense NAPLs in heterogeneous subsurface systems.
The model is derived for the three-phase flow of water, NAPL, and air in porous media. The basic governing equations are based
upon the mass conservation of the constitutents within the phases. The descretization chosen to transform the governing equations
into the approximating equations, although logically regular, is very general. The approximating equations are a set of simultaneous
coupled nonlinear equations which are solved by the Newton-Raphson method. The linear system solutions needed for the Newton-Raphson
method are obtained using a matrix of preconditioner/accelerator iterative methods.
Because of the special way the governing equations are implemented, the model is capable of simulating many of the phenomena
considered necessary for the sucessful simulation of field problems including entry pressure phenomena, entrapment, and preferential
flow paths. The model is verified by comparing it with several exact analytic test solutions and three soil flume experiments
involving the introduction and movement of light nonaqueous-phase liquid (LNAPL) or dense nonaqueous-phase liquid (DNAPL)
in heterogeneous sand containing a watertable.
This revised version was published online in August 2006 with corrections to the Cover Date. 相似文献
18.
The surface subsidence above a compacting saturated oil reservoir is the main topic of this paper. From a literature review, it is obvious that extensive efforts have been conducted for investigating this phenomenon in various situations. Herein, a numerical model, based on the finite element method, was used for simulating three-dimensional three-phase fluid flow in a deforming saturated oil reservoir. The mathematical formulation describes a fully coupled governing equation system which consists of the equilibrium and continuity equations for three immiscible fluids flowing in a porous media. An elastoplastic soil model, based on a Mohr Coulomb yield surface, was utilized. The finite element method was applied to obtain simultaneous solutions to the governing equations where the displacements and the fluid pressures are the primary unknowns. The final discretized equations are solved by a direct solver using fully implicit procedures. A linear analysis was used to study the stability conditions of the present model. Finally, a series of simulations were conducted to indicate the validity and the utility of the developed model. 相似文献
19.
We present a method to transform the governing equations of multispecies reactive transport in porous media. The reformulation leads to a smaller problem size by decoupling of equations and by elimination of unknowns, which increases the efficiency of numerical simulations. The reformulation presented here is a generalization of earlier works. In fact, a whole class of transformations is now presented. This class is parametrized by the choice of certain transformation matrices. For specific choices, some known formulations of reactive transport can be retrieved. Hence, the software based on the presented transformation can be used to obtain efficiency comparisons of different solution approaches. For our efficiency tests, we use the MoMaS benchmark problem on reactive transport. 相似文献
20.
A numerical simulation is presented for three-dimensional three-phase fluid flow in a deforming saturated oil reservoir. The mathematical formulation describes a fully coupled governing equation systen which consists of the equilibrium and continuity equations for three immiscible fluids flowing in a porous medium. An elastoplastic soil model, based on a Mohr–Coulomb yield surface, is used. The finite element method is applied to obtain simultaneous solutions to the governing equations where displacement and fluid pressures are the primary unknowns. The final discretized equations are solved by a direct solver using fully implicit procedures. The developed model is used to investigate the problems of three-phase fluid flow in a deforming saturated oil reservoir. 相似文献