波动方程正演在海拉尔盆地的应用.pdfVIP

波动方程正演在海拉尔盆地的应用 吴清岭1李来林1陈志德1毕民1董良国2 l(大庆油田有限责任公司勘探开发研究院，163712) 2(同济大学) 序言 已知地下地质模型，那么它的地震响应是怎样的，用计算机把结果计算出来，这就是地 球物理正演模型的研究内容。因此，利用正演模型，地震解人员可以方便地在计算机上验证 自自己的解释成果，研究地下地层参数情况、构造情况。因此可以直接为找油、气勘探服务。 随着计算机技术的发展和叠前深度偏移技术在实际资料处理中的应用，研究其理论模型正演 记录的叠前偏移效果，具有实际的知道意义。利用多次覆盖采集技术，虽然提高了信噪比， 但也降低了分辨率，耗损了许多有用的信息。我们可以通过正演模型，让其最佳接近实际地 震资料的野外采集，这样可以弥补对时间剖面的认识不足。因此正演模型是地震勘探中一项 重要的工作，也是地震反演和叠前深度偏移的基础。由于模拟计算震源和检波器的位置可任 意设置，因此，正演模型还可以用于勘探方面其它研究工作。 在计算机上进行地震波正演数值模拟，国内八十年代开始进行大量研究。由于当时计算 机条件的限制，只能进行小规模的理论研究。地震数值模拟按照原理的不同，可以分为射线 方法和波动方法。按照几何地震学原理，进行射线追踪，来模拟地震波的传播，就是所谓的 射线方法。这种方法参数输入比较复杂，计算的地震响应不能保持地震波的动力学特征，特 别是复杂地质构造的射线追踪，会出现射线的盲区，但是计算速度快，所得到的地震波到达 时间比较准确。按照波动理论，通过求解波动方程来模拟地震波的传播过程，就是波动方法。 根据地下介质的不同，可以采用不同的波动方程进行近似。目前声学介质的声波方程可以模 拟野外的炸药震源激发，地震纵波检波器接收，与实际地震记录有很好的相似性，因此在实 际工作中得到了广泛的应用，得到令人满意的效果。求解波动方程的方法有多种，目前主要 富立叶方法、有限元法和有限差分法。有限差分法是一种快速有效的波动方程数值计算方法， 它可以对任意复杂地质构造，进行数值模拟，得到地质构造的正演单炮记录和自激自收时间 剖面。对于—个给定大小的地质模型，一旦差分网格大小确定，运算时间和模型的复杂程度 无关。 理论与方法 变速变密度声波方程 在二维情况下，地震波的传播可以用变密度声波方程可以用下列方程表示： 害一c2【(睾+》一吉(芸‘芸+鲁+》】“(X--XO,Z--,K0)．，(，) 这里c=如z)表示地震波的传播速度：H=M伍二f)表示声波压力；p=P(葺z)表示地下 介质的密度，6(X-Xo，z—zo)·，(f)表示在地下坐标为(勘zD)处有—个点震源贝f)。 边界条件和初始条件 自由表面边界上应力为零 102 仃琵Iz卸=0 两种介质的分界面上，法向应力连续 仃鬣‘1’(工，z)I铲^=仃杞‘2’(工，z)I：曲 两种介质的分界面上，法向位移连续 l印，，Jz，^2砧f2JI：；^ 初始位移为零 引 ∥O卸O 初始速度为零 = 塑孤 ‰蚓 I|^ 数值计算公式 在x,z方向上等间隔进行离散采样，对求解区域剖分，设X方向上节点编号为m，z方向 上的节点编号为n，t方向上的节点编号为j，空间采样间隔为／NX=AZ=h，时间采样间隔为dt。 对于二阶偏导数，采用二阶中心差分格式；对于一阶偏导数，采用一阶中心差分格式，因此 波动方程的有限差分数值计算公式为 u(m，n，J+1)=一u(m，，l，七一1)+(2-4p)u(m，n，jf) 州l-也掣州斛咖册 +p(1一兰堕!!芝掣)“(m，n+l，jf) +p(1-丛掣Mm-l,n,D 州l-丛铲Mm,n-1,j『) +(5(m-mO，万-nO)·f(J·出) 这里，p=兰器三，mD，一。是用网格节点表示的震源坐标。 差分公式的稳定条件dt、h／c 为了避免计算时出现网格频散现象，要求