蘭冬東，李 冕，許 妍，鮑晨光，于春艷，梁 斌，馬明輝

（國(guó)家海洋環(huán)境監(jiān)測(cè)中心，遼寧 大連 116023）

?

渤海水質(zhì)監(jiān)測(cè)站位優(yōu)化研究

蘭冬東，李 冕，許 妍，鮑晨光，于春艷，梁 斌，馬明輝

（國(guó)家海洋環(huán)境監(jiān)測(cè)中心，遼寧 大連 116023）

摘要：以2000年-2010年無(wú)機(jī)氮、活性磷酸鹽和化學(xué)需氧量這三個(gè)主要污染因子趨勢(shì)性監(jiān)測(cè)數(shù)據(jù)為基礎(chǔ)，采用水質(zhì)空間穩(wěn)定性評(píng)價(jià)、監(jiān)測(cè)站位時(shí)空重要性評(píng)價(jià)以及站位網(wǎng)絡(luò)優(yōu)化設(shè)計(jì)等技術(shù)方法，對(duì)渤海2010年原有的海洋環(huán)境趨勢(shì)性水質(zhì)監(jiān)測(cè)站位進(jìn)行優(yōu)化調(diào)整，為渤海海域水質(zhì)的時(shí)空變化趨勢(shì)評(píng)價(jià)，提供更為有效、客觀的支撐。

關(guān)鍵詞：海水質(zhì)量監(jiān)測(cè)； 站位評(píng)價(jià)； 站位優(yōu)化； 渤海

環(huán)渤海地區(qū)是我國(guó)快速發(fā)展的重要經(jīng)濟(jì)集聚區(qū)和基礎(chǔ)產(chǎn)業(yè)分布區(qū)之一，為我國(guó)經(jīng)濟(jì)社會(huì)的發(fā)展提供了重要的支撐作用。隨著經(jīng)濟(jì)的快速發(fā)展，渤海面臨的壓力也在持續(xù)增大，近年來(lái)，渤海海域污染嚴(yán)重，污染范圍逐漸擴(kuò)大且出現(xiàn)新的污染區(qū)域。1974年我國(guó)在渤海首次開(kāi)展海水水質(zhì)監(jiān)測(cè)，當(dāng)年渤海全域水質(zhì)符合一類(lèi)海水水質(zhì)標(biāo)準(zhǔn)，2010年渤海一類(lèi)水質(zhì)水域面積僅占渤海總面積的 45%，劣四類(lèi)水域面積占渤?？偯娣e的比例高達(dá)18%。近年來(lái)，渤海海域水質(zhì)污染范圍逐漸擴(kuò)大，為準(zhǔn)確評(píng)價(jià)渤海水質(zhì)現(xiàn)狀及變化趨勢(shì)，在20世紀(jì)70年代初期和90年代末期兩次海洋污染基線調(diào)查的基礎(chǔ)上，渤海水質(zhì)監(jiān)測(cè)站位也在不斷調(diào)整和增加，到2010年水質(zhì)監(jiān)測(cè)站位已經(jīng)達(dá)到 151個(gè)。渤?，F(xiàn)有水質(zhì)監(jiān)測(cè)站位網(wǎng)絡(luò)是否合理，在監(jiān)測(cè)經(jīng)費(fèi)投入一定的情況下，如何設(shè)計(jì)優(yōu)化監(jiān)測(cè)站位網(wǎng)絡(luò)使得監(jiān)測(cè)評(píng)價(jià)結(jié)果更準(zhǔn)確，成為海洋環(huán)境監(jiān)測(cè)設(shè)計(jì)的重要課題。目前針對(duì)海水水質(zhì)站位優(yōu)化的研究尚屬空白，本研究對(duì)渤海水質(zhì)監(jiān)測(cè)站位的科學(xué)合理優(yōu)化，能夠?yàn)椴澈：Ｓ蛩|(zhì)的時(shí)空變化趨勢(shì)評(píng)價(jià)提供有力支撐。

1 研究方法

1.1 IDW空間插值

