CN113107474B - 利用短期关井数据预测低渗气藏的储层参数的方法及装置 - Google Patents

Abstract

本公开提供了一种利用短期关井数据预测低渗气藏的储层参数的方法及装置，属于低渗气藏开发技术领域。本公开通过获取气藏的地层压力、渗透率、井底静压以及关井时间之间的对应关系；基于关井后的较短时间内的多个关井时间对应的井底静压，获取的该气井的井底静压随关井时间变化的函数关系；基于上述对应关系和函数关系，计算得到关井时间较长时的储层参数，避免了为了获取储层参数而进行长时间的关井，提高了气藏的经济效益。

Description

利用短期关井数据预测低渗气藏的储层参数的方法及装置

技术领域

本公开涉及低渗气藏开发技术领域，特别涉及一种利用短期关井数据预测低渗气藏的储层参数的方法及装置。

背景技术

在低渗气藏的开发过程中，准确获取生产井目前的储层参数，对于生产井后期的产能规划十分重要，其中，储层参数包括渗透率及地层压力等。

目前通常使用的利用短期关井数据预测低渗气藏的储层参数的方法包括：在关井一段时间后，井底压力恢复到真实的地层压力时获取，关井时间越长，获取到的储层参数越准确，但是长时间的关井，不能正常生产，极大地影响了气藏的经济效益。

发明内容

本公开实施例提供了一种利用短期关井数据预测低渗气藏的储层参数的方法及装置，能够解决目前通常使用的利用短期关井数据预测低渗气藏的储层参数的方法需要进行长时间的关井，影响生产的问题。该技术方案如下：

一方面，提供了一种利用短期关井数据预测低渗气藏的储层参数的方法，该方法包括：

基于平面径向流渗流数学模型以及气藏特性，采用假想注气井法，获取该气藏的地层压力、渗透率、井底静压以及关井时间之间的对应关系；

在该气藏中的气井关井后，获取多个关井时间对应的井口油压，每个该关井时间小于预设时长；

基于该关井时间及相应的井口油压，获取该多个关井时间的井底静压；

基于该多个关井时间的井底静压，确定该气井的井底静压随关井时间变化的函数关系；

基于该对应关系以及该气井的井底静压随关井时间的函数关系，获取该气藏的储层参数，该储层参数包括：渗透率以及地层压力。

在一种可能实现方式中，该基于该关井时间及相应的井口油压，获取该多个关井时间的井底静压，包括：

基于该气井的气层中部垂深、井筒静气柱平均温度、稳定生产产量以及气井关井前累积生产时间，该气井生产的气体在储层条件下的相对密度、井筒静气柱平均偏差因子以及气体粘度，在该气井关井后，每隔预设时长获取到的关井时间及相应的井口油压，以及下述公式(1)，获取每个该关井时间对应的井底静压：

式中，P_w表示关井时间为Δt时的井底静压，MPa；P_wh表示关井时间为Δt时的井口油压，MPa；γ_g表示气体的相对密度，无因次；H表示气井的气层中部垂深，m；表示井筒静气柱平均温度，K；/>表示井筒静气柱平均偏差因子，无因次。

在一种可能实现方式中，该基于该多个关井时间的井底静压，确定该气井的井底静压随关井时间变化的函数关系，包括：

基于每个该关井时间以及对应的井底静压，以为横坐标，以/>为纵坐标，绘制关系曲线，其中，Δt表示该关井时间，t_p表示气井关井前累积生产时间，P_w表示关井时间为Δt时的井底静压；

截取该关系曲线中的直线段；

对该直线段进行拟合，得到与/>之间的函数关系。

在一种可能实现方式中，基于下述公式(2)，获取该渗透率：

式中，k表示渗透率，μm²；q_sc表示气井稳定生产产量，10⁴m³/d；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；μ表示在储层条件下的气体粘度，mPa·s；m为该直线段的斜率的绝对值，无因次。

在一种可能实现方式中，基于下述公式(3)，获取该地层压力：

