王家彪　雷曉輝　王浩

摘要：天然河流水污染事件頻發(fā)，如何運用水庫應(yīng)急調(diào)度降低污染損失具有重要意義。以突發(fā)水污染事件應(yīng)急處置為目標，分別構(gòu)建了水污染溯源、濃度預測和水庫應(yīng)急調(diào)度模型，并通過數(shù)值模擬的方式對水庫調(diào)度方案進行優(yōu)選，最終建立了一套完善的水庫應(yīng)急調(diào)度技術(shù)體系。將技術(shù)運用于賀江2013年水污染事件應(yīng)急處置分析，結(jié)果表明，模型重構(gòu)的監(jiān)測斷面污染物濃度過程與實測過程較為接近，優(yōu)選的水庫調(diào)度方案及其處置效果也與實際情況基本吻合。結(jié)果說明本文建立的水庫應(yīng)急調(diào)度關(guān)鍵技術(shù)體系可運用于指導河流突發(fā)水污染事件應(yīng)急處置。

關(guān)鍵詞：突發(fā)水污染；水庫調(diào)度；應(yīng)急；溯源

中圖分類號：X522文獻標志碼：A文章編號：

16721683（2018）02000106

Abstract：

Water pollution incidents occur frequently in the natural river systems.How to operate the reservoirs on rivers to minimize the pollution loss is surely of great significance to the emergency disposal.In this paper，we established several corresponding models for the sequential processes of source identification，prediction of contaminant concentration，and reservoir operation.Using numerical simulation method，we optimized the reservoir operation scheme.Eventually，we proposed a relatively complete technical system for emergency disposal and then applied it to the case of Hejiang pollution incident in 2013.Results showed that the simulated time series of contaminant concentration were similar to the measured series.The optimized scheme and its disposal result generally fit the actual situation.The results indicated that the proposed technical system of emergency disposal with reservoir operation is feasible and efficient for real sudden water pollution incidents in natural rivers.

Key words：sudden water pollution；reservoir operation；emergency；source identification

突發(fā)水污染事件具有不確定性、擴散性、危害性、處理艱巨性和影響長期性等特點，且難以從根本上杜絕，其污染物排放也無固定途徑[13]。對于突發(fā)水污染的應(yīng)急處置通?？煞譃楣こ绦源胧┖头枪こ绦源胧﹥纱箢怺4]。在天然河流中，若采用工程措施進行處置（吸附、混凝等），不僅投入成本大，而且一[JP2]些化學藥物的運用還存在負面影響。而通過水庫應(yīng)急調(diào)度的方式不僅能稀釋河流污染物濃度，還可加快污染團的運移、弱化擴散作用，有[JP]效縮減水污染影響范圍。如2005年廣東北江鎘污染事件中采用了“加大上游水庫排量以稀釋水體污染物、利用人工小洪峰加快污染物運移到下游處置區(qū)”的應(yīng)急措施有效地控制了污染事故的惡化。2012年廣西龍江發(fā)生的鎘污染事件中，通過調(diào)度柳江上游大埔、麻石等水庫水量進行沖污稀釋，保證了下游柳州市供水安全。

關(guān)于水庫應(yīng)急調(diào)度的研究一直是流域管理的熱點和難點。蘇友華[5]研究了崇左市各縣區(qū)突發(fā)水污染時如何利用上游水庫進行應(yīng)急調(diào)度的調(diào)水方案與實施辦法。辛小康等[6]探討了三峽水庫應(yīng)急調(diào)度措施對長江宜昌段水污染事故處置的有效性和可行性，研究表明水庫調(diào)度對瞬排型水污染事故的處置作用明顯。陶亞等[7]探討了包括工程應(yīng)急調(diào)度、吸附攔截等在內(nèi)的多種污染物應(yīng)急處置措施其應(yīng)用原理和處置效果。余真真等[8]研究了小浪底水庫應(yīng)急調(diào)度對下游水污染事故的處置情況，結(jié)果表明水庫應(yīng)急調(diào)度降低了下游一定范圍內(nèi)的污染程度。

