收藏本站
《中国科学技术大学》 2010年
收藏 | 手机打开
二维码
手机客户端打开本文

地球电磁三维数值模拟的多重网格方法及其应用研究

鲁晶津  
【摘要】: 地球电磁勘探方法在能源及矿产资源勘探开发、水文工程地质勘察以及地壳上地幔深部研究中均发挥着重要作用。随着地质勘探目标日益复杂,三维精细结构探测成为迫切需求,导致复杂地下结构的地球物理三维数值模拟在计算效率上面临着越来越大的挑战。就目前地球电磁三维数值模拟的有限差分、有限元法而言,大规模非均匀网格、大的电性差异以及复杂的电性结构等因素,使得最后形成的大型线性方程组的系数矩阵条件数很差,严重影响着现有迭代解法的收敛速度(如不完全Cholesky共轭梯度法,即ICCG法)。多重网格(MG)方法为近年来计算数学领域发展的新的高效算法,在数学上被证明是求解线性椭圆型偏微分方程的最优解法,其收敛速度与网格节点数无关。但由于诸多困难,在地球电磁三维数值模拟中的应用尚极少见。 论文首先针对地球物理中常见的Poisson方程及Helmholtz方程进行了几何多重网格法数值模拟,并与目前应用于电磁三维数值模拟中较好的ICCG迭代算法进行比较。结果表明,多重网格算法收敛速度均与网格节点数或网格尺度无关;而ICCG法的收敛速度则随网格节点数增加迅速减慢。几何多重网格法的计算速度比ICCG法要快数倍,证明了多重网格法三维数值模拟的高效性。 然而以上几何多重网格(GMG)法应用于直流电阻率三维数值模拟遇到很大困难。点源直流电阻率法的数值模拟是一个开域的问题,三维正演中的网格高度非均匀,同时大的电性差异(不连续电性界面)、点电源的奇异性等问题都使得几何多重网格法无法获得其固有的高效性。Zhebel(2006)在其数学系博士论文中改进了传统的几何多重网格法,应用基于矩阵的多重网格法求解带电性差异的直流电阻率三维正演问题,但收效甚微。 代数多重网格(AMG)法近10年发展非常快,在计算流体力学领域的应用表明该方法在非均匀网格和物性系数不连续的情况下依旧能保持其固有的收敛性。论文对直流电阻率三维数值模拟的有限差分格式进行求解,针对大规模非均匀网格、大的电性差异等主要问题进行了研究,实现了电阻率三维数值模拟的代数多重网格算法,并与ICCG法进行系统的比较。大量的模型计算结果表明,AMG算法收敛速度与网格节点数、计算区域、网格精细度、电性差异及模型尺寸等因素基本无关,迭代收敛次数均在7至8次左右,其计算时间随网格节点数增加基本呈较低的线性增长;ICCG法收敛曲线则呈振荡型,受上述各种因素的影响,其迭代收敛次数从数百次变化至数千次,计算时间随网格节点数增加呈非线性快速增长。一般大规模网格情况下,单个点源的三维数值模拟,AMG法的计算速度约为ICCG法的3至4倍;对多个点源进行求解时,如21个点电源,AMG方法的计算速度约为ICCG法的9倍,优势将更加明显。 在电阻率三维正演的代数多重网格快速算法基础上,将其应用于巷道超前探测的三维数值模拟研究中。近年来,国内煤矿屡次发生透水矿难,造成巨大的人员伤亡和财产损失,直流电阻率巷道超前探测面临着迫切的生产需求,然而由于可利用的观测空间有限、巷道掘进面前方异常信号弱以及巷道空腔和旁侧异常体的影响等诸多问题,目前仍是国际上的研究难点。本文在大量地下巷道模型的电阻率三维数值模拟基础上,总结分析了巷道掘进面前方异常体真实位置与观测得到的视电阻率曲线极值点位置之间的关系,首次给出了定量估计巷道掘进面前方异常体真实位置的预测模型。文中还分析了异常体尺度对视电阻率异常曲线的影响,据此,可以参照理论预测模型,定性推断前方异常体的大小。同时,对巷道旁侧异常体的影响及定位也进行了细致的研究。研究结果对巷道超前探测的生产实践有重要的理论指导意义。 多相介质的有效电导率研究在测井资料定量解释中非常重要,也是建立壳幔深部电导率模型的基础。本文应用电阻率三维数值模拟代数多重网格快速算法,发展了一种新的多相介质有效电导率反演计算方法,在两相介质电导率差异变化很宽(达108)的情况下,大量模型的有效电导率计算结果与有限元数值模拟ElecFEM3D的计算结果基本一致,表明本文的计算结果是可靠的。 上述巷道超前探测和多相介质有效电导率计算的应用研究中,均涉及大规模三维网格、大电性差异的大量模型的电阻率三维数值模拟,工作量巨大。因此,电阻率三维数值模拟代数多重网格快速算法发挥了重要作用。 论文最后研究了大型复系数矩阵的代数多重网格算法并探索其在大地电磁三维数值模拟中的应用。结果表明,代数多重网格(AMG)算法在解大型复数方程组时可以获得了多重网格算法的最优效率,但应用于大地电磁三维数值模拟还需要对AMG更深层次的数学问题进行研究才有可能获得更好的效果。
【学位授予单位】:

