收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

基于MPI的大地电磁三维正反演并行算法研究

李焱  
【摘要】:大地电磁测深是利用天然交变电磁场为场源研究地球电性结构的一种地球物理方法,是深部地球物理探测必不可少的重要手段,在资源及矿产勘查方面有广泛的应用。当前,大地电磁数据主要使用一二维反演解释,但是实际的地电断面电性分布往往很复杂,难以满足一二维条件的假设。对这样的实测资料做一二维反演处理时,得不到可靠的结果。为了解决这种问题,需要进行三维大地电磁反演。三维大地电磁与一维和二维大地电磁最大的区别在于对地电模型假设的不同,三维大地电磁认为地下电阻率沿横向、纵向、垂向都有变化,这样就完全符合地下空间的真实状态。三维大地电磁反演研究的难点在于正演计算时间长、反演雅克比偏导数矩阵计算时间长、反演的非唯一性严重。因此,三维大地电磁反演的一个主要问题是计算量太大,特别是频点数和测点数较多、模型网格剖分较密、采用全张量阻抗元素计算灵敏度矩阵时,普通微机是很难承受的,运算速度更是让人难以接受,这也是目前制约三维大地电磁反演推广使用的最大障碍,并行计算的出现使全三维反演真正实用化成为可能。并行计算是指同时使用多种计算资源解决计算问题的过程,它能快速解决大型且复杂的计算问题,克服单个计算机上存在的存储器限制,既可以提高计算速度,义可以节省内存开支。因此,将并行处理技术引入三维大地电磁反演领域,为从根本上解决三维大地电磁反演计算的难题提供了一种有效途径。论文主要内容包括以下六个部分。 第一部分概述大地电磁反演的研究现状,介绍大地电磁一维、二维、三维反演方法以及并行计算在电磁领域的应用现状。 第二部分简要介绍MPI,包括MPI的特点、两种编程模式、通信方式以及MPI程序的基本设计。 第三部分介绍大地电磁三维正反演基本理论。三维正演使用的是交错网格有限差分数值模拟,详细给出其理论公式,包括边界条件、方程求解,以及张量阻抗、视电阻率相位的求取。三维反演使用的是数据空间算法,该算法是基于Occam方法的空间版本。通过数学推导将灵敏度矩阵从原先的模型空间形式(MxM)计算和存储转换到以数据空间形式(NxN)计算和存储,由于通常情况下模型参数量M大于甚至远大于数据量N,因此该方法可以大大降低运算量和所需要的内存。最后给出了数据空间反演方法的计算过程。 第四部分给出大地电磁三维正演并行算法,针对不同频率值对应的电磁响应是相互独立这一特点,我们采用分频并行思想。为了检验开发的三维正演并行算法正确性和有效性,我们设计了四个模型分别进行试算。详细给出了四个模型视电阻率相位的响应断面图以及频率对应的切片图,并对其分析。最后通过和串行算法对比,证明了我们开发的三维正演并行程序是正确的、高效的。 第五部分给出大地电磁三维反演并行算法。通过三维反演程序计算分析和测试,我们发现,三维反演计算时间长、需要存储多的主要是灵敏度矩阵、多次正演以及义积矩阵。针对这一特点,我们开发了大地电磁三维反演并行程序,并给出其计算流程图。为了验证开发的三维反演并行程序,我们对第四章中模型一和模型四分别进行试算,结果表明,开发的三维反演并行程序是正确的、高效的。另外我们还开发了三维正反演程序的接口可视化程序,通过直接操作该程序可以很方便的准备正反演需要的数据以及结果成图。 最后第六部分总结了本论文取得的主要成果,并指出若干不足及未来需要改进的地方。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 李焱;胡祥云;吴桂桔;叶益信;廖国忠;;基于MPI的二维大地电磁正演的并行计算[J];地震地质;2010年03期
2 L.B.Pedersen ,李然;大地电磁数据反演:非线性最小二乘法[J];地球物理学进展;1990年01期
3 徐世浙,刘斌;大地电磁一维连续介质反演的曲线对比法[J];地球物理学报;1995年05期
4 陈东敬,于鹏;大地电磁与地震联合反演技术在西亚某区的综合应用[J];石油物探;2005年03期
5 王书明,王家映,林长佑;大地电磁在地震预报中的应用[J];地质科技情报;2002年02期
6 魏胜,王家映;大地电磁解释工作站[J];地球科学-中国地质大学学报;1993年05期
7 林长佑,武玉霞,杨长福,陈军营;利用大地电磁测深法识别EDA地震电性变化前兆[J];西北地震学报;1995年04期
8 王连新,韩德波,党英;大地电磁测深法在大地构造研究中的一个应用实例[J];矿产与地质;1996年04期
9 周红,宋维琪,赵冲;一种快速二维大地电磁正演方法[J];石油物探;1998年02期
10 严良俊,胡文宝;大地电磁测深资料的二次函数逼近非线性反演[J];地球物理学报;2004年05期
11 王书明,王家映;高阶统计量对改进大地电磁测深资料处理方法的作用[J];石油地球物理勘探;2004年S1期
12 林昌洪;谭捍东;佟拓;;大地电磁全信息资料三维共轭梯度反演研究(英文)[J];Applied Geophysics;2011年01期
13 赵国泽,刘国栋,王家映;北京—柏各庄剖面连续介质大地电磁反演解释[J];地震学报;1986年02期
14 苏鸿尧,何展翔;表层电性不均匀对大地电磁测深曲线的畸变研究[J];地质科技情报;2000年03期
15 邓明,杜刚,张启升,余平,邓靖武;海洋大地电磁场的特征与测量技术[J];仪器仪表学报;2004年06期
16 张全胜,王家映;大地电磁测深资料的去噪方法[J];石油地球物理勘探;2004年S1期
17 于鹏,王家林,吴健生;大地电磁全频段偏移成像结合反演技术[J];同济大学学报(自然科学版);2005年04期
18 陈乐寿;大地电磁方法——能源勘探中的一种新手段[J];物探与化探;1984年04期
19 J.M.Travassos;D.Beamish;冰清;;大地电磁数据处理——分析一例[J];地球物理学进展;1988年04期
20 李予国,徐世浙,刘斌,于涛;电导率分块连续变化的二维MT有限元模拟(Ⅱ)[J];高校地质学报;1996年04期
中国重要会议论文全文数据库 前10条
1 李焱;胡祥云;吴桂桔;廖国忠;;基于MPI的三维大地电磁正反演的并行算法研究[A];中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集[C];2010年
2 刘鹏茂;柳建新;刘文劼;杨杰;周安发;;基于MPI大地电磁二维正则化反演并行算法研究[A];中国地球物理学会第二十七届年会论文集[C];2011年
3 徐光晶;黄清华;李墩柱;汤吉;;天然电磁干扰对大地电磁反演的可能影响[A];中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集[C];2010年
4 管然浩;景建恩;魏文博;韩杰;;大地电磁数据常见噪声及处理方法研究[A];中国地球物理学会第二十七届年会论文集[C];2011年
5 陈小斌;;地形对大地电磁观测资料的影响及二维带地形反演校正[A];中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集[C];2010年
6 王立凤;詹艳;赵国泽;汤吉;陈小斌;王继军;肖骑彬;;汶川强余震大地电磁观测试验和有关问题的思考[A];2011中国地震学会地震电磁学专业委员会换届暨学术研讨会摘要集[C];2011年
7 鲍光淑;鲁来玉;张碧星;;三维大地电磁模拟[A];2000年中国地球物理学会年刊——中国地球物理学会第十六届年会论文集[C];2000年
8 陈小斌;叶涛;;大地电磁多重网格二维反演方法的应用实践[A];中国地球物理学会第二十七届年会论文集[C];2011年
9 于鹏;王家林;吴健生;;提高大地电磁曲线反演分辨率的新方法[A];中国地球物理.2003——中国地球物理学会第十九届年会论文集[C];2003年
10 张秉政;徐义贤;;大别山腹地大地电磁WAL分析及二维反演[A];中国地球物理学会第二十七届年会论文集[C];2011年
中国博士学位论文全文数据库 前10条
1 蔡剑华;基于Hilbert-Huang变换的大地电磁信号处理方法与应用研究[D];中南大学;2010年
2 彭淼;大地电磁与天然地震数据联合反演研究[D];中国地质大学(北京);2012年
3 刘云;起伏地形大地电磁、时间域瞬变电磁二维数值模拟及直接反演法[D];成都理工大学;2012年
4 张大海;3-D天然场源地球电磁场响应的分析研究[D];中国科学院研究生院(上海天文台);2003年
5 杨生;大地电磁测深法环境噪声抑制研究及其应用[D];中南大学;2004年
6 陈小斌;大地电磁正反演新算法研究及资料处理与解释的可视化集成系统开发[D];中国地震局地质研究所;2003年
7 张帆;基于MPI和GPU直流电法和大地电磁法三维正演的并行算法研究[D];中国地质大学(北京);2011年
8 邱殿明;黑龙江省东部岩石圈演化特征[D];吉林大学;2005年
9 田越;cCTA与MPI异机融合图像在冠状动脉慢性完全闭塞病变中的应用研究[D];第四军医大学;2010年
10 霍光谱;二维大地电磁各向异性研究[D];中国地质大学;2012年
中国硕士学位论文全文数据库 前10条
1 李焱;基于MPI的大地电磁三维正反演并行算法研究[D];中国地质大学;2011年
2 刘鹏茂;基于MPI大地电磁二维正则化反演并行算法研究[D];中南大学;2012年
3 任伟;三维大地电磁数据反演效果研究[D];中国地质大学(北京);2011年
4 刘谋荣;低噪声低功耗大地电磁数据采集电路的设计与实现[D];哈尔滨工程大学;2012年
5 徐志敏;庐枞大地电磁干扰噪声研究[D];中南大学;2012年
6 倪鹏;大地电磁环境噪声抑制研究及其在冷水坑矿集区中的应用[D];东华理工大学;2012年
7 南东辰;淮南—溧阳地质廊带大地电磁数据处理与反演解释研究[D];中南大学;2012年
8 荣莹;基于MPI和Linux的集群系统的构建和应用[D];成都理工大学;2005年
9 汪睿;KD60平台MPI通信库优化设计[D];中国科学技术大学;2011年
10 赵少林;广域网协议在HPC集群上的应用分析[D];郑州大学;2005年
中国重要报纸全文数据库 前10条
1 张丽华 章力 (本报记者 张丽华 );“入地”先锋:点亮我们脚下的“黑狱”[N];中国矿业报;2010年
2 记者 廖文根;世界屋脊在向东漂移[N];人民日报;2005年
3 郭韶华 耿文明;福建局一九六队引进V8系统[N];中煤地质报;2007年
4 王俐 编译;越南加强吸引外资力度[N];国际商报;2002年
5 康维海;青海局将投资2000万元引进先进设备[N];中国矿业报;2007年
6 张素红;晋山区地下水勘测技术难题破解[N];中国矿业报;2006年
7 康维海;省地矿局投资2000万购置先进勘探设备[N];青海日报;2007年
8 丁丁;德国邮政面面观[N];中国邮政报;2001年
9 记者 史延廷;迪拜卓美亚集团与MPI签署全球合作协议[N];中国旅游报;2010年
10 刘浩;韩国服装业发展速度渐缓[N];国际经贸消息;2002年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978