CN103645506B - 一种检测地层中裂缝发育程度的方法 - Google Patents

Abstract

本发明提供一种检测地层中裂缝发育程度的方法，所述方法利用横波在裂缝发育地层中传播时发生横波分裂现象，导致井下三分量检波器的两个水平分量接收到的横波初至波的波形随着地层中裂缝的发育程度变化而变化，通过计算两个水平分量上的横波初至波互相关误差系数来定量检测地层裂缝的发育程度，为降低勘探开发的风险，提高油气藏开采效率提供参考数据。

Description

一种检测地层中裂缝发育程度的方法

技术领域

本发明涉及地震勘探技术领域，是一种通过对三维三分量VSP中井下三分量检波器的两个水平分量记录的地震波场的横波初至能量的比值进行计算，从而对地层中裂缝发育程度进行定量检测的分析方法。本发明主要应用于地震勘探领域，为裂缝型油气藏勘探与开发提供地层裂缝发育程度信息，为降低勘探开发的风险，提高油气藏开采效率提供参考数据。

背景技术

在三维三分量VSP地震勘探数据中，当横波可控震源在井旁激发时，震源所产生的横波在地层中传播，穿透地层到达井中的三分量检波器所在位置，并在该点引起质点振动，该点的振动被放置于井中的三分量检波器记录下来。如果地层中裂缝发育，可控震源所激发的横波在穿透地层传播时会发生横波分裂现象，即横波会分裂为沿着裂缝走向偏振的快横波和垂直裂缝走向偏振的慢横波。通过计算可控震源激发、井中三分量检波器接收的三维三分量VSP地震数据中井下三分量检波器的水平分量的横波初至波互相关误差系数的大小，可以定量检测地层中裂缝发育的程度。

现有技术手段是通过地层的地质露头进行观测、或者对钻井中提取的岩心样本在实验室中观测、或者对测井资料进行分析来获得地层裂缝发育程度的相关信息。由于上述方法资料采集的局限性，只能够获得地层裂缝的局部信息。

现有技术手段也可以通过宽方位角纵波地震勘探获得纵波地震传播的方位各向异性属性来检测地层裂缝的走向方位角度信息，但是由于地震纵波对地层裂缝的发育情况不够敏感，因此通过纵波地震方位各向异性属性分析来对储层裂缝的发育特征只能够做到定性的或者是误差较大的半定量描述。

总之，现有的技术手段由于受到技术方法本身的限制，对地层裂缝发育程度的检测能力不能满足油田的生产需要。

发明内容

本发明通过计算三维三分量VSP地震勘探中井下三分量检波器水平分量的横波初至的互相关误差系数来定量检测地层裂缝发育程度的分析方法与流程。

当地层中发育有裂缝时，横波在地层中传播时将发生横波分裂现象，横波分裂导致井下三分量检波器的两个水平分量上接收到的横波初至波的能量不同，地层中裂缝发育程度越高，井下三分量检波器的两个水平分量上接收到的横波初至波的能量大小就越接近。通过计算井下三分量检波器的两个水平分量上所记录的地震波场中横波初至波互相关误差系数,可以定量检测地层中裂缝的发育程度。

图1是当横波穿过呈平行排列的垂直裂缝介质时发生横波分裂的原理示意图。当由横波可控震源激发的入射横波S穿过裂缝发育的地层时，发生了横波分裂，入射横波S分解为快横波S₁和慢横波S₂，其中快横波S₁的偏振方向平行于裂缝走向，慢横波S₂的偏振方向垂直于裂缝走向。

为了便于分析横波分裂现象，我们将三分量检波器的径向方向记录的地震波场记为R(t)、切向分量记录的地震波场记为T(t)，三分量检波器的两个水平分量同时记录下了快横波S₁(t)和慢横波S₂(t)在观测点引起的质点振动，三分量检波器的两个水平分量R(t)与T(t)可以写成如下形式：

$\left[\begin{array}{c}R\left(t\right)\\ T\left(t\right)\end{array}\right]=\left[\begin{array}{cc}\cos \mathrm{\θ}& \sin \mathrm{\θ}\\ \sin \mathrm{\θ}& -\cos \mathrm{\θ}\end{array}\right]\left[\begin{array}{c}{S}_1\left(t\right)\\ S_2\left(t\right)\end{array}\right]+\left[\begin{array}{c}{N}_R\left(t\right)\\ N_T\left(t\right)\end{array}\right]---\left(1\right)$

