郁斢蘭，王 諾，劉忠波，李 慧 （大連海事大學(xué)交通運(yùn)輸管理學(xué)院，遼寧 大連 116026）

自凈作用下渤海海域COD環(huán)境容量研究

郁斢蘭，王 諾*，劉忠波，李 慧 （大連海事大學(xué)交通運(yùn)輸管理學(xué)院，遼寧 大連 116026）

為測(cè)算自凈作用下我國(guó)渤海COD環(huán)境容量，按地理特征將渤海水域劃分為萊州灣、渤海灣、遼東灣以及渤海中部4個(gè)分區(qū)，從水動(dòng)力交換擴(kuò)散和化學(xué)降解兩個(gè)方面入手，一方面利用各分區(qū)水半交換周期分別求得COD的擴(kuò)散能力，另一方面利用COD衰減規(guī)律得到其降解能力，綜合后得到各分區(qū)在不同水質(zhì)標(biāo)準(zhǔn)下COD的環(huán)境容量.研究表明，渤海各灣水域COD濃度和環(huán)境容量各不相同，其中，萊州灣水質(zhì)屬于較清潔，在此區(qū)間內(nèi)還有26%的環(huán)境余量；渤海灣水質(zhì)屬于輕度污染，在此區(qū)間內(nèi)還有5%的環(huán)境余量，遼東灣水質(zhì)屬于較清潔，在此區(qū)間內(nèi)還有88%的環(huán)境余量.

海洋；環(huán)境；容量；COD；渤海

海洋對(duì)污染物本身具有巨大的自凈能力，這一能力主要來(lái)自兩個(gè)方面，一是水體會(huì)持續(xù)不斷從沿岸向遠(yuǎn)海流動(dòng)交換，使得水中污染物得以擴(kuò)散，濃度逐漸被稀釋；再是大自然具有生物化學(xué)降解作用，部分污染物隨著時(shí)間推移會(huì)逐漸自行降解.海洋的這種自凈能力，使得其使所接納的污染物濃度始終處于動(dòng)態(tài)變化中.

目前，關(guān)于海域水交換問(wèn)題已有一些研究成果，如引入半交換時(shí)間的概念，采用類似于放射性同位素的半衰期來(lái)研究海域通過(guò)對(duì)流擴(kuò)散作用稀釋為初始濃度一半時(shí)所需的時(shí)間［1］；采用不同水交換數(shù)值模擬的計(jì)算方法［2-5］；研究不同海域水交換能力和特性等［6-8］.關(guān)于COD生物化學(xué)降解問(wèn)題，邱巍等［9］通過(guò)大量實(shí)驗(yàn)研究了長(zhǎng)江口竹園排污口的COD在各種條件下的化學(xué)降解過(guò)程，得出降解系數(shù)主要跟COD濃度和水體溫度有關(guān)；郭棟鵬［10］通過(guò)實(shí)驗(yàn)發(fā)現(xiàn)COD降解規(guī)律符合一級(jí)反應(yīng)動(dòng)力學(xué)過(guò)程.在海洋環(huán)境容量方面，Zhao等［11］利用三維對(duì)流擴(kuò)散模型并結(jié)合降解過(guò)程模擬了渤海COD濃度分布，推算了渤海水域的COD環(huán)境容量.綜上分析發(fā)現(xiàn)，因不同海域在地理位置、水域面積、水交換能力以及岸線與地形復(fù)雜程度上的差異，對(duì)面積較大水域僅以平均指標(biāo)描述水交換過(guò)程及降解效率并不能有效反映各局部水域的水質(zhì)情況，整體上水質(zhì)標(biāo)準(zhǔn)清潔，其局部水域的污染可能已到了相當(dāng)嚴(yán)重的程度.因此， 對(duì)于大范圍海域的研究， 對(duì)于大范圍海域的研究，只有具體分析各局部水域的污染狀況和環(huán)境容量，才能夠有針對(duì)性的對(duì)海洋環(huán)境進(jìn)行管理.

