共查询到19条相似文献，搜索用时 218 毫秒
1.
基于双程波动方程的叠前逆时偏移技术是解决复杂构造成像难题的强有力手段.本文首先对该技术的基本原理和实现方法进行简要介绍,该技术主要包括震源波场与检波波场构建、对构建的源、检波场应用互相关成像条件以及成像过程中或成像后低波数噪声压制三个主要环节.在此基础上对其实用化过程中涉及的高性能计算策略、低存储方案和高效低波数噪声压制方法等关键环节所需考虑的因素进行分析并给出相应的解决方案,从而给出了一套存储消耗较低且计算效率较高的逆时偏移实现流程.然后以Sigsbee2A模型和深水盐丘数据为例,与工业上标准的克希霍夫叠前深度偏移方法成像结果进行对比从而验证方法的正确性及对复杂构造的成像优势.最后结合理论分析与实际测试结果对该技术的特点和发展方向进行探讨与总结. 相似文献
2.
3.
最小二乘逆时偏移（LSRTM）相对于常规逆时偏移（RTM）具有分辨率更高、振幅更准确、噪音更少等优势,可以对复杂的地质构造进行有效的成像.这种迭代更新反演成像方法十分依赖目标函数的梯度质量和计算效率.当地质模型中存在强反射界面或者记录中存在折射波时,基于常规互相关成像条件（CCC）的最小二乘逆时偏移梯度会包含很强的低频噪音,从而使反演的收敛速度和成像质量降低.为此,本文在最小二乘逆时偏移的梯度中引进了逆散射成像条件来压制这种低频噪音,并以此提出基于逆散射成像条件（ISC）的最小二乘逆时偏移方法.数值模拟结果表明,两者计算耗时基本一致,但逆散射成像条件能高效压制梯度中的低频噪音,从而使反演过程中收敛加速,成像质量得到显著提高. 相似文献
4.
地震波场模拟是叠前逆时偏移方法的核心计算单元,因此,边界条件和存储方案的选择是高精度以及高效率逆时偏移算法不可回避的问题.本文就随机边界条件和吸收边界条件在逆时偏移算法中的应用效果以及相应的四种存储策略展开具体分析,给出了不同存储策略下的计算成本和存储量需求.文中优选随机边界条件和吸收边界条件的有效边界存储方案,引入GPU并行加速技术对SEG起伏地表模型进行测试,详细讨论分析了两种方法的计算效率和计算精度.测试表明,随机边界存储的计算成本要低于有效边界存储的计算成本,而受边界散射的影响,随机边界逆时偏移在浅层伴有噪音干扰.针对具体数据情况,合理地运用GPU加速技术、选择有效的边界条件以及存储方法,是改善逆时偏移成像精度和效率的有效途径. 相似文献
5.
本文深入分析了多次波逆时偏移的成像原理和串扰假象产生机制,并提出了多次波分阶逆时偏移策略,即首先应用自由界面多次波衰减方法对原始炮集记录中的多次波进行剔除得到一次波记录,然后应用一次波记录预测获得各阶多次波记录,最后将各阶多次波记录分别进行逆时偏移成像.模型实验结果表明,多次波分阶逆时偏移其各阶多次波成像剖面均可对真实界面正确成像且其能够有效压制串扰假象,其中一阶多次波逆时偏移剖面的成像精度最高,其深层成像质量明显优于常规多次波逆时偏移. 相似文献
6.
7.
相比地面地震勘探,垂直地震剖面(VSP)资料信噪比高,从VSP偏移结果中能够准确地识别井旁小的地质构造.在VSP偏移方法中,VSP逆时偏移方法基于双程波动方程,具有成像精确且无倾角限制、适应任意复杂速度模型等诸多优点,受到人们广泛关注.传统的VSP逆时偏移方法会在成像剖面上产生低频噪声,同时井中检波点分布位置附近会产生井轨迹噪声,严重影响了成像剖面的质量.在本文中,我们利用优化时空域频散关系求解差分系数,提高数值模拟精度;采用混合吸收边界条件压制边界反射,减少边界反射对成像的影响;通过Hilbert变换来构建复数波场,利用波场分解成像条件进行VSP逆时偏移成像,来压制低频噪声和井轨迹噪声.数值模型和实际资料测试表明,基于复数波场分解成像方法能够实现高精度的VSP逆时偏移,并且能够有效地压制VSP逆时偏移中的低频成像噪声和井轨迹噪声;相比地面地震逆时偏移,本文方法能够对井旁构造和盐丘下构造取得更清晰的成像结果. 相似文献
8.
叠前逆时偏移假象去除方法 总被引：3，自引：0，他引：3
叠前逆时偏移(RTM)求解全波动方程进行波场延拓,可以精确描述地震波的传播,在陡倾界面和复杂构造成像方面具有显著优势.但在实际应用中,互相关成像条件会带来低波数高振幅的假象.本文分析了假象产生的原因,实现并比较了Laplace滤波、坡印廷矢量法和上下、左右行波分解这三种假象压制方法,针对不同的计算平台分析其优劣.数值试验表明,坡印廷矢量法当波场简单时较为有效;而上下、左右行波分解法和Laplace滤波可以应对较复杂的波场.就去假象结果而言,使用波场分解成像条件并施加Laplace滤波效果最佳;然而就计算效率而言,CPU平台应用波场分解成像条件只增加少量的计算量,但GPU平台的实现则需要对数据的读写付出较大的代价.因此,综合考虑,对于CPU平台上述两种方案的组合为最佳策略;而基于GPU平台的逆时偏移采用Laplace滤波去噪仍然是目前最经济的选择.Marmousi模型和SEG/HESS VTI模型试验验证了上述结论. 相似文献
9.
与地面地震资料相比,VSP资料具有分辨率高、环境噪声小及能更好地反映井旁信息等优点.常规VSP偏移主要对上行反射波进行成像,存在照明度低、成像范围受限等问题.为了增加照明度、拓宽成像范围、提高成像精度,本文采用直达波除外的所有声波波场数据(全波),包括一次反射波、多次反射波等进行叠前逆时偏移成像.针对逆时偏移中的四个关键问题,即波场延拓、吸收边界条件、成像条件及低频噪声的压制,本文分别采用自适应变空间差分算子长度的优化有限差分方法(自适应优化有限差分方法)求解二维声波波动方程以实现高精度、高效率的波场延拓,采用混合吸收边界条件压制因计算区域有限所引起的人工边界反射,采用震源归一化零延迟互相关成像条件进行成像,采用拉普拉斯滤波方法压制逆时偏移中产生的低频噪声.本文对VSP模型数据的逆时偏移成像进行了分析,结果表明:自适应优化有限差分方法比传统有限差分方法具有更高的模拟精度与计算效率,适用于VSP逆时偏移成像;全波场VSP逆时偏移成像比上行波VSP逆时偏移的成像范围大、成像效果好;相对于反褶积成像条件,震源归一化零延迟互相关成像条件具有稳定性好、计算效率高等优点.将本文方法应用于某实际VSP资料的逆时偏移成像,进一步验证了本文方法的正确性和有效性. 相似文献
10.
《地球物理学进展》2016,(1)
为克服各向同性和VTI介质逆时偏移方法对复杂地质构造成像的局限,研究了TTI介质拟声波逆时偏移方法.首先从精确的TTI介质频散关系出发,引入一个各向异性控制参数σ,推导了新的二阶耦合TTI介质拟声波方程,以保证波场延拓的稳定性;然后引入波场的伪速度分量,推导了等价的一阶拟声波方程.相比于规则网格有限差分法,交错网格有限差分(SGFD)法能够有效地压制数值频散,模拟精度更高;因此利用高阶SGFD法求解TTI介质一阶拟声波方程,构建逆时偏移所需的正向和逆时波场延拓算子,并应用归一化互相关成像条件实现精确的TTI介质逆时偏移成像.最后,简单讨论了伪横波的产生机制,并给出了伪横波的联合压制策略.模型试验结果验证了方法的有效性和稳定性. 相似文献
11.
研究不同偏移方法的成像机理是实现复杂构造高精度成像的前提,研究了两类共炮点域的相关型叠前深度偏移成像方法：基于单程波动方程的叠前深度成像方法和基于双程波动方程的叠前深度成像方法,同时对比了它们在计算效率、数据存储量、成像精度、成像机理、速度敏感性等方面的差异及其共性。以复杂构造模型为例,采用了傅里叶有限差分法（FFD）和逆时成像法（RTM）,这两种方法实现了121个共炮点道集的叠前深度偏移成像处理。计算结果表明,当速度准确时,两种深度域波动方程成像方法均可以恢复出各个地质反射界面,其中逆时偏移对陡倾角成像效果显著,当速度存在百分比误差和随机扰动情况时,逆时成像结果要差于单程波方法,因此,逆时偏移方法对速度的敏感性较大,且低频噪声较为严重。 相似文献
12.
叠前逆时偏移技术是解决地震成像问题的有力工具,但由于计算量大、成像噪音以及存储量大等原因没有得到广泛的应用.本文给出了逆时偏移的实现过程,分析了高阶有限差分格式的稳定性与频散关系.针对叠前逆时偏移计算量大的问题,使用图形处理器(Graphic Processing Unit,简称GPU)实现算法加速,比传统的CPU计算速度提高了一个数量级.文中对理论模型进行了计算,并与单程波偏移方法做比较,结果表明:叠前逆时偏移有效突破了成像倾角限制,对垂直断层、盐丘空腔内幕等特殊构造成像效果均有显著提高.本文尚未涉及成像噪音去除以及存储量等问题,笔者将另文阐述. 相似文献
13.
逆时偏移被认为是对地下复杂构造进行成像的精确偏移方法,尤其是能够有效地对回转波、绕射波、多次波等各种波动现象进行成像.近几年来随着并行计算机和存储设备的快速发展,逆时偏移方法备受关注.本文采用伪谱法实现了VSP逆时偏移,该方法不仅实施简便,而且计算效率高,精度好.并运用反周期扩展法来消除伪谱法中特殊的周期性边界效应问题.对VSP绕射点模型进行试算,分析了因覆盖次数不足在近井区产生的假象问题.对地堑模型和半圆隆起模型也获得了较好的VSP逆时偏移成像效果.并分别对VSP全波波场及分离出的上行波场进行了逆时偏移成像,可明显发现直达波在炮点和检波点位置处收敛成像,也产生了很强的成像噪声.最后对某地区实际观测的VSP资料进行了逆时偏移成像,并与Kirchhoff法VSP偏移结果和地面地震偏移结果进行了对比,显示了VSP逆时偏移在近井区成像上的优势. 相似文献
14.
声反射成像测井中常用的基于射线理论的Kirchhoff积分偏移算法和基于单程波理论的F-K偏移算法均可实现井旁缝洞反射体的快速偏移成像,但其仅适用于地层垂向变化较弱的速度场和高陡角裂缝的偏移成像,无法实现低角度反射体的准确偏移归位,产生偏移假象误导测井解释.逆时偏移基于全波动方程,可适应强垂向变化速度场,实现近似水平反射体的偏移成像.本文详细分析了将逆时偏移应用于声反射成像测井时存在的数据准备、时间采样间隔匹配和成像条件改进等若干问题,通过设置多组理论模型来说明算法对井旁不同反射体的识别能力.模拟资料和实际资料处理结果证实,较F-K偏移算法,逆时偏移算法成像精度更高、收敛性更好,可有效实现近似水平构造偏移归位.改进的归一化互相关成像条件可解决深部地层的远井壁成像衰减问题,降低测井解释的多解性.逆时偏移将成为声反射成像测井高精度偏移技术的发展方向. 相似文献
15.
逆时偏移成像(RTM)常用来处理复杂速度模型,包括陡倾角及横向速度变化剧烈的模型.与常规偏移成像方法(如Kirchhoff偏移)相比,逆时偏移成像能提供更好的偏移成像结果,近些年逆时偏移成像越来越广泛地应用到勘探地震中,它逐渐成为石油地震勘探中的一种行业标准.电磁波和弹性波在动力学和运动学上存在相似性,故本文开发了基于麦克斯韦方程组的电磁波逆时偏移成像算法,并将其应用到探地雷达数据处理中.时间域有限差分(FDTD)用于模拟电磁波正向和逆向传播过程,互相关成像条件用于获得最终偏移结果.逆时偏移成像算法中,偏移成像结果受初始模型影响较大,而其中决定电磁波传播速度的介电常数的影响尤为重要.本文基于时间域全波形反演(FWI)算法反演获得了更为精确的地下介电常数模型,并将其反演结果作为逆时偏移成像的初始介电常数模型.为了验证此算法的有效性,首先构建了一个复杂地质结构模型,合成了共偏移距及共炮点探地雷达数据,分别应用常规Kirchhoff偏移算法及逆时偏移成像算法进行偏移处理,成像结果显示由逆时偏移成像算法得到的偏移结果与实际模型具有较高的一致性;此外本文在室内沙槽中进行了相关的物理模拟实验,采集了共偏移距及共炮点探地雷达数据,分别应用Kirchhoff和叠前逆时偏移成像算法进行处理,结果表明叠前逆时偏移成像在实际应用中能获得更好的成像效果. 相似文献
16.
在常规偏移方法中一般都需要压制地震数据中的多次波,仅利用一次波信息成像,把自由表面反射的多次波视为噪声,但是在多次波中也包含着地下结构信息,应该将其充分利用到成像中来.事实上,已经有不少成像方法试图利用多次波信息,但是大部分方法都需要对多次波进行预测.本文提出了基于傅里叶有限差分偏移算子的数据自相关偏移方法.在这种偏移方法中,对含有一次波和多次波的地震数据,分别进行下行和上行延拓,然后直接利用常规的互相关成像条件成像.由于波场延拓采用了傅里叶有限差分算子,其计算效率高,能够很好地对复杂介质中的地震数据进行延拓.在数值试验中,使用了一个含散射点的三层模型和Marmousi模型.合成数据测试结果表明,这种方法可以对更大范围的地下构造成像,比常规的只利用一次波的傅里叶有限差分法照明度更好,并且在浅层可以提供更高的分辨率.我们提出的数据自相关策略易于实现且避免了繁杂的多次波预测,这对于复杂地下构造成像可能有着重大意义. 相似文献
17.
It is important to include the viscous effect in seismic numerical modelling and seismic migration due to the ubiquitous viscosity in an actual subsurface medium. Prestack reverse‐time migration (RTM) is currently one of the most accurate methods for seismic imaging. One of the key steps of RTM is wavefield forward and backward extrapolation and how to solve the wave equation fast and accurately is the essence of this process. In this paper, we apply the time‐space domain dispersion‐relation‐based finite‐difference (FD) method for visco‐acoustic wave numerical modelling. Dispersion analysis and numerical modelling results demonstrate that the time‐space domain FD method has great accuracy and can effectively suppress numerical dispersion. Also, we use the time‐space domain FD method to solve the visco‐acoustic wave equation in wavefield extrapolation of RTM and apply the source‐normalized cross‐correlation imaging condition in migration. Improved imaging has been obtained in both synthetic and real data tests. The migration result of the visco‐acoustic wave RTM is clearer and more accurate than that of acoustic wave RTM. In addition, in the process of wavefield forward and backward extrapolation, we adopt adaptive variable‐length spatial operators to compute spatial derivatives to significantly decrease computing costs without reducing the accuracy of the numerical solution. 相似文献
18.
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. 相似文献
19.
Jia-Jia Yang Xi-Wu Luan Bing-Shou He Gang Fang Jun Pan Wei-Min Ran Tao Jiang 《应用地球物理》2017,14(4):492-504
Angle-domain common-image gathers (ADCIGs) transformed from the shotdomain common-offset gathers are input to migration velocity analysis (MVA) and prestack inversion. ADCIGs are non-illusion prestack inversion gathers, and thus, accurate. We studied the extraction of elastic-wave ADCIGs based on amplitude-preserving elastic-wave reversetime migration for calculating the incidence angle of P-and S-waves at each image point and for different source locations. The P-and S-waves share the same incident angle, namely the incident angle of the source P-waves. The angle of incidence of the source P-wavefield was the difference between the source P-wave propagation angle and the reflector dips. The propagation angle of the source P-waves was obtained from the polarization vector of the decomposed P-waves. The reflectors’ normal direction angle was obtained using the complex wavenumber of the stacked reverse-time migration (RTM) images. The ADCIGs of P-and S-waves were obtained by rearranging the common-shot migration gathers based on the incident angle. We used a horizontally layered model, the graben medium model, and part of the Marmousi-II elastic model and field data to test the proposed algorithm. The results suggested that the proposed method can efficiently extract the P-and S-wave ADCIGs of the elastic-wave reverse-time migration, the P-and S-wave incident angle, and the angle-gather amplitude fidelity, and improve the MVA and prestack inversion. 相似文献