3D scattering of obliquely incident plane SV waves by an alluvial valley embedded in a fluid-saturated, poroelastic layered half-space

The indirect boundary element method is used to study the 3D dynamic response of an infinitely long alluvial valley embedded in a saturated layered half-space for obli- quely incident SV waves. A wave-number transform is first applied along the valley＇s axis to reduce a 3D problem to a 2D plane strain problem. The problem is then solved in the section perpendicular to the axis of the valley. Finally, the 3D dynamic responses of the valley are obtained by an inverse wave-number transform. The validity of the method is con- firmed by comparison with relevant results. The differences between the responses around the valley embedded in dry and in saturated poroelastic medium are studied, and the effects of drainage conditions, porosity, soil layer stiffness, and soil layer thickness on the dynamic response are dis- cussed in detail resulting in some conclusions.

Dynamic Stress and Pore Pressure Around Circular Cavity in Saturated Porous Elastic Half-space Under Harmonic Plane Dilatational Waves

Dynamic stress concentration and pore pressure concentration around an infinitely long cylindrical cavity of circular cross-section subjected to harmonic plane dilatational waves in fluid-saturated porous elastic half-space were obtained by a complex function method based on potential function and multi-polar coordinate. The steady state Biot’s dynamic field equations of porous elastic solid with a viscous liquid were uncoupled into Helmholtz equations via given potential functions. A circular cavity with large radius is used to replace the straight boundary of the saturated porous elastic half-space. The stresses and pore pressures were obtained by using complex functions in multi-polar coordinates with certain boundary conditions of the solid matrix and the fluid matrix. The approximate solutions were compared to existing numerical solutions. Then the variations of the coefficients of dynamic stress concentration and the pore pressures concentration on boundaries of the cavity were discussed with different parameter conditions. The results of the given numerical example indicate that the method used is useful and efficient to the scattering and dynamic stress concentration of plane dilatational waves in saturated porous elastic half-space.

Vertical vibrations of elastic foundation resting on saturated half-space

This paper is mainly concerned with the dynamic response of an elastic foundation of finite height bounded to the surface of a saturated half-space. The foundation is subjected to time-harmonic vertical loadings. First, the transform solutions for the governing equations of the saturated media are obtained. Then, based on the assumption that the contact between the foundation and the half-space is fully relaxed and the halfspace is completely pervious or impervious, this dynamic mixed boundary-value problem can lead to dual integral equations, which can be further reduced to the Predhohn integral equations of the second kind and solved by numerical procedures. In the numerical extortples, the dynamic colnpliances, displacements and pore pressure are developed for a wide range of frequencies and material/geometrical properties of the saturated soil-foundation system. In most of the cases, the dynamic behavior of an elastic foundation resting on the saturated media significantly differs from that of a rigid disc on the saturated half-space. The solutions obtained can be used to study a variety of wave propagation problems and dynamic soil-structure interactions.

The reflection and transmission of elastic waves at a interface between a transversely isotropi celastic soils and saturated soils

According to Biot′s wave equation of transversely isotropic saturated soil, this paper deduces the general equation of the reflection coefficients and transmission coefficients when qP 1 wave goes through from saturated soil to elastic media. The effects of anisotropies and boundary drainage condition on reflection coefficients and transmission coefficients are analyzed by numerical method. The idea of this paper can be applied to the case when qSV wave or qP 2 wave goes through from saturated soils to elastic soils.

Non-axisymmetrical vibration of elastic circular plate on layered transversely isotropic saturated ground

The non-axisymmetrical vibration of elastic circular plate resting on a layered transversely isotropic saturated ground was studied. First, the 3-d dynamic equations in cylindrical coordinate for transversely isotropic saturated soils were transformed into a group of governing differential equations with 1-order by the technique of Fourier expanding with respect to azimuth, and the state equation is established by Hankel integral transform method, furthermore the transfer matrixes within layered media are derived based on the solutions of the state equation. Secondly, by the transfer matrixes, the general solutions of dynamic response for layered transversely isotropic saturated ground excited by an arbitrary harmonic force were established under the boundary conditions, drainage conditions on the surface of ground as well as the contact conditions. Thirdly, the problem was led to a pair of dual integral equations describing the mixed boundaryvalue problem which can be reduced to the Fredholm integral equations of the second kind solved by numerical procedure easily. At the end of this paper, a numerical result concerning vertical and radical displacements both the surface of saturated ground and plate is evaluated.

Scattering of elastic wave by a cylindrical shell deeply embeded in saturated soils

