Inversion of shear wave interval velocity on prestack converted wave data

Seismic velocity is important to migration of seismic data,interpretation of lithology and lithofacies as well asprediction of reservoir.The information of shear wave velocity is required to reduce the uncertainty for discrimi-nating lithology,identifying fluid type in porous material and calculating gas saturation in reservoir prediction.Based on Zoeppritz equations,a numeral and scanning method was proposed in this paper.Shear wave velocitycan be calculated with prestack converted wave data.The effects were demonstrated by inversion of theoreticaland real seismic data.     

The study and analysis of floating datum selection for rugged terrain

We analyze the characteristics of different floating datums for static corrections and discuss the methods for determining them. The effect of different floating datum corrections was studied using theoretical model experiments, resulting in the conclusion that the velocity obtained after the floating datum correction with the minimum static correction errors depends on the velocity of the layer below the low velocity layer （LVL） lower boundary and is not related to topographic relief and LVL structure. For the real data processing case, wave equation numerical model experiments were conducted which resulted in a new method for calculating objective functions based on the waveform and modifications to the calculation equation for minimum static correction errors to make the method suitable for real data static correction processing using inhomogeneous velocity models with lower velocity boundary relief. Real data processing results demonstrate the method＇s superiority.     

Propagation characteristics of converted refracted wave and its application in static correction of converted wave

Three-component seismic exploration through P-wave source and three-component geophone is an effective technique used in complicated reservoir exploration. In three-component seismic exploration data processing,one of the difficulties is static correction of converted wave. This paper analyzes propagation characteristics of non-converted and converted refracted waves,and discovers a favor-able condition for the formation of converted refracted wave,i.e. the velocity of overlaying medium S wave is much lower than that of underlying medium S wave. In addition,the paper proposes the static correction method of converted wave based on PPS converted refracted wave,and processes the real three-component seismic data with better results of static correction of converted wave.     

Migration velocity modeling based on common reflection surface gather

The common-reflection-surface (CRS) stacking is a new seismic imaging method, which only depends on seismic three parameters and near-surface velocity instead of macro-velocity model. According to optimized three parameters obtained by CRS stacking, we derived an analytical relationship between three parameters and migration velocity field, and put forward CRS gather migration velocity modeling method, which realize velocity estimation by optimizing three parameters in CRS gather. The test of a sag model proved that this method is more effective and adaptable for velocity modeling of a complex geological body, and the accuracy of velocity analysis depends on the precision of optimized three parameters.     

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.     

基于三角网格的有限差分法叠后逆时偏移

Compared with other migration methods, reverse-time migration is based on a precise wave equation, not an approximation, and performs extrapolation in the depth domain rather than the time domain. It is highly accurate and not affected by strong subsurface structure complexity and horizontal velocity variations. The difference method based on triangular grids maintains the simplicity of the difference method and the precision of the finite element method. It can be used directly for forward modeling on models with complex top surfaces and migration without statics preprocessing. We apply a finite difference method based on triangular grids for post-stack reverse-time migration for the first time. Tests on model data verify that the combination of the two methods can achieve near-perfect results in application.     

Converted-wave Seismology in Anisotropic Media Revisited, Part I： Basic Theory

We have developed new basic theories for calculating the conversion point and the travel time of the P-SV converted wave (C-wave) in anisotropic, inhomogeneous media. This enables the use of conventional procedures such as semblance analysis, Dix-type model building and Kirchhoff summation, to implement anisotropic processing, and makes anisotropic processing affordable. Here we present these new developments in two parts: basic theory and application to velocity analysis and parameter estimation. This part deals with the basic theory, including both conversion-point calculation and moveout analysis.Existing equations for calculating the PS-wave (C-wave) conversion point in layered media with vertical transverse isotropy (VTI) are strictly limited to offsets about half the reflector depth (an offset-depth ratio, x/z, of 0.5), and those for calculating the C-wave traveltimes are limited to offsets equal to the reflector depth (x/z=1.0). In contrast, the new equations for calculating the conversion-point extend into offsets about three-times the reflector depth (x/z=3.0), those for calculating the C-wave traveltimes extend into offsets twice the reflector depth (x/z=2.0). With the improved accuracy, the equations can help in C-wave data processing and parameter estimation in anisotropic, inhomogeneous media.     

