王 竹，劉 揚，王 芳，朱士江

(1. 流域水循環(huán)模擬與調(diào)控國家重點實驗室，北京 100038；2. 中國水利水電科學研究院 水資源所，北京 100038； 3. 三峽大學 水利與環(huán)境學院，湖北宜昌443002)

水作為人類賴以生存和發(fā)展的基本物質(zhì)條件，維系著生態(tài)系統(tǒng)自然功能的同時，也支撐著社會經(jīng)濟系統(tǒng)的發(fā)展。隨著經(jīng)濟社會的快速發(fā)展，水環(huán)境污染造成水生態(tài)系統(tǒng)惡化，進而引發(fā)用水緊張、用水安全等問題。因此評價水環(huán)境質(zhì)量、分析水體主要污染物的時空變化特征，識別造成水環(huán)境現(xiàn)狀的主要原因，對保護當?shù)氐乃h(huán)境質(zhì)量、維護良好的生態(tài)系統(tǒng)具有至關重要的作用[1,2]。

目前，眾多學者普遍認為水質(zhì)污染的影響因素主要分為自然因素和社會因素，其中自然因素有降水、流量等[3-5]，社會因素主要有人口規(guī)模、區(qū)域經(jīng)濟發(fā)展、土地利用結(jié)構等[6-8]。王建華等[9]在全球水資源面臨人口增長、土地利用變化、氣候變化等壓力下，提出了水量水質(zhì)聯(lián)合模擬與評價模型，通過模擬計算灤河流域2006年5項水質(zhì)指標，發(fā)現(xiàn)灤河干流河道水體質(zhì)量以Ⅲ類為主，劣Ⅴ類次之。夏星輝等[10,11]從水質(zhì)水量相結(jié)合的角度評價黃河的水資源狀況，認為大約47%的水資源短缺是由水體污染所致，因此對沿岸排污進行合理的時空優(yōu)化可以幫助緩解黃河水資源短缺的局面。石代軍[12]研究了氣溫、降水、光照、風速等氣象因子對地表水環(huán)境質(zhì)量的變化，發(fā)現(xiàn)人類活動引起的氣候變化會加重氣象因子對水環(huán)境質(zhì)量的影響。黃奕龍等[13]通過觀測深圳市近20年的水質(zhì)數(shù)據(jù)，發(fā)現(xiàn)工業(yè)化、城市化和土地利用格局變化是水質(zhì)退化的主要驅(qū)動機制。張瑞等[14]綜合土地利用指標分析了雅魯藏布江35年間土地利用/覆被時空變化特征，發(fā)現(xiàn)氣候和水文要素是該流域土地利用的主要驅(qū)動因子，流域內(nèi)經(jīng)濟的發(fā)展及人口增長是城鄉(xiāng)建設用地和耕地發(fā)生變化的驅(qū)動力。經(jīng)濟發(fā)展往往伴隨著水環(huán)境質(zhì)量的惡化，因此在社會經(jīng)濟發(fā)展的同時做好污染源頭控制和清潔生產(chǎn)顯得尤為重要。

近些年來，環(huán)渤海區(qū)域經(jīng)濟快速發(fā)展、人口日益密集，各種污染源排入河流和海洋，造成該區(qū)域水質(zhì)惡化，環(huán)渤海區(qū)域已經(jīng)成為我國三大經(jīng)濟帶中水體污染最為嚴重的區(qū)域[15]。作為中國首批沿海開放城市和世界最大的能源輸出港，秦皇島市工農(nóng)業(yè)發(fā)展迅速，人口日益增多，直接導致了秦皇島市及其他環(huán)渤海區(qū)港口城市陸域入海河流和渤海海域水生態(tài)環(huán)境狀況的惡化。因此，本文選擇單因子評價與主成分分析相結(jié)合的方法評價研究區(qū)水質(zhì)現(xiàn)狀，旨在識別不同季節(jié)、不同時期的水質(zhì)時空變化特征，分析土地利用與水質(zhì)污染的相關關系，為評價秦皇島市流域水質(zhì)狀況、識別流域污染主要原因提供參考，并為進一步制定該市水資源保護和管理措施提供科學依據(jù)。

1 研究區(qū)概況

