《電子技術(shù)應(yīng)用》
您所在的位置:首頁 > 模擬設(shè)計(jì) > 設(shè)計(jì)應(yīng)用 > 高軌SAR信號(hào)對(duì)GNSS接收機(jī)性能影響評(píng)估方法
高軌SAR信號(hào)對(duì)GNSS接收機(jī)性能影響評(píng)估方法
2020年電子技術(shù)應(yīng)用第3期
趙冠先,趙思浩,崔曉偉
清華大學(xué) 電子工程系,北京100084
摘要: 由于高軌SAR衛(wèi)星信號(hào)設(shè)計(jì)頻段可能與多種已投入使用的GNSS信號(hào)所在頻段共用或重疊,導(dǎo)致地面GNSS接收機(jī)受到由高軌SAR信號(hào)產(chǎn)生的脈沖射頻干擾。為評(píng)估高軌SAR信號(hào)對(duì)GNSS地面接收機(jī)性能的影響,提出了一種基于ITU相關(guān)標(biāo)準(zhǔn)的參數(shù)化理論模型。該模型使用“源-路徑-接收機(jī)”的分析方法,分別從高軌SAR信號(hào)、信號(hào)傳播路徑以及用于干擾評(píng)估的接收機(jī)模型三個(gè)方面提取參數(shù),然后根據(jù)ITU相關(guān)標(biāo)準(zhǔn),評(píng)估高軌SAR信號(hào)作為加性脈沖射頻干擾時(shí),對(duì)GNSS接收機(jī)性能產(chǎn)生的影響。以具有代表性的B3I信號(hào)為例,評(píng)估在設(shè)定參數(shù)條件下,高軌SAR信號(hào)對(duì)不同類型的B3I接收機(jī)性能的影響。結(jié)果表明,在當(dāng)前設(shè)定的高軌SAR信號(hào)干擾存在情況下,接收機(jī)性能略有下降,但仍然能夠正常工作。
中圖分類號(hào): TN965.5
文獻(xiàn)標(biāo)識(shí)碼: A
DOI:10.16157/j.issn.0258-7998.191049
中文引用格式: 趙冠先,趙思浩,崔曉偉. 高軌SAR信號(hào)對(duì)GNSS接收機(jī)性能影響評(píng)估方法[J].電子技術(shù)應(yīng)用,2020,46(3):14-18.
英文引用格式: Zhao Guanxian,Zhao Sihao,Cui Xiaowei. Evaluation method for the performance of GNSS receiver influenced by GEO SAR signal[J]. Application of Electronic Technique,2020,46(3):14-18.
Evaluation method for the performance of GNSS receiver influenced by GEO SAR signal
Zhao Guanxian,Zhao Sihao,Cui Xiaowei
Department of Electronic Engineering, Tsinghua University,Beijing 100084,China
Abstract: Because the designed frequency band of Geosynchronous synthetic aperture radar(GEO SAR) satellite signal possibly shares or overlaps with the frequency band of various Global Navigation Satellite System(GNSS) signals which have been put into use, the ground GNSS receiver may be subject to the pulse radio frequency interference generated by GEO SAR signal. In order to evaluate the influence of GEO SAR signal on the performance of GNSS ground receiver, a parametric theoretical model based on International Telecommunication Union(ITU) related standards is proposed in this paper. The model uses the "source-path-receiver" analysis method to extract parameters from the three aspects of GEO SAR signal, signal propagation path and receiver model for interference evaluation. Then, according to ITU standards, the influence of high-orbit SAR signal as additive pulse radio frequency interference on the performance of GNSS receiver is evaluated. Taking the representative B3I signal as an example, this paper evaluates the performance of different types of B3I receivers for GEO SAR signal under the condition of the assumed design parameters. The results show that the performance of the receiver decreases slightly in the presence of high-orbit SAR signal interference, but it can still work normally.
Key words : satellite navigation;GNSS receiver;radio frequency interference

