多波多分量高斯束叠前深度偏移

本文对基于弹性波动理论的多波多分量高斯束偏移进行了完整且详细的分析和公式推导,实现了3D空间多分量(矢量)波场的直接成像.由于当前多数基于弹性波动方程的偏移方法只是假设应力边界条件为自由地表边界条件,这种假设不符合垂直地震剖面(VSP)和海底电缆(OBC)等地震数据.为此本文详细分析了实际应用中常见的三种弹性各向同性介质模型的应力边界条件:自由空间、海底和自由地表模型.在上行传播假设情况下,获得了应力边界条件与位移边界条件的关系式.在此基础上,准确推导了3D多波多分量高斯束波场延拓和偏移成像公式,并在偏移过程中实现了完整的多波型自动分离.由于常规的互相关成像条件不适用于矢量波场成像,本文引用了散度/旋度互相关成像条件.通过约定PS转换波的正向旋转方向解决了3D空间PS成像极性翻转问题.利用2D和3D模型数据偏移成像验证了我们所提出的多波多分量高斯束偏移方法的可行性.     

基于海底电缆数据的声弹耦合介质多参数最小二乘逆时偏移（英文）

在海洋地震勘探中,海底电缆技术采集的多分量地震数据,涉及到多参数的反演与成像问题,本文针对海底电缆多分量地震数据提出了一种弹性波多参数最小二乘逆时偏移方法。该方法的波场延拓算子为混合方程,即在海水介质中采用声波方程进行波场计算,而海底固体介质的波场由纵横波形分离的矢量弹性波方程得到。在海底界面中采用声弹耦合控制方程将两种类型的方程结合起来。通过推导纵横波形分离的弹性波偏移算子、反偏移算子和梯度公式实现基于纵横波形分离的弹性波最小二乘逆时偏移方法。通过模型试算证明了该方法能够得到高质量的纵波速度和横波速度分量的成像剖面,相比与传统弹性波最小二乘逆时偏移方法,多参数耦合造成的成像串扰噪音得到了很好地压制。     

一种多波联合预测油气技术研究及应用

常规纵波(P波)油气检测利用单一纵波波场信息进行预测,预测结果多解性较强,而本文实现了多种波场信息综合预测,可减少预测的多解性,提高预测的准确性.多波多分量地震勘探比常规纵波地震勘探可获得到更多的波场信息,包括纵波(P波)波场信息和转换横波(P-SV波、P-SH波)波场信息,这为多波场信息综合运用减少勘探多解性提供了较好的资料基础.首先,对P波进行地震层位解释,可得到P波目的层顶、底时间;其次,对P-SV波进行地震层位解释,可得到P-SV波目的层顶、底时间;再次,对P波做P波AVO分析,可求得P波AVO截距;再次,对P-SV波做P-SV波AVO分析,可求得P-SV波AVO截距;最后,计算机编程研制多波流体因子计算程序,同时输入上述6个参数,计算得到多波多分量流体因子数据体,据此可进行多波多分量油气预测.根据实际工区应用,取得了较好的实际应用效果.实际应用结果表明:本方法技术是可行的、可靠的,可提高油气预测的精度和可靠性.     

低信噪比海底多分量地震数据波场分离方法

以多分量地震观测为基础,联合纵波和转换横波数据能更有效地估计地下介质的弹性和物性参数,提升地质构造成像与油气储层描述的精度.在海底多分量地震数据处理过程中,观测记录的上-下行波分解和P/S波分离可压制水层鸣震以及P与S波之间的串扰,对偏移成像和纵横波速度建模至关重要.但受海底环境、仪器与观测因素共同影响,许多海底多分量地震资料都无法基于现有的海底波场分离方法与流程取得合理的结果.本文以海底声波场与弹性波场分离基本原理为基础,通过对方法流程的修正,摆脱常规流程对中小偏移距直达波信号的依赖性.借助模拟数据实验讨论了波场分离对海底介质参数、噪声的敏感性.结合东海YQ探区海底多分量地震资料上-下行P/S波分离及其叠前深度偏移处理,验证了本文方法流程的可行性.     

多分量地震资料处理解释技术研究

