大地电磁二维CG-OCCAM反演研究

发布时间：2020-05-28 15:43

【摘要】：计算效率和对计算机的要求(比如内存)是反演方法需要考虑的关键问题。本文在融合CG(共轭梯度)和OCCAM二种反演方法优点的基础上,提出了一种新的CG-OCCAM反演方法。CG-OCCAM反演方法以OCCAM反演流程为基础框架,采用共轭梯度反演的“拟正演”思想解决雅可比矩阵的计算问题,不需要计算与存储雅可比矩阵的每个元素,直接得到反演模型的更新量,大大的减少了内存的使用,提高了计算效率。同时,本文对OCCAM最优正则化权重因子的搜索流程进行了改进,得到了一套新的反演流程。二维正演采用有限单元法,二维反演采用上述CG-OCCAM方法,本文实现了大地电磁二维CG-OCCAM反演算法,并使用Fortran语言编程开发了大地电磁二维CG-OCCAM反演程序代码。本文设计了低阻异常体与高阻异常体二种理论模型,理论模型合成数据的试算结果验证了所实现的大地电磁二维CG-OCCAM反演的正确性和可靠性。
【图文】：

× h = j + d 两个方程： e = b × h = j + d 式j = δe，b = ，d = ，并取时间因子为e iωt，则电、 × E = iωμH × H = δ iω H 是介质的磁导率，，δ 是电导率， 是介电常数， 为角频率究的是地下二维电性结构，取其走向为图 2-1 中所示的 x ，y 轴水平并与走向垂直。本文所有的公式推导、模型计

使用的边界条件如表 2-1 所示：方法元方法进行大地电磁正演数值模拟计算。格剖分。TE 模式 T空气层，顶边界，u=1 地表， u n= 0,u = Exp, Hxp u n= u n+ ku = 0 u n+ 表 2-1 两种极化模式下正演使用边界条件列表
【学位授予单位】：中国地质大学(北京)
【学位级别】：硕士
【学位授予年份】：2018
【分类号】：P631.325

【参考文献】

相关期刊论文 前10条

1 林昌洪;谭捍东;舒晴;佟拓;谭嘉言;;可控源音频大地电磁三维共轭梯度反演研究[J];地球物理学报;2012年11期

2 林昌洪;谭捍东;佟拓;;倾子资料三维共轭梯度反演研究[J];地球物理学报;2011年04期

3 张罗磊;于鹏;王家林;吴健生;;光滑模型与尖锐边界结合的MT二维反演方法[J];地球物理学报;2009年06期

4 陈小斌,赵国泽,汤吉,詹艳,王继军;大地电磁自适应正则化反演算法[J];地球物理学报;2005年04期

5 陈小斌,赵国泽;基本结构有限元算法及大地电磁测深一维连续介质正演[J];地球物理学报;2004年03期

6 谭捍东,魏文博,邓明,金胜;大地电磁法张量阻抗通用计算公式[J];石油地球物理勘探;2004年01期

7 谭捍东,余钦范,John Booker,魏文博;大地电磁法三维快速松弛反演[J];地球物理学报;2003年06期

8 谭捍东,余钦范,John Booker,魏文博;大地电磁法三维交错采样有限差分数值模拟[J];地球物理学报;2003年05期

9 刘光鼎;回顾与展望——21世纪的固体地球物理[J];地球物理学进展;2002年02期

10 魏文博;我国大地电磁测深新进展及瞻望[J];地球物理学进展;2002年02期

相关博士学位论文 前1条

1 林昌洪;大地电磁张量阻抗三维共轭梯度反演研究[D];中国地质大学（北京）;2009年

相关硕士学位论文 前2条

1 李晋津;可控源音频大地电磁法二维Occam反演研究[D];中国地质大学（北京）;2013年

2 曾国;大地电磁二维有限元正演数值模拟[D];中南大学;2008年

本文编号：2685465