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

    基于路徑估計的衛(wèi)星雙向傳播時延差估算方法

    2014-11-30 05:47:24楊文可占建偉龔航朱祥維孫廣富
    關(guān)鍵詞:電離層雙向時延

    楊文可,占建偉,龔航,朱祥維,孫廣富

    (國防科學(xué)技術(shù)大學(xué) 電子科學(xué)與工程學(xué)院,湖南 長沙,410073)

    采用靜地軌道(GEO)衛(wèi)星的衛(wèi)星雙向時間頻率傳遞(TWSTFT)是國際上最重要的實現(xiàn)遠(yuǎn)程時間比對的方法之一,而包括地球自轉(zhuǎn)效應(yīng)(Sagnac效應(yīng))修正、電離層時延差、對流層時延差等在內(nèi)的雙向空間傳播時延差是造成 TWSTFT所測量的鐘差存在偏差的主要因素[1]。實際上,由于各種攝動力的存在[2],GEO衛(wèi)星相對于地面站存在相對運動。國際權(quán)度局(BIPM)給出的TWSTFT數(shù)據(jù)處理指導(dǎo)[1]以理想GEO衛(wèi)星為基礎(chǔ),給出了計算Sagnac效應(yīng)修正量、電離層時延差的方法,該方法是國際上一般處理TWSTFT數(shù)據(jù)的方法,本文簡稱其為BIPM法。雖然其中提到了衛(wèi)星運動對雙向Sagnac效應(yīng)修正量、雙向幾何距離時延差的影響量級,但未給出估算方法。研究衛(wèi)星運動對TWSTFT影響的文獻(xiàn)中,一般只分析了衛(wèi)星運動對 Sagnac效應(yīng)修正量的影響[3?6]及雙向幾何距離差[3?4,7?8],而忽略了衛(wèi)星運動對電離層時延差、對流層時延差的影響。其中,Sagnac效應(yīng)修正量的分析一般以BIPM法中的Sagnac效應(yīng)修正公式為基礎(chǔ)[3,5?6];分析雙向幾何距離差時,一般是以兩地面站與衛(wèi)星之間的距離差除以光速估算得到地面站信號到達(dá)衛(wèi)星的時差,再結(jié)合衛(wèi)星速度進(jìn)行估算[7?8]。以上分析方法未考慮信號真實傳播路徑,未綜合考慮信號傳播方向上經(jīng)歷的Sagnac效應(yīng)、幾何路徑時延、電離層時延、對流層時延等,并且一般以模型[3]或精度為幾十千米的軌道數(shù)據(jù)[4,8]來描述衛(wèi)星運動,存在誤差。在北斗導(dǎo)航系統(tǒng)中,基于GEO導(dǎo)航衛(wèi)星的TWSTFT實現(xiàn)站間時間同步[9],且該GEO衛(wèi)星具有2 h擬合精度好于0.1 m、徑向誤差為厘米量級的廣播星歷[10]。不同于前述文獻(xiàn)所提方法,本文作者以廣播星歷數(shù)據(jù)為依托,對衛(wèi)星或地面站的坐標(biāo)進(jìn)行旋轉(zhuǎn)來修正Sagnac效應(yīng),獲知真實信號傳播路徑,并綜合考慮幾何路徑時延、電離層時延修正、對流層時延修正等對衛(wèi)星與地面站之間的單向空間傳播時延建模,提出一種利用迭代法分別計算上行、下行信號時延,并精確計算TWSTFT中雙向空間傳播時延差的方法。

    1 衛(wèi)星雙向時延模型

    當(dāng)考慮衛(wèi)星相對于地面站的相對運動時,衛(wèi)星雙向時間頻率傳遞的時延模型如圖1所示。

    圖1中,τUi為地面站i(i=1, 2)到衛(wèi)星的上行信號空間傳播時延,包含4部分:由地面站i到衛(wèi)星的偽幾何距離時延τUi|SatP、地球自轉(zhuǎn)效應(yīng)(Sagnac效應(yīng))時延修正τSCUi、電離層時延修正τUi|Ion、對流層時延修正τUi|Trop;τDi為衛(wèi)星到地面站i的下行信號空間傳播時延,也包含4部分:由衛(wèi)星到地面站i的偽幾何距離時延τDi|SatP、Sagnac效應(yīng)時延修正τSCDi、電離層時延修正τDi|Ion、對流層時延修正τDi|Trop;τTi和τRi分別為地面站i發(fā)射信號時延和接收信號時延;τSij(i=1, 2,j=1, 2,i≠j)為衛(wèi)星由地面站i向地面站j轉(zhuǎn)發(fā)信號時延。

    其中,之所以稱τUi|SatP和τDi|SatP為偽幾何距離時延,是由于該幾何距離時延未考慮 Sagnac效應(yīng)的影響,不同于將偽幾何距離時延經(jīng)過Sagnac效應(yīng)修正得到的地面站到衛(wèi)星(或衛(wèi)星到地面站)真實傳播幾何路徑時延,記為τUi|Sat(或τDi|Sat)。

    圖1 衛(wèi)星雙向時間頻率傳遞的時延模型示意圖Fig. 1 Schematic diagram of two-way satellite time and frequency transfer with satellite in motion

    當(dāng)記以鐘i為參考測量所得的地面站j的信號時延為PRij(i=1, 2,j=1, 2,i≠j)時,鐘1與鐘2之間的時差,記為T1?T2,可表示為[1]:

    式中:0.5((τT1?τR1)?(τT2?τR2))為雙向地面站發(fā)射時延與接收時延差;0.5(τS12?τS21)為雙向衛(wèi)星轉(zhuǎn)發(fā)時延差,這兩者可通過標(biāo)定獲得[1,11]。因此,為精確計算雙向時延差,即式(1)等式右邊除雙向測量時差0.5(PR21?PR12)以外,還需確定雙向空間傳播時延差 0.5((τU1?τD1)?(τU2?τD2)),記為τUD,且可分解為

    2 算法原理

    對下行或上行信號傳播時延的估算是一個基于信號路徑估計的迭代計算的過程,基本流程如圖2所示。

    圖2 單向信號傳播時延迭代計算流程Fig. 2 One-way signal propagation delay iterative calculation flow

    2.1 下行信號傳播時延模型

    如圖1所示,以衛(wèi)星到地面站1的下行鏈路為例,下行信號時延τDn1,可表示為

    其中,地面站1信號接收時延可通過標(biāo)校獲得[11]。以下分析除τR1以外的時延計算模型,即

    2.1.1 Sagnac效應(yīng)修正

    在式(2)中對下行信號時延的 Sagnac效應(yīng)的修正是以時延修正量τSCD1的方式表示的。然而,本質(zhì)上,下行信號時延的 Sagnac效應(yīng)修正是為了獲得衛(wèi)星到地面站方向的信號傳播路徑。為此,本文對衛(wèi)星位置坐標(biāo)進(jìn)行修正,獲得信號實際傳播路徑,具體過程如下。

    (1) 假設(shè)下行信號時延為τDn1,地面站1接收信號時刻為tr1,則可知衛(wèi)星轉(zhuǎn)發(fā)信號時刻ttr21為

    (2) 對衛(wèi)星廣播星歷進(jìn)行插值,獲得在ttr21時刻的地心地固(ECEF)坐標(biāo)系下衛(wèi)星位置{xs,ys,zs},記為S;

    (4) 地面站1在tr1的ECEF坐標(biāo)系下的坐標(biāo){x1,y1,z1},記為GS1,與衛(wèi)星位置S′之間的連線方向是實際信號傳播路徑。其中,GS1為已知量。另外,據(jù)式(6),S′由tr1,ttr21和S共同決定。已知tr1,同時,當(dāng)星歷已知時,S實際上由ttr21決定。由此,S′是ttr21的函數(shù),且τR1已知,則據(jù)式(5)和式(4)可知S′實際上是τD1的函數(shù),記為

    2.1.2 幾何距離時延

    當(dāng)經(jīng)過 Sagnac效應(yīng)修正后,偽幾何距離時延τD1|SatP修正為幾何距離時延τD1|Sat,即

    其中:c為光在真空中傳播的速度,為299 792 458 m/s。

    可見:幾何距離時延τD1|Sat是地面站1坐標(biāo)GS1與經(jīng) Sagnac效應(yīng)修正的轉(zhuǎn)發(fā)信號時刻衛(wèi)星坐標(biāo)S′的函數(shù)。若隱去已知量,結(jié)合式(7),則τD1|Sat也可簡記為τD1的函數(shù):

    2.1.3 電離層時延修正

    由于電離層折射引起的偽碼相位時延修正量τD1|Ion可寫為[1]

    其中:s為傾斜因子,由信號傳播路徑?jīng)Q定,是經(jīng)Sagnac效應(yīng)修正的轉(zhuǎn)發(fā)信號時刻ttr21衛(wèi)星坐標(biāo)S′和地面站 1坐標(biāo)GS1的函數(shù);TEC為地面站1頂端垂直方向上總電子數(shù),可采用IGS提供的電離層網(wǎng)格參數(shù)[13]進(jìn)行時間和空間上的線性插值計算得到,是衛(wèi)星轉(zhuǎn)發(fā)信號時刻ttr21和地面站1坐標(biāo)GS1的函數(shù);fd為下行信號載波頻率,為已知量。

    若隱去已知量,結(jié)合式(6),則τD1|Ion同樣可簡記為τD1的函數(shù):

    2.1.4 對流層時延修正

    對流層時延修正模型有多種,本文使用 Collins模型,則τD1|Trop可寫為[14]

    其中:h為地面站1高度;ε為地面站1對衛(wèi)星的仰角。

    由于地面站1高度h可由其ECEF坐標(biāo)GS1轉(zhuǎn)換得到,地面站1對衛(wèi)星仰角ε是GS1和S′的函數(shù)。則隱去已知量,結(jié)合式(7),τD1|Trop可簡寫為

    2.2 上行信號空間傳播時延

    仍以地面站2到衛(wèi)星到地面站1的鏈路為例,說明當(dāng)計算得到衛(wèi)星到地面站1的下行信號空間傳播時延τD1后,如何計算地面站 2到衛(wèi)星的上行信號空間傳播時延,記為τU2。

    上行信號時延τUp2為

    其中,地面站2信號發(fā)送時延可通過標(biāo)校獲得[11]。以下,分析除τT2以外的時延計算模型,即

    當(dāng)忽略衛(wèi)星轉(zhuǎn)發(fā)信號時延,及期間的衛(wèi)星位置改變,認(rèn)為衛(wèi)星轉(zhuǎn)發(fā)信號時刻的位置即為衛(wèi)星接收信號時刻的位置,則已知衛(wèi)星轉(zhuǎn)發(fā)信號時刻ttr21和該時刻ECEF坐標(biāo)系下的坐標(biāo)為S,則可按照類似下行單向時延的模型對上行單向時延進(jìn)行建模。

    區(qū)別僅在于,修正Sagnac效應(yīng)時,需修正地面站2(信號發(fā)射方)的位置,獲知信號實際傳播路徑,其過程如下。

    (1) 已知衛(wèi)星轉(zhuǎn)發(fā)信號時刻ttr21,則可知地面站2發(fā)射信號時刻tt2為:

    (2) 地球旋轉(zhuǎn)速率設(shè)為,將tt2時刻地面站2在ECEF坐標(biāo)系下的坐標(biāo){x2,y2,z2},記為GS2,修正至衛(wèi)星轉(zhuǎn)發(fā)信號時刻ttr21的ECEF坐標(biāo)系下坐標(biāo)},記為:

    (3) 衛(wèi)星轉(zhuǎn)發(fā)信號時刻ttr21的ECEF坐標(biāo)系下地面站2的坐標(biāo)到衛(wèi)星位置坐標(biāo)S之間的連線方向為實際信號傳播路徑。其中,據(jù)式(21),由ttr21,tt2和GS2共同決定。而GS2為已知量,ttr21已知,tt2由τUp2和ttr21共同決定。隱去已知量,GS2 ′實際上是τU2的函數(shù)。

    按照類似于下行單向時延的建模和推導(dǎo)過程,可得到:

    其中:

    即上行信號幾何距離時延τU2|Sat,電離層時延修正τU2|Ion,對流層時延修正τU2|Trop都是τU2的函數(shù)。

    2.3 迭代法計算

    將式(9),(11),(13)代入式(4),得到:

    另外,式(23)與式(21)具有類似的形式,可通過迭代法[15]求等式(21)和(23)的解τD1和τU2。以求解式(21)的解τD1的過程進(jìn)行說明:

    (1) 設(shè)定初始解τD1(0)=130 ms:

    (2) 第j(j≥1)次迭代計算過程如下:

    (3) 設(shè)置門限λ(λ>0),收斂條件為:

    其中,門限選取為時延有效數(shù)據(jù)位量級的1/10。若有效數(shù)據(jù)位量級為1 ps,則λ=0.1 ps。

    (4) 當(dāng)滿足式(16)的收斂條件時,停止迭代計算,則

    當(dāng)計算得到地面站2到衛(wèi)星再到地面站1以及地面站1到衛(wèi)星再到地面站2的雙向鏈路的下行單向時延τDn1和τDn2以及上行單向時延τUp1和τUp2時,可求得雙向空間傳播時延差τUD。

    3 實驗結(jié)果

    本文以北斗導(dǎo)航系統(tǒng)中相距3 376 km的地面站1、地面站2之間通過位于東經(jīng)140°的GEO衛(wèi)星建立的C波段TWSTFT鏈路的2 d實測數(shù)據(jù)為例,采用本文方法和BIPM法計算雙向傳播時延差,并分析兩者的差異。

    圖3所示為GEO衛(wèi)星的星下點軌跡。由圖3可見:GEO衛(wèi)星的星下點軌跡具有以d為單位的周期波動(以下簡稱日波動),經(jīng)度變化幅度為0.09°,緯度變化幅度為3.18°。

    圖3 GEO衛(wèi)星的星下點軌跡Fig. 3 Sub-satellite point track of GEO satellite

    雙向幾何距離時延差τUD|Sat包含偽幾何距離時延差τUD|SatP和Sagnac效應(yīng)修正量τUD|SC。如圖4所示,使用BIPM法,即認(rèn)為GEO衛(wèi)星相對地面站無相對運動時,τUD|Sat為定值?84.149 ns;使用本文方法,即考慮衛(wèi)星運動以及信號傳播路徑時,發(fā)現(xiàn)τUD|Sat呈現(xiàn)幅度為0.412 ns的日波動,均值為?84.266 ns,與BIPM法結(jié)果偏差77 ps。

    圖4 雙向幾何距離時延差(包含雙向Sagnac效應(yīng)修正)τUD|SatFig. 4 Two-way geometry signal path difference (including two-way Sagnac effect correction) τUD|Sat

    圖5 雙向電離層時延的雙差(本文方法與BIPM法計算所得雙向電離層時延差τUD|Ion 的差值)Fig. 5 Double difference of two-way ionospheric corrections(the difference of two-way ionospheric correction differenceτUD| Ion calculated by method introduced in this work and BIPM method)

    圖 5所示為本文方法計算得雙向電離層時延差τUD|Ion與BIPM法計算得雙向電離層時延差τUD|Ion的差值。由于衛(wèi)星運動造成地面站到衛(wèi)星,以及衛(wèi)星到地面站的雙向信號路徑改變,使得雙向電離層時延修正量較衛(wèi)星相對地面站靜止的情況有所改變。并且,由于電離層電子濃度在白天較大,在夜晚較小,使得該差值呈現(xiàn)日波動,且幅值白天較大,夜晚較小,最大幅度為0.120 ns,均值為31 ps。

    雙向?qū)α鞑顣r延差如圖6所示。從圖6可見:當(dāng)不考慮衛(wèi)星相對地面站運動時(BIPM法),雙向?qū)α鲗訒r延差τUD|Trop幾乎為0;使用本文方法計算,發(fā)現(xiàn)D|Trop也呈現(xiàn)日波動,幅度為0.3 ps,均值為?0.9 ps。

    同時考慮τUD|Sat,τUD|Ion和τUD|Trop,得到衛(wèi)星雙向傳播時延差τUD。衛(wèi)星雙向傳播時延的雙差(本文方法與 BIPM 法計算所得衛(wèi)星雙向傳播時延差τUD的差值)如圖 7所示。從圖 7可見:對比本文方法結(jié)果與BIPM 法結(jié)果,衛(wèi)星星下點軌跡的日波動引入了τUD幅度為0.406 ns,均值為?0.087 ns的日波動。

    圖6 雙向?qū)α鲗訒r延差τUD|TropFig. 6 Two-way tropospheric correction difference τUD| Trop

    圖7 衛(wèi)星雙向傳播時延的雙差(本文方法與BIPM法計算所得衛(wèi)星雙向傳播時延差τUD的差值)Fig. 7 Double difference of two-way propagation delay (the difference of non-reciprocity of satellite two-way propagation delay τUD calculated by method introduced in this work and BIPM method)

    4 結(jié)論

    (1) 提出了一種估計信號傳播路徑,結(jié)合衛(wèi)星星歷和迭代法精確計算衛(wèi)星雙向時間頻率傳遞(TWSTFT)中雙向傳播時延差的新方法。該方法的特點在于:通過修正衛(wèi)星位置或者地面站位置來進(jìn)行Sagnac效應(yīng)修正,并由此獲得信號真實傳播路徑,并以此為基礎(chǔ),綜合考慮地面站到衛(wèi)星的幾何距離時延、Sagnac效應(yīng)、電離層時延修正、對流層時延修正,計算雙向空間傳播時延差。該方法只依賴于衛(wèi)星星歷、兩地面站坐標(biāo)和接收信號時刻,可推廣應(yīng)用于采用IGSO衛(wèi)星等非靜地軌道衛(wèi)星的TWSTFT中。另外,還可繼續(xù)研究兩地面站時間與衛(wèi)星星歷的參考時間之差造成的影響。

    (2) 本文方法能夠精確估計衛(wèi)星雙向傳播時延差,反映衛(wèi)星運動對TWSTFT的影響;并發(fā)現(xiàn)衛(wèi)星運動不但引起雙向幾何距離時延差(包含 Sagnac效應(yīng)修正)幅度達(dá)0.412 ns的日波動和均值達(dá)77 ps的偏差,對于實驗中的C波段TWSTFT鏈路而言,還造成了雙向電離層時延差幅度達(dá)0.120 ns的日波動和均值達(dá)31 ps的偏差,必須從雙向測量結(jié)果中扣除,否則將在鐘差中疊加具有偏差的日波動,影響對站間鐘差的評估。

    [1]Bauch A, Piester D, Fujieda M, et al. Directive for operational use and data handling two-way satellite time and frequency transfer (TWSTFT)[R]. Paris: BIPM, 2011: 6?8.

    [2]周善石. 基于區(qū)域監(jiān)測網(wǎng)的衛(wèi)星導(dǎo)航系統(tǒng)精密定軌方法研究[D]. 上海: 中國科學(xué)院研究生院, 2011: 31?36.ZHOU Shanshi. Studies on precise orbit determination theory and application for satellite navigation system with regional tracking network[D]. Shanghai: Graduate University of Chinese Academy of Sciences, 2011: 31?36.

    [3]Hongwei S, Imae M, Gotoh T. Impact of satellite motion on two-way satellite time and frequency transfer[J]. Electronics Letters, 2003, 39(5): 482?483.

    [4]Piester D, Bauch A, Fujieda M, et al. Studies on instabilities in long-baseline two-way satellite time and frequency transfer(TWSTFT) including a troposphere delay model[C]// The 39th Annual Precise Time and Time Interval Systems and Applications Meeting. Long Beach, 2007: 211?222.

    [5]武文俊, 李志剛, 楊旭海, 等. 衛(wèi)星雙向時間頻率傳遞中的Sagnac效應(yīng)[J]. 宇航學(xué)報, 2012, 33(7): 936?941.WU Wenjun, LI Zhigang, YANG Xuhai, et al. Sagnac effect on two-way satellite time and frequency transfer[J]. Journal of Astronautics, 2012, 33(7): 936?941.

    [6]Tseng W H, Feng K M, Lin S Y, et al. Sagnac effect and diurnal correction on two-way satellite time transfer[J]. IEEE Transactions on Instrumentation and Measurement, 2011, 60(7):2298?2303.

    [7]Parker T E, Zhang V. Sources of instabilities in two-way satellite time transfer[C]// Proceedings of the 2005 IEEE International Frequency Control Symposium and Exposition. Vancouver, 2005:745?751.

    [8]楊旭海, 李志剛, 李偉超, 等. 基于小傾角 IGSO衛(wèi)星的雙向時間傳遞試驗[C]// 2007全國時間頻率學(xué)術(shù)會議論文集. 宜昌,2007: 372?376.YANG Xuhai, LI Zhigang, LI Weichao, et al. Two-way time transfer experiment based on small obliquity IGSO satellite[C]//2007 National Time and Frequency Research Conference Proceeding. Yichang, 2007: 372?376.

    [9]HAN Chunhao, YANG Yuanxi, CAI Zhiwu. Beidou navigation satellite system and its time scales[J]. Metrologia, 2011, 48(4):s213?s218.

    [10]黃勇, 胡小工, 王小亞, 等. 中高軌衛(wèi)星廣播星歷精度分析[J].天文學(xué)進(jìn)展, 2006, 24(1): 81?87.HUANG Yong, HU Xiaogong, WANG Xiaoya, et al. Precision analysis of broadcast ephemeris for medium and high orbit satellites[J]. Progress in Astronomy, 2006, 24(1): 81?87.

    [11]李星, 耿淑敏,李垣陵. 雙向時間同步系統(tǒng)的設(shè)備時延校準(zhǔn)技術(shù)研究[J]. 中國空間科學(xué)技術(shù), 2011, 31(4): 23?28.LI Xing, GENG Shumin, LI Yuanling. Equipment delay calibration technique for two-way time synchronization system[J]. Chinese Space Science and Technology, 2011, 31(4):23?28.

    [12]Kaplan E D, Hegerty C J. GPS原理與應(yīng)用[M]. 2版. 寇艷紅,譯. 北京: 電子工業(yè)出版社, 2008: 229.Kaplan E D, Hegerty C J. Understanding GPS principles and applications[M]. 2nd ed. KOU Yanhong, trans. Beijing:Publishing House of Electronics Industry, 2008: 229.

    [13]Kouba J. A guide to using International GNSS service (IGS)products[EB/OL]. [2009?05?13]. http://igscb.jpl.nasa.gov/igscb/resource/pubs/UsingIGSProductsVer21.pdf

    [14]Farrell J, Barth M. The global positioning system and inertial navigation[M]. New York: McGraw Hill, 1998: 156.

    [15]葉其孝, 沈永歡. 實用數(shù)學(xué)手冊[M]. 北京: 科學(xué)出版社, 2010:740.YE Qixiao, SHEN Yonghuan. Hand book of applied mathematics[M]. Beijing: Science Press, 2010: 740.

    猜你喜歡
    電離層雙向時延
    雙向度的成長與自我實現(xiàn)
    出版人(2022年11期)2022-11-15 04:30:18
    一種電離層TEC格點預(yù)測模型
    Kalman濾波估算電離層延遲的一種優(yōu)化方法
    基于GCC-nearest時延估計的室內(nèi)聲源定位
    電子制作(2019年23期)2019-02-23 13:21:12
    基于改進(jìn)二次相關(guān)算法的TDOA時延估計
    電離層對中高軌SAR影響機理研究
    FRFT在水聲信道時延頻移聯(lián)合估計中的應(yīng)用
    基于分段CEEMD降噪的時延估計研究
    一種軟開關(guān)的交錯并聯(lián)Buck/Boost雙向DC/DC變換器
    一種工作頻率可變的雙向DC-DC變換器
    亚洲av成人精品一区久久| 色综合婷婷激情| 69av精品久久久久久| 老司机午夜福利在线观看视频| 18+在线观看网站| 一级毛片久久久久久久久女| 99国产极品粉嫩在线观看| 国产精品一区www在线观看 | 国产高清激情床上av| 简卡轻食公司| 成人高潮视频无遮挡免费网站| 久久人人精品亚洲av| 天天躁日日操中文字幕| 日韩国内少妇激情av| 午夜福利在线观看吧| 动漫黄色视频在线观看| 一本一本综合久久| av福利片在线观看| 性色avwww在线观看| 国产探花在线观看一区二区| 色噜噜av男人的天堂激情| 亚洲一区高清亚洲精品| 欧美日韩乱码在线| 99国产极品粉嫩在线观看| xxxwww97欧美| 成人特级黄色片久久久久久久| 观看美女的网站| 国产精品久久久久久亚洲av鲁大| 亚洲欧美精品综合久久99| 九色成人免费人妻av| 国产伦在线观看视频一区| 久久精品综合一区二区三区| 日本成人三级电影网站| 全区人妻精品视频| 午夜精品在线福利| av.在线天堂| 国产精品久久久久久久久免| 国产视频内射| 精品无人区乱码1区二区| 一个人观看的视频www高清免费观看| 亚洲成人久久性| 亚洲午夜理论影院| 性欧美人与动物交配| 最近在线观看免费完整版| 亚洲狠狠婷婷综合久久图片| 亚洲美女视频黄频| 久久欧美精品欧美久久欧美| 国产免费男女视频| 91在线观看av| 国产免费男女视频| 亚洲最大成人中文| 中国美白少妇内射xxxbb| 国产午夜精品论理片| 乱码一卡2卡4卡精品| 成熟少妇高潮喷水视频| 内地一区二区视频在线| 亚洲天堂国产精品一区在线| 色av中文字幕| 成年女人毛片免费观看观看9| 亚洲av一区综合| 久久精品影院6| 国产精品免费一区二区三区在线| 亚洲国产高清在线一区二区三| 真实男女啪啪啪动态图| 亚洲天堂国产精品一区在线| 国产精品久久久久久久久免| 亚洲成人久久爱视频| 欧美一区二区亚洲| 99久久成人亚洲精品观看| 国产欧美日韩一区二区精品| 国产精品爽爽va在线观看网站| 真人做人爱边吃奶动态| 亚洲av不卡在线观看| 中出人妻视频一区二区| 色视频www国产| 国产日本99.免费观看| 国产成人aa在线观看| 日本 欧美在线| 国内精品宾馆在线| 欧美一区二区亚洲| 日韩高清综合在线| 国产精品爽爽va在线观看网站| 嫩草影视91久久| 国产一区二区在线观看日韩| 欧美黑人欧美精品刺激| 女人十人毛片免费观看3o分钟| 桃红色精品国产亚洲av| 成年女人毛片免费观看观看9| 最近最新中文字幕大全电影3| 午夜免费男女啪啪视频观看 | 成人综合一区亚洲| 精品福利观看| 午夜免费成人在线视频| 国产精品免费一区二区三区在线| 丰满乱子伦码专区| 97热精品久久久久久| 97超级碰碰碰精品色视频在线观看| 九九爱精品视频在线观看| 人人妻人人看人人澡| 九九在线视频观看精品| 91av网一区二区| 国产一区二区在线av高清观看| 丝袜美腿在线中文| 99热6这里只有精品| 男女那种视频在线观看| 日韩国内少妇激情av| 国产精品亚洲一级av第二区| 波多野结衣高清无吗| 久久久久久久久久黄片| 亚洲av第一区精品v没综合| 简卡轻食公司| 最近最新中文字幕大全电影3| 国产午夜精品论理片| ponron亚洲| 国产蜜桃级精品一区二区三区| 噜噜噜噜噜久久久久久91| 久久午夜亚洲精品久久| 精品久久久久久成人av| 亚洲成av人片在线播放无| 精品一区二区三区视频在线观看免费| 深夜a级毛片| 九九热线精品视视频播放| 中国美女看黄片| 欧美3d第一页| 成人一区二区视频在线观看| av福利片在线观看| 精品久久久久久久久亚洲 | 免费在线观看日本一区| 变态另类成人亚洲欧美熟女| 亚洲精品成人久久久久久| 免费高清视频大片| 一本久久中文字幕| 亚洲av不卡在线观看| 色5月婷婷丁香| 中文字幕av成人在线电影| 五月玫瑰六月丁香| 亚洲va在线va天堂va国产| 人妻久久中文字幕网| 欧美日本视频| 两个人的视频大全免费| 99久久中文字幕三级久久日本| 搡老妇女老女人老熟妇| 亚洲成a人片在线一区二区| 精品国内亚洲2022精品成人| 女同久久另类99精品国产91| 久久亚洲精品不卡| 男女之事视频高清在线观看| 免费在线观看成人毛片| av专区在线播放| 免费观看在线日韩| 91麻豆精品激情在线观看国产| 天美传媒精品一区二区| 99热这里只有是精品在线观看| 精品久久久久久久久久免费视频| 日本熟妇午夜| www.www免费av| 美女高潮喷水抽搐中文字幕| 国产伦人伦偷精品视频| 亚洲人成伊人成综合网2020| 日日夜夜操网爽| 一级黄色大片毛片| 国内揄拍国产精品人妻在线| av视频在线观看入口| 国产单亲对白刺激| 国产亚洲精品综合一区在线观看| 在线观看66精品国产| 少妇丰满av| 最近中文字幕高清免费大全6 | 欧美日本视频| 国产精品久久久久久亚洲av鲁大| 日本精品一区二区三区蜜桃| 日韩欧美精品免费久久| 国内久久婷婷六月综合欲色啪| 久久精品国产亚洲av涩爱 | 国产色婷婷99| 亚洲一区二区三区色噜噜| 久久久久久久久久黄片| 亚洲欧美清纯卡通| 天堂影院成人在线观看| 久久久久性生活片| 国产色婷婷99| av在线蜜桃| АⅤ资源中文在线天堂| 亚洲欧美日韩高清专用| 国产高清三级在线| 夜夜看夜夜爽夜夜摸| 99久久中文字幕三级久久日本| 亚洲aⅴ乱码一区二区在线播放| 桃色一区二区三区在线观看| 91在线观看av| 亚洲精品一区av在线观看| 亚洲最大成人av| bbb黄色大片| 精品久久久久久久人妻蜜臀av| 国产午夜福利久久久久久| 国产成年人精品一区二区| 美女cb高潮喷水在线观看| 国产色爽女视频免费观看| 夜夜看夜夜爽夜夜摸| 88av欧美| 免费观看的影片在线观看| 久久久久久久亚洲中文字幕| 亚洲中文字幕日韩| 88av欧美| 精品一区二区三区视频在线| 大又大粗又爽又黄少妇毛片口| 成人特级黄色片久久久久久久| 久久亚洲精品不卡| 亚洲精华国产精华精| 日本三级黄在线观看| 一级黄色大片毛片| 精品乱码久久久久久99久播| 欧美成人一区二区免费高清观看| 琪琪午夜伦伦电影理论片6080| 久久精品久久久久久噜噜老黄 | 老司机福利观看| 伊人久久精品亚洲午夜| 五月伊人婷婷丁香| 欧美zozozo另类| 欧美国产日韩亚洲一区| 日本成人三级电影网站| 精品乱码久久久久久99久播| 久久久色成人| 国产精品野战在线观看| 神马国产精品三级电影在线观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲欧美日韩高清专用| 国产精品野战在线观看| 欧美又色又爽又黄视频| 美女xxoo啪啪120秒动态图| 久久草成人影院| 99久久九九国产精品国产免费| 亚洲欧美日韩无卡精品| 久久久久久久精品吃奶| 日本在线视频免费播放| 在线天堂最新版资源| 国产精品人妻久久久影院| 99在线视频只有这里精品首页| 亚洲在线自拍视频| 欧美一区二区精品小视频在线| 国产精品三级大全| 两人在一起打扑克的视频| 伊人久久精品亚洲午夜| 久久这里只有精品中国| 俺也久久电影网| 99精品在免费线老司机午夜| 两个人的视频大全免费| 精品久久久久久久久av| 精品人妻熟女av久视频| 欧美高清成人免费视频www| 亚洲自拍偷在线| 欧美激情在线99| 中文在线观看免费www的网站| 国产高清不卡午夜福利| 亚洲熟妇中文字幕五十中出| 国产中年淑女户外野战色| 特级一级黄色大片| 男女下面进入的视频免费午夜| 欧美日韩乱码在线| 成人二区视频| 欧美绝顶高潮抽搐喷水| 久久国产乱子免费精品| 欧美在线一区亚洲| 午夜福利成人在线免费观看| 内射极品少妇av片p| 精品人妻视频免费看| 两个人视频免费观看高清| 乱码一卡2卡4卡精品| 国产精品一区二区性色av| a级毛片免费高清观看在线播放| 99热这里只有精品一区| 自拍偷自拍亚洲精品老妇| 久久久色成人| 熟妇人妻久久中文字幕3abv| 有码 亚洲区| 级片在线观看| 国产免费av片在线观看野外av| 18禁黄网站禁片午夜丰满| 麻豆成人av在线观看| 午夜亚洲福利在线播放| 国产午夜精品论理片| 男女边吃奶边做爰视频| 国产色婷婷99| 亚洲国产欧洲综合997久久,| 99热这里只有是精品50| 在线免费观看的www视频| 亚洲男人的天堂狠狠| 人人妻,人人澡人人爽秒播| 午夜精品在线福利| 国产亚洲精品久久久久久毛片| 少妇裸体淫交视频免费看高清| 亚洲在线观看片| 国产伦精品一区二区三区视频9| 99热6这里只有精品| 国产一区二区三区视频了| 亚洲国产精品成人综合色| 久久久久免费精品人妻一区二区| 国产一区二区在线av高清观看| 亚洲精品乱码久久久v下载方式| 又粗又爽又猛毛片免费看| 一本久久中文字幕| 国产精品亚洲一级av第二区| videossex国产| 亚洲va日本ⅴa欧美va伊人久久| 有码 亚洲区| 亚洲av成人精品一区久久| 免费av毛片视频| 91麻豆av在线| 亚洲成人中文字幕在线播放| 中出人妻视频一区二区| 精华霜和精华液先用哪个| 夜夜夜夜夜久久久久| 91午夜精品亚洲一区二区三区 | 日韩欧美国产一区二区入口| 国产美女午夜福利| 久久精品人妻少妇| 国产麻豆成人av免费视频| 人妻夜夜爽99麻豆av| 春色校园在线视频观看| 日韩中文字幕欧美一区二区| 国产精华一区二区三区| 一个人观看的视频www高清免费观看| 日日夜夜操网爽| 亚洲综合色惰| 又爽又黄a免费视频| 久久精品夜夜夜夜夜久久蜜豆| 麻豆精品久久久久久蜜桃| 在线观看av片永久免费下载| 91午夜精品亚洲一区二区三区 | 伦精品一区二区三区| 在线a可以看的网站| 国产精品一及| 天天躁日日操中文字幕| 一夜夜www| 乱码一卡2卡4卡精品| 日本-黄色视频高清免费观看| 哪里可以看免费的av片| 老司机午夜福利在线观看视频| 免费无遮挡裸体视频| 天堂√8在线中文| 91麻豆av在线| 亚洲成人久久性| 亚洲aⅴ乱码一区二区在线播放| 又黄又爽又刺激的免费视频.| 亚洲精品亚洲一区二区| 日韩av在线大香蕉| 免费一级毛片在线播放高清视频| 国产真实伦视频高清在线观看 | 日韩欧美在线乱码| 日本 av在线| 亚洲成人久久爱视频| 精品福利观看| 丰满乱子伦码专区| 99久久成人亚洲精品观看| av天堂在线播放| 男人的好看免费观看在线视频| 中文字幕av在线有码专区| www.色视频.com| 桃色一区二区三区在线观看| 亚洲国产高清在线一区二区三| 久久久久久久午夜电影| 亚洲图色成人| 男女啪啪激烈高潮av片| 最近最新免费中文字幕在线| 少妇的逼水好多| 香蕉av资源在线| 美女高潮的动态| 国产亚洲精品久久久久久毛片| 成人特级av手机在线观看| 级片在线观看| 婷婷色综合大香蕉| 国产成人影院久久av| 日日摸夜夜添夜夜添小说| 在线观看一区二区三区| 国产欧美日韩一区二区精品| 成人特级黄色片久久久久久久| 久久国产精品人妻蜜桃| ponron亚洲| 最近最新免费中文字幕在线| 国产精品人妻久久久久久| 久久久成人免费电影| 在线看三级毛片| 亚洲精华国产精华液的使用体验 | 中亚洲国语对白在线视频| 色5月婷婷丁香| 99九九线精品视频在线观看视频| 欧美高清成人免费视频www| 亚洲最大成人av| 中文资源天堂在线| 俄罗斯特黄特色一大片| 国内精品一区二区在线观看| 日本 欧美在线| 国产老妇女一区| 欧美zozozo另类| 国产精品不卡视频一区二区| 亚洲av.av天堂| 欧美中文日本在线观看视频| 一级黄色大片毛片| 欧美日韩乱码在线| 午夜视频国产福利| 国产精品人妻久久久影院| 国产69精品久久久久777片| 一区二区三区激情视频| 亚洲七黄色美女视频| 免费av观看视频| 国产高清视频在线观看网站| 久久久久性生活片| 国产视频一区二区在线看| 乱人视频在线观看| 精品久久久噜噜| 五月伊人婷婷丁香| 国产麻豆成人av免费视频| 欧美区成人在线视频| 精品免费久久久久久久清纯| 99久久精品一区二区三区| 免费av毛片视频| 色av中文字幕| 午夜免费成人在线视频| 亚洲最大成人手机在线| 成年版毛片免费区| 99九九线精品视频在线观看视频| 国产又黄又爽又无遮挡在线| 18禁黄网站禁片免费观看直播| 成人性生交大片免费视频hd| 国产av在哪里看| 91麻豆av在线| 精品久久久噜噜| 少妇高潮的动态图| 美女cb高潮喷水在线观看| 国产一级毛片七仙女欲春2| 久久九九热精品免费| 白带黄色成豆腐渣| av在线天堂中文字幕| a级一级毛片免费在线观看| 色综合亚洲欧美另类图片| 亚洲一区高清亚洲精品| 国产女主播在线喷水免费视频网站 | 亚洲黑人精品在线| 韩国av一区二区三区四区| 免费不卡的大黄色大毛片视频在线观看 | 别揉我奶头~嗯~啊~动态视频| 亚洲第一区二区三区不卡| 可以在线观看的亚洲视频| 久久欧美精品欧美久久欧美| 欧美日韩亚洲国产一区二区在线观看| 久久这里只有精品中国| 国内精品久久久久精免费| 琪琪午夜伦伦电影理论片6080| 18+在线观看网站| 日本a在线网址| av天堂中文字幕网| 在线免费十八禁| 久久欧美精品欧美久久欧美| 国内精品久久久久久久电影| 久久精品人妻少妇| 午夜爱爱视频在线播放| 精品一区二区三区视频在线观看免费| 亚洲四区av| 亚洲一区二区三区色噜噜| 成熟少妇高潮喷水视频| 色综合亚洲欧美另类图片| 中文字幕久久专区| 亚洲av一区综合| 国产精品98久久久久久宅男小说| 午夜老司机福利剧场| 三级国产精品欧美在线观看| 久久精品国产亚洲网站| 国产高清激情床上av| 亚洲av成人精品一区久久| 淫秽高清视频在线观看| 免费无遮挡裸体视频| 一本一本综合久久| 村上凉子中文字幕在线| 尤物成人国产欧美一区二区三区| 国产欧美日韩精品亚洲av| 亚洲国产日韩欧美精品在线观看| 欧美zozozo另类| 精品一区二区三区视频在线观看免费| 熟女人妻精品中文字幕| 18禁黄网站禁片免费观看直播| 亚洲不卡免费看| 久久精品影院6| 国产成人一区二区在线| www.www免费av| 又粗又爽又猛毛片免费看| 特大巨黑吊av在线直播| 亚洲国产欧洲综合997久久,| 亚洲,欧美,日韩| 欧美xxxx黑人xx丫x性爽| ponron亚洲| 一级av片app| 国产精品,欧美在线| 中文字幕av在线有码专区| 丝袜美腿在线中文| 国产高清视频在线播放一区| 特大巨黑吊av在线直播| 免费高清视频大片| 免费看光身美女| 男插女下体视频免费在线播放| 一区二区三区高清视频在线| 婷婷丁香在线五月| 亚洲18禁久久av| 国产精品亚洲美女久久久| www.www免费av| 波多野结衣高清无吗| 国产久久久一区二区三区| 十八禁网站免费在线| 国产免费男女视频| 国产亚洲精品久久久com| 国产三级在线视频| 老师上课跳d突然被开到最大视频| 国产成人福利小说| 伦精品一区二区三区| 久久天躁狠狠躁夜夜2o2o| 国内少妇人妻偷人精品xxx网站| 亚洲一区高清亚洲精品| 国产精品,欧美在线| 在线观看舔阴道视频| 成人无遮挡网站| 亚洲 国产 在线| 亚洲美女搞黄在线观看 | 免费观看在线日韩| 亚洲最大成人av| 国产精品98久久久久久宅男小说| 真人一进一出gif抽搐免费| 精品人妻1区二区| 在线播放国产精品三级| 91久久精品国产一区二区三区| 欧美激情国产日韩精品一区| 给我免费播放毛片高清在线观看| 亚洲av中文av极速乱 | 日日撸夜夜添| netflix在线观看网站| 少妇猛男粗大的猛烈进出视频 | 亚洲av二区三区四区| 18禁黄网站禁片午夜丰满| 男女下面进入的视频免费午夜| 亚洲在线自拍视频| 一夜夜www| 精品一区二区三区人妻视频| 亚洲av不卡在线观看| 床上黄色一级片| 人妻制服诱惑在线中文字幕| 亚洲欧美日韩高清专用| 少妇丰满av| 看片在线看免费视频| 亚洲五月天丁香| 午夜福利在线在线| 国产综合懂色| 一个人免费在线观看电影| 色哟哟哟哟哟哟| 色av中文字幕| 成人特级av手机在线观看| av天堂中文字幕网| 国产亚洲欧美98| 亚洲18禁久久av| 国产毛片a区久久久久| 能在线免费观看的黄片| 久久99热这里只有精品18| 午夜日韩欧美国产| 成人无遮挡网站| 亚洲美女黄片视频| 成人综合一区亚洲| 久久亚洲真实| 国产av在哪里看| 国产成人a区在线观看| 亚洲av一区综合| 亚洲成人久久性| 欧洲精品卡2卡3卡4卡5卡区| 中文字幕熟女人妻在线| 免费观看在线日韩| 成熟少妇高潮喷水视频| 成人特级av手机在线观看| 精品不卡国产一区二区三区| 最新在线观看一区二区三区| 波多野结衣高清作品| 色综合色国产| 免费在线观看成人毛片| 深夜精品福利| 老熟妇乱子伦视频在线观看| 国内精品久久久久久久电影| 深夜精品福利| 成年版毛片免费区| 97超级碰碰碰精品色视频在线观看| 香蕉av资源在线| 欧美绝顶高潮抽搐喷水| 精品人妻熟女av久视频| 国产精品一及| 人妻制服诱惑在线中文字幕| 三级毛片av免费| 亚洲一级一片aⅴ在线观看| 国产午夜精品论理片| 精品一区二区三区视频在线| 九九久久精品国产亚洲av麻豆| 波多野结衣巨乳人妻| 午夜免费激情av| 国产一区二区在线观看日韩| 午夜日韩欧美国产| 国内精品宾馆在线| 热99在线观看视频| 国产成人一区二区在线| 国产精品一区www在线观看 | 精品国内亚洲2022精品成人| 亚洲五月天丁香| 亚洲美女黄片视频| 可以在线观看毛片的网站| 最后的刺客免费高清国语| 成人三级黄色视频| 美女被艹到高潮喷水动态| 在线观看午夜福利视频| 亚洲成人中文字幕在线播放| 少妇的逼水好多| 日本欧美国产在线视频| 婷婷六月久久综合丁香| 国产精品亚洲美女久久久|