曹 維，仵彥卿，唐德富，楊耀芳，沈繼平，費(fèi)岳軍

（1.國(guó)家海洋局 寧波海洋環(huán)境監(jiān)測(cè)中心站，浙江 寧波 315012；2.上海交通大學(xué) 環(huán)境科學(xué)與工程學(xué)院，上海 200240）

0 引言

象山港位于寧波市東南部，穿山半島與象山半島之間，東臨大目洋，是一個(gè)由東北向西南深入內(nèi)陸的狹長(zhǎng)型半封閉海灣，全港縱深60km，水深為10～15 m，港灣內(nèi)有西滬港、黃墩港和鐵港3個(gè)支港。象山港屬內(nèi)陸海港，為半封閉型海灣，盡管有潮汐的作用，但其海水交換能力差，不利于污染物的降解和擴(kuò)散，環(huán)境容量較?。?］。

20世紀(jì)70年代初期，寧波市衛(wèi)生防疫站負(fù)責(zé)象山港水質(zhì)監(jiān)測(cè)任務(wù)［2］；1976年起，寧波市環(huán)境監(jiān)測(cè)中心站承擔(dān)了象山港水質(zhì)監(jiān)測(cè)任務(wù)，設(shè)置了5個(gè)站位（吉奇、黃墩港、下山、黃避岙和野龍山）；1983年和1985年在上述5個(gè)站位增測(cè)了沉積物樣品；2002年，國(guó)家海洋局在象山港海域設(shè)置了赤潮監(jiān)控區(qū)［3］。受監(jiān)測(cè)能力和監(jiān)測(cè)經(jīng)費(fèi)的影響，象山港監(jiān)測(cè)站位數(shù)量有限，站位分布也不具代表性。隨著我國(guó)對(duì)海洋環(huán)境保護(hù)工作的重視，港灣海域監(jiān)測(cè)工作有待進(jìn)一步規(guī)范和深化，因此，有必要對(duì)監(jiān)測(cè)站位進(jìn)行優(yōu)化。本研究利用2011年度海洋環(huán)境評(píng)價(jià)項(xiàng)目《象山港海域海洋環(huán)境質(zhì)量綜合評(píng)價(jià)方法 DOMEP（MEA）－03－02》專題調(diào)查監(jiān)測(cè)數(shù)據(jù)，結(jié)合同期赤潮監(jiān)控區(qū)和趨勢(shì)性監(jiān)測(cè)資料①，進(jìn)行站位優(yōu)化和驗(yàn)證分析，希望通過優(yōu)化后的監(jiān)測(cè)站位全面掌握港灣環(huán)境質(zhì)量，更好地為海洋環(huán)境保護(hù)管理提供決策依據(jù)。

1 站位優(yōu)化方法

目前關(guān)于監(jiān)測(cè)站位優(yōu)化的新方法，主要有物元分析法［4］、灰色局勢(shì)決策分析法［5］、序貫數(shù)論法［6］、模糊優(yōu)選法［7］、密切值法［8］、模糊聚類法［9］和人工神經(jīng)網(wǎng)絡(luò)優(yōu)選法［10］等，但上述方法在一定程度上存在缺陷。

克立金方法（Kriging）［11］是 G.Matheron教授以南非礦山地質(zhì)工程師D.G.Krige的名字命名的一種方法。從數(shù)學(xué)上講，克立金法是一種對(duì)空間分布數(shù)據(jù)

①寧波市海洋與漁業(yè)局.2011年寧波市海洋與漁業(yè)環(huán)境監(jiān)測(cè)成果匯編（內(nèi)部資料）［G］.寧波：寧波市海洋與漁業(yè)局，2012.求最優(yōu)、線性、無偏內(nèi)插估計(jì)量的方法。從海洋環(huán)境研究角度講，它是根據(jù)已知監(jiān)測(cè)站位點(diǎn)上的環(huán)境要素變量，如某種污染物濃度的實(shí)測(cè)數(shù)據(jù)，對(duì)環(huán)境要素變量進(jìn)行結(jié)構(gòu)性分析（變差函數(shù)模型的確定）之后，為了對(duì)待估點(diǎn)作出一種線性、無偏、最小方差的估計(jì)，而對(duì)周圍已知站位點(diǎn)的測(cè)量值賦予一定權(quán)系數(shù)，進(jìn)行加權(quán)平均來估計(jì)待估點(diǎn)環(huán)境要素變量的方法。

