1. 引言
随着石油钻探认识的不断增加,具有工业性开采价值的潜山油气藏越来越引起人们的重视。渤海JZ区块发现的大型变质岩潜山油气藏与典型的碳酸盐潜山油气藏不同,其主要靠风化壳和裂缝储集油气。在钻进变质岩潜山的作业过程中,常常会在潜山界面之上钻遇一层较薄的砂砾岩和砂岩。该岩层特征与潜山及其上覆岩性具有差异性,由于其主要分布于潜山顶面缓坡区,故称之为“坡积岩”。坡积岩层横向发育无规律,预测难度大;纵向上较薄,由于随钻测井测量工具和钻头的距离差,无法及时判断是否钻进坡积岩,在潜山界面、取心位置和下套管位置的判断上等存在较大困难,造成钻井风险增大。因此,为了降低钻井风险,获取更为精确的地质资料,笔者在分析坡积岩层的岩石特征和录井工程参数响应特征的基础上,利用钻井工程参数建立功指数和机械比能模型,实现了对潜山坡积岩层的随钻判别,达到了降低钻井地质作业风险的目的。
2. 工区地质背景
渤海JZ区块变质岩潜山位于辽东湾辽西低凸起中北段,西侧以辽西大断层与辽西凹陷相邻,东侧向辽中凹陷下倾,毗邻辽中凹陷中、北洼,处于油气富集的有利位置 [1] 。自上而下主要发育新近系明化镇组(N1m)和馆陶组(N1g)、古近系东营组(E3d)、沙河街组一~三段(E3s1、E3s2、E3s3),基底为上太古界,缺失古近系沙河街组四段、孔店组,中生界,古生界和元古界 [2] 。变质岩潜山油藏主力产层为太古宇变质岩,盖层为上覆E3s3泥岩,属风化体块状储集层地貌潜山 [3] 。
根据岩性特征和裂缝发育特征,变质岩潜山可分为半风化壳上段、半风化壳下段、基岩等。在研究区内,潜山上部的风化黏土层通常被剥蚀掉,在半风化壳之上直接沉积E3s3 [4] 。E3s3沉积特征复杂、相变快、岩性类型多样,尤其是E3s3底部局部发育有坡积岩层。坡积岩层以下为花岗质变质岩,受长期风化影响,变质岩上部裂缝及溶洞十分发育,部分裂缝被充填。
3. 坡积岩层的岩石特征
尽管在油田实施开发井钻井作业中,受作业手段限制,缺少坡积岩层的测井等资料,但通过现场岩屑和岩心资料及探井试验资料,可对坡积岩层进行较为系统的分析。
分别对研究区的坡积段砂岩、E3s2砂岩和E3s3砂岩开展了岩石薄片分析鉴定,通过三角图(图1)投点发现,E3s2砂岩主要落在V和VI区域,属于岩屑长石砂岩和长石岩屑砂岩;E3s3砂岩主要落在VI区域,属于长石岩屑砂岩;坡积岩层砂岩分为2种,一种是长石含量较多的长石砂岩,另一种是岩屑含量较多的岩屑砂岩。

Figure 1. The delta pitch point diagram of sandstones in the studied area
图1. 研究区砂岩三角投点图
研究区的E3s2砂岩、E3s3砂岩以及坡积岩层砂岩的岩石薄片体积分数具有较大差异(图2)。E3s2砂岩和E3s3砂岩的成分及其含量较为一致,以石英、长石为主;相比之下,坡积岩层的长石砂岩和岩屑砂岩的含量差异较大;坡积岩层岩屑砂岩的火成岩岩块体积分数较高,达到了66%,使得石英、长石的体积分数较低;坡积岩层长石砂岩不含有火成岩岩块,其钾长石比其他3种岩石的体积分数平均值都要高。

Figure 2. The histogram of sandstone composition and its average volumetric fraction (rock slice identification) in the studied area
图2. 研究区砂岩成分及其平均体积分数柱状图(岩石薄片鉴定)
图3为E3s2砂岩、E3s3砂岩以及坡积岩层砂岩的粒度分布图,分析认为坡积岩层岩屑砂岩和长石砂岩的粒度中值要大于E3s2砂岩、E3s3砂岩的粒度中值。坡积岩层岩屑砂岩的磨圆度差,有粗粒也有细粒,且分布范围较广;坡积岩层长石砂岩的粒径主要分布在1 mm左右,最大2.6 mm,比岩屑砂岩的分布范围小;E3s2砂岩、E3s3砂岩粒度较小,在1 mm以下。

