基于NAD算子和FCT技术改进的Runge-Kutta方法及其波场模拟.pdfVIP

２０１６ ３１ ３ 　 Ｐｒｏｇｒｅｓｓ ｉｎ Ｇｅｏｐｈｙｓｉｃｓ ＩＳＳＮ １００４２９０３ 年第卷第期 地球物理学进展 ２０１６，３１（３）：１２７４１２８０ ｈｔｔｐ／／ ｗｗｗ．ｐｒｏｇｅｏｐｈｙｓ．ｃｎ ＣＮ １１２９８２／ Ｐ ．２０１６． ＮＡＤ ＦＣＴ ＲｕｎｇｅＫｕｔｔａ ． ３１ ３ １２７４１２８０ ｄｏｉ １０． ６０３８／ 张朝元 基于 算子和 技术改进的 方法及其波场模拟地球物理学进展， （）： ， ： ｐ ＺＨＡＮＧ Ｃｈａｏｙｕａｎ． ２０１６． ＩｍｐｒｏｖｅｄＲｕｎｇｅＫｕｔｔａｍｅｔｈｏｄａｎｄｉｔｓｗａｖｅｆｉｅｌｄｓｉｍｕｌａｔｉｏｎｂａｓｅｄｏｎＮＡＤｏｐｅｒａｔｏｒａｎｄＦＣＴｔｅｃｈｎｉｑｕｅ．Ｐｒｏｇｒｅｓｓｉｎ Ｇｅｏｐｈｙｓｉｃｓ ｉｎ Ｃｈｉｎｅｓｅ ３１ ３ １２７４１２８０ ｄｏｉ １０．６０３８／ ｐ （ ）， （）： ， ： 基于ＮＡＤ算子和ＦＣＴ技术改进的ＲｕｎｇｅＫｕｔｔａ方法及其 波场模拟 Ｉｍｐｒｏｖｅｄ ＲｕｎｇｅＫｕｔｔａ ｍｅｔｈｏｄ ａｎｄ ｉｔｓ ｗａｖｅｆｉｅｌｄ ｓｉｍｕｌａｔｉｏｎ ｂａｓｅｄ ｏｎ ＮＡＤ ｏｐｅｒａｔｏｒ ａｎｄ ＦＣＴ ｔｅｃｈｎｉｑｕｅ 张朝元 ＺＨＡＮＧ Ｃｈａｏｙｕａｎ 　 ６７１００３ 大理大学数学与计算机学院，大理 ， ， ， Ｃｏｌｌｅｇｅ ｏｆ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｃｏｍｐｕｔｅｒ Ｄａｌｉ Ｕｎｉｖｅｒｓｉｔｙ Ｄａｌｉ６７１００３ Ｃｈｉｎａ 摘　 要　 地震波正演数值模拟中所产生的数值频散严重地影 Ａｂｓｔｒａｃｔ　 Ｔｈｅ ｎｕｍｅｒｉｃａｌ ｄｉｓｐｅｒｓｉｏｎ，ｅｘｉｓｔｉｎｇ ｉｎ ｔｈｅ ｆｏｒｗａｒｄ ｎｕｍｅｒｉｃａｌ ｓｉｍｕｌａｔｉｏｎ ｏｆ ｓｅｉｓｍｉｃｗａｖｅ ｉｎｔｅｒｆｅｒｅｓｗｉｔｈ ｔｈｅｗａｖｅｆｉｅｌｄ 响波场模拟的精度和分辨率基于声波方程，以阶精度的近 ， ． ８ ｍｏｄｅｌｉｎｇｓｅｒｉｏｕｓｌｙａｎｄｄｅｃｒｅａｓｅｓｔｈｅｐｒｅｃｉｓｉｏｎｏｆ ｓｉｍｕｌａｔｉｏｎｇｒｅａｔｌｙ． 似解析离散算子来校正差分算子对空间偏导数进行离散，以三 Ｂａｓｅｄ ｏｎ ｔｈｅ ａｃｏｕｓｔｉｃ ｗａｖｅ ｅｑｕａｔｉｏｎ，ａｎ ｉｍｐｒｏｖｅｄ ＲｕｎｇｅＫｕｔｔａ 阶 方法对时间导数进行离散，并结合通量校正传 ＲｕｎｇｅＫｕｔｔａ ｍｅｔｈｏｄ ｉｓ ｐｒｏｐｏｓｅｄ ｂｙ ｃｏｍｂｉｎｉｎｇ ＮＡＤ ｏｐｅｒａｔｏｒ ａｎｄ ＦＣＴ 输技术以抑制数值频散，从而建立基于 算子和 技术 ｔｅｃｈｎｉｑｕｅ． Ｔｈｅ ｎｅｗ ｍｅｔｈｏｄ ｕｓｅｓ ａｎ ｅｉｇｈｔｈｏｒｄｅｒ ａｃｃｕｒａｔｅ ｎｅａｒｌｙ ＮＡＤ ＴＣＦ ａｎａｌｙｔｉｃ ｄｉｓｃｒｅｔｅ ｏｐｅｒａｔｏｒ ｔｏ ｄｉｓｃｒｅｔｉｚｅ ｔｈｅ ｈｉｇｈｏｒｄｅｒ ｓｐａｔｉａｌ 改进的ＲｕｎｇｅＫｕｔｔａ方法．波场模拟结果表明，该方法在抑制数