基于局部斜率属性的VSP波场分离研究

基于垂直地震剖面(VSP)资料中上、下行波视速度的差异,利用地震剖面同相轴局部斜率属性参数,提出了一种分离上、下行波场的新方法.首先利用Fourier变换初步分离上、下行波场,然后利用平面波分解滤波器(Plane Wave Destruction (PWD) Filter)技术估计初始分离波场的同相轴局部斜率属性参数,在此基础上对VSP原始资料波场分离.该方法是一种时间域最小平方优化分离波场的方法,不存在其他滤波方法阈值滤波器边界的影响,减少了因镶边问题带来的假象.模拟和实际资料处理结果表明,该方法与传统方法相比,分离出的上、下行波噪声假象少,振幅保持好,更好地消除了上、下行波的相互影响.     

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

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

一种稳定的VSP反Q滤波方法

VSP-垂直地震剖面法勘探技术能够提供准确的速度参数、层位标定、井孔周围的构造、岩性及储层的分布范围,为利用地面地震反射信息进行构造精细解释、储层横向预测和油藏描述提供可靠的资料依据.但是地震波在地下介质中传播会产生吸收衰减,降低VSP资料的信噪比和分辨率.本文提出了一种针对VSP资料的稳定的反Q滤波方法.先分离波场,然后对下行波和上行波分别进行反Q滤波补偿.用模型数据和实际VSP资料进行试算,都取得了很好的效果.     

相对保幅的角度域VSP逆时偏移(英文)

本文介绍了一种改进的角度域VSP逆时偏移方法。对VSP逆时偏移中的逆推公式进行了改进,为方便数值计算出相对保幅的角度域共成像点道集(ADCIGs)。此外VSP记录到的波场信息丰富,包括上行波场、下行波场和直达波场等,本文分析了这些波场的响应特征,发现直达波和下行波在角度域共成像点道集(ADCIGs)上都产生了成像噪音,直达波产生的噪音尤为严重。把该方法用于我国西部地区实际观测的VSP资料,不仅获得相对保幅角度域共成像点道集(ADCIGs),而且压制了成像噪音。通过数值模型试算,实际资料的应用验证了该方法的实用性与有效性,从而为VSP偏移速度分析、VSP AVA/AVO分析和反演等提供可靠的基础资料。     

奇异值分解(SVD)实现地震波场分离与去噪新思路

依据不同性质的地震信号（反射波、折射波、直达波、面波、VSP上\下行波、多次波、随机干扰等）之间在运动学、视速度和相干性上的差异,借助某种数学变换、SVD分解与重构联合、时域与频域结合的方式,通过这些间接的处理手段,把要提取的目标信号或要剔除的干扰信号转换到一种相干性更好的空间域中,再进行SVD分解与重构,最充分利用SVD滤波技术特点,实现地震波场分离与去噪,而不是直接对信号进行SVD分解与重构来实现地震波场分离与去噪,这样做可有效地避免以往对SVD波场分离与去噪技术应用空间狭窄、有效信号损失严重等缺陷性.     

SVD地震波场分离与去噪技术在鄂尔多斯盆地地震资料处理中的应用

针对鄂尔多斯盆地地区地震资料自身的特点,本文依据地震信号传播的特性和波场之间的差异,通过多域(多种线性变换域和频域)奇异值分解(SVD),然后提取目标信号的奇异值重构地震信号的方法,实现地震波场分离与去噪处理.与传统SVD地震波场分离与去噪技术相比较,该方法的目的性更强,直接针对感兴趣的地震信号成分进行SVD波场分离与去噪,在提高地震资料信噪比的同时,确保了信号的高保真度和分辨率;同时,避免了以往SVD技术应用空间狭窄,有效信号损失严重等缺陷性.从实际资料处理结果来看,取得了较好的效果.     

多方向正交多项式变换压制多次波

提出一种基于Radon 变换和正交多项式变换的多方向正交多项式变换压制多次波方法.抛物Radon变换对不同曲率方向的同相轴叠加,根据速度差异区分一次波和多次波,但Radon反变换会损伤振幅特性,不利于AVO分析.多方向正交多项式变换在Radon变换(某一曲率方向的零阶特性)的基础上,利用正交多项式变换进一步分析同相轴的高阶多项式分布特性,用正交多项式谱表征同相轴AVO特性;根据一次波和多次波速度差异和同相轴能量分布特征实现多次波压制.该方法的优点是仅用一个曲率参数就可描述同相轴剩余时差参数,提高了一次波和多次波的剩余时差分辨率.实验结果表明,该方法可以有效压制多次波并保留一次波AVO特性.     

基于Marchenko理论一步法压制层间多次波

