氣候變化條件下降雨-融水型泥石流流量預(yù)測(cè)方法

1、氣候變化條件下降雨-融水型泥石流流量預(yù)測(cè)方法楊宗佶付校龍游勇Summary：全球氣候變化背景下，目前唯一的冰川泥石流流量經(jīng)驗(yàn)?zāi)Ｐ蜔o法定量 反映雨熱耦合變化的影響。根據(jù)度日模型冰雪融水當(dāng)量和泥石流流量雨洪計(jì)算 方法，提出了考慮雨熱耦合條件下降雨-融水型泥石流流量的計(jì)算模型。通過統(tǒng) 計(jì)51 a降雨和溫度觀測(cè)資料的變化趨勢(shì)，建立了帕隆藏布流域降雨-融水型泥 石流流量計(jì)算重現(xiàn)期標(biāo)準(zhǔn)，并通過對(duì)典型降雨-融水型泥石流流量的實(shí)測(cè)資料進(jìn) 行驗(yàn)證，說明該方法準(zhǔn)確性相比經(jīng)驗(yàn)?zāi)Ｐ兔黠@增加。進(jìn)一步以帕隆藏布流域15 條典型降雨-融水型泥石流為例，開展氣候變化條件下泥石流流量的計(jì)算和預(yù) 測(cè)。結(jié)果顯示：該方法可量化雨熱耦

2、合條件對(duì)泥石流流量的影響，可為氣候變 化條件下冰川泥石流的防治關(guān)鍵參數(shù)計(jì)算和預(yù)測(cè)提供參考依據(jù)。關(guān)鍵詞：降雨-融水型泥石流；泥石流流量；度日模型；雨熱耦合變化；氣候 變化；帕隆藏布流域中圖法分類號(hào)：P642. 23文獻(xiàn)標(biāo)志碼：ADOI： 10. 16232/j. cnki. 1001-4179. 2021. 12. 0060引言冰川泥石流廣泛分布于中國青藏高原1-2及世界各地的高山地區(qū)3-5。隨著 全球氣候變暖，冰川消融加劇6,將導(dǎo)致冰川泥石流頻繁發(fā)生。國內(nèi)外對(duì)冰川 泥石流展開的相關(guān)研究主要是成災(zāi)條件和機(jī)理方面。對(duì)于冰川泥石流的形成過 程，局部學(xué)者認(rèn)為冰植物在其中發(fā)揮了重要作用7-8,還有學(xué)者認(rèn)

3、為冰川泥石 流的形成是復(fù)合過程9,并將其啟動(dòng)過程總結(jié)為侵蝕-滑移型10。對(duì)于冰川 泥石流的形成條件，前人認(rèn)為其爆發(fā)與水熱組合關(guān)系密切，且多發(fā)生在濕熱環(huán) 境下1, 11;還有局部學(xué)者研究認(rèn)為其起動(dòng)和起動(dòng)方式與氣溫上升下冰川融水 12-13,黏粒含量14以及初始含水率15有關(guān)。但是目前有關(guān)冰川泥石流流 量的計(jì)算方法的研究較少，只有一套中科院成都山地所提出的經(jīng)驗(yàn)公式16, 無法反映溫度和降雨兩個(gè)關(guān)鍵因素對(duì)冰川泥石流流量的影響。因此本研究以帕 隆藏布流域?yàn)檠芯繀^(qū)域，以降雨與氣溫17為關(guān)鍵參數(shù)，構(gòu)建了疊加小流域暴 雨洪水計(jì)算公式和改進(jìn)度日模型18-20的降雨-融水型泥石流流量的計(jì)算模型 和預(yù)測(cè)方法，為氣

4、候變化背景下冰川泥石流的工程防治關(guān)鍵參數(shù)的選取和預(yù)測(cè) 預(yù)報(bào)提供了新的方法。1研究方法雨熱耦合條件下降雨-融水型泥石流流量的計(jì)算方法分別考慮全流域的降雨量和 冰雪流域的冰雪融水當(dāng)量。改進(jìn)冰雪融水當(dāng)量計(jì)算方法全流域24 h內(nèi)冰雪的消融水當(dāng)量Ml根據(jù)冰雪消融與氣溫之間關(guān)系建立的度日 模型18-20進(jìn)行計(jì)算：M1=DDF - T24 (S1F) (1)式中：T24為一天的平均氣溫，,考慮冰雪區(qū)域不同海拔位置的氣溫梯度， T24=T- (ELE-2 736) X0. 4721 ； S1為冰雪融水面積,km2；DDF為冰川或雪的度日因子，mm/ (d , ) o改進(jìn)的度日模型是在劉金平等22擬合的公式基礎(chǔ)

