龍門(mén)山北段地貌發(fā)育特征及其形成因素分析
——以通口河治城段流域?yàn)槔?/h1>
2018-04-28 02:21:09馮金煒邱占林劉玉法
馮金煒, 馬 超, 李 勇, 陳 歡, 邱占林, 劉玉法, 張 威, 趙 重
( 1. 成都理工大學(xué) 油氣藏地質(zhì)及開(kāi)發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室, 四川 成都 610059； 2. 四川省地震局 工程地震研究院, 四川 成都 610041； 3. 中國(guó)地質(zhì)科學(xué)院 探礦工藝研究所, 四川 成都 610031； 4. 龍巖學(xué)院 資源工程學(xué)院, 福建 龍巖 364012)
地貌是不斷發(fā)展變化的,地貌發(fā)展變化的物質(zhì)過(guò)程稱(chēng)地貌過(guò)程,內(nèi)力和外力是塑造地貌的2種營(yíng)力,地貌是內(nèi)力過(guò)程與外力過(guò)程相互統(tǒng)一的產(chǎn)物.龍門(mén)山是青藏高原山脈中落差高度變化最大的山脈,構(gòu)造作用強(qiáng)烈,李勇等[1]通過(guò)研究龍門(mén)山構(gòu)造的地貌標(biāo)志定量計(jì)算了龍門(mén)山主斷裂的逆沖速率與走滑速率的關(guān)系,進(jìn)一步分析了龍門(mén)山周緣地貌演化趨勢(shì).Ouchi[2]認(rèn)為斷層及褶皺影響河流形態(tài)及搬運(yùn)沉積過(guò)程,而河流的多樣性則影響著周緣的地貌的演化.李亞林等[3]通過(guò)研究構(gòu)造運(yùn)動(dòng)對(duì)河流地貌特征的響應(yīng)指出：早期逆滑和走滑斷裂控制水系形態(tài),而晚期的正斷層和地塹構(gòu)造控制河流地貌的發(fā)育演變.2008年“5·12”地震導(dǎo)致映秀—北川斷裂發(fā)生地表破裂,導(dǎo)致龍門(mén)山地區(qū)基巖的同震隆升變形[1,3],觸發(fā)大規(guī)模的同震滑坡、崩塌等地質(zhì)災(zāi)害發(fā)生,形成大面積松散物質(zhì)堆積,在氣候的作用下進(jìn)一步改造龍門(mén)山區(qū)域的地貌[4-6].因此研究區(qū)域性的構(gòu)造抬升及河流侵蝕對(duì)未來(lái)地貌演化方向的影響十分重要.本文基于ArcGis技術(shù)提取龍門(mén)山北段通口河治城段的地形因子及河流剖面等參數(shù),進(jìn)一步探討構(gòu)造、氣候等因素對(duì)治城段地貌演化發(fā)展趨勢(shì)的影響.
1 資料與處理方法
通口河系屬龍門(mén)山中段山前河流,是涪江的重要支流之一,起源于九頂山南麓,西北—南東流向,在彭州市境內(nèi)穿越龍門(mén)山,在丹景山鎮(zhèn)進(jìn)入成都平原,是涪江流域的一級(jí)支流.本文研究區(qū)域通口河段挾持于后龍門(mén)山深斷裂與江油—灌縣大斷裂之間,最高海拔2 260 m,最低海拔600 m(圖1),落差1 660 m,地勢(shì)呈現(xiàn)自北西往南東降低(圖2).
本次研究通過(guò)無(wú)人機(jī)勘察及野外實(shí)測(cè)獲得高程點(diǎn),由高程點(diǎn)繪制成等高線(xiàn).利用ArcGis 10.2軟件將等高線(xiàn)矢量數(shù)據(jù)生成TIN,再由TIN轉(zhuǎn)成DEM,經(jīng)過(guò)DEM柵格數(shù)據(jù)多次處理獲得通口河流域治城段地形特征和水系網(wǎng)絡(luò)特征.
1.1地形坡度的獲取利用ArcGis 10.2軟件,將獲得的高分辨率DEM數(shù)據(jù)(1∶1萬(wàn),柵格大小5 m),經(jīng)過(guò)ArcToolbox中3D Analyst工具的柵格表面計(jì)算地面坡度的平均值,獲得治城段平均坡度分布區(qū)間為0°～69.1°,通過(guò)圖層屬性將平均坡度范圍分為4個(gè)等級(jí),分別為：0°～15.6°、15.6°～30.9°、30.9°～44.5°、44.5°～69.1°.4個(gè)坡度區(qū)間的分布范圍如圖2所示.
1.2水系網(wǎng)絡(luò)的獲取利用ArcGis 10.2軟件,將得到的DEM進(jìn)行填洼獲得無(wú)洼DEM,經(jīng)過(guò)水文分析工具得到流向和流量,將閾值設(shè)置為5 000獲得流水累計(jì)柵格,以此得到通口河流域水系分布圖(圖3).通口河流域主要發(fā)育4個(gè)較大支流：小園河支流(支流1)、白草河支流(支流2)、白坭溝支流(支流3)、玉溪河支流(支流4).下面將通過(guò)對(duì)通口河支流的河流縱剖面和橫剖面特征進(jìn)一步分析.
圖 1 通口河流域縱剖面圖
圖 2 治城段流域平均坡度分布圖
圖 3 治城段流域地形與水系圖Fig. 3 Topographic and water system of Zhicheng river basin
2 通口河治城段流域發(fā)育特征
2.1通口河治城段流域的地形特征通過(guò)DEM計(jì)算出來(lái)的通口河流域治城段平均坡度分布圖(圖2)可宏觀(guān)表達(dá)該區(qū)域地貌特征.北川老縣城附近(即距北川—映秀斷裂北西側(cè)4 km區(qū)域內(nèi))平均坡度較大,多介于30.9°～62.3°；玉溪河和通口河河流區(qū)域坡度較小,多為10°～30°；研究區(qū)北西段(以禹里鄉(xiāng)和漩坪鄉(xiāng)連線(xiàn)為界,下同)劉家窩—茶園坪—白巖坪—燈草坪斷斷裂破碎帶區(qū)域平均坡度較大,多介于30.9°～62.3°；挾于各斷裂帶之間的區(qū)域區(qū)域均屬于小角度坡度：16.3°～26.0°.
不同地層巖性(圖4)可體現(xiàn)出不同坡度分布,結(jié)合地質(zhì)圖與坡度圖可進(jìn)一步分析通口河流域的地貌特征.北川映秀斷裂北西向(4 km)為志寒武系油房組,巖性以中、粗粒巖屑砂巖、粉砂巖互層,南東側(cè)以泥盆系和石炭系灰?guī)r為主,北西向一側(cè)坡度明顯較大.研究區(qū)中部以志留系茂縣群第二組灰?guī)r環(huán)繞茂縣群第一組千枚巖,灰?guī)r區(qū)域坡度大于千枚巖區(qū)域.位于研究區(qū)北段劉家窩—茶園坪—白巖坪—燈草坪斷層附近巖性以千枚巖夾灰?guī)r、砂巖為主,坡度大于單純千枚巖地層.
通過(guò)斷裂帶與巖性對(duì)坡度分布影響可知,斷裂(劉家窩—茶園坪—白巖坪—燈草坪斷裂、北川斷裂)上下盤(pán)周?chē)钠露染哂幸欢ǖ牟町愋?斷層上盤(pán)4～6 km范圍內(nèi)坡度大,說(shuō)明2條逆沖斷裂上盤(pán)處于隆升狀態(tài),斷裂上盤(pán)的河流侵蝕下切作用將相對(duì)較快.通過(guò)對(duì)比地形坡度(圖2)及地層巖性(圖4)可知巖性的“軟硬”程度同樣對(duì)坡度的分布起一定的作用：研究表明砂巖和灰?guī)r區(qū)域坡度明顯大于千枚巖區(qū)域坡度.
2.2通口河治城段流域的水系特征本文基于ArcGis 10.2軟件,將得到的DEM進(jìn)行填洼或得無(wú)洼DEM,經(jīng)過(guò)水文分析工具得到流向和流量,將閾值設(shè)置為5 000獲得流水累計(jì)柵格,以此得到通口河流域水系分布圖(圖2).通口河流域主要發(fā)育4個(gè)較大支流：小園河支流(支流1)、白草河支流(支流2)、白坭溝支流(支流3)、玉溪河支流(支流4).下面通過(guò)對(duì)通口河支流的河流縱剖面和橫剖面特征進(jìn)一步分析.
圖4治城段流域地質(zhì)構(gòu)造簡(jiǎn)圖
Fig.4ThesimplifiedgeologicalandstructuralmapofZhicehngriverbasin
2.2.1河流縱剖面 受到地質(zhì)構(gòu)造、流域基部巖性特性、氣候降雨等諸多因素影響,河流縱剖面呈現(xiàn)不同形態(tài).研究區(qū)域較小,氣候條件幾乎相同；基部巖性以志留系茂縣群千枚巖為主,4條支流流向大致垂直于巖層走向,因而該區(qū)域河流剖面形態(tài)主受構(gòu)造影響.通過(guò)對(duì)河流縱剖面形態(tài)特征的分析可有效反映通口河流域的構(gòu)造差異性.
關(guān)于河流剖面與地貌演化的關(guān)系前人做出了大量的研究,陳彥杰等[7]通過(guò)對(duì)河流縱剖面模型進(jìn)行了總結(jié)：1)地形面在受到地體抬升后構(gòu)造運(yùn)動(dòng)即停止,則初始形成的河流其縱剖面的下凹程度較小或近乎直線(xiàn), 可以用線(xiàn)性剖面來(lái)擬合；2)隨著中上游侵蝕作用的進(jìn)行,河流將侵蝕下來(lái)的物質(zhì)搬運(yùn)至下游堆積,使得河流縱剖面的下凹程度逐漸變大而逐漸演化成指數(shù)剖面；3)隨著下游物質(zhì)的堆積,河流中-上游的侵蝕,河流縱剖面下凹程度進(jìn)一步加大,可以用對(duì)數(shù)函數(shù)來(lái)擬合河流縱剖面;4)再進(jìn)一步發(fā)展,如果流量增加或由于構(gòu)造運(yùn)動(dòng)等形成破碎帶導(dǎo)致輸沙量增加,使得河流縱剖面的下凹程度快速增大而呈現(xiàn)出乘冪函數(shù)剖面.河流縱剖面的演化趨勢(shì)為:線(xiàn)性→指數(shù)→對(duì)數(shù)→乘冪.
通口河流域治城段圖幅較小,在考慮到流域內(nèi)各條支流的氣候和巖性類(lèi)似的情況下,構(gòu)造運(yùn)動(dòng)是影響河流縱剖面的主控因素.通過(guò)對(duì)研究區(qū)4條支流(圖2、圖5)對(duì)比分析：支流1、支流2、支流3下凹程度低,呈近線(xiàn)性指數(shù)擬合,即考慮該3條支流較“年輕”,在氣候和基底巖性相似的情況下,支流1靠近劉家窩—茶園坪—白巖坪—燈草坪斷層；支流2穿過(guò)劉家窩—茶園坪—白巖坪—燈草坪斷層、盛家院—密地坪—楊柳松斷層及茂汶斷裂；支流3位于盛家院—密地坪—楊柳松斷層附近；很可能由于構(gòu)造導(dǎo)致山體相對(duì)抬升,使得這3條支流的縱剖面接近線(xiàn)性擬合.支流4和通口河主干接近對(duì)數(shù)剖面,說(shuō)明這2條河流演化速度慢受地質(zhì)構(gòu)造作用相對(duì)較少,河流進(jìn)入“中、老年化”,根據(jù)地質(zhì)資料及實(shí)地路線(xiàn)觀(guān)察,該2條河流附近斷層地質(zhì)作用較弱,處于相對(duì)穩(wěn)定狀態(tài).
圖 5 治城流域主要支流縱剖面
2.2.2河流橫剖面 通過(guò)研究區(qū)橫剖面線(xiàn)可得各支流情況(圖6)：支流1、支流2、支流3河谷呈深“V”型,說(shuō)明河流正進(jìn)行快速下切侵蝕作用,也進(jìn)一步反映這3條支流范圍內(nèi)抬升作用較強(qiáng).支流3和通口河(湔江)主干河谷為U型谷,說(shuō)明這2條河流有側(cè)向侵蝕作用.根據(jù)河流發(fā)育特征,河流先以縱向下切作用為主,當(dāng)河流縱剖面達(dá)到平衡狀態(tài)時(shí),河流再以側(cè)向侵蝕為主.通過(guò)5條河谷對(duì)比分析：支流4和通口河主流相對(duì)于支流1、2、3下切侵蝕速率慢,構(gòu)造相對(duì)穩(wěn)定,山體抬升不明顯,地貌不會(huì)發(fā)生太劇烈的演化.
圖 6 流域治城段河流橫剖面分析圖
3 通口河治城段流域地貌發(fā)育的主要影響因素
3.1構(gòu)造活動(dòng)對(duì)地貌演化的影響2008年汶川地震的發(fā)生導(dǎo)致震區(qū)地質(zhì)災(zāi)害大面積爆發(fā),據(jù)野外核查,研究區(qū)構(gòu)造活動(dòng)的較明顯(斷層附近)地區(qū)地質(zhì)災(zāi)害成面狀分布.在地質(zhì)構(gòu)造的發(fā)育的地區(qū),構(gòu)造面發(fā)育導(dǎo)致周?chē)鷰r體破碎,穩(wěn)定性差,極易發(fā)生地質(zhì)災(zāi)害.構(gòu)造對(duì)地質(zhì)災(zāi)害的控制主要表現(xiàn)在以下幾方面：1) 斷層均以北東向排布,控制著研究區(qū)內(nèi)山脈及河流的展布.區(qū)域內(nèi)多條小斷層交匯,交匯處形成大面積“X”節(jié)理,使得巖體失穩(wěn),形成規(guī)模較大的地質(zhì)災(zāi)害.斷裂帶附近巖體較其它地區(qū)破碎,研究區(qū)西北段無(wú)根密集斷層、中段茂汶斷裂、南西段北川斷裂比較典型.2) 現(xiàn)代新構(gòu)造運(yùn)動(dòng),特別是地震活動(dòng),直接破壞了巖土體的結(jié)構(gòu),使巖土體的物理力學(xué)強(qiáng)度降低.研究區(qū)主要以絹云千枚巖為主,巖性較為軟弱,在斷層作用下,節(jié)理面發(fā)育,巖體破碎,當(dāng)?shù)乇硭疂B入時(shí)使巖層軟化形成活動(dòng)帶,當(dāng)遇大量降雨滲入土體,增加了土體自重或下滑力,使土體的粘聚力和磨擦力減小,便產(chǎn)生了滑坡、崩塌,也為泥石流提供大量物源.3) 構(gòu)造帶中多期次構(gòu)造運(yùn)動(dòng)的疊加,使巖體完整性降低,更易于被風(fēng)化,同時(shí)由于破碎帶的影響,其風(fēng)化深度更大,研究區(qū)多條區(qū)域性斷裂為后期在北川斷層上的疊加.4) 構(gòu)造裂隙及斷層面常與卸荷裂隙組成向臨空面的不利組合.
2008年特大地震后不僅造成龍門(mén)山斷裂帶急劇的隆升,也進(jìn)一步觸發(fā)了大面積的地質(zhì)災(zāi)害的爆發(fā).根據(jù)2016年9月野外勘察,共查明研究區(qū)域內(nèi)崩塌12處、滑坡56處、泥石流5處,地質(zhì)災(zāi)害面積達(dá)到5.6 km2,松散堆積物體積約0.019 km3.2010年,對(duì)龍門(mén)山北川—映秀斷裂帶同震巖石隆升研究得出的數(shù)據(jù)中,該研究區(qū)域增加的山脈體積估算值約為0.176 km3,大于松散堆積物體積.Robert等[8]也指出,沿?cái)鄬臃较蛏系牟糠謪^(qū)域同震巖石隆升增加的山脈體積大于地震導(dǎo)致的堆積物體積,即山脈體積增加.研究區(qū)北西段發(fā)育小型密集“無(wú)根斷層”,地震后巖石隆升作用較為劇烈(圖7),雖然同震地質(zhì)災(zāi)害產(chǎn)生了松散堆積體,但研究區(qū)基巖隆升體積大于同震松散物質(zhì)體積,即研究區(qū)山體體積并未“虧損”,而是處于抬升狀態(tài).
圖 7 汶川地震觸發(fā)的同震抬升圖(據(jù)文獻(xiàn)[8]修改)
3.2氣候?qū)Φ孛惭莼挠绊懲诤恿饔蛑纬嵌螝夂驅(qū)儆趤啛釒貪?rùn)疾風(fēng)氣候,降雨量充足,1967年為最大年降雨量,達(dá)到2 340 mm,日最大降雨量達(dá)到101 mm,年平均降雨量可達(dá)1 399.1 mm.強(qiáng)降雨時(shí)期主要集中在7、8、9月份(圖8),達(dá)到全年降雨量的70%以上.2008年汶川后大面積地質(zhì)災(zāi)害爆發(fā),2013年“7.9”洪災(zāi)使得一些潛在地質(zhì)災(zāi)害發(fā)生,暴雨誘發(fā)的泥石流、崩塌、滑坡共31處.暴雨導(dǎo)致的山洪峰值較大,攜帶溝谷內(nèi)崩塌和滑坡物質(zhì)進(jìn)入河流.研究區(qū)域較大的泥石流有倒開(kāi)門(mén)泥石流和羊盤(pán)溝泥石流,每年雨季均攜帶大量松散物質(zhì)分別匯入通口河和白草河,地震導(dǎo)致的崩塌物質(zhì)是主要的物源.許強(qiáng)通過(guò)實(shí)驗(yàn),表明：地震災(zāi)區(qū)發(fā)生大型泥石流的臨界雨降為35～40 mm/h,最低僅為15 mm/h[9].
根據(jù)實(shí)地野外核查,泥石流位置均靠近河流,由于構(gòu)造(主要地震)產(chǎn)生的大面積松散物質(zhì)通過(guò)河流作用搬運(yùn)至下游地區(qū).由于山脈物質(zhì)的遷移,影響地貌的改變.氣候(以降水為主)影響地貌演化的最直接表現(xiàn)便是河流中輸砂率的變化.
研究區(qū)松散物通過(guò)洪流由白草河、白泥溝、玉溪河等支流攜帶搬運(yùn)至通口河,根據(jù)通口河下游將軍水文站近十年數(shù)據(jù)表明,降雨量直接影響輸砂量,7～9月份為研究區(qū)雨季,輸砂量野達(dá)到峰值.圖8體現(xiàn)出降雨量和輸砂率呈正相關(guān)關(guān)系.
圖 8 1-12月通口河流域平均輸砂率與北川縣的月降水量對(duì)比分析
前人通過(guò)流域輸砂量的計(jì)算獲來(lái)得剝蝕速率,本文通過(guò)輸砂率估算出剝蝕速率,以此研究流域內(nèi)地貌的演化發(fā)展.據(jù)綿陽(yáng)將軍石水文站實(shí)測(cè)懸砂數(shù)據(jù)統(tǒng)計(jì)[10]：多年平均輸砂率171 kg/s,多年平均輸砂量1.8 kg/m3,多年平均汛期含砂量為2.25 kg/m3,最大輸砂量為320萬(wàn)t(1964年),最小年輸砂量96萬(wàn)t(1974年).據(jù)此比擬推求通口河流域年平均輸砂量為2.27×106t,沉積物密度取1.50 t/m3,即1.51×106m3.通口河流域面積為2 808 km2,計(jì)算出通口河流域剝蝕速率為0.537 mm/a.侵蝕速率應(yīng)包括風(fēng)力搬運(yùn)作用、化學(xué)作用等,因此通口河流域的剝蝕速率應(yīng)遠(yuǎn)大于0.537 mm/a.研究數(shù)據(jù)表明該區(qū)域平均隆升幅度為0.456 m,即在未來(lái)850年內(nèi)將剝蝕掉地震成因的山脈增加量. 在地質(zhì)構(gòu)造不太強(qiáng)雷的情況下,在氣候因素的影響下,山脈物質(zhì)增加量將被河流遷移至下游.
前人對(duì)龍門(mén)山在非地震期間的隆升速率已經(jīng)進(jìn)行了大量研究.劉樹(shù)根[11]根據(jù)測(cè)量地表的變形機(jī)制獲得龍門(mén)山逆沖推覆帶的年均隆升速率在0.3～0.4 mm之間；江娃利等[12]研究指出龍門(mén)山構(gòu)造帶可達(dá)0.2～2.5 m的垂直位移量；邵崇建等[13]等根據(jù)多年(1970—1997年)野外實(shí)測(cè)資料得出在理縣附近的隆升速率最大,可達(dá)到 3. 5 mm/a.龍門(mén)山處于非地震階段的垂直隆升速率目前數(shù)據(jù)值較多,大致處于0.2～3.5 mm/a范圍內(nèi).對(duì)比通口河的侵蝕速率(>0.537 mm/a)與隆升速率(0.2～3.5 mm/a),區(qū)域總體呈現(xiàn)出抬升作用,海拔處于上升階段.
4 地貌演化分區(qū)
綜合上述坡度、河網(wǎng)、構(gòu)造、氣候等要素對(duì)地貌演化的影響進(jìn)行分析,對(duì)通口河流域治城段劃分為兩個(gè)區(qū)域,以此討論地貌演化過(guò)程.
4.1Ⅰ區(qū)(緩慢抬升區(qū)) Ⅰ區(qū)位于北川斷裂上盤(pán)附近,汶川地震使得該區(qū)山脈抬升.通過(guò)各支流縱坡面(圖5)可知位于該處的通口河段具有破折點(diǎn)和匯入通口河的玉溪河(支流4)具有明顯的破折點(diǎn),說(shuō)明該處正處于抬升狀態(tài),并且正經(jīng)歷河流下切作用,但由河流橫剖面(圖6)反應(yīng)出這2條河流河谷較為寬闊,即河流下切作用弱.此外,該處Ⅰ區(qū)巖性由志留系千枚巖過(guò)渡到較為堅(jiān)硬的砂巖,進(jìn)一步制約了河流剝蝕速率.Ⅰ區(qū)堰塞湖位置也可能淤積上流搬運(yùn)下來(lái)的河流物質(zhì),使得松散物質(zhì)堆積量增加.Ⅰ區(qū)剝蝕速率小于基巖抬升速率,可以推斷Ⅰ區(qū)可能為緩慢隆升的地貌演化趨勢(shì).
4.2Ⅱ區(qū)(快速抬升區(qū)) Ⅱ區(qū)為研究區(qū)北西段,該處發(fā)育大量“無(wú)根密集”斷層,巖性較為破碎,表現(xiàn)為構(gòu)造作用強(qiáng)烈.支流1、2、3縱坡面呈深“V”型,縱剖面為線(xiàn)性指數(shù),即該處河流下切速度快,進(jìn)一步說(shuō)明該區(qū)屬于隆升階段.此外,由地形坡度(圖3)可知,位于支流1、2、3附近坡度大,河流兩側(cè)的將形成陡峭的河谷地形.雖然在汶川地震觸發(fā)下形成大面積的崩塌、滑坡等松散堆積物質(zhì),在河流作用下將被搬運(yùn)到下流,但通過(guò)對(duì)比隆升幅度及侵蝕速率,山體并未“虧損”.Ⅱ區(qū)基巖抬升速率大于河流侵蝕作用,且在較為劇烈構(gòu)造運(yùn)動(dòng)的情況下,地貌演化將呈快速抬升趨勢(shì).
圖 9 治城段流域地貌演化分區(qū)
5 結(jié)束語(yǔ)
通過(guò)對(duì)龍門(mén)山北段通口河流域地貌、構(gòu)造、氣候等特征進(jìn)行分析,初步得下以下幾點(diǎn)認(rèn)識(shí).
1) 斷層上盤(pán)附近坡度普遍較大,區(qū)域處于抬升狀態(tài),河流侵蝕下切作用強(qiáng),使得地形變得陡峭.巖性影響地形特征,砂巖和灰?guī)r區(qū)域地形坡度明顯大于千枚巖區(qū)域坡度.
2) 以如今通口河流域治城段的構(gòu)造作用和氣候特征,大部分河流處于“年輕”狀態(tài),下切剝蝕速率較快,搬運(yùn)能力較強(qiáng),上游沉積物極易搬運(yùn)到下游,使得流域局部地貌往削低趨勢(shì)發(fā)展.
3) 南東段(北川斷裂上盤(pán)附近)受構(gòu)造作用影響處于抬升階段,該位置河流河谷較為寬闊,剝蝕山體速率較慢,受堰塞湖等因素影響將繼續(xù)積蓄沉積物,山脈將以緩慢速度抬升.
4) 北西段發(fā)育“無(wú)根密集”斷層,巖體破碎,構(gòu)造強(qiáng)烈,受構(gòu)造影響基巖將處于抬升階段.支流1、2、3河谷呈表現(xiàn)為“深V”下切作用侵蝕作用,河流兩側(cè)將形成陡峭河谷地形,即進(jìn)一步證明基巖的抬升,區(qū)域呈現(xiàn)出快速隆升狀態(tài).
致謝中國(guó)地質(zhì)科學(xué)院探礦工藝研究所為本文提供野外地質(zhì)災(zāi)害信息,審稿專(zhuān)家及編輯對(duì)本文提出寶貴的修改意見(jiàn),謹(jǐn)致謝意.
[1] 李勇,周榮軍, DENSMORE A L,等. 青藏高原東緣大陸動(dòng)力學(xué)過(guò)程與地質(zhì)響應(yīng)[M]. 北京:地質(zhì)出版社,2006:10-16,106-116,125-129.
[2] OUCHI S. Response of alluvial rivers to slow active tectonic movement[J]. Geological Society of America Bulletin,1985,96(4):504-515.
[3] 李亞林,王成善,王謀,等. 藏北長(zhǎng)江源地區(qū)河流地貌特征及其對(duì)新構(gòu)造運(yùn)動(dòng)的響應(yīng)[J]. 中國(guó)地質(zhì),2006,33(2):374-382.
[4] 顏照坤,李勇,黃潤(rùn)秋,等. 龍門(mén)山北段平通河流域地貌演化過(guò)程[J]. 山地學(xué)報(bào),2012,30(2)：136-146.
[5] 李勇,黃潤(rùn)秋,周榮軍,等. 龍門(mén)山地震帶的地質(zhì)背景與汶川地震的地表破裂[J]. 工程地質(zhì)學(xué)報(bào),2009,17(1):3-16.
[6] 趙國(guó)華,李勇,閆亮,等. 汶川Ms8.0級(jí)地震后龍門(mén)山構(gòu)造地貌及地表侵蝕過(guò)程研究：以湔江海子段河為例[J]. 四川地震,2014(1):15-23.
[7] 陳彥杰,鄭光佑,宋國(guó)城. 面積尺度與空間分布對(duì)流域面積高程積分及其地質(zhì)意義的影響[J]. 地理學(xué)報(bào)(臺(tái)灣),2005,39:53-59.
[8] ROBERT N P, ALEXANDER L D, NICHOLAS J R, et al. Mass wasting triggered by the 2008 Wenchuan earthquake is greater than orogenic growth[J]. Nature Geoscience,2011,4(7):449-452.
[9] 許強(qiáng). 四川省“8·13”特大泥石流災(zāi)害特點(diǎn)、成因與啟示[J]. 工程地質(zhì)學(xué)報(bào),2010,18(5):596-608.
[10] 四川省水利廳. 都江堰水利詞典[M]. 北京:科學(xué)出版社,2004:254-255.
[11] 劉樹(shù)根. 龍門(mén)山?jīng)_斷帶與川西前陸盆地的形成演化[M]. 成都:成都科技大學(xué)出版社,1993.
[12] 江娃利,謝新生. 龍門(mén)山后山斷裂汶川Ms8.0地震地表破裂帶[J]. 國(guó)際地震動(dòng)態(tài),2009(4):4.
[13] 邵崇建,李勇,顏照坤,等. 基于DEM數(shù)據(jù)的龍門(mén)山流域地形起伏度研究[J]. 四川師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,38(5):766-773.
[14] 顏照坤,李勇,黃潤(rùn)秋,等. 汶川地震對(duì)龍門(mén)山地區(qū)平通河河流系統(tǒng)沉積物輸移的影響[J]. 自然雜志,2011,33(6):337-339.
[15] 陳浩,陳朝鎮(zhèn). 汶川地震驅(qū)動(dòng)的滑坡災(zāi)害空間發(fā)育規(guī)律研究:以四川省安縣為例[J]. 四川師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2016,39(3):427-431.
[16] 曾洪揚(yáng),鄧斌,李勇. 汶川大地震對(duì)岷江上游水資源與水環(huán)境的影響[J]. 四川師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,32(1):134-138.
[17] 李勇,周榮軍,ALEXANDER L,等. 映秀—北川斷裂的地表破裂與變形特征[J]. 地質(zhì)學(xué)報(bào),2008,82(12):1688-1706.
[18] 黃潤(rùn)秋. 汶川8.0級(jí)地震觸發(fā)崩滑災(zāi)害機(jī)制及其地質(zhì)力學(xué)模式[J]. 巖石力學(xué)與工程學(xué)報(bào),2009,28(6):1239-1249.
馮金煒, 馬 超, 李 勇, 陳 歡, 邱占林, 劉玉法, 張 威, 趙 重
( 1. 成都理工大學(xué) 油氣藏地質(zhì)及開(kāi)發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室, 四川 成都 610059； 2. 四川省地震局 工程地震研究院, 四川 成都 610041； 3. 中國(guó)地質(zhì)科學(xué)院 探礦工藝研究所, 四川 成都 610031； 4. 龍巖學(xué)院 資源工程學(xué)院, 福建 龍巖 364012)
地貌是不斷發(fā)展變化的,地貌發(fā)展變化的物質(zhì)過(guò)程稱(chēng)地貌過(guò)程,內(nèi)力和外力是塑造地貌的2種營(yíng)力,地貌是內(nèi)力過(guò)程與外力過(guò)程相互統(tǒng)一的產(chǎn)物.龍門(mén)山是青藏高原山脈中落差高度變化最大的山脈,構(gòu)造作用強(qiáng)烈,李勇等[1]通過(guò)研究龍門(mén)山構(gòu)造的地貌標(biāo)志定量計(jì)算了龍門(mén)山主斷裂的逆沖速率與走滑速率的關(guān)系,進(jìn)一步分析了龍門(mén)山周緣地貌演化趨勢(shì).Ouchi[2]認(rèn)為斷層及褶皺影響河流形態(tài)及搬運(yùn)沉積過(guò)程,而河流的多樣性則影響著周緣的地貌的演化.李亞林等[3]通過(guò)研究構(gòu)造運(yùn)動(dòng)對(duì)河流地貌特征的響應(yīng)指出：早期逆滑和走滑斷裂控制水系形態(tài),而晚期的正斷層和地塹構(gòu)造控制河流地貌的發(fā)育演變.2008年“5·12”地震導(dǎo)致映秀—北川斷裂發(fā)生地表破裂,導(dǎo)致龍門(mén)山地區(qū)基巖的同震隆升變形[1,3],觸發(fā)大規(guī)模的同震滑坡、崩塌等地質(zhì)災(zāi)害發(fā)生,形成大面積松散物質(zhì)堆積,在氣候的作用下進(jìn)一步改造龍門(mén)山區(qū)域的地貌[4-6].因此研究區(qū)域性的構(gòu)造抬升及河流侵蝕對(duì)未來(lái)地貌演化方向的影響十分重要.本文基于ArcGis技術(shù)提取龍門(mén)山北段通口河治城段的地形因子及河流剖面等參數(shù),進(jìn)一步探討構(gòu)造、氣候等因素對(duì)治城段地貌演化發(fā)展趨勢(shì)的影響.
1 資料與處理方法
通口河系屬龍門(mén)山中段山前河流,是涪江的重要支流之一,起源于九頂山南麓,西北—南東流向,在彭州市境內(nèi)穿越龍門(mén)山,在丹景山鎮(zhèn)進(jìn)入成都平原,是涪江流域的一級(jí)支流.本文研究區(qū)域通口河段挾持于后龍門(mén)山深斷裂與江油—灌縣大斷裂之間,最高海拔2 260 m,最低海拔600 m(圖1),落差1 660 m,地勢(shì)呈現(xiàn)自北西往南東降低(圖2).
本次研究通過(guò)無(wú)人機(jī)勘察及野外實(shí)測(cè)獲得高程點(diǎn),由高程點(diǎn)繪制成等高線(xiàn).利用ArcGis 10.2軟件將等高線(xiàn)矢量數(shù)據(jù)生成TIN,再由TIN轉(zhuǎn)成DEM,經(jīng)過(guò)DEM柵格數(shù)據(jù)多次處理獲得通口河流域治城段地形特征和水系網(wǎng)絡(luò)特征.
1.1地形坡度的獲取利用ArcGis 10.2軟件,將獲得的高分辨率DEM數(shù)據(jù)(1∶1萬(wàn),柵格大小5 m),經(jīng)過(guò)ArcToolbox中3D Analyst工具的柵格表面計(jì)算地面坡度的平均值,獲得治城段平均坡度分布區(qū)間為0°～69.1°,通過(guò)圖層屬性將平均坡度范圍分為4個(gè)等級(jí),分別為：0°～15.6°、15.6°～30.9°、30.9°～44.5°、44.5°～69.1°.4個(gè)坡度區(qū)間的分布范圍如圖2所示.
1.2水系網(wǎng)絡(luò)的獲取利用ArcGis 10.2軟件,將得到的DEM進(jìn)行填洼獲得無(wú)洼DEM,經(jīng)過(guò)水文分析工具得到流向和流量,將閾值設(shè)置為5 000獲得流水累計(jì)柵格,以此得到通口河流域水系分布圖(圖3).通口河流域主要發(fā)育4個(gè)較大支流：小園河支流(支流1)、白草河支流(支流2)、白坭溝支流(支流3)、玉溪河支流(支流4).下面將通過(guò)對(duì)通口河支流的河流縱剖面和橫剖面特征進(jìn)一步分析.
圖 1 通口河流域縱剖面圖
圖 2 治城段流域平均坡度分布圖
圖 3 治城段流域地形與水系圖Fig. 3 Topographic and water system of Zhicheng river basin
2 通口河治城段流域發(fā)育特征
2.1通口河治城段流域的地形特征通過(guò)DEM計(jì)算出來(lái)的通口河流域治城段平均坡度分布圖(圖2)可宏觀(guān)表達(dá)該區(qū)域地貌特征.北川老縣城附近(即距北川—映秀斷裂北西側(cè)4 km區(qū)域內(nèi))平均坡度較大,多介于30.9°～62.3°；玉溪河和通口河河流區(qū)域坡度較小,多為10°～30°；研究區(qū)北西段(以禹里鄉(xiāng)和漩坪鄉(xiāng)連線(xiàn)為界,下同)劉家窩—茶園坪—白巖坪—燈草坪斷斷裂破碎帶區(qū)域平均坡度較大,多介于30.9°～62.3°；挾于各斷裂帶之間的區(qū)域區(qū)域均屬于小角度坡度：16.3°～26.0°.
不同地層巖性(圖4)可體現(xiàn)出不同坡度分布,結(jié)合地質(zhì)圖與坡度圖可進(jìn)一步分析通口河流域的地貌特征.北川映秀斷裂北西向(4 km)為志寒武系油房組,巖性以中、粗粒巖屑砂巖、粉砂巖互層,南東側(cè)以泥盆系和石炭系灰?guī)r為主,北西向一側(cè)坡度明顯較大.研究區(qū)中部以志留系茂縣群第二組灰?guī)r環(huán)繞茂縣群第一組千枚巖,灰?guī)r區(qū)域坡度大于千枚巖區(qū)域.位于研究區(qū)北段劉家窩—茶園坪—白巖坪—燈草坪斷層附近巖性以千枚巖夾灰?guī)r、砂巖為主,坡度大于單純千枚巖地層.
通過(guò)斷裂帶與巖性對(duì)坡度分布影響可知,斷裂(劉家窩—茶園坪—白巖坪—燈草坪斷裂、北川斷裂)上下盤(pán)周?chē)钠露染哂幸欢ǖ牟町愋?斷層上盤(pán)4～6 km范圍內(nèi)坡度大,說(shuō)明2條逆沖斷裂上盤(pán)處于隆升狀態(tài),斷裂上盤(pán)的河流侵蝕下切作用將相對(duì)較快.通過(guò)對(duì)比地形坡度(圖2)及地層巖性(圖4)可知巖性的“軟硬”程度同樣對(duì)坡度的分布起一定的作用：研究表明砂巖和灰?guī)r區(qū)域坡度明顯大于千枚巖區(qū)域坡度.
2.2通口河治城段流域的水系特征本文基于ArcGis 10.2軟件,將得到的DEM進(jìn)行填洼或得無(wú)洼DEM,經(jīng)過(guò)水文分析工具得到流向和流量,將閾值設(shè)置為5 000獲得流水累計(jì)柵格,以此得到通口河流域水系分布圖(圖2).通口河流域主要發(fā)育4個(gè)較大支流：小園河支流(支流1)、白草河支流(支流2)、白坭溝支流(支流3)、玉溪河支流(支流4).下面通過(guò)對(duì)通口河支流的河流縱剖面和橫剖面特征進(jìn)一步分析.
圖4治城段流域地質(zhì)構(gòu)造簡(jiǎn)圖
Fig.4ThesimplifiedgeologicalandstructuralmapofZhicehngriverbasin
2.2.1河流縱剖面 受到地質(zhì)構(gòu)造、流域基部巖性特性、氣候降雨等諸多因素影響,河流縱剖面呈現(xiàn)不同形態(tài).研究區(qū)域較小,氣候條件幾乎相同；基部巖性以志留系茂縣群千枚巖為主,4條支流流向大致垂直于巖層走向,因而該區(qū)域河流剖面形態(tài)主受構(gòu)造影響.通過(guò)對(duì)河流縱剖面形態(tài)特征的分析可有效反映通口河流域的構(gòu)造差異性.
關(guān)于河流剖面與地貌演化的關(guān)系前人做出了大量的研究,陳彥杰等[7]通過(guò)對(duì)河流縱剖面模型進(jìn)行了總結(jié)：1)地形面在受到地體抬升后構(gòu)造運(yùn)動(dòng)即停止,則初始形成的河流其縱剖面的下凹程度較小或近乎直線(xiàn), 可以用線(xiàn)性剖面來(lái)擬合；2)隨著中上游侵蝕作用的進(jìn)行,河流將侵蝕下來(lái)的物質(zhì)搬運(yùn)至下游堆積,使得河流縱剖面的下凹程度逐漸變大而逐漸演化成指數(shù)剖面；3)隨著下游物質(zhì)的堆積,河流中-上游的侵蝕,河流縱剖面下凹程度進(jìn)一步加大,可以用對(duì)數(shù)函數(shù)來(lái)擬合河流縱剖面;4)再進(jìn)一步發(fā)展,如果流量增加或由于構(gòu)造運(yùn)動(dòng)等形成破碎帶導(dǎo)致輸沙量增加,使得河流縱剖面的下凹程度快速增大而呈現(xiàn)出乘冪函數(shù)剖面.河流縱剖面的演化趨勢(shì)為:線(xiàn)性→指數(shù)→對(duì)數(shù)→乘冪.
通口河流域治城段圖幅較小,在考慮到流域內(nèi)各條支流的氣候和巖性類(lèi)似的情況下,構(gòu)造運(yùn)動(dòng)是影響河流縱剖面的主控因素.通過(guò)對(duì)研究區(qū)4條支流(圖2、圖5)對(duì)比分析：支流1、支流2、支流3下凹程度低,呈近線(xiàn)性指數(shù)擬合,即考慮該3條支流較“年輕”,在氣候和基底巖性相似的情況下,支流1靠近劉家窩—茶園坪—白巖坪—燈草坪斷層；支流2穿過(guò)劉家窩—茶園坪—白巖坪—燈草坪斷層、盛家院—密地坪—楊柳松斷層及茂汶斷裂；支流3位于盛家院—密地坪—楊柳松斷層附近；很可能由于構(gòu)造導(dǎo)致山體相對(duì)抬升,使得這3條支流的縱剖面接近線(xiàn)性擬合.支流4和通口河主干接近對(duì)數(shù)剖面,說(shuō)明這2條河流演化速度慢受地質(zhì)構(gòu)造作用相對(duì)較少,河流進(jìn)入“中、老年化”,根據(jù)地質(zhì)資料及實(shí)地路線(xiàn)觀(guān)察,該2條河流附近斷層地質(zhì)作用較弱,處于相對(duì)穩(wěn)定狀態(tài).
圖 5 治城流域主要支流縱剖面
2.2.2河流橫剖面 通過(guò)研究區(qū)橫剖面線(xiàn)可得各支流情況(圖6)：支流1、支流2、支流3河谷呈深“V”型,說(shuō)明河流正進(jìn)行快速下切侵蝕作用,也進(jìn)一步反映這3條支流范圍內(nèi)抬升作用較強(qiáng).支流3和通口河(湔江)主干河谷為U型谷,說(shuō)明這2條河流有側(cè)向侵蝕作用.根據(jù)河流發(fā)育特征,河流先以縱向下切作用為主,當(dāng)河流縱剖面達(dá)到平衡狀態(tài)時(shí),河流再以側(cè)向侵蝕為主.通過(guò)5條河谷對(duì)比分析：支流4和通口河主流相對(duì)于支流1、2、3下切侵蝕速率慢,構(gòu)造相對(duì)穩(wěn)定,山體抬升不明顯,地貌不會(huì)發(fā)生太劇烈的演化.
圖 6 流域治城段河流橫剖面分析圖
3 通口河治城段流域地貌發(fā)育的主要影響因素
3.1構(gòu)造活動(dòng)對(duì)地貌演化的影響2008年汶川地震的發(fā)生導(dǎo)致震區(qū)地質(zhì)災(zāi)害大面積爆發(fā),據(jù)野外核查,研究區(qū)構(gòu)造活動(dòng)的較明顯(斷層附近)地區(qū)地質(zhì)災(zāi)害成面狀分布.在地質(zhì)構(gòu)造的發(fā)育的地區(qū),構(gòu)造面發(fā)育導(dǎo)致周?chē)鷰r體破碎,穩(wěn)定性差,極易發(fā)生地質(zhì)災(zāi)害.構(gòu)造對(duì)地質(zhì)災(zāi)害的控制主要表現(xiàn)在以下幾方面：1) 斷層均以北東向排布,控制著研究區(qū)內(nèi)山脈及河流的展布.區(qū)域內(nèi)多條小斷層交匯,交匯處形成大面積“X”節(jié)理,使得巖體失穩(wěn),形成規(guī)模較大的地質(zhì)災(zāi)害.斷裂帶附近巖體較其它地區(qū)破碎,研究區(qū)西北段無(wú)根密集斷層、中段茂汶斷裂、南西段北川斷裂比較典型.2) 現(xiàn)代新構(gòu)造運(yùn)動(dòng),特別是地震活動(dòng),直接破壞了巖土體的結(jié)構(gòu),使巖土體的物理力學(xué)強(qiáng)度降低.研究區(qū)主要以絹云千枚巖為主,巖性較為軟弱,在斷層作用下,節(jié)理面發(fā)育,巖體破碎,當(dāng)?shù)乇硭疂B入時(shí)使巖層軟化形成活動(dòng)帶,當(dāng)遇大量降雨滲入土體,增加了土體自重或下滑力,使土體的粘聚力和磨擦力減小,便產(chǎn)生了滑坡、崩塌,也為泥石流提供大量物源.3) 構(gòu)造帶中多期次構(gòu)造運(yùn)動(dòng)的疊加,使巖體完整性降低,更易于被風(fēng)化,同時(shí)由于破碎帶的影響,其風(fēng)化深度更大,研究區(qū)多條區(qū)域性斷裂為后期在北川斷層上的疊加.4) 構(gòu)造裂隙及斷層面常與卸荷裂隙組成向臨空面的不利組合.
2008年特大地震后不僅造成龍門(mén)山斷裂帶急劇的隆升,也進(jìn)一步觸發(fā)了大面積的地質(zhì)災(zāi)害的爆發(fā).根據(jù)2016年9月野外勘察,共查明研究區(qū)域內(nèi)崩塌12處、滑坡56處、泥石流5處,地質(zhì)災(zāi)害面積達(dá)到5.6 km2,松散堆積物體積約0.019 km3.2010年,對(duì)龍門(mén)山北川—映秀斷裂帶同震巖石隆升研究得出的數(shù)據(jù)中,該研究區(qū)域增加的山脈體積估算值約為0.176 km3,大于松散堆積物體積.Robert等[8]也指出,沿?cái)鄬臃较蛏系牟糠謪^(qū)域同震巖石隆升增加的山脈體積大于地震導(dǎo)致的堆積物體積,即山脈體積增加.研究區(qū)北西段發(fā)育小型密集“無(wú)根斷層”,地震后巖石隆升作用較為劇烈(圖7),雖然同震地質(zhì)災(zāi)害產(chǎn)生了松散堆積體,但研究區(qū)基巖隆升體積大于同震松散物質(zhì)體積,即研究區(qū)山體體積并未“虧損”,而是處于抬升狀態(tài).
圖 7 汶川地震觸發(fā)的同震抬升圖(據(jù)文獻(xiàn)[8]修改)
3.2氣候?qū)Φ孛惭莼挠绊懲诤恿饔蛑纬嵌螝夂驅(qū)儆趤啛釒貪?rùn)疾風(fēng)氣候,降雨量充足,1967年為最大年降雨量,達(dá)到2 340 mm,日最大降雨量達(dá)到101 mm,年平均降雨量可達(dá)1 399.1 mm.強(qiáng)降雨時(shí)期主要集中在7、8、9月份(圖8),達(dá)到全年降雨量的70%以上.2008年汶川后大面積地質(zhì)災(zāi)害爆發(fā),2013年“7.9”洪災(zāi)使得一些潛在地質(zhì)災(zāi)害發(fā)生,暴雨誘發(fā)的泥石流、崩塌、滑坡共31處.暴雨導(dǎo)致的山洪峰值較大,攜帶溝谷內(nèi)崩塌和滑坡物質(zhì)進(jìn)入河流.研究區(qū)域較大的泥石流有倒開(kāi)門(mén)泥石流和羊盤(pán)溝泥石流,每年雨季均攜帶大量松散物質(zhì)分別匯入通口河和白草河,地震導(dǎo)致的崩塌物質(zhì)是主要的物源.許強(qiáng)通過(guò)實(shí)驗(yàn),表明：地震災(zāi)區(qū)發(fā)生大型泥石流的臨界雨降為35～40 mm/h,最低僅為15 mm/h[9].
根據(jù)實(shí)地野外核查,泥石流位置均靠近河流,由于構(gòu)造(主要地震)產(chǎn)生的大面積松散物質(zhì)通過(guò)河流作用搬運(yùn)至下游地區(qū).由于山脈物質(zhì)的遷移,影響地貌的改變.氣候(以降水為主)影響地貌演化的最直接表現(xiàn)便是河流中輸砂率的變化.
研究區(qū)松散物通過(guò)洪流由白草河、白泥溝、玉溪河等支流攜帶搬運(yùn)至通口河,根據(jù)通口河下游將軍水文站近十年數(shù)據(jù)表明,降雨量直接影響輸砂量,7～9月份為研究區(qū)雨季,輸砂量野達(dá)到峰值.圖8體現(xiàn)出降雨量和輸砂率呈正相關(guān)關(guān)系.
圖 8 1-12月通口河流域平均輸砂率與北川縣的月降水量對(duì)比分析
前人通過(guò)流域輸砂量的計(jì)算獲來(lái)得剝蝕速率,本文通過(guò)輸砂率估算出剝蝕速率,以此研究流域內(nèi)地貌的演化發(fā)展.據(jù)綿陽(yáng)將軍石水文站實(shí)測(cè)懸砂數(shù)據(jù)統(tǒng)計(jì)[10]：多年平均輸砂率171 kg/s,多年平均輸砂量1.8 kg/m3,多年平均汛期含砂量為2.25 kg/m3,最大輸砂量為320萬(wàn)t(1964年),最小年輸砂量96萬(wàn)t(1974年).據(jù)此比擬推求通口河流域年平均輸砂量為2.27×106t,沉積物密度取1.50 t/m3,即1.51×106m3.通口河流域面積為2 808 km2,計(jì)算出通口河流域剝蝕速率為0.537 mm/a.侵蝕速率應(yīng)包括風(fēng)力搬運(yùn)作用、化學(xué)作用等,因此通口河流域的剝蝕速率應(yīng)遠(yuǎn)大于0.537 mm/a.研究數(shù)據(jù)表明該區(qū)域平均隆升幅度為0.456 m,即在未來(lái)850年內(nèi)將剝蝕掉地震成因的山脈增加量. 在地質(zhì)構(gòu)造不太強(qiáng)雷的情況下,在氣候因素的影響下,山脈物質(zhì)增加量將被河流遷移至下游.
前人對(duì)龍門(mén)山在非地震期間的隆升速率已經(jīng)進(jìn)行了大量研究.劉樹(shù)根[11]根據(jù)測(cè)量地表的變形機(jī)制獲得龍門(mén)山逆沖推覆帶的年均隆升速率在0.3～0.4 mm之間；江娃利等[12]研究指出龍門(mén)山構(gòu)造帶可達(dá)0.2～2.5 m的垂直位移量；邵崇建等[13]等根據(jù)多年(1970—1997年)野外實(shí)測(cè)資料得出在理縣附近的隆升速率最大,可達(dá)到 3. 5 mm/a.龍門(mén)山處于非地震階段的垂直隆升速率目前數(shù)據(jù)值較多,大致處于0.2～3.5 mm/a范圍內(nèi).對(duì)比通口河的侵蝕速率(>0.537 mm/a)與隆升速率(0.2～3.5 mm/a),區(qū)域總體呈現(xiàn)出抬升作用,海拔處于上升階段.
4 地貌演化分區(qū)
綜合上述坡度、河網(wǎng)、構(gòu)造、氣候等要素對(duì)地貌演化的影響進(jìn)行分析,對(duì)通口河流域治城段劃分為兩個(gè)區(qū)域,以此討論地貌演化過(guò)程.
4.1Ⅰ區(qū)(緩慢抬升區(qū)) Ⅰ區(qū)位于北川斷裂上盤(pán)附近,汶川地震使得該區(qū)山脈抬升.通過(guò)各支流縱坡面(圖5)可知位于該處的通口河段具有破折點(diǎn)和匯入通口河的玉溪河(支流4)具有明顯的破折點(diǎn),說(shuō)明該處正處于抬升狀態(tài),并且正經(jīng)歷河流下切作用,但由河流橫剖面(圖6)反應(yīng)出這2條河流河谷較為寬闊,即河流下切作用弱.此外,該處Ⅰ區(qū)巖性由志留系千枚巖過(guò)渡到較為堅(jiān)硬的砂巖,進(jìn)一步制約了河流剝蝕速率.Ⅰ區(qū)堰塞湖位置也可能淤積上流搬運(yùn)下來(lái)的河流物質(zhì),使得松散物質(zhì)堆積量增加.Ⅰ區(qū)剝蝕速率小于基巖抬升速率,可以推斷Ⅰ區(qū)可能為緩慢隆升的地貌演化趨勢(shì).
4.2Ⅱ區(qū)(快速抬升區(qū)) Ⅱ區(qū)為研究區(qū)北西段,該處發(fā)育大量“無(wú)根密集”斷層,巖性較為破碎,表現(xiàn)為構(gòu)造作用強(qiáng)烈.支流1、2、3縱坡面呈深“V”型,縱剖面為線(xiàn)性指數(shù),即該處河流下切速度快,進(jìn)一步說(shuō)明該區(qū)屬于隆升階段.此外,由地形坡度(圖3)可知,位于支流1、2、3附近坡度大,河流兩側(cè)的將形成陡峭的河谷地形.雖然在汶川地震觸發(fā)下形成大面積的崩塌、滑坡等松散堆積物質(zhì),在河流作用下將被搬運(yùn)到下流,但通過(guò)對(duì)比隆升幅度及侵蝕速率,山體并未“虧損”.Ⅱ區(qū)基巖抬升速率大于河流侵蝕作用,且在較為劇烈構(gòu)造運(yùn)動(dòng)的情況下,地貌演化將呈快速抬升趨勢(shì).
圖 9 治城段流域地貌演化分區(qū)
5 結(jié)束語(yǔ)
通過(guò)對(duì)龍門(mén)山北段通口河流域地貌、構(gòu)造、氣候等特征進(jìn)行分析,初步得下以下幾點(diǎn)認(rèn)識(shí).
1) 斷層上盤(pán)附近坡度普遍較大,區(qū)域處于抬升狀態(tài),河流侵蝕下切作用強(qiáng),使得地形變得陡峭.巖性影響地形特征,砂巖和灰?guī)r區(qū)域地形坡度明顯大于千枚巖區(qū)域坡度.
2) 以如今通口河流域治城段的構(gòu)造作用和氣候特征,大部分河流處于“年輕”狀態(tài),下切剝蝕速率較快,搬運(yùn)能力較強(qiáng),上游沉積物極易搬運(yùn)到下游,使得流域局部地貌往削低趨勢(shì)發(fā)展.
3) 南東段(北川斷裂上盤(pán)附近)受構(gòu)造作用影響處于抬升階段,該位置河流河谷較為寬闊,剝蝕山體速率較慢,受堰塞湖等因素影響將繼續(xù)積蓄沉積物,山脈將以緩慢速度抬升.
4) 北西段發(fā)育“無(wú)根密集”斷層,巖體破碎,構(gòu)造強(qiáng)烈,受構(gòu)造影響基巖將處于抬升階段.支流1、2、3河谷呈表現(xiàn)為“深V”下切作用侵蝕作用,河流兩側(cè)將形成陡峭河谷地形,即進(jìn)一步證明基巖的抬升,區(qū)域呈現(xiàn)出快速隆升狀態(tài).
致謝中國(guó)地質(zhì)科學(xué)院探礦工藝研究所為本文提供野外地質(zhì)災(zāi)害信息,審稿專(zhuān)家及編輯對(duì)本文提出寶貴的修改意見(jiàn),謹(jǐn)致謝意.
[1] 李勇,周榮軍, DENSMORE A L,等. 青藏高原東緣大陸動(dòng)力學(xué)過(guò)程與地質(zhì)響應(yīng)[M]. 北京:地質(zhì)出版社,2006:10-16,106-116,125-129.
[2] OUCHI S. Response of alluvial rivers to slow active tectonic movement[J]. Geological Society of America Bulletin,1985,96(4):504-515.
[3] 李亞林,王成善,王謀,等. 藏北長(zhǎng)江源地區(qū)河流地貌特征及其對(duì)新構(gòu)造運(yùn)動(dòng)的響應(yīng)[J]. 中國(guó)地質(zhì),2006,33(2):374-382.
[4] 顏照坤,李勇,黃潤(rùn)秋,等. 龍門(mén)山北段平通河流域地貌演化過(guò)程[J]. 山地學(xué)報(bào),2012,30(2)：136-146.
[5] 李勇,黃潤(rùn)秋,周榮軍,等. 龍門(mén)山地震帶的地質(zhì)背景與汶川地震的地表破裂[J]. 工程地質(zhì)學(xué)報(bào),2009,17(1):3-16.
[6] 趙國(guó)華,李勇,閆亮,等. 汶川Ms8.0級(jí)地震后龍門(mén)山構(gòu)造地貌及地表侵蝕過(guò)程研究：以湔江海子段河為例[J]. 四川地震,2014(1):15-23.
[7] 陳彥杰,鄭光佑,宋國(guó)城. 面積尺度與空間分布對(duì)流域面積高程積分及其地質(zhì)意義的影響[J]. 地理學(xué)報(bào)(臺(tái)灣),2005,39:53-59.
[8] ROBERT N P, ALEXANDER L D, NICHOLAS J R, et al. Mass wasting triggered by the 2008 Wenchuan earthquake is greater than orogenic growth[J]. Nature Geoscience,2011,4(7):449-452.
[9] 許強(qiáng). 四川省“8·13”特大泥石流災(zāi)害特點(diǎn)、成因與啟示[J]. 工程地質(zhì)學(xué)報(bào),2010,18(5):596-608.
[10] 四川省水利廳. 都江堰水利詞典[M]. 北京:科學(xué)出版社,2004:254-255.
[11] 劉樹(shù)根. 龍門(mén)山?jīng)_斷帶與川西前陸盆地的形成演化[M]. 成都:成都科技大學(xué)出版社,1993.
[12] 江娃利,謝新生. 龍門(mén)山后山斷裂汶川Ms8.0地震地表破裂帶[J]. 國(guó)際地震動(dòng)態(tài),2009(4):4.
[13] 邵崇建,李勇,顏照坤,等. 基于DEM數(shù)據(jù)的龍門(mén)山流域地形起伏度研究[J]. 四川師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,38(5):766-773.
[14] 顏照坤,李勇,黃潤(rùn)秋,等. 汶川地震對(duì)龍門(mén)山地區(qū)平通河河流系統(tǒng)沉積物輸移的影響[J]. 自然雜志,2011,33(6):337-339.
[15] 陳浩,陳朝鎮(zhèn). 汶川地震驅(qū)動(dòng)的滑坡災(zāi)害空間發(fā)育規(guī)律研究:以四川省安縣為例[J]. 四川師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2016,39(3):427-431.
[16] 曾洪揚(yáng),鄧斌,李勇. 汶川大地震對(duì)岷江上游水資源與水環(huán)境的影響[J]. 四川師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,32(1):134-138.
[17] 李勇,周榮軍,ALEXANDER L,等. 映秀—北川斷裂的地表破裂與變形特征[J]. 地質(zhì)學(xué)報(bào),2008,82(12):1688-1706.
[18] 黃潤(rùn)秋. 汶川8.0級(jí)地震觸發(fā)崩滑災(zāi)害機(jī)制及其地質(zhì)力學(xué)模式[J]. 巖石力學(xué)與工程學(xué)報(bào),2009,28(6):1239-1249.
