当前中国石油公司开发业务优质资源接替面临较大挑战,新增储量以低渗透、稠油和特殊岩性油藏等低品位资源为主,产能建设投资巨大、开发效益变差。而以中高渗透为主体的老油田相对是优质资源,储量规模大,目前已进入高含水、高采出程度的“双高”开发阶段。高含水老油田的产量对中石油的总产量起决定性作用,进一步挖掘高含水油田的潜力,提高油田采收率对中石油的长期可持续发展具有重要的战略意义[1-2]。但对于注水开发层状砂岩油藏而言,储层非均质性比较严重,在油田进入高含水后期,层间矛盾、层内矛盾、平面矛盾日益突出,呈现了注采关系不清、水流优势通道更为突显、水驱控制程度低等问题,对油田的稳油控水造成了极大困难[3]。通过对渤海湾中深层油藏储层的精细认识[4]及实际生产动态资料的分析,发现与油层相通的干层具有一定吸水能力,利用这样的有效干层完善注采关系、改变水流方向能够起到较好改善开发效果的作用。
目前,与油层相通的干层有效性识别还没有专门深入系统地研究,对于干层的利用也只是针对某个油田一些实际生产现状的分析,并没有系统地研究有效干层的判别标准及如何系统地分析利用干层对注采关系的影响。本文针对渤海湾中深层油藏特点,采用主成分分析统计方法与模糊贴近度方法相结合,在如何识别有效干层和利用干层完善注采井网调整方面提供了一种方法,以指导渤海湾中深层油藏在高含水期更好的完善注采关系,最终提高水驱采收率。
1 有效干层对注采关系的影响在陆相油田储层常出现纵向上砂泥岩频繁互层,加之厚、薄砂体连通关系的复杂性以及单成因砂体自身的岩性和物性的空间差异,导致了储层复杂的结构非均质性强,致使部分储层砂体孔喉结构细,砂体致密,在生产上表现出储层砂体出油能力差异大,甚至部分砂体在试油或生产时不出油,也就是所称谓的干砂层[5]。这些所谓的干层在油田生产过程中,受油田多层合注合采,层间非均质性的影响,抑制了部分物性相对较差储层砂体自身能力的发挥,被认为是干储层[6-7]。但是,对于注水井就不尽然,通过大量油田实际吸水测试资料可以看出,部分与油层相通的干砂层具有一定的吸水能力,这样的干层称之为有效干层。开发时可以利用有效干层以达到提高水驱控制程度的效果。
在生产过程中,受储层非均质性的影响,注入水会始终趋向砂体比较厚、物性较好的层段,对那些储层物性相对较差的砂体可能就会有大量的剩余油,并且长期对储层好的层段进行注水开发后,容易出现水淹,形成水流优势通道。因此对于这种储层非均质性严重的油田,充分利用有效干层是很有必要的,通过对与油层相通干层的有效性识别,在有效干层层位注水井进行补孔注水,能够提高井网对砂体的控制程度,同时大大改善注采连通受效关系(图 1),以最大可能挖潜剩余油,达到提高石油采收率的目的[8-9]。
判别干层是否吸水的影响因素多而复杂,首先,水井干层能否吸水的前提条件是该层是否与相应的油井相连;其次,在此基础上,需进一步判别相连通的干层是否能够满足注采驱替。传统的干层识别方法在实际应用过程中常遇到输入参数过多或输入参数过少等问题,并且各参数之间还存在相互关联,导致信息重叠与数据冗余,影响预测的精度。因此,对干层性质影响因素进行优化是预测精度的重要保证。
一般识别储层的参数包括孔隙度、渗透率、含油饱和度、有效厚度、孔喉半径、微电位、声波时差、电阻率、自然伽马、泥质含量、产能系数、储能系数、可动水指数等等,在这些参数之间存在一定的相关性。通过分析这些参数可以识别储层,获得储层物性情况。在实际生产过程中,储层对油水井的贡献受多重因素影响,特别是一些相对较差的储层,因此,对注采的影响难以给出定量判断。鉴于此,为了更精细地进行高含水期注采调控,完善注采关系[10-11],结合实际生产过程情况,应用主成分分析方法对影响因素进行筛选,得出综合影响因子,进而为干层有效利用的判别奠定基础,也为进一步完善注采关系提供依据和条件[12-13]。
上述的影响因素具有不同的量纲,对于主成分分析方法来说,方差起到至关重要的影响,因此,为消除量纲不同可能带来的一些不合理影响,采用变量标准化的方法对参数进行处理。假设评价干层的有效性影响因素有
$ {E(}X_{i}^{*})=\mu _{i}^{*} $ | (1) |
$ {\rm Var}(X_{i}^{*})={{(\sigma _{i}^{*})}^{2}} $ | (2) |
标准化处理后的参数为
$ {{X}_{i}}=\dfrac{X_{i}^{*}-\mu _{i}^{*}}{\sigma _{i}^{*}}, {\kern 10pt}i=1, \cdots , p $ | (3) |
即处理后的判别参数为
$ \left\{ \begin{matrix}{{Z}_{1}}={{{{a}}}_{1}'}X={{a}_{11}}{{X}_{1}}+{{a}_{21}}{{X}_{2}}+L+{{a}_{p1}}{{X}_{p}} \\{{Z}_{2}}={{{{a}}}_{2}'}X={{a}_{12}}{{X}_{1}}+{{a}_{22}}{{X}_{2}}+L+{{a}_{p2}}{{X}_{p}} \\\cdots \\{{Z}_{p}}={{{{a}}}_{p}'}X={{a}_{1p}}{{X}_{1}}+{{a}_{2p}}{{X}_{2}}+L+{{a}_{pp}}{{X}_{p}} \\\end{matrix} \right. $ | (4) |
可得
$\left\{ \begin{array}{l} {\rm Var}({{Z}_{i}})={{{{a}}}_{i}'}\sum {{a}_{i}} {\kern 28pt}i=1, 2, \cdots, p \\ {\rm Cov}({{Z}_{i}}, {{Z}_{j}})={{{{a}}}_{i}'}\sum {{a}_{j}} {\kern 10pt} i, j=1, 2, \cdots, p \\ \end{array} \right. $ | (5) |
若希望用
$ \left\{ \begin{array}{l} {{{{a}}}_{i}'}{{a}_{i}}=1, {\kern 10pt} i=1, 2, \cdots, p \\{{{{a}'}}_{i}}\sum {{a}_{j}}=0, {\kern 10pt} {i>1}, j=1, \cdots, i-1 \\ {\rm Var(}{{Z}_{i}}{\rm )=}\underset{{{{{a}'}}_{i}}{{a}_{i}}=1, {{{{a}'}}_{i}}\sum {{a}_{j}}=0(j=1, \cdots, i-1)}{\mathop{{\rm max}}}\, {\rm Var(}{a}'X{\rm )} \\\end{array} \right. $ | (6) |
则称
$ R=\dfrac{1}{n-1}{X}'X $ | (7) |
由(3)式可知,求取主成分就是一个条件极值问题,可应用拉格朗日乘子法求解,得出
$ \left\{ \begin{array}{l} Z=XA(i=1, 2, \cdots , p) \\ { D(}R{\rm )=diag}({{\lambda }_{1}}, {{\lambda }_{2}}, \cdots , {{\lambda }_{p}}) \\ \end{array} \right . $ | (8) |
主成分分析的目的是简化数据结构,故在实际应用中一般不会用
由于有效干层的识别属于非线性预测问题,加之地质条件复杂,因此采用模糊数学的方法来对干层有效性进行识别可以提高预测精度[14-20]。采用主成分分析法求解出对判别影响较重要的几个因素作为模糊贴近度的输入变量。假设对给定的
首先,经过模糊运算得出决策对象
贴近度决策法是通过
$\left\{ \begin{array}{l} \sigma{:F(U)}\times {F(U)}\!\to \! \left[ {\rm 0, 1} \right]{ (\underset{\scriptscriptstyle\thicksim}{C}\!, \!\underset{\scriptscriptstyle\thicksim}{D})}\mapsto \sigma{ (\underset{\scriptscriptstyle\thicksim}{C}, \underset{\scriptscriptstyle\thicksim}{D})} \\ \sigma{ (\underset{\scriptscriptstyle\thicksim}{C}, \underset{\scriptscriptstyle\thicksim}{C})=1;} \\ {( }\varPhi{, \bigcup)=0;} \\ \sigma{ (\underset{\scriptscriptstyle\thicksim}{C}, \underset{\scriptscriptstyle\thicksim}{D})= }\sigma{ (\underset{\scriptscriptstyle\thicksim}{D}, \underset{\scriptscriptstyle\thicksim}{C})} \\ {\underset{\scriptscriptstyle\thicksim}{C}}\subset { \underset{\scriptscriptstyle\thicksim}{D}}\subset { \underset{\scriptscriptstyle\thicksim}{E}}\Rightarrow \sigma{ (\underset{\scriptscriptstyle\thicksim}{C}\!, \!\underset{\scriptscriptstyle\thicksim}{E})}\leqslant \sigma{ (\underset{\scriptscriptstyle\thicksim}{C}\!, \!\underset{\scriptscriptstyle\thicksim}{D})}\wedge \sigma{ (\underset{\scriptscriptstyle\thicksim}{C}\!, \!\underset{\scriptscriptstyle\thicksim}{E})} \\ \end{array} \right. $ | (9) |
采用非对称模糊称贴近度法
$ \sigma (\underset{\scriptscriptstyle\thicksim}{C}, \underset{\scriptscriptstyle\thicksim}{D})=1-\dfrac{1}{n}\sum\limits_{i=1}^{n}{\left| {{\mu }_{{\underset{\scriptscriptstyle\thicksim}{C}}}}({{u}_{i}})-{{\mu }_{{\underset{\scriptscriptstyle\thicksim}{D}}}}({{u}_{i}}) \right|} $ | (10) |
其中,
根据冀东油田G5断块的实际生产数据,采用主成分分析方法及模糊贴近度方法对可利用干层进行了判别,首先采用孔隙度、渗透率、有效厚度、孔喉半径、排驱压力作为判别参数,统计某一口井在测井解释中解释为干层的小层,具体参数值见表 1。
首先,对判别参数进行标准化处理求出各参数的均值及方差;然后,应用式(1)
将数据代入,求出各类参数的特征向量及特征值为42.77,13.38,4.16,0.47,0.08,得到累计贡献率为0.700 0,0.920 0,0.990 0,0.999 8,1.000 0。
一般主成分分析取累计贡献率达到70%以上的数据,为了达到更加精确的目的,选取累计贡献率达到90%以上的数据,因此选取
由
$ \left\{ \begin{array}{l} {Z_1} = - 0.4922{X_1} - 0.5116{X_2} - 0.2746{X_3} - \\{\kern 40pt}0.5136{X_4} + 0.3961{X_5} \\ {Z_2} = - 0.3619{X_1} + 0.187{X_2} + 0.7547{X_3} + \\{\kern 40pt}0.1389{X_4} + 0.4952{X_5} \end{array} \right. $ | (11) |
优选的影响因素见表 3。
选取影响因子Z的隶属度函数为
$ \mu ({Z_i}) = \left\{ {\begin{array}{l} {1 - \dfrac{1}{2}{{\rm{e}}^{ - {\left(\frac{{{Z_{ij}} - {{\bar Z}_i}}}{{{S_i}}}\right)^2}}}, {\kern 5pt}{Z_{ij}} \leqslant {{\bar Z}_i}}\\[6pt] {{{\rm{e}}^{ - {\left(\frac{{{Z_{ij}} - {{\bar Z}_i}}}{{{S_i}}}\right)^2}}}, {\kern 30pt} {{\bar Z}_i} + {S_i} > {Z_{ij}} \geqslant {{\bar Z}_i} }\\ 0, {\kern 60pt} {Z_{ij}} > {{\bar Z}_i} + {S_i} \end{array}} \right. $ | (12) |
其中,
将
由
当
判别发现,15、25、36小层为没有射孔的有效干层,因此,可对15、25、36号小层进行注水对应油井射孔生产,以提高注采井的水驱动用程度。应用数值模拟进行计算,与不利用有效干层对比,利用有效干层的方案整体采出程度提高了约2
以25号小层剩余油分析结果为例,图 4为解释为有干层的水井没有射孔的波及面积图,图 5为解释为有效干层的水井打开射孔后充分利用有效干层的水驱油波及面积图,从图 4和图 5对比可以看出,图 5的水驱油波及面积明显大于图 4的水驱油波及面积,表明充分利用有效干层对剩余油的挖潜有一定的作用,能够增大水驱油面积,提高水驱效率(图 4,图 5)。
(1) 渤海湾中深层油藏纵向上砂泥岩频繁交互,厚、薄砂体连通方式复杂,以及单成因砂体自身的岩性和物性的空间差异导致了储层非均质性较强,注采系统难以完善,利用模糊隶属度判别函数的识别方法,能够识别出有效干层,以便加以利用。
(2) 通过对识别的出有效干层进行射孔、补孔等措施进而完善注采关系,提高水驱波及面积,能够有效的提高水驱动用程度。
(3) 通过对冀东油田高5断块干层有效性的验证,应用本文计算结果与实际吸水剖面测试数据基本一致,符合率达到85%以上。
[1] |
胡文瑞. 论老油田实施二次开发工程的必要性与可行性[J]. 石油勘探与开发, 2008, 35(1): 1-5. HU Wenrui. Necessity and feasibility of petrochina mature field redevelopment[J]. Petroleum Exploration and Development, 2008, 35(1): 1-5. doi: 10.3321/j.issn:1000-0747.-2008.01.001 |
[2] |
胡文瑞. 中国石油二次开发技术综述[J]. 特种油气藏, 2007, 14(6): 1-4. HU Wenrui. The summarize of the PetroChina mature field redevelopment technologies[J]. Special Oil and Gas Reservoir, 2007, 14(6): 1-4. doi: 10.3969/j.issn.1006-6535.2007.06.001 |
[3] |
韩大匡. 关于高含水油田二次开发理念、对策和技术路线的探讨[J]. 石油勘探与开发, 2010, 37(5): 583-591. HAN Dakuang. Discussions on concepts, countermeasures and technical routes for the redevelopment of high water-cut oilfields[J]. Petroleum Exploration and Development, 2010, 37(5): 583-591. |
[4] |
冉启佐. 高尚堡深层北区沙三段二、三亚段沉积体系与沉积相[J]. 石油与天然气地质, 1997, 18(1): 13-16. RAN Qizuo. Gao Shangpu deep and north Es32+3 sedimentary system and sedimentary facies[J]. Oil & Gas Geology, 1997, 18(1): 13-16. |
[5] |
叶庆全. 油气田开发地质[M]. 北京: 石油工业出版社, 1999. YE Qingquan. Development of oil and gas geology[M]. Beijing: Petroleum Industry Press, 1999. |
[6] |
穆龙新. 储层井细研究方法[M]. 北京: 石油工业出版社, 2000. MU Longxin. Reservoir sophisticated research methods[M]. Beijing: Petroleum Industry Press, 2000. |
[7] |
陈程, 孙义梅, 邓宏文. 油田开发后期扇三角洲前缘微相分析及应用[J]. 现代地质, 2001, 15(1): 88-93. CHEN Cheng, SUN Yimei, DENG Hongwen. Microfacies analysis of fan delta front and its application in Shuanghe Oil Field[J]. Geoscience, 2001, 15(1): 88-93. doi: 10.3969/j.issn.1000-8527.2001.01.016 |
[8] |
付国民, 李永军, 石京平. 特高含水期扇三角洲储集层剩余油分布及挖潜途径[J]. 成都理工大学学报(自然科学版), 2003, 30(2): 178-183. FU Guomin, LI Yongjun, SHI Jingping. Residual oil distribution and the approach of tapping the latent power of the fan delta reservoir during the high water cut stage[J]. Journal Of Chengdu University Of Technology (Science & Technology Edition), 2003, 30(2): 178-183. doi: 10.3969/j.issn.1671-9727.2003.02.010 |
[9] |
王友净. 高尚堡深层北区沙三段二、三亚段沉积特征与砂体结构[J]. 油气地质与采收率, 2010, 17(2): 14-16. WANG Youjing. Gao Shangpu deep and north Es32+3 sedimentary characteristics and sand body structure[J]. Petroleum Geology and Recovery Efficiency, 2010, 17(2): 14-16. doi: 10.3969/j.issn.1009-9603.2010.02.004 |
[10] |
秦同洛. 实用油藏工程方法[M]. 北京: 石油工业出版社, 1989. QIN Tongluo. Practical reservoir engineering method[M]. Beijing: Petroleum Industry Press, 1989. |
[11] |
孙国. 胜坨油田特高含水期井网重组技术优化研究[J]. 油气地质与采收率, 2005, 12(3): 48-50. SUN Guo. The well pattern reorganization of Sheng Tuo Oilfield in extra-high water-cut period[J]. Petroleum Geology and Recovery Efficiency, 2005, 12(3): 48-50. doi: 10.3969/j.issn.1009-9603.2005.03.015 |
[12] |
辛治国. 应用模糊综合评判法评价油层水淹状况[J]. 油气地质与采收率, 2007, 14(6): 88-90. XIN Zhiguo. Evaluation of watered out behavior of reservoirs using fuzzy comprehensive assessment[J]. Petroleum Geology and Recovery Efficiency, 2007, 14(6): 88-90. doi: 10.3969/j.issn.1009-9603.2007.06.027 |
[13] |
何琰. 基于模糊综合评判与层次分析的储层定量评价——以包界地区须家河组为例[J]. 油气地质与采收率, 2011, 18(1): 23-25, 29. HE Yan. A quantitative reservoir evaluation method based on fuzzy comprehensive appraisal and analytical hierarchy process -Case of Xujiahe Formation, Baojie Area[J]. Petroleum Geology and Recovery Efficiency, 2011, 18(1): 23-25, 29. doi: 10.3969/j.issn.1009-9603.2011.01.007 |
[14] |
刘正锋, 燕军. 模糊识别方法在储层识别中的应用[J]. 西南石油学院学报, 1998, 20(3): 4-6. LIU Zhengfeng, YAN Jun. The Application of fuzzy technique in reservoir identification[J]. Journal of Southwest Petroleum Institute, 1998, 20(3): 4-6. doi: 10.3863/j.issn.-1674-5086.1998.03.002 |
[15] |
冯国庆, 李允, 谈德辉. 模糊贴近度在储层识别中的应用[J]. 西南石油学院学报, 1999, 21(4): 46-49. FENG Guoqing, LI Yun, TAN Dehui. The application of the fuzzy degree in reservoir identification[J]. Journal of Southwest Petroleum Institute, 1999, 21(4): 46-49. doi: 10.3863/j.issn.1674-5086.1999.04.013 |
[16] |
刘丽. 基于密闭取心井和模糊评判分析埕岛油田水淹状况[J]. 西北大学学报(自然科学版), 2012, 42(6): 995-999. LIU Li. Water out behavior analysis of Chengdao Oilfield by sealing core drilling and fuzzy comprehensive evaluation[J]. Journal of Northwest University (Natural Science Edition), 2012, 42(6): 995-999. doi: 10.3969/j.issn.1000-274X.2012.06.025 |
[17] |
夏崇双, 刘林清, 张理, 等. 四川盆地老气田二次开发优化技术及应用[J]. 天然气工业, 2016, 36(9): 80-89. XIA Chongshuang, LIU Linqing, ZHANG Li, et al. Optimization for the secondary development of old gas fields in the Sichuan Basin and its application[J]. Natural Gas Industry, 2016, 36(9): 80-89. doi: 10.3787/j.issn.1000——0976.2016.09.009 |
[18] |
查娜. 大庆油田二次开发技术策略及实践[J]. 科技和产业, 2013, 13(4): 151-154. CHA Na. The technical policy and practice of the secondary development in Daqing Oilfield[J]. Science Technology and Industry, 2013, 13(4): 151-154. doi: 10.-3969/j.issn.1671-1807.2013.04.035 |
[19] |
何江川. 中国石油二次开发技术与实践(20082010年)[M]. 北京: 石油工业出版社, 2012. HE Jiangchuan. The second development technology and practice of china petroleum(20082010)[M]. Beijing: Petroleum Industry Press, 2012. |
[20] |
马贵阳, 宫清君, 潘振, 等. 基于支持向量机结合遗传算法的天然气水合物相平衡研究[J]. 天然气工业, 2017, 37(5): 46-52. MA Guiyang, GONG Qingjun, PAN Zhen, et al. GA-SVM based study on natural gas hydrate phase equilibrium[J]. Natural Gas Industry, 2017, 37(5): 46-52. doi: 10.3787/j.-issn.1000-0976.2017.05.006 |