GeoExpl第讲培训讲义.ppt

电法数据处理与正反演计算 中国地质调查局发展研究中心 2005年10月 北京 一维直流电阻率和极化率测深正反演程序简介 该程序能对二极电位、对称四极和轴向偶极测深装置进行一维层状介质的视电阻率和视极化率正演理论计算或对实测视电阻率和视极化率自动反演层参数。 该程序设计的二极电位、对称四极(MN0)、对称四极(MN＝0)、轴向偶极测深(MN0)和轴向偶极测深(MN＝0)等五种正反演装置实质上仅为三种装置，但是它们的用途不同。 测量电极距MN0与MN＝0在算法上有很大差别，理论上对梯度型装置一般计算电场强度。在野外无法测量电场强度，只能测M，N点电位值△V，用△V/MN来逼近电场强度。本程序依此设计，以备不同使用目的。对MN0，需要M、N两点的电位求差△V，而对MN＝0的装置，必须计算电场强度E，依据ρs＝K*E/I计算。 中国地质调查局发展研究中心 GeoExpl 第十一讲 一维直流电阻率和极化率测深正反演程序 选择当前数据分析模块 选择当前菜单，进入“一维直流电阻率和极化率测深正反演”模块 一维直流电阻率和极化率测深正反演程序之正演计算 用鼠标点击“正演数据输入” ，弹出左下对话框。 在“计算参数选择”选项里选“视电阻率和视极化率”。 在“装置类型选择”选项里选“偶极装置(MN0)”。 下面以“偶极装置(MN0)”为例，介绍正演流程。 对于其他装置的操作步骤类似。 一维直流电阻率和极化率测深正反演程序之正演计算 选择完“装置类型”和“计算参数”，点击确定，弹出如下文件输入对话框； 选择相应的正演数据文件名；弹出“正演参数输入” 对话框； 输入要正演计算的电性层层数，弹出“层参数输入” 对话框； 输入第一层的层参数； 输入第二层的层参数； 输入第三层的层参数； 点击“确定”完成正演计算参数的输入。 由于方法和算法等方 面的局限性程序对层参 数输入的输入有一些规 定： 电阻率参数大于零， 小于10000； 极化率参数大于零， 小于100； 如输入的电性层层数 为n层，则n层的厚度之 和不能大于100米。 一维直流电阻率和极化率测深正反演程序之正演计算 点击“正演计算”按钮，弹出“正演结果文件保存”对话框； 选择“是”，输入相应的正演结果数据文件名；点击“确定” 退出文件保存对话框；此时正演计算很快完成并保存了计算结果； 点击“正演曲线图示”按钮，可以查看正演曲线。 点击“图形保存”按钮，以位图格式保存正演曲线。 点击“退出正演计算”按钮，系统将清除屏幕图形，以方便用户进行下一步操作，至此整个正演过程结束。 一维直流电阻率和极化率测深正反演程序之反演计算 用鼠标点击“反演数据输入” ，弹出左下对话框。在“计算参数选择”选项里选“视电阻率和视极化率”，在“装置类型选择”选项里选“对称四极装置(MN0)”。下面以“对称四极装置(MN0)”为例，介绍反演流程。 依次输入“反演数据文件名”和“工区位置数据文件名”； 依次输入“反演点号”反演的“迭代次数”、“迭代误差”，完成反演数据输入过程。 取对话框中的默认值 一维直流电阻率和极化率测深正反演程序之反演计算 用鼠标点击“实测曲线图示” ，程序将显示实测曲线。 如果实测曲线上有畸变点，可以点击“曲线畸点校正”按钮，纠正或园滑实测视电阻率，这一步不是必须的，若用户认为曲线已足够园滑无需校正则可以跳过； 对于二极和对称四极装置的视电阻率反演，程序提供了简易方便的图上选择初值方法，极大地方便了用户使用。用鼠标点击“实测曲线变换” 可以执行这一功能。 使用方法是：用鼠标左键单击需要校正的实测点，然后在认为图上的正确位置双击鼠标左键，曲线将自动移到正确位置。程序允许进行同一个多次重复点或进行多个畸点校正操作。 一维直流电阻率和极化率测深正反演程序之反演计算 经过或跳过“曲线变换”这一步，用鼠标点击“输入反演初值” 按钮，程序将弹出“层数输入” 对话框。 对于二极和对称四极装置，均可使用变换曲线方法直接从图上选取初始参数。其使用方法为：在变换曲线上依曲线特征判定层数，在变换曲线各层特征点位置点附近点击鼠标右键，点击点位置的纵横坐标即为该层的初始层参数进行分层。经过变换，在变换曲线上，各层分层处电阻率均有较大变化(梯度较大)，选择其稍后处即为层的厚度分界处，如下图所示。 由于目前尚无办法获得偶极装置的初值选择曲线，只能由用户凭经验选择输入初始层参数值。操作步骤为：通过对话框输入层数和初始层参数值。对于偶数装置依次出现层数，各层初始层参数输入对话框，请用户逐一输入即可。 一维直流电阻率和极化率测深正反演程序之反演计算 用鼠标点击“反演计算” 按钮，程序将显示反演结果图形和拟合均方差。若不满意反演结果，可点击“否”按钮，再取初值反演。 如点击“是”按钮，程序继续提示“是否保存反演结果”，点击“