CN101366041A - 利用溶气驱重油进行储层模拟的方法、系统和程序存储设备 - Google Patents
利用溶气驱重油进行储层模拟的方法、系统和程序存储设备 Download PDFInfo
- Publication number
- CN101366041A CN101366041A CNA2006800525037A CN200680052503A CN101366041A CN 101366041 A CN101366041 A CN 101366041A CN A2006800525037 A CNA2006800525037 A CN A2006800525037A CN 200680052503 A CN200680052503 A CN 200680052503A CN 101366041 A CN101366041 A CN 101366041A
- Authority
- CN
- China
- Prior art keywords
- correlationship
- reservoir
- capillary
- gas
- relative permeability
- 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
- 238000004088 simulation Methods 0.000 title claims abstract description 68
- 238000000034 method Methods 0.000 title claims abstract description 49
- 239000000295 fuel oil Substances 0.000 title claims abstract description 39
- 230000035699 permeability Effects 0.000 claims abstract description 100
- 239000012530 fluid Substances 0.000 claims abstract description 45
- 210000001736 capillary Anatomy 0.000 claims description 75
- 238000012360 testing method Methods 0.000 claims description 39
- 230000008569 process Effects 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 238000007620 mathematical function Methods 0.000 claims 1
- 239000007789 gas Substances 0.000 description 121
- 239000003921 oil Substances 0.000 description 71
- 239000011435 rock Substances 0.000 description 40
- 230000035508 accumulation Effects 0.000 description 15
- 238000009825 accumulation Methods 0.000 description 15
- 229920006395 saturated elastomer Polymers 0.000 description 10
- 239000003350 kerosene Substances 0.000 description 9
- 230000008859 change Effects 0.000 description 6
- 238000006073 displacement reaction Methods 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 6
- 238000005755 formation reaction Methods 0.000 description 5
- 239000011148 porous material Substances 0.000 description 5
- 230000009467 reduction Effects 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 238000002591 computed tomography Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 230000000670 limiting effect Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 239000004576 sand Substances 0.000 description 3
- 238000010998 test method Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 239000010779 crude oil Substances 0.000 description 2
- 238000007701 flash-distillation Methods 0.000 description 2
- 238000009533 lab test Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000005381 potential energy Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 229920002449 FKM Polymers 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 239000006260 foam Substances 0.000 description 1
- 239000001307 helium Substances 0.000 description 1
- 229910052734 helium Inorganic materials 0.000 description 1
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000006911 nucleation Effects 0.000 description 1
- 238000010899 nucleation Methods 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 238000009418 renovation Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000010729 system oil Substances 0.000 description 1
- 238000003325 tomography 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
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种用于预测流体的性质的方法、系统和程序存储设备,例如用于预测来自包含含气重油的地下储层中的流体生产量。所述方法包括得到气体相对渗透率krg相对于气体饱和度Sg的基线相关关系。与毛细数相关的相关关系基于多个衰竭速度通过获取临界气体饱和度Sge与毛细数Nca之间以及气体相对渗透率krgro与毛细数Nca之间至少一个的关系来确定。为代表地下储层的储层模型中的多个单元计算毛细数Nc。然后调整基线相关关系以与根据与毛细数相关的相关关系选择的Sgc和krgro中的至少一个一致,从而产生多个相应的经过调整的基线相关关系。用于多个单元的气体相关渗透率krg根据多个相应的经过调整的基线相关关系选择。然后利用选定的相对渗透率krg来运行储层模拟以预测在包含含气重油的地下储层中的至少一种流体的性质。
Description
技术领域
本发明总体涉及用于预测地下储层中的流体流的储层模拟的方法和系统,更特别地,涉及通过考虑由于溶气驱重油(heavy oil solutiongas drive)对流体流的影响来提高储层性能预测。
背景技术
储层模拟用于预测地下储层中的流体流。流体流可以包括油、气和水。这种储层预测对于储层管理和估计储层的潜在开采量很重要。
储层模拟在石油工业和科学文献中众所周知。关于储层模拟的原理的一本好的入门书是K.Aziz和A.Settari的Petroleum ReserviorSimulation,Elsevier Applied Science Publishers,London(1979)。另外,在Watts III等人的美国专利No.6,052,520中描述了一般如何进行储层模拟的说明。这些参考文献的全部内容在此引入作为参考。
下列是在传统的储层模拟中所采取的一般步骤。首先,选择要对其岩石和流体性质进行建模和模拟的储层。储层被建模并被离散化成多个单元。一般以有限差分方程的形式为每一个单元构造非线性控制方程,其表示储层中的岩石和流体的性质。岩石性质的实例包括孔隙度、毛细压力和流体(油、水、气)的每一个相的相对渗透率。流体性质的实例包括油的粘度、油地层因数(B0)、以及每一个单元的压力、温度、和饱和度。这些方程中的非线性项被线性化以便形成用于模拟的每一个时间步长的一组线性方程。然后这些线性方程可被求解以估算未知项(比如单元中的压力和饱和度)的解值。从压力和饱和度的这些值可以估算出其他性质,其包括每个时间步长的来自储层的油、气和水的总生产量。重复进行前述步骤,经过许多这种时间步长来模拟储层中相对于时间的流体流。
储层模拟中所需的关键性质中的一个是待流经的岩石的渗透率。绝对渗透率K是岩石传送流的能力的量度,而且可能在储层和周围地层各处变化很大。当气、油和水移动穿过多孔岩石时,它们不会等速运动。而是流体彼此竞相运动。当存在多于一种流体时,相对渗透率kr是有效渗透率ke与绝对渗透率K的比率。有效渗透率ke是当存在另一个流体时,多孔介质对一种流体所测量的渗透率。相对渗透率kr和饱和度S之间的关系取决于储层岩石和流体,而且在不同地层之间可以改变。此外,相对渗透率kr取决于存在的流体的相对比例,即流体饱和度。
图1说明了气体的典型的相对渗透率krg相对于饱和度Sg的曲线。直到气体饱和度达到最小阈值之前,气体不能以任何可感知的速度流动。参照图1,该阈值被称为临界气体饱和度并且起始于近似0.03或大约3%饱和度的值。曲线的另一端是端点相对渗透率其为保留在岩石中的残余油的运动最小时的气体相对渗透率值krg。由于储层岩石总是包含最小量的残余油,气体饱和度不会达到100%。饱和度的总百分比必须加起来达100%。在这种情况下,存在最大的76%气体饱和度Sg和24%的残余油饱和度Sorg。从图1中可看出,最大的相对渗透率发生在接近0.76的饱和度处,这时kr=0.40。和的这些值将被称为气体饱和度Sg和相对渗透率krg的端点基线值。
理想地,通过对取自将进行储层模拟的储层的岩心试样进行实验室试验来得到相对渗透率曲线。例如,驱替试验可以用于得到相对渗透率krg相对于饱和度Sg的曲线。这种试验众所周知。特别地,在文献E.F.Johnson、D.P.Bossier和V.O.Nauniann的Calculations ofRelative Permeability from Displacement Experiments,Trans.Am.Inst.Mining Engineers,Volume 216,1959,pp.370-378以及S.C.Jones和W.O.Roszelle的Graphical Techniques for DeterminingRelative Permeability from Displacement Experiments,Journal ofPetroleum Engineering,Volume30,pp.807-817(1978)中描述了公知的驱替试验程序。由于通常认为渗透率曲线基本与气体流动穿过储层岩石的速度多快无关,因此这些驱替试验通常在低衰竭速度下进行。
替代地,如果得不到岩心试样,可以从理论上建立相对渗透率krg相对于饱和度Sg的曲线。例如,曲线可以根据可比的类似储层得到。
一旦已经获得相对渗透率krg相对于饱和度Sg曲线,假定储层模型单元中的饱和度Sg已知,则从这些曲线就能简单地获得将在储层模拟中使用的相对渗透率krg。饱和度Sg一般从在模拟的起始点处建立的初始条件、从模拟的最后的时间步长或者从时间步长的迭代中的计算得知。
重油的生产量最初主要由油的压力驱动。可以认为重油包括具有200或更小的API重力值的油。当在高储层压力下时,大量的气体常常存在于重油内。在从储层生产出足够量的重油之后,储层部分的压力可以下降到低于泡点压力。在该压力,气体容易从重油溶液中逸出。一旦足够的气体已经从油中释放出来,则认为气体形成连续相,而且气体可以流动穿过储层,而且气体的生产率显著提高。如上所述,气体流动开始处的饱和度Sg称为临界气体饱和度或者Sgc。图11显示了从岩心试样产生的累积气体相对于时间(分钟)的曲线图。曲线中所示的断点表示Sgc。
试验表明来自重油储层的油开采量取决于储层的衰竭速度。较高的衰竭速度常常将导致整体提高的油开采量。由于不能很好地理解溶气驱重油的机制,储层模拟器通常利用静态气体相对渗透率krg相对于饱和度Sg的曲线,比如从图1中看到的曲线,其与流体流或衰竭速度不相关。一旦用于待建模的相应类型的岩石的这些曲线建立,不管穿过储层单元的流速如何,曲线将在整个储层模拟过程中保持相同(即和的端点保持固定)。对于包含非重油的含烃地下地层中的一般储层模拟而言渗透率曲线是稳定的假设一般能得到满足。
然而,在重油的情况下,不平衡的溶气驱油(“泡沫油”)是影响临界气体饱和度Sgc和油开采量的重要生产机制。当前,对于所有规模(气孔、岩心和矿区)的溶气驱重油的理解是有限的。传统的储层模拟器在预测储层中的流体流时不能准确地考虑这种溶气驱油。这是常常导致低估重油产量的预测的显著缺点。本发明通过考虑了溶气驱重油的影响来克服这一缺点。
发明内容
本发明公开了一种用于预测在包含含气重油的地下储层中的至少一种流体的性质的方法。例如,该性质可以包括来自储层的流体(即油、气和水)的总生产量。使用储层模拟器进行预测,该储层模拟器使用具有代表储层的多个单元的储层模型。对于至少一些单元以及对于至少一些储层模拟的迭代,气体相对渗透率krg取决于单元中的局部流体速度vo。
在本方法的一个优选实施例中,通常基于在慢速的衰竭速度下进行的驱替试验,来得到气体相对渗透率krg相对于气体饱和度Sg的基线相关关系。接着,在临界气体饱和度Sgc与毛细数Nca之间和气体相对渗透率端点Krgro与毛细数Nca之间的至少一个得到与毛细数Nca相关的相关关系,而且最优选在临界气体饱和度Sgc与毛细数Nca之间和气体相对渗透率端点Krgro与毛细数Nca之间都得到与毛细数Nca相关的相关关系。可如何表示这种相关关系的非限制性实例包括(举例来说但不局限于)使用描述曲线的数学方程或者建立相应的查阅表。
然后,可以使用这些通过实验方法获得的与毛细数Nca相关的相关关系结合储层模拟来获得溶气驱重油和衰竭速度对重油和含在重油中的气体的生产的影响。为代表其流体性质待模拟的地下储层的储层模型中的多个单元计算毛细数Nc。Sgc和/或Krgro的值基于为单元计算的毛细数Nc根据与毛细数相关的相关关系选择。然后得到经过调整的基线相关关系。例如,基线曲线的初始端点,即和用新的与毛细数相关的Sgc和Krgro值替代,其间的曲线比如通过线性比例缩放而调整。图2建议经过调整的基线曲线可以通过将初始端点值和改变为其他的Sgc和Krgro值而得到,其部分地基于油流动穿过单元的速度vo。
多个单元的气体相对渗透率krg根据相应的经过调整的基线相关关系选择。然后,在储层模拟中使用这些相对渗透率krg来预测在包含含气重油的地下储层中的至少一种流体的性质。这种预测的性质可以是油、水或气的生产量。优选地,一旦单元中的饱和度Sg等于或大于临界气体饱和度Sgc的水平,该单元的当前经过调整的基线相关关系在剩余的模拟时间步长中被固定。一旦气体开始流动,该经过调整的基线相关关系的固定有助于在求解用于为储层建模的方程组的过程中保持稳定性。
可以使用Sgc或Krgro的与毛细数相关的相关关系中的一个或两个来调整基线相关关系以提出经过调整的基线相关关系。这些经过调整的基线相关关系,通过使用毛细数Nc,获得衰竭速度/流体速度流和粘度在溶气驱重油作用下的重油生产期间对相对渗透率的影响。优选地,在不同衰竭速度下进行衰竭试验以便得到Sgc和Krgro的与毛细数相关的相关关系。然而,如果需要,可通过理论方法预测所述与毛细数相关的相关关系应该是什么样。
可以选择在储层模拟中的时间步长开始处所计算的与毛细数相关的相对渗透率krg。替代地,毛细数Nc可以在时间步长中的迭代期间一直重复计算,以便在模拟期间提供相对渗透率曲线的不断更新。此外,一旦单元中的饱和度Sg在模拟期间保持在临界气体饱和度Sgc处或之上,则优选停止单元的相对渗透率曲线的毛细数Nc的这种更新。
本发明的一个目的是通过比传统的储层模拟器更好地考虑由于溶气驱重油对流体流的影响来提高储层性能预测,从而改进地下地层中涉及重油流的储层模拟的预测能力,这可能导致改进储层管理策略。
另一个目的是用实验方法在许多不同的衰竭速度下确定岩心试样的临界气体饱和度Sgc值和/或气体相对渗透率端点Krgro的值,以及使这些值与毛细数Nca相关而得到与毛细数相关的相关关系。这些与毛细数Nca相关的相关关系可以结合储层模型和在储层模拟过程中所计算的算出的毛细数Nc一起使用,以便更准确地估计在重油的储层模拟中要使用的相对渗透率krg。
附图说明
根据下列说明、所附的权利要求书和附图,将更好地理解本发明的这些和其他目的、特征和优点,其中:
图1显示了传统的气体相对渗透率krg相对于饱和度Sg的曲线;
图3显示了在本发明的一个优选实施例中进行储层模拟所采取的步骤的流程图,其利用了取决于单元中的流体流的局部速度vo的气体相对渗透率krg;
图4显示了用于确定岩心和填砂柱(sandpack)试样的气体饱和度Sg的试验装置的示意图;
图5描绘了在快速衰竭试验中横跨填砂柱试样的平均填砂柱压力和压差相对于时间的曲线图;
图6示出了在图5的快速衰竭试验中所产生的累积的油和气的曲线图;
图7显示了慢速衰竭试验的平均填砂柱试样压力和流出物(effluent)密度的曲线图;
图8描绘了慢速衰竭试验的平均填砂柱试样压力和所产生的累积的油的曲线图;
图9显示了油开采量作为在衰竭速度为0.3和0.03cc/min下所进行的填砂柱试验的平均孔隙压力的函数的曲线图;
图10是油开采量作为在衰竭速度分别为0.082、0.08和0002cc/min下所进行的岩心试验的平均孔隙压力的函数的曲线图;以及
图11是产生的累积气体(测得的)和产生的累积溶气(计算的)相对于时间的曲线图。
具体实施方式
I.引言
本发明考虑了溶气驱重油的影响,而且更特别地考虑了流体衰竭速度对重油生产的影响。在重油储层模拟中利用与速度或衰竭速度相关的相对渗透率krg值来提供比用传统的储层模拟所获得的储层模拟预测更准确的储层模拟预测。
在一个优选实施例中,为储层单元计算与油速度vo相关的毛细数Nc。这些毛细数Nc用于调整基线相对渗透率相关关系,以考虑速度或衰竭速度对相对渗透率krg的影响。在本优选实施例中,首先,优选基于实验室试验,得到与毛细数Nca相关的临界气体饱和度Sgc和/或端点相对渗透率krgro的相关关系。然后,利用相应于为单元计算的毛细数Nc的Sgc和/或krgro的值来调整该单元的基线相对渗透率的相关关系。然后基于单元中的饱和度Sg由这些毛细数调整的基线相对渗透率的相关关系选择相对渗透率krg的值。
图3提供了可以用于进行本发明的溶气驱重油储层模拟的步骤的示例性流程图。在步骤100中,得到krg和Sg之间的基线相关关系。然后在步骤110中得到Sgc和Nca和/或krgro和Nca之间的相关关系。在步骤120中为储层模型中的许多单元计算毛细数Nc。在步骤130中,为这些单元中的每一个建立krg和Sg之间的经过调整的基线相关关系,其取决于Nc和在步骤110中得到的相关关系。然后在步骤140中使用来自单元的饱和度Sg值,从krg和Sg之间的经过调整的基线相关关系为每一个单元选择气体相对渗透率krg。然后在步骤150中,在储层模拟中使用与这些毛细数相关的渗透率krg来预测储层模型中的流体流的性质。
将描述用于建立Sgc和Nca之间以及krgro和Nca之间的相关关系的示例性试验方法的说明。然后,将描述对传统的储层模拟器所进行的修改,通过所述修改来在进行储层模拟时结合衰竭速度/与毛细数相关的Sgc和/或krgro相关关系来选择相对渗透率krg。
II.建立相关关系
A.基线气体相对渗透率k rg 相对于饱和度Sg 的相关关系
建立气体相对渗透率krg与饱和度Sg之间的相关关系,以使得基于储层模型单元中的已知的饱和度Sg值,相对渗透率值k rg 可以被储层模拟器利用。理想地,这些相关关系由来自将进行储层模拟的储层的岩心试样通过试验方法得到。替代地,也可以利用典型的填砂柱和/或合成油来得到相关关系。用于建立这些基线相关关系的优选方法是Johnson、Bossier和Naumann的方法,或者Jones和Roszell的方法,它们在上述的背景技术部分中被引用并且被建立渗透率曲线的技术人员所公知。替代地,存在用于为储层岩石和流体建立气体相对渗透率krg相对于饱和度Sg曲线的许多其他公知的方案。通常,krg将成为Sg的函数。由于实际原因,人们常常将用于krg相关关系中的气体饱和度归一化。在方程12中描述了这种归一化。这种归一化允许模拟器容易地估算用于改变端点(例如Sgc和Sorg)的krg。
如果得不到岩心试样,则相对渗透率krg与饱和度Sg之间的相关关系可以用理论方法估计。作为非限制性的实例,可利用相似地层来最初建立基线曲线。相关关系的非限制性的实例可以采用几种形式,比如曲线、数学表达式、查阅表等。
图1是气体相对渗透率krg相对于饱和度Sg的示例性基线曲线或相关关系。所示的的基线值位于大约0.03或3%处。
在该值之上,期望气体开始自由流动而不是主要被截留在多孔介质内。最大气体饱和度Sg为大约76%,并且残余油饱和度Sorg为24%饱和度。本实例假设几乎不存在水。在最大气体饱和度Sg=76%处,最大气体相对渗透率为近似0.4%。
B.Sgc相对于Nca和krgro相对于Nca之间的相关关系
在不同衰竭速度下进行实验室试验以便建立Sgc相对于Nca和krgro相对于Nca的相关关系。Sgc在下述方法中获得。使用下面的方程(8)计算Nca。从试验和使用对岩心或填砂柱试样的储层模拟的历史匹配(history matching)获得每一个衰竭速度的Sgc、krgro和Nca的值。然后,通过曲线拟合Sgc、krgro和Nca的数据而获得Sgc与Nca和krgro与Nca之间的相关关系。岩心试样的生产数据的历史匹配可以用于提高相关关系的准确度。
1.含气石油的制备
通过使未过滤的脱气石油(dead oil)与甲烷组合来制备含气石油(live oil)。油的水含量可忽略。PVT(压力、体积和温度)数据:油气比(Rs)、油层体积因子(Bo)和气层体积因子(Bg)通过试验(恒定的成分膨胀、闪蒸(flash)、密度测量)与状态方程的调整相组合来确定。含气石油的粘度在储层温度下在毛细管粘度计(ID=0.05in)中进行测量。表1列出了在178℉下的含气石油的相关关系质。
表1 原油的性质
泡点压力(Psia) | 1350 |
溶液的油气比(cc/cc) | 20 |
在泡点压力下的Bo | 1.0918 |
含气石油的粘度(cp) | 240 |
脱气石油的粘度(cp) | 1300 |
2.衰竭试验
在恒定衰竭速度下在水平80cm长的填砂柱或者29cm的水平合成岩心(4塞)中进行衰竭试验。填砂柱试验中使用的砂是尺寸范围从75到125微米的清洁的Ottawa砂。砂被包装在配有压力端口的特制维托套管(Viton sleeve)中。用氦孔隙度仪测量填砂柱和合成岩心的孔隙度。填砂柱和合成岩心的性质列在表2中:
表2 填砂柱和合成岩心的性质
表2 | 填砂柱 | 合成岩心 |
温度,℉ | 178 | 178 |
长度,cm | 80 | 29 |
直径,cm | 5.04 | 5.04 |
上覆岩层压力,psia | 2050 | 2050 |
孔隙度 | 0.33 | 0.27 |
孔隙体积,cm3 | 560 | 162 |
含气石油的渗透率,md | 2000 | 1850 |
衰竭速度范围,cm3/min | 0.002到0.3 | 0.0003到0.03 |
使用以再填充模式操作的一个或两个ISCO泵来控制衰竭速度。图4显示了试验装置的示意图。在衰竭期间,监控压力(入口、出口和沿着岩心的几个点处)、油和气的生产量以及流出物的密度。将岩心夹持器放置在Siemens Somatom HiQ CT扫描仪上以监控空间的和时间的气体饱和度。
3.程序
干填砂柱首先在储层条件(即在上覆岩层应力和温度下)下被进行CT(计算机X射线层析摄影法)扫描。然后岩心在~1600psia的背压下被用CO2闪蒸、被抽空并用煤油饱和。填砂柱(或合成岩心)的渗透率在不同的流速下用煤油测量。充满煤油的填砂柱也被进行CT扫描。使用湿的和干的空气和煤油的CT扫描和CT数来计算填砂柱的孔隙度。然后将含气石油慢慢注入到岩心中来替代煤油。也在若干流速下用含气石油测量填砂柱渗透率。然后减少含气石油的注入率使得横跨岩心的差压小于2psi。
被含气石油饱和的填砂柱被进行CT-扫描以便记录初始状态。在~1500-1700psi(在泡点压力之上大约150-350psi)的压力下开始衰竭。关闭入口阀,并且以恒定的抽出速度操作下游的Isco泵A。在给定衰竭时间之后,转换泵并且Isco泵B抽出流体而Isco泵A将油和气传送到收集系统中。重复进行泵循环直到出口压力减小至大约200psia。使用传统的数据采集软件连续记录收集系统中的压力、温度和流体积累。使用在线密度计连续测量产生的流体的密度。还周期地扫描填砂柱以直接确定作为时间和位置的函数的气体饱和度Sc。
4.CT-扫描
Siemens Somatom HiQ CT扫描仪用于监控空间的和时间的气体饱和度。这种第三代的CT-扫描仪具有768个固定的探测器和旋转的X射线源。以133kV进行扫描并且扫描时间为2.7秒。三维像素尺寸对于10mm的扫描厚度来说为近似0.625mm3而且饱和度测量中的不确定度为+/-1.5饱和度单位。获得了10mm和/或5mm的扫描厚度。
5.结果
在试验的过程中,获得了沿着岩心和在封闭的岩心入口和开放的岩心出口处的压力信息、产生的油和气的量、流出物密度和气体饱和度(经由CT-扫描仪)。在试验期间观察到的典型响应如图5和图6所示。图5显示了在快速衰竭试验期间的平均填砂柱压力和横跨填砂柱的压差。图6示出了在快速衰竭试验期间产生的累积的油和气。
虽然不希望依据特别的理论,相信在早期仅通过油和地层膨胀进行生产(系统中没有游离气体),而且压力迅速降低。在(明显的)泡点压力下气泡开始成核。当压力减小到泡点压力之下时,气泡的尺寸缓慢增长而且油的生产由气体膨胀控制。如从图5可看到的,压力减小率显著下降。油是唯一的运动相,而且收集的气体是由收集系统中溶解的气体所释放的。在临界气体饱和度Sgc处,气泡在填砂柱各处连接而且气体开始自由流动。注意气体的生产量有显著增加而油的生产量渐渐减少(见累积的气生产量图中在~270分钟时的急剧突变)。
对于填砂柱中的较慢衰竭速度试验和对于岩心试验,也测量流出物密度。图7和8显示了在这种仪器中观察到的典型响应。图7示出了慢速衰竭试验的平均填砂柱压力和流出物密度。图8描绘了慢速衰竭试验的平均填砂柱压力、产生的累积的油(收集在分离器中以及基于流出物密度推断的)。
6.速度影响
在衰竭试验期间观察到的主要影响是油开采量对衰竭速度非常敏感。大填砂柱试验(图9)和小岩心试验(图10)都观察到这种现象。图9示出了油开采量为平均孔隙压力(填砂柱试验-速度=0.3和0.03cc/min)的函数。图10显示了油开采量作为平均孔隙压力(岩心试验-速度=0.082、0.08和0.002cc/min)的函数。
除了速度影响之外,注意在这些试验中观察的总油开采量非常大(达~30%OOIP)。如此高的开采量和这种对衰竭速度的相关关系不容易用传统的物理学来解释。此外,这种现象用当前商用模拟器不能进行正确建模。
7.数据分析-S g 和S gc 的确定
临界气体饱和度Sgc是在产生的累积气体开始显著增加处的饱和度。图11显示了产生的累积气(测量的)和产生的累积溶气(计算的)相对于时间的关系。临界气体饱和度Sgc也可以基于流出物密度来确定。
用图4所描述的装置,有几种方式来确定气体饱和度:
(1)直接在原地用CT-扫描仪测量;
(2)使用收集在收集系统中的流体量的物质平衡;以及
(3)使用流出物流的密度的物质平衡。
方法2和3需要使用PVT数据(即地层体积因子和作为压力函数的密度)。
物质平衡:
SR=1-So (1)
其中:N是试验开始的适当位置处而且在压力Pi下的油(stb);Np是在压力P下产生的累积油(stb)(Np是用收集系统测量的);Bo和Boi分别是在P和Pi下的油层体积因子;以及cf是岩石或填砂柱的可压缩性(1/psi)。
在泡点之上,油仅通过油和地层膨胀产生。即
其中,油的可压缩性给定为
由于co已知,填砂柱和合成岩心的可压缩性使用方程(3)来计算。
如上所指出,Np是通过收集系统来测量的。替代地,产生的油的量可基于流出物密度,ρeff:
孔隙度和气体饱和度都可使用CT-扫描仪来计算。孔隙度被给定为
其中:CTsaturated_core是被煤油饱和的填砂柱(在初始压力下)的CT数,CTdry_core是被气体饱和的填砂柱的CT数。CTliq和CTgas分别是煤油和空气的CT数。
类似地,气体饱和度通过下列方程获得:
其中:CTp是在衰竭期间(在压力P下)测得的CT数.CTsaturated_core是被含气石油饱和(在初始压力下)的填砂柱的CT数,以及CTdry_core是被空气饱和的填砂柱在初始压力下的CT数。
8.数据分析-毛细数的计算
对于每一个试验,使用在衰竭期间记录的压差来计算平均毛细数(Nca)。毛细数可用几种方式来计算。在本优选实施例中,使用下列公式:
其中:K是岩心或填砂柱的渗透率,σ是气油表面张力(在试验中使用的油被估算为80dyn/cm),L是填砂柱的长度,以及ΔP是在气体变得可动之前所测得的压差。
9.数据分析—作为N ca 的函数的S gc 和k rgro
基于上述分析,对于所有可用试验,Sgc被绘制成N ca 的函数的曲线。然后对数据进行曲线拟合,优选地,使用指数函数(方程(9))来内插/外插缺少的数据。系数“a”和指数“b”的值针对每一个油/岩石系统而异。
举例来说式,但不局限于此,优选的作为Nca的函数的Sgc和krgro之间的数学相关关系性如下:
以及
在不同衰竭速度下对岩心试样进行的储层模拟用于确定krgro的值。对于每一个模拟操作,临界气体饱和度Sgc是已知的,因此气体相对渗透率krg相对于饱和度Sg曲线上的该端点是已知的。对于曲线krgro的另一个端点,做出多个估算值。然后,使用反复试验方法来确定在特定的衰竭速度下krgro的估算值中的哪一个值与由岩心试样产出的试验生产量相匹配。使用利用多次模拟得到的试验生产量结果的历史记录匹配来在多个衰竭速度下确定krgro,其与Nca值对应。然后这些krgro与Nca的值被进行曲线拟合以便得到与毛细数相关的相关关系。最优选地,该相关关系为带有确定的“c”和“d”值的方程(10)形式。
III.利用溶气驱重油进行储层模拟
从试验数据获得的Sgc和krgro相对于Nca的函数形式在本示例性实施例中执行,优选在储层模拟器中使用修正的隐式算法执行。举例来说,但不局限于此,Sgc和krgro的优选形式是使用上述的方程(9)和(10)作为Nca函数的输入量。参数a、b、c和d是使用者输入到储层模拟器的输入量。应注意的是,在图2中Sgc是a、b和毛细数Nca的函数。类似地,krgro是c、d和Nca的函数。在本发明的优选实施例中,下列是默认值:a=l04;b=1.0;c=l04;以及d=1。理想地,所计算的Sgc和krgro的值分别受使用者指定的最大值和最小值的限制。例如,可使用最大值Sgc=0.1和最小值krgro=10-4。由于Nc是定向的,为每个单元面计算Sgc和krgro,而且从而它们也是定向的。
为了减少振荡和收敛问题,执行优选实施例的修正的隐式算法来计算Sgc和krgro。当气相是不可动的,即饱和度Sg≤Sgc时,例如分别使用方程(9)和(10)来计算Sgc和krgro。当气相流动时,Sgc和krgro变成不变量—既不增加也不减小。它们的值使用当气相变得可动在时间步长的起点处的毛细数Nc来计算,并且对于所有剩余的时间步长都固定。
A.单元毛细数Nc的计算
在优选示例性实施例中,毛细数Nc的修正表达式优选地使用下列表达式来并入储层模拟器中:
其中:σog为油气界面张力,K为岩石渗透率,φo为油相势能,Po是横跨单元面的压力变化,ρo=油的密度,g=重力常数,以及D=离基准点的深度变化。
Nc的该修正定义在方程中省去油相对渗透率。由于Nc理想地以隐式计算,这在雅可比行列式(Jacobian)生成过程中大大简化了作为原始变量的函数的气体相对渗透率的导数的计算。同样,Nc计算中的势梯度是定向的,并且基于横跨两个相邻网格块的面的梯度。对于每一次牛顿迭代,为每一个网格块的面计算毛细数Nc。在三维模型中,每一个网格块将有六个方向的Nc。每一个Nc对应于单元面处的六个值中的一个。定向的Nc的使用导致了雅可比行列式可容易地通过传统的线性方程解算器进行求解。在本优选实施例中,计算了所有网格块的面的平均Nc。
B.调整基线相对渗透率相关关系
为每一个单元分配特定的岩石类型或相。这些岩石类型或相的每一个对应于特定基线气体相对渗透率krg相对于饱和度Sg的曲线,比如图1中所示的曲线。为每一个相应的单元调整这些相应的基线曲线。这通过由使用方程(9)和(10)以及为每一个相应单元计算的特定毛细数Nc而计算得到的与毛细数相关的Sgc和krgro值来替换和的初始值来为每一个单元实现。连接这些端点的基线曲线优选通过按比例缩放而进行调整。可使用几种方法确定相对渗透率的比例。方程(12)显示了这样一种方法,即:
方程(12)简单陈述了krg是Sg、Sgc和Sorg的函数。(气体饱和度大于Sorg时,油相是不可动的—即,Kro=0)。函数F可为(但不局限于)一个简单的幂函数:
另外,如果气体相对渗透率krgro的端点相对于基线曲线的初始减小10%,则相关关系或曲线上的所有气体相对渗透率值将减小10%。所属领域的技术人员应该理解的是可使用许多其他方式来调整基线曲线以便反映端点Sgc和/或krgro的更新值的变化,而且也在本发明范围内。
C.选择用于并入储层模拟器中的气体相对渗透率k rg
饱和度值Sg可来自当储层模拟第一次开始时的初始条件,来自前一个时间步长,或来自在时间步长内迭代计算的值。然后,检查每一个储层单元的饱和度Sg,并且由经过调整的基线相关关系选择相应的相对渗透率krg。如上所述,如果Sg≥Sgc,则使用来自先前计算的曲线相关关系来确定krg。
D.使用选定的气体相对渗透率k rg 来运行储层模拟
求解有限差分方程来确定未知量,比如压力P或饱和度Sg。这些有限差分方程取决于最近更新的相对渗透率kr,包括储层单元的与毛细数相关的气体渗透率krg。这种有限差分方程对于储层模拟领域的技术人员来说是公知的。该方程公知的求解方法的实例包括:(1)完全显式的;(2)隐式压力、显式饱和度(IMPES);(3)完全隐式的;(4)顺序隐式的(SEQ)、适应性隐式的(AIM);以及(阶式)Cascade。在优选实施例中,使用完全隐式方法来求解这些方程。
如果状态变量(即压力或饱和度)的解在迭代期间在满意的公差范围内,则将为一个时间步长建立最终的流体性质。在时间步长期间的气、水和油的生产容积可根据这些流体性质建立,如使用储层模拟器的通常做法一样。然后储层模拟器可以运行经历更多的时间步长直到满足预定的时间长度。这些时间步长累积的生产量提供了对地下地层生产的估计。
本发明还包括用于使用取决于衰竭速度/流体速度和原油的粘度的相对渗透率krg进行上述储层模拟的系统。此外,本发明还包括带有用于使用与流体速度相关的相对渗透率来进行该储层模拟的指令的程序存储设备。
尽管在前述说明中,已经结合某些优选实施例描述了本发明,并且为了示例说明目的已经陈述了许多详细资料,但对所属领域的技术人员而言显而易见的是,本发明容许改变而且在此描述的某些其他详细资料可显著改变而不偏离本发明的基本原理。
术语
a=用于计算Sgc的系数;
b=用于计算Sgc的指数;
Boj=在压力Pi下的油层体积因子;
Bo=在压力P下的油层体积因子;
Bg=在压力P下的气层体积因子;
c=用于计算krgro的系数;
cf=岩石或岩心试样的可压缩性(1/psi);
co=油试样的可压缩性(1/psi);
CTdry_core=被气体饱和的试样的CT数;
CTsaturated_core=被煤油饱和(在初始压力下)的试样的CT数;
CTp=在压力P下衰竭期间测量的CT数;
CTliq=煤油的CT数;
CTgas=空气的CT数;
d=用于计算krgro的指数;
D=离基准点的深度的变化;
g=重力常数;
ke=有效渗透率;
kr=相对渗透率,无量纲;
krg=气体相对渗透率,无量纲;
krgro=在最小残余油情况下的气体相对渗透率,无量纲;
kro=油相对渗透率,无量纲;
K=岩石渗透率;
K=溶气曲线的斜率,psi-1;
Nc=为储层模型的特定单元计算的毛细数;
Nca=毛细数;
ΔP=压力变化(psi);
L=试验腔的长度(英寸);
N=在初始条件下适当位置的油(stb);
Np=压力P(cm3)下产生的累积的油(stb);
φo=油相势能;
Pi=在时间i时的压力,psi;
ΔPo=横跨面的压力变化;
ρeff=有效密度;
ρg=气的密度;
ρo=油的密度;
Rs=气油比;
S=饱和度,无量纲;
Sg=气体饱和度,无量纲;
Sgc=临界气体饱和度,无量纲;
So=油饱和度,无量纲;
Sorg=特定岩石区的残余油对气体的饱和度,无量纲;
stb=油罐桶数;
σ=界面张力;
σog=油气的界面张力;以及
vo=油的速度。
Claims (16)
1.一种用于预测在包含含气重油的地下储层中的至少一种流体的性质的方法,所述方法包括:
a)在储层模拟器中利用气体相对渗透率krg相对于气体饱和度Sg的基线相关关系;
b)在储层模拟器中利用临界气体饱和度Sgc与毛细数Nca之间以及气体相对渗透率端点krgro与毛细数Nca之间至少一个的与毛细数相关的相关关系;
c)为代表地下储层的储层模型中的多个单元计算毛细数Nc;
d)使用在步骤c)中计算的毛细数Nc来调整基线相关关系以与根据步骤b)的与毛细数相关的相关关系选择的Sgc和krgro中的至少一个一致,从而产生多个相应的经过调整的基线相关关系;
e)从步骤d)的相应的经过调整的基线相关关系中为多个单元选择相对渗透率krg;以及
f)利用步骤(e)的选定的相对渗透率krg来运行储层模拟以预测在包含含气重油的地下储层中的至少一种流体的性质。
2.如权利要求1所述的方法,其中
利用与毛细数相关的相关关系的步骤包括利用Sgc与Nca之间的相关关系。
3.如权利要求1所述的方法,其中
利用与毛细数相关的相关关系的步骤包括利用krgro与Nca之间的相关关系。
4.如权利要求1所述的方法,其中
利用与毛细数相关的相关关系的步骤包括利用Sgc与Nca之间的相关关系以及利用krgro与Nca之间的相关关系。
5.如权利要求4所述的方法,其中
调整基线相关关系的步骤使得与根据Sgc与Nca之间的相关关系选择的Sgc的值以及根据krgro与Nca之间的相关关系选择的krgro的值一致。
6.如权利要求1所述的方法,其中
利用与毛细数相关的相关关系的步骤包括对来自地下储层的岩心试样进行衰竭试验。
7.如权利要求1所述的方法,其中
利用与毛细数相关的相关关系的步骤包括预测相关关系而不对来自地下储层的岩心试样进行衰竭试验。
8.如权利要求1所述的方法,其中
得到与毛细数相关的相关关系的步骤包括对填砂柱试样进行衰竭试验。
9.如权利要求1所述的方法,其中
为多个单元计算的毛细数Nca在时间步长迭代中保持固定。
10.如权利要求1所述的方法,其中
为多个单元计算的毛细数在储层模拟中进行的时间步长的迭代期间被更新。
11.如权利要求1所述的方法,其中
与毛细数相关的相关关系是查阅表。
12.如权利要求1所述的方法,其中
与毛细数相关的相关关系是数学函数。
13.如权利要求1所述的方法,其中
储层模拟器使用完全隐式方法来求解方程。
14.一种用于模拟地下储层中的重油流的方法,所述方法包括:
创建代表地下储层的储层模型,所述地下储层的流体流待模拟,所述储层模型包括多个单元;
为储层单元确定与速度相关的相对渗透率;以及
利用与速度相关的渗透率运行储层模拟以便模拟地下储层的重油流。
15.一种程序存储设备,所述设备包含用于进行储层模拟方法的指令,所述方法包括:
带有用于储层模拟方法的计算机指令的有形介质,包括以下步骤:
a)在储层模拟器中利用气体相对渗透率krg相对于气体饱和度Sg的基线相关关系;
b)在储层模拟器中利用临界气体饱和度Sgc与毛细数Nca之间以及气体相对渗透率krgro与毛细数Nca之间至少一个的与毛细数相关的相关关系;
c)为代表地下储层的储层模型中的多个单元计算毛细数Nc;
d)使用在步骤c)中计算的毛细数Nc来调整基线相关关系以与根据步骤b)的与毛细数相关的相关关系选择的Sgc和krgro中的至少一个一致,从而产生多个相应的经过调整的基线相关关系;
e)从步骤d)的相应的经过调整的基线相关关系中为多个单元选择相对渗透率krg;以及
f)利用步骤(e)的选定的相对渗透率krg来运行储层模拟以预测在包含含气重油的地下储层中的至少一种流体的性质。
16.一种用于预测在包含含气重油的地下储层中的至少一种流体性质的系统,所述系统包括:
a)用于将气体相对渗透率krg相对于气体饱和度Sg的基线相关关系输入储层模拟器的装置;
b)用于将基于多个衰竭速度的临界气体饱和度Sgc与毛细数Nca之间以及气体相对渗透率端点krgro与毛细数Nca之间至少一个的与毛细数相关的相关关系输入储层模拟器的装置;
c)用于为代表地下储层的储层模型中的多个单元计算毛细数Nc的装置;
d)用于调整基线相关关系以与根据步骤b)的与毛细数相关的相关关系选择的Sgc和krgro中的至少一个相一致的装置,所述Sgc和krgro中的至少一个和在步骤c)中计算的毛细数Nc相对应从而产生多个相应的经过调整的基线相关关系;
e)用于根据步骤d)的相应的经过调整的基线相关关系为多个单元选择相对渗透率krg的装置;以及
f)用于利用步骤e)的选定的相对渗透率krg来运行储层模拟以预测在包含含气重油的地下储层中的至少一种流体的性质的装置。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US75390005P | 2005-12-22 | 2005-12-22 | |
US60/753,900 | 2005-12-22 | ||
PCT/US2006/049095 WO2007076044A2 (en) | 2005-12-22 | 2006-12-20 | Method, system and program storage device for reservoir simulation utilizing heavy oil solution gas drive |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101366041A true CN101366041A (zh) | 2009-02-11 |
CN101366041B CN101366041B (zh) | 2011-10-12 |
Family
ID=38218675
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2006800525037A Expired - Fee Related CN101366041B (zh) | 2005-12-22 | 2006-12-20 | 利用溶气驱重油进行储层模拟的方法 |
Country Status (4)
Country | Link |
---|---|
CN (1) | CN101366041B (zh) |
BR (1) | BRPI0620170A2 (zh) |
CA (1) | CA2634757C (zh) |
WO (1) | WO2007076044A2 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102576374A (zh) * | 2009-09-02 | 2012-07-11 | 界标制图有限公司 | 用于油气形成建模的系统和方法 |
CN102713683A (zh) * | 2010-02-26 | 2012-10-03 | 雪佛龙美国公司 | 用于利用多点统计模拟来对储层特性走向建模的方法和系统 |
CN103590807A (zh) * | 2012-08-13 | 2014-02-19 | 中国石油天然气集团公司 | 一种确定稠油热采蒸汽腔形态的方法 |
CN104481504A (zh) * | 2014-11-11 | 2015-04-01 | 中国石油天然气股份有限公司 | 稠油溶解气驱开发模拟实验系统及其方法 |
CN107101773A (zh) * | 2017-04-14 | 2017-08-29 | 中国石油天然气股份有限公司 | 泡点压力值和拟泡点压力值的测试方法、及测试装置 |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2717353C (en) | 2008-04-22 | 2016-06-21 | Exxonmobil Upstream Research Company | Functional-based knowledge analysis in a 2d and 3d visual environment |
US8892407B2 (en) | 2008-10-01 | 2014-11-18 | Exxonmobil Upstream Research Company | Robust well trajectory planning |
US8849640B2 (en) | 2008-11-06 | 2014-09-30 | Exxonmobil Upstream Research Company | System and method for planning a drilling operation |
US8931580B2 (en) | 2010-02-03 | 2015-01-13 | Exxonmobil Upstream Research Company | Method for using dynamic target region for well path/drill center optimization |
CA2786584C (en) | 2010-03-12 | 2017-07-18 | Exxonmobil Upstream Research Company | Dynamic grouping of domain objects via smart groups |
US9595129B2 (en) | 2012-05-08 | 2017-03-14 | Exxonmobil Upstream Research Company | Canvas control for 3D data volume processing |
US10584570B2 (en) | 2013-06-10 | 2020-03-10 | Exxonmobil Upstream Research Company | Interactively planning a well site |
CN104481523B (zh) * | 2014-11-11 | 2018-07-13 | 中国石油天然气股份有限公司 | 稠油溶解气驱开发模拟实验系统和方法 |
FR3034894B1 (fr) * | 2015-04-09 | 2018-08-10 | Services Petroliers Schlumberger | Reservoir de saturation de gisement de petrole et modelisation de permeabilite |
FR3036820B1 (fr) | 2015-06-01 | 2021-12-31 | Services Petroliers Schlumberger | Modelisation de la saturation et permeabilite de reservoir de champ petrolifere |
FR3038408B1 (fr) | 2015-06-30 | 2017-08-25 | Services Petroliers Schlumrberger | Modelisation de la saturation et permeabilite de reservoir de champ petrolifere |
CN105443080A (zh) * | 2015-12-04 | 2016-03-30 | 中国石油天然气股份有限公司 | 边底水含气油藏衰竭开发的模拟装置 |
CN107944599B (zh) * | 2017-10-31 | 2020-10-09 | 中国石油天然气股份有限公司 | 油气水平井产量的预测方法 |
CN109281662B (zh) * | 2018-08-16 | 2022-03-29 | 中国石油天然气股份有限公司 | 水平井不同节段位置储层相对可压性指数的确定方法 |
US11492895B2 (en) | 2018-11-13 | 2022-11-08 | Saudi Arabian Oil Company | Relative permeability ratio from wellbore drilling data |
CN110671084B (zh) * | 2019-09-04 | 2021-07-06 | 中国石油化工股份有限公司 | 普通稠油水溶性降粘剂驱数值模拟的方法 |
CN110863806B (zh) * | 2019-11-28 | 2021-07-23 | 西安石油大学 | 一种二氧化碳驱气体前缘动态变化预测方法 |
US20240218769A1 (en) * | 2022-12-30 | 2024-07-04 | Landmark Graphics Corporation | Adjusting fluid injection into a wellbore based on relative permeability of a formation |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6052520A (en) * | 1998-02-10 | 2000-04-18 | Exxon Production Research Company | Process for predicting behavior of a subterranean formation |
US6227296B1 (en) * | 1998-11-03 | 2001-05-08 | Exxonmobil Upstream Research Company | Method to reduce water saturation in near-well region |
FR2837572B1 (fr) * | 2002-03-20 | 2004-05-28 | Inst Francais Du Petrole | Methode pour modeliser la production d'hydrocarbures par un gisement souterrain soumis a une depletion |
US6810332B2 (en) * | 2003-01-31 | 2004-10-26 | Chevron U.S.A. Inc. | Method for computing complexity, confidence and technical maturity indices for reservoir evaluations |
CN1302386C (zh) * | 2003-10-17 | 2007-02-28 | 大庆油田有限责任公司 | 低浓度表面活性剂与相态结合的三元复合驱计算机仿真方法 |
-
2006
- 2006-12-20 CN CN2006800525037A patent/CN101366041B/zh not_active Expired - Fee Related
- 2006-12-20 WO PCT/US2006/049095 patent/WO2007076044A2/en active Application Filing
- 2006-12-20 CA CA2634757A patent/CA2634757C/en active Active
- 2006-12-20 BR BRPI0620170-9A patent/BRPI0620170A2/pt not_active IP Right Cessation
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102576374A (zh) * | 2009-09-02 | 2012-07-11 | 界标制图有限公司 | 用于油气形成建模的系统和方法 |
CN102576374B (zh) * | 2009-09-02 | 2015-04-08 | 界标制图有限公司 | 用于油气形成建模的系统和方法 |
US9390207B2 (en) | 2009-09-02 | 2016-07-12 | Landmark Graphics Corporation | System and method of hydrocarbon formation modeling |
CN102713683A (zh) * | 2010-02-26 | 2012-10-03 | 雪佛龙美国公司 | 用于利用多点统计模拟来对储层特性走向建模的方法和系统 |
CN102713683B (zh) * | 2010-02-26 | 2016-04-06 | 雪佛龙美国公司 | 用于利用多点统计模拟来对储层特性走向建模的方法和系统 |
CN103590807A (zh) * | 2012-08-13 | 2014-02-19 | 中国石油天然气集团公司 | 一种确定稠油热采蒸汽腔形态的方法 |
CN103590807B (zh) * | 2012-08-13 | 2016-06-08 | 中国石油天然气集团公司 | 一种确定稠油热采蒸汽腔形态的方法 |
CN104481504A (zh) * | 2014-11-11 | 2015-04-01 | 中国石油天然气股份有限公司 | 稠油溶解气驱开发模拟实验系统及其方法 |
CN104481504B (zh) * | 2014-11-11 | 2017-09-01 | 中国石油天然气股份有限公司 | 稠油溶解气驱开发模拟实验系统及其方法 |
CN107101773A (zh) * | 2017-04-14 | 2017-08-29 | 中国石油天然气股份有限公司 | 泡点压力值和拟泡点压力值的测试方法、及测试装置 |
CN107101773B (zh) * | 2017-04-14 | 2019-07-09 | 中国石油天然气股份有限公司 | 泡点压力值和拟泡点压力值的测试方法、及测试装置 |
Also Published As
Publication number | Publication date |
---|---|
CN101366041B (zh) | 2011-10-12 |
CA2634757A1 (en) | 2007-07-05 |
WO2007076044A2 (en) | 2007-07-05 |
WO2007076044A3 (en) | 2007-12-27 |
BRPI0620170A2 (pt) | 2011-11-01 |
CA2634757C (en) | 2014-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101366041B (zh) | 利用溶气驱重油进行储层模拟的方法 | |
US9020793B2 (en) | Method, system and program storage device for reservoir simulation utilizing heavy oil solution gas drive | |
Sheng et al. | Modelling foamy oil flow in porous media | |
Behmanesh et al. | Treatment of rate-transient analysis during boundary-dominated flow | |
CN107563899A (zh) | 油气井产能预测方法及装置 | |
CA2888205A1 (en) | Predicting performance of gas condensate reservoirs | |
US20100185424A1 (en) | Method, Program and Computer System for Conciliating Hydrocarbon Reservoir Model Data | |
Nagarajan et al. | The Critical Role of Rock and Fluid in Unconventional Shale Reservoir Performance | |
Maini | Is it futile to measure relative permeability for heavy oil reservoirs? | |
Lenormand et al. | Comparison of four numerical simulators for SCAL experiments | |
Lala | A novel model for reservoir rock tortuosity estimation | |
US20150226061A1 (en) | System and method for estimating flow capacity of a reservoir | |
Zhang et al. | Using digital rock modeling to estimate permeability and capillary pressure from NMR and geochemical logs | |
US20190196060A1 (en) | Method for modeling a sedimentary basin | |
Ali | Developments in measurement and interpretation techniques in coreflood tests to determine relative permeabilities | |
US20190196059A1 (en) | Method for modeling a sedimentary basin | |
Lucia | Petrophysical rock properties | |
Elsharkawy | MB solution for high pressure gas reservoirs | |
Kosova et al. | Uncertainty in Oil Reservoir Properties Deterministic and Stochastic Methods of Reserves Estimation | |
Taghavi et al. | The Impact of Autonomous Inflow Control Valve on Improved Oil Recovery in a Thin-Oil-Rim Reservoir | |
RU2794707C1 (ru) | Способ прогнозирования объемов добычи углеводородов из месторождений нефти и газа с использованием компьютерного моделирования | |
Li | Optimization of multistage hydraulic fracturing treatment for maximization of the tight gas productivity | |
Bejan et al. | Methane hydrates in porous layers: Gas formation and convection | |
RU2819121C1 (ru) | Способ определения гидродинамической связи между участками продуктивного пласта и фильтрационно-емкостных свойств межскважинного пространства сеноманской залежи при запуске промысла после остановок по результатам интегрального гидропрослушивания на скважинах | |
Savioli et al. | Modeling multiphase fluid flow in unconventional reservoirs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20111012 Termination date: 20141220 |
|
EXPY | Termination of patent right or utility model |