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.

The source vector and static displacement field by elastic dislocation on the two-phase saturated medium

Based on the Biot＇s theory about two-phase saturated medium, according to the character of δ function, the Green function on two-phase saturated medium by the point source under concentrated force can be derived. By the Betti＇s theorem for the two-phase saturated medium field, the source vector and static displacement field by elastic dislocation on the two-phase saturated medium were comprehensively discussed.

The numerical simulation for a 3D two-phase anisotropic medium based on BISQ model

Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.

Physical mechanism of seismic attenuation in a two-phase medium

High-frequency seismic attenuation is conventionally attributed to anelastic absorption. In this paper, I present three studies on high-frequency seismic attenuation and propose that the physical mechanism results from the interference of elastic microscopic multiple scattering waves. First, I propose a new theory on wave propagation in a two-phase medium which is based on the concept that the basic unit for wave propagation is a nano- mass point. As a result of the elasticity variations of pore fluid and rock framework, micro multiple scattering waves would emerge at the wavelength of the seismic waves passing through the two-phase medium and their interference and overlap would generate high- frequency seismic attenuation. Second, I present a study of the frequency response of seismic transmitted waves by modeling thin-layers with thicknesses no larger than pore diameters. Results indicate that high-frequency seismic waves attenuate slightly in a near-surface water zone but decay significantly in a near-surface gas zone. Third, I analyze the seismic attenuation characteristics in near-surface water and gas zones using dual-well shots in the Songliao Basin, and demonstrate that the high-frequency seismic waves attenuate slightly in water zones but in gas zones the 160-1600 Hz propagating waves decay significantly. The seismic attenuation characteristics from field observations coincide with the modeling results. Conclusions drawn from these studies theoretically support seismic attenuation recovery.

TIME DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD FOR DYNAMIC ANALYSES IN SATURATED PORO-ELASTO-PLASTIC MEDIUM

A time-discontinuous Galerkin finite element method for dynamic analyses in saturated poro-elasto-plastic medium is proposed.As compared with the existing discontinuous Galerkin finite element methods,the distinct feature of the proposed method is that the continuity of the displacement vector at each discrete time instant is automatically ensured,whereas the discontinuity of the velocity vector at the discrete time levels still remains.The computational cost is then obviously reduced, particularly,for material non-linear problems.Both the implicit and explicit algorithms to solve the derived formulations for material non-linear problems are developed.Numerical results show a good performance of the present method in eliminating spurious numerical oscillations and providing with much more accurate solutions over the traditional Galerkin finite element method using the Newmark algorithm in the time domain.

A Formulation of the Porous Medium Equation with Time-Dependent Porosity: A Priori Estimates and Regularity Results

We consider a generalized form of the porous medium equation where the porosity ϕis a function of time t: ϕ=ϕ(x,t): ∂(ϕS)∂t−∇⋅(k(S)∇S)=Q(S).In many works, the porosity ϕis either assumed to be independent of (or to depend very little of) the time variable t. In this work, we want to study the case where it does depend on t(and xas well). For this purpose, we make a change of unknown function V=ϕSin order to obtain a saturation-like (advection-diffusion) equation. A priori estimates and regularity results are established for the new equation based in part on what is known from the saturation equation, when ϕis independent of the time t. These results are then extended to the full saturation equation with time-dependent porosity ϕ=ϕ(x,t). In this analysis, we make explicitly the dependence of the various constants in the estimates on the porosity ϕby the introduced transport vector w, through the change of unknown function. Also we do not assume zero-flux boundary, but we carry the analysis for the case Q≡0.

Gas hydrate saturation from NGHP 02 LWD data in the Mahanadi Basin

During the Indian National Gas Hydrate Program(NGHP)Expedition 02,Logging-while-drilling(LWD)logs were acquired at three sites(NGHP-02-11,NGHP-02-12,and NGHP-02-13)across the Mahanadi Basin in area A.We applied rock physics theory to available sonic velocity logs to know the distribution of gas hydrate at site NGHP-02-11 and NGHP-02-13.Rock physics modeling using sonic velocity at well location shows that gas hydrate is distributed mainly within the depth intervals of 150-265 m and 100 -215 mbsf at site NGHP-02-11 and NGHP-02-13,respectively,with an average saturation of about 4%of the pore space and the maximum concentration of about 40%of the pore space at 250 m depth at site NGHP-02-11,and at site NGHP-02-13 an average saturation of about 2%of the pore space and the maximum concentration of about 20%of the pore space at 246 m depth,as gas hydrate is distributed mainly within 100-246 mbsf at this site.Saturation of gas hydrate estimated from the electrical resistivity method using density derived porosity and electrical resistivity logs from Archie's empirical formula shows high saturation compared to that from the sonic log.However,estimates of hydrate saturation based on sonic P-wave velocity may differ significantly from that based on resistivity,because gas and hydrate have higher resistivity than conductive pore fluid and sonic P-wave velocity shows strong effect on gas hydrate as a small amount of gas reduces the velocity significantly while increasing velocity due to the presence of hydrate.At site NGHP-02-11,gas hydrate saturation is in the range of 15%e30%,in two zones between 150-180 and 245-265 mbsf.Site NGHP-02-012 shows a gas hydrate saturation of 20%e30%in the zone between 100 and 207 mbsf.Site NGHP-02-13 shows a gas hydrate saturation up to 30%in the zone between 215 and 246 mbsf.Combined observations from rock physics modeling and Archie’s approximation show the gas hydrate concentrations are relatively low(<4%of the pore space)at the sites of the Mahanadi Basin in the turbidite channel system.

