南水北調(diào)中線工程突發(fā)水污染下游供水方式研究

摘 要：南水北調(diào)中線工程以明渠調(diào)水的形式進(jìn)行輸水，具有投資小、費(fèi)用低、流量大的優(yōu)勢，但敞開的輸水線路也帶來了水源保護(hù)難度大、容易引發(fā)突發(fā)水污染事故的風(fēng)險(xiǎn)。系統(tǒng)總結(jié)國內(nèi)外研究成果，以南水北調(diào)中線工程為研究區(qū)，以蘭河節(jié)制閘—潁河節(jié)制閘段發(fā)生水污染事件為例，重點(diǎn)開展突發(fā)水污染事故期事故段下游供水時(shí)間的研究，基于明渠調(diào)水的常規(guī)供水方式，采用常規(guī)和優(yōu)化兩種供水方式對事故段下游供水時(shí)間進(jìn)行模擬，提出了延長事故段下游供水時(shí)間的優(yōu)化分區(qū)供水方式。結(jié)果表明，采用優(yōu)化分區(qū)供水方式在面對渠池下游大流量分水情況下可以維持下游渠池長時(shí)間的供水穩(wěn)定，有效地延長下游不利渠池的供水時(shí)間，保障區(qū)域供水安全。

關(guān)鍵詞：南水北調(diào)中線工程；突發(fā)水污染；明渠調(diào)水；優(yōu)化分區(qū)；供水安全

中圖分類號：ＴＶ２１３．９ 文獻(xiàn)標(biāo)志碼：Ａｄｏｉ：１０．３９６９／ ｊ．ｉｓｓｎ．１０００－１３７９．２０２３．０４．０１１

引用格式：張乃豐，任紅蕾，林飛．南水北調(diào)中線工程突發(fā)水污染下游供水方式研究［Ｊ］．人民黃河，２０２３，４５（４）：６３－６９．

１ 引 言

明渠調(diào)水是大型調(diào)水工程中的常用輸水方式，以天然河道或人工渠道為主要輸水載體，具有投資小、運(yùn)行費(fèi)用低、輸水流量大的特點(diǎn)［１］ 。明渠調(diào)水工程中，水質(zhì)安全是調(diào)水工程發(fā)揮經(jīng)濟(jì)效益和社會效益的重要保障。但是近些年來突發(fā)水污染事故頻繁發(fā)生，不僅給環(huán)境帶來不可估量的影響，還給社會和經(jīng)濟(jì)發(fā)展帶來威脅，引發(fā)人們對明渠調(diào)水工程用水安全的擔(dān)憂［２］ 。突發(fā)水污染事故是指突然發(fā)生的且由人為因素或自然因素等引起大量污染物進(jìn)入水體事故［３］ ，從而導(dǎo)致水質(zhì)惡化，影響生活生產(chǎn)用水安全，造成巨大經(jīng)濟(jì)損失，并且危害生態(tài)環(huán)境、帶來惡劣的社會影響［４］ ?，F(xiàn)階段針對明渠調(diào)水工程突發(fā)水污染事故，主要通過構(gòu)建機(jī)理模型，預(yù)測分析突發(fā)水污染事故后污染物的遷移過程，結(jié)合渠道的水工建筑物進(jìn)行聯(lián)合調(diào)控，從而實(shí)現(xiàn)對突發(fā)水污染事故的應(yīng)急調(diào)控［５］ 。Ｇｒｉｆｏｕｌｌ等［６］ 在巴塞羅那港建立三維水質(zhì)模型，結(jié)合港內(nèi)的節(jié)制閘進(jìn)行調(diào)控，實(shí)現(xiàn)對突發(fā)水污染事件的應(yīng)急處理。Ｓａａｄａｔｐｏｕｒ 等［７］ 在伊拉姆水庫建立二維水質(zhì)模型，利用水庫內(nèi)的退水閘實(shí)現(xiàn)對污染水體的精準(zhǔn)退水。Ｈｏｕ等［８］ 建立可模擬突發(fā)水污染事故的實(shí)時(shí)動(dòng)態(tài)預(yù)警模型，模擬了新安江苯酚泄漏事故后苯酚在錢塘江的擴(kuò)散過程。吳輝明等［９］ 建立淮河干流一維水動(dòng)力－水質(zhì)模型，分析了在支流不同流量情況下干流突發(fā)水污染事件時(shí)污染物濃度隨時(shí)間、空間變化規(guī)律，據(jù)此并結(jié)合邊界退水閘可實(shí)現(xiàn)突發(fā)水污染的應(yīng)急調(diào)控。Ｔａｎｇ等［１０］ 基于ＭＩＫＥ１１ 建立一維水動(dòng)力－水質(zhì)模型，分析了南水北調(diào)工程中線干渠突發(fā)水污染事故后的污染物擴(kuò)散過程，研究３ 種流量工況下的４ 種污染物（磷肥、氰化物、石油、鉻）的擴(kuò)散規(guī)律，結(jié)合水工建筑物實(shí)現(xiàn)對４ 種污染物精準(zhǔn)控制。周超等［１１］ 建立南水北調(diào)中線邢臺和石家莊交界處至古運(yùn)河南段一維水質(zhì)模型，通過分析突發(fā)水污染事件后污染物的擴(kuò)散規(guī)律和退水閘的響應(yīng)時(shí)間，建立了污染物溯源應(yīng)急決策系統(tǒng)。王昭亮等［１２］ 選取淮河流域沙潁河下游及淮河干流作為研究區(qū)域，考慮沙潁河潁上閘的調(diào)度調(diào)控作用，建立以目標(biāo)節(jié)點(diǎn)幾種典型污染物濃度最小為優(yōu)化目標(biāo)、考慮閘壩利用和調(diào)度約束的優(yōu)化調(diào)度模型，并對比在同樣入河污染物荷載的情況下常規(guī)調(diào)度和優(yōu)化調(diào)度模型得到的調(diào)度結(jié)果，結(jié)果表明，經(jīng)過優(yōu)化后潁上閘的調(diào)度調(diào)控能夠有效地降低目標(biāo)節(jié)點(diǎn)污染物濃度，且污染物峰值濃度明顯降低。