IDW 空間插值方法能夠合理地展示海域水質(zhì)的空間變化趨勢(shì)和主要熱點(diǎn)區(qū)域，尤其是在出現(xiàn)極高值和低值的區(qū)域，已被廣泛應(yīng)用到空間要素分析［1-4］、降雨量預(yù)測(cè)［5-6］、污染物預(yù)測(cè)［7-9］等諸多領(lǐng)域，以及海洋水質(zhì)污染評(píng)價(jià)與預(yù)測(cè)方面［10-12］，模型為：

式中，Z*（B）為待插值點(diǎn)的值，Z（xi）為插值結(jié)點(diǎn)的實(shí)測(cè)值，λi為插值結(jié)點(diǎn)xi所占的權(quán)重，其中λi的計(jì)算方法為：

式中，di為估值點(diǎn)與已知點(diǎn)間的距離，k為正數(shù)。

1.2 水質(zhì)穩(wěn)定性評(píng)價(jià)

單要素水質(zhì)穩(wěn)定性評(píng)價(jià)是對(duì)每一柵格統(tǒng)計(jì)各類(lèi)水質(zhì)出現(xiàn)的次數(shù)，出現(xiàn)次數(shù)最高的即為最可能出現(xiàn)的水質(zhì)等級(jí)，水質(zhì)穩(wěn)定概率Pi為第i類(lèi)水質(zhì)出現(xiàn)的次數(shù)除以總的統(tǒng)計(jì)次數(shù)，即：

其中，iN為水質(zhì)等級(jí)i類(lèi)出現(xiàn)的次數(shù)。

綜合水質(zhì)穩(wěn)定性評(píng)價(jià)是先分別統(tǒng)計(jì)各污染要素最可能出現(xiàn)的水質(zhì)等級(jí)及相應(yīng)概率，然后取污染等級(jí)最高的概率作為水質(zhì)綜合穩(wěn)定的概率，若污染等

1.3 站位時(shí)空重要性評(píng)價(jià)方法

1.3.1 歷史重要性站位篩選

海洋環(huán)境監(jiān)測(cè)站位的連續(xù)性對(duì)評(píng)價(jià)海洋環(huán)境變化趨勢(shì)十分重要，因此可作為確定站位重要性的標(biāo)準(zhǔn)。以2010年渤海水質(zhì)趨勢(shì)性監(jiān)測(cè)站位為基礎(chǔ)站位，篩選自2000年以來(lái)的監(jiān)測(cè)站位，通過(guò)相近站位的空間比對(duì)，確定監(jiān)測(cè)年限不低于9年的、具有歷史延續(xù)性的時(shí)間重要站位。再對(duì)距離過(guò)近的站位進(jìn)行篩選，0.1度范圍內(nèi)只保留一個(gè)站位作為歷史重要性站位。

1.3.2 站位貢獻(xiàn)率評(píng)價(jià)

綜合水質(zhì)站位貢獻(xiàn)率為站位各分指標(biāo)貢獻(xiàn)率的和。設(shè)X1是站位集X的一個(gè)子集，則X1在X中的貢獻(xiàn)率計(jì)算公式為：

X1的貢獻(xiàn)率=X的水質(zhì)評(píng)價(jià)-X/X1的水質(zhì)評(píng)價(jià)，其中X/X1=｛x| x∈X & x?X1｝，“-”表示水質(zhì)評(píng)價(jià)的差異性，即不同水質(zhì)等級(jí)的柵格數(shù)量。

1.3.3 站位綜合評(píng)價(jià)

設(shè)監(jiān)測(cè)站位總數(shù)為n個(gè)，歷史重要性站位數(shù)為m。

（1） n個(gè)監(jiān)測(cè)站位的數(shù)據(jù)水質(zhì)評(píng)價(jià)結(jié)果作為原始結(jié)果，記錄各類(lèi)水質(zhì)的柵格數(shù)量及空間位置，站位貢獻(xiàn)率計(jì)算均以此為基礎(chǔ)。

