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

    基于協(xié)方差分析研究VGOS不同頻段組合對(duì)群時(shí)延精度的影響

    2022-01-06 01:59:24黃逸丹舒逢春鄭為民
    天文學(xué)進(jìn)展 2021年4期
    關(guān)鍵詞:群時(shí)延基線頻段

    黃逸丹,舒逢春,吳 德,鄭為民

    (1.中國科學(xué)院 上海天文臺(tái),上海200030; 2.中國科學(xué)院大學(xué),北京100049; 3.國家基礎(chǔ)學(xué)科公共科學(xué)數(shù)據(jù)中心,北京100190; 4.中國科學(xué)院 射電天文重點(diǎn)實(shí)驗(yàn)室,南京210008; 5.上海市導(dǎo)航定位重點(diǎn)實(shí)驗(yàn)室,上海200030)

    1 引言

    20世紀(jì)60年代末,射電天文學(xué)家提出甚長基線干涉測(cè)量(very long baseline interferometry,VLBI)技術(shù)[1]。經(jīng)過數(shù)十年的發(fā)展,VLBI臺(tái)站定位精度從米量級(jí)提高到亞厘米水平[2]。VLBI精度的提高主要?dú)w功于帶寬綜合技術(shù)[3]、用于儀器時(shí)延改正的校正系統(tǒng)[4]和用于消除電離層影響的雙頻觀測(cè)策略等。

    隨著電子設(shè)備、全球光纖網(wǎng)絡(luò)和低成本天線等領(lǐng)域的技術(shù)發(fā)展,21世紀(jì)初國際VLBI大地測(cè)量與天體測(cè)量服務(wù)組織(International VLBI Service for Geodesy and Astrometry,IVS[5,6])提出了下一代測(cè)地VLBI系統(tǒng)的構(gòu)想[7],2012年正式命名為VLBI全球觀測(cè)系統(tǒng)(VLBI Global Observing System,VGOS),以適應(yīng)未來天體測(cè)量與大地測(cè)量的高精度需求。

    2003年,IVS成立了專門的工作組WG3對(duì)下一代測(cè)地VLBI系統(tǒng)的需求進(jìn)行調(diào)研[2],2005年9月,WG3提出VGOS的目標(biāo)為:(1)全球尺度實(shí)現(xiàn)臺(tái)站位置和EOP連續(xù)監(jiān)測(cè);(2)24 h快速數(shù)據(jù)處理;(3)臺(tái)站位置精度達(dá)到1 mm[8]。為了實(shí)現(xiàn)上述目標(biāo),VGOS采用一系列措施減少時(shí)延觀測(cè)量的系統(tǒng)誤差和隨機(jī)誤差影響。為了減少隨機(jī)誤差,VGOS采用口徑約為12 m的天線,以縮短單個(gè)射電源的觀測(cè)時(shí)間,增加參考源的全天區(qū)覆蓋率;在2~14 GHz頻率范圍內(nèi)可自由選擇最多4個(gè)頻段開展觀測(cè),以增加觀測(cè)的有效帶寬。新的寬帶群時(shí)延精度預(yù)期為4 ps[7,9]。

    測(cè)地VLBI觀測(cè)中,群時(shí)延精度與有效帶寬和信噪比有關(guān),而信噪比取決于通道數(shù)、單通道帶寬、觀測(cè)時(shí)長、臺(tái)站的系統(tǒng)等效流量密度(SEFD)以及射電源的相關(guān)流量密度等因素。對(duì)于特定的基線、目標(biāo)射電源和觀測(cè)時(shí)長來說,頻率設(shè)置是影響群時(shí)延觀測(cè)精度的主要因素。由于超寬頻帶的使用,還需解決色散效應(yīng)帶來的電離層影響,在得到群時(shí)延觀測(cè)量的同時(shí)擬合基線的差分總電子含量(dT E C)[10,11]。寬帶群時(shí)延與dT E C具有很強(qiáng)的相關(guān)性[12],在分析VGOS的寬帶群時(shí)延精度時(shí)需考慮dT E C的影響。

    通過在干涉測(cè)量相位模型中引入dT E C參數(shù),我們建立了一種基于協(xié)方差分析對(duì)群時(shí)延精度進(jìn)行計(jì)算的方法,并利用國際VGOS實(shí)測(cè)數(shù)據(jù)進(jìn)行了群時(shí)延精度的統(tǒng)計(jì)分析和檢驗(yàn)。該方法可在一系列頻率設(shè)置預(yù)選方案中,比較不同頻率組合的優(yōu)缺點(diǎn),為VGOS觀測(cè)的頻率選擇和設(shè)置提供參考。

    2 干涉測(cè)量相位模型中引入dT E C參數(shù)的群時(shí)延精度計(jì)算方法

    VLBI觀測(cè)中,一般記錄多個(gè)頻率通道的數(shù)據(jù),將多通道數(shù)據(jù)組合在一起進(jìn)行帶寬綜合解算,得到多通道的基線群時(shí)延觀測(cè)量。對(duì)于VGOS的超寬帶頻率觀測(cè),不能忽略電離層的色散效應(yīng)影響,因此在擬合群時(shí)延的同時(shí),引入新的觀測(cè)量dT E C,干涉測(cè)量相位模型可表示為[12]:

    其中,φ為相位(單位為周),f為頻率(單位為GHz),f0為參考頻率(單位為GHz),τg為群時(shí)延(單位為ns),φ0為參考頻率處的相位(單位為周),dT E C為差分電離層總電子含量(單位為TECU)。

    根據(jù)線性最小二乘原理,參數(shù)解算系數(shù)矩陣Aij可表示為:

    其中,k表示通道號(hào),ρk表示第k通道的權(quán)重,參數(shù)β1,β2和β3分別表示τg,φ0和dT E C,它們的偏導(dǎo)數(shù)為:

    其中,σ1,σ2和σ3分別表示τg,φ0和dT E C的中誤差。

    3 理論計(jì)算分析

    在利用實(shí)際觀測(cè)數(shù)據(jù)驗(yàn)證本文方法的理論計(jì)算中,我們選擇代碼為vo0009的實(shí)際VGOS觀測(cè)的4個(gè)頻段(3 GHz,5 GHz,6 GHz和10 GHz)組合成5種頻率設(shè)置模式。每個(gè)頻段有8個(gè)頻率通道,每個(gè)通道帶寬為32 MHz,天空頻率如表1所示。天空頻率表示從觀測(cè)綱要文件(vex文件)中讀取的射頻頻率,VGOS觀測(cè)中采用下邊帶(LSB),以第1通道天空頻率3 032.40 MHz為例,表示該通道的頻率范圍為3 000.40~3 032.40 MHz。本文用A,B,C和D分別代表3 GHz,5 GHz,6 GHz和10 GHz頻段。

    表1 天空頻率設(shè)置 MHz

    5種頻率設(shè)置組合模式如表2所示,模式1包含了A,B,C和D共4個(gè)頻段;模式2包含B,C和D共3個(gè)頻段,缺失A頻段;模式3包含A,C和D共3個(gè)頻段,缺失B頻段;模式4包含A,B和D共3個(gè)頻段,缺失C頻段;模式5包含A,B和C共3個(gè)頻段,缺失D頻段。其中模式1共32個(gè)頻率通道,模式2~5皆為24個(gè)頻率通道。

    近年來,瑞豐生態(tài)在集團(tuán)化發(fā)展過程中積極實(shí)施國際化戰(zhàn)略,不斷加強(qiáng)與國際知名科研機(jī)構(gòu)和優(yōu)秀跨國企業(yè)的技術(shù)與戰(zhàn)略合作,從技術(shù)引進(jìn)、產(chǎn)品開發(fā)、國際貿(mào)易合作、新型服務(wù)模式打造等多方面深度融合發(fā)展,實(shí)現(xiàn)了集團(tuán)企業(yè)的多業(yè)態(tài)組合,形成了更加具有國際競(jìng)爭(zhēng)力的多元發(fā)展格局。此次安徽宣城50萬噸新型肥料基地的投產(chǎn)建設(shè)正是基于“碳能結(jié)晶”工藝技術(shù)的應(yīng)用,是瑞豐生態(tài)與西班牙海拉全面戰(zhàn)略合作的第一個(gè)落地項(xiàng)目。

    表2 5種模式的頻率設(shè)置組合情況

    式(6)中的權(quán)值主要與該通道信噪比有關(guān),單通道信噪比與射電源的相關(guān)流量密度、兩個(gè)測(cè)站的系統(tǒng)等效流量密度(SEFD)、單通道帶寬以及積分時(shí)間有關(guān)[16]。對(duì)同一基線觀測(cè)同一顆射電源,單通道帶寬與積分時(shí)間相同,射電源的相關(guān)流量密度和測(cè)站SEFD對(duì)不同頻率通道的影響復(fù)雜,很難量化,因此本文在假設(shè)單通道信噪比相同的前提下,將公式(6)中的權(quán)ρk設(shè)為相同值(等權(quán))。根據(jù)公式(6)和(7),我們?cè)诘葯?quán)情況下計(jì)算了5種頻率設(shè)置模式中參數(shù)τg,φ0和dT E C的中誤差,見表3。表3第2列為5種頻率設(shè)置模式下的理論有效帶寬,第3列為通道數(shù),第4列為群時(shí)延中誤差,第5列為參考相位中誤差,第6列為dT E C中誤差,第7列為群時(shí)延與dT E C的相關(guān)系數(shù)。有效帶寬計(jì)算公式如下:

    表3 由相位模型估計(jì)的群時(shí)延、參考相位及dT E C的中誤差

    其中N為通道數(shù),fi為第i個(gè)通道的天空頻率,ˉf為天空頻率的平均值。

    從表3中可以看出,不同頻率設(shè)置模式的群時(shí)延中誤差是有差異的,即不同的觀測(cè)頻率設(shè)置會(huì)影響群時(shí)延精度。其中模式1的群時(shí)延誤差最小,群時(shí)延擬合精度最高,接下來的群時(shí)延誤差從小到大的順序是:模式4、模式3、模式2及模式5,也就是說,在VGOS觀測(cè)中,對(duì)群時(shí)延精度影響最大的頻段分別為:D,A,B和C。模式3和4雖然有效帶寬比模式1略大,但是考慮到通道數(shù)量減少帶來的信噪比損失,這兩種模式的群時(shí)延中誤差仍然比模式1稍大。第7列統(tǒng)計(jì)了寬帶群時(shí)延與dT E C的相關(guān)性,當(dāng)缺失A頻段(模式2)時(shí),寬帶群時(shí)延與dT E C的相關(guān)性最大;其次是缺失D頻段(模式5);相關(guān)性最小的是4頻段組合(模式1)。

    值得注意的是,本文假設(shè)每個(gè)通道信噪比皆相同,將式(6)不同通道的權(quán)值設(shè)為相等值,因此,群時(shí)延精度的理論估計(jì)值,僅用于不同頻率設(shè)置模式下的相對(duì)比較,不能對(duì)群時(shí)延精度進(jìn)行絕對(duì)評(píng)估。其目的旨在分析不同頻率設(shè)置模式對(duì)群時(shí)延精度影響的相對(duì)大小,用于觀測(cè)頻率的選擇。但第4章實(shí)測(cè)數(shù)據(jù)處理的結(jié)果基于實(shí)測(cè)信噪比,包含了通道數(shù)、帶寬、積分時(shí)間、測(cè)站天線的系統(tǒng)等效流量密度、射電源相關(guān)流量密度等所有信息,可以用于計(jì)算群時(shí)延的實(shí)際測(cè)量精度。

    4 實(shí)測(cè)數(shù)據(jù)處理

    vo0009是一次24 h的VGOS觀測(cè),觀測(cè)時(shí)間從2020年1月9日UTC 18:00至2020年1月10日18:00,共觀測(cè)了71個(gè)射電源。參與觀測(cè)的6個(gè)臺(tái)站分別為美國的GGAO12M(Gs)站、Westford(Wf)站[14,15]、德國的WETTZ13S(Ws)站[17]、日本的ISHIOKA(Is)站[18,19]、瑞典的ONSA13NE(Oe)和ONSA13SW(Ow)[20]站,見圖1。

    圖1 vo0009觀測(cè)臺(tái)站分布圖

    本文利用中國科學(xué)院上海天文臺(tái)的相關(guān)處理機(jī)平臺(tái)對(duì)原始數(shù)據(jù)進(jìn)行了相關(guān)處理,用HOPS[21]軟件對(duì)本次觀測(cè)所有基線的相關(guān)數(shù)據(jù)進(jìn)行了后處理,得到每條基線所有scan的群時(shí)延觀測(cè)量及誤差。本次觀測(cè)4個(gè)頻段的頻率設(shè)置與表1一致,每個(gè)通道的帶寬為32 MHz。

    我們處理了表2中5種頻率設(shè)置模式的數(shù)據(jù),考慮到某些通道信號(hào)幅值過低、相位校正信號(hào)(PCAL)質(zhì)量不好以及條紋較差等因素的影響,本文僅統(tǒng)計(jì)了條紋擬合過程中質(zhì)量因子高于5的觀測(cè)量。

    群時(shí)延精度與射電源的相關(guān)流量密度有關(guān),而射電源的相關(guān)流量密度又與射電源的總流量密度和基線的長短有關(guān)?;€越長,射電源的相關(guān)流量密度越小,信噪比也越小,群時(shí)延精度越低;反之,基線越短,射電源的相關(guān)流量密度則越大,信噪比也越大,群時(shí)延精度越高。我們處理了vo0009觀測(cè)所有基線的實(shí)測(cè)數(shù)據(jù)。為說明不同基線長度的情況,本文展示了最長基線Gs-Is(9 593 km)以及短基線Gs-Wf(600 km)所有scan的群時(shí)延誤差的統(tǒng)計(jì)結(jié)果。其他基線的結(jié)果由于篇幅所限,本文未做展示。

    4.1 以最長基線Gs-Is為例

    表4給出了Gs-Is基線5種模式對(duì)群時(shí)延誤差影響的統(tǒng)計(jì)情況:第3列為群時(shí)延誤差小于4 ps的觀測(cè)量百分比;第4列為群時(shí)延誤差中值;第5列為以模式1為基準(zhǔn),其他模式群時(shí)延誤差與模式1群時(shí)延誤差比值的時(shí)間序列均值;第6列為以模式1為基準(zhǔn),其他模式群時(shí)延誤差與模式1群時(shí)延誤差比值的時(shí)間序列標(biāo)準(zhǔn)差。

    表4 Gs-Is(9 593 km)基線不同頻率設(shè)置情況對(duì)群時(shí)延誤差影響統(tǒng)計(jì)

    從表中可以看出,當(dāng)基線長度大于9 000 km時(shí),基于模式1的觀測(cè)有64%的群時(shí)延誤差在4 ps以內(nèi),群時(shí)延誤差中值為3.4 ps,滿足VGOS群時(shí)延的設(shè)計(jì)目標(biāo)。基于模式5的觀測(cè),僅2%的scan群時(shí)延誤差在4 ps以內(nèi),群時(shí)延誤差中值為14.4 ps,遠(yuǎn)遠(yuǎn)低于VGOS的預(yù)期;基于模式5的觀測(cè),11%的scan群時(shí)延誤差在4 ps以內(nèi),群時(shí)延誤差中值為10.2 ps,亦不能滿足VGOS的預(yù)期;而基于模式3和4的觀測(cè),群時(shí)延誤差在4 ps以內(nèi)的scan百分比分別為50%和54%,群時(shí)延誤差中值分別為4.0 ps和3.9 ps,可滿足VGOS的預(yù)期。圖2較直觀地顯示了5種頻率設(shè)置模式在不同時(shí)延誤差區(qū)間的觀測(cè)數(shù),橫軸為群時(shí)延誤差,縱軸為群時(shí)延誤差區(qū)間的累計(jì)數(shù)目。

    圖2 Gs-Is(9 593 km)基線不同頻率設(shè)置模式的群時(shí)延誤差累計(jì)數(shù)目分布

    從表4和圖2對(duì)實(shí)際數(shù)據(jù)的處理情況分析可以看出,5種模式的群時(shí)延誤差由小到大的順序?yàn)椋耗J?、模式4、模式3、模式2和模式5;對(duì)應(yīng)4個(gè)頻段對(duì)群時(shí)延精度影響從大到小的順序?yàn)椋篋,A,B和C。與第3章的理論分析一致。

    4.2 以短基線Gs-Wf為例

    同樣地,我們處理了vo0009觀測(cè)中短基線Gs-Wf(600 km)的相關(guān)數(shù)據(jù),得到了5種不同頻率設(shè)置模式下所有scan的群時(shí)延誤差值,統(tǒng)計(jì)情況見表5,表中每列含義與表4相同。Gs-Wf基線在模式1、模式2、模式3、模式4以及模式5中的有效觀測(cè)量分別為684,680,682,681以及686個(gè)。為便于對(duì)比,采用了5種模式的674個(gè)共有觀測(cè)量進(jìn)行統(tǒng)計(jì)分析。

    表5 Gs-Wf(600 km)基線不同頻率設(shè)置情況對(duì)群時(shí)延誤差影響統(tǒng)計(jì)

    表5顯示,本次觀測(cè)中的基線Gs-Wf基于模式1有99%的群時(shí)延誤差在4 ps以內(nèi),群時(shí)延誤差中值為1.5 ps,滿足VGOS群時(shí)延的設(shè)計(jì)目標(biāo)?;谀J?和4,群時(shí)延誤差在4 ps以內(nèi)的scan百分比分別為98%和96%,群時(shí)延誤差中值分別為1.6 ps和1.8 ps,與模式1相比,群時(shí)延誤差在4 ps以內(nèi)的觀測(cè)數(shù)略小,群時(shí)延誤差略大?;谀J?,僅45%的scan群時(shí)延誤差在4 ps以內(nèi),群時(shí)延誤差中值為4.3 ps,接近VGOS的預(yù)期;基于模式5,僅30%的scan群時(shí)延誤差在4 ps以內(nèi),群時(shí)延誤差中值為5.5 ps,未達(dá)到VGOS的預(yù)期。圖3顯示了Gs-Wf基線在5種頻率設(shè)置模式下不同時(shí)延誤差區(qū)間的觀測(cè)數(shù),橫軸為群時(shí)延誤差,縱軸為群時(shí)延誤差區(qū)間的累計(jì)數(shù)目。

    圖3 Gs-Wf基線不同頻率設(shè)置模式的群時(shí)延誤差累計(jì)數(shù)目分布

    從表5和圖3可以看出,5種模式的群時(shí)延誤差由小到大的順序?yàn)椋耗J?、模式4、模式3、模式2和模式5,對(duì)應(yīng)4個(gè)頻段對(duì)群時(shí)延精度影響從大到小的順序?yàn)椋篋,A,B和C。規(guī)律與Gs-Is基線的統(tǒng)計(jì)結(jié)果一致,且與第3章的理論分析相符。

    4.3 其他基線

    我們還處理了vo0009觀測(cè)的其他基線,由于篇幅所限,具體結(jié)果本文未做展示。值得指出的是,其他基線的統(tǒng)計(jì)結(jié)果與Gs-Is,Gs-Wf基線的規(guī)律一致,且與第3章的理論分析相符。

    5 結(jié)論與討論

    傳統(tǒng)測(cè)地VLBI采用S/X雙頻觀測(cè)消除電離層效應(yīng),VGOS觀測(cè)由于頻率范圍跨度很大,必須根據(jù)多頻段干涉條紋同時(shí)提取群時(shí)延觀測(cè)量和差分電離層時(shí)延。本文采用了一種基于VGOS超寬帶頻率序列的協(xié)方差分析方法,通過在干涉測(cè)量相位模型中引入電離層參數(shù),研究了不同頻率組合對(duì)群時(shí)延精度的影響。

    基于代碼為vo0009的國際VGOS觀測(cè),設(shè)計(jì)了5種頻率組合模式,從理論上計(jì)算不同模式對(duì)寬帶群時(shí)延精度的相對(duì)影響。接著處理了本次觀測(cè)所有基線的實(shí)際觀測(cè)數(shù)據(jù),由于篇幅所限,本文僅展示了長基線和短基線的實(shí)測(cè)數(shù)據(jù)分析結(jié)果。我們統(tǒng)計(jì)了5種頻率組合模式下寬帶群時(shí)延誤差的實(shí)測(cè)結(jié)果,與理論分析一致,驗(yàn)證了本文方法的可靠性。

    VGOS觀測(cè)通常采用4個(gè)頻段和固定的頻率設(shè)置,但有時(shí)候也需要針對(duì)實(shí)際觀測(cè)條件進(jìn)行觀測(cè)頻率調(diào)整,例如,有的參加臺(tái)站需要避免特定頻段的無線電干擾,有的臺(tái)站只能記錄3個(gè)頻段,有時(shí)候VGOS臺(tái)站與傳統(tǒng)VLBI臺(tái)站開展混合模式觀測(cè)。上述情況下,都需要重新設(shè)計(jì)觀測(cè)頻率,本文方法可為VGOS觀測(cè)頻率的選擇和設(shè)置提供參考。

    每個(gè)通道信噪比與射電源的相關(guān)流量密度、臺(tái)站SEFD等因素有關(guān),由于這些影響因素的復(fù)雜性,無法具體量化,因此本文的理論計(jì)算忽略了上述因素對(duì)不同頻率通道條紋信噪比的影響。此外,本文未考慮時(shí)延分辨率函數(shù)第2旁瓣峰值對(duì)群時(shí)延模糊度的影響。后續(xù)將考慮上述因素,進(jìn)一步優(yōu)化VGOS頻率的設(shè)置方法。

    致謝

    感謝美國Haystack天文臺(tái)提供VGOS實(shí)測(cè)數(shù)據(jù),為本文開展VGOS數(shù)據(jù)的相關(guān)后處理和群時(shí)延精度計(jì)算提供了便利。

    猜你喜歡
    群時(shí)延基線頻段
    群時(shí)延傅里葉分解模型及其估計(jì)方法*
    適用于MAUV的變基線定位系統(tǒng)
    gPhone重力儀的面波頻段響應(yīng)實(shí)測(cè)研究
    地震研究(2021年1期)2021-04-13 01:04:56
    航天技術(shù)與甚長基線陣的結(jié)合探索
    科學(xué)(2020年5期)2020-11-26 08:19:14
    群時(shí)延對(duì)混合擴(kuò)頻信號(hào)接收的影響分析
    電子科技(2018年4期)2018-04-08 02:06:40
    一種改進(jìn)的干涉儀測(cè)向基線設(shè)計(jì)方法
    推擠的5GHz頻段
    CHIP新電腦(2016年3期)2016-03-10 14:07:52
    TD—LTE在D頻段和F頻段的覆蓋能力差異
    中國新通信(2015年1期)2015-05-30 10:30:46
    強(qiáng)干擾下群時(shí)延失真對(duì)偽碼測(cè)距性能的影響研究
    技術(shù)狀態(tài)管理——對(duì)基線更改的控制
    航天器工程(2014年5期)2014-03-11 16:35:50
    亚洲国产高清在线一区二区三 | 最新在线观看一区二区三区| 热re99久久国产66热| 在线观看免费视频日本深夜| 两个人免费观看高清视频| 久久中文字幕一级| 国产日本99.免费观看| 可以在线观看毛片的网站| 欧美成人午夜精品| 巨乳人妻的诱惑在线观看| 久久99热这里只有精品18| 搡老妇女老女人老熟妇| 国产v大片淫在线免费观看| 国产亚洲精品综合一区在线观看 | 国产成年人精品一区二区| 国产高清有码在线观看视频 | 国产成人啪精品午夜网站| 免费在线观看视频国产中文字幕亚洲| 日韩成人在线观看一区二区三区| 少妇裸体淫交视频免费看高清 | 无人区码免费观看不卡| 在线国产一区二区在线| 男人舔奶头视频| 长腿黑丝高跟| 国产国语露脸激情在线看| 怎么达到女性高潮| 天天躁夜夜躁狠狠躁躁| 哪里可以看免费的av片| 国产精品美女特级片免费视频播放器 | 99国产极品粉嫩在线观看| 亚洲国产欧洲综合997久久, | 久久久久亚洲av毛片大全| 久久久水蜜桃国产精品网| 亚洲成av人片免费观看| 亚洲 欧美一区二区三区| 久久婷婷成人综合色麻豆| 国产精品免费一区二区三区在线| tocl精华| 一级毛片精品| 男人舔女人的私密视频| 国产黄a三级三级三级人| 精品欧美国产一区二区三| 男女那种视频在线观看| 成人av一区二区三区在线看| 亚洲无线在线观看| 亚洲免费av在线视频| 国产亚洲精品综合一区在线观看 | av在线天堂中文字幕| 久久这里只有精品19| 国产激情偷乱视频一区二区| 亚洲精品一区av在线观看| 日韩视频一区二区在线观看| 老司机午夜福利在线观看视频| 国产成+人综合+亚洲专区| 2021天堂中文幕一二区在线观 | 一级黄色大片毛片| 美女高潮喷水抽搐中文字幕| 日韩有码中文字幕| 国产激情偷乱视频一区二区| 精品久久久久久,| 午夜亚洲福利在线播放| 国产激情久久老熟女| 婷婷亚洲欧美| 久久99热这里只有精品18| 欧美黑人巨大hd| 午夜成年电影在线免费观看| 日韩中文字幕欧美一区二区| 国产午夜精品久久久久久| 别揉我奶头~嗯~啊~动态视频| 高潮久久久久久久久久久不卡| 欧美三级亚洲精品| 99在线人妻在线中文字幕| 99在线视频只有这里精品首页| 黄色a级毛片大全视频| 精品一区二区三区视频在线观看免费| 美女高潮到喷水免费观看| 日韩有码中文字幕| 又黄又粗又硬又大视频| 国产av一区二区精品久久| 18美女黄网站色大片免费观看| 天堂影院成人在线观看| 999久久久精品免费观看国产| 久久久久久久久久黄片| 99精品久久久久人妻精品| 日日摸夜夜添夜夜添小说| 999久久久国产精品视频| 色综合站精品国产| 中出人妻视频一区二区| 精品福利观看| 日本精品一区二区三区蜜桃| 精品熟女少妇八av免费久了| 欧美绝顶高潮抽搐喷水| 欧美午夜高清在线| 久99久视频精品免费| 久久国产亚洲av麻豆专区| 波多野结衣高清无吗| 国产精品av久久久久免费| 成人国语在线视频| 国产精品99久久99久久久不卡| 丰满的人妻完整版| 欧美日韩福利视频一区二区| 看片在线看免费视频| 国产精品影院久久| 69av精品久久久久久| 国产精品 欧美亚洲| 一二三四在线观看免费中文在| 国产私拍福利视频在线观看| 制服诱惑二区| 亚洲男人的天堂狠狠| 午夜日韩欧美国产| 女性被躁到高潮视频| 国内少妇人妻偷人精品xxx网站 | 夜夜躁狠狠躁天天躁| 久久精品国产综合久久久| 久热爱精品视频在线9| 神马国产精品三级电影在线观看 | 国产高清激情床上av| 十八禁网站免费在线| 听说在线观看完整版免费高清| 黄频高清免费视频| 在线观看日韩欧美| 亚洲第一欧美日韩一区二区三区| 午夜免费鲁丝| 亚洲国产精品合色在线| 国产爱豆传媒在线观看 | 亚洲中文av在线| 精品欧美国产一区二区三| 精品高清国产在线一区| 日韩精品免费视频一区二区三区| 男人的好看免费观看在线视频 | 精品熟女少妇八av免费久了| 欧美精品啪啪一区二区三区| 在线免费观看的www视频| 高清毛片免费观看视频网站| 黑人操中国人逼视频| 波多野结衣av一区二区av| 99国产精品一区二区蜜桃av| 看免费av毛片| 18禁黄网站禁片免费观看直播| 男女视频在线观看网站免费 | 黄色成人免费大全| 亚洲中文av在线| 久久久国产成人精品二区| 亚洲午夜精品一区,二区,三区| 免费电影在线观看免费观看| 精品一区二区三区四区五区乱码| 国产精品99久久99久久久不卡| 国产成年人精品一区二区| 人人妻人人澡人人看| 无限看片的www在线观看| 欧美黑人精品巨大| 99久久国产精品久久久| 99国产精品一区二区三区| 国产蜜桃级精品一区二区三区| 成年人黄色毛片网站| 中文字幕高清在线视频| 99re在线观看精品视频| 在线观看日韩欧美| 国产乱人伦免费视频| 亚洲久久久国产精品| x7x7x7水蜜桃| 熟女电影av网| 一夜夜www| 天堂影院成人在线观看| 亚洲成人精品中文字幕电影| 午夜精品在线福利| 女性被躁到高潮视频| 亚洲成人国产一区在线观看| 黄色片一级片一级黄色片| 不卡av一区二区三区| 午夜免费观看网址| 久久久国产成人精品二区| 精品不卡国产一区二区三区| 亚洲九九香蕉| 亚洲欧美日韩高清在线视频| 美女午夜性视频免费| 黄片大片在线免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 久久99热这里只有精品18| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品99久久99久久久不卡| 老汉色av国产亚洲站长工具| 久久草成人影院| 亚洲中文字幕日韩| 极品教师在线免费播放| 精品国产乱子伦一区二区三区| 国产精品影院久久| a级毛片在线看网站| 99精品在免费线老司机午夜| 久久精品国产亚洲av香蕉五月| 国产真实乱freesex| 又黄又粗又硬又大视频| 女生性感内裤真人,穿戴方法视频| 亚洲成国产人片在线观看| 一级黄色大片毛片| 精品免费久久久久久久清纯| 亚洲av中文字字幕乱码综合 | 精品国产乱子伦一区二区三区| 久久国产乱子伦精品免费另类| 免费人成视频x8x8入口观看| 免费一级毛片在线播放高清视频| 午夜免费观看网址| 国产v大片淫在线免费观看| 精品一区二区三区视频在线观看免费| 国产精品av久久久久免费| 最新美女视频免费是黄的| 午夜福利在线在线| 麻豆国产av国片精品| 无限看片的www在线观看| 欧美黄色片欧美黄色片| 国产亚洲欧美精品永久| 神马国产精品三级电影在线观看 | 中文字幕最新亚洲高清| 一级a爱片免费观看的视频| 国产精华一区二区三区| 美女国产高潮福利片在线看| 亚洲狠狠婷婷综合久久图片| 怎么达到女性高潮| 久久精品国产清高在天天线| 国产亚洲欧美精品永久| 国产三级在线视频| 精品熟女少妇八av免费久了| 夜夜看夜夜爽夜夜摸| 亚洲性夜色夜夜综合| 日韩欧美国产在线观看| 欧美在线黄色| 日韩一卡2卡3卡4卡2021年| 又紧又爽又黄一区二区| 两个人视频免费观看高清| 精品国产乱子伦一区二区三区| 18禁黄网站禁片免费观看直播| 久久香蕉激情| 亚洲欧美精品综合一区二区三区| 国产一区二区三区在线臀色熟女| 免费观看精品视频网站| 他把我摸到了高潮在线观看| 国产日本99.免费观看| 一区二区三区激情视频| 亚洲五月天丁香| 亚洲国产日韩欧美精品在线观看 | 国产亚洲精品av在线| 免费在线观看完整版高清| 国产区一区二久久| 悠悠久久av| 啪啪无遮挡十八禁网站| netflix在线观看网站| 久久久久久久久久黄片| 久久久久精品国产欧美久久久| 久久久久久九九精品二区国产 | 免费电影在线观看免费观看| 亚洲黑人精品在线| 国产国语露脸激情在线看| 国产精品野战在线观看| 欧美三级亚洲精品| 午夜免费观看网址| 一边摸一边抽搐一进一小说| 老司机福利观看| 精品一区二区三区av网在线观看| 熟女电影av网| 欧美激情 高清一区二区三区| 两个人视频免费观看高清| 亚洲人成伊人成综合网2020| 亚洲av成人一区二区三| 免费看a级黄色片| 两性夫妻黄色片| 岛国在线观看网站| 免费看美女性在线毛片视频| 99精品在免费线老司机午夜| 天堂动漫精品| 成人18禁高潮啪啪吃奶动态图| 一级a爱片免费观看的视频| 18禁黄网站禁片午夜丰满| 日本成人三级电影网站| 夜夜爽天天搞| 欧美乱色亚洲激情| 亚洲精品国产精品久久久不卡| 欧美另类亚洲清纯唯美| 男人舔奶头视频| 欧美乱妇无乱码| 老司机在亚洲福利影院| 亚洲av五月六月丁香网| 嫩草影院精品99| 我的亚洲天堂| 美女高潮喷水抽搐中文字幕| 观看免费一级毛片| 色综合亚洲欧美另类图片| 可以免费在线观看a视频的电影网站| 久久久久久人人人人人| 国产欧美日韩一区二区精品| 国产精品亚洲一级av第二区| 满18在线观看网站| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久久大精品| 中文字幕av电影在线播放| 久久国产精品男人的天堂亚洲| 美女 人体艺术 gogo| 精品一区二区三区视频在线观看免费| 桃红色精品国产亚洲av| 熟女电影av网| 男人的好看免费观看在线视频 | 在线永久观看黄色视频| 狂野欧美激情性xxxx| 久久欧美精品欧美久久欧美| 人人妻人人澡人人看| 欧美日韩黄片免| 亚洲五月天丁香| 在线十欧美十亚洲十日本专区| or卡值多少钱| 91在线观看av| 国产精品爽爽va在线观看网站 | 国产欧美日韩一区二区三| 成人三级做爰电影| 国产又色又爽无遮挡免费看| 中出人妻视频一区二区| 国产在线观看jvid| 91麻豆av在线| 国产一区在线观看成人免费| 欧美丝袜亚洲另类 | 亚洲三区欧美一区| 免费看十八禁软件| 日本成人三级电影网站| 99riav亚洲国产免费| 欧美精品啪啪一区二区三区| av视频在线观看入口| 国产精品亚洲美女久久久| 欧美日韩中文字幕国产精品一区二区三区| 日日夜夜操网爽| 深夜精品福利| 一级毛片女人18水好多| 国产午夜精品久久久久久| 国产成人av激情在线播放| 免费电影在线观看免费观看| 男女那种视频在线观看| 欧美黑人欧美精品刺激| 日韩大尺度精品在线看网址| 亚洲熟妇熟女久久| 看片在线看免费视频| 日韩三级视频一区二区三区| 国产1区2区3区精品| 波多野结衣高清作品| 最近最新中文字幕大全电影3 | 99久久综合精品五月天人人| 久久国产精品男人的天堂亚洲| 欧美中文日本在线观看视频| 久久国产乱子伦精品免费另类| 97超级碰碰碰精品色视频在线观看| 自线自在国产av| 精品熟女少妇八av免费久了| 伊人久久大香线蕉亚洲五| 国产欧美日韩一区二区精品| 国产爱豆传媒在线观看 | 99国产精品99久久久久| or卡值多少钱| 国产伦在线观看视频一区| 久久久久久国产a免费观看| 国产激情欧美一区二区| 久久久久久大精品| 好看av亚洲va欧美ⅴa在| 听说在线观看完整版免费高清| 色综合婷婷激情| 久久久久久大精品| 香蕉久久夜色| 国产亚洲精品久久久久5区| 韩国av一区二区三区四区| 国产一区二区三区在线臀色熟女| 国产伦在线观看视频一区| 叶爱在线成人免费视频播放| 亚洲va日本ⅴa欧美va伊人久久| 日韩欧美免费精品| 色老头精品视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 两个人免费观看高清视频| 法律面前人人平等表现在哪些方面| 久久国产精品男人的天堂亚洲| 国产v大片淫在线免费观看| 97人妻精品一区二区三区麻豆 | 成人国产一区最新在线观看| x7x7x7水蜜桃| 欧美av亚洲av综合av国产av| 成熟少妇高潮喷水视频| 色综合亚洲欧美另类图片| 伦理电影免费视频| 国产成年人精品一区二区| 欧美一级a爱片免费观看看 | 99久久精品国产亚洲精品| 国产精品影院久久| 免费看日本二区| 美女免费视频网站| 成人18禁在线播放| 99精品欧美一区二区三区四区| 亚洲人成网站高清观看| 国产黄a三级三级三级人| 看片在线看免费视频| 90打野战视频偷拍视频| 搞女人的毛片| 一级片免费观看大全| 国产精品98久久久久久宅男小说| a在线观看视频网站| 国产精品免费视频内射| 两性夫妻黄色片| 老司机靠b影院| 国产精品美女特级片免费视频播放器 | 成人一区二区视频在线观看| 国产一区二区三区在线臀色熟女| 不卡av一区二区三区| 日日干狠狠操夜夜爽| 美女国产高潮福利片在线看| 欧美成人一区二区免费高清观看 | 99热6这里只有精品| 中文在线观看免费www的网站 | 欧美精品亚洲一区二区| 久久精品国产清高在天天线| 久久这里只有精品19| 亚洲五月色婷婷综合| 国产精品久久视频播放| 成人精品一区二区免费| 欧美黄色片欧美黄色片| 国产激情久久老熟女| 美女大奶头视频| 精品无人区乱码1区二区| 免费在线观看影片大全网站| 欧美日韩福利视频一区二区| 亚洲 国产 在线| aaaaa片日本免费| 欧美av亚洲av综合av国产av| 校园春色视频在线观看| 久久99热这里只有精品18| xxx96com| 19禁男女啪啪无遮挡网站| 99久久99久久久精品蜜桃| 国产精品av久久久久免费| 亚洲最大成人中文| 青草久久国产| 国产成人欧美| bbb黄色大片| 国产三级在线视频| 国产亚洲欧美精品永久| 国产精品久久久久久精品电影 | 听说在线观看完整版免费高清| 日韩国内少妇激情av| 中亚洲国语对白在线视频| 亚洲国产高清在线一区二区三 | 一级黄色大片毛片| 成年人黄色毛片网站| 亚洲 国产 在线| 热99re8久久精品国产| 国产亚洲精品第一综合不卡| 中出人妻视频一区二区| 色综合亚洲欧美另类图片| 精品国产乱子伦一区二区三区| 91成年电影在线观看| 欧美乱色亚洲激情| av超薄肉色丝袜交足视频| 精品国产亚洲在线| 国产国语露脸激情在线看| 亚洲成人精品中文字幕电影| 法律面前人人平等表现在哪些方面| 精品乱码久久久久久99久播| 欧美最黄视频在线播放免费| 嫩草影院精品99| 色老头精品视频在线观看| 十八禁人妻一区二区| 一进一出抽搐gif免费好疼| 夜夜夜夜夜久久久久| 国产真人三级小视频在线观看| 亚洲欧洲精品一区二区精品久久久| 波多野结衣高清作品| 国产亚洲精品综合一区在线观看 | 亚洲avbb在线观看| 国产一卡二卡三卡精品| 精品国内亚洲2022精品成人| 精品久久久久久久毛片微露脸| 757午夜福利合集在线观看| a在线观看视频网站| 999久久久国产精品视频| 精品一区二区三区四区五区乱码| 国产三级在线视频| 亚洲男人的天堂狠狠| 日本精品一区二区三区蜜桃| 午夜日韩欧美国产| 99国产精品一区二区蜜桃av| 精品国产超薄肉色丝袜足j| 热99re8久久精品国产| www日本黄色视频网| 级片在线观看| 久久久久久久午夜电影| 亚洲自偷自拍图片 自拍| 18美女黄网站色大片免费观看| 一进一出抽搐gif免费好疼| 女人被狂操c到高潮| 久久天堂一区二区三区四区| 久久热在线av| 免费高清在线观看日韩| 白带黄色成豆腐渣| 久久九九热精品免费| 国产男靠女视频免费网站| 久久伊人香网站| 久久精品影院6| 搞女人的毛片| 一本久久中文字幕| 欧美+亚洲+日韩+国产| 一卡2卡三卡四卡精品乱码亚洲| 欧美黄色片欧美黄色片| 91在线观看av| 久久精品91无色码中文字幕| 日韩免费av在线播放| 色av中文字幕| 日本精品一区二区三区蜜桃| 国产99久久九九免费精品| 视频区欧美日本亚洲| 激情在线观看视频在线高清| 国产色视频综合| 亚洲精品一区av在线观看| www.熟女人妻精品国产| 19禁男女啪啪无遮挡网站| 欧美精品啪啪一区二区三区| 久久午夜亚洲精品久久| 日韩成人在线观看一区二区三区| 成人18禁高潮啪啪吃奶动态图| 亚洲狠狠婷婷综合久久图片| 亚洲专区国产一区二区| 美女 人体艺术 gogo| 两个人免费观看高清视频| 国产激情偷乱视频一区二区| 亚洲aⅴ乱码一区二区在线播放 | 亚洲男人天堂网一区| 韩国av一区二区三区四区| 麻豆av在线久日| 亚洲精品在线美女| 国产精品亚洲av一区麻豆| 久久天堂一区二区三区四区| 亚洲午夜理论影院| a级毛片在线看网站| 男人舔女人的私密视频| 国产成人啪精品午夜网站| 国产亚洲欧美在线一区二区| 成人国产综合亚洲| 亚洲人成电影免费在线| 亚洲精品久久成人aⅴ小说| 黄色毛片三级朝国网站| 午夜亚洲福利在线播放| 亚洲精品国产区一区二| 久久久久国产一级毛片高清牌| 男女床上黄色一级片免费看| 欧美 亚洲 国产 日韩一| 欧美成狂野欧美在线观看| 国产成人系列免费观看| 淫妇啪啪啪对白视频| 视频区欧美日本亚洲| 久热爱精品视频在线9| 免费在线观看亚洲国产| 91麻豆精品激情在线观看国产| 欧美乱码精品一区二区三区| 高清毛片免费观看视频网站| 波多野结衣av一区二区av| 午夜免费鲁丝| 欧美在线一区亚洲| 天天躁夜夜躁狠狠躁躁| 丰满人妻熟妇乱又伦精品不卡| 亚洲狠狠婷婷综合久久图片| 黄网站色视频无遮挡免费观看| 两人在一起打扑克的视频| 国产91精品成人一区二区三区| 一级毛片高清免费大全| 麻豆av在线久日| 国产亚洲精品第一综合不卡| 我的亚洲天堂| 69av精品久久久久久| 自线自在国产av| 午夜久久久久精精品| 少妇裸体淫交视频免费看高清 | 美女大奶头视频| 免费在线观看日本一区| 国产成人精品无人区| 好男人电影高清在线观看| 成人特级黄色片久久久久久久| 性欧美人与动物交配| 亚洲久久久国产精品| 婷婷亚洲欧美| 精品无人区乱码1区二区| 99热只有精品国产| 亚洲国产精品合色在线| 99久久国产精品久久久| 免费电影在线观看免费观看| 国产精品乱码一区二三区的特点| 18禁黄网站禁片午夜丰满| 91麻豆精品激情在线观看国产| 国产免费男女视频| 午夜精品久久久久久毛片777| 一区二区日韩欧美中文字幕| 国产熟女xx| 成年人黄色毛片网站| avwww免费| 精品国产国语对白av| 欧美成狂野欧美在线观看| av欧美777| 欧美乱码精品一区二区三区| 真人一进一出gif抽搐免费| 日韩欧美一区视频在线观看| 国产亚洲精品一区二区www| 国产av又大| 久久99热这里只有精品18| 成人永久免费在线观看视频| 中文字幕人妻熟女乱码| 久久天堂一区二区三区四区| 亚洲va日本ⅴa欧美va伊人久久| 国产精品美女特级片免费视频播放器 | 欧美性长视频在线观看| 日本一区二区免费在线视频| 国产精品,欧美在线| 99国产精品一区二区三区| 一进一出抽搐动态| 亚洲第一青青草原| 亚洲国产精品合色在线| 日韩一卡2卡3卡4卡2021年|