盡管如此，運用水庫調(diào)度進行河流突發(fā)水污染應(yīng)急處置仍存在諸多問題[9]，包括水污染事件信息掌握不足、水庫調(diào)蓄功能發(fā)揮不充分、應(yīng)急處置理論和技術(shù)不完善等。本文以天然河流突發(fā)水污染應(yīng)急處置為目標，研究如何運用水庫進行合理調(diào)度以實現(xiàn)突發(fā)水污染的快速處置、最大可能降低污染損失，完善水庫應(yīng)急調(diào)度理論和技術(shù)，并通過實例分析水庫應(yīng)急調(diào)度的可行性，以期為我國河流突發(fā)水污染應(yīng)急處置提供理論指導和技術(shù)借鑒。

1基本原理與模型構(gòu)建

運用水庫調(diào)度進行突發(fā)水污染應(yīng)急處置主要通過污染物溯源、濃度預測和水庫調(diào)度方案優(yōu)選三個關(guān)鍵技術(shù)環(huán)節(jié)來實現(xiàn)。本文對這三個環(huán)節(jié)分別構(gòu)建相應(yīng)的數(shù)值計算模型。

[BT3]1.1河道水流計算模型

河道中污染物濃度預測離不開水流的計算。在重力、水壓力、沿程阻力和慣性力作用下，由動量守恒原理推導水流運動方程見式（1），并在不可壓縮假定下得到水流連續(xù)性方程見式（2），方程組可采用四點偏心隱格式（Preimann格式）差分方法求解[10]。

[SX（]A[]t[SX）]+[SX（]Q[]x[SX）]=q[JY]（1）

[SX（]Q[]t[SX）]+[SX（][]x[SX）][JB（（][SX（]Q2[]A[SX）][JB））]+gA[SX（]Z[]x[SX）]+g[SX（]n2|u|[]R4/3[SX）]Q=0[JY]（2）

式中：x為沿河道距離（m）；t為時間（s）；A為斷面過流面積（m2）；Q為斷面過流流量（m3/s）；u為斷面平均流速（m/s）；q為支流入流流量（出流為負） （m3/（s·m）；Z為河道水位（m）；g為重力加速度（m/s2）；n為曼寧糙率系數(shù)；R為水力半徑（m）。

[BT3]1.2污染物濃度預測模型

河道中污染物的遷移擴散是三維非穩(wěn)態(tài)問題，其過程不僅與水流變化有關(guān)，還受水溫、風速等因素影響。實際計算中，模擬水質(zhì)全過程并不可行，因此在構(gòu)建水質(zhì)模型時可對問題進行適當簡化。當河道中污染物沿水平和垂向混合時間短、能快速混合至橫向均勻時，近似認為污染物只隨水流沿縱向發(fā)生對流擴散?？紤]一級反應(yīng)下的污染物對流擴散方程見式（3）。

[SX（]AC（x，t）[]t[SX）]+[SX（]QC（x，t）[]x[SX）]=[SX（][]x[SX）][JB（（]DA[SX（]C（x，t）[]x[SX）][JB））]-[WTB1X]k[WTBX]AC（x，t）+∑S[JY]（3）

[JP3]式中：C（x，t）為所計算的斷面時刻污染物濃度（mg/L）；[JP+1]D為包括彌散作用在內(nèi)的縱向擴散系數(shù)[JP]（m2/s）；[WTB1X]k[WTBZ]為反應(yīng)速率常數(shù)（s1）；[JP]S為源匯項（g/（m·s），當河道中有污染物匯入時，此項不為0；其余同式（1）和（2）。

在水流計算基礎(chǔ)上，方程（3）采用差分方法進行求解。特別的，在均勻紊流條件下可求得方程（3）解析解[11]，見式（4）。

C（x，t）=[SX（]M0[]A[KF（]4πDt[KF）][SX）]exp[JB（（]-[SX（]（x-ut）2[]4Dt[SX）]-[WTB1X]k[WTBX]t[JB））][JY]（4）

式中：u為斷面平均流速（m/s）；m0為進入水體污染物強度（kg）。