（2） 對(duì)n-m中的每一站位，計(jì)算其貢獻(xiàn)率，按升序排序，去除貢獻(xiàn)率最小的站位，記錄各站位對(duì)應(yīng)的貢獻(xiàn)率。

（3） 對(duì)剩下的 n-m-1個(gè)站位中的首站位，重新計(jì)算貢獻(xiàn)率并升序排序，若順序不變，則去除，若改變，則按排序后重新計(jì)算首站位貢獻(xiàn)率；

（4） 剩下的n-m-2個(gè)站位按步驟三依次循環(huán)。

1.4 站位重新優(yōu)化布設(shè)方法

三角剖分已經(jīng)被廣泛應(yīng)用到物理，計(jì)算機(jī)等領(lǐng)域［13-17］，考慮到海洋環(huán)境的特殊性，本文采用的是改進(jìn)的 Delaunay三角剖分進(jìn)行站位布設(shè)，即加入跨島判斷，如果兩個(gè)站位之間的連線經(jīng)過(guò)某一大島嶼或陸地，則認(rèn)為該兩站位之間無(wú)關(guān)系，刪除同時(shí)以它們倆為頂點(diǎn)的三角面。

增加的站位主要布設(shè)在污染變化劇烈的A區(qū)域以及污染變化的B區(qū)域，清潔穩(wěn)定的C區(qū)域不增加站位。先對(duì)站位進(jìn)行 Delaunay三角剖分，記錄所得三角形的邊長(zhǎng)，找出經(jīng)過(guò) A區(qū)且邊長(zhǎng)最長(zhǎng)的邊。如圖1所示，若該邊為某兩個(gè)三角形的公共邊，則新加站位為這兩個(gè)三角形四個(gè)頂點(diǎn)的加權(quán)（公式 4），若該邊不是公共邊，則新加站位為該邊中點(diǎn)。每新增一個(gè)站位則重新進(jìn)行一次Delaunay三角剖分。若新增站位落在 C區(qū)，則重新尋找次長(zhǎng)邊進(jìn)行站位添加，如此循環(huán)。B區(qū)域站位增加同A區(qū)域站位增加過(guò)程。

圖1 四頂點(diǎn)加權(quán)圖Fig. 1 Four-vertex weighted graph

2 結(jié)果與討論

2.1 渤海水質(zhì)污染與穩(wěn)定性

對(duì)2000年-2010年渤海水質(zhì)趨勢(shì)性監(jiān)測(cè)數(shù)據(jù)進(jìn)行分析整理，選取無(wú)機(jī)氮、活性磷酸鹽和化學(xué)需氧量三個(gè)渤海區(qū)域主要的污染因子進(jìn)行評(píng)價(jià)，采用表、中、底三層的平均值，地圖?。?.01度）^2為一柵格數(shù)據(jù)，采用 IDW 空間插值方法，得到歷年渤海海域水質(zhì)狀況的空間分布。如圖2所示，近年來(lái)渤海水質(zhì)污染成加重趨勢(shì)，無(wú)機(jī)氮污染增長(zhǎng)趨勢(shì)明顯，污染范圍持續(xù)擴(kuò)大。如圖3所示，渤海中部一直穩(wěn)定在一類(lèi)水質(zhì)，無(wú)機(jī)氮污染主要圍繞遼東灣，渤海灣和萊州灣變動(dòng)。活性磷酸鹽污染主要在遼東灣區(qū)域變動(dòng)，萊州灣基本穩(wěn)定在一類(lèi)水質(zhì)，渤海灣水質(zhì)等級(jí)和污染區(qū)域均變動(dòng)較大。渤海區(qū)域COD基本穩(wěn)定在一類(lèi)水質(zhì)，僅遼東灣上部、萊州灣西部等區(qū)域會(huì)出現(xiàn)變動(dòng)，在某些年份會(huì)出現(xiàn)污染。站位布設(shè)和監(jiān)測(cè)應(yīng)重點(diǎn)關(guān)注污染嚴(yán)重和污染變化劇烈區(qū)域。

