官壩河泥石流動力學特征參數

李 俊，陳寧生

（1．四川大學 水力學與山區河流開發保護國家重點實驗室，成都610065；2．中國科學院 水利部 成都山地災害與環境研究所，成都610041）

25萬km2的青藏高原東南緣內發育有11個構造斷陷湖［1］，這些構造斷陷湖都不同程度地受到泥石流侵蝕帶來的泥沙淤積的影響。滇池、瀘沽湖、洱海、邛海均地處構造活躍區，地震頻繁，人為生態破壞嚴重，使得這些區域的山洪和泥石流發育［2－9］。山洪和泥石流攜帶大量的泥沙進入湖泊，使其呈現嚴重萎縮的趨勢。目前人們已經認識到僅修建河堤或排導槽并不能解決青藏高原東南緣湖泊的泥沙淤積問題，還必須進行有效的工程治理。青藏高原東南緣眾多的高原構造斷陷湖中，邛海的泥沙淤積問題最具有代表性。日益嚴峻的泥沙淤積使邛海面臨著生死存亡的憂患［10－11］。據調查，西昌市邛海北岸的官壩河平均每年流入邛海泥沙量為15萬m3，并且泥沙每年平均向邛海推進的距離為55m。同時，邛海南岸的鵝掌河將大量泥沙帶入邛海［12］，鵝掌河輸入邛海的泥沙有4．96萬m3，并且泥沙每年平均向邛海推進距離為2m［6］。因此，本文通過總結國內外泥石流防治工程的研究現狀［13］，選取治理難度較大的典型泥石流防治工程，如邛海官壩河泥石流防治工程，采用合理的方法計算其泥石流動力學特征參數，并給出官壩河泥石流防治建議。以期為治理青藏高原東南緣構造斷陷湖區域內的山洪和泥石流提供參考。

1 官壩河流域概況

官壩河系邛海的一級支流，位于西昌市川興鎮東北方，屬青藏高原東南部螺髻山系，東連昭覺縣，北接喜德縣，西與著名的瀘山風景區遙相對應。早更新世末，沿川滇菱形斷塊的東側邊界斷裂曾發生強烈的拉張性斷陷活動，并在西昌附近形成了一個北北西向的地塹地壘系，其中有大興單斷同生斷陷盆地［14］，受川滇南北向構造作用影響，官壩河流域地形地貌較為復雜。海拔2 700m以上區域主要為中高山剝蝕地貌，約占流域總面積的27．7%；海拔高程為1 700～2 700m區域主要為深切峽谷，約占流域總面積的59．9%，溝道多呈“V”字型，滑坡、崩塌和坡面泥石流發育；海拔高程1 500～1 700m地帶斷陷盆地發育。經過現場調查，推測張巴寺河上游處可能有一斷層，走向呈北西—南東向，傾角50°左右，該斷層經過地區巖體破碎，崩塌滑坡發育，為泥石流的發生提供了大量物源。官壩河基本情況見表1。

官壩河物源豐富，總物源為1．2×107m3，物源動儲量為1．16×106m3。動物源主要集中于張巴寺河上游。官壩河流域上游裸地區產生了大量的細顆粒，經分析認為這些細顆粒以山洪泥石流方式進入邛海。

表1 官壩河基本參數

2 官壩河泥石流活動特征和邛海泥沙淤積特征

2．1 1998年官壩河泥石流特征值參數

2．1．1 容重 2010年5月調查了官壩河流域的泥石流歷史活動（表2）。通過調查，發現三處明顯的1998－07－06官壩河泥石流洪痕，分別為1—1′，2—2′，3—3′斷面。但由于后期水流改造，難以采集到泥石流原樣，根據現場顆分實驗結果（圖1），采用基于黏粒含量的泥石流容重計算方法［15］和根據泥石流沉積物計算泥石流容重計算方法［16］，取這兩種計算結果的平均值作為官壩河泥石流容重值。計算結果表明，張巴寺河沿程容重逐漸降低，泥石流漿體運動至斷面2—2′處以后，容重值由2．33 g／cm3降低至1．69g／cm3，官壩河泥石流由黏性過渡到稀性。官壩河其他支流如麻雞窩河和新仁寺河容重值分別為1．79g／cm3和1．93g／cm3，屬于過渡性泥石流，而大蘿卜溝泥石流容重值為1．65g／cm3，屬于稀性泥石流。官壩河流域泥石流總體屬于過渡性—黏性泥石流［17］。

表2 官壩河泥石流歷史活動情況和容重

2．1．2 流量 采用雨洪法［18］計算各斷面處的泥石流洪峰流量（表3），然后用泥石流形態調查法［18］進行校核（表4）。

