淤砂层模拟方法对上层理想流体动力反应的影响——平面P波入射

采用解析法,以弹性半空间-淤砂层-理想流体层为对象,分别将淤砂层模拟为固体介质和流体介质,分析了平面P波入射时上层理想流体的动力反应特性。通过算例分析,比较了将淤砂层模拟为弹性固体、黏弹性固体、理想流体和黏性流体介质时与将淤砂层模拟为两相多孔介质时计算结果的差异。研究表明当渗透系数很小时,多孔介质可以近似简化为黏弹性固体,甚至线弹性固体。这一分析研究对于在高坝-库水-淤砂-地基系统地震反应分析的数值模型中采用合理的淤砂层模型,提高计算效率具有重要参考意义,而且还可以用于校核数值模型的合理性。     

Steady‐state analytical solutions of flow and deformation coupling due to a point sink within a finite fluid‐saturated poroelastic layer

An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.     

Poroelastic model for pile–soil interaction in a half‐space porous medium due to seismic waves

In this paper, frequency domain dynamic response of a pile embedded in a half‐space porous medium and subjected to P, SV seismic waves is investigated. According to the fictitious pile methodology, the problem is decomposed into an extended poroelastic half‐space and a fictitious pile. The extended porous half‐space is described by Biot's theory, while the fictitious pile is treated as a bar and a beam and described by the conventional 1‐D structure vibration theory. Using the Hankel transformation method, the fundamental solutions for a half‐space porous medium subjected to a vertical or a horizontal circular patch load are established. Based on the obtained fundamental solutions and free wave fields, the second kind of Fredholm integral equations describing the vertical and the horizontal interaction between the pile and the poroelastic half‐space are established. Solution of the integral equations yields the dynamic response of the pile to plane P, SV waves. Numerical results show the parameters of the porous medium, the pile and incident waves have direct influences on the dynamic response of the pile–half‐space system. Significant differences between conventional single‐phase elastic model and the poroelastic model for the surrounding medium of the pile are found. Copyright © 2007 John Wiley & Sons, Ltd.     

Lateral dynamic response of a pipe pile in saturated soil layer

This paper presents an analytical solution for the lateral dynamic response of a pipe pile in a saturated soil layer. The wave propagations in the saturated soil and the pipe pile are simulated by Biot's three‐dimensional poroelastic theory and one‐dimensional elastic theory, respectively. The governing equations of soil are solved directly without introducing potential functions. The displacement response and dynamic impedances of the pipe pile are obtained based on the continuous conditions between the pipe pile and both the outer and inner soil. A comparison with an existing solution is performed to verify the proposed solution. Selected numerical results for the lateral dynamic responses and impedances of the pipe pile are presented to reveal the lateral vibration characteristics of the pile‐soil system. Copyright © 2015 John Wiley & Sons, Ltd.     

Prediction of long-term settlement on soft clay using shear wave velocity and damping characteristics

This paper presents a method of calculating long-term settlement of a loaded pad on soft clay at Klang, Selangor, Malaysia where the soil model is treated as an anelastic material of viscoelastic property. Initially, an elastic shear modulus (G) value from shear wave velocity profiles of the seismic tests from spectral analysis of surface wave (SASW) and continuous surface wave (CSW) tests was obtained. A value of damping (D) at an equivalent elastic strain is then calculated from the hysteresis of the plate load tests corrected to equivalent strain using the Damping–Strain formula. The calculated elastic settlement and its equivalent damping are then used to calculate the long-term settlement by applying the generalised viscoelastic formula. Comparisons to traditional methods of settlement predictions were made and the viscoelastic formula has shown better agreement to the observed settlement. Further modification of the settlement formula is introduced to improve the settlement accuracy to 10%.     

半空间饱和土中圆形衬砌对弹性压缩波的散射

借助于Biot 波动理论和弹性波的传播理论,采用复变函数和多级坐标法,对半空间饱和土中圆形衬砌结构对弹性稳态压缩波的散射问题进行求解和分析。利用一个半径很大的圆弧来逼近半空间直边界,将待解问题转化为稳态弹性压缩波在一个大圆孔和一个弹性衬砌结构的散射问题。通过引入势函数,将饱和土的Biot波动方程和衬砌的弹性波动方程解耦成Helmholtz 方程,借助复变函数级数展开便可以预先写出该组Helmholtz方程的通解。然后,通过引用复变量,把饱和土和衬砌结构中的应力、位移及孔压用设定的势函数表示出来,再利用半空间饱和土和衬砌结构的连续性条件和近似直边界的圆弧边界和衬砌内边界的边界条件求解出该组势函数的特解。最后,利用势函数的特解,得到饱和土中的位移,应力和孔压及衬砌结构的位移和应力;变换不同的参数求解衬砌结构内外边界的动应力和孔压的集中系数,通过对算例结果的分析得出一系列有益的结论。     

