共查询到19条相似文献，搜索用时 234 毫秒
1.
2.
航空瞬变电磁数据一维Occam反演 总被引：1，自引：0，他引：1
分析了奥克姆(Occam)反演法的原理和计算方法,并将其应用于时间域航空电磁一维反演中,给出了一种简便的正则化参数搜索方法,对迭代过程进行了改进,大大提高了反演速度,5次迭代约20s。在反演初始模型中,电阻率数据由原始感应电动势数据通过自动迭代法转换而来,深度数据由初始电阻率计算得出。经算例表明,收敛过程中满足拟合差最小的拉格朗日乘子μ值随迭代次数增加而递减;反演过程非常稳定,一般5次迭代至10次迭代就能达到收敛,对低阻层的反演深度和电阻率值都接近真实值,但对高阻层的反演深度和电阻率值误差较大。 相似文献
3.
二维视电阻率断面的快速最小二乘反演 总被引：12，自引：0，他引：12
快速最小二乘反演是以平滑限定的最小二乘方法为基础,是对二维视电阻率断面进行反演的一种方法。反演过程不需要提供初始模型,在首次迭代时使用一均匀介质地下模型作为初始模型,该模型的视电阻率偏导数值可以用解析法得到。在后面的迭代中,使用了拟牛顿法去修改每一次迭代的偏导数矩阵,避免了偏导数矩阵的直接计算,从而减少了计算时间和存储空间。同时运用牛顿矩阵校正技术解最小二乘方程组也减少了大量的计算时间。总之,该方 相似文献
4.
有限单元法三维电阻率最小二乘反演中存在问题的研究 总被引：6，自引：1，他引：5
针对有限单元法在三维电阻率最小二乘反演中关于深部单元的反演能力、反演稳定性和速度问题进行了讨论与改进.首先构造了全新的"体积因子"作为先验信息加入到电阻率反演中,使地质体赋存区域的电阻率变化更符合实际情况,提高了深度较大的地质体的分辨能力.反演计算中对实测电阻率数据和模型正演数据分别取对数,并改进迭代修正量的计算方法,将以前的加减修正量改为与修正量指数之积,这样,不但可以改进电阻率模型出现负值的缺陷,提高反演计算的稳定性,而且还可减少迭代次数;针对反演时间过长的缺陷,将三维全域反演改为三维局域反演,使每次迭代的时间由4 h减少到目前的15 min;对模型的反演计算表明,以上改进是非常有效的. 相似文献
5.
利用传输线理论、Sommerfeld积分快速计算以及最小平方拟合技术研究建立多分量感应测井数据的一种新的快速参数化迭代反演算法,同时重构水平层状横向同性地层的纵、横向电阻率以及水平层界面深度。首先,通过Fourier变换与传输线理论给出频率波数域中电磁场并矢Green函数在各个地层中的解析解,并利用三次样条插值和贝塞尔函数递推公式建立Sommerfeld积分的半解析算法,快速计算多分量感应的测井响应。然后在此基础上,利用摄动理论建立磁场并矢Green函数与模型向量间变化关系的摄动方程,并将摄动方程中各个积分转化为Sommerfeld积分,实现正演模拟的同时用半解析算法快速确定多分量感应测井响应的Fréchet导数。最后,利用归一化处理和奇异值分解技术,同时反演所有地层的纵、横向电阻率和层界面深度,实现输入数据和反演模型的模拟数据优化拟合。理论模型的数值结果验证了该反演算法的有效性及抗噪性。 相似文献
6.
7.
通过一维层状介质电磁法数据的时域和频域正演结果,对不同电阻率大小的电阻层进行了反演计算,说明了在这种情况下电磁法数据时频联合反演的合理性.正演计算了一维层状模型下水平电偶极频率域水平电场响应及利用G-S变换计算时间域垂直磁场导数响应;利用最光滑模型约束反演方法建立了时频联合反演迭代公式.通过理论模型反演试算表明,该反演效果较好. 相似文献
8.
9.
10.
基于磁导率的频率域航空电磁法双频反演方法 总被引：1，自引：0，他引：1
以往的频率域航空电磁法数据反演方法中,通常在模型中以真空磁导率代替介质磁导率,忽略了磁性层影响。但在高磁性区域,这样做将会导致计算电阻率值偏高。针对这一问题,首先根据层状介质模型研究了考虑磁导率和电阻率因素条件下的电磁场正演计算方法,并且分析了不同磁性介质对频率域电磁响应的影响规律。在此基础上,提出了磁导率和电阻率的双频反演方法,电阻率是通过2个频率的虚分量之比计算得到的,计算结果更加可靠,提高了磁性区域视电阻率填图的精度。通过对模型计算和实测数据进行反演计算,表明这种方法具有较高的反演精度,为强磁性区域航空电磁法视电阻率填图提了供一种行之有效的反演方法。 相似文献
11.
三维地电断面激发极化法有限元数值模拟 总被引：5，自引：0，他引：5
用有限单元法求解三维地电断面激发极化法正演模拟算法.首先给出了三维构造中点源电场异常电位的边值问题与变分问题, 简化无穷远边界上的边界条件以提高计算速度及计算精度.以此为基础计算视电阻率对模型电阻率的偏导数矩阵, 并进行三维地电断面激发极化法正演模拟, 与等效电阻率法相比节约了计算时间.对几例较典型的地电模型进行计算, 结果表明本方法是正确可行的. 相似文献
12.
电阻率三维复杂结构的快速反演 总被引：3，自引：1，他引：2
运用共轭梯度迭代算法解三维反问题中的线性方程组,并结合求Jacobi矩阵 G 的Rodi算法,则每次反演迭代仅需一次正演计算,大大加快了计算速度,实现了直流电阻率三维快速反演。另外,由于避免了存储 G 和 G T G 所需的庞大存储量,以及在三维反演中加入光滑约束,有利于精细网格下的复杂模型反演。 相似文献
13.
讨论了应用Ｋｏｅｆｏｅｄ法自动直接解释电测深曲线中的数据采集、归化采样值的计算、反演计算、正演计算、图形显示以及计算精度等问题，并介绍了由作者开发的一套完整的电测深自动解释软件。 相似文献
14.
边界积分方程法二维电阻率层析成像 总被引：1，自引：1，他引：0
本文基于作者在文献［１］中提出的边界积分方程,建立了一种新的二维电阻率层析成像方案。与目前国内外用有限元法和α中心法的成像方法相比,该方法无需计算雅可比矩阵,也不必做线性化近似。这不仅节省了计算时间,而且使成像迭代过程稳定,可以对实际上广泛存在的大电性对比度介质获得满意的成像效果,从而大大提高了电阻率层析成像技术的实用性。文中给出的数值模拟结果显示了该方法的有效性。 相似文献
15.
三维电阻率成像新方法及应用 总被引：4，自引：0，他引：4
三维电阻率成像法包括近似三维反演因子和综合子空间解释器,三维反演因子用于正向线性电阻率近似解译,综合子空间法是共轭梯度法的变换,是求解大型最优化问题的有效方法.三维成像是线性反演问题,不需正演模拟或敏感度修改.在非线性三维反演中,三维电阻率成像技术可直接获取三维电阻率的分布信息,或提供中间过程的模式修改.野外资料验证表明,三维成像技术可提供电阻率空间的分布信息. 相似文献
16.
这里引入在工程电磁场（EEF）中使用的新型等效源法（NESM）,并给出了新型等效源法模拟感应测井响应的理论推导及处理过程,模拟了感应测井在不同地质条件下的响应。新型等效源法是在广义多级技术（GMT）的基础上发展起来的,具有理论推导简单,计算速度快,计算精度高等优点。新型等效源法是半解析解法,占用计算机内存少。另外,解的可导性和光滑性比有限元法更胜一筹。采用新型等效源法对感应测井的响应进行数值模拟,并把数值模拟的结果和数值模式匹配法（NMM）进行对比,验证了新型等效源法可行性,并且精度可以随着等效源和级数项的增加而提高。新型等效源法形成的矩阵条件数比较大,应用双精度求解可以增加稳定性。 相似文献
17.
Richards方程常用于非饱和土渗流问题,并且应用广泛。在数值求解中,对Richards方程线性化,进而采用有限差分法进行数值离散以及迭代计算。其中传统的迭代法比如Jacobi迭代、Gauss-Seidel迭代法（GS）和连续超松驰迭代法（successive over-relaxation method,简称SOR）迭代收敛率较慢,尤其在离散空间步长较小以及离散时间步长较大时。因此,采用整体校正法以及多步预处理法对传统迭代法进行改进,提出一种基于整体校正法的多步预处理Gauss-Seidel迭代法（improved Gauss-Seidel iterative method with multistep preconditioner based on the integral correction method,简称ICMP(m)-GS）求解Richards方程导出的线性方程组。通过非饱和渗流算例,并与传统迭代法和解析解对比,对改进算法的收敛率和加速效果进行了验证。结果表明,提出的ICMP(m)-GS可以很大程度地改善线性方程组的病态性,相较于常规方法GS,SOR以及单一改进方法,ICMP(m)-GS具有更快的收敛率,更高的计算效率和计算精度。该方法可以为非饱和土渗流的数值模拟提供一定参考。 相似文献
18.
利用一维正则化反演进行大地电磁测深数据拟二维反演解释 总被引：1，自引：0，他引：1
大地电磁测深的反演问题是不适定的,其反演结果不稳定,且具有非唯一性。通过在目标函数中采用正则化方法,可以使得不适定反演问题具有稳定的反演结果,并改善解的稳定性和非唯一性问题。为了提高野外大地电磁测深数据的处理效率和初步解释的精度,提出了大地电磁测深数据的一维正则化反演进行拟二维反演解释方法。这里所述的大地电磁测深一维反演解释,与以往的解释方法不同,其思路首先用Bostick反演的深度来控制层参数,使反演计算的模型参数仅存在电阻率;最后采用阻尼高斯-牛顿算法进行反演计算,并将Bostick反演结果作为反演计算的初始模型。通过模型试算,结果表明其处理速度快、解释直观,对野外大地电磁测深数据进行初步反演解释是可行的。 相似文献
19.
为提高起伏地形三维激电数据的反演解释精度和处理效率,笔者开展了连续介质模型快速反演方法研究。针对起伏地形连续电性介质模型,提出了四面体单元电性参数分块线性连续变化的参数化方法。在反演中对模型参数施加光滑和背景约束信息,以提高反演的稳定性和分辨率。通过采用降低反问题的维数、压缩存储线性反演方程各矩阵的元素以及将互换原理和拟牛顿法相结合的方式计算偏导数矩阵等手段,可有效地加快反演的计算速度。最后,为验证反演方法的有效性,对2例地电模型进行反演试算,计算结果表明：反演耗费时间较少,仅迭代6次拟合差便趋于稳定,反演结果能较好地刻画异常体形态,编制的反演解释软件可用于实际生产。 相似文献