近岸水域環(huán)境容量計(jì)算及實(shí)例應(yīng)用分析研究

摘要：提出了適用于近岸水域的環(huán)境容量計(jì)算方法，通過構(gòu)建三維水動(dòng)力、三維生態(tài)動(dòng)力模型對(duì)近岸水域的水動(dòng)力、水質(zhì)過程進(jìn)行了模擬驗(yàn)證，結(jié)合水質(zhì)功能分區(qū)，將容量計(jì)算歸結(jié)為線性規(guī)劃問題，求解在規(guī)劃的水污染控制單元?jiǎng)澐趾团盼劭谖恢脳l件下的水環(huán)境容量。以大亞灣水域?yàn)閷?shí)例求解與分析了CODMn環(huán)境容量，驗(yàn)證了計(jì)算方法的適用性和有效性。

關(guān)鍵詞：近岸水域；環(huán)境容量；水動(dòng)力；生態(tài)動(dòng)力；線性規(guī)劃

中圖分類號(hào)： X26

文獻(xiàn)標(biāo)識(shí)碼： A 文章編號(hào)： 16749944（2015）06025104

1 引言

環(huán)境容量概念最早由日本環(huán)境學(xué)界在20世紀(jì)60年代提出[1]，國內(nèi)外的研究者對(duì)水環(huán)境容量提出過多種表述的定義[2～4]，具體可描述為：在一定的排污口位置下，對(duì)象水域在水環(huán)境質(zhì)量達(dá)到環(huán)境標(biāo)準(zhǔn)的前提下所能容納污染物的最大量。水環(huán)境容量的影響因素除了排污口的位置和環(huán)境標(biāo)準(zhǔn)之外，還有支配污染物輸移、稀釋和擴(kuò)散的水動(dòng)力條件以及污染物本身的降解速率等。按照計(jì)算條件的不同，環(huán)境容量可區(qū)分為兩類，一是天然環(huán)境容量，就是在假設(shè)各類污染物背景濃度為0時(shí)（即假設(shè)水體完全不受污染），所求出的各排污口總排放最大負(fù)荷量；二是背景環(huán)境容量，是指在只考慮外海邊界污染物輸入時(shí)，所求出的各排污口總排放最大負(fù)荷量。水環(huán)境容量評(píng)估工作對(duì)近岸區(qū)域的可持續(xù)發(fā)展具有重要意義，本研究具體提出了適用于近岸水域環(huán)境容量的評(píng)估計(jì)算方法，并結(jié)合實(shí)例應(yīng)用開展研究分析。

2 計(jì)算方法

根據(jù)研究水域的地形、水動(dòng)力特征，建立起反映研究水域?qū)嶋H動(dòng)力情況的三維水動(dòng)力，并根據(jù)觀測數(shù)據(jù)對(duì)建立起來的水動(dòng)力過程進(jìn)行模擬，然后根據(jù)研究水域的水質(zhì)及生態(tài)現(xiàn)狀，建立起反映水質(zhì)情況的三維生態(tài)動(dòng)力模型，并根據(jù)觀測數(shù)據(jù)對(duì)建立起來的生態(tài)動(dòng)力模型進(jìn)行水質(zhì)過程模擬，并驗(yàn)證模型應(yīng)用于近岸水域的可靠性，最后根據(jù)研究水域現(xiàn)有的水質(zhì)功能分區(qū)，求解在規(guī)劃的水污染控制單元?jiǎng)澐趾团盼劭谖恢脳l件下的水環(huán)境容量。

2.1 水動(dòng)力與生態(tài)動(dòng)力過程模擬

采用ECOMSED（Estuaries and Coastal Ocean Model with Sediment Module）模型[5]進(jìn)行水動(dòng)力過程模擬。ECOMSED為三維數(shù)學(xué)模型，模型主要模擬變量包括三維流速、溫度、鹽度、湍動(dòng)能、湍流長度等。模型的動(dòng)力方程為非線性的，考慮由于溫、鹽導(dǎo)致流體密度差異所產(chǎn)生的上升和下降流，以及流體的水平對(duì)流。在進(jìn)行水動(dòng)力模型的構(gòu)建和驗(yàn)證時(shí)，首先根據(jù)近岸水域的地形、水動(dòng)力特征，建立起反映近岸水域?qū)嶋H動(dòng)力情況的三維水動(dòng)力；接著則根據(jù)現(xiàn)場觀測數(shù)據(jù)對(duì)建立起來的水動(dòng)力過程進(jìn)行模擬，并對(duì)結(jié)果進(jìn)行驗(yàn)證，以檢驗(yàn)?zāi)Ｐ偷目煽啃浴ＤＰ偷睦字Z動(dòng)力方程為：