Solution of a Green's function for a saturated porous medium in a half-space subjected to a torsional force

Based on the solutions of the Green＇s function for a saturated porous medium obtained by the authors, and using transformation of axisymmetric coordinates, Sommerfeld integrals and superposition of the influence field on a free surface, the authors have obtained displacement solutions of a saturated porous medium subjected to a torsional force in a half-space. The relationship curves of the displacement solutions and various parameters （permeability, frequency, etc.） under action of a unit of torque are also given in this paper. The results are consistent with previous Reissner＇s solutions, where a two-phase medium decays to a single-phase medium. The solution is useful in solving relevant dynamic problems of a two- phase saturated medium in engineering.

On the higher-order thermal vibrations of FG saturated porous cylindrical micro-shells integrated with nanocomposite skins in viscoelastic medium

Since the multi-layered structures are widely used nowadays, and due to interesting applications of cylindrical shells, this study is dedicated to analyzing free vibrational behaviors of functionally graded saturated porous micro cylindrical shells with two nanocomposite skins. Based on Biot's assumptions, constitutive relations for the core are presented and effective properties of the skins are determined via the rule of mixture. A sinusoidal theory is used to capture the shear deformation effects, and to account for the scale effects, the modified couple stress theory is employed which suggests a material length-scale parameter for predicting the results in small-dimension. With the aid of extended form of Hamilton's principle for dynamic systems, differential equations of motion are extracted. Fourier series functions are used to obtain natural frequencies and after validating them, a set of parametric studies are carried out. The results show the significant effects of porosity and Skempton coefficient, pores placement patterns, CNTs addition and distribution patterns, temperature variations, material length-scale parameter and viscoelastic medium on the natural frequencies of the microstructure. The outcomes of this work could be used to design and manufacture more reliable micro cylindrical structures in thermo-dynamical environments.

Pressure/Saturation System for Immiscible Two-Phase Flow: Uniqueness Revisited

We give a sufficient condition for uniqueness for the pressure/saturation system. We establish this condition through analytic arguments, and then construct 'mobilities' (or mobility-like functions) that satisfy the new condition (when the parameter is 2). For the constructed 'mobilities', we do graphical experiments that show, empirically, that this condition could be satisfied for other values of . These empirical experiments indicate that the usual smoothness condition on the fractional flow function (and on the total mobility), for uniqueness and convergence, might not be necessary. This condition is also sufficient for the convergence of a family of perturbed problems to the original pressure/saturation problem.

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.

Effects of Solid Matrix and Porosity of Porous Medium on Heat Transfer of Marangoni Boundary Layer Flow Saturated with Power-Law Nanofluids

The effect of the solid matrix and porosity of the porous medium are first introduced to the study of power-law nanofluids, and the Marangoni boundary layer flow with heat generation is investigated. Two cases of solid matrix of porous medium including glass balls and aluminum foam are considered. The governing partial differential equations are simplified by dimensionless variables and similarity transformations, and are solved numerically by using a shooting method with the fourth-fifth-order Runge-Kutta integration technique. It is indicated that the increase of the porosity leads to the enhancement of heat transfer in the surface of the Marangoni boundary layer flow.

Deformation of a two-phase medium due to a long buried strike-slip fault

The aim of the present paper is to obtain the two-dimensional deformation of a two-phase elastic medium consisting of half-spaces of different ri- gidities in welded contact due to a buried long strike-slip fault. The solution is valid for arbitrary values of the fault-depth and the dip angle. The effect of fault-depth on the displacement and stress fields for different values of dip angle has been studied numerically. It is found that the displacement field varies significantly for a buried fault from the corresponding displacement field for an interface-breaking fault. The contour maps showing the stress field for various dip angles for buried and interface-breaking fault have been plotted. It has been observed that the stress field varies significantly for a buried fault from the corresponding stress field for an interface-breaking fault.

