利用伪谱法模拟横向各向同性介质中的波(英文)

SIMULATING SEISMIC WAVE PROPAGATION IN TRANSVERSELY ISOTROPIC MEDIA BY USING PSEUDO-SPECTRAL METHOD

介质的弹性常数为三维四阶张量的分量，共有８１个，由于应力张量和应变张量的对称性及能量密度是应变的二次函数，一般各向异常性介质的独立弹性常数可减为２１个，如果介质具有较高的对称性，独立弹性常数的数目会更少。 对于地壳和上地幔，具有５个独立弹性常数的横向各向同性介质是一个非常好的近似，本研究中横向各向同性介质的对称轴方向可以是任意的（即对称轴可以不平等于铅直方向），在此情况下，需要进行坐标变换，如果已知介质在某一坐标系（其坐标轴平行或垂直于介质的对称轴）中的弹性常数，我们能够容易地利用变换公式得到变换后新坐标系中的弹性常数。 本文提出了一种方案，利用伪谱法既能模拟横向各向同性介质中的平面波，也能模拟点源激发的波场。在勘探地球物理和地震学中，模拟横向各向同性介拮中传播的平面波及区域源产生的波是最重要的研究课题之一。然而在一般各向异性介质中，很难或不可能确定弹性波的相速度和偏振方向，但在横向各向同性介质中，则可以通过坐标变换来实现，这里我们所提出的方法可以用于横向各向同性介质中弹性波的模拟。

The elastic constant is a component of a three-dimensional fourth-rank tensor, hav- ing 81 components, in all. According to the symmetry of both the stress and strain ten- sors and the existence of a density energy function, which is quadratic in strain, the number of independent constants is 21 for general anisotropic media. The number of in- dependent elastic constants can be reduced still mord if media have higher symmetry. Transversely isotropic medium, which has only five indepindent constants, is a good ap- proximation of rocks in the crust and the upper mantle of the earth.In this paper, we are concerned about transversely isotropic media with an arbitiary direction of symmetrical axis(i. e., the symmetrical axis may not be parallel to the verticalaxis) . In this case we need to change coordinates from one system to another. If we know the elastic constants in one particular coordinate system, for example, whose axes are parallel or perpendicular to the symmetrical axis of the midia, we can easily obtain these elastic constants in new coordinate system by using the transformation formula. In this paper we present an approach for modeling wave-fields excited by not only a source but also a plane-wave incidence in transversely isotiopic media mentioned above by the pseudo-spectral method. Modeling of plane waves propagating in transversely isotropic media is one of the most important subjects as well as that of waves emitted from a localized source in exploration geophsics and seismology. While it is deffcult or even impossible to determine the phase velocity and the polarization direction of plane waves in general an anisotropic media, in the case of transversely isotropic media we can achieve this purpose through coordinate transformation. We here develop a scheme that can be used for plane-wave modeling in travsversely isotropic media,