再论各向异性介质转换波地震学,第二部分:参数计算研究

在具有垂直对称轴横向各向同性介质中，利用四种参数来确定中间至远偏移距转换波（C-波）动校正。它们是C-波叠加速度Vc2，垂直速度比和有效速度比γ0和γeff以及各向异性参数χeff。我们将这四种参数作为C波叠加速度模型。C-波速度分析的目的就是确定这种叠加速度模型。C-波叠加速度模型Vc2，γ0，γeff，和χeff可以由P-波和C-波反射动校正资料获得。然而错误的传播是C-波反射动校正反演中的严重问题。当前短排列叠加速度由于是从双曲线动校正推算而得，因而其精度不足以为各向异性参数提供有意义的反演值。中间偏移非双曲线动校正不再被人们所勿略，而是可以用一个背景γ加以量化。非双曲线分析通过中间偏移距的γ校正量可以产生Vc2，若数据不含燥音，其误差小于1％。方法稳健，允许γ启始假定值的误差达20％。该方法也适用垂直非均匀各向异性介质。精度的提高使能够用4分量地震资料计算各向异性参数。为此提出了两种工作流程：双扫描和单扫描流程。理论数据和实际数据的应用表明这两种流程得出的结果其精度相似，但是单扫描流程比双扫描更有效。     

Converted-wave seismology in anisotropic media Revisited,Part II: Application to parameter estimation

In transversely isotropic media with a vertical symmetry axis (VTI), the converted-wave (C-wave) moveout over intermediate-to-far offsets is determined by four parameters. These are the C-wave stacking velocity V _C2, the vertical and effective velocity ratios γ ₀and γ _eff, and the anisotropic parameter X _eff. We refer to the four parameters as the C-wave stacking velocity model. The purpose of C-wave velocity analysis is to determine this stacking velocity model. The C-wave stacking velocity model V _C2, γ ₀, γ _geff, and X _eff can be determined from P- and C-wave reflection moveout data. However, error propagation is a severe problem in C-wave reflection-moveout inversion. The current short-spread stacking velocity as deduced from hyperbolic moveout does not provide sufficient accuracy to yield meaningful inverted values for the anisotropic parameters. The non-hyperbolic moveout over intermediate-offsets (x/z from 1.0 to 1.5) is no longer negligible and can be quantified using a background γ. Non-hyperbolic analysis with a γ correction over the intermediate offsets can yield V _C2 with errors less than 1% for noise free data. The procedure is very robust, allowing initial guesses of γ with up to 20% errors. It is also applicable for vertically inhomogeneous anisotropic media. This improved accuracy makes it possible to estimate anisotropic parameters using 4C seismic data. Two practical work flows are presented for this purpose: the double-scanning flow and the single-scanning flow. Applications to synthetic and real data show that the two flows yield results with similar accuracy but the single-scanning flow is more efficient than the double-scanning flow. This work is funded by the Edinburgh Anisotropy Project of the British Geological Survey. First Author Li Xiangyang, he is currently a professorial research seismologist (Grade 6) and technical director of the Edinburgh Anisotropy Project in the British Geological Survey. He also holds a honorary professorship multicomponent seismology at the School of Geosciences, University of Edinburgh. He received his BSc(1982) in Geophysics from Changchun Geological Institute, China, an MSc (1984) in applied geophysics from East China Petroleum Institute (now known as the China University of Petroleum), and a PhD (1992) in seismology from the University of Edinburgh. During 1984–1987, he worked as a lecturer with the East China Petroleum Institute. Since 1991, he has been employed by the British Geological Survey. His research interests include seismic anisotropy and multicomponent seismology.     

再论各向异性介质转换波地震学,第一部分:理论

我们业已研发了计算各向异性、非均质介质中P- SV转换波(C-波)的转换点和旅行时的新理论。据此 可以利用诸如相似性分析、迪克斯模型建模、克契 霍夫求和等常规方法来完成各向异性的处理和各向 异性处理,并使各向异性的处理成为可能。这里将 我们的新发展分作两部分来介绍。第一部分为理 论,第二部分为对速度分析和参数计算的应用。第 一部分理论包括转换点的计算和动校正的分析。     

CONVERSION POINTS AND TRAVELTIMES OF CONVERTED WAVES IN PARALLEL DIPPING LAYERS1