5、上進(jìn)行了改進(jìn)，考慮了坡度 的影響，由于冰雪的消融和其與空氣的接觸面積密切相關(guān)23,不同坡度投影 面積下冰雪區(qū)域與空氣的接觸面積不同，因而改進(jìn)公式如下：DDF= (0. 009XELE-O. 934XLAT-8. 1) Xcos 0 (2)式中：ELE代表海拔高度，m；LAT代表緯度，( ) ； 0為冰雪融水區(qū)域以上范 圍的平均坡度，( ) O降雨-融水型泥石流流量計(jì)算(1)改進(jìn)算法。改進(jìn)四川省中小流域暴雨洪水計(jì)算手冊(cè)中的計(jì)算公式，疊加融水和降雨 量，求得全流域在降雨-融水條件下的洪水洪峰流量。Q=0. 278 (H24+M1) 24n-l tr-hF (3)式中：Q為洪水洪峰流量，m3/s；H

6、24是全流域降雨量，采用24 h降雨量的統(tǒng)計(jì) 數(shù)值，mm；Ml為24 h融水量，mm；n為暴雨參數(shù)；t為流域匯流時(shí)間，h, T=0.278LmJl/3Ql/424；m為匯流參數(shù);J為流域平均縱比降;p為平均下滲強(qiáng) 度,mm/h； F為全流域面積，km2o本文采用配方法計(jì)算泥石流的峰值流量，其計(jì)算公式為Qc= (l+6c) QDu (4)4)c= Y c- Y w v s- v c (5)式中：Qc為泥石流洪峰流量，m3/s； 6c為泥石流洪峰流量修正系數(shù);Du為堵塞 系數(shù)；YC為泥石流容重，t/m3; Y w為清水比重，t/m3; Y s為泥石流中固體物 質(zhì)比重，t/m3o(2)傳統(tǒng)算法。目前

7、唯一能考慮冰川泥石流流量計(jì)算的經(jīng)驗(yàn)?zāi)Ｐ褪怯芍锌圃撼啥忌降厮岢龅?經(jīng)驗(yàn)公式16：Qc= (Q2+Q0) (l+4)c) d (6)式中：Q0為流域內(nèi)非冰川區(qū)的洪峰流量，m3/s；Q2為降雨型冰川消融流量，m3/s, Q2=F1 (0.05H+2. 1) , H為降雨量,mm；d為冰川消融洪峰系數(shù),d=l+7. 6 (F1/F) +0.05。O；F1為冰川面積，km2；0 0為冰川坡度，( )。2數(shù)據(jù)來源 帕隆藏布流域是中國最典型的海洋性冰川區(qū)和雅魯藏布江降水中心且具有雨熱 同期的特點(diǎn)，升溫和暴雨成為誘發(fā)大型泥石流的有利條件。本文對(duì)波密縣19612011年51 a的降雨和氣溫?cái)?shù)據(jù)通過皮爾遜TH型

8、曲線進(jìn)行概率分析統(tǒng)計(jì) 計(jì)算（見圖12）,在皮爾遜-III型曲線的分析中得到在不同頻率下年最高日 平均氣溫以及日降雨量如表1所列。由于波密縣的氣溫呈現(xiàn)整體上升的趨勢(shì)，因此符這51 a的日最高平均氣溫?cái)M合 如圖2所示，在對(duì)20n年以前的泥石流實(shí)例驗(yàn)證中的溫度指標(biāo)減去了這個(gè)升溫 趨勢(shì)，公式為TPa=TP-0.03612X （2011-a） ；TPa為第a年頻率為P的年最高日 平均氣溫，。C；TP為波密縣統(tǒng)計(jì)不同頻率年最高日平均氣溫，。C；a為年份。對(duì) 以后泥石流的預(yù)測(cè)那么是在不同情景條件下疊加升溫趨勢(shì)，公式為TPn=TP+0. 03612 Xn, TPn為n年一遇的年最高日平均氣溫，。C。3模型驗(yàn)證

