穆紅波 徐燕飛 殷夢(mèng)杰

摘 要："珠海一號(hào)"衛(wèi)星高光譜數(shù)據(jù)具有高時(shí)空分辨率的特點(diǎn)，為采煤沉陷區(qū)水質(zhì)時(shí)空監(jiān)測(cè)提供了更高精度的數(shù)據(jù)源。葉綠素a濃度是評(píng)價(jià)水質(zhì)狀況的重要參數(shù)。為探究淮南采煤沉陷區(qū)葉綠素a濃度反演模型，以安徽淮南顧橋采煤沉陷區(qū)為研究對(duì)象，基于"珠海一號(hào)"高光譜影像和水樣測(cè)試數(shù)據(jù)，對(duì)采煤沉陷水域葉綠素a濃度反演模型進(jìn)行構(gòu)建，反演并分析了顧橋采煤沉陷區(qū)葉綠素a濃度空間分布特征及其成因。結(jié)果表明：基于珠海一號(hào)高光譜影像（b10+b19-b2）/（b10+b19+b2）波段組合方式構(gòu)建的三次函數(shù)模型精度和穩(wěn)定性較高，可以用于顧橋采煤沉陷區(qū)水質(zhì)的動(dòng)態(tài)監(jiān)測(cè);顧橋沉陷區(qū)水域總體葉綠素a濃度較高，高值區(qū)主要位于水域中部，且南部水域高于北部;農(nóng)田施肥和工礦生產(chǎn)、居民生活污水的排放是造成沉陷區(qū)內(nèi)葉綠素a濃度上升的主要原因。

關(guān)鍵詞：水質(zhì)遙感;珠海一號(hào)衛(wèi)星;高光譜影像;葉綠素a;采煤沉陷區(qū)

Abstract：The hyperspectral data OHS-1 satellite，with the high spatial and temporal resolution， provides a more accurate data source for the spatial and temporal monitoring of water quality in mining subsidence area. Chlorophyll-a concentration is an important parameter to evaluate water quality. To explore the inversion model of chlorophyll-a concentration in Huainan coal mining subsidence area， the inversion model of chlorophyll-a concentration in the coal mining subsidence area was constructed， with Guqiao coal mining subsidence area in Huainan， Anhui province as the research object， and the spatial distribution characteristics of chlorophyll-a and its origin in the subsidence area were analyzed based on the hyperspectral image of OHS-1 satellite and water sampling data. The results show that the cubic function model based on （b10+b19-b2）/（b10+b19+b2）band combination from hyperspectral image of OHS-1 satellite has high inversion accuracy and stability， and can be used for dynamic monitoring of water quality in Guqiao coal mining subsidence area. The overall chlorophyll-a concentration in the water of subsidence area is high， mainly located in the middle part of the area， which in the southern part is higher than in the northern part. The main reasons for the increase of chlorophyll-a concentration in the subsidence area are the fertilization of farmland， the discharge of the industrial and mining production and the domestic sewage.

Key words：water quality remote sensing; OHS-1 satellite; hyperspectral image; Chlorophyll-a;coal mining subsidence

我國(guó)東部高潛水位礦區(qū)，常年大規(guī)模地下煤炭的開采形成了大面積的沉陷水域，嚴(yán)重影響了當(dāng)?shù)氐霓r(nóng)業(yè)生產(chǎn)和居民生活，沉陷水域的水環(huán)境問題成為限制礦區(qū)生態(tài)環(huán)境修復(fù)和區(qū)域經(jīng)濟(jì)可持續(xù)發(fā)展的重要因素[1]85。水體水質(zhì)參數(shù)的有效監(jiān)測(cè)和研究是開展水環(huán)境治理和修復(fù)的重要前提。傳統(tǒng)的通過采集水樣進(jìn)行水質(zhì)測(cè)試分析的方法盡管精度較高，但由于受采樣點(diǎn)數(shù)量和代表性的限制，難以實(shí)現(xiàn)水質(zhì)參數(shù)空間化、時(shí)效性、代表性的監(jiān)測(cè)和時(shí)空變化研究[2-3]。相比于傳統(tǒng)監(jiān)測(cè)分析方法，遙感技術(shù)具有成本低、范圍廣、時(shí)效性強(qiáng)的優(yōu)勢(shì)，利用遙感影像進(jìn)行水域水質(zhì)參數(shù)的反演，可以實(shí)現(xiàn)水質(zhì)全區(qū)域的空間分布及其變化分析。