传统基于Marchenko理论自聚焦法压制层间多次波技术需要初始下行聚焦函数的估计,再进行多次波压制的自适应相减.论文基于传统方法的研究成果,提出了一种无需先验信息和自适应相减步骤的层间多次波压制方法.当聚焦点在实际反射层的反射点上时,上行格林函数为聚焦点上的脉冲源在表面所产生的上行波场,该反射层的反射将是上行格林函数的第一个同相轴.通过将格林函数与下行直达波场进行褶积,可将反射层处的所有焦点重新计算到地表,得到只含一次反射波的波场.利用初始下行聚焦函数和初始直达下行格林函数互为逆函数的关系,引入delta函数,将特定深度上所有聚焦点的聚焦函数投影到采集面上的检波器位置,可只利用地表获取的地震记录一步计算得到一次反射波.水平层状介质模型和复杂的含高速薄透镜地质体的介质模型中,进行了多次波压制的数值试验,通过与传统基于Marchenko理论自聚焦方法的压制效果进行对比分析,显示论文给出的一步法压制层间多次波的效果更好,进而验证了所提出方法的有效性和高精度的特点.     

不同参数对地震干涉实现VSP向SWP波场重构的影响

地震干涉是近几年勘探地球物理领域研究的一个热点方向,它是波场重构的有力工具.然而,地震干涉往往引入虚假同相轴,影响波场重构的质量.为进一步分析虚假同相轴产生的原因并改善波场重构的质量,本文基于稳相分析,详细探讨子波主频、激发炮数、检波器埋深、检波器间距和地层倾角五个因素对波场重构的影响.模型结果表明,震检排列方式和地层倾角等因素通过改变有限震检范围内稳相点的位置来影响波场重构的质量.将垂直地震剖面VSP中的下行直达波与下行反射波进行地震干涉处理,可有效重构来自高陡反射面的下行反射波,将传统的VSP转化为单井地震剖面SWP.直接利用重构SWP波场进行成像,不仅扩大传统VSP的成像范围,而且避免常规勘探面临的静校正和近地表速度建模问题,为高陡构造成像提供一种新方法.     

基于零偏VSP三种Q值反演方法对比分析及应用

准确提取Q值是研究地层吸收衰减特性的关键.对三种Q值反演方法(质心频率法、振幅衰减法、频谱比法)进行对比分析,旨在为衰减特性的求取提供参考.针对零偏VSP数据进行计算,对比总结薄层、频带宽度及低衰减层、波场成分、界面干扰等条件下Q值反演与层位揭示的准确性及差异性:对于相对较薄的层位,质心频率法几乎能准确揭示所有地层,而其他两种方法在薄层分界面处出现异常,误差超过200%;高频成分对频率域方法影响较大,尤其是对于低衰减层的反演;只有下/上行波场时三种方法计算结果相似,全波场时质心频率法效果较好;反射界面的存在会对三种方法造成干扰,反演值在界面处出现跳跃.对比分析三种方法在实际VSP数据应用可知:Q值分层比速度曲线更为敏感,质心频率法Q值反演曲线与地质分层吻合程度最好.     

基于改进线性Radon变换的井孔地震上下行波场分离方法

反射波场分离是井孔地震资料处理中极其重要的一个环节,波场分离的质量直接影响成像结果的精度.不管是VSP还是井间地震资料,其反射波时距曲线都近似直线型,根据这一特征,本文提出一种改进的线性Radon变换方法来进行井孔资料的反射波上下行波场分离.该方法基于频率域线性Radon变换,通过引入一个新的变量λ来消除变换算子对频率的依赖性,避免了求取每一频率分量对应的不同变换算子,显著降低了计算成本;文中在求解该方法对应的最小二乘问题时,引入了发展较为成熟的高分辨率Radon变换技术来进一步提高波场分离的精度.采用本文方法进行井孔地震资料的上下行波场分离可以在保证分离精度的前提下有效地提高计算效率.根据上下行波在λ-f域内分布的特殊性,设计简单的滤波算子就可实现上下行波场的分离.最后通过合成数据试算以及实际资料处理（VSP数据和井间地震数据）验证了该方法的可行性和有效性.     

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

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

VERTICAL SEISMIC PROFILING IN COAL*

The intellection of seismic wave propagation in coal measures demands direct observation of the wavefield progression. Two vertical seismic profiles with high spatial and temporal sampling, were recently recorded in the Sydney Basin coalfields as part of an experimental coal seismic program. Static corrections and interval velocities were obtained by an automated system to determine first kicks and pulse rise times. Upgoing and downgoing waves were separated in the f—k-plane using a novel technique of contour slice filtering. The isolated upgoing waves clearly display reflections from the major coal seams within the stratigraphic sequence. The downgoing wave spectra were subjected to attenuation analysis. The deduced specific quality factor Q for Permian coal measure rocks lies in the range 20–70. Similar estimates were obtained in the time domain from measurements of pulse broadening. Synthetic VSP seismograms, computed using an exact recursive formulation, are an indispensable aid to interpretation. They illustrate the filtering effects of coal seams and sequences, and the effects of the contribution of internal and free-surface multiple reflections in the recorded wavetrains.     