式中，P_r表示关井时间为Δt时的地层压力，MPa；P_w表示关井时间为Δt时的井底静压，MPa；q_sc表示气井稳定生产产量，10⁴m³/d；μ表示在储层条件下的气体粘度，mPa·s；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；k表示气藏渗透率，μm²；h表示气层有效厚度，m；t_p表示气井关井前累积生产时间，h；Δt表示关井时间，h；

当该等于1时，将所得的P_r的值获取为该地层压力。

一方面，提供了一种利用短期关井数据预测低渗气藏的储层参数的装置，该装置包括：

对应关系获取模块，用于基于平面径向流渗流数学模型以及气藏特性，采用假想注气井法，获取该气藏的地层压力、渗透率、井底静压以及关井时间之间的对应关系；

井口油压获取模块，用于在该气藏中的气井关井后，获取多个关井时间对应的井口油压，每个该关井时间小于预设时长；

井底静压获取模块，用于基于该关井时间及相应的井口油压，获取该多个关井时间的井底静压；

函数关系获取模块，用于基于该多个关井时间的井底静压，确定该气井的井底静压随关井时间变化的函数关系；

储层参数获取模块，用于基于该对应关系以及该气井的井底静压随关井时间的函数关系，获取该气藏的储层参数，该储层参数包括：渗透率以及地层压力。

在一种可能实现方式中，该井底静压获取模块，用于基于该气井的气层中部垂深、井筒静气柱平均温度、稳定生产产量以及气井关井前累积生产时间，该气井生产的气体在储层条件下的相对密度、井筒静气柱平均偏差因子以及气体粘度，在该气井关井后，每隔预设时长获取到的关井时间及相应的井口油压，以及下述公式(1)，获取每个该关井时间对应的井底静压：

式中，P_w表示关井时间为Δt时的井底静压，MPa；P_wh表示关井时间为Δt时的井口油压，MPa；γ_g表示气体的相对密度，无因次；H表示气井的气层中部垂深，m；表示井筒静气柱平均温度，K；/>表示井筒静气柱平均偏差因子，无因次；

函数关系获取模块，用于基于每个该关井时间以及对应的井底静压，拟合得到与/>之间的函数关系；

储层参数获取模块，用于基于该对应关系以及该气井的井底静压随关井时间的函数关系，获取该气藏的储层参数，该储层参数包括：渗透率以及地层压力。

在一种可能实现方式中，该函数关系获取模块，用于：

基于每个该关井时间以及对应的井底静压，以为横坐标，以/>为纵坐标，绘制关系曲线，其中，Δt表示该关井时间，t_p表示气井关井前累积生产时间，P_w表示关井时间为Δt时的井底静压；

截取该关系曲线中的直线段；

对该直线段进行拟合，得到与/>之间的函数关系。

在一种可能实现方式中，基于下述公式(2)，获取该渗透率：

式中，k表示渗透率，μm²；q_sc表示气井稳定生产产量，10⁴m³/d；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；μ表示在储层条件下的气体粘度，mPa·s；m为该直线段的斜率的绝对值，无因次。

在一种可能实现方式中，基于下述公式(3)，获取该地层压力：

式中，P_r表示关井时间为Δt时的地层压力，MPa；P_w表示关井时间为Δt时的井底静压，MPa；q_sc表示气井稳定生产产量，10⁴m³/d；μ表示在储层条件下的气体粘度，mPa·s；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；k表示气藏渗透率，μm²；h表示气层有效厚度，m；t_p表示气井关井前累积生产时间，h；Δt表示关井时间，h；

当该等于1时，将所得的P_r的值获取为该地层压力。

本公开实施例提供的方法，通过获取气藏的地层压力、渗透率、井底静压以及关井时间之间的对应关系；基于关井后的较短时间内的多个关井时间对应的井底静压，获取的该气井的井底静压随关井时间变化的函数关系；基于上述对应关系和函数关系，计算得到关井时间较长时的储层参数，避免了为了获取储层参数而进行长时间的关井，提高了气藏的经济效益。

附图说明

为了更清楚地说明本公开实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本公开的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本公开实施例提供的一种利用短期关井数据预测低渗气藏的储层参数的方法的流程图；

图2是本公开实施例提供的一种利用短期关井数据预测低渗气藏的储层参数的方法的流程图；

图3是本公开实施例提供的一种关井时间和井底静压的关系曲线图；

