WAVE FIELD EXTRAPOLATION TECHNIQUES FOR INHOMOGENEOUS MEDIA WHICH INCLUDE CRITICAL ANGLE EVENTS. PART I: METHODS USING THE ONE-WAY WAVE EQUATIONS*

Part I of this series starts with a brief review of the fundamental principles underlying wave field extrapolation. Next, the total wave field is split into downgoing and upgoing waves, described by a set of coupled one-way wave equations. In cases of limited propagation angles and weak inhomogeneities these one-way wave equations can be decoupled, describing primary waves only. For large propagation angles (up to and including 90°) an alternative choice of sub-division into downgoing and upgoing waves is presented. It is shown that this approach is well suited for modeling as well as migration and inversion schemes for seismic data which include critical angle events.     

利用矩阵分解理论的双程波方程叠前深度偏移方法

叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践.     

Numerical study of seismic scattering and waveguide excitation in faulted coal seams

Finite‐difference P‐SV simulations of seismic scattering characteristics of faulted coal‐seam models have been undertaken for near‐surface P‐ and S‐wave sources in an attempt to understand the efficiency of body‐wave to channel‐wave mode conversion and how it depends on the elastic parameters of the structure. The synthetic seismograms clearly show the groups of channel waves generated at the fault: one by the downgoing P‐wave and the other by the downgoing S‐wave. These modes travel horizontally in the seam at velocities less than the S‐wavespeed of the rock. A strong Airy phase is generated for the fundamental mode. The velocity contrast between the coal and the host rock is a more important parameter than the density contrast in controlling the amplitude of the channel waves. The optimal coupling from body‐wave energy to channel‐wave energy occurs at a velocity contrast of 1.5. Strong guided waves are produced by the incident S‐sources for source angles of 75° to 90° (close to the near‐side face of the fault). As the fault throw increases, the amplitude of the channel wave also increases. The presence of a lower‐velocity clay layer within the coal‐seam sequence affects the waveguiding characteristics. The displacement amplitude distribution is shifted more towards the lower‐wavespeed layer. The presence of a ‘washout’ zone or a brecciated zone surrounding the fault also results in greater forward scattering and channel‐wave capture by the coal seam.     

TTI介质的交错网格伪P波正演方法

研究了三维弱各向异性近似下,利用伪P波(伪纵波)模拟弹性波场P分量在倾斜对称轴的横向各向同性(TTI)介质中的传播过程,并对比了分别基于弹性Hooke定律、弹性波投影和运动学色散方程所建立的三种二阶差分伪P波方程的正演特点.目前这些伪P波方程数值计算主要采用规则网格差分,但是规则网格在TTI模拟中有低效率、低精度以及不稳定的缺点.为了提高计算的精度,本文构建出相应方程的交错网格有限差分格式.通过对比伪P波方程在三维TTI介质中不同的数值模拟的表达形式,本文认为基于色散方程所建立的伪P波方程在模拟弹性波中P波传播的过程中具有最小的噪声.本文分析不同的各向同性对称轴空间角度的频散特征,并引入适当的横波速度维持计算的稳定.二维模型算例表明,本文提出的交错网格正演算法可以得到稳定光滑的伪P波正演波场.使用本文交错网格算法对二维BP TTI模型的逆时偏移也具有较稳定的偏移结果.     

VSP上下行反射波联合成像方法研究

VSP资料上下行波场发育丰富.本文在分析VSP直达波、上行反射波、下行反射波传播路径及其照明范围的基础上,指出了常规VSP波动方程偏移方法缺陷,进而通过修改波场延拓方式,提出了上下行反射波联合成像方法,并在高频近似下分析了该方法的成像原理.该方法不需要进行VSP上下行反射波场分离,能够同时对VSP资料中的一次反射波、自由表面多次波、层间多次波进行成像,比常规成像剖面具有更宽的成像范围和更好的成像效果.该方法能够对下行一次反射波进行成像,从而可以实现常规偏移方法难以处理的高陡倾角构造成像.模拟资料和实际资料处理证明了本文方法的正确性.     

角度域叠前时间偏移:非均匀介质中的单程波方法

从单程波方程出发,推导出角度域叠前时间偏移的走时、入射波与反射波夹角、成像幅值计算方法;构建了可直接生成角度域成像道集的叠前时间偏移方法与偏移流程.文中定量分析了速度梯度对走时、角度、幅值的影响,给出了可更好考虑介质非均匀性的角度、幅值计算方法.理论模型以及实际数据验证了本文方法的有效性.     

基于波场传播方向的波动方程偏移有效去噪方法(英文)

