Reflection and transmission of seismic waves at an interface betwen two saturated soils

ＲｅｆｌｅｃｔｉｏｎａｎｄｔｒａｎｓｍｉｓｉｏｎｏｆｓｅｉｓｍｉｃｗａｖｅｓａｔａｎｉｎｔｅｒｆａｃｅｂｅｔｗｅｎｔｗｏｓａｔｕｒａｔｅｄｓｏｉｌｓＪＵＮＹＡＮＧ（杨峻）ＳＨＩＭＩＮＧＷＵ（吴世明）ＺｈｅｊｉａｎｇＵｎｉｖｅｒｓｉｔｙ，Ｈａｎｇｚ...     

Seismic reflection and transmission coefficients at an air-water interface of saturated porous soil

Based on the modified Biot＇s theory of two-phase porous media, a study was presented on seismic reflection and transmission coefficients at an air-water interface of saturated porous soil media. The major differences between air-saturated soils and water-saturated soils were theoretically discussed, and the theoretical formulas of reflection and transmission coefficients at an air-water interface were derived. The characteristics of propagation and attenuation of elastic waves in air-saturated soils were given and the relations among the frequency, the angle of incidence and the reflection, transmission coefficients were analyzed by using numerical methods. Numerical results show that the propagation characteristic of the wave in air-saturated soils is great different from that in water-saturated soils. The frequency and the angle of incidence can have great influences on the reflection and transmission coefficients at interface. Some new cognition about the wave propagation is obtained and the study suggests that we may carefully pay attention to the influence of air on the dynamic analysis of seismic wave.     

Calculation of reflection and transmission coefficients for qP waves incident on a planar interface between isotropic and triclinic media

In the paper by Chattopadhyay and Rajneesh (2006, “Reflection and refraction of waves at the interface of an isotropic medium over a highly anisotropic medium’, Acta Geophysica, vol. 54, no. 3, pp. 239–249), the authors proposed a process to calculate R/T (reflection and transmission) coefficients at the interface between isotropic and triclinic half-spaces, with incident qP waves in triclinic media. Unfortunately, besides several misprints, the authors made a fatal assumption that there is no transmitted SH wave generated in isotropic media, which led the successive analytical derivations and numerical calculations thoroughly wrong. In this paper, the errors are analyzed at length and corrections are given. Then an alternative approach to solve the problem is proposed and numerical results are shown and discussed.     

Dynamic analysis of slab track on multi-layered transversely isotropic saturated soils subjected to train loads

The dynamic responses of a slab track on transversely isotropic saturated soils subjected to moving train loads are investigated by a semi-analytical approach. The track model is described as an upper Euler beam to simulate the rails and a lower Euler beam to model the slab. Rail pads between the rails and slab are represented by a continuous layer of springs and dashpots. A series of point loads are formulated to describe the moving train loads. The governing equations of track-ground systems are solved using the double Fourier transform, and the dynamic responses in the time domain are obtained by the inverse Fourier transform. The results show that a train load with high velocity will generate a larger response in transversely isotropic saturated soil than the lower velocity load, and special attention should be paid on the pore pressure in the vicinity of the ground surface. The anisotropic parameters of a surface soil layer will have greater influence on the displacement and excess pore water pressure than those of the subsoil layer. The traditional design method taking ground soil as homogeneous isotropic soil is unsafe for the case of RE 1 and RG 1, so a transversely isotropic foundation model is of great significance to the design for high train velocities.     

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.     

二维横各向同性弹性随机介质中的波场特征

本文通过交错网格有限差分正演．模拟了平面地震波在二维横各向同性弹性随机介质模型中的传播及其自激自收时间记录．为研究横各向同性弹性随机介质模型中的波场特征,我们在五个不同的时间区段上,分别计算剖面的三个统计特征(横向中心频率、纵向中心频率、波场能量相对值)．这样,对应每一个横各向同性弹性随机介质模型．均可计算得到15个不同的波场特征量．我们通过在二维横各向同性弹性随机介质中的正演模拟．研究当自相关长度以及介质的各向异性系数变化时,对应的上述波场特征量的变化特点．证实了在随机介质模型中．各向异性系数的变化会引起波场记录上的某些统计特征的变化,归纳得出了若干结论．     

