基于分部求和(Summation By Parts)方法和同时逼近项(Simultaneous Approximation Terms)技术建立的有限差分方法,具有更高的精度和稳定性。同时在介质几何不连续、参数突变条件具有较大的优势。国内对SBP-SAT方法的相关研究目前较少,...基于分部求和(Summation By Parts)方法和同时逼近项(Simultaneous Approximation Terms)技术建立的有限差分方法,具有更高的精度和稳定性。同时在介质几何不连续、参数突变条件具有较大的优势。国内对SBP-SAT方法的相关研究目前较少,论文对该方法的研究背景,方法发展过程进行了介绍并基于SBP-SAT方法和弹性波动理论,结合初边值条件,推导出曲线网格条件下的弹性波动SBP-SAT离散方程。最后,通过数值模拟实现地震波传播过程,介绍该方法在地震数值模拟领域中的应用价值和前景。展开更多
文摘基于分部求和(Summation By Parts)方法和同时逼近项(Simultaneous Approximation Terms)技术建立的有限差分方法,具有更高的精度和稳定性。同时在介质几何不连续、参数突变条件具有较大的优势。国内对SBP-SAT方法的相关研究目前较少,论文对该方法的研究背景,方法发展过程进行了介绍并基于SBP-SAT方法和弹性波动理论,结合初边值条件,推导出曲线网格条件下的弹性波动SBP-SAT离散方程。最后,通过数值模拟实现地震波传播过程,介绍该方法在地震数值模拟领域中的应用价值和前景。
文摘基于速度-应力形式的弹性波动方程,采用分部求和和同时逼近项技术建立的SBP-SAT方法,推导了横向各向同性(transversely isotropy,TI)介质的矩阵对称型(symmetric matrix form,SMF)弹性波动方程离散形式,并通过能量法进行了稳定性分析。将该方法应用于倾斜横向各向同性(tilted transverse isotropic,TTI)介质模型、垂直横向各向同性(transverse isotropy with a vertical axis of symmetry,VTI)介质和含裂缝及曲线域的复杂介质模型,对所得的速度幅值和单炮记录分析并总结规律;对不同时间步长、单元网格数的结果进行对比,得出计算效率并验证该方法在求解P-SV波传播问题上的正确性。数值模拟结果表明,该方法模拟精度高,适用性好,在地震数值模拟领域有很好的应用价值和前景。