图4是本公开实施例提供的一种与Δt的关系曲线图；

图5是本公开实施例提供的一种试井解释图；

图6是本公开实施例提供的一种利用短期关井数据预测低渗气藏的储层参数的装置的结构示意图；

图7是本公开实施例提供的一种计算机设备的结构示意图。

具体实施方式

为使本公开的目的、技术方案和优点更加清楚，下面将结合附图对本公开实施方式作进一步地详细描述。

图1是本公开实施例提供的一种利用短期关井数据预测低渗气藏的储层参数的方法的流程图，请参见图1，该方法可以应用于计算机设备，该方法包括：

101、基于平面径向流渗流数学模型以及气藏特性，采用假想注气井法，获取该气藏的地层压力、渗透率、井底静压以及关井时间之间的对应关系。

102、在该气藏中的气井关井后，获取多个关井时间对应的井口油压，每个该关井时间小于预设时长。

103、基于该关井时间及相应的井口油压，获取该多个关井时间的井底静压。

104、基于该多个关井时间的井底静压，确定该气井的井底静压随关井时间变化的函数关系。

105、基于该对应关系以及该气井的井底静压随关井时间的函数关系，获取该气藏的储层参数，该储层参数包括：渗透率以及地层压力。

本公开实施例提供的方法，通过获取气藏的地层压力、渗透率、井底静压以及关井时间之间的对应关系；基于关井后的较短时间内的多个关井时间对应的井底静压，获取的该气井的井底静压随关井时间变化的函数关系；基于上述对应关系和函数关系，计算得到关井时间较长时的储层参数，避免了为了获取储层参数而进行长时间的关井，提高了气藏的经济效益。

在一种可能实现方式中，该基于该关井时间及相应的井口油压，获取该多个关井时间的井底静压，包括：

基于该气井的气层中部垂深、井筒静气柱平均温度、稳定生产产量以及气井关井前累积生产时间，该气井生产的气体在储层条件下的相对密度、井筒静气柱平均偏差因子以及气体粘度，在该气井关井后，每隔预设时长获取到的关井时间及相应的井口油压，以及下述公式(1)，获取每个该关井时间对应的井底静压：

式中，P_w表示关井时间为Δt时的井底静压，MPa；P_wh表示关井时间为Δt时的井口油压，MPa；γ_g表示气体的相对密度，无因次；H表示气井的气层中部垂深，m；表示井筒静气柱平均温度，K；/>表示井筒静气柱平均偏差因子，无因次。

在一种可能实现方式中，该基于该多个关井时间的井底静压，确定该气井的井底静压随关井时间变化的函数关系，包括：

基于每个该关井时间以及对应的井底静压，以为横坐标，以/>为纵坐标，绘制关系曲线，其中，Δt表示该关井时间，t_p表示气井关井前累积生产时间，P_w表示关井时间为Δt时的井底静压；

截取该关系曲线中的直线段；

对该直线段进行拟合，得到与/>之间的函数关系。

在一种可能实现方式中，基于下述公式(2)，获取该渗透率：

式中，k表示渗透率，μm²；q_sc表示气井稳定生产产量，10⁴m³/d；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；μ表示在储层条件下的气体粘度，mPa·s；m为该直线段的斜率的绝对值，无因次。

在一种可能实现方式中，基于下述公式(3)，获取该地层压力：

式中，P_r表示关井时间为Δt时的地层压力，MPa；P_w表示关井时间为Δt时的井底静压，MPa；q_sc表示气井稳定生产产量，10⁴m³/d；μ表示在储层条件下的气体粘度，mPa·s；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；k表示气藏渗透率，μm²；h表示气层有效厚度，m；t_p表示气井关井前累积生产时间，h；Δt表示关井时间，h。

当该等于1时，将所得的P_r的值获取为该地层压力。

上述所有可选技术方案，可以采用任意结合形成本公开的可选实施例，在此不再一一赘述。

图2是本公开实施例提供的一种利用短期关井数据预测低渗气藏的储层参数的方法的流程图，请参见图2，该方法可以应用于计算机设备，该方法包括：

201、基于平面径向流渗流数学模型以及气藏特性，采用假想注气井法，获取该气藏的地层压力、渗透率、井底静压以及关井时间之间的对应关系。

