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

    JPL行星歷表的發(fā)展及比較

    2022-08-30 05:55:36劉婉逸鄒賢才衷路萍
    大地測量與地球動力學(xué) 2022年9期
    關(guān)鍵詞:質(zhì)心月球探測器

    劉婉逸 鄒賢才 衷路萍

    1 武漢大學(xué)測繪學(xué)院,武漢市珞喻路129號,430079

    行星歷表廣泛應(yīng)用于航天器導(dǎo)航、對行星及其衛(wèi)星等的空間觀測、天文年歷編制、地球物理及引力理論驗證等領(lǐng)域[1]。目前行星和月球歷表編制主要有分析與半分析法及數(shù)值積分法兩大類[2],分析與半分析歷表主要有法國巴黎天文臺VSOP系列歷表,數(shù)值歷表主要有美國噴氣推進實驗室(JPL)發(fā)布的DE系列歷表[1,3]、歐空局(ESA)支持編制的INPOP歷表[4]以及俄羅斯應(yīng)用天文研究所發(fā)布的EPM系列歷表[5]。我國紫金山天文臺也自行編制了PMOE系列歷表,但目前其精度不及上述高精度歷表[6-8]。

    DE系列歷表發(fā)展于20世紀60年代,是以運動方程的數(shù)值積分為基礎(chǔ)推出的歷表[9],在國際上廣泛應(yīng)用。DE系列歷表中比較著名的有1984年以來廣泛使用的DE200、國際地球自轉(zhuǎn)與參考系服務(wù)(IERS)在2003、2010規(guī)范中推薦的DE405以及DE421歷表。DE405可被用于天體視位置計算,DE421歷表中月球天平動角可被用于GRAIL任務(wù)中月心慣性系和月固系間的坐標轉(zhuǎn)換。目前使用比較廣泛的為DE430,于2013-04發(fā)布,時間跨度為1550~2650年,采用IAU1980章動模型,包含月球天平動[10],曾被用于天體物理學(xué)研究。目前最新版本為DE440和DE441,后者假定月球液體核和固體幔間無阻尼,時間跨度更長[11]。

    由于觀測數(shù)據(jù)和動力學(xué)模型不同,DE系列模型的精度存在差異。本文以DE405、DE421、DE430、DE440為例給出不同行星歷表擬合深空數(shù)據(jù)適用范圍的差異,同時比較不同歷表中大行星在地心慣性系及太陽系質(zhì)心慣性系中的位置差異,最后分析各歷表用于月球月固系和慣性坐標系間坐標和速度轉(zhuǎn)換的精度。

    1 基本行星歷表

    近年來,隨著科學(xué)技術(shù)的不斷提升以及深空探測任務(wù)的不斷增多,行星歷表的更新也越來越頻繁,從而出現(xiàn)很多為特定飛行任務(wù)服務(wù)的臨時星歷。觀測數(shù)據(jù)的增多也使得歷表精度不斷提高。

    DE405[12]星歷對內(nèi)部坐標系指向進行改進,主要通過對繞金星的麥哲倫號探測器進行甚長基線干涉(very long baseline interferometry,VLBI)觀測,使其指向國際天球參考框架(ICRF),整體精度為1 mas,距離相當于1 km。DE421[13]星歷中動力學(xué)參考框架同ICRF的對準由包括火星偵察軌道器MRO、火星奧德賽(ODY)、火星全球探測者MGS在內(nèi)的火星探測器VLBI測量完成,精度為0.25 mas[1]。DE430星歷通過對火星軌道航天器進行VLBI觀測來實現(xiàn)內(nèi)部行星軌道與ICRF2的對準,精度為0.2 mas[10]。DE440星歷通過VLBI測量火星軌道航天器相對于ICRF中位置已知的類星體來實現(xiàn)同ICRF3的對準,整體精度為0.2 mas[11]。DE440星歷和DE430星歷相比,主要有以下改進:Lense-Thirring效應(yīng)、柯伊伯帶(Kuiper belt)新模型、Vondrak歲差模型、考慮地心進動對月球軌道的影響、考慮太陽輻射壓對地月系統(tǒng)軌道的影響、內(nèi)部行星動力學(xué)框架同ICRF的連接由ICRF2[14]變?yōu)镮CRF3[15]、2013~2020年新的航天器和LLR數(shù)據(jù)。

    2 擬合行星歷表的高精度觀測數(shù)據(jù)

    擬合行星歷表采用的觀測資料主要包括行星的光學(xué)觀測、地面雷達測距、行星探測器測量以及激光測月4類數(shù)據(jù)。自DE102星歷以來,DE系列歷表已擬合越來越多的深空探測器測量數(shù)據(jù),數(shù)據(jù)類型覆蓋更全面,同時動力學(xué)模型也不斷精化,這都有利于歷表精度的提高??紤]到歷表使用的廣泛性及國際標準性,本文選取4個星歷,對比其月球、火星、木星軌道擬合使用的觀測數(shù)據(jù)。具體來講,相比于DE405,DE421放棄使用部分天文觀測數(shù)據(jù),數(shù)據(jù)時間跨度增長、航天器觀測數(shù)據(jù)增多。DE430剔除了DE421中部分精度較低的觀測數(shù)據(jù),同時使用更高精度的數(shù)據(jù),通過信使號的距離測量以及2008~2009年3次飛掠的無線電跟蹤測量來改進水星軌道;對卡西尼(Cassini)探測器進行VLBI測量以及測距使得土星軌道精度有所提升。DE430中探測器觀測數(shù)據(jù)時間跨度的增加也可提升金星、地球和火星軌道精度。相比于DE430,DE440擬合朱諾號(Juno)探測器無線電測距和超長基線陣列(VLBA)測量,使得木星軌道顯著提升。得益于卡西尼探測器的VLBA測量,土星軌道也有所改進。7 a跨度的新數(shù)據(jù)使得月球、金星、水星、火星軌道精度均有所提升。利用高精度觀測數(shù)據(jù)得到行星位置和速度的過程,實際上是動力學(xué)法確定行星狀態(tài)矢量及其他參數(shù)的過程。在建立DE歷表時,考慮的因素主要有星體的點質(zhì)量、非球形引力場的影響、地球潮汐對月球加速度的影響、月球天平動、相對論加速度以及太陽輻射壓等。對于動力學(xué)模型,4種歷表均基于參數(shù)化后牛頓N體引力模型進行數(shù)值積分(PPN-EIH)。DE421、DE430、DE440模型中所有參數(shù)個數(shù)均為228個,DE405為156個。對于太陽系中的小行星,DE440除包含DE421和DE430中考慮的谷神星的牛頓力以及342顆主帶小行星外,還包含30顆柯伊伯帶星體,其余的柯伊伯帶星體被當作一個環(huán)。

    3 行星歷表比較

    3.1 不同歷表大行星位置比較

    在進行深空探測過程中,對探測器進行觀測的射電或光學(xué)望遠鏡多位于地球,觀測數(shù)據(jù)具有顯著的地基屬性,數(shù)據(jù)處理也需要使用地心天球坐標系。作為基本空間參考系之一,太陽系質(zhì)心天球坐標系也常用于深空探測。因此,比較不同星歷間大行星相對于太陽系質(zhì)心及地心的位置分量差異對于深空探測及導(dǎo)航等具有重要意義。本文計算4個不同歷表間大行星(包括太陽、冥王星、月球)在太陽系質(zhì)心以及地心慣性系中位置差異的最大值(max)及均方根(RMS),數(shù)據(jù)時間跨度為1960~2020年,時間間隔為1 d,結(jié)果見表1和表2(單位km)。本文精度分析均以DE440星歷作為參考,將DE405、DE421、DE430與之作差。Park等[11]給出DE440星歷中各天體的部分觀測數(shù)據(jù)殘差以表征星歷精度,如月球軌道的LLR數(shù)據(jù)殘差為20 cm(早期)和1.3 cm(近期),水星軌道的Messenger測距殘差為0.7 m,金星、火星、木星、土星的觀測數(shù)據(jù)殘差分別為8 m、0.7 m、13 m、3 m。Folkner等[10]也基于各種觀測數(shù)據(jù)的殘差對DE430星歷精度進行分析。但需要注意的是,上述精度結(jié)果應(yīng)乘以2倍尺度因子來合理評估軌道變化[16]。盡管本文的差值并不能代表歷表用于測量行星位置的絕對誤差,但在一定程度上可以反映歷表的相對精度。

    由表1可知,月球在所有天體中位置精度最高,DE421和DE430星歷精度為1~2 m,DE405歷表精度略差,為7 m。月球歷表的高精度得益于月球LLR測量數(shù)據(jù)的大體量和高精度。對于太陽、水星、金星、火星,DE405歷表精度均在103m量級,DE421和DE430歷表精度在102m量級(DE421星歷中水星位置精度為2 km除外),可提高一個數(shù)量級。其中,除金星外,DE430相對于DE421在精度上至少可提高50%,這是因為水星、火星的航天器測量數(shù)據(jù)時間跨度變長以及觀測精度提高。對于木星和土星,DE405歷表精度在105m量級,比DE421和DE430歷表精度高1~2個量級。其中,DE430歷表精度相對于DE421依然有50%的提高。由于卡西尼號測量,

    表1 1960~2020年大行星相對地心位置差異的最大值max和均方根RMS

    表2 1960~2020年大行星相對太陽系質(zhì)心位置差異的最大值max和均方根RMS

    土星軌道精度比木星軌道精度高。天王星、海王星歷表精度從DE405的106m量級提高到DE421和DE430的105m量級。對于冥王星,DE405和DE421歷表精度均在106m量級,DE430歷表精度有一個量級的提高。天王星、海王星和冥王星的軌道主要由天文測量觀測確定,受制于地球大氣測量以及恒星目錄的不確定性,DE430之前位置精度均限制在105~106m量級。

    由表2可知,相對于太陽系質(zhì)心的位置偏差,地球、水星、金星和火星的DE405、DE421和DE430同DE440歷表的差異均為數(shù)10 km到100 km,而DE405、DE421和DE430之間的差異在數(shù)km。以金星為例,時間跨度為1960~2020年,步長為1 d,結(jié)果如圖1所示。這是由于DE440歷表在計算太陽系質(zhì)心時加入了30個柯伊伯帶星體(KBOs)的擾動,同時將其他柯伊伯帶星體當作一個圓環(huán),其結(jié)果與DE430歷表相比,太陽系質(zhì)心移動100 km左右。天王星、冥王星和海王星相對于太陽系質(zhì)心的位置差異隨著新星歷的發(fā)布也越來越小。

    圖1 不同歷表間金星相對于太陽系質(zhì)心的位置差異Fig.1 The position difference of Venus relative to the mass center of solar system based on different ephemerides

    DE405與其他星歷的差距明顯,可能由于DE405、DE421、DE430、DE440均采用參數(shù)化后牛頓N體運動方程描述太陽系天體運動,其中DE405僅包含156個參數(shù),而DE421、DE430、DE440均包含228個參數(shù)。

    3.2 不同歷表中月球相對于地心的位置差異

    月球是地球唯一的天然衛(wèi)星,同時也是離地球最近的天體,月球超高真空、無磁場的環(huán)境可為科學(xué)研究提供良好條件,因此月球探測一直是深空探測的重要部分。隨著激光測月數(shù)據(jù)量的增多,月球軌道的測量精度越來越高。同時,隨著探月任務(wù)的發(fā)展,對月球探測器位置的精度要求也越來越高。設(shè)一個繞月探測器的地心位置矢量為rES,月心位置矢量為rLS,月球的地心位置矢量為rEL,則存在以下關(guān)系式:

    rES=rEL+rLS

    (1)

    由式(1)可知,歷表造成的月球地心慣性系下的位置和速度差異會影響繞月探測器在地心慣性系下的位置和速度。本文基于不同星歷計算和比較月球在地心天球坐標系下的位置及速度差異,分析其對繞月探測器的位置影響。

    基于不同歷表計算1960~2020年月球在J2000.0地心天球坐標系下的位置和速度,以DE440結(jié)果為參考值,研究其他不同歷表(DE405、DE421、DE430)引起的月球在地心天球坐標系下的位置和速度誤差,結(jié)果如圖2、3所示。由圖可知,DE405位置矢量誤差為7 m,DE421和DE430誤差分別為1.5 m、1.3 m。DE405速度誤差在0.02 mm/s左右,DE421和DE430速度誤差在0.005 mm/s以內(nèi)。DE430在X和Y方向的位置誤差均值比DE421小0.25 m左右??梢园l(fā)現(xiàn),相比于DE405歷表,DE421和DE430的誤差明顯減小,這是因為DE421比DE405所用數(shù)據(jù)跨度更長、精度更高。相比于DE421歷表,DE430歷表使用了更高精度的重力場數(shù)據(jù)以及更長跨度的LLR數(shù)據(jù)。在DE430基礎(chǔ)上,DE440歷表又使用7 a的LLR數(shù)據(jù),同時LLR測距數(shù)據(jù)的殘差為1.3 cm,較DE430有所提升。因此,在繞月探測任務(wù)中,若要求探測器位置精度在亞米級,則需要采用DE421和DE430,甚至DE440星歷。

    圖2 不同歷表間月球的地球質(zhì)心位置差異Fig.2 The position difference of Moon relative to the mass center of Earth based on different ephemerides

    圖3 不同歷表間月球的地球質(zhì)心速度差異Fig.3 The velocity difference of Moon relative to the mass center of Earth based on different ephemerides

    3.3 不同歷表對月固坐標系坐標的影響

    在月球參考系中建立運動方程時,通常選擇月心慣性系,而月球非球形引力作為最主要的攝動力,在月固系下通常采用球諧函數(shù)展開形式來表示,因此需要進行月固系與月心慣性系間的轉(zhuǎn)換,以建立環(huán)月探測器的運動方程。JPL中定義的月固坐標系原點在月球質(zhì)心,3個軸由月球的3個慣性主軸確定。本文提及的4個歷表中均已給出月球天平動,月球天平動的3個歐拉角可以用來實現(xiàn)從月心慣性坐標系到月固坐標系的轉(zhuǎn)換,具體轉(zhuǎn)換過程可參考文獻[17]。

    本文利用美國國家航空航天局(NASA)提供的重力恢復(fù)數(shù)據(jù)和內(nèi)部結(jié)構(gòu)實驗室GRAIL的Level 1B數(shù)據(jù)研究星歷對月固系坐標的影響。其中,GNI1B數(shù)據(jù)為GRAIL A/B探測器的月心慣性系坐標。選取2012-03-20的GRAIL B探測器GNI1B數(shù)據(jù),將其轉(zhuǎn)換為月固坐標系坐標,時間間隔為1 min。將基于不同行星歷表得到的月固系位置與速度矢量分別同DE440歷表轉(zhuǎn)換得到的結(jié)果作差,并將差值視為利用不同歷表進行坐標轉(zhuǎn)換時的位置和速度矢量誤差,結(jié)果見表3。由表可知,對于位置矢量,DE405歷表誤差高達30 m,均值偏離也達m級,DE421和DE430歷表誤差分別為1.3 m、1 m,兩者均值也較小,為cm級;對于速度矢量,DE405歷表誤差為3 cm/s,DE421歷表誤差為1.2 mm/s,DE430歷表誤差為0.9 mm/s。DE405在Y方向精度最差,而DE421在Y方向精度最好,甚至比DE430精度略高,但DE430歷表在Z方向精度有較大提高。綜上可知,在進行月球探測器定軌及其他涉及月球坐標系轉(zhuǎn)換任務(wù)時,當精度要求在m級時,不建議使用DE405歷表。若要求導(dǎo)航精度在亞米級別時,需要選用較新的DE430或DE440歷表。

    表3 利用不同歷表進行坐標及速度轉(zhuǎn)換時誤差最大值、均值及RMS

    4 結(jié) 語

    本文選取DE405、DE421、DE430、DE440星歷計算大行星相對于地球及太陽系質(zhì)心位置,分析各大行星的星歷位置精度。對于地心位置,由于受觀測數(shù)據(jù)等因素影響,各大行星的位置精度不同,從最高的月球m級精度,到太陽、水星、金星和火星的103m級精度,到木星、土星的105m量級,到天王星、海王星和冥王星的106m量級。總體而言,DE421和DE430星歷相對DE405在行星位置精度方面有1~2個量級的提高,DE430星歷精度相對DE421提高50%。對于太陽系質(zhì)心位置,DE440歷表由于加入30個柯伊伯帶星體的擾動,和其他3個歷表存在幾百到幾千km量級的差異。同時,針對我國正在進行的“嫦娥”系列探月工程以及其他月球科學(xué)任務(wù),比較不同歷表對于月球地心位置以及月固和月慣坐標系間轉(zhuǎn)換的影響。對于月球地心位置和速度,DE405歷表精度為7 m和0.02 mm/s,DE421歷表精度約為1.5 m和0.004 mm/s,DE430歷表精度為1.3 m和0.003 5 mm/s。利用DE405歷表進行坐標轉(zhuǎn)換時位置和速度誤差高達30 m和3 cm/s,DE421和DE430歷表位置和速度誤差分別為1.3 m、1.2 mm/s和1 m、0.9 mm/s。因此對于高精度的探月導(dǎo)航及科學(xué)研究,推薦使用高精度的DE430或DE440歷表。

    隨著我國“嫦娥”計劃的持續(xù)推進,以及地月激光測距技術(shù)的提升,將提供類型更多、精度更高的測月數(shù)據(jù),有利于我國建立自主的月球歷表,也可為我國開展月基觀測及導(dǎo)航奠定基礎(chǔ)。

    猜你喜歡
    質(zhì)心月球探測器
    重型半掛汽車質(zhì)量與質(zhì)心位置估計
    到月球上“飆車”
    軍事文摘(2022年12期)2022-07-13 03:11:44
    陪我去月球
    飛碟探索(2022年1期)2022-03-12 14:11:20
    月球上的另一個我
    軍事文摘(2021年18期)2021-12-02 01:28:54
    基于GNSS測量的天宮二號質(zhì)心確定
    第二章 探測器有反應(yīng)
    EN菌的引力波探測器
    第二章 探測器有反應(yīng)
    快拉住月球
    有7顆彗星已經(jīng)被探測器造訪過
    太空探索(2014年9期)2014-07-10 13:06:30
    www.熟女人妻精品国产| 亚洲成人手机| 建设人人有责人人尽责人人享有的| 国产日韩欧美亚洲二区| 在线观看人妻少妇| 亚洲欧美色中文字幕在线| 在线观看人妻少妇| 午夜福利视频在线观看免费| 国产日韩欧美亚洲二区| 国产深夜福利视频在线观看| 极品少妇高潮喷水抽搐| 午夜免费成人在线视频| 亚洲全国av大片| 男女午夜视频在线观看| 午夜免费成人在线视频| 亚洲国产看品久久| 日韩精品免费视频一区二区三区| 精品乱码久久久久久99久播| 国产成人精品无人区| 一本色道久久久久久精品综合| 99久久综合免费| 国产欧美日韩一区二区三区在线| 国产免费现黄频在线看| 亚洲国产av新网站| 91九色精品人成在线观看| 国产日韩欧美视频二区| 一本一本久久a久久精品综合妖精| 午夜福利在线观看吧| 成人三级做爰电影| 欧美激情久久久久久爽电影 | 俄罗斯特黄特色一大片| 爱豆传媒免费全集在线观看| 欧美乱码精品一区二区三区| videos熟女内射| 午夜福利免费观看在线| 久久久久久亚洲精品国产蜜桃av| 超碰成人久久| tube8黄色片| 老司机影院毛片| 18禁黄网站禁片午夜丰满| 人人妻人人添人人爽欧美一区卜| 男人舔女人的私密视频| 国产av精品麻豆| 极品人妻少妇av视频| 丁香六月欧美| av天堂在线播放| 人人妻人人澡人人爽人人夜夜| 日韩有码中文字幕| 女人爽到高潮嗷嗷叫在线视频| 99久久人妻综合| 性高湖久久久久久久久免费观看| av片东京热男人的天堂| 亚洲欧美一区二区三区久久| 日本欧美视频一区| 一级,二级,三级黄色视频| 亚洲欧洲日产国产| 啦啦啦啦在线视频资源| 人妻一区二区av| 青草久久国产| 色视频在线一区二区三区| 久久久久久免费高清国产稀缺| 亚洲精品一二三| 这个男人来自地球电影免费观看| av国产精品久久久久影院| 久久久久久人人人人人| 美女高潮到喷水免费观看| 欧美日韩福利视频一区二区| 国产亚洲精品一区二区www | 久久久久久久国产电影| 精品国产一区二区三区久久久樱花| 极品人妻少妇av视频| 亚洲欧美一区二区三区久久| 无遮挡黄片免费观看| 少妇裸体淫交视频免费看高清 | 999久久久国产精品视频| 老汉色av国产亚洲站长工具| 后天国语完整版免费观看| 午夜免费鲁丝| 人人妻,人人澡人人爽秒播| 国产男女超爽视频在线观看| 极品人妻少妇av视频| 91九色精品人成在线观看| tube8黄色片| 男女下面插进去视频免费观看| 成人三级做爰电影| 热re99久久精品国产66热6| 一本一本久久a久久精品综合妖精| 日本av免费视频播放| 一级毛片女人18水好多| 999久久久国产精品视频| 亚洲av电影在线进入| 精品久久久久久电影网| 国产在视频线精品| 另类亚洲欧美激情| 久久香蕉激情| 欧美 亚洲 国产 日韩一| 淫妇啪啪啪对白视频 | 亚洲激情五月婷婷啪啪| av国产精品久久久久影院| 欧美+亚洲+日韩+国产| 国产精品自产拍在线观看55亚洲 | 99香蕉大伊视频| 丁香六月欧美| 精品亚洲成国产av| 久久天躁狠狠躁夜夜2o2o| 中国国产av一级| 久久久久久免费高清国产稀缺| 日韩,欧美,国产一区二区三区| 香蕉国产在线看| 啦啦啦在线免费观看视频4| 国产在线一区二区三区精| 国内毛片毛片毛片毛片毛片| 午夜激情久久久久久久| 12—13女人毛片做爰片一| 亚洲中文av在线| 18禁黄网站禁片午夜丰满| 91九色精品人成在线观看| 国产真人三级小视频在线观看| 母亲3免费完整高清在线观看| av一本久久久久| 999久久久精品免费观看国产| 岛国毛片在线播放| 2018国产大陆天天弄谢| 少妇粗大呻吟视频| 欧美激情高清一区二区三区| 成人18禁高潮啪啪吃奶动态图| 国产精品久久久久久精品古装| 美女高潮到喷水免费观看| 两人在一起打扑克的视频| www日本在线高清视频| 欧美变态另类bdsm刘玥| 国产成人一区二区三区免费视频网站| 美女国产高潮福利片在线看| 老汉色av国产亚洲站长工具| 午夜激情久久久久久久| 中文字幕制服av| 黑丝袜美女国产一区| 首页视频小说图片口味搜索| 另类精品久久| 肉色欧美久久久久久久蜜桃| 人人澡人人妻人| 91字幕亚洲| 肉色欧美久久久久久久蜜桃| 淫妇啪啪啪对白视频 | 国产欧美日韩一区二区精品| 美国免费a级毛片| 午夜福利乱码中文字幕| 日本wwww免费看| 高清欧美精品videossex| 久久精品熟女亚洲av麻豆精品| a在线观看视频网站| 啦啦啦在线免费观看视频4| 国产野战对白在线观看| 亚洲av欧美aⅴ国产| 久久国产精品人妻蜜桃| 国产成人系列免费观看| 国产熟女午夜一区二区三区| 欧美性长视频在线观看| 欧美午夜高清在线| 国产区一区二久久| 热re99久久精品国产66热6| 人人妻人人澡人人爽人人夜夜| 国产男人的电影天堂91| 妹子高潮喷水视频| 久久人妻福利社区极品人妻图片| 青青草视频在线视频观看| 国产深夜福利视频在线观看| 亚洲人成电影观看| 色播在线永久视频| 亚洲av美国av| 国产亚洲欧美精品永久| 天天躁夜夜躁狠狠躁躁| 久久精品成人免费网站| 亚洲精品国产av成人精品| 欧美黑人精品巨大| 性少妇av在线| 丰满饥渴人妻一区二区三| 91成年电影在线观看| 亚洲精品国产一区二区精华液| 亚洲精品在线美女| 最近最新中文字幕大全免费视频| 国产日韩一区二区三区精品不卡| 啦啦啦免费观看视频1| 在线观看免费日韩欧美大片| av不卡在线播放| 久久久欧美国产精品| 一进一出抽搐动态| 亚洲欧美日韩另类电影网站| 99精国产麻豆久久婷婷| 青青草视频在线视频观看| 午夜免费成人在线视频| 香蕉国产在线看| 夜夜骑夜夜射夜夜干| 久久中文看片网| 精品国产一区二区久久| 欧美日韩av久久| 一二三四社区在线视频社区8| 日本撒尿小便嘘嘘汇集6| 欧美精品亚洲一区二区| 亚洲中文av在线| 亚洲精品国产一区二区精华液| 黄色 视频免费看| 少妇裸体淫交视频免费看高清 | 青春草视频在线免费观看| 中文欧美无线码| 精品久久久精品久久久| 精品人妻熟女毛片av久久网站| 老司机午夜福利在线观看视频 | 十八禁高潮呻吟视频| 亚洲激情五月婷婷啪啪| 超碰成人久久| 国产一区二区三区在线臀色熟女 | 黄色视频不卡| 大片电影免费在线观看免费| av天堂在线播放| 免费高清在线观看视频在线观看| 老司机亚洲免费影院| 亚洲成av片中文字幕在线观看| 女人精品久久久久毛片| 18禁观看日本| 一区二区三区精品91| 午夜福利影视在线免费观看| 日本精品一区二区三区蜜桃| 丰满人妻熟妇乱又伦精品不卡| 久久人人97超碰香蕉20202| 亚洲国产欧美一区二区综合| 高清欧美精品videossex| 91精品三级在线观看| 一本色道久久久久久精品综合| 国产无遮挡羞羞视频在线观看| 亚洲七黄色美女视频| 久久99热这里只频精品6学生| 免费少妇av软件| 国产熟女午夜一区二区三区| 国产日韩欧美视频二区| 美女大奶头黄色视频| 中文字幕高清在线视频| 男女床上黄色一级片免费看| 国产精品欧美亚洲77777| 国产在线一区二区三区精| 欧美亚洲 丝袜 人妻 在线| 国产成+人综合+亚洲专区| 午夜福利,免费看| 黑人欧美特级aaaaaa片| 日韩精品免费视频一区二区三区| 欧美黄色淫秽网站| 国产深夜福利视频在线观看| 久久久国产成人免费| 精品福利观看| 天堂中文最新版在线下载| 久久久久久久精品精品| 欧美老熟妇乱子伦牲交| svipshipincom国产片| 极品少妇高潮喷水抽搐| 人人妻人人澡人人看| 亚洲欧美激情在线| 亚洲av欧美aⅴ国产| 9191精品国产免费久久| 国产精品九九99| 中文精品一卡2卡3卡4更新| 悠悠久久av| 天天操日日干夜夜撸| 久久国产亚洲av麻豆专区| 久久精品aⅴ一区二区三区四区| 热99re8久久精品国产| 亚洲成人免费电影在线观看| av超薄肉色丝袜交足视频| 亚洲国产欧美网| 欧美黄色片欧美黄色片| 国产亚洲欧美精品永久| 精品亚洲成a人片在线观看| 久久九九热精品免费| 热99国产精品久久久久久7| 国产精品.久久久| 欧美中文综合在线视频| 免费少妇av软件| 午夜老司机福利片| 欧美 日韩 精品 国产| 日韩人妻精品一区2区三区| 日韩免费高清中文字幕av| 亚洲精品美女久久久久99蜜臀| 大片电影免费在线观看免费| 欧美精品高潮呻吟av久久| 国产极品粉嫩免费观看在线| 老司机靠b影院| 中国美女看黄片| 欧美人与性动交α欧美软件| 国产一区二区三区在线臀色熟女 | 中文字幕另类日韩欧美亚洲嫩草| 精品卡一卡二卡四卡免费| 高清视频免费观看一区二区| 青春草亚洲视频在线观看| 日韩精品免费视频一区二区三区| 国产男人的电影天堂91| 在线观看人妻少妇| 国产成人免费观看mmmm| 女人高潮潮喷娇喘18禁视频| bbb黄色大片| 建设人人有责人人尽责人人享有的| 天天操日日干夜夜撸| av电影中文网址| 国产亚洲欧美在线一区二区| 脱女人内裤的视频| 免费在线观看完整版高清| 国产精品欧美亚洲77777| 午夜影院在线不卡| 久久久久国产一级毛片高清牌| 好男人电影高清在线观看| 电影成人av| 亚洲精品久久午夜乱码| 亚洲精品一区蜜桃| 别揉我奶头~嗯~啊~动态视频 | 日本一区二区免费在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲专区中文字幕在线| 在线观看免费午夜福利视频| 三上悠亚av全集在线观看| 一区二区三区精品91| 一二三四在线观看免费中文在| 高清黄色对白视频在线免费看| netflix在线观看网站| 国产三级黄色录像| 久久 成人 亚洲| 深夜精品福利| 黑人欧美特级aaaaaa片| 亚洲成av片中文字幕在线观看| 搡老乐熟女国产| h视频一区二区三区| 一本—道久久a久久精品蜜桃钙片| 国产深夜福利视频在线观看| 精品国产乱码久久久久久小说| 亚洲精品在线美女| 成人黄色视频免费在线看| 亚洲熟女毛片儿| 一级,二级,三级黄色视频| 一边摸一边抽搐一进一出视频| 男女之事视频高清在线观看| 中文字幕精品免费在线观看视频| 日韩欧美免费精品| 免费日韩欧美在线观看| 亚洲欧美一区二区三区黑人| 五月天丁香电影| 精品免费久久久久久久清纯 | 在线观看免费高清a一片| 国产在线视频一区二区| 啦啦啦视频在线资源免费观看| 国产精品国产av在线观看| 69精品国产乱码久久久| 亚洲成人免费电影在线观看| 欧美成狂野欧美在线观看| 99九九在线精品视频| 久久精品国产亚洲av香蕉五月 | e午夜精品久久久久久久| www.自偷自拍.com| 性少妇av在线| 少妇精品久久久久久久| 欧美一级毛片孕妇| 亚洲国产看品久久| 一级毛片电影观看| 黄片小视频在线播放| 正在播放国产对白刺激| 欧美黄色片欧美黄色片| 久久青草综合色| 十八禁网站网址无遮挡| 久久国产亚洲av麻豆专区| 国产三级黄色录像| 9191精品国产免费久久| 免费观看人在逋| 高清欧美精品videossex| 91av网站免费观看| 女性被躁到高潮视频| 最近最新中文字幕大全免费视频| 丝袜脚勾引网站| 老鸭窝网址在线观看| 国产精品熟女久久久久浪| 亚洲七黄色美女视频| 岛国在线观看网站| 最近最新免费中文字幕在线| 建设人人有责人人尽责人人享有的| 桃红色精品国产亚洲av| 不卡一级毛片| 中文字幕人妻熟女乱码| 亚洲伊人色综图| 性少妇av在线| 搡老岳熟女国产| 老司机福利观看| 巨乳人妻的诱惑在线观看| 国产亚洲欧美精品永久| 久久天堂一区二区三区四区| 久久综合国产亚洲精品| 国产三级黄色录像| av电影中文网址| 天堂中文最新版在线下载| 老司机影院成人| 人成视频在线观看免费观看| 国产精品1区2区在线观看. | 亚洲精品美女久久av网站| 国产av国产精品国产| 日韩熟女老妇一区二区性免费视频| 黄色毛片三级朝国网站| 老司机影院毛片| 日韩免费高清中文字幕av| 精品免费久久久久久久清纯 | 亚洲精品国产av成人精品| 久久人妻熟女aⅴ| 久久久久网色| 汤姆久久久久久久影院中文字幕| 蜜桃国产av成人99| 国产成人一区二区三区免费视频网站| 丝袜喷水一区| 欧美国产精品一级二级三级| 亚洲,欧美精品.| 色视频在线一区二区三区| 岛国毛片在线播放| 亚洲专区字幕在线| 精品福利永久在线观看| 国产视频一区二区在线看| 欧美成狂野欧美在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 免费一级毛片在线播放高清视频 | 亚洲免费av在线视频| 国产日韩欧美在线精品| 国产高清videossex| 亚洲成人国产一区在线观看| 中国国产av一级| 老汉色av国产亚洲站长工具| 黑人操中国人逼视频| 国产伦人伦偷精品视频| 国产精品1区2区在线观看. | 五月开心婷婷网| 欧美乱码精品一区二区三区| 交换朋友夫妻互换小说| 欧美日韩精品网址| 欧美精品一区二区免费开放| 自线自在国产av| 十分钟在线观看高清视频www| 日本欧美视频一区| 青春草亚洲视频在线观看| 欧美成人午夜精品| 亚洲美女黄色视频免费看| 亚洲精品国产av蜜桃| 两人在一起打扑克的视频| 精品一区二区三卡| 欧美久久黑人一区二区| 90打野战视频偷拍视频| 黄色a级毛片大全视频| 国产又色又爽无遮挡免| 精品久久蜜臀av无| 久久精品国产a三级三级三级| 成人黄色视频免费在线看| 亚洲五月色婷婷综合| 亚洲欧美一区二区三区黑人| 黑人操中国人逼视频| 国产精品.久久久| 精品国产乱子伦一区二区三区 | 伦理电影免费视频| 日韩电影二区| 可以免费在线观看a视频的电影网站| 亚洲精品一区蜜桃| 免费在线观看黄色视频的| 婷婷成人精品国产| 十八禁人妻一区二区| 97人妻天天添夜夜摸| 亚洲成av片中文字幕在线观看| 黄片播放在线免费| 国产成人精品在线电影| 丝袜脚勾引网站| 欧美成狂野欧美在线观看| 乱人伦中国视频| 国产亚洲欧美精品永久| 亚洲久久久国产精品| netflix在线观看网站| 天天操日日干夜夜撸| 亚洲精品在线美女| 两人在一起打扑克的视频| av欧美777| 黄色怎么调成土黄色| 俄罗斯特黄特色一大片| 久热爱精品视频在线9| 日本av免费视频播放| 日本五十路高清| 欧美另类亚洲清纯唯美| 日本精品一区二区三区蜜桃| av电影中文网址| 中国美女看黄片| 国产欧美日韩一区二区三区在线| 欧美一级毛片孕妇| 中文字幕人妻丝袜一区二区| 国产成+人综合+亚洲专区| 亚洲一区二区三区欧美精品| 午夜精品国产一区二区电影| 久久久国产一区二区| 日本精品一区二区三区蜜桃| 午夜福利影视在线免费观看| 日韩大码丰满熟妇| 国产在线一区二区三区精| www.精华液| 欧美日本中文国产一区发布| 中文字幕精品免费在线观看视频| 嫁个100分男人电影在线观看| 黑人巨大精品欧美一区二区蜜桃| 美女扒开内裤让男人捅视频| 国产一区二区三区综合在线观看| 久久精品亚洲熟妇少妇任你| 国产一区二区三区av在线| 精品亚洲成a人片在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产三级黄色录像| 9热在线视频观看99| 亚洲美女黄色视频免费看| 新久久久久国产一级毛片| 伦理电影免费视频| 久久国产精品大桥未久av| a 毛片基地| 精品熟女少妇八av免费久了| 国产淫语在线视频| 亚洲 欧美一区二区三区| 成年美女黄网站色视频大全免费| 999久久久精品免费观看国产| 亚洲免费av在线视频| 1024香蕉在线观看| 国产免费一区二区三区四区乱码| 亚洲一区二区三区欧美精品| 亚洲精品一区蜜桃| 搡老乐熟女国产| 天天添夜夜摸| 少妇 在线观看| 亚洲国产av影院在线观看| 日韩欧美一区二区三区在线观看 | 国产精品秋霞免费鲁丝片| 精品视频人人做人人爽| 免费一级毛片在线播放高清视频 | 欧美另类一区| 亚洲人成77777在线视频| 久久国产精品影院| 久久久精品国产亚洲av高清涩受| 国产精品久久久人人做人人爽| 国产免费视频播放在线视频| 欧美在线一区亚洲| 91字幕亚洲| 国产激情久久老熟女| 成人免费观看视频高清| 久久久久国产精品人妻一区二区| 视频区图区小说| 国产av一区二区精品久久| 午夜日韩欧美国产| 国产黄色免费在线视频| 嫩草影视91久久| 久久久精品94久久精品| 成在线人永久免费视频| 天堂俺去俺来也www色官网| 一本色道久久久久久精品综合| 深夜精品福利| 五月天丁香电影| 欧美日韩福利视频一区二区| 久久精品国产亚洲av香蕉五月 | 久久青草综合色| 99国产精品一区二区蜜桃av | 国产精品一二三区在线看| 各种免费的搞黄视频| 十八禁网站免费在线| 欧美激情久久久久久爽电影 | 国产免费一区二区三区四区乱码| 汤姆久久久久久久影院中文字幕| 亚洲成人免费av在线播放| 50天的宝宝边吃奶边哭怎么回事| 日本a在线网址| 纯流量卡能插随身wifi吗| 国产亚洲午夜精品一区二区久久| 一个人免费在线观看的高清视频 | 国产成人精品久久二区二区91| 亚洲国产成人一精品久久久| 国产激情久久老熟女| 亚洲av电影在线进入| 国产亚洲一区二区精品| av又黄又爽大尺度在线免费看| 国产欧美亚洲国产| 亚洲专区国产一区二区| 久久久久久亚洲精品国产蜜桃av| 十八禁网站免费在线| 9热在线视频观看99| 日本vs欧美在线观看视频| 欧美黑人欧美精品刺激| 男人爽女人下面视频在线观看| 热99re8久久精品国产| 欧美老熟妇乱子伦牲交| 久久精品亚洲熟妇少妇任你| 99热国产这里只有精品6| 久久天堂一区二区三区四区| 欧美 亚洲 国产 日韩一| 亚洲人成电影免费在线| 热99久久久久精品小说推荐| 国产成+人综合+亚洲专区| 久久香蕉激情| 国产又色又爽无遮挡免| 超碰成人久久| 国产成人av教育| 国产91精品成人一区二区三区 | 久久精品熟女亚洲av麻豆精品| 大片免费播放器 马上看| 中文欧美无线码| 精品一品国产午夜福利视频| 午夜日韩欧美国产| 日本av免费视频播放| 日韩欧美一区视频在线观看| 国产免费一区二区三区四区乱码| 亚洲美女黄色视频免费看| 亚洲精品久久成人aⅴ小说| 咕卡用的链子| 亚洲色图 男人天堂 中文字幕| 国产淫语在线视频| 欧美乱码精品一区二区三区| 午夜精品久久久久久毛片777| 国产精品久久久人人做人人爽| 免费在线观看视频国产中文字幕亚洲 |