縱觀國內(nèi)外眾多學(xué)者的研究成果，運(yùn)用模型模擬突發(fā)水污染事件污染物的遷移規(guī)律、實(shí)現(xiàn)污染渠池的精準(zhǔn)退水技術(shù)已十分完善，但對突發(fā)水污染事故段下游不同供水方式的供水時(shí)間研究較少。調(diào)水工程多服務(wù)于沿線城市的供水［１３］ ，當(dāng)渠池發(fā)生突發(fā)水污染事件時(shí)，維持事故段下游渠池的供水時(shí)間顯得尤為重要。故本文針對明渠調(diào)水工程渠池突發(fā)水污染事件時(shí)如何保障事故段下游供水時(shí)間進(jìn)行研究，提出基于不同分水等級的優(yōu)化分區(qū)方法，從而延長事故段下游渠池的供水時(shí)間。

２ 研究區(qū)域與方法

２．１ 研究區(qū)域概況

南水北調(diào)中線工程從漢江中上游的丹江口水庫引水，途經(jīng)河南、河北兩省，橫穿海河、黃河、淮河和長江四大流域［１４］ ，自流輸水將“南水”送至北京的團(tuán)城湖和天津的外環(huán)河［１５］ 。工程全程長１ ４３２ ｋｍ，沿線共布置６１ 座節(jié)制閘、９５ 座分水口、１ 座泵站以及倒虹吸、渡槽、無壓涵洞等眾多類型建筑物［１６－１７］ 。南水北調(diào)中線工程于２０１４ 年１２ 月１２ 日正式通水，迄今已安全輸水近３ ０００ ｄ，累計(jì)輸水超過４４１ 億ｍ３，使沿線６ ０００ 萬人從中獲益［１８］ 。本次研究的區(qū)域?yàn)槟纤闭{(diào)中線工程蘭河節(jié)制閘—北拒馬河節(jié)制閘段，區(qū)域內(nèi)共有４５ 座節(jié)制閘、５８ 座分水口、４１ 座退水口。

２．２ 研究方法

２．２．１ 水動(dòng)力模型構(gòu)建

南水北調(diào)中線工程總干渠被６１ 座節(jié)制閘劃分為６０ 個(gè)渠池，每一個(gè)渠池水面線受設(shè)計(jì)水深和當(dāng)前流量影響、也受相鄰渠池間的耦合影響［１９］ 。在推求渠池蓄水量時(shí)，假定每一個(gè)渠池的水面線為恒定的，然后對這個(gè)渠池進(jìn)行恒定流計(jì)算。

式中： Ｅｓｄ 為欲求斷面的斷面比能； Ｅｓｕ 為已知斷面的斷面比能； ΔＳ 為欲求斷面和已知斷面之間的間距； Ｊ為欲求斷面和已知斷面之間的水力坡度平均值。以已知斷面為控制斷面，可根據(jù)式（３）推導(dǎo)出欲求斷面。連續(xù)推求則水面線便可確定。

２．２．２ 分水口等級劃分

對于突發(fā)水污染事故下游段，在事故發(fā)生后由于進(jìn)口的供水被中斷，處于一種利用渠池內(nèi)蓄水量給分水口門供水的工況。因此，在事故發(fā)生后，給不同的分水口門持續(xù)供水，事故段下游須采用不同的應(yīng)急調(diào)控策略，首先應(yīng)對事故段下游的分水口門持續(xù)供水情況進(jìn)行等級劃分。根據(jù)沿線供水對象的重要程度，將分水口分為４ 個(gè)等級，依次對應(yīng)北京供水口門（Ⅰ級）、省會城市供水口門（Ⅱ級）、普通地級市供水口門（Ⅲ級）和其他城市供水口門（Ⅳ級）。受限于篇幅，本文只列舉Ⅰ級分水口門、Ⅱ級分水口門以及Ⅲ級分水口門的具體信息，見表１、表２ 和表３。

按照突發(fā)水污染事故期事故段下游分水口門供水對象的不同，可分為Ⅰ級分水口門供水、Ⅰ＋Ⅱ級分水口門供水、Ⅰ＋Ⅱ＋Ⅲ級分水口門供水、下游所有分水口門供水工況。Ⅰ級供水是將中線的水全部用于北京的城市供水；Ⅰ＋Ⅱ級供水主要保障北京和其他地區(qū)省會（包括天津）的供水；Ⅰ＋Ⅱ＋Ⅲ級供水主要保障北京、省會和地級市城市的供水。這種等級劃分可服務(wù)于調(diào)度人員，在不同的應(yīng)急工況以及供水矛盾情況下，快速合理地決策事故段下游的分水供給情況。

２．２．３ 優(yōu)化分區(qū)方法

明渠調(diào)水工程正常輸水狀態(tài)通常采用局部分區(qū)方法（見圖１）對渠池下游城市進(jìn)行供水，局部分區(qū)的供水方法是指將當(dāng)前渠池上游節(jié)制閘與下游臨近第１ 個(gè)分水口的節(jié)制閘作為１ 個(gè)供水分區(qū)，當(dāng)前渠池的水量僅供給當(dāng)前渠池或臨近下游分水口。局部分區(qū)在面向突發(fā)水污染事件時(shí)，上游長時(shí)間不能提供來水的情況下，就可能出現(xiàn)渠池的蓄水量不能滿足城市供水的問題。

面對明渠突發(fā)水污染事件，為延長不利渠池的供水時(shí)間，提出了一種通過合理分區(qū)從而延長不利渠池的持續(xù)供水時(shí)間的優(yōu)化分區(qū)方法（見圖２）。首先運(yùn)用水動(dòng)力模型計(jì)算研究區(qū)段中每１ 個(gè)渠池的蓄水量，利用式（４）計(jì)算當(dāng)前渠池向下游依次累加后的供水時(shí)間，再按照式（５）找出不利渠池；然后關(guān)閉不利渠池下游的節(jié)制閘，將不利渠池及以上所有的渠池當(dāng)作１ 個(gè)供水分區(qū)，不利渠池的下１ 個(gè)渠池將作為新的第１ 個(gè)渠池，重復(fù)以上的工作，重新識別不利渠池，直至識別出不利渠池為最后１ 個(gè)渠池，優(yōu)化分區(qū)工作完成。

式中： Ｔ１ 為第１ 個(gè)渠池供水時(shí)間； Ｔ１＋２ 為第１ 個(gè)渠池和第２ 個(gè)渠池累加后供水時(shí)間； Ｔ１＋２＋···＋Ｎ 為第１ 個(gè)渠池依次累加至第Ｎ 個(gè)渠池后供水時(shí)間； Ｖ１ 為第１ 個(gè)渠池蓄水量； Ｑ１ 為第１ 個(gè)渠池分水量； Ｖ２ 為第２ 個(gè)渠池蓄水量； Ｑ２ 為第２ 個(gè)渠池分水量； Ｖｉ 為第ｉ 個(gè)渠池蓄水量；Ｑｉ 為第ｉ 個(gè)渠池分水量。Ｔ１＋２＋…＋（Ｎ－１） ＜ Ｔ１＋２＋…＋Ｎ ＜ Ｔ１＋２＋…＋（Ｎ＋１） （５）式中：Ｔ１＋２＋…＋Ｎ 為從第１ 個(gè)渠池到下游第Ｎ 個(gè)渠池的累計(jì)供水時(shí)間。

上述方法計(jì)算的供水時(shí)間為下游所有渠池的平均供水時(shí)間，由于從第１ 個(gè)渠池累加至第Ｎ ＋ １ 個(gè)渠池的供水時(shí)間大于從第１ 個(gè)渠池累加至第Ｎ 個(gè)渠池的供水時(shí)間，因此從以下兩點(diǎn)進(jìn)行分析：① 從蓄水量的角度分析，當(dāng)累加到第Ｎ ＋ １ 個(gè)渠池時(shí)，整個(gè)系統(tǒng)的供水時(shí)間增加，說明第Ｎ ＋ １個(gè)渠池的蓄水量被平均分配到上游的渠池中，致使上游所有渠池的平均供水時(shí)間增加，因?yàn)橹芯€調(diào)水工程渠池的水量無法倒流，所以應(yīng)該將第Ｎ 個(gè)渠池下游的節(jié)制閘關(guān)閉；② 從分水量的角度分析，當(dāng)發(fā)生上述情況時(shí)，說明第Ｎ 個(gè)渠池的分水量較大，對整個(gè)供水系統(tǒng)影響較大，無法繼續(xù)向下游供水，此時(shí)可關(guān)閉第Ｎ 個(gè)渠池下游的節(jié)制閘，從而達(dá)到延長下游供水時(shí)間的效果。

