郭榮幸　王超梁　陳濟(jì)民等

關(guān)鍵詞：多光譜遙感；水質(zhì)反演；Ｓｅｎｔｉｎｅｌ－２；反演模型；小浪底水庫(kù)

中圖分類號(hào)：Ｘ８２３；ＴＶ２１１．１＋ １ 文獻(xiàn)標(biāo)志碼：Ａ ｄｏｉ：１０．３９６９／ ｊ．ｉｓｓｎ．１０００－１３７９．２０２４．０１．０１６

引用格式：郭榮幸，王超梁，陳濟(jì)民，等．基于Ｓｅｎｔｉｎｅｌ－２ 多光譜遙感影像的小浪底水質(zhì)反演［Ｊ］．人民黃河，２０２４，４６（１２）：９３－９６，１０２．

水是人類賴以生存的重要自然資源。黃河是中華民族的母親河，也是我國(guó)北方九省（區(qū)）主要的工業(yè)、農(nóng)業(yè)和生活水源［１］ 。近年來(lái)隨著城市化的快速發(fā)展，工業(yè)廢水和生活污水的排放造成黃河流域水環(huán)境污染嚴(yán)重，已成為制約區(qū)域高質(zhì)量發(fā)展的主要因素，而快速高效的水質(zhì)監(jiān)測(cè)是黃河流域生態(tài)保護(hù)的關(guān)鍵。

傳統(tǒng)的水質(zhì)監(jiān)測(cè)以實(shí)地采樣為主，不僅周期長(zhǎng)、成本高，而且只能獲取采樣點(diǎn)的水質(zhì)信息，很難進(jìn)行長(zhǎng)時(shí)序、大范圍的監(jiān)測(cè)［２］ 。多光譜遙感技術(shù)可根據(jù)遙感波段信息反演水質(zhì)參數(shù)，降低監(jiān)測(cè)成本，提高監(jiān)測(cè)速度和質(zhì)量，為大范圍水質(zhì)監(jiān)測(cè)提供了一種新的方法［３］ 。章佩麗等［４］ 基于無(wú)人機(jī)多光譜影像進(jìn)行城市河道水質(zhì)參數(shù)反演，建立了水質(zhì)參數(shù)的反演模型并驗(yàn)證了模型的精確性和適用性；朱熹等［５］ 采用衛(wèi)星多光譜遙感數(shù)據(jù)及河道斷面采樣點(diǎn)數(shù)據(jù)，反演了上海市浦東新區(qū)河道的總磷、氨氮、溶解氧等水質(zhì)參數(shù)；Ｍｃｅｌｉｅｃｅ 等［６］ 基于多光譜影像對(duì)近岸海域葉綠素ａ 濃度與濁度進(jìn)行了反演；張立福等［７］ 采用結(jié)合多源遙感數(shù)據(jù)的星地協(xié)同水質(zhì)光譜在線監(jiān)測(cè)技術(shù)，提升了長(zhǎng)時(shí)序、大范圍水域水質(zhì)參數(shù)反演的精確度和穩(wěn)定性。本文以黃河小浪底水庫(kù)為研究對(duì)象，通過(guò)Ｓｅｎｔｉｎｅｌ－２ 獲取小浪底水庫(kù)多光譜遙感影像數(shù)據(jù)，結(jié)合采樣點(diǎn)實(shí)測(cè)水質(zhì)數(shù)據(jù)建立了ＣＯＤ、ＴＰ、ＴＮ、ＮＨ３ －Ｎ 等參數(shù)的反演模型，并驗(yàn)證了模型的精確度和穩(wěn)定性，進(jìn)而反演各水質(zhì)參數(shù)的空間分布，以期為黃河流域快速高效的水質(zhì)監(jiān)測(cè)提供參考。

１ 材料與方法

１．１ 研究區(qū)概況

