时空移动成像条件及偏移速度分析

首先比较了深度聚焦速度分析和剩余曲率速度分析中的成像条件,然后通过时空移动成像条件得到了时移偏移距域共成像点道集和时移角度域共成像点道集.基于时移角度域共成像点道集,统一了偏移速度分析中通常应用的两个偏移速度判断准则：深度聚焦准则和成像道集拉平准则.最后基于时移角度域共成像点道集,推导了速度更新公式,并设计了速度分析流程.合成数据和实际地震资料上的测试证明了方法的可行性和有效性.     

库车坳陷复杂高陡构造地震成像研究

复杂构造地震成像主要取决于叠前地震数据品质、偏移速度可靠性和偏移算子成像精度. 库车坳陷异常复杂的近地表条件导致极低信噪比的地震采集数据. 该区逆冲推覆高陡构造刺穿盐体大面积分布, 盐层厚度变化大、顶底面形态复杂, 盐下断裂带破碎、小断块发育, 形成异常复杂的地震成像问题. 本文重点研究三个关键环节:(1)精细的叠前地震预处理研究: 根据该区地震地质复杂性和地震资料特征, 采用一些新的方法技术和技术组合从振幅与时移的大、中、小尺度变化三个层次来解决资料信噪比问题, 重建深部反射信号; (2)三级偏移速度分析研究:利用库车坳陷盐刺穿逆冲推覆构造建模理论及变速成图配套技术解决叠前时间偏移速度场时深转换问题,利用井约束低频速度地震迭代反演技术解决连井层速度场与偏移速度场的融合问题,实现从DMO速度分析、叠前时间偏移速度分析到叠前深度偏移速度分析的有机衔接,建立拓扑结构相对保持的叠前深度偏移速度模型;(3)基于退化Fourier偏移算子的半解析波动方程叠前时间和深度偏移研究, 极大地改善了地震偏移过程中高波数波的成像问题. 通过对库车坳陷大北、博孜、却勒、西秋4和西秋10等复杂高陡构造的叠前时间和深度偏移地震成像处理,取得了较好的应用效果.     

基于共聚焦点道集的叠前深度偏移

共聚焦点(CFP)偏移技术是一种基于等时原理,将Kirchhoff积分法的一步偏移分两步聚焦（即激发聚焦和检波聚焦）来完成的叠前地震成像方法.该方法借助于逆时聚集算子和共聚焦点道集来实现叠前偏移成像.基于共聚焦点道集的叠前深度偏移是把基于共炮集的深度偏移的算法引入到CFP技术上来,基于波场延拓的理论来实现偏移成像,该方法首先生成共聚焦点道集,然后基于面炮合成的理论合成聚焦震源,最终通过相关成像来实现叠前偏移成像.该方法选取较少的聚焦点就可以实现对于地下构造的偏移成像,和炮域波动方程偏移相比,其计算效率得到了提高.通过模型试算和实际资料的试处理,验证了该方法在实现叠前深度偏移成像上的正确性和有效性.     

双平方根方程三维叠前深度偏移

从双平方根（DSR）形式的波动方程出发，基于沉降观测概念和地震波扰动理论，介绍了深度域的DSR全偏移算子及共成像道集的生成方法. 根据三维地震数据的方位角特征，通过对全偏移算子的稳相近似，依次导出了适应于零方位角道集、Cross line共偏移距道集以及共偏移距矢量道集的偏移算子. 理论分析与合成数据的数值试验表明，DSR全偏移算子、共方位角偏移算子对介质速度变化的适应性很强，而其余两种偏移算子仅适用于缓变速情况.     

地震波探测地质构造复杂性的定量分析方法

