• <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

    久久99精品国语久久久| 成人毛片a级毛片在线播放| 亚洲av电影在线观看一区二区三区 | 丰满人妻一区二区三区视频av| 色播亚洲综合网| 天堂网av新在线| a级一级毛片免费在线观看| 永久网站在线| or卡值多少钱| 精品久久久噜噜| 亚洲欧美成人综合另类久久久 | 久久国内精品自在自线图片| 久久精品熟女亚洲av麻豆精品 | 日韩av在线大香蕉| 国产精品99久久久久久久久| 听说在线观看完整版免费高清| 欧美潮喷喷水| 长腿黑丝高跟| 亚洲人成网站高清观看| 男女下面进入的视频免费午夜| 干丝袜人妻中文字幕| av在线观看视频网站免费| 亚洲人与动物交配视频| 乱码一卡2卡4卡精品| 日韩一本色道免费dvd| 99九九线精品视频在线观看视频| 精品人妻一区二区三区麻豆| 超碰av人人做人人爽久久| 久久久午夜欧美精品| 精品久久久久久久人妻蜜臀av| 欧美色视频一区免费| av又黄又爽大尺度在线免费看 | 欧美xxxx黑人xx丫x性爽| 91狼人影院| 国产成年人精品一区二区| 成人亚洲精品av一区二区| 国产爱豆传媒在线观看| 岛国在线免费视频观看| 可以在线观看毛片的网站| 国内揄拍国产精品人妻在线| 中文在线观看免费www的网站| 一区二区三区四区激情视频| 中文欧美无线码| 国产伦一二天堂av在线观看| 久99久视频精品免费| 97超视频在线观看视频| 中文在线观看免费www的网站| 男人舔女人下体高潮全视频| 日本黄色片子视频| 亚洲欧美精品自产自拍| 91在线精品国自产拍蜜月| 啦啦啦韩国在线观看视频| 久久久久精品久久久久真实原创| 成人性生交大片免费视频hd| 夜夜看夜夜爽夜夜摸| 午夜亚洲福利在线播放| 色尼玛亚洲综合影院| 亚洲精品日韩av片在线观看| 五月玫瑰六月丁香| 麻豆av噜噜一区二区三区| 国产高潮美女av| 中文精品一卡2卡3卡4更新| 精品99又大又爽又粗少妇毛片| 2022亚洲国产成人精品| 午夜福利视频1000在线观看| 亚洲欧美日韩东京热| 九色成人免费人妻av| 天堂影院成人在线观看| 最新中文字幕久久久久| 久久6这里有精品| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久久久久久久av| 国产日韩欧美在线精品| 午夜日本视频在线| 大香蕉97超碰在线| 老司机影院毛片| 精品免费久久久久久久清纯| 欧美变态另类bdsm刘玥| 爱豆传媒免费全集在线观看| 1024手机看黄色片| www.av在线官网国产| 中国国产av一级| 中文字幕久久专区| 亚洲在线自拍视频| 国产精品无大码| 日本免费在线观看一区| 日韩国内少妇激情av| 亚洲精品乱码久久久v下载方式| 中文乱码字字幕精品一区二区三区 | 成人国产麻豆网| 成年av动漫网址| 97人妻精品一区二区三区麻豆| 国产亚洲一区二区精品| 国产三级中文精品| 成人性生交大片免费视频hd| 亚洲精品自拍成人| 在线免费十八禁| 亚洲18禁久久av| 最后的刺客免费高清国语| 纵有疾风起免费观看全集完整版 | 成人毛片a级毛片在线播放| 国产一级毛片在线| 最近最新中文字幕免费大全7| 免费av观看视频| 欧美成人a在线观看| ponron亚洲| 国产精品熟女久久久久浪| 国产av码专区亚洲av| 亚洲aⅴ乱码一区二区在线播放| 人人妻人人看人人澡| 美女黄网站色视频| 日韩在线高清观看一区二区三区| 女人十人毛片免费观看3o分钟| 高清午夜精品一区二区三区| 岛国毛片在线播放| 高清午夜精品一区二区三区| 成年女人看的毛片在线观看| 91狼人影院| 日韩视频在线欧美| 亚洲欧美成人综合另类久久久 | 欧美高清成人免费视频www| 亚洲精品影视一区二区三区av| 精品人妻偷拍中文字幕| 尤物成人国产欧美一区二区三区| 深爱激情五月婷婷| 1000部很黄的大片| 久久精品国产亚洲av天美| 国产成人福利小说| 美女脱内裤让男人舔精品视频| 成人鲁丝片一二三区免费| 国产av不卡久久| 大香蕉久久网| 国产中年淑女户外野战色| 免费不卡的大黄色大毛片视频在线观看 | 少妇被粗大猛烈的视频| 国产激情偷乱视频一区二区| 十八禁国产超污无遮挡网站| 午夜日本视频在线| 欧美一级a爱片免费观看看| 国产成人a∨麻豆精品| 国产午夜精品久久久久久一区二区三区| 精品久久久噜噜| 一区二区三区免费毛片| 日本午夜av视频| 亚洲成av人片在线播放无| 1024手机看黄色片| 国产私拍福利视频在线观看| 国产免费男女视频| 免费黄色在线免费观看| 长腿黑丝高跟| 夜夜看夜夜爽夜夜摸| 亚洲美女搞黄在线观看| av在线天堂中文字幕| 一区二区三区高清视频在线| 99久久人妻综合| 校园人妻丝袜中文字幕| 亚洲精品亚洲一区二区| 日韩视频在线欧美| 18禁裸乳无遮挡免费网站照片| 3wmmmm亚洲av在线观看| 波多野结衣巨乳人妻| 能在线免费看毛片的网站| 国产一区二区在线观看日韩| 91狼人影院| 久久久久性生活片| 国产午夜精品论理片| 免费看美女性在线毛片视频| 欧美一区二区精品小视频在线| 老女人水多毛片| 久久午夜福利片| 热99在线观看视频| 免费看光身美女| 欧美3d第一页| 亚洲av电影在线观看一区二区三区 | 亚洲在久久综合| 韩国高清视频一区二区三区| 亚洲欧洲日产国产| 91久久精品电影网| 国产成人a区在线观看| av国产久精品久网站免费入址| 成人国产麻豆网| 国产精品1区2区在线观看.| 特大巨黑吊av在线直播| 免费大片18禁| 国产亚洲5aaaaa淫片| 精品久久久久久久人妻蜜臀av| 99在线视频只有这里精品首页| 亚洲av成人av| 国产一级毛片七仙女欲春2| 中文在线观看免费www的网站| 亚洲国产最新在线播放| 亚洲怡红院男人天堂| 亚洲成人av在线免费| 日本-黄色视频高清免费观看| 搡女人真爽免费视频火全软件| 草草在线视频免费看| 久久精品夜夜夜夜夜久久蜜豆| 97人妻精品一区二区三区麻豆| 婷婷色综合大香蕉| 热99在线观看视频| 两个人的视频大全免费| 小蜜桃在线观看免费完整版高清| 秋霞在线观看毛片| 一级二级三级毛片免费看| 欧美日韩在线观看h| 欧美人与善性xxx| 欧美成人免费av一区二区三区| 91精品伊人久久大香线蕉| 国产精品久久久久久精品电影| 日本黄色视频三级网站网址| 国产精品电影一区二区三区| 精品99又大又爽又粗少妇毛片| 国产免费又黄又爽又色| 亚洲天堂国产精品一区在线| 国产一区有黄有色的免费视频 | 亚洲国产最新在线播放| 免费在线观看成人毛片| 日韩在线高清观看一区二区三区| 日本熟妇午夜| h日本视频在线播放| 白带黄色成豆腐渣| 欧美最新免费一区二区三区| av.在线天堂| 国产精品野战在线观看| 如何舔出高潮| 嫩草影院入口| 国产精品国产三级国产av玫瑰| 日韩人妻高清精品专区| 亚洲性久久影院| 精品久久久久久久久久久久久| 99久久精品国产国产毛片| 精品久久国产蜜桃| 黄片无遮挡物在线观看| 女人久久www免费人成看片 | 中国美白少妇内射xxxbb| 精品久久久久久成人av| 午夜精品在线福利| 听说在线观看完整版免费高清| 在线观看美女被高潮喷水网站| a级毛片免费高清观看在线播放| 嘟嘟电影网在线观看| 日韩av不卡免费在线播放| a级一级毛片免费在线观看| 国产单亲对白刺激| 国产又色又爽无遮挡免| 国产黄片美女视频| 69人妻影院| 激情 狠狠 欧美| 日韩三级伦理在线观看| 日韩制服骚丝袜av| 成人三级黄色视频| 国产av在哪里看| 久久99蜜桃精品久久| 国产白丝娇喘喷水9色精品| 看片在线看免费视频| a级一级毛片免费在线观看| 人人妻人人看人人澡| 级片在线观看| 美女xxoo啪啪120秒动态图| 久久99蜜桃精品久久| 村上凉子中文字幕在线| 波多野结衣高清无吗| 99久久九九国产精品国产免费| 国产伦在线观看视频一区| 久久久久久久久中文| 国产 一区 欧美 日韩| 国产v大片淫在线免费观看| 久久精品国产亚洲网站| 国产麻豆成人av免费视频| 国产精品久久久久久精品电影小说 | av免费在线看不卡| 国产成人精品婷婷| 国产男人的电影天堂91| 国产精品一及| 亚洲精品日韩av片在线观看| 国产午夜精品久久久久久一区二区三区| 精品人妻视频免费看| 神马国产精品三级电影在线观看| 久久鲁丝午夜福利片| 国产探花极品一区二区| 天天一区二区日本电影三级| 久久久欧美国产精品| 国产激情偷乱视频一区二区| 欧美高清成人免费视频www| 精品久久久久久久人妻蜜臀av| 精品久久久久久久人妻蜜臀av| 伦精品一区二区三区| kizo精华| 91狼人影院| 久久久久免费精品人妻一区二区| 天堂av国产一区二区熟女人妻| 97超碰精品成人国产| 国产成人一区二区在线| 伦精品一区二区三区| 国产熟女欧美一区二区| 精品熟女少妇av免费看| 成年女人看的毛片在线观看| 国产免费一级a男人的天堂| 少妇高潮的动态图| 少妇高潮的动态图| 亚洲中文字幕日韩| 免费无遮挡裸体视频| 国产老妇女一区| 你懂的网址亚洲精品在线观看 | 亚洲自偷自拍三级| 男人舔女人下体高潮全视频| 七月丁香在线播放| 亚州av有码| 亚洲18禁久久av| 色噜噜av男人的天堂激情| 色噜噜av男人的天堂激情| 亚州av有码| АⅤ资源中文在线天堂| 91久久精品国产一区二区三区| 日韩强制内射视频| 丝袜美腿在线中文| 欧美成人a在线观看| 好男人视频免费观看在线| 美女cb高潮喷水在线观看| 久99久视频精品免费| 久热久热在线精品观看| 免费看av在线观看网站| 日韩视频在线欧美| 夜夜爽夜夜爽视频| 人人妻人人看人人澡| 成人综合一区亚洲| 国产精品乱码一区二三区的特点| 黄色配什么色好看| 成人欧美大片| 99久久九九国产精品国产免费| 国产在视频线在精品| 日韩视频在线欧美| 国产又色又爽无遮挡免| 国产探花极品一区二区| 亚洲成av人片在线播放无| 在线播放无遮挡| 久久久久久久久中文| 黄色配什么色好看| 亚洲欧美精品自产自拍| 久久久亚洲精品成人影院| 亚洲欧洲国产日韩| 淫秽高清视频在线观看| 亚洲国产欧美在线一区| 亚洲av电影在线观看一区二区三区 | 亚洲高清免费不卡视频| 午夜久久久久精精品| 波野结衣二区三区在线| 精品久久久久久久末码| 国产精品久久电影中文字幕| 伦理电影大哥的女人| 色噜噜av男人的天堂激情| 亚洲婷婷狠狠爱综合网| 日韩欧美国产在线观看| 麻豆久久精品国产亚洲av| 久久婷婷人人爽人人干人人爱| 亚洲欧美成人精品一区二区| 亚洲国产最新在线播放| 久久久久久久国产电影| 少妇高潮的动态图| 国产真实伦视频高清在线观看| av在线天堂中文字幕| 美女高潮的动态| 天天一区二区日本电影三级| 欧美bdsm另类| 乱人视频在线观看| 国产69精品久久久久777片| 波多野结衣巨乳人妻| 国产高潮美女av| 亚洲精品乱久久久久久| av国产久精品久网站免费入址| 国产精品美女特级片免费视频播放器| 久久久久久久国产电影| 日韩成人伦理影院| 午夜日本视频在线| 国产黄a三级三级三级人| 亚洲国产色片| 99热这里只有精品一区| 亚洲在久久综合| 亚洲av.av天堂| 亚洲自偷自拍三级| 免费黄色在线免费观看| 日韩成人伦理影院| 久久这里只有精品中国| 久久久国产成人免费| 少妇的逼水好多| 2021天堂中文幕一二区在线观| 久久韩国三级中文字幕| 春色校园在线视频观看| 成人午夜精彩视频在线观看| 五月伊人婷婷丁香| 亚洲aⅴ乱码一区二区在线播放| 中文字幕免费在线视频6| 小说图片视频综合网站| 亚洲国产精品专区欧美| 女人十人毛片免费观看3o分钟| 亚洲精品456在线播放app| 国产精品国产高清国产av| 久久精品久久久久久噜噜老黄 | 两个人视频免费观看高清| 国产真实乱freesex| 中文字幕久久专区| 精品人妻偷拍中文字幕| 欧美性猛交黑人性爽| 久久精品国产亚洲网站| 高清视频免费观看一区二区 | 国产探花在线观看一区二区| 国产女主播在线喷水免费视频网站 | 久久久a久久爽久久v久久| 国产高潮美女av| 国产精品日韩av在线免费观看| 精品人妻熟女av久视频| 97超碰精品成人国产| 免费看a级黄色片| 成人午夜高清在线视频| 国产av一区在线观看免费| 婷婷色麻豆天堂久久 | 成人午夜精彩视频在线观看| 午夜a级毛片| 午夜激情欧美在线| 欧美+日韩+精品| 九九久久精品国产亚洲av麻豆| 国产片特级美女逼逼视频| 不卡视频在线观看欧美| 精品久久久久久久久av| 91精品国产九色| 久久99热这里只有精品18| 国产精品精品国产色婷婷| 99久国产av精品| 免费av不卡在线播放| 免费av毛片视频| 男人和女人高潮做爰伦理| 国产色婷婷99| 亚洲精品自拍成人| 视频中文字幕在线观看| 国产91av在线免费观看| 黄色配什么色好看| 精品久久久久久电影网 | 欧美另类亚洲清纯唯美| 国产在线男女| 亚洲自拍偷在线| 性插视频无遮挡在线免费观看| 九九爱精品视频在线观看| 精品酒店卫生间| 搡女人真爽免费视频火全软件| 国产精品久久久久久av不卡| 国产大屁股一区二区在线视频| 成人毛片60女人毛片免费| 久久久精品欧美日韩精品| 久久综合国产亚洲精品| 观看美女的网站| 亚洲久久久久久中文字幕| 亚洲第一区二区三区不卡| 亚洲av福利一区| 国产成人aa在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 色综合亚洲欧美另类图片| 男人的好看免费观看在线视频| 亚洲av中文字字幕乱码综合| 国产精品福利在线免费观看| 美女内射精品一级片tv| 黄色欧美视频在线观看| 欧美又色又爽又黄视频| 天堂√8在线中文| 午夜激情福利司机影院| 日日摸夜夜添夜夜添av毛片| 久久国产乱子免费精品| 午夜免费激情av| 天堂av国产一区二区熟女人妻| 久久久午夜欧美精品| 久久韩国三级中文字幕| 国产色婷婷99| av播播在线观看一区| 欧美成人一区二区免费高清观看| 可以在线观看毛片的网站| 午夜爱爱视频在线播放| 免费观看精品视频网站| 黄片无遮挡物在线观看| 国语对白做爰xxxⅹ性视频网站| 国产精品电影一区二区三区| 汤姆久久久久久久影院中文字幕 | 亚洲在线自拍视频| 成年女人看的毛片在线观看| 男女下面进入的视频免费午夜| 春色校园在线视频观看| 麻豆一二三区av精品| 亚洲欧美成人综合另类久久久 | 特大巨黑吊av在线直播| .国产精品久久| 久久久精品94久久精品| 亚洲国产精品久久男人天堂| 精品人妻偷拍中文字幕| 69人妻影院| 啦啦啦观看免费观看视频高清| 国产午夜福利久久久久久| 久久久久久久午夜电影| 国产v大片淫在线免费观看| 你懂的网址亚洲精品在线观看 | 国产在视频线在精品| 看十八女毛片水多多多| 国产一区亚洲一区在线观看| 青春草国产在线视频| 欧美最新免费一区二区三区| 日本猛色少妇xxxxx猛交久久| 九九爱精品视频在线观看| 欧美区成人在线视频| 亚洲欧洲日产国产| 99久国产av精品| 久久精品91蜜桃| 深爱激情五月婷婷| 尤物成人国产欧美一区二区三区| 亚洲av中文av极速乱| 亚洲五月天丁香| 啦啦啦韩国在线观看视频| 欧美精品国产亚洲| 国产人妻一区二区三区在| 国产探花极品一区二区| 黄片无遮挡物在线观看| 久久国内精品自在自线图片| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲欧洲国产日韩| 男人和女人高潮做爰伦理| 少妇熟女aⅴ在线视频| 日韩大片免费观看网站 | 欧美成人一区二区免费高清观看| 国产久久久一区二区三区| 国产av在哪里看| 国内精品宾馆在线| 国产精品国产高清国产av| 亚洲国产欧洲综合997久久,| 国产v大片淫在线免费观看| 少妇高潮的动态图| 久久精品综合一区二区三区| 内地一区二区视频在线| 国产爱豆传媒在线观看| 亚洲五月天丁香| 91精品一卡2卡3卡4卡| 日韩精品有码人妻一区| 日韩高清综合在线| 18+在线观看网站| 在线免费观看的www视频| 国产精品久久久久久久久免| 中文字幕av成人在线电影| 日韩欧美精品免费久久| 国产精品99久久久久久久久| 少妇的逼好多水| 尤物成人国产欧美一区二区三区| 蜜臀久久99精品久久宅男| 久久久色成人| 你懂的网址亚洲精品在线观看 | 亚洲天堂国产精品一区在线| 看黄色毛片网站| 一边亲一边摸免费视频| 丝袜喷水一区| 美女黄网站色视频| 日韩欧美国产在线观看| 国内揄拍国产精品人妻在线| 一级av片app| 少妇人妻一区二区三区视频| 一级av片app| 免费观看a级毛片全部| 中文字幕亚洲精品专区| 日本午夜av视频| 久久草成人影院| 精品久久久久久久久亚洲| 久久草成人影院| www日本黄色视频网| 精品99又大又爽又粗少妇毛片| 亚洲精品成人久久久久久| 18禁动态无遮挡网站| 男人舔女人下体高潮全视频| 免费黄网站久久成人精品| 国产一区二区在线观看日韩| 国产久久久一区二区三区| 黑人高潮一二区| 国产成人a区在线观看| 男女边吃奶边做爰视频| 国产午夜精品一二区理论片| 免费看美女性在线毛片视频| 国产伦理片在线播放av一区| av视频在线观看入口| 国产高清不卡午夜福利| 一边摸一边抽搐一进一小说| .国产精品久久| 亚洲精品国产av成人精品| av.在线天堂| 国产精品熟女久久久久浪| 大话2 男鬼变身卡| 久久久久网色| 在线观看av片永久免费下载| a级毛色黄片| av在线观看视频网站免费| 极品教师在线视频| 日韩制服骚丝袜av| 中文字幕久久专区| 狂野欧美激情性xxxx在线观看| 欧美+日韩+精品| av在线播放精品| 少妇高潮的动态图| 精品久久久久久久末码| h日本视频在线播放| 亚洲国产最新在线播放| 亚洲av日韩在线播放| 中文字幕精品亚洲无线码一区| 久久99蜜桃精品久久| 三级国产精品片| 自拍偷自拍亚洲精品老妇| 在线免费观看不下载黄p国产| 男女啪啪激烈高潮av片| 最近2019中文字幕mv第一页| 蜜臀久久99精品久久宅男| 女人被狂操c到高潮| 色噜噜av男人的天堂激情| 男女那种视频在线观看| 又粗又硬又长又爽又黄的视频| 午夜免费男女啪啪视频观看|