In this paper,we propose a multiphysics finite element method for a nonlinear poroelasticity model with nonlinear stress-strain relation.Firstly,we reformulate the original problem into a new coupled fluid system-a ge...In this paper,we propose a multiphysics finite element method for a nonlinear poroelasticity model with nonlinear stress-strain relation.Firstly,we reformulate the original problem into a new coupled fluid system-a generalized nonlinear Stokes problem of displacement vector field related to pseudo pressure and a diffusion problem of other pseudo pressure fields.Secondly,a fully discrete multiphysics finite element method is performed to solve the reformulated system numerically.Thirdly,existence and uniqueness of the weak solution of the reformulated model and stability analysis and optimal convergence order for the multiphysics finite element method are proven theoretically.Lastly,numerical tests are given to verify the theoretical results.展开更多
In this paper,we design a new error estimator and give a posteriori error analysis for a poroelasticity model.To better overcome“locking phenomenon”on pressure and displacement,we proposed a new error estimators bas...In this paper,we design a new error estimator and give a posteriori error analysis for a poroelasticity model.To better overcome“locking phenomenon”on pressure and displacement,we proposed a new error estimators based on multiphysics discontinuous Galerkin method for the poroelasticity model.And we prove the upper and lower bound of the proposed error estimators,which are numerically demonstrated to be computationally very efficient.Finally,we present numerical examples to verify and validate the efficiency of the proposed error estimators,which show that the adaptive scheme can overcome“locking phenomenon”and greatly reduce the computation cost.展开更多
An elastomeric gel is a mixture of a polymer network and a solvent. In response to changes in mechanical forces and in the chemical potential of the solvent in the environment, the gel evolves by two concurrent molecu...An elastomeric gel is a mixture of a polymer network and a solvent. In response to changes in mechanical forces and in the chemical potential of the solvent in the environment, the gel evolves by two concurrent molecular processes： the conformational change of the network, and the migration of the solvent. The two processes result in viscoelasticity and poroelasticity, and are characterized by two material-specific properties： the time of viscoelastic relaxation and the effective diffusivity of the solvent through the network. The two properties define a material- specific length. The material-specific time and length enable us to discuss macroscopic observations made over different lengths and times, and identify limiting conditions in which viscoelastic and poroelastic relaxations have either completed or yet started. We formulate a model of homogeneous deformation, and use several examples to illustrate viscoelasticity-limited solvent migration, where the migration of the solvent is pronounced, but the size of the gel is so small that the rate of change is limited by viscoelasticity. We further describe a theory that evolves a gel through inhomogeneous states. Both infinitesimal and finite deformation are considered.展开更多
Seismic amplitude variation with offset(AVO) inversion is an important approach for quantitative prediction of rock elasticity,lithology and fluid properties.With Biot-Gassmann's poroelasticity,an improved statist...Seismic amplitude variation with offset(AVO) inversion is an important approach for quantitative prediction of rock elasticity,lithology and fluid properties.With Biot-Gassmann's poroelasticity,an improved statistical AVO inversion approach is proposed.To distinguish the influence of rock porosity and pore fluid modulus on AVO reflection coefficients,the AVO equation of reflection coefficients parameterized by porosity,rock-matrix moduli,density and fluid modulus is initially derived from Gassmann equation and critical porosity model.From the analysis of the influences of model parameters on the proposed AVO equation,rock porosity has the greatest influences,followed by rock-matrix moduli and density,and fluid modulus has the least influences among these model parameters.Furthermore,a statistical AVO stepwise inversion method is implemented to the simultaneous estimation of rock porosity,rock-matrix modulus,density and fluid modulus.Besides,the Laplace probability model and differential evolution,Markov chain Monte Carlo algorithm is utilized for the stochastic simulation within Bayesian framework.Models and field data examples demonstrate that the simultaneous optimizations of multiple Markov chains can achieve the efficient simulation of the posterior probability density distribution of model parameters,which is helpful for the uncertainty analysis of the inversion and sets a theoretical fundament for reservoir characterization and fluid discrimination.展开更多
Abstract Hydrogel can swell to many times of its dry volume, resulting in large deformation which is vital for its function. The swelling process is regulated by many physical and chemical mechanisms, and can, to some...Abstract Hydrogel can swell to many times of its dry volume, resulting in large deformation which is vital for its function. The swelling process is regulated by many physical and chemical mechanisms, and can, to some extent, be fairly described by the poroelasticity theory. Implementation of the poroelastieity theory in the framework of finite element method would aid the design and optimization of hydrogel-based soft devices. Choosing chemical potential and displacement as two field variables, we present the implementation of poroelastieity tailored for hydrogel swelling dynamics, detail the normalization of physical parameters and the treatment of boundary conditions. Several examples are presented to demonstrate the feasibility and correctness of the proposed strategy.展开更多
In this paper,we develop a fully discrete virtual element scheme based on the local pressure projection stabilization for a three-field poroelasticity problem with a storage coefficient c00.We not only provide the wel...In this paper,we develop a fully discrete virtual element scheme based on the local pressure projection stabilization for a three-field poroelasticity problem with a storage coefficient c00.We not only provide the well-posedness of the proposed scheme by proving a weaker form of the discrete inf-sup condition,but also show optimal error estimates for all unknowns,whose generic constants are independent of the Lam´e coefficient.Moreover,our proposed scheme avoids pressure oscillation and applies to general polygonal elements,including hanging-node elements.Finally,we numerically validate the good performance of our virtual element scheme.展开更多
In this paper,we will investigate a multigrid algorithm for poroelasticity problem by a new finite element method with homogeneous boundary conditions in two dimensional space.We choose N´ed´elec edge elemen...In this paper,we will investigate a multigrid algorithm for poroelasticity problem by a new finite element method with homogeneous boundary conditions in two dimensional space.We choose N´ed´elec edge element for the displacement variable and piecewise continuous polynomials for the pressure variable in the model problem.In constructing multigrid algorithm,a distributive Gauss-Seidel iteration method is applied.Numerical experiments shows that the finite element method achieves optimal convergence order and the multigrid algorithm is almost uniformly convergent to mesh size h and parameter dt on regular meshes.展开更多
We consider a strongly heterogeneous medium saturated by an incompressible viscous fluid as it appears in geomechanical modeling.This poroelasticity problem suffers from rapidly oscillating material parameters,which c...We consider a strongly heterogeneous medium saturated by an incompressible viscous fluid as it appears in geomechanical modeling.This poroelasticity problem suffers from rapidly oscillating material parameters,which calls for a thorough numerical treatment.In this paper,we propose a method based on the local orthogonal decomposition technique and motivated by a similar approach used for linear thermoelasticity.Therein,local corrector problems are constructed in line with the static equations,whereas we propose to consider the full system.This allows to benefit from the given saddle point structure and results in two decoupled corrector problems for the displacement and the pressure.We prove the optimal first-order convergence of this method and verify the result by numerical experiments.展开更多
In this paper,a multirate time iterative scheme with multiphysics finite element method is proposed and analyzed for the nonlinear poroelasticity model.The original problem is reformulated into a generalized nonlinear...In this paper,a multirate time iterative scheme with multiphysics finite element method is proposed and analyzed for the nonlinear poroelasticity model.The original problem is reformulated into a generalized nonlinear Stokes problem coupled with a diffusion problem of a pseudo pressure field by a new multiphysics approach.A multiphysics finite element method is adopted for the spatial discretization,and the generalized nonlinear Stokes problem is solved in a coarse time step and the diffusion problem is solved in a finer time step.The proposed algorithm is a decoupled algorithm,which is easily implemented in computation and reduces greatly computation cost.The stability analysis and the convergence analysis for the multirate iterative scheme with multiphysics finite element method are given.Some numerical tests are shown to demonstrate and validate the analysis results.展开更多
Understanding the effects of point liquid loading on transversely isotropic poroelastic media is crucial for advancing geomechanics and biomechanics, where precise modeling of fluid-structure interactions is essential...Understanding the effects of point liquid loading on transversely isotropic poroelastic media is crucial for advancing geomechanics and biomechanics, where precise modeling of fluid-structure interactions is essential. This paper presents a comprehensive analysis of infinite transversely isotropic poroelasticity under a fluid source, based on Biot's theory, aiming to uncover new and previously unexplored insights in the literature. We begin our study by deriving a general solution for fluid-saturated, transversely isotropic poroelastic materials in terms of harmonic functions that satisfy sixth-order homogeneous partial differential equations, using potential theory and Almansi's theorem. Based on these general solutions and potential functions, we construct a Green's function for a point fluid source, introducing three new harmonic functions with undetermined constants. These constants are determined by enforcing continuity and equilibrium conditions. Substituting these into the general solution yields fundamental solutions for poroelasticity that provide crucial support for a wide range of project problems. Numerical results and comparisons with existing literature are provided to illustrate physical mechanisms through contour plots. Our observations reveal that all components tend to zero in the far field and become singular at the concentrated source. Additionally, the contours exhibit rapid changes near the point fluid source but display gradual variations at a distance from it. These findings highlight the intricate behavior of the system under point liquid loading, offering valuable insights for further research and practical applications.展开更多
In order to quantify the poroelastic mechanical signals conduction and evaluate the biomechanical effectiveness of functional units(osteocyte processes,canaliculi and lacuna)in lacunar-canalicular system(LCS),a multis...In order to quantify the poroelastic mechanical signals conduction and evaluate the biomechanical effectiveness of functional units(osteocyte processes,canaliculi and lacuna)in lacunar-canalicular system(LCS),a multiscale poroelastic finite element model was established by using the Comsol Multiphysics software.The poroelastic mechanical signals(pore pressure,fluid velocity,von-Mises stress,strain)were analyzed inside the osteon-osteocyte system.The effects of osteocyte(OCY)’s shape(ellipse and circle),long axis directions(horizontal and vertical)and mechanical properties(Elastic modulus and permeability)on its poroelastic responses were examined.It is found that the OCY processes is the best mechanosensor compared with the OCY body,lacunae and canaliculi.The mechanotransduction ability of the elliptic shaped OCY is stronger than that of circular shaped.The pore pressure and flow velocity around OCYs increase as the elastic modulus and permeability of OCY increase.The established model can be used for studying the mechanism of bone mechanotransduction at the multiscale level.展开更多
Connecting earthquake nucleation in basement rock to fluid injection in basal,sedimentary reservoirs,depends heavily on choices related to the poroelastic properties of the fluid-rock system,thermo-chemical effects no...Connecting earthquake nucleation in basement rock to fluid injection in basal,sedimentary reservoirs,depends heavily on choices related to the poroelastic properties of the fluid-rock system,thermo-chemical effects notwithstanding.Direct constraints on these parameters outside of laboratory settings are rare,and it is commonly assumed that the rock layers are isotropic.With the Arbuckle wastewater disposal reservoir in Osage County,Oklahoma,high-frequency formation pressure changes and collocated broadband ground velocities measured during the passing of large teleseismic waves show a poroelastic response of the reservoir that is both azimuthally variable and anisotropic;this includes evidence of static shifts in pressure that presumably relate to changes in local permeability.The azimuthal dependence in both the static response and shear coupling appears related to tectonic stress and strain indicators such as the orientations of the maximum horizontal stress and faults and fractures.Using dynamic strains from a nearby borehole strainmeter,we show that the ratio of shear to volumetric strain coupling is~0.41 which implies a mean Skempton's coefficient of A=0.24 over the plausible range of the undrained Poisson's ratio.Since these observations are made at relatively low confining pressure and differential stress,we suggest that the hydraulically conductive fracture network is a primary control on the coupling between pore pressure diffusion and elastic stresses in response to natural or anthropogenic sources.展开更多
From 2009 to 2017,parts of Central America experienced marked increase in the number of small to moderate-sized earthquakes.For example,three significant earthquakes(~Mw 5)occurred near Prague,Oklahoma,in the U.S.in 2...From 2009 to 2017,parts of Central America experienced marked increase in the number of small to moderate-sized earthquakes.For example,three significant earthquakes(~Mw 5)occurred near Prague,Oklahoma,in the U.S.in 2011.On 6 Nov 2011,an Mw 5.7 earthquake occurred in Prague,central Oklahoma with a sequence of aftershocks.The seismic activity has been attributed to slip on the Wilzetta fault system.This study provides a 3 D fully coupled poroelastic analysis(using FLAC3 D)of the Wilzetta fault system and its response to saltwater injection in the underpressured subsurface layers,especially the Arbuckle group and the basement,to evaluate the conditions that might have led to the increased seismicity.Given the data-limited nature of the problem,we have considered multiple plausible scenarios,and use the available data to evaluate the hydromechanical response of the faults of interest in the study area.Numerical simulations show that the injection of large volumes of fluid into the Arbuckle group tends to bring the part of the Wilzetta faults in Arbuckle group and basement into near-critical conditions.展开更多
A linear viscoporoelastic model is developed to describe the problem of reflection and transmission of an obliquely incident plane P-wave at the interface between an elastic solid and an unsaturated poroelastic medium...A linear viscoporoelastic model is developed to describe the problem of reflection and transmission of an obliquely incident plane P-wave at the interface between an elastic solid and an unsaturated poroelastic medium, in which the solid matrix is filled with two weakly coupled fluids （liquid and gas）. The expressions for the amplitude reflection coefficients and the amplitude transmission coefficients are derived by using the potential method. The present derivation is subsequently applied to study the energy conversions among the incident, reflected, and transmitted wave modes. It is found that the reflection and transmission coefficients in the forms of amplitude ratios and energy ratios are functions of the incident angle, the liquid saturation, the frequency of the incident wave, and the elastic constants of the upper and lower media. Numerical results are presented graphically. The effects of the incident angle, the frequency, and the liquid saturation on the amplitude and the energy reflection and transmission coefficients are discussed. It is verified that in the transmission process, there is no energy dissipation at the interface.展开更多
A mathematical formulation is presented for the dynamic stress intensity factor (mode I) of a finite permeable crack subjected to a time-harmonic propagating longitudinal wave in an infinite poroelastic solid. In part...A mathematical formulation is presented for the dynamic stress intensity factor (mode I) of a finite permeable crack subjected to a time-harmonic propagating longitudinal wave in an infinite poroelastic solid. In particular, the effect of the wave-induced fluid flow due to the presence of a liquid-saturated crack on the dynamic stress intensity factor is analyzed. Fourier sine and cosine integral transforms in conjunction with Helmholtz potential theory are used to formulate the mixed boundary-value problem as dual integral equations in the frequency domain. The dual integral equations are reduced to a Fredholm integral equation of the second kind. It is found that the stress intensity factor monotonically decreases with increasing frequency, decreasing the fastest when the crack width and the slow wave wavelength are of the same order. The characteristic frequency at which the stress intensity factor decays the fastest shifts to higher frequency values when the crack width decreases.展开更多
Geological sequestration (GS) of carbon dioxide (CO2) is considered as one of the most promising technologies to reduce the amount of anthropogenic CO2 emission in the atmosphere. To ensure success of CO2 GS, monitori...Geological sequestration (GS) of carbon dioxide (CO2) is considered as one of the most promising technologies to reduce the amount of anthropogenic CO2 emission in the atmosphere. To ensure success of CO2 GS, monitoring is essential on ascertaining movement, volumes and locations of injected CO2 in the sequestration reservoir. One technique is to use time-lapsed seismic survey mapping to provide spatial distribution of seismic wave velocity as an indicator of CO2 migration and volumes in a storage reservoir with time. To examine the use of time-lapsed seismic survey mapping as a monitoring tool for CO2 sequestration, this paper presents mathematical and experimental studies of the effects of supercritical CO2 injection on the seismic velocity of sandstone initially saturated with saline water. The mathematical model is based on poroelasticity theory, particularly the application of the Biot-Gassmann substitution theory in the modeling of the acoustic velocity of porous rocks containing two-phase immiscible pore fluids. The experimental study uses a high pressure and high temperature triaxial cell to clarify the seismic response of a sample of Berea sandstone to supercritical CO2 injection under deep saline aquifer conditions. Measured ultrasonic wave velocity changes during CO2 injection in the sandstone sample show the effects of pore fluid distribution in the seismic velocity of porous rocks. CO2 injection was shown to decrease the P-wave velocity with increasing CO2 saturation whereas the S-wave velocity was almost constant. The results confirm that the Biot-Gassmann theory can be used to model the changes in the acoustic P-wave velocity of sandstone containing different mixtures of supercritical CO2 and saline water provided the distribution of the two fluids in the sandstone pore space is accounted for in the calculation of the pore fluid bulk modulus. The empirical relation of Brie et al. for the bulk modulus of mixtures of two-phase immiscible fluids, in combination with the Biot-Gassmann theory, was found to satisfactorily represent the pore-fluid dependent acoustic P-wave velocity of sandstone.展开更多
Numerical modeling of seepage-induced consolidation process usually encounters significant uncertainty in the properties of geotechnical materials.Assessing the effect of uncertain parameters on the performance variab...Numerical modeling of seepage-induced consolidation process usually encounters significant uncertainty in the properties of geotechnical materials.Assessing the effect of uncertain parameters on the performance variability of the seepage consolidation model is of critical importance to the simulation and tests of this process.To this end,the uncertainty and sensitivity analyses are performed on a seepage consolidation model in a fractured porous medium using the Bayesian sparse polynomial chaos expansion(SPCE)method.Five uncertain parameters including Young’s modulus,Poisson’s ratio,and the permeability of the porous matrix,the permeability within the fracture,and Biot’s constant are studied.Bayesian SPCE models for displacement,flow velocity magnitude,and fluid pressure at several reference points are constructed to represent the input-output relationship of the numerical model.Based on these SPCE models,the total and first-order Sobol’indices are computed to quantify the contribution of each uncertain input parameter to the uncertainty of model responses.The results show that at different locations of the porous domain,the uncertain parameters show different effects on the output quantities.At the beginning of the seepage consolidation process,the hydraulic parameters make major contributions to the uncertainty of the model responses.As the process progresses,the effect of hydraulic parameters decreases and is gradually surpassed by the mechanical parameters.This work demonstrates the feasibility to apply Bayesian SPCE approach to the uncertainty and sensitivity analyses of seepage-induced consolidation problems and provides guidelines to the numerical modelling and experimental testing of such problems.展开更多
The effect of the fracture distribution on CO 2 injection into coal seams was studied with a heterogeneous model having dual porosity to represent both the primary medium (the coal matrix) and the secondary medium (th...The effect of the fracture distribution on CO 2 injection into coal seams was studied with a heterogeneous model having dual porosity to represent both the primary medium (the coal matrix) and the secondary medium (the fractures) under variable stress conditions. A numerical generation method and a digital image processing method were used to model the heterogeneous fracture distribution in the coal. The model solutions demonstrate that: (1) the fractures are the main channel for gas flow and their distribution has an important impact on the gas injection rate; (2) the fractures only affect the injection rate of CO 2 into the coal but not the final storage amount; (3) when gas is injected into coal the fractures will first expand and then close due to the changing effective stresses and the adsorption induced swelling of coal grains. This fully coupled dual-porosity model with a heterogeneous fracture distribution provides a way to predict the CO 2 injection into a coal seam.展开更多
An extended and reasonable stress boundary condition at an osteon exte- rior wall is presented to solve the model proposed by Remond and Naili. The obtained pressure and fluid velocity solutions are used to investigat...An extended and reasonable stress boundary condition at an osteon exte- rior wall is presented to solve the model proposed by Remond and Naili. The obtained pressure and fluid velocity solutions are used to investigate the osteonal poroelastic behaviors. The following results are obtained. （i） Both the fluid pressure and the velocity amplitudes are proportional to the strain amplitude and the loading frequency. （ii） In the physiological loading state, the key role governing the poroelastic behaviors of the osteon is the strain rate. （iii） At the osteon scale, the pressure is strongly affected by the permeability variations, whereas the fluid velocity is not.展开更多
Pre-stack seismic inversion is an important method for fluid identification and reservoir characterization in exploration geophysics. In this study, an effective fluid factor is initially established based on Biot por...Pre-stack seismic inversion is an important method for fluid identification and reservoir characterization in exploration geophysics. In this study, an effective fluid factor is initially established based on Biot poroelastic theory, and a pre-stack seismic inversion method based on Bayesian framework is used to implement the fluid identification. Compared with conventional elastic parameters, fluid factors are more sensitive to oil and gas. However, the coupling effect between rock porosity and fluid content is not considered in conventional fluid factors, which may lead to fuzzy fluid identification results. In addition,existing fluid factors do not adequately consider the physical mechanisms of fluid content, such as squirt flow between cracks and pores. Therefore, we propose a squirt fluid factor(SFF) that minimizes the fluid and pore mixing effects and takes into account the squirt flow. On this basis, a novel P-wave reflection coefficient equation is derived, and the squirt fluid factor is estimated by amplitude variation with offset(AVO) inversion method. The new reflection coefficient equation has sufficient accuracy and can be utilized to estimate the parameters. The effectiveness and superiority of the proposed method in fluid identification are verified by the synthetic and field examples.展开更多
基金Supported by the National Natural Science Foundation of China(Grant Nos.12371393,11971150 and 11801143)Natural Science Foundation of Henan Province(Grant No.242300421047).
文摘In this paper,we propose a multiphysics finite element method for a nonlinear poroelasticity model with nonlinear stress-strain relation.Firstly,we reformulate the original problem into a new coupled fluid system-a generalized nonlinear Stokes problem of displacement vector field related to pseudo pressure and a diffusion problem of other pseudo pressure fields.Secondly,a fully discrete multiphysics finite element method is performed to solve the reformulated system numerically.Thirdly,existence and uniqueness of the weak solution of the reformulated model and stability analysis and optimal convergence order for the multiphysics finite element method are proven theoretically.Lastly,numerical tests are given to verify the theoretical results.
基金supported by the National Natural Science Foundation of China(Grant Nos.12371393 and 11971150)Natural Science Foundation of Henan(Grant No.242300421047).
文摘In this paper,we design a new error estimator and give a posteriori error analysis for a poroelasticity model.To better overcome“locking phenomenon”on pressure and displacement,we proposed a new error estimators based on multiphysics discontinuous Galerkin method for the poroelasticity model.And we prove the upper and lower bound of the proposed error estimators,which are numerically demonstrated to be computationally very efficient.Finally,we present numerical examples to verify and validate the efficiency of the proposed error estimators,which show that the adaptive scheme can overcome“locking phenomenon”and greatly reduce the computation cost.
基金supported by the National Science Foundation (NSF) (CMMI-0800161)Multidisciplinary University Research Initiative (MURI) (W911NF-09-1-0476)the Materials Research Science and Engineering Center at Harvard University (DMR-0820484)
文摘An elastomeric gel is a mixture of a polymer network and a solvent. In response to changes in mechanical forces and in the chemical potential of the solvent in the environment, the gel evolves by two concurrent molecular processes： the conformational change of the network, and the migration of the solvent. The two processes result in viscoelasticity and poroelasticity, and are characterized by two material-specific properties： the time of viscoelastic relaxation and the effective diffusivity of the solvent through the network. The two properties define a material- specific length. The material-specific time and length enable us to discuss macroscopic observations made over different lengths and times, and identify limiting conditions in which viscoelastic and poroelastic relaxations have either completed or yet started. We formulate a model of homogeneous deformation, and use several examples to illustrate viscoelasticity-limited solvent migration, where the migration of the solvent is pronounced, but the size of the gel is so small that the rate of change is limited by viscoelasticity. We further describe a theory that evolves a gel through inhomogeneous states. Both infinitesimal and finite deformation are considered.
基金the sponsorship of National Grand Project for Science and Technology(2016ZX05024004,2017ZX05009001,2017ZX05032003)the Fundamental Research Funds for the Central Universities(20CX06036A)+1 种基金the Postdoctoral Applied Research Project of Qingdao(QDYY20190040)the Science Foundation from SINOPEC Key Laboratory of Geophysics(wtyjy-wx2019-01-04)。
文摘Seismic amplitude variation with offset(AVO) inversion is an important approach for quantitative prediction of rock elasticity,lithology and fluid properties.With Biot-Gassmann's poroelasticity,an improved statistical AVO inversion approach is proposed.To distinguish the influence of rock porosity and pore fluid modulus on AVO reflection coefficients,the AVO equation of reflection coefficients parameterized by porosity,rock-matrix moduli,density and fluid modulus is initially derived from Gassmann equation and critical porosity model.From the analysis of the influences of model parameters on the proposed AVO equation,rock porosity has the greatest influences,followed by rock-matrix moduli and density,and fluid modulus has the least influences among these model parameters.Furthermore,a statistical AVO stepwise inversion method is implemented to the simultaneous estimation of rock porosity,rock-matrix modulus,density and fluid modulus.Besides,the Laplace probability model and differential evolution,Markov chain Monte Carlo algorithm is utilized for the stochastic simulation within Bayesian framework.Models and field data examples demonstrate that the simultaneous optimizations of multiple Markov chains can achieve the efficient simulation of the posterior probability density distribution of model parameters,which is helpful for the uncertainty analysis of the inversion and sets a theoretical fundament for reservoir characterization and fluid discrimination.
基金supported by the National Natural Science Foundation of China(11072185,11372239,and 11021202)
文摘Abstract Hydrogel can swell to many times of its dry volume, resulting in large deformation which is vital for its function. The swelling process is regulated by many physical and chemical mechanisms, and can, to some extent, be fairly described by the poroelasticity theory. Implementation of the poroelastieity theory in the framework of finite element method would aid the design and optimization of hydrogel-based soft devices. Choosing chemical potential and displacement as two field variables, we present the implementation of poroelastieity tailored for hydrogel swelling dynamics, detail the normalization of physical parameters and the treatment of boundary conditions. Several examples are presented to demonstrate the feasibility and correctness of the proposed strategy.
基金supported in part by the National Natural Science Foundation of China(Grant No.12371405).
文摘In this paper,we develop a fully discrete virtual element scheme based on the local pressure projection stabilization for a three-field poroelasticity problem with a storage coefficient c00.We not only provide the well-posedness of the proposed scheme by proving a weaker form of the discrete inf-sup condition,but also show optimal error estimates for all unknowns,whose generic constants are independent of the Lam´e coefficient.Moreover,our proposed scheme avoids pressure oscillation and applies to general polygonal elements,including hanging-node elements.Finally,we numerically validate the good performance of our virtual element scheme.
基金The first author is supported by the National Natural Science Foundation of China under Grant No.11501473,No.11426189the Fundamental Research Funds for the Central Universities of China(No.2682016CX108)The second author is supported by the National Natural Science Foundation of China under Grant No.11671157.
文摘In this paper,we will investigate a multigrid algorithm for poroelasticity problem by a new finite element method with homogeneous boundary conditions in two dimensional space.We choose N´ed´elec edge element for the displacement variable and piecewise continuous polynomials for the pressure variable in the model problem.In constructing multigrid algorithm,a distributive Gauss-Seidel iteration method is applied.Numerical experiments shows that the finite element method achieves optimal convergence order and the multigrid algorithm is almost uniformly convergent to mesh size h and parameter dt on regular meshes.
文摘We consider a strongly heterogeneous medium saturated by an incompressible viscous fluid as it appears in geomechanical modeling.This poroelasticity problem suffers from rapidly oscillating material parameters,which calls for a thorough numerical treatment.In this paper,we propose a method based on the local orthogonal decomposition technique and motivated by a similar approach used for linear thermoelasticity.Therein,local corrector problems are constructed in line with the static equations,whereas we propose to consider the full system.This allows to benefit from the given saddle point structure and results in two decoupled corrector problems for the displacement and the pressure.We prove the optimal first-order convergence of this method and verify the result by numerical experiments.
基金supported by the National Natural Science Foundation of China(Grant No.11971150)partially by the National Natural Science Foundation of China(Grant No.11801143).
文摘In this paper,a multirate time iterative scheme with multiphysics finite element method is proposed and analyzed for the nonlinear poroelasticity model.The original problem is reformulated into a generalized nonlinear Stokes problem coupled with a diffusion problem of a pseudo pressure field by a new multiphysics approach.A multiphysics finite element method is adopted for the spatial discretization,and the generalized nonlinear Stokes problem is solved in a coarse time step and the diffusion problem is solved in a finer time step.The proposed algorithm is a decoupled algorithm,which is easily implemented in computation and reduces greatly computation cost.The stability analysis and the convergence analysis for the multirate iterative scheme with multiphysics finite element method are given.Some numerical tests are shown to demonstrate and validate the analysis results.
基金supported by the National Natural Science Foundation of China (Grant Nos. 12272269, 11972257,11832014 and 11472193)the Shanghai Pilot Program for Basic Researchthe Shanghai Gaofeng Project for University Academic Program Development。
文摘Understanding the effects of point liquid loading on transversely isotropic poroelastic media is crucial for advancing geomechanics and biomechanics, where precise modeling of fluid-structure interactions is essential. This paper presents a comprehensive analysis of infinite transversely isotropic poroelasticity under a fluid source, based on Biot's theory, aiming to uncover new and previously unexplored insights in the literature. We begin our study by deriving a general solution for fluid-saturated, transversely isotropic poroelastic materials in terms of harmonic functions that satisfy sixth-order homogeneous partial differential equations, using potential theory and Almansi's theorem. Based on these general solutions and potential functions, we construct a Green's function for a point fluid source, introducing three new harmonic functions with undetermined constants. These constants are determined by enforcing continuity and equilibrium conditions. Substituting these into the general solution yields fundamental solutions for poroelasticity that provide crucial support for a wide range of project problems. Numerical results and comparisons with existing literature are provided to illustrate physical mechanisms through contour plots. Our observations reveal that all components tend to zero in the far field and become singular at the concentrated source. Additionally, the contours exhibit rapid changes near the point fluid source but display gradual variations at a distance from it. These findings highlight the intricate behavior of the system under point liquid loading, offering valuable insights for further research and practical applications.
基金supported by the National Natural Science Foundation of China(Grants 11972242,11702183,11632013,and 11572213)the Scientific and Technological Innovation Projects of Colleges and Universities in Shanxi Province(Grant 2017135)Philosophy and Social Sciences Research of Higher Learning Institutions of Shanxi(Grant 2017313).
文摘In order to quantify the poroelastic mechanical signals conduction and evaluate the biomechanical effectiveness of functional units(osteocyte processes,canaliculi and lacuna)in lacunar-canalicular system(LCS),a multiscale poroelastic finite element model was established by using the Comsol Multiphysics software.The poroelastic mechanical signals(pore pressure,fluid velocity,von-Mises stress,strain)were analyzed inside the osteon-osteocyte system.The effects of osteocyte(OCY)’s shape(ellipse and circle),long axis directions(horizontal and vertical)and mechanical properties(Elastic modulus and permeability)on its poroelastic responses were examined.It is found that the OCY processes is the best mechanosensor compared with the OCY body,lacunae and canaliculi.The mechanotransduction ability of the elliptic shaped OCY is stronger than that of circular shaped.The pore pressure and flow velocity around OCYs increase as the elastic modulus and permeability of OCY increase.The established model can be used for studying the mechanism of bone mechanotransduction at the multiscale level.
文摘Connecting earthquake nucleation in basement rock to fluid injection in basal,sedimentary reservoirs,depends heavily on choices related to the poroelastic properties of the fluid-rock system,thermo-chemical effects notwithstanding.Direct constraints on these parameters outside of laboratory settings are rare,and it is commonly assumed that the rock layers are isotropic.With the Arbuckle wastewater disposal reservoir in Osage County,Oklahoma,high-frequency formation pressure changes and collocated broadband ground velocities measured during the passing of large teleseismic waves show a poroelastic response of the reservoir that is both azimuthally variable and anisotropic;this includes evidence of static shifts in pressure that presumably relate to changes in local permeability.The azimuthal dependence in both the static response and shear coupling appears related to tectonic stress and strain indicators such as the orientations of the maximum horizontal stress and faults and fractures.Using dynamic strains from a nearby borehole strainmeter,we show that the ratio of shear to volumetric strain coupling is~0.41 which implies a mean Skempton's coefficient of A=0.24 over the plausible range of the undrained Poisson's ratio.Since these observations are made at relatively low confining pressure and differential stress,we suggest that the hydraulically conductive fracture network is a primary control on the coupling between pore pressure diffusion and elastic stresses in response to natural or anthropogenic sources.
文摘From 2009 to 2017,parts of Central America experienced marked increase in the number of small to moderate-sized earthquakes.For example,three significant earthquakes(~Mw 5)occurred near Prague,Oklahoma,in the U.S.in 2011.On 6 Nov 2011,an Mw 5.7 earthquake occurred in Prague,central Oklahoma with a sequence of aftershocks.The seismic activity has been attributed to slip on the Wilzetta fault system.This study provides a 3 D fully coupled poroelastic analysis(using FLAC3 D)of the Wilzetta fault system and its response to saltwater injection in the underpressured subsurface layers,especially the Arbuckle group and the basement,to evaluate the conditions that might have led to the increased seismicity.Given the data-limited nature of the problem,we have considered multiple plausible scenarios,and use the available data to evaluate the hydromechanical response of the faults of interest in the study area.Numerical simulations show that the injection of large volumes of fluid into the Arbuckle group tends to bring the part of the Wilzetta faults in Arbuckle group and basement into near-critical conditions.
文摘A linear viscoporoelastic model is developed to describe the problem of reflection and transmission of an obliquely incident plane P-wave at the interface between an elastic solid and an unsaturated poroelastic medium, in which the solid matrix is filled with two weakly coupled fluids （liquid and gas）. The expressions for the amplitude reflection coefficients and the amplitude transmission coefficients are derived by using the potential method. The present derivation is subsequently applied to study the energy conversions among the incident, reflected, and transmitted wave modes. It is found that the reflection and transmission coefficients in the forms of amplitude ratios and energy ratios are functions of the incident angle, the liquid saturation, the frequency of the incident wave, and the elastic constants of the upper and lower media. Numerical results are presented graphically. The effects of the incident angle, the frequency, and the liquid saturation on the amplitude and the energy reflection and transmission coefficients are discussed. It is verified that in the transmission process, there is no energy dissipation at the interface.
基金supported by the National Natural Science Foundation of China (Grant 11372091)China Scholarship Council (Grant 201406120086)
文摘A mathematical formulation is presented for the dynamic stress intensity factor (mode I) of a finite permeable crack subjected to a time-harmonic propagating longitudinal wave in an infinite poroelastic solid. In particular, the effect of the wave-induced fluid flow due to the presence of a liquid-saturated crack on the dynamic stress intensity factor is analyzed. Fourier sine and cosine integral transforms in conjunction with Helmholtz potential theory are used to formulate the mixed boundary-value problem as dual integral equations in the frequency domain. The dual integral equations are reduced to a Fredholm integral equation of the second kind. It is found that the stress intensity factor monotonically decreases with increasing frequency, decreasing the fastest when the crack width and the slow wave wavelength are of the same order. The characteristic frequency at which the stress intensity factor decays the fastest shifts to higher frequency values when the crack width decreases.
文摘Geological sequestration (GS) of carbon dioxide (CO2) is considered as one of the most promising technologies to reduce the amount of anthropogenic CO2 emission in the atmosphere. To ensure success of CO2 GS, monitoring is essential on ascertaining movement, volumes and locations of injected CO2 in the sequestration reservoir. One technique is to use time-lapsed seismic survey mapping to provide spatial distribution of seismic wave velocity as an indicator of CO2 migration and volumes in a storage reservoir with time. To examine the use of time-lapsed seismic survey mapping as a monitoring tool for CO2 sequestration, this paper presents mathematical and experimental studies of the effects of supercritical CO2 injection on the seismic velocity of sandstone initially saturated with saline water. The mathematical model is based on poroelasticity theory, particularly the application of the Biot-Gassmann substitution theory in the modeling of the acoustic velocity of porous rocks containing two-phase immiscible pore fluids. The experimental study uses a high pressure and high temperature triaxial cell to clarify the seismic response of a sample of Berea sandstone to supercritical CO2 injection under deep saline aquifer conditions. Measured ultrasonic wave velocity changes during CO2 injection in the sandstone sample show the effects of pore fluid distribution in the seismic velocity of porous rocks. CO2 injection was shown to decrease the P-wave velocity with increasing CO2 saturation whereas the S-wave velocity was almost constant. The results confirm that the Biot-Gassmann theory can be used to model the changes in the acoustic P-wave velocity of sandstone containing different mixtures of supercritical CO2 and saline water provided the distribution of the two fluids in the sandstone pore space is accounted for in the calculation of the pore fluid bulk modulus. The empirical relation of Brie et al. for the bulk modulus of mixtures of two-phase immiscible fluids, in combination with the Biot-Gassmann theory, was found to satisfactorily represent the pore-fluid dependent acoustic P-wave velocity of sandstone.
基金National Natural Science Foundation of China(11702199)。
文摘Numerical modeling of seepage-induced consolidation process usually encounters significant uncertainty in the properties of geotechnical materials.Assessing the effect of uncertain parameters on the performance variability of the seepage consolidation model is of critical importance to the simulation and tests of this process.To this end,the uncertainty and sensitivity analyses are performed on a seepage consolidation model in a fractured porous medium using the Bayesian sparse polynomial chaos expansion(SPCE)method.Five uncertain parameters including Young’s modulus,Poisson’s ratio,and the permeability of the porous matrix,the permeability within the fracture,and Biot’s constant are studied.Bayesian SPCE models for displacement,flow velocity magnitude,and fluid pressure at several reference points are constructed to represent the input-output relationship of the numerical model.Based on these SPCE models,the total and first-order Sobol’indices are computed to quantify the contribution of each uncertain input parameter to the uncertainty of model responses.The results show that at different locations of the porous domain,the uncertain parameters show different effects on the output quantities.At the beginning of the seepage consolidation process,the hydraulic parameters make major contributions to the uncertainty of the model responses.As the process progresses,the effect of hydraulic parameters decreases and is gradually surpassed by the mechanical parameters.This work demonstrates the feasibility to apply Bayesian SPCE approach to the uncertainty and sensitivity analyses of seepage-induced consolidation problems and provides guidelines to the numerical modelling and experimental testing of such problems.
基金supported by the Chinese National Science Foundation (51104147)the Fundamental Research Funds for the Central Universities (2011QNA17)+1 种基金National Basic Research Program of China (2010CB226800)State Key Laboratory for Geo-mechanics and Deep Underground Engineering in China
文摘The effect of the fracture distribution on CO 2 injection into coal seams was studied with a heterogeneous model having dual porosity to represent both the primary medium (the coal matrix) and the secondary medium (the fractures) under variable stress conditions. A numerical generation method and a digital image processing method were used to model the heterogeneous fracture distribution in the coal. The model solutions demonstrate that: (1) the fractures are the main channel for gas flow and their distribution has an important impact on the gas injection rate; (2) the fractures only affect the injection rate of CO 2 into the coal but not the final storage amount; (3) when gas is injected into coal the fractures will first expand and then close due to the changing effective stresses and the adsorption induced swelling of coal grains. This fully coupled dual-porosity model with a heterogeneous fracture distribution provides a way to predict the CO 2 injection into a coal seam.
基金Project supported by the National Natural Science Foundation of China (No. 11032008)
文摘An extended and reasonable stress boundary condition at an osteon exte- rior wall is presented to solve the model proposed by Remond and Naili. The obtained pressure and fluid velocity solutions are used to investigate the osteonal poroelastic behaviors. The following results are obtained. （i） Both the fluid pressure and the velocity amplitudes are proportional to the strain amplitude and the loading frequency. （ii） In the physiological loading state, the key role governing the poroelastic behaviors of the osteon is the strain rate. （iii） At the osteon scale, the pressure is strongly affected by the permeability variations, whereas the fluid velocity is not.
基金the sponsorship of National Natural Science Foundation of China (41974119, 42030103)Science Foundation from Innovation and Technology Support Program for Young Scientists in Colleges of Shandong Province and Ministry of Science and Technology of China (2019RA2136)Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology (Qingdao) (2021QNLM020001-6)。
文摘Pre-stack seismic inversion is an important method for fluid identification and reservoir characterization in exploration geophysics. In this study, an effective fluid factor is initially established based on Biot poroelastic theory, and a pre-stack seismic inversion method based on Bayesian framework is used to implement the fluid identification. Compared with conventional elastic parameters, fluid factors are more sensitive to oil and gas. However, the coupling effect between rock porosity and fluid content is not considered in conventional fluid factors, which may lead to fuzzy fluid identification results. In addition,existing fluid factors do not adequately consider the physical mechanisms of fluid content, such as squirt flow between cracks and pores. Therefore, we propose a squirt fluid factor(SFF) that minimizes the fluid and pore mixing effects and takes into account the squirt flow. On this basis, a novel P-wave reflection coefficient equation is derived, and the squirt fluid factor is estimated by amplitude variation with offset(AVO) inversion method. The new reflection coefficient equation has sufficient accuracy and can be utilized to estimate the parameters. The effectiveness and superiority of the proposed method in fluid identification are verified by the synthetic and field examples.
