第十五讲非线性反演演示教学.pptx

第七章 非线性反演计算地球物理学讲义非线性反演概论地球物理学家研究地球所依据的物性参数不同，方法各异，但就工作程序而言，一般都可分为数据采集，资料处理和反演解释等三个阶段。数据采集就是按照一定的观测系统、一定的测线、测网布置,在现场获得第一手、真实可靠的原始资料。所以数据采集是地球物理工作的基础，是获得高质量地质成果的前提和条件。资料处理的目的是通过各种手段，去粗取精，去伪存真，压制干扰，提高信噪比，使解释人员能从经过处理的资料(异常,或响应) 中、较准确的提取出测区的地质、地球物理信息。资料处理是从原始观测数据到地球物理模型之间的必不可少的手段和过渡阶段。反演解释的目的，用地球物理的术语来说，就是实现从地球物理异常(或响应)到地球物理模型的映射，使解释人员能从经过处理的地球物理资料(异常或响应)中提取出获得最接近真实情况的地质、地球物理模型，圆满的完成提出的地质任务。引自：王家映，地球物理反演问题概述，工程地球物理学报，4（1），2007非线性反演概论资料处理的基础是时间（空间）序列分析。反演解释的基础是反演理论。地球物理观测数据和地球物理模型参数之间的数学表达式，只有线性和非线性两大类。以d表示观测数据向量，m表示模型参数向量，f表示联系d和m的函数或泛函表达式。线性反演问题满足：不满足上面两个条件之一的所有问题都统称为非线性问题。引自：王家映，地球物理反演问题概述，工程地球物理学报，4（1），2007非线性反演概论从模型参数m 到观测数据d 的变换(或映射) ，统称为正演或正问题,并记为：由观测数据d，反推模型参数m的变换(或映射) ，叫反演或反问题。引自：王家映，地球物理反演问题概述，工程地球物理学报，4（1），2007非线性反演概论引自：王家映，地球物理反演问题概述，工程地球物理学报，4（1），2007解决非线性问题有两种办法，一是把非线性问题线性化，按线性问题解决，然后通过迭代的办法，逐次逼近，求得反演问题的近似解；一是不走线性化的路子，而按非线性的办法进行反演。非线性问题线性化的办法，简单、易行，在许多情况下也可以取得较好的结果，但是在目标函数具多极值的情况下，在反演迭代中，容易陷入局部极小，而且反演结果极大的取决于初始模型，也可能使反演出现不稳定，甚至无解。非线性反演概论一般说来，地球物理反演的目标函数都是高次非线性函数，有多个极值。反演时，如初始模型选取不当，不靠近目标函数的全局极小，因而在迭代过程中,只能在局部极小点附近搜索，很难跳出局部极小的阱，这时，只能获得局部最优解；有的反演方法，由于初始模型选取和参数修改的随机性，在搜索的过程中，可能在模型空间全局寻优，其解就是全局最优解。引自：王家映，地球物理反演问题概述，工程地球物理学报，4（1），2007非线性反演概论反演理论家R. L . Parker 在其著名的论文Understanding Inverse Theory中提出，地球物理反演理论必须回答以下4 个问题：1) 解的存在性（Existence）。即给定一组地球物理观测数据： di ， i = 1 , 2 , 3 , ?, M 之后，是否存在一个能拟和观测数据的地球物理模型m？2) 解的非唯一性（Non-uniqueness）。如能求得能拟合观测数据的地球物理模型，解是唯一的还是非唯一的？3) 模型构制（Model Construction）。如何求得能拟和观测数据的一个地球物理模型？4) 解的评价（Appraisal） 。既然解是非唯一，地球物理反演所获得的任一解又有何意义？引自：王家映，地球物理反演问题概述，工程地球物理学报，4（1），2007非线性反演概论理论严格证明，给定一组地球物理观测数据以后，总可以找到一个能拟合它的地球物理模型。由于观测数据的个数并非无限，不构成一个数据的完备群，加之每一个观测数据均有误差，这就决定了地球物理反问题的解是非惟一的。虽然，反演问题的解是非惟一的，但这个非惟一解仍然是有意义的。引自：王家映，地球物理反演问题概述，工程地球物理学报，4（1），20071、物理类 1) 原子跃迁法 (Atomic Transition Algorithm, ATA) 2) 模拟退火法 （Simulated Annealing, SA） 3) 量子退火法 (Quantum SA) 4) 量子遗传算法 (Quantum GA) 2、数学类1) 贝叶斯法 (Bayesian Inversion) 2) 同伦算法 (Homotopy Algorithm) 3) 多尺度非线性反演法 (Multi-Scale Algorithm) 4) 蒙特卡罗法 (Monte-Carlo)3、仿生类1) 粒子群算法 (Particle Swarm Optimization, PSO) 2) 免疫算法