Walkaway VSP数据的联合相干成像(英文)

本文将相干成像用于WVSP数据的一次反射波和下行地表多次波的成像中,通过对叠加次数分布的分析,发现一次反射波相干成像的成像范围小但在井位置处的叠加次数高;下行地表多次波相干成像的成像范围大而在井位置处的叠加次数少。基于此,本文尝试利用一种WVSP的联合相干成像方法来解决这一问题。该方法利用WVSP的直达波与下行地表多次波和直达波与一次反射波来进行联合相干成像。这种联合相干成像方法在有效拓宽成像面积的同时增加了覆盖次数(特别是在井旁),使干涉记录中的虚假同相轴得到了有效地压制,提高了成像剖面的信噪比,使成像结果更加可靠。正演模拟得到的WVSP合成记录和野外实际采集到的WVSP数据处理结果均表明这种方法的有效性。     

VSP地震干涉成像及应用研究

受井中检波器串级数局限,垂直地震剖面(VSP)反射波成像范围窄,且不能对井中最浅接收点上方区域有效成像.虽然多次波成像能扩大成像范围,但在实际应用中尚有诸多困难和挑战.本文根据Wapenaar的地震干涉理论,基于上下行波场分离技术,研发了VSP地震干涉成像方法.该方法将VSP自由表面多次波重建为在地表震源位置激发(虚震源)接收的拟地面地震反射波,然后偏移成像,以达到对多次波间接成像的目的.通过数值模型实验,测试了VSP干涉成像的极限分辨率,并讨论了主要采集参数的影响,结果表明:该方法的垂向和水平极限分辨率分别达约10m和20m,且能分辨深度达6500m处的50m×100m溶洞;采用12至24道井中检波器串采集的VSP资料,其干涉成像结果显著优于VSP反射波成像,与相应的地面地震成像效果相当.将本文方法应用于新疆地区采集的VSP资料,结果表明:与VSP反射波成像相比,成像同相轴更加连续,成像范围显著扩大;与地面地震成像相比,成像结果相当,尤其在浅中部甚至更好.新方法不仅无需进行井中接收点静校正,且能显著增大成像范围,有利于成像同相轴的追踪对比、地震属性提取和地质解释,尤其对中国新疆地区深部缝洞型储层的成像,具有广泛的实际应用前景.     

Estimation of Thomsen＇s anisotropic parameters from walkaway VSP and applications

Estimation of Thomsen＇s anisotropic parameters is very important for accuratetime-to-depth conversion and depth migration data processing. Compared with othermethods, it is much easier and more reliable to estimate anisotropic parameters that arerequired for surface seismic depth imaging from vertical seismic profile （VSP） data, becausethe first arrivals of VSP data can be picked with much higher accuracy. In this study, wedeveloped a method for estimating Thomsen＇s P-wave anisotropic parameters in VTImedia using the first arrivals from walkaway VSP data. Model first-arrival travel times arecalculated on the basis of the near-offset normal moveout correction velocity in VTI mediaand ray tracing using Thomsen＇s P-wave velocity approximation. Then, the anisotropicparameters 0 and e are determined by minimizing the difference between the calculatedand observed travel times for the near and far offsets. Numerical forward modeling, usingthe proposed method indicates that errors between the estimated and measured anisotropicparameters are small. Using field data from an eight-azimuth walkaway VSP in TarimBasin, we estimated the parameters 0 and e and built an anisotropic depth-velocity modelfor prestack depth migration processing of surface 3D seismic data. The results showimprovement in imaging the carbonate reservoirs and minimizing the depth errors of thegeological targets.     

双井Walkaway VSP技术在辽河坳陷曙光高陡构造地区的应用研究

针对辽河坳陷西部凹陷曙光潜山高陡构造带顶界面断裂成像精度差、内幕地层识别难的问题,开展双井Walkaway VSP应用研究.利用大阵列检波器接收、炸药与震源联合采集方法,获得了高品质采集数据,是高陡构造成像的基础;子波一致性校正、基于射线追踪的矢量波场分离、VSP数据与地面地震数据联合各向异性速度建模、双井Walkaway VSP成像方法对比等针对性处理技术和波阻抗反演、地震属性提取等综合研究方法的应用,精细刻画了潜山顶界面,新识别出井旁小断层,落实了大倾角潜山构造形态,分析了内幕地层的接触关系,对板岩隔夹层及储层分布有进一步认识,展现了Walkaway VSP技术在高陡构造地区的应用前景.研究得到以下认识:(1)逆时偏移成像方法对高陡潜山成像更具适用性;(2)增加检波器进入潜山内幕的深度,避开潜山高速顶界面的能量屏蔽及滑行波、转换波等干扰,对潜山勘探具有重要意义;(3)双井Walkaway VSP联采时,两井距离要控制在一定范围内,防止两井间距离过大对深层偏移成像造成影响.     

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

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