横观各向同性饱和地基中埋置荷载的非轴对称瞬态响应

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

Surface wave propagation in a transversely isotropic elastic layer overlying a liquid saturated porous solid half-space and lying under the uniform layer of liquid

Dispersion of Rayleigh-type surface wave is studied in a homogeneous transversely isotropic elastic layer overlying a nondissipative liquid-saturated porous solid half-space and lying under a uniform layer of homogeneous liquid. The frequency equation in the form of ninth-order determinant is obtained.Special cases have been deduced by reducing the depth of the layers to zero and by changing the transverse isotropic layer to an isotropic layer. Dispersion curves for the phase velocity have been plotted for a particular model.     

横向各向同性地层(VTI)井孔声弹效应对弯曲波的影响

本文发展了建立在地层参考状态为各向同性介质假定下的现行井孔声弹性理论,就井外为横向各向同性面与井轴垂直的、具有9个独立三阶弹性模量的横向各向同性介质(VTI井况),水平面内受双轴应力作用下给出了一个简洁的与井内压力、应力差、应力和以及多极源偏振方位角有关的井孔弯曲波声弹公式,并且导出了平面纵、横波速度的声弹公式.数值考察了弯曲波速度之改变量的灵敏系数随频率的变化、受井外水平双轴应力作用时两种偏振的偶极弯曲波频散曲线以及对应不同方位径向偏振的横波速度.研究结果表明弯曲波声弹公式与5个二阶弹性系数以及7个独立的三阶弹性模量有关;而且由内压引起的井孔弯曲波声弹性公式中的三阶弹性模量仅与6个独立的弹性模量有关.横向各向同性介质井孔弯曲波速度的交叉现象仍是判断地应力存在的标志;一个重要的认识是受双轴应力作用的弯曲波速度变化在低频区主要与c₁₄₄和c₁₅₅两个三阶弹性模量有关,而且此认识与径向偏振的平面横波一致.在缺乏足够实验条件情况下,对VTI情况,以c₁₄₄,c₁₅₅和c₁₂₃三个独立的量进行测量,然后可暂不考虑三阶弹性模量的各向异性,建立简化的应力反演公式.反之,如果已知地层的地应力信息,由简化的声弹公式可以反演三阶弹性模量c₁₄₄,c₁₅₅和c₁₂₃.     

双相介质分界面上弹性波的反射与透射

本文基于Biot理论,推导出Zoeppritz形式的双相介质分界面上弹性波的反射与透射公式,对单相Zoeppritz公式与双相反射系数公式进行了比较,对双相介质含油、水、气不同流体时的反射规律以及孔隙度、渗透率、饱和度等储层参数对纵波反射的影响进行了研究.数值模拟分析表明,双相与单相反射公式的主要差异在于双相介质反射公式中考虑了液相、固液耦合相弹性模量的影响;油砂、水砂、气砂岩的慢P波反射差异明显;快P波反射对孔隙度的变化敏感,饱和度次之,对渗透率和频率的变化不敏感.     

含黏性流体孔隙介质分界面上反透射系数特征研究

地下岩石由岩石骨架和孔隙流体组成,通常流体含黏性.地震波在地下介质中传播时受岩石骨架和黏性流体的影响会呈现出复杂的变化.本文将流、固体位移和应力连续作为边界条件,推导出含黏性流体孔隙介质分界面上反透射系数方程;通过建立上层为饱油、下层为饱盐水的砂岩孔隙介质模型,开展反透射系数特征研究,分别分析不同频率、不同黏滞系数条件下,含黏性流体孔隙介质分界面上反透射系数随入射角的变化.研究表明,孔隙介质分界面上和等效介质分界面上的反透射系数分别随入射角的变化趋势基本一致,说明方程推导和数值计算的正确性;快纵波反透射系数受频率、流体黏性的影响较小,而快横波反透射系数在一定入射角范围内受频率、流体黏性的影响比较大;由于黏性孔隙流体的作用,慢纵波和慢横波的反透射系数受入射角、频率及流体黏性的影响都很大.

     

