1. 引言
蓬莱9-1油田位于庙西北凸起之上,毗邻渤东凹陷和庙西北洼,油源充足,油气成藏条件比较优越。新近系(明化镇组和馆陶组)直接披覆在元古宇古潜山之上,具有圈闭面积大,埋藏浅等特点。受辽河油田潜山成功勘探的影响[1] -[3] ,认为该区潜山具备油气富集的有利条件,是勘探的一个重点目标。蓬莱9-1潜山勘探钻遇大套花岗岩体,并且油气显示较好,而元古宇浅变质的围岩无油气显示或显示极差。元古宇变质石英砂岩、碳酸盐岩以及太古宇的高级变质岩和晚期侵入到太古宇中的中、酸性侵入体形成的油气藏在国内较为常见[4] -[6] ,而单纯的侵入到中新元古界地层中的花岗岩侵入体型油藏类型较为罕见。花岗岩体成藏的主控因素、岩体的空间展布以及油气富集规律等成为蓬莱9-1潜山油气藏勘探开发的重点。本文主要利用钻井取心、井壁取心、岩屑等,运用岩石铸体薄片鉴定、单颗粒锆石定年、物性分析等实验测试手段获得的基础数据资料为基础,结合测井、地震等技术,对蓬莱9-1潜山岩性特征、储层以及岩体的空间展布等进行了详实的论述,并重点对花岗岩侵入体岩石类型的划分、岩石相变及储层成因类型、物性特征等进行剖析,明确了蓬莱9-1油田潜山勘探方向,为花岗岩体储量申报奠定了扎实的基础。同时,本文研究形成的技术方法及研究结果,对类似地区的花岗岩储层以及结晶基岩储层研究,具有一定的借鉴意义。
2. 地质概况
蓬莱9-1油田位于渤海中东部的庙西北凸起上,紧邻富生烃的渤东凹陷和庙西北洼陷,具有双凹供烃、凸起聚烃的优越成藏背景[7] (图1)。
渤海湾盆地前古生界基底形成之后,先后经历了加里东、海西、印制、燕山和喜马拉雅运动等多期构造及改造作用[8] 。早古生代的加里东运动结束了稳定的克拉通发展阶段,开始了晚奥陶世–早石炭世长达140 Ma的构造抬升、风化剥蚀。三叠纪末华北地区处于南北挤压应力场之中,庙西北凸起初步形成。中生代晚三叠世至早、中侏罗世内陆盆地发展阶段,华北板块受扬子板块碰撞构造效应影响结束,继而受太平洋板块俯冲、碰撞作用影响,华北东部地壳增厚,局部地区在一定深度和温度条件下壳幔重熔形成岩浆,就近沿张性断裂呈整合型“岩盖”式侵入到元古界变质岩地层之中,即庙西北凸起鞍部的花岗岩侵入。中生代晚侏罗至白垩纪裂陷盆地发育阶段,进入受太平洋构造域控制新时期,同时诱发NNE向郯庐断裂发生大规模的左旋走滑[9] -[11] ,进入岩石圈减薄的峰期[12] ,同时出现地幔柱活动,伴随着软流圈的隆升。在上述两种机制的联合作用下,中国东部进入大规模的裂陷或断陷盆地发育阶段,沉积了

