层析反演中的正则化方法研究

合集下载

四川盆地西北部深层低幅度构造建模及成像

四川盆地西北部深层低幅度构造建模及成像

2021年2月第56卷第1期•偏移成像•文章编号:1000-7210(2021)01-0109-09四川盆地西北部深层低幅度构造建模及成像王艳香*苏勤乐幸福张军舵刘威袁焕(中国石油勘探开发研究院西北分院,甘肃兰州730020)王艳香,苏勤,乐幸福,张军舵,刘威,袁焕.四川盆地西北部深层低幅度构造建模及成像.石油地球物理勘探,2021,56(1):109-117.摘要为落实四川盆地西北部龙门山前带构造圈闭细节,开展了低幅度构造成像、网格层析速度建模和构造物理模拟实验研究。

龙门山前带地表地下构造复杂,采用构造正演模拟的方法研究其动力学成因;利用自适应微测井约束初至层析静校正技术解决该区影响低幅度构造成像的静校正问题:应用各向异性动校正技术解决在大炮检距处产生的动校过量问题,提高了浅部陡倾地层及深部低幅度构造的速度分析精度;将约束层析的近地表速度模型和沿层层析获得的中深层速度模型融合,充分利用钻丼数据等多信息指导初始速度模型的建立,利用网格层析优化速度模型并偏移,为地震资料解释提供了准确可靠的数据。

关键词四川盆地西北部深部低幅度构造自适应撖测丼约束层析静校正层控网格层析构造物理模拟中图分类号:P631文献标识码:A doi:10. 13810/j. cnki. issn. 1000-7210. 2021. 01. 013〇引言低幅度构造是指构造相对平缓、闭合幅度只有 10〜20m的地质体1_14],又被称为小幅度构造、微幅 度构造等。

近年来,低幅度构造油气藏逐渐被重视,成为油气勘探的一个重要接替领域,因此得到了人 们的广泛关注[3]。

由于构造幅度较低,因此地震资 料上表现为反射同相轴平直而变化幅度很小,不易 被识别。

低幅度构造的解释精度取决于地震数据的 采集、处理和解释等各个环节,处理人员就低幅度构 造圈闭地震资料处理中的关键点、横向速度变化对 构造成像影响进行了大量研究[46],解释人员在低幅 度构造解释技术、识别技术、储层横向预测等方面也 取得了较大进展[78]。

地球物理反演理论

地球物理反演理论

地球物理反演理论一、解释下列概念1.分辨矩阵数据分辨矩阵描述了使用估计的模型参数得到的数据预测值与数据观测值的拟合程度,可以表示为[][]pre est g obs g obs obs d Gm G G d GG d Nd --====,其中,方阵g N GG -=称为数据分辨矩阵。

它不是数据的函数, 而仅仅是数据核G (它体现了模型及实验的几何特征)以及对问题所施加的任何先验信息的函数。

模型分辨矩阵是数据核和对问题所附加的先验信息的函数,与数据的真实值无关,可以表示为()()est g obs g true g ture ture m G d G Gm G G m Rm ---====,其中R 称为模型分辨矩阵。

2.协方差模型参数的协方差取决于数据的协方差以及由数据误差映射成模型参数误差的方式。

其映射只是数据核和其广义逆的函数, 而与数据本身无关。

在地球物理反演问题中,许多问题属于混定形式。

在这种情况下,既要保证模型参数的高分辨率, 又要得到很小的模型协方差是不可能的,两者不可兼得,只 有采取折衷的办法。

可以通过选择一个使分辨率展布与方差大小加权之和取极小的广义逆来研究这一问题:()(1)(cov )u aspread R size m α+-如果令加权参数α接近1,那么广义逆的模型分辨矩阵将具有很小的展布,但是模型参数将具有很大的方差。

而如果令α接近0,那么模型参数将具有相对较小的方差, 但是其分辨率将具有很大的展布。

3.适定与不适定问题适定问题是指满足下列三个要求的问题:①解是存在的;②解是惟一的;③解连续依赖于定解条件。

这三个要求中,只要有一个不满足,则称之为不适定问题4.正则化用一组与原不适定问题相“邻近”的适定问题的解去逼近原问题的解,这种方法称为正则化方法。

对于方程c Gm d =,若其是不稳定的,则可以表述为()T T c G G I m G d α+=,其中α称为正则参数,其正则解为1()T T c m G G I G d α-=+。

地球物理反演中的正则化技术分析

地球物理反演中的正则化技术分析

地球物理反演中的正则化技术分析地球物理反演是一种通过观测地球上各种现象和数据,来推断地球内部结构和物质分布的方法。

在地球物理反演中,由于观测数据的不完整性和不精确性,常常需要借助正则化技术来提高反演结果的可靠性和准确性。

正则化技术是一种以一定规则限制解的优化方法。

通过在反演过程中引入附加信息或者假设,正则化技术可以帮助减小反演问题的不确定性,提高解的稳定性和可靠性。

在地球物理反演中,正则化技术有多种应用。

下面将介绍几种常见的正则化技术,并对其进行分析和比较。

1. Tikhonov正则化Tikhonov正则化是一种基本的正则化技术,它通过在目标函数中加入一个范数约束来限制解的空间。

常见的约束可以是L1范数和L2范数。

L1范数可以使解具有稀疏性,即解中的大部分分量为零,适用于具有稀疏特性的反演问题。

而L2范数可以使解具有平滑性,适用于具有平滑特性的反演问题。

2. 主成分分析正则化主成分分析正则化是一种通过将反演问题映射到低维空间来减小问题的维度的正则化技术。

它可以通过选择重要的主成分来实现数据降维,从而减少反演问题的不确定性。

主成分分析正则化在处理高维数据时可以提高反演的效率和精度。

3. 奇异值正则化奇异值正则化是一种基于奇异值分解的正则化技术。

通过对反演问题进行奇异值分解,可以将问题分解为多个低维子问题,从而减小高维问题的不确定性。

奇异值正则化适用于非线性反演问题,可以提高反演结果的稳定性和可靠性。

4. 稀疏表示正则化稀疏表示正则化是一种基于稀疏表示理论的正则化技术。

它通过将反演问题转化为对系数矩阵的优化问题,并引入L1范数约束,使得解具有稀疏性。

稀疏表示正则化适用于信号重构和图像恢复等问题,并在地震勘探和地球成像中有广泛应用。

在选择正则化技术时,需要考虑问题的特性和数据的特点。

不同的正则化技术适用于不同的问题,并且各自具有一些优势和限制。

