張俊芝，張子元，王慶明

(1.河北省水文勘測(cè)研究中心，河北 石家莊 050031；2.河北省水土保持工作總站，河北 石家莊 050031；3.中國(guó)水利水電科學(xué)研究院 流域水循環(huán)模擬與調(diào)控國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100038)

1 研究背景

準(zhǔn)確計(jì)算降雨入滲補(bǔ)給量是評(píng)價(jià)平原區(qū)地下水資源量的基礎(chǔ)[1]，尤其對(duì)于河北省中東部平原區(qū)來(lái)說(shuō)，確定降雨入滲補(bǔ)給量是合理制定區(qū)域地下水超采綜合治理方案的關(guān)鍵環(huán)節(jié)[2]。目前常用的年補(bǔ)給系數(shù)法[3]，即根據(jù)區(qū)域的年降雨量乘以年補(bǔ)給系數(shù)得出降雨入滲補(bǔ)給量，在實(shí)際應(yīng)用中難以解決兩個(gè)問(wèn)題：一是每年的降雨頻次、分布特征均不相同，采用固定的年降雨入滲補(bǔ)給系數(shù)是否導(dǎo)致地下水資源量評(píng)價(jià)結(jié)果的誤差[4]，二是前期土壤含水量、地下水埋深等因素時(shí)刻變化，同等的降雨未必形成同等的入滲補(bǔ)給量，如何確定單次降雨的入滲補(bǔ)給量[5]。

準(zhǔn)確評(píng)價(jià)降雨入滲補(bǔ)給量是地下水水文學(xué)的傳統(tǒng)研究方向，前人利用不同方法開(kāi)展過(guò)大量的研究工作。申豪勇等[6]利用氯離子質(zhì)量守恒的原理估算了裸露區(qū)中等巖溶發(fā)育區(qū)和覆蓋區(qū)的降雨入滲系數(shù)，取得的結(jié)果與該地區(qū)已有的成果基本吻合，為同類(lèi)型地區(qū)降雨入滲系數(shù)的準(zhǔn)確估算提供了新的方法。劉瑋等[7-8]綜合全國(guó)各流域片內(nèi)的分析成果，列出不同巖性在不同降雨量年份條件下的平均年降雨入滲補(bǔ)給系數(shù)的取值范圍。宋秋波等[9]采用水均衡方法，以年降雨作為一次降雨過(guò)程，進(jìn)一步改進(jìn)了計(jì)算降雨入滲補(bǔ)給系數(shù)的方法。王潤(rùn)冬等[10]利用MODCECLE(An Object Oriented Model for Basin Scale Water Cycle Simulation)模型在衡水野外試驗(yàn)田田間尺度下，分析研究了不同灌溉條件的地下水補(bǔ)給規(guī)律和補(bǔ)給量。李金柱[11]分析得出不同巖性、不同年降雨量在不同地下水埋深條件下的降雨入滲補(bǔ)給系數(shù)，并對(duì)降雨入滲補(bǔ)給的最佳地下水埋深和大埋深穩(wěn)定點(diǎn)存在的機(jī)理進(jìn)行了探討。崔浩浩等[12]基于地質(zhì)地貌特征，針對(duì)降雨入滲機(jī)理及地質(zhì)結(jié)構(gòu)對(duì)降雨入滲過(guò)程和速率影響進(jìn)行了研究。還有一些學(xué)者總結(jié)了降雨入滲補(bǔ)給的模式，綜合考慮降雨入滲及其滯后、降雨入滲優(yōu)先補(bǔ)給等，研究了土壤不同巖性、不同埋深、不同降雨時(shí)段與不同土壤含水量條件下的降雨入滲補(bǔ)給過(guò)程[13-14]?？偟膩?lái)說(shuō)，自然資源部門(mén)多集中于包氣帶巖性對(duì)年降雨入滲補(bǔ)給系數(shù)的影響[15]，水利部門(mén)多采用試驗(yàn)站法、動(dòng)態(tài)分析法、數(shù)值模型法等方法研究年降雨入滲補(bǔ)給系數(shù)的規(guī)律[16-17]，但目前專(zhuān)門(mén)針對(duì)大區(qū)域尺度的次降雨入滲補(bǔ)給過(guò)程分析研究相對(duì)較少，無(wú)法支撐決策部門(mén)掌握次降雨入滲補(bǔ)給量的實(shí)際工作。鑒于此，本研究在總結(jié)前人成果的基礎(chǔ)上，重點(diǎn)研究次降雨入滲補(bǔ)給規(guī)律，分析降雨間隔劃分、影響因素識(shí)別、地下水位影響等問(wèn)題，為更精確評(píng)估河北中東部平原降雨入滲補(bǔ)給量提供方法的探索和規(guī)律認(rèn)知的深化。

