杨攀新 胡伟华 胡朝忠 孙鑫喆黄帅堂 郭春杉 陈志华 任金卫
1)中国北京 100036 中国地震局地震预测研究所
2)中国乌鲁木齐 830011 新疆维吾尔自治区地震局
3)中国西安 710077 中国煤炭科工集团西安研究院有限公司
4)中国广州 510800 广东省地质物探工程勘察院
印度板块和欧亚板块的碰撞使得欧亚大陆中部产生剧烈变形,从而导致形成于古生代末的天山褶皱带在晚新生代以后再次强烈隆升成再生造山带(Allenet al,1993;邓起东等,2000).由于青藏高原西构造结节向北推挤和塔里木地块顺时针旋转导致了西天山和哈萨克斯坦地块内部发育一组NW向右旋走滑断层(图1).塔拉斯—费尔干纳断裂是中亚地区的巨型右旋走滑断裂带,始于哈萨克斯坦境内,终至中国喀什以北,全长约为1 000 km,累计位移达200 km (Burtmanet al,1996),研究表明该断裂上距今5 000年内发生了三次古地震事件,全新世以来右旋滑动走滑速率为13—14 mm/a (Korjenkovetal,2010),现今GPS观察到的右旋滑动速率为1—4 mm/a (杨少敏等,2008; Zubovichet al,2010).而博罗可努—阿齐克库都克断裂(以下简称博阿断裂)作为中天山与北天山地块的分界断裂,起始于巴尔喀什湖东侧,经阿拉湖南缘、阿拉山口、艾比湖到达吐鲁番盆地南缘,止于鄯善以南的觉罗塔格山山体内,是一条发育在古缝合带上的活动断裂,全长超过1 000 km,走向NW (邓起东等,2000),对于博阿断裂的滑动速率测算得出了不同的结果.杨晓平和沈军(2000)得到的精河县东南段晚更新世以来的走滑速率约为5 mm/a;Shen 等(2011)根据下切河流阶地的形成年龄与河流位错的关系,认为博阿断裂自西向东其走滑速率逐渐减小,阿拉山口附近晚更新世以来的走滑速率约为5 mm/a,吐鲁番盆地西缘晚更新世以来的走滑速率仅约为1.4 mm/a;Campbell 等(2013)得到晚第四纪以来在阿拉山口以北哈萨克斯坦境内的博阿断裂的右旋走滑速率约为2.2 mm/a;基于GPS数据估算的博阿断裂的右旋走滑速率为1.7—2.4 mm/a (李杰等,2010),对冲沟位错及洪积扇年代学的研究表明,博阿断裂精河段的活动速率为(3.1±0.3)mm/a (胡宗凯等,2019).尽管前人对该断裂已进行了较深入的研究,但主要集中在断裂展布和活动性方面,对古地震和最新活动时代均未深入分析;且该区域内所有NW向断裂规模宏大,但尚未发生过M≥7.5地震,最大发震能力值得商榷.前人关于形变带长度和单次位错研究认为,博阿断裂的阿拉山口—精河东南段为一整体,具有发生M7.5—7.8地震危险性(杨晓平,沈军,2000;Shenet al, 2011),但并未对断裂活动性分段予以分析,此外,1765年精河M61/2地震导致古城墙坍塌,房屋倒塌,但其发震构造尚待确定.
图1 青藏高原西构造结节和天山地区的构造及强震分布图F1: 塔拉斯—费尔干纳断裂; F2: 肯迪克塔什断裂; F3: 博罗可努-阿齐克库都克断裂;F4:温泉—博乐断裂;F5:四棵树—古尔图南断裂Fig. 1 Structure of western Syntaxis of Tibet Plateau and Tianshan area and distribution of earthquakesF1:Talas-Fergana fault; F2:Kindyktash fault; F3:Bolokenu-Aqikekuduk fault;F4:Wenquan-Bole fault;F5:Sikeshu-Guertunan fault
针对以上问题,本文在高分遥感解译的基础上,对断裂沿线进行地质考察,针对重点区段进行无人机航测,获取断裂位错量信息.并选择在地貌和沉积典型位置开挖探槽,对第四纪沉积层和堆积楔进行分析,之后结合新年代学测试结果,讨论断裂不同段落的古地震期次,厘定各段最新活动时间,并根据离逝时间探讨断裂各段地震危险性,为断裂各段地震危险性判定提供可靠地质依据.
天山以北博阿断裂主要沿NW−NNW方向展布,遥感解译其总长度约为370 km,其中巴尔喀什湖东缘—阿拉湖南总长约为250 km,以断裂走向和活动性质为依据将其划分为三段:阿拉湖以北段(长约100 km),活动地貌表现不明显,地貌宏观线性较好,迹线较平直,推测其为右旋走滑性质,未见近期活动变形带;阿拉湖段,长度约为100 km,以正断为主,变形带非常清晰,以陡坎地貌发育为主,未见明显的走滑活动性质,该段为阿拉湖盆地南缘控制断裂,严格控制了洪积扇与山体的分界,总体走向为285°,断裂沿山体边缘展布,局部段落呈弧形展布;阿拉山口—巴音阿门段,长度约为160 km,自阿拉湖东侧开始,断裂走向发生转变,由西侧走向315°进一步向NW偏转25°,总体走向为340°,以右旋走滑运动为主,兼具正断及逆冲分量.本文野外研究重点区域为中国境内的精河段,该段长约120 km (图2).以精河县城为界,由于其它断裂截切,将该段分为艾比湖段(阿拉山口—精河县)和精河以东段(精河县—巴音阿门)两段(图1,2).
图2 博阿断裂天山以北段的展布影像解译及活动性质(a)博阿断裂总体展布图;(b)艾比湖西岸断裂影像;(c)100棵树附近断裂影像;(d)精河县城断裂影像.图中黄色框为艾比湖段;蓝色框为精河以东段;红色剪头指示形变带.影像来自Google EarthFig. 2 The distribution and activity of the northern Tianshan mountain segment of Bolokenu-Aqikekuduk fault(a)The trace of Bolokenu−Aqikekuduk fault;(b)The trace of Aibi lake segment;(c)The trace of 100 trees segment;(d)The trace of Jinghe county. The yellow frame is Aibi lake segment,blue frame is east Jinghe segment,red arrows indicate deformation trace of the fault. The figures are from Google Earth
博阿断裂艾比湖段始于阿拉山口市西侧山麓戈壁滩,由于断裂北东盘上升,南侧山体所降雨水未能继续北流进入低洼河谷,而是沿断裂自艾比湖南侧注入该断裂控制的大冲沟,断裂在湖相阶地上呈线性展布,最新活动迹象最为明显,也是本文研究的重点区段(图2a,b).随后断裂延入艾比湖,湖内长度约为30 km,并在湖东侧的河湖三角冲积漫滩出露,延伸至精河县城西北,由于耕地改造等原因,行迹不清.而精河县城北部的达流村可见断裂控制耕地边界,呈断续线性展布,断裂南侧为耕地,而北侧则土壤沙砾化严重,几无植被发育(图2d),断裂在此处对植被控制的现象与艾比湖西侧几乎一致,因此推测其活动性质具有同样的正断性质.艾比湖段总长度约为60 km.
在艾比湖西岸湖相阶地上发育线性断裂地貌,其严格的控制了该区域的植被发育,由于断裂北侧抬升,山前冲沟的水补充减少,所以其植被较少;而断裂南侧由于地下水相对丰沛,因此,沙漠植物梭梭草生长较为密集,也相对高大.该段形变带长约2 km,走向为340°(图2b,3).艾比湖由于晚更新世以来水位剧烈退缩,发育有多级古湖岸线,其原始形态为一天然坝体,宽数米,由于湖水波浪的冲击作用,坝体顶面一般高出其对应湖面2—3 m,坝体后缘发育小型水体,也叫障壁塘.博阿断裂(红色剪头之间)在Ⅱ级湖岸线或阶地面以上线性明显,表现为正断陡坎,陡坎高约1 m,Ⅱ级湖岸线高出湖面约18 m,宽约200 m,前人研究其形成年代距今4.5—6 ka (吴敬禄,1995),表明其形成于全新世,但由于断裂南侧阶地被冲沟所破坏,位错标志不能获得,因此只能判定该断裂为全新世活动断裂,尚未获得其右旋位错距离(图3).
图3 艾比湖西北岸断裂无人机航测合成影像(a)和高程图(b)(图中红色箭头指向断层滑动面)Fig. 3 Unmanned aerial vehicle aerial composite image (a)and digital elevation model in the northwest bank of Aibi lake (b)where red arrows indicate fault scarp trace
在Ⅱ与Ⅲ级湖岸线之间,垂直线性陡坎处开挖探槽,探槽长为20 m,最大深度为3.5 m(图3,4),当挖至陡坎南侧时,冲沟一侧便有水涌出,探槽南端开始坍塌,无法继续开挖,断裂发育在陡坎中部位置,由于张性滑动和充填作用,断裂面两侧多充填细沙,其中砾石均发生沿滑动面的定向排列,形成厚约为20 cm的断层破碎带,局部发育多个次级滑动面(图4c).由于细砂的渗透系数为0.001—0.005 cm/s,而砂砾石为0.01—0.1 cm/s,渗透系数相差20倍以上(毛昶熙,2009),因此,细砂层具有明显的阻隔作用,而地下潜水多沿渗透系数小的冲沟一侧流走,水量极少渗透到断裂北侧,加之断裂北侧抬升,与南侧存在明显高差,导致其潜水层埋深过大,地上植被很难生存,从而致使上述断裂位置和植物线性重合的地貌现象产生.
图4 探槽东壁局部照片图中红色箭头指向断层迹线,E1和E2为地震崩积楔体Fig. 4 Splicing photos of east wall of the trenchThe red arrows indicate fault slip face,E1 and E2 are earthquake colluvial wedge
探槽东西壁均揭露多套地层,但断裂两盘地层可对比性差,其中两盘最顶部均有棕黑色砾石土层发育.上盘厚约60 cm,探槽西壁该土层底部存在砾石定向排列现象,并发育纵向裂缝;下盘厚度小于30 cm,层理不发育,根系发育,底部弱钙质胶结;而探槽东壁该砾石土层无明显错动现象,发育崩积楔E1,断裂错动面上达地表附近.仅下盘发育的土层为: 黑色砾石粗砂层,厚度约为60 cm,水平层理发育; 黑褐色砂砾石层,水平层理发育,偶夹薄砂透镜体;黄褐色砾石粗砂层,水平层理发育,可见厚度大于2 m,未见底.仅上盘发育的土层有:灰黄色细砾粗砂层,层理发育,钙质胶结,断裂一侧明显错动,裂缝发育;褐色钙质胶结砾石夹砂层,近断裂一侧发生明显错动,与上层共同形成三角形楔体E2,该崩积楔由砂砾石组成,砾石定向和裂缝发育,推测在楔体形成之后,后期地震活动,在断层面错动摩擦牵引作用下,楔体内发生不均匀轻微错动形成(图5,6).
探槽西壁可见断裂发育于陡坎中部,产状为220°∠68°,为正断性质,北盘上升,南盘下降.断裂带内存在砾石定向,充填的砂脉中发育平行断裂的次级裂缝(图5c).断层下盘为洪积物,分选较差,层理较清晰,均为水平的粗砂砾石层,未见变动.上盘为砾石和砂层混杂堆积,层理有掀斜,并有崩积楔发育,在探槽西壁崩积楔附近采集6个光释光(optically stimulated luminescence,缩写为OSL)年代样品,测试结果其年代均介于1.95—3.04 ka之间(表1,图5b,6),存在年代倒序现象,可能是由于样品要求为粉细砂,所采样品可能是沿裂隙贯入的沉积物,或者还可能存在个别样品年龄测量误差,但从大多数样品年代可基本推测楔体形成在1.95—3.04 ka之间,也即代表在此之间存在一次地震事件,以楔体周边样品年代推测其年代更可能为2 ka.而楔体内发育砾石定向和纵向裂缝,且探槽东壁砾石层①下部还发育崩积楔E1,且错动面接近地表,表明断裂在2 ka之后还发生另一次事件,根据该层沉积物异常松散,其可能对应该区域的历史地震事件,下文将详细分析.而E2楔体底部距下盘地表覆盖层①底的距离为0.9—1.1 m,此距离应为两次地震事件的累积位移(图4),因此综合判定单次正断量在0.5 m左右.
表1 艾比湖探槽光释光年代测试结果Table 1 Optically stimulated luminescence dating results of Aibi lake trench
图5 探槽西壁拼合照片及局部断裂细节图(a)探槽西壁拼合;(b)探槽西壁采样分布;(c)断层带细节.红色剪头指示断层面Fig. 5 Splicing photos of west wall of the trench and details of faults section(a)Splicing photos of west wall of the trench;(b)The aging sample sites;(c)The details of the fault section. Red arrows indicate fault slip face
图6 探槽西壁剖面素描图及样品位置Fig. 6 Profile of the western wall of the trench and the location of samples
断裂在精河县城东北山麓冲/洪积扇连续出露(图2c),经100棵树到天山山麓的巴音阿门北西,断裂穿过洪积扇,呈反向坎,其错动了一系列洪积扇上发育的冲沟,由于冲沟发育较混乱,可参照的标志较为混乱,因此根据无人机航测结果,仅选择最可靠的一条冲沟为错断标志,其位错距离约为40 m (图7),由于其切割深度较大,发育年代应在洪积扇形成之后或同期.距离该处5 km的同期洪积扇探槽中接近地表的砂砾石层样品作为该洪积扇的形成年代,测年结果为(10.67±0.38)ka,而洪积扇形成之后风成黄土沙层最底部年代为(7.16±0.53)ka,因此洪积扇形成年代介于7.16—10.67 ka之间(图8,表2),判定其右旋活动速率应大于等于4 mm/a,此结果与近期根据天山洪积扇年代对比所得速率相当(胡宗凯等,2019).
图7 100棵树附近无人机航测DEM及断裂冲沟位错(白色箭头之间为断裂)Fig. 7 DEM and dislocation gully by fault near 100 trees (The section between the white arrow is fault trace)
距离100棵树西约5 km处(位置见图2c),同期洪积扇垂直陡坎开挖大型探槽,深度约为8.5 m,探槽壁大致可分为两套沉积物,其顶部为风成黄土沙层,断裂上盘薄,厚度仅为30 cm,而断裂下盘厚度可大90 cm,最上层层理较连续,未发现错动证据,采集OSL样品,其结果为(3.70±0.28)ka,可判断此后未发生地震错动事件.下部为总厚度达6 m的洪积扇砾石粗沙互层沉积,偶夹细沙透镜体,层理发育清晰,根据年测定结果,其年龄介于(10.26±0.32)—(18.27±1.2)ka之间,表明此洪积扇应形成于末次盛冰期之后,断裂带表现为地层轻微逆冲错动,各层均发生明显的“Z”字型膝折,由于该断裂以右旋走滑为主,其垂向位错量不大,不易进行细致的古地震分析,但从剖面中膝折位错大小至少可以分辨出4期事件.其中,Ⅰ期为风成沙土层,中下层样品年代分别为(4.86±0.35)ka和(7.16±0.53)ka,JH-TC1-02样品层位有明显错动,有下部的砾石层混入,其中两处有明显定向,且错动处黄土颜色呈暗褐色,而两侧为浅黄色;而该错动现象距离主断裂约1 m,推测可能是由于断裂压扭作用在断层下盘形成的次级断裂活动所致,由于位错较小,其下方的砾石混杂沉积层中不易观察,而主断裂面由于地表水流侵蚀作用而未能在黄土层中保持错动面(图8b),所以认为该段在3.70—4.86 ka之间发生了最新错动,根据黄土层底层和砾石顶层位错距离,其位错量约为0.5 m.Ⅱ期在距离地表3—4.5 m处,洪积层位错约为1.5 m,对应的年代为(10.26±0.32)ka.Ⅲ期在距离地表4.5—6.5 m处,膝折位错近2 m,对应的年代为(18.27±1.2)ka.Ⅳ期在距离地表6.5—8.5 m处,膝折只揭露一半,可见位错近2 m,真实位错应大于该距离(图8,表2).显然,这肯定不是连续的四次地震所致,其中肯定有多次事件无法分辨,但基本表明其晚更新世以来断裂一直活动.由顶部风成沙土层未错动,表明其离逝时间至少有3.7 ka.而根据Ⅰ期与Ⅱ期的时间间隔约为3 ka,如果以此为断裂强震复发周期且为准周期活动,判断该断裂离逝时间已经超过复发周期,其发震危险性不容忽视.近期其它研究组在巴音阿门附近开挖了两个探槽,结果表明精河段古地震发震周期为1—3 ka,但最近一次古地震事件的离逝时间为6.2—6.9 ka (胡宗凯等,2020).该结果与本文相近,但离逝时间相较本研究差别较大,超出其判定的发震周期两倍多,推测可能是由于全新世以来沉积均为黄土沉积,沉积速度慢厚度小,且容易愈合或被侵蚀,导致最近事件被遗漏.
图8 (a)探槽2照片拼接图及OSL年代测试结果;(b)黄土层中错动细节Fig. 8 (a)Splicing photos of trench 2 and dating results of OSL;(b)The location detail in loess
表2 探槽2样品测试结果Table 2 Test results of trench 2 sample
由以上两个探槽和地貌分析,博阿断裂自阿拉山口市—天山山麓带大致可分为艾比湖段和精河以东段两个活动段.
艾比湖段西端分界根据前人在阿拉山口西南对该断裂地表形变带进行了测量,据断层陡坎的坡向,可将该断层陡坎带分为南段和北段,其中南段断层陡坎坡向南西,北段坡向北东(杨晓平,沈军,2000).表明在阿拉山口西南断裂可能发生转变,而南段与艾比湖西岸陡坎倾向一致,转变位置距离艾比湖探槽约8 km.该段断裂具有张扭性质,倾向南西,表现为明显的正断性质,2 ka以来发生过两次强震事件.
精河以东段,长度约为40 km,精河北—100棵树段具轻微压扭性质,倾向南西,100棵树—天山山麓带呈为反向坎,南盘上升,北盘下降(沈军,杨晓平,1998;杨晓平,沈军,2000),其最新活动时代约3.7 ka以前.地震构造图显示在分段位置附近,正好有两条断裂与其相交,分别为西南侧的温泉—博乐断裂和东北侧的准噶尔盆地南缘断裂西段四棵树—古尔图南断裂(图1),由于这两条断裂研究较少,尤其是前者,到精河附近处于隐伏状态.而后者控制了天山和准噶尔盆地的分界断裂,断裂在距离精河50 km的古尔图南,切割基岩山体达8 km,地貌显示明显,而到精河附近,断裂南侧为陡立的基岩山体,北侧为山麓冲洪积及黄土沉积,但影像未解译出明显活动现象,其延伸恰好至精河县城东侧,而博阿断裂在该处两侧地貌表现也发生明显变化.因此,推测四棵树—古尔图断裂可能对断裂的活动性分段起控制作用,当然也不能排除温泉—博乐断裂的控制作用,还需要进一步研究予以确定.
1765年历史地震记载精河发生M61/2地震,历史记载当时灾情如下:在呈报皇帝的奏折和军队奏报均有记载.其中满文奏折记载“正月二十日自卯时,地震不宁,至未时方休”.伊犁将军明瑞奏报:“本年正月二十日卯时,伊犁地方时震时止约二刻时间,闻听泾水(精河)地方地震严重,城垣房屋倒塌,居民受伤.” (新疆维吾尔自治区地震局,1981).而在清代官修的编年体史料长编《清实录》(清高宗实录,卷七二九,页八)不仅有更详细的破坏记录,并有赈灾措施:“据伍弥泰奏称,晶河(今精河)地方.正月二十日地震.城垣、兵房、仓厫及商民房屋俱倒坏.兵丁王哲昭压毙.城垣兵房仓厫,陆续修补.商民房屋,听其量力自修等语,晶河地震,压毙兵丁,着加恩照兵丁失足落水例,赏恤,其商民倒坏房屋,令其量力自修,未免拮据,着伍弥泰查酌,如系征租房间,量免数月半年,以纾商民之力,若官地商民自盖,征房身地租者,宽免一年.而作为处理军队事务最高决策机构的军机处,详细记录了破坏的具体数字.“军机处满文月褶档,乾隆三十年二月二十日,据晶河管庄都司张文等呈称,正月二十日自卯时地震不宁,至未时方休,我至晶河查看城墙、门楼、仓厫、营房,其城垣倒塌七丈五尺,门楼、仓厫、营房倒塌69间,仓厫、墙垣破裂二十四间,商人房屋墙垣倒裂五十四间.除兵丁王哲昭伤势过重殒命外,官兵商人皆安(谢毓寿,蔡美彪,1987).
有关安阜城的有以下记录:清于乾隆二十四年(1759年)在精河修筑安阜城.古城位于精河县城西北部约一千米处,史料所载内容十分简单: “乾隆二十三年(1758年),定边将军兆惠率兵在精河境内全歼准葛尔叛军阿睦尔撒纳残部,活捉首领布尔古特吉珲.第二年,清政府在精河东岸筑建一城,城南北长166米,东西宽122米,呈矩形,取名安阜,设典吏一名,隶属乌鲁木齐都统管辖.当时此城驻扎几百丁士,再安置百十官绅、商贾,其城也足可居矣!可随着人口和商贸不断发展,城区太小的鄙陋就显现出来,乾隆四十八年(1783年),清政府又于城东二公里处再筑新城,仍称安阜城”(赵世勇,2015).
新城位于精河县城内,由于后期建设被破坏,几无遗迹保存.旧城现作为遗址仍保存较好,城墙形态基本保存完好,Google Earth影像显示,其大致呈边长140 m的正方形(图2d).四周城墙大半犹存,黄土夯筑,高约4 m.南墙中段辟城门,宽约4.5 m,高4 m.城墙底部宽度约5 m,顶部宽约3 m,至今依然坚固,只有表面有轻微风化.
由于1765年地震时,新城还未兴建,因此以上破坏记录均为旧安阜城,由于城墙作为军事堡垒,每隔一定距离有加宽矩形墩台(即马面),非常坚固,破坏不很严重,因而只有少部倒塌,而城内房屋、仓库均被破坏倒塌,其烈度可能大于历史地震所记Ⅷ度,可能接近Ⅸ度,其震级有极大可能超过历史记录震级为M61/2,原因为:① 其发震构造为艾比湖段,其震中位置可能更偏向湖盆中心或西侧,以艾比湖西岸至精河北为地震破裂带,总长度约为60 km,依据适用的矩震级−地表破裂长度经验关系式M=5.16+1.12lgL(Wellset al,1994)其震级为M7.15;② 由于1765年,当地除安阜城之外很少居民点分布,因此当地震害记录很不完备,安阜城位于活动段尾端,且古城距离断裂距离约4 km,加之城堡坚固,破坏性明显减少,但其烈度可达Ⅸ度,导致震级低估.
对博阿断裂中国境内的艾比湖—巴音阿门段进行高分遥感解译及沿线调查,并在重点区段进行无人机航测和探槽古地震研究,结果表明该断裂为全新世活动断裂,并根据活动性质差异和断裂交切关系,将断裂分为艾比湖段及精河段,其中艾比湖段的探槽开挖及测年结果显示,距今约2000年发生了两次错动事件,最新一次错动至地表,其地表破裂带可延至精河县城以北,推断极有可能为1765年精河M61/2地震的发震构造,表现为张扭活动性质,因其破裂长度达60 km,历史记载该地震破坏性强,其震级有被低估的可能,震级可能达到7级.而精河以东段已超过其复发周期,其地震危险性需要关注.本研究将古地震研究与历史地震灾害记录考证相结合,首次提出该断裂为历史地震发震构造.但由于工作程度不足,其中部分结论还有待进一步证实,如古地震研究还基于一个探槽研究结果,尚需要更多的材料予以证实,如在艾比湖东南岸或精河县城北侧开挖新探槽,为断裂最新活动活动时代和古地震期次提供更强的约束.另外,断裂活动性分段除断裂截切原因之外,造成断裂正断活动分量的动力机制,以及断裂上下盘在百公里之内发生转变构造成因也需要更深入的研究.