[发明专利]裂缝流体因子预测方法及系统

权利要求书

1.一种裂缝流体因子预测方法，其特征在于，包括：

计算各向异性截距阻抗和梯度阻抗数据体，计算HTI介质中的零阶扩展方位傅里叶系数，确定预测流体因子的最佳旋转角，通过线性拟合得到HTI介质中的零阶扩展方位傅里叶系数与裂缝流体因子的关系式，结合HTI介质中的零阶扩展方位傅里叶系数与裂缝流体因子的关系式，得到裂缝流体因子；

所述HTI介质中的零阶扩展方位傅里叶系数为：

其中，AI_iso＝αρ、分别表示各向同性截距阻抗和梯度阻抗；AI₀＝α₀ρ₀表示纵波阻抗均值；α、β、ρ分别表示各向同性背景介质的纵、横波速度及密度，α₀、β₀、ρ₀分别表示各向同性背景介质的纵、横波速度及密度的平均值；g＝β²/α²表示横波速度平方与纵波速度平方之比；k＝1-2g；δ_N、δ_T分别表示法向及切向弱度；χ表示旋转角。

2.根据权利要求1所述的裂缝流体因子预测方法，其特征在于，基于方位部分角度叠加地震数据、方位角度地震子波以及测井曲线插值得到各向异性截距阻抗和梯度阻抗低频模型，通过分方位叠前地震贝叶斯反演得到各向异性截距阻抗和梯度阻抗数据体。

3.根据权利要求2所述的裂缝流体因子预测方法，其特征在于，利用各向异性截距阻抗和梯度阻抗反演结果，计算不同旋转角下的扩展方位弹性阻抗，对扩展方位弹性阻抗进行傅里叶级数展开，得到HTI介质中的零阶扩展方位傅里叶系数。

4.根据权利要求1所述的裂缝流体因子预测方法，其特征在于，计算HTI介质中的零阶扩展方位傅里叶系数与裂缝流体因子的相关系数，确定预测流体因子的最佳的旋转角。

5.根据权利要求4所述的裂缝流体因子预测方法，其特征在于，结合各向异性截距阻抗和梯度阻抗数据体以及最佳旋转角，计算得到扩展方位弹性阻抗数据体，对扩展方位弹性阻抗数据体进行傅里叶级数展开，得到HTI介质中的零阶扩展方位傅里叶系数数据体，利用HTI介质中的零阶扩展方位傅里叶系数与裂缝流体因子的关系式，将HTI介质中的零阶扩展方位傅里叶系数数据体转化为裂缝流体因子。

6.根据权利要求2所述的裂缝流体因子预测方法，其特征在于，利用测井曲线插值得到各向异性截距阻抗和梯度阻抗低频模型。

7.一种基于如权利要求1-6任一项所述的裂缝流体因子预测方法的裂缝流体因子预测系统，其特征在于，包括：

第一计算模块，用于计算各向异性截距阻抗和梯度阻抗数据体；

第二计算模块，用于计算HTI介质中的零阶扩展方位傅里叶系数；

确定模块，用于确定预测流体因子的最佳旋转角；

拟合模块，用于将HTI介质中的零阶扩展方位傅里叶系数与裂缝流体因子进行线性拟合，得到HTI介质中的零阶扩展方位傅里叶系数与裂缝流体因子的关系式；

转化模块，用于结合各向异性截距阻抗和梯度阻抗数据体、最佳旋转角以及HTI介质中的零阶扩展方位傅里叶系数与裂缝流体因子的关系式，将HTI介质中的零阶扩展方位傅里叶系数转化为裂缝流体因子；

所述HTI介质中的零阶扩展方位傅里叶系数为：

其中，AI_iso＝αρ、分别表示各向同性截距阻抗和梯度阻抗；AI₀＝α₀ρ₀表示纵波阻抗均值；α、β、ρ分别表示各向同性背景介质的纵、横波速度及密度，α₀、β₀、ρ₀分别表示各向同性背景介质的纵、横波速度及密度的平均值；g＝β²/α²表示横波速度平方与纵波速度平方之比；k＝1-2g；δ_N、δ_T分别表示法向及切向弱度；χ表示旋转角。