３ 案例分析

渠池發(fā)生突發(fā)水污染事件時(shí)，發(fā)生污染的渠池被稱為事故段，渠池從發(fā)生污染至污染物從退水口全部退出這段時(shí)間稱之為事故期。當(dāng)發(fā)生突發(fā)水污染事件后事故段渠池下游節(jié)制閘關(guān)閉，將明渠調(diào)水工程分成事故段上游、事故段以及事故段下游，分段示意見圖３。

本文以南水北調(diào)中線干渠工程蘭河—潁河節(jié)制閘段發(fā)生突發(fā)水污染事件為例，將中線干渠所有的分水口按照重要等級分成Ⅰ、Ⅱ、Ⅲ、Ⅳ級。為了體現(xiàn)優(yōu)化分區(qū)適應(yīng)于不同的供水模式，當(dāng)發(fā)生突發(fā)水污染事件時(shí)，本文選用４ 種不同的供水模式來展示優(yōu)化分區(qū)相較于局部分區(qū)的供水方法可以延長不利渠池供水時(shí)間的效果。對于事故段，無論下游采用何種供水模式，事故段都應(yīng)當(dāng)關(guān)閉其下游的第１ 個(gè)節(jié)制閘（即潁河節(jié)制閘），使得污染團(tuán)在事故段全部排出，降低污染物擴(kuò)散的范圍；對于事故段下游，首先按照分級將供水方式分為Ⅰ級供水、Ⅰ＋Ⅱ級供水、Ⅰ＋Ⅱ＋Ⅲ級供水以及全線供水４ 種供水模式，再按照不同的等級串聯(lián)渠池，利用渠池自身的蓄水量，向相應(yīng)等級的供水口門進(jìn)行供水，從而生成４ 種供水方案。本文選擇流量工況為渠首（陶岔渠首工程）３５０ ｍ３ ／ ｓ，此工況下事故段下游渠池蓄水量見表４。