本次站位優(yōu)化設(shè)計(jì)運(yùn)用普通Kriging方差分析法，采用下式計(jì)算研究區(qū)海洋環(huán)境要素的估計(jì)誤差標(biāo)準(zhǔn)差，同時(shí)結(jié)合研究區(qū)污染現(xiàn)狀，進(jìn)行海洋環(huán)境監(jiān)測(cè)站位優(yōu)化分析：

式中：σ為估計(jì)誤差的標(biāo)準(zhǔn)差；λi（i＝1，2，...，N）為克立金權(quán)系數(shù)；N為已知監(jiān)測(cè)點(diǎn)總數(shù)；μ是拉格朗日乘子；xi＝（xi，yi），為已知監(jiān)測(cè)i點(diǎn)的坐標(biāo)；x0＝（x0，y0），為待插值點(diǎn)的坐標(biāo)；γ（xi，x0）＝γ（h）為變差函數(shù)，h為（x0，y0）到（xi，yi）之間的距離矢量，于是，γ（h）可用下列公式計(jì)算：

求得變差函數(shù)后，代入下式可以計(jì)算權(quán)系數(shù)λi（i＝1，2，...，N）和拉格朗日乘子μ：

將式（3）求得λi（i＝1，2，...，N）和μ代入到式（1），可以得到所有未知點(diǎn)的方差值。該方差的大小可反映已有監(jiān)測(cè)站位質(zhì)量的優(yōu)劣。

利用方差減少法，進(jìn)行現(xiàn)有監(jiān)測(cè)站位的優(yōu)化分析，站位優(yōu)化調(diào)整的原則為：

（1）污染物擴(kuò)散濃度梯度大的區(qū)域，加大監(jiān)測(cè)站位密度，否則，減少站位密度；

（2）考慮到海洋水動(dòng)力條件，在流速梯度大的區(qū)域，加大監(jiān)測(cè)站位密度，否則，減少站位密度；

（3）污染源分布的地區(qū)，需要加密站位；

（4）對(duì)于有長(zhǎng)序列的站位，盡可能保留；

（5）對(duì)海洋近岸地區(qū)，加大站位密度，遠(yuǎn)海區(qū)域減少站位密度；

（6）監(jiān)測(cè)站位優(yōu)化是一個(gè)過程，與經(jīng)費(fèi)投入和環(huán)境信息提取精度有關(guān)，因此，在站位優(yōu)化時(shí)，還要考慮目前的海洋環(huán)境研究狀況和經(jīng)費(fèi)投入。

2 站位優(yōu)化實(shí)例

2.1 數(shù)據(jù)來源

研究采用的水質(zhì)、沉積物質(zhì)量、生物生態(tài)數(shù)據(jù)主要來源于2011年度海洋環(huán)境評(píng)價(jià)項(xiàng)目專題監(jiān)測(cè)數(shù)據(jù)，結(jié)合同期象山港赤潮監(jiān)控區(qū)和趨勢(shì)性監(jiān)測(cè)數(shù)據(jù)；水文（水溫和鹽度）、沉積物類型采用2011年度海洋環(huán)境評(píng)價(jià)項(xiàng)目專題調(diào)查數(shù)據(jù)。

2.2 現(xiàn)有監(jiān)測(cè)站位

象山港現(xiàn)有13個(gè)赤潮監(jiān)控區(qū)站位、6個(gè)趨勢(shì)性監(jiān)測(cè)站位、31個(gè)專題調(diào)查站位（港口7個(gè)、港中11個(gè)、港底13個(gè)），共有50個(gè)監(jiān)測(cè)站位（圖1）。

