高分辨率过井地震层析成像的级联算法及其应用.pdfVIP

第34卷第6期 物 探 与 化 探 Vol_34．No．6 2010年 12月 GE0PHYSICAL＆ GE0CHEMICALEXPLORAT10N Dec．．2010 高分辨率过井地震层析成像的级联算法及其应用 瞿辰 ，杨文采，于常青 (中国地质科学院 地质研究所，北京 100037) 摘 要：将过井地震层析成像的级连算法应用到垦71区井间地震资料处理中，得到了高分辨率、精确的井问纵波速 度结构。首先根据初至出射角，旋转合成得到清晰的纵波波场，再精确拾取初至到时。初始速度模型是参考测井 资料的分层模型，第一级反演算法为分辨率较低但稳定的BPT反演算法，为第二级高分辨率反演算法——阻尼 LSQR迭代反演算法提供模型。层析结果显示出清晰的泥、砂互层，与测井参数构造解释对比，高、低速带与砂 ／泥 岩层基本吻合，最高能分辨 1／4地震波波长的薄层。研究证明级连算法能很好地兼顾高分辨率和高精度 ，大大提 高了井间地震走时层析成像的质量，可以广泛地应用于油田的勘探和检测。 关键词：过井地震层析成像 ；级联算法；高分辨率；阻尼LSQR迭代反演 中图分类号：P631．4 文献标识码：A 文章编号：1000—8918(2010)06—0814—07 过井地震走时层析成像因其简便、高效的优点， 用前一级的输出作为输入模型，在解空间对解估计 近年来在众多油田已经得到了广泛的应用，但 目前 的搜索沿着一条正确的路线进行，大大地提高了数 的过井地震走时层析成像使用的反演方法大多是低 值计算的稳定性。 分辨率反演方法，如联合迭代重建法 (SIRT)¨ ， 使用的一级反演方法为修正的BPT反演 ；二 共轭梯度最小二乘法 (CGLS)E3]等，这些方法的最 级反演方法为 LSQR阻尼迭代反演算法 ；射线 佳分辨率宽度约为震源和接收器密度的3倍 J。 追踪方法是 Langan等人 1985年提出的精确和快速 过井地震层析成像的成败关键是垂向分辨率提 射线追踪方法 。 高到 1／4地震波波长，才能满足勘探开发的要求 。 越是分辨率低的反演算法稳定性越好，要求输 只有高分辨率井间地震层析成像，才能作为油 田寻 入的初始信息越少，而且计算成本越低。由于层析 找利于油气聚集岩性圈闭、对油 田储集层描述和动 反演问题的非线性，成像分辨率高的方法要求接近 态检测的工具。跨孔地震层析成像级联方法 能 真实地球模型的准确背景波速场 ，否则反演算法可 够将分辨宽度提高到 1／4地震波波长，满足深入研 能不收敛，或者得到方差极大的解估计。因此使用 究地震资料的要求。 这样的思路：利用分辨率低而稳定性好的算法为分 跨孔地震层析成像级联方法 的基本思想是 辨率高的算法提供比较准确的初始模型，以保证在 在反演 中由低分辨率 向高分辨力逐步逼近 ，将具有 解空间中搜索解估计的过程沿一条正确的路线进 不同分辨率的算法 串连起来 ，通过非线性迭代逐步 行。所 以层析反演第一级使用分辨率低但稳定 的 提高成像的分辨率，同时抑制 由噪声带来的人造假 BPT算法为下一级反演提供初始模型。 象。将该方法应用在胜利油 田垦71区井间地震勘 井间地震走时层析成像是一个非线性问题，但 探中，并使用一套精确的井间地震数据处理方法，得 因为要同时确定未知的射线路径和未知速度结构， 到了该井间高分辨率的纵波速度结构 。将得到的层 一 般都使用线性反演方法求取速度模型。LSQR反 析成像结果与测井参数对 比分