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

    陜西連續(xù)GPS基準(zhǔn)站坐標(biāo)時(shí)間序列分析*

    2014-09-20 08:05:54蘇利娜丁曉光張彥芬張永奇
    關(guān)鍵詞:共模漫步基準(zhǔn)

    蘇利娜 丁曉光 張彥芬 張永奇

    1)中國(guó)地震局地質(zhì)研究所,北京 100029

    2)陜西省地震局,西安 710068

    3)山西省測(cè)繪工程院,太原 030002

    陜西連續(xù)GPS基準(zhǔn)站坐標(biāo)時(shí)間序列分析*

    蘇利娜1,2)丁曉光1,2)張彥芬3)張永奇1)

    1)中國(guó)地震局地質(zhì)研究所,北京 100029

    2)陜西省地震局,西安 710068

    3)山西省測(cè)繪工程院,太原 030002

    對(duì)陜西省GPS連續(xù)基準(zhǔn)站2011-01-01~2013-05-30坐標(biāo)時(shí)間序列進(jìn)行分析,利用區(qū)域疊加濾波分離出共模誤差,濾波后站坐標(biāo)時(shí)間序列N、E、U方向的平均均方根分別減小69%、60%和41%。通過(guò)譜指數(shù)計(jì)算發(fā)現(xiàn),空間濾波后的時(shí)間序列含有白噪聲和有色噪聲;利用CATS軟件對(duì)各噪聲分量估計(jì)發(fā)現(xiàn),垂直方向噪聲分量大于水平方向;空間濾波使閃爍噪聲大幅減小,部分站出現(xiàn)隨機(jī)漫步噪聲。

    陜西省連續(xù)GPS基準(zhǔn)站;時(shí)間序列;區(qū)域疊加濾波;共模誤差;噪聲特征

    利用GNSS觀測(cè)技術(shù)研究地殼運(yùn)動(dòng)和構(gòu)造變形,是地球動(dòng)力學(xué)研究領(lǐng)域的熱點(diǎn)之一。然而,GPS位置時(shí)間序列中還包含了時(shí)、空相關(guān)的噪聲,與構(gòu)造變形引起的位移混疊在一起。時(shí)間序列的研究有助于去除GPS數(shù)據(jù)中的噪聲影響和非構(gòu)造形變,得到對(duì)研究地震有用的構(gòu)造形變信息[1-2]。

    GPS坐標(biāo)時(shí)間序列中的噪聲研究主要分為時(shí)域分析和空域分析。時(shí)域分析是確定坐標(biāo)時(shí)間序列在時(shí)間上的相關(guān)性特征,其方法有譜分析和最大似然估計(jì)(MLE,maximum likelihood estimator),用以確定時(shí)間相關(guān)噪聲的類型和強(qiáng)度。大量研究表明,坐標(biāo)時(shí)間序列中包含有白噪聲和與時(shí)間相關(guān)的有色噪聲[3-10]。空域分析是研究區(qū)域GPS網(wǎng)站點(diǎn)坐標(biāo)時(shí)間序列之間的空間相關(guān)性,通過(guò)剔除共模誤差(CME,common mode error)可以有效提高站坐標(biāo)精度[5,7],常用方法有區(qū)域疊加濾波、主成分分析和 KL變換空間濾波。本文首先去除坐標(biāo)時(shí)間序列中的地殼構(gòu)造形變信息,利用區(qū)域疊加濾波分離出坐標(biāo)殘差時(shí)間序列中的共模誤差,比較空間疊加濾波前后基準(zhǔn)站坐標(biāo)平均均方根。然后,利用 CATS軟件[9]分析空間濾波后的基準(zhǔn)站坐標(biāo)時(shí)間序列的頻譜特征和噪聲特性。

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

    首先采用GAMIT軟件對(duì)陜西連續(xù)GPS基準(zhǔn)站2011-01-01~2013-05-30的數(shù)據(jù)進(jìn)行基線解算。解算策略為BASELINE,觀測(cè)值類型為L(zhǎng)C_AUTCLN,高度角為10°,采樣率為30 s,每隔2 h對(duì)對(duì)流層進(jìn)行估計(jì),并估計(jì)大氣梯度等。然后,通過(guò)約束9個(gè)IGS參考站在ITRF08框架下的坐標(biāo)和速度,利用GLOBK軟件平差得到ITRF08框架下的陜西GPS連續(xù)觀測(cè)站的坐標(biāo)時(shí)間序列。各個(gè)坐標(biāo)時(shí)間序列水平分量的重復(fù)性(加權(quán)均方根誤差)優(yōu)于4 mm,垂直分量的重復(fù)性優(yōu)于6 mm。

    2 共模誤差

    Wdowinski等[11]研究表明,不同 GPS基準(zhǔn)站之間存在共模誤差,但誤差來(lái)源尚未明確[2]。為了獲得可靠的GPS數(shù)據(jù),需要提高GPS觀測(cè)結(jié)果的信噪比,并通過(guò)一些方法降低噪聲的影響。對(duì)于區(qū)域GPS坐標(biāo)時(shí)間序列,區(qū)域疊加濾波、主成分分析、K-L濾波等都能很好地去除共模誤差,提高信噪比[12]。

    區(qū)域疊加濾波是將每個(gè)歷元各站殘差的加權(quán)平均值當(dāng)作該歷元的共模誤差,算法簡(jiǎn)單易實(shí)現(xiàn)。本文采用區(qū)域疊加濾波法來(lái)分離共模誤差。首先利用擬合函數(shù)去除趨勢(shì)項(xiàng)和周期項(xiàng),獲得殘差坐標(biāo)時(shí)間序列;然后以GPS單日解中誤差定權(quán),利用區(qū)域疊加濾波提取出共模誤差(圖1);最后在原始時(shí)間序列中剔除共模誤差,時(shí)間序列中的周期性和線性趨勢(shì)更加明顯(圖2),站點(diǎn)精度也有所提高(表1)。

    從濾波前后23個(gè)站點(diǎn)的精度可以看出,利用區(qū)域疊加濾波消除共模誤差后的均方根明顯降低,濾波前最大均方差為5.56 mm,過(guò)濾后為3.9 mm;濾波后北、東和高程方向平均均方根分別為0.78、0.56和2.89 mm,相對(duì)于過(guò)濾前分別減小了69%、60%和41%,空間濾波后陜西省GPS基準(zhǔn)站坐標(biāo)精確性和可靠性得到了較大提高。

    表1 濾波前后部分站點(diǎn)均方根比較(單位:mm)Tab.1 Comparison of RMS at some sites before and after filtering(unit:mm)

    圖1 區(qū)域疊加濾波提取的共模誤差Fig.1 Common mode error extracted by regional stacking

    3 空間濾波后時(shí)間序列的噪聲特征分析

    GPS坐標(biāo)時(shí)間序列中的噪聲功率譜可表示為[13]:

    其中,Px(f)是功率譜密度,f是頻率,P0和f0是待求參數(shù),κ是譜指數(shù)。在時(shí)域分析中,譜指數(shù)k能有效判斷噪聲特性,當(dāng)-3<k<-1時(shí)為非平穩(wěn)過(guò)程;當(dāng)-1<k<1時(shí)為平穩(wěn)過(guò)程,也稱為分形白噪聲。其中,k=0時(shí)為白噪聲;k=-1時(shí)為閃爍噪聲;k=-2時(shí)為隨機(jī)漫步噪聲。通過(guò)CATS軟件計(jì)算空間濾波后的陜西省GPS基準(zhǔn)站坐標(biāo)時(shí)間序列的譜指數(shù)如表2所示,23個(gè)基準(zhǔn)站各坐標(biāo)分量的噪聲譜指數(shù)基本介于-1~0,SNAK站E方向略小于-1。說(shuō)明這些基準(zhǔn)站各坐標(biāo)分量不僅有白噪聲,也有有色噪聲。

    表2 空間濾波后23個(gè)測(cè)站坐標(biāo)時(shí)間序列譜指數(shù)(單位:mm)Tab.2 Results of spectral index of 23 sites after filtering(unit:mm)

    大量研究表明,GPS坐標(biāo)時(shí)間序列噪聲一般都存在白噪聲(WN)和閃爍噪聲(FN)[3-10],有的存在隨機(jī)漫步噪聲(RN)[10,14]。將噪聲模型分為白噪聲、白噪聲+閃爍噪聲、白噪聲+隨機(jī)漫步噪聲、白噪聲+閃爍噪聲+隨機(jī)漫步噪聲4類,利用CATS軟件估計(jì)不同噪聲模型下最大似然值,以確定最佳噪聲模型。比較23個(gè)站N、E、U方向4種不同噪聲模型計(jì)算的MLE值發(fā)現(xiàn)(N方向MLE值見(jiàn)圖3),各站略有不同,但大部分站4個(gè)噪聲模型下的MLE值相近。為了估計(jì)各類噪聲大小,采用“白噪聲+閃爍噪聲+隨機(jī)漫步噪聲(WN+FN+RN)”來(lái)模擬時(shí)間序列中的噪聲。

    圖3 N方向4種不同噪聲模型下的最大似然值Fig.3 MLE of 4 noise model in N direction

    限于篇幅,表3僅列出10個(gè)站空間濾波前后基于WN+FN+RN模型的坐標(biāo)時(shí)間序列在北、東、高3個(gè)分量上噪聲參數(shù)估值。從表3可以看出:

    1)由于受到不同噪聲信號(hào)的影響,陜西省連續(xù)GPS基準(zhǔn)站水平和垂直方向表現(xiàn)出不同的噪聲性質(zhì)。空間濾波前絕大部分站表現(xiàn)為白噪聲+閃爍噪聲性質(zhì),空間濾波后各站點(diǎn)3個(gè)分量39.1%含有白噪聲+閃爍噪聲+隨機(jī)漫步噪聲,49.3%含有白噪聲+閃爍噪聲。

    2)三個(gè)坐標(biāo)分量具有不同的噪聲性質(zhì),垂直分量噪聲明顯大于水平分量,這與一般認(rèn)為的垂直分量精度要低于水平分量的結(jié)果一致。空間濾波前N方向的噪聲分量高于E方向(特別是閃爍噪聲),原因尚無(wú)法確定,可能與框架定義有關(guān)[2,15]。閃爍噪聲歸因于空間相關(guān)的物理性偏差[8],空間濾波可以分離空間相關(guān)噪聲,主要削弱閃爍噪聲[6,8],經(jīng)過(guò)濾波后N、E方向噪聲分量基本相當(dāng)。

    3)當(dāng)閃爍噪聲占主導(dǎo)地位時(shí),隨機(jī)漫步噪聲易被掩蓋[2,8],但原因不明。經(jīng)過(guò)空間濾波后的 GPS位置時(shí)間序列中閃爍噪聲明顯減小,部分站出現(xiàn)隨機(jī)漫步噪聲。一般認(rèn)為,隨機(jī)漫步噪聲源自天線墩的運(yùn)動(dòng),由氣候、土壤水分等的季節(jié)變化造成。天線墩通常由混凝土澆筑而成,在太陽(yáng)照射下可能產(chǎn)生日變化,同時(shí)可能還受到外部溫度、氣壓、濕度等環(huán)境因素的影響[8]。

    4)陜西省連續(xù)GPS基準(zhǔn)站均采用水泥柱天線墩,大部分建立在基巖之上,SNDL、SNHX、SNJY、SNYL由于地表較厚,采用20 m埋深水泥柱。比較土層水泥墩和基巖水泥墩的噪聲大小發(fā)現(xiàn),兩者噪聲水平無(wú)明顯差別,表明埋深20 m的水泥天線墩基本達(dá)到基巖水泥墩的穩(wěn)定性。

    表3 空間濾波前后基于WN+FN+RN模型的測(cè)站坐標(biāo)時(shí)間序列噪聲參數(shù)估值(單位:mm)Tab.3 Results of estimated parameters of noise property based on WN+FN+RN model before and after filtering(unit:mm)

    4 結(jié)論

    1)空間相關(guān)分析表明,在1 000 km之內(nèi),坐標(biāo)序列具有較高的相關(guān)性[8]。陜西省南北長(zhǎng)880 km,東西長(zhǎng)160~490 km,空間相關(guān)性高。區(qū)域疊加濾波去除共性粗差后,N、E、U方向的坐標(biāo)平均均方根相比濾波前分別減小了69%、60%和41%,基準(zhǔn)站坐標(biāo)精確性和可靠性得到很大提高。水平分量比垂向分量空間相關(guān)性更強(qiáng),共模誤差不僅含有季節(jié)性變化等影響(對(duì)垂向分量影響顯著),可能還含有衛(wèi)星軌道誤差、解算策略的影響,具體原因有待進(jìn)一步分析。

    2)譜指數(shù)計(jì)算發(fā)現(xiàn),基準(zhǔn)站坐標(biāo)時(shí)間序列不僅包括白噪聲,也包括有色噪聲。利用最大似然值估計(jì)各噪聲分量,表明空間濾波后各站點(diǎn)3個(gè)分量39.1%含有白噪聲+閃爍噪聲+隨機(jī)漫步噪聲,49.3% 含有白噪聲+閃爍噪聲,垂直分量噪聲明顯大于水平分量,東西分量與南北分量噪聲水平相當(dāng)。經(jīng)過(guò)空間濾波,白噪聲和閃爍噪聲(尤其是閃爍噪聲)分量減小,部分站出現(xiàn)隨機(jī)漫步噪聲。

    3)空間濾波前,水平分量各站點(diǎn)噪聲水平相當(dāng),垂向分量略有差異,表明水平方向比垂向噪聲表現(xiàn)更一致。經(jīng)濾波后,各分量還殘留一些未能被剔除的噪聲,特別是垂向分量還殘留較大的噪聲,一部分可能源于濾波方法未能徹底剔除的共模誤差,另一部分可能由于各站點(diǎn)自有的不具有共性的噪聲。

    1 王敏,沈正康,董大南.非構(gòu)造形變對(duì)GPS連續(xù)站位置時(shí)間序列的影響和修正[J].地球物理學(xué)報(bào),2005,48(5):1 045 - 1 052.(Wang Min,Shen Zhengkang,Dong Danan.Effects of non-tectonic crustal deformation on continuous GPS position time series and correction to them[J].Chinese J Geophy,2005,48(5):1 045 -1 052)

    2 田云鋒,沈正康.GPS坐標(biāo)時(shí)間序列中非構(gòu)造噪聲的剔除方法研究進(jìn)展[J].地震學(xué)報(bào),2009,31(1):68 -81.(Tian Yunfeng,Shen Zhengkang.Progress on reduction of non-tectonic noise in GPS position time series[J].Acta Seismologica Sinica,2009,31(1):68 -81)

    3 符養(yǎng).中國(guó)大陸現(xiàn)今地殼形變與GPS坐標(biāo)時(shí)間序列分析[D].上海:中國(guó)科學(xué)院上海天文臺(tái),2002.(Fu Yang.Present-day crustal deformation in China and GPS-derived coordinate time series analysis[D].Shanghai:Shanghai Astronomical Observatory,CAS,2002)

    4 朱文耀,符養(yǎng),李彥.GPS高程導(dǎo)出的全球高程振蕩運(yùn)動(dòng)及季節(jié)變化[J].中國(guó)科學(xué) D 輯,2003,33(5):470 -481.(Zhu Wenyao,F(xiàn)u Yang,Li Yan.Global height vibration and its seasonal variation induced by GPS height[J].Science in China:Series D,2003,33(5):470 -481)

    5 Mao A,Harrision C G A,Dixon T H.Noise in GPS coordinate time series[J].J Geophys Res,1999,104(B2):2 797-2 816.

    6 Nikolaidis R.Observation of geodetic and seismic deformation with the global positioning system[D].San Digeo:University of California,2002.

    7 Dong D.Spatiotemporal filtering using principal component analysis and Karhunen 2 loeve expansion approaches for regional GPS network analysis[J].Journal of Geophysical Research,2006,111(B3):3 405 -3 421.

    8 Williams S D P.Error analysis of continuous GPS position time series[J].Journal of Geophysical Research,2004,109(B3):412-430.

    9 Williams S D P.CATS:GPS coordinate time series analysis software[J].GPS Solutions,2008,12(2):147 -153.

    10 黃立人.GPS基準(zhǔn)站時(shí)間序列的噪聲特性分析[J].大地測(cè)量與地球動(dòng)力學(xué),2006(2):31-33.(Huang Liren.Noise properties in time series of coordinate component at GPS fiducial stations[J].Journal of Geodesy and Geodynamics,2006(2):31 -33)

    11 Wdowinski S.Southern California permanent GPS geodetic array:Spatial filtering of daily positions for estimating coseismic and postseismic displacements induced by the 1992 Landers earthquake[J].J Geophys Res,1997,102(B8):18 057-18 070.

    12 胡守超,伍吉倉(cāng),孫亞峰.區(qū)域GPS網(wǎng)三種時(shí)空濾波方法的比較[J].大地測(cè)量與地球動(dòng)力學(xué),2009(3):95-99.(Hu Shouchao,Wu Jicang,Sun Yafeng.Comparison among three spatiotemporal filtering methods for regional GPS networks analysis[J].Journal of Geodesy and Geodynamics,2009(3):95 -99)

    13 Agnew D C.The time-domain behavior of power-law noises[J].Geophys Res Lett,1992,19(4):333 -336.

    14 姜衛(wèi)平.利用連續(xù)GPS觀測(cè)數(shù)據(jù)分析水庫(kù)長(zhǎng)期變形[J].測(cè)繪學(xué)報(bào),2012,41(5):682 -689.(Jiang Weiping.Analysis of long-term deformation of reservoir using continuous GPS observations[J].Acta Geodaetica et Cartographica Sinica,2012,41(5):682 -689)

    15 蔣志浩.顧及有色噪聲影響的CGCS2000下我國(guó)CORS站速度估計(jì)[J].測(cè)繪學(xué)報(bào),2010,39(4):355-363.(Jiang Zhihao.Velocity estimation on the colored noise properties ofCORS network in Chinabased on the CGCS2000 frame[J].Acta Geodaetica et Cartographica Sinica,2010,39(4):355 -363)

    STUDY ON COORDINATE TIME SERIES OF SHAANXI CONTINUOUS GPS REFRENCE STATIONS

    Su Lina1,2),Ding Xiaoguang1,2),Zhang Yanfen3)and Zhang Yongqi1)
    1)Institute of Geology,CEA,Beijing 100029
    2)Earthquake Administration of Shaanxi Province,Xi’an 710068
    3)Surveying and Mapping Engineering Institute of Shanxi Province,Taiyuan030002

    The data of the Shaanxi continuous GPS reference stations from Jan 1st,2011 to May 30th,2013 are computed and analyzed in this article.After common mode error being removed by the regional stacking,the mean RMS of coordinate time series are declined 69%,60%and 41%respectively in N,E,U directions.The spectral indexes analysis states that the after-filtered coordinate time series contain the white noise and the color noise.The CATS software analysis indicates that the noise components of vertical direction are much larger than those of horizontal direction,and random noise appear in some stations along with flicker noise being decreased by the spatial filtering.

    Shaanxi continuous GPS reference stations;time series;regional stacking;common mode error;noise characteristics

    P228.41

    A

    1671-5942(2014)05-0106-04

    2013-11-11

    陜西省地震局啟航與創(chuàng)新基金項(xiàng)目(201305)。

    蘇利娜,女,1986年生,助工,主要從事GPS高精度解算和地殼形變分析研究。E-mail:sulinasgg@163.com。

    猜你喜歡
    共模漫步基準(zhǔn)
    海邊漫步
    漫步春天
    月下漫步
    小讀者(2019年24期)2019-01-10 23:00:37
    關(guān)于差模和共模干擾的研究
    憶中傘
    明基準(zhǔn)講方法??待R
    滑落還是攀爬
    非隔離型光伏并網(wǎng)逆變器共模電流分析
    單相逆變器共模電磁干擾特性研究
    巧用基準(zhǔn)變換實(shí)現(xiàn)裝配檢測(cè)
    河南科技(2014年15期)2014-02-27 14:12:35
    各种免费的搞黄视频| 国产精品女同一区二区软件| 大片免费播放器 马上看| 国产精品一二三区在线看| 精品久久久噜噜| 亚洲无线观看免费| av有码第一页| 蜜臀久久99精品久久宅男| 2018国产大陆天天弄谢| 丝袜脚勾引网站| 久久精品国产鲁丝片午夜精品| 欧美人与善性xxx| 男女免费视频国产| 九色成人免费人妻av| 春色校园在线视频观看| 一级爰片在线观看| 国产免费一级a男人的天堂| 美女国产视频在线观看| 亚洲高清免费不卡视频| 国产日韩欧美亚洲二区| 国产精品一区二区三区四区免费观看| 免费播放大片免费观看视频在线观看| 美女内射精品一级片tv| 高清在线视频一区二区三区| 免费av不卡在线播放| 婷婷色综合大香蕉| 亚洲精品色激情综合| 国产精品久久久久久av不卡| av线在线观看网站| 最近中文字幕2019免费版| 成人黄色视频免费在线看| 我要看日韩黄色一级片| 国产亚洲精品久久久com| 蜜桃久久精品国产亚洲av| 国产探花极品一区二区| 国产极品天堂在线| 偷拍熟女少妇极品色| 国产黄色视频一区二区在线观看| kizo精华| 亚洲熟女精品中文字幕| 91久久精品国产一区二区成人| 国产69精品久久久久777片| 亚洲人成网站在线观看播放| 国产精品成人在线| 丰满迷人的少妇在线观看| 亚洲高清免费不卡视频| 丝袜脚勾引网站| 夜夜爽夜夜爽视频| 久久久亚洲精品成人影院| 色哟哟·www| 男人爽女人下面视频在线观看| 亚洲不卡免费看| 人妻一区二区av| 2022亚洲国产成人精品| 色视频在线一区二区三区| 韩国高清视频一区二区三区| 国产在线视频一区二区| 国产精品秋霞免费鲁丝片| 中文字幕av电影在线播放| 国产69精品久久久久777片| 亚洲国产精品999| 亚洲精品日韩在线中文字幕| 国产综合精华液| 韩国高清视频一区二区三区| 国产在线免费精品| 色吧在线观看| 久久久久精品性色| 亚洲无线观看免费| 少妇人妻 视频| 美女福利国产在线| 国产精品麻豆人妻色哟哟久久| 欧美激情极品国产一区二区三区 | 精品亚洲成国产av| 亚洲性久久影院| 亚洲欧美精品自产自拍| 亚洲精品,欧美精品| 有码 亚洲区| 亚洲欧美中文字幕日韩二区| 蜜桃在线观看..| 日日啪夜夜撸| 久久精品久久精品一区二区三区| 国产爽快片一区二区三区| a级片在线免费高清观看视频| 国产有黄有色有爽视频| 有码 亚洲区| 少妇人妻精品综合一区二区| 亚洲久久久国产精品| 这个男人来自地球电影免费观看 | 美女内射精品一级片tv| 久久狼人影院| 两个人的视频大全免费| 99九九线精品视频在线观看视频| 久久久久久久大尺度免费视频| 国产成人精品一,二区| 女人久久www免费人成看片| 边亲边吃奶的免费视频| 麻豆成人av视频| 伊人亚洲综合成人网| 亚洲欧洲日产国产| 亚洲精品中文字幕在线视频 | 国产一区亚洲一区在线观看| 国产黄片视频在线免费观看| 国产免费一区二区三区四区乱码| 99久久综合免费| 一本久久精品| 妹子高潮喷水视频| av黄色大香蕉| 久久久精品免费免费高清| 国产伦在线观看视频一区| 黄色配什么色好看| 你懂的网址亚洲精品在线观看| a级毛片在线看网站| 精品久久久久久电影网| 91久久精品国产一区二区三区| 久久国产亚洲av麻豆专区| 美女cb高潮喷水在线观看| 热re99久久精品国产66热6| 人妻 亚洲 视频| 免费人妻精品一区二区三区视频| 国产高清有码在线观看视频| 日日啪夜夜爽| 男女啪啪激烈高潮av片| 91久久精品国产一区二区成人| 伊人久久国产一区二区| 天天躁夜夜躁狠狠久久av| 亚洲精品色激情综合| 最后的刺客免费高清国语| tube8黄色片| 免费观看av网站的网址| 久久午夜综合久久蜜桃| 久久久久久久久久久丰满| 91精品一卡2卡3卡4卡| 纵有疾风起免费观看全集完整版| 亚洲精品久久午夜乱码| 一个人免费看片子| 人妻制服诱惑在线中文字幕| 日日撸夜夜添| 熟女人妻精品中文字幕| 国产精品.久久久| 七月丁香在线播放| 久久99热6这里只有精品| 久久鲁丝午夜福利片| 亚洲av成人精品一二三区| 国产欧美日韩精品一区二区| 伊人亚洲综合成人网| 成人漫画全彩无遮挡| 午夜精品国产一区二区电影| 三级国产精品欧美在线观看| 久久99热这里只频精品6学生| 黑人巨大精品欧美一区二区蜜桃 | 国产伦在线观看视频一区| 亚洲美女视频黄频| 久久久久网色| 大香蕉久久网| 2018国产大陆天天弄谢| 国产爽快片一区二区三区| 亚洲av二区三区四区| 爱豆传媒免费全集在线观看| 久热久热在线精品观看| 永久免费av网站大全| 丰满迷人的少妇在线观看| 亚洲av综合色区一区| 精品一区二区三卡| 永久免费av网站大全| 亚洲情色 制服丝袜| 久久久久久久大尺度免费视频| 在线天堂最新版资源| 夫妻午夜视频| 久久精品国产鲁丝片午夜精品| 天天操日日干夜夜撸| 91久久精品国产一区二区三区| 大片免费播放器 马上看| 欧美 日韩 精品 国产| 一个人看视频在线观看www免费| 亚洲欧美成人综合另类久久久| 色婷婷久久久亚洲欧美| 看免费成人av毛片| 亚洲欧洲日产国产| av黄色大香蕉| av播播在线观看一区| av国产久精品久网站免费入址| 国产欧美日韩精品一区二区| 国产成人91sexporn| 黄色视频在线播放观看不卡| 三级经典国产精品| 狂野欧美激情性xxxx在线观看| 国产爽快片一区二区三区| 人妻夜夜爽99麻豆av| 美女xxoo啪啪120秒动态图| 久久久久久人妻| 高清黄色对白视频在线免费看 | 国产精品一区二区在线观看99| 啦啦啦在线观看免费高清www| 人妻系列 视频| 国产男女超爽视频在线观看| 9色porny在线观看| 丁香六月天网| 中文字幕av电影在线播放| 国产精品久久久久久久久免| av网站免费在线观看视频| 久久久久久久精品精品| 黄色欧美视频在线观看| 亚洲欧洲国产日韩| 国产美女午夜福利| 国产熟女欧美一区二区| 精品少妇久久久久久888优播| 新久久久久国产一级毛片| 日韩av在线免费看完整版不卡| 国产成人免费无遮挡视频| 成年女人在线观看亚洲视频| av免费观看日本| 2021少妇久久久久久久久久久| 熟女av电影| 亚洲伊人久久精品综合| 国产无遮挡羞羞视频在线观看| 国产黄片美女视频| 久久久久久久久大av| 一区二区三区乱码不卡18| 精品亚洲成国产av| 久久久久久久久久久久大奶| 欧美国产精品一级二级三级 | 97超视频在线观看视频| 国产精品嫩草影院av在线观看| 久热这里只有精品99| 国内揄拍国产精品人妻在线| 午夜福利视频精品| 日韩欧美一区视频在线观看 | 精品久久久噜噜| 免费观看a级毛片全部| 一本一本综合久久| 日韩,欧美,国产一区二区三区| 亚洲精品乱久久久久久| 麻豆成人午夜福利视频| 久久久欧美国产精品| 成人黄色视频免费在线看| 亚洲成人手机| 亚洲av成人精品一区久久| 汤姆久久久久久久影院中文字幕| 欧美性感艳星| 女的被弄到高潮叫床怎么办| 色94色欧美一区二区| www.色视频.com| 成人无遮挡网站| 美女中出高潮动态图| 日日啪夜夜撸| 免费观看无遮挡的男女| 欧美成人午夜免费资源| 亚洲av免费高清在线观看| 高清毛片免费看| 国产精品久久久久久久电影| 搡女人真爽免费视频火全软件| 啦啦啦中文免费视频观看日本| 国产精品国产三级国产专区5o| 亚洲性久久影院| 日韩制服骚丝袜av| 国产成人aa在线观看| 岛国毛片在线播放| 亚洲av男天堂| 亚洲欧美清纯卡通| 极品少妇高潮喷水抽搐| 久久人妻熟女aⅴ| 最近最新中文字幕免费大全7| 久久6这里有精品| 亚洲成人一二三区av| 国产成人午夜福利电影在线观看| 亚洲精品乱码久久久久久按摩| 亚洲精品国产色婷婷电影| 日韩免费高清中文字幕av| 99久国产av精品国产电影| 人人澡人人妻人| 极品少妇高潮喷水抽搐| 美女大奶头黄色视频| 亚洲精品久久午夜乱码| 精品亚洲乱码少妇综合久久| 中文字幕人妻丝袜制服| 多毛熟女@视频| 啦啦啦啦在线视频资源| 高清午夜精品一区二区三区| 久久久国产一区二区| 狂野欧美白嫩少妇大欣赏| 亚洲无线观看免费| 如日韩欧美国产精品一区二区三区 | 国国产精品蜜臀av免费| 另类精品久久| 久久久国产精品麻豆| xxx大片免费视频| 啦啦啦啦在线视频资源| 国产精品一区二区性色av| 一级二级三级毛片免费看| 丝袜在线中文字幕| 一级毛片黄色毛片免费观看视频| 伦理电影免费视频| 啦啦啦啦在线视频资源| 最近中文字幕2019免费版| 午夜激情久久久久久久| 我要看黄色一级片免费的| 国产免费又黄又爽又色| 国产 精品1| 一级毛片aaaaaa免费看小| 国产成人精品无人区| 大片电影免费在线观看免费| 老司机影院成人| a级一级毛片免费在线观看| 国产精品国产三级专区第一集| 在线观看美女被高潮喷水网站| 亚洲欧美一区二区三区国产| 日日摸夜夜添夜夜添av毛片| 日韩伦理黄色片| 欧美日韩av久久| 黄色配什么色好看| 久久青草综合色| 亚洲精品一二三| 国产精品免费大片| h视频一区二区三区| 97超视频在线观看视频| 久久久久久久亚洲中文字幕| 免费不卡的大黄色大毛片视频在线观看| 内地一区二区视频在线| 国产精品秋霞免费鲁丝片| 久久久久久久大尺度免费视频| 日本与韩国留学比较| 少妇人妻精品综合一区二区| 亚洲av国产av综合av卡| 少妇人妻精品综合一区二区| freevideosex欧美| 狂野欧美白嫩少妇大欣赏| 国产女主播在线喷水免费视频网站| 日韩免费高清中文字幕av| www.av在线官网国产| 色哟哟·www| 精品亚洲乱码少妇综合久久| 日本vs欧美在线观看视频 | 日本vs欧美在线观看视频 | 一级a做视频免费观看| 中国国产av一级| 亚洲av在线观看美女高潮| 久久99精品国语久久久| 韩国高清视频一区二区三区| freevideosex欧美| 久久精品久久久久久久性| 亚洲精品日韩在线中文字幕| 国产精品国产三级国产专区5o| 男女啪啪激烈高潮av片| 少妇熟女欧美另类| 久久久午夜欧美精品| 午夜激情福利司机影院| 自拍偷自拍亚洲精品老妇| 免费黄色在线免费观看| 日韩人妻高清精品专区| 精品卡一卡二卡四卡免费| 99re6热这里在线精品视频| 9色porny在线观看| 精品人妻熟女毛片av久久网站| 国产亚洲午夜精品一区二区久久| 国产精品久久久久久久久免| 亚洲国产精品一区二区三区在线| 在线观看美女被高潮喷水网站| 亚洲国产毛片av蜜桃av| 午夜av观看不卡| 日韩欧美精品免费久久| 久久这里有精品视频免费| 精品久久久久久电影网| 秋霞在线观看毛片| 人体艺术视频欧美日本| 亚洲国产最新在线播放| 午夜福利在线观看免费完整高清在| 一本大道久久a久久精品| 在线天堂最新版资源| 亚洲成人手机| 人妻一区二区av| 国产国拍精品亚洲av在线观看| 亚洲精品国产av蜜桃| 日本黄大片高清| 久久 成人 亚洲| 免费人成在线观看视频色| 看非洲黑人一级黄片| 一本大道久久a久久精品| 国产精品不卡视频一区二区| 一级毛片 在线播放| 最近2019中文字幕mv第一页| 久久久国产精品麻豆| 激情五月婷婷亚洲| 亚洲人成网站在线播| 日韩av在线免费看完整版不卡| 亚洲高清免费不卡视频| 亚洲,欧美,日韩| 日韩成人av中文字幕在线观看| 内射极品少妇av片p| √禁漫天堂资源中文www| 少妇精品久久久久久久| 免费久久久久久久精品成人欧美视频 | 我的女老师完整版在线观看| 99久久人妻综合| 亚洲美女黄色视频免费看| 午夜激情福利司机影院| 如日韩欧美国产精品一区二区三区 | 久久精品熟女亚洲av麻豆精品| 国产国拍精品亚洲av在线观看| 成年女人在线观看亚洲视频| av天堂久久9| 国产伦理片在线播放av一区| 日日摸夜夜添夜夜爱| 亚洲成色77777| 日本爱情动作片www.在线观看| 在线观看人妻少妇| 自拍偷自拍亚洲精品老妇| 日本欧美国产在线视频| 九草在线视频观看| 中国三级夫妇交换| 久久久久国产精品人妻一区二区| 精品少妇黑人巨大在线播放| 一本色道久久久久久精品综合| 最近的中文字幕免费完整| 人人妻人人澡人人爽人人夜夜| 国产亚洲午夜精品一区二区久久| 欧美3d第一页| av天堂中文字幕网| 黑人猛操日本美女一级片| 内地一区二区视频在线| 亚洲一区二区三区欧美精品| 嫩草影院入口| 少妇被粗大的猛进出69影院 | 极品教师在线视频| 免费观看av网站的网址| av卡一久久| 成人亚洲欧美一区二区av| 国内精品宾馆在线| 亚洲精品久久久久久婷婷小说| 国产精品久久久久久精品古装| 少妇丰满av| 中文在线观看免费www的网站| 777米奇影视久久| 秋霞伦理黄片| 日本av手机在线免费观看| 国产在线免费精品| 免费观看的影片在线观看| 亚洲精品aⅴ在线观看| 韩国av在线不卡| 欧美日韩视频精品一区| 十八禁高潮呻吟视频 | 国产伦在线观看视频一区| 亚洲,一卡二卡三卡| 免费av不卡在线播放| 国产片特级美女逼逼视频| 我的老师免费观看完整版| 狂野欧美白嫩少妇大欣赏| 国产精品无大码| 亚洲不卡免费看| 五月玫瑰六月丁香| 久久亚洲国产成人精品v| 一区二区三区四区激情视频| 全区人妻精品视频| 在线免费观看不下载黄p国产| 国产 一区精品| 99久久人妻综合| 欧美精品国产亚洲| 麻豆成人av视频| 国产成人a∨麻豆精品| 国产av码专区亚洲av| 成年人午夜在线观看视频| 汤姆久久久久久久影院中文字幕| 亚洲人成网站在线播| 各种免费的搞黄视频| 国产av一区二区精品久久| 久久综合国产亚洲精品| 国产高清不卡午夜福利| 水蜜桃什么品种好| 日韩人妻高清精品专区| 欧美日韩av久久| 国产永久视频网站| 国产淫片久久久久久久久| 欧美三级亚洲精品| 热re99久久国产66热| 人妻一区二区av| av天堂久久9| 欧美变态另类bdsm刘玥| 一边亲一边摸免费视频| 在线观看人妻少妇| www.色视频.com| 老熟女久久久| 人人妻人人澡人人看| 久久久久久久国产电影| 日韩中字成人| 香蕉精品网在线| 免费观看a级毛片全部| 国产精品久久久久久精品古装| av播播在线观看一区| 麻豆精品久久久久久蜜桃| 亚洲精品视频女| 亚洲人与动物交配视频| 日韩一本色道免费dvd| 熟女人妻精品中文字幕| 日韩av不卡免费在线播放| 黄色视频在线播放观看不卡| 日韩欧美 国产精品| 亚洲国产色片| 国产免费福利视频在线观看| 国产免费一区二区三区四区乱码| av在线播放精品| 国产黄片美女视频| 黄色视频在线播放观看不卡| 少妇被粗大猛烈的视频| 在线观看一区二区三区激情| 一级二级三级毛片免费看| 欧美最新免费一区二区三区| 最近的中文字幕免费完整| 少妇丰满av| 精品一区二区三卡| 久久99热这里只频精品6学生| 肉色欧美久久久久久久蜜桃| 麻豆精品久久久久久蜜桃| 国产一区二区在线观看av| 国产伦精品一区二区三区四那| 国产av国产精品国产| 日本-黄色视频高清免费观看| 99热6这里只有精品| 久久精品久久久久久噜噜老黄| 亚洲三级黄色毛片| 啦啦啦中文免费视频观看日本| 午夜福利在线观看免费完整高清在| 国产视频内射| 一本—道久久a久久精品蜜桃钙片| 日韩在线高清观看一区二区三区| 99久久综合免费| 日韩人妻高清精品专区| 高清毛片免费看| 国产精品三级大全| 高清av免费在线| 伊人久久精品亚洲午夜| 亚洲精品日本国产第一区| 日本色播在线视频| 亚洲三级黄色毛片| 亚洲欧美日韩卡通动漫| 成人毛片60女人毛片免费| 国产一级毛片在线| 中文精品一卡2卡3卡4更新| 不卡视频在线观看欧美| 日韩亚洲欧美综合| 美女主播在线视频| 国产欧美日韩精品一区二区| 美女福利国产在线| 亚洲精品亚洲一区二区| 又爽又黄a免费视频| 日韩强制内射视频| 人体艺术视频欧美日本| 性色av一级| 国产欧美另类精品又又久久亚洲欧美| 亚洲在久久综合| 一本一本综合久久| 国产一区二区在线观看日韩| 久久这里有精品视频免费| 国产成人a∨麻豆精品| 亚洲精品自拍成人| 97在线人人人人妻| 一级二级三级毛片免费看| 夜夜看夜夜爽夜夜摸| 亚洲av综合色区一区| 亚洲av福利一区| 国产欧美亚洲国产| 最近2019中文字幕mv第一页| 日本vs欧美在线观看视频 | 男男h啪啪无遮挡| 亚洲激情五月婷婷啪啪| 中文字幕人妻熟人妻熟丝袜美| 偷拍熟女少妇极品色| 制服丝袜香蕉在线| 国产91av在线免费观看| 91久久精品国产一区二区三区| 人体艺术视频欧美日本| 国产精品久久久久久精品古装| 午夜激情久久久久久久| 看十八女毛片水多多多| 51国产日韩欧美| 久久久久久久久大av| 麻豆成人av视频| 成年人午夜在线观看视频| 国内揄拍国产精品人妻在线| 国产精品人妻久久久久久| 春色校园在线视频观看| 日韩中字成人| 久久毛片免费看一区二区三区| 亚洲综合精品二区| 日本色播在线视频| 成人午夜精彩视频在线观看| 国产免费一区二区三区四区乱码| av播播在线观看一区| 美女中出高潮动态图| 精品少妇内射三级| 久久99蜜桃精品久久| 男人添女人高潮全过程视频| 亚洲高清免费不卡视频| 日本与韩国留学比较| 欧美97在线视频| 亚洲婷婷狠狠爱综合网| 亚洲美女黄色视频免费看| 免费人成在线观看视频色| 精品国产国语对白av| 国产真实伦视频高清在线观看| 男人添女人高潮全过程视频| 嫩草影院入口| 91aial.com中文字幕在线观看| 99久久人妻综合| 色哟哟·www| 国产黄片美女视频| 欧美丝袜亚洲另类| 国产精品三级大全| 一级毛片黄色毛片免费观看视频| 如何舔出高潮| 成人影院久久| 亚洲av日韩在线播放| av线在线观看网站| 色视频在线一区二区三区| a级毛片在线看网站| 亚洲国产毛片av蜜桃av| 久久韩国三级中文字幕|