[发明专利]一种硫沉积储层伤害时空分布特征预测方法

权利要求书

1.一种硫沉积储层伤害时空分布特征预测方法，其特征在于，包括以下步骤：

获得气井控制储层的单质硫流体相图，并从所述单质硫流体相图获取所述气井控制储层温度条件下的硫析出压力P_s及析出状态；

测量所述气井控制储层的气样硫溶解量和地层饱和单质硫含量；

基于所述单质硫流体相图和所述气井控制储层的气样硫溶解量和地层饱和单质硫含量，建立所述气井控制储层的元素硫溶解模型；

通过单质硫岩心驱替实验获得单质硫的硫饱和度S_s，并根据岩芯驱替实验的参数，获得不同参数条件下获得的硫饱和度S_s对应的气相相对渗透率K_rg；

基于封闭定容气藏物质平衡方程获得所述气井控制储层的平均地层压力P与所述气井控制储层的气井累积产气量Q的关系；

基于元素硫溶解模型、不同参数条件下获得的硫饱和度S_s对应的气相相对渗透率K_rg和所述气井控制储层的平均地层压力P与所述气井控制储层的气井累积产气量G_p的关系，建立所述气井控制储层的硫沉积储层伤害时空变化规律预测模型；

根据所述硫沉积储层伤害时空变化规律预测模型预测所述气井控制储层出现硫沉积的时间和位置、硫沉积程度。

2.根据权利要求1所述的硫沉积储层伤害时空分布特征预测方法，其特征在于，所述气井控制储层的单质硫流体相图通过单质硫在所述气井的地层条件下的单质硫析出曲线、单质硫凝固曲线、单质硫露点线、单质硫泡点线获得。

3.根据权利要求2所述的硫沉积储层伤害时空分布特征预测方法，其特征在于，所述元素硫溶解模型为：

式中：T—单质硫流体温度，K；

ρ—单质硫流体密度，kg/m³；

C_r—硫溶解度，g/m³；

K，M，N—拟合回归常数；

P—平均地层压力，MPa；

P_s—气井控制储层温度条件下元素硫开始析出压力，MPa；

C_i—气井控制储层条件下硫组分初始浓度，g/m³。

4.根据权利要求1所述的硫沉积储层伤害时空分布特征预测方法，其特征在于，所述岩芯驱替实验的参数至少包括岩心驱替压差、岩心长度、气体粘度、气体渗流速度、岩芯孔隙度中的任一参数。

5.根据权利要求4所述的硫沉积储层伤害时空分布特征预测方法，其特征在于，所述硫饱和度S_s与气相相对渗透率K_rg的对应关系满足以下表达式：

K_rg＝k_rgiexp(αS_s) (2)

式中：α—回归参数；

k_rgi—气相相对渗透率最大端点值。

6.根据权利要求5所述的硫沉积储层伤害时空分布特征预测方法，其特征在于，所述硫饱和度S_s与所述岩芯驱替实验的参数的关系式为：

其中，p表示井筒内压力，MPa；q表示产气量，m³/d；h表示储层厚度，m；r表示径向半径，即气井控制半径以内的储层到气井井底的距离，m；μ表示流体黏度，mPa.s；k表示储层渗透率，m²；表示受紊流影响的系数，m^-1；ρ表示流体密度，kg/m³；γ_g表示气体相对密度，无量纲；R表示摩尔气体常数，无量纲；Z表示气体偏差因子，无量纲；T表示气体绝对温度，K；q_sc表示标准状态下的气井产气量；s_wi表示束缚水饱和度，无量纲；ρ_s表示液相硫的密度，kg/m³；B_g表示气体体积系数，无量纲。