共查询到20条相似文献，搜索用时 671 毫秒
1.
针对利用地震道进行相对波阻抗反演中遇到的横向连续性难以保持、初始子波容错度差以及随机噪声干扰影响反演结果等问题,提出了一种基于矩阵Toeplitz稀疏分解的相对波阻抗反演方法.该方法将地震数据剖面的Toeplitz稀疏分解问题分解为两个子反演问题,其一以Toeplitz子波矩阵元素为待反演的参数,用Fused Lasso方法求解,可保证子波具有紧支集且是光滑的;其二以稀疏反射系数矩阵元素为待反演参数,用基于回溯的快速萎缩阈值迭代算法求解,大大降低了目标函数中参数选择的难度.通过交替迭代求解上述两个子反演问题可将地震数据剖面因式分解为一个Toeplitz子波矩阵和一个稀疏反射系数矩阵;然后由反射系数矩阵递推反演可以得到高分辨率的相对波阻抗剖面;利用测井资料加入低频分量后,也可得到高分辨率的绝对波阻抗剖面.Marmousi2模型生成的合成记录算例和实际地震资料算例均表明:本文方法可以从带限地震数据中有效地反演相对波阻抗,反演结果分辨率高并且能够很好地保持地震数据的横向连续性;即使在初始估计子波存在误差和地震数据被随机噪声污染的情况下也能取得较好的效果. 相似文献
2.
时间域航空电磁数据的反演 总被引：3，自引：0，他引：3
本文用广义逆矩阵理论讨论了时间域航空电磁数据的反演方法.以水平二层大地和球体模型为例进行了计算,结果表明这种方法对时间域电磁数据的反演是行之有效的.对于理论数据,一般仅需迭代数次即可稳定地收敛到真值.通过奇异值分析,可以确定一个模型中的相对重要参数和无关紧要参数;对于求解某一特定参数,可确定哪些数据起决定性的作用.由奇异值分解而得出的信息密度矩阵、模型分辨矩阵,分别给出了模型响应拟合数据和模型参数分辨率的度量.这些信息对于时间域航空电磁系统的设计、野外测量和资料解释都具有重要的参考价值.最后给出了实例. 相似文献
3.
基于知识的分阶段不确定性多角度遥感反演 总被引：12，自引：4，他引：8
针对多角度遥感的反演提出了一种基于知识的分阶段不确定性反演方法 .反演以Tarantola的信息理论为基础将人类的先验知识、BRDF模型以及实验测量数据结合在一起 ,在反演决策上使用了多阶段的目标决策方案 .最后以SAIL模型为例 ,用室内实验测量数据对该反演方法进行了反演试验 ,反演结果表明这种反演方法较传统的确定性反演方法更为合理、可靠 . 相似文献
4.
地球表面的各向异性特性可以用地表二向反射函数（BRDF）恰当地描述.BRDF的核心是利用线性核驱动模型,数学上表述为各向同性核、体散射核和几何光学核的线性组合.随着多角度遥感领域的发展,BRDF模型越来越被看作是可以反演重要的有关地表生物的或气候的参数,比如说叶面积指数和地表反照率.一个线性逼近的核驱动BRDF模型通常可以写成下述形式（Roujean等,1992）：fiso＋kvol（ti,tv,φ）fvol＋kgeo（ti,tv,φ）fgeo=r（ti,tv,φ）,其中r表示地表的二向反射;kvol和kgeo为通常所说的核函数,即为已知的入射和观测几何特性的函数,分别描述了体散射和几何散射（包括折射和反射）;ti是太阳方向天顶角,tv是观测方向天顶角;φ表示太阳-观测方向的相对方位角;fiso,fvol和fgeo为未知的待反演参数,可以用来拟合观测.计算过程的稳定性是由核矩阵的代数算子特征谱和观测噪音/误差来刻画的.因此为了计算地表反照率,成功反演模型参数是至关重要的环节.我们首先考虑了为计算BRDF模型反演的光滑解方法.业已知道,这是一个不适定的反问题.不适定性是由线性核驱动BRDF模型的欠定性表征的,比如说观测严重不足或观测方向范围有限,或者是观测数据高度线性相关以及噪音的污染等.例如,一次单角度观测可以导致一个欠定的系统（核算子的零空间含有非零向量）或者系统无解（系数矩阵的秩不等于增广矩阵的秩）.因此,光滑性或正则化技巧应当加以利用来压制不适定性.Li等（2001）应用先验知识把原始模型转换为一个超定的模型并求得最小二乘解.Pokrovsky等（2002）应用QR分解反演BRDF模型.Wang等（2007）考虑到了反演的正则化策略并提出了不适定地表参数反演的一个完整的正则化理论.在文中,强调从不同的空间添加先验信息于反演模型中.首先从数学物理的观点,第一次提出了? 相似文献
5.
AVO反演可以获得地层岩性和流体信息,而叠前反演问题都是高维的和非适定的,因此获得可靠稳定的解对叠前反演至关重要. 本文给出了一种基于贝叶斯理论的纵波和转换波联合反演密度比和模量比的方法. 鉴于剪切模量比、体积模量比可以更好地指示油气,基于岩石物理中速度比与模量比之间的关系,将此关系式代入Zoeppritz方程的近似形式Aki-Richards公式中,得到与模量比有关的反射系数近似公式. 联合纵波和转换波,利用最小二乘准则构建目标函数,最终反演出密度比、剪切模量比、体积模量比三个参数. 在反演过程中引入贝叶斯理论,假定先验信息服从高斯分布,待求参数服从改进的Cauchy分布,并去除待求参数之间的相关性. 利用模型数据和实际数据对本文方法进行测试,并与常规的单独利用纵波数据来反演方法进行比较,结果表明联合反演稳定性更好、精度更高、抗噪音能力更强,验证了本文方法的可行性和有效性. 相似文献
6.
采用弹性波全波形反演方法精确重建深部金属矿多参数模型,建模过程采用基于地震照明的反演策略.首先给出基于照明理论的观测系统可视性定义,利用可视性分析构建新的目标函数,对反演目标可视性较高的炮检对接收到的地震记录在波场匹配时占有更高的权重,确保了参与反演计算中的地震数据的有效性;其次将给定观测系统对地下介质的弹性波场照明强度作为优化因子,根据地震波在波阻抗界面处的能量分配特点,自适应补偿波场能量分布和优化速度梯度,以提高弹性波全波形反演过程的稳定性和反演结果的精度.理论模型和金属矿模型反演试验结果表明,基于可视性分析和能量补偿的反演策略可以使弹性波全波形反演更快地收敛到目标函数的全局极小值,获得适用于金属矿高分辨率地震偏移成像的多参数模型. 相似文献
7.
在实际的AVO反演问题中，叠前数据体中的噪声或其他因素严重影响了AVO反演问题的适定性，而采用先验地质信息作为AVO反演问题的约束条件是解决AVO反演问题不适定的一种可行方法. 文中的似然函数采用了[WTBX]ι[WTBX]p范数的解，并用Cauchy分布表示先验模型参数的分布. 以此为基础，在反演中建立了测井数据的参数协方差矩阵对反演过程进行约束，并采用了共轭梯度算法实现多参数非线性的反演过程. 同时，为了提高反演精度，避免动校正拉伸及依赖于炮检距的调谐效应对参数估计的影响，反演采用动校前地震数据进行参数估计. 从应用效果分析来看，即使叠前道集的信噪比不高，反演的结果也能较好地与实际情况相匹配，为识别储层流体性质提供了新的手段. 相似文献
8.
常规AVO三参数反演是通过Zoeppritz方程的近似公式来建立AVO正演模拟的过程,然而在P波入射角过临界角和弹性参数在纵向上变化剧烈的情况下,Zoeppritz方程近似公式精度有限.针对这种情况,可以使用精确的Zoeppritz方程来构建反演目标函数,由于精确Zoeppritz方程中P波反射系数和弹性参数之间是一种复杂的非线性关系,通常解决途径是利用非线性的优化算法来进行数值计算,但是非线性优化算法的缺点是计算量过大;另外一种途径是利用广义线性反演的方法,通过泰勒一阶展开式将P波反射振幅展开后,用线性关系近似表达非线性关系,经过几次迭代后,在理论上可以达到很高的精度,但是广义线性反演算法的核心部分--Jacobian矩阵由于矩阵条件数过大,往往会造成反演算法的不稳定,其应用范围得到了限制.贝叶斯反演方法是通过引入模型参数的先验分布结合噪声的似然函数,生成模型参数的后验分布,通过求取模型参数的最大后验概率分布来得到模型参数的反演解,由于引入模型参数的先验分布信息,可以有效的降低反演的不适定问题.本文将两种反演算法的思想相结合,利用广义线性反演算法的思想,构建AVO正演模拟的过程来提高大角度地震数据反演的精度,同时结合贝叶斯理论,通过引入模型参数的先验分布信息构建反演目标函数的正则化项,可以有效降低由于Jacob矩阵条件数过大带来的反演不适定问题,该算法假设模型参数服从三变量柯西分布. 相似文献
9.
10.
三维密度反演是地球物理领域的研究热点,而在大数据及人工智能发展的时代背景下如何快速高效地实现重力数据反演显得更为重要.传统反演方法通常需要存储大型系数矩阵,内存占用大,耗费时间长,同时为约束反演结果而加入的正则化约束项参数难以确定;深度学习可以不依赖先验信息,也不需要计算及存储系数矩阵,使得计算效率大大提高.基于此,本文提出了一种基于UNet++网络的重力异常反演方法.首先将UNet++网络中部分参数进行更改：选择在输入数据绝对值较大时梯度更稳定的LeakyReLU作为激活函数;加入了Batch Normalization层,增强了网络的收敛速度及稳定性.然后为了提高网络的全局最优化能力,引入了基于余弦退火的学习率更新策略,使用梯度的一阶以及二阶矩估计的Adam最优化算法,利用数据集与标签集进行网络训练,实现了重力异常的三维密度反演.通过实验验证了UNet++、LeakyReLU更快速稳定的收敛能力,而余弦退火学习率更新策略具有更强的全局寻优能力.含噪模型实验及实际数据反演结果进一步证明该方法的正确性和有效性,及其良好的泛化能力与抗噪能力. 相似文献
11.
叠前反演的目的是基于弹性波理论从地震数据中获得地层参数的可靠估计,进而用于描述地层的流体和岩性特征.然而叠前反演问题都是高维的和非适定的,并且容易受各种噪声和采集过程中不确定因素的影响,因此,为了获得稳定可靠的解必需对反演过程加以合理的约束.本文提出了一种基于非线性二次规划的叠前三参数反演方法.首先基于贝叶斯参数估计理论,假设似然函数服从高斯分布,并使待反演的参数服从于改进的Cauchy分布,从而提高了反演结果的分辨率;其次用协方差矩阵来描述参数间的相关程度,进一步提高了反演结果的稳定性;最后将问题转化为一个非线性二次规划的求解问题,并在多种约束下得到问题的解.仿真实验和实际应用皆已表明,本文提出的反演方法运算速度快捷,既使在信噪比很低的情况下也可获得较好的反演结果,为储层的进一步识别提供更多的物性参数. 相似文献
12.
地球介质的黏度是地球动力学研究中的重要参数. 在考虑震后震源区介质物性变化的情况下,本文提出利用GPS观测的时间序列资料反演地壳、地幔黏度的方法. 计算中,将遗传算法与Marquardt法相结合,利用遗传算法搜索结果作为Marquardt法反演的初始模型,运用微扰法形成Jacobian矩阵. 根据余震分布及地震波速度变化确定震源区的范围及介质的物性变化情况,利用黏弹性有限单元模型（Maxwell体）及台湾集集地震后GPS观测的时间序列数据,反演了台湾地区地壳、地幔的黏性系数. 反演计算及数值实验表明,该方法反演黏度的效率高,稳定性好. 反演结果初步显示台湾地区的地壳、地幔的黏度分别为12×1018Pa·s和36×1019 Pa·s. 相似文献
13.
随着重力和重力梯度测量技术的日趋成熟,基于重力和重力梯度数据的反演技术得到了广泛关注.针对反演多解性严重、计算效率低和内存消耗大等难点问题,本文开展了三维重力和重力梯度数据的联合反演研究,该方法结合重力和重力梯度两种数据,将L0范数正则化项加入到目标函数中,并在数据空间下采用改进的共轭梯度算法求解反演最优化问题.同时,本文摒弃了依赖先验信息的深度加权函数,引入了自适应模型积分灵敏度矩阵,用来克服因重力和重力梯度数据核函数随深度增加而衰减引起的趋肤效应问题.为了提高反演计算效率,本文又推导出基于规则网格化的重力和重力梯度快速正演计算方法.模拟试算表明,改进的共轭梯度法可以降低反演的迭代次数,提高反演的收敛速度;自适应模型积分灵敏度矩阵,可以有效解决趋肤效应,提高反演纵向分辨能力;数据空间和改进的共轭梯度算法结合,可以更好地降低反演求解方程的维度,避免存储灵敏度矩阵,有效地降低反演计算时间和内存消耗量.野外实例表明,该算法可以在普通计算机下快速地获得地下密度分布模型,表现出较强的稳定性和适用性. 相似文献
14.
重力反演是恢复地下密度空间分布的有效工具,而选择合理的密度模型约束方法是提升重力反演分辨率和可靠性的关键.常规约束方法大多是从剖分网格空间中的密度模型出发,通过调整光滑或稀疏约束权重来匹配反演目标,但当地质体类型多样、异常分离不准确及网格剖分方案不合理时,模型约束的合理性与灵活性难以得到有效保证.为此,本文提出了一种基于密度模型稀疏表征的重力反演方法.首先假设待反演的密度模型表征为模型特征矩阵和稀疏分解系数的线性组合,之后重新推导了重力反演目标函数,并给出了分解系数的稀疏求解过程.相比现有重力反演方法,用于构建模型特征矩阵的特征模型可包含不同类型地质体的先验几何信息,分解系数的稀疏性保证了待反演目标来自于最典型的地质模式组合.最后,通过模型试验及实际资料验证了基于密度模型稀疏表征的重力反演方法的有效性. 相似文献
15.
慢度平方三角网格模型VTI介质二维QP波各向异性立体层析反演 总被引：1，自引：1，他引：0
基于慢度平方三角网格剖分模型下的各向异性射线扰动理论,建立了具有垂直对称轴的横向各向同性(VTI)介质中二维QP波立体层析所需的FRECHET导数矩阵,实现了慢度平方三角网格模型的VTI介质二维QP波各向异性立体层析.考虑到各向异性参数是影响地震波运动学特征的次一级因素,制定了先反演慢度再反演各向异性参数的反演策略.由于慢度平方三角网格模型的稀疏性,使得数据空间对模型空间的FRECHET导数矩阵规模被大幅压缩,在降低了计算成本的同时也很好的保证了界面位置、背景慢度、各向异性等各个参数的反演精度.理论数据算例证实了 FRECHET导数的正确性和反演策略的合理性. 相似文献
16.
基于交叉梯度结构约束的不同数据的联合反演可以提高地球物理成像的可靠度,但是由于不同观测数据对地下介质的灵敏度不同以及不同反演系统灵敏度矩阵元素的数值可能存在较大的差异,导致构建和求解联合反演系统存在很大的挑战.针对以上问题,本文提出一种新的基于单独反演模型更新量与交叉梯度结构约束相结合的联合反演策略.该策略利用单独反演系统分别确定出各个模型的更新量,然后利用它们约束交叉梯度系统的反演,得出新的模型更新量.通过这样的方式,有效实现了数据拟合与结构约束的平衡,实现了基于交叉梯度联合反演的目标.新的联合反演策略不需要对原来单独反演程序修改即可实现联合反演,减小了联合反演实现的难度,极大地提高了联合反演的易实现性,而且避免了联合反演矩阵存储及结构过于复杂难以求解的问题.基于新的策略,本文首次实现了基于交叉梯度结构约束的三维地震走时与直流电阻率联合反演.合成模型测试表明,与单一成像相比,联合成像减少了地震走时反演中出现的干扰异常并提高了电阻率反演的分辨率. 相似文献
17.
《地球物理学进展》2015,(5)
在海域天然气水合物勘探过程中,全波形反演可以在缺少井资料情况下提供水合物及游离气区域精细速度结构,与AVO分析相比,全波形反演可以有效降低薄层效应影响.概述了全波形反演的研究进展及其在海域天然气水合物参数结构分析中的应用现状;讨论了天然气水合物全波形反演策略,包括基于层状介质假设的全局优化与局部优化相结合反演策略和基于有限差分、有限元正演模拟的多参数结构直接反演策略等;探讨了天然气水合物全波形反演中一些具有代表性的技术特点,包括地震子波估计方式、海底多分量地震数据的应用、多参数同时反演;最后给出天然气水合物全波形反演发展趋势和可能改进建议.随着三维全波形反演算法的不断改进,海域天然气水合物勘探精度将进一步提高,为水合物商业开发提供依据. 相似文献
18.
有别于传统基于梯度信息的反演方法在正则化约束中用总梯度逼近海塞逆矩阵的技术,本文将正则化约束问题的数据拟合项和模型光滑项分开考虑,只利用数据拟合函数的梯度信息对数据拟合项的海塞矩阵进行逼近,通过求解类高斯牛顿下降方向方程得到不依赖前几次迭代正则化因子的更精确下降方向,在求解当前迭代下降方向的过程中,通过保证右端项中两个向量的二范数在同一数量级的原则,实现了正则化因子的自动更新.对理论模型的试算表明这种自适应正则化反演方案可以在拟牛顿反演框架下基本达到OCCAM的算法稳定性,反演结果对初始模型依赖性较小,同时又无需在一次迭代中多次搜索最佳正则化因子.本文还基于此算法讨论了大地电磁各参数对于反演结果的影响,由于本文的反演结果能得到充分的正则化约束,因而在此框架下讨论阻抗和倾子在反演中的作用相对更为客观. 相似文献
19.
海面浮油膜是海洋石油污染监测与海洋油气资源勘探的共同关注对象,其厚度的遥感定量反演是研究难点之一.随着油膜厚度的变化,入射光在浮油膜中的辐射传输有所不同,是海面浮油膜厚度光学遥感反演的理论依据.本研究将均匀分布的浮油膜视为一个平行平板,分析其可见光范围内的辐射传输过程,探讨油膜反射光的双光束干涉,构建以油膜厚度及其光谱反射率为变量的油膜厚度遥感定量反演模型,并结合菲涅耳公式进行模型参数的物理意义分析,揭示油膜厚度遥感定量反演的关键在于油膜的消光特性.利用油膜光谱响应实验数据,开展理论模型的分析与验证,结果表明：油膜厚度遥感定量反演的理论模型不仅具有明确的光学物理意义,同时也具有较好的模拟精度;浮油膜光谱反射率随厚度变化响应的主要因素是油膜对入射光的分子散射作用,该散射作用与入射光波长具有显著的关系,因此油膜光学遥感散射特性研究将有助于油膜种类识别与厚度的遥感定量反演;模型参数的理论意义与数据分析表明,建立典型油种油膜的模型参数查找表,将有效提高海面浮油膜厚度遥感定量反演的效率. 相似文献
20.
在界面两侧地层的弹性参数弱反差的假设难以成立的情况下,本文提出用平均入射角道集进行PP波与PS波的联合反演.首先,在PP波与PS波AVA(amplitude versus angle,振幅随入射角变化)道集的基础上,分别选择小入射角范围与大入射角范围的AVA道集进行局部加权叠加,以获得由两个角度组成的平均入射角道集,并作为后续反演的输入数据.然后,再通过最小二乘原理建立了PP波与PS波联合反演目标函数,推导了模型修改量的向量公式,建立了平均入射角道集联合反演的流程.模型数据与实际数据的测试结果表明:在信噪比较低、地层弹性参数反差较大、层厚较薄的情况下,该反演方法的精度在很大程度上超过了基于近似反射系数的反演方法,为复杂油气藏勘探提供了新的思路. 相似文献