大連灣海域排污狀況分析及防治

趙 輝

（大連市生態(tài)環(huán)境事務服務中心，遼寧 大連 116023）

1 研究背景

海洋占地球表面積的71%，是生態(tài)系統(tǒng)的重要組成成分，對人類的發(fā)展至關(guān)重要。然而海洋資源開發(fā)程度高的同時，排污問題也會隨之而來。隨著我國沿海城市的經(jīng)濟發(fā)展和城鎮(zhèn)化進程的推進，濱海城市工業(yè)廢水和生活污水所攜帶的污染物相應增加，近岸海域環(huán)境質(zhì)量污染呈加劇趨勢[1]。排污是否超量，如何防治海域污染成為許多臨海城市尋求可持續(xù)發(fā)展的重要課題。大連灣位于遼東半島最南端，其岸線長125 km[2]，是北方最大的水產(chǎn)品交易中心，海域資源高度開發(fā)，有較高的生態(tài)敏感性；大連市有近600萬人口，排污量大；大連灣為半封閉性海域，海水交換周期長，水體內(nèi)無機氮，無機磷等含量較高[3]，這些都給大連灣的污染防治工作帶來極大的考驗。該文對大連灣附近的排污狀況進行了量化分析，并通過數(shù)學模型等手段分析污染物的演變規(guī)律，這些模型可以合理推斷污染物對相關(guān)海域的影響，以及對相應防治策略的有效性進行預測。

2 研究對象與研究思路

2.1 研究范圍

該文主要研究了大連灣海域的排污及防治措施。大連灣海域隸屬于大連市甘井子區(qū)，根據(jù)《海灣志》對大連灣海域的劃分，結(jié)合環(huán)境管理的相關(guān)標準，我們研究的區(qū)域范圍是東經(jīng) 121°36'—121°52'，北緯 38°55'—39°4'，總研究海域面積為161 km2。

2.2 研究方法

該文認為，科學防治大連灣海域污染大致分為3個步驟。1）基于數(shù)據(jù)分析大連灣排污現(xiàn)狀，根據(jù)污染狀況的比例，找出重要的污染源。2）建立數(shù)學模型，預測大連灣水質(zhì)的演化規(guī)律，對無干預情況下大連灣水域的發(fā)展演化做出預測。3）提出科學防治的措施，結(jié)合數(shù)學模型分析為什么相應的措施可以有效解決污染狀況。該文的研究思路框圖如圖1所示。

圖1量化分析研究策略

3 大連灣污染物數(shù)據(jù)分析

根據(jù)2016年監(jiān)測統(tǒng)計數(shù)據(jù)，大連灣廢水排量為50 325萬噸，其中18 892萬噸來自工業(yè)源排放，25 510萬噸來自生活源，剩余5 923萬噸來自降雨徑流排放。各種廢水排放的源頭比例情況如圖2所示。

圖2 大連灣陸源廢水排放量分布

大連灣總氮排放量為9 293 t，其中627 t為工業(yè)源排放，7 352 t為生活源排放，1 314 t位降雨徑流排放，各源頭總氮排放量如圖3所示。

我們也對其他幾種污染物進行了數(shù)據(jù)分析，大連灣化學需氧量（COD）排放為53 951 t，其中4 264 t為工業(yè)源排放，37 466 t為生活源排放，12 321噸則來自于降雨徑流排放;氨氮排放總量為4 702 t，其中111 t來自工業(yè)源排放，3 814 t來自生活源排放，777 t為降雨徑流排放;總磷排放量為613 t，其中27 t來自于工業(yè)排放，489 t來自于生活排放，而97 t來自于降雨徑流排放。

圖3 大連灣陸源總氮排放量分布

從大連灣各種污染物排放源頭的數(shù)據(jù)來看，生活源是大部分污染物的主要來源，其次是降雨徑流，然后是工業(yè)源，但是對于廢水排放來說，工業(yè)源排放遠大于降雨徑流排放。從污染物來源的數(shù)據(jù)分析可以看出，生活源是污染防治的重點，然而對于廢水污染的防治，工業(yè)源的污染減少也十分重要。

4 污染物傳播模型與數(shù)值模擬

