Analogue experiments on the piggyback propagation in northwestern Sichuan and latest propagation in its deeps
-
摘要: 近年来油气勘探揭示了川西前陆北段油气资源的巨大潜力,但因受控于龙门山—川西前陆系统多期构造活动和山前带深埋藏等,川西前陆山前隐伏构造带含油气构造模型解释具有明显的多解性和不确定性。基于川西北前陆原型与实验室砂箱构造物理模型间几何学—运动学—动力学相似性原理,通过两组川西前陆北段扩展变形构造物理模型对比实验(即标准模型实验和断坪—断坡模型实验),揭示川西冲断褶皱带—前陆盆地受控于中、下三叠统膏盐和下寒武统泥岩2套主滑脱层系,晚中生代—新生代呈现出分层式前陆扩展变形特征,基底断坪—断坡结构对川西北前陆深层冲断—冲起构造的控制影响作用最为明显。进一步结合川西北前陆双鱼石地区三维地震解释资料,揭示川西北前展式扩展晚期变形的古生界隐伏冲断—冲起构造构成了川西北深层主要的勘探潜力区。Abstract: Recent explorations reveal a great potential for oil and gas resources in the northwestern Sichuan foreland basin. However, constrained by multi-phase tectonic activities of the Longmenshan-West Sichuan foreland system and the deep burial of the piedmont belt, there remain many uncertainties in the interpretation of petroleum structure models in the foreland concealed tectonic belt. Based on the geometry-kinematics-dynamic similarity theory between the foreland prototype of the northwestern Sichuan and the sandbox analogue model, two groups of controlled experiments were carried out (including standard scaled experiments and analogue experiments with a ramp-flat structure). The fold-thrust belt and foreland basin of the western Sichuan are controlled by two sets of main detachments of the Middle-Lower Triassic gypsum salt and the Lower Cambrian mudstone, leading the layered-style propagation to the foreland in the Late-Middle Cenozoic. The ramp-flat structure plays a profound influence on deep thrust and pop-up structures in the northwestern Sichuan foreland basin. Analogue results combined with seismic interpretation of the northwestern Sichuan illustrates that the main exploration potential may rely on the Paleozoic blind thrust and pop-up structures caused by the latest propagation to the northwestern Sichuan foreland basin.
-
基于几何学—运动学相似性和“实验—实例”互证过程的砂箱构造物理模型模拟,能够最大限度地再现自然界地质结构中潜在的复杂性和不确定性,尤其是多边界条件的物理模型,如:结构机制层模型、应力—应变叠加模型和同构造剥蚀—沉积模型等[1-2]。砂箱构造物理模型模拟和断层相关褶皱端元模型等方法相结合,已广泛地应用于冲断褶皱带—前陆盆地系统等结构—构造特征研究,它们能够有效解释自然界冲断带—前陆盆地系统浅表构造变形过程及其相关含油气系统过程[3-5]。
扬子板块西缘晚三叠世以来形成的龙门山造山带—川西前陆盆地系统既是中国大陆西部中—新生代冲断构造和陆内前陆盆地的典型代表,又是四川盆地天然气勘探最重要的区块之一。前陆盆地形成演化是四川盆地油气基础地质研究的重要内容[6-9],从造山带向前陆盆地的前展式扩展变形特征与川西油气聚集有紧密联系[10]。但限于地层埋深大、钻探难等复杂条件,川西油气勘探成果主要集中于上三叠统须家河组和侏罗系等几套浅层陆相储集体中,如:中坝气田和平落坝气田。前陆盆地陆相地层天然气探明储量超过1.0×1012 m3 [11],近年川西北双鱼石地区深层下二叠统栖霞组(完钻7口井、获工业气井5口、平均单井测试产量56.67×104 m3/d)和泥盆系观雾山组的突破[12],揭示川西北段深层海相层系巨大的油气勘探潜力。但当前川西北段油气勘探仍存在诸多问题,如:多相态的油气显示特征与构造分带性、多期构造演化与油气保存条件、复杂构造与圈闭落实等。
因此,本文基于构造物理模型模拟手段,以川西北段造山带—前陆盆地系统为研究目标,模拟川西前陆盆地山前带隐伏带构造形成演化过程,为川西北山前带双鱼石地区古生界层系含油气构造模型解释及其油气成藏过程研究提供模型参考。
1. 川西北段区域地质特征
晚三叠世松潘—甘孜褶皱带SE向逆冲推覆与川西前陆盆地形成了龙门山造山带—川西前陆盆地系统,导致龙门山冲断带自北西向南东发育4条主干断裂:茂县—汶川断裂、北川—映秀断裂、安县—灌县断裂和广元—大邑(隐伏)断裂。通常以北川—映秀断裂为界将龙门山冲断带划分为龙门山后山带和前山带,将安县—灌县—马角坝断裂和广元—大邑(隐伏)断裂之间称为前陆扩展变形带或山前带。龙门山—川西前陆系统沿主干断裂倾向由北西向南东具明显构造变形层次渐浅、变形强度减弱、卷入层位变新的趋势,总体上呈从冲断带到前陆方向的扩展分带性特征;走向上基底性质及展布、地层发育和活动构造等具差异性走向分带性特征[13-15]。川西北段通常指川西北川—安县一带以北至广元地区的龙门山—川西前陆盆地系统的北段(图 1)。
川西北段通常由凉水—南坝断裂和马角坝断裂分为以具浅变质特征的志留系—元古代层系后山带(常简称为轿子顶推覆体)和以古生界碳酸盐岩层系为主的前山带(常简称为唐王寨推覆体)。川西北段山前隐伏断层常被上覆侏罗系红层覆盖而断续出露,分割前陆扩展变形带(山前带)和前陆盆地;前者为寒武系—三叠系组成、发育强烈断裂和褶皱体系,明显区别于前陆盆地地区的近水平展布的侏罗系—白垩系陆相红层。受控于龙门山—川西前陆盆地系统晚三叠世和新生代两期NW向SE扩展变形过程[16-17],在山前带—前陆盆地区形成沿NE—SW走向展布的冲断褶皱构造,如:矿山梁、天井山、青林口、中坝和双鱼石等构造。由于前陆冲断带多期活动性,沿山前带不同层系出露沥青(如:矿山梁、碾子坝等下寒武统和天井山泥盆系黑色沥青)、液态油(如:天井山泥盆系、青林口侏罗系等油砂)和浅层凝析油—稠油(如:长江沟飞仙关组原油、厚坝侏罗系稠油)等。
川西前陆北段乃至大部分四川盆地区,受加里东期构造运动相关剥蚀作用导致区域泥盆系、石炭系和志留系总体上大量缺失,其余古生界—中生界层系在本地区出露较全,尤其是近期基于盆缘地表—盆内地球物理资料联合研究揭示出的川西早古生界拉张槽及其相关下寒武统巨厚深水沉积成因的筇竹寺组黑色泥岩层系[18]。早—中三叠世受控于碳酸盐岩台地相沉积,川西地区沉积了巨厚的泥岩—膏盐层系(即飞仙关组泥岩、嘉陵江—雷口坡组膏盐层系);随后川西前陆盆地进入陆相前陆盆地演化阶段。受控于上述2套滑脱层系和冲断带前展式扩展变形过程,川西前陆盆地发育复杂构造样式,如:隐伏堆垛背斜构造/叠置背斜、叠瓦状冲断、断层转折/传播褶皱和双重构造等,且存在深、浅部滑脱分割形成的不同构造样式。川西北双鱼石地区三维地震资料解释,揭示出本地区具巨厚下寒武统筇竹寺组黑色页岩层系和中—下三叠统膏盐滑脱层系,构造样式上具有典型的前展式冲断构造与冲起构造样式,尤其古生界层系发育典型的反向冲断层,它们与前展式冲断层局部构成断层三角带等。
2. 砂箱物理模拟模型设置
砂箱物理模型基于相似性原理模拟冲断褶皱带—前陆盆地系统浅表构造变形作用过程,即它们普遍被认为遵循Mohr-Coulomb破裂准则[19-21]。因此实验中采用石英砂建模,它们普遍具有中等磨圆、较好分选性,平均粒度和密度分别为0.2~0.4 mm和1.55 g/cm3,内摩擦角和摩擦系数分别为29°~31° 和0.58[22-23],且石英砂普遍具有近似线性的库伦—摩尔破裂变形行为和近似为零的内聚力值(20~100 Pa),能够有效模拟0~10 km浅表地壳构造变形过程[24-26]。实验过程中,分别使用玻璃珠和硅胶代表自然界泥岩(下寒武统筇竹寺组泥岩、下志留统龙马溪组泥页岩)和膏盐滑脱层系(中—下三叠统膏盐层系)[27]。为了获得合理的实验模拟结果,需要遵循实验室模型和自然界实际构造模型几何学、运动学和动力学的相似性[2, 24]。因此,需结合川西前陆地区实际地层展布特征和区域结构—构造特征设计具体物理模型垂向物质比例系数和结构边界条件(图 2)。
本次共设计2组川西北段砂箱物理模型模拟实验,即(1)实际地层比例模型标准实验和(2)断坪—断坡结构模拟实验。地层比例模型实验为标准化实验基础,以川西北双鱼石地区平均地层厚度结构为几何学标准,采用的几何学比例系数为5×10-6,即砂箱模型中1 cm代表自然界中2 km厚度,模型底部、中部和上部分别加入1~2 mm、2~3 mm、3~4 mm和2~3 mm玻璃珠或硅胶代替垂向地层中滑脱层系。断坪—断坡结构模型实验中与标准化实验仅2个边界条件不同,(1)基底加入4~5 mm厚度玻璃珠代替现今下寒武统牛蹄塘组页岩层系和(2)基底厚度4~5 mm的树脂面板代替现今下寒武统底部上超界面“陡坡带”(为基底滑脱变形的应变集中带),断坡角为国际常用的45°[28-29]。实验都在成都理工大学“油气藏地质与开发工程”国家重点实验室构造物理模拟实验室开展。均质石英砂物质通过匀速筛网自动铺设并分步骤间序加入玻璃珠和硅胶,形成长宽高为850 mm×345 mm×39 mm砂体,且石英砂中铺设约1 mm厚彩色石英砂作为标志层以便于构造变形过程监测和量化对比等(图 2)。为减小砂箱物理模型中石英砂均质体周围玻璃隔板的边界效应,铺设石英砂前使用有机去泡剂润滑模型周围玻璃隔板。
重力条件下应变比例系数能够通过物理模型与自然界模型比例系数来表达,应变比例系数σ *=ρ *L *g *, 其中ρ *、L *、g *分别为物理模型和自然模型的密度、长度和重力加速度间比例系数,因此计算得到,σ *=6.6×10-6。若自然界岩石剪切强度为5~20 MPa,物理模型内聚力强度约为20~100 Pa,它们普遍适用于正常重力条件下比例砂箱物理模型。本次实验模型中采用挤压速率为0.003 mm/s,代表自然界中大致缩短速率为7.1 km/Myr,进行模型模拟对比实验。
表 1 砂箱物理模型模拟实验比例系数特征Table 1. Scale coefficient for sandbox physical model of simulation experiments参数 物理模型(M) 自然实例(N) 比例系数(M/N) 长度(L) /m 8.5×10-1 1.7×105 L*=5×10-6 缩短量(D)/m 0.30 6.0×104 D*=5×10-6 内聚力 20~100 Pa 5~20 MPa 重力加速度(g)/(m·s-2) 9.81 9.81 g* =1.0 密度(石英砂)/(kg·m-3) 1 550 2 200 ρ*1=0.70 密度(玻璃珠)/(kg·m-3) 1 370 2 400 ρ*2=0.57 密度(硅胶)/(kg·m-3) 960 2 300 ρ*3=0.42 应力(σ) σ* =ρ*g*L* =6.60×10-6 动力黏度(η)/(Pa·s) 2.1×104 1.7×1019 η* =1.24×10-15 时间 1×105 s tN=tM/t*=8.44 Ma t* =η*/σ* =3.76×10-10 速率 0.003 mm/s 7.1 km/Myr 注:tN,tM分别代表自然实际和实验模型中的时间跨度,t*为实验—实例时间比例系数,η*为动力黏度比例系数。表中部分参数来源于文献[23],计算公式参考文献[4]。 每组砂箱物理模型实验中,通过砂箱装置左侧活动挡板恒定速率挤压砂箱物质产生缩短变形。挤压缩短变形过程中,使用同步数码相机每间隔1.0 mm缩短量记录物质变形过程,挤压缩短完成后使用图像编辑软件每间隔5 mm缩短量量化测量楔形体几何参数并进行对比研究,如:楔高、楔长、楔顶角等。
3. 砂箱物理模型模拟结果
3.1 川西北段地层比例模型标准实验结果
持续挤压过程中,初始挤压缩短量为40 mm时,砂箱物质发生缩短变形、形成典型的冲起结构,由2条前展式逆冲断层及其相关反向冲断层构成。伴随挤压缩短量增大,逐渐形成多条新生的前展式逆冲断层,如:T3、T4逆断层(图 3a中的②和③),且楔形体后缘早期断层倾角伴随缩短发生大量的旋转变形。在缩短量达到140 mm时,受后缘挤压缩短变形砂箱物质沿浅部滑脱层系发生明显的滑脱扩展变形,形成浅部逆断层T5a,而深部石英砂物质未发生明显的冲断变形特征(图 3a中的④)。进一步缩短变形至缩短量为155 mm时,深部石英砂物质发生逆冲变形、形成逆断层T5b,其具有明显的下部断坡、中部断坪和上部断坡结构。上下部断坡结构分别切割砂箱物质,中部断坪沿硅胶滑脱层系发生明显的顺层滑脱变形,最终贯通连接形成T5断层及其相关双重构造,双重构造底板断层为基底断层、顶板断层为滑脱逆冲断层构成。随后挤压缩短变形过程导致浅部滑脱层系(即硅胶滑脱层)进一步顺层缩短扩展变形,依次形成浅部前展式逆断层T6a、T7a。当缩短量达到230 mm时,深部石英砂物质发生明显的逆冲断层变形、形成T6b及其相应的下部断坡、中部滑脱层系断坪、上部断坡结构组合样式,最终形成T6逆冲断层与双重构造。至最后挤压缩短过程中,与前述过程相似砂箱物质依次形成浅部滑脱逆冲断层T8a(缩短量为260 mm时)、和深部T7b逆冲断层及其相关变形结构,且由于后缘强烈挤压缩短变形导致多个断片叠置,形成背斜状叠置双重构造;同时,楔形体后缘早期断层由于挤压缩短相关的物质旋转形成“似花状结构”。
伴随持续挤压缩短过程,砂箱物质变形过程可以分为2个典型阶段,即前展式扩展变形阶段和双重构造扩展变形阶段。前展式扩展变形阶段中,楔形体后缘挤压缩短导致砂箱物质浅部和深部发生一致的构造变形,以从冲断带到前陆盆地的冲断扩展变形过程为主;楔形体楔高和楔长伴随砂箱物质前缘新生断层形成而快速生长,且楔长具有与新生断层形成相一致的跃迁性和阶段性(图 4),每当新生断层形成后楔长伴随挤压缩短具有逐渐减小的特征,与砂箱物质明显的平行层挤压缩短变形特征相一致,如:当缩短量为40 mm时形成新生T2逆断层、楔长从75 mm快速生长至95 mm,随后挤压缩短过程中楔长逐渐减小至约80 mm(对应缩短量约为65 mm时)。双重构造扩展变形阶段中,伴随挤压缩短砂箱物质首先沿浅部滑脱层系发生前展式逆冲与扩展变形,形成浅部滑脱逆冲断层,随后发生深部逆冲断层及其相关断坪—断坡结构,形成典型的双重构造。由于砂箱物质浅部和深部先后分别发生前展式扩展过程,虽然本阶段楔形体楔长生长仍然与浅部滑脱逆冲断层新生具有一致性,呈现出阶段性生长过程,但其楔形体高度生长速率明显较早期阶段缓慢(图 4)。
3.2 川西北段断坪—断坡结构模型模拟结果
川西北断坪—断坡结构模型挤压缩短变形过程大致也可分为2个演化阶段,即前展式扩展变形阶段和双重构造扩展变形阶段(图 3b)。伴随逐渐挤压缩短变形,当缩短量为50 mm时楔形体物质形成第一个冲起构造,由2条前展式逆冲断层及其相关反向冲断层构成。伴随持续挤压,形成第三条前展式逆断层T3。随后进入双重构造扩展变形阶段,伴随挤压缩短过程首先在砂箱浅表沿滑脱层系形成典型的前展式滑脱逆冲断层,如:D=90 mm时,形成T4a逆断层;D=140 mm时,形成T5a逆断层等(图 4)。随后进一步挤压缩短导致砂箱物质深部石英砂物质发生逆冲缩短构造变形作用,如:D=70 mm时,形成T3b逆断层;D=110 mm时,形成T4b逆断层;D=170 mm时,形成T5b逆断层等(图 3b的④和⑤)。同时深部逆冲断层构成下部断坡结构,连通中部沿滑脱层系的断坪结构和上部断坡结构(即早期的滑脱逆冲断层),构成典型的双重构造。
当挤压缩短量为190 mm时,砂箱浅表物质发生滑脱缩短变形过程、形成浅部T6a逆冲断层。楔形体后缘进一步挤压缩短导致深部石英砂物质形成形态宽缓的冲起构造,该构造分别由深部T6b前展式逆断层和T6b'反冲逆断层构造构成。深部冲起构造及其反向逆冲断层的形成,可能与深部物质变形前缘带的断坪—断坡结构控制形成反向挤压应力密切相关。随后进一步挤压缩短导致冲起构造发生构造叠加变形,导致其早期冲起结构发生强烈叠加改造并卷入了新的双重构造断片中(图 3b)。
与标准实验模型模拟结果相似,断坪—断坡结构模型中楔形体高度生长大致也分为早期快速生长阶段和晚期缓慢生长阶段,同时楔形体长度(即楔长)伴随楔形体前缘新生逆冲断层形成呈现出阶段性跃迁生长特征(图 4)。至最终挤压缩短量为300 mm时,标准模型和断坪—断坡结构模型楔长、楔形体高度和楔顶角等特征总体一致,揭示出断坪—断坡结构虽然导致砂箱物质深部变形过程与构造特征存在重要差别,但对楔形体浅表几何学特征作用不明显。
4. 川西北前陆冲断构造特征及其物理模拟意义
4.1 分层式扩展变形过程及其双重构造特征
冲断褶皱带挤压变形过程中通常形成3种典型的扩展变形过程模式(图 5),即“背驮式”扩展变形模式(或前展式扩展变形模式)、后陆式扩展变形模式和介于二者之间的无序冲断扩展变形模式[30],常形成典型的叠瓦状冲断岩片。前展式扩展变形模式强调受造山带向盆地方向挤压变形作用控制,新生断层常发育于早期断层下盘,成“背驮式”依次向前陆盆地冲断扩展,导致前陆盆地方向断层形成活动年龄变新。后陆式扩展变形模式则与前者相反,强调新生断层通常在早期断层的上盘形成,逐渐向后陆地区冲断扩展、导致往盆地方向断层年龄变老,常常发育在以重力扩展成因为主的构造变形地区[31],它通常会导致基底主滑脱层及其前缘层系构造变形停止;而无序冲断扩展变形过程则介于上述两种端元之间,没有明显的扩展变形活动序列过程。
川西前陆盆地构造物理模型中(即标准实验模型和断坪—断坡结构模型中),伴随楔形体后缘持续挤压缩短过程,前展式扩展变形过程分别发生在基底滑脱层系和浅部滑脱层系,从而体现出分层式扩展冲断变形特征。在楔形体挤压生长后期(临界楔后期)通常率先发生浅部滑脱层系的前展式滑脱扩展冲断变形(形成浅部叠瓦状冲断结构),随后再发生深部层系的前展式滑脱扩展冲断变形(图 5d),从而在冲断层前缘深部(如:山前隐伏断层带地区等)形成典型的双重构造,其前缘断层普遍存在多期活动特征,而明显区别于上述3种扩展变形模式及其相关叠瓦状冲断带特征。结合川西北地区地震剖面结构解析,也揭示出前陆盆地山前隐伏带深层区域存在广泛的双重构造,如:矿山梁双重构造、天井山双重构造等[9, 17, 32],其顶底板逆冲断层应分别沿中—下三叠统膏盐层系和下寒武统筇竹寺组泥岩层系滑脱。BONINI通过多滑脱层系构造物理模型模拟,强调受控于基底、上部滑脱层系和砂箱物质间强度比及其厚度[33-34],通常会形成(1)从冲断带到前陆盆地扩展传播的双重构造和(2)被动顶板冲断双重构造样式,而前陆盆地区同构造沉积作用也进一步加剧形成被动顶板双重构造的趋势。与之相似,NILFOUROUSHAN等[35]通过双层滑脱属性数字模型模拟实验中也得到相似的分层式扩展变形演化过程及其相关双重构造,但他更强调物质材料本身属性(即内聚力大小)对于多层滑脱模型与变形过程的控制作用。
分层式前陆扩展变形模式及其形成演化过程也表明砂箱物质和自然界实例中垂向空间格架上可能存在构造变形强度、变形样式和缩短变形量等差异。如:标准模型试验中缩短量为300 mm时其深部简单的双重构造与上覆层系滑脱冲断变形(图 3a的⑨),断坪—断坡构造模型中缩短量为205 mm时,深部简单冲起构造与上覆层系滑脱冲断变形(图 3b的⑦)等。垂向空间变形特征差异性主要受控于浅部滑脱层系的顺层滑脱作用,从而导致山前带深部层系构造变形明显弱于浅部强构造变形,因此具有更加有利的油气保存条件。
4.2 川西北构造带新生代扩展变形构造特征
分层式前陆扩展变形演化过程强调冲断褶皱带—前陆盆地中浅部滑脱冲断变形过程明显早于深部逆冲断层构造变形活动过程,冲断带后缘持续挤压过程导致深部层系发生前展式冲断扩展变形,形成典型深部断坡结构,最终连接顺滑脱层的断坪和浅部断坡结构,形成典型的双重构造;导致冲断带系统持续向前陆盆地扩展变形生长,其砂箱物质深部层系构造样式通常为楔形体前展式扩展最晚期构造变形过程中形成。晚三叠世以来,龙门山造山带—川西前陆盆地系统在印支期、燕山期和新生代喜玛拉雅期持续发生多期挤压构造变形,导致冲断带持续由NW向SE扩展变形[7, 14, 36]。因此,川西前陆地区山前隐伏带中—下三叠统滑脱层之下的古生界层系构造样式应为龙门山冲断带晚期前展式扩展变形过程中形成,结合龙门山地区新生代显著的构造缩短特征和现今强烈的地震活动等[37],我们认为川西北双鱼石构造带深层古生界结构样式及其相关变形为晚期新生代构造成因。
标准实验模型和断坪—断坡结构模型实验对比,揭示出后者深部砂箱物质受控于基底滑脱层和断坡结构的反向作用过程,从而形成典型的深部冲起构造和相关冲断层等(图 3b)。DENG等[23]通过砂箱物理模型模拟,揭示出低速挤压缩短变形过程及其相关典型的对称性冲起构造特征。同时NILFOUROUSHAN等[35]也强调深部冲起结构带发育的位置及其规模受控于浅部滑脱层系与下伏层系间的脱耦作用,伴随砂箱物质属性内聚力减小,深部层系更加容易形成冲起构造带。川西前陆盆地早寒武世受控于兴凯地裂运动,沉积了巨厚的下寒武统筇竹寺组泥岩层系,它既是川西地区最优质的烃源岩之一,又是典型基底滑脱层系[18]。受控于龙门山冲断带向前陆盆地扩展变形过程,沿基底滑脱层系发生晚期滑脱扩展变形过程,滑脱层系越厚、滑脱扩展变形作用越强,其上覆地层卷入变形形成冲起构造规模影响越显著。对川西北前陆区的地震测线构造解释,揭示出在下伏包含巨厚筇竹寺组泥岩层系的区域普遍发育典型的冲起构造(图 6),且部分冲起构造主断层滑脱归并于下伏基底滑脱面,它们充分表明区域晚期扩展变形作用的重要性。值得一提的是,川西北双鱼石构造冲断前缘深部隐伏了绵阳—长宁拉张槽的东部边界,现今的构造解释表明,拉张槽边界(几何形态可充当基底断坡)对上覆地层中—新生代以来的前展式扩展变形有重要控制作用,主要体现在断坡位置附近发育由NW向前展式逆冲断层和SE向反冲断层构成的冲起构造和构造三角带,不同于两侧仅发育单向的逆冲断层,钻获高产工业气流的ST1井和ST3井均钻于该区构造高点(图 6)。
图 6 川西北前陆区地震测线及其深部冲起结构构造解释对比测线位置见图 1。Figure 6. Seismic lines of the foreland zone in northwestern Sichuan Basin and the interpretation of their thrust and pop-up structures4.3 山前带古生界深层油气勘探潜力
川西北前陆地区古生界层系主要为海相地层,具有良好的源—储油气地质条件,前陆盆地山前带古生界层系不仅地表出露大量古油藏和油苗,且潜伏构造钻探也有较好的油气显示,如:天井山、矿山梁和河湾场等构造。由于龙门山冲断带多期前展式冲断扩展变形,导致山前带隐伏断层上盘地层(包括古生界—中生界层系)普遍发生强烈缩短变形、抬升剥蚀出露,如:香水断层等,封盖保存条件是成藏关键因素。砂箱构造物理模型模拟揭示出,冲断带前缘断层受控于分层式扩展变形演化过程,普遍具有多期活动性,如:T5a和T5b、T6a和T6b逆断层等,因而山前带主断层上盘目标层系难以具有良好的封盖保存条件,这与现今川西北浅层大量钻井和地表沥青、油苗显示等特征相一致。同时,山前隐伏带半原地系统由于扩展变形可能形成典型的双重构造,如图 3a和3b的⑨,半原地系统中双重构造,其古生界—中生界目的层系则相对具有利的油气聚集和封盖保存条件,如:中坝构造等。下伏深层层系的晚期构造变形常常导致山前带主断层上盘层系构造改造,从而使早期古油气藏遭破坏,发生油气再运聚,它是川西北段山前带大量古油藏和油苗出露的根本原因。
受龙门山造山带新生代扩展变形作用控制,川西北前陆盆地区古生界层系主体构造形成期较晚,同时其上覆中—下三叠统膏盐层系作为盆地区域盖层具较好封盖保存条件。因此,川西北前陆区深层古生界圈闭成藏的关键是能否有效捕获和聚集油气。由于川西地区生油气层较多,如:陡山沱组黑色泥岩、筇竹寺组黑色泥岩等,我们认为川西前陆区晚期形成的原地构造系统具有较佳的油气勘探前景;尤其是古生界层系中发育大量构造圈闭,部分断层能够垂向沟通多个有效烃源岩层,具有形成大中型油气聚集带的天然优势条件。目前川西北双鱼石地区下二叠统栖霞组多口探井平均测试产量高达56.67×104 m3/d,初步展示了川西前陆区原地系统古生界层系的良好勘探前景。
5. 结论
(1) 龙门山冲断带—川西前陆盆地系统砂箱构造物理模型模拟揭示,伴随持续挤压缩短过程砂箱物质变形可以大致分为前展式扩展变形阶段和双重构造扩展变形两阶段。早期阶段以砂箱物质浅部和深部发生一致的前展式冲断扩展变形为主;晚期双重构造扩展变形阶段,首先沿浅部滑脱层系发生前展式扩展变形与逆冲变形,形成浅部叠瓦状滑脱逆冲断层,随后再发生深部逆冲断层并发育相关断坪—断坡结构,最终形成深部隐伏双重构造,整体上表现为分层式扩展变形特征。
(2) 分层式扩展变形模式揭示出龙门山川西前陆区原地系统隐伏构造晚期新生代构造扩展变形成因机制及其相对较弱的构造变形特征。深部物质受控于基底滑脱和断坪—断坡结构形成前陆盆地深部古生界层系冲起构造和构造三角带,它们构成了川西北深层古生界最有利的油气圈闭,是目前川西北深层勘探的潜力区。
-
图 6 川西北前陆区地震测线及其深部冲起结构构造解释对比
测线位置见图 1。
Figure 6. Seismic lines of the foreland zone in northwestern Sichuan Basin and the interpretation of their thrust and pop-up structures
表 1 砂箱物理模型模拟实验比例系数特征
Table 1. Scale coefficient for sandbox physical model of simulation experiments
参数 物理模型(M) 自然实例(N) 比例系数(M/N) 长度(L) /m 8.5×10-1 1.7×105 L*=5×10-6 缩短量(D)/m 0.30 6.0×104 D*=5×10-6 内聚力 20~100 Pa 5~20 MPa 重力加速度(g)/(m·s-2) 9.81 9.81 g* =1.0 密度(石英砂)/(kg·m-3) 1 550 2 200 ρ*1=0.70 密度(玻璃珠)/(kg·m-3) 1 370 2 400 ρ*2=0.57 密度(硅胶)/(kg·m-3) 960 2 300 ρ*3=0.42 应力(σ) σ* =ρ*g*L* =6.60×10-6 动力黏度(η)/(Pa·s) 2.1×104 1.7×1019 η* =1.24×10-15 时间 1×105 s tN=tM/t*=8.44 Ma t* =η*/σ* =3.76×10-10 速率 0.003 mm/s 7.1 km/Myr 注:tN,tM分别代表自然实际和实验模型中的时间跨度,t*为实验—实例时间比例系数,η*为动力黏度比例系数。表中部分参数来源于文献[23],计算公式参考文献[4]。 -
[1] BUTLER R W H, BOND C E, COOPER M A, et al. Interpreting structural geometry in fold-thrust belts: why style matters[J]. Journal of Structural Geology, 2018, 114: 251-273. doi: 10.1016/j.jsg.2018.06.019 [2] GRAVELEAU F, MALAVIEILLE J, DOMINGUEZ S. Experimental modelling of orogenic wedges: a review[J]. Tectonophysics, 2012, 538-540: 1-66. doi: 10.1016/j.tecto.2012.01.027 [3] DUERTO L, MCCLAY K. The role of syntectonic sedimentation in the evolution of doubly vergent thrust wedges and foreland folds[J]. Marine and Petroleum Geology, 2009, 26(7): 1051-1069. doi: 10.1016/j.marpetgeo.2008.07.004 [4] DOOLEY T P, JACKSON M P A, HUDEC M R. Breakout of squeezed stocks: dispersal of roof fragments, source of extrusive salt and interaction with regional thrust faults[J]. Basin Research, 2015, 27(1): 3-25. doi: 10.1111/bre.12056 [5] 邓宾, 黄瑞, 马华灵, 等. "从源到汇": 褶皱冲断带-前陆盆地系统砂箱物理模型浅表作用研究[J]. 大地构造与成矿学, 2018, 42(3): 431-444. https://www.cnki.com.cn/Article/CJFDTOTAL-DGYK201803002.htmDENG Bin, HUANG Rui, MA Hualing, et al. "Source to sink": a review of tectonics-erosion-sedimentation interactions in the analogue experiments of accretionary wedge[J]. Geotectonica et Metallogenia, 2018, 42(3): 431-444. https://www.cnki.com.cn/Article/CJFDTOTAL-DGYK201803002.htm [6] 刘树根, 罗志立, 赵锡奎, 等. 试论中国西部陆内俯冲型前陆盆地的基本特征[J]. 石油与天然气地质, 2005, 26(1): 37-48. https://www.cnki.com.cn/Article/CJFDTOTAL-SYYT200501006.htmLIU Shugen, LUO Zhili, ZHAO Xikui, et al. Discussion on essential characteristics of intracontinental-subduction type foreland basins in western China[J]. Oil & Gas Geology, 2005, 26(1): 37-48. https://www.cnki.com.cn/Article/CJFDTOTAL-SYYT200501006.htm [7] 刘和甫, 梁慧社, 蔡立国, 等. 川西龙门山冲断系构造样式与前陆盆地演化[J]. 地质学报, 1994, 68(2): 101-118. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXE199402000.htmLIU Hefu, LIANG Huishe, CAI Liguo, et al. Structural styles of the Longmenshan thrust belt and evolution of the foreland basin in western Sichuan Province, China[J]. Acta Geologica Sinica, 1994, 68(2): 101-118. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXE199402000.htm [8] 杨克明. 川西坳陷须家河组天然气成藏模式探讨[J]. 石油与天然气地质, 2006, 27(6): 786-793. https://www.cnki.com.cn/Article/CJFDTOTAL-SYYT200606010.htmYANG Keming. Gas reservoiring mode in Xujiahe Formation of Western Sichuan Depression[J]. Oil & Gas Geology, 2006, 27(6): 786-793. https://www.cnki.com.cn/Article/CJFDTOTAL-SYYT200606010.htm [9] JIA Dong, WEI Guoqi, CHEN Zhuxin, et al. Longmen Shan fold-thrust belt and its relation to the western Sichuan Basin in central China: new insights from hydrocarbon exploration[J]. AAPG Bulletin, 2006, 90(9): 1425-1447. doi: 10.1306/03230605076 [10] 徐旭辉, 方成名, 刘金连, 等. 中国中西部山前构造变形结构分带模式与油气[J]. 石油实验地质, 2019, 41(6): 779-790. doi: 10.11781/sysydz201906779XU Xuhui, FANG Chengming, LIU Jinlian, et al. Deformation zoning model of piedmont thrust, western China, and its petroleum response[J]. Petroleum Geology & Experiment, 2019, 41(6): 779-790. doi: 10.11781/sysydz201906779 [11] 刘树根, 孙玮, 王国芝, 等. 四川叠合盆地油气富集原因剖析[J]. 成都理工大学学报(自然科学版), 2013, 40(5): 481-497. https://www.cnki.com.cn/Article/CJFDTOTAL-CDLG201305001.htmLIU Shugen, SUN Wei, WANG Guozhi, et al. Analysis of causes of oil and gas accumulation in superimposed Sichuan Basin of China[J]. Journal of Chengdu University of Technology (Science & Technology Edition), 2013, 40(5): 481-497. https://www.cnki.com.cn/Article/CJFDTOTAL-CDLG201305001.htm [12] 张本健, 谢继容, 尹宏, 等. 四川盆地西部龙门山地区中二叠统碳酸盐岩储层特征及勘探方向[J]. 天然气工业, 2018, 38(2): 33-42. https://www.cnki.com.cn/Article/CJFDTOTAL-TRQG201802006.htmZHANG Benjian, XIE Jirong, YIN Hong, et al. Characteristics and exploration direction of the Middle Permian carbonate reservoirs in the Longmenshan Mountain areas, western Sichuan Basin[J]. Natural Gas Industry, 2018, 38(2): 33-42. https://www.cnki.com.cn/Article/CJFDTOTAL-TRQG201802006.htm [13] 刘树根, 李智武, 曹俊兴, 等. 龙门山陆内复合造山带的四维结构构造特征[J]. 地质科学, 2009, 44(4): 1151-1180. https://www.cnki.com.cn/Article/CJFDTOTAL-DZKX200904012.htmLIU Shugen, LI Zhiwu, CAO Junxing, et al. 4-D textural and structural characteristics of Longmen intracontinental composite orogenic belt, southwest, China[J]. Chinese Journal of Geology, 2009, 44(4): 1151-1180. https://www.cnki.com.cn/Article/CJFDTOTAL-DZKX200904012.htm [14] LIU Shugen, DENG Bin, LI Zhiwu, et al. Architecture of basin-mountain systems and their influences on gas distribution: a case study from the Sichuan Basin, South China[J]. Journal of Asian Earth Sciences, 2012, 47: 204-215. doi: 10.1016/j.jseaes.2011.10.012 [15] DENG Bin, LIU Shugen, JANSA L, et al. Sedimentary record of Late Triassic transpressional tectonics of the Longmenshan thrust belt, SW China[J]. Journal of Asian Earth Sciences, 2012, 48: 43-55. doi: 10.1016/j.jseaes.2011.12.019 [16] 胡烨, 陈迎宾, 王彦青, 等. 川西坳陷回龙构造雷口坡组天然气成藏条件[J]. 特种油气藏, 2018, 25(1): 46-51. https://www.cnki.com.cn/Article/CJFDTOTAL-TZCZ201801010.htmHU Ye, CHENYingbin, WANG Yanqing, et al. Naturalgas accumulation conditions in LeiKoupo Formation of Huilongstructure, Western Sichuan Depression[J]. Special Oil & Gas Reservoirs, 2018, 25(1): 46-51. https://www.cnki.com.cn/Article/CJFDTOTAL-TZCZ201801010.htm [17] 陈竹新, 贾东, 魏国齐, 等. 龙门山北段冲断前锋构造带特征[J]. 石油学报, 2008, 29(5): 657-662. https://www.cnki.com.cn/Article/CJFDTOTAL-SYXB200805006.htmCHEN Zhuxin, JIA Dong, WEI Guoqi, et al. Characteristics of thrust structures in the northern Longmenshan front belt[J]. Acta Petrolei Sinica, 2008, 29(5): 657-662. https://www.cnki.com.cn/Article/CJFDTOTAL-SYXB200805006.htm [18] LIU Shugen, DENG Bin, JANSA L, et al. The Early Cambrian Mianyang-Changning intracratonic sag and its control on petroleum accumulation in the Sichuan Basin, China[J]. Geofluids, 2017, 2017: 6740892. [19] CHAPPLE W M. Mechanics of thin-skinned fold-and-thrust belts[J]. GSA Bulletin, 1978, 89(8): 1189-1198. [20] DAVIS D, SUPPE J, DAHLEN F A. Mechanics of fold-and-thrust belts and accretionary wedges[J]. Journal of Geophysical Research: Solid Earth, 1983, 88(B2): 1153-1172. [21] DAHLEN F A, SUPPE J, DAVIS D. Mechanics of fold-and-thrust belts and accretionary wedges: cohesive coulomb theory[J]. Journal of Geophysical Research: Solid Earth, 1984, 89(B12): 10087-10101. [22] 邓宾, 赵高平, 万元博, 等. 褶皱冲断带构造砂箱物理模型研究进展[J]. 大地构造与成矿学, 2016, 40(3): 446-464. https://www.cnki.com.cn/Article/CJFDTOTAL-DGYK201603004.htmDENG Bin, ZHAO Gaoping, WAN Yuanbo, et al. A review of tecto-nic sandbox modeling of fold-and-thrust belt[J]. Geotectonica et Metallogenia, 2016, 40(3): 446-464. https://www.cnki.com.cn/Article/CJFDTOTAL-DGYK201603004.htm [23] DENG Bin, JIANG Lei, ZHAO Gaoping, et al. Insights into the velocity-dependent geometry and internal strain in accretionary wedges from analogue models[J]. Geological Magazine, 2018, 155(5): 1089-1104. [24] BUITER S J H. A review of brittle compressional wedge models[J]. Tectonophysics, 2012, 530-531: 1-17. [25] WEIJERMARS R, JACKSON M P A, VENDEVILLE B. Rheological and tectonic modeling of salt provinces[J]. Tectonophysics, 1993, 217(1/2): 143-174. [26] LOHRMANN J, KUKOWSKI N, ADAM J, et al. The impact of analogue material properties on the geometry, kinematics, and dynamics of convergent sand wedges[J]. Journal of Structural Geology, 2003, 25(10): 1691-1711. [27] KOYI H A, HESSAMI K, TEIXELL A. Epicenter distribution and magnitude of earthquakes in fold-thrust belts: insights from Sandbox Models[J]. Geophysical Research Letters, 2000, 27(2): 273-276. [28] KOYI H A, MAILLOT B. Tectonic thickening of hanging-wall units over a ramp[J]. Journal of Structural Geology, 2007, 29(6): 924-932. [29] BONINI M, SOKOUTIS D, MULUGETA G, et al. Modelling hanging wall accommodation above rigid thrust ramps[J]. Journal of Structural Geology, 2000, 22(8): 1165-1179. [30] MCCLAY K R. Glossary of thrust tectonic terms[M]//Thrust Tectonics. London: Chapman and Hall, 1992: 419-433. [31] ALSOP G I, WEINBERGER R, MARCO S. Distinguishing thrust sequences in gravity-driven fold and thrust belts[J]. Journal of Structural Geology, 2018, 109: 99-119. [32] 管树巍, 李本亮, 何登发, 等. 构造楔形体的识别与勘探: 以准噶尔盆地南缘为例[J]. 地学前缘, 2009, 16(3): 129-137. https://www.cnki.com.cn/Article/CJFDTOTAL-DXQY201901017.htmGUAN Shuwei, LI Benliang, HE Dengfa, et al. Recognition and exploration of structural wedges: a case study of the southern margin of Junggar Basin[J]. Earth Science Frontier, 2009, 16(3): 129-137. https://www.cnki.com.cn/Article/CJFDTOTAL-DXQY201901017.htm [33] BONINI M. Passive roof thrusting and forelandward fold propagation in scaled brittle-ductile physical models of thrust wedges[J]. Journal of Geophysical Research: Solid Earth, 2001, 106(B2): 2291-2311. [34] BONINI M. Deformation patterns and structural vergence in brittle-ductile thrust wedges: an additional analogue modelling perspective[J]. Journal of Structural Geology, 2007, 29(1): 141-158. [35] NILFOUROUSHAN F, PYSKLYWEC R, CRUDEN A. Sensitivity analysis of numerical scaled models of fold-and-thrust belts to granular material cohesion variation and comparison with analog experiments[J]. Tectonophysics, 2012, 526-529: 196-206. [36] 邓涛, 李勇. 龙门山前陆盆地南段须五段页岩气富集条件及有利区评价[J]. 特种油气藏, 2018, 25(4): 51-56. https://www.cnki.com.cn/Article/CJFDTOTAL-TZCZ201804010.htmDENG Tao, LIYong. Shalegas enrichment condition and favorable area evaluation of the Xu5 Member in the southern section of Longmenshan Foreland Basin[J]. Special Oil & Gas Reservoirs, 2018, 25(4): 51-56. https://www.cnki.com.cn/Article/CJFDTOTAL-TZCZ201804010.htm [37] ZHANG Peizhen, WEN Xueze, SHEN Zhengkang, et al. Oblique, high-angle, listric-reverse faulting and associated development of strain: the Wenchuan Earthquake of May 12, 2008, Sichuan, China[J]. Annual Review of Earth and Planetary Sciences, 2010, 38(1): 353-382. 期刊类型引用(7)
1. 金文正,白万奎,王俊鹏,叶治续. 川西北碧口地块中新生代隆升及构造扩展特征. 华南地质. 2023(01): 24-36 . 百度学术
2. 张欢,陈应涛,陶威,陈涛,余文鑫,艾卉卉. 不同拉伸方式和速度下的伸展构造砂箱物理模拟实验研究. 西北地质. 2023(02): 327-336 . 百度学术
3. 王耀华,甘军,梁刚,李兴,杨威,高之业,熊书苓. 断裂-砂体-潜山复式天然气输导体系及成藏模式——以琼东南盆地深水区为例. 断块油气田. 2022(03): 319-324 . 百度学术
4. 徐宏远,孙玮,郭虹兵,邵红君,田小彬. 米市盆地构造物理模拟研究. 中国石油大学胜利学院学报. 2021(01): 11-15 . 百度学术
5. 张连进,黄家强,罗强,何宇,于豪,杨荣军,梁霄,周政,邓宾. 川西北前陆双鱼石地区砂箱物理模拟及其油气地质意义. 地质科技通报. 2021(02): 156-166 . 百度学术
6. 邓宾,何宇,黄家强,罗强,杨荣军,于豪,张静,刘树根. 前陆盆地形成与演化砂箱物理模拟启示——以四川盆地西部龙门山为例. 石油与天然气地质. 2021(02): 401-415 . 百度学术
7. 魏子博,陈雷,冯亮,陈鑫,蔡家兰,闫梦楠,谭秀成. 川东奉节地区长兴组沉积相及储层发育特征. 断块油气田. 2021(06): 816-822 . 百度学术
其他类型引用(1)
-