2 研究區(qū)域及方法

2.1 研究區(qū)域研究區(qū)域?yàn)楹颖敝袞|部平原兩個(gè)典型小流域，下博小流域和滄縣運(yùn)西小流域。下博流域位于衡水市，流域形狀近似長(zhǎng)方形，集水面積150 km2，流域相對(duì)封閉，地勢(shì)西高東低，流域地面平均高程23.30 m，平均坡度0.45‰。流域內(nèi)主要河流為朱家河，全長(zhǎng)28.5 km。2011—2015年下博流域的平均降水量為466 mm，其中汛期降水量為329 mm，占全年降水量的71%左右；平均蒸發(fā)量為815 mm，其中汛期蒸發(fā)量為397 mm，占全年總蒸發(fā)的49%左右，多年平均淺層地下水埋深為7.9～ 8.6 m。滄縣運(yùn)西小流域位于滄州市中部，集水面積582 km2，整個(gè)地形西南高東北低，平均坡度為0.1‰。2011—2015年平均降水量為550 mm，其中汛期降水量為422 mm，占全年降水量的77%左右，多年平均淺層地下水埋深為1.7～3.6 m。兩個(gè)小流域地下徑流微弱，地下水主要以降雨補(bǔ)給為主，消耗項(xiàng)主要是潛水蒸發(fā)和開(kāi)采，在一次降雨的短時(shí)間內(nèi)，水平排泄量和蒸發(fā)消耗量都很小，可以忽略不計(jì)[6]。研究區(qū)位置及站點(diǎn)分布圖見(jiàn)圖1。

圖1 典型小流域位置及監(jiān)測(cè)站點(diǎn)分布圖Fig.1 Location of typical small watershed and distribution of monitoring stations

2.2 計(jì)算方法

(1)動(dòng)態(tài)分析法計(jì)算次降雨入滲補(bǔ)給量

動(dòng)態(tài)分析法是在地下水水平排泄微弱的平原地區(qū)，考慮次降雨特征和前期土壤墑情，分次累計(jì)得到時(shí)段內(nèi)的降雨入滲補(bǔ)給量方法，相較于傳統(tǒng)年補(bǔ)給系數(shù)法，動(dòng)態(tài)分析法的次降雨入滲補(bǔ)給過(guò)程分析更科學(xué)，結(jié)果也更可靠[18-19]。

動(dòng)態(tài)分析法首先要對(duì)次降雨進(jìn)行劃分，根據(jù)以往研究成果分析[20-21]，降雨間隔選取小于3 d，降雨引起的地下水位動(dòng)態(tài)變化不易劃分，降雨間隔選取大于14 d，則可能將多次降雨混合，弱化了前期土壤含水量的影響，故本次選取3～14 d無(wú)雨作為次降雨劃分的時(shí)間間隔，采取試算的方法，確定最優(yōu)的次降雨劃分時(shí)間間隔。

在確定了次降雨量劃分間隔，根據(jù)地下水位變化幅度ΔH和水位變動(dòng)帶給水度μ(見(jiàn)表1，給水度來(lái)源于河北省第二次水資源評(píng)價(jià)成果，也是目前的通用值)，確定次降雨的入滲補(bǔ)給量Pr。當(dāng)計(jì)算時(shí)段內(nèi)有數(shù)次降雨，則將每次降雨引起的地下水位上升幅度相加，再乘以給水度，得到時(shí)段內(nèi)的降雨入滲補(bǔ)給量，公式如下。

(1)

式中：Pr為時(shí)段內(nèi)入滲補(bǔ)給量，m3；n為時(shí)段內(nèi)降雨總次數(shù)；i為次降雨；ΔH為次降雨后地下水位上升幅度，m；μ為給水度；A為研究區(qū)面積，m2。

表1 河北省平原區(qū)淺層地下水位變化帶給水度表

(2)水均衡法驗(yàn)證