2.2 站位時(shí)空重要性評(píng)價(jià)結(jié)果

根據(jù)渤海 2000年-2010年的監(jiān)測(cè)數(shù)據(jù)，選定連續(xù) 9年的監(jiān)測(cè)站位作為站位重要性的選擇標(biāo)準(zhǔn)，結(jié)果表明連續(xù)9年以上的站位為77個(gè)。再對(duì)77個(gè)站位中距離過(guò)近（小于0.1度）的站位進(jìn)一步篩選，最終確定59個(gè)站位具有歷史延續(xù)性。

圖2 渤海主要污染物濃度分布圖（mg/L）Fig. 2 Concentration distribution of main pollutants in the Bohai Sea（mg/L）

采用站位重要性評(píng)價(jià)方法評(píng)價(jià)92個(gè)站位在影響渤海環(huán)境質(zhì)量的主要指標(biāo)無(wú)機(jī)氮、活性磷酸鹽和COD中的重要性，最終得出92個(gè)站位對(duì)水質(zhì)等級(jí)變化影響的排序，根據(jù)這三項(xiàng)指標(biāo)評(píng)價(jià)的最大值來(lái)確定站位最終的空間重要性。選擇水質(zhì)類(lèi)別改變大于或等于 2%的站位在水質(zhì)空間評(píng)價(jià)中具有重要作用，評(píng)估結(jié)果表明歷史延續(xù)重要性評(píng)估確定的59個(gè)以外的92個(gè)站位中，有61個(gè)站位對(duì)渤海水質(zhì)空間評(píng)價(jià)具有重要作用（圖4）。2010年渤海151個(gè)監(jiān)測(cè)站位的評(píng)估結(jié)果如圖5所示，其中59個(gè)站位具有歷史延續(xù)性，61個(gè)站位對(duì)渤海水質(zhì)空間評(píng)價(jià)具有重要作用，其他站位為31個(gè)。

圖3 2000年-2010年渤海水質(zhì)穩(wěn)定性分布圖Fig. 3 Stability distribution of water environment in the Bohai Sea from 2000 to 2010

圖4 去除站位數(shù)量與評(píng)價(jià)結(jié)果改變的關(guān)系Fig. 4 Relationship between nember of stations removed and changes in evaluation results

圖5 2010年渤海海域監(jiān)測(cè)站位的評(píng)估結(jié)果Fig. 5 Evaluation results from monitoring stations in the Bohai Sea in 2010

2.3 渤海站位重新優(yōu)化布設(shè)結(jié)果

2010年渤海趨勢(shì)性監(jiān)測(cè)站位優(yōu)化布設(shè)的結(jié)果如圖6所示，保持原有151個(gè)站位總數(shù)不變，調(diào)整的31個(gè)站位在污染變化劇烈和較大的A區(qū)和B區(qū)，數(shù)量分配分別為9個(gè)和22個(gè)，主要分布在遼東灣、渤海灣和萊州灣及其周邊海域，這些區(qū)域污染面積和污染區(qū)域一直處于不穩(wěn)定狀態(tài)，是監(jiān)測(cè)需要重點(diǎn)關(guān)注的區(qū)域，清潔穩(wěn)定的渤海中部沒(méi)有增設(shè)站位。

2.4 站位優(yōu)化的效益分析

站位優(yōu)化的效益主要體現(xiàn)在經(jīng)濟(jì)效益和監(jiān)測(cè)的精確度上。站位監(jiān)測(cè)的經(jīng)濟(jì)投入主要與出海天數(shù)及監(jiān)測(cè)站位數(shù)量有關(guān)，站位優(yōu)化前后監(jiān)測(cè)海域的范圍及站位數(shù)量均無(wú)變化，所以站位監(jiān)測(cè)的經(jīng)濟(jì)投入不變。監(jiān)測(cè)效果體現(xiàn)在水質(zhì)評(píng)價(jià)的結(jié)果上，水質(zhì)評(píng)價(jià)采用的是 IDW 插值方法，插值誤差與水質(zhì)變化率及站位布設(shè)疏密程度相關(guān)，水質(zhì)變化率越低，站位布設(shè)越密，則插值誤差越小。站位優(yōu)化后，水質(zhì)污染嚴(yán)重和污染變化劇烈的區(qū)域整體站位增多且布設(shè)更加均勻，有效的避免了該區(qū)域出現(xiàn)較大空白而使插值誤差過(guò)大，從而使得監(jiān)測(cè)結(jié)果的精確度更佳。

