一种根据实际钻井轨迹校正VSP下行波的处理方法

摘要

本发明属于地震勘探和开发领域，尤其涉及一种根据实际钻井轨迹校正VSP下行波的处理方法。所述方法包括步骤1：探测井眼轨迹数据和地震资料，步骤2：根据输入的地震数据，求取VSP检波器在深度h时与中垂线的距离d，走时钻井轨迹数据s和s1；步骤3：计算在深度为h时的走时路径误差步骤4：误差分析及校正步骤；步骤5：输出步骤。本发明的校正方法，解决了VSP下行波由路程误差进而引起由时间求得的速度误差，提高了VSP勘探的精度。该校正方法是在各向同性介质的条件下进行校正，利用钻井轨迹的方位与偏离数据校正VSP下行波资料，精确程度比较高。在各向异性非均匀介质的假设条件下使用则存在一定误差。

主权项

一种根据实际钻井轨迹校正VSP下行波的处理方法，其特征在于，所述方法包括如下步骤，步骤1：探测井眼轨迹数据和地震资料，探测得到钻井轨迹数据，并输入(1)钻井深度：h；(2)井斜角：钻井与垂直线的夹角；(3)方位角：在钻井深度h时钻井轨迹点与中心垂直线的方位角；探测得到VSP数据，并输入：(4)偏移距X：钻井井口到爆炸点的距离；(5)θ：以井口为原点，井口与爆炸点方位角；步骤2：根据输入的地震数据，求取VSP检波器在深度h时与中垂线的距离d，s为地面震源激发点到理论垂直钻井中VSP检波器波的走时路径和s1为地面震源激发点到实际钻井中VSP检波器波的走时路径；(1)求取d，d为钻井中VSP检波器在深度h时与中垂线的距离，d＝h·tgγ(2)求取S， s = x 2 + h 2 (3)求取S1， s 1 = x 2 + d 2 - 2 xd cos θ + h 2 步骤3：计算在深度为h时的VSP下行波的走时路径误差Δs： ΔS = s 1 - s = x 2 + d 2 - 2 xd cos θ + h 2 - x 2 + h 2 其中，x：为井口到地面震源激发点距离；h：为VSP检波器在井中接收的深度；s：为地面震源激发点到理论垂直钻井中VSP检波器波的走时路径；s1：为地面震源激发点到实际钻井中VSP检波器波的走时路径；d：钻井中VSP检波器波偏离钻井井口中垂线的距离；θ：以井口为圆点，VSP检波器波与震源投影在水平面上的夹角；步骤4：误差分析及校正步骤根据步骤3得到的误差，进行分析校正：对于走时路径误差进行相对误差分析和绝对误差分析，并且进行校正；把实际钻井中VSP检波器波的位置平移d到井口中心垂直线上，其中d为钻井中VSP检波器在深度h时与中垂线的距离，深度h不变；并且求出地面震源激发点与在钻井井口中心垂直线上VSP检波器的走时。由于VSP检波器在钻井中接收到的是地面激发波的振幅，波的振幅最先到达时间t0称为初至(t0为原始观测数据)；而这一段时间与走时路径的比V＝s1/t0为平均速度；步骤包括，(1)求由地面震源激发点与在钻井中VSP检波器的走时路径与初至时间t0的比V＝s1/t0，该比为平均速度；(2)把实际钻井中VSP检波器波的位置平移d(d为钻井中VSP检波器在深度h时与中垂线的距离)到井口中心垂直线上(深度h不变)；(3)由求得的钻井中心垂直线上VSP检波器波的位置与地面震源激发点距离s(s为到理论垂直钻井中VSP检波器波的走时路径)，用V＝s1/t0为平均速度求得钻井中心垂直线上VSP检波器波的位置与地面震源激发点走时t0o，t0o＝s/v；(4)求由地面震源激发点与在钻井井口中心垂直线上VSP检波器的走时t0o与初至时间t0(t0为原始观测数据)的差Δt，Δt＝t0‑t0o；(5)对原始观测数据t0进行校正，即t0o＝t0‑Δt，也就是求得在所有深度的中地面震源激发点与在钻井井口中心垂直线上VSP检波器的精确走时t0o；(6)误差分析过程：1)计算误差走时路径绝对误差： ΔS = s 1 - s = x 2 + d 2 - 2 xd cos θ + h 2 - x 2 + h 2 , 计算走时路径相对误差：Δs/s1计算走时绝对误差：Δt＝t0‑t0o计算走时相对误差：Δt/t02)分析误差：a，分析走时路径绝对误差，找出最大绝对误差所在深度，最大数值，对实际生产的影响的程度。也可以分析一段所在深度绝对误差比较大的对实际生产的影响的程度；b，分析走时路径相对误差，找出最大走时路径相对误差所在深度，最大数值，对实际生产的影响的程度；也可以分析一段所在深度走时路径相对误差比较大的对实际生产的影响的程度；步骤5：输出结果输出深度(m)，井斜角(度)，方位角(度)，地面震源激发点到中心垂直线上VSP检波器波走时路径(m)，地面震源激发点到实际钻井中VSP检波器波的走时路径(用钻井轨迹)(m)，走时路径误差(m)，走时路径相对误差(％)，地面震源激发点到中心垂直线上VSP检波器波走时(ms)，地面震源激发点到实际钻井中VSP检波器波的走时(用钻井轨迹)(ms)，走时误差(ms)，相对时间误差(％)。