高保真VSP波场分离方法研究

三维VSP资料中,各种不同类型的波混杂一起形成复杂的波场.因此,波场分离是三维VSP数据处理关键的第一步.从不同波场的偏振方向和传播方向之差异着手,提出了一种高保真的VSP波场分离方法.首先通过射线追踪和偏振滤波的结合,把复杂波场(分解为简单波场;然后根据简单波场中不同波的传播方向截然相反的特点,进行方向滤波,达到波场分离的目的.实际数据处理表明,与常规波场分离方法相比,本方法大大降低了混波作用以及由此而生的波形畸变.     

平面网格位场数据的空间域非线性曲率滤波方法

剖面测线位场数据的空间域非线性曲率涯波方法能够较好地分离提取剖面到线位场的趋势和剩余异常，本文首次将该方法推广应用到平面网格位场数据情况，为平面网格位场的趋势和剩余异常的分离提取提供了一种新的实用方法。在新疆某地区的实例计算表明，本文提出的平面网格位场数据的空间域非线性曲率滤波方法的应用效果好、计算速度快。     

双相各向同性介质弹性波高精度波场分离数值模拟方法

为了便于研究双相介质固流相混合弹性波场中纵横波波场的传播规律,提出了基于交错网格的Biot双相各向同性介质弹性波动方程高精度波场分离正演数值模拟方法.采用高阶交错网格有限差分法来构建一阶双曲型双相各向同性介质弹性波动方程正演算子实现波场正演,并在每一步递推过程中,分别计算出同相和流相分量相应的散度场(纯纵波场)和旋度场...     

扩散滤波方法在地震资料处理中的应用研究

在地震勘探中,由于地下结构错综复杂,多尺度非均匀的地质体常会形成包含反射波、散射波等在内的复杂的地震波场。传统的成像方法一般只考虑反射波场,忽略了散射波场,这使得细小结构无法准确成像,从而影响对复杂构造的识别。为了对小尺度构造进行准确的地震成像,要将散射波从地震波场中分离出来。在众多波场分离算法中,基于滤波的波场分离方法可以准确提取散射波,提高成像分辨率。本文调研和归纳多种基于滤波处理的地震散射波分离方法,围绕国内外学者在滤波处理波场分离方面的研究成果,总结各种方法的研究进展,并对比和分析各方法的分离效果,最后结合人工智能深度学习的研究趋势,对未来滤波处理散射波分离的发展方向进行展望。

     

多分量地震波波场分解研究

提出了一种有效的多分量地震波波场分解方法, 该方法假设在相邻道集波场的振幅、相位变化不大的情况下, 根据纵波和转换波在传播中相对于不同的极化矢量和慢度得到的垂向和水平向分量, 在最小方差约束下从多分量输入数据中获取纵波和转换波的偏振方向和振幅并重构纵波和转换波, 从而分离纵波和转换波. 这种方法精度较高, 不仅能有效分解一次反射波的多分量波场, 多次反射的纵波和转换波也能被恰当地分解. 该方法不需要输入地层速度等参数信息, 减少了人为误差的影响. 通过对合成资料的处理, 验证了方法的可行性和正确性.     

井间透射-反射波场分离方法及应用

近十几年来,井间地震在石油勘探领域发挥着重要作用.然而井间地震波极其复杂,给同相轴拾取、层析成像、偏移成像等工作带来了困难.因此,波场分离成为井间地震数据处理的一项重要内容.目前,不少学者研究了井间上下行波分离,纵横波分离等多种方法,但尚未出现相关文献研究透射波和反射波的分离方法.本文根据井间透射波和反射波在不同道集(共炮点集或共接收点集)中视速度极性的差异,首次提出了一种"透射.反射"波场分离方法.该方法首先在不同道集中利用τ-p变换进行视速度分离,得到反射上行波、反射下行波、透射上行波、透射下行波四个波场,然后选择对应的波场叠加,从而得到透射波和反射波.其中,透射波为井间层析成像提供了相对清晰的初至同相轴,尤其是横渡初至同相轴,以便进行到时提取;所得反射波为井中偏移成像提供了可靠的反射波数据,比人工切除透射波方法更为准确有效.本文通过理论模型试验和实际数据处理,证明了该方法的有效性和实用性.此外,在本文基础上,再进行纵横波分离处理可以得到更多的波场分量,为井间成像提供更有效的数据.     

航磁水平分量数据的均衡滤波位场边缘增强研究

航磁矢量数据包含异常的强弱与方向信息,能够有效减少反演的多解性,从而提高数据处理与解释的精度和可靠性.因此,我们提出了航磁水平分量数据的均衡滤波位场边缘增强技术,提高了场源边界识别的分辨率.组合模型试验证明此方法不仅能够有效均衡不同埋深的地质体的异常信息,而且能提高对于相近或叠加地质体的识别能力.将此方法应用于新疆启鑫地区实测航磁矢量数据处理,结果表明,此方法使研究区弱异常得到了有效的均衡与增强,获取了更加清晰的位场边缘细节信息,结合研究区地质概况推断出29条断裂构造线,为研究区构造演化的推断与解释提供了更加丰富的参考依据.因此,该方法在应用航磁矢量数据进行位场边缘研究方面具有重要意义,为航磁矢量数据处理与解释提供了新的技术支撑.

     

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

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

     

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

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

用地震波参数预测地震的实践

１９９４年８月２３日及１９９５年５月１５日在新疆巴楚县的东北和西北部两地相继发生了ＭＳ５．２、５．３地震。依据以往地震波参数震例研究结果及其异常指标，笔者在上述两次地震震前对其作了短临预报，实践检验令人鼓舞。     

带阻尼的一维地震反射波成像

本文是文献［８］的继续，文中提出了地震勘探的带阻尼波动方程反演的新方法，应用此反演算法可以通过微机识别地层的波阻抗，这里把地层分成１００层，地震数据含有１０％的误差，且反演计算时间仅仅是正演计算时间的三倍．文献［１］－［４］的计算方法是基于波动方程的混合问题，而文献［８］所提出的反演算法是基于波动方程的哥西问题，许多例子充分说明后者更有效．     

非测震学前兆异常的震情窗口效应

通过对１９７６年唐山７．８级、１９８３年菏泽５．９极和１９８９年大同６．１级．１９９１年大同５．８级四次地震的非测震学前兆异常（短水准、地电、水氡含量、水位等）的研究，发现在这些地震以前，邢台地区的非测震学前兆观测值均有异常出现，显示出该地区不仅有地震活动性的窗口效应，而且有非测震学前兆异常的震情窗口效应。其它地区也有相类似的情况存在。这对区域应力场研究和前兆场特征的研究具有重要意义，对强震的预报     

LINEAR FILTERING OF VLF DIP-ANGLE MEASUREMENTS*

The suggested interpretation technique is based on discrete linear filtering of VLF data. The output of the described filtering results is expressed in terms of an equivalent current density at a specified depth that would cause the measured magnetic field. The most practical six-point filter gives an accuracy of 8%. The filter is an extension of the commonly used Fraser filter to process VLF dip-angle data. Filtering the same data set for various depths gives an idea about the change of current density with depth. Areas with high current-density correspond to good conductors. The conductor dip can also be determined. The use of the technique is illustrated on theoretical and field examples. In all cases a good correlation with original models and other types of geophysical measurements was obtained. As shown in the last example, the filtering technique is also applicable in interpretation of other electromagnetic methods.     

WAVE SCATTERING DECONVOLUTION BY SEISMIC INVERSION*

We propose a wave scattering approach to the problem of deconvolution by the inversion of the reflection seismogram. Rather than using the least-squares approach, we study the full wave solution of the one-dimensional wave equation for deconvolution. Randomness of the reflectivity is not a necessary assumption in this method. Both the reflectivity and the section multiple train can be predicted from the boundary data (the reflection seismogram). This is in contrast to the usual statistical approach in which reflectivity is unpredictable and random, and the section multiple train is the only predictable component of the seismogram. The proposed scattering approach also differs from Claerbout's method based on the Kunetz equation. The coupled first-order hyperbolic wave equations have been obtained from the equation of motion and the law of elasticity. These equations have been transformed in terms of characteristics. A finite-difference numerical scheme for the downward continuation of the free-surface reflection seismogram has been developed. The discrete causal solutions for forward and inverse problems have been obtained. The computer algorithm recursively solves for the pressure and particle velocity response and the impedance log. The method accomplishes deconvolution and impedance log reconstruction. We have tested the method by computer model experiments and obtained satisfactory results using noise-free synthetic data. Further study is recommended for the method's application to real data.     

TRANSFORMATION OF DIPOLE RESISTIVITY SOUNDING MEASUREMENTS OVER LAYERED EARTH BY LINEAR DIGITAL FILTERING *

The technique of linear digital filtering as developed for the direct interpretation of Schlumberger and Wenner soundings (Ghosh 1971) has been applied here for the derivation of the resistivity transform function from the field dipole measurements as the first step in directly interpreting dipole data. Filter coefficients for this transformation have been worked out for the radial-polar, perpendicular and parallel (30°) arrays of dipole sounding. The procedure combines speed with accuracy.     

EFFICIENT MULTICHANNEL FILTERING OF SEISMIC DATA1

Multichannel filters are used to eliminate coherent noise from surface seismic data, for wavefield separation from VSP stacks, and for signal enhancement. Their success generally depends on the choice of the filter parameters and the domain of application. Multichannel filters can be applied to shots (monitors), common-receiver traces, CDP traces and stacked sections. Cascaded applications in these domains are currently performed in the seismic industry for better noise suppression and for signal enhancement. One-step shot-domain filtering is adequate for some applications. However, in practice, cascaded applications in shot-and common-receiver domains usually give better results when the S/N ratio is low. Multichannel filtering after stacking (especially after repeated applications in shot and/or receiver domains) may create undesirable results such as artificial continuations, or smearing and smoothing of small features such as small throw faults and fine stratigraphic details. Consequently, multichannel filtering after stacking must be undertaken with the utmost care and occasionally only as a last resort. Multichannel filters with fan-shaped responses (linear moveout filters) should be applied after NMO correction. These are the filters commonly used in the seismic industry where they have such names as velocity filters, moveout filters, f-k filters and coherency filters. Filtering before NMO correction may result in break-up and flattening especially of those shallow reflection events with relatively higher curvatures and diffractions. NMO correction is needed prior to wavefield separation from VSP stacks for the same practical reasons outlined above whenever source-receiver offsets are involved. Creation of artificial lineup and smearing at the outputs of multichannel filters is presently the common practical concern. Optimum multichannel filters with well-defined pass, reject and transition bands overcome the latter problems when applied before stacking and after NMO correction. The trace dimension of these filters must be kept small to avoid such lineups and the smoothing of small structures. Good results can be obtained with only five traces, but seven traces seems to be a better compromise both in surface and well seismic applications. The so-called f-k filtering and τ-p domain filtering are no exceptions to the above practical considerations. Residual static computations after multichannel filtering also need special consideration. Since multichannel filtering improves spatial continuity, residual static algorithms using local correlation, i.e. nonsurface-consistent algorithms, may be impractical especially after multichannel filtering.     

地震反射波成像的“三基色”方法

基于动态情形介质参数是频率的复函数,利用复小波变换处理地震勘探反射波数据,求出每时刻瞬时谱曲线,它刻划了该时刻所对应的地层介质性质．本文提供一种全新的由反射波数据绘制地下剖面图以取代反褶积的方法,提出用“三基色”方法来提取这条谱线的特征,合成彩色地下剖面图,并利用特制的小波软件对油田实际资料进行研究．     

地震反射波成像的“三基色”方法

基于动态情形介质参数是频率的复函数,利用复小波变换处理地震勘探反射波数据,求出每时刻瞬时谱曲线,它刻划了该时刻所对应的地层介质性质．本文提供一种全新的由反射波数据绘制地下剖面图以取代反褶积的方法,提出用“三基色”方法来提取这条谱线的特征,合成彩色地下剖面图,并利用特制的小波软件对油田实际资料进行研究．     

P^—波定位方法的研究及软件编制

Ｐ＾－波定位方法是针对网内及网缘发生较大震级地震时，地震波记录限幅Ｓ＾＝难以辨认的情况，而编；与的单震相定位程序。该程序可由数字化仪输入地震波到时数据，然后进行定位，定位计算分两步进行，即初定和修定，初定仅根据区域台风提供直达波到时，列出走时方程，将走时方程线性化，解出四个震源参数，为了防止输入的个别错误数据影响定位，采取了图形显示的方法。