知网文化
【相似文献】
中国期刊全文数据库 前20条
1 谢伟,陈材侃,程尔升;用于曲线坐标生成的多重网格与有限分析法[J];华中理工大学学报;1995年05期
2 刘之行,封卫兵;三维区域上的多重网格算法[J];西安交通大学学报;1998年12期
3 卫加宁,王伟沧,皮新明,章社生;一类虚拟边界预条件多重网格并行算法[J];武汉理工大学学报(交通科学与工程版);2001年01期
4 葛永斌,吴文权,田振夫;二维波动方程的加权平均隐格式及多重网格算法[J];上海理工大学学报;2002年03期
5 欧阳洁;求解跨音速势流的自适应并行多重网格算法[J];计算机工程与科学;1995年04期
6 葛永斌,吴文权,卢曦;解二维扩散方程的高精度多重网格方法[J];工程热物理学报;2002年S1期
7 刘星,卞恩荣,童谷生;非结构网格上二维欧拉方程的多重网格解[J];华东交通大学学报;1999年03期
8 葛永斌,吴文权,卢曦;求解二维扩散方程的加权平均隐式多重网格方法[J];内蒙古大学学报(自然科学版);2003年05期
9 田小波,吴庆举,曾融生;弹性波场数值模拟的隐式差分多重网格算法[J];地球物理学报;2004年01期
10 冯树常,彭爱和;一个在球面半规则网格上的求解Helmholtz方程的多重网格法[J];气象科技;1994年02期
11 徐正权;多重网格区域分裂分布式计算[J];数值计算与计算机应用;1996年01期
12 纪欢;王宁宁;;二维抛物方程初边值问题的拟多重网格预处理迭代法[J];科技信息;2011年05期
13 李桦,王承尧,王磊;用多重网格TVD方法求解三维高超音速粘性流场[J];国防科技大学学报;1998年02期
14 石钟慈,许学军;Wilson非协调元的V循环多重网格法[J];中国科学A辑;1999年10期
15 石钟慈,许学军;非对称不定问题的非协调多重网格法[J];计算数学;1999年04期
16 刘学强,伍贻兆,夏健;多重网格法在非结构网格中的应用[J];计算物理;2002年04期
17 曾晓清,吴雄华,王云飞;有限体积法与多重网格法用于提高激波捕捉质量及加速收敛[J];应用数学与计算数学学报;1995年01期
18 曾晓清,吴雄华,王云飞;多重网格法与有限体积法在流体力学中的应用[J];同济大学学报(自然科学版);1995年05期
19 黎明,张政;计算流体力学与传热中的多重网格块修正算法[J];北京化工大学学报;1999年03期
20 葛永斌,吴文权,卢曦;基于二维泊松方程六阶紧致格式的多重网格方法[J];上海理工大学学报;2002年04期
中国重要会议论文全文数据库 前10条
1 周磊;谭伟伟;;聚合多重网格的体网格可视化技术[A];探索 创新 交流(第4集)——第四届中国航空学会青年科技论坛文集[C];2010年
2 袁礼;;三维不可压Navier-Stokes方程的多重网格求解[A];第十届全国计算流体力学会议论文集[C];2000年
3 王青平;白武明;王洪亮;;多重网格在二维泊松方程有限元分析中的应用[A];中国科学院地质与地球物理研究所第11届(2011年度)学术年会论文集(上)[C];2012年
4 王斐;赵杰煜;;基于多重网格的水平集图像分割方法[A];第十四届全国图象图形学学术会议论文集[C];2008年
5 刘刚;刘伟;牟斌;肖中云;;涡流发生器绕流模拟的数值方法研究[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集(下)[C];2007年
6 葛永斌;田振夫;吴文权;;三维非定常对流扩散方程的高精度全隐式多重网格方法[A];计算流体力学研究进展——第十二届全国计算流体力学会议论文集[C];2004年
7 朱培烨;;基于非结构网格的欧拉方程聚合多重网格法[A];计算流体力学研究进展——第十二届全国计算流体力学会议论文集[C];2004年
8 李笑郁;毛士艺;;基于相位信息加权的合成多重网格干涉SAR相位展开算法[A];中国体视学学会图像分析专业、中国体视学学会仿真与虚拟现实专业、中国航空学会信号与信息处理专业第一届联合学术会议论文集[C];2000年
9 吴爱弟;牟永光;杨文杰;曹思远;符力耘;;微分方程的小波解法[A];1995年中国地球物理学会第十一届学术年会论文集[C];1995年
10 陈明;柯涛;丁大志;樊振宏;陈如山;;电磁散射问题中的等级基二重网格法[A];2007年全国微波毫米波会议论文集(上册)[C];2007年
中国重要报纸全文数据库 前1条
1 张伟;管理出“妙方” 治理全覆盖[N];泸州日报;2009年
中国博士学位论文全文数据库 前10条
1 刘学强;基于混合网格和多重网格上的N-S方程求解及应用研究[D];南京航空航天大学;2002年
2 阳述林;高维中子输运方程的离散格式与并行算法研究[D];中国工程物理研究院;2003年
3 冯国峰;波动方程反问题的多尺度—信赖域反演方法[D];哈尔滨工业大学;2006年
4 谭勤学;p型多重网格间断有限元及其在燃机冷却系统中的应用[D];清华大学;2013年
5 张挺;基于多点地质统计的多孔介质重构方法及实现[D];中国科学技术大学;2009年
6 李壮;探地雷达反问题的同伦算法研究[D];哈尔滨工业大学;2007年
7 付朝江;集群MPI环境下有限元结构分析并行计算研究[D];上海大学;2006年
8 葛长江;基于交错、同位和非正交网格求解流场的半隐式快速算法研究[D];吉林大学;2010年
9 李文军;块代数多重网格算法的研究及其实现技术在在油藏数值模拟软件中的应用[D];中国科学院软件研究所;2001年
10 牛海峰;流体及非线性最优控制问题的有限元方法:状态受限与超收敛分析[D];山东大学;2011年
中国硕士学位论文全文数据库 前10条
1 李林;一维抛物型方程拟多重网格预处理迭代法及数值试验[D];渤海大学;2013年
2 杨磊;基于PC机群环境MPI的多重网格并行算法研究[D];武汉理工大学;2004年
3 刘鑫;基于PDE的几种图像处理方法及其数值解法[D];大连理工大学;2005年
4 杨艳南;二维抛物型方程初边值问题拟多重网格预处理迭代法[D];渤海大学;2013年
5 赵振宇;求解第一类算子方程的多重网格算法[D];河北工业大学;2004年
6 陈飞;基于多重网格的图像恢复算法[D];福州大学;2006年
7 万建军;椭圆边值问题非协调有限元多重网格法分析[D];郑州大学;2004年
8 陆华才;基于摩擦学理论的缸套激光表面微造型参数优化的研究[D];江苏大学;2005年
9 张杨;图形处理器并行计算应用研究[D];西南交通大学;2006年
10 王军华;多尺度计算方法的初探及应用[D];西北工业大学;2006年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978