波在传播过程中与非均匀介质的相互作用是多年来地球物理研究的理论问题之一。用地震波探测地下地质结构时,介质的复杂性是相对于地震波长而言的。地震成像的效果与地下速度横向变化和地层的倾角密切相关。基于地震成像算子,研究地下非均匀介质复杂性的定量分析方法,实现定量表征地震成像过程中成像算子的角谱与地下介质地质非均质谱的相干作用。首先,将复杂地质构造表示为慢度非均质谱和倾角非均质谱,以此量化速度横向变化和地层倾角变化分布。其次,通过频散方程分析建立地震成像算子的角谱函数,以此描述其成像精度随慢度和传播角度的变化规律。最后,通过地下复杂构造的地质非均质谱与地震成像算子角谱的点积运算来定义该成像算子对给定地区复杂构造介质的成像效率碾其值越大,说明地震波的探测能力越强,复杂构造的地质复杂性就越小。因此,与地震成像效率相对应的非均匀介质变化的复杂系数可定义为φ=1－η,从而实现对地下复杂介质地震探测复杂性的定量评估。     

利用单程波算子进行叠前时间偏移走时计算(英文)

传统的叠前时间偏移公式都是基于水平层状模型的假设而推导而来,在速度横向变化剧烈的介质中则不能够达到理想的聚焦效果。本文基于李代数积分方程以及拟微分算子等理论推导了非对称走时公式,计算公式由于包含有速度横向导数项,在速度横向变化大的介质中,也能有较高的聚焦效果,而且也为走时计算提供了精确的相对振幅保持权系数。本文对推导的方法进行模型测试并进行实际数据的试算,其结果证明非对称走时方法的成像精度远高于对称走时计算方法。     

面向目标的合成曲面波叠前深度偏移及其应用

面向目标的控制照明叠前深度偏移属于面炮偏移方法的一种特殊情况.通过对目标区域波场照射方向的控制,实现波场在目的层的最佳入射,提高目标区域成像质量.将目的层定义的目标合成算子外推到地表形成地表合成算子,应用地表合成算子进行波面源和记录的合成,最后通过地表合成算子的一系列旋转,实现旋转控制照明叠前深度偏移.本文系统阐述了面向目标的合成曲面波控制照明叠前深度偏移方法的基本原理和实现方法,通过数值模型试算验证了该方法的有效性,最后将其应用到某探区实际资料的处理中,并进行了实际资料的应用效果分析,总结了其在实际资料应用中的特点.     

基于波动方程的广义屏叠前深度偏移

地震波传播算子的计算效率和精度是制约三维叠前深度偏移的关键因素. 广义屏传播算子(GSP, Generalized Screen Propagator)是一种在双域中实现的广角单程波传播算子. 这一方法略去了在非均匀体之间发生的交混回响,但它可以正确处理包括聚焦、衍射、折射和干涉在内的各种多次前向散射现象. 通过背景速度下的相移和扰动速度下的陡倾角校正,广义屏算子能够适应地层速度的强烈横向变化. 这种算子可以直接应用于炮集叠前偏移,通过将广义屏算子作用于双平方根方程,还可以获得一种高效率、高精度的炮检距域叠前深度偏移方法,用于二维共炮检距道集和三维共方位角道集的深度域成像. 本文首先简述了炮检距域广义屏传播算子的理论,进而讨论了共照射角成像(CAI, Common Angle Imaging)条件,由此给出各个不同照射角(炮检距射线参数)下的成像结果,进而得到共照射角像集. 由于照射角和炮检距的对应关系,共照射角像集又为偏移速度分析和AVO(振幅随炮检距变化)分析等提供了有力工具.     

非稳态相移法叠前深度偏移

介绍一种能够适应介质速度横向变化的非稳态相移算子及其叠前深度偏移方法．为了克服常规相移偏移算法中要求速度横向不变的缺点,出现了基于非稳态滤波器理论的非稳态相移算子,即PSPI算子、NSPS算子和SNPS算子,其中SNPS算子是将前二者结合起来的一种对称的非稳态相移算子,它比前二者具有更高的精度和稳定性．为了提高运算速度,基于非稳态相移算子的叠前深度偏移算法采取了分片均匀近似的策略,Marmousi模型的叠前深度偏移结果证明了该算法的可行性和有效性。     

非零偏VSP弹性波叠前逆时深度偏移技术探讨