0 引言

    地球同步軌道合成孔徑雷達(dá)(Geosynchronous Synthetic Aperture Radar,GEO SAR)即高軌SAR,是一種新體制星載SAR,能夠?qū)⒌厍蛲杰壍赖能壩粌?yōu)勢(shì)與SAR衛(wèi)星穿透能力強(qiáng)、不受氣象與光照條件影響的優(yōu)勢(shì)結(jié)合起來,實(shí)現(xiàn)對(duì)地全天候、大范圍、高時(shí)間分辨率的觀測(cè)[1]。考慮到信號(hào)衰減及穿透能力,高軌SAR衛(wèi)星信號(hào)計(jì)劃使用L頻段,該頻段上現(xiàn)已存在多種全球?qū)Ш叫l(wèi)星系統(tǒng)(Global Navigation Satellite System,GNSS)業(yè)務(wù),高軌SAR計(jì)劃頻段與多種GNSS信號(hào)頻段重疊或鄰近。由于衛(wèi)星導(dǎo)航信號(hào)自身的特點(diǎn),其到達(dá)地面的功率非常小,非常容易受到各種無線電干擾[2]。為保證地面GNSS接收機(jī)在高軌SAR信號(hào)存在時(shí)仍能夠正常工作,有必要對(duì)高軌SAR信號(hào)對(duì)GNSS地面接收機(jī)的性能影響進(jìn)行分析與評(píng)估,結(jié)果可作為高軌SAR項(xiàng)目立項(xiàng)及未來開展國內(nèi)外頻率協(xié)調(diào)工作的理論依據(jù)。

    對(duì)于GNSS接收機(jī)來說,高軌SAR信號(hào)可看作是脈沖射頻干擾(Radio Frequency Interference,RFI)。近年來,已有許多有關(guān)脈沖射頻干擾對(duì)導(dǎo)航接收機(jī)影響的分析方法,國際電聯(lián)(International Telecommunication Union,ITU)給出了多個(gè)相關(guān)的標(biāo)準(zhǔn)可供借鑒:其中ITU-R M.1902給出了工作在1 215~1 300 MHz頻段的地面GNSS接收機(jī)的工作特性和保護(hù)標(biāo)準(zhǔn),其中特定的接收機(jī)工作特性在分析高軌SAR信號(hào)干擾時(shí)具有參考意義[3]。ITU-R M.2220包含了用于干擾評(píng)估、計(jì)算綜合脈沖干擾參數(shù)的背景知識(shí)和基本計(jì)算方法,舉例計(jì)算了在多種射頻干擾存在的復(fù)雜環(huán)境下,不同類型的GNSS接收機(jī)受干擾情況的評(píng)估方法和結(jié)果[4]。ITU-R M.2030針對(duì)兩種不同類型的基本GNSS接收機(jī)給出了等效載噪比衰減公式,提出評(píng)估加性脈沖RFI影響的方法[5]。ITU-R RS.2311包含了主動(dòng)地球探測(cè)衛(wèi)星業(yè)務(wù)(Earth Exploration Satellite Service,EESS)對(duì)衛(wèi)星導(dǎo)航接收機(jī)影響的兩個(gè)獨(dú)立研究結(jié)果[6]。本文整理了上述文獻(xiàn),提出用于評(píng)估高軌SAR信號(hào)對(duì)GNSS接收機(jī)性能的影響的方法。

1 干擾場(chǎng)景及評(píng)估方法

    本文采用典型的“源-路徑-接收機(jī)”方法評(píng)估高軌SAR信號(hào)對(duì)GNSS接收機(jī)性能影響,這種方法需要生成(仿真)以下三個(gè)基本的且具有相關(guān)性的關(guān)鍵因素,或從已有的標(biāo)準(zhǔn)中收集相關(guān)信息:

    (1)高軌SAR信號(hào)的位置及相關(guān)參數(shù);

    (2)高軌SAR信號(hào)及GNSS信號(hào)在干擾場(chǎng)景下的傳播路徑中的參數(shù);

    (3)用于接收機(jī)性能評(píng)估的GNSS接收機(jī)模型。

    干擾場(chǎng)景構(gòu)成了這三個(gè)關(guān)鍵因素的基礎(chǔ),本文分析的干擾場(chǎng)景為高軌SAR衛(wèi)星載荷工作時(shí)對(duì)地面接收機(jī)進(jìn)行持續(xù)干擾。為了便于分析,將其他脈沖干擾源以及連續(xù)干擾源對(duì)GNSS接收機(jī)性能產(chǎn)生的影響作為干擾基準(zhǔn),單顆高軌SAR衛(wèi)星信號(hào)作為加性干擾源。高軌SAR信號(hào)的參數(shù)包括衛(wèi)星軌道、天線特性、發(fā)射波形、功率電平以及頻率;傳播路徑考慮了各種信號(hào)傳播過程當(dāng)中的衰減,提供了計(jì)算到達(dá)接收機(jī)的信號(hào)功率的方法;而接收機(jī)模型提供了用于確定高軌SAR信號(hào)對(duì)于接收機(jī)影響的性能參數(shù)。

    評(píng)估流程如下:首先將高軌SAR信號(hào)、傳播路徑、接收機(jī)模型參數(shù)匯總,計(jì)算到達(dá)接收機(jī)天線的接收功率;接著將接收功率與接收機(jī)生存電平與壓縮電平進(jìn)行比較,若均不超過門限值,則進(jìn)行等效載噪比計(jì)算;最后通過等效載噪比評(píng)估高軌SAR信號(hào)對(duì)GNSS地面接收機(jī)的影響,評(píng)估流程圖如圖1所示。

