[发明专利]基于方位弹性阻抗傅里叶级数展开的裂缝预测方法及系统

权利要求书

1.一种基于方位弹性阻抗傅里叶级数展开的裂缝预测方法，包括：

根据原始道集，获得包含方位角信息的角道集；

根据方位角与入射角将所述角道集进行角度划分，获得不同方位角、不同入射角的部分叠加地震数据；

根据测井的纵波速度、横波速度、密度、各向异性信息，计算与所述部分叠加地震数据对应的每一个方位角、每一个入射角的方位弹性阻抗曲线；

根据每一个方位角、每一个入射角的方位弹性阻抗曲线开展正演，与对应的部分叠加地震数据进行匹配标定，提取每一个方位角、每一个入射角对应的地震子波，进而建立每一个方位角、每一个入射角对应的弹性阻抗模型；

根据所述部分叠加地震资料、所述弹性阻抗模型、所述地震子波，计算每一个方位角、每一个入射角对应的方位弹性阻抗体，并对方位弹性阻抗体进行傅里叶级数展开，实现裂缝预测。

2.根据权利要求1所述的基于方位弹性阻抗傅里叶级数展开的裂缝预测方法，其中，通过公式(1)计算方位弹性阻抗曲线：

其中，EI(θ，φ)为方位弹性阻抗曲线，g为横波速度与纵波速度的比值的平方，V_p为纵波速度，V_s为横波速度，ρ为密度，ε^(v)、δ^(v)、γ分别表示为Thomsen弱各向异性参数，θ表示地震波入射角，φ表示为测线方向与裂缝倾向的夹角，V_p0、V_s0、ρ₀分别表示为纵波速度、横波速度、密度的均值。

3.根据权利要求1所述的基于方位弹性阻抗傅里叶级数展开的裂缝预测方法，其中，通过公式(2)计算方位弹性阻抗体：

其中，F(m)为目标函数，G为正演算子，m为模型参数向量，d为观测数据向量，i＝1,2，…，N，N为采样点数量，X_n为地震数据的协方差矩阵，C为积分矩阵，ξ为相对弹性阻抗，表示第i组模型参数向量，表示第i组模型参数向量的协方差矩阵，X_Cm-ξ表示模型误差的协方差矩阵。

4.根据权利要求1所述的基于方位弹性阻抗傅里叶级数展开的裂缝预测方法，其中，通过公式(3)进行傅里叶级数展开：

ln(EI_A(θ，φ))＝A₀+A₂cos2φ+A₄cos4φ (3)

其中，A₀、A₂与A₄分别为零阶、二阶、四阶傅里叶系数，表达式分别为：

其中，EI_A(θ,φ)表示归一化后的方位弹性阻抗，g为横波速度与纵波速度的比值的平方，V_p为纵波速度，V_s为横波速度，ρ为密度，ε^(v)、δ^(v)、分别表示为Thomsen弱各向异性参数，φ表示为测线方向与裂缝倾向的夹角，θ表示地震波入射角。

5.根据权利要求1所述的基于方位弹性阻抗傅里叶级数展开的裂缝预测方法，其中，所述部分叠加地震数据的方位角与入射角分别为划分区域的方位角与入射角的平均值。