小浪底水庫(kù)位于河南省洛陽(yáng)市與濟(jì)源市之間的黃河干流，水域面積２７２ ｋｍ^２，是黃河干流三門(mén)峽以下唯一能取得較大庫(kù)容的控制性工程，主要用于控制洪水、減緩下游河道淤積，兼顧供水、灌溉和發(fā)電。近年來(lái)，小浪底水庫(kù)水質(zhì)為地表水Ⅱ類和Ⅲ類［８］ 。

１．２ 采樣點(diǎn)分布和水質(zhì)參數(shù)分析

２０２２ 年５ 月５ 日在小浪底水庫(kù)出水口附近進(jìn)行水質(zhì)采樣，采樣點(diǎn)布置見(jiàn)圖１。采樣點(diǎn)間距約１００ ｍ，在采樣點(diǎn)用取水器取表層０．５ ｍ 深處水樣５００ ｍＬ，分別采用檢測(cè)試劑和儀器進(jìn)行水質(zhì)分析［９］ ，共獲取了１８?jìng)€(gè)采樣點(diǎn)的水質(zhì)樣本，采樣點(diǎn)ＣＯＤ、ＴＮ、ＴＰ、ＮＨ_３－Ｎ 質(zhì)量濃度見(jiàn)圖２。

１．３ 多光譜影像獲取和數(shù)據(jù)處理

遙感影像數(shù)據(jù)來(lái)源為搭載多光譜成像儀（ＭＳＩ）的高分辨率多光譜成像衛(wèi)星Ｓｅｎｔｉｎｅｌ－２，由軌道彼此相差１８０°的Ｓｅｎｔｉｎｅｌ－２Ａ 和Ｓｅｎｔｉｎｅｌ－２Ｂ 兩顆衛(wèi)星組成，主要包含從可見(jiàn)光到短波紅外的１３ 個(gè)工作波段，圖像采集周期為５ ｄ［１０］ 。本文采用從歐空局下載的２０２２年５ 月５ 日Ｓｅｎｔｉｎｅｌ－２Ａ 經(jīng)過(guò)輻射校正、大氣校正的具有更高亮度和對(duì)比度的Ｌ２Ａ 級(jí)遙感數(shù)據(jù)，并使用ＳＮＡＰ 軟件對(duì)遙感數(shù)據(jù)進(jìn)行重采樣、裁剪等處理。

１．４ 水質(zhì)參數(shù)反演模型構(gòu)建流程

水質(zhì)參數(shù)反演模型主要通過(guò)以下步驟建立：對(duì)Ｓｅｎｔｉｎｅｌ－２ 遙感數(shù)據(jù)進(jìn)行預(yù)處理，提取光譜反射率數(shù)據(jù)；分析實(shí)測(cè)水質(zhì)參數(shù)和光譜反射率數(shù)據(jù)之間的相關(guān)性，獲取相應(yīng)水質(zhì)參數(shù)的敏感波段及其組合；根據(jù)敏感波段及其組合，選擇經(jīng)驗(yàn)統(tǒng)計(jì)模型進(jìn)行水質(zhì)參數(shù)反演，建立各水質(zhì)參數(shù)的遙感反演定量模型；通過(guò)驗(yàn)證集測(cè)試反演模型，對(duì)模型精確度和穩(wěn)定性進(jìn)行評(píng)價(jià)，并反演水質(zhì)參數(shù)空間分布，模型構(gòu)建流程見(jiàn)圖３。

１．５ 反演模型評(píng)價(jià)標(biāo)準(zhǔn)

采用含量梯度法將１８ 個(gè)采樣點(diǎn)按２ ∶ １ 的比例劃分為建模集和驗(yàn)證集，水質(zhì)參數(shù)反演模型采用確定系數(shù)Ｒ２、均方根誤差ＲＭＳＥ、相對(duì)分析誤差ＲＰＤ 進(jìn)行精度評(píng)價(jià)，模型的準(zhǔn)確性越高，則相應(yīng)的Ｒ２ 和ＲＰＤ 值越大、ＲＭＳＥ 越小。模型穩(wěn)定且準(zhǔn)確時(shí)ＲＰＤ 的值大于２．０，模型穩(wěn)定性一般時(shí)ＲＰＤ 的值介于１．４ 和２．０ 之間，否則模型的預(yù)測(cè)能力較差［１１］ 。

