CN108252688A - 致密油藏二氧化碳吞吐影响因素敏感性分析方法及其应用 - Google Patents
致密油藏二氧化碳吞吐影响因素敏感性分析方法及其应用 Download PDFInfo
- Publication number
- CN108252688A CN108252688A CN201810017037.3A CN201810017037A CN108252688A CN 108252688 A CN108252688 A CN 108252688A CN 201810017037 A CN201810017037 A CN 201810017037A CN 108252688 A CN108252688 A CN 108252688A
- Authority
- CN
- China
- Prior art keywords
- carbon dioxide
- reservoir
- oil reservoir
- tight
- influence
- 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.)
- Pending
Links
- 238000010206 sensitivity analysis Methods 0.000 title claims abstract description 21
- 230000000638 stimulation Effects 0.000 title abstract description 8
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims abstract description 292
- 229910002092 carbon dioxide Inorganic materials 0.000 claims abstract description 146
- 239000001569 carbon dioxide Substances 0.000 claims abstract description 146
- 230000000694 effects Effects 0.000 claims abstract description 111
- 238000000034 method Methods 0.000 claims abstract description 60
- 238000009792 diffusion process Methods 0.000 claims abstract description 51
- 238000011084 recovery Methods 0.000 claims abstract description 47
- 230000035699 permeability Effects 0.000 claims abstract description 38
- 239000011159 matrix material Substances 0.000 claims abstract description 36
- 230000008569 process Effects 0.000 claims abstract description 31
- 239000012071 phase Substances 0.000 claims description 54
- 239000011148 porous material Substances 0.000 claims description 26
- 238000004364 calculation method Methods 0.000 claims description 20
- 230000008859 change Effects 0.000 claims description 18
- 239000007791 liquid phase Substances 0.000 claims description 18
- 239000012530 fluid Substances 0.000 claims description 15
- 230000035945 sensitivity Effects 0.000 claims description 14
- 238000009826 distribution Methods 0.000 claims description 11
- 238000013178 mathematical model Methods 0.000 claims description 9
- 238000002474 experimental method Methods 0.000 claims description 6
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 3
- 238000005094 computer simulation Methods 0.000 claims 1
- 238000013461 design Methods 0.000 abstract description 2
- 239000003921 oil Substances 0.000 description 86
- 239000007789 gas Substances 0.000 description 25
- 238000004458 analytical method Methods 0.000 description 12
- 239000000243 solution Substances 0.000 description 11
- 238000002347 injection Methods 0.000 description 10
- 239000007924 injection Substances 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 238000006073 displacement reaction Methods 0.000 description 8
- 238000011161 development Methods 0.000 description 7
- 238000004519 manufacturing process Methods 0.000 description 7
- 238000007906 compression Methods 0.000 description 6
- 239000000203 mixture Substances 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 5
- 230000006835 compression Effects 0.000 description 5
- 239000003129 oil well Substances 0.000 description 4
- 230000009897 systematic effect Effects 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000010779 crude oil Substances 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 230000009747 swallowing Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000005431 greenhouse gas Substances 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
- E21B43/164—Injecting CO2 or carbonated water
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Physics & Mathematics (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)
Abstract
本发明属于致密油藏技术领域，涉及致密油藏二氧化碳吞吐影响因素敏感性分析方法及其应用。本发明提供的致密油藏二氧化碳吞吐影响因素敏感性分析方法，包括如下步骤：(a)建立油藏模型；(b)建立致密油藏相态理论模型；(c)利用油藏组分模型模拟器，对影响二氧化碳吞吐过程的因素包括基质渗透率、毛细管力效应和二氧化碳的扩散效应进行敏感性分析；(d)根据敏感性分析的结果，明确提高油藏采收率的措施或方向。本发明系统性的分析了基质渗透率、毛细管力效应、二氧化碳的扩散效应等因素对致密油藏的二氧化碳吞吐过程的影响，为致密油藏二氧化碳吞吐过程的设计提供更为有力的理论依据和指导。
Description
技术领域
本发明属于致密油藏领域，尤其涉及一种致密油藏二氧化碳吞吐 影响因素敏感性分析方法及其应用。
背景技术
随着世界石油供需矛盾的日益突出和优质石油资源的逐渐匮乏， 以及常规油藏的大量开采，非常规油气资源例如致密油气已经成为重 要的油气资源接替区。近年来，非常规油气藏的勘探开发日益增多。 由于非常规油气藏具有低孔、低渗，储集层中纳米级孔隙发育的特点， 常采用水平井多段压裂等手段进行开发。根据已有实验和数值模拟研 究显示，注二氧化碳驱是提高致密油藏最终采收率的有效方法之一。 在二氧化碳注入的过程中，二氧化碳与原油接触，从油相中抽提轻质 组分及中间组分，从而降低原油粘度及界面张力。当界面张力为零时， 二氧化碳与地层流体达到混相，从而实现二氧化碳的混相驱替过程。 二氧化碳驱不但有效解决了工业生产、人民生活温室气体排放和储存 等重大问题，而且由于其独特有效的驱替机制，使其成为气驱提高油 藏采收率的一项关键技术。
近年来，已有众多国内外科研工作者对致密油藏的二氧化碳驱油 过程进行研究。例如，Yu W等研究发现，二氧化碳的扩散效应是影 响二氧化碳驱油效果的重要因素之一；Alharthy N等基于室内实验 和数值模拟，分析了Bakken致密油藏的提高采收率过程；Zuloaga－Molero P等分析了天然裂缝的存在的复杂裂缝网络对二氧 化碳驱油过程的影响；Kim T H等研究了应力敏感效应对二氧化碳吞 吐采油过程的影响。但上述研究均未对影响二氧化碳吞吐效果的因素 进行系统的敏感性分析。
鉴于此，特提出本发明。
发明内容
本发明的第一目的在于提供一种致密油藏二氧化碳吞吐影响因 素敏感性分析方法，以解决上述背景技术中提到的缺乏系统性的问 题。
本发明的第二目的在于提供一种包含上述致密油藏二氧化碳吞 吐影响因素敏感性分析方法在提高油藏采收率中的应用。
为实现上述目的，本发明采用的技术方案为：
根据本发明的一个方面，本发明提供一种致密油藏二氧化碳吞吐 影响因素敏感性分析方法，包括如下步骤：
(a)建立油藏模型；
(b)建立致密油藏相态理论模型；
(c)利用油藏组分模型模拟器，对影响二氧化碳吞吐过程的因 素包括基质渗透率、毛细管力效应和二氧化碳的扩散效应进行敏感性 分析；
(d)根据敏感性分析的结果，明确提高油藏采收率的措施或方 向。
作为进一步优选技术方案，步骤(a)中，统计分析致密油藏的 典型特征，建立相应的油藏三维和/或平面模型，并确定油藏及裂缝 的相关参数。
作为进一步优选技术方案，步骤(a)中，所述的油藏为Bakken 致密油藏。
作为进一步优选技术方案，步骤(b)中，为考虑致密油藏流体 各组分临界性质的变化，构建数学模型，并利用蒙特卡罗模拟方法对 临界温度和临界压力进行计算，所述的数学模型为：
其中，rp为孔喉半径，ΔTc *和ΔPc *分别表示毛细管力效应影响的临 界温度和临界压力的相对变化量，σLJ为Lennard－Jones参数。
作为进一步优选技术方案，步骤(b)中，考虑毛细管力效应的 相平衡计算模型为：
当油相和气相中各组分的逸度相等时，体系达到相平衡，即：
根据质量守恒定律，可得：
其中，fi为液相、气相中组分i的逸度；PV和PL分别为气相和液 相的压力；zi为体系中组分i的总摩尔分数；xi和yi分别为组分i在 液相和气相中的摩尔分数；F为总摩尔数；L和V分别为液相和气相 的摩尔数；Nc为体系中的总组分数；
采用Young－Laplace方程来计算毛细管力Pcap：
PV-PL＝Pcap (7)
其中，θ为接触角，σ为界面张力。
作为进一步优选技术方案，步骤(c)中，考虑毛细管力效应的 影响时，设定一系列孔隙半径值，利用所述的数学模型和相平衡计算 模型，考察致密油藏流体各组分临界性质的变化，并根据实验测得的 实际孔隙分布数据，预测致密油藏的最终采收率；
通过对比考虑毛细管力效应和忽略毛细管力效应时的油藏最终 采收率，明确毛细管力效应对二氧化碳吞吐的影响程度。
作为进一步优选技术方案，步骤(c)中，考虑基质渗透率的影 响时，保持其他参数不变，设定一系列基质渗透率值，分析油藏基质 渗透率的变化对油藏最终采收率的影响，明确基质渗透率对二氧化碳 吞吐的影响程度。
作为进一步优选技术方案，步骤(c)中，考虑二氧化碳的扩散 效应的影响时，通过对比考虑二氧化碳的扩散效应和忽略二氧化碳的 扩散效应时的二氧化碳的分布情况以及对油藏最终采收率的影响，明 确二氧化碳的扩散效应对二氧化碳吞吐的影响程度。
作为进一步优选技术方案，步骤(d)中，根据敏感性分析的结 果，确定考虑基质渗透率、毛细管力效应和/或二氧化碳的扩散效应 时对二氧化碳吞吐过程的影响程度，以明确提高油藏采收率的措施或 方向。
根据本发明的另一个方面，本发明还提供一种包含上述的致密油 藏二氧化碳吞吐影响因素敏感性分析方法在提高油藏采收率中的应 用。
与现有技术相比，本发明的有益效果在于：
1、本发明提供的致密油藏二氧化碳吞吐影响因素敏感性分析方 法，系统性的分析了基质渗透率、毛细管力效应、二氧化碳的扩散效 应等因素对致密油藏的二氧化碳吞吐过程的影响，为致密油藏二氧化 碳吞吐过程的设计提供更为有力的理论依据和指导。
2、通过本发明的分析方法可知，二氧化碳的扩散效应是影响二 氧化碳吞吐过程的重要参数之一，考虑二氧化碳的扩散效应时，最终 采收率提高4.3％；较低的渗透率及纳米孔隙中的毛细管力效应提高 了最终采收率，对致密油藏的二氧化碳吞吐过程是有益的；考虑注二 氧化碳过程中的毛细管力效应时，二氧化碳的扩散效应更加显著。
3、本发明明确了影响致密油藏二氧化碳驱开发效果的主控因素， 提高判断影响因素的准确性，不仅能够为提高致密油藏油井产量提供 有力的技术支持，还可以为致密油藏合理开发提供了理论依据及参 考，具有重大的理论和实际应用意义。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方 案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简 单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式， 对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可 以根据这些附图获得其他的附图。
图1为本发明一种实施方式提供的基准算例的平面示意图；
图2为本发明一种实施方式提供的不同孔隙半径下Bakken致密 油藏注二氧化碳驱油过程中最小混相压力的变化示意图；
图3为本发明一种实施方式提供的基质渗透率对二氧化碳吞吐效 果的影响示意图；
图4为本发明一种实施方式提供的二氧化碳的扩散效应对最终采 收率的影响示意图；
图5为本发明一种实施方式提供的忽略二氧化碳扩散效应时二氧 化碳分布示意图(图中数值表示该网格处二氧化碳的浓度)；
图6为本发明一种实施方式提供的考虑二氧化碳扩散效应时二氧 化碳分布示意图(图中数值表示该网格处二氧化碳的浓度)；
图7为本发明一种实施方式提供的Bakken油藏流体各组分临界 性质的变化示意图；图7(a)和图7(b)分别为临界温度和临界压 力的变化示意图；
图8为本发明一种实施方式提供的毛细管力效应对二氧化碳吞吐 效果的影响示意图；
图9为本发明一种实施方式提供的二氧化碳吞吐过程对最终采收 率的影响的实例分析。
具体实施方式
下面将结合实施方式和实施例对本发明的实施方案进行详细描 述，但是本领域技术人员将会理解，下列实施方式和实施例仅用于说 明本发明，而不应视为限制本发明的范围。基于本发明中的实施例， 本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其 他实施例，都属于本发明保护的范围。实施例中未注明具体条件者， 按照常规条件或制造商建议的条件进行。
第一方面，在至少一个实施例中提供一种致密油藏二氧化碳吞吐 影响因素敏感性分析方法，包括如下步骤：
(a)建立油藏模型；
(b)建立致密油藏相态理论模型；
(c)利用油藏组分模型模拟器，对影响二氧化碳吞吐过程的因 素包括基质渗透率、毛细管力效应和二氧化碳的扩散效应进行敏感性 分析；
(d)根据敏感性分析的结果，明确提高油藏采收率的措施或方 向。
本发明基于热力学理论，建立了考察致密油藏储层毛细管力效应 的相态平衡模型，并与油藏组分模型模拟器相结合，考察了基质渗透 率、毛细管力效应和二氧化碳的扩散效应等因素对致密油藏二氧化碳 吞吐过程的影响，为提高油藏采收率提供了方向和理论依据。
可以理解的是，所述的油藏组分模型模拟器为本领域技术人员公 知的模拟软件。
这里需要说明的是，根据本发明所建立的模型以及本发明的分析 方法，还可以对致密油藏二氧化碳吞吐过程的其他影响因素进行分 析，不仅限于基质渗透率、毛细管力效应和二氧化碳的扩散效应的敏 感性分析，本发明只是以上述三种影响因素为例进行了进一步的详细 说明。
需要说明的是，尽管上述步骤是以(a)、(b)、(c)和(d)顺序 描述了本发明的分析方法，但是，这并非要求或者暗示必须按照该特 定顺序来执行这些操作。步骤(a)、(b)、(c)和(d)的分析可以同 时或者任意先后执行。例如，可以同时建立油藏模型和致密油藏相态理论模型，也可以先建立油藏模型后建立致密油藏相态理论模型，还 可以先建立致密油藏相态理论模型后建立油藏模型。
[关于油藏模型]
建立油藏模型时，可以先统计分析致密油藏的典型特征，建立相 应的油藏三维和/或平面模型，并确定油藏及裂缝的相关参数。
优选地，以典型的Bakken致密油藏为例进行进一步的分析。
在一种优选的实施方式中，建立油藏模型的基准算例的平面示意 图，如图1所示，其长、宽和高分别为104m、384m和12m，包含 4条人工裂缝；裂缝半长为64m，裂缝间距为24m。
油藏及裂缝的主要性质及相关参数如表1所示；其中，油藏的初 始压力为51.7MPa，井底流压为10.4MPa。
表1 Bakken致密油藏及裂缝的相关参数
| 参数 | 值 |
| 油藏初始压力 | 51.7MPa |
| 油藏温度 | 116℃ |
| 生产时间 | 10年 |
| 基质孔隙度 | 0.056 |
| 基质渗透率 | 0.01md |
| 初始水饱和度 | 0.49 |
| 岩石压缩系数 | 1.5×10－15MPa－1 |
| 裂缝半长 | 64m |
| 裂缝间距 | 24m |
数值模拟中考虑二氧化碳的扩散效应，并将扩散系数定为 0.0005cm2/s。生产过程中，水平井首先作为生产井生产，后转为注 入井，以一定的注入速率(如2.8×104m3/day的注入速率)注入二 氧化碳，停止注入后，焖井一段时间(如一个月)开井生产，至此， 一次二氧化碳吞吐结束。油井生产一年后，进行新一轮的吞吐作业； 该井可经历几次(如两次、三次、四次或五次)吞吐过程。通过所建 立的理论计算模型与油藏组分模型模拟器相结合，对基质渗透率、毛 细管力效应等敏感性因素进行系统性分析。
[关于致密油藏相态理论模型]
1)纳米孔隙中的高毛细管力效应会显著影响流体的临界性质和 相态平衡。当孔隙半径小于10nm时，由于壁表面及分子间的相互作 用，组分的临界性质变化较大。为考虑致密油藏流体各组分临界性质 的变化，构建数学模型，并利用蒙特卡罗模拟方法对临界温度和临界 压力进行计算，所述的数学模型为：
其中，rp为孔喉半径，和分别表示毛细管力效应影响的临 界温度和临界压力的相对变化量，σLJ为Lennard－Jones参数。
这里需要说明的是，在目前的研究中，常将rp近似为孔隙半径。
2)考虑毛细管力效应的相平衡计算模型为：
当油相和气相中各组分的逸度相等时，体系达到相平衡，即：
根据质量守恒定律，可得：
其中，fi为液相、气相中组分i的逸度；PV和PL分别为气相和液 相的压力；zi为体系中组分i的总摩尔分数；xi和yi分别为组分i在 液相和气相中的摩尔分数；F为总摩尔数；L和V分别为液相和气相 的摩尔数；Nc为体系中的总组分数。
采用Young－Laplace方程来计算毛细管力Pcap：
PV-PL＝Pcap (7)
其中，θ为接触角，σ为界面张力。
界面张力σ的计算可利用式(9)：
其中，和分别为液相、气相的密度。[P]i为组分i的等张 比容。纯组分的性质可通过查阅化工手册得到。
逸度系数可通过Peng－Robinson状态方程求得：
其中，Vm为组分i的摩尔体积，R为普适气体常数。参数a和b 利用van der Waals混合定律来计算。由于毛细管力的存在，气、 液两相的压力不相等，因此，利用式(10)对气、液相的压缩因子分别 求解。
对于液相，式(10)可转化为：
(ZL)3-(1-BL)(ZL)2+(AL-2BL-3(BL)2)ZL-(ALBL-(BL)2-(BL)3)＝0 (11)
对于气相，式(10)可转化为：
(ZV)3-(1-BV)(ZV)2+(AV-2BV-3(BV)2)ZV-(AVBV-(BV)2-(BV)3)＝0 (12)
其中，ZL和ZV分别 为液相和气相的压缩因子。
为确定液相和气相的压缩因子，采用Gibbs自由能最小化原理确 定式(11)和式(12)的解。即：在P＝PL条件下求解式(11)，方程三个 解中的最小值为液相的压缩因子；在P＝PV条件下求解式(12)，三个 解的最大值为气相的压缩因子。采用Newton－Raphson迭代法求解非 线性方程组，进而求得组分i在气、液相中的逸度系数。
应当理解的是，上述公式中未特别说明的参数或计算方式按照本 领域中常规的理解、应用、计算或查阅即可。
3)致密油藏相态理论模型的验证
当界面张力为零时，此时对应的压力即为最小混相压力(Minimum MiscibilityPressure，简称MMP)。利用该原理，结合上述理论模 型计算公式计算了Bakken致密油藏注二氧化碳驱油过程中最小混相 压力的变化。其中，孔隙半径为20nm，油藏温度为116℃。同时为 进一步探究毛细管力对最小混相压力的影响，分别计算了孔隙半径为 10nm和5nm时，最小混相压力的变化。计算结果如图2所示。
图2的计算结果表明，二氧化碳与Bakken致密油藏流体的最小 混相压力为23.1MPa，与Kurtoglu等利用升泡仪法测得的22.8MPa 相对误差仅为2％，计算结果符合较好，进而验证了采用上述致密油 藏相态理论模型算法的准确性。
此外，当孔隙半径为10nm时，最小混相压力为21.9MPa；当 孔隙半径降至5nm时，最小混相压力为20MPa。因此，孔隙半径的 减小会使体系的最小混相压力降低，即二氧化碳与油藏流体可以在较 低的压力下达到混相，有利于提高二氧化碳的驱油效果。
[基质渗透率的影响]
考虑基质渗透率的影响时，保持其他参数不变，设定一系列基质 渗透率值，分析油藏基质渗透率的变化对油藏最终采收率的影响，明 确基质渗透率对二氧化碳吞吐的影响程度。
在一种优选的实施方式中，一系列基质渗透率值可以为0.001 md、0.01md和0.1md，保持其他参数不变，在上述油藏模型、理论 计算模型以及油藏组分模型模拟器的结合下，进行模拟计算，分析油 藏基质渗透率的变化对油藏最终采收率的影响。
图3显示了基质渗透率对二氧化碳吞吐效果的影响，由图3可以 看出，当基质渗透率分别为0.001md、0.01md和0.1md时，采用 二氧化碳吞吐方式，会使油藏的最终采收率分别提高7.5％、6.4％和 5.9％。即基质渗透率越低，二氧化碳的吞吐效果越明显。这是因为一 次采油结束时，渗透率较低的储层具有较大的残油饱和度，与渗透率 较高的储层相比，二氧化碳的扩散效应比对流效应更加显著。因此， 较低的渗透率对于二氧化碳驱油过程是非常有益的。
[二氧化碳扩散效应的影响]
考虑二氧化碳的扩散效应的影响时，通过对比考虑二氧化碳的扩 散效应和忽略二氧化碳的扩散效应时的二氧化碳的分布情况以及对 油藏最终采收率的影响，明确二氧化碳的扩散效应对二氧化碳吞吐的 影响程度。
在一种优选的实施方式中，在上述油藏模型、理论计算模型以及 油藏组分模型模拟器的结合下，进行模拟计算，分析考虑二氧化碳扩 散效应和忽略二氧化碳的扩散效应时对油藏最终采收率的影响。
图4显示了二氧化碳的扩散效应对最终采收率的影响。由图4 可以看出，二氧化碳的扩散效应是影响二氧化碳吞吐效果的一个重要 因素，考虑二氧化碳的扩散效应时，油藏生产10年时的最终采收率 比未考虑扩散效应的最终采收率高出4.3％。
进一步地，图5和图6分别显示了忽略二氧化碳扩散效应时二氧 化碳的分布示意图和考虑二氧化碳扩散效应时二氧化碳的分布示意 图。由图5和图6的对比分析可以看出，未考虑二氧化碳的扩散效应 时，由于基质渗透率较低，大量的二氧化碳分子主要集中在人工裂缝 周围；而考虑二氧化碳的扩散效应时，更多的二氧化碳分子从裂缝运 移到基质中，与地层流体接触并逐渐达到混相，从而有效提高了油藏 的最终采收率。因此，研究二氧化碳的吞吐过程时，不应忽略二氧化 碳扩散效应的影响。
[纳米孔隙中毛细管力效应的影响]
考虑毛细管力效应的影响时，设定一系列孔隙半径值，利用所述 的数学模型和相平衡计算模型，考察致密油藏流体各组分临界性质的 变化，并根据实验测得的实际孔隙分布数据，预测致密油藏的最终采 收率；通过对比考虑毛细管力效应和忽略毛细管力效应时的油藏最终 采收率，明确毛细管力效应对二氧化碳吞吐的影响程度。
在一种优选的实施方式中，一系列孔隙半径值可以为20nm、10 nm和5nm，在上述油藏模型、理论计算模型以及油藏组分模型模拟 器的结合下，进行模拟计算，分析毛细管力效应对油藏最终采收率的 影响。
由于低渗透储层中纳米级孔隙发育，毛细管力效应显著，因此流 体的相态平衡和流动过程较常规储层差别较大。图7显示了孔隙半径 值为20nm、10nm和5nm时油藏流体各组分临界性质的变化；图7 (a)和图7(b)分别为临界温度和临界压力的变化。从图7可以看出，随着孔隙半径的减小，组分的临界温度和临界压力均降低，重组 分性质的变化更显著。
根据实验测得的实际孔隙分布数据，将孔隙半径划分为五个区 域：小于10nm(27％)、10－20nm(26％)、20－30nm(30％)、30－50nm(13％) 及大于50nm(4％)的区域。不同区域的流体具有不同的PVT性质， 并将其随机分布于油藏模型中，从而能够更加合理的预测致密油藏的 最终采收率。
进一步地，为探究毛细管力效应对二氧化碳吞吐效果的影响，对 比了油藏生产10年后的最终采收率，如图8所示。从图8可以看出， 忽略毛细管力效应对相平衡的影响时，采用二氧化碳吞吐会使最终采 收率提高7.1％，而考虑毛细管力效应时，最终采收率继续提高2.5％。 考虑毛细管力效应时，二氧化碳的扩散更加显著。二氧化碳从裂缝逐 渐运移到基质中，从而使基质中二氧化碳的浓度升高，有利于二氧化 碳的驱油过程。因此，毛细管力效应可提高二氧化碳吞吐的效果。研 究致密油藏的二氧化碳吞吐过程时，不应忽略毛细管力效应的影响。
第二方面，本实施方式还提供一种包含上述的致密油藏二氧化碳 吞吐影响因素敏感性分析方法在提高油藏采收率中的应用。
通过本发明提供的理论模型以及系统性的分析，明确了影响致密 油藏二氧化碳驱开发效果的主控因素，提高判断影响因素的准确性， 不仅能够为提高致密油藏油井产量提供有力的技术支持，还可以为致 密油藏合理开发提供了理论依据及参考，具有重大的理论和实际应用 意义。
应用例
对Bakken致密油藏二氧化碳的吞吐过程进行研究，同时考虑二 氧化碳的扩散效应及毛细管力对相平衡的影响，分析吞吐过程提高采 收率的效果。油藏模型的尺寸为3219m×805m×15m，平面上网格数 为264×66。油藏中两口水平井进行生产，每口水平井伴有30条人工 裂缝。油藏及裂缝的参数如表1所示，在上述油藏模型、本发明提供 的理论模型以及油藏组分模型模拟器的结合下，对比有无二氧化碳注 入的油藏最终采收率的变化，结果如图9所示。从图9可以看出，考 虑二氧化碳的扩散效应及毛细管力对相平衡的影响时，采用二氧化碳 吞吐过程可使最终采收率提高4.8％，说明二氧化碳吞吐能够提高致 密储层的动用程度，是提高致密油藏油井产量的有效手段之一。
综上，通过本发明的分析可以得知，(1)二氧化碳的扩散效应是 影响二氧化碳吞吐过程的重要参数之一，考虑二氧化碳的扩散效应 时，最终采收率提高4.3％；(2)较低的渗透率及纳米孔隙中的毛细 管力效应提高了最终采收率，对致密油藏的二氧化碳吞吐过程是有益 的；(3)考虑注二氧化碳过程中的毛细管力效应时，二氧化碳的扩散 效应更加显著。实例分析表明，Bakken致密油藏二氧化碳吞吐，最 终采收率提高4.8％。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案， 而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明， 本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载 的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替 换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各 实施例技术方案的范围。
Claims (10)
1.一种致密油藏二氧化碳吞吐影响因素敏感性分析方法，其特征在于，包括如下步骤：
(a)建立油藏模型；
(b)建立致密油藏相态理论模型；
(c)利用油藏组分模型模拟器，对影响二氧化碳吞吐过程的因素包括基质渗透率、毛细管力效应和二氧化碳的扩散效应进行敏感性分析；
(d)根据敏感性分析的结果，明确提高油藏采收率的措施或方向。
2.根据权利要求1所述的致密油藏二氧化碳吞吐影响因素敏感性分析方法，其特征在于，步骤(a)中，统计分析致密油藏的典型特征，建立相应的油藏三维和/或平面模型，并确定油藏及裂缝的相关参数。
3.根据权利要求2所述的致密油藏二氧化碳吞吐影响因素敏感性分析方法，其特征在于，步骤(a)中，所述的油藏为Bakken致密油藏。
4.根据权利要求1所述的致密油藏二氧化碳吞吐影响因素敏感性分析方法，其特征在于，步骤(b)中，为考虑致密油藏流体各组分临界性质的变化，构建数学模型，并利用蒙特卡罗模拟方法对临界温度和临界压力进行计算，所述的数学模型为：
其中，rp为孔喉半径，和分别表示毛细管力效应影响的临界温度和临界压力的相对变化量，σLJ为Lennard－Jones参数。
5.根据权利要求4所述的致密油藏二氧化碳吞吐影响因素敏感性分析方法，其特征在于，步骤(b)中，考虑毛细管力效应的相平衡计算模型为：
当油相和气相中各组分的逸度相等时，体系达到相平衡，即：
根据质量守恒定律，可得：
Fzi＝xiL+yiV,i＝1,...,Nc， (6)
其中，fi为液相、气相中组分i的逸度；PV和PL分别为气相和液相的压力；zi为体系中组分i的总摩尔分数；xi和yi分别为组分i在液相和气相中的摩尔分数；F为总摩尔数；L和V分别为液相和气相的摩尔数；Nc为体系中的总组分数；
采用Young－Laplace方程来计算毛细管力Pcap：
PV-PL＝Pcap (7)
其中，θ为接触角，σ为界面张力。
6.根据权利要求5所述的致密油藏二氧化碳吞吐影响因素敏感性分析方法，其特征在于，步骤(c)中，考虑毛细管力效应的影响时，设定一系列孔隙半径值，利用所述的数学模型和相平衡计算模型，考察致密油藏流体各组分临界性质的变化，并根据实验测得的实际孔隙分布数据，预测致密油藏的最终采收率；
通过对比考虑毛细管力效应和忽略毛细管力效应时的油藏最终采收率，明确毛细管力效应对二氧化碳吞吐的影响程度。
7.根据权利要求5所述的致密油藏二氧化碳吞吐影响因素敏感性分析方法，其特征在于，步骤(c)中，考虑基质渗透率的影响时，保持其他参数不变，设定一系列基质渗透率值，分析油藏基质渗透率的变化对油藏最终采收率的影响，明确基质渗透率对二氧化碳吞吐的影响程度。
8.根据权利要求5所述的致密油藏二氧化碳吞吐影响因素敏感性分析方法，其特征在于，步骤(c)中，考虑二氧化碳的扩散效应的影响时，通过对比考虑二氧化碳的扩散效应和忽略二氧化碳的扩散效应时的二氧化碳的分布情况以及对油藏最终采收率的影响，明确二氧化碳的扩散效应对二氧化碳吞吐的影响程度。
9.根据权利要求1～8任一项所述的致密油藏二氧化碳吞吐影响因素敏感性分析方法，其特征在于，步骤(d)中，根据敏感性分析的结果，确定考虑基质渗透率、毛细管力效应和/或二氧化碳的扩散效应时对二氧化碳吞吐过程的影响程度，以明确提高油藏采收率的措施或方向。
10.包含权利要求1～9任一项所述的致密油藏二氧化碳吞吐影响因素敏感性分析方法在提高油藏采收率中的应用。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810017037.3A CN108252688A (zh) | 2018-01-08 | 2018-01-08 | 致密油藏二氧化碳吞吐影响因素敏感性分析方法及其应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810017037.3A CN108252688A (zh) | 2018-01-08 | 2018-01-08 | 致密油藏二氧化碳吞吐影响因素敏感性分析方法及其应用 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108252688A true CN108252688A (zh) | 2018-07-06 |
Family
ID=62726094
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810017037.3A Pending CN108252688A (zh) | 2018-01-08 | 2018-01-08 | 致密油藏二氧化碳吞吐影响因素敏感性分析方法及其应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108252688A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113241127A (zh) * | 2021-03-19 | 2021-08-10 | 中国石油大学(北京) | 相平衡模型构建方法、装置、设备及存储介质 |
| CN113592194A (zh) * | 2021-08-23 | 2021-11-02 | 北京科技大学 | 建立co2吞吐效果预测模型的方法以及co2吞吐效果评价方法 |
| CN113919448A (zh) * | 2021-12-14 | 2022-01-11 | 武汉大学 | 一种任意时空位置二氧化碳浓度预测影响因素分析方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4217955A (en) * | 1978-05-15 | 1980-08-19 | Petroleum Recovery Institute | Oil recovery process |
| CN107130960A (zh) * | 2017-06-26 | 2017-09-05 | 中海石油(中国)有限公司 | 一种考虑致密砂岩油藏非线性渗流特征的数值模拟方法 |
| CN107130949A (zh) * | 2016-02-25 | 2017-09-05 | 中国石油化工股份有限公司 | 深层稠油二氧化碳吞吐注入过程相态变化测试方法 |
| CN107366526A (zh) * | 2017-07-25 | 2017-11-21 | 中国石油天然气股份有限公司 | 一种适合二氧化碳混相驱的致密油藏的筛选方法 |
-
2018
- 2018-01-08 CN CN201810017037.3A patent/CN108252688A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4217955A (en) * | 1978-05-15 | 1980-08-19 | Petroleum Recovery Institute | Oil recovery process |
| CN107130949A (zh) * | 2016-02-25 | 2017-09-05 | 中国石油化工股份有限公司 | 深层稠油二氧化碳吞吐注入过程相态变化测试方法 |
| CN107130960A (zh) * | 2017-06-26 | 2017-09-05 | 中海石油(中国)有限公司 | 一种考虑致密砂岩油藏非线性渗流特征的数值模拟方法 |
| CN107366526A (zh) * | 2017-07-25 | 2017-11-21 | 中国石油天然气股份有限公司 | 一种适合二氧化碳混相驱的致密油藏的筛选方法 |
Non-Patent Citations (3)
| Title |
|---|
| 唐杰: "《致密油注二氧化碳数值模拟研究》", 《中国优秀硕士学位论文全文数据库（电子期刊）》 * |
| 孙仁远 等: "《吸附对页岩中气体临界参数的影响》", 《科学技术与工程》 * |
| 张园 等: "《纳米孔隙中毛管力效应对致密油藏产能的影响》", 《浙江科技学院学报》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113241127A (zh) * | 2021-03-19 | 2021-08-10 | 中国石油大学(北京) | 相平衡模型构建方法、装置、设备及存储介质 |
| CN113241127B (zh) * | 2021-03-19 | 2023-06-02 | 中国石油大学(北京) | 相平衡模型构建方法、装置、设备及存储介质 |
| CN113592194A (zh) * | 2021-08-23 | 2021-11-02 | 北京科技大学 | 建立co2吞吐效果预测模型的方法以及co2吞吐效果评价方法 |
| CN113919448A (zh) * | 2021-12-14 | 2022-01-11 | 武汉大学 | 一种任意时空位置二氧化碳浓度预测影响因素分析方法 |
| CN113919448B (zh) * | 2021-12-14 | 2022-04-01 | 武汉大学 | 一种任意时空位置二氧化碳浓度预测影响因素分析方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2732868C1 (ru) | Способ расчета продуктивности горизонтальных скважин в залежах сланцевого газа при гидроразрыве в условиях нестационарной диффузии | |
| Jing et al. | Mechanism simulation of oil displacement by imbibition in fractured reservoirs | |
| Meng et al. | Influence of effective stress on gas slippage effect of different rank coals | |
| Yu-Long et al. | Performance analysis of fractured wells with stimulated reservoir volume in coal seam reservoirs | |
| CN108133080B (zh) | 一种考虑非稳态吸附的非均质裂缝页岩气藏数值模拟方法 | |
| CN106481332A (zh) | 用于确定页岩气多段压裂水平井内外区动态储量的方法 | |
| CN108252688A (zh) | 致密油藏二氧化碳吞吐影响因素敏感性分析方法及其应用 | |
| Fu et al. | Productivity model with mechanisms of multiple seepage in tight gas reservoir | |
| Shabib-Asl et al. | Evaluation of nanopore confinement during CO 2 huff and puff process in liquid-rich shale formations | |
| CN113406307A (zh) | 一种泥质砂岩储层电阻率指数与相对渗透率的转换方法 | |
| You et al. | Production data analysis of shale gas using fractal model and fuzzy theory: Evaluating fracturing heterogeneity | |
| Lei et al. | Characterization of multiphase flow in shale oil reservoirs considering multiscale porous media by high-resolution numerical simulation | |
| Wang et al. | Monitoring and evaluating the volume fracturing effect of horizontal well | |
| Kong et al. | Experimental evaluation and optimization of improved conductivity of gelling acid etched fractures in deep, low-permeability reservoirs | |
| Wang et al. | The dynamic convection-diffusion effect of CO2 Huff-n-Puff on Pre-Darcy flow behavior | |
| CN115965107A (zh) | 基于岩石蠕变效应的深层页岩裂缝长期导流能力计算方法 | |
| CN114372398B (zh) | 一种考虑储层裂缝闭合效应的聚合物驱试井解释方法 | |
| CN112031719A (zh) | 一种基于流动系数下启动压力的油藏开发方式优选方法 | |
| Wang et al. | Semi-analytical model-based research on influencing factors of the spacing of volume fracturing horizontal well in tight oil reservoirs via consideration of imbibition | |
| CN107481158B (zh) | 一种计算页岩气有效储层改造体积的方法 | |
| CN107247816B (zh) | 一种稠油油藏开展冷采化学降粘最远施工半径的判定方法 | |
| CN106884633B (zh) | 断块油藏注气形成人工气顶合理注气量预测方法 | |
| CN104712299B (zh) | 适合气井控水增气压裂的设计方法 | |
| Wu et al. | Pressure analysis for volume fracturing well in tight reservoirs considering dynamic threshold pressure gradient and stress sensitivity | |
| Guo et al. | Study on gas-liquid relative permeability experiments of fractured-porous reservoirs |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180706 |
|
| RJ01 | Rejection of invention patent application after publication |