上式中，N_R(t)与N_T(t)分别表示三分量检波器的径向分量与横向分量记录的环境噪声；θ是地层中裂缝走向方位角度β与三分量检波器的径向分量指向的方位角度α的夹角，见图2中所示。

本发明的技术方案为：

一种检测地层中裂缝发育程度的方法，所述方法包括以下步骤：

(1)根据所欲观测的井所在的现场地形，按照如下要求设计三维三分量VSP观测方案：从井口所在位置O点指向可控震源所在位置S点的观测方位角度取值α∈{α_{j,j＝1,2,......,M}}，α_j∈[0,360)，α的M个取值将圆周均匀分割为M等份，根据对观测精度的要求，M取值从集合{6、8、10、12、16、18、36}中选取一个合适的数值；井下三分量检波器所放置的位置G与井口位置O的距离为z₀米，z₀的取值保证井下三分量检波器的正常工作并尽量接近要观测的地层；可控震源所在位置S与井口位置O之间的距离为x₀米，x₀的取值满足当激发可控震源时，井中三分量检波器能够接收到较高信噪比的有效信号；设定井下三分量检波器的离散采样间隔dt与离散采样点数N，dt与N的取值保证记录到完整的可控震源引发的位置G处的质点振动；

(2)对于α∈{α_{j,j＝1,2,......,M}}的每一个取值α_j，将可控震源放置在与井口O距离为x₀米方位角度为α_j处；将井下三分量检波器放置于与井口位置O的距离为z₀米的井下，并使井下三分量检波器的R分量指向的方位角度与α_j保持一致；激发可控震源产生沿着α_j方位角度偏振的横波，并同时以井下三分量检波器记录可控震源引发的井下三分量检波器所在位置G处的质点振动，当α取值为α＝α_j时井下三分量检波器两个水平分量R、T记录到的质点振动记为R_j(t_k),T_j(t_k),k＝1,2,......,N；

(3)对步骤(2)中得到的所有观测方位的数据记录{R_j(t_k),T_j(t_k)},j＝1,2,......,M,k＝1,2,......,N，编程计算直达波初至的到达时间{P_j}，j＝1,2,......,M，

编程的具体算法为：

对于α＝α_j的每一个取值得到的井下三分量检波器的两个水平分量的数据记录R_j(t_k),T_j(t_k),k＝1,2,......,N，扫描计算{r_k＝F_k/E_k,k＝L+1,L+2,......,N-L}中的最大值所对应的k的取值，则得到P_j＝kdt；在公式{r_k＝F_k/E_k,k＝L+1,L+2,......,N-L}中， $E_k=\underset{l=1}{\overset{L}{\mathrm{\Σ}}}(R_j^2\left(t_{k-l}\right)+T_j^2\left(t_{k-l}\right)),$ L的取值由可控震源激发的地震子波的延续时间t0决定，L＝t0/dt，L＜N；

(4)根据{P_j}的值，求取在所有观测方位角度观测的R分量上的横波初至与T分量上横波初至的互相关误差系数 $C_j=1-\frac{\left|\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{R}_j\left(t_{k+1}\right)T_j\left(t_{k+1}\right)\right|}{\sqrt{\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{R}_j^2\left(t_{k+1}\right)}\sqrt{\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{T}_j^2\left(t_{k+1}\right)}},k=P_j/,$ L的取值同上所述，j=1,2,……,M；

(5)计算所有C_j，j=1,2,……,M的平均值则sigma值的大小用于定量描述所观测地层中的裂缝发育程度，sigma的取值范围为sigma∈（0,1），当sigma值增大时，代表地层中的裂缝发育程度提高，当sigma值减小时代表地层中的裂缝发育程度降低；

(6)流程结束。

本发明的技术效果：利用横波在裂缝发育地层中传播时发生横波分裂现象，导致井下三分量检波器的两个水平分量接收到的横波初至波能量比值随着地层中裂缝的发育程度变化而变化，通过计算两个水平分量上的横波初至波互相关误差系数来定量检测地层裂缝的发育程度。

附图说明

图1为横波穿过呈平行排列的垂直裂缝介质时发生横波分裂的示意图。