本文研究分析了双程波波动方程偏移成像中广泛存在的三种主要噪声,特别是针对过去研究中没能很好解决的存在于高速盐丘悬垂边界附近的射线状噪声,提出了基于优化成像条件的有效去噪方法。射线状噪声主要来自于震源一侧波场的下行透射波分量和接收阵列一侧波场的上行散射波互相关成像。这一部分能量具有较强的互相关性,但并不携带真实的反射面信息。它广泛存在叠前偏移成像中,与信号在强度上同量级。多数情况下偏移成像中的相关噪声由方向性传播的波场能量产生。利用波场梯度得到的波场传播角度,可以分离出噪声对应的波场能量,并在成像条件中减去。采用这一方法可以有效地去除多种噪声,包括直达波噪声、散射波噪声和射线状噪声。该去噪方法不依赖波场外推算子,在需要时可以方便地运用到几乎所有的波动方程偏移中去。并且该去噪方法针对噪声的物理根源,对信号的损害很小。对去噪后的偏移成像结果额外地进行波数域滤波处理,可以进一步提高叠加图像的质量。这一去噪方法在超广角单程波偏移成像中取得良好效果,我们同时期待其在其他双程波波动方程偏移特别是逆时偏移(RTM)中的成功运用。     

The discrete wave number formulation of boundary integral equations and boundary element methods: A review with applications to the simulation of seismic wave propagation in complex geological structures

We review the application of the discrete wave number method to problems of scattering of seismic waves formulated in terms of boundary integral equation and boundary element methods. The approach is based on the representation of the diffracting surfaces and interfaces of the medium by surface distributions of sources or by boundary source elements, the radiation from which is equivalent to the scattered wave field produced by the diffracting boundaries. The Green's functions are evaluated by the discrete wave number method, and the boundary conditions yield a linear system of equations. The inversion of this system allows the calculation of the full wave field in the medium. We investigate the accuracy of the method and we present applications to the simulation of surface seismic surveys, to the diffraction of elastic waves by fractures, to regional crustal wave propagation and to topographic scattering.     

标量波动方程全倾角有限差分偏移

传统的有限差分法偏移是建立在单程波波动方程的不同阶数的近似式基础上的。因此,各阶近似式的偏移方法都存在一个倾角极限问题。克服倾角极限和提高极限倾角是八十年代以来有限差分偏移技术的研究目标。在这方面取得了显著的进展。本文从标量波动方程出发,通过函数替代,将它化为一个等价的方程组,用有限差分法解此方程组进行波场外推,实现全倾角偏移成象的目的。与国外已有的使用波动方程全式进行偏移的方法相比较,除方法原理本身不同外,该方法具有稳定性好,计算工作量较少的优点。     

High angle prestack depth migration with absorption compensation

The absorption effect of actual subsurface media can weaken wavefield energy, decrease the dominating frequency, and further lead to reduced resolution. In migration, some actions can be taken to compensate for the absorption effect and enhance the resolution. In this paper, we derive a one-way wave equation with an attenuation term based on the timespace domain high angle one-way wave equation. A complicated geological model is then designed and synthetic shot gathers are simulated with acoustic wave equations without and with an absorbing term. The derived one-way wave equation is applied to the migration of the synthetic gathers without and with attenuation compensation for the simulated shot gathers. Three migration profiles are obtained. The first and second profiles are from the shot gathers without and with attenuation using the migration method without compensation, the third one is from the shot gathers with attenuation using the migration method with compensation. The first and third profiles are almost the same, and the second profile is different from the others below the absorptive layers. The amplitudes of the interfaces below the absorptive layers are weak because of their absorption. This method is also applied to field data. It is concluded from the migration examples that the migration method discussed in this paper is feasible.     

VSP时空域高角度单程波方程偏移及其应用研究

偏移成像是VSP数据处理中的一个重要环节,常规的VSP成像方法通常利用VSP-CDP转换或Kirchhoff偏移,均存在保幅性差及成像精度低等问题,而波动方程叠前深度偏移被认为是对地下复杂构造进行成像的精确偏移方法.任意广角波动方程作为一种高精度的空间域单程波波动方程,同时由于只含有二阶偏导数项,易于数值实现,与其他单程波波动方程相比,具有更大的成像倾角,因此是偏移成像的有力工具之一.本文将AWWE推广应用到VSP数据成像中,实现了VSP时空域高角度单程波方程偏移.首先从三维标量任意广角波动方程出发,推导了完全匹配层吸收边界条件,在基本不增加计算量的前提下有效地压制了边界反射成像噪音,同时利用非线性反演算法优选参考速度来提高平方根算子的近似程度,从而提高高角度地层的成像精度.模型数值模拟实验验证了该方法的有效性,同时表明该方法在陡倾角构造情况下能取得很好的成像效果.最后对某地区实际观测的VSP资料进行了偏移成像,并与地面地震偏移结果进行了对比,显示出VSP波动方程偏移在成像分辨率上的优势.     

