1. 引言
随着油气勘探程度的深入，基岩油气藏也越来受到广泛的重视。目前世界范围内的基岩潜山油气田中，40%产于花岗岩潜山储层中，储量约占基岩油气田储量的75% [1] [2] ，以越南的白虎油田和锦州25-1S油田较为著名。国内外勘探开发实践表明，花岗岩潜山储层非均质性强，且内幕变化复杂，储层的纵横向分布对油气的富集具有明显的控制作用。花岗岩原岩致密，只有经历后期风化淋滤以及构造裂缝改造才能形成有效储层，因此从成因上来说该类油藏多为风化壳型油藏。
国内外对碳酸盐岩风化壳型储层的研究起步较早且较为成熟 [3] [4] ，而针对花岗岩类风化壳型储层分带性尚处于起步阶段。钻井取心和野外露头观察表明花岗岩类风化壳型储层同样具有垂直分带特征，其优势储层集中在风化壳附近100 m以内，特别是在50 m左右最好。国内外不同学者根据风化壳的结构、岩矿特征、化学指标等对这些带进行了定性或半定量的划分 [5] ，主要是依据风化程度的强弱划分为极强风化带、强风化带、次风化带、弱风化带，本文中沿用以上分类标准。但在实际油田操作中，不同风化带界面识别存在困难，潜山内幕地震成像质量差，且测井响应不明显，而又无系统取芯资料，由此制约了储层地质模式的建立和优质储层的预测。
2. 研究区概况
本次研究区(渤海PL油田)位于渤海庙西北凸起之上，整体构造为依附于凸起东界大断层发育的大型半背斜，潜山南北高点为元古界变质岩潜山，中间宽缓鞍部为中生界花岗岩侵入体。PL油田为国内海上第一个亿吨级花岗岩潜山油田，发育裂缝孔隙型储层，潜山岩性较为单一、缺乏常规沉积岩旋回特征，纵向分带难度大。本文综合利用PL油田钻录井、测井、分析化验等方面的资料，以及XRF元素录井和INPEFA频谱属性分析等技术手段，形成了一套新的潜山分带界面识别技术方法，在实际应用中取得了良好的效果，有效支撑了潜山地质模式的建立。
3. 花岗岩潜山地层纵向分带识别技术
3.1. 钻录井分带特征综合识别技术
3.1.1. 岩心、薄片识别
通过大量的岩心和薄片观察，认为潜山纵向上主要分为极强风化带、强风化带、次风化带和弱风化带，其中极强风化带因遭受风化作用最为严重，花岗岩中长石多风化为高岭土等粘土矿物，孔隙多为泥质充填，无储集性能(图1(A))；强风化带岩心多呈散砂状和碎块状，铸体薄片见大量溶蚀孔隙，微裂缝也极其发育(图1(B))；次风化带取心整体致密，仅在裂缝发育段储层得到改善，裂缝面含油，见少量溶蚀孔隙(图1(C))。
(A) (B) (C)(A) 1284 m~1288 m花岗岩遭强烈风化，松散呈土壤状。(B) (a) 风化较重，疏松，裂缝发育，裂缝面含油；(b) 铸体薄片(单偏光25×)孔缝发育，见角闪石溶孔和溶解构造缝；(c) 1605.79 m，荧光薄片，未充填裂缝发中亮-亮绿色荧光。(C) 1495 m~1495.85 m，基质致密，发育高角度裂缝。
Figure 1. Identification characteristics of cores and thin sections in each weathering zone
图1. 各风化带岩心、薄片识别特征
3.1.2. 钻石曲线识别
钻时反映的是岩石可钻性，可钻性与岩石本身的致密程度密切相关，而岩石的致密程度主要是由岩性和物性决定的，因此，钻时能够反映岩石的岩性和物性，而且不受岩石本身的电性、磁性或放射性等因素的影响 [6] 。而本区储层纵向上岩性变化不大，在钻头、钻压、转盘转数等参数稳定时钻时基本上反应了岩石物性纵向上的差异。
图2中PL11井钻时曲线在强风化带和次风化的分界处可见明显的钻时异常。
Figure 2. Integrated logging map of well PL11
图2. PL11井综合录井图in well PL5
3.2. X射线元素录井(XRF)界面识别技术
X射线荧光元素录井(通常用英文缩写XRF表示)技术是近年来兴起的一项新技术，通过X射线脉冲数来记录元素产生的荧光能量和强度，定性地指示元素含量的高低。XRF技术目前尚处于应用研究阶段，目前国内很少开展X射线元素录井地层划分的工作，仅有少数油田进行了探索性试验与研究 [7] ，PL潜山多数评价井均采用XRF技术开展录井。
XRF技术的原理在于各种矿物的抗风化能力不同，各种元素的迁移能力不同，在不同的风化作用阶段，元素的富集情况和形成的矿物也不一样。但风化作用最终会造成不稳定矿物(铁镁矿物、长石类矿物等)向稳定矿物(石英和黏土矿物)的转变；易迁移元素(Cl、S、Ca、Mg、Na、Mn、硅酸盐中的Si、P等元素)的流失和弱迁移元素(Fe、Al、Ti等)、不迁移元素(Si)的富集。因此，从理论上能够通过元素的不同及含量的变化来进行地层的划分与对比。
以PL5井为例，不迁移元素SI的富集以及易迁移元素Ca，P，Na，Al的流失，很好的说明了元素在垂向上规律性的变化(图3)。
Figure 3. Vertical distribution of rock components
图3. PL5井岩石组分纵向分布图
3.3. 基于“测井频谱分析”的测井识别技术
从PL潜山15口评价井测井曲线特征来看，个别井纵向上可以通过曲线基值和形态的变化识别出各个风化带，尤其是对于极强风化带顶底界的识别。但就强风化带、次风化带底界却很难形成一套统一的标准。
为此引入了INPEFA测井频谱分析技术，INPEFA测井频谱分析技术是荷兰沉积学家S.D.Nio研发的一种以最大熵分析方法为基础的测井曲线频谱分析方法，该方法通过一系列的数学运算把使曲线由深度域转化为频率域，使隐藏在测井曲线中的地层特征更加显现化。其原理主要是已经依据“求异对比”的理念—原始的花岗岩在未经风化构造改造这一因素影响时是不会有频谱曲线纵向上的差异的，相似的风化构造改造程度有相似的频谱曲线响应特征，因此只需要把这一因素导致的测井响应差异放大，形态相似性好的曲线其地层发育也相似，就可以实现地层多井对比 [8] 。
以PL2和PL4井为例，分别对GR，DT，DEN三条测井曲线做了整体INPEFA频谱分析，利用该技术可以将原始测井曲线中的趋势或是差异放大，使地层界面易于识别(图4)，风化带底界DT呈明显的负转折点，DEN则呈正转折点。
Figure 4. Stratigraphic division correlation map by log spectrum analysis
图4. 测井频谱分析划分对比图
4. 潜山地层纵向分带综合识别图版
在多种潜山纵向分带技术应用的基础上，最终确定了PL油潜山地层划分图版(图5)，各风化带综合识别特征如下
Figure 5. Vertical zoning identification chart of granite buried hill in PL oilfield
图5. PL油田花岗岩潜山纵向分带识别图版
极强风化带位于花岗岩侵入体顶部，录井表现为泥质含量高，低钻时特征，测井响应则为高GR、低电阻、低速的特征。PL2井在极强风化带内的取心证实，本带风化非常强烈，部分呈粉砂质泥状结构，难以形成有效储层。地球化学元素特征上以富Si、Fe等难迁移元素，贫Ca、Na等易迁移元素为特征，岩矿特征上则以高石英含量，低长石含量为特征，该带风化成熟度最高，化学蚀变程度也最强。
强风化带位于花岗岩侵入体的上部，经过强烈风化、淋滤作用，溶蚀孔、洞，裂缝均较发育，钻时较下部明显偏低，测井响应表现为低GR，较高电阻，较高波速的特征，INPEFA频谱曲线上DEN为负转折点，DT为正转折点。综合解释证实，该带储层发育，横向分布稳定，连通性较好。本带岩心，多见散砂状和碎块状储层，铸体薄片见大量的颗粒溶蚀孔和部分粒间孔，构造裂缝发育，形成网状缝，沟通溶孔，形成很好的储集空间，以孔隙型储层和裂缝-孔隙型储层为主。岩石地化以及岩矿特征上较极强风化带SI含量降低，Ca、Na等含量升高，石英含量降低，长石含量升高，该带化学风化成熟度降低，但较下部次风化带和弱风化带仍较高。
次风化带的风化、淋滤作用变弱，溶蚀孔、洞较少，以裂缝为主，钻时曲线基值较上部明显增加，表现为低GR、高电阻、高速的测井响应特征，INPEFA频谱曲线上DEN为负转折点，DT为正转折点。致密段较多，储层发育变差，横向分布不稳定，连通性较差。取心证实，本带存在有效储层，但整体比较致密。
弱风化带钻遇井较少，该带裂缝较发育，致密带厚度增加，整体物性变差，钻时基值较上部增加，但在裂缝发育段钻时降低明显，测井曲线上表现为低GR、高电阻、较高速的特征。铸体薄片观察仅在构造缝周围发育有少量的溶蚀孔隙。
5. 应用效果分析
该潜山纵向分带综合识别图版在实际应用显示出较高的技术优势，一方面在PL18和PL19评价井钻进中通过钻时曲线和XRF技术准确、快速的确定潜山界面；另一方面各井分带的技术成果结合潜山古地貌分析，形成了由高部位向鞍部中心及北西方向环带状依次减薄的“簸箕圈”模式，模式对于后期有利储层的预测，开发井位的优选起到了很好的指示作用(图6)。
Figure 6. Reservoir distribution model in PL buried hill
图6. PL潜山储层发育模式
6. 结论
1) 渤海PL花岗岩潜山具纵向分带性，依据风化程度的不同将其分为极强风化带、强风化带、次风化带和弱风化带；
2) 综合利用PL油田钻录井、分析化验、测井等方面的资料，引入XRF元素录井和INPEFA频谱属性分析等技术手段，形成潜山纵向分带综合识别图版，能弥补不同资料间地层划分对比当中存在的不足，在实际应用中具有较高的技术优势。