GREEN FUNCTION ON TWO-PHASE SATURATED MEDIUM BY CONCENTRATED FORCE IN TWO-DIMENSIONAL DISPLACEMENT FIELD

The Green function on two-phase saturated medium by concentrated force has a broad and important use In seismology, seismic engineering, soil mechanics, geophysics, dynamic foundation theory and so on. According to the Green function on two-phase saturated medium by concentrated force in three-dimentional displacement field obtained by Ding Bo-yang et al., it gives out the Green function in two-dimensional displacement field by infinite integral method along x(3)-direction derived by De Hoop and Manolis. The method adopted in the thesis is simpler. The result will be simplified to the boundary element method of dynamic problem.

Solutions of Green's function for Lamb's problem of a two-phase saturated medium

The solutions of Green’s function are significant for simplification of problem on a two-phase saturated medium.Using transformation of axisymmetric cylindrical coordinate and Sommerfeld’s integral,superposition of the influence field on a free surface,authors obtained the solutions of a two-phase saturated medium subjected to a concentrated force on the semi-space.

Digital image processing of saturation for two-phase flow in planar porous media model

In this paper, the accuracy of estimating stained non-wetting phase saturation using digital image processing is examined, and a novel post-processing approach for calculating threshold is presented. In order to remove the effect of the background noise of images and to enhance the high-frequency component of the original image, image smoothing and image sharpening methods are introduced. Depending on the correct threshold, the image binarization processing is particularly useful for estimating stained non-wetting phase saturation. Calculated saturation data are compared with the measured saturation data during the two-phase flow experiment in an artificial steel planar porous media model. The results show that the calculated saturation data agree with the measured ones. With the help of an artificial steel planar porous media model, digital image processing is an accurate and simple method for obtaining the stained non-wetting phase saturation.

AVO forwarding modeling in two-phase media: multiconstrained matrix mineral modulus inversion

AVO forward modeling is based on two-phase medium theory and is considered an effective method for describing reservoir rocks and fluids. However, the method depends on the input matrix mineral bulk modulus and the rationality of the two-phase medium model. We used the matrix mineral bulk modulus inversion method and multiple constraints to obtain a two-phase medium model with physical meaning. The proposed method guarantees the reliability of the obtained AVO characteristicsin two-phase media. By the comparative analysis of different lithology of the core sample, the advantages and accuracy of the inversion method can be illustrated. Also, the inversion method can be applied in LH area, and the AVO characteristics can be obtained when the porosity, fluid saturation, and other important lithology parameters are changed. In particular, the reflection coefficient amplitude difference between the fast P wave and S wave as a function of porosity at the same incidence angle, and the difference in the incidence angle threshold can be used to decipher porosity.

Reflection and transmission of bottom simulating reflectors in gas hydrate-bearing sediments: Two-phase media models

The bottom simulating reflector （BSR） in gas hydrate-bearing sediments is a physical interface which is composed of solid, gas, and liquid and is influenced by temperature and pressure. Deep sea floor sediment is a porous, unconsolidated, fluid saturated media. Therefore, the reflection and transmission coefficients computed by the Zoeppritz equation based on elastic media do not match reality. In this paper, a two-phase media model is applied to study the reflection and transmission at the bottom simulating reflector in order to find an accurate wave propagation energy distribution and the relationship between reflection and transmission and fluid saturation on the BSR. The numerical experiments show that the type I compressional （fast） and shear waves are not sensitive to frequency variation and the velocities change slowly over the whole frequency range. However, type II compressional （slow） waves are more sensitive to frequency variation and the velocities change over a large range. We find that reflection and transmission coefficients change with the amount of hydrate and free gas. Frequency, pore fluid saturation, and incident angle have different impacts on the reflection and transmission coefficients. We can use these characteristics to estimate gas hydrate saturation or detect lithological variations in the gas hydrate-bearing sediments.

A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

Numerical simulation of two-phase flow in fractured porous media using streamline simulation and IMPES methods and comparing results with a commercial software

Streamline simulation is developed to simulate waterflooding in fractured reservoirs. Conventional reservoir simulation methods for fluid flow simulation in large and complex reservoirs are very costly and time consuming. In streamline method, transport equations are solved on one-dimensional streamlines to reduce the computation time with less memory for simulation. First, pressure equation is solved on an Eulerian grid and streamlines are traced. Defining the "time of flight", saturation equations are mapped and solved on streamlines. Finally, the results are mapped back on Eulerian grid and the process is repeated until the simulation end time. The waterflooding process is considered in a fractured reservoir using the dual porosity model. Afterwards, a computational code is developed to solve the same problem by the IMPES method and the results of streamline simulation are compared to those of the IMPES and a commercial software. Finally, the accuracy and efficiency of streamline simulator for simulation of two-phase flow in fractured reservoirs has been proved.

