获取由水淹气藏改建成的储气库的库容量的方法和装置制造方法
【专利摘要】本发明提供了一种获取由水淹气藏改建成的储气库的库容量的方法和装置，该方法包括：获取储气库的气水界面随储层压力变化的变化关系；根据变化关系将储气库的气藏动态储量占据的孔隙体积划分为气区、气水交互区和水区三个储气库流体区域，并计算各个储气库流体区域对气藏动态储量贡献的孔隙体积；根据各个储气库流体区域对气藏动态储量贡献的孔隙体积对应扣除各个储气库流体区域内不能用于气体储存的孔隙体积；在扣除各个储气库流体区域内不能用于气体储存的孔隙体积基础上，计算储气库的库容量。本发明的计算结果准确度较高，而且随着气藏水体能力的增强，其计算精度会进一步提高，从而可为储气库前期评价提供准确的库容量技术参数。
【专利说明】获取由水淹气藏改建成的储气库的库容量的方法和装置
【技术领域】
[0001]本发明涉天然气地下储存及【技术领域】，尤其是涉及一种获取由水淹气藏改建成的储气库的库容量的方法和装置。
【背景技术】
[0002]储气库是调节天然气管网供求平衡的有力手段，衡量储气库调节能力的重要指标是库容量，准确计算库容量是困扰相关科研人员的技术难题。
[0003]目前，水淹气藏改建储气库后的库容量一般取动态储量数值，(如文献2012年12月第31卷第12期《油气储运》的文章“低渗岩性气藏改建地下储气库工作气量的确定”，和2013年5月第20卷第3期《断块油气田》的文章“不封闭期藏建库库容设计”介绍的计算方法)。该类计算方法认为，动态储量所占据的孔隙体积改建储气库后将完全用于储气库气体储存。但在现有的水淹气藏改建储气库的工程中，没有考虑水侵对储层的影响，错误地将一部分动态储量占据的孔隙体积计算为建库后用于气体储存的储气库孔隙体积，导致设计库容量偏大，而且这种偏差将随着水侵量的增加而扩大。
【发明内容】
[0004]本发明的目的在于提供一种获取由水淹气藏改建成的储气库的库容量的方法和装置，以获得更为准确的由水淹气藏改建成的储气库的库容量。
[0005]为达到上述目的，一方面，本发明提供了一种获取由水淹气藏改建成的储气库的库容量的方法，包括以下步骤:
[0006]获取储气库的气水界面随储层压力变化的变化关系；
[0007]根据所述变化关系将所述储气库的气藏动态储量占据的孔隙体积划分为气区、气水交互区和水区三个储气库流体区域，并计算各个储气库流体区域对所述气藏动态储量贡献的孔隙体积；
[0008]根据各个储气库流体区域对所述气藏动态储量贡献的孔隙体积对应扣除各个储气库流体区域内不能用于气体储存的孔隙体积；
[0009]在扣除所述各个储气库流体区域内不能用于气体储存的孔隙体积基础上，计算所述储气库的库容量。
[0010]本发明的获取由水淹气藏改建成的储气库的库容量的方法，所述获取储气库的气水界面随储层压力变化的变化关系，包括:
[0011]根据公式W = f(P)获取气水界面的净水侵量与储层压力的关系；
[0012]根据公式V = GX Bg获取储气库的动态储量对应的孔隙体积V ；
[0013]其中，W为压力P对应的净水侵量，f (P)为净水侵量和储层压力满足的函数关系，G为气藏动态储量，BgS天然气在原始压力时的体积系数。
[0014]本发明的获取由水淹气藏改建成的储气库的库容量的方法，所述根据所述变化关系将所述储气库的气藏动态储量占据的孔隙体积划分为气区、气水交互区和水区三个储气库流体区域，并计算各个储气库流体区域对所述气藏动态储量贡献的孔隙体积，包括:
[0015]利用储气库下限压力对应的净水侵量Wmin和改建储气库时的原始储层压力对应的净水侵量wab，计算气区净水侵量W1= W ab-wmin;
[0016]根据所述气区净水侵量^计算所述气区对所述气藏动态储量贡献的孔隙体积V !
= va-swcl-sgrl);
[0017]利用储气库下限压力对应的净水侵量Wmin和储气库上限压力对应的净水侵量Wmax，计算气水交互区净水侵量W2= Wmax-Wmin;
[0018]根据所述气水交互区净水侵量W2计算所述气水交互区对所述气藏动态储量贡献的孔隙体积 V2= W 2/(1-Swc2-Sgr2);
[0019]水区净水侵量W3等于储气库上限压力对应的净水侵量Wmax;
[0020]根据水区净水侵量W3计算所述水区对所述气藏动态储量贡献的孔隙体积V3= W3/
(I Swc3 Sgr3)；
[0021]其中，Swc;1、Sirc2和S wc;3分别为气区、气水交互区和水区的束缚术饱和度，S grl, Sgrf和Sgri分别为气区、气水交互区和水区因水侵而产生的残余气体饱和度。
[0022]本发明的获取由水淹气藏改建成的储气库的库容量的方法，所述根据各个储气库流体区域对所述气藏动态储量贡献的孔隙体积对应扣除各个储气库流体区域内不能用于气体储存的孔隙体积，包括:
[0023]根据公式AV1= (Swcln-SwJSgrtn)V1获得气区内不能用于气体储存的孔隙体积AV1;
[0024]根据公式AV2= (Swc;2n_Swc;2+Sgun)V^得气水交互区内不能用于气体储存的孔隙体积AV2;
[0025]根据公式AV3= (1-SwJV3获得水区内不能用于气体储存的孔隙体积AV3;
[0026]其中，Swcln, Swc;2dP Sire3l^别为储气库稳定运行时气区、气水交互区和水区的束缚术饱和度。
[0027]本发明的获取由水淹气藏改建成的储气库的库容量的方法，所述Swc;1、Sirc2和S wc3,Sgrn Sgrt和S gr3,以及Swc;ln、Swc2n^P S κ3η通过以下方式得到:
[0028]根据各个储气库流体区域的特征分别设计对应的储气库渗流实验，并通过所述储气库渗流实验测得各个储气库流体区域的束缚术饱和度swc;1、Sire2和Sire3，残余气体饱和度Sgrp Sgr2和S gr3,以及储气库稳定运行时的束缚术饱和度swc;ln、Swc2n^P S wc;3n。
[0029]本发明的获取由水淹气藏改建成的储气库的库容量的方法，所述在扣除所述各个储气库流体区域内不能用于气体储存的孔隙体积基础上，计算所述储气库的库容量，包括:
[0030]根据公式Vv= V- Λ V厂Λ V2- Δ %获得储气库用于气体储存的孔隙体积V ν;
[0031]根据公式Gg= V V/Bg获得的库容量G g;
[0032]其中，BgS天然气在原始压力时的体积系数。
[0033]另一方面，本发明还提供了一种获取由水淹气藏改建成的储气库的库容量的装置，包括:
[0034]变化关系获取模块，用于获取储气库的气水界面随储层压力变化的变化关系；
[0035]孔隙体积划分模块，用于根据所述变化关系将所述储气库的气藏动态储量占据的孔隙体积划分为气区、气水交互区和水区三个储气库流体区域，并计算各个储气库流体区域对所述气藏动态储量贡献的孔隙体积；
[0036]体积扣除模块，用于根据各个储气库流体区域对所述气藏动态储量贡献的孔隙体积对应扣除各个储气库流体区域内不能用于气体储存的孔隙体积；
[0037]库容量计算模块，用于在扣除所述各个储气库流体区域内不能用于气体储存的孔隙体积基础上，计算所述储气库的库容量。
[0038]本发明的获取由水淹气藏改建成的储气库的库容量的装置，所述获取储气库的气水界面随储层压力变化的变化关系，包括:
[0039]根据公式W = f (P)获取气水界面的净水侵量与储层压力的关系；
[0040]根据公式V = GX Bg获取储气库的动态储量对应的孔隙体积V ；
[0041]其中，W为压力P对应的净水侵量，f (P)为净水侵量和储层压力满足的函数关系，G为气藏动态储量，BgS天然气在原始压力时的体积系数。
[0042]本发明的获取由水淹气藏改建成的储气库的库容量的装置，所述根据所述变化关系将所述储气库的气藏动态储量占据的孔隙体积划分为气区、气水交互区和水区三个储气库流体区域，并计算各个储气库流体区域对所述气藏动态储量贡献的孔隙体积，包括:
[0043]利用储气库下限压力对应的净水侵量Wmin和改建储气库时的原始储层压力对应的净水侵量wab，计算气区净水侵量W1= W ab-wmin;
[0044]根据所述气区净水侵量WJf算所述气区对所述气藏动态储量贡献的孔隙体积V !
= va-swcl-sgrl);
[0045]利用储气库下限压力对应的净水侵量Wmin和储气库上限压力对应的净水侵量Wmax，计算气水交互区净水侵量W2= Wmax-Wmin;
[0046]根据所述气水交互区净水侵量W2计算所述气水交互区对所述气藏动态储量贡献的孔隙体积 V2= W 2/(1-Swc2-Sgr2);
[0047]水区净水侵量胃3等于储气库上限压力对应的净水侵量W max;
[0048]根据水区净水侵量W3计算所述水区对所述气藏动态储量贡献的孔隙体积V3= W3/
(I Swc3 Sgr3)；
[0049]其中，Swc;1、Sirc2和S wc;3分别为气区、气水交互区和水区的束缚术饱和度，S grl, Sgrf和Sgri分别为气区、气水交互区和水区因水侵而产生的残余气体饱和度。
[0050]本发明的获取由水淹气藏改建成的储气库的库容量的装置，所述根据各个储气库流体区域对所述气藏动态储量贡献的孔隙体积对应扣除各个储气库流体区域内不能用于气体储存的孔隙体积，包括:
[0051]根据公式AV1= ^-、。--^、获得气区内不能用于气体储存的孔隙体积AV1;
[0052]根据公式AV2= (Swe2n-Swc;2+S#2n)V^得气水交互区内不能用于气体储存的孔隙体积AV2;
[0053]根据公式AV3= (1-SwJV3获得水区内不能用于气体储存的孔隙体积AV3;
[0054]其中，Swcln, Swc;2dP Sire3l^别为储气库稳定运行时气区、气水交互区和水区的束缚术饱和度。
[0055]本发明的获取由水淹气藏改建成的储气库的库容量的装置，所述Swc;1、Sirc2和S wc3,Sgrn Sgrt和S gr3,以及Swc;ln、Swc2n^P S K3n通过以下方式得到:
[0056]根据各个储气库流体区域的特征分别设计对应的储气库渗流实验，并通过所述储气库渗流实验测得各个储气库流体区域的束缚术饱和度swc;1、Sire2和Sire3，残余气体饱和度Sgrp Sgr2和S gr3,以及储气库稳定运行时的束缚术饱和度swc;ln、Swc2n^P S wc;3n。
[0057]本发明的获取由水淹气藏改建成的储气库的库容量的装置，所述在扣除所述各个储气库流体区域内不能用于气体储存的孔隙体积基础上，计算所述储气库的库容量，包括:
[0058]根据公式Vv= V- Λ V厂Λ V2- Δ %获得储气库用于气体储存的孔隙体积V ν;
[0059]根据公式Gg= V V/Bg获得的库容量G g;
[0060]其中，BgS天然气在原始压力时的体积系数。
[0061]本发明实施例在计算库容量时，没有像现有技术那样直接取动态储量的数值，而是将动态储量占据的孔隙体积分为气区、气水交互区和水区这三个储气库流体区域，分别扣除各个储气库流体区域的动态储量改建储气库后不能用于储气库气体储存的孔隙体积，因此，本发明实施例的计算结果准确度大幅度提高，而且随着气藏水体能力的增强，本发明实施例的计算精度会进一步提高，技术优势凸显，从而可为储气库前期评价提供准确的库容量技术参数。
【专利附图】
【附图说明】
[0062]此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中:
[0063]图1为本发明实施例的获取由水淹气藏改建成的储气库的库容量的方法流程图；
[0064]图2为本发明实施例的获取由水淹气藏改建成的储气库的库容量的方法中储气库孔隙体积分割示意图；
[0065]图3为本发明实施例的获取由水淹气藏改建成的储气库的库容量的结构框图。
【具体实施方式】
[0066]为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
[0067]下面结合附图，对本发明的【具体实施方式】作进一步的详细说明。
[0068]参考图1所示，本发明实施例的获取由水淹气藏改建成的储气库的库容量的方法包括以下步骤:
[0069]步骤S101、获取储气库的气水界面随储层压力变化的变化关系。具体的，根据公式W = f(p)获取气水界面的净水侵量与储层压力的关系；根据公式V = GXBg获取储气库的动态储量对应的孔隙体积V。其中，其中，W为-储层平均压力P对应的净水侵量，可由气藏工程方法或数值模拟方法得到；f (P)为净水侵量和储层压力满足的函数关系(例如净水侵量和储层平均压力回归成函数关系，以便计算任一储层压力的净水侵量)；G为气藏动态储量由气藏工程方法或数值模拟方法得到，以便计算动态储量对应的孔隙体积；Bg为天然气在原始压力时的体积系数。
[0070]本发明实施例中，储气库的气水界面与常规气藏开发的平面形状不同，本发明实施例中气水界面可以是任意曲面，通过净水侵量来表示:
[0071]Wfflax= f (Pffl J
[0072]Wfflin= f (Pfflin)
[0073]Wab= f (Pab)
[0074]式中，Pmaxm储气库上限压力，MPa ；P minm储气库下限压力，MPa ；P &为改建储气库时原始气藏压力，MPa。
[0075]步骤S102、根据所述变化关系将所述储气库的气藏动态储量占据的孔隙体积划分为气区、气水交互区和水区三个储气库流体区域(如图2所示)，并计算各个储气库流体区域对气藏动态储量贡献的孔隙体积。具体的:
[0076]气区位于储气库下限压力对应的气水界面以上的储层，储气库运行阶段一直处于气区状态；利用储气库下限压力对应的净水侵量Wmin和改建储气库时的原始储层压力对应的净水侵量wab，计算气区净水侵量W1 = W ab-wmin;
[0077]根据气区净水侵量W1计算所述气区对气藏动态储量贡献的孔隙体积V1= W1/(I Swci Sgri);
[0078]气水交互区位于储气库下限压力与下限压力对应的气水界面之间的储层，储气库注满气时为气区、采气末期为水区；利用储气库下限压力对应的净水侵量Wmin和储气库上限压力对应的净水侵量wmax，计算气水交互区净水侵量W2= W max-wmin；
[0079]根据所述气水交互区净水侵量W2计算气水交互区对所述气藏动态储量贡献的孔隙体积 V2= W 2/ (1-Swc2-Sgr2)；
[0080]水区位于储气库上限压力对应的气水界面以下的储层，储气库运行阶段一直处于水淹状态；水区净水侵量W3等于储气库上限压力对应的净水侵量Wmax;
[0081]根据水区净水侵量W3计算所述水区对所述气藏动态储量贡献的孔隙体积V3= W3/
(I Swc3 Sgr3)；
[0082]其中，Swel、Swe2和S we3分别为气区、气水交互区和水区的束缚术饱和度，S grl, Sgrf和Sgri分别为气区、气水交互区和水区因水侵而产生的残余气体饱和度。
[0083]步骤S103、根据各个储气库流体区域对气藏动态储量贡献的孔隙体积对应扣除各个储气库流体区域内不能用于气体储存的孔隙体积。具体的:
[0084]根据公式AV1= (Swcln-SwJSgrtn)V1获得气区内不能用于气体储存的孔隙体积AV1;
[0085]根据公式AV2= (3￥。211-5￥。2+5#211)￥2获得气水交互区内不能用于气体储存的孔隙体积AV2;
[0086]根据公式AV3= (1-SwJV3获得水区内不能用于气体储存的孔隙体积AV3;
[0087]其中，Swcln, Swc;2, S ^分别为储气库稳定运行时气区、气水交互区和水区的束缚术饱和度。
[0088]步骤S104、在扣除各个储气库流体区域内不能用于气体储存的孔隙体积基础上，计算储气库的库容量。具体的:
[0089]根据公式Vv= V- Λ V厂Λ V2- Δ %获得储气库用于气体储存的孔隙体积V ν;
[0090]根据公式Gg= V V/Bg获得的库容量G g;
[0091]其中，BgS天然气在原始压力时的体积系数。
[0092]本发明实施例中，上述Swc;1、Sire2和 S wc3, Sgrt、Sgi2和 S gr3,以及 Swc;ln、Swc2n^P S wc;3n可通过以下方式得到:
[0093]根据各个储气库流体区域的特征分别设计对应的储气库渗流实验，并通过储气库渗流实验测得各个储气库流体区域的束缚术饱和度swc;1、Sirc2和Sire3，残余气体饱和度sgrt、Sgrt和S gr3,以及储气库稳定运行时的束缚术饱和度Swc;ln、Swc2n^P S wc3no具体的:
[0094]①气驱水实验，描述气区注气驱水的物理过程。得到气区的端点值饱和度:得到气区束缚术饱和度Swca，气区因为水侵而产生的残余气饱和度Swcl，储气库稳定运行气区束缚术饱和度Swcln，储气库稳定运行气区因为水侵而产生的残余气饱和度swc;ln。
[0095]②气驱水-水驱气实验，描述气水交互区气水互驱替的物理过程，得到气水交互区的端点值饱和度:气水交互区束缚术饱和度Swe2，气水交互区因为水侵而产生的残余气饱和度Swe2，储气库稳定运行气水交互区束缚术饱和度Swe2n，储气库稳定运行气水交互区因为水侵而产生的残余气饱和度Swc;2n。
[0096]③水驱气实验，描述水区水驱气的物理过程，得到水区的端点值饱和度:水区束缚术饱和度Sirc3，水区因为水侵而产生的残余气饱和度Sire3，储气库稳定运行水区束缚术饱和度Sire3n，储气库稳定运行水区因为水侵而产生的残余气饱和度swc;3n。
[0097]与本发明实施例上述的方法对应，如图3所示，本发明实施例的获取由水淹气藏改建成的储气库的库容量的装置包括变化关系获取模块31、孔隙体积划分模块32、体积扣除模块33和库容量计算模块34。其中:
[0098]变化关系获取模块31，用于获取储气库的气水界面随储层压力变化的变化关系。具体的，根据公式W = f(P)获取气水界面的净水侵量与储层压力的关系；根据公式V =GXBg获取储气库的动态储量对应的孔隙体积V。其中，其中，W为-储层平均压力P对应的净水侵量，可由气藏工程方法或数值模拟方法得到；f(P)为净水侵量和储层压力满足的函数关系(例如净水侵量和储层平均压力回归成函数关系，以便计算任一储层压力的净水侵量)；G为气藏动态储量由气藏工程方法或数值模拟方法得到，以便计算动态储量对应的孔隙体积；Bg为天然气在原始压力时的体积系数。
[0099]本发明实施例中，储气库的气水界面与常规气藏开发的平面形状不同，本发明实施例中气水界面可以是任意曲面，通过净水侵量来表示:
[0100]Wfflax= f (Pffl J
[0101]Wfflin= f (Pfflin)
[0102]Wab= f (Pab)
[0103]式中，Pmaxm储气库上限压力，MPa ;P minm储气库下限压力，MPa ；P &为改建储气库时原始气藏压力，MPa。
[0104]孔隙体积划分模块32，用于根据所述变化关系将所述储气库的气藏动态储量占据的孔隙体积划分为气区、气水交互区和水区三个储气库流体区域，并计算各个储气库流体区域对所述气藏动态储量贡献的孔隙体积。具体的:
[0105]气区位于储气库下限压力对应的气水界面以上的储层，储气库运行阶段一直处于气区状态；利用储气库下限压力对应的净水侵量Wmin和改建储气库时的原始储层压力对应的净水侵量wab，计算气区净水侵量W1 = W ab-wmin;
[0106]根据气区净水侵量W1计算所述气区对气藏动态储量贡献的孔隙体积V1= W1/(I Swci Sgri);
[0107]气水交互区位于储气库下限压力与下限压力对应的气水界面之间的储层，储气库注满气时为气区、采气末期为水区；利用储气库下限压力对应的净水侵量Wmin和储气库上限压力对应的净水侵量wmax，计算气水交互区净水侵量W2= W max-wmin；
[0108]根据所述气水交互区净水侵量W2计算气水交互区对所述气藏动态储量贡献的孔隙体积 V2= W 2/ (1-Swc2-Sgr2)；
[0109]水区位于储气库上限压力对应的气水界面以下的储层，储气库运行阶段一直处于水淹状态；水区净水侵量W3等于储气库上限压力对应的净水侵量Wmax;
[0110]根据水区净水侵量W3计算所述水区对所述气藏动态储量贡献的孔隙体积V3= W3/
(I Swc3 Sgr3)；
[0111]其中，Swel、Swe2和Swe3分别为气区、气水交互区和水区的束缚术饱和度，S grt、sgrf和Sgri分别为气区、气水交互区和水区因水侵而产生的残余气体饱和度。
[0112]体积扣除模块33，用于根据各个储气库流体区域对所述气藏动态储量贡献的孔隙体积对应扣除各个储气库流体区域内不能用于气体储存的孔隙体积。具体的:
[0113]根据公式AV1= (Sweln-Swc^sgrtn)V1获得气区内不能用于气体储存的孔隙体积AV1;
[0114]根据公式AV2= 61^-51?。2+5^1)￥2获得气水交互区内不能用于气体储存的孔隙体积AV2;
[0115]根据公式AV3= (1-SwJV3获得水区内不能用于气体储存的孔隙体积AV3;
[0116]其中，Swcln, Swe2l^ Swe3l^别为储气库稳定运行时气区、气水交互区和水区的束缚术饱和度。
[0117]库容量计算模块34，用于在扣除所述各个储气库流体区域内不能用于气体储存的孔隙体积基础上，计算所述储气库的库容量。具体的:
[0118]根据公式Vv= V- Λ V厂Λ V2- Δ %获得储气库用于气体储存的孔隙体积V ν;
[0119]根据公式Gg= V V/Bg获得的库容量Gg;
[0120]其中，BgS天然气在原始压力时的体积系数。
[0121]本发明实施例中，上述Swc;1、Swc^P S wc;3，Sgrt、S@和 S gr3,以及 Swc;ln、swc;2l^p Sire3n可通过上述方法实施例的方式得到，在次不再赘述。
[0122]由此可见，本发明实施例在计算库容量时，没有像现有技术那样直接取动态储量的数值，而是将动态储量占据的孔隙体积分为气区、气水交互区和水区这三个储气库流体区域，分别扣除各个储气库流体区域的动态储量改建储气库后不能用于储气库气体储存的孔隙体积，因此，本发明实施例的计算结果准确度大幅度提高，而且随着气藏水体能力的增强，本发明实施例的计算精度会进一步提高，技术优势凸显，从而可为储气库前期评价提供准确的库容量技术参数。
[0123]本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块、单元和步骤可以通过硬件、软件或两者的结合来实现。至于是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用，可以使用各种方法实现所述的功能，但这种实现不应被理解为超出本发明实施例保护的范围。
[0124]本发明实施例中所描述的各种说明性的逻辑块，或单元都可以通过通用处理器，数字信号处理器，专用集成电路0310，现场可编程门阵列或其它可编程逻辑装置，离散门或晶体管逻辑，离散硬件部件，或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器，可选地，该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现，例如数字信号处理器和微处理器，多个微处理器，一个或多个微处理器联合一个数字信号处理器核，或任何其它类似的配置来实现。
[0125]本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于狀1存储器、闪存、1？01存储器、存储器、存储器、寄存器、硬盘、可移动磁盘、⑶-801或本领域中其它任意形式的存储媒介中。示例性地，存储媒介可以与处理器连接，以使得处理器可以从存储媒介中读取信息，并可以向存储媒介存写信息。可选地，存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于…X中，八31(:可以设置于用户终端中。可选地，处理器和存储媒介也可以设置于用户终端中的不同的部件中。
[0126]在一个或多个示例性的设计中，本发明实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的媒介上，或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但不限于狀1、801、22？1？01、0)-801或其它光盘存储、磁盘存储或其它磁性存储装置，或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外，任何连接都可以被适当地定义为电脑可读媒介，例如，如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(01)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片((114)和磁盘((11%)包括压缩磁盘、镭射盘、光盘、1^0、软盘和蓝光光盘，磁盘通常以磁性复制数据，而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。
[0127]以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
【权利要求】
1.一种获取由水淹气藏改建成的储气库的库容量的方法，其特征在于，包括以下步骤: 获取储气库的气水界面随储层压力变化的变化关系； 根据所述变化关系将所述储气库的气藏动态储量占据的孔隙体积划分为气区、气水交互区和水区三个储气库流体区域，并计算各个储气库流体区域对所述气藏动态储量贡献的孔隙体积； 根据各个储气库流体区域对所述气藏动态储量贡献的孔隙体积对应扣除各个储气库流体区域内不能用于气体储存的孔隙体积； 在扣除所述各个储气库流体区域内不能用于气体储存的孔隙体积基础上，计算所述储气库的库容量。
2.根据权利要求1所述的获取由水淹气藏改建成的储气库的库容量的方法，其特征在于，所述获取储气库的气水界面随储层压力变化的变化关系，包括: 根据公式W = f(P)获取气水界面的净水侵量与储层压力的关系； 根据公式V = GXBg获取储气库的动态储量对应的孔隙体积V ； 其中，W为压力P对应的净水侵量，f (P)为净水侵量和储层压力满足的函数关系，G为气藏动态储量，BgS天然气在原始压力时的体积系数。
3.根据权利要求2所述的获取由水淹气藏改建成的储气库的库容量的方法，其特征在于，所述根据所述变化关系将所述储气库的气藏动态储量占据的孔隙体积划分为气区、气水交互区和水区三个储气库流体区域，并计算各个储气库流体区域对所述气藏动态储量贡献的孔隙体积，包括: 利用储气库下限压力对应的净水侵量Wmin和改建储气库时的原始储层压力对应的净水侵量wab，计算气区净水侵量W1= W ab-wmin; 根据所述气区净水侵量W1计算所述气区对所述气藏动态储量贡献的孔隙体积V i = W i/(I Swci Sgri); 利用储气库下限压力对应的净水侵量Wmin和储气库上限压力对应的净水侵量Wmax，计算气水交互区净水侵量W2= W max-wmin； 根据所述气水交互区净水侵量胃2计算所述气水交互区对所述气藏动态储量贡献的孔隙体积 V2= W 2/ (1-Swc2-Sgr2)； 水区净水侵量胃3等于储气库上限压力对应的净水侵量W max; 根据水区净水侵量W3计算所述水区对所述气藏动态储量贡献的孔隙体积V 3= W 3/(I Swc3 Sgr3)； 其中，Swcl, Sire2和Sire3分别为气区、气水交互区和水区的束缚术饱和度，S grl, Sgrt和S gr3分别为气区、气水交互区和水区因水侵而产生的残余气体饱和度。
4.根据权利要求3所述的获取由水淹气藏改建成的储气库的库容量的方法，其特征在于，所述根据各个储气库流体区域对所述气藏动态储量贡献的孔隙体积对应扣除各个储气库流体区域内不能用于气体储存的孔隙体积，包括: 根据公式AV1= (Swcln-Swc^sgrtn)V1获得气区内不能用于气体储存的孔隙体积AV1; 根据公式AV2= (Swc;2n-Swc;2+Sgm)V2获得气水交互区内不能用于气体储存的孔隙体积AV2; 根据公式AV3= (1-SwJV3获得水区内不能用于气体储存的孔隙体积AV3; 其中，Swc;ln、Swc;2dP Sire3n*别为储气库稳定运行时气区、气水交互区和水区的束缚术饱和度。
5.根据权利要求4所述的获取由水淹气藏改建成的储气库的库容量的方法，其特征在于，所述 swc;1、Sire2和 S wc3, Sgrl, Sgrf和 S gr3,以及 Swc;ln、Swc2n^P S wc;3n通过以下方式得到: 根据各个储气库流体区域的特征分别设计对应的储气库渗流实验，并通过所述储气库渗流实验测得各个储气库流体区域的束缚术饱和度swc;1、Sirc2和Sire3，残余气体饱和度sgrt、Sgrt和S gr3,以及储气库稳定运行时的束缚术饱和度Swc;ln、Swc2n^P S wc;3n。
6.根据权利要求5所述的获取由水淹气藏改建成的储气库的库容量的方法，其特征在于，所述在扣除所述各个储气库流体区域内不能用于气体储存的孔隙体积基础上，计算所述储气库的库容量，包括: 根据公式Vv= V-AV1-A V2- Δ %获得储气库用于气体储存的孔隙体积V V； 根据公式Gg= V V/Bg获得的库容量G g; 其中，BgS天然气在原始压力时的体积系数。
7.一种获取由水淹气藏改建成的储气库的库容量的装置，其特征在于，包括: 变化关系获取模块，用于获取储气库的气水界面随储层压力变化的变化关系； 孔隙体积划分模块，用于根据所述变化关系将所述储气库的气藏动态储量占据的孔隙体积划分为气区、气水交互区和水区三个储气库流体区域，并计算各个储气库流体区域对所述气藏动态储量贡献的孔隙体积； 体积扣除模块，用于根据各个储气库流体区域对所述气藏动态储量贡献的孔隙体积对应扣除各个储气库流体区域内不能用于气体储存的孔隙体积； 库容量计算模块，用于在扣除所述各个储气库流体区域内不能用于气体储存的孔隙体积基础上，计算所述储气库的库容量。
8.根据权利要求7所述的获取由水淹气藏改建成的储气库的库容量的装置，其特征在于，所述获取储气库的气水界面随储层压力变化的变化关系，包括: 根据公式W = f(P)获取气水界面的净水侵量与储层压力的关系； 根据公式V = GXBg获取储气库的动态储量对应的孔隙体积V ； 其中，W为压力P对应的净水侵量，f (P)为净水侵量和储层压力满足的函数关系，G为气藏动态储量，BgS天然气在原始压力时的体积系数。
9.根据权利要求8所述的获取由水淹气藏改建成的储气库的库容量的装置，其特征在于，所述根据所述变化关系将所述储气库的气藏动态储量占据的孔隙体积划分为气区、气水交互区和水区三个储气库流体区域，并计算各个储气库流体区域对所述气藏动态储量贡献的孔隙体积，包括: 利用储气库下限压力对应的净水侵量Wmin和改建储气库时的原始储层压力对应的净水侵量wab，计算气区净水侵量W1= W ab-wmin; 根据所述气区净水侵量W1计算所述气区对所述气藏动态储量贡献的孔隙体积V i = W i/(I Swci Sgri); 利用储气库下限压力对应的净水侵量Wmin和储气库上限压力对应的净水侵量Wmax，计算气水交互区净水侵量W2= W max-wmin； 根据所述气水交互区净水侵量胃2计算所述气水交互区对所述气藏动态储量贡献的孔隙体积 V2= W 2/ (1-Swc2-Sgr2)； 水区净水侵量胃3等于储气库上限压力对应的净水侵量W max; 根据水区净水侵量W3计算所述水区对所述气藏动态储量贡献的孔隙体积V 3= W 3/(I Swc3 Sgr3)； 其中，Swcl, Sire2和Sire3分别为气区、气水交互区和水区的束缚术饱和度，S grl, Sgrt和S gr3分别为气区、气水交互区和水区因水侵而产生的残余气体饱和度。
10.根据权利要求9所述的获取由水淹气藏改建成的储气库的库容量的装置，其特征在于，所述根据各个储气库流体区域对所述气藏动态储量贡献的孔隙体积对应扣除各个储气库流体区域内不能用于气体储存的孔隙体积，包括: 根据公式AV1= (Swcln-Swc^sgrtn)V1获得气区内不能用于气体储存的孔隙体积AV1; 根据公式AV2= (Swc;2n_Swc;2+Sgm)V2获得气水交互区内不能用于气体储存的孔隙体积AV2; 根据公式AV3= (1-SwJV3获得水区内不能用于气体储存的孔隙体积AV3; 其中，Swc;ln、Swc;2dP Sire3n*别为储气库稳定运行时气区、气水交互区和水区的束缚术饱和度。
11.根据权利要求10所述的获取由水淹气藏改建成的储气库的库容量的装置，其特征在于，所述 swc;1、Sire2和 S wc3, Sgrl, Sgi2和 S gr3,以及 Swc;ln、Swc2n^P S wc;3n通过以下方式得到: 根据各个储气库流体区域的特征分别设计对应的储气库渗流实验，并通过所述储气库渗流实验测得各个储气库流体区域的束缚术饱和度swc;1、Sirc2和Sire3，残余气体饱和度sgrt、Sgrt和S gr3,以及储气库稳定运行时的束缚术饱和度Swc;ln、Swc2n^P S wc;3n。
12.根据权利要求11所述的获取由水淹气藏改建成的储气库的库容量的装置，其特征在于，所述在扣除所述各个储气库流体区域内不能用于气体储存的孔隙体积基础上，计算所述储气库的库容量，包括: 根据公式Vv= V-AV1-A V2- Δ %获得储气库用于气体储存的孔隙体积V V； 根据公式Gg= V V/Bg获得的库容量G g; 其中，BgS天然气在原始压力时的体积系数。
【文档编号】G06F19/00GK104484550SQ201410636368
【公开日】2015年4月1日 申请日期:2014年11月6日 优先权日:2014年11月6日
【发明者】王皆明, 唐立根, 胥洪成, 李春 申请人:中国石油天然气股份有限公司