Errors Due to the Common Ray Approximations of the Coupling Ray Theory

The common ray approximation considerably simplifies the numerical algorithm of the coupling ray theory for S waves, but may introduce errors in travel times due to the perturbation from the common reference ray. These travel-time errors can deteriorate the coupling-ray-theory solution at high frequencies. It is thus of principal importance for numerical applications to estimate the errors due to the common ray approximation.We derive the equations for estimating the travel-time errors due to the isotropic and anisotropic common ray approximations of the coupling ray theory. These equations represent the main result of the paper. The derivation is based on the general equations for the second-order perturbations of travel time. The accuracy of the anisotropic common ray approximation can be studied along the isotropic common rays, without tracing the anisotropic common rays.The derived equations are numerically tested in three 1-D models of differing degree of anisotropy. The first-order and second-order perturbation expansions of travel time from the isotropic common rays to anisotropic-ray-theory rays are compared with the anisotropic-ray-theory travel times. The errors due to the isotropic common ray approximation and due to the anisotropic common ray approximation are estimated. In the numerical example, the errors of the anisotropic common ray approximation are considerably smaller than the errors of the isotropic common ray approximation.The effect of the isotropic common ray approximation on the coupling-ray-theory synthetic seismograms is demonstrated graphically. For comparison, the effects of the quasi-isotropic projection of the Green tensor, of the quasi-isotropic approximation of the Christoffel matrix, and of the quasi-isotropic perturbation of travel times on the coupling-ray-theory synthetic seismograms are also shown. The projection of the travel-time errors on the relative errors of the time-harmonic Green tensor is briefly presented.     

地震波在垂向不均匀介质中传播问题的某些探讨

本文计算了含有高速夹层介质中首波的理论地震图。通过分析得到,当高速夹层薄到一定程度时,就会产生干涉型首波,从而从一个侧面证明了射线理论的局限性。通过对地震波反射—折射系数能量守恒关系的分析,探讨了反射—折射系数大于1的可能性。最后,介绍了一种计算垂向不均匀介质中拉梅问题理论地震图的数值方法——有限差分法。     

浅层有限频率面波成像中的3D灵敏度核分析

本文利用面波散射的模式耦合方法,基于波恩近似和远场假设,研究了有限频率面波三维灵敏度核,针对面波在工程应用中常遇到的水平分层的背景介质模型,计算了介质扰动引起的面波相位和幅度扰动的三维灵敏度核,分析了模式耦合对三维灵敏度核的影响.结果表明,仅考虑模式自身耦合的JWKB近似,介质密度和波速扰动引起的三维灵敏度核可以蜕化为...     

声波介质一次散射波场高斯束Born正演

Born正演是一种常用的地震波场正演模拟方法,也是线性化地震反演的理论基础.在实际应用时,Born正演通常结合常规的地震射线方法进行实现.为了克服常规地震射线方法的弊端,并且保证地震波场的模拟精度和计算效率,本文提出了一种基于高斯束的一阶散射波场Born正演方法.该方法分为两个环节:首先,我们利用高斯束的走时和振幅信息将地下散射点处的反射率映射为地表束中心位置处的局部平面波;然后,我们利用逆倾斜叠加将局部平面波转化为接收点处的时空域散射波场.在具体的实施过程中,我们提出一种以wavelet-bank方式实现的局部平面波合成方法,同现有的算法相比,可以在保持计算精度的同时,大大减少计算时间;此外,我们还利用最速下降法优化了高斯束的迭代循环过程,进一步提高了Born正演的计算效率.两个模型的应用效果证明,本文所提出的高斯束Born正演方法可以精确、高效的实现声波介质一次散射波场的正演模拟,为三维大规模地震波场的正演问题提供了一种切实可行的实现方案.     

Grid dispersion in generating finite-differences synthetic seismograms

Ｇｒｉｄｄｉｓｐｅｒｓｉｏｎｉｎｇｅｎｅｒａｔｉｎｇｆｉｎｉｔｅ－ｄｉｆｆｅｒｅｎｃｅｓｓｙｎｔｈｅｔｉｃｓｅｉｓｍｏｇｒａｍｓＡｂｄｏｌｒａｈｉｍＪａｖａｈｅｒｉａｎ（ＩｎｓｔｉｔｕｔｅｏｆＧｅｏｐｈｙｓｉｃｓ，ｔｈｅＵｎｉｖｅｒｓｉｔｙｏｆＴｅｈ...     

Application of the edge wave superposition method

