区域地震范围内地壳介质的低频响应

本文对于计算P_nl波格林函数的广义地震射线方法进行了改进。包括:1.简化了地震射线的描述方法;2.将射线按其运动学相似性进行分组计算,从而较大地提高了计算速度;3.进一步推广用于计算SH波部分的响应。对于1983年山东菏泽地震和1985年云南禄劝地震,用综合地震图方法模拟了763型长周期地震仪的P_nl波记录,求出两个地震的断层面解答,其结果与有关研究所得该地区构造应力场的图象基本一致。     

区域地震范围内地壳介质的低频响应

本文对于计算Pnl波格林函数的广义地震射线方法进行了改进。包括:1.简化了地震射线的描述方法;2.将射线按其运动学相似性进行分组计算,从而较大地提高了计算速度;3.进一步推广用于计算SH波部分的响应。对于1983年山东菏泽地震和1985年云南禄劝地震,用综合地震图方法模拟了763型长周期地震仪的Pnl波记录,求出两个地震的断层面解答,其结果与有关研究所得该地区构造应力场的图象基本一致。     

Space-time ray method for seismic wave fields

Summary The space-time ray method can be applied to the evaluation and continuation (extrapolation) of the complete seismic wave field in laterally inhomogeneous media with curved interfaces. The wave field propagates along certain space-time curves, called space-time rays. Their space projections correspond to standard rays. Examples of possible applications of the space-time ray method, where the standard ray method fails, are as follows: a) The propagation of seismic waves in slightly dissipative media, b) The computation of seismic wave fields generated by seismic sources with direction-dependent source-time variations. c) Downward continuation of the seismic wave field (actual seismograms) measured at the Earth's surface.     

多频段组合菲涅尔带走时层析成像方法研究

针对传统射线层析存在的种种局限性,菲涅尔带走时层析成像摒弃了传统的数学射线,考虑到地震信号具有一定的频带宽度,中央射线附近的介质对地震波的传播产生不同程度的影响。本文提出了多频段组合菲涅尔带走时层析成像方法。该方法以频率域波动方程Born和Rytov近似为基础,推导出建立在带限地震波理论基础上的波动方程Rytov近似走时敏感核函数,实现第一菲涅尔带约束下的波动方程走时层析反演方法。同时由于多个频段的引入,充分利用低频段和高频段的特有优势,从而兼顾菲涅尔带层析的计算效率与分辨率。模型试算结果证明了本方法的有效性和稳定性。     

Spatial dispersion of surface waves at the free surface of a general anisotropic elastic medium

A new technique relates the wave velocity of the surface waves in anisotropic elastic medium to its elastic constants. Anisotropic propagation of surface waves is studied in a half-space occupied by a general anisotropic elastic solid. The phase velocity expressions of quasi-waves, in three-dimensional space, are used to derive the secular equation of surface waves. The complex secular equation is resolved, analytically, into real and imaginary parts and is then solved, numerically, for phase velocity along a given phase direction on the surface. The complete procedure is thus analogous to the one used for conventional Rayleigh waves in isotropic medium. A non-linear equation relates the ray direction of the surface waves to its phase direction on the (plane) surface of the medium. The analytical differentiation of secular equation yields the directional derivative of phase velocity. This derivative is used to calculate the wave velocity of surface waves. Spatial variations of phase velocity, wave velocity and ray direction over the free plane surface are plotted for the numerical models of crustal rocks with orthorhombic, monoclinic and triclinic anisotropies.     

Recursive seismic ray modelling: applications in inversion and VSP

The recursive nature of rays in blocky models can be exploited to solve some difficult problems in seismic modelling. Each segment of a ray travels from an initial point up to a reflecting interface, where it is split into reflected and transmitted ray segments, which each continue in a similar way. The tree structure that thus emanates is conveniently handled by a recursive scheme. Recursion allows an automatic generation of all phases on a seismogram, together with all information necessary to analyse or select them. By operating recursively with a ray cell, bounded by a pair of vicinal rays in 2D, or a triplet of vicinal rays in 3D, and two successive isochrons, the two-point ray-tracing problem is reduced to a simple interpolation. Also, the cellular approach allows for a stable and robust evaluation of dynamic ray quantities without any paraxial tracing, which is cumbersome in blocky models of realistic complexity. Geometric shadows are filled by recursively generated diffractions. The recursive ray tracer has found applications in the fast computation of Green's functions in target-oriented inversion and in phase identification in VSP.     

Common-ray tracing and dynamic ray tracing for S waves in a smooth elastic anisotropic medium

