• <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ā)展對策
    欧美成人免费av一区二区三区| www.999成人在线观看| 午夜亚洲福利在线播放| 色视频www国产| 怎么达到女性高潮| 99视频精品全部免费 在线| 亚洲成人精品中文字幕电影| 极品教师在线免费播放| 久久久久性生活片| 欧美日本视频| 日韩成人在线观看一区二区三区| 亚洲av一区综合| 简卡轻食公司| 窝窝影院91人妻| 午夜两性在线视频| 亚洲,欧美精品.| 日本a在线网址| 国产日本99.免费观看| 美女被艹到高潮喷水动态| 国产一级毛片七仙女欲春2| 黄色女人牲交| 亚洲精华国产精华精| 精品午夜福利视频在线观看一区| 757午夜福利合集在线观看| 无遮挡黄片免费观看| 亚洲最大成人av| 亚洲成人久久爱视频| 尤物成人国产欧美一区二区三区| 女同久久另类99精品国产91| 亚洲精品色激情综合| 18美女黄网站色大片免费观看| 亚洲av不卡在线观看| 好看av亚洲va欧美ⅴa在| 91字幕亚洲| 色哟哟·www| 1024手机看黄色片| 小蜜桃在线观看免费完整版高清| 大型黄色视频在线免费观看| 我要看日韩黄色一级片| 99热这里只有是精品在线观看 | 我要搜黄色片| 色av中文字幕| 亚洲国产精品合色在线| 久久99热这里只有精品18| 亚洲国产精品成人综合色| 女人被狂操c到高潮| 好看av亚洲va欧美ⅴa在| 亚洲精品一区av在线观看| 露出奶头的视频| 国产麻豆成人av免费视频| 欧美最黄视频在线播放免费| 午夜福利高清视频| 黄色视频,在线免费观看| 三级毛片av免费| 欧美最黄视频在线播放免费| 日韩人妻高清精品专区| 精品久久国产蜜桃| 国产成人av教育| 亚洲av中文字字幕乱码综合| 日本五十路高清| 一个人免费在线观看的高清视频| 亚洲第一欧美日韩一区二区三区| 黄色日韩在线| 天堂动漫精品| 欧美色视频一区免费| 啦啦啦观看免费观看视频高清| 动漫黄色视频在线观看| 成人特级av手机在线观看| 亚洲成人久久爱视频| 国产在线男女| 亚洲国产精品999在线| 少妇裸体淫交视频免费看高清| 色视频www国产| 国产免费男女视频| 国产一区二区在线观看日韩| 最后的刺客免费高清国语| 国产又黄又爽又无遮挡在线| 午夜福利高清视频| 男女之事视频高清在线观看| 亚洲第一欧美日韩一区二区三区| 国产精品98久久久久久宅男小说| 久久久久久久久大av| 十八禁人妻一区二区| 国产综合懂色| ponron亚洲| 欧美一区二区国产精品久久精品| 性欧美人与动物交配| 欧美激情国产日韩精品一区| 美女大奶头视频| 搡老妇女老女人老熟妇| av黄色大香蕉| 大型黄色视频在线免费观看| 亚洲欧美日韩卡通动漫| 一本精品99久久精品77| 亚洲国产高清在线一区二区三| 日本精品一区二区三区蜜桃| 国产日本99.免费观看| 18+在线观看网站| 最好的美女福利视频网| 亚洲黑人精品在线| 中文字幕人妻熟人妻熟丝袜美| 亚洲电影在线观看av| 一二三四社区在线视频社区8| 国产欧美日韩一区二区精品| 成人永久免费在线观看视频| 欧美bdsm另类| 欧美成狂野欧美在线观看| 看免费av毛片| 最好的美女福利视频网| 在线观看午夜福利视频| a级一级毛片免费在线观看| 嫩草影视91久久| 99国产综合亚洲精品| 女人被狂操c到高潮| 欧美不卡视频在线免费观看| 久久精品影院6| 国产麻豆成人av免费视频| 91狼人影院| 精品久久久久久久久亚洲 | 一卡2卡三卡四卡精品乱码亚洲| 免费无遮挡裸体视频| 国产精品,欧美在线| 国产淫片久久久久久久久 | 欧美精品国产亚洲| 麻豆av噜噜一区二区三区| 少妇熟女aⅴ在线视频| 精品国产三级普通话版| 一级作爱视频免费观看| 欧美成人免费av一区二区三区| 波野结衣二区三区在线| 亚洲自偷自拍三级| 欧美中文日本在线观看视频| 黄色配什么色好看| 国产国拍精品亚洲av在线观看| 美女 人体艺术 gogo| 好看av亚洲va欧美ⅴa在| 一a级毛片在线观看| 老鸭窝网址在线观看| 最近最新中文字幕大全电影3| 身体一侧抽搐| 搡老妇女老女人老熟妇| 怎么达到女性高潮| 国产不卡一卡二| 在线观看午夜福利视频| 99久久精品一区二区三区| 亚洲av成人不卡在线观看播放网| 狠狠狠狠99中文字幕| 蜜桃久久精品国产亚洲av| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久久久久黄片| 无人区码免费观看不卡| 亚洲经典国产精华液单 | 真实男女啪啪啪动态图| 亚洲av成人不卡在线观看播放网| 亚洲av熟女| 中文字幕人妻熟人妻熟丝袜美| 国产爱豆传媒在线观看| 国产精品女同一区二区软件 | 国语自产精品视频在线第100页| 亚洲,欧美精品.| 午夜激情欧美在线| 一个人看的www免费观看视频| 97超视频在线观看视频| 久久99热6这里只有精品| 激情在线观看视频在线高清| 国产三级黄色录像| 夜夜夜夜夜久久久久| 白带黄色成豆腐渣| 成年女人永久免费观看视频| 蜜桃久久精品国产亚洲av| 色哟哟哟哟哟哟| 国产精品自产拍在线观看55亚洲| 别揉我奶头 嗯啊视频| 欧美日韩黄片免| 噜噜噜噜噜久久久久久91| 欧美极品一区二区三区四区| 亚洲精品在线美女| 日本a在线网址| 国产又黄又爽又无遮挡在线| 97人妻精品一区二区三区麻豆| 日韩高清综合在线| 亚洲,欧美,日韩| 久久中文看片网| 热99re8久久精品国产| 757午夜福利合集在线观看| 亚州av有码| 桃色一区二区三区在线观看| 国内精品美女久久久久久| 听说在线观看完整版免费高清| 亚洲avbb在线观看| 色哟哟哟哟哟哟| 在线观看美女被高潮喷水网站 | 中文字幕人成人乱码亚洲影| 给我免费播放毛片高清在线观看| 日本熟妇午夜| 欧美性猛交黑人性爽| 婷婷丁香在线五月| 又粗又爽又猛毛片免费看| 久久久久久久精品吃奶| 人人妻人人看人人澡| 久久亚洲精品不卡| 久久99热6这里只有精品| 欧美日韩综合久久久久久 | 亚洲av熟女| 免费观看精品视频网站| 国内揄拍国产精品人妻在线| 狂野欧美白嫩少妇大欣赏| 99久久九九国产精品国产免费| 九色成人免费人妻av| 精品午夜福利视频在线观看一区| 国产精品久久久久久亚洲av鲁大| 亚洲精品影视一区二区三区av| 如何舔出高潮| 成年女人看的毛片在线观看| 欧美日韩乱码在线| 欧美极品一区二区三区四区| 国产精品一区二区性色av| 日本免费a在线| 99热这里只有是精品50| eeuss影院久久| 欧美色视频一区免费| 欧美激情久久久久久爽电影| 成年女人永久免费观看视频| 国产成人a区在线观看| 美女大奶头视频| 国产淫片久久久久久久久 | 99国产综合亚洲精品| 老女人水多毛片| 亚洲精品一区av在线观看| 桃红色精品国产亚洲av| 日韩精品中文字幕看吧| 99国产综合亚洲精品| 国产在线精品亚洲第一网站| 久久精品国产亚洲av香蕉五月| 18禁黄网站禁片午夜丰满| 亚洲经典国产精华液单 | 亚洲色图av天堂| 麻豆国产97在线/欧美| 国产久久久一区二区三区| 嫩草影视91久久| 俺也久久电影网| 9191精品国产免费久久| 一二三四社区在线视频社区8| 亚洲av电影不卡..在线观看| 欧美日本亚洲视频在线播放| 国产麻豆成人av免费视频| 精品不卡国产一区二区三区| 欧美zozozo另类| 麻豆一二三区av精品| 九九久久精品国产亚洲av麻豆| 很黄的视频免费| 亚洲电影在线观看av| а√天堂www在线а√下载| 国产探花在线观看一区二区| 亚洲av美国av| 亚洲av电影不卡..在线观看| 亚洲国产精品sss在线观看| 毛片女人毛片| 国产精品乱码一区二三区的特点| 18禁在线播放成人免费| 天堂av国产一区二区熟女人妻| 免费av毛片视频| 最新在线观看一区二区三区| 日本成人三级电影网站| 欧美激情在线99| 久久久久久久久中文| 一进一出好大好爽视频| 老司机福利观看| 老熟妇仑乱视频hdxx| 一区二区三区免费毛片| 黄色配什么色好看| 超碰av人人做人人爽久久| 免费在线观看日本一区| 国产私拍福利视频在线观看| 桃色一区二区三区在线观看| 大型黄色视频在线免费观看| 成人性生交大片免费视频hd| 免费在线观看影片大全网站| 欧美成狂野欧美在线观看| 国产黄片美女视频| 亚洲欧美日韩卡通动漫| 国产亚洲精品av在线| 亚州av有码| 欧美高清性xxxxhd video| 亚洲乱码一区二区免费版| 久久天躁狠狠躁夜夜2o2o| 91九色精品人成在线观看| 欧美黑人欧美精品刺激| 一本精品99久久精品77| 免费av毛片视频| 国产v大片淫在线免费观看| 91字幕亚洲| 欧美乱妇无乱码| 一个人看的www免费观看视频| 免费大片18禁| 赤兔流量卡办理| 亚洲av成人av| 日本一本二区三区精品| 久久精品久久久久久噜噜老黄 | 国产伦精品一区二区三区四那| 高清日韩中文字幕在线| 级片在线观看| 国产精品久久久久久精品电影| 中出人妻视频一区二区| 日韩欧美精品v在线| 欧美色视频一区免费| 熟女电影av网| 久久天躁狠狠躁夜夜2o2o| 特级一级黄色大片| 听说在线观看完整版免费高清| 99热这里只有是精品50| 国产私拍福利视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲自拍偷在线| bbb黄色大片| 欧美日韩黄片免| 在线观看午夜福利视频| 熟妇人妻久久中文字幕3abv| 一个人免费在线观看的高清视频| 可以在线观看毛片的网站| 高清日韩中文字幕在线| 老女人水多毛片| 亚洲国产精品sss在线观看| 18禁黄网站禁片午夜丰满| 亚洲精品在线美女| 午夜免费激情av| 婷婷精品国产亚洲av| 亚洲最大成人av| 真人做人爱边吃奶动态| 90打野战视频偷拍视频| 免费人成在线观看视频色| 亚洲自偷自拍三级| 啦啦啦观看免费观看视频高清| 国产精品,欧美在线| 在线播放国产精品三级| 国产伦精品一区二区三区视频9| 国产精品伦人一区二区| 国产真实伦视频高清在线观看 | 国产亚洲精品综合一区在线观看| 欧美性猛交黑人性爽| av专区在线播放| 亚洲国产欧洲综合997久久,| 国内精品久久久久精免费| 国产精品,欧美在线| 级片在线观看| 亚州av有码| 国产大屁股一区二区在线视频| 精品国产三级普通话版| 欧洲精品卡2卡3卡4卡5卡区| 热99在线观看视频| 亚洲第一电影网av| 免费av毛片视频| 丰满人妻熟妇乱又伦精品不卡| 久久精品夜夜夜夜夜久久蜜豆| 亚洲第一欧美日韩一区二区三区| 午夜福利欧美成人| 中文字幕免费在线视频6| 精品人妻视频免费看| 99热这里只有是精品在线观看 | 欧洲精品卡2卡3卡4卡5卡区| 免费av毛片视频| www.熟女人妻精品国产| 国产欧美日韩一区二区三| 日本在线视频免费播放| 国产精品久久久久久精品电影| 色尼玛亚洲综合影院| 最新中文字幕久久久久| 亚洲不卡免费看| 亚洲美女视频黄频| 国产av一区在线观看免费| 精品国产三级普通话版| 日韩欧美国产在线观看| 亚洲18禁久久av| 高潮久久久久久久久久久不卡| 丰满人妻一区二区三区视频av| 在线播放国产精品三级| 久久精品国产清高在天天线| 夜夜躁狠狠躁天天躁| 精品午夜福利在线看| 亚洲人成网站高清观看| 少妇的逼好多水| 久久久久久久亚洲中文字幕 | 成人精品一区二区免费| 九九热线精品视视频播放| 少妇的逼水好多| 久久久久久久午夜电影| 久久久久免费精品人妻一区二区| 少妇裸体淫交视频免费看高清| 欧美绝顶高潮抽搐喷水| 亚洲黑人精品在线| 熟女电影av网| 99精品久久久久人妻精品| 每晚都被弄得嗷嗷叫到高潮| 真人一进一出gif抽搐免费| h日本视频在线播放| 啪啪无遮挡十八禁网站| 欧美区成人在线视频| 国产精华一区二区三区| 日本黄大片高清| 真人一进一出gif抽搐免费| 91久久精品电影网| 非洲黑人性xxxx精品又粗又长| 精品人妻1区二区| 精品无人区乱码1区二区| 亚洲精品一区av在线观看| 欧美黄色淫秽网站| eeuss影院久久| 日本三级黄在线观看| www.色视频.com| 午夜久久久久精精品| 黄片小视频在线播放| 日本一本二区三区精品| 国产精品日韩av在线免费观看| 真人做人爱边吃奶动态| 日本免费a在线| 99久久九九国产精品国产免费| 久久天躁狠狠躁夜夜2o2o| 在线天堂最新版资源| 亚洲成av人片免费观看| 日韩人妻高清精品专区| 性欧美人与动物交配| 精品久久久久久久久亚洲 | 变态另类成人亚洲欧美熟女| 国产精品嫩草影院av在线观看 | 精华霜和精华液先用哪个| 亚洲美女搞黄在线观看 | 日本 av在线| 高清在线国产一区| 国产高清视频在线观看网站| 18+在线观看网站| 男女之事视频高清在线观看| av天堂中文字幕网| 真人一进一出gif抽搐免费| 91字幕亚洲| 最新在线观看一区二区三区| 搡女人真爽免费视频火全软件 | 好男人在线观看高清免费视频| 精品人妻偷拍中文字幕| 成年免费大片在线观看| 搡女人真爽免费视频火全软件 | 久9热在线精品视频| 91午夜精品亚洲一区二区三区 | 黄色配什么色好看| 嫁个100分男人电影在线观看| 欧美成人a在线观看| 欧美日韩综合久久久久久 | 国产成人欧美在线观看| 午夜a级毛片| 级片在线观看| 琪琪午夜伦伦电影理论片6080| 久久久久免费精品人妻一区二区| 毛片一级片免费看久久久久 | 久久精品人妻少妇| 老司机午夜福利在线观看视频| 九九久久精品国产亚洲av麻豆| 亚洲av成人精品一区久久| 丰满乱子伦码专区| 欧美激情久久久久久爽电影| 国模一区二区三区四区视频| 在线天堂最新版资源| 别揉我奶头~嗯~啊~动态视频| 国产av在哪里看| 精品午夜福利视频在线观看一区| 亚洲,欧美精品.| 精品久久久久久久人妻蜜臀av| 欧美日韩乱码在线| 久久香蕉精品热| 一级作爱视频免费观看| 成人三级黄色视频| 久久国产精品人妻蜜桃| 噜噜噜噜噜久久久久久91| 国内精品一区二区在线观看| 男插女下体视频免费在线播放| 日本免费一区二区三区高清不卡| 亚洲一区高清亚洲精品| 黄色女人牲交| 免费一级毛片在线播放高清视频| 天美传媒精品一区二区| 国产精品国产高清国产av| 午夜福利视频1000在线观看| 午夜免费成人在线视频| 久久伊人香网站| 亚洲18禁久久av| 韩国av一区二区三区四区| 91在线精品国自产拍蜜月| 日本成人三级电影网站| 无人区码免费观看不卡| 一区福利在线观看| 欧美潮喷喷水| 亚洲七黄色美女视频| 女同久久另类99精品国产91| 午夜精品在线福利| 日韩成人在线观看一区二区三区| 老女人水多毛片| 亚洲av电影在线进入| 国产精品国产高清国产av| 青草久久国产| 精品人妻1区二区| 亚洲国产精品久久男人天堂| 国产成人福利小说| 日韩精品中文字幕看吧| 国产淫片久久久久久久久 | 热99re8久久精品国产| 少妇被粗大猛烈的视频| 精品久久久久久久久久免费视频| av在线天堂中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 精品熟女少妇八av免费久了| 午夜福利欧美成人| 99精品在免费线老司机午夜| 性色av乱码一区二区三区2| 99视频精品全部免费 在线| 狂野欧美白嫩少妇大欣赏| 一a级毛片在线观看| 中文字幕久久专区| 天美传媒精品一区二区| 999久久久精品免费观看国产| 久久久久性生活片| 国产精品一及| 国产午夜精品论理片| 国产成人影院久久av| 国产三级中文精品| 久久精品国产亚洲av天美| 国产三级黄色录像| 国语自产精品视频在线第100页| 久99久视频精品免费| 少妇丰满av| 国产av麻豆久久久久久久| 欧美激情在线99| 91久久精品电影网| 校园春色视频在线观看| 天堂√8在线中文| 日韩国内少妇激情av| 男插女下体视频免费在线播放| av在线观看视频网站免费| 久久久久精品国产欧美久久久| or卡值多少钱| 在线观看66精品国产| 亚洲人与动物交配视频| 国产蜜桃级精品一区二区三区| a级一级毛片免费在线观看| 啪啪无遮挡十八禁网站| 中文字幕av成人在线电影| 免费黄网站久久成人精品 | 搡老熟女国产l中国老女人| 国产成人影院久久av| 2021天堂中文幕一二区在线观| 国产aⅴ精品一区二区三区波| 久久精品综合一区二区三区| 国语自产精品视频在线第100页| 51国产日韩欧美| 丰满的人妻完整版| 亚洲午夜理论影院| 亚洲五月婷婷丁香| 午夜福利18| 国内精品久久久久精免费| 亚洲欧美日韩卡通动漫| 大型黄色视频在线免费观看| www.熟女人妻精品国产| 美女xxoo啪啪120秒动态图 | 在线看三级毛片| 丰满人妻一区二区三区视频av| 特级一级黄色大片| av女优亚洲男人天堂| 在线观看一区二区三区| 看片在线看免费视频| 毛片一级片免费看久久久久 | 男女做爰动态图高潮gif福利片| 久久久久久久久中文| 日韩高清综合在线| 少妇人妻一区二区三区视频| 免费人成在线观看视频色| 欧美黑人巨大hd| 日本五十路高清| 99热这里只有是精品50| 99国产精品一区二区蜜桃av| 亚洲专区中文字幕在线| 久久久久精品国产欧美久久久| 免费大片18禁| 亚洲成人久久性| 日本三级黄在线观看| 极品教师在线免费播放| 精品一区二区三区人妻视频| 老司机福利观看| 欧美性猛交黑人性爽| 国产精品久久久久久久电影| 蜜桃亚洲精品一区二区三区| 国产一区二区在线av高清观看| 波野结衣二区三区在线| 欧美丝袜亚洲另类 | 最近在线观看免费完整版| 男女下面进入的视频免费午夜| 18+在线观看网站| 欧美乱妇无乱码| 久久亚洲精品不卡| 好看av亚洲va欧美ⅴa在| 1024手机看黄色片| 99久久精品热视频| 国产伦在线观看视频一区| 国产精品嫩草影院av在线观看 | 成年女人永久免费观看视频| 在线观看免费视频日本深夜| 老司机午夜福利在线观看视频| 国产精品三级大全| 久久香蕉精品热| 色综合欧美亚洲国产小说| 国产一区二区亚洲精品在线观看| 国产亚洲精品久久久久久毛片| 色综合欧美亚洲国产小说| 亚洲av电影在线进入| 身体一侧抽搐| 国内久久婷婷六月综合欲色啪| 欧美性感艳星| 欧美一区二区精品小视频在线| 国产大屁股一区二区在线视频|