２ 結(jié)果與分析

２．１ 小浪底水庫(kù)水質(zhì)光譜特征

采用ＳＮＡＰ 軟件對(duì)遙感數(shù)據(jù)進(jìn)行處理，獲取小浪底水庫(kù)采樣點(diǎn)光譜反射率曲線，見(jiàn)圖４。從圖４ 中可以看出，波長(zhǎng)在４４０～８７５ ｎｍ 范圍內(nèi)水體的光譜反射率呈波動(dòng)變化。４４０～５６０ ｎｍ 處反射峰出現(xiàn)是葉綠素吸收較弱以及水體中懸浮物的反射作用造成的，６６５ ｎｍ處反射波谷是葉綠素的強(qiáng)吸收作用引起的，７０５ ｎｍ 處反射峰是藻類植物的反射造成的，７８０ ｎｍ 處反射峰是懸浮物等物質(zhì)對(duì)光的反射作用造成的［１２］ 。綜合分析采樣點(diǎn)水質(zhì)光譜反射曲線可知，不同區(qū)域光譜變化趨勢(shì)總體一致，由于水質(zhì)參數(shù)濃度不同，因此反射曲線的峰谷高低和變化快慢不同。

２．２ 水質(zhì)參數(shù)和反射率相關(guān)性分析

將小浪底水庫(kù)各采樣點(diǎn)的水質(zhì)參數(shù)分別與光譜反射率進(jìn)行Ｐｅａｒｓｏｎ 相關(guān)性分析?？紤]水體遙感反演敏感波段分布和噪聲的影響，采用１ ０００ ｎｍ 以下波段進(jìn)行相關(guān)性分析，水質(zhì)參數(shù)不同波長(zhǎng)下與反射率的相關(guān)系數(shù)見(jiàn)圖５。由圖５ 可知，在波長(zhǎng)４５０ ～ ９００ ｎｍ 范圍內(nèi)，ＮＨ_３－Ｎ 與光譜各波段反射率正相關(guān)，相關(guān)系數(shù)在８４２ ｎｍ 處達(dá)到最大值０．４２，為中等相關(guān)性（０．４＜ｒ ＜０．６）［１３］ ；ＣＯＤ 與各波段反射率存在正相關(guān)和負(fù)相關(guān)，相關(guān)系數(shù)在７０５ ｎｍ 處達(dá)到最大值０．７９，為強(qiáng)相關(guān)性（０．６＜ｒ＜０．８）；ＴＰ 在６８０～７６０ ｎｍ 范圍內(nèi)與各波段反射率負(fù)相關(guān)，在其余波段正相關(guān)，相關(guān)系數(shù)在７８３ ｎｍ 處達(dá)到最大值０．３１，為較弱相關(guān)性（０．２＜ｒ＜０．４）；ＴＮ 與各波段反射率存在正相關(guān)和負(fù)相關(guān)，相關(guān)系數(shù)在７４０ ｎｍ處達(dá)到最大值０．５１，為中等相關(guān)性。

２．３ 水質(zhì)參數(shù)反演模型

分析水質(zhì)參數(shù)和光譜各波段反射率的相關(guān)性，篩選出各水質(zhì)參數(shù)相關(guān)性峰值對(duì)應(yīng)的光譜波段，ＣＯＤ、ＴＰ、ＴＮ 和ＮＨ_３－Ｎ 相關(guān)性最高的波段分別為７０５ ｎｍ（Ｂ５）、７８３ ｎｍ（Ｂ７）、７４０ ｎｍ（Ｂ６）和８４２ ｎｍ（Ｂ８）波段，建立水質(zhì)參數(shù)與所選波段之間反演模型（見(jiàn)表１），在各反演模型中選取確定系數(shù)最大的反演模型（見(jiàn)圖６）。從圖６ 可以看出，建模集數(shù)據(jù)點(diǎn)均勻分布在擬合曲線兩側(cè)，ＣＯＤ、ＴＰ 和ＴＮ 的反演模型確定系數(shù)Ｒ２ 均高于ＮＨ_３－Ｎ 反演模型的。

