共查询到20条相似文献，搜索用时 765 毫秒
1.
抛物Radon变换法(Parabolic Radon Transform)在地震资料处理中有广泛的应用。PRT可对不同频率的地震数据解耦处理,这一特点使得抛物Radon变换的计算效率比双曲Radon变换有数量级上的提高。在频率域求解时,需要对每一个频率成份求解同样大小的线性方程组。求解抛物Radon正变换的计算方法主要有Levinson递推法、共轭梯度法、Cholesky分解法和直接矩阵求逆法。最小平方抛物Radon正变换所形成的矩阵具有Toeplitz结构,可采用Levinson递推法进行计算。高分辨率抛物Radon正变换所形成矩阵的Toeplitz结构被破坏,一般采用共轭梯度法或Cholesky分解法进行求解。这里详细推导了复Toeplitz矩阵的Levinson递推算法,并分别对求解方程的四种方法进行了讨论,最后给出抛物Radon正变换求解的数值算例,并对所给出的四种方程求解方法的计算效率及计算精度进行了对比。 相似文献
2.
《物探与化探》2017,(2)
为优化二维各向同性介质中弹性波频率域正演时阻抗矩阵的结构,减小正演所需内存,提高正演效率,在25点差分格式的基础上进行适当的简化,得到了二维弹性波频率域17点差分格式。该格式重新计算了弹性波中偏微分项和加速项的差分算子,减少了计算过程中的网格节点需求,构造了优化阻抗矩阵后的频率域正演矩阵方程;推导了纵波和横波相速度的频散公式,给出了不同泊松比条件下的频散曲线,得到了相速度误差控制范围±1%时每一横波波长内网格数需求。通过对比频散曲线和数值模拟时得到的波场快照及检波点处U、V分量,验证了17点差分格式与25点差分格式相比,具有稍严格的网格间距需求、相当的计算精度、略少的计算时间和更小的阻抗矩阵带宽等特点。 相似文献
3.
实现稀疏反褶积的预条件双共轭梯度法 总被引：9，自引：3，他引：9
地震勘探稀疏反褶积计算一般要导出一个Toeplitz矩阵的线性系统，通常可以用矩阵求逆、Levison递推及共轭梯度等方法直接求解。当Toeplitz矩阵的条件数很大时，数值稳定性差，甚至无法求解。使用共轭梯度法，在矩阵的对角元素上加入规则化因子，可以改善这种情况，但不能彻底解决数值稳定性和精度问题。若求解最小二乘问题的原始问题，结果会好些。线性系统形式的细微改变，将导致不同的数值计算特性。在规则化策略基础上，可巧妙地构造稀疏反褶积的问题原型，引入预条件，采用双共轭梯度法求解，从而实现稀疏反褶积，获得较好结果。数值算例表明，预条件双共轭梯度法比直接稀疏反褶积方法收敛快、精度高。 相似文献
4.
5.
6.
Matlab环境下瑞利波有限差分正演与曲线绘制 总被引：1，自引：0，他引：1
以瑞利波频散方程为出发点的Abo-Zena传递矩阵等方法,只能研究层状介质中瑞利波的传播特性,对于非层状介质,传递矩阵无能为力。因此,不得不考虑有限元、有限差分等方法。根据弹性动力学方程,采用交错网格有限差分方法对均匀弹性半空间介质进行全波场模拟,并在Matlab环境下实现编程计算,再现了瑞利波在近地表的传播状态,实现了地震剖面曲线的正振幅充填问题,从而获得了类似野外地震勘探的剖面记录。模拟结果表明,采用一阶差分格式所得到的地震记录存在较严重的数值频散,在采用有限差分法研究瑞利波"之"字形成因时,须尽量压制数值频散现象,才能获得更接近实际情况的频散曲线。 相似文献
7.
殷文 《吉林大学学报(地球科学版)》2008,38(1):144-151
在弹性波频率空间域有限差分数值模拟方面,差分网格及边界条件是影响弹性波模拟成功与否的关键,为了压制数值模拟中的网格频散,采用25点有限差分算子,建立了有限差分矩阵方程,且借鉴匹配层衰减边界条件思想,设计了弹性波频率空间域有限差分数值模拟算法。由于采用高阶有限差分法来提高差分格式的精度,将会导致计算量显著增加,为此,对频率空间域有限差分弹性波数值模拟方法,采用流水线技术与分治策略进行了并行算法研究,提高了计算效率,使得在合理的计算时间内更精确地模拟弹性波在弹性介质中的传播过程。 相似文献
8.
9.
10.
三维电阻率法对反演的精度和速度的要求越来越高,而正演是反演的基础,因此直流电阻率三维正演计算的速度和精度是三维电阻率反演实用化的关键。这里利用对称超松弛预条件共轭梯度法(SSOR-PCG),求解有限差分法离散生成的大型稀疏线性方程组,预条件矩阵的选择大大降低了系数矩阵的条件数,结合矩阵的一维非零元素压缩存储模式,使得正演计算速度得以提高,而内存占用量明显减小。在直流电阻率三维正演中采用异常场法,提高了电源点附近的解的精度。利用编制的有限差分正演程序,对两层模型、垂直接触带模型和低阻异常体模型进行了数值模拟,计算结果表明该算法是可行的,且可以明显提高正演计算的速度和精度。 相似文献
11.
利用小波变换快速反演大型磁数据 总被引：4，自引：0，他引：4
我们应用小波变换反演大型磁数据,恢复三维磁化率分布问题,快速小波变换通过给出小波系数的门槛值来产生灵敏度矩阵表达式,在波域进行正演模拟计算,通过小波域内稀疏矩阵矢量的乘积运算,和压缩后的矩阵来实施快速正演模拟,所用ＣＰＵ的时间的减少与矩阵压缩在正比,而且,我们采用对数阻尼法,用共轭梯度解算器来获得反演结果,这种求解手段的组合,直接减少了正演模拟的时间,进而加快了反演速度,产生了更的三维磁反演算法。 相似文献
12.
13.
面波频散能量谱成像精度是近地表横波速度计算的重要影响因素,有效提取面波频散曲线是后续反演地层结构及介质物性参数的重要环节。文中论述了τ-p变换法（τ为垂直波慢度,p为水平波慢度）、高分辨率线性Radon变换法、频率分解法、相移法、f-k变换法（f为频率,k为波数）、SFK变换法和矢量波数变换法等7种不同面波频散能量谱计算方法的基本原理;通过上述方法,对三层速度递增水平层状介质模型正演模拟合成Rayleigh波记录进行频散能量谱计算,并与理论频散曲线进行比较;根据所得频散能量谱,分析不同方法在数据处理过程中的特点与应用差异。结果表明：在提取频散曲线过程中,高分辨率线性Radon变换法和矢量波数变换法提取精度较高,f-k变换法提取精度最低,低频信息与高频信息无法有效识别;SFK变换法、矢量波数变换法、高分辨率线性Radon变换法抗噪性能强;τ-p变换法、f-k变换法、相移法计算速度优于其他方法。 相似文献
14.
波动方程数值模拟是研究地震波传播机理的重要工具,有限差分求解波动方程是当前地震波数值模拟的主要方法之一。当地下介质中的地震波速度较低或地震波高频成分丰富时,常规有限差分技术常常产生严重的数值频散误差,这种误差会降低数值模拟的精度,影响对地震波传播机理的分析。为压制地震波数值模拟时产生的数值频散误差,提高波场模拟精度,提出了基于NAD算子的时间四阶精度波动方程差分格式。根据对应的差分格式,分析了该差分格式的数值频散关系。与常规四阶精度差分算法的频散曲线相比,基于NAD时间四阶精度差分方法不但能够实现时间频散的有效压制,同时其基于更多网格点的位移分量和位移梯度分量空间微分求解方法还能够实现空间频散的有效压制。另外在相同模型条件下,基于NAD算法的声波方程时间四阶差分解法可采用大网格对模拟空间进行差分离散,减少网格数,提高计算效率。 相似文献
15.
抛物线Radon变换是地震资料处理中常用的一种多次波压制方法,常规Radon变换由于采用最小二乘算法的原因,对不同的曲率值采用相同的加权值进行求解,分辨率较低。为提高Radon域数据的分辨率,应采用对不同曲率值应用不同的加权值的方法。一般的思路是利用前一次迭代的结果得到加权矩阵,该算法需要迭代进行,运算量较大。这里拟采用低频约束的方式提高Radon变换的分辨率,即利用前一个频率的运算结果对下一个频率的计算进行约束,该算法分辨率高,计算速度快。模拟数据和实际数据的测试表明,本文方法在多次波压制处理中,可以快速有效地完成多次波去除。 相似文献
16.
研制了利用共轭梯度松驰法的三维直流电阻率正演模拟与反演算法。在正演网络的模拟计算中，将一种用于预处理的不完全Ｃｈｏｌｅｓｋｙ分解法与稀疏矩阵程序结合起来，便可得到一种快速有效的算法。侧面和底部的边界条件用阻抗条件来标度，而阻抗条件考虑了由于电源流的一些配置而在边界上产生的电流。 相似文献
17.
TEM正演响应计算的几种频-时域转换方法对比 总被引：1，自引：0，他引：1
在瞬变电磁的一维正演响应模拟中,常用的方法是先在频率域中求解,之后将结果转换到时间域。但该方法在晚期的计算精度通常不高,因此,使用5种频-时域转换方法(正、余弦变换的数值滤波算法,G-S逆拉普拉斯算法、正、余弦变换的折线逼近法)进行了计算,与解析解对比,得出余弦变换的数值滤波算法在晚期计算中精度最高的结论,并对这五种转换方法产生误差问题的原因进行了讨论和分析。本研究有利于瞬变电磁一维正演响应的高精度计算,使其在多维计算中得到更好的应用。 相似文献
18.
TTI介质弹性波频率-空间域有限差分数值模拟 总被引：3，自引：0，他引：3
由周期性薄互层引起的VTI介质是研究比较广泛的一类各向异性介质。当VTI介质对称轴偏离垂向,本构坐标系与观测坐标系不重合时,会形成观测坐标系下的TTI介质。引入25点优化差分算子,推导出二维TTI介质频率域弹性波动方程;为压制边界反射,采用完全匹配层法吸收边界条件,并计算出优化差分系数;最后采用集中力源,模拟了弹性波在TTI介质中的传播过程。从波场快照和地面共炮记录可以看出,笔者采用的数值模拟算法能有效压制数值频散。TTI介质中的波场传播比较复杂,纵波传播相对稳定,横波波前的三分叉现象比较明显,并存在振幅奇异性。当VTI介质的对称轴偏转后,还会增加地面地震记录的复杂性。 相似文献
19.
《物探化探计算技术》2016,(3)
井中激电是寻找深部金属矿产的有效方法之一,这里基于有限差分和异常电位算法实现了时间域井中激发极化法的三维数值模拟。首先直接给出点源三维异常电位场满足的基本微分方程;然后采用一维非零元素行压缩存储模式对所形成的大型稀疏矩阵进行存储,节约了计算所需内存空间;同时引入不完全Cholesky分解稳定化双共轭梯度(ICBG)法求解有限差分线性方程组,提高了求解的效率。构建典型的地电模型,采用地-井、井-地、井-井的观测模式分别进行正演计算,结果表明,时间域井中激电视极化率异常响应较视电阻率异常响应灵敏,综合考虑这两种视参量的空间分布特征,在一定程度上可以提高对相关异常体的分辨率,为今后野外的实际数据资料的解释提供依据。 相似文献
20.
ZTEM(Z轴倾子电磁法)是一种天然场源的频率域航空电磁法,其特点是磁场垂直分量在空中机载平台测量,磁场水平分量在地面的固定基站测量,具有勘探深度大、速度快、成本低、覆盖面积大等技术优势。本文实现了ZTEM二维有限差分正演和二维非线性共轭梯度(NLCG)反演算法。研究对象是倾子资料,反演过程中通过解"拟正演"问题来避开雅克比矩阵的直接计算。通过理论模型合成数据反演试算,验证了ZTEM倾子资料二维NLCG反演算法的稳定性与可靠性。与大地电磁(MT) TE模式阻抗资料反演结果进行对比,发现在异常体横向边界的约束方面,ZTEM倾子反演比MT阻抗反演更具优越性。 相似文献