因此,根据问题的具体要求和数据的特征,选择合适的正则化技术可以提高反演结果的可靠性和准确性。

速度层析成像正反演技术

速度层析成像正反演技术

速度层析成像正反演技术杨晓弘;何继善;谢冬琪【摘要】利用层析成像的基本原理对地质异常体进行速度层析成像正、反演研究,计算结果可以确定出地质模型中速度异常体的大小、位置和慢度值等物性参数.首先通过建立正演模型,将含有异常体的模型区域剖分成规则的矩形网格,在区域两边的钻孔中分别安置发射和接受装置,采用多点发射、多点接收的方法获得已知地质条件下的每条速度射线的走时值,而后根据线性回归的原则对所得的时间数据加入20%的噪声,采用医学中成熟的ART算法进行模型的反演计算,计算结果表明,反演结果与正演模型十分逼近,计算精度和速度都能达到满意的效果.【期刊名称】《物探与化探》【年(卷),期】2009(033)002【总页数】4页(P217-219,223)【关键词】速度层析成像;正反演;地质模型;慢度【作者】杨晓弘;何继善;谢冬琪【作者单位】中南大学,信息物理工程学院,湖南,长沙,410083;中南大学,信息物理工程学院,湖南,长沙,410083;中南大学,信息物理工程学院,湖南,长沙,410083【正文语种】中文【中图分类】P631.4层析成像技术(computerized tomography)是一种用数学方法把许多射线路径得到的信息组合成射线在其中传播的介质图像的技术,医学上利用层析成像技术可清晰地重现身体内部的构造[1-2]。

目前,国内外关于层析成像方法的研究已经取得了丰硕的成果,地球物理中的层析成像技术则是利用对象的各种物理性质和物性参数来重建地质体的内部图像,为其他的地球物理资料处理方法提供精确的速度模型[3]。

自20世纪70年代以来,速度层析成像技术在油气田勘探开发、地质灾害预报、无损检测等领域中的应用日益广泛,已从直射线层析成像发展到弯曲射线层析成像,目前存在的主要问题仍然是计算效率和精度相对较低[4-6]。

根据Backus-Gilbert理论[7-8],分辨率和精度之间存在折中关系,在增加射线数目的情况下,正演和反演的关键就变成了求解一个大型稀疏线性方程组的问题。

地球物理反演的理论基础与方法研究

地球物理反演的理论基础与方法研究

地球物理反演的理论基础与方法研究地球物理反演是研究地球内部结构和性质的一种重要方法。

它通过利用地球表面或近地表的观测数据,推断地球内部的物理参数分布。

地球物理反演的理论基础与方法是支撑反演技术的关键,下面将重点介绍地球物理反演的理论基础和常用方法。

1. 理论基础地球物理反演的理论基础主要涉及地球内部物理参数与观测数据之间的关系。

常用的理论基础包括地球物理学原理、数学方法、统计学方法等。

(1)地球物理学原理:地球物理学原理是地球物理反演的基础。

它包括重力学、磁力学、地震学、电磁学等学科的原理,通过分析这些物理过程的规律,可以推断地下介质的性质和结构。

(2)数学方法:数学方法是地球物理反演中处理观测数据和求解反演问题的重要工具。

常用的数学方法包括线性与非线性最小二乘方法、正则化方法、优化算法等。

这些方法可以将观测数据与地下介质的参数之间建立数学模型,通过数值计算来求解最优解。

(3)统计学方法:统计学方法在地球物理反演中的应用越来越广泛。

它可以解决一些非唯一性问题,通过统计分析建立多个可能的模型,提供多个可能的解释。

统计学方法还可以对反演结果进行可靠性评估,提供不确定性估计。

2. 常用方法地球物理反演的方法多种多样,根据不同的物理量和观测方法可以分为地震反演、重磁反演、电磁反演等。

(1)地震反演:地震反演是利用地震波在地下传播的特性,通过分析地震波的传播速度、振幅等信息,推断地下介质的密度、泊松比、剪切模量等物理参数。

常用的地震反演方法有全波形反演、层析成像、声波全息等。

(2)重磁反演:重磁反演是利用地球重力场和地球磁场的观测数据,推断地下介质的密度、磁化率等物理参数。

常用的重磁反演方法有静态反演、动态反演、傅立叶反演等。

(3)电磁反演:电磁反演是利用电磁场的观测数据,推断地下介质的电导率、介电常数等物理参数。

常用的电磁反演方法有研究地电场、研究磁场、研究电磁场构造等。

此外,还有多物理场反演、岩石物理反演、非线性反演等方法,可以根据不同的需求和观测数据选择合适的方法进行反演。

强偏折场的光偏折层析重建技术

强偏折场的光偏折层析重建技术

