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

    基于GNSS的武漢區(qū)域電離層TEC建模

    2022-03-09 07:16:34吳偉銓周晨禹智斌左平兵趙家奇趙正予
    電波科學(xué)學(xué)報 2022年1期
    關(guān)鍵詞:電離層級數(shù)建模

    吳偉銓 周晨 禹智斌 左平兵 趙家奇 趙正予

    (1. 武漢大學(xué)電子信息學(xué)院空間物理系,武漢 430072;2. 哈爾濱工業(yè)大學(xué)空間科學(xué)與應(yīng)用技術(shù)研究院,深圳 518055)

    引 言

    電離層天頂方向總電子含量(vertical total electron content, VTEC)是衛(wèi)星導(dǎo)航系統(tǒng)的重要產(chǎn)品. 近年來不少優(yōu)秀的學(xué)者利用經(jīng)典的電離層模型,對VTEC進(jìn)行了建模研究. Lanyi等[1]提出利用2階或3階多項(xiàng)式模型計算單站VTEC;Schaer等[2]利用球諧級數(shù)模型進(jìn)行全球電離層建模;Yuan等[3-4]提出廣義三角級數(shù)函數(shù)模型,該模型在小區(qū)域VTEC 建模中有很好的效果. 在各種電離層模型被提出之后,又有不少專家學(xué)者對這些模型進(jìn)行了比較分析. 文獻(xiàn)[5]對多項(xiàng)式模型、三角級數(shù)模型和低階球諧級數(shù)模型之間的一致性和參數(shù)設(shè)置局限性進(jìn)行了分析,指出了它們的優(yōu)缺點(diǎn). 文獻(xiàn)[6]根據(jù)應(yīng)用環(huán)境對電離層模型進(jìn)行了分類,并提出了模型的一些問題和改進(jìn)措施. 文獻(xiàn)[7]指出了多項(xiàng)式模型和三角級數(shù)模型不適用于大區(qū)域的電離層建模. 文獻(xiàn)[8-9]表明在模型參數(shù)選取適當(dāng)?shù)那闆r下,球諧級數(shù)模型適用于各種大小區(qū)域的電離層建模. 文獻(xiàn)[10]基于球諧級數(shù)模型對山東區(qū)域進(jìn)行電離層建模,并采用單頻精密單點(diǎn)定位驗(yàn)證了該區(qū)域電離層模型的有效性. 電離層在不同高度上的電子密度會有很大差異, 文獻(xiàn)[11]提出了一種三層的電離層模型,可以反映不同高度區(qū)域的VTEC情況. 在電離層建模中,通常要采用一定的估計策略來對模型系數(shù)進(jìn)行估計. 常采用的是分段常數(shù)(piece-wise constant, PWC)估計策略和分段線性(piece-wise linear, PWL)估計策略. 文獻(xiàn)[12]指出,只有在低緯度地區(qū)的電離層模型改正精度才會對估計策略比較敏感,而這種影響是隨緯度升高而減小的. 文獻(xiàn)[13]在總結(jié)全球?qū)Ш叫l(wèi)星系統(tǒng)(Global Navigation Satellite System, GNSS)電離層延遲修正成果的基礎(chǔ)上,從建模及誤差分析等方面介紹了我國在這一領(lǐng)域的研究進(jìn)展. 但是,上述文獻(xiàn)都沒有對電離層建模方法和估計策略進(jìn)行綜合對比分析,也沒有研究地磁活動對電離層建模的影響. 因此,本文在上述研究的基礎(chǔ)上,對武漢地區(qū)進(jìn)行電離層建模,目的是比較不同建模方法和估計策略在該區(qū)域中的適用性,以及探究地磁活動和電離層建模的相關(guān)性.

    國際GNSS服務(wù)(International GNSS Service,IGS)的電離層聯(lián)合分析中心(Ionosphere Associate Analysis Centers, IAACs)可以提供高精度的電離層產(chǎn)品,主要包括歐洲軌道確定中心(Center for Orbit Determination in Europe, CODE)[8]、美國噴氣推進(jìn)實(shí)驗(yàn)室(Jet Propulsion Laboratory, JPL)[14]、歐洲太空局(European Space Agency/European Space Operations Center, ESA/ESOC)[15]、西班牙加泰羅尼亞理工大學(xué)(Universitat Politècnica de Catalunya, UPC)[16]、加拿大自然資源部門(Natural Resources Canada, NRCan)、中國科學(xué)院測量與地球物理研究所和光電研究院(Chinese Academy of Science, CAS)[17]及武漢大學(xué)衛(wèi)星導(dǎo)航與定位技術(shù)研究中心(Wuhan University, WHU)[18].本文選擇CODE提供的電離層產(chǎn)品作為參考源,其采用15階球諧級數(shù)進(jìn)行全球的電離層建模,能提供時間分辨率為1 h、空間分辨率為2.5°×5°的全球電離層格網(wǎng)地圖(global ionospheric map, GIM)[19]. 選取以武漢為中心的十七個連續(xù)運(yùn)行參考站(continuous operation reference system, CORS)在不同地磁活動情況下的數(shù)據(jù)進(jìn)行試驗(yàn),并與CODE公布的GIM作對比,對于分析建立武漢地區(qū)的電離層模型、削弱電離層延遲誤差具有重要意義[10].

    傳統(tǒng)的二維電離層模型是區(qū)域電離層建模的經(jīng)典手段,但是,它們只能根據(jù)經(jīng)度和緯度擬合區(qū)域內(nèi)的VTEC,卻不能反映高度上的電子密度變化情況.而通過三層電離層模型可以從經(jīng)度、緯度、高度三個維度上對電離層進(jìn)行建模,并且可以有效提高電離層地圖的擬合精度[11]. 此外,多層電離層模型對空間目標(biāo)監(jiān)視和低軌衛(wèi)星航天測控具有重要應(yīng)用價值.本文通過研究三種二維電離層模型和估計策略對VTEC計算精度的影響,探索適用于武漢區(qū)域的電離層模型、模型參數(shù)和估計策略,以及分析模型精度和地磁活動的相關(guān)性. 在此研究的基礎(chǔ)上,進(jìn)一步通過三層電離層模型,來分析不同高度區(qū)域的VTEC情況.

    1 電離層模型

    二維電離層模型包括多項(xiàng)式模型、三角級數(shù)模型、球諧級數(shù)模型、樣條函數(shù)模型、廣義三角級數(shù)模型. 下面主要介紹多項(xiàng)式模型和球諧級數(shù)模型,以及在二維電離層模型基礎(chǔ)上提出的三層電離層模型.

    1.1 多項(xiàng)式模型

    多項(xiàng)式模型[1]將VTEC看作如下表達(dá)式:

    式中:Tv(φ,S)是電離層穿刺點(diǎn)(ionosphere piece point,IPP)處的VTEC;m和n是 模型階數(shù);Eij是模型系數(shù);φ是地理緯度; φ0是IPP覆蓋中心的地理緯度;S是太陽時角;S0是IPP覆蓋中心的太陽時角.

    太陽時角S具有如下形式:

    式中:θ是地理經(jīng)度,單位是(°);t是UTC時間,單位是h.

    通常求解一個多項(xiàng)式模型用2~3 h,所以要每2~3 h計算一組模型系數(shù).

    1.2 球諧級數(shù)模型

    球諧級數(shù)模型[2]的表達(dá)式如下:

    式中:n是模型階數(shù);(sinφ) 是i度j階的締和勒讓德函數(shù);φ是IPP的緯度;s是IPP的經(jīng)度;Aij和Bij是模型系數(shù).

    求解球諧級數(shù)模型也是使用2~3 h. 區(qū)域電離建模時,可以采用低階的球諧級數(shù)模型.

    1.3 三層電離層模型

    三層電離層模型[11]把總電子含量(total electron content, TEC)看成與經(jīng)度θ、緯度φ和高度h都相關(guān)的函數(shù),其表達(dá)式如下:

    式中:T是發(fā)送機(jī)到接收機(jī)路徑上的TEC;Rtra和Rrec分別是發(fā)送機(jī)和接收機(jī)的位置;Ne(θ,φ,h)是電子密度函數(shù). 即TEC是電子密度沿發(fā)送機(jī)和接收機(jī)路徑的積分. 電子密度函數(shù)可看成水平基函數(shù)和垂直基函數(shù)的線性組合,表達(dá)式如下:

    式中:cij是 參數(shù);Hi是水平方向的電子密度基函數(shù),可沿用在二維電離層模型中使用的函數(shù);Vj是垂直方向的電子密度基函數(shù),可采用國際參考電離層(international reference ionosphere, IRI)模型[20]輸出的電子密度作為該函數(shù);NH和NV是兩個基函數(shù)的階數(shù).

    三層電離層模型把VTEC看作如下表達(dá)式:

    式中:Tv(φ,S,h1,h2)代 表緯度φ、 太陽時角S、高度h1到h2之 間的VTEC;Tv(φ,S,0,∞)是二維電離層模型計算的結(jié)果;cj是參數(shù). 即通過垂直方向電子密度函數(shù)的積分可以得到與高度、緯度、經(jīng)度都相關(guān)的VTEC.

    2 武漢地區(qū)電離層建模

    本文對武漢周邊17個CORS 2016-01-01—05的數(shù)據(jù)進(jìn)行建模,如圖1所示. 圖2是2015-12-31—2016-01-07的地磁活動指數(shù)(Dst和Kp). 可以看到,01-01處在磁暴恢復(fù)相期,而01-05處在地磁平靜時期,選取這兩天數(shù)據(jù)有利于研究不同地磁活動水平對區(qū)域電離層建模的影響.

    圖1 武漢周邊17個CORS分布圖Fig. 1 Distribution map of 17 continuously operating reference stations around Wuhan

    圖2 2015-12-31—2016-01-07 Dst指數(shù)和Kp指數(shù)Fig. 2 Dst index and Kp index from Dec. 31st,2015 to Jan. 7th, 2016

    2.1 建模方法

    本文基于GNSS觀測數(shù)據(jù),使用載波相位平滑偽距法[14]提高偽距觀測精度,衛(wèi)星截止仰角設(shè)置為15°,投影函數(shù)采用標(biāo)準(zhǔn)投影函數(shù),時間采樣率為30 s,電離層薄層高度設(shè)置為450 km.

    本文利用參考文獻(xiàn)[1]中提出的2階多項(xiàng)式模型、文獻(xiàn)[5]中使用的4階球諧級數(shù)模型,以及文獻(xiàn)[21]中使用的4×3階多項(xiàng)式模型這三種模型進(jìn)行區(qū)域VTEC建模. 估計策略分別采用PWC估計策略和PWL估計策略[12]. PWC假定VTEC模型系數(shù)在適用時間內(nèi)不變;PWL假定模型系數(shù)在這段時間線性變化,通常需要兩組系數(shù),分別代表該時間段開始和結(jié)束時的模型系數(shù),則任意時刻的模型系數(shù)都能由這兩組系數(shù)線性表示. 將VTEC模型系數(shù)和硬件延遲偏差(differential code biases, DCB)聯(lián)合建立法方程[22],采用最小二乘法,每2 h計算一組系數(shù). 得到聯(lián)合DCB后,通過假設(shè)所有衛(wèi)星的DCB之和為零來分離接收機(jī)和衛(wèi)星DCB. 按照0.1°的間隔劃分網(wǎng)格點(diǎn),求得模型系數(shù)后再反代入網(wǎng)格點(diǎn),計算出各網(wǎng)格點(diǎn)處VTEC.

    2.2 結(jié)果和分析

    對2016-01-01的數(shù)據(jù)進(jìn)行建模,采用PWC估計策略時,4×3階多項(xiàng)式模型、2階多項(xiàng)式模型、4階球諧級數(shù)模型的VTEC建模結(jié)果如圖3所示.

    圖3 2016-01-01三種模型采用PWC估計策略的VTEC建模結(jié)果Fig. 3 Modeling diagram of three models and piece-wise constant strategy on Jan. 1st, 2016

    對2016-01-05的數(shù)據(jù)進(jìn)行建模,采用PWC估計策略時,4×3階多項(xiàng)式模型、2階多項(xiàng)式模型、4階球諧級數(shù)模型的VTEC建模結(jié)果如圖4所示.

    圖4 2016-01-05三種模型采用PWC估計策略的VTEC建模結(jié)果Fig. 4 Modeling diagram of three models and piece-wise constant strategy on Jan. 5th, 2016

    從圖3~4的二維電離層模型結(jié)果可以看出:在0:00—6:00UT時,電離層電子濃度逐漸增大,并在6:00UT時達(dá)到峰值,而在16:00—24:00UT電子濃度都偏??;低緯度地區(qū)的電子濃度比高緯度地區(qū)高;發(fā)生磁暴時,電離層電子濃度顯著增大,說明地磁活動會影響電離層的電子濃度.

    鑒于CODE在硬件延遲及電離層解算方面的優(yōu)越性,采用CODE公布的GIM作為外部對比源. CODE提供的該區(qū)域2016-01-01和01-05的VTEC結(jié)果如圖5所示.

    圖5 2016-01-01和01-05 CODE公布的VTEC建模圖Fig. 5 Modeling diagram published by CODE on Jan. 1st and Jan. 5th, 2016

    通過三種模型的VTEC建模圖和CODE公布的建模圖進(jìn)行對比可以發(fā)現(xiàn),三種模型都能達(dá)到較好的擬合效果. 為了進(jìn)一步分析這些模型在地磁平靜時期和地磁擾動時期的建模精度變化,對武漢單站VTEC在PWC估計策略下的計算結(jié)果和CODE公布值進(jìn)行了比較,結(jié)果如圖6所示.

    圖6 2016-01-01和01-05武漢單站VTEC建模結(jié)果和CODE公布結(jié)果對比Fig. 6 Comparison of VTEC modeling results and code results of Wuhan single station on Jan. 1st and Jan. 5th, 2016

    從圖6可以看出,三種模型的建模結(jié)果幾乎重疊在一起,且和CODE公布的VTEC變化趨勢基本一致,但VTEC大小依然存在一定的誤差. 原因是CODE進(jìn)行的是全球電離層擬合,在中國區(qū)域用的站不到十個,在小區(qū)域范圍內(nèi)擬合精度不夠.

    鑒于此原因,本文將VTEC通過電離層投影函數(shù)轉(zhuǎn)化成斜向總電子含量(slant total electron content,STEC),與實(shí)測值比較計算30天的殘差值均方根(root mean square, RMS),以此來分析模型和估計策略的精度. 其中殘差值的表達(dá)式如下:

    式中:R是殘差值;TS是STEC的觀測值;fm是電離層投影函數(shù),定義為STEC和VTEC的比;TV是VTEC;D是DCB.

    統(tǒng)計2016-01-01—30的STEC殘差值RMS,在PWL和PWC估計策略下的結(jié)果如圖7所示.

    圖7 兩種估計策略下的RMS分布Fig. 7 RMS distribution with two estimation strategies

    從圖7的殘差RMS結(jié)果可以看出:無論采用哪種估計策略,效果最好的是4階球諧級數(shù)模型,其次是4×3階多項(xiàng)式模型,而2階多項(xiàng)式模型的效果最差,這與其模型參數(shù)較少有關(guān);三種模型在PWL估計策略下的效果要優(yōu)于PWC估計策略,說明PWL估計策略更適用于區(qū)域的VTEC建模. 當(dāng)使用PWL估計策略和4階球諧級數(shù)模型時,RMS集中在0~1.5 TECU,區(qū)域VTEC建模達(dá)到較高的精度.

    對比01-01和01-05的殘差,發(fā)現(xiàn)01-05三種模型的RMS都比01-01對應(yīng)的RMS要小. 另外,從圖6可以看出,同一時刻,模型的計算結(jié)果會比CODE公布值小,除了在01-01T6:00UT,此時VTEC也恰好達(dá)到峰值. 從以上的分析可以推斷出,地磁擾動會對VTEC計算結(jié)果產(chǎn)生干擾. 一方面,地磁活動強(qiáng)的時候GNSS觀測誤差會增大;另一方面,VTEC模型的計算精度也會降低.

    谷子生育期為80~140 d,適合小麥、大麥等收后復(fù)種。筆者在內(nèi)蒙古農(nóng)牧業(yè)科學(xué)院試驗(yàn)地曾進(jìn)行飼用谷子不同生長時期全株?duì)I養(yǎng)物質(zhì)測定試驗(yàn),結(jié)果表明:初花期預(yù)干全株的蛋白質(zhì)含量要明顯高于成熟期的蛋白質(zhì)含量,而初花期預(yù)干全株的纖維含量要明顯低于成熟期纖維含量,初花期預(yù)干全株的干物質(zhì)量與成熟期相當(dāng),所以飼用谷子在初花期就可以進(jìn)行收獲,大大縮短了生長周期,是短日期作物收獲后復(fù)種的優(yōu)良作物之一。

    以上都是基于經(jīng)度和緯度的二維電離層建模.然而,二維電離層模型不能反映VTEC與高度之間的關(guān)系. 為了進(jìn)一步研究底部電離層、峰高電離層和頂部電離層的VTEC情況,本文在上述研究的基礎(chǔ)上利用三層電離層模型對武漢地區(qū)進(jìn)行建模. 參考文獻(xiàn)[11]中對電離層高度的分層方法,將電離層分為三個區(qū)域:90~350 km,350~600 km,600~1 500 km. 利用IRI模型計算出每個IPP三個高度區(qū)域的VTEC比例,再代入4階球諧級數(shù)模型在PWC估計策略下的計算結(jié)果中,構(gòu)建出三層電離層模型.2016-01-01和01-05的建模結(jié)果如圖8所示.

    圖8 2016-01-01和01-05三層電離層模型建模Fig. 8 Modeling diagram of three layer ionosphere model on Jan. 5th, 2016

    從圖8可以看出:90~350 km高度區(qū)間的電離層TEC最大,其次是350~600 km,而在600 km以上,TEC會變得很稀少. 在得出該結(jié)果后,本文進(jìn)一步研究三層電離層VTEC隨時間的變化規(guī)律. 利用2015-12-31—2016-01-06七天的數(shù)據(jù),對武漢上空的VTEC進(jìn)行了計算,得出三個高度區(qū)間以及總的VTEC結(jié)果如圖9所示.

    從圖9可以看出:在每天VTEC達(dá)到峰值時刻前后,底部電離層的VTEC要明顯大于峰高電離層和頂部電離層的VTEC. 當(dāng)VTEC在谷值時,底部VTEC和峰高VTEC較為接近. 同時,不同高度的VTEC對地磁活動變化的響應(yīng)也有明顯區(qū)別,底部TEC和峰高TEC受地磁活動的影響較大.

    圖9 2015-12-31—2016-01-06武漢上空三層電離層VTEC結(jié)果Fig. 9 VTEC results of three layers over Wuhan from Dec.31st, 2015 to Jan. 6th, 2016

    3 結(jié) 論

    本文針對區(qū)域電離層建模的三種模型和兩種估計策略,分別選取地磁平靜和地磁擾動的兩天數(shù)據(jù)進(jìn)行計算,并以CODE公布的GIM和STEC實(shí)測值作為對比源進(jìn)行分析,得出一些結(jié)論.

    1)在區(qū)域VTEC建模中,三種二維電離層模型的建模結(jié)果和CODE公布的GIM基本一致,達(dá)到較好的擬合精度.

    2)地磁擾動會對模型的擬合結(jié)果產(chǎn)生干擾,降低模型的擬合精度.

    4)當(dāng)使用PWC估計策略時,所有模型的擬合效果都會變差. 因此,在該區(qū)域VTEC建模中建議使用PWL估計策略.

    5)從三層電離層模型結(jié)果可以看出,電離層的電子含量集中在90~600 km高度區(qū)間,不同高度的VTEC對地磁活動變化的響應(yīng)有很大區(qū)別.

    文中的試驗(yàn)是基于武漢周邊地區(qū)兩天的數(shù)據(jù)來進(jìn)行的,沒有考慮不同的地區(qū)和更長時間的情況. 而且,模型的精度受實(shí)際應(yīng)用場景和觀測數(shù)據(jù)密集程度等多種因素的影響,為了達(dá)到詳細(xì)的分析對比,需要進(jìn)行更加完善的試驗(yàn)設(shè)計.

    猜你喜歡
    電離層級數(shù)建模
    一種電離層TEC格點(diǎn)預(yù)測模型
    Kalman濾波估算電離層延遲的一種優(yōu)化方法
    聯(lián)想等效,拓展建?!浴皫щ娦∏蛟诘刃鲋凶鰣A周運(yùn)動”為例
    Dirichlet級數(shù)及其Dirichlet-Hadamard乘積的增長性
    基于PSS/E的風(fēng)電場建模與動態(tài)分析
    電子制作(2018年17期)2018-09-28 01:56:44
    不對稱半橋變換器的建模與仿真
    電離層對中高軌SAR影響機(jī)理研究
    幾個常數(shù)項(xiàng)級數(shù)的和
    p級數(shù)求和的兩種方法
    Linux Shell語言在電離層解算中的應(yīng)用
    香蕉国产在线看| 亚洲精品美女久久av网站| 波多野结衣av一区二区av| 男人舔女人的私密视频| 性色av一级| 少妇人妻精品综合一区二区| 女性被躁到高潮视频| 一本久久精品| 日韩 亚洲 欧美在线| 黄色怎么调成土黄色| 伦理电影免费视频| 天天躁日日躁夜夜躁夜夜| 天堂俺去俺来也www色官网| av在线老鸭窝| 久久人妻熟女aⅴ| 99热国产这里只有精品6| 制服人妻中文乱码| 久久久久精品性色| 日韩精品免费视频一区二区三区| 亚洲精品aⅴ在线观看| 少妇人妻 视频| 久久青草综合色| 人妻少妇偷人精品九色| 国产免费又黄又爽又色| 亚洲伊人色综图| 熟女少妇亚洲综合色aaa.| 中文字幕av电影在线播放| 成人国语在线视频| av免费在线看不卡| 午夜福利在线观看免费完整高清在| 91午夜精品亚洲一区二区三区| 99久久中文字幕三级久久日本| xxxhd国产人妻xxx| 久久久精品区二区三区| 亚洲,欧美精品.| 交换朋友夫妻互换小说| 日韩中文字幕视频在线看片| 国产精品久久久久久av不卡| 涩涩av久久男人的天堂| 久久精品人人爽人人爽视色| 在线观看免费日韩欧美大片| 精品一区在线观看国产| 免费在线观看视频国产中文字幕亚洲 | 18禁观看日本| 99re6热这里在线精品视频| 国产极品天堂在线| 伦理电影免费视频| 人妻系列 视频| 另类亚洲欧美激情| 街头女战士在线观看网站| 可以免费在线观看a视频的电影网站 | 中国国产av一级| 日韩一区二区三区影片| 精品一区二区三卡| 久久99一区二区三区| 91国产中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 国产人伦9x9x在线观看 | 亚洲av欧美aⅴ国产| av又黄又爽大尺度在线免费看| 欧美日韩国产mv在线观看视频| 久久人人爽人人片av| 青春草视频在线免费观看| 国产熟女欧美一区二区| 高清av免费在线| videossex国产| 婷婷色av中文字幕| 1024香蕉在线观看| 国产精品成人在线| 少妇猛男粗大的猛烈进出视频| 婷婷成人精品国产| 黄色配什么色好看| 日本午夜av视频| 青春草亚洲视频在线观看| 国产精品嫩草影院av在线观看| 国产亚洲av片在线观看秒播厂| 成年人免费黄色播放视频| 两性夫妻黄色片| 国产精品 欧美亚洲| 久久午夜福利片| 美女视频免费永久观看网站| 尾随美女入室| 一区二区三区激情视频| 99热全是精品| 综合色丁香网| 亚洲精品久久午夜乱码| 免费久久久久久久精品成人欧美视频| 国产毛片在线视频| 一级黄片播放器| 最近2019中文字幕mv第一页| 国产精品一二三区在线看| 精品99又大又爽又粗少妇毛片| 天天躁狠狠躁夜夜躁狠狠躁| 午夜福利一区二区在线看| 看十八女毛片水多多多| 国产女主播在线喷水免费视频网站| 制服诱惑二区| 亚洲综合色网址| 黄色毛片三级朝国网站| 97人妻天天添夜夜摸| 日日啪夜夜爽| 免费黄色在线免费观看| 在线观看一区二区三区激情| 永久免费av网站大全| 日日摸夜夜添夜夜爱| 免费日韩欧美在线观看| 久久久国产欧美日韩av| 亚洲av欧美aⅴ国产| 日韩av免费高清视频| 热99国产精品久久久久久7| 国产亚洲av片在线观看秒播厂| 亚洲欧美一区二区三区黑人 | 999精品在线视频| 99re6热这里在线精品视频| 五月伊人婷婷丁香| 精品少妇内射三级| 国产一区有黄有色的免费视频| 亚洲第一区二区三区不卡| 日韩大片免费观看网站| 国产 精品1| 国产又色又爽无遮挡免| 最新中文字幕久久久久| 丝袜脚勾引网站| 热99久久久久精品小说推荐| 精品国产乱码久久久久久小说| 中文字幕av电影在线播放| 十八禁高潮呻吟视频| 日本免费在线观看一区| 成人影院久久| 日韩制服丝袜自拍偷拍| 久久免费观看电影| 啦啦啦中文免费视频观看日本| 亚洲美女搞黄在线观看| av网站免费在线观看视频| www.熟女人妻精品国产| 晚上一个人看的免费电影| 亚洲经典国产精华液单| 午夜激情久久久久久久| 国产成人欧美| 日本-黄色视频高清免费观看| av网站在线播放免费| 一级a爱视频在线免费观看| 亚洲国产精品999| 69精品国产乱码久久久| 国产精品久久久久久av不卡| av视频免费观看在线观看| 男人操女人黄网站| 日日爽夜夜爽网站| 日韩电影二区| 最黄视频免费看| 亚洲欧美一区二区三区久久| 久久精品熟女亚洲av麻豆精品| 日韩av不卡免费在线播放| 久久午夜综合久久蜜桃| 日日啪夜夜爽| 麻豆av在线久日| 久久久精品94久久精品| 中文字幕人妻熟女乱码| 国产精品免费视频内射| 日韩大片免费观看网站| 香蕉国产在线看| 国产 一区精品| 伦理电影免费视频| www日本在线高清视频| 久久久久久久大尺度免费视频| 亚洲欧美色中文字幕在线| 成人免费观看视频高清| 亚洲精品,欧美精品| 男女高潮啪啪啪动态图| 丰满少妇做爰视频| 欧美亚洲日本最大视频资源| 80岁老熟妇乱子伦牲交| 免费观看在线日韩| 日本vs欧美在线观看视频| 色吧在线观看| 超碰成人久久| 亚洲精华国产精华液的使用体验| 青草久久国产| 欧美日韩国产mv在线观看视频| 欧美97在线视频| 高清黄色对白视频在线免费看| 成人影院久久| 丝瓜视频免费看黄片| 精品福利永久在线观看| 在线观看免费高清a一片| 午夜影院在线不卡| 日韩人妻精品一区2区三区| 国产视频首页在线观看| 国产精品久久久久久精品古装| 国产白丝娇喘喷水9色精品| 国产一区二区在线观看av| 人妻少妇偷人精品九色| 亚洲经典国产精华液单| 中文字幕最新亚洲高清| 国产爽快片一区二区三区| 久久人妻熟女aⅴ| 国产精品欧美亚洲77777| 国产片特级美女逼逼视频| 一个人免费看片子| 高清黄色对白视频在线免费看| 午夜免费鲁丝| av在线app专区| 久久99一区二区三区| 色视频在线一区二区三区| 美女视频免费永久观看网站| 亚洲中文av在线| 日韩一卡2卡3卡4卡2021年| 超碰成人久久| 中文字幕人妻丝袜一区二区 | 国产亚洲最大av| 成年av动漫网址| 人体艺术视频欧美日本| 人人妻人人添人人爽欧美一区卜| 成年女人在线观看亚洲视频| 色网站视频免费| 国产精品久久久久久精品古装| 久久97久久精品| √禁漫天堂资源中文www| 国产男人的电影天堂91| xxxhd国产人妻xxx| videossex国产| 欧美精品亚洲一区二区| 精品少妇一区二区三区视频日本电影 | 美女中出高潮动态图| 日韩熟女老妇一区二区性免费视频| 69精品国产乱码久久久| 国精品久久久久久国模美| 免费看av在线观看网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一级毛片电影观看| av女优亚洲男人天堂| 午夜福利乱码中文字幕| 晚上一个人看的免费电影| 两个人免费观看高清视频| 日韩在线高清观看一区二区三区| 欧美日韩综合久久久久久| 亚洲av.av天堂| 亚洲精品在线美女| 大话2 男鬼变身卡| 日韩欧美精品免费久久| 亚洲久久久国产精品| 欧美日韩一级在线毛片| 水蜜桃什么品种好| 久久久国产精品麻豆| av网站在线播放免费| 在线观看www视频免费| 成年美女黄网站色视频大全免费| 王馨瑶露胸无遮挡在线观看| 午夜福利一区二区在线看| 中文字幕人妻熟女乱码| 侵犯人妻中文字幕一二三四区| 90打野战视频偷拍视频| 日韩成人av中文字幕在线观看| 国产精品久久久av美女十八| 免费观看a级毛片全部| 欧美日韩精品网址| 女性生殖器流出的白浆| 亚洲伊人久久精品综合| 秋霞在线观看毛片| 黄频高清免费视频| 大片免费播放器 马上看| 国产伦理片在线播放av一区| 亚洲精品乱久久久久久| 国产精品久久久久久久久免| 美女大奶头黄色视频| 亚洲婷婷狠狠爱综合网| kizo精华| 男女国产视频网站| 亚洲欧美色中文字幕在线| 99热国产这里只有精品6| 亚洲精品成人av观看孕妇| 免费高清在线观看日韩| 国产日韩一区二区三区精品不卡| 一级爰片在线观看| 免费观看a级毛片全部| 国产黄色视频一区二区在线观看| 亚洲欧美一区二区三区黑人 | 纵有疾风起免费观看全集完整版| 一级毛片我不卡| 久久亚洲国产成人精品v| 可以免费在线观看a视频的电影网站 | av在线老鸭窝| 99香蕉大伊视频| 午夜激情av网站| 精品国产露脸久久av麻豆| 我的亚洲天堂| 1024香蕉在线观看| 激情五月婷婷亚洲| 日韩电影二区| 大片电影免费在线观看免费| 日本91视频免费播放| 一区在线观看完整版| 亚洲欧美日韩另类电影网站| 一级,二级,三级黄色视频| 国语对白做爰xxxⅹ性视频网站| 亚洲精品美女久久久久99蜜臀 | 国产综合精华液| 99久久综合免费| 九九爱精品视频在线观看| 考比视频在线观看| 亚洲精品,欧美精品| 亚洲伊人久久精品综合| 人人妻人人添人人爽欧美一区卜| 婷婷色综合www| 午夜福利乱码中文字幕| av在线观看视频网站免费| av国产精品久久久久影院| 国产野战对白在线观看| 亚洲激情五月婷婷啪啪| 我的亚洲天堂| 亚洲伊人色综图| 纵有疾风起免费观看全集完整版| 看免费av毛片| 老熟女久久久| 不卡视频在线观看欧美| 亚洲国产av新网站| 老司机影院成人| 少妇的逼水好多| 欧美日韩精品成人综合77777| 9191精品国产免费久久| 成年女人在线观看亚洲视频| 亚洲色图 男人天堂 中文字幕| 好男人视频免费观看在线| 久久婷婷青草| 欧美日韩一级在线毛片| 又黄又粗又硬又大视频| 成人亚洲精品一区在线观看| 99久久综合免费| 欧美激情极品国产一区二区三区| 亚洲精品久久午夜乱码| 国产在线视频一区二区| 亚洲,一卡二卡三卡| 99re6热这里在线精品视频| 久久久久人妻精品一区果冻| 久久亚洲国产成人精品v| 久久久精品94久久精品| 亚洲精华国产精华液的使用体验| 激情五月婷婷亚洲| 欧美变态另类bdsm刘玥| 男女高潮啪啪啪动态图| 1024香蕉在线观看| 成人漫画全彩无遮挡| 人妻 亚洲 视频| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品国产亚洲av天美| 免费播放大片免费观看视频在线观看| 日韩制服丝袜自拍偷拍| 色婷婷av一区二区三区视频| 香蕉精品网在线| 亚洲熟女精品中文字幕| 中文字幕亚洲精品专区| 精品国产乱码久久久久久小说| 久热久热在线精品观看| 亚洲天堂av无毛| 成人二区视频| 欧美bdsm另类| 亚洲精品一区蜜桃| 成人国语在线视频| 国产免费又黄又爽又色| 色哟哟·www| h视频一区二区三区| 亚洲国产欧美在线一区| 少妇精品久久久久久久| 精品久久蜜臀av无| 日韩视频在线欧美| 高清av免费在线| 国产精品香港三级国产av潘金莲 | 成人黄色视频免费在线看| 热re99久久精品国产66热6| 在线观看免费日韩欧美大片| 黄色 视频免费看| 如日韩欧美国产精品一区二区三区| 亚洲五月色婷婷综合| 亚洲精品日韩在线中文字幕| 人妻系列 视频| 超色免费av| 激情视频va一区二区三区| 成人手机av| 国产无遮挡羞羞视频在线观看| 亚洲av中文av极速乱| 亚洲精品成人av观看孕妇| 18在线观看网站| 一区二区三区激情视频| 观看av在线不卡| 三上悠亚av全集在线观看| 91aial.com中文字幕在线观看| 国产亚洲一区二区精品| 制服诱惑二区| 九九爱精品视频在线观看| 日韩制服骚丝袜av| 精品亚洲乱码少妇综合久久| 精品久久蜜臀av无| 两个人免费观看高清视频| 精品一区在线观看国产| 99re6热这里在线精品视频| 欧美日韩视频高清一区二区三区二| 亚洲精品美女久久久久99蜜臀 | 国产精品二区激情视频| 国产日韩一区二区三区精品不卡| 高清视频免费观看一区二区| 卡戴珊不雅视频在线播放| 国产欧美亚洲国产| 一级毛片我不卡| 黄频高清免费视频| 亚洲精品日本国产第一区| 中文天堂在线官网| 建设人人有责人人尽责人人享有的| 国产极品天堂在线| 久久久久久久大尺度免费视频| 视频区图区小说| 九九爱精品视频在线观看| 极品人妻少妇av视频| 老女人水多毛片| 午夜激情久久久久久久| 宅男免费午夜| 黄色毛片三级朝国网站| 成人亚洲欧美一区二区av| 我要看黄色一级片免费的| 天堂中文最新版在线下载| 久久久久精品性色| 在线观看免费日韩欧美大片| 亚洲国产欧美日韩在线播放| 久久这里只有精品19| 老女人水多毛片| 在线观看三级黄色| 日韩在线高清观看一区二区三区| 色视频在线一区二区三区| 伦理电影大哥的女人| 日本av手机在线免费观看| 亚洲久久久国产精品| 成人影院久久| 美女高潮到喷水免费观看| 丰满少妇做爰视频| 国产福利在线免费观看视频| 少妇人妻精品综合一区二区| 超色免费av| 啦啦啦啦在线视频资源| 久久久精品国产亚洲av高清涩受| av在线老鸭窝| 国产精品 欧美亚洲| 亚洲精品一二三| 我的亚洲天堂| 精品国产一区二区三区久久久樱花| 午夜日韩欧美国产| 寂寞人妻少妇视频99o| 国产片特级美女逼逼视频| 一级,二级,三级黄色视频| 亚洲五月色婷婷综合| 免费观看a级毛片全部| 亚洲精品一区蜜桃| 国产激情久久老熟女| 久久久欧美国产精品| 国产综合精华液| 90打野战视频偷拍视频| 欧美激情极品国产一区二区三区| 91精品国产国语对白视频| 亚洲人成77777在线视频| 在现免费观看毛片| 丝袜人妻中文字幕| 免费黄网站久久成人精品| 国产一区二区在线观看av| 国产色婷婷99| 成人18禁高潮啪啪吃奶动态图| 黄色配什么色好看| 免费日韩欧美在线观看| 大香蕉久久成人网| 伊人久久大香线蕉亚洲五| 精品久久蜜臀av无| 亚洲,欧美精品.| 久久 成人 亚洲| 黄频高清免费视频| 欧美精品亚洲一区二区| 日韩中字成人| 久久国产精品男人的天堂亚洲| 欧美精品av麻豆av| 美女高潮到喷水免费观看| 在线观看免费高清a一片| 亚洲中文av在线| 久久久久国产一级毛片高清牌| av又黄又爽大尺度在线免费看| 天天躁夜夜躁狠狠躁躁| 国产 一区精品| 国产黄色视频一区二区在线观看| 你懂的网址亚洲精品在线观看| 欧美日韩国产mv在线观看视频| 久久久久精品人妻al黑| 黄片播放在线免费| 99re6热这里在线精品视频| 亚洲国产精品一区三区| 1024香蕉在线观看| 欧美精品av麻豆av| 人妻一区二区av| 国产精品一区二区在线不卡| 日韩电影二区| 久久久精品94久久精品| 赤兔流量卡办理| 国产日韩一区二区三区精品不卡| 尾随美女入室| 在线观看国产h片| 亚洲第一青青草原| 在线观看美女被高潮喷水网站| 99re6热这里在线精品视频| 免费观看无遮挡的男女| 亚洲国产精品999| 日韩成人av中文字幕在线观看| 大香蕉久久网| 91精品伊人久久大香线蕉| 国产精品av久久久久免费| 美国免费a级毛片| 成人免费观看视频高清| 嫩草影院入口| 熟女电影av网| 中文字幕精品免费在线观看视频| 国产有黄有色有爽视频| 色婷婷久久久亚洲欧美| 一级,二级,三级黄色视频| 亚洲三区欧美一区| 青草久久国产| 蜜桃国产av成人99| 观看av在线不卡| 999久久久国产精品视频| 日韩精品有码人妻一区| 国产爽快片一区二区三区| 亚洲欧洲国产日韩| 黑人巨大精品欧美一区二区蜜桃| av福利片在线| 韩国精品一区二区三区| 国产一区二区在线观看av| 极品少妇高潮喷水抽搐| 日韩制服骚丝袜av| 高清在线视频一区二区三区| 久久这里只有精品19| 丝瓜视频免费看黄片| 亚洲国产av影院在线观看| 亚洲国产精品国产精品| 亚洲av成人精品一二三区| 三上悠亚av全集在线观看| 亚洲av成人精品一二三区| 久久精品国产亚洲av天美| 少妇人妻久久综合中文| 波多野结衣av一区二区av| 日本猛色少妇xxxxx猛交久久| 国产午夜精品一二区理论片| 欧美精品av麻豆av| 精品少妇黑人巨大在线播放| 日本av手机在线免费观看| 看免费成人av毛片| 人人妻人人添人人爽欧美一区卜| 久久av网站| 丰满迷人的少妇在线观看| 久久毛片免费看一区二区三区| 免费av中文字幕在线| 赤兔流量卡办理| 青青草视频在线视频观看| 精品少妇黑人巨大在线播放| 久久99热这里只频精品6学生| 黄色配什么色好看| 精品国产国语对白av| 人体艺术视频欧美日本| 老司机影院成人| 亚洲精品自拍成人| 一边摸一边做爽爽视频免费| 成人手机av| 久久久久久久大尺度免费视频| 看非洲黑人一级黄片| 精品亚洲成a人片在线观看| xxx大片免费视频| 国产色婷婷99| 久久ye,这里只有精品| 老汉色av国产亚洲站长工具| 高清不卡的av网站| 人体艺术视频欧美日本| 成年人午夜在线观看视频| 久久久亚洲精品成人影院| 麻豆精品久久久久久蜜桃| 波多野结衣一区麻豆| 在线观看人妻少妇| 看免费av毛片| 日韩人妻精品一区2区三区| 精品人妻一区二区三区麻豆| 我的亚洲天堂| 美女大奶头黄色视频| 十八禁高潮呻吟视频| 在线观看美女被高潮喷水网站| 久久99一区二区三区| 亚洲美女搞黄在线观看| 国产一区二区 视频在线| 亚洲第一区二区三区不卡| 纯流量卡能插随身wifi吗| 欧美精品av麻豆av| 欧美xxⅹ黑人| 欧美成人午夜免费资源| 两个人看的免费小视频| 永久免费av网站大全| 亚洲精品一区蜜桃| 色视频在线一区二区三区| 伊人亚洲综合成人网| 涩涩av久久男人的天堂| 久久精品久久精品一区二区三区| 亚洲激情五月婷婷啪啪| 波多野结衣av一区二区av| 中文字幕av电影在线播放| 国产精品国产三级专区第一集| 午夜影院在线不卡| 男男h啪啪无遮挡| 人人澡人人妻人| 在线天堂最新版资源| 老汉色∧v一级毛片| 欧美少妇被猛烈插入视频| 狠狠婷婷综合久久久久久88av| 纯流量卡能插随身wifi吗| 久久影院123| 人成视频在线观看免费观看| 老熟女久久久| 久久久久久久精品精品|