共查询到20条相似文献，搜索用时 515 毫秒
1.
叠前逆时偏移将是应用前景最为广泛的地震偏移技术,但存在计算量大和存储单元要求多的瓶颈问题.随机边界条件思想的提出较好地解决了其中的存储问题,但同时也大大增加了计算成本.为了减小边界宽度,降低宽边界带来的计算负担,同时又不影响随机边界的随机效果,本文首先提出了一种多形式随机边界函数,并从随机介质散射理论出发,研究了随机边界的散射特性,并将随机边界归纳为增强模式和衰减模式两类,深入分析推导了不同模式边界的随机函数形式选择的最优判据,最后形成了一套能够使得在边界宽度有限情况下获得最佳散射效果的随机边界构建算法.理论分析和数值模拟实验结果表明:在边界宽度较窄的情况下,利用本文提出的随机边界构建策略,能够最大程度地发挥有限宽随机边界散射性能,压制相干能量的产生.这为提高随机边界的随机性能,减小随机边界的宽度,提高逆时偏移的计算效率提供了一种有效的方法. 相似文献
2.
《地球物理学进展》2017,(6)
TTI介质逆时偏移是复杂构造区地震成像的一种重要有效技术.然而,其3D数据实际应用仍然存在一些挑战.本文以TTI拟声波方程为基础,系统梳理了TTI逆时偏移计算中的几个关键问题并给出解决方案.首先,P波与SV波耦合问题.本文系统分析了各向异性参数匹配方法伪横波压制效果以及其对波传播的影响;其次,为了提高计算效率,本文基于CUDA平台实现GPU加速计算方案.主要解决了两个方面的问题,其一为扩展了随机边界在TTI介质波场传播中的应用,确定了TTI介质速度和其他四个参数对应的"均匀各向异性"随机边界;其二,本文将单GPU拓展到多GPU实现三维TTI介质逆时偏移,并且采用了一种异步流技术,大大提高了GPU计算的高并发策略,使双GPU实现的加速比与单GPU相比接近理想的2∶1.最后通过3D模型试验证明了方法的效率和准确性. 相似文献
3.
基于双程波动方程的逆时偏移被认为是目前最好的偏移成像技术,更适合于复杂构造成像.然而,大计算量和大存储量使得逆时偏移的计算成本很高而无法用于大数据量的地震成像.本文分析了目前常用的存储策略,并分别在空间和时间上对存储策略进行了研究:空间上,根据有限差分格式,在边界存储策略的基础上通过修改波场逆向传播的边界条件,提出了有效边界存储策略.该策略可在不增加任何计算量的情况下大幅降低逆时偏移对存储量的需求;在时间上,使用checkpointing技术对有效边界存储策略进行了改进,使叠前逆时偏移在增加少量计算量的情况下进一步降低存储量需求.Marmousi模型测试结果表明了有效边界存储策略的有效性和优越性. 相似文献
4.
《地球物理学进展》2015,(4)
基于双程波动方程的叠前逆时偏移方法成像精度高,而且无地层倾角限制,较适合于复杂地下构造成像.但是,由于地下介质的各向异性广泛存在,基于各向同性的正演算法,尚难准确描述真实的地下波场传播,逆时偏移的成像精度也因此受到限制.鉴于此,本文研究了各向异性VTI介质逆时偏移方法,首先根据VTI介质一阶准P波方程推导出了炮点和检波点的逆时延拓的交错网格高阶差分格式,针对算法计算量和存储量大的问题,文中研究了一种改进的基于GPU加速的有效PML边界存储策略.本文建议的方法只需增加少量的额外计算,就可降低大幅度的存储成本,进而实现高精度和高效率的各向异性逆时偏移.Hess 2DVTI模型测试表明,本文提出的方法不需要存储全部历史时刻的波场,可以实现高效率高精度的VTI介质叠前逆时偏移成像. 相似文献
5.
在地震资料的处理中,逆时偏移方法可以实现复杂构造高精度成像,但计算量和存储量两大问题影响了该方法的实际应用.逆时偏移算法中的边界处理方式的优化可以大大减少存储量,本文讨论了记录波场边界信息的吸收边界和随机边界这两种方法的计算效率和成像效果,提出了吸收边界中只记录边界范围内未衰减单层波场的方法并通过数值实验验证了其可行性,选择了不同的边界策略应用于模型数据和实际资料的处理.结果表明:记录单层波场边界信息的吸收边界的成像效果同传统存储波场历史的方式几乎无差别,但要额外存储每个时刻的波场边界信息;随机边界不需要额外存储波场信息,但会带来边界漫反射影响和计算区域的增加;记录单层波场边界可以明显减少存储量,并且不影响成像效果. 相似文献
6.
地震波走时层析是主要的近地表速度建模技术.对于变速模型,尤其是复杂近地表模型,波动方程能够更加精确地描述地震波的一阶绕射效应,理论上反演精度高于射线(束)类层析.然而,波动方程的数值求解的计算量较射线理论增加了几个数量级,严重制约了波动方程走时层析技术的三维实用化.为解决波动方程走时层析中遇到的计算和存储问题,本文通过引入随机边界条件以及离散Fourier积分,提出了基于单频梯度反演策略的实用化三维波动方程初至层析方法.本文方法的优势在于：在计算和存储方面,相比于基于波场重构算法(需要三次波传播)的梯度计算策略,本文提出的单频梯度计算方法对内存需求极低(相对于波动方程正演,仅需要增加两个单频波场以及一个成像体),且计算量至少减少1/3(仅需要两次波传播且无需处理数值吸收边界).此外,随机边界的引入使得正演算法大幅简化因而更适用于GPU(Graphics Processing Unit,图形处理器)加速.在理论层面,由于影响地震波一阶散射效应的主要区域集中在连接炮检的中心射线邻域的第一菲涅尔带内,随机边界条件的引入降低了随机边界内产生的散射波的空间相关性,使得第一菲涅尔带内的单频梯度相干... 相似文献
7.
叠前逆时偏移是当前最为准确的地震成像方法,由于计算量大、存储量大等原因需要合适的实现策略和高效的计算平台.本文以高阶有限差分逆时偏移为基础,重点讨论了在GPU上实现需要解决的显存不足问题和人工边界问题.利用区域分解技术可以在当前GPU上高效地实现任意生产规模的三维逆时偏移成像,不会受到GPU显存规模的制约.常规最佳匹配层边界条件边界区域控制方程与内部区域差异较大,不适于GPU高速运算.本文在GPU上实现近似最佳匹配层(NPML)边界条件,使得高阶有限差分计算不需要分支判断,边界区域辅助波场的存储量也较低,保证了在GPU上进行波场传播的高效性.三维理论数据和实际资料成像结果表明了本文方法的正确性. 相似文献
8.
地震叠前逆时偏移是当前公认的地震成像的有效途径,然而它面临着计算量甚巨,低频成像噪音以及存储量大等问题,因此,业内科研工作者对其研究乐此不疲.借助GPU/CPU协同计算可以有效解决计算量的难点,笔者已在另文中阐述,本文着重探讨成像噪音抑制以及存储问题.文中分析了叠前逆时偏移产生成像噪音的机制,据此提出在叠前地震资料中先对数据进行相位与振幅校正,进而在成像后运用拉普拉斯算子滤波法消除成像噪音,从而有效去除成像所产生的低频噪音;针对存储量,采用随机边界,用计算换存储,并借助GPU实现,节省了GPU与CPU之间的数据通讯,数值实验结果表明,采用随机边界方法的逆时偏移结果与直接存储波场的方法得到的结果差别甚小. 相似文献
9.
《应用地球物理》2019,(4)
本文从反射声波测井逆时偏移成像方法出发,考虑逆时偏移本身存在海量波场数据存储问题,探讨了波场存储优化策略,提出并实现了激发振幅成像条件逆时偏移和随机边界法逆时偏移两种避免波场数据存储的优化方法。首先,激发振幅成像条件方法利用计算激发时间作为成像时间判断,只需存储小部分波场,如几十个时间点的波场数据,通过实例证明甚至仅需两个时间点波场即可完成成像。传统逆时偏移通常需存储数千时间点的波场数据,与之相比数据存储和读写量可缩小几十倍甚至数千倍;其次,随机边界方法利用波场在随机介质中发生散射而非反射进行震源声场反传重构,从而避免正传波场数据的存储。数值计算表明两种方法均可有效减少反射声波逆时偏移的波场数据存储,与其它偏移算法和常规逆时偏移相比,在保证成像精度的同时提高了数据处理效率,为井旁缝洞体高效、高精度成像提供了重要的技术手段。 相似文献
10.
地震数据TTI介质逆时偏移计算量巨大,阻碍了其深度应用。NVIDIA? GPU以及其通用计算具有高度并行特性,为加速该类高密度计算提速提供了可能,可解决由于地震处理成像空间增加所带来的计算量大难题。本文提出了一种高效率的多GPU并行计算策略,以解决TTI介质逆时偏移庞大的数据处理问题。该策略流程,以GPU及其统一设备架构CUDA为前提,拥有多核心同时计算、作业流、点对点GPU直接传输等一些列运算特点。其核心是将GPU间边界数据的数据传输时间和不同区域的差分计算时间重叠。由于逆时偏移计算强度主要与差分计算空间相关。重叠计算后,数据传输时间可以忽略,因而计算效率随GPU数量的增加呈线性提高。并用于TTI逆时偏移成像的处理,以验证本文提出的高计算效率的多GPU并行计算策略的正确性。对比试验表明,利用本文提出的多GPU策略可大大提供高密度数据成像计算的效率及实现多GPU计算时效率呈线性增加,提升了该计算的延展特性。 相似文献
11.
一种改进的二阶弹性波动方程的最佳匹配层吸收边界条件 总被引：3，自引：0，他引：3
在分析前人有关二阶波动方程的最佳匹配层(PML)吸收边界条件的构建方式的基础上,讨论了这些PML吸收边界条件在计算效率和计算精度方面的不足,并给出了一种新的改进PML吸收边界条件.本文的核心思想是在频率域引入一个中间变量,避免了褶积运算,对该变量做傅里叶逆变换得到时间域的精确解.通过模型试算,把本文的改进算法与前人的算法的计算精度和计算时间进行比较.对比分析表明,本文提出的PML吸收边界条件计算量小且精度较高,并且是一种稳定的算法. 相似文献
12.
三维交错网格有限差分地震波模拟的GPU集群实现 总被引：1，自引：1，他引：0
有限差分实现简单、速度快,作为地震波场模拟一种有效数值方法,被广泛用于正演计算密集的波形反演和逆时偏移中.三维地震波正演模拟计算量大,一直以来制约着三维叠前逆时偏移和反演的工业化应用,GPU通用计算技术的产生及其内在的数据并行性有望改变这一现状.本文通过分析三维交错网格有限差分方法在GPU上的实施,利用片内共享存储器实现了三维地震波数值模拟的高效算法,取得了较单核CPU快79x～108x的加速比;通过区域分解技术将单GPU上不能计算的地质体模型沿Z轴方向进行粗粒度分解,采用消息传递接口交换边界数据,运用MPI+CUDA的方式实现了大尺度三维地震波场模拟,并着重分析了影响GPU并行计算效率的一些关键因素.大尺度三维地震波场模拟的加速实现,为促进叠前逆时偏移和波形反演技术的工业化转化提供了可能,因此具有重要的研究意义. 相似文献
13.
叠前逆时偏移等基于波场互相关原理的地球物理方法存在极大的计算与存储需求,因此采用合适的波场重构方法显得尤为重要.常规的随机边界法容易产生成像噪声,而有效边界法在三维情况仍难以实现,检查点技术具有内存要求小的特点,但存在较高的重算率,因此本文提出了插值原理的检查点技术波场重构方法.在满足Nyquist采样定理的前提下对相邻检查点间的波场进行规则抽样,将抽样波场作为插值节点,运用多项式插值算法重构任意时刻的波场,从而避免优化检查点技术反复递推造成的计算效率问题.数值实验表明：插值检查点重构算法能有效的恢复波场,其中三次样条插值重构精度最高,而牛顿法插值法计算代价较小适合于快速重构.经Sigsbee模型的叠前逆时偏移证明了插值算法的可行性,并且极大的提高了波场重构的计算效率.三维模型分析得出在增加少量存储的情况下插值重构法的重算率大幅度降低,存储量减少为有效边界法的7.1%,对于三维尺度的叠前逆时偏移有实际意义. 相似文献
14.
相较于单程波偏移算法而言,逆时偏移成像方法以其物理基础为依托优势,几十年来一直备受国内外地球物理学家的青睐.目前的逆时偏移(RTM)若直接采用双程波动方程进行延拓,尽管可以回避上下行波的分离处理,然就已有算法而言,其计算量和I/O(输入/输出)量却是最大的.针对此问题,本文在分析现行逆时偏移的多种算法基础上,提出利用CPU/GPU(中央处理器/图形处理器)作为数值计算核心,建立随机边界模型,从而克服存储I/O难题和提高计算效率.在实际的数据测试中,本文的方法可以大幅度的提高计算效率和减少存储单元,从而促使其高效地应用于生产实际. 相似文献
15.
《地球物理学进展》2017,(6)
地震资料逆时偏移方法在计算过程中需要保存大量的波场信息,这种巨大的存储需求限制了该方法的实际应用.因此,研究能够显著降低波场存储量的存储策略对于推进逆时偏移方法的实用化具有重要的理论意义和实际价值.本文分别从时间和空间的角度提出了相应的波场存储策略:(1)提出变阶数有限差分算子,在计算波场左、右和下边界处时,将差分阶数由高阶递减至二阶,结合有效边界存储策略,可在不损失波场重建精度的同时大幅降低波场存储量;(2)借鉴checkpoint策略和缓冲区思想,提出了基于拉格朗日乘子法的参数最优化方案,在时间方向上压缩波场存储量.通过对二维盐丘模型进行震源波场重建和逆时偏移存储量测试,结果表明:本文方法可实现对震源波场的高精度重建,与传统方法相比,在不损失成像精度的前提下,新的存储策略可大幅降低逆时偏移方法的波场存储需求. 相似文献
16.
传统有限差分系数是通过泰勒级数展开求取的,这样导致所计算的频散曲线在大波数区域会产生较强的数值误差.针对二阶空间偏导数的显式有限差分离散,本文发展了一种新的优化差分系数方法：首先将泰勒级数展开与多点采样方法结合应用于空间频散关系,基于最大范数建立直观有效的优化目标函数,采用Remez算法求解该目标函数,从而获得最优化差分系数.利用优化有限差分方法求解三维垂直对称轴横向各向同性（VTI）介质中的声波和弹性波方程.另外,本文将二维混合吸收边界条件推广到三维VTI介质中,用于吸收人工截断边界反射;基于各向异性特征,合理调整了边界区域的速度值来提高吸收效果.考虑到三维情况下计算效率的问题,本文波场外推过程中采用图形处理器（GPU）取代传统的中央处理器（CPU）.数值精度分析表明,相比较于传统的泰勒级数展开方法,优化有限差分方法在大波数区域对频散误差的压制效果更明显.在三维均匀和修改的Hess VTI模型中的数值模拟实验证明了本文方法具有更高的精度与效率,混合吸收边界条件在三维VTI介质中具有良好的边界吸收效果. 相似文献
17.
叠前地震逆时偏移是目前公认地震成像的有效途径.声学介质逆时偏移已经成熟,弹性波逆时偏移研究相对较少,本文研究了各向同性介质弹性波逆时偏移.传统逆时偏移由于存储量的巨大不便于实际应用.本文针对逆时偏移存储量大的缺点,采用随机边界条件,用计算换存储减少了存储量.首先从波动方程出发,推导了各向同性介质弹性波逆时偏移高阶有限差分算子;其次运用随机边界逆时偏移的流程进行逆时偏移;最后采用传统的拉普拉斯滤波的方法消除低频噪音.模型试算表明,采用随机边界得到逆时偏移结果与常规有限差分得到的逆时偏移结果差别较小,证明了该方法的正确性. 相似文献
18.
在地震波传播数值模拟的过程中,需要使用吸收边界条件从而达到衰减人为边界 反射的目的. 本文基于傍轴近似法提出了计算三维弹性波方程的吸收边界条件公式,表示了 各边界面、边棱和角点处波场所满足的单程波方程,并在三维弹性波数值模拟中进行了应用 . 理论模型及三维盐丘地质模型波场切片快照试算结果表明,该吸收边界条件可以有效地吸 收人为边界反射,适用于较大入射角情况,从而消除了边界有效波信息的干扰. 由于采用四 阶近似方程,在保证计算精度的前提下,该方法具有节省计算工作量和易于实现的特点. 相似文献
19.
相对于传统的逆时偏移,最小二乘逆时偏移具有更高的成像质量,这种改善是通过迭代反演来获得的,另外其精度与效率依赖于求解波动方程算法的精度与效率.本文给出了基于非结构化网格的最小二乘逆时偏移的方法,该方法能够充分结合最小二乘算法与非结构化网格精细刻画地下界面以及随速度自适应剖分的优点;并采用带补偿的拉普拉斯滤波算法,来消除梯度计算中的低频噪声,从而加速目标函数的收敛速度.通过简单倾斜模型以及复杂Marmousi模型测试,显示了该方法的有效性和潜力. 相似文献
20.
叠前逆时偏移技术是解决地震成像问题的有力工具,但由于计算量大、成像噪音以及存储量大等原因没有得到广泛的应用.本文给出了逆时偏移的实现过程,分析了高阶有限差分格式的稳定性与频散关系.针对叠前逆时偏移计算量大的问题,使用图形处理器(Graphic Processing Unit,简称GPU)实现算法加速,比传统的CPU计算速度提高了一个数量级.文中对理论模型进行了计算,并与单程波偏移方法做比较,结果表明:叠前逆时偏移有效突破了成像倾角限制,对垂直断层、盐丘空腔内幕等特殊构造成像效果均有显著提高.本文尚未涉及成像噪音去除以及存储量等问题,笔者将另文阐述. 相似文献