共查询到20条相似文献,搜索用时 390 毫秒
1.
可控源音频大地电磁法在资源勘探等领域中发挥着重要的作用.我们把有限差分数值模拟方法用于可控源音频大地电磁三维正演,结合正则化反演方案和共轭梯度反演的思路,将反演中的雅可比矩阵计算问题转为求解两次"拟正演"问题,得到模型参数的更新步长,形成反演迭代,实现了可控源音频大地电磁三维共轭梯度反演算法.该反演算法可用于对有限长度电偶源激发下采集到的可控源音频大地电磁全区(近区、过渡区和远区)视电阻率和相位资料进行三维反演定量解释,获得地下三维模型的电阻率结构.理论模型合成数据的反演算例验证了所实现的可控源音频大地电磁三维共轭梯度反演算法的有效性和稳定性. 相似文献
2.
最近开发了一种针对二维大地电磁野外数据进行处理解释的新反演方法。该方法以加入阻抗相位信息的一维大地电磁连续介质曲线对比法为基础 ,把一维反演得到的电阻率和相位的数据集作为二维反演的初始模型 ,使用二维有限单元法做正演模拟。在程序的后继迭代中 ,深度方向上用一维反演修改模型的电阻率和深度值 ,沿测线方向由二维有限元作修改 ,反演结束可得到一个接近真实电性分布的电阻率数据集 ,并绘制成电阻率断面图。对模型的反演实验结果显示 ,该反演方法能够较真实地反映地下电性分布 ,而且避免了偏导数矩阵的计算 ,其原理简单 ,计算速度快 ,表明该反演方案是可行的 相似文献
3.
航空电磁数据的三维解释由于数据量大需要有高效的反演算法作为支撑.本文利用两种目前主流的数值优化技术(非线性共轭梯度和有限内存的BFGS法)实现了三维频率域航空电磁反演,并进一步比较了两种方法的有效性和运算效率.在反演过程中,为了更好地反演异常体的空间位置,模型方差矩阵中的光滑系数在反演起始阶段取值较大;当数据拟合差下降趋于平缓时,再利用较小的光滑因子约束反演过程来实现聚焦和获得精确的反演结果.理论数据反演表明这两种优化策略具有相似的内存需求,但是有限内存的BFGS技术比非线性共轭梯度法在计算时间和模型反演分辨率上具有一定的优越性,因此有限内存BFGS法更适合于求解大规模三维反演问题. 模型试验进一步表明目前主流的迭代法求解技术不适合大规模航空电磁数据反演,未来移动平台多源电磁数据快速正反演可通过引入矩阵分解技术来实现. 相似文献
4.
推导出频率域有限元声波正演方程,为了消除边界反射,将Clayton-Engquist旁轴波动方程吸收边界条件引入频率域,并对有限元刚度矩阵和质量矩阵进行压缩存储,利用广义共轭梯度法求解有限元方程获得正演解.在此基础上,推导出在某一频率下波场数据残差δU与单元物性参数修改量δλ之间关系的Jacobi矩阵,反演方法允许利用地面二维炮集全波场资料与给出初始模型参数的正演值的差值δU,迭代求得δλ.由于计算机内存的限制,方法计算不允许有过多数目的未知数个数,因此还提出了对同一介质物性单元的Jacobi矩阵元素进行压缩组装的措施,从而使反演的未知量个数减少,结合采用共轭梯度迭代法,使得只需利用有效波频段的少数一些频率即可进行迭代反演.正演和反演理论模型的数值模拟结果表明方法是有效的. 相似文献
5.
提出了非线性共轭梯度法大地电磁拟三维反演.该方法选取共轭梯度反演算法为拟三维反演的核心.在计算灵敏度(Jacobian)矩阵时,吸取近似灵敏度矩阵思想,采用一维灵敏度矩阵来代替三维灵敏度矩阵,并对非测点的灵敏度元素提出一种近似方法.在第一次反演之后,采用拟牛顿法更新灵敏度矩阵.拟三维反演法在很大程度上节省了计算时间,并且理论模型和实际资料的反演试算结果表明大地电磁拟三维反演法具有一定的实用价值. 相似文献
6.
推导出频率域有限元声波正演方程,为了消除边界反射,将Clayton-Engquist旁轴波动方程吸收边界条件引入频率域,并对有限元刚度矩阵和质量矩阵进行压缩存储,利用广义共轭梯度法求解有限元方程获得正演解.在此基础上,推导出在某一频率下波场数据残差δU与单元物性参数修改量δλ之间关系的Jacobi矩阵,反演方法允许利用地面二维炮集全波场资料与给出初始模型参数的正演值的差值δU,迭代求得δλ.由于计算机内存的限制,方法计算不允许有过多数目的未知数个数,因此还提出了对同一介质物性单元的Jacobi矩阵元素进行压缩组装的措施,从而使反演的未知量个数减少,结合采用共轭梯度迭代法,使得只需利用有效波频段的少数一些频率即可进行迭代反演.正演和反演理论模型的数值模拟结果表明方法是有效的. 相似文献
7.
应用地震资料约束下大地电磁资料反演古潜山或基底内部结构的新方法. 首先利用地震及井资料,反演得到古潜山或基底之上各层的厚度及电阻率,建立了古潜山或基底以上的地电模型;并利用大地电磁一维反演资料,建立了古潜山或基底及其以下的地电模型. 在初始模型建立的基础上,应用高效快速的二维大地电磁正演算法——模式匹配方法,进行正演计算. 利用共扼梯度迭代方法,进行正反演迭代计算. 通过对胜利油田南北618剖面反演, 基底以下在地震资料没有显示的构造信息,在反演结果中得到了较好地反映. 相似文献
8.
本文用广义逆矩阵理论讨论了时间域航空电磁数据的反演方法.以水平二层大地和球体模型为例进行了计算,结果表明这种方法对时间域电磁数据的反演是行之有效的.对于理论数据,一般仅需迭代数次即可稳定地收敛到真值.通过奇异值分析,可以确定一个模型中的相对重要参数和无关紧要参数;对于求解某一特定参数,可确定哪些数据起决定性的作用.由奇异值分解而得出的信息密度矩阵、模型分辨矩阵,分别给出了模型响应拟合数据和模型参数分辨率的度量.这些信息对于时间域航空电磁系统的设计、野外测量和资料解释都具有重要的参考价值.最后给出了实例. 相似文献
9.
本文用广义逆矩阵理论讨论了时间域航空电磁数据的反演方法.以水平二层大地和球体模型为例进行了计算,结果表明这种方法对时间域电磁数据的反演是行之有效的.对于理论数据,一般仅需迭代数次即可稳定地收敛到真值.通过奇异值分析,可以确定一个模型中的相对重要参数和无关紧要参数;对于求解某一特定参数,可确定哪些数据起决定性的作用.由奇异值分解而得出的信息密度矩阵、模型分辨矩阵,分别给出了模型响应拟合数据和模型参数分辨率的度量.这些信息对于时间域航空电磁系统的设计、野外测量和资料解释都具有重要的参考价值.最后给出了实例. 相似文献
10.
三维反演解释是电磁法勘探发展的重要趋势,而如何提高三维反演的可靠性、稳定性和计算效率是算法开发者们目前的研究重点.本文实现了一种频率域可控源电磁(CSEM)三维反演算法.其中正演基于拟态有限体积法离散化,利用直接矩阵分解技术来求解大型线性系统方程,不仅准确、稳定,而且特别有利于含有大量发射场源位置的CSEM勘探情况;对目标函数的最优化采用高斯牛顿法(GN),具有近似二次的收敛性;使用预条件共轭梯度法(PCG)求解每次GN迭代所得到的法方程,避免了显式求解和存储灵敏度矩阵,减小了计算量.以上这些方法的结合应用,使得本文的三维反演算法准确、稳定且高效.通过陆地和海洋CSEM勘探场景中的典型理论模型的反演测试,验证了本文算法的有效性. 相似文献
11.
常规的三维时间域航空电磁模拟通常采用隐式步长方法进行时间离散,需要几次矩阵分解和上百次右端源项回带,计算效率较低.为了提高正演计算效率,本文提出使用有理Krylov方法求解时间域电场扩散方程.首先使用非结构四面体网格进行空间离散,采用Nédélec矢量基函数近似四面体单元内的电场;然后基于有限元离散给出矩阵指数和矢量乘积表示的电场显式解;最后采用有理Arnoldi算法构造Krylov子空间内的正交基函数并进一步求解矩阵指数与矢量的乘积,直接得到任意时刻的电场解向量,避免步长离散过程.此外,本文还提出一种指数加权偏移参数优化方法,使得有理Arnoldi近似在瞬变衰减晚期具备更高的精度,从而降低Krylov子空间阶数并提高计算效率.通过和层状模型解析解的对比验证了有理Krylov方法的精度.针对三维异常体模型使用全局网格和局部网格剖分并和其他数值方法比较,进一步说明了有理Krylov方法的有效性. 相似文献
12.
位场向下延拓是位场数据处理和反演中的重要运算,但是它的不稳定性影响了它在许多处理和反演方法技术中的应用.本文通过把位场向下延拓视为向上延拓的反问题,得到向下延拓的褶积型线性积分方程,再利用Fourier变换矩阵的正交对称特性,并结合矩阵的奇异值分解和广义逆原理,提出了一种稳定的不需要进行求逆运算的位场向下延拓广义逆方法——波数域广义逆算法,解决了位场大深度向下延拓的不稳定性问题.把这种方法用于三维理论模型数据和实际磁场数据的向下延拓获得了理想的结果. 相似文献
13.
To improve the inversion accuracy of time-domain airborne electromagnetic data, we propose a parallel 3D inversion algorithm for airborne EM data based on the direct Gauss–Newton optimization. Forward modeling is performed in the frequency domain based on the scattered secondary electrical field. Then, the inverse Fourier transform and convolution of the transmitting waveform are used to calculate the EM responses and the sensitivity matrix in the time domain for arbitrary transmitting waves. To optimize the computational time and memory requirements, we use the EM “footprint” concept to reduce the model size and obtain the sparse sensitivity matrix. To improve the 3D inversion, we use the OpenMP library and parallel computing. We test the proposed 3D parallel inversion code using two synthetic datasets and a field dataset. The time-domain airborne EM inversion results suggest that the proposed algorithm is effective, efficient, and practical. 相似文献
14.
Although waveform inversion has been intensively studied in an effort to properly delineate the Earth's structures since the early 1980s, most of the time‐ and frequency‐domain waveform inversion algorithms still have critical limitations in their applications to field data. This may be attributed to the highly non‐linear objective function and the unreliable low‐frequency components. To overcome the weaknesses of conventional waveform inversion algorithms, the acoustic Laplace‐domain waveform inversion has been proposed. The Laplace‐domain waveform inversion has been known to provide a long‐wavelength velocity model even for field data, which may be because it employs the zero‐frequency component of the damped wavefield and a well‐behaved logarithmic objective function. However, its applications have been confined to 2D acoustic media. We extend the Laplace‐domain waveform inversion algorithm to a 2D acoustic‐elastic coupled medium, which is encountered in marine exploration environments. In 2D acoustic‐elastic coupled media, the Laplace‐domain pressures behave differently from those of 2D acoustic media, although the overall features are similar to each other. The main differences are that the pressure wavefields for acoustic‐elastic coupled media show negative values even for simple geological structures unlike in acoustic media, when the Laplace damping constant is small and the water depth is shallow. The negative values may result from more complicated wave propagation in elastic media and at fluid‐solid interfaces. Our Laplace‐domain waveform inversion algorithm is also based on the finite‐element method and logarithmic wavefields. To compute gradient direction, we apply the back‐propagation technique. Under the assumption that density is fixed, P‐ and S‐wave velocity models are inverted from the pressure data. We applied our inversion algorithm to the SEG/EAGE salt model and the numerical results showed that the Laplace‐domain waveform inversion successfully recovers the long‐wavelength structures of the P‐ and S‐wave velocity models from the noise‐free data. The models inverted by the Laplace‐domain waveform inversion were able to be successfully used as initial models in the subsequent frequency‐domain waveform inversion, which is performed to describe the short‐wavelength structures of the true models. 相似文献
15.
The dynamic response of a double infinite beam system connected by a viscoelastic foundation under the harmonic line load is studied. The double infinite beam system consists of two identical and parallel beams, and the two beams are infinite elastic homogeneous and isotropic. A viscoelastic layer connects the two beams continuously. To decouple the two coupled equations governing the response of the double infinite beam system, a variable substitution method is introduced. The frequency domain solutions of the decoupled equations are obtained by using Fourier transforms as well as Laplace transforms successively. The time domain solution in the generalized integral form are then obtained by employing the corresponding inverse transforms, i.e. Fourier transform and inverse Laplace transform. The solution is verified by numerical examples, and the effects of parameters on the response are also investigated. 相似文献
16.
In this paper we propose a 3D acoustic full waveform inversion algorithm in the Laplace domain. The partial differential equation for the 3D acoustic wave equation in the Laplace domain is reformulated as a linear system of algebraic equations using the finite element method and the resulting linear system is solved by a preconditioned conjugate gradient method. The numerical solutions obtained by our modelling algorithm are verified through a comparison with the corresponding analytical solutions and the appropriate dispersion analysis. In the Laplace‐domain waveform inversion, the logarithm of the Laplace transformed wavefields mainly contains long‐wavelength information about the underlying velocity model. As a result, the algorithm smoothes a small‐scale structure but roughly identifies large‐scale features within a certain depth determined by the range of offsets and Laplace damping constants employed. Our algorithm thus provides a useful complementary process to time‐ or frequency‐domain waveform inversion, which cannot recover a large‐scale structure when low‐frequency signals are weak or absent. The algorithm is demonstrated on a synthetic example: the SEG/EAGE 3D salt‐dome model. The numerical test is limited to a Laplace‐domain synthetic data set for the inversion. In order to verify the usefulness of the inverted velocity model, we perform the 3D reverse time migration. The migration results show that our inversion results can be used as an initial model for the subsequent high‐resolution waveform inversion. Further studies are needed to perform the inversion using time‐domain synthetic data with noise or real data, thereby investigating robustness to noise. 相似文献
17.
时间域航空电磁作为一种高效地球物理勘探技术特别适合我国地形复杂地区(沙漠、高山、湖泊、沼泽等)资源勘查.然而,这些地区地形起伏较大,对航空电磁响应有严重影响,忽略地形影响会给航空电磁数据解释造成很大误差.到目前为止人们对航空电磁地形效应特征研究十分有限.本文提出了基于非结构化网格的有限元法模拟带地形时间域航空电磁系统响应.该方法与基于结构化网格的有限差分相比能更好地模拟地形.首先通过傅里叶变换将2.5维问题转化成二维问题,利用伽辽金方法对二维问题进行离散.通过使用MUMPS求解器,得到波数域电磁响应.利用反傅里叶变换将波数域电磁响应变换到空间域,并利用正弦变换将其变换到时间域,得到2.5维时间域航空电磁响应.通过将本文的计算结果与半空间模型解析解及其他已发表的结果进行对比,检验了本文算法的精度.最后,我们系统分析了山峰和山谷地形对航空响应的影响特征.本文研究结果对航空电磁地形效应的识别和校正具有指导意义. 相似文献
18.
目前,瞬变电磁法(TEM)数据基本都是基于各向同性模型进行反演解释,这对于存在明显电性各向异性的勘探区域会产生较大的反演解释误差.为分析电各向异性对回线源瞬变电磁信号的影响方式与程度,本文通过求解离散化的全张量电导率时间域Helmholtz方程,实现了基于有限体积法的TEM任意各向异性的三维正演算法.该算法采用基于交错网格的拟态有限体积法(MFV)对时域Maxwell方程组进行空间域离散,并利用后退欧拉算法(Backward Euler Method)进行时间域离散.为提高时域电磁场的求解精度与效率,该算法将时间分段等步长算法与方程直接求解法相结合.通过对一维各向异性模型以及三维复杂各向同性模型进行测试,验证了本算法对于回线源瞬变电磁响应计算的正确性及有效性.最后,通过对几类典型电各向异性介质中大回线源瞬变电磁信号响应的分析,总结了不同电各向异性类型对TEM电磁信号的影响模式,结果表明,主轴各向异性情况下TEM信号主要受水平方向电导率的影响,倾斜各向异性对TEM信号的影响程度远大于水平各向异性,而通过水平各向异性信号能较清晰判断出各向异性主轴方向. 相似文献
19.
以二维声波方程为模型,在时间域深入研究了全波形速度反演.全波形反演要解一个非线性的最小二乘问题,是一个极小化模拟数据与已知数据之间残量的过程.针对全波形反演易陷入局部极值的困难,本文提出了基于不同尺度的频率数据的"逐级反演"策略,即先基于低频尺度的波场信息进行反演,得出一个合理的初始模型,然后再利用其他不同尺度频率的波场进行反演,并且用前一尺度的迭代反演结果作为下一尺度反演的初始模型,这样逐级进行反演.文中详细阐述和推导了理论方法及公式,包括有限差分正演模拟、速度模型修正、梯度计算和算法描述,并以Marmousi复杂构造模型为例,进行了MPI并行全波形反演数值计算,得到了较好的反演结果,验证了方法的有效性和稳健性. 相似文献
20.
利用钻孔进行超前探测地质构造及含水体是地下开挖工程中的常规手段,如何利用这些钻孔进行钻孔瞬变电磁法扫描探测,从而实现钻孔孔壁外围地质异常体的精细探测,对实现地下工程地质透明化具有重要的指导意义.本文提出钻孔瞬变电磁法扫描探测2.5D反演的数据解译方法,首先针对随机性反演算法时效性低,易陷入局部最优解,而确定性反演算法依赖初始模型的问题,提出了组合策略的量子粒子群优化算法用来随机搜索最优初始模型.在此基础上,利用Levenberg-Marquarat方法求解Occam反演的目标函数,形成了RCQPSO-LMO组合算法进行2.5D反演,通过对比组合算法和单一算法,验证了组合算法具有更精确的反演结果.其次结合屏蔽条件下扫描探测,对比分析了有无屏蔽的2.5D反演结果,通过设定屏蔽系数对非探测方向信号进行部分压制,可以较好地解决钻孔径向扫描探测中对非探测方向信号部分屏蔽下的反演及成像.最后建立三组理论模型进行组合算法2.5D反演,结果表明:组合算法反演结果与理论模型的一致性较好,对低阻异常体的反演精度较高,验证了组合算法对钻孔孔壁外围低阻异常体具有较高的反演精度和分辨能力. 相似文献
|