共查询到20条相似文献，搜索用时 419 毫秒
1.
通过建立二维物理模型,研究非饱和层状非均质多孔介质中轻非水相液体LNAPLs(Light Non-a-queous Phase Liquids)的入渗机制与变化特征,建立LNAPLs运移与分布的锋面扩展模型,探讨入渗阶段油流锋面扩展速率的变化规律。结果发现:模型中砂土的渗透系数、孔隙度及油的相对渗透系数对LNAPLs运移的锋面扩展速率影响较大,而油的饱和度与压头变化对其影响较小。当LNAPLs由中砂进入细砂与细砂倾斜透镜体时,砂的渗透系数、孔隙度及油的相对渗透系数减小,这种介质结构面的突变改变了LNAPLs锋面扩展速率,其在中砂中锋面扩展速率快而在细砂与细砂倾斜透镜体中锋面扩展速率慢。此外,当LNAPLs进入干湿界面后,由于毛细作用的增强与含水量的加大,其横向比垂向运移的平均锋面扩展速率大。 相似文献
2.
将二氧化碳注入到深部咸水层中,形成复杂的多组分、多相流系统。二氧化碳在压力梯度、浓度差作用下不断扩散,逐步带走盐溶液中的水分,导致各组分的相态变化,盐结晶析出,阻塞了咸水层孔隙通道,从而降低了二氧化碳的注入效率,研究该干化效应的影响因素并为工程选址提供依据具有重要意义。采用二维径向模型建立多相流体的流动方程,并结合相对渗透率和毛细压力方程探讨二氧化碳注入速率、咸水层含盐量、毛细压力的特征参数对干化效应的影响,干化效应可用固体饱和度值进行定量描述。结果表明：二氧化碳运移分3个区域：干涸区、气液相混合区及液相咸水区,干化效应主要发生在井周的干涸区。在毛细作用下固体饱和度随注入速率的减小而增大,随咸水层含盐量增大而增大,随毛细作用增大而增大。因此,提高二氧化碳的注入速率,向咸水层中注水稀释含盐量或选择粒径较大的均质咸水层减小毛细作用,均可降低盐结晶对孔隙通道的阻塞,提高注入效率。 相似文献
3.
In this paper, a numerical model is developed for the fully coupled hydro‐mechanical analysis of deformable, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non‐wetting pore fluids, in which the coupling between various processes is taken into account. The governing equations involving the coupled solid skeleton deformation and two‐phase fluid flow in partially saturated porous media including cohesive cracks are derived within the framework of the generalized Biot theory. The fluid flow within the crack is simulated using the Darcy law in which the permeability variation with porosity because of the cracking of the solid skeleton is accounted. The cohesive crack model is integrated into the numerical modeling by means of which the nonlinear fracture processes occurring along the fracture process zone are simulated. The solid phase displacement, the wetting phase pressure and the capillary pressure are taken as the primary variables of the three‐phase formulation. The other variables are incorporated into the model via the experimentally determined functions, which specify the relationship between the hydraulic properties of the fracturing porous medium, that is saturation, permeability and capillary pressure. The spatial discretization is implemented by employing the extended finite element method, and the time domain discretization is performed using the generalized Newmark scheme to derive the final system of fully coupled nonlinear equations of the hydro‐mechanical problem. It is illustrated that by allowing for the interaction between various processes, that is the solid skeleton deformation, the wetting and the non‐wetting pore fluid flow and the cohesive crack propagation, the effect of the presence of the geomechanical discontinuity can be completely captured. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
4.
When fluid flows in porous media under subsurface conditions, significant deformation can occur. Such deformation is dependent on structural and phase characteristics. In this paper, we investigate the effect of multiphase flow on the deformation of porous media at the pore scale by implementing a strongly coupled partitioned solver discretized with finite volume (FV) technique. Specifically, the role of capillary forces on grain deformation in porous media is investigated. The fluid and solid subdomains are meshed using unstructured independent grids. The model is applied for solving multiphase coupled equations and is capable of capturing pore scale physics during primary drainage by solving the Navier-Stokes equation and advecting fluid indicator function using volume of fluid (VOF) while the fluid is interacting with a nonlinear elastic solid matrix. The convergence of the coupled solver is accelerated by Aitken underrelaxation. We also reproduce geomechanical stress conditions, at the pore scale, by applying uniaxial stress on the solid while simultaneously solving the multiphase fluid-solid interaction problem to investigate the effect of external stress on fluid occupancy, velocity-field distribution, and relative permeability. We observe that the solid matrix exhibits elasto-capillary behavior during the drainage sequence. Relative permeability endpoints are shifted on the basis of the external stress exerted. 相似文献
5.
《Physics and Chemistry of the Earth, Part A: Solid Earth and Geodesy》2000,25(2):169-174
A simple cutoff approach based on the capillary bundle model has become an industrial standard method to quickly obtain the irreducible water saturation from NMR T2 distribution. However this approach is not always valid. To overcome the shortcoming of the capillary bundle model that ignores pore-pore connectivity, we have conducted a two-phase flow simulation in a CT-based pore network. The CT -based pore network is a representation of a real rock pore structure that is described by a binary X-ray tomographic data set. Simulation on such network mimics a process of the porous plate measurement. The generated capillary curve is quite reasonable. The oil accession distributions at different water saturations plotted as a function of pore size provide an insight for the immiscible displacement process in the real rock pore structure. Water is trapped not only in dead ends but also in the well-connected pores due to a pore level by-pass mechanism. At the capillary end point pressure, a plot of the trapped water distribution as a function of pore size has the same lognormal distribution as the pore size distribution, which is much different from what the simple capillary bundle model suggests. 相似文献
6.
Steam drive recovery of oil is an economical way of producing oil even in times of low oil prices and is used worldwide. This paper focuses on the one-dimensional setting, where steam is injected into a core initially containing oil and connate water while oil and water are produced at the other end. A three-phase (oil, water, steam) hot zone develops, which is abruptly separated from the two-phase (oil + water) cold zone by the steam condensation front. The oil, water and energy balance equations (Rankine–Hugoniot conditions) cannot uniquely solve the system of equations at the steam condensation front. In a previous study, we showed that two additional constraints follow from an analysis of the traveling wave equation representing the shock; however, within the shock, we assumed local thermodynamic equilibrium. Here we extend the previous study and include finite condensation rates; using that appropriate scaling requires that the Peclet number and the Damkohler number are of the same order of magnitude. We give a numerical proof, using a color-coding technique, that, given the capillary diffusion behavior and the rate equation, a unique solution can be obtained. It is proven analytically that the solution for large condensation rates tends to the solution obtained assuming local thermodynamic equilibrium. Computations with realistic values to describe the viscous and capillary effects show that the condensation rate can have a significant effect on the global saturation profile, e.g. the oil saturation just upstream of the steam condensation front. 相似文献
7.
In this paper, a fully coupled numerical model is presented for the finite element analysis of the deforming porous medium interacting with the flow of two immiscible compressible wetting and non-wetting pore fluids. The governing equations involving coupled fluid flow and deformation processes in unsaturated soils are derived within the framework of the generalized Biot theory. The displacements of the solid phase, the pressure of the wetting phase and the capillary pressure are taken as the primary unknowns of the present formulation. The other variables are incorporated into the model using the experimentally determined functions that define the relationship between the hydraulic properties of the porous medium, i.e. saturation, relative permeability and capillary pressure. It is worth mentioning that the imposition of various boundary conditions is feasible notwithstanding the choice of the primary variables. The modified Pastor–Zienkiewicz generalized constitutive model is introduced into the mathematical formulation to simulate the mechanical behavior of the unsaturated soil. The accuracy of the proposed mathematical model for analyzing coupled fluid flows in porous media is verified by the resolution of several numerical examples for which previous solutions are known. Finally, the performance of the computational algorithm in modeling of large-scale porous media problems including the large elasto-plastic deformations is demonstrated through the fully coupled analysis of the failure of two earth and rockfill dams. Furthermore, the three-phase model is compared to its simplified one which simulates the unsaturated porous medium as a two-phase one with static air phase. The paper illustrates the shortcomings of the commonly used simplified approach in the context of seismic analysis of two earth and rockfill dams. It is shown that accounting the pore air as an independent phase significantly influences the unsaturated soil behavior. 相似文献
8.
多孔介质中的流动问题,与孔隙介质的特性,含水量状态以及含水量的变化历史密切相关。基于毛细循环滞回理论模型,考虑含水量变化历史对土水特征关系的影响,在开发的U-DYSAC2有限元程序中进行了相应的数值实施。在试验给定的初边值条件下进行了非饱和渗流模拟分析,并将模拟结果与实测数据比较,表明在压力边界条件反复变化下,考虑滞回效应能获得更接近实测的结果,证实该模型在模拟各种循环变化条件下非饱和土渗流初边值问题的适用性与必要性。对入渗重分布反复变化条件下非饱和土柱流动的数值模拟表明,考虑滞回与不考虑滞回条件下,含水量、孔隙水压力和湿峰的迁移的预测在入渗后的重分布过程差异较大。考虑滞回效应时,土柱上部的脱湿速率、下部的吸湿速率比不考虑滞回时要低。从而证实了非饱和多孔介质中的土水状态依赖于含水量变化,而且强烈依赖于土体的水力路径变化。因此,循环边界条件变化下,毛细滞回效应在非饱和渗流模拟中的影响显著,必须加以考虑。 相似文献
9.
Air sparging (AS) is an in situ soil/groundwater remediation technology, which involves the injection of pressurized air/oxygen through an air sparging well below the zone of contamination. Characterizing the mechanisms governing movement of air through saturated porous media is critical for the design of an effective cleanup treatment system. In this research, micromechanical investigation was performed to understand the physics of air migration and subsequent spatial distribution of air at pore scale during air sparging. The void space in the porous medium was first characterized by pore network consisting of connected pore bodies and bonds. The biconical abscissa asymmetric concentric bond was used to describe the connection between two adjacent pore bodies. Then a rule‐based dynamic two‐phase flow model was developed and applied to the pore network model. A forward integration of time was performed using the Euler scheme. For each time step, the effective viscosity of the fluid was calculated based on fractions of two phases in each bond, and capillary pressures across the menisci was considered to compute the pressure field. The developed dynamic model was used to study the rate‐dependent drainage during air sparging. The effect of the capillary number and geometrical properties of the network on the dynamic flow properties of two‐phase flow including residual saturation, spatial distribution of air and water, dynamic phase transitions, and relative permeability‐capillary pressure curves were systematically investigated. Results showed that all the above information for describing the air water two‐phase flow are not intrinsic properties of the porous medium but are affected by the two‐phase flow dynamics and spatial distribution of each phase, providing new insight to air sparging. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
10.
In this paper we study a one-dimensional model for oil recovery by steamdrive. This model consists of two parts: a (global)
interface model and a (local) steam condensation/capillary diffusion model. In the interface model a steam condensation front
(SCF) is present as an internal boundary between the hot steam zone (containing water, oil and steam) and the cold liquid
zone (containing only water and oil). Disregarding capillary pressure away from the SCF, a 2× 2 hyperbolic system arises for
the water and steam saturation. This system cannot be solved uniquely without additional conditions at the SCF. To find such
conditions we blow up the SCF and consider a parabolic transition model, including capillary diffusion. We study in detail
the existence conditions for traveling wave solutions. These conditions provide the missing matching conditions at the SCF
in the hyperbolic limit. We show that different transition models yield different matching conditions, and thus different
solutions of the interface model. We also give a relatively straightforward approximation and investigate its validity for
certain ranges of model parameters.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
11.
One of the driving forces in porous media flow is the capillary pressure. In standard models, it is given depending on the
saturation. However, recent experiments have shown disagreement between measurements and numerical solutions using such simple
models. Hence, we consider in this paper two extensions to standard capillary pressure relationships. Firstly, to correct
the nonphysical behavior, we use a recently established saturation-dependent retardation term. Secondly, in the case of heterogeneous
porous media, we apply a model with a capillary threshold pressure that controls the penetration process. Mathematically,
we rewrite this model as inequality constraint at the interfaces, which allows discontinuities in the saturation and pressure.
For the standard model, often finite-volume schemes resulting in a nonlinear system for the saturation are applied. To handle
the enhanced model at the interfaces correctly, we apply a mortar discretization method on nonmatching meshes. Introducing
the flux as a new variable allows us to solve the inequality constraint efficiently. This method can be applied to both the
standard and the enhanced capillary model. As nonlinear solver, we use an active set strategy combined with a Newton method.
Several numerical examples demonstrate the efficiency and flexibility of the new algorithm in 2D and 3D and show the influence
of the retardation term.
This work was supported in part by IRTG NUPUS. 相似文献
12.
蒸汽注射法是一种新近发展起来并对包气带非水相液体（NAPLs）污染修复具有较好应用前景的修复方法。为了探讨热蒸汽在非饱和多孔介质中的运移特征和蒸汽注射技术修复氯苯污染土壤的效果,本文通过一维模拟柱分别开展了热蒸汽在多孔介质中的迁移规律实验和蒸汽注射修复氯苯实验。研究表明：模拟柱中给定点处温度随时间变化可分为环境温度段、升温段和饱和蒸汽温度段3个阶段,随着温度锋面向下迁移,其迁移速度逐渐减小;模拟柱中的压力分布在蒸汽覆盖区域且空间分布呈线性规律,给定点处的压力随时间先增加后趋于平稳;实验结束后测得的水饱和度随蒸汽迁移距离的增加逐渐增加。在注入蒸汽质量流量为0.3 kg/h和氯苯初始质量分数为56.8 mg/kg的条件下,热蒸汽修复细砂中氯苯3.5 h后去除率达98.0%。模拟柱底部因蒸汽冷凝水达到饱和而产生出流液,出流10 min时出流液中氯苯最大质量浓度达到152.98 mg/L,70 min后出流液中氯苯质量浓度低于7.00 mg/L;介质中氯苯质量分数随距离的增加而增加,最大残留量为0.36 mg/kg。蒸汽注射法有效地改善了SVE（土壤气相抽提）法"拖尾"现象。 相似文献
13.
M D Sharma 《Journal of Earth System Science》2018,127(1):7
Propagation of surface waves is discussed in a cylindrical borehole through a liquid-saturated porous solid of infinite extent. The porous medium is assumed to be a continuum consisting of a solid skeletal with connected void space occupied by a mixture of two immiscible inviscid fluids. This model also represents the partial saturation when liquid fills only a part of the pore space and gas bubbles span the remaining void space. In this isotropic medium, potential functions identify the existence of three dilatational waves coupled with a shear wave. For propagation of plane harmonic waves along the axially-symmetric borehole, these potentials decay into the porous medium. Boundary conditions are chosen to disallow the discharge of liquid into the borehole through its impervious porous walls. A dispersion equation is derived for the propagation of surface waves along the curved walls of no-liquid (all gas) borehole. A numerical example is studied to explore the existence of cylindrical waves in a particular model of the porous sandstone. True surface waves do not propagate along the walls of borehole when the supporting medium is partially saturated. Such waves propagate only beyond a certain frequency when the medium is fully-saturated porous or an elastic one. Dispersion in the velocity of pseudo surface waves is analysed through the changes in consolidation, saturation degree, capillary pressure or porosity. 相似文献
14.
水力特性是非饱和土力学理论与工程的重要课题之一,土样水力特性的变化过程本质上可反映为土颗粒间毛细力的演化规律。为此,以粗颗粒为研究对象,将其简化为一对不等径球体颗粒,而其间的水分形态可视为形如圆环的液桥,不考虑颗粒重力和浮力的影响。以Young-Laplace方程为基础,先推得计算液桥毛细力的控制方程组,再结合液桥的无量纲体积最大、最小值及其外曲率半径的割线迭代算法提出求解毛细力的动态计算方法,进而研究毛细力与颗粒间距、颗粒半径比和液桥体积的无量纲关系,结果表明：当颗粒间距一定时,液桥的毛细力随其体积和颗粒半径比增大均呈递增趋势;当液桥体积一定时,其抗拉刚度随颗粒间距增大呈递减趋势。最后,利用已有文献中液桥毛细力与颗粒间距的实测关系,验证了该动态计算方法在表征液桥从形成至断裂时毛细力演化规律的有效性。 相似文献
15.
冰点下水合物在多孔介质的生成是一个复杂的多相转化过程,为了研究冰点下多孔介质中水合物生成过程的水相转化率、气体消耗量与稳定压力等生成特性,在定容条件下,进行不同孔径与粒径的多孔介质中甲烷水合物在冰点下的生成实验。所使用的多孔介质平均孔径为1295 nm、1796 nm和3320 nm。研究结果表明：水合物生成结束时水的转化率随着初始生成压力的增大而增大,随着温度的升高而降低,随着孔径的增大而增大;在相同的孔径下,多孔介质粒径的增大降低了水合物的生成速率但对最终气体消耗量没有影响;在相同的温度下,随着初始生成压力的增大,实验最终压力、气体消耗量与最终水的转化率均随之升高;温度越高,不同的生成初始生成压力下体系的最终稳定压力与水的转化率相差越大;在多孔介质的毛细管作用力与结合水的共同作用下,冰点下水合物生成的水的转化率会大大地降低。在本实验条件下,水相转化为水合物的比例最高为32.39%。 相似文献
16.
饱和多孔介质的动力响应研究在众多工程领域具有重要意义。充分考虑孔隙率的变化规律与影响因素,有利于合理揭示饱和多孔介质的相关力学行为。为此,将动态孔隙率模型与用于表征饱和多孔介质动力特性的u-U-p型方程结合,构建相应的非线性力学模型,利用Comsol Multiphysis PDE求取相应的数值解,以此研究不同透水条件下,受谐波载荷激励的二维饱和土体的孔隙率、变形量及孔隙水压力的变化规律。结果表明：孔隙率的变化与土骨架的体应变及孔隙水压力直接相关,土体压缩过程中,孔隙率相应减小,土骨架与孔隙流体的相互作用增强,土体运动时所受阻力增大,其无量纲竖向位移小于孔隙率被视为常数时的情况,在此条件下,由于土体的变形量减小,其孔隙水压力也相对减小。故充分考虑动态孔隙率,有利于更加精确地研究等饱土体和多孔介质的相关力学行为。此外,土体上表面透水条件下,孔隙流体可以从土体表面自由排出,土骨架承受的载荷更大,与不透水条件相比,土体孔隙率、竖向位移、孔隙水压力等变化更为显著。 相似文献
17.
We present a network flow model to compute transport, through a pore network, of a compositional fluid consisting of water with a dissolved hydrocarbon gas. The model captures single-phase flow (below local bubble point conditions) as well as the genesis and migration of the gas phase when bubble point conditions are achieved locally. Constant temperature computational tests were run on simulated 2D and 3D micro-networks near bubble point pressure conditions. In the 2D simulations which employed a homogeneous network, negligible capillary pressure, and linear relative permeability relations, the observed concentration of CO2 dissolved in the liquid phase throughout the medium was linearly related to the liquid pressure. In the case of no gravity, the saturation of the gas phase throughout the medium was also linearly related to the liquid pressure; under gravity, the relationship became nonlinear in regions where buoyancy forces were significant. The 3D heterogeneous network model had nonnegligible capillary pressure and nonlinear relative permeability functions. While 100 % of the CO2 entered the 3D network dissolved in the liquid phase, 25 % of the void space was occupied by gas phase and 47 % of the CO2 exiting the outlet face did so via the gaseous phase after 500 s of simulation time. 相似文献
18.
Chongbin Zhao B. E. Hobbs P. Hornby A. Ord Shenglin Peng Liangming Liu 《国际地质力学数值与分析法杂志》2008,32(9):1107-1130
The chemical‐dissolution front propagation problem exists ubiquitously in many scientific and engineering fields. To solve this problem, it is necessary to deal with a coupled system between porosity, pore‐fluid pressure and reactive chemical‐species transport in fluid‐saturated porous media. Because there was confusion between the average linear velocity and the Darcy velocity in the previous study, the governing equations and related solutions of the problem are re‐derived to correct this confusion in this paper. Owing to the morphological instability of a chemical‐dissolution front, a numerical procedure, which is a combination of the finite element and finite difference methods, is also proposed to solve this problem. In order to verify the proposed numerical procedure, a set of analytical solutions has been derived for a benchmark problem under a special condition where the ratio of the equilibrium concentration to the solid molar density of the concerned chemical species is very small. Not only can the derived analytical solutions be used to verify any numerical method before it is used to solve this kind of chemical‐dissolution front propagation problem but they can also be used to understand the fundamental mechanisms behind the morphological instability of a chemical‐dissolution front during its propagation within fluid‐saturated porous media. The related numerical examples have demonstrated the usefulness and applicability of the proposed numerical procedure for dealing with the chemical‐dissolution front instability problem within a fluid‐saturated porous medium. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
19.
冻土离心模型试验相似准则分析 总被引：1，自引：0，他引：1
在研究土壤冻融问题时,土工离心模型试验在时间效应和应力全等型模拟方面更具优势。然而,目前关于冻土离心模型试验的研究成果较少,其中的相似准则也不够完善。针对这一现状,在充分考虑土体的水-热-力耦合作用过程的基础上,采用Butterfield量纲分析法确定了控制饱和土冻融变形性状的无量纲项,建立了冻土离心模型试验中孔隙压力、热扩散、未冻水迁移、融土固结以及冻融变形的相似准则。分析结果表明,冻土离心模型试验中未冻水迁移、融土固结及热扩散效应的非稳态时间具有统一的比尺,即离心模型是原型的1/N2倍,而未冻水迁移流速的比尺为离心模型是原型的N倍。此外,利用冻土离心模拟装置完成了一则对渠道基土冻融作用的离心模拟。 相似文献
20.
Jianhui Yang Igor Bondino Mohamed Regaieg Arthur Moncorgé 《Computational Geosciences》2017,21(5-6):849-862
Pore network modelling (PNM) has been widely used to study the multiphase flow and transport in porous media. Although a number of recent papers discussed the PNM validation on core-scale parameters such as permeability, relative permeability and capillary pressure; quantitative predictive potential of PNM on pore by pore basis has rarely been studied. The aim of this paper is to present a direct comparison between PNM simulations and corresponding micro-model experiments at the same scale and the same geometry. A number of well-defined and constrained two-phase flow in porous medium experimental scenarios were utilized to validate the physics solving part in PNM (filling rules, capillary and viscous pressure). This work validates that a dynamic pore network flow solver can predict two-phase flow displacements for these experiments for drainage situations at both pore and plug scales. A glass-etched micro-model is used to quantify the accuracy of a dynamic PNM solver on pore and core levels. Two-phase drainage micro fluidic experiments at different flow conditions are performed on micro-models. PNM simulations are performed on the same pattern and flow conditions as used in micro-model experiments. The two-phase distribution extracted from experiment images is registered onto rsults of PNM simulations for direct pore to pore comparison. Pore-scale matching level is found at around 75 % for all three test cases. The matching level of core-scale parameters such as S w c and oil-phase permeability varies from case to case; the relative error to micro-model experiment measurements varies from 15 to 60 %. Possible reasons leading to discrepancies on core-scale parameters are discussed: missing considerations during validation of the combination of uncertainty in both simulator input parameters and experiments are seen as the principal factors. 相似文献