共查询到20条相似文献，搜索用时 125 毫秒
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.
针对弹性参数满足一种特殊形式的二维各向异性介质，结合边界条件，求解波动方程组，得到了该种介质中传播的准瑞利波的波速度满足的方程，进而探讨了该种介质中传播的准瑞利面波的传播规律。 相似文献
7.
这里实现了在2.5维各向异性介质中地震波传播的数值模拟。首先从2.5维一阶速度-应力弹性各向异性波动方程出发,得出了该方程的拟谱法数值解法;然后通过对Marmousi速度模型进行数值模拟,表明了拟谱法的有效性;最后计算并分析了均匀和混合各向异性介质中波的传播快照和理论模拟的地震记录,进一步认识了波在各向异性介质中的传播规律。 相似文献
8.
9.
起伏地形频率域可控源电磁二维快速正反演 总被引：1，自引：0，他引：1
起伏地形频率域可控源电磁正反演研究受到众多学者的关注和重视,目前普遍实用的为未考虑地形的二维正反演算法程序。本文基于积分方程法正演和多场源、多频率、对比源反演算法,开展起伏地形条件下频率域可控源电磁二维快速正反演研究。通过引入层状参考模型、将地形与目标体整体作为异常场剖分区域,实现起伏地形可控源电磁正反演计算。针对包含地形在内的大尺度剖分区域引起的大型计算代价问题,采用快速傅里叶算法提高正反演计算效率。通过与现有正反演算法进行模型算例对比,验证了本文采用方法的可行性与有效性。 相似文献
10.
《物探化探计算技术》2017,(1)
地震波场频率域正演是频率域全波形反演的基础。针对反演计算量巨大的问题,利用频率域二维波形正演算法中频率和炮点计算的独立性,开发出粗细粒度结合的MPI并行算法:粗粒度为频率并行,细粒度为单个频率解方程并行。实现方法是将正演频率分组放入不同的MPI通讯域内,每个通讯域内单个频率求解方程过程,采用基于MPI的MUMPS(多波前大规模并行稀疏直接解法器)软件包并行加速。模型测试结果表明:MPI多通讯域并行算法计算结果正确,计算效率显著提高,加速效果稳定。 相似文献
11.
针对起伏地表条件下的地震波数值模拟问题,提出将起伏的地表映射到一个规则的长方形网格坐标系中,并在此基础上推导出变换域中的波动方程。根据导出来的波动方程应用有限余弦变换有限差分方法进行地震波场的数值模拟,以解决起伏地表条件下的数值模拟。 相似文献
12.
基于势流理论的数值水池可以快速计算波浪的传播及其对建筑物的作用,但是势流理论是基于波浪的无黏性假设的,而在工程中,通常需要在固体边界及波浪破碎的区域考虑黏性效应。针对基于求解Navier-Stokes方程的黏性水池计算量较大、速度较慢的缺点,采用耦合的方法模拟多向不规则波浪的传播,即在外域通过基于势流理论的数值水池产生多向不规则波浪,内域采用求解Navier-Stokes方程和流体体积法(VOF方法)对自由表面进行追踪,通过外域所提供的边界波浪,内域计算可以在较小区域进行计算,从而达到减少计算量、提高计算效率的目的。 相似文献
13.
针对地基隔振控制,提出了一类新型的空沟?波阻板联合隔振屏障,并对其隔振性能进行了数值分析。首先,利用复伸展坐标变换,在频域内建立了完全匹配层（perfect matched layer, PML）吸收边界的控制方程;其次,利用Galerkin近似技术,给出了以位移为基本未知量的二阶非分裂格式PML的频域有限元计算列式;最后,通过数值算例分析了空沟?波阻板联合隔振屏障的物性参数（地基与波阻板的模量比）、几何参数（空沟深度、波阻板深度）以及载荷参数（振动波频率）等对其隔振性能的影响规律。结果表明,空沟?波阻板联合隔振屏障结合了空沟和波阻板各自的优势,可以有效地控制不同频率振源引起的地基振动。 相似文献
14.
Acoustic imaging and sensor modeling are processes that require repeated solution of the acoustic wave equation. Solution of the wave equation can be computationally expensive and memory intensive for large simulation domains. One scheme for speeding up solution of the wave equation is the operator-based upscaling method. The algorithm proceeds in two steps. First, the wave equation is solved for fine grid unknowns internal to coarse blocks assuming the coarse blocks do not need to communicate with neighboring blocks in parallel. Second, these fine grid solutions are used to form a new problem which is solved on the coarse grid. Accurate and efficient wave propagation schemes also must avoid artificial reflections off of the computational domain edges. One popular method for preventing artificial reflections is the nearly perfectly matched layer (NPML) method. In this paper, we discuss applying NPML to operator upscaling for the wave equation. We show that although we only apply NPML to the first step of this two step algorithm (directly affecting the fine grid unknowns only), we still see a significant reduction of reflections back into the domain. We describe three numerical experiments (one homogeneous medium experiment and two heterogeneous media examples) in which we validate that the solution of the wave equation exponentially decays in the NPML regions. Numerical experiments of acoustic wave propagation in two dimensions with a reasonable absorbing layer thickness resulted in a maximum pressure reflection of 3–8%. While the coarse grid acceleration is not explicitly damped in our algorithm, the tight coupling between the two steps of the algorithm results in only 0.1–1% of acceleration reflecting back into the computational domain. 相似文献
15.
假定各物理量沿径向不变化,建立了低应变瞬态集中荷载作用下现浇大直径管桩振动响应的计算模型和波动方程。采用Laplace变换法,求得了波动方程的频域解析解,采用Fourier逆变换求得了时域响应。将文中解计算结果与三维频域解析解进行了对比分析,文中解的入射峰-反射峰时间差与三维频域解差别很小,说明对于PCC桩这种大直径薄壁管桩,采用径向不变假定对计算结果几乎没有影响。因此,在PCC桩低应变检测波形分析时,采用基于径向不变假定的二维解是完全合理的。将文中二维频域解的计算结果与二维时域解的结果进行了对比分析,结果表明,2种解在入射波、第一个桩底反射波、第二个桩底反射波峰值大小和到达时间非常吻合,这说明考虑环向位移与否对计算结果没有太大的影响,研究PCC桩低应变动力响应的问题时忽略水平方向的位移是完全可以的。 相似文献
16.
17.
TTI介质弹性波频率-空间域有限差分数值模拟 总被引：3，自引：0，他引：3
由周期性薄互层引起的VTI介质是研究比较广泛的一类各向异性介质。当VTI介质对称轴偏离垂向,本构坐标系与观测坐标系不重合时,会形成观测坐标系下的TTI介质。引入25点优化差分算子,推导出二维TTI介质频率域弹性波动方程;为压制边界反射,采用完全匹配层法吸收边界条件,并计算出优化差分系数;最后采用集中力源,模拟了弹性波在TTI介质中的传播过程。从波场快照和地面共炮记录可以看出,笔者采用的数值模拟算法能有效压制数值频散。TTI介质中的波场传播比较复杂,纵波传播相对稳定,横波波前的三分叉现象比较明显,并存在振幅奇异性。当VTI介质的对称轴偏转后,还会增加地面地震记录的复杂性。 相似文献
18.
波动方程的系数反演是一种用于识别地下介质物理力学参数的重要方法。过去这方面的研究一般都建立在弹性模型的基础上,本文则提出了粘弹性反演模型,这种模型更能真实地反映波在地下传播的实际情形。在粘弹性反演模型的基础上,我们还提出了一整套行之有效的数值反演方法,并在频域中完成了反演计算。最后的数值计算结果是令人满意的,证明了本文的模型和方法的合理性。 相似文献
19.
动力响应问题的求解对于地基在外荷载作用下引起的弹性波动问题研究有重要的意义。本文提出了一种求解横观各向同性层状地基在施加时间简谐荷载作用下任意点的应力响应的算法。此算法利用傅里叶变换将广义平面应变问题频率-空间域的动力方程转化到频率-波数域内,结合对偶变量的引入,利用高精度的精细积分算法对状态方程进行求解,在得到频率-波数域内的位移响应的基础上,利用傅里叶逆变换得到任意点的动应力响应。简谐荷载不仅可以施加在地基表面,而且可以施加在地基内部。对比算例验证了本文算法的准确性,同时对地基各向异性特性、激励频率和阻尼比对动应力响应的影响进行了参数分析,为工程实际提供可靠的数值依据。 相似文献
20.
正演数值模拟是反演成像的基础。为了实现磁法勘探精细化反演成像与定量解释,本文利用把一个大问题分解为多个小问题的思路,提出一种基于泊松方程的高效、高精度空间波数混合域二度体磁异常数值模拟方法。该方法利用傅里叶变换把二度体磁位偏微分方程转换为一维常微分方程,并采用基于二次插值的一维有限单元法求解该方程,进而通过反傅里叶变换得到空间域磁异常。在模型算例中,分别设计截面为矩形的常磁化率和变磁化率二度体模型,针对本文算法的计算精度和计算效率进行了验证。模型算例结果表明：该算法计算精度高,相对误差绝对值均小于1%;计算速度快,网格节点剖分2 501×2 501的模型模拟时间为4.18 s;适用于任意复杂地形模型。 相似文献