史帅雨,余一欣,殷进垠,邬长武,刘静静,刘延莉,王 波
[1.中国石油大学(北京) 油气资源与探测国家重点实验室,北京 102249; 2.中国石油大学(北京) 地球科学学院,北京102249;3.中国石化 石油勘探开发研究院,北京100083; 4.中海油能源发展股份有限公司 工程技术分公司 中海油实验中心,天津300452]
含盐盆地是全球油气资源最丰富的一类盆地[1],含盐盆地中已探明的油、气储量分别占全球已探明油、气储量的89%和80%[2-3]。近些年来,位于西非被动大陆边缘的下刚果盆地发现一系列与盐构造相关的大型油气田,使其成为了全球油气勘探的热点地区[3]。前人已针对下刚果盆地的构造演化和油气聚集成藏特征等进行了研究[4-14],对盐构造的变形和演化特征也有所讨论[15-28],同时还完成了一些物理模拟实验研究[16,30-33]。但目前还缺少对盐构造变形期次及主控因素的详细分析,而且前人的物理模拟实验较少涉及对比模型,从而限制了对盐构造变形影响因素的深入分析。本文主要在前人工作基础上,利用地震剖面解释成果分析下刚果盆地不同区带盐构造变形特征,并通过平衡剖面复原和构造物理模拟实验厘定盐构造变形期次,对比分析重力滑动和沉积差异负载作用对盐构造变形与演化的影响,明确不同基底倾斜程度对盐构造发育的影响,为下刚果盆地盐构造分析以及有利油气勘探区带评价提供借鉴和参考。
下刚果盆地位于西非被动大陆边缘,北部以Mayumba隆起与加蓬盆地相接,南部以Ambriz隆起与宽扎盆地相邻,东至出露地表的前寒武系基底,西至非洲大陆坡,面积约为68 000 km2,其中海上面积约占90%(图1)[4-6]。下刚果盆地的形成演化与西非海岸其他盆地相似,经历了冈瓦纳大陆分离、裂谷发育和板块漂移等过程,属于典型的被动大陆边缘叠合盆地,可划分为前裂谷期、同裂谷期、过渡期和被动大陆边缘期等4个阶段[7-11]。
下刚果盆地沉积地层主要包括侏罗系—新生界(图1),纵向上可划分为盐下层序,盐岩层序和盐上层序[6]。盐下层序不整合覆盖在前寒武系基底之上,以晚侏罗世、早白垩世尼欧克姆和巴雷姆期的河流-湖泊-三角洲沉积体系为主。盐岩层序在早白垩世晚阿普特期沉积,厚几十米到上千米,整体上厚度由南向北变薄。盐上层序是一套裂谷后的被动边缘沉积,主要发育三角洲与深海沉积体系[12-14]。
根据下刚果盆地盐层及其上覆地层构造变形特征的差异性,从大陆架到深海平原方向可以整体划分为后缘伸展带、中间过渡带和前缘挤压带(图2)。各构造带盐相关构造样式差异明显,从陆地向海洋方向盐构造发育成熟度依次增大。
图1 下刚果盆地区域地质背景(a)及地层综合柱状图(b)Fig.1 Regional geological map and stratigraphic column of the Lower Congo Basin
后缘伸展带位于大陆架及部分上斜坡,主要发育盐滚、筏状构造和铲式正断层等。盐滚和盐筏构造是下刚果盆地典型的构造样式(图3a)。盐滚构造幅度较低,两翼不对称,缓倾一翼与上覆地层整合接触,陡倾一翼与上覆地层以正断层接触,是伸展应力作用下盐岩向下倾方向滑动的结果[20-21]。盐岩滑脱早期,盐上地层发育一系列向盆地内倾斜的铲式断层;晚期上盘滑脱速率大于下盘,上盘与下盘完全分离,不再通过断面接触,形成了盐筏构造。根据盐筏构造发育的年代,下刚果盆地的盐筏构造又可细分为白垩系盐筏构造和新近系盐筏构造。
中间过渡带位于斜坡区之上,主要发育各类不同形态的盐底辟以及龟背构造和微盆等(图3a)。该区带受拉张应力和挤压应力的共同作用,盐岩主要以垂向运动为主,呈直立状,断裂相对不发育。向陆一侧盐底辟构造幅度比盐滚构造幅度大,两翼较为对称,与上覆地层呈整合接触;向海一侧发育的盐底辟构造幅度更大,呈直立状刺穿上覆地层。盐底辟的活动对上覆地层有很大的影响,导致上覆地层通常发育龟背构造、披覆背斜和微盆地等构造样式[22-23]。
图2 下刚果盆地区域地质剖面(剖面位置见图1)[25]Fig.2 Regional geological section of the Lower Congo Basin (see Fig.1 for the location)[25]
图3 下刚果盆地盐构造地质剖面(剖面位置见图1)Fig.3 Geological section of salt structures in the Lower Congo Basin (see Fig.1 for the location)
前缘挤压带位于坡脚及局部深海平原区,以盐蓬、逆冲断层和厚层盐岩等挤压构造为主要特征(图3b)。盐蓬和盐盖构造位于挤压区向陆一侧,在挤压应力的作用下下层盐岩侵入上覆地层,并向下倾方向挤入、展布,通过补给通道与盐岩层相连,在强烈挤压作用下会伴生一些逆冲推覆构造[27-28]。
厚层盐岩(图3b)位于挤压区向海一侧,是上坡原始沉积的盐岩在重力和沉积载荷作用下向下坡滑动并充注、聚集的结果,厚度可能达到2~3 km。盆地向海前缘的区域盐岩厚度都较大,且分布均匀,形成现今仍在活动的陡坡构造。上覆地层沉积厚度较薄,以逆冲断层和挤压褶皱等构造为主。
下刚果盆地盐岩的塑性流动和变形影响了上覆地层的沉积和变形,反过来,上覆地层的沉积厚度变化和变形特征也可以反映盐构造的变形过程。本文主要利用分层复原法[29,34]编制下刚果盆地盐构造平衡演化剖面,进而分析盐构造形成演化期次。结果表明,下刚果盆地盐构造形成演化主要经历了阿尔比期—晚白垩纪初始活动、晚渐新世—上新世末强烈活动和上新世末—现今微弱活动3个阶段(图4)。
阿尔布早期地层沉积之后,在重力和上覆层差异载荷的作用下,盐岩开始发生早期的滑脱变形。阿尔布期末碳酸盐岩沉积后,热沉降作用导致陆缘向海掀斜,后缘伸展区在拉张作用下形成了一系列规模较小的断层,伴生一些沿滑脱面分布的孤立断块及同沉积地层;中部过渡区和前缘挤压区发生局部增厚和减薄,形成一系列小规模的盐枕,导致上覆地层变形形成小型凹陷和突起。在晚渐新世—上新世末,由于非洲板块构造掀斜导致大量陆缘沉积物供给,盐构造变形剧烈,伸展区盐岩向海方向发生大规模滑脱,形成盐滚和盐焊接;过渡区早期形成的盐枕构造进一步演化成为小型底辟构造,盐底辟顶部地堑开始发育;挤压区盐岩发生强烈刺穿,发育刺穿型底辟。由于伸展区盐岩大规模滑脱,导致大陆斜坡位置形成滑脱较为严重的铲式断层及反向调节断层。上新世末—现今,盐岩变形相对较弱,没有新的盐构造发育,早期形成的盐底辟缓慢、持续发育,挤压区的盐底辟在重力扩展作用下有向海方向挤入的趋势。
从下刚果盆地的构造演化特征及其平衡复原剖面来看,控制盐构造发育的动力学机制主要是重力滑动作用和上覆地层的差异负载作用。而在不同构造带内,重力滑动作用和重力扩展作用所起的程度是不一样的。在下刚果盆地伸展区和过渡区,由于大陆斜坡坡度较大,重力滑动作用占主导地位,而在挤压区重力扩展作用完全被挤压作用吸收,所以只有重力滑动作用起作用。同时由于上覆沉积物厚度向海方向减薄,在差异负载作用下,盐岩在重力滑脱变形过程中一方面向海流动,聚集在挤压变形区,同时在不同构造区内形成多种类型的盐构造。
本文在分析下刚果盆地主要盐构造变形特征和形成演化的基础上,设计了比例化物理模拟实验,通过再现下刚果盆地盐岩塑性流动及其变形过程,探讨重力滑动作用、上覆层差异负载作用和不同基底倾斜角度对盐构造发育的影响。
同前人对盐构造开展的盐构造物理模拟实验类似,本次模拟实验用硅胶来模拟具有牛顿流体特征的塑性盐层,硅胶的密度约为987 kg/m3,粘度系度约为5×104Pa·s[30,31]。用干燥松散的石英砂来模拟沉积岩,石英砂的密度约为1 500 kg/m3,粘结力约200 Pa。在自然重力场条件下,石英砂的形变遵循莫尔—库仑破坏准则,破裂内摩擦角为25°~30°,非常接近地壳浅部(< 10~15 km)沉积岩层的脆性形变行为[32,33]。
本次物理模拟实验的模型尺寸和实验材料都满足相似性原则,表1列出了下刚果盆地模拟实验模型的参数。实验中的长度比例为1×10-5,即模型中的1 cm代表自然界中的1 km,实验材料的密度约为地质原型岩石密度的一半,硅胶与盐岩的粘度比约为5×10-13,由此可以得出时间比为1×10-8,即实验中模拟1 h相当于自然界中11.4×103a。
以下刚果盆地实际地震剖面(图3a)为模拟对象,同时考虑到下刚果盆地所处大陆边缘的初始倾斜程度[17],设计了基底斜坡倾角为3°的模型1(图5a)。模型尺寸为80 cm×20 cm,模型两侧为固定端,在靠近模型右侧固定端处使用光滑的PVC板设置一个长度为80 cm的基底斜坡,其坡度为3°。硅胶的流动变形完全依靠上覆沉积物差异负载作用以及重力滑动作用来实现。在斜坡上方铺厚约2 cm的硅胶,静置约10 h,以等待硅胶自然流动使其表面水平,并排出气泡。然后开始实验,以类似大陆边缘砂体进积的形式向模型左侧铺砂。实验结束后加砂、浇水,将模型固定,最后沿长轴方向切开模型,以观察模型内部的硅胶和石英砂层的变形特征。
图4 下刚果盆地平衡复原剖面(剖面位置见图1)Fig.4 Cross-section restoration of the Lower Congo Basin (see Fig.1 for the location)
表1 下刚果盆地物理模拟实验参数Table 1 Scaling parameters of physical modeling of the Lower Congo Basin
图5 实验模型1(a)和实验模型2(b)Fig.5 Physical Model 1 (a) and Model 2 (b)
为了对比分析基底倾斜程度对盐构造发育的影响,设计了对比模型2(图5b)。与模型1相比,仅是将基底斜坡坡角改为5°,其他实验条件均不改变。
在模型1中,当实验进行约3 h后(图6b),在上覆砂层差异负载作用和重力滑动作用下,硅胶拖拽上覆层沿斜坡向下流动,砂层开始下陷和伸展,上覆砂层中也仅仅发育了两条向下倾方向伸展的断层。随着上覆砂层向下坡方向逐渐进积,硅胶的流动更为明显,同时断层持续发育,规模变大,数量也开始增多。7 h后(图6c),盐岩滑脱较为严重,斜坡上倾方向上的硅胶发生了明显的局部增厚和减薄现象,早期发育的断层被砂岩覆盖后,逐渐停止活动。至实验结束时,已发育较明显的盐滚和盐筏构造(图6d)。
在模型2中,随着上覆地层的沉积,当实验进行3 h后(图7b),硅胶整体流动不明显,出现局部增厚和减薄现象,上覆砂层中发育了两条向下倾方向伸展的正断层。7 h后(图7c),硅胶的流动更为明显,同时上覆砂层发育了几条规模较大的断层。至实验结束时(图7d),已发育较为明显的盐滚和盐底辟构造。与模型1相比,硅胶整体流动特征比较相似,但上覆砂层中的断层数量明显减少。
对比模型1和模型2的内部变形特征,可以发现在模型1中,盐岩在重力滑动作用和上覆地层差异负载作用下,形成了一系列盐滚和盐筏构造,向斜坡下倾方向,筏状构造的发育成为逐渐变新。上覆地层形成一系列向下倾方向伸展的断层,后期伴生一些反向调节断层(图8)。而在模型2中,整体变形特征与模型1比较相似,上坡方向在正断层下盘发育的盐滚数量较少,上覆砂层中形成了比较明显的筏状构造,而且在斜坡下倾方向发育规模较大、隆起幅度较高的盐底辟构造。两个底辟中间形成明显的龟背构造(图9)。两个模型有比较明显差异的是,模型1形成的盐滚和盐底辟的数量要比模型2多,但其规模明显要小一些。
图6 模型1模拟过程及结果Fig.6 Simulation process and results of Model 1a.初始沉积;b.硅胶开始流动,上覆地层形成两条断层;c.断层持续发育,规模变大;d.硅胶发生明显增厚和减薄,形成盐滚、盐筏构造
图7 模型2模拟过程及结果Fig.7 Simulation process and results of Model 2a.初始沉积;b.硅胶开始流动,上覆砂层形成两条断层;c.硅胶在斜坡中下部发生明显增厚;d.形成盐滚、盐筏和盐底辟构造
图8 模型1内部变形特征Fig.8 Internal deformation characteristics of Model 1
图9 模型2内部变形特征Fig.9 Internal deformation characteristics of Model 2
该组模拟实验结果表明,较大的基底坡度有利于重力滑动作用使盐体向海洋方向发生快速流动,从而导致盐底辟快速发育,其规模较大,但数量较少,盐构造之间的单个盐焊接面积也更大,构造分带特征更为明显。而较小的基底坡度有利于盐体缓慢流动,盐滚构造充分发育,其规模较小,但数量较多,上覆层中同时发育大量断层。
与下刚果盆地实际地质剖面(图3a)对比,发现基底斜坡倾角更大的模型2与其有更好的相似性。这可能也表明西非被动大陆边缘盐岩初始沉积在较缓的基底斜坡背景下,而后受渐新世非洲板块隆升作用影响,下刚果盆地盐下基底斜坡逐渐变陡,盐岩则发生强烈的塑性流动变形,从而形成了一系列不同类型的盐构造样式,如图9反映的在较陡基底斜坡条件下的盐构造变形特征。后缘伸展区盐岩较薄,形成小规模的盐滚和盐筏构造,上覆层中铲式正断层较为发育。中部过渡区发育规模较大的三角形盐底辟,盐底辟之间盐岩几乎流失殆尽,形成盐焊接,上覆地层形成龟背构造。前缘挤压带的盐岩发生了大规模聚集,形成厚度较大的盐岩聚集区。
1) 下刚果盆地从大陆架到深海平原下刚果盆地可划分为后缘伸展带、中间过渡带和前缘挤压带。伸展带以盐滚、盐筏、铲式断层为主要特征,过渡带以不同形态的盐底辟、龟背构造为主要特征,挤压带以厚层盐岩、盐蓬和逆冲推覆构造为主要特征。
2) 下刚果盆地盐构造的形成演化主要经历了阿尔比期—晚白垩世初始活动、晚渐新世—上新世末强烈活动和上新世—现今微弱活动3个阶段。
3) 比例化物理模拟实验证实,重力滑动作用和上覆地层差异负载作用是导致下刚果盆地盐岩发生塑性流动的主要因素,而基底倾斜程度也对盐构造的变形有重要影响。当基底斜坡较大时,有利于盐体快速流动,引起盐构造快速发育,形成的规模大而数量少,上覆地层断层数量较少;当基底斜坡较小时,有利于盐体缓慢流动,形成盐构造数量多而规模小,上覆地层发育一系列断层。西非板块的隆升作用导致下刚果盆地基底斜坡较陡,进而影响了盐构造变形。