图2为三分量检波器径向分量指向的方位角度α、裂缝走向β与θ的关系图。

图3为三维三分量VSP地震观测方法示意图。

图4为α取值（M=16）及其与井下检波器水平分量指向方位关系示意图。

图5与图6分别是理论合成的三维三分量VSP观测数据的水平分量R与T图。模型设计的裂缝方位为NE60°，地层速度100m/s，x=500m，z=500m。观测方位角度α的取值为：NE10°,NE20°,……,NE360°。井下三分量检波器的方向采用右手坐标系,径向分量与观测方位角度α保持一致，切向分量垂直于观测方位角度α，位于径向分量顺时针旋转90°的方位。合成数据地震子波采用的是40Hz的零相位雷克子波。

图7和图8是对图5和图6中的合成地震记录水平分量进行横波初至时间拾取的结果图。

图9是对图7和图8中的数据计算每个观测方位RT分量互相关误差系数数值图，从图9中可知所有观测方位RT分量互相关误差系数的平均值为sigma=0.62，表明地层中裂缝的发育程度为0.62。

具体实施方式

一种检测地层中裂缝发育程度的方法，所述方法包括以下步骤：

(1)根据所欲观测的井所在的现场地形，按照如下要求设计三维三分量VSP观测方案：参考图2所示，从井口所在位置O点指向可控震源（炮点）所在位置S点的观测方位角度取值α∈{α_{j,j＝1,2,......,M}}，α_j∈[0,360)，α的M个取值将圆周均匀分割为M等份，根据对观测精度的要求，M取值可从集合{6、8、10、12、16、18、36}中选取一个合适的数值；井下三分量检波器所放置的位置G与井口位置O的距离为z₀米，z₀的取值应能够保证井下三分量检波器的正常工作并尽量接近要观测的地层；可控震源所在位置S与井口位置O之间的距离为x₀米，x₀的取值应满足当激发可控震源时，井中三分量检波器能够接收到较高信噪比的有效信号；设定井下三分量检波器的离散采样间隔dt与离散采样点数N，dt与N的取值应能保证记录到完整的可控震源引发的G点处的质点振动。

(2)对于α∈{α_{j,j＝1,2,......,M}}的每一个取值α_j，将可控震源放置在与井口O距离为x₀米方位角度为α_j处；将井下三分量检波器放置于与井口位置O的距离为z₀米的井下，并使井下三分量检波器的R分量指向的方位角度与α_j保持一致（参考图3所示）；激发可控震源产生沿着α_j方位角度偏振的横波，并同时以井下三分量检波器记录可控震源引发的井下三分量检波器所在位置G处的质点振动。当α取值为α＝α_j时井下三分量检波器两个水平分量R、T记录到的质点振动记为R_j(t_k),T_j(t_k),k＝1,2,......,N。

(3)对(2)中得到的所有观测方位的数据记录{R_j(t_k),T_j(t_k)},j＝1,2,......,M,k＝1,2,......,N，编程计算直达波初至的到达时间{P_j}，j＝1,2,......,M。

编程的具体算法为：

对于α＝α_j的每一个取值得到的井下三分量检波器的两个水平分量的数据记录R_j(t_k),T_j(t_k),k＝1,2,......,N，扫描计算{r_k＝F_k/E_k,k＝L+1,L+2,......,N-L}中的最大值所对应的k的取值，则得到P_j＝kdt；在公式{r_k＝F_k/E_k,k＝L+1,L+2,......,N-L}中， $E_k=\underset{l=1}{\overset{L}{\mathrm{\Σ}}}(R_j^2\left(t_{k-l}\right)+T_j^2\left(t_{k-l}\right)),$ L的取值由可控震源激发的地震子波的延续时间t0决定，L＝t0/dt，L＜N；

(4)根据{P_j}的值，求取在所有观测方位角度观测的R分量上的横波初至与T分量上横波初至的互相关误差系数 $C_j=1-\frac{|\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{R}_j\left(t_{k+l}\right)T_j\left(t_{k+l}\right)}{\sqrt{\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{R}_j^2\left(t_{k+l}\right)}\sqrt{\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{T}_j^2\left(t_{k+l}\right)}},k=P_j/\mathrm{dt},$ L的取值同上所述，j=1,2,……,M；