３．１ Ⅰ級供水

采用Ⅰ級供水方式目的是保障北京的用水安全，不考慮沿線其他省市的供水問題，所以沿線的全部分水口門都關(guān)閉，只打開北拒馬河節(jié)制閘用于向北京供水。由于Ⅰ級供水只有１ 個(gè)節(jié)制閘進(jìn)行分水，所以優(yōu)化分區(qū)和局部分區(qū)的分區(qū)方法一致，供水時(shí)間都為３２．８１ ｄ。

３．２ Ⅰ＋Ⅱ級供水

Ⅰ＋Ⅱ級供水方式通過打開北京和沿線省會城市的分水口，保障北京和沿線省會城市的供水。分水口門的流量采用實(shí)測分水流量，北拒馬河節(jié)制閘采用實(shí)際過閘流量。按照局部分區(qū)方法可將事故段下游渠池分為潁河—古運(yùn)河、古運(yùn)河—磁河、磁河—西黑山、西黑山—北拒馬河４ 個(gè)分區(qū)，分別計(jì)算每個(gè)分區(qū)的供水時(shí)間；按照優(yōu)化分區(qū)方法分別計(jì)算潁河—古運(yùn)河、潁河—磁河、潁河—西黑山、潁河—北拒馬河４ 個(gè)分區(qū)的供水時(shí)間。