秦皇島市地處河北省和華北平原東北部，北倚燕山，南襟渤海，北緯39°24′～40°37′，東經(jīng)118°33′～119°51′，位于環(huán)渤海經(jīng)濟圈的中心地帶，總面積7 812.4 km2。秦皇島市內(nèi)的河流分屬灤河水系和冀東沿海諸河，其中灤河在境內(nèi)的流域面積為3 773.7 km2，沿海諸河自西向東有石河、湯河、戴河、洋河、飲馬河等。隨著周邊地區(qū)近些年來工農(nóng)業(yè)的快速發(fā)展，秦皇島市水污染日益加重，生態(tài)環(huán)境受到了一定程度的破壞，也進一步威脅著渤海海洋生態(tài)系統(tǒng)。

2 材料與方法

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

圖1 秦皇島市水系及監(jiān)測斷面分布Fig.1 Water system and monitoring section distribution in Qinhuangdao

(2)土地利用數(shù)據(jù)來源。本研究選取2014年TM遙感影像(空間分辨率30 m)，采用ENVI 4.8軟件，進行遙感解譯和影像配準，采用最大似然分類與交互式目視解譯相結(jié)合的方法，得到2014年土地利用柵格數(shù)據(jù)。

(3)子流域劃分。研究選取地理空間數(shù)據(jù)云(www.gscloud.cn)90 m的DEM為原始數(shù)據(jù)，利用ArcGIS 10.2中的Hydrology工具提取水系，并運用水文分析中的Basin工具將秦皇島市劃分為7個子流域(圖2)。

2.2 研究方法

2.2.1 水質(zhì)評價方法

(1)單因子評價法。單因子評價法是通過比較評價因子的實測濃度與評價標準值的大小來判斷水質(zhì)級別的方法。

(1)

表1 監(jiān)測斷面地理信息Tab.1 Geographical information of monitoring sections

圖2 秦皇島市子流域劃分Fig.2 Sub-basin division of Qinhuangdao

式中：Pi是某項指標的相對污染值；ci是某指標的最大實測濃度(溶解氧取最小值)；c0是某指標的評價標準值。當Pi>1時，表示水質(zhì)超標，反之則達標。

(2)主成分分析法。主成分分析法可通過建立水質(zhì)綜合評價指標，對各斷面水質(zhì)進行排序，并分析各項指標在綜合評價中的作用，識別主要污染因子，分析水質(zhì)污染原因。主要步驟為：

①建立原始變量矩陣X，由m個樣本的n個因子組成：

(2)

②對原始變量矩陣X進行標準化處理，采用Z-Score進行標準化：

(3)

(4)

(5)

③計算標準化數(shù)據(jù)的相關系數(shù)矩陣及其特征根，確定主成分個數(shù)P；

④確定主成分Fi(i=1,2,…,p)的表達式如下：

(6)

⑤確定綜合評價函數(shù)F：

(7)

式中：a1m,a2m,…,anm為原始變量矩陣X的協(xié)方差矩陣∑的特征值對應的特征向量；ZX1,ZX2,…,ZXn為原始變量矩陣X經(jīng)過Z-Score標準化處理后的值；λ1，λ2，…，λp為矩陣ZX的特征值；m為樣本個數(shù)；n為評價因子的個數(shù)；p為主成分個數(shù)。監(jiān)測斷面水質(zhì)的時間變化趨勢采用Excel和SPSS 22軟件進行分析。

2.2.2 土地利用與水質(zhì)相關分析

根據(jù)ArcGIS 10.2提取各點所在子流域的土地利用數(shù)據(jù)，采用Spearman相關分析法在SPSS 22軟件中建立水質(zhì)參數(shù)與土地利用類型面積比的關系，以定性評價土地利用類型對水質(zhì)的影響。

3 結(jié)果與分析

3.1 水質(zhì)指標時間變化特征

將2017年1-12月分為春季(3-5月)、夏季(6-8月)、秋季(9-11月)、冬季(1-2月、12月)、汛期(6-9月)和非汛期(1-5月、10-12月)的不同季節(jié)和時期進行水質(zhì)指標評價，具體統(tǒng)計值見表2，各項評價指標的變化趨勢見圖3～圖6。

3.1.1 季節(jié)變化特征

