首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 577 毫秒
1.
饱和多孔微极介质的波动方程及其势函数方程   总被引:1,自引:0,他引:1       下载免费PDF全文
胡亚元 《地球物理学报》2005,48(5):1132-1140
土是由一定尺寸大小颗粒所构成的多孔介质,具有明显的颗粒特性,当土颗粒间的孔隙被流体(如水或油)充满时则成为饱和土.利用微极理论和Biot波动理论的研究成果,把饱和土中多孔固体骨架部分近似地视为微极介质,孔隙中的流体部分视为质点介质,获得饱和多孔微极介质的弹性波动方程.借鉴Greetsma理论,建立了饱和多孔微极介质弹性本构方程力学参数与相应单相介质弹性参数的相互关系,使饱和多孔微极介质弹性波动方程中的物理参数具有明确的物理意义,易于在试验中确定.运用场论理论把饱和多孔微极介质的波动方程简化为势函数方程,建立了饱和多孔微极介质中五种弹性波的弥散方程,数值分析了五种简谐体波在无限饱和多孔微极介质中的传播特性. 结果表明,P1波、P2波和剪切S1波的波速弥散曲线与经典饱和多孔介质基本相同,当频率小于临界频率ω0时旋转纵波θ波和横波S2波不存在,当频率大于临界频率ω0时,θ波和S2波的传播速度随频率增加而减小.  相似文献   

2.
本文基于Biot的饱和多孔介质本构方程,考察具有辐射阻尼的外行球面波,推导了饱和多孔介质三维黏弹性人工边界的法向和切向边界方程;在已有的饱和多孔介质二维显式有限元数值计算方法基础上,提出该理论的三维方法,并开发了实现该三维方法的有限元程序.算例表明饱和多孔介质三维时域黏弹性人工边界与动力反应分析的显式有限元法具有较好的精度和稳定性.  相似文献   

3.
流体饱和多孔介质黏弹性动力人工边界   总被引:1,自引:0,他引:1  
基于Biot流体饱和多孔介质本构方程,采用平面波和远场散射波经验叠加来反映外行波传播,以经验参数反映人工边界外行波动的衰减和多角度透射特性。在人工边界处分别施加反映固相和液相介质传播效应的弹簧及阻尼来模拟人工边界以外的无限域介质对来自有限域的外行波的能量的吸收作用。从而形成一种流体饱和多孔介质的黏弹性动力人工边界。数值算例表明:边界的精度和稳定性高于现有的黏性边界、黏弹性人工边界及一阶透射边界。  相似文献   

4.
圆弧形凹陷饱和土场地对平面P波散射问题的解析解   总被引:13,自引:7,他引:13       下载免费PDF全文
以Biot 饱和多孔介质动力学理论为基础,利用Fourier Bessel级数展开法, 得到饱和多孔介质半空间中圆弧形凹陷地形对平面P波的散射问题的解析解. 数值计算给出地表位移幅值,分析了入射波波长、入射角、圆弧高宽比对地表位移幅值的影响,并与现存的在单相介质情况下得到的结论进行对比.  相似文献   

5.
用Biot饱和多孔介质动力学理论模拟半球形两相饱和土沉积谷,用单相介质弹性动力学理论模拟周围半空间场地.利用Fourier Bessel 级数展开法,在频域内给出了半空间中半球形饱和土沉积谷场地在平面Rayleigh波入射下三维散射问题的解析解.利用这一解析解计算分析了入射波频率、场地特征(包括孔隙比)对地表位移幅值的影响,并与已有的三维半球形沉积谷场地在单相介质中的散射问题进行了对比分析.  相似文献   