The compression of soil grain and pore fluid as well as viscid coupling of pore fluid and soil skeleton is considered, the scattering problem of incident plane P1 wave （fast compressional wave） by an infinite cylindrical shell deeply embedded in isotropic saturated soils is studied by adopting the amended Biot model, amplitude equations about potential functions of scattering and refracting fields are obtained, and the effect of dimensionless frequencies and shell thickness on the back-scattering spectra and dynamic stress concentration factors of two types of cylindrical shells with high and low rigidity are numerically computed and analyzed.

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 ILLUSTRATION CALCULATIONS OF SECOND ORDER EFFECTS IN ELASTIC HALF-SPACE ACTED UPON BY A NON-UNIFORM SHEAR LOAD

This paper is a continuation of [1]. An example is discussed in derail to illustrate the second order effects. Numerical calculations for the second order elastic material for the z-direction displacement and the stress t(rz) are carried out. It is found that the second order effect is to reduce z-direction displacement and to decrease t(rz)inside the circle but to increase its value outside the circle.

THE GENERAL SOLUTION OF SECOND ORDER EFFECTS IN ELASTIC HALF-SPACE ACTED UPON BY A NON-UNIFORM SHEAR LOAD

This paper is a continuation of [1]. A closed form solution to the second order elasticity problem, when an isotropic compressible elastic half-space undergoes a deformation owing to a non-uniformly distributed shear load, is presented. The method of integral transform is employed to determine the solutions.

DYNAMIC INTERACTION BETWEEN ELASTIC THICK CIRCULAR PLATE AND TRANSVERSELY ISOTROPIC SATURATED SOIL GROUND

A study of the dynamic interaction between foundation and the underlying soil has been presented in a recent paper based on the assumption of saturated ground and elastic circular plate excited by the axisymmetrical harmonic source. However, the assumption may not always be valid. The work is extended to the case of a circular plate resting on transversely isotropic saturated soil and subjected to a non-axisymmetrical harmonic force. The analysis is based on the theory of elastic wave in transversely isotropic saturated poroelastic media established. By the technique of Fourier expansion and Hankel transform, the governing difference equations for transversely isotropic saturated soil are easily solved and the cooresponding Hankel transformed stress and displacement solutions are obtained. Then, under the contact conditions, the problem leads to a pair of dual integral equations which describe the mixed boundary-value problem. Furthermore, the dual integral equations can be reduced to the Fredholm integral equations of the second kind solved by numerical procedure. At the end, a numerical result is presented which indicates that on a certain frequency range, the displacement amplitude of the surface of the foundation increases with the increase of the frequency of the exciting force, and decreases in vibration form with the increase of the distance.

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.

PROPAGATION VELOCITIES OF ELASTIC WAVES IN SATURATED SOILS

Bused on the wave equations established by the authors, the characteristics of propagation velocities of elastic vaves in saturated soils arc analyzed and verified by ultrasonic test in laboratory and seismic survey in the field. The results provide theoretical basis for the determination of physical and mechanical parameters of saturated soils using propagation velocities of elastic waves. especially P-wave Velocity.

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.

Influence of thermal contact resistance on dynamic response of bilayered saturated porous strata

Porous materials can be found in a variety of geophysical and engineering applications.The existence of thermal contact resistance at the interface between bilayered saturated porous strata would result in a significant temperature difference at the interface.An attempt is made to study the thermo-hydro-mechanical coupling dynamic response of bilayered saturated porous strata with thermal contact resistance and elastic wave impedance.The corresponding analytical solutions for the dynamic response of bilayered saturated porous strata under a harmonic thermal load are derived by the operator decomposition method,and their rationality is verified by comparing them with existing solutions.The influences of thermal contact resistance,thermal conductivity ratio,and porosity ratio on the dynamic response of bilayered saturated porous strata are systematically investigated.Outcomes disclose that with the increase of thermal contact resistance,the displacement,pore water pressure and stress decrease gradually,and the temperature jump at the interface between two saturated porous strata increases.

Fluids discrimination by ray-path elastic impedance inversion: A successful case from Sulige tight gas field

Existing seismic prediction methods struggle to effectively discriminate between fluids in tight gas reservoirs,such as those in the Sulige gas field in the Ordos Basin,where porosity and permeability are extremely low and the relationship between gas and water is complicated.In this paper,we have proposed a comprehensive seismic fluid identification method that combines ray-path elastic impedance(REI)inversion with fluid substitution for tight reservoirs.This approach is grounded in geophysical theory,forward modeling,and real data applications.We used geophysics experiments in tight gas reservoirs to determine that Brie's model is better suited to calculate the elastic parameters of mixed fluids than the conventional Wood’s model.This yielded a more reasonable and accurate fluid substitution model for tight gas reservoirs.We developed a forward model and carried out inversion of REI.which reduced the non-uniqueness problem that has plagued elastic impedance inversion in the angle domain.Our well logging forward model in the ray-path domain with different fluid saturations based on a fluid substitution model proved that REI identifies fluids more accurately when the ray parameters are large.The distribution of gas saturation can be distinguished from the crossplot of REI(p=0.10)and porosity.The inverted ray-path elastic impedance profile was further used to predict the porosity and gas saturation profile.Our new method achieved good results in the application of 2D seismic data in the western Sulige gas field.