wxdh4-t1.gif

2 參數(shù)匯總

2.1 高軌SAR信號(hào)特性及參數(shù)

    SAR一般通過發(fā)射矩形包絡(luò)的線性調(diào)頻信號(hào)獲得距離向的高分辨率,矩形包絡(luò)且幅度歸一化線性調(diào)頻信號(hào)的復(fù)指數(shù)形式為[7]

wxdh4-gs1.gif

wxdh4-gs1-x1.gif

wxdh4-b1.gif

    高軌SAR設(shè)計(jì)頻段與當(dāng)前已經(jīng)存在的多種GNSS信號(hào)存在頻率重疊或相近的情況,如圖2所示。

wxdh4-t2.gif

    由圖2可以看出,B3I信號(hào)頻段完全被高軌SAR信號(hào)頻段覆蓋,具有代表性,故本文選取北斗B3I信號(hào)通用接收機(jī)進(jìn)行分析。

2.2 傳播路徑分析

    在理想情況下,高軌SAR及GNSS信號(hào)在自由空間內(nèi)進(jìn)行傳播。電磁波在自由空間中傳播時(shí)僅存在因信號(hào)能量擴(kuò)散引發(fā)的衰減,不存在任何其他形式的損耗。傳播模型是在Friis傳輸公式的基礎(chǔ)上經(jīng)過一定簡(jiǎn)化、處理得來的,該模型的對(duì)數(shù)形式如下[8]

    wxdh4-gs2.gif

其中,PRmax為接收最大功率(dBW);PTmax為發(fā)射最大功率(dBW);GTmax為發(fā)射天線最大增益(dB);GRmax為接收天線最大增益(dB);Lpolar為極化失配損耗(dB),當(dāng)SAR為線極化而B3I接收機(jī)采用圓極化時(shí),該值取1.46 dB;LA為大氣損耗,約為2 dB;Lpath為路徑損耗(dB),其計(jì)算公式如下:

    wxdh4-gs3.gif

其中,fc為信號(hào)頻率(MHz);D為傳輸距離(km),即衛(wèi)星到接收機(jī)的距離。

    根據(jù)自由空間傳播模型,在已知高軌SAR信號(hào)及GNSS接收機(jī)相關(guān)參數(shù)的情況下,能夠計(jì)算出信號(hào)到達(dá)接收機(jī)的功率電平,以此與GNSS接收機(jī)的生存電平、輸入壓縮電平作比較判斷高軌SAR信號(hào)是否對(duì)GNSS接收機(jī)產(chǎn)生影響。

2.3 用于脈沖射頻干擾評(píng)估的接收機(jī)模型

    本文主要分析高軌SAR信號(hào)對(duì)B3I通用接收機(jī)的影響,通用接收機(jī)設(shè)計(jì)用于車輛導(dǎo)航、步行導(dǎo)航、一般定位等。根據(jù)文獻(xiàn)[3],給出B3I信號(hào)通用接收機(jī)參數(shù)如表2所示。

wxdh4-b2.gif

3 高軌SAR信號(hào)對(duì)GNSS接收機(jī)影響評(píng)估模型

