The Influence of Fractures in Shale Oil Reservoirs on CO2 Huff and Puff and Its Pore Production Characteristics
-
摘要: 为了分析页岩水力压裂后产生的裂缝对CO2吞吐效果的影响，采用低场核磁共振测试技术，开展了不同渗透率级别页岩CO2吞吐试验，研究了裂缝对不同渗透率储层CO2吞吐效果的影响。研究发现，裂缝显著提高了CO2吞吐初期的采油速度和采收率；但随着渗透率升高和吞吐次数增多，裂缝对采收率的影响程度逐渐降低；渗透率对裂缝岩心的吞吐效果影响明显小于无裂缝岩心，表明裂缝能够降低渗透率对CO2吞吐采收率的影响；随着吞吐次数增多，裂缝提高大孔隙原油采出程度的幅度减小，而提高小孔隙原油采出程度的幅度增大，但大孔隙仍是原油的主要产出部位。研究结果表明，大孔隙原油的产出主要靠体积膨胀和溶解气驱，速度快且产量大；而小孔隙中原油的产出主要靠抽提和传质方式，过程缓慢且产量小。研究结果为评价裂缝性油藏的产油特征、改善生产动态提供了理论依据。Abstract: To study the influence of fractures from shale hydraulic fracturing on the CO2 huff and puff, low-field nuclear magnetic resonance was adopted in CO2 huff and puff experiments with shales of different permeabilities, and the influence of fractures on the CO2 huff and puff effect in reservoirs of different permeabilities was investigated. Experimental results show that fractures significantly improve the oil production rate and recovery at the initial stage of the CO2 huff and puff operations. However, the influence of fractures on the recovery gradually decreases as huff and puff continues and the permeability increases. The influence of permeability on huff and puff of fractured cores is much lower than that of unfractured cores, which indicates that fractures may reduce the influence of permeability on recovery by CO2 huff and puff. With continuing huff and puff, the increase rate in the degree of reserve recovery of crude oil from macropores contributed by fractures declines while that of crude oil from micropores increases. Despite this finding, macropores continued to be the main contributor of crude oil production. The results indicate that the quick and high production of crude oil from macropores mainly depends on volume expansion and dissolved gas drive. In contrast, the slow and low production of crude oil from micropores is primarily brought by extraction and mass transfer. The results of this study have provided a theoretical basis for evaluating the characteristics of oil production from fractured reservoirs and improving production performance.
-
Keywords:
- shale oil /
- fracture /
- CO2 huff and puff /
- production characteristics /
- nuclear magnetic resonance
-
浊积岩油藏属于典型的深海油藏圈闭，由深水环境中的浊流沉积而成[1-3]。根据沉积特性，浊积岩油藏同一平面储层的孔隙度、渗透率等物性参数具有较大差异，平面非均质性强[4-7]，在实际生产中易出现注入流体窜逸，影响油藏开发效果[8-9]。
气水交替驱和氮气泡沫驱是提高原油采收率的有效方法。气水交替驱可通过气、水段塞的交替注入，在扩大水驱波及系数的同时，实现气体流度控制，从而延缓气体窜逸[10-13]；氮气泡沫驱兼具改善水油流度比、降低油水界面张力、选择性封堵等驱油特性[14-17]，可增大波及体积、提高储量动用程度，达到提高采收率的目的。国内外现有研究主要集中于不同提高采收率方法在纵向非均质储层的适应性方面[18-23]，针对气水交替驱或氮气泡沫驱的效果评价的研究较少[24-25]。张丽娟等人[26]针对纵向非均质储层，通过筛选注入气体种类，提高了气水交替驱的采油效果；M. M. Salehi等人[27]对比了气水交替驱和泡沫驱两种方法，发现泡沫驱可以减缓油水前缘黏性指进；元福卿等人[28]优化了氮气泡沫驱所用泡沫配方，提高了泡沫调堵、分流性能，从而改善了地层纵向非均质性。
不同于纵向非均质性，平面非均质性主要为同一深度储层在平面上由于沉积条件不同而产生的渗透率差异，一般不涉及流体重力差异对驱替介质波及效率的影响。因此，对于纵向非均质储层，可以利用凝胶调剖实现对注入井近井地带的封堵；但对于平面非均质储层，其条带间复杂的连通性导致难以采用近井调剖措施，而需考虑采用具有深度驱替特性的注入介质。从理论角度分析，气水交替驱或氮气泡沫驱的注入介质为气、水/泡沫液段塞，地层注入性优于凝胶类调剖堵水剂，具有深度驱替和整体波及特性，更适于平面非均质油藏。但是，目前针对平面非均质储层调驱效果的研究较少，渗透率级差、高渗条带宽度占比等非均质特性对气水交替驱或氮气泡沫驱效果的影响机制尚不清楚。因此，笔者结合室内试验和数值模拟，开展不同平面非均质条件下的气水交替驱与氮气泡沫驱的适用性及提高采收率机制研究，为浊积岩油藏开发中后期延缓水窜、稳油控水提供了新的技术思路。
1. 油藏地质背景
X油藏为非洲西海岸安哥拉典型浊积岩油藏，属于新近–古近系中新统、渐新统深海沉积储层，含油砂体为海底水道形成的浊积砂体和席状砂，储层含少量泥质胶结物，岩性为长石、石英砂岩，油藏构造和砂体共同控制油气聚集。X油藏由海岸沿斜坡向下（北—南）发育3条浊积水道及席状砂体系。河道储层由2部分砂体构成，中部高净毛比砂岩构成高渗条带，与边缘砂体呈现显著的物性差异，储层平面非均质性强。水道中心轴线以中粗砂为主，物性较好，渗透率达800～3 000 mD；水道边缘天然堤沉积则以砂泥互层为主，电性特征为锯齿状，渗透率通常低于500 mD，中心轴线、边缘沉积及周围深水泥质沉积在平面上构成了3层结构。地质模型计算及前期地质勘探结果表明，高渗条带宽度占砂体整体宽度的1/6～1/3（见图1）。
X油藏平均孔隙度为26.0%，平均气测渗透率为408 mD。地层压力为33 MPa，地层温度为85 ℃，地层原油黏度为0.52～0.79 mPa·s，地面脱气原油密度为0.85～0.86 kg/L。目前采用5注6采井网开发，其中一口注水井X-IF井已转为注气井（见图1），当前油井产出液综合含水率为59%。
2. 岩心驱替试验
基于X油藏中部渗透率较高、边缘渗透率较低和高渗条带宽度占比变化大等地质特征，设计制作代表不同平面非均质性的岩心模型，进行了室内岩心驱替试验，研究了渗透率级差和高渗条带宽度占比对X油藏提高采收率效果与驱油动态特征的影响。
2.1 试验材料与试验设备
由于深海浊积岩油藏天然岩心钻取困难，设计制作了标准尺寸人造方岩心（45 mm×45 mm×300 mm）。根据上述X油藏储层物性和分布特征，设计了2类人造方岩心（见图2），平均孔隙度均为26.0%。考虑X油藏高渗条带宽度占砂体整体宽度的1/6～1/3，2类平面非均质岩心的高渗条带宽度占比分别取上、下极值。其中，X-1型岩心渗透率级差为2.08，模拟X油藏的中部或西部浊积水道储层非均质性较弱的情况（岩心平均渗透率为408 mD，高渗条带的宽度占比为1/3，条带渗透率为624 mD；两侧低渗条带的宽度占比均为1/3，条带渗透率均为300 mD）；X-2型岩心渗透率级差为3.16，模拟X油藏的东部浊积水道储层非均质性较强的情况（岩心模型平均渗透率为408 mD，高渗条带的宽度占比为1/6，条带渗透率为948 mD；两侧低渗条带的宽度占比均为5/12，条带渗透率均为300 mD）。
2类非均质岩心模型的各条带渗透率和孔隙度均符合油藏实际地质认识。人造方岩心选用不同粒径的天然石英砂胶结而成，各条带要满足浊积岩油藏不同区域的渗透率；同时，不同条带间直接接触且相互连通，以模拟平面非均质油藏的条带间复杂连通性。
采用85 ℃下黏度为0.472 mPa·s的原油进行模拟试验，其黏度与地层原油黏度一致；采用由蒸馏水和氯化钠配制而成、矿化度为133 000 mg/L的地层水进行试验。氮气泡沫驱所用的泡沫液采用蒸馏水配制而成，发泡剂选用十二烷基硫酸钠（SDS），质量浓度为4 000 mg/L；稳泡剂选用部分水解聚丙烯酰胺（HPAM，相对分子质量1.2×107，水解度17%，工业品），质量浓度为800 mg/L。气水交替驱和氮气泡沫驱所用氮气纯度为99.95%。气水交替驱或氮气泡沫驱的岩心驱替试验装置主要由ISCO泵（驱替液体）、LF485-FD型气体质量流量控制器（控制注气流速）、压差变送器及数据采集模块、DHZ-50-180型自控恒温箱、高压中间容器、方岩心夹持器（适用岩心规格45 mm×45 mm×300 mm）、气液分离装置、六通阀、试管、量筒和若干管线组成（见图3）。
![]()
图 3 气水交替驱和氮气泡沫驱试验装置Figure 3. Experimental set-up of water and gas alternating flooding and nitrogen foam flooding2.2 试验方法
进行气水交替驱和氮气泡沫驱试验前，根据X油藏注水井的实际注入量，计算油水井中部的线性渗流速率[29-31]，得到试验条件下的注入体积流量。X油藏中部典型井P井的日注入量为4 770 m3/d，油水井距为2 000 m，油层厚度为7.30 m，注入水在油水井中部的线性渗流速率为0.104 m/d，利用式（1）计算出岩心尺度下注入体积流量为0.146 mL/min。
qV=Qi14.4πhDL2 (1) 式中：qV为岩心尺度的注入体积流量，mL/min；Qi为P井的实际注入量，m3/d；h为油层厚度，m；D为单井控制直径，即油水井距，m；L为人造方岩心截面宽度，cm。
2.2.1 气水交替驱试验
1)岩心抽真空、饱和水。将岩心放入抽真空、饱和水的密封钢筒中，用真空泵将压力降至–0.1 MPa并持续抽真空24 h；将模拟地层水注入钢筒中，待压力稳定后，利用手摇泵继续向钢筒加注模拟地层水，直至压力达到10 MPa，饱和24 h。
2）测定岩心孔隙度。根据岩心抽真空、饱和水前后的质量差以地层水密度，计算得到岩心孔隙度。
3）测定饱和油及含油饱和度。用双缸恒流泵以0.05～0.30 mL/min的变流速从岩心夹持器两端反复注入配制好的模拟油，直至出口端产油率达到100%且模拟油注入量达到10倍孔隙体积以上。根据驱替出的水相体积，计算岩心含油饱和度。
4）开始注水驱替，注入体积流量为0.146 mL/min，驱替至岩心出口端含水率达到X油藏当前综合含水率（即59%）后转为气水交替驱；气水交替驱时，气、水注入体积流量均为0.146 mL/min，每个交替注入轮次的注气段塞为0.3倍孔隙体积、注水段塞为0.1倍孔隙体积。
5）岩心出口端含水率达到95%时，停止试验。
2.2.2 氮气泡沫驱试验
1）岩心抽真空、饱和水及饱和油，具体操作过程与气水交替驱试验相同，在完全饱和油后进行后续试验。
2）以注入体积流量0.146 mL/min注水驱替，待岩心出口端含水率达到59%后，转为气–泡沫液交替注入的氮气泡沫驱；氮气泡沫驱时，气、泡沫液的注入体积流量均为0.146 mL/min，每个交替注入轮次注气段塞0.3倍孔隙体积、注泡沫液段塞0.1倍孔隙体积。
3）岩心出口端含水率达到95%时，停止试验。
气水交替驱和氮气泡沫驱试验过程中，实时记录各时间段的产液量、产油量和驱替压差等产出参数，并计算瞬时含水率与原油采出程度。
2.3 试验结果分析
2.3.1 弱非均质条件对驱替效果的影响
选用X-1型岩心进行弱非均质条件下的气水交替驱和氮气泡沫驱试验，岩心渗透率级差为2.08，岩心编号、注入方式及驱油结果见表1。
表 1 不同注入方式X-1型岩心提高采出程度结果Table 1. Oil recovery percent from different injection methods of X-1 core岩心
编号初始含油饱和度，
%采出程度，% 注入方式 前期水驱 最终 X-1-1 59.65 45.75 57.20 气水交替驱 X-1-2 59.40 44.16 62.42 氮气泡沫驱 分析表1可知，采用气水交替驱和氮气泡沫驱均可提高采出程度，且氮气泡沫驱效果略优于气水交替驱，相比提高了5.22百分点。对于气水交替驱，注入水后，孔隙中的水相饱和度增加，气相相对渗透率降低，可以在孔隙尺度减缓气体窜逸，同时注入氮气段塞可以在水驱基础上进一步增大波及体积；由于泡沫具有洗油和调驱双重功能，氮气泡沫驱既能提高洗油效率，又能有效封堵高渗条带并提高波及系数，采出程度提高幅度较大。
X-1型岩心气水交替驱和氮气泡沫驱试验的采出程度、驱替压差及出口端含水率随注入体积的变化曲线如图4—图6所示。
![]()
图 4 不同注采方式下X-1型岩心采出程度随注入体积的变化关系Figure 4. Relationship of oil recovery percentage and injected volume in different injection methods of X-1 core![]()
图 6 不同注采方式下X-1型岩心出口端含水率随注入体积的变化关系Figure 6. Relationship of water cut and injected volume in different injection methods of X-1 core由图4可知，在前期水驱结束后，气水交替驱和氮气泡沫驱在注入体积小于0.7倍孔隙体积时，采出程度曲线几乎一致，且上升较为缓慢。这是由于二者在第一轮次（将一个注气段塞与一个注水或泡沫液段塞的组合视为一个注入轮次）中均为先注入氮气，再注入水或泡沫液。第一轮次注入氮气过程中，气相为连续相，与水驱过程相近，采出程度上升幅度较小；但注入0.7～0.9倍孔隙体积流体时，二者的采出程度出现明显跃升。研究表明，随着注入体积增大，气水交替驱和氮气泡沫驱的贾敏效应提高了相对低渗条带的原油动用程度，该部分原油逐渐被驱替至出口端，产油量明显上升。由于氮气泡沫驱的贾敏效应更强、且兼具表面活性剂的洗油作用，其采出程度的上升幅度更为明显。
由图5可知，气水交替驱过程中，注气时驱替压差先上升后迅速下降，注水时驱替压差保持平稳，与前期水驱接近。这是由于注气时，贾敏效应的存在使驱替压差大幅度增大，扩大了驱替剂的波及范围；但随着气体继续注入，形成气流通道，气体变为连续相，贾敏效应减弱，出现气窜现象，驱替压差迅速下降。后续注水时，注入水主要流入气窜大孔道，填补地层能量，驱替压差变化较小。氮气泡沫驱过程中，驱替压差随着注入轮次增加而逐渐上升。注入泡沫液后，由于其黏度较高，驱替压差上升，随后注气在岩心中形成氮气泡沫，封堵高渗条带。此外，第二轮次和第三轮次注入泡沫液后，后续气驱驱替压差均明显高于气水交替驱，表明氮气泡沫封堵高渗条带的效果优于气水交替驱。
![]()
图 5 不同注采方式下X-1型岩心驱替压差随注入体积的变化关系Figure 5. Relationship of displacement pressure and injected volume in different injection methods of X-1 core由图6可知，气水交替驱和氮气泡沫驱过程中，岩心出口端含水率曲线的整体变化趋势基本一致，且在注入体积小于0.7倍孔隙体积时，出口端含水率继续上升；注入0.7～0.9倍孔隙体积流体时，含水率明显下降，曲线呈漏斗形。结合图5中的驱替压差分析认为，对于气水交替驱，由于驱替压差增大，非均质岩心中相对低渗条带的部分原油被动用，产油量增加，导致出口端含水率降低；但在后续注入轮次中，随着驱替压差下降，相对低渗条带的原油动用逐渐困难，出口端含水率逐渐上升。对于氮气泡沫驱，其驱替压差上升幅度较大，调驱效果明显，且泡沫可以提高洗油效率，使产油量大幅度上升，出口端含水率明显降低。但由于平面非均质性较弱，氮气泡沫驱过程中所形成的泡沫沿前缘地带均匀推进，未能有效封堵高渗条带，后续注入轮次中氮气突破泡沫封锁，发生气窜，出口端含水率逐渐上升。
气水交替驱和氮气泡沫驱的前2个注入轮次中，氮气泡沫驱稳定出口端含水率的效果明显优于气水交替驱。但随着注入轮次增加，2种方法在第三轮次的效果均明显下降，表明弱平面非均质条件下气水交替驱和氮气泡沫驱的增油期主要集中在前2个注入轮次。
2.3.2 强非均质条件对驱替效果的影响
选用X-2型岩心模型开展了强非均质条件下的气水交替驱和氮气泡沫驱试验，渗透率级差为3.16，所用岩心编号、注入方式和驱油结果见表2。
表 2 不同注入方式下X-2型岩心提高采出程度结果Table 2. Oil recovery percentage from different injection methods of X-2 core岩心
编号初始含油饱和度，
%采出程度，% 注入方式 前期水驱 最终 X-2-1 52.91 28.66 41.46 气水交替驱 X-2-2 53.27 29.25 74.33 氮气泡沫驱 分析表2可知，氮气泡沫驱提高原油采出程度的效果明显优于气水交替驱，相比提高了31.87百分点。这是因为，X-2型岩心高渗条带宽度较窄且渗透率更高，平面非均质性更强，气水交替驱虽然能够提高波及系数，但其调整储层非均质的能力有限，难以有效控制强平面非均质条件下的气体窜逸，使得最终采出程度大幅度降低。然而，氮气泡沫驱对于非均质油藏的适用性更强，可以延缓强非均质条件下高渗条带中的气体窜逸，使低渗条带的原油得到有效动用，同时泡沫可以提高洗油效率，使最终采出程度显著提高。对比表1和表2可以看出，气水交替驱较适用于弱非均质条件，在强非均质条件下的适用性较差；而氮气泡沫驱在弱非均质和强非均质条件下都适用。
X-2型岩心气水交替驱和氮气泡沫驱试验的采出程度、驱替压差及出口端含水率随注入体积的变化曲线如图7—图9所示。
![]()
图 7 不同注采方式下X-2型岩心采出程度随注入体积的变化关系Figure 7. Relationship of oil recovery percentage and injected volume in different injection methods of X-2 core![]()
图 8 不同注采方式下X-2型岩心驱替压差随注入体积的变化关系Figure 8. Relationship of displacement pressure and injected volume in different injection methods of X-2 core![]()
图 9 不同注采方式下X-2型岩心出口端含水率随注入体积的变化关系Figure 9. Relationship of water cut and injected volume in different injection methods of X-2 core由图7可知，注入体积小于0.6倍孔隙体积时，气水交替驱和氮气泡沫驱的采出程度曲线变化几乎一致，其原因与图4相同。后续注入轮次中，气水交替驱的采出程度仅在第二轮次有明显提升，整体提高幅度较小；氮气泡沫驱每个注入轮次的采出程度均有明显提高，驱替效果显著。平面非均质性较强时，气水交替驱在高渗条带中的气窜现象较为明显，提高采出程度的效果降低；由于强非均质条件下高渗条带的孔隙半径较大，氮气泡沫的稳定性增强[32-33]，延缓窜逸效果更优，使得低渗条带得到有效动用，因此采出程度明显提高。
由图8可知，气水交替驱过程中，驱替压差波动较小，且在2个注入轮次后迅速下降，表明在强平面非均质条件下，气水交替驱的波及范围主要集中在高渗条带，低渗条带难以被波及，因此提高采出程度的效果较差；氮气泡沫驱过程中，驱替压差随着注入轮次增加而逐渐上升，表明在强平面非均质条件下氮气泡沫可以封堵高渗条带的大孔隙，迫使后续注入流体转向并波及低渗条带，从而提高低渗条带的原油动用程度。
由图9可知，气水交替驱和氮气泡沫驱的出口端含水率曲线呈现明显差异。气水交替驱的前2个注入轮次中，岩心出口端含水率先明显降低后迅速上升，表明强非均质条件下气水交替驱可以降低出口端含水率，但稳定性较差，仅能维持一个注入轮次；氮气泡沫驱过程中，出口端含水率先下降至0，后稳定在50%附近，低含水阶段能够持续4个注入轮次，说明氮气泡沫驱降低并稳定岩心出口端含水率的效果显著。这是由于随着注入体积增大，岩心中泡沫含量增加，封堵高渗条带大孔道的效果提升，氮气泡沫驱波及范围逐渐从高渗条带延伸到低渗条带，从而降低并稳定了岩心出口端含水率。
3. 气水交替驱/氮气泡沫驱数值模拟研究
由于岩心驱替试验难以可视化描述注入介质在岩心中的波及规律，因此在岩心试验基础上进行了气水交替驱和氮气泡沫驱的数值模拟研究。以强非均质岩心为例，建立了岩心尺度的数值模拟模型，其中的模型尺寸、孔渗物性以及初始含油、含水饱和度均与岩心试验参数保持一致；通过调整油–水、油–气相对渗透率曲线和毛管力曲线等，实现气水交替驱或氮气泡沫驱试验结果的历史拟合，进而通过分析驱替过程中不同渗透率条带的含油饱和度的变化特征，揭示不同提高采收率方法提高采收率的机制。
利用CMG数值模拟软件的STARS模块进行岩心尺度数值模拟，单位制选用LAB试验单位制。模型选用正交网格，i，j和z方向网格数分别为32，3和1个，i方向网格尺寸为1.000 cm，j方向3个网格尺寸分别为18.75，7.50和18.75 mm，z方向网格尺寸为45.00 mm。各条带渗透率分别为300，948和300 mD。设定模型左侧为注入端，设置一口注水井、一口注气井，右侧为出口端，设置一口生产井；岩心数值模型尺寸、孔隙度、含油/含水饱和度、渗透率分布均与试验参数保持一致。模型初始压力和生产井井底流压均设定为33 MPa，以模拟岩石试验的围压和出口端回压条件。
选用机理法进行氮气泡沫驱数值模拟，利用反应式表示泡沫的生成和破灭（S表示表面活性剂，L表示液膜）：
0.99993660682H2O+6.33931845830×10−5S+1.0N2→1.0L+1.0N2 (2) 1.0L+1.0N2→0.99993660682H2O+6.33931845830×10−5S+1.0N2 (3) 式（2）中，左侧表示水、表面活性剂与氮气的反应，右侧表示反应形成液膜（泡沫）。其中，式（2）左右两侧均含有氮气，左侧氮气表示氮气作为反应物参与反应，右侧氮气则表示生成的泡沫具有氮气的组分特性。
根据气水交替驱和氮气泡沫驱的岩心试验步骤，先注水驱替300 min，然后开始气水交替驱或氮气泡沫驱。每个注入轮次均先注375 min的氮气再注125 min的水或表面活性剂溶液，以此循环注入，驱替至2175 min停止模拟运算。气水交替驱和氮气泡沫驱采出程度的数值模拟历史拟合结果如图10（a）、图10（b）所示。由图10可知，氮气泡沫驱的开发效果远优于气水交替驱。
![]()
图 10 气水交替驱和氮气泡沫驱采出程度数值模拟历史拟合结果Figure 10. History matching of recovery percentage from water and gas alternating flooding and nitrogen foam flooding由10（a）可知，对于气水交替驱，前2个轮次的采出程度试验数据与模拟数据拟合效果较好，表明可以用数值模拟结果研究气水交替驱的提高采收率机制。利用该岩心尺度数值模型，对比开展注水驱替数值模拟。水驱和气水交替驱过程中，不同渗透率条带的含油饱和度变化特征对比如图11所示。水驱过程中，注入水主要波及高渗条带（见图11（a）所示的岩心模型中间层网格），导致含油饱和度降低，但即使继续注水驱替，低渗条带（见图11（a）所示的岩心模型上、下两层网格）的含油饱和度未发生明显变化，表明注入水主要沿高渗条带低效循环，需要进行液流转向措施。对于气水交替驱，气水交替注入的初始时刻（t=300 min），高渗条带的含油饱和度较低，剩余油主要赋存在低渗条带；随着气水交替注入轮次增多，低渗条带含油饱和度逐渐降低，且降低幅度比水驱方案大，表明气水交替驱在平面非均质岩心模型中发挥了扩大波及体积的作用，能够驱动剩余油。
![]()
图 11 水驱与气水交替驱过程中含油饱和度变化特征对比Figure 11. Comparison of oil saturation distribution change during water flooding and water and gas alternating flooding由图10（b）可知，氮气泡沫驱驱替过程中的采出程度试验数据与模拟数据拟合效果较好，表明可以用数值模拟结果研究氮气泡沫驱的驱替规律。对比注水驱替数值模拟，氮气泡沫驱过程中不同渗透率条带的含油饱和度变化特征如图12所示。初始时刻（t=300 min），高渗条带（见图12（b）所示的岩心模型中间层网格）含油饱和度较低，剩余油主要赋存于低渗条带（见图12（b）所示的岩心模型上、下两层网格）。随着氮气泡沫驱的进行，高渗条带含油饱和度基本不变，表明氮气泡沫发挥了堵大不堵小、堵水不堵油的性能，迫使后续注入介质进入低渗条带。分析低渗条带含油饱和度变化特征可知，氮气泡沫在低渗条带近乎呈活塞式驱替，将低渗条带中的剩余油采出，因此能够大幅度提高采收率，这与岩心驱替试验的分析结果基本一致，表明氮气泡沫驱可以作为平面非均质油藏剩余油挖潜的有效手段。
![]()
图 12 水驱与氮气泡沫驱过程中含油饱和度变化特征对比Figure 12. Comparison of oil saturation distribution change during water flooding and nitrogen foam flooding4. 结 论
1）驱替试验结果表明，氮气泡沫驱效果优于气水交替驱。气水交替驱的主要增油期为前1～2个轮次，而氮气泡沫驱增油期可以持续2～4个轮次，稳油控水效果更优。
2）岩心模型渗透率级差较小时，气水交替驱表现出良好的提高波及系数和降低出口端含水率的能力，但当渗透率级差较大时，该方法控制气体流度的能力降低，驱替压差上升幅度较小，提高采收率效果变差。
3）岩心模型渗透率级差较大时，氮气泡沫驱仍可发挥泡沫堵大不堵小、堵水不堵油及表面活性剂洗油的多重特性，抑制高渗条带中的流体窜逸，使得氮气泡沫在低渗条带近似呈活塞式驱替，从而实现深部调驱，大幅度提高采收率。
4）气水交替驱及氮气泡沫驱数值模拟结果与试验结果拟合度较好，综合岩心驱替试验和数值模拟结果可知，气水交替驱较适用于弱非均质条件，但在强非均质条件下适用性较差；氮气泡沫驱可同时适用于弱非均质和强非均质条件，有效提高深水浊积岩油藏采收率，为该类油藏的经济高效开发提供技术支持。
-
![]()
图 2 不同CO2溶解比例下原油膨胀系数随平衡压力的变化
Figure 2. Variation of the expansion coefficient of crude oil with different proportions of dissolved CO2 with equilibrium pressure
![]()
图 4 第一轮降压衰竭中1#、2#岩心采油速度及采出程度对比
Figure 4. Comparison of oil production rate and degree of reserve recovery between Core 1# and Core 2# in the first cycle of depressurization
![]()
图 5 1#、2#岩心饱和油及每轮吞吐后的T2谱分布
Figure 5. T2 spectrum distribution of Core 1# and Core 2# saturated oil and after each huff and puff cycle
![]()
图 6 1#、2#岩心每轮吞吐后不同孔径孔隙原油采出程度对比
Figure 6. Comparison of degree of reserve recovery of crude oil for Core 1# and Core 2# after each huff and puff cycle from pores in different sizes
![]()
图 7 3#和5#岩心饱和油及每轮吞吐后的T2谱分布
Figure 7. T2 spectrum distribution of Core 3# and Core 5# saturated oil and after each huff and puff cycle
![]()
图 8 含有裂缝的不同渗透率岩心大孔隙、小孔隙采出程度对比
Figure 8. Comparison between degree of reserve recovery from macropores and micropores in fractured cores with different permeability
表 1 试验岩心基本参数
Table 1 Basic parameters of experimental cores
岩心编号 长度/mm 直径/mm 孔隙度, % 渗透率/mD 有机碳含量,% 平均孔隙半径/nm 裂缝横向长度/mm 裂缝纵向长度/mm 1# 80.3 25.1 9.45 0.007 5 3.77 11.72 80.3 21.6 2# 81.3 25.2 8.76 0.005 2 3.59 8.13 无裂缝 无裂缝 3# 78.4 25.1 8.95 0.027 0 3.08 16.81 78.4 22.4 4# 79.2 25.2 9.57 0.044 0 2.74 19.24 无裂缝 无裂缝 5# 81.1 25.1 10.22 0.086 0 1.66 27.43 81.1 22.1 6# 80.6 25.1 11.38 0.094 0 2.14 24.86 无裂缝 无裂缝 
下载: 导出CSV
表 2 试验岩心CO2吞吐采收率
Table 2 CO2 huff and puff recovery of experimental cores
岩心编号 吞吐采收率,% 累计吞吐
采收率，%第一轮 第二轮 第三轮 1# 20.6 8.8 2.4 31.8 2# 11.4 2.8 1.4 15.6 3# 22.4 11.6 4.3 38.3 4# 16.2 6.7 2.2 25.1 5# 25.8 14.1 4.7 44.6 6# 19.6 8.4 3.2 31.2
下载: 导出CSV
表 3 不同吞吐轮次中裂缝对大、小孔隙采出程度提高幅度的影响
Table 3 Influence on the increase rate in degree of reserve recovery from macropores and micropores contributed by fractures in each huff and puff cycle
吞吐轮次 孔隙类型 采出程度，% 采出程度提高
幅度，%1#岩心 2#岩心 第一轮 小孔隙 6.8 4.4 54.5 大孔隙 37.5 24.6 52.4 第二轮 小孔隙 4.3 1.8 138.9 大孔隙 19.6 13.8 42.0 第三轮 小孔隙 1.8 0.6 200.0 大孔隙 10.5 8.7 20.7
下载: 导出CSV
-
[1] 张林晔,李钜源,李政,等. 北美页岩油气研究进展及对中国陆相页岩油气勘探的思考[J]. 地球科学进展,2014,29(6):700–711. doi: 10.11867/j.issn.1001-8166.2014.06.0700 ZHANG Linye, LI Juyuan, LI Zheng, et al. Advances in shale oil/gas research in North America and considerations on exploration for continental shale oil/gas in China[J]. Advances in Earth Science, 2014, 29(6): 700–711. doi: 10.11867/j.issn.1001-8166.2014.06.0700
[2] 卢双舫,李俊乾,张鹏飞,等. 页岩油储集层微观孔喉分类与分级评价[J]. 石油勘探与开发,2018,45(3):436–444. LU Shuangfang, LI Junqian, ZHANG Pengfei, et al. Classification of microscopic pore-throats and the grading evaluation on shale oil reservoirs[J]. Petroleum Exploration and Development, 2018, 45(3): 436–444.
[3] 陈作,刘红磊,李英杰,等. 国内外页岩油储层改造技术现状及发展建议[J]. 石油钻探技术,2021,49(4):1–7. doi: 10.11911/syztjs.2021081 CHEN Zuo, LIU Honglei, LI Yingjie, et al. The current status and development suggestions for shale oil reservoir stimulation at home and abroad[J]. Petroleum Drilling Techniques, 2021, 49(4): 1–7. doi: 10.11911/syztjs.2021081
[4] 胡伟,吕成远,王锐,等. 水驱转CO2混相驱渗流机理及传质特征[J]. 石油学报,2018,39(2):201–207. doi: 10.7623/syxb201802008 HU Wei, LYU Chengyuan, WANG Rui, et al. Porous flow mechanisms and mass transfer characteristics of CO2 miscible flooding after water flooding[J]. Acta Petrolei Sinica, 2018, 39(2): 201–207. doi: 10.7623/syxb201802008
[5] WANG Haitao, LUN Zengmin, LYU Chengyuan, et al. Nuclear-magnetic-resonance study on mechanisms of oil mobilization in tight sandstone reservoir exposed to carbon dioxide[J]. SPE Journal, 2018, 23(3): 750–761. doi: 10.2118/179554-PA
[6] TORABI F, QAZVINI FIROUZ A, KAVOUSI A, et al. Comparative evaluation of immiscible, near miscible and miscible CO2 huff-n-puff to enhance oil recovery from a single matrix-fracture system (experimental and simulation studies)[J]. Fuel, 2012, 93: 443–453. doi: 10.1016/j.fuel.2011.08.037
[7] 张怿赫,盛家平,李情霞,等. CO2吞吐技术应用进展[J]. 特种油气藏,2021,28(6):1–10. ZHANG Yihe, SHENG Jiaping, Li Qingxia, et al. Advances in the application of CO2 stimulation technology[J]. Special Oil & Gas Reservoirs, 2021, 28(6): 1–10.
[8] 杨正明,刘学伟,张仲宏,等. 致密油藏分段压裂水平井注二氧化碳吞吐物理模拟[J]. 石油学报,2015,36(6):724–729. doi: 10.7623/syxb201506009 YANG Zhengming, LIU Xuewei, ZHANG Zhonghong, et al. Physical simulation of staged-fracturing horizontal wells using CO2 huff and puff in tight oil reservoirs[J]. Acta Petrolei Sinica, 2015, 36(6): 724–729. doi: 10.7623/syxb201506009
[9] 钱坤,杨胜来,马轩,等. 超低渗透油藏CO2吞吐利用率实验研究[J]. 石油钻探技术,2018,46(6):77–81. QIAN Kun, YANG Shenglai, MA Xuan, et al. CO2 utilization ratio simulation during a CO2 huff-and-puff process in ultra-low permeability oil reservoirs[J]. Petroleum Drilling Techniques, 2018, 46(6): 77–81.
[10] QIAN Kun, YANG Shenglai, DOU Hongen, et al. Experimental investigation on microscopic residual oil distribution during CO2 huff-and-puff process in tight oil reservoirs[J]. Energies, 2018, 11(10): 2843. doi: 10.3390/en11102843
[11] 黄小亮,贾新峰,周翔,等. 延长油田低渗油藏长岩心二氧化碳吞吐参数优化[J]. 新疆石油地质,2015,36(3):313–316. HUANG Xiaoliang, JIA Xinfeng, ZHOU Xiang, et al. Optimization of long core CO2 huff-n-puff experimental parameters in low permeability reservoir, Yanchang Oilfield[J]. Xinjiang Petroleum Geology, 2015, 36(3): 313–316.
[12] XIAO Pufu, YANG Zhengming, WANG Xuewu, et al. Experimental investigation on CO2 injection in the Daqing extra/ultra-low permeability reservoir[J]. Journal of Petroleum Science and Engineering, 2017, 149: 765–771. doi: 10.1016/j.petrol.2016.11.020
[13] 霍刚,范潇. 混注烟道气辅助蒸汽吞吐驱替机理数值模拟研究[J]. 油气地质与采收率,2012,19(4):59–61,65. doi: 10.3969/j.issn.1009-9603.2012.04.017 HUO Gang, FAN Xiao. Numerical simulation study on displacement mechanism of flue gas assisting steam huff and puff[J]. Petroleum Geology and Recovery Efficiency, 2012, 19(4): 59–61,65. doi: 10.3969/j.issn.1009-9603.2012.04.017
[14] ABEDINI A, TORABI F. Oil recovery performance of immiscible and miscible CO2 huff-and-puff processes[J]. Energy & Fuels, 2014, 28(2): 774–784.
[15] 李二党,韩作为,高祥瑞,等. 不同注气介质驱替致密油藏微观孔隙动用特征研究[J]. 石油钻探技术,2020,48(5):85–91. doi: 10.11911/syztjs.2020078 LI Erdang, HAN Zuowei, GAO Xiangrui, et al. Research on the microscopic pore producing characteristics of tight reservoirs displaced by different gas injection media[J]. Petroleum Drilling Techniques, 2020, 48(5): 85–91. doi: 10.11911/syztjs.2020078
[16] WEI Bing, GAO Ke, SONG Tao, et al. Nuclear-magnetic-resonance monitoring of mass exchange in a low-permeability matrix/fracture model during CO2 cyclic injection: a mechanistic study[J]. SPE Journal, 2020, 25(1): 440–450. doi: 10.2118/199345-PA
[17] 唐人选,梁珀,吴公益,等. 苏北复杂断块油藏二氧化碳驱油效果影响因素分析及认识[J]. 石油钻探技术,2020,48(1):98–103. doi: 10.11911/syztjs.2019125 TANG Renxuan, LIANG Po, WU Gongyi, et al. Analyzing and understanding the influencing factors of CO2 flooding in the Subei complex fault block reservoirs[J]. Petroleum Drilling Techniques, 2020, 48(1): 98–103. doi: 10.11911/syztjs.2019125
[18] WEI Bing, LU Laiming, PU Wanfen, et al. Production dynamics of CO2 cyclic injection and CO2 sequestration in tight porous media of Lucaogou Formation in Jimsar Sag[J]. Journal Petroleum Science and Engineer, 2017, 157: 1084–1094. doi: 10.1016/j.petrol.2017.08.023
[19] 贾瑞轩,孙灵辉,苏致新,等. 二氧化碳吞吐致密油藏的可动用性[J]. 断块油气田,2020,27(4):504–508. JIA Ruixuan, SUN Linghui, SU Zhixin, et al. Availability of CO2 huff and puff in tight reservoir[J]. Fault-Block Oil & Gas Field, 2020, 27(4): 504–508.
[20] 巩联浩,刘继梓,武兴,等. 裂缝性致密油藏二氧化碳吞吐基质-裂缝间流体渗流特征研究[J]. 特种油气藏,2021,28(1):118–124. doi: 10.3969/j.issn.1006-6535.2021.01.017 GONG Lianhao, LIU Jizi, WU Xing, et al. Study on seepage characteristics of fluid between matrix and fracture in CO2 huff-puff process in fractured tight reservoirs[J]. Special Oil & Gas Reservoirs, 2021, 28(1): 118–124. doi: 10.3969/j.issn.1006-6535.2021.01.017
-
期刊类型引用(4)
1. 聂法健，郭志华，吴罕，侯丽莉，李真，刘振洋，王艳平. 浊积岩油藏优势渗流通道精准识别——以胜利油区牛庄洼陷牛20区块为例. 石油地质与工程. 2024(03): 9-13 . 
百度学术
2. 李舒展，杨进，朱国倞，李中，王宁. 深水钻井井口吸力锚最小下入深度预测方法. 石油钻探技术. 2023(02): 29-36 . 本站查看
3. 葛罗. 大庆油田萨北区块中渗透砂岩油藏凝胶调剖剂运移吸附试验研究. 石油钻探技术. 2023(03): 119-125 . 本站查看
4. 苑登御. N_2泡沫/CO_2复合吞吐提高采收率三维物理模拟试验研究. 石油钻探技术. 2022(06): 126-132 . 本站查看
其他类型引用(1)
计量
- 文章访问数: 507
- HTML全文浏览量: 274
- PDF下载量: 85
- 被引次数: 5