結(jié)合表2及圖3～6可知，四項指標的濃度在不同季節(jié)差異較大，且變化趨勢各不相同。具體來看，COD濃度：春季>冬季>夏季>秋季，其中秋季最優(yōu)，為Ⅱ類水標準，其他3個季節(jié)為Ⅲ類水，春季最差。這主要是因為隨著溫度的升高，春季積雪開始融化，水流攜帶大量混合著有機物的泥沙進入河道，導致COD濃度升高[17]。從全年來看，飲馬河-楊古泊村北斷面9月COD濃度達71 mg/L，超出水功能區(qū)水質(zhì)目標1.37倍，為研究區(qū)全年COD濃度峰值。

高錳酸鹽指數(shù)濃度：春季>秋季>夏季>冬季，在夏秋冬三季的濃度均較低，為Ⅱ類水標準。可能原因是春耕導致了農(nóng)業(yè)面源污染的增加，或存在部分工業(yè)企業(yè)偷排等現(xiàn)象[19-21]。高錳酸鹽指數(shù)年度峰值出現(xiàn)在9月的減河-364省道橋斷面，超水功能區(qū)水質(zhì)目標0.72倍。

表2 不同時期監(jiān)測斷面水質(zhì)指標結(jié)果統(tǒng)計Tab.2 Statistics of water quality indicators of monitoring sections in different periods

3.1.2 汛期/非汛期變化特征

圖3 監(jiān)測斷面污染濃度變化(COD)Fig.3 Concentration change of pollutants in monitoring sections (COD)

圖4 監(jiān)測斷面污染濃度變化Fig.4 Concentration change of pollutants in monitoring sections

圖5 監(jiān)測斷面污染濃度變化(TP)Fig.5 Concentration change of pollutants in monitoring sections (TP)

圖6 監(jiān)測斷面污染濃度變化(高錳酸鹽指數(shù))Fig.6 Concentration change of pollutants in monitoring sections (permanganate index)

圖7 監(jiān)測斷面污染濃度變化(汛期)Fig.7 Concentration change of pollutants in monitoring sections (flood season)

圖8 監(jiān)測斷面污染濃度變化(非汛期)Fig.8 Concentration change of pollutants in monitoring sections (non-flood season)

3.2 水質(zhì)指標空間變化特征

表3 監(jiān)測斷面年均值評價結(jié)果Tab.3 Monthly evaluation results of monitoring sections

采用主成分分析法分析4項指標全年的均值，確定了一個主要包含與高錳酸鹽指數(shù)、TP信息有關的主成分，并根據(jù)成分矩陣計算所得綜合得分F值對各斷面進行排序，F(xiàn)值越小表示斷面水質(zhì)越好。由表3結(jié)果可知，沙河-冷口(S5)、青龍河大匯合(S4)、青龍河界嶺口(S8)等斷面水質(zhì)偏好，水質(zhì)較差的斷面主要有減河-364省道橋(S15)、飲馬河-楊古泊村北(S14)、飲馬河-昌黎(S13)等。

兩種評價方法的結(jié)果稍有差異，如主成分分析法認為沙河-冷口(S5)斷面水質(zhì)相對最好，其次為青龍河大匯合(S4)、洋河青龍界嶺口(S8)等，而單因子評價法結(jié)果顯示洋河水庫壩上(S7)、石河水庫壩上(S12)水質(zhì)最好，其次為南河-青龍鎮(zhèn)(S6)等斷面。由于主成分分析法相較單因子評價考慮了多項指標的貢獻率，評價結(jié)果的綜合性更強，但兩種評價方法得出污染嚴重的斷面較為一致，主要為飲馬河、減河上的S13、S14、S15斷面。整體來看，研究區(qū)內(nèi)上游水質(zhì)較下游好，北部山區(qū)水質(zhì)優(yōu)于南部沿海諸河。主要原因有兩個方面，其一是各河流沿程建設攔河壩等水利工程，導致中下游水量較少，致使河流水體的水環(huán)境容量明顯下降；其二是工業(yè)和生活排水大部分都集中在河流中下游地區(qū)，導致污水排放較為集中，超過河流水體的自凈能力。

3.3 土地利用與水質(zhì)的相關性分析

