高密度电阻率法2.5维反演演示程序使用说明c要点讲解.doc

高密度电阻率法2.5维反演演示程序 hdrminv_ys_c使用说明书 中国地质大学电法科研组 2005年7月 hdrminv_ys是用以说明高密度电阻率法2.5维反演程序使用方法和反演效果的演示程序。它使用Fortran语言设计，在QuickWin工程环境下编译为执行程序，可以在任何IBM PC或其兼容机上运行。 一、输入文件 反演所需的实际观测条件参数和实测视电阻率数据，采用数据文件输入。演示程序hdrminv_ys，只能用来对给定的三个输入数据文件demo1.dat, demo2.dat和demo3.dat做反演。其中，demo1.dat是一个对称四极（施伦贝格尔）装置的正演计算理论模型数据；demo2.dat是浙江地球物理研究所采用(装置测得的实测数据；而demo3.dat是物化探研究所在秦始皇陵探测中，用双向单极梯度装置的实测结果。虽然我们的高密度电阻率法2.5维反演程序，可以用来反演以下九种常用电极装置的测量数据： 装置类型1((双向单极梯度装置NMAMN，MN = a ，AM = n·a，B((。 其中，a(相邻电极间距，即测点距；n(电极隔离系数。 装置类型2((正向单极梯度装置（三极装置）AMN，MN = a ，AM = n·a， B((。 装置类型3((反向单极梯度装置（三极装置）NMB，MN = a ， MB = n·a， A((。 装置类型4((对称四极(施伦贝格尔)装置AMNB，MN = a，AM = NB = n·a。 装置类型5((偶极-偶极装置ABMN，AB=MN= a，BM = n·a。 装置类型6(( 电位装置（二极装置）AM，AM= n·a， B((和N((。 装置类型7(( (装置（温纳装置）AMNB，AM=MN=NB= n·a。 装置类型8(( (装置ABMN，AB=BM=MN= n·a。 装置类型9(( (装置AMBN，AM=MB=BN= n·a。 但这里提供的演示程序只能反演上述三个数据。这些输入数据文件的共同格式是： LINE Narray,Nd,AES,Nmin,Nmax H(i), i = 1,Nd Is0,dIs,Xs0 Ps0 其中： 第一行LINE为说明测区和测线信息的字符串，字符数不超过60个。 第二行输入观测装置的有关参数。其中： 整型数Narray为电极装置类型号。根据所用电极装置类型，取值为1(9之间。 整型数Nd取值为被反演的实测剖面上布置的电极总数。约定剖面上第一个电极的位置序号为1，以后依此类推，直到最后一个电极的位置序号为Nd。 实型数AES =a取值为观测剖面上相邻电极的间距，也等于测点距，单位米。 整型数 Nmin和Nmax分别为电极装置的最小和最大间隔系数，所谓间隔系数是指供电电极和相邻测量电极间距与a的比值n。 第三行输入布置电极的各测点的高程（实型数组），共Nd个，单位为米。 第四行输入三个整型或实型数：Is0(第1个电极（测点）的编号，整型数；dIs(相邻电极（测点）的编号增量，整型数和Xs0(第1个电极（测点）的X坐标（米），实型数。 从第五行开始以数组形式输入实测视电阻率 Ps0。在第一个演示数据demo1.dat中，Ps0按重庆地质仪器厂所生产仪器的观测数据格式排列；在第二个演示数据demo2.dat中，Ps0按重庆奔腾仪器厂所生产仪器的观测数据格式排列；而在第三个演示数据demo3.dat中，Ps0按物化探研究所的观测数据格式排列。 二、程序的运行 双击执行程序hdrminv_ys.exe便可开始运行本演示程序。运行时，屏幕出现类似WINDOWS的QuickWin用户界面。为便于观察，用户需将其中的“高密度电阻率法2.5维成像系统……”窗口最大化；然后，按其上屏幕提问输入相应数据： 屏幕显示“请键入输入文件名”((用户键入demo1.dat或demo2.dat或demo3.dat，然后按“Enter”键。 2）“需要对实测数据做偶然误差处理吗？（Y/N）”。本程序设计了对实测视电阻率数据做“压制过失误差”和“压制偶然误差”的预处理子段。前者压制单点突变性质的过失误差，效果甚佳，且不会影响有用异常，故列为必做内容。后者采用“移动窗口光滑”算法压制具随机性质的偶然误差，有一定效果，但可能会降低对有用异常的分辨率，故列为供用户选择的处理内容。用户可根据自己的意愿，选用（键入‘Y’或‘y’）或不选用（键入‘N’或‘n’）。 3）屏幕显示“请键入下面内容，这些内容将写入输出文件中：” 1．屏幕显示“完成计算工作的单位:” ((用户键入相应字符，然后按“Enter” 键；也可不输入任何信息，直接按“Enter”键。 2．屏幕显示“计算者:” ((用户键入