具体地，若气层均质等厚，内边界有一口定产量气井，外边界无限大，原始地层条件下的流体处于压力静平衡状态。当气井投产后，地层中的流体流动为平面径向流。例如，平面径向流渗流数学模型可以是如公式(4)所示；

式中，ψ表示地层拟压力，MPa²/(mPa·s)；

ψ_i表示原始条件下的地层拟压力，MPa²/(mPa·s)；

r表示径向断面半径，m；

μ表示气层中部压力和温度下的气体粘度，mPa·s；

Z表示气层中部压力和温度下的气体偏差因子；

T表示气层中部温度，K；

k表示气藏渗透率，μm²；

h表示气层有效厚度，m；

q_sc表示气井稳定生产产量，10⁴m³/d；

p_sc表示地面标准压力，0.10MPa；

T_sc表示地面标准温度，为293.15K；

Z_sc表示地面标准状态下的气体偏差因子，无因次。

其中，该气藏具有气压较低等特性。在涉及压力恢复问题的计算过程中，一般是通过原生产井点处钻一口假想注气井来实现的，注气井的流量与生产井相同，但符号相反。压力恢复时，可认为生产井并未关井，以产量q_sc继续生产(t_p+Δt)时间，而注气井从t_p开始注气，注气时间为Δt。气藏中的产量为两口井的产量代数和，井点处地层压力降落为两口井单独生产时所产生的压力降落之和。

例如，该对应关系可以是如公式(3)所示：

式中，P_r表示关井时间为Δt时的地层压力，MPa；P_w表示关井时间为Δt时的井底静压，MPa；q_sc表示气井稳定生产产量，10⁴m³/d；μ表示在储层条件下的气体粘度，mPa·s；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；k表示气藏渗透率，μm²；h表示气层有效厚度，m；t_p表示气井关井前累积生产时间，h；Δt表示关井时间，h。

202、在该气藏中的气井关井后，获取多个关井时间对应的井口油压，每个该关井时间小于预设时长。

具体地，可以是每隔一段时间获取井口油压，形成多个关井时间对应的井底油压。可以根据气井的生产情况设定该预设时长，该预设时长为一个相对较短的时间。从而可以保证该气井无需停井过长时间，即可恢复生产，保证了气井的经济效益。例如，该预设时长可以为10天。

例如，大深XX井是四川盆地大兴场潜伏构造西南翼的一口开发井，生产层位为茅口组，本实施例以该井为例进行了现场分析，并和同一关井时间的压力恢复试井解释成果进行了对比验证。对于该井，可以是在连续生产155天后，开始进行压力恢复测试，关井10天，记录关井时间与关井后的井口油压数据如下表1所示。

表1

关井时间Δt(h)	井口油压Pwh(MPa)
		24	33
48	35
		72	36.2
96	36.7
		120	37.2
144	37.5
		168	37.8
192	37.9
		216	38.1
240	38.2

基于上述表1可得，关井10天后，井口油压仍然保持上涨趋势。

203、获取气井的各项生产参数。

具体地，获取气井的气层中部垂深、井筒静气柱平均温度、稳定生产产量以及气井关井前累积生产时间，该气井生产的气体在储层条件下的相对密度、井筒静气柱平均偏差因子以及气体粘度。

上述参数可以是在气井的关井后短期内获取到的，基于上述数据进行后续的计算过程，使得关井时间较短，对生产影响较小。具体地，可以根据生产需要和计算需要设计生产时间和关井时间。对于上述大深XX井，生产层位储层具有典型的低孔低渗特征，平均孔隙度在1％左右，测井及岩心分析的渗透率在0.01～3.08mD之间，储层温度为417.9K，储层厚度在15～20m。该井的日产气量为10*10⁴m³。

204、获取每个关井时间对应的井底静压。

具体地，基于步骤202和步骤203中获取的数据以及下述公式(1)，获取每个该关井时间对应的井底静压：

式中，P_w表示关井时间为Δt时的井底静压，MPa；P_wh表示关井时间为Δt时的井口油压，MPa；γ_g表示气体的相对密度，无因次；H表示气井的气层中部垂深，m；表示井筒静气柱平均温度，K；/>表示井筒静气柱平均偏差因子，无因次。