在上一節(jié)中，我們統(tǒng)計數(shù)據(jù)分析了水中污染物的來源，可以幫助我們從源頭處防治，減少污染物排放。然而污染程度不僅僅取決于排放狀況，相關(guān)環(huán)保指標依賴于排放污染物在海域內(nèi)濃度，因此了解污染物在水中擴散規(guī)律至關(guān)重要。該文提出了基于數(shù)學模型來模擬各種污染物在大連灣內(nèi)擴散狀況，根據(jù)數(shù)學模型，我們可以對海域內(nèi)污染物的濃度及擴散趨勢進行數(shù)據(jù)模擬，掌握科學防治的重點，也可以嘗試采取不同措施對濃度的變化進行預測，從而判斷措施的有效性。

污染物擴散模型基于三維流場模型，在近岸海域內(nèi)，污染物會隨著海水輸運擴散，在運輸過程中自然衰減，同時源頭處又不斷有新的污染物輸入，使得整個過程變得非常復雜。物理上，我們用偏微分方程來描述以上污染物擴散規(guī)律。根據(jù)水質(zhì)模型，在近岸處污染物的排放滿足以下方程：

式中：x，y為二維平面坐標，t為時間，S為任意時刻污染物的濃度，S0為污染物排放濃度，ζ表示海面水體上升的高度，H為靜水面下水深，D、w在以三角形節(jié)點為中心的控制體上計算，而u和v則基于三角形單位控制體。Kh和Ah為波流聯(lián)合作用下的垂向紊動系數(shù)和水平湍流擴散系數(shù)，可以根據(jù)疊加法得到。

結(jié)合邊界條件，我們可以對偏微分方程離散化，使用有限差分或者有限元方法來得到迭代格式，然后通過Matlab進行數(shù)值模擬。數(shù)值模擬的結(jié)果是，我們可以得到任意時刻，污染物濃度在海域內(nèi)不同位置的數(shù)值，了解不同時間段海域內(nèi)污染物的濃度，判斷污染程度S。基于水動力模擬模型，我們以無機氮濃度為例進行數(shù)值模擬，預測不同時間段大連灣海域附近的無機氮濃度，結(jié)果如圖4所示。

經(jīng)過對比，數(shù)值模擬的結(jié)果和真實監(jiān)測的結(jié)果類似。我們可用同樣的方法，對磷、氨氮等污染物進行類似的預測，根據(jù)邊界條件、污染物排入、海域特點等預測海灣各處污染物濃度。

5 大連灣海域排污防治建議

在前面兩節(jié)，我們提出數(shù)據(jù)分析的方法和數(shù)學模型方法來分析污染物的源頭以及污染物在大連灣海域內(nèi)的擴散規(guī)律，可以根據(jù)這些量化手段來合理防治。

5.1 污染物入海排污源頭總量控制

首先可以對污染入海的源頭進行控制，針對每種污染物，根據(jù)污染物來源統(tǒng)計，對每種源頭的總量進行確定，例如對于無機氮控制，應主要考察并限制生活來源的氮排放，而對于廢水污染，同時考察和限制生活來源廢水排放和工業(yè)源廢水。首先根據(jù)大連地區(qū)排污口分布，對大連各地區(qū)排污總量加以限制?？蓪⒄{(diào)整口的排污量代入上一節(jié)的擴散模型中，檢驗海域內(nèi)污染物是否超標。

圖4 大連夏季雨后灣海域無機氮濃度分布數(shù)值結(jié)果

5.2 排污分配方法

在大連灣的排污總量合理的前提下，可以綜合考慮經(jīng)濟效益和社會效益，來分配入海污染物允許排放量和消減量。依據(jù)“陸海聯(lián)動、海陸統(tǒng)籌”的原則、重視歷史和現(xiàn)狀的原則、兼顧公平和效益的原則等，從目前廣泛應用的等比例分配法、動態(tài)層次分析法和基于容量總量分配法中選取等比例分配法作為分配方法。確定分配方法后，可以用上一節(jié)的方程模型來預測相應分配策略下海域各處污染物濃度是否達標。

我們建議可以指定討論減排方案，這些方案會影響模型中的邊界條件和污染輸入項，從而影響模擬的結(jié)果，對相應措施的效果進行分析。

6 結(jié)論

大連灣是大連市工業(yè)廢水和生活污水的主要受納海域[4]，其環(huán)境治理極其重要。量化分析是多個學科的基石，通過量化分析，我們往往可以找出解決問題的關(guān)鍵，并且預測相應措施的效果。該文用數(shù)據(jù)分析來找出不同污染物的源頭，并提出用數(shù)學模型來預測防治措施的有效性，可以大大降低項目成本，高效篩選出行之有效的策略。這些研究方法不僅適用于大連灣海域，也適合推廣到其他海域。