[JP3]在水污染應(yīng)急調(diào)度過程中，可根據(jù)式（4）對不同調(diào)水方案下的河道污染物濃度變化情況進行快速預測。[JP]

[BT3]1.3污染物溯源模型

突發(fā)水污染溯源技術(shù)通過研究污染物在河渠中遷移轉(zhuǎn)化規(guī)律，依據(jù)所觀測的污染物濃度過程推測出污染物排放位置、排放時間以及排放強度。作為污染物濃度預測的反問題[12]，污染物溯源具有非線性和不適定性的特征[13]，計算內(nèi)容更為復雜。本次研究考慮了一種基于概率密度函數(shù)的溯源方法，方法通過對污染物逆向位置概率密度與正向濃度之間關(guān)系進行回歸分析，構(gòu)建一個以污染物排放位置、排放時間和排放強度為參數(shù)的優(yōu)化模型。

[JP3]由∫[DD（][]x[DD）]C（x，t）dx=M0，可對C（x，t）歸一化，見下式：[JP]

c（x，t）=C（x，t）/M0[JY]（5）

式中：c（x，t）為對應(yīng)于C（x，t）的單位質(zhì)量污染物濃度值，具有m1的量綱，表示了污染物質(zhì)t時刻出現(xiàn)在x斷面的概率。

[JP2]在河道中，以P（x，t′）表示由觀測斷面xd判定的t′時刻污染源在xs處的概率（即污染物質(zhì)由xs斷面經(jīng)時間td-t′輸運到xd斷面的概率），則P（xs，t′）滿足對流擴散方程（3）的伴隨狀態(tài)方程以及歸一化條件[12]（P（x，t′）也具有m-1的量綱），見式（6）和式（7）[JP]。

-[SX（]P（xs，t′）[]t[SX）]+[SX（]（QP（xs，t′）/A）[]x[SX）]+D[SX（]2P（xs，t′）[]2x[SX）]=0[JY]（6）

P（xd，td）=1[JY]（7）

式中：t′為逆向計算時間點；td為污染物濃度觀測時間點（s）；式（7）表示的是：污染物質(zhì)未發(fā)生輸運（t′=td）而出現(xiàn)在觀測斷面時，污染源只能是在觀測斷面處。

類似于濃度預測，在均勻紊流條件下，可得到方程（6）解析解，如下：

P（x，t′）=[SX（]1[][KF（]4πD（td-t′）[KF）][SX）]·exp[JB（（]-[SX（]（xd-xs-u（td-t′））2[]4D（td-t′）[SX）][JB））][JY]（8）

結(jié)合式（5），分別比較式（3）和（6）與式（4）和（8）可看出，P（xs，t′）與c（x，t）形式完全一致。事實上，P（xs，t′）與c（x，t）兩者關(guān)系都可由圖1進行確定，圖中箭頭表示輸運（移流項）方向。

[JP2]圖1中，（a）和（b）表示初始時間t0和經(jīng)過Δt時間后的污染物濃度分布，（c）和（d）表示初始監(jiān)測時間td以及往前追溯Δt的位置概率密度分布。當t-t0=td-t′時，污染物質(zhì)從源x0經(jīng)時間t-t0運動到斷面xd處的概率c（xd，t），與觀測者位于斷面xd處判斷污染物由斷面x0[JP]經(jīng)時間td-t′運動到斷面xd處的概率相等。即當t-t0=td-t′時，式（9）成立。

P（x0，t′）=c（xd，t）[JY]（9）

由圖1和式（9）可看出，正向濃度過程與逆向位置概率過程具有高度耦合性，兩者除了計算時間方向相反外，其余完全一致。因此可基于這種耦合關(guān)系由P（xs，t′）計算c（xd，t），構(gòu)建優(yōu)化模型實現(xiàn)溯源計算。

首先以計算濃度和觀測濃度誤差最小確定出模型目標函數(shù)如下[1416]：

min∑（M0·Pi（x0，t0）-Ci）2）[JY]（10）

式中：x0、t0和M0即為所求污染物排放位置、時間和強度三參數(shù)；Ci代表觀測系列濃度值（mg/L）。

