基于時間序列的三沙灣水質狀況變化趨勢分析

王臻

(福建省海洋環境與漁業資源監測中心，福建福州 350003)

?

基于時間序列的三沙灣水質狀況變化趨勢分析

王臻

(福建省海洋環境與漁業資源監測中心，福建福州 350003)

根據2007—2016年福建省寧德市三沙灣海水水質狀況調查資料，采用時間序列模型對三沙灣海域水體的富營養化和劣四類海水水質面積比例變化趨勢進行分析。同時，利用資源與環境承載力評價模型對三沙灣海域空間開發強度進行分析。結果表明，三沙灣海域開發強度指數為0.33，開發強度較高。2017年海水富營養化指數預測為4.5，變化不大，依然處于中度富營養化狀態。劣四類水質面積比例為36.77%，較2016年有所惡化。三沙灣水質的惡化可能與其開發強度較高有關。

三沙灣；水質變化趨勢；時間序列

三沙灣位于福建省東北部沿海，是我國天然良港，漁業資源十分豐富。海灣由一澳(三都澳)、三港(盧門港、白馬港、鹽田港)、三洋(東吾洋、官井洋、福鼎洋)等次一級海灣匯集而成[1]，是福建省重要的港灣之一。作為許多專家學者的研究對象，對三沙灣水質狀況的研究報告也頗為眾多[2-6]。但是對水質狀況變化趨勢的預測較少，行政管理部門在對三沙灣的開發與保護的決策上，缺少科學依據。本文利用2007—2016年三沙灣的調查結果，從水體污染總體狀況出發，利用SPSS軟件中的時間序列模型預測三沙灣水質變化趨勢，為相關部門做出決策提供科學支撐。

1 材料與方法

1.1 水樣采集與監測方法

水樣的采集依據《海洋監測規范 第3部分樣品采集、貯存和運輸》(GB 17378.3—2007)執行。實驗室分析方法分按照《海洋監測規范 第四部分：海水分析》(GB 17378.4—2007)執行，其中COD采用堿性高錳酸鉀法，石油類采用紫外分光光度法，PO4-P采用磷鉬藍分光光度法，NO2-N采用重氮-偶氮分光光度法，NO3-N采用鋅-鎘還原法，NH4-N采用次溴酸鈉氧化法進行分析。

1.2 評價依據和方法

1.2.1 評價依據

水質評價按照《海水水質標準》(GB 3097—1997)。

1.2.2 空間插值法

劣四類水質面積采用距離反比例法進行空間插值。

1.2.3 富營養化評價法

海水中富營養化狀況評價采用富營養化指數法，見公式：

E=(CCOD×CDIN×CDIP×106)/4500

式中，E為富營養化指數；CCOD為化學需氧量濃度，mg/L；CDIN為無機氮濃度，即亞硝酸鹽氮(NO2--N)、硝酸鹽氮(NO3--N)、氨氮(NH4+-N)的總和，mg/L；CDIP為活性磷酸鹽濃度，mg/L。

1.3 預測模型

本文主要采用SPSS軟件中的B-J模型來預測水質趨勢，基本模型有3種：自回歸(AR)模型、移動平均(MA)模型和自回歸綜合移動(ARIMA)模型。

1.3.1 自回歸模型

一般形式為：xt=Φ1xt-1+Φ2xt-2+…+Φpxt-p+εt