含两项不混合/单项黏性流体孔隙介质分界面上反透射系数特征研究

弹性孔隙介质分界面上的反透射系数特征,在岩性划分、流体识别、储层边界判识等方面有重要的应用.本文研究上层为含两项不混合黏性流体孔隙介质、下层为含单项黏性流体孔隙介质分界面上的反透射理论.首先根据两种孔隙介质分界面上的能量守恒得到边界条件,再将波函数、位移、应力与应变关系代入边界条件,推导出完全连通孔隙情况下,第一类纵波入射到孔隙介质分界面上的反透射系数方程.通过建立砂岩孔隙介质模型,分别分析不同孔隙流体类型、不同含油饱和度及不同入射角情况下,各类波的反透射系数特征.研究表明,第二、三类纵波反透射系数数值比第一类纵波小多个数量级,且两者对入射角的变化不敏感,但对孔隙流体性质、含油饱和度的变化较敏感,而横波反透射系数特征恰好与此相反;第一类纵波反透射系数特征比较复杂,入射角、孔隙流体的性质及含油饱和度的变化都对其产生影响.不同孔隙流体弹性物性的差异、孔隙介质中含油饱和度的变化及不同入射角引起垂向和切向应力分量的变化都会影响各类波的反透射系数特征,分析这些特征可以为研究储层含油气性提供理论基础.

     

Axisymmetric vibration of an elastic circular plate bonded on a transversely isotropic multilayered half-space

Based on the analytical layer-element method, an analytical solution is proposed to determine the dynamic interaction between the elastic circular plate and transversely isotropic multilayered half-space. The dynamic response of the elastic circular plate is governed by the classical thin-plate theory with the assumption that the contact surface between the plate and soil is frictionless. The total stiffness matrix of the transversely isotropic multilayered half-space is acquired by assembling the analytical layer-element of each soil layer with the aid of the continuity conditions between adjacent layers. According to the displacement condition of coordination between the plate and soil, the dynamic interaction problem is reduced to that of multilayered transversely isotropic half-space subjected to axisymmetric harmonic vertical loading. Some numerical examples are given to study the vertical vibration of the plate, and the results indicate that the dynamic response of elastic circular plate depends strongly on the material properties of the soils, the rigidity of the plate, the frequency of excitation and the external load form.     

流体-孔隙介质界面震电反射和透射系数的简化表达式

震电测井和海洋地震勘探中均遇到流体-孔隙介质界面.对于平面声波入射平面界面模型, 本文利用前人给出的震电反射和透射电磁波的电场式, 导出了震电反射和透射系数的表达式, 其中声波与电磁波的物理参数是解耦的.在给定介质参数和入射角且已知三类弹性体波的透射系数时, 就能通过该表达式计算震电反射和透射系数.该表达式对于所有入射角都是成立的, 当入射角为复数时对应非均匀平面波入射的情况.利用电磁波数远小于界面切向波数的前提条件, 将震电反射和透射系数的表达式化简为了简化表达式.简化式中与快纵波、慢纵波和横波有关的量均以透射系数与渗-固位移比的乘积形式出现, 而动电耦合系数、电导率和渗透率以乘积形式作为一个整体出现.简化式表明, 孔隙介质中的快纵波、慢纵波和横波通过渗流位移势影响震电反射和透射系数.数值算例表明, 当平面声波入射角不是极其小时, 震电反射和透射系数的简化式是适用的; 快纵波对震电反射和透射系数的影响可以忽略.震电反射和透射系数的简化表达式不仅简化了震电波场的计算, 而且明确表达了三类弹性体波引起反射和透射电磁信号的物理因素, 为地下和海底储层的参数反演提供了思路.

     

