三维复杂介质中转换波走时快速计算

复杂介质中转换波走时计算是多波勘探地震学中重要内容之一.本项研究利用惠更斯原理和费玛原理,获得了三维复杂介质中转换波快速计算的改进型最小走时树方法.其中,在保证精度的条件下,为了提高三维转换波走时计算效率,首先对初至波最小走时树基本算法进行了改进.本方法通过将转换波分为上、下行波分别进行射线追踪以实现三维转换波走时的快速计算.模型计算表明,方法的计算速度快,而且稳定性强,对多波地震勘探具有较大的应用价值.     

VTI介质转换波动态聚焦束偏移成像方法研究

采用多分量地震数据进行处理,可以得到高质量的转换波成像结果,进而为流体识别和裂缝预测等提供有效支持.本文在前人研究的基础上,通过求解各向异性介质qPSV转换波运动学和动力学射线追踪方程,获得转换波中心射线走时、路径及动力学信息;同时,通过修改动态聚焦型传播算子,构建出各向异性介质动态聚焦束,提出了一种各向异性介质转换波...     

北京地震台远震Pn、PnPn波分析

在北京地震台NCDSN宽频带地震图集编辑工作中,发现Pn波在远至震中距16．8°的地震中仍然作为初至波被清晰记录,而PnPn波则延续至震中距30°左右仍以非常明显的震相特征被记录下来。     

山西地区远震P波到时残差分布

收集2008—2016年山西地震台网记录的震中距30°—90°范围内1 253个远震事件波形,拾取7 600余条高质量P波初至到时,使用IASP91模型计算相对到时残差,分析残差水平分布特征,结果显示:①以山西地区中部的山西断裂带为界,西部地震台站记录的P波初至主要表现为早到时,东部位于大同火山区的地震台站记录则主要表现为晚到时;②位于山西断裂带内部的地震台站记录的P波初至主要表现为早到时,残差水平显著低于西部地震台站;③研究区P波到时整体呈现自西向东逐渐由早到晚的分布特征。推测山西断裂带西部地区下方可能存在高速异常结构,山西断裂带内部及大同火山区下方可能存在低速异常结构。     

山西地区远震P波到时残差分布

收集2008-2016年山西地震台网记录的震中距30°-90°范围内1 253个远震事件波形,拾取7 600余条高质量P波初至到时,使用IASP91模型计算相对到时残差,分析残差水平分布特征,结果显示：①以山西地区中部的山西断裂带为界,西部地震台站记录的P波初至主要表现为早到时,东部位于大同火山区的地震台站记录则主要表现为晚到时;②位于山西断裂带内部的地震台站记录的P波初至主要表现为早到时,残差水平显著低于西部地震台站;③研究区P波到时整体呈现自西向东逐渐由早到晚的分布特征。推测山西断裂带西部地区下方可能存在高速异常结构,山西断裂带内部及大同火山区下方可能存在低速异常结构。     

转换波和非转换波与射线参数p的关系

本文对Aki固体-固体分界面的反射和透射公式进行适当变换,使之成为一组适合不同坐标系下,固体、液体、气体和真空之间所有可能形成的分界面上,入射P、SV和SH波所形成的反射和透射公式,从而把不同边界条件的反射和透射方程统一起来,并且阐明了转换波和非转换波动力学特征的差异.指出:在不同弹性介质分界面上入射平面弹性波,转换波的反射系数和透射系数是射线参数p的奇函数;非转换波的反射系数和透射系数是射线参数p的偶函数.根据奇函数和偶函数的性质可见:垂直入射时不存在转换波.同时也为简化Zoeppricz方程,开展AVO分析打开新思路.     

利用近震Sp转换波到时和远震接收函数P波峰值延迟研究华北盆地浅部结构

收集2017—2020年中国地震科学探测台阵在华北地区布设观测宽频带流动地震观测台采集的高信噪比近震远震地震波形,采用远震接收函数P波峰值延迟和近震Sp转换波到时对华北沉积层厚度进行了估算.由远震接收函数P波峰值延迟,近震Sp转换波和S波到时差方法推断的华北地区沉积层厚度在6～7 km,燕山造山带区域沉积层厚度较小,结晶基底埋深一般小于1 km,华北盆地中部和东南部盆地结晶基底埋深较厚,厚度普遍大于3 km,局部厚度可达7 km.盆地内地垒式隆起区邢衡隆起有3 km厚度的沉积层,冀中坳陷、黄骅坳陷厚度6～7 km;沉积层厚度与地质构造相对应,坳陷区沉积层厚度较大,隆起区沉积层厚度较小,体现了沉积盆地内部不同二级块体的沉降差异.震相拾取误差、震源方位角和震中距对沉积厚度计算结果影响不明显;进行时深转换采用的一维速度模型准确性对计算结果存在一定影响.总体来说,远震接收函数P波峰值延迟和近震Sp转换波到时可以研究沉积层厚度,能够可靠、有效地确定沉积层厚度结构的特征,为地震风险评估提供基础信息.     

