共查询到19条相似文献，搜索用时 328 毫秒
1.
探讨了井巷隧道工程反射波法超前地质预报成像问题.从一阶速度－应力弹性波方程出发,推导了二维各向同性介质情况下弹性波逆时传播的高阶差分格式,实现了弹性波在数值空间中的逆时延拓.构建反射波法隧道超前探测中的断层、软弱夹层等介质模型,以反射波探测的正演记录作为初始条件,并从程函方程出发,采用逆时差分格式求取介质模型网格空间中各点的直达波旅行时作为弹性波逆时偏移的成像条件,实现多波多分量资料的逆时偏移.偏移结果表明,逆时偏移能够使隧道壁接收到的波场准确归位,提高隧道反射波超前探测的资料处理的精度. 相似文献
2.
各向异性弹性波动方程多分量联合叠后逆时偏移 总被引：6，自引：1，他引：5
为了更好地描述地下介质的各向异性特性,实现准确的波场成像,开展了各向异性介质逆时偏移研究。采用高阶交错网格有限差分法推导出各向异性弹性波动方程叠后逆时深度偏移方程,证明各向异性介质向各向同性介质转化的条件,并研究Thom son参数对弹性波场的影响,合成了各向同性介质和各向异性介质多分量零偏移距剖面,实现了多分量模拟记录联合叠后逆时偏移。计算结果表明,能够准确地实现多分量波场叠后逆时成像,使地质层位中的断层、断点等目标成像清晰准确,且偏移成像精度较高,因此各向异性逆时偏移法是一种高精度的深度域偏移成像方法,可以指导实际天然地震资料和人工地震资料的数据处理。 相似文献
3.
《地球物理学进展》2015,(4)
基于双程波动方程的叠前逆时偏移方法成像精度高,而且无地层倾角限制,较适合于复杂地下构造成像.但是,由于地下介质的各向异性广泛存在,基于各向同性的正演算法,尚难准确描述真实的地下波场传播,逆时偏移的成像精度也因此受到限制.鉴于此,本文研究了各向异性VTI介质逆时偏移方法,首先根据VTI介质一阶准P波方程推导出了炮点和检波点的逆时延拓的交错网格高阶差分格式,针对算法计算量和存储量大的问题,文中研究了一种改进的基于GPU加速的有效PML边界存储策略.本文建议的方法只需增加少量的额外计算,就可降低大幅度的存储成本,进而实现高精度和高效率的各向异性逆时偏移.Hess 2DVTI模型测试表明,本文提出的方法不需要存储全部历史时刻的波场,可以实现高效率高精度的VTI介质叠前逆时偏移成像. 相似文献
4.
研究了三维弱各向异性近似下,利用伪P波(伪纵波)模拟弹性波场P分量在倾斜对称轴的横向各向同性(TTI)介质中的传播过程,并对比了分别基于弹性Hooke定律、弹性波投影和运动学色散方程所建立的三种二阶差分伪P波方程的正演特点.目前这些伪P波方程数值计算主要采用规则网格差分,但是规则网格在TTI模拟中有低效率、低精度以及不稳定的缺点.为了提高计算的精度,本文构建出相应方程的交错网格有限差分格式.通过对比伪P波方程在三维TTI介质中不同的数值模拟的表达形式,本文认为基于色散方程所建立的伪P波方程在模拟弹性波中P波传播的过程中具有最小的噪声.本文分析不同的各向同性对称轴空间角度的频散特征,并引入适当的横波速度维持计算的稳定.二维模型算例表明,本文提出的交错网格正演算法可以得到稳定光滑的伪P波正演波场.使用本文交错网格算法对二维BP TTI模型的逆时偏移也具有较稳定的偏移结果. 相似文献
5.
《地球物理学进展》2017,(4)
弹性波逆时偏移是当前多分量地震资料相对准确的偏移算法,它能够形成多波模式的成像剖面,从而减少纵波勘探的多解性.本文首先依据各向同性介质中矢量分离纵横波场的速度-应力方程组,利用高阶交错网格有限差分数值方法求解弹性波方程,进而构建矢量的纵波和横波波场,不同于散度和旋度算子分离纵横波场的传统方法,文中提出的矢量分离纵横波场方法保持了原始波场的振幅和相位特征.文中也提出将震源归一化的内积成像条件应用于分离后的纯纵波和横波矢量场,由此得到的转换波成像避免了传统弹性波成像方法中出现的极性反转.水平层状和复杂构造模型测试表明,文中提出的基于矢量分离纵横波场的弹性波逆时偏移方法成像精度高,转换波成像PS和SP极性无反转,所形成的多种模式纯波剖面能够准确地对复杂地下构造成像. 相似文献
6.
多分量联合逆时偏移最佳匹配层吸收边界 总被引：3，自引：2，他引：1
有限空间内的波动方程逆时偏移需要利用有效的边界处理技术用以消除人工截断对偏移结果产生的影响。本文以横向各向同性介质弹性波速度-应力方程为基础,依据传统分裂式最佳匹配层（Perfect Matched Layer,PML）吸收边界技术的思想,推导了应用于逆时偏移的完全匹配层波动方程,并给出了其高阶交错网格有限差分格式。针对由边界处向计算区域内传播的＂反射波＂,以及地震记录排列两端地震同相轴突变对计算区域的影响这两方面问题,本文给出了逆时偏移中吸收层的布设方式。模型和实际资料的弹性波叠前多分量联合逆时深度偏移结果表明本文的边界处理方法取得了较好的吸收效果,获得了好的联合偏移成像结果。 相似文献
7.
根据具有垂直对称轴的横向各向同性（VTI）介质中的一阶准P波方程,导出了该方程在交错网格中逆时延拓的高阶有限差分格式,给出了其稳定性条件,采用完全匹配层吸收边界条件解决边界反射问题,分别应用下行波最大能量法和归一化互相关成像条件, 实现了VTI介质中准P波方程的叠前逆时深度偏移．各向异性Marmousi模型的试算结果表明,VTI介质准P波方程叠前逆时深度偏移算法不受地下构造倾角和介质横向速度变化的限制,对复杂模型具有良好的成像能力;应用归一化互相关成像条件能得到更好的成像效果．对比该模型的各向异性和各向同性逆时偏移剖面表明,在各向异性地区采集的纵波数据用各向异性偏移算法理论上能得到更好的成像结果． 相似文献
8.
《应用地球物理》2019,(3)
在海洋地震勘探中,海底电缆技术采集的多分量地震数据,涉及到多参数的反演与成像问题,本文针对海底电缆多分量地震数据提出了一种弹性波多参数最小二乘逆时偏移方法。该方法的波场延拓算子为混合方程,即在海水介质中采用声波方程进行波场计算,而海底固体介质的波场由纵横波形分离的矢量弹性波方程得到。在海底界面中采用声弹耦合控制方程将两种类型的方程结合起来。通过推导纵横波形分离的弹性波偏移算子、反偏移算子和梯度公式实现基于纵横波形分离的弹性波最小二乘逆时偏移方法。通过模型试算证明了该方法能够得到高质量的纵波速度和横波速度分量的成像剖面,相比与传统弹性波最小二乘逆时偏移方法,多参数耦合造成的成像串扰噪音得到了很好地压制。 相似文献
9.
《地球物理学进展》2016,(1)
为克服各向同性和VTI介质逆时偏移方法对复杂地质构造成像的局限,研究了TTI介质拟声波逆时偏移方法.首先从精确的TTI介质频散关系出发,引入一个各向异性控制参数σ,推导了新的二阶耦合TTI介质拟声波方程,以保证波场延拓的稳定性;然后引入波场的伪速度分量,推导了等价的一阶拟声波方程.相比于规则网格有限差分法,交错网格有限差分(SGFD)法能够有效地压制数值频散,模拟精度更高;因此利用高阶SGFD法求解TTI介质一阶拟声波方程,构建逆时偏移所需的正向和逆时波场延拓算子,并应用归一化互相关成像条件实现精确的TTI介质逆时偏移成像.最后,简单讨论了伪横波的产生机制,并给出了伪横波的联合压制策略.模型试验结果验证了方法的有效性和稳定性. 相似文献
10.
逆时偏移作为一种高精度偏移方法已成为复杂构造成像的重要技术,描述纵波独立传播的延拓方程是各向异性介质逆时偏移的一个关键问题.在对VTI介质几个经典相速度近似公式回顾的基础上,针对常用于描述纯P波的Harlan近似公式在各向异性参数ε较大情况下近似精度较低的问题,本文对Harlan公式中的非椭圆项进行了修正,在非椭圆项前添加了一个与各向异性参数ε有关的修正系数,得到了三种改进型Harlan公式,并以近似精度最高的改进式为基础,推导了TTI介质纯P波方程.针对该伪微分方程,本文利用伪谱法和有限差分法联合实现波场延拓,对于常密度二阶方程,基于中心网格实现;对于一阶应力-速度方程则基于旋转交错网格实现.通过数值试验分析了TTI介质纯P波一阶应力-速度方程的近似精度,并以一阶纯P波方程为基础进行了TTI介质逆时偏移数值模拟试验.结果表明,本文给出的方法能够较准确地描述TTI介质纯P波波场特征,可以应用至各向异性介质逆时偏移. 相似文献
11.
When treating the forward full waveform case, a fast and accurate algorithm for modelling seismic wave propagation in anisotropic inhomogeneous media is of considerable value in current exploration seismology. Synthetic seismograms were computed for P-SV wave propagation in transversely isotropic media. Among the various techniques available for seismic modelling, the finite-difference method possesses both the power and flexibility to model wave propagation accurately in anisotropic inhomogeneous media bounded by irregular interfaces. We have developed a fast high-order vectorized finite-difference algorithm adapted for the vector supercomputer. The algorithm is based on the fourth-order accurate MacCormack-type splitting scheme. Solving the equivalent first-order hyperbolic system of equations, instead of the second-order wave equation, avoids computation of the spatial derivatives of the medium's anisotropic elastic parameters. Examples indicate that anisotropy plays an important role in modelling the kinematic and the dynamic properties of the wave propagation and should be taken into account when necessary. 相似文献
12.
针对地下工程领域隧道超前预报地震波波场传播与成像中存在的问题,通过数值模拟,构建二维含低速异常的隧道介质模型,研究隧道弹性波场传播规律和异常体边界成像准确性.首先,利用一阶速度-应力波动方程和高阶交错网格有限差分计算方法,导出隧道超前预报数值模拟的稳定性条件和边界条件,对上述隧道模型进行数值模拟,识别波场特征;其次,利用叠前逆时偏移成像方法,对压制噪音干扰后的波场在互相关成像条件下,对隧道模型中的异常体边界进行逆时偏移成像.研究结果表明：采用高阶交错网格有限差分正演获得异常体边界清晰的反射波和角点产生的散射波;逆时偏移算法获得隧道内异常体准确成像结果,从而大大提高隧道超前预报的分辨率与准确性;靠近掌子面单一震源、多道接收观测系统对异常体成像效果最佳,为隧道内高效数据采集提供理论依据. 相似文献
13.
Viscoacoustic prestack reverse time migration based on the optimal time-space domain high-order finite-difference method 总被引：3，自引：1，他引：2
Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution. 相似文献
14.
提出了一种新的叠前时间偏移方法和流程,可不必应用野外静校正,直接对起伏地表采集的地震数据进行叠前时间偏移.本文采用输入道成像方式,通过基于稳相点原理给出单道数据的走时和振幅计算方法,发展了一个表驱动的叠前时间偏移算法.偏移方法可依据同相轴是否被拉平确定叠加速度和修正近地表速度模型,也可依据拟成像的构造倾角,自适应地确定偏移孔径;后者既减少了偏移计算量,也压制了偏移噪声.文中用二维起伏地表的断陷盆地模型的理论数据验证了所发展方法的成像效果. 相似文献
15.
转换波偏移可以利用纵横波波场信息,得到高分辨率的成像结果,从而为油藏描述提供高质量的地震资料.目前的研究主要是利用纵波波场信息进行偏移成像,然而,传统的纵波方法在复杂探区成像时具有一定的局限性.为此,本文在各向异性介质声波射线追踪算法的基础上,推导出各向异性介质转换波射线追踪方程,发展了一种转换波射线追踪算法;并将研究的追踪算法应用到偏移成像中,提出了一种各向异性VTI介质角度域转换波高斯束偏移成像方法.通过各向异性VTI介质断块模型和复杂构造模型试算,说明了本文方法的正确性和有效性.模型试算的结果表明,在考虑地下各向异性时,本文研究的方法具有更好的成像效果,提取的角道集结果可以为偏移速度分析提供依据. 相似文献
16.
各向异性介质纯P波方程完全不受横波的干扰,在一定程度上可以减缓由于介质各向异性引起的数值不稳定,本文推导了具有垂直对称轴的横向各向同性(VTI)介质纯P波一阶速度-应力方程.由于纯P波方程存在一个分数形式的伪微分算子,无法直接采用有限差分法求解.针对该问题,本文采用伪谱法和高阶有限差分法联合求解波动方程,重点分析了混合法求解纯P波一阶速度-应力方程的稳定性问题,并给出了混合法求解纯P波方程的稳定性条件.数值模拟结果表明纯P波方程伪谱法和高阶有限差分混合法能够进行复杂介质的正演模拟,在强变速度、变密度的地球介质中仍然具有较好的稳定性. 相似文献
17.
经过多年的研究发展,各向异性叠前深度偏移算法已经趋于完善.然而,在地震资料处理过程中导致成像结果不理想的主要原因还是由于建立的地层参数场不够精确.当地层参数接近其真实值时,基于波动方程的剩余曲率建模方法由于不受构造的影响,能够在各向异性和横向变速介质中进行速度分析,所以得到了广泛的研究.本文从偏移结果中抽取共成像道集,然后通过交互运用叠前深度偏移和参数更新实现各向异性偏移速度建模.对理论模型和实际资料进行的试算表明,该方法具有较强的适应性,能极大改善VTI介质反射界面成像效果和分辨率. 相似文献
18.
Based on arbitrarily wide-angle wave equations,a reverse-time propagation scheme is developed by substituting the partial derivatives of depth and time with central differences. The partial derivative of horizontal direction is replaced with high order difference. The imaging condition is computed by solving the eikonal equations. On the basis of above techniques,a prestack reverse-time depth migration algorithm is developed. The processing exam-ples of synthetic data show that the method can remove unwanted internal reflections and decrease the migration noise. The method also has the advantage of fidelity and is applicable of dip angle reflector imaging. 相似文献
19.
We present a layer-stripping method of migration for irregularly layered media in which first-order velocity discontinuities separate regions of constant or smoothly varying velocity. We use the reverse-time method to migrate seismic data layer by layer, from the surface downwards. As part of the migration of a given layer, the bottom boundary of the layer is defined based on power in the migrated signal, and a seismic section is collected along it. This new section serves as the boundary condition for migration in the next layer. This procedure is repeated for each layer, with the final image formed from the individual layer images. Layer-stripping migration consists of three steps: (1) layer definition, (2) wavefield extrapolation and imaging, and (3) boundary determination. The migration scheme when used with reverse-time extrapolation is similar to datuming with an imaging condition. The reverse-time method uses an explicit fourth-order time, tenth-order space, finite-difference approximation to the scalar wave equation. The advantages of layer-stripping reverse-time migration are: (1) it preserves the benefits of the reverse-time method by handling strong velocity contrasts between layers and steeply dipping structures; (2) it reduces computer memory and saves computation time in high-velocity layers, and (3) it allows interpretational control of the image. Post-stack layer-stripping reverse-time migration is illustrated with a synthetic CMP data example. Prestack migration is illustrated with a synthetic data set and with a marine seismic reflection profile across the Santa Maria Basin and the Hosgri Fault in central California. 相似文献