計算結果表明，雨洪法和形態調查法計算的泥石流峰值流量較為接近，官壩河泥石流的最大洪峰流量出現在1—1′斷面處為813．7m3／s，比斷面7—7′的泥石流洪峰流量大。分析認為，由于1—1′斷面以下的河流為寬淺泥石流灘地，沖出出山口后斷面逐漸展寬，泥石流流速降低，泥石流漿體沿程停淤，最后導致泥石流在入?？谔幜髁拷档?。

圖1 官壩河各斷面處沉積物顆分曲線

2．1．3 一次泥石流固體物質總量 目前國內外還沒有好的方法計算泥石流固體總量，因此采用泥石流五邊形法［18］，粗略計算了1998年7月6日官壩河斷面7—7′處一次泥石流固體物質總量（表5）。

通過大量的實地走訪，確認1998年7月6日官壩河泥石流為百年一遇的泥石流。計算結果表明，1998年7月6日官壩河入?？谔幰淮文嗍鞴腆w物質總量46．4萬m3，另外通過現場實地調查，官壩河入?？谔?998年泥石流堆積物為69．0萬m3，故可判斷剩余22．6萬m3泥石流體來源于同年的其他小規模山洪泥石流和常流水期間被水流所攜帶的懸移質。說明1998年官壩河進入邛海的泥沙以山洪泥石流方式運移的泥沙總量最大，因此官壩河泥石流的工程治理對減少邛海泥沙淤積問題具有重要意義。

表3 雨洪法計算的泥石流洪峰流量

表4 形態調查法計算的泥石流洪峰流量

表5 官壩河泥石流斷面7－7′的一次泥石流總量和一次泥石流固體物質總量

2．2 1952—2003年邛海泥沙淤積特征

據調查與測量，邛海水面面積從1952年的31 km2減少到2003年的27．4km2，按50a計算，則邛海平均每年約減少7．2萬m2，50a邛海面積萎縮11．6%。邛海的庫容由1952年的3．2億m3減少到2003年的2．93億m3（表6），最大水深和平均水深也由原來的34．0m和14．0m降低至2003年的18．32 m 和10．95m［12］。

表6 邛海典型年份淤積變化

據表6，邛海面積和庫容變化率為0．16%，水深和平均水深變化率分別為0．9%和0．43%，面積、庫容和水深變化都說明邛海正在不斷地萎縮，按此變化率計算，則邛海的壽命將只有600a左右。

據調查，2010年的邛海面積為27km2，泥沙淤積現狀嚴重（圖2）。2003—2010年，邛海水面面積平均每年約減少5萬m2，8a內邛海水面面積萎縮了1．45%。近年來邛海流域在小箐溝、鵝掌河實施的水土保持工程對入海泥沙起到了一定的控制作用，其他未進行工程治理的流域依然十分嚴重。2008—2009年官壩河河口泥沙推進90m，高倉河河口泥沙前進20m。通過對邛海1952—2003年泥沙淤積特征的分析和官壩河98年百年一遇泥石流的動力學特征參數計算結果可知，僅通過修建河堤或排導槽并不能解決邛海的泥沙淤積問題，還須依賴于有效的工程治理。

2．3 官壩河泥石流動壓力和沖擊力

依據舟曲泥石流的大石塊沖擊力計算實例［19］，基于1998年7月6日官壩河泥石流動壓力和大石塊沖擊力，分別采用鐵二院公式［20］和基于彈塑性理論的泥石流大石塊沖擊力計算方法［21］（表7），計算擬建混凝土攔砂壩的動壓力和大石塊沖擊力。現場調查1998年官壩河泥石流最大石塊直徑為6m。官壩河溝道堆積大石塊巖石類型為紅砂巖，大石塊密度取2．31g／cm3。

圖2 2010年官壩河入海口處泥沙現狀

表7 官壩河泥石流動壓力和大石塊沖擊力

3 官壩河泥石流防治建議

總結了官壩河泥石流的活動歷史、動力學特征參數和邛海泥沙淤積特征，建議官壩河泥石流的防治原則應以保護官壩河上游居民和保護邛海為主，盡量減少官壩河上游滑坡和崩塌物產生的細顆粒進入邛海，結合官壩河溝道地形地貌特征、物源量和動力學特征參數，建議在官壩河主溝及各支溝下游處設置以攔擋細顆粒為主的攔砂壩和修建排導槽并定時清淤，以減少進入邛海的泥沙。同時在主溝及各支溝上游有施工條件的地方盡量多修建谷坊，以形成谷坊壩群效應，起到穩固溝床，攔擋大顆粒作用，在保護官壩河上游居民的安全，同時減輕下游攔砂壩的攔擋壓力。另外考慮了地震和極端降雨情況，建議攔砂壩和谷坊采用高強度混凝土，確保其施工質量。

