基于加权平均导数的频率-空间域正演模拟及其GPU实现.pdfVIP

第 卷 第 期 地 球 物 理 学 报 ， ５８ ４ Ｖｏｌ．５８ Ｎｏ．４ 　 　 　 　 　 年 月 ， ２０１５ ４ ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＧＥＯＰＨＹＳＩＣＳ Ａ ｒ．２０１５ 　 　 　 ｐ 唐祥德，刘洪，张衡 基于加权平均导数的频率 空间域正演模拟及 实现 地球物理学报， （）： ， ： ．２０１５．  ＧＰＵ ． ５８４ １３４１１３５４ｄｏｉ１０． ／ ６０３８ ｊｇ ， ， Ｔａｎ ＸＤ ＬｉｕＨ Ｚｈａｎ Ｈ．２０１５．Ｆｒｅｕｅｎｃｓａｃｅｄｏｍａｉｎａｃｏｕｓｔｉｃｍｏｄｅｌｉｎ ｂａｓｅｄｏｎａｎａｖｅｒａｅｄｅｒｉｖａｔｉｖｅａｎｄＧＰＵ ｇ ｇ ｑ ｙ ｐ ｇ ｇ （ ）， （）： ， ： ／ ｉｍ ｌｅｍｅｎｔａｔｉｏｎ．犆犺犻狀犲狊犲犑．犌犲狅犺狊． ｉｎＣｈｉｎｅｓｅ ５８４ １３４１１３５４ｄｏｉ１０．６０３８ ｐ 狆狔 ｊｇ 基于加权平均导数的频率 空间域正演模拟  及犌犘犝实现 ， １ ２ １ ３ 唐祥德 ，刘洪 ，张衡 １中国科学院地质与地球物理研究所，中国科学院油气资源研究重点实验室，北京　１０００２９ ２中国科学院大学，北京 １０００４９ 　 ３国土资源部海底矿产资源重点实验室，广州海洋地质调查局，广州　５１００７５ 摘要 传统基于旋转坐标系的频率 空间域正演模拟方法仅适用于方形网格，而实际生产中矩形网格广泛存在，本  　 文提出一种适用性广的正演差分算子，不仅适用于方形网格而且适用于矩形网格 通过综合运用平均导数法、加速 ． 项加权平均、模拟退火法压制频散和减少单个波长所需网格点数，从而提高算法精度和减少计算量 在该方法的基 ． 础上采用不完全 ＬＵ分解作为求解 Ｈｅｌｍｈｏｌｔｚ方程的预条件，并利用图形处理器加速计算速度，很大程度上提高了 频率域正演的效率． 关键词 频率 空间域正演；矩形网格；不完全 分解； 加速  ＬＵ ＧＰＵ 　 ： ／ 中图分类