共查询到20条相似文献，搜索用时 15 毫秒
1.
Fatiha Bouchelaghem 《Mathematical Geosciences》2011,43(1):23-53
This paper presents an analytical study of fluid flow in a porous medium presenting pores of two different length scales: at the smallest or microscopic scale, the presence of connected voids confers a porous medium structure to the material investigated, while at the upper or mesoscopic scale, macro-pores are present. This microstructure is employed to represent the progressive opening of inter-aggregate pore spaces observed in natural compacted montmorillonites polluted by heavy metal ions. Three-dimensional analytical expressions are rigorously derived for pore fluid velocity and excess pore fluid pressure within the porous matrix, around an occluded ellipsoidal inter-aggregate void. The eccentricity ratio is employed to characterize the geometrical shape of the ellipsoidal void, while its orientation with respect to the inflow in the far field is determined by the dip angle θ. As an application, we investigate the flow focusing effect for varying eccentricity ratios and dip angles. 相似文献
2.
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. 相似文献
3.
4.
A three-phase hydro-mechanical model for hydraulic fracturing is proposed. Three phases include: porous solid, fracturing fluid and host fluid. Discontinuity is handled using extended finite element method (XFEM) while cohesive crack model is used as fracturing criterion. Flow through fracture is defined as one-dimensional laminar flow, and flow through porous medium (host reservoir) is defined as two-dimensional Darcy flow. Coupling between two fluids in each space, fracture and pore, is captured through capillary pressure–saturation relationship, while the identical fluids in fracture and pore are coupled through a so-called leak-off mass transfer term. Coupling between fluids and deformation is captured through compatibility of volumetric strain of fluids within fracture and pore, and volumetric strain of the matrix. Spatial and temporal discretisation is achieved using the standard Galerkin method and the finite difference technique, respectively. The model is verified against analytical solutions available from literature. The leaking of fracturing fluid into the medium and suction of porous fluid into the fracture around the tip, are investigated. Sensitivity analyses are carried out for cases with slow and fast injection rates. It is shown that the results by single-phase flow may underestimate the leak-off. 相似文献
5.
Frictionless contact formulation for dynamic analysis of nonlinear saturated porous media based on the mortar method 
下载免费PDF全文
下载免费PDF全文 A finite element algorithm for frictionless contact problems in a two‐phase saturated porous medium, considering finite deformation and inertia effects, has been formulated and implemented in a finite element programme. The mechanical behaviour of the saturated porous medium is predicted using mixture theory, which models the dynamic advection of fluids through a fully saturated porous solid matrix. The resulting mixed formulation predicts all field variables including the solid displacement, pore fluid pressure and Darcy velocity of the pore fluid. The contact constraints arising from the requirement for continuity of the contact traction, as well as the fluid flow across the contact interface, are enforced using a penalty approach that is regularised with an augmented Lagrangian method. The contact formulation is based on a mortar segment‐to‐segment scheme that allows the interpolation functions of the contact elements to be of order N. The main thrust of this paper is therefore how to deal with contact interfaces in problems that involve both dynamics and consolidation and possibly large deformations of porous media. The numerical algorithm is first verified using several illustrative examples. This algorithm is then employed to solve a pipe‐seabed interaction problem, involving large deformations and dynamic effects, and the results of the analysis are also compared with those obtained using a node‐to‐segment contact algorithm. The results of this study indicate that the proposed method is able to solve the highly nonlinear problem of dynamic soil–structure interaction when coupled with pore water pressures and Darcy velocity. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
6.
This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek−1 is 𝒪(1). These variations are estimated to be less than 5% if Ek−1<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
7.
A program for the simulation of two‐dimensional (2‐D) fluid flow at the microstructural level of a saturated anisotropic granular medium is presented. The program provides a numerical solution to the complete set of Navier–Stokes equations without a priori assumptions on the viscous or convection components. This is especially suited for the simulation of the flow of fluids with different density and viscosity values and for a wide range of granular material porosity. The analytical solution for fluid flow in a simple microstructure of porous medium is used to verify the computer program. Subsequently, the flow field is computed within microscopic images of granular material that differ in porosity, particle size and particle shape. The computed flow fields are shown to follow certain paths depending on air void size and connectivity. The permeability tensor coefficients are derived from the flow fields, and their values are shown to compare well with laboratory experimental data on glass beads, Ottawa sand and silica sands. The directional distribution of permeability is expressed in a functional form and its anisotropy is quantified. Permeability anisotropy is found to be more pronounced in the silica sand medium that consists of elongated particles. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
8.
Zhendong Shan Daosheng Ling Zhinan Xie Liping Jing Yongqiang Li 《国际地质力学数值与分析法杂志》2019,43(13):2184-2199
One-dimensional transient wave propagation in a saturated single-layer porous medium with a fluid surface layer is studied in this paper. An analytical solution for a special case with a dynamic permeability coefficient kf → ∞ and a semianalytical solution for a general case with an arbitrary dynamic permeability coefficient are presented. The eigenfunction expansion and precise time step integration methods are employed. The solution is presented in series form, and thus, the long-term dynamic responses of saturated porous media with small permeability coefficients can be easily computed. We first transform the nonhomogeneous boundary conditions into homogeneous boundary conditions, and then we obtain the eigenvalues and orthogonal eigenfunctions of the fluid–solid system. Finally, the solutions in the time domain are developed. As the model is one dimensional, geometric attenuation is absent, and only the attenuation in the saturated porous medium is considered. We can apply this model to analyse the influences of different seabed types on the propagation of acoustic waves in the fluid layer, which is very important in ocean acoustics and ocean seismic. This solution can also be employed to validate the accuracies of various numerical methods. 相似文献
9.
There are many expressions proposed for the permeability of isotropic media based on flow channel and pore size distribution concepts, but there are no such expressions for anisotropic media. In this paper an expression for the permeability of an anisotropic medium is proposed, which has been verified in the laboratory. The mechanism behind fluid flow through soil was investigated using microscopic computer simulations to propose an expression for macroscopic permeability. The soil was assumed to be a spatially periodic porous medium, and the Navier-Stokes equation was solved using the FEM with appropriate boundary conditions for several different arrangements of the porous medium. The basic variables influencing flow through soil at the microscopic level were identified as specific surface area, void ratio, particle shape, material heterogeneity and the arrangement of particles in a porous medium. A sensitivity analysis was carried out to obtain an expression for the permeability in terms of the above variables. The corresponding macroscopic variables for the above microscopic variables are average specific surface area, average void ratio, anisotropy, tortuosity due to material heterogeneity, and the arrangement of particles respectively. An expression for the directional permeability is proposed in terms of these variables for the most common occurrence of particles in a porous medium. For the verification of the proposed equation, the permeability values of a fine-grained sand were measured at different void ratios and were compared with those predicted by the proposed equation. The results show that the predicted permeability values from the proposed equation are very close to the measured values. 相似文献
10.
Upscaling electroosmosis in porous media is a challenge due to the complexity and scale-dependent nonlinearities of this coupled phenomenon. “Pore-network modeling” for upscaling electroosmosis from pore scale to Darcy scale can be considered as a promising approach. However, this method requires analytical solutions for flow and transport at pore scale. This study concentrates on the development of analytical solutions of flow and transport in a single rectangular channel under combined effects of electrohydrodynamic forces. These relations will be used in future works for pore-network modeling. The analytical solutions are valid for all regimes of overlapping electrical double layers and have the potential to be extended to nonlinear Boltzmann distribution. The innovative aspects of this study are (a) contribution of overlapping of electrical double layers to the Stokes flow as well as Nernst–Planck transport has been carefully included in the analytical solutions. (b) All important transport mechanisms including advection, diffusion, and electromigration have been included in the analytical solutions. (c) Fully algebraic relations developed in this study can be easily employed to upscale electroosmosis to Darcy scale using pore-network modeling. 相似文献
11.
Visualizing complex pore structure and fluid flow in porous media using 3D printing technology and LBM simulation 下载免费PDF全文
下载免费PDF全文 Pore structures of porous media and properties of fluid flow are key factors for the study of non-Darcy groundwater flow. However, it is difficult to directly observe pore structures and flow properties, resulting in a “black box” problem of porous media. This problem has hindered the in-depth study of the groundwater flow mechanism at the pore scale. In recent years, 3D rapid prototyping technology has seen tremendous development. 3D printing provides digital models and printing models of porous media with clear internal structure. Thus, Lattice Boltzmann Method can be used to simulate the flow processes at the pore scale based on real pore structures. In this study, 3D printing cores and Lattice Boltzmann Method were coupled to conduct both laboratory and numerical experiments in spherical porous media with different sphere diameters and periodic arrays. The LBM simulation results show a good agreement with laboratory experimental results. With the advantages of LBM and 3D printing, this approach provides a visualization of the complex pore structure and fluid flow in pores, which is a promising method for studies of non-Darcy groundwater flow at the pore scale. 相似文献
12.
The coupled heat-fluid-stress problem of circular wellbore or spherical cavity subjected to a constant temperature change and a constant fluid flow rate is considered. Transient analytical solutions for temperature, pore pressure and stress are developed by coupling conductive heat transfer with Darcy fluid flow in a poroelastic medium. They are applicable to low permeability porous media suitable for liquid-waste disposal and also simulating reservoir for enhanced oil recovery, where conduction dominates the heat transfer process. A full range of solutions is presented showing separately the effects of temperature and fluid flow on pore pressure and stress development. It is shown that injection of warm fluid can be used to restrict fracture development around wellbores and cavities and generally to optimize a fluid injection operation. Both the limitations of the solutions and the convective flow effect are addressed. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
13.
深入探究孔隙尺度下的流体流动特性和溶质运移规律对石油开采、农田养分管理、地下水污染修复有着重要意义。以人工构建的多孔介质结构和同步辐射X射线显微CT扫描的土壤团聚体(分辨率3.7μm)为研究对象,在空间节点数多达64 000 000的情况下,基于格子Boltzmann模型和GPU并行技术计算得到多孔介质流体运动和溶质运移过程的关键参数,并据此探究多孔介质空间异质性对水力学特性的影响。通过对3组不同结构的多孔介质比较发现,结构复杂程度最高的土壤样品和不规则堆叠的圆球结构的渗透率在100 mD(即10^-13m^2)量级,远低于规则堆叠的圆球结构(>20 000 mD);土壤的迂曲度为1.40~1.60,明显高于规则堆叠的圆球结构。研究结果表明,渗透率大的样品具有较小的迂曲度,这与结构的空间异质性有较强的关系;土壤的渗透率和迂曲度呈现各向异性;在水力梯度一定的前提下,渗透率较大的样品,纵向弥散系数也较大;同时,结构的异质性也会影响溶质的穿透曲线。本研究提出的模拟方法可在土壤结构中进行高效的水流运动和溶质运移模拟,可用于土壤多孔介质在孔隙尺度下的水力学特性研究。 相似文献
14.
We consider a stationary flow of an incompressible non-Newtonian flow through a porous medium, induced by an injection velocity when inertial effects are negligible. At the pore scale, the governing equations are based on a nonlinear relation between the stress and the rate of deformation. In such a situation, the limit problem obtained when the pore size tends to zero, is called the homogenized problem that leads to the filtration law. This filtration law is given by a non-linear system coupling a local problem on a typical cell of the porous medium to a global problem at the scale of the whole porous medium. We propose, in this work, a numerical method to solve this homogenized problem and apply this method when the velocity dependent viscosity is given by the power law. Finally, we propose some numerical experiments to illustrate our approach. 相似文献
15.
Mathematical modelling of the ascent of free fluid through relatively strong rock, deep in the Earth's mantle, presents a challenge in geomechanics. Here the medium is considered as fluid-saturated, porous, elastic and bounded, and the fluid enters at a point source. An explicit finite difference method is developed for the numerical solution to the problem of the dilatation of a fluid-saturated porous elastic sphere due to a point fluid source of constant strength at the centre of the sphere. A cubic spline interpolant is used to evaluate a definite integral which occurs in the boundary condition for the pore fluid pressure at the surface of the sphere. The numerical solutions for the dilatation and pore fluid pressure are compared with analytical solutions and the absolute and relative errors of the numerical solutions are calculated. When the fluid source is switched on, the pore fluid pressure starts to decrease, reaches a minimum value and then steadily increases. The initial time rate of decrease of the pore fluid pressure is independent of the radial distance from the source. It decreases as the radius of the sphere increases and vanishes for a point fluid source in an infinite porous elastic medium. 相似文献
16.
深入探究孔隙尺度下的流体流动特性和溶质运移规律对石油开采、农田养分管理、地下水污染修复有着重要意义。以人工构建的多孔介质结构和同步辐射X射线显微CT扫描的土壤团聚体（分辨率3.7 μm）为研究对象，在空间节点数多达64 000 000的情况下，基于格子Boltzmann模型和GPU并行技术计算得到多孔介质流体运动和溶质运移过程的关键参数，并据此探究多孔介质空间异质性对水力学特性的影响。通过对3组不同结构的多孔介质比较发现，结构复杂程度最高的土壤样品和不规则堆叠的圆球结构的渗透率在100 mD（即10-13m2）量级，远低于规则堆叠的圆球结构（>20 000 mD）；土壤的迂曲度为1.40~1.60，明显高于规则堆叠的圆球结构。研究结果表明，渗透率大的样品具有较小的迂曲度，这与结构的空间异质性有较强的关系；土壤的渗透率和迂曲度呈现各向异性；在水力梯度一定的前提下，渗透率较大的样品，纵向弥散系数也较大；同时，结构的异质性也会影响溶质的穿透曲线。本研究提出的模拟方法可在土壤结构中进行高效的水流运动和溶质运移模拟，可用于土壤多孔介质在孔隙尺度下的水力学特性研究。 相似文献
17.
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. 相似文献
18.
Guillermo A. Narsilio Olivier Buzzi Stephen Fityus Tae Sup Yun David W. Smith 《Computers and Geotechnics》2009,36(7):1200-1206
The accurate estimation of hydraulic conductivity is important for many geotechnical engineering applications, as the presence of fluids affects all aspects of soil behaviour, including its strength. Darcy’s law is the key experimental (or phenomenological) equation employed to model ground water flow. Yet, this phenomenological equation can be linked to a more fundamental microscale model of flow through the pore spaces of the porous material. This paper provides an experimental verification of the relationships between Darcy’s law (macroscale) and the Navier–Stokes equations (microscale) for actual complex pore geometries of a granular material. The pore geometries are experimentally obtained through state-of-the-art X-ray computer assisted micro-tomography. From the numerical modelling of the microscale flow based on actual pore geometries, it is possible to quantify and visualize the development of pore-scale fluid preferential flow-paths through the porous material, and to assess the importance of pore connectivity in soil transport properties. 相似文献
19.
20.
This paper presents a novel analytical solution to the transient, z‐dependent, and asymmetric problem of an infinite wellbore drilled into a fluid‐saturated porous medium. The formulations are based on Biot's linear theory of poroelasticity, in which the dependency of poroelastic field variables to spatial coordinates as well as time domain is considered in the most general form. This gives flexibility to the solution in cases that cannot be analyzed using the conventional plane strain or symmetric models. One such case is when calculating the stress variations around an inclined wellbore where the far‐field stresses are acting over a finite vertical section. The results of our solution to this case with a three‐dimensional state of far‐field stress are used to analyze the stability of inclined wellbores passing through abnormally stressed formations. The presented solution is capable of finding expressions for fundamental solutions with stress or flow boundary conditions at the wellbore. These solutions are here adopted to analyze the pressure disturbances generated by multiprobe formation tester, a standard wireline device that is designed for downhole fluid sampling as well as estimating the directional permeabilities of subsurface earth formations. A comparison with the conventional solution for the relevant pressure diffusion equation indicates that the poroelastic effect is relatively significant in relation to the transient response of the pore pressure. Further, it is shown that the finite dimensions of sink probe would, to a great extent, contribute to the formation's pore pressure variations at its immediate proximity. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献