CN113404485A - 一种凝析气顶油藏的油气界面移动速度的确定方法及装置 - Google Patents
一种凝析气顶油藏的油气界面移动速度的确定方法及装置 Download PDFInfo
- Publication number
- CN113404485A CN113404485A CN202010186156.9A CN202010186156A CN113404485A CN 113404485 A CN113404485 A CN 113404485A CN 202010186156 A CN202010186156 A CN 202010186156A CN 113404485 A CN113404485 A CN 113404485A
- Authority
- CN
- China
- Prior art keywords
- gas
- oil
- gas cap
- volume
- condensate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 99
- 239000011148 porous material Substances 0.000 claims abstract description 97
- 230000008859 change Effects 0.000 claims abstract description 73
- 230000009545 invasion Effects 0.000 claims abstract description 41
- 239000011435 rock Substances 0.000 claims abstract description 40
- 239000000463 material Substances 0.000 claims abstract description 35
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 239000007789 gas Substances 0.000 claims description 477
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 37
- 238000004519 manufacturing process Methods 0.000 claims description 27
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 26
- 238000011161 development Methods 0.000 claims description 22
- 238000004590 computer program Methods 0.000 claims description 17
- 239000003345 natural gas Substances 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 10
- 230000033001 locomotion Effects 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 239000003921 oil Substances 0.000 description 207
- 238000010586 diagram Methods 0.000 description 18
- 238000000605 extraction Methods 0.000 description 9
- 239000010779 crude oil Substances 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 238000004891 communication Methods 0.000 description 6
- 239000012530 fluid Substances 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 239000011800 void material Substances 0.000 description 3
- 230000005465 channeling Effects 0.000 description 2
- 238000009833 condensation Methods 0.000 description 2
- 230000005494 condensation Effects 0.000 description 2
- 238000007872 degassing Methods 0.000 description 2
- 238000001704 evaporation Methods 0.000 description 2
- 230000008020 evaporation Effects 0.000 description 2
- 239000008398 formation water Substances 0.000 description 2
- 239000008239 natural water Substances 0.000 description 2
- 239000003129 oil well Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 206010021137 Hypovolaemia Diseases 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000004090 dissolution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/10—Locating fluid leaks, intrusions or movements
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
Abstract
本发明提供一种凝析气顶油藏的油气界面移动速度的确定方法及装置，所述方法包括：获取气顶区孔隙变化体积和油环区孔隙变化体积；根据所述气顶区孔隙变化体积、所述油环区孔隙变化体积、油气藏内岩石的总膨胀体积以及凝析气顶油藏物质平衡方程，获得凝析气顶油藏的当前地层压力；其中，所述凝析气顶油藏物质平衡方程是预设的；根据所述当前地层压力和气顶侵入体积，获得所述凝析气顶油藏的油气界面移动速度。所述装置用于执行上述方法。本发明实施例提供的凝析气顶油藏的油气界面移动速度的确定方法及装置，提高了凝析气顶油藏衰竭开采方式下的油气界面移动速度的计算效率。
Description
技术领域
本发明涉及石油开采技术领域，具体涉及一种凝析气顶油藏的油气界面移动速度的确定方法及装置。
背景技术
在自然界中，具有自由气和原油共存的油气藏中，大部分是凝析气顶油藏。在深层高压条件下，凝析气顶油藏的原油多数具有油质轻的特点。在凝析气顶油藏中，一般存在油气水三相或油气两相可动流体共存，油气水按重力差异，它们之间有油气和油水两个界面存在。
凝析气顶油藏投入开发时，由于油和气都处于饱和状态，一旦地层压力下降，就会在凝析气顶中发生反凝析现象，造成凝析油损失；在油环中出现原油脱气现象，增加原油粘度，加大油区开发难度。在凝析气顶油藏开发过程中处理好气顶与油环之间的相互作用、相互影响，防止油气互窜，保持油气界面稳定或缓慢移动，是实现合理开发该类油气藏的关键因素。因此，凝析气顶油藏的开发难度比单纯的气藏或油藏要高得多。当凝析气顶压力高于油环压力时，气顶中的凝析气容易向下侵入油环，从而导致靠近油气界面的油井发生气窜，并降低油井的产油能力；而当油环压力高于气顶压力时，油环中的原油则会向上侵入到气顶中，由于侵入的原油处于分散状态而不易被采出到地面，因此造成了部分原油资源的损失。因此，如何准确识别凝析气顶油藏油气界面的移动情况，对于该类油气藏的开发调整具有重要意义。
发明内容
针对现有技术中的问题，本发明实施例提供一种凝析气顶油藏的油气界面移动速度的确定方法及装置，能够至少部分地解决现有技术中存在的问题。
一方面，本发明提出一种凝析气顶油藏的油气界面移动速度的确定方法，包括：
获取气顶区孔隙变化体积和油环区孔隙变化体积；
根据所述气顶区孔隙变化体积、所述油环区孔隙变化体积、油气藏内岩石的总膨胀体积以及凝析气顶油藏物质平衡方程，获得凝析气顶油藏的当前地层压力；其中，所述凝析气顶油藏物质平衡方程是预设的；
根据所述当前地层压力和气顶侵入体积，获得所述凝析气顶油藏的油气界面移动速度。
另一方面，本发明提供一种凝析气顶油藏的油气界面移动速度的确定装置，包括：
获取单元，用于获取气顶区孔隙变化体积和油环区孔隙变化体积；
第一获得单元，用于根据所述气顶区孔隙变化体积、所述油环区孔隙变化体积、油气藏内岩石的总膨胀体积以及凝析气顶油藏物质平衡方程，获得凝析气顶油藏的当前地层压力；其中，所述凝析气顶油藏物质平衡方程是预设的；
第二获得单元，用于根据所述当前地层压力和气顶侵入体积，获得所述凝析气顶油藏的油气界面移动速度。
再一方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述任一实施例所述凝析气顶油藏的油气界面移动速度的确定方法的步骤。
又一方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述任一实施例所述凝析气顶油藏的油气界面移动速度的确定方法的步骤。
本发明实施例提供的凝析气顶油藏的油气界面移动速度的确定方法及装置，能够获取气顶区孔隙变化体积和油环区孔隙变化体积，根据气顶区孔隙变化体积、油环区孔隙变化体积、油气藏内岩石的总膨胀体积以及凝析气顶油藏物质平衡方程，获得凝析气顶油藏的当前地层压力，根据当前地层压力和气顶侵入体积，获得凝析气顶油藏的油气界面移动速度，提高了凝析气顶油藏衰竭开采方式下的油气界面移动速度的计算效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1a是本发明一实施例提供的凝析气顶油藏的开发前的流体分布示意图；
图1b是本发明一实施例提供的凝析气顶油藏的开发后的流体分布示意图；
图2是本发明一实施例提供的凝析气顶油藏的油气界面移动速度的确定方法的流程示意图。
图3是本发明另一实施例提供的凝析气顶油藏的油气界面移动速度的确定方法的流程示意图。
图4是本发明一实施例提供的衰竭开采方式下采油速度、采气速度与油气界面移动速度关系示意图。
图5是本发明一实施例提供的凝析气顶油藏的油气界面移动速度的确定装置的结构示意图。
图6是本发明另一实施例提供的凝析气顶油藏的油气界面移动速度的确定装置的结构示意图。
图7是本发明一实施例提供的电子设备的实体结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。
为了便于理解本申请提供的技术方案，下面先对本申请技术方案的相关内容进行说明。在本发明实施例中，假设凝析气顶油藏中凝析气顶与底部油环处于同一压力***内，且油气界面均匀推进；原始地层压力高于凝析气体的露点压力；凝析气顶与油环均有束缚水存在；忽略凝析气顶气在油环中的溶解以及油环逸出溶解气窜入凝析气顶；忽略气体在地层的吸附现象。本发明实施例提供的凝析气顶油藏的油气界面移动速度的确定方法的执行主体包括但不限于计算机。
图1a是本发明一实施例提供的凝析气顶油藏的开发前的流体分布示意图，如图1a所示，凝析气顶油藏在开发前存在油、气、水三相，层次分明。图1b是本发明一实施例提供的凝析气顶油藏的开发后的流体分布示意图，如图1b所示，经过一段时间的开发，当凝析气顶压力高于油环压力时，凝析气顶中发生反凝析现象，形成凝析油；当油环压力高于气顶压力时，在油环中出现原油脱气现象，逸出溶解气。
图2是本发明一实施例提供的凝析气顶油藏的油气界面移动速度的确定方法的流程示意图，如图2所示，本发明实施例提供的凝析气顶油藏的油气界面移动速度的确定方法包括：
S201、获取气顶区孔隙变化体积和油环区孔隙变化体积；
具体地，由于在凝析气顶油藏的开发过程中，随着地层压力的不断下降，凝析气顶中的凝析油会不断析出，同时气顶中的原生水也会不断蒸发。因此，考虑凝析油析出和地层水蒸发等因素的影响，当地层压力降至当前地层压力时，气顶区孔隙变化体积可以根据如下公式计算获得：
其中，ΔVG为所述气顶区孔隙变化体积，单位可以采用m3，G为所述凝析气顶油藏的地层原始天然气储量，单位可以采用m3，Gp为气顶区域内累积天然气采出量，单位可以采用m3，Bg为当前气顶气体的体积系数，Bgi为原始条件下气顶气的体积系数，yw为当前地层压力下凝析气相中水蒸汽百分数，Sco为气顶体积内凝析油的饱和度，SwcG为气顶体积内的束缚水饱和度，ywi为原始条件下气相中水蒸汽含量，以小数表示。
在凝析气顶油藏的开发过程中，随着地层压力的持续降低，油环中的溶解气会不断逸出，继而成为油藏的游离气。其中一部分游离气会流入井筒而被采出至地面，而另外一部分游离气会滞留在地层中。因此，考虑油环溶解气逸出的影响，地层压力降至当前地层压力时的油环孔隙体积的变化为：
其中，ΔVO为油环区孔隙变化体积，单位为m3，m为原始条件下油环孔隙体积与气顶孔隙体积之比，G为所述凝析气顶油藏的地层原始天然气储量，Bgi为原始条件下气顶气的体积系数，SwcO为油环体积内的束缚水饱和度，Bo为油环油在当前地层压力下的体积系数，ywi为原始条件下气相中水蒸汽含量，SwcG为气顶体积内的束缚水饱和度，Boi为油环油在原始压力下的体积系数，Np为采出油环油的地面体积，ΔSw为油环区含水饱和度增量，SgO为当前油环区内的含气饱和度。其中，所述原始条件是指凝析气顶油藏未投入开发时的地层条件(压力及温度等)。
当前油环区内的含气饱和度SgO，可以通过下式计算获得：
其中，SwcO为油环体积内的束缚水饱和度，Rsi为原始条件下的溶解气油比，Rs为当前条件下的溶解气油比，Bg为当前气顶气体的体积系数，Boi为油环油在原始压力下的体积系数，Np为采出油环油的地面体积，Rp为当前条件下的生产气油比ywi为原始条件下气相中水蒸汽含量，SwcG为气顶体积内的束缚水饱和度，m为原始条件下油环孔隙体积与气顶孔隙体积之比，G为所述凝析气顶油藏的地层原始天然气储量，Bgi为原始条件下气顶气的体积系数。
考虑油藏天然水侵的影响，油环区含水饱和度增量ΔSw可通过下式求得：
其中，We为天然水侵量，单位可以采用m3，Wp为当前累积产水量，Bw为水的体积系数，ywi为原始条件下气相中水蒸汽含量，SwcG为气顶体积内的束缚水饱和度，m为原始条件下油环孔隙体积与气顶孔隙体积之比，G为所述凝析气顶油藏的地层原始天然气储量，Bgi为原始条件下气顶气的体积系数。
S202、根据所述气顶区孔隙变化体积、所述油环区孔隙变化体积、油气藏内岩石的总膨胀体积以及凝析气顶油藏物质平衡方程，获得凝析气顶油藏的当前地层压力；其中，所述凝析气顶油藏物质平衡方程是预设的；
具体地，在获得所述气顶区孔隙变化体积和所述油环区孔隙变化体积之后，可以根据所述气顶区孔隙变化体积、所述油环区孔隙变化体积、油气藏内岩石的总膨胀体积以及凝析气顶油藏物质平衡方程，获得凝析气顶油藏的当前地层压力，即将所述气顶区孔隙变化体积和所述油环区孔隙变化体积带入到凝析气顶油藏物质平衡方程中，并所述油气藏内岩石的总膨胀体积通过所述凝析气顶油藏的当前地层压力表示，获得关于所述凝析气顶油藏的当前地层压力的方程，求解该方程即可获得凝析气顶油藏的当前地层压力。其中，所述凝析气顶油藏物质平衡方程是预设的。
S203、根据所述当前地层压力和气顶侵入体积，获得所述凝析气顶油藏的油气界面移动速度。
具体地，在获得所述当前地层压力之后，根据所述当前地层压力和气顶侵入体积，可以获得所述凝析气顶油藏的油气界面移动速度。
例如，根据所述当前地层压力以及第一预设公式，计算获得所述气顶侵入体积，然后根据所述气顶侵入体积以及第二预设公式，计算获得所述油气界面移动速度。其中，所述第一预设公式和所述第二预设公式是预先设定的。
本发明实施例提供的凝析气顶油藏的油气界面移动速度的确定方法，能够获取气顶区孔隙变化体积和油环区孔隙变化体积，根据气顶区孔隙变化体积、油环区孔隙变化体积、油气藏内岩石的总膨胀体积以及凝析气顶油藏物质平衡方程，获得凝析气顶油藏的当前地层压力，根据当前地层压力和气顶侵入体积，获得凝析气顶油藏的油气界面移动速度，提高了凝析气顶油藏衰竭开采方式下的油气界面移动速度的计算效率。
在上述各实施例的基础上，进一步地，所述凝析气顶油藏物质平衡方程为：
ΔVG+ΔVO+ΔVTf＝0
其中，ΔVG表示所述气顶区孔隙变化体积，ΔVO表示所述油环区孔隙变化体积，ΔVTf表示所述油气藏内岩石的膨胀体积，所述油气藏内岩石的总膨胀体积根据如下公式计算获得：
其中，
为所述凝析气顶油藏的平均岩石压缩系数，pi为所述凝析气顶油藏的原始地层压力，m为原始条件下油环孔隙体积与气顶孔隙体积之比，p为所述凝析气顶油藏的当前地层压力，G为所述凝析气顶油藏的地层原始天然气储量，Bgi为原始条件下气顶气的体积系数，ywi为原始条件下气相中水蒸汽含量，SwcG为气顶体积内的束缚水饱和度。
具体地，分别考虑气顶区和油环区内岩石的膨胀作用，当所述凝析气顶油藏的原始地层压力降至当前地层压力时，油气藏内岩石的总膨胀体积为：
其中，ΔVTf为油气藏内岩石的总膨胀体积，单位可以采用m3，ΔVGf为气顶内岩石的总膨胀体积，ΔVOf为油环内岩石的总膨胀体积量，
为所述凝析气顶油藏的平均岩石压缩系数。
根据凝析气顶油藏开发过程中的体积守恒原理，即原始油气藏孔隙体积等于当前油气藏孔隙体积与油气藏内岩石的总膨胀体积之和，可以获得所述凝析气顶油藏物质平衡方程为：
ΔVG+ΔVO+ΔVTf＝0 (6)
将公式(5)带入到公式(6)中，可以获得关于凝析气顶油藏的当前地层压力的方程，通过迭代计算的方法便可以求解获得所述凝析气顶油藏的当前地层压力p。具体的过程为，假定当前油藏压力下降Δp，即当前地层压力p＝pi-Δp，通过公式(5)计算获得此时的油气藏内岩石的膨胀体积，并计算此时的气顶区孔隙变化体积和油环区孔隙变化体积，然后代入所述凝析气顶油藏物质平衡方程(6)中，对比方程(6)的左边和右边是否相等，如果相等，则当前地层压力p＝pi-Δp，如果不相等，则需要重新假设一个地层压力的变化量Δp，并重复上述计算过程，直到获得满足计算精度要求的当前地层压力p为止。其中，计算精度要求根据实际需要进行设置，本发明实施例不做限定。
图3是本发明另一实施例提供的凝析气顶油藏的油气界面移动速度的确定方法的流程示意图，如图3所示，在上述各实施例的基础上，进一步地，所述根据所述当前地层压力和气顶侵入体积，获得所述凝析气顶油藏的油气界面移动速度包括：
S2031、根据所述当前地层压力以及第一预设公式，计算获得所述气顶侵入体积；
具体地，在获得所述当前地层压力之后，可以将当前地层压力代入到第一预设公式中，计算获得所述气顶侵入体积。
S2032、根据所述气顶侵入体积以及第二预设公式，计算获得所述油气界面移动速度。
具体地，在计算获得所述气顶侵入体积之后，将所述气顶侵入体积代入到第二预设公式中，计算获得所述油气界面移动速度。
在上述各实施例的基础上，进一步地，所述第一预设公式为：
其中，VGinvade为所述气顶侵入体积，G为所述凝析气顶油藏的地层原始天然气储量，Gp为气顶区域内累积天然气采出量，Bg为当前气顶气体的体积系数，yw为当前地层压力下凝析气相中水蒸汽百分数，Sco为气顶体积内凝析油的饱和度，SwcG为气顶体积内的束缚水饱和度，Bgi为原始条件下气顶气的体积系数，
为所述凝析气顶油藏的平均岩石压缩系数，pi为所述凝析气顶油藏的原始地层压力，p为所述凝析气顶油藏的当前地层压力。
具体地，在凝析气顶与油环协同开发条件下，假设油气界面向油环方向移动，可以得到气顶侵入体积VGinvade的计算公式为：
其中，VG为气顶区当前条件下的孔隙体积，ΔVGf为气顶内岩石的总膨胀体积，VGi为气顶区原始条件下的孔隙体积。
在上述各实施例的基础上，进一步地，所述第二预设公式为：
其中，vog为所述油气界面移动速度，VGinvade为所述气顶侵入体积，t为开发时间，L为油气界面内外边界间的宽度，W为油藏宽度，φ为地层孔隙度，SwcO为油环体积内的束缚水饱和度，Sor为残余油饱和度，α为地层倾角。
具体地，根据气顶侵入体积，利用容积法可以计算出油气界面的移动速度，假设内、外油气界面移动速度相等，即油气界面平行下移，气顶侵入体积VGinvade可以表示为：
VGinvade＝LhWφ(1-SwcO-Sor)＝LWφ(1-SwcO-Sor)xsinα (8)
其中，L为油气界面内外边界间的宽度，单位可以采用m，h为油气界面纵向移动的距离，单位可以采用m，W为油藏宽度，x为油气界面横向移动的距离，单位可以采用m，φ为地层孔隙度，SwcO为油环体积内的束缚水饱和度，Sor为残余油饱和度，α为地层倾角，单位可以采用度。
油气界面横向移动速度vog可以表示为：
vog＝x/t (9)
其中，vog为油气界面横向移动速度，单位可以采用米/年，x为油气界面横向移动的距离，单位可以采用m，t为开发时间，单位可以采用年。此处的油气界面横向移动速度即上述各实施例所述的油气界面移动速度。
根据公式(8)和公式(9)，可以获得所述第二预设公式为：
利用衰竭开采方式开展凝析气顶与油环协同开采，油气界面的移动方向受气顶与油环的体积亏空差来决定。当气顶亏空体积大于油环时，油气界面向气顶方向移动；而当油环亏空体积大于气顶时，油气界面便向油环方向移动。因此，衰竭式油气协同开采时，影响油气界面稳定的主要开发因素为气顶采气速度与油环采油速度。
以让纳若尔油气田的Г北气顶油藏为研究对象，分别计算出不同油环采油速度、不同气顶采气速度下的油气界面移动速度，建立采油速度、采气速度与油气界面移动速度的关系，如图4所示，其中采气速度设定为0％、1％、2％、4％、6％，采油速度设定为0％、0.3％、0.6％、0.9％、1.2％。在图4中，以横坐标轴为界，横轴以上代表油环区域，横轴以下代表气顶区域。从图4还可以看出，不同的采油速度、采气速度组合，油气界面的移动方向是不同的，由此也印证了我们对衰竭式油气同采时影响油气界面稳定的主控因素的推断。
从图4所示的衰竭开采方式下采油速度、采气速度与油气界面移动速度关系可以看出：
(1)当气顶亏空大于油环亏空时，油气界面偏向气区方向移动。而且采气速度相同时，采油速度越大，油环的亏空速度加快导致气顶油环间的压力差变小，油气界面移动速度越小；当采油速度相同时，采气速度越大导致气顶压力下降越快，气顶油环间的压力差越来越大，油气界面的移动速度变得越来越大。
(2)当油环亏空大于气顶亏空时，油气界面偏向油区方向移动。而且相同采气速度下，采油速度越大，油环压力下降越快，气顶油环间的压力差越来越大，油气界面向油区移动速度越大；相同采油速度下，采气速度越大，气顶压力下降越快，气顶油环间的压力差逐渐降低，油气界面移动速度也变得越来越小。
(3)对于某一油环采油速度，均存在一个合理的采气速度，实现油气界面的相对稳定(即油气界面移动速度为零)；而且油气界面稳定时，合理的采气速度、采油速度成正比例关系。由此可以看出，衰竭式气顶油环协同开采时，只要合理地匹配采油速度和采气速度的大小关系，保持油气界面相对稳定，就可以有效防止油侵或气侵现象的发生。
图5是本发明一实施例提供的凝析气顶油藏的油气界面移动速度的确定装置的结构示意图，如图5所示，本发明实施例提供的凝析气顶油藏的油气界面移动速度的确定装置包括获取单元501、第一获得单元502和第二获得单元503，其中：
获取单元501用于获取气顶区孔隙变化体积和油环区孔隙变化体积；第一获得单元502用于根据所述气顶区孔隙变化体积、所述油环区孔隙变化体积、油气藏内岩石的总膨胀体积以及凝析气顶油藏物质平衡方程，获得凝析气顶油藏的当前地层压力；其中，所述凝析气顶油藏物质平衡方程是预设的；第二获得单元503用于根据所述当前地层压力和气顶侵入体积，获得所述凝析气顶油藏的油气界面移动速度。
具体地，由于在凝析气顶油藏的开发过程中，随着地层压力的不断下降，凝析气顶中的凝析油会不断析出，同时气顶中的原生水也会不断蒸发。因此，考虑凝析油析出和地层水蒸发等因素的影响，当地层压力降至当前地层压力时，获取单元501可以根据公式(1)计算获得气顶区孔隙变化体积。
在凝析气顶油藏的开发过程中，随着地层压力的持续降低，油环中的溶解气会不断逸出，继而成为油藏的游离气。其中一部分游离气会流入井筒而被采出至地面，而另外一部分游离气会滞留在地层中。因此，考虑油环溶解气逸出的影响，地层压力降至当前地层压力时，获取单元501可以根据公式(2)计算获得油环孔隙变化体积。
在获得所述气顶区孔隙变化体积和所述油环区孔隙变化体积之后，第一获得单元502可以根据所述气顶区孔隙变化体积、所述油环区孔隙变化体积、油气藏内岩石的总膨胀体积以及凝析气顶油藏物质平衡方程，获得凝析气顶油藏的当前地层压力，即将所述气顶区孔隙变化体积和所述油环区孔隙变化体积带入到凝析气顶油藏物质平衡方程中，并所述油气藏内岩石的总膨胀体积通过所述凝析气顶油藏的当前地层压力表示，获得关于所述凝析气顶油藏的当前地层压力的方程，求解该方程即可获得凝析气顶油藏的当前地层压力。其中，所述凝析气顶油藏物质平衡方程是预设的。
在获得所述当前地层压力之后，第二获得单元503根据所述当前地层压力和气顶侵入体积，可以获得所述凝析气顶油藏的油气界面移动速度。
本发明实施例提供的凝析气顶油藏的油气界面移动速度的确定装置，能够获取气顶区孔隙变化体积和油环区孔隙变化体积，根据气顶区孔隙变化体积、油环区孔隙变化体积、油气藏内岩石的总膨胀体积以及凝析气顶油藏物质平衡方程，获得凝析气顶油藏的当前地层压力，根据当前地层压力和气顶侵入体积，获得凝析气顶油藏的油气界面移动速度，提高了凝析气顶油藏衰竭开采方式下的油气界面移动速度的计算效率。
在上述各实施例的基础上，进一步地，所述凝析气顶油藏物质平衡方程为：
ΔVG+ΔVO+ΔVTf＝0
其中，ΔVG表示所述气顶区孔隙变化体积，ΔVO表示所述油环区孔隙变化体积，ΔVTf表示所述油气藏内岩石的总膨胀体积，所述油气藏内岩石的总膨胀体积根据如下公式计算获得：
其中，
为所述凝析气顶油藏的平均岩石压缩系数，pi为所述凝析气顶油藏的原始地层压力，m为原始条件下油环孔隙体积与气顶孔隙体积之比，p为所述凝析气顶油藏的当前地层压力，G为所述凝析气顶油藏的地层原始天然气储量，Bgi为原始条件下气顶气的体积系数，ywi为原始条件下气相中水蒸汽含量，SwcG为气顶体积内的束缚水饱和度。
图6是本发明另一实施例提供的凝析气顶油藏的油气界面移动速度的确定装置的结构示意图，如图6所示，在上述各实施例的基础上，进一步地，第二获得单元503包括第一计算子单元5031和第二计算子单元5032，其中：
第一计算子单元5031用于根据所述当前地层压力以及第一预设公式，计算获得所述气顶侵入体积；第二计算子单元5032用于根据所述气顶侵入体积以及第二预设公式，计算获得所述油气界面移动速度。
具体地，在获得所述当前地层压力之后，第一计算子单元5031可以将当前地层压力代入到第一预设公式中，计算获得所述气顶侵入体积。
在计算获得所述气顶侵入体积之后，第二计算子单元5032将所述气顶侵入体积代入到第二预设公式中，计算获得所述油气界面移动速度。
本发明实施例提供的装置的实施例具体可以用于执行上述各方法实施例的处理流程，其功能在此不再赘述，可以参照上述方法实施例的详细描述。
图7是本发明一实施例提供的电子设备的实体结构示意图，如图7所示，该电子设备可以包括：处理器(processor)701、通信接口(Communications Interface)702、存储器(memory)703和通信总线704，其中，处理器701，通信接口702，存储器703通过通信总线704完成相互间的通信。处理器701可以调用存储器703中的逻辑指令，以执行如下方法：获取气顶区孔隙变化体积和油环区孔隙变化体积；根据所述气顶区孔隙变化体积、所述油环区孔隙变化体积、油气藏内岩石的总膨胀体积以及凝析气顶油藏物质平衡方程，获得凝析气顶油藏的当前地层压力；其中，所述凝析气顶油藏物质平衡方程是预设的；根据所述当前地层压力和气顶侵入体积，获得所述凝析气顶油藏的油气界面移动速度。
此外，上述的存储器703中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
本实施例公开一种计算机程序产品，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，计算机能够执行上述各方法实施例所提供的方法，例如包括：获取气顶区孔隙变化体积和油环区孔隙变化体积；根据所述气顶区孔隙变化体积、所述油环区孔隙变化体积、油气藏内岩石的总膨胀体积以及凝析气顶油藏物质平衡方程，获得凝析气顶油藏的当前地层压力；其中，所述凝析气顶油藏物质平衡方程是预设的；根据所述当前地层压力和气顶侵入体积，获得所述凝析气顶油藏的油气界面移动速度。
本实施例提供一种计算机可读存储介质，所述计算机可读存储介质存储计算机程序，所述计算机程序使所述计算机执行上述各方法实施例所提供的方法，例如包括：获取气顶区孔隙变化体积和油环区孔隙变化体积；根据所述气顶区孔隙变化体积、所述油环区孔隙变化体积、油气藏内岩石的总膨胀体积以及凝析气顶油藏物质平衡方程，获得凝析气顶油藏的当前地层压力；其中，所述凝析气顶油藏物质平衡方程是预设的；根据所述当前地层压力和气顶侵入体积，获得所述凝析气顶油藏的油气界面移动速度。
本领域内的技术人员应明白，本发明的实施例可提供为方法、***、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在本说明书的描述中，参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种凝析气顶油藏的油气界面移动速度的确定方法，其特征在于，包括：
获取气顶区孔隙变化体积和油环区孔隙变化体积；
根据所述气顶区孔隙变化体积、所述油环区孔隙变化体积、油气藏内岩石的总膨胀体积以及凝析气顶油藏物质平衡方程，获得凝析气顶油藏的当前地层压力；其中，所述凝析气顶油藏物质平衡方程是预设的；
根据所述当前地层压力和气顶侵入体积，获得所述凝析气顶油藏的油气界面移动速度。
2.根据权利要求1所述的方法，其特征在于，所述凝析气顶油藏物质平衡方程为：
ΔVG+ΔVO+ΔVTf＝0
其中，ΔVG表示所述气顶区孔隙变化体积，ΔVO表示所述油环区孔隙变化体积，ΔVTf表示所述油气藏内岩石的总膨胀体积，所述油气藏内岩石的总膨胀体积根据如下公式计算获得：
其中，
为所述凝析气顶油藏的平均岩石压缩系数，pi为所述凝析气顶油藏的原始地层压力，m为原始条件下油环孔隙体积与气顶孔隙体积之比，p为所述凝析气顶油藏的当前地层压力，G为所述凝析气顶油藏的地层原始天然气储量，Bgi为原始条件下气顶气的体积系数，ywi为原始条件下气相中水蒸汽含量，SwcG为气顶体积内的束缚水饱和度。
3.根据权利要求1或2所述的方法，其特征在于，所述根据所述当前地层压力和气顶侵入体积，获得所述凝析气顶油藏的油气界面移动速度包括：
根据所述当前地层压力以及第一预设公式，计算获得所述气顶侵入体积；
根据所述气顶侵入体积以及第二预设公式，计算获得所述油气界面移动速度。
4.根据权利要求3所述的方法，其特征在于，所述第一预设公式为：
其中，VGinvade为所述气顶侵入体积，G为所述凝析气顶油藏的地层原始天然气储量，Gp为气顶区域内累积天然气采出量，Bg为当前气顶气体的体积系数，yw为当前地层压力下凝析气相中水蒸汽百分数，Sco为气顶体积内凝析油的饱和度，SwcG为气顶体积内的束缚水饱和度，Bgi为原始条件下气顶气的体积系数，
为所述凝析气顶油藏的平均岩石压缩系数，pi为所述凝析气顶油藏的原始地层压力，p为所述凝析气顶油藏的当前地层压力。
5.根据权利要求3所述的方法，其特征在于，所述第二预设公式为：
其中，vog为所述油气界面移动速度，VGinvade为所述气顶侵入体积，t为开发时间，L为油气界面内外边界间的宽度，W为油藏宽度，φ为地层孔隙度，SwcO为油环体积内的束缚水饱和度，Sor为残余油饱和度，α为地层倾角。
6.一种凝析气顶油藏的油气界面移动速度的确定装置，其特征在于，包括：
获取单元，用于获取气顶区孔隙变化体积和油环区孔隙变化体积；
第一获得单元，用于根据所述气顶区孔隙变化体积、所述油环区孔隙变化体积、油气藏内岩石的总膨胀体积以及凝析气顶油藏物质平衡方程，获得凝析气顶油藏的当前地层压力；其中，所述凝析气顶油藏物质平衡方程是预设的；
第二获得单元，用于根据所述当前地层压力和气顶侵入体积，获得所述凝析气顶油藏的油气界面移动速度。
7.根据权利要求6所述的装置，其特征在于，所述凝析气顶油藏物质平衡方程为：
ΔVG+ΔVO+ΔVTf＝0
其中，ΔVG表示所述气顶区孔隙变化体积，ΔVO表示所述油环区孔隙变化体积，ΔVTf表示所述油气藏内岩石的总膨胀体积，所述油气藏内岩石的总膨胀体积根据如下公式计算获得：
其中，
为所述凝析气顶油藏的平均岩石压缩系数，pi为所述凝析气顶油藏的原始地层压力，m为原始条件下油环孔隙体积与气顶孔隙体积之比，p为所述凝析气顶油藏的当前地层压力，G为所述凝析气顶油藏的地层原始天然气储量，Bgi为原始条件下气顶气的体积系数，ywi为原始条件下气相中水蒸汽含量，SwcG为气顶体积内的束缚水饱和度。
8.根据权利要求6或7所述的装置，其特征在于，所述第二获得单元包括：
第一计算子单元，用于根据所述当前地层压力以及第一预设公式，计算获得所述气顶侵入体积；
第二计算子单元，用于根据所述气顶侵入体积以及第二预设公式，计算获得所述油气界面移动速度。
9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5任一项所述方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至5任一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010186156.9A CN113404485B (zh) | 2020-03-17 | 2020-03-17 | 一种凝析气顶油藏的油气界面移动速度的确定方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010186156.9A CN113404485B (zh) | 2020-03-17 | 2020-03-17 | 一种凝析气顶油藏的油气界面移动速度的确定方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113404485A true CN113404485A (zh) | 2021-09-17 |
| CN113404485B CN113404485B (zh) | 2023-08-22 |
Family
ID=77677140
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010186156.9A Active CN113404485B (zh) | 2020-03-17 | 2020-03-17 | 一种凝析气顶油藏的油气界面移动速度的确定方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113404485B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118008229A (zh) * | 2024-03-04 | 2024-05-10 | 西南石油大学 | 一种稠油注蒸汽热采油气界面前缘动态位置的确定方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102953717A (zh) * | 2011-08-26 | 2013-03-06 | 中国石油天然气股份有限公司 | 废弃凝析气藏注水开发方法 |
| WO2014074474A2 (en) * | 2012-11-09 | 2014-05-15 | Saudi Arabian Oil Company | Predicting performance of gas condensate reservoirs |
| CN104389592A (zh) * | 2014-10-08 | 2015-03-04 | 西南石油大学 | 底水带油环凝析气藏水淹层油损失评价实验测试方法 |
-
2020
- 2020-03-17 CN CN202010186156.9A patent/CN113404485B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102953717A (zh) * | 2011-08-26 | 2013-03-06 | 中国石油天然气股份有限公司 | 废弃凝析气藏注水开发方法 |
| WO2014074474A2 (en) * | 2012-11-09 | 2014-05-15 | Saudi Arabian Oil Company | Predicting performance of gas condensate reservoirs |
| CN104389592A (zh) * | 2014-10-08 | 2015-03-04 | 西南石油大学 | 底水带油环凝析气藏水淹层油损失评价实验测试方法 |
Non-Patent Citations (2)
| Title |
|---|
| 刘佳等: "气顶油环协同开发下油气界面运移规律研究", 西南石油大学学报(自然科学版), vol. 37, no. 05, pages 99 - 105 * |
| 张安刚等: "计算凝析气顶油藏油气界面移动距离的新方法", 油气地质与采收率, vol. 22, no. 01, pages 102 - 105 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118008229A (zh) * | 2024-03-04 | 2024-05-10 | 西南石油大学 | 一种稠油注蒸汽热采油气界面前缘动态位置的确定方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113404485B (zh) | 2023-08-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wardlaw et al. | Mercury porosimetry and the interpretation of pore geometry in sedimentary rocks and artificial models | |
| Amaefule et al. | The effect of interfacial tensions on relative oil/water permeabilities of consolidated porous media | |
| CN105160146B (zh) | 水驱特征关系图版生成方法及装置 | |
| CN109884269B (zh) | 岩心自发渗吸采收率预测方法及*** | |
| Cusini et al. | Algebraic dynamic multilevel method for compositional flow in heterogeneous porous media | |
| CN113404485A (zh) | 一种凝析气顶油藏的油气界面移动速度的确定方法及装置 | |
| Dutt et al. | Predicting concentrations of solutes in water percolated through a column of soil | |
| CN103775078A (zh) | 幂律流体在多孔介质中的有效粘度的预测方法 | |
| Pentland et al. | Capillary trapping in water-wet sandstones: coreflooding experiments and pore-network modeling | |
| Singh et al. | Pore-scale Ostwald ripening of gas bubbles in the presence of oil and water in porous media | |
| US20240069577A1 (en) | Methods and devices for dynamic pore network modeling of two-phase flow in mixed-wet porous media | |
| CN110242261B (zh) | 垂向注气孔隙油气水微观渗流规律的预测方法及*** | |
| CN113404486B (zh) | 一种凝析气顶油藏的流体界面移动速度确定方法及装置 | |
| CN110674463A (zh) | 一种确定油水两相相对渗透率曲线的方法和装置 | |
| Al-Shalabi | A new insight into modeling of polymer flooding in carbonate reservoirs | |
| CN112112639A (zh) | 一种凝析气藏循环注气条件下地层压力确定方法及*** | |
| CN107977732A (zh) | 边水气藏水平井的见水时间预测方法及装置 | |
| Larsen et al. | Comparing hysteresis models for relative permeability in WAG studies | |
| CN111434889A (zh) | 带凝析气顶油藏中注入水配比的获取方法及*** | |
| Suekane et al. | Non-wetting phase saturation after drainage from the wetting-phase-filled porous medium | |
| Nomeli et al. | A data driven model for the impact of IFT and density variations on CO2 storage capacity in geologic formations | |
| CN110287545B (zh) | 一种计算油藏内部任意部位原始含水饱和度的方法及*** | |
| CN107989581A (zh) | 一种提高高含水致密凝析气藏采收率的方法及其验证方法 | |
| CN113468476A (zh) | 一种聚合物驱增油量预测方法及装置 | |
| CN112487594A (zh) | 油藏水体倍数计算方法及装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |