共查询到17条相似文献，搜索用时 125 毫秒
1.
曲线坐标系下的完全匹配层吸收边界条件 总被引：1，自引：0，他引：1
在地震波数值模拟中,需要采用吸收边界条件以吸收人为边界反射。本文针对曲线坐标系下的二阶弹性波方程提出了一种完全匹配层(PML)吸收边界条件。与直角坐标系下的PML吸收边界条件类似,曲线坐标系下的PML吸收边界条件是一种在频率域中给出的人工边界条件,由相应的复坐标变换得到。在变换到时间域后,完全匹配层中将出现复杂的卷积运算。为了避免这些卷积运算,引入了4个中间变量。为了简化自由边界条件,采用正交贴体网格对起伏地表模型进行网格剖分。数值算例表明,该方法可以有效消除人为边界反射。 相似文献
2.
地震数值模拟中,完全匹配层(PML)边界能有效地吸收衰减地震波,得到无边界反射干扰的波场快照和地震记录。在前人研究的基础上,进一步分析了分裂完全匹配层(SPML)的衰减机制,通过波动方程模拟,分析了震源主频、空间网格间距、介质速度等参数对SPML边界吸收衰减特征的影响。得到了不同条件下,PML边界对地震波的吸收效果。通过分析对比,得出了震源主频对PML吸收效果无直接影响,空间网格间距与PML吸收效果成反比,高速层PML吸收效果缓慢等结论。并通过Marmousi2模型测算,对结论进行了验证。最后给出了复杂模型、不同尺度数值模拟中SPML参数的选取方法,为地震数值模拟中SPML边界参数的定量选取奠定了基础。 相似文献
3.
从一阶速度—应力弹性波方程出发,基于旋转交错网格,推导了时间二阶精度空间2M阶精度的有限差分离散格式。阐述了递归卷积复频移完全匹配层(CPML)边界条件的原理,建立了一阶速度—应力弹性波高阶差分CPML边界条件的递推公式。开展了CPML边界中关键参数m、κ和α的选取实验,通过分析反射误差分布图,选取了CPML边界条件中最优参数。全局反射误差与波场快照都说明,CPML较PML对隐失波具有更优的吸收性能。基于Matlab平台,编写了基于CPML边界的旋转交错网格弹性波正演模拟程序,应用该程序对各向异性介质及随机介质进行了模拟,得到了弹性波正演剖面记录及波场快照,通过对正演剖面记录及波场快照的分析,可以更清楚地了解弹性波在各向异性介质及随机介质的传播特性,指导非均匀介质中地震勘探资料解释。 相似文献
4.
为了更好地研究近岸海域波生沿岸流,建立了基于高阶Boussinesq水波方程的波生沿岸流时域数值模型。控制方程在中等水深范围内具有较好的色散性和变浅作用性能,同时具有二阶完全非线性特征,适合描述近岸区域波浪强非线性运动。通过采用松弛造波方法实现了非线性波浪的无反射入射,采用周期性侧边界条件模拟开敞边界。通过数值试验,讨论了模型中主要参数对数值结果的影响。利用率定后的参数模拟了均匀坡度海岸上产生的沿岸流,通过和实验数据的对比验证了模型的准确性和适用性。利用模型数值模拟了不同波浪入射条件(包括周期、波高和波浪入射角度)对波生沿岸流的影响。 相似文献
5.
通过对声波方程采用二阶有限差分与四阶有限差分近似,进行数值模拟及波场分析,边界条件利用吸收边界和增加衰减带混合使用的方法。计算结果表明,四阶有限差分的精度明显高于二阶有限差分,并且边界吸收效果好。对于单炮记录难于模拟复杂模型的情况下,实现了单次叠加剖面和爆炸反射面的零炮检距剖面,并进行了偏移处理。地震剖面更加清晰地反映了复杂构造的信息,更好地实现了对复杂地质模型的正演模拟。 相似文献
6.
完全匹配层吸收边界（PML）已经被证明是非常有效的边界吸收技术,对体波和面波的吸收都具有非常好的效果,已经被广泛应用于弹性波的数值模拟中。但是在某些情况下传统的PML技术还是存在一定的问题,比如对掠射情况下的体波和窄区域自由表面条件下的面波的吸收等等。在坐标变换中采用复频移拉伸函数的复频移PML可以有效地改善PML边界条件的吸收性能。基于弹性波一阶速度-应力方程,推导了复频移PML的递推卷积实现方法,并采用交错网格高阶有限差分法对其进行了数值模拟,与传统的PML进行了对比。结果表明：传统的PML对掠射情况下的体波和窄区域自由表面条件下的面波吸收不足,会产生虚假反射,影响真实波场;而基于递推卷积的复频移PML算法能够有效地改善困难情况下的吸收效果,并且在实现过程中不用分裂变量,应用更加方便简单。计算卷积时采用递推的形式,推导过程直观易懂,易于编程,而且不会增加计算量,存储量也没有太大的变化。 相似文献
7.
采用颗粒离散单元法进行动力计算时,人工截断边界上需设置吸收边界条件,以防止波的反射。鉴于颗粒离散单元数值计算模型的人工边界上颗粒单元半径大小不一、边界面凸凹不平,在连续介质的黏性、黏弹性、自由场边界条件方程基础之上,推导出适用于离散介质的等效方程。在离散介质的黏性边界条件等效方程中引入微调系数,提出比值迭代法以快速确定其最优值,以实现对波的最佳吸收。采用二维颗粒离散单元计算软件PFC2D,分别建立黏性、黏弹性、自由场边界条件相关数值分析模型,探讨颗粒分布模式对黏性边界上颗粒单元半径、速度分布及比值迭代过程的影响;采用外源波动算例及经典Lamb问题算例验证黏弹性边界设置方法的正确性;通过隧洞算例检验提出的自由场边界条件设置方法的正确性。 相似文献
8.
地震数值模拟中选取的有限计算区域产生的边界反射会干扰正常波场模拟结果,因此引入人工边界条件来降低边界反射的影响。本文针对PML（perfectly matched layer）边界条件的解耦与非解耦差分形式,应用不同空间差分阶数进行地震波场数值模拟。空间差分阶数提高后,非解耦PML差分形式在计算效率和实现方式上均更具优势。针对CE（Clayton Engquist）边界条件受入射波入射角度限制、边界处精度低吸收效果不好等问题,将2阶CE边界条件和PML边界条件组合成一种新的组合边界条件,在保证吸收效果的同时减少衰减带厚度,从而达到提高计算效率的目的。数值模拟结果验证了算法的有效性。 相似文献
9.
弹性波正演模拟中PML吸收边界条件的改进 总被引：4，自引：0，他引：4
在弹性波有限差分正演模拟中, 完全匹配层(PML) 吸收边界条件是使用广泛、吸收效果最好的吸收边界条件.在目前的两种PML实现方法中, 分裂形式的完全匹配层(SPML) 方法计算存储量大、编程实现复杂; 非分裂形式的完全匹配层(NPML) 方法计算效率低、计算过程复杂.针对传统PML吸收边界条件在实现过程中存在的问题, 推导出了一种简洁有效的非卷积实现的NPML吸收边界条件, 既不需要对场分量进行分裂, 也不需要做复杂的卷积运算.分析结果表明, 本文实现的NPML吸收边界条件不仅具有良好的吸收衰减性能, 而且计算方程简单, 编程实现容易, 占有内存更小. 相似文献
10.
TTI介质qP波数值模拟方法因为考虑了倾角因素,可以比VTI介质qP波数值模拟方法更加准确地描述各向异性介质中地震波场的传播规律。文中用拟声波方程对TTI介质中的地震波场进行了高阶有限差分数值模拟,在改进衰减函数分布方式后,通过坐标变换,利用改进的完全匹配层(perfectly matched layer,PML)边界控制方程对波场边界进行吸收处理,取得了良好的效果;然后分析了拟声波方程数值模拟中的稳定性问题,并对波场中的伪横波进行压制。通过对不同模型的数值模拟,验证了文中使用的TTI介质拟声波波动方程的稳定性以及所采用的PML边界控制方程的可靠性和适用性。 相似文献
11.
研究了二阶声波方程频域PML边界条件和频域变网格并行计算技术。PML边界是一种较为理想的吸收边界方法,多用在求解时域应力速度方程中,但对于频域声波正演,二阶位移方程更常用。从一阶声波方程PML吸收边界条件导出频域二阶位移方程PML边界条件,模拟算例得到的频率切片、时间切片和地震记录对比都说明该边界条件吸收效果很好。频域单炮正演不同频率间是独立的,据此低频部分采用大网格计算,高频采用小网格,实现变网格步长计算技术,这是较时间域正演的一个优势,在保证模拟质量的同时,减少计算量和内存消耗。 相似文献
12.
针对起伏地表条件下的地震波数值模拟问题,提出将起伏的地表映射到一个规则的长方形网格坐标系中,并在此基础上推导出变换域中的波动方程。根据导出来的波动方程应用有限余弦变换有限差分方法进行地震波场的数值模拟,以解决起伏地表条件下的数值模拟。 相似文献
13.
14.
影响完全匹配层方法吸收效果的主要因素有吸收层厚度和衰减系数,笔者通过地震波场数值模拟,讨论了吸收层厚度和衰减系数对吸收效果的影响。研究表明,在衰减系数一定时,吸收层厚度越大,边界反射越弱。吸收层的厚度一般在12～20道较为适当,在吸收层的厚度比较小时,随着衰减系数的增大,边界反射逐渐减弱甚至完全没有边界反射,但是若再进一步增大衰减系数则又会逐步出现边界反射。因此,在实际应用中应该注意衰减系数的选择,衰减系数一般在500～2 000之间,衰减系数较大反而影响其效果。 相似文献
15.
16.
17.
反转构造模型的地震波模拟及成像 总被引：1，自引：0，他引：1
为了检验由地震剖面解释得出的反转构造的可靠性,探讨反转构造地区在地震数据采集,处理中应采用的有效方法,设计了一个典型的反转构造模型－突起构造。以该模型为基础,针对复杂构造地震波模拟中常见起源频散严重、人为边界反射吸收不理想两个问题,应用高阶差分法以及一种新的边界吸收条件,对反转模型进行了地震波场数值模拟和深度偏移成像处理。采用中心放炮两边接收的观测方式,共模拟了185炮地震记录。研究表明,波动方程高阶差分法是模拟复杂构造模型中地震波传播的有效方法,除了近于直立的断层因难以接收到其反射波而较难成像外,其余的反转断层及其组合在理论上是可以通过合适方法（如Fourier有限差分深度偏移）较好地成像的。 相似文献