目前，水質(zhì)遙感涉及到的參數(shù)主要包括葉綠素a（Chlorophyll-a，Chla）、透明度、濁度、總氮、總磷、富營(yíng)養(yǎng)化指數(shù)等直接和間接相關(guān)參數(shù)。所使用的數(shù)據(jù)源從多光譜遙感數(shù)據(jù)、非成像高光譜數(shù)據(jù)到高光譜成像光譜數(shù)據(jù)，發(fā)展了經(jīng)驗(yàn)?zāi)Ｐ?、半?jīng)驗(yàn)/半分析模型和分析模型三種反演方法。針對(duì)內(nèi)陸地表水水質(zhì)遙感，國(guó)內(nèi)外學(xué)者建立了眾多的水質(zhì)參數(shù)遙感反演模型，然而，由于內(nèi)陸水體光學(xué)特征的復(fù)雜性，目前還沒有普適性的水質(zhì)遙感反演模型[4-6]。近年來(lái)，一些學(xué)者在對(duì)內(nèi)陸湖泊水體水質(zhì)參數(shù)遙感反演的基礎(chǔ)上，對(duì)我國(guó)東部采煤沉陷水域水質(zhì)參數(shù)的遙感反演進(jìn)行了研究。文獻(xiàn)[7]以兗州礦濟(jì)寧礦區(qū)為研究區(qū)，基于Hyperion高光譜數(shù)據(jù)，通過對(duì)礦業(yè)擾動(dòng)區(qū)水體組分敏感波段的分析，構(gòu)建礦業(yè)擾動(dòng)的反演模型，實(shí)現(xiàn)了礦業(yè)擾動(dòng)區(qū)地表水體的遙感反演。文獻(xiàn)[8]綜合應(yīng)用地面測(cè)量光譜和無(wú)人機(jī)高光譜數(shù)據(jù)，對(duì)工礦城市尾礦庫(kù)、河流、湖泊等典型區(qū)域的水質(zhì)參數(shù)進(jìn)行了反演。文獻(xiàn)[9]以淮南顧橋采煤沉陷水域?yàn)檠芯繉?duì)象，利用HJ-1 A/B衛(wèi)星的和的影像數(shù)據(jù)，結(jié)合同步水體波譜信息和現(xiàn)場(chǎng)水質(zhì)監(jiān)測(cè)數(shù)據(jù)，建立了針對(duì)淮南釆煤沉陷水域水環(huán)境特點(diǎn)的遙感反演模型，并對(duì)采煤沉陷水域富營(yíng)養(yǎng)化現(xiàn)狀進(jìn)行了評(píng)價(jià)。

文獻(xiàn)[1]89利用HJ-1衛(wèi)星數(shù)據(jù)，分別對(duì)淮南顧橋和潘集兩個(gè)不同類型的沉陷區(qū)水體進(jìn)行了相關(guān)水質(zhì)指標(biāo)的遙感反演并進(jìn)行了富營(yíng)養(yǎng)化評(píng)價(jià)，研究表明顧橋沉陷水域富營(yíng)養(yǎng)化程度明顯高于潘集沉陷水域。盡管前人已經(jīng)在顧橋采煤沉陷區(qū)水質(zhì)遙感監(jiān)測(cè)方面取得了一定的成果，但上述研究主要是基于中等分辨率遙感衛(wèi)星數(shù)據(jù)開展的，相較于大型湖泊、河流的水質(zhì)反演，采煤沉陷區(qū)水體光譜特征更為復(fù)雜，且水域面積較小，鄰近像元效應(yīng)影響較為嚴(yán)重，更加需要高分辨率的衛(wèi)星遙感數(shù)據(jù)。2018年12月17日，“珠海一號(hào)”商業(yè)航天高光譜衛(wèi)星數(shù)據(jù)首次對(duì)外發(fā)布，其數(shù)據(jù)具有高光譜、高空間、高時(shí)間分辨等特點(diǎn)，開啟了我國(guó)商業(yè)高光譜遙感定量新時(shí)代，為小范圍高分辨率水域水質(zhì)時(shí)空監(jiān)測(cè)提供了數(shù)據(jù)源，在內(nèi)陸水體遙感監(jiān)測(cè)中具有較大潛力[10-11]。