渤海是我國(guó)水域面積最大的半封閉內(nèi)海，水域面積約75600km2，平均水深18m，東面以寬59n mile的渤海海峽與黃海相通.近幾十年來(lái)，隨著渤海沿海地區(qū)經(jīng)濟(jì)建設(shè)的飛速發(fā)展，海洋污染情況日益嚴(yán)重，降低陸源污染投資代價(jià)巨大.為了有效制定環(huán)境保護(hù)政策，避免水質(zhì)污染狀況繼續(xù)惡化，結(jié)合自凈能力分析計(jì)算渤海中各海灣在不同污染標(biāo)準(zhǔn)下COD環(huán)境容量，已成為當(dāng)務(wù)之急. 文中按地理特征將渤海水域劃分為萊州灣、渤海灣、遼東灣以及渤海中部4個(gè)分區(qū)，從水動(dòng)力交換擴(kuò)散和化學(xué)降解兩個(gè)方面入手，研究各分區(qū)在不同水質(zhì)標(biāo)準(zhǔn)下COD的環(huán)境容量.

1 計(jì)算模型

1.1 水動(dòng)力與水質(zhì)模型

海洋水交換的能力主要取決于水半交換的時(shí)間.根據(jù)水流形成機(jī)理，構(gòu)建渤海水半交換時(shí)間的數(shù)值模型由兩部分組成:即在風(fēng)、潮汐和熱鹽作用下瞬時(shí)流動(dòng)的水動(dòng)力模型和計(jì)算保守物質(zhì)濃度輸運(yùn)的水質(zhì)模型.前一模型的分析結(jié)果可為后一模型提供對(duì)流輸運(yùn)場(chǎng)和湍流擴(kuò)散系數(shù)，通過(guò)對(duì)水質(zhì)點(diǎn)的追蹤以及考慮風(fēng)、潮汐和水密度梯度的作用，便可計(jì)算出水半交換模型的特征值.

本文采用斜壓模型作為水動(dòng)力模型，該模型已多次用于渤海海域的水動(dòng)力分析［8］.設(shè)岸界水質(zhì)點(diǎn)的法向流速為零，利用渤海海峽水文站大長(zhǎng)山島水文站（39°16'N，122°35'E）和雞鳴島水文站（37°7'N，122°29'E）實(shí)測(cè)結(jié)果，采用HAMSOM三維斜壓分層模型，以5'×5'水平網(wǎng)格，垂直分10層建立水動(dòng)力模型.水質(zhì)變量的質(zhì)量守恒控制方程為:

式中:C表示水質(zhì)狀態(tài)變量濃度；u、v和w分別為x、y和z方向的速度；Ax、Ay和Az分別為x、y和z方向的擴(kuò)散系數(shù)；Sc為單位體積源匯項(xiàng).

由于式（1）包含了物理輸運(yùn)、平流擴(kuò)散和動(dòng)力學(xué)過(guò)程，因而在求解時(shí)，需將動(dòng)力學(xué)項(xiàng)與物理輸運(yùn)項(xiàng)脫耦.若對(duì)物理輸運(yùn)的質(zhì)量守恒方程與鹽度方程采用相同形式，則對(duì)流擴(kuò)散輸運(yùn)方程為:

當(dāng)C表示保守物質(zhì)的濃度時(shí)，則只考慮物理輸運(yùn)，可用式（2）計(jì)算，其數(shù)值解法和水動(dòng)力模型中鹽度的質(zhì)量守恒方程相同.岸界和垂直邊界采用無(wú)通量條件，初始濃度假定北黃海（渤海的邊界）濃度為0.5單位，渤海全海區(qū)為1單位分布，水質(zhì)模型與水動(dòng)力模型設(shè)置網(wǎng)格相同，時(shí)間步長(zhǎng)為900s.

1.2 水交換作用下COD累積量計(jì)算模型

排入海中污染物會(huì)通過(guò)水體對(duì)流產(chǎn)生輸運(yùn)和稀釋擴(kuò)散，進(jìn)而不斷與外部海域進(jìn)行交換，其間不斷有剩余的污染物累積下來(lái).如果污染物排海量逐年減少，那么這一累積量便會(huì)減少，反之，則會(huì)逐年增加.

對(duì)于渤海，水交換是由沿岸向中部水域進(jìn)行的（圖1），根據(jù)水半交換時(shí)間的定義，可構(gòu)建出污染物累積量數(shù)學(xué)模型.由水交換的遞進(jìn)過(guò)程，該模型包含兩個(gè)部分:一部分是陸源污染物進(jìn)入渤海沿岸水域的累積量，見(jiàn)式（3）；另一部分是污染物經(jīng)渤海中部與渤海外部實(shí)現(xiàn)水交換的累積量，見(jiàn)式（4）.