进一步地，可以基于上述数据绘制关井时间和井底静压的关系曲线，以便获取井底静压随关井时间的变化规律。

例如，该大深XX井的关井时间和井底静压的关系曲线如图3所示，图3是本公开实施例提供的一种关井时间和井底静压的关系曲线图，在图3中，横坐标为关井时间，h，纵坐标为井底静压，MPa。由图3可知，井底静压随关井时间的延长而增加，基于上述方式获取的井底静压更为准确，相比于仅仅使用关井末期的油压数据进行井底静压的转换计算，本实施例得到的井底静压为连续的多个数据，基于此，能够获取更为准确的储层参数。

上述步骤获取到的井底静压为一系列的数据点，基于该数据点，可以进行一系列的拟合，进而获取到关井后井底静压的变化趋势，从而对长时间关井后的井底静压进行预测，以便进行后续的计算。

205、基于该多个关井时间的井底静压，确定该气井的井底静压随关井时间变化的函数关系。

具体地，该步骤可以是基于每个该关井时间以及对应的井底静压，以为横坐标，以/>为纵坐标，绘制关系曲线。其中，Δt表示该关井时间，t_p表示气井关井前累积生产时间，P_w表示关井时间为Δt时的井底静压。

具体地，基于步骤202中获取的数据，可以绘制井底静压与关井时间之间的关系曲线，例如，对于上述大深XX井，获取到的关系曲线如图4所示，图4是本公开实施例提供的一种与/>的关系曲线图，在图4中，横坐标为/>纵坐标为/>由图4可知，该曲线的斜率的绝对值随/>逐渐变化，但该曲线中包括直线段。

206、截取该关系曲线中的直线段。

该直线段能够更准确地反映关井足够长的时间后，该井底静压随时间的变化规律。

207、对该直线段进行拟合，得到与/>之间的函数关系。

将该通过实际数据获取的函数关系代入步骤203中获取的公式(3)中，即可获得各项所需的储层参数。

例如，对于上述大深XX井，该函数关系为y＝-331.43x+2808.2。

208、基于该对应关系以及该气井的井底静压随关井时间的函数关系，获取该气藏的储层参数，该储层参数包括：渗透率以及地层压力。

对于公式(3)，当以为自变量，以/>为因变量的函数的图像为一条直线时，则该直线的斜率m的表达式为下述公式(5)：

式中，m表示斜率，无因次；k表示渗透率，μm²；q_sc表示气井稳定生产产量，10⁴m³/d；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；μ表示在储层条件下的气体粘度，mPa·s；m为该直线段的斜率的绝对值，无因次。

具体地，对上述公式(5)进行变形，可以得到下述公式(6)，基于下述公式(6)，获取该地层系数：

式中，kh表示地层系数；q_sc表示气井稳定生产产量，10⁴m³/d；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；μ表示在储层条件下的气体粘度，mPa·s；m为该直线段的斜率的绝对值，无因次。

例如，对于上述步骤206中的函数关系，可得其斜率的绝对值m为331.43，将该m代入上述公式(6)，可以得到该大深XX井的地层系数kh＝6.11mD·m。

进一步地，对上述公式(6)进行变形，可以得到下述公式(2)，基于下述公式(2)，获取该渗透率：

式中，k表示渗透率，μm²；q_sc表示气井稳定生产产量，10⁴m³/d；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；μ表示在储层条件下的气体粘度，mPa·s；m为该直线段的斜率的绝对值，无因次。

例如，对于上述大深XX井，其储层有效厚度h＝9012m，将其代入上述公式(2)，可得该大深XX井的渗透率k＝0.66mD。

在一种可能实现方式中，基于下述公式(3)，获取该地层压力：

式中，P_r表示关井时间为Δt时的地层压力，MPa；P_w表示关井时间为Δt时的井底静压，MPa；q_sc表示气井稳定生产产量，10⁴m³/d；μ表示在储层条件下的气体粘度，mPa·s；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；k表示气藏渗透率，μm²；h表示气层有效厚度，m；t_p表示气井关井前累积生产时间，h；Δt表示关井时间，h。

