复杂构造网格化及高精度地震波场谱元法数值模拟.pdfVIP

第 53卷第 4期 石 油 物 探 Vo1．53，No．4 2014年 7月 GE0PHYSICALPR0SPECTING FDRPETR0LEUM Ju1．，2014 文章编号：1000-1441(2014)04—0375—09 复杂构造网格化及高精度地震波场谱元法数值模拟 李洪建，韩立国，巩向博 (吉林大学地球探测科学与技术学院，吉林长春 13OO26) 摘要 ：深部复杂隐伏构造油气藏具有复杂的地质结构和地层条件，模拟并分析其地震波场的传播特征尤为重要。 研究并提出了复杂构造条件下的网格化及其高精度谱元数值模拟方法。首先基于Delaunay三角网格剖分方法 将模拟区域剖分为若干三角形网格，进而采用波前法将所有三角形网格重组 ，得到谱元法计算所需的四边形网 格。对比常规网格剖分方法 ，复杂构造网格剖分方法能够灵活适应速度边界的变化且网格均匀 ，质量较高。应 用基于复杂构造网格剖分的高精度谱元法对 Marmousi速度模型和二维随机介质模型进行数值模拟，并与常规 谱元法的模拟结果进行对比，可以看出高精度谱元法数值模拟结果同相轴连续，数值频散小，信噪比和分辨率较 高，取得了较好的应用效果。 关键词：谱元法；复杂构造；网格剖分；Delaunay算法 DOI：10．3969／j．issn．1000—1441．2014．04．001 ‘ ． 中图分类号：P631．4 文献标识码：A High precision spectralelementmethodbased0ngriddiscretization of complicated structureforseismicwavefieldnumericalsimulation LiHongjian，HanLiguo，GongXiangbo (CollegeofGeo-ExplorationScienceandTechnology，JilinUniversity，Changchun130026，China) Abstract：Forseismicexplorationofthedeep subtlestructuralreservoirswithcomplicatedgeologicalstructure andformationconditions，itisparticularlyimportantto simulateandanalyzetheirseismicwavefieldcharacteris— tics．Inthispaperwestudythegriddiscretizationforcomplicatedstructureandpropo seahighprecisionspectral elementmethod．Inordertogetthequadrilateralgridsusedforcalculatingspectralelementmethod，wedividethe targetareaintoseveraltriangulargridsbyDelaunaygriddiscretizationalgorithm，andthenallthetriangulargrids arereorganizedinto quadrilateralgridswithwavefrontmethod．Comparedwith conventionalgrid discretization method，thegriddiscretizationmethodforcomplicatedstructureisflexibleandadaptstothechang eofvelocity boundary，with fairly quality．NumericalsimulationofM armousivelocitymodeland2D random medium model show thattheres