周嘉旺
(福建省水利水電勘測設計研究院,福建 福州 350001)
九龍江北溪橋閘位于九龍江北溪河口,庫區作為福建省重要的水源地,涉及漳州第二水廠水源保護區、漳州市福糖水廠水源保護區、北溪廈門市引水水源保護區、龍海市自來水廠水源保護區等4個水源保護區,承擔著給廈門、漳州供水重擔,其水質安全極其重要。但受到九龍江北溪流域天然來水量較小、流域內人口較多等因素影響,北溪橋閘庫區的水質安全隱患較大,甚至偶發水華現象。因此,研究北溪橋閘庫區污染分布對水源地保護、流域綜合治理都具有重要意義。
CODMn(高錳酸鹽指數)是反映水體中有機和無機可氧化物質污染的重要指標,常用在地表水、飲用水的監測中,同時CODMn也是富營養化評價的重要指標。
CDOM即有色可溶性有機物,是光學上可測量的水中溶解性有機物的,主要由氨基酸、糖、脂肪酸類、酚類等組成,在紫外和藍光范圍具有強烈的吸收,在黃色波段吸收較小,隨著CDOM的增長,水體顏色逐漸由藍色過渡到綠色、黃綠色、黃色。根據相關研究表明,在天然水體中,CDOM與CODMn具有較好的相關性,使CODMn的遙感監測具有可行性。
OLI影像來自landsat8遙感衛星,其波段主要包括海岸波段、紅、綠、藍、近紅外、短波紅外等7個波段以及全色波段和卷云波段。除全色波段外,其它波段空間分辨率為30 m。通過輻射定標以及大氣校正后的OLI影像,可從各波段讀取實測點的反射率,獲取各實測點的光譜信息,用于定量分析。
通過分析各實測點的CODMn值和遙感光譜反射率,找到適合研究區域的波段組合方式,并建立二者的定量模型,從而實現遙感監測。
九龍江北溪為九龍江干流,是福建省第二大江,全長約274 km,流經龍巖、漳州、廈門等市。北溪橋閘位于九龍江北溪河口,是九龍江北溪干流末端閘壩,控制流域面積約9640 km2,庫區內涉及4個水源保護區,是福建省重要水源地。為了了解北溪橋閘庫區內水質情況,2018年3月6日進行水樣監測,共選取7個監測點位,詳見表1。

表1 北溪橋閘庫區監測點位分布
OLI影像從地理空間數據云獲取??紤]到landsat8衛星一個運行周期為16 d,且衛星影像質量受大氣影響較大,在選取影像的時候,應考慮成像日期以及云量的影響。本次OLI影像的數據標識為LC81190432018070LGN00,其成像日期為2018年3月11日,云量為2.98%,影像數據較好。
所獲取的影像投影坐標系為WGS84 UTM 50,已經過系統性的幾何校正。查看影像像元值為整型,為DN值,需通過輻射定標及大氣校正轉為反射率。采用ENVI5.3中Radiometric Calibration工具對所需影像及波段進行輻射定標,并進一步進行FLAASH大氣校正:影像平均地面高程設為0.05 km,大氣模型參數選擇Tropical(低緯度),氣溶膠模型選擇urban,氣溶膠反演方法選擇2-band(K-T)法。
在處理后的影像上繼續提取研究區域水系。影像水系提取主要有MNDWI(改進型水體歸一化指數)以及NDWI(水體歸一化指數)兩種方法。相比起NDWI,MNDWI采用綠色和短波紅外來增強無冰水面要素,同時降低了人工建筑的響應??紤]到研究區域河岸有較長的防洪堤、且水面架設橋梁較多,人工建筑影響較多,采用MNDWI來提取水系更加合適,MNDWI計算公式如下。
(1)
最后在arcgis Pro中采用柵格裁剪工具,對影像提取的水系進行進一步的修剪,獲得最終的研究范圍。


表2 各水樣CODMn監測值與反射率

表3 波段及波段組合與監測值的相關性(節選)

圖1 庫區高錳酸鹽指數預測分布
將驗證樣本的4、6、7波段反射率帶入訓練樣本的線性擬合公式,獲得驗證樣本CODMn濃度預測值,并將預測值與實測值進行對比,發現相對誤差較小,用本方法預測庫區CODMn濃度分布結果較可靠。驗證樣本預測結果分析見表4,線性擬合如下:

(2)
采用arcgis Pro的柵格計算器,將研究范圍內的各波段通過上述公式組合,獲得九龍江北溪橋閘庫區的CODMn濃度分布,可發現庫區內CODMn平均值約1.71 mg/L,且大部分地區可達到地表水Ⅱ類標準,庫區CODMn較高的區域主要分布在龍津溪口,此外受沿岸農田面源散排以及部分河岸水土流失的影響,部分沿岸區域CODMn也較高。
本文基于OLI影像及實測水質數據,利用遙感影像對污染物特征光譜的識別,反演了九龍江北溪橋閘庫區內的高錳酸鹽指數分布(圖1),并判斷出庫區內主要污染分布及來源。雖然實測和預測數據表明,龍津溪口水質未超過Ⅱ類標準,但結合相關工程經驗,當CODMn較高時(>3 mg/L),總磷總氮容易超標。此外,郭坑大橋至北溪橋閘區間,CODMn分布呈增大趨勢,說明沿岸農田面源的散排也是庫區水質的重要影響因素。因此,為進一步保障九龍江北溪橋閘水源地的供水水質安全,建議可針對龍津溪開展小流域綜合整治,并對庫區沿岸農田采取節水灌溉、測土配方肥等手段進行治理。