具体地，当关井时间无限长时，Δt趋近于无穷大，井底压力应恢复到地层真实压力，该地层真实压力能够反映储层的储量等信息，可以为后续的生产提供可靠的依据。

基于此，在上述公式(3)中，当该等于1时，将所得的P_r的值获取为该地层压力。

例如，对于该大深XX井，其地层压力为52.9345。

上述步骤201～步骤208，说明了如何根据关井后短期内获取到的参数来获取储层参数的过程。该方法已经得到实际检验，例如，在获取大深XX井的储层参数的过程的同时，对同一关井时间的压力恢复试井解释成果进行对比验证。两者得到的储层参数的对比结果如下表2和图5所示：

表2

解释参数	试井解释结果	本实施例结果
			地层系数kh(mD·m)	5.42	6.11
渗透率(mD)	0.59	0.66
			地层压力(MPa)	52.9223	52.9345

该试井解释的结果如图5所示，图5是本公开实施例提供的一种试井解释图，在图5中，横坐标表示等效时间，纵坐标表示累积时间，图中上方的曲线表示压力恢复曲线，下方的曲线表示压力时间双对数曲线，通过拟合曲线，可以得到储层渗流相关参数和压力，由图可见，渗透率为表2中所示的0.59。

由表2和图5可知，本实施例提供的方法获取到的储层参数与经过长时间试井得到的试井解释结果比较接近，表明上述方法具有较强的可靠性和操作性。

上述所有可选技术方案，可以采用任意结合形成本公开的可选实施例，在此不再一一赘述。

本公开实施例提供的方法包括：获取气藏的地层压力、渗透率、井底静压以及关井时间之间的对应关系；基于关井后的较短时间内的多个关井时间对应的井底静压，获取的该气井的井底静压随关井时间变化的函数关系；基于上述对应关系和函数关系，计算得到关井时间较长时的储层参数，避免了为了获取储层参数而进行长时间的关井，提高了气藏的经济效益。

图6是本公开实施例提供的一种利用短期关井数据预测低渗气藏的储层参数的装置的结构示意图，该装置包括：

对应关系获取模块601，用于基于平面径向流渗流数学模型以及气藏特性，采用假想注气井法，获取该气藏的地层压力、渗透率、井底静压以及关井时间之间的对应关系；

井口油压获取模块602，用于在该气藏中的气井关井后，获取多个关井时间对应的井口油压，每个该关井时间小于预设时长；

井底静压获取模块603，用于基于该关井时间及相应的井口油压，获取该多个关井时间的井底静压；

函数关系获取模块604，用于基于该多个关井时间的井底静压，确定该气井的井底静压随关井时间变化的函数关系；

储层参数获取模块605，用于基于该对应关系以及该气井的井底静压随关井时间的函数关系，获取该气藏的储层参数，该储层参数包括：渗透率以及地层压力。

在一种可能实现方式中，该井底静压获取模块603，用于基于该气井的气层中部垂深、井筒静气柱平均温度、稳定生产产量以及气井关井前累积生产时间，该气井生产的气体在储层条件下的相对密度、井筒静气柱平均偏差因子以及气体粘度，在该气井关井后，每隔预设时长获取到的关井时间及相应的井口油压，以及下述公式(1)，获取每个该关井时间对应的井底静压：

/>

式中，P_w表示关井时间为Δt时的井底静压，MPa；P_wh表示关井时间为Δt时的井口油压，MPa；γ_g表示气体的相对密度，无因次；H表示气井的气层中部垂深，m；表示井筒静气柱平均温度，K；/>表示井筒静气柱平均偏差因子，无因次；

在一种可能实现方式中，该函数关系获取模块604，用于：

基于每个该关井时间以及对应的井底静压，以为横坐标，以/>为纵坐标，绘制关系曲线，其中，Δt表示该关井时间，t_p表示气井关井前累积生产时间，P_w表示关井时间为Δt时的井底静压；

截取该关系曲线中的直线段；

对该直线段进行拟合，得到与/>之间的函数关系。

在一种可能实现方式中，基于下述公式(2)，获取该渗透率：

