[发明专利]一种储气库有效储气空间计算方法

权利要求书

1.一种储气库有效储气空间计算方法，其特征在于，按照以下步骤进行：

S1、建立储气库三维地质模型；

根据油气藏的地质静态资料和油气水纵向分布规律，建立储气库三维地质模型；

S2、网格系统粗化；

在储气库三维地质模型的基础上，进行网格及属性数据体粗化；

S3、准备数值模拟基础参数；

准备所需的数值模拟基础参数，该数值模拟基础参数包括储气库储层参数与流体参数；

S4、获取相对渗透率曲线；

根据从实验样品测试得到油气相对渗透率数据和油水相对渗透率数据，通过归一化处理，得到油气相对渗透率曲线和油水相对渗透率曲线；

S5、获取PVT参数；

利用生产井的井流物高压物性实验数据，建立地层流体状态方程，拟合得到PVT参数，并确定临界温度和临界压力；

S6、建立储气库数值模拟模型；

基于储气库三维地质模型、数值模拟基础参数、相对渗透率曲线和PVT参数，建立储气库数值模拟模型，并对储气库数值模拟模型进行初始化设置；

S7、对储气库数值模拟模型进行全过程历史拟合；

根据储气库生产历史动态特征及流体分布特征，对储气库数值模拟模型进行全过程历史拟合；

S8、基于临界饱和度参数对储气库进行区带划分；

按照流体饱和度参数的不同，将储气库的各渗透区带分别划分为纯气区、气油过渡带、纯油区和水淹区；

S9、提取储气库的孔隙体积；

将区带划分结果以数据流的形式导入储气库数值模拟模型中，提取储气库的各渗透区带的孔隙体积，该孔隙体积包括含气孔隙体积、含油孔隙体积、含水孔隙体积和含烃孔隙体积；

S10、确定储气库各渗透区带的动用效率；

根据多周期注采模拟实验结果，确定储气库各渗透区带的动用效率；

S11、确定储气库的有效储气空间；

根据储气库的各渗透区带的孔隙体积以及动用效率，计算得到储气库的有效储气空间。

2.根据权利要求1所述的一种储气库有效储气空间计算方法，其特征在于，所述步骤S1按照以下步骤进行：

S1-1、设定模型的网格单元和尺度；

S1-2、建立储气库构造模型；

S1-3、建立储气库沉积微相模型；

S1-4、建立储气库岩相模型；

S1-5、建立储气库属性模型，该储气库属性模型包括孔隙度模型、渗透率模型和饱和度模型；

S1-6、检验储气库构造模型、储气库沉积微相模型、储气库岩相模型和储气库属性模型的精度是否合格：否，返回步骤S1-1，并调整参数；是，进入步骤S2。

3.根据权利要求2所述的一种储气库有效储气空间计算方法，其特征在于：所述步骤S1-2中，基于储气库区域的单井基础资料、地震解释资料和区域地质资料建立储气库的构造模型，其中，储气库由多个叠合在一起的小层组成，将每个小层的层面模型叠合在一起组成储气库构造模型。

4.根据权利要求3所述的一种储气库有效储气空间计算方法，其特征在于：所述步骤S1-3中，对各个小层的二维沉积微相分布图进行边界数字化，采用赋值方法建立储气库沉积微相模型。

5.根据权利要求4所述的一种储气库有效储气空间计算方法，其特征在于：所述步骤S1-4中，在完成各井点单砂体精细识别的基础上，在等时的小层内部，以各单井识别的单砂体顶底数据为约束，确定各单砂体顶部和底部的等时砂体包络面，以此包络面为约束，建立单砂体岩相模型以及砂体间的隔夹层泥相模型，从而由各单砂体岩相模型和各隔夹层泥相模型组成储气库岩相模型。

6.根据权利要求5所述的一种储气库有效储气空间计算方法，其特征在于：所述步骤S1-5中，在储气库沉积微相模型和储气库岩相模型的基础上，以沉积相为约束，首先分析各个沉积相内的储层物性相关性，再利用随机模拟方法对测井数据进行插值计算的同时，通过协克里金方法，对模拟计算进行加权和条件约束，使测井数值的插值与井间分布相贴近，可实现确定性建模和随机性建模相结合，分别建立了孔隙度模型、渗透率模型和饱和度模型。

7.根据权利要求1所述的一种储气库有效储气空间计算方法，其特征在于：所述步骤S4中，归一化处理采用以下公式进行：

油水两相系统：

式(1)和式(2)中，K_ro表示油的相对渗透率；表示束缚水饱和度下的油相相对渗透率，mD；S_w表示含水饱和度；S_orw表示油水两相中残余油饱和度，％；S_wc表示束缚水饱和度；K_rw表示水的相对渗透率；表示残余油饱和度下水相相对渗透率；n_o表示油相相对渗透率曲线指数，无因次；n_w表示水相对渗透率曲线指数，无因次；

油气两相系统：

式(3)和式(4)中，表示残余气饱和度下的油相相对渗透率，mD；S_g表示含气饱和度；S_lc表示总临界液体饱和度，％，并且，S_lc＝S_wc+S_org；S_org表示油气两相中残余油饱和度，％；S_gc表示临界气饱和度，％；K_rg表示气的相对渗透率；表示残余油饱和度下气相相对渗透率；S_gc表示残余气饱和度，％；n_g表示气相对渗透率曲线指数，无因次；n_go表示气油两相中油的相对渗透率曲线指数，无因次。