約束條件x0、t0和M0的取值范圍由先驗信息給定，一般是現(xiàn)場調(diào)查或由已有資料分析得出，見式（11）和至（13）：

x0min≤x0≤x0min[JY]（11）

t0min≤t0≤t0min[JY]（12）

M0min≤M0≤M0min[JY]（13）

通過求解上述優(yōu)化模型，可實現(xiàn)污染物的快速溯源。

[BT3]1.4水庫應(yīng)急調(diào)度模型

基于水庫調(diào)度的突發(fā)水污染應(yīng)急處置其本質(zhì)上是水庫出庫流量的分配問題。若以水庫水資源量損失最小和處置歷時最短為目標進行調(diào)度，則目標函數(shù)如下：

G1=min[JB（{]∑[DD（]n[]i=1[DD）]QiT[JB）}][JY]（14）

G2=min[JB（{]T[JB）}][JY]（15）

式中：G1為目標調(diào)水量（m3）；G2為目標調(diào)水歷時（h）；n為參與應(yīng)急調(diào)度的水庫數(shù)目；Qi為流域第i個水庫出庫流量（m3/s）；T為調(diào)水時間。

水庫調(diào)度過程中，不僅受水庫可用水量和河道水流演進的約束，還要求控制斷面處水質(zhì)必須達標。

水庫供水量約束：QiT≤Vimax[JY]（16）

[JP4]水流水質(zhì)聯(lián)動約束：C[WTB1X]k[WTBX]（t）=f（x0，t0，M，Q1，…，Qn）[JY]（17）[JP]

控制斷面水質(zhì)濃度約束：C[WTB1X]k[WTBX]（t）≤C[WTB1X]k[WTBZ]max[JY]（18）

式中：Vimax表示第i個水庫最大可供水量（m3）；[WTB1X]k[WTBX]為流域水質(zhì)控制斷面編號，取值1，2，…；C[WTB1X]k[WTBX]（t）表示第[WTB1X]k[WTBZ]個控制斷面處污染物計算濃度值（mg/L）；f（x0，t0，M，Q1，…，Qn，t，[WTB1X]k[WTBX]）表示水流水質(zhì)聯(lián)動約束下第[WTB1X]k[WTBX]個控制斷面處的濃度過程，是濃度值與污染源參數(shù)（x0，t0，M）及調(diào)度水量Q1，…，Qn和變量t、[WTB1X]k[WTBX]的一種抽象映射；C[WTB1X]k[WTBZ]max表示第[WTB1X]k[WTBZ]個控制斷面處允許的最大污染物濃度值（mg/L）。

由此建立的優(yōu)化調(diào)度模型其約束條件都具有非線性，其中式（17）中水流水質(zhì)聯(lián)動約束很難通過具體表達式進行簡化。因此，模型可采用遺傳算法（GA）或微分進化算法（DEA）進行求解[14，17]。盡管如此，流域多水庫調(diào)度時通常難以找出或無最優(yōu)解，此時可通過多組方案比選的方式尋找滿意解[18]，實現(xiàn)水庫調(diào)度方案的優(yōu)選。

[BT2]2關(guān)鍵技術(shù)體系

流域突發(fā)水污染后，需第一時間啟動應(yīng)急監(jiān)測方案，對污染事件進行診斷。在確定污染物超標河段后，根據(jù)監(jiān)測濃度和初步診斷結(jié)果，對水污染進行快速追蹤溯源，并依據(jù)溯源結(jié)果到現(xiàn)場進行排查確認。然后，運用水污染快速預測模型預估污染可能波及的河段和控制斷面處污染濃物度變化過程，對整個污染事件進行重構(gòu)。最后，根據(jù)水污染應(yīng)急處置目標擬定多組水庫調(diào)度方案并進行模擬計算，對不同方案的處置效果進行對比分析。在分析調(diào)度方案合理性時，還需對方案的可能影響進行評價，如發(fā)電效益損失。另外，在水污染應(yīng)急處置中，流域啟動水量調(diào)度時一般會配合以工程方法，如攔截吸附、混凝沉降等，在應(yīng)急調(diào)度方案擬定時應(yīng)予以考慮。

分別將所構(gòu)建的計算模型運用于突發(fā)水污染應(yīng)急處置各個環(huán)節(jié)，構(gòu)建出一套完整的水庫應(yīng)急調(diào)度技術(shù)體系，并確立基于水庫調(diào)度的突發(fā)水污染應(yīng)急處置實現(xiàn)流程如圖2所示。由于應(yīng)急調(diào)度需實時決策，因此實現(xiàn)流程中存在動態(tài)調(diào)整與反饋修正的閉環(huán)過程，即根據(jù)控制斷面處計算的濃度過程不斷調(diào)整應(yīng)急調(diào)水方案，直到斷面處濃度達標。

3應(yīng)用實例

2013年7月，廣西賀江流域發(fā)生鎘污染事件，賀江干流近110 km河段受到污染波及，干流上合面獅水庫整個庫區(qū)水體被污染，應(yīng)急處理經(jīng)費高達幾千萬元。本文以賀江2013年水污染事件為應(yīng)用實例，分析如何應(yīng)用所構(gòu)建的水庫應(yīng)急調(diào)度技術(shù)體系實現(xiàn)河流突發(fā)水污染的應(yīng)急處置。

3.1污染物溯源分析

根據(jù)實測資料，2013年7月7日6時監(jiān)測到賀江支流馬尾河河口處鎘濃度嚴重超標，同時在賀江干流廈島斷面處（馬尾河匯入口上游約2 km）并未檢出鎘超標情況。事發(fā)時具體監(jiān)測斷面布設(shè)見圖3。

由圖4可看出，通過污染物溯源模型和濃度預測模型還原的馬尾河口污染物濃度過程與實際過程基本一致，無論是濃度峰值還是峰現(xiàn)時間都很接近。盡管計算的濃度過程未能反映出實測過程中的第二峰值，但兩者具有相同的變化趨勢?？紤]到應(yīng)急調(diào)度方案的擬定主要受峰值和峰現(xiàn)時間以及濃度變化趨勢影響，模型溯源結(jié)果可用于調(diào)度方案分析。

3.2應(yīng)急調(diào)度方案優(yōu)選

根據(jù)《賀江重金屬污染處置應(yīng)急調(diào)水方案》，擬定合面獅水庫按300 m3/s的流量下泄污水，并通過聯(lián)合調(diào)用龜石水庫2.484億m3和爽島水庫1.173億m3有效庫容進行補水和稀釋，盡可能降低下游水體污染程度，控制下游信都超標2.5倍（標準0005 mg/L）以內(nèi)，封開江口鎮(zhèn)（賀江河口）超標1.5倍以內(nèi)。為此，需通過模型分析確定出滿足合面獅污水處置要求的龜石水庫和爽島水庫調(diào)水量及調(diào)水時長（圖5）。受篇幅限制，本文僅探討上游龜石水庫按180 m3/s的固定流量進行放水時下游爽島水庫的最優(yōu)調(diào)水方案。

由于合面獅庫區(qū)為狹長型河道，水流和水質(zhì)計算都可采用一維模型。假定應(yīng)急處置前合面獅庫區(qū)污染物濃度已混合均勻。根據(jù)圖5中所示信息，計算出調(diào)水處置10 d（1旬）后合面獅庫區(qū)范圍內(nèi)濃度變化見圖6。

從圖6可看出，實施應(yīng)急處置10 d后水庫壩前（圖6中區(qū)域右邊界最上端）濃度已稀釋到0012 mg/L，超標約14倍，庫區(qū)范圍有一半以上河段濃度降到了0005 mg/L。可見上游水庫的補水有力稀釋了庫區(qū)污染水體，降低了出庫污水的濃度。

根據(jù)得到的出庫水流濃度過程，分別計算爽島水庫調(diào)水100 m3/s、150 m3/s、180 m3/s 以及250 m3/s、300 m3/s時合面獅下游信都和江口鎮(zhèn)10 d后對應(yīng)濃度過程見圖7和圖8。其中，由于東安江匯入斷面在信都下游，因此信都斷面濃度過程不受爽島水庫調(diào)水量的影響。

從圖7和圖8中可看出，上游龜石水庫按180 m3/s的流量放水75 d后可滿足信都超標25倍以內(nèi)的控制要求。而爽島水庫若按同樣流量進行調(diào)水，則大約需10 d的時間才能滿足江口超標1.5倍以內(nèi)的控制要求，但在調(diào)水流量為250 m3/s時可將調(diào)水時間縮短到75 d。因此，初步確定調(diào)水方案為龜石水庫放水180 m3/s，爽島水庫對應(yīng)調(diào)水250 m3/s，應(yīng)急處置75 d后能滿足下游信都和江口斷面水質(zhì)控制要求。

為進一步確定各水庫間調(diào)水關(guān)系，盡可能節(jié)約水資源，對爽島水庫不同調(diào)水流量下允許的調(diào)水時長和調(diào)水10 d后的江口水質(zhì)進行分析，并計算出對應(yīng)調(diào)水方案下所需的處置時間和耗水量，結(jié)果見表1

從表1中可看出，在爽島調(diào)水流量為300 m3/s時，應(yīng)急處置所需時間和有效庫容允許的調(diào)水時間最為接近，但此時允許調(diào)水時間只有45 d，而根據(jù)圖7，在龜石水庫補水180 m3/s情況下信都斷面至少要75 d才能滿足水質(zhì)控制要求，加之流量過大有可能帶來洪水問題，300 m3/s的調(diào)水方案并不可行。相比而言，在爽島調(diào)水流量為180 m3/s時，雖然應(yīng)急處置所需時間比有效庫容允許的調(diào)水時間多23 d，但分析圖8可發(fā)現(xiàn)，在合面獅污水開始下排的64 h以內(nèi)江口濃度并沒有超控制線，而這64 h的時長正好與調(diào)水180 m3/s方案超出的允許調(diào)水時長接近，這說明如果爽島水庫在起調(diào)時間上適當滯后的話，采用180 m3/s調(diào)水方案是可行的。事實上可給出一種可行方案：控制爽島水庫起調(diào)后的出庫水流剛好滯后合面獅開始排污后25 d流入賀江干流，按180 m3/s的方案調(diào)度爽島水庫有效庫容用完后可結(jié)束合面獅水庫應(yīng)急處置任務(wù)。

另外，根據(jù)水利廳下發(fā)的《賀江重金屬污染處置應(yīng)急調(diào)水方案》，事發(fā)后擬定的調(diào)水方案之一是：合面獅按300 m3/s下泄污水，龜石按180 m3/s放水，而爽島水庫調(diào)則按181 m3/s調(diào)水。根據(jù)水利廳《信息簡報》第8期，截止7月14日18時賀江流域鎘指標全線達標。從資料中可看出，本文研究的調(diào)度方案和處置效果都與實際情況較為吻合。

4結(jié)論

本文從我國日益嚴峻的河流水污染問題出發(fā)，研究了基于水庫調(diào)度的河流突發(fā)水污染應(yīng)急處置問題。論文首先構(gòu)建了污染物溯源、預測和水庫調(diào)度模型，建立了一套較為完善的水庫應(yīng)急調(diào)度技術(shù)體系。將技術(shù)運用于2013年賀江水污染事件中合面獅水庫的污水應(yīng)急處置，經(jīng)分析得出了在“合面獅水庫按300 m3/s流量下泄污水、龜石水庫按180 m3/s下放清水” 情況下的爽島水庫調(diào)水方案：控制爽島起調(diào)后的出庫水流剛好滯后合面獅開始排污后25 d流入賀江干流，按180 m3/s調(diào)度爽島水庫1173億m3有效庫容用完。所得出的調(diào)度方案和處置效果與實際情況較為吻合。通過實例應(yīng)用，一方面論證了天然河流運用水庫調(diào)度進行突發(fā)水污染處置的可行性，另一方面也說明本文研究的技術(shù)體系可用于指導實際突發(fā)水污染事件的應(yīng)急處置。