式中，k表示渗透率，μm²；q_sc表示气井稳定生产产量，10⁴m³/d；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；μ表示在储层条件下的气体粘度，mPa·s；m为该直线段的斜率的绝对值，无因次。

在一种可能实现方式中，基于下述公式(3)，获取该地层压力：

式中，P_r表示关井时间为Δt时的地层压力，MPa；P_w表示关井时间为Δt时的井底静压，MPa；q_sc表示气井稳定生产产量，10⁴m³/d；μ表示在储层条件下的气体粘度，mPa·s；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；k表示气藏渗透率，μm²；h表示气层有效厚度，m；t_p表示气井关井前累积生产时间，h；Δt表示关井时间，h；

当该等于1时，将所得的P_r的值获取为该地层压力。

需要说明的是：上述实施例提供的利用短期关井数据预测低渗气藏的储层参数的装置在利用短期关井数据预测低渗气藏的储层参数时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将装置的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的利用短期关井数据预测低渗气藏的储层参数的装置与利用短期关井数据预测低渗气藏的储层参数的方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。

本公开实施例提供的装置，通过获取气藏的地层压力、渗透率、井底静压以及关井时间之间的对应关系；基于关井后的较短时间内的多个关井时间对应的井底静压，获取的该气井的井底静压随关井时间变化的函数关系；基于上述对应关系和函数关系，计算得到关井时间较长时的储层参数，避免了为了获取储层参数而进行长时间的关井，提高了气藏的经济效益。

图7是本公开实施例提供的一种计算机设备的结构示意图，该计算机设备700可因配置或性能不同而产生比较大的差异，可以包括一个或一个以上处理器(centralprocessing units，CPU)701和一个或一个以上的存储器702，其中，上述存储器702中存储有至少一条指令，上述至少一条指令由上述处理器701加载并执行以实现上述各个方法实施例提供的利用短期关井数据预测低渗气藏的储层参数的方法。当然，该计算机设备还可以具有有线或无线网络接口、键盘以及输入输出接口等部件，以便进行输入输出，该计算机设备还可以包括其他用于实现设备功能的部件，在此不做赘述。

在示例性实施例中，还提供了一种计算机可读存储介质，例如包括指令的存储器，上述指令可由计算机设备中的处理器执行以完成上述实施例中的利用短期关井数据预测低渗气藏的储层参数方法。例如，该计算机可读存储介质可以是只读存储器(Read-OnlyMemory,ROM)、随机存取存储器(Random Access Memory，RAM)、只读光盘(Compact DiscRead-Only Memory，CD-ROM)、磁带、软盘和光数据存储设备等。

本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，上述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。

上述仅为本公开的可选实施例，并不用以限制本公开，凡在本公开的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本公开的保护范围之内。

Claims (8)

1.一种利用短期关井数据预测低渗气藏的储层参数的方法，其特征在于，所述方法包括：

基于平面径向流渗流数学模型以及气藏特性，采用假想注气井法，获取所述气藏的地层压力、渗透率、井底静压以及关井时间之间的对应关系；

在所述气藏中的气井关井后，获取多个关井时间对应的井口油压，每个所述关井时间小于预设时长；

基于所述关井时间及相应的井口油压，获取所述多个关井时间的井底静压；

基于每个所述关井时间以及对应的井底静压，以为横坐标，以/>为纵坐标，绘制关系曲线，其中，Δt表示所述关井时间，t_p表示气井关井前累积生产时间，P_w表示关井时间为Δt时的井底静压；

截取所述关系曲线中的直线段；

对所述直线段进行拟合，得到与/>之间的函数关系；

基于所述对应关系以及所述与/>之间的函数关系，获取所述气藏的储层参数，所述储层参数包括：渗透率以及地层压力。

2.根据权利要求1所述的方法，其特征在于，所述基于所述关井时间及相应的井口油压，获取所述多个关井时间的井底静压，包括：

基于所述气井的气层中部垂深、井筒静气柱平均温度、稳定生产产量以及气井关井前累积生产时间，所述气井生产的气体在储层条件下的相对密度、井筒静气柱平均偏差因子以及气体粘度，在所述气井关井后，每隔预设时长获取到的关井时间及相应的井口油压，以及下述公式(1)，获取每个所述关井时间对应的井底静压：