随着油田勘探开发的不断深入,多波多分量地震资料的应用也越来越广泛.本文结合实际资料,对多波多分量地震资料的处理和解释方法进行了探讨和研究,主要针对多波资料的静校正、叠加成像、转换波资料的共转换点的求取方法等技术进行了分析研究,开发了综合利用多波资料进行地震属性和各向异性研究的方法,形成了一套较完善的多分量地震资料的处理、解释技术系列.经过在胜利油田罗家地区的初步应用,表明了预测结果与钻井资料相吻合.     

应用三分量浅层地震反射方法探测隐伏活动断裂

多波多分量地震勘探方法在油气田勘探开发中日益受到重视,但在工程和活动断裂探测中过去还没有开展多波多分量地震勘探的先例.本文简要介绍了在北京黄庄—高丽营断裂上开展的三分量浅层地震反射试验数据采集和资料处理方法,以及取得的初步结果.试验结果表明,采用纵波震源或横波震源激发、三分量接收的浅层地震勘探方法,可获得较高信噪比的浅层地震纵波、横波和转换波剖面,特别是水平分量采集的横波和转换波可提供很有意义的地下结构和构造信息.综合利用纵波、横波和转换波剖面特征进行的地质分层和断层解释结果取得了与跨断层的钻孔联合地质剖面结果较好的一致性.     

一种多波宽频地震数据采集设计与仪器实现关键技术

随着多波多分量技术的发展,宽频数据采集逐渐成为地震勘探的发展趋势.通过高保真获取全波场信息,并尽可能保留高频有效成分,进而利用多种波场信息相互验证和补充,是提高地震勘探精度的有效方法.受仪器带宽限制,常规设备在高频段很难同时保持高分辨率和大动态范围特征,在多波地震采集时,极容易出现分辨率不足和限幅失真等现象,极大影响多波多分量方法的实际应用效果.鉴于此,本文给出一种利用积分器、高频ADC、DAC和抽取滤波器构成的低成本宽频数据采集方法.通过构造环路滤波器(调制器),对高频ADC输出数据中的量化噪声二次整形,并配合抽取滤波,对带外噪声进行有效压制.实际测试表明,通过精细电路结构设计,采集站有效分辨率可达24位(采样率为2 kHz),动态范围超过140 dB,SNR达142 dB.相比传统设备,其采集数据在5～25 Hz的低频段和100 Hz以上的高频段都具有更加丰富的有效信息,非常适合多波多分量数据采集.     

二维弹性多波时空域高斯束偏移方法

随着我国勘探开发难度逐步增大,勘探目标开始向裂缝油气藏、岩性油气藏等复杂探区转移,研究高精度、适应性强的多波多分量深度偏移算法在后续的地震解释、属性分析及储层预测中具有重要意义.针对多波多分量地震数据,本文提出了一种二维弹性波时空域高斯束偏移方法.时空域高斯束沿中心射线传播时能够面向成像目标描述局部波场,且对振幅和频率可调制的Gabor基函数有天然的适应性,因而将基于Gabor分解的子波重构方法应用于震源波场构建,从而得到任意点源函数产生的时空域高斯束波场.该方法由于直接在时间域进行计算,可以避开频率域中出现的假频和边缘截断效应等问题.基于各向同性弹性波动方程的Kirchhoff-Helmholtz积分解,利用矢量时空域高斯束传播算子构建格林函数和格林位移张量,并结合上行射线追踪策略,实现了检波点波场的反向延拓.针对矢量波成像问题,本文借鉴弹性波逆时偏移方法从矢量延拓波场中分离出纯纵波分量和纯横波分量,进而采用修改后的内积成像条件产生具有明确物理意义的PP、PS成像结果,避免了转换波成像的极性反转问题.最后利用简单两层模型和不含盐体构造的部分Sigsbee2a模型的成像结果,并将其与应用近似纵横波成像条件、标量和矢量势成像条件的偏移剖面进行对比,验证了本文方法的正确性和有效性.     

多波多分量地震探测技术

随着油气勘探程度的提高,越来越多地显示出地下介质的复杂性,多波多分量地震勘探可以提供更多的信息。本文主要论述多波多分量地震勘探在各向异性介质中应用时的效应,可以解决的实际问题和各向异性介质中地震波传播特征,并给出实例说明多波多分量地震勘探的应用范围。     

带有多道相关的抛物线Radon变换法分离P-P、P-SV波