分別用局部分區(qū)方法和優(yōu)化分區(qū)方法計(jì)算Ⅰ＋Ⅱ級供水方式的下游供水時(shí)間，計(jì)算結(jié)果見表５。

表５ 為采用兩種分區(qū)方法的供水時(shí)間，第１ 次分區(qū)，兩種方法都是將潁河—古運(yùn)河作為第１ 分區(qū)；第２次分區(qū)，由于優(yōu)化分區(qū)涵蓋潁河—古運(yùn)河的蓄水量，分區(qū)的蓄水量遠(yuǎn)遠(yuǎn)大于該分區(qū)內(nèi)分水口的分水量，采用優(yōu)化分區(qū)方法較局部分區(qū)方法延長下游渠池供水時(shí)間１１．８５ ｄ；第３ 次分區(qū)，從局部分區(qū)方法來看，由于磁河—西黑山分區(qū)只串聯(lián)了７ 個(gè)渠池，但是西黑山分水口分水流量較大，所以局部分區(qū)的供水時(shí)間為８．１７ ｄ，但是采用優(yōu)化分區(qū)方法，串聯(lián)了４０ 個(gè)渠池的蓄水量，大大減緩了西黑山大流量分水對整個(gè)分區(qū)蓄水量的影響，所以供水時(shí)間遠(yuǎn)大于采用局部分區(qū)方法的供水時(shí)間；第４ 次分區(qū)，由于北拒馬河節(jié)制閘分水流量最大，所以采用局部分區(qū)方法只能維持１．３４ ｄ，但是采用優(yōu)化分區(qū)方法供水時(shí)間可達(dá)３７．２１ ｄ，遠(yuǎn)大于局部分區(qū)方法的供水時(shí)間。相較于第３ 次優(yōu)化分區(qū)的供水時(shí)間，第４ 次優(yōu)化分區(qū)的供水時(shí)間減少了０．３６ ｄ，說明優(yōu)化分區(qū)方法在大流量分水情況下仍然能穩(wěn)定維持下游的供水時(shí)間。

通過優(yōu)化分區(qū)方法識別出事故段下游不利渠池為西黑山—北拒馬河節(jié)制閘，根據(jù)計(jì)算，得到選用優(yōu)化分區(qū)方法可以使北拒馬河節(jié)制閘持續(xù)供水３７．２１ ｄ，但是采用局部分區(qū)方法僅使北拒馬河節(jié)制閘持續(xù)供水１．３４ｄ，所以通過優(yōu)化分區(qū)的供水方法可延長不利渠池的供水時(shí)間為３５．８７ ｄ。

３．３ Ⅰ＋Ⅱ＋Ⅲ級供水

采用Ⅰ＋Ⅱ＋Ⅲ級供水方式的目的是保障全線所有地級市的供水安全，打開所有地級市的分水口門和北拒馬河節(jié)制閘，分水口門的流量采用實(shí)際分水流量，節(jié)制閘的流量采用實(shí)際過閘流量。

依據(jù)上述局部分區(qū)和優(yōu)化分區(qū)兩種方法，渠池分區(qū)與渠池的供水時(shí)間見表６。

表６ 為在Ⅰ＋Ⅱ＋Ⅲ級供水工況下局部分區(qū)和優(yōu)化分區(qū)的供水時(shí)間計(jì)算結(jié)果，可發(fā)現(xiàn)分區(qū)下游段為古運(yùn)河節(jié)制閘、西黑山節(jié)制閘、北拒馬河節(jié)制閘的分區(qū)供水時(shí)間都特別短，原因是３ 個(gè)渠池的分水流量較大，遠(yuǎn)大于事故段下游其他渠池的分水流量，對渠池的蓄水量影響較大。由表６ 局部分區(qū)供水時(shí)間可知，在經(jīng)歷大流量分水的情況下，七里河—古運(yùn)河渠池、崗頭隧洞—西黑山渠池、西黑山—北拒馬河渠池供水時(shí)間分別為１４．９５、０．６４、１．３４ ｄ，而優(yōu)化分區(qū)方法在同樣經(jīng)歷了大流量分水的情況下，西黑山節(jié)制閘可以持續(xù)供水３３．６８ ｄ，北拒馬河節(jié)制閘可以持續(xù)供水１７．６３ ｄ。優(yōu)化分區(qū)方法相較于局部分區(qū)方法，同樣在大流量分水情況下可以維持更長時(shí)間的穩(wěn)定供水。

