譚光超, 李姍姍, 李智民, 柳曉晨, 李 業(yè)
(湖北省地質局 水文地質工程地質大隊,湖北 荊州 434020)
基流是補給河水的地下徑流,是河流中比較穩(wěn)定的徑流部分。通常認為地下徑流是枯水季節(jié)河流的基本流量,故稱其為基流[1]?;鞯难芯砍晒麖V泛地應用于流域水文模擬、水資源評價與管理、生態(tài)環(huán)境保護等多個學科領域。
國內外學者大量的科學研究及實際應用表明,數值模擬法計算的基流值滿足基流所應具備的特征,具有客觀、可重復性、易操作等特點[2-4],成為目前基流分割中最常采用的方法。數值模擬基流分割法包括數字濾波法、HYSEP法和BFI法,由于不同研究區(qū)的產匯流條件存在較大的差異,各種基流分割方法在不同研究區(qū)的穩(wěn)定性、適應性并無絕對性。
隨縣位于長江中下游北部,屬于長江流域二級地下水流系統(tǒng),該區(qū)地貌形態(tài)包括構造侵蝕剝蝕低山、構造侵蝕剝蝕丘陵和弱侵蝕堆積崗波狀平原,屬于典型的丘陵山區(qū)地貌區(qū)[5],河川基流量是丘陵山區(qū)地下水的主要排泄項,在河川徑流中占有一定比例。隨縣水文地質研究程度相對較低,1∶50萬湖北省地下水資源調查評價報告采用徑流模數法和入滲系數法評價了多年平均天然資源量,1∶20萬隨縣幅水文地質調查未開展地下水資源評價工作,探索適合于鄂北丘陵山區(qū)隨縣各河流的基流分割方法,是研究該區(qū)河川基流量及其動態(tài)變化的基礎,對于提高該區(qū)地下水資源評價精度,加強地下水資源管理具有重要意義。
以設立于隨縣境內的各支流水文站實測逐日徑流量資料為基礎,采用數字濾波法、HYSEP法和BFI法等基流分割法對平水年、枯水年、特枯年3個代表年的實測數據進行基流分割,選取計算結果最合理、最穩(wěn)定的方法,為后續(xù)開展該區(qū)域地下水資源評價工作奠定基礎,并為水文地質條件類似的流域河川基流量動態(tài)變化研究提供參考。
隨縣位于長江中下游北部,屬北亞熱帶季風氣候,四季分明,春秋季節(jié)短、夏冬季節(jié)長,年平均氣溫15.6 ℃,多年平均年降雨量960 mm,降雨集中在5—8月,約占年降水量的70%,年蒸發(fā)量1 450~1 520 mm。流經隨縣的長江支流有府河、漢水、淮河,其中隨縣境內府河流域面積為4 539.85 km2,占縣域面積的80.03%,漢水流域面積437.98 km2,占縣域面積的7.72%,淮河流域面積695.17 km2,占縣域面積的12.25%。
根據隨縣境內流域分布及水文站的布設情況,隨縣境內漢水、淮河流域支流均未布設水文站,府河流域的涢水源、均水、涢水、漂水均布設了水文站,并且4個水文站均有多年實測逐日徑流量資料。考慮到府河流域占據隨縣縣域面積的80.03%,且隨縣境內府河流域與漢水及淮河流域水文地質條件、產匯流條件類似,地下水類型以變質巖裂隙水為主,主要賦存于風化裂隙和構造裂隙中,局部富水帶多沿斷裂構造和侵入巖與變質巖接觸帶展布,地下水排泄以分散徑流型為主,地下水天然動態(tài)類型均屬于入滲—徑流型[注]李智民、卞學軍、譚光超等,鄂北丘陵山區(qū)嚴重缺水地區(qū)隨縣、應城市(北部)水文地質調查成果報告,湖北省地質環(huán)境總站,2016。。因此將府河流域各水文站實測數據作為本次研究的數據基礎具有相當程度的代表性,其研究結果可以類比用于整個隨縣。各水文站基本情況如表1所示。
表1 水文站基本情況Table 1 Basic information of the hydrological station
數值模擬基流分割方法包括數字濾波法、HYSEP法和BFI法,根據計算公式的不同或參數選取的方法不同,這3類數值模擬基流分割法可劃分為9種。
數字濾波法由學者Nathan和McMahon[6]于1990年引入到水文計算中,并由Chapman、Maxwell[7]進行了多次改進,該方法源于信號分析,其主要原理是將日徑流量視為高頻信號的地表徑流和低頻信號的河川基流組成,利用信號處理技術,將高頻信號和低頻信號分離,從而從逐日徑流量過程中分割出基流[8-9],該方法利用長歷時逐日徑流量資料計算,具有較好的客觀性和可重復性。根據基流分割方程的不同,可以分為F1、F2、F3、F4四種。
F1計算方法的基流分割方程為:
Qdt=Qd(t-1)·α+[Qt-Q(t-1)]·(1+α)/2
Qbt=Qt-Qdt
F2計算方法的基流分割方程為:
Qdt=Qd(t-1)·(3α-1)/(3-α)
+[Qt-α·Q(t-1)]·2/(3-α)
Qbt=Qt-Qdt
F3計算方法的基流分割方程為:
Qbt=Qb(t-1)·β/(2-β)+Qt·(1-β)/(2-β)
F4計算方法的基流分割方程為:
Qbt= Qb(t-1)·k/(1+k)+Qt·k/(1+k)
式中:Qdt和Qd(t-1)分別為第t時刻和第t-1時刻的地表徑流;Qt和Q(t-1)分別為第t時刻和第t-1時刻的河川徑流總量;Qbt為第t時刻的基流;α為濾波參數,一般取值0.95;β為退水參數,一般取值0.95;k為固定值,一般取值0.15。
HYSEP法由PettyJohn 和Henning[10]于1979 年提出,是美國地質調查局推薦使用的方法。原理是利用經驗公式計算直接徑流的持續(xù)時間,選定相應的時間間隔,并在此基礎上計算基流,根據選取的時間間隔內基流取值的選取方法不同,HYSEP法分為固定步長法(Fixed)H1、滑動步長法(Slide)H2、局部最小值法(Lealman)H3,3種方法的經驗公式為:
N=A0.2
式中:A為流域面積;N為直接徑流的持續(xù)時間,時間間隔通常3—11 d,選擇與2N最接近的奇數作為時間間隔,在此基礎上進行基流計算。
BFI法由英國水文所學者Wels[11]于1980年提出,是以基流指數為權系數計算基流量的方法。根據拐點檢驗因子f和k經驗值取值不同,BFI法分為標準BFI(F)法B1、改進BFI(K)法B2。該方法基于“簡單平滑分割法”,其計算原理是將每年以N天為一時段劃分成365/N個時段,確定每一時段內的最小流量值,如果某時段內的最小流量值與拐點檢驗因子的積小于相鄰左右時段內的最小流量值,則確定其為拐點,將所有拐點直線連接,即可得到基流過程線,過程線以下的面積即為該年的基流量。
根據隨縣雨量站1957—2015年的逐年降雨量實測數據,采用水文頻率計算適線法,選取頻率P=50%、P=75%、P=95%對應的降雨量作為平水年、枯水年、特枯年的設計值,并按照“年內分配最不利”原則確定2002年、2004年、2011年為平、枯、特枯年。
2002年4個水文站實測逐日徑流量的基流分割結果如圖1-圖4所示。
圖1 2002年澴潭站基流分割結果Fig.1 Base flow segmentation results of Huantan Station in 2002
圖2 2002年均川站基流分割結果Fig.2 Base flow segmentation results of Junchuan Station in 2002
圖3 2002年隨州站基流分割結果Fig.3 Base flow segmentation results of Suizhou Station in 2002
圖4 2002年萬店站基流分割結果Fig.4 Base flow segmentation results of Wandian Station in 2002
圖5 2004年澴潭站基流分割結果Fig.5 Base flow segmentation results of Huantan Station in 2004
圖6 2004年均川站基流分割結果Fig.6 Base flow segmentation results of Junchuan Station in 2004
2004年各水文站實測逐日徑流量的基流分割結果如圖5-圖8所示。
圖7 2004年隨州站基流分割結果Fig.7 Base flow segmentation results of Suizhou Station in 2004
圖8 2004年萬店站基流分割結果Fig.8 Base flow segmentation results of Wandian Station in 2004
2011年各水文站實測逐日徑流量的基流分割結果如圖9-圖12所示。
圖9 2011年澴潭站基流分割結果Fig.9 Base flow segmentation results of Huantan Station in 2011
圖10 2011年均川站基流分割結果Fig.10 Base flow segmentation results of Junchuan Station in 2011
圖11 2011年隨州站基流分割結果Fig.11 Base flow segmentation results of Suizhou Station in 2011
圖12 2011年萬店站基流分割結果Fig.12 Base flow segmentation results of Wandian Station in 2011
從圖1-圖12中可以看出,各基流分割方法的計算結果存在較為明顯的差異,主要體現(xiàn)為三個方面:
(1) 不同方法之間的結果差異較大,9種方法均可實現(xiàn)基流過程的分割,得到的基流過程趨勢基本一致,但是9種基流過程曲線重合程度較低、差異較大。
(2) 相同方法、不同月份數據得到的結果差異較大,在降雨量少、河川徑流量較小且動態(tài)變化幅度較小的1月初—4日底、9月初—12月底的兩個階段,基流過程線與徑流量過程線均呈現(xiàn)較為平緩的曲線,曲線重合程度較高;在降雨量多而集中、河川徑流量較大的5月初—8月底,基流、徑流曲線整體動態(tài)變化特征基本一致,但是基流過程線始終位于徑流過程曲線下方,能夠較好的反映出流域產流的阻尼效應。
(3) 相同月份數據、不同方法得到的結果差異較大,濾波分析法的F1、F2、F3基流過程曲線整體較為平緩,對于不同年份或不同月份的基流指數值辨識能力均較弱,HYSEP法H1、H2、H3基流過程曲線拐點較多,曲線整體呈鋸齒狀且與河川徑流量過程曲線呈現(xiàn)高度的同步性,峰值出現(xiàn)的時間點高度重合,未反映出流域產匯流的遲滯效應;濾波分析法F4、BFI法B1和B2基流過程曲線整體較為平滑,反映出了徑流的漲水、退水過程,且體現(xiàn)了徑流過程的遲滯效應,基流分割計算結果較為合理。
基流指數是基流量與年徑流量的比值,依據收集的實測數據,利用9種數值模擬基流分割法得到的各河流不同代表年的基流指數計算結果見表2。
從表2可以看出,9種數值模擬基流分割方法的計算結果差異巨大,基流指數結果最小值F3為0.006 2、最大值H1為0.630 6,差距約100倍;基流指數結果平均值最小值F2為0.017 1、平均值最大值H2為0.476 6,差距約28倍。數字濾波法的基流指數值計算結果較分散,平均值最小值0.017 1、平均值最大值0.339 4,尤其是F2、F3方法計算基流指數結果不足0.02,即基流占比不足河川徑流量的2%,計算結果過于保守,不符合鄂北丘陵山區(qū)地下水排泄特征;HYSEP法的計算結果差異不大,基流指數值計算結果相對較集中,平均值在0.4~0.5;BFI法的計算結果差異也不大,基流指數值計算結果相對較集中,平均值在0.25~0.3。
表2 基流指數計算結果Table 2 Base flow index calculation result
上述基流分割過程對比結果得知濾波分析法F4、BFI法B1和B2計算結果較為合理,而基流指數計算結果對比得知濾波分析法F2、F3方法不適用與本研究區(qū),其余方法的合理性有待進一步對比研究。為綜合對比分析不同基流分割方法計算結果的穩(wěn)定性,對9種數值模擬基流分割法得到的計算結果進一步統(tǒng)計分析,結果如表3。
從表3可以看出,除去濾波分析法F2、F3以外,濾波分析法F4基流指數統(tǒng)計值極值比最小(1.961 1)、方差最小(0.003 5),不同水文站年際基流指數值在0.20~0.42,反映出該方法計算結果年際變化小、基流成果最穩(wěn)定,說明濾波分析法F4是這9種方法中最為可靠、穩(wěn)定的基流分割方法。
本文以設立于隨縣境內的涢水源、均水、涢水、漂水4個水文站實測逐日徑流量資料為數據基礎,采用3類9種數值模擬基流分割法對平水年、枯水年、特枯年3個代表年的實測數據進行基流分割,得到如下結論:
相同研究區(qū)、相同數據,9種基流分割方法的計算結果差異巨大。通過基流分割過程對比,濾波分析法F4、BFI法B1和B2基流過程曲線最為平滑,較好的反映出了降雨—徑流過程的阻尼效應和遲滯效應,基流分割過程曲線較為合理;通過進一步分析,根據基流指數計算結果對比、基流指數結果統(tǒng)計值對比,濾波分析法F4基流指數統(tǒng)計值極值比最小(1.961 1)、方差最小(0.003 5),不同水文站年際基流指數值在0.20~0.42,反映出該方法計算結果年際變化最小、基流成果最穩(wěn)定。
表3 基流指數值計算結果的統(tǒng)計值Table 3 Statistical value of the calculation result of the base flow index value
綜上所述,依據基流過程的年際變化和年內變化特征,濾波分析法F4是這9種方法中最為可靠、穩(wěn)定的基流分割方法??紤]到本次研究的數據基礎具有相當程度的代表性,其研究結果可以類比用于整個隨縣,因此濾波分析法F4可以適用于鄂北丘陵山區(qū)隨縣境內河流的基流分割,從而為后續(xù)隨縣地下水資源評價工作提供依據,并可以為類似的變質巖區(qū)河川基流量及其動態(tài)變化的研究提供重要參考。結合其他學者研究成果,充分表明不同基流分割方法在不同地區(qū)表現(xiàn)的適宜性不同,適用其他地區(qū)的具體基流分割方法尚需開展專項研究。