叠前逆时偏移影响因素分析

反射地震勘探中的偏移成像技术是获取地下介质构造形态最有效的手段之一.在叠前深度域偏移方法中,目前工业界采用的方法包括基于射线理论的波动方程积分解法和基于波动理论的微分波动方程单程波解法,这两类方法难以处理地震波横向速度变化剧烈的高陡倾角构造成像问题.近年来勘探地震学研究领域发展起来的叠前逆时偏移采用了双程波求解微分波动方程的算法,这种方法具有相位准确、不受介质横向速度变化和高陡倾角构造的影响、成像精度高、可以利用回转波正确成像等优点,从理论上弥补了当前工业界常规地震偏移所面临的成像缺陷.然而,叠前逆时偏移成像方法从理论走向实用尚需解决如下问题:计算速度和数据存储空间的节省、初始速度模型的建立、震源子波的选择、数值模型边界条件的定义和假像的消除等等.对于计算速度和存储量大的问题,随着计算机硬件的快速发展,将会不断得到改善,同时可以采取一些计算技术和存储策略来加以缓解.本文主要针对初始速度模型的建立、震源子波的选择、数值模型边界条件的定义和假像的消除这些因素,利用简单模型进行了分析.对于反射波造成的传播路径上的假像,给出了一种振幅补偿滤波方法.对勘探地球物理学界给出的SEG/EAGE二维盐丘模型、Marmousi模型和本研究设计的崎岖海底模型进行了叠前逆时偏移成像,均取得了较好的成像效果.     

基于波场垂向和水平方向外推的高陡构造偏移

常规的单程波波动方程偏移成像方法对大角度的高陡构造偏移成像存在内在的限制.根据波动方程在各个空间方向的数学特性和高陡构造反射地震波的传播特征,通过把地震波分解为垂向的上下行波、水平方向的前后行波和左右行波,提出基于波场垂向外推和水平方向外推相结合的单程波波动方程高陡构造偏移成像方法,即用波场垂向外推的单程波波动方程偏移成像方法解决中低角度平缓构造的偏移成像,用波场水平方向外推的单程波波动方程偏移成像方法解决中高角度陡倾构造的偏移成像.这种基于波场垂向和水平方向外推相结合的高陡构造偏移成像方法是常规单程波波动方程叠前深度偏移成像方法的补充和改进,它相对基于全波方程的逆时偏移具有计算效率上的优势.     

三维复杂构造中地震波模拟的单程波方法

复杂构造中单程波与双程波方法模拟结果的比较表明,就地震勘探中主要关心的一次反射波而言,单程波算法已具有足够的精度. 使用单程波方程将极大地减少数值计算的计算量,同时对介质的几何和物理参数建模也降低了要求. 单程波算法可视为深度偏移的“逆运算”,这样可以很好地借用已知的深度偏移方法及其程序系统. 基于计算效率和计算精度的双重考虑,本文在介质速度结构较复杂时采用显式短算子波场延拓方法,而在介质速度结构相对简单时采用分裂步相移法. 反射系数的计算中考虑了其随入射角的变化.     

伴随等离子体密度下降的磁声波与辐射带电子的波粒相互作用及其散射效应

利用范阿伦卫星的高质量观测数据,我们报道了伴随等离子体密度下降的磁声波现象.通过选取分别发生于2013年7月26日(事件A)和2013年9月19日(事件B)的两个相应事件进行细致分析,我们开展试验粒子模拟计算了磁声波对辐射带电子的散射系数,并求解二维福克-普朗克扩散方程量化了磁声波散射导致的辐射带电子动态变化.结果表明,事件A中的磁声波的散射作用主要发生于投掷角范围为60°～80°、能量范围为20～200keV的辐射带电子,而事件B中的磁声波的散射作用主要发生于投掷角范围为50°～80°、能量范围为20～400keV的辐射带电子;两个事件中的磁声波均能导致辐射带电子的蝴蝶状投掷角分布,但是由于事件B的磁声波幅度更强,形成的电子蝴蝶状分布更明显.     

单程波算子地震波入射角计算

基于单程波深度延拓方法,发展了一种地震波入射角度计算方法.入射角度的计算仅利用简谐波场,可得到整个成像区域内所有点的入射波波前面方向.该方法具有较高的计算效率,可服务于合成角道集等深度偏移方法;与偏移算法相比,其计算量几乎可以忽略.与射线法或基于走时梯度的入射角度计算方法相比,本文方法更稳健,避免了速度场的微小变化导致的入射角较大变化,因此更适用于实际偏移速度模型,也与波动方程深度偏移方法更匹配.数值算例表明,本文方法既有较高的计算效率又有很好的精度,且有很好的稳定性.     

