ENGINEERING EFFECT OF THE HIMALAYAN OROGEN AND ENGINEERING GEOLOGICAL ZONING OF CHINA-NEPAL RAILWAY
-
摘要: 造山带内动力作用的工程效应是地学与工程领域基础科学研究的前沿问题。结合GIS技术和相关地学理论,对喜马拉雅造山带由构造划分的高喜马拉雅、低喜马拉雅和特提斯喜马拉雅3个地块进行研究,揭示不同地块工程地质特性和灾害效应受内动力制约的普适性规律。提出:在挤压碰撞造山机制作用下快速隆升的高、低喜马拉雅,地震以逆冲型为主,地震强度大、频率高,水平应力相对较大,主应力方向近NE-WS方向,地形向大高差发展、河流下切强烈,山地灾害严重;而属于拆离地系,处于相对沉陷状态的特提斯喜马拉雅,地震以正断型为主,地震活动性相对较弱,水平应力相对较小,主应力方向近E-W方向,地形演化向着减弱地势的趋向发展,雪崩灾害严重;此外,高喜马拉雅特有的海洋性冰川地貌、冰湖和冰川泥石流,可能是控制跨喜马拉雅山铁路线路方案的重要问题。基于上述各地块工程效应存在显著差异的认识,提出以构造划分作为铁路工程地质分区的建议,并以拟建中尼铁路交通廊道为例,绘制了工程地质分区图。研究有助于将造山带理论推进到工程应用层面,为铁路大范围方案比选阶段,广域、高效、低成本地获取信息提供了新途径。Abstract: The engineering effect of endogenic processes in Himalayan orogen is a frontier problem of basic scientific research in the field of geosciences and engineering. In the longitude direction, based on the tectonic and lithological differences, the Himalayan orogenic belt is divided into 3 parts: Tethys Himalaya, Higher Himalaya and Lesser Himalaya. These studies reveal the universal law that the engineering geological characteristics and disaster effects of the different parts are restricted by endogenic processes. This paper presents the following arguments. Firstly, the rapid uplift of Higher and Lesser Himalaya under the mechanism of compressional collision and orogeny, the earthquakes are mainly thrust fault type and have high earthquake intensity and high frequency. The dominant direction of geostress is close to NE-WS. The terrain develops to a large height difference. The river cuts strongly. The Mountain disaster is serious. Secondly, the Tethys Himalaya belongs to the detachment system and is in a relatively subsidence state. The earthquakes are mainly normal fault type. The seismicity is relatively weak. The geostress direction is close to the E-S. The terrain evolution is trending towards weakening the terrain. And the avalanche disaster is serious. Lastly, the marine glacial landforms, glacial lakes and glacial debris flows are unique to the Higher Himalaya. These may be important issues in controlling the trans-Himalayan railway line scheme. According to the understanding that there are significant differences in the engineering effects of the above parts, this paper puts forward the suggestion of taking the structural division as the railway engineering geological zoning. Taking the traffic corridor of the proposed China-Nepal railway as an example, the engineering geological zoning map is drawn. The research is helpful to advance the theory of Himalayan orogen to the level of engineering application. At the stage of railway large-scale scheme comparison and selection, it provides a new way to obtain information in wide area, high efficiency and low cost.
-
0. 引言
自始新世以来,印度板块向北俯冲,沿雅鲁藏布江缝合带向亚洲大陆南缘俯冲挤压,褶皱隆升和断块挤压上升共同作用使喜马拉雅山快速隆起,成为世界上海拔最高的山脉(Gansser,1964; Li et al., 1979; Harrion et al., 1992; Yin,2006)。它不仅是东亚大陆与南亚次大陆的天然界山,也是目前中国与印度、尼泊尔、巴基斯坦和不丹诸国国界的基本轮廓。“一带一路”倡议的优先领域是铁路的互联互通,在我国2016年7月颁布的《中长期铁路网规划》中,面向“一带一路”国际通道的喀什—红其拉甫(中巴铁路国内段)、日喀则—吉隆(中尼铁路国内段)、日喀则—亚东等口岸铁路均列入其中。在喜马拉雅造山带修建铁路将是人类铁路建设史上新的挑战,国、内外均无经验可借鉴,前期研究工作亟需开展(姚令侃等,2012)。
在干线工程的可行性研究阶段,为了查明铁路所经地区的工程地质条件,一般是通过收集与研究已有资料,如地形图、区域地质调查报告、地震烈度区划图等,概略了解拟建交通走廊的地形、地貌、地质、水文、气象、地震等自然条件,进行沿线工程地质分区。但在复杂艰险、构造活动强烈山区,地貌可能在铁路百年服务期内就发生显著变化,加之大量深埋长隧工程的修建,使得铁路越来越多的希望在建设前期就获得区域动力地质现象、地貌演变规律、地应力状态等信息。而目前常规资料难以满足这种需求,甚至,对跨境铁路而言,有些资料境外路段根本就没有。
造山带作为独特的构造单元,一直是国内外研究地壳结构和地球动力学的热点地区,长期以来许多学者和专家在此完成了多方面的课题研究,取得了大量的研究成果。然而,造山带地区的工程地质研究远远滞后于造山带理论的进展。随着我国大型工程建设不断向地质条件复杂地区延伸,不少学者如王思敬(2013)、黄润秋等(2017)、许强(2012)纷纷呼吁应将新的地学基础理论引入到工程地质领域,更进一步,部分学者已开展了一些先导性工作。彭建兵(2001a, 2001b)在其提出的区域稳定动力学中,论述了3种类型造山带(逆冲推覆型、伸展型和走滑型)各构造部位的区域稳定性特征,建议作为工程选址的原则;张永双等(2004)以秦岭造山带为例,提出造山带形成演化过程对工程地质问题的制约作用,介绍了造山带常见的地质灾害类型和工程地质问题。但总体上看,现有文献仍多以提出命题或倡导为主,将造山带理论深入到工程实践的研究甚少。
喜马拉雅山是陆陆碰撞的典型地区,造山带地质演化和壳幔深部动力学过程与表面地壳运动之间的相互作用,无疑会在地表工程圈得到体现。若能利用造山带理论直接获取工程地质特性,不仅可望为铁路大范围方案比选阶段,高效低成本地获取信息提供新途径,还有可能获得常规方法难以获取的信息。本文即从造山机制出发,结合GIS技术和相关地学理论,对喜马拉雅山由构造划分的3个地块进行对比分析,揭示不同地块工程地质特性和灾害效应受内动力制约的普适性规律,并以中尼铁路为例,提出一种以构造线作为分区界线的工程地质分区方法。
1. 喜马拉雅造山带主要地块及内动力机制
从广义上讲,喜马拉雅造山带由喜马拉雅弧和东、西构造结-挤压转换带组成(许志琴等,2007),其中喜马拉雅弧是造山带的主体,长达2000多公里,是碰撞造山作用最强烈的区域(图 1),为本文的主要研究区。
虽然构造地质学家对印度与欧亚大陆碰撞方式和过程有诸多争议,但普遍达成共识是,印度板块与欧亚板块发生碰撞后,在印度河—雅鲁藏布缝合带(ITSZ)以南,由南向北依次形成了一系列的逆冲型断层,主要有主前锋断裂(MFT)、主边界断裂(MBT)、主中央断裂(MCT);此外,还有一条以拆离张拉作用为主的藏南拆离断裂(STDS)(Burchfiel et al., 1985; Tapponnier et al., 1986; Yin,2006)(图 2)。
这些构造线将整个喜马拉雅地体分割成为特提斯喜马拉雅、高喜马拉雅、低喜马拉雅和次喜马拉雅4个部分(Brookfield,1993; Lefort,1996; Ding et al., 2017)(图 3)(构造线的确定主要依据我国青藏高原1 ︰ 25万的区域地质报告,并参考了尼泊尔1 ︰ 25万地质图)。其中次喜马拉雅基本属于丘陵和平原地貌区,不作为本文研究对象。
特提斯喜马拉雅位于ITSZ和STDS之间。STDS是世界上最大的正断层体系(Burg et al., 1984; Burchfield et al., 1992),在拆离滑脱作用下,高喜马拉雅(下盘)和特提斯喜马拉雅(上盘)的最大差异性隆升达3 mm ·a-1以上,因此,特提斯喜马拉雅处于相对沉陷状态。岩片主要形成于地壳中、上部的浅层冷环境,多是特提斯洋消亡后的海相沉积层,并在STDS附近存在低级变质作用,主要岩性为砂岩、砂泥岩、白云岩、灰岩和页岩。
高喜马拉雅位于MCT和STDS之间。由于沿着MCT的逆冲和沿着STDS的拆离同时发生,促使高喜马拉雅快速挤出,形成以高级变质岩、普遍的深熔反应以及高温韧性变形为主要特征的“热碰撞造山”特点(Godin et al., 2006)。现今高喜马拉雅的隆升速度约为4~6 mm ·a-1(Jackson et al., 1994)。岩片来自于地壳下部,受高温高压变质作用影响,岩石以片麻岩、花岗岩和片岩为主。
低喜马拉雅位于MBT和MCT之间,系印度板块俯冲至欧亚板块之下,印度板块表层物质被欧亚板块刮擦剥离而成,并伴随着低级变质作用,表现出“冷碰撞造山”特点(Law et al., 2006)。现今局部隆升速度约为2~3 mm ·a-1(Jackson et al., 1994),岩片来自于地壳中、上部,以页岩、泥岩、石英砂岩及其变质岩等堆积而成的混杂岩为主(许志琴等,2013)。
从宏观角度,人们习惯于将喜马拉雅造山带统一按挤压隆升区域处理,但若按构造划分的3个地块细化分析,各地块的内动力机制截然不同,并直接表现出各地块地层岩性的差异。我们认为,内动力机制不同的构造单元,工程地质特性也应存在差异,以下将对比分析各地块的其他工程地质特性和灾害效应。
2. 各地块工程效应分析
2.1 地震分析
地震是岩石圈表层的构造动力学效应主要表现之一。喜马拉雅造山带是世界上地震活动最强烈、大地震频繁发生的地震带之一。3个地块内动力机制的差异应能直接在地震上得到体现。
2.1.1 地震类型
根据中国地震局地壳应力研究所提供的中、强地震震源机制解,将喜马拉雅造山带1947年7月~2007年5月期间,震源深度0~100 km的80个岩石圈地震按地块分别统计地震类型(图 4,表 1)。
表 1 各地块地震类型占比Table 1. Proportion of earthquake types in different blocks地块 正断型/% 正断-走滑型/% 走滑型/% 逆冲型/% 逆冲-走滑型/% 特提斯喜马拉雅 45.5 9.1 40.9 0 0 高喜马拉雅 0 0 34.4 59.4 0 低喜马拉雅 4.0 0 28.0 64.0 4.0 由表 1可知,特提斯喜马拉雅所发生的地震大部分为正断型,其次为走滑型,未见逆冲型;而高、低喜马拉雅绝大部分为逆冲型,其次为走滑型。
2.1.2 地震活动性
按照我国现行概率地震危险性评价规范,b值和v4是反映地震活动性的主要参数(高孟潭,2015)。b值是统计区地震震级-频度关系(G-R模型)中的斜率,反映了大、小地震发生的比例关系。b值越大,表示大震占的比例越小,反之,b值越小,大震占的比例越大(Gutenberg et al., 1956)。v4是地震年平均发生率,其值越大,代表这个地区4.0级以上地震越频繁,地震活动性水平越高。
除了以上两个参数,我们还定义了地震能量密度(统计区内历史地震的总能量除以面积)作为表征地震强度的指标;此外,震源深度是体现地震对工程影响程度的重要指标,震源深度越浅,对工程影响越不利。根据USGS(美国地质调查局)公布的历史地震目录,对喜马拉雅造山带1958~2018年期间,震源深度为0~100 km的1545个3.0级以上的岩石圈地震,分地块统计上述参数列于表 2~表 4。
表 2 各地块地震活动性参数Table 2. Seismic activity parameters of different blocks地块 地震数目 b值 v4 地震活动性排序 特提斯喜马拉雅 313 1.23 4.15 3 高喜马拉雅 651 1.05 9.18 1 低喜马拉雅 581 1.12 8.79 2 表 3 各地块地震能量密度Table 3. Seismic energy density of different blocks地块 面积/km2 总能量/面积/J·km-2 地震能量密度排序 特提斯喜马拉雅 244883.3 60841.67 3 高喜马拉雅 185251.6 430864.20 1 低喜马拉雅 96102.3 423588.49 2 表 4 各地块4级及以上震源深度Table 4. Focal depths of earthquakes with 4 magnitude and above in different blocks地块 平均震源深度/km < 20 km地震比例/% >20 km地震比例/% 震源深度排序 特提斯喜马拉雅 30.58 26.09 73.91 3 高喜马拉雅 28.19 35.71 64.29 2 低喜马拉雅 21.92 50.75 49.25 1 由表 2~表 4可看出,从地震活动性和地震能量密度来看,高喜马拉雅最高,其次是低喜马拉雅,最低是特提斯喜马拉雅;从震源深度来看,低喜马拉雅最浅,其次是高喜马拉雅,最深是特提斯喜马拉雅。
造山机制决定了地震类型,地震类型又决定了地震强度,以上分析表明,3个地块内动力机制的差异确实在地震指标上得到体现。此外,根据造山带理论,在地壳隆起区一般为高水平应力区,而现代沉降区则为低应力水平区(陈彭年等,1987)。高喜马拉雅隆升速度最快,低喜马拉雅其次,而特提斯喜马拉雅处于相对沉陷,由此可推测高喜马拉雅相对水平地应力最大,低喜马拉雅其次,特提斯喜马拉雅水平地应力最小。更进一步,b值越大,说明该地区岩体非均匀性和破裂程度越高(Mogi,1962;杨百存等,2017)。按照表 2中3个地块b值大小的排序,可推测高喜马拉雅岩体均匀性最好,破裂程度较轻,其次是低喜马拉雅,而特提斯喜马拉雅岩体的非均匀性和破裂程度都最高。该推论在有关3个地块岩性与成因的基础上,又增加了岩体相对破裂程度的信息。
2.2 构造应力方向分析
在造山带区域,地应力方向也是选线与隧道工程设计需要获取的重要参数,一般通过物理钻探手段获取,在构造地貌显著地区,用地质证据分析地应力的方法具有优势(许忠淮,1990)。喜马拉雅造山带是世界上内动力塑造地貌效应最强烈的地区,若能由地貌反推内动力特征,无疑是获取地应力信息的一种低成本的手段。
20世纪70年代末,瑞士学者Scheidegger提出了地表发育形态的对抗性原理,成为地球表面动力学的基本理论之一(Scheidegger,1979)。根据对抗性原理,在新构造应力场的作用下,岩体内部会形成一对共轭的剪裂面(X型破裂,抗侵蚀能力弱)。此剪裂面在外营力的风化或侵蚀作用下,易破坏形成断裂带,而河流多沿构造节理及破裂面发育。因此沿节理破裂发育的水系优势方向可以反映新构造运动以来的应力分布状态。根据莫尔强度理论,新构造应力场的主压应力方向即为共轭剪节理的锐角角平分线。若将区域内的河流看作一组共轭剪节理,则节理的锐角角平分线则为构造应力场的主应力方向(Scheidegger,2004)。
首先提取研究区河网,用直线连接河流的自由端或交汇点,湖泊则取其长轴方位,得到3个地块水系折线图(图 5)。然后取5°为一个单元,统计河段总长度,并绘制各地块水系方向玫瑰图(图 6)。再采用密集度方法(余庆余等,1983),得到第一优势走向(图 6中红线),在第一优势方向两侧30°以外的区域寻找第二优势走向(图 6中蓝线)。最后以两个优势方向相交的锐角等分线(图 6红色箭头)作为主应力方向。
由图 6可知,高喜马拉雅第一优势方向为82.5°,第二优势方向为47.5°,主应力方向为65°;同理,低喜马拉雅主应力方向为72.5°;特提斯喜马拉雅主应力方向为115°。从宏观角度,受印度板块向北东向挤压欧亚板块的影响,喜马拉雅造山带的广大地区,均在一个广阔的北东—南西向挤压应力场的作用范围内,根据上述分析,按照水系推求应力场方向,高、低喜马拉雅主应力方向基本与这一结论相符;但在特提斯喜马拉雅,主应力方向为东偏南25°,与宏观结论不一致,说明这一地块的岩石圈应力场受近东西向的水平扩张应力控制。
2.3 地貌特征
地壳挤压作用越强烈,山地隆升的速度越快,形成的褶皱山、断块山发育越快,向着增强地势的趋向发展;同时河流下切作用越强烈,切割深度越大,河流纵断面越陡峻。反之,当山地上升减弱或趋于稳定,风化、剥蚀、搬运与沉积作用等相对占优势,山地就会剥蚀降低,向着减弱地势的趋向发展。高喜马拉雅隆升最快,低喜马拉雅其次,特提斯喜马拉雅相对沉陷,这种隆升状态的差异也应在地貌上得到反映。下面按3个地块对地貌进行统计,列于表 5。
表 5 各地块地貌指标Table 5. Landform indices of different blocks地块 平均海拔/m 地形起伏度/m 地表切割深度/m 坡度/(°) 特提斯喜马拉雅 4785 647.6 292.9 16.5 高喜马拉雅 3394 1060.8 510.5 25.7 低喜马拉雅 1675 873.9 416.0 24.2 从表 5可看出,地形起伏度、地表切割深度和坡度这3个表征地形陡峻程度指标的排序与3个地块隆升速度的排序完全一致。
河流纵坡也是对构造活动响应最敏感的地貌参数。分布在喜马拉雅山脉南北向的主要河流有16条,其中源于特提斯喜马拉雅,而后流经高喜马拉雅和低喜马拉雅3个地块的河流有9条(图 7),为了方便起见,将9条河流由西至东,分别以R1~R9命名,分别是阿勒格嫩达河(R1)、Mahakali河(R2)、Bheri河(R3)、Kali-Gangdaki(R4)、Marsyangdi河(R5)、Budhi-Gandaki河(R6)、吉隆藏布河(R7)、波曲-Sun koshi河(R8)和朋曲—阿伦河(R9)。提取出这9条主要河流的纵剖面,并计算每条河流位于不同地块河段的河流纵坡列于表 6。
表 6 河流纵坡Table 6. Longitudinal slopes of rivers河流名称 河段地块 特提斯喜马拉雅 高喜马拉雅 低喜马拉雅 R1 0.0140 0.0233 0.0038 R2 0.0552 0.0632 0.0073 R3 0.0284 0.0319 0.0090 R4 0.0228 0.0088 0.0089 R5 0.0500 0.0350 0.0066 R6 0.0475 0.0344 0.0087 R7 0.0206 0.0351 0.0152 R8 0.0199 0.0551 0.0054 R9 0.0021 0.0228 0.0131 平均值 0.0289 0.0344 0.0087 由表 6可看出,高喜马拉雅河流纵坡最陡,其次是特提斯喜马拉雅,最后是低喜马拉雅。一般情况下,河流纵坡应为上游陡下游缓,但在喜马拉雅造山带,处于中游的高喜马拉雅河段河流纵坡最陡,我们认为3个地块的差异性隆升是导致这种反常现象的主要原因。
造山作用过程是漫长的、缓变的、相对稳态的地质构造动力学演化过程,时间单位以百万年、千万年乃至亿年计。各地块地貌特征差异的现象是3个地块长期演变的积累效应,可预见的是,起码在百年至千年时间尺度内,高、低喜马拉雅地势将继续向着增大的方向发展,而特提斯喜马拉雅地势向着减小的方向发展,可以作为铁路服务期内地貌和孕灾环境演变预测的基础。例如,高、低喜马拉雅的河流会继续下切变陡,而特提斯喜马拉雅的河流会持续上涨变缓,可为铁路选择线路高程和防洪设计提供依据。
2.4 气候与冰川分布特征
造山带的形成演化过程与环境灾害密切相关。由于青藏高原隆升,喜马拉雅山脉异峰突起,在冬季阻挡来自北方的大陆冷空气流入印度,同时迫使西南季风在穿越山脉向北移动之前抛弃大部分水分,从而造成印度一侧的巨大降雨量和西藏的干旱状况,根据印度洋暖湿气流可以到达的前锋线(降雨量500~1000 mm ·a-1,海拔3500~4000 m),可大致划分为干旱寒冷气候和温带-亚热带季风气候两个气候区;在喜马拉雅山南坡,海拔高度急剧下降,气候则呈从寒温带、温带至亚热带的垂直变化,这就是喜马拉雅造山带气候的基本格局(林振耀等,1981)。
冰川按照物理性质,可以分为大陆型冰川和海洋性冰川两大类。大陆型冰川是在前锋线以北的大陆气候条件下形成的,活动性小,积累与消融均比较缓慢,冰湖和冰川泥石流也不发育;在前锋线以南,冰川得到印度洋季风丰富的降水补给,发育了我国罕有的海洋性冰川类型(施雅风等,2000)。海洋性冰川活动性强,冰川的补给和消融水平高,容易形成冰湖溃决、冰川泥石流等大型灾害。因此,由前锋线划分出的气候分区也决定了冰川灾害分布。
根据我国青藏高原海洋性冰川和大陆型冰川的分布界线(中国科学院青藏高原综合科学考察队,1983),我们发现,在相当一部分地区,这条分界线与STDS在空间上接近甚至重合(图 8),呈现出特提斯喜马拉雅处于干旱寒冷气候区,而高、低喜马拉雅处于温带至亚热带季风气候区的现象。因此,在这些区域选择STDS作为区划线,不仅可以直接体现岩性、气候的差异,还可以反映地貌和灾害类型等更加丰富的信息。
2.5 山地灾害分布特征
现对3个地块山地灾害的形成条件分析如下:
(1) 特提斯喜马拉雅以沉积岩为主,具备崩塌滑坡发生的地形条件;处于寒冷的半干旱高原气候区,常年少雨,冬季降雪量较大,不利于泥石流形成,平均海拔4785 m,山峰终年积雪,在外力作用下容易发生雪崩,成为该地块的主要灾害类型之一。地块处于相对沉陷状态,河流下切速率减缓,铁路百年服务期内山地灾害孕灾环境变化不大。
(2) 高喜马拉雅以硬岩为主,地形最为陡峻,滑坡,特别是崩塌容易发生,处于温带-亚热带半温润季风气候区,具备泥石流形成的降雨条件,其中的海洋性冰川覆盖区容易引发特殊冰川灾害。受地块强烈隆升的影响,河流下切速率加大,崩塌滑坡灾势加剧;尤其在全球气候变暖的背景下,铁路百年服务期内要考虑冰川萎缩导致冰川泥石流、冰湖溃决灾害加剧的问题。
(3) 低喜马拉雅以混杂岩为主,岩性最差,地形陡峻,崩塌,特别是滑坡容易发生;处于亚热带半温润季风气候区,降水量最大,有利于降雨泥石流发生。受地块隆升的影响,山地灾害有增强趋势。
据此得到各地块灾害的主要类型及发展趋势见表 7。
表 7 各地块山地灾害的主要类型及发展趋势Table 7. Main types and trend of mountain hazards in different blocks地块 主要灾害 特殊灾害 灾害发展趋势 特提斯喜马拉雅 崩滑、滑坡 雪崩 平缓 高喜马拉雅 崩塌、滑坡、泥石流 冰川泥石流、冰湖溃决、雪崩 加剧 低喜马拉雅 滑坡、崩塌、降雨泥石流 — 增强 3. 中尼铁路交通廊道工程地质分区
利用对喜马拉雅造山带各地块工程效应的宏观分析结论,以中尼铁路为例,参照铁路大范围方案比选阶段工程地质分区的技术要求,提出一种以构造线作为分区界线的工程地质分区方法。
中尼铁路是“一带一路”倡议中的重要工程,线路起点由我国拉日铁路日喀则站引出,沿国道G318向西布线,经拉孜后折向南,途经萨迦、定结、定日县向南,过吉隆县至吉隆口岸,出境后至尼泊尔首都加德满都,线路全长约666 km,采用吉隆藏布廊道方案穿越喜马拉雅山(姚志勇等,2019)(图 9)。线路沿吉隆藏布河行走,从青藏高原急促下降至恒河平原,是典型的单面大下坡段。现以STDS和MCT作为工程地质分区的界线,将中尼铁路吉隆藏布廊道划分为3个分区类型,并给出分区常用的7个指标,工程地质分区图见图 10,工程地质分区表见表 8。
表 8 中尼铁路吉隆藏布廊道工程地质分区表Table 8. Engineering geology zoning table of China-Nepal Railway along Gyirong river工程地质分段范围及所属地块 地貌特征 地质构造 地层岩性 地震 气象及水文地质 地质灾害 主要工程地质问题 Ⅰ区:吉隆县—隆目特提斯喜马拉雅 高山峡谷区,地形较为陡峻,河流下切程度较强,河流纵坡较陡 构造活动表现为相对下降运动,水平地应力相对较小,主应力方向近东西向 砂岩、砂泥岩、白云岩、灰岩和页岩等沉积岩为主。岩体破裂程度相对较高 地震类型以正断型为主,属高烈度地震区 干旱寒冷气候区。主要接受大气降水和冰川融水补给,石灰岩和白云岩分布区存在岩溶裂隙水 崩塌、滑坡、雪崩 岩性复杂,对南北向的隧道工程,地应力方向与隧道轴线大角度相交,对隧道工程不利 Ⅱ区:隆目—Syabru Bensi高喜马拉雅 极高山峡谷区,地形极为陡峻,河流强烈下切,河流纵坡陡峻 构造活动表现为整体上升运动,水平地应力相对较大,主应力方向近北东向,岩体主要承受挤压应力 以片麻岩、花岗岩和片岩为主。岩体破裂程度相对较低 地震类型以逆冲型为主,地震活动强度大、频率高,对工程不利影响严重 温带气候区,海洋性冰川覆盖区。主要为大气降水和降雪,因地形起伏,切割深度大,一般无法形成大片的补给径流区,高山峡谷地形有利于排水,地下水不丰富 崩塌、滑坡、雪崩;冰川泥石流、冰湖溃决 地震活动性强;冰川泥石流、高地应力、高地温等都可能成为影响线路方案的控制性问题。对南北向的隧道工程,地应力方向与隧道轴线大角度相交,对隧道工程不利;岩石以硬岩为主,岩层受强烈的挤压作用,隧道工程需注意岩爆问题 Ⅲ-1区:SyabruBensi-Betrawati低喜马拉雅 中-高山峡谷区,地形陡峻,河流下切相对较弱,河流纵坡相对较缓 构造活动表现为整体上升运动,水平地应力相对中等,地应力方向近北东向,岩体主要承受挤压应力 页岩、泥岩、石英砂岩及其变质岩等堆积而成的混杂岩。岩体破裂程度相对中等 地震类型以逆冲型为主,地震活动强度大、频率高,对工程不利影响严重 温带至亚热带气候区。主要为大气降水,混杂岩地区构造裂隙含水层、风化带裂隙含水层发育 滑坡、崩塌、降雨泥石流 地震活动性强;高地应力、高地温都可能对隧道工程造成影响。对南北向的隧道工程,地应力方向与隧道轴线大角度相交,对隧道工程不利;软岩分布较多,隧道工程需注意围岩大变形问题 Ⅲ-2区:SyabruBensi-加德满都高喜马拉雅 同Ⅱ区 从表 8来看,3个分区的指标确实反映了同一分区自然条件相近,不同分区差异明显的原则。
根据吉隆藏布流域地貌指标的计算,第Ⅰ区平均起伏度787 m,地表切割深度362 m,坡度21.3°,河流纵坡20.6‰;第Ⅱ区平均起伏度1141 m,地表切割深度525 m,坡度27.3°,河流纵坡35.1‰;第Ⅲ-1区平均起伏度831 m,地表切割深度385 m,坡度22.5°,河流纵坡15.2‰,与2.3节归纳的各地块地貌特征相符。根据中国科学院水利部成都山地灾害与环境研究所的调查资料,第Ⅰ区有崩塌6处,泥石流2条,雪崩10处;第Ⅱ区有崩塌9处,滑坡2处,泥石流8条(其中冰川泥石流4条),雪崩2处;第Ⅲ-1区有崩塌10处,滑坡9处,泥石流7条。根据中国科学院寒区旱区研究所和国际山地综合开发中心(ICIMOD)提供的冰川、冰湖数据,第Ⅰ区有冰湖42个,第Ⅱ区有冰湖212个。这些数据均与表 7归纳的各地块灾害的主要类型相符。可以认为本文提出以构造划分作为工程地质分区的方法具有合理性,可为铁路选线设计提供依据。如可提出3个地块的选线原则:特提斯喜马拉雅地块基本没有对线路方案起控制作用的灾害类型,选线时可尽量利用有利地形,多走明线,降低桥隧比,节省工程造价;高喜马拉雅地块只能采用高桥长隧克服地形障碍,且需注意规避冰川泥石流、冰湖溃决等冰川引发灾害;低喜马拉雅地块气候地质条件均有利于山地灾害发育,必须开展加深地质勘察工作,以规避大型山地灾害为原则选择线路方案。
4. 结论
(1) 本文从造山机制出发,分析了喜马拉雅造山带内动力机制不同的3个地块地层岩性、地震类型、地震活动性、地应力、地貌、灾害分布等工程地质特性。提出以构造划分作为工程地质分区的建议,以中尼铁路为案例的研究表明,该方法具有合理性。
(2) 造山带理论研究的最小单位是由构造划分的地块单元,这代表了造山带理论能够解释的最小空间尺度。本文以构造线作为分区界线,对中尼铁路的分区尺度为几十公里,这种精度能够满足铁路方案比选阶段工程地质分区的要求,表明将造山带理论应用于工程实践,可望取得有意义的成果。
(3) 利用造山带理论可获得常规方法不能获取的信息。如铁路百年服务期内河流演变的趋势、山地灾害的发育趋势,在内动力描述方面,可以判定地震类型、地震相对活动性以及地应力的相对大小等特殊信息。因此,本文方法在铁路领域更具应用前景。
最后,对于新生的造山带,构造作用塑造的地貌占主导成分。若能根据内动力与地貌之间的因果关系,由地貌现象反推内动力状态,获得不使用钻探手段就能得到的地应力信息,是值得进一步深化研究的方向。
致谢: 感谢中国地震局地壳应力研究所、中国科学院寒区旱区环境与工程研究所和国际山地综合开发中心(ICIMOD)在数据方面提供的帮助和支持。感谢中国科学院水利部成都山地灾害与环境研究所的苏立君研究员、张伟峰博士和伍纯昊博士的热情帮助和指教,在此一并感谢。 -
表 1 各地块地震类型占比
Table 1 Proportion of earthquake types in different blocks
地块 正断型/% 正断-走滑型/% 走滑型/% 逆冲型/% 逆冲-走滑型/% 特提斯喜马拉雅 45.5 9.1 40.9 0 0 高喜马拉雅 0 0 34.4 59.4 0 低喜马拉雅 4.0 0 28.0 64.0 4.0 表 2 各地块地震活动性参数
Table 2 Seismic activity parameters of different blocks
地块 地震数目 b值 v4 地震活动性排序 特提斯喜马拉雅 313 1.23 4.15 3 高喜马拉雅 651 1.05 9.18 1 低喜马拉雅 581 1.12 8.79 2 表 3 各地块地震能量密度
Table 3 Seismic energy density of different blocks
地块 面积/km2 总能量/面积/J·km-2 地震能量密度排序 特提斯喜马拉雅 244883.3 60841.67 3 高喜马拉雅 185251.6 430864.20 1 低喜马拉雅 96102.3 423588.49 2 表 4 各地块4级及以上震源深度
Table 4 Focal depths of earthquakes with 4 magnitude and above in different blocks
地块 平均震源深度/km < 20 km地震比例/% >20 km地震比例/% 震源深度排序 特提斯喜马拉雅 30.58 26.09 73.91 3 高喜马拉雅 28.19 35.71 64.29 2 低喜马拉雅 21.92 50.75 49.25 1 表 5 各地块地貌指标
Table 5 Landform indices of different blocks
地块 平均海拔/m 地形起伏度/m 地表切割深度/m 坡度/(°) 特提斯喜马拉雅 4785 647.6 292.9 16.5 高喜马拉雅 3394 1060.8 510.5 25.7 低喜马拉雅 1675 873.9 416.0 24.2 表 6 河流纵坡
Table 6 Longitudinal slopes of rivers
河流名称 河段地块 特提斯喜马拉雅 高喜马拉雅 低喜马拉雅 R1 0.0140 0.0233 0.0038 R2 0.0552 0.0632 0.0073 R3 0.0284 0.0319 0.0090 R4 0.0228 0.0088 0.0089 R5 0.0500 0.0350 0.0066 R6 0.0475 0.0344 0.0087 R7 0.0206 0.0351 0.0152 R8 0.0199 0.0551 0.0054 R9 0.0021 0.0228 0.0131 平均值 0.0289 0.0344 0.0087 表 7 各地块山地灾害的主要类型及发展趋势
Table 7 Main types and trend of mountain hazards in different blocks
地块 主要灾害 特殊灾害 灾害发展趋势 特提斯喜马拉雅 崩滑、滑坡 雪崩 平缓 高喜马拉雅 崩塌、滑坡、泥石流 冰川泥石流、冰湖溃决、雪崩 加剧 低喜马拉雅 滑坡、崩塌、降雨泥石流 — 增强 表 8 中尼铁路吉隆藏布廊道工程地质分区表
Table 8 Engineering geology zoning table of China-Nepal Railway along Gyirong river
工程地质分段范围及所属地块 地貌特征 地质构造 地层岩性 地震 气象及水文地质 地质灾害 主要工程地质问题 Ⅰ区:吉隆县—隆目特提斯喜马拉雅 高山峡谷区,地形较为陡峻,河流下切程度较强,河流纵坡较陡 构造活动表现为相对下降运动,水平地应力相对较小,主应力方向近东西向 砂岩、砂泥岩、白云岩、灰岩和页岩等沉积岩为主。岩体破裂程度相对较高 地震类型以正断型为主,属高烈度地震区 干旱寒冷气候区。主要接受大气降水和冰川融水补给,石灰岩和白云岩分布区存在岩溶裂隙水 崩塌、滑坡、雪崩 岩性复杂,对南北向的隧道工程,地应力方向与隧道轴线大角度相交,对隧道工程不利 Ⅱ区:隆目—Syabru Bensi高喜马拉雅 极高山峡谷区,地形极为陡峻,河流强烈下切,河流纵坡陡峻 构造活动表现为整体上升运动,水平地应力相对较大,主应力方向近北东向,岩体主要承受挤压应力 以片麻岩、花岗岩和片岩为主。岩体破裂程度相对较低 地震类型以逆冲型为主,地震活动强度大、频率高,对工程不利影响严重 温带气候区,海洋性冰川覆盖区。主要为大气降水和降雪,因地形起伏,切割深度大,一般无法形成大片的补给径流区,高山峡谷地形有利于排水,地下水不丰富 崩塌、滑坡、雪崩;冰川泥石流、冰湖溃决 地震活动性强;冰川泥石流、高地应力、高地温等都可能成为影响线路方案的控制性问题。对南北向的隧道工程,地应力方向与隧道轴线大角度相交,对隧道工程不利;岩石以硬岩为主,岩层受强烈的挤压作用,隧道工程需注意岩爆问题 Ⅲ-1区:SyabruBensi-Betrawati低喜马拉雅 中-高山峡谷区,地形陡峻,河流下切相对较弱,河流纵坡相对较缓 构造活动表现为整体上升运动,水平地应力相对中等,地应力方向近北东向,岩体主要承受挤压应力 页岩、泥岩、石英砂岩及其变质岩等堆积而成的混杂岩。岩体破裂程度相对中等 地震类型以逆冲型为主,地震活动强度大、频率高,对工程不利影响严重 温带至亚热带气候区。主要为大气降水,混杂岩地区构造裂隙含水层、风化带裂隙含水层发育 滑坡、崩塌、降雨泥石流 地震活动性强;高地应力、高地温都可能对隧道工程造成影响。对南北向的隧道工程,地应力方向与隧道轴线大角度相交,对隧道工程不利;软岩分布较多,隧道工程需注意围岩大变形问题 Ⅲ-2区:SyabruBensi-加德满都高喜马拉雅 同Ⅱ区 -
Brookfield M E. 1993. The Himalayan passive margin from Precambrian to Cretaceous times[J]. Sedimentary Geology, 84(1-4): 1-35. DOI: 10.1016/0037-0738(93)90042-4
Burchfiel B C, Royden L H. 1985. North-south extension within the convergent Himalayan region[J]. Geology, 13(10): 679-682. DOI: 10.1130/0091-7613(1985)13<679:NEWTCH>2.0.CO;2
Burchfield B C, Chen Z, Hodges K V, et al. 1992. The South Tibetan Detachement System, Himalayan orogen: extension contemporaneous with and parallel to shortening in a collisional mountain belt[J]. Geological Society of America Special Paper, 269(1): 1-41.
Burg J P, Chen G M. 1984. Tectonics and structural formation of southern Tibet, China[J]. Nature, 311(1): 219-223. DOI: 10.1038-311219a0/
Chen P N, Gao L Q. 1987. Characteristics of horizontal stress distribution and their relation to seismicity[J]. Earthquake Research in China, 3 (S): 85-89.
Ding L, Maksatbek S, Cai F L, et al. 2017. Processes of initial collision and suturing between India and Asia[J]. Science China Earth Sciences: 47(3)293-309.
Gansser A. 1964. Geology of the Himalayas[M]. New York: Wiley Inter Science: 289.
Gao M T. 2015. Seismic parameter zonation map of China(GB 18306-2015)[S]. Beijing: Standards Press of China.
Godin L, Grujic D, Law R D, et al. 2006. Channel flow, ductile extrusion and exhumation in continental collision zones[M]. London: Geological Society Special Publication, 268: 1-23.
Gutenberg B, Richter C F. 1956. Magnitude and energy of earthquakes[J]. Annali di Geofisica, 9(1): 1-15.
Harrion T M, Copeland P, Kidd W S, et al. 1992. Raising Tibet[J]. Science, 255(5052): 1663-1670. DOI: 10.1126/science.255.5052.1663
Huang R Q, Qi S W. 2017. Engineering geology: Review and prospect of past ten years in China[J]. Journal of Engineering Geology, 25(2): 257-276.
Jackson M, Bilham R. 1994. Constraints on Himalayan deformation inferred from vertical velocity fields in Nepal and Tibet[J]. Journal of Geophysical Research, 99 (B7): 13897-13912. DOI: 10.1029/94JB00714
Lefort P. 1996. Evolution of the Himalaya, in the tectonics of Asia[M]. New York: Cambridge University Press: 95-106.
Li J J, Wen S X, Zhang Q S, et al. 1979. A discussion on the period, amplitude and type of the uplift of the Qinghai-Xizang plateau[J]. Scientia Sinica, 22(11): 1314-1328.
Lin Z Y, Wu X D. 1981. Climatic regionalization of the Qinghai-Xizang Plateau[J]. Acta Geographica Sinica, 36(1): 22-32.
Mogi K. 1962. Study of elastic shocks caused by the fracture of heterogeneous materials and its relations to earthquake phenomena[J]. Bulletin of the Earthquake Research Institute, 40(1): 125-173.
Peng J B. 2001. Theory-method system in study of dynamics of the regional stability[J]. Journal of Engineering Geology, 9(1): 3-11.
Peng J B. 2001. Fundamental principles for study on dynamics of regional stability[J]. Journal of Engineering Geology, 9(1): 12-16.
Peng J B, Mao Y L, Fan Wen. 2001. Study of dynamics on regional stability[M]. Beijing: Science Press.
Scheidegger A E. 1979. The principle of antagonism in the Earth's evolution[J]. Tectonophysics, 55:T7-10. DOI: 10.1016/0040-1951(79)90177-X
Scheidegger A E. 2004. Morphotectonics[M]. Springer: New York: 111-118.
Shi Y F, Huang M H, Yao T D, et al. 2000. Glaciers and their environments in China[M]. Beijing: Science Press: 368-370.
Tapponnie P, Peltzer G, Armijo R. 1986. On the mechanics of the collision between India and Asia[J]. Geological Society of London Special Publications, 19(1): 113-157. DOI: 10.1144/GSL.SP.1986.019.01.07
The Comprehensive Scientific Expedition to the Qinghai-Xizang Plateau, Chinese Academy of Sciences. 1983. Glaciers of Xizang(Tibet)[M]. Beijing: Science Press: 6-36.
Wang S J. 2013. Century evolution of engineering geo-sciences and further prospect[J]. Journal of Engineering Geology, 21(1): 1-5.
Xu Q. 2012. New development tendency of engineering geology in China—summary report of the 9th national congress of engineering geology[J]. Journal of Engineering Geology, 20(6): 1087-1095.
Xu Z Q, Wang Q, Zeng L S, et al. 2013. Three-dimensional extrusion model of the Great Himalaya slice[J]. Geology in China, 40(3): 671-680.
Xu Z Q, Yang J S, Li H B, et al. 2007. Orogenic plateaux—Terrane amalgamation, collision and uplift in the Qinghai-Tibet Plateau[M]. Beijing: Geological Pubishing House: 259-263.
Xu Z H. 1990. Current status and prospects of geostress[J]. Advances in Earth Science, (5): 27-34.
Yao L K, Qiu Y L, Wei Y X. 2012. Challenges in construction of railway and highway from Sichuan to Tibet through eastern margin of Tibetan Plateau[J]. Journal of Southwest Jiaotong University, 47(5): 719-734.
Yu Q Y, Jiang Z Z, Ai N S. 1983. A new method for calculating water system direction data[J]. Journal of Lanzhou University, 19(3): 164-165.
Yao Z Y, Meng X L, Li X, et al. 2019. Research on engineering geological environment characteristics and engineering geological problems of China-Nepal cross-border railway[J]. Railway Standard Design, 63(7): 28-34.
Yang B C, Qin S Q, Xue L, et al. 2017. A physical self-similarity law describing the accelerated failure behavior of rocks[J]. Chinese Journal of Geophysics, 60(5): 1746-1760.
Yin A. 2006. Cenozoic tectonic evolution of the Himalayan orogen as constrained by along-strike variation of structural geometry, exhumation history, and foreland sedimentation[J]. Earth-Science Reviews, 76(1-2): 1-131. DOI: 10.1016/j.earscirev.2005.05.004
Yin A. 2006. Cenozoic tectonic evolution of the Himalayan orogen as constrained by along-strike variation of structural geometry, exhumation history, and foreland sedimentation[J]. Earth Science Frontiers, 13(5): 416-515.
Zhang Y S, Qu Y X, Wu S R, et al. 2004. The research introduction of engineering geology in Qinling orogenic belt[J]. Journal of Engineering Geology, 12 (S): 44-49.
陈彭年, 高莉青. 1987.地应力水平分布特征及其与地震的关系[J].中国地震, 3 (S): 85-89. 高孟潭. 2015. 《中国地震动参数区划图》宣贯教材(GB 18306-2015)[S].北京: 中国标准出版社. 黄润秋, 祁生文. 2017.工程地质:十年回顾与展望[J].工程地质学报, 25(2): 257-276. DOI: 10.13544/j.cnki.jeg.2017.02.001 林振耀, 吴祥定. 1981.青藏高原气候区划[J].地理学报, 36(1): 22-32. DOI: 10.3321/j.issn:0375-5444.1981.01.003 彭建兵. 2001a.区域稳定动力学研究(一)理论与方法体系[J].工程地质学报, 9(1): 3-11. 彭建兵. 2001b.区域稳定动力学研究(二)基本原则[J].工程地质学报, 9(1): 12-16. 彭建兵, 毛彦龙, 范文. 2001.区域稳定动力学研究[M].北京:科学出版社: 78-89. 施雅风, 黄茂桓, 姚檀栋, 等. 2000.中国冰川与环境[M].北京:科学出版社: 368-370. 王思敬. 2013.工程地质学科的世纪演化与前景[J].工程地质学报, 21(1): 1-5. DOI: 10.3969/j.issn.1004-9665.2013.01.001 许强. 2012.工程地质学科发展的新趋势-第九届全国工程地质大会学术总结[J].工程地质学报, 20(6): 1087-1095. DOI: 10.3969/j.issn.1004-9665.2012.06.025 许志琴, 杨经绥, 李海兵, 等. 2007.造山的高原-青藏高原的地体拼合、碰撞造山及隆升机制[M].北京:地质出版社: 259-263. 许志琴, 王勤, 曾令森, 等. 2013.高喜马拉雅的三维挤出模式[J].中国地质, 40(3): 671-680. DOI: 10.3969/j.issn.1000-3657.2013.03.002 许忠淮. 1990.地应力研究现状与展望[J].地球科学进展, (5): 27-34. 姚令侃, 邱燕玲, 魏永幸. 2012.青藏高原东缘进藏高等级道路面临的挑战[J].西南交通大学学报, 47(5): 719-734. DOI: 10.3969/j.issn.0258-2724.2012.05.001 余庆余, 蒋柱中, 艾南山. 1983.水系方向资料计算的一种新方法[J].兰州大学学报(自然科学版), 19(3): 164-165. 姚志勇, 孟祥连, 李响, 等. 2019.中尼跨境铁路工程地质环境特征及主要工程地质问题研究[J].铁道标准设计, 63(7): 28-34. 杨百存, 秦四清, 薛雷, 等. 2017.岩石加速破裂行为的物理自相似律[J].地球物理学报, 60(5): 1746-1760. 张永双, 曲永新, 吴树仁, 等. 2004.秦岭造山带工程地质研究导论[J].工程地质学报, 12(增刊): 44-49. 中国科学院青藏高原综合科学考察队. 1983.西藏冰川[M].北京:科学出版社: 6-36.