水均衡一般指評(píng)價(jià)期內(nèi)多年平均地下水總補(bǔ)給量、總排泄量和淺層地下水蓄變量三者之間的均衡關(guān)系，如式2所示，不同方法計(jì)算出的水均衡差越小，說(shuō)明方法的精度越高，結(jié)果越可靠。水均衡法[22]適用于降雨、地下水位、補(bǔ)給量、排泄量等監(jiān)測(cè)資料以及水文地質(zhì)參數(shù)數(shù)據(jù)完整豐富的地區(qū)，研究選擇資料豐富、且有較多研究成果的衡水市下博和滄州滄縣運(yùn)西兩個(gè)小流域作為驗(yàn)證流域，分別采用動(dòng)態(tài)分析法和傳統(tǒng)年補(bǔ)給系數(shù)法計(jì)算水均衡分項(xiàng)，對(duì)比兩種方法計(jì)算的水均衡差。

水均衡差需分析多年平均地下水總補(bǔ)給量、總排泄量和淺層地下水蓄變量三者之間的均衡關(guān)系，如下式：

Q總補(bǔ)-Q總排±ΔW=X

(2)

δ=X/Q總補(bǔ)

(3)

式中：Q總補(bǔ)為淺層地下水總補(bǔ)給量，萬(wàn)m3；Q總排為淺層地下水總排泄量，萬(wàn)m3；ΔW為淺層地下水蓄變量，萬(wàn)m3；X為絕對(duì)均衡差，萬(wàn)m3；δ為相對(duì)均衡差。

在河北省中東部平原區(qū)，公式表述如下：

Q總補(bǔ)=Qh(q)+Qg+Pr

(4)

Q總排=Qk+E+Qy

(5)

其中：Qh(q)為河渠滲漏補(bǔ)給量，萬(wàn)m3；Qg為農(nóng)田灌溉入滲補(bǔ)給量，萬(wàn)m3；Pr為降雨入滲補(bǔ)給量，萬(wàn)m3，可采用動(dòng)態(tài)分析法和年入滲補(bǔ)給系數(shù)法等計(jì)算；Qk為淺層地下水開(kāi)采量，萬(wàn)m3；Qy越流排泄量，萬(wàn)m3；E為潛水蒸發(fā)量，萬(wàn)m3。其中，除降雨入滲補(bǔ)給量外，其他各分項(xiàng)在兩個(gè)小流域均有成熟的監(jiān)測(cè)或計(jì)算方法，本研究不在贅述，在對(duì)比水均衡差時(shí)，除降雨入滲補(bǔ)給量，其他各分項(xiàng)采用相同的值。

3 結(jié)果與分析

3.1 基于動(dòng)態(tài)分析法的次降雨入滲補(bǔ)給量計(jì)算

3.1.1 次降雨間隔的劃分 為精確的劃分次降雨時(shí)段長(zhǎng)度，選取雨量站、地下水位觀測(cè)井距離較近，數(shù)據(jù)關(guān)系較好的測(cè)站，作為劃分場(chǎng)次降雨的數(shù)據(jù)來(lái)源，本次研究選取兩個(gè)小流域內(nèi)5個(gè)典型站點(diǎn)的次降雨量和次降雨入滲補(bǔ)給量進(jìn)行分析，見(jiàn)表2所示。選取3～14 d無(wú)雨日作為次降雨劃分時(shí)段間隔，通過(guò)分析次降雨入滲補(bǔ)給量與次降雨量擬合度確定最優(yōu)劃分時(shí)段，結(jié)果如圖2所示。選取2003—2015年的127次降雨，在不同劃分時(shí)間間隔下，次降雨入滲補(bǔ)給量(Pr)與次降雨量(P)關(guān)系較好，其中G-001、G-003、G-004三個(gè)站點(diǎn)7 d無(wú)雨時(shí)Pr～P相關(guān)系數(shù)最高分別為0.79、0.83和0.86，G-002和G-005兩個(gè)站點(diǎn)8 d無(wú)雨日時(shí)Pr～P相關(guān)系數(shù)最高，分別為0.84和0.70，綜合來(lái)說(shuō)，可選取7 d無(wú)雨日作為次降雨劃分適合河北中部平原地區(qū)的降雨特征。

表2 地下水監(jiān)測(cè)井與土壤墑情監(jiān)測(cè)站對(duì)應(yīng)情況