2.3 站位優(yōu)化結(jié)果

根據(jù)普通Kriging方差分析法，計(jì)算得到象山港現(xiàn)有監(jiān)測(cè)站位水質(zhì)要素估計(jì)誤差標(biāo)準(zhǔn)差平均值為0.190（圖2），灣底和灣中部的估計(jì)誤差標(biāo)準(zhǔn)差較小，而在灣口區(qū)域的估計(jì)誤差標(biāo)準(zhǔn)差較大。此外，在觀測(cè)點(diǎn)較密的區(qū)域，估計(jì)誤差標(biāo)準(zhǔn)差??；反之，估計(jì)誤差標(biāo)準(zhǔn)差大。

同時(shí)，通過象山港海域水質(zhì)各監(jiān)測(cè)要素值的標(biāo)準(zhǔn)化和歸一化處理，進(jìn)行現(xiàn)有水質(zhì)污染狀況分析?？紤]到不同時(shí)期數(shù)據(jù)的差異性，采用31個(gè)專題調(diào)查站位的水質(zhì)監(jiān)測(cè)數(shù)據(jù)進(jìn)行水質(zhì)污染狀況分析（圖3）?？梢缘贸?，象山港灣底污染程度較嚴(yán)重，灣中部區(qū)域次嚴(yán)重，灣口區(qū)域較輕。

結(jié)合象山港現(xiàn)有監(jiān)測(cè)站位水質(zhì)估計(jì)誤差標(biāo)準(zhǔn)差和象山港現(xiàn)有水質(zhì)監(jiān)測(cè)污染程度（圖2和圖3），對(duì)象山港水質(zhì)監(jiān)測(cè)站位進(jìn)行優(yōu)化。優(yōu)化后象山港有38個(gè)水質(zhì)測(cè)點(diǎn)（圖4），刪減了12個(gè)近岸站點(diǎn)，優(yōu)化后的水質(zhì)要素平均估計(jì)誤差標(biāo)準(zhǔn)差為0.168（圖5），增加了11.6%，刪減的監(jiān)測(cè)站位包括5個(gè)赤潮監(jiān)測(cè)站位（XS05、XS07、XS08、XS09和XS12）、2個(gè)趨勢(shì)性監(jiān)測(cè)站位（D33ZQ037和D33ZQ038）以及5個(gè)專題調(diào)查站監(jiān)測(cè)站位（QS2、QS6、QS10、QS22和 QS27）。

水質(zhì)站位設(shè)置：根據(jù)象山港豐水期和枯水期水體中各類污染物分布特點(diǎn)（圖6），充分考慮能監(jiān)測(cè)各類污染物重點(diǎn)分布區(qū)域，同時(shí)統(tǒng)籌考慮象山港水團(tuán)（圖7）和水動(dòng)力［12］（圖8）分布區(qū)域，優(yōu)化后得到38個(gè)水質(zhì)監(jiān)測(cè)站位（圖4）。

圖1 象山港現(xiàn)有監(jiān)測(cè)站位Fig.1 Present monitoring stations in Xiangshangang Bay

圖2 現(xiàn)有站位水質(zhì)估計(jì)誤差標(biāo)準(zhǔn)差Fig.2 Water quality estimated standard deviation of present stations

圖3 現(xiàn)有水質(zhì)污染狀況Fig.3 Present pollution condition of water quality

圖4 優(yōu)化后的監(jiān)測(cè)站位Fig.4 Optimized monitoring stations

圖5 優(yōu)化站位水質(zhì)估計(jì)誤差標(biāo)準(zhǔn)差Fig.5 Water quality estimated standard deviation of optimized stations

圖6 水質(zhì)污染物分布Fig.6 Water pollutant distribution

圖7 水團(tuán)分布Fig.7 Water mass distribution

圖8 水體半交換時(shí)間Fig.8 Half exchange time of water