Surface Wave Group Velocity Tomography Imaging from Ambient Noise for Fujian Province and Its Adjacent Areas

Two-month continuous waveforms of 108 broadband seismic stations in Fujian Province and its adjacent areas are used to compute noise cross-correlation function (NCF). The signal quality of NCF is improved via the application of time-frequency phase weighted stacking. The Rayleigh and Love waves group velocities between 1s-20s are measured on the symmetrical component of the NCF with the multiple filter method. More than 5,000 Rayleigh wave dispersion curves and about 4,000 Love wave dispersion curves are obtained and used to invert for group velocity maps. This data set provides about 50km resolution that is demonstrated with checkerboard tests. Considering the off great circle effect in inhomogeneous medium, the ray path is traced based on the travel time field computed with a finite difference method. The inverted group velocity maps show good correlation with the geological features in the upper and middle crust. The Fuzhou basin and Zhangzhou basin showed low velocity on the short period group velocity maps. On the long period group velocity maps, the low velocity anomaly in the high heat flow region near Zhangzhou and clear velocity contrast across the Zhenghe-Dapu faults, which suggests that the Zhenghe-Dapu fault might be a deep fault.     

The analysis and application of simplified two-parameter moveout equation for C-waves in VTI anisotropy media

Several parameters are needed to describe the converted-wave （C-wave） moveout in processing multi-component seismic data, because of asymmetric raypaths and anisotropy. As the number of parameters increases, the converted wave data processing and analysis becomes more complex. This paper develops a new moveout equation with two parameters for C-waves in vertical transverse isotropy （VTI） media. The two parameters are the C-wave stacking velocity （Vc2） and the squared velocity ratio （7v,i） between the horizontal P-wave velocity and C-wave stacking velocity. The new equation has fewer parameters, but retains the same applicability as previous ones. The applicability of the new equation and the accuracy of the parameter estimation are checked using model and real data. The form of the new equation is the same as that for layered isotropic media. The new equation can simplify the procedure for C-wave processing and parameter estimation in VTI media, and can be applied to real C-wave processing and interpretation. Accurate Vc2 and Yvti can be deduced from C-wave data alone using the double-scanning method, and the velocity ratio model is suitable for event matching between P- and C-wave data.     

用于AVO分析的振幅保真平面波叠前时间偏移

To support amplitude variation with offset （AVO） analysis in complex structure areas, we introduce an amplitude-preserving plane-wave prestack time migration approach based on the double-square-root wave equation in media with little lateral velocity variation. In its implementation, a data mapping algorithm is used to obtain offset-plane-wave data sets from the common-midpoint gathers followed by a non-recursive phase-shift solution with amplitude correction to generate common-image gathers in offset-ray-parameter domain and a structural image. Theoretical model tests and a real data example show that our prestack time migration approach is helpful for AVO analysis in complex geological environments.     

火山岩地震屏蔽层的转换波叠前时间偏移成像

在反射地震转换波资料处理中,准确求取共转换点一直是一个难题,采用叠前时间偏移技术能避免共转换点道集的抽取,而且能够使转换波归位到真正的反射点上,实现准确成像.本文针对火山岩地震屏蔽层的转换波成像问题,通过对转换波共近似转换点道集进行速度分析,建立了转换波叠前时间偏移的初始速度场,通过速度扫描和纵、横波速度比值扫描确定最佳的偏移速度场和纵、横波速度比值,实现了在火山岩高速层覆盖区域的转换波偏移成像.实际资料的成像结果表明,本文采用的近似转换点计算以及转换波叠前时间偏移方法是有效的.     

黏声介质平面波有限差分叠前深度偏移及在莺歌海盆地的应用