２．４ 水質(zhì)參數(shù)反演模型驗(yàn)證及評(píng)價(jià)

通過(guò)遙感反演模型對(duì)小浪底水質(zhì)參數(shù)進(jìn)行反演，用實(shí)測(cè)數(shù)據(jù)中的６ 組驗(yàn)證集數(shù)據(jù)與模型反演值進(jìn)行曲線擬合，擬合結(jié)果表明，各水質(zhì)參數(shù)反演值和驗(yàn)證集實(shí)測(cè)值均勻分布在擬合曲線兩側(cè)，擬合曲線確定系數(shù)Ｒ^２均高于反演模型的，且曲線斜率均在１ 附近，說(shuō)明各反演模型可用于水質(zhì)參數(shù)反演，各水質(zhì)參數(shù)反演模型精度評(píng)價(jià)見(jiàn)表２。從表２ 可以看出，ＣＯＤ 反演模型的確定系數(shù)Ｒ^２最大，ＲＰＤ 大于２．０，反演模型精確度和穩(wěn)定性最好；ＴＰ 和ＴＮ 的反演模型精確度和穩(wěn)定性較好；ＮＨ_３－Ｎ 的確定系數(shù)Ｒ^２較小，ＲＰＤ 小于２．０，模型精確度和穩(wěn)定性較差。可嘗試通過(guò)提高多光譜數(shù)據(jù)獲取的實(shí)時(shí)性和增加采樣點(diǎn)數(shù)量來(lái)提高反演精度。

２．５ 水質(zhì)參數(shù)質(zhì)量濃度分布圖反演

根據(jù)建立的ＣＯＤ、ＴＰ、ＴＮ 和ＮＨ_３－Ｎ 反演模型，對(duì)獲取的多光譜圖像各像元進(jìn)行估算，反演小浪底水質(zhì)參數(shù)質(zhì)量濃度分布，結(jié)果見(jiàn)圖７。由圖７ 可知，水庫(kù)出水口和部分邊界處ＣＯＤ 質(zhì)量濃度較高，往水庫(kù)內(nèi)部延伸ＴＮ、ＴＰ 和ＮＨ_３－Ｎ 質(zhì)量濃度高于邊緣處。

３ 結(jié)束語(yǔ)

以小浪底水庫(kù)Ｓｅｎｔｉｎｅｌ－２ 多光譜遙感影像和采樣點(diǎn)實(shí)測(cè)水質(zhì)參數(shù)為數(shù)據(jù)源，通過(guò)相關(guān)性分析篩選各水質(zhì)參數(shù)的敏感波段，選擇相關(guān)性好的波段進(jìn)行水質(zhì)參數(shù)反演，建立了反演模型，并通過(guò)驗(yàn)證集驗(yàn)證了各模型的精確度和穩(wěn)定性，進(jìn)而反演了各水質(zhì)參數(shù)的空間分布，結(jié)果表明：小浪底水庫(kù)ＣＯＤ 反演模型精確度和穩(wěn)定性較好，ＴＰ 和ＴＮ 反演模型次之，ＮＨ_３－Ｎ 的性能最差，模型可以用來(lái)反演小浪底水庫(kù)ＣＯＤ、ＴＰ 和ＴＮ 的濃度和空間分布，用于小浪底水庫(kù)和黃河流域水生態(tài)環(huán)境水質(zhì)參數(shù)大范圍快速監(jiān)測(cè)。

【責(zé)任編輯 呂艷梅】