因为多波多分量地震勘探中P-P波和P-SV波通常混杂在一起,所以较好地分离P-P波和P-SV波能够提高数据处理和解释的质量.抛物线Radon变换法在分离P-P波和P-SV波时取得了一定效果,但是在离散叠加的计算过程中会带来假频,这些假频会干扰波场分离.本文针对这一问题,将多道相关算法引入抛物线Radon变换,发展了带有多道相关的抛物线Radon变换法.该方法利用叠加信号具有相似性的特点,依据多道相关中衡量多道信号相似性的能量比标准,对叠加过程加以控制,压制变换中出现的假频.本文用该方法对合成地震记录进行了波场分离,取得了较好的效果.     

VTI介质P波非双曲时差分析

具有垂直对称轴的横向各向同性介质模型(VTI)是目前各向异性理论研究和多波多分量地震资料叠前成像处理中最常用的一种各向异性模型.VTI介质中反射 P波时距曲线一般不再是双曲线.基于不同的相速度近似公式会得到不同的时距关系式.文中对几种典型的非双曲时距曲线与射线追踪得到的准确时距曲线在不同各向异性强度下进行了对比，结果表明Muir等和Stovas等提出的非双曲时距公式由于过高地考虑了横波垂直速度的影响与精确的时距曲线有很大偏差；Tsvankin等提出的弱各向异性非双曲时距公式在ε-δ<0时误差增大；Alkhalifah等提出的非双曲时距公式在大炮检距任意各向异性强度下都具有较高的精度，适于在实际资料处理中应用.     

Optimization approach to automatic first arrival picking for three‐component three‐dimensional vertical seismic profiling data

A new, adaptive multi‐criteria method for accurate estimation of three‐component three‐dimensional vertical seismic profiling of first breaks is proposed. Initially, we manually pick first breaks for the first gather of the three‐dimensional borehole set and adjust several coefficients to approximate the first breaks wave‐shape parameters. We then predict the first breaks for the next source point using the previous one, assuming the same average velocity. We follow this by calculating an objective function for a moving trace window to minimize it with respect to time shift and slope. This function combines four main properties that characterize first breaks on three‐component borehole data: linear polarization, signal/noise ratio, similarity in wave shapes for close shots and their stability in the time interval after the first break. We then adjust the coefficients by combining current and previous values. This approach uses adaptive parameters to follow smooth wave‐shape changes. Finally, we average the first breaks after they are determined in the overlapping windows. The method utilizes three components to calculate the objective function for the direct compressional wave projection. An adaptive multi‐criteria optimization approach with multi three‐component traces makes this method very robust, even for data contaminated with high noise. An example using actual data demonstrates the stability of this method.     

Application of 3D vertical seismic profile multi‐component data to tight gas sands‡

Due to the complicated geophysical character of tight gas sands in the Sulige gasfield of China, conventional surface seismic has faced great challenges in reservoir delineation. In order to improve this situation, a large‐scale 3D‐3C vertical seismic profiling (VSP) survey (more than 15 000 shots) was conducted simultaneously with 3D‐3C surface seismic data acquisition in this area in 2005. This paper presents a case study on the delineation of tight gas sands by use of multi‐component 3D VSP technology. Two imaging volumes (PP compressional wave; PSv converted wave) were generated with 3D‐3C VSP data processing. By comparison, the dominant frequencies of the 3D VSP images were 10–15 Hz higher than that of surface seismic images. Delineation of the tight gas sands is achieved by using the multi‐component information in the VSP data leading to reduce uncertainties in data interpretation. We performed a routine data interpretation on these images and developed a new attribute titled ‘Centroid Frequency Ratio of PSv and PP Waves’ for indication of the tight gas sands. The results demonstrated that the new attribute was sensitive to this type of reservoir. By combining geologic, drilling and log data, a comprehensive evaluation based on the 3D VSP data was conducted and a new well location for drilling was proposed. The major results in this paper tell us that successful application of 3D‐3C VSP technologies are only accomplished through a synthesis of many disciplines. We need detailed analysis to evaluate each step in planning, acquisition, processing and interpretation to achieve our objectives. High resolution, successful processing of multi‐component information, combination of PP and PSv volumes to extract useful attributes, receiver depth information and offset/ azimuth‐dependent anisotropy in the 3D VSP data are the major accomplishments derived from our attention to detail in the above steps.     

Receiver function method in reflection seismology