式中：（U，V）為水平流速矢量，為水平梯度運(yùn)算符，ρ0為海水密度常數(shù)，ρ是實(shí)際密度，g是重力加速度，P是壓力，KM是垂向湍動(dòng)能混合的渦粘系數(shù)，f是隨緯度變化的科氏力參數(shù)。

采用RCA生態(tài)動(dòng)力模型（Row-Column Aesop）[6]進(jìn)行水質(zhì)過程模擬。RCA模型是從Aesop箱式生態(tài)模型及美國環(huán)保署的WASP水質(zhì)模型發(fā)展而來，RCA包括溶解氧、碳循環(huán)（其中包括了浮游植物的生長）等五個(gè)相互作用系統(tǒng)。RCA最明顯的特征是耦合了一個(gè)底泥營養(yǎng)鹽通量模型，可以根據(jù)每天沉積到底泥中的顆粒態(tài)有機(jī)碳等形態(tài)的污染物通量及底泥中原有的顆粒態(tài)有機(jī)碳等形態(tài)污染物的含量，動(dòng)態(tài)地模擬研究區(qū)域污染物底通量及底泥耗氧量的時(shí)空變化過程。RCA水質(zhì)模型描述各水質(zhì)因子物理、生化過程的方程的一般形式可用下式表示：

式（4）中：c為水質(zhì)因子濃度（mg/L），t為時(shí)間（s），E為由于潮汐、密度或速度梯度所引起的梯度導(dǎo)致的擴(kuò)散系數(shù)（m2/s），U為對(duì)流速度（m/s），S為各水質(zhì)因子內(nèi)部源匯項(xiàng)（mg/（L·s）），W為外部源匯項(xiàng)（mg/（L·s））。

用于構(gòu)建水動(dòng)力及生態(tài)動(dòng)力模型的數(shù)據(jù)主要分為四類，分別為模型驗(yàn)證數(shù)據(jù)、邊界條件數(shù)據(jù)、環(huán)境容量水質(zhì)目標(biāo)數(shù)據(jù)及地形數(shù)據(jù)。模型驗(yàn)證數(shù)據(jù)主要是用于水動(dòng)力模型及水質(zhì)模型中各參數(shù)率定及計(jì)算結(jié)果驗(yàn)證的各類水動(dòng)力及水質(zhì)數(shù)據(jù)；邊界條件數(shù)據(jù)是用于驅(qū)動(dòng)水動(dòng)力模型及生態(tài)動(dòng)力模型所需的在外海邊界、河流邊界及海表邊界所需要的水位、流量、污染物濃度、氣象條件等數(shù)據(jù)；環(huán)境容量水質(zhì)目標(biāo)依據(jù)是研究水域執(zhí)行的水環(huán)境質(zhì)量標(biāo)準(zhǔn)的有關(guān)規(guī)定，近岸水域水質(zhì)目標(biāo)參考《中華人民共和國海水水質(zhì)標(biāo)準(zhǔn)》（GB3097-1997）；地形數(shù)據(jù)是研究水域的水底地形數(shù)據(jù)。

2.2 容量求解

為了使環(huán)境容量的求解更加科學(xué)化與合理化，同時(shí)為了避免傳統(tǒng)的試算法解環(huán)境容量所帶來的巨大工作量，環(huán)境容量的計(jì)算問題可轉(zhuǎn)化為線性規(guī)劃問題。

根據(jù)定義，水環(huán)境容量問題可表述為：在選定的一組水質(zhì)控制斷面的指標(biāo)污染物濃度不超過其各自對(duì)應(yīng)的環(huán)境標(biāo)準(zhǔn)值的前提下，使各排污口的污染負(fù)荷排放量之和最大，即

式中：L為對(duì)象水域所有排污口的總排放負(fù)荷量；x為某個(gè)排污口的排放負(fù)荷量；i為水質(zhì)控制斷面編號(hào)；m為水質(zhì)控制斷面數(shù)目；j為排污口編號(hào)；n為排污口數(shù)目；aij為第j個(gè)排污口的單位負(fù)荷量對(duì)第i個(gè)水質(zhì)控制斷面的污染貢獻(xiàn)度系數(shù)；Cbi為水質(zhì)控制斷面處的污染背景濃度；為水質(zhì)控制斷面處的環(huán)境標(biāo)準(zhǔn)值。

