VSP多波旅行时反演探究.pdfVIP

　２０１１ 年 １１ 月 第 ４６ 卷 　 增刊 １ 　 ＊ 北京市昌平区府学路 １８ 号中国石油大学（北京）地质地球物理综合研究中心， １０２２４９ 。 Ｅｍａｉｌ ： ｓａｍｚｄｓｕｎ ＠ｙ ａｈｏｏ．ｃｏｍ ， ｓｚｄ ＠ ｃｕ ｐ ．ｅｄｕ．ｃｎ 本文于 ２０１１ 年 ３ 月 １８ 日收到，最终修改稿于 ２０１１ 年 ９ 月 ３０ 日收到。 本文研究由国家重点基础研究发展规划 ９７３ 项目（ ２０１１ＣＢ２０１１０３ ）“中国西部叠合盆地深部有效碳酸盐岩储层形成机制与分布预测”与国家 油气重大专项（ ２０１１Ｚ）“海相碳酸盐岩储层地震描述与油气藏有效预测技术研究”资助。 · ＶＳＰ 技术· 文章编号： １０００ － ７２１０ （ ２０１１ ）增刊 １ － ００４０ － ０７ ＶＳＰ 多波旅行时反演 谢春辉①② 　 孙赞东＊① 　 王学军③ 　 杨午阳② （ ① 中国石油大学（北京）地质地球物理综合研究中心，北京昌平 １０２２４９ ； ② 中国石油勘探 开发研究院西北分院，甘肃兰州 ７３００２０ ； ③ 东方地球物理公司研究院，河北涿州 ０７２７５１ ） 谢春辉，孙赞东，王学军，杨午阳 ．ＶＳＰ 多波旅行时反演 ． 石油地球物理勘探， ２０１１ ， ４６ （增刊 １ ）： ４０ ～ ４６ 摘要 　 本文基于几何地震学原理，分别利用 ＶＳＰ 直达波、反射纵波及上行转换横波旅行时反演层速度。直达波 旅行时反演只能获得检波点所在井段地层的速度；上行纵波旅行时反演能获得最深检波器以下地层的速度，即 井底以下地层的速度；上行转换横波旅行时反演则能获得横波速度。本文主要研究利用 ＶＳＰ 直 达 波 和 上 行 反 射波求取井段及井底以下地层的速度，最终根据三维 ＶＳＰ 资料获取三维速度体。 关键词 　 速度分析 　 旅行时反演 　 多波 中图分类号： Ｐ６３１ 　　 文献标识码： Ａ １ 　 引言 通过对零井源距 ＶＳＰ 资料进行初至拾取，可以 直接获得时间—深度关系，由此计算地震平均速度、 层速度等参数；也可以进行大套速度界面的划分及 地质层位的划分。除零井源距 ＶＳＰ 外，人们还希望 从非零井源距及三维 ＶＳＰ 资料中获取速度信息。 有关应用 ＶＳＰ 数据进 行 速 度 分 析 的 研 究 由 来 已久， Ｄｉｘ ［ １ ］在 １９４５ 年讨论了单一震源的直井 ＶＳＰ 计 算 层 速 度 的 方 法； Ｓｔｅｗａｒｔ ［ ２ ～ ４ ］ 、 Ｌｉｎｅｓ ［ ５ ］ 、 Ｓａｌｏ 等［６］及 Ｌｅｅ ［ ７ ］对旅行时反演进行过系统的研究。贺 元凯［８］研究了不同井源距情况的速度求取，并且分 析了初至拾取误差对结果的影响；王成礼等［９］根据 折射定律，给出了关于 ＶＳＰ 直达 波 的 方 程 组，已 知 检波点深度和初至时间求解速度；周熙襄等［１０］研究 了用旅行时反演确定层速度的横向变化；陈信平［１１］ 研究了 ＶＳＰ 初至时间反演层速度的算法 和 误 差 估 计；庄东海［１２］提 出 了 射 线 追 踪 折 线 法，逐 层 迭 代 求 解地层速 度；邬 世 英 等［１３］改 进 雅 克 比 矩 阵 求 取 办 法，使得反演过程更加稳定、收敛速度更快。 从严格意义上说，基于直达波初至时间计算速 度的方法，只是获得每两个检波点之间的平均速度。 当地下接收 点 距 较 大，相 应 地 层 厚 度 较 薄 时，求 得 的速度并不能真正反映地层的速度；当接收点距较 小时，所 计 算 出 来 的 速 度 受 初 至 时 间 的 影 响 又 较 大，初至时间的较小变化可能引起层速度计算上的 较大误差。利用直达波信息，只能解决最深检波点 以上地层的相对准确的速度，最深检波点以下的地 层速度一般采用插值外推的算法求得，存在先天的 缺陷。 随着物探技术的发展和应用领域的扩大，人们 求解速度的视野也逐 渐 变 宽。朱 键 等［１４］利 用 ＶＳＰ 上行波进行旅行时反演；徐仲达等［１５］、张继国［１６］研 究了 ＶＳＰ 横波反演的可行性。为 弥 补 井 底 以 下 地 层速度分析的盲区，本文将旅行时反演由直达波推 广至上行反射波和上行转换横波，其中确定反射同 相轴的位置、拾取旅行时是反演中至关重要的步骤。 ２ 　 旅行时反演原理 ＶＳＰ 各波的传播路径如图 １ 所示，绿色代表直 达纵波，蓝色 表 示 上 行 纵 波，红 色 表 示 上 行 转 换 横 波。已知井源距 和 第 ｉ 个 检 波 点 的 观 测 旅 行 时 ｔ 　 ｉ ０ ， 以上三种波的旅行时信息均能反演相应深度段的地 层速度。 　第 ４６ 卷 　 增刊 １ 谢春辉等： ＶＳＰ 多波旅行时反演 ４１ 　 　 图 １ 　 水平层状介质直达波、反射纵波及 转换横波传播路径