9、本文選取了帕隆藏布流域前期研究和實(shí)測(cè)數(shù)據(jù)相對(duì)較豐富的古鄉(xiāng)溝和培龍溝兩 條典型冰川泥石流開展實(shí)例驗(yàn)證（見圖3）。歷史統(tǒng)計(jì)數(shù)據(jù)如表2所列。通過皮爾遜TH型曲線計(jì)算不同頻率下的溫度TP以及TPa得到1983, 1984年和 2005年的不同頻率下的最高日平均氣溫如表3所列。通過改進(jìn)算法和傳統(tǒng)算法計(jì)算古鄉(xiāng)溝與培龍溝不同頻率下的泥石流數(shù)據(jù)如表 45所列。改進(jìn)算法和傳統(tǒng)算法在不同頻率下計(jì)算得到的古鄉(xiāng)溝和培龍溝泥石流流量與實(shí) 測(cè)數(shù)據(jù)比照方圖45所示。雨熱耦合條件下降雨-融水型泥石流流量的計(jì)算 方法分別考慮全流域的降雨量和冰雪流域的冰雪融水 當(dāng)量。改進(jìn)冰雪融水當(dāng)量計(jì)算方法全流域24 h內(nèi)冰雪的消融水當(dāng)量M根據(jù)

10、冰雪消 融與氣溫之間關(guān)系建立的度日模型8 - 2。進(jìn)行計(jì)算：-S.以二DDF-心5)r式中廳24為一天的平均氣溫，工，考慮冰雪區(qū)域不同海 拔位置的氣溫梯度，= T - (ELE - 2 736 ) x 0.47；S為冰雪融水面積雪m?；DOF為冰川或雪的 度日因子，mm/(d工)。改進(jìn)的度日模型是在劉金平等口擬合的公式基礎(chǔ) 上進(jìn)行了改進(jìn)，考慮了坡度的影響，由于冰雪的消融和其如表23所列，以溫度為衡量指標(biāo)，古鄉(xiāng)溝2005年約為27 a 一遇的泥石流， 由改進(jìn)算法求得的流量為831.33 m3/s,與前人25計(jì)算的33 a一遇的泥石流 規(guī)模849. 65 m3/s比擬符合。以溫度為衡量指標(biāo)，培龍溝

11、1983年和1984年的 泥石流頻率分別為36 a一遇和154 a一遇，由改進(jìn)算法求得的流量分別為2 808 m3/s與5 259 m3/s,與歷史記載流量2 950 m3/s與5 245 m3/s相似度較 高。而由圖45和表45比照分析可知，傳統(tǒng)算法的計(jì)算結(jié)果普遍偏大，改 進(jìn)算法對(duì)古鄉(xiāng)溝和培龍溝的計(jì)算準(zhǔn)確性相比傳統(tǒng)方法分別提高47. 26%與 21. 12%32. 5%o根據(jù)歷史資料揭示，古鄉(xiāng)溝泥石流的規(guī)模在1953年以后在逐漸減小，且20世 紀(jì)70年代以后泥石流流量都普遍小于1 000 m3/s28,而傳統(tǒng)經(jīng)驗(yàn)算法計(jì)算 古鄉(xiāng)溝20 a一遇流量達(dá)1 156 m3/s,嚴(yán)重偏大，相比之下改進(jìn)算

12、法精度更 高。4氣候變化條件下泥石流預(yù)測(cè)1波密地區(qū)泥石流溝基本數(shù)據(jù)本文選擇了 15條典型降雨-融水型泥石流溝進(jìn)行氣候變化條件下泥石流流量預(yù) 測(cè)（見圖3）,基礎(chǔ)數(shù)據(jù)如表6所列。波密地區(qū)泥石流冰雪消融量波密地區(qū)的泥石流冰雪融水當(dāng)量采用度日模型公式（1）進(jìn)行計(jì)算。以表1的 不同頻率下的年最高日平均氣溫為標(biāo)準(zhǔn)，疊加最近51 a升溫趨勢(shì)下的日平均氣 溫與降雨量，如表7所列。計(jì)算得到的全流域24 h冰雪融水當(dāng)量如表8所列。波密地區(qū)泥石流溝計(jì)算結(jié)果 結(jié)合全流域24 h冰雪融水量和降雨量，通過傳統(tǒng)算法（式（6）和改進(jìn)算法（式（3）、（4）計(jì)算不同頻率下泥石流溝的泥石流流量的結(jié)果如表9和圖6所示。將上述15條泥