式中，P_w表示关井时间为Δt时的井底静压，MPa；P_wh表示关井时间为Δt时的井口油压，MPa；γ_g表示气体的相对密度，无因次；H表示气井的气层中部垂深，m；表示井筒静气柱平均温度，K；/>表示井筒静气柱平均偏差因子，无因次。

3.根据权利要求1所述的方法，其特征在于，基于下述公式(2)，获取所述渗透率：

式中，k表示渗透率，μm²；q_sc表示气井稳定生产产量，10⁴m³/d；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；μ表示在储层条件下的气体粘度，mPa·s；m为所述直线段的斜率的绝对值，无因次。

4.根据权利要求1所述的方法，其特征在于，基于下述公式(3)，获取所述地层压力：

式中，P_r表示关井时间为Δt时的地层压力，MPa；P_w表示关井时间为Δt时的井底静压，MPa；q_sc表示气井稳定生产产量，10⁴m³/d；μ表示在储层条件下的气体粘度，mPa·s；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；k表示气藏渗透率，μm²；h表示气层有效厚度，m；t_p表示气井关井前累积生产时间，h；Δt表示关井时间，h；

当等于1时，将所得的P_r的值获取为所述地层压力。

5.一种利用短期关井数据预测低渗气藏的储层参数的装置，其特征在于，所述装置包括：

对应关系获取模块，用于基于平面径向流渗流数学模型以及气藏特性，采用假想注气井法，获取所述气藏的地层压力、渗透率、井底静压以及关井时间之间的对应关系；

井口油压获取模块，用于在所述气藏中的气井关井后，获取多个关井时间对应的井口油压，每个所述关井时间小于预设时长；

井底静压获取模块，用于基于所述关井时间及相应的井口油压，获取所述多个关井时间的井底静压；

函数关系获取模块，用于基于每个所述关井时间以及对应的井底静压，以为横坐标，以/>为纵坐标，绘制关系曲线，其中，Δt表示所述关井时间，t_p表示气井关井前累积生产时间，P_w表示关井时间为Δt时的井底静压；截取所述关系曲线中的直线段；对所述直线段进行拟合，得到/>与/>之间的函数关系；

储层参数获取模块，用于基于所述对应关系以及所述与/>之间的函数关系，获取所述气藏的储层参数，所述储层参数包括：渗透率以及地层压力。

6.根据权利要求5所述的装置，其特征在于，所述井底静压获取模块，用于基于所述气井的气层中部垂深、井筒静气柱平均温度、稳定生产产量以及气井关井前累积生产时间，所述气井生产的气体在储层条件下的相对密度、井筒静气柱平均偏差因子以及气体粘度，在所述气井关井后，每隔预设时长获取到的关井时间及相应的井口油压，以及下述公式(1)，获取每个所述关井时间对应的井底静压：

式中，P_w表示关井时间为Δt时的井底静压，MPa；P_wh表示关井时间为Δt时的井口油压，MPa；γ_g表示气体的相对密度，无因次；H表示气井的气层中部垂深，m；表示井筒静气柱平均温度，K；/>表示井筒静气柱平均偏差因子，无因次。

7.根据权利要求5所述的装置，其特征在于，基于下述公式(2)，获取所述渗透率：

式中，k表示渗透率，μm²；q_sc表示气井稳定生产产量，10⁴m³/d；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；μ表示在储层条件下的气体粘度，mPa·s；m为所述直线段的斜率的绝对值，无因次。

8.根据权利要求5所述的装置，其特征在于，基于下述公式(3)，获取所述地层压力：

式中，P_r表示关井时间为Δt时的地层压力，MPa；P_w表示关井时间为Δt时的井底静压，MPa；q_sc表示气井稳定生产产量，10⁴m³/d；μ表示在储层条件下的气体粘度，mPa·s；Z表示气层中部压力和温度下的气体偏差因子，无因次；T表示气层中部温度，K；k表示气藏渗透率，μm²；h表示气层有效厚度，m；t_p表示气井关井前累积生产时间，h；Δt表示关井时间，h；

当等于1时，将所得的P_r的值获取为所述地层压力。