共查询到20条相似文献，搜索用时 703 毫秒
1.
起伏地形频率域可控源电磁二维快速正反演 总被引：1，自引：0，他引：1
起伏地形频率域可控源电磁正反演研究受到众多学者的关注和重视,目前普遍实用的为未考虑地形的二维正反演算法程序。本文基于积分方程法正演和多场源、多频率、对比源反演算法,开展起伏地形条件下频率域可控源电磁二维快速正反演研究。通过引入层状参考模型、将地形与目标体整体作为异常场剖分区域,实现起伏地形可控源电磁正反演计算。针对包含地形在内的大尺度剖分区域引起的大型计算代价问题,采用快速傅里叶算法提高正反演计算效率。通过与现有正反演算法进行模型算例对比,验证了本文采用方法的可行性与有效性。 相似文献
2.
三维磁异常傅里叶正演,能够计算整个场源区域内与磁化率分布网格一致的三维磁异常。通过三维傅里叶变换推导了长方体三维磁谱表达式,当源体埋藏较深或者位于反演区域边缘时,标准FFT正演异常场由于强制周期化,边界震荡等原因,使得正演结果发生较大的畸变,为了减少标准FFT算法引起的误差,引入了3DGauss-FFT技术用于三维磁异常频率域正演。通过简单的模型正演验证,从计算时间、计算精度以及内存需求上与空间域算法及标准3DFFT算法进行比较,结果表明,3DGauss-FFT磁异常正演相比于标准FFT算法在计算精度上有很大提升,显著降低了标准3DFFT由于自身缺陷引起的误差,且在运行时间上,4点3DGauss-FFT磁异常正演算法相对于空间域算法降低了三个数量级,但内存需求有所增加,证明了3DGauss-FFT正演方法在磁异常正演方面的高效性以及准确性。 相似文献
3.
《物探化探计算技术》2017,(1)
地震波场频率域正演是频率域全波形反演的基础。针对反演计算量巨大的问题,利用频率域二维波形正演算法中频率和炮点计算的独立性,开发出粗细粒度结合的MPI并行算法:粗粒度为频率并行,细粒度为单个频率解方程并行。实现方法是将正演频率分组放入不同的MPI通讯域内,每个通讯域内单个频率求解方程过程,采用基于MPI的MUMPS(多波前大规模并行稀疏直接解法器)软件包并行加速。模型测试结果表明:MPI多通讯域并行算法计算结果正确,计算效率显著提高,加速效果稳定。 相似文献
4.
《物探与化探》2021,(3)
Z轴倾子电磁法(ZTEM)是一种新型频率域航空天然场源电磁法,该方法采用倾子作为研究参数,可用于起伏地形下的大规模地质勘探。本文在对TEM三维有限差分正演和数据空间OCCAM反演算法研究的基础上,考虑地形起伏的影响,研究了带地形的频率域三维ZTEM正反演算法。首先,对起伏地表下ZTEM正演算法的准确性进行了验证,计算和分析了纯地形的ZTEM三维异常响应特征。其次,通过对山峰和山谷地形下低阻棱柱体合成算例的反演分析,表明带地形的ZTEM数据空间OCCAM反演算法能够获得比较接近真实导电性结构的地电模型,尤其对地下目标体的横向边界具有理想的约束效果。最后,与不带地形的ZTEM反演结果对比,验证了所开发的带地形ZTEM倾子资料反演算法的有效性。 相似文献
5.
随着传感器技术的发展,重力场与重力张量场测量技术发展迅速,为实现地下密度分布精细反演提供了数据保障。正演是反演的基础,解决任意密度分布复杂地质体重力场与重力张量正演高效、高精度计算问题,是实现重力高效、精细反演、人机交互反演解释的关键。针对起伏地形和任意密度分布这种复杂条件下二维重力场及重力张量场高效高精度正演问题,这里提出了一种空间波数混合域正演算法,其关键环节包括:①结合新的矩形二度体组合模型波数域表达式和一维Gauss-FFT算法,提出了一种任意密度分布和起伏地形下重力场及重力张量高效、高精度正演算法;②采用新的二维正演算法,计算观测最高点和最低点之间多个不同高度水平网格重力场及重力张量,结合三次样条插值方法,实现了起伏地形上重力场及重力张量场高效、高精度正演。模型算例结果表明,新方法具有高效、高精度的显著特点。 相似文献
6.
在地震勘探中,地震正演模拟是非常重要的技术。与时间域正演相比,频率域正演速度快,计算效率高。如何高效准确地完成频率域正演计算是目前该领域的一个重要问题。数值频散问题和如何提高计算效率降低求解分解阻抗内存占用量一直是频率域正演所需要解决的问题。与传统的直接法求解阻抗矩阵的频率域正演方法不同,本文采用最小二乘共轭梯度法(LSCG法)求解阻抗矩阵进行频率域正演,并提出了一种波数补偿的表达式来压制数值频散现象。经过简单模型和复杂模型的数值测试,采用最小二乘共轭梯度法(LSCG法)求解阻抗矩阵进行频率域正演能够有效降低计算时间,且采用波数补偿的频率域正演方法能够有效压制数值频散现象,提高波场模拟精度。 相似文献
7.
8.
《物探与化探》2017,(2)
为优化二维各向同性介质中弹性波频率域正演时阻抗矩阵的结构,减小正演所需内存,提高正演效率,在25点差分格式的基础上进行适当的简化,得到了二维弹性波频率域17点差分格式。该格式重新计算了弹性波中偏微分项和加速项的差分算子,减少了计算过程中的网格节点需求,构造了优化阻抗矩阵后的频率域正演矩阵方程;推导了纵波和横波相速度的频散公式,给出了不同泊松比条件下的频散曲线,得到了相速度误差控制范围±1%时每一横波波长内网格数需求。通过对比频散曲线和数值模拟时得到的波场快照及检波点处U、V分量,验证了17点差分格式与25点差分格式相比,具有稍严格的网格间距需求、相当的计算精度、略少的计算时间和更小的阻抗矩阵带宽等特点。 相似文献
9.
频率域正演是频率域全波形反演的基础,高效准确地完成频率域正演计算是目前该领域的一个热门问题。这里利用优化17点差分格式和PML完全匹配层构建波场,并使用LSQR算法求解波场,通过在简单层状模型和Overthrust模型上计算验证可知,LSQR算法可大幅提高频率域正演的计算速度和模拟精度,从而为此方向的研究提供一个可行的技术手段。 相似文献
10.
穆石敏 《物探化探计算技术》1981,(4)
正反演是重磁解释中很重要的一个步骤,过去这部分工作都是在空间域中进行的,如各种规则形体的正演公式,点元线元或面元的正演程序,最优化选择法(如马奎特方法)反演程序等,这些在实际工作已取得一定效果。近一年来,频率域位场转换工作已在数据处理中广为利用,同时,频率域正反演工作也在逐步开展,如组合斜顶面棱柱体正演程序,功率谱计算基底深度程序等已在开始使用,此外对频率域其他反演方法也在研究与试用。考虑 相似文献
11.
全波形反演是一种全新的地震成像方法,主要利用全波形信息反演地下介质参数,通过非线性优化波场理论值和观测值的残差实现波形反演。基于时间域声波方程,建立了时间域波场残差目标函数,以分层模型为例,分别从波场精度、目标函数收敛性和运行时间3个方面,比较了共轭梯度（Conjugate Gradient,CG）算法和拟牛顿算法（Broyden-Fletcher-Goldfarb-Shanno,BFGS）反演的效果。同时,应用共轭梯度法对正、逆断层模型和Marmousi模型进行了速度结构反演。反演结果表明:共轭梯度法计算效率较高,反演得到的速度模型精度更高,反演效果较好,是一种有效的波形反演方法。 相似文献
12.
研究了二阶声波方程频域PML边界条件和频域变网格并行计算技术。PML边界是一种较为理想的吸收边界方法,多用在求解时域应力速度方程中,但对于频域声波正演,二阶位移方程更常用。从一阶声波方程PML吸收边界条件导出频域二阶位移方程PML边界条件,模拟算例得到的频率切片、时间切片和地震记录对比都说明该边界条件吸收效果很好。频域单炮正演不同频率间是独立的,据此低频部分采用大网格计算,高频采用小网格,实现变网格步长计算技术,这是较时间域正演的一个优势,在保证模拟质量的同时,减少计算量和内存消耗。 相似文献
13.
斜缆采集资料中鬼波的陷频特征限制了地震记录的频带宽度和分辨率,给地震资料反演、解释带来困难。通过压制鬼波,可获得高分辨率的宽频数据。基于拉东域的水平拖缆鬼波压制方法,结合斜缆鬼波随偏移距变化的特点,推导出一次波和鬼波在频率-拉东域的逆变换算子,建立检波器处总波场和海表面处上行波场之间新的关系式,利用最小二乘反演精确求解获得海表面处上行波场,并延拓得到拖缆处鬼波压制后的记录。通过考虑鬼波延迟时间受偏移距和出射角的影响,弥补鬼波延迟时间估计中存在的误差,无需进行反演迭代求取最优鬼波延迟时间,提高了计算效率。合成数据及海上实际斜缆数据测试结果表明,研究方法能较好地压制鬼波,达到拓宽地震记录频带的目的。 相似文献
14.
全波形反演方法利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而,巨大的计算量是阻碍其发展的一个瓶颈问题。为此,研究者们提出了震源编码技术来减少计算量,但是此方法在模型更新过程中会引进随机串扰噪声,降低反演结果准确性。所以,在保证计算精度的情况下,本文提出了采用逐减随机震源采样的方法来高效计算全波形反演问题。笔者将此方法应用于频率域二维黏滞声波波动方程全波形反演,开始了在频率域进行随机震源采样类方法的研究,计算过程中共使用了依次增大的8个频率段;并应用Overthrust模型来验证此类随机震源采样法的正确性。实验结果表明：基于逐减随机震源采样法的反演结果与实际Overthrust模型的拟合误差为0.065 65,而应用基于全部震源的全波形反演方法得到的反演结果与实际Overthrust模型的拟合误差为0.064 64,两者差别不大;但计算用时由740 min减少到291.2 min,即计算效率提高了2.54倍。为了更好地确定方法的有效性,将其应用于Marmousi模型进行试算。模型试算结果表明：基于逐减随机震源和基于全部震源得到的反演结果与实际Marmousi模型的拟合误差分别为0.080 12和0.078 97,相差不大;但计算用时由1 218.9 min减少到274.4 min,计算效率提高了4.44倍。综上,在保证反演精度的情况下,基于逐减随机震源采样法的频率域全波形反演方法大大减少了计算量,具有不可替代的计算优势,并且没有引进随机串扰噪声。 相似文献
15.
16.
全波形反演可以利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而,庞大的计算量和存储空间需求,限制了全波形反演的发展。在频率多尺度全波形反演中将L BFGS数值优化算法与同时激发震源技术相结合的方法来改善这一现状。首先,对Marmousi模型进行了速度反演：在计算过程中明显发现对计算机内存的占用减少,最终反演结果与实际Marmousi模型的拟合误差为0.095 9,较小;采用10个频带单炮震源正演384炮所需时间约为32 640 s,而采用同时激发震源（384炮）正演一次所需时间仅约为700 s。然后,基于高速楔形体模型进行了抗噪能力研究：原始含噪地震记录信噪比为11.147 3 dB;对反演得到的速度模型进行正演,其地震记录信噪比为22.251 8 dB。最后,基于逆冲断层模型进行了反演速度扰动能力研究,反演得到的最终模型很清晰,与具有速度扰动特性的实际模型非常接近,拟合误差仅为0.036 0。数值模拟试验结果表明：此方法反演精度高,内存开销较小,能够显著提高计算效率,并且具有良好的抗噪能力,能够反演出具有速度扰动特性的介质。 相似文献
17.
构建近地表横波速度模型是煤田多分量地震资料处理的重要环节。相较于面波多道分析法,全波形反演在构建近地表横波模型中具有更高的分辨率。然而,在基于梯度的全波形反演中,由于地震记录频带有限、波场的非均匀覆盖以及双重散射等原因导致梯度算子不随深度的增加而缩放,模型深部参数得不到明显更新。目标函数的Hessian算子包含曲率信息,可清晰预测梯度算子中的焦散现象及双重散射产生的伪影,因此,逆Hessian算子则可作为反卷积算子实现对梯度的预处理,加强对模型深部的照明能力。然而Hessian算子具有巨大维度,对其显式计算十分困难。基于此,借鉴逆散射理论的思想,给出勒夫波全波形反演目标函数的拟Hessian算子的表达式,并提出一种梯度预处理的全波形反演方法。将该方法分别应用于断层模型、凹陷模型以及起伏界面模型的重构试验,反演结果表明:与传统的共轭梯度全波形反演方法相比,基于拟Hessian算子的预处理共轭梯度方法可加快收敛速度,提升成像质量。 相似文献
18.
19.
利用高阶交错网格有限差分模拟Kelvin-Voigt黏弹性介质中传播的地震波,同时将完全匹配层吸收边界条件引入到其边界处理中。数值模拟结果表明,完全匹配层吸收边界效果好,高阶有限差分能模拟得到的黏弹性介质波场精度较高。对模拟的黏弹性波场进行分析,表明介质的粘滞性使地震反射波的能量变弱,高频衰减明显,并比低频衰减得快,主频向低频方向移动,有效频带变窄,即降低了地震波的分辨率;并且反射转换波比反射纵波要衰减得快;而且还随着传播距离的增加,其峰值频率也逐渐降低。通过数值模拟分析具有不同的粘滞系数介质对地震波的吸收和衰减,结果表明随着粘滞系数的增大,地下介质对地震波的吸收衰减更明显。 相似文献
20.
THE HIGH RESOLUTION SEISMIC TOMOGRAPHIC IMAGE IN QINGHAI—TIBET PLATEAU AND ITS DYNAMIC IMPLICATIONS 总被引：1，自引：1，他引：0
THE HIGH RESOLUTION SEISMIC TOMOGRAPHIC IMAGE IN QINGHAI—TIBET PLATEAU AND ITS DYNAMIC IMPLICATIONSeasthenospherehadbe 相似文献