转换波和非转换波与射线参数p的关系

本文对Aki固体-固体分界面的反射和透射公式进行适当变换,使之成为一组适合不同坐标系下,固体、液体、气体和真空之间所有可能形成的分界面上,入射P、SV和SH波所形成的反射和透射公式,从而把不同边界条件的反射和透射方程统一起来,并且阐明了转换波和非转换波动力学特征的差异.指出:在不同弹性介质分界面上入射平面弹性波,转换波的反射系数和透射系数是射线参数p的奇函数;非转换波的反射系数和透射系数是射线参数p的偶函数.根据奇函数和偶函数的性质可见:垂直入射时不存在转换波.同时也为简化Zoeppricz方程,开展AVO分析打开新思路.     

地震转换波对叠前同时反演的影响研究

叠前同时反演是油气探测的一种有效工具.其理论基础是平面P波Zoeppritz方程计算的反射系数的近似,是入射角的函数.叠前同时反演可以利用三项或两项Fatti方程进行反演分析.本文针对实际油田的测井数据,利用反射率法模拟了仅包含P波一次反射记录,包含P波一次反射和P波层间多次波记录以及全波场地震记录,再利用叠前同时反演对合成地震记录进行反演研究.研究结果表明,在大偏移距处P波主要反射受到其它模式波的污染,从而影响了叠前同时反演结果的精度.对于薄互层介质当转换波影响严重时,使用小角度数据的两项AVO反演比使用大角度数据的三项AVO反演更合理可靠.     

多层倾斜层状介质P-SV转换波转换点计算的最佳角度搜索法

转换波转换点的位置对转换波道集的抽取和叠加都非常重要.目前对于单一倾斜反射界面或水平层状介质模型,已经能较好地计算转换波转换点的位置.本文针对多层倾斜层状介质,提出最佳角度搜索法求取P-SV转换波转换点的位置,此方法是通过搜索最符合Snell定律的入射角和反射角来确定最佳的P-SV转换波转换点的位置.在搜索最佳转换点位置的同时,该方法也能给出最佳的P-SV波的传播路径.本文在单层和双层倾斜层状介质模型中测试了该方法,获得了精确的P-SV转换波转换点的位置和P-SV波的传播路径.     

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.     

Ray path of head waves with irregular interfaces

Head waves are usually considered to be the refracted waves propagating along flat interfaces with an underlying higher velocity. However, the path that the rays travel along in media with irregular interfaces is not clear. Here we study the problem by simulation using a new approach of the spectral-element method with some overlapped elements (SEMO) that can accurately evaluate waves traveling along an irregular interface. Consequently, the head waves are separated from interface waves by a time window. Thus, their energy and arrival time changes can be analyzed independently. These analyses demonstrate that, contrary to the case for head waves propagating along a flat interface, there are two mechanisms for head waves traveling along an irregular interface: a refraction mechanism and transmission mechanism. That is, the head waves may be refracted waves propagating along the interface or transmitted waves induced by the waves propagating in the higher-velocity media. Such knowledge will be helpful in constructing a more accurate inversion method, such as head wave travel-time tomography, and in obtaining a more accurate model of subsurface structure which is very important for understanding the formation mechanism of some special areas, such as the Tibetan Plateau.     

Stratigraphic filtering of seismic waves: The influence of mode conversion multiples onP-wave andS-wave phase and attenuation

