[发明专利]一种应力敏感储层产水气井动态反演方法

权利要求书

1.一种应力敏感储层产水气井动态反演方法，其特征在于，包括如下步骤：

(1)、日常生产数据或试井测试数据、压裂施工及测井数据、室内实验测试数据相关基础数据的整理及筛选；

(2)、基于以下关系式，计算考虑相渗曲线应力敏感、绝对渗透率应力敏感及动态滑脱效应等气-水两相复杂渗流特征的改进拟压力及改进拟时间或简化拟压力及简化拟时间，其中，简化的拟参数在考虑气-水两相复杂渗流特征的同时结合孔隙水产水特征对拟参数进行了相应的简化；

其中，ψ_two为考虑应力敏感储层气-水两相复杂渗流特征的改进的拟压力MPa/cp；t_two为考虑应力敏感储层气-水两相复杂渗流特征的改进的拟时间d；ψ_two-g为考虑应力敏感储层气-水两相复杂渗流特征及孔隙水产水特征的简化拟压力MPa²/cp；t_two-g为考虑应力敏感储层气-水两相复杂渗流特征的简化的拟时间d；p为压力MPa；p_a为参考压力MPa；ρ_g为气相地下密度kg/m³；t为真实时间d；t_a为参考时间d；ρ_w为水的地下密度kg/m³；ρ_wsc为水的地面标准密度kg/m³；为平均压力所对应的气相地下密度kg/m³；为平均压力所对应的水相地面标准密度kg/m³；k_rgE为气相等效相对渗透率；k_rwE为水相等效相对渗透率；为平均压力所对应的气相等效相对渗透率；为平均压力所对应的水相等效相对渗透率；μ_g为气体粘度cp；μ_w为水的粘度cp；μ_gi为原始地层压力下的气体粘度cp；为平均压力所对应的气体粘度cp；为平均压力所对应的水的粘度cp；C_t-twoi为原始地层压力下的气水两相综合压缩系数MPa^-1；为平均压力所对应的气水两相综合压缩系数MPa^-1；z_g为天然气压缩因子C_t-two-gi、分别为原始地层压力及地层平均压力所对应的简化气水两相综合压缩系数MPa^-1；C_t-two-g＝C_p+S_gC_g+S_wC_w，C_t-two-g为简化的气水两相综合压缩系数MPa^-1；k_i为储层原始渗透率mD；

(3)、利用考虑了气-水两相复杂渗流特征的流动阶段分析方程绘制相应的分析曲线对流动阶段进行识别；

(4)、利用相应的流动阶段分析曲线，结合参数与曲线之间的定量关系，初步确定相应的储层及裂缝参数；

(5)、利用考虑复杂渗流特征的数值模型对实际双对数曲线进行拟合，从而验证流动阶段分析方法所获得参数的正确性，若拟合效果较好，可进入下一步骤，若拟合效果较差，则返回步骤(4)重新确定储层及裂缝参数初始值；

(6)、利用考虑复杂渗流特征的数值模型，以实际井底流压作为输入工作制度，从而模拟得到产气量及产水量数据，并与实际数据在直角坐标系下进行生产动态数据历史拟合，从而进一步验证以上步骤所确定参数的准确性，若拟合效果较好，则说明步骤(4)～(5)所得到的储层及裂缝参数较为准确，否则需要对以上步骤的分析过程进行检查，直至得到较好的拟合效果。

2.根据权利要求1所述的应力敏感储层产水气井动态反演方法，其特征在于：在步骤(1)中需要收集整理室内试验数据包括：不同压力条件下的气-水相对渗透率曲线或可用于计算该曲线的基础数据。

3.根据权利要求1所述的应力敏感储层产水气井动态反演方法，其特征在于：在步骤(5)及(6)中的数值模型综合考虑了相渗曲线应力敏感、储层绝对渗透率应力敏感及动态气体滑脱效应等气-水两相复杂渗流特征，其中气相及水相的流速方程分别为：

其中，v_g为气体流速m/s；v_w为水相流速m/s；p为压力MPa；k<p>表示为压力函数的储层绝对渗透率md；μ_g<p>表示为压力函数的气体粘度cp；μ_w<p>表示为压力函数的水相粘度cp；b<p,S_w>表示同时为压力及饱和度函数的气体滑脱因子MPa；k_rg<p,S_w>表示同时为压力及饱和度函数的气相相对渗透率；k_rw<p,S_w>表示同时为压力及饱和度函数的水相相对渗透率。