共查询到20条相似文献，搜索用时 78 毫秒
1.
孔隙介质弹性波传播理论在地球物理勘探、地震工程和岩土动力学等领域有着广泛的应用.而孔隙介质中的弹性波受孔隙度、渗透率、流体黏滞系数等参数的影响,因此研究波场的传播特征将有助于分析和提取这些信息.本文在Biot理论的基础上,针对三维层状孔隙介质模型,利用在合成理论地震图的研究中已经被证实具有稳定、高效且适用范围较广的Luco-Apsel-Chen（LAC）广义反透射方法,给出了弹性波场的一种积分形式的半解析解,可通过数值方法高效、准确地计算层状孔隙介质中的理论波场,所以该积分形式的半解析解可为三维层状孔隙介质波场传播特征的理论数值模拟研究提供一种新的途径和手段. 相似文献
2.
B.O. Ruud 《Studia Geophysica et Geodaetica》2006,50(3):479-498
The calculation of reflection and transmission coefficients of plane waves at a plane interface between two homogeneous anelastic
media may become ambiguous because it is not always obvious how to determine the sign of the vertical component of the slowness
vector of the scattered waves. For elastic media, the sign is determined by applying so-called radiation condition when the
slowness vector is complex-valued, but it has long been known that this approach does not work satisfactorily for anelastic
media. Other approaches have been suggested, e.g., by requiring that the reflection and transmission coefficients should vary
continuously with increasing incident angles, or by relating the sign to the direction of the energy flux. In the present
paper, it is shown that these approaches may give different results, and that the results can be inconsistent with the elastic
case even for weak attenuation. Instead, it is demonstrated that the ambiguity in the reflection coefficient can be resolved
by expressing the seismic response of a point source over an interface as a superposition of plane waves and their reflection
coefficients, and solving the resulting integral by the saddle point approximation. Although the saddle point itself (point
of stationary phase) does not provide new insight, the ambiguity is removed by considering the steepest descent path through
the point. Ray synthetic seismograms computed by this method compare well with synthetics computed by the reflectivity method,
which does not suffer from the above-mentioned ambiguity since the integration path is taken along the real axis.
This paper concentrates on the isotropic case, but it is discussed how the result may be extended to layered transversely
isotropic media. The suggested approach, derived for a point source and plane layers, does not directly apply to 2-D or 3-D
laterally inhomogeneous media, or to media of general anisotropy. A generalization of the result found is that the sign of
the vertical slowness components should be chosen according to the energy flux direction for subcritical incidence and according
to the radiation condition for supercritical incidence, even if this creates a discontinuity in the coefficients at the critical
incidence angle. Such a discontinuity is sometimes necessary to get results which are consistent with the elastic case. It
is discussed how the generalized result can be obtained by applying certain continuity criteria for the sub-and supercritical
angle intervals, but the validity of this approach for general models remains to be proved. 相似文献
3.
地震波在各向异性介质中以一个准P波(qP)和两个准S波(qS1和qS2)的形式传播.研究三种波的相速度、群速度以及偏振方向等传播性质能够为各向异性介质中的正反演问题提供有效支撑.具有比横向各向同性(TI)介质更一般对称性的正交各向异性介质通常需要9个独立参数对其进行描述,这使得对传播特征的计算更为复杂.当两个准S波速度相近时具有耦合性,从而令慢度的计算产生奇异性.因此,奇异点(慢度面的鞍点和交叉点)附近的反射与透射(R/T)系数的求解不稳定,会导致波场振幅不准确.本文首次通过结合耦合S波射线理论和基于迭代的各向异性相速度与偏振矢量的高阶近似解,得到了适用于正交各向异性介质以qP波入射所产生的二阶R/T系数的计算方法.与基于一阶近似的结果相比,基于二阶近似的方法提高了qP波R/T系数的精度,能得到一阶耦合近似无法表达的准确的qP-qS转换波的R/T系数解,且方法适用于较强的各向异性介质. 相似文献
4.
本文采用传播矩阵技术研究并建立了层状正交各向异性地层中多分量感应测井响应的有效算法.首先通过Fourier变换将频率空间域中的Maxwell方程组求解问题转化为频率波数域中关于电磁场水平分量常微分方程组的定解问题.利用该方程组系数矩阵的本征值和归一化本征向量将电磁场分解成上行波和下行波模式的组合,推导出均匀正交各向异性介质中由任意方向磁偶极子产生的电磁波模式解析表达式;在此基础上,利用叠加原理和边界条件研究了电磁波在层状正交各向异性地层中的反射和透射,给出各个界面上的广义反射系数和不同地层中电磁波振幅的递推公式,进而得到电磁波模式的解析解.为了有效确定频率空间域中的电磁场,采用二维Patterson自适应求积算法结合有限连分式展开技术计算傅氏逆变换.最后通过数值模拟结果证明了该算法的有效性,考察了不同各向异性系数、不同井眼倾角以及仪器长度和工作频率变化等情况下的多分量感应测井响应特征. 相似文献
5.
弹性孔隙介质分界面上的反透射系数特征,在岩性划分、流体识别、储层边界判识等方面有重要的应用.本文研究上层为含两项不混合黏性流体孔隙介质、下层为含单项黏性流体孔隙介质分界面上的反透射理论.首先根据两种孔隙介质分界面上的能量守恒得到边界条件,再将波函数、位移、应力与应变关系代入边界条件,推导出完全连通孔隙情况下,第一类纵波入射到孔隙介质分界面上的反透射系数方程.通过建立砂岩孔隙介质模型,分别分析不同孔隙流体类型、不同含油饱和度及不同入射角情况下,各类波的反透射系数特征.研究表明,第二、三类纵波反透射系数数值比第一类纵波小多个数量级,且两者对入射角的变化不敏感,但对孔隙流体性质、含油饱和度的变化较敏感,而横波反透射系数特征恰好与此相反;第一类纵波反透射系数特征比较复杂,入射角、孔隙流体的性质及含油饱和度的变化都对其产生影响.不同孔隙流体弹性物性的差异、孔隙介质中含油饱和度的变化及不同入射角引起垂向和切向应力分量的变化都会影响各类波的反透射系数特征,分析这些特征可以为研究储层含油气性提供理论基础. 相似文献
6.
黏弹各向异性介质中传播不均匀波,其反射、透射模式不仅与介质分界面两侧速度对比有关,还与品质因子Q的对比有关. 用伪谱技术模拟黏弹各向异性介质分界面上波的反射、透射,并与弹性各向异性介质、黏弹各向同性介质和弹性各向同性介质的模拟结果做比较. 计算平面波的反射、透射系数,分析介质的黏弹性和各向异性对反射、透射系数的影响. 数值模拟了一个三层介质模型中的波场,分析两个分界面上产生的反射波的特征. 黏弹各向异性介质中,qS波比qP波衰减程度大. 相似文献
7.
8.
A spectral-domain method, for the solution of the two-dimensional electromagnetic plane-wave scattering by a finite set of perfectly-conducting or dielectric cylinders buried in a dielectric half-space, has been developed. The scattered field is represented in terms of a superposition of cylindrical waves, and use is made of the plane-wave spectrum to take into account the reflection and transmission of such waves by the interface. The problem is solved for both the near- and the far-field regions, for TM and TE polarizations. In this work we briefly resume the theoretical basis of our approach. For configurations in which more obstacles are buried in the ground, and they are near to one another, we give details about the convergence rate of our method, and about the properties of our algorithms for the integration of cylindrical functions. With our technique it is possible to simulate two-dimensional buried obstacles of general shape, by means of a suitable set of circular-section cylinders: in this paper we show preliminary results of simulations carried out using arrays of same-radius circular cylinders, and of different-radius circular cylinders. 相似文献
9.
The bottom simulating reflector （BSR） in gas hydrate-bearing sediments is a physical interface which is composed of solid, gas, and liquid and is influenced by temperature and pressure. Deep sea floor sediment is a porous, unconsolidated, fluid saturated media. Therefore, the reflection and transmission coefficients computed by the Zoeppritz equation based on elastic media do not match reality. In this paper, a two-phase media model is applied to study the reflection and transmission at the bottom simulating reflector in order to find an accurate wave propagation energy distribution and the relationship between reflection and transmission and fluid saturation on the BSR. The numerical experiments show that the type I compressional （fast） and shear waves are not sensitive to frequency variation and the velocities change slowly over the whole frequency range. However, type II compressional （slow） waves are more sensitive to frequency variation and the velocities change over a large range. We find that reflection and transmission coefficients change with the amount of hydrate and free gas. Frequency, pore fluid saturation, and incident angle have different impacts on the reflection and transmission coefficients. We can use these characteristics to estimate gas hydrate saturation or detect lithological variations in the gas hydrate-bearing sediments. 相似文献
10.
单程波近似实际上是一种多次前向散射和单次后向散射近似.利用单程波近似来描述波传播可以极大地节省地震数值模拟的计算时间和内存,实现地震波长距离传播模拟和三维地震模拟快速计算.本文基于单程波近似和波动积分方程的分离变量逼近,从广义Lippmann-Schwinger波动积分方程推导出耦合反射/透射系数的单程波传播算子.该算子由两部分构成：分离变量Fourier单程波传播算子和薄板间的反射/透射系数表达.前者将常规的Fourier分裂步单程波传播算子(SSF)推广适应横向强速度变化介质和大角度传播波场.后者是利用垂直波数来表示反射/透射系数,自然耦合到波场传播的计算过程中,其为地质界面倾角的隐式表达,精确描述振幅随入射角的变化,能适应任意复杂的模型.通过两个数值算例和一个实际地质模型的计算,本文将该方法和边界元法进行了比较,结果表明：在算例给出的介质横向速度变化情况下,本文提出的方法在相位和振幅方面与全波数值方法基本吻合. 相似文献
11.
12.
R.-G. FERBER 《Geophysical Prospecting》1988,36(8):857-877
A new time-domain method is introduced for the calculation of theoretical seismograms which include frequency dependent effects like absorption. To incorporate these effects the reflection and transmission coefficients become convolutionary operators. The method is based on the communication theory approach and is applicable to non-normal incidence plane waves in flat layered elastic media. Wave propagation is simulated by tracking the wave amplitudes through a storage vector inside the computer memory representing a Goupillaud earth model discretized by equal vertical transit times. Arbitrary numbers of sources and receivers can be placed at arbitrary depth positions, while the computational effort is independent of that number. Therefore, the computation of a whole plane-wave vertical seismic profile is possible with no extra effort compared to the computation of the surface seismogram. The new method can be used as an aid to the interpretation of plane-wave decomposed reflection data where the whole synthetic vertical seismic profile readily gives the interpreter the correct depth position of reflection events. 相似文献
13.
岩石圈各向异性主要由上地幔矿物晶格优势排列方向和上地壳内裂缝、裂隙的定向分布造成.在各向异性特征显著区域,利用SKS震相剪切波分裂获得的延迟时间高达1.5s以上.本文根据方位各向异性,利用广义反射透射系数矩阵方法正演S波接收函数,研究各向异性对不同反方位角接收函数转换震相走时的影响.我们发展了基于HTI模型各向异性走时校正的方法,成功的在单层和多层(快轴方向相同或不同)的各向异性介质中对齐不同反方位角接收函数的Moho面和LAB的转换震相走时.我们将该方法应用于在青藏高原东北缘的流动台站,试图在各向异性强度较大区域对实测数据转换震相走时的校正效果进行测试.结果表明:各向异性走时校正能够加强单台接收函数转换震相的可追踪性,能量增强的叠后转换震相在时深转换后更利于对界面深度的识别与判断;在考虑三维成像的情况下,我们的各向异性校正方法对提高成像结果的准确性有重要意义. 相似文献
14.
15.
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. 相似文献
16.
进一步提高瞬变电磁法对地探测的解释精度,提出了回线源瞬变电磁成像原理及数值计算方法. 讨论了频率域中水平层状介质中瞬变电磁响应,得到一个以波阻抗为积分核的双重积分式;然后对水平层介质下电磁场的解进行上、下行波分离,得到含有以反射系数序列为未知的线性方程组,并给出了求取波阻抗和反射系数的数值解法：对实测磁场值进行域的变换,以均匀半空间下的等效波阻抗代替积分核函数,经过线性数字滤波后,在频率域求出等效波阻抗;把频率域中的波阻抗转换到时间域,以此为参数,构建方程组,在时间域用线性规划法求出反射系数序列. 最终以反射系数为参数进行成像. 对理论模型的数值计算结果表明,用本文提出的成像方法可以增强瞬变电磁法识别地下电性分界面的能力. 相似文献
17.
地下岩石由岩石骨架和孔隙流体组成,通常流体含黏性.地震波在地下介质中传播时受岩石骨架和黏性流体的影响会呈现出复杂的变化.本文将流、固体位移和应力连续作为边界条件,推导出含黏性流体孔隙介质分界面上反透射系数方程;通过建立上层为饱油、下层为饱盐水的砂岩孔隙介质模型,开展反透射系数特征研究,分别分析不同频率、不同黏滞系数条件下,含黏性流体孔隙介质分界面上反透射系数随入射角的变化.研究表明,孔隙介质分界面上和等效介质分界面上的反透射系数分别随入射角的变化趋势基本一致,说明方程推导和数值计算的正确性;快纵波反透射系数受频率、流体黏性的影响较小,而快横波反透射系数在一定入射角范围内受频率、流体黏性的影响比较大;由于黏性孔隙流体的作用,慢纵波和慢横波的反透射系数受入射角、频率及流体黏性的影响都很大. 相似文献
18.
In this article, we have derived a new and more general formulation of Love waves in arbitrarily irregular multi-layered media by using the global generalized reflection/transmission (abbreviated to R/T thereafter) matrices method developed earlier by Chen [17~20]. From the basic principle that the modal solutions are the non-trivial solutions of the free elastodynamic equation under appropriate boundary conditions, we naturally derived the characteristic frequencies and the corresponding distorted modes of Love wave in irregular multi-layered media. Moreover, we have derived the corresponding excitation formulation of Love waves in such laterally heterogeneous media by using the general solution of elastodynamic equation [17~20]. Similar to the result for laterally homogeneous layered structure, the Love waves radiated from a point source in irregular multi-layered media can be expressed as a superposition of distorted modes. Since the structure model used here is quite arbitrary, it can be used for so 相似文献
19.
自从发现震电现象以来,众多学者进行了相关研究.其中,Pride提出了一套描述流体饱和孔隙介质中震电波场的耦合与传播的宏观控制方程组,该方程组后来被广泛地应用到相关的震电研究中.Chen发展了一套广义反透射系数方法并将其应用到层状介质合成地震图的研究当中,该方法数值计算效率高并且可以处理带弯曲界面层状介质这种复杂模型.基于Pride的震电波场控制方程组,我们将Chen的广义反透射系数方法推广应用到层状孔隙介质中震电波场的数值模拟研究中,但是在数值计算过程中发现,当含源层的厚度相对于地震波波长较大时（即高频情况）,会出现数值计算的不稳定,此即为高频不稳定性问题.针对高频不稳定性问题,一种自然的处理方法就是在原来的含源层中插入两个虚拟界面,构造出一个新的含源薄层,但是这会带来一些额外的计算量,此外,由于虚拟含源薄层的厚度是有限的,必须针对具体模型参数设定一个合适的厚度值.高频不稳定性问题同样存在于层状介质合成地震图的数值计算过程中,Chen提出了一种解析的处理方法,即在原含源层内引入一无限薄的虚拟含源薄层,通过解析的方法解决高频不稳定性问题,该方法不会降低计算效率且适用于任意参数模型.本文首先对层状孔隙介质中的震电波场数值计算公式进行分析,指出源项积分中的指数增长因子是导致高频不稳定性问题的根本原因;其次将Chen在合成地震图数值模拟研究中采用的解析处理方法推广到震电波场研究中,得到了适用于数值计算的公式;然后给出数值算例,并针对一个含源层过厚的模型,比较了自然处理方法和解析处理方法,两种方法得到的结果具有相当好的一致性,而解析处理方法计算效率更高,证实了本文给出的解析处理方法在解决层状孔隙介质震电波场数值模拟的高频不稳定性问题方面的有效性. 相似文献
20.
Summary On the basis of Alsop's method (1966), approximated reflection and transmission coefficients are obtained for surface waves obliquely incident at a vertical discontinuity. Rayleigh and Love waves are arranged in such a way that they form a homogeneous eigenfunction system and their mutual conversion into each other is compatible with this theory. The structure of the eigenvalue equation is responsible for the mathematical complications in many surface wave problems. On the basis of a computer routine, first result concerning the angular dependence of the reflection and transmission coefficients including mode conversion are presented. The problem of phase jump is discussed. A comparison of theoretical and experimental results shows the applicability of Alsop's method for 90°-corners, too. 相似文献