共查询到20条相似文献,搜索用时 156 毫秒
1.
在地震勘探中,地震正演模拟是非常重要的技术。与时间域正演相比,频率域正演速度快,计算效率高。如何高效准确地完成频率域正演计算是目前该领域的一个重要问题。数值频散问题和如何提高计算效率降低求解分解阻抗内存占用量一直是频率域正演所需要解决的问题。与传统的直接法求解阻抗矩阵的频率域正演方法不同,本文采用最小二乘共轭梯度法(LSCG法)求解阻抗矩阵进行频率域正演,并提出了一种波数补偿的表达式来压制数值频散现象。经过简单模型和复杂模型的数值测试,采用最小二乘共轭梯度法(LSCG法)求解阻抗矩阵进行频率域正演能够有效降低计算时间,且采用波数补偿的频率域正演方法能够有效压制数值频散现象,提高波场模拟精度。 相似文献
2.
电磁场从频率域转换到时间域的几种算法比较 总被引:8,自引:1,他引:8
研究了电磁场从频率域向时间域转换的三种具体算法,即余弦变换多项式近似法。G-S法和数值线性滤波法;用这三种变换法对电偶源电场的瞬变响应E_(x2)(t)进行了计算;将计算结果同解析解进行了比较,对比了三种算法的计算精度。 相似文献
3.
研究了提高频率域中导数计算结果的精度以及局部波数反演结果的稳定性和精度,对比了余弦扩边方法和最小曲率扩边方法在频率域中计算导数的精度。通过对比表明,最小曲率扩边结果的精度高于余弦扩边结果的精度,提出了向上延拓某一高度来进行局部波数的反演技术。此技术不但能够解决高频干扰,提高导数计算精度,而且利用不同延拓高度的反演结果能够得到一个平稳的反演结果,从而提高了局部波数反演的精度。理论模型试算验证了方法技术的精度,实际资料处理验证了方法技术的稳定性,达到了满意的效果。 相似文献
4.
5.
通常的叠前深度偏移都是通过偏移距域共成像道集和炮域共成像道集来实现的,这里则通过对G-D波场局部平面波的分解,提出了局部角度域共成像点道集方法,从而改善了偏移成像速度及效果,为研究地下构造,岩性分析以及偏移速度分析提供新的途径。 相似文献
6.
三维反演在磁数据定量解释中具有重要作用。常用的空间域三维反演方法通常需要大量的正演和反演计算,因此对大规模数据的反演效率较低。三维成像是另一种定性和定量解释磁数据的重要方法。文中给出了一种磁异常与梯度三维成像的频率域迭代方法,该方法可以提高成像效率,适用于大规模数据的快速成像。笔者推导了磁总场异常和异常梯度频率域正演公式和成像公式,并将一种深度尺度因子引入成像公式中,提高了深度精度;笔者采用了迭代优化算法,减小了拟合误差,进一步提高了成像精度。通过理论模型数据试验和中国新疆某金属矿床实测数据,验证了本文方法的有效性、可行性。 相似文献
7.
垂向导数广泛应用于位场数据处理和解释中,常规的空间域求导公式多以台劳级数为基础,这不可避免会引入截断误差,且不便于求取高阶导数。另一种应用较多的波数域算法。其求导算子在高波数成份产生振荡失真,很不稳定,在实际应用中很难取得理想结果。这里提出了一种计算位场数据垂向导数的稳定算法,在波数域求垂向积分,在空间域计算二个二阶水平方向导数,借助拉普拉斯方程,求取位场异常的各阶垂向导数。算法只在求垂向积分时引入了一次付里叶变换,其它计算全部在空间域中进行,从而有效地避免了吉布斯效应的重复累加传递,有效降低了常规波数域求导算子在高波数成份产生的振荡失真。通过理论模型对比,证明了这里所提出的算法较常规的波数域方法具有良好的稳定性,特别是在计算高阶垂向导数时有明显优势,能有效地提高与垂向导数有关的位场定量解释方法的准确性。 相似文献
8.
均匀介质中的F-K反偏移是实现一般介质下F-K反偏移的基础,而均匀介质中的F-K反偏移的实质是一个插值映射过程。根据插值对象为复函数的特点,采用了实部虚部分别插值、去振荡因子法、模和相位分别插值、直接频率域插值4种插值方法对点脉冲模型进行了试算。结果表明:截断的直接频率域插值在适应性和精度综合衡量下为最佳;去振荡因子法精度虽高,但可应用的模型太具理想,不适合一般模型中的应用;实部虚部分别插值与模和相位分别插值两者精度较低。 相似文献
9.
10.
角度域照明可以更加细致地对地下复杂地质体进行照明分析。提出了一种适用于起伏地表条件下的角度域照明分析方法,其主要过程是利用局部倾斜叠加技术分解波场,对起伏地表处采用加权函数判定局部时窗内的振幅包络位置,以区分其位于地表之上或之下;并引入相似系数函数计算反假频局部平面波场,根据局部斜率与平面波场入射角关系,得到局部角度域地震波场;根据帕斯瓦尔等式下的时间-角度域照明关系式,计算得出角度域照明强度。通过起伏地表高速侵入矿体模型试算,得到给定观测系统下角度域照明强度与地下分析点照明随炮点、波场入射角的变化关系。此方法可用于起伏地表条件的观测系统优化、叠前AVA分析等方面。 相似文献
11.
12.
面波频散能量谱成像精度是近地表横波速度计算的重要影响因素,有效提取面波频散曲线是后续反演地层结构及介质物性参数的重要环节。文中论述了τ-p变换法(τ为垂直波慢度,p为水平波慢度)、高分辨率线性Radon变换法、频率分解法、相移法、f-k变换法(f为频率,k为波数)、SFK变换法和矢量波数变换法等7种不同面波频散能量谱计算方法的基本原理;通过上述方法,对三层速度递增水平层状介质模型正演模拟合成Rayleigh波记录进行频散能量谱计算,并与理论频散曲线进行比较;根据所得频散能量谱,分析不同方法在数据处理过程中的特点与应用差异。结果表明:在提取频散曲线过程中,高分辨率线性Radon变换法和矢量波数变换法提取精度较高,f-k变换法提取精度最低,低频信息与高频信息无法有效识别;SFK变换法、矢量波数变换法、高分辨率线性Radon变换法抗噪性能强;τ-p变换法、f-k变换法、相移法计算速度优于其他方法。 相似文献
13.
The method of fundamental solution for 3‐D wave scattering in a fluid‐saturated poroelastic infinite domain 下载免费PDF全文
Zhongxian Liu Zhikun Wang Alexander H.D. Cheng Jianwen Liang Chuchu Wang 《国际地质力学数值与分析法杂志》2018,42(15):1866-1889
By using a complete set of poroelastodynamic spherical wave potentials (SWPs) representing a fast compressional wave PI, a slow compressional wave PII, and a shear wave S with 3 vectorial potentials (not all are independent), a solution scheme based on the method of fundamental solution (MFS) is devised to solve 3‐D wave scattering and dynamic stress concentration problems due to inhomogeneous inclusions and cavities embedded in an infinite poroelastic domain. The method is verified by comparing the result with the elastic analytical solution, which is a degenerated case, as well as with poroelastic solution obtained using other numerical methods. The accuracy and stability of the SWP‐MFS are also demonstrated. The displacement, hoop stress, and fluid pore pressure around spherical cavity and poroelastic inclusion with permeable and impermeable boundary are investigated for incident plane PI and SV waves. The scattering characteristics are examined for a range of material properties, such as porosity and shear modulus contrast, over a range of frequency. Compared with other boundary‐based numerical strategy, such as the boundary element method and the indirect boundary integral equation method, the current SWP‐MFS is a meshless method that does not need elements to approximate the geometry and is free from the treatment of singularities. The SWP‐MFS is a highly accurate and efficient solution methodology for wave scattering problems of arbitrary geometry, particularly when a part of the domain extends to infinity. 相似文献
14.
The authors propose a semi‐analytical approach to studying wave propagation in multilayered poroviscoelastic grounds due to transient loads. The theoretical development is based on the exact stiffness matrix method for the Biot theory coupled with a matrix conditioning technique. It is developed in the wavenumber frequency domain after a Fourier transform on the surface space variables and the time variable. The usual methods yield a poorly conditioned numerical system. This is due in particular to the presence of mismatched exponential terms. In this article, increasing exponential terms are eliminated and only decreasing exponential terms remain. Consequently, the method can be applied to a large field of configurations without restriction concerning high frequencies, large Fourier transform parameters or large layer thicknesses. Validation and efficiency of the method are discussed. Effects of layering show that the layer impedance influence on solid and fluid displacements. Moreover, this approach can be of interest for the validation of numerical tools. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
15.
合成平面波数据的Rytov近似叠前深度偏移 总被引:6,自引:0,他引:6
作者把Berkhout(1992)提出的合成平面波技术应用于线性Rytov近似叠前深度偏移,提高了波动方程叠前深度偏移的计算效率。通过在Marmousi模型数据试验,我们认为在达到基本相同成像效果的条件下,利用合成平面波技术的波动方程叠前深度偏移相对于基本共炮集数据的波动方程叠前深度偏移可提高计算效率约10倍。 相似文献
16.
频率-空间域正演模拟是频率域及Laplace-Fourier域全波形反演的基础,起伏地表条件下波形反演算法的关键是正演算法中考虑起伏地表的影响。基于带PML吸收边界的声波波动方程,在已有最优9点有限差分正演算法的基础上构建了起伏地表条件下频率-空间域正演算法。通过应用变网格技术,进一步提高算法的计算效率、降低内存开销,使得大规模起伏地表模型的频率域正反演问题成为可能。理论分析及数值测试表明:通过对近地表区域进行局部网格加密,可有效地压制由于矩形网格离散引起的角点散射;结合变网格技术可较易获得5倍以上计算效率的提高及内存占用的降低,且随着模型尺度的增加及地表起伏高程差的减小,倍数将显著增加;在细网格与粗网格交界处产生的虚假反射振幅幅值控制在原始波场的2%以内,满足地震波场正反演的需求。 相似文献
17.
波动方程角度域共成像道集 总被引:1,自引:0,他引:1
针对基于波场深度递推的波动方程叠前深度偏移成像方法难以准确地给出如Kirchhoff积分叠前深度偏移成像方法的偏移距域共成像道集的不足,在分析研究现有的各种波动方程偏移成像共成像道集方法的基础上,利用波场外推的单平方根算子和波场的窗口Fourier框架展开与重构方法,提出一种局部角度域共成像道集方法——成像点处反射角共成像道集方法.把这种角度域共成像道集方法应用于国际标准的Marmousi模型数据和一条实际二维地震数据都取得了理想的结果.提出的波动方程反射角度域共成像道集方法可为进一步的叠前偏移数据振幅随角度变化分析和偏移速度分析工作提供基础. 相似文献
18.
偏移距域/角度域共成像点道集与偏移速度的关系 总被引:3,自引:0,他引:3
为了得到波动方程偏移距域共成像点道集(ODCIGs) 和角度域共成像点道集(ADCIGs) 与偏移速度的关系, 在匀速单层水平反射情况下, 对ODCIGs和ADCIGs随偏移速度的变化进行了定量的推导, 得到偏移速度偏小时ODCIGs与半偏移距成双曲关系, 反之成椭圆关系; 结合深度聚焦分析, 得到偏移速度偏大时, ADCIGs中视入射角大于真入射角, 反之, 小于真入射角; 在速度偏小时, ADCIGs上剩余时差(RMO) 与视入射角的正切值成椭圆关系, 反之成双曲关系.在地震有效入射角度范围内, ADCIGs对于偏大的偏移速度更加敏感.在偏移速度分析(MVA) 中, 从稍大的速度开始分析会对速度分析更加有利. 相似文献
19.
井间地震资料中蕴含丰富的地下储层岩性和物性信息,利用井间地震直达波初至信息,通过层析反演可以得到两井之间的速度剖面。通过对井间地震直达波振幅衰减信息的处理,可以实现井间的吸收衰减层析成像。基于井间地震直达波射线路径的一致性,提出了井间地震时域直达波走时层析反演和频域质心频率衰减层析同步反演方法,增加了约束条件,增强了抗干扰能力。对模型和实际资料的处理结果表明,速度和吸收衰减同步反演方法,提高了速度层析反演精度,得到了对油气更为敏感的衰减属性,增强了井间地震在储层描述和油气预测中的应用能力。 相似文献