For converted waves, stacking as well as AVO analysis requires a true common reflection point gather which, in this case, is also a common conversion point (CCP) gather. The coordinates of the conversion points for PS or SP waves, in a single homogeneous layer can be calculated exactly as a function of the offset, the reflector depth and the ratio v_p/v_s. An approximation of the conversion point on a dipping interface as well as for a stack of parallel dipping layers is given. Numerical tests show that the approximation can be used for offsets smaller than the depth of the reflector under consideration. The traveltime of converted waves in horizontal layers can be expanded into a power series. For small offsets a two-term truncation of the series yields a good approximation. This approximation can also be used in the case of dipping reflectors if a correction is applied to the traveltimes. This correction can be calculated from the approximated conversion point coordinates.     

VTI介质中P-SV波转换点与各向异性参数关系

对于转换波地震勘探中的转换点位置这个重要问题,提出转换点位置不仅与纵横波速度比,偏移距深度比以及源检距有关,还与地下介质的各向异性的性质有关,计算了忽略地下介质的各向异性影响对转换点的确定带来的严重误差,从而影响地下地质体的精确成像.通过对层状VTI介质中的转换点近似方程的推导过程,引入该方程不同于传统方程的导出是对层状各向同性介质而言,该方程通过引入各向异性参数,使我们对转换波可以有进一步的认识,拓展了转换波处理中各向异性的应用.该方程对于偏移距深度比小于3.0的情况是比较准确的,这对于大偏移距转换波勘探具有实际意义.     

Moveout approximation for vertical seismic profile geometry in a 2D model with anisotropic layers

I introduce a new explicit form of vertical seismic profile (VSP) traveltime approximation for a 2D model with non‐horizontal boundaries and anisotropic layers. The goal of the new approximation is to dramatically decrease the cost of time calculations by reducing the number of calculated rays in a complex multi‐layered anisotropic model for VSP walkaway data with many sources. This traveltime approximation extends the generalized moveout approximation proposed by Fomel and Stovas. The new equation is designed for borehole seismic geometry where the receivers are placed in a well while the sources are on the surface. For this, the time‐offset function is presented as a sum of odd and even functions. Coefficients in this approximation are determined by calculating the traveltime and its first‐ and second‐order derivatives at five specific rays. Once these coefficients are determined, the traveltimes at other rays are calculated by this approximation. Testing this new approximation on a 2D anisotropic model with dipping boundaries shows its very high accuracy for offsets three times the reflector depths. The new approximation can be used for 2D anisotropic models with tilted symmetry axes for practical VSP geometry calculations. The new explicit approximation eliminates the need of massive ray tracing in a complicated velocity model for multi‐source VSP surveys. This method is designed not for NMO correction but for replacing conventional ray tracing for time calculations.     

CDP mapping in tilted transversely isotropic (TTI) media. Part I: Method and effectiveness

Transverse isotropy with a tilted axis of symmetry (TTI) causes image distortion if isotropic models are assumed during data processing. A simple anisotropic migration approach needs long computational times and is sensitive to the signal-to-noise ratio. This paper presents an efficient, general approach to common-depth-point (CDP) mapping to image the subsurface in TTI media from qP-wave seismic data by adding anisotropic and dip parameters to the velocity model. The method consists of three steps: (i) calculating traveltimes and positions of the CDP points; (ii) determining CDP trajectories; (iii) CDP imaging. A crucial step is the rapid computation of traveltimes and raypaths in the TTI media, which is achieved by the Fermat method, specially adapted for anisotropic layered media. The algorithm can image the subsurface of a given model quickly and accurately, and is suitable for application to a bending reflector. The effectiveness of the method is demonstrated by comparing the raypaths, the traveltimes and the results of CDP mapping, when assuming isotropic media, transversely isotropic media with a vertical axis of symmetry (TIV), and TTI media.     

3D description and inversion of reflection moveout of PS‐waves in anisotropic media