GURTIN VARIATIONAL PRINCIPLE AND FINITE ELEMENT SIMULATION FOR DYNAMICAL PROBLEMS OF FLUID-SATURATED POROUS MEDIA

Based on the theory of porous media, a general Gurtin variational principle for the initial boundary value problem of dynamical response of fluid-saturated elastic porous media is developed by assuming infinitesimal deformation and incompressible constituents of the solid and fluid phase. The finite element formulation based on this variational principle is also derived. As the functional of the variational principle is a spatial integral of the convolution formulation, the general finite element discretization in space results in symmetrical differential-integral equations in the time domain. In some situations, the differential-integral equations can be reduced to symmetrical differential equations and, as a numerical example, it is employed to analyze the reflection of one-dimensional longitudinal wave in a fluid-saturated porous solid. The numerical results can provide further understanding of the wave propagation in porous media.

THE NON-AXISYMMETRICAL DYNAMIC RESPONSE OF TRANSVERSELY ISOTROPIC SATURATED POROELASTIC MEDIA

The Blot's wave equations of transversely isotropic saturated poroelastic media excited hy non-axisymmetrical harmonic source were solved by means of Fourier expansion and Hankel transform. Then the components of total stress in porous media are expressed with the solutions of Biot's wave equations. The method of research on non-axisymmetrical dynamic response of saturated porous media is discussed, and a numerical result is presented.

A control volume based finite element method for simulating incompressible two-phase flow in heterogeneous porous media and its application to reservoir engineering

Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.

Response of saturated porous media subjected to local thermal loading on the surface of semi-infinite space

Heat source function method is adopted in the present paper to derive elementary solutions of coupled thermo-hydro-mechanical consolidation for saturated porous media under conjunct actions of instantaneous point heat source, instantaneous point fluid source and constant volume force. By using the so-called fictitious heat source method and images method, the solutions of a semi-infinite saturated porous medium subjected to a local heat source with time-varied intensity on its free surface are developed from elementary solutions. The numerical integral methods for calculating the unsteady temperature, pore pressure and displacement fields are given. The thermomechanical response are analyzed for the case of a circular planar heat source. Besides, the thermal consolidation characteristics of a saturated porous medium subjected to a harmonic thermal loading are also given, and the fluctuation processes of the field variables located below the center of heat source are analyzed.

THE NON-AXISYMMETRICAL DYNAMIC RESPONSE OF TRANSVERSELY ISOTROPIC SATURATED POROELASTIC MEDIA

The Biot’s wave equations of transversely isotropic saturated poroelastic media excited by non_axisymmetrical harmonic source were solved by means of Fourier expansion and Hankel transform. Then the components of total stress in porous media are expressed with the solutions of Biot’s wave equations. The method of research on non_axisymmetrical dynamic response of saturated porous media is discussed, and a numerical result is presented.

A lattice Boltzmann exploration of two-phase displacement in 2D porous media under various pressure boundary conditions

While experimental designs developed in recent decades have contributed to research on dynamic nonequilibrium effects in transient two-phase flow in porous media,this problem has been seldom investigated using direct numerical simulation(DNS).Only a few studies have sought to numerically solve Navier—Stokes equations with level-set(LS)or volume-of-fluid(VoF)methods,each of which has constraints in terms of meniscus dynamics for various flow velocities in the control volume(CV)domain.The Shan—Chen multiphase multicomponent lattice Boltzmann method(SC-LBM)has a fundamental mechanism to separate immiscible fluid phases in the density domain without these limitations.Therefore,this study applied it to explore two-phase displacement in a single representative elementary volume(REV)of two-dimensional(2D)porous media.As a continuation of a previous investigation into one-step inflow/outflow in 2D porous media,this work seeks to identify dynamic nonequilibrium effects on capillary pressure—saturation relationship(P_(c)—S)for quasi-steady-state flow and multistep inflow/outflow under various pressure boundary conditions.The simulation outcomes show that P_(c),S and specific interfacial area(a_(nw))had multistep-wise dynamic effects corresponding to the multistep-wise pressure boundary conditions.With finer adjustments to the increase in pressure over more steps,dynamic nonequilibrium effects were significantly alleviated and even finally disappeared to achieve quasisteady-state inflow/outflow conditions.Furthermore,triangular wave-formed pressure boundary conditions were applied in different periods to investigate dynamic nonequilibrium effects for hysteretical Pc—S.The results showed overshoot and undershoot of P_(c)to S in loops of the nonequilibrium hysteresis.In addition,the flow regimes of multistep-wise dynamic effects were analyzed in terms of Reynolds and capillary numbers(Re and Ca).The analysis of REV-scale flow regimes showed higher Re(1<Re<10)for more significant dynamic nonequilibrium effects.This indicates that inertia is critical for transient twophase flow in porous media under dynamic nonequilibrium conditions.