Prediction of elastic and acoustic behaviors of calcarenite used for construction of historical monuments of Rabat, Morocco

Natural materials(e.g. rocks and soils) are porous media, whose microstructures present a wide diversity.They generally consist of a heterogeneous solid phase and a porous phase which may be fully or partially saturated with one or more fluids. The prediction of elastic and acoustic properties of porous materials is very important in many fields, such as physics of rocks, reservoir geophysics, civil engineering, construction field and study of the behavior of historical monuments. The aim of this work is to predict the elastic and acoustic behaviors of isotropic porous materials of a solid matrix containing dry, saturated and partially saturated spherical pores. For this, a homogenization technique based on the Morie Tanaka model is presented to connect the elastic and acoustic properties to porosity and degree of water saturation. Non-destructive ultrasonic technique is used to determine the elastic properties from measurements of P-wave velocities. The results obtained show the influence of porosity and degree of water saturation on the effective properties. The various predictions of Morie Tanaka model are then compared with experimental results for the elastic and acoustic properties of calcarenite.

SECOND ORDER EFFECTS IN AN ELASTIC HALF－SPACE ACTED UPON BY A NON－UNIFORM NORMAL LOAD

A closed form solution to the second order elasticity problem ,when an isotropiccompressible elastic half-space undergoes a deformation owing to a non-uniformlydistributed normal load,is presented,The method of integral transform is employedand the case when loading is distributed,in accordance with Hertz'x law ,is discussed.The limiting solution for incompressible isotropic elastic material is also derived.Numerical calculations for the second order elastic material for the displacement and the normal stress in the z-direction are carried out .It is found that,in comparison to the linear elastic case,the displacement increases and the normal stress decreases in the second order elastic material.

Laboratory investigation on the mechanical behavior of ice-saturated frozen loess

To investigate the mechanical properties of ice-saturated frozen soil, a series of triaxial tests under various confining pressures（0.5 to 9.0 MPa） on ice-saturated frozen loess with ice content of 23.7% were carried out at a temperature-6 °C, and at 1.25 mm/min of loading rate. The triaxial tests include two loading modes, one with monotonic loading（i.e., triaxial compression）, and another with static cycle loading. The test results under triaxial compression show that the strength and deformation behaviors of ice-saturated frozen loess are affected by confining pressure. According to the test results of triaxial loading-unloading cycle test, the elastic modulus evolution with the number of cycles under different confining pressures are analyzed.

双重孔隙流体饱和介质弹性波散射二维IBIEM模拟

基于平面波势函数,采用间接边界积分方程法(indirect boundary integral equation method, IBIEM)研究了双孔隙流体饱和介质中弹性波入射下二维孔洞的散射特性。推导得到了双孔隙介质中全空间二维线源动力格林函数,并给出了各散射波的位移场和应力场。在数值精度验证的基础上,以双孔隙二维饱和全空间中孔洞为例,解决了平面P、SV波入射下的地震波散射问题。数值结果表明:双重孔隙介质中的位移幅值、环向应力幅值、孔隙压力变化规律与不同入射波形,入射频率,孔隙率和边界排水条件密切相关,位移幅值在低频(无量纲频率η≤2)入射时出现峰值。环向应力幅值与干土条件相比更为复杂,基质孔压与裂缝孔压的存在增大了双重孔隙饱和介质的能量效应,总体震动趋势大于干土条件,环向应力放大可达62%。

快压缩P_(1)波在饱和冻土与弹性基岩分界面上的能量传输特性

基于弹性波在冻结饱和多孔介质与单相弹性介质中的传播理论,选取了饱和冻土中传播速度最大的快压缩P 1波入射在饱和冻土与弹性基岩分界面上的能量传输问题。根据分界面上的边界条件,推导出了快压缩P 1波从饱和冻土介质入射到弹性基岩分界面上透反射振幅比和能量率的解析表达式。研究了快压缩P_(1)波入射在饱和冻土与弹性基岩分界面上的能量与入射角度,入射频率,温度(含冰量),孔隙率,胶结参数以及接触参数的关系。研究结果表明:当入射角度为0°时仅存在压缩波,达到临界角后透射P波消失;各种波在达到临界角时出现不同程度的脉冲,其中反射P_(1)波最为显著;随着胶结参数,孔隙率,接触参数的增大,临界角越早出现;入射频率仅对反射P_(2),P_(3)和S_(2)波的能量反射率影响较大;当温度和含冰量较低或较高时,均不利于反射S_(2)波的能量产生。