[发明专利]有限方位角条件下压噪叠前纵波各向异性裂缝检测方法

权利要求书

有限方位角条件下压噪叠前纵波各向异性裂缝检测方法

1.有效的有限方位角条件下压噪叠前纵波各向异性裂缝检测方法，其特征在于，包括如下步骤：

步骤1：对地震资料进行分方位处理(包括保幅偏移处理)，抽取分方位共反射点(CRP)道集；

步骤2：将分方位的CRP道集进行数据合并，并为下一步作准备；

步骤3：用本发明所采用的算法进行裂缝反演，算出每个CRP点上每个时间采样点上的裂缝密度和裂缝方位角值；

步骤4：对反演数据进行质量控制；

步骤5：输出裂缝反演结果即裂缝密度数据体和方位角数据体；

步骤6：对输出后的数据体提取层位数据，并进行地质意义分析。

2.根据权利要求1所述的有限方位角条件下压噪叠前纵波各向异性裂缝检测方法，其特征在于，所述的步骤2中的算法，从世界著名的Ruger反射系数公式出发，当覆盖次数为n时，一般的裂缝反演问题可归结为求解线性方程组AX＝B；其中A是n行4列的矩阵；X是四个未知量构成的列向量；B是不含噪音的地震反射振幅；其不足之处是，只能使用宽方位采集的地震资料进行反演，以克服噪音对反演结果的不利影响；当使用有限(窄)方位角地震资料时，含有噪音的资料会使反演的结果失效或不准确；

本发明的改进之处在于通过添加噪音因子N，使要求解的方程转化为AX＝B+N

根据奇异值分解的方法求解A的广义逆，进而解上述方程可得：

X＝(A^TA)^-1A^T(B+N)

 ＝V^TW^-1U^TB+V^TW^-1U^TN

式中U是n行4列正交矩阵；W是4行4列对角矩阵，对角元素是正值或0；V是4行4列正交矩阵；通过以上方程，将地震数据的有效信号和噪音项分开；具体方式是将N乘上一个很小的数或除以一个很大的数，从而使噪音项趋近于0；

本算法的特色是可以利用含有噪音的有限方位角资料，准确地反演出裂缝的各向异性特征，极大地节约采集所需的成本。

3.根据权利要求1所述的有限方位角条件下压噪叠前纵波各向异性裂缝检测方法，其特征在于，裂缝探测的适用范围如下：所述的步骤2中反演的裂缝密度，是指有效裂缝引起的各向异性强度；所谓有效裂缝是指开启并充填了流体的裂缝；有效裂缝是油、气、水的有效储集空间；当一定深度范围内有效裂缝走向大致相同时，各向异性程度高，反之，若在一定的深度范围内，裂缝的走向呈现出不同的方向，地震波的各向异性特征会被削弱，甚至会显示出各向同性的特征，无法通过反演获得地下裂缝的信息。