基于HMA方法的宜昌站生態徑流演變規律分析

胡 楊，林凱榮，劉 攀

(1.中山大學水資源與環境系，廣州 510275; 2.華南地區水循環與水安全廣東省普通高校重點實驗室，廣州 510275;3 武漢大學水利水電學院， 武漢 430072)

1 研究區概況

長江發源于青海省唐古拉山脈各拉丹冬峰西南側,河源名為沱沱河,從河源到四川宜賓稱為金沙江,岷江匯入之后始稱川江。湖北宜昌以上河段為長江上游，橫跨中國陸地地勢的一、二級階梯。該段河道全長4 529 km，占長江主河道全長的70%，上游沿岸有數個重要城市，包括宜賓、瀘州、重慶、宜昌等，主要的一級支流包括岷江、嘉陵江、烏江等。

宜昌水文站始于1877年設立的宜昌海關，此后開始有水位記載。1946年3月正式設立宜昌水文站，為長江上游出口控制站[1]。宜昌站位于西陵峽口的宜昌市，宜昌因水而興旺發達，但是也飽受水害的摧殘，屬亞熱帶季風性濕潤氣候，四季分明，水熱同季，寒旱同季的氣候特征，多年平均降水量1 215.6 mm。地處長江上中游的咽喉之地，宜昌站控制流域面積100萬km2，占全流域面積的55%，為下游防洪以及葛洲壩、三峽工程等大型工程提供了詳盡的水文資料，也為宜昌市的城市建設和經濟建設起到了舉足輕重的作用。上游數公里處的葛洲壩水電站，是長江干流上第一座大型水電站，葛洲壩水利樞紐于1970年開始動工，宜昌水文站的徑流量資料可以體現出葛洲壩水電站的建設對長江天然徑流量的影響程度。因此本文采用Mann-Kendall法和基于IHA指標體系的HMA法，分析宜昌站徑流量年際變化的規律和趨勢，為水資源管理部門和工程建設部門的規劃提供依據，便于水庫調度的決策，保護上中游流域生態環境。

2 研究方法

2.1 Mann-Kendall檢驗

Mann-Kendall檢驗法[2,3]是世界氣象組織推薦并已廣泛使用的非參數檢驗方法，不需要遵從一定的分布，也不受少數異常值的干擾，適用于水文、氣象等非正態分布的數據，計算比較簡便[4]。

2.2 水文變化指標法(IHA法)

Richter等人(1996年, 1998年)提出的受自然保護協會(TNC)支持的IHA統計包[5,6]計算水文變異。IHA方法在研究人類活動如大壩和地下水取水等對水文情勢的影響方面非常有用。IHA法[5]是依據河流的日測水文流量資料，計算32個影響生態環境意義的關鍵水文特征值，并計算集中量數(如中位數、平均值)和離散量數(如標準偏差、變異系數)，以對受人類活動影響的河流進行干擾前和干擾后的對比分析[7]。32個水文關鍵值可以分為5組，分別是月徑流量大小、年極端流量大小、極端流量出現時間、高低流量的頻率和延時以及徑流變化的改變率和頻率[8]。流量平均增加率、減少率指相鄰兩天的流量平均增加率、減少率，逆轉次數是指河流日流量由增加變為減少或者由減少變為增加的次數[9]。

2.3 直方圖匹配法(HMA法)

HMA法是Shiau and Wu[10]提出的一種以直方圖差異度作為衡量標準評價水文情勢改變度的方法。HMA的核心思想是：若兩個徑流序列的32個IHA指標的頻率直方圖相互類似，則這兩個徑流序列的水文情勢也類似，這種相似性可以用變異前后頻率直方圖的“統計距離”來衡量。因此，水文情勢的改變可以用基于距離的差異度來度量。HMA的過程如下：

(1)

式中：Dm是第m個指標的直方圖差異度；dm是第m個指標的二次方間距；max(dm)是每個dm的最大值。

(2)

式中：h=(h1，h2，…，hnc)T，k=(k1，k2，…，knc)T，分別為變異前后直方圖H、K的頻率向量， |h-k|是統計距離向量，nc是徑流序列的分類數；A=[aij]為相似矩陣，aij是H的第i類和K的第j類之間的相似性：