6.
具有饱和土沉积层的充水河谷对平面波的散射   总被引:6,自引:0,他引:6       下载免费PDF全文
把波函数展开方法用于地震波散射问题的研究中,首次在频域内给出了具有饱和土沉积层的圆弧形充水河谷对平面P波和SV散射问题的解析解答. 其中半空间场地用单相介质弹性动力理论模拟,河谷中的饱和土沉积层用饱和多孔介质的Biot动力学理论模拟,河谷中的水用无黏性流体(理想流体)介质模拟. 文中还给出算例,计算了不同高宽比的河谷谷底的位移幅值,分析表明河谷地形的存在使得半空间介质表面的位移幅值随着观察点位置的变化变化较大.  相似文献   

7.
三维半球形凹陷饱和土场地对平面P波散射   总被引:8,自引:5,他引:8       下载免费PDF全文
以Biot饱和多孔介质动力学理论为基础,利用Fourier Bessel 级数展开法,得到饱和多孔介质半空间中半球形凹陷地形对入射平面P波三维散射问题的解析解.利用这一解析解数值计算给出地表位移幅值,分析入射角、入射波频率对地表位移幅值的影响,并与现存的单相介质情况下三维半球形凹陷问题和两相介质情况下二维圆弧形凹陷问题得到的结论进行对比.  相似文献   

8.
王小岗 《地球物理学报》2009,52(8):2084-2092
基于孔隙介质的Biot理论,首先利用Laplace变换,给出圆柱坐标系下横观各向同性饱和弹性多孔介质在变换域上的波动方程;将波动方程解耦后,根据方位角的Fourier展开和径向Hankel变换,求解了Biot波动方程,得到以土骨架位移、孔隙水压力和土介质总应力分量的积分形式的一般解;借助一般解,建立了有限厚度饱和土层和饱和半空间的精确动力刚度矩阵,并由土层的层间界面连续条件建立三维非轴对称层状饱和地基的总刚度方程;在此基础上,系统研究了横观各向同性饱和半空间体在内部集中荷载激励下的动力响应,并给出了问题的瞬态解答.该研究为运用边界元法求解饱和地基动力响应奠定了理论基础.  相似文献   

9.
针对饱和多孔介质中热弹性波的传播特性问题,基于多孔介质理论和广义的热弹性模型,研究平面S波在饱和多孔热弹性介质边界上的反射问题。以考虑流-固耦合的饱和多孔介质波动方程和热-弹耦合的广义热弹性基本方程出发,建立饱和多孔介质的热-流-固耦合弹性波动模型。通过引入势函数并考虑自由透水和绝热的边界条件,经过理论推导最终给出在饱和多孔热弹性介质边界上的四种反射波的振幅反射率的理论表达式。在此基础上进行数值计算,分别讨论平面S波的入射频率、入射角和热膨胀系数等参数对四种反射波的振幅反射率的影响情况。结果表明:各反射波的振幅反射率分别随频率和热膨胀系数的增大而增大,同时也受到平面S波入射角变化的影响。该结论对于土动力学的理论研究及其相关的工程勘探具有一定的指导意义。  相似文献   

10.
饱和土中浅圆弧状埋置基础对平面P波散射解答   总被引:1,自引:0,他引:1  
李伟华 《华南地震》2008,28(1):11-20
在Biot饱和多孔介质动力学理论的基础上,利用Fourier—Bessel级数展开法,首次得到饱和土中浅圆弧状埋置基础对平面波散射问题的解答.并利用该解对有埋置基础的饱和土场地地面运动的幅值及其分布进行分析,研究了基础弹性对平面P波作用下饱和土地基与基础动力相互作用的影响。  相似文献   

11.
A model of wave propagation in fluid-saturated porous media is developed where the principal fluid/solid interaction mode affecting the propagation of the acoustic wave results from the conjunction of the Biot and the Squirt flow mechanism. The difference between the original Biot/Squirt (BISQ) flow theory and the new theory, which we call the reformulated BISQ, is that the average fluid pressure term appearing in the dynamic equation for a two component solid/fluid continuum is independent of squirt flow length. P-velocity and attenuation relate to measurable rock physical parameters: the Biot's poroelastic constants, porosity, permeability, pore fluid compressibility and viscosity. Modelling shows that velocity and attenuation dispersion obtained using the reformulated BISQ theory are of the same order of magnitude as those obtained using the original BISQ theory. Investigation on permeability effect on velocity and attenuation dispersion indicate that the transition zone in velocity and attenuation peak, occurring both at the relaxation frequency, shifts toward high frequency when permeability decreases. This behaviour agrees with Biot's theory prediction.  相似文献   