局域化相空间中的VSP偏移成像方法(英文)

对VSP资料进行偏移成像可提高井附近地下结构的成像分辨率。本文给出了一种基于局域化相空间波场分解的VSP偏移成像方法。此方法采用了基于Gabor-Daubechies紧标架的延拓算子(G-D延拓算子)及其高频渐近形式对相空间波场进行延拓;基于局部平面假设,提出了一种局部角度域相关成像条件。合成和实际VSP资料的偏移成像结果表明,在满足渐近展开的条件下利用G-D延拓算子的高频近似式能够有效的减少计算时间;采用局部角度域相关成像条件能够在不增加计算量的同时,有效减弱VSP成像剖面上的偏移假象。     

Migration scheme for imaging offset VSP data within local phase space

The imaging of offset VSP data in local phase space can improve the image of the subsurface structure near the well. In this paper, we present a migration scheme for imaging VSP data in a local phase space, which uses the Gabor-Daubechies tight framebased extrapolator (G-D extrapolator) and its high-frequency asymptotic expansion to extrapolate wavefields and also delineates an improved correlation imaging condition in the local angle domain. The results for migrating synthetic and real VSP data demonstrate that the application of the high-frequency G-D extrapolator asymptotic expansion can effectively decrease computational complexity. The local angle domain correlation imaging condition can be used to weaken migration artifacts without increasing computation.     

基于伪谱法的VSP逆时偏移及其应用研究

逆时偏移被认为是对地下复杂构造进行成像的精确偏移方法,尤其是能够有效地对回转波、绕射波、多次波等各种波动现象进行成像.近几年来随着并行计算机和存储设备的快速发展,逆时偏移方法备受关注.本文采用伪谱法实现了VSP逆时偏移,该方法不仅实施简便,而且计算效率高,精度好.并运用反周期扩展法来消除伪谱法中特殊的周期性边界效应问题.对VSP绕射点模型进行试算,分析了因覆盖次数不足在近井区产生的假象问题.对地堑模型和半圆隆起模型也获得了较好的VSP逆时偏移成像效果.并分别对VSP全波波场及分离出的上行波场进行了逆时偏移成像,可明显发现直达波在炮点和检波点位置处收敛成像,也产生了很强的成像噪声.最后对某地区实际观测的VSP资料进行了逆时偏移成像,并与Kirchhoff法VSP偏移结果和地面地震偏移结果进行了对比,显示了VSP逆时偏移在近井区成像上的优势.     

基于walkaway VSP下行P波的TTI各向异性参数反演

本文基于弱各向异性（WA）介质的正反演公式和qP波的坐标变换,推导了利用qP波反演任意倾斜对称轴的横向各向同性（TTI）介质的各向异性参数和对称轴方向的公式.理论和数值实验表明,利用2个相互正交的变井源距垂直地震剖面（walkaway VSP）可以完全确定钻井中TTI介质qP波的3个WA参数和对称轴的2个方向参数.我们完成了几个由不同数量剖面组成的walkaway VSP模拟实验,使用TTI模型和一般各向异性模型对模拟数据进行了反演,证明了反演公式的正确性和可靠性.使用这些公式,对来自Java Sea的由3条剖面组成的walkaway VSP观测数据进行了各向异性反演,获得了钻井中接收点处介质的WA参数.     

CDP mapping and image reconstruction by using offset VSP data

Because zero-offset VSP (Vertical Seismic Profile) data can only provide the information of rock properties and structure in the area around the Fresnel zone within the well, the scheme of VSP with offset was developed to acquire the reflection information away from the borehole in order to widen the range of VSP survey and to improve the precision of imaging.In this paper, we present a new CDP (Common Depth Point) mapping approach to image the reflecting structure by using offset VSP data. For the processing of offset VSP data, we firstly separated the up-going and down-going wave-fields from VSP data by means of F-K filtering technique, and we can calculate the mapping conditions (position and reflecting traveltime for CDP point) in homogeneous media, and then reconstruct the inner structure of the earth.This method is tested by using the offset VSP data which are used to simulate the case of super-deep borehole by means of finite-difference method. The imaged structure matches the real model very well. The results show that the method present here could accurately image the inner structure of the earth if the deviation of initial velocity model from the true model is less than 10%. Finally, we presented the imaged results for the real offset data by using this method.     