13、石流溝按冰雪融水面積占全流域面積大于1/4和小于1/4分為兩 類（見圖7和圖8）。由上述圖表可知，傳統(tǒng)經(jīng)驗(yàn)算法的計(jì)算結(jié)果很大程度上取決于冰川區(qū)域，在冰 川面積占比大于1/4的區(qū)域（見圖7）,改進(jìn)算法的計(jì)算結(jié)果稍低于偏向保守 的傳統(tǒng)算法。而在冰川面積占比小于1/4的區(qū)域（見圖8）,傳統(tǒng)算法計(jì)算的 結(jié)果那么遠(yuǎn)小于改進(jìn)算法，因而在冰川占比擬小的區(qū)域，傳統(tǒng)算法的結(jié)果那么會(huì)過 于偏小。而現(xiàn)如今冰川的面積在逐年減少29,所以相比傳統(tǒng)依賴于冰川面積 的經(jīng)驗(yàn)算法，改進(jìn)算法在今后預(yù)測(cè)氣候變化條件下泥石流流量的實(shí)用性和準(zhǔn)確 性更強(qiáng)。5結(jié)論和討論（1）建立了降雨-融水型泥石流溫度頻率確定標(biāo)準(zhǔn)，構(gòu)建了波密地區(qū)氣候變化

14、 條件下的升溫趨勢(shì)，提出了基于改進(jìn)的度日模型和雨洪法的改進(jìn)降雨-融水型泥 石流流量算法。（2）以古鄉(xiāng)溝和培龍溝為例，改進(jìn)算法的結(jié)果與歷史數(shù)據(jù)比照，誤差分別為2. 24%與0. 27%4. 85%,準(zhǔn)確性比傳統(tǒng)經(jīng)驗(yàn)?zāi)Ｐ吞岣吡?47. 26%與21. 12% 32. 5%o（3）對(duì)于冰川占比大于1/4的4條溝，傳統(tǒng)算法與改進(jìn)算法誤差在3%20%之間，由于傳統(tǒng)算法較為保守，因此傳統(tǒng)算法的結(jié)果均大于本文的改進(jìn)算法, 所以傳統(tǒng)的經(jīng)驗(yàn)算法對(duì)于冰川面積占比大于1/4的泥石流流域仍然具有不錯(cuò)的 參考價(jià)值。但是對(duì)于冰川面積占比擬小的泥石流溝，計(jì)算結(jié)果那么會(huì)偏低甚至無 法作為參考，因此改進(jìn)算法適用性更廣。(4)本

15、研究在用度日模型計(jì)算冰雪融水的時(shí)候疊加了波密縣近51 a來的升溫 趨勢(shì)，使預(yù)測(cè)結(jié)果更加符合氣候變化的趨勢(shì)。Reference：1呂儒仁，李德基.西藏波密冬茹弄巴的冰雪融水泥石流J.冰川凍土， 1989, 11 (2) ： 148-160.GE Y G, CUI P, SU F H, et al. Case history of the disastrous debris flows of Tianmo Watershed in Bomi County, Tibet, China： some mitigation suggestionsJ. Journal of Mountain Science

16、, 2014, 11(5) ： 1253-1265.JACKSON L E, HUNGR 0, GARDNER J S, et al. Cathedral Mountain debris flows, CanadaJ. Bulletin of Engineering Geology and the Environment, 1989, 40 (1) ： 3554.WILKERSON F D, SCHMID G L. Distribution of debris flows in Glacier National Park, Montana, U. S. A. J. Journal of Mou

17、ntain Science, 2008, 5 (4) ： 318-326.HUGGEL C, HAEBERLI W, KAAB A, et al. An assessment procedure for glacial hazards in the Swiss AlpsJ. Canadian Geotechnical Journal, 2004, 41 (6) ： 1068-1083.6何秋樂，匡星星，梁四海，等.19662015年長江源冰川融水變化及其對(duì)徑流的影響：以冬克瑪?shù)缀恿饔驗(yàn)槔齁.人民長江，2020, 51 (2) ： 77-85, 130.7施雅風(fēng)，楊宗輝，謝自楚，等.西藏古鄉(xiāng)地區(qū)

18、的冰川泥石流J.科學(xué)通報(bào)， 1964, 9 (6) ： 542-544.8杜榕桓，章書成.西藏高原東南部冰川泥石流的特征J.冰川凍土，1981, 3 (3) ： 10-16.9 FLEISCHMAN C M. Debris FlowM. Beijing： Science Press, 1985.10 RUOFU S. Basic characteristics and measuring method of debris flowM. Chongqing： Science and Technology Literature Press Chongqing Branch, 1986.11齊云龍，

