电法可视化正反演报告.docx

一维电测深正反演报告院 （系）：地空学院专业年级: 地球物理学2010级 作 者: 061102班 第6组 韩冰凯、黄凯、金焕、凌晓 时间： 2013年4月25日目录第一章前言1第二章一维电测深原理2第三章一维电测深正演6第四章一维电测深反演8第五章总结与思考10参考资料10附录：相关数据处理代码前言（一）本次作业题目为：一维电测深正反演作业要求：学习掌握一维电测深原理，具体模型为：一维水平层状大地电测深模型学习使用Matlab编程，并设计一维电测深正反演程序在设计正反演程序的基础上实现可视化编程软件基本操作和演示（三）小组成员：061102班第6组韩冰凯、黄凯、金焕、凌晓小组分工：美工设计凌晓，小组每人掌握图形界面创建工具GUIDE中各个模块的使用（四）进展与成果：课堂学习一维电测深基本原理，并利用课余时间做拓展学习学习使用matlab编写调试一维水平层状大地电测深程序在3）的基础上实现可视化编程一维电测深原理水平层状大地对称四极电阻率测深曲线多层水平地层上的视电阻率表达式1.水平地层上地面点电流源的电场水平地层上地面点电流源电位的定解问题如图所示，假定地面水平，在地下有n层状地层，各层的电阻率；厚度分别；每层底面到地面的距离为其中，在A点有一电流源供电，电流强度为I用柱坐标系，将原点设在O点，z轴垂直向下。由于问题的解具有轴对称性，与无关，因此电位分布满足以下形式的拉普拉斯方程：及如下边界条件：2.水平地层上地面点电流源电位表达式其中，第一式的解为第一类零阶贝塞尔函数，第二式的解为由此得到（n-1）元方程组（二）水平地层上视电阻率表达式及滤波算法水平地层上视电阻率表达式为：水平地层视电阻率滤波算法，使用一套滤波系数CK，如下表：第三章 一维电测深正演打开水平地电断面反演界面：该界面主要功能有：输入模型的地电参数，依次为地层层数n,层电阻率res和层厚度thick；绘制电测深曲线：根据设置的参数绘制电测深曲线；二、流程图如下：一维电测深正演程序的设计编程采用自顶而底的思想，首先编写伪代码，设计程序构架，然后逐步实现程序函数的功能。流程图中重点展示的是程序如何实现读入数据、保存数据的过程。开始选择要输入的数据类型(层数、厚度)进行正演绘图导出数据结束关闭文件判断指针是否指向最后一个字符建立新文件（*.dat） N Y输出所指数据文件指针加1三、具体操作步骤如下：1.登陆界面，设置地电参数：2.绘制电测深曲线：第四章 一维电测深反演打开水平地电断面反演界面：该界面主要功能有：1.导入实测数据：可将.dat文件中的数据（即实测视电阻率或者正演结果数据）导入matlab中以进行运算；2.绘制实测电测深曲线：根据导入的数据绘制曲线；3.进行正演计算：正演的结果会在绘图框中显示出来；4.进行反演计算：通过设置界面中的电性参数，将反演结果以电测深曲线的形 式表现出来；5.计算相对误差：可以计算实测视电阻率和模型反演视电阻率的相对误差。二、具体操作步骤如下：1.导入实测数据：2.绘制实测曲线：3.进行正演计算：4.修改数据，使正演结果接近观测结果：通过修改数据，然后进行正演，得到误差范围内的模型参数第五章 总结与思考通过老师上课的讲解，我们小组成员在课下进行学习，最终完成了本次作业，对于这一过程的收获和反思总结如下：1.基本掌握了一维电测深的正反演程序及实现可视化编程，达到每个人都能够独立操作，这也是我们最初设计的最低要求；2.对于老师给的源程序，我们对核心代码部分做了研究，并且结合一维电测深的原理，核心程序做了创新和改动，增加核函数的表达方式；3.掌握了可视化编程的基本操作，希望成为今后编程的基础。特别是对于图形界面创建工具GUIDE的使用，从启动到完成甄格格界面的设计都有了深刻的认识。其中，还经过图形的界面设计和修改。4.我们在实践过程中认为，最重要的是实现GUIDE和M代码的联系，在这一过程中我们查找了不少资料，对其编程实现的原理有了系统的认识；5.总结我们的缺点和下一步继续编程的目标：可视化界面需要增加更多的功能和文件菜单，而且希望进一步对界面进行简化，调试软件，编写软件说明书。 计划实现界面主要参考资料：[1]MATLAB从入门到精通 周建兴[2]电法勘探教程 程志平[3]其他网络资料