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

    舟山近海域懸浮泥沙濃度遙感反演研究

    2016-03-16 06:29:49錢揚(yáng)風(fēng)郭碧云羅鳳云韓錄維
    安徽農(nóng)業(yè)科學(xué) 2016年3期

    錢揚(yáng)風(fēng), 郭碧云, 黃 東, 羅鳳云, 韓錄維

    (浙江海洋學(xué)院, 浙江舟山 31600)

    ?

    舟山近海域懸浮泥沙濃度遙感反演研究

    錢揚(yáng)風(fēng), 郭碧云*, 黃 東, 羅鳳云, 韓錄維

    (浙江海洋學(xué)院, 浙江舟山 31600)

    摘要[目的]采用衛(wèi)星遙感方法反演舟山海域懸浮泥沙的分布,并進(jìn)一步探索適合的反演模型,為以后舟山海域懸浮泥沙遙感研究提供參考依據(jù)。[方法]利用實(shí)測數(shù)據(jù)、遙感數(shù)據(jù)和實(shí)驗(yàn)室分析的方法,對(duì)舟山海域不同水體進(jìn)行光譜分析,建立水體遙感反射率與懸浮泥沙濃度之間3種不同的反演模型。[結(jié)果]遙感參數(shù)反演模型的反演效果最好;對(duì)數(shù)反演模型次之,對(duì)數(shù)模型下TM4波段的相關(guān)系數(shù)最大,TM3波段的相關(guān)系數(shù)最??;線性反演模型最差,線性模型下TM2波段的相關(guān)系數(shù)最大,TM3波段的相關(guān)系數(shù)最小。[結(jié)論]在舟山海域采用遙感參數(shù)模型可以獲得較高的反演精度。

    關(guān)鍵詞舟山海域;懸浮泥沙;遙感反射率;反演模型

    Remote Sensing Retrieval of Suspended Sediment Concentration in Near Sea Area of Zhoushan

    QIAN Yang-feng, GUO Bi-yun*, HUANG Dong et al(Zhejiang Ocean University, Zhoushan, Zhejiang 31600)

    Abstract[Objective] The distribution of suspended sediment in Zhoushan sea area was retrieved by using satellite remote sensing method, and the suitable inversion model was further explored, which will provide reference for remote sensing study of suspended sediment in Zhoushan sea area.[Method] Using field data, remote sensing data and laboratory analysis method, spectral analysis of different water bodies in Zhoushan sea area was carried out, and three kinds of inversion models between water body remote sensing reflectance and suspended sediment concentration were established.[Result] The results showed that the inversion effect of remote sensing parameters model was the best.The second was Logarithmic inversion model.The correlation coefficient was the largest in TM4 band and the smallest in TM3.The linear inversion model was the worst.The correlation coefficient was the largest in TM2 band and the smallest in TM3.[Conclusion] The higher inversion accuracy can be obtained by using the remote sensing parameters model in Zhoushan sea area.

    Key wordsZhoushan sea area; Suspended sediment; Remote sensing reflectance; Inversion model

    水體中的懸浮泥沙含量是重要的水質(zhì)參數(shù)之一。懸浮泥沙濃度對(duì)研究水流挾沙能力和河口水流泥沙運(yùn)動(dòng)規(guī)律等具有重要作用,對(duì)研究海岸帶的水質(zhì)、生態(tài)、環(huán)境以及海岸工程、港口建設(shè)等具有重要意義[1]。傳統(tǒng)懸浮泥沙調(diào)查方法是利用船只在相關(guān)海域進(jìn)行逐點(diǎn)采樣分析,速度慢,周期長,只能獲得在時(shí)間和空間分布上都很離散的少量數(shù)據(jù),而且這種傳統(tǒng)的調(diào)查方法存在著很多實(shí)際操作上的困難。隨著衛(wèi)星遙感技術(shù)的發(fā)展,衛(wèi)星遙感監(jiān)測逐漸從定性發(fā)展到定量,遙感精度逐漸提高。利用衛(wèi)星遙感大范圍、周期短的優(yōu)點(diǎn),可以實(shí)現(xiàn)大面積水域懸浮泥沙動(dòng)態(tài)、連續(xù)的觀測,為近海海域懸浮泥沙研究提供了新的有效手段。

    關(guān)于我國海區(qū)懸浮泥沙遙感的研究較多,但比較集中在一些重點(diǎn)河口和港灣,如長江口、珠江口、杭州灣、渤海灣等,對(duì)舟山海域懸浮泥沙遙感研究較少,建立的反演模型不夠完善。舟山海域緊鄰長江入??冢瑣u嶼星羅棋布,近岸海流運(yùn)動(dòng)復(fù)雜,影響泥沙運(yùn)動(dòng)的水動(dòng)力機(jī)制多樣,其特殊的地理位置、島嶼的緊密分布以及人類活動(dòng)對(duì)舟山海域懸浮泥沙分布有較大影響。對(duì)舟山海域懸浮泥沙分布的研究可以為舟山未來的海洋生態(tài)建設(shè)、海洋工程建設(shè)與海洋資源的開發(fā)利用提供參考依據(jù)。筆者在前人研究的基礎(chǔ)上,探索適合舟山海域的懸浮泥沙反演模型。

    1研究區(qū)概況

    舟山群島位于浙江省東北部,長江口南側(cè),杭州灣外緣的東海洋面上,地理坐標(biāo)為121°31′~123°25′E、29°32′~31°04′N,總面積22 000 km2,其中海域面積20 800 km2,陸域面積1 440 km2。舟山群島由1 390個(gè)大小島嶼組成,是全國第一大群島(圖1)。舟山群島為海島丘陵地貌,群島呈西南—東北走向,南部大島較多,海拔較高,排列密集;北部以小島為主,地勢漸低,分布稀疏。舟山群島屬北亞熱帶南緣季風(fēng)海洋型氣候,四季分明,冬暖夏涼,年溫適中,光熱充裕,蒸發(fā)量大,降雨偏少[2]。2011年6月30日,國務(wù)院正式批準(zhǔn)設(shè)立浙江舟山群島新區(qū),也是首個(gè)以海洋經(jīng)濟(jì)為主題的國家級(jí)新區(qū)。研究舟山近海域懸沙分布對(duì)舟山建設(shè)國家級(jí)新區(qū),發(fā)展海洋經(jīng)濟(jì)具有重要意義。

    2數(shù)據(jù)來源與處理

    2.1實(shí)測數(shù)據(jù) 該研究采用的測量數(shù)據(jù)是2015年7月2日搭乘浙江海洋學(xué)院浙???號(hào)科考船出海測量所得,觀測點(diǎn)共6個(gè)(圖2)。采用合肥儀思特光電技術(shù)有限公司生產(chǎn)的ISI921VF系列野外地物光譜輻射計(jì)對(duì)不同海域水體的光譜值進(jìn)行測量,獲得該海域水體的各種光譜數(shù)據(jù),同時(shí)利用采水器獲得不同觀測點(diǎn)位置的水樣。

    圖2 測量點(diǎn)位置示意Fig.2 Position of observation point

    2.2數(shù)據(jù)處理

    2.2.1等效遙感反射率計(jì)算。建立懸浮泥沙實(shí)際濃度與遙感反射率之間的關(guān)系,需將各種實(shí)地測量的光譜值轉(zhuǎn)化為遙感反射率。由遙感反射率Rs的定義可知:

    Rs=Lw/Es

    (1)

    式中,Lw為離水輻亮度;Es為水面入射輻照度。這2個(gè)值都不是直接的測量值,需對(duì)測量數(shù)據(jù)處理分析才能得到。海水的離水輻亮度Lw與海水輻射率Lseawater、天空漫散射光信號(hào)Lsky的關(guān)系為:

    Lseawater= Lw+ rLsky+ Lg

    (2)

    式中,海水輻射率Lseawater、天空漫散射光信號(hào)Lsky是直接的測量值,r是水-氣界面反射率, 也稱為菲涅耳反射系數(shù)。Lg為太陽的鏡面輻射率,通常情況下這一項(xiàng)很難測量,但如果采用適當(dāng)?shù)臏y量方位角和儀器觀測角度,此項(xiàng)的作用可以避免。在避開太陽直接反射、忽略或避開白帽的情況下,可以忽略太陽鏡面反射率Lg。式(2)化簡可得到海水離水輻射率:

    Lw=Lseawater- rLsky

    (3)

    水面入射輻照度Es與標(biāo)準(zhǔn)板的反射信號(hào)Lp及標(biāo)準(zhǔn)板的反射率ρP有關(guān),兩者均是直接的測量值,計(jì)算公式為:

    Es= π Lp/ρP

    (4)

    推算遙感反射率的目的是建立其與懸浮泥沙濃度的相關(guān)關(guān)系,但是反演泥沙分布使用的是TM衛(wèi)星影像,所以在進(jìn)行分析時(shí),要實(shí)現(xiàn)光譜儀測得的光譜值與TM數(shù)據(jù)之間的等效遙感反射率的計(jì)算[3]。TM衛(wèi)星影像指美國陸地衛(wèi)星Landsat4-5號(hào)專題制圖儀(thematicmapper)所獲取的多波段掃描影像,分為7個(gè)波段:TM1(450~520nm)藍(lán)光波段、TM2(520~600nm)綠光波段、TM3(630~0.69nm)紅光波段、TM4(760~900nm)近紅外波段、TM5(1 550~1 750nm)中紅外波段、TM6(1 040~1 250nm)熱紅外波段、TM7(2 080~2 350nm)中遠(yuǎn)紅外波段。該研究選用前4個(gè)波段,將遙感反射率等效計(jì)算為TM影像上1~4波段的反射率,結(jié)果見表1。

    表1 等效計(jì)算各TM波段下對(duì)應(yīng)的遙感反射率

    2.2.2懸浮泥沙濃度測量。在實(shí)驗(yàn)室采用過濾稱重法對(duì)在不同觀測點(diǎn)采集的水樣泥沙濃度進(jìn)行測量。先定量量取200 mL搖勻水樣并置于漏斗中,用抽濾機(jī)將漏斗中的海水快速抽離,得到覆蓋懸浮泥沙的濾紙片。然后將覆蓋懸浮泥沙的濾紙片中的水分烘干,用電子天平稱量得到烘干后濾紙與懸浮泥沙總質(zhì)量,減去提前測出的濾紙片質(zhì)量即可得到懸浮泥沙質(zhì)量,最后除以水樣體積即為懸浮泥沙含量。懸浮泥沙含量計(jì)算公式為:

    懸沙含量=懸沙質(zhì)量/水樣體積

    (5)

    2015年7月2日各取樣點(diǎn)表層含沙量見表2。

    3結(jié)果與分析

    3.1懸浮泥沙水體的光譜特征分析運(yùn)用MATLAB繪制出遙感反射率光譜曲線,結(jié)果見圖3。不同波段電磁波對(duì)海水的穿透能力存在差異,不同含沙量水體對(duì)太陽輻射光譜的反射率明顯不同。從圖3可以看出,隨著水體含沙量的增加,各波段的反射率都普遍增大,但其增幅不同。懸浮泥沙光譜反射率普遍具有雙峰特征,第一反射峰位置在500~600 nm波段,第二反射峰位置在750~850 nm波段。當(dāng)水體中懸浮泥沙含量增加時(shí),波譜上的反射峰由短波向長波方向移動(dòng),即所謂的“紅移現(xiàn)象”。隨著水體中懸浮泥沙濃度改變,各波段上反射率的改變存在差異。有的波段改變大,說明該波段對(duì)懸浮泥沙變化敏感,懸浮泥沙信息比較豐富;而有的波段基本不改變,即該波段基本不含有懸浮泥沙信息,因此在反演過程中要選擇對(duì)懸浮泥沙濃度敏感的波段[4-5]。對(duì)此,考慮建立不同的反演模型來探求舟山海域懸浮泥沙的適合反演模型。

    表2不同觀測點(diǎn)表層水體含沙量

    Table 2The sediment concentration of surface water in various observation sites

    樣點(diǎn)標(biāo)號(hào)Locationlabel濃度Concentrationmg/L經(jīng)緯度Latitudeandlongitude1112.75122°20.463'E、29°53.908'N2108.25122°11.680'E、29°55.099'N3108.00122°15.191'E、29°51.302'N495.75122°18.132'E、29°54.248'N578.50122°17.880'E、29°50.638'N678.25122°13.547'E、29°52.513'N

    圖3 不同濃度懸浮泥沙水體反射波譜曲線Fig.3 Reflection spectrum curve of suspended sediment in different concentration

    3.2懸浮泥沙反演模型的建立

    3.2.1線性模型。不同TM波段對(duì)懸浮泥沙的敏感程度不同,得到的反演模型精度也不同,據(jù)此建立不同TM單波段與遙感反射率的線性相關(guān)模型。首先在整個(gè)波段,如圖3所示,每個(gè)固定波長下存在6個(gè)觀測點(diǎn)遙感反射率,將這6個(gè)觀測點(diǎn)的遙感反射率與其對(duì)應(yīng)的6個(gè)觀測點(diǎn)的懸浮泥沙濃度比較,對(duì)遙感反射率與懸浮泥沙濃度兩組數(shù)作線性回歸分析得出該波長下的線性相關(guān)系數(shù),進(jìn)而推廣至整個(gè)波段,計(jì)算出每個(gè)波長下采樣點(diǎn)遙感反射率與懸浮泥沙濃度的線性相關(guān)系數(shù),結(jié)果見圖4。

    圖4 波長與線性相關(guān)系數(shù)對(duì)應(yīng)曲線Fig.4 The corresponding curves of wavelength and linear correlation coefficient

    從圖4可以看出,在400~900 nm波段下相關(guān)系數(shù)均非常大,平均值在0.9以上。在900 nm之后的波段出現(xiàn)不規(guī)則大幅波動(dòng),平均相關(guān)系數(shù)不大,故不考慮900 nm之后的TM5、TM6、TM7波段。分別在TM1、TM2、TM3、TM4波段上,對(duì)不同TM波段等效遙感反射率與泥沙濃度做線性相關(guān)分析,并得到線性回歸方程與相關(guān)系數(shù)(表3)。從表3可以看出,TM2波段下的相關(guān)系數(shù)最大,TM1與TM4波段次之,TM3波段下的相關(guān)系數(shù)最小。

    表3 線性回歸方程與相關(guān)系數(shù)

    注:參數(shù)X代表遙感反射率;Y代表懸浮泥沙濃度。

    Note: ParameterXstands for remote sensing reflectance;Ydenotes suspended sediment concentration.

    3.2.2對(duì)數(shù)模型。2000年朱立俊等[6]根據(jù)不同海區(qū)的陸地衛(wèi)星資料及地面同步實(shí)測數(shù)據(jù),建立了水體懸浮泥沙對(duì)數(shù)反演模型,其表達(dá)式為:

    X= - 5 714 + 0.595 2 log10S

    (6)

    式中,X為陸地衛(wèi)星TM3波段與TM2波段(或MSS5波段與MSS4波段)反射率比值;S為水體懸浮泥沙濃度[6-8]。參考建立此對(duì)數(shù)模型的方法,建立適合舟山海域的對(duì)數(shù)模型。采用“3.2.1”同樣的數(shù)據(jù)分析方法,得出整個(gè)波長下遙感反射率與懸浮泥沙濃度的對(duì)數(shù)相關(guān)系數(shù),結(jié)果見圖5。從圖5中可以得出與表3線性模型相似的結(jié)論,在400~920 nm波段下相關(guān)系數(shù)均非常高,平均相關(guān)系數(shù)在0.9以上。在920 nm之后的波段出現(xiàn)不規(guī)則大幅波動(dòng),不適宜取TM5、TM6、TM7波段,故略去。據(jù)此在TM1、TM2、TM3、TM4波段建立回歸方程,得到對(duì)應(yīng)的相關(guān)系數(shù)(表4)。從表4可以看出,總體上4個(gè)波段相關(guān)系數(shù)都比較大,不同波段之間差別不明顯。其中TM4波段下的相關(guān)系數(shù)最大,TM1、TM2波段次之,TM3波段下的相關(guān)系數(shù)最小。

    圖5 波長與對(duì)數(shù)相關(guān)系數(shù)對(duì)應(yīng)曲線Fig.5 The corresponding curves of wavelength and logarithm correlation coefficient

    Table 4 Logarithmic regression equation and correlation coefficient

    波段Band∥nm回歸方程Regressionequation相關(guān)系數(shù)CorrelationcoefficientTM1450~520Y=22.435lnX+57.160.9814TM2520~600Y=25.548lnX+44.110.9832TM3620~690Y=23.091lnX+53.860.9814TM4760~900Y=12.648lnX+87.530.9867

    注:參數(shù)X代表遙感反射率;Y代表懸浮泥沙濃度。

    Note: Parameter X stands for remote sensing reflectance; Y denotes suspended sediment concentration.

    3.2.3遙感參數(shù)模型。550 nm波段對(duì)葉綠素濃度較為敏感,而490 nm波段對(duì)葉綠素濃度也較為敏感。通過兩者的相關(guān)性分析,發(fā)現(xiàn)490 nm處的遙感反射率和550 nm處的遙感反射率高度相關(guān),兩者相除,可以去除葉綠素對(duì)懸浮泥沙遙感信息的干擾,于是韓震等[9]提出采用遙感參數(shù)XS來建立與泥沙濃度的反演模型,遙感參數(shù)XS即:

    XS=[R(λ5)+R(λ6) /R(λ4)]/R(λ5)

    (7)

    式中,R為各個(gè)波段下的遙感反射率,λ4、λ5、λ6分別對(duì)應(yīng)490、550、670 nm波段。R(λ5)+R(λ6)可以較好地反映出懸浮泥沙的特征,R(λ4)/R(λ5)可以去除葉綠素對(duì)懸浮泥沙遙感信息的干擾[9-10]。利用公式(7)計(jì)算得到6個(gè)觀測點(diǎn)的遙感參數(shù)XS(表5),與對(duì)應(yīng)的懸浮泥沙濃度做相關(guān)分析,結(jié)果發(fā)現(xiàn)當(dāng)遙感參數(shù)XS與懸浮泥沙濃度呈對(duì)數(shù)關(guān)系時(shí)相關(guān)系數(shù)最大。圖6為遙感參數(shù)XS與懸浮泥沙濃度的擬合曲線。圖6中反演模型的相關(guān)系數(shù)為0.993 4,是3個(gè)模型中最高的,可見遙感參數(shù)模型對(duì)懸浮泥沙濃度遙感反演結(jié)果最好。

    表5 不同采樣點(diǎn)的遙感參數(shù)XS

    4結(jié)論

    懸浮泥沙分布研究對(duì)舟山海域的生態(tài)建設(shè)與經(jīng)濟(jì)建設(shè)具有重要意義。該研究基于遙感技術(shù)在該海域不同海區(qū)進(jìn)行現(xiàn)場測量,建立遙感反射率與懸浮泥沙濃度之間3種不同的反演模型,結(jié)果表明:

    (1)線性模型是3個(gè)模型中相關(guān)系數(shù)最小的,在TM1~TM4波段中相關(guān)系數(shù)最大的在TM2波段,最小在TM3波段。

    圖6 泥沙遙感參數(shù)與實(shí)際含沙量的相關(guān)性分析曲線Fig.5  Correlation analysis curve between remote sensing parameters and actual sediment

    從整個(gè)波段上來看,平均相關(guān)系數(shù)最好的波段是420~600 nm和730~830 nm,兩者的平均相關(guān)系數(shù)均在0.94以上。

    (2)對(duì)數(shù)模型的相關(guān)性在3者中處于中間位置,相關(guān)系數(shù)最大的是在TM4波段,最小的是在TM3波段。

    (3)遙感參數(shù)模型是3個(gè)模型中相關(guān)系數(shù)最大的,相關(guān)系數(shù)高達(dá)0.993 4。在以后對(duì)舟山海域懸浮泥沙遙感反演研究中可以首先考慮遙感參數(shù)模型,以取得較高的精度。

    參考文獻(xiàn)

    [1] 汪小欽,陳崇成. 遙感在近岸海洋環(huán)境監(jiān)測中的應(yīng)用[J].海洋環(huán)境科學(xué),2000,19(4):72-76.

    [2] 宋亞民.舟山群島水文特性[J]. 水文,2001,21(6):59-62.

    [3] 鐘凱文,劉旭,攏解靚. 基于遙感方法反演珠江三角洲西江干流懸浮泥沙分布研究[J] 遙感應(yīng)用,2009,12(3):49-59.

    [4] 李四海,惲才興.河口表層懸浮泥沙氣象衛(wèi)星遙感定量模式研究[J].遙感學(xué)報(bào),2001,5(2):154-161.

    [5] 李炎,李京.基于海面遙感器光譜反射率料率傳遞現(xiàn)象的懸浮泥沙遙感算法[J]. 科學(xué)通報(bào),1999,44(17):1892-1897.

    [6] 朱立俊,尤玉明.遼東灣綏中海岸演變及懸沙分布特征的遙感分析[J].海洋工程,2000,18(1):66-69.

    [7] 李洪靈,張鷹,姜杰. 基于遙感方法反演懸浮泥沙分布[J ].水科學(xué)進(jìn)展,2006 ,17(2):242-2451

    [8] 鄧孺孺,何執(zhí)兼,陳曉翔. 珠江口水域水污染遙感定量分析[J ].中山大學(xué)學(xué)報(bào),2002 ,41(3):99-103

    [9] 韓震,惲才興,蔣雪中.懸浮泥沙反射光譜特性實(shí)驗(yàn)研究[J].水利學(xué)報(bào),2003(12):118-122.

    [10] 黎夏 .遙感定量的統(tǒng)一模式及其在珠江口中的應(yīng)用[J ].環(huán)境遙感,1992,7(2): 106-113.

    收稿日期2015-12-16

    作者簡介錢揚(yáng)風(fēng)(1994- ),男,安徽舒城人,本科生,專業(yè):海洋科學(xué)。*通訊作者,副教授,博士,從事海洋遙感和地理信息科學(xué)研究。

    基金項(xiàng)目浙江省大學(xué)生科技創(chuàng)新活動(dòng)計(jì)劃暨新苗人才計(jì)劃項(xiàng)目(2014R4 11036);國家自然科學(xué)基金項(xiàng)目(51479179);浙江海洋學(xué)院大學(xué)生科技創(chuàng)新項(xiàng)目;浙江海洋學(xué)院科研啟動(dòng)項(xiàng)目。

    中圖分類號(hào)S 951.4

    文獻(xiàn)標(biāo)識(shí)碼A

    文章編號(hào)0517-6611(2016)03-294-04

    av免费观看日本| 国产爱豆传媒在线观看| 中文字幕久久专区| 国产精品伦人一区二区| 国产精品久久久久久精品电影小说 | 日韩欧美在线乱码| 日本一本二区三区精品| 蜜桃久久精品国产亚洲av| 免费看av在线观看网站| 国产av不卡久久| 大型黄色视频在线免费观看| 久久久久久九九精品二区国产| 成人毛片a级毛片在线播放| 岛国毛片在线播放| 亚洲最大成人中文| 热99re8久久精品国产| 美女xxoo啪啪120秒动态图| 国产一区二区激情短视频| 亚洲欧美日韩无卡精品| 亚洲成人久久性| 全区人妻精品视频| 欧美+日韩+精品| 男女那种视频在线观看| 禁无遮挡网站| 国模一区二区三区四区视频| 国产老妇伦熟女老妇高清| 性欧美人与动物交配| 一区二区三区高清视频在线| 国产探花极品一区二区| 国产精品人妻久久久久久| 又黄又爽又刺激的免费视频.| 久99久视频精品免费| 一个人观看的视频www高清免费观看| 一级毛片久久久久久久久女| 亚洲国产欧美人成| 国产不卡一卡二| 草草在线视频免费看| 久久欧美精品欧美久久欧美| 白带黄色成豆腐渣| 精品99又大又爽又粗少妇毛片| 亚洲乱码一区二区免费版| 色5月婷婷丁香| 久久精品久久久久久久性| 国产一区亚洲一区在线观看| 日韩精品有码人妻一区| 黄片无遮挡物在线观看| 国产极品精品免费视频能看的| 日韩强制内射视频| 成人av在线播放网站| av免费观看日本| 99久久久亚洲精品蜜臀av| 在线观看免费视频日本深夜| 欧美色欧美亚洲另类二区| 国模一区二区三区四区视频| 亚洲精品国产av成人精品| 国产一区二区三区在线臀色熟女| 日本免费a在线| 亚洲av熟女| 国产一区亚洲一区在线观看| 亚洲自拍偷在线| 精品国产三级普通话版| 99国产极品粉嫩在线观看| 国产三级在线视频| 成人永久免费在线观看视频| 在线观看午夜福利视频| 成人国产麻豆网| 一个人看的www免费观看视频| 久久精品久久久久久久性| 观看美女的网站| 成人国产麻豆网| 青春草亚洲视频在线观看| 春色校园在线视频观看| 最近最新中文字幕大全电影3| 成人高潮视频无遮挡免费网站| 色哟哟哟哟哟哟| 国产蜜桃级精品一区二区三区| 少妇裸体淫交视频免费看高清| 精品人妻熟女av久视频| 国产欧美日韩精品一区二区| 亚洲精品乱码久久久久久按摩| 亚洲自偷自拍三级| 欧美色视频一区免费| 99久久人妻综合| 少妇熟女aⅴ在线视频| 神马国产精品三级电影在线观看| 狂野欧美激情性xxxx在线观看| 亚洲av免费高清在线观看| 国内久久婷婷六月综合欲色啪| 天天一区二区日本电影三级| 最新中文字幕久久久久| 少妇猛男粗大的猛烈进出视频 | 国产精品麻豆人妻色哟哟久久 | 日日摸夜夜添夜夜添av毛片| 美女内射精品一级片tv| 久久精品夜色国产| 小说图片视频综合网站| 午夜福利高清视频| 午夜精品国产一区二区电影 | 亚洲欧美成人综合另类久久久 | 三级经典国产精品| 91精品一卡2卡3卡4卡| 国产乱人偷精品视频| 99久国产av精品| 国产单亲对白刺激| 久久久久性生活片| 午夜精品在线福利| 美女内射精品一级片tv| 久久久久久久亚洲中文字幕| 日韩一本色道免费dvd| 国产一区二区在线av高清观看| 国产精品嫩草影院av在线观看| 午夜久久久久精精品| av.在线天堂| 日日撸夜夜添| 狠狠狠狠99中文字幕| 欧美3d第一页| 国产精品国产高清国产av| 国产 一区精品| 日本爱情动作片www.在线观看| 男人狂女人下面高潮的视频| .国产精品久久| www.色视频.com| 成人无遮挡网站| 在线播放国产精品三级| 日韩欧美 国产精品| 哪里可以看免费的av片| 麻豆精品久久久久久蜜桃| 97超碰精品成人国产| 亚洲真实伦在线观看| 噜噜噜噜噜久久久久久91| 久久人人爽人人片av| 欧美日本视频| 少妇丰满av| 六月丁香七月| 啦啦啦韩国在线观看视频| 一区二区三区高清视频在线| 国产成人福利小说| 亚洲人与动物交配视频| 日韩,欧美,国产一区二区三区 | 精品国产三级普通话版| 久久国内精品自在自线图片| 午夜福利视频1000在线观看| 国产午夜精品一二区理论片| 久久亚洲精品不卡| 久久久久网色| 欧美人与善性xxx| 91av网一区二区| 亚洲高清免费不卡视频| 国产精品电影一区二区三区| 国产 一区精品| 99国产精品一区二区蜜桃av| 国产精品电影一区二区三区| 欧美bdsm另类| 国产午夜精品一二区理论片| 在线观看av片永久免费下载| 狠狠狠狠99中文字幕| 国产乱人偷精品视频| 国产成人精品久久久久久| 亚洲精品粉嫩美女一区| 欧美成人精品欧美一级黄| 人妻制服诱惑在线中文字幕| 亚洲精品色激情综合| 老司机影院成人| 亚洲av中文字字幕乱码综合| 在线观看一区二区三区| 可以在线观看的亚洲视频| 国产亚洲精品av在线| 亚洲欧美精品专区久久| 一本久久精品| 亚洲欧美日韩高清在线视频| 国产真实伦视频高清在线观看| а√天堂www在线а√下载| 欧美成人免费av一区二区三区| 变态另类成人亚洲欧美熟女| 国产黄a三级三级三级人| 久久精品国产亚洲网站| 变态另类成人亚洲欧美熟女| 丰满乱子伦码专区| 人妻久久中文字幕网| 看十八女毛片水多多多| 久久久久久国产a免费观看| 久久精品久久久久久久性| 久久久久国产网址| 丰满人妻一区二区三区视频av| 亚洲va在线va天堂va国产| 黄色欧美视频在线观看| 男人舔奶头视频| 国产黄片视频在线免费观看| 日韩视频在线欧美| 免费观看a级毛片全部| 国产精品久久视频播放| 国产成人午夜福利电影在线观看| www.色视频.com| 能在线免费观看的黄片| 久久久久久大精品| 一本精品99久久精品77| 自拍偷自拍亚洲精品老妇| 色播亚洲综合网| 久久久久久久亚洲中文字幕| eeuss影院久久| 一级毛片我不卡| 欧美潮喷喷水| 热99在线观看视频| 人人妻人人看人人澡| 国产高清三级在线| 欧美+日韩+精品| 一区二区三区免费毛片| 亚洲av中文字字幕乱码综合| 校园春色视频在线观看| 18+在线观看网站| 欧美精品一区二区大全| 色吧在线观看| 国产成人a∨麻豆精品| 蜜桃亚洲精品一区二区三区| 亚洲av中文字字幕乱码综合| 久久热精品热| 九九久久精品国产亚洲av麻豆| 免费人成视频x8x8入口观看| 麻豆乱淫一区二区| 亚洲丝袜综合中文字幕| 午夜福利在线观看免费完整高清在 | 国产极品天堂在线| 变态另类丝袜制服| 亚洲国产色片| 99热这里只有是精品50| 在线观看一区二区三区| 干丝袜人妻中文字幕| 亚洲国产高清在线一区二区三| 青春草亚洲视频在线观看| 尤物成人国产欧美一区二区三区| 99国产极品粉嫩在线观看| 男插女下体视频免费在线播放| 成人午夜精彩视频在线观看| 日韩精品有码人妻一区| 国产亚洲欧美98| 91精品一卡2卡3卡4卡| 亚洲一级一片aⅴ在线观看| 色视频www国产| 夜夜看夜夜爽夜夜摸| 免费大片18禁| 欧美高清成人免费视频www| 国产蜜桃级精品一区二区三区| 日韩一区二区三区影片| 两个人的视频大全免费| 成人无遮挡网站| 熟妇人妻久久中文字幕3abv| 一级毛片电影观看 | 国产精品电影一区二区三区| 91麻豆精品激情在线观看国产| a级毛色黄片| 欧美日韩一区二区视频在线观看视频在线 | 99热6这里只有精品| 赤兔流量卡办理| 欧美精品一区二区大全| www日本黄色视频网| АⅤ资源中文在线天堂| 国产蜜桃级精品一区二区三区| 色哟哟·www| 一区福利在线观看| 在线观看一区二区三区| 我的女老师完整版在线观看| 69人妻影院| 熟女人妻精品中文字幕| 好男人在线观看高清免费视频| 国产69精品久久久久777片| 亚洲欧美精品综合久久99| 18+在线观看网站| 国产极品天堂在线| 亚洲精品粉嫩美女一区| 国产真实乱freesex| 联通29元200g的流量卡| 国产免费一级a男人的天堂| 床上黄色一级片| 一个人看的www免费观看视频| 毛片一级片免费看久久久久| 国产视频内射| 久久久久久久久久成人| 中文欧美无线码| 黑人高潮一二区| 国产精品美女特级片免费视频播放器| 国产精品野战在线观看| 在线观看av片永久免费下载| 97在线视频观看| 晚上一个人看的免费电影| 一级毛片久久久久久久久女| 美女 人体艺术 gogo| 亚洲,欧美,日韩| 精品人妻偷拍中文字幕| 国产av不卡久久| 偷拍熟女少妇极品色| 国产黄色小视频在线观看| 亚洲国产精品成人久久小说 | 麻豆成人av视频| 我的老师免费观看完整版| 看非洲黑人一级黄片| 日本免费a在线| 久久精品久久久久久久性| 少妇熟女aⅴ在线视频| 国产亚洲5aaaaa淫片| 狠狠狠狠99中文字幕| 久久精品国产亚洲av天美| 亚洲第一区二区三区不卡| 久久精品国产亚洲av香蕉五月| 毛片女人毛片| 乱码一卡2卡4卡精品| 91午夜精品亚洲一区二区三区| 中文精品一卡2卡3卡4更新| 亚洲欧美中文字幕日韩二区| 观看免费一级毛片| 亚洲一区二区三区色噜噜| 少妇的逼好多水| 中出人妻视频一区二区| 国产精品99久久久久久久久| 精华霜和精华液先用哪个| 丝袜喷水一区| 黄色一级大片看看| 中文字幕av成人在线电影| 2021天堂中文幕一二区在线观| 久久久久久伊人网av| 国产精品免费一区二区三区在线| 狂野欧美激情性xxxx在线观看| 久久精品久久久久久久性| 久久久久久久久久成人| 久久人人精品亚洲av| 亚洲aⅴ乱码一区二区在线播放| 人妻久久中文字幕网| 97超碰精品成人国产| 97人妻精品一区二区三区麻豆| 91在线精品国自产拍蜜月| 日本一本二区三区精品| 2022亚洲国产成人精品| 最近中文字幕高清免费大全6| 99热6这里只有精品| 精品久久久噜噜| 亚洲最大成人中文| 变态另类成人亚洲欧美熟女| 久久这里只有精品中国| 精品欧美国产一区二区三| 可以在线观看毛片的网站| 久久热精品热| 欧美日韩乱码在线| 亚洲天堂国产精品一区在线| 日产精品乱码卡一卡2卡三| а√天堂www在线а√下载| 精品久久久久久成人av| 黄色配什么色好看| 国产精品久久久久久亚洲av鲁大| 菩萨蛮人人尽说江南好唐韦庄 | 成人亚洲欧美一区二区av| 久久精品国产亚洲av香蕉五月| 97热精品久久久久久| 日韩大尺度精品在线看网址| 22中文网久久字幕| 日日摸夜夜添夜夜爱| 老师上课跳d突然被开到最大视频| 一进一出抽搐gif免费好疼| 一个人免费在线观看电影| 精品99又大又爽又粗少妇毛片| 中文字幕久久专区| 国产精品一区二区在线观看99 | 三级国产精品欧美在线观看| av国产免费在线观看| 嫩草影院入口| 尤物成人国产欧美一区二区三区| av福利片在线观看| 日产精品乱码卡一卡2卡三| 成人国产麻豆网| 亚洲欧美日韩卡通动漫| 男女那种视频在线观看| 亚洲av二区三区四区| 免费不卡的大黄色大毛片视频在线观看 | 女的被弄到高潮叫床怎么办| 夜夜看夜夜爽夜夜摸| 偷拍熟女少妇极品色| 成熟少妇高潮喷水视频| 亚洲电影在线观看av| 午夜亚洲福利在线播放| 亚洲精品乱码久久久v下载方式| 悠悠久久av| 69av精品久久久久久| 日本一本二区三区精品| 3wmmmm亚洲av在线观看| 亚洲人与动物交配视频| av.在线天堂| 日本av手机在线免费观看| 老女人水多毛片| 国产真实伦视频高清在线观看| 插阴视频在线观看视频| 老师上课跳d突然被开到最大视频| h日本视频在线播放| 亚洲精品色激情综合| 午夜爱爱视频在线播放| 岛国毛片在线播放| 国产精品人妻久久久影院| 日韩一区二区三区影片| 变态另类成人亚洲欧美熟女| 五月玫瑰六月丁香| 国产黄色视频一区二区在线观看 | 国语自产精品视频在线第100页| 亚洲无线在线观看| 国产乱人偷精品视频| 国产午夜精品久久久久久一区二区三区| 亚洲婷婷狠狠爱综合网| 国产成人a区在线观看| 亚洲成人久久性| 有码 亚洲区| 欧美成人免费av一区二区三区| 欧美三级亚洲精品| 免费看美女性在线毛片视频| 欧美极品一区二区三区四区| av.在线天堂| 日日摸夜夜添夜夜爱| 联通29元200g的流量卡| 久久婷婷人人爽人人干人人爱| 热99re8久久精品国产| 亚洲人成网站在线播放欧美日韩| 中文字幕熟女人妻在线| av黄色大香蕉| 波多野结衣高清无吗| 卡戴珊不雅视频在线播放| 亚洲人成网站在线播| 你懂的网址亚洲精品在线观看 | 成人国产麻豆网| 能在线免费观看的黄片| 日韩国内少妇激情av| 九九久久精品国产亚洲av麻豆| 欧美色视频一区免费| 成人亚洲欧美一区二区av| 国国产精品蜜臀av免费| 国产综合懂色| 日韩av不卡免费在线播放| 欧美xxxx黑人xx丫x性爽| 免费看光身美女| 久久久久久久午夜电影| 校园春色视频在线观看| 日韩av不卡免费在线播放| 91久久精品国产一区二区成人| 少妇裸体淫交视频免费看高清| 欧美zozozo另类| 成人欧美大片| 91麻豆精品激情在线观看国产| 十八禁国产超污无遮挡网站| 九色成人免费人妻av| 成人午夜精彩视频在线观看| 免费观看人在逋| 久久久欧美国产精品| 精品午夜福利在线看| 中文字幕制服av| 欧美精品一区二区大全| 狂野欧美激情性xxxx在线观看| 国产精品久久电影中文字幕| 悠悠久久av| 亚洲国产精品成人综合色| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品久久久久久婷婷小说 | 日韩人妻高清精品专区| 国产亚洲av嫩草精品影院| 黄色一级大片看看| 免费大片18禁| 亚洲欧美日韩高清在线视频| 亚洲欧美精品专区久久| 97热精品久久久久久| 亚洲天堂国产精品一区在线| 久久久久久九九精品二区国产| 国产精品久久视频播放| 不卡视频在线观看欧美| 精品久久久久久久久av| 国产精品久久电影中文字幕| 深夜a级毛片| 夜夜夜夜夜久久久久| 亚洲成人久久性| 一本久久精品| 亚洲在久久综合| 国产男人的电影天堂91| 免费看av在线观看网站| av福利片在线观看| 亚洲国产欧美人成| 国产乱人偷精品视频| 深夜a级毛片| 亚洲最大成人手机在线| 白带黄色成豆腐渣| 18禁在线播放成人免费| 亚洲成av人片在线播放无| 久久久久久国产a免费观看| 国产成人精品久久久久久| 国产精品女同一区二区软件| 亚洲三级黄色毛片| 1000部很黄的大片| 高清毛片免费观看视频网站| 一个人看视频在线观看www免费| 成人漫画全彩无遮挡| 亚洲av中文字字幕乱码综合| 黑人高潮一二区| 欧美色欧美亚洲另类二区| 免费观看在线日韩| 久久人人爽人人爽人人片va| 久久久久性生活片| 欧美日韩综合久久久久久| 成人漫画全彩无遮挡| 国产成人精品一,二区 | 日韩高清综合在线| 亚洲av二区三区四区| 天天躁日日操中文字幕| 少妇高潮的动态图| av在线播放精品| 在线播放无遮挡| 免费电影在线观看免费观看| 国产精品久久久久久精品电影小说 | 国产午夜精品论理片| 看非洲黑人一级黄片| 免费看av在线观看网站| 边亲边吃奶的免费视频| 亚洲最大成人手机在线| 成年女人看的毛片在线观看| 99在线视频只有这里精品首页| 91久久精品电影网| 日本欧美国产在线视频| 精品熟女少妇av免费看| 中国国产av一级| 亚洲av不卡在线观看| 麻豆久久精品国产亚洲av| 91在线精品国自产拍蜜月| 人人妻人人看人人澡| 麻豆国产97在线/欧美| 国产精品麻豆人妻色哟哟久久 | 国产精品一二三区在线看| 国产av在哪里看| 啦啦啦啦在线视频资源| 黄片无遮挡物在线观看| 国产麻豆成人av免费视频| 女同久久另类99精品国产91| 亚洲国产精品成人综合色| 欧美+日韩+精品| 欧美性猛交黑人性爽| 午夜精品一区二区三区免费看| 久久99热这里只有精品18| 午夜精品国产一区二区电影 | 中文亚洲av片在线观看爽| 在线观看午夜福利视频| 久久鲁丝午夜福利片| 久久久精品欧美日韩精品| 波多野结衣巨乳人妻| 99久久中文字幕三级久久日本| 成年av动漫网址| 成人三级黄色视频| 中出人妻视频一区二区| 欧美日韩精品成人综合77777| 老司机福利观看| 91av网一区二区| 日本黄色片子视频| 成人三级黄色视频| 大又大粗又爽又黄少妇毛片口| 日本av手机在线免费观看| av在线蜜桃| 一卡2卡三卡四卡精品乱码亚洲| 有码 亚洲区| 亚洲经典国产精华液单| 亚洲va在线va天堂va国产| 亚洲成人久久性| 午夜激情福利司机影院| 91久久精品电影网| 成人av在线播放网站| 亚洲欧美精品自产自拍| 亚洲婷婷狠狠爱综合网| 看免费成人av毛片| 免费看美女性在线毛片视频| 国产伦在线观看视频一区| 成人漫画全彩无遮挡| 日韩在线高清观看一区二区三区| 18禁在线播放成人免费| 亚洲国产精品合色在线| 卡戴珊不雅视频在线播放| 午夜福利高清视频| 干丝袜人妻中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 成年免费大片在线观看| 91午夜精品亚洲一区二区三区| 成人av在线播放网站| 国产女主播在线喷水免费视频网站 | 国产一级毛片在线| av免费观看日本| 欧美高清性xxxxhd video| 一进一出抽搐动态| 你懂的网址亚洲精品在线观看 | 国产精华一区二区三区| 亚洲第一区二区三区不卡| 人妻少妇偷人精品九色| av在线观看视频网站免费| 日本撒尿小便嘘嘘汇集6| 国产精品福利在线免费观看| 一级av片app| 欧美日韩精品成人综合77777| 成熟少妇高潮喷水视频| avwww免费| 国产av一区在线观看免费| 菩萨蛮人人尽说江南好唐韦庄 | 女人十人毛片免费观看3o分钟| 精品不卡国产一区二区三区| 国产片特级美女逼逼视频| 黄片无遮挡物在线观看| 国产久久久一区二区三区| 五月玫瑰六月丁香| 国产在线男女| 国产精品三级大全| 在线国产一区二区在线| 偷拍熟女少妇极品色| 特级一级黄色大片| 久久6这里有精品| 国产毛片a区久久久久| 嫩草影院入口| 永久网站在线| 国产精品一二三区在线看| 99久久精品热视频| 国产成人精品一,二区 | 中文字幕人妻熟人妻熟丝袜美|