We demonstrate how multiples, generated at the interfaces of plane parallel beds, modify the propagation characteristics of an originally coherent seismic wave. For waves propagating at an angle to the bedding plane we find that theSV andP-waves couple so that neither is a pure mode. TheSH-wave, while modified in its propagation characteristics by multiples, remains a pure mode. The coupling ofSV-multiples into the quasi-P-mode appears weaker than the coupling ofP-wave multiples into the quasi-SV mode; at least this is so for the two simple cases of (a) density fluctuations only and (b) correlatedV _p andV _s fluctuations which conserve Poisson's ratio.We also find that the coupling is sensitive to both the angle of propagation and frequency. In addition there is a cut-off angle forP-wave multiples influencing the quasi-SV mode. Propagation angles larger than the cut-off permit theP-multiples to modify the phase of the quasi-SV mode, but not its effective attenuation. No such cut-off effect is found for SV-multiples influencing the quasi-P mode, whose angle-dependent and frequency-dependent phase distortion and effective attenuation are influenced both byP-wave multiples andSV-multiples.In view of the mathematical complexity of the expressions describing the phase, and effective attenuation of modes when allowance is made forP-andS-wave multiples, we strongly advocate numerical coding of the major mathematical formulae. By so doing a systematic study can be undertaken of the frequency and offset dependence of seismic waves as a function of seismic source input and power spectral behavior of the fluctuations in density and elastic constants of beds. It is our opinion that the full mathematical expressions are too involved to permit an analytic, systematic investigation to be given of the phase and attenuation of seismic waves with any degree of sophistication or generality.     

横向不均匀介质中Kirchhoff-Helmholtz积分合成地震图

本文将Kirchhoff-Helmholtz积分方法应用于横向非均匀介质的合成地震图计算.与反射率方法及二维有限差分方法进行了比较.表明KH积分方法能很好的模拟反射波震相,且精度较高.KH 积分方法能够计算复杂界面的反射波,且其计算效率明显大于有限差分.     

Ray theory of modulated body waves in inhomogeneous viscoelastic media

Ray theory based upon real rays is presented for high frequencyP andS waves in continuously inhomogeneous isotropic linear media. Simple explicit formulae for the polarization anomalies described by higher order terms of the ray asymptotic expansions are considered.     

In-Plane Ray Tracing with Calculation of Out-of-Plane Geometrical Spreading in Anisotropic Media

The objective is to provide, in one single paper, a complete collection of equations governing kinematic and dynamic ray tracing related to a symmetry plane of an anisotropic medium. Well known systems for kinematic ray tracing and in-plane dynamic ray tracing are reformulated for the purpose of clarity, by taking advantage of a vector representation of the Christoffel matrix elements and related quantities. A generalized formula is derived for the integrand in out-of-plane dynamic ray tracing, pertaining to a monoclinic medium. Integrands corresponding to non-tilted orthorhombic and transversely isotropic media are obtained as special cases.     

Ray-centred coordinate systems in anisotropic media

Whereas the ray-centred coordinates for isotropic media by Popov and Pšenčík are uniquely defined by the selection of the basis vectors at one point along the ray, there is considerable freedom in selecting the ray-centred coordinates for anisotropic media. We describe the properties common to all ray-centred coordinate systems for anisotropic media and general conditions, which may be imposed on the basis vectors. We then discuss six different particular choices of ray-centred coordinates in an anisotropic medium. This overview may be useful in choosing the ray-centred coordinates best suited for a particular application. The equations are derived for a general homogeneous Hamiltonian of an arbitrary degree and are thus applicable both to the anisotropic-ray-theory rays and anisotropic common S-wave rays.     

Method for correction of the prestack migration amplitude of converted waves

The paper presents a method for correction of amplitude of prestack migration using the reflectivity function. The solution of wave equations for heterogeneous media expressed in the form of Neumann series for converted waves was used to find the reflectivity function. The performance of the proposed method was verified using synthetic models of wavefields. The synthetic models of wavefields were also used to compare amplitude correction levels required by the presented method and the standard amplitude correction methods, i.e., the methods using the spherical divergence and extention factor.     

基于成像射线的偏移剖面和速度时深转换评述

深度域地震波速度模型是地震勘探中最重要的参数之一,是获得精确深度域偏移剖面的基础.本文基于成像射线的概念,阐述了时间偏移域和深度域之间的关系,明确了时间偏移后的像点是通过成像射线与深度域的散射点联系起来的.在有横向变速的情况下,时间偏移剖面需要经过时深转换才能真实反映深度域层位,同样,从时间偏移剖面上提取的速度也需要经过时深转换才能得到深度域的速度.与前者不同的是,在速度的时深转换中,不仅速度的位置需要改变,且其数值也需要进行校正.该校正因子有几种等价的描述形式：度规量、速度扩散因子和几何扩散因子.