基于两步SVD变换的零偏VSP资料上下行波场分离方法

垂直地震剖面(Vertical Seismic Profiling,VSP)资料处理中波场分离是关键问题之一.随着属性提取技术的发展,新的属性参数(例如Q值)提取技术对波场分离的保真性要求越来越高.本文改进了传统奇异值分解(Singular Value Decomposition,SVD)法,给出了一种对波场的动力学特征具有更好的保真性,可以作为Q值提取的预处理步骤的零偏VSP资料上下行波场分离方法.该方法通过两步奇异值分解变换实现:第一步,排齐下行波同相轴,利用SVD变换压制部分下行波能量;第二步,在剩余波场中排齐上行波同相轴,使用SVD变换提取上行波场.在该方法的实现过程中,压制部分下行波能量后的剩余波场中仍然存在较强的下行波干扰,使得上行波同相轴的排齐比较困难.本文给出了一种通过极大化多道数据线性相关程度(Maximize Coherence,MC)排齐同相轴的算法,在一定程度上解决了低信噪比下排齐同相轴的问题.将本文提出的方法用于合成数据和实际资料的处理,并与传统SVD法的处理结果进行对比,结果表明本文提出的波场分离方法具有良好的保真性,得到波场的质量明显优于传统SVD法.通过对本文方法和传统SVD法处理合成数据得到的下行波场提取Q值,然后进行对比可知,本文方法可以有效提高所提取Q值的准确性,适合作为Q值提取的预处理步骤.     

The One-Way Wave Equation: A Full-Waveform Tool for Modeling Seismic Body Wave Phenomena

The study of seismic body waves is an integral aspect in global, exploration and engineering scale seismology, where the forward modeling of waves is an essential component in seismic interpretation. Forward modeling represents the kernel of both migration and inversion algorithms as the Green’s function for wavefield propagation and is also an important diagnostic tool that provides insight into the physics of wave propagation and a means of testing hypotheses inferred from observational data. This paper introduces the one-way wave equation method for modeling seismic wave phenomena and specifically focuses on the so-called operator-root one-way wave equations. To provide some motivation for this approach, this review first summarizes the various approaches in deriving one-way approximations and subsequently discusses several alternative matrix narrow-angle and wide-angle formulations. To demonstrate the key strengths of the one-way approach, results from waveform simulation for global scale shear-wave splitting modeling, reservoir-scale frequency-dependent shear-wave splitting modeling and acoustic waveform modeling in random heterogeneous media are shown. These results highlight the main feature of the one-way wave equation approach in terms of its ability to model gradual vector (for the elastic case) and scalar (for the acoustic case) waveform evolution along the underlying wavefront. Although not strictly an exact solution, the one-way wave equation shows significant advantages (e.g., computational efficiency) for a range of transmitted wave three-dimensional global, exploration and engineering scale applications.     

Q-LOG DETERMINATION ON DOWNGOING WAVELETS AND TUBE WAVE ANALYSIS IN VERTICAL SEISMIC PROFILES1

Seismic attenuation introduces modifications in the wavelet shape in vertical seismic profiles. These modifications can be quantified by measuring particular signal attributes such as rise-time, period and shape index. Use of signal attributes leads to estimations of a seismic-attenuation log (Q-log). To obtain accurate signal attributes it is important to minimize noise influence and eliminate local interference between upgoing and downgoing waves at each probe location. When tube waves are present it is necessary to eliminate them before performing separation of upgoing and downgoing events. We used a trace-by-trace Wiener filter to minimize the influence of tube waves. The separation of upgoing and downgoing waves was then performed in the frequency domain using a trace-pair filter. We used three possible methods based on signal attribute measurements to obtain g-log from the extracted downgoing wavefield. The first one uses a minimum phasing filter and the arrival time of the first extremum. The two other methods determine the Q-factor from simple relations between the amplitudes of the first extrema and the pseudo-periods of the down-going wavelet. The relations determined between a signal attribute and traveltime over quality factor were then calibrated using field source signature and constant-Q models computed by Ganley's method. Q-logs thus obtained from real data are discussed and compared with geological information, specifically at reservoir level. Analysis of the tube wave arrivals at the level of the reservoir showed a tube wave attenuation that could not be explained by simple transmission effects. There was also a loss of signal coherence. This could be interpreted as tube wave diffusion in the porous reservoir, followed by dispersion. If this interpretation can be verified, tube wave analysis could lead to further characterization of porous permeable zones.