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

    基于高光譜分析的長蕩湖水質(zhì)監(jiān)測技術(shù)研究

    2022-08-27 07:34:12蔣志昊梁文廣石一凡王冬梅唐大偉宋亞君
    江蘇水利 2022年8期
    關(guān)鍵詞:長蕩湖懸浮物濁度

    蔣志昊,梁文廣,石一凡,王冬梅,唐大偉,宋亞君

    (1.江蘇省水利科學(xué)研究院,江蘇南京 210000;2.江蘇省水文水資源勘測局常州分局,江蘇常州 213300)

    1 研究背景

    目前,利用衛(wèi)星遙感反演水質(zhì)參數(shù)和人工實測這2 種常規(guī)監(jiān)測方法為主要的湖泊水質(zhì)監(jiān)測手段。常規(guī)的水質(zhì)監(jiān)測技術(shù)耗時耗力,投入產(chǎn)出比低,難以覆蓋區(qū)域范圍內(nèi)所有地點水質(zhì)數(shù)據(jù),不能滿足多樣化的監(jiān)測需求。當前,多光譜遙感技術(shù)具有時態(tài)多、范圍廣的動態(tài)監(jiān)測優(yōu)勢,廣泛應(yīng)用于水體參數(shù)的反演研究中,但光譜波段少、光譜分辨率較低[1],影響水質(zhì)參數(shù)反演精度。高光譜技術(shù)具有光譜波段多和連續(xù)等優(yōu)勢,海量的光譜信息可反映真實水表光譜信息,是常規(guī)水質(zhì)監(jiān)測和遙感水質(zhì)監(jiān)測的重要補充,為區(qū)域性水環(huán)境監(jiān)測提供新的技術(shù)手段。曹引等[2]基于南四湖實測的高光譜數(shù)據(jù)和同步水質(zhì)采樣分析數(shù)據(jù)構(gòu)建了反演南四湖水體葉綠素a濃度的多種數(shù)學(xué)模型;鞏彩蘭等[3]利用地物光譜測量及同步配套水質(zhì)采樣分析實驗,對黃浦江全河段水體進行調(diào)查研究,并以TN 和懸浮物為例得出最佳反射率反演組合;包起帆等[4]利用星地同步監(jiān)測實驗,建立長江口含沙量濃度高光譜反演模型;吳廷寬等[5]以貴陽市百花湖為例,對水體進行光譜實測,探討湖泊水質(zhì)高光譜反射率與水質(zhì)參數(shù)濃度之間的定量關(guān)系,對湖泊富營養(yǎng)化進行監(jiān)測評價。由上述研究成果可知,不同湖泊之間的下墊面條件有所差別,反演模型普適性較差,對典型湖泊光譜特性和水質(zhì)參數(shù)之間定量關(guān)系需單獨分析。

    根據(jù)常州水環(huán)境監(jiān)測中心2008—2018 年水質(zhì)監(jiān)測成果,長蕩湖整體水質(zhì)基本處于劣Ⅴ類??紤]到長蕩湖在西太湖流域中的重要過水作用和對周邊水系的影響,提高長蕩湖水質(zhì)監(jiān)測精度水平,實現(xiàn)全湖泊、多時段常態(tài)化管理,對改善水環(huán)境、實行長效化管理意義重大。因此,本研究以常州市長蕩湖為研究對象,依照《地表水和污水監(jiān)測技術(shù)規(guī)范》進行水質(zhì)采樣,獲取采樣點的水質(zhì)參數(shù),利用高光譜掃描儀提取對應(yīng)點位水表光譜信息。通過水質(zhì)參數(shù)與高光譜的相關(guān)性分析,分析比較不同水質(zhì)指標的光譜,篩選出對光譜反射率敏感的水質(zhì)參數(shù),建立可見光至近紅外波長范圍內(nèi)多種水質(zhì)參數(shù)-高光譜反演模型,并驗證分析模型精度和適用性,確定適用于長蕩湖的水質(zhì)高光譜監(jiān)測模型。

    2 材料與方法

    2.1 研究區(qū)概況

    長蕩湖(31°33′~31°40′N,119°30′~119°37′E)位于江蘇省常州市境內(nèi),地處亞熱帶季風氣候帶,湖泊南北長13.6 km,東西寬9.3 km,現(xiàn)有水域面積約84 km2,平均水深1.1~1.2 m?,F(xiàn)狀湖盆地形平坦,無顯著起伏,北半部湖盆水稍深,南半部水淺,多沼澤性蘆葦淺灘,淤積嚴重。

    2.2 水質(zhì)采樣及處理

    為保證實驗效果,要求采樣當日天氣晴朗,湖泊水體流速不宜過高,以提供較為穩(wěn)定的實驗環(huán)境。研究區(qū)域選定后,在湖泊的進出水區(qū)、中心敞水區(qū)、濱岸帶等各類型典型區(qū)域設(shè)置一定數(shù)量的采樣點,于2020 年5 月4 日利用YSI 公司水質(zhì)多參數(shù)分析儀和實驗室水質(zhì)采樣化驗分析的方式,獲取水溫、DO、ORP、電導(dǎo)率、透明度、NH3-N、NO3-N、TP、TN、CODcr、CODMn、懸浮物SS、濁度、Chl-a 和pH 等15 項水質(zhì)參數(shù)。本次采樣共設(shè)20 個點位,取水表0.5 m表層水,點位具體分布情況見表1。

    表1 長蕩湖布設(shè)點分布

    通過傳感器現(xiàn)場測量以及實驗室內(nèi)水樣測定,長蕩湖5月水質(zhì)統(tǒng)計結(jié)果如表2所示。結(jié)果顯示調(diào)研區(qū)域內(nèi)各水質(zhì)參數(shù)具有顯著空間差異。根據(jù)國家《地表水環(huán)境質(zhì)量標準》(GB3838—2001)和現(xiàn)場采樣實驗室分析結(jié)果,當日長蕩湖CODMn、COD 全湖平均達到Ⅱ類水體水質(zhì)標準,南部水質(zhì)劣于北部,大部分達到Ⅲ類標準。TP、TN 質(zhì)量濃度均低于Ⅴ類地表水標準,湖泊水質(zhì)為劣Ⅴ類,溧陽市域內(nèi)圍網(wǎng)養(yǎng)殖嚴重,湖泊富營養(yǎng)化程度嚴重。

    表2 長蕩湖水質(zhì)參數(shù)統(tǒng)計結(jié)果

    2.3 現(xiàn)場水體高光譜數(shù)據(jù)

    高光譜采集與水質(zhì)參數(shù)測量同步進行,根據(jù)布設(shè)點位的序號對各點依次進行測量。本次實驗使用的是美國ASD 公司生產(chǎn)的FieldSpec?4 Std-Res野外便攜式高光譜儀,波長范圍350~2 500nm,能夠捕獲可見和近紅外光譜。在350~1 000 nm范圍內(nèi),光譜采樣間隔為1.4 nm,光譜分辨率為3 nm,而在1 000~2 500 nm 范圍內(nèi),光譜采樣間隔為2 nm,光譜分辨率為10nm。獲取的原始數(shù)據(jù)經(jīng)過數(shù)據(jù)篩選,反射率計算和異常值剔除處理后,與水質(zhì)參數(shù)進行擬合分析。

    2.4 光譜參數(shù)及相關(guān)性分析

    將經(jīng)實驗室檢測15 種水質(zhì)參數(shù)和經(jīng)預(yù)處理的光譜反射率進行相關(guān)性分析。分別采樣單波段[6-7](Single Band,SB)、波段差值[8](Difference of Band,DB)和波段比值(Ratio of Band,RB)[9-11](式1~2),比較350~1 000 nm 波長范圍內(nèi)光譜反射率和各參數(shù)間的相關(guān)關(guān)系,篩選出對光譜敏感的水質(zhì)參數(shù)。式中:BD和BR分別為波段差值和波段比值;Ri和Rj為對應(yīng)i和j波段的光譜反射率。

    2.5 模型構(gòu)建

    利用檢測到的建模樣本,使用各種算法建立光譜數(shù)據(jù)與通過傳統(tǒng)方法檢測的屬性值之間的線性或非線性回歸模型。本次研究根據(jù)光譜相關(guān)性分析結(jié)果,參照水質(zhì)參數(shù)與光譜反射率之間關(guān)系所得規(guī)律選擇最優(yōu)波段組合,分別建立指數(shù)、線性、多項式、冪函數(shù)和對數(shù)模型。

    2.6 模型驗證

    從所有采樣樣本中隨機選取2/3的樣本用于構(gòu)建數(shù)學(xué)模型,其他1/3 的樣本用來檢驗?zāi)P蛿M合精度。選擇相關(guān)性最高的數(shù)學(xué)模型進行水質(zhì)反演,并采用F檢驗法對模型顯著性進行驗證,模型的精度評價指標為絕對誤差(Absloute Error,AE)、均方根誤差(Root Mean Squared Error,RMSE)、相對誤差(Relative Error,RE)和平均相對誤差(Mean Relative Error,MRE)來表示進行評價。通過R2、RSE、Std.Res(殘差標準差)將驗證組樣本真實值與模型模擬值進行對比,模型的評價標準為R2越高,且RMSE、RE以及MRE越低,則表示模型的擬合效果越好。實測值與模擬值散點越接近1∶1 線,代表反演的精度越高,最終確定反演效果最佳的模型。

    3 結(jié)果與分析

    3.1 水體光譜特征

    本次實測光譜曲線在810nm 開始下降,而由于水分子在近紅外波段對光有強吸收作用,波長大于1 000 nm 以后反射率迅速下降并接近于0,所以選擇350~1 000 nm 波長范圍進行樣本分析。本次實驗最終選取的20 個采樣點光譜處理結(jié)果如圖1 所示。根據(jù)光譜實測結(jié)果,長蕩湖水體在可見光范圍內(nèi)共有3個反射峰值區(qū),分別集中于550 nm、700 nm和800 nm 波長附近,分別對應(yīng)可見光綠光、紅光波長。其中,在550 nm 處的反射率對應(yīng)最高峰值,在700 nm附近出現(xiàn)2個強吸收帶,反射率陡降,800 nm波長以外的反射率逐漸下降。

    圖1 長蕩湖采樣點光譜反射率響應(yīng)曲線

    2.2 水質(zhì)參數(shù)與高光譜波段相關(guān)性分析

    將2020 年5 月4 日同步獲取的上述15 項水質(zhì)參數(shù)和實測光譜數(shù)據(jù)進行Pearson相關(guān)分析,獲取水質(zhì)參數(shù)相關(guān)系數(shù)最佳的波段或波段組合反射率值(表3)。單波段相關(guān)性分析優(yōu)點在于模型構(gòu)建簡單且直觀,能排除波段自相關(guān)影響,并且具有豐富的先驗知識和理論基礎(chǔ)作為支撐;雙波段相關(guān)性分析可以一定程度上消除測量時代入人的噪聲、測量環(huán)境變化帶來的誤差和大氣的影響,達到信息增強與壓縮的目的,以利于更好建立水表反射率與水質(zhì)參數(shù)化驗指標之間的關(guān)系。本文分別采用波段差值和比值的方法進行相關(guān)性分析。

    表3 長蕩湖水質(zhì)參數(shù)與350~1 000 nm范圍內(nèi)光譜反射率相關(guān)性分析

    根據(jù)最優(yōu)波段相關(guān)性表,雙波段組合與光譜反射率之間的相關(guān)關(guān)系整體優(yōu)于單波段,其中波段比值與水溫之間具有最高的相關(guān)性,R=0.96;波段差值對懸浮物SS、電導(dǎo)率、水溫、濁度相關(guān)性最強,R=0.94;單波段對懸浮物SS 具有最高的相關(guān)性,在波長為816 nm處取得,R=0.84。

    將各種波段組合下最優(yōu)相關(guān)性大于0.7的水質(zhì)參數(shù)看作光譜反射率的敏感參數(shù),進行建模分析和對比驗證,最終選取TP、懸浮物SS、濁度、水溫和透明度等5 項參數(shù)。圖2~6 為上述5 種參數(shù)分別與單波段、波段差值、波段比值任意組合的相關(guān)關(guān)系分布圖。5 種參數(shù)均在近紅外波段的波長處與反射率具有較好的相關(guān)性,TP、懸浮物SS、濁度、水溫和透明度在單波段的最大相關(guān)系數(shù)分別為0.815、0.84、0.78、0.74、-0.74。圖2~6 為上述5 項參數(shù)與350~1 000 nm 波長范圍內(nèi),3 種波段反射率組合的相關(guān)性分布圖。單波段光譜反射率相關(guān)性最高的波段集中位于近紅外波段附近,900 nm之后伴隨強水汽吸收帶,相關(guān)性迅速減弱。波段差值組合的相關(guān)性分布圖具有一致的特征規(guī)律,雙波段在均位于500 nm 波長范圍內(nèi)的組合與各參數(shù)相關(guān)性較弱,對應(yīng)反射率曲線為Chl-a 和黃色物質(zhì)的強烈吸收帶,此區(qū)間各波段反射率較低。500 nm以上時相關(guān)性較為顯著,且波長差在200 nm 以內(nèi),波段反射率差值與參數(shù)相關(guān)性較強,這與顆粒懸浮物的散射作用有關(guān)。波段相除組合具有相似的特征,但是當波長超過930 nm后相關(guān)性突然下降。

    圖2 TP與光譜反射率相關(guān)性分析

    圖3 懸浮物SS與光譜反射率相關(guān)性分析

    圖4 濁度與光譜反射率相關(guān)性分析

    圖5 水溫與光譜反射率相關(guān)性分析

    圖6 透明度與光譜反射率相關(guān)性分析

    2.3 模型構(gòu)建

    經(jīng)驗?zāi)P停?2-15]是通過對觀測對象表層光譜特性和水質(zhì)參數(shù)同步實測,并進行適當?shù)慕y(tǒng)計分析得到水質(zhì)參數(shù)的估測算法。本次實驗以單波段、波段差值、波段比值為自變量,分別利用指數(shù)函數(shù)、線性函數(shù)、多項式函數(shù)、冪函數(shù)和對數(shù)函數(shù)對長蕩湖TP、懸浮物SS、濁度、水溫和透明度這5 項參數(shù)進行模擬(表4~8)。

    表4 長蕩湖水體TP高光譜擬合模型

    利用單波段、波段差值和波段比值組合,利用多種數(shù)學(xué)模型,建立光譜反射率和TP質(zhì)量濃度之間的定量關(guān)系。結(jié)果顯示,比值模型整體相關(guān)性最優(yōu),各類函數(shù)擬合均保持R2>0.8,其中將864 nm 和914 nm處的反射率比值作為自變量,建立的一元二次多項式方程模擬的TP 質(zhì)量濃度具有最高的相關(guān)性,R2為0.833,并且通過F檢驗。

    利用單波段、波段差值和波段比值組合,結(jié)合多種數(shù)學(xué)模型,建立光譜反射率和懸浮物SS之間的定量關(guān)系。結(jié)果顯示,用波段差值和波段比值進行擬合均能取得較好的效果,其中將928 nm和964 nm處的反射率比值作為自變量,建立的一元二次多項式方程模擬的懸浮物質(zhì)量濃度具有最高的相關(guān)性,R2為0.926,并且通過F檢驗。

    利用單波段、波段差值和波段比值組合,利用多種數(shù)學(xué)模型,建立光譜反射率和濁度之間的定量關(guān)系。結(jié)果顯示,單波段反射率與濁度在可見光范圍內(nèi)雖然有較好的相關(guān)性,但不適宜用本文提供的經(jīng)驗?zāi)P头ㄟM行擬合。雙波段進行擬合的效果明顯優(yōu)于單波段,并且以909 nm 和923 nm 處的反射率差值作為自變量,建立的線性方程模擬的水體濁度具有最高的相關(guān)性,R2為0.922。

    利用單波段、波段差值和波段比值組合,利用多種數(shù)學(xué)模型,建立光譜反射率和水溫之間的定量關(guān)系。結(jié)果顯示,雙波段模擬效果明顯優(yōu)于單波段,將779 nm 和740 nm 處的反射率比值作為自變量,建立的線性方程和多相式方程模擬的pH 值相關(guān)性最高,R2均為0.941,并且通過F檢驗??紤]到線性方程式模型更加簡易,最終選擇線性模型作為反演模型。

    表5 長蕩湖水體懸浮物SS高光譜擬合模型

    表6 長蕩湖水體濁度高光譜擬合模型

    利用單波段、波段差值和波段比值組合,利用多種數(shù)學(xué)模型,建立光譜反射率和透明度之間的定量關(guān)系。結(jié)果顯示,利用差值組合擬合的效果最佳,其中將971 nm 和930 nm 處的反射率差值作為自變量,建立的指數(shù)方程模擬的透明度具有最高的相關(guān)性,R2為0.959,并且通過F檢驗。

    對上述5種參數(shù)最佳波段組合下的最優(yōu)模型進行整體評價,除TP 的擬合采用的是波段比值,懸浮物、濁度、水溫和透明度均宜采用波段差值的組合。對透明度模擬的R值最高,對水溫進行模擬后MRE取得最小值(表9)。

    表9 高光譜擬合模型整體評價

    2.4 模型評價

    根據(jù)最優(yōu)半經(jīng)驗?zāi)P皖A(yù)測值與TP、懸浮物SS、濁度、水溫和透明度實測值的對比,根據(jù)本文模型預(yù)測的長蕩湖TP、濁度、透明度均與實測值高度吻合。TPR2=0.97,RES=0.005,std.Res=0.015;透明度R2=0.97,RES=0.069,std.Res=1.204;透明度R2=0.9,RES=-0.971,std.Res=3.69。預(yù)測值與實測值均勻分布在X=Y線的兩側(cè)且緊鄰分布。水溫和懸浮物SS的模擬相關(guān)性弱于上述3 項,但R2均大于0.74。整體而言,本文建立模型對長蕩湖TP、懸浮物、濁度、水溫和透明度的預(yù)測均取得較好的效果。圖7為擬合效果對比圖。

    圖7 長蕩湖參數(shù)反演模型精度評價

    3 結(jié)論與討論

    本文根據(jù)實測長蕩湖的光譜特征,利用同步獲取20個點位的15項常規(guī)水質(zhì)參數(shù)。通過光譜反射率和水質(zhì)參數(shù)相關(guān)性分析,篩選出TP、懸浮物SS、濁度、水溫和透明度等5項相關(guān)性強的水質(zhì)參數(shù),建立適用于長蕩湖的多種水質(zhì)-高光譜反演模型,并進行模型驗證和適用性分析等,得出以下結(jié)論:

    (1)根據(jù)2021年5月對長蕩湖15項水質(zhì)特征參數(shù)的采樣結(jié)果分析,當日長蕩湖CODMn、COD全湖平均達到Ⅱ類水體水質(zhì)標準,南部水質(zhì)劣于北部,大都達到Ⅲ類標準。TP、TN 質(zhì)量濃度均低于Ⅴ類地表水標準,湖泊水質(zhì)為劣Ⅴ類。受溧陽市域內(nèi)大量圍網(wǎng)養(yǎng)殖的影響,全湖各位置懸浮物SS、濁度、Chl-a 和透明度差異較大,南部湖泊富營養(yǎng)化程度嚴重。

    表7 長蕩湖水體水溫高光譜擬合模型

    表8 長蕩湖水體透明度高光譜擬合模型

    (2)根據(jù)光譜實測結(jié)果,長蕩湖水體在可見光范圍內(nèi)共有3 個反射峰值區(qū),分別集中于550 nm、700 nm 和800 nm 波長附近,分別對應(yīng)可見光綠光、紅光波長。其中,在550 nm處的反射率對應(yīng)最高峰值,在700 nm附近出現(xiàn)2個強吸收帶,反射率陡降,800 nm波長以外的反射率逐漸下降。

    (3)在實測水表光譜的基礎(chǔ)上,分別利用單波段、雙波段差值、雙波段比值3 種組合方式與15 項水質(zhì)參數(shù)進行相關(guān)性分析。結(jié)果顯示,雙波段普遍比單波段具有更強相關(guān)性。單波段則受限于單邊自變量,相關(guān)性在給定波長內(nèi)表現(xiàn)出正負性。TP、懸浮物SS、濁度、水溫和透明度在不同波段反射率組合條件下Rmax均大于0.7,表現(xiàn)出較強相關(guān)性,對光輻射敏感性較高。受自身作用機制和被其他水質(zhì)參數(shù)掩蓋特征的影響,另外10種水質(zhì)參數(shù)輻射特性表現(xiàn)不太顯著。

    (4)基于半經(jīng)驗?zāi)P?,分別利用指數(shù)、線性、多項式、冪函數(shù)和對數(shù)函數(shù),篩選出最優(yōu)波段組合下的最佳反演模型。結(jié)果顯示,多項式模型對TP、懸浮物SS和透明度的反演效果最佳,相關(guān)系數(shù)平方分別達到R2864/914=0.92,R2928-964=0.93,R2971-930=0.96。濁度和水溫更適合用線性函數(shù)反演,兩者反演的相關(guān)系數(shù)平方分別為R2909-923=0.92,R2779-740=0.94。將上述模型帶入檢驗集中進行精度驗證,均表現(xiàn)出于真實值相近的結(jié)果。TP 預(yù)測值和真實值之間R2=0.97,RES=0.005;懸浮物SSR2=0.74,RES=1.014;濁度R2=0.9,RES=0.971;水溫R2=0.80,RES=0.144;透明度R2=0.97,RES=0.069。

    (5)高光譜技術(shù)在湖泊水質(zhì)參數(shù)研究中發(fā)揮重要作用,利用實測光譜和水質(zhì)參數(shù)建立的數(shù)學(xué)模型有利于提高遙感反演精度,為實際應(yīng)用和推廣建立理論依據(jù)。同時,長蕩湖內(nèi)構(gòu)筑物分布復(fù)雜,具有飲用水源地、水上景觀、餐飲船聚集點和圍網(wǎng)養(yǎng)殖等多種功能區(qū)?,F(xiàn)已在全湖開展退圩還湖工作,高精度高光譜反演模型結(jié)合衛(wèi)星遙感可有效監(jiān)測工程實施過程中全湖水質(zhì)改善狀況,并對各功能區(qū)的分配決策給出指導(dǎo)性意見。

    猜你喜歡
    長蕩湖懸浮物濁度
    豐水期長蕩湖水質(zhì)評價與微生物群落結(jié)構(gòu)特征研究
    懸浮物對水質(zhì)氨氮測定的影響
    化工管理(2022年14期)2022-12-02 11:43:52
    丙烯酰胺強化混凝去除黑河原水濁度的研究
    動態(tài)濁度補償技術(shù)在總磷在線自動監(jiān)測儀上的應(yīng)用
    云南化工(2021年6期)2021-12-21 07:31:06
    11°角應(yīng)用于啤酒過濾濁度測量
    『長蕩湖』八鮮宴
    美食(2017年9期)2017-09-15 08:52:47
    廢水懸浮物國標檢測方法的修改建議
    壓裂返排液中懸浮物去除的室內(nèi)研究
    雷州灣海域懸浮物濃度遙感監(jiān)測研究
    長蕩湖河蟹產(chǎn)業(yè)現(xiàn)狀及發(fā)展對策
    99在线视频只有这里精品首页| 亚洲精品国产精品久久久不卡| 在线观看美女被高潮喷水网站 | 一本久久中文字幕| 在线看三级毛片| 两个人看的免费小视频| 欧美区成人在线视频| 神马国产精品三级电影在线观看| 有码 亚洲区| 18禁裸乳无遮挡免费网站照片| 首页视频小说图片口味搜索| 欧美乱色亚洲激情| 亚洲精品一区av在线观看| 久久中文看片网| 丝袜美腿在线中文| 国产97色在线日韩免费| 国产亚洲精品综合一区在线观看| 亚洲精品亚洲一区二区| 久久久国产精品麻豆| 给我免费播放毛片高清在线观看| 国产久久久一区二区三区| 久久精品亚洲精品国产色婷小说| 精品久久久久久,| 床上黄色一级片| 欧美成人a在线观看| 女人被狂操c到高潮| 午夜久久久久精精品| 校园春色视频在线观看| 欧美成人免费av一区二区三区| 日韩高清综合在线| 色精品久久人妻99蜜桃| 国产高清有码在线观看视频| 色综合站精品国产| АⅤ资源中文在线天堂| 久久精品国产清高在天天线| 国产成年人精品一区二区| 一级a爱片免费观看的视频| 国产男靠女视频免费网站| 久久精品91蜜桃| 国产精品亚洲一级av第二区| 日韩精品中文字幕看吧| 国产免费男女视频| 在线天堂最新版资源| 精品久久久久久,| 国产精品精品国产色婷婷| 午夜精品久久久久久毛片777| 国产av一区在线观看免费| 久久久久免费精品人妻一区二区| 国产真实伦视频高清在线观看 | 国产精品久久久人人做人人爽| 亚洲欧美日韩卡通动漫| 久久伊人香网站| www日本在线高清视频| 久久久久亚洲av毛片大全| 亚洲精品粉嫩美女一区| 国产成人aa在线观看| 人妻丰满熟妇av一区二区三区| 国产精品亚洲一级av第二区| 哪里可以看免费的av片| 国产美女午夜福利| 色综合婷婷激情| 久久久久久九九精品二区国产| 欧美黑人欧美精品刺激| 熟女人妻精品中文字幕| 九色成人免费人妻av| 青草久久国产| 国产精品电影一区二区三区| 国产精品久久久久久亚洲av鲁大| 久久久国产精品麻豆| 一进一出抽搐动态| 好看av亚洲va欧美ⅴa在| 此物有八面人人有两片| 国产精品女同一区二区软件 | 黄片大片在线免费观看| 可以在线观看毛片的网站| 国产极品精品免费视频能看的| 国产精品一及| 一区二区三区免费毛片| 国产成人a区在线观看| 国产视频一区二区在线看| 黄色女人牲交| 亚洲成av人片免费观看| 免费一级毛片在线播放高清视频| 国产av在哪里看| 亚洲狠狠婷婷综合久久图片| 少妇的丰满在线观看| 亚洲最大成人中文| 国产主播在线观看一区二区| 久久久久国内视频| 最新美女视频免费是黄的| 亚洲人成网站在线播放欧美日韩| 国产成年人精品一区二区| 最近在线观看免费完整版| 久久久精品大字幕| 好男人电影高清在线观看| 露出奶头的视频| 三级毛片av免费| 成人精品一区二区免费| 日韩 欧美 亚洲 中文字幕| 欧美成人性av电影在线观看| 国内久久婷婷六月综合欲色啪| 精品福利观看| 88av欧美| 一级黄片播放器| 亚洲午夜理论影院| www日本黄色视频网| 日日干狠狠操夜夜爽| 精品日产1卡2卡| 亚洲最大成人中文| 久久中文看片网| 成年人黄色毛片网站| 午夜精品一区二区三区免费看| 精品一区二区三区人妻视频| 亚洲av电影在线进入| 老司机午夜十八禁免费视频| 久久久久国产精品人妻aⅴ院| 亚洲成人精品中文字幕电影| 母亲3免费完整高清在线观看| 国产精品久久电影中文字幕| 动漫黄色视频在线观看| 淫秽高清视频在线观看| 亚洲av电影不卡..在线观看| 午夜两性在线视频| 麻豆国产97在线/欧美| 国产精品美女特级片免费视频播放器| 熟妇人妻久久中文字幕3abv| 丰满的人妻完整版| 国产成人av教育| 91在线观看av| 五月伊人婷婷丁香| 桃红色精品国产亚洲av| netflix在线观看网站| 男女午夜视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 精品99又大又爽又粗少妇毛片 | 国产伦在线观看视频一区| 少妇人妻一区二区三区视频| 日韩高清综合在线| 9191精品国产免费久久| 丰满的人妻完整版| bbb黄色大片| 国产精品亚洲美女久久久| av国产免费在线观看| 亚洲黑人精品在线| 国内精品久久久久精免费| 日韩免费av在线播放| 嫩草影院精品99| 在线看三级毛片| 色老头精品视频在线观看| 午夜免费观看网址| 无人区码免费观看不卡| 老熟妇乱子伦视频在线观看| 国产精品一区二区三区四区久久| 最新在线观看一区二区三区| 午夜两性在线视频| 香蕉丝袜av| 色av中文字幕| 我的老师免费观看完整版| 欧美性猛交黑人性爽| 色视频www国产| 成人国产综合亚洲| АⅤ资源中文在线天堂| 精品99又大又爽又粗少妇毛片 | 欧美不卡视频在线免费观看| 波多野结衣高清作品| 亚洲人与动物交配视频| 国产精品1区2区在线观看.| 国产精品免费一区二区三区在线| aaaaa片日本免费| 国产探花极品一区二区| 亚洲av免费在线观看| 国产成人aa在线观看| 啦啦啦观看免费观看视频高清| 两性午夜刺激爽爽歪歪视频在线观看| 一本综合久久免费| 两性午夜刺激爽爽歪歪视频在线观看| 人妻夜夜爽99麻豆av| 国内少妇人妻偷人精品xxx网站| 亚洲国产精品999在线| 日韩欧美三级三区| 嫩草影院精品99| 性色avwww在线观看| 非洲黑人性xxxx精品又粗又长| 99久久成人亚洲精品观看| 亚洲av成人精品一区久久| 日韩欧美在线乱码| 国产老妇女一区| 午夜免费成人在线视频| 婷婷精品国产亚洲av在线| 搡女人真爽免费视频火全软件 | 日韩欧美精品v在线| 久久6这里有精品| 久久久久性生活片| 日本黄大片高清| 在线a可以看的网站| 日韩精品青青久久久久久| 久久精品夜夜夜夜夜久久蜜豆| 国产黄色小视频在线观看| 国产亚洲精品久久久com| 中文字幕精品亚洲无线码一区| 亚洲成a人片在线一区二区| 亚洲国产精品久久男人天堂| 亚洲在线自拍视频| 久久精品国产综合久久久| 日本黄大片高清| av欧美777| 51午夜福利影视在线观看| 好男人在线观看高清免费视频| 老司机午夜福利在线观看视频| 免费av观看视频| 国产爱豆传媒在线观看| 日韩欧美精品v在线| 亚洲av日韩精品久久久久久密| 午夜精品在线福利| 制服人妻中文乱码| 国产午夜福利久久久久久| 国产精品一区二区三区四区久久| 99riav亚洲国产免费| 窝窝影院91人妻| 亚洲成人中文字幕在线播放| 熟女人妻精品中文字幕| 麻豆成人av在线观看| 亚洲精品美女久久久久99蜜臀| av片东京热男人的天堂| 不卡一级毛片| 夜夜爽天天搞| 午夜视频国产福利| 黄片小视频在线播放| 在线看三级毛片| 久久香蕉精品热| 国产国拍精品亚洲av在线观看 | 老熟妇仑乱视频hdxx| 一个人免费在线观看电影| 精品久久久久久久久久免费视频| 国产一区二区在线观看日韩 | 婷婷精品国产亚洲av| 麻豆国产97在线/欧美| 精品99又大又爽又粗少妇毛片 | 三级男女做爰猛烈吃奶摸视频| 90打野战视频偷拍视频| 男女那种视频在线观看| 国产精品香港三级国产av潘金莲| 国产亚洲欧美98| 熟女电影av网| 久久亚洲真实| 国产精品亚洲美女久久久| 美女黄网站色视频| 亚洲国产中文字幕在线视频| 久久国产乱子伦精品免费另类| 亚洲精品久久国产高清桃花| 日韩精品中文字幕看吧| 俄罗斯特黄特色一大片| 真人做人爱边吃奶动态| 日韩欧美三级三区| 国产av一区在线观看免费| 亚洲美女黄片视频| 国产成人av教育| 久久香蕉精品热| 两个人的视频大全免费| 精品人妻一区二区三区麻豆 | 午夜福利成人在线免费观看| 九色国产91popny在线| 国产精品国产高清国产av| 久久久久久九九精品二区国产| 免费电影在线观看免费观看| 色综合站精品国产| 操出白浆在线播放| 国产伦精品一区二区三区四那| 亚洲一区二区三区色噜噜| 九九久久精品国产亚洲av麻豆| 久久精品91无色码中文字幕| 久久草成人影院| 亚洲国产日韩欧美精品在线观看 | 中文字幕精品亚洲无线码一区| 亚洲无线观看免费| 97人妻精品一区二区三区麻豆| 欧美成人a在线观看| 久久精品国产综合久久久| 深夜精品福利| 最新美女视频免费是黄的| www日本黄色视频网| 午夜激情福利司机影院| 一进一出抽搐gif免费好疼| 免费搜索国产男女视频| 老司机深夜福利视频在线观看| 久久精品影院6| 欧美中文综合在线视频| 人人妻人人澡欧美一区二区| 99精品久久久久人妻精品| 小说图片视频综合网站| 最近最新中文字幕大全免费视频| 久久精品国产综合久久久| 岛国视频午夜一区免费看| 色av中文字幕| av国产免费在线观看| 非洲黑人性xxxx精品又粗又长| 在线天堂最新版资源| 一个人免费在线观看电影| 中文在线观看免费www的网站| 免费在线观看影片大全网站| 免费高清视频大片| 亚洲欧美日韩高清专用| 国产在视频线在精品| www日本在线高清视频| 午夜老司机福利剧场| 国产在线精品亚洲第一网站| av天堂在线播放| 久久天躁狠狠躁夜夜2o2o| 中文字幕高清在线视频| 亚洲乱码一区二区免费版| 久久久久国内视频| 成人精品一区二区免费| 日本 av在线| 18禁黄网站禁片午夜丰满| 婷婷亚洲欧美| 免费高清视频大片| 久久精品国产99精品国产亚洲性色| 免费无遮挡裸体视频| 婷婷六月久久综合丁香| 亚洲国产色片| 韩国av一区二区三区四区| 一区福利在线观看| 欧美乱妇无乱码| 国产免费男女视频| 色综合亚洲欧美另类图片| 麻豆成人午夜福利视频| 无限看片的www在线观看| 久久这里只有精品中国| 日韩亚洲欧美综合| 国产伦在线观看视频一区| 99久久精品热视频| av女优亚洲男人天堂| 国产午夜精品论理片| 搡女人真爽免费视频火全软件 | 桃红色精品国产亚洲av| 在线观看舔阴道视频| 国产私拍福利视频在线观看| 免费在线观看亚洲国产| 麻豆成人av在线观看| 欧美性猛交黑人性爽| 国产精品久久视频播放| 最近视频中文字幕2019在线8| 在线免费观看不下载黄p国产 | 日韩欧美免费精品| 高清在线国产一区| 麻豆久久精品国产亚洲av| 俄罗斯特黄特色一大片| 蜜桃亚洲精品一区二区三区| 一个人看的www免费观看视频| 精品乱码久久久久久99久播| 亚洲人成伊人成综合网2020| 99国产精品一区二区三区| 色吧在线观看| x7x7x7水蜜桃| 99久久九九国产精品国产免费| 亚洲精品美女久久久久99蜜臀| 女人被狂操c到高潮| 国产精品久久久久久久电影 | a级一级毛片免费在线观看| 免费看a级黄色片| 国产成人影院久久av| 搡老熟女国产l中国老女人| 99久久精品热视频| 51国产日韩欧美| 亚洲va日本ⅴa欧美va伊人久久| 色吧在线观看| 99久久精品国产亚洲精品| 欧美最新免费一区二区三区 | 亚洲美女黄片视频| 欧美中文日本在线观看视频| 亚洲国产精品sss在线观看| 97超级碰碰碰精品色视频在线观看| 91字幕亚洲| 亚洲欧美日韩高清专用| 日韩人妻高清精品专区| 在线观看日韩欧美| 国产高清视频在线播放一区| 久99久视频精品免费| 欧美三级亚洲精品| 99在线视频只有这里精品首页| 国产一区二区三区视频了| 日韩亚洲欧美综合| 好男人电影高清在线观看| 日韩有码中文字幕| 嫩草影视91久久| 精品无人区乱码1区二区| 精品人妻偷拍中文字幕| 国产视频内射| 啦啦啦免费观看视频1| 午夜福利成人在线免费观看| 亚洲片人在线观看| 在线播放国产精品三级| 宅男免费午夜| 啦啦啦韩国在线观看视频| 天堂√8在线中文| 精品久久久久久久久久免费视频| 亚洲精品美女久久久久99蜜臀| 国产成人a区在线观看| 亚洲精品色激情综合| 久久久久久久亚洲中文字幕 | 日本一二三区视频观看| 一区二区三区激情视频| 亚洲最大成人手机在线| 两个人的视频大全免费| 一进一出抽搐动态| 亚洲自拍偷在线| 国产单亲对白刺激| 国产精品香港三级国产av潘金莲| 天堂影院成人在线观看| 黄片小视频在线播放| 两个人看的免费小视频| 国产三级黄色录像| netflix在线观看网站| 国产精品女同一区二区软件 | 丁香六月欧美| 琪琪午夜伦伦电影理论片6080| 真人做人爱边吃奶动态| 欧美性感艳星| 757午夜福利合集在线观看| 亚洲欧美日韩卡通动漫| 国产精品亚洲一级av第二区| 欧美成狂野欧美在线观看| 国产高清视频在线播放一区| e午夜精品久久久久久久| 亚洲国产精品999在线| 婷婷丁香在线五月| 久久精品国产清高在天天线| 高潮久久久久久久久久久不卡| 老汉色av国产亚洲站长工具| 亚洲午夜理论影院| 久久6这里有精品| 国产乱人伦免费视频| 国产成人影院久久av| 三级国产精品欧美在线观看| 国产成人av教育| 亚洲乱码一区二区免费版| 亚洲av不卡在线观看| 久久亚洲精品不卡| 国产在视频线在精品| 欧美黄色片欧美黄色片| 中文字幕高清在线视频| 91在线精品国自产拍蜜月 | 欧美色视频一区免费| 人妻夜夜爽99麻豆av| 非洲黑人性xxxx精品又粗又长| 午夜福利18| 老司机福利观看| 欧美午夜高清在线| 免费一级毛片在线播放高清视频| 哪里可以看免费的av片| av视频在线观看入口| 亚洲成av人片在线播放无| 国产av不卡久久| 日韩欧美国产在线观看| 五月伊人婷婷丁香| 性欧美人与动物交配| 成人av一区二区三区在线看| 禁无遮挡网站| 国产老妇女一区| 亚洲av第一区精品v没综合| 一区二区三区国产精品乱码| 欧美黑人欧美精品刺激| 黄色丝袜av网址大全| 热99在线观看视频| 三级毛片av免费| 国产精品亚洲一级av第二区| 18禁在线播放成人免费| 国产色婷婷99| 中文字幕人成人乱码亚洲影| 黄色片一级片一级黄色片| 激情在线观看视频在线高清| 久久精品综合一区二区三区| 久久人妻av系列| 亚洲精品成人久久久久久| 一本一本综合久久| 成年版毛片免费区| 听说在线观看完整版免费高清| 日韩欧美在线二视频| 精品久久久久久,| 国产成人av教育| 麻豆久久精品国产亚洲av| 中国美女看黄片| 国产精品一区二区三区四区免费观看 | 看片在线看免费视频| 亚洲熟妇中文字幕五十中出| 国产 一区 欧美 日韩| avwww免费| 在线看三级毛片| 精品国产亚洲在线| 桃红色精品国产亚洲av| 精品不卡国产一区二区三区| 欧美一区二区亚洲| av黄色大香蕉| 男女做爰动态图高潮gif福利片| 91麻豆av在线| 色在线成人网| 亚洲人成伊人成综合网2020| 国产男靠女视频免费网站| 在线免费观看不下载黄p国产 | 亚洲精品乱码久久久v下载方式 | 很黄的视频免费| 久久久久亚洲av毛片大全| 亚洲国产色片| 国产一级毛片七仙女欲春2| 狠狠狠狠99中文字幕| 啦啦啦韩国在线观看视频| 日韩欧美一区二区三区在线观看| 精品久久久久久成人av| 国产精品香港三级国产av潘金莲| 国产精品久久久久久久久免 | 变态另类成人亚洲欧美熟女| 亚洲七黄色美女视频| 99国产精品一区二区三区| 国产精品女同一区二区软件 | 国产成人av教育| 欧美+亚洲+日韩+国产| 欧美日韩福利视频一区二区| 成人av一区二区三区在线看| 国产成人aa在线观看| 91字幕亚洲| 国产亚洲av嫩草精品影院| 中文字幕av在线有码专区| 精品乱码久久久久久99久播| 法律面前人人平等表现在哪些方面| 身体一侧抽搐| 99国产精品一区二区蜜桃av| 日韩国内少妇激情av| 亚洲av第一区精品v没综合| 99久久99久久久精品蜜桃| 日本免费a在线| 中文字幕高清在线视频| 99热只有精品国产| 日本一本二区三区精品| 亚洲精品色激情综合| 国产探花在线观看一区二区| 美女黄网站色视频| 国产在视频线在精品| 一a级毛片在线观看| 精品99又大又爽又粗少妇毛片 | 91九色精品人成在线观看| 51午夜福利影视在线观看| 亚洲无线观看免费| 欧美性猛交╳xxx乱大交人| 色老头精品视频在线观看| 岛国在线免费视频观看| 香蕉久久夜色| 免费大片18禁| 免费人成在线观看视频色| 久久精品综合一区二区三区| 91在线观看av| 免费大片18禁| 搞女人的毛片| 国产精品女同一区二区软件 | 亚洲人与动物交配视频| 两个人的视频大全免费| 国产成人aa在线观看| 男人舔奶头视频| 国产亚洲av嫩草精品影院| 全区人妻精品视频| 亚洲成av人片免费观看| 两个人的视频大全免费| 99精品久久久久人妻精品| 午夜福利在线在线| 色视频www国产| 精品国产三级普通话版| 国产又黄又爽又无遮挡在线| 91麻豆精品激情在线观看国产| 国产在线精品亚洲第一网站| 国产精品日韩av在线免费观看| 欧美+日韩+精品| 黄色视频,在线免费观看| 无限看片的www在线观看| 18美女黄网站色大片免费观看| 动漫黄色视频在线观看| 18美女黄网站色大片免费观看| 久久精品91蜜桃| 久久久久性生活片| 91av网一区二区| 99热这里只有精品一区| 婷婷六月久久综合丁香| 欧美一级a爱片免费观看看| 亚洲国产欧美人成| 亚洲午夜理论影院| 欧美在线一区亚洲| 国产午夜福利久久久久久| 熟女人妻精品中文字幕| 日本与韩国留学比较| 99久久成人亚洲精品观看| 久久久久久久久久黄片| 国产探花在线观看一区二区| 国产一区二区在线av高清观看| 免费av观看视频| 日韩av在线大香蕉| av在线蜜桃| 亚洲成av人片免费观看| 两个人的视频大全免费| 国产美女午夜福利| 日韩欧美国产一区二区入口| 国产97色在线日韩免费| 欧美av亚洲av综合av国产av| 国产97色在线日韩免费| 亚洲精品一区av在线观看| 老司机福利观看| 亚洲激情在线av| 99热这里只有是精品50| 亚洲内射少妇av| 在线国产一区二区在线| 亚洲欧美一区二区三区黑人| 黄色日韩在线| 亚洲精品久久国产高清桃花| 国产极品精品免费视频能看的| 国产视频内射| 亚洲美女黄片视频| 国产成人啪精品午夜网站| 久久亚洲精品不卡|