The receiver function method was originally developed to analyse earthquake data recorded by multicomponent (3C) sensors and consists in deconvolving the horizontal component by the vertical component. The deconvolution process removes travel path effects from the source to the base of the target as well as the earthquake source signature. In addition, it provides the possibility of separating the emergent P and PS waves based on adaptive subtraction between recorded components if plane waves of constant ray parameters are considered. The resulting receiver function signal is the local PS wave's impulse response generated at impedance contrasts below the 3C receiver.We propose to adapt this technique to the wide‐angle multi‐component reflection acquisition geometry. We focus on the simplest case of land data reflection acquisition. Our adapted version of the receiver function approach consists in a multi‐step procedure that first removes the P wavefield recorded on the horizontal component and next removes the source signature. The separation step is performed in the τ?p domain while the source designature can be achieved in either the τ?p or the t?x domain. Our technique does not require any a priori knowledge of the subsurface. The resulting receiver function is a pure PS‐wave reflectivity response, which can be used for amplitude versus slowness or offset analysis. Stack of the receiver function leads to a high‐quality S wave image.     

三分量数字检波器在ZY油田应用效果分析

在水平分量的炮集记录中,数字三分量检波器直达波的椭圆率明显呈现出线性偏振的形态,并具有较高的分辨率,而模拟三分量检波器直达波的线性偏振不清晰,抗干扰能力差.在频谱分析和能量曲线分析的对比中,数字检波器的性能也高于模拟检波器.在PP波弱波阻抗界面的沙河组沙三12段,PS波表现出强反射界面,而且构造幅度比PP波剖面明显.在PS波剖面上,剖面深处两边的构造联系清晰可见,而在PP波剖面上很难解释剖面两边构造之间的关系,表明PS波资料对深层界面的成像能力有时不比PP波能力差,多分量资料在构造勘探上具有互补性.在沙四层底钻遇气层的马62井和马11井的测线CDP位置上,PP波和PS波的振幅比剖面上表现出低振幅比异常区,与钻井结果吻合.用多分量衰减系数或衰减比剖面,能够识别砂岩含气/水产生所产生的高频衰减异常.多分量地震属性及其比值剖面是储层含气异常指示的有效参数.     

局部和全局边界条件的组合

时域局部人工边界是以行波概念为基础的,原则上不能用于模拟零频分量（消减运动）。Ｗｏｌｆ提出了能同时模拟零频分量和行波的双渐近－多向透射边界。本文指出,在这一组合中的双渐近边界的高频渐近项是不合理的,还查能降低模拟精度。据此,本文提出一种概念清楚、简便实用的组合方案－多次射透－静力边界,并通过数值试验对此方案进行检验。     

Unified multi‐depth‐level field decomposition

Wavefield decomposition forms an important ingredient of various geophysical methods. An example of wavefield decomposition is the decomposition into upgoing and downgoing wavefields and simultaneous decomposition into different wave/field types. The multi‐component field decomposition scheme makes use of the recordings of different field quantities (such as particle velocity and pressure). In practice, different recordings can be obscured by different sensor characteristics, requiring calibration with an unknown calibration factor. Not all field quantities required for multi‐component field decomposition might be available, or they can suffer from different noise levels. The multi‐depth‐level decomposition approach makes use of field quantities recorded at multiple depth levels, e.g., two horizontal boreholes closely separated from each other, a combination of a single receiver array combined with free‐surface boundary conditions, or acquisition geometries with a high‐density of vertical boreholes. We theoretically describe the multi‐depth‐level decomposition approach in a unified form, showing that it can be applied to different kinds of fields in dissipative, inhomogeneous, anisotropic media, e.g., acoustic, electromagnetic, elastodynamic, poroelastic, and seismoelectric fields. We express the one‐way fields at one depth level in terms of the observed fields at multiple depth levels, using extrapolation operators that are dependent on the medium parameters between the two depth levels. Lateral invariance at the depth level of decomposition allows us to carry out the multi‐depth‐level decomposition in the horizontal wavenumber–frequency domain. We illustrate the multi‐depth‐level decomposition scheme using two synthetic elastodynamic examples. The first example uses particle velocity recordings at two depth levels, whereas the second example combines recordings at one depth level with the Dirichlet free‐surface boundary condition of zero traction. Comparison with multi‐component decomposed fields shows a perfect match in both amplitude and phase for both cases. The multi‐depth‐level decomposition scheme is fully customizable to the desired acquisition geometry. The decomposition problem is in principle an inverse problem. Notches may occur at certain frequencies, causing the multi‐depth‐level composition matrix to become uninvertible, requiring additional notch filters. We can add multi‐depth‐level free‐surface boundary conditions as extra equations to the multi‐component composition matrix, thereby overdetermining this inverse problem. The combined multi‐component–multi‐depth‐level decomposition on a land data set clearly shows improvements in the decomposition results, compared with the performance of the multi‐component decomposition scheme.     