A wavelet finite-difference method for numerical simulation of wave propagation in fluid-saturated porous media

In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.

Finite element equations and numerical simulation of elastic wave propagation in two-phase anisotropic media

Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical solution of finite element equations is given.Finally,Properties of elastic wave propagation are observed and analyzed through FEM modeling.

MATHEMATICAL MODEL OF TWO-PHASE FLUID NONLINEAR FLOW IN LOW-PERMEABILITY POROUS MEDIA WITH APPLICATIONS

A mathematical model of two-phase fluid nonlinear flow in the direction of normal of ellipse through low-permeability porous media was established according to a nonlinear flow law expressed in a continuous function with three parameters, a mass conservation law and a concept of turbulent ellipses. A solution to the model was obtained by using a finite difference method and an extrapolation method. Formulas of calculating development index not only before but also after water breaks through an oil well in the condition of two-phase fluid nonlinear flow in the media were derived. An example was discussed. Water saturation distribution was presented. The moving law of drainage front was found. Laws of change of pressure difference with time were recognized. Results show that there is much difference of water saturation distribution between nonlinear flow and linear flow; that drainage front by water moves faster, water breaks through sooner and the index gets worse because of the nonlinear flow; and that dimensionless pressure difference gets larger at the same dimensionless time and difficulty of oil development becomes bigger by the nonlinear flow. Thus, it is necessary that influence of nonlinear flow on development indexes of the oil fields be taken into account. The results provide water-flooding development of the oilfields with scientific basis.

Wave field in two-phase media by the convolutional differentiator method

This paper applies the convolutional differentiator method, based on generalized Forsyte orthogonal polynomial （CFPD）, to simulate the seismic wave propagation in two-phase media. From the numerical results we can see that three types of waves, fast P-waves, S-waves and slow P-waves, can be observed in the seismic wave field. The experiments on anisotropic models demonstrate that the wavefront is elliptic instead of circular and S-wave splitting occurs in anisotropic two-phase media. The research has confirmed that the rules of elastic wave propagation in fluid-saturated porous media are controlled by Biot＇s theory. Experiment on a layered fault model shows the wavefield generated by the interface and the fault very well, indicating the effectiveness of CFPD method on the wavefield modeling for real layered media in the Earth. This research has potential applications to the investigation of Earth＇s deep structure and oil/gas exploration.

PARAMETERS INVERSION OF FLUID-SATURATED POROUS MEDIA BASED ON NEURAL NETWORKS

The multi- layers feedforward neural network is used for inversion ofmaterial constants of fluid-saturated porous media. The direct analysis of fluid-saturated porousmedia is carried out with the boundary element method. The dynamic displacement responses obtainedfrom direct analysis for prescribed material parameters constitute the sample sets training neuralnetwork. By virtue of the effective L-M training algorithm and the Tikhonov regularization method aswell as the GCV method for an appropriate selection of regu-larization parameter, the inversemapping from dynamic displacement responses to material constants is performed. Numerical examplesdemonstrate the validity of the neural network method.

Petroleum Hydrocarbon Transport through Saturated-unsaturated Media: A Numerical Model and Practice

The mathematical model of migration of total petroleum hydrocarbons in unsaturated media was described,including convection,molecular diffusion,mechanical dispersion and adsorption,and chemical reactions.By finite element method,a numerical model of evaluating petroleum hydrocarbon migration through contaminated soils was created and applied to the environmental investigations of a relocated mechanical factory in Shanghai.The model consisted of three compacted soil layers:plain fill,sandy silt and silty clay.The results showed that pollutants in the sandy silt traveled faster than that in the plain fill and silty clay.The same decreasing trend of migration velocity was observed in all of the three soil layers.After 180 d,the concentrations of pollutants in the sandy silt can be as low as 40% of the original maximum,while its counterpart in the silty clay is 64%.

Vibrating Liquefaction Experiment and Mechanism Study in Saturated Granular Media

By the vibrating liquefaction experiment of tailings and fine-ores of iron, it is observed and noted that the change of pore water pressure when the vibrating liquefaction takes place. Based on relevant suppositions, the equation of wave propagation in saturated granular media is obtained. This paper postulates the potential vector equation and the velocity expression of three kinds of body waves under normal conditions. Utilizing the wave theory and the experimental results, the influence of three body waves on pore water pressure and granules has been analyzed in detail. This revealed the rapid increment mechanism of pore water pressure and the wave mechanism of vibrating liquefaction.