地球电磁场的连续与间断有限元三维数值模拟

发布时间：2018-11-21 10:01

【摘要】：电磁法是地球物理勘探中的重要分支。该方法以地球导电媒质的物性差异为基础,通过观测和分析在其中传播的电磁场,从而预测测区的地下电性结构。经过一百多年的发展,电磁法勘探技术己广泛应用于各个勘探领域。随着油气和固体矿产资源的日益减少,勘探深度越来越大,地形地质条件越来越复杂,勘探难度也越来越大。在工程应用中,施工工艺和质量要求的提高也对电磁法勘察提出了更高的要求。除了需要有更高精度的电磁法勘探仪器外,也要求正演算法能应对更复杂的地质模型,能提供更高的计算精度和效率。本文研究的正是地球电磁场的高效、高精度算法,具体研究内容包括直流电场、频率域可控源电磁场和时间域电磁场的三维正演。采用的数值计算方法是基于非结构四面体网格的有限单元法,包括连续和间断有限元方法,二者的本质区别在于前者采用连续的有限维函数空间逼近真实解,而后者的解空间不连续。点电源直流电场的正反演技术已相当成熟,而长电极源电场的正演计算却依然存在许多难点。其中最大的难点是长电极源穿过非均匀媒质时其表面的电流密度难以准确计算。此外,在大测区中构建穿过非均匀媒质的细长电极模型是比较困难的,而任意形状和倾角的长电极则更难构造。同时,长电极的电阻率远远低于围岩电阻率,若考虑长电极的电阻率进行正演计算,将会带来一个非常病态的系统,影响求解的收敛速度和数值解的精度。我们在长电极源直流电场正演中将长电极源近似为电阻率可忽略不计的有限长线电流源,既避免了悬殊的电阻率差异带来的困难,也免去了构建长电极网格模型的麻烦。利用电磁场基本原理推导的长电极源表面电流密度近似公式允许正演模型中长电极穿过非均匀媒质。此部分研究采用基于节点基函数的连续有限元方法实现了任意形状多段长电极源的直流电场,基于总电位和二次电位的正演算法均由解析解验证了数值精度。我们首次提出了任意长电极源观测装置的视电阻率公式,分析了长电极源和点电源电场的异同。通过地面和井中的数值模拟算例发现长电极源装置对地下异常体的响应比点电源装置更敏感,具有更强的勘探能力。文中还提出了混合双极供电源装置(即长电极和点电极配合供电)以弥补测区中长电极不足的情况,数值模拟结果显示混合供电极同样能获得比点电极装置更强的异常响应。最后还进行了油田水力压裂模型的长电极直流电法正演,预示了油田套管可用作长电极源进行水力压裂过程的监测。基于有限元方法的频率域电磁场正演难点在于电场法向分量不连续且散度为零,尽管节点有限元方法的应用相当成熟,却不能直接用于求解非均匀媒质中的电磁场问题,而只能先求解连续的矢量势和标量势,再转换成电磁场,而且在求解矢量势和标量势时依然要做散度校正才能得到可靠的数值解。考虑到以上问题,我们在频率域可控源电磁场的三维正演中采用的是基于非结构网格的矢量有限元方法。矢量有限元方法将求解的自由度设置在剖分单元的棱边上,只要求待求物理量切向连续,同时矢量基函数自动满足散度为零的条件,因此是直接求解非均匀媒质中电磁场的极佳选择。此部分正演中以电场双旋度方程作为控制方程,由广义变分原理得到对应的泛函形式再利用矢量有限元离散和求解。截断边界采用了索末菲辐射边界条件,可以得到比齐次边界条件更好的截断效果。通过与一维层状模型的解析解和二维模型的有限元数值解的对比,验证了矢量有限元算法的稳定性和计算精度。我们还研究了海底地形对海洋可控源电磁勘探的影响,数值模拟结果显示,海底的起伏地形会造成较大的异常,可能掩盖由海底油气藏产生的异常,在数据处理时,应考虑海底地形的影响。由于采用的是非结构的四面体网格,我们的频率域可控源电磁正演算法可进行复杂海底模型中低频电磁场的正演计算。间断有限元方法是最近二十年来非常热门的数值计算方法,该方法易于实现高阶插值,具有天然的并行性,具有处理间断问题的能力,在计算电磁学中广受欢迎,但该方法目前还没有引入到低频地球电磁场的计算中。我们完成的第三部分研究正是利用间断有限元方法实现三维时间域高频和低频电磁场的正演,因此既可用于探地雷达正演,也能用于瞬变电磁场的计算,可实现任意电流输入波形的电偶极子、磁偶极子和任意形状大回线源的瞬变电磁正演模拟。在边界的处理上,采用的是各向异性的完全匹配层。引入完全匹配层后出现了频率依赖关系,利用12个辅助变量可以将依赖关系解耦,从而推导了完全匹配层中的时域麦克斯韦方程组。完全匹配层的存在可以高效的吸收外传电磁波而不产生明显反射,因此不需要像常规的齐次边界那样需要将计算域设得非常大。空间离散采用的是层级型的高阶正交基函数,而时间离散则采用高阶龙格-库塔格式。该基函数的正交性使得单元矩阵变成单元对角阵,而高阶的刚度矩阵依然是良态矩阵,计算复杂度大大降低。数值模拟结果表明数值解的误差随着所采用的基函数的阶数增加而指数降低。当采用k阶基函数时,数值解具有接近k+1阶的数值精度。为了将间断有限元的优势最大化,还利用OpenMP库实现了并行计算,测试结果显示其并行效率非常高。
[Abstract]:......
【学位授予单位】：中国科学技术大学
【学位级别】：博士
【学位授予年份】：2016
【分类号】：P631.325

