郎宏磊
(吉林省水文水資源局白城分局,吉林白城137000)
烏拉特中旗狼山地區(qū)地表水資源量評價研究
郎宏磊
(吉林省水文水資源局白城分局,吉林白城137000)
通過對烏拉特中旗狼山地區(qū)水文氣象資料收集和分析處理,對地表水資源量進行了評價計算。旨在說明其地表水資源數(shù)量和地表水資源可利用量,為實現(xiàn)該地區(qū)的水資源的供需平衡及水資源合理開發(fā)、利用、節(jié)約、保護提供依據(jù)、對策和建議。
水資源評價;烏拉特中旗狼山地區(qū);內蒙古
以內蒙古烏拉特中旗狼山地區(qū)為評價區(qū),總面積為9960.50km2,屬大陸性干旱氣候區(qū)。評價區(qū)北部山區(qū)屬于陰山山脈西段,包括前、后狼山各一部分。山脈走向近東西,主要是中低山、低山丘陵、沖狀態(tài)高平原和河谷平原地形。
烏拉特中旗以狼山山地丘陵與渣爾泰山余脈為分水嶺,劃成兩個水系。南部為黃河水系,從黃河最大年引水量為2.22億m3,多年平均引水量為1.76億m3;北部為內陸河水系。狼山地區(qū)內山地、丘陵,溝谷河系較發(fā)育,河床較寬,地形坡度在1/30~1/300之間,出口處形成較大的沖洪積扇。狼山地區(qū)內河流、水系不發(fā)育,在山丘區(qū)多以南北向干谷為主,僅布爾嘎斯溝、東烏蓋溝、呼勒斯太溝、前達門溝、罕烏拉溝、敖布拉格、木倫河、扎拉格河等常年有水,其它溝谷均屬季節(jié)性河流。
2.1數(shù)據(jù)收集
針對評價區(qū),共收集了狼山地區(qū)及其周邊地區(qū)3個水文站、12個降水量站和5個氣象站的水文氣象資料。影響地表水資源數(shù)量評價主要為降水、蒸發(fā)(E601數(shù)據(jù))和徑流3個因子,其它次要因子(氣溫、相對濕度、風速、日照時數(shù)等)對計算出的結果影響不大,在這里不詳細敘述。
2.2水資源計算分區(qū)
根據(jù)《水資源評價導則》及《全國水資源綜合規(guī)劃技術大綱》要求,水資源評價先以流域為基礎,再套合匯總到各行政分區(qū)上。因此,此次水資源評價按照河流水系流入同一個流域劃分為一個區(qū),共分了14個二級區(qū),在二級區(qū)基礎上,我們又劃分了18個三級區(qū)。
2.3數(shù)據(jù)的處理及初步分析
2.3.1降水資料
對于需要插補、延長的某一站點來說,在對參證站的降水資料代表性、一致性、連續(xù)性分析的基礎上,采用單變量或雙變量把其與參證站進行相關分析,進而建立相關方程,進行插補延展。針對插補的測站而言,在選擇用于插補的參證站時,盡量選擇鄰近的降水量測站,其次選擇緯度相近的降水量測站,最后選擇經度相近的降水量測站。當有多個測站可供選擇,無法確認哪一個作為參證站時,分別對插補站、各參證站的降水量進行單變量回歸分析,最后選擇相關最顯著的參證站。如果有多個測站的相關關系都較好,相關系數(shù)又相差不大,則均選作參證站。此次降水量回歸分析時,參證站選擇控制在2個以內。利用插補延展后的降水資料通過距離反比法計算流域各分區(qū)的平均降水量。
2.3.2徑流量資料
1)利用流域平均年降水量與年徑流量、測站年徑流量與年徑流量之間關系進行插補。在流域平均年降水量與年徑流量的相關分析中,只有金門相關性達到相關系數(shù)檢驗值要求,即r(金門)=0.801>R臨界(2,17,0.01)=0.606,回歸方程通過F檢驗,方程穩(wěn)定,可以用于回歸模擬。
2)將其余各水文站流域平均月降水量與月徑流量、測站月徑流量與月徑流量系列作相關分析,相關系數(shù)均接近相關系數(shù)檢驗值,但未達到檢驗值。因此,嘗試同時采用月降水和月徑流分段插補月徑流量資料。例如海流圖站月徑流量序列中,11—1月徑流量序列、2—4月徑流量序列與阿塔山對應月序列相關性較好,而5—10月徑流量序列與該站對應月降水量序列相關性較好,滿足相關檢驗臨界值要求,同時回歸方程通過F檢驗,方程穩(wěn)定,可以用于回歸模擬。將插補延展后得到各徑流站的年徑流量進行一致性分析,才能用于天然徑流量的還原計算。
2.3.3蒸發(fā)
作為水循環(huán)的重要環(huán)節(jié),蒸發(fā)是水量平衡的重要要素。用折算系數(shù)法計算水面蒸發(fā),公式為:
E0=KE601(1)式中:K為折算系數(shù),目前許多地區(qū)尚無實測的折算系數(shù),實際計算中往往是用臨近地區(qū)數(shù)值進行計算。E601為E601型蒸發(fā)器,mm。E0為蒸發(fā)量,mm。計算已有站點蒸發(fā)數(shù)據(jù)的干旱系數(shù),見表1。
表1 各時間段干旱系數(shù)
根據(jù)烏拉特中旗氣象站多年平均E601各月月均水面蒸發(fā)量的統(tǒng)計分析,蒸發(fā)量的月分配受氣候的月際變化影響較大,1月蒸發(fā)量最小,其值為15.55mm,蒸發(fā)量自2—3月開始迅速增大,到6月達最大,其值為273.94mm,6—10月緩慢下降,在10—11月間又迅速下降,而后才緩慢下降到1月的最小值。蒸發(fā)量多年平均最大連續(xù)4個月(5—8月)水面為968.4mm,占全年1546.5mm的62.6%。
陸面蒸發(fā)(土壤蒸發(fā)和植物散發(fā))難于直接觀測,采用水量平衡估算,即流域年降水量減去年徑流量得到流域蒸發(fā)量。計算可得評價區(qū)內多年平均陸面蒸發(fā):1956—1979年為184.9mm;1980— 2006年為185.4mm;1956—2006年為185.1mm。
3.1流域區(qū)域徑流量計算
用水文比擬法(降水、面積倍比法)計算評價區(qū)域徑流量。按已知流域的流域面積、流域平均年降水量和流域天然徑流量,換算出未知流域天然年徑流量,換算公式為:
式中:Wy為未知流域天然年徑流量,萬m3;Fy為未知流域產流區(qū)面積,km2;Fx為已知流域產流區(qū)面積,km2為未知流域平均年降水量,mm;為已知流域平均年降水量,mm;Wx為已知流域天然年徑流量,萬m3。
從3個角度對研究區(qū)域的徑流量進行了計算:1)將3個已知水文站按照地理位置,分配給相鄰的分區(qū),從而計算研究區(qū)域的水資源量;2)由于阿塔山距離所有分區(qū)較遠,在計算過程中,只選取海流圖和金門來計算研究區(qū)域水資源量;3)阿塔山的徑流資料最完整,它沒有進行插補延展,計算出的數(shù)據(jù)誤差較小,故只用該站計算研究區(qū)水資源量。
3.2結果與分析
阿塔山站計算的水資源量最大,用海流圖、金門和阿塔山3站計算量次之,而用海流圖和金門兩站計算量最?。挥嬎憬Y果顯示1981年之后水資源量有明顯的減少。
原因分析:1)由阿塔山計算出的水資源量相對偏大,因為阿塔山站地處研究區(qū)域東部,降水規(guī)律為自西向東降水量逐漸增加,因此該站測出的降水量、徑流量都相對較大,由它計算出的研究區(qū)水資源量也會隨之較大。綜合3個方法計算的水資源量,最后選擇由阿塔山、海流圖、金門3個水文站計算出的研究區(qū)水資源量。2)在計算兩時間段內的區(qū)域蒸發(fā)量時發(fā)現(xiàn),區(qū)域年蒸發(fā)量變化不大;而區(qū)域水資源量與降水量相關性密切,而年平均降雨量由 1981—2006年時段內的174.9 mm下降到1956—1980年的165.7mm,減少了5.3%,可見降雨量的減少可能是地表水資源量較少的主導因素。
3.3地表水資源可利用量估算模擬
鑒于海流圖站以上流域從各方面條件均可近似代表評價區(qū)的實際情況,因此,選擇海流圖站以上流域地表水資源可利用率作為評價區(qū)地表水資源可利用率。
海流圖站以上流域地表水資源可利用率計算過程:
1)計算海流圖站地表水資源量的多年平均值。據(jù)1956—2006年海流圖站天然年徑流量,統(tǒng)計出海流圖站地表水資源量的多年平均值為443.844萬m3。
2)計算維持河道基本功能的生態(tài)需水量。多年平均河道最小生態(tài)環(huán)境需水量可取年徑流量的10%~20%,根據(jù)該水系的實際情況,取多年平均年徑流量的15%作為平均河道最小生態(tài)環(huán)境需水量。其值為W生態(tài)=443.844×0.15=66.577萬m3。
3)計算地表水資源可利用量,其值為地表水資源量減去河道內最小生態(tài)環(huán)境需水量后,再減去汛期難于控制利用的洪水量。
海流圖站下游有德嶺山水庫,該水庫除控制海流圖河外,另還控制博格爾高勒溝和烏敦熱格高勒溝,兩溝匯合后又與海流圖河匯合進入德嶺山水庫。德嶺山水庫在除險加固后調節(jié)庫容為488萬m3,由于該水庫是該流域主要利用地表水的水源工程,因此,我們可采用其對海流圖河的一次最大調節(jié)水量作為控制汛期洪水下泄的水量。經過對另兩溝天然徑流量的求取,按照三溝同時發(fā)生洪水的原則,德嶺山水庫對單一海流圖河的一次最大調節(jié)水量為231.52萬m3。
由于洪水的隨機性以及評價區(qū)現(xiàn)狀地表水資源開發(fā)利用程度較低,尤其是流域上游地區(qū)。因此,流域汛期可以控制利用的最大洪水量要比231.52億m3還大??紤]全國水資源綜合評價規(guī)范要求,并結合當?shù)氐乇硭Y源開發(fā)利用現(xiàn)狀,可以控制利用的洪水量上限取實測值的1.1倍,即本地區(qū)為231.52×1.1=254.672萬m3。
根據(jù)以上確定的汛期可以控制利用的最大洪水量,應用海流圖站1956—2006年51年汛期洪水量(7—9月天然)系列,逐年計算汛期下泄洪水量。汛期洪水量中大于的部分作為難于控制利用的洪水量,汛期洪水量小于或等于,則下泄洪水量為0。根據(jù)算出的下泄洪水量系列,計算多年平均汛期難于控制利用的洪水量。
式中:W泄為多年平均汛期難于控制利用的洪水量,萬m3;W1d為海流圖站第i年汛期(7—9月)天然徑流量,萬m3;Wm為汛期可以控制利用的最大洪水量,萬m3;n為系列年數(shù)(51年)。
歷年汛期難于控制利用的洪水量51年的多年平均值為127.68萬m3。
4)計算地表水資源可利用率,其值為地表水資源可利用量除以地表水資源量。
根據(jù)以上計算的多年平均最小生態(tài)環(huán)境需水量66.58萬m3和汛期難于控制利用的洪水量127.68萬m3,計算得出海流圖站以上流域多年平均地表水資源可利用量為249.58萬m3。則得出地表水資源可利用率為56.23%。
研究區(qū)地表水資源可利用量估算,由公式W可= 0.56·W總直接計算可利用資源量。1956—2006年,狼山地區(qū)的地表水資源量為5599.20萬m3,可利用量為3135.55萬m3。
水資源作為一種動態(tài)資源,隨人為干預及氣候環(huán)境的變化而變化。以水文氣象數(shù)據(jù)為依據(jù),估算了烏拉特中旗黃河流域狼山地區(qū)的地表水資源量及可利用地表水資源。1956—2006年,狼山地區(qū)的地表水資源量為5599.20萬m3。可利用量為3135.55萬m3。
針對狼山地區(qū)內徑流泥沙數(shù)據(jù)較少的現(xiàn)狀,建議增加數(shù)據(jù)測站,為研究區(qū)內水資源的節(jié)約、開發(fā)利用、配置及治理工作提供足夠的數(shù)據(jù)支持。由于研究區(qū)內降雨量年內分配極其不均,汛期降水集中,加之水利工程少,易造成階段性水資源匱乏,建議增加水利工程。另外,由于未收集到最新數(shù)據(jù)資料,只能對已有數(shù)據(jù)年限水資源量研究分析,而近十年隨著地區(qū)政策施行及水利工程措施的進步,水資源量計算是會出現(xiàn)量值上的差異,而此研究可作為對比研究的依據(jù)。
X832
A
1002-0624(2015)10-0035-03
2015-07-15