根據(jù)優(yōu)化分區(qū)方法，可以識別出不利渠池潁河—古運(yùn)河段供水時(shí)間為５０．８０ ｄ、潁河—北拒馬河段供水時(shí)間為１７．６３ ｄ。由于優(yōu)化分區(qū)方法計(jì)算的是分區(qū)內(nèi)各子渠池的供水時(shí)間，因此在優(yōu)化分區(qū)供水方法下，上述渠段內(nèi)對應(yīng)的子渠池分別為七里河—古運(yùn)河、西黑山—北拒馬河，但在同樣條件下，采用局部分區(qū)方法計(jì)算出七里河—古運(yùn)河渠池、西黑山—北拒馬河渠池的持續(xù)供水時(shí)間分別為１４．９５、１．３４ ｄ，通過對比分析，以上兩個(gè)不利渠池的供水時(shí)間分別延長了３５． ８５、１６．２９ ｄ。

３．４ 全部口門供水

采用全部分水口門供水方式的目的是保障全線所有沿線城市供水，打開所有的分水口門和北拒馬河節(jié)制閘，分水口門的流量采用實(shí)際分水流量，北拒馬河節(jié)制閘的流量采用實(shí)際過閘流量。分別用局部分區(qū)方法和優(yōu)化分區(qū)方法計(jì)算全部口門供水方式的供水時(shí)間，計(jì)算結(jié)果如圖４、圖５ 所示。

圖４ 表示局部分區(qū)方法事故段下游供水時(shí)間，從圖中可以看出局部分區(qū)下游為西黑山、北拒馬河節(jié)制閘的分區(qū)供水時(shí)間較短，分別為０．６４、０．６９ ｄ。圖５ 為優(yōu)化分區(qū)方法事故段下游供水時(shí)間，下游為西黑山節(jié)制閘的分區(qū)，雖然相較于Ⅰ＋Ⅱ＋Ⅲ供水方式減少了１３．０１ ｄ，但是也達(dá)到了２０．６７ ｄ，相較于局部分區(qū)方法延長供水時(shí)間為２０ ｄ 左右；同理，下游為北拒馬河節(jié)制閘的分區(qū)，相較于較于Ⅰ＋Ⅱ＋Ⅲ供水方式減少了４．３７ ｄ，但是比局部分區(qū)方法多向城市供水１３ ｄ。根據(jù)優(yōu)化分區(qū)方法，可以識別出不利渠池潁河—須水河、潁河—穿黃、潁河—淇河、潁河—沁河、潁河—南沙河、潁河—古運(yùn)河、潁河—放水河、潁河—西黑山、潁河—北拒馬河的供水時(shí)間分別為１７． ４８、１９． ７４、２９．７０、１５．８７、３３．２７、２６．５６、２６．４０、２０．６７、１３．２６ ｄ；根據(jù)優(yōu)化分區(qū)方法的原理，上述渠段對應(yīng)的子渠池分別為金水河—須水河、索河—穿黃、香泉河—淇河、牯牛河—沁河、洺河—南沙河、洨河—古運(yùn)河、漠道溝—放水河、崗頭隧洞—西黑山、北易水—北拒馬河，但采用局部分區(qū)方法以上各不利渠池的供水時(shí)間分別為３．８０、０．６９、２５．８７、１５．８５、２５．６４、２．４６、１７．７５、０．６４、０．６９ｄ，通過對比分析，優(yōu)化分區(qū)方法可以大大延長不利渠池的供水時(shí)間。

４ 結(jié) 論

國內(nèi)外為應(yīng)對突發(fā)水污染事件，開展了突發(fā)水污染事件事故段下游不同供水方式下供水時(shí)間的研究。本文基于不同分水等級，通過對南水北調(diào)中線干渠蘭河—潁河突發(fā)水污染事件進(jìn)行研究，提出了優(yōu)化分區(qū)的供水方法，相比較中線傳統(tǒng)局部分區(qū)的供水方法可以有效地延長下游供水時(shí)間。因此，通過４ 種不同工況的驗(yàn)證，在下游分水口較大流量分水的情況下，優(yōu)化分區(qū)的供水方法相比較于傳統(tǒng)局部分區(qū)的供水方法可以較好地緩解大流量分水對渠池蓄水量的影響，保障大流量分水情況下渠道下游供水的穩(wěn)定。優(yōu)化分區(qū)的供水方法識別出不利渠池為墳莊河—北拒馬河渠段，相較于局部分區(qū)的供水方法，優(yōu)化分區(qū)大大地延長不利渠池的供水時(shí)間。

