共查询到20条相似文献,搜索用时 31 毫秒
1.
2.
如何得到快速稳定的反演结果和更清晰的地质体分界面等问题仍然是当前MT反演研究的一个重点.为了解决反演结果不能得到清晰的电性分界面的问题,本文在前人研究的基础上,基于OCCAM反演以及突出尖锐边界反演的思想,结合最小支撑梯度泛函,构建了新的反演目标函数,并利用共轭梯度法优化目标函数,实现了一种光滑模型与尖锐边界结合的MT二维反演方法.模型实验证明了该方法的准确性,通过与已发表的相关反演方法的结果进行比较,证明了该方法在光滑稳定的基础上可突出对尖锐电性边界的刻画.对广东徐闻地区的实测MT资料进行了处理,表明了该方法的适用性和效果. 相似文献
3.
4.
尖锐边界反演是大地电磁及其他地球物理反演中的一个较难解决的问题,目前有不少专家正对这方面进行研究.我们在分析了Zhdanov(2004)提出的基于最小支撑泛函聚焦反演方法的基础上,在目标函数中引入对角梯度支撑以改善倾斜电性分界面的反演效果.同时,为了提高计算效率,本文还使用互易定义计算灵敏度矩阵及使用正则化共轭梯度法优化目标函数.通过模型试验并与未引入对角梯度支撑的反演结果进行了比较,表明本文的方法反演效果更好而且算法稳定和实用. 相似文献
5.
6.
电磁波走时层析可用于空洞和孤石等离散块体异常的探测中.在传统电磁波走时层析中,通常采用一个低阶差分泛函来稳定反演过程,而这类泛函都具有光滑约束的性质,在成像结果中难以对此类离散块体异常进行准确的解释.本文基于Tikhonov正则化思想,提出了一种反余切泛函,该泛函具有紧凑约束的性质且无需选择一个聚焦因子.结合弯曲胖射线层析理论和重加权正则化共轭梯度反演算法,对两个理论模型进行了成像效果研究.与低阶差分泛函和最小支撑泛函的成像结果相比,该方法能够更好地突出异常的边界,成像结果中的虚假异常也要更少.同时,还分析讨论了激发电磁波的中心频率以及数据噪音对成像结果的影响.此外,针对一个地下连续墙墙体内部缺陷的现场模型,利用该方法取得了理想的成像结果,良好地揭示了缺陷的形态和位置. 相似文献
7.
为推进大地电磁三维反演的实用化,本文实现了基于L-BFGS算法的带地形大地电磁三维反演.首先推导了大地电磁法三维反演的Tikhonov正则化目标函数以及Hessian矩阵逆矩阵近似表达式和计算方法,然后设计了一种既能保证空气电阻率固定不变又能保证模型平滑约束的协方差矩阵统一表达式,解决带地形反演问题.在反演算法中采用正则化因子冷却法以及基于Wolf条件的步长搜索策略,提升了反演的稳定性.利用开发的算法对多个带地形地电模型(山峰地形下的单个异常模型、峰-谷地形下的棋盘模型)的合成数据进行了三维反演,并与已有大地电磁三维反演程序(ModEM)进行对比,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了华南某山区大地电磁实测数据,得到该区三维电性结构,揭示了研究区以高阻介质为基底,中间以低阻不整合面和相对低阻介质连续分布,浅部覆盖高阻介质的电性结构特征,进一步验证了本文算法的实用性. 相似文献
8.
9.
10.
11.
《应用地球物理》2015,(2)
<正>则化反演通过引入模型约束和正则化因子求解病态的地球物理反演问题,但该方法存在正则化因子选取困难和初始模型依赖的问题。针对该问题,本文提出多目标粒子群反演算法。该算法反演中不需要目标函数梯度信息和正则化因子,先同时求数据拟合和模型约束的多目标反演解集,再权衡两者的相对重要程度,最后从反演解集中优选出最终反演结果,从而起到正则化因子的作用。以二维磁测数据反演为例,进行理论模型反演试验,试验结果表明,多目标粒子群反演算法能尽可能多地保留可行解,得到反演解集;通过分析反演解集,既能深入的理解反演过程,又能灵活地从数据拟合和模型约束两方面进行权衡与选择,得到比正则化反演更合理的反演结果;该算法能同时解决正则化因子选取困难和初始模型依赖问题。 相似文献
12.
《地球物理学进展》2016,(3)
航空电磁探测数据量大,二维、三维反演算法复杂、计算速度慢,通常采用一维反演,利用层状模型拼接描述地下复杂结构,但航空电磁数据信噪比低,容易引起一维反演结果横向连续性差等问题.本文针对上述问题,基于一维反演算法,通过整合测线观测数据,建立了测线数据整体的目标函数,并根据Tikhonov正则化反演理论,引入包含空间粗糙度和先验信息的模型参数约束项,确定了拟二维整体反演的目标函数,推导了反演迭代方程组,利用超松弛共轭梯度算法,求得由于整条测线整体反演所致的大型稀疏矩阵的极小化解,实现了对整条测线数据同时反演的固定翼航空电磁数据的拟二维整体反演算法.在反演迭代过程中,正则化因子采用线性搜索自适应迭代的方法自动选取,提高了反演结果的稳定性.对比分析了仿真数据的一维反演与拟二维整体反演结果,得出拟二维整体反演算法横向连续性较好,对高导覆盖层下的导体分辨率优于一维反演,同时受高斯噪声的影响较小.最后,将直升机飞行实测噪声加入仿真数据中,拟二维整体反演结果平均相对误差较一维反演结果降低了31.6%,进一步验证了拟二维整体反演算法的有效性. 相似文献
13.
14.
可控源音频大地电磁数据的反演方法 总被引:18,自引:7,他引:18
从反演方程、构造目标函数和求解三方面对用于可控源音频大地电磁法(CSAMT)的实用反演方法中的四种进行了描述.水平层状地层CSAMT法资料的直接反演法首次尝试了一维空间的全资料CSAMT反演,效果较好,但该方法尚难应用于2D、3D复杂介质中;奥克姆反演方法既考虑了横向的光滑函数,又考虑了纵向的光滑函数,得到比较光滑的横向、纵向变化的背景电性结果,但有可能把一些小构造光滑掉.快速松驰反演算法和共轭梯度算法由于计算速度快,占内存少而被用于三维反演中,二者相比,快速松驰算法在求解雅可比矩阵时只做一次正演计算,在更新模型时解小型方程组,所以在速度上更胜一筹.在后三种算法中,由于复杂电性结构无解析解,正演计算都采用数值计算.数值计算的可靠性、速度影响着反演算法的有效性,这方面的研究也将是2D、3D复杂电性结构反演的研究方向之一. 相似文献
15.
16.
《应用地球物理》2015,(3)
频率域全波形反演充分利用全波场的振幅、相位以及频率信息,采用较少的频率便能反演得到精度很高的速度模型。本文以有限单元法为基础,对起伏地形条件下二维声波频率域全波形反演进行了研究。在正演算法中,针对截断边界问题,并考虑多频率联合反演中计算区域采用同一套剖分网格的需求,提出了一种适用于起伏地形的衰减边界条件算法。该算法的核心思想是在控制方程波数项中引入衰减因子,通过一定方式调节衰减因子使得声波在衰减层中充分衰减,达到压制截断边界影响的目的。根据指数衰减规律,文中推导出了一种新的衰减因子计算公式,并给出了不同频率条件下衰减层厚度计算公式;在反演算法中,采用共轭梯度法求解高斯牛顿反演迭代方程组,避免直接求解雅克比矩阵和HeSSian矩阵带来的巨额计算量,并采用相同的反演模型,对比分析了不同初始模型和频率组合对全波形反演结果的影响。起伏地形模型数值模拟和全波形反演数值试验表明,本文提出的指数衰减边界条件算法和基于该算法的全波形反演算法具有很好的应用效果。 相似文献
17.
18.
为了克服时间域航空电磁数据单点反演结果中常见的电阻率或层厚度横向突变造成数据难以解释的问题,通过引入双向约束实现航空电磁拟三维空间约束反演.除考虑沿测线方向相邻测点之间的横向约束外,同时还考虑了垂直测线方向测点在空间上的相互约束.为此,首先设计拟三维模型中固定层厚和可变层厚两种空间约束反演方案,然后通过在目标函数中引入沿测线和垂直测线方向上的模型参数约束矩阵,并使用L-BFGS算法使目标函数最小化,获得最优拟三维模型空间反演解.基于理论模型和实测数据反演,对单点反演与两种空间约束反演方案的有效性进行比较,证明本文空间约束反演算法对于噪声的压制效果好,反演的界面连续光滑,同时内存需求和反演时间少,是一种快速有效的反演策略. 相似文献
19.
研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法(BICG)计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次“拟正演”问题,进而将PARDISO的快速正演方案应用于“拟正演”问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性. 相似文献