3.2.1 土地利用結(jié)構分析

將研究區(qū)2014年土地利用柵格數(shù)據(jù)按居民用地、未利用土地、耕地、草地、水域等(一級)土地利用類型進行合并，研究區(qū)2014年土地利用狀況見圖9。利用ArcGIS中的疊加工具統(tǒng)計各土地利用類型在所劃分的7個子流域中的比例，具體結(jié)果見圖10。

圖9 秦皇島2014年土地利用狀況Fig.9 Land use patterns of Qinhuangdao in 2014

圖10 各子流域土地利用類型比例Fig.10 Proportion of land use patterns in Sub-basins

根據(jù)各子流域中土地利用結(jié)構可知，秦皇島市內(nèi)最主要的土地利用類型為耕地，占比42.2%，其次為居民用地和林地，分別為22.2%和22.3%，草地、水域及未利用土地占比較少，分別為10.0%、3.3%、0.05%。各子流域的主要土地利用類型均為耕地，耕地占比：W2>W3>W5>W7>W6>W1>W4，其中W2比例最高達71.5%；草地占比：W1>W4>W3>W5>W6>W2>W7，其中W1比例最高達28.2%；居民用地占比：W6>W7>W5>W4>W2>W3>W1，其中W6比例最高達39.0%；林地占比：W4>W1>W3>W6>W7>W2>W5，其中W4比例最高達47.1%；水域占比：W7>W1>W3>W6>W2>W5>W4，其中W7比例最高達6.6%；未利用土地占比：W1>W4>W3，最高比例為0.27%，W2、W5、W6和W7內(nèi)均無未利用土地。秦皇島境內(nèi)的灤河上游以林地、草地為主，下游以耕地為主，冀東沿海諸河以耕地和居民用地為主，可見研究區(qū)內(nèi)河流受人類活動影響較大。

3.2.2 相關性分析

分析監(jiān)測斷面4項水質(zhì)指標的平均值與各子流域2014年土地利用類型面積比例的相關關系，結(jié)果見表4。由于本次數(shù)據(jù)來源中，在W6子流域中未設置監(jiān)測斷面，因此該子流域不參與相關性分析。

表4 土地流域類型面積比例與水質(zhì)指標的相關性分析Tab.4 Analysis of the correlation between land basin type area ratio and water quality index

注：*表示顯著性水平為0.05，**表示顯著性水平為0.01。

以上分析可知，污染嚴重的河流如減河、飲馬河位于以耕地為主的W2、W5子流域，湯河位于以居民用地為主的W4子流域，而水質(zhì)較好的斷面如灤河上游、洋河水庫等位于以林草地為主的W1子流域，說明人類活動對水質(zhì)影響程度巨大。耕地、居民用地等人類用地代表著人類工農(nóng)業(yè)活動強度的增加，產(chǎn)生各種污水排入水體中，成為水環(huán)境污染的“源頭”；林地、草地等自然用地往往代表著豐富植被、微生物的存在，此類區(qū)域生態(tài)系統(tǒng)的生物地球化學過程效率往往較高[26]，能夠更有效地降解有機物、吸收轉(zhuǎn)化氮磷等營養(yǎng)物質(zhì)，以上結(jié)果與相關研究[27,28]結(jié)論一致。因此，秦皇島市內(nèi)林地、草地、未利用土地的管理和保護，以及對居民用地、耕地的控制是減少水質(zhì)污染的防治關鍵。

4 結(jié) 論

通過分析研究區(qū)的水質(zhì)時空變化特征發(fā)現(xiàn)：

(3)研究區(qū)耕地、居民用地等人類用地面積與污染指標濃度呈正相關，林地、草地、水域及未利用土地等自然用地與水質(zhì)指標呈負相關。表明農(nóng)業(yè)種植、養(yǎng)殖、城鄉(xiāng)生活等人類活動是該區(qū)域水質(zhì)污染的主要來源。因此，科學耕種、節(jié)水減肥、發(fā)展集約化農(nóng)業(yè)、設施農(nóng)業(yè)，提高污水回收和處理率、增加林地、草地等自然用地是改善水質(zhì)惡化現(xiàn)狀、提升水生態(tài)環(huán)境狀況的必要措施。

□