圖6 渤海站位優(yōu)化結(jié)果Fig. 6 Optimization results from Bohai Sea monitoring station

3 結(jié)論

以2000年-2012年夏季無(wú)機(jī)氮、活性磷酸鹽和化學(xué)需氧量這三個(gè)主要污染因子趨勢(shì)性監(jiān)測(cè)數(shù)據(jù)為基礎(chǔ)，采用水質(zhì)空間穩(wěn)定性評(píng)價(jià)、監(jiān)測(cè)站位時(shí)空重要性評(píng)價(jià)以及站位優(yōu)化布設(shè)等方法，對(duì)渤海2010年趨勢(shì)性水質(zhì)監(jiān)測(cè)站位進(jìn)行優(yōu)化布設(shè)。在保持原有 151個(gè)站位總數(shù)不變的情況下，將距離過(guò)近或者位于清潔穩(wěn)定區(qū)的31個(gè)站位重新布設(shè)在污染變化較大的區(qū)域，重新布設(shè)的站位主要分布在遼東灣、渤海灣和萊州灣及其周邊海域，這些區(qū)域污染面積和污染區(qū)域一直處于不穩(wěn)定狀態(tài)。在監(jiān)測(cè)經(jīng)費(fèi)投入不變的情況下，這些區(qū)域增加站位能夠?yàn)楦鼫?zhǔn)確、更客觀地監(jiān)測(cè)海洋環(huán)境質(zhì)量提供保障和基礎(chǔ)，為渤海海域水質(zhì)的時(shí)空變化趨勢(shì)評(píng)價(jià)，提供更為有效的支撐。

參考文獻(xiàn)：

［1］ Jonathan A G，Carlos R，Erin L H，et al. Least cost distance analysis for spatial interpolation［J］. Computers & Geosciences，2011，37（2）： 272-276.

［2］ Wolfram Rühaak. A Java application for quality weighted 3-d interpolation［J］. Computers& Geosciences，2006，32（1）： 43-51.

［3］ Raymond D W，Niall G B，Christopher B J. The use of medium point density LIDAR elevation data to determine plant community types in Baltic coastal wetlands［J］. Ecological Indicators，2013，33： 96-104.

［4］ Vincent C，F(xiàn)rédéric D，Hocine B. Accuracy of interpolation techniques for the derivation of digital elevation models in relation to landform types and data density［J］. Geomorphology，2006，77（1-2）： 126-141.

［5］ Teegavarapu R S V，Chandramouli V. Improved weighting methods，deterministic and stochastic data-driven models for estimation of missing precipitation records［J］. Journal of Hydrology，2005，312： 191-206.

［6］ Michael K，Ida W，Allan R et al. Estimating areal rainfall over Lake Victoria and its basin using groundbased and satellite data［J］. Journal of Hydrology，2012，464-465： 401-411.

［7］ Louis de Mesnard. Pollution models and inverse distance weighting： Some critical remarks［J］. Computers & Geosciences，2013，52： 459-469.

［8］ Richard R，Mark B，Gregory M. Integrating forest biomass and distance from channel to develop an indicator of riparian condition［J］. Ecological Indicators，2012，23： 46-55.

［9］ Yunfeng Xie，Tong-bin Chen，Mei Lei，et al. Spatial distribution of soil heavy metal pollution estimated by different interpolation methods： Accuracy and uncertainty analysis［J］. Chemosphere，2011，82（3）： 468-476.