莺琼盆地诸探区中存在底劈现象,深层气源产生的气体沿底劈产生的裂隙通道向上漫溢.漫溢过程中,一是充填在遇到的砂体中形成气藏;二是弥漫在上溢通道中,使得通道中的纵波速度发生变化,进而纵波波阻抗差异变小,反射变弱.另外,通道中气体的存在,会加强地层的吸收衰减,使得地震波振幅变弱、高频成分损失导致同相轴分辨率降低.利用OBC数据进行多波地震勘探和利用黏声介质的叠前深度偏移都是改善模糊区成像质量的重要方法技术.为此,本文提出用黏声介质平面波有限差分法叠前深度偏移成像方法改善气体充填区域的成像质量.黏声介质成像目的是补偿地震波的吸收衰减;平面波偏移成像目的是适应海上单炮数量巨大,提高波动方程叠前深度偏移成像的效率;有限差分法叠前偏移的目的是适应该区浅层气分布局域性极强、Q值的空间变化大的情况.在莺歌海某探区的实际数据上的黏声介质平面波有限差分叠前深度偏移试验证明,本方法是改善模糊区成像质量的较为有效的途径.     

Velocity anisotropy characteristics and pre-stack time migration imaging of the marine sedimentary strata in the South Yellow Sea Basin

The South Yellow Sea is a superimposed basin overlying Mesozoic-Cenozoic continental sediments, which in turn overlie Paleozoic-Mesozoic marine deposits that are now the target of hydrocarbon exploration. Strongly modified by multiple tectonic events, the marine sediments feature a large tectonic relief, with obvious horizontal anisotropy in seismic velocity, which significantly affects the seismic image quality. In this study, the sedimentary velocity anisotropy and its influence on image quality were analyzed using an analytical theory method, assuming transversely isotropic medium with vertical axis of symmetry (VTI), and using seismic and well-log data. Additionally, an anisotropic prestack time migration was used for the imaging of the field data. The results showed that the anisotropic pre-stack time migration processing could be used to significantly improve the accuracy of the seismic images in areas with distinct faults, offering clear images of accurately located fault planes and fault edges, thereby improving the lateral resolution of the seismic data and its signal-to-noise ratio.     

两种类型共炮点域相关型叠前深度偏移方法分析

研究不同偏移方法的成像机理是实现复杂构造高精度成像的前提,研究了两类共炮点域的相关型叠前深度偏移成像方法：基于单程波动方程的叠前深度成像方法和基于双程波动方程的叠前深度成像方法,同时对比了它们在计算效率、数据存储量、成像精度、成像机理、速度敏感性等方面的差异及其共性。以复杂构造模型为例,采用了傅里叶有限差分法（FFD）和逆时成像法（RTM）,这两种方法实现了121个共炮点道集的叠前深度偏移成像处理。计算结果表明,当速度准确时,两种深度域波动方程成像方法均可以恢复出各个地质反射界面,其中逆时偏移对陡倾角成像效果显著,当速度存在百分比误差和随机扰动情况时,逆时成像结果要差于单程波方法,因此,逆时偏移方法对速度的敏感性较大,且低频噪声较为严重。     

Velocity model updating in prestack Kirchhoff time migration for PS converted waves: Part I – Theory

A velocity model updating approach is developed based on moveout analysis of the diffraction curve of PS converted waves in prestack Kirchhoff time migration. The diffraction curve can be expressed as a product of two factors: one factor depending on the PS converted‐wave velocity only, and the other factor depending on all parameters. The velocity‐dependent factor represents the hyperbolic behaviour of the moveout and the other is a scale factor that represents the non‐hyperbolic behaviour of the moveout. This non‐hyperbolic behaviour of the moveout can be corrected in prestack Kirchhoff time migration to form an inverse normal‐moveout common‐image‐point gather in which only the hyperbolic moveout is retained. This hyperbolic moveout is the moveout that would be obtained in an isotropic equivalent medium. A hyperbolic velocity is then estimated from this gather by applying hyperbolic moveout analysis. Theoretical analysis shows that for any given initial velocity, the estimated hyperbolic velocity converges by an iterative procedure to the optimal velocity if the velocity ratio is optimal or to a value closer to the optimal velocity if the velocity ratio is not optimal. The velocity ratio (V_P/V_S) has little effect on the estimation of the velocity. Applying this technique to a synthetic seismic data set confirms the theoretical findings. This work provides a practical method to obtain the velocity model for prestack Kirchhoff time migration.     

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

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

基于全波形反演的探地雷达数据逆时偏移成像

