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

    多星座組合RAIM技術(shù)在列車定位中的應(yīng)用

    2022-06-10 03:49:10鐵凱博米根鎖王彥快
    導(dǎo)航定位學(xué)報(bào) 2022年3期
    關(guān)鍵詞:偽距可用性星座

    鐵凱博,米根鎖,王彥快

    (蘭州交通大學(xué), 自動化與電氣工程學(xué)院, 蘭州 730070)

    0 引言

    列車定位精度直接影響列車的運(yùn)行安全,提高列車定位的精度和可靠性已經(jīng)成為不可忽視的重要目標(biāo)[1]。全球衛(wèi)星導(dǎo)航系統(tǒng)(global navigation satellite system, GNSS)在鐵路領(lǐng)域已廣泛應(yīng)用,但在列車實(shí)際運(yùn)行過程中,衛(wèi)星出現(xiàn)故障,會導(dǎo)致定位結(jié)果不準(zhǔn)確[2],這就要求導(dǎo)航系統(tǒng)能對所用衛(wèi)星進(jìn)行實(shí)時(shí)檢測,剔除出現(xiàn)故障的衛(wèi)星。因此,有學(xué)者提出將接收機(jī)自主完好性監(jiān)測(receiver autonomous integrity monitoring,RAIM)技術(shù)[3]應(yīng)用于列車定位。文獻(xiàn)[4]提出了RAIM技術(shù)在列車定位中的可用性預(yù)測方法,并對該技術(shù)在鐵路領(lǐng)域內(nèi)的可用性做出了判斷;文獻(xiàn)[5]將RAIM技術(shù)應(yīng)用于北斗衛(wèi)星導(dǎo)航系統(tǒng)(BeiDou navigation satellite system,BDS)的列車定位中,提出沿股道水平保護(hù)級(along track protect level,ATPL)算法,提高了RAIM技術(shù)在列車定位中的可用性;文獻(xiàn)[6]分析了BDS完好性問題在列車定位中的風(fēng)險(xiǎn),并設(shè)計(jì)了BDS RAIM可用性預(yù)測系統(tǒng),保障了RAIM預(yù)測結(jié)果的可靠性;文獻(xiàn)[7]在故障衛(wèi)星檢測環(huán)節(jié)中引入加權(quán)因子,降低檢測過程中的漏警率,提高了列車定位精度。文獻(xiàn)[8]解決了GNSS性能指標(biāo)與鐵路可用性指標(biāo)的映射問題,為基于GNSS的 RAIM 算法在鐵路領(lǐng)域中的應(yīng)用提供了參考。通過分析文獻(xiàn)可知:目前主要針對單星座RAIM 在列車定位中的應(yīng)用進(jìn)行了研究,單星座RAIM應(yīng)于列車定位時(shí),常會出現(xiàn)因可見星不足或幾何分布較差而降低算法可用性的情況,進(jìn)而影響算法的故障檢測率和故障識別率。由于可見星數(shù)目的限制,單星座 RAIM 在某些高海拔鐵路定位精度并不高[9],甚至不能使用,RAIM算法在鐵路領(lǐng)域內(nèi)的研究還處于起步階段。

    鑒于此,本文提出了將多星座RAIM技術(shù)應(yīng)用于列車定位,通過增加可見星數(shù)目提高RAIM算法的可用性。在構(gòu)造多星座系統(tǒng)觀測方程的基礎(chǔ)上,采用奇偶矢量法對導(dǎo)航衛(wèi)星進(jìn)行故障檢測。通過青藏鐵路列車定位數(shù)據(jù)對全球定位系統(tǒng)(global positioning system, GPS)、伽利略衛(wèi)星導(dǎo)航系統(tǒng)(Galileo satellite navigation system, Galileo)和 BDS多星座組合RAIM和GPS單星座RAIM進(jìn)行對比仿真,分析多星座RAIM在列車定位中的檢測性能。

    1 多星座系統(tǒng)列車定位原理

    1.1 時(shí)空基準(zhǔn)統(tǒng)一

    為構(gòu)成基于列車定位的多星座導(dǎo)航系統(tǒng),需將不同星座的時(shí)間系統(tǒng)進(jìn)行統(tǒng)一,坐標(biāo)系統(tǒng)進(jìn)行轉(zhuǎn)換,采用相同的時(shí)間和坐標(biāo)標(biāo)準(zhǔn),實(shí)現(xiàn)多星座系統(tǒng)的數(shù)據(jù)融合[10-11]。本文將多星座時(shí)間、坐標(biāo)系統(tǒng)統(tǒng)一轉(zhuǎn)換為以GPS為基準(zhǔn)。轉(zhuǎn)換過程如下:

    1)將Galileo系統(tǒng)時(shí)(Galileo system time, GST)轉(zhuǎn)化為GPS時(shí)(GPS time, GPST)的計(jì)算公式為

    式中:TGPST為 GPS時(shí);TGST為 Galileo系統(tǒng)時(shí);TGGTO為GPST與GST間的偏差;A0G、A1G、t0G和SWNOG分別為時(shí)間偏差的常數(shù)項(xiàng)、變化率、參考時(shí)間和參考星期個(gè)數(shù);TTOW為星期內(nèi)的時(shí)間;T是常數(shù),為604 800;SWN為Galileo的星期個(gè)數(shù);mod為求余函數(shù)。

    2)將北斗時(shí)(BDS time, BDT)轉(zhuǎn)化為GPST的計(jì)算公式為

    式中,TBDT為北斗系統(tǒng)時(shí)。

    3)將格洛納斯衛(wèi)星導(dǎo)航系統(tǒng)(global navigation satellite system, GLONASS)時(shí)(GLONASS time,GLONASST)轉(zhuǎn)化為GPST的計(jì)算公式為

    式中:TGLONASST為GLONASS時(shí);n為協(xié)調(diào)世界時(shí)(coordinated universal time, UTC)與國際原子時(shí)(international atomic time, TAI)之間不斷調(diào)整的參數(shù),由國際地球自轉(zhuǎn)服務(wù)組織(International Earth Rotation Service, IERS)提供。

    GPS采用1984世界大地坐標(biāo)系(world geodetic coordinate system 1984, WGS84);Galileo采用伽利略地球參考框架(Galileo terrestrial reference frame,GTRF);BDS采用北斗坐標(biāo)系(BeiDou coordinate system,BDCS);GLONASS采用俄羅斯大地坐標(biāo)框架(Parametry Zelmy1990, PZ90)。由于列車偽距定位的精度為米級,BDCS坐標(biāo)和GTRF坐標(biāo)轉(zhuǎn)換為WGS84坐標(biāo)產(chǎn)生的誤差可以忽略不計(jì)[12],只需將PZ90坐標(biāo)轉(zhuǎn)換為WGS84坐標(biāo)即可,采用布爾薩-沃爾夫(Bursa-Wolf)模型[13]進(jìn)行PZ90和WGS84間的坐標(biāo)轉(zhuǎn)換,Bursa-Wolf模型亦稱為黑爾默特(Helmert)7參數(shù)模型,其計(jì)算公式為

    1.2 列車定位原理

    在列車運(yùn)行過程中,可見星連續(xù)不斷地向車載接收機(jī)發(fā)送星歷參數(shù)和時(shí)間信息,車載計(jì)算機(jī)計(jì)算出與N顆可見星之間的偽距(ρ1,ρ2,ρ3,…,ρN),通過聯(lián)立方程進(jìn)行列車位置的定位解算,即

    式中:ρi(i=1,2,...,N)為第i號衛(wèi)星到車載接收機(jī)的偽距,N為可見星數(shù)目;(xi,yi,zi)為第i號衛(wèi)星的空間位置;(xu,yu,zu)為列車位置;c為光速;tu為車載接收機(jī)與衛(wèi)星間的鐘差,其定位原理如圖1所示。

    圖1 列車偽距定位原理

    若參與解算的可見星發(fā)生故障,會造成ρi的測量值出現(xiàn)偏差,導(dǎo)致列車定位解算結(jié)果不準(zhǔn)確。因此,需在列車定位解算之前,采用RAIM技術(shù)對衛(wèi)星的健康狀況進(jìn)行判斷,剔除出現(xiàn)故障的衛(wèi)星。

    2 基于列車定位的多星座RAIM算法

    在列車定位中,RAIM算法采用冗余技術(shù)對車載接收機(jī)接收到的導(dǎo)航信息進(jìn)行一致性校驗(yàn),其本質(zhì)是通過增加可見星數(shù)目來增加信息冗余,實(shí)現(xiàn)對衛(wèi)星的故障檢測[14]。當(dāng)列車導(dǎo)航系統(tǒng)至少監(jiān)測到4顆衛(wèi)星時(shí),可進(jìn)行用戶位置的定位解算;至少監(jiān)測到5顆衛(wèi)星時(shí),可實(shí)現(xiàn)故障衛(wèi)星的檢測;至少監(jiān)測到6顆衛(wèi)星時(shí),可對故障衛(wèi)星進(jìn)行識別。因此,可見星數(shù)目對 RAIM 算法的可用性有很大影響,通過采用多星座系統(tǒng)增加可見星數(shù)目進(jìn)而改善幾何分布情況對 RAIM 算法的可用性會有相應(yīng)提升[15]。

    在鐵路沿線中,海拔較高時(shí),會使部分衛(wèi)星的檢測效果變差[16],影響算法的可用性。多星座組合系統(tǒng)通過增加線路中可見衛(wèi)星數(shù)目增加冗余觀測信息,進(jìn)而提高 RAIM 算法在列車定位中的可用性。當(dāng)車載接收機(jī)監(jiān)測到的可見衛(wèi)星數(shù)N>n時(shí)(n=4+m,m為星座數(shù)),可進(jìn)行故障衛(wèi)星的檢測與識別。對于單星座系統(tǒng)至少監(jiān)測到n=4+1=5顆衛(wèi)星;對于雙星座系統(tǒng)至少監(jiān)測到n=4+2=6顆衛(wèi)星;對于三星座系統(tǒng)至少監(jiān)測到n=4+3=7顆衛(wèi)星;以此類推。圖2為多星座 RAIM算法的程序流程圖。圖2中n為故障檢測所需的最少可見星數(shù)。

    圖2 多星座RAIM算法程序流程

    2.1 構(gòu)造多星座系統(tǒng)觀測方程

    假設(shè)在某一時(shí)刻共有N顆可見星,對式(5)進(jìn)行泰勒展開,并考慮觀測誤差,線性化后的多星座系統(tǒng)觀測方程為

    式中:y?RN×1為觀測偽距與計(jì)算偽距的差值;x?Rl×1為列車的三維位置信息和m個(gè)星座對應(yīng)的時(shí)鐘誤差參數(shù),l=m+3;ε?RN×1為觀測偽距噪聲矢量,服從均值為 0,方差為的正態(tài)分布。H?RN?l為N ?l維的觀測矩陣,其表達(dá)式為

    式中:A為列車與衛(wèi)星之間的方向余弦值;B為對應(yīng)的各星座車載接收機(jī)鐘差;Aij(i=1,…,m,j=1,…,3)為對應(yīng)第i個(gè)星座觀測矩陣Ai的第j列元素,Ai可表示為

    2.2 基于奇偶矢量法的故障衛(wèi)星檢測識別

    假設(shè)偽距誤差ε=O,則系統(tǒng)的線性化觀測方程為

    將系數(shù)矩陣H分解成

    式中:Q是N?N階的正交矩陣;R為N?l階的上三角矩陣;O為零矩陣。

    將式(10)代入式(9)并左乘QT,可得

    式中,QT和R可分別表示為

    式中:Qx為QT的前l(fā)行,Q為QT的后N-l行;Rx為R的前l(fā)行,O是R的后N-l行,為零矩陣。

    將式(12)代入式(11)中可得

    可得最小二乘解

    若考慮到觀測誤差ε,即y=Hx+ε ,則有

    式中:p為奇偶空間上的奇偶矢量;Q為奇偶空間矩陣;奇偶空間矢量p是由ε在奇偶空間矩陣Q上投影所得,可反映故障衛(wèi)星偏差量等信息。因此,可基于奇偶矢量p構(gòu)造檢測統(tǒng)計(jì)量,進(jìn)行故障衛(wèi)星的檢測識別。

    衛(wèi)星偽距殘差向量平方和為v=pTp。根據(jù)統(tǒng)計(jì)分布理論,當(dāng)偽距噪聲矢量ε中的各個(gè)分量相互獨(dú)立,且服從均值為 0,方差為的正態(tài)分布時(shí),v/服從自由度為N-l的中心X2分布;均值不為0時(shí),v/服從自由度為N?l的非中心χ2分布,非中心參數(shù)為λ,即

    構(gòu)造檢驗(yàn)統(tǒng)計(jì)量ri為

    式中:ri為第i顆衛(wèi)星的檢驗(yàn)統(tǒng)計(jì)量;Qi為第i顆衛(wèi)星的奇偶空間矩陣。

    已知誤警率PFA,則有

    式中:P(ri>Tr)為檢測統(tǒng)計(jì)量大于檢測門限的概率,為標(biāo)準(zhǔn)正態(tài)分布的概率密度函數(shù)。

    根據(jù)式(18)可求出檢測門限Tr,將每個(gè)檢測統(tǒng)計(jì)量ri與Tr進(jìn)行比較,若ri>Tr,說明第i顆衛(wèi)星出現(xiàn)故障。

    2.3 多星座RAIM算法可用性判斷

    為保障列車定位結(jié)果的可靠性,在可見星數(shù)目達(dá)到要求后,還需對沿線可見星的空間幾何分布進(jìn)行預(yù)測,對RAIM算法的可用性進(jìn)行判斷。

    本文采用水平保護(hù)級(horizontal protect level,HPL)THPL來判斷RAIM可用性。根據(jù)車載接收機(jī)接收到的列車定位數(shù)據(jù)計(jì)算THPL值,當(dāng)THPL小于水平告警閾值(horizontal alert limit,HAL)THAL時(shí),算法可用,反之則無效,車載需重新接收衛(wèi)星信號進(jìn)行可用性判斷。

    利用最大特征斜率Lmax乘以臨界偏差Pbias計(jì)算出THPL的值。首先,計(jì)算每顆衛(wèi)星的特征斜率為

    式中:Li為第i顆衛(wèi)星的特征斜率;G1i、G2i分別為矩陣G中元素,G= (HTH)?1HT;S為矩陣ii S的對角元素,S=I? (HTH)?1HT。

    然后根據(jù)檢測門限Tr,基于漏檢概率PMD計(jì)算出臨界偏差Pbias=σ0., λ的計(jì)算方法為

    水平保護(hù)限值可由式(21)計(jì)算出,即

    式中,SHPL為計(jì)算出的水平保護(hù)限值。

    將SHPL與HAL進(jìn)行比較,若SHPL

    3 多星座RAIM算法仿真分析

    利用矩陣實(shí)驗(yàn)室(matrix laboratory, MATLAB)軟件對算法進(jìn)行仿真分析,本次實(shí)驗(yàn)選用 GPS、Galileo和BDS三星座系統(tǒng)。首先根據(jù)3個(gè)星座系統(tǒng)的衛(wèi)星及軌道信息參數(shù),對星座進(jìn)行模擬,模擬中選擇24顆GPS衛(wèi)星,24顆Galileo衛(wèi)星和14顆BDS衛(wèi)星,選擇列車用戶位置在青藏鐵路某段線路的固定位置處,用戶觀測衛(wèi)星的高度截至角為15°。

    3.1 多星座RAIM算法可用性驗(yàn)證

    首先對多星座系統(tǒng)和單星座系統(tǒng)下的可見星數(shù)量進(jìn)行仿真對比。圖3為GPS/Galileo/BDS多星座系統(tǒng)與GPS單星座系統(tǒng)在某天同時(shí)刻下可見星數(shù)量的對比情況。由圖3可知,在24 h內(nèi)的任意時(shí)刻,多星座系統(tǒng)的可見星數(shù)量均高于單星座系統(tǒng),由于多星座系統(tǒng)增加了衛(wèi)星數(shù)量,所以在相同路況下,可見星的數(shù)量也隨之增多。隨著可見星數(shù)量增多,更容易滿足 RAIM 算法對冗余觀測信息的需求,在列車運(yùn)行中的某些特殊路段,如城市、山川,多星座 RAIM 能夠更好地保障導(dǎo)航信息的連續(xù)性和完好性。

    圖3 可見星數(shù)量對比

    圖4為 GPS/Galileo/BDS多星座系統(tǒng)和 GPS單星座系統(tǒng)在 24 h內(nèi)對應(yīng)的位置精度因子(position dilution of precision, PDOP)曲線。

    圖4 位置精度因子曲線對比

    由圖4可知,多星座系統(tǒng)各時(shí)刻的位置精度因子均小于GPS單星座系統(tǒng),說明多星座系統(tǒng)能夠有效改善可見星幾何分布的情況,進(jìn)而提高算法的可用性。

    圖5表示青藏鐵路某段中密度線路在GPS/Galileo/BDS多星座系統(tǒng)和 GPS單星座系統(tǒng)下的SHPL值對比情況。

    圖5 SHPL對比

    根據(jù)鐵路領(lǐng)域完好性水平要求,中密度線路的水平告警閾值為20 m[17]。由圖5可知,GPS單星座系統(tǒng)SHPL值的波動范圍為11.2~27.0 m,多星座系統(tǒng)SHPL值的波動范圍為6.0~18.7 m,在34 min之后,GPS單星座系統(tǒng)的SHPL開始出現(xiàn)大于水平告警閾值的情況,發(fā)生虛警現(xiàn)象。多星座系統(tǒng)的SHPL值始終小于水平告警閾值,未出現(xiàn)虛警,且在同一時(shí)刻,多星座系統(tǒng)的SHPL值均小于 GPS單星座系統(tǒng),說明列車在相同運(yùn)行路況下,多星座組合RAIM算法的可用性比GPS單星座系統(tǒng)更好。

    3.2 多星座RAIM算法故障檢測識別

    本次實(shí)驗(yàn)所用數(shù)據(jù)為無故障衛(wèi)星時(shí)的列車定位數(shù)據(jù),需通過人為給衛(wèi)星注入偽距偏差來模擬衛(wèi)星出現(xiàn)故障,進(jìn)而對算法進(jìn)行驗(yàn)證。選用線路數(shù)據(jù)中幾何構(gòu)型較好的21顆無故障可見星,依次編號為1號至21號,其中:1號至8號為GPS衛(wèi)星;9號至16號為Galileo衛(wèi)星;17號至21號為BDS衛(wèi)星。設(shè)定誤警率P= 1 × 1 0?6,漏警率P= 1 × 1 0?6,F(xiàn)AMD假設(shè)各導(dǎo)航系統(tǒng)的偽距觀測噪聲水平相同,取σ0=5 m 。

    在第80—120秒的時(shí)間間隔內(nèi),給9號衛(wèi)星注入60 m的持續(xù)故障,其結(jié)果如圖6所示。

    圖6 持續(xù)故障檢測

    由圖6可知,在第80—120秒的時(shí)間間隔內(nèi),檢測統(tǒng)計(jì)量超出了預(yù)設(shè)閾值,算法檢測出了衛(wèi)星出現(xiàn)持續(xù)偽距偏差。

    在第20—80秒的時(shí)間間隔內(nèi),給第17號衛(wèi)星注入0~60 m的慢變斜坡故障,其結(jié)果如圖7所示。

    圖7 慢變斜坡故障檢測

    由圖7可知,在第32秒時(shí),檢測統(tǒng)計(jì)量超出閾值,檢測到有故障發(fā)生。隨著偽距偏差的增大,檢測統(tǒng)計(jì)量也不斷增大,算法檢測出了衛(wèi)星出現(xiàn)慢變斜坡偽距偏差。

    通過故障識別算法得到兩種故障情況下各衛(wèi)星的檢測統(tǒng)計(jì)量如表1所示,檢測統(tǒng)計(jì)量最大的衛(wèi)星即為故障星。

    表1 各星座系統(tǒng)衛(wèi)星檢測統(tǒng)計(jì)量的值

    由表1可知,在持續(xù)故障情況下,9號衛(wèi)星的檢測統(tǒng)計(jì)量最大,為算法識別出的故障星。在慢變斜坡故障情況下,17號衛(wèi)星的檢測統(tǒng)計(jì)量最大,為算法識別出的故障星。算法有效識別出了在持續(xù)故障和慢變斜坡故障兩種情況下預(yù)設(shè)的故障星。

    3.3 多星座RAIM算法檢測性能分析

    為檢驗(yàn)多星座 RAIM 算法的檢測性能,利用3.2節(jié)列車定位數(shù)據(jù)對算法的故障檢測率和故障識別率進(jìn)行仿真分析。

    給某顆衛(wèi)星注入 0~200 m的偽距偏差,采樣時(shí)間為24 h,采樣間隔1 min,仿真結(jié)果如圖8和圖9所示。

    圖8 故障檢測率對比

    圖9 故障識別率對比

    由仿真結(jié)果可知,當(dāng)偽距偏差達(dá)到100 m時(shí),多星座系統(tǒng)的故障檢測率和故障識別率均可達(dá) 100%,而對于GPS單星座系統(tǒng)當(dāng)偽距偏差達(dá)到136 m時(shí),故障檢測率才達(dá)到100%,偽距偏差達(dá)到160 m時(shí),故障識別率才達(dá)到94.1%。表明多星座組合RAIM的故障檢測性能優(yōu)于GPS單星座RAIM。

    4 結(jié)束語

    本文利用多星座組合RAIM技術(shù),實(shí)時(shí)檢測接收機(jī)自身以及衛(wèi)星故障,進(jìn)而提高列車運(yùn)行過程中定位的精度。在列車定位過程中,采用多星座系統(tǒng)可在相同路況下獲得更多的可見星,從而提高RAIM算法的可用性及信息冗余度。在故障檢測部分,利用奇偶矢量法計(jì)算量少、檢測速度快的優(yōu)勢,可實(shí)現(xiàn)故障衛(wèi)星的快速檢測。結(jié)合青藏鐵路列車定位數(shù)據(jù)分別從算法可用性、故障檢測識別和檢測性能三個(gè)方面進(jìn)行了仿真分析。結(jié)果表明,多星座 RAIM 算法的可用性較單星座系統(tǒng)更好,且能迅速地檢測出定位過程中出現(xiàn)的故障衛(wèi)星,并予以識別。由于信息冗余度的增加,多星座系統(tǒng)的故障檢測率和故障識別率較單星座系統(tǒng)均有所提高,位置精度因子較單星座系統(tǒng)有明顯的降低。多星座RAIM技術(shù)可有效提高列車定位精度,保證鐵路運(yùn)輸行車安全。

    猜你喜歡
    偽距可用性星座
    基于文獻(xiàn)計(jì)量學(xué)的界面設(shè)計(jì)可用性中外對比研究
    包裝工程(2023年24期)2023-12-27 09:18:26
    基于輻射傳輸模型的GOCI晨昏時(shí)段數(shù)據(jù)的可用性分析
    北斗偽距觀測值精度分析
    星座
    健康女性(2017年3期)2017-04-27 21:34:31
    12星座之我愛洗澡
    GNSS偽距粗差的開窗探測及修復(fù)
    星座
    都市麗人(2016年12期)2016-12-14 09:54:15
    星座
    都市麗人(2016年7期)2016-08-03 10:06:05
    聯(lián)合碼偽距和載波寬巷組合的相對定位技術(shù)研究
    空客A320模擬機(jī)FD1+2可用性的討論
    河南科技(2015年7期)2015-03-11 16:23:13
    久久久久免费精品人妻一区二区| 国产三级在线视频| 久久中文看片网| 性色av乱码一区二区三区2| 99精品久久久久人妻精品| 老司机福利观看| 小蜜桃在线观看免费完整版高清| 国产高潮美女av| 亚洲激情在线av| 少妇裸体淫交视频免费看高清| 一进一出抽搐动态| 国内精品久久久久精免费| 亚洲人成网站在线播放欧美日韩| 国产精品av视频在线免费观看| 亚洲性夜色夜夜综合| 国产精品美女特级片免费视频播放器| 麻豆一二三区av精品| 日本黄色视频三级网站网址| 99精品欧美一区二区三区四区| 精品免费久久久久久久清纯| 又爽又黄无遮挡网站| 在线播放无遮挡| 身体一侧抽搐| 国产真实伦视频高清在线观看 | 一个人看视频在线观看www免费 | 日韩 欧美 亚洲 中文字幕| 亚洲中文字幕日韩| av福利片在线观看| 在线观看免费视频日本深夜| 亚洲精品一卡2卡三卡4卡5卡| 内射极品少妇av片p| 国产精品乱码一区二三区的特点| 亚洲va日本ⅴa欧美va伊人久久| 真实男女啪啪啪动态图| 亚洲avbb在线观看| 天堂网av新在线| 精品久久久久久成人av| 欧美av亚洲av综合av国产av| 动漫黄色视频在线观看| 亚洲中文字幕日韩| 久久国产精品影院| 两个人的视频大全免费| 成人特级黄色片久久久久久久| 国产精品久久久久久久电影 | 波多野结衣高清作品| 国产精品一及| 村上凉子中文字幕在线| 国产精品免费一区二区三区在线| 亚洲av免费高清在线观看| 亚洲天堂国产精品一区在线| 日韩欧美 国产精品| 99在线人妻在线中文字幕| 欧美成狂野欧美在线观看| 亚洲在线自拍视频| 欧美+亚洲+日韩+国产| 亚洲精品国产精品久久久不卡| 亚洲欧美日韩高清专用| 中文字幕高清在线视频| 村上凉子中文字幕在线| 嫩草影院入口| 精品久久久久久久久久免费视频| 日韩免费av在线播放| 男人舔女人下体高潮全视频| 一区福利在线观看| 91麻豆av在线| 99久久精品国产亚洲精品| 国产欧美日韩一区二区精品| 欧美日韩福利视频一区二区| 久久久色成人| 老司机午夜十八禁免费视频| 久久午夜亚洲精品久久| 久久久久九九精品影院| 男人和女人高潮做爰伦理| 国产美女午夜福利| 性欧美人与动物交配| 欧美中文日本在线观看视频| 成人鲁丝片一二三区免费| bbb黄色大片| 日韩中文字幕欧美一区二区| 国产一区在线观看成人免费| 在线观看午夜福利视频| 精品久久久久久久久久免费视频| 国内少妇人妻偷人精品xxx网站| 女生性感内裤真人,穿戴方法视频| tocl精华| 国产成+人综合+亚洲专区| 无遮挡黄片免费观看| 久久6这里有精品| 一卡2卡三卡四卡精品乱码亚洲| 白带黄色成豆腐渣| 国产精品,欧美在线| 久久性视频一级片| 18+在线观看网站| 嫩草影视91久久| 亚洲精品国产精品久久久不卡| 少妇人妻精品综合一区二区 | 久久久精品大字幕| 久久精品91蜜桃| 天堂√8在线中文| 18禁黄网站禁片午夜丰满| 午夜免费男女啪啪视频观看 | 欧美另类亚洲清纯唯美| 国产极品精品免费视频能看的| 欧美黑人欧美精品刺激| 欧美日韩福利视频一区二区| 欧美绝顶高潮抽搐喷水| 日本 av在线| 日韩高清综合在线| 色综合婷婷激情| 免费看a级黄色片| 国产黄色小视频在线观看| 校园春色视频在线观看| 免费在线观看影片大全网站| 午夜亚洲福利在线播放| 国产成人av激情在线播放| 亚洲专区中文字幕在线| 久久亚洲真实| 成人精品一区二区免费| 国内揄拍国产精品人妻在线| 精品欧美国产一区二区三| 看片在线看免费视频| 午夜影院日韩av| 日韩av在线大香蕉| 美女被艹到高潮喷水动态| 人人妻人人澡欧美一区二区| 最近最新免费中文字幕在线| 精品久久久久久久人妻蜜臀av| 国产三级在线视频| 日本一二三区视频观看| 中文资源天堂在线| 精品国产美女av久久久久小说| tocl精华| 又黄又爽又免费观看的视频| 国产av不卡久久| 日本撒尿小便嘘嘘汇集6| 老司机深夜福利视频在线观看| 亚洲片人在线观看| 色精品久久人妻99蜜桃| 国产v大片淫在线免费观看| 亚洲精品色激情综合| 色吧在线观看| 久久久久久久午夜电影| 国产精品女同一区二区软件 | 中文亚洲av片在线观看爽| 成人国产一区最新在线观看| 少妇丰满av| 一二三四社区在线视频社区8| 女人十人毛片免费观看3o分钟| 久久久久国内视频| 国产69精品久久久久777片| 日本一本二区三区精品| 老熟妇乱子伦视频在线观看| 亚洲成人精品中文字幕电影| 国产午夜福利久久久久久| 亚洲在线自拍视频| 欧美日韩国产亚洲二区| 国产乱人伦免费视频| 亚洲中文字幕日韩| 欧美又色又爽又黄视频| 久久婷婷人人爽人人干人人爱| 久久精品综合一区二区三区| 免费看美女性在线毛片视频| 午夜精品在线福利| 国产v大片淫在线免费观看| 午夜精品一区二区三区免费看| 国产av麻豆久久久久久久| 国产精品亚洲一级av第二区| 88av欧美| 99久久综合精品五月天人人| 国产高清视频在线播放一区| 亚洲欧美一区二区三区黑人| 日日摸夜夜添夜夜添小说| 午夜福利18| 在线观看午夜福利视频| 一个人看视频在线观看www免费 | 亚洲av成人精品一区久久| 国产三级中文精品| 国产成人影院久久av| 看片在线看免费视频| 欧美bdsm另类| 国产美女午夜福利| 又紧又爽又黄一区二区| 噜噜噜噜噜久久久久久91| 欧洲精品卡2卡3卡4卡5卡区| 国产三级在线视频| 日韩欧美精品免费久久 | 日韩免费av在线播放| 狠狠狠狠99中文字幕| 国产 一区 欧美 日韩| 国产三级中文精品| 日本一二三区视频观看| 少妇熟女aⅴ在线视频| 亚洲av免费在线观看| 老司机午夜十八禁免费视频| 精品不卡国产一区二区三区| 舔av片在线| 亚洲第一电影网av| 亚洲成人精品中文字幕电影| 一级作爱视频免费观看| 又紧又爽又黄一区二区| av在线天堂中文字幕| 日韩欧美国产在线观看| 国产成人系列免费观看| 热99re8久久精品国产| 亚洲av二区三区四区| 国产av不卡久久| 特大巨黑吊av在线直播| 亚洲精品在线观看二区| 欧美乱妇无乱码| 国产探花极品一区二区| 精品国产超薄肉色丝袜足j| 长腿黑丝高跟| 尤物成人国产欧美一区二区三区| 伊人久久大香线蕉亚洲五| 99久久无色码亚洲精品果冻| 国产成+人综合+亚洲专区| 97超视频在线观看视频| 一卡2卡三卡四卡精品乱码亚洲| 国产精品综合久久久久久久免费| 18美女黄网站色大片免费观看| 成人精品一区二区免费| 亚洲,欧美精品.| 欧美黄色淫秽网站| 国产爱豆传媒在线观看| 欧美成人一区二区免费高清观看| 精品久久久久久久久久久久久| 亚洲一区二区三区色噜噜| 久久久久久久午夜电影| 亚洲欧美日韩高清在线视频| 日韩欧美一区二区三区在线观看| 特大巨黑吊av在线直播| 成人国产一区最新在线观看| 午夜精品在线福利| 手机成人av网站| 亚洲国产欧美人成| 日本一本二区三区精品| 精华霜和精华液先用哪个| 国产欧美日韩一区二区精品| 在线天堂最新版资源| 欧美性感艳星| 午夜免费男女啪啪视频观看 | 精品久久久久久久人妻蜜臀av| 99久久精品热视频| 色视频www国产| 成年女人看的毛片在线观看| 国模一区二区三区四区视频| 少妇熟女aⅴ在线视频| 亚洲熟妇熟女久久| 亚洲人成网站高清观看| 久久这里只有精品中国| 99久久精品国产亚洲精品| 亚洲 欧美 日韩 在线 免费| 精品久久久久久久末码| 久久精品亚洲精品国产色婷小说| 国产一区二区在线av高清观看| 两人在一起打扑克的视频| 夜夜看夜夜爽夜夜摸| 香蕉丝袜av| 日韩大尺度精品在线看网址| 久久婷婷人人爽人人干人人爱| 桃色一区二区三区在线观看| 亚洲 国产 在线| 在线观看午夜福利视频| 日本熟妇午夜| 听说在线观看完整版免费高清| 国产综合懂色| 欧美bdsm另类| 欧美一级毛片孕妇| 亚洲欧美日韩无卡精品| aaaaa片日本免费| 不卡一级毛片| 一边摸一边抽搐一进一小说| 少妇高潮的动态图| а√天堂www在线а√下载| 国产激情偷乱视频一区二区| 欧美黄色片欧美黄色片| 尤物成人国产欧美一区二区三区| 国产一区在线观看成人免费| 欧美色欧美亚洲另类二区| av在线蜜桃| 亚洲精品影视一区二区三区av| 久久精品国产综合久久久| av欧美777| 免费人成视频x8x8入口观看| 在线观看一区二区三区| 久久天躁狠狠躁夜夜2o2o| 窝窝影院91人妻| 亚洲真实伦在线观看| 制服丝袜大香蕉在线| 久久草成人影院| 激情在线观看视频在线高清| 婷婷六月久久综合丁香| 母亲3免费完整高清在线观看| 99热这里只有精品一区| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久久人人做人人爽| 欧美成人一区二区免费高清观看| 亚洲天堂国产精品一区在线| 亚洲国产高清在线一区二区三| 国产欧美日韩一区二区三| 国产亚洲av嫩草精品影院| 亚洲av成人精品一区久久| 12—13女人毛片做爰片一| 听说在线观看完整版免费高清| 午夜福利成人在线免费观看| 亚洲国产精品久久男人天堂| 欧美色欧美亚洲另类二区| 看免费av毛片| 波野结衣二区三区在线 | 国产精品女同一区二区软件 | 久久久成人免费电影| 最近视频中文字幕2019在线8| 亚洲国产色片| 日韩欧美在线乱码| 最新中文字幕久久久久| 免费观看的影片在线观看| 婷婷亚洲欧美| 夜夜爽天天搞| 国产精品98久久久久久宅男小说| 丝袜美腿在线中文| 亚洲色图av天堂| 性色av乱码一区二区三区2| 欧美高清成人免费视频www| 亚洲av成人av| 淫妇啪啪啪对白视频| 久久99热这里只有精品18| 欧美日韩亚洲国产一区二区在线观看| 男女下面进入的视频免费午夜| 国产亚洲欧美98| 欧美+亚洲+日韩+国产| 99精品久久久久人妻精品| 不卡一级毛片| 精品欧美国产一区二区三| 夜夜看夜夜爽夜夜摸| 久久久久国内视频| 日韩中文字幕欧美一区二区| 久久九九热精品免费| 中文亚洲av片在线观看爽| 午夜福利在线观看吧| 色老头精品视频在线观看| 超碰av人人做人人爽久久 | 一二三四社区在线视频社区8| 国产色爽女视频免费观看| 免费av毛片视频| 亚洲久久久久久中文字幕| 欧美av亚洲av综合av国产av| 特大巨黑吊av在线直播| 国内揄拍国产精品人妻在线| 欧美国产日韩亚洲一区| 18禁国产床啪视频网站| 级片在线观看| 搡老岳熟女国产| 美女高潮的动态| 男女午夜视频在线观看| 久久精品国产99精品国产亚洲性色| 大型黄色视频在线免费观看| 久久久久性生活片| 女人高潮潮喷娇喘18禁视频| av片东京热男人的天堂| 亚洲精品日韩av片在线观看 | 久久欧美精品欧美久久欧美| 国内精品久久久久久久电影| 亚洲乱码一区二区免费版| 97超级碰碰碰精品色视频在线观看| 神马国产精品三级电影在线观看| 国内精品久久久久久久电影| 国产午夜精品久久久久久一区二区三区 | 精品无人区乱码1区二区| 91字幕亚洲| 热99re8久久精品国产| 又粗又爽又猛毛片免费看| 成人性生交大片免费视频hd| 在线观看av片永久免费下载| 亚洲人成伊人成综合网2020| 亚洲自拍偷在线| 最新在线观看一区二区三区| 3wmmmm亚洲av在线观看| 99久国产av精品| 美女大奶头视频| 久久精品国产99精品国产亚洲性色| 在线视频色国产色| 亚洲国产中文字幕在线视频| 99久久精品热视频| 黑人欧美特级aaaaaa片| 天堂动漫精品| 99久久精品一区二区三区| 亚洲精品国产精品久久久不卡| 国产激情偷乱视频一区二区| 99国产综合亚洲精品| 一区二区三区国产精品乱码| 变态另类丝袜制服| 免费av观看视频| or卡值多少钱| 亚洲最大成人中文| 天天躁日日操中文字幕| 少妇高潮的动态图| 18美女黄网站色大片免费观看| 久久久久国产精品人妻aⅴ院| av视频在线观看入口| 亚洲色图av天堂| 欧美成人一区二区免费高清观看| 日本黄大片高清| 国产淫片久久久久久久久 | 一级a爱片免费观看的视频| 免费看日本二区| 国产成人系列免费观看| 日本三级黄在线观看| 他把我摸到了高潮在线观看| xxx96com| 久久久久久久午夜电影| 无限看片的www在线观看| 我要搜黄色片| 一个人看的www免费观看视频| 亚洲精品影视一区二区三区av| 99视频精品全部免费 在线| 国产亚洲精品久久久com| 国产免费一级a男人的天堂| 国产成人aa在线观看| 99国产极品粉嫩在线观看| 桃色一区二区三区在线观看| 亚洲欧美日韩无卡精品| 欧美国产日韩亚洲一区| www.999成人在线观看| 国产激情欧美一区二区| 深爱激情五月婷婷| h日本视频在线播放| 午夜老司机福利剧场| 久久精品91蜜桃| 欧美日韩一级在线毛片| 嫩草影视91久久| 夜夜夜夜夜久久久久| 成人三级黄色视频| 亚洲国产精品久久男人天堂| 日本精品一区二区三区蜜桃| 窝窝影院91人妻| 最近视频中文字幕2019在线8| 久久精品国产亚洲av涩爱 | 高清在线国产一区| 99在线人妻在线中文字幕| 日韩欧美在线乱码| 18禁美女被吸乳视频| 亚洲欧美日韩东京热| 国内精品久久久久精免费| 免费高清视频大片| 亚洲国产日韩欧美精品在线观看 | 99视频精品全部免费 在线| 最近最新免费中文字幕在线| 我的老师免费观看完整版| www.www免费av| 变态另类成人亚洲欧美熟女| 久久久色成人| 国产三级在线视频| 色尼玛亚洲综合影院| 97超视频在线观看视频| 18+在线观看网站| 亚洲欧美一区二区三区黑人| 亚洲成人免费电影在线观看| 国产精品美女特级片免费视频播放器| 亚洲国产精品999在线| 国产69精品久久久久777片| 神马国产精品三级电影在线观看| 好男人在线观看高清免费视频| 黄色成人免费大全| 国内毛片毛片毛片毛片毛片| 精品99又大又爽又粗少妇毛片 | 老司机午夜十八禁免费视频| 国内精品美女久久久久久| 午夜亚洲福利在线播放| 人人妻人人看人人澡| 国产国拍精品亚洲av在线观看 | av国产免费在线观看| av片东京热男人的天堂| 男人舔奶头视频| 在线观看日韩欧美| 国产免费一级a男人的天堂| 亚洲最大成人手机在线| 夜夜看夜夜爽夜夜摸| 十八禁人妻一区二区| 黄色日韩在线| 日韩av在线大香蕉| 国产一区二区三区视频了| 老司机午夜十八禁免费视频| 搡女人真爽免费视频火全软件 | 精品国产三级普通话版| 精品久久久久久久久久免费视频| av在线天堂中文字幕| 国产精品99久久99久久久不卡| 1024手机看黄色片| 在线视频色国产色| 黑人欧美特级aaaaaa片| av在线天堂中文字幕| 婷婷精品国产亚洲av| 久久久色成人| 成人三级黄色视频| 亚洲人成网站高清观看| 香蕉丝袜av| 欧美区成人在线视频| 狂野欧美白嫩少妇大欣赏| 91av网一区二区| 真实男女啪啪啪动态图| av在线天堂中文字幕| 亚洲国产精品合色在线| 国产精品 国内视频| 黄色成人免费大全| 高清毛片免费观看视频网站| www国产在线视频色| 久久精品91蜜桃| 俺也久久电影网| 国产午夜福利久久久久久| 18禁裸乳无遮挡免费网站照片| 成年女人看的毛片在线观看| 天美传媒精品一区二区| 色老头精品视频在线观看| 免费看十八禁软件| 亚洲中文字幕一区二区三区有码在线看| 久久久色成人| 色哟哟哟哟哟哟| tocl精华| 日本成人三级电影网站| 精品久久久久久久末码| 国产亚洲精品久久久久久毛片| 俺也久久电影网| 国产99白浆流出| 亚洲人成电影免费在线| 大型黄色视频在线免费观看| 一级黄片播放器| 1000部很黄的大片| 精品熟女少妇八av免费久了| 熟女电影av网| 黄片大片在线免费观看| 变态另类成人亚洲欧美熟女| 中亚洲国语对白在线视频| 一个人看视频在线观看www免费 | 午夜激情福利司机影院| 亚洲一区二区三区色噜噜| 久久久久九九精品影院| 变态另类成人亚洲欧美熟女| 欧美在线一区亚洲| 亚洲欧美日韩东京热| 少妇裸体淫交视频免费看高清| 好看av亚洲va欧美ⅴa在| 天天添夜夜摸| 午夜免费男女啪啪视频观看 | 麻豆成人av在线观看| 日本一二三区视频观看| 麻豆久久精品国产亚洲av| 亚洲成人久久爱视频| 国产av在哪里看| 国产精品久久久久久久久免 | 久久久久久九九精品二区国产| 18禁黄网站禁片午夜丰满| 精品人妻1区二区| 国产97色在线日韩免费| 国产一区在线观看成人免费| 国产久久久一区二区三区| 精品乱码久久久久久99久播| 啦啦啦观看免费观看视频高清| 国产精品三级大全| 国产伦精品一区二区三区四那| 长腿黑丝高跟| 可以在线观看毛片的网站| 亚洲精品国产精品久久久不卡| 9191精品国产免费久久| 黄片大片在线免费观看| 成人三级黄色视频| 99精品在免费线老司机午夜| 欧美黑人巨大hd| 男人舔奶头视频| 国产综合懂色| 51午夜福利影视在线观看| 亚洲成人久久性| 国产成人啪精品午夜网站| 一个人免费在线观看的高清视频| 中文亚洲av片在线观看爽| 少妇的丰满在线观看| 亚洲内射少妇av| 欧美不卡视频在线免费观看| 非洲黑人性xxxx精品又粗又长| 18禁国产床啪视频网站| 伊人久久精品亚洲午夜| 欧美日韩福利视频一区二区| www.色视频.com| 国产伦在线观看视频一区| 最近在线观看免费完整版| 日韩 欧美 亚洲 中文字幕| 少妇的逼水好多| 黄片大片在线免费观看| 国产一区在线观看成人免费| 亚洲av中文字字幕乱码综合| 亚洲精品久久国产高清桃花| 久久人人精品亚洲av| 国产色婷婷99| 亚洲精华国产精华精| 亚洲美女黄片视频| 免费人成在线观看视频色| 母亲3免费完整高清在线观看| 精品人妻1区二区| 午夜精品一区二区三区免费看| 啪啪无遮挡十八禁网站| 国产男靠女视频免费网站| 国产成人a区在线观看| 免费搜索国产男女视频| 人人妻人人澡欧美一区二区| 看免费av毛片| 欧美成人一区二区免费高清观看| 亚洲欧美日韩无卡精品| 国产伦精品一区二区三区四那| 99riav亚洲国产免费| 91av网一区二区| 中文字幕高清在线视频| 国产亚洲精品久久久com| avwww免费| 91在线观看av| 亚洲国产色片| 亚洲激情在线av|