The reflection of plane waves in a poroelastic half-space saturated with inviscid fluid

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.     

Reflection and refraction of waves at the interface of an isotropic medium over a highly anisotropic medium

In this paper, we have considered the reflection and refraction of a plane wave at an interface between two half-spaces. The lower half-spaces is composed of highly anisotropic triclinic crystalline material and the upper half-space is homogeneous and isotropic. It has been assumed that due to incidence of a plane quasi-P (qP) wave, three types of waves, namely, quasi-P (qP), quasi-SV (qSV) and quasi-SH (qSH), will be generated in the lower half space whereas P and S waves will be generated in the upper half space. The phase velocities of all the quasi waves have been calculated. It has been assumed that the direction of particle motion is neither parallel nor perpendicular to the direction of propagation. Some specific relations have been established between directions of motion and propagation, respectively. The expressions for reflection coefficients of qP, qSV, qSH and refracted coefficients of P and SV waves are obtained. Results of reflection and refraction coefficients are presented.     

非均匀P\|偏振电磁波在导电界面的反射系数曲线

基于非均匀电磁波在导电媒质中传播时其相移常数和振幅衰减常数方向的不一致性，利用电磁波在导电媒质界面的边界条件，导出了非均匀P\|偏振电磁波的反射系数. 反射系数随入射角的变化曲线显示：当电磁波由电导率大的介质射向电导率小的介质时，在相移常数临界角和衰减常数临界角附近存在峰值，且反射系数的数值小于1，这一结果与全反射光的结果有明显的差异.     

横向各向同性地层随钻多极子测井声场的分波分析与特性研究

随钻声场在各向异性地层中的传播一直是随钻声波测井的重要课题.该问题中钻铤与井孔和地层的相互作用极大地影响着声场,特别是系统中的模式波的传播特征.因此,本文采用分波计算法来详细地研究在横向各向同性（VTI）地层中随钻多极子声波的传播规律,将声场全波分解成来自波数域中的奇异点（包括支点和极点）的分波的贡献,以此模拟分析了声场中各个分波的传播特征,包括模式波在频域中的频散、激发响应及时域中的分波波形.结果表明：在快速VTI地层中,钻铤的存在使得单极伪瑞利波和二阶四极地层螺旋波的频散大为降低,可以此来测量地层垂直向横波速度;慢速VTI地层中仅存在单一四极地层螺旋波,其低频速度为地层垂直向横波速度.值得注意的是,在各向异性参数满足δ>ε+c₄₄/2c₃₃的地层中,上述模式波的低频截止速度均小于地层的垂向横波速度,这一结果与已知的裸眼井情况一致.

     

The transmission of waves through a layered voigt solid sandwiched between isotropic homogeneous elastic media of same nature

Summary The paper is considered to show that the effect of solid viscosity on the transmission of waves through a visco-elastic layer sandwiched between isotropic homogeneous elastic media increases with decreasing wave lengths.     

Reflection of waves from the boundary of a random elastic semi-infinite medium

In this paper the smooth perturbation technique is employed to investigate the problem of reflection of waves incident on the plane boundary of a semi-infinite elastic medium with randomly varying inhomogeneities. Amplitude ratios have been obtained for various types of incident and reflected waves. It has been shown that an incidentSH orSV type of wave gives rise to reflectedSH, P andSV waves, the main components beingSH andP, SV in the respective cases. The reflected amplitudes have been calculated depending upon the randomness of the medium to the square of the small quantity , where measures the deviation of the medium from homogeneity. An incidentP-type wave produces mainly aP component and also a weakSH component to the order of ². The reflected amplitudes obtainable for elastic media are also altered by terms of the same order. The direction of the reflected wave is influenced by randomness in some cases.