式中:i=1，2，…，N-1表示渤海沿岸水域分區(qū)；Ai，n表示沿岸水域分區(qū)i第n個(gè)單位時(shí)間的污染物累積量；ai，n表示沿岸水域分區(qū)i第n個(gè)單位時(shí)間污染物排放量；Hi表示沿岸水域分區(qū)i的水半交換時(shí)間；t表示積分變量.

式中: AN，n表示渤海中部（即分區(qū)N）第n個(gè)時(shí)間單位的污染物累積量；HN表示渤海中部水半交換時(shí)間.

1.3 COD降解計(jì)算模型

根據(jù)渤海灣模擬實(shí)驗(yàn)研究成果［12］，COD的降解過(guò)程基本符合一級(jí)反應(yīng)動(dòng)力學(xué)規(guī)律，其衰減速率方程為:

式中:C表示t時(shí)刻COD濃度，mg/L；t表示反應(yīng)時(shí)間，d；C0表示最初時(shí)刻COD濃度，mg/L； k表示衰減速率，d-1，即單位生物量在單位時(shí)間內(nèi)完成降解的污染量.

1.4 水交換和降解過(guò)程共同作用下COD殘余量計(jì)算

結(jié)合式（3）、式（4）和式（5），可推導(dǎo)出渤海各分區(qū)水交換和降解共同作用下COD殘余量計(jì)算公式.設(shè)有:

式中:Ei，n表示分區(qū)i第n個(gè)單位時(shí)間在水交換作用下COD經(jīng)交換減少的量；Di，n表示分區(qū)i第n個(gè)單位時(shí)間在生化作用下COD降解的量；Si，n表示分區(qū)i第n個(gè)單位時(shí)間水交換和降解共同作用下COD殘余量；ai，n表示分區(qū)i第n個(gè)單位時(shí)間污染物排放量；Hi表示分區(qū)i半交換時(shí)間；t表示積分變量.

將式（12）～（13）累加后，即得到渤海COD的殘余量

式中:Sn表示渤海COD殘余量.計(jì)算時(shí)，可以假定開(kāi)始COD排放量為零，這樣處理，在經(jīng)過(guò)多個(gè)半交換周期以后，對(duì)最終殘余量的計(jì)算結(jié)果并無(wú)較大影響.

2 模型求解

2.1 水半交換時(shí)間

根據(jù)渤海岸線的輪廓、地形和水域的氣候條件及其獨(dú)特的水流循環(huán)結(jié)構(gòu)，將渤海劃分為萊州灣、渤海灣、遼東灣和渤海中部4個(gè)分區(qū)，應(yīng)用計(jì)算機(jī)軟件MIKE21進(jìn)行深海潮汐及海流的模擬，并與實(shí)測(cè)的潮汐和潮流資料作比較后，得到渤海水域的流場(chǎng)見(jiàn)圖1.自圖1可知，渤海4個(gè)分區(qū)水交換過(guò)程是先由萊州灣、渤海灣和遼東灣交換到渤海中部，然后再由渤海中部交換到黃海（圖2）.各分區(qū)水體參數(shù)見(jiàn)表1.

根據(jù)魏浩［8］的結(jié)果，以平均氣象條件和協(xié)振潮、外部流驅(qū)動(dòng)得到的環(huán)流場(chǎng)輸入水質(zhì)模型，同時(shí)考慮風(fēng)、潮汐和水密度梯度的作用，積分3年分別計(jì)算出各灣水域平均濃度為初始濃度一半的時(shí)間，即萊州灣水體完成一次半交換的時(shí)間需0.52a；渤海灣需0.83a；遼東灣需2.26a；渤海中部需1.38a.

2.2 衰減速率求解

為得到渤海COD衰減速率，先根據(jù)渤海近10年COD的排放量的統(tǒng)計(jì)值，利用式（8）～（14）試算出不同k值下渤海整體的累積濃度，將計(jì)算結(jié)果與渤海實(shí)測(cè)的COD平均濃度進(jìn)行對(duì)比，其誤差最小時(shí)的k值即為所求.

圖1 渤海海區(qū)流場(chǎng)模擬Fig.1 Simulation of tidal flow field in the Yellow Sea and Bohai Sea

圖2 渤海水域分區(qū)及水交換方向Fig.2 Partitions of Bohai Sea and the illusion of the water exchange

表1 渤海水域各分區(qū)參數(shù)Table 1 Parameters of the zones of the Bohai Sea

