趙 巖,孟興民,2,慶 豐,劉林通,郭富赟,3
(1.蘭州大學(xué) 資源環(huán)境學(xué)院 西部環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,甘肅 蘭州 730000;2.甘肅省環(huán)境地質(zhì)與災(zāi)害防治工程研究中心,甘肅 蘭州 730000;3.甘肅省地質(zhì)環(huán)境監(jiān)測院,甘肅 蘭州 730050)
泥石流發(fā)源于地貌學(xué)[1],說其特殊因?yàn)樗且环N主要由水和重力作用引發(fā)的強(qiáng)烈的地表流水侵蝕、搬運(yùn)與堆積的過程[2];而說其不特殊是因?yàn)槟嗍骰顒?dòng)屬于地貌演化過程中眾多夷平作用過程的其中一種。既然泥石流活動(dòng)屬于地貌演化的一種形式,而地貌演化是分階段的,所以泥石流溝的發(fā)育階段也可以通過地貌參數(shù)來定量劃分。本文通過分析與泥石流發(fā)育階段相關(guān)的地貌參數(shù),將泥石流溝的發(fā)育階段進(jìn)行定量劃分。
地貌信息熵(H)與超熵(?xP)被應(yīng)用到了許多領(lǐng)域的研究中。劉小鳳等[3]將H應(yīng)用在了構(gòu)造活動(dòng)和區(qū)域穩(wěn)定性評價(jià);鐵永波等[4]將H應(yīng)用在泥石流溝谷的危險(xiǎn)性評價(jià)研究中;孫然好等[5]將H應(yīng)用在山體演化階段的判定;王曉朋等[6]利用H來判斷泥石流的發(fā)育程度,用于定量評價(jià)泥石流危險(xiǎn)性;李雅輝等[7]基于H對泥石流敏感性進(jìn)行定量評價(jià);岳俊生等[8]將?xP應(yīng)用在泥石流發(fā)育現(xiàn)狀及其危險(xiǎn)度評價(jià);鄒強(qiáng)等[9]將?xP應(yīng)用在泥石流活躍度的定量分析;王鈞等[10]將H應(yīng)用在震后泥石流危險(xiǎn)性評價(jià);劉麗娜等[11]利用H進(jìn)行泥石流危險(xiǎn)性評價(jià);張寶軍等[12]將H應(yīng)用在沖溝溝頭活躍度的判定。H和?xP的廣泛應(yīng)用,有必要對其來源及其與泥石流發(fā)育階段的關(guān)系進(jìn)行分析與討論。
1899年,Davis[13]提出地貌循環(huán)理論,并將地貌演化階段分為幼年期、壯年期和老年期三個(gè)階段,1952年Strahler[14]提出侵蝕流域的面積高程分析方法,并用曲線積分值S的大小來劃分侵蝕流域地貌演化階段:S>0.6時(shí)為幼年期,0.35≤S≤0.6時(shí)為壯年期,S<0.35是為老年期,從而將戴維斯的地貌發(fā)育模型定量化。面積高程積分值S的計(jì)算方法如下。
設(shè)A表示流域面積,a表示該流域內(nèi)某條等高線上的面積,h表示該等高線與流域內(nèi)最低點(diǎn)的高差,ΔH表示流域內(nèi)最大相對高差,則對直角坐標(biāo)軸上一系列X=a/A,Y=h/ΔH組成的點(diǎn)進(jìn)行擬合,得到擬合效果最好的函數(shù)y=f(x)即為面積高程積分曲線,而其積分即曲線包圍的面積,稱為面積高程積分值:
(1)
由于地貌系統(tǒng)演化過程的數(shù)學(xué)描述與熱傳導(dǎo)方程相似,1962年里奧普和拜里把熵的概念引入地貌學(xué)中。1970年Leopold與Langbein也提出過地貌學(xué)熵的定義,熱力場是由溫度T(x,y)和熱量Q來表現(xiàn),地貌場類似可以用高度h(x,y)和質(zhì)量M表現(xiàn),進(jìn)而提出地貌學(xué)熵的定義:
ds=dQ/T?dM/h。
(2)
1987年我國學(xué)者艾南山[15]以Davis地貌循環(huán)理論、Penck的內(nèi)外營力同時(shí)作用的地貌發(fā)育模型以及Shannon的信息熵理論,將熵的概念和計(jì)算方法引入到Strahler曲線和積分中,提出了侵蝕流域系統(tǒng)的地貌信息熵理論并推導(dǎo)出其值H的計(jì)算公式,推導(dǎo)過程簡述如下。
1948年Shannon[16]提出能夠定量描述系統(tǒng)演化階段的的信息熵概念:
(3)
后有人將其推廣到連續(xù)情形(p(x)為密度函數(shù)):
(4)
用面積高程積分曲線f(x)和面積高程積分值S構(gòu)造密度函數(shù):
(5)
將(1)式、(5)式帶入(4)式便可得到地貌信息熵:
H=S-1-lnS。
(6)
在之后的研究中又將信息熵由侵蝕流域引入到一般流域,利用Иванов的河流縱剖面方程(式中,h為縱剖面上某點(diǎn)與河口的高差,l為縱剖面上某點(diǎn)距河口的距離,H和L分別為河源與河口之間的高差和水平距離,N為形態(tài)指數(shù),且0
h=H(l/L)N。
(7)
使用同樣的方法推導(dǎo)出利用縱剖面形態(tài)指數(shù)N表示的一般流域的地貌信息熵公式[17]:
H=ln (1+N)-N/(1+N)。
(8)
在1989年為使其可以應(yīng)用在一般的遠(yuǎn)離平衡態(tài)的耗散結(jié)構(gòu)的流域系統(tǒng)內(nèi),引入了超熵的概念,并推導(dǎo)出公式(α、β為參數(shù))[18]:
?xP=[-βα3(α2-1)]/[α(-β-1)-3]。
(9)
1992年蔣忠信[19]提出以溝谷縱剖面形態(tài)指數(shù)N表達(dá)超熵?cái)?shù)學(xué)式,并試用于泥石流溝的判別與活動(dòng)性評價(jià)。
從地貌信息熵和超熵的理論來源來看,其理論基礎(chǔ)分別是Davis的侵蝕循環(huán)理論和Иванов的河流縱剖面方程,對應(yīng)的定量化公式分別為(1)式和(7)式來計(jì)算面積高程積分值(以下簡稱S)和縱剖面形態(tài)指數(shù)(以下簡稱N)??梢钥闯鯯和N分別從侵蝕流域的演化和河谷縱剖面的演化兩種不同的角度對地貌的演化過程進(jìn)行定量描述。
而拋開推導(dǎo)過程中所包含的意義,由S(或N)到H、?xP的推導(dǎo)中關(guān)鍵的一步是通過構(gòu)建含有S(或N)的密度公式(公式5)帶入到熵的表達(dá)式(公式4)中實(shí)現(xiàn)的,這種變換從數(shù)學(xué)角度來看只是對S和N做了一定的拉伸和轉(zhuǎn)換使其表達(dá)方式符合了熵的表現(xiàn)形式(即將原本隨地貌侵蝕演化過程遞減的函數(shù)S轉(zhuǎn)換為隨地貌侵蝕演化過程遞增的函數(shù)H),而H、?xP對地貌發(fā)育階段的劃分完全依據(jù)S和N對應(yīng)公式計(jì)算出的H、?xP值來劃分的。所以地貌熵和超熵的公式只是對S和N做的一種變換,而之后的很多研究只關(guān)注到地貌熵和超熵的應(yīng)用,卻忽視了其來源S和N的應(yīng)用。所以,對于泥石流溝發(fā)育階段的研究,從應(yīng)用的角度來看直接使用S和N值更為簡便。
S、N、H和?xP被廣泛應(yīng)用的根本依據(jù)是這些參數(shù)能夠反映泥石流溝的發(fā)育階段,且不同階段泥石流的活躍程度不同。而利用這些參數(shù)對泥石流溝的發(fā)育階段進(jìn)行定量的劃分方面的研究還比較少。劉希林[21]討論了泥石流溝谷與其橫、縱剖面的演化關(guān)系;蔣忠信[19]通過?xP的三個(gè)特征值將泥石流的發(fā)育階段劃分為五個(gè)階段(表1);蔣忠信[22]研究了泥石流溝谷縱剖面的形態(tài)特征與泥石流活動(dòng)性的關(guān)系;呂學(xué)軍等[23]根據(jù)面積高程積分值將泥石流發(fā)育階段劃分為五個(gè)階段;李泳等[24]用比降曲線和面積高程曲線討論了流域可能發(fā)生泥石流的曲線形態(tài)和演化趨勢。倪華勇[25]采用面積高程分析法對泥石流流域地貌所處的發(fā)育階段進(jìn)行了判別。然而大多研究對泥石流溝發(fā)育階段的判定都是根據(jù)Strahler的三個(gè)地貌階段劃分的。另外,蔣忠信對N的劃分是根據(jù)數(shù)學(xué)意義上的特征點(diǎn)來確定的界限值,呂學(xué)軍等對S的劃分也是根據(jù)經(jīng)驗(yàn)來確定的(表1)。
表1 各學(xué)者對泥石流發(fā)育階段的定量劃分的界限值
首先討論面積高程積分值S,由式(1)可知S實(shí)際上代表著坐標(biāo)軸上一些列x=a/A,y=h/ΔH點(diǎn)的擬合曲線的積分值,這種利用定義計(jì)算S的方法稱為積分曲線法。由積分定義可知,將x分為n等份則(其中V剩表示流域現(xiàn)有體積,V總表示未侵蝕時(shí)流域的總體積(圖1)):
(12)
圖1 S、N值計(jì)算各參數(shù)示意圖
故S值的實(shí)際意義表示該流域經(jīng)侵蝕后所剩體積占未經(jīng)侵蝕時(shí)流域總體積的百分比,這種計(jì)算方法稱為體積比法。Pike等[26]推算出用起伏比法來計(jì)算S(其中Hmean、Hmax和Hmin分別代表流域內(nèi)平均高程、最大高程和最小高程):
(13)
由于平均高程與最小高程之差表示流域的平均高差,最大高差與最小高差之差表示流域最大高差,故由式(11)也可以推出式(13)(其中A為流域面積):
(14)
常直楊等[27]對體積比法和起伏比法計(jì)算方法的準(zhǔn)確性也進(jìn)行了驗(yàn)證,結(jié)果為三種方法計(jì)算出的S值幾乎一致。
對于縱剖面形態(tài)指數(shù)N,其實(shí)是從主溝道的侵蝕狀態(tài)來判斷流域發(fā)育的所處階段的。由S和N表達(dá)的地貌信息熵式(6)和式(8)讓S和N建立了一定的聯(lián)系,據(jù)蔣忠信的推導(dǎo)為:
S=1/(N+1)。
(15)
而這種關(guān)系是從兩種地貌熵的形式變換得到的,實(shí)際中S和N到底其關(guān)系是不是如此還需要實(shí)例來進(jìn)行驗(yàn)證,為此本文以白龍江流域?yàn)檠芯繀^(qū),通過區(qū)內(nèi)DEM進(jìn)行水文分析共提取出1 984條流域,通過歷史災(zāi)害資料的收集以及實(shí)地的調(diào)查共收錄有泥石流災(zāi)害記錄的泥石流流域714條(圖2),從中挑選流域內(nèi)典型的泥石流溝對其流域的面積高程積分值和流域縱剖面形態(tài)進(jìn)行對比。
圖2 白龍江泥石流溝流域分布圖
名稱H/mL/mS值N值R21/(N+1)誤差仁家溝157168860.690.7680.880.560.13馬泉溝146445930.650.8490.960.540.11陰山溝116422210.60.8780.970.530.07三眼峪溝243586700.551.5210.990.400.15羅家峪溝237677890.511.3180.990.430.07舊房溝189861240.492.4030.970.290.2紅土灣溝168969610.422.2500.940.310.11梨坪溝99523930.381.6160.980.380聯(lián)豐溝149751220.322.5310.990.280.04段河壩溝117745370.233.3990.990.230王石壩溝84035050.213.9130.990.200.01
從表2和圖3的結(jié)果我們可以看出,只有在縱剖面形態(tài)為比較規(guī)則的上凸或下凹曲線時(shí),利用公式擬合的效果才比較好,得到的N值也與S值有較好的一致性,比如陰山溝、羅家峪溝、梨坪溝、聯(lián)豐溝、段河壩溝、王石壩溝,而縱剖面的形態(tài)指數(shù)N會(huì)受到DEM精度、主溝道縱剖面精確性、縱剖面形態(tài)凹凸度、擬合公式局限性等眾多因素的影響且計(jì)算繁瑣,而利用DEM通過式(12)和式(13)來計(jì)算S值時(shí)受到其他因素的影響較小,結(jié)果較為可靠且計(jì)算便捷。
通過以上分析,使用S對泥石流溝發(fā)育階段進(jìn)行定量的劃分是比較便捷和可靠的。S劃分的幼年期、壯年期和老年期是針對侵蝕循環(huán)理論的,而對于泥石流溝的發(fā)育階段的劃分則要從泥石流發(fā)育特點(diǎn)來進(jìn)行重新劃分。從已有的研究來看泥石流的發(fā)育階段劃分為五個(gè)階段更為詳細(xì)和合理,分別為發(fā)育期、發(fā)展期、旺盛期、衰減期和消退期。
這里說的泥石流溝的發(fā)育時(shí)期應(yīng)該是泥石流在同等條件(即拋開隨機(jī)性較大的降雨等激發(fā)因素)下的泥石流的活躍程度或者說易發(fā)程度,所以并不能代表泥石流實(shí)際的活躍程度和規(guī)模大小,因?yàn)樵趯?shí)際條件下,泥石流的活躍程度受到物質(zhì)條件的動(dòng)態(tài)性、激發(fā)條件的隨機(jī)性等的影響,其活躍程度會(huì)有著較大的區(qū)別。
然而現(xiàn)實(shí)條件下無法找到同等條件下且發(fā)育在不同時(shí)期的泥石流溝,所以我們選取泥石流多發(fā)的白龍江流域[28]內(nèi)714條有泥石流記錄的溝進(jìn)行調(diào)查,得到550條泥石流溝的發(fā)生頻率資料(圖4),通過統(tǒng)計(jì)分析得到石流發(fā)生頻率與S值的關(guān)系(圖5)。
圖3 各泥石流溝溝谷縱剖面及擬合曲線
圖4 流域泥石流頻率圖圖5 泥石流溝的發(fā)生頻率與S值的關(guān)系圖
表3 泥石流發(fā)育階段劃分及其對應(yīng)泥石流發(fā)育特征
白龍江內(nèi)所有溝的S值范圍為0.15~0.8,有泥石流記錄的溝的S值范圍為0.3~0.7,且92%的溝S值集中在0.35~0.6之間,這也和地貌演化階段的壯年期相吻合,說明泥石流作為地貌演化中的一個(gè)過程是遵循地貌演化規(guī)律的,且發(fā)育狀態(tài)也和地貌演化狀態(tài)有著較好的一致性。泥石流溝相對較少的S的值域?yàn)?.3~0.4和0.6~0.7之間,據(jù)這幾個(gè)特征值我們將其劃分了五個(gè)發(fā)育階段,這五個(gè)階段及其對應(yīng)泥石流發(fā)育特征總結(jié)如表3所示。
本文總結(jié)了侵蝕流域系統(tǒng)的地貌信息熵與超熵的應(yīng)用現(xiàn)狀,通過分析地貌信息熵與超熵理論來源,認(rèn)為這些應(yīng)用的根本依據(jù)是其能夠反映泥石流溝的發(fā)育階段。通過S和N的計(jì)算方法以及其實(shí)際意義的分析認(rèn)為使用S來對泥石流發(fā)育階段進(jìn)行定量的劃分是比較便捷和可靠。通過泥石流發(fā)育頻率與S值的統(tǒng)計(jì)關(guān)系,初步將泥石流的發(fā)育階段劃分為發(fā)育期、發(fā)展期、旺盛期、衰減期和消退期五個(gè)階段,并分析了各個(gè)階段對應(yīng)泥石流的發(fā)育特征。
本文對泥石流的發(fā)育階段的研究只是一個(gè)初步的研究,但從一定程度上抓住了泥石流溝演化的一些規(guī)律,并在以后的研究中會(huì)逐漸深入和優(yōu)化。泥石流的發(fā)育階段的定量化研究以及各個(gè)發(fā)育階段的泥石流特征的分析,能為泥石流在各個(gè)發(fā)育階段的防治對策提供理論支持,同時(shí)也對潛在泥石流的判定、泥石流危險(xiǎn)性評價(jià)以及泥石流發(fā)生頻率的深入分析等研究提供支持。