Effects of Near-Surface Absorption on Reflection Characteristics of Continental Interbedded Strata: the Dagang Oilfield as an Example

<正>Due to the effects of seismic wave field interference,the reflection events generated from interbedded and superposed sand and shale strata no longer have an explicit corresponding relationship with the geological interface.The absorption of the near-surface layer decreases the resolution of the seismic wavelet,intensifies the interference of seismic reflections from different sand bodies,and makes seismic data interpretation of thin interbedded strata more complex and difficult. In order to concretely investigate and analyze the effects of the near-surface absorption on seismic reflection characteristics of interbedded strata,and to make clear the ability of current technologies to compensate the near-surface absorption,a geological model of continental interbedded strata with near-surface absorption was designed,and the prestack seismic wave field was numerically simulated with wave equations.Then,the simulated wave field was processed by the prestack tune migration, the effects of near-surface absorption on prestack and poststack reflection characteristics were analyzed,and the near-surface absorption was compensated for by inverse Q-filtering.The model test shows that:(1) the reliability of prediction and delineation of a continental reservoir with AVO inversion is degraded due to the lateral variation of the near-surface structure;(2) the corresponding relationships between seismic reflection events and geological interfaces are further weakened as a result of near-surface absorption;and(3) the current technology of absorption compensation probably results in false geological structure and anomaly.Based on the model experiment,the real seismic data of the Dagang Oil Field were analyzed and processed.The seismic reflection characteristics of continental interbedded strata were improved,and the reliability of geological interpretation from seismic data was enhanced.     

Three‐dimensional time‐harmonic fundamental solutions for a fluid‐saturated poroelastic half‐space with partially permeable free surface

By virtue of a pair of scalar potentials for the displacement of the solid skeleton and the pore fluid pressure field of a saturated poroelastic medium, an alternative solution method to the Helmholtz decomposition is developed for the wave propagation problems in the framework of Biot's theory. As an application, a comprehensive solution for three‐dimensional response of an isotropic poroelastic half‐space with a partially permeable hydraulic free surface under an arbitrarily distributed time‐harmonic internal force field and fluid sources is developed. The Green's functions for the poroelastic fields, corresponding to point, ring, and disk loads, are reduced to semi‐infinite complex‐valued integrals that can be evaluated numerically by an appropriate quadrature scheme. Analytical and numerical comparisons are made with existing elastic and poroelastic solutions to illustrate the quality and features of the solution. Copyright © 2016 John Wiley & Sons, Ltd.     

Response of poroelastic halfspace to steady-state harmonic surface tractions

Traditionally, most formulations of dynamic halfspace problems have represented the material as either an elastic or a viscoelastic solid. Herein the counterpart of Lamb's elastodynamic problem is reformulated and solved for a liquid-saturated poroelastic halfspace using Biot's theory of poroelasticity. The responses of the solid and fluid phases are evaluated due to steady-state harmonic concentrated loads applied to each phase at the surface. The solutions are presented over a broad range of permeabilities and are compared to solutions to Lamb's problem for equivalent drained and undrained solids. Methodology is then introduced by which these results are treated as Green functions for the solution of a mixed boundary-value problem. namely, the response of the poroelastic halfspace to steady-state harmonic vertical motion of a rigid. massless plate. It is observed that small differences exist among overall compliance functions for a drained solid, an undrained solid, and a liquid-saturated porous, halfspace. However, use of the poroelastic model permits the distribution between effective skeletal normal stresses and fluid stresses to be determined.     

运动荷载附近有限层厚软土地基的振动研究

基于Biot多孔弹性介质的波动理论,研究了运动荷载附近软土地基的振动问题。假设一条形均布荷载作用在地基表面,则该模型可视为平面应变问题进行分析。通过引入4个势函数和Helmholtz原理,并利用Fourier变换及逆变换技术,获得了运动荷载作用下软土地基的应力、位移和孔隙水压力的解答。利用离散Fourier逆变换得到数值计算结果,分析了荷载速度,频率以及软土的渗透系数及多孔弹性参数对运动荷载作用下地表竖向位移及土体中任一点孔隙水压力分布的影响。     

Dynamic responses of a poroelastic half‐space from moving trains caused by vertical track irregularities