(3)

其中，dij=|Vi-Vj|是第i，j類之間的距離，Vi、Vj分別為i、j類的均值；dmax=max(dij)=|V1-Vnc|；α為(1，∞)，本文采用α=1的線性相似函數。

(4)

式中：r為實測數據中最大值與最小值之差；n為實測數據長度；riq為四分位數間距。

(5)

式中：nc,m為第m個IHA指標的分類值。

因Dm的個別值可以表現出不同程度的差異性，所以用一個綜合指數來評價“水文情勢的整體變化程度”，即：

(6)

式中：D′o最小值被視為等同于保護最好的天然徑流量，用以替代生態系統需水量。

3 結論與分析

3.1 Kendall趨勢分析

由于人類活動和氣候變化的因素,使得宜昌站的年徑流數據產生變化,其中上游的大型水利設施對徑流的影響非常顯著。圖1可以體現宜昌站年徑流量的變化趨勢和變異特點，從圖中可知宜昌站的年平均流量呈降低的趨勢，宜昌站的變異點在1970年左右和1996年左右，徑流量的變異應該和長江上游1971年開工建設的龔嘴水庫、黑龍灘水庫，特別是葛洲壩水利樞紐等工程有關。而這也與秦年秀、曾娟等人的研究結果相近，并符合實際狀況。與此同時，由于經濟的增長，對灌溉用水和工業用水的需求相應地大規模增加，此外湖泊的開墾、氣候變化等都會對年徑流量產生不同程度的影響。因此將1878-1970年作為長江上游徑流量未受擾動時間段，將1971-2010年作為由于上游水利設施而產生擾動的時間段。

圖1 宜昌站徑流變化和突變分析圖

3.2 IHA指標變化分析

運用IHA軟件將各站逐日實測徑流資料分為月平均流量、年極端流量、極端流量發生時間、高低流量的頻率和延時以及徑流變化的改變率和頻率5組共32個指標。通過每一個指標在變異前后的年平均值分析變異趨勢，結果顯示，除了4月的月平均流量略微增加外，其他月的月平均流量都有不同程度的減少；年最大和最小極值流量都有小幅度下降的趨勢；每年發生低流量的次數和延時增加，每年發生高流量的次數和延時減少；每年的流量逆轉次數增加。從上述這些指標的改變趨勢可以看出，包括葛洲壩在內地一系列在長江上游修建的水利設施確實對天然徑流量產生了作用，這主要是因為上游水庫的蓄水和蒸發使得徑流量減少，而關閘蓄水和開閘放水的使得逆轉次數增多。

3.3 用HMA法作徑流量變異分析

3.3.1月平均流量

我們定義HMA方法中，改變度高于67%為高度改變，改變度在33%和67%之間為中度改變，改變度小于33%為低度改變。中度改變的月平均流量指標有8月、10月和11月，見圖2。HMA法這3個月的差異度分別為40.39%、33.77%和37.15%。

在8月份，天然徑流序列小于13 308 m3/s的頻率值為0，在13 308～17 025 m3/s的頻率為0.02，在17 025～20 743 m3/s的頻率為0.12，在20 743～28 178 m3/s的頻率為0.33，變異后，徑流小于9 590 m3/s的頻率為0.07，在9 590～13 308 m3/s的頻率為0.07，在17 025～20 743 m3/s的頻率為0.2，在20 743～28 178 m3/s的頻率上升為0.71。

在10月份，天然徑流序列流量11 882～14 064 m3/s的頻率值為0.03，14 064～16 245 m3/s的頻率為0.15，小于16 245 m3/s的頻率為0.23。變異后月平均徑流序列在11 882～14 064 m3/s的頻率值升高至0.26，在14 064～16 245 m3/s的頻率升高為0.45，小于16 245 m3/s的頻率升高為0.75，沒有高于27 155 m3/s的流量。

在11月份，變異前流量小于8 912 m3/s的頻率為0.38，變異后卻漲至0.88。

上述結果表明，高要站8月、10月、11月的低流量頻率變異后高于天然徑流量的頻率。

3.3.2流量極值和延時變化