［10］ Dimitra K，Michael K. Categorical mapping of marine eutrophication based on ecological indices［J］. Science of The Total Environment，2000，255： 113-127.

［11］ Dilip K J，Anup K D，Saravanane N，et al. Sensitivity of GIS-based interpolation techniques in assessing water quality parameters of Port Blair bay，Andaman［J］. Journal of the Marine Biological Association of India，2010，52（1）： 55-61.

［12］ 蘭冬東，李冕，梁斌等.渤海水環(huán)境無(wú)機(jī)氮污染風(fēng)險(xiǎn)研究［J］.海洋科學(xué)，2013，37（11）59-64. Lan Dongdong，Li Mian，Liang Bin，et al. Study on the pollution risk of inorganic nitrogen in water environment of Bohai Sea［J］. Marine Sciences，2013，37（11）：59-64.

［13］ Chun-Hsien Wu，Kuo-Chuan Lee，Yeh-Ching Chung. A Delaunay Triangulation based method for wireless sensor network deployment［J］. Computer Communications，2007，30（14-15）： 2744-2752.

［14］ Weatherill N P. Delaunay triangulation in computational fluid dynamics［J］. Computers & Mathematics with Applications，1992，24（5-6）： 129-150.

［15］ Lertrattanapanich S，Bose N K. High resolution image formation from low resolution frames using Delaunay triangulation［J］. IEEE Trans Image Process，2002，11（12）： 1427-41.

［16］ Gopi M，Krishnan S，Silva C T. Surface Reconstruction based on Lower Dimensional Localized Delaunay Triangulation［J］. Computer Graphics Forum，2000，19（3）：467-478.

［17］ Sloan S W. A fast algorithm for constructing Delaunay triangulations in the plane［J］. Advances in engineering software，1987，9（1）： 34-55.

（本文編輯： 康亦兼）

中圖分類(lèi)號(hào)：X55

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1000-3096（2016）04-0088-06

doi：10.11759/hykx20150608001

收稿日期：2015-06-08； 修回日期： 2015-10-25

基金項(xiàng)目：國(guó)家自然科學(xué)基金青年科學(xué)基金（41306098）； 海洋公益性行業(yè)科研專(zhuān)項(xiàng)項(xiàng)目（201305023）； 國(guó)家海洋局近岸海域生態(tài)環(huán)境重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金（201312）

［Foundation item：The National Natural Science Foundation of China （No.41306098）； Public science and technology research funds projects of ocean （No.201305023）； The Key Laboratory for Ecological Environment in Coastal Areas，State Oceanic Administration（No.201312）］

作者簡(jiǎn)介：蘭冬東（1984-），女，遼寧大連人，碩士，助理研究員，主要從事海洋環(huán)境評(píng)價(jià)研究，E-mail： ddlan@nmemc.org.cn； 李冕，通信作者，助理研究員，主要從事海洋環(huán)境評(píng)價(jià)研究E-mail： mli@nmemc.org.cn級(jí)相同，則取概率值較大的。水質(zhì)發(fā)生概率大于75%為穩(wěn)定區(qū)域，50%～75%為變化區(qū)域，0～50%為變化劇烈區(qū)域。

Layout optimization of water quality monitoring stations in the Bohai Sea

LAN Dong-dong，LI Mian，XU Yan，BAO Chen-guang，YU Chun-yan，LIANG Bin，MA Ming-hui
（National Marine Environmental Monitoring Center，Dalian 116023，China）

Received： Jun. 8，2015

Key words：water quality monitoring； station assessment； station optimization； Bohai Sea

Abstract：The trend of water quality monitoring stations in the Bohai Sea was optimized and adjusted using a stability assessment method in relation to the water quality. An assessment of the spatiotemporal importance of monitoring stations and an optimal design method for the stations was determined on the basis of the trend of monitoring data between 2000 and 2010 for three main pollution factors，namely，inorganic nitrogen，active phosphate，and COD. Optimization of the layout of water quality monitoring stations could provide effective and objective support for the assessment of temporal trends in the water quality of the Bohai Sea.