3.1 接收機(jī)功率輸入分析

    根據(jù)對(duì)高軌SAR信號(hào)參數(shù)、傳播路徑以及接收機(jī)參數(shù),首先判斷接收到的SAR信號(hào)峰值功率是否超過接收機(jī)的生存電平。由第2節(jié)提到的自由空間傳播模型,將表1中的高軌SAR信號(hào)參數(shù)及表2中的接收機(jī)參數(shù)帶入公式(2),得到接收機(jī)接收功率:PRmax=-97.36 dBW。未超過接收機(jī)生存電平-20 dBW及輸入壓縮電平-70 dBW,進(jìn)行下一步分析。

3.2 接收機(jī)等效噪聲功率變化分析[4,5]

    前述各標(biāo)準(zhǔn)中采用的分析指標(biāo)均為等效噪聲功率的變化值。由于通用接收機(jī)實(shí)現(xiàn)方式多樣,本文分別對(duì)使用脈沖消隱方式和使用脈沖飽和方式的接收機(jī)進(jìn)行等效噪聲功率變化分析。

3.2.1 脈沖消隱接收機(jī)

    有些接收機(jī)會(huì)使用快速數(shù)字脈沖消隱器減少脈沖射頻干擾的影響。脈沖消隱是指當(dāng)接收到的信號(hào)功率高于消隱閾值時(shí),將信號(hào)、干擾及噪聲全部歸零;而低于消隱閾值的信號(hào)、干擾及噪聲則正常通過。這種接收機(jī)稱為脈沖消隱接收機(jī)。

    對(duì)于這種類型的接收機(jī),其相關(guān)器輸出處的等效噪聲功率譜密度為:

    wxdh4-gs4.gif

其中,N0是接收機(jī)系統(tǒng)熱噪聲功率譜密度(W/Hz);I0,WB是所有連續(xù)RFI干擾的功率譜密度,本文主要分析高軌SAR信號(hào)的強(qiáng)脈沖干擾,因此將這一項(xiàng)歸納到場(chǎng)景基準(zhǔn)RFI當(dāng)中,不做詳細(xì)討論;RI是低于消隱器閾值的平均脈沖干擾功率譜密度與接收機(jī)噪聲溫度之比,與上一項(xiàng)相同,不做考慮;PDCB是所有超出消隱門限的強(qiáng)脈沖RFI的凈占空比,對(duì)于有多個(gè)脈沖RFI干擾源的情況,PDCB的計(jì)算公式如下:

    wxdh4-gs5.gif

其中,PDCB,j是第j個(gè)干擾源超出消隱門限的強(qiáng)脈沖RFI占空比。

3.2.2 脈沖飽和接收機(jī)

    另有一些接收機(jī)沒有采用脈沖消隱的方式,其射頻前端將會(huì)被干擾源的強(qiáng)脈沖造成短暫飽和,稱為脈沖飽和接收機(jī)。

    對(duì)于這一類接收機(jī),其相關(guān)器輸出處的等效噪聲功率譜密度為:

    wxdh4-gs6.gif

其中,N0、I0,WB、PDCLIM、RI與脈沖消隱的情況基本相同,對(duì)于RI的閾值變?yōu)轱柡烷撝?對(duì)應(yīng)于消隱閾值),對(duì)于PDCLIM變?yōu)檩斎腼柡碗娖?對(duì)應(yīng)于輸入壓縮電平);NLIM是接收機(jī)A/D飽和電平與AGC的1 σ噪聲電壓之比,該值與A/D實(shí)現(xiàn)方式有關(guān),對(duì)于硬限制的1 bit接收機(jī)來說,NLIM=1。

3.2.3 將高軌SAR作為加性脈沖RFI

    如前所述,本文將高軌SAR信號(hào)作為加性干擾引入干擾基準(zhǔn)場(chǎng)景當(dāng)中,當(dāng)高軌SAR信號(hào)存在時(shí),對(duì)于脈沖消隱接收機(jī),等效載噪比的變化為:

     wxdh4-gs7.gif

其中,N0,EFF+SAB是添加了高軌SAR信號(hào)后的等效載噪比,PDCSAR是高軌SAR信號(hào)等效占空比,RSAR是低于消隱器閾值的高軌SAR信號(hào)功率譜密度與接收機(jī)噪聲溫度之比。PDCSAR與RSAR都是以接收機(jī)輸入壓縮電平作為功率參考點(diǎn)進(jìn)行估算的。

    同理,對(duì)于脈沖飽和接收機(jī),載噪比的變化為:

wxdh4-gs8-12.gif

其中,τPW是發(fā)射脈寬(s);BWoverlap是高軌SAR信號(hào)帶寬與GNSS接收機(jī)RF/IF濾波器帶寬重合部分的帶寬(MHz);BWchirp是高軌SAR信號(hào)發(fā)射帶寬(MHz)。當(dāng)公式(11)中的τPW,EFF為零時(shí),PDCSAR=0。

    將表1和表2中的相關(guān)參數(shù)帶入公式(9)及公式(10)可得接收機(jī)等效噪聲功率的變化值。不同脈沖重復(fù)頻率情況下,兩種接收機(jī)的等效噪聲功率變化值如圖3所示。

wxdh4-t3.gif

    從圖3中可以看出,脈沖消隱接收機(jī)的等效噪聲功率譜密度的變化為相同參數(shù)時(shí)脈沖飽和接收機(jī)的一半,且兩者變化量都不大,在0.4 dB以內(nèi)。

    根據(jù)得到的等效噪聲功率變化值,在信號(hào)功率不變的情況下,可以得到增加了高軌SAR信號(hào)后的等效載噪比:

    wxdh4-gs13.gif

其中,C/N0是未加入高軌SAR信號(hào)情況下的載噪比,ΔN0,EFF根據(jù)接收機(jī)類型進(jìn)行選擇。

    射頻干擾對(duì)GNSS接收機(jī)測(cè)距精度產(chǎn)生影響,主要體現(xiàn)在碼跟蹤環(huán)路的熱噪聲顫動(dòng),而通過C/N0和C/N0,EFF即可計(jì)算出加入高軌SAR信號(hào)前后碼環(huán)熱噪聲均方誤差值,并進(jìn)一步確定測(cè)距精度變化。由于本例中ΔN0,EFF較小,未超過門限,對(duì)測(cè)距精度產(chǎn)生的影響也較小,故不作詳細(xì)討論。

4 軟件實(shí)現(xiàn)

    為了方便進(jìn)行參數(shù)修改,對(duì)不同參數(shù)的高軌SAR信號(hào)、不同參數(shù)的接收機(jī)進(jìn)行分析,明確每個(gè)參數(shù)對(duì)接收機(jī)影響程度的大小,根據(jù)上述分析流程編寫軟件,界面如圖4所示。

wxdh4-t4.gif

    通過軟件界面可以對(duì)高軌SAR參數(shù)、接收機(jī)參數(shù)進(jìn)行配置,給定相應(yīng)的門限,即可計(jì)算出當(dāng)前參數(shù)配置是否超過門限。

5 結(jié)論

    本文將高軌SAR信號(hào)作為GNSS接收機(jī)的脈沖射頻干擾,參考ITU相關(guān)標(biāo)準(zhǔn),使用“源-路徑-接收機(jī)”的分析方法,從干擾源、傳播路徑、接收機(jī)模型三個(gè)方面設(shè)置了干擾評(píng)估相關(guān)參數(shù),通過計(jì)算接收機(jī)功率、等效噪聲變化值評(píng)估高軌SAR信號(hào)對(duì)GNSS地面接收機(jī)的性能影響;最終通過軟件實(shí)現(xiàn)了多種不同參數(shù)的SAR信號(hào)對(duì)不同類型GNSS接收機(jī)影響的分析。

    理論推導(dǎo)及計(jì)算結(jié)果表明:第一,在脈沖的消隱/飽和閾值相同時(shí),脈沖干擾對(duì)使用脈沖消隱方式接收機(jī)的影響為其對(duì)使用脈沖飽和方式接收機(jī)影響的一半(dB);第二,按照當(dāng)前參數(shù)設(shè)計(jì)的高軌SAR信號(hào)對(duì)B3I接收機(jī)會(huì)造成一定程度的干擾,但對(duì)性能的影響不大,不會(huì)影響B(tài)3I接收機(jī)的正常使用。

參考文獻(xiàn)

[1] 張薇,楊思全,范一大,等.高軌SAR衛(wèi)星在綜合減災(zāi)中的應(yīng)用潛力和工作模式需求[J].航天器工程,2017(1):127-131.

[2] 何瑞珠,劉彥明.通用導(dǎo)航接收機(jī)的脈沖干擾評(píng)估方法[J].移動(dòng)通信,2017(11):40-43.