圖2 不同無(wú)雨日間隔下Pr～P相關(guān)性分析Fig.2 Correlation analysis of Pr～P at different rainless day intervals

3.1.2 動(dòng)態(tài)分析法計(jì)算結(jié)果及對(duì)比驗(yàn)證 選取2011年、2012年和2014年分別代表平水年、豐水年和枯水年，分別采用動(dòng)態(tài)分析法和傳統(tǒng)年補(bǔ)給系數(shù)法計(jì)算降雨入滲補(bǔ)給量，與排泄量和地下水蓄變量對(duì)比求水均衡差，評(píng)價(jià)兩種計(jì)算方法結(jié)果的精度。動(dòng)態(tài)分析法由公式(1)計(jì)算得出，年補(bǔ)給系數(shù)法即根據(jù)年降雨量乘以年補(bǔ)給系數(shù)得出，年補(bǔ)給系數(shù)來(lái)源于河北省水資源評(píng)價(jià)結(jié)果。下博流域水均衡法各分項(xiàng)數(shù)據(jù)如表3所示，2011年、2012年和2014年下博流域采用動(dòng)態(tài)分析法計(jì)算的地下水總補(bǔ)給量分別為555萬(wàn)m3、810萬(wàn)m3和150萬(wàn)m3，采用傳統(tǒng)年補(bǔ)給系數(shù)法計(jì)算的地下水總補(bǔ)給量分別為1245萬(wàn)m3、1830萬(wàn)m3和465萬(wàn)m3，兩種方法計(jì)算結(jié)果相差較大，在2012年甚至相差一倍以上，三個(gè)年份平均，年補(bǔ)給系數(shù)法較動(dòng)態(tài)分析法計(jì)算的降雨入滲補(bǔ)給量多675萬(wàn)m3。三個(gè)水平年平均總排泄量為866萬(wàn)m3，蓄變量減少120萬(wàn)m3，其中降雨是地下水的最大補(bǔ)給項(xiàng)，人工開(kāi)采是最大的地下水排泄項(xiàng)。

表3 下博小流域水均衡計(jì)算結(jié)果表

滄縣運(yùn)西地區(qū)水均衡法計(jì)算各分項(xiàng)如表4所示，采用動(dòng)態(tài)分析法得到的降雨入滲補(bǔ)給量分別為3110萬(wàn)m3、3648萬(wàn)m3和897萬(wàn)m3，采用年補(bǔ)給系數(shù)法得到的降雨入滲補(bǔ)給量分別為4246萬(wàn)m3、6877萬(wàn)m3和1375萬(wàn)m3，兩種方法的計(jì)算結(jié)果同樣相差巨大，三個(gè)水平年平均降雨入滲補(bǔ)給量為2552萬(wàn)m3和4087萬(wàn)m3，相差超過(guò)60%?？偱判沽繛?213萬(wàn)m3，蓄變量增加767萬(wàn)m3，其中降雨入滲是最大的補(bǔ)給項(xiàng)，潛水蒸發(fā)是最大的排泄項(xiàng)。

表4 滄縣運(yùn)西小流域水均衡計(jì)算結(jié)果表

將動(dòng)態(tài)分析法和年補(bǔ)給系數(shù)法的結(jié)果進(jìn)行對(duì)比，結(jié)果如表5所示，在下博小流域動(dòng)態(tài)分析法計(jì)算的降雨入滲補(bǔ)給量帶入水均衡法得出的水均衡差為1%～15%，年補(bǔ)給系數(shù)法得到的水均衡差為43%～52%，動(dòng)態(tài)分析法平均水均衡差精度提高40%；在滄縣運(yùn)西小流域動(dòng)態(tài)分析法的水均衡差在-1%～-13%，而年補(bǔ)給系數(shù)法得到的水均衡差20%～39%，動(dòng)態(tài)分析法平均水均衡差精度提高22%。從水均衡差上看動(dòng)態(tài)分析法的結(jié)果更準(zhǔn)確，補(bǔ)給項(xiàng)與排泄項(xiàng)平衡結(jié)果更好。動(dòng)態(tài)分析法根據(jù)次降雨入滲補(bǔ)給量累積得出年入滲補(bǔ)給量，能夠考慮雨型和前期土壤水量的影響，而年補(bǔ)給系數(shù)法，是根據(jù)歷史經(jīng)驗(yàn)的入滲補(bǔ)給系數(shù)乘以年降雨量得出，誤差較大。另外，年補(bǔ)給系數(shù)法計(jì)算的降雨入滲補(bǔ)給量遠(yuǎn)高于動(dòng)態(tài)分析法計(jì)算結(jié)果，意味著河北中東部平原的真實(shí)入滲補(bǔ)給量可能長(zhǎng)期被高估。