Figure 1. The explore array map for Penglai 9-1 structure
图1. 蓬莱9-1构造勘探部署图
较厚的J3-K1地层,并伴有强烈的火山活动。晚白垩世,整体处于挤压抬升构造环境,庙西北凸起区出露地表,由于花岗岩体的侵入,在侵入体的上部产生张性裂缝,加上岩性不同而差异风化剥蚀,花岗岩体出露地表。古近纪强烈的断陷作用,渤海地区形成了一系列的箕状半地堑,构成了复杂的的多重“盆–岭”构造。庙西北凸起一直处于风化剥蚀之中,其中鞍部花岗岩体风化成良好的储层。新近纪进入岩石圈热沉降期,新近系地层披覆于庙西北凸起之上。
3. 潜山地层及花岗岩体形成时代和储层成因
3.1. 潜山地层及花岗岩体形成时代
蓬莱9-1含油气构造带自上而下钻遇的地层依次为:第四系平原组,新近系明化镇组(上段、下段)、新近系馆陶组和元古宇(未穿)。上覆盖层和基底潜山岩性分界明显,潜山上部新生代盖层为砂泥岩沉积,潜山岩性为中低级变质岩(变质石英砂岩、云母片岩等)和花岗岩。在电性上,两套地层具有明显的差别,下部地层由于岩石成岩强,达到中低级变质,以及大套岩浆侵入岩,在电性上表现为具有三高一低的特点,即具有较高的岩石密度、补偿中子、自然伽玛和较低的声波时差。
对围岩和侵入体岩心分别进行了单颗粒锆石年代分析。选取了云母石英片岩围岩(PL
9-1-3
井1435.0 m)和二长花岗岩侵入体((PL9-1-4, 1405 m)分别5 Kg左右,送到了河北地质局廊坊实验室进行了锆石分离,随后将分离好的锆石在西北大学大陆动力学国家重点实验室进行锆石年代分析。首先是在双目镜下挑出具有代表性的锆石颗粒,镶嵌在环氧树脂中并抛光至锆石颗粒的一半,然后进行锆石的光学、阴极发光(CL)显微图像及LA-ICP-MS分析。其中,投射光、反射光照相在该实验室Nikon显微镜下完成;CL图像分析在该实验室的Gatan阴极发光MonoCL3+及Quanta 400 FEG热场发射环境扫描电子显微镜下完成;锆石微区U-Pb年龄测定在西北大学大陆动力学国家重点实验室最新引进的Hewlett packard公司最新一代带有Shield Torch的Agilient
7500a
ICP-MS和德国Lambda Physik公司的ComPex102 Excimer激光器(工作物质ArF,波长193 nm)以及MicroLas公司的GeoLas
200M
光学系统的联机上进行。激光束斑直径为30 μm,激光剥蚀样品的深度为20~40 μm,详细分析步骤和数据处理方法参见文献[13] 。云母石英片岩获得了三个年龄集中段(表1、图2),2685 Ma、2423 Ma和1793 Ma (由于变质程度低没有获得变质年龄),说明该套岩石的沉积年龄应晚于1793 Ma,即晚于古元古代。花岗岩体锆石年代分析(表2、图3),19个点的年龄结果集中,并都在协和线上,花岗岩的侵位年龄为164.6 Ma,形成于中侏罗世。与吴福元等[14] ,研究的中国东部花岗质岩浆作用的划分的第二个阶段(180~156 Ma)吻合。该时期东亚多板块拼贴运动学发生重大调整,构造体制发生重大转换,形成了以内陆俯冲和内陆造山为特征东亚汇聚构造体系等[12] 。
3.2. 花岗岩储层成因机理
花岗岩体在漫长的地质历史时期,经历了岩体结晶收缩、风化剥蚀、溶蚀淋滤、构造破碎等改造,储集空间较发育。尤其是潜山在新生代新近系馆陶组时期才被覆盖,长期风化淋滤,物理和化学成因的储集空间发育,岩溶作用相当强烈。李振宏[15] 等研究裸露期风化壳岩溶模式,从上到下大致可分五个带:残积带、垂直淋滤带、水平深流溶蚀带、溶蚀沉积共生带和滞留带,该种模式虽然是在研究碳酸盐岩储层时建立起来的,但同样适用于火成岩储层。火成岩岩体易形成收缩缝,在构造应力作用下,易形成构造裂缝,经历较长时期的风化淋滤,形成的风化淋滤带厚度仅次于碳酸盐岩,但远大于砂质岩。王艳忠[16] 等人对不整合面上下岩性特征进行了研究,认为不整合应包括不整合面上下岩石三部分,不整合面之上的叫残积带,主要为底砾岩和水进砂体;不整合面之下的叫风化粘土层和半风化岩石,风化粘土层位于风化壳最上部,主要是细粒残积物,是识别不整合面的重要标志,在上覆沉积物压实下岩性较致密,具有良好的封盖能力;半风化岩石,在不整合形成过程中,褶皱及断裂作用可形成大量构造裂缝[17] ,上

Table 1. LA-ICP-MS U-Pb analytical results of the mica quartz
表1. 云母石英片岩样品LA-ICP-MS U-Pb同位素分析结果

Table 2. LA-ICP-MS U-Pb analytical results of monzogranite
表2. 二长花岗岩样品同位素分析结果

Figure 2. U-Pb concordia diagram of zircons in the mica quartz
图2. 云母石英片岩锆石U-Pb谐和图

Figure 3. U-Pb concordia diagram of zircons in monzogranite
图3. 二长花岗岩锆石U-Pb谐和图
覆地层的剥蚀产生大量的卸荷裂缝,大气淡水沿早期形成的裂缝下渗,是下覆岩层发生岩溶形成大量的风化裂缝和溶蚀孔洞[18] ,因此,储层主要分布在半风化岩石中,所以这类储层从高部位向低部位半风化岩层厚度逐渐变薄。
4. 花岗岩体的岩性识别
4.1. 岩相学鉴定
蓬莱9-1构造在项目研究的2012年,共钻井13口,其中4口井(PL
9-1-3
、PL9-1-6、PL9-1-9、PL9-1-12)揭示的为元古宙中低级副变质岩,其余9口井均钻遇大套酸性侵入体,并且9口井受到强烈的风化剥蚀,其中PL9-1-4、PL9-1-5、PL9-1-7、PL9-1-8、PL9-1-13、PL9-1-14井元古宙云母石英片岩和石英岩再沉积形成了残积层,其余PL9-1-1、PL9-1-2、PL9-1-6井新近系地层直接覆盖在侵入体之上。潜山目的层取心5口井,4口井取到酸性侵入体(即PL9-1-2、PL9-1-5、PL9-1-11、PL9-1-13),同时有三口井旋转井壁取心取到酸性侵入体(PL9-1-2、PL9-1-4、PL9-1-5),为进一步认识花岗岩体奠定了基础。由于岩浆的结晶分析作用,广义的花岗岩体在不同井段和不同井在矿物组成和含量上存在差异性,反映在测井曲线上表现为自然伽玛、岩石密度、补偿中子的差异性。岩矿鉴定主要岩石类型为花岗闪长岩、二长花岗岩和石英二长闪长岩。三类岩石以灰白色、浅粉色为主,均具有半自形粒状结构,在矿物成分上表现为石英、斜长石和钾长石相对含量的差别[19] [20] 。
4.2. 岩性的测井识别
变质岩和岩浆岩岩石类型的测井识别有人做过研究[6] ,并初步建立了岩浆岩和变质岩大类的测井识别方法。测井曲线识别岩性的原理主要是依据每种岩性的矿物组合及元素组成的不同,反映到常规测井上尤其是自然伽玛、岩石密度和补偿中子等数值及形态特殊的差异性,建立起岩性与测井曲线的对应关系,进而达到测井曲线识别岩性的目的。无论是变质岩还是岩浆侵入体等结晶基岩,主要造岩矿物可分为硅铝矿物和铁镁矿物两个系列。硅铝矿物矿物中SiO2与Al2O3的含量高,不含FeO和MgO,包括石英、钾长石和斜长石等,它们基本不含色素原子,颜色较浅,又称浅色矿物;铁镁矿物矿物中FeO、MgO含量较高,包括橄榄石、辉石、角闪石、黑云母等,这些矿物的颜色一般较深,又称暗色矿物。岩浆侵入体随着岩石类型由基性向酸性过度,主要造岩矿物发生明显变化,铁镁矿物含量降低,硅铝矿物矿物含量增加,反映在测井曲线上自然伽玛质增大,原因是根据鲍文反应序列[9] ,Th、U、K等放射性元素主要进入岩浆晚期结晶的矿物,如钾长石、黑云母等,对于早期结晶的矿物斜长石、辉石、角闪石等自然加码值低。利用测井曲线识别岩性,需要搞清测井参数所反应的岩石物理特征,每一种矿物有特征测井相应值。通过矿物的测井相应特征,来研究每一种岩石的测井曲线特征,建立岩性与测井曲线的对应关系,从而实现利用测井曲线恢复巨厚的基岩岩性剖面,最终有效划分储集层和非储集层。
以往的研究成果只做到了利用测井曲线识别岩浆岩基性、中性和酸性大类,但还没有达到对具体岩石类型的细分。本次对花岗岩侵入体的研究,建立了测井曲线与具体岩石类型(如二长花岗岩、花岗闪长岩等)的对应关系,同时还建立了不同围岩的测井响应特征(表3)。通过对花岗岩岩体不同类型的岩石的测井识别,以及在单井中与围岩接触关系,以此来研究岩浆演化特点以及与地震结合确定花岗岩侵入体的侵入方向、分布面积等,为储量计算提供依据。
潜山中酸性侵入体,包括二长花岗岩、花岗闪长岩、石英二长闪长岩;元古宇围岩变质石英砂岩、石英岩和云母石英片岩等的岩石密度和补偿中子交汇特征及自然伽玛曲线特征见图4。
二长花岗岩,自然伽玛曲线平直状,值在90 API~110 API;岩石密度和补偿中子交汇曲线一般补偿中子曲线在右,岩石密度曲线在左,“小的正差异”或“绞合状”,岩石密度
2.45 g
/cm3~2.65 g/cm3,补偿中子小于6%。
花岗闪长岩,自然伽玛曲线平直状,值在75 API~90 API;岩石密度和补偿中子交汇曲线一般补偿中子曲线在右,岩石密度曲线在左,“小的正差异”或“绞合状”,岩石密度
2.45 g
/cm3~2.65 g/cm3,补偿中子小于6%。
石英二长闪长岩,自然伽玛曲线平直状,值在60 API~75 API;岩石密度和补偿中子交汇曲线呈“绞合状”,岩石密度大于2.65 g/cm3,补偿中子小于6%~12%。
变质石英砂岩,测井曲线上表现为:中等密度、较大中子、较大声波时差,较高自然伽玛等特征,自然伽玛曲线呈“锯齿状”,岩石密度和补偿中子交汇曲线一般补偿中子曲线在左,岩石密度曲线在右,“小的负差异”或“绞合状”。
石英岩,测井曲线上表现为:中等密度、较低中子、低自然伽玛等特征,自然伽玛曲线呈“低–平直”状,岩石密度与补偿中子较平直的“绞合状”。
云母石英片岩,测井曲线上表现为:较高密度、较大中子、较高自然伽玛等特征,自然伽玛曲线呈“高–平直”状,岩石密度和补偿中子呈较大的“负差异”。

Table 3. The logging identification standard of oilfield PL9-1
表3. PL9-1油田岩性测井识别标准

Figure 4. Logging response characteristics of major rock types
图4. 主要岩石类型测井响应特征
4.3. 花岗岩体的空间展布
花岗岩侵入体在侵入过程中由于岩浆的分异作用,发生分离结晶,岩浆分离结晶的总趋势是愈到晚期,岩浆愈向富硅、富碱和铁代替镁的方向演化。花岗岩体岩性分布的特点是由东南向西北由二长花岗岩(PL
9-1-1
、PL9-1-2、PL9-1-4)演变为花岗闪长岩(PL9-1-5、PL9-1-11、PL9-1-13等),二长花岗岩为岩浆分离结晶晚期的产物。根据鲍文反应原理,矿物结晶顺序:橄榄石→Ca质斜长石(Mg辉石)→Ca-碱质斜长石(Mg-Ca辉石)→碱质-Ca斜长石(角闪石)→碱质斜长石(黑云母)→钾长石→白云母→石英,同样是越晚期形成的岩石钾长石和石英越多。在流动分离结晶过程中,早期结晶的矿物容易随着流体带走,也就是岩浆在地下深处早期可形成少量斜长石等晶体,在侵入到地层中随着岩浆的流动被带走,形成的岩石斜长石等含量较高,同时早期向前侵入的岩浆本身偏中性,形成的岩体会在岩体侵入的前端方向较末端盆钾长石等。因此,根据岩浆流动分异作用的特点,推断,岩浆由东南部大断裂上升,向西北部侵入,造成东南部为二长花岗岩,西北部为花岗闪长岩(图5),同时岩浆在纵向上也具有结晶分析特点,PL9-1-13、PL9-1-5井等下部以石英二长闪长岩、石英闪长岩为主,向上过渡到花岗闪长岩和二长花岗岩。根据地震剖面反映的该岩体的形态、大小及与围岩的关系等,确定该岩体产状为岩盖。此花岗岩体岩盖整体为一旋转45˚的正方形形态,只是在西南角向外有一定的延伸,本岩盖东南–西北轴长13.6 km,东北–西南轴长12.7 km,最长处17.1 km,总面积166 km2。
5. 储集空间类型及物性特征
本区潜山岩性主要为一套中新元古代中低级副变质岩和中生代岩浆侵入体,4口井钻遇中新元古代围岩(云母石英片岩、变质石英砂岩和石英岩),无油气显示;11口井钻遇中生代花岗岩体,钻井取心和井壁取心具有较好的油气显示,因此,花岗岩体是本区的主要储集岩。
5.1. 储集空间类型
岩心观察、铸体薄片图像分析、以及毛细管压力曲线等资料均表明花岗岩体储集空间具有双重介质的特点,储集空间分为孔隙型和裂隙型两大类,主要的储集空间类型:构造裂缝、构造–溶解缝、收缩缝、破碎粒间孔及溶孔(图6)。研究发现储集空间组合类型以裂缝–孔隙型为主,次为裂缝型。

Figure 5. Cross-well seismic profile of Penglai
9-1- 13
and Penglai 9-1-2
图5. 蓬莱
9-1-13
-蓬莱9-1-2联井地震剖面

(a) (b)
(c) (d)
Figure 6. The reservoir space types in the granite. (a) broken intergranular pores (PL
9-1-2
well, 1285.0 m, Single-Mode polarizing), (b) dissolution pores (PL
9-1-13
well, 1514 m, Single-Mode polarizing), (c) tectonic fractures (PL9-1-13 well, 1503.50 m, Single-Mode polarizing), (d) shrink-dissolve fractures (PL9-1-11 well, 1605.0 m, Single-Mode polarizing)
图6. 花岗岩体储集空间特征。(a) 破碎粒间孔(PL
9-1-2
井,
1285.0 m
,单偏光),(b) 溶孔(PL
9-1-13
井,
1514 m
,单偏光),(c) 构造裂缝(PL
9-1-13
井,1503.5 m,单片光),(d) 收缩-溶解缝(PL
9-1-11
井,
1605.0 m
,单偏光)
5.2. 物性特征
对8口井,106块样品物性分析结果表明,花岗岩体储层孔隙度0.39%~25.4%,平均值为5.49%。孔隙度1%以下占12.3%,1%~5%占51.9%,5%~10%占20.8%,10%以上的占15.1%;渗透率0.07 mD~94.2 mD。其中,大于50 mD占4.9%,(10~50) mD占7.3%,(1~10) × 1 mD占4.9%,1 mD以下的占82.9%。根据孔隙度、渗透率划分标准[21] ,孔隙度以Ⅲ类为主,次为I、II类,渗透率以IV为主。孔隙度、渗透率具有良好的相关性。
6. 储层有利区带分析
根据试油资料结合测井、钻井取心、旋转井壁取心、气测录井资料分析论证,本区花岗岩侵入体潜山油藏整体含油,油层主要分布在潜山半风化岩石中,受构造、风化程度等因素的控制,在潜山上部、靠近大断裂等构造高部位,由于风化程度较强,油层较发育,如PL
1-1-1
井、PL1-1-2井、PL1-1-5井等井有效厚度都在百米以上;而对于潜山下部、远离大断裂等构造低部位来说,由于风化程度较弱,甚至没有风化,油层不发育,有效厚度较薄,如PL1-1-13井,该井处于构造较低位置,花岗岩侵入体几乎没有风化,仅在潜山中部,由于受到构造及垂直淋滤等作用的影响,发育部分裂缝型油层。油层从构造高部位向低部位半风化岩层厚度逐渐变薄。花岗岩体潜山油藏以风化壳型为主,次位构造内幕型,因此,构造高部位或靠近大断裂为有利勘探区域。
7. 结论
1) 蓬莱9-1潜山花岗岩体单颗粒锆石U-Pb同位素年龄为164.6 Ma,为中生代侏罗纪中侏罗世时期侵入到元古宇地层中的侵入体。
2) 花岗岩侵入体主要岩石类型为二长花岗岩、花岗闪长岩和石英二长闪长岩,岩体具有横向结晶分析和纵向结晶分析等特点,在自然伽玛测井曲线上有很好的反映。
3) 花岗岩体在本区形成了很好的储层,与岩体自身矿物组成密切相关,花岗岩体主要由石英、斜长石、碱性长石等浅色矿物组成,而且长石含量达到55%~70%,岩石脆性大,抗风化能力弱易形成储层;而围岩主要矿物成分为石英、云母类矿物,石英抗风化能力强,而黑云母则蚀变形成绿泥石等柔性矿物,因此,不易形成有效储集空间。
4) 储集空间类型主要为风化淋滤、构造破碎等形成的破碎颗粒粒间孔、溶孔、收缩缝、溶解缝和构造裂缝等,储集空间组合类型以裂缝+孔隙型为主,次为裂缝型。
5) 花岗岩体油藏以风化壳型为主,次为构造型,有利勘探区域为构造高部位和大断裂附近。
基金项目
国家重大科技专项(2011ZX05023-002,2011ZX05023-006)资助。