地球物理反演理论综述

《地球物理反演理论综述》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、目录摘要1一、反演问题基本概念1二、线性反演问题2三、线性反演问题的求解43.1适定和超定问题43.2欠定问题43.3混定问题4四、非线性反演方法54.1线性化迭代算法54.2最速下降法54.3共轭梯度法64.4遗传算法74.5模拟退货法74.6人工神经网络法8总结：89地球物理反演理论综述摘要在地球物理学中，其核心问题就是如何根据地面上的观测信号推测地球内部与信号有关部分的物理状态。不同的地球物理问题，其数学物理是不同的；同一个物理问题，应为观测方式不同，也会有不同的物理模型。在地球物理学中，大多数的观测数据核模型参数之间是不满足线性关系的。但是在一定

2、近似条件想均可简化或近似简化为线性关系。因此线性反演是地球物理的关键问题。关键词：反演；线性反演；非线性反演一、反演问题基本概念把数据模型中的一个点定义为m，把数据空间中的一个点定义为d，两者的关系可以成：式中，G为模型空间M到数据空间D的一个映射，也称反函数算子，反应了模型m与数据d之间的物理规律从空间映射来看，如果存在一个映射A，使得则A为有数据空间D到模型空间M的映射，即A为G的逆映射，称逆算子。也可以写成我们把给定模型m求解数据d的过程称为正演；把给定数据d求解模型参数m的过程称为繁衍问题。图1.1模型空间域数据空间之间的映射关系示意图9反演问题

3、的研究归纳为四个方面的问题：1)解的存在性：给定数据d，按照物理定律，能否找到满足要求的模型参数m；2)模型构制：若解存在，如何让构制问题的数学模型使得反演问题的解能迅速而准确地确定；3)解的非唯一性：若解存在，其是否唯一；4)解的评价：若解的非唯一性的，如何从非唯一解中获取真实解的信息。关于上述四方面问题的研究就构成了地球物理反演的基本理论。二、线性反演问题为了使问题简单明了而又不失一般性，我们在此讨论一维问题。设有积分方程式中，。在观测数据数目有限的情况下，为便于书写，我们把各参量表示成如下形式由于与线性无关，则式(2-2)可以表示成内积形式我们先用

4、核函数构造另一组正交函数，即以为系数对观测数据作一个线性组合，并令其为，则由此可见，是m在正交基轴上的投影。9这里是Hilbert空间的任意坐标基，可以正交，也可以是不正交。若将其分成两部分，并取为其他任意坐标基则式(2-8)可写成可以证明。因为所以有。第二项是无限维空间中一个向量投影之和，且该向量在M维正交基中的投影为零，则对于问题中的模型m，它可视为零向量，即故即9讨论:1)对于给定的观测数据总是能找到与之对应的数据模型，即解是存在的；2)模型的构制本质上就是对线性无关的核函数实行正交变换，求得相应的新正交坐标基及模型在这个正交坐标基上投影的过程；3

5、)反演问题是在特征解上加以任何零化子向量所得的模型，都可拟合观测数据，所以姐是非唯一的。三、线性反演问题的求解3.1适定和超定问题在线性反演问题中，如果观测数据的个数多于模型参数的个数，我们想得到一组与观测系统之间误差平方和最小的观测数据所对应的模型参数，也就是适用最小二乘法。3.2欠定问题欠定问题中假设方程数比未知的模型参数少，则可以找到很多的最小方差解。即，虽然数据能提供有关模型参数的信息，但是由于信息不足所以不能唯一确定模型参数。为了唯一确定解，可以把某些为引入的信息附加到该问题上，这些附加信息称为先验信息。它是不依赖实际数据使解以某种定量的形式出

6、现。3.3混定问题混定问题是一种混合模式，观测数据个数多于模型参数的个数，但特征值接近或等于零，具有欠定性质。混定问题可以引入求解，取决于预测误差E与模型长度L在极小化过程中的相对重要性，称为阻尼因子或加权因子。9四、非线性反演方法4.1线性化迭代算法我们所遇到的问题中都是求解多远函数的最优化问题，即寻找目标函数极小点或极大点所对应的数学实现过程。迭代算法在给定一个初始点b0后按照一定的规则产生一个新的点b1，如此迭代产生k个点，形成一个序列{bk}，并使得bk不断逼近极值点b*，最终得到最优化问题的解。对于迭代酸度最重要的是极小化序列的收敛性和收敛速度

7、问题。线性化迭代反应过程中包含有分辨率的快速、慢速提高，方差的急增和混沌的相态，数据误差的大小会影响相态的转变速度。在解估计方差急增之前停止迭代就可以得到分辨率高而方差不大的反问题。线性迭代算法的主要问题是可能陷入空间的局部极小区，因此建议把所有可能的模型都作为初始模型输入，然后再反演结果中找出拟合误差最小的解估计。4.2最速下降法最速下降法是一种运用梯度与极值的性质，综合数值计算方法寻找局部极值。基本思想：任一点的负梯度方向是函数值在该点下降最快的方向。将维问题转化为一系列沿负梯度方向用一维搜索方法寻优的问题，利用负梯度作为搜索方向，故称最速下降法。图

8、4.1：最快下降法具体步骤：Step1给定初始点，允许误差，令。Step2计算搜