基于SHAW模型的山西沁河流域水資源配置方法研究

趙電雷

(晉城市水務(wù)局水利水電事務(wù)中心，山西 晉城 048000)

0 前 言

目前常用的水資源配置方法有基于改進果蠅算法[4]和基于地下水“雙控”2種水資源配置方法[5]。但上述方法計算得到流域產(chǎn)流損失數(shù)值較小，易導(dǎo)致水資源出現(xiàn)供需失衡的情況。SHAW模型最初用于模擬土壤凍融過程，研究凍融過程水分的遷移，經(jīng)過研究人員不斷地改進，SHAW模型可實現(xiàn)對循環(huán)系統(tǒng)中多種物理指標的模擬，并整合多指標的物理性質(zhì)形成垂向的一維剖面，在該剖面內(nèi)建立一個生物分層系統(tǒng)，描述剖面的水通量，從而能夠更加準確得到流域周圍土壤成分中水分的變化。

山西沁河流域發(fā)源于長治市沁源縣霍山南麓的二郎神溝，屬于黃河的一級支流，位于太行山脈西側(cè)，地形北高南低，大部分流域為海拔700 m左右及以上的海拔高度的山區(qū)，支流數(shù)量眾多[1]。由于地勢原因，沁河的徑流大多由降水形成，降水一部分形成地表水匯入河道中，另外一部分滲入地表下形成地下水[2-3]。為此，本文研究利用SHAW模型，結(jié)合沁河流域的自然環(huán)境，提出1種山西沁河流域水資源配置方法，因地制宜的調(diào)度沁河流域的水資源，合理配置水資源，改善沁河的水生態(tài)環(huán)境。

1 水資源配置方法研究

1.1 SHAW模型處理流域下邊界水熱通量

由于沁河流域的特殊地理位置及特殊的土壤剖面，故采用SHAW模型處理流域下邊界水熱通量。考慮到沁河流域海拔過高，存在部分土壤凍結(jié)融化的過程，首先計算該過程中土壤內(nèi)的水通量，計算過程如下：

(1)

式中：θl為土壤體積內(nèi)含水率變化；θi為土壤體積內(nèi)含冰率變化；t為時間；ρi表示進入土層的液態(tài)水通量；ρl為進入土層的氣體通量；z為土壤深度；h為土壤基質(zhì)勢；k為土壤導(dǎo)水率；qv為土壤內(nèi)水分蒸發(fā)所需的熱量；U為根系吸水源系數(shù)。在此基礎(chǔ)上，根據(jù)水通量計算結(jié)果，采用BC模型描述流域周圍土壤水分特征，描述公式可表示為：

(2)

式中：he為進氣值;θsat為孔隙含水量;b為土壤孔徑分布參數(shù)。當沁河流域土壤凍結(jié)的情況下，假定該部分凍土水分位于非飽和土壤中，土壤孔隙的基質(zhì)勢不變，導(dǎo)水率受到冰量的影響，導(dǎo)水率逐漸減小直至為零，此時土壤內(nèi)的有效孔隙率為0.13[6]。但土壤孔隙中存在部分的溶質(zhì)勢，孔隙內(nèi)溫度在凍結(jié)溫度之下，土壤孔隙內(nèi)部的液態(tài)水與冰處于平衡狀態(tài)，故得到土壤總土水勢與溫度間的關(guān)系，如下：

(3)

式中：E為土壤總土水勢;R為溶質(zhì)勢，Lf為溶質(zhì)水的凍結(jié)潛熱，g為重力加速度，Tk為冰量轉(zhuǎn)化為溶質(zhì)的時間，其余參數(shù)含義保持不變。根據(jù)上述數(shù)值結(jié)果，計算得到溶質(zhì)勢數(shù)量關(guān)系：

(4)

式中：c為土壤水質(zhì)中溶質(zhì)濃度;A為氣體常數(shù)。根據(jù)上述得到的土壤液態(tài)水含量，假設(shè)沁河流域土壤內(nèi)大于含水量部分的水分以冰形態(tài)存在[7]，考慮到冰在融化時吸收的熱量，沁水流域周圍的溫度分布就可表示為：

