2. 中国地质大学(北京)能源学院, 北京 100083
2. School of Energy Resources, China University of Geosciences (Beijing), Beijing, 100083, China
CO2驱油技术是一种有效的提高采收率技术，提高采收率幅度可达15%[1-5]，其中CO2混相驱具有更好的驱油效果[6-9]。对于CO2混相驱，国内外的研究主要集中在最小混相压力[10-12]、驱替机理[13-15]和驱替效率等方面的室内试验[16-17]以及注采参数优化数值模拟等方面[18-19]。油层压力是决定CO2驱油效率的关键因素之一，目前油层压力主要通过井口压力和静液面测试结果来估算，不能实现井底压力实时测定，而且精度不高。文献调研表明，目前尚无矿场注CO2驱油中油层压力连续实时监测方面的系统研究和报道。
胜利油田樊142块于2013年7月开始注CO2驱油，但对该区块的油层压力状况并不十分清楚。为实时了解油层压力特征并实现混相驱油，胜利油田选取了F142-7-X4注采井组，对该井组中的采油井采取了临时关井措施，并在其井底安装了压力计，通过数据线将压力数据传送至地面，实现了井底压力的连续实时监测。笔者根据井底压力实时监测结果，综合运用CO2驱油渗流理论和油藏数值模拟技术，准确地验证了油藏数值模拟得到的井底压力，在此基础上揭示油层内部的压力分布特征和规律，确定CO2混相前缘的推进距离和混相驱面积，为制定最佳注采井措施和开发方案提供依据，从而提高CO2驱油采收率。
1 樊142块概况胜利油田樊142块为低渗透油藏，油层的平均孔隙度为12.5%，平均渗透率为1.89 mD。地面原油密度0.877 1 kg/L，地层原油密度0.746 3 kg/L；地面原油黏度20.4 mPa·s，地下黏度1.18 mPa·s，含硫0.21%，凝固点35.6 ℃，为低黏、低含硫、高凝固点原油。原始油层压力为44.87 MPa，含油饱和度为0.72。其中，F142-7-X4井组为CO2驱油矿场试验井组，包括一口注气井(F142-7-X4井)和6口采油井(见图 1)。该井组的油气开发经历了2个阶段：第一阶段(2007.03-2013.06)，利用天然能量开采；第二阶段(2013.07-2015.08)，F142-7-X4井注CO2，日注入量达20~30 t，为了提高油层压力，其余6口井临时关井。
|
| 图 1 F142-7-X4井组油层顶部深度及井位分布 Fig.1 Depths of oil layer and well location of Well Cluster F142-7-X4 |
通过室内试验研究以及与类似区块的对比分析，认为该区块具备CO2驱油的基本条件，注CO2驱油是可行的。此外，胜利油田已经利用胜利发电厂分离的CO2，在高89块进行了矿场注气试验，为该区块实施注气开采积累了较丰富的矿场经验。
樊142块油层注CO2后能否达到理想的驱油效果关键取决于地层压力，如果不能实现实时压力监测，就难以揭示地层压力的变化规律，无法判断达到混相驱油的时机。因此，需要通过对注采井进行实时压力监测研究油层压力动态变化情况，确定樊142块油层CO2驱油混相区域面积和混相前缘的推进规律。
2 单井井底压力与油层压力监测 2.1 依靠天然能量开采阶段为了准确进行压力监测，采用了矿场井底压力测试和油藏数值模拟相结合的研究手段。前者只能监测到注采井的井底压力，无法直接监测油层内部不同位置的压力，后者不仅可以模拟出井底压力，而且能够模拟出油层各部位压力随时间连续变化的情况。尽管如此，数值模拟出的井底压力必须与监测的井底压力进行对比来验证其可靠性，在油井关闭阶段，井底压力能够代表井底周围油层的压力。为此，笔者利用樊142块的动静态资料，建立了油藏数值模拟地质模型。网格尺寸，x方向为30 m，y方向为30 m，z方向平均为0.92 m；网格数为66×35×4=9 240个。利用多相多组分CO2驱油数值模拟软件Tough-CO2[20]进行了油藏动态历史拟合(该软件是胜利油田勘探开发研究院联合国内外2家研究机构，根据胜利油田油藏和流体的性质联合攻关研发的油藏模拟器)，拟合时间段为井组衰竭式生产阶段与注气井底压力恢复阶段。其中，生产阶段为依靠天然能量开采阶段。天然能量开采末的压力拟合结果如图 2所示。
|
| 图 2 天然能量开采阶段末各生产井井底压力实测值与拟合值的对比 Fig.2 Measured bottom-hole pressures and simulation results in different producers at the end of natural recovery stage |
由图 2可知，依靠天然能量开采的末段，数值模拟的井底压力接近实测井底压力，压力拟合结果较好。在该阶段结束时，6口井的井底压力在24 MPa左右，低于CO2驱油最小混相压力(31 MPa)。与油藏投产时相比，井底压力比原始油层压力低20 MPa左右。由此可见，油层压力低是制约该油藏生产的主控因素之一。
2.2 压力恢复阶段压力恢复阶段仅有F142-7-X4井注气，其余6口生产井均停产，为了跟踪压力变化过程，在6口生产井井底安装了压力计，进行井底压力实时监测。该阶段结束时6口井的井底压力的实测值与拟合结果如图 3所示。
|
| 图 3 压力恢复阶段末井底压力实测值与拟合值的对比 Fig.3 Simulated bottom-hole pressures at the end of wellbore pressure build-up phase |
由图 3可知，压力恢复阶段末，6口井的数值模拟井底压力接近实测井底压力，压力拟合结果较好。压力恢复阶段结束时，F141-1井、F142-6-3井、F142-8-X4井和F142-6-2井等4口井的井底压力高于30 MPa，其中F142-6-2井的井底压力最高，已经高于最小混相压力(31 MPa)，进一步注气可在这4口生产井和注入井之间的油层内实现混相驱；F142-8-3井和F142-7-3井等2口井的井底压力仍然低于25 MPa，储层物性和压力测试分析结果表明，导致这2口井井底压力低的原因是其与注气井间的储层物性差，砂体连通性不好。
2.3 单井压力拟合与对比为对樊142块的开发动态和油层压力进行跟踪研究，利用Tough-CO2软件，对F142-7-X4井组自投产以来的开发历史进行了数值模拟，根据模拟结果，绘制了6口生产井井底压力随开发时间的动态曲线，并添加压力恢复阶段的实测压力数据，结果如图 4所示。
|
| 图 4 6口井的井底压力拟合值与实测值的对比 Fig.4 Simulated and measured bottom-hole pressures in 6 wells |
从图 4可以看出，在依靠天然能量开发阶段，6口井的井底压力下降较快，除F142-8-X4井外，其余5口井在第1年内井底压力急剧下降，其中F142-6-2井和F142-6-3井在生产1个月后，井底压力下降幅度达20 MPa。油层压力经过剧烈下降后，进入相对缓慢的递减阶段，明显体现出低渗透率油藏的特征。
在注CO2开发阶段，关井后油井的井底压力开始恢复，但是F142-8-3井和F142-7-3井实测和模拟的井底压力恢复速度都非常缓慢。根据该区块的地质资料和试井分析资料，可以判断这2口井井底压力恢复慢的主要原因是其与注入井F142-7-X4井之间的地层连通性不好。其余4口井的地层压力恢复较快，其中距注入井F142-7-X4井较近的F142-6-2井和F141-1井的井底压力上升最快，至2015年3月，这2口井的井底压力均高于最小混相压力(31 MPa)。从这4口井井底压力随时间的变化曲线看，注CO2开发后，井底压力呈近似线性关系升高，这是由CO2注入速度保持相对稳定造成的。从拟合效果看，实测和数值模拟的井底压力基本一致，恢复阶段末期井底压力模拟值接近实测值，达到拟合的精度要求，可以在此基础上进行开发指标预测。
3 压力特征及最小混相压力前缘分析 3.1 平面压力分布及前缘推进特征分别选取由数值模拟计算得到的不同时间点的部分地层压力数据，以地质建模的第11-14层为模拟层，下面以主力层第14号模拟层为例，绘制该层8个时间阶段的压力等值线(见图 5)，分析压力恢复阶段压力前缘的运移特征。
|
| 图 5 14号模拟层压力恢复阶段油层压力前缘运移特征 Fig.5 Characteristics of pressure front migration during pressure build-up stage in simulation Layer 14 |
由图 5可知：连续注CO2 14个月后，该层仍然无法达到混相驱油；注CO217个月后，CO2混相压力前缘已传至F142-6-2井，该井和注气井间实现了混相驱油；注CO225个月后，除F142-8-3井外，其余各井的油层压力均已达到或高于混相压力，14号模拟层基本实现了CO2混相驱油。F142-8-3井的试井分析解释结果表明，该井钻遇地层与其他井钻遇的地层连通性差，故该井地层压力上升缓慢，仍未达到最小混相压力。
3.2 垂向剖面油层压力分布及前缘推进特征图 6为F142-7-X4井与F142-6-3井连线方向垂向剖面油层压力前缘推进特征(其中，p为油层压力，MPa；pmm为最小混相压力，MPa)。
|
| 图 6 F142-7-X4井与F142-6-3井连线方向垂向剖面油层压力变化 Fig.6 Changes in formation pressures vertical to profile between Well F142-7-X4 and Well F142-6-3 |
由图 6可知：压力前缘传播快，注气初期便已到达采油井；注气450 d后，注气井附近直径80 m范围内油层压力p已经高于最小混相压力pmm；随着注气时间延长，高于最小混相压力的区域迅速增大；累计注气540 d后，高于最小混相压力的区域已经达到注采井之间的1/3处；累计注气660 d后，油层垂向剖面上的压力均大于最小混相压力，表明CO2在该剖面上波及到位置的驱油属于混相驱；从垂向剖面看，油层的上部压力稍高于下部位置，这主要是油层厚度较小，层内存在不明显的垂向流动所致。
4 混相区域确定方法混相前缘定义为在注气过程中油藏内原油与CO2形成混相区域的最前缘位置，在平面上表现为混相压力前缘与CO2浓度前缘所围区域交集的外边界。形成混相区需同时满足2个条件：区域内压力达到混相压力与注气端注入的CO2已到达该区域。混相区域定义为高于最小混相压力的区域与CO2波及区的交集区域。因此，为了确定混相区域，将CO2浓度分布图和压力分布图叠合，得到的交集区域为CO2混相区域，交集的边界线为混相前缘。其中最小混相压力采用室内细管试验测定值，为31.0 MPa，压力分布是油藏数值模拟的结果，通过对最小混相压力与模拟的油层压力进行对比，并综合考虑CO2浓度分布来确定混相前缘。
图 7所示为压力恢复阶段混相前缘运移特征(其中，物质的量分数为混合物中某一种组分的物质的量与各组分物质的量之和的比)。
|
| 图 7 压力恢复阶段混相前缘运移特征 Fig.7 Characteristics of miscible front migration during pressure build-up stage |
从图 7可以清晰地看出不同时刻混相前缘的推进位置，即混相压力前缘与CO2浓度前缘所围区域交集的外边界。
5 结论1)胜利油田樊142块采取采油井关井和注CO2驱油措施后，油层压力迅速提高，注气15个月后注气井附近就达到了混相驱油条件。
2)应用新建立的CO2混相区域确定方法，判定目前樊142块油层已经具备混相驱油条件，达到油井开井生产的时机。
3)在特低渗透油藏采取注气驱油措施，运用井底压力实时监测和油藏数值模拟相结合的方法，实现了矿场CO2驱油压力实时监测及混相驱油区域的识别，对同类油田实施CO2驱油动态监测有借鉴和指导作用。
4)低渗透油层中CO2驱油最小混相压力前缘监测是目前面临的难点，建议今后开展地球物理监测和设计观察井测试分析方面的研究。
| [1] |
王维波, 师庆三, 余华贵, 等.
二氧化碳驱油注入方式优选实验[J]. 断块油气田 , 2015, 22 (4) : 497–500.
WANG Weibo, SHI Qingsan, YU Huagui, et al. Optimization experiment for carbon dioxide flooding injection mode[J]. Fault-Block Oil & Gas Field , 2015, 22 (4) : 497–500. |
| [2] |
杜勇.
底水稠油油藏水平井复合增产技术[J]. 石油钻探技术 , 2016, 44 (1) : 67–72.
DU Yong. Compound stimulation techniques for heavy oil reservoirs with bottom water[J]. Petroleum Drilling Techniques , 2016, 44 (1) : 67–72. |
| [3] | AZZOLINA N A, NAKLES D V, GORECKI C D, et al. CO2 storage associated with CO2 enhanced oil recovery:a statistical analysis of historical operations[J]. International Journal of Greenhouse Gas Control , 2015, 37 (3) : 384–397. |
| [4] |
沈平平, 廖新维.
二氧化碳地质埋存与提高石油采收率技术[M]. 北京: 石油工业出版社 ,2009 : 151 -168.
SHEN Pingping, LIAO Xinwei. The technology of CO2 geological storage and enhanced oil recovery[M]. Beijing: Petroleum Industry Press , 2009 : 151 -168. |
| [5] |
江怀友, 沈平平, 卢颖, 等.
CO2提高世界油气资源采收率现状研究[J]. 特种油气藏 , 2010, 17 (2) : 5–10.
JIANG Huaiyou, SHEN Pingping, LU Ying, et al. Present situation of enhancing hydrocarbon recovery factor by CO2[J]. Special Oil & Gas Reservoirs , 2010, 17 (2) : 5–10. |
| [6] |
秦积舜, 韩海水, 刘晓蕾.
美国CO2驱油技术应用及启示[J]. 石油勘探与开发 , 2015, 42 (2) : 209–216.
QIN Jishun, HAN Haishui, LIU Xiaolei. Application and enlightenment of carbon dioxide flooding in the United States of America[J]. Petroleum Exploration and Development , 2015, 42 (2) : 209–216. |
| [7] | KOOTTUNGAL L. 2014 worldwide EOR survey[J]. Oil & Gas Journal , 2014, 112 (4) : 79–91. |
| [8] | KOOTTUNGAL L. 2012 worldwide EOR survey[J]. Oil & Gas Journal , 2012, 110 (4) : 57–69. |
| [9] | KOOTTUNGAL L. 2010 worldwide EOR survey[J]. Oil & Gas Journal , 2010, 108 (14) : 41–53. |
| [10] |
鞠斌山, 秦积舜, 李治平, 等.
二氧化碳-原油体系最小混相压力预测模型[J]. 石油学报 , 2012, 33 (2) : 274–277.
JU Binshan, QIN Jishun, LI Zhiping, et al. A prediction model for the minimum miscibility pressure of the CO2-crude oil system[J]. Acta Petrolei Sinica , 2012, 33 (2) : 274–277. |
| [11] |
孙业恒, 吕广忠, 王延芳, 等.
确定CO2最小混相压力的状态方程法[J]. 油气地质与采收率 , 2006, 13 (1) : 82–84.
SUN Yeheng, LYU Guangzhong, WANG Yanfang, et al. A method of state equation for determining minimum miscible pressure of CO2[J]. Petroleum Geology and Recovery Efficiency , 2006, 13 (1) : 82–84. |
| [12] | EISSA M. CO2-oil minimum miscibility pressure model for impure and pure CO2 streams[J]. Journal of Petroleum Science and Engineering , 2006, 58 (1/2) : 173–185. |
| [13] |
杨承志, 岳清山, 沈平平.
混相驱提高石油采收率[M]. 北京: 石油工业出版社 ,1991 : 25 -34.
YANG Chengzhi, YUE Qingshan, SHEN Pingping. Miscible flooding to enhance oil recovery[M]. Beijing: Petroleum Industry Press , 1991 : 25 -34. |
| [14] |
谈士海, 张文正.
非混相CO2驱油在油田增产中的应用[J]. 石油钻探技术 , 2001, 29 (2) : 58–60.
TAN Shihai, ZHANG Wenzheng. Applications of non-mixed phase CO2 for oil displacement in field stimulations[J]. Petroleum Drilling Techniques , 2001, 29 (2) : 58–60. |
| [15] |
沈平平, 江怀友.
温室气体提高采收率的资源化利用及地下埋存[J]. 中国工程科学 , 2009, 11 (5) : 54–59.
SHEN Pingping, JIANG Huaiyou. Utilization of greenhouse gas as resource in EOR and storage it underground[J]. Engineering Science , 2009, 11 (5) : 54–59. |
| [16] |
郝永卯, 薄启炜, 陈月明.
CO2驱油实验研究[J]. 石油勘探与开发 , 2005, 32 (2) : 110–112.
HAO Yongmao, BO Qiwei, CHEN Yueming. Laboratory investigation of CO2 flooding[J]. Petroleum Exploration and Development , 2005, 32 (2) : 110–112. |
| [17] |
李向良, 李振泉, 郭平, 等.
二氧化碳混相驱的长岩心物理模拟[J]. 石油勘探与开发 , 2004, 31 (5) : 102–104.
LI Xiangliang, LI Zhenquan, GUO Ping, et al. Long core physical simulation for CO2 miscible displacement[J]. Petroleum Exploration and Development , 2004, 31 (5) : 102–104. |
| [18] | JU Binshan, FAN Tailiang, JIANG Zaixing. Modeling asphaltene precipitation and flow behavior in the processes of CO2 flood for enhanced oil recovery[J]. Journal of Petroleum Science and Engineering , 2013, 109 (3) : 144–154. |
| [19] |
王杰, 谭保国, 吕广忠.
一种通过数值模拟手段划分CO2驱替相带的新方法:以高89块油藏为例[J]. 科技导报 , 2013, 31 (9) : 46–49.
WANG Jie, TAN Baoguo, LYU Guangzhong. Numerical simulation for division of phase Zone displacement:with Gao 89 as an example[J]. Science & Technology Review , 2013, 31 (9) : 46–49. |
| [20] | JU Binshan, WU Yushu, QIN Jishun. Computer modeling of the displacement behavior of carbon dioxide in undersaturated oil reservoirs[J]. Oil & Gas Science and Technology , 2015, 70 (6) : 951–965. |
