首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 601 毫秒
1.
首先以频点为并行粒度,完成了MT三维正演并行计算,然后将该算法引入共轭梯度法MT三维反演过程中的正演求数据的残差、反演方程右端项的求取、共轭梯度求解模型增量这三个计算环节中,使得MT三维反演中的主要计算量实现并行化。对编制的MT三维正、反演并行计算程序,在多核工作站上通过理论地电模型进行了正演试算和反演验证,并在串行和并行两种计算方式下进行了计算效率对比。  相似文献   

2.
孙慧秋  韩立国  徐杨杨  高晗  周岩  张盼 《世界地质》2015,34(4):1106-1112
估计地层品质因子Q对描述地下构造分布和油气预测都有着重要的意义。考虑到地震波在介质中传播时的黏滞吸收作用,笔者选择2D黏滞性声波方程全波形反演,将观测波场与正演模拟计算得到的波场构建一个目标函数,引入L--BFGS算法反演出地层品质因子。将该方法应用于异常体模型和Marmousi模型测试,结果均表明反演的模型与理论模型的相对误差较小,计算精度高,反演效果好。  相似文献   

3.
地震波场频率域正演是频率域全波形反演的基础。针对反演计算量巨大的问题,利用频率域二维波形正演算法中频率和炮点计算的独立性,开发出粗细粒度结合的MPI并行算法:粗粒度为频率并行,细粒度为单个频率解方程并行。实现方法是将正演频率分组放入不同的MPI通讯域内,每个通讯域内单个频率求解方程过程,采用基于MPI的MUMPS(多波前大规模并行稀疏直接解法器)软件包并行加速。模型测试结果表明:MPI多通讯域并行算法计算结果正确,计算效率显著提高,加速效果稳定。  相似文献   

4.
近年来2.5维复电阻率反演取得了一些理论和应用方面的成果,但计算过程复杂、计算量巨大,特别是正演计算过程中计算量大导致了反演速度较慢。这里在2.5维复电阻率反演计算基础上,将反演算法中的正演过程的排列和频率循环分割成多个规模较小的程序段,并发执行正演过程中的每个排列和频率循环的求解计算过程,并根据并行计算的节点数目或者计算核心数目,动态调整每个节点或者计算核心的计算量,从而提高了算法的计算效率和适应性。  相似文献   

5.
在地震勘探方面,地震数值模拟技术在地球物理反演及观测系统设计领域有着广泛地应用。这里主要提出并实现了利用C++AMP并行加速内存/计算密集型的三维弹性波的正演模拟,使用C++AMP编写的代码可以有效地在图形处理器(GPU)等支持数据并行的硬件上运行。在重点研究C++AMP三维弹性波数值模拟的并行算法实现的基础上,结合CPU串行及CPU多核并行进行加速对比,三维盐丘实验的结果说明,利用GPU的强大并行计算能力的C++AMP,可以极大地提高数值模拟的效率。此外,C++AMP在Windows平台具有良好的可移植性,同时兼容NVIDA和AMD的设备,极大地提升了开发效率。  相似文献   

6.
王观石  熊鹏  胡世丽  孟世明  龙平  谭谈 《岩土力学》2018,39(6):2175-2183
考虑节理质量对应力波传播的影响,运用波的位移位函数推导了谐波在厚黏弹性节理的透、反射系数计算公式,采用波形相关系数描述子波穿过黏弹性节理的波形变化,讨论具有一定厚度的黏弹性节理简化位移不连续模型的适用条件。设厚黏弹性节理模型和位移不连续模型的透射波波形相关系数为0.9时对应的节理厚度为临界厚度,岩体与节理的阻抗比对临界厚度的影响很小;临界厚度随子波中心频率增大呈负指数减小,入射角越大,临界厚度随中心频率减小得越慢。试验数据分析表明:当节理厚度为0.03 m时,采用位移不连续模型和厚黏弹性模型计算得到的节理力学参数非常接近,随节理厚度和子波中心频率增加,运用位移不连续模型的计算结果偏差越大,试验结果与理论分析是一致的。  相似文献   

7.
利用快速标量传递算法进行了一维层状模型的地震面波频散曲线正演计算,得到了基阶频散曲线;利用传统方法对一维层状模型进行了大地电磁正演计算,得到大地电磁的视电阻率曲线数据和相位数据。利用遗传算法对合成数据分别进行了单独反演和联合反演,得到拟合差分布图来分析两种数据的兼容性并研究联合反演的可行性。通过对比单独反演和联合反演的结果,可以看出地震面波频散曲线和大地电磁数据的联合反演结果较之单方法的反演结果要更加精确,收敛速度更快,两种数据的耦合情况较好,同时联合反演能够有效地限制多解性。  相似文献   