沉積物站位設(shè)置：原則上沉積物站位數(shù)為優(yōu)化后水質(zhì)監(jiān)測(cè)站位數(shù)的50%，合理布局，均勻分布。同時(shí)，充分考慮布點(diǎn)覆蓋所有沉積物類型（圖9），且各沉積物類型上設(shè)點(diǎn)密度基本保持一致，尤其是獅子口北面的QS26號(hào)站位因其僅有的貝殼－砂沉積物類型，故予以保留。對(duì)于有長(zhǎng)序列的站位，盡可能保留，在站位設(shè)置中保留了具有多年沉積物監(jiān)測(cè)數(shù)據(jù)的趨勢(shì)性監(jiān)測(cè)站位（D33ZQ035）。優(yōu)化后共得到19個(gè)沉積物監(jiān)測(cè)站位（圖4）。

圖9 沉積物類型分布Fig.9 Distribution of sediment types

生物生態(tài)站位設(shè)置：充分考慮象山港水團(tuán)分布特征，浮游生物和底棲生物的群落分布特征，使站位覆蓋各類型水團(tuán)和各生物群落類型。根據(jù)水質(zhì)優(yōu)化站位來看，已覆蓋所有類型的水團(tuán)和浮游生物、底棲生物群落，生物生態(tài)站位數(shù)在取60%水質(zhì)站位數(shù)的過程中，做到均勻、有序地刪減。根據(jù)浮游植物優(yōu)勢(shì)種類的生態(tài)類型和分布特點(diǎn)，將象山港浮游植物進(jìn)行生態(tài)類型劃分。夏季象山港可劃分為沿岸種分布區(qū)、近岸溫帶種分布區(qū)和半咸水種分布區(qū)3個(gè)主分布區(qū)，在獅子口海域有少量熱帶種出現(xiàn)（圖10a）。冬季象山港可分為外洋廣溫種分布區(qū)、近岸暖海種分布區(qū)、廣溫廣布種分布區(qū)和半咸水種分布區(qū)4個(gè)區(qū)域（圖10b）。底棲生物站位應(yīng)充分考慮布點(diǎn)覆蓋所有沉積物類型，且各沉積物類型上設(shè)點(diǎn)密度基本保持一致，根據(jù)此原則，在獅子口北面的貝殼－砂沉積物類型中，保留QS26號(hào)站位，以充分覆蓋和均勻分布各底質(zhì)類型。為了盡可能保留長(zhǎng)序列的站位，在站位設(shè)置中保留了具有多年生物監(jiān)測(cè)數(shù)據(jù)的趨勢(shì)性監(jiān)測(cè)站位（D33ZQ035）。優(yōu)化后共為23個(gè)生物生態(tài)監(jiān)測(cè)站位（圖4）。

圖10 浮游植物類型分布Fig.10 Distribution of phytoplankton communities

3 驗(yàn)證分析

為了驗(yàn)證優(yōu)化后監(jiān)測(cè)站位所反映的生態(tài)環(huán)境質(zhì)量與現(xiàn)有監(jiān)測(cè)站位所反映的生態(tài)環(huán)境質(zhì)量之間是否存在差異，需對(duì)2組數(shù)據(jù)進(jìn)行顯著性檢驗(yàn)。在數(shù)理統(tǒng)計(jì)中一般以概率（P）0.05作為顯著評(píng)定標(biāo)準(zhǔn)，即在100次試驗(yàn)中，由于偶然因素造成差異的可能性在5次以上，其差異被認(rèn)為是不顯著。如果兩者差異在概率為5%的范圍內(nèi)，出現(xiàn)這樣概率的機(jī)會(huì)非常小而出現(xiàn)了，那么我們就認(rèn)為此差數(shù)具有顯著差異程度。