式（6）左邊第一項(xiàng)表示個(gè)排污口對(duì)控制點(diǎn)的水質(zhì)濃度貢獻(xiàn)量的總和，亦即區(qū)域內(nèi)排污引起的控制點(diǎn)處的濃度增值；第二項(xiàng)是該控制點(diǎn)的背景濃度，這是在沒有污染源的條件下，僅僅由邊界水質(zhì)影響而產(chǎn)生的濃度。因此式（6）左邊兩項(xiàng)之和實(shí)際上是控制點(diǎn)處的濃度。由于水質(zhì)擴(kuò)散方程是線性的，濃度有可迭加性，所以用線性迭加的方法來求解某一點(diǎn)的濃度是可行的。因此，求解海域環(huán)境容量問題可歸結(jié)到求解式（5）～式（7）所表達(dá)的線性規(guī)劃問題。線性規(guī)劃問題用單純型法（Simplex Method）求解。

對(duì)于排污口的設(shè)置，主要按照三個(gè)方面進(jìn)行統(tǒng)籌劃分，這三個(gè)方面分別是沿岸排污口調(diào)查結(jié)果、集水區(qū)域分布（面源集污單元）及近岸水域功能區(qū)劃。沿岸排污口調(diào)查主要是指通過對(duì)沿岸工廠污水、城鎮(zhèn)居民生活污水、水產(chǎn)養(yǎng)殖及入海河流等直接進(jìn)入近岸水域的污染源進(jìn)行野外調(diào)查。排污口設(shè)置還應(yīng)考慮近岸水域功能區(qū)劃，因?yàn)楸狙芯凯h(huán)境容量的計(jì)算，不僅要求出現(xiàn)存排污口的環(huán)境容量，還要求出各類未來將設(shè)置的排污口的環(huán)境容量，這樣才能為未來的排污決策提供科學(xué)依據(jù)。綜上所述，近岸水域環(huán)境容量計(jì)算排污口設(shè)置的原則總結(jié)為以下4點(diǎn)。

（1）現(xiàn)存沿岸排污口的位置應(yīng)設(shè)置排污口；

（2）每一個(gè)集水區(qū)域至少要設(shè)置一個(gè)排污口，集水區(qū)域的設(shè)置依據(jù)各集水區(qū)域河流入海口所處位置而定；

（3）每一個(gè)水質(zhì)功能區(qū)至少要設(shè)有一個(gè)排污口；

（4）近岸水域功能區(qū)劃規(guī)劃的排污口應(yīng)考慮設(shè)置排污口；

對(duì)于水質(zhì)控制點(diǎn)的選定，根據(jù)《海水水質(zhì)標(biāo)準(zhǔn)》（GB3097-1997）[7]，海水水質(zhì)分為4類，為了控制各功能區(qū)的水質(zhì)達(dá)到其對(duì)應(yīng)的環(huán)境標(biāo)準(zhǔn)（圖1），對(duì)水質(zhì)控制點(diǎn)的選取須滿足以下要求。

（1）對(duì)處在采用海水水質(zhì)第三類標(biāo)準(zhǔn)功能區(qū)的排污口，在其周邊設(shè)第三類水質(zhì)標(biāo)準(zhǔn)控制點(diǎn)；

（2）對(duì)處在采用海水水質(zhì)第二類標(biāo)準(zhǔn)功能區(qū)的排污口，在其周邊設(shè)第二類水質(zhì)標(biāo)準(zhǔn)控制點(diǎn)；

（3）對(duì)處在采用海水水質(zhì)第一類標(biāo)準(zhǔn)功能區(qū)的排污口，在其周邊設(shè)第一類水質(zhì)標(biāo)準(zhǔn)控制點(diǎn)；

（4）在采用第三類與第二類海水水質(zhì)標(biāo)準(zhǔn)的2種水域的交界線上設(shè)第二類水質(zhì)控制點(diǎn)；

（5）在采用第一類與第二類海水水質(zhì)標(biāo)準(zhǔn)的2種水域的交界線上設(shè)第一類水質(zhì)控制點(diǎn)，使這些點(diǎn)的濃度不得超過第一類海水水質(zhì)標(biāo)準(zhǔn)。

