基于叠前全波形反演储层流体识别技术.pdfVIP

·620· 中国地球物理2011 (18)储层地球物理 基于叠前全波形反演的储层流体识别技术 Reservoirfluiddiscriminationbasedon full—waveforminversion pre-stack 陈章+彭真明黎殿来 Chen DianlaiLi Zhang ZhenmingPeng 电子科技大学光电信息学院 成都610054 地震资料流体识别技术研究已经走向纵深发展的阶段，各种条件下的流体因子研究成为勘探地球物 理、油藏地球物理的热点。 现有的基于地震资料的储层流体识别主要基于叠前、叠后地震反演及属性分析技术。目前，地震反演 方法主要足基于褶积模型，存在其局限性。从描述地震波在地下传播规律的波动方程出发的全波形反演自 动考虑了地震波的反射、透射、绕射等行为，可从地震数据中直接反演出介质的弹性参数和密度；并充分 利用接收到的地震信号信息，可为基于地震资料的流体识别提供可靠的地球物理参数。因此，利用叠前全 波形反演直接进行储层流体识别近年来受到勘探地球物理界越来越多的学者的关注，成为新的研究热点。 本文主要介绍利用叠前全波形反演进行储层流体识别的关键技术，包括以下几个方面。 1．基于声波波动方程的模型正演 在声波近似的情况下，地下介质可用密度P和体积模量K完全描述，波动方程形式为： 厂 ， ．， 厂， ， 、] f南紊“沁l南grad肛刈)引(x,t) (1) 在众多的波场正演模拟方法中，有限差分法使用方便灵活，计算速度快、计算量小。本文采用高阶(时 间4阶空问10阶)交错网格有限差分法进行地震波场正演模拟。 在震源模拟上，选用Ricker子波作为震源子波，空间衰减形式采用指数形衰减方式，表达式为： h(x)=exp(一d2Itx-xoll2)，其中，X0为离散网格中震源中心位置，a为衰减系数。 数值频散是差分方程固有的本质特征，其程度直接决定地震波数值模拟效果。单纯依靠提高差分精度 并不能达到很好压制效果，所以本文采用通量校正传输技术(FCT)来压制数值频散现象，取得较好的效果。 有限差分另一个问题是用有界的计算区域来模拟地震波在无界区域中的传播，会使向区域外传播的波 在边界上产生反射，因而需要引入一定的人工边界条件来吸收人为边界的反射。本文采用完全匹配层(PML) 作为人工边界条件消除反射。 2．波场回传与叠前参数反演 本文采用的是最小二乘意义_F的时间域全波形反演。在反演方法的选择上，从全局最优的角度来说采 用随机搜索(如模拟退火、遗传算法等)的反演方法能够取得最好的效果，然而这些方法是以牺牲CPU 时间为代价的。而全波形反演呕待解决的重要问题就是计算效率问题，所以本文选用局部最优的共轭梯度 算法作为伞波形反演方法。 常大，所以本文采用Tarantola提出的波场回传技术，利用正演得到的波场与回传的剩余波场(野外观测记 录与合成记录的残差)的互相关构造梯度方向，避开了求解Frech∈t导数的过程。 要归一化。本文采用对数形式进行归一化，即令 (2) 足’X)=111(XolK(x))，P’(x)=111(po／p(x)) Ko、P。为任意的体积模量和密度值，模型参数为：J，l=【置’p]1。 此时，声波方程形式变为： 中国地球物理2011 ·62l· [掣簧“V[掣鲫d肛∽叫川 模型参数的扰动变为： 训小(删硝小一孚掣莓』d，鲥p(引；对鲥出，，；薯) (4a) 种(小(删硝小一掣掣军p吲刈蚓帅，；气)