起伏地表条件下双变网格正演模拟.pdfVIP

·586· 中国地球物理2011 (17)油气田与煤田地球物理勘探 起伏地表条件下双变网格正演模拟 in finite conditionsdual-variabledifference free-surface Irregular modeling 黄建平3 谢金平1，+韩文功2李振春1 Xie Han LiZhenchun JinpingWengong HuangJianping 1．中国石油大学(华东)地球资源与信息学院青岛266555 2．中石化胜利油田分公司物探研究院 东营257000 1．引言 石油勘探正向复杂构造、复杂地表等地区发展，消除地表地形对地震资料的影响一直都是地震勘探中 的一个难以克服的问题。研究复杂近地表地震波传播问题数值模拟方法，对于认识复杂地表下地震波传播 的规律，弄清地震波在复杂地表条件下的传播特征，消除地表地形对地震资料的影响，以及研究基于起伏地 表的地震资料处理、采集和解释方法都有非常重要的意义【1|。 2．方法原理 在起伏地表条件下的地震波传播数值模拟方面，两个关键问题是如何选择适合起伏地表的数值计算方 法以及如何更好地实现自由边界条件。 通常起伏地表以下为低速介质，对于低速介质地震波模拟时，有限差分法需要较密的网格。而对于深 部速度较高的介质，则可以采用较疏的网格。因此可以根据介质的速度结构，在不同的计算区域，采用不 同的网格密度，这就是变网格技术。当空间网格变化时，时间稳定性仍然满足最短波长的原则，即时间采 样需要满足精细网格采样对应的稳定性条件。因此在变网格的基础上，又发展了局部变时间来配合变网格 技术，在全局时间层内，对需要精细划分的空间区域进行局部时间层精细采样，即双变网格高阶差分技术。 从弹性波的一阶速度一应力方程出发，其差分近似式如下式所示： ．．譬r- ． 、1 (1) O,g(x，z)一∑【．c2I-|g(x+△2H，z)+c2ig(x—As∥)J i-I 式中，Ci为待求的差分系数，△，是空间差分算子，它是变网格的步长出的函数。给定地质模型，在确定 网格剖分后就可求出各网格节点的差分算子供递推计算时调用。 合理选取自由边界是处理起伏地表的关键技术。实际地表是介质与空气之间的物性突变界面，称为自 由界面，其边界条件为垂直地表的应力为零，即： ％=0，f一工，z (2) 的广义虚像法的基础上实现的。不规则起伏表面可以分解成许多水平或竖直的表面单元和它们形成的顶角。 其特点是：自由边界通过切应力的采样位置，即在自由边界上直接处理剪应力。应力分量在水平和竖直两个 方向上做镜像处理。所谓虚像处理即自由边界处的应力在高阶差分算子的长度上～定是关于z一0反对称的。 3．结论 针对一系列起伏地表模型如水平界面、洼陷、隆起及光滑起伏地表等模型正演试算，结论是近地表速 度结构对深层反射波的影响很大，在单炮记录上呈现出复杂的同相轴，降低了信噪比：地表起伏可以引起 强散射体波和面波以及对地震波能量的调制和衰减作用，并造成反射波同相轴严重畸变、扭曲，甚至淹没， 在成像过程中可能产生构造假象。通过以上结论可知，正演起伏地表条件下的复杂模型对于认识复杂波场 特征及对后续研究有重大指导意义。 本课题由国家油气重大专项课题(2008ZX05014—001—008HZ)