式中，εt假設為白噪聲序列，且和t時刻之前的原始序列xk(k

1.3.2 移動平均模型

一般形式為：xt=εt+θ1εt-1+θ2εt-2+…+θqεt-q

式中，εt為白噪聲序列，說明時間序列xt能表示是為若干個白噪聲的加權平均和。此式稱為q階移動平均模型，記為MA(q)。

1.3.3 自回歸移動平均模型

一般形式為：xt=Φ1xt-1+Φ2xt-2+…+Φpxt-p+εt+θ1εt-1+θ2εt-2+…+θqεt-q

式中，εt為白噪聲序列，且和t時刻之前的原始序列xk(k

2 結果與討論

2.1 統計結果

從表1中可以看出，十年來，活性磷酸鹽、無機氮、石油類含量和COD基本變化不大。2013年后富營養化程度有所減輕，劣四類水質面積逐漸減少。由于2013年劣四類水質比例面積變化太大，對預測結果可能產生較大影響，因此在預測時應作特殊處理，如圖1所示。

圖1 2007—2016年水質變化趨勢Fig.1 The variation trend of water quality from 2007 to 2016

圖2 富營養化指數自相關函數和偏相關函數值Fig.2 The ACF and PACF of eutrophication index

年份活性磷酸鹽/(mg/L)無機氮/(mg/L)石油類/(mg/L)COD/(mg/L)富營養化劣四類水質面積比例/%20070.0320.3760.0320.902.439.120080.0230.2750.0140.550.845.620090.0570.4280.0140.573.151.420100.0480.4120.0130.924.045.120110.0430.3460.0120.872.934.420120.0350.3750.0120.992.931.420130.0560.3800.0141.245.999.920140.0470.5110.0160.864.658.820150.0390.4860.0170.803.434.820160.0370.4540.0180.823.014.4

2.2 預測與分析

2.2.1 富營養化指數

從圖2可以看出富營養化指數的ACF(自相關函數值)和PACF(偏自相關函數值)完全落入上下置信區間，判斷序列為平穩序列。通過表2看出，自相關函數拖尾，偏自相關函數5階截尾，因此采用AR(5)模型進行預測，得到2017年富營養化指數預測值為4.5，處于中度富營養化狀態。

表2 富營養化指數偏相關和偏自相關函數值

圖3可以看出AR(5)在2009年以后觀測值與預測值可以較好地擬合。由表3看出，5階滯后t統計量顯著性(Sig.)小于0.05，判斷該模型較為恰當[7]。根據圖4殘差自相關和偏相關函數值均在置信區間以內，表明殘差統計值滿足白噪聲的兩個基本標準[8-9]，由此說明AR(5)可以較好地模擬富營養化指數變化趨勢，2017年預測值基本可信。

圖3 富營養化指數擬合圖Fig.3 The fitting chart of eutrophication index

估計SEtSig.常數3.3410.22414.9230.000AR滯后10.3910.1402.8030.049滯后2-0.0900.202-0.4440.680滯后3-0.0510.176-0.2920.785滯后40.3710.1852.0060.115滯后5-0.9460.151-6.2490.003

圖4 殘差自相關和偏相關函數值Fig.4 The ACF and PACF of residual

2.2.2 劣四類水質面積比例

在剔除了2013年的異常數據后，在SPSS軟件中采用線性插值法替換缺失值，替換值為45.1%，并對新的數據組進行自相關分析。

表4 劣四類水質面積比例偏相關和偏自相關函數值

圖5 劣四類水質面積比例自相關函數和偏相關函數Fig.5 The ACF and PACF of the water area worse than Grade IV

從圖5可以看出，劣四類水質面積比例的ACF(自相關函數值)和PACF(偏自相關函數值)基本落入上下置信區間，判斷序列為平穩序列。通過表4看出，自相關函數拖尾，偏自相關函數2階截尾，因此采用AR(2)模型進行預測，得到2017年劣四類水質面積預測值為36.77%，見圖6。

由圖7看出AR(5)觀測值與預測值可以較好地擬合。2階滯后t統計量顯著性(Sig.)小于0.05，判斷該模型較為恰當(見表5)。由圖7看出，AR(2)模型殘差自相關和偏自相關函數值均在置信區間以內，預測值基本可信。

圖6 劣四類水質面積比例擬合圖Fig.6 The fitting chart of the proportion of water area worse than Grade IV

估計SEtSig.常數41.0701.64125.0280.000AR滯后10.3740.2631.4200.199滯后2-0.9040.324-2.7880.027

圖7 殘差自相關和偏自相關函數值Fig.7 The ACF and PACF of residual

2.3 海域開發強度評價

鑒于三沙灣的富營養化指數和劣四類水質面積比例較2016年均有所提高可能是由于海域開發造成的，因此以2015年的海域調查資料為基礎，對三沙灣海域開發強度進行評價。海洋資源環境承載力的模型眾多[10-12]，本文采用國家海洋局海洋領域研究組編制的《海洋資源環境承載能力監測預警技術方法

(2016)》進行評價。

2.3.1 海域開發資源效應指數

三沙灣海域使用類型有漁業基礎設施用海、圍海養殖用海、開放式養殖用海、港口用海、路橋用海、臨海工業用海、旅游基礎設施用海、城鎮建設用海和海岸防護工程用海,如圖8所示

海域開發資源效應指數(PE)采用如下計算公式：

式中，PE為海域開發資源效應指數，n為海域使用類型數，Si為第i種用海類型的面積，S為省級海洋功能區劃數據獲得的海域總面積，li為第i種用海類型的資源耗用系數。通過遙感影像及現場補充調查，統計得到各類型用海面積，見表6。根據表6，計算得三沙灣海域開發資源效應指數PE=0.15。

圖8 三沙灣海域開發利用現狀Fig.8 The current status of development and utilization of the sea area in Sansha Bay

2.3.2 海域空間開發利用標準

以《福建省海洋功能區劃(2011—2020年)》為基礎得出三沙灣海域各類功能區面積，并計算該海域空間開發利用標準。海域空間開發利用標準計算公式如下：

式中，PM0為海洋空間開發利用標準;ai為第i類海洋功能區面積;hi為第i類海洋功能區允許的海洋開發程度，賦值方法見表7;S為區域海洋功能區劃面積總和。經計算，三沙灣海域海洋空間開發利用標準PM0=0.46。

2.3.3 三沙灣海域開發強度評估

根據海域開發資源效應指數與海洋空間開發利用標準之比，得到三沙灣海域開發強度S2，其計算公式如下：

經計算，三沙灣海域開發資源效應指數PE=0.15。海洋空間開發利用標準PM0=0.46，因此，海域開發強度S2=0.33。

表6 海域使用類型資源耗用指數

由于當S2>0.3時，海域開發強度較高；當S2介于0.15～0.3時，海域開發強度臨界；當S2≤0.15時，海域開發強度適宜。因此，三沙灣海域開發強度較高。

表7 主要海洋功能區海洋開發對海域空間資源的影響

3 結論

(1)時間序列模型預測水質變化，能在一定程度上反映排除人為因素干擾下，水體變化情況。預測過程需要排除異常值的影響。異常值不能直接剔除，需要通過替換缺失值的方式進行補齊，否則會導致原有時間序列周期性的破壞，而無法得到正確的預測分析結果。

(2)三沙灣的主要污染物為無機氮和活性磷酸鹽。近年來，各參數變化幅度不大。2013年的劣四類水質面積上升幅度較大，之后幾年回落到正常水平。2013年水質狀況較為異常，需要結合當年多方面因素對數據作進一步分析判斷。

(3)2017年富營養化指數預測值為4.5，較2016年有所升高，但變化幅度不大，屬于基本維持現狀，依然處于中度富營養化狀態。劣四類水質面積比例預測值為36.77%，〗較2016年來說，水質有所惡化。三沙灣海域開發強度為0.33，開發強度較高，水質惡化原因可能與三沙灣海域開發強度較高有關。相關行政主管部門應對此高度重視。

[1] 劉修德. 福建省海灣數模與環境研究——三沙灣[M]. 北京: 海洋出版社, 2009, 8- 9.

[2] 劉家富, 鄭欽華, 陳洪清, 等. 三沙灣的水質狀況[J]. 臺灣海峽, 2003, 22(2): 201- 204.

[3] 蔡清海. 福建三沙灣海洋生態環境研究[J]. 中國環境監測, 2007, 23(6): 101- 105.

[4] 鄭欽華, 王曉清, 何湘蓉, 等. 三沙灣海水中無機氮的變化特征[J]. 河北漁業, 2008(12): 43- 46.

[5] 馬祖友, 夏永健, 石志洲, 等. 2011年三沙灣增養殖區水環境質量評價[J]. 海洋開發與管理, 2013(7): 75- 78.

[6] 邵留, 于克鋒, 吳海龍, 等. 三沙灣海域水質周年變化分析與評價[J]. 上海海洋大學學報, 2014, 23(2): 228- 236.

[7] 杜強, 賈麗艷. SPSS統計分析從入門到精通[M]. 北京: 人民郵電出版社, 2013: 131.

[8] 黃杰, 錢瑾, 殷浩文, 等. ARIMA模型在蘇州河總磷預測中的應用[J]. 環境工程學報, 2007, 1(1): 139- 143.

[9] 楊位欽, 顧嵐. 時間序列分析與動態數據建模[M]. 北京: 北京工業學院出版社, 1986: 361.

[10] 劉佳, 萬榮, 陳曉文. 山東省藍色經濟區海洋資源承載力測評[J]. 海洋環境科學, 2013, 32(4): 619- 624.

[11] 劉蕊. 海洋資源承載力指標體系的設計與評價[J]. 廣東海洋大學學報, 2009, 29(5): 6- 9.

[12] 鄧宗成, 孫英蘭, 周皓, 等. 沿海地區海洋生態環境承載力定量化研究[J]. 海洋環境科學, 2009, 28(4): 438- 441.

Analysis of Variation Trend of Water Quality Based on Time Series in Sansha Bay

WANG Zhen

(Ocean Environment and Fishery Resources Monitoring Center of Fujian Province, Fuzhou 350003, China)

According to the survey data of water quality in Sansha Bay, Ningde City, Fujian Province from 2007 to 2016, the time series model is used to analyze the variation trend of the eutrophication and proportion of the water area which was worse than Grade IV. At the same time, the paper analyzes the spatial development intensity of the area in Sansha Bay based on the resource and environment bearing capacity model. The results showed that the development intensity index was 0.33 in Sansha Bay, which indicated high development intensity. In 2017, the eutrophication index was forecasted to be 4.5, still a moderate eutrophication. The proportion of the water area worse than Grade IV was 36.77%, which was worse than that in 2016. The deterioration of water quality in Sansha Bay may be related to its high development intensity.

Sansha Bay; variation trend of water quality; time series

2017-06-08

福建省海洋與漁業廳資源與環境承載力監測專項(閩海漁〔2015〕289號)

王臻(1982—)，男，福建福州人，工程師，本科，主要從事海洋環境監測和評價工作，E-mail：wangcom20@126.com

10.14068/j.ceia.2017.04.017

X824

A

2095-6444(2017)04-0076-06