共查询到10条相似文献,搜索用时 15 毫秒
1.
地球物理反问题线性化处理之后, 各种反演算法归结为对病态线性方程组的求解. 为了快速准确地计算出地球物理参数, 本文提出了一种全新的基于LSQR算法的混合差分进化算法(Hybrid Differential Evolution Algorithm, HDE). 该算法利用LSQR算法给出DE算法的初始种群, 提高DE算法的计算速度和稳定性. 在不同噪声水平下, 对四种正则化方法Tikhonov、TSVD、LSQR和HDE的反演结果进行详细比较. 理论模型和实际数据反演的结果都表明: 改进的HDE算法应用于地球物理反问题的求解是成功的: 反演结果与原设定模型具有较高的相关性, 在稳定性和准确性上较常规的反演算法都具有一定的优势; 而且不需要给定正则化参数, 具有更强的实用性. 相似文献
2.
本文将Tikhonov正则化方法与active-set算法相结合,利用双频电磁波电导率成像原理,求解其反演成像方程.不仅对现有算法进行了改进,也促进了算法的实际应用.本文研究了在双频电磁波电导率成像方程建立后,如何根据其严重病态性质,选择合适的算法求解矩阵成像方程.针对电导率非负的特性,引入正则化参数,将问题转化为一个非负最小二乘问题,并用active-set算法求解.采用改进后的迭代算法对理论模型进行了数值模拟计算,验证了该方法的有效性.应用到实际电导率成像反演,与常规的LSQR、SP-LSQR、Tikhonov正则化等算法进行比较,取得了满意的结果. 相似文献
3.
4.
5.
6.
本文将Bayes算法应用到Abel积分方程中,利用End-to-End Generic Occultation Performance Simulation and Processing System(EGOPS)软件仿真出的弯角数据来反演大气折射指数,并与Tikhonov正则化的反演结果进行对比.当直接利用仿真所得的弯角数据(认为它是不含误差的)时,Bayes算法与Tikhonov正则化的反演结果具有很好的一致性,二者的均方根误差均为2.5470×10~(-8);但在实际观测过程中,由于电离层以及水汽等方面的影响,所得到的弯角数据不可避免地会含有误差,可能产生高频成分,甚至存在间断的现象,因此我们在仿真的弯角数据中加入了满足高斯型分布的随机噪声,然后进行折射指数的反演试验.结果表明,相比于Tikhonov正则化技术,Bayes算法具有更高的反演精度. 相似文献
7.
8.
电磁波走时层析可用于空洞和孤石等离散块体异常的探测中.在传统电磁波走时层析中,通常采用一个低阶差分泛函来稳定反演过程,而这类泛函都具有光滑约束的性质,在成像结果中难以对此类离散块体异常进行准确的解释.本文基于Tikhonov正则化思想,提出了一种反余切泛函,该泛函具有紧凑约束的性质且无需选择一个聚焦因子.结合弯曲胖射线层析理论和重加权正则化共轭梯度反演算法,对两个理论模型进行了成像效果研究.与低阶差分泛函和最小支撑泛函的成像结果相比,该方法能够更好地突出异常的边界,成像结果中的虚假异常也要更少.同时,还分析讨论了激发电磁波的中心频率以及数据噪音对成像结果的影响.此外,针对一个地下连续墙墙体内部缺陷的现场模型,利用该方法取得了理想的成像结果,良好地揭示了缺陷的形态和位置. 相似文献
9.
10.
《地球物理学进展》2016,(1)
本文基于重力梯度张量密度反演基本理论,建立了模型约束正则化密度反演矩阵方程.分析了代数重构算法(ART)中迭代初始值、松弛因子和终止条件三个关键参数的影响;与最小二乘求逆法对应比较分析了算法的时间和精度.结果表明:在地震、地质等地球物理手段提供初值、边界等约束较多的情况下,ART可以克服方程的不适定进行直接求解,并且合理的松弛因子和终止条件可有效提高反演效率.当初始信息不足时,添加光滑假设、深度加权等模型约束,正则化方程可以提高反演结果的可靠性.ART的行迭代可有效避免观测误差的积累和矩阵求逆的计算,从而使计算精度和速度提高数倍.最后基于GOCE地球重力场模型所得重力梯度,以地震层析成像所得速度模型为约束,对华北克拉通密度结构进行了反演,并与该区已有密度研究结果进行了对比.结果表明:利用GOCE重力场系数计算重力梯度扰动,以速度模型为约束,基于代数重构算法进行重力梯度反演所得密度模型与重力-地震联合反演所得密度模型具有很好的对应性.ART算法为重力梯度张量反演中大规模复杂问题的快速计算提供了又一种有效手段. 相似文献