[发明专利]基于枯竭气藏储气库地应力反演的临界出砂压差预测方法

权利要求书

1.一种基于枯竭气藏储气库地应力反演的临界出砂压差预测方法，其特征在于，包括：

获取不同深度的内聚力C和内摩擦角φ；

基于测井曲线、所述内聚力C和所述内摩擦角φ，建立单井弹性力学参数剖面，以初始应力场结果为基础，开展三维应力场数值模拟，结合不同时间点的应力场分布特征，进行地应力反演，得到不同深度的三向地应力反演结果；

通过公式计算得到x、y、z方向的正应力分量σ_x、σ_y、σ_z和xy面、yz面、xz面上的应力分量σ_xy、σ_yz、σ_xz；其中，S_H为最大水平地应力，S_h为最小水平地应力，S_v为垂向地应力，β为方位角，θ为井斜角；

通过公式计算得到沿井周柱坐标系的径向应力σ_R，沿井周柱坐标系的切向应力σ_α，沿井周柱坐标系的轴向应力σ_z’，井周的切应力τ_αz；其中，P_w为井底流压，σ_z为z方向的正应力分量，a为井周角，μ为泊松比；

通过公式计算得到柱坐标R方向的有效主应力σ_R’、除R向外的另外两个有效主应力σ_1m和σ_2m；其中，ε为有效应力系数，P_p为孔隙压力；

通过公式计算得到井周的最大主应力σ₁、中间主应力σ₂及最小主应力σ₃；

通过公式计算得到井壁骨架岩石失稳指数F；其中，σ_m，2为正应力平均值，τ_otc为八面体剪切应力；

当F＝0时，通过公式CDP＝P_p-P_w计算得到临界生产压差CDP。

2.如权利要求1所述的基于枯竭气藏储气库地应力反演的临界出砂压差预测方法，其特征在于，所述获取不同深度的内聚力C和内摩擦角φ，包括：

取芯制样后，分别进行室内常规三轴压缩实验、不同原始含水饱和状态下的剪切特性实验和周期性循环加载岩石力学实验，得到不同沉积微相的内聚力C和内摩擦角φ。

3.如权利要求1或2所述的基于枯竭气藏储气库地应力反演的临界出砂压差预测方法，其特征在于，在所述计算得到临界生产压差CDP之后，还包括：

建立枯竭气藏储气库储层不同沉积微相临界生产压差表，同时形成基于测井数据的单井不同沉积微相临界出砂压差剖面分布模式图。