[3] Recommendation ITU-R M.1902.Characteristics and protection criteria for receiving earth stations in the radionavigation-satellite service (space-to-Earth) operating in the band 1215-1300 MHz[R].International Telecommunications Union,2012.

[4] Recommendation ITU-R M.2220-0.Calculation method to determine aggregate interference parameters of pulsed RF systems operating in and near the bands 1164-1215 MHz and 1215-1300 MHz that may impact radionavigation-satellite service airborne and ground-based receivers operating in those frequency bands[R].International Telecommunications Union,2011.

[5] Recommendation ITU-R M.2030.Evaluation method for pulsed interference from relevant radio sources other than in the radionavigation-satellite service to the radionavigation-satellite service systems and networks operating in the 1164-1215 MHz, 1215-1300 MHz and 1559-1610 MHz frequency bands[R].International Telecommunications Union,2012.

[6] Recommendation ITU-R RS.2311.Pulsed radio frequency signal impact measurements and possible mitigation techniques between Earth exploration-satellite service(active) systems and RNSS systems and networks in the band 1215-1300 MHz[R].International Telecommunications Union,2014.

[7] 汪海波.大時(shí)帶寬線性調(diào)頻信號(hào)產(chǎn)生方法與系統(tǒng)實(shí)現(xiàn)技術(shù)研究[D].長(zhǎng)沙:中國人民解放軍國防科學(xué)技術(shù)大學(xué),2002.

[8] 謝鋼.GPS原理與接收機(jī)設(shè)計(jì)[M].北京:電子工業(yè)出版社,2009.



作者信息:

趙冠先,趙思浩,崔曉偉

(清華大學(xué) 電子工程系,北京100084)