年極端流量中相對改變最大的是年1日最小流量，從變異前2 770 m3/s到變異后2 470 m3/s，減少了300 m3/s。年極端流量和出現時間的HMA法結果只有年1天最小流量和最小極端徑流時間改變度在33%和67%之間，為中度改變，其他均屬于低度改變。差異度最小的是基流指數，為17.55%，最大的是年1日最小流量發生時間，為39.44%。跟月平均流量一樣，各項指標的低流量頻率變異后高于天然徑流量的頻率，見圖3。

圖2 月平均流量序列變異前后頻率分布圖

3.3.3流量變化的頻率

HMA法計算得到的變異前后宜昌站流量平均上升率、下降率和逆轉次數改變度分別為51.68%、34.17%、63.80%，都屬于中度改變、低度。變異后流量逆轉次數改變較大，從逆轉次數頻率分布圖中(圖4)可以明顯看出變異后逆轉次數較大區間的頻率明顯增加，逆轉次數均大于66，大于101的頻率為0.80，而變異前大于101的頻率為0。變異之后上升率改變也較大，從圖4中可知變異前上升率小于260%的頻率值為0，在260%到338%之間的頻率為0.07，在338%到417%之間的頻率為0.08，在417%到495%之間的頻率為0.22，變異后，上升率小于260%的頻率為0.23，在260%到338%之間的頻率為0.39，在338%到417%之間的頻率為0.19，在417%到495%之間的頻率上升為0.32。上升率較小區間的頻率明顯增加。

圖3 流量極值和延時變化變異前后頻率分布圖

圖4 流量變化的頻率變異前后頻率分布圖

4 結 語

(1)通過宜昌站水文情勢趨勢分析，年平均凈流量在1970年左右發生變異，這與上游的葛洲壩工程和黑龍灘水庫等水利設施的建成時間較為符合。

(2)宜昌站32個水文變化指標大部分指標變化趨勢較明顯，其中月平均流量指標普遍減少；年極端流量普遍減小，年低流量發生次數也增加，低流量平均延時同樣增加；流量平均上升率和流量逆轉次數顯著增加。

(3)宜昌站徑流趨勢的減少對生態環境問題，例如水質的變化、水生生物的生存等都有重大影響，同時也會影響到航運、下游人民用水、輸沙、地下水位等一系列問題，如何保證所有方面的問題都考慮周全，需要水資源管理部門的保護和調控。

□

[1] 藺秋生, 范北林, 黃 莉. 宜昌水文站年徑流量演變多時間尺度分析[J].長江科學院院報, 2009,26(4):1-3.

[2] Aziz O I A, Burn D H. Trends and variability in the hydrological regime of the Mackenzie River Basin[J]. Journal of Hydrology,2006,319(1-4):282-294.

[3] Kendall MG. Rank Correlation Methods.Griffin[M]:London, UK, 1975.

[4] 陳 華,閆寶偉,郭生練,等. 漢江流域徑流時空變化趨勢分析[J]. 南水北調與水利科技,2008,(3):49-53.

[5] Richter B D, Baumgartner J V, Powell J, et al. A method for assessing hydrologic alteration within ecosystems[J]. Conservation Biology,1996,10(4):1 163-1 174.

[6] Richter B D, Baumgartner J V, Braun D P, et al. A spatial assessment of hydrologic alteration within a river network[J]. Regulated Rivers: Research & Management,1998,14(4):329-340.

[7] 張 鑫, 丁志宏, 謝國權, 等. 水庫運用對河流水文情勢影響的IHA法評價——以伊河陸渾水庫為例[J]. 水利與建筑工程學報,2012,02:79-83.

[8] Poff N L, J D Allan, M B Bain,et al. The natural flow regime: A paradigm for river conservation and restoration[J]. BioScience,1997,47(11):769-784.

[9] 陳棟為, 陳曉宏, 李 翀, 等. 基于RAV 法的水利工程對河流水文情勢改變的累積效應研究——以東江流域為例[J].水文, 2011,31(2):54-57.

[10] Shiau J T, Wu F C. A histogram matching approach for assessment of flow regime alteration: application to environmental flow optimization[J]. River Research and Applications 24: published online, 2008,24:914-928.