參考文獻（References）：[HJ1.9mm]

[1]徐興東.流域突發(fā)性水污染事故風險應(yīng)急防范系統(tǒng)研究[D].蘭州：蘭州大學，2008.（XU X D.Study on the risk prevention systems for basin sudden water pollution accidents[D].Lanzhou：Lanzhou University，2008.（in Chinese））

[2]李青云，趙良元，林莉，等.突發(fā)性水污染事故應(yīng)急處理技術(shù)研究進展[J].長江科學院院報，2014，31（4）：611.（LI Q Y，ZHAO L Y，LIN L，et al.Review of emergency treatment technologies for sudden water pollution accidents[J].Journal of Yangtze River Scientific Research Institute，2014，31（4）：611.（in Chinese）） DOI：10.3969/j.issn.10015485.2014.04.002.

[3]陸曦，梅凱.突發(fā)性水污染事故的應(yīng)急處理[J].中國給水排水，2007，23（8）：1418.（LU X，MEI K.Emergency treatment of sudden water pollution accident[J].China Water&Wastewater;，2007，23（8）：1418.（in Chinese）） DOI：10.3969/j.issn.10015485.2014.04.002.

[4]徐月華.南水北調(diào)東線一期工程南四湖突發(fā)水污染仿真模擬及應(yīng)急處置研究[D].濟南：山東大學，2014.（XU Y H.Simulation and emergency disposal research on sudden water pollution of Nansi Lake for the firststage of the Eastern Route of the SouthtoNorth water Diversion Project[D].Jinan：Shandong University，2014.（in Chinese））

[5]蘇友華.崇左市突發(fā)性水污染事件應(yīng)急調(diào)水分析[J].企業(yè)科技與發(fā)展月刊，2011，20：115117.（SU Y H.The Analysis of emergent water dispatch for sudden water contamination events in Chongzuo[J].Enterprise Science And Technology & Development，2011，20：115117.（in Chinese）） DOI：10.3969/j.issn.16740688.2011.20.042.

[6]辛小康，葉閩，尹煒.長江宜昌江段水污染事故的水庫調(diào)度措施研究[J].水電能源科學，2011，29（6）：4648.（XIN X K，YE M，YIN W.Research on reservoir operation measure of water pollution accident in Yichang Segment of Yangtze River[J].Water Resources and Power，2011，29（6）：4648.（in Chinese）） DOI：10.3969/j.issn.10007709.2011.06.015.

[7]陶亞，任華堂，夏建新.突發(fā)水污染事故不同應(yīng)對措施處置效果模擬[J].應(yīng)用基礎(chǔ)與工程科學學報，2013，21（2）：203213.（TAO Y，REN H T，XIA J X.Effect analysis of different emergency measures for accidental water pollution[J].2013，21（2）：203213.（in Chinese）） DOI：10.3969/j.issn.10050930.2013.02.002.

[8]余真真，張建軍，馬秀梅，等.小浪底水庫應(yīng)急調(diào)度對下游水污染事件的調(diào)控[J].人民黃河，2014，36（8）：7375，100.（YU Z Z，ZHANG J J，MA X M，et al.Control ranges of emergency dispatches in Xiaolangdi reservoir for downstream water pollution incidents[J]..Yellow River，2014，36（8）：7375，100.（in Chinese）） DOI：10.3969/j.issn.10001379.2014.08.022.

[9]郝麗娟，綦中躍.應(yīng)對突發(fā)性水污染事故存在的問題和建議[J].北京水務(wù)，2007（2）：3435.（HAO L J，QI Z Y.Question and advice of facing accident water pollution[J].Beijing Water，2007（2）：3435.（in Chinese）） DOI：10.3969/j.issn.16734637.2007.02.016.[ZK）]

[10][ZK（#]王船海，李光熾.實用河網(wǎng)水流計算[M].南京：河海大學出版社，2003，9.（WANG C H，LI G C.Practical flow calculation of river network[M].Nanjing：Hohai University Press，2003，9.（in Chinese））