(5)

式中：Cs為體積熱容量;T為水域土壤溫度;ql為液態(tài)水通量;cl為水的比熱容，其余參數(shù)含義不變。部分土壤中水溶質(zhì)存在部分的植物吸收，計算溶質(zhì)在流域土壤內(nèi)的吸收系數(shù)，計算過程如下：

(6)

式中：Dw為溶質(zhì)的濃度;τ為土壤迂曲度參數(shù)。綜合上述SHAW模型處理過程，流域土壤含水量及土壤溫度可根據(jù)不同沁河流域的季節(jié)改變，模型可自行對流域土壤含水量及溫度進行線性處理，得到不同季節(jié)時間點下沁水流域下邊界熱量變化[8]。

根據(jù)下邊界數(shù)值結(jié)果，采用SHAW模型模擬得到整條沁水流域含有水力梯度的重力梯度，進而得到沁水流域土壤底部溫度[9]。根據(jù)溫度變化值，得到沁水流域外部降水部分轉(zhuǎn)化為地下水的含量，便于合理化配置沁水流域水資源。

1.2 山西沁水流域可配置水資源調(diào)度量計算

(7)

(8)

式中：Kp為模比系數(shù);Hp為點暴雨均值。計算得到定點雨量的計算結(jié)果如表1所示。

參照上表所得到的各項數(shù)值，得到面暴雨的計算公式，如下：

(9)

式中：Hp,A(tb)為設(shè)計面暴雨；ηp(A,tb)為點面折減系數(shù)，可按下式計算：

(10)

式中：A為流域面積;C，N均為經(jīng)驗參數(shù)。整理得到設(shè)計面的暴雨計算公式：

表1 設(shè)計定點雨量計算結(jié)果表

(11)

式中：Sp為設(shè)計雨力，即1 h的設(shè)計雨量;λ為經(jīng)驗參數(shù);ns為設(shè)計面雨量修正值。在此基礎(chǔ)上，計算山西沁河流域產(chǎn)流，計算過程如下：

(12)

式中：tanh為雙曲正切運算符;tZ為設(shè)計暴雨的主雨歷時;Rp為設(shè)計洪水凈雨深;FA(tZ)為主雨歷時內(nèi)的流域可能損失。根據(jù)該部分流域的可能流失，得到該時段下凈雨量，計算公式如下：

Δhp,j=hp(tj-1+Δt)-hp(j-1)

(13)

式中：hp(t)為設(shè)計時段凈雨深;Δt為雨量計算時段;j為雨量流失時間段。綜合沁水流域內(nèi)的調(diào)度量作為總配置量[11]，實現(xiàn)山西沁河流域水資源的的合理配置。

1.3 實現(xiàn)水資源配置

應(yīng)用上述得到的可配置調(diào)度量，采用綜合瞬時單位線法進行匯流計算，量算各匯流地類面積占流域面積的權(quán)重，根據(jù)權(quán)重值計算沁水流域洪水過程線，過程線計算公式如下：

(14)

式中：J為河道比降;C1,A為復(fù)合地類匯流參數(shù);C1,i為單地類匯流參數(shù);β為經(jīng)驗性指數(shù);α為面積權(quán)重。根據(jù)上述計算出的數(shù)值，得到書院河段的洪峰流量及最大24 h洪量值，結(jié)果如表2所示。

表2 洪峰流量及最大24 h洪量值表

以沁水流域洪水過程線內(nèi)的水資源作為分配指標，預(yù)測山西省對沁水流域的需水情況時，以沁河流域人口分布情況作為配置依據(jù)[12]，見表3所示。

表3 沁河流域現(xiàn)狀人口分布情況表