表5 兩種方法計(jì)算結(jié)果均衡差對(duì)比結(jié)果

3.2 小流域次降雨入滲補(bǔ)給規(guī)律

3.2.1 前期土壤墑情與降雨入滲補(bǔ)給關(guān)系 降雨入滲補(bǔ)給量與前期土壤墑情密切相關(guān)，通常前期雨量(Pa)越大、前期土壤含水量(SW)越高，土壤吸收的水量越少，形成的降雨入滲補(bǔ)給量越多，入滲補(bǔ)給系數(shù)越大。本研究將入滲補(bǔ)給量(Pr)分別和次降雨量(P)、次降雨量加前期影響雨量(P+Pa)以及次降雨量加前期土壤含水量(P+SW)建立關(guān)系曲線，確定最佳擬合關(guān)系。如圖3所示，下博小流域和滄縣運(yùn)西小流域點(diǎn)的Pr～P、Pr～(P+Pa)、Pr～(P+SW)均符合線性關(guān)系，考慮前期影響雨量或前期土壤含水量后，與降雨入滲補(bǔ)給量的相關(guān)性均有所提高，Pr～P決定系數(shù)為0.64，Pr～(P+Pa)和Pr～(P+SW)的決定系數(shù)分別為0.70和0.76，降雨量加前期土壤含水量與次降雨入滲補(bǔ)給量的相關(guān)關(guān)系更強(qiáng)。故在評(píng)價(jià)地下水補(bǔ)給量時(shí)，應(yīng)結(jié)合土壤墑情站，考慮前期土壤墑情才能得到相對(duì)準(zhǔn)確的結(jié)果。

圖3 降雨入滲與土壤墑情的關(guān)系Fig.3 Correlation between rainfall infiltration and soil moisture

3.2.2 補(bǔ)給地下水的降雨閾值 下博小流域地下水埋深約為7.9～8.6 m左右，如圖4所示，2011年、2012年和2014年下博小流域共有44場(chǎng)降雨，其中有19場(chǎng)降雨監(jiān)測(cè)到地下水位波動(dòng)，即次降雨補(bǔ)給到地下水，其它次降雨均未引發(fā)地下水位波動(dòng)，最高次降雨量為364 mm(P+SW為451 mm)，入滲補(bǔ)給地下水量為42 mm，平均(P+SW)超過(guò)165 mm才會(huì)入滲補(bǔ)給到地下水，三年平均次降雨入滲補(bǔ)給系數(shù)為0.025(Pr/(P+SW))。滄縣運(yùn)西小流域地下水埋深約為1.7～3.6 m，如圖5所示，在滄縣運(yùn)西小流域三個(gè)水平年共產(chǎn)生45次降雨，有33次降雨補(bǔ)給到地下水，最高次降雨量為436 mm(P+SW為609 mm)，入滲補(bǔ)給地下水量為30 mm，平均(P+SW)超過(guò)94 mm才會(huì)入滲補(bǔ)給到地下水，三年平均次降雨入滲補(bǔ)給系數(shù)為0.038(Pr/(P+SW))。

圖4 下博小流域動(dòng)態(tài)分析法入滲補(bǔ)給量計(jì)算結(jié)果Fig.4 Calculation results of infiltration recharge by dynamic analysis method in Xiabo watershed

圖5 滄縣運(yùn)西小流域動(dòng)態(tài)分析法入滲補(bǔ)給量計(jì)算結(jié)果Fig.5 Calculation results of infiltration recharge by dynamic analytical method in Cangxian Yunxi small watershed

