大地电磁测深有限单元法正演与混合遗传算法正则化反演研究

来源 :中南大学 | 被引量 : 0次 | 上传用户:xiestephen
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
大地电磁测深是以天然电磁场为场源研究地球内部电性结构的一种重要的地球物理手段。自上世纪五十年代提出以来,无论是仪器野外采集、资料处理与解释还是在理论方法的基础研究方面都有了长足的进展。然而,随着理论研究的深入和实际应用的推广,出现了许多新的问题需要解决或进一步研究完善,这些问题涉及到正演的数值模拟技术和反演算法研究。本文主要讨论了有限单元法的大地电磁测深正演算法和实数编码混合遗传算法的大地电磁测深正则化反演算法。从麦克斯韦方程组出发,推导出二维大地电磁测深的边值问题、变分问题、有限单元法的插值函数以及单元刚度矩阵。结合散度条件,采用广义变分原理推导出三维大地电磁测深的有限元方程。有限单元法最后归结于解线性方程组Ax=b,这里的A是大型稀疏的带状对称复系数矩阵,并且其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难。不完全LU分解预处理的BICGSTAB算法可用于该线性方程组的求解,并且具有速度快、精度高和稳定性好等优点。为了验证大地电磁测深二维、三维正演算法的正确性,用二维、三维正演模拟结果和一维大地电磁测深正演结果进行对比。另外,对COMMEMI组织提供的模型进行正演计算,将模拟的结果与其对比研究,也进一步对二维、三维的正演算法进行评估。正演模拟结果表明:本文设计的三维大地电磁测深正演算法是非常有效的,该方法具有很多优点:速度快、精度高、解题过程比较规范、能自动满足内部边界条件、适合于各种复杂的物性分布问题。大地电磁测深的反演问题是不适定的,因此反演结果是不稳定的且具有非唯一性,这也就意味着不同的地电模型对观测数据的拟合具有同样的精度。在目标函数中采用正则化方法,可以使得不适定反演问题具有稳定的反演结果。设计了基于正则化反演的实数编码混合遗传算法,避免了非线性问题的线性化困难,并较好的解决了大地电磁测深反演的不适定性和非唯一性。通过反演模拟研究,验证了本文反演算法在大地电磁测深资料反演的有效性和稳定性。另一方面也指出,大地电磁测深反演问题是一个高度非线性的不适定问题,基于正则化的实数编码混合遗传反演方法本身并不能解决反演解的非唯一性,因为这种非唯一性是由于观测数据的缺乏,物理问题本身所固有的。然而该反演方法是非常稳健的,通过反演模拟研究,我们认为这种全新的反演方法给大地电磁测深反演问题的解决带来了新的希望,值得深入研究。
其他文献
对于一类具有等时中心的二次系统x=-y+x^2-y^2,y=x+2 xy,应用Abel积分和Picard-Fuchs方程,得到了该等时系统在n次实多项式扰动下从相应周期环域中分支出极限环个数的上界.
目的:初步分析DNA测序结果.方法:采用Visual Basic 6.0编程语言编写了DNA碱基序列分析软件.结果:此软件具备计数碱基、自动查找开放阅读框架(open reading frame,ORF)、翻译
目的:探讨急性脑血管病患者应激状态下外周血皮质醇、淋巴细胞和T细胞亚群的变化及临床意义.方法:将70例急性脑血管患者分为有应激反应组(A组,34例)和无应激反应组(B组,36例)