基于紧凑存储的任意四边形网格有限元法地震波场数值模拟.pdfVIP

第24卷 第 5 期 CT理论与应用研究 Vol.24, No.5 2015年9月（667-680） CT Theory and Applications Sep., 2015 刘瑞合, 周建科, 周学锋, 等. 基于紧凑存储的任意四边形网格有限元法地震波场数值模拟[J]. CT 理论与应用研究, 2015, 24(5): 667-680. doi:10.15953/j.1004-4140.2015.24.05.04. Liu RH, Zhou JK, Zhou XF, et al. Finite element method with arbitrary quadrilateral meshes for numerical modeling of seismic wave based on compact storage[J]. CT Theory and Applications, 2015, 24(5): 667-680. (in Chinese). doi:10.15953/j.1004-4140.2015.24.05.04. 基于紧凑存储的任意四边形网格 有限元法地震波场数值模拟 1,2 3 3 1 1 刘瑞合 ，周建科 ，周学锋 ，印兴耀 ，曹丹平 1.中国石油大学（华东）地球科学与技术学院，山东 青岛266580 2.胜利油田分公司油气勘探管理中心，山东 东营 257000 3.中海石油（中国）有限公司天津分公司渤海石油研究院，天津 300452 摘要：采用有限元法精确模拟复杂介质条件下的地震波场，一般采用三角单元，但在节点数相 同的情况下，三角单元的计算精度不如矩形单元高，采用三角单元模拟复杂界面时，编制确定 结构刚度矩阵非零元素位置的程序也较为麻烦。采用矩形单元离散含有倾斜或起伏界面的地质 模型时，无法避免绕射噪声，若加密网格又会增加计算量。为此，本文基于任意四边形单元模 拟声波的传播，在倾斜或起伏界面条件下，可以有效避免因离散引起的“阶梯状”界面，在不 增加计算量以及内存占用的前提下，有效地消除离散绕射噪声。采用对角的集中质量矩阵代替 一致质量矩阵，避免矩阵的求逆运算，从而提高显式有限元法的计算效率；对结构刚度矩阵采 用紧凑存储格式，每一行需要存储的元素最多为 5 个，同时零元素不参与运算, 既减少内存的 占用，又极大地提高计算效率。 关键词：有限元法；数值模拟；集中质量矩阵 doi:10.15953/j.1004-4140.2015.24.05.04 中图分类号：O242.21，P315 文献标志码：A 波动方程数值模拟对认识地震波传播规律，进行地震属性研究，地震资料地质解释以 及储层评价等具有重要的理论和实际意义[1－2]。目前有限元法是最精确的一种正演方法，具 有可模拟任意复杂结构和易于进行边界处理等优点[3－4]。采用有限元法精确模拟起伏地形、 复杂构造和复杂介质条件下的地震波场，以前一般采用三角网格单元[5－8]。虽然三角单元具 有很好的“适应性”，几乎任何复杂边界的弹性体总可以划分为三角形，并且三角单元计算 公式简单，但在相同节点的情况下，矩形单元的计算精度比三角单元高许多[4]。采用常规的 矩形单元离散模型，不能适应曲线边界和非正交的直线边界，而是以“阶梯状”的界面来 替代实际的界面，其结果是产生了非物理的绕射波，干扰了对有效波场的识别，降低了数 值模拟的精度。 为此，Ke 等[9]联合使用三角单元和矩形单元求解二维起伏地表条件下标量波动方程的 响应，即采用三角单元拟合地形的起伏，而在其他区域采用矩形单元。此外，近年来也尝 收稿日期：2015-06-26。 基金项目：国家重点基础研究发展计划“973”项目（2013CB228604）；山东省自然科学基金（ZR2014DM009）。 668 CT理论与应用研究 24卷 试把有限元法与有限差分法的优