共查询到20条相似文献，搜索用时 31 毫秒
1.
《地球物理学进展》2017,(5)
本文将二维时间空间域和频率空间域声波全波形速度反演方法分别应用到Marmousi模型,进行数值试验.两种方法均采用相同的观测系统和其他的参数,理论模型的数值试验结果证实了:使用较多的计算集群的CPU进行二维频率空间域直接法声波全波形反演时,其加速有限(正演数值模拟的计算量主要用于稀疏矩阵的LU分解,炮点计算波场时为线性关系).二维时间空间域声波全波形反演计算时更灵活,多炮同时计算时,可以多倍提高其计算效率;二维声波全波形速度反演时,直接法求解频率空间域的计算速度远快于时间空间域,所需要的计算机内存也比时间空间域少.二维声波全波形速度反演时,相比较于时间空间域的方法,频率空间域直接法声波全波形反演具有计算速度快和节省计算机内存需求的优势. 相似文献
2.
3.
起伏地形和关断时间对地面瞬变电磁响应影响严重,这给传统基于水平地表模型和理论阶跃波形的瞬变电磁数据解释技术带来很大困难.为此,本文开展考虑起伏地形和关断时间的地面瞬变电磁三维反演算法研究.正演采用基于非结构网格和后退欧拉隐式时间离散格式的时间域有限元算法,快速模拟起伏地表模型瞬变电磁响应.反演采用L-BFGS算法,减少每次反演迭代的计算量.针对发射线圈随地表起伏变化的特点,利用基于偶极子离散的场源处理技术模拟发射源的实际形状,采用瞬时电流脉冲技术实现考虑关断时间的地面瞬变电磁三维正演模拟.我们首先将本文开发的三维反演算法应用于理论模型的反演计算中,检验本文算法的可靠性,并分析地形和关断时间对反演结果的影响特征.在此基础上,进一步将本文算法应用于实测数据反演,验证本文算法的实用性. 相似文献
4.
《地球物理学进展》2017,(2)
<正>演是反演的基础,有效正确的正演差分格式可以保证反演结果的精度和效率.本文通过声波波动方程域间转换提出标量地震波Laplace-Fourier域数值模拟方法,并推导了同时引入衰减因子和频率的Laplace-Fourier域标量波方程的9点法有限差分格式和Laplace-Fourier域对应的加入PML(perfectly matched layer)吸收边界条件的差分格式,并通过模型试算验证了本文提出的Laplace-Fourier域正演方法的有效性和准确性,通过与时间域正演方法得到的地震记录比较,可以看出该方法能满足正演数值模拟的要求,为下一步进行Laplace-Fourier域标量波全波形反演奠定了基础. 相似文献
5.
《地球物理学进展》2017,(6)
我们采用区域分解(物理上分割模型,使用基于MPI的分布式存储架构的计算集群,从而节约单个CPU内核的内存使用量,加快正演数值模拟的计算速度)和炮并行(能够加快计算速度)的双并行算法,进行L-BFGS算法的二维时间-空间域声波全波形速度反演.我们采用多尺度的策略,只需要使用三个离散频率(5 Hz,8 Hz,12 Hz),从数据的低频成分开始反演,将低频的反演结果作为高频反演时的初始速度模型,依次反演数据的高频成分,进行Marmousi理论模型的全波形反演数值试验.数值试验反演所恢复得到的速度证实了:二维时间-空间域声波全波形速度反演方法计算灵活,可以适用于任何的采集观测系统,对地震数据可以方便地加时窗;使用多个计算节点同时计算多炮时,能够多倍提高数值计算效率.二维时间-空间域声波全波形速度反演所恢复得到的速度模型的分辨率较高. 相似文献
6.
研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法(BICG)计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次“拟正演”问题,进而将PARDISO的快速正演方案应用于“拟正演”问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性. 相似文献
7.
8.
9.
从各向异性弹性波的有限元正演方程出发 ,导出了反问题中时间域雅可比矩阵求解的计算公式 .它具有与时间域有限元正演方程相同的表达形式 ,故可通过有限元正演计算来获得雅可比矩阵 .研究了有限元正演算法的效率和精度、吸收边界条件等方面的问题 ,以提高反演系统的效率和精度 .在此基础上 ,实现了叠前全波场各向异性弹性参数反演 .计算表明 ,在初始模型偏离真实模型较大的情况下 ,层状模型和横向不均匀模型的反演结果均能准确地收敛到真实模型上 . 相似文献
10.
为推进大地电磁三维反演的实用化,本文实现了基于L-BFGS算法的带地形大地电磁三维反演.首先推导了大地电磁法三维反演的Tikhonov正则化目标函数以及Hessian矩阵逆矩阵近似表达式和计算方法,然后设计了一种既能保证空气电阻率固定不变又能保证模型平滑约束的协方差矩阵统一表达式,解决带地形反演问题.在反演算法中采用正则化因子冷却法以及基于Wolf条件的步长搜索策略,提升了反演的稳定性.利用开发的算法对多个带地形地电模型（山峰地形下的单个异常模型、峰-谷地形下的棋盘模型）的合成数据进行了三维反演,并与已有大地电磁三维反演程序（ModEM）进行对比,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了华南某山区大地电磁实测数据,得到该区三维电性结构,揭示了研究区以高阻介质为基底,中间以低阻不整合面和相对低阻介质连续分布,浅部覆盖高阻介质的电性结构特征,进一步验证了本文算法的实用性. 相似文献
11.
在利用地震波数据进行地球物理反演时,地层对地震波的吸收衰减效应会对地层物性参数的准确反演产生较大的影响,因此利用黏弹性声波方程进行反演更符合实际情形.本文在考虑地层衰减效应进行频率空间域正演模拟的基础上,提出基于黏弹性声波方程的频率域逆散射反演算法并对地震波传播速度进行反演重建,在反演过程中分别用地震波传播复速度和实速度来表征是否考虑地层吸收衰减效应.基于反演参数总变差的正则化处理使反演更加稳定,在反演中将低频反演速度模型作为高频反演的背景模型进行逐频反演,由于单频反演过程中背景模型保持不变,故该方法不需要在每次迭代中重新构造正演算子,具有较高的反演效率;此外本文在反演过程中采用了基于MPI的并行计算策略,进一步提高了反演计算的效率.在二维算例中分别对是否考虑地层吸收衰减效应进行了地震波速度反演,反演结果表明考虑衰减效应可以得到与真实模型更加接近的速度分布结果,相反则无法得到正确的地震波速度重建结果.本文算法对复杂地质模型中浅层可以反演得到分辨率较高的速度模型,为其他地震数据处理提供比较准确的速度信息,在地层深部由于地震波能量衰减导致反演分辨率不太理想. 相似文献
12.
13.
为更好地处理与解释复杂海底地形条件下测得的海洋可控源电磁数据,本文提出了一种基于非结构网格剖分的频率域海洋可控源电磁数据三维正则化反演方法.该方法首先对海洋地电模型以非结构四面体单元进行离散,然后基于矢量有限元方法获得海洋可控源电磁响应和灵敏度信息,最后采用共轭梯度法求解高斯-牛顿反演方程计算模型修正量.为提高反演的稳定性,通过在反演过程中采用对数转换方法实现反演模型参数的上下限约束.本文分别测试了单测线水平海底地形反演算例和面积性测量的起伏海底地形反演算例.反演结果表明,本文提出的频率域海洋可控源电磁三维反演能够准确地恢复高阻储油层的位置和电阻率信息,且计算效率较高,可用于实测海洋电磁资料的处理与解释. 相似文献
14.
采用新近提出的多次波射线追踪正演算法,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,分析讨论了利用多震相走时资料进行联合反演成像的方法及技术.考虑到不同震相走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数; 由于同时反演速度模型和反射界面起伏中不同模型参数变化对走时影响程度的不同, Jacobi偏导矩阵元素中引入了不同参数的归一化因子; 另外,为了克服射线密度过大(或过小)区域速度模型的过度(或欠)更新问题,反演算法中引入了等权射线密度的概念.几种数值模拟实例表明(含噪声敏感性试验): 多波走时的联合或同时反演成像技术是一种提高走时成像空间分辨率,进而降低重建模型失真度行之有效的方法. 相似文献
15.
《应用地球物理》2015,(4)
弹性波场的复杂性使得反演问题非线性增强,容易陷入局部极值,需要采用合理的多尺度反演策略降低非线性。在逐频组多尺度反演的基础上引入第二级别的反演策略,即基于阻尼波场的层剥离方法,能够改善反演过程的稳定性。针对频域全波形反演计算效率低、内存占用大的问题,采用双级并行算法:(i)利用多波前大规模并行直接解法(MUMPS)软件包,多节点并行实现波场正演;(ii)基于MPI实现频组内各频率并行计算梯度和步长等,使得多尺度反演算法在提高精度的前提下,保证了计算效率,提高了算法的实用性。Overthrust模型的数值实验表明,本文反演算法能够在有效改善反演稳定性的前提下,高效地获得精度较高的反演结果。 相似文献
16.
频率域数值模拟是频率域全波形反演的基础,在地震波场数值模拟中占有重要地位.相对于时间域数值模拟,频率域数值模拟具有两个明显的优势:没有时间累计误差,适合于并行计算.然而,严重的数值频散和巨大的内存损耗是阻碍其应用的两大瓶颈.为解决这两个问题,基于有限差分方法,学者提出了多种差分格式,如优化9点、15点、17点以及25点差分格式.本文从频散关系、计算效率和存储量三个方面,对比、分析了以上四种差分方法.基于2D声波方程,通过在均匀模型、层状模型以及Marmousi模型上的应用效果,对每种方法的优缺点进行了总结,为高精度数值模拟和声波频率域全波形反演提供方法选择上的参考. 相似文献
17.
波动方程数值模拟是深入研究地震波传播规律的有效方法.有限差分法因其方法简单、精度高而得到了广泛的应用.但其缺点是不能准确模拟具有复杂几何形态的物性界面.因而当遇到起伏地表或复杂构造时,求解精度低.为了准确模拟起伏地形、复杂构造和复杂介质条件下的地震波场,本文采用有限元法模拟二维声波方程.用三角形单元模拟地形和速度界面;把单元内的场和波速均看作单元上的线性函数,以适应复杂介质压制边角散射;采用吸收边界条件去除来自截断边界上的反射;采用集中质量矩阵和集中阻尼矩阵使得显式时间递推无需对矩阵求逆,提高了计算效率.对模型的计算表明该方法正确有效. 相似文献
18.
频率域全波形反演是重要的地震成像方法,而频率域波动方程数值模拟是频率域全波形反演的基础.对于大规模的问题,由于受存储和计算量的限制,基于LU分解的直接方法一般不再适用,而是采用迭代方法.基于多重网格预条件的双共轭梯度稳定化方法是一种重要的迭代方法.本文重点讨论了多重网格预条件求解过程中的松弛因子选择方法,研究结果表明,(1)对于一般选取的松弛因子,随模型复杂性的增加,所能计算的重数逐渐下降,方法的实用性也随之下降;(2)对于复杂模型,采用局部模式分析方法选取松弛因子,提高了所能计算的重数,保证了多重网格方法的收敛性和实用性.这些研究成果对基于多重网格预条件的迭代算法的实际应用具有重要意义. 相似文献
19.
层状介质的声波波动方程反演 总被引：4，自引：3，他引：1
基于广义反射透射系数矩阵正演方法 ,讨论了层状介质的声波波动方程反问题 .推导出波数频率域中的雅可比矩阵的解析表达式 ,其计算在正演过程中求出 .采用最小二方法可得到层介质参数 .数值结果表明反演方法的正确有效性 . 相似文献
20.
全波形反演已被广泛应用于获取地下速度结构.而反演问题与正演方法密切联系,针对特定反演问题,合适的正演方法能极大提高反演效率和精度.本文首先验证谱元法在含起伏界面模型数值模拟方面的优势,在此基础上将谱元法作为正演引擎应用于全波形反演,并为克服未知子波的影响,采用一种归一化的频率域目标泛函.结果表明,起伏地表情况下,基于谱元法的全波形反演相比于基于传统有限差分法反演,具有更高的反演精度.进一步,本文将基于谱元法的波形反演方法应用于OBS观测系统的理论合成数据和野外采集数据.谱元法非结构化网格剖分自然满足自由边界条件,能很好地适应不规则海床并模拟多次波.理论实验表明,即使在OBS观测系统很稀疏的情况下,基于谱元法的全波形反演仍能获得海底以下正确的高波数速度结构.在处理实际OBS数据时,本文采用分频策略以减少反演非线性,初始模型成功更新,其结果揭示了西沙海槽海底以下更多的细节信息. 相似文献