王忠國(guó)，宋果毅，張文韜

目前，山東省正在進(jìn)行第三次水資源調(diào)查評(píng)價(jià)工作。由于人類活動(dòng)對(duì)河川徑流的影響越來(lái)越大，所以正確估算天然河川徑流量，使年徑流統(tǒng)計(jì)系列具有統(tǒng)一基礎(chǔ)，使年內(nèi)分配符合天然過(guò)程，就必須進(jìn)行還原計(jì)算。還原計(jì)算過(guò)程中，年月天然河川徑流量出現(xiàn)了負(fù)值，對(duì)此分析如下：

1 單站天然河川徑流量的計(jì)算

1.1 計(jì)算要求

按照第三次《山東省水資源評(píng)價(jià)技術(shù)細(xì)則（印發(fā)稿）》（以下簡(jiǎn)稱《細(xì)則》）以及山東省水文局地表水資源評(píng)價(jià)階段成果審查及技術(shù)研討會(huì)的要求，結(jié)合高崖水庫(kù)站的流域特性和資料情況，推求天然徑流量。對(duì)于控制面積內(nèi)不存在蓄水、引水、提水及河道分洪或堤防決口的水文站，實(shí)測(cè)河川徑流量即為天然河川徑流量；對(duì)于控制面積內(nèi)存在蓄水、引水、提水及分洪或決口的水文站，進(jìn)行逐月、逐年的實(shí)測(cè)河川徑流量進(jìn)行還原計(jì)算。高崖水庫(kù)站按逐年、逐月進(jìn)行還原計(jì)算。當(dāng)還原后的天然月河川徑流量出現(xiàn)負(fù)值時(shí)，對(duì)各項(xiàng)逐月還原水量進(jìn)行具體分析。

1.2 計(jì)算方法

1）高崖水庫(kù)站選為水資源評(píng)價(jià)的水文站，對(duì)單站年實(shí)測(cè)河川徑流系列進(jìn)行了一致性還原計(jì)算。做出了歷年逐月的天然河川徑流量。還原計(jì)算采用全面收集資料和典型調(diào)查分析相結(jié)合的方法，按照評(píng)價(jià)要求逐年逐月進(jìn)行。

2）按《細(xì)則》要求，對(duì)地表水利用的耗損量進(jìn)行還原。還原的主要項(xiàng)目包括：農(nóng)業(yè)灌溉、工業(yè)和生活用水的耗損量（含蒸發(fā)消耗和入滲損失），跨流域引入、引出水量，河道分洪水量，水庫(kù)蓄水變量等。第三次水資源評(píng)價(jià)時(shí)段內(nèi)，高崖水庫(kù)站還原的主要項(xiàng)目是農(nóng)業(yè)還原水量、工業(yè)和生活還原水量、水庫(kù)蓄水變量；本次評(píng)價(jià)時(shí)段內(nèi)無(wú)跨流域引入、引出水量，河道分洪水量。

3）還原計(jì)算時(shí)段內(nèi)天然徑流量的計(jì)算公式：

式中：W為天然河川徑流量，m3；W1為實(shí)測(cè)河川徑流量，m3；W2為農(nóng)業(yè)灌溉耗損量，m3；W3為工業(yè)用水耗損量，m3；W4為城鎮(zhèn)生活用水耗損量，m3；W5為跨流域（或跨區(qū)間）引水量，引出為正，引入為負(fù)，m3；W6為河道分洪不能回歸后的水量，分出為正，分入為負(fù)，m3；W7為大中型水庫(kù)蓄水變量，增加為正，減少為負(fù)，m3。

對(duì)于高崖水庫(kù)站式（1）簡(jiǎn)化為式（2）：

式中符號(hào)意義同上。

式（2）中僅列出了對(duì)水文站實(shí)測(cè)徑流量影響較大的還原項(xiàng)目，《細(xì)則》要求各地可根據(jù)具體情況增減項(xiàng)目。

天然河川徑流量W的計(jì)算公式（1）中，W2、W3、W4都是耗水量，均為正值；W1為實(shí)測(cè)河川徑流量有正有負(fù)；W7為大中型水庫(kù)蓄水變量，增加為正，減少為負(fù)。還原計(jì)算過(guò)程中，濰河區(qū)地表水資源評(píng)價(jià)年月天然河川徑流量出現(xiàn)了負(fù)值，計(jì)算結(jié)果見(jiàn)表1。

表1 濰河區(qū)高崖水庫(kù)站歷年逐月河川徑流量表

2 天然河川徑流量負(fù)值分析

2.1 天然河川徑流量概念

在第三次水資源評(píng)價(jià)中，地表水資源量是指河流、湖泊等地表水體中由當(dāng)?shù)亟邓纬傻?、可以逐年更新的?dòng)態(tài)水量，用天然河川徑流量表示。地表水資源量循環(huán)的時(shí)間單位是一年，由當(dāng)?shù)亟邓纬?，是正值?/p>

