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

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

?

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

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

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

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

關鍵詞：海水質量監(jiān)測； 站位評價； 站位優(yōu)化； 渤海

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

1 研究方法

1.1 IDW空間插值

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

式中，Z*（B）為待插值點的值，Z（xi）為插值結點的實測值，λi為插值結點xi所占的權重，其中λi的計算方法為：

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

1.2 水質穩(wěn)定性評價

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

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

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

1.3 站位時空重要性評價方法

1.3.1 歷史重要性站位篩選

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

1.3.2 站位貢獻率評價

綜合水質站位貢獻率為站位各分指標貢獻率的和。設X1是站位集X的一個子集，則X1在X中的貢獻率計算公式為：

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

1.3.3 站位綜合評價

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

（1） n個監(jiān)測站位的數(shù)據(jù)水質評價結果作為原始結果，記錄各類水質的柵格數(shù)量及空間位置，站位貢獻率計算均以此為基礎。

（2） 對n-m中的每一站位，計算其貢獻率，按升序排序，去除貢獻率最小的站位，記錄各站位對應的貢獻率。

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

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

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

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

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

圖1 四頂點加權圖Fig. 1 Four-vertex weighted graph

2 結果與討論

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

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

2.2 站位時空重要性評價結果

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

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

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

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

圖4 去除站位數(shù)量與評價結果改變的關系Fig. 4 Relationship between nember of stations removed and changes in evaluation results

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

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

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

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

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

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

3 結論

以2000年-2012年夏季無機氮、活性磷酸鹽和化學需氧量這三個主要污染因子趨勢性監(jiān)測數(shù)據(jù)為基礎，采用水質空間穩(wěn)定性評價、監(jiān)測站位時空重要性評價以及站位優(yōu)化布設等方法，對渤海2010年趨勢性水質監(jiān)測站位進行優(yōu)化布設。在保持原有 151個站位總數(shù)不變的情況下，將距離過近或者位于清潔穩(wěn)定區(qū)的31個站位重新布設在污染變化較大的區(qū)域，重新布設的站位主要分布在遼東灣、渤海灣和萊州灣及其周邊海域，這些區(qū)域污染面積和污染區(qū)域一直處于不穩(wěn)定狀態(tài)。在監(jiān)測經(jīng)費投入不變的情況下，這些區(qū)域增加站位能夠為更準確、更客觀地監(jiān)測海洋環(huá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)境無機氮污染風險研究［J］.海洋科學，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.

（本文編輯： 康亦兼）

中圖分類號：X55

文獻標識碼：A

文章編號：1000-3096（2016）04-0088-06

doi：10.11759/hykx20150608001

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

基金項目：國家自然科學基金青年科學基金（41306098）； 海洋公益性行業(yè)科研專項項目（201305023）； 國家海洋局近岸海域生態(tài)環(huán)境重點實驗室開放基金（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）］

作者簡介：蘭冬東（1984-），女，遼寧大連人，碩士，助理研究員，主要從事海洋環(huán)境評價研究，E-mail： ddlan@nmemc.org.cn； 李冕，通信作者，助理研究員，主要從事海洋環(huán)境評價研究E-mail： mli@nmemc.org.cn級相同，則取概率值較大的。水質發(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.