t ele a nni a n es n ii lry t cn to s d s n d fr cl l in o a r etr n o raz o l er ivri .A dgt a r ig meh d i ei e o a ua o f ry tjcoy i i n o a a g c t a
第 3 卷第 6 7 期
21 0 0年 6月
光 电工 程
Op o Elc r n cEn i e rn t — e to i g n e i g
V l3 . . 0 _ - 7 NO 6 J n , 01 ue2 0
文章编 号 :1 0 9 0
结合光线追迹算法提 出基于曲线路径反 演的正则化修正重建技术。 该技术克服 了现有重建技术 中光路 直进的假设 ,
基 于几何光学原理设 计 了光线通过 强偏折场 的积分路径追踪算法。使 用代数 重建技术和正则化重建技术相 结合 的
思想 ,通过偏折 角投 影数据 的迭代修正 实现非线性反演。使用该技 术对 强偏折场进行数值重建 ,证 实了其在有 限
ZHANG n , SONG ng , HE Bi Ya An. hi z 2
( . ol e fEe t meh ncl n i ei , ig a nvri cec n c n l y Qi d o2 6 6 , 1 C lg lc o ca i gn r g Qn doU i syo i e dT h oo , n a 6 0 e o r aE e n e t fS n a e g g 1
Ab t a t sr c :W h n t mo r p y i e f r d wi  ̄o g y r fa t g f l , e a t n ma a s u c e tb n i g o e e o g a h sp r o me t s n l e c i e d r f c i y c u e s f in e d n f h h r n i r o i t p o i g r y . d n r e o s u t n a g r h r b n s Or i ay r c n t ci l o t m, a e n t ea s mp i n o r ih a s d e o il c u ae r s l . a r o i b s d o su t fs a g t y , o sn t e d a c r t e u t h o t r y s A e a p o c o c r e -a o g a h cr c n tu t n o r n l e a t g r fa t e i d x d sr u i n i r p s d n w p r a h t u v d- y t mo r p i e o sr c i f to g y r f c i er c i ・ e it b to sp o o e r o s r n v ・ n i

地震波阻抗反演方法综述

地震波阻抗反演方法综述

地震波阻抗反演方法综述地震波阻抗反演方法可以分为直接方法和间接方法。

直接方法是指直接根据地震波观测数据反演地下结构的方法,常见的直接方法有全波形反演。

间接方法是指通过建立模型和计算地震波传播路径来反演地下结构的方法,常见的间接方法有层析成像、正则化反演和遗传算法等。

全波形反演是一种直接方法,它利用完整的地震波观测数据来反演地下结构。

全波形反演的核心是通过比较实际观测数据和模拟数据的差异来优化模型参数。

全波形反演可以获取高分辨率的地下结构信息,但由于计算复杂度高、非线性程度强等因素,全波形反演面临着一些挑战。

层析成像是一种常用的间接方法,它通过在空间上离散化模型并计算地震波在传播路径上的传播时间与振幅的差异来重建地下结构。

层析成像的原理是建立了地震波传播路径上的散射模型,通过优化模型参数使计算值与实际观测值吻合。

层析成像具有分辨率高、计算效率高等优点,适用于复杂地质环境的反演。

正则化反演是一种常用的间接方法,它通过在反演过程中引入先验信息来约束模型的解。

正则化反演的核心是将反演问题构建成最优化问题,并添加正则化项以保证解的稳定性。

常见的正则化方法有Tikhonov正则化、L1正则化和全变差正则化等。

正则化反演可以提高反演结果的稳定性,但其分辨率相对较低。

遗传算法是一种通过模拟进化过程来求解最优问题的优化方法。

在地震波阻抗反演中,遗传算法可通过定义模型参数的染色体编码、适应度函数以及遗传操作等步骤来最优解。

遗传算法能够全局,适用于非线性、多峰反演问题,但也存在计算复杂度高、空间维度大等问题。

除了上述的方法,还有一些其他地震波阻抗反演方法,如基于人工神经网络的反演、基于模糊数学的反演等。

这些方法各有特点,适用于不同的反演问题。

地震波阻抗反演方法在地球物理勘探、地震灾害预测等领域有着广泛的应用。

不同的反演方法具有不同的优点和缺点,需要根据具体问题的需求选择合适的方法。

未来地震波阻抗反演方法的发展方向将是提高反演的分辨率和稳定性,减少计算复杂度,开展多物理场的耦合反演研究。

全波形反演正则化方法对比

全波形反演正则化方法对比

全波形反演正则化方法对比
李昕洁;王维红;郭雪豹;张庭俊
【期刊名称】《石油地球物理勘探》
【年(卷),期】2022(57)1
【摘要】正则化是缓解反演不适定性、约束解特征的重要方式。

Tikhonov正则化、全变分(TV)正则化是全波形反演中常用的两种正则化方法,分别具有压制高波数和
保护地层边缘的特点。

双参数整形正则化、混合双参数正则化和稀疏结构约束正则化是在二者的基础上发展而来,并具备各自优势。

为系统论证不同正则化方法特点,
对五种正则化方法的全波形反演进行了对比分析。

背斜—超覆模型、Marmousi
模型测试表明,不同的正则化方法均对反演结果有不同程度的改进。

双参数整形正
则化兼具Tikhonov正则化和TV正则化的优势,并可有效提高深部反演精度。

混合双参数正则化能进一步提高浅层反演精度。

相较于其他方法,稀疏结构约束正则化
无论在地层连续性,还是在边缘结构的刻画上均有明显优势。

【总页数】13页(P129-139)
【作者】李昕洁;王维红;郭雪豹;张庭俊
【作者单位】东北石油大学地球科学学院
【正文语种】中文
【中图分类】P631
【相关文献】
1.井间地震资料全变差正则化波形反演
2.基于改进全变差正则化的GPR多尺度全波形双参数同步反演
3.隧道环境下频率域声波全波形反演优化方法对比
4.多尺度全波形反演的正则化思想及速度模型的多尺度分解
5.弹性波全波形反演中的四种优化方法对比
因版权原因,仅展示原文概要,查看原文内容请购买。

Tikhonov正则化参数的选取及两类反问题的研究的开题报告

Tikhonov正则化参数的选取及两类反问题的研究的开题报告

Tikhonov正则化参数的选取及两类反问题的研究的开题报告题目:Tikhonov正则化参数的选取及两类反问题的研究一、研究背景和意义:随着科学技术的进步,反问题研究成为了最热门的研究领域之一。

反问题的研究涉及到的学科领域非常广泛,其中数学、物理和工程等领域是最为重要的。

反问题包括了许多子领域,如参数反问题、区域反问题、混合反问题等等。

其中参数反问题是最为基础和重要的子领域之一。

Tikhonov正则化方法在参数反问题中得到了广泛应用,因为它可以通过降低噪声波动和提高解的光滑性来改进问题的稳定性。

然而,在应用Tikhonov正则化方法时,如何选取正则化参数是一个非常重要的问题,因为不同的正则化参数会影响到结果的精度和稳定性。

此外,不同类型的反问题需要对正则化参数作出不同的选择,这也是一个需要进一步探究的问题。

因此,我们需要对Tikhonov正则化参数的选取以及在不同类型的反问题中的应用进行深入的研究。

二、研究内容和目标:本文将主要研究Tikhonov正则化参数的选取方法,探讨其在参数反问题和区域反问题中的应用。

具体研究内容包括以下几个方面:1. 对Tikhonov正则化方法的优化算法进行研究,包括最小二乘方法、正交匹配迭代算法等。

2. 针对参数反问题,研究不同类型的Tikhonov正则化方法与对应的正则化参数的选取方法,并比较其性能和精度。

3. 针对区域反问题,研究不同类型的Tikhonov正则化方法与对应的正则化参数的选取方法,并比较其性能和精度。

4. 开发相应的计算程序,实现研究结果的数值验证和实际应用。

通过以上研究,本文旨在实现以下目标:1. 系统性地总结不同类型的Tikhonov正则化方法与对应的正则化参数的选取方法,并探讨其适用范围和局限性。

2. 比较不同类型的Tikhonov正则化方法及其选取的正则化参数在参数反问题和区域反问题中的应用效果,提出相应改进措施,提高解的稳定性和精度。

3. 开发相应的计算程序,实现研究结果的数值验证和实际应用,为相关领域的研究提供参考。

地球物理学研究中的反演方法

地球物理学研究中的反演方法

地球物理学研究中的反演方法地球物理学研究是一门涉及地球内部结构和物质组成的学科,从事这项研究需要掌握一定的物理知识和专业技能,而反演方法则是地球物理学研究的重要工具之一。

反演方法是指根据测量得到的地球物理数据,推算出地球内部结构和物质组成的过程,是一种重要的物理数学分析手段。

在地球物理学研究中,常用的反演方法包括地震层析成像、电磁场反演、地磁场反演、重力反演等。

本文将就地球物理学研究中的反演方法进行阐述。

一、地震层析成像方法地震层析成像方法是一种通过地震波传播路径来推断地球的三维结构的方法。

地震波可以沿着曲折的路径穿过地球中的各种物质,而当地震波沿着不同的路径传播时,它们会受到不同的影响,如反射、折射、散射、压缩等,根据这些影响就可以推断地球内部横截面的结构。

地震层析成像方法主要包括射线追踪、全波形反演和双向波路径方法等。

二、电磁场反演方法电磁场反演方法是一种通过测量地球表面或近表面电磁场的变化来推断地下物质电导率的分布状况的方法。

电磁场反演方法主要包括电阻率层析成像、磁化率层析成像、电场、磁场重力反演等。

三、地磁场反演方法地磁场反演方法是一种通过测量地球表面或近表面磁场的变化来推断地下物质磁性的分布状况的方法。

地磁场反演方法主要包括磁性层析成像、重力反演等。

四、重力反演方法重力反演方法是一种通过测量地球表面或近表面重力值的变化来推断地下物质密度分布状况的方法。

重力反演方法主要包括引力异常反演、引力梯度反演、重力谱反演等。

总之,地球物理学研究中的反演方法是一个复杂的科学体系,需要将物理学、数学、计算机科学等多个学科融合在一起,才能够高效地推算出地球内部结构的分布情况。

虽然反演方法在地球物理学研究中起到了重要的作用,但是它也存在一定的局限性。

例如测量误差、相位问题、非唯一性等问题都会影响到反演结果的准确性。

因此,在进行地球物理学研究的过程中,需要结合多种反演方法,将不同的地球物理数据综合起来,才能获得更加准确和完整的地球内部结构信息,为地球科学研究提供更加可靠的数据支撑。

正则化和反问题

正则化和反问题

正则化和反问题正则化(regularization)在线性代数理论中,不适定问题通常是由一组线性代数方程定义的,而且这组方程组通常来源于有着很大的条件数的不适定反问题。

大条件数意味着舍入误差或其它误差会严重地影响问题的结果。

反问题有两种形式。

最普遍的形式是已知系统和输出求输入,另一种系统未知的情况通常也被视为反问题。

许多反问题很难被解决,但是其他反问题却很容易得到答案。

显然,易于解决的问题不会比很难解决的问题更能引起人们的兴趣,我们直接解决它们就可以了。

那些很难被解决的问题则被称为不适定的。

一个不适定问题通常是病态的,并且不论是简单地还是复杂地改变问题本身的形式都不会显著地改善病态问题。

另一方面,病态问题不一定是不适定的,因为通过改变问题的形式往往可以改善病态问题。

在严格的数学意义上,我们通常不可能对不适定问题进行求解并得到准确解答。

然而,通过使用我们的先验知识,我们通常有希望能够得到一个接近准确解答的答案。

求解不适定问题的普遍方法是:用一族与原不适定问题相"邻近"的适定问题的解去逼近原问题的解,这种方法称为正则化方法。

如何建立有效的正则化方法是反问题领域中不适定问题研究的重要内容。

通常的正则化方法有基于变分原理的Tikhonov正则化、各种迭代方法以及其它的一些改进方法,这些方法都是求解不适定问题的有效方法,在各类反问题的研究中被广泛采用,并得到深入研究。

正则化:Normalization,代数几何中的一个概念。

通俗来说,就是给平面不可约代数曲线以某种形式的全纯参数表示。

即对于PC^2中的不可约代数曲线C,寻找一个紧Riemann面C*和一个全纯映射σ:C*→PC^2,使得σ(C*)=C严格的定义如下:设C是不可约平面代数曲线,S是C的奇点的集合。

如果存在紧Riemann 面C*及全纯映射σ:C*→PC^2,使得(1)σ(C*)=C(2)σ^(-1)(S)是有限点集(3)σ:C*\σ^(-1)(S)→C\S是一对一的映射则称(C*,σ)为C的正则化。

一种新的TTI介质多参数联合层析反演方法

一种新的TTI介质多参数联合层析反演方法

一种新的TTI介质多参数联合层析反演方法郭恺;杨林【摘要】倾斜横向各向同性(TTI)介质具有构造复杂、速度和各向异性参数的数量级差距大、对旅行时同时产生影响等特性,因而层析反演时不稳定且存在多解性.常规的单参数顺序反演方法计算量大,存在误差累积现象,难以建立准确的速度和各向异性参数模型.在常规TTI介质层析反演的基础上,采用理论公式将数量级相差巨大的Thomsen参数转换为数量级一致的速度参数,实现了等效参数同时反演.提出了先采用稳定性较高的顺序反演获得一定精度的参数模型,再使用精度较高的同时反演进一步提高模型精度的联合反演策略.SEG标准TTI模型数据试算和实际资料处理结果表明,联合反演策略精度高、稳定性强,具有较强的有效性和实用性,适用于复杂构造地区大规模实际资料的生产应用.【期刊名称】《石油物探》【年(卷),期】2019(058)003【总页数】8页(P412-418,426)【关键词】TTI介质;层析反演;Thomsen参数;等效参数;顺序反演;同时反演;联合反演;参数耦合【作者】郭恺;杨林【作者单位】中国石油化工股份有限公司石油物探技术研究院,江苏南京 211103;中国石油化工股份有限公司西北油田分公司研究院,新疆乌鲁木齐 830011【正文语种】中文【中图分类】P631传统地震学和地震勘探主要以地球介质具有完全弹性和各向同性的物理假设为基础,由于早期的地震勘探方位较窄、成像方法简单、硬件设施相对落后,因而地震数据体现不出各向异性特性,采用各向同性处理技术能够取得较好的效果。

近几年来,为了获得高品质的地震数据,宽方位地震采集技术得到了广泛应用,地下介质的各向异性问题日益突出;另外,一些高精度地震成像新技术,例如逆时偏移(RTM)成像技术,克服了偏移孔径和偏移倾角的限制,引入了更多的大偏移距数据,因而必须考虑各向异性因素的影响;同时,计算机的发展使各向异性复杂介质高精度处理成为可能。

因此,为了更精确地刻画地质构造,开展地震各向异性处理方法研究是高精度地震成像技术的必然发展趋势。

理论地球物理学的地震层析成像方法

理论地球物理学的地震层析成像方法

理论地球物理学的地震层析成像方法引言地震层析成像是一种利用地震数据推断地下结构的方法,它在地球物理学研究中具有重要的理论和实际意义。

理论地球物理学的地震层析成像方法是基于地震波传播理论和信号处理原理,通过对地震数据进行处理和解释,得到地球内部结构的信息。

本文将介绍理论地球物理学的地震层析成像方法的基本原理、算法和应用。

地震波传播理论地震波是地表上发生的地震源产生的机械波动力。

根据波动方向的不同,地震波可分为纵波(P波)和横波(S波)。

P波是一种有压缩和扩张性的波动,其传播速度较快;S波是一种只能沿垂直于波动方向传播且传播速度较慢的波动。

地震波在地下的传播受到地球结构的影响,由此可以推断地球内部的物理性质和结构。

地震层析成像的基本原理地震层析成像方法基于地震波的传播特性,通过对地震波数据的采集和处理,推断出地下结构的信息。

其基本原理是利用地震波的反射、透射、散射等现象,将地震数据的波形分析和解释,定量地反映地下介质的速度、密度和衰减等特性。

地震层析成像算法地震层析成像算法是将地震数据通过一系列的数学和物理方法进行处理和分析,从而得到地下结构的信息。

常用的地震层析成像算法包括正演算法、反演算法、匹配滤波算法等。

正演算法正演算法是一种将地下结构和初始条件作为输入,通过对地震波方程进行求解,得到地震波的传播情况的方法。

常用的正演算法有有限差分法、波动方程正演法等。

反演算法反演算法是将地震数据作为输入,通过对地震波反问题的求解,推断出地下结构的方法。

常用的反演算法有共轭梯度法、正则化反演法、全波形反演等。

匹配滤波算法匹配滤波算法是一种基于地震数据的频率和波形特征进行分析和处理的方法。

它通过与地下结构的响应进行匹配,提取出地下介质的特征信息。

地震层析成像的应用地震层析成像方法在地球物理学的研究和实践中具有广泛的应用。

以下是地震层析成像在不同领域的应用示例。

石油勘探地震层析成像方法在石油勘探中得到广泛应用。

通过分析地震数据,确定石油或天然气藏的位置、形状和分布,指导油气勘探与开发。

Abel变换数值反演的离散正则化方法

Abel变换数值反演的离散正则化方法

Abel变换数值反演的离散正则化方法
肖庭延;宋金来
【期刊名称】《计算物理》
【年(卷),期】2000(17)6
【摘要】基于Tikhonov的正则化思想 ,将Abel变换的理论反演公式与对数值求导的离散正则化处理以及带权的Gauss型积分相结合 ,给出了Abel变换数值反演的一个新算法 ,并进行了理论分析与数值实验。

结果表明 :该算法具有精度高。

【总页数】9页(P602-610)
【关键词】Abel变换;数值反演;离散正则化;Gauss积分公式
【作者】肖庭延;宋金来
【作者单位】河北工业大学文理学院
【正文语种】中文
【中图分类】O241.8
【相关文献】
1.用正则化方法求一类Abel型积分方程数值解 [J], 罗兴钧
2.Abel变换数值反演的积分算子方法 [J], 吕小红;吴传生
3.Abel变换的数值反演 [J], 赵振宇;贺国强;刘伟
4.数值求解Abel积分方程的正则化方法 [J], 蒋慧琴;许雅静
5.磨光法在Abel变换数值反演中的应用 [J], 刘增锐;郭献洲;张相梅
因版权原因,仅展示原文概要,查看原文内容请购买。

正则化

正则化

正则化(regularization)在线性代数理论中,不适定问题通常是由一组线性代数方程定义的,而且这组方程组通常来源于有着很大的条件数的不适定反问题。

大条件数意味着舍入误差或其它误差会严重地影响问题的结果。

反问题有两种形式。

最普遍的形式是已知系统和输出求输入,另一种系统未知的情况通常也被视为反问题。

许多反问题很难被解决,但是其他反问题却很容易得到答案。

显然,易于解决的问题不会比很难解决的问题更能引起人们的兴趣,我们直接解决它们就可以了。

那些很难被解决的问题则被称为不适定的。

一个不适定问题通常是病态的,并且不论是简单地还是复杂地改变问题本身的形式都不会显著地改善病态问题。

另一方面,病态问题不一定是不适定的,因为通过改变问题的形式往往可以改善病态问题。

在严格的数学意义上,我们通常不可能对不适定问题进行求解并得到准确解答。

然而,通过使用我们的先验知识,我们通常有希望能够得到一个接近准确解答的答案。

求解不适定问题的普遍方法是:用一族与原不适定问题相“邻近”的适定问题的解去逼近原问题的解,这种方法称为正则化方法。

如何建立有效的正则化方法是反问题领域中不适定问题研究的重要内容。

通常的正则化方法有基于变分原理的Tikhonov 正则化、各种迭代方法以及其它的一些改进方法,这些方法都是求解不适定问题的有效方法,在各类反问题的研究中被广泛采用,并得到深入研究。

正则化:Normalization,代数几何中的一个概念。

通俗来说,就是给平面不可约代数曲线以某种形式的全纯参数表示。

即对于PC^2中的不可约代数曲线C,寻找一个紧Riemann面C*和一个全纯映射σ:C*→PC^2,使得σ(C*)=C严格的定义如下:设C是不可约平面代数曲线,S是C的奇点的集合。

如果存在紧Riemann面C*及全纯映射σ:C*→PC^2,使得(1) σ(C*)=C (2) σ^(-1)(S)是有限点集(3) σ:C*\σ^(-1)(S)→C\S是一对一的映射则称(C*,σ)为C的正则化。

初至旅行时层析反演近地表模型精度分析

初至旅行时层析反演近地表模型精度分析

初至旅行时层析反演近地表模型精度分析李宇;杨德义;邓辉;严锐;陈凤英【摘要】为了提高初至旅行时层析反演近地表速度模型的精度,改善静校正效果,本文研究了激发井深以及约束速度对反演精度的影响.通过建立典型低降速层近地表模型,首先对比分析零井深、井深穿过低速层和井深穿过降速层三种情况下,射线穿过网格的特性及激发井深在低速层以下时,最小约束速度对低速层反演效果的影响;其次利用井口时间确定最小约束速度,并提出了辅助点激发确定最小约束速度的新方法;最后应用于合成与实际地震数据,提高了反演精度.【期刊名称】《石油地球物理勘探》【年(卷),期】2011(046)002【总页数】5页(P221-225)【关键词】初至旅行时;层析反演;近地表模型;激发井深;约束速度【作者】李宇;杨德义;邓辉;严锐;陈凤英【作者单位】中国矿业大学(北京)煤炭资源与安全开采国家重点实验室,北京,100083;太原理工大学矿业工程学院,山西,太原,030024;有色金属华东地质勘查局资源调查与评价研究院,江苏,南京,210007;中国矿业大学(北京)煤炭资源与安全开采国家重点实验室,北京,100083;中国矿业大学(北京)煤炭资源与安全开采国家重点实验室,北京,100083【正文语种】中文【中图分类】P631在复杂近地表地区,应用初至旅行时层析反演近地表模型能较好地解决静校正问题,因此该方法在实际中得到广泛的应用。

随着地震勘探精度要求的提高,需要进一步提高近地表模型的精度来改善静校正效果。

同时在复杂近地表地区,认为射线在近地表垂直传播会带来较大误差,基于地表一致性假设的静校正方法已不能满足精度要求,需要运用波场延拓和动态静校正,如波动方程基准面校正和射线追踪静校正[1],这两类方法是解决复杂近地表地区静校正问题最为有效的方法,而采用这类方法需要建立精确的近地表模型,这对初至旅行时层析反演模型的精度提出了更高的要求。

为了提高初至旅行时层析反演近地表模型的精度,人们对初始速度模型和反演算法等进行了大量研究[2~6],却忽略了陆地煤田、油田激发井深参数对反演精度的影响。

地球物理反演中的正则化技术及应用

地球物理反演中的正则化技术及应用

地球物理反演中的正则化技术及应用地球物理反演是研究地球内部结构和性质的重要手段。

正则化技术是地球物理反演过程中常用的数学方法,旨在解决反问题的不适定性和非唯一性,提高反演结果的稳定性和可靠性。

本文将介绍正则化技术的基本原理和常用方法,并探讨其在地球物理反演中的应用。

地球物理反演是根据观测数据推断地下地球结构和性质的过程。

然而,由于地球介质的复杂性以及观测数据的不完备性和噪声污染等因素的影响,地球物理反演问题往往是不适定的反问题。

也就是说,同一个目标模型可以对应多个不同的观测数据解释,使得反演结果存在非唯一性。

此外,反演过程中还可能存在数值不稳定性,即小的扰动可能导致大的误差放大。

为了克服不适定性和非唯一性问题,正则化技术在地球物理反演中得到广泛应用。

正则化通过向反演问题中引入附加信息,对反演结果进行约束,从而提高反演结果的稳定性和可靠性。

正则化的基本思想是在目标函数中同时考虑拟合数据的残差项和模型的先验信息项,通过调节两者之间的平衡,得到最优的反演结果。

在正则化技术中,最常用的方法是Tikhonov正则化。

Tikhonov正则化通过在目标函数中引入二阶范数(L2范数)惩罚项,对模型进行光滑约束,使得反演结果具有空间平滑性。

这种正则化方法在稳定性和保真性之间取得了良好的平衡,常用于地震波速度、重力场和磁场等反演问题。

除了Tikhonov正则化,还有一些其他常用的正则化技术,如L1范数正则化、TV正则化和脉冲响应正则化等。

L1范数正则化通过在目标函数中引入一阶范数惩罚项,将反演结果稀疏化,适用于具有稀疏结构的反演问题。

TV正则化是对图像进行处理的一种方法,通过将图像的梯度惩罚项添加到目标函数中,实现对反演结果的边缘保持和去噪。

脉冲响应正则化是将目标函数中的滤波器参数视作反演模型的参数,通过滤波器设计对反演结果进行约束。

这些正则化技术在地球物理反演中具有各自的优势和适用范围,根据具体的反演问题选择合适的正则化方法进行应用。

地震-地质导向技术在盆缘复杂构造带的应用——以南川区为例

地震-地质导向技术在盆缘复杂构造带的应用——以南川区为例

1751 南川区块概况南川区块位于渝东南盆缘复杂构造带,面积约为1604km 2 ,构造改造作用相对较强,整体北东向展布,呈“四隆四凹”构造格局,南部志留系页岩出露,页岩分布面积525km 2,资源量5500×108m 3。

工区内页岩气主要开发构造单元为平桥、东胜、阳春沟构造带。

开发的目的层段为上奥陶统五峰组-下志留统龙马溪组地层,主力气层段(①-⑤小层)厚33~33.7m。

2 地震-地质导向技术地震数据作为覆盖了整个产建区的唯一三维立体数据,蕴含丰富的地质信息,具有钻井、测井、录井这些点、线数据无法比拟的优势[1]。

在钻头钻进过程中,应用少量高精度的点、线信息迭代校正地震解释成果,及时更新地质模型,地震、地质、测井、录井多专业协同工作,能够极大程度的排除各类地质异常对钻井工程的影响,提高靶窗钻遇率[2]。

2.1 速度模型构建为实现工区复杂高陡构造准确成像,减小井震误差,处理-解释一体化合作,交互作业,形成了层控、断控、井控的“三控”特色速度建模技术。

(1)构造层位约束的局部速度模型构建技地震-地质导向技术在盆缘复杂构造带的应用——以南川区为例倪锋中国石化华东油气分公司勘探开发研究院 江苏 南京 210000摘要:南川页岩气田位于川东南盆缘构造带,目标层系五峰组—龙马溪组经过多期构造运动改造,地层变形强烈,断裂、褶皱发育,地层产状变化剧烈,水平井的地质导向工作难度大。

为保证水平井顺利入窗,提高靶窗钻遇率,本文以地震-地质一体化导向技术为研究目标。

建立了层控、断控、井控的“三控”特色速度建模技术,得到较为可靠的深度域数据体。

以深度域地震数据为基础,建立工区构造模型、地质模型,立体把控目的层空间展布形态,优化了靶点深度和水平段地层产状的预测方法。

综合分析认为,地震-地质导向技术在水平井精准着窗、靶窗平稳穿行等导向关键环节均呈现了良好的指示意义。

是复杂构造带油气效益勘探开发行之有效的技术方法,值得在同类型区块推广应用。

井间地震资料全变差正则化波形反演

井间地震资料全变差正则化波形反演

井间地震资料全变差正则化波形反演王守东【期刊名称】《石油物探》【年(卷),期】2011(050)003【摘要】针对井间地震资料走时层析成像分辨率低、稳定性差的问题,研究了一种井间地震资料成像新方法--全变差正则化波形反演方法(TVWI).该方法基于Gauss-Newton波动方程反演,将井间地震资料的速度求取问题归结为波动方程反演问题,通过极小化合成地震记录与观测地震记录的差实现波动方程反演;利用局部线性化逐步迭代的方法实现非线性泛函优化问题的求解;将波动方程反演的目标泛函加入全变差约束项,提高反演过程的稳定性和反演结果的分辨率.设计了一个复杂的速度模型验证了TVWI方法的有效性,结果表明,在给定不同速度初值的情况下,反演结果都具有较高的精度.选取与真实速度相差较大的初始速度模型讨论了TVWI方法对初值选取的依赖性:反演结果对初值选取的依赖性不强,但初值选取不好,会增加迭代次数和计算成本.最后进行了抗噪能力试验,结果表明,TVWI方法具有较强的抗噪能力.【总页数】7页(P234-240)【作者】王守东【作者单位】中国石油大学(北京)地球物理与信息工程学院CNPC物探重点实验室,北京102249【正文语种】中文【中图分类】P631.4【相关文献】1.不连续介质反演的原对偶牛顿法和全变差正则化 [J], 冯立新;李媛;张磊2.一种反演抛物型方程未知系数的全变差正则化方法 [J], 李照兴;张雷;胡志东3.井间地震资料的声波波动方程走时和波形反演 [J], Zhou,CX;张恒山4.基于改进全变差正则化的GPR多尺度全波形双参数同步反演 [J], 王珣;冯德山;王向宇5.不连续介质反演的原对偶牛顿法和全变差正则化 [J], 李娜;冯立新因版权原因,仅展示原文概要,查看原文内容请购买。

环状流ERT图像重建算法研究

环状流ERT图像重建算法研究

环状流ERT图像重建算法研究环状流ERT(Electrical Resistivity Tomography,电阻率层析成像)是一种非侵入式地下勘探方法,它利用地下介质的电阻率差异反演出地下结构的分布情况。

环状流ERT算法是一种改进的ERT方法,它利用环形电极的布置和交替通电方式提高了测量精度和空间分辨率。

本文将重点介绍环状流ERT图像重建算法的研究进展及其存在的问题和挑战。

一、环状流ERT图像重建算法研究现状目前,环状流ERT图像重建算法主要可以分为以下几类:(1)直接使用常规ERT算法进行反演;(2)利用环状流ERT的特殊性质设计新的数据处理方法;(3)利用约束反演等逆问题求解方法进行反演;(4)结合其他地球物理探测方法进行反演。

其中,第一类方法基本是将环状流ERT看作一种特殊的ERT方式,直接采用常规的ERT算法求解反演问题。

这种方法的优点是简单易于实现,但由于未能充分利用环状流ERT的特殊性质,精度和空间分辨率都较低。

第二类方法则是根据环状流ERT的特殊性质,设计了一些特殊的数据处理方法。

例如Yang等人(2019)提出了一种基于自回归移动平均(ARMA)模型的数据拟合方法,在保证数据平滑性的同时,可以自动修正电极不对称性和非理想等效电路的影响,提高了反演精度。

Shi等人(2020)则提出了一种基于多方向分段梯度的反演方法,能够提高反演效率和空间分辨率。

第三类方法则是采用逆问题求解方法,利用各种约束条件(如先验信息、正则化等)对反演问题进行求解。

刘国阳等人(2018)采用了一种基于反演先验信息的逆问题求解方法,在提高反演精度的同时保证了图像的平滑性和连通性。

陈彬等人(2019)则利用了约束最小二乘(CLS)方法进行反演,在充分考虑多种数据约束条件的情况下,提高了反演精度和稳定性。

第四类方法则是将环状流ERT与其他地球物理探测方法相结合,如采用声波、磁法或重力等方法共同反演地下结构。

徐琛等人(2020)采用了环状流ERT与声波波动方程联合反演的方法,在提高图像分辨率和减小模型误差的同时,得到了更为准确的地下结构信息。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

李辉,王华忠,张兵.层析反演中的正则化方法研究[J].石油物探,2015,54(5):569 -581Li Hui,Wang Huazhong,Zhang Bing.The study of regularization in tomography[J].Geophysical Prospecting for Petroleum,2015,54(5):569 -581收稿日期:2014-11-24;改回日期:2015-02-26。

作者简介:李辉(1985—),男,博士,现从事射线类偏移与反演的研究工作。

基金项目:国家自然科学基金(41374117)、国家重点基础研究发展计划(973计划)项目(2011CB201002)、国家科技重大专项项目(2011ZX05003-003,2011ZX05005-005-008HZ,2011ZX05006-002)和中国石化地球物理重点实验室开放基金项目(33550006-14-FW2099-0026)共同资助。

层析反演中的正则化方法研究李 辉1,2,王华忠1,张 兵1,3(1.同济大学海洋与地球科学学院波现象与反演成像研究组,上海200092;2.青凤致远应用地球物理研究所,上海200093;3.中国石油化工股份有限公司石油物探技术研究院,江苏南京211103)摘要:正则化可显著降低层析反演解的非唯一性,提高层析反演结果的质量。

主要研究了模型参数正则化和数据正则化。

地下介质参数之间的关联性如何加入模型正则化是讨论的问题之一;观测数据之间的关联性加入数据正则化的方法则是另一个主要议题。

此外,讨论了Tikhonov正则化和预条件两种模型正则化实现策略,指出前者理论比较直观,后者计算效率更高,并证明了两者在理论上的等价性。

模型正则化通过构造各向异性光滑算子加入地质构造特征,数据正则化则通过在层析矩阵中加入预先构造的数据预条件矩阵来实现。

通过层析偏移速度分析给出了模型正则化和数据正则化的具体实现策略。

理论分析和层析偏移速度分析的数值实验说明本文的模型正则化和数据正则化可显著提高层析反演的质量。

关键词:层析偏移速度分析;模型正则化;数据正则化;预条件;地质构造约束中图分类号:P631文献标识码:A文章编号:1000-1441(2015)05-0569-13 DOI:10.3969/j.issn.1000-1441.2015.05.010The study of regularization in tomographyLi Hui 1,Wang Huazhong1,Zhang Bing1,2(1.Wave Phenomena and Inversion Imaging Group(WPI),Tongji University,Shanghai 200092,China;2.Qingfeng-zhiyuan Applied Geophysics Institute,Shanghai 200093,China;3.Sinopec Geophysical Research Institute,Nanjing211103,China)Abstract:Regularization in tomography is able to weaken the non-uniqueness of tomography to improve the inversion result-The discussion of regularization in this paper includes model-regularization and data-regularizationModel parameters are not i-solated,how to add the relationship of these parameters into tomography is one of the missions hereSimilarly,considering da-tum relationship in tomography is another problemThe so-called“straightforward regularization”and the“precondition regu-larization”are focused,and we achieve that the former is intuitionistic and the latter is more efficiencyAlso,we point out thatthe above two algorithms are equivalent to each other,and this will be shown in this paperThe geological structure character-istics of the medium can be integrated into the tomography using the model-regularization with anisotropic smooth matrix.The data-regularization is realized with another smooth operator which will be integrated into the tomographic matrix.Themodel-regularization and data-regularization are tested with tomographic migration velocity analysis(MVA)algorithm.Theresults of theory and numerical experiments with tomographic MVA show that the proposed model-regularization and the da-ta-regularization are both able to improve the quality of tomography obviously.Keywords:tomographic MVA,model-regularization,data-regularization,precondition,geological structure constraint 随着勘探地震技术的发展以及石油工业需求的提高,叠前深度偏移逐渐成为工业应用中偏移技965第54卷第5期2015年9月石 油 物 探GEOPHYSICAL PROSPECTING FOR PETROLEUMVol.54,No.5Sep.,2015术的主流。

在观测数据品质有保障的前提下,叠前深度偏移在实际应用中成功与否最重要的前提是速度模型的准确程度。

深度域速度建模的方法包括分析类方法(如扫描速度分析、剩余曲率分析)、层析方法、波动方程类反演(如全波形反演)方法等。

其中,分析类方法存在较强的理论假设,导致估计的速度精度及分辨率非常低;波动方程类反演方法的效率非常低,且严重依赖于初始速度模型;基于射线理论的层析方法是目前工业界应用最广泛的深度域速度估计方法[1-3],但存在明显的缺点,如:只能反演光滑的背景速度、建立的反问题病态性较明显等,所幸的是,正则化可显著解决射线层析中的上述问题[4-5]。

尽管正则化的本质思想十分明确,但具体实现方法及种类繁多,本文针对模型参数之间的关联性和数据之间的关联性讨论层析反演中加入正则化的策略。

正则化是地球物理反问题中非常重要的环节,可以显著降低层析反演解的非唯一性,提高层析反演结果的质量。

地震层析反演中的模型参数是地下介质参数(通常是地震波速度),数据是观测到的波场信息(通常是旅行时)。

地下不同空间区域的介质参数之间有一定关联性;同样,不同观测点的数据也相互关联。

此外,测量的地震波数据良莠不齐,可靠性存在明显差异。

如何把这些信息加入层析过程以提高反演质量是一个值得研究的课题。

在地震层析反演中,利用地下构造信息约束模型参数的空间分布特征,把模型参数在空间中的相关特征通过构造信息提取出来,结合Tikhonov模型正则化或预条件思想把此信息加入层析反演中,可显著改善估计的模型参数,这里称之为模型正则化。

Zhou等[6]、Zhou[7]和Zdraveva等[8]把光滑算子加入灵敏度核函数中对反演的模型直接进行正则化约束;Clapp[9]和Clapp等[10]利用预条件思想实现对反演模型的预条件约束。

事实上,Tik-honov模型正则化和预条件模型正则化等价,虽然两者有不同的表达形式及理论基础。

类似于模型参数,利用协方差矩阵可提取出观测数据之间的相关性,将数据相关性作为数据先验信息加入反演过程是贝叶斯反演思想[4]的一个关键环节。

此外,测量的数据有优劣之分,反演中如何通过正则化突出质量较高的数据也是本文的研究内容。

层析中加入数据关联性和突出高质量的数据称为数据正则化。

1 模型正则化理论及实现策略层析反演中针对模型施加正则化是把我们对模型参数的认识加入层析反演中。

地下介质的模型参数化后,不同参数之间存在一定的关联性,因为所有参数组成的地下介质遵循一定的地质规律。

这里研究的模型正则化试图在层析过程中加入地质构造约束,我们通过在层析矩阵中加入光滑矩阵来实现。

光滑矩阵可约束模型中不同方向的突变成分,令垂直于反射界面与平行于反射界面的方向光滑强度不同来实现对模型的各向异性光滑约束。

模型正则化的具体实现策略包括Tikhonov模型正则化和预条件模型正则化两种。

无论是Tikhonov正则化还是预条件模型正则化,光滑矩阵的构建是重要的一环。

1.1 Tikhonov模型正则化贝叶斯框架下层析反演等价于优化如下的目标泛函[4]:J(m)=12{[A(m)-d]TC-1D[A(m)-d]+ε(m-m0)TC-1M(m-m0)}(1)式中:m和d是模型参数与观测数据;m0是参考模型;A(m)是模拟数据;C-1D和C-1M分别是数据协方差矩阵和模型协方差矩阵的逆矩阵;ε是阻尼因子。

层析反演通过线性化的迭代实现对上述目标函数的优化,一次迭代的线性方程为:(ATC-1DA+εC-1M)Δm=ATC-1DΔd(2)式中:矩阵A是A(m)的线性化算子;Δm=m-m0是当前迭代中模型的更新量;m0是参考模型;Δd=Am0-d是当前迭代中的数据残差。

目标泛函((1)式)和方程(2)中的协方差矩阵及阻尼因子实现了对层析反演施加正则化。

不考虑正则化,仅用数据残差进行层析反演,即协方差矩阵为单位矩阵且ε=0时,层析方程(2)退化成一般的层析方程,即:AΔm=Δd(3) 令方程(2)中的数据协方差矩阵为单位矩阵,只考虑模型正则化时的层析方程为:ATAΔm+εC-1MΔm=ATΔd(4) 模型协方差矩阵的逆矩阵不能直接构造,一般利用其它矩阵代替,从而实现对模型的正则化约束。

相关文档
最新文档