摘要：基于信息熵理論，計算了大通水文站的流量歷時曲線。利用兩個簡單約束方程對流量歷時曲線進行建模，運用概率總和與平均流量以及約束方程推導最大熵獲得曲線參數，最終確定了流量歷時曲線。通過對大通控制站觀測到的流量數據進行驗證，發現基于信息熵理論得到的流量歷時曲線與觀測到的曲線一致。自2003年三峽工程建設投運以來，大通水文站觀測到的均值流量和最大流量趨勢發生了變化，卻未發現均值與最大值的比值存在顯著差異。在此期間，其表征分布形態的熵參數也并未發生顯著變化。

　　關鍵詞：流量歷時曲線;信息熵;熵參數;大通水文站

　　推薦閱讀：水利科技水利論文投稿期刊

　　1 研究背景

　　大通水文站位于安徽省池州市梅隴，為長江下游的控制水文站，距長江口約700 km，控制面積約170萬km2。該站于1922年10月設立，1949年10月1日后由長江水利委員會管轄，觀測資料包括水位、流量、泥沙、雨量等。

　　流量歷時曲線代表水文控制站不同頻率對應的流量由低到高排列，一般在某一時間尺度用0-100%的數值刻畫。流量歷時曲線可看作某一站點某一時段的流量累積頻率曲線，而非時間序列。由于流量歷時曲線是通過觀測范圍內所有觀測到的流量判定的，考慮到流域的綜合水文條件，可用來比對不同河流的情況，經常用來預測未來供水流量的分布[1]、水力發電量[2]、沉積物負荷[3]及污染指數[4]。

　　一般來說，可通過每日平均流量構建流量歷時曲線，也可根據需要建立周、月和季步長上的流量歷時曲線，但隨著時間間隔增加，流量變化的細節變得模糊。對于多數河流，月流量的歷時曲線并不能充分反映流量的變化情況，而基于年平均流量建立的流量歷時曲線由于其變化范圍較小，幾乎沒有價值。H.A. Foster[5]研究表明，基于日流量確定的流量歷時曲線與基于月流量確定曲線的差異可能高達35%。對于不同的河道，時間間隔變化的影響也不盡相同。對于一些日流量幾乎均勻的河道，周流量的歷時曲線可能與日流量的歷時曲線幾乎相同。而對于突然發生洪水的河道，其日流量和周流量的歷時曲線將存在較大差異。通常根據經驗分布構建流量歷時曲線。本文采用2014年V.P.Singh等[6]研發的信息熵模型，結合最大熵理論，確定大通水文站流量歷時曲線，并與其觀測值進行比較驗證。

　　2 研究方法

　　為利用信息熵理論計算流量歷時曲線，需假設大通站觀測到的日流量Q為獨立隨機變量，服從未知概率分布f(Q)。

　　V.P.Singh等[6]基于信息熵理論的流量歷時曲線計算步驟包括：①定義信息熵;②確定流量約束條件;③利用拉格朗日方程求最大熵并獲得流量概率分布;④計算拉格朗日參數并求得流量歷時曲線。

　　2.1 定義流量的Shannon熵

　　考慮觀測到的流量Q服從未知分布f(Q)，利用Shannon熵[7]，流量的信息熵H(Q)可定義為式中，Qmax和Qmax分別為上下限，表征最大及最小流量。H表征流量概率分布f(Q)所包含的不確定性，代表樣本Q所包含信息量的期望。通過式(1)，根據最大熵理論求得最不偏的概率分布f(Q)。最大熵理論由E.T.Jaynes[8]的研究成果發展而來，只需明確流量約束條件即可求得。根據最大熵理論可知，對僅知部分信息的變量或關系，對其最不偏的判斷滿足熵最大的情況。

　　2.2 約束條件

　　首先，作為Q的概率分布f(Q)必須滿足其總和等于1的條件：

　　能量守恒定律也能作為約束條件，然而D.E.Barbe等[9]的研究表明，對于河道流速的分布，利用信息熵理論時僅需上述兩個約束條件就已足夠。為不增加計算的復雜性，該研究中僅采用上述兩個約束條件進行計算。

　　2.3最 大化

　　利用拉格朗日算子，可以將式(1)所定義的信息熵最大化，并求得熵最大時的概率分布：

　　2.4 流量歷時曲線

　　流量歷時曲線代表了流量和不小于該流量時長之間的關系，因此可假設觀測到的時均流量為在其最大值Qmax與最小值Qmin之間的概率分布。這符合實際情況和Laplacian不充分原則。為此，流量的累積概率F(Q)可表示為

　　F(Q)=1-t/T=1-τ，τ=t/T

　　(6)式中，t為流量不小于Q的時長(天數)，T為考慮的總時長(365 d)。值得注意的是，F(Q)為關于流量Q的函數，t/T是關于τ的函數。上述線性累積分布假設是對流量Q與時長百分比τ的線性假設，相似的假設已應用于不少研究[10-12]。相應地，也可假設流量與時長百分比間具有非線性關系。為簡化計算步驟，本文采用上述線性假設關系。

　　參照文獻[6]，對通過最大化信息熵獲得的概率分布，在Qmax和Qmax之間進行積分，即可得到流量歷時曲線：

　　3 計算結果

　　3.1 大通站日流量

　　采用大通水文控制站1951-2013年日均流量進行計算，以2003年三峽工程建成投運為界，分為2003年前后兩個序列。表1總結了2003年前后大通站年均流量、年最大流量、年最小流量的歷史平均值(Qmax，Qmax，Qmax)及其觀測到的歷史最大流量Qpeak。由表1可以看出，三峽工程建設投運以來，年均流量變化幅度不顯著，年最大流量減少，年最小流量增加，其流量變化范圍減小，然而該變化并不顯著。

　　為更準確地表明2003年前后的流量變化，圖1顯示了2003年前后兩組觀測流量的概率分布。由圖1可以看出，2003年之后的流量概率分布相比2003年前更為集中，流量在12 000-20 000m3/s之間的情況明顯增加，無流量小于9 000m3/s或大于65 000 m3/S的觀測值。對2003年前后的概率進行配線并對比，發現兩組序列所配的概率分布存在明顯差異，通過t檢驗得知2003年前后的流量不服從同一分布。