基于河網(wǎng)水流水質(zhì)耦合模型的塔河干流引水方案優(yōu)選

張小清

（新疆塔里木河流域干流管理局，新疆 庫爾勒 841000）

0 引言

塔里木河流域是我國第一大內(nèi)流性河流，其干流全長1231 km，自身不產(chǎn)流，屬于純耗散型流域，其水資源完全依靠外部補給。目前主要靠上游阿克蘇河、和田河及葉爾羌河進行流域干流水資源供給，自1976年通過連接庫爾勒和塔河的干渠進行孔雀河水源調(diào)入后，便逐漸形成“四源一干”的水量供給格局[1]。近年來，塔河干流來自源流的補給水量逐年減少，再加上干流上中游無序低效的水土資源開發(fā)，導(dǎo)致塔河干流水生態(tài)環(huán)境日益惡化，尾閭羅布泊和臺特瑪湖面臨干涸的危機。2001年，國家出臺的《塔里木河流域近期綜合治理規(guī)劃》基于戰(zhàn)略高度提出了塔河流域水資源開發(fā)利用、治理、節(jié)約、保護的總體思路、布局及實施意見，將塔河干流水資源可持續(xù)開發(fā)利用提升到戰(zhàn)略高度。2010年，塔河流域管理局基于所建設(shè)的《塔河流域水量調(diào)度管理系統(tǒng)》進行干流水量調(diào)度方案的實施。本文正是基于這一背景，對塔里木河干流三種水量調(diào)度方案進行比選，并進行了水量調(diào)度對水環(huán)境可能影響的量化分析，對于流域干流生態(tài)環(huán)境的恢復(fù)、水資源節(jié)約及綜合開發(fā)利用具有現(xiàn)實意義。

1 研究方法

1.1 河網(wǎng)水流水質(zhì)耦合模型

此處選用流域水量調(diào)度規(guī)劃管理中常用的由河海大學陳文輝教授所開發(fā)的水文水動力學與水質(zhì)耦合模型[2]，進行河網(wǎng)水流水質(zhì)耦合模擬。該模型以塔里木河干流流域為計算范圍，包括降水產(chǎn)匯流、河網(wǎng)動力模型和水質(zhì)模型等部分。

1.1.1 降水產(chǎn)匯流

（1）產(chǎn)流模型。分別計算不同下墊面的水面、水田、旱地、道路的產(chǎn)流情況:水面產(chǎn)流是日降水量與蒸發(fā)量之差；水田產(chǎn)流則結(jié)合作物生長所需水量變化、水田水深限度及耐淹水深等按天推求而得；旱地產(chǎn)流則考慮運用新安江蓄滿產(chǎn)流模型；道路產(chǎn)流根據(jù)下墊面情況，分別推求透水層、填挖不透水層、無填挖不透水層等情況下的產(chǎn)流過程。

（2）匯流模型。為進行干流匯流計算，假設(shè)日凈雨量按照40%、40%和20%過程進行分配，歷時3日分別匯入河網(wǎng)，還需同時考慮排澇模數(shù)等的制約。

1.1.2 河網(wǎng)水動力模型

以降水徑流模型計算結(jié)果為依據(jù)并結(jié)合流域水負荷的面與點數(shù)據(jù)，進行河網(wǎng)水質(zhì)成果的耦合模擬，再考慮流域內(nèi)引排水工程的實際作用，進行河網(wǎng)水流運動態(tài)勢的模擬，并推求各斷面流量與水位。河網(wǎng)一維非恒定流模型如下:

其中:q為流域旁側(cè)入流；Q為流域河道斷面流量，m3/s；A為河道過水面積，m2；B 為河寬，m；Z 為水位，m；VX為沿水流方向的流域旁側(cè)入流流速分配，m/s；K為流量模數(shù)；α為動量調(diào)整系數(shù)，反映河道斷面流速分布均勻程度。

1.1.3 水質(zhì)模型

