• <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在线大香蕉| 男人狂女人下面高潮的视频| 熟妇人妻久久中文字幕3abv| 亚洲av中文av极速乱| 国产精品一及| 亚洲精品影视一区二区三区av| 久久人人爽人人片av| 黄色欧美视频在线观看| 1000部很黄的大片| 深夜精品福利| 久久九九热精品免费| 内地一区二区视频在线| 国产黄片美女视频| 亚洲欧洲国产日韩| 日韩av在线大香蕉| av黄色大香蕉| 免费av不卡在线播放| 91av网一区二区| 激情 狠狠 欧美| 熟女人妻精品中文字幕| 小蜜桃在线观看免费完整版高清| 日韩精品青青久久久久久| 久久精品影院6| 一本久久精品| 亚洲在线自拍视频| 一进一出抽搐gif免费好疼| 色吧在线观看| 简卡轻食公司| 国产午夜精品一二区理论片| av天堂在线播放| 国产精品日韩av在线免费观看| 精品欧美国产一区二区三| av黄色大香蕉| 亚洲成人久久性| 亚洲自偷自拍三级| 国产精品一区www在线观看| 亚洲欧美日韩高清在线视频| 精品人妻偷拍中文字幕| 青春草亚洲视频在线观看| 亚洲成人久久爱视频| 国产在线精品亚洲第一网站| 边亲边吃奶的免费视频| 国产精品爽爽va在线观看网站| 麻豆国产av国片精品| 久久久久久久亚洲中文字幕| 永久网站在线| 久久中文看片网| 日韩欧美精品免费久久| 黄色欧美视频在线观看| 欧美又色又爽又黄视频| 校园人妻丝袜中文字幕| 日本三级黄在线观看| 看十八女毛片水多多多| 爱豆传媒免费全集在线观看| 在线免费十八禁| 联通29元200g的流量卡| 亚洲真实伦在线观看| 亚洲最大成人av| 国产黄色视频一区二区在线观看 | 两个人的视频大全免费| 噜噜噜噜噜久久久久久91| 少妇熟女欧美另类| 国产探花极品一区二区| 亚洲成人精品中文字幕电影| 国产精品精品国产色婷婷| 亚洲人与动物交配视频| 91狼人影院| 国产亚洲5aaaaa淫片| 狂野欧美激情性xxxx在线观看| 亚洲人与动物交配视频| 亚洲最大成人中文| 少妇丰满av| 毛片一级片免费看久久久久| 十八禁国产超污无遮挡网站| 日韩人妻高清精品专区| 嫩草影院精品99| 久久99热6这里只有精品| 日韩在线高清观看一区二区三区| 看非洲黑人一级黄片| 91久久精品国产一区二区成人| 人妻久久中文字幕网| 两性午夜刺激爽爽歪歪视频在线观看| 国产伦精品一区二区三区视频9| 亚洲精品影视一区二区三区av| 尤物成人国产欧美一区二区三区| av在线蜜桃| 国产av不卡久久| 国内精品宾馆在线| 久久亚洲国产成人精品v| 日本与韩国留学比较| 一个人免费在线观看电影| 国产一级毛片在线| 高清午夜精品一区二区三区 | 男女那种视频在线观看| 国产精品乱码一区二三区的特点| 亚洲成人久久性| 国产成人91sexporn| 美女xxoo啪啪120秒动态图| 麻豆成人午夜福利视频| 在线播放无遮挡| 亚洲美女视频黄频| 日日干狠狠操夜夜爽| 国产精品蜜桃在线观看 | 国产一区二区三区在线臀色熟女| 午夜福利视频1000在线观看| 特大巨黑吊av在线直播| 国产精品电影一区二区三区| 国产精品久久久久久久电影| 亚洲熟妇中文字幕五十中出| 干丝袜人妻中文字幕| 天天一区二区日本电影三级| 日日撸夜夜添| 简卡轻食公司| 欧美激情久久久久久爽电影| 黑人高潮一二区| 性插视频无遮挡在线免费观看| 干丝袜人妻中文字幕| 岛国毛片在线播放| 日韩精品青青久久久久久| 久久久精品欧美日韩精品| 国产v大片淫在线免费观看| 婷婷精品国产亚洲av| 国产私拍福利视频在线观看| 国产中年淑女户外野战色| 国产午夜精品论理片| 国产一区二区在线观看日韩| 1024手机看黄色片| 一个人看视频在线观看www免费| 哪个播放器可以免费观看大片| 欧美成人a在线观看| 亚洲人成网站高清观看| 欧美日韩在线观看h| 免费观看人在逋| 中文字幕免费在线视频6| 国产熟女欧美一区二区| or卡值多少钱| 中文字幕人妻熟人妻熟丝袜美| 两个人的视频大全免费| 欧美性猛交╳xxx乱大交人| 亚洲av不卡在线观看| 免费av不卡在线播放| 中文精品一卡2卡3卡4更新| 欧美xxxx黑人xx丫x性爽| 亚洲人成网站在线播| av.在线天堂| 床上黄色一级片| 亚洲一区高清亚洲精品| 国产极品精品免费视频能看的| 一级毛片电影观看 | 国产v大片淫在线免费观看| 三级男女做爰猛烈吃奶摸视频| 欧美精品国产亚洲| 亚洲精品456在线播放app| 免费无遮挡裸体视频| av又黄又爽大尺度在线免费看 | 成人亚洲精品av一区二区| 在线观看66精品国产| 一边亲一边摸免费视频| 男人狂女人下面高潮的视频| 1000部很黄的大片| 色播亚洲综合网| av在线老鸭窝| www日本黄色视频网| 禁无遮挡网站| 人妻制服诱惑在线中文字幕| 哪里可以看免费的av片| 蜜桃久久精品国产亚洲av| 国产精品久久久久久久电影| 日本av手机在线免费观看| 99久久成人亚洲精品观看| 国产成人精品久久久久久| 日韩强制内射视频| 国产在视频线在精品| 日韩欧美 国产精品| 国产精品精品国产色婷婷| 99热这里只有是精品在线观看| 伦理电影大哥的女人| 深夜精品福利| 成人亚洲精品av一区二区| 天堂中文最新版在线下载 | 免费一级毛片在线播放高清视频| 成人国产麻豆网| 麻豆国产av国片精品| 国产精品99久久久久久久久| 变态另类丝袜制服| 99久久人妻综合| 成人毛片a级毛片在线播放| 国产黄片视频在线免费观看| 国产欧美日韩精品一区二区| 搡女人真爽免费视频火全软件| 国产真实乱freesex| 亚洲欧美日韩东京热| 99久久中文字幕三级久久日本| 51国产日韩欧美| 网址你懂的国产日韩在线| 天堂√8在线中文| 干丝袜人妻中文字幕| 你懂的网址亚洲精品在线观看 | 免费看av在线观看网站| 99久久久亚洲精品蜜臀av| 国产黄片视频在线免费观看| 久久久精品大字幕| 午夜a级毛片| 在现免费观看毛片| 久久亚洲国产成人精品v| 久久久久久伊人网av| 麻豆乱淫一区二区| 亚洲最大成人中文| 搡女人真爽免费视频火全软件| 能在线免费观看的黄片| 国产白丝娇喘喷水9色精品| 床上黄色一级片| 精品久久久久久成人av| 国产极品天堂在线| av女优亚洲男人天堂| 免费看光身美女| 黑人高潮一二区| 男人狂女人下面高潮的视频| 亚洲精品亚洲一区二区| 久久久色成人| 看十八女毛片水多多多| 亚洲av一区综合| 天堂网av新在线| 村上凉子中文字幕在线| 成人av在线播放网站| 少妇高潮的动态图| 天堂√8在线中文| 麻豆乱淫一区二区| 99久久久亚洲精品蜜臀av| 亚洲国产精品成人综合色| 免费观看精品视频网站| 免费电影在线观看免费观看| 国产91av在线免费观看| 日韩欧美精品v在线| 男插女下体视频免费在线播放| 亚洲国产日韩欧美精品在线观看| 五月玫瑰六月丁香| 色尼玛亚洲综合影院| 国产亚洲5aaaaa淫片| 美女cb高潮喷水在线观看| 欧美激情国产日韩精品一区| or卡值多少钱| av专区在线播放| 国内揄拍国产精品人妻在线| 日本黄色片子视频| 亚洲精品日韩在线中文字幕 | 看免费成人av毛片| 少妇人妻精品综合一区二区 | 亚洲av第一区精品v没综合| 欧美最黄视频在线播放免费| 男女做爰动态图高潮gif福利片| 狂野欧美白嫩少妇大欣赏| 国产亚洲精品久久久久久毛片| 欧美性感艳星| 18+在线观看网站| 国产精品免费一区二区三区在线| 亚洲图色成人| 看非洲黑人一级黄片| 非洲黑人性xxxx精品又粗又长| 日本黄色视频三级网站网址| av卡一久久| 91狼人影院| 变态另类丝袜制服| 国产色婷婷99| 国产成年人精品一区二区| 极品教师在线视频| 国产一区二区三区av在线 | 欧美bdsm另类| 久久精品国产99精品国产亚洲性色| 男的添女的下面高潮视频| 男女视频在线观看网站免费| 91精品一卡2卡3卡4卡| 国产亚洲91精品色在线| 亚洲av成人精品一区久久| 最新中文字幕久久久久| 一级av片app| 欧美+亚洲+日韩+国产| 99riav亚洲国产免费| 国产成人freesex在线| 91久久精品国产一区二区三区| 亚洲天堂国产精品一区在线| 久久99热这里只有精品18| 一级av片app| 欧美+亚洲+日韩+国产| 99久国产av精品国产电影| 国产精品1区2区在线观看.| 午夜亚洲福利在线播放| 尾随美女入室| 青春草亚洲视频在线观看| 99热精品在线国产| 爱豆传媒免费全集在线观看| 国产又黄又爽又无遮挡在线| 久久精品国产清高在天天线| 别揉我奶头 嗯啊视频| 村上凉子中文字幕在线| 久久精品久久久久久久性| 日本免费一区二区三区高清不卡| 日本一二三区视频观看| 欧美日韩一区二区视频在线观看视频在线 | 舔av片在线| 联通29元200g的流量卡| 国产av不卡久久| 尤物成人国产欧美一区二区三区| 国产蜜桃级精品一区二区三区| 国产69精品久久久久777片| 日韩一本色道免费dvd| 久久九九热精品免费| 成人亚洲欧美一区二区av| 亚洲欧美日韩卡通动漫| 久久精品国产亚洲av涩爱 | 黄色日韩在线| 色播亚洲综合网| 不卡视频在线观看欧美| 我要搜黄色片| 一个人看视频在线观看www免费| 日韩强制内射视频| 赤兔流量卡办理| 亚洲18禁久久av| 老师上课跳d突然被开到最大视频| 少妇丰满av| 久久草成人影院| 国产av在哪里看| 午夜激情欧美在线| 日产精品乱码卡一卡2卡三| 国产片特级美女逼逼视频| 女同久久另类99精品国产91| 美女内射精品一级片tv| 在线a可以看的网站| 少妇猛男粗大的猛烈进出视频 | 亚洲欧美日韩高清专用| 人妻夜夜爽99麻豆av| 插阴视频在线观看视频| av.在线天堂| 日韩高清综合在线| 高清午夜精品一区二区三区 | 此物有八面人人有两片| 免费在线观看成人毛片| 在线观看av片永久免费下载| 非洲黑人性xxxx精品又粗又长| 日本免费一区二区三区高清不卡| 最近的中文字幕免费完整| 悠悠久久av| 久久国内精品自在自线图片| 国内精品美女久久久久久| 毛片女人毛片| 12—13女人毛片做爰片一| 1024手机看黄色片| 久久热精品热| 欧美xxxx性猛交bbbb| 国内少妇人妻偷人精品xxx网站| 日本欧美国产在线视频| 不卡一级毛片| 可以在线观看的亚洲视频| 身体一侧抽搐| 国产精品美女特级片免费视频播放器| 日本免费一区二区三区高清不卡| 中文字幕精品亚洲无线码一区| 欧美三级亚洲精品| 久久久国产成人免费| 国产成人精品久久久久久| 亚洲最大成人手机在线| 身体一侧抽搐| 亚洲人与动物交配视频| 国产亚洲5aaaaa淫片| 人妻系列 视频| 亚洲精品乱码久久久久久按摩| 国产真实乱freesex| 91久久精品电影网| 亚洲人与动物交配视频| 免费看a级黄色片| 夜夜夜夜夜久久久久| 欧美又色又爽又黄视频| 性欧美人与动物交配| 亚洲欧美日韩无卡精品| 亚洲av不卡在线观看| 色噜噜av男人的天堂激情| 中文字幕免费在线视频6| 色噜噜av男人的天堂激情| 精品99又大又爽又粗少妇毛片| 久久精品国产自在天天线| 成人亚洲精品av一区二区| 在线观看一区二区三区| 亚洲精品乱码久久久v下载方式| 久久亚洲精品不卡| 久久人妻av系列| 五月玫瑰六月丁香| 午夜视频国产福利| 我要看日韩黄色一级片| 久久久久久久久久成人| 69av精品久久久久久| 亚洲经典国产精华液单| 禁无遮挡网站| 国内少妇人妻偷人精品xxx网站| 国产麻豆成人av免费视频| 亚洲精品久久久久久婷婷小说 | 又爽又黄a免费视频| 99久久中文字幕三级久久日本| 真实男女啪啪啪动态图| 麻豆国产97在线/欧美| 国产午夜精品一二区理论片| 免费看美女性在线毛片视频| 婷婷亚洲欧美| 国产成人影院久久av| 国产一级毛片七仙女欲春2| 国产老妇伦熟女老妇高清| 一区二区三区高清视频在线| 欧美色欧美亚洲另类二区| 国产一区二区三区av在线 | 国产成人精品婷婷| 18+在线观看网站| 青青草视频在线视频观看| 一个人看视频在线观看www免费| 成人午夜高清在线视频| 黄色日韩在线| 日韩亚洲欧美综合| 亚洲精品久久久久久婷婷小说 | 丝袜喷水一区| 国产av不卡久久| 国产在线男女| 超碰av人人做人人爽久久| 看免费成人av毛片| 一本久久精品| 日韩人妻高清精品专区| 一级av片app| 精品国产三级普通话版| 精品久久久久久久久久免费视频| 欧美性猛交黑人性爽| 国产午夜精品论理片| 日日摸夜夜添夜夜爱| 欧美日韩综合久久久久久| 中文字幕熟女人妻在线| 看免费成人av毛片| 一本久久精品| 91午夜精品亚洲一区二区三区| 成年免费大片在线观看| 日韩中字成人| 最近最新中文字幕大全电影3| 国产一区亚洲一区在线观看| 免费一级毛片在线播放高清视频| 观看美女的网站| 亚洲经典国产精华液单| 1024手机看黄色片| 国产精品福利在线免费观看| 欧美高清性xxxxhd video| 免费av不卡在线播放| 国产精品久久视频播放| 色哟哟哟哟哟哟| 秋霞在线观看毛片| 久久精品夜夜夜夜夜久久蜜豆| 久久精品国产鲁丝片午夜精品| 亚洲国产欧美在线一区| 国产亚洲精品av在线| 丝袜美腿在线中文| 身体一侧抽搐| 啦啦啦啦在线视频资源| 久久久成人免费电影| 国内精品久久久久精免费| 狂野欧美激情性xxxx在线观看| 日韩,欧美,国产一区二区三区 | 最近最新中文字幕大全电影3| 有码 亚洲区| 久久亚洲国产成人精品v| 亚洲aⅴ乱码一区二区在线播放| 中文亚洲av片在线观看爽| 熟女人妻精品中文字幕| 亚洲精品成人久久久久久| 亚洲国产精品国产精品| 麻豆国产97在线/欧美| 国产伦一二天堂av在线观看| 18禁在线无遮挡免费观看视频| 亚洲国产高清在线一区二区三| 成人性生交大片免费视频hd| 精品人妻视频免费看| 亚洲国产精品国产精品| 国产精品日韩av在线免费观看| 久久精品国产亚洲av涩爱 | 波多野结衣高清无吗| 亚洲欧美精品综合久久99| 国产精品免费一区二区三区在线| 免费av观看视频| av在线播放精品| 在线观看免费视频日本深夜| 成人av在线播放网站| 国产真实伦视频高清在线观看| 国产av一区在线观看免费| 精品久久国产蜜桃| 欧美激情久久久久久爽电影| 1024手机看黄色片| 久久国产乱子免费精品| 99热这里只有是精品50| av在线天堂中文字幕| av福利片在线观看| 热99re8久久精品国产| 26uuu在线亚洲综合色| 天堂网av新在线| 超碰av人人做人人爽久久| 小说图片视频综合网站| 欧美xxxx性猛交bbbb| 少妇丰满av| 日韩大尺度精品在线看网址| 色综合站精品国产| 人人妻人人看人人澡| 亚洲最大成人av| 久久久久久久久久久免费av| 久久国内精品自在自线图片| 精品一区二区三区人妻视频| 亚洲欧美日韩高清专用| 国产真实伦视频高清在线观看| 中文在线观看免费www的网站| 国产视频首页在线观看| 三级男女做爰猛烈吃奶摸视频| 最近2019中文字幕mv第一页| 五月伊人婷婷丁香| 日韩大尺度精品在线看网址| 69人妻影院| 永久网站在线| 中出人妻视频一区二区| 中文精品一卡2卡3卡4更新| 精品少妇黑人巨大在线播放 | 联通29元200g的流量卡| av在线播放精品| 麻豆久久精品国产亚洲av| 91精品国产九色| 亚洲国产欧美在线一区| 亚洲最大成人中文| 在线免费十八禁| 久久精品国产自在天天线| 99精品在免费线老司机午夜| 欧美激情在线99| 欧美一区二区精品小视频在线| 中出人妻视频一区二区| 变态另类丝袜制服| 久久久久久久久久久免费av| 少妇人妻精品综合一区二区 | 亚洲欧美精品综合久久99| 欧美精品国产亚洲| 欧美激情久久久久久爽电影| 午夜爱爱视频在线播放| 亚洲欧洲国产日韩| 午夜免费男女啪啪视频观看| 亚洲精品久久国产高清桃花| 国产亚洲5aaaaa淫片| 国产 一区 欧美 日韩| 激情 狠狠 欧美| 国产熟女欧美一区二区| 性欧美人与动物交配| 亚洲国产欧洲综合997久久,| 国产精品国产高清国产av| 可以在线观看的亚洲视频| 国产精品av视频在线免费观看| 看黄色毛片网站| 欧美xxxx黑人xx丫x性爽| 国产av在哪里看| 日本一本二区三区精品| 91在线精品国自产拍蜜月| 日韩欧美 国产精品| 欧美高清性xxxxhd video| 91久久精品国产一区二区三区| 欧美日韩综合久久久久久| 日本黄色片子视频| 亚洲av成人av| 伦理电影大哥的女人| 午夜激情福利司机影院| 久久中文看片网| 伊人久久精品亚洲午夜| 能在线免费观看的黄片| 色综合亚洲欧美另类图片| 欧美高清成人免费视频www| 最近手机中文字幕大全| 18禁黄网站禁片免费观看直播| 亚洲丝袜综合中文字幕| 精品久久久久久久末码| 99国产极品粉嫩在线观看| 又黄又爽又刺激的免费视频.| 欧美日韩乱码在线| 国产精品综合久久久久久久免费| 亚洲久久久久久中文字幕| 国产一区二区三区在线臀色熟女| 亚洲电影在线观看av| 如何舔出高潮| 一级毛片电影观看 | 免费看a级黄色片| 国产日韩欧美在线精品| 国产高清有码在线观看视频| 国产女主播在线喷水免费视频网站 | av在线播放精品| 日韩av不卡免费在线播放| 亚洲天堂国产精品一区在线| 日本爱情动作片www.在线观看| 22中文网久久字幕| 亚洲国产精品国产精品| 午夜亚洲福利在线播放| ponron亚洲| 老师上课跳d突然被开到最大视频| 嫩草影院新地址| 男女啪啪激烈高潮av片| 级片在线观看| 国产高清有码在线观看视频| 男人舔奶头视频| 国产在线精品亚洲第一网站| 卡戴珊不雅视频在线播放| 亚洲av中文av极速乱| 亚洲在线观看片| 丝袜美腿在线中文| 人人妻人人澡欧美一区二区| 啦啦啦韩国在线观看视频| 午夜老司机福利剧场| 91av网一区二区| 久久久欧美国产精品| 一本一本综合久久| 久久久成人免费电影| 天天一区二区日本电影三级|