根據(jù)2011年象山港生態(tài)環(huán)境專題調(diào)查結(jié)果，結(jié)合同年在象山港開展的赤潮監(jiān)控區(qū)、趨勢(shì)性監(jiān)測(cè)等相關(guān)資料，采用統(tǒng)計(jì)學(xué)中的T檢驗(yàn)方法，對(duì)站位優(yōu)化前、后象山港海域的水質(zhì)、沉積物和生物生態(tài)監(jiān)測(cè)數(shù)據(jù)進(jìn)行顯著性檢驗(yàn)。其中，水質(zhì)采用溶解氧、化學(xué)需氧量、活性磷酸鹽、無機(jī)氮、石油類、重金屬（汞、砷、銅、鉛、鋅、鎘、鉻）等12個(gè)指標(biāo)的監(jiān)測(cè)值；沉積物采用石油類、滴滴涕、多氯聯(lián)苯、重金屬（汞、砷、銅、鉛、鋅、鎘、鉻）等10個(gè)指標(biāo)的監(jiān)測(cè)數(shù)據(jù)；生物生態(tài)采用浮游植物、浮游動(dòng)物和底棲生物的多樣性指數(shù)。經(jīng)統(tǒng)計(jì)檢驗(yàn)（表1），站位優(yōu)化前、后的2組數(shù)據(jù)之間無顯著性差異，說明優(yōu)化后的監(jiān)測(cè)站位（38個(gè)水質(zhì)站位、19個(gè)沉積物站位、23個(gè)生物生態(tài)站位）所表征的監(jiān)測(cè)信息量與現(xiàn)有50個(gè)監(jiān)測(cè)站位的信息量等效，可作為今后的常規(guī)監(jiān)測(cè)站位。

表1 站位優(yōu)化前、后監(jiān)測(cè)數(shù)據(jù)顯著性檢驗(yàn)Tab.1 Significant testing of original and optimized monitoring data

4 小結(jié)

根據(jù)2011年象山港海域生態(tài)環(huán)境監(jiān)測(cè)資料，運(yùn)用普通Kriging方差分析法進(jìn)行海洋環(huán)境監(jiān)測(cè)站位優(yōu)化，再采用統(tǒng)計(jì)學(xué)中的T檢驗(yàn)方法，對(duì)站位優(yōu)化前、后象山港海域生態(tài)監(jiān)測(cè)數(shù)據(jù)進(jìn)行方差分析，經(jīng)統(tǒng)計(jì)檢驗(yàn)，2組數(shù)據(jù)之間無顯著性差異，說明優(yōu)化后的38個(gè)監(jiān)測(cè)站位能較好地代表象山港海域的生態(tài)環(huán)境質(zhì)量，其所表征的監(jiān)測(cè)信息量與現(xiàn)有50個(gè)監(jiān)測(cè)站位的信息量等效，故可作為今后的常規(guī)監(jiān)測(cè)站位。

（References）：

［1］LI Nai－fang，LI Ying－bao.Discussion on resource characteristics and development mode of Xiangshan port［J］.Ningbo Teachers College：Social Science Edition，1986（4）：101－104.

李乃芳，李英保.試論象山港的資源特征和開發(fā)方式［J］.寧波師范學(xué)院：社會(huì)科學(xué)版，1986（4）：101－104.

［2］HU Wen－xiang，CHEN Tie－rong，XIN Ying，et al.The assessment on the environment monitoring of XiangSanPort［J］.Marine Environmental Science，1995，14（4）：57－63.

胡文翔，陳鐵熔，忻穎，等.象山港環(huán)境監(jiān)測(cè)總結(jié)評(píng)價(jià)［J］.海洋環(huán)境科學(xué)，1995，14（4）：57－63.

［3］ZHANG Li－xu，JIANG Xiao－shan，CAI Yan－h(huán)ong.Variations on nutrients and the characteristic of their structure in the red tidemonitoring area of Xiangshan Habor in recent four years［J］.Marine Sciences Bulletin，2006，25（6）：1－9.

張麗旭，蔣曉山，蔡燕紅.近4年來象山港赤潮監(jiān)控區(qū)營(yíng)養(yǎng)鹽變化及其結(jié)構(gòu)特征［J］.海洋通報(bào)，2006，25（6）：1－9.