水質(zhì)模型主要是調(diào)蓄節(jié)點、河網(wǎng)及塔河干流水質(zhì)模型，通過進行水量耦合模型與水質(zhì)模型的聯(lián)算，進行數(shù)值離散程度的分析。水質(zhì)模型為:

其中:A為斷面面積，m2；C為水質(zhì)指標濃度，mg/L；t為觀測時間，s；Ex為縱向分散系數(shù)，m2/s；U 為斷面流速均值，m/s；S 為具體水質(zhì)指標生物化學反映項，g/（m3·d）；Sw為具體水質(zhì)指標外部源項，g/s；αe為系數(shù)，取 αe=0.01；C0為謝才系數(shù)；θ為流域斷面寬深比；q為斷面單寬流量均值，m3/s。

1.2 研究時段及模型參數(shù)

1.2.1 研究時段

在1980年～2010年的數(shù)據(jù)序列中，2010年塔河干流年降水量、流域水位保證率接近50%，降水過程與多年平均降水量接近，能反映和代表流域平水年基本情況。

1.2.2 模型參數(shù)

為保證研究結(jié)論在流域水量調(diào)度方面的適用性，以干流相關(guān)規(guī)劃及具體進展情況為依據(jù)，采用《新疆塔里木河流域水環(huán)境綜合治理方案》所規(guī)定的流域近期工況。流域干流下墊面采用綜合治理方案中所預(yù)測的下墊面資料為邊界條件，污染源邊界則沿用年度實際污染源數(shù)據(jù)資料。

（1）水動力學參數(shù)

采用塔河干流阿拉爾斷面和新其滿斷面逐日實測流量及水位過程進行水動力學驗證，水位驗證結(jié)果見圖1和圖2。由圖可知，水位計算值和實測值擬合程度高，兩個斷面絕對誤差均在0.55 m以內(nèi)，本文所提出的模型對于干流水位擬合精度較高。

圖1 阿拉爾斷面水位計算值和實測值擬合程度

圖2 新其滿斷面水位計算值和實測值擬合程度

限于篇幅所限，英巴扎典型段、烏斯?jié)M樞紐、恰拉水庫、大西海子水庫、臺特瑪湖斷面水位計算值與實測值擬合程度對比圖略。所選取的塔河干流八個斷面流量計算值變化規(guī)律與實測值變化規(guī)律一致性較好，模型精度高。

（2）水質(zhì)參數(shù)

采用2017年逐月實測水質(zhì)數(shù)據(jù)對塔河干流各斷面水質(zhì)模型進行驗證，并將模型計算結(jié)果與實測值進行對比，從表1各斷面氨氮、COD、TN、TP等水質(zhì)指標實測值與計算值相對誤差均值結(jié)果可以看出:氨氮指標方面阿拉爾斷面、帕滿水庫、烏斯?jié)M樞紐相對誤差較大；COD指標方面，英巴扎典型段、大西海子水庫和帕滿水庫、新其滿斷面誤差均值較大；英巴扎典型段、阿拉爾斷面TN指標相對誤差均值較大；帕滿水庫、恰拉水庫、大西海子水庫和臺特瑪湖TP指標相對誤差均值較大。但是上述誤差均值取值均在30%以內(nèi)（除帕滿水庫TP指標相對誤差均值為30.45%外），表明本文所提出模型的模擬值與實測值擬合程度較好，能基本滿足塔河干流水質(zhì)模擬的要求。

表1 各斷面水質(zhì)指標相對誤差均值 單位:%

2 水量調(diào)度方案管理研究

2.1 水量調(diào)度方案

常規(guī)調(diào)度模式難以形成水量循環(huán)，即便增加水量也無法改善干流水環(huán)境，難以從根源上解決流域水質(zhì)型缺水問題。為了充分利用干流現(xiàn)有水利工程，優(yōu)化調(diào)度方案，和塔河干流河網(wǎng)有網(wǎng)有綱、縱橫交錯、四通八達以及較好的引排能力和區(qū)位優(yōu)勢，構(gòu)建既引又排的中循環(huán)調(diào)水線路。