（1） 渤海各分區(qū)COD排放量統(tǒng)計(jì)

按照地理劃分，萊州灣沿岸有山東省的東營(yíng)、煙臺(tái)、濰坊和濱州4個(gè)城市；渤海灣沿岸有滄州、天津和唐山3個(gè)城市；遼東灣沿岸有葫蘆島、錦州、盤錦、營(yíng)口、大連和秦皇島6個(gè)城市.設(shè)排入海域污染物主要以陸源污染統(tǒng)計(jì)為基礎(chǔ)，2004～2013年統(tǒng)計(jì)值數(shù)據(jù)來(lái)源于我國(guó)環(huán)境統(tǒng)計(jì)公報(bào)和海洋統(tǒng)計(jì)年鑒，當(dāng)個(gè)別數(shù)據(jù)缺失時(shí)，采用內(nèi)插值或以多年平均值補(bǔ)齊［13-15］.由于自2011年開(kāi)始，我國(guó)環(huán)境保護(hù)部對(duì)統(tǒng)計(jì)方法及相關(guān)技術(shù)規(guī)定等進(jìn)行了修訂，COD的統(tǒng)計(jì)范圍擴(kuò)展為工業(yè)源、農(nóng)業(yè)源和城鎮(zhèn)生活源.為便于分析，將2004～2013年的COD排放統(tǒng)計(jì)值按新的口徑統(tǒng)一進(jìn)行換算，結(jié)果見(jiàn)表2.

（2） 渤海COD濃度實(shí)際觀測(cè)

目前，我國(guó)在渤海共設(shè)環(huán)境觀測(cè)站160個(gè)，其中，國(guó)家環(huán)保局140個(gè)，國(guó)家海洋局有20個(gè)，根據(jù)我國(guó)海洋統(tǒng)計(jì)年鑒統(tǒng)計(jì)，2004～2007年以及2009年渤海COD的實(shí)測(cè)平均濃度依次為1.04、1.18、1.21、1.50、1.36mg/L.

（3） 渤海COD衰減速率計(jì)算

利用COD殘余量計(jì)算結(jié)果與實(shí)測(cè)COD濃度進(jìn)行比較，便可確定渤海水域COD衰減速率k值.經(jīng)試算，當(dāng)衰減速率k取為0.061時(shí)，兩者誤差為4.86%，此點(diǎn)為誤差最小值（表3）.因此，本文選取k=0.061作為渤海各分區(qū)COD濃度分析的計(jì)算參數(shù).不同k值下COD濃度實(shí)測(cè)值與計(jì)算值的平均相對(duì)誤差見(jiàn)圖3.

表2 2004～2013年渤海各水域COD排放量（×104t）Table 2 COD emissions of the bays in the Bohai Sea from 2004 to 2013 （×104t）

表3 渤海COD濃度實(shí)測(cè)值與計(jì)算值比較（k=0.061）Table 3 The comparison between the calculated value and the measured COD concentrations （k=0.061）

圖3 不同k值下COD濃度實(shí)測(cè)值與計(jì)算值的相對(duì)誤差Fig.3 Relative error between the computed COD concentration and the measured data for different degradation coefficient

2.3 COD濃度計(jì)算

根據(jù)渤海近岸各灣COD排放量（表2），當(dāng)k=0.061時(shí)，由式（8）～（14）可得到2004～2013年渤海各分區(qū)COD濃度（表4）.

表4 2004～2013年渤海各水域COD累計(jì)濃度（mg/L）Table 4 Accumulated COD concentrations of the bays in the Bohai Sea from 2004 to 2013 （mg/L）

3 不同水質(zhì)標(biāo)準(zhǔn)下各灣的COD容納余量

3.1 COD污染狀況

根據(jù)我國(guó)國(guó)家海水水質(zhì)標(biāo)準(zhǔn)（GB3097-1997），COD含量小于等于2mg/L屬于清潔水質(zhì)；大于2mg/L，小于等于3mg/L屬于較清潔水質(zhì)；大于3mg/L，小于等于4mg/L屬于輕度污染水質(zhì)；大于4mg/L，小于等于5mg/L屬于中度污染水質(zhì)；大于5mg/L屬于重度污染水質(zhì).對(duì)比表3，可以判斷出在2004年至2013年間，萊州灣COD污染程度屬于較清潔，但呈上升趨勢(shì)；渤海灣COD濃度由較清潔一度上升為中度污染，后又降到輕度污染，總體上升幅度較大；遼東灣COD濃度由清潔上升為較清潔，呈直線上升趨勢(shì)，詳見(jiàn)圖4.