8.
利用余弦变换的高精度优势,采用并行计算来弥补余弦变换在计算速度方面的不足。深入分析原串行程序中的外层时间节点循环和内层余弦变换这种多层次循环结构的特点,给出了3种并行方案,在WINDOWS系统上使用FORTRAN和MPICH2相结合的开发工具实现了一维瞬变电磁并行正演程序。通过理论模型并行试算,与串行程序进行对比,验证了并行程序的可行性和正确性。3种并行方案均有效地提高了正演计算速度,在相同进程数目的情况下,外层循环时刻点之间的粗粒度并行效率最高,为高维的正反演并行提供了设计思路及研究基础。  相似文献   

9.
频率域正演是频率域全波形反演的基础,高效准确地完成频率域正演计算是目前该领域的一个热门问题。这里利用优化17点差分格式和PML完全匹配层构建波场,并使用LSQR算法求解波场,通过在简单层状模型和Overthrust模型上计算验证可知,LSQR算法可大幅提高频率域正演的计算速度和模拟精度,从而为此方向的研究提供一个可行的技术手段。  相似文献   

10.
周雪  于平  翁爱华  陈瑞鼎 《世界地质》2018,37(3):897-904
为提高重力及重力梯度数据的正演效率,笔者引入基于MPI (Message Passing Interface)和OpenMP (Open Multi-Processing)的并行计算,通过对比分析不同数据规模的网格数和模型体个数对并行效率和加速比的影响,得出随着正演数据规模的增加,并行效率和加速比均得到提高。同时对比了基于MPI和基于OpenMP的两种并行方式的性能,结果表明,重力及重力梯度数据正演的并行计算中MPI的提速能力优于OpenMP,且在较大规模数据的正演计算中基于MPI的并行效率优于基于OpenMP的并行效率。  相似文献   

11.
蔡伟  宋先海  袁士川  胡莹 《地球科学》2017,42(9):1608-1622
反演瑞雷波频散曲线能有效地获取横波速度和地层厚度,传统的多模式瑞雷波频散曲线反演需要正确的模式判别.然而,当地层中含有低速软弱夹层或高速硬夹层等复杂结构时,瑞雷波可能会出现"模式接吻"和"模式跳跃"等现象,这些现象极易造成模式误判,进而导致错误的反演结果;同时,传统的频散曲线反演方法需要进行求根运算,进而导致现有的瑞雷波非线性反演速度慢,运算时间长.鉴于此,对传统的Haskell-Thomson频散曲线正演模拟算法进行了改进,提出了一种新颖有效的目标函数.该目标函数直接利用实测频散曲线与迭代更新模型频散函数表面形状进行最佳拟合,无需将多模式频散数据归于特定的模式,可有效避免多模式瑞雷波频散曲线反演模式误识别;同时,该目标函数不需要求根运算,进而大大加快了非线性反演速度.基于粒子群优化算法,利用实际工作中经常遇到的3种典型理论地质模型和某一高速公路路基实测资料进行了理论模型试算和实例分析,检验了本文提出的瑞雷波多模式频散曲线反演新方法的有效性和实用性.   相似文献   

12.
全波形反演不仅利用相位和振幅信息,还利用波形的细节变化,具有刻画模型精确细节的能力.在对稀疏矩阵直接LU分解求解的基础上,采用梯度预处理方法对声波介质速度模型进行了反射波全波形反演.采用误差反向传播算法计算目标函数梯度以及伪Hessian矩阵的对角线元素来做梯度预处理.数值模型的实验结果表明,利用有效的频率段便能反演出分辨率较高的速度结构,用低频反演出的结果作为高频反演的初始模型,减少了解的非唯一性.二维高斯光滑初始模型提供了有利的低频信息,得到较好的反演结果.伪Hessian矩阵的预处理吸收了高斯牛顿法的二次收敛优势,在不增加计算量的前提下,加快收敛速度.   相似文献   