葉綠素a作為水體光合生物諸如藍(lán)藻、綠藻等生物的重要組成成分，是衡量水體初級(jí)生產(chǎn)力和富營(yíng)養(yǎng)化程度最基本的指標(biāo)，也是水質(zhì)狀態(tài)評(píng)價(jià)的關(guān)鍵參數(shù)之一，成為國(guó)內(nèi)外學(xué)者水質(zhì)遙感反演研究的主要目標(biāo)[12-16]。本文以安徽淮南顧橋采煤沉陷水域?yàn)檠芯繉?duì)象，基于“珠海一號(hào)”高光譜數(shù)據(jù)，結(jié)合野外實(shí)測(cè)光譜和水樣測(cè)試分析，通過對(duì)葉綠素a濃度敏感性波段及其組合的相關(guān)分析，構(gòu)建遙感反演模型，對(duì)研究區(qū)葉綠素a空間分布特征及其成因加以分析，為該區(qū)域水環(huán)境的監(jiān)測(cè)和治理提供科學(xué)指導(dǎo)。

1 研究區(qū)概況

顧橋礦位于安徽省淮南市鳳臺(tái)縣以東約20km，地處潘謝礦區(qū)中西部，東與丁集礦為鄰，西與張集礦相接，北與顧北礦相通，井田面積約為106km2（見圖1）。礦區(qū)內(nèi)地下水位平均埋深1～2m，長(zhǎng)期的煤炭開采形成了大面積沉陷積水區(qū)，平均水深約為4m。經(jīng)遙感影像解譯，截止2020a，整個(gè)積水區(qū)域面積達(dá)到11.6km2。

礦區(qū)內(nèi)主要土地利用類型為耕地，采煤沉陷水域的主要利用方式為水產(chǎn)養(yǎng)殖和光伏發(fā)電。由于沉陷水域主要由農(nóng)田、村莊塌陷而來(lái)，周圍有發(fā)電廠、煤矸石堆積，各種農(nóng)業(yè)污水、生活及工業(yè)廢水進(jìn)入積水區(qū)域，造成了一定程度的污染。同時(shí)，不合理的水產(chǎn)養(yǎng)殖也不同程度地影響著水體水質(zhì)狀況[1]87-88。

2 數(shù)據(jù)與方法

（1）樣品采集與測(cè)試

2019年5月23日對(duì)顧橋沉陷區(qū)水質(zhì)進(jìn)行了實(shí)地采樣，采樣點(diǎn)分布如圖2所示，采樣當(dāng)天天氣晴朗，水面平靜。此次采樣共采集有效水體樣本20個(gè)，同步利用ASD野外光譜輻射儀，采用水體水面以上光譜測(cè)量方法[17]對(duì)各樣點(diǎn)的輻射數(shù)據(jù)進(jìn)行測(cè)量，共獲得水體光譜曲線1 400條。水體樣本按照相關(guān)要求保存在棕色瓶中帶回實(shí)驗(yàn)室進(jìn)行測(cè)試分析。其中，2/3的測(cè)試數(shù)據(jù)用來(lái)進(jìn)行模型的構(gòu)建，1/3的數(shù)據(jù)進(jìn)行模型反演精度的評(píng)價(jià)。

（2）遙感數(shù)據(jù)處理

本文收集了與采樣時(shí)間接近的2019年5月4日的一景成像質(zhì)量?jī)?yōu)良的“珠海一號(hào)”高光譜衛(wèi)星數(shù)據(jù)， 該數(shù)據(jù)具有256個(gè)波段（有效成像波段數(shù)32）， 波長(zhǎng)范圍為400nm～1 000nm，光譜分辨率為2.5nm，空間分辨率為10m，幅寬為150km，重訪周期為2d。