19、鄧明楓.川藏公路波密段九絨溝泥石流形成機(jī)制研究J.災(zāi)害 學(xué),2019, 34 (3) ： 123-127.12陳曉清，崔鵬，楊忠，等.聶拉木縣沖堆普2002年泥石流成因分析及防治 對(duì)策J.冰川凍土, 2006, 28 (5) ： 776-781.ZIMMERMANN M, HAEBERLI W. Climatic change and debris flow activity in high-mountain areas： a case study in the Swiss Alps J. Catena. Supplement, 1992, 22 (1) ： 59-72.14陳寧生，王政，田樹

20、峰，等.冰磺土體起動(dòng)泥石流的特征研究J.第四紀(jì)研究,2019, 39 (5) ： 1235-1245.15潘蕾，魏學(xué)利，張遠(yuǎn)芳，等,初始含水率對(duì)冰川泥石流的起動(dòng)影響分析 J.水土保持學(xué)報(bào)，2017, 31 (6) ： 116-122.16羅德富，毛濟(jì)周.川藏公路南線(西藏境內(nèi))山地災(zāi)害及防治對(duì)策國.北 京：科學(xué)出版社，社95.17屈永平，肖進(jìn)，潘義為.藏東南地區(qū)冰川泥石流形成條件初步分析：以天 摩溝冰川泥石流為例J.水利水電技術(shù)，2018, 49 (12) ： 180-187.SINGH P, KUMAR N, ARORA M. Degree-day factors for snow and

21、ice for Dokriani Glacier, Garhwal HimalayasJ. Journal of Hydrology, 2000, 235 (1) ： 1-11.HOCK R. A distributed temperature-index ice-and snowmelt model including potential direct solar radiationJ. Journal of Glaciology, 1999, 45 (149) ： 101-111.HOCK R. Temperature index melt modelling in mountain ar

22、easJ. Journal of Hydrology, 2003, 282 (1/4) ： 104-115.21登K明楓.帕隆藏布流域冰川降雨型泥石流頻發(fā)成因分析D.北京：中國科 學(xué)院大學(xué)，2019.22劉金平，張萬昌.雅魯藏布江流域度日因子空間變化J.中國科學(xué)院大學(xué) 學(xué)報(bào),2018, 35 (5) ： 131-138.23 HOCK R. Glacier melt： a review of processes and their modellingJ. Progress in Physical Geography, 2005, 29 (3) ： 162-391.24熊文華，孫立霞，李誠，等,

23、瀾滄江與怒江流域洪峰流量公式推求J .人 民長江,2020, 51 (11) ： 82-86.25魯安新，鄧曉峰，趙尚學(xué)，等.2005年西藏波密古鄉(xiāng)溝泥石流爆發(fā)成因分 析J,冰川凍土, 2006, 28 (6) ： 956-960.26程尊蘭，吳積善.西藏東南部培龍溝泥石流堵塞壩的形成機(jī)理A.第八屆 海峽兩岸山地災(zāi)害與環(huán)境保育學(xué)術(shù)研討會(huì)C 烏魯木齊：中國水土保持學(xué)會(huì)， 2011： 39-45.27國家防汛總指揮部辦公室，中國科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究 所.山洪泥石流滑坡災(zāi)害及防治M.北京：科學(xué)出版社，1994.28朱平一，羅德富.西藏古鄉(xiāng)溝泥石流開展趨勢(shì)J.山地學(xué)報(bào)，1997, 15

24、(4) ： 296-299.29劉時(shí)銀，姚曉軍，郭萬欽，等.基于第二次冰川編目的中國冰川現(xiàn)狀J. 地理學(xué)報(bào),2015, 70 (1) ： 3-16.(編輯：劉媛)Discharge prediction method for rainfall-melting water triggered glacial debris flows under climate change scenariosYANG Zongjil, FU Xiaolongl, 2, YOU Yongl(1. Institute of Mountain Hazards and Environment, CAS, Chengdu

25、 610041, China； 2. University of Chinese Academy of Sciences, Beijing 100049, China)Abstract：Under the background of global warming, current empirical model for the glacial debris flow discharge calculation cannot quantitatively consider the coupling effect of rainfall water and the snow/ice melting

26、 water. We proposed a discharge calculation model that considers coupling effect of both rain and melting water for glacial debris flows by studying ice-snow melt equivalent in the improved degree-day model as well as the current storm flood calculation method for rainfall triggered debris flows. The criteria of return period for the glacial debris flows discharge prediction was established based on the rainfall and temperature change trends fitted by the 51 years long period observation data in the Palong Zangbo watershed. Then this prop