馮宇鵬，鄭永路，叢娜，高誠

（中水北方勘測設計研究有限責任公司，天津 300222）

1 基本情況

清漳河屬漳河兩大支流中的北支，流域總面積5 142 km2，干流上有劉家莊、匡門口2個水文站。地表水供水工程主要有戀思水庫、九京供水工程、石匣水庫、漳北渠、漳西渠、漳南渠和大躍峰渠等。地下水水井約1 670眼，供水能力為1.06億m3，主要包括：山西省和順縣自來水廠若干水井，左權縣自來水廠若干水井、鑫源電廠水井以及若干煤礦水井；河北省涉縣自來水廠若干水井以及天鐵集團、龍山電廠、崇利冶鋼、焦化廠和一些中小型企業(yè)的自備井。再生水源主要來自3座污水處理廠，其中和順縣污水處理廠處理能力為1.0萬m3/d，潛在供水對象為潞安星光煤電公司，目前暫未供水；左權縣污水處理廠處理能力為1.0萬m3/d，供水對象為左權縣華能電廠；涉縣污水處理廠處理能力為5.0萬m3/d，供水對象為涉縣龍山電廠。流域內供、用水量2.14億m3，供水工程供水量情況詳見表1。

表1 清漳河流域供水工程供水量 億m3

海河水利委員會制定的《清漳河水資源配置方案》［1］中對清漳河上劉家莊、匡門口水文站斷面的下泄水量做出了要求，詳見表2。

表2 清漳河水資源配置方案中斷面下泄水量要求億m3

2 水文年景的判別方法——直接法

為了管理需要，需對當年的水文年景進行預判。直接法首先分析劉家莊水文站、匡門口水文站的還原計算方案，提出主要還原項，給出天然徑流量計算公式。之后，通過分析相關關系確定該方法的預判期，并給出截至預判期的天然徑流量與全年天然徑流量相關關系，即由預判期天然徑流量計算全年天然徑流量的公式。該方法可由預判期內獲得的諸數(shù)據(jù)預測該年天然徑流量，以此獲得水文年景，由水文年景來確定應下泄水量。

根據(jù)對流域的實際調查情況，建立了流域供需分析模型，模型中率定了各用水戶的回歸系數(shù)。結合調查與模型，給出劉家莊水文站、匡門口水文站徑流量還原計算公式：

式中：h為農業(yè)灌溉回歸系數(shù)，取0.1；如進行全年徑流 量 還 原 計 算，W劉家莊以上農業(yè)耗水量取0.23億m3，W劉家莊至匡門口區(qū)間農業(yè)耗水量取0.33億m3；如計算劉家莊水文站與匡門口水文站截至某一月份的天然徑流量，則蓄變量還原項中的末期蓄水量為截至該月底水庫蓄水量，供水量為截至該月份累計供水量；W劉家莊以上農業(yè)耗水量和W劉家莊至匡門口區(qū)間農業(yè)耗水量根據(jù)實際情況打折采用。

預判期的研究方法是逐一建立全年徑流量與截至某一月份的徑流量相關關系，經過分析，全年徑流量與1—8月的相關關系尚可（見圖1），與1—7月及以前的相關關系均不好，因此預判期定為8月底。選取1—8月天然徑流量來預測全年徑流量，其計算公式如下：

圖1 劉家莊和匡門口水文站1—8月徑流量與全年徑流量相關關系

根據(jù)灌溉制度，1—8月灌溉用水量占全年的約70%，因此W劉家莊以上農業(yè)耗水量取0.23億m3的70%，為0.16億m3；W劉家莊至匡門口區(qū)間農業(yè)耗水量取0.33億m3的70%，為0.23億m3。

3 水文年景的判別方法——間接法

由于人類活動和下墊面因素變化等方面的影響，地表徑流量在某些情況下不能很好地對下泄水量做出客觀的預判，而降水從降落、觀測到形成地表徑流的過程中受外界因素影響小，因此可作為水量下泄目標的指示性因子。間接法（降水徑流法）通過流域內截止某一月底降水資料等推算劉家莊水文站、匡門口水文站全年天然徑流量。

影響地表徑流量的因子主要為降水，汛期雨量充沛，是徑流形成的主要時期。而短期連續(xù)降水的集中強度和降水量級等往往決定地表徑流量的大小，日降水大于20 mm的累積降水量和連續(xù)10 d降水量能夠很好地表現(xiàn)這一特征；同時，上一年的徑流退水量和汛期降水量對次年的水量也有一定的影響關系。本次選取清漳河流域1—8月降水量、最大10 d降水量、主汛期7—8月降水量和上一年5—9月降水量作為影響因子進行初步分析。

由于影響因子的復雜性，在一般情況下只考慮一個因子是不夠的，必須考慮多個因子對預測對象的影響。這類問題就是多元回歸分析所要解決的問題。分別對全年徑流與上述諸因子中的兩者或兩者以上進行多元線性回歸分析。其預報模型為：