按照表3所示的人口成分，綜合區(qū)域內(nèi)城鎮(zhèn)化的情況，將人均用水量定額為200 L/人·d，參照城鎮(zhèn)與農(nóng)村人口的數(shù)量預(yù)測得到沁河流域配置的生活需水量[13]。利用表3中農(nóng)村人口數(shù)量，認定當前山西沁河流域周圍的耕地面積不變。根據(jù)沁河流域產(chǎn)業(yè)分布情況，確定3項產(chǎn)業(yè)的水資源需求以及產(chǎn)業(yè)分布情況，結(jié)果如表4所示。

使用表4所示的產(chǎn)業(yè)分布，在考慮農(nóng)業(yè)節(jié)水灌溉技術(shù)水平的同時，利用第一產(chǎn)業(yè)占比數(shù)值，查閱歷史資料的農(nóng)業(yè)用水來配置沁河的水資源[14]。由于工業(yè)、建筑業(yè)需水受到行業(yè)結(jié)構(gòu)、生產(chǎn)技術(shù)等外部因素的影響，該行業(yè)配置水資源存在較大不確定性，故將一個工作單位作為計算單位，每個單位分配可用于一周日常工作活動的預(yù)算水量[15]。同理配置得到第三產(chǎn)業(yè)配置得到的水資源。

表4 沁河流域產(chǎn)業(yè)分布表

2 仿真實驗與結(jié)果分析

為驗證基于SHAW模型的山西沁河流域水資源配置方法的應(yīng)用效果，采用仿真實驗加以驗證。

2.1 實驗準備

實驗在MATLAB仿真平臺中完成。選用沁河流經(jīng)晉城市內(nèi)的書院河作為實驗對象。書院河道總長為4.61 km，流域面積5.95 km2，平均縱坡26.8‰。書院河的地理位置如圖1所示。

圖1 實驗流域?qū)ο蟮乩砦恢脠D

在此基礎(chǔ)上，結(jié)合書院河流域周圍的產(chǎn)業(yè)分布及人口分布，設(shè)定不同降水保證率下工業(yè)、農(nóng)業(yè)及生活用水情況，結(jié)果如表5所示。

表5 實驗設(shè)定的水資源使用情景表

按照表5所示的設(shè)定情景，分別采用傳統(tǒng)的基于改進果蠅算法、基于地下水“雙控”與基于SHAW模型的水資源配置方法完成水資源模擬配置，并對比3種配置方法的應(yīng)用性能。

2.2 實驗結(jié)果及分析

根據(jù)上述實驗準備，以配水量為計算指標，定義水資源的供需平衡系數(shù)計算公式如下：

(15)

式中：ε為供需平衡系數(shù)；Q為設(shè)定情景內(nèi)的需水量；Q′為配備方法配置水量。一般來說，ε的數(shù)值越接近于1，則表示水資源配置方法的配比結(jié)果更佳。

匯總3種水資源配置方法的供需平衡系數(shù)，結(jié)果如表6所示。

表6 三種配備方法供需結(jié)果表

由表6所示的供需平衡系數(shù)可知，基于改進果蠅算法的方法得到的供需平衡系數(shù)小于1，表明該方法并超出書院河所有產(chǎn)業(yè)的用水需求的預(yù)期，配置的水量過多，存在水資源浪費的情況，不利于流域內(nèi)的水資源循環(huán)。基于地下水“雙控”的方法得到的供需平衡系數(shù)在1.20上下，大于1，表明該方法所配置的水量無法滿足該支流下產(chǎn)業(yè)的日常運營，不適合在沁河流域中使用。而應(yīng)用基于SHAW模型的水資源配置方法得到的供需平衡系數(shù)接近1，表明該方法能夠在滿足產(chǎn)業(yè)的日常用水需求外，有效避免了配置的水資源過量。

綜上所述，基于SHAW模型的水資源配置方法更加適合在沁水流域水資源配置中應(yīng)用，符合水資源供需平衡的實際需求。

3 結(jié) 語

隨著流域水資源與社會用水供需矛盾的日趨激化，對水資源進行合理分配成為了相關(guān)領(lǐng)域的研究熱點。本研究在SHAW模型的基礎(chǔ)上，實現(xiàn)了對山西沁河流域水資源的合理配置，有效減少了配置過程中水資源的浪費。