依據(jù)《塔里木河防洪規(guī)劃報告》以及上述對塔河干流水動力學、水質(zhì)模擬的結(jié)果，針對塔里木河干流提出以下三種可行的引水方案，方案一是從阿拉爾斷面、新其滿斷面及帕滿水庫引水；方案二是從干流中游烏斯?jié)M樞紐和英巴扎典型段引水，方案三是從恰拉水庫、大西海子生態(tài)水庫和臺特瑪湖引水。具體見表2。

表2 塔里木河干流水量調(diào)度循環(huán)方案設(shè)計

2.2 方案對比

當前，塔河干流引水以滿足水資源區(qū)域性供給和改善水環(huán)境為主要目的，不同的水量調(diào)度方案和線路對水資源、區(qū)域水環(huán)境和防洪具有不同的影響。

方案一雖然引水量大幅增加，但是引水過程中，存在水資源大量流失的循環(huán)短路。另外在托什干河、庫瑪拉克河聯(lián)合排水時托什干河沿程水位降低將導(dǎo)致污水流入塔河干流，造成水量進一步流失，拉低河網(wǎng)整體水位，比其余方案河網(wǎng)水位平均降低0.11 m。根據(jù)《流域水循環(huán)質(zhì)量標準》，上游引水區(qū)域大部分河道COD濃度為18 mg/L～28 mg/L以內(nèi)，屬于Ⅳ類水，僅阿拉爾地區(qū)屬于Ⅴ類水。而NH3-N濃度含量多數(shù)河道均超過2 mg/L，整體屬于劣Ⅴ類水，污染嚴重。該方案聯(lián)合排水有助于水位降低從而形成水位差和水體循環(huán)線路，但是河道各斷面的COD和NH3-N濃度并無明顯降低趨勢，水體沒有很大改善。

方案二和常規(guī)調(diào)度水位較為接近，引排水量均有較大幅度增加，但其并未導(dǎo)致河網(wǎng)水位降低，表明，該方案雖然在本區(qū)域內(nèi)形成了水循環(huán)回路，但是卻導(dǎo)致水質(zhì)惡化，水環(huán)境明顯變差，主要是因為該方案下排水河道所排出的污染物相當有限。

方案三具有明顯的水質(zhì)改善效果，通過引水渠引水控制線進行控制，大部分污染物排走，該調(diào)度方案下水量比常規(guī)調(diào)度有所降低，本方案所形成的水量循環(huán)回路還能有效改善各斷面水質(zhì)狀況，即便是水環(huán)境污染較為嚴重的下游荒漠區(qū)域也形成了水體流動態(tài)勢。

根據(jù)式（1）和式（2）分別進行塔河干流不同調(diào)度方案下斷面流量水位以及水質(zhì)情況的計算，三個水量調(diào)度管理方案的水位、進出水量及水質(zhì)計算成果見表3和表4。

表3 不同調(diào)度方案的進出水量成果統(tǒng)計

表4 不同調(diào)度方案主要斷面COD、NH3-N濃度統(tǒng)計表 單位:mg/L

塔河干流水量調(diào)度管理以改善水環(huán)境為主要目的，所以水質(zhì)改善程度是評判的主要標準，河網(wǎng)水位和引排水量則是次要的評判因子。上述三個方案中，方案三通過加強干流控制有效防止了高片污水進入干流低片地勢低洼區(qū)域，使得干流水環(huán)境顯著改善，河網(wǎng)整體水位明顯提高[3]，所以以方案三為推薦方案。

3 結(jié)論

基于河網(wǎng)水流水質(zhì)耦合模型，求出不同引水方案的進出水量、水質(zhì)情況，再根據(jù)這些內(nèi)容對三個引水方案的分析和比選，分析的具體結(jié)果是:方案三為最優(yōu)引水方案。以帶動流域水環(huán)境改善為目的所進行的干流水資源優(yōu)化配置的引水方案的研究，對于改善水質(zhì)、提升流域水量管理調(diào)度水平具有積極意義。