共查询到19条相似文献，搜索用时 187 毫秒
1.
由于传统的阻尼最小二乘法只适合于模型较少的简单模型,因此当介质的层数较多时,反演就会受到多解性的影响,有时甚至出现不收敛的情况,并且反演十分耗时。为此,使用正则化思想引入模型约束进行反演,且正则化因子通过计算每次迭代的数据目标函数和模型目标函数自适应得到,使反演能够稳定地进行;引入拟牛顿法来更新雅可比矩阵,大大缩短反演所需要的时间,通过典型的3层与多层理论模型的反演试算,证明了拟牛顿法自适应正则化反演算法对初始模型的要求不高,拟合效果好,收敛速度快,适应性强,体现了良好的稳定性和可靠性。 相似文献
2.
《物探与化探》2021,(5)
长导线源半航空瞬变电磁正则化反演正则项通常采用L2范数,其拟合结果较光滑,不能有效刻画层界面信息。针对层状介质陡变模型实现正则项为L1范数的反演算法,采用迭代重加权最小二乘法将原问题转化为L2正则化子问题求解,解决L1范数存在不可导问题;采用OpenMP技术对雅可比矩阵并行计算,提高了反演速度;对自适应正则化因子分段迭代法的调整策略进行分析并改进,改进后的自适应正则化因子调整策略更适合半航空瞬变电磁L1正则反演算法。最后对电阻率进行反演并与Occam反演结果作比较,结果表明L1正则反演充分迭代后能够突出符合真实模型的电性界面,反演电阻率与模型真实值更接近。 相似文献
3.
二维大地电磁正则化共轭梯度法反演算法 总被引：2，自引：0，他引：2
针对大地电磁二维反演中目标函数收敛速度慢而且解的稳定性较差等问题,提出了大地电磁数据的正则化共轭梯度法反演算法(Regularized Conjugate Gradient Algorithm,RCGA)。此算法在构建目标函数时引入正则化的思想,利用共轭梯度法求解最优化问题。在每次迭代过程中根据目标的收敛情况更新正则化因子,有效地解决了迭代时目标函数发散的问题。最后对两个典型模型进行了试算,对其收敛速度及解的稳定性进行了分析,并与传统反演方法的计算结果进行了对比,结果都表明RCGA具有收敛速度快、精度高、结果稳定等优点。 相似文献
4.
5.
层析反演的本质是解病态方程组,且系数矩阵通常为大型稀疏矩阵。正则化可提高解的稳定性,LSQR法可减小存储空间、加快计算效率,因此将二者结合提出了基于光滑约束正则化的层析反演算法。正则化参数在反演中有着至关重要的作用,采用L曲线法选择最佳正则化参数。列举了基于走时层析成像的数值计算实例,由最佳正则化参数下的速度剖面图可清晰地识别异常体特征,证明了基于光滑约束正则化的层析反演算法的有效性。 相似文献
6.
为了减小地球物理反演的多解性,通常采用累加型目标函数,即数据拟合差项加上正则化项。基于累加型目标函数的反演存在如何优选正则化因子的问题,这个过程通常需要做大量的反演计算。本研究系统分析和实现了基于乘积型目标函数(数据拟合差项和正则化项相乘)的电阻率法三维反演,乘积型目标函数的反演不存在优选正则化因子的问题。三维正演采用非结构化有限单元法,三维反演采用有限内存拟牛顿方法。使用理论模型合成数据进行了三维反演,检验了基于乘积型目标函数的电阻率法三维反演的可行性与有效性。反演算例结果表明：基于乘积型目标函数的反演方法能够可靠恢复异常体的电阻率值、形态和位置。 相似文献
7.
波阻抗约束反演的几种方法 总被引：2，自引：0，他引：2
介绍三种理论与实际效果比较好的波阻抗反演方法,即宽带约束反演(BCI)、快速模拟退火算法(FSA)和正则化方法。BCI是一种线性迭代方法,FSA是一种具有全局寻优特点的随机搜索方法,正则化方法是通过对目标函数施加先验约束,来解决反演的病态问题———多解性和不稳定性。最后,使用正则化方法并结合快速模拟退火算法,进行了理论模型试算和实际资料处理,实用效果比较好。 相似文献
8.
《物探化探计算技术》2015,(3)
针对目前地面直流电阻率法所面临的勘探深度有限的问题,开展了基于TV井地2.5D直流电阻率正则化反演研究。首先从井地直流电满足的边值问题出发,采用三角单元二次插值有限元法推导了2.5D井地电阻率正演公式。为了提高计算效率,对边界采用近似处理、以及利用图论理论的矩阵重排与填入元分析方法,实现了稀疏矩阵的直接解法,大大提高了计算效率。为了能够提高反演结果的稳定以及对异常边界分辨能力,采用L曲线法求取正则化因子、TV的稳定因子以及共轭梯度方法进行反演目标函数的求解,研究表明,反演结果高效、准确。最后对不同井地观测方式的反演结果进行对比分析,得出井地直流电阻率异常特征和分布规律。 相似文献
9.
冰川内部温度与冰面温度变化过程密切相关,因此可由实测冰温-深度剖面重建冰面温度变化。耦合热传导-冰流物理模型及其相关反演算法是基于冰温-深度剖面重建冰面温度变化的理论基础和关键。通过构造理想条件下的数值实验,分析比较了最小二乘法、Tikhonov正则化方法和蒙特卡罗算法重建百年尺度冰面温度变化的实验结果。同时,以我国藏北高原腹地的马兰冰川钻孔资料为例,结合气象资料构造真实钻孔数据模拟实验。分别给出了三种反演算法重建10 a和40 a尺度的冰面温度变化结果,并讨论了不同算法的优劣性和适用性。两组模拟实验结果均表明:Tikhonov正则化方法是目前求解该反问题的最优方法。与其他算法相比,Tikhonov正则化方法在各个时间尺度的重建结果与真实冰面温度变化吻合更好。另外,分别对重建百年和40 a尺度的冰温-深度剖面添加±0.02℃和±0.01℃的随机误差,研究发现Tikhonov正则化方法能够有效降低噪声干扰,相较其他算法得到的解更加稳定。这在一定程度上解决了该反问题中解的稳定性问题。 相似文献
10.
滑坡致灾范围的预测研究一直是滑坡研究的重点难点之一。以陕西泾阳南塬滑坡为研究对象,选取滑坡高度、体积、滑源区长度以及宽度为影响因子,采用蝙蝠算法对最小二乘支持向量机中的正则化参数γ和σ2进行寻优计算,建立BA-LSSVM滑坡致灾范围预测模型,并于多元线性回归模型进行对比。结果表明,该模型具有较高的预测精度和效果,可作为该地区防灾减灾依据。 相似文献
11.
12.
为评价采用抗滑桩支护后的边坡稳定性,从设计和验算两个角度着手,基于Morgenstern-Price（M-P）法建立了抗滑桩支护边坡的分析模型,进而得到了抗滑桩下滑力和边坡安全系数的表达式;通过引入自适应遗传优化算法,建立边坡稳定性分析优化模型,搜索采用抗滑桩支护边坡的非圆弧最危险滑动面;从确定安全系数求抗滑桩所受下滑力和确定抗滑桩抗力求边坡最小安全系数两种情况出发,探讨了抗滑桩位置对边坡稳定性的影响。结果表明：该方法可搜索出更加符合实际情况的边坡非圆弧最危险滑动面,并且能够得到抗滑桩的下滑力或边坡的最小安全系数,相同条件下抗滑桩应设置在边坡中部较为适合,才能最大发挥抗滑桩的加固作用。 相似文献
13.
线性方程组的正则化解法及其应用 总被引：3，自引：0，他引：3
在重力资料反演解释中，正则化法解线性方程组是一个行之有效的方法。这种方法的特点是引入正则化因子，求泛函极小和构组一族正则化近似并从中寻求最佳解。在求解过程中，要求对所求的解答预先有一个初步的推断解释。为此，要求对所解释的问题先进行一番研究，并对综合利用各种方法和手段所获得的先验资料进行对比，提出比较合适的基望值。本文给出了一个参数为已知的理论模型体的横剖面，与正则化法和其它线性方程解法的结果作对比，说明正则化法的有效性。文章也给出了利用正则化法对横穿吕宋岛、马尼拉海沟的一段地质地球物理剖面以及中国南黄海海域内─地质地球物理调查剖面进行反演解释的应用效果。 相似文献
14.
在重力资料反演解释中,正则化法解线性方程组是一个行之有效的方法。这种方法的特点是引入正则化因子,求泛函极小和构组一族正则化近似并从中寻求最佳解。在求解过程中,要求对所求的解答预先有一个初步的推断解释。为此,要求对所解释的问题先进行一番研究,并对综合利用各种方法和手段所获得的先验资料进行对比,提出比较合适的基望值。本文给出了一个参数为已知的理论模型体的横剖面,与正则化法和其它线性方程解法的结果作对比,说明正则化法的有效性。文章也给出了利用正则化法对横穿吕宋岛、马尼拉海沟的一段地质地球物理剖面以及中国南黄海海域内─地质地球物理调查剖面进行反演解释的应用效果。 相似文献
15.
为使传统的维纳滤波方法在最小均方准则下的解更准确、稳定,采用正则化思想,提出导数算子约束下的最小均方准则;在这种准则下,推导了维纳滤波因子及其频率响应.通过在不同频段取不同的正则化函数值,控制滤波结果.理论及实际地震勘探资料的分析表明:正则化维纳滤波去噪效果比传统维纳滤波好,且对有效波的损失较小,表明了文中所取正则项的可行性及有效性.最后讨论了通过建立正则化函数关系式,方便地提供正则化函数值以及圆滑振幅谱对改进滤波效果的可能性. 相似文献
16.
LSQR法具有计算效率高、对计算机内存要求低的优点,适合于大规模问题的求解。为探讨其应用于位场反演的稳定性和可靠性,笔者以加入不同噪声的两个合成模型数据为实验对象,比较分析了Tikhonov正则化与LSQR法求解结果,显示直接利用LSQR法求解位场反问题能够得到满意的正则化解,其解模型相对Tikhonov正则化,最大相对误差仅为0.36%,说明直接利用LSQR法求解位场反问题是可行的。将其应用于四川盆地雅安地区重力三维反演,极大地降低计算成本,获取了区内沉积盆及主要断裂分布情况,为页岩气靶区优选提供了有力支撑。 相似文献
17.
为消除TEM定源回线装置增大接收区域后明显的边框效应,且提高TEM的分层能力,达到提高对不同电性地层精细分辨的目的,采用带约束的正则化反演方法。在反演算法中,拉格朗日因子采用CMD正则化因子的自适应调节方案,自主调整。为保证反演过程的高效、稳定且结果收敛性,使反演模型每一层的厚度按照等对数间隔进行离散,减少矩阵方程的病态性;同时对反演电阻率值采用可行的上下限约束,使迭代的结果不至于出现极值或错值,降低反演的多解性。通过对多层理论模型的计算,证明了约束反演算法的正确性和可靠性,将其用于实测瞬变电磁法勘探数据的反演计算中,反演电阻率纵向分布规律与已知地层电性资料吻合,解释的异常区得到了钻探验证,结果可靠。 相似文献
18.
The iteratively re-weighted least squares (IRLS) is a commonly used algorithm which has received significant attention in geophysics and other fields of scientific computing for regularization of discrete ill-posed problems. The IRLS replaces a difficult optimization problem by a sequence of weighted linear systems. The optimum solution of the original problem is usually determined by computing the solution for various regularization parameters λ, each needing several re-weighted iterations (usually 10–15). In this paper, in order to decrease the required computation time (iterations) while maintaining good properties of the algorithm such as edge-preserving, the IRLS is augmented with a refinement strategy and the value of λ is progressively updated in a geometrical form during the iterations. The new algorithm, called iteratively re-weighted and refined least squares (IRRLS), can be interpreted as a Landweber iteration with a non-stationary shaping matrix which is updated based on the solution obtained from previous iteration. Two main properties of IRRLS are (1) the regularization parameter is the stopping iteration and (2) it is equipped with a tuning parameter which makes it flexible for recovering models with different smoothness. We show numerically that both the residual and regularization norms are monotone functions of iteration and hence well behaved for automatic determination of stopping parameter. The Stain’s unbiased risk estimate (SURE), generalized cross validation (GCV), L-curve analysis, and discrepancy principle (DCP) techniques are employed for automatic determination of optimum iteration. Experimental results from seismic deconvolution and seismic tomography are included showing that the proposed methodology outperforms the conventional IRLS with significantly lower computational burden. 相似文献
19.
瞬变电磁法（TEM）以其对低阻体反映灵敏而被广泛用于煤矿水文地质勘探中,当煤层埋藏较深,上覆地层纵向上有多种电性层存在时,瞬变电磁法晚期经验公式计算结果已不能反映实际情况。通过深入研究OCCAM反演算法,认为选择使误差目标函数达到最小的正则化参数,在反演迭代步长过大时进行压缩迭代步长的控制方法,以及最光滑的模型约束,可以提高反演分辨率。在瞬变电磁数据处理中使用加入了上述参数、控制和约束的反演技术,经对理论模型数据进行反演,反演结果与模型的电性分布规律吻合,表明算法正确;在华北煤田深部采场的瞬变电磁法探测应用中,OCCAM反演分辨出了煤层上覆多种电性地层,效果较好。 相似文献