此內(nèi)容為AET網(wǎng)站原創(chuàng),未經(jīng)授權(quán)禁止轉(zhuǎn)載。
亚洲一区二区欧美_亚洲丝袜一区_99re亚洲国产精品_日韩亚洲一区二区
欧美日韩在线精品| 国产一区二区电影在线观看| 欧美一级久久久| 中日韩美女免费视频网站在线观看| 亚洲国产精品第一区二区| 欧美一区二区三区四区视频 | 亚洲久久视频| 亚洲精品国久久99热| 亚洲国产精品久久人人爱蜜臀| 久久国产精品久久久| 久久精品青青大伊人av| 亚洲大胆人体在线| 亚洲欧洲视频在线| 99ri日韩精品视频| 一本色道久久99精品综合| av成人免费| 亚洲午夜电影| 午夜在线观看欧美| 久久成人精品一区二区三区| 欧美中文字幕在线| 久久这里只有| 欧美黄污视频| 国产精品成人一区二区网站软件| 欧美性开放视频| 国产精品夜夜夜| 国产综合婷婷| 亚洲国产精品一区制服丝袜| 亚洲日本电影| 99视频+国产日韩欧美| 亚洲亚洲精品三区日韩精品在线视频| 亚洲一区美女视频在线观看免费| 欧美一级午夜免费电影| 亚洲黄色大片| 亚洲视频一区在线观看| 欧美一区二区视频在线| 久久夜色精品| 欧美日韩一区二区三区在线看| 国产精品日本一区二区| 好看的日韩av电影| 亚洲日本国产| 午夜国产精品影院在线观看 | 亚洲日本精品国产第一区| 一区二区日韩欧美| 欧美在线高清| 欧美韩日一区| 国产精品天天看| 亚洲成色www8888| 一区二区免费看| 久久精品免视看| 亚洲手机成人高清视频| 欧美在线高清| 欧美乱大交xxxxx| 国产亚洲精久久久久久| 亚洲激情成人| 午夜精品影院| 日韩一区二区精品在线观看| 午夜视频在线观看一区二区三区| 久久天天躁狠狠躁夜夜av| 欧美日韩国产大片| 狠狠久久亚洲欧美专区| 一区二区久久久久| 亚洲高清不卡在线| 亚洲免费婷婷| 欧美国产91| 国产亚洲人成a一在线v站| 亚洲日韩欧美视频| 欧美一二三视频| 在线中文字幕日韩| 另类av一区二区| 国产精品视频一区二区三区| 91久久精品国产| 久久成人这里只有精品| 亚洲一区二区高清视频| 麻豆成人在线播放| 国产精品视频网站| 亚洲老司机av| 最新日韩在线| 久久九九久久九九| 国产精品狼人久久影院观看方式| 尤物视频一区二区| 午夜精品一区二区三区在线视| 99天天综合性| 男女精品视频| 国产一区二区在线免费观看| 亚洲视频www| 日韩视频在线免费观看| 久久精品午夜| 国产精品日韩在线观看| 日韩视频免费在线| 亚洲乱码国产乱码精品精可以看 | 午夜国产一区| 亚洲欧美日韩国产成人| 欧美日韩二区三区| 亚洲福利视频在线| 久久精品99国产精品酒店日本| 亚洲欧美激情在线视频| 欧美日韩国产经典色站一区二区三区| 在线激情影院一区| 国产精品影音先锋| 亚洲精品视频一区| 亚洲国产高清在线观看视频| 欧美一区二区免费| 欧美视频四区| 亚洲国产精品久久久久秋霞不卡| 久久成人免费视频| 欧美亚洲综合另类| 欧美午夜精品久久久久久浪潮| 亚洲激情视频在线| 亚洲韩国精品一区| 久久一区二区三区国产精品| 国产视频一区在线观看| 亚洲欧美色婷婷| 欧美亚洲综合网| 国产精品夜色7777狼人| 亚洲一区日韩| 欧美一区二区国产| 免费成人美女女| 欧美日韩卡一卡二| 亚洲欧洲日韩综合二区| 亚洲国产一区二区三区在线播| 久久久久国产精品一区三寸| 国产精品少妇自拍| 中文精品在线| 亚洲欧美一区二区三区极速播放| 国产精品mm| 亚洲自拍都市欧美小说| 香蕉亚洲视频| 国产欧美日韩视频| 午夜在线视频观看日韩17c| 久久av红桃一区二区小说| 国产欧美日韩视频一区二区三区| 亚洲欧美视频在线观看| 久久精品91久久香蕉加勒比| 国产午夜一区二区三区| 午夜在线a亚洲v天堂网2018| 久久精品国产亚洲精品| 狠狠色香婷婷久久亚洲精品| 久久精品国产第一区二区三区最新章节| 久久久青草婷婷精品综合日韩 | 欧美成人一区二区在线| 亚洲人在线视频| 亚洲天堂av综合网| 国产精品视频一| 久久国产精品亚洲77777| 美女视频网站黄色亚洲| 亚洲日本黄色| 亚洲视屏一区| 国产欧美一区二区色老头 | 在线观看视频免费一区二区三区| 亚洲乱码国产乱码精品精天堂 | 亚洲韩国一区二区三区| 中日韩美女免费视频网址在线观看| 国产精品久久一区二区三区| 午夜亚洲精品| 欧美激情一区二区三区在线视频 | 午夜精品福利一区二区蜜股av| 久久综合给合| 日韩午夜在线电影| 久久gogo国模裸体人体| 亚洲国产精品尤物yw在线观看 | 欧美精品国产一区二区| 亚洲深夜福利网站| 久久久亚洲精品一区二区三区| 亚洲高清精品中出| 亚洲欧美日韩精品| 狠狠色综合网| 亚洲午夜国产成人av电影男同| 国产精品视频一区二区高潮| 亚洲高清视频在线| 欧美视频在线观看视频极品| 欧美一区二区三区免费观看视频 | 欧美亚洲一区二区在线| 欧美高清视频在线播放| 亚洲综合国产激情另类一区| 久久婷婷久久| 一本色道综合亚洲| 久久综合婷婷| 亚洲一区二区在线免费观看| 六月丁香综合| 亚洲一区国产视频| 欧美黑人国产人伦爽爽爽| 亚洲欧美春色| 欧美日韩mv| 欧美在线播放一区| 国产精品av免费在线观看| 亚洲福利视频一区二区| 国产精品地址| 亚洲精品国产精品久久清纯直播 | 欧美日韩精品免费| 欧美呦呦网站| 国产精品久久999| 亚洲激情影视| 国产丝袜美腿一区二区三区| 99热精品在线| 精品1区2区3区4区| 欧美一级片久久久久久久| 亚洲免费高清| 乱中年女人伦av一区二区| 亚洲在线一区二区三区| 欧美日韩高清在线|