對(duì)于貢獻(xiàn)度系數(shù)的計(jì)算，由于海洋受潮汐影響，一個(gè)連續(xù)恒定的污染源對(duì)某一控制點(diǎn)的影響是隨時(shí)間變化的，所以貢獻(xiàn)度系數(shù)也是時(shí)間的函數(shù)。如果將貢獻(xiàn)度系數(shù)的時(shí)間變化也加以考慮，就會(huì)更加實(shí)際，但同時(shí)也使得容量計(jì)算量劇增。而且如果給出一個(gè)潮周期內(nèi)隨時(shí)間變化的水環(huán)境容量，在現(xiàn)階段根據(jù)環(huán)境容量變化來實(shí)施管理操作也具有難度。因此貢獻(xiàn)度系數(shù)以及環(huán)境容量都取一個(gè)月內(nèi)的平均值。

3 計(jì)算實(shí)例應(yīng)用分析

3.1 水環(huán)境容量計(jì)算因子的確定

容量計(jì)算因子的選擇考慮以下幾方面的原則：

（1）選擇海洋管理部門監(jiān)測管理和制定海洋污染治理措施一般考慮的控制因子；

（2）選擇符合研究海域環(huán)境污染特征，能夠反映研究海域污染水平，決定水環(huán)境質(zhì)量的主要超標(biāo)因子或定類因子；

（3）選擇能夠反映水域污染物排放現(xiàn)狀和趨勢，決定水污染特征的主要污染負(fù)荷因子；

（4）選擇各行業(yè)污染源主要監(jiān)控的污染因子；

（5）選擇被收入在《海水水質(zhì)標(biāo)準(zhǔn)》（GB3097-1997）和《污水綜合排放標(biāo)準(zhǔn)》（GB8978-1996）之內(nèi)的因子；

（6）選擇具有可靠的監(jiān)測計(jì)量手段和在線監(jiān)控方法的控制因子；

（7）選擇具有切實(shí)可行的污染控制和削減措施的控制因子。

根據(jù)上述原則要求，重點(diǎn)考慮海洋管理部門需要、研究水域?qū)嶋H情況和掌握資料的充分性三方面因素，本研究選擇CODMn作為計(jì)算因子。

3.2 計(jì)算實(shí)例與結(jié)果分析

根據(jù)掌握的數(shù)據(jù)資料，以大亞灣水域作為應(yīng)用實(shí)例進(jìn)行分析。大亞灣是廣東省沿岸最大的海灣之一，位于珠江口東側(cè)（114°29′42″～114°49′42″E，22°31′12″～22°50′00″N），正在大力發(fā)展海洋經(jīng)濟(jì)，石化等產(chǎn)業(yè)已成為區(qū)域經(jīng)濟(jì)發(fā)展支柱，而大亞灣又是國家級(jí)海龜自然保護(hù)區(qū)和省級(jí)水產(chǎn)資源自然保護(hù)區(qū)所在地，經(jīng)濟(jì)發(fā)展與環(huán)境保護(hù)在并行中博弈。水環(huán)境容量對(duì)于大亞灣區(qū)域的發(fā)展至關(guān)重要，而摸清大亞灣水域的環(huán)境容量是大亞灣區(qū)域?qū)嵤┛沙掷m(xù)發(fā)展戰(zhàn)略的前提條件。

根據(jù)《大亞灣海洋環(huán)境質(zhì)量現(xiàn)狀調(diào)查報(bào)告》關(guān)于污染源和水環(huán)境質(zhì)量現(xiàn)狀調(diào)查的基礎(chǔ)數(shù)據(jù)[8]，結(jié)合大亞灣海域總體規(guī)劃及社會(huì)經(jīng)濟(jì)發(fā)展與環(huán)境保護(hù)的需要，以CODMn為指標(biāo)，在大亞灣水體功能區(qū)劃的基礎(chǔ)上，運(yùn)用經(jīng)過驗(yàn)證的水動(dòng)力模型和生態(tài)動(dòng)力模型求出大亞灣沿岸深圳龍崗區(qū)、惠州大亞灣區(qū)及惠東縣在規(guī)劃的水污染控制單元?jiǎng)澐趾团盼劭谖恢脳l件下的水環(huán)境容量。用單純性法求解大亞灣水域環(huán)境容量解線性規(guī)劃問題式（1）～式（3），得到總排放負(fù)荷量最大時(shí)的各排污口允許排放CODMn環(huán)境容量，結(jié)果見表1（排污口對(duì)應(yīng)位置見圖2）。計(jì)算結(jié)果包括兩種類型的容量，天然容量和背景容量。

