• <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
    少妇裸体淫交视频免费看高清| 人妻夜夜爽99麻豆av| 亚洲真实伦在线观看| 有码 亚洲区| 美女cb高潮喷水在线观看| 免费看日本二区| 日韩一区二区视频免费看| 精品熟女少妇av免费看| 最新中文字幕久久久久| 亚洲婷婷狠狠爱综合网| 晚上一个人看的免费电影| 在线播放无遮挡| 国内少妇人妻偷人精品xxx网站| 日本wwww免费看| 国产精品精品国产色婷婷| 一本久久精品| 国产色爽女视频免费观看| 人人妻人人澡人人爽人人夜夜| 97人妻精品一区二区三区麻豆| 亚洲国产欧美人成| 一级片'在线观看视频| 搞女人的毛片| 免费大片黄手机在线观看| 免费av不卡在线播放| 伊人久久精品亚洲午夜| 欧美老熟妇乱子伦牲交| 亚洲国产日韩一区二区| 看十八女毛片水多多多| 大香蕉久久网| 插逼视频在线观看| 一个人看的www免费观看视频| 你懂的网址亚洲精品在线观看| 免费观看av网站的网址| 欧美精品人与动牲交sv欧美| 国产真实伦视频高清在线观看| 一本色道久久久久久精品综合| 超碰av人人做人人爽久久| 亚洲色图综合在线观看| 亚洲精品日本国产第一区| 国产69精品久久久久777片| 别揉我奶头 嗯啊视频| 直男gayav资源| 国产日韩欧美在线精品| 国产精品国产av在线观看| 日韩,欧美,国产一区二区三区| 亚洲国产精品999| 精品酒店卫生间| 你懂的网址亚洲精品在线观看| 亚洲一区二区三区欧美精品 | 91精品伊人久久大香线蕉| 亚洲av免费高清在线观看| av黄色大香蕉| 精品久久久久久久人妻蜜臀av| 麻豆久久精品国产亚洲av| 亚洲自偷自拍三级| 免费黄色在线免费观看| 搞女人的毛片| 免费av毛片视频| 亚洲高清免费不卡视频| av黄色大香蕉| 国产永久视频网站| 亚洲精品国产色婷婷电影| 久久精品国产亚洲网站| 色综合色国产| 免费不卡的大黄色大毛片视频在线观看| 亚洲va在线va天堂va国产| 亚洲不卡免费看| 三级经典国产精品| 一本色道久久久久久精品综合| 国产亚洲av嫩草精品影院| 免费黄频网站在线观看国产| 国产午夜福利久久久久久| 亚洲av免费在线观看| 亚洲av免费高清在线观看| 成人综合一区亚洲| 国内精品美女久久久久久| 3wmmmm亚洲av在线观看| 高清av免费在线| a级毛片免费高清观看在线播放| 丰满人妻一区二区三区视频av| 一本色道久久久久久精品综合| 高清午夜精品一区二区三区| 国产黄色免费在线视频| 自拍偷自拍亚洲精品老妇| 少妇人妻 视频| 自拍欧美九色日韩亚洲蝌蚪91 | 中文欧美无线码| 99久久中文字幕三级久久日本| 一本久久精品| 欧美一级a爱片免费观看看| 国产精品一区二区性色av| 有码 亚洲区| 97精品久久久久久久久久精品| 欧美高清成人免费视频www| 国产爱豆传媒在线观看| 久久人人爽av亚洲精品天堂 | 亚洲欧美日韩东京热| 又大又黄又爽视频免费| 禁无遮挡网站| 国产成人午夜福利电影在线观看| 国内少妇人妻偷人精品xxx网站| 视频区图区小说| 听说在线观看完整版免费高清| 国产亚洲一区二区精品| av天堂中文字幕网| 亚州av有码| 亚洲国产精品成人久久小说| 中文资源天堂在线| 哪个播放器可以免费观看大片| 久久97久久精品| a级毛色黄片| 久久精品熟女亚洲av麻豆精品| 看十八女毛片水多多多| 国产亚洲5aaaaa淫片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美一区二区亚洲| 亚洲高清免费不卡视频| 欧美xxxx性猛交bbbb| 欧美成人午夜免费资源| 波多野结衣巨乳人妻| av福利片在线观看| 久久精品人妻少妇| 视频区图区小说| 2018国产大陆天天弄谢| 中国国产av一级| 欧美成人一区二区免费高清观看| 国产v大片淫在线免费观看| 大香蕉久久网| 黑人高潮一二区| 51国产日韩欧美| 伦精品一区二区三区| 精品国产乱码久久久久久小说| 亚洲精品影视一区二区三区av| 国产高清三级在线| 国产精品无大码| 国产免费视频播放在线视频| 国产色爽女视频免费观看| 国产爱豆传媒在线观看| 一个人观看的视频www高清免费观看| 最近最新中文字幕大全电影3| 三级国产精品欧美在线观看| 国产精品精品国产色婷婷| 亚洲精品国产av蜜桃| 插逼视频在线观看| av在线天堂中文字幕| 又爽又黄无遮挡网站| 久久女婷五月综合色啪小说 | 久久精品国产鲁丝片午夜精品| av在线蜜桃| 久久影院123| 亚洲国产成人一精品久久久| 成年女人看的毛片在线观看| 在线观看国产h片| 女的被弄到高潮叫床怎么办| 肉色欧美久久久久久久蜜桃 | 国产精品一区二区三区四区免费观看| 啦啦啦在线观看免费高清www| 国产人妻一区二区三区在| 亚洲性久久影院| 日韩制服骚丝袜av| 人妻少妇偷人精品九色| 国产精品一区二区在线观看99| av又黄又爽大尺度在线免费看| 免费黄色在线免费观看| 九九在线视频观看精品| 久久99热这里只频精品6学生| av黄色大香蕉| 久久99精品国语久久久| 免费观看的影片在线观看| 18+在线观看网站| 少妇丰满av| 一本一本综合久久| 99久久九九国产精品国产免费| 国产精品国产三级国产av玫瑰| 97精品久久久久久久久久精品| 国产精品一区二区三区四区免费观看| 观看免费一级毛片| 精品人妻视频免费看| 亚洲国产最新在线播放| 色5月婷婷丁香| 精品久久久久久电影网| 久久久久性生活片| 日本黄大片高清| 久久综合国产亚洲精品| 国产日韩欧美在线精品| 久久精品久久久久久噜噜老黄| 亚洲av免费高清在线观看| av在线亚洲专区| 菩萨蛮人人尽说江南好唐韦庄| av女优亚洲男人天堂| 人妻制服诱惑在线中文字幕| 中文精品一卡2卡3卡4更新| 麻豆精品久久久久久蜜桃| h日本视频在线播放| 2018国产大陆天天弄谢| 久久精品国产鲁丝片午夜精品| 久久久亚洲精品成人影院| 女人十人毛片免费观看3o分钟| 一级爰片在线观看| 久久精品国产自在天天线| 日本一二三区视频观看| 亚洲在线观看片| 精品国产露脸久久av麻豆| 少妇人妻 视频| 毛片一级片免费看久久久久| 91精品国产九色| 国产成人精品久久久久久| 青青草视频在线视频观看| 在线观看人妻少妇| 最新中文字幕久久久久| 国产精品一二三区在线看| 99视频精品全部免费 在线| a级毛片免费高清观看在线播放| 亚洲久久久久久中文字幕| 交换朋友夫妻互换小说| 欧美成人一区二区免费高清观看| 日韩成人伦理影院| 国产精品人妻久久久影院| 国产精品伦人一区二区| 亚洲激情五月婷婷啪啪| 日本欧美国产在线视频| 黄色日韩在线| 听说在线观看完整版免费高清| 乱码一卡2卡4卡精品| 亚洲欧美精品专区久久| 国产黄色视频一区二区在线观看| 欧美区成人在线视频| 国产精品女同一区二区软件| 在线观看一区二区三区激情| 欧美bdsm另类| 热re99久久精品国产66热6| 人体艺术视频欧美日本| 日本av手机在线免费观看| 国产一区有黄有色的免费视频| 国产精品久久久久久久电影| 亚洲伊人久久精品综合| 丝袜美腿在线中文| 晚上一个人看的免费电影| 内射极品少妇av片p| 国产永久视频网站| 最近最新中文字幕免费大全7| 黄色怎么调成土黄色| 超碰97精品在线观看| 久久精品熟女亚洲av麻豆精品| av.在线天堂| 日韩精品有码人妻一区| 免费观看性生交大片5| 99热这里只有是精品在线观看| 人人妻人人看人人澡| 久久久久精品久久久久真实原创| 午夜福利视频精品| 老师上课跳d突然被开到最大视频| 国产精品一区二区在线观看99| 精品久久国产蜜桃| 精品久久久噜噜| 日韩伦理黄色片| 免费黄色在线免费观看| av网站免费在线观看视频| 欧美性猛交╳xxx乱大交人| 国产av国产精品国产| 精品国产乱码久久久久久小说| 亚洲,一卡二卡三卡| 免费黄频网站在线观看国产| 国内少妇人妻偷人精品xxx网站| 一级毛片电影观看| 女的被弄到高潮叫床怎么办| 熟女电影av网| 亚洲精品色激情综合| 亚洲欧美一区二区三区国产| 在线精品无人区一区二区三 | 国产精品av视频在线免费观看| 99久国产av精品国产电影| 国产一区有黄有色的免费视频| av网站免费在线观看视频| 黄色怎么调成土黄色| 国产日韩欧美在线精品| 亚洲欧美日韩无卡精品| 2021少妇久久久久久久久久久| 男女无遮挡免费网站观看| 国内精品美女久久久久久| 黄色怎么调成土黄色| 国产淫语在线视频| 亚洲真实伦在线观看| 欧美激情在线99| 成人欧美大片| 免费看a级黄色片| 男男h啪啪无遮挡| 男人添女人高潮全过程视频| 自拍偷自拍亚洲精品老妇| 成人无遮挡网站| 久久久久久久亚洲中文字幕| 成人毛片a级毛片在线播放| 亚洲精品国产成人久久av| 欧美日韩亚洲高清精品| 婷婷色麻豆天堂久久| 毛片女人毛片| 熟女电影av网| 国产在视频线精品| 高清午夜精品一区二区三区| videossex国产| 国产成人精品福利久久| 一区二区三区四区激情视频| 国产精品99久久久久久久久| 18禁在线播放成人免费| 人妻一区二区av| 成人特级av手机在线观看| 国产色婷婷99| 五月玫瑰六月丁香| 天堂网av新在线| 亚洲伊人久久精品综合| 伊人久久国产一区二区| 久久99热这里只有精品18| 国产精品蜜桃在线观看| 亚洲一区二区三区欧美精品 | 97在线人人人人妻| 久久精品久久久久久久性| 久久鲁丝午夜福利片| 精品一区二区免费观看| 国产老妇伦熟女老妇高清| 91狼人影院| 中文欧美无线码| 在线观看美女被高潮喷水网站| 国产精品久久久久久精品古装| 天天一区二区日本电影三级| 欧美bdsm另类| 亚洲在线观看片| 午夜福利网站1000一区二区三区| 激情五月婷婷亚洲| 干丝袜人妻中文字幕| 亚洲成色77777| 色吧在线观看| 精品一区二区三区视频在线| 国产伦理片在线播放av一区| 亚洲成人中文字幕在线播放| 亚洲va在线va天堂va国产| 亚洲最大成人中文| 欧美激情在线99| 国产一区二区三区综合在线观看 | 亚洲欧美清纯卡通| 天堂中文最新版在线下载 | 国产毛片a区久久久久| 我要看日韩黄色一级片| 欧美日本视频| 18禁在线播放成人免费| 久久精品国产亚洲av涩爱| 久久综合国产亚洲精品| 国产一区二区三区综合在线观看 | 内地一区二区视频在线| 日本wwww免费看| 在线亚洲精品国产二区图片欧美 | 精品视频人人做人人爽| 精品久久久久久久久av| 人人妻人人看人人澡| 成人特级av手机在线观看| 亚洲av中文字字幕乱码综合| 一级毛片电影观看| 久久久欧美国产精品| 亚洲欧美成人精品一区二区| 一边亲一边摸免费视频| 秋霞伦理黄片| 亚洲最大成人av| 国产又色又爽无遮挡免| 欧美性感艳星| 内射极品少妇av片p| 日本av手机在线免费观看| 在线精品无人区一区二区三 | 免费看日本二区| 欧美日韩视频精品一区| 国产精品蜜桃在线观看| 草草在线视频免费看| 久久久久久久国产电影| 国产成人91sexporn| 尤物成人国产欧美一区二区三区| 免费电影在线观看免费观看| 欧美一级a爱片免费观看看| 亚洲av福利一区| 国产探花在线观看一区二区| av一本久久久久| kizo精华| 久久精品久久精品一区二区三区| 三级经典国产精品| 亚洲va在线va天堂va国产| 欧美日韩一区二区视频在线观看视频在线 | 亚洲四区av| 精品熟女少妇av免费看| 精品人妻熟女av久视频| 亚洲真实伦在线观看| 成人鲁丝片一二三区免费| 亚洲自拍偷在线| 香蕉精品网在线| 黄片wwwwww| 亚洲精品国产av蜜桃| 视频区图区小说| 亚洲av在线观看美女高潮| 久久久久久伊人网av| 新久久久久国产一级毛片| 中国三级夫妇交换| 男人添女人高潮全过程视频| 日韩强制内射视频| 内射极品少妇av片p| 久久6这里有精品| 亚洲成色77777| 亚洲va在线va天堂va国产| 永久网站在线| 少妇人妻精品综合一区二区| 91aial.com中文字幕在线观看| 欧美成人精品欧美一级黄| 26uuu在线亚洲综合色| 别揉我奶头 嗯啊视频| 最近最新中文字幕免费大全7| 天堂网av新在线| 精品久久国产蜜桃| 伦精品一区二区三区| 欧美激情在线99| 夫妻性生交免费视频一级片| 热99国产精品久久久久久7| 国产免费一级a男人的天堂| 国产黄色视频一区二区在线观看| 街头女战士在线观看网站| 中文字幕制服av| 亚洲伊人久久精品综合| 青春草视频在线免费观看| 男女国产视频网站| 国产一区二区三区综合在线观看 | 国产色爽女视频免费观看| 又粗又硬又长又爽又黄的视频| 久久久精品免费免费高清| 亚洲av电影在线观看一区二区三区 | 亚洲av福利一区| 亚洲av中文字字幕乱码综合| 校园人妻丝袜中文字幕| 国产老妇女一区| 大香蕉97超碰在线| 91午夜精品亚洲一区二区三区| 日本黄色片子视频| 亚洲怡红院男人天堂| 一级爰片在线观看| 舔av片在线| 九草在线视频观看| 国产亚洲精品久久久com| 午夜福利高清视频| 干丝袜人妻中文字幕| 精品国产三级普通话版| 国产成人精品一,二区| 成人特级av手机在线观看| 搞女人的毛片| 性色avwww在线观看| 成人毛片60女人毛片免费| 亚洲精品第二区| 亚洲成人精品中文字幕电影| 伊人久久国产一区二区| 人妻一区二区av| 午夜激情久久久久久久| 一级毛片aaaaaa免费看小| av在线天堂中文字幕| 伦理电影大哥的女人| 性色avwww在线观看| 国产高潮美女av| 免费在线观看成人毛片| 亚洲欧洲国产日韩| 男人添女人高潮全过程视频| 在线观看美女被高潮喷水网站| 亚洲成色77777| 老司机影院成人| 亚洲av国产av综合av卡| 熟妇人妻不卡中文字幕| 自拍欧美九色日韩亚洲蝌蚪91 | 可以在线观看毛片的网站| 亚洲av电影在线观看一区二区三区 | 91久久精品国产一区二区三区| 特级一级黄色大片| 亚洲精品国产av成人精品| 国内少妇人妻偷人精品xxx网站| 女人被狂操c到高潮| 亚洲怡红院男人天堂| 大又大粗又爽又黄少妇毛片口| 国产精品精品国产色婷婷| 三级经典国产精品| 亚洲欧美中文字幕日韩二区| 91aial.com中文字幕在线观看| 亚洲图色成人| 一边亲一边摸免费视频| 伊人久久精品亚洲午夜| 2021少妇久久久久久久久久久| 午夜日本视频在线| 欧美成人午夜免费资源| 国产精品三级大全| 看免费成人av毛片| 一本一本综合久久| 日韩欧美精品v在线| 国产欧美日韩一区二区三区在线 | 免费大片18禁| 人人妻人人看人人澡| 久久久精品94久久精品| 午夜精品一区二区三区免费看| 99热国产这里只有精品6| 热re99久久精品国产66热6| 亚洲精品成人久久久久久| 熟妇人妻不卡中文字幕| 国产精品久久久久久精品电影| 在线亚洲精品国产二区图片欧美 | 六月丁香七月| 亚洲色图av天堂| 午夜亚洲福利在线播放| 久久久久久久久久久免费av| videossex国产| 插逼视频在线观看| 久久久久九九精品影院| 黄色一级大片看看| 国产高清有码在线观看视频| 天天一区二区日本电影三级| 欧美xxxx黑人xx丫x性爽| 中文天堂在线官网| 国产欧美日韩一区二区三区在线 | 国产一区二区三区综合在线观看 | 少妇的逼水好多| 国产成人a∨麻豆精品| 男的添女的下面高潮视频| 亚洲精品456在线播放app| 国产成年人精品一区二区| 久久久久久久久久成人| 春色校园在线视频观看| 国产黄频视频在线观看| 涩涩av久久男人的天堂| 日韩强制内射视频| 国产男女超爽视频在线观看| 狂野欧美激情性bbbbbb| 在线观看一区二区三区| 成年女人在线观看亚洲视频 | 制服丝袜香蕉在线| 又爽又黄无遮挡网站| 日韩大片免费观看网站| 人人妻人人爽人人添夜夜欢视频| 美国免费a级毛片| 丝袜人妻中文字幕| 十八禁高潮呻吟视频| 欧美精品av麻豆av| 欧美日韩综合久久久久久| 亚洲精品成人av观看孕妇| 成人国产麻豆网| 一区二区日韩欧美中文字幕| 卡戴珊不雅视频在线播放| 人成视频在线观看免费观看| 亚洲国产欧美日韩在线播放| 午夜激情久久久久久久| www.av在线官网国产| 美女国产高潮福利片在线看| 在线观看免费午夜福利视频| 久久婷婷青草| 亚洲成人手机| 日韩av不卡免费在线播放| www.自偷自拍.com| 国产精品三级大全| 午夜福利免费观看在线| 亚洲五月色婷婷综合| 一区二区av电影网| 久久天堂一区二区三区四区| 欧美乱码精品一区二区三区| 日韩大码丰满熟妇| 国产野战对白在线观看| 青青草视频在线视频观看| 黄频高清免费视频| 欧美成人精品欧美一级黄| 一级片'在线观看视频| 国产在线免费精品| 九色亚洲精品在线播放| 日韩制服骚丝袜av| 黄片小视频在线播放| 高清欧美精品videossex| 午夜av观看不卡| 日韩精品免费视频一区二区三区| 国产成人a∨麻豆精品| 日韩大码丰满熟妇| 午夜影院在线不卡| 高清不卡的av网站| 这个男人来自地球电影免费观看 | 亚洲av电影在线观看一区二区三区| 免费观看人在逋| 女人被躁到高潮嗷嗷叫费观| 欧美少妇被猛烈插入视频| 美女脱内裤让男人舔精品视频| 亚洲精品国产区一区二| 亚洲国产中文字幕在线视频| 老司机影院毛片| 国产无遮挡羞羞视频在线观看| 午夜激情久久久久久久| 波野结衣二区三区在线| 亚洲精品在线美女| 国产激情久久老熟女| 国产高清不卡午夜福利| 亚洲成人一二三区av| 欧美中文综合在线视频| 99国产综合亚洲精品| 黄色 视频免费看| 欧美国产精品一级二级三级| 丝袜美腿诱惑在线| 欧美成人精品欧美一级黄| 男女国产视频网站| 欧美日韩一区二区视频在线观看视频在线| 如何舔出高潮| 国产亚洲最大av| 国产免费福利视频在线观看| 九九爱精品视频在线观看| 黑人欧美特级aaaaaa片| 亚洲av欧美aⅴ国产| 精品亚洲成a人片在线观看| 久久99热这里只频精品6学生| 一级,二级,三级黄色视频| 国产一区有黄有色的免费视频| 国产免费视频播放在线视频| 免费人妻精品一区二区三区视频| 日韩电影二区| 一区二区三区乱码不卡18| 欧美日韩一级在线毛片| 国产一卡二卡三卡精品 |