(5)计算所有C_j，j=1,2,……,M的平均值则sigma值的大小可用于定量描述所观测地层中的裂缝发育程度，sigma的取值范围为sigma∈（0,1），当sigma值增大时，代表地层中的裂缝发育程度提高，当sigma值减小时降低；

(6)流程结束。

本发明利用横波在裂缝发育地层中传播时发生横波分裂现象，导致井下三分量检波器的两个水平分量接收到的横波初至波能量比值随着地层中裂缝的发育程度变化而变化，通过计算两个水平分量上的横波初至波互相关误差系数来定量检测地层裂缝的发育程度。

Claims (1)

1.一种检测地层中裂缝发育程度的方法，所述方法包括以下步骤：

(1)根据所欲观测的井所在的现场地形，按照如下要求设计三维三分量VSP观测方案：从井口所在位置O点指向可控震源所在位置S点的观测方位角度取值α∈{α_{j,j＝1,2,......,M}}，α_j∈[0,360)，α的M个取值将圆周均匀分割为M等份，根据对观测精度的要求，M取值从集合{6、8、10、12、16、18、36}中选取一个合适的数值；井下三分量检波器所放置的位置G与井口位置O的距离为z₀米，z₀的取值保证井下三分量检波器的正常工作并尽量接近要观测的地层；可控震源所在位置S与井口位置O之间的距离为x₀米，x₀的取值满足当激发可控震源时，井中三分量检波器能够接收到较高信噪比的有效信号；设定井下三分量检波器的离散采样间隔dt与离散采样点数N，dt与N的取值保证记录到完整的可控震源引发的位置G处的质点振动；

(2)对于α∈{α_{j,j＝1,2,......,M}}的每一个取值α_j，将可控震源放置在与井口O距离为x₀米方位角度为α_j处；将井下三分量检波器放置于与井口位置O的距离为z₀米的井下，并使井下三分量检波器的R分量指向的方位角度与α_j保持一致；激发可控震源产生沿着α_j方位角度偏振的横波，并同时以井下三分量检波器记录可控震源引发的井下三分量检波器所在位置G处的质点振动，当α取值为α＝α_j时井下三分量检波器两个水平分量R、T记录到的质点振动记为R_j(t_k),T_j(t_k),k＝1,2,......,N；

(3)对步骤(2)中得到的所有观测方位的数据记录{R_j(t_k),T_j(t_k)},j＝1,2,......,M,k＝1,2,......,N，编程计算直达波初至的到达时间{P_j}，j＝1,2,......,M，

编程的具体算法为：

对于α＝α_j的每一个取值得到的井下三分量检波器的两个水平分量的数据记录R_j(t_k),T_j(t_k),k＝1,2,......,N，扫描计算{r_k＝F_k/E_k,k＝L+1,L+2,......,N-L}中的最大值所对应的k的取值，则得到P_j＝kdt；在公式{r_k＝F_k/E_k,k＝L+1,L+2,......,N-L}中， $E_k=\underset{l=1}{\overset{L}{\mathrm{\Σ}}}(R_j^2\left(t_{k-l}\right)+T_j^2\left(t_{k-l}\right)),$ $F_k=\underset{l=1}{\overset{L}{\mathrm{\Σ}}}(R_j^2\left(t_{k+l}\right)+T_j^2\left(t_{k+l}\right)),$ L的取值由可控震源激发的地震子波的延续时间t0决定，L＝t0/dt，L＜N；

(4)根据{P_j}的值，求取在所有观测方位角度观测的R分量上的横波初至与T分量上横波初至的互相关误差系数 $C_j=1-\frac{\left|\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{R}_j\left(t_{k+l}\right)T_j\left(t_{k+l}\right)\right|}{\sqrt{\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{R}_j^2\left(t_{k+1}\right)}\sqrt{\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{T}_j^2\left(t_{k+1}\right)}},k=P_j/,$ L的取值同上所述，j=1,2,……,M；

(5)计算所有C_j，j=1,2,……,M的平均值则sigma值的大小用于定量描述所观测地层中的裂缝发育程度，sigma的取值范围为sigma∈（0,1），当sigma值增大时，代表地层中的裂缝发育程度提高，当sigma值减小时代表地层中的裂缝发育程度降低；

(6)流程结束。