3.2 COD環(huán)境容納余量

環(huán)境容納余量是指在一定標(biāo)準(zhǔn)下，某一環(huán)境容納污染物后還剩余的負(fù)荷容量，反映了環(huán)境容納污染物的能力.根據(jù)定義，不同海水水質(zhì)標(biāo)準(zhǔn)下渤海各分區(qū)COD環(huán)境容納余量計(jì)算公式為:

式中:i=1，2，3分別表示萊州灣、渤海灣和遼東灣；j=1，2，3，4分別表示清潔、較清潔、輕度污染和中度污染；Qi是分區(qū)i的COD環(huán)境余量是水質(zhì)標(biāo)準(zhǔn)j下的COD濃度上限，如=2.00，以此類推；是分區(qū)i的COD濃度值.

圖4 渤海各分區(qū)COD濃度的年度變化Fig.4 Annual COD concentrations evolution process for different bays in the Bohai Sea

利用式（15），可得到2013年萊州灣、渤海灣和遼東灣在相應(yīng)水質(zhì)標(biāo)準(zhǔn)區(qū)間內(nèi)COD的環(huán)境容納余量，即在較清潔區(qū)間內(nèi)萊州灣環(huán)境容納余量為26%；在輕度污染區(qū)間內(nèi)渤海灣的環(huán)境容納余量為5%；在較清潔區(qū)間內(nèi)遼東灣的環(huán)境容納余量為88%.

4 結(jié)語(yǔ)

以渤海沿岸的萊州灣、渤海灣和遼東灣以及渤海中部為不同水域分區(qū)，以水動(dòng)力交換和降解作用兩個(gè)方面分別分析計(jì)算，即利用不同水域分區(qū)水半交換周期構(gòu)建出COD累積量計(jì)算模型，計(jì)算出COD的擴(kuò)散能力；同時(shí)，利用COD衰減速率方程計(jì)算渤海水域的COD降解能力，綜合后得到渤海沿岸各分區(qū)在不同水質(zhì)標(biāo)準(zhǔn)下COD的環(huán)境容納余量.結(jié)果發(fā)現(xiàn)，盡管萊州灣、渤海灣及遼東灣同屬渤海，但COD污染情況卻相差較大.到2013年，萊州灣水質(zhì)的COD污染程度屬于較清潔，在此標(biāo)準(zhǔn)下還剩余26%的環(huán)境余量；渤海灣屬于輕度污染，在此標(biāo)準(zhǔn)下還剩余5%的環(huán)境余量；遼東灣屬于較清潔，在此標(biāo)準(zhǔn)下還剩余88%的環(huán)境余量.

上述分析結(jié)果對(duì)于有效制定環(huán)境保護(hù)措施具有重要參考價(jià)值.在治理渤海水域COD污染時(shí)，對(duì)萊州灣、渤海灣和遼東灣應(yīng)實(shí)行區(qū)別對(duì)待，尤其是近期應(yīng)嚴(yán)格控制渤海灣的COD排放量.分析中還發(fā)現(xiàn)，由于渤海水域面積較大，對(duì)沿岸海域僅按3個(gè)分區(qū)開(kāi)展污染容量研究還顯不夠，尤其是遼東灣面積過(guò)大，所以對(duì)水質(zhì)污染程度的分析結(jié)果相對(duì)該灣北部營(yíng)口一帶水域的實(shí)際污染狀況過(guò)于樂(lè)觀.因此，還可以按本文思路將研究水域范圍進(jìn)一步縮小，以便得到更貼近實(shí)際污染狀況的結(jié)果，這可作為下一步深入研究的內(nèi)容.

［1］Luff R， Pohlmann T. Calculation of water exchange times in the ICES-boxes with a eulerian dispersion model using a half-life time approach ［J］. Deutsche Hydrografische Zeitschrift，1995，47（4）:287-299.

［2］董禮先，蘇紀(jì)蘭.象山港水交換數(shù)值研究Ⅰ.對(duì)流-擴(kuò)散型的水交換模式 ［J］. 海洋與湖沼， 1999，30（4）:410-415.

［3］姚炎明，彭 輝，杜雅杰，等.象山港分區(qū)水交換數(shù)值研究 ［J］. 海洋學(xué)報(bào)（中文版）， 2014，36（1）:126-130.