Figure 3. The distribution of sandstone grain size in the studied area
图3. 研究区砂岩粒度的分布图
通过上述资料可以看出,坡积岩层分为岩屑砂岩和长石砂岩2种类型的砂岩;在成分和粒度对比上,岩屑砂岩更趋近于花岗岩,长石砂岩更趋近于正常砂岩,表明该坡积岩层应为2种不同成因形成的砂岩,其中岩屑砂岩的形成与花岗岩存在一定关系。
4. 钻井工程参数识别模型
钻井工程参数具有实时性、连续性的优点,根据钻井工程参数建立的功指数模型和机械比能模型,两者能从不同角度反映地层的致密程度,从侧面反映岩性的变化,结合岩屑分析结果,可实现对坡积岩层的随钻识别。
4.1. 功指数模型
钻头钻进过程实际是其做功的过程,不同强度的地层其做功也有不同。要计算钻头每破碎单位体积岩石所做的功是十分困难的,目前国内外尚无较好的计算方法。因此需要引入经过变形简化后的计算模型,该模型的计算结果即为钻头的做功指数,即“功指数”。功指数与钻头做功大小呈正相关性,从而可以从侧面反映钻头的做功情况。
钻头破岩做功可简要表述为垂直纵向力做功与水平剪切力做功之和,纵向力分为静钻压和冲击钻压,水平剪切力为扭矩,则钻头每进尺一个单位深度(1 m)所做功为:
(1)
式中:Wz为钻头每进尺一个深度单位所做功大小,kJ;YJ为静钻压,kN;YC为冲击钻压,kN;R为转速,r/min;Z为钻时,min;N为扭矩,kN∙m。
YC目前尚无法获取,其余参数均可通过综合录井系统监测得到,因此需对式(1)进行变形,使之能够
直接计算并具应用价值。变形后的功指数表达式为 [5] :
(2)
式中:Wm为变形后的功指数,kJ;a为区域经验静钻压值,kN;b为区域经验转速值,r/min;c为区域经验系数,1。
4.2. 机械比能模型
机械比能理论是作为一种用来描述钻头性能的概念而被提出,它提供了一种实时评价钻井性能的工具。开展机械比能的研究对随钻监测井下工作状态,随钻岩性分析,提前预测而避免发生钻井事故,及时做出合理的钻进参数调整,提高钻井效率,降低钻井成本等方面都具有重要意义。在国外,机械比能理论已广泛应用于钻进过程的监测与预测、钻井工程设计优化、随钻岩性评价等方面,并取得了很好的应用效果。国内对机械比能的研究相对较少,主要应用于钻头选型、监测钻头磨损等。
一种完善的机械比能模型应具备的条件是:① 在岩性不同的地层内钻进,最小机械比能值应该约等于岩石的抗压强度;② 模型参数易于测量和计算;③ 适用范围广,能很好地适用于各类钻头与各种井型。
通过分析不同模型的优劣,并结合渤海JZ地区的特点,研究区使用樊洪海模型 [6] :
(3)
式中:EMS为机械比能,MPa;Ab为钻头面积,mm2;vROP为机械钻速,m/h;Db为钻头直径,mm。
4.3. 模型识别规律
通过对渤海JZ区块功指数和机械比能曲线响应特征的分析及研究,确定模型系数,计算识别特征值范围,从而达到随钻识别的效果。渤海JZ地区的功指数模型系数,E3d为:a = 145 kN,b = 125 r/min,c = 0.8;E3s (沙河街组)为:a = 180 kN,b = 127 r/min,c = 0.8;太古界为:a = 212 kN,b = 103 r/min,c = 0.8。
总结了研究区E3s3泥岩、坡积岩层、太古界花岗岩的功指数、机械比能曲线变化规律(见表1)。该区E3s3泥岩功指数较低,一般小于2 kJ;太古界花岗岩较高,一般大于8 kJ;坡积岩层功指数则呈现出高低值跳跃变化,一般分布于1~8 kJ之间。机械比能与功指数曲线类似,将其放置同一图道中,赋予相同刻度,则E3s3泥岩段呈填充特征,太古界花岗岩段呈重叠特征,坡积岩层段则呈现出曲线填充、重叠交替变化特征。

Table 1. The characteristics of work index and mechanical specific energy identification
表1. 研究区功指数、机械比能识别特征表
5. 实例应用分析
以JZ-X1井为例,从图4中可以看出,E3s3泥岩与太古界花岗岩功指数特征差异明显,E3s3泥岩功指

Figure 4. The method of model identification in Well JZ-X1
图4. JZ-X1井模型识别方法图
数较低,平均在0.6左右;太古界花岗岩潜山功指数则较高,平均在100左右;而坡积岩功指数则出现高、低值跳动的现象,最高达8,最低仅为1。分析认为,功指数高、低变化值主要与坡积岩自身的致密程度相关,坡积岩越致密,其功指数越高,坡积岩越疏松,则功指数越低。将功指数与机械比能曲线放置在同一图道,赋予相同刻度,机械比能与功指数曲线形态呈类似特征。机械比能与功指数曲线在E3s3泥岩段呈填充特征,在太古界花岗岩段呈重叠特征,在坡积岩层段却呈现出曲线填充、重叠交替变化特征。
应用该方法对研究区21口实际井资料的处理与验证,随钻识别率高达86.7%,表明该方法效果较好,可用于研究区的坡积岩层随钻识别应用。
6. 结论
1) 渤海JZ区块变质花岗岩潜山中的坡积岩层位于E3s3泥岩与太古界变质花岗岩之间,分为2种类型砂岩,一种是长石含量较高的长石砂岩,其特征与正常沉积的砂岩相似,另一种是岩屑含量较高的岩屑砂岩,其特征与下覆花岗岩特征相似。表明该坡积岩层应为2种不同成因形成的砂岩。
2) 功指数、机械比能曲线组合能够较好地对坡积岩层进行随钻识别。在E3s3泥岩段,功指数呈低值,功指数与机械比能曲线呈填充特征;在太古界花岗岩段,功指数呈高值,功指数与机械比能曲线呈重叠特征;在坡积岩层段,功指数呈高、低值跳跃变化,且功指数与机械比能曲线亦呈填充、重叠交替变化特征。
3) 功指数模型、机械比能模型构建参数简单、易得,有利于后期的推广及应用。在潜山界面识别、潜山内部储层储集性能快速识别等方面具有广泛的应用价值。