[发明专利]致密油气藏脆性的预测方法及装置

权利要求书

1.一种致密油气藏脆性的预测方法，其特征在于，包括以下步骤：

基于获取的测井曲线得到横波速度与纵波速度曲线、密度曲线、GR曲线、电阻率曲线之间的模型关系，利用所述模型关系再建立横波速度预测模型；

基于所述横波速度预测模型预测得到横波速度曲线；

根据所述横波速度曲线、所述密度曲线和纵波速度曲线计算弹性参数曲线；

基于石英含量测井曲线和叠后地震数据体得到石英含量数据体；

基于所述石英含量数据体得到脆性指数，具体计算公式如下：

其中，B_i表示脆性系指数，E_i表示地层的杨氏模量，E_min表示最小杨氏模量，E_max表示最大杨氏模量，S_i表示岩石的石英含量；

基于叠前地震资料、纵波速度曲线、密度曲线得到时深曲线和地震子波；

基于所述时深曲线、所述密度曲线和所述弹性参数曲线建立初始模型，根据所述初始模型与地震子波褶积得到不同角度的合成角度道集，通过所述合成角度道集与叠前角道集建立反演方程进而反演得到杨氏模量的数据体，基于所述时深曲线、所述密度曲线和所述弹性参数曲线中的杨氏模量曲线和泊松比曲线得到的层位数据，通过插值外推的方法建立杨氏模量、泊松比和密度的初始模型，根据所述初始模型计算得到不同角度的反射系数，再与地震子波褶积得到不同角度的合成角度道集，通过所述合成角度道集与叠前角道集进行比较使两者的误差达到预设条件时建立反演方程，利用最小二乘的反演算法反演得到杨氏模量的数据体；

基于所述石英含量数据体和所述杨氏模量的数据体得到脆性指数的数据体。

2.根据权利要求1所述的致密油气藏脆性的预测方法，其特征在于，在所述基于叠前地震资料、纵波速度曲线、密度曲线得到时深曲线的步骤中具体包括：

基于叠前地震资料得到叠后地震资料；

基于所述叠后地震资料、纵波速度曲线、密度曲线得到时深曲线。

3.根据权利要求1所述的致密油气藏脆性的预测方法，其特征在于，在所述基于获取的测井曲线得到横波速度与纵波速度曲线、密度曲线、GR曲线、电阻率曲线之间的模型关系，利用所述模型关系再建立横波速度预测模型的步骤中，具体为：对密度曲线、GR曲线、电阻率曲线进行环境校正和标准化处理，根据获取的测井曲线得到横波速度，构建所述横波速度和纵波速度曲线、密度曲线、GR曲线、电阻率曲线之间的模型关系，利用所述模型关系建立横波速度预测模型。

4.根据权利要求1所述的致密油气藏脆性的预测方法，其特征在于，在所述根据所述横波速度曲线、所述密度曲线和所述纵波速度曲线计算弹性参数曲线的步骤中，弹性参数曲线中至少包括拉梅常数、剪切模量、体积模量、杨氏模量和泊松比其中之一。

5.根据权利要求1所述的致密油气藏脆性的预测方法，其特征在于，在所述基于石英含量测井曲线和叠后地震数据体得到石英含量数据体的步骤中，具体为利用石英含量测井曲线和叠后地震数据体之间的非线性关系，运用AVF反演得到石英含量数据体。

6.根据权利要求1所述的致密油气藏脆性的预测方法，其特征在于，在所述基于叠前地震资料、纵波速度曲线、密度曲线得到时深曲线和地震子波的步骤中，具体为运用纵波速度曲线和密度曲线，依托由叠前地震资料获得的叠后地震资料进行标定，使得实际井旁道和合成地震记录达到匹配，此时的时深曲线和地震子波为所求的时深曲线和地震子波。

7.根据权利要求1所述的致密油气藏脆性的预测方法，其特征在于，在所述根据所述初始模型计算得到不同角度的反射系数的步骤中，具体计算公式如下：

经转化得到，

其中，θ表示入射角，k表示横纵波速度比，E、σ、ρ分别表示杨氏模量的数据体、泊松比的数据体和裂缝密度的数据体，R_pp表示反射系数。