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

    GPT2w模型在南極地區(qū)精度分析

    2018-10-26 08:09:32姚宜斌單路路王澤民
    測繪學報 2018年10期
    關(guān)鍵詞:探空對流層水汽

    孔 建,姚宜斌,單路路,王澤民

    1. 武漢大學中國南極測繪研究中心,湖北 武漢 430079; 2. 武漢大學測繪學院,湖北 武漢 430079; 3. 武漢大學地球空間環(huán)境與大地測量教育部重點實驗室,湖北 武漢 430079

    GNSS定位過程中由于非色散效應影響,對流層延遲不能采用頻率組合的方式消除,因此其延遲量是影響定位精度的主要誤差源之一[1]。天頂方向?qū)α鲗友舆t影響一般能達到2~3 m左右,當高度角在10°以下時,傾斜對流層誤差影響可達20 m以上[2-3]。目前常采用模型改正的方法來削弱對流層延遲影響,即將信號傳播路徑上的對流層延遲通過映射函數(shù)投影到天頂方向,然后對天頂對流層延遲(zenith tropospheric delay,ZTD)進行建模[4-5]。

    傳統(tǒng)的對流層氣象參數(shù)模型如文獻[6—8]提出的模型改正精度可達厘米級,但模型計算要求輸入實測氣象參數(shù)。由于實測氣象參數(shù)難以獲得,又發(fā)展出一系列不依賴實測氣象參數(shù)的對流層經(jīng)驗模型。文獻[9—10]建立了UNB系列模型,UNB3m模型精度已達到4.9 cm;文獻[11]指出,歐盟星基廣域增強系統(tǒng)(The European Geostationary Navigation Overlay Service,EGNOS)采用EGNOS模型,平均精度約為5 cm;文獻[12]構(gòu)建了IGGtrop模型,雖然精度有所改進,但是參數(shù)數(shù)量太多;文獻[13]建立了TropGrid模型,以格網(wǎng)形式提供氣象元素的時序變化特征參數(shù);而后文獻[14]對該模型進行了簡化,提出了TropGrid2格網(wǎng)模型;文獻[15]建立了基于球諧函數(shù)表達的全球非氣象參數(shù)模型GZTD,計算簡便且精度較高。文獻[16]在文獻[17]提出的GPT2模型的基礎(chǔ)上,增加了大氣加權(quán)平均溫度Tm和水汽壓垂直梯度λ兩個參數(shù),建立了GPT2w模型,可同時獲得計算ZHD和ZWD所需的氣象元素,是目前精度最高、使用較為廣泛的對流層延遲模型[18]。GPT2w模型采用ECMWF(European Center for Medium-range Weather Forecast)提供的2001—2010年共10年的高質(zhì)量的ERA-37數(shù)據(jù)進行建模,氣象參數(shù)可由式(1)[18]計算得到

    (1)

    式中,r(t)為氣象參數(shù);DOY為年積日;A0為平均值;(A1,B1)為年周期振幅;(A2,B2)為半年周期振幅。

    式(1)得到的氣象參數(shù)是相對地面測站高度的,其他高度上是通過地面數(shù)值進行高度歸算得到的,如式(2)所示[19](http:∥ggosatm.hg.tuwien.ac.at/DELAY/SOURCE/GPT2w/)

    (2)

    式中,T0、P0、e0分別為格網(wǎng)點處溫度、氣壓、水汽壓;dT為溫度遞減率;dh為高度;T、P、e分別為格網(wǎng)點對應高度為dh處的溫度、氣壓、水汽壓;gm為重力,gm=9.806 65 m/s2;dMtr為干空氣的摩爾質(zhì)量;dMtr=28.965×10-3kg/mol;Rg為氣體常數(shù),Rg=8.314 3 J/K/mol;O為比濕;λ為水汽壓垂直梯度。

    對流層延遲由天頂靜力學延遲(zenith hydrostatic delay,ZHD)和天頂濕延遲(zenith wet delay,ZWD)組成,一般而言,ZHD可利用Saastamoinen模型[20]精確計算得出

    (3)

    式中,P、φ、H分別代表格網(wǎng)點處總氣壓、緯度和正高。

    而GPT2w模型計算ZWD可由加權(quán)平均溫度、水汽壓及其垂直梯度[16]通過式(4)得到

    (4)

    進而可得到對流層延遲

    ZTD=ZHD+ZWD

    (5)

    從式(1)可以看到,GPT2w模型氣象參數(shù)變化周期項僅考慮到年周期項和半年周期項,沒有考慮更短周期變化,如季節(jié)周期,日周期項。文獻[21]指出,GPT2w模型提供的氣壓、溫度、水汽壓在全球范圍內(nèi)均具有很高的精度和穩(wěn)定性,但在南北極地區(qū)精度仍然較低。文獻[16]利用全球341個IGS站對GPT2w模型進行了精度評估,發(fā)現(xiàn)其在全球范圍具有很好的精度,但是統(tǒng)計過程中南極區(qū)域IGS測站分布較少且不均勻。文獻[22—24]對比了目前的對流層延遲模型,指出GPT2w模型較其他模型精度高,與IGS站得到的對流層延遲最為接近。文獻[25]指出,GPT2w模型得到的氣壓等參數(shù)由于沒有考慮其日周期項,在寒冷季節(jié)存在較大誤差,可能更適用于熱帶地區(qū)。

    南極地區(qū)地理位置特殊,氣象條件不同于地球其他區(qū)域,目前專門研究GPT2w在南極區(qū)域精度評價的文章較少,而隨著后勤保障的巨大改善,在南極地區(qū)開展的大地測量科學研究日趨增多,因此有必要對GPT2w模型在南極地區(qū)的精度進行評估。另外,由于GPT2w模型沒有考慮日周期項,本文利用中國第33次南極科考實測探空數(shù)據(jù)進一步檢驗了模型在短時間尺度下的精度狀況。

    1 數(shù)據(jù)來源及數(shù)據(jù)處理方法

    無線電探空數(shù)據(jù)可以提供垂直氣象元素(如氣壓、溫度、水汽壓等)的剖面信息,且垂直分辨率和精度較高,經(jīng)常被作為參考(真值),用來評估其他方式獲取的氣象元素的精度[26-27]。本文選取兩種類型的探空數(shù)據(jù)來對GPT2w模型提供的氣壓、溫度、水汽壓進行精度評估:一是全球站點無線電探空資料數(shù)據(jù)集(the integrated global radiosonde archive,IGRA)中南極探空站數(shù)據(jù),二是中國第33次南極科考期間施放的探空氣球?qū)崪y數(shù)據(jù)。

    美國國家氣候數(shù)據(jù)中心(NCDC)研制的IGRA包含全球1538個探空站點自1938年至今的多層次、多要素逐日探空資料。該數(shù)據(jù)是目前全球范圍內(nèi)時空密度最高、資料最為完整的一套高質(zhì)量探空站點資料數(shù)據(jù)集[28]。探空站施放探空氣球的時間為世界時每天的00時和12時,但由于天氣、設備故障等原因,會導致數(shù)據(jù)缺失。

    IGRA探空站能直接提供氣壓、溫度、水汽壓等氣象參數(shù)。南極科考期間施放的探空氣球數(shù)據(jù)只能提供氣壓和溫度,不能直接獲得水汽壓,可通過露點溫度對水汽壓進行求解[29],具體公式如下

    (6)

    式中,a、b為常數(shù):溫度高于0℃時,a=17.26,b=237.29;溫度低于-40℃時,a=21.87,b=265.49;若溫度介于兩者之間,a、b用線性內(nèi)插獲得。

    為了對GPT2w模型在高度上的精度進行評估,本文按照氣球上升高度進行分層,以探究模型隨高度的精度變化。由于對流層及平流層沒有明顯的分層特性,結(jié)合數(shù)據(jù)的垂直分布特點采用不均勻分層[30],共分為9層,分層高度見表1,其中高度表示離地面高度。

    表1 高度分層信息

    采用偏差(Bias)和均方根誤差(RMS)作為精度評定標準,分別對分層數(shù)據(jù)進行評估。Bias和RMS分別由式(7)和式(8)計算得到

    (7)

    (8)

    式中,xmodel,i表示模型計算值;xreal,i表示探空數(shù)據(jù)氣象元素(本文視為真值);n為數(shù)據(jù)總個數(shù)。

    2 GPT2w模型南極地區(qū)精度檢驗及分析

    2.1 利用IGRA探空站數(shù)據(jù)進行模型氣象參數(shù)精度分析

    取南極地區(qū)探空站2016年全年數(shù)據(jù)作為參考值,共選取14個具有可用數(shù)據(jù)的探空站,與GPT2w模型計算的氣壓、溫度、水汽壓進行分層對比。探空站坐標信息見表2。由于南極特殊的地理位置,全年僅分為寒、暖兩季,4月—10月是寒季,11月—次年3月是暖季。取1月和7月分別作為暖季、寒季代表月份,對模型精度的季節(jié)性特征進行初步分析,并分層分季節(jié)統(tǒng)計Bias和RMS。

    表2 探空站坐標信息

    表3和表4為GPT2w模型氣壓、溫度、水汽壓的分層分季節(jié)Bias和RMS統(tǒng)計結(jié)果。從表中可以看出,隨著層數(shù)(高度)的增加,三者的Bias和RMS均呈現(xiàn)增大的趨勢;當?shù)竭_第7、8層左右,隨著高度增加,氣壓的Bias有減小趨勢,其Bias最大出現(xiàn)在第8層(10~15 km)左右,數(shù)值為7月的41 mb左右;溫度的Bias逐漸減小為負數(shù),但其絕對值仍遠大于地面附近數(shù)值,其Bias絕對值最大出現(xiàn)在第9層(15 km以上),數(shù)值為7月的-45 K左右;水汽壓的Bias和RMS也有增加趨勢,但是沒有氣壓和溫度變化顯著,其Bias最大出現(xiàn)在第9層(15 km以上),數(shù)值為7月的38 mb。對比發(fā)現(xiàn),在全年時間尺度上,GPT2w在地面附近精度較高,但是隨高度增加,精度大致呈下降趨勢。GPT2w模型一般提供地面高度的氣象參數(shù)初值,其他高度氣象參數(shù)的計算都是從地面進行歸算的,隨著高度增加,誤差會累積,從而導致精度降低。

    表3 1月GPT2w模型氣壓、溫度、水汽壓的精度統(tǒng)計

    表4 7月GPT2w模型氣壓、溫度、水汽壓的精度統(tǒng)計

    對比1月與7月統(tǒng)計結(jié)果,氣壓、溫度和水汽壓精度整體1月精度數(shù)值小于7月,說明1月精度較7月高,即GPT2w模型精度在南極區(qū)域呈現(xiàn)出季節(jié)性特征。

    2.2 利用南極探空氣球?qū)崪y數(shù)據(jù)進行模型氣象參數(shù)精度分析

    本文選取2016年11月到2017年4月第33次南極考察期間施放的環(huán)南極18顆探空氣球獲得的實測探空氣象數(shù)據(jù),通過高度分層對比其與GPT2w模型計算的氣溫、氣壓、水汽壓,對模型在南極地區(qū)的短時間尺度精度進行進一步分析。由于探空氣球是在航行過程中施放的,每個特定位置施放一個,達到既定高度氣球會自行爆炸,因此每顆氣球只能記錄當天的一組數(shù)據(jù),共有18組數(shù)據(jù),分布在不同時間和地點。表5為施放探空氣球的坐標信息。

    表5 探空氣球坐標信息

    2.2.1 利用ECMWF氣壓分層數(shù)據(jù)驗證氣球?qū)崪y數(shù)據(jù)

    歐洲中期天氣預報中心(European Centre for Medium-Range Weather Forecasts,ECMWF)是全球數(shù)值天氣預報組織,可提供高精度的ERA-Interim產(chǎn)品,在全球范圍內(nèi)包括南北極地區(qū)都具有很高的精度[20]。本文選取ECMWF氣壓分層數(shù)據(jù),與探空氣球?qū)崪y氣壓數(shù)據(jù)進行對比,來驗證氣球氣象數(shù)據(jù)的可靠性,以15185729號探空氣球數(shù)據(jù)為例,選取與其時間及地點對應的ECMWF氣壓分層數(shù)據(jù),兩者隨高度的變化如圖1所示??梢钥吹?,兩者符合較好,說明南極探空氣球?qū)崪y數(shù)據(jù)具有較高的可靠性,可以用來對GPT2w模型精度進行評估。

    2.2.2 利用探空氣球?qū)崪y數(shù)據(jù)對模型氣象參數(shù)進行精度分析

    表6為Bias和RMS統(tǒng)計結(jié)果。可以看出,對于氣壓和溫度而言,其趨勢與探空站統(tǒng)計結(jié)果類似,表現(xiàn)為地面精度最高,隨高度增加,精度下降。當?shù)竭_第7、8層左右時,氣壓的Bias有減小趨勢,Bias最大出現(xiàn)在第8層(10~15 km),數(shù)值為31 mb左右;溫度的Bias逐漸減小為負數(shù),但其絕對值仍遠大于地面附近數(shù)值,Bias絕對值最大出現(xiàn)在第9層(15 km以上),數(shù)值為-60 K左右;而水汽壓的Bias和RMS都比較小,在0附近波動,保持在基本穩(wěn)定狀態(tài),這是由于實測數(shù)據(jù)與GPT2w模型的水汽壓均是利用指數(shù)模型計算得來的,兩者之間差異較小。

    圖1 ECMWF氣壓分層數(shù)據(jù)與氣球?qū)崪y氣壓數(shù)據(jù)對比Fig.1 Comparison between ECMWF pressure data and measured pressure data from sounding balloons

    表6 利用實測數(shù)據(jù)對GPT2w模型氣壓、溫度、水汽壓的精度統(tǒng)計

    與2.1節(jié)相比,雖然GPT2w沒有考慮日周期項,但是溫度和氣壓的天內(nèi)精度與月統(tǒng)計精度在地面附近精度相當。但是垂直高度上,天內(nèi)精度與月統(tǒng)計精度相比有所下降,特別是與1月精度對比,隨著高度的增加精度在下降,可見未考慮日周期變化項對模型精度有所影響。

    圖2、圖3、圖4是GPT2w模型分別與探空站1月、7月及氣球?qū)崪y數(shù)據(jù)對比得到的氣壓、溫度、水汽壓的Bias,圖中給出了各層的Bias分布。對于氣壓和溫度,三者的趨勢有一致性,表現(xiàn)為地面附近精度較高,隨著高度增加精度下降;本文使用的實測數(shù)據(jù)采集于2016年11月至2017年4月,對應于2.1節(jié)統(tǒng)計的暖季區(qū)間,整體而言,GPT2w模型與實測數(shù)據(jù)的精度對比結(jié)果差于探空站1月對比精度,而優(yōu)于7月對比精度,與2.1節(jié)得到的暖季精度略優(yōu)于寒季精度的結(jié)果相符。而與探空站水汽壓增幅較大不同,氣球?qū)崪y數(shù)據(jù)的水汽壓較小且波動不大。

    圖2 GPT2w模型與探空站1月、7月及氣球?qū)崪y氣壓數(shù)據(jù)對比的BiasFig.2 Comparison of pressure between GPT2w and the radiosonde stations in January and July and sounding balloons,respectively

    圖3 GPT2w模型與探空站1月、7月及氣球?qū)崪y溫度數(shù)據(jù)對比的BiasFig.3 Comparison of temperature between GPT2w and the radiosonde stations in January and July and sounding balloons,respectively

    2.3 利用探空數(shù)據(jù)進行模型對流層延遲精度分析

    2.3.1 利用探空數(shù)據(jù)計算對流層延遲的方法及原理

    探空數(shù)據(jù)計算對流層延遲是根據(jù)其定義進行求解的。GPS衛(wèi)星信號在對流層中的延遲為[31]

    (9)

    圖4 GPT2w模型與探空站1月、7月及氣球?qū)崪y水汽壓數(shù)據(jù)對比的BiasFig.4 Comparison of water vapor pressure between GPT2w and the radiosonde stations in January and July and sounding balloons,respectively

    因此,ΔL可以進一步表示為

    (10)

    (11)

    式中,N=(n-1)·106表示大氣折射率。

    大氣折射率計算公式為[32]

    Nd+Nw

    (12)

    (13)

    1.75×10-4t2+1.44×10-6t3)

    (14)

    將對流層延遲繼續(xù)沿高度進行積分,則得到GPS信號穿過對流層天頂方向的總延遲

    (15)

    2.3.2 利用探空數(shù)據(jù)對模型對流層延遲進行精度分析

    利用上述計算方法,分別計算GPT2w模型與探空數(shù)據(jù)計算得到的對流層延遲,結(jié)果見表7、表8。

    表7 GPT2w模型與探空站計算結(jié)果對比

    數(shù)據(jù)對比時段:2017-01-16

    表8 GPT2w模型與實測探空數(shù)據(jù)對比

    由兩表可以看出,由GPT2w模型及探空數(shù)據(jù)計算得到的ZTD差距約為幾個厘米,這與其他學者驗證得到的GPT2w模型對流層延遲在全球范圍內(nèi)精度為厘米級的結(jié)論是一致的[16]。同時可以看出,ZWD數(shù)值較小,僅為幾個厘米,一般來說,在大氣中,ZHD占90%左右,ZWD占10%左右,而南極處于高緯度地區(qū),水汽含量很少,故ZWD相對來說更小。且兩者ZWD相差最多在氣球編號為15188344處,約為5.7 cm,數(shù)值較小。

    在探空數(shù)據(jù)計算結(jié)果作為真值的情況下,統(tǒng)計模型ZWD和ZTD的Bias及RMS,結(jié)果見表9??梢钥闯?,GPT2w得到的對流層延遲與探空數(shù)據(jù)計算得到的對流層延遲較為接近,Bias及RMS均處于厘米級,說明該模型可用于一般定位等對流層延遲的求解。但在精密定位中,仍需對該模型進行精化。

    表9 ZWD及ZTD的Bias、RMS統(tǒng)計

    3 結(jié)果討論及分析

    3.1 結(jié)果分析

    一般認為,探空氣球在上升過程中,只有高度發(fā)生變化,經(jīng)緯度不發(fā)生變化。

    (1) 由上述分析可知,在南極地區(qū),Bias和RMS隨著高度的增加而發(fā)生變化,GPT2w模型在地面附近精度最高;用ECMWF氣壓分層數(shù)據(jù)對實測氣壓數(shù)據(jù)進行檢驗,兩者符合較好,說明了實測數(shù)據(jù)具有可靠的測量精度,Bias和RMS隨高度變差主要是由于GPT2w模型本身誤差造成的。

    (2) GPT2w模型中,氣壓和溫度的計算均與高度有關(guān),氣壓采用高度的指數(shù)函數(shù)擬合,溫度采用高度的線性函數(shù)擬合得到,而水汽壓的計算與氣壓相關(guān)。因此,高度的變化,會導致模型精度的變化。GPT2w模型一般用來提供地面高度的氣象參數(shù)初值,其他高度的參數(shù)值模型內(nèi)部進行了函數(shù)改正,即高度歸算,而這種改正本身存在一定誤差。且會隨著高度增加,空氣成分及其他因素的變化,進一步導致模型精度下降,從而導致模型值與實測值出現(xiàn)較大不符。

    3.2 殘差擬合

    GPT2w模型計算的氣象參數(shù)值與實測數(shù)據(jù)在高度上存在較大的差異,本節(jié)對兩者之間的差異作進一步的分析。以15187632號探空氣球數(shù)據(jù)為例,圖5給出的是實測氣壓與GPT2w模型計算氣壓隨時間變化對比圖。從圖中可以看出,在地面附近,兩者相差不大,隨著高度的增加,兩者之間的差異在增大,在最高處略有下降。從兩者的殘差圖(圖6)也可以看出這一現(xiàn)象。圖6中開始時段異常是由于探空氣球是在走航期間施放的,開始時段氣球運動處于不穩(wěn)定的上下浮動狀態(tài)。對兩者殘差進行初步擬合,如圖7所示,可以看到,用三次多項式就可以較好的擬合兩者之間的差異。通過對溫度和水汽壓進行類似計算,可以得到相同的結(jié)論,如圖8、圖9所示。

    因此,下一步可以對GPT2w模型在高度上與實測數(shù)據(jù)的差異進行模型擬合改正,進而對該模型在南極地區(qū)乃至全球垂直方向上的精度進行改進。

    圖5 氣球?qū)崪y氣壓數(shù)據(jù)與模型計算氣壓對比Fig.5 Comparison of pressure between results calculated by GPT2w and results measured by sounding balloons

    圖6 氣球?qū)崪y氣壓數(shù)據(jù)與模型計算氣壓殘差Fig.6 Residual of pressure between results calculated by GPT2w and results measured by sounding balloons

    圖7 氣壓殘差三階多項式擬合 Fig.7 Three-order polynomial fitting of pressure’s residual

    圖8 溫度殘差三階多項式擬合Fig.8 Three-order polynomial fitting of temperature’s residual

    圖9 水汽壓殘差三階多項式擬合Fig.9 Three-order polynomial fitting of water vapor pressure’s residual

    4 結(jié) 論

    本文通過探空站和探空氣球提供的氣象數(shù)據(jù)對GPT2w模型在南極區(qū)域計算的氣壓、溫度、水汽壓進行精度檢驗,結(jié)果表明,該模型在地面附近精度可達到與全球其他區(qū)域精度一致的水平,但隨著高度的增加,Bias和RMS有不同程度的增加,精度有所下降。由模型計算得到的對流層延遲可達到厘米級精度,可用于一般定位等方面的對流層延遲改正。得到如下結(jié)論:

    (1) GPT2w一般用于地面氣象參數(shù)的計算,本文分析顯示,在地面附近,模型在南極地區(qū)精度與全球其他位置精度一致;且在南極區(qū)域,模型精度表現(xiàn)出季節(jié)性。

    (2) 通過與實測數(shù)據(jù)對比發(fā)現(xiàn),GPT2w模型雖然沒有考慮日周期項,但是對于南極區(qū)域而言,在地面高度上,計算的氣象參數(shù)天內(nèi)精度與月平均精度相當。但是隨著高度的增加,精度較1月精度有所下滑,而略優(yōu)于7月精度。

    (3) 隨著高度的增加,模型值與探空數(shù)據(jù)差異增加,除了氣球升空過程中產(chǎn)生漂移等原因外,主要原因是模型本身計算氣象元素時所使用的高度歸算函數(shù)存在一定誤差。同時可以發(fā)現(xiàn),殘差規(guī)律性較強,因此,可以在考慮季節(jié)性因素的基礎(chǔ)上,對其GPT2w模型加以改進,使其在垂直高度上精度得以改善。

    (4) 總體而言,由GPT2w模型計算的ZTD在南極地區(qū)有較高的精度,Bias和RMS均在厘米級,與全球其他地區(qū)得到的ZTD精度較為一致。

    猜你喜歡
    探空對流層水汽
    探空氣球為高空氣象“把脈”
    青藏高原上空平流層水汽的時空演變特征
    郴州地區(qū)對流層頂氣候概況
    實時干涉測量中對流層延遲與鐘差精修正建模
    載人航天(2016年4期)2016-12-01 06:56:24
    TK-2GPS人影火箭探空數(shù)據(jù)與L波段探空數(shù)據(jù)對比分析
    1979~2011年間平流層溫度及平流層水汽的演變趨勢
    成都地區(qū)2005~2015年對流層NO2柱濃度趨勢與時空分布
    深圳“5·11”特大暴雨過程的水汽輸送特征分析
    皖北地區(qū)對流層頂氣象特征分析
    淺談凈舉力對探空氣球升速及施放高度的影響
    河南科技(2014年14期)2014-02-27 14:12:08
    亚洲一区高清亚洲精品| 日本与韩国留学比较| 十八禁网站免费在线| 日本精品一区二区三区蜜桃| 久久韩国三级中文字幕| 国产国拍精品亚洲av在线观看| 有码 亚洲区| 男女之事视频高清在线观看| 毛片一级片免费看久久久久| 亚洲精品456在线播放app| av在线老鸭窝| 国产老妇女一区| 亚洲乱码一区二区免费版| 丰满人妻一区二区三区视频av| 久久久久久大精品| 在线免费十八禁| 午夜精品在线福利| 男女视频在线观看网站免费| 国产精品一区www在线观看| 精品99又大又爽又粗少妇毛片| 日韩一区二区视频免费看| 伦理电影大哥的女人| 91午夜精品亚洲一区二区三区| 国产精品,欧美在线| 一级毛片久久久久久久久女| 桃色一区二区三区在线观看| 免费av观看视频| 日本-黄色视频高清免费观看| 日韩av不卡免费在线播放| 国产私拍福利视频在线观看| 午夜福利在线观看吧| 又爽又黄无遮挡网站| 精品人妻一区二区三区麻豆 | 久久久久久九九精品二区国产| 国内精品美女久久久久久| 久久久久精品国产欧美久久久| 久久久久免费精品人妻一区二区| 男女边吃奶边做爰视频| 国产成人freesex在线 | 亚洲欧美日韩东京热| 国产在线精品亚洲第一网站| 99视频精品全部免费 在线| 亚洲av中文av极速乱| 国产成年人精品一区二区| 丝袜美腿在线中文| 内地一区二区视频在线| 黄色日韩在线| 精品一区二区免费观看| 欧美xxxx黑人xx丫x性爽| 成人二区视频| 人妻制服诱惑在线中文字幕| 精品久久久久久久久亚洲| 嫩草影院精品99| 99在线视频只有这里精品首页| 级片在线观看| 亚洲中文字幕日韩| 十八禁网站免费在线| 午夜福利视频1000在线观看| 久久久久久久亚洲中文字幕| 三级国产精品欧美在线观看| 九色成人免费人妻av| 久久九九热精品免费| 日韩欧美国产在线观看| 全区人妻精品视频| 国产亚洲91精品色在线| 一进一出抽搐动态| 91狼人影院| 亚洲精品亚洲一区二区| 我要看日韩黄色一级片| 国产精品伦人一区二区| 国产精品一区二区三区四区免费观看 | 色噜噜av男人的天堂激情| 欧美性感艳星| 又爽又黄无遮挡网站| 国产午夜精品论理片| av在线观看视频网站免费| 国产一区二区在线av高清观看| 老司机影院成人| 日本在线视频免费播放| 欧美区成人在线视频| 久久亚洲精品不卡| 国产高清有码在线观看视频| 精品久久久久久成人av| 一级黄色大片毛片| 精品午夜福利在线看| 亚洲精品国产av成人精品 | 国产成人福利小说| 热99re8久久精品国产| 成人一区二区视频在线观看| 国产三级中文精品| 美女xxoo啪啪120秒动态图| 午夜激情福利司机影院| 欧美人与善性xxx| ponron亚洲| 香蕉av资源在线| 一区二区三区免费毛片| 精品一区二区三区av网在线观看| 成人av一区二区三区在线看| 久久久久久大精品| 国产v大片淫在线免费观看| 国产免费男女视频| 男人舔女人下体高潮全视频| 亚洲精品国产av成人精品 | 听说在线观看完整版免费高清| 少妇人妻精品综合一区二区 | 大香蕉久久网| 久久这里只有精品中国| 亚洲av免费高清在线观看| 国产aⅴ精品一区二区三区波| 欧美日韩综合久久久久久| 99热网站在线观看| 亚洲国产精品成人久久小说 | 亚洲性久久影院| 国产乱人视频| 欧美性猛交╳xxx乱大交人| 午夜激情福利司机影院| 亚洲四区av| 最近视频中文字幕2019在线8| 国产精品亚洲美女久久久| 青春草视频在线免费观看| 又爽又黄a免费视频| 干丝袜人妻中文字幕| 秋霞在线观看毛片| 久久草成人影院| 女的被弄到高潮叫床怎么办| 高清毛片免费看| 老师上课跳d突然被开到最大视频| 亚洲人成网站在线播| 一级av片app| 久久久国产成人免费| 国产一区二区三区在线臀色熟女| 欧美日韩在线观看h| 国产午夜福利久久久久久| 不卡一级毛片| 高清午夜精品一区二区三区 | 波多野结衣高清作品| 99热网站在线观看| 看非洲黑人一级黄片| 热99在线观看视频| 国产亚洲精品久久久com| 国产色婷婷99| 91av网一区二区| av在线播放精品| 蜜桃亚洲精品一区二区三区| 欧美区成人在线视频| 国产精品国产三级国产av玫瑰| 午夜福利视频1000在线观看| 免费av观看视频| 成人三级黄色视频| 最后的刺客免费高清国语| 丝袜喷水一区| 最好的美女福利视频网| av卡一久久| 日韩欧美国产在线观看| 波多野结衣高清无吗| 俄罗斯特黄特色一大片| 亚洲图色成人| 午夜福利高清视频| 最近中文字幕高清免费大全6| 成年av动漫网址| 三级男女做爰猛烈吃奶摸视频| 黑人高潮一二区| 午夜精品在线福利| 国产又黄又爽又无遮挡在线| 精品一区二区三区人妻视频| 日本a在线网址| 直男gayav资源| 啦啦啦韩国在线观看视频| 亚洲高清免费不卡视频| 一本精品99久久精品77| 97热精品久久久久久| 性插视频无遮挡在线免费观看| 在现免费观看毛片| 一本久久中文字幕| 久久久久九九精品影院| 欧美色欧美亚洲另类二区| 天堂影院成人在线观看| 国产精品不卡视频一区二区| 国产极品精品免费视频能看的| 久久精品国产亚洲av香蕉五月| or卡值多少钱| 少妇裸体淫交视频免费看高清| 久久6这里有精品| 校园人妻丝袜中文字幕| 一个人观看的视频www高清免费观看| 国产精品嫩草影院av在线观看| 欧美另类亚洲清纯唯美| 日日摸夜夜添夜夜添小说| 中文亚洲av片在线观看爽| 午夜福利在线在线| 国产亚洲精品久久久com| 中国美白少妇内射xxxbb| 日本 av在线| 日本一二三区视频观看| 欧美绝顶高潮抽搐喷水| 男女下面进入的视频免费午夜| 99热网站在线观看| 亚洲精品乱码久久久v下载方式| 男人狂女人下面高潮的视频| 51国产日韩欧美| 免费观看的影片在线观看| 日韩 亚洲 欧美在线| 成人鲁丝片一二三区免费| 啦啦啦啦在线视频资源| 丝袜喷水一区| 插逼视频在线观看| 黄片wwwwww| 久久久久久大精品| 最后的刺客免费高清国语| 成人精品一区二区免费| 人人妻人人澡欧美一区二区| 久久综合国产亚洲精品| 国产探花极品一区二区| 精品人妻一区二区三区麻豆 | 久久久精品大字幕| 天堂网av新在线| 中文资源天堂在线| 搞女人的毛片| 国产精品久久久久久久电影| 色5月婷婷丁香| 亚洲欧美精品自产自拍| 亚洲成人中文字幕在线播放| 久久欧美精品欧美久久欧美| 亚洲欧美日韩无卡精品| 精品国产三级普通话版| 变态另类丝袜制服| 卡戴珊不雅视频在线播放| 三级毛片av免费| 男女那种视频在线观看| 国产亚洲精品av在线| 成年女人看的毛片在线观看| 日本与韩国留学比较| 亚洲自偷自拍三级| 色播亚洲综合网| 美女 人体艺术 gogo| 亚洲人成网站在线观看播放| 欧美区成人在线视频| 亚洲自拍偷在线| 国产精品久久久久久精品电影| 亚洲内射少妇av| 精品久久久久久久末码| 啦啦啦啦在线视频资源| 国产毛片a区久久久久| 国产不卡一卡二| 亚洲自偷自拍三级| 日韩欧美精品免费久久| 天堂影院成人在线观看| 亚洲中文字幕一区二区三区有码在线看| 12—13女人毛片做爰片一| 深夜精品福利| av视频在线观看入口| 亚洲av免费高清在线观看| 欧美成人免费av一区二区三区| 成人亚洲欧美一区二区av| av在线天堂中文字幕| 久久久精品大字幕| 美女黄网站色视频| 可以在线观看的亚洲视频| 亚洲国产日韩欧美精品在线观看| 久久久国产成人免费| 国产亚洲欧美98| 国产亚洲精品综合一区在线观看| 婷婷精品国产亚洲av在线| 国产精品日韩av在线免费观看| 国产一区二区激情短视频| 国产91av在线免费观看| 91av网一区二区| 最近2019中文字幕mv第一页| 成年女人看的毛片在线观看| 国产成人影院久久av| 成人高潮视频无遮挡免费网站| 变态另类丝袜制服| 国产 一区 欧美 日韩| av国产免费在线观看| 精品无人区乱码1区二区| 黄色日韩在线| 九色成人免费人妻av| 久久午夜亚洲精品久久| 嫩草影院入口| 真实男女啪啪啪动态图| 中文字幕av成人在线电影| 亚州av有码| 久久精品夜色国产| 国产伦一二天堂av在线观看| 日韩中字成人| 欧美成人a在线观看| 久久久久精品国产欧美久久久| 夜夜看夜夜爽夜夜摸| 五月玫瑰六月丁香| 国产欧美日韩一区二区精品| 九九在线视频观看精品| 中国美女看黄片| 伦精品一区二区三区| 国产伦精品一区二区三区视频9| 国产淫片久久久久久久久| 国产高清三级在线| 亚洲激情五月婷婷啪啪| av视频在线观看入口| 色视频www国产| 久久久久久久久久黄片| 少妇被粗大猛烈的视频| 美女 人体艺术 gogo| 精品一区二区三区人妻视频| 成人无遮挡网站| 欧美绝顶高潮抽搐喷水| 中文字幕熟女人妻在线| 日韩大尺度精品在线看网址| 久久精品国产清高在天天线| 老师上课跳d突然被开到最大视频| 男女下面进入的视频免费午夜| 三级毛片av免费| 国语自产精品视频在线第100页| 国国产精品蜜臀av免费| 九九热线精品视视频播放| 91av网一区二区| 男女之事视频高清在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲av不卡在线观看| 亚洲av五月六月丁香网| 国产午夜精品久久久久久一区二区三区 | 欧美日韩综合久久久久久| 亚洲av美国av| 淫妇啪啪啪对白视频| eeuss影院久久| 久久99热6这里只有精品| 亚洲精品亚洲一区二区| 99热全是精品| 国产三级在线视频| 性欧美人与动物交配| 亚洲精品日韩在线中文字幕 | 免费一级毛片在线播放高清视频| 国产一区二区在线av高清观看| 精品一区二区三区av网在线观看| 十八禁网站免费在线| 久久久a久久爽久久v久久| 国产不卡一卡二| 国产精品久久久久久久久免| 欧美一区二区精品小视频在线| 久久久久九九精品影院| 99视频精品全部免费 在线| 精品久久久久久久末码| 国产视频一区二区在线看| 深爱激情五月婷婷| 精品少妇黑人巨大在线播放 | 色哟哟·www| 中文字幕久久专区| 日韩欧美三级三区| 国产色婷婷99| 亚洲国产欧美人成| 国产激情偷乱视频一区二区| 淫秽高清视频在线观看| av在线亚洲专区| 露出奶头的视频| 赤兔流量卡办理| 波多野结衣高清作品| 精品久久久久久久久久免费视频| 国产欧美日韩精品亚洲av| 亚洲av不卡在线观看| 日韩中字成人| 久久久久久久久久久丰满| 久久婷婷人人爽人人干人人爱| 免费黄网站久久成人精品| 五月伊人婷婷丁香| 久久精品国产清高在天天线| 观看免费一级毛片| 天堂动漫精品| 国产激情偷乱视频一区二区| 舔av片在线| 中文资源天堂在线| 69人妻影院| 两性午夜刺激爽爽歪歪视频在线观看| av视频在线观看入口| 黄色配什么色好看| 久久6这里有精品| 99久久九九国产精品国产免费| 深爱激情五月婷婷| 黄色配什么色好看| 中国美女看黄片| 欧美区成人在线视频| 俺也久久电影网| 久久精品91蜜桃| 成人无遮挡网站| 国产亚洲精品久久久久久毛片| 久久午夜福利片| 九九在线视频观看精品| 好男人在线观看高清免费视频| 欧美不卡视频在线免费观看| 久久久久国内视频| 麻豆久久精品国产亚洲av| 久久久久国内视频| av在线蜜桃| 国产成人影院久久av| 成年女人看的毛片在线观看| 能在线免费观看的黄片| 男女啪啪激烈高潮av片| 国产69精品久久久久777片| 久久国产乱子免费精品| 亚洲图色成人| 五月玫瑰六月丁香| 少妇的逼水好多| 中国国产av一级| 特级一级黄色大片| 99久久无色码亚洲精品果冻| 你懂的网址亚洲精品在线观看 | 嫩草影院新地址| 国产精品国产高清国产av| 嫩草影院入口| 身体一侧抽搐| aaaaa片日本免费| 亚洲精品在线观看二区| 久久精品国产99精品国产亚洲性色| 人人妻人人澡欧美一区二区| 国产成人aa在线观看| 在线a可以看的网站| 夜夜爽天天搞| 极品教师在线视频| 欧美xxxx性猛交bbbb| 国产精品99久久久久久久久| 色哟哟·www| 又爽又黄a免费视频| 日韩中字成人| 特大巨黑吊av在线直播| 午夜福利在线在线| 狠狠狠狠99中文字幕| 如何舔出高潮| 男人舔奶头视频| 国产私拍福利视频在线观看| 51国产日韩欧美| 高清午夜精品一区二区三区 | 免费高清视频大片| 一本精品99久久精品77| 亚洲精品成人久久久久久| 一进一出抽搐动态| 欧美三级亚洲精品| 亚洲七黄色美女视频| 18+在线观看网站| 黑人高潮一二区| 亚洲在线观看片| 国产精品一区二区三区四区久久| 成人二区视频| 国产精品av视频在线免费观看| av在线老鸭窝| 狂野欧美激情性xxxx在线观看| 在线国产一区二区在线| 美女xxoo啪啪120秒动态图| 国产精品国产高清国产av| 国产成人影院久久av| 国产高清视频在线播放一区| 又爽又黄无遮挡网站| 男女做爰动态图高潮gif福利片| 日本一本二区三区精品| 国产一区亚洲一区在线观看| 美女被艹到高潮喷水动态| 欧美一区二区亚洲| 亚洲专区国产一区二区| 男插女下体视频免费在线播放| 亚洲中文字幕一区二区三区有码在线看| 亚洲人与动物交配视频| 亚洲精品在线观看二区| 免费看av在线观看网站| 久久婷婷人人爽人人干人人爱| 国产成人a区在线观看| 欧美一区二区国产精品久久精品| 成人一区二区视频在线观看| 十八禁网站免费在线| 国产精品三级大全| 久久精品国产亚洲av天美| avwww免费| 免费黄网站久久成人精品| 91午夜精品亚洲一区二区三区| 精品欧美国产一区二区三| 亚洲性久久影院| 国产男人的电影天堂91| 99久国产av精品| 国内精品美女久久久久久| 日韩高清综合在线| 国产午夜精品论理片| 五月玫瑰六月丁香| 尤物成人国产欧美一区二区三区| 国产一区二区三区av在线 | 国产蜜桃级精品一区二区三区| 免费看日本二区| 精华霜和精华液先用哪个| 老司机午夜福利在线观看视频| 欧美日本亚洲视频在线播放| 少妇熟女欧美另类| ponron亚洲| 亚洲欧美中文字幕日韩二区| 两个人的视频大全免费| 在线免费十八禁| 91午夜精品亚洲一区二区三区| 亚洲欧美中文字幕日韩二区| av福利片在线观看| 免费搜索国产男女视频| 欧美激情国产日韩精品一区| 舔av片在线| 亚洲成av人片在线播放无| 亚洲精品粉嫩美女一区| 国产av麻豆久久久久久久| 精品久久久久久久久av| 联通29元200g的流量卡| 五月玫瑰六月丁香| 精品久久久久久久久久免费视频| 国产成年人精品一区二区| 国产成人a区在线观看| 国产高清视频在线观看网站| 大又大粗又爽又黄少妇毛片口| 免费高清视频大片| 亚洲人成网站在线观看播放| 亚洲性夜色夜夜综合| 日韩国内少妇激情av| 色播亚洲综合网| 国产av麻豆久久久久久久| 日本黄色视频三级网站网址| 一边摸一边抽搐一进一小说| 精品不卡国产一区二区三区| 亚洲欧美清纯卡通| 免费看光身美女| 亚洲综合色惰| 中国国产av一级| 欧美日韩国产亚洲二区| 18禁在线播放成人免费| 亚洲国产日韩欧美精品在线观看| 狂野欧美激情性xxxx在线观看| 日韩一区二区视频免费看| or卡值多少钱| 亚洲无线观看免费| 两个人视频免费观看高清| 亚洲专区国产一区二区| 99久久精品国产国产毛片| 一区二区三区高清视频在线| 国产真实乱freesex| 大型黄色视频在线免费观看| 麻豆乱淫一区二区| 亚洲欧美日韩无卡精品| 成人二区视频| 91麻豆精品激情在线观看国产| 亚洲三级黄色毛片| 日韩强制内射视频| 国产伦一二天堂av在线观看| 最近中文字幕高清免费大全6| 婷婷亚洲欧美| 免费人成视频x8x8入口观看| 国产高清不卡午夜福利| 春色校园在线视频观看| 永久网站在线| 国产精品一区www在线观看| 久久久午夜欧美精品| 最新在线观看一区二区三区| 国产精品无大码| 亚洲久久久久久中文字幕| 香蕉av资源在线| 人妻丰满熟妇av一区二区三区| 97在线视频观看| 亚洲,欧美,日韩| 少妇猛男粗大的猛烈进出视频 | 别揉我奶头~嗯~啊~动态视频| 日韩人妻高清精品专区| 亚洲国产色片| 久久久午夜欧美精品| 成年女人毛片免费观看观看9| 久久亚洲精品不卡| 少妇裸体淫交视频免费看高清| 熟女人妻精品中文字幕| 欧美+日韩+精品| 亚洲最大成人av| 给我免费播放毛片高清在线观看| 直男gayav资源| 长腿黑丝高跟| 最近最新中文字幕大全电影3| 少妇的逼好多水| 午夜日韩欧美国产| 男女之事视频高清在线观看| 国语自产精品视频在线第100页| 午夜老司机福利剧场| .国产精品久久| 亚洲七黄色美女视频| 高清午夜精品一区二区三区 | 国产精品美女特级片免费视频播放器| 免费人成在线观看视频色| 免费看美女性在线毛片视频| 亚洲国产精品合色在线| а√天堂www在线а√下载| av天堂中文字幕网| 中国美女看黄片| 亚洲专区国产一区二区| 日韩亚洲欧美综合| 日日啪夜夜撸| 男人狂女人下面高潮的视频| 国产精品av视频在线免费观看| 国产色爽女视频免费观看| 国产国拍精品亚洲av在线观看| 一级毛片aaaaaa免费看小| 听说在线观看完整版免费高清| 在线观看66精品国产| 国产中年淑女户外野战色| 国产高清激情床上av| 国产乱人偷精品视频| 国产真实伦视频高清在线观看| 亚洲精品456在线播放app| 热99在线观看视频| 日韩大尺度精品在线看网址| 国产一区二区三区在线臀色熟女| 人妻丰满熟妇av一区二区三区| 欧美激情久久久久久爽电影| 一级毛片我不卡| 永久网站在线| 国产视频内射| 国产高清激情床上av| 99热精品在线国产| 美女黄网站色视频| 亚洲婷婷狠狠爱综合网| 国产女主播在线喷水免费视频网站 | 国产精品女同一区二区软件| 国产成人a区在线观看|