• <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
    日韩不卡一区二区三区视频在线| 精品酒店卫生间| 高清午夜精品一区二区三区| 哪个播放器可以免费观看大片| 久久久国产一区二区| 久久精品久久精品一区二区三区| 狂野欧美白嫩少妇大欣赏| 观看免费一级毛片| 全区人妻精品视频| 草草在线视频免费看| 美女xxoo啪啪120秒动态图| 国产在视频线精品| 超碰97精品在线观看| 国产午夜精品一二区理论片| 91久久精品电影网| 麻豆国产97在线/欧美| 春色校园在线视频观看| 亚洲av成人精品一区久久| 久久女婷五月综合色啪小说 | 精品一区二区免费观看| av国产免费在线观看| 又黄又爽又刺激的免费视频.| 久久精品久久久久久久性| 在线免费十八禁| 久久午夜福利片| 精品久久久精品久久久| 天天躁日日操中文字幕| 成人特级av手机在线观看| 深爱激情五月婷婷| 亚州av有码| 视频区图区小说| 99热国产这里只有精品6| 国产男人的电影天堂91| 久久久精品94久久精品| 美女xxoo啪啪120秒动态图| 午夜福利视频1000在线观看| 天美传媒精品一区二区| 搡老乐熟女国产| 丰满乱子伦码专区| 亚洲欧美日韩卡通动漫| 听说在线观看完整版免费高清| 国产亚洲av嫩草精品影院| 亚洲精品成人av观看孕妇| 国产精品蜜桃在线观看| 国产欧美日韩精品一区二区| 九草在线视频观看| 国产成人午夜福利电影在线观看| 黄色欧美视频在线观看| 最后的刺客免费高清国语| 日韩一区二区三区影片| 亚洲天堂国产精品一区在线| 免费大片黄手机在线观看| 天天躁夜夜躁狠狠久久av| 日本黄大片高清| 亚洲国产日韩一区二区| 不卡视频在线观看欧美| 国产精品99久久久久久久久| 麻豆成人av视频| 国产黄片美女视频| 成人亚洲精品一区在线观看 | 精品久久国产蜜桃| 小蜜桃在线观看免费完整版高清| 在线观看免费高清a一片| 性插视频无遮挡在线免费观看| 国产精品熟女久久久久浪| 亚洲综合精品二区| 成人鲁丝片一二三区免费| 麻豆乱淫一区二区| 久久精品国产亚洲av天美| 亚洲av不卡在线观看| 别揉我奶头 嗯啊视频| 免费少妇av软件| 日韩成人伦理影院| 国产白丝娇喘喷水9色精品| 高清午夜精品一区二区三区| 纵有疾风起免费观看全集完整版| 亚洲欧美日韩无卡精品| 久久久成人免费电影| 男女边吃奶边做爰视频| 国产精品女同一区二区软件| 一本久久精品| 亚洲综合精品二区| 国产探花在线观看一区二区| 久久精品久久久久久久性| 日韩亚洲欧美综合| 能在线免费看毛片的网站| 国产成人aa在线观看| 又爽又黄无遮挡网站| 一级爰片在线观看| 内射极品少妇av片p| 不卡视频在线观看欧美| 身体一侧抽搐| 国产老妇女一区| 欧美亚洲 丝袜 人妻 在线| 亚洲成人中文字幕在线播放| 少妇人妻久久综合中文| 久久综合国产亚洲精品| 国语对白做爰xxxⅹ性视频网站| 美女cb高潮喷水在线观看| 久久久a久久爽久久v久久| 在线天堂最新版资源| 18禁裸乳无遮挡动漫免费视频 | av天堂中文字幕网| 男女边吃奶边做爰视频| 免费看av在线观看网站| 在线a可以看的网站| 我的女老师完整版在线观看| 韩国av在线不卡| 人人妻人人看人人澡| 日日啪夜夜撸| 在线观看美女被高潮喷水网站| 国产精品一二三区在线看| 亚洲精华国产精华液的使用体验| 成人黄色视频免费在线看| 日韩成人av中文字幕在线观看| 特大巨黑吊av在线直播| 欧美区成人在线视频| 国产精品秋霞免费鲁丝片| 精品久久久精品久久久| 97超碰精品成人国产| 欧美亚洲 丝袜 人妻 在线| 国产黄片视频在线免费观看| 国产爽快片一区二区三区| 精品国产露脸久久av麻豆| 久久久国产一区二区| 亚洲欧美日韩另类电影网站 | 久久久久久久久久久免费av| 国产男女内射视频| 女人被狂操c到高潮| 国产黄片美女视频| 欧美另类一区| 美女视频免费永久观看网站| 欧美最新免费一区二区三区| 国产精品国产三级国产av玫瑰| 插阴视频在线观看视频| 可以在线观看毛片的网站| 亚洲精品日韩在线中文字幕| 免费观看性生交大片5| 男人狂女人下面高潮的视频| 大话2 男鬼变身卡| 色视频在线一区二区三区| 欧美日韩亚洲高清精品| 国产成人精品久久久久久| 一级毛片 在线播放| 国产精品99久久99久久久不卡 | 亚洲最大成人手机在线| 日韩,欧美,国产一区二区三区| 99久久精品热视频| 国产淫片久久久久久久久| 少妇的逼水好多| 亚洲成人中文字幕在线播放| 亚洲国产精品成人久久小说| 另类亚洲欧美激情| 欧美bdsm另类| 亚洲精品日本国产第一区| 亚洲精品一二三| 日韩精品有码人妻一区| av女优亚洲男人天堂| 亚洲精华国产精华液的使用体验| 国产一区二区在线观看日韩| 久久久久精品性色| 人人妻人人澡人人爽人人夜夜| 国产成人精品婷婷| 蜜桃久久精品国产亚洲av| 三级国产精品欧美在线观看| 99九九线精品视频在线观看视频| 舔av片在线| 女的被弄到高潮叫床怎么办| 国内精品美女久久久久久| 日韩制服骚丝袜av| 深爱激情五月婷婷| 精品久久久久久久末码| 久久久久网色| 伊人久久国产一区二区| 欧美激情久久久久久爽电影| 91aial.com中文字幕在线观看| a级毛片免费高清观看在线播放| 国产精品熟女久久久久浪| 男女那种视频在线观看| 国产精品久久久久久久电影| 日韩不卡一区二区三区视频在线| 日日摸夜夜添夜夜爱| 国产精品久久久久久精品古装| 亚洲精华国产精华液的使用体验| 久久久久久久久久人人人人人人| 日韩欧美精品v在线| 国产成年人精品一区二区| 久久综合国产亚洲精品| 久久人人爽人人片av| 国产精品麻豆人妻色哟哟久久| 日韩国内少妇激情av| 亚洲精品,欧美精品| 成人特级av手机在线观看| 国产精品久久久久久精品古装| 国产成人一区二区在线| 国产黄片美女视频| 精品久久久久久久久亚洲| 成人漫画全彩无遮挡| 日韩免费高清中文字幕av| 国产精品久久久久久精品古装| 亚洲国产精品999| 2022亚洲国产成人精品| 国产黄a三级三级三级人| av卡一久久| 九九爱精品视频在线观看| 国产成人精品福利久久| 边亲边吃奶的免费视频| 久久99热这里只有精品18| 精品久久国产蜜桃| 一区二区三区乱码不卡18| 在现免费观看毛片| 久久精品久久久久久久性| 老女人水多毛片| 欧美日本视频| 国产黄片美女视频| 少妇人妻一区二区三区视频| 久久ye,这里只有精品| 美女cb高潮喷水在线观看| av黄色大香蕉| 一个人看的www免费观看视频| 日韩中字成人| 黄色怎么调成土黄色| 日韩一区二区三区影片| 久久久欧美国产精品| 自拍偷自拍亚洲精品老妇| 午夜激情久久久久久久| 王馨瑶露胸无遮挡在线观看| 中文字幕免费在线视频6| 国产毛片在线视频| 久久久午夜欧美精品| 欧美激情在线99| 精品一区二区三卡| 美女主播在线视频| 午夜爱爱视频在线播放| 精品国产一区二区三区久久久樱花 | 欧美潮喷喷水| 一级毛片 在线播放| 啦啦啦啦在线视频资源| 中文天堂在线官网| 18禁裸乳无遮挡免费网站照片| 人人妻人人看人人澡| 日日撸夜夜添| 九九爱精品视频在线观看| freevideosex欧美| 亚洲精品日本国产第一区| 99久久九九国产精品国产免费| 一级毛片电影观看| 国产日韩欧美在线精品| 日韩不卡一区二区三区视频在线| 国产精品久久久久久久久免| 亚洲一区二区三区欧美精品 | 亚洲精品视频女| 国产在视频线精品| av天堂中文字幕网| 免费观看性生交大片5| 久久久精品免费免费高清| 久久久久久久久大av| 国产视频内射| 国产永久视频网站| 男人爽女人下面视频在线观看| 天天一区二区日本电影三级| 嫩草影院精品99| 亚洲天堂国产精品一区在线| 精品国产乱码久久久久久小说| 亚洲欧美日韩另类电影网站 | av在线观看视频网站免费| 午夜老司机福利剧场| 亚洲性久久影院| 2018国产大陆天天弄谢| 日本一本二区三区精品| 七月丁香在线播放| 国产在线一区二区三区精| 精品人妻视频免费看| 久久久久九九精品影院| 国产又色又爽无遮挡免| 1000部很黄的大片| 国产伦精品一区二区三区视频9| 亚洲高清免费不卡视频| 性色av一级| 99热全是精品| 久热这里只有精品99| 亚洲精品日韩在线中文字幕| 在线观看国产h片| 久久久久性生活片| 特级一级黄色大片| 精品久久久久久久久亚洲| 国产在视频线精品| 久久精品国产a三级三级三级| 亚洲成人一二三区av| 成人免费观看视频高清| 精华霜和精华液先用哪个| 91精品伊人久久大香线蕉| 日韩欧美一区视频在线观看 | 九色成人免费人妻av| 免费黄色在线免费观看| 成人高潮视频无遮挡免费网站| av国产免费在线观看| 午夜亚洲福利在线播放| 九九爱精品视频在线观看| 人体艺术视频欧美日本| 777米奇影视久久| 最近的中文字幕免费完整| 亚洲国产精品999| 一区二区av电影网| 老司机影院成人| 国产精品国产av在线观看| 欧美bdsm另类| 大香蕉97超碰在线| 欧美老熟妇乱子伦牲交| 国产伦精品一区二区三区视频9| 亚洲国产精品成人久久小说| 欧美激情久久久久久爽电影| 亚洲内射少妇av| 又爽又黄无遮挡网站| 国产黄片美女视频| 高清av免费在线| 一区二区三区免费毛片| 日本一本二区三区精品| 久久久欧美国产精品| av网站免费在线观看视频| 一级毛片久久久久久久久女| 国产乱人偷精品视频| 你懂的网址亚洲精品在线观看| 高清av免费在线| 一级av片app| av福利片在线观看| 又粗又硬又长又爽又黄的视频| 亚洲成人一二三区av| 亚洲第一区二区三区不卡| 九色成人免费人妻av| 久热久热在线精品观看| 国产av码专区亚洲av| 日韩强制内射视频| 亚洲精品456在线播放app| 欧美+日韩+精品| 草草在线视频免费看| 国产成人aa在线观看| 大香蕉久久网| 99热网站在线观看| 下体分泌物呈黄色| 亚洲,欧美,日韩| 国产精品av视频在线免费观看| 好男人在线观看高清免费视频| 色综合色国产| 九色成人免费人妻av| 欧美成人午夜免费资源| 高清午夜精品一区二区三区| 日韩成人av中文字幕在线观看| 午夜免费观看性视频| av免费观看日本| 成人无遮挡网站| 欧美老熟妇乱子伦牲交| 国产伦在线观看视频一区| 丝袜喷水一区| av.在线天堂| 黄色一级大片看看| 2022亚洲国产成人精品| 久久人人爽av亚洲精品天堂 | 亚洲美女视频黄频| 国产男女超爽视频在线观看| 亚洲av日韩在线播放| 色视频在线一区二区三区| 久久精品国产a三级三级三级| 国产免费一级a男人的天堂| 中文精品一卡2卡3卡4更新| 一本久久精品| 天天一区二区日本电影三级| 国产成人精品婷婷| 极品教师在线视频| 欧美区成人在线视频| 听说在线观看完整版免费高清| 中文在线观看免费www的网站| 国产亚洲5aaaaa淫片| 99久久精品国产国产毛片| 国产探花极品一区二区| 久久人人爽人人爽人人片va| 中文字幕亚洲精品专区| 亚洲欧美精品自产自拍| 久久久久久久精品精品| 国产综合懂色| 国产精品福利在线免费观看| 男插女下体视频免费在线播放| 伦精品一区二区三区| 中国三级夫妇交换| 天堂俺去俺来也www色官网| 亚洲精品自拍成人| 韩国av在线不卡| 免费看a级黄色片| 美女高潮的动态| 亚洲精品,欧美精品| 免费看av在线观看网站| 久久久久精品久久久久真实原创| 极品少妇高潮喷水抽搐| 啦啦啦在线观看免费高清www| 国产免费又黄又爽又色| 色哟哟·www| 亚洲人成网站在线观看播放| 欧美性猛交╳xxx乱大交人| 久久久久久久久久成人| 午夜福利视频精品| 欧美日韩一区二区视频在线观看视频在线 | 国产亚洲av嫩草精品影院| 身体一侧抽搐| 免费黄色在线免费观看| 少妇人妻一区二区三区视频| 九九爱精品视频在线观看| 少妇人妻久久综合中文| 精品久久久噜噜| 一级a做视频免费观看| 秋霞伦理黄片| 一个人看视频在线观看www免费| 最近最新中文字幕免费大全7| 国产亚洲av片在线观看秒播厂| 大码成人一级视频| 人妻 亚洲 视频| 男女国产视频网站| 高清视频免费观看一区二区| 男男h啪啪无遮挡| 99视频精品全部免费 在线| 好男人视频免费观看在线| 久久久久久久久久久免费av| 2021少妇久久久久久久久久久| av.在线天堂| 简卡轻食公司| 国产精品国产三级专区第一集| 男人狂女人下面高潮的视频| 久久久成人免费电影| 国产高清三级在线| 精品一区在线观看国产| 成人高潮视频无遮挡免费网站| 亚洲婷婷狠狠爱综合网| av在线蜜桃| 大香蕉97超碰在线| 亚洲一区二区三区欧美精品 | 久久久久久伊人网av| 亚洲在线观看片| 久久久久久久午夜电影| 伦精品一区二区三区| 国产免费视频播放在线视频| 亚洲最大成人中文| 青青草视频在线视频观看| 亚洲国产成人一精品久久久| 亚洲aⅴ乱码一区二区在线播放| 可以在线观看毛片的网站| 国产成人免费观看mmmm| 丝瓜视频免费看黄片| 国产乱人偷精品视频| 只有这里有精品99| 人妻系列 视频| 自拍偷自拍亚洲精品老妇| 国产淫片久久久久久久久| 美女脱内裤让男人舔精品视频| 欧美zozozo另类| 搞女人的毛片| 噜噜噜噜噜久久久久久91| 中文在线观看免费www的网站| 肉色欧美久久久久久久蜜桃 | 国产成人a区在线观看| 成人国产麻豆网| 欧美高清性xxxxhd video| 日本色播在线视频| 99久久精品国产国产毛片| 人妻制服诱惑在线中文字幕| 日韩成人av中文字幕在线观看| 日韩国内少妇激情av| 久久精品久久精品一区二区三区| 国产伦理片在线播放av一区| 我要看日韩黄色一级片| 欧美变态另类bdsm刘玥| 国产一区亚洲一区在线观看| 日本一本二区三区精品| 又粗又硬又长又爽又黄的视频| 国产色婷婷99| 日韩av免费高清视频| 日韩视频在线欧美| 七月丁香在线播放| 免费黄网站久久成人精品| 亚洲电影在线观看av| 有码 亚洲区| 亚洲三级黄色毛片| 国产精品无大码| 啦啦啦啦在线视频资源| 婷婷色麻豆天堂久久| 2021少妇久久久久久久久久久| 在线a可以看的网站| 国产av国产精品国产| 日日摸夜夜添夜夜添av毛片| 夫妻午夜视频| 看黄色毛片网站| 国产男女内射视频| 国产欧美日韩精品一区二区| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲色图综合在线观看| 亚洲在久久综合| 中国国产av一级| 欧美激情国产日韩精品一区| 久久久精品94久久精品| 国产精品国产av在线观看| 在线观看国产h片| 国产在线男女| 亚洲国产欧美人成| 日韩一区二区视频免费看| 少妇丰满av| 啦啦啦啦在线视频资源| 在线a可以看的网站| 国产视频首页在线观看| 人妻 亚洲 视频| 一级片'在线观看视频| 晚上一个人看的免费电影| 国产 一区精品| 美女高潮的动态| 国产精品成人在线| 天堂俺去俺来也www色官网| 久久久午夜欧美精品| 亚洲成人久久爱视频| 中文天堂在线官网| 2021少妇久久久久久久久久久| 大陆偷拍与自拍| av一本久久久久| 可以在线观看毛片的网站| 麻豆国产97在线/欧美| 精品一区二区三卡| 亚洲自偷自拍三级| av在线老鸭窝| 美女国产视频在线观看| 欧美极品一区二区三区四区| 3wmmmm亚洲av在线观看| 亚洲精品456在线播放app| 2021天堂中文幕一二区在线观| 国产老妇伦熟女老妇高清| 大话2 男鬼变身卡| 热99国产精品久久久久久7| 亚洲国产高清在线一区二区三| 国语对白做爰xxxⅹ性视频网站| 人人妻人人澡人人爽人人夜夜| 少妇猛男粗大的猛烈进出视频 | 激情五月婷婷亚洲| 男女无遮挡免费网站观看| 久久精品国产自在天天线| 欧美 日韩 精品 国产| 在线观看av片永久免费下载| 国产精品嫩草影院av在线观看| 亚洲av欧美aⅴ国产| 国产精品99久久99久久久不卡 | www.色视频.com| 乱码一卡2卡4卡精品| 少妇人妻 视频| 亚洲av在线观看美女高潮| 一级毛片 在线播放| 日韩电影二区| 久久ye,这里只有精品| 欧美+日韩+精品| 久久6这里有精品| 白带黄色成豆腐渣| 性插视频无遮挡在线免费观看| 精品国产露脸久久av麻豆| 日本一二三区视频观看| 亚洲av日韩在线播放| 麻豆成人午夜福利视频| 男女边摸边吃奶| 美女高潮的动态| 五月伊人婷婷丁香| 九九爱精品视频在线观看| 搡女人真爽免费视频火全软件| 中文字幕制服av| 国内揄拍国产精品人妻在线| 国产v大片淫在线免费观看| 欧美bdsm另类| 亚洲图色成人| 美女主播在线视频| 高清毛片免费看| 精品一区二区免费观看| 国产日韩欧美亚洲二区| 黑人高潮一二区| av专区在线播放| 日韩av免费高清视频| 久久久欧美国产精品| 偷拍熟女少妇极品色| 久久精品国产亚洲av涩爱| 男男h啪啪无遮挡| 欧美日韩一区二区视频在线观看视频在线 | 天天躁日日操中文字幕| 国产av码专区亚洲av| 国产淫片久久久久久久久| 亚洲一区二区三区欧美精品 | 久久99蜜桃精品久久| 欧美日韩亚洲高清精品| 国产精品秋霞免费鲁丝片| videos熟女内射| 成人国产麻豆网| 久久久a久久爽久久v久久| 亚洲国产欧美在线一区| 美女被艹到高潮喷水动态| 亚洲精品国产色婷婷电影| 内射极品少妇av片p| 性插视频无遮挡在线免费观看| 国产免费一级a男人的天堂| 久久精品国产亚洲网站| 大片免费播放器 马上看| 欧美最新免费一区二区三区| 最近2019中文字幕mv第一页| 国产精品成人在线| av黄色大香蕉| 国产亚洲午夜精品一区二区久久 | 99久久精品热视频| 日日啪夜夜爽| 国产成人精品婷婷| 亚洲,欧美,日韩| 一区二区三区免费毛片| 成人毛片a级毛片在线播放| 最近2019中文字幕mv第一页| 男女边摸边吃奶| 麻豆精品久久久久久蜜桃| 亚洲欧美精品专区久久| 国产一级毛片在线|