4 討 論

官壩河泥石流容重計算采用了兩種截然不同的方法，這兩種方法在計算官壩河泥石流最大容重值時存在較大波動，而計算其他支溝泥石流容重偏差較小。這是由于目前基于物體顆粒級配的泥石流容重計算的不成熟和青藏高原東南緣構造斷陷湖泥石流溝的特殊性引起的結果，即泥石流搬運的粗大顆粒堆積在溝床中，而泥石流演化成高含沙水流以后細顆粒隨洪水進入湖中，調查中無法獲得更全面的顆分資料。因而進一步研究青藏高原東南緣構造斷陷湖的泥石流溝容重值變化過程存在困難。另外1998年官壩河向邛海以山洪泥石流方式運移的泥沙總量占當年邛海泥沙淤積總量的67．2%。說明泥石流對邛海可持續利用非常不利。最后總結了官壩河百年一遇的清水流量497m3／s與集水面積106．64km2成5倍關系，其他斷面的百年一遇清水流量與面積之比在4～7倍之間浮動，可為計算川滇地區其他構造斷陷湖的泥石流溝提供經驗估算公式。

［1］ 張永雙，曲永新，趙希濤，等．青藏高原東南緣第四紀工程地質概論［M］．北京：地質出版社，2010．

［2］ 周靜，王蘇民，呂靜．洱海地區一萬多年以來氣候環境演化的湖泊沉積記錄［J］．湖泊科學，2003，15（2）：104－111．

［3］ 張振克，沈潔，羊向東，等．近8ka來云南洱海湖泊沉積記錄的氣候變化與夏季印度季風強弱變化的關系［J］．亞熱帶資源與環境學報，2008，3（3）：1－6．

［4］ 張振克，吳瑞金，沈潔，等．近1800年云南洱海流域氣候變化與人類活動的湖泊沉積記錄［J］．湖泊科學，2000，12（4）：297－303．

［5］ 姚遠，張恩樓，沈潔，等．云南屬都湖流域人類活動的湖泊沉積響應［J］．海洋地質與第四紀，2007，27（5）：115－120．

［6］ 余斌，王士革，章書成，等．鵝掌河泥石流對四川邛海影響的初步研究［J］．湖泊科學，2006，18（1）：57－62．

［7］ 朱海虹．云南斷陷湖泊的形成和晚新生代的沉積及其演化［J］．海洋與湖沼，1991，22（6）：509－516．

［8］ 朱海虹，姚秉衡．云南斷陷湖泊三角州沉積及其在石油勘探中的意義［J］．石油與天然氣地質，1989，10（2）：95－106．

［9］ 王小雷，楊浩，趙其國，等．云南蕪仙湖近現代環境變化的沉積物粒度記錄［J］．沉積學報，2010，28（4）：776－782．

［10］ 袁國林，李昌侯，米鴻雁，等．邛海湖盆形態研究［J］．海洋湖盆通報，2006（2）：18－26．

［11］ Liu Yong，Zhou Feng，Guo Huaicheng，et al．Biotic condition assessment and the implication for lake fish conservation：a case study of Lake Qionghai，China［J］．Water and Environment Journal，2009，23（3）：189－199．

［12］ 余斌，章書成，王士革．四川西昌邛海的濁流沉積初探［J］．沉積學報，2005，2（4）：559－565．

［13］ 蔣樹，文寶萍．國內外泥石流活動關鍵指標估算方法之比較［J］．水文地質工程地質，2012，39（3）：86－96．

［14］ 聞學澤，吳迪忠．四川西昌邛海第四紀盆地成因的探討［J］．地質評論，1985，31（4）：360－366．

［15］ 陳寧生，崔鵬，劉中港，等．基于黏土顆粒含量的泥石流容重計算［J］．中 國科學：E 輯，2003，33（增刊）：164－174．

［16］ 余斌．根據泥石流沉積物計算泥石流容重的方法研究［J］．沉積學報，2008，26（5）：789－796．

［17］ 崔鵬，莊建琦，陳興長，等．汶川地震區震后泥石流活動特征與防治對策［J］．四川大學學報：工程科學版，2010，42（5）：10－19．

［18］ 陳寧生，楊成林，周偉，等．泥石流勘查技術［M］．北京：科學出版社，2011．

［19］ 鄧虎，陳寧生，胡桂勝，等．甘肅舟曲三眼峪溝泥石流動力學特征參數計算［J］．重慶交通大學學報：自然科學版，2011，30（4）：833－838．

［20］ 周必凡，李德基，羅德富．泥石流防治指南［M］．北京：科學出版社，1991．

［21］ 何思明，吳永，沈均．泥石流大塊石沖擊力的簡化計算［J］．自然災害學報，2009，18（5）：51－56．