Monte-Carlo Bayesian look-ahead inversion of walkaway vertical seismic profiles

We describe a method to invert a walkaway vertical seismic profile (VSP) and predict elastic properties (P‐wave velocity, S‐wave velocity and density) in a layered model looking ahead of the deepest receiver. Starting from Bayes's rule, we define a posterior distribution of layered models that combines prior information (on the overall variability of and correlations among the elastic properties observed in well logs) with information provided by the VSP data. This posterior distribution of layered models is sampled by a Monte‐Carlo method. The sampled layered models agree with prior information and fit the VSP data, and their overall variability defines the uncertainty in the predicted elastic properties. We apply this technique first to a zero‐offset VSP data set, and show that uncertainty in the long‐wavelength P‐wave velocity structure results in a sizable uncertainty in the predicted elastic properties. We then use walkaway VSP data, which contain information on the long‐wavelength P‐wave velocity (in the reflection moveout) and on S‐wave velocity and density contrasts (in the change of reflectivity with offset). The uncertainty of the look‐ahead prediction is considerably decreased compared with the zero‐offset VSP, and the predicted elastic properties are in good agreement with well‐log measurements.     

基于虚源估计的复杂上覆地层下地震相干成像

在上覆地层比较复杂的情况下,常规地震勘探方法常常难以得到好的成像.本文研究了基于地震相干避开复杂上覆地层对地震波的影响,利用VSP数据估计地震虚源直接对目的地层进行成像的方法.在地震相干成像过程中,震源子波对分辨率有比较大的影响,尤其是存在薄层的条件下,两个非常近的反射同相轴将无法辨认.利用估计出虚源地震子波的性质,对...     

变井源距垂直地震剖面各向异性参数反演

利用微扰理论推导了弱各向异性(WA)介质参数正反演计算的基本公式;给出了在已知慢度矢量的一个分量和偏振矢量情况下确定WA参数的方法.如果这一方法被用于单一的变井源距垂直地震剖面(walkawayVSP)资料,可以获得9个WA参数.这9个WA参数完全决定了qP波和两个qS波在由剖面和钻井所决定的平面内传播的特性.对单条walkawayVSP观测系统数据的产生和WA参数的反演进行了数值模拟计算,对所能确定的WA参数及其可靠性进行了详细的讨论.     

OFFSET SOURCE VSP SURVEYS AND THEIR IMAGE RECONSTRUCTION*

Zero-offset-source VSP surveys provide information about the subsurface only within the Fresnel zone centered at the well. Offsetting the source location moves the reflection zones away from the well thus providing lateral cover. Conventional processing of this type of data gives rise to a distorted image of the subsurface. Using a simple ray-tracing scheme, this image may be reconstructed into the more familiar coordinate system of the surface seismic section. This simple data-independent mapping is based on the assumption of horizontal layering and requires a vertical velocity profile. The technique of placing the source away from the borehole was first applied to the single-offset-source VSP survey. However, data from any survey geometry (such as deviated well with rig source, walkaway VSP, etc.) can be mapped to the coordinate system defined by the appropriate seismic section. To obtain the best results from this type of survey the target area must be defined and simple modeling techniques used to optimize the source location(s). These pre-survey modeling methods may also be used to anticipate—and hence avoid a number of problem areas which experience has highlighted. The data from any VSP survey is the result of a realizable experiment and as such obeys the wave equation. This implies that the wave equation may be used to migrate the data to its true subsurface location. Theoretically, such a process is more secure than ray-tracing techniques, although its practice presents many difficulties.     

True amplitude elastic Gaussian beam imaging of multicomponent walkaway vertical seismic profiling data