式中：b0，b1，b2，…，bm為回歸系數(shù)，根據(jù)Y與每個X過去的實測資料確定。本次計算置信度α=0.05。

以匡門口水文站1961—2000年天然徑流系列為預測對象，以諸因子不同組合建立多元線性回歸模型進行模擬計算。通過兩組共計24個方案的比選，得出以1—8月降水量、最大旬降水量、主汛期7—8月降水量和上一年5—9月降水量為組合的影響因子，計算結果相對誤差較小，模擬效果最好。

4 水文年景判別方法的取舍

從客觀情況來看，直接法所需諸參數(shù)中，每一參數(shù)的獲取環(huán)節(jié)均存在一定的誤差，而間接法主要使用降水資料進行計算，降水資料屬于較為客觀的數(shù)據(jù)，誤差較小。同時，從管理機構的可操作性來講[2]，間接法有一定優(yōu)勢，推薦采用間接法，直接法為備選方法，待實踐后進一步驗證2種方法的優(yōu)劣。間接法計算公式如下：

式中：W匡門口為匡門口水文站當年天然徑流量（億m3）；P1—8為當年1—8月降水量（mm）；P10為最大10 d降水量（mm）；P7—8為主汛期7—8月降水量（mm）；P上年5—9為上一年5—9月降水量（mm）。

由于匡門口水文站與劉家莊水文站年徑流相關系數(shù)達99%，因此劉家莊水文站的預測水量可根據(jù)兩者相關關系公式進行推求：

5 流域內可調度水量預測方法

清漳河流域內具有調蓄功能的供水工程主要為山西省和順縣的戀思水庫、左權縣的石匣水庫及澤城西安水電站（二期）工程[3]，其調節(jié)庫容分別為戀思水庫調節(jié)庫容1 234萬m3、石匣水庫調節(jié)庫容1 732萬m3、澤城西安水電站（二期）調節(jié)庫容2 099萬m3。根據(jù)前述的水文年型預判方法，預見期設置在8月底，因此需預測兩工程9—12月的徑流量（澤城西安水電站無供水任務，不必要預測其徑流量）。戀思水庫、石匣水庫上游基本無引水工程，開發(fā)利用程度很低，因此將戀思水庫、石匣水庫天然徑流量作為入庫水量，應用劉家莊水文站水文比擬法計算。

劉家莊水文站全年天然徑流量系列與9—12月天然徑流量系列的相關關系較好，因此該站全年徑流量計算公式為：

式中：W戀思9—12月徑流量為戀思水庫9—12月預測徑流量（億m3）；W石匣9—12月徑流量為石匣水庫9—12月預測徑流量（億m3）；W劉家莊9—12月徑流量為劉家莊水文站9—12月預測徑流量（億m3）；F戀思，F(xiàn)石匣分別為戀思水庫，石匣水庫流域面積（km2）；F劉家莊為劉家莊水文站以上流域面積（km2）；P戀思，P石匣為戀思水庫，石匣水庫多年平均面雨量（mm）；P劉家莊為劉家莊水文站多年平均面雨量（mm）。

劉家莊水文站、戀思水庫和石匣水庫相關參數(shù)，詳見表3。

表3 劉家莊水文站、戀思水庫和石匣水庫相關參數(shù)

可將戀思水庫、石匣水庫9—12月徑流量計算公式整理如下：

式中：W9—12月可調度水量為9—12月可調度水量（億m3）；W戀思蓄水量為 戀思水庫8月末的蓄水量（億m3）；W石匣蓄水量為 石 匣水庫8月 末的蓄水量（億m3）；W澤城西安水電站（二期）蓄水量為澤城西安水電站（二期）8月末蓄水量（億m3）；W戀思9—12月徑流量為戀思水庫9—12月預測徑流量（億m3）；W石匣9—12月徑流量為石匣水庫9—12月預測徑流量（億m3）；W戀思9—12月供水量為戀思水庫9—12月預計供水量（億m3）；W石匣9—12月供水量為石匣水庫9—12月預計供水量（億m3）。

6 水量調度方案

通過以上水文年景、可調度水量的預判方法的技術準備，可以制定可操作性強的清漳河水量調度方案，相關參數(shù)詳見表4—5。①預判年景，在預判期（8月底9月初）提供流域內雨量站的降水量資料。②確定應下泄水量。通過預測的年景對照水量分配的要求表格，得到當年全年應下泄水量。③判斷豐欠程度。將應下泄水量與截至8月底的已下泄水量（1—8月累計實測徑流量）進行對比，判斷豐欠程度。④可調度水量計算，按表4—5進行。⑤行政干預。通過流域機構建立的監(jiān)督管理辦法，完成水量下泄目標。從前述可調度水量計算公式可見，需要建立供水計劃申報制度，流域內的供水工程管理單位需在年初向流域管理機構申報全年的供水計劃。

表4 水文年景計算相關參數(shù)

表5 可調度水量計算相關參數(shù) 億m3