Anisotropic common S-wave rays are traced using the averaged Hamiltonian of both S-wave polarizations. They represent very practical reference rays for calculating S waves by means of the coupling ray theory. They eliminate problems with anisotropic-ray-theory ray tracing through some S-wave slowness-surface singularities and also considerably simplify the numerical algorithm of the coupling ray theory for S waves. The equations required for anisotropic-common-ray tracing for S waves in a smooth elastic anisotropic medium, and for corresponding dynamic ray tracing in Cartesian or ray-centred coordinates, are presented. The equations, for the most part generally known, are summarized in a form which represents a complete algorithm suitable for coding and numerical applications.     

PP amplitude bias caused by interface scattering: are diffracted waves guilty?

This paper is concerned with the problem of interpretation of anomalous seismic amplitudes, induced by the amplitude‐scattering phenomenon. This phenomenon occurs in the vicinity of a crack distribution at the interface between elastic layers. The purpose of this work is to obtain a better understanding of the physics of this distinctive phenomenon, in order to interpret correctly the amplitudes of the reflected events. By analogy with studies in optics and in acoustics, we suggest that diffraction is widely involved in the amplitude‐scattering phenomenon. Analytical evaluation of the amount of energy carried by the reflected and the diffracted waves shows that neglecting diffraction in numerical models leads to local underestimation of the amplitude of waves reflected at interfaces with gas‐filled crack distribution.     

基于GSLS模型TI介质衰减拟声波方程

随着计算机硬件技术的发展以及高分辨率勘探需求的增加,我们希望能够更准确地模拟地下介质,得到更丰富的地层信息.然而,传统的声学假设并不能描述实际地层所存在各向异性和黏滞性,使得成像分辨率较低.为了实现深部储层的高精度成像,本文同时考虑了介质的各向异性和黏滞性,从TI介质弹性波的基本理论出发,结合各向异性GSLS理论,并通过声学近似方法导出基于GSLS模型的各向异性衰减拟声波方程.数值模拟表明该方程既能准确地描述各向异性介质下的准P波运动学规律,又能体现地层的吸收衰减效应;模型逆时偏移结果表明,在实现成像过程中考虑各向异性和黏滞性的影响,能对高陡构造清晰成像,且剖面振幅相对均衡,分辨率较高.     

固体介质中孤立波的传播及演化特征

由于地球介质中广泛存在断裂、微裂缝等地质现象,实际地震资料中会出现形似孤立波的非线性地震现象。因此,对固体介质中孤立波的研究有利于解释这些非线性地震现象的形成机制。本文基于KdV方程,以雷克子波作为初始条件,采用伸缩子机理构建体力模型,利用有限差分的方法模拟了孤立波的演化过程。理论结果表明,非线性地震纵波可以从雷克子波逐渐演化成孤立波,而且地震波的初始振幅和频散系数对模拟结果也有重要影响。通过与实际资料的对比也能说明这种演化的可能性。同时根据方程系数矩阵中元素带状分布的特征,采用稀疏矩阵的存取方法,可以减小计算内存,提高计算效率。     

2.5维地震波场褶积微分算子法数值模拟

早期的褶积微分算子都是基于正反傅立叶变换而实现的,其精度比四阶有限差分的精度稍高,本文将计算数学中的Forsyte广义正交多项式微分算子与褶积算子相结合,构建了一个新的快速、高精度褶积微分算子,其计算结果非常接近实验函数微分的精确值,精度与16阶有限差分的精度相当,远优于错格伪谱法的精确度.另外,2.5维数值模拟比二维模拟可以更真实地模拟三维介质的臬个剖面的波场,并且2.5维地震波模拟的计算量比三维模拟的计算量及计算耗时要大大减少.本文利用基于Forsyte广义正交多项式褶积微分算子法计算2.5维非均匀介质地震波场,模拟结果表明,该算法的计算速度快,计算精度高,能够直观、高效地反映复杂介质中波场的传播规律,并且2.5维波场数值模拟具有更高的计算效率,是一种非常值得深入研究并广泛应用的方法.     

Eigenrays in 3D heterogeneous anisotropic media,Part II: Dynamics

