• / 10
  • 下载费用:1 下载币  

地球物理反演理论综述

关 键 词:
物探 地震资料解释 地震处理 反演
资源描述:
目录摘要 .................................................................................................................................................演问题基本概念 .................................................................................................................性反演问题 .........................................................................................................................性反演问题的求解 ............................................................................................................. 适定和超定问题 .............................................................................................................. 欠定问题 .......................................................................................................................... 混定问题 ..........................................................................................................................线性反演方法 ..................................................................................................................... 线性化迭代算法 .............................................................................................................. 最速下降法 ...................................................................................................................... 共轭梯度法 ..................................................................................................................... 遗传算法 .......................................................................................................................... 模拟退货法 ...................................................................................................................... 人工神经网络法 ...........................................................................................................................................................................................................................................................核心问题就是如何根据地面上的观测信号推测地球内部与信号有关部分的物理状态。不同的地球物理问题,其数学物理是不同的;同一个物理问题,应为观测方式不同,也会有不同的物理模型。在地球物理学中,大多数的观测数据核模型参数之间是不满足线性关系的。但是在一定近似条件想均可简化或近似简化为线性关系。因此线性反演是地球物理的关键问题。关键词: 反演;线性反演;非线性反演一、反演问题基本概念把数据模型中的一个点定义为 m,把数据空间中的一个点定义为 d,两者的关系可以成: 中,G 为模型空间 M 到数据空间 D 的一个映射,也称反函数算子,反应了模型 m 与数据 d 之间的物理规律从空间映射来看,如果存在一个映射 A,使得A则 A 为有数据空间 D 到模型空间 M 的映射,即 A 为 G 的逆映射,称逆算子。也可以写成 我们把给定模型 m 求解数据 d 的过程称为正演;把给定数据 d 求解模型参数 m 的过程称为繁衍问题。m 型空间域数据空间之间的映射关系示意图2反演问题的研究归纳为四个方面的问题:1) 解的存在性:给定数据 d,按照物理定律,能否找到满足要求的模型参数 m;2) 模型构制:若解存在,如何让构制问题的数学模型使得反演问题的解能迅速而准确地确定;3) 解的非唯一性:若解存在,其是否唯一;4) 解的评价:若解的非唯一性的,如何从非唯一解中获取真实解的信息。关于上述四方面问题的研究就构成了地球物理反演的基本理论。二、线性反演问题为了使问题简单明了而又不失一般性,我们在此讨论一维问题。设有积分方程式中, 。在观测数据数目有限的情况下,为便于书写,我们把各参量表示成如下形式x, m,2,1由于 与 线性无关,则式(2以表示成内积形式,j,,我们先用核函数 构造另一组正交函数,即jM,2,1以 为系数对观测数据 作一个线性组合,并令其为 ,则,,111 由此可见, 是 m 在正交基 轴上的投影。11这里 是 间的任意坐标基,可以正交,也可以是不正交。若将其分成两部分,并取kM,2,1为其他任意坐标基 k则式(2写成可以证明 。11 因为11,1l 10所以有 。第二项 是无限维空间中一个向量投影之和,且该1向量在 M 维正交基 中的投影为零,则对于问题中的模型 m,它可视为零向量,即10故10即412110,,讨论:1) 对于给定的观测数据总是能找到与之对应的数据模型,即解是存在的;2) 模型的构制本质上就是对线性无关的核函数实行正交变换,求得相应的新正交坐标基及模型在这个正交坐标基上投影的过程;3) 反演问题是在特征解上加以任何零化子向量所得的模型,都可拟合观测数据,所以姐是非唯一的。 三、果观测数据的个数多于模型参数的个数,我们想得到一组与观测系统之间误差平方和最小的观测数据所对应的模型参数,也就是适用最小二乘法。可以找到很多的最小方差解。即,虽然数据能提供有关模型参数的信息,但是由于信息不足所以不能唯一确定模型参数。为了唯一确定解,可以把某些为引入的信息附加到该问题上,这些附加信息称为先验信息。它是不依赖实际数据使解以某种定量的形式出现。测数据个数多于模型参数的个数,但特征值5接近或等于零,具有欠定性质。混定问题可以引入 求解, 取决于预测误差22E 与模型长度 L 在极小化过程中的相对重要性,称为阻尼因子或加权因子。四、寻找目标函数极小点或极大点所对应的数学实现过程。迭代算法在给定一个初始点 按照一定的规则产生一个新的点 如此迭代产生 k 个点,形成一个序列{b k} ,并使得 断逼近极值点 b*,最终得到最优化问题的解。对于迭代酸度最重要的是极小化序列的收敛性和收敛速度问题。线性化迭代反应过程中包含有分辨率的快速、慢速提高,方差的急增和混沌的相态,数据误差的大小会影响相态的转变速度。在解估计方差急增之前停止迭代就可以得到分辨率高而方差不大的反问题。线性迭代算法的主要问题是可能陷入空间的局部极小区,因此建议把所有可能的模型都作为初始模型输入,然后再反演结果中找出拟合误差最小的解估计。合数值计算方法寻找局部极值。基本思想:任一点的负梯度方向是函数值在该点下降最快的方向。将 用负梯度作为搜索方向,故称最速下降法。6图 快下降法具体步骤: 给定初始点 0允许误差 0,令 1k。 计算搜索方向 ()。 若 kA,则 则,求解最优步长 k,使得()。 令11最速下降方向是反映了目标函数的局部性质,它只是局部目标函数值下降最快的方向。轭梯度法基本思想:将共轭性与最速下降法相结合利用已知迭代点的梯度方向构造一组共轭方向,并沿此方向搜索,求出函数的极小值。例如: 12其中, 0 是对称正定矩阵具体步骤: 取初始点 (0)x,取第一次搜索方向为 (0)(0)。 设已求得 (1)k,若 (1),令 (1)则下一个搜索方向 (1)()1由于 (1)k与 ()共轭,所以给(1)两边同时乘以 :()(1)()()()10解得:()() 搜索步长的确定,已知迭代点 ()搜索方向 ()定步长 k,即:()()7解得:()共轭梯度法是对最速下降法的一种改进,减少了迭代次数从而提高了程序运行效率。设置迭代参数确定进化代数 随机产生规模为 且满足约束条件的群体,0对群体 的个体进行评价,如果个体 满足约束条件,则随机生成一个满足约束条件的个体来替换,并保存当前最好的个体采用比例选择算子对当前群体进行选择操作,选择群体规模为 的对群体中的个体进行两两随机配对。以概率 交换部分基因。对群体中个体的每一个基因以概率 进行变异,未变异的个体直变异后的新群体取代上一代群体,并对当前群体 中的个体进行1若当前群体中最大适应值大于上一代群体中的最大适应值,则当前群体中的最好个体取代上一代保留下的最好个体,否则不取代。条件判断,如满足条件,则停止迭代。否则, ,转回1传算法的核心是由“繁殖” , “杂交” , “变异”三步组成的转移过程。同时计算繁殖杂交概率与变异概率用于控制迭代过程并使数据进一步拟合。热力学的理论套用到统计学上,将搜寻空间内每一点想像成空气内的分子;分子的能量,就是它本身的动能;而搜寻空间内的每一点,也像空气分子一样带有“能量” ,以表示该点对命题的合适程度。演算法先以搜寻空间内一个任意点作起始:每一步先选择一个“邻居” ,然后再计算从现8有位置到达“ 邻居” 的概率。模拟退火的基本思想:(1) 初始化:初始温度 T(充分大),初始解状态 S(是算法迭代的起点), 每个 T 值的迭代次数 L(2) 对 k=1, ……,L 做第(3)至第(6)步:(3) 产生新解 S′(4) 计算增量 Δt′=C(S′)),其中 C(S)为评价函数(5) 若 Δt′0,然后转第 2 步。模拟退火法是用于求解非线性地球物理反问题的一种好方法,对于模型参数较小的情况下计算很方便,但实际应用上其反演结果还是依赖于初始参数选择和降温方式选择等问题。种映射既可以用模式识别的方式实现、也可以用数字的方式定量的实现。人工神经网络在反演中的应用,大多数是基于模式识别的原理。和常规反演一样,基于 传理论建立起来的反演方法 ,也要有目标函数。此时的目标函数定义为所有输入模式对上输出单元之希望输出与实际输出之误差的平方和。和常规的反演不同,基于 传理论 ,可将反演过程分为学习(训练)和反演( 测试 )两个阶段。先学习后反演。学习就是对网络进行训练。神经网络的基本特征:巨量并行,运用多方面知识和经验同事并发,迅速做出解答;信息存储和信息处理和在一起;自组织,自学习的功能。总结:目前地球物理反演的主要应用是在三维结构的研究,地球物理迭代反应的9线性化迭代方法具有较好的数学物理基础随着计算机软硬件的快速发展,还将会出现更多的新变化。遗传算法只要繁殖杂交概率和变异概率选择合适总能收敛到较优化的解估计。但是某些地球物理反问题不适宜会出现传代太快,不能收敛全局的极值。所以要选择适宜的方法解决反问题。人工神经网络算法是发展比较迅速的信息处理技术,但在物探方面还局限于模式识别,开发人工神经网络的潜能石下一阶段反演研究的重要方向。参考文献[1]]质出版社, 1997.[2]]国地质大学出版社 3]]等教育出版社 2002.[4]周明, 孙树栋 . 遗传算法原理及其应用[M]. 北京: 国防工业出, 1999[5]]995.[6]]995[7] 侯格贤, 吴成柯. 遗传算法的性能分析 [J]. 控制与决策 , 1999.[8]秦浩宇. 遗传算法的改进研究及其在酵母扩培系统中的应用[D]. 中南大学, 2010.[9]]992.[10]栾文贵 ]11]杨文采 J]995.
展开阅读全文
  石油文库所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

关于本文
本文标题:地球物理反演理论综述
链接地址:http://www.oilwenku.com/p-63032.html
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2016-2020 石油文库网站版权所有
经营许可证编号:川B2-20120048,ICP备案号:蜀ICP备11026253号-10号
收起
展开