由于在遙感圖像中，水體表現(xiàn)為弱信息，因此在對(duì)沉陷區(qū)水質(zhì)參數(shù)進(jìn)行反演時(shí)，對(duì)遙感圖像的預(yù)處理就顯得尤為重要[1]87，主要包括遙感圖像的幾何校正、輻射定標(biāo)、大氣校正和圖像裁剪。

此外，在進(jìn)行水質(zhì)相關(guān)指標(biāo)反演之前還需對(duì)研究區(qū)水體范圍進(jìn)行提取。因水體在藍(lán)光波段對(duì)入射能量具有較強(qiáng)反射峰，在近紅外、短波紅外幾乎吸收全部能量，國(guó)內(nèi)相關(guān)學(xué)者根據(jù)上述水體光譜特征先后提出了單波段閾值法、多波段譜間關(guān)系法、歸一化水體指數(shù)法、改進(jìn)的歸一化水體指數(shù)法等來(lái)進(jìn)行水體信息的提取[18-19]。根據(jù)在該地區(qū)的實(shí)驗(yàn)，使用歸一化水體指數(shù)方法能很好的提取出沉陷積水區(qū)范圍。

（3）反演模型構(gòu)建

利用實(shí)測(cè)采樣點(diǎn)，對(duì)珠海一號(hào)高光譜數(shù)據(jù)每一波段反射率值和葉綠素a濃度測(cè)試值進(jìn)行相關(guān)性分析，結(jié)果如圖3所示。

由圖3可知，葉綠素a濃度與波段b2、b10、b19相關(guān)性較強(qiáng)，所以可以使用上述波段進(jìn)行葉綠素a濃度的反演。由于單波段與葉綠素a濃度進(jìn)行相關(guān)性分析時(shí)其相關(guān)性可能較弱，所以一般方法是將幾個(gè)較強(qiáng)相關(guān)性的波段進(jìn)行組合，以提高與葉綠素a濃度的相關(guān)性，進(jìn)而增加反演模型的精度，常用的增強(qiáng)方式有波段間比值、波段間歸一化、波段間求和、波段間差值等。據(jù)前人研究結(jié)果表明，通過上述增強(qiáng)方法后的波段運(yùn)算結(jié)果在一定程度上能增加與葉綠素a濃度的相關(guān)性。因此，這里將相關(guān)性較強(qiáng)的波段b2、波段b10、波段b19進(jìn)行上述組合，最后將組合結(jié)果與采樣點(diǎn)葉綠素a濃度進(jìn)行相關(guān)性分析，結(jié)果如表1所示。

由表1可知，葉綠素a濃度與（b10+b19-b2）/（b10+b19+b2）組合方式相關(guān)性最高，所以選用此波段組合進(jìn)行葉素綠a濃度反演模型的構(gòu)建。分別采用常用的線性函數(shù)、指數(shù)函數(shù)、二次函數(shù)、三次函數(shù)模型對(duì)每個(gè)采樣點(diǎn)的波段組合值和葉綠素a濃度值進(jìn)行曲線擬合，結(jié)果如圖4所示。

擬合效果的評(píng)價(jià)一般使用R2與調(diào)整后R2來(lái)衡量，R2越接近1，表示效果越好，調(diào)整后R2越接近1模型穩(wěn)定性越高。由圖4可知，三次函數(shù)模型的R2最高，調(diào)整后R2也達(dá)到了0.603，因此選擇三次函數(shù)模型作為本研究區(qū)葉綠素a濃度的最佳反演模型。

（4）反演模型評(píng)價(jià)