This paper is the second in a sequel of two papers and dedicated to the computation of paraxial rays and dynamic characteristics along the stationary rays obtained in the first paper. We start by formulating the linear, second‐order, Jacobi dynamic ray tracing equation. We then apply a similar finite‐element solver, as used for the kinematic ray tracing, to compute the dynamic characteristics between the source and any point along the ray. The dynamic characteristics in our study include the relative geometric spreading and the phase correction due to caustics (i.e. the amplitude and the phase of the asymptotic form of the Green's function for waves propagating in 3D heterogeneous general anisotropic elastic media). The basic solution of the Jacobi equation is a shift vector of a paraxial ray in the plane normal to the ray direction at each point along the central ray. A general paraxial ray is defined by a linear combination of up to four basic vector solutions, each corresponds to specific initial conditions related to the ray coordinates at the source. We define the four basic solutions with two pairs of initial condition sets: point–source and plane‐wave. For the proposed point–source ray coordinates and initial conditions, we derive the ray Jacobian and relate it to the relative geometric spreading for general anisotropy. Finally, we introduce a new dynamic parameter, similar to the endpoint complexity factor, presented in the first paper, used to define the measure of complexity of the propagated wave/ray phenomena. The new weighted propagation complexity accounts for the normalized relative geometric spreading not only at the receiver point, but along the whole stationary ray path. We propose a criterion based on this parameter as a qualifying factor associated with the given ray solution. To demonstrate the implementation of the proposed method, we use several isotropic and anisotropic benchmark models. For all the examples, we first compute the stationary ray paths, and then compute the geometric spreading and analyse these trajectories for possible caustics. Our primary aim is to emphasize the advantages, transparency and simplicity of the proposed approach.     

Efficient and accurate computation of seismic traveltimes and amplitudes

We describe two practicable approaches for an efficient computation of seismic traveltimes and amplitudes. The first approach is based on a combined finite‐difference solution of the eikonal equation and the transport equation (the ‘FD approach’). These equations are formulated as hyperbolic conservation laws; the eikonal equation is solved numerically by a third‐order ENO–Godunov scheme for the traveltimes whereas the transport equation is solved by a first‐order upwind scheme for the amplitudes. The schemes are implemented in 2D using polar coordinates. The results are first‐arrival traveltimes and the corresponding amplitudes. The second approach uses ray tracing (the ‘ray approach’) and employs a wavefront construction (WFC) method to calculate the traveltimes. Geometrical spreading factors are then computed from these traveltimes via the ray propagator without the need for dynamic ray tracing or numerical differentiation. With this procedure it is also possible to obtain multivalued traveltimes and the corresponding geometrical spreading factors. Both methods are compared using the Marmousi model. The results show that the FD eikonal traveltimes are highly accurate and perfectly match the WFC traveltimes. The resulting FD amplitudes are smooth and consistent with the geometrical spreading factors obtained from the ray approach. Hence, both approaches can be used for fast and reliable computation of seismic first‐arrival traveltimes and amplitudes in complex models. In addition, the capabilities of the ray approach for computing traveltimes and spreading factors of later arrivals are demonstrated with the help of the Shell benchmark model.     

Complex Rays and Wave Packets for Decaying Signals in Inhomogeneous,Anisotropic and Anelastic Media

Diffraction and anelasticity problems involving decaying, evanescent or inhomogeneous waves can be studied and modelled using the notion of complex rays. The wavefront or eikonal equation for such waves is in general complex and leads to rays in complex position-slowness space. Initial conditions must be specified in that domain: for example, even for a wave originating in a perfectly elastic region, the ray to a real receiver in a neighbouring anelastic region generally departs from a complex point on the initial-values surface. Complex ray theory is the formal extension of the usual Hamilton equations to complex domains. Liouville's phase-space-incompressibility theorem and Fermat's stationary-time principle are formally unchanged. However, an infinity of paths exists between two fixed points in complex space all of which give the same final slowness, travel time, amplitude, etc. This does not contradict the fact that for a given receiver position there is a unique point on the initial-values surface from which this infinite complex ray family emanates. In perfectly elastic media complex rays are associated with, for example, evanescent waves in the shadow of a caustic. More generally, caustics in anelastic media may lie just outside the real coordinate subspace and one must trace complex rays around the complex caustic in order to obtain accurate waveforms nearby or the turning waves at greater distances into the lit region. The complex extension of the Maslov method for computing such waveforms is described. It uses the complex extension of the Legendre transformation and the extra freedom of complex rays makes pseudocaustics avoidable. There is no need to introduce a Maslov/KMAH index to account for caustics in the geometrical ray approximation, the complex amplitude being generally continuous. Other singular ray problems, such as the strong coupling around acoustic axes in anisotropic media, may also be addressed using complex rays. Complex rays are insightful and practical for simple models (e.g. homogeneous layers). For more complicated numerical work, though, it would be desirable to confine attention to real position coordinates. Furthermore, anelasticity implies dispersion so that complex rays are generally frequency dependent. The concept of group velocity as the velocity of a spatial or temporal maximum of a narrow-band wave packet does lead to real ray/Hamilton equations. However, envelope-maximum tracking does not itself yield enough information to compute synthetic seismograms. For anelasticity which is weak in certain precise senses, one can set up a theory of real, dispersive wave-packet tracking suitable for synthetic seismogram calculations in linearly visco-elastic media. The seismologically-accepiable constant-Q rheology of Liu et al. (1976), for example, satisfies the requirements of this wave-packet theory, which is adapted from electromagnetics and presented as a reasonable physical and mathematical basis for ray modelling in inhomogeneous, anisotropic, anelastic media. Dispersion means that one may need to do more work than for elastic media. However, one can envisage perturbation analyses based on the ray theory presented here, as well as extensions like Maslov's which are based on the Hamiltonian properties.     

利用偏移进行视反射率估计的初步研究

视反射率估计是地震数据处理解释中的一项重要内容,通常采用反演的方法得到.本文以地震偏移和地震线性反演理论相结合为基础,并利用保幅单程波传播算子和保幅波动方程叠前偏移算法以及成像空间中的角度域波动方程偏移成像和照明补偿等方法技术,提出了一种利用单程波波动方程偏移进行地下反射面视反射率估计方法,并进行了理论模型的数值试验.这种估计方法得到的视反射率估计是一种近法向入射的小角度反射率.     

基于精细积分法的三维弹性波数值模拟(英文)

波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。     

地球介质自组织性对地震波走时和振幅的影响

本文采用随机介质描述地球内部大尺度背景场上存在的小尺度不均匀性和自组织结构;文中没有采用传统的层状网格结构,提出了块状结构来描述复杂的二维自组织结构的方法,分别以高斯型、指数型和von Karman型等自相关函数描述各向同性和各向异性非均匀分布自组织特征;采用射线追踪分析了不同分布特征自组织结构对地震波运动学和动力学特征参数的影响;结果表明,由于地球内部介质的自组织性存在,地震波射线轨迹可能发生明显的畸变; 不同偏移距处,地震反射振幅减弱或增强; 自组织结构从高斯型到von Karman型,在小尺度上表现更大的非均匀性,因此走时和振幅表征依次更强的平均效应.     

近似常Q模型的黏声各向异性纯qP波逆时偏移

地下介质中普遍存在各向异性及黏滞性, 各向异性会使地震波走时发生变化, 黏滞性会使地震波发生振幅衰减和相位畸变.在进行地震资料偏移成像时, 忽略介质各向异性及黏滞性的影响会导致绕射波不收敛、构造位置不准确, 能量不均衡, 偏移剖面频带变窄等问题, 大幅降低偏移剖面质量.本文基于TI介质纯qP波动方程, 结合常Q模型的松弛函数, 推导出新的黏声各向异性纯qP波动方程.在该方程中, 振幅衰减项与相位频散项是解耦的, 因此利用该方程进行成像时, 可通过改变振幅衰减项的符号进行衰减补偿.多个数值算例结果表明文中所提出的方程可以正确描述地震波在各向异性衰减介质中的传播特征, 应用本文方程进行成像时能够有效校正各向异性及黏滞性对偏移剖面的影响, 提高成像剖面的分辨率及精度.

     

复杂地表边界元-体积元波动方程数值模拟

复杂近地表引起来自深部构造的地震反射信号振幅和相位的异常变化,是影响复杂近地表地区地震资料品质的主要原因.本文采用边界元-体积元方法,通过求解含复杂地表的波动积分方程,来模拟地震波在复杂近地表构造中的传播.其中,边界元法模拟地形起伏和表层地质结构对地震波传播的影响;体积元法模拟起伏地表下非均质低降速层的影响.与其他数值模拟方法比较,其主要优点为几何上精确描述不规则地表界面,实现精确模拟自由表面对地震波的边界散射;显式应用近地表地层界面的连续边界条件,实现半解析的数值模拟;分区处理近地表复杂结构,有效模拟复杂地表下非均匀介质对地震波场的体散射.数值试验结果表明了该方法的实用性和有效性.     

倾斜界面中反射波走时计算

对于平层界面和倾斜界面的反射地震波,基于射线理论,采用波动方程进行数值模拟;使用Matlab软件,编写2种界面不同倾角、不同震中距时反射地震波走时程序,绘制走时曲线,分析地震波的传播特点,为研究地震波在地球内部的传播路径提供参考。