基于格林函数理论的波场预测和鬼波压制方法

作为一种特殊的噪声,鬼波对一次波的波形及频带宽度产生极大的影响,鬼波压制是提高海上地震资料分辨率及保真度的重要因素.以格林公式为基础,详细论述了基于格林函数理论的鬼波压制方法,在不需要地下介质信息的条件下,进行地震数据驱动鬼波压制,并根据"Double Dirichlet"(双狄利克雷)边界条件,预测压力波场和垂直速度波场.建立了基于格林函数理论鬼波压制的处理流程,数值模拟和实际资料处理结果表明,基于格林函数理论鬼波压制方法在很好地去除鬼波的同时极大地拓宽了地震资料的频带,尤其提升了低频端能量,有利于后续资料的处理解释.     

Marine seismic wavefield measurement to remove sea-surface multiples

We propose a new method for removing sea-surface multiples from marine seismic reflection data in which, in essence, the reflection response of the earth, referred to a plane just above the sea-floor, is computed as the ratio of the plane-wave components of the upgoing wave and the downgoing wave. Using source measurements of the wavefield made during data acquisition, three problems associated with earlier work are solved: (i) the method accommodates source arrays, rather than point sources; (ii) the incident field is removed without simultaneously removing part of the scattered field; and (iii) the minimum-energy criterion to find a wavelet is eliminated. Pressure measurements are made in a horizontal plane in the water. The source can be a conventional array of airguns, but must have both in-line and cross-line symmetry, and its wavefield must be measured and be repeatable from shot to shot. The problem is formulated for multiple shots in a two-dimensional configuration for each receiver, and for multiple receivers in a two-dimensional configuration for each shot. The scattered field is obtained from the measurements by subtracting the incident field, known from measurements at the source. The scattered field response to a single incident plane wave at a single receiver is obtained by transforming the common-receiver gather to the frequency–wavenumber domain, and a single component of this response is obtained by Fourier transforming over all receiver coordinates. Each scattered field component is separated into an upgoing wave and a downgoing wave using the zero-pressure condition at the water-surface. The upgoing wave may then be expressed as a reflection coefficient multiplied by the incident downgoing wave plus a sum of scattered downgoing plane waves, each multiplied by the corresponding reflection coefficient. Keeping the upgoing scattered wave fixed, and using all possible incident plane waves for a given frequency, yields a set of linear simultaneous equations for the reflection coefficients which are solved for each plane wave and for each frequency. To create the shot records that would have been measured if the sea-surface had been absent, each reflection coefficient is multiplied by complex amplitude and phase factors, for source and receiver terms, before the five-dimensional Fourier transformation back to the space–time domain.     

用于零偏移距VSP资料的自适应波形反演方法研究

提出一种利用零偏移距VSP资料初至下行波(即直达波)反演介质品质因子Q及层速度V等参数的方法, 称为自适应时域波形反演法(ATWI). 为了充分地利用有效信息, 该方法根据实际VSP资料的信噪比及直达波与上行波干涉的程度,自适应最大限度地选取未受干扰的初至波片段,并用该片段构造目标函数; 通过恰当地构造数据加权矩阵提高目标函数对Q值变化的敏感性;为克服非线性反演的病态问题,采用近来发展的乘性正则化方法,并通过约束条件限制待求参数的取值范围;文中推导出了雅可比矩阵各元素的解析表达式,从而减小了反问题的计算量.合成数据反演结果表明,与谱比值法和子波包络峰值瞬时频率法相比较,ATWI法受上行波影响相对较小、抗噪性能更强.实际资料算例进一步证明了ATWI方法的有效性.     

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.     

Joint viscoacoustic waveform inversion with the separated upgoing and downgoing wavefields in zero-offset vertical seismic profile data

The attenuation of seismic waves propagating in reservoirs can be obtained accurately from the data analysis of vertical seismic profile in terms of the quality-factor Q. The common methods usually use the downgoing wavefields in vertical seismic profile data. However, the downgoing wavefields consist of more than 90% energy of the spectrum of the vertical seismic profile data, making it difficult to estimate the viscoacoustic parameters accurately. Thus, a joint viscoacoustic waveform inversion of velocity and quality-factor is proposed based on the multi-objective functions and analysis of the difference between the results inverted from the separated upgoing and downgoing wavefields. A simple separating step is accomplished by the reflectivity method to obtain the individual wavefields in vertical seismic profile data, and then a joint inversion is carried out to make full use of the information of the individual wavefields and improve the convergence of viscoacoustic full-waveform inversion. The sensitivity analysis of the different wavefields to the velocity and quality-factor shows that the upgoing and downgoing wavefields contribute differently to the viscoacoustic parameters. A numerical example validates our method can improve the accuracy of viscoacoustic parameters compared with the direct inversion using full wavefield and the separate inversion using upgoing or downgoing wavefield. The application on real field data indicates our method can recover a reliable viscoacoustic model, which helps reservoir appraisal.     

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

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