[发明专利]一种砂砾岩储层纵向连续脆性指数预测方法

摘要

本发明属于砂砾岩油气储层水力压裂技术领域，涉及一种砂砾岩储层纵向连续脆性指数预测方法。步骤如下：一、根据测井解释数据获取砂砾岩储层纵向连续的动态力学参数、孔隙度和砾石粒径及含量；二、建立弹性参数动静态转化关系；三、建立砂砾岩储层强度参数‑测井数据之间的解释模型；四、建立数值岩心模型，从而建立数值岩心宏细观参数转化关系；五、建立不同纵向深度的储层数值岩心并开展单轴压缩数值试验，获取应力、应变数据和全应力应变曲线；六、结合静态弹性模量计算脆性指数。本发明实现了对全井段砂砾岩储层岩心的数值化表征，可为砂砾岩储层可压性评价和射孔压裂层段优选提供指导，有助于促进油气藏经济高效开发。

主权项

1.一种砂砾岩储层纵向连续脆性指数预测方法，其特征在于，包括如下步骤：/n第一步：根据测井解释数据获取砂砾岩储层纵向连续的动态力学参数和孔隙度；/n第二步：开展岩心单轴压缩测试，获取纵向上非连续的静态力学参数，结合由测井数据得到的动态弹性参数和由室内测试得到的静态弹性参数，采用最小二乘回归方法建立砂砾岩储层弹性参数的动静态转化关系，见公式(1)和(2)，由此获得纵向连续的静态弹性参数分布；/nE_s＝f₁(E_d) (1)/nν_s＝f₂(ν_d) (2)/n式中：E_s和E_d分别为静态和动态弹性模量，GPa，ν_s和ν_d别为静态和动态泊松比，无量纲；/n第三步：利用砂砾岩岩心静态强度参数对由测井解释得到的动态强度参数进行归位校核，建立砂砾岩储层强度参数-测井数据之间的解释模型，见公式(3)和(4)；/nσ_c＝F(C₁) (3)/nθ＝F(C₂,Φ) (4)/n式中：σ_c为校核后的单轴抗压强度；θ为校核后的残余强度系数，θ＝σ_r/σ_c，σ_r为残余强度；Φ为测井解释得到的岩心孔隙度；C₁和C₂分别为不同类型的测井数据；/n利用砂砾岩储层弹性参数的动静态转化关系和强度参数-测井数据解释模型，对测井解释力学参数进行校核，得到砂砾岩储层纵向连续的精细的力学参数表征；/n第四步：根据真实物理岩心的尺寸以及砾石含量和粒径，建立储层岩心数值有限元模型，对模型划分网格和施加轴向位移荷载并开展多次单轴压缩数值试验；利用物理岩心的力学参数和应力应变曲线对相应数值岩心基质的细观力学参数进行标定，保持砾石和界面材料的细观力学参数不变，从而建立数值岩心宏细观参数转化关系，见公式(6)和(7)；/n /n /n式中：σ_c和σ_c0分别为岩心宏观和细观单轴抗压强度，MPa；E_s和E_s0分别为Weibull分布赋值时的岩心宏观和细观弹性模量均值，单位GPa；m为均值度系数，用以表征岩心力学参数的非均质性，m值越小则非均质性越强；/n第五步：根据步骤第二步、第三步中得到的岩心静态宏观力学参数和第四步中得到的宏细观力学参数转化关系，可对任意层位储层对应的数值岩心细观力学参数进行标定，从而建立不同纵向深度的储层数值岩心并开展单轴压缩数值试验，获取应力、应变数据和全应力应变曲线；/n第六步：根据不同纵向深度储层对应的数值岩心应力应变曲线，获取对应的脆性破坏残余强度特征参数σ_r和ε_r，进而根据式(8)计算脆性指数BI，由此得到纵向连续的脆性指数预测值；/n /n式中，σ_r为残余强度，ε_r为残余应变。/n