參考文獻(xiàn)：

［１］ 鄧銘江，王遠(yuǎn)超，李江，等．沙漠地區(qū)大型供水工程建設(shè)的關(guān)鍵技術(shù)與方法［Ｊ］． 巖石力學(xué)與工程學(xué)報(bào)，２０１７，３６（１２）：３０８５－３０９４．

［２］ 蘭才有，儀修堂．南水北調(diào)工程部分線段應(yīng)采用管道輸水［Ｊ］．南水北調(diào)與水利科技，２００４，２（１）：２２－２４．

［３］ ＹＡＮＧ Ｈ Ｄ，ＳＨＡＯ Ｄ Ｇ，ＬＩＵ Ｂ Ｙ，ｅｔ ａｌ．Ｍｕｌｔｉ?Ｐｏｉｎｔ ＳｏｕｒｃｅＩｄｅｎｔｉｆｉｃａｔｉｏｎ ｏｆ Ｓｕｄｄｅｎ Ｗａｔｅｒ Ｐｏｌｌｕｔｉｏｎ Ａｃｃｉｄｅｎｔｓ ｉｎ ＳｕｒｆａｃｅＷａｔｅｒｓ Ｂａｓｅｄ ｏｎ Ｄｉｆｆｅｒｅｎｔｉａｌ Ｅｖｏｌｕｔｉｏｎ ａｎｄ Ｍｅｔｒｏｐｏｌｉｓ?Ｈａｓｔｉｎｇ Ｍａｒｋｏｖ Ｃｈａｉｎ Ｍｏｎｔｅ Ｃａｒｌｏ［Ｊ］．Ｓｔｏｃｋ Ｅｎｖｉｒｏｎ Ｒｅ?ｓｅａｒｃｈ ａｎｄ Ｒｉｓｋ Ａｓｓｅｓｓ，２０１６，３０（２）：５０７－５２２．

［４］ 王浩，賈仰文．變化中的流域“自然－社會”二元水循環(huán)理論與研究方法［Ｊ］．水利學(xué)報(bào)，２０１６，４７（１０）：１２１９－１２２６．

［５］ 鄭和震．南水北調(diào)中線干渠突發(fā)水污染擴(kuò)散預(yù)測與應(yīng)急調(diào)度［Ｄ］．杭州：浙江大學(xué)，２０１８：３２－４２．

［６］ ＧＲＩＦＯＵＬＬ Ｍ，ＪＯＲＤà Ｇ，ＥＳＰＩＮＯ Ｍ，ｅｔ ａｌ．Ａ ＭａｎａｇｅｍｅｎｔＳｙｓｔｅｍ ｆｏｒ Ａｃｃｉｄｅｎｔａｌ Ｗａｔｅｒ Ｐｏｌｌｕｔｉｏｎ Ｒｉｓｋ ｉｎ ａ Ｈａｒｂｏｕｒ：Ｔｈｅ Ｂａｒｃｅｌｏｎａ Ｃａｓｅ Ｓｔｕｄｙ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ｍａｒｉｎｅ Ｓｙｓｔｅｍｓ，２０１１，８８（１）： ６０－７３．

［７］ ＳＡＡＤＡＴＰＯＵＲ Ｍ，ＡＦＳＨＡＲ Ａ．Ｍｕｌｔｉ Ｏｂｊｅｃｔｉｖｅ Ｓｉｍｕｌａｔｉｏｎ? Ｏｐｔｉｍｉｚａｔｉｏｎ Ａｐｐｒｏａｃｈ ｉｎ Ｐｏｌｌｕｔｉｏｎ Ｓｐｉｌｌ Ｒｅｓｐｏｎｓｅ Ｍａｎａｇｅ? ｍｅｎｔ Ｍｏｄｅｌ ｉｎ Ｒｅｓｅｒｖｏｉｒｓ ［ Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｗａｔｅｒ ＲｅｓｏｕｒｃｅＭａｎａｇｅｍｅｎｔ，２０１３，２７：１８５１－１８６５．

