共查询到20条相似文献，搜索用时 15 毫秒
1.
完全匹配层吸收边界（PML）已经被证明是非常有效的边界吸收技术,对体波和面波的吸收都具有非常好的效果,已经被广泛应用于弹性波的数值模拟中。但是在某些情况下传统的PML技术还是存在一定的问题,比如对掠射情况下的体波和窄区域自由表面条件下的面波的吸收等等。在坐标变换中采用复频移拉伸函数的复频移PML可以有效地改善PML边界条件的吸收性能。基于弹性波一阶速度-应力方程,推导了复频移PML的递推卷积实现方法,并采用交错网格高阶有限差分法对其进行了数值模拟,与传统的PML进行了对比。结果表明：传统的PML对掠射情况下的体波和窄区域自由表面条件下的面波吸收不足,会产生虚假反射,影响真实波场;而基于递推卷积的复频移PML算法能够有效地改善困难情况下的吸收效果,并且在实现过程中不用分裂变量,应用更加方便简单。计算卷积时采用递推的形式,推导过程直观易懂,易于编程,而且不会增加计算量,存储量也没有太大的变化。 相似文献
2.
弹性波正演模拟中PML吸收边界条件的改进 总被引：4，自引：0，他引：4
在弹性波有限差分正演模拟中, 完全匹配层(PML) 吸收边界条件是使用广泛、吸收效果最好的吸收边界条件.在目前的两种PML实现方法中, 分裂形式的完全匹配层(SPML) 方法计算存储量大、编程实现复杂; 非分裂形式的完全匹配层(NPML) 方法计算效率低、计算过程复杂.针对传统PML吸收边界条件在实现过程中存在的问题, 推导出了一种简洁有效的非卷积实现的NPML吸收边界条件, 既不需要对场分量进行分裂, 也不需要做复杂的卷积运算.分析结果表明, 本文实现的NPML吸收边界条件不仅具有良好的吸收衰减性能, 而且计算方程简单, 编程实现容易, 占有内存更小. 相似文献
3.
曲线坐标系下的完全匹配层吸收边界条件 总被引：1，自引：0，他引：1
在地震波数值模拟中,需要采用吸收边界条件以吸收人为边界反射。本文针对曲线坐标系下的二阶弹性波方程提出了一种完全匹配层(PML)吸收边界条件。与直角坐标系下的PML吸收边界条件类似,曲线坐标系下的PML吸收边界条件是一种在频率域中给出的人工边界条件,由相应的复坐标变换得到。在变换到时间域后,完全匹配层中将出现复杂的卷积运算。为了避免这些卷积运算,引入了4个中间变量。为了简化自由边界条件,采用正交贴体网格对起伏地表模型进行网格剖分。数值算例表明,该方法可以有效消除人为边界反射。 相似文献
4.
多个激励源无延时发射(多源并发)相同中心频率脉冲会形成平面波束信号,增强数据记录质量。本文通过数值模拟对比分析在多源并发情况下,非分裂递归卷积完全匹配层作为吸收边界条件和Mur二阶吸收边界条件对电磁波的吸收效果。其研究结果表明,传统的Mur二阶吸收边界条件对多源并发、多角度掠射情况下电磁波的吸收效果不佳,在大偏移距下会造成波形畸变和形成虚假反射。而在多源并发情况下采用非分裂递归卷积完全匹配层作为吸收边界条件,将坐标伸缩因子引进时域有限差分算法中。通过傅里叶逆变换将频率域坐标伸缩变换PML方程转换到时域,对电场和磁场值在离散状态下进行递归卷积运算求解。从而避免了直接对卷积进行数值求解的复杂计算,在保证计算准确性的同时,节约了内存空间,提高了计算效率。在不分裂波场情况下,改善了网格截断位置对电磁波的吸收效果。 相似文献
5.
地震数值模拟中,完全匹配层(PML)边界能有效地吸收衰减地震波,得到无边界反射干扰的波场快照和地震记录。在前人研究的基础上,进一步分析了分裂完全匹配层(SPML)的衰减机制,通过波动方程模拟,分析了震源主频、空间网格间距、介质速度等参数对SPML边界吸收衰减特征的影响。得到了不同条件下,PML边界对地震波的吸收效果。通过分析对比,得出了震源主频对PML吸收效果无直接影响,空间网格间距与PML吸收效果成反比,高速层PML吸收效果缓慢等结论。并通过Marmousi2模型测算,对结论进行了验证。最后给出了复杂模型、不同尺度数值模拟中SPML参数的选取方法,为地震数值模拟中SPML边界参数的定量选取奠定了基础。 相似文献
6.
地质雷达数值模拟中有损耗介质吸收边界条件的实现 总被引：5，自引：1，他引：4
利用时间域有限差分(FDTD)法将麦克斯韦方程进行离散化,可以对地质雷达进行数值模拟。利用完全匹配层(PML)作为吸收边界条件可以有效地吸收向外的电磁波,从而大大提高了计算效率。对于有损耗介质的情况,采用扩张坐标系下改正的麦克斯韦方程,只要扩张变量满足一定的关系,同样可有效地吸收向外的电磁波, 这种吸收边界条件称为通用完全匹配层。 相似文献
7.
针对地基隔振控制，提出了一类新型的空沟?波阻板联合隔振屏障，并对其隔振性能进行了数值分析。首先，利用复伸展坐标变换，在频域内建立了完全匹配层（perfect matched layer, PML）吸收边界的控制方程；其次，利用Galerkin近似技术，给出了以位移为基本未知量的二阶非分裂格式PML的频域有限元计算列式；最后，通过数值算例分析了空沟?波阻板联合隔振屏障的物性参数（地基与波阻板的模量比）、几何参数（空沟深度、波阻板深度）以及载荷参数（振动波频率）等对其隔振性能的影响规律。结果表明，空沟?波阻板联合隔振屏障结合了空沟和波阻板各自的优势，可以有效地控制不同频率振源引起的地基振动。 相似文献
8.
从一阶速度—应力弹性波方程出发,基于旋转交错网格,推导了时间二阶精度空间2M阶精度的有限差分离散格式。阐述了递归卷积复频移完全匹配层(CPML)边界条件的原理,建立了一阶速度—应力弹性波高阶差分CPML边界条件的递推公式。开展了CPML边界中关键参数m、κ和α的选取实验,通过分析反射误差分布图,选取了CPML边界条件中最优参数。全局反射误差与波场快照都说明,CPML较PML对隐失波具有更优的吸收性能。基于Matlab平台,编写了基于CPML边界的旋转交错网格弹性波正演模拟程序,应用该程序对各向异性介质及随机介质进行了模拟,得到了弹性波正演剖面记录及波场快照,通过对正演剖面记录及波场快照的分析,可以更清楚地了解弹性波在各向异性介质及随机介质的传播特性,指导非均匀介质中地震勘探资料解释。 相似文献
9.
研究井间地震波场的形成过程以及波场的传播机理、规律,对于指导实际井间地震勘探有着重要的意义.基于具有垂直对称轴的横向各向同性(VTI)介质中的一阶准P波方程,应用具有无条件稳定性质的紧致交错网格隐式差分格式求解该方程.重点研究了紧致交错网格求解该方程的完全匹配层(PML)吸收边界条件,在此基础上实现了VTI介质中一阶准P波方程的井间地震波场的正演模拟.数值算例表明:紧致交错网格能精准模拟VTI介质中准P波的传播过程,得到高精度的正演结果.一阶准P波方程能以足够的精度描述VTI介质中准P波特征.完全匹配层吸收边界能有效地解决人工边界问题,是一种高效的边界吸收算法. 相似文献
10.
11.
基于一阶速度-应力波动方程,采用高阶交错网格有限差分数值模拟方法,对弹性及黏弹性TTI介质进行正演数值模拟。模拟时采用完全匹配层吸收边界条件(PML)消除边界反射。同时设计了层状介质模型、断层模型,通过模型的正演计算,得到了不同时刻的地震波波场快照及合成地震记录,分析其波场运动学及动力学特征。模拟结果表明,交错网格有限差分法可以很好地完成对复杂介质的波场模拟,具有较高的精度和可靠性。 相似文献
12.
褶积微分算子法是一种全新的数值模拟方法,已被广泛应用于复杂介质的地震波场数值模拟,但是其边界反射问题一直没有解决。这里将最佳匹配层(PML)吸收边界条件引入到褶积微分算子法中,此方法是在研究区域的边界上加入吸收层,使边界上传入吸收层的波,随传播距离按指数规律衰减,不产生任何反射,以达到消除边界反射的目的。构造不同的模型,通过对比分析证明,PML吸收边界条件能比较好地解决褶积微分算子法的边界问题,从而验证了完全匹配层吸收效果的优越性。 相似文献
13.
基于散度和旋度的弹性波波场分离数值模拟方法 总被引：4，自引：1，他引：3
弹性波波场分离数值模拟方法是研究混合波场中纯纵波和纯横波波场传播规律的一种重要手段,通常采用算子分裂来实现,这不但增加了计
算量,而且还浪费了大量计算时间。因此提出了基于散度和旋度的弹性波波场分离数值模拟新方法,仅依赖于交错网格来计算出散度场（纯纵波）
和旋度场（纯横波）,从而提高计算效率。在人工截断边界处采用完全匹配层吸收边界条件（PML）。数值模拟结果表明,该方法能够从混合波场中
准确地分离出纯纵横波波场,且精度高,边界吸收效果好,计算速度快,计算过程稳定,可用于指导实际多波多分量地震资料处理,同时指出本文
算法对相互耦合的波型（如面波）无法实现准确分离。 相似文献
14.
以具有倾斜对称轴的横向各向同性(Titled Transverse Isotropic,TTI)介质中纯准P波二阶方程为基础,通过引入辅助波场实现了方程的降阶,推导出了TTI介质纯准P波一阶压力-速度方程,并利用旋转交错网格高阶有限差分方法对方程进行数值求解,给出了完全匹配层(Perfectly Matched Layers,PML)吸收边界条件在该方程中应用的方法。理论分析及数值模拟结果表明:相对于二阶形式的准P波方程,该TTI介质纯准P波一阶压力-速度方程在计算精度相同的情况下,具有更高的计算效率;同时能更准确地描述对称轴参数变化剧烈的TTI介质中准P波的传播规律。 相似文献
15.
16.
从双程声波方程出发,在交错网格空间中推导了地震波逆时延拓的高阶有限差分算子,依据最佳匹配层(PML)的方程分裂思路,得到了一阶声波方程的PML边界条件及其高阶差分格式,采用零时间成像条件和上行、下行波场互相关成像条件,实现了声波方程的叠后与叠前逆时深度偏移。逆时偏移对sigsbee_2b模型理论数据的偏移成像得到了满意效果。 相似文献
17.
地震数值模拟中选取的有限计算区域产生的边界反射会干扰正常波场模拟结果,因此引入人工边界条件来降低边界反射的影响。本文针对PML（perfectly matched layer）边界条件的解耦与非解耦差分形式,应用不同空间差分阶数进行地震波场数值模拟。空间差分阶数提高后,非解耦PML差分形式在计算效率和实现方式上均更具优势。针对CE（Clayton Engquist）边界条件受入射波入射角度限制、边界处精度低吸收效果不好等问题,将2阶CE边界条件和PML边界条件组合成一种新的组合边界条件,在保证吸收效果的同时减少衰减带厚度,从而达到提高计算效率的目的。数值模拟结果验证了算法的有效性。 相似文献
18.
波动方程数值模拟是研究地震波传播机理的重要工具,有限差分求解波动方程是当前地震波数值模拟的主要方法之一。当地下介质中的地震波速度较低或地震波高频成分丰富时,常规有限差分技术常常产生严重的数值频散误差,这种误差会降低数值模拟的精度,影响对地震波传播机理的分析。为压制地震波数值模拟时产生的数值频散误差,提高波场模拟精度,提出了基于NAD算子的时间四阶精度波动方程差分格式。根据对应的差分格式,分析了该差分格式的数值频散关系。与常规四阶精度差分算法的频散曲线相比,基于NAD时间四阶精度差分方法不但能够实现时间频散的有效压制,同时其基于更多网格点的位移分量和位移梯度分量空间微分求解方法还能够实现空间频散的有效压制。另外在相同模型条件下,基于NAD算法的声波方程时间四阶差分解法可采用大网格对模拟空间进行差分离散,减少网格数,提高计算效率。 相似文献
19.
三维频散介质中地质雷达信号的FDTD数值模拟 总被引：10，自引：0，他引：10
考虑频散介质的电磁波传播,引入随频率变化的电位移矢量D,并对电场强度E和电位移D进行标准化。开发出一种计算介电常数随频率变化满足Debye关系的频散介质中麦克斯韦方程的时间域有限差分（FDTD）解法,并编写了相应的计算原代码。在边界部分引入假象的介电常数和磁导率,实现了边界的无反射吸收,避免了Berenger完全匹配层中场分裂过程,提高了计算效率。计算实例可以看出频散介质中波的传播情况。 相似文献
20.
TTI介质qP波数值模拟方法因为考虑了倾角因素,可以比VTI介质qP波数值模拟方法更加准确地描述各向异性介质中地震波场的传播规律。文中用拟声波方程对TTI介质中的地震波场进行了高阶有限差分数值模拟,在改进衰减函数分布方式后,通过坐标变换,利用改进的完全匹配层(perfectly matched layer,PML)边界控制方程对波场边界进行吸收处理,取得了良好的效果;然后分析了拟声波方程数值模拟中的稳定性问题,并对波场中的伪横波进行压制。通过对不同模型的数值模拟,验证了文中使用的TTI介质拟声波波动方程的稳定性以及所采用的PML边界控制方程的可靠性和适用性。 相似文献