Common‐midpoint moveout of converted waves is generally asymmetric with respect to zero offset and cannot be described by the traveltime series t²(x²) conventionally used for pure modes. Here, we present concise parametric expressions for both common‐midpoint (CMP) and common‐conversion‐point (CCP) gathers of PS‐waves for arbitrary anisotropic, horizontally layered media above a plane dipping reflector. This analytic representation can be used to model 3D (multi‐azimuth) CMP gathers without time‐consuming two‐point ray tracing and to compute attributes of PS moveout such as the slope of the traveltime surface at zero offset and the coordinates of the moveout minimum. In addition to providing an efficient tool for forward modelling, our formalism helps to carry out joint inversion of P and PS data for transverse isotropy with a vertical symmetry axis (VTI media). If the medium above the reflector is laterally homogeneous, P‐wave reflection moveout cannot constrain the depth scale of the model needed for depth migration. Extending our previous results for a single VTI layer, we show that the interval vertical velocities of the P‐ and S‐waves (V_P0 and V_S0) and the Thomsen parameters ε and δ can be found from surface data alone by combining P‐wave moveout with the traveltimes of the converted PS(PSV)‐wave. If the data are acquired only on the dip line (i.e. in 2D), stable parameter estimation requires including the moveout of P‐ and PS‐waves from both a horizontal and a dipping interface. At the first stage of the velocity‐analysis procedure, we build an initial anisotropic model by applying a layer‐stripping algorithm to CMP moveout of P‐ and PS‐waves. To overcome the distorting influence of conversion‐point dispersal on CMP gathers, the interval VTI parameters are refined by collecting the PS data into CCP gathers and repeating the inversion. For 3D surveys with a sufficiently wide range of source–receiver azimuths, it is possible to estimate all four relevant parameters (V_P0, V_S0, ε and δ) using reflections from a single mildly dipping interface. In this case, the P‐wave NMO ellipse determined by 3D (azimuthal) velocity analysis is combined with azimuthally dependent traveltimes of the PS‐wave. On the whole, the joint inversion of P and PS data yields a VTI model suitable for depth migration of P‐waves, as well as processing (e.g. transformation to zero offset) of converted waves.     

三维复杂介质中转换波走时快速计算

复杂介质中转换波走时计算是多波勘探地震学中重要内容之一.本项研究利用惠更斯原理和费玛原理,获得了三维复杂介质中转换波快速计算的改进型最小走时树方法.其中,在保证精度的条件下,为了提高三维转换波走时计算效率,首先对初至波最小走时树基本算法进行了改进.本方法通过将转换波分为上、下行波分别进行射线追踪以实现三维转换波走时的快速计算.模型计算表明,方法的计算速度快,而且稳定性强,对多波地震勘探具有较大的应用价值.     

层状各向异性介质转换波克希霍夫叠前时间偏移

在克希霍夫叠前时间偏移处理中,地震波走时的计算方法是决定大偏移距地震资料成像品质的重要因素.在常规的三维转换波各向异性叠前时间偏移公式中,走时的计算是基于等效单层各向异性介质的非双曲线方法.用这种方法处理的成像道集,在偏移/深度比超过一定阈值后,成像道集中的反射同相轴将出现过偏现象,这种偏移不平的同相轴将影响偏移叠加的最佳响应,使得偏移成像波组呈低频化特征,最终降低三维转换波偏移成像质量.我们采用层状介质的走时计算方法代替常规算法,并且利用了常规方法的转换波各向异性偏移速度模型.基于层状介质的算法能够提高大偏移距转换波走时计算精度,克服中浅地层大偏移距远道成像道集中反射同相轴逐渐上翘的问题.两个地区的三维转换波资料处理结果证实,基于层状各向异性介质的转换波克希霍夫叠前时间偏移方法,明显改善了反射成像剖面的连续性和分辨率,提高成像剖面构造的可解释性.     

多震相层析成像与在介质中的多波型射线跟踪

本文介绍了多震相的层析成像的思路和算法，通过穿透和反射走时可以同时作出２维和３维慢度（速度的倒数）重建。我们分析了在穿透和反射数据中确定速度和深度的不确定性，并认识到深度扰动对反射走时异常比慢度扰动更敏感。由不同波类型所提供的对速度和深度的约束，这个算法实际上减少了在一般反射层析成像在速度和深度之间的不确定性，并且也避免了在穿透层析成像中的不确定问题。线性化反演是通过从反射界面深度由分离速度参数迭代进行的。使用一个快速的２－Ｄ和３－Ｄ射线跟踪算法来计算穿透和反射走时和对幔度及反射界面深度的偏导数。深度和速度都用立方Ｂ样条函数来进行参数化。合成例子表明，当同时考虑穿透和反射时间，层析成像的结果得到改进。这个方法也应用到英国煤炭测量局（ＢｒｉｔｉｓｈＣｏａｌＭｅａｓｕｒｅｓ）沿跨线排列所记录的逆ＶＳＰ数据组。通过使用波形配合技术，用同时确定时间延迟和叠加权，可以自动拾取旅行时间。所观察到的逆ＶＳＰ层析成像可比周围介质具有较低速度的两个断层区域成像。断层的位置由附近的反射测线所确定。本文还讨论了在复杂２－Ｄ和３－Ｄ非均匀各向同性介质中地震射线跟踪方法。界面的几何形状和水平速度场都通过使用非均匀步长立方Ｂ－样条节点     