非零偏VSP地震资料是一种多分量资料,处理非零偏VSP资料,弹性波叠前逆时深度偏移技术无疑是最适合的处理技术.本文从二维各向同性介质的弹性波波动方程出发,研究了对非零偏VSP资料进行叠前逆时深度偏移的偏移算法,讨论了逆时传播过程中的边值问题和数值频散问题及其相应的解决方案;采用求解程函方程计算得到地下各点的地震波初至时间作为成像时间,实现了非零偏VSP资料的叠前逆时深度偏移.最后进行了模型试算和非零偏VSP地震资料的试处理,结果表明该方法不受地层倾角限制,较适用于高陡构造地区或介质横向速度变化较大地区的非零偏VSP地震资料处理.     

An integrated method for resolving the seismic complex near-surface problem

We describe an integrated method for solving the complex near‐surface problem in land seismic imaging. This solution is based on an imaging approach and is obtained without deriving a complex near‐surface velocity model. We start by obtaining from the data the kinematics of the one‐way focusing operators (i.e. time‐reversed Green's functions) that describe propagation between the acquisition surface and a chosen datum reflector using the common‐focus‐point technology. The conventional statics solutions obtained from prior information about the near surface are integrated in the initial estimates of the focusing operators. The focusing operators are updated iteratively until the imaging principle of equal traveltime is fulfilled for each subsurface gridpoint of the datum reflector. Therefore, the seismic data is left intact without any application of time shifts, which makes this method an uncommitted statics solution. The focusing operators can be used directly for wave‐equation redatuming to the respective reflector or for prestack imaging if determined for multiple reflecting boundaries. The underlying velocity model is determined by tomographic inversion of the focusing operators while also integrating any hard prior information (e.g. well information). This velocity model can be used to perform prestack depth imaging or to calculate the depth of the new datum level. We demonstrate this approach on 2D seismic data acquired in Saudi Arabia in an area characterized by rugged topography and complex near‐surface geology.     

Subsurface offset behaviour in velocity analysis with extended reflectivity images

Migration velocity analysis with the constant‐density acoustic wave equation can be accomplished by the focusing of extended migration images, obtained by introducing a subsurface shift in the imaging condition. A reflector in a wrong velocity model will show up as a curve in the extended image. In the correct model, it should collapse to a point. The usual approach to obtain a focused image involves a cost functional that penalizes energy in the extended image at non‐zero shift. Its minimization by a gradient‐based method should then produce the correct velocity model. Here, asymptotic analysis and numerical examples show that this method may be too sensitive to amplitude peaks at large shifts at the wrong depth and to artefacts. A more robust alternative is proposed that can be interpreted as a generalization of stack power and maximizes the energy at zero‐subsurface shift. A real‐data example is included.     

Obstacles to long‐term soil moisture monitoring with heated distributed temperature sensing

The field deployment of a heated distributed temperature sensor (DTS) for over three years has revealed two obstacles to estimating soil moisture (θ) that may hamper subsurface DTS applications as well as use of other subsurface thermal probes. The first observed obstacle was a hysteretic response of the DTS sensor. The relationship between θ and the temperature response (?T) within the cable was not only dependent on θ of the soil, but also on the previous wetting and drying cycles leading to that state. The second observed obstacle was soil structure healing. Soil structure healing causes the relationship between ?T and θ to evolve through time; this calibration curve becomes flatter, or less sensitive, as the surrounding soil makes better contact with the cable. Effects of the hysteretic response of the instrument and soil structure healing are largely the result of small gaps between the cable and soil. These small gaps can be approximated by a contact resistance between the cable and soil. In this article we characterize the occurrence of hysteretic and soil structure healing effects from field data and parameterize contact resistance by simulating heat transfer using a numerical modelling approach Copyright © 2015 John Wiley & Sons, Ltd.     

Wave‐equation migration velocity analysis with time‐shift imaging

Wave‐equation migration velocity analysis is a technique designed to extract and update velocity information from migrated images. The velocity model is updated through the process of optimizing the coherence of images migrated with the known background velocity model. The capacity for handling multi‐pathing of the technique makes it appropriate in complex subsurface regions characterized by strong velocity variation. Wave‐equation migration velocity analysis operates by establishing a linear relation between a slowness perturbation and a corresponding image perturbation. The linear relationship and the corresponding linearized operator are derived from conventional extrapolation operators and the linearized operator inherits the main properties of frequency‐domain wavefield extrapolation. A key step in the implementation is to design an appropriate procedure for constructing an image perturbation relative to a reference image that represents the difference between the current image and a true, or more correct image of the subsurface geology. The target of the inversion is to minimize such an image perturbation by optimizing the velocity model. Using time‐shift common‐image gathers, one can characterize the imperfections of migrated images by defining the focusing error as the shift of the focus of reflections along the time‐shift axis. The focusing error is then transformed into an image perturbation by focusing analysis under the linear approximation. As the focusing error is caused by the incorrect velocity model, the resulting image perturbation can be considered as a mapping of the velocity model error in the image space. Such an approach for constructing the image perturbation is computationally efficient and simple to implement. The technique also provides a new alternative for using focusing information in wavefield‐based velocity model building. Synthetic examples demonstrate the successful application of our method to a layered model and a subsalt velocity update problem.     

Anisotropic migration velocity analysis: Application to a data set from West Africa

Although it is widely recognized that anisotropy can have a significant influence on the focusing and positioning of migrated reflection events, conventional depth imaging methods still operate with isotropic velocity fields. Here, we present an application of a 2D migration velocity analysis (MVA) algorithm, designed for factorized v(x, z) VTI (transversely isotropic with a vertical symmetry axis) media, to an offshore data set from West Africa. By approximating the subsurface with factorized VTI blocks, it is possible to decouple the spatial variations in the vertical velocity from the anisotropic parameters with minimal a priori information. Since our method accounts for lateral velocity variation, it produces more accurate estimates of the anisotropic parameters than those previously obtained with time‐domain techniques. The values of the anellipticity parameter η found for the massive shales exceed 0.2, which confirms that ignoring anisotropy in the study area can lead to substantial imaging distortions, such as mis‐stacking and mispositioning of dipping events. While some of these distortions can be removed by using anisotropic time processing, further marked improvement in image quality is achieved by prestack depth migration with the estimated factorized VTI model. In particular, many fault planes, including antithetic faults in the shallow part of the section, are better focused by the anisotropic depth‐migration algorithm and appear more continuous. Anisotropic depth migration facilitates structural interpretation by eliminating false dips at the bottom of the section and improving the images of a number of gently dipping features. One of the main difficulties in anisotropic MVA is the need to use a priori information for constraining the vertical velocity. In this case study, we successfully reconstructed the time–depth curve from reflection data by assuming that the vertical velocity is a continuous function of depth and estimating the vertical and lateral velocity gradients in each factorized block. If the subsurface contains strong boundaries with jumps in velocity, knowledge of the vertical velocity at a single point in a layer is sufficient for our algorithm to determine all relevant layer parameters.     

Utilization of multiple scattering: the next big step forward in seismic imaging

Surface removal and internal multiple removal are explained by recursively separating the primary and multiple responses at each depth level with the aid of wavefield prediction error filtering. This causal removal process is referred to as “data linearization.” The linearized output (primaries only) is suitable for linear migration algorithms. Next, a summary is given on the migration of full wavefields (primaries + multiples) by using the concept of secondary sources in each subsurface gridpoint. These secondary sources are two‐way and contain the gridpoint reflection and the gridpoint transmission properties. In full wavefield migration, a local inversion process replaces the traditional linear imaging conditions. Finally, Marchenko redatuming is explained by iteratively separating the full wavefield response from above a new datum and the full wavefield response from below a new datum. The redatuming output is available for linear migration (Marchenko imaging) or, even better, for full wavefield migration. Linear migration, full wavefield migration, and Marchenko imaging are compared with each other. The principal conclusion of this essay is that multiples should not be removed, but they should be utilized, yielding two major advantages: (i) illumination is enhanced, particularly in the situation of low signal‐to‐noise primaries; and (ii) both the upper side and the lower side of reflectors are imaged. It is also concluded that multiple scattering algorithms are more transparent if they are formulated in a recursive depth manner. In addition to transparency, a recursive depth algorithm has the flexibility to enrich the imaging process by inserting prior geological knowledge or by removing numerical artefacts at each depth level. Finally, it is concluded that nonlinear migration algorithms must have a closed‐loop architecture to allow successful imaging of incomplete seismic data volumes (reality of field data).     

Enabling affordable omnidirectional subsurface extended image volumes via probing

Image gathers as a function of subsurface offset are an important tool for the inference of rock properties and velocity analysis in areas of complex geology. Traditionally, these gathers are thought of as multidimensional correlations of the source and receiver wavefields. The bottleneck in computing these gathers lies in the fact that one needs to store, compute, and correlate these wavefields for all shots in order to obtain the desired image gathers. Therefore, the image gathers are typically only computed for a limited number of subsurface points and for a limited range of subsurface offsets, which may cause problems in complex geological areas with large geologic dips. We overcome increasing computational and storage costs of extended image volumes by introducing a formulation that avoids explicit storage and removes the customary and expensive loop over shots found in conventional extended imaging. As a result, we end up with a matrix–vector formulation from which different image gathers can be formed and with which amplitude‐versus‐angle and wave‐equation migration velocity analysis can be performed without requiring prior information on the geologic dips. Aside from demonstrating the formation of two‐way extended image gathers for different purposes and at greatly reduced costs, we also present a new approach to conduct automatic wave‐equation‐based migration‐velocity analysis. Instead of focusing in particular offset directions and preselected subsets of subsurface points, our method focuses every subsurface point for all subsurface offset directions using a randomized probing technique. As a consequence, we obtain good velocity models at low cost for complex models without the need to provide information on the geologic dips.     

呼图壁气枪数据中面波振幅变化研究

利用新疆呼图壁气枪数据,在引入小波传播子方法测量到时延迟的基础上,进一步引入面波振幅作为监测地下介质变化的一个重要参数。结果表明:相对于到时的突跳,面波振幅能克服气枪源震源强度突变带来的测量不稳定;同时,通过面波到时延迟和振幅变化特征的模拟,发现面波振幅和到时延迟测量存在较好的正相关关系,且可以用介质速度变化所导致的聚焦和散焦效应来解释。面波振幅也可以作为较为稳定的参数来衡量地下介质的波速变化。     

Anisotropic migration of coincident VSP and cross-hole seismic reflection surveys

Seismic depth migration may result in false reflector positioning and destructive interference when an incorrect velocity field is used to convert from time to depth. The assumption of isotropy to describe anisotropic rocks is one major source of error in the velocity model, although individual survey images may not be impaired by such an approximation. When different survey types such as VSP and cross-hole reflection seismics have coincident illumination of the subsurface, it is important not only to produce consistent images upon depth migration, but also to determine a consistent velocity model. Using real data sets as examples, both objectives are successfully achieved when anisotropy is incorporated into the velocity model.     

Near-surface models in Saudi Arabia

A single‐layer model of the near surface throughout the Kingdom of Saudi Arabia is available. While this simple model suffices for most areas and large subsurface structures, it fails in situations where the surface topography is complex, the base of weathering is below the datum, or where the time structural closure is less than the uncertainty in the static correction. In such cases, multiple‐layered models that incorporate velocities derived from analysis of first arrivals picked from seismic shot records have proved to be successful in defining the lateral heterogeneity of the near surface. The additional velocity information obtained from this first‐arrival analysis (direct as well as refracted arrivals) vastly improves the velocity–depth model of the near surface, regardless of the topography. Static corrections computed from these detailed near‐surface velocity models have significantly enhanced subsurface image focusing, thereby reducing the uncertainty in the closure of target structures. Other non‐seismic methods have been used either to confirm qualitatively or to enhance the layer models previously mentioned. Gravity data may be particularly useful in sandy areas to confirm general structure, while geostatistical modelling of vibrator base‐plate attributes has yielded information that enhances the velocity field. In the global context, exploration targets of the oil and gas industry are seeking smaller and lower relief‐time structures. Thus, near‐surface models will need to enhance and integrate these methods, particularly in areas where the assumption of flat‐lying near‐surface layers cannot be met.