The vibrations of railway tracks on a poroelastic half‐space generated by moving trains are investigated through a vehicle–track–ground coupling model. The theoretical model incorporates a vehicle, a track, and a fully saturated poroelastic half‐space soil medium. The source of vibration excitation is divided into two components: the quasi‐static loads and the dynamic loads. The quasi‐static loads are related to the static component of the axle loads, whereas the dynamic loads are due to the dynamic wheel–rail interaction. A linear Hertizian contact spring is introduced between each wheelset and the rail to consider the dynamic loads. Biot's dynamic theory is used to characterize the poroelastic half‐space soil medium. Using the Fourier transform, the governing equations for the track–ground system are solved and the numerical results are presented for a single axle vehicle model. The different dynamic characteristics of the elastic soil medium and the saturated poroelastic medium are investigated. In addition, the different roles of the moving axle loads and the roughness‐induced dynamic loads are identified. It is concluded that the vibration level of the free field off the track predicted by the poroelastic soil medium is smaller than that predicted by the elastic soil medium for vehicle speed below the Rayleigh wave speed of the poroelastic half‐space, whereas it is larger for vehicle speed above the Rayleigh wave speed. The dynamic loads play an important role in the dynamic responses of the track–ground system. Copyright © 2010 John Wiley & Sons, Ltd.     

淤砂层模型对上层理想流体动力反应的影响——平面SV波入射

采用解析方法,分析了平面SV波斜入射时弹性半空间-淤砂层-理想流体层系统的动力反应,其中将淤砂层分别模拟为线弹性固体、黏弹性固体、理想流体和黏性流体等不同特性的介质,并与将淤砂层模拟为两相多孔介质的计算结果进行了比较。研究表明在平面SV波斜入射的情况下,当渗透系数很小时,多孔介质可以近似简化为黏弹性固体。     

Quasi-static deformation due to two-dimensional seismic sources embedded in an elastic half-space in welded contact with a poroelastic half-space

The Biot linearized theory of fluid saturated porous materials is used to study the plane strain deformation of a two-phase medium consisting of a homogeneous, isotropic, poroelastic half-space in welded contact with a homogeneous, isotropic, perfectly elastic half-space caused by a two-dimensional source in the elastic half-space. The integral expressions for the displacements and stresses in the two half-spaces in welded contact are obtained from the corresponding expressions for an unbounded elastic medium by applying suitable boundary conditions at the interface. The case of a long dip-slip fault is discussed in detail. The integrals for this source are solved analytically for two limiting cases: (i) undrained conditions in the high frequency limit, and (ii) steady state drained conditions as the frequency approaches zero. It has been verified that the solution for the drained case (ω → 0) coincides with the known elastic solution. The drained and undrained displacements and stresses are compared graphically. Diffusion of the pore pressure with time is also studied.     

Surface load dynamic solution of saturated transversely isotropic multilayer half-space

This paper treats the dynamic response of a multilayered transversely isotropic fluid saturated poroelastic half-space under surface time-harmonic traction. The governing system of partial differential equations is uncoupled with the use of a set of physically meaningful and complete potential functions that decompose different body waves in a saturated poroelastic transversely isotropic medium. After expressing the equations in the Hankel-Fourier domain, a proper algebraic factorization is applied to generate reflection and transmission matrices for decomposed waves. All responses including displacements, stresses, and pore fluid pressure for both general patch load and point load are presented in the form of semi-infinite line integrals. The verification of the method is confirmed with the degeneration of the solutions presented here to the existing solutions for dried both homogeneous and multilayered elastic half-spaces as well as poroelastic half-space. Selected numerical results are depicted to investigate the effects of layering and pore pressure on responses of a transversely isotropic poroelastic medium. The load distribution effects are studied by comparison of the patch and point load responses. Also, resonance notion and effective parameters on this phenomenon such as layering system and anisotropy contrast are discussed. Significant influence of materials and layering configuration on number and amplitude of resonances depicted through the numerical evaluation.     

Deformation of a poroelastic layer overlying an elastic half‐space due to dip‐slip faulting

An analytical solution of the plane strain problem of the deformation of a homogeneous, isotropic, poroelastic layer of uniform thickness overlying a homogeneous, isotropic, elastic half‐space due to two‐dimensional seismic sources buried in the elastic half‐space has been obtained. The integral expressions for the displacements, stresses and pore pressure have been obtained using the stress function approach by applying suitable boundary conditions at the free surface and the interface. The solution obtained is in the Laplace–Fourier transform domain. The case of a vertical dip‐slip line dislocation for the oceanic crust model of Earth is studied in detail. Schapery's formula is used for the Laplace inversion and the extended Simpson's formula for the Fourier inversion. Diffusion of pore pressure in the layer is studied numerically. Contour maps showing the pore pressure in the poroelastic layer have been plotted. The effect of the compressibility of the solid and fluid constituents on pore pressure has also been studied. Copyright © 2015 John Wiley & Sons, Ltd.     