從計(jì)算結(jié)果可以看出，惠州市惠東縣占了大亞灣沿岸三個(gè)主要行政區(qū)天然環(huán)境容量的49%，深圳龍崗區(qū)和惠州市大亞灣區(qū)的CODMn天然環(huán)境容量分別占大亞灣各污染物天然環(huán)境容量的19%和32%。背景環(huán)境容量與天然環(huán)境容量有著類似的空間分布特點(diǎn)。環(huán)境容量的計(jì)算結(jié)果與排污口位置、控制點(diǎn)水質(zhì)標(biāo)準(zhǔn)及水體交換能力有關(guān)。由于排污口位置及控制點(diǎn)水質(zhì)標(biāo)準(zhǔn)是確定的，所以使環(huán)境容量存在區(qū)域差異的主導(dǎo)因素應(yīng)該是水體交換能力。

由于惠東的海岸線占據(jù)了整個(gè)大亞灣東岸，而大亞灣東岸是大亞灣海域水體交換能力最強(qiáng)的區(qū)域，因此，惠東占有將近50%的天然環(huán)境容量。大亞灣區(qū)及龍崗區(qū)由于所占有的區(qū)域是大亞灣海域中水體交換能力最差的西岸，導(dǎo)致其天然環(huán)境容量相對(duì)較小。大亞灣沿岸各排污口的天然和背景環(huán)境容量空間分布基本與水體交換能力強(qiáng)弱相關(guān)。另一方面，結(jié)合大亞灣水域的監(jiān)測數(shù)據(jù)分析，大亞灣北部范和港及大鵬澳等近岸水域多在三類水或差于三類水，而這些區(qū)域的天然環(huán)境容量和背景環(huán)境容量計(jì)算結(jié)果都較少，表明本研究對(duì)大亞灣環(huán)境容量的計(jì)算是合理的。

4 結(jié)語

根據(jù)近岸水域的地形、水動(dòng)力特征，利用ECOMSED建立起反映近岸水域?qū)嶋H動(dòng)力情況的三維水動(dòng)力模型并根據(jù)現(xiàn)場觀測數(shù)據(jù)對(duì)建立起來的水動(dòng)力過程進(jìn)行模擬，采用RCA三維生態(tài)動(dòng)力模型對(duì)近岸水域的水質(zhì)過程進(jìn)行模擬，并驗(yàn)證模型應(yīng)用于近岸水域的可靠性，結(jié)合水域水質(zhì)功能分區(qū)，將容量計(jì)算歸結(jié)為線性規(guī)劃問題，從而求解在規(guī)劃的水污染控制單元?jiǎng)澐趾团盼劭谖恢脳l件下的水環(huán)境容量。以大亞灣水域?yàn)閷?shí)例研究對(duì)象，對(duì)大亞灣水域的CODMn環(huán)境容量進(jìn)行了求解與分析，計(jì)算結(jié)果與實(shí)際狀況相符合理，表明計(jì)算方法具有較強(qiáng)的適用性和有效性，為污染物的總量控制提供科學(xué)依據(jù)，也為近岸水域的環(huán)境管理提供決策依據(jù)和技術(shù)指導(dǎo)。

參考文獻(xiàn)：

[1]張永良， 劉培哲. 水環(huán)境容量綜合手冊[M].北京：清華大學(xué)出版社，1991：574～575.

[2]張永良，洪繼華，夏 青.我國水環(huán)境容量研究與展望[J].環(huán)境科學(xué)研究，1988，1（1）：73～81.

[3]Rauch W.， Henze M.， Koncsos L. et al. River quality modeling： I. State of the art[J].Water Science and Technology，1998，38（11）： 237～244.

[4]Litak B.G. Environmental Engineers Handbook， Vol.1， Water Pollution[R].ChiltonBook Company， 1997（12）：34～37.

[5]Blumberg A.F. A Primer for ECOMSED Version 1.3 Users Manual[R].Mahwah， N.J. USA： HydroQual， Inc.，2002.

[6]Fitzpatrick J.J. USER'S GUIDE FOR RCA（Release 3.0）[R].Mahwah， N.J. USA： HydroQual， Inc.，2004.

[7]國家環(huán)境保護(hù)局.海水水質(zhì)標(biāo)準(zhǔn)GB3097-1997[S].北京：國家環(huán)境保護(hù)局，1997.

[8]李耀初，何智遠(yuǎn).大亞灣海洋環(huán)境質(zhì)量現(xiàn)狀調(diào)查報(bào)告[R].中山大學(xué)，廣東省海洋與漁業(yè)局等.