Large-offset approximation to seismic reflection traveltimes

Conventional approximations of reflection traveltimes assume a small offset-to-depth ratio, and their accuracy decreases with increasing offset-to-depth ratio. Hence, they are not suitable for velocity analysis and stacking of long-offset reflection seismic data. Assuming that the offset is large, rather than small, we present a new traveltime approximation which is exact at infinite offset and has a decreasing accuracy with decreasing offset-to-depth ratio. This approximation has the form of a series containing powers of the offset from 1 to −∞. It is particularly accurate in the presence of a thin high-velocity layer above the reflector, i.e. in a situation where the accuracy of the Taner and Koehler series is poor. This new series can be used to gain insight into the velocity information contained in reflection traveltimes at large offsets, and possibly to improve velocity analysis and stacking of long-offset reflection seismic data.     

转换波方位各向异性裂缝检测技术研究及应用

HTI裂缝各向异性介质中,转换波随方位角的变化比较复杂,目前还没有解析公式可以表达其变化特征,只能通过物理实验或数值模拟来分析其应用的可能性.数值模拟结果表明,转换波在裂缝各向异性介质中传播时,其R分量和T分量的振幅属性都具有明显的方位各向异性特征,R分量振幅方位各向异性拟合椭圆的长轴方向指示裂缝方位,这与纵波方位各向异性特征相似;根据P波AVAZ方位各向异性分析原理,对转换波R分量振幅方位各向异性曲线进行方位椭圆拟合,寻找椭圆的长轴方向,即裂缝主方位,再由P波AVAZ技术中振幅响应与炮检方向和裂缝走向之间的夹角关系式得到裂缝的发育密度,从而构建转换波方位各向异性AVAZ裂缝检测技术.该技术已用于川西新场气田某区块的裂缝储层预测,取得了较好的应用效果.     

Background velocity estimation using non-linear optimization for reflection tomography and migration misfit

We show that it is possible to estimate the background velocity for prestack depth migration in 2D laterally varying media using a non-linear optimization technique called very fast simulated annealing (VFSA). We use cubic splines in the velocity model parametrization and make use of either successive pairs of shot gathers or several constant-offset sections as input data for the inversion. A Kirchhoff summation scheme based on first-arrival traveltimes is used to migrate/model the input data during the velocity analysis. We evaluate and compare two different measures of error. The first is defined in the recorded data or (x,t) domain and is based on a reflection-tomography criterion. The second is defined in the migrated data or (x,z) domain and is based on a migration-misfit criterion. Depth relaxation is used to improve the convergence and quality of the velocity analysis while simultaneously reducing the computational cost. Further, we show that by coarse sampling in the offset domain the method is still robust. Our non-linear optimization approach to migration velocity analysis is evaluated for both synthetic and real seismic data. For the velocity-analysis method based on the reflection-tomography criterion, traveltimes do not have to be picked. Similarly, the migration-misfit criterion does not require that depth images be manually compared. Interpreter intervention is required only to restrict the search space used in the velocity-analysis problem. Extension of the proposed schemes to 3D models is straightforward but practical only for the fastest available computers.     

AN ADVANCED CMP RAY-TRACING1

The ray-tracing algorithm presented in this paper is based on formulae derived for the common reflecting element (CRE) stacking method. A 2D, smooth, laterally-varying media is assumed where offset rays and traveltimes are evaluated from normal-incidence (central) rays. The method uses a second-order asymmetrical approximation for rays and an additional oblique spherical approximation of the central wavefronts for calculating offset traveltimes. In order to solve the two-point ray-tracing problem for the common midpoint (CMP) configuration of source-receiver pairs located symmetrically around the CMP stations, the central rays are perturbed to satisfy the above-mentioned asymmetrical distribution. Although the accuracy of the calculations is limited for far offsets, it is still good for distances of the order of the reflecting depths. Since only a few normal-incidence rays are traced through the medium, the method is very fast and is found to be most attractive for iterative inversions in macromodel estimation.     

