• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于GF-1號衛(wèi)星WFV數(shù)據(jù)的太湖水質遙感監(jiān)測

    2015-03-09 08:47:48朱利李云梅趙少華郭宇龍
    自然資源遙感 2015年1期
    關鍵詞:懸浮物富營養(yǎng)化太湖

    朱利, 李云梅, 趙少華, 郭宇龍

    (1.環(huán)境保護部衛(wèi)星環(huán)境應用中心,北京 100094;2.國家環(huán)境保護衛(wèi)星遙感重點實驗室,北京 100094;3.南京師范大學地理科學學院教育部虛擬地理環(huán)境重點實驗室,南京 210023)

    基于GF-1號衛(wèi)星WFV數(shù)據(jù)的太湖水質遙感監(jiān)測

    朱利1,2, 李云梅3, 趙少華1,2, 郭宇龍3

    (1.環(huán)境保護部衛(wèi)星環(huán)境應用中心,北京 100094;2.國家環(huán)境保護衛(wèi)星遙感重點實驗室,北京 100094;3.南京師范大學地理科學學院教育部虛擬地理環(huán)境重點實驗室,南京 210023)

    為研究高分一號(GF-1)衛(wèi)星數(shù)據(jù)監(jiān)測太湖水質的可行性,基于新發(fā)射的GF-1衛(wèi)星16 m分辨率的多光譜寬覆蓋(wide field of view,WFV)相機和HJ-1A CCD數(shù)據(jù),對太湖的葉綠素a、懸浮物、透明度和富營養(yǎng)化狀況進行遙感監(jiān)測,以評價GF-1 WFV相機的應用潛力。研究結果表明: GF-1 WFV與HJ-1A CCD數(shù)據(jù)對水質參數(shù)的反演結果具有一致性,可有效反映葉綠素a濃度、懸浮物濃度、透明度和富營養(yǎng)化指數(shù)的空間變化規(guī)律。其中,太湖西北部分布有少量水華藍藻,在大面積藍藻爆發(fā)區(qū)域附近,葉綠素a濃度明顯高于其他區(qū)域的水體,平均濃度為62.46 mg·m-3;懸浮物濃度以竺山灣及西部沿岸湖區(qū)較大,沿西北向東南方向遞減,平均濃度為26.07 mg·L-1;透明度整體從西北向東南遞增,與懸浮物濃度的分布趨勢相反,平均值為22.1 cm;富營養(yǎng)化指數(shù)整體從西北向東南遞減,與葉綠素濃度的分布趨勢相同,平均值為69.62。遙感監(jiān)測指標的結果均符合常規(guī)監(jiān)測規(guī)律。

    葉綠素a濃度;懸浮物濃度;透明度;富營養(yǎng)化指數(shù)

    0 引言

    太湖是我國第三大淡水湖,位于經濟發(fā)達的長江三角洲南緣,屬于典型的富營養(yǎng)化淺水湖泊,是國家確定的“三河三湖”水污染重點治理區(qū)域之一。水體富營養(yǎng)化是太湖水體污染的主要表征之一。葉綠素a濃度和懸浮物濃度等是表征水體富營養(yǎng)化狀況的主要指標。葉綠素a是浮游植物或藻類植物中最豐富的色素,是浮游生物體的重要組成成分之一,其含量的高低與該水體中藻類的種類和數(shù)量密切相關,是反映水體營養(yǎng)狀態(tài)的一個客觀生物學指標,也是水體營養(yǎng)狀態(tài)評價和富營養(yǎng)化評價的重要參數(shù)之一;其含量的大小在一定程度上反映了水體富營養(yǎng)化的程度。懸浮物是指懸浮在水中的固體物質,包括不溶于水的無機物、有機物、泥沙、粘土和微生物等,是我國近岸海水環(huán)境和湖泊水環(huán)境質量檢測的重要參數(shù)之一。懸浮物的存在減少了進入水體的光能量,影響了水生植被的生長,還可以通過輸送磷酸鹽、銨鹽、重金屬和一些致命的細菌等來影響水質;其濃度的大小直接影響著真光層深度、水體透明度及水色等光學性質,從而進一步影響水體的初級生產力和生物量。透明度是描述水體光學特性的基本參數(shù)之一,與水體中的懸浮物、葉綠素、黃色物質的含量和成分密切相關;利用透明度可以間接地反演漫衰減系數(shù),進而獲取真光層深度,進行水下廣場分布的研究。水體富營養(yǎng)化是指湖泊等水體中氮和磷等營養(yǎng)鹽含量過高,導致藻類以及其他水生生物過量繁殖,使得水體透明度下降,溶解氧含量降低,造成湖泊水質惡化從而使生態(tài)系統(tǒng)和水體功能受到損害和破壞,嚴重時產生水華現(xiàn)象。衛(wèi)星遙感監(jiān)測具有大面積、實時性和經濟性等優(yōu)勢,是水質和富營養(yǎng)化監(jiān)測的重要技術手段。通過衛(wèi)星遙感獲取一些常規(guī)監(jiān)測方法難以揭示的污染源和污染物遷移特征,可為污染源的分布和污染物質的擴散規(guī)律等研究提供基礎數(shù)據(jù)。

    2013年4月26日,我國自主研制的第一顆民用高分辨率衛(wèi)星高分一號(GF-1)成功發(fā)射,為我國水環(huán)境的遙感監(jiān)測帶來了巨大機遇。該衛(wèi)星攜載了2 m分辨率全色、8 m分辨率多光譜CCD相機和16 m分辨率多光譜CCD相機,其中,2 m分辨率全色、8 m分辨率多光譜CCD數(shù)據(jù)的幅寬優(yōu)于60 km,由2臺相機拼接而成;16 m分辨率多光譜CCD數(shù)據(jù)幅寬優(yōu)于800 km,由4臺相機拼接而成。GF-1衛(wèi)星不僅可以用于精細探測局部的環(huán)境污染狀況,還可用于水環(huán)境、大氣環(huán)境和生態(tài)環(huán)境質量等大范圍的宏觀監(jiān)測與評價。作為新發(fā)射的衛(wèi)星,并且是我國高分辨率衛(wèi)星系列的首發(fā)星,開展GF-1衛(wèi)星載荷應用的在軌測試工作非常重要,環(huán)保部作為GF-1衛(wèi)星的主要用戶,承擔了其環(huán)境應用的在軌測試工作。

    對于水污染嚴重的湖泊,國內外已有不少關于水污染的遙感監(jiān)測研究,例如: 呂恒等[1]總結了國內外用于湖泊水質參數(shù)反演的3種常用方法,并分析了3種模型的優(yōu)缺點;馬榮華等[2]系統(tǒng)分析了湖泊水色遙感的發(fā)展現(xiàn)狀,表明湖泊水體物質組分的復雜性以及衛(wèi)星傳感器與實際需求的矛盾決定了湖泊水色遙感的難度;李俊生等[3]應用緊密型高分辨率光譜儀(compact high resolution imaging spectrometer,CHRIS)圖像反演了太湖梅梁灣的葉綠素濃度分布圖,取得了較好的結果;閻福禮等[4]對比了環(huán)境小衛(wèi)星和Hyperion高光譜監(jiān)測太湖水質參數(shù)的反演精度,展示了星載高光譜監(jiān)測湖泊水體水質的應用潛力;王橋等[5]介紹了環(huán)境一號衛(wèi)星在水環(huán)境、大氣環(huán)境和生態(tài)遙感監(jiān)測與評價中的應用方法和范例;王彥飛等[6]針對HJ-1A HSI數(shù)據(jù)在巢湖水質監(jiān)測應用方面的適宜性進行了評價,表明HSI水體圖像在530~900 nm波段范圍內數(shù)據(jù)質量較為真實可靠;潘邦龍等[7]利用HJ-1A HIS高光譜數(shù)據(jù),采用回歸克里格方法建立了總氮濃度的定量反演模型;朱利等[8]基于環(huán)境一號衛(wèi)星多光譜數(shù)據(jù),建立了分地區(qū)、分季節(jié)的經驗模型反演葉綠素a濃度;徐雯佳等[9]利用MODIS的第一波段數(shù)據(jù)建立了河北省海域的葉綠素a濃度的遙感模型,取得了較好的效果;陳蕾等[10]提出一種基于水質類型的TM圖像水體信息提取方法,可以快速、準確、有效地提取TM圖像中的水體分布信息。但由于GF-1衛(wèi)星是新發(fā)射衛(wèi)星,目前還沒有相關應用的研究報道,所以本文根據(jù)其在軌測試的結果,評價GF-1衛(wèi)星在水環(huán)境遙感監(jiān)測方面的應用潛力,以期為后續(xù)GF-1衛(wèi)星的水環(huán)境遙感監(jiān)測提供參考。因太湖面積較大,故本文主要探討GF-1寬覆蓋(wide field of view,WFV)相機對太湖的葉綠素a濃度、懸浮物濃度、透明度和富營養(yǎng)化狀況的遙感監(jiān)測,并分析上述水質參數(shù)的空間分布特征,進而綜合評價GF-1 WFV相機對水質參數(shù)的監(jiān)測能力。

    1 研究區(qū)概況與數(shù)據(jù)源

    1.1 研究區(qū)概況

    本文選擇太湖作為研究區(qū),其地理范圍在E119°55′~120°40′,N30°55′ ~31°30′之間;湖水面積約2 338 km2,平均水深約2 m;水體富營養(yǎng)化較嚴重。

    1.2 遙感圖像預處理

    使用2013年8月9日獲取的GF-1衛(wèi)星16 m分辨率多光譜CCD圖像,傳感器為WFV2。GF-1 WFV2圖像的預處理主要包括輻射定標、幾何糾正和大氣校正。

    1)輻射定標利用ENVI軟件,將GF-1 WFV2圖像灰度值Le通過絕對定標系數(shù)和輻射定標公式進行定標,即

    Le=Gain·DN+Bias,

    (1)

    式中:Gain為波段增益,即定標斜率,W·m-2·sr-1·μm-1;DN為衛(wèi)星載荷觀測值;Bias為波段偏置值,即定標截距,W·m-2·sr-1·μm-1。GF-1號衛(wèi)星WFV圖像的4個波段的增益和偏置值見表1。

    表1 GF-1號衛(wèi)星WFV相機絕對輻射定標系數(shù)Tab.1 Absolute radiation calibrated coefficients of WFV cameras of GF-1 satellite

    2)幾何糾正利用ERDAS軟件中的Georeferencing模塊,以研究區(qū)經過幾何精糾正的TM圖像為基準,通過在研究區(qū)選取控制點完成對GF-1 WFV2圖像的幾何精糾正。

    3)水色大氣校正則利用太湖站CE318實測數(shù)據(jù),輸入大氣輻射傳輸模型6S,獲得大氣氣溶膠散射和瑞利散射,得到各個波段圖像的反射率[6,9]。

    4)由于水質參數(shù)反演模型只針對水體進行計算,因而要根據(jù)水陸分布圖對預處理后的GF-1 WFV2圖像數(shù)據(jù)進行掩模處理,得到水體圖像。

    1.3 野外光譜實測

    太湖野外實測光譜的采集使用ASD公司生產的ASD FieldSpec Pro便攜式光譜輻射計,其波譜范圍為350~1 050 nm。為減少水體鏡面反射和船體自身陰影的影響,測量采用唐軍武等[11]提出的陸上Ⅱ類水體水面以上光譜測量的方法。

    2 研究方法

    目前,利用遙感器觀測得到的光譜數(shù)據(jù)反演水質參數(shù)通常有3種方法,即經驗方法、半經驗-半分析方法和分析方法。經驗方法根據(jù)統(tǒng)計分析理論,在大量現(xiàn)場實測數(shù)據(jù)的基礎上建立遙感反射率與水質參數(shù)濃度之間的統(tǒng)計關系,進而外推水體水質參數(shù)濃度,常用的模型主要有單波段模型、一階微分模型、波段比值模型和多波段組合模型等;半經驗-半分析方法借助于對水體主要水色要素固有光學特性的研究,采用分析和經驗模式獲得各種水體成分的吸收和后向散射系數(shù),進而估算水質參數(shù)濃度,常用的模型主要有3波段模型和4波段模型等;分析方法則利用生物光學模型描述水體組分與輻照度比之間的關系,并模擬水體的輻照度比,水體中各組分用其單位吸收系數(shù)和單位散射系數(shù)表示,用輻射傳輸模型模擬光在大氣和水體中的傳播過程,然后根據(jù)遙感數(shù)據(jù)反演水質參數(shù)濃度。此外,一些新型的經驗方法,如神經網(wǎng)絡、偏最小二乘及支持向量機等也被引入到水色遙感中,用以反演水質參數(shù)濃度。由于沒有同步觀測的數(shù)據(jù)和缺少散射系數(shù)、后向散射系數(shù)等固有光學數(shù)據(jù),暫無法建立以生物光學模型為基礎的分析模型,因此本研究的水質遙感監(jiān)測模型主要以經驗模型和半經驗-半分析模型為主。

    盡管目前尚缺乏與GF-1衛(wèi)星過頂同步野外采樣數(shù)據(jù),難以根據(jù)實測數(shù)據(jù)建立估算模型;但考慮到GF-1數(shù)據(jù)與HJ-1A CCD數(shù)據(jù)有相似的波段設置,且目前已經有許多學者利用HJ-1A CCD數(shù)據(jù)進行了水色參數(shù)估算研究[8,11,12-15],并取得了較好的效果;因此,本文首先根據(jù)不同傳感器的波段響應函數(shù)與野外ASD實測的太湖大量高光譜數(shù)據(jù),對GF-1 WFV與HJ-1A CCD圖像的波段設置進行對比。

    有效波長范圍內的等效反射率計算方法為

    (2)

    式中:f(λ)為光譜響應函數(shù);λ為觀測波長;λmin和λmax分別為波段范圍最小和最大波長;Rb為b波段范圍內的等效反射率。

    進而參考前人的研究[8,12],選擇葉綠素a、懸浮物、透明度和富營養(yǎng)化指數(shù)的估算模型分別計算相應的水質參數(shù)。

    1)葉綠素a濃度的反演模型為

    Cchl-a=4.089(b4/b3)2-0.746(b4/b3)+29.733 ,

    (3)

    式中:Cchl-a為葉綠素a的濃度,mg·m-3;b3和b4分別為GF-1 WFV2圖像經過輻射校正和大氣校正后的第3波段和第4波段圖像像元亮度值。

    2)懸浮物濃度的反演模型為

    CTSS=119.62(b3/b2)6.082 3,

    (4)

    式中:CTSS為總懸浮物濃度,mg·L-1;b2和b3分別為GF-1 WFV2圖像經過輻射校正和大氣校正后的第2波段和第3波段圖像像元亮度值。

    3)透明度的反演模型為

    (5)

    式中Zsd為透明度,cm。前人研究表明,太湖區(qū)域的透明度主要受懸浮物影響[14]。

    4)富營養(yǎng)化指數(shù)的估算模型為

    TLI=0.54TLIchl-a+0.46TLITSS,

    (6)

    TLIchl-a=10(2.5+1.086lnCchl-a),

    (7)

    TLITSS=10(5.118+0.49lnCTSS),

    (8)

    式(6)—(8)中:TLI為綜合遙感富營養(yǎng)化指數(shù),無量綱,可用于對湖泊營養(yǎng)狀態(tài)進行連續(xù)數(shù)值分級評價;TLIchl-a為葉綠素富營養(yǎng)化指數(shù);TLITSS為懸浮物富營養(yǎng)化指數(shù)。

    3 結果與討論

    3.1 遙感圖像的一致性分析

    GF-1 WFV圖像與HJ-1A CCD圖像的波段響應函數(shù)如圖1所示,其中虛線和實線分別表示HJ-1A CCD圖像和GF-1 WFV圖像4個波段的波段響應函數(shù)。

    圖1 GF-1 WFV與HJ-1A CCD波譜設置對比Fig.1 Comparison of spectrum settings between sensers of GF-1 WFV and HJ-1A CCD

    從圖1可以看出: ① 2種傳感器4個波段的有效波譜覆蓋范圍比較一致;② 2種傳感器不同設置中的強響應波段位置較為一致。說明GF-1 WFV和HJ-1A CCD這2種傳感器的理論波段設置比較相似。

    為了進一步對比2種傳感器的波段設置差異,利用在太湖野外實測得到的345組高光譜反射率,分別對2種傳感器進行模擬,并逐波段對比2種傳感器的模擬反射率,結果如圖2所示。

    (a) 第1波段(b) 第2波段(c) 第3波段(d) 第4波段

    圖2 GF-1 WFV與HJ-1A CCD的4個波段模擬反射率

    Fig.2 Simulated reflectivities of 4 bands of GF-1 WFV and HJ-1A CCD

    從圖2可以看出,在4個波段上,2種傳感器都呈現(xiàn)出很高的一致性: 線性擬合的決定系數(shù)都在0.99以上,斜率都在1±0.04之間,截距也都接近0。

    上述分析表明,GF-1 WFV圖像與HJ-1A CCD圖像在波段設置上具有很強的一致性,因而在理論上可以采取通用的估算模型進行水質參數(shù)的計算。本文選擇2013年8月9日在太湖過頂?shù)腉F-1 WFV與HJ-1A CCD圖像數(shù)據(jù)(表2),二者過頂時間相差46 min,水質參數(shù)的空間變化不是太大。

    表2 同步獲取的GF-1 WFV和HJ-1A CCD參數(shù)對比Tab.2 Comparison of parameters between GF-1 WFV and HJ-1A CCD acquired synchronously

    3.2 葉綠素a濃度的反演結果

    利用上述數(shù)據(jù)和方法,得出太湖葉綠素a濃度的空間分布情況圖(圖3)和統(tǒng)計表(表3)。

    (a) GF-1 WFV2 (b) HJ-1A CCD1

    圖3 太湖葉綠素a濃度(mg·m-3)分布圖

    Fig.3 Distribution of chlorophyll-a(mg·m-3) in Taihu Lake

    表3 太湖葉綠素a濃度統(tǒng)計值Tab.3 Statistical value of chlorophyll-a concentration in Taihu Lake (mg·cm-3)

    圖3中的黃色區(qū)域為藍藻水華分布區(qū)域,基于歸一化植被指數(shù)(NDVI),結合灰度直方圖確定閾值來區(qū)分水華和正常水體區(qū)域(由于篇幅所限,本文不展開討論)。進行水質參數(shù)反演時,將對水華分布區(qū)域做掩模,并和陸地信息一并去除,使二者不參與水質參數(shù)反演。為了間接驗證反演結果的正確性,對同步獲取的HJ-1A CCD數(shù)據(jù)進行同樣的處理,并對比2組數(shù)據(jù)的反演結果。

    由圖3(a)可以看出,2013年8月9日太湖葉綠素a濃度的空間分布平緩,高值區(qū)分布在竺山灣、梅梁灣大部分湖區(qū)以及西部湖區(qū),濃度值集中在60~110 mg·m-3之間;東部湖區(qū)及貢湖灣大部分區(qū)域葉綠素a濃度值較低,在20~40 mg·m-3之間。低值區(qū)面積較小,零星分布于太湖南部,葉綠素a濃度在20 mg·m-3以下。

    從表3可以看出,由GF-1 WFV反演的全湖區(qū)葉綠素a濃度變化范圍較大。同時,2種傳感器反演得到的葉綠素濃度的最大值和最小值的差距較大,這可能是由于不同傳感器的觀測角度、大氣狀況及衛(wèi)星過頂時間的差異所導致;但從均值和空間分布情況來看,2種傳感器的監(jiān)測結果具有較好的一致性。

    3.3 懸浮物濃度的反演結果

    利用上述數(shù)據(jù)和方法,得出太湖懸浮物濃度的空間分布情況圖(圖4)和統(tǒng)計表(表4)。

    (a) GF-1 WFV2 (b) HJ-1A CCD1

    圖4 太湖懸浮物濃度(mg·L-1)分布圖

    Fig.4 Distribution of suspended solids concentration(mg·L-1) in Taihu Lake

    表4 太湖懸浮物濃度統(tǒng)計值Tab.4 Statistical value of suspended solids concentration in Taihu Lake (mg·L-1)

    從圖4(a)可以看出,2013年8月9日太湖懸浮物濃度的空間變化劇烈,高值區(qū)集中在梅梁灣沿岸湖區(qū)以及太湖西部入湖河流河口附近,懸浮物濃度在110 mg·L-1以上;從貢湖灣南到太湖南部沿岸為低值區(qū),濃度都在55 mg·L-1以下,甚至有部分區(qū)域小于25 mg·L-1。

    從表4可以看出,由GF-1 WFV反演的全湖區(qū)懸浮物濃度變化范圍較大。與葉綠素濃度類似,在懸浮物濃度的統(tǒng)計值中,2種傳感器的最大值和最小值也呈現(xiàn)出一定的偏差,但均值仍然比較接近,體現(xiàn)出較好的一致性。

    3.4 透明度的反演結果

    利用上述數(shù)據(jù)和方法,得出太湖透明度的空間分布情況圖(圖5)和統(tǒng)計表(表5)。

    (a) GF-1 WFV2 (b) HJ-1A CCD1

    圖5 太湖透明度(cm)分布圖

    Fig.5 Distribution of transparency(cm) in Taihu Lake

    表5 太湖透明度統(tǒng)計值Tab.5 Statistical value of transparency in Taihu Lake (cm)

    從圖5(a)可以看出,2013年8月9日太湖透明度的分布情況較為簡單,與懸浮物濃度分布有相反的趨勢,在梅梁灣及太湖西部湖區(qū),透明度值較低,都在18 cm以下;而從貢湖灣延伸到太湖南部,這些區(qū)域透明度值都較高,在22 cm以上。

    從表5可以看出,全湖最大透明度為62.19 cm,最小透明度僅為3.6 cm,均值為22.17 cm。

    3.5 富營養(yǎng)化狀況的反演結果

    利用上述數(shù)據(jù)和方法,得出太湖富營養(yǎng)化的空間分布情況圖(圖6)和統(tǒng)計表(表6)。

    (a) GF-1 WFV2 (b) HJ-1A CCD1

    圖6 太湖富營養(yǎng)化指數(shù)分布圖

    Fig.6 Distribution of eutrophication index in Taihu Lake

    表6 太湖富營養(yǎng)化指數(shù)統(tǒng)計值Tab.6 Statistical value of eutrophication index in Taihu Lake

    從圖6(a)可以看出,2013年8月9日太湖富營養(yǎng)化狀況以梅梁灣和西北沿岸湖區(qū)最為嚴重,富營養(yǎng)化指數(shù)達到69以上,湖心區(qū)也有部分高值區(qū);其次為西部大面積湖區(qū),富營養(yǎng)化指數(shù)在69左右;南部湖區(qū)和東部湖區(qū)富營養(yǎng)化水平相對較低,富營養(yǎng)化指數(shù)在65左右。

    從表6可以看出,由GF-1 WFV反演的全湖區(qū)最大富營養(yǎng)化指數(shù)為87.90,最小富營養(yǎng)化指數(shù)為49.45,均值為69.62。與其余3個指標相比,2種傳感器估算得到的富營養(yǎng)化狀態(tài)在3項統(tǒng)計指標上都體現(xiàn)出較好的一致性。

    4 結論

    本文利用GF-1號衛(wèi)星16 m空間分辨率多光譜寬覆蓋(WFV)數(shù)據(jù),對太湖的葉綠素a、懸浮物、透明度和富營養(yǎng)化狀況開展在軌測試遙感監(jiān)測,進而評價GF-1號衛(wèi)星WFV相機的應用潛力,得出如下結論:

    1)在軌測試的初步評價表明,GF-1號衛(wèi)星WFV相機獲取的圖像可以反映葉綠素a濃度、懸浮物濃度、透明度和富營養(yǎng)化指數(shù)的空間變化規(guī)律,能夠有效地用于水質等環(huán)境監(jiān)測。

    2)通過GF-1號衛(wèi)星的監(jiān)測發(fā)現(xiàn),太湖西北部分布有少量水華藍藻;在大面積藍藻爆發(fā)區(qū)域附近,葉綠素a濃度明顯高于其他區(qū)域的水體。懸浮物濃度以竺山灣及西部沿岸湖區(qū)較大,沿西北向東南方向遞減。透明度在整體上從西北到東南遞增,與懸浮物濃度的分布趨勢相反,符合常規(guī)監(jiān)測規(guī)律。富營養(yǎng)化指數(shù)在整體上從西北到東南遞減,與葉綠素濃度的分布趨勢相同,也符合常規(guī)監(jiān)測規(guī)律。

    3)實驗表明,一些針對HJ-1A CCD圖像數(shù)據(jù)建立的水色要素估算模型,也能較好地應用在GF-1號衛(wèi)星16 m分辨率的圖像上。對比GF-1 WFV傳感器與HJ-1A CCD傳感器,發(fā)現(xiàn)二者具有較為一致的波段設置,且可以得到近似的水色要素監(jiān)測結果,間接表明GF-1 WFV傳感器數(shù)據(jù)反演水質參數(shù)結果的可靠性。

    4)由于在軌測試時間緊張,覆蓋太湖研究區(qū)的無云或少云圖像很少,并缺乏同步實測數(shù)據(jù),無法客觀地驗證有關模型的精度。今后,一方面要采用更多的遙感圖像和開展同步的湖面觀測試驗;另一方面可以考慮采用普適性更強的半分析模型進行葉綠素a、懸浮物濃度和透明度的反演,或引進更多技術指標進行湖水富營養(yǎng)化狀態(tài)的評價。

    [1] 呂恒,江南,李新國.內陸湖泊的水質遙感監(jiān)測研究[J].地球科學進展,2005,20(2):185-192. Lü H,Jiang N,Li X G.The study on water quality of inland lake monitoring by remote sensing[J].Advances in Earth Science,2005,20(2):185-192.

    [2] 馬榮華,唐軍武,段洪濤,等.湖泊水色遙感研究進展[J].湖泊科學,2009,21(2):143-158. Ma R H,Tang J W,Duan H T,et al.Progress in lake water color remote sensing[J].Journal of Lake Sciences,2009,21(2):143-158.

    [3] 李俊生,張兵,申茜,等.航天成像光譜儀CHRIS在內陸水質監(jiān)測中的應用[J].遙感技術與應用,2007,22(5):593-597. Li J S,Zhang B,Shen Q.Application of spaceborne imaging spectrometry CHRIS in monitoring of inland water quality[J].Remote Sensing Technology and Application,2007,22(5):593-597.

    [4] 閻福禮,王世新,周藝,等.利用Hyperion星載高光譜傳感器監(jiān)測太湖水質的研究[J].紅外與毫米波學報,2006,25(6):460-464. Yan F L,Wang S X,Zhou Y,et al.Monitoring the water quality of Taihu Lake by using Hyperion hyperspectral data[J].Journal of Infrared and Millimeter Waves,2006,25(6):460-464.

    [5] 王橋,魏斌,王昌佐,等.基于環(huán)境一號衛(wèi)星的生態(tài)環(huán)境遙感監(jiān)測[M].北京:科學出版社,2010. Wang Q,Wei B,Wang C Z,et al.Environmental Monitoring Based on HJ-1 Satellite[M].Beijing:Science Press,2010.

    [6] 王彥飛,李云梅,呂恒,等.環(huán)境一號衛(wèi)星高光譜遙感數(shù)據(jù)的內陸水質監(jiān)測適宜性——以巢湖為例[J].湖泊科學,2011,23(5):789-795. Wang Y F,Li Y M,Lü H,et al.Suitability assessment of lake water quality monitoring on waterbody images acquired by HJ-1A hyperspectral imager:A case study of Lake Chaohu[J].Journal of Lake Sciences,2011,23(5):789-795.

    [7] 潘邦龍,易維寧,王先華,等.基于環(huán)境一號衛(wèi)星超光譜數(shù)據(jù)的多元回歸克里格模型反演湖泊總氮濃度的研究[J].光譜學與光譜分析,2011,31(7):1884-1888. Pan B L,Yi W N,Wang X H,et al.Inversion of the lake total nitrogen concentration by multiple regression Kriging model based on hyperspectral data of HJ-1A[J].Spectroscopy and Spectral Analysis,2011,31(7):1884-1888.

    [8] 朱利,姚延娟,吳傳慶,等.基于環(huán)境一號衛(wèi)星的內陸水體水質多光譜遙感監(jiān)測[J].地理與地理信息科學,2010,26(2):81-84. Zhu L,Yao Y J,Wu C Q,et al.Remote sensing monitoring on water quality of interior water body based on the HJ-1 satellite multi-spectral data[J].Geography and Geo-Information Science,2010,26(2):81-84.

    [9] 徐雯佳,楊斌,田力,等.應用MODIS數(shù)據(jù)反演河北省海域葉綠素a濃度[J].國土資源遙感,2012,24(4):152-156. Xu W J,Yang B,Tian L,et al.Retrieval of chlorophyll-a concentration by using MODIS data in Hebei sea area[J].Remote Sensing for Land and Resources,2012,24(4):152-156.

    [10]陳蕾,鄧孺孺,陳啟東,等.基于水質類型的TM圖像水體信息提取[J].國土資源遙感,2012,24(1):90-94. Chen L,Deng R R,Chen Q D,et al.The extraction of water body information from TM imagery based on water quality types[J].Remote Sensing for Land and Resources,2012,24(1):90-94.

    [11]唐軍武,田國良,汪小勇,等.水體光譜測量與分析I:水面以上測量法[J].遙感學報,2004,8(1):37-44. Tang J W,Tian G L,Wang X Y,et al.The methods of water spectra measurement and analysis I:Above-water method[J].Journal of Remote Sensing,2004,8(1):37-44.

    [12]徐祎凡,李云梅,王橋,等.基于環(huán)境一號衛(wèi)星多光譜影像數(shù)據(jù)的三湖一庫富營養(yǎng)化狀態(tài)評價[J].環(huán)境科學學報,2011,31(1):81-93. Xu Y F,Li Y M,Wang Q,et al.Eutrophication evaluation of three lakes and one reservoir using CCD images from the HJ-1 satellite[J].Acta Scientiae Circumstantiae,2011,31(1):81-93.

    [13]郭宇龍,李云梅,朱利,等.基于HJ-1A CCD數(shù)據(jù)的高光譜影像重構研究[J].環(huán)境科學,2013,34(1):69-76. Guo Y L,Li Y M,Zhu L,et al.Research of hyperspectral reconstruction based on HJ-1A CCD data[J].Environmental Science,2013,34(1):69-76.

    [14]肖艷芳,趙文吉,朱琳.近海水體懸浮物HJ-1號小衛(wèi)星CCD定量反演研究[J].海洋科學,2012,36(8):59-63. Xiao Y F,Zhao W J,Zhu L.Quantitative retrieval model of suspended sediment concentration in estuary based on HJ-1 CCD image[J].Marine Sciences,2012,36(8):59-63.

    [15]曠達,韓秀珍,劉翔,等.基于環(huán)境一號衛(wèi)星的太湖葉綠素a濃度提取[J].中國環(huán)境科學,2010,30(9):1268-1273. Kuang D,Han X Z,Liu X,et al.Quantitative estimation of Taihu chlorophyll-a concentration using HJ-1A and 1B CCD imagery[J].China Environmental Science,2010,30(9):1268-1273.

    (責任編輯: 邢宇)

    Remote sensing monitoring of Taihu Lake water quality by using GF-1 satellite WFV data

    ZHU Li1,2, LI Yunmei3, ZHAO Shaohua1,2, GUO Yulong3

    (1.EnvironmentalSatelliteCenter,MinistryofEnvironmentalProtection,Beijing100094,China;2.StateEnvironmentalProtectionKeyLaboratoryofSatelliteRemoteSensing,Beijing100094,China;3.SchoolofGeographyScience,TheKeyLaboratoryofVirtualGeographicEnvironment,MinistryofEducation,NanjingNormalUniversity,Nanjing210023,China)

    The chlorophyll-a concentration, suspended solids concentration, transparency and eutrophication index for the Taihu Lake were obtained by using GF-1 satellite wide field of view (WFV) camera and HJ-1A satellite CCD camera data with the purpose of evaluating the GF-1 satellite’s application capability in environmental monitoring. The results show that the WFV camera data, consistent with the HJ-1A inversion results, could be used to monitor the distribution of the water quality parameters including chlorophyll-a concentration, suspended solids concentration, transparency and eutrophication index. It was shown that the algal blooming areas were found in the northwest of the Taihu Lake, the chlorophyll-a concentration was the highest near the algal blooming area, and the average of chlorophyll-a concentration in the whole lake was 62.46 mg·m-3. The suspended solids concentration was higher near the Zhusan Valley and the west coast area and decreased in the NW-SE direction, with the whole average value being 26.07mg·L-1. The transparency decreased in the NW-SE direction in contrast to the distribution of suspended solids concentration, with the whole average value being 22.1 cm. The eutrophication index decreased southeastward, being the same as the distribution of chlorophyll-a concentration, with the whole average value being 69.62. These results are consistent with the results of the routine surveying.

    chlorophyll-a concentration;suspended solids concentration;transparency;eutrophication index

    2013-11-11;

    2014-01-13

    國家自然科學基金項目“基于四波段CCD數(shù)據(jù)的藍藻水華遙感識別指數(shù)構建與藻密度定量反演研究”(編號: 41101378)和國家高分辨率對地觀測重大專項“環(huán)境保護遙感動態(tài)監(jiān)測信息服務系統(tǒng)”(編號: E05-Y30B02-9001-13/15)共同資助。

    10.6046/gtzyyg.2015.01.18

    朱利,李云梅,趙少華,等.基于GF-1號衛(wèi)星WFV數(shù)據(jù)的太湖水質遙感監(jiān)測[J].國土資源遙感,2015,27(1):113-120.(Zhu L,Li Y M,Zhao S H,et al.Remote sensing monitoring of Taihu Lake water quality by using GF-1 satellite WFV data[J].Remote Sensing for Land and Resources,2015,27(1):113-120.)

    TP 79; X 87

    A

    1001-070X(2015)01-0113-08

    朱利(1978-),男,博士,高級工程師,主要從事水環(huán)境定量遙感研究。Email: zlpinan@163.com。

    郭宇龍(1988-),男,博士,主要從事水環(huán)境定量遙感研究。Email: gyl.18@163.com。

    猜你喜歡
    懸浮物富營養(yǎng)化太湖
    懸浮物對水質氨氮測定的影響
    化工管理(2022年14期)2022-12-02 11:43:52
    基于臨界點的杭州灣水體富營養(yǎng)化多年變化研究
    海洋通報(2022年4期)2022-10-10 07:41:48
    洪口水庫近年富營養(yǎng)化程度時間分布的研究
    太湖思變2017
    玩具世界(2017年4期)2017-07-21 13:27:24
    太湖攬春
    寶藏(2017年2期)2017-03-20 13:16:42
    洞庭湖典型垸內溝渠水體富營養(yǎng)化評價
    太湖
    中亞信息(2016年3期)2016-12-01 06:08:24
    富營養(yǎng)化藻的特性與水熱液化成油的研究
    廢水懸浮物國標檢測方法的修改建議
    壓裂返排液中懸浮物去除的室內研究
    久久鲁丝午夜福利片| 免费观看a级毛片全部| 97超视频在线观看视频| 2022亚洲国产成人精品| 99九九线精品视频在线观看视频| 成人性生交大片免费视频hd| 欧美成人午夜免费资源| 看十八女毛片水多多多| www.色视频.com| 日日干狠狠操夜夜爽| 狂野欧美白嫩少妇大欣赏| 别揉我奶头 嗯啊视频| 国产成人freesex在线| 午夜福利在线在线| 一级二级三级毛片免费看| 午夜精品在线福利| 天堂网av新在线| 成人亚洲欧美一区二区av| av播播在线观看一区| 日韩强制内射视频| 麻豆精品久久久久久蜜桃| 亚洲国产精品国产精品| 欧美日韩综合久久久久久| 国产一区二区在线观看日韩| av在线播放精品| 亚洲成人一二三区av| 国产精品不卡视频一区二区| 亚洲av国产av综合av卡| 国产精品女同一区二区软件| 韩国av在线不卡| 日韩一区二区视频免费看| 精品不卡国产一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 99久久精品热视频| 国模一区二区三区四区视频| 亚洲人成网站高清观看| 国产精品一区二区在线观看99 | 男人舔女人下体高潮全视频| 天堂av国产一区二区熟女人妻| 久久亚洲国产成人精品v| 18+在线观看网站| 尾随美女入室| 插逼视频在线观看| 91aial.com中文字幕在线观看| 美女主播在线视频| 国产成年人精品一区二区| 久久久精品欧美日韩精品| 欧美xxxx黑人xx丫x性爽| 最后的刺客免费高清国语| 日本黄色片子视频| 极品教师在线视频| 天堂影院成人在线观看| 亚洲精华国产精华液的使用体验| 高清午夜精品一区二区三区| 免费观看a级毛片全部| 精品午夜福利在线看| 69人妻影院| 一区二区三区乱码不卡18| 人妻少妇偷人精品九色| 乱系列少妇在线播放| 免费看日本二区| 亚洲天堂国产精品一区在线| 欧美日本视频| 少妇熟女aⅴ在线视频| 国产伦一二天堂av在线观看| 久久久久免费精品人妻一区二区| 一级毛片aaaaaa免费看小| 国产大屁股一区二区在线视频| 别揉我奶头 嗯啊视频| 国产大屁股一区二区在线视频| 免费黄频网站在线观看国产| 乱码一卡2卡4卡精品| 成人综合一区亚洲| 亚洲国产日韩欧美精品在线观看| 黄色一级大片看看| 色网站视频免费| 在线观看av片永久免费下载| 久久国产乱子免费精品| 最近最新中文字幕免费大全7| 亚洲国产最新在线播放| 99久久精品一区二区三区| 亚洲精品aⅴ在线观看| 中文精品一卡2卡3卡4更新| 亚洲av免费高清在线观看| 亚洲在线观看片| 国产成人精品久久久久久| 五月玫瑰六月丁香| 水蜜桃什么品种好| 久久鲁丝午夜福利片| 日韩不卡一区二区三区视频在线| 亚洲精品国产av成人精品| 国产综合懂色| 神马国产精品三级电影在线观看| 两个人视频免费观看高清| 亚洲人成网站在线播| 人人妻人人澡人人爽人人夜夜 | 搡老妇女老女人老熟妇| 丝袜喷水一区| 国产精品久久久久久久电影| 亚洲熟妇中文字幕五十中出| 简卡轻食公司| 免费看不卡的av| 日本与韩国留学比较| 少妇裸体淫交视频免费看高清| 免费黄色在线免费观看| 免费av不卡在线播放| 久久国产乱子免费精品| 日本一二三区视频观看| 亚洲av中文字字幕乱码综合| 久久久午夜欧美精品| 欧美性猛交╳xxx乱大交人| 又爽又黄无遮挡网站| 伊人久久精品亚洲午夜| 亚洲综合色惰| 韩国av在线不卡| 特级一级黄色大片| 国产伦精品一区二区三区视频9| 18禁裸乳无遮挡免费网站照片| 最近最新中文字幕大全电影3| 美女脱内裤让男人舔精品视频| 国产一区亚洲一区在线观看| 在线免费观看不下载黄p国产| 亚洲欧美成人精品一区二区| 美女主播在线视频| 午夜福利视频1000在线观看| 日韩在线高清观看一区二区三区| 免费播放大片免费观看视频在线观看| 亚洲欧美清纯卡通| kizo精华| 免费av毛片视频| 久久久色成人| 高清视频免费观看一区二区 | 观看免费一级毛片| 婷婷色麻豆天堂久久| 一级毛片久久久久久久久女| 国产乱人偷精品视频| 亚洲欧美一区二区三区黑人 | 亚洲国产日韩欧美精品在线观看| 久久久久久久久大av| 久久久色成人| 久久鲁丝午夜福利片| 国产成人免费观看mmmm| 非洲黑人性xxxx精品又粗又长| 不卡视频在线观看欧美| 亚洲美女视频黄频| 国产美女午夜福利| 伦精品一区二区三区| 草草在线视频免费看| 男插女下体视频免费在线播放| 精品国内亚洲2022精品成人| 精品久久久久久久久亚洲| 精品久久久久久久末码| 黄色欧美视频在线观看| 欧美bdsm另类| 水蜜桃什么品种好| 亚洲经典国产精华液单| 国产亚洲最大av| 久久草成人影院| 精品国产露脸久久av麻豆 | 国产在线一区二区三区精| 国内精品美女久久久久久| 狂野欧美激情性xxxx在线观看| 亚洲精品456在线播放app| 国产男女超爽视频在线观看| 淫秽高清视频在线观看| 乱码一卡2卡4卡精品| 久久久久久九九精品二区国产| 日本免费a在线| 亚洲久久久久久中文字幕| 身体一侧抽搐| 69av精品久久久久久| 丝袜喷水一区| 少妇人妻一区二区三区视频| 能在线免费看毛片的网站| 久久久精品免费免费高清| 日韩中字成人| 免费黄色在线免费观看| 欧美日韩一区二区视频在线观看视频在线 | 小蜜桃在线观看免费完整版高清| 热99在线观看视频| 午夜激情福利司机影院| 成人亚洲欧美一区二区av| 97在线视频观看| 极品教师在线视频| 国产一级毛片七仙女欲春2| 国产精品精品国产色婷婷| av线在线观看网站| 男女啪啪激烈高潮av片| 欧美xxxx黑人xx丫x性爽| 日韩一区二区视频免费看| 五月伊人婷婷丁香| 男人和女人高潮做爰伦理| 亚洲欧美成人精品一区二区| 国产精品一区二区三区四区久久| 中文字幕制服av| 欧美xxⅹ黑人| 亚洲aⅴ乱码一区二区在线播放| 老女人水多毛片| 久久精品人妻少妇| 最后的刺客免费高清国语| 免费看不卡的av| 久久99热这里只有精品18| 久久久久精品性色| 国产老妇女一区| 久久久精品94久久精品| 极品少妇高潮喷水抽搐| 婷婷色综合大香蕉| 最近的中文字幕免费完整| 少妇高潮的动态图| 噜噜噜噜噜久久久久久91| 丝瓜视频免费看黄片| 日韩av不卡免费在线播放| 毛片一级片免费看久久久久| 国产亚洲精品av在线| 久久久成人免费电影| 国产成人aa在线观看| 国产精品一区二区三区四区免费观看| 亚洲av成人av| 亚洲精品一区蜜桃| 欧美日韩精品成人综合77777| 2021少妇久久久久久久久久久| 国产精品不卡视频一区二区| 日韩av免费高清视频| 欧美 日韩 精品 国产| 亚洲欧美精品自产自拍| 亚洲高清免费不卡视频| 精品少妇黑人巨大在线播放| 亚洲自偷自拍三级| 亚洲不卡免费看| 午夜日本视频在线| 午夜精品在线福利| 欧美日韩综合久久久久久| 一级二级三级毛片免费看| 国产 一区 欧美 日韩| 国产亚洲5aaaaa淫片| 插阴视频在线观看视频| 大香蕉久久网| 中国国产av一级| 国产一区有黄有色的免费视频 | 亚洲欧美中文字幕日韩二区| 一级毛片黄色毛片免费观看视频| 91久久精品国产一区二区成人| 一级毛片aaaaaa免费看小| 亚洲欧美成人综合另类久久久| 水蜜桃什么品种好| 大话2 男鬼变身卡| 午夜老司机福利剧场| 国产伦理片在线播放av一区| 国产亚洲午夜精品一区二区久久 | h日本视频在线播放| 国产亚洲午夜精品一区二区久久 | 国产黄色小视频在线观看| 亚洲人成网站高清观看| 亚洲综合色惰| 熟妇人妻不卡中文字幕| av免费观看日本| 中文乱码字字幕精品一区二区三区 | 国产伦精品一区二区三区视频9| 少妇熟女欧美另类| 亚洲电影在线观看av| 成人午夜高清在线视频| 国产精品人妻久久久久久| 在线观看一区二区三区| 一级黄片播放器| 日本免费在线观看一区| 在线观看免费高清a一片| 亚洲人成网站在线观看播放| 777米奇影视久久| 欧美成人午夜免费资源| 国产一区二区在线观看日韩| 国产单亲对白刺激| 国产三级在线视频| 久久久久性生活片| 亚洲精品视频女| 丰满少妇做爰视频| 看十八女毛片水多多多| 欧美激情在线99| 18+在线观看网站| xxx大片免费视频| 国产久久久一区二区三区| 蜜臀久久99精品久久宅男| 亚洲四区av| 人人妻人人澡人人爽人人夜夜 | 国产欧美日韩精品一区二区| 日韩欧美三级三区| 国产综合懂色| 国产成人免费观看mmmm| 国产乱人偷精品视频| 日本色播在线视频| 成人二区视频| 亚洲欧美精品专区久久| 三级国产精品片| 国产免费一级a男人的天堂| 男女视频在线观看网站免费| 国产精品一区二区三区四区免费观看| 一级毛片久久久久久久久女| a级一级毛片免费在线观看| 久久久久久久久久黄片| 真实男女啪啪啪动态图| 久久精品夜夜夜夜夜久久蜜豆| 99热这里只有是精品在线观看| 边亲边吃奶的免费视频| 99久国产av精品国产电影| 亚洲色图av天堂| 免费大片18禁| 亚洲无线观看免费| 偷拍熟女少妇极品色| 精品久久久久久久人妻蜜臀av| 国产一区亚洲一区在线观看| 一级爰片在线观看| 国产成年人精品一区二区| 国产在线一区二区三区精| 美女主播在线视频| 中文字幕制服av| 中文欧美无线码| 国产中年淑女户外野战色| 两个人的视频大全免费| 男人爽女人下面视频在线观看| 久久久久久久久久人人人人人人| 久久韩国三级中文字幕| 大香蕉久久网| 午夜福利视频精品| 欧美不卡视频在线免费观看| 成人午夜高清在线视频| 精品酒店卫生间| 观看美女的网站| 亚洲人成网站在线观看播放| 91久久精品国产一区二区三区| 成人综合一区亚洲| 2018国产大陆天天弄谢| 国产色婷婷99| 最近最新中文字幕免费大全7| 丝瓜视频免费看黄片| 欧美bdsm另类| 九九爱精品视频在线观看| 日韩欧美精品v在线| a级毛色黄片| 男女啪啪激烈高潮av片| 最近的中文字幕免费完整| 纵有疾风起免费观看全集完整版 | 国内精品一区二区在线观看| 成人午夜高清在线视频| 看十八女毛片水多多多| 精品久久久久久久末码| 精品国产露脸久久av麻豆 | 久久精品国产亚洲av涩爱| 美女xxoo啪啪120秒动态图| 啦啦啦韩国在线观看视频| 亚洲国产欧美在线一区| 午夜精品一区二区三区免费看| 七月丁香在线播放| 日本一本二区三区精品| 又粗又硬又长又爽又黄的视频| 国产 一区 欧美 日韩| 草草在线视频免费看| 久久精品久久久久久久性| 亚洲第一区二区三区不卡| 亚洲内射少妇av| 狂野欧美激情性xxxx在线观看| 日韩成人av中文字幕在线观看| 观看免费一级毛片| 99久久精品热视频| 嫩草影院入口| 国产伦理片在线播放av一区| 美女大奶头视频| 成人性生交大片免费视频hd| 麻豆成人av视频| 成年女人看的毛片在线观看| 国产综合精华液| 日本一本二区三区精品| www.av在线官网国产| 午夜福利在线在线| 久久久久网色| 国产成人freesex在线| 97超视频在线观看视频| 色尼玛亚洲综合影院| av在线老鸭窝| 欧美日韩在线观看h| 国产成人aa在线观看| 国产精品一区二区性色av| 综合色丁香网| 精品人妻一区二区三区麻豆| av网站免费在线观看视频 | 国产v大片淫在线免费观看| 欧美精品一区二区大全| 欧美成人一区二区免费高清观看| 一个人看视频在线观看www免费| 日本猛色少妇xxxxx猛交久久| 乱码一卡2卡4卡精品| 91av网一区二区| 好男人在线观看高清免费视频| 亚洲欧美日韩东京热| 日韩欧美三级三区| 亚洲精品中文字幕在线视频 | av网站免费在线观看视频 | 99热这里只有是精品50| 少妇被粗大猛烈的视频| 水蜜桃什么品种好| 国产有黄有色有爽视频| 搡老妇女老女人老熟妇| 天堂av国产一区二区熟女人妻| 日日撸夜夜添| 啦啦啦啦在线视频资源| 国产爱豆传媒在线观看| 亚洲高清免费不卡视频| 国产探花极品一区二区| 在线免费观看的www视频| 中文字幕制服av| 蜜臀久久99精品久久宅男| 麻豆成人午夜福利视频| 伦理电影大哥的女人| 久久99热这里只频精品6学生| 久久鲁丝午夜福利片| av卡一久久| 国产高清不卡午夜福利| 久久国内精品自在自线图片| 久久久亚洲精品成人影院| 亚洲av电影在线观看一区二区三区 | 国产欧美日韩精品一区二区| 听说在线观看完整版免费高清| 国产精品99久久久久久久久| 嫩草影院入口| 亚洲欧美一区二区三区黑人 | av国产免费在线观看| 午夜福利成人在线免费观看| 美女国产视频在线观看| 欧美97在线视频| 国产精品三级大全| 日韩在线高清观看一区二区三区| 国精品久久久久久国模美| 三级经典国产精品| 99久久精品国产国产毛片| 人妻夜夜爽99麻豆av| 午夜激情欧美在线| 国产 一区精品| 午夜福利成人在线免费观看| 美女黄网站色视频| 国产欧美日韩精品一区二区| 一级二级三级毛片免费看| 国内精品一区二区在线观看| 亚洲综合精品二区| av网站免费在线观看视频 | 两个人的视频大全免费| 水蜜桃什么品种好| 久久99精品国语久久久| 亚洲av电影不卡..在线观看| 国产熟女欧美一区二区| 美女cb高潮喷水在线观看| 我的老师免费观看完整版| 亚洲成人久久爱视频| 亚洲欧美精品专区久久| 久久久久久久久久久丰满| av又黄又爽大尺度在线免费看| 美女被艹到高潮喷水动态| av一本久久久久| 国产成人freesex在线| av国产久精品久网站免费入址| 国内精品美女久久久久久| 亚洲精品第二区| 免费高清在线观看视频在线观看| 精品99又大又爽又粗少妇毛片| 日日撸夜夜添| 亚洲av中文字字幕乱码综合| 69人妻影院| 精品不卡国产一区二区三区| 国产美女午夜福利| 国产一区二区在线观看日韩| 国产精品麻豆人妻色哟哟久久 | 国产探花极品一区二区| 一区二区三区乱码不卡18| av网站免费在线观看视频 | 22中文网久久字幕| 韩国高清视频一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 真实男女啪啪啪动态图| 国产高清国产精品国产三级 | 日韩一区二区视频免费看| 免费看美女性在线毛片视频| 最近中文字幕2019免费版| 精品一区二区三区视频在线| 国产真实伦视频高清在线观看| 亚洲国产精品国产精品| 国国产精品蜜臀av免费| 亚洲最大成人中文| 最近手机中文字幕大全| 国产精品国产三级专区第一集| 精品久久久久久成人av| 又粗又硬又长又爽又黄的视频| 色5月婷婷丁香| 美女内射精品一级片tv| 国产精品综合久久久久久久免费| 晚上一个人看的免费电影| 美女黄网站色视频| 欧美日本视频| 丰满少妇做爰视频| 亚洲av中文av极速乱| 18禁动态无遮挡网站| 97超视频在线观看视频| 国产伦理片在线播放av一区| 亚洲国产最新在线播放| 亚洲精品日韩在线中文字幕| 日本三级黄在线观看| 少妇的逼好多水| 久久精品国产亚洲av涩爱| 亚洲综合精品二区| 2021天堂中文幕一二区在线观| 国产在线一区二区三区精| 日本猛色少妇xxxxx猛交久久| 十八禁网站网址无遮挡 | 国产单亲对白刺激| 久久这里有精品视频免费| 2018国产大陆天天弄谢| 国产日韩欧美在线精品| 在线免费观看的www视频| 一级毛片aaaaaa免费看小| 久久久精品欧美日韩精品| 亚洲在久久综合| 69av精品久久久久久| 精品不卡国产一区二区三区| 国内精品宾馆在线| 久久热精品热| 国产精品一二三区在线看| 3wmmmm亚洲av在线观看| 精品亚洲乱码少妇综合久久| 99热网站在线观看| 男女那种视频在线观看| 国产精品一区二区三区四区久久| 国产在线男女| 国产视频首页在线观看| 搞女人的毛片| 欧美成人a在线观看| 国产av码专区亚洲av| 成年人午夜在线观看视频 | 久久久成人免费电影| 国产成人a区在线观看| 97在线视频观看| 亚洲在久久综合| 国产精品久久视频播放| 寂寞人妻少妇视频99o| 国产中年淑女户外野战色| 国产 一区精品| 精品久久久久久久久久久久久| 久久久久久久久久黄片| 小蜜桃在线观看免费完整版高清| 成人漫画全彩无遮挡| 少妇高潮的动态图| 亚洲欧美清纯卡通| 日韩欧美 国产精品| 麻豆av噜噜一区二区三区| 男女下面进入的视频免费午夜| 日本三级黄在线观看| 五月玫瑰六月丁香| 可以在线观看毛片的网站| 成人漫画全彩无遮挡| 看十八女毛片水多多多| videossex国产| 男的添女的下面高潮视频| 国产69精品久久久久777片| 久久久久精品久久久久真实原创| 天堂影院成人在线观看| 国产av不卡久久| 久久亚洲国产成人精品v| 国产精品久久久久久av不卡| 欧美xxⅹ黑人| 麻豆乱淫一区二区| 少妇人妻一区二区三区视频| 精品亚洲乱码少妇综合久久| 国产黄a三级三级三级人| 久久久久免费精品人妻一区二区| 日韩av在线大香蕉| 日本免费在线观看一区| 在线 av 中文字幕| 免费高清在线观看视频在线观看| 狠狠精品人妻久久久久久综合| 日韩大片免费观看网站| 亚洲在线自拍视频| 国产黄a三级三级三级人| 国产伦理片在线播放av一区| 91精品伊人久久大香线蕉| 国产免费一级a男人的天堂| 天堂俺去俺来也www色官网 | 国产精品久久视频播放| av在线蜜桃| 免费在线观看成人毛片| 哪个播放器可以免费观看大片| 天堂影院成人在线观看| 日韩av免费高清视频| 内射极品少妇av片p| 亚洲高清免费不卡视频| 国产免费一级a男人的天堂| 国产亚洲一区二区精品| 少妇被粗大猛烈的视频| 又黄又爽又刺激的免费视频.| 国产欧美另类精品又又久久亚洲欧美| 久久精品熟女亚洲av麻豆精品 | 久久草成人影院| 黑人高潮一二区| 中文字幕av在线有码专区| 一夜夜www| 人体艺术视频欧美日本| 亚洲成人久久爱视频| av播播在线观看一区| 久久草成人影院| 丰满人妻一区二区三区视频av| av专区在线播放| 一夜夜www| 九九在线视频观看精品| 日韩亚洲欧美综合| 2022亚洲国产成人精品| 99久国产av精品| 亚洲国产精品专区欧美| 欧美潮喷喷水| 国产真实伦视频高清在线观看| 久久久久久久国产电影| 亚州av有码|