学位论文 > 优秀研究生学位论文题录展示
电性参数分块连续的大地电磁二维有限元数值模拟
作 者: 尤淼
导 师: 王绪本
学 校: 成都理工大学
专 业: 地球探测与信息技术
关键词: 大地电磁测深法 电性参数分块连续 二维正演模拟 地形起伏 有限单元法 线性插值 矩形三角形网格剖分 变带宽矩阵存储
分类号: P631.325
类 型: 硕士论文
年 份: 2011年
下 载: 71次
引 用: 0次
阅 读: 论文下载
内容摘要
大地电磁测深法(MT)属于电磁法勘探中的频率域方法,利用天然交变场源对地球深部岩石的电性参数进行研究。目前,MT法的二维正演问题已基本解决,在使用有限单元法,有限差分法和积分方程法等数值模拟方法解决二维MT正演问题方面有广泛的研究结果,但在网格剖分方式,电性参数设定,辅助场定义和起伏地形模拟等方面,仍有改进之处。本文研究电性参数分块连续变化情况下的MT二维正演问题,推导出二维水平和起伏地形条件下大地电磁法有限元数值模拟算法,以提高正演模拟精度,为反演解释工作的进行提供良好的基础。本文以频率域Maxwell方程组为基础,并结合大地电磁测深法的边界条件探讨,推导出二维MT法的变分方程。在有限单元法网格剖分方式上,采用矩形网格内剖分三角形网格的方案,这种剖分方式便于对倾斜异常体和倾斜地形的模拟,以适应各种水平或起伏地形情况。考虑到实际地层中的岩石,矿物体等在水平方向和垂直方向上电性参数是连续变化的,而在一些反演方法中,反演结果的电性参数也是连续变化,故将网格单元内的电性参数设定为线性变化。根据单元节点主场值和线性插值形函数间的关系,计算出单元节点辅助场值。在方程组的求解方面,采用变带宽存储解决含有大量零元素的大型稀疏矩阵的存储和方程组的求解问题,以节约内存使用量和提高计算速度。根据起伏地形情况下实测电磁场分量的特征,定义TE和TM两种模式下的视电阻率和阻抗相位计算公式。编制一套二维大地电磁法的正演程序,实现了使用高程信息自动生成网格和使用不同的模型参数计算指定频率的视电阻率和相位值。通过对多种模型的验证,不同地电断面的水平地形和起伏地形的正演模拟结果与前人模拟结果一致,模型参数基本吻合。结合两种极化模式的视电阻率和相位信息的横向和纵向分辨率特点,对地下异常体的深度定位,规模大小和倾斜方位判定表现出良好的效果。山峰、山谷和斜坡等非水平地形情况下的正演响应结果也与正演模型基本符合,验证了本文方法的正确性和有效性。
|
全文目录
摘要 4-6 Abstract 6-10 第1章 引言 10-13 1.1 研究背景 10 1.2 国内外发展现状和趋势 10-11 1.3 选题依据和研究意义 11-12 1.3.1 选题依据 11 1.3.2 研究意义 11-12 1.4 研究目标和主要研究内容 12-13 1.4.1 研究目标 12 1.4.2 主要研究内容 12-13 第2章 大地电磁二维变分问题推导 13-21 2.1 基本方程 13-15 2.2 边界条件 15-18 2.2.1 外边界条件 15-17 2.2.2 内边界条件 17-18 2.2.3 边值问题方程组 18 2.3 变分问题 18-20 2.4 本章小结 20-21 第3章 有限单元法 21-34 3.1 网格剖分 21-22 3.2 三角形单元线性插值 22-26 3.2.1 面积坐标 23-25 3.2.2 线性插值 25-26 3.3 单元分析 26-29 3.4 总体合成 29 3.5 求变分 29-30 3.6 简化系数矩阵 30-31 3.7 代入边界条件 31-32 3.8 解线性方程组 32 3.9 本章小结 32-34 第4章 辅助场和视电阻率计算 34-37 4.1 计算辅助场 34 4.2 视电阻率和阻抗相位计算 34-36 4.2.1 水平地形 35 4.2.2 起伏地形 35-36 4.3 本章小结 36-37 第5章 正演程序编制 37-44 5.1 Fortran90 语言 37 5.2 程序输入和输出 37-40 5.2.1 程序输入部分 37-39 5.2.2 输出部分 39-40 5.3 程序流程图 40-42 5.3.1 主程序流程图 41-42 5.3.2 有限元子程序流程图 42 5.4 程序源代码 42-43 5.5 本章小结 43-44 第6章 模型验证 44-69 6.1 一维水平层状模型 45-48 6.1.1 模型一 45-47 6.1.2 模型二 47-48 6.2 二维水平地形模拟 48-63 6.2.1 模型三 49-51 6.2.2 模型四 51-54 6.2.3 模型五 54-55 6.2.4 模型六 55-57 6.2.5 模型七 57-60 6.2.6 模型八 60-63 6.3 起伏地形模拟 63-68 6.3.1 模型九 64-65 6.3.2 模型十 65-66 6.3.3 模型十一 66-68 6.4 本章小结 68-69 结论 69-70 致谢 70-72 参考文献 72-74 攻读学位期间取得学术成果 74-75 附录 75-80
|
相似论文
- 便携式风速风向测量仪,P414.7
- 人脸自动美化算法研究,TP391.41
- 攀钢钢渣作为大体积地基回填料的数值模拟研究,TU521.4
- 基于ANSYS二次开发的高填方拱涵的受力性状分析,U449
- 单梁双吊钩门式起重机结构分析与优化设计,TH213.5
- 基于NUSS手术的漏斗胸矫正手术方案研究,R655
- 弹性波有限差分数值模拟及井间地震逆时偏移成像研究,P631.4
- 里伍铜矿采空区处理技术研究,TD862.1
- 寒冷地区开敞式水闸的温度应力分析,TV66
- 满堂红桥加固提载的数值模拟分析,U445.72
- 椒江大桥拓宽改造方案研究,U445.6
- 板式无砟道岔岔区板力学特性分析,U213.6
- 城轨车辆切削式专用吸能装置研究,U270.8
- 基于无线移动技术的交互式多媒体E-leaning系统设计与实现,TP391.6
- 有限单元法在跨断层水准变化机理研究中的应用,P315.2
- 考虑土—结构相互作用的建筑结构水平地震作用研究,TU435
- 全自动捆钢机关键零部件设计研究,TG96
- 一种5KW离网小型风力发电机组的工作可靠性研究,TM315
- 增压锅炉用SiC耐火砖墙的工程热分析,TK425
- 超长高层混凝土结构的等效综合温差效应研究,TU375
- 既有隧道衬砌病害评估与治理决策,U457.2
中图分类: > 天文学、地球科学 > 地质学 > 地质、矿产普查与勘探 > 地球物理勘探 > 电法勘探 > 方法 > 电磁法
© 2012 www.xueweilunwen.com
|