We develop the true‐amplitude prestack migration of multicomponent data based on the use of elastic Gaussian beams for walkaway vertical seismic profile (VSP) acquisition systems. It consists in a weighted summation of multishot data with specific weights, computed by tracing elastic Gaussian beams from each imaging point of the target area towards the sources and receivers. Each pair of beams may be connected with either a pair of P‐rays (PP‐image) or the P‐ray towards sources and the S‐ray to receivers (PS‐image) and is uniquely determined by dip (the angle of the bisector between the rays and the vertical direction) and opening (the angle between the rays) angles. Shooting from the bottom towards the acquisition system helps to avoid well‐known troubles, in particular multipathing for the imaging conditions in complex velocity models. The ability to fix the dip angle and implement summation over opening angles leads to the so‐called selective images that contain mostly interfaces with desired slopes. On the other hand, a set of images computed for a range of opening angles by summation over all available dip angles is used as input of an AVO‐like inversion procedure for the recovery of elastic parameters. The feasibility of this imaging procedure is verified by synthetic data for 2D realistic elastic models.     

多测线变偏移距VSP地震各向异性反演

本文从一般弱各向异性介质参数反演中得到了使用两条相互正交的变偏VSP测线计算弱各向异性(WA)参数的反演公式. 如果仅仅使用qP波, 则可以确定9个独立的WA参数, 这9个WA参数可以完全地描述井中接收点在两个剖面内介质的各向异性性质. 通过对走时曲线进行最小二乘条件下的三次样条光滑, 可以获得慢度矢量的垂直分量和水平分量. 如果介质是横向非均匀介质, 则水平分量仅仅被用作反演时的约束条件. 为了获得偏振矢量, 本文引入质心计算方法, 该方法计算简单、 稳定, 而且不存在解的奇异问题. 在寻找与一般各向异性介质最接近的高对称性正交各向异性介质和TTI各向异性介质时, 使用qP波各向异性坐标变换方法和最小二乘求解方法, 得到了与一般各向异性介质最接近的正交各向异性和TTI各向异性参数及其对称轴方向参数的计算公式. 使用这些方法, 对瓜哇海地区布设的两条相互正交的变偏VSP测线数据进行各向异性反演, 获得了井中10个接收点处介质的WA参数. 数值计算和实际资料反演表明, 本文所使用的反演方法能够准确地得到VSP井中接收点处介质的WA参数, 这是地震勘探中研究地壳介质各向异性性质最直接和最可靠的方法.      

克希霍夫法VSP多波联合成像

VSP 转换波跟VSP 纵波或常规地面转换波相比,具有较高的分辨率和信噪比,但传统的VSP成像方法只利用了反射P波信息,而把转换波（反射S波、透射S波）以及透射P波当作影响成像质量的噪音.本文给出了一种VSP共炮点道集多分量地震资料克希霍夫法偏移成像的方法.本方法充分利用了多波（反射P波、反射S波、透射P波、透射S波）信息,根据转换点处四种波同时起跳,能量叠加最大的原理,从接收点分别用向绕射点延拓它们的能量,并将其叠加起来,求得的和最大的一点即反射点.通过模型试算和实际资料处理表明,此法成像精度高,信噪比高,且有利于改善剖面的频率特性.     

Decoupled elastic prestack depth migration

This paper presents a new decoupled form of the formula for common-shot or common-receiver amplitude-preserving elastic prestack depth migration (PreSDM), which can be used for estimating angle-dependent elastic reflection coefficients in laterally inhomogeneous anisotropic media. The multi-shot or multi-receiver extension of this formula is suitable for automated prestack amplitude-versus-angle (AVA) elastic inversion of ocean-bottom cable (OBC), walkaway VSP (WVSP) or standard towed-cable data at any subsurface location. The essence of the theory is a systematic application of the stationary-phase principle and high-frequency approximations to the basic elastic Green's theorem. This leads to nonheuristic explicit wave mode decoupling and scalarization of vector PreSDM. Used in combination, ray-trace and finite-difference (FD) eikonal solvers create a useful tool to calculate accurate Green's function travel time and amplitude maps. Examples of synthetic OBC data and applications to field WVSP data show that the new imaging technique can produce a clear multi-mode elastic image.     

Estimation of effective source signatures from marine VSP data1

We present a simple method for estimating an effective source wavelet from the first arrival in marine vertical seismic profiling (VSP) data. The method, which utilizes the free-space Green's function of the Helmholz equation, is simple and very computer efficient. We show examples from synthetic and real offset and walkaway VSP data. In the synthetic examples, we show that data modelled with the estimated wavelet give small residuals when subtracted from the reference data. In the real data examples, we show that when modelling with the wavelet estimated from the real data, in a smooth macromodel, we obtain a good fit between the first arrivals in the real and modelled data.