商建霞,张乔良,李华,孔令辉,罗佼,叶苗
中海石油(中国)有限公司湛江分公司,广东湛江 524057
南海西部海域大部分在生产油田已进入“高含水、高采出程度”的双高阶段,加之地下储层非均质性及开采不均匀性,剩余油分布十分复杂,挖潜难度愈来愈大。储层内部不同层次的复杂非均质性是导致高含水期油藏剩余油高度分散的主要地质因素,因此,必须系统地厘清储层的非均质体系,描述砂体内部结构并有效评价其渗流品质,为改善高含水期油藏开发效果提供有效的技术支持。
当前,储层流动单元刻画是描述油藏非均质性及评价储层渗流性能行之有效的手段之一,它是指在垂向上和横向上连续的储集带且内部具有相似的岩性及岩石物理性质[1]。自1984年Hearnet al.[1]首次提出流动单元的概念以来,国内外学者提出了多种流动单元划分的方法:冯晓宏等[2]、姚光庆等[3]、吕晓光等[4]提出利用流动层段指标FZI方法,吴胜和等[5]提出了以地质研究为主的储层层次分析法,彭仕宓等[6]、姚合法等[7]提出了采用粒度中值、孔隙度及流动层段指标FZI方法等多参数判别法,每种方法都有自身的适用条件和优缺点,国内目前比较常用的是利用FZI的方式进行流动单元划分,但是该方法主要反映岩石的孔隙度和渗透率,据此划分精度不够,划分结果与生产动态吻合度不高,无法满足高含水油田精细研究的需要,本次以研究区珠江口盆地C油田珠海组二段一油组(简称“ZH2Ⅰ油组”)为例,在构型研究基础上,结合毛管压力曲线等微观孔喉参数提出了一种全新的流动单元划分方法——基于构型的渗流单元刻画新技术,有效评价了储层的渗流品质,并据此将研究区划分60个渗流单元,提高了剩余油预测精度,为高含水油田调整挖潜提供关键的技术支持。
文昌C油田位于南海珠江口盆地西部珠三坳陷的文昌B凹陷中部东侧(图1)。油田所在的珠三坳陷发育了完整的新生界地层。油田内所钻遇油层属新近系中新统韩江组、珠江组和古近系渐新统珠海组、恩平组地层,其中ZH2Ⅰ油组是油田主力油组之一。珠海组含油范围内断层较为发育,主断层断距大,延伸较长,且次一级小断层相对发育,对珠海组的油水分布具有控制作用。ZH2Ⅰ油组属于扇三角洲前缘沉积,扇三角洲前缘亚相,以水下分流河道微相为主,油藏类型属于受断层遮挡的边水油藏。由于纵向多期砂体相互叠置,内部结构复杂、非均质性强,经过近十年的有效开发,目前已进入高含水、高采出程度的“双高”开发阶段,剩余油分布十分复杂。
图1 珠江口盆地西部构造区划图及研究区位置[8]Fig.1 Geotectonic map of the western Pearl River Mouth Ba⁃sin and location of the study area[8]
储层构型研究是渗流单元研究的基础,通过开展构型研究可以明确连通体及渗流屏障的分布关系,进而在连通体内部开展渗流单元的划分,构型单元界面级次与储层渗流屏障的分级也有较好的对应关系[9⁃10]。
研究区ZH2Ⅰ油组为扇三角洲前缘沉积,构型单元包括前缘水下分流河道主体、分下分流河道侧缘、河口坝、席状砂以及分流河道间5种,沉积微相与储层构型具有成因联系。物源来自东部的神狐隆起,受物源供给及古地形影响,水下分流河道砂体主要在研究区中东部连片分布,沿物源方向过渡为河口坝及席状砂,沉积微相控制储层构型平面分布特征(图2)。纵向剖面上水下分流河道发育区由于多期河道侧向迁移,呈相互切割叠置接触,砂体之间连通性好;远离物源方向席状砂砂体处呈孤立分布特点,剖面沉积微相控制储层构型剖面分布特征(图3)。
图2 研究区ZH2Ⅰ⁃2⁃1⁃2单层构型平面分布Fig.2 Planar distribution of two small⁃layer/one single⁃layer configuration of ZH2 I oil group
图3 研究区ZH2Ⅰ⁃2⁃1⁃2单层构型剖面分布Fig.3 Vertical distribution of two small⁃layer/one single⁃layer configuration of ZH2 I oil group
依据Miall[11⁃13]提出的界面划分原则,结合用岩心、地震、测井等资料,采用井震结合、相控约束、模式指导的方法[14⁃17],将目的层构型界面按正序由大到小从7级划分到3级构型界面(表1、图4)。
图4 研究区ZH2Ⅰ油组构型模式Fig.4 Area configurations of the ZH2 I oil group in the study area
表1 ZH2Ⅰ油组构型界面划分Table 1 Structural interface division of ZH2 I oil group
7级构型界面主要为扇三角洲前缘复合朵叶体之间的界面,对应7级渗流屏障界面,主要发育厚层泥岩,厚度在2 m左右,分布相对稳定,具有较强的阻隔作用,该级次界面对应的是油层组界面。6级构型界面为扇三角洲前缘不同朵叶体之间界面,对应6级渗流屏障界面,成因类型主要为几个复合水下分流河道顶部的泥质隔层,厚度大于1.5 m。从地震反射特征上,可以看出垂直物源方向三期朵叶砂体形态均呈顶平底凸,自下而上,砂体延伸范围变小,表明随着水体的加深,扇三角洲发育规模缩小,该级次界面对应的是砂层组界面。5级构型界面为多期水下分流河道复合体的顶界面,对应5级渗流屏障界面,成因类型为水下分流河道叠置体间泥岩,厚度约1.5 m左右,与小层组的界面相对应。4级构型界面为单一水下分流河道砂体的顶界面,在单井上表现为泥岩夹层或者较大程度的曲线回返面,对应4级渗流屏障界面,在岩心上对应厚度较薄的砂泥过渡段,物性差,且延伸短,厚度大于1 m。3级界面为单一河道内部增生体的分界面,成因类型有两种,一是河道内顶部静水沉积形成泥质夹层,二是底部滞留沉积形成砂砾混杂物性差的夹层[18⁃19],对应3级渗流屏障界面,厚度较薄,大多小于1 m。
在构型界面划分的基础上,进行构型单元的解剖,结果表明ZH2Ⅰ油组自下而上由于水下分流河道的不断迁移摆动,形成了纵向上多期河道相互切割叠置,平面上连片分布的空间格局。油组内部进一步分为7个复合砂体、26个单砂体,单砂体宽度100~450 m,厚度1.5~6 m,宽厚比为40~80。
流动单元与渗流单元既有联系又有区别,实际上,它们都是不同尺度储层非均质单元,主要用于油田不同开发阶段;二者的区别在于:前者以孔隙结构、孔隙度和渗透率等静态特征为主,结合流体流动的动态特征进行划分;而后者则以储层内部动态渗流性能为主,强调渗流动用状态,结合储层构型静态特征进行划分。由于研究区已进入开发中后期,因此本次在储层构型划分的基础上,在最小构型单元内部,开展了基于构型的渗流单元的研究。
(1)参数优选
传统的流动单元划分往往采用流动分层指标(FZI)进行划分。理论上,FZI是把结构和矿物地质特征、孔喉特征结合起来判定孔隙几何相的一个参数,实际上FZI计算方法过于理想化,弱化了孔隙结构非均质性的影响,实际上决定因素仅有孔隙度和渗透率两个参数,根据FZI划分流动单元精度不够,在同一单元内部仍然存在水驱不均的现象,因此,仅用FZI一个参数划分流动单元的方法无法满足开发后期对渗流单元划分精细的需要[20⁃22]。考虑到C油田层内非均质性强,流动单元划分难度较大,仅仅依靠孔渗参数不能反映流动单元的全部特征,因此,本次在充分考虑研究区储层地质特征、流体特征基础上,综合考虑构型成果(构型要素、渗流屏障、砂体叠置关系等)、岩性参数(岩性、粒度中值、泥质含量)、岩心物性参数(孔隙度、渗透率)、孔喉参数(毛管压力曲线、启动压力、中值压力、均质半径)四大类评价参数。
(2)取心井渗流单元分类
依据四类评价参数将目的层分为优质、次优、差三类渗流单元,具体的划分标准见图5所示。
图5 ZH2Ⅰ油组渗流单元分类依据表Fig.5 Classifications of seepage unit of ZH2 I oil formation
①优质渗流单元此类构型要素主要为水下分流河道主体,砂体叠置样式为叠置式、渗流屏障发育特点以无夹层式为主。岩性以含砾中—细砂岩为主,粒度中值为9.92µm,平均泥质含量为2.26%。岩心孔隙度平均26%,岩心渗透率平均为1 560×10-3µm2。毛管压力曲线以I类为主,启动压力为0.003 MPa,中值压力小于0.1 MPa,均质半径大于20µm。
②次优渗流单元此类构型要素为水下分流河道侧缘、河口坝,砂体叠置样式为对接式为主、渗流屏障发育特点以半遮挡型夹层为主。岩性以中—细砂岩为主,粒度中值为1.03µm,平均泥质含量为9.1%。岩心孔隙度平均23%,岩心渗透率平均为564×10-3µm2。毛管压力曲线以II类为主,启动压力为0.003~0.007 MPa,中值压力0.1~0.15 MPa,均质半径大于15~20µm。
③差渗流单元此类构型要素为席状砂、溢岸砂,砂体叠置样式以孤立式为主,渗流屏障发育特点以全遮挡型夹层为主。岩性以粉砂岩为主,粒度中值为0.05µm,平均泥质含量为10.46%。岩心孔隙度平均19%,岩心渗透率平均为115×10-3µm2。毛管压力曲线以III类为主,启动压力为0.07~0.02 MPa,中值压力0.2~0.5 MPa,均质半径大于5~15µm。
(3)取心井渗流单元刻画
结合上述分类标准对取心井C8井进行渗流单元的刻画,建立取心井不同类别渗流单元的电性参数区间(表2),比传统方法更加精细(图6),与生产动态更加匹配,优质渗流单元水淹严重,而次优、差的渗流单元弱水淹—未水淹。
表2 取心井渗流单元的电性参数分布Table 2 Distribution of electrical parameters of percolation unit in coring well
假设G个总体t1,t2,…,tg是由m个变量a1,a2,…,am组成,假设已知每个总体的密度函数为fg(a)和先验概率为Pg。如果把t1,t2,…,tg视为总体样本空间的一个划分,依据Bayes公式可以求得样品a属于tg(t=1,2,…,G)的条件概率:
式中:P(ti/a)为a的条件概率;Pg为每个总体的先验概率;fg(a)为每个样品的密度函数;Pi为第i个变量概率值;fi(a)为第i个变量a的函数。
研究区共划分优质、次优、差三类渗流单元,可视为三个总体G=1,2,3。选取了六条测井曲线为变量,经过标准化之后对非取心井进行判别分析,建立判别函数如下:
式中:GR为伽马值(API);CNC为中子值(%);DEN为密度值(g/cm3);VSH为泥质含量(%);Por为孔隙度值(%);Perm为渗透率(10-3µm2)。
依据判别函数划分单井渗流单元划分结果与动态监测资料匹配分析结果表明:渗流单元划分结果与实际水淹情况具有较好的相关性(图7),优质和次优渗流单元产出高(17.1%~28.4%)、含水率高(82%~97.6%),而次优、差储层产出少(4.7%~10.3%),含水率低(47.6%)。
图7 渗流单元划分与动态监测对比图Fig.7 Comparison of seepage unit division and dynamic monitoring
在储层构型单元内部,依据由单井渗流单元划分—平面渗流单元划分—剖面渗流单元划分的多维互动研究思路,完成对渗流单元分布特征研究,结果表明,横向上,靠近物源区,以优质和次优储层为主,远离物源区,次优和差储层比例增加,优势渗流单元主要分布于水下分流河道主体处(图8)。纵向上,席状砂、水下分流河道侧缘以及水下分流河道顶部及底部次优、差储层为主,水下分流河道主体中部以优质储层为主,从而说明了构型单元的展布对渗流单元分布具有一定的控制作用(图9)。
图8 研究区ZH2Ⅰ⁃2⁃1⁃2单层渗流单元平面分布Fig.8 Planar seepage cell division distribution of two small⁃layer/one single⁃layer configuration of ZH2 I oil group
图9 研究区ZH2Ⅰ⁃2⁃1⁃2单层渗流单元剖面分布Fig.9 Vertical seepage cell division distribution of two small⁃layer/one single⁃layer configuration of ZH2 I oil group
结合构型成果及渗流单元划分,对整个ZH2Ⅰ油组进行了渗流单元的划分,纵向上15个单层、平面5个井区,一共划分75个单元,其中位于油水界面以上含油的是60个单元。其中优质单元29个,储量占55.6%;次优单元20个,储量占30.2%;差单元11个,储量占14.2%。
根据动态监测结果,将产量劈分到所在的单元,得到各单元的采出程度(图10),采出程度相对较低的渗流单元即剩余油富集区,通过分析总结控油模式有两种类型:渗流单元差异平面控油模式、渗流屏障遮挡垂向控油模式。
图10 不同渗流单元采出程度三维立体柱状图Fig.10 Three⁃dimensional histogram of recovery percent of different seepage units
(1)渗流单元差异平面控油模式
由于看似连片的砂体内部存在渗流差异性,边水在驱替过程中会优先选择物性好的优质、次优渗流单元推进,而物性差的差渗流单元储层难以波及,水驱程度弱,剩余油相对富集。研究区内优质、次优渗流单元砂体成因类型多为水下分流河道、河口坝砂体,储层物性好,内部往往水淹严重,剩余油潜力小;而差渗流单元例如席状砂、溢岸砂等砂体,水驱程度弱,剩余油富集,比如ZH2Ⅰ⁃1⁃1⁃1单层内部5个席状砂成因的差渗流单元采出程度低,剩余潜力大。
(2)渗流屏障遮挡垂向控油模式
渗流屏障也即夹层,对油水垂向运移起到一定的阻挡作用。垂向上,在渗流单元间或渗流单元内部夹层发育时,往往造成界面上下储层水淹程度不同,导致夹层下部单元或该单元内部剩余油富集。例如研究区四级构型单元河道内部存在两种夹层类型,分别是顶部静水细粒沉积物与底部滞留沉积物,若该类型夹层发育会导致河道内部水洗程度较轻;流动单元间隔层的存在同样对流体垂向运移起到明显的阻隔作用,在研究区分布均比较稳定夹层主要为三级或四级界面,厚度大于1 m,在研究区分布均比较稳定,造成界面上下储层水驱程度的不均匀,剩余油富集。该类型剩余油主要分布在B3井区、B5井区ZH2Ⅰ⁃2⁃2~ZH2Ⅰ⁃3⁃2共10个单元。
针对渗流差异控制的剩余油富集区,由于储层物性相对较差,本次暂不调整;受渗流屏障控制的剩余油富集区,其储层品质主要为优质、次优,潜力相对较大,因此作为本次挖潜主要目标,分别实施调整井、补孔共七项措施,平均单井增油3.7万方,累积增油26万方,提高采收率5%。
(1)流动单元与渗流单元属于不同尺度储层非均质单元,适用于不同开发阶段。前者以储层静态特征为主,结合流体流动的动态特征进行划分;而后者则以储层内部动态渗流性能为主,结合储层构型静态特征进行划分。文昌C油田ZH2Ⅰ油组的研究成果应用表明,油田进入开发中后期之后,在成因单元内部开展渗流单元研究能有效指导剩余油分布预测以及油田开发调整策略。
(2)依据储层构型理论,将研究区构型界面按正序由大到小从7级划分至3级构型界面,7级界面主要为扇三角洲前缘复合朵叶体之间界面,6级界面为不同朵叶体之间界面,5级界面为水下分流河道复合体的顶界面,4级界面为单一水下分流河道砂体的顶界面,3级界面为单一河道内部增生体的分界面,其中4、5级界面对开发有较大影响。在沉积微相认知基础上,识别出5种构型单元。
(3)在构型单元内部,选取构型成果、岩性参数、岩心物性参数、孔喉参数四大类评价参数,将取心井储层划分为优质、次优、差三类渗流单元,并运用神经网络的方式运用到非取心井进行渗流单元的划分。优质渗流单元主要分布在水下分流河道主体,次优渗流单元主要分布于水下分流河道侧缘、河口坝处,而差渗流单元主要集中于水下分流河道顶部漫溢细粒沉积和底部滞留沉积段、席状砂、溢岸砂处。
(4)受构型及渗流品质的影响,剩余油富集模式有两类:渗流单元差异平面控油模式、渗流屏障遮挡垂向控油模式。渗流单元与剩余油分布具有较好的相关关系,无渗流屏障遮挡的优质渗流单元水淹严重,次优渗流单元、差渗流单元、受渗流屏障遮挡的下部优质单元水淹程度弱,剩余油富集,本次调整挖潜选取受渗流屏障控制的优质、次优剩余油富集单元进行挖潜,累积增油26万方,提高采收率5%。