2.2 流域下墊面條件的影響

從水量交換的角度也可把全球陸海水量平衡的區(qū)域劃分為4個(gè)自然系統(tǒng)，即分為大氣系統(tǒng)，流域系統(tǒng)，土壤系統(tǒng)和地下水系統(tǒng)。

其中流域系統(tǒng)水量平衡方程式為：

式中：±△S為流域蓄水量，m3；P為降水量，m3；R 為徑流量，m3；E 為蒸發(fā)量，m3。

由于人類活動(dòng)改變了流域下墊面條件，導(dǎo)致入滲、徑流、蒸發(fā)等水平衡要素發(fā)生一定的變化，從而造成徑流的減少（或增加）。下墊面變化對(duì)產(chǎn)流的影響非常復(fù)雜，而在半干旱半濕潤(rùn)地區(qū)，許多流域的徑流因下墊面變化而衰減的現(xiàn)象已經(jīng)非常明顯。地表層下墊面的變化反映在流域蓄水量△S的變化上，是由于人類活動(dòng)影響造成的，如農(nóng)業(yè)灌溉、工業(yè)用水、城市化增加城鎮(zhèn)生活用水、跨流域（或跨區(qū)間）引出引進(jìn)水量、河道分洪、大中型水庫(kù)數(shù)量增加、庫(kù)湖串聯(lián)、河道梯級(jí)攔河閘壩等水利工程建設(shè)，造成流域蒸散發(fā)和滲入地下水量的變化。在干旱年份或枯水季節(jié)，公式（3）中，降水量P很?。ㄈ绺哐滤畮?kù)站2006年是枯水年，流域平均降水量412.7mm），農(nóng)業(yè)灌溉等人類活動(dòng)用水量就增大，加大河道、水庫(kù)、湖泊的取水量，流域內(nèi)蒸散發(fā)E也隨之增大，流域內(nèi)水庫(kù)、湖泊蓄水量減少，△S成為負(fù)值，水量平衡方程式中徑流量R變?yōu)樨?fù)值。即年月天然河川徑流量成為負(fù)值。對(duì)于非閉合流域或處于流域的上游的單站還原計(jì)算，下墊面變化對(duì)天然河川徑流量影響更加復(fù)雜，若產(chǎn)生的影響是增大蒸散發(fā)量或增大滲入地下水的補(bǔ)給量，年月天然河川徑流量成為負(fù)值就具有了合理性。

2.3 水庫(kù)蒸發(fā)量和滲漏量分析

水庫(kù)蒸發(fā)量和滲漏量作為下墊面因素，徑流量還原計(jì)算不考慮。第三次水資源評(píng)價(jià)，《細(xì)則》中年月天然河川徑流量計(jì)算公式（1），根據(jù)高崖水庫(kù)站流域特性和資料情況，推求年月天然河川徑流量公式（2）中：W1是高崖水庫(kù)站流入河道、流入渠道、城市供水和本站蓄水變量之和。其中，該站流入河道、流入渠道和城市供水都是正值；W2、W3、W4是高崖水庫(kù)站流域內(nèi)農(nóng)業(yè)灌溉耗損量、工業(yè)用水耗損量和城鎮(zhèn)生活用水耗損量，也都是正值；在干旱年份或枯水季節(jié)，降水量小，農(nóng)業(yè)灌溉等人類活動(dòng)用水量就增大，本站水庫(kù)蓄水量和流域內(nèi)大中型水庫(kù)蓄水量就將大量減小，大中型水庫(kù)蓄水變量是負(fù)值，且比較大。

同時(shí)，流域內(nèi)蒸散發(fā)量將增大。根據(jù)式（1），高崖水庫(kù)站本次水資源評(píng)價(jià)不考慮水庫(kù)蒸發(fā)損失量還原、流域內(nèi)蒸散發(fā)量還原和地表水滲入地下水的補(bǔ)給量，在干旱年份或枯水季節(jié)，其年月天然徑流量就出現(xiàn)了負(fù)值。

3 結(jié)語(yǔ)

當(dāng)前，進(jìn)行第三次山東省水資源調(diào)查評(píng)價(jià)工作，要在單站還原計(jì)算的基礎(chǔ)上，對(duì)其同步期逐年天然河川徑流量進(jìn)行系列一致性分析。同量級(jí)降水條件下2000年以后降雨徑流點(diǎn)據(jù)明顯偏離2000年以前的降雨徑流點(diǎn)據(jù)時(shí)，則表明下墊面變化對(duì)徑流影響較大，對(duì)原1956—2000年天然年河川徑流系列要進(jìn)行一致性修正。