下载文档原格式
基于偏导数和统计学方法的电阻率测深二维反演程勃;底青云【摘要】为了获得详细直观的视电阻率测深数据的反演结果，本文研究了利用视电阻率测深数据的偏导数变化和统计学参数，设计初始模型的反演方法。
该方法先判断地层电性结构类型，再圈定目标体位置。
利用包含地下地层参数的偏导数信息和一维反演结果的统计学特征，设计出有层参数和地层结构的初始模型。
在此基础上使用遗传算法修改模型，配合二维有限元法正演达到拟合误差最小，从而完成电阻率测深数据的二维反演。
此法的优点在于初始模型与真实的地层结构有很大的相似性；不需要解大型线性方程组，降低计算量；反演过程中独立计算总体拟合误差和异常体影响区域拟合误差，指导优先变异方向。
反演后给出直观的解释结果。
对理论模型及实测数据试算证明该方法可以获得较实用的反演效果。
【期刊名称】《石油地球物理勘探》【年(卷),期】2012(047)006【总页数】8页(P1006-1013)【关键词】电阻率测深;二维反演;偏导数;统计学;遗传算法【作者】程勃;底青云【作者单位】中国科学院地质与地球物理研究所中国科学院工程地质力学重点实验室,北京100029 中国科学院研究生院,北京100049;中国科学院地质与地球物理研究所中国科学院工程地质力学重点实验室,北京100029【正文语种】中文【中图分类】P6311 概述电阻率测深方法广泛应用于大地分层及探测地质构造等问题。
电阻率测深数据定量解释是人们非常感兴趣的研究课题。
多年来人们研究了多种正反演方法，也取得了很好的效果。
尤其是正演模拟的有限元计算方法已经获得了成功[1～6]，在普通的个人计算机上就可以完成二维、三维地质模型的计算。
在二维、三维电阻率测深数据反演问题中，初始模型设计往往成为困扰反演算法的问题，设计模型太简单，可能会无法接近拟合误差极小，而设计模型复杂人为干扰增多[7～11]。
如果可以了解地层结构的类型，进而设计出最接近实际地层结构的初始模型，反演计算的速度和效果会有较大改善，并且反演结果准确、可信。
高密度电阻率法二维和三维有限差分正演计算高密度电阻率法（High-Resistivity Method，简称HRM）是一种用于精确描述电流在较大空间内散射和复杂反射情况下的分布和阻抗、电磁场和电动势场强度的方法。
它可以精确地计算出沿一个给定路径上电流和电场强度的分布，从而获得某一个时刻电阻地层的特点。
有别于PT和MT方法，HRM有明显优势，可以处理电磁参数和地质参数的复杂变化，适应复杂的反射现象，以及更好地解析各种地表以下界面及其结构形成的脉动反射特征。
HRM二维和三维正演计算都基于相同的数学模型，有一个主要的区别，即边界条件的处理不同。
二维正演计算是基于一维电磁方程组和单边带边界条件，而三维正演计算基于3维电磁方程组，具有等张场及大气层的沿着深度的变化。
HRM的数值正演模拟具有以下特点：（1）快速收敛：基于电磁方程组，通过精确描述电极间耦合时，就可以获得快速收敛。
（2）反射损失小：以解耦形式处理反射，它们的损失更小，从而可以更好地获得精确的结果。
（3）地质参数可调：由于地表以下空间有不同的地质参数，如电导率、磁导率等，只需简单地修改这些参数，就可以很好地模拟不同的深部地质结构。
HRM正演模拟分为2步完成：第一步建立模型，根据真实的地质情况设置正演路径，并确定正演模型的地表，地表以下的地质学参数及边界条件；第二步进行数值正演模拟，根据建立的模型，使用有限差分法的方法解决电磁方程，获得强度场及电场、磁场，并与水文地质诸要素参数建立联系以便反演出水文地质参数。
HRM正演模拟应用非常广泛，能够获得较精确的电流分布和阻抗、电磁场和电动势数据集，对地质勘查提供了有效的技术手段，例如辅助高精度的地质资料的获取和评价，预测浅层以及深部复杂地质结构非常有用，可帮助开发者识别和定位储层分布，有助于钻井现场安全和高效地进行。
点源二维直流电阻率法正演模拟吴曲波;张志勇;柯丹;范利飞【摘要】从点源二维电场边值问题出发,利用傅里叶变换将三维问题转换为二维问题,而后使用有限元法求解正演变分问题得节点电位.对三种常用解大型有限元线性方程组解法的计算机时进行试算统计,认为共轭梯度法计算速度较快.最后,通过对比均匀半空间模型有限元数值解和解析解的误差,以及对四个模型分别进行正演模拟试算,算例结果表明本文采用算法具备较好的模拟效果.【期刊名称】《工程地球物理学报》【年(卷),期】2014(011)001【总页数】6页(P44-49)【关键词】点源二维直流电阻率法;正演模拟;有限元法【作者】吴曲波;张志勇;柯丹;范利飞【作者单位】核工业北京地质研究院中核集团铀资源勘查与评价技术重点实验室,北京 100029;中南大学地球科学与信息物理学院,湖南长沙 410083;核工业北京地质研究院中核集团铀资源勘查与评价技术重点实验室,北京 100029;劳雷仪器有限公司,北京 100020【正文语种】中文【中图分类】P631.31 引言点源二维电阻率法作为电法勘探技术的分支，广泛应用于水、工、环等地球物理勘探领域。
开展点源二维电阻率法正演研究，能够提高反演计算精度、地质解释准确性，更好地指导实际生产。
电法正演模拟方法分三类：解析法、物理模拟法和数值模拟法［1］。
由于解析法和物理模拟法不适合场源和物性分布复杂的正演情况，数值模拟法已成为电阻率法正演最主要的方法。
电法数值模拟方法有：有限差分法，有限元法［2，3］，积分方程法［4］，边界单元法［5］。
有限差分法不涉及计算耗时的偏导数运算，计算简便速度较快，但当物性参数分布复杂或场域的几何特征不规则时，适应性比较差。
积分方程法和边界单元法由于其降维的特点具有较高的计算效率，但只适用简单模型的正演计算。
有限元法具有较高的求解精度且理论较为成熟，能够适应任意模型形态。
因此，有限元法现已是解电法正演问题的主流方法［6，7］，本文利用有限元法正演模拟点源二维电场，采用三角形剖分和线性插值方式，取得了较好的模拟效果。
2D电阻率正反演成像在水平和垂直模型上的异常响应研究2D电阻率正反演成像在水平和垂直模型上的异常响应研究2D电阻率成像技术近年来被广泛应用于工程、水文、环境和矿产等领域,在实际应用中它具有多种的装置类型,不同的装置类型对特定的地质情况有不同的应用效果.本文通过采用Wenner、Wenner-Schlumberger和dipole-dipole三种排列装置在一个水平和一个垂直模型上正演模拟和块反演,揭示了不同的排列装置在水平结构和垂直结构的异常响应,并对高阻体和低阻体进行了模拟.结果显示当采用Wenner,Wenner-Schlumberger和dipole-dipole数据采集技术时,不同的模型结构异常响应有明显的不同,三种排列类型对地下低阻体的分辨能力均高于高阻体,Wenner装置和dipole-dipole装置分别对水平层状结构和垂直结构有较好的分辨力,wenner-Schlumberger装置对水平层状结构有好的分辨力,对垂直结构有中等的分辨能力.作者：李天成牛滨华孙春岩赵丽 LI Tian-cheng NIU Bin-hua SUN Chun-yan ZHAO Li 作者单位：李天成,LI Tian-cheng(中国地质大学(北京)地球物理与信息技术学院,北京,100083;有色金属矿产地质调查中心,北京,100012)牛滨华,NIU Bin-hua(中国地质大学(北京)地球物理与信息技术学院,北京,100083)孙春岩,SUN Chun-yan(中国地质大学(北京)工程技术学院,北京,100083)赵丽,ZHAO Li(中国煤炭地质总局物探研究院,涿州,720750)刊名：地球物理学进展ISTIC PKU英文刊名：PROGRESS IN GEOPHYSICS 年，卷(期)：2007 22(3) 分类号：P631 关键词：2D电阻率勘探装置排列类型数值模拟平滑约束反演块反演。
电阻率层析成像的二维改进粒子群优化算法反演张倩;王玲;江沸菠【摘要】Particle swarm optimization ( PSO) is a global random search algorithm put forward by simulating the flock foraging in the process of social behavior based on swarm intelligence. Researchers have proved that PSO algorithm is an effective geophysical inversion method, and it does not rely on the initial model. Because the conventional PSO is easy to be stuck in relative extremum, slow conver⁃gence speed in the late and the inversion accuracy is not high, this paper presented an improved fully chaotic oscillations particle swarm optimization algorithm based on same conventional PSO theory. It improved the formula of updating speed, made the particles getting the difference between the current global best position quickly, enhanced the learning ability of particles. The paper did a two⁃dimen⁃sional numerical test on ERT data in matlab2012b programming environment,the results show that this algorithm inversion is not de⁃pendent on the initial model, increases the search space,and have higher inversion in accuracy than the standard PSO, and the image quality is better than that of Levenberg⁃Marquardt method.%粒子群优化算法（ PSO）是通过模拟鸟群觅食过程中的社会行为而提出的一种基于群体智能的全局随机搜索算法，已有研究学者证明PSO算法是一种有效的地球物理反演方法，不依赖初始模型。
1、程序简介RES2DINV是一种能自动确定电子成象测量资料的地下二维电阻率模型的最小二乘法计算机反演计算程序，适用于二维电阻率&激发极化资料快速反演，可用于约25～650个电极采集的大型数据（约100～5000个数据点）资料反演。
本反演计算程序除适用于电极布设于地面的正常勘查外，还可用于水下及跨孔高密度电阻率法勘查。
由于它是基于WINDOWS的程序，能支持任何与WINDOWS兼容的图形卡或打印机。
程序在1600X1200象素、256色的显示器上测试通过。
2、计算机系统需求由于反演计算需要运行二维逼近模拟和最小二乘子程序，本程序设计运行于80386及其以上的IBM PC兼容微机、操作系统为Windows 3.1、Windows 95或Windows NT。
在奔腾机上，进行50个电极的数据资料反演，仅需几分钟。
最低配置：48MB硬盘自由空间、16MB RAM、640×480 SVGA彩色图形系统、Windows3.1 或95操作系统、80386/387 80486DX Pentium, Pentium Pro or Pentium II CPU（或兼容CPU）。
建议使用800×600分辨率（对于14英寸或15英寸监视器）或1024×768分辨率（对于17英寸或21英寸监视器）、256色彩色SVGA图形模式。
如果需要处理多于300根电极、2000个数据点的资料，建议使用32MB以上内存的系统。
3、程序的安装使用在Windows环境下先运行安装盘上的SETUP.EXE安装程序，然后再运行硬盘res2dinv子目录下的JACOBWIN.EXE程序，完成系统安装。
此时，主程序RES2DINV.EXE以及支持文件（GRADWEN, GRADTWO and GRADDIP）、示例文件将被安装在硬盘的res2dinv子目录下。
双击RES2DINV图标，便可运行电阻率反演程序。
《电法资料处理与解释》实验七二维MT正演实验专业名称：地球物理学学生姓名：学生学号：指导老师：冯兵、周建美完成时间：2017.01.05一、实验目的：1、掌握大地电磁法的二维有限元正演理论。
2、明确TE极化模式和TM极化模式特征。
3、编制二维正演程序（参考徐世浙有限元与边界元）。
4、掌握二维反演软件的应用，理解网格剖分的规则，合理的进网格剖分和模型设计。
5、设计不同的二维地电模型，包括一维层状模型进行程序精度的检验，分别就TE和TM模式进行正演，并根据结果分析两种模式的异同。
二、实验内容：1、均匀半空间及层状介质模型模拟，并和一维正演程序进行比较，分析误差，以检验二维正演软件的正确性。
2、设计一个低阻体和高阻体分别就TE和TM模式进行正演。
3、设计一个倾斜断层模型。
4、设计一个地堑地电模型。
5、其他模型。
三、实验要求：1、能够正确应用软件，掌握各参数的设置。
2、了解软件设计的基本思路，能够自己调整修改相关输出输入参数。
3、对设计模型进行计算，并绘制相关误差分析剖面曲线，进行误差统计，并对模拟结果做以分析。
4、完成实验报告，word排版，图标清晰，分析合理。
四、结果及分析（一）、均匀半空间1、均匀半空间TM二维正演模型其中电阻率为50，如图1：图表1均匀半空间正演模型2、均匀半空间TE模式视电阻率等值线图（图2）：（以下等值线图的纵坐标为f 的对数，横坐标为模型测量点位号，值为视电阻率）图表1均匀半空间TE模式视电阻率等值线图均匀半空间TM模式视电阻率等值线图，如图3：图表3均匀半空间TM模式视电阻率等值线图（二）、断层模型1、断层模型TM二维正演模型，左上侧电阻率200，右下侧电阻率50.如图4：图表2断层TM二维正演模型2、断层模型TE模式视电阻率等值线图（图5）：（以下等值线图的纵坐标为f的对数，横坐标为模型测量点位号，值为视电阻率）图表3断层TE模式视电阻率等值线图3、断层模型TM模式视电阻率等值线图（图6）：图表4断层TM模式视电阻率等值线图（三）低阻模型1、低阻模型TM二维正演模型其中，低阻体电阻10，背景电阻率50，如图7：图表5低阻模型TM二维正演模型2、低阻模型TE模式视电阻率等值线图（图8）：（以下等值线图的纵坐标为f的对数，横坐标为模型测量点位号，值为视电阻率）图表6低阻模型TE模式视电阻率等值线图3、低阻模型TM模式视电阻率等值线图（图9）：图表7低阻模型TM模式视电阻率等值线图（四）、高阻模型1、高阻模型TM二维正演模型其中，高阻体电阻率100 ，背景电阻率50，如图10：图表8高阻模型TM二维正演模型2、高阻模型TE模式视电阻率等值线图（图11）：（以下等值线图的纵坐标为f 的对数，横坐标为模型测量点位号，值为视电阻率）图表9高阻模型TE模式视电阻率等值线图3、高阻模型TM模式视电阻率等值线图（图12）：图表10高阻模型TM模式视电阻率等值线图（五）、地堑模型1、地堑模型TM二维正演模型其中，上部电阻率，下部电阻率，如图16：图表13地堑模型TM二维正演模型2、地堑模型TE模式视电阻率等值线图（图17）：（以下等值线图的纵坐标为f 的对数，横坐标为模型测量点位号，值为视电阻率）图表14地堑模型TE模式视电阻率等值线图3、地堑模型TM模式视电阻率等值线图（图18）：图表15 地堑模型TM模式视电阻率等值线图五、结论及分析（一）均匀半空间均匀半空间模型的二维MT正演试验的TE和TM模式的视电阻率计算结果在50欧米左右，符合模型结果。
基于全区视电阻率的线源FCSEM二维正演模拟
李杰;柳建新;童孝忠;刘颖
【期刊名称】《地球物理学进展》
【年(卷),期】2012(27)1
【摘要】频率域可控源电磁是在大地电磁测深的基础上发展起来的一种人工源电磁测深法,其二维电磁响应的计算须采用数值模拟方法.本文以Matlab为程序编译工具,采用双二次插值的有限单元法,推导出相应的计算公式.为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分.在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间.针对频率域可控源电磁法中卡尼亚电阻率在过渡区和近区畸变的问题,给出了全区视电阻率的迭代公式,并对典型的一维层状模型以及简单二维模型进行了计算.过渡区和近区数据经过校正后,可以正确反映出模型的地电特征,证明了线源下近区勘探的可能性.【总页数】6页(P343-348)
【关键词】线源;频率域可控源电磁法;有限单元法;正演;全区视电阻率
【作者】李杰;柳建新;童孝忠;刘颖
【作者单位】中南大学信息物理工程学院
【正文语种】中文
【中图分类】P631
【相关文献】
1.基于双二次插值的
2.5维FCSEM有限元正演模拟 [J], 柳建新;汤文武;童孝忠
2.CSAMT全区视电阻率电场正演迭代拟合近场校正方法 [J], 詹少金;钱美平;冯戋戋
3.随钻电磁波电阻率测井视电阻率提取及正演模拟分析 [J], 王滨涛;吴锡令;苑娜;丁柱;王晓星
4.基于全区视电阻率瞬变电磁一维正演中滤波系数的选取 [J], 郭嵩巍;郑凯;刘小畔;王成;王淼
因版权原因，仅展示原文概要，查看原文内容请购买。