12.
This paper examines stresses and excess pore fluid pressure that are induced in a saturated poroelastic soil of halfspace extent by a concentrated line load. The line load is moving at a constant velocity along the surface of the poroelastic halfspace. The governing equations for the proposed analysis are based on the Biot's theory of dynamics in saturated poroelastic soils. The governing partial differential equations are solved using Fourier transforms. The solutions for the stresses and excess pore pressure are expressed in the forms of inverse Fourier transforms. The numerical results are obtained by performing the numerical inversion of the transform integrals. A parametric study is presented to illustrate the influences of the velocity of moving load and the poroelastic material parameters on the stresses and excess pore pressure. At a high velocity, the maximum values of the stresses in a poroelastic halfspace are smaller than those in an elastic solid, whilst at a low velocity the stresses in a poroelastic halfspace are larger than those in an elastic halfspace. The potential of diffusivity has an important influence on the stresses and excess pore pressure.  相似文献   

13.
This paper discusses surface displacements, surface strain, rocking, and energy partitioning during reflection-of-plane waves in a fluid-saturated poroelastic half-space. The medium is modeled by Biot's theory, and is assumed to be saturated with inviscid fluid. A linear porosity-modulus relation based on experimental data on sandstones is used to determine the material parameters for Biot's model. Numerical results in terms of angle of incident waves and Poisson's ratio are illustrated for various porosities and degrees of solid frame stiffness. The results show that the amount of solid frame stiffness controls the response of a fluid-saturated porous system. A poroelastic medium with essentially dry-frame stiffness behaves like an elastic medium, and the influence of pore fluid increases as dry-frame stiffness is reduced. The effects of a second P-wave become noticeable in poroelastic media with low dry-frame stiffness.  相似文献   

14.
Heterogeneous wave equations are more complicated numerically than homogeneous wave equations, but are necessary for physical validity. A wide variety of numerical solutions of seismic wave equations is available, but most produce strong numerical artefacts and local instabilities where model parameters change rapidly. Accuracy and stability of heterogeneous equations is achieved through staggered-grid formulations. A new pseudospectral staggered-grid algorithm is developed for the poroelastic (Biot) equations. The algorithm may be reduced to handle the elastic and acoustic limits of the Biot equations. Comparisons of results from poroelastic, elastic, acoustic and scalar computations for a 2D model show that porous medium parameters may affect amplitudes significantly. The use of homogeneous wave equations for modelling of a heterogeneous medium, or of a centred rather than a staggered grid, or of simplified (e.g. acoustic) wave equations when elastic or poroelastic media are synthesized, may produce erroneous or ambiguous interpretations.  相似文献   

15.
Tight oil/gas medium is a special porous medium, which plays a significant role in oil and gas exploration. This paper is devoted to the derivation of wave equations in such a media, which take a much simpler form compared to the general equations in the poroelasticity theory and can be employed for parameter inversion from seismic data. We start with the fluid and solid motion equations at a pore scale, and deduce the complete Biot's equations by applying the volume averaging technique.The underlying assumptions are carefully clarified. Moreover, time dependence of the permeability in tight oil/gas media is discussed based on available results from rock physical experiments. Leveraging the Kozeny-Carman equation, time dependence of the porosity is theoretically investigated. We derive the wave equations in tight oil/gas media based on the complete Biot's equations under some reasonable assumptions on the media. The derived wave equations have the similar form as the diffusiveviscous wave equations. A comparison of the two sets of wave equations reveals explicit relations between the coefficients in diffusive-viscous wave equations and the measurable parameters for the tight oil/gas media. The derived equations are validated by numerical results. Based on the derived equations, reflection and transmission properties for a single tight interlayer are investigated. The numerical results demonstrate that the reflection and transmission of the seismic waves are affected by the thickness and attenuation of the interlayer, which is of great significance for the exploration of oil and gas.  相似文献   