［８］ ＨＯＵ Ｄ Ｂ，ＧＥ Ｘ Ｆ，ＨＵＡＮＧ Ｐ Ｊ，ｅｔ ａｌ．Ａ Ｒｅａｌ?Ｔｉｍｅ ＤｙｎａｍｉｃＥａｒｌｙ?Ｗａｒｎｉｎｇ Ｍｏｄｅｌ Ｂａｓｅｄ ｏｎ Ｕｎｃｅｒｔａｉｎｔｙ Ａｎａｌｙｓｉｓ ａｎｄＲｉｓｋ Ａｓｓｅｓｓｍｅｎｔ ｆｏｒ Ｓｕｄｄｅｎ Ｗａｔｅｒ Ｐｏｌｌｕｔｉｏｎ Ａｃｃｉｄｅｎｔｓ［Ｊ］．Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｓｃｉｅｎｃｅ ａｎｄ Ｐｏｌｌｕｔｉｏｎ Ｒｅｓｅａｒｃｈ Ｉｎｔｅｒｎａｔｉｏｎａｌ，２０１４，２１（１４）：８８７８－８８９２．

［９］ 吳輝明，雷曉輝，廖衛(wèi)紅，等．淮河干流突發(fā)性水污染事故預(yù)測模擬研究［Ｊ］．人民黃河，２０１６，３８（１）：７５－７８，８４．

［１０］ ＴＡＮＧ Ｃ Ｈ，ＹＩ Ｙ Ｊ，ＹＡＮＧ Ｚ Ｆ，ｅｔ ａｌ．Ｗａｔｅｒ Ｐｏｌｌｕｔｉｏｎ ＲｉｓｋＳｉｍｕｌａｔｉｏｎ ａｎｄ Ｐｒｅｄｉｃｔｉｏｎ ｉｎ ｔｈｅ Ｍａｉｎ Ｃａｎａｌ ｏｆ ｔｈｅ Ｓｏｕｔｈ?ｔｏ?Ｎｏｒｔｈ Ｗａｔｅｒ Ｔｒａｎｓｆｅｒ Ｐｒｏｊｅｃｔ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ｈｙｄｒｏｌｏｇｙ，２０１４，５１９：２１１１－２１２０．

［１１］ 周超，陳政，蔣婷婷，等．突發(fā)水污染事故污染云團(tuán)快速追蹤實(shí)驗(yàn)研究［Ｊ］．水資源與水工程學(xué)報(bào)，２０１４，２５（２）：２００－２０５．

［１２］ 王昭亮，高仕春，艾澤．閘壩對河流水質(zhì)的調(diào)控作用初步分析［Ｊ］．水利科技與經(jīng)濟(jì)，２０１０，１６（１２）：１３３９－１３４０，１３４７．

［１３］ 黃會勇．南水北調(diào)中線總干渠水量調(diào)度模型研究及系統(tǒng)開發(fā)［Ｄ］．北京：中國水利水電科學(xué)研究院，２０１３：１９－３７．

［１４］ 王浩，王建華．中國水資源與可持續(xù)發(fā)展［Ｊ］．中國科學(xué)院院刊，２０１２，２７（３）：３５２－３５８，３３１．

［１５］ 蔣云鐘，趙紅莉，董延軍，等．南水北調(diào)中線水資源調(diào)度關(guān)鍵技術(shù)研究［Ｊ］．南水北調(diào)與水利科技，２００７，５（４）：１－５．

［１６］ 王浩，雷曉輝，尚毅梓．南水北調(diào)中線工程智能調(diào)控與應(yīng)急調(diào)度關(guān)鍵技術(shù)［Ｊ］．南水北調(diào)與水利科技，２０１７，１５（２）：１－８．

［１７］ 王頔．“南水北調(diào)”的工程學(xué)意義［Ｊ］．河南水利與南水北調(diào)，２０１６，４５（８）：２７－２８，４９．

［１８］ 陳翔．南水北調(diào)中線干線工程應(yīng)急調(diào)控與應(yīng)急響應(yīng)系統(tǒng)研究［Ｄ］．北京：中國水利水電科學(xué)研究院，２０１５：７７－７８．

［１９］ 張大偉．南水北調(diào)中線干線水質(zhì)水量聯(lián)合調(diào)控關(guān)鍵技術(shù)研究［Ｄ］．上海：東華大學(xué)，２０１３：２３－２７．

【責(zé)任編輯 簡 群】