13.
全波形反演可以利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而,庞大的计算量和存储空间需求,限制了全波形反演的发展。在频率多尺度全波形反演中将L BFGS数值优化算法与同时激发震源技术相结合的方法来改善这一现状。首先,对Marmousi模型进行了速度反演:在计算过程中明显发现对计算机内存的占用减少,最终反演结果与实际Marmousi模型的拟合误差为0.095 9,较小;采用10个频带单炮震源正演384炮所需时间约为32 640 s,而采用同时激发震源(384炮)正演一次所需时间仅约为700 s。然后,基于高速楔形体模型进行了抗噪能力研究:原始含噪地震记录信噪比为11.147 3 dB;对反演得到的速度模型进行正演,其地震记录信噪比为22.251 8 dB。最后,基于逆冲断层模型进行了反演速度扰动能力研究,反演得到的最终模型很清晰,与具有速度扰动特性的实际模型非常接近,拟合误差仅为0.036 0。数值模拟试验结果表明:此方法反演精度高,内存开销较小,能够显著提高计算效率,并且具有良好的抗噪能力,能够反演出具有速度扰动特性的介质。  相似文献   

14.
在一般的地震波计算中,地球介质可以当作各向同性的完全弹性体来对待。而对地表土壤和精细观测,粘弹性介质模型比完全弹性模型更符合实际。可控震源是近年来发展起来的精确可控、安全有效的新型震源,但国内对可控震源传播动力学的研究较少。对可控震源激发的地震波在粘弹性介质中的传播,用有限差分方法进行了数值模拟,给出了Kelvin一阶粘弹性波动方程组、吸收边界条件和震源处理。提出了检波器接收信号加速度和传播距离的三次多项式指数关系式,比线性指数关系更准确。使用该公式可以预测采集系统有效接收距离,指导测线布置。   相似文献   

15.
混凝土是一种含有骨料、水泥浆、钢筋、波纹管包裹的预应力孔道及空洞、孔隙的多相非均匀复合材料,作为一种粘弹性介质,超声在混凝土中传播机理比较复杂,常规的弹性声波方程不能确切地模拟超声传播方式。这里从粘弹性介质材料性质出发,推导了基于Kelvin粘弹性介质模型的高阶有限差分方程,通过建立几种粘弹性混凝土介质数值模型,模拟了超声波在混凝土中的传播过程,分析了含空洞、裂隙等不同缺陷的混凝土介质中超声波的波形特性,讨论了缺陷部位波形的变化规律,探讨了空洞、裂隙等不同缺陷对超声传播机理的影响以及能量分布的差异,为准确判断识别混凝土内部不同结构的缺陷提供科学借鉴。  相似文献   

16.
简谐波在黏弹性介质中波动传播的应力场研究   总被引:1,自引:0,他引:1  
得到了简谐波在具有Kelvin模型黏弹性介质中的外行扩散波与内行入射波应力场表达式。将波动在介质中的传播分为外行扩散波与内行入射波两部分,先得到简谐波在弹性介质中波动传播时外行扩散波与内行入射波各自的应力场,结合弹性波理论与黏弹性理论,分别得到了波动在Kelvin模型黏弹性介质中外行扩散波与内行入射波两部分的应力场,并将外行扩散波与内行入射波两部分的应力场分别表示为外行扩散波与内行入射波位移与速度场的表达式,其结果有利于波动理论在黏弹性介质中的应用  相似文献   

17.
重磁异常数据三维人机联作模拟   总被引:5,自引:0,他引:5  
在研究三角形多面体模型重、磁异常三维正演和反演技术的基础上,吸取人机交互正演模拟的优点,实现了三角形多面体模型重、磁异常数据三维人机联作反演。通过研究三角形多面体模型节点偏导数的计算方法,对目标函数进行线性化处理,形成了计算机自动迭代修改模型体的技术。利用计算机图形技术,在三维空间显示重、磁场和模型体,开发了模型的交互修改技术,使数据解释过程中,可以结合已知信息及人的推断和经验,完成重、磁异常数据的三维模拟,减少了数据解释结果的不确定性。  相似文献   

18.
张平  韩立国 《世界地质》2017,36(1):255-265
基于谱元法计算理论,通过坐标旋转算法实现任意各向异性黏弹性介质中的地震波场计算,研究复合介质中地震波场特征。根据计算结果,qSV波表现为剪切滞弹性形变,qP波表现为膨胀滞弹性形变;品质因子不仅影响地震波振幅,还影响其传播频率;黏弹性介质吸收高频快,在一定频率范围内,黏弹性介质没有频散现象。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号