【相似文献】

相关期刊论文 前10条

1 陈二云;马大为;赵改平;周海军;;燃气自由射流的高精度间断有限元数值模拟[J];弹道学报;2009年01期

2 杜明笙,冯铁凯,傅连祥,曹常恕,刘玉兰;解二维输运问题的三角网间断有限元法[J];计算物理;1984年01期

3 李宏;高分辨率间断有限元方法[J];计算物理;2004年04期

4 蔚喜军,周铁;流体力学方程的间断有限元方法[J];计算物理;2005年02期

5 黄灿;陈传淼;;一阶线性双曲组的时空全间断有限元的收敛性[J];湖南师范大学自然科学学报;2006年01期

6 张冬云;杨永;郭永恒;李喜乐;;二维间断有限元混合网格混合算法研究[J];航空计算技术;2009年01期

7 吴迪;蔚喜军;;自适应间断有限元方法求解三维欧拉方程[J];计算物理;2010年04期

8 夏轶栋;伍贻兆;吕宏强;宋江勇;;高阶间断有限元法的并行计算研究[J];空气动力学学报;2011年05期

9 张铁;李铮;;一阶双曲问题间断有限元的后验误差分析[J];计算数学;2012年02期

10 刘兴强;吕宏强;黄文超;;基于高阶间断有限元法的气动噪声数值模拟[J];航空计算技术;2012年06期

相关会议论文 前8条

1 王博;李文成;邓子辰;;基于间断有限元方法的高精度数值计算[A];中国计算力学大会'2010（CCCM2010）暨第八届南方计算力学学术会议（SCCM8）论文集[C];2010年

2 王言金;冯其京;张树道;郝鹏程;;基于动网格的欧拉守恒方程的间断有限元方法研究:1D[A];第十四届全国激波与激波管学术会议论文集（下册）[C];2010年

3 蔚喜军;;间断有限元方法求解流体力学方程组[A];中国工程物理研究院科技年报（2001）[C];2001年

4 程剑;卢耀文;刘铁钢;;基于间断有限元方法与加权本质非振荡格式的多区域耦合算法[A];第十六届全国流体力学数值方法研讨会2013论文集[C];2013年

5 洪振英;袁光伟;;粒子输运方程中的线性间断有限元方法[A];全国计算物理学会第六届年会和学术交流会论文摘要集[C];2007年

6 张磊;袁礼;;求解含源项守恒律组的间断有限元方法及在反应流计算中的应用[A];中国力学学会学术大会'2009论文摘要集[C];2009年

7 赵张益;张庆河;;三维悬沙运动的间断有限元模型[A];中国海洋湖沼学会第十次全国会员代表大会暨学术研讨会论文集[C];2012年

8 王言金;;基于动网格多介质流场的欧拉数值方法研究[A];第七届全国流体力学学术会议论文摘要集[C];2012年

相关博士学位论文 前10条

1 王疆兴;电磁场方程的间断有限元方法[D];湖南师范大学;2015年

2 贺茜君;求解波动方程的间断有限元方法及其波场模拟[D];清华大学;2015年

3 杨军;地球电磁场的连续与间断有限元三维数值模拟[D];中国科学技术大学;2016年

4 陈韵骋;线性弹性力学问题的间断有限元方法研究[D];上海交通大学;2011年

5 季良月;非线性方程间断有限元方法误差分析以及后处理技术[D];中国科学技术大学;2012年

6 冯涛;自适应间断有限元方法求解三维可压缩欧拉方程[D];中国科学技术大学;2013年

7 夏银华;非线性问题的间断有限元方法及其时间离散[D];中国科学技术大学;2008年

8 张福涛;基于间断有限元的等几何分析及应用[D];中国科学技术大学;2015年

9 杨继明;渗流驱动问题间断有限元高效数值方法研究[D];湘潭大学;2007年

10 吴迪;自适应间断有限元方法求解三维欧拉方程的研究[D];中国工程物理研究院;2010年

相关硕士学位论文 前10条

1 郭江;求解矢量辐射传输方程的间断有限元法[D];哈尔滨工业大学;2015年

2 李佳;热传导方程的新型间断有限元方法[D];哈尔滨工业大学;2015年

3 罗娟;可对称化双曲守恒律方程组的龙格库塔间断有限元方法光滑解的先验误差估计[D];南京大学;2013年

4 张峰;二维非协调网格自适应局部间断有限元[D];南京大学;2014年

5 赵京昌;自适应间断有限元方法在中子输运问题中的应用[D];华北电力大学(北京);2016年

6 李灿华;常微分方程的一类强间断有限元法[D];湖南师范大学;2006年

7 郑玉军;非线性微分方程的间断有限元法[D];湖南师范大学;2007年

8 温建;一类稳定混合型间断有限元方法[D];郑州大学;2006年

9 刘彦忠;流动与传热问题的间断有限元法数值模拟[D];上海交通大学;2012年

10 黄灿;一阶双曲组的时空全间断有限元的收敛性[D];湖南师范大学;2006年

，

本文编号：2346637