Elimination of the water-layer response from multi-component source and receiver marine electromagnetic data

This paper presents the theory to eliminate from the recorded multi‐component source, multi‐component receiver marine electromagnetic measurements the effect of the physical source radiation pattern and the scattering response of the water‐layer. The multi‐component sources are assumed to be orthogonally aligned above the receivers at the seabottom. Other than the position of the sources, no source characteristics are required. The integral equation method, which for short is denoted by Lorentz water‐layer elimination, follows from Lorentz' reciprocity theorem. It requires information only of the electromagnetic parameters at the receiver level to decompose the electromagnetic measurements into upgoing and downgoing constituents. Lorentz water‐layer elimination replaces the water layer with a homogeneous half‐space with properties equal to those of the sea‐bed. The source is redatumed to the receiver depth. When the subsurface is arbitrary anisotropic but horizontally layered, the Lorentz water‐layer elimination scheme greatly simplifies and can be implemented as deterministic multi‐component source, multi‐component receiver multidimensional deconvolution of common source gathers. The Lorentz deconvolved data can be further decomposed into scattering responses that would be recorded from idealized transverse electric and transverse magnetic mode sources and receivers. This combined electromagnetic field decomposition on the source and receiver side gives data equivalent to data from a hypothetical survey with the water‐layer absent, with idealized single component transverse electric and transverse magnetic mode sources and idealized single component transverse electric and transverse magnetic mode receivers. When the subsurface is isotropic or transverse isotropic and horizontally layered, the Lorentz deconvolution decouples into pure transverse electric and transverse magnetic mode data processing problems, where a scalar field formulation of the multidimensional Lorentz deconvolution is sufficient. In this case single‐component source data are sufficient to eliminate the water‐layer effect. We demonstrate the Lorentz deconvolution by using numerically modeled data over a simple isotropic layered model illustrating controlled‐source electromagnetic hydrocarbon exploration. In shallow water there is a decrease in controlled‐source electromagnetic sensitivity to thin resistors at depth. The Lorentz deconvolution scheme is designed to overcome this effect by eliminating the water‐layer scattering, including the field's interaction with air.     

Elastic full waveform inversion of multi‐component OBC seismic data

Elastic full waveform inversion of seismic reflection data represents a data‐driven form of analysis leading to quantification of sub‐surface parameters in depth. In previous studies attention has been given to P‐wave data recorded in the marine environment, using either acoustic or elastic inversion schemes. In this paper we exploit both P‐waves and mode‐converted S‐waves in the marine environment in the inversion for both P‐ and S‐wave velocities by using wide‐angle, multi‐component, ocean‐bottom cable seismic data. An elastic waveform inversion scheme operating in the time domain was used, allowing accurate modelling of the full wavefield, including the elastic amplitude variation with offset response of reflected arrivals and mode‐converted events. A series of one‐ and two‐dimensional synthetic examples are presented, demonstrating the ability to invert for and thereby to quantify both P‐ and S‐wave velocities for different velocity models. In particular, for more realistic low velocity models, including a typically soft seabed, an effective strategy for inversion is proposed to exploit both P‐ and mode‐converted PS‐waves. Whilst P‐wave events are exploited for inversion for P‐wave velocity, examples show the contribution of both P‐ and PS‐waves to the successful recovery of S‐wave velocity.     

YC地区转换波数据处理及裂隙预测

裂隙介质的横波响应能提供对裂隙走向与密度一种直接测量.为了探测YC地区裂隙型储层特征,在该区进行了纵波源三分量数据采集.经过精细地处理,YC地区转换波成像质量得到了改善.以此为基础,对该区的裂隙特征进行了预测.预测结果与本区的应力场分析和钻井揭示的裂隙特征比较吻合.本文所介绍的一些关键转换波资料处理技术、解释方法和分析结论对今后陆地转换波勘探有一定借鉴作用.