［4］GAO Ming－h(huán)ui.Reseach on optimized selection in water quality environmental monitoring by substance elemant analysis［J］.Environmental Science and Technology，1997，20（4）：26－27.

高明慧.物元分析用于水質(zhì)環(huán)境監(jiān)測(cè)優(yōu)化布點(diǎn)的研究［J］.環(huán)境科學(xué)與技術(shù)，1997，20（4）：26－27.

［5］ZHANG Yong，HUANG Xing－qun，YI Zhi－wei，et al.Optimized selection in river monitoring points by grey situation substance elemant analysis［J］.Environmental Monitoring in China，1996，15（3）：24－26.

張勇，黃醒群，易志偉，等.灰色局勢(shì)物元分析進(jìn)行河流優(yōu)化布點(diǎn)［J］.中國(guó)環(huán)境監(jiān)測(cè)，1996，15（3）：24－26.

［6］LI Feng，LI Yong－gan.Optimized selection in river monitoring sites by sequential number－theoretic optimization method and silhouette graph［J］.Environmental Monitoring in China，2000，16（4）：15－17.

李峰，李永干.序貫數(shù)論優(yōu)化法和側(cè)影圖進(jìn)行河流優(yōu)化布點(diǎn)［J］.中國(guó)環(huán)境監(jiān)測(cè)，2000，16（4）：15－17.

［7］DOU Su－zhen，WANG Li－h(huán)ong，HOU Cun－dong，et al.A fuzzy optimization model of preferential location in urban air environmental monitoring and case study［J］.Ubran Environment ＆ Ubran Ecology，2001，14（5）：55－58.

竇素珍，王黎虹，侯存東，等.城市大氣環(huán)境監(jiān)測(cè)優(yōu)化布點(diǎn)模糊優(yōu)選模型及應(yīng)用實(shí)例［J］.城市環(huán)境與城市生態(tài)，2001，14（5）：55－58.

［8］XU Li－zhong，ZHANG Jiang－shan.Optimization of atmospheric environmental monitoring sites with modified intimate value method［J］.Environmental Engineering，2000，18（2）：50－53.

許麗忠，張江山.改進(jìn)密切值法優(yōu)化大氣環(huán)境監(jiān)測(cè)布點(diǎn)［J］.環(huán)境工程.2000，18（2）：50－53.

［9］LIANG Wei－zhen，YE Jin－run，YANG Jing.Optimized selection in water quality monitoring by fuzzy cluster analysis method［J］.Environmental Monitoring Management and Technology，2002，14（3）：6－7.

梁偉臻，葉錦潤(rùn)，楊靜.模糊聚類分析法優(yōu)化城市河涌水質(zhì)監(jiān)測(cè)點(diǎn)［J］.環(huán)境監(jiān)測(cè)管理與技術(shù)，2002，14（3）：6－7.

［10］LI Zou－yong，XU Ting－ting，DING Jing.Application of artificial neural network to optimize groundwater quality monitoring sites［J］.Ubran Environment ＆ Ubran Ecology，2003，16（6）：169－171.

李祚泳，徐婷婷，丁晶.地下水環(huán)境監(jiān)測(cè)優(yōu)化布點(diǎn)的人工神經(jīng)網(wǎng)絡(luò)模型［J］.城市環(huán)境與城市生態(tài)，2003，16（6）：169－171.

［11］KANG Yong－shang，SHEN Jin－song，CHEN Zhuo－h(huán)eng.Modern mathematics geology［M］.Beijing：Petroleum Industry Press，2005：114－131.

康永尚，沈金松，諶卓恒.現(xiàn)代數(shù)學(xué)地質(zhì)［M］.北京：石油工業(yè)出版社，2005：114－131.

［12］PENG Hui，YAO Yan－ming，LIU Lian.Study on the features of water exchange in Xiangshangang Bay［J］.Journal of Marine Sciences，2012，30（4）：1－12.

彭輝，姚炎明，劉蓮.象山港水體交換特征研究［J］.海洋學(xué)研究，2012，30（4）：1－12.