速度模型精度对vsp资料偏移成像质量的影响.docxVIP

速度模型精度对vsp资料偏移成像质量的影响 vsp检测具有信噪比高、分辨率高的优点，受到了业界的高度重视。VSP资料波场信息丰富 (如:直达下行波、上行反射波、大角度出射波以及多次波等) , 有利于陡倾角成像、多次波成像和多波成像。常规VSP资料偏移成像包括VSP-CDP变换方法、Kirchhoff积分法和波动方程偏移法等。然而, 在介质复杂的情况下, 这些方法都存在成像精度不高、不能对多次波准确成像以及保幅性较差等问题。逆时偏移利用双程波方程进行波场延拓, 其偏移结果不受地层倾角的限制, 并且能够对多次波进行成像 然而, 在复杂构造情况下 (特别是复杂上覆构造) , 受VSP资料速度分析精度的限制, 很难得到较为准确的速度模型, 因此, VSP资料逆时偏移质量也受到极大的影响。鉴于复杂上覆构造的影响, 可以转向另一种处理思路, 即避开上覆复杂构造的影响, 专门针对井周区域进行局部成像。局部成像方法受偏移速度精度的影响较小, 对提高VSP资料的成像精度有较大贡献。 很多学者对VSP资料局部成像方法进行了研究。一类是基于干涉成像 我们对不受上覆复杂构造影响的VSP资料局部逆时偏移方法进行研究, 重点讨论直达波旅行时重构和基于激发时间成像条件局部逆时偏移两项关键技术。借鉴了He等 1 方法原理 1.1 逆时偏移研究 传统意义上, 偏移成像包括两步———利用波场外推算子重构波场, 然后利用成像条件将反传到散射点处的能量作为“像”提取出来。Claerbout 基于激发时间成像条件的逆时偏移可表达为 式中:t 基于激发时间成像条件的逆时偏移需要计算震源到各散射点的旅行时场, 一般的处理思路是基于波的正向传播, 从震源向散射点计算旅行时。然而, 当上覆地层比较复杂时, 速度模型的精度不会很高, 计算出的旅行时会有较大误差, 在很大程度上影响偏移成像的质量。对于VSP资料来说, 检波器更加靠近目的层, 这就给另外一种思路提供了便利 VSP资料局部逆时偏移的技术路线为 (图2) :首先, 数据准备———从VSP资料中分离出上行P波、拾取直达波初至以及建立观测井段深度范围的速度模型;接着, 利用直达波初至时间重构出直达波的旅行时场;然后, 对分离出的上行P波数据进行反向时间延拓, 利用激发时间成像条件进行偏移成像;最后, 通过去噪和切除得到最终成像结果。 1.2 源点旅行时重构 逆向旅行时计算是从检波点向震源点方向计算, 离震源点越近旅行时将越小。逆向旅行时计算依然需要满足程函方程, 未知点的旅行时为所有已知点的旅行时减去两者间局部走时后的最大值 (图3) 。逆向旅行时计算与正向旅行时计算相类似, 只是计算方向不同, 因此, 可以借鉴正向旅行时计算的相关原理。关于正向旅行时计算的研究已经非常深入, 可以利用射线追踪 如图4所示, 已知 (x 任意一点 (x, z) 到源点的旅行时关系可表达为 其中, 源点位于 (0, 0) ;S 在波传播距离较远的情况下, (x 利用线性关系, (x 假设, 过点 (x 其中, S为当前网格点的慢度。利用 (7) 式对z dt/dz 逆向旅行时计算过程中, 必须按照波传播过程的因果律, 制定合适的计算策略。违反因果律的任何计算方法或计算过程都不可能得到正确的旅行时。假设已知第一列的旅行时, 从左向右计算旅行时, 在实际处理过程中, 我们并不清楚射线的方向, 计算每一点的旅行时时必须将所有可能的射线方向都考虑到 (图5) 。对地震射线可能的方向分别计算旅行时, 然后取最大旅行时动态更新当前计算点的旅行时 (可参见文献[15]中的图3) 。 为了验证上述直达波旅行时重构方法的效果, 对常梯度模型进行了测试。速度模型如图6a所示, 源点位于右上角;图6b为正演的旅行时;图6c为利用最左侧记录的旅行时重构出的直达波旅行时;图6d为两者间旅行时间的差异。从图6中可以看出, 本文方法可以很好地重构直达波旅行时。但是, 旅行时重构也存在不准确区域, 该区域是连接源点和检波点的直达波地震射线不经过的区域。由于直达波初至时间不包含该误差区域的信息, 此处旅行时不能被准确重构。那么, VSP局部逆时偏移方法成像区域将有所限制, 只能对井周围有限区域成像。所以, 对于仅针对目标区域成像来说, 本文方法是有效的。 2 主要逆时偏移结果 首先利用模拟资料验证本文方法的有效性。设计一个6层速度模型 (图7) , 利用声波方程有限差分法进行VSP资料模拟。共模拟10炮数据, 炮点偏移距从400m到960m, 偏移间隔为40m;井口位于0处, 井中观测井段300~2 000m, 观测点距为5m。图8为炮点偏移距为400m时的VSP波场模拟记录。 对图8所示模拟资料进行直达波初至拾取和上行波分离;然后分别进行常规逆时偏移 (基于相关成像条件)