對(duì)利用高光譜數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)建立的各水質(zhì)指標(biāo)反演模型，使用剩余1/3樣本進(jìn)行模型精度評(píng)價(jià)，計(jì)算其平均絕對(duì)誤差（Mean Absolute Error，MAE）、均方根誤差（Root Mean Square Error，RMSE），并作出其誤差圖（見圖5）、實(shí)測(cè)值與模型估算值對(duì)比圖（見圖6）。平均絕對(duì)誤差（MAE）是將各次測(cè)量的誤差取絕對(duì)值后再求平均值，可以避免誤差相互抵消的問題，可以準(zhǔn)確反映實(shí)際估算誤差的大小。均方根誤差（RMSE）是模型預(yù)測(cè)值與真實(shí)觀測(cè)值偏差的平方與觀測(cè)次數(shù)的n比值的平方根，主要用來(lái)評(píng)價(jià)模型估算值與真實(shí)值之間的偏差，同時(shí)在一定程度上也能代表估算模型的穩(wěn)定性。一般來(lái)說(shuō)，均方根誤差越小，模型精度和穩(wěn)定性越好。由圖5可知，葉綠素a濃度MAE為1.132，RMSE為1.231，因此，該模型可以用于顧橋沉陷積水區(qū)葉綠素a濃度的反演。

3 結(jié)果分析

基于上述葉綠素a濃度反演模型得到研究區(qū)葉綠素濃度空間分布，并將其劃分為五級(jí)，為便于闡述將沉陷水域標(biāo)記為五個(gè)小區(qū)域，沉陷水域西北部分因大量光伏電池的鋪設(shè)，未參與反演，如圖7所示。

由圖3.1可知，葉綠素a濃度在整個(gè)水域濃度都比較高，尤其是水域中部，且南部Ⅲ、Ⅳ區(qū)域濃度相對(duì)北部高值區(qū)范圍更大，北部I、Ⅱ區(qū)域以及西部V區(qū)域濃度相對(duì)較低、范圍較小，與裴文明等研究結(jié)果相一致[1]88-89。這主要與當(dāng)?shù)氐霓r(nóng)業(yè)生產(chǎn)有關(guān)，5月份為農(nóng)田播種季節(jié)，大量農(nóng)業(yè)化肥使用過程中，氮磷等元素通過地表徑流匯集到沉陷積水區(qū)內(nèi)，并且其南部靠近發(fā)電廠、礦區(qū)和居民區(qū)，大量工業(yè)廢水和生活污水的流入也導(dǎo)致富營(yíng)養(yǎng)化程度加重。此外鄰近沉陷水域中煤矸石填充的道路南側(cè)，葉綠素a濃度也相對(duì)較高。

4 結(jié)論

本文以淮南顧橋采煤沉陷區(qū)為研究對(duì)象，基于“珠海一號(hào)”高光譜數(shù)據(jù)各波段反射率值與葉綠素a濃度實(shí)測(cè)值的相關(guān)性分析，確立最佳波段組合，建立了研究區(qū)葉綠素a濃度的反演模型，對(duì)研究區(qū)葉綠素a濃度進(jìn)行了空間反演和分析。

（1）利用“珠海一號(hào)”高光譜數(shù)據(jù)能夠得到更高分辨率的水質(zhì)反演結(jié)果，為采煤沉陷區(qū)水質(zhì)監(jiān)測(cè)和時(shí)空變化研究提供了有力支持。對(duì)于顧橋采煤沉陷區(qū)來(lái)說(shuō)，葉綠素a濃度與“珠海一號(hào)”高光譜數(shù)據(jù)b2、b10、b19波段相關(guān)性較強(qiáng)，基于（b10+b19-b2）/（b10+b19+b2）組合所構(gòu)建的三次函數(shù)模型反演效果最佳，反演模型精度和穩(wěn)定性較高。但模型在整個(gè)礦區(qū)的適用性以及季節(jié)性差異還有待進(jìn)一步驗(yàn)證。

（2）整個(gè)顧橋采煤沉陷區(qū)葉綠素a濃度都比較高，說(shuō)明水體的富營(yíng)養(yǎng)化程度較高，高值區(qū)主要位于水域中部及煤矸石充填的道路南側(cè)，且南部水域高于北部。這一方面與農(nóng)業(yè)生產(chǎn)過程中化肥使用有關(guān)，另一方面與煤矸石堆積、工礦生產(chǎn)廢水和生活污水的排放有關(guān)。

參考文獻(xiàn)：

