[发明专利]一种变应力场中非常规储层气-水相渗曲线计算方法

权利要求书

1.一种变应力场中非常规储层气-水相渗曲线计算方法，其特征是：具有如下步骤：

(a)、首先利用实验方法测定5个以上不同应力条件(需包含原始地层压力所对应的应力条件)下的孔隙度及相渗曲线各端点值，并将实验室应力条件折算为所对应的地层压力；

(b)、基于原始地层压力所对应的应力条件下的实验结果，拟合以下关系式获得毛管分布指数λ及毛管弯曲系数η；

其中，k_rgi为原始地层压力所对应的应力条件下的气相相对渗透率；k_rgendi为原始地层压力下的气相相渗端点值；S_w为含水饱和度；S_wci为原始地层压力下的束缚水饱和度；S_gri为原始地层压力下的残余气饱和度；λ为毛管分布指数；η为毛管弯曲系数；k_rwi为原始地层压力所对应的应力条件下的水相相对渗透率；k_rwendi为原始地层压力下的水相相渗端点值；

(c)、基于不同应力条件下的实验数据，拟合孔隙度与地层压力变化的关系得到孔隙度应力敏感系数γ；并基于实验结果，分别根据以下关系式拟合获得相渗曲线各端点值与应力之间的定量表征系数C、D、E、F：

利用所获得的相渗曲线各端点值与应力之间的定量表征系数，代入以下关系式，从而构建相渗曲线端点值随应力场的动态变化模型：

其中，φ为指定应力条件下的储层孔隙度；φ_i为原始地层压力所对应的应力条件下的储层孔隙度；k_rgend-test、k_rwend-test、S_wc-test、S_gr-test分别为实验压力下的气相相渗端点值、水相相渗端点值、束缚水饱和度、残余气饱和度测试结果；k_rgendi、k_rwendi、S_wci、S_gri分别为原始地层压力所对应的应力条件下的气相相渗端点值、水相相渗端点值、束缚水饱和度、残余气饱和度；k_rgend、k_rwend、S_wc、S_gr分别为指定地层压力条件下的气相相渗端点值、水相相渗端点值、束缚水饱和度、残余气饱和度；γ为孔隙度应力敏感系数，单位为MPa^-1；C为气相相渗端点值应力敏感系数；D为水相相渗端点值应力敏感系数；E为束缚水饱和度应力敏感系数；F为残余气饱和度应力敏感系数；p为指定地层压力，单位为MPa；p_i为原始地层压力，单位为MPa；

(d)、利用所构建的相渗曲线端点值随应力场的动态变化模型，进一步建立考虑相渗曲线端点值随应力变化特征的气-水相对渗透率计算模型，并计算不同应力条件下的气-水两相相对渗透率；

(e)、利用计算获得的不同应力条件下的气-水两相相对渗透率结果，绘制变应力场中的气-水相对渗透率曲线。

2.如权利要求1所述的变应力场中非常规储层气-水相渗曲线计算方法，其特征是：所述步骤(d)包括以下步骤：

(d1)、基于孔隙度应力敏感系数γ、气相相渗端点值应力敏感系数C，水相相渗端点值应力敏感系数D，利用下式计算不同压力条件下的气-水相对渗透率端点值；

(d2)、基于孔隙度应力敏感系数γ、束缚水饱和度应力敏感系数E及残余气饱和度应力敏感系数F，利用下式计算不同压力条件下的标准化水相饱和度；

其中，为标准化水相饱和度；

(d3)、基于标准化水相饱和度、孔隙度应力敏感系数及相渗曲线各端点值的应力敏感系数，利用下式，计算不同压力及饱和度条件下的气-水相对渗透率；

其中，k_rg为气相相对渗透率；k_rw为水相相对渗透率。