[发明专利]一种基于叠前地震数据的煤层气储层脆性指数反演方法

权利要求书

1.一种基于叠前地震数据的煤层气储层脆性指数反演方法，其特征在于，包括以下步骤：

一、叠前三维地震数据处理

针对叠前地震数据的CDP道集，抽取16个面元叠加形成超道集，进一步结合声波时差测井数据建立的速度场信息生成角度道集d(θ)；

利用声波时差和密度测井曲线制作合成记录，通过与井旁地震记录做统计分析，提取得到地震子波矩阵W；

二、测井数据处理

利用P波速度、S波速度和密度测井曲线，依据公式计算得到动弹性参数(E_d和υ_d)测井曲线，

式中，V_P和V_S分别表示纵波速度和横波速度，ρ为密度，E_d和υ_d分别表示动杨氏模量和动泊松比；

三、反演初始模型建立

1)、对步骤二中计算得到的动弹性参数、密度测井曲线进行低通滤波，得到反演参数的低频趋势，即反演的初始模型；

2)、计算反演方程中的系数因子:为正则化系数，为m的偏导数算子；

3)、结合提取的地震子波矩阵W和系数矩阵A(θ)得到G(θ)＝WA(θ)，

式中，σ_n为噪声的标准方差，σ_m为反演参数m的标准方差，为横波速度与纵波速度的比，θ表示入射角度；

四、动弹性参数反演

将反演初始模型的参数m、系数因子与处理后的角道集数据d(θ)带入到反演方程中，利用共轭梯度法求解，通过更新反演参数m和系数因子迭代5次后得到动弹性参数(E_d和υ_d)的反演结果；

五、全区目标煤层气储层脆性指数计算

根据全区的动弹性参数测井曲线统计结果，确定动弹性参数的最大和最小值(E_dmax、E_dmin、υ_dmax和υ_dmin)，基于脆性指数计算公式，利用动弹性参数反演结果计算得到全区的脆性指数分布。

2.根据权利要求1所述的一种基于叠前地震数据的煤层气储层脆性指数反演方法，其特征在于，所述步骤二中依据的公式为：

3.根据权利要求1所述的一种基于叠前地震数据的煤层气储层脆性指数反演方法，其特征在于，所述步骤三中的反演方程为：(G(θ)^TG(θ)+δQ(m))m＝G(θ)^Td(θ)。

4.根据权利要求1所述的一种基于叠前地震数据的煤层气储层脆性指数反演方法，其特征在于，所述步骤三中的系数矩阵为：

5.根据权利要求1所述的一种基于叠前地震数据的煤层气储层脆性指数反演方法，其特征在于，所述步骤五中的脆性指数计算方程为：