逆时偏移成像(RTM)常用来处理复杂速度模型,包括陡倾角及横向速度变化剧烈的模型.与常规偏移成像方法(如Kirchhoff偏移)相比,逆时偏移成像能提供更好的偏移成像结果,近些年逆时偏移成像越来越广泛地应用到勘探地震中,它逐渐成为石油地震勘探中的一种行业标准.电磁波和弹性波在动力学和运动学上存在相似性,故本文开发了基于麦克斯韦方程组的电磁波逆时偏移成像算法,并将其应用到探地雷达数据处理中.时间域有限差分(FDTD)用于模拟电磁波正向和逆向传播过程,互相关成像条件用于获得最终偏移结果.逆时偏移成像算法中,偏移成像结果受初始模型影响较大,而其中决定电磁波传播速度的介电常数的影响尤为重要.本文基于时间域全波形反演(FWI)算法反演获得了更为精确的地下介电常数模型,并将其反演结果作为逆时偏移成像的初始介电常数模型.为了验证此算法的有效性,首先构建了一个复杂地质结构模型,合成了共偏移距及共炮点探地雷达数据,分别应用常规Kirchhoff偏移算法及逆时偏移成像算法进行偏移处理,成像结果显示由逆时偏移成像算法得到的偏移结果与实际模型具有较高的一致性;此外本文在室内沙槽中进行了相关的物理模拟实验,采集了共偏移距及共炮点探地雷达数据,分别应用Kirchhoff和叠前逆时偏移成像算法进行处理,结果表明叠前逆时偏移成像在实际应用中能获得更好的成像效果.     

结合基准面重建的叠前时间偏移方法

提出了一种结合虚拟界面、瑞利积分和相移法的混合的基准面重建方法.通过与叠前时间偏移方法结合,形成了针对起伏地表采集数据的叠前时间偏移方法和新流程.该方法能正确考虑波在近地表传播的实际路径,克服了高速层出露时静校正方法的误差;它也能自己确定虚拟层速度,避免了现行基于波场延拓的基准面重建方法需要准确近地表速度的困难.文中分别用近地表存在明显低速层和近地表有高速层出露这两类模型的理论数据,验证了所发展方法和流程的有效性.     

适于Kirchhoff叠前深度偏移的地震走时李代数积分算法

本文基于拟微分算子理论和李代数积分法,根据程函方程和波场坐标变换,提出一种新的适于横向变速介质Kirchhoff叠前深度偏移的地震波走时算法.该算法与Kirchhoff叠前时间偏移所用李代数时间积分表达相比,差异在于增加了波数一次项,且二次项的系数在求积时亦需进行修正.针对单平方根算子象征、李代数积分、指数映射和走时多项式的求解而言,皆需对以往Kirchhoff叠前时间偏移中所用算法进行深化调整.文中数值算例对比了本文李代数积分表达与时间积分的区别,本算法计算结果与线性横向变速介质中的理论值相当吻合.通过走时多项式中各项对结果的影响分析,可知非对称项使计算精度得到了进一步提高.数值试验表明,本算法对横向变速介质中走时求取是可行的,且不需要存储海量走时表,有利于提高Kirchhof叠前深度偏移的精度和效率.     

Study and application of PS-wave pre-stack migration in HTI media and an anisotropic correction method

Anisotropy correction is necessary during the processing of converted PSwave seismic data to achieve accurate structural imaging, reservoir prediction, and fracture detection. To effectively eliminate the adverse effects of S-wave splitting and to improve PSwave imaging quality, we tested methods for pre-stack migration imaging and anisotropic correction of PS-wave data. We based this on the propagation rules of seismic waves in a horizontal transverse isotropy medium, which is a fractured medium model that reflects likely subsurface conditions in the field. We used the radial (R) and transverse (T) components of PS-wave data to separate the fast and slow S-wave components, after which their propagation moveout was effectively extracted. Meanwhile, corrections for the energies and propagation moveouts of the R and T components were implemented using mathematical rotation. The PS-wave imaging quality was distinctly improved, and we demonstrated the reliability of our methods through numerical simulations. Applying our methods to three-dimensional and three-component seismic field data from the Xinchang-Hexingchang region of the Western Sichuan Depression in China, we obtained high-quality seismic imaging with continuous reflection wave groups, distinct structural features, and specific stratigraphic contact relationships. This study provides an effective and reliable approach for data processing that will improve the exploration of complex, hidden lithologic gas reservoirs.