16.
The nonsplitting perfectly matched layer (NPML) absorbing boundary condition (ABC) was first provided by Wang and Tang (2003) for the finite-difference simulation of elastic wave propagation in solids. In this paper, the method is developed to extend the NPML to simulating elastic wave propagation in poroelastic media. Biot's equations are discretized and approximated to a staggered-grid by applying a fourth-order accurate central difference in space and a second-order accurate central difference in time. A cylindrical twolayer seismic model and a borehole model are chosen to validate the effectiveness of the NPML. The results show that the numerical solutions agree well with the solutions of the discrete wavenumber (DW) method.  相似文献   

17.
Numerical simulation in coupled elastic and poroelastic media is important in oil and gas exploration. However, the interface between elastic and poroelastic media is a challenge to handle. In order to deal with the coupled model, the first-order velocity–stress wave equations are used to unify the elastic and poroelastic wave equations. In addition, an arbitrary high-order discontinuous Galerkin method is used to simulate the wave propagation in coupled elastic–poroelastic media, which achieves same order accuracy in time and space domain simultaneously. The interfaces between the two media are explicitly tackled by the Godunov numerical flux. The proposed forms of numerical flux can be used efficiently and conveniently to simulate the wave propagation at the interfaces of the coupled model and handle the absorbing boundary conditions properly. Numerical results on coupled elastic–poroelastic media with straight and curved interfaces are compared with those from a software that is based on finite element method and the interfaces are handled by boundary conditions, demonstrating the feasibility of the proposed scheme in dealing with coupled elastic–poroelastic media. In addition, the proposed method is used to simulate a more complex coupled model. The numerical results show that the proposed method is feasible to simulate the wave propagation in such a media and is easy to implement.  相似文献   

18.
Three-dimensional time-harmonic response of a poroelastic half space subjected to an arbitrary buried loading is investigated. The analysis starts with the field equations in cylindrical coordinates based on Biot's general theory of poroelasticity. General solutions for the displacements are first derived using the Fourier expansions and Hankel integral transform with respect to the circumferential and radial coordinates, respectively. The transformed-domain solutions are obtained in explicit form. The physical-domain displacements and stress components are then obtained numerically by inverse integral transform. Comparisons illustrating the accuracy of the developed approach are made with existing solutions for an elastic half space, which is reduced directly from the general solution developed in the paper. Numerical results are presented for the displacements of a saturated soil subjected to a horizontal internal excitation.  相似文献   

19.
This review article consists of two parts. The first part concerns the admissibility of four contributions characteristic for Biot's model of poroelastic materials: coupling of partial stresses, presence of relative accelerations in equations of motion, dependence of permeability on frequency and changes of porosity. The second part is devoted to the demonstration of those contributions in the analysis of acoustic waves. Propagation of fronts, monochromatic waves and surface waves are considered.  相似文献   

20.
平面P波在饱和半空间中洞室周围的散射(I):解析解   总被引:1,自引:0,他引:1  
利用波函数展开法给出了入射平面P波在饱和半空间中圆形洞室周围散射问题的一个解析解。半空间假定为无粘性流体饱和介质,满足Biot理论。采用一种基于实验数据的孔隙率和模量之间的线性关系来确定Biot模型中的介质参数。解答考虑了透水边界和非透水边界两种情况。对边界条件进行了数值检验,结果表明,随着级数截断项数的增大,边界残量衰减很快。解答为进一步研究入射波频率和角度、边界渗透条件、孔隙率、泊松比等参数对散射的影响奠定了基础。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号