Vertical seismic profile at Pike's Peak, Saskatchewan, Canada: turning rays and velocity anisotropy

First-arrival traveltimes from a multi-offset vertical seismic profile (VSP) were used to estimate velocity anisotropy in the presence of a vertical velocity gradient. A numerical model consisting of two layers with vertical velocity gradients of 3.1 and 1.2 s⁻¹, respectively, and global anisotropy parameters of =0.12±0.02 and δ=0.30±0.06 yielded first-arrival traveltimes that matched the observed traveltimes well. Shallow receivers were found to be crucial for constraining the vertical velocity field and for determining the parameters of anisotropy at depth.     

Finite-difference calculation of traveltimes based on rectangular grid

Introduction The calculation of seismic wave traveltimes is a basic and the most important step in tomo-graphy, seismic wave forward modeling and Kirchhoff prestack depth migration. Limitations withtraditional ray tracing fall into four categories. a) Analytical methods can only realize ray tracingfor simply varying velocity fields, so they have relative small applied-range; b) Shooting methodsof ray tracing can cause shadow zones. When the shadow zones exist the method will invalid; c)…     

Offset dependent amplitudes: Fault curvature effects and dip move-out corrections

Equations have been presented previously which predict that reflector curvature can affect significantly seismic reflection amplitudes at both zero and nonzero source-receiver offsets. Here the fact that faults are generally concave-upward is used to examine the curvature effect for compaction-driven faults which the sediments have both exponential and logarithmic porosity decreases with increasing depth.The curvature effect is generally larger over upwardly concave reflectors where amplitudes can either increase with offset (exposed focus) or decrease with offset (buried focus). The magnitude of the effect depends on the ratio between the depth to the structure and the radius of curvature of the structure. A phase change of 90° also occurs at a critical offset in the case of an exposed focus concave-upward reflector, with decreasing amplitude at offsets larger than the critical value. Dip move-out (DMO) dominantly removes the amplitude variation with offset due to curvature. These results suggest that when looking for amplitude variations with offset in a fault prospect, DMO should be applied as a preprocessing step. Compaction-driven faults have an exposed focus and, for an exponential variation of porosity with depth, there is a maximum curvature effect at a depth roughly the same as the scaling depth for the porosity. Logarithmic porosity with depth variations suggest a continued increase in the curvature effect of faults with increasing depth, to a maximum amplitude increase of about 35% relative to a plane interface.     

P‐wave stacking‐velocity tomography for VTI media

A major complication caused by anisotropy in velocity analysis and imaging is the uncertainty in estimating the vertical velocity and depth scale of the model from surface data. For laterally homogeneous VTI (transversely isotropic with a vertical symmetry axis) media above the target reflector, P‐wave moveout has to be combined with other information (e.g. borehole data or converted waves) to build velocity models for depth imaging. The presence of lateral heterogeneity in the overburden creates the dependence of P‐wave reflection data on all three relevant parameters (the vertical velocity V_P0 and the Thomsen coefficients ε and δ) and, therefore, may help to determine the depth scale of the velocity field. Here, we propose a tomographic algorithm designed to invert NMO ellipses (obtained from azimuthally varying stacking velocities) and zero‐offset traveltimes of P‐waves for the parameters of homogeneous VTI layers separated by either plane dipping or curved interfaces. For plane non‐intersecting layer boundaries, the interval parameters cannot be recovered from P‐wave moveout in a unique way. Nonetheless, if the reflectors have sufficiently different azimuths, a priori knowledge of any single interval parameter makes it possible to reconstruct the whole model in depth. For example, the parameter estimation becomes unique if the subsurface layer is known to be isotropic. In the case of 2D inversion on the dip line of co‐orientated reflectors, it is necessary to specify one parameter (e.g. the vertical velocity) per layer. Despite the higher complexity of models with curved interfaces, the increased angle coverage of reflected rays helps to resolve the trade‐offs between the medium parameters. Singular value decomposition (SVD) shows that in the presence of sufficient interface curvature all parameters needed for anisotropic depth processing can be obtained solely from conventional‐spread P‐wave moveout. By performing tests on noise‐contaminated data we demonstrate that the tomographic inversion procedure reconstructs both the interfaces and the VTI parameters with high accuracy. Both SVD analysis and moveout inversion are implemented using an efficient modelling technique based on the theory of NMO‐velocity surfaces generalized for wave propagation through curved interfaces.