共查询到20条相似文献，搜索用时 890 毫秒
1.
反射系数与子波同时迭代反演的预条件共轭梯度法 总被引：4，自引：0，他引：4
在反射系数白噪、子波最小相位的假设下,研究基于线性反演的地震反射系数和子波同时估计问题。在Cauchy准则稀疏反演求解中,应用预条件共轭梯度法实现反射系数和子波同时迭代反演。在迭代求解正则化方程时,用共轭梯度法求解相应的原问题,初猜子波求解也使用该策略。模型数据试算与比较,表明了该算法正确而有效。用实际数据检验算法的实用性,经研究表明,预条件共轭梯度法计算的反射系数和子波,要比直接稀疏反演精度高,而且收敛较快,数值稳定,实用性强。 相似文献
2.
常规磁性基底反演方法往往忽略基底磁性变化采用常磁化强度模型,这里分析了其中的局限性,提出了先构建相关搜索重磁三维定量反演技术的变磁化强度模型,然后通过Parker界面反演算法进行基底反演.其中,对反演磁化强度所需磁异常分离采用了小波变换方法,并分析了该方法的关键技术和措施.通过利用该方法对某工区进行了实际资料处理后,得到较好效果. 相似文献
3.
《物探化探计算技术》2015,(6)
利用阻尼型高斯牛顿法求解反演目标函数,并将粗糙度矩阵约束和最小支撑泛函约束以及参考模型和模型参数界限约束引入到目标函数中,同时结合预条件共轭梯度法进行反演迭代,实现了激发极化数据三维聚焦反演快速计算。通过对几个典型模型的试算,并与传统光滑模型约束的反演结果进行比较,表明该反演算法的可靠性和稳定性较好,并且算法速度快,占用内存低且易于并行化。 相似文献
4.
基于稀疏离散τ-p变换的叠后地震道外推 总被引：2，自引：0，他引：2
这里研究的是一种基于稀疏离散τ-p变换的地震道缺失重建方法。该方法根据叠后零偏移距剖面,在局部时窗内可以看作是一系列线性同相轴的组合,使用稀疏离散τ-p和预条件双共轭梯度算法进行地震道外推,使空间方向的缺失道得到恢复。τ-p变换由于信息不足,有限的孔径和离散等因素,导致在τ-p域的结果存在假象、不准确,而在τ-p域进行稀疏对于重建缺失的信息是非常有用和必要的。同时,还可使用迭代的预条件双共轭梯度算法进行计算。 相似文献
5.
全波形反演是一种全新的地震成像方法,主要利用全波形信息反演地下介质参数,通过非线性优化波场理论值和观测值的残差实现波形反演。基于时间域声波方程,建立了时间域波场残差目标函数,以分层模型为例,分别从波场精度、目标函数收敛性和运行时间3个方面,比较了共轭梯度（Conjugate Gradient,CG）算法和拟牛顿算法（Broyden-Fletcher-Goldfarb-Shanno,BFGS）反演的效果。同时,应用共轭梯度法对正、逆断层模型和Marmousi模型进行了速度结构反演。反演结果表明:共轭梯度法计算效率较高,反演得到的速度模型精度更高,反演效果较好,是一种有效的波形反演方法。 相似文献
6.
随着地球物理设备和探测技术的不断发展,快速处理大规模地球物理数据的需求也随之增长。为了解决三维重力数据密度反演的耗时问题,提出一种并行的预处理共轭梯度算法来提高计算效率。本文分别采用两种不同的预处理算子通过组合模型数据反演进行测试比较,并利用迭代残差和计算用时共同评价其加速效果。结果表明：对称逐次超松弛预处理方法比对角预处理方法反演计算速度快,密度结果更贴近实际模型;与传统串行的共轭梯度算法相比,本文并行预处理快速算法可以获得近19倍的加速比。将该算法应用于美国Vinton盐丘的实测重力数据中,反演结果能够很好地圈定出岩体的位置,验证了本文并行预处理共轭梯度法在三维重力数据快速反演中的高效性和可行性。 相似文献
7.
通过将地下介质划分为一定规模的网格单元,并构造一个与磁性体中心网格位置有关的指数函数模拟地下介质磁性分布特征,利用粒子群算法反演该函数参数,根据求得的函数表达式,推算各网格单元磁化强度,由此实现磁化强度的快速成像。理论模型试验显示:该方法不仅能将一个求解20×10网格规模的物性分布问题简化为求解指数函数的10个参数值的问题;而且,当在不同模型地质体磁异常中加入20%的随机噪声时,也能取得理想的效果,且需要调整的参数少,算法通用性强。利用该算法处理老挝爬立山铁矿区某磁测剖面,得到厚度近140 m的强磁性体;钻孔证实,实际钻遇矿体累积最大厚度达130.6 m,且矿体中心与成像中心相对应:由此证明该方法应用于实际资料处理与解释的可行性和有效性。 相似文献
8.
起伏地形电导率连续变化的三维激电数据有限元数值模拟 总被引：2，自引：1，他引：1
为了模拟起伏地形条件下的地电体,这里给出了电导率连续变化的三维激电数据有限元数值模拟方法.在模拟过程中,采用超松弛预条件共轭梯度法解线性方程组,并探讨了模拟的计算精度和效率问题.最后,对起伏地形下的地电体进行模拟,验证了本文算法的正确性和可行性. 相似文献
9.
为了适应卫星磁测资料解释的需要,作者在磁性界面和磁性层磁异常正反演方法研究的基础上,首次提出了球壳形磁性层磁异常正反演理论和方法。 作者首先阐述了垂直磁化球壳型磁性层模型,导出了场垂直分量Z_(λο)的球谐系数正演公式。然后,给出了磁性层上,下层面埋深变化值与磁化强度变化值的球谐系数反演公式。进而导出了迭代反演收敛条件。本文的磁化强度反演方法不同于目前国内外已有的方法,其反演对象是磁化强度变化值△J_(λθ),磁化强度平均值J是必须给定的定解条件。这种反演方法能给出尽可能满足人们解释要求的解,克服了前人反演方法总存在负值的困难。 为了增加球壳形磁性层反演方法的实用性,作者还在类似平面频率城反演方法所作假设的前提下,把这种利用球谐分析理论在全球面导出的方法,严格地推广应用于局部区域。这也是球谐分析理论在局部区域应用的一个实例,在理沦上具有一定的意义。 相似文献
10.
新疆色尔特能地区视密度和视磁性的反演 总被引：1，自引：0，他引：1
将起伏地形上测量的重力异常和磁异常,延拓至平均水平面;对平面上的场进行不同深度层源的切割分离;再将各层的场向下延拓至相应的深度,并反演出该深度层的密度和磁化强度结构。以此得到的密度和磁化强度,反映该深度层密度和磁性的近似分布,称为视密度和视磁性反演。本方法反演速度很快,在主频1.73 GHz的微机上,完成71×81数据点反演的计算时间为10s。这里以新疆色尔特能地区的实例说明这种方法的反演效果。 相似文献
11.
电阻率三维反演稳定性和可靠性研究 总被引：1，自引：0，他引：1
电阻率三维反演参数太多,将导致反演不稳定性和高度非唯一性等诸多方面的问题。比较传统的最小二乘反演方法,利用共轭梯度的电阻率三维最小构造反演显示了更好的效果。讨论了初始模型选择及数据误差对三维反演结果的影响,说明了反演方法的稳定性和可靠性。 相似文献
12.
遗传算法是近些年来产生和发展的一种模拟生物进化过程的自适应启发式全局优化的搜索算法。它不完全依赖于初始猜测,且具有全局收敛的特点,可以被用来解决各种复杂的实际问题,如工程优化设计,人工智能和决策系统,以及地球物理反演等。尽管遗传算法是一种效率很高的全局优化算法,但许多仿真结果表明,它具有计算时间长,局部搜索能力弱的缺点。而共轭梯度法属于非启发式全局优化搜索方法,收敛速度快,但容易陷入局部极值,且严重依赖初始猜测。根据遗传算法和共轭梯度法的特点,这里提出了一种混合遗传算法,用来进行地球物理反演。该算法既具有遗传算法的全局收敛性,又有共轭梯度法的快速收敛性,经实际应用,取得了良好的效果。 相似文献
13.
用改进的共轭梯度算法解波阻抗反演问题 总被引：2，自引：0，他引：2
研究快速有效的波阻抗反演算法,是地震反演中一直追求的目标。这里从目标函数出发,推导的反演方程,在常规共轭梯度算法的基础上,对算法搜索的共轭梯度方向提出了具体的改进措施,并在阻尼项及松弛因子的给定上实现了自适应化。另外,在整个算法过程中实现了向量运算,无需额外的内存。通过合成记录及实测资料试算对比发现,该算法具有精度高,运算速度快及抗病态能力强的特点。 相似文献
14.
《物探化探计算技术》2020,(2)
磁化强度矢量反演方法在磁法勘探中具有现实意义,能够获取磁场强度大小和磁化方向。磁化强度矢量反演方法通常是基于笛卡尔直角坐标系进行计算,在野外和实验室测得的岩石磁性数据都是以磁化强度和磁化方向的形式给出,将其变换到笛卡尔直角坐标系中进行计算,显然增加了计算的冗余性。这里首先通过模型实验证明球坐标系下的磁化强度矢量反演方法具有两个优点:①减少了未知数的个数,因而可以减少计算时间;②反演结果更加聚焦,深度分辨率更高。然后将球坐标系下的磁化强度矢量反演方法应用于卡拉塔格红海-黄土坡勘查区,对该区的实测航磁数据进行处理,计算得到了勘查区的磁化率模型,为深部找矿提供地球物理借鉴。 相似文献
15.
实现稀疏反褶积的预条件双共轭梯度法 总被引：9，自引：3，他引：9
地震勘探稀疏反褶积计算一般要导出一个Toeplitz矩阵的线性系统，通常可以用矩阵求逆、Levison递推及共轭梯度等方法直接求解。当Toeplitz矩阵的条件数很大时，数值稳定性差，甚至无法求解。使用共轭梯度法，在矩阵的对角元素上加入规则化因子，可以改善这种情况，但不能彻底解决数值稳定性和精度问题。若求解最小二乘问题的原始问题，结果会好些。线性系统形式的细微改变，将导致不同的数值计算特性。在规则化策略基础上，可巧妙地构造稀疏反褶积的问题原型，引入预条件，采用双共轭梯度法求解，从而实现稀疏反褶积，获得较好结果。数值算例表明，预条件双共轭梯度法比直接稀疏反褶积方法收敛快、精度高。 相似文献
16.
17.
18.
针对充电电场反演计算参数较多以及复形调优算法处理约束条件的直接有效性,作者就均匀充电矩形薄板反演问题编制程序,计算地质体理深、倾角及大小。尽管迭代次数较高,计算结果颇为理想,取得了较好的效果。 相似文献
19.
遗传算法是模拟自然界生物进化思想所得到的一种全局非线性最优化算法,该算法可有效地避免局部最优解。这里介绍遗传算法的原理以及基于其算法的磁异常反演算例,对不同迭代次数所得理论模型反演的磁性体空间位置与磁性等参数的相对误差进行了对比。据此,对四川会理黎溪地区某磁铁矿区外围8号剖面上实测△T磁异常进行遗传算法,并迭代拟合反演,求其反问题,获得了与后续工程揭露及磁性测定结果相吻合的结果。 相似文献
20.
构建近地表横波速度模型是煤田多分量地震资料处理的重要环节。相较于面波多道分析法,全波形反演在构建近地表横波模型中具有更高的分辨率。然而,在基于梯度的全波形反演中,由于地震记录频带有限、波场的非均匀覆盖以及双重散射等原因导致梯度算子不随深度的增加而缩放,模型深部参数得不到明显更新。目标函数的Hessian算子包含曲率信息,可清晰预测梯度算子中的焦散现象及双重散射产生的伪影,因此,逆Hessian算子则可作为反卷积算子实现对梯度的预处理,加强对模型深部的照明能力。然而Hessian算子具有巨大维度,对其显式计算十分困难。基于此,借鉴逆散射理论的思想,给出勒夫波全波形反演目标函数的拟Hessian算子的表达式,并提出一种梯度预处理的全波形反演方法。将该方法分别应用于断层模型、凹陷模型以及起伏界面模型的重构试验,反演结果表明:与传统的共轭梯度全波形反演方法相比,基于拟Hessian算子的预处理共轭梯度方法可加快收敛速度,提升成像质量。 相似文献