地震道非线性约束反演.pdfVIP

地震道的非线性约束反演 姜 亮l 黄捍东2 魏修成2 贺振华3 (1清华大学；2石油大学(北京)；3成都理工大学) 摘要对于复杂波场储层研究，非线性反演在解空间的性质、状态方面比线性反演更具 有优势。基于非线性最优化理论本文提出了构造一测井约束下的地震资料非线性反演方法，该 方法综合了测井约束反演和宽带约束反演等线性反演方法的优点。实际地震资料的处理结果表 明，该方法运算速度快、稳定、精度高，在油气储层研究中具有良好的应用前景。 1 引言 在复杂波场油气储层研究中，波阻抗是识别地下介质岩性最有效的参数之一。20世纪 80年代末，合成声波测井、宽带约束反演、波动方程反演等方法有了较大的发展，成为我 们进行目标层段储层预测的一个重要手段。 但是，随着我们对地层岩性尖灭、透镜体和多断层小断块等复杂油气储层的研究发现， 基于线性反演方法的地震约束反演已不能适应复杂波场储层的进一步研究，这是因为地震道 的反射应看作一个非线性系统。对于这样的非线性问题，采用相应的非线性反演，其解空间 的性质、状态将优于线性反演。在本文的研究中，我们吸收了宽带约束反演、测井约束反演 和模型法反演等线性反演的优势，在进行非线性反演运算的同时，采用了构造一测井共同约 束下的井外推波阻抗反演。实际资料处理结果表明，该方法可以得到与实际测井资料对比的 波阻抗剖面，反演出的地震参数具有较高的可信度和明确的地质意义，可用于油气勘探开发 中的储层或油藏预测。 2方法原理 2．1地震道线性反演方法 目前地震道线性反演的原理基本相同，即从褶积模型出发，利用拉格朗日乘子法解以下 线性方程组： s(t)=乏叫(t—r)r(r) (1) 求出反射系数系列r，，其中S(t)为地震道，叫(t—r)为具有不同延迟的子波。 r=，．o+(硼T删+ai)一1训T(S—So) (2) 式中，A为阻尼因子；ro为反射系数向量的初始猜测，so=叫·，．o。 为了得到波阻抗还需将反射系数序列用递归法进行转换： 2 ㈥¨) ri乏2z1掣llriJ Zi+I=篝铲z。星[毒] 式中，z，为给定的地表波阻抗值。 在上述反演过程中反射系数序列与地震记录之间的关系是线性的，因而该方法称为褶积 ·244· 模型的线性反演。 2．2地震道非线性反演方法 假设反射系数序列模型r可表示为： ，．(t)=∑≥生■劫[t—iAt]i=1，2，…，L (4) 扣lZi+1十Zi 式中，△￡为采样间隔。 地震道可表示为： 一 (5) s(￡)=叫·．三型坤[￡一iAt]+挖(t) {=lZl+1十Zi 式中，地震道与波阻抗的关系是非线性的，因而称为非线性反演。 在解上述非线性问题时，人们为便于求解常采用线性反演方法来求解。这样就大大地降 低了反演的精度。基于非线性最优化理论，我们提出了地震道非线性最优化反演的思想，其 目标函数为： s—d II—mJn (6) f(z)=lI 其中，z为波阻抗；s为模型响应由式(5)得到；d为实际地震记录。 利用泰勒级数将式(8)在给定的波阻抗初始值zo附近展开得：