共转换点道集S变换谱相似性的多波地震数据横波波场提取技术

横波速度动校正后的共转换点(CCP)道集内,同时刻的各道横波信号S变换(ST)谱与其叠加道ST谱具有相似关系.因此,可基于这种相似关系设计自适应滤波器来提取多波地震数据中的横波波场.首先对共中心点(CMP)道集应用纵波速度动校正并在各道减去叠加道来去除数据中的纵波波场;然后在CCP道集应用横波速度动校正,将地震道振幅水平调整至叠加道振幅水平并做S变换,以叠加道ST谱为参考对地震道ST谱进行自适应滤波,去除数据中的残余纵波和噪声;最后,将滤波结果的振幅水平恢复至滤波前振幅水平.理论和实际数据试算表明,本文方法可有效提取多波地震数据中的横波波场,为多波多分量横波数据处理提供新思路.     

WEIGHTED STACKING FOR ROCK PROPERTY ESTIMATION AND DETECTION OF GAS*

Amplitude versus offset concepts can be used to generate weighted stacking schemes (here called geo-stack) which can be used in an otherwise standard seismic data processing sequence to display information about rock properties. The Zoeppritz equations can be simplified and several different approximations appear in the literature. They describe the variation of P-wave reflection coefficients with the angle of incidence of a P-wave as a function of the P-wave velocities, the S-wave velocities and the densities above and below an interface. Using a smooth, representative interval velocity model (from boreholes or velocity analyses) and assuming no dip, the angle of incidence can be found as a function of time and offset by iterative ray tracing. In particular, the angle of incidence can be computed for each sample in a normal moveout corrected CMP gather. The approximated Zoeppritz equation can then be fitted to the amplitudes of all the traces at each time sample of the gather, and certain rock properties can be estimated. The estimation of the rock properties is achieved by the application of time- and offset-variant weights to the data samples before stacking. The properties which can be displayed by geo-stack are: P-wave reflectivity (or true zero-offset reflectivity), S-wave reflectivity, and the reflectivity of P-wave velocity divided by S-wave velocity (or ‘pseudo-Poisson's ratio reflectivity’). If assumptions are made about the relation between P-wave velocity and S-wave velocity for water-bearing clastic silicate rocks, then it is possible to create a display which highlights the presence of gas.     

厚风化层覆盖区转换波静校正方法

P-SV转换波处理与传统的P-P波处理有很大的不同,如S波静校正、CCP叠加、P-SV速度分析和偏移等,其中最大的难题就是S波静校正问题.S波速度基本不受潜水面的影响,与纵波静校正没有直接相关性,有时横波静校正量能达到纵波静校正量的十倍,用纵波静校正量乘以比例系数来解决横波静校正问题将导致较大误差.同一接收点X和Z分量存在一定的初至时差,该时差代表了P波和S波在低降速带的走时差,可以利用该时差和近地表纵横波速度比信息去除低降速带对横波的影响,得到准确的静校正量.本文利用多分量初至时差推导了较为精确的横波静校正公式,再结合共检波点叠加求取剩余静校正量的方法,形成了完整的转换波静校正配套方法.利用该方法对苏里格气田二维及三维多波地震资料进行了实际处理,数据处理结果证明了该方法的有效性,该方法尤其适用于其他方法难以奏效的风化层较厚地区的横波静校正量求解,该方法也同时考虑了长波长横波静校正问题.     

天水-礼县地区地壳速度结构

1984—1985年,利用厂坝铅锌矿工业爆破,在天水—礼县地区布设测线进行了大范围的地震测深工作。对该地区的地壳速度结构的研究结果表明,该地区沉积层平均厚度为2.5km,速度为4.0km/s(P波);地壳平均厚度为43.68km,平均速度为6.20km/s;徽县—礼县地壳速度剖面可分为5层,其中在24—29km深处有一低速层,基底深度变化较大,在礼县地壳浅部发现一断层。对天水—礼县地区还进行了P波、S波联合反演,获得了该区P波与S波速度结构,其地壳范围内的平均波速比为1.73。     

The bound weighted average method (BWAM) for predicting S-wave velocity

The shear-wave velocity is a very important parameter in oil and gas seismic exploration, and vital in prestack elastic-parameters inversion and seismic attribute analysis. However, shearing-velocity logging is seldom carried out because it is expensive. This paper presents a simple method for predicting S-wave velocity which covers the basic factors that influence seismic wave propagation velocity in rocks. The elastic modulus of a rock is expressed here as a weighted arithmetic average between Voigt and Reuss bounds, where the weighting factor, w, is a measurement of the geometric details of the pore space and mineral grains. The S-wave velocity can be estimated from w, which is derived from the P-wave modulus. The method is applied to process well-logging data for a carbonate reservoir in Sichuan Basin, and shows the predicted S-wave velocities agree well with the measured S-wave velocities.     

Linearized elastic parameter sections1

Contrasts in elastic parameters can be estimated directly from seismic data using offset-dependent information in the PP reflection coefficient. A linear approximation to the PP reflection coefficient including three coefficients is fitted to the data, and relative contrasts in various elastic parameters are obtained from linear combinations of the estimated coefficients. Linearized elastic parameter sections for the contrasts in P-wave impedance, P-wave velocity, density, plane-wave modulus, and the change in bulk modulus and shear modulus normalized with the plane-wave modulus are estimated. If the average P- to S-wave velocity ratio is known, linearized parameter sections including the contrast in the average P- to S-wave velocity ratio and a fluid factor section can be computed. Applied to synthetic data, visual comparison of the estimated and true elastic parameter sections agree qualitatively, and the results are confirmed by an analysis of the standard deviation of the estimated parameters. The parameter sections obtained by inversion of a shallow seismic anomaly in the Barents Sea are promising, but the reliability is uncertain because neither well data nor regional trends are available.     

三维地震与OBS联合勘探揭示的神狐海域含水合物地层声波速度特征

以三维高分辨地震与海底高频地震仪（OBS）联合勘探数据为基础,获得海底之下沉积层的地震反射成像剖面及多波信息,并以此确定研究区含天然气水合物沉积层的纵、横波速度的变化特征.根据走时反演获得的横波速度与纵波速度对比分析发现,研究区海底之下500 m深度范围内的某些沉积层具有较高的纵横波速度,这一纵波速度升高区域与水合物稳定带对应,而纵波速度下降并且横波速度变化较小的区域,可能与游离气的存在相关.游离气的可能存在与基于这一区域2007年钻探测井结果的普遍认识不完全相符.     

CO2 saturation estimates at Sleipner (North Sea) from seismic tomography and rock physics inversion

CO₂ saturations are estimated at Sleipner using a two-step imaging workflow. The workflow combines seismic tomography (full-waveform inversion) and rock physics inversion and is applied to a two-dimensional seismic line located near the injection point at Sleipner. We use baseline data (1994 vintage, before CO₂ injection) and monitor data that was acquired after 12 years of CO₂ injection (2008 vintage). P-wave velocity models are generated using the Full waveform inversion technology and then, we invert selected rock physics parameters using an rock physics inversion methodology. Full waveform inversion provides high-resolution P-wave velocity models both for baseline and monitor data. The physical relations between rock physics properties and acoustic wave velocities in the Utsira unconsolidated sandstone (reservoir formation) are defined using a dynamic rock physics model based on well-known Biot–Gassmann theories. For data prior to injection, rock frame properties (porosity, bulk and shear dry moduli) are estimated using rock physics inversion that allows deriving physically consistent properties with related uncertainty. We show that the uncertainty related to limited input data (only P-wave velocity) is not an issue because the mean values of parameters are correct. These rock frame properties are then used as a priori constraint in the monitor case. For monitor data, the Full waveform inversion results show nicely resolved thin layers of CO₂–brine saturated sandstones under intra-reservoir shale layers. The CO₂ saturation estimation is carried out by plugging an effective fluid phase in the rock physics model. Calculating the effective fluid bulk modulus of the brine–CO₂ mixture (using Brie equation in our study) is shown to be the key factor to link P-wave velocity to CO₂ saturation. The inversion tests are done with several values of Brie/patchiness exponent and show that the CO₂ saturation estimates are varying between 0.30 and 0.90 depending on the rock physics model and the location in the reservoir. The uncertainty in CO₂ saturation estimation is usually lower than 0.20. When the patchiness exponent is considered as unknown, the inversion is less constrained and we end up with values of exponent varying between 5 and 20 and up to 33 in specific reservoir areas. These estimations tend to show that the CO₂–brine mixing is between uniform and patchy mixing and variable throughout the reservoir.     

Calculation of Rayleigh-wave phase velocities due to models with a high-velocity surface layer

Rayleigh-wave phase velocities have been utilized to determine shear (S)-wave velocities in near-surface geophysics since early 1980s. One of the key steps is to calculate theoretical dispersion curves of an earth model. When the S-wave velocity of the surface layer is higher than some of the layers below, however, the Rayleigh-wave phase velocity in a high-frequency range calculated by existing algorithms approaches the lowest S-wave velocity among the layers above the half-space, rather than a value related to the S-wave velocity of the surface layer. According to our numerical modeling results based on wave equation, trends of the Rayleigh-wave dispersive energy approach about a 91% of the S-wave velocity of the surface layer at a high-frequency range when its wavelength is much shorter than the thickness of the surface layer, which cannot be fitted by a dispersion curve calculated by existing algorithms. We propose a method to calculate Rayleigh-wave phase velocities of models with a high-velocity surface layer by considering its penetration depth. We build a substituted model that only contains the layer with the lowest S-wave velocity among the layers above the half-space and the layers above it. We use the substituted model to replace the original model to calculate phase velocities when the Rayleigh-wave wavelength is not long enough to penetrate the lowest S-wave velocity layer. Several synthetic models are used to verify fitness between the dispersion curve calculated by our proposed method and the trend of the highest dispersive energy. Examples of inversion also demonstrate high accuracy of using our method as the forward calculation method during the inversions.     

Joint inversion of high-frequency surface waves with fundamental and higher modes

Joint inversion of multimode surface waves for estimating the shear (S)-wave velocity has received much attention in recent years. In this paper, we first analyze sensitivity of phase velocities of multimodes of surface waves for a six-layer earth model, and then we invert surface-wave dispersion curves of the theoretical model and a real-world example. Sensitivity analysis shows that fundamental mode data are more sensitive to the S-wave velocities of shallow layers and are concentrated on a very narrow frequency band, while higher mode data are more sensitive to the parameters of relatively deeper layers and are distributed over a wider frequency band. These properties provide a foundation of using a multimode joint inversion to define S-wave velocities. Inversion results of both synthetic data and a real-world example demonstrate that joint inversion with the damped least-square method and the singular-value decomposition technique to invert high-frequency surface waves with fundamental and higher mode data simultaneously can effectively reduce the ambiguity and improve the accuracy of S-wave velocities.     

Joint AVO inversion, wavelet estimation and noise-level estimation using a spatially coupled hierarchical Bayesian model

The main objective of the AVO inversion is to obtain posterior distributions for P-wave velocity, S-wave velocity and density from specified prior distributions, seismic data and well-log data. The inversion problem also involves estimation of a seismic wavelet and the seismic-noise level. The noise model is represented by a zero mean Gaussian distribution specified by a covariance matrix. A method for joint AVO inversion, wavelet estimation and estimation of the noise level is developed in a Bayesian framework. The stochastic model includes uncertainty of both the elastic parameters, the wavelet, and the seismic and well-log data. The posterior distribution is explored by Markov-chain Monte-Carlo simulation using the Gibbs' sampler algorithm. The inversion algorithm has been tested on a seismic line from the Heidrun Field with two wells located on the line. The use of a coloured seismic-noise model resulted in about 10% lower uncertainties for the P-wave velocity, S-wave velocity and density compared with a white-noise model. The uncertainty of the estimated wavelet is low. In the Heidrun example, the effect of including uncertainty of the wavelet and the noise level was marginal with respect to the AVO inversion results.     

北京地区地震参数与速度结构的联合测定

为提高北京遥测地震台网的地震参数测定能力,对北京地区（38°-41°N,114°-119°E）建立了一个新的地震波速度模型MDBJ81。这一模型有四个水平层置于均匀半空间之上,第一层分成三块,以反映该地区浅部结构的横向差异及地形起伏效应。在联合测定地震参数和速度结构过程中,其可调整速度模型参数包括:各层中的平均P波速度,各层层顶深度,第一层三块厚度以及平均P波S波速度比。通过分离参数,避免了求解大型方程组。 最后求解中使用了北京台网1979年记录的43个地震,共836个P波及S波震相到时。结果表明,北京周围地区的确存在速度结构的横向差异,第一层底深度自西北向东南逐块增加。得到的新震中位置一般比原台网测定位置更靠近表面断层线。     

Inversion of Scholte wave dispersion and waveform modeling for shallow structure of the Ninetyeast Ridge

The construction of S-wave velocity models of marine sediments down to hundreds of meters below the seafloor is important in a number of disciplines. One of the most significant trends in marine geophysics is to use interface waves to estimate shallow shear velocities which play an important role in determining the shallow crustal structure. In marine settings, the waves trapped near the fluid–solid interface are called Scholte waves, and this is the subject of the study. In 1998, there were experiments on the Ninetyeast Ridge (Central Indian Ocean) to study the shallow seismic structure at the drilled site. The data were acquired by both ocean bottom seismometer and ocean bottom hydrophone. A new type of seafloor implosion sources has been used in this experiment, which successfully excited fast and high frequency (>500 Hz) body waves and slow, intermediate frequency (<20 Hz) Scholte waves. The fundamental and first higher mode Scholte waves have both been excited by the implosion source. Here, the Scholte waves are investigated with a full waveform modeling and a group velocity inversion approach. Shear wave velocities for the uppermost layers of the region are inferred and results from the different methods are compared. We find that the full waveform modeling is important to understand the intrinsic attenuation of the Scholte waves between 1 and 20 Hz. The modeling shows that the S-wave velocity varies from 195 to 350 m/s in the first 16 m of the uppermost layer. Depths levels of high S-wave impedance contrasts compare well to the layer depth derived from a P-wave analysis as well as from drilling data. As expected, the P- to S-wave velocity ratio is very high in the uppermost 16 m of the seafloor and the Poisson ratio is nearly 0.5. Depth levels of high S-wave impedance contrasts are comparable to the layer depth derived from drilling data.     

A physical model for shear-wave velocity prediction1

The clay-sand mixture model of Xu and White is shown to simulate observed relationships between S-wave velocity (or transit time), porosity and clay content. In general, neither S-wave velocity nor S-wave transit time is a linear function of porosity and clay content. For practical purposes, clay content is approximated by shale volume in well-log applications. In principle, the model can predict S-wave velocity from lithology and any pair of P-wave velocity, porosity and shale volume. Although the predictions should be the same if all measurements are error free, comparison of predictions with laboratory and logging measurements show that predictions using P-wave velocity are the most reliable. The robust relationship between S- and P-wave velocities is due to the fact that both are similarly affected by porosity, clay content and lithology. Moreover, errors in the measured P-wave velocity are normally smaller than those in porosity and shale volume, both of which are subject to errors introduced by imperfect models and imperfect parameters when estimated from logs. Because the model evaluates the bulk and shear moduli of the dry rock frame by a combination of Kuster and Toksöz’ theory and differential effective medium theory, using pore aspect ratios to characterize the compliances of the sand and clay components, the relationship between P- and S-wave velocities is explicit and consistent. Consequently the model sidesteps problems and assumptions that arise from the lack of knowledge of these moduli when applying Gassmann's theory to this relationship, making it a very flexible tool for investigating how the v_P-v_s relationship is affected by lithology, porosity, clay content and water saturation. Numerical results from the model are confirmed by laboratory and logging data and demonstrate, for example, how the presence of gas has a more pronounced effect on P-wave velocity in shaly sands than in less compliant cleaner sandstones.     

Cauchy prior distribution-based AVO elastic parameter estimation via weakly nonlinear waveform inversion

Cauchy priori distribution-based Bayesian AVO reflectivity inversion may lead to sparse estimates that are sensitive to large reflectivities. For the inversion, the computation of the covariance matrix and regularized terms requires prior estimation of model parameters, which makes the iterative inversion weakly nonlinear. At the same time, the relations among the model parameters are assumed linear. Furthermore, the reflectivities, the results of the inversion, or the elastic parameters with cumulative error recovered by integrating reflectivities are not well suited for detecting hydrocarbons and fuids. In contrast, in Bayesian linear AVO inversion, the elastic parameters can be directly extracted from prestack seismic data without linear assumptions for the model parameters. Considering the advantages of the abovementioned methods, the Bayesian AVO reflectivity inversion process is modified and Cauchy distribution is explored as a prior probability distribution and the time-variant covariance is also considered. Finally, we propose a new method for the weakly nonlinear AVO waveform inversion. Furthermore, the linear assumptions are abandoned and elastic parameters, such as P-wave velocity, S-wave velocity, and density, can be directly recovered from seismic data especially for interfaces with large reflectivities. Numerical analysis demonstrates that all the elastic parameters can be estimated from prestack seismic data even when the signal-to-noise ratio of the seismic data is low.     

Constraining the anisotropy structure of the crust by joint inversion of seismic reflection travel times and wave polarizations

In the context of wide-angle seismic profiling, the determination of the physical properties of the Earth crust, such as the elastic layer depth and seismic velocity, is often performed by inversion of P- and/or S-phases propagation data supplying the geometry of the medium (reflector depths) or any other structural parameter (P- or S-wave velocity, density...). Moreover, the inversion for velocity structure and interfaces is commonly performed using only seismic reflection travel times and/or crustal phase amplitudes in isotropic media. But it is very important to utilize more available information to constrain the non-uniqueness of the solution. In this paper, we present a simultaneous inversion method of seismic reflection travel times and polarizations data of transient elastic waves in stratified media to reconstruct not only layer depth and vertical P-wave velocity but also the anisotropy feature of the crust based on the estimation of the Thomsen’s parameters. We carry out a checking with synthetic data, comparing the inversion results obtained by anisotropic travel-time inversion to the results derived by joint inversion of seismic reflection travel times and polarizations data. The comparison proves that the first procedure leads to biased anisotropic models, while the second one fits nearly the real model. This makes the joint inversion method feasible. Finally, we investigate the geometry, P-wave velocity structure and anisotropy of the crust beneath Southeastern China by applying the proposed inversion method to previously acquired wide-angle seismic data. In this case, the anisotropy signature provides clear evidence that the Jiangshan-Shaoxing fault is the natural boundary between the Yangtze and Cathaysia blocks.     

基阶与高阶瑞利波联合反演研究

研究了六层层状介质模型瑞利波基阶和高模式波相速度对横波速度、深度的敏感性,结果表明：基阶波较高模式波对7 m以内浅部地层的横波速度更敏感,敏感性频带在10～25 Hz范围内,峰值频带集中在18 Hz左右;高模式波较基阶波对深部地层的横波速度更敏感,敏感性频带宽,峰值分散.基阶波对浅层的敏感性和高模式波穿透深度更深的特点为近地表岩土层二维横波速度结构的联合反演提供了前提条件.利用阻尼最小二乘SVD（Singular Value Decomposition）算法联合基阶与高模式波对理论模型和实例数据进行横波速度反演,反演结果表明联合反演增强了反演的稳定性,提高了反演的精度.     

Seismic inversion for underground fractures detection based on effective anisotropy and fluid substitution

Underground fractures play an important role in the storage and movement of hydrocarbon fluid. Fracture rock physics has been the useful bridge between fracture parameters and seismic response. In this paper, we aim to use seismic data to predict subsurface fractures based on rock physics. We begin with the construction of fracture rock physics model. Using the model, we may estimate P-wave velocity, S-wave velocity and fracture rock physics parameters. Then we derive a new approximate formula for the analysis of the relationship between fracture rock physics parameters and seismic response, and we also propose the method which uses seismic data to invert the elastic and rock physics parameters of fractured rock. We end with the method verification, which includes using well-logging data to confirm the reliability of fracture rock physics effective model and utilizing real seismic data to validate the applicability of the inversion method. Tests show that the fracture rock physics effective model may be used to estimate velocities and fracture rock physics parameters reliably, and the inversion method is resultful even when the seismic data is added with random noise. Real data test also indicates the inversion method can be applied into the estimation of the elastic and fracture weaknesses parameters in the target area.     

纵横波弹性阻抗联立反演在GD地区的应用

在GD油田复杂油藏描述中应用叠前纵横波弹性阻抗反演,精确地进行了油藏岩性的划分。利用三个或三个以上部分叠加数据,进行纵横波弹性阻抗联立反演,既克服了因叠后地震反演结果单一而不能满足复杂储层描述的需求,又避免了由于叠前道集信噪比低造成反演结果不稳定的缺陷。本文论述了叠前弹性波阻抗反演的基本原理,结合GD地区实际资料,对反演过程中涉及的角道集子波提取、层位标定、横波速度预测、弹性参数提取与解释等关键步骤进行了详细研究,指出基于测井资料分析的多种弹性参数综合解释是提高叠前地震反演应用效果的关键。