Numerical examples of high-frequency synthetic seismograms of body waves in a 2-D layered medium with complex interfaces (faults, wedges, curvilinear, corrugated) are presented. The wave field modeling algorithm combines the possibilities of the ray method and the edge wave superposition method. This approach preserves all advantages of the ray method and eliminates restrictions related to diffraction by boundary edges and to caustic effects in singular regions. The method does not require two-point ray tracing (source-to-receiver), and the position of the source, as well as the type of source, and the position of receivers can be chosen arbitrarily. The memory and the time required for synthetic seismogram computation are similar to ray synthetic seismograms. The computation of the volume of the medium (the Fresnel volume or Fresnel zones), which gives the essential contribution to the wave field, is included in the modeling program package. In the case of complicated irregular interface (or a layered medium with a regular ray field at the last interface), the method displays a high accuracy of wave field computation. Otherwise, the method can be considered a modification of the ray method with regularization by the superposition of edge waves.     

二阶Born近似有限频率走时敏感核

有限频率层析成像考虑了非均匀介质中波的散射、衍射、波前愈合等物理性质,使得其对速度异常体的分辨能力远大于射线层析成像.推导和计算有限频率敏感核是进行有限频率层析成像的关键,当前推导有限频率敏感核多借助一阶Born近似,但这只适用于弱散射介质的情况.本文基于二阶Born近似并利用傅里叶变换推导了三维均匀介质情况下有限频率敏感核的解析表达式,并将其推广到非均匀介质中得到了三维非均匀介质中有限频率敏感核.研究表明:当介质中速度扰动小于2%时,基于二阶Born近似的有限频率敏感核与基于一阶Born近似的有限频率敏感核差别很小,可近似认为相同;当介质中速度扰动大于5%时,基于二阶Born近似的有限频率敏感核与基于一阶Born近似的有限频率敏感核有较大不同,表明此时已不能忽略二次散射.     

偶然干扰对地震脉冲震相对比的影响

在地震勘探中,有效地震波是在干扰背景上进行記录的,在記录上識別有效波一直是地震勘探的基本問題。本文中,我們假定地震脉冲是雷克对称形式的波漣,用随机过程分析法,討論了地震脉冲在平稳正态分布随机干扰影响下的幅度和相角分布函数及其它主要的統計特点,并指出地震脉冲波同相軸能够予以識别的条件,同相軸的可靠程度,能够予以識別的同相軸对应要求的最少地震脉冲波瞬时强度--門限值,并分析了用組合方法控制門限值的方法。     

Synthetic seismograms for finite sources in spherically symmetric Earth using normal-mode summation

Normal-mode summation is the most rapidly used method in calculating synthetic seismograms. However, normal-mode summation is mostly applied to point sources. For earthquakes triggered by faults extending for as long as several 100 km, the seismic waves are usually simulated by point source summation. In this paper, we attempt to follow a different route, i.e., directly calculate the excitation of each mode, and use normal-mode summation to obtain the seismogram. Furthermore, we assume the finite source to be a ‘‘line source' and numerically calculate the transverse component of synthetic seismograms for vertical strike-slip faults. Finally, we analyze the features in the Love waves excited by finite faults.     

TheoreticalPKP seismograms near the caustic calculated by sum of reflected waves

Summary This paper gives some examples of theoretical seismograms of PKP waves near the caustic. Seismograms of refracted waves for the original medium are compared with seismograms composed as a sum of the reflected waves, generated at boundaries of a substitute medium. All seismograms are calculated by zero approximation of the ray theory. The influence of some parameters of the source function and of the substitute medium on the results is shown.     

直达Ps与接收函数Ps联合分析地壳各向异性

本文尝试直接利用地震图记录中的Ps转换波,而非传统的接收函数方法,检测地壳介质的各向异性.虽然远震记录中的P波对Ps转换波会有干扰,但P波与Ps转换波具有不同的偏振特性,最小能量法可以将Ps转换波的各向异性参数检测出来,对时间窗的选择也很强壮.本文利用IC台网（New China Digital Seismograph Network,中国新数字地震台网）中的KMI,ENH,LSA,BJT,XAN,WMQ六个台站接收的宽频远震记录进行了测试,并且对比了接收函数方法的结果.研究表明,利用最小能量法直接测量与利用旋转相关法测量接收函数中Ps转换波的结果基本一致.两种方法测量所得的快波偏振方向与地壳运动方向,最大水平主压应力方向,震源机制解压缩轴方向,最大主压应变率方向以及已有的各向异性研究结果基本一致.     

弱各向异性介质弹性波的准各向同性近似正演模拟

