我卡实现反向时间偏移技术,PDE(偏微分方程)是一个实现声波方程,所使用的数值方法是一个FD(有限差分法)使用SBP分部(求和)属性,以确保稳定域,即。这变化系数,FD用附近的边界和边界点,除了坐(同时近似计算),模拟连续能量法是如何实现的,这将保证长时间稳定模型。我需要帮助与实现坐2 d网格,每次运行的模型我实现它发散坐在最后一个时间步长。有什么著名acoustic-2D-wave方程弱强制边界条件,包括坐吗? ?
1回答
\ begingroup美元
\ endgroup美元
2
与SBP-SAT方法我没有经验,但是你可能尝试了和通过引用(和潜在联系的作者)https://arxiv.org/pdf/1802.06123.pdf。
就我个人而言,我不认为严格需要RTM使用这种复杂的造型方法。在理想情况下,您有deghosted(即。向上移动的)记录数据和平滑的速度模型。然后,使用标准的低阶有限差分模型在一个统一的网格应该已经提供相当不错的结果相当小的代价。可能,你会取得更好的结果通过优化数据和速度模型比改善数值代码之外的标准。这个简单的MATLAB代码https://csim.kaust.edu.sa/files/ErSE328.2013/LAB/WE_LABS/index.html已经做了杰出的工作使用二阶精度在空间时间和八阶精度,使用一个简单的FD方法。
当然,如果你所做的最先进的(学术)研究或处理一个非常典型的速度模型与一个极端的动态范围,这些都是有趣的话题来调查,我希望有人可以帮助你……但我希望我的回答可以帮助你。
-
\ begingroup美元 埃里克,我真的很想谢谢你接触和回答这个问题,我有很多困难在表达我做的工作和寻找可以联系的人。然而,我完全理解你的观点,这个框架的主要边缘(SBP-SAT)增加了RTM multi-grid技术,我们可以使用它。每一块可以有自己的空间和时间间隔不破坏整体的稳定域和数值方法。这可以帮助很多复杂的速度域可以实现巨大的增长在保持精度计算方面。 \ endgroup美元- - - - - -卡里姆Morad2019年5月20日,在球队的
-
\ begingroup美元 卡里姆,那么我认为你最好联系相关论文的作者如上图。其他作者必须经历过最典型的错误可以体验当建模算法。如果你幸运的话,他们可能有一些共享代码. . ! \ endgroup美元- - - - - -埃里克2019年5月21日9:53