[1] 裴文明，張慧，姚素平，等.淮南礦區(qū)不同類型沉陷水域水質(zhì)遙感反演和時(shí)空變化分析[J].煤田地質(zhì)與勘探，2018，46（3）：85-90.

[2] 徐鵬飛，毛峰，金平斌，等.基于GF1_WFV的千島湖水體葉綠素a濃度時(shí)空變化[J].中國(guó)環(huán)境科學(xué)，2020，40（10）：4 580-4 588.

[3] 侯鵬飛.城鎮(zhèn)化背景下滇池葉綠素a濃度時(shí)空變化特征研究[D].昆明：云南師范大學(xué)，2020.

[4] 段洪濤，羅菊花，曹志剛，等.流域水環(huán)境遙感研究進(jìn)展與思考[J].地理科學(xué)進(jìn)展，2019，38（8）：1 182-1 195.

[5] 曹引，冶運(yùn)濤，趙紅莉，等.內(nèi)陸水體水質(zhì)參數(shù)遙感反演集合建模方法[J].中國(guó)環(huán)境科學(xué)， 2017， 37（10）：3 940-3 951.

[6] DANIEL O，ANATOLY G，VITTORIO E B，et al.Review of constituent retrieval in optically deep and complex waters from satellite imagery[J].Remote Sensing of Environment，2012，118（4）：116-126.

[7] 韓留生.礦業(yè)擾動(dòng)區(qū)水質(zhì)參數(shù)高光譜遙感反演[D].青島：山東科技大學(xué)，2011.

[8] 彭令，梅軍軍，王娜，等.工礦業(yè)城市區(qū)域水質(zhì)參數(shù)高光譜定量反演[J].光譜學(xué)與光譜分析，2019，39（9）：2 921-2 928.

[9] 葉圓圓.基于RS淮南采煤沉陷水域水質(zhì)實(shí)時(shí)監(jiān)測(cè)技術(shù)研究[D].淮南：安徽理工大學(xué)，2014.

[10] 洪韜.珠海一號(hào)高光譜衛(wèi)星在內(nèi)陸湖泊監(jiān)測(cè)中的應(yīng)用[J].衛(wèi)星應(yīng)用，2019，（8）：19-22.

[11] 尹慧新，樊彥國(guó).基于“珠海一號(hào)”高光譜衛(wèi)星的洱海水體提取方法研究[J].航天返回與遙感，2020（4）：118-127.

[12] 陶然，彭金嬋，張豪，等.內(nèi)陸水體葉綠素a濃度遙感監(jiān)測(cè)方法研究進(jìn)展[J].地理信息世界，2019，26（4）：44-53.

[13] 潘應(yīng)陽(yáng)，國(guó)巧真，孫金華.水體葉綠素a濃度遙感反演方法研究進(jìn)展[J].測(cè)繪科學(xué)，2017，42（1）：43-48.

[14] 宋挺，周文鱗，劉軍志，等.利用高光譜反演模型評(píng)估太湖水體葉綠素a濃度分布[J].環(huán)境科學(xué)學(xué)報(bào)，2017，37（3）：888-899.

[15] 封紅娥，李家國(guó)，朱云芳，等.GF-1與Landsat8水體葉綠素a濃度協(xié)同反演——以太湖為例[J].國(guó)土資源遙感，2019，31（4）：182-189.

[16] WANG T，HUANG M T，WANG Z J.Estimation of chlorophyll-a Concentration of lakes based on SVM algorithm and Landsat 8 OLI images[J].Environmental Science and Pollution Research，2020，27：14 977-14 990.

[17] 唐軍武，田國(guó)良，汪小勇，等.水體光譜測(cè)量與分析Ⅰ：水面以上測(cè)量法[J].遙感學(xué)報(bào)，2004，8（1）：37-44.

[18] 朱鶴.遙感技術(shù)在地表水源地水體監(jiān)測(cè)中的應(yīng)用研究[D].北京：中國(guó)水利水電科學(xué)研究院，2013.

[19] 王航，秦奮.遙感影像水體提取研究綜述[J].測(cè)繪科學(xué)，2018，43（5）：23-32.

（責(zé)任編輯：李 麗，范 君）