高密度电法软件用于电测深剖面反演

各位代表，大家好！,高密度电法软件用于电测深剖面反演高建东中国冶金地质总局山东正元地质勘查院2010.3,主要内容,一、普通电测深剖面的高密度软件反演二、高密度电法正反演软件GDF介绍三、反演实例,（一）开展电测深剖面反演的主要优点（二）普通电测深与高密度电法的主要异同点（三）高密度电法剖面自动拟合反演定量解释与普通电测深点曲线量板法定量解释的主要异同点（四）野外工作中需注意的问题（五）反演数据格式（六）编排反演数据的几个注意点（七）反演参数设置的注意点,一、普通电测深剖面的高密度软件反演,1、可将人为规定的测深点位视电阻率转换为地电断面模型的“真”电阻率和“真”极化率；把纵坐标从AB/n转换为深度，将成果图件从拟断面图转变为断面图；2、可消除二维地形影响，减少地形假异常误导；3、可有效使用单边多道三极、双边多道三极或其它效率更高的电极排列，提高野外测深工作效率。,（一）开展电测深剖面反演的主要优点,单边多道三极测深电极排列示意图,三极四极混合测深电极排列示意图,B,（二）普通电测深与高密度电法的主要异同点,1、电极布设方式不同2、数据点密度不同3、数据采集方式不同4、自动化程度不同5、基本原理相同6、电极排列装置相同,1、反演的维空间不同2、对测深曲线左支的要求不同3、对测深曲线数据点的密度要求不同4、所依据的基本理论相同,（三）高密度电法剖面自动拟合反演定量解释与普通电测深点曲线量板法定量解释的主要异同点：,1、不要求测深曲线左支达到1渐近线，最小AB极距可适当放大；2、为了提高深部反演效果，大极距的测深数据不宜太稀疏；3、与单点测深曲线反演不同，参加剖面反演的测深点距和极距设计应综合考虑，尽量取整。一般情况下，测深点距、AB距和MN距最好都设计为10m的整数倍；,（四）野外工作中需注意的问题,4、为了便于野外施工，设计的测深极距宜使所有AB电极均处在剖面的测深点位上；5、MN点位必须与单位电极距相匹配，否则读数据时将出现错误提示；6、为了便于反演程序图像显示，整条剖面各测深点的电极排列必须一致。7、宜使用测绘仪器测放剖面及延长线上的测深点位，使用测绳或皮尺测放MN点位，,例：对称四极测深装置，测深点距50m，最小供电极距AB取100m，最大供电极距AB取4000m，设计a=10m，电测深供电极距与测量极距可选择如下：等比排列，n=AM/MN=2：非等比排列,n=AM/MN常数：,电测深极距参考排列,RES2DINV和GDF软件的输入输出数据均采用ASCII码文本格式，GDF软件可以直接读入重庆奔腾公司和中国地质大学（武汉）的高密度电阻率仪器数据。数据文件可用任何文本编辑程序（建议使用UltraEdit）进行编辑。数据项之间可用逗号、空格、或换行/回车符分隔。反演数据文本格式文件后缀为“.DAT”，程序要求DAT数据文件必须严格按规定的顺序排列，不能随意改变，否则数据读入过程将会出错。,（五）反演数据格式(以RES2DINV和GDF软件为例),XX矿区XX号激电测深剖面20727511充电率ms0.20.02-200401.5319.207166718.9-2001201.51246.18288237.52-2002001.51588.53807741.5218008801.5455.840186442.1180010001.5437.190594160.16128-2200146.5485121080160.1278541160162.7756571720159.029591800147.8030643800147.80306410000,测线名称。限制255个字符以内,单位电极距,排列类型标识（温纳=1，二极=2，偶极-偶极=3，温纳=4，温纳=5，三极(单极-偶极)=6，对称四极=7，赤道偶极-偶极=8）,数据点的总点数,数据点的记录位置类型:0=记录点位于排列的第一个电极处;1=记录点位于排列的中点,激发极化资料标志:（1=有电阻率和极化率；0仅有电阻率,后续的3行不必输入）,IP数据类型,可以使用：充电率、相位角、金属因子、百分频率效应,IP单位，可以使用：ms、mrad、ms/ohm.m、0.01FE,IP积分延时，IP充电率积分时长,点位，MN电极距，(使用偶极-偶极、单极-偶极和温纳-对称四极排列时，增加电极间隔系数n)，实测视电阻率值，实测充电率值。,地形数据标志（0=无地形数据,下面的地形点数据不必输入；1=地形点和物探点都有水平和垂直坐标；2=地形点有水平和垂直坐标，但是物探点在测线方向的距离是沿地表的斜距）,剖面地形点数,地形点的x,z,剖面第一个电极位置的地形点序号,最后以4个0作为结束标志,1、数据点点数要与实际的数据量匹配一致；2、舍弃误差大、不准确的数据时，应该将该数据点的点位、极距、场值同时删除；3、排列(剖面)的首电极是指位于剖面最小点位的电极；4、数据点的第一个电极是指采集该数据所使用的电极中相对位于最小点位的电极；5、数据中的重复观测数据要事先整理，每个观测极距只能保留1组数据。,（六）编排反演数据的几个注意点,点距或极距不匹配的典型错误提示,错误原因：野外施工时没有从整条剖面反演的角度安排极距，致使该数据点电极排列中的某个电极不在以“a=10”为单位的位置上。改正措施：修改a值或舍弃该数据点或在数据中微调电极位置。预防措施：野外施工时从整条剖面反演的角度安排极距，务必使所有的电极均位于a的整数倍位置。,1、首先使用删除坏点功能检查数据质量，删除突变点；2、一般应限制模型电阻率的上下限（默认不限制）；3、一般应设置降低边缘单元影响选项（默认不降低）；4、一般应选择模型电阻率圆滑选项（默认圆滑）,（七）反演参数设置的注意点,二、高密度电法正反演软件GDF介绍,（一）概述（二）GDF软件特点（三）GDF正反演模型（四）模型计算结果,（一）概述,GDF高密度电法及电测深剖面正反演软件是一种全中文的电法模型正演和异常自动反演计算程序。程序采用了有限元法、有限差分法和准牛顿最优化非线性圆滑约束最小二乘算法，在大数据量下的计算速度较常规最小二乘法快10倍以上且占用内存较少。程序采用VC+编制，运行环境为Windows系列操作系统，它不但可以完成各种排列的高密度电法剖面电阻率和极化率模型正演及异常反演，还可以完成普通电测深剖面的自动拟合反演。一次可以处理2000个电极、8000个数据点的高密度电法剖面，允许电极排列的最大电极距可以达到最小电极距的47倍。,（二）GDF软件特点,1)可以直接读取多种高密度电法仪器采集的数据、完成各种电极排列在起伏地形条件下二维模型视电阻率和激发极化数据的反演计算；2)可以选择使用有限差分法或有限元法，以满足对计算时间和计算精度的不同需求；3)具有多种地形模拟网格可供选择：可以选用均匀变形有限元网格、阻尼变形有限元网络或变形有限元变换网格；4)具有多种反演方法可供选择：可以单独选用快速的准牛顿最优化非线性最小二乘新算法或常规高斯牛顿法，也可以混合使用高斯牛顿法和准牛顿法；,5)可使用具有相同尺度的等截面模型单元或随深度增大而逐渐扩大截面积的变截面模型单元，程序允许设置的模型层数可多达27层。6)输入输出数据文件均采用txt文本字符格式，便于使用各种仪器输出的数据；输出的成果图像还可以选择普通图像文件格式和SURFER等值线文件格式。7)具有较强的坏点删除，测线数据拼接、拆分，测线反向、数据编辑等数据预处理功能；8)可控制调节的参数多，便于提高反演计算效果。可调的参数主要有：阻尼系数、阻尼深度系数、阻尼系数优化方式、模型单元电阻率范围、垂向平滑滤波(fz)与水平平滑滤波(fx)的阻尼系数比、模型单元尺度、反演收敛极限、截止均方误差、迭代次数等。,（二）GDF软件特点,（三）GDF正反演模型,递增厚度变截面模型,等截面模型,（三）GDF正反演模型,（四）模型计算结果,模型参数：低阻体从第19测点开始(x=95m)，延续到第24测点(x=120m)。高阻体从第30测点开始(x=150m)，延续到第35测点(x=175m)。,模型正反演结果简表（PIV1.8G/256MB笔记本计算机）,1A双向单极梯度装置正演成果,拟合误差0.103，耗时12.4分,2A正向单极梯度装置(三极装置AMN)正演成果,拟合误差0.146，耗时6.7分,3A反向单极梯度装置(三极装置MNA)正演成果,3B反向单极梯度装置(三极装置MNA)反演成果,拟合误差0.152，耗时6.7分,4A对称四极装置(AMNB)正演成果,4B对称四极装置(AMNB)反演成果,拟合误差0.197，耗时5.0分,5A偶极偶极装置(ABMN)正演成果,5B偶极偶极装置(ABMN)反演成果,拟合误差0.155，耗时9.7分,6A电位装置二极AM正演成果,拟合误差0.184，耗时8.7分,7A温纳装置(AMNB)正演成果,7B温纳装置(AMNB)反演成果,拟合误差0.177，耗时3.9分,8A温纳装置(ABMN)正演成果,拟合断面,原始断面,拟合误差0.193，耗时6.5分,8B温纳装置(ABMN)反演成果,9A温纳装置(AMBN)正演成果,拟合误差0.231，耗时4.4分,9B温纳装置(AMBN)反演成果,三、反演实例,秦始皇陵2-2线双向单极梯度(NMAMN)反演成果(bmp格式),秦始皇陵2-2线双向单极梯度反演成果(1)及解释,低阻区：推测为墓室的保护材料,次高阻区：推测为墓道位置,高阻区：推测为未塌落、未充水、比较干燥的墓室位置,高阻区：推测为墓道位置,秦始皇陵2-2线双向单极梯度反演成果(2)及解释,秦始皇陵2-2线双向单极梯度反演成果(3)解释,迭代10次拟合误差4.54耗时17.4分,山东省XX金矿区1