电阻率／激发极化法测深数据的一维最优化反演方法.pdfVIP

维普资讯 第 19卷 第4期 桂 林 工 学 院 学 报 Ⅷ l9 No4 JOURNALOFGUILIN INSTIT15I~ OFTECHNOLOGY Oct 1999 1999年 10月 电阻率／激发极化法测深数据的一维 弓 ／一 最优化反演方法 阮 百 尧 (桂韩工学院贾菰写 j 541004) 摘 要 采用差分法求偏导数矩阵．同时对各地层中的电阻率和厚度参数进行几何修 改．以及用奇异值分解法解线性方程组，介绍了一种收敛速度快、对初始模型要求不 高、反演效果好的电阻率测深数据的一维反演方法。在此基础上，保持各层厚度参数 不变，用塞格尔理论．完成对极化率的反演。 关键词 !!墨 避垄壑i鸯 些些星 法 奇异值分解 分类号 1631．324；P631，322 1 引 言 电阻率测深数据的反演方法可归结为使下面 目标 函数趋于极小： = ∑ 一 (M)l。 (1) 这里，NS是供电极距数； 是第 个极距的实测视电阻率；M (M ．j=1．2，…，NM)是模型参数 (地层的电阻率或厚度)．NM是预测模型参数个数； 是由预测模型正演计算所得的第 个极距 ， 的理论视电阻率 ；n是范数，当n=2时即为最小二乘法。 将在预测模型 M 处展开，并忽略二次项 以上的项 ，得到 = =I∑；L0ff一一阳阳一一∑，lu々嘶’。【 ㈤(2) 上式要趋于极小，则对于各供电极距 (i=1．2．…，NS)要满足下面的线性方程 ： 器= 阳 (3) 解上述线性方程组 ．就可得到预钡I模型 M 的修改量△M ．于是可得到新的预测模型。计算 新模型的理论视电阻率，与实测视电阻率进行对 比，如果精度满足要求，则新的预测模型即为 反演结果 ；否则重新展开，计算模型修改量．直到精度满足要求。 上述反演过程中，有 3个fⅡ1题需要解决 ：第一 ．偏导数的求解 ；第二，模型参数的处理。 1999—04一∞ 收稿。 作者简介：阮百尧．男．1962年出生．博士．教授．应用地球物理专业．￡一Ⅱmd—rby@ gliteed帅 *国家 自然科学基金资助项 目 和广西自然科学基金匹配项 目成果之一。 维普资讯 桂 林 工 学 院 学 报 1999妊 模型参数中有不 同量纲的电阻率和厚度．尤其对电阻率参数，变化范围很大，如果直接由上述 方法求解．不但会导致方程 (3)严重病态，而且 电阻率和厚度参数的修改也不会正确，从而 会导致反演方法不收敛 ；第三，线性方程组 (3)的求解。线性方程组 (3)往往是不相容方程 组，求解方法不同 收敛速度和收敛性也不同。前人也 曾对这些问题进行过探讨，并提出了许 多解决方法 ，本文将在此基础上，首先介绍对这 3个 问题新 的处理方法．然后再根据塞格 尔理论，完成对极化率的反演。 2 偏导数的计算 电阻率测深 中。视电阻率对模型参数 的偏导 卸 ／ 通常是从数值滤波正演过程 中直接给 出 。本文考虑到对模型参数的生睫 ．用差分方法来求计算．取 △f =0．1M ．则 ／a = [ ( ，M2，…，1，1 ，…，Mm)一pe(M,，M2，…，M ，…， M)]／(0．IMp)