共查询到20条相似文献，搜索用时 281 毫秒
1.
为了实现三维倾斜界面真正的共反射点叠加,获得更加精确的P波速度模型,进行可靠的AVO/AVA分析,形成与三维多分量PS转换波共转换点道集配套的技术,在常速介质共反射点精确走时的基础上,推导并给出三维倾斜界面P波共反射点道集近似双曲时距关系及NMO速度;认识到近似时距在偏移距与深度比小于2时具有较高的拟合精度;指出P波共反射点道集NMO速度是P波速度、界面倾角、界面倾向和测线方位的函数,且随测线方位变化而呈现椭圆特征.与P波共中心点道集NMO速度不同.P波共反射点道集NMO速度小于介质速度.且随界面倾角增大而减小;椭圆长轴方向为倾斜界面走向方位,短轴方向为倾斜界面倾向方位. 相似文献
2.
三维三分量(3D3C)陆地反射PS转换波共中心点(CMP)叠加成像方法,虽然抽道集简单,但是对实际资料处理结果往往不理想.尤其当反射界面为三维倾斜界面时,其成像质量较差.本文提出有三个主要因素影响其成像质量:第一,转换点离散.运用实例计算得出,转换点离散度随着纵横波速度比、偏移距和界面倾角的增大而增大.相同界面倾角,不同测线方位的转换点离散度不同,视倾角的绝对值越大离散度也越大;第二,道集内静校正量差异增大.CMP道集中,由于转换点离散使得转换点横向跨度较大,经倾斜界面反射转换的S波出射到近地表地层时的角度差异也较大,导致静校突出;第三,加大动校叠加复杂性.三维倾斜界面PS波CMP道集近炮检距时距方程可表示为双曲形式,但是曲线的顶点位置和动校速度同时随测线方位变化,使得CMP道集同相轴很难校平,动校叠加过程很复杂. 相似文献
3.
对于三维反射地震,文中将三维弯曲界面视为倾向和倾角变化的三维倾斜界面的包络,并依据三维倾斜界面反射提出计算三维弯曲界面P波反射时距的一种新方法。该方法基于倾角CDP理论,从界面反射点出发,给出地表CDP网格点,由检波点位置计算炮点位置,得到三维地震的炮点网格和检波点网格的分布关系,最终计算P波的三维反射时距。炮集数值模拟的结果表明,该方法有较高的计算效率,便于分析反射时距随炮点位置的变化,特别能揭示界面反射点位置和时距的关系。三维弯曲界面和三维倾斜界面的反射时距特征对比表明,仅从同一测线方位的时距难以将这两种类型的界面反射区分开来,但比较不同方位的时距差异可以解决这一问题。 相似文献
4.
《地球物理学进展》2016,(1)
近年来,由于勘探区域构造的逐渐复杂化,要求我们对复杂界面反射的时距特征有更进一步认识.本文在已有对复杂界面研究的基础上,结合实际3D复杂地质构造特征,对复杂条件做了两个限定:一是反射界面有界,二是复杂界面有多个不同性质界面组合而成.在这两条限定条件下,从界面反射点出发,进行炮集P波反射时距特征研究.结果表明,复杂界面反射会生产多条时距曲线,反射时距特征随测线方位的变化而变化,不同方位的时距对应的三维反射路径、反射点范围、反射点位置以及界面特征不同.其次时距曲线之间存在着比较复杂的位置关系,包括交叉、断续和无信号道等特征,这些复杂特征与组合模型中倾斜界面的倾向、倾角、模型范围和深度有关.最后对复杂模型的理论数值分析得出了反射点范围与接收点范围的对应关系,并理论证明了复杂条件下部分反射特征以及反射盲区的存在,这与实际三维地震勘探相符合.此研究有助于实际生产的观测设计及数据处理解释,对于提高三维复杂构造反射界面的认识具有理论意义. 相似文献
5.
系统研究了台站下方存在倾斜界面时,远震体波接收函数的波场特征及多方位接收函数对台站下方介质横向非均匀特征的反映能力．我们利用合成三维横向非均匀介质接收函数的Maslov方法,具体模拟了台站下方存在倾斜界面时的接收函数响应．结果表明,当台站下方存在倾向一致的倾斜界面时,远震体波接收函数的径向分量和切向分量分别随震源方位角的变化呈现对称和反对称的规律性变化．利用不同震源方位角接收函数径向分量和切向分量的变化规律,可以估计台站下方界面的倾向和倾角．当台站下方各界面的倾向不同时,随方位角的变化,接收函数只能直观给出界面的整体倾向．实际观测数据的分析结果表明,对于利用单个台站接收函数研究台站下方介质的横向非均匀特征来说,简单倾斜界面可以看作较好的一级近似． 相似文献
6.
自然界的裂隙排列并非都是直立的,在复杂地壳应力作用下裂隙排列很可能是倾斜的.本文从各向异性弹性动力学基本方程出发,在特征矩阵方法基础上提出了计算任意倾斜裂隙介质水平界面上的反射、透射系数的简单方法,并利用该方法计算了倾斜裂隙介质P-P波、P-SV波和P-SH波的反透射系数.结果表明,P-P波、P-SV波和P-SH波的反透射系数与倾斜裂隙的倾角、测线方位以及偏移距这些物理参数之间存在一定的依赖关系.此外我们计算了C波的零入射相角时的反射系数,发现该反射系数跟裂隙倾角、裂隙的方位角有着明显的对应规律,且各向异性越强时,该反射系数越大.这些研究对于深入了解倾斜裂隙的动力学特征具有一定的意义. 相似文献
7.
非双曲（远偏移距）时差为各向异性介质中正、反演研究,特别是各向异性参数估计提供了重要信息.本文在任意空间取向TI（ATI）介质水平界面同类反射波四次时差系数(A4)精确解的基础上,进一步讨论我们推导得出的ATI介质中四次时差系数解析近似解,比较随CMP测线方位变化的近似解与精确解之间的差别,为利用近似解来解析研究ATI介质中非双曲时距以及参数反演提供有价值的信息．结合实际岩性资料的数值研究表明,ATI条件下四次时差系数近似解与精确解之间存在差别,不仅表现在A4数的大小及符号特征上,更突出地表现在A4系数随方位的变化特征上;在强各向异性条件下,近似解相比精确解存在较大误差．但在各向异性参数满足0<ε-δ<0.15、｜δ｜<0.20的情况下,对于TI对称轴的特殊倾角范围(75°～80°),近似解与精确解的差别很小,可用近似解进行各向异性观测解释及参数反演． 相似文献
8.
本文利用各向异性反射率技术计算理论地震图,提出海底高速薄层会产生沿高速层水平传播的波（简称径向波）,这种波在水层中作为P波,在固液界面激发下行横波,该均匀横波以临界角入射高速薄层,在层内作为超临界角的非均匀横波水平传播,再以临界角转换为上行传播的均匀横波,最终在固液界面上行透射转换为水层中P波.高速薄层传播的径向波不同于界面折射波,也不同于具有频散的面波和通道波.理论地震图的研究表明,径向波具有线性时距,能与海底强反射具有同等振幅水平;径向波有其振幅、时距位置和斜率这些观测记录参数,分别对应高速层的厚度、深度和近似的横波速度;径向波可以克服折射波解释中遇到的振幅强弱和高速层速度等困难.径向波可作为探测海底高速薄层的有力工具,对于研究高速层屏蔽、海底反射类型的多样性和相应的资料处理解释有重要意义. 相似文献
9.
相速度和偏振方向是研究地震波传播规律和描述介质特性的重要参数,在理论研究和实际应用中有重要作用.本文假定倾斜横向各向同性(TTI)介质对称轴位于观测坐标系XOZ面内,在此观测坐标系下直接推导了TTI介质弹性波相速度和偏振方向的解析表达式,再进一步利用Thomsen弱各向异性理论,推导了弱各向异性近似条件下弹性波相速度以及qP波和qSV波偏振方向表达式.理论分析和数值试例表明,在相速度方面,随着各向异性介质参数改变,qP波和qSH波速度变化较为平缓,qSV波速度变化较为剧烈.弹性波相速度近似式误差均较小,能较好地近似精确相速度.在偏振方向方面,SH波偏振方向只是传播方向和对称轴倾角的函数,而与各向异性参数无关,SH波偏振方向既垂直于传播方向,又垂直于TTI介质对称轴方向.除特定方向外,qP波和qSV波的偏振方向与传播方向均成一定角度,并且随TTI介质对称轴倾角的改变而改变;在精确和近似情况下,qP波和qSV波的偏振方向始终垂直;在精度允许范围内,偏振方向的弱各向异性近似式与理论解析式吻合较好. 相似文献
10.
11.
Stacking velocity V C2, vertical velocity ratio γ 0, effective velocity ratio γ eff, and anisotropic parameter χ eff are correlated in the PS-converted-wave (PS-wave) anisotropic prestack Kirchhoff time migration (PKTM) velocity model and are thus difficult to independently determine. We extended the simplified two-parameter (stacking velocity V C2 and anisotropic parameter k eff) moveout equation from stacking velocity analysis to PKTM velocity model updating and formed a new four-parameter (stacking velocity V C2, vertical velocity ratio γ 0, effective velocity ratio γ eff, and anisotropic parameter k eff) PS-wave anisotropic PKTM velocity model updating and process flow based on the simplified two-parameter moveout equation. In the proposed method, first, the PS-wave two-parameter stacking velocity is analyzed to obtain the anisotropic PKTM initial velocity and anisotropic parameters; then, the velocity and anisotropic parameters are corrected by analyzing the residual moveout on common imaging point gathers after prestack time migration. The vertical velocity ratio γ 0 of the prestack time migration velocity model is obtained with an appropriate method utilizing the P- and PS-wave stacked sections after level calibration. The initial effective velocity ratio γ eff is calculated using the Thomsen (1999) equation in combination with the P-wave velocity analysis; ultimately, the final velocity model of the effective velocity ratio γ eff is obtained by percentage scanning migration. This method simplifies the PS-wave parameter estimation in high-quality imaging, reduces the uncertainty of multiparameter estimations, and obtains good imaging results in practice. 相似文献
12.
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. 相似文献
13.
Transversely isotropic models with a tilted symmetry axis have become standard for imaging beneath dipping shale formations and in active tectonic areas. Here, we develop a methodology of wave-equation-based image-domain tomography for acoustic tilted transversely isotropic media. We obtain the gradients of the objective function using an integral wave-equation operator based on a separable dispersion relation that takes the symmetry-axis tilt into account. In contrast to the more conventional differential solutions, the integral operator produces only the P-wavefield without shear-wave artefacts, which facilitates both imaging and velocity analysis. The model is parameterized by the P-wave zero-dip normal-moveout velocity, the Thomsen parameter δ, anellipticity coefficient η and the symmetry-axis tilt θ. Assuming that the symmetry axis is orthogonal to reflectors, we study the influence of parameter errors on energy focusing in extended (space-lag) common-image gathers. Distortions in the anellipticity coefficient η introduce weak linear defocusing regardless of reflector dip, whereas δ influences both the energy focusing and depth scale of the migrated section. These results, which are consistent with the properties of the P-wave time-domain reflection moveout in tilted transversely isotropic media, provide important insights for implementation of velocity model-building in the image-domain. Then the algorithm is tested on a modified anticline section of the BP 2007 benchmark model. 相似文献
14.
It has been shown in the past that the interval-NMO velocity and the non-ellipticity parameter largely control the P-wave reflection time moveout of VTI media. To invert for these two parameters, one needs either reasonably large offsets, or some structure in the subsurface in combination with relatively mild lateral velocity variation.This paper deals with a simulation of an inversion approach, building on the assumption that accurately measured V
NMO, as defined by small offset asymptotics for a particular reflector, were available. Instead of such measurements we take synthetically computed data. First, an isotropic model is constructed which explains these V
NMO. Subsequently, residual moveout in common image gathers is modelled by ray tracing (replacing real data), along with its sensitivity for changes in the interval-NMO velocity and the non-ellipticity parameter under the constraint that V
NMO is preserved. This enables iterative updating of the non-ellipticity parameter and the interval-NMO velocity in a layer that can be laterally inhomogeneous.This approach is successfully applied for a mildly dipping reflector at the bottom of a layer with laterally varying medium parameters. With the exact V
NMO assumed to be given, lateral inhomogeneity and anisotropy can be distinguished for such a situation. However, for another example with a homogeneous VTI layer overlying a curved reflector with dip up to 30°, there appears to be an ambiguity which can be understood by theoretical analysis. Consistently with existing theory using the NMO-ellipse, the presented approach is successfully applied to the latter example if V
NMO in the strike direction is combined with residual moveout in dip direction. 相似文献
15.
16.
Despite the complexity of wave propagation in anisotropic media, reflection moveout on conventional common-midpoint (CMP) spreads is usually well described by the normal-moveout (NMO) velocity defined in the zero-offset limit. In their recent work, Grechka and Tsvankin showed that the azimuthal variation of NMO velocity around a fixed CMP location generally has an elliptical form (i.e. plotting the NMO velocity in each azimuthal direction produces an ellipse) and is determined by the spatial derivatives of the slowness vector evaluated at the CMP location. This formalism is used here to develop exact solutions for the NMO velocity in anisotropic media of arbitrary symmetry. For the model of a single homogeneous layer above a dipping reflector, we obtain an explicit NMO expression valid for all pure modes and any orientation of the CMP line with respect to the reflector strike. The contribution of anisotropy to NMO velocity is contained in the slowness components of the zero-offset ray (along with the derivatives of the vertical slowness with respect to the horizontal slownesses) — quantities that can be found in a straightforward way from the Christoffel equation. If the medium above a dipping reflector is horizontally stratified, the effective NMO velocity is determined through a Dix-type average of the matrices responsible for the ‘interval’ NMO ellipses in the individual layers. This generalized Dix equation provides an analytic basis for moveout inversion in vertically inhomogeneous, arbitrarily anisotropic media. For models with a throughgoing vertical symmetry plane (i.e. if the dip plane of the reflector coincides with a symmetry plane of the overburden), the semi-axes of the NMO ellipse are found by the more conventional rms averaging of the interval NMO velocities in the dip and strike directions. Modelling of normal moveout in general heterogeneous anisotropic media requires dynamic ray tracing of only one (zero-offset) ray. Remarkably, the expressions for geometrical spreading along the zero-offset ray contain all the components necessary to build the NMO ellipse. This method is orders of magnitude faster than multi-azimuth, multi-offset ray tracing and, therefore, can be used efficiently in traveltime inversion and in devising fast dip-moveout (DMO) processing algorithms for anisotropic media. This technique becomes especially efficient if the model consists of homogeneous layers or blocks separated by smooth interfaces. The high accuracy of our NMO expressions is illustrated by comparison with ray-traced reflection traveltimes in piecewise-homogeneous, azimuthally anisotropic models. We also apply the generalized Dix equation to field data collected over a fractured reservoir and show that P-wave moveout can be used to find the depth-dependent fracture orientation and to evaluate the magnitude of azimuthal anisotropy. 相似文献
17.
Common‐midpoint moveout of converted waves is generally asymmetric with respect to zero offset and cannot be described by the traveltime series t2(x2) 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 (VP0 and VS0) 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 (VP0, VS0, ε 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. 相似文献
18.
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. 相似文献
19.
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 VP0 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. 相似文献
20.
Proper stacking of three-dimensional seismic CDP-data generally requires the knowledge of normal moveout velocities in all source-receiver directions contributing to a CDP-gather. The azimuthal variation of the stacking velocities mainly depends on the dip of the seismic interfaces. For a single dipping plane a simple relation exists between the dip and the azimuthal variation of NMO-velocity. Varying strike and dip of subsequent reflectors, however, result in a complex dependency of the seismic parameters. Reliable information on the spatial distribution of the normal moveout (NMO)-velocity can be derived from a wavefront curvature estimation using a 3-D ray-tracing technique. These procedures require additional information, e.g. reflection time gradients or depth maps to show interval velocities between leading interfaces. Moreover, their application to an extended 3-D data volume is restricted by high costs. The need for a routine 3-D procedure resulted in a special data selection to create pseudo 2-D profiles and to apply existing velocity estimation routines to these profiles. At least three estimates in different directions are necessary to derive the full azimuthal velocity variation, characterized by the large and the small main axis and the orientation of the velocity ellipse. Errors are estimated by means of computer models. Stacking velocities obtained by mathematical routines (least-squares fit) and by seismic standard routines (NMO-correction and correlation) are compared. Finally, a general 3-D velocity procedure using cross-correlation of preliminarily NMO-corrected traces is proposed. 相似文献