CSC频率—空间域波动方程数值模拟

详细信息   

• 责任者：吕晓春^c中国地质大学(武汉)构造与油气资源教育部重点实验室^d1985^e博士研究生 ; 顾汉明 ; 成景旺 ; 周丽
• 刊名：石油地球物理勘探
• 出版年：2014
• 卷期：49^b2
• 页码：288～294
• 索书号：P613.06/167-2
• 图表：^c参8插图
• 文摘期号：2014/09

摘要

针对频率‐空间域波动方程数值模拟需要巨大内存空间的现状，提出了利用列索引压缩存储(CSC)技术存储大型稀疏非对称复数型的矩阵系数。CSC技术将系数矩阵转化为三个一维数组来存储，分别存储系数非零元素、非零元素对应所在的行以及每列起始非零元素所在位置。经CSC技术压缩存储后显著减少了内存空间及计算量，在计算时只有少许的非零元素参加计算，且根据三个一维数组可以简便地找到对应的非零元素，进而采用LU分解快速而精确地求解。该文基于Jo等提出的最优化9点差分方法，应用CSC技术在频率—空间域进行二维声波方程数值模拟。通过对Corner‐edge模型和二维Marmousi模型进行试算，可以显著节省内存需求，明显提高计算速度，进而得到精度较高的正演结果。