史建忠 才巨宏 张玲 杨英珍 田莹摘要 碳酸盐岩、花岗片麻岩潜山油藏的主要特点是岩性复杂、储集空间类型多、非均质性强,储集层评价比较困难。文章以埕岛油田埕北30潜山为例,对该类储集层定量解释中的几个关键参数进行了深入细致的探讨,初步形成一套适用于复杂岩性潜山油藏的储集层测井评价技术。关键词 济阳坳陷 埕岛油田 潜山 碳酸盐岩 变质岩 储集层 测井评价一、引言埕岛油田埕北30块油藏类型为潜山内幕型易挥发轻质油藏,储集层为古生界的碳酸盐岩和太古宇的花岗片麻岩,储集空间类型有裂缝、溶蚀孔洞、晶簇孔等多种,其基质也具有一定的储集能力。由于岩性复杂、储集空间类型多、储集层非均质强,给测井定量评价工作带来了很大困难。为此,结合该区实际情况,借助新的测井方法、新的测井解释软件,初步建立了一套基于复杂岩性储集层的测井定量评价方法,实际应用效果比较显著。二、孔隙度解释技术总孔隙度埕北30潜山原生孔隙不发育,对油气富集高产起决定作用的是次生孔隙,具有缝、洞、孔三大类,另外,荧光分析发现,其基质也具有一定的储油能力。埕北30潜山孔隙度解释主要使用由 Schumberger公司引进的Petrophysics软件包进行解释,核心程序是ELAN,解释时需结合岩心分析和核磁共振测井资料。ELAN软件的基本思路是以实际测井值为基础,根据地层矿物组分建立合适的解释模型和测井响应方程,通过合理选择解释参数,反算相应的理论测井值,并与实际测井值比较,按非线性加权最小二乘法原理建立目标函数,不断调整未知储集层参数,使目标函数达到极小值。其优点是充分利用所有测井信息,采用最优化技术使解释结果最为合理。单井处理过程包括填写参数卡、初步解释、解释结果与岩心分析对比、修改参数卡、再解释等五个步骤。在对各井进行处理时,首先根据测井曲线及地区地质经验,填写参数卡进行初步解释,然后将解释结果与岩心分析进行对比,如果处理井段没有岩心分析数据,则根据反算的理论曲线和实测曲线的拟合情况适当修改参数卡,直到与岩心分析数据吻合或理论曲线与实测曲线拟合较好为止。为了更好地利用好井眼段的核磁共振测井资料,做了好井眼井段的核磁孔隙度与声波、密度、中子三种测井视孔隙度的关系研究,发现相互对应关系均比较好(图1)。图1 埕北302井古生界核磁孔隙度与补偿中子关系图视孔隙度求取公式为:胜利油区勘探开发论文集式中:ΦD——视密度孔隙度,%;ΦS——视声波孔隙度,%;ρb——补偿密度测井值,g/cm3;t——补偿声波测井值,μs/m。那么,在好井眼段用核磁孔隙度;在坏井眼段核磁共振测井失真,密度未受影响时用视密度孔隙度与核磁孔隙度的关系求得核磁共振孔隙度;如果核磁共振、密度测井都不可信,则用视声波孔隙度与核磁共振孔隙度的关系求取核磁共振孔隙度,得出一条综合的“核磁共振孔隙度”曲线,将其作为一条输入曲线,参加ELAN的反演。这样,既利用了核磁孔隙度的准确性,又充分利用了其他测井曲线,提高了ELAN解释结果的准确性和可靠性。为检验解释结果的可靠性,进行了精度分析。从图2中可看出,在致密段,因为岩心分析代表的是总孔隙度,测井解释与岩心分析吻合较好;在储集层段,由于缝、洞的存在,测井解释孔隙度大于岩心分析孔隙度,也是比较合理的。裂缝孔隙度埕北30潜山油藏具有双重孔隙结构特征,油田开发中裂缝孔隙度是一个重要参数。根据专业文献资料,裂缝孔隙度一般不超过1%,考虑有与裂缝连通的溶洞的存在,包括缝洞的裂缝系统孔隙通常低于2%。裂缝孔隙度通常根据双侧向测井资料求得,AMSibbit和QFaivre提出的利用双侧向电阻率计算裂缝孔隙度公式为:[1]图2 埕北303井太古宇测井解释与岩心分析孔隙度交会图油气层胜利油区勘探开发论文集水层胜利油区勘探开发论文集式中:mf——裂缝孔隙度指数;Rm——泥浆电阻率,Ω·m;Rth——岩块电阻率,Ω·m;Rlls——浅侧向电阻率,Ω·m;Kr——双侧向畸变系数,低角度缝取2,斜交缝取1,垂直缝取0;Rw——地层水电阻率,Ω·m。对于进行了岩心分析的井段,可以认为岩心分析为岩块系统孔隙度,测井解释为总孔隙度,用测井解释孔隙度减去岩心分析孔隙度后可得该井段裂缝孔隙度,然后以此对mf和Kr进行刻度,也可根据成像资料或录井资料确定mf和Kr的值。埕北30潜山解释4口井,平均裂缝孔隙度为44%,其中埕北303井解释裂缝孔隙度为15%,岩块孔隙度为75%,这与试井解释的裂缝孔隙度2%、岩块孔隙度8%对应较好,说明裂缝孔隙度解释比较可靠,用岩心刻度法求取裂缝参数是可行的。三、渗透率解释技术在双重孔隙结构的裂缝性地层中,渗透率为岩块渗透率和裂缝渗透率的综合反映,由于岩块系统渗透率非常低,大都小于1×10-3μm2,因此储集层渗透率主要为裂缝渗透率的反映。经验建模法根据该区全直径岩心分析资料,建立了孔隙度和渗透率的经验关系模型(图3),由于全直径分析样品比较少,这种方法计算的渗透率代表性较差,仅供参考。图3 埕北30潜山全直径岩心分析 孔隙度、渗透率关系图核磁共振测井解释核磁共振测井解释渗透率为:胜利油区勘探开发论文集式中:k——渗透率,10-3μm2;Φnmr——核磁测井有效孔隙度,小数;T2g——T2几何平均值,ms;C、m、n——经验系数。根据埕北302井古生界、太古宇 14块岩心样品的核磁测试数据,对上式中的经验系数进行刻度,古生界 6块岩样的C、m、n平均值分别为639、711、531,太古宇8块岩样的C、m、n平均值分别为632、524、089,T2g根据核磁测试古生界、太古宇平均值分别为88ms和71ms,用(5)式分别对埕北302、303两口井好井眼段进行了解释。由于公式中的各项参数均经过岩心刻度,且好井眼段核磁测量孔隙度是可靠的,用该式解释的渗透率基本代表井眼的实际情况。ELAN软件解释ELAN软件解释渗透率为一种地球化学算法,公式如下:胜利油区勘探开发论文集式中:Φt——总孔隙度,小数;胜利油区勘探开发论文集N——地层中矿物总数;Fi——第i种矿物的渗透率因子;Wi——第 i种矿物的重量百分比,%。这种算法既考虑了总孔隙度,又考虑了各种矿物组分及其百分含量,是一种比较合理的渗透率解释方法。试井解释该区对埕北301、302、303井进行了试井,并用试井解释软件进行解释,利用压力恢复典型曲线拟合分析解释了渗透率各参数(表1)。表1 试井解释成果表对比上述四种方法解释结果(表2),经验公式法、ELAN、试井解释三种方法解释结果比较接近,再将ELAN和核磁共振解释结果进行了对比(图4),两种解释结果吻合较好。结合地质、油藏方面的研究成果,认为经验公式、核磁、ELAN、试井等四种方法解释结果符合地下实际情况,由于经验公式及试井解释的局限性,最终结果以ELAN和核磁共振解释为准。表2 渗透率解释对比表四、含油饱和度解释技术针对埕北30潜山没有进行密闭取心和油基泥浆取心分析,以及裂缝性油藏其裂缝的发育程度和分布是多变的,采用以下几种方法获取油藏原始含油饱和度资料。图4 埕北302井古生界 ELAN解释渗透率与核磁共振解释渗透率交会图阿尔奇方程阿尔奇方程是建立在均匀孔隙基础上的饱和度解释方程[2],即:胜利油区勘探开发论文集式中:Sw——含水饱和度,小数;Φ——孔隙度,小数;m——胶结指数;Rw——饱水电阻率,Ω·m;Rt——岩块电阻率,Ω·m;n——饱和度指数;a、b——岩电系数,一般取1。由于该区无法做岩电实验分析,式中 m、n等参数均根据理论值选取,m=n=2,a=b=1。在裂缝性地层中,泥浆侵入深度大,而且侵入深度的变化范围也很大,求得的饱和度值是在侵入带至原状地层之间变化。对于缝、洞不发育的孔隙性储集层,用该方程解释的饱和度基本反映原状地层情况。ELAN软件ELAN软件采用双水模型,由于缝、洞的影响可能使部分层解释的含油饱和度偏低。压汞资料处理对有代表性的岩心样品,经J函数处理后转换成含油高度与含油饱和度的关系,依据油藏的平均含油高度可确定油藏的含油饱和度。核磁共振解释核磁共振测井可以求得地层可动流体和束缚流体孔隙度,由于本地区几口井均未见到明显油水界面,因此地层中的可动流体应为油,所以可用核磁测井资料解释含油饱和度胜利油区勘探开发论文集式中:So——含油饱和度,%;MBVM——可动流体孔隙度,%;MPHI——核磁共振测井总孔隙度,%。这种方法的关键是求准 T2截止值,根据岩心样品的核磁共振实验分析,古生界的T2截止值平均为1ms,太古宇的T2截止值平均为4ms。利用核磁测井资料和 T2截止值可求出每口井的含油饱和度。以上各种方法求得的含油饱和度具有不同的含义,油藏的含油饱和度选值应综合考虑。五、结论和认识埕岛油田埕北30潜山具有岩性复杂、储集类型多、非均质强的特点。利用新的测井方法——核磁共振测井,结合取心统计、试井解释、压汞处理对孔隙度、渗透率、饱和度等参数进行分析,建立了储集层参数的解释模型,并论证了参数的解释精度。利用该方法处理埕岛油田埕北30潜山5口探井,均达到较好的应用效果。该套方法也适用于类似的复杂岩性、裂缝型油藏。主要参考文献[1]周文裂缝性油气储集层评价方法成都:四川科学技术出版社,[2]柏松章等碳酸盐岩潜山油田开发北京:石油工业出版社,