基于改進(jìn)新安江模型的巖溶地區(qū)徑流過程模擬
——以廟溝巖溶流域?yàn)槔?/h1>
2020-06-29 05:28李玉坤燕子琪王紀(jì)元
安全與環(huán)境工程 2020年3期
李玉坤，燕子琪，王紀(jì)元，陳 麗，周 宏
(中國地質(zhì)大學(xué)(武漢)地質(zhì)調(diào)查研究院，湖北 武漢 430074)
流域水文模型是在流域產(chǎn)匯流結(jié)構(gòu)概化的基礎(chǔ)上，以一個實(shí)體或數(shù)學(xué)模型來模擬流域降雨-徑流的形成過程，是一種常用的徑流過程模擬方法。流域水文模型通常分為集總式水文模型如新安江模型、SAC模型和水箱模型，分布式水文模型如SWAT模型、SHE模型，半分布式水文模型如TOPMODEL模型等。
巖溶地區(qū)含水介質(zhì)內(nèi)部結(jié)構(gòu)復(fù)雜多樣，實(shí)測資料獲取難度較大，使得流域水文模型在該地區(qū)的應(yīng)用成為難點(diǎn)。分布式水文模型基于水動力學(xué)機(jī)制來描述和模擬流域內(nèi)的水文過程，但對資料要求較高，因而在缺乏資料地區(qū)的應(yīng)用受到限制；黑箱模型又稱經(jīng)驗(yàn)?zāi)Ｐ?，?yīng)用門檻低、適用性強(qiáng)，但模型不能反映含水系統(tǒng)內(nèi)部的水文過程，也無法用于評價(jià)和管理人類活動對泉流量的影響[1]。新安江模型是一種概念模型，由于采用張力水和自由水蓄水容量曲線來描述張力水和自由水蓄水容量分布的空間異質(zhì)性，因此已成為一種具有分布意義的概念性流域水文模型。近年來，隨著GIS和RS技術(shù)以及計(jì)算機(jī)的發(fā)展，新安江模型參數(shù)的精度和計(jì)算分區(qū)的精度已得到顯著提升，如李致家等[2]將新安江模型應(yīng)用于柵格之上，與高精度雷達(dá)降雨資料耦合進(jìn)行了洪水預(yù)報(bào)；井立陽等[3]根據(jù)地形因子分布曲線推求了自由水蓄水容量;Yuan等[4]應(yīng)用雙源蒸散發(fā)模型改進(jìn)了新安江模型結(jié)構(gòu)，模擬了在植被參與的水循環(huán)過程中的徑流過程；Zhao等[5]將新安江模型應(yīng)用于包括土壤雨水侵蝕、土地資源利用等土壤資源評價(jià)中。目前新安江模型已由功能單一的洪水預(yù)報(bào)模型發(fā)展為適用于多領(lǐng)域、多功能的流域水文模型。
在巖溶地區(qū)，新安江模型的應(yīng)用也得到了發(fā)展。如莊一鸰等[6]在總結(jié)南方巖溶地區(qū)水文特性的基礎(chǔ)上，對新安江模型的穩(wěn)定下滲、水源劃分和匯流模塊進(jìn)行了改進(jìn)，并取得了較好的模擬效果；郝慶慶等[7]在新安江模型原有自由水庫的基礎(chǔ)上，再增加一個線性水庫用來模擬深層地下水的出流過程，使模擬精度得到了顯著提高。但目前新安江模型在巖溶地區(qū)的應(yīng)用研究仍較為有限。為此，本文選取廟溝巖溶流域?yàn)榈湫蛶r溶流域，以降雨入滲為切入點(diǎn)，對新安江模型進(jìn)行了改進(jìn)，以期提高新安江模型對巖溶地區(qū)徑流過程模擬的精度。
1 研究區(qū)概況
廟溝巖溶流域位于湖北省宜昌市興山縣黃糧鎮(zhèn)高陽鎮(zhèn)境內(nèi),為高嵐河流域的子流域，整體上處于中國大陸二級階地向三級階地過渡地帶。該流域面積約為25 km2，平均海拔為1 100 m，地貌類型為溶蝕侵蝕中山溝谷地貌。該流域西側(cè)為斜坡，巖性以志留系粉砂巖為主；東側(cè)為巖溶洼地地貌，巖性以寒武-奧陶系碳酸鹽巖為主，區(qū)內(nèi)落水洞、巖溶漏斗等巖溶形態(tài)十分發(fā)育,見圖1。
圖1 廟溝巖溶流域水系概況圖Fig.1 Drainage map of Miaogou karst basin
2 研究方法
2.1 新安江模型及預(yù)模擬
新安江模型以蓄滿產(chǎn)流為理論基礎(chǔ)，該模型結(jié)構(gòu)為分散性結(jié)構(gòu)。如圖2所示，降雨信號輸入模型后轉(zhuǎn)化為張力水，此時(shí)一部分水量耗于蒸散發(fā)，采用三層蒸發(fā)模型進(jìn)行蒸散發(fā)計(jì)算，另一部分水量轉(zhuǎn)化為自由水，并將水源劃分為地表徑流、壤中流和地下徑流三部分；最后采用馬斯京根法進(jìn)行河道匯流計(jì)算。
圖2 新安江模型結(jié)構(gòu)簡圖Fig.2 Structure diagram of Xin’anjiang model
為了驗(yàn)證原始新安江模型對巖溶流域徑流過程的模擬效果，采用原始新安江模型對廟溝巖溶流域進(jìn)行了一次預(yù)模擬，選取的時(shí)間序列為2014年3月2日到2014年9月2日，期間共有4次小洪峰，1次大洪峰。原始新安江模型預(yù)模擬采用遺傳算法對模型參數(shù)進(jìn)行率定，預(yù)模擬參數(shù)率定值見表1，預(yù)模擬的效果圖見圖3。
表1 原始新安江模型預(yù)模擬參數(shù)率定值[8]Table 1 Pre-simulation parameters of the original Xin’anjiang model[8]
圖3 原始新安江模型預(yù)模擬的效果圖Fig.3 Pre-simulation results of the original Xin’anjiang model
從原始新安江模型預(yù)模擬的結(jié)果可以發(fā)現(xiàn)以下現(xiàn)象：
(1) 根據(jù)賈云飛等[9]對金沙江流域張力水蓄水容量的研究，一般流域平均張力水蓄水容量WM取值范圍在90～150 mm之間，而預(yù)模擬的WM值為31.88 mm，與經(jīng)驗(yàn)值相比偏小。
(2) 模擬流量對降雨的響應(yīng)過于敏感，多大于實(shí)際流量，即使無效降雨，模擬值也會有所響應(yīng)。
(3) 模擬洪峰流量較為滯后，起峰過程長，退水過程短，峰形較寬緩，拖尾現(xiàn)象不明顯。
結(jié)合廟溝巖溶流域的水文地質(zhì)條件，對上述現(xiàn)象分析如下：
非巖溶地區(qū)降雨降至地面之后，經(jīng)過緩慢入滲至地下水，這一過程在原始新安江模型中體現(xiàn)為三層蒸散發(fā)模型和表層、下層、深層3個串聯(lián)張力水水庫的計(jì)算，因此該模型對降雨信號的過濾能力強(qiáng)，較好地模擬了非巖溶地區(qū)含水介質(zhì)對降雨的調(diào)蓄功能。而在巖溶地區(qū)，降雨經(jīng)過地表產(chǎn)流通過巖溶管道直接匯入地下，相對于非巖溶區(qū)，巖溶區(qū)降雨在地表和淺層土壤滯留的時(shí)間更短，匯入地下轉(zhuǎn)化為地下水的速度更快，效率更高，所以巖溶地區(qū)地下水對降雨的響應(yīng)更為迅速，退水也更為迅速。因此，將新安江模型直接應(yīng)用于巖溶地區(qū)徑流過程的模擬會造成以下兩個結(jié)果：①加長了洪峰的響應(yīng)時(shí)間，導(dǎo)致模擬洪峰流量較為滯后；②減小了流域平均張力水蓄水容量WM值和自由水蓄水容量SM值，減小了產(chǎn)流閾值，縮短了模擬值的響應(yīng)時(shí)間。這將導(dǎo)致模擬結(jié)果對降雨的響應(yīng)過于敏感、拖尾現(xiàn)象不明顯。
2.2 新安江模型的改進(jìn)
Hartmann等[10]認(rèn)為，流域的巖溶發(fā)育程度越高，滲透性能越好。而巖溶地區(qū)降雨主要通過洼地、落水洞等巖溶地貌直接灌入并補(bǔ)給地下水，因此等效于降雨中的一部分受包氣帶調(diào)蓄，參與蒸發(fā)過程，另一部分直接補(bǔ)給地下水(見圖4)。故本文根據(jù)新安江模型預(yù)模擬的結(jié)論，將直接補(bǔ)給地下水的一部分降雨量定義為P1，總降雨量為P2，則參與蒸發(fā)過程的降雨量為P2-P1，并定義了降雨分配系數(shù)φ：
(1)
郭緒磊[11]、羅明明[12]的研究表明，南方巖溶地區(qū)地下水降雨入滲系數(shù)多在0.07～0.7之間，故本文中φ的取值范圍為0～0.7。
圖4 降雨-自由水轉(zhuǎn)換結(jié)構(gòu)示意圖Fig.4 Schematic diagram of rainfall-free water transition
另外，在地下徑流(RG)出流之前，增加一個線性水庫用以模擬深層溶洞水(快速徑流)與裂隙水(慢速徑流)的水量交換過程(見圖5)，進(jìn)入快速流水庫和慢速流水庫的水量由巖溶管道和裂隙的發(fā)育情況決定[7]。
圖5 溶洞快速和裂隙慢速徑流結(jié)構(gòu)示意圖Fig.5 Schematic diagram of fast runoff in karst caves and slow runoff in fracture
溶洞快速徑流流量和裂隙慢速徑流流量的計(jì)算公式如下:
(2)
式中：Q出，1、Q出，2分別代表快速徑流和慢速徑流的出庫流量(m3/s)；Q入，1、Q入，2分別代表快速徑流和慢速徑流的入庫流量(m3/s)；α1、α2分別為溶洞水和裂隙水的地下水衰退系數(shù)。
2.3 模型參數(shù)的敏感性分析
原始新安江模型改進(jìn)后共涉及17個參數(shù)，其中φ、α1和α2為新增參數(shù)。本文首先采用One-factor-At-a-Time(OAT)法對模型參數(shù)進(jìn)行局部敏感性分析。參數(shù)的敏感系數(shù)Si表達(dá)式為
(3)
式中:x為參數(shù)變化量;y為因變量模擬流量；Si為參數(shù)的敏感系數(shù)即敏感度，其大小代表敏感性的強(qiáng)弱。
改進(jìn)的新安江模型各參數(shù)取值同表1，模型參數(shù)的敏感度Si計(jì)算結(jié)果及其排序見表2[13]，模型參數(shù)敏感性分類見表3。
表2 改進(jìn)的新安江模型參數(shù)敏感度及其排序Table 2 Sensitivity of parameters of the improved Xin’anjiang model and the ranking
注：敏感度較低的模型參數(shù)未列出。
表3 改進(jìn)的新安江模型參數(shù)敏感性分類Table 3 Classification of the sensitivity of parameters of the improved Xin’anjiang model
結(jié)合表2和表3可知，改進(jìn)的新安江模型參數(shù)SM、WM、KG和φ為極敏感參數(shù)，KC、KI和CI為敏感參數(shù)。其中，張力水蓄水容量WM和自由水蓄水容量SM均屬于模型中線型水庫的出流閾值，WM控制產(chǎn)流，SM控制自由水水量，影響水源劃分，兩者直接影響著模型中水量的轉(zhuǎn)換，參數(shù)的敏感性較高；地下水出流系數(shù)KG控制著自由水S中地下水的出流比例，對于Nash-Sutcliffe效率系數(shù)DC有著十分顯著的影響[14]，因而參數(shù)的敏感性較高；降雨分配系數(shù)φ描述了通過巖溶管道直接灌入地下轉(zhuǎn)化為自由水的部分降雨，實(shí)質(zhì)上該參數(shù)減少了水分蒸發(fā)量，并將減少的水分直接轉(zhuǎn)化為自由水，減小了模型對降雨的調(diào)蓄作用，使得模擬徑流對降雨的響應(yīng)更加迅速，而這一點(diǎn)恰恰與巖溶地區(qū)的產(chǎn)流特點(diǎn)相吻合。
2.4 φ值對巖流流域徑流過程模擬的影響分析
為了進(jìn)一步探究參數(shù)φ值對巖溶流域徑流過程模擬效果的影響，本文針對廟溝巖溶流域140831號洪峰流量，利用改進(jìn)的新安江模型對其進(jìn)行了模擬，分析不同φ值對巖溶流域徑流過程模擬效果的影響，并與參數(shù)SM值和WM值進(jìn)行比較，其模擬結(jié)果見圖6。其中圖6(a)、(b)、(c)分別為只改變φ、SM、WM參數(shù)時(shí)，廟溝巖溶流域洪峰流量模擬曲線與實(shí)際流量曲線的對比圖。
由圖6(b)、(c)可見，隨著參數(shù)SM、WM的增大，模擬洪峰流量整體上逐漸減小，強(qiáng)降雨和弱降雨的模擬流量的“平均流量水平”逐漸降低。
表4為不同敏感參數(shù)取值時(shí)模擬得到的廟溝巖溶流域洪水信息統(tǒng)計(jì)表。
由表4可知，當(dāng)WM值逐漸增大時(shí)，平均流量、洪峰流量和確定性系數(shù)DC值逐漸減小，流量誤差逐漸增大；當(dāng)SM逐漸增大時(shí)，平均流量、洪峰流量和DC值逐漸增大，流量誤差逐漸減小，但當(dāng)SM取值大于80 mm時(shí)，DC值和流量誤差分別保持在96.4%和19.78%；在參數(shù)φ值逐漸增大的過程中，各項(xiàng)指標(biāo)均符合先增大后減小的規(guī)律，并且當(dāng)參數(shù)φ取值為0.5時(shí)，巖溶流域徑流過程的模擬效果達(dá)到最優(yōu)。
圖6 不同敏感參數(shù)取值時(shí)廟溝巖溶流域洪峰流量模擬曲線與實(shí)際流量曲線的對比圖Fig.6 Comparison between simulated and actual curves of runoff process in Miaogou karst basin with different sensitive parameter values
表4 不同敏感參數(shù)取值時(shí)模擬得到的廟溝巖溶流域洪水信息統(tǒng)計(jì)表Table 4 Statistical table of simulated flood information of Miaogou karst basin with different sensitivity parameter values
由此可見，WM值控制模擬流量的“平均流量水平”，影響模擬流量的縱向“高度”，決定著模型模擬效果的整體水平；SM為自由水蓄水容量，向下由地下水日出流系數(shù)KG控制地下水徑流，壤中流出流系數(shù)RI控制壤中流，當(dāng)自由水量大于SM值時(shí)，出流成為地表徑流RS，而巖溶地區(qū)地表徑流成分較小，因此當(dāng)自由水量大于SM值時(shí)，SM值的大小不會改變模型的模擬效果，DC值和流量誤差保持穩(wěn)定。
結(jié)合圖6(a)可以發(fā)現(xiàn)，當(dāng)參數(shù)φ值較小時(shí)，洪水峰值更多落在圖中Ⅰ區(qū)，較小洪峰流量的模擬效果較好，較大洪峰流量的模擬效果較差，改進(jìn)新安江模型對強(qiáng)降雨的響應(yīng)效率較低，模擬效果更接近于原始新安江模型；隨著參數(shù)φ值的增大，弱降雨模擬洪峰流量逐漸減小，強(qiáng)降雨模擬洪峰流量逐漸增大，洪水峰值逐漸由Ⅰ區(qū)向Ⅱ區(qū)過渡；當(dāng)參數(shù)φ值較大時(shí)，洪水峰值主要落在Ⅱ區(qū)，較小洪峰流量的模擬效果較差，較大洪峰流量的模擬效果較好，說明改進(jìn)新安江模型對強(qiáng)降雨的響應(yīng)效率更高。整體上模擬值的“平均流量水平”變化幅度相對較小。表4中數(shù)據(jù)也表明，參數(shù)φ值在整個取值域中流量誤差的變化范圍在3%以內(nèi)，而SM值和WM值在整個取值范圍內(nèi)誤差的最大變化范圍分別為7.5%和48.5%。
3 模型應(yīng)用與分析
3.1 模型參數(shù)率定
南方地區(qū)蒸散發(fā)折算系數(shù)KC一般取0.6，上層初始含水量WU、下層初始含水量WL和深層初始含水量WD是比較重要的指標(biāo)，其值的大小可以看作是整個流域內(nèi)蓄水量變化的起點(diǎn)，直接影響后期模型參數(shù)率定的準(zhǔn)確度以及模型模擬的效果[15]。
流域平均張力水蓄水容量WM可用實(shí)測資料分析，即選擇前期特別干旱的一次全流域產(chǎn)流過程進(jìn)行分析，根據(jù)廟溝巖溶流域?qū)嶋H降雨資料，確定WM取值為50 mm。利用2014年3月2日——2014年9月15日降雨徑流數(shù)據(jù)作為改進(jìn)的新安江模型參數(shù)率定期，采用遺傳算法對模型參數(shù)進(jìn)行率定，最后確定了Ⅰ區(qū)和Ⅱ區(qū)改進(jìn)新安江模型的參數(shù)，見表5。
表5 改進(jìn)的新安江模型參數(shù)率定值Table 5 Parameter values of the improved Xin’anjiang model
3.2 模擬結(jié)果與分析
根據(jù)參數(shù)率定表確定的改進(jìn)新安江模型參數(shù)，對Ⅰ區(qū)和Ⅱ區(qū)改進(jìn)的新安江模型分別賦值并計(jì)算2014年3月1日—2014年12月1日的降雨徑流，其模擬結(jié)果見圖7,模擬結(jié)果的精度見表6。
圖7 模擬徑流量與實(shí)際徑流量的關(guān)系圖Fig.7 Diagram of the relationship between simulated flow and actual flow
表6 新安江模型模擬精度評定表Table 6 Simulation accuracy evaluation table for Xin’anjiang model
由圖7可見，改進(jìn)的新安江模型的模擬效果較好，所模擬的洪峰起峰時(shí)間及峰值均與實(shí)際洪峰過程較為吻合。由表6可知，改進(jìn)后新安江模型模擬結(jié)果的確定性系數(shù)DC由0.592增長至0.901，平均絕對誤差和相對誤差均有較大幅度的減小?？傮w上看，改進(jìn)的新安江模型模擬值偏小，從而導(dǎo)致水量平衡偏小。
為了探究改進(jìn)的新安江模型對廟溝巖溶流域洪峰過程的模擬效果，利用改進(jìn)的新安江模型針對廟溝巖溶流域2014年4月10日、2014年7月14日、2014年8月13日以及2014年8月31日4次小洪峰過程進(jìn)行模擬，并分析模型參數(shù)對改進(jìn)新安江模型模擬效果的影響，其模擬結(jié)果見圖8和表7。
圖8 廟溝巖溶流域4次洪峰徑流過程模擬結(jié)果與實(shí)際流量的對比圖Fig.8 Comparison between the simulated and actual runoff during the runoff process of four flood peak periods in Miaogou karst basin
表7 廟溝巖溶流域4次洪峰徑流過程模擬值與實(shí)際值信息表Table 7 Information table of simulated runoff and actual runoff in the runoff process of four flood peak periods
注：a為原始新安江模型；b為改進(jìn)的新安江模型。
由圖8和表7可見，改進(jìn)新安江模型的確定性系數(shù)明顯優(yōu)于原始新安江模型，140410號洪峰可能是由于降雨-流量的相關(guān)關(guān)系較差，故而其模擬的效果較差。
4 結(jié) 論
(1) 原始新安江模型直接應(yīng)用于巖溶地區(qū)徑流過程的模擬效果較差，本文根據(jù)巖溶地區(qū)降水灌入并補(bǔ)給地下水的特點(diǎn)，對原始新安江模型入滲模塊結(jié)構(gòu)進(jìn)行了調(diào)整，增加了降雨-自由水直接轉(zhuǎn)換環(huán)節(jié)，改進(jìn)后新安江模型模擬結(jié)果的確定性系數(shù)由0.592增為0.901，單次徑流過程模擬精度最高可達(dá)0.993。
(2) 降雨分配系數(shù)φ實(shí)質(zhì)上將部分降雨直接轉(zhuǎn)化為自由水，減少了包氣帶對降雨的調(diào)蓄時(shí)間，縮短了模擬徑流對降雨的響應(yīng)時(shí)間，提高了新安江模型降雨-徑流的轉(zhuǎn)化效率。
(3) 降雨分配系數(shù)φ的敏感度為1.991，在新安江模型眾多參數(shù)中排在第4位，僅次于SM、WM和KG參數(shù)，且該參數(shù)更側(cè)重于改變確定性系數(shù)DC值，對模擬流量的影響較小。
綜上所述，改進(jìn)的新安江模型對廟溝巖溶流域徑流過程的模擬效果有顯著提升。
猜你喜歡
格陵蘭島積雪區(qū)地表徑流增加研究水利水電快報(bào)(2022年8期)2022-11-23 基于SWAT模型的布爾哈通河流域徑流模擬研究黑龍江大學(xué)自然科學(xué)學(xué)報(bào)(2022年1期)2022-03-29 雅魯藏布江河川徑流變化的季節(jié)性規(guī)律探索品牌研究(2020年32期)2020-08-09 新安江時(shí)代文學(xué)·上半月(2019年6期)2019-12-13 近40年來蒲河流域徑流變化及影響因素分析人民珠江(2019年4期)2019-04-20 無定河流域洪峰流量的空間變化統(tǒng)計(jì)分析城市地理(2017年10期)2017-11-04 鐵力水文站水文特性分析科技創(chuàng)新與應(yīng)用(2017年25期)2017-09-09 清流河滁縣站歷年洪峰水位洪峰流量趨勢分析及應(yīng)對措施南方農(nóng)業(yè)·下旬(2017年5期)2017-08-16 浙江省數(shù)學(xué)高考考法剖析科學(xué)與財(cái)富(2017年17期)2017-06-16 加強(qiáng)水生生物資源養(yǎng)護(hù)促進(jìn)新安江水環(huán)境生態(tài)平衡安徽農(nóng)學(xué)通報(bào)(2015年10期)2015-06-15
李玉坤，燕子琪，王紀(jì)元，陳 麗，周 宏
(中國地質(zhì)大學(xué)(武漢)地質(zhì)調(diào)查研究院，湖北 武漢 430074)
流域水文模型是在流域產(chǎn)匯流結(jié)構(gòu)概化的基礎(chǔ)上，以一個實(shí)體或數(shù)學(xué)模型來模擬流域降雨-徑流的形成過程，是一種常用的徑流過程模擬方法。流域水文模型通常分為集總式水文模型如新安江模型、SAC模型和水箱模型，分布式水文模型如SWAT模型、SHE模型，半分布式水文模型如TOPMODEL模型等。
巖溶地區(qū)含水介質(zhì)內(nèi)部結(jié)構(gòu)復(fù)雜多樣，實(shí)測資料獲取難度較大，使得流域水文模型在該地區(qū)的應(yīng)用成為難點(diǎn)。分布式水文模型基于水動力學(xué)機(jī)制來描述和模擬流域內(nèi)的水文過程，但對資料要求較高，因而在缺乏資料地區(qū)的應(yīng)用受到限制；黑箱模型又稱經(jīng)驗(yàn)?zāi)Ｐ?，?yīng)用門檻低、適用性強(qiáng)，但模型不能反映含水系統(tǒng)內(nèi)部的水文過程，也無法用于評價(jià)和管理人類活動對泉流量的影響[1]。新安江模型是一種概念模型，由于采用張力水和自由水蓄水容量曲線來描述張力水和自由水蓄水容量分布的空間異質(zhì)性，因此已成為一種具有分布意義的概念性流域水文模型。近年來，隨著GIS和RS技術(shù)以及計(jì)算機(jī)的發(fā)展，新安江模型參數(shù)的精度和計(jì)算分區(qū)的精度已得到顯著提升，如李致家等[2]將新安江模型應(yīng)用于柵格之上，與高精度雷達(dá)降雨資料耦合進(jìn)行了洪水預(yù)報(bào)；井立陽等[3]根據(jù)地形因子分布曲線推求了自由水蓄水容量;Yuan等[4]應(yīng)用雙源蒸散發(fā)模型改進(jìn)了新安江模型結(jié)構(gòu)，模擬了在植被參與的水循環(huán)過程中的徑流過程；Zhao等[5]將新安江模型應(yīng)用于包括土壤雨水侵蝕、土地資源利用等土壤資源評價(jià)中。目前新安江模型已由功能單一的洪水預(yù)報(bào)模型發(fā)展為適用于多領(lǐng)域、多功能的流域水文模型。
在巖溶地區(qū)，新安江模型的應(yīng)用也得到了發(fā)展。如莊一鸰等[6]在總結(jié)南方巖溶地區(qū)水文特性的基礎(chǔ)上，對新安江模型的穩(wěn)定下滲、水源劃分和匯流模塊進(jìn)行了改進(jìn)，并取得了較好的模擬效果；郝慶慶等[7]在新安江模型原有自由水庫的基礎(chǔ)上，再增加一個線性水庫用來模擬深層地下水的出流過程，使模擬精度得到了顯著提高。但目前新安江模型在巖溶地區(qū)的應(yīng)用研究仍較為有限。為此，本文選取廟溝巖溶流域?yàn)榈湫蛶r溶流域，以降雨入滲為切入點(diǎn)，對新安江模型進(jìn)行了改進(jìn)，以期提高新安江模型對巖溶地區(qū)徑流過程模擬的精度。
1 研究區(qū)概況
廟溝巖溶流域位于湖北省宜昌市興山縣黃糧鎮(zhèn)高陽鎮(zhèn)境內(nèi),為高嵐河流域的子流域，整體上處于中國大陸二級階地向三級階地過渡地帶。該流域面積約為25 km2，平均海拔為1 100 m，地貌類型為溶蝕侵蝕中山溝谷地貌。該流域西側(cè)為斜坡，巖性以志留系粉砂巖為主；東側(cè)為巖溶洼地地貌，巖性以寒武-奧陶系碳酸鹽巖為主，區(qū)內(nèi)落水洞、巖溶漏斗等巖溶形態(tài)十分發(fā)育,見圖1。
圖1 廟溝巖溶流域水系概況圖Fig.1 Drainage map of Miaogou karst basin
2 研究方法
2.1 新安江模型及預(yù)模擬
新安江模型以蓄滿產(chǎn)流為理論基礎(chǔ)，該模型結(jié)構(gòu)為分散性結(jié)構(gòu)。如圖2所示，降雨信號輸入模型后轉(zhuǎn)化為張力水，此時(shí)一部分水量耗于蒸散發(fā)，采用三層蒸發(fā)模型進(jìn)行蒸散發(fā)計(jì)算，另一部分水量轉(zhuǎn)化為自由水，并將水源劃分為地表徑流、壤中流和地下徑流三部分；最后采用馬斯京根法進(jìn)行河道匯流計(jì)算。
圖2 新安江模型結(jié)構(gòu)簡圖Fig.2 Structure diagram of Xin’anjiang model
為了驗(yàn)證原始新安江模型對巖溶流域徑流過程的模擬效果，采用原始新安江模型對廟溝巖溶流域進(jìn)行了一次預(yù)模擬，選取的時(shí)間序列為2014年3月2日到2014年9月2日，期間共有4次小洪峰，1次大洪峰。原始新安江模型預(yù)模擬采用遺傳算法對模型參數(shù)進(jìn)行率定，預(yù)模擬參數(shù)率定值見表1，預(yù)模擬的效果圖見圖3。
表1 原始新安江模型預(yù)模擬參數(shù)率定值[8]Table 1 Pre-simulation parameters of the original Xin’anjiang model[8]
圖3 原始新安江模型預(yù)模擬的效果圖Fig.3 Pre-simulation results of the original Xin’anjiang model
從原始新安江模型預(yù)模擬的結(jié)果可以發(fā)現(xiàn)以下現(xiàn)象：
(1) 根據(jù)賈云飛等[9]對金沙江流域張力水蓄水容量的研究，一般流域平均張力水蓄水容量WM取值范圍在90～150 mm之間，而預(yù)模擬的WM值為31.88 mm，與經(jīng)驗(yàn)值相比偏小。
(2) 模擬流量對降雨的響應(yīng)過于敏感，多大于實(shí)際流量，即使無效降雨，模擬值也會有所響應(yīng)。
(3) 模擬洪峰流量較為滯后，起峰過程長，退水過程短，峰形較寬緩，拖尾現(xiàn)象不明顯。
結(jié)合廟溝巖溶流域的水文地質(zhì)條件，對上述現(xiàn)象分析如下：
非巖溶地區(qū)降雨降至地面之后，經(jīng)過緩慢入滲至地下水，這一過程在原始新安江模型中體現(xiàn)為三層蒸散發(fā)模型和表層、下層、深層3個串聯(lián)張力水水庫的計(jì)算，因此該模型對降雨信號的過濾能力強(qiáng)，較好地模擬了非巖溶地區(qū)含水介質(zhì)對降雨的調(diào)蓄功能。而在巖溶地區(qū)，降雨經(jīng)過地表產(chǎn)流通過巖溶管道直接匯入地下，相對于非巖溶區(qū)，巖溶區(qū)降雨在地表和淺層土壤滯留的時(shí)間更短，匯入地下轉(zhuǎn)化為地下水的速度更快，效率更高，所以巖溶地區(qū)地下水對降雨的響應(yīng)更為迅速，退水也更為迅速。因此，將新安江模型直接應(yīng)用于巖溶地區(qū)徑流過程的模擬會造成以下兩個結(jié)果：①加長了洪峰的響應(yīng)時(shí)間，導(dǎo)致模擬洪峰流量較為滯后；②減小了流域平均張力水蓄水容量WM值和自由水蓄水容量SM值，減小了產(chǎn)流閾值，縮短了模擬值的響應(yīng)時(shí)間。這將導(dǎo)致模擬結(jié)果對降雨的響應(yīng)過于敏感、拖尾現(xiàn)象不明顯。
2.2 新安江模型的改進(jìn)
Hartmann等[10]認(rèn)為，流域的巖溶發(fā)育程度越高，滲透性能越好。而巖溶地區(qū)降雨主要通過洼地、落水洞等巖溶地貌直接灌入并補(bǔ)給地下水，因此等效于降雨中的一部分受包氣帶調(diào)蓄，參與蒸發(fā)過程，另一部分直接補(bǔ)給地下水(見圖4)。故本文根據(jù)新安江模型預(yù)模擬的結(jié)論，將直接補(bǔ)給地下水的一部分降雨量定義為P1，總降雨量為P2，則參與蒸發(fā)過程的降雨量為P2-P1，并定義了降雨分配系數(shù)φ：
(1)
郭緒磊[11]、羅明明[12]的研究表明，南方巖溶地區(qū)地下水降雨入滲系數(shù)多在0.07～0.7之間，故本文中φ的取值范圍為0～0.7。
圖4 降雨-自由水轉(zhuǎn)換結(jié)構(gòu)示意圖Fig.4 Schematic diagram of rainfall-free water transition
另外，在地下徑流(RG)出流之前，增加一個線性水庫用以模擬深層溶洞水(快速徑流)與裂隙水(慢速徑流)的水量交換過程(見圖5)，進(jìn)入快速流水庫和慢速流水庫的水量由巖溶管道和裂隙的發(fā)育情況決定[7]。
圖5 溶洞快速和裂隙慢速徑流結(jié)構(gòu)示意圖Fig.5 Schematic diagram of fast runoff in karst caves and slow runoff in fracture
溶洞快速徑流流量和裂隙慢速徑流流量的計(jì)算公式如下:
(2)
式中：Q出，1、Q出，2分別代表快速徑流和慢速徑流的出庫流量(m3/s)；Q入，1、Q入，2分別代表快速徑流和慢速徑流的入庫流量(m3/s)；α1、α2分別為溶洞水和裂隙水的地下水衰退系數(shù)。
2.3 模型參數(shù)的敏感性分析
原始新安江模型改進(jìn)后共涉及17個參數(shù)，其中φ、α1和α2為新增參數(shù)。本文首先采用One-factor-At-a-Time(OAT)法對模型參數(shù)進(jìn)行局部敏感性分析。參數(shù)的敏感系數(shù)Si表達(dá)式為
(3)
式中:x為參數(shù)變化量;y為因變量模擬流量；Si為參數(shù)的敏感系數(shù)即敏感度，其大小代表敏感性的強(qiáng)弱。
改進(jìn)的新安江模型各參數(shù)取值同表1，模型參數(shù)的敏感度Si計(jì)算結(jié)果及其排序見表2[13]，模型參數(shù)敏感性分類見表3。
表2 改進(jìn)的新安江模型參數(shù)敏感度及其排序Table 2 Sensitivity of parameters of the improved Xin’anjiang model and the ranking
注：敏感度較低的模型參數(shù)未列出。
表3 改進(jìn)的新安江模型參數(shù)敏感性分類Table 3 Classification of the sensitivity of parameters of the improved Xin’anjiang model
結(jié)合表2和表3可知，改進(jìn)的新安江模型參數(shù)SM、WM、KG和φ為極敏感參數(shù)，KC、KI和CI為敏感參數(shù)。其中，張力水蓄水容量WM和自由水蓄水容量SM均屬于模型中線型水庫的出流閾值，WM控制產(chǎn)流，SM控制自由水水量，影響水源劃分，兩者直接影響著模型中水量的轉(zhuǎn)換，參數(shù)的敏感性較高；地下水出流系數(shù)KG控制著自由水S中地下水的出流比例，對于Nash-Sutcliffe效率系數(shù)DC有著十分顯著的影響[14]，因而參數(shù)的敏感性較高；降雨分配系數(shù)φ描述了通過巖溶管道直接灌入地下轉(zhuǎn)化為自由水的部分降雨，實(shí)質(zhì)上該參數(shù)減少了水分蒸發(fā)量，并將減少的水分直接轉(zhuǎn)化為自由水，減小了模型對降雨的調(diào)蓄作用，使得模擬徑流對降雨的響應(yīng)更加迅速，而這一點(diǎn)恰恰與巖溶地區(qū)的產(chǎn)流特點(diǎn)相吻合。
2.4 φ值對巖流流域徑流過程模擬的影響分析
為了進(jìn)一步探究參數(shù)φ值對巖溶流域徑流過程模擬效果的影響，本文針對廟溝巖溶流域140831號洪峰流量，利用改進(jìn)的新安江模型對其進(jìn)行了模擬，分析不同φ值對巖溶流域徑流過程模擬效果的影響，并與參數(shù)SM值和WM值進(jìn)行比較，其模擬結(jié)果見圖6。其中圖6(a)、(b)、(c)分別為只改變φ、SM、WM參數(shù)時(shí)，廟溝巖溶流域洪峰流量模擬曲線與實(shí)際流量曲線的對比圖。
由圖6(b)、(c)可見，隨著參數(shù)SM、WM的增大，模擬洪峰流量整體上逐漸減小，強(qiáng)降雨和弱降雨的模擬流量的“平均流量水平”逐漸降低。
表4為不同敏感參數(shù)取值時(shí)模擬得到的廟溝巖溶流域洪水信息統(tǒng)計(jì)表。
由表4可知，當(dāng)WM值逐漸增大時(shí)，平均流量、洪峰流量和確定性系數(shù)DC值逐漸減小，流量誤差逐漸增大；當(dāng)SM逐漸增大時(shí)，平均流量、洪峰流量和DC值逐漸增大，流量誤差逐漸減小，但當(dāng)SM取值大于80 mm時(shí)，DC值和流量誤差分別保持在96.4%和19.78%；在參數(shù)φ值逐漸增大的過程中，各項(xiàng)指標(biāo)均符合先增大后減小的規(guī)律，并且當(dāng)參數(shù)φ取值為0.5時(shí)，巖溶流域徑流過程的模擬效果達(dá)到最優(yōu)。
圖6 不同敏感參數(shù)取值時(shí)廟溝巖溶流域洪峰流量模擬曲線與實(shí)際流量曲線的對比圖Fig.6 Comparison between simulated and actual curves of runoff process in Miaogou karst basin with different sensitive parameter values
表4 不同敏感參數(shù)取值時(shí)模擬得到的廟溝巖溶流域洪水信息統(tǒng)計(jì)表Table 4 Statistical table of simulated flood information of Miaogou karst basin with different sensitivity parameter values
由此可見，WM值控制模擬流量的“平均流量水平”，影響模擬流量的縱向“高度”，決定著模型模擬效果的整體水平；SM為自由水蓄水容量，向下由地下水日出流系數(shù)KG控制地下水徑流，壤中流出流系數(shù)RI控制壤中流，當(dāng)自由水量大于SM值時(shí)，出流成為地表徑流RS，而巖溶地區(qū)地表徑流成分較小，因此當(dāng)自由水量大于SM值時(shí)，SM值的大小不會改變模型的模擬效果，DC值和流量誤差保持穩(wěn)定。
結(jié)合圖6(a)可以發(fā)現(xiàn)，當(dāng)參數(shù)φ值較小時(shí)，洪水峰值更多落在圖中Ⅰ區(qū)，較小洪峰流量的模擬效果較好，較大洪峰流量的模擬效果較差，改進(jìn)新安江模型對強(qiáng)降雨的響應(yīng)效率較低，模擬效果更接近于原始新安江模型；隨著參數(shù)φ值的增大，弱降雨模擬洪峰流量逐漸減小，強(qiáng)降雨模擬洪峰流量逐漸增大，洪水峰值逐漸由Ⅰ區(qū)向Ⅱ區(qū)過渡；當(dāng)參數(shù)φ值較大時(shí)，洪水峰值主要落在Ⅱ區(qū)，較小洪峰流量的模擬效果較差，較大洪峰流量的模擬效果較好，說明改進(jìn)新安江模型對強(qiáng)降雨的響應(yīng)效率更高。整體上模擬值的“平均流量水平”變化幅度相對較小。表4中數(shù)據(jù)也表明，參數(shù)φ值在整個取值域中流量誤差的變化范圍在3%以內(nèi)，而SM值和WM值在整個取值范圍內(nèi)誤差的最大變化范圍分別為7.5%和48.5%。
3 模型應(yīng)用與分析
3.1 模型參數(shù)率定
南方地區(qū)蒸散發(fā)折算系數(shù)KC一般取0.6，上層初始含水量WU、下層初始含水量WL和深層初始含水量WD是比較重要的指標(biāo)，其值的大小可以看作是整個流域內(nèi)蓄水量變化的起點(diǎn)，直接影響后期模型參數(shù)率定的準(zhǔn)確度以及模型模擬的效果[15]。
流域平均張力水蓄水容量WM可用實(shí)測資料分析，即選擇前期特別干旱的一次全流域產(chǎn)流過程進(jìn)行分析，根據(jù)廟溝巖溶流域?qū)嶋H降雨資料，確定WM取值為50 mm。利用2014年3月2日——2014年9月15日降雨徑流數(shù)據(jù)作為改進(jìn)的新安江模型參數(shù)率定期，采用遺傳算法對模型參數(shù)進(jìn)行率定，最后確定了Ⅰ區(qū)和Ⅱ區(qū)改進(jìn)新安江模型的參數(shù)，見表5。
表5 改進(jìn)的新安江模型參數(shù)率定值Table 5 Parameter values of the improved Xin’anjiang model
3.2 模擬結(jié)果與分析
根據(jù)參數(shù)率定表確定的改進(jìn)新安江模型參數(shù)，對Ⅰ區(qū)和Ⅱ區(qū)改進(jìn)的新安江模型分別賦值并計(jì)算2014年3月1日—2014年12月1日的降雨徑流，其模擬結(jié)果見圖7,模擬結(jié)果的精度見表6。
圖7 模擬徑流量與實(shí)際徑流量的關(guān)系圖Fig.7 Diagram of the relationship between simulated flow and actual flow
表6 新安江模型模擬精度評定表Table 6 Simulation accuracy evaluation table for Xin’anjiang model
由圖7可見，改進(jìn)的新安江模型的模擬效果較好，所模擬的洪峰起峰時(shí)間及峰值均與實(shí)際洪峰過程較為吻合。由表6可知，改進(jìn)后新安江模型模擬結(jié)果的確定性系數(shù)DC由0.592增長至0.901，平均絕對誤差和相對誤差均有較大幅度的減小?？傮w上看，改進(jìn)的新安江模型模擬值偏小，從而導(dǎo)致水量平衡偏小。
為了探究改進(jìn)的新安江模型對廟溝巖溶流域洪峰過程的模擬效果，利用改進(jìn)的新安江模型針對廟溝巖溶流域2014年4月10日、2014年7月14日、2014年8月13日以及2014年8月31日4次小洪峰過程進(jìn)行模擬，并分析模型參數(shù)對改進(jìn)新安江模型模擬效果的影響，其模擬結(jié)果見圖8和表7。
圖8 廟溝巖溶流域4次洪峰徑流過程模擬結(jié)果與實(shí)際流量的對比圖Fig.8 Comparison between the simulated and actual runoff during the runoff process of four flood peak periods in Miaogou karst basin
表7 廟溝巖溶流域4次洪峰徑流過程模擬值與實(shí)際值信息表Table 7 Information table of simulated runoff and actual runoff in the runoff process of four flood peak periods
注：a為原始新安江模型；b為改進(jìn)的新安江模型。
由圖8和表7可見，改進(jìn)新安江模型的確定性系數(shù)明顯優(yōu)于原始新安江模型，140410號洪峰可能是由于降雨-流量的相關(guān)關(guān)系較差，故而其模擬的效果較差。
4 結(jié) 論
(1) 原始新安江模型直接應(yīng)用于巖溶地區(qū)徑流過程的模擬效果較差，本文根據(jù)巖溶地區(qū)降水灌入并補(bǔ)給地下水的特點(diǎn)，對原始新安江模型入滲模塊結(jié)構(gòu)進(jìn)行了調(diào)整，增加了降雨-自由水直接轉(zhuǎn)換環(huán)節(jié)，改進(jìn)后新安江模型模擬結(jié)果的確定性系數(shù)由0.592增為0.901，單次徑流過程模擬精度最高可達(dá)0.993。
(2) 降雨分配系數(shù)φ實(shí)質(zhì)上將部分降雨直接轉(zhuǎn)化為自由水，減少了包氣帶對降雨的調(diào)蓄時(shí)間，縮短了模擬徑流對降雨的響應(yīng)時(shí)間，提高了新安江模型降雨-徑流的轉(zhuǎn)化效率。
(3) 降雨分配系數(shù)φ的敏感度為1.991，在新安江模型眾多參數(shù)中排在第4位，僅次于SM、WM和KG參數(shù)，且該參數(shù)更側(cè)重于改變確定性系數(shù)DC值，對模擬流量的影響較小。
綜上所述，改進(jìn)的新安江模型對廟溝巖溶流域徑流過程的模擬效果有顯著提升。
猜你喜歡
水利水電快報(bào)(2022年8期)2022-11-23
黑龍江大學(xué)自然科學(xué)學(xué)報(bào)(2022年1期)2022-03-29
品牌研究(2020年32期)2020-08-09
時(shí)代文學(xué)·上半月(2019年6期)2019-12-13
人民珠江(2019年4期)2019-04-20
城市地理(2017年10期)2017-11-04
科技創(chuàng)新與應(yīng)用(2017年25期)2017-09-09
南方農(nóng)業(yè)·下旬(2017年5期)2017-08-16
科學(xué)與財(cái)富(2017年17期)2017-06-16
安徽農(nóng)學(xué)通報(bào)(2015年10期)2015-06-15