Generalized poroelastic wellbore problem

This paper presents a novel analytical solution to the transient, z‐dependent, and asymmetric problem of an infinite wellbore drilled into a fluid‐saturated porous medium. The formulations are based on Biot's linear theory of poroelasticity, in which the dependency of poroelastic field variables to spatial coordinates as well as time domain is considered in the most general form. This gives flexibility to the solution in cases that cannot be analyzed using the conventional plane strain or symmetric models. One such case is when calculating the stress variations around an inclined wellbore where the far‐field stresses are acting over a finite vertical section. The results of our solution to this case with a three‐dimensional state of far‐field stress are used to analyze the stability of inclined wellbores passing through abnormally stressed formations. The presented solution is capable of finding expressions for fundamental solutions with stress or flow boundary conditions at the wellbore. These solutions are here adopted to analyze the pressure disturbances generated by multiprobe formation tester, a standard wireline device that is designed for downhole fluid sampling as well as estimating the directional permeabilities of subsurface earth formations. A comparison with the conventional solution for the relevant pressure diffusion equation indicates that the poroelastic effect is relatively significant in relation to the transient response of the pore pressure. Further, it is shown that the finite dimensions of sink probe would, to a great extent, contribute to the formation's pore pressure variations at its immediate proximity. Copyright © 2013 John Wiley & Sons, Ltd.     

可控震源粘弹性波动方程有限差分模拟

在一般的地震波计算中,地球介质可以当作各向同性的完全弹性体来对待。而对地表土壤和精细观测,粘弹性介质模型比完全弹性模型更符合实际。可控震源是近年来发展起来的精确可控、安全有效的新型震源,但国内对可控震源传播动力学的研究较少。对可控震源激发的地震波在粘弹性介质中的传播,用有限差分方法进行了数值模拟,给出了Kelvin一阶粘弹性波动方程组、吸收边界条件和震源处理。提出了检波器接收信号加速度和传播距离的三次多项式指数关系式,比线性指数关系更准确。使用该公式可以预测采集系统有效接收距离,指导测线布置。      

横观各向同性饱和地基竖向振动的衰减特性

基于提出的横观各向同性饱和多孔介质Biot波动方程的一般解,研究了饱和半空间地基在竖向点源简谐激振荷载作用下地表振动的衰减特征,分析了激振频率以及横观各向同性饱和土介质的各向异性参数和孔隙渗透系数对地表振动特征的影响。计算结果表明,低频和高频激振时,地表位移衰减特性存在明显差异;在饱和土的各向异性参数中,纵向和水平方向动态渗透系数比值和刚度系数比值对地表位移衰减影响最大,这也说明采用各向同性饱和介质的动力学模型不能准确地描述具有明显各向异性特性的饱和土地基的动力特性。     

粘弹性波动方程的系数反演

波动方程的系数反演是一种用于识别地下介质物理力学参数的重要方法。过去这方面的研究一般都建立在弹性模型的基础上,本文则提出了粘弹性反演模型,这种模型更能真实地反映波在地下传播的实际情形。在粘弹性反演模型的基础上,我们还提出了一整套行之有效的数值反演方法,并在频域中完成了反演计算。最后的数值计算结果是令人满意的,证明了本文的模型和方法的合理性。     

2008年汶川大地震孕震、同震及震后变形和应力演化全过程的数值模拟

2008年M_S 8.0级汶川大地震发生在具有复杂的地质构造背景、强烈的地表起伏、不均匀的弹性和黏性结构的龙门山断裂带上。由于震前地震活动性不够强烈且地表构造变形较小,龙门山断裂带的地震危险性在汶川地震之前被低估。从数值模拟的角度,建立黏弹性有限元模型,考虑了初始地形、重力、构造加载、黏弹性松弛等因素对2008年汶川大地震的孕震、同震及震后150年变形全过程的影响,定量研究了映秀-北川断裂带的同震及震后变形,分析了弹性层、黏弹性层的应力积累、释放、调整的特点,模拟得到地表同震和震后位移与大地测量资料较为吻合,对汶川大地震的余震分布进行了力学上的解释,模拟得到震前、同震及震后的应力变化有助于深入分析大地震的动力学成因及其对周围区域的地震危险性影响。