准各向同性(QI)近似可用于弱各向异性介质的正演模拟.本文通过运用QI方法的零阶和一阶近似,计算了VTI介质模型的地震记录.得出的地震记录与标准各向同性射线理论(IRT)和基于伪谱法的三维地震正演模拟得出的地震记录作了比较,可以认为是精确的合成地震记录.     

WAVEFRONT DIVERGENCE,MULTIPLES, AND CONVERTED WAVES IN SYNTHETIC SEISMOGRAMS*

Synthetic seismograms can be very useful in aiding understanding of wave propagation through models of real media, verification of geologic models derived from interpretation of field seismic data, and understanding the nature and complexity of wave phenomena. If meaningful results are to be obtained from synthetic seismograms, the method of their computation must, in general, include three-dimensional geometrical spreading of wavefronts associated with highly concentrated (i.e., point) sources. The method should also adequately represent the seismic response of solid-layered media by including enough primaries, multiples, and converted phases to accurately approximate the total wavefield. In addition to these features, it is also very helpful, although not always essential, if the method of seismogram computation provides for explicit identification of wave type and ray path for each arrival. Various seismograms, computed via asymptotic ray theory and an automatic ray generation scheme, are presented for a highly simplified North Sea velocity structure. This is done to illustrate the importance of the above features and to demonstrate the inadequacy of the plane-wave synthesis method of seismogram computation for point sources and the limitations of acoustic models of solid-layered media.     

Theoretical ray seismograms of short periodPKP waves

Summary This paper shows some examples of theoretical seismograms of short period PKP waves calculated by the zero approximation of the ray theory. The influence of the epicentral distance, source time function and the Earth's crust on the form of the seismograms is shown.     

STABILIZATION OF NORMAL-INCIDENCE SEISMOGRAM INVERSION REMOVING THE NOISE-INDUCED BIAS*

A main problem in computing reflection coefficients from seismograms is the instability of the inversion procedure due to noise. This problem is attacked for two well-known inversion schemes for normal-incidence reflection seismograms. The crustal model consists of a stack of elastic, laterally homogeneous layers between two elastic half-spaces. The first method, which directly computes the reflection coefficients from the seismogram is called “Dynamic Deconvolution”. The second method, here called “Inversion Filtering”, is a two-stage procedure. The first stage is the construction of a causal filter by factorization of the spectral function via Levinson-recursion. Filtering the seismogram is the second stage. The filtered seismogram is a good approximation for the reflection coefficients sequence (unless the coefficients are too large). In the non-linear terms of dynamic deconvolution and Levinson-recursion the noise could play havoc with the computation. In order to stabilize the algorithms, the bias of these terms is estimated and removed. Additionally incorporated is a statistical test for the reflection coefficients in dynamic deconvolution and the partial correlation coefficients in Levinson-recursion, which are set to zero if they are not significantly different from noise. The result of stabilization is demonstrated on synthetic seismograms. For unit spike source pulse and white noise, dynamic deconvolution outperforms inversion filtering due to its exact nature and lesser computational burden. On the other hand, especially in the more realistic bandlimited case, inversion filtering has the great advantage that the second stage acts linearly on the seismogram, which allows the calculation of the effect of the inversion procedure on the wavelet shape and the noise spectrum.     

有限频率线性理论的波恩近似佯谬

对有限频率层析成像线性理论的波恩近似问题进行梳理, 用数值方法统计分析其适用范围, 结果表明波恩近似要求最大速度扰动不超过1%; 然后对相关走时一阶近似进行统计分析, 结果表明它也只适用于最大速度扰动在1%以内的情形. 然而, 结合波恩近似和相关走时一阶近似而得到的有限频率线性理论, 其适用的速度扰动范围最大可达10%. 这个表面上的逻辑悖论, 称为“波恩近似佯谬”. 此佯谬是由于不恰当地使用波恩近似造成的. 本文摒弃波恩近似, 使用泛函的Fréchet微分和隐函数定理推导得到有限频率线性理论, 圆满解释了波恩近似佯谬. 由于有限频率非线性理论早已摒弃了波恩近似, 因此波恩近似概念在有限频率层析成像理论中完全没有必要.      

Interpolation of the coupling-ray-theory Green function within ray cells