Torsional vibration of a pipe pile in transversely isotropic saturated soil

This study considers the torsional vibration of a pipe pile in a transversely isotropic saturated soil layer. Based on Biot＇s poroelastic theory and the constitutive relations of the transversely isotropic medium, the dynamic governing equations of the outer and inner transversely isotropic saturated soil layers are derived. The Laplace transform is used to solve the governing equations of the outer and inner soil layers. The dynamic torsional response of the pipe pile in the frequency domain is derived utilizing 1D elastic theory and the continuous conditions at the interfaces between the pipe pile and the soils. The time domain solution is obtained by Fourier inverse transform. A parametric study is conducted to demonstrate the influence of the anisotropies of the outer and inner soil on the torsional dynamic response of the pipe pile.

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.

Axisymmetrical analytical solution for vertical vibration of end-bearing pile in saturated viscoelastic soil layer

Based on elasticity and the theory of saturated porous media, and regarding the pile and the soil as a single phase elastic and a saturated viscoelastic media, respectively, the dynamical behavior of vertical vibration of an end-bearing pile in a saturated viscoelastic soil layer is investigated in the frequency domain using the Helmholtz decomposition and variable separation method. The axisymmetrical analytical solutions for vertical vibrations of the pile in a saturated viscoelastic soil layer are obtained, and the analytical expression of the dynamical complex stiffness of the pile top is presented. Responses of dynamic stiffness factor and equivalent damping of pile top with respect to the frequency are shown in figures using a numerical method. Effects of the saturated soil parameters, modulus ratio of the pile to soil, slenderness ratio of pile and pile＇s Poisson ratio, etc. on the stiffness factor and damping are examined. It is shown that, due to the effect of the transversal deformation of the pile and the radial force of the saturated viscoelastic soil acting on the pile, the dynamic stiffness factor and the damping derived from the axisymmetrical solution are greatly different from those derived from the classical Euler-Bernoulli rod model, especially at some specific excitation frequencies. Therefore, there are limitations on applicability of the Euler-Bernoulli rod model in analyzing verticai vibration of the pile. More accurate analysis should be based on a three-dimensional model.

3-D DYNAMIC RESPONSE OF TRANSVERSELY ISOTROPIC SATURATED SOILS

A study on dynamic response of transversely isotropic saturated poroelastic media under a circular non-axisymmetrical harmonic source has been presented by Huang Yi et al. using the technique of Fourier expansion and Hankel transform. However, the method may not always be valid. The work is extended to the general case being in the rectangular coordinate. The purpose is to study the 3-d dynamic response of transversely isotropic saturated soils under a general source distributing in arbitrary rectangular zoon on the medium surface. Based on Biot＇s theory for fluid- saturated porous media, the 3-d wave motion equations in rectangular coordinate for transversely isotropic saturated poroelastic media were transformed into the two uncoupling governing differential equations of 6-order and 2-order respectively by means of the displacement functions. Then, using the technique of double Fourier transform, the governing differential equations were easily solved. Integral solutions of soil skeleton displacements and pore pressure as well as the total stresses for poroelastic media were obtained. Furthermore, a systematic study on half-space problem in saturated soils was performed. Integral solutions for surface displacements under the general harmonic source distributing on arbitrary surface zone, considering both case of drained surface and undrained surface, were presented.

Exact solutions to one-dimensional transient response of incompressible fluid-saturated single-layer porous media

Based on the Biot theory of porous media, the exact solutions to one- dimensional transient response of incompressible saturated single-layer porous media un- der four types of boundary conditions are developed. In the procedure, a relation between the solid displacement u and the relative displacement ＇w is derived, and the well-posed initial conditions and boundary conditions are proposed. The derivation of the solution for one type of boundary condition is then illustrated in detail. The exact solutions for the other three types of boundary conditions are given directly. The propagation of the compressional wave is investigated through numerical examples. It is verified that only one type of compressional wave exists in the incompressible saturated porous media.

Vertical vibration of a floating pile in a saturated viscoelastic soil layer overlaying bedrock

An axisymmetrical analytical solution is developed to investigate the vertical time-harmonic vibration of a floating pile in a saturated viscoelastic soil layer overlaying bedrock. The soil is described by porous medium model established by Boer, while the pile is described by a beam vibration theory. By using separation theory of differential operator and variables to solve the dynamic governing equations for the soil, the fundamental solutions for the soil reactions on side and bottom of the pile are obtained. The dynamic impedance of the pile head is then derived by solving the vibration equation for the pile according to the compatibility condition between the pile and the soil. The proposed model is validated by comparing special cases of our model with the existing results. Numerical examples are presented to analyze the vibration characteristics of the pile.