摘 要:懸浮泥沙濃度(SSC)是重要的水質(zhì)監(jiān)測參數(shù)。以寧夏清水河王團(tuán)段為研究對象,基于2022 年8 月26 日至11 月5 日原位高光譜數(shù)據(jù)與SSC 構(gòu)建了適用于極高含沙量水體的SSC 反演模型,選用單波段模型、波段差值模型、波段比值模型、二元線性模型4 種經(jīng)驗(yàn)?zāi)P停瑢Ρ绕渚?。結(jié)果表明:針對極高含沙量水體,單波段模型擬合效果不佳,模型決定系數(shù)(R2 ) 不足0.25。波段差值模型、波段比值模型、二元線性模型一定程度上可以降低噪聲的影響,提高模型R2。整體上,波段差值模型效果較好,基于650~720nm 波段反射率與560~700 nm 波段反射率構(gòu)成的波段差值模型的R2大于0.40,其中R687-R685 模型的R2最大,為0.76。
關(guān)鍵詞:原位高光譜;懸浮泥沙濃度;反演模型;寧夏清水河王團(tuán)段
中圖分類號:TV145+ .3 文獻(xiàn)標(biāo)志碼:A doi:10.3969/ j.issn.1000-1379.2025.04.005
引用格式:侯崢,侯英卓,邢前國,等.基于原位高光譜的極高懸浮泥沙濃度監(jiān)測模型研究[J].人民黃河,2025,47(4):28-31,37.
0 引言
受黃土高原水土流失與降水沖刷等影響,黃河懸浮泥沙濃度(SSC)極高。2022 年《黃河泥沙公報(bào)》顯示,黃河最大含沙量超過450 g/ L。SSC 是重要的水質(zhì)監(jiān)測參數(shù),其與水體濁度密切相關(guān),高SSC 水體的濁度通常較高[1] 。懸浮泥沙是污染物的重要載體[2] ,水體SSC 越高,污染物含量越高,太陽光在水體中的衰減速度越快,從而對水體浮游生物與底棲生物的生存產(chǎn)生不利影響[2-3] 。因此,開展黃河水體SSC 監(jiān)測,對提高黃河水質(zhì)分析與生態(tài)管理效率至關(guān)重要。
水體SSC 監(jiān)測方法主要有野外調(diào)查與遙感監(jiān)測兩種,其中遙感監(jiān)測方法具有監(jiān)測時(shí)序長、省時(shí)省力等優(yōu)點(diǎn),已被廣泛應(yīng)用于海水、湖泊及河流等多種水體的SSC 監(jiān)測[4-6] 。目前已有眾多SSC 遙感監(jiān)測模型被提出,主要包括經(jīng)驗(yàn)?zāi)P?、半分析模型與神經(jīng)網(wǎng)絡(luò)模型等[7-10] 。Qiu 等[11] 利用經(jīng)驗(yàn)?zāi)P蛯?shí)現(xiàn)了1986—2022年黃河水體SSC 反演,反演結(jié)果的平均相對誤差為39%,約97%河段的SSC 呈下降趨勢;Kong 等[9] 利用半分析模型對渤海灣水體SSC 進(jìn)行了反演研究,發(fā)現(xiàn)針對較渾濁水體SSC 的反演精度較高;Jiang 等[10] 對中山市水體SSC 進(jìn)行了反演研究,發(fā)現(xiàn)相比經(jīng)驗(yàn)?zāi)P?,在反演渾濁水體SSC 時(shí)神經(jīng)網(wǎng)絡(luò)模型沒有顯著優(yōu)勢。值得注意的是,目前被提出的眾多模型僅適用于低或高含沙量水體的SSC 反演,例如:章敏超等[12] 研究的舟山近岸水體SSC 范圍為0~1.2 g/ L,Hou 等[13]研究的長江中下游湖泊SSC 最高值為0.3 g/ L,Qiu等[11] 研究的黃河SSC 最高值為31 g/ L。目前,針對極高含沙量(SSC>100 g/ L)水體的SSC 遙感反演模型相對較少。原位高光譜數(shù)據(jù)是近地面采集的數(shù)據(jù),具有較多的波段信息,其數(shù)據(jù)質(zhì)量較少受大氣的影響,可反映不同水體的含沙量[14-17] 。因此,本文基于野外實(shí)測的極高SSC 數(shù)據(jù)與同步的原位高光譜數(shù)據(jù)構(gòu)建極高SSC 反演模型,實(shí)現(xiàn)對極高含沙量水體SSC 的監(jiān)測,以期為黃河泥沙預(yù)警與水沙管理提供支撐。
1 研究區(qū)概況與數(shù)據(jù)來源
1.1 研究區(qū)概況
寧夏位于我國西北部黃河中上游地區(qū),是黃土高原與內(nèi)蒙古高原的過渡帶,其西北和東北部與騰格里沙漠、烏蘭布和沙漠及毛烏素沙地相鄰[18] 。寧夏氣候類型為溫帶大陸性干旱半干旱氣候,常年干旱少雨,土地沙漠化較嚴(yán)重。清水河是寧夏境內(nèi)流入黃河的流域面積最大、長度最長的支流。清水河流域水土流失嚴(yán)重,降水時(shí)空分布不均。清水河含沙量非常高,實(shí)測多年平均含沙量達(dá)216 g/ L[19] ,平均輸沙模數(shù)為3 160 t/ km2。清水河流域年輸沙量分布極不均勻,高輸沙量多出現(xiàn)在汛期,全年80%以上的輸沙量集中在7—8 月,這與流域降水年內(nèi)分配密切相關(guān)。本研究選取寧夏清水河王團(tuán)段為研究區(qū)域,該區(qū)域水體平均含沙量超過100 g/L,是典型的極高SSC 水體。
1.2 數(shù)據(jù)來源
所用數(shù)據(jù)包括實(shí)測水體原位高光譜數(shù)據(jù)與實(shí)測SSC數(shù)據(jù)。原位高光譜數(shù)據(jù)通過ZK-UVIR-I 型原位光譜水質(zhì)在線監(jiān)測儀獲取,光譜測量范圍為350~950 nm,光譜分辨率為0.5 nm,采集時(shí)段為2022 年8 月26 日至11月5 日。水體SSC 數(shù)據(jù)通過人工野外采樣與室內(nèi)測量獲得,其采樣時(shí)段同原位高光譜數(shù)據(jù)。人工野外采樣采用水邊一線水面一點(diǎn)法,在河流平水期(含沙量變化平緩)每日8 時(shí)采樣1 次,在河流洪水期(含沙量變化劇烈)每日采樣7~10 次。室內(nèi)SSC 測量采用過濾法,具體步驟如下:將水樣注入量筒,測得水樣體積;待泥沙沉淀后,倒出上部清水,將沉淀的泥沙倒在濾紙上;水分下滲后,在100~110 ℃烘箱中烘干泥沙2 h;稱取泥沙質(zhì)量,泥沙質(zhì)量與水樣體積之比即為水體SSC。
經(jīng)過質(zhì)量控制,共獲得47 組水體原位高光譜數(shù)據(jù)(反射率)與SSC 數(shù)據(jù)。不同水體的光譜響應(yīng)存在很大差異,部分水體光譜反射率在短波波段與長波波段存在零值,故選取波長在400~900 nm 范圍內(nèi)的數(shù)據(jù)。圖1為水體不同SSC 對應(yīng)的光譜曲線。SSC 為0~5 g/ L 時(shí),SSC 越高各波段反射率越高;SSC>5 g/ L 時(shí),SSC 越高各波段反射率越低。
2 研究方法
基于水體原位高光譜與SSC 的關(guān)系,構(gòu)建極高SSC 反演模型。水體原位高光譜數(shù)據(jù)主要采用光譜分辨率重采樣至1、5、10 nm 的數(shù)據(jù),所用模型包括單波段模型、波段差值模型、波段比值模型、二元線性模型。單波段模型基于單波段反射率與SSC 構(gòu)建,波段差值模型基于單波段反射率兩兩差值與SSC 構(gòu)建,波段比值模型基于單波段反射率兩兩比值與SSC 構(gòu)建,二元線性模型由兩個(gè)單波段反射率與ln SSC 構(gòu)建,其中模型構(gòu)建所用的單波段反射率為任一波段反射率。單波段模型、波段差值模型與波段比值模型的函數(shù)形式見式(1),二元線性模型的函數(shù)形式見式(2)。
y1 = f(x) (1)
式中:y1 為SSC;x 為單波段反射率或單波段反射率的兩兩差值或單波段反射率的兩兩比值;f 為擬合函數(shù),包括一次多項(xiàng)式、指數(shù)函數(shù)、冪函數(shù)與對數(shù)函數(shù)。
y2 = ax1 + bx2 + c (2)
式中:y2為ln SSC,x1與x2均為單波段反射率,a、b、c 為擬合系數(shù)。
采用決定系數(shù)(R2)評估模型精度,計(jì)算公式為
式中:y^ i為預(yù)測值,yi 為實(shí)測值,y- i 為平均值,n 為樣本總數(shù)。
3 結(jié)果分析
3.1 單波段模型反演精度
不同光譜分辨率的單波段模型的反演精度見圖2。1 nm 光譜分辨率的單波段模型在401 nm 波長處R2最大,為0.23。5 nm 與10 nm 光譜分辨率的模型均在400 nm 波長處R2最大,為0.22。3 種光譜分辨率的單波段模型的R2在400~630 nm 波長范圍內(nèi)整體呈減小趨勢,在700~900 nm 波長范圍內(nèi)呈增大趨勢。原因是極高SSC 水體在可見光波段處反射信號易達(dá)到飽和,而在長波、短波波段處不易飽和,這與周媛等[7]得到的結(jié)論基本一致。整體上,單波段模型的反演效果較差,R2不足0.25。
3.2 波段差值模型反演精度
不同光譜分辨率的波段差值模型的反演精度見圖3,反演精度最高時(shí)模型擬合曲線見圖4[其中子圖(a)橫坐標(biāo)軸標(biāo)值為687 nm 波段與685 nm 波段的反射率差值,其他同理]。分析圖3 可知,波段差值模型的R2大多小于0.40,而由650~720 nm 波段反射率與560~700 nm 波段反射率構(gòu)成的波段差值模型的R2 大于0.40。光譜分辨率為1 nm 由687 nm 波段反射率與685 nm 波段反射率構(gòu)成的波段差值模型(記為R687-R685 模型)的R2最大,為0.76;光譜分辨率為5 nm 由695 nm 波段反射率與685 nm 波段反射率構(gòu)成的波段差值模型(記為R695-R685 模型)的R2 最大,為0.68;光譜分辨率為10 nm 由690 nm 波段反射率與680 nm 波段反射率構(gòu)成的波段差值模型(記為R690-R680 模型)的R2最大,為0.67。
分析圖4 可知,3 種光譜分辨率的R2最大時(shí),模型函數(shù)形式均為冪函數(shù),且SSC 越高反射率差值越小。
3.3 波段比值模型反演精度
不同光譜分辨率的波段比值模型的反演精度見圖5(其中橫坐標(biāo)軸標(biāo)值為比值模型的分母,縱坐標(biāo)軸標(biāo)值為比值模型的分子),反演精度最高時(shí)模型擬合曲線見圖6[其中子圖(a)橫坐標(biāo)軸標(biāo)值為430 nm 波段與411 nm 波段的反射率比值,其他同理]。分析圖5可知,與波段差值模型類似,波段比值模型的R2 大部分也小于0.40。由420~460 nm 或730~900 nm 波段反射率作分子與400~420 nm 波段反射率作分母構(gòu)建的比值模型,以及由650~720 nm 波段反射率作分子與400~700 nm 波段反射率作分母構(gòu)建的比值模型的R2多數(shù)大于0.40。相比波段差值模型,波段比值模型中R2大于0.40 的模型數(shù)量更多,但R2 最大值小于波段差值模型的R2最大值。光譜分辨率為1 nm 由430nm 波段反射率與411 nm 波段反射率構(gòu)成的波段比值模型(記為R430/ R411 模型)的R2最大,為0.67;光譜分辨率為5 nm 由685 nm 波段反射率與660 nm 波段反射率構(gòu)成的波段比值模型(記為R685/ R660 模型)的R2最大,為0.53;光譜分辨率為10 nm 由430 nm 波段反射率與410 nm 波段反射率構(gòu)成的波段比值模型(記為R430/ R410 模型)的R2最大,為0.53。分析圖6可知,波段比值模型的函數(shù)形式均為指數(shù)函數(shù),SSC 越高反射率比值越大。
3.4 二元線性模型反演精度
不同光譜分辨率的二元線性模型的反演精度見圖7(其中橫坐標(biāo)軸代表模型的x1,縱坐標(biāo)軸代表模型的x2),反演精度最高的模型對ln SSC 的預(yù)測值與實(shí)際值對比見圖8。分析圖7 可知,二元線性模型的R2分布與波段比值模型相似,R2 大多數(shù)也不足0.40。由400 ~420 nm波段反射率作x1與400~460 nm 或730~900 nm波段反射率作x2構(gòu)建的模型,以及由400~700 nm 波段反射率作x1與650~720 nm 波段反射率作x2構(gòu)建的模型的R2大于0.40。整體上,二元線性模型的R2大于波段比值模型的。光譜分辨率為1、5、10 nm 時(shí)R2 最大的模型分別為[R411,R434]模型(基于411 nm 波段反射率作x1、434 nm 波段反射率作x2構(gòu)建的二元線性模型,以下同理)、[R665,R670] 模型、[R660,R670]模型,其R2依次為0.69、0.61、0.60。
綜上,波段差值模型、波段比值模型與二元線性模型的R2整體大于單波段模型的,原因是單波段模型易受噪聲影響,差值運(yùn)算、比值運(yùn)算與二元線性組合在一定程度上可以降低噪聲的干擾。針對3 種光譜分辨率數(shù)據(jù),最優(yōu)的SSC 反演模型均為波段差值模型,其R2均在0.65 以上。
4 結(jié)論
基于原位高光譜數(shù)據(jù)與SSC 數(shù)據(jù)構(gòu)建了極高SSC反演模型(包括單波段模型、波段差值模型、波段比值模型、二元線性模型)。研究發(fā)現(xiàn),針對極高含沙量水體的SSC 監(jiān)測,波段差值模型效果相對較好,最大R2為0.76。針對4 類反演模型,1 nm 光譜分辨率的模型在同類中精度最高,最優(yōu)模型的組成波段波長集中在400~440 nm 或560~720 nm。
本文針對2022 年8—11 月監(jiān)測數(shù)據(jù)進(jìn)行建模,在研究區(qū)氣溫、水位、水流條件下模型適用性高,若要實(shí)現(xiàn)模型在全年全周期的普適性,需要采集更加全面的原位高光譜數(shù)據(jù)與河流SSC 值,進(jìn)行模型優(yōu)化。在此基礎(chǔ)上還可以構(gòu)建基于葉綠素a 濃度與透明度等水質(zhì)參數(shù)的監(jiān)測模型,綜合評估極高含沙量水體的水質(zhì)情況。
參考文獻(xiàn):
[1] LAVIGNE H,DOGLIOTTI A,DOXARAN D,et al.The Hyper?maq Dataset:Bio?Optical Properties of Moderately to ExtremelyTurbid Waters[J].Earth System Science Data,2022,14(11):4935-4947.
[2] CUMMINGS V J,BEAUMONT J,MOBILIA V,et al.Responsesof a Common New Zealand Coastal Sponge to Elevated Sus?pended Sediments:Indications of Resilience[J].Mar EnvironRes,2020,155:104886.
[3] OLIVEIRA M L S,DOTTO G L,PINTO D,et al.Nanoparticlesas Vectors of Other Contaminants in Estuarine SuspendedSediments:Natural and Real Conditions[J].Mar Pollut Bull,2021,168:112429.
[4] 段夢偉,李如仁,劉東,等.河流水體懸浮泥沙遙感研究進(jìn)展與展望[J].地球科學(xué)進(jìn)展,2023,38(7):675-687.
[5] 王思夢,秦伯強(qiáng).湖泊水質(zhì)參數(shù)遙感監(jiān)測研究進(jìn)展[J].環(huán)境科學(xué),2023,44(3):1228-1243.
[6] HAN B,LOISEL H,VANTREPOTTE V,et al.Development of aSemi?Analytical Algorithm for the Retrieval of Suspended Par?ticulate Matter from Remote Sensing over Clear to VeryTurbid Waters[J].Remote Sensing,2016,8(3):211.
[7] 周媛,郝艷玲,劉東偉,等.基于Landsat 8 影像的黃河口懸浮物質(zhì)量濃度遙感反演[J].海洋學(xué)研究,2018,36(1):35-45.
[8] 褚言皓,吳文娟,李鵬,等.黃河口懸浮泥沙時(shí)空動態(tài)及其驅(qū)動機(jī)制[J].海洋學(xué)報(bào),2022,44(6):150-163.
[9] KONG J L,SUN X M,WONG D W,et al.A Semi?AnalyticalModel for Remote Sensing Retrieval of Suspended SedimentConcentration in the Gulf of Bohai,China[J].Remote Sensing,2015,7(5):5373-5397.
[10] JIANG B,LIU H L,XING Q G,et al.Evaluating TraditionalEmpirical Models and Bpnn Models in Monitoring the Con?centrations of Chlorophyll?A and Total Suspended Particulateof Eutrophic and Turbid Waters[J].Water,2021,13(5):650.
[11] QIU Z Q,LIU D,DUAN M W,et al.Four?Decades of Sedi?ment Transport Variations in the Yellow River on the LoessPlateau Using Landsat Imagery[J].Remote Sensing of En?vironment,2024,306:114147.
[12] 章敏超,郭碧云.基于GF-1 號的舟山近岸懸浮泥沙濃度反演研究[J].海洋開發(fā)與管理,2018,35(1):126-131.
[13] HOU X J,FENG L,DUAN H T,et al.Fifteen?Year Monitoringof the Turbidity Dynamics in Large Lakes and Reservoirs inthe Middle and Lower Basin of the Yangtze River,China[J].Remote Sensing of Environment,2017,190:107-121.
[14] 唐軍武,田國良,汪小勇,等.水體光譜測量與分析Ⅰ:水面以上測量法[J].遙感學(xué)報(bào),2004,8(1):37-44.
[15] 劉茜,ROSSITER D G.基于高光譜數(shù)據(jù)和MODIS 影像的鄱陽湖懸浮泥沙濃度估算[J].遙感技術(shù)與應(yīng)用,2008,
23(1):7-11,119.[16] SIYOON K,WON S I,SIWAN L.Investigating Mixing Pat?terns of Suspended Sediment in a River Confluence UsingHigh?Resolution Hyperspectral Imagery[J].Journal of Hy?drology,2023,620(PB):129505.
[17] HARRINGMEYER J P,GHOSH N,WEISER M W,et al.A Hy?perspectral View of the Nearshore Mississippi River Delta:Characterizing Suspended Particles in Coastal WetlandsUsing Imaging Spectroscopy[J].Remote Sensing of Envi?ronment,2024,301:113943.
[18] 段麗婷,馬永剛.寧夏60 年泥沙特性和變化規(guī)律分析[J].寧夏工程技術(shù),2021,20(1):34-38.
[19] 艾成,丁環(huán).寧夏清水河流域水文特性分析[J].寧夏農(nóng)林科技,2010,51(3):71-72.
【責(zé)任編輯 栗 銘】
基金項(xiàng)目:山東省創(chuàng)新技術(shù)引導(dǎo)計(jì)劃項(xiàng)目(YDZX2024021);寧夏回族自治區(qū)重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2023BEG02052);中國科學(xué)院重點(diǎn)儀器研發(fā)項(xiàng)目(YJKYYQ20170048)