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 th...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.展开更多
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 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.展开更多
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...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.展开更多
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 wa...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.展开更多
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 fro...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.展开更多
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 fea...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.展开更多
Conversion of labile Pb species into chloropyromorphite(CPY)using phosphorus-bearing amendments was considered to be an ideal strategy in soil passivation remediation.However,the fate and transport of CPY in the soil ...Conversion of labile Pb species into chloropyromorphite(CPY)using phosphorus-bearing amendments was considered to be an ideal strategy in soil passivation remediation.However,the fate and transport of CPY in the soil are poorly understood.This study aims to fill the knowledge gap by evaluating the fate and transport of CPY under environmentally relevant conditions of humic acid(HA),pH,electrolyte concentration,and species through the saturated sandy medium.Results showed that bare CPY colloids are basically immobile in sandy porous media while the co-existence of HA made the transport of CPY improved by30%-93.5%.Facilitated transport of CPY was attributed to the increased stability of CPY and the repulsive interaction between CPY particles and sands due to HA adsorption.The mobility of CPY was also increased with increasing pH from 5.0 to 9.0.When the pH was 9 with a10 mmol/L NaCl background solution,the stronger energy barrier between CPY and sand led to enhanced transport behavior.The divalent Ca^(2+)had a more dramatic effect than monovalent Na^(+)on the aggregation and sedimentation of CPY colloids due to its effectivescreening of the surface charge of CPY and bridging interaction with CPY particles.Derjaguin-LandauVerwey-Overbeek theory and attachment efficiency calculation indicated that high energy barriers were responsible for the high mobility of CPY colloids,while the retention of CPY in sands was mainly caused by secondary energy minimum and physically straining.The findings of this work can help to evaluate the fate of soil passivation remediation products in natural water and soil.展开更多
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 fiel...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.展开更多
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 ...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.展开更多
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 satis...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.展开更多
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 remov...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.展开更多
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 soli...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.展开更多
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...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.展开更多
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, th...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.展开更多
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 consum...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.展开更多
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 vibratio...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.展开更多
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 ...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.展开更多
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 th...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.展开更多
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 medi...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.展开更多
This paper examines the two-phase flow for a horizontal well penetrating a naturally fractured reservoir with edge water injection by means of a fixed streamlinemodel,The mathematical model of the verical two-dimensio...This paper examines the two-phase flow for a horizontal well penetrating a naturally fractured reservoir with edge water injection by means of a fixed streamlinemodel,The mathematical model of the verical two-dimensional flow or oil-water for a horizontal well in a medium with double-porosity is established ,and whose accuratesolutions are obtained by using the characteristic method .The saturation distributionsin the fractured system and the matrix system as well as the formula of the time of water free production are presented .All these results provide a theoretical basis and a computing method for oil displacement by edge water from naturally fracturedreservirs.展开更多
文摘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.
基金supported by the National Natural Science Foundation of China (10572129)
文摘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.
文摘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.
基金National Natural Science Foundation (Project number 40604013).
文摘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.
文摘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.
基金The project supported by the National Natural Science Foundation of China(19832010,50278012,10272027)the National Key Basic Research and Development Program(973 Program,2002CB412709)
文摘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.
基金supported by the National Natural Science Foundation of China(No.41931292)the Priority Academic Program Development of Jiangsu Higher Education Institutions(No.164320H101)。
文摘Conversion of labile Pb species into chloropyromorphite(CPY)using phosphorus-bearing amendments was considered to be an ideal strategy in soil passivation remediation.However,the fate and transport of CPY in the soil are poorly understood.This study aims to fill the knowledge gap by evaluating the fate and transport of CPY under environmentally relevant conditions of humic acid(HA),pH,electrolyte concentration,and species through the saturated sandy medium.Results showed that bare CPY colloids are basically immobile in sandy porous media while the co-existence of HA made the transport of CPY improved by30%-93.5%.Facilitated transport of CPY was attributed to the increased stability of CPY and the repulsive interaction between CPY particles and sands due to HA adsorption.The mobility of CPY was also increased with increasing pH from 5.0 to 9.0.When the pH was 9 with a10 mmol/L NaCl background solution,the stronger energy barrier between CPY and sand led to enhanced transport behavior.The divalent Ca^(2+)had a more dramatic effect than monovalent Na^(+)on the aggregation and sedimentation of CPY colloids due to its effectivescreening of the surface charge of CPY and bridging interaction with CPY particles.Derjaguin-LandauVerwey-Overbeek theory and attachment efficiency calculation indicated that high energy barriers were responsible for the high mobility of CPY colloids,while the retention of CPY in sands was mainly caused by secondary energy minimum and physically straining.The findings of this work can help to evaluate the fate of soil passivation remediation products in natural water and soil.
基金National Natural Science Foundation of China Under Grant No.11172268
文摘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.
文摘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.
文摘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.
基金supported by the National Natural Science Foundation of China(Grant No51079043)the Special Fund for Public Welfare Industry of Ministry of Water Resources of China(Grants No200901064 and 201001020)the Research Innovation Program for College Graduates of Jiangsu Province(Grant No CXZZ11_0450)
文摘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.
基金Supported by the National Natural Science Foundation of China under Grant No 51305080
文摘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.
文摘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.
基金The 111 Project under Grant No.B13024the National Natural Science Foundation of China under Grant Nos.U1134207 and 51378177the Program for New Century Excellent Talents in University under Grant No.NCET-12-0843
文摘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.
文摘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.
基金Project supported by the National Natural Science Foundation of China (No. 10872124)
文摘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.
文摘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.
基金Project supported by the Earthquake Administration Foundation for Seismological Researches of China(No.200808022)the National Natural Science Foundation of China(Nos.50778163 and 50708095)the National Basic Research Program of China(No.2007CB714200)
文摘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.
基金Projects(50809009,51578100) supported by the National Natural Science Foundation of ChinaProjects(3132014326,3132015095) supported by the Fundamental Research Funds for the Central Universities of China
文摘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.
文摘This paper examines the two-phase flow for a horizontal well penetrating a naturally fractured reservoir with edge water injection by means of a fixed streamlinemodel,The mathematical model of the verical two-dimensional flow or oil-water for a horizontal well in a medium with double-porosity is established ,and whose accuratesolutions are obtained by using the characteristic method .The saturation distributionsin the fractured system and the matrix system as well as the formula of the time of water free production are presented .All these results provide a theoretical basis and a computing method for oil displacement by edge water from naturally fracturedreservirs.
