共查询到20条相似文献，搜索用时 390 毫秒
1.
运用有限差分法研究频率域中二维弹性波动方程组的正演算法及理论合成记录的制作。频率域中二维弹性波动方程组的建立还未见报道，且此种算法要比时间域中正演算法计算量小得多，大大提高了计算速度，并且为弹性波动方程组的多参数反演提供了基本保证。理论合成记录的制作是获取波场信息的重要途径，可以全面地反映地震波在地下介质中的分布与传播情况 相似文献
2.
3.
频率-空间域正演模拟是频率域及Laplace-Fourier域全波形反演的基础,起伏地表条件下波形反演算法的关键是正演算法中考虑起伏地表的影响。基于带PML吸收边界的声波波动方程,在已有最优9点有限差分正演算法的基础上构建了起伏地表条件下频率-空间域正演算法。通过应用变网格技术,进一步提高算法的计算效率、降低内存开销,使得大规模起伏地表模型的频率域正反演问题成为可能。理论分析及数值测试表明：通过对近地表区域进行局部网格加密,可有效地压制由于矩形网格离散引起的角点散射;结合变网格技术可较易获得5倍以上计算效率的提高及内存占用的降低,且随着模型尺度的增加及地表起伏高程差的减小,倍数将显著增加;在细网格与粗网格交界处产生的虚假反射振幅幅值控制在原始波场的2%以内,满足地震波场正反演的需求。 相似文献
4.
5.
《物探与化探》2021,(3)
Z轴倾子电磁法(ZTEM)是一种新型频率域航空天然场源电磁法,该方法采用倾子作为研究参数,可用于起伏地形下的大规模地质勘探。本文在对TEM三维有限差分正演和数据空间OCCAM反演算法研究的基础上,考虑地形起伏的影响,研究了带地形的频率域三维ZTEM正反演算法。首先,对起伏地表下ZTEM正演算法的准确性进行了验证,计算和分析了纯地形的ZTEM三维异常响应特征。其次,通过对山峰和山谷地形下低阻棱柱体合成算例的反演分析,表明带地形的ZTEM数据空间OCCAM反演算法能够获得比较接近真实导电性结构的地电模型,尤其对地下目标体的横向边界具有理想的约束效果。最后,与不带地形的ZTEM反演结果对比,验证了所开发的带地形ZTEM倾子资料反演算法的有效性。 相似文献
6.
起伏地形频率域可控源电磁二维快速正反演 总被引：1，自引：0，他引：1
起伏地形频率域可控源电磁正反演研究受到众多学者的关注和重视,目前普遍实用的为未考虑地形的二维正反演算法程序。本文基于积分方程法正演和多场源、多频率、对比源反演算法,开展起伏地形条件下频率域可控源电磁二维快速正反演研究。通过引入层状参考模型、将地形与目标体整体作为异常场剖分区域,实现起伏地形可控源电磁正反演计算。针对包含地形在内的大尺度剖分区域引起的大型计算代价问题,采用快速傅里叶算法提高正反演计算效率。通过与现有正反演算法进行模型算例对比,验证了本文采用方法的可行性与有效性。 相似文献
7.
《物探化探计算技术》2017,(1)
地震波场频率域正演是频率域全波形反演的基础。针对反演计算量巨大的问题,利用频率域二维波形正演算法中频率和炮点计算的独立性,开发出粗细粒度结合的MPI并行算法:粗粒度为频率并行,细粒度为单个频率解方程并行。实现方法是将正演频率分组放入不同的MPI通讯域内,每个通讯域内单个频率求解方程过程,采用基于MPI的MUMPS(多波前大规模并行稀疏直接解法器)软件包并行加速。模型测试结果表明:MPI多通讯域并行算法计算结果正确,计算效率显著提高,加速效果稳定。 相似文献
8.
9.
各向异性介质中地震波场计算的研究仍较少。本文根据Ｆ一Ｋ域格林函数的解析解，采用玻恩近似及积分法，计算了二维层状、横向各向同性介质中的合成记录。 相似文献
10.
ZTEM(Z轴倾子电磁法)是一种天然场源的频率域航空电磁法,其特点是磁场垂直分量在空中机载平台测量,磁场水平分量在地面的固定基站测量,具有勘探深度大、速度快、成本低、覆盖面积大等技术优势。本文实现了ZTEM二维有限差分正演和二维非线性共轭梯度(NLCG)反演算法。研究对象是倾子资料,反演过程中通过解"拟正演"问题来避开雅克比矩阵的直接计算。通过理论模型合成数据反演试算,验证了ZTEM倾子资料二维NLCG反演算法的稳定性与可靠性。与大地电磁(MT) TE模式阻抗资料反演结果进行对比,发现在异常体横向边界的约束方面,ZTEM倾子反演比MT阻抗反演更具优越性。 相似文献
11.
《物探化探计算技术》2016,(3)
对三维场源二维地电模型的正演计算称为2.5维数值模拟。这里从麦克斯韦方程组出发,1分别求解电磁场的一次场和二次场,将三维场源降为二维;2利用傅氏变换将空间域方程转化为波数域,应用有限元求解波数域电磁场方程,引入无限元解决无穷远边界的收敛问题;3根据电磁场值实虚部的曲线特征,按对数等间隔选取21个波数。编写代码计算均匀半空间与解析解结果对比,电阻率的均方相对误差均小于0.5%,证明其有效性。计算了三种地电模型的电磁场响应,对单一高、低阻体和高低阻组合体的模拟效果真实,异常中心位置基本吻合。结果证明,这里的CSAMT2.5-D正演算法可以模拟较为复杂的地电模型,并取得良好的效果。 相似文献
12.
研究了二维横向各向同性介质中关于位移分量u_x,u_y与u_z的偏微分波动方程组进行三分量弹性波场有限差分法模拟的稳定性问题,分别获得了y向分量与x、z向分量正演模拟时所必须满足的稳定性条件,提供了三分量记录数字仿真稳定所需的时空离散间隔约束。 相似文献
13.
通过ｚ－τ变换，把（ｘ，ｚ，ｔ）域的波动方程变换为（ｘ，τ，ｔ）域的波动方程，以分离出波动方程中速度（或弹性系数）的变化项，并利用有限差分法来计算横向各向同位介质中的合成记录。差分格式采用二阶中心差商及一阶向前差商，边界条件为吸能边界。计算结果显示了横向各向同位介质中的ｑＰ、ｑＳＶ及ｑＳＨ波；在有界面的情形下，还较清楚地反映了反射波的存在。 相似文献
14.
二维横向各向同性介质中的差分法地震波正演 总被引：1，自引：0，他引：1
通过ｚ－τ变换，把（ｘ，ｚ，ｔ）域的波动方程变换为（ｘ，τ，ｔ）域的波动方程，以分离出波动方程中速度（或弹性系数）的变化项，并利用有限差分法来计算横向各向同性介质中的合成记录。差分格式采用二阶中心差商及一阶向前差商，边界条件为吸能边界。计算结果显示了横向各向同性介质中的ｑＰ、ｑＳＶ及ｑＳＨ波；在有界面的情形下，还较清楚地反映了反射波的存在。 相似文献
15.
针对弹性参数满足一种特殊形式的二维各向异性介质，结合边界条件，求解波动方程组，得到了该种介质中传播的准瑞利波的波速度满足的方程，进而探讨了该种介质中传播的准瑞利面波的传播规律。 相似文献
16.
零偏移距地震记录是叠后偏移的基础,常用射线追踪法来制作,但当介质复杂时,射线会出现交叉现象,对三维情况计算量更大。从二维标量波动方程出发,用伪谱法研究了零偏移距地震记录的合成,给出了不同情况下的数值计算公式,并对边界吸收问题作了有效处理。数值计算结果表明了方法的正确性和有效性。该方法可用于偏移中正演模拟的快速计算,且易推广到三维情形。 相似文献
17.
《物探化探计算技术》2015,(4)
近年来2.5维复电阻率反演取得了一些理论和应用方面的成果,但计算过程复杂、计算量巨大,特别是正演计算过程中计算量大导致了反演速度较慢。这里在2.5维复电阻率反演计算基础上,将反演算法中的正演过程的排列和频率循环分割成多个规模较小的程序段,并发执行正演过程中的每个排列和频率循环的求解计算过程,并根据并行计算的节点数目或者计算核心数目,动态调整每个节点或者计算核心的计算量,从而提高了算法的计算效率和适应性。 相似文献
18.
《物探与化探》2017,(2)
为优化二维各向同性介质中弹性波频率域正演时阻抗矩阵的结构,减小正演所需内存,提高正演效率,在25点差分格式的基础上进行适当的简化,得到了二维弹性波频率域17点差分格式。该格式重新计算了弹性波中偏微分项和加速项的差分算子,减少了计算过程中的网格节点需求,构造了优化阻抗矩阵后的频率域正演矩阵方程;推导了纵波和横波相速度的频散公式,给出了不同泊松比条件下的频散曲线,得到了相速度误差控制范围±1%时每一横波波长内网格数需求。通过对比频散曲线和数值模拟时得到的波场快照及检波点处U、V分量,验证了17点差分格式与25点差分格式相比,具有稍严格的网格间距需求、相当的计算精度、略少的计算时间和更小的阻抗矩阵带宽等特点。 相似文献
19.
20.
这里实现了在2.5维各向异性介质中地震波传播的数值模拟。首先从2.5维一阶速度-应力弹性各向异性波动方程出发,得出了该方程的拟谱法数值解法;然后通过对Marmousi速度模型进行数值模拟,表明了拟谱法的有效性;最后计算并分析了均匀和混合各向异性介质中波的传播快照和理论模拟的地震记录,进一步认识了波在各向异性介质中的传播规律。 相似文献