由大地电磁数据同时反演电阻率和磁化率参数的方法研究

论文题目: 由大地电磁数据同时反演电阻率和磁化率参数的方法研究

论文类型: 博士论文

论文专业: 地球探测与信息技术

作者: 李学民

导师: 曹俊兴

关键词: 大地电磁法,同时反演,同相分量,异相分量,灵敏度,电阻率,磁化率,非线性共轭梯度法,大杨树盆地

文献来源: 成都理工大学

发表年度: 2005

论文摘要: 大地电磁法(MT)是重要的深部地球物理探测手段,在油气、地热勘探,深部探测等领域有着广泛的应用。传统的MT数据处理通常假定地下介质无磁性,仅对介质的电阻率参数进行反演。本文以国家自然科学基金重大项目为依托,在分析计算岩石磁化率参数对大地电磁数据影响的基础上,提出了由MT数据同时反演介质电阻率和磁化率参数分布的方法,并用所发展的方法技术成功地进行了实测资料的处理。论文的主要内容与成果如下: (1) 推导了电阻率和磁化率同为参数时的二维大地电磁正演计算公式。利用Madden和Mackie等所提出的基于传输网络的有限差分算法,实现了其正演计算,并详细地阐述了电阻率和磁化率对大地电磁响应的影响,证实了磁化率主要是影响磁场信号的同相分量,而对磁场信号的异相分量和电场信号的影响比较弱。当磁化率k≥0.01SI时,在某些低频点处视电阻率值可以增加5Ωm左右,并且随着磁化率的增加,这种增加的幅度显著增大;而阻抗相位的变化则比较小。通过对比研究,认为在大地电磁勘探中,起决定作用的是岩石的电阻率;但是在高磁性物质存在的地区进行大地电磁观测时,还必需重视磁化率的影响。 (2) 提出了由大地电磁实测数据同时反演电阻率和磁化率分布的方法技术。考虑电阻率和磁化率同为参数变量,推导了由MT数据同时反演电阻率和磁化率参数分布的反演计算公式。在同时反演过程中,观测数据和物性参数之间以及电阻率和磁化率物性参数本身之间是一种复杂的非线性关系。本文对电阻率和磁化率之间的非线性关系利用一个线性函数来处理,重点是采用了Rodi和Mackie所发展的非线性共轭梯度算法来处理观测数据和物性参数之间的非线性关系,完成了新反演方法的程序代码编写。 (3) 推导了二维大地电磁法中TM和TE两种极化模式的观测数据关于模型参数的灵敏度表达式。利用Mora提出的计算灵敏度的思想,根据Rodi、Mackie和Madden发展的算法,给出了观测数据关于模型参数灵敏度的计算表达式和求解过程。 (4) 通过数值实验和实际观测资料实验,证明了本文研究发展的算法和计算程序的正确性与价值。使用所研发的由MT数据同时反演电阻率和磁化率参数分布的方法技术对内蒙古大杨树盆地9条MT测线的观测资料进行了重新处理,厘定了电阻率分布图像,并新获得了磁化率的分布图像。反演结果表明,在火山岩分布地区,由MT数据获得的磁化率分布信息更有助于推断火山岩的分布,并能为在火成岩覆盖地区寻找沉积岩的分布提供一种新的依据。利用新发展的方法技术同时还可以解决在某些大地电磁观测数据中出现负的同相分量值问题,使得这些出现在低频段的负的同相分量值不再是一种“污染源”,而是作为一种提供磁化率分布很重要的信息源。

论文目录:

摘要

ABSTRACT

目录

序言

1 电阻率和磁化率同为参数的二维MT正演问题

1.1 二维大地电磁正演问题的计算公式

1.2 边界条件的确定

1.3 数值模拟实验

1.4 小结

2 岩石电阻率和磁化率对大地电磁响应的影响

2.1 岩石的电阻率和磁化率

2.2 大地电磁响应的计算

2.3 数值模拟实验

2.3.1 岩石磁化率对大地电磁响应的影响

2.3.2 岩石电阻率对大地电磁响应的影响

2.4 小结

3 电阻率和磁化率同为参数的二维MT反演方法

3.1 常用非线性反演方法

3.1.1 梯度法

3.1.2 牛顿法

3.1.3 共轭梯度法

3.2 大地电磁数据的二维反演方法

3.3 求目标函数极小值的算法

3.3.1 高斯—牛顿迭代法

3.3.2 Mackie—Madden算法

3.3.3 非线性共轭梯度法

3.4 预处理因子

3.5 模型函数的转换

3.6 小结

4 计算观测数据关于参数的灵敏度

4.1 计算Fréchet导数的基本原理

4.2 计算灵敏度的方法

4.2.1 扰动法

4.2.2 灵敏度方程法

4.2.3 伴随方程法

4.3 计算灵敏度的解析解

4.4 计算灵敏度的数值解

4.5 计算雅克比矩阵的过程

4.6 小结

5 数值实验与实测资料的处理解释

5.1 正则化因子和加权系数的选择

5.2 数值实验

5.3 实测数据处理与解释

5.3.1 工区的概况

5.3.2 观测数据的定性解释

5.3.3 观测数据的博斯蒂克(Bostick)反演

5.3.4 观测数据的二维反演

5.4 小结

结论与建议

致谢

参考文献

发布时间: 2005-09-09

参考文献

• [1].基于自适应矢量有限元法的三维大地电磁正反演研究[D]. 秦策.成都理工大学2018
• [2].大地电磁正则化反演方法研究[D]. 张昊.中国地质大学(北京)2018
• [3].中国东北地区新生代火山深部电性结构及流变分析[D]. 李建平.吉林大学2018
• [4].长周期大地电磁数据处理方法研究[D]. 张刚.成都理工大学2015
• [5].大地电磁信号处理理论及方法研究[D]. 严家斌.中南大学2003
• [6].大地电磁数据非线性反演方法研究[D]. 尹彬.中国地质大学2017
• [7].大地电磁测深静态位移数据的反演方法研究[D]. 周军.成都理工大学2016
• [8].井震约束的大地电磁-重力联合反演研究[D]. 胡祖志.中国地质大学2017
• [9].基于Hilbert-Huang变换的大地电磁信号处理方法与应用研究[D]. 蔡剑华.中南大学2010
• [10].三维大地电磁有限内存拟牛顿反演[D]. 阮帅.成都理工大学2015

相关论文

• [1].大地电磁测深有限单元法正演与混合遗传算法正则化反演研究[D]. 童孝忠.中南大学2008
• [2].大地电磁聚焦反演成像方法研究[D]. 刘小军.同济大学2007
• [3].地球物理非线性联合反演方法研究[D]. 敬荣中.中南大学2002
• [4].三维电阻率/极化率有限元正演模拟与反演成像[D]. 黄俊革.中南大学2003
• [5].复电阻率三维电磁场正反演研究[D]. 张辉.吉林大学2006

标签：大地电磁法论文;  同时反演论文;  同相分量论文;  异相分量论文;  灵敏度论文;  电阻率论文;  磁化率论文;  非线性共轭梯度法论文;  大杨树盆地论文;