共查询到20条相似文献，搜索用时 343 毫秒
1.
Anisotropic reverse-time migration for tilted TI media 总被引：1，自引：0，他引：1
Seismic anisotropy in dipping shales results in imaging and positioning problems for underlying structures. We develop an anisotropic reverse‐time depth migration approach for P‐wave and SV‐wave seismic data in transversely isotropic (TI) media with a tilted axis of symmetry normal to bedding. Based on an accurate phase velocity formula and dispersion relationships for weak anisotropy, we derive the wave equation for P‐wave and SV‐wave propagation in tilted transversely isotropic (TTI) media. The accuracy of the P‐wave equation and the SV‐wave equation is analyzed and compared with other acoustic wave equations for TTI media. Using this analysis and the pseudo‐spectral method, we apply reverse‐time migration to numerical and physical‐model data. According to the comparison between the isotropic and anisotropic migration results, the anisotropic reverse‐time depth migration offers significant improvements in positioning and reflector continuity over those obtained using isotropic algorithms. 相似文献
2.
地球介质相对于地震波波长尺度的定向非均匀性会导致波速的各向异性,进而影响地震波场的运动学与动力学特征．各向异性弹性波动方程是描述该类介质波场传播的基本工具,在正演模拟、偏移成像与参数反演中起着关键作用．为了面向实际应用构建灵活、简便的各向异性波场传播算子,人们一直在寻求简化的各向异性波动方程．本文借鉴各向异性弹性波波型分离思想,通过对平面波形式的弹性波方程(即Christoffel方程)实施一种代表向波矢量方向投影的相似变换,推导出了一种适应任意各向异性介质、运动学上与原始弹性波方程完全等价,在动力学上突出qP波的新方程,即qP波伪纯模式波动方程．文中以横向各向同性(TI)介质为例,给出了相应的qP波伪纯模式波动方程及其声学与各向同性近似,并在此基础上开展了正演模拟和逆时偏移试验,展示了这种描述各向异性波场传播的新方程的特点与优势． 相似文献
3.
常规的单程波波动方程偏移成像方法对大角度的高陡构造偏移成像存在内在的限制.根据波动方程在各个空间方向的数学特性和高陡构造反射地震波的传播特征,通过把地震波分解为垂向的上下行波、水平方向的前后行波和左右行波,提出基于波场垂向外推和水平方向外推相结合的单程波波动方程高陡构造偏移成像方法,即用波场垂向外推的单程波波动方程偏移成像方法解决中低角度平缓构造的偏移成像,用波场水平方向外推的单程波波动方程偏移成像方法解决中高角度陡倾构造的偏移成像.这种基于波场垂向和水平方向外推相结合的高陡构造偏移成像方法是常规单程波波动方程叠前深度偏移成像方法的补充和改进,它相对基于全波方程的逆时偏移具有计算效率上的优势. 相似文献
4.
各向异性介质纯P波方程完全不受横波的干扰,在一定程度上可以减缓由于介质各向异性引起的数值不稳定,本文推导了具有垂直对称轴的横向各向同性(VTI)介质纯P波一阶速度-应力方程.由于纯P波方程存在一个分数形式的伪微分算子,无法直接采用有限差分法求解.针对该问题,本文采用伪谱法和高阶有限差分法联合求解波动方程,重点分析了混合法求解纯P波一阶速度-应力方程的稳定性问题,并给出了混合法求解纯P波方程的稳定性条件.数值模拟结果表明纯P波方程伪谱法和高阶有限差分混合法能够进行复杂介质的正演模拟,在强变速度、变密度的地球介质中仍然具有较好的稳定性. 相似文献
5.
3D P‐wave traveltime computation in transversely isotropic media using layer‐by‐layer wavefront marching 
下载免费PDF全文
下载免费PDF全文 Jiangtao Hu Junxing Cao Huazhong Wang Xingjian Wang Renfei Tian 《Geophysical Prospecting》2018,66(7):1303-1314
Subsurface rocks (e.g. shale) may induce seismic anisotropy, such as transverse isotropy. Traveltime computation is an essential component of depth imaging and tomography in transversely isotropic media. It is natural to compute the traveltime using the wavefront marching method. However, tracking the 3D wavefront is expensive, especially in anisotropic media. Besides, the wavefront marching method usually computes the traveltime using the eikonal equation. However, the anisotropic eikonal equation is highly non‐linear and it is challenging to solve. To address these issues, we present a layer‐by‐layer wavefront marching method to compute the P‐wave traveltime in 3D transversely isotropic media. To simplify the wavefront tracking, it uses the traveltime of the previous depth as the boundary condition to compute that of the next depth based on the wavefront marching. A strategy of traveltime computation is designed to guarantee the causality of wave propagation. To avoid solving the non‐linear eikonal equation, it updates traveltime along the expanding wavefront by Fermat's principle. To compute the traveltime using Fermat's principle, an approximate group velocity with high accuracy in transversely isotropic media is adopted to describe the ray propagation. Numerical examples on 3D vertical transverse isotropy and tilted transverse isotropy models show that the proposed method computes the traveltime with high accuracy. It can find applications in modelling and depth migration. 相似文献
6.
Depth migration consists of two different steps: wavefield extrapolation and imaging. The wave propagation is firmly founded on a mathematical frame-work, and is simulated by solving different types of wave equations, dependent on the physical model under investigation. In contrast, the imaging part of migration is usually based on ad hoc‘principles’, rather than on a physical model with an associated mathematical expression. The imaging is usually performed using the U/D concept of Claerbout (1971), which states that reflectors exist at points in the subsurface where the first arrival of the downgoing wave is time-coincident with the upgoing wave. Inversion can, as with migration, be divided into the two steps of wavefield extrapolation and imaging. In contrast to the imaging principle in migration, imaging in inversion follows from the mathematical formulation of the problem. The image with respect to the bulk modulus (or velocity) perturbations is proportional to the correlation between the time derivatives of a forward-propagated field and a backward-propagated residual field (Lailly 1984; Tarantola 1984). We assume a physical model in which the wave propagation is governed by the 2D acoustic wave equation. The wave equation is solved numerically using an efficient finite-difference scheme, making simulations in realistically sized models feasible. The two imaging concepts of migration and inversion are tested and compared in depth imaging from a synthetic offset vertical seismic profile section. In order to test the velocity sensitivity of the algorithms, two erroneous input velocity models are tested. We find that the algorithm founded on inverse theory is less sensitive to velocity errors than depth migration using the more ad hoc U/D imaging principle. 相似文献
7.
8.
Václav Bucha 《Studia Geophysica et Geodaetica》2012,56(2):533-552
We use Kirchhoff prestack depth migration to calculate migrated sections in 3-D simple anisotropic homogeneous velocity models in order to demonstrate the impact of anisotropy on migrated images. The recorded wave field is generated in models composed of two homogeneous layers separated by one either non-inclined or inclined curved interface. The anisotropy in the upper layer is triclinic. We apply Kirchhoff prestack depth migration to velocity models with different types of anisotropy: a triclinic anisotropic medium, an isotropic medium, transversely isotropic media with a horizontal (HTI) and vertical (VTI) symmetry axis. We observe asymmetry in migration caused by triclinic anisotropy and we show the errors of the migrated interface caused by inaccurate velocity models used for migration. The study is limited to P-waves. 相似文献
9.
Estimating azimuthal stress‐induced P‐wave anisotropy from S‐wave anisotropy using sonic log or vertical seismic profile data 下载免费PDF全文
下载免费PDF全文 Most sedimentary rocks are anisotropic, yet it is often difficult to accurately incorporate anisotropy into seismic workflows because analysis of anisotropy requires knowledge of a number of parameters that are difficult to estimate from standard seismic data. In this study, we provide a methodology to infer azimuthal P‐wave anisotropy from S‐wave anisotropy calculated from log or vertical seismic profile data. This methodology involves a number of steps. First, we compute the azimuthal P‐wave anisotropy in the dry medium as a function of the azimuthal S‐wave anisotropy using a rock physics model, which accounts for the stress dependency of seismic wave velocities in dry isotropic elastic media subjected to triaxial compression. Once the P‐wave anisotropy in the dry medium is known, we use the anisotropic Gassmann equations to estimate the anisotropy of the saturated medium. We test this workflow on the log data acquired in the North West Shelf of Australia, where azimuthal anisotropy is likely caused by large differences between minimum and maximum horizontal stresses. The obtained results are compared to azimuthal P‐wave anisotropy obtained via orthorhombic tomography in the same area. In the clean sandstone layers, anisotropy parameters obtained by both methods are fairly consistent. In the shale and shaly sandstone layers, however, there is a significant discrepancy between results since the stress‐induced anisotropy model we use is not applicable to rocks exhibiting intrinsic anisotropy. This methodology could be useful for building the initial anisotropic velocity model for imaging, which is to be refined through migration velocity analysis. 相似文献
10.
各向异性研究对地下介质精确成像有着重要的意义,在当前计算机硬件迅速发展及宽方位地震数据采集日益普遍的情况下,成像必须考虑介质的各向异性.逆时偏移是基于双程波动方程的较为精确的数值解的成像方法,所以相对于其他地震成像方法,它具有很大的优势,譬如不受反射界面的倾角限制、偏移速度结构合适时能够使回转波及多次波正确成像.在各向同性介质中,可使用标量波方程来模拟波场.而在各向异性介质中,P波和SV波是相互耦合的,即不存在单纯的标量波传播,通常利用能代表耦合波场中P波分量运动学特征的拟声波(qP波)进行偏移成像.本文中,我们推导出了TTI介质下qP波控制方程.该方程可采用显式有限差分格式进行求解.通过声学近似,若沿对称轴方向的剪切波速度为零,对于对称轴方向不变且ε≥δ的模型来说,可得到稳定的数值解.但对于TTI介质来说,由于沿对称轴方向各向异性参数是变化的,声学近似会引起波场传播及数值计算的不稳定.因此,我们提出了正则化有限横波的方法,很好地解决了这一问题.最后,给出了Foothill模型的测试结果及某探区实际资料试算结果,展示了采用这个方程进行复杂TTI模型正演和高质量逆时偏移成像结果,证实了该方法的正确性和实际资料应用中的有效性. 相似文献
11.
Characterizing the expressions of seismic waves in elastic anisotropic media depends on multiparameters. To reduce the complexity, decomposing the P-mode wave from elastic seismic data is an effective way to describe the considerably accurate kinematics with fewer parameters. The acoustic approximation for transversely isotropic media is widely used to obtain P-mode wave by setting the axial S-wave phase velocity to zero. However, the separated pure P-wave of this approach is coupled with undesired S-wave in anisotropic media called S-wave artefacts. To eliminate the S-wave artefacts in acoustic waves for anisotropic media, we set the vertical S-wave phase velocity as a function related to propagation directions. Then, we derive a pure P-wave equation in transversely isotropic media with a horizontal symmetry axis by introducing the expression of vertical S-wave phase velocity. The differential form of new expression for pure P-wave is reduced to second-order by inserting the expression of S-wave phase velocity as an auxiliary operator. The results of numerical simulation examples by finite difference illustrate the stability and accuracy of the derived pure P-wave equation. 相似文献
12.
经过多年的研究发展,各向异性叠前深度偏移算法已经趋于完善.然而,在地震资料处理过程中导致成像结果不理想的主要原因还是由于建立的地层参数场不够精确.当地层参数接近其真实值时,基于波动方程的剩余曲率建模方法由于不受构造的影响,能够在各向异性和横向变速介质中进行速度分析,所以得到了广泛的研究.本文从偏移结果中抽取共成像道集,然后通过交互运用叠前深度偏移和参数更新实现各向异性偏移速度建模.对理论模型和实际资料进行的试算表明,该方法具有较强的适应性,能极大改善VTI介质反射界面成像效果和分辨率. 相似文献
13.
Pure-mode wave propagation is important for applications ranging from imaging to avoiding parameter tradeoff in waveform inversion. Although seismic anisotropy is an elastic phenomenon, pseudo-acoustic approximations are routinely used to avoid the high computational cost and difficulty in decoupling wave modes to obtain interpretable seismic images. However, such approximations may result in inaccuracies in characterizing anisotropic wave propagation. We propose new pure-mode equations for P- and S-waves resulting in an artefact-free solution in transversely isotropic medium with a vertical symmetry axis. Our approximations are more accurate than other known approximations as they are not based on weak anisotropy assumptions. Therefore, the S-wave approximation can reproduce the group velocity triplications in strongly anisotropic media. The proposed approximations can be used for accurate modelling and imaging of pure P- and S-waves in transversely isotropic media. 相似文献
14.
Di Luo Zhiqiang Wu Tianjiu Zeng Xunhua Zhang Weina Zhao Huaishan Liu Guolin Xiao 《Studia Geophysica et Geodaetica》2018,62(3):466-484
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. 相似文献
15.
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. 相似文献
16.
由所建立的三维qP波相速度表示式出发，导出并解析求解各向异性介质中的频散方程，得到三维各向异性介质中的相移算子，进而将以相移算子为基础的对称非平稳相移方法推广到各向异性介质，发展了一个三维各向异性介质的深度偏移方法. 文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致，将各向同性、弱各向异性及强各向异性统一在一个模型中. 所建立的各向异性介质对称非平稳相移波场延拓算子可以同时适应速度及各向异性参数横向变化；文中给出的算例虽然是针对二维VTI介质的，但所提出的算法同样适用于三维TI介质. 相似文献
17.
地震各向异性研究进展 总被引：28，自引：21，他引：28
简要介绍了各向异性介质中地震波传播理论的最新研究进展，总结了地震各向异性在地球动力学，地震灾害监测以及资源探查等方面的应用研究新进展。目前，主要通过数值模拟和物理模拟来研究地球内部典型各向异性介质中波的传播规律；利用横波分裂，Pn波的方位变化，面波成像方法研究地球内部的结构和物质交换，应力场的变化以及地震的孕育和发生过程。在资源探查中，地震各向异性理论主要用于解决测线不闭合，中长排列动校，时深转换以及裂缝分布等问题。 相似文献
18.
Pure acoustic wave propagation in transversely isotropic media by the pseudospectral method 总被引：4，自引：0，他引：4
Seismic wave propagation in transversely isotropic (TI) media is commonly described by a set of coupled partial differential equations, derived from the acoustic approximation. These equations produce pure P‐wave responses in elliptically anisotropic media but generate undesired shear‐wave components for more general TI anisotropy. Furthermore, these equations suffer from instabilities when the anisotropy parameter ε is less than δ. One solution to both problems is to use pure acoustic anisotropic wave equations, which can produce pure P‐waves without any shear‐wave contaminations in both elliptical and anelliptical TI media. In this paper, we propose a new pure acoustic transversely isotropic wave equation, which can be conveniently solved using the pseudospectral method. Like most other pure acoustic anisotropic wave equations, our equation involves complicated pseudo‐differential operators in space which are difficult to handle using the finite difference method. The advantage of our equation is that all of its model parameters are separable from the spatial differential and pseudo‐differential operators; therefore, the pseudospectral method can be directly applied. We use phase velocity analysis to show that our equation, expressed in a summation form, can be properly truncated to achieve the desired accuracy according to anisotropy strength. This flexibility allows us to save computational time by choosing the right number of summation terms for a given model. We use numerical examples to demonstrate that this new pure acoustic wave equation can produce highly accurate results, completely free from shear‐wave artefacts. This equation can be straightforwardly generalized to tilted TI media. 相似文献
19.
利用地震波正向传播方程对属于波形线性反演问题近似求解方法的地震数据偏移成像进行重新推导,得到了适合散射地震数据的散射偏移成像方法和适合反射地震数据的反射偏移成像方法.以地震波传播的散射理论为出发点,首先根据描述一次散射波正向传播的线性方程研究建立散射地震数据的偏移成像方法理论;利用高频近似对产生散射波场的地下速度扰动函数的空间变化进行近似,推导出地下反射率函数,再由散射波传播方程推导出基于反射率函数的反射波传播方程,然后根据描述一次反射波正向传播的线性方程研究建立反射地震数据的偏移成像方法理论.本文指出和修正了Claerbout偏移成像方法中的不足,提出的地震数据偏移成像方法是对当前偏移成像方法理论的完善,使反射地震数据偏移成像具有了更坚实的数学物理理论基础,得到的偏移成像结果相位正确、位置准确、分辨率提高. 相似文献
20.
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. 相似文献