The coupling–ray–theory tensor Green function for electromagnetic waves or elastic S waves is frequency dependent, and is usually calculated for many frequencies. This frequency dependence represents no problem in calculating the Green function, but may pose a significant challenge in storing the Green function at the nodes of dense grids, typical for applications such as the Born approximation or non–linear source determination. Storing the Green function at the nodes of dense grids for too many frequencies may be impractical or even unrealistic. We have already proposed the approximation of the coupling–ray–theory tensor Green function, in the vicinity of a given prevailing frequency, by two coupling–ray–theory dyadic Green functions described by their coupling–ray–theory travel times and their coupling–ray–theory amplitudes. The above mentioned prevailing–frequency approximation of the coupling ray theory enables us to interpolate the coupling–ray–theory dyadic Green functions within ray cells, and to calculate them at the nodes of dense grids. For the interpolation within ray cells, we need to separate the pairs of prevailing–frequency coupling–ray–theory dyadic Green functions so that both the first Green function and the second Green function are continuous along rays and within ray cells. We describe the current progress in this field and outline the basic algorithms. The proposed method is equally applicable to both electromagnetic waves and elastic S waves. We demonstrate the preliminary numerical results using the coupling–ray–theory travel times of elastic S waves.     

A NONITERATIVE PROCEDURE FOR INVERTING PLANE-WAVE REFLECTION DATA AT SEVERAL ANGLES OF INCIDENCE USING THE RICCATI EQUATION*

Various exact methods of inverting the complete waveform of vertical seismic reflection data to produce acoustic impedance profiles have been suggested. These inverse methods generally remain valid for nonvertical, plane-wave data, provided total reflection does not occur. Thus, in principle, the “seismogram” at each ray parameter in a slant stack can be interpreted separately. Rather than invert each plane-wave seismogram separately, they can all be interpreted simultaneously and an “average” model thus obtained. Inversion for both the velocity and the density also becomes possible when two or more plane-wave seismograms are simultaneously inverted. The theory for a noniterative inversion method, based on the time-domain Riccati equation, is discussed. Numerical examples of inversions using this technique on synthetic data demonstrate its numerical stability and the advantage of simultaneous inversion of several seismograms to reduce the effect of noise in the data and increase the stability of the inversion process.     

Determining the focal mechanisms and depths of relatively small earthquakes using a few stations by full‐waveform modelling

Determining the focal mechanism of earthquakes helps us to better define faults and understand the stress regime. This technique can be helpful in the oil and gas industry where it can be applied to microseismic events. The objective of this paper is to find double couple focal mechanisms, excluding scalar seismic moments, and the depths of small earthquakes using data from relatively few local stations. This objective is met by generating three‐component synthetic seismograms to match the observed normalized velocity seismograms. We first calculate Green's functions given an initial estimate of the earthquake's hypocentre, the locations of the seismic recording stations and a 1D velocity model of the region for a series of depths. Then, we calculate the moment tensor for different combinations of strikes, dips and rakes for each depth. These moment tensors are combined with the Green's functions and then convolved with a source time function to produce synthetic seismograms. We use a grid search to find the synthetic seismogram with the largest objective function that best fits all three components of the observed velocity seismogram. These parameters define the focal mechanism solution of an earthquake. We tested the method using three earthquakes in Southern California with moment magnitudes of 5.0, 5.1 and 4.4 using the frequency range 0.1–2.0 Hz. The source mechanisms of the events were determined independently using data from a multitude of stations. Our results obtained, from as few as three stations, generally match those obtained by the Southern California Earthquake Data Center. The main advantage of this method is that we use relatively high‐frequency full‐waveforms, including those from short‐period instruments, which makes it possible to find the focal mechanism and depth of earthquakes using as few as three stations when the velocity structure is known.     

A central difference method with low numerical dispersion for solving the scalar wave equation

In this paper, we propose a nearly‐analytic central difference method, which is an improved version of the central difference method. The new method is fourth‐order accurate with respect to both space and time but uses only three grid points in spatial directions. The stability criteria and numerical dispersion for the new scheme are analysed in detail. We also apply the nearly‐analytic central difference method to 1D and 2D cases to compute synthetic seismograms. For comparison, the fourth‐order Lax‐Wendroff correction scheme and the fourth‐order staggered‐grid finite‐difference method are used to model acoustic wavefields. Numerical results indicate that the nearly‐analytic central difference method can be used to solve large‐scale problems because it effectively suppresses numerical dispersion caused by discretizing the scalar wave equation when too coarse grids are used. Meanwhile, numerical results show that the minimum sampling rate of the nearly‐analytic central difference method is about 2.5 points per minimal wavelength for eliminating numerical dispersion, resulting that the nearly‐analytic central difference method can save greatly both computational costs and storage space as contrasted to other high‐order finite‐difference methods such as the fourth‐order Lax‐Wendroff correction scheme and the fourth‐order staggered‐grid finite‐difference method.