［4］張宇銘，宋朝陽(yáng)，吳克儉，等.環(huán)渤海排污口臨近海域水交換能力研究 ［J］. 中國(guó)海洋大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2014，44（5）:1-7.

［5］Sun J， Tao J H. Relation matrix of water exchange for sea bays and its application ［J］. China Ocean Engineering， 2006，20（4）: 529-544.

［6］Wei H， Zhao L， Feng S. Comparison of the tide-induced Lagrangian and Eulerian mean circulation in the Bohai Sea ［J］. Chinese Limin and Oceanogra， 2001，19（2）:121-132.

［7］Wei H， Hainbucher D， Pohlmann T， et al. Tidal-induced Lagrangian and Eulerian mean circulation in the Bohai Sea ［J］. Journal of Marine Systems， 2004，44（4）:141-151.

［8］魏 皓，田 恬，周 鋒，等.渤海水交換的數(shù)值研究-水質(zhì)模型對(duì)半交換時(shí)間的模擬 ［J］. 青島海洋大學(xué)學(xué)報(bào)（自然科學(xué)版），2002，32（4）:519-525.

［9］邱 巍.長(zhǎng)江口竹園排污區(qū)COD降解系數(shù)的測(cè)試與分析 ［J］.上海水利， 1996，（4）:33-36+12.

［10］郭棟鵬，徐明德.黃海南部海水中COD降解規(guī)律的研究 ［J］. 太原理工大學(xué)學(xué)報(bào)， 2008，39（4）:358-361.

［11］Zhao X X， Wang X L， Shi X Y， et al. Environmental capacity of 18chemical oxygen demand in the Bohai Sea: modeling and calculation ［J］. Chinese Journal of 19Oceanology and Limnology，2011，29（1）:46-52.

［12］喬璐璐，劉容子，鮑獻(xiàn)文，等.經(jīng)濟(jì)增長(zhǎng)下的渤海環(huán)境容量預(yù)測(cè)［J］. 中國(guó)人口·資源與環(huán)境， 2008，18（2）:76-81.

［13］蘇一兵，雷 坤，孟 偉.陸域活動(dòng)對(duì)渤海海岸帶的影響 ［J］. 中國(guó)水利， 2003，（3）:78-80.

［14］王修林，崔正國(guó)，李克強(qiáng).渤海COD入海通量估算及其分配容量?jī)?yōu)化研究 ［J］. 海洋環(huán)境科學(xué)， 2009，28（5）:497-500.

［15］渤海環(huán)境保護(hù)總體規(guī)劃（2008-2020） ［Z］.

COD environmental capacity of different zones of the Bohai Sea with the consideration of self-purification.

YU Tiao-lan， WANG Nuo*， LIU Zhong-bo， LI Hui （College of Transportation Management， Dalian Maritime University，Dalian 116026， China）. China Environmental Science， 2015，35（5）：1579～1585

To determine the partial water quality， it is essential to divide large water area into smaller partitions according to the geographical structure. In present paper， China's Bohai Sea was divided into four zones （i.e.， Laizhou Bay， Bohai Bay， Liaodong Bay and the center of Bohai Sea）. The hydrodynamic induced mass convection and diffusion and the chemical degradation were considered， the former was used to compute the diffusion capability of cumulative chemical oxygen demand （COD） in different zones， and the latter was applied to solve COD decreasing property， thus the degradation ability could be estimated. After the comprehensive consideration of these two parts， the COD environmental capacity was obtained with different water quality standards considered. The result shows that COD capacity is different in these four areas. For Laizhou Bay， the COD concentrations could satisfy the standard of “relatively clean”， and the environmental capacity has additional 26% to reach the standard； For Laizhou Bay， the COD concentrations was very close to the standard of “l(fā)ight pollution”， and the environmental capacity only has 5% left. For Liaodong Bay， COD concentrations met the standard of “relatively clean”， and the environmental capacity remained 88%.

sea；environment；capacity；COD；Bohai Sea

X145

A

100-6923（2015）05-1579-07

郁斢蘭（1988-），女，江西吉安人，大連海事大學(xué)博士研究生，主要研究方向?yàn)榻煌üこ膛c環(huán)境.

2014-10-05

教育部人文社會(huì)科學(xué)研究規(guī)劃基金項(xiàng)目（12YJA630128）

* 責(zé)任作者， 教授， wangnuodl@126.com