3.2.3 不同降雨量對(duì)入滲補(bǔ)給量的影響 兩個(gè)小流域不同降雨量條件下Pr～(P+SW)關(guān)系，如圖6所示，滄縣運(yùn)西小流域淺層地下水埋深較淺，同等降雨量條件，入滲補(bǔ)給量多于下博小流域，并且隨著降雨量增多差距越來(lái)越大。當(dāng)埋深確定的情況下，(P+SW)～Pr成S型曲線關(guān)系，當(dāng)(P+SW)<250 mm時(shí)，Pr隨著(P+SW)增大緩慢增加，此時(shí)下博小流域的降雨入滲補(bǔ)給系數(shù)為0.017，滄縣運(yùn)西小流域?yàn)?.029；當(dāng)250 mm≤(P+SW)≤500 mm時(shí)，Pr隨著(P+SW)增大增速變快，此時(shí)下博小流域的降雨入滲補(bǔ)給系數(shù)為0.043，滄縣運(yùn)西小流域?yàn)?.075；當(dāng)(P+SW)>500 mm時(shí)，Pr隨著(P+SW)增大增速變緩，此時(shí)下博小流域的降雨入滲補(bǔ)給系數(shù)為0.066，滄縣運(yùn)西小流域?yàn)?.11。(P+SW)較小時(shí)，土壤初始含水量低，次降雨首先被土壤吸收，不會(huì)產(chǎn)生入滲補(bǔ)給量，隨著(P+SW)增大，土壤逐漸達(dá)到田間持水率，土壤吸力減小，入滲補(bǔ)給量增大，直至達(dá)到穩(wěn)定下滲量，此時(shí)，(P+SW)再增大，入滲補(bǔ)給速率是穩(wěn)定的，多余的降雨將形成地表產(chǎn)流。

圖6 不同地下水埋深條件下Pr～(P+SW)擬合曲線Fig.6 Fitting curves of Pr～(P+SW) under different groundwater depths

由于地下水補(bǔ)給過(guò)程十分復(fù)雜，受評(píng)價(jià)時(shí)段、降雨量、地下水埋深、數(shù)據(jù)觀測(cè)精度等多種因素影響，即使在相近的區(qū)域不同學(xué)者評(píng)估出的地下水補(bǔ)給量也有相當(dāng)大的差異[23-25]。目前，采用水均衡法校正評(píng)價(jià)結(jié)果的精度是常用的方法，本研究的結(jié)果說(shuō)明在研究時(shí)段內(nèi)，動(dòng)態(tài)分析法的評(píng)價(jià)精度是令人滿意的，但是在不同時(shí)段、不同地下水埋深等條件下，動(dòng)態(tài)分析法也應(yīng)有相應(yīng)的調(diào)整。動(dòng)態(tài)分析法所需要地下水位觀測(cè)數(shù)據(jù)，目前在部分地區(qū)推廣還存在困難，但隨著河北省地下水監(jiān)測(cè)站網(wǎng)的逐步完善，未來(lái)評(píng)價(jià)地下水補(bǔ)給量的精度將逐步提高，技術(shù)方法也會(huì)逐步更新發(fā)展，動(dòng)態(tài)分析法的應(yīng)用前景也會(huì)更加廣泛。

4 結(jié)論

本研究采用動(dòng)態(tài)分析法評(píng)價(jià)了河北平原兩個(gè)典型小流域2011、2012和2014年地下水補(bǔ)給量，通過(guò)水均衡差驗(yàn)證了動(dòng)態(tài)分析法的適用性，并分析了降雨間隔、前期土壤墑情和地下水埋深等因素對(duì)地下水補(bǔ)給的影響，主要結(jié)論如下：

(1)河北中東部平原適合以7 d無(wú)雨日作為劃分場(chǎng)次降雨的時(shí)間間隔標(biāo)準(zhǔn)，將次降雨量(P)加上前期土壤含水量(SW)后，能夠提高單純以次降雨量推求地下水補(bǔ)給量的精度。

(2)地下水埋深是影響地下水補(bǔ)給的重要因素，下博小流域地下水埋深約為7.9～8.6 m，(P+SW)超過(guò)165 mm才會(huì)入滲補(bǔ)給到地下水；滄縣運(yùn)西小流域地下水埋深約為1.7～3.6 m，平均(P+SW)超過(guò)94 mm才入滲補(bǔ)給到地下水。

(3)(P+SW)與地下水補(bǔ)給量(Pr)呈S型曲線關(guān)系，當(dāng)(P+SW)<250 mm時(shí)，Pr隨著(P+SW)增大緩慢增加，當(dāng)250 mm≤(P+SW)≤500 mm時(shí)，Pr隨著(P+SW)增大增速變快，(P+SW)>500 mm時(shí)，Pr隨著(P+SW)增大增速變緩。