• <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
    成人亚洲精品一区在线观看| 精品99又大又爽又粗少妇毛片| 免费在线观看视频国产中文字幕亚洲 | 又大又黄又爽视频免费| 色网站视频免费| 久久久久国产网址| 亚洲,欧美,日韩| 一区福利在线观看| 久久久a久久爽久久v久久| 成人二区视频| 建设人人有责人人尽责人人享有的| 国产精品熟女久久久久浪| 久久99热这里只频精品6学生| 国产日韩一区二区三区精品不卡| 久久久久国产一级毛片高清牌| 男人操女人黄网站| 麻豆av在线久日| a级毛片在线看网站| 边亲边吃奶的免费视频| 人人妻人人澡人人爽人人夜夜| 老汉色av国产亚洲站长工具| 国产精品不卡视频一区二区| 久久久久久久亚洲中文字幕| 一个人免费看片子| 国产亚洲av片在线观看秒播厂| 97在线人人人人妻| 男人舔女人的私密视频| 久久99精品国语久久久| 免费在线观看完整版高清| 日韩在线高清观看一区二区三区| 日韩一区二区视频免费看| 在线观看一区二区三区激情| 男女午夜视频在线观看| 亚洲av在线观看美女高潮| 亚洲美女搞黄在线观看| 欧美 日韩 精品 国产| 建设人人有责人人尽责人人享有的| 麻豆乱淫一区二区| 精品国产超薄肉色丝袜足j| 一级黄片播放器| 一二三四在线观看免费中文在| 亚洲欧美日韩另类电影网站| 亚洲一区中文字幕在线| 亚洲四区av| 欧美 亚洲 国产 日韩一| 亚洲av免费高清在线观看| 在现免费观看毛片| 久久久久网色| 一级a爱视频在线免费观看| 久久人人爽人人片av| 少妇人妻精品综合一区二区| 人妻少妇偷人精品九色| 97人妻天天添夜夜摸| 97人妻天天添夜夜摸| 最近2019中文字幕mv第一页| 曰老女人黄片| 日韩 亚洲 欧美在线| 久久精品夜色国产| 秋霞在线观看毛片| 日韩一区二区三区影片| 三级国产精品片| 一级黄片播放器| 啦啦啦啦在线视频资源| 女人被躁到高潮嗷嗷叫费观| www.自偷自拍.com| 亚洲图色成人| 久久久国产欧美日韩av| 美女中出高潮动态图| 国产爽快片一区二区三区| 999久久久国产精品视频| 免费黄网站久久成人精品| 如何舔出高潮| 欧美 日韩 精品 国产| 国产一级毛片在线| 亚洲精品国产av蜜桃| 国产成人欧美| 999久久久国产精品视频| 伊人久久大香线蕉亚洲五| 三上悠亚av全集在线观看| 久久人人97超碰香蕉20202| 久久精品国产亚洲av天美| 久久久久久久大尺度免费视频| 伊人久久国产一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 国产淫语在线视频| 人妻一区二区av| 深夜精品福利| 亚洲中文av在线| 国产精品.久久久| 亚洲成av片中文字幕在线观看 | 在线天堂最新版资源| 欧美最新免费一区二区三区| 寂寞人妻少妇视频99o| 精品人妻偷拍中文字幕| 日韩av不卡免费在线播放| 在线观看免费日韩欧美大片| 中文欧美无线码| 日本av免费视频播放| 波多野结衣一区麻豆| 免费在线观看黄色视频的| 菩萨蛮人人尽说江南好唐韦庄| 日日撸夜夜添| 国产一区亚洲一区在线观看| 99热网站在线观看| 一本—道久久a久久精品蜜桃钙片| 一本—道久久a久久精品蜜桃钙片| 777米奇影视久久| 国产精品欧美亚洲77777| 激情五月婷婷亚洲| av国产久精品久网站免费入址| 美国免费a级毛片| 国产成人a∨麻豆精品| 久久国内精品自在自线图片| 最近2019中文字幕mv第一页| 狂野欧美激情性bbbbbb| 国产精品99久久99久久久不卡 | 午夜老司机福利剧场| 2021少妇久久久久久久久久久| 色婷婷久久久亚洲欧美| 少妇熟女欧美另类| 免费黄网站久久成人精品| 欧美在线黄色| 亚洲欧美成人精品一区二区| 九草在线视频观看| 国产免费现黄频在线看| 亚洲第一区二区三区不卡| 汤姆久久久久久久影院中文字幕| 飞空精品影院首页| 婷婷成人精品国产| 久久精品久久久久久久性| 自线自在国产av| 亚洲三级黄色毛片| 午夜免费男女啪啪视频观看| 欧美激情极品国产一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 纵有疾风起免费观看全集完整版| 1024视频免费在线观看| 少妇被粗大猛烈的视频| 99热网站在线观看| 国产精品欧美亚洲77777| 国产成人av激情在线播放| 欧美 亚洲 国产 日韩一| 国产精品一二三区在线看| www.自偷自拍.com| av视频免费观看在线观看| 九草在线视频观看| 亚洲av欧美aⅴ国产| av在线观看视频网站免费| av免费观看日本| 国产一区亚洲一区在线观看| 天堂8中文在线网| 人妻人人澡人人爽人人| 人妻系列 视频| 久久久久久人妻| 在线观看免费视频网站a站| 国产成人一区二区在线| 亚洲人成77777在线视频| 午夜日本视频在线| 精品一区二区免费观看| av又黄又爽大尺度在线免费看| 亚洲情色 制服丝袜| 又黄又粗又硬又大视频| 亚洲五月色婷婷综合| 在线观看人妻少妇| 国产又色又爽无遮挡免| 亚洲中文av在线| 国产成人精品久久二区二区91 | 在线天堂中文资源库| 国产有黄有色有爽视频| 丝袜人妻中文字幕| 一区二区av电影网| 亚洲av福利一区| 日本av手机在线免费观看| 一级毛片我不卡| 下体分泌物呈黄色| 国产成人欧美| 精品国产一区二区三区四区第35| 亚洲,欧美精品.| 叶爱在线成人免费视频播放| 母亲3免费完整高清在线观看 | 多毛熟女@视频| 久久这里只有精品19| 亚洲精品乱久久久久久| 一区在线观看完整版| 女人精品久久久久毛片| 久久人人爽av亚洲精品天堂| 久久国产亚洲av麻豆专区| 午夜激情av网站| 久久人人97超碰香蕉20202| 80岁老熟妇乱子伦牲交| 久久精品人人爽人人爽视色| 欧美精品一区二区大全| 一区在线观看完整版| 欧美精品国产亚洲| 丰满迷人的少妇在线观看| 亚洲 欧美一区二区三区| 婷婷色综合www| 亚洲国产欧美日韩在线播放| 一级爰片在线观看| 欧美激情极品国产一区二区三区| 黄频高清免费视频| 日本午夜av视频| 国产精品欧美亚洲77777| 王馨瑶露胸无遮挡在线观看| 丰满饥渴人妻一区二区三| 国产亚洲精品第一综合不卡| 高清视频免费观看一区二区| 成年美女黄网站色视频大全免费| 麻豆精品久久久久久蜜桃| 在线免费观看不下载黄p国产| 精品亚洲成a人片在线观看| av在线老鸭窝| 日本黄色日本黄色录像| 亚洲精品视频女| 国产精品成人在线| 亚洲av男天堂| xxxhd国产人妻xxx| 亚洲国产色片| 亚洲国产成人一精品久久久| 久久婷婷青草| 久久女婷五月综合色啪小说| 日韩av不卡免费在线播放| 波多野结衣一区麻豆| 亚洲av福利一区| 国精品久久久久久国模美| 少妇人妻精品综合一区二区| 国产日韩欧美在线精品| www.av在线官网国产| 男女边摸边吃奶| 久久这里只有精品19| 激情五月婷婷亚洲| av国产久精品久网站免费入址| 热re99久久精品国产66热6| 美女中出高潮动态图| 不卡av一区二区三区| a级毛片黄视频| 伊人久久国产一区二区| 精品少妇黑人巨大在线播放| 天堂中文最新版在线下载| 日本午夜av视频| 伦理电影免费视频| 亚洲婷婷狠狠爱综合网| 国产亚洲一区二区精品| 在线观看人妻少妇| 国产精品.久久久| 午夜av观看不卡| 亚洲综合色惰| 亚洲av电影在线观看一区二区三区| 午夜福利影视在线免费观看| 久久久精品区二区三区| 精品国产一区二区三区四区第35| 国产免费现黄频在线看| 成人手机av| 亚洲精品国产色婷婷电影| 一本大道久久a久久精品| 丁香六月天网| 捣出白浆h1v1| 久久99一区二区三区| 国产成人91sexporn| 日本av手机在线免费观看| 久久久久国产精品人妻一区二区| 国产深夜福利视频在线观看| 亚洲成av片中文字幕在线观看 | 桃花免费在线播放| 亚洲色图综合在线观看| 青春草国产在线视频| 久久这里有精品视频免费| 午夜福利在线观看免费完整高清在| 免费观看a级毛片全部| www.熟女人妻精品国产| 亚洲婷婷狠狠爱综合网| 久热久热在线精品观看| 国产午夜精品一二区理论片| 香蕉丝袜av| 一区二区av电影网| 国产精品女同一区二区软件| 欧美日韩精品网址| 亚洲欧美色中文字幕在线| 国产免费现黄频在线看| 日韩中字成人| 最近中文字幕高清免费大全6| 丰满饥渴人妻一区二区三| 亚洲国产看品久久| 欧美xxⅹ黑人| 又粗又硬又长又爽又黄的视频| 欧美人与性动交α欧美软件| 少妇精品久久久久久久| 亚洲成人av在线免费| 亚洲美女视频黄频| 亚洲婷婷狠狠爱综合网| www.自偷自拍.com| 国产免费现黄频在线看| 欧美bdsm另类| av国产精品久久久久影院| 亚洲成人一二三区av| 中文字幕人妻丝袜一区二区 | 男人添女人高潮全过程视频| 日韩电影二区| 免费观看在线日韩| 中国三级夫妇交换| 夜夜骑夜夜射夜夜干| 在线观看人妻少妇| 一级毛片黄色毛片免费观看视频| 少妇的丰满在线观看| 2021少妇久久久久久久久久久| 波多野结衣av一区二区av| 女人久久www免费人成看片| 日本色播在线视频| 卡戴珊不雅视频在线播放| 欧美成人精品欧美一级黄| 久久久久精品人妻al黑| 欧美日韩精品成人综合77777| 高清不卡的av网站| 国产欧美日韩一区二区三区在线| 亚洲成av片中文字幕在线观看 | 一区二区三区四区激情视频| 美女视频免费永久观看网站| 看十八女毛片水多多多| 亚洲欧洲国产日韩| 色播在线永久视频| 精品人妻一区二区三区麻豆| 亚洲av日韩在线播放| 免费av中文字幕在线| 亚洲五月色婷婷综合| 亚洲,欧美精品.| 五月伊人婷婷丁香| 综合色丁香网| 国产爽快片一区二区三区| 午夜福利在线免费观看网站| 久久久久精品久久久久真实原创| 久久午夜福利片| 亚洲少妇的诱惑av| 日本-黄色视频高清免费观看| 一边亲一边摸免费视频| 少妇 在线观看| 亚洲中文av在线| 日韩大片免费观看网站| 在线 av 中文字幕| 国产av一区二区精品久久| 美女大奶头黄色视频| 中文字幕最新亚洲高清| 国产精品蜜桃在线观看| 国产在线免费精品| 免费高清在线观看视频在线观看| 国产 一区精品| 永久网站在线| 国产亚洲av片在线观看秒播厂| 日韩大片免费观看网站| 飞空精品影院首页| 观看美女的网站| 九色亚洲精品在线播放| 久久婷婷青草| 麻豆av在线久日| 免费观看性生交大片5| av在线老鸭窝| 涩涩av久久男人的天堂| 99国产综合亚洲精品| 国产亚洲欧美精品永久| 韩国av在线不卡| 国产欧美亚洲国产| 少妇人妻久久综合中文| 欧美人与性动交α欧美软件| 国产精品香港三级国产av潘金莲 | 狂野欧美激情性bbbbbb| 考比视频在线观看| 青草久久国产| 咕卡用的链子| 精品一区二区三卡| 欧美日韩成人在线一区二区| 国产1区2区3区精品| 亚洲国产毛片av蜜桃av| 国产精品 国内视频| 欧美精品av麻豆av| 久久久久久久国产电影| 国产有黄有色有爽视频| 国产精品 欧美亚洲| 蜜桃国产av成人99| 看十八女毛片水多多多| 色婷婷av一区二区三区视频| 中文乱码字字幕精品一区二区三区| 欧美人与善性xxx| 亚洲国产av影院在线观看| 高清欧美精品videossex| 午夜免费观看性视频| 在线观看免费日韩欧美大片| 亚洲成人手机| 日本欧美视频一区| 国产 一区精品| 丰满少妇做爰视频| 精品一品国产午夜福利视频| 制服诱惑二区| 精品久久久精品久久久| 日韩一区二区视频免费看| 美女主播在线视频| 欧美人与性动交α欧美软件| 1024香蕉在线观看| 极品人妻少妇av视频| 久久久久网色| 久久久久精品人妻al黑| 99久久中文字幕三级久久日本| 欧美成人午夜精品| 国产男女超爽视频在线观看| 国产精品欧美亚洲77777| 亚洲精品国产av蜜桃| 国产成人精品久久二区二区91 | 欧美精品一区二区大全| 免费观看性生交大片5| 日本欧美视频一区| 成人手机av| 丰满少妇做爰视频| 国产又色又爽无遮挡免| 精品一区二区免费观看| 亚洲精品成人av观看孕妇| 亚洲中文av在线| 国产男人的电影天堂91| 黄色 视频免费看| 国产成人精品无人区| 久久热在线av| 涩涩av久久男人的天堂| 午夜日本视频在线| 性色avwww在线观看| 两个人免费观看高清视频| 最新中文字幕久久久久| 热99国产精品久久久久久7| 伦精品一区二区三区| 国产成人91sexporn| 一级黄片播放器| 国产成人精品在线电影| 欧美中文综合在线视频| 伊人久久大香线蕉亚洲五| 纵有疾风起免费观看全集完整版| 国产1区2区3区精品| 两个人免费观看高清视频| 免费播放大片免费观看视频在线观看| 夜夜骑夜夜射夜夜干| 亚洲,欧美精品.| 在线观看三级黄色| 久久久久久人人人人人| 老鸭窝网址在线观看| 亚洲婷婷狠狠爱综合网| 考比视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 97精品久久久久久久久久精品| 国产男女超爽视频在线观看| 激情五月婷婷亚洲| 18禁观看日本| 欧美精品人与动牲交sv欧美| 高清不卡的av网站| 国产又爽黄色视频| 99热国产这里只有精品6| 热99国产精品久久久久久7| 大香蕉久久网| 日韩av免费高清视频| www.自偷自拍.com| 久久99精品国语久久久| 午夜福利影视在线免费观看| 1024香蕉在线观看| 久久久久国产一级毛片高清牌| 老汉色∧v一级毛片| 国产成人精品无人区| av在线播放精品| 国产不卡av网站在线观看| 国产97色在线日韩免费| 国产97色在线日韩免费| 91精品伊人久久大香线蕉| 欧美+日韩+精品| 丰满乱子伦码专区| 97在线人人人人妻| 18在线观看网站| av不卡在线播放| 国产免费现黄频在线看| 国产成人av激情在线播放| 在线观看一区二区三区激情| 99香蕉大伊视频| 黄色视频在线播放观看不卡| 久久久久久久久久久免费av| 中文字幕最新亚洲高清| 久久久国产欧美日韩av| 看非洲黑人一级黄片| 国产精品免费视频内射| 中文字幕人妻丝袜制服| 91精品国产国语对白视频| 午夜福利,免费看| 国产亚洲欧美精品永久| 亚洲av在线观看美女高潮| 99香蕉大伊视频| 亚洲精品第二区| 久久99一区二区三区| 午夜福利影视在线免费观看| 国产一区二区激情短视频 | av网站在线播放免费| 国产深夜福利视频在线观看| 欧美日韩精品成人综合77777| 免费在线观看视频国产中文字幕亚洲 | 高清不卡的av网站| 宅男免费午夜| 大码成人一级视频| 伊人久久国产一区二区| 哪个播放器可以免费观看大片| 一级黄片播放器| 久久99一区二区三区| 欧美 亚洲 国产 日韩一| 黄色一级大片看看| 伦理电影免费视频| 一级爰片在线观看| 少妇人妻久久综合中文| 晚上一个人看的免费电影| 国产极品天堂在线| 亚洲国产精品成人久久小说| 熟女电影av网| 日日摸夜夜添夜夜爱| 国产探花极品一区二区| 高清av免费在线| 巨乳人妻的诱惑在线观看| 久久久精品国产亚洲av高清涩受| 纯流量卡能插随身wifi吗| 五月伊人婷婷丁香| 欧美成人午夜精品| 亚洲一区二区三区欧美精品| 五月天丁香电影| 你懂的网址亚洲精品在线观看| 大片免费播放器 马上看| 精品少妇黑人巨大在线播放| 国产亚洲最大av| 成人午夜精彩视频在线观看| 亚洲综合色网址| 久久久国产欧美日韩av| 少妇人妻久久综合中文| 男女午夜视频在线观看| 日本wwww免费看| 乱人伦中国视频| 哪个播放器可以免费观看大片| 亚洲精品在线美女| 日产精品乱码卡一卡2卡三| 亚洲国产精品成人久久小说| 亚洲av日韩在线播放| av国产久精品久网站免费入址| 女性生殖器流出的白浆| 国产成人av激情在线播放| 亚洲男人天堂网一区| 成人18禁高潮啪啪吃奶动态图| 蜜桃在线观看..| 性少妇av在线| 高清在线视频一区二区三区| 五月伊人婷婷丁香| 国产精品一二三区在线看| 欧美日韩成人在线一区二区| 韩国av在线不卡| 这个男人来自地球电影免费观看 | 我要看黄色一级片免费的| 99国产综合亚洲精品| 丰满少妇做爰视频| 国产熟女欧美一区二区| av电影中文网址| 你懂的网址亚洲精品在线观看| 丝袜美足系列| 美女高潮到喷水免费观看| 日韩成人av中文字幕在线观看| 人妻一区二区av| 我的亚洲天堂| 国产探花极品一区二区| 青青草视频在线视频观看| 黄片小视频在线播放| 国产成人av激情在线播放| 伦精品一区二区三区| 亚洲情色 制服丝袜| 深夜精品福利| 观看美女的网站| 精品人妻熟女毛片av久久网站| 天美传媒精品一区二区| 久久鲁丝午夜福利片| 校园人妻丝袜中文字幕| 老汉色av国产亚洲站长工具| 啦啦啦在线免费观看视频4| av有码第一页| 亚洲,一卡二卡三卡| 美女高潮到喷水免费观看| 国产精品久久久久久av不卡| 国产高清不卡午夜福利| 美国免费a级毛片| 免费观看av网站的网址| 99久国产av精品国产电影| 亚洲av欧美aⅴ国产| 国产精品熟女久久久久浪| 亚洲国产精品一区二区三区在线| 777久久人妻少妇嫩草av网站| 人人妻人人澡人人爽人人夜夜| 久久综合国产亚洲精品| 亚洲第一青青草原| 中国国产av一级| 国产色婷婷99| 少妇人妻久久综合中文| 丝袜喷水一区| 亚洲综合色网址| 少妇被粗大的猛进出69影院| 两个人看的免费小视频| 成年人免费黄色播放视频| 色视频在线一区二区三区| 女人久久www免费人成看片| 国产精品久久久久久av不卡| 久久国产精品大桥未久av| 免费在线观看视频国产中文字幕亚洲 | 久久精品国产鲁丝片午夜精品| 亚洲伊人久久精品综合| av在线老鸭窝| 亚洲内射少妇av| 少妇的逼水好多| 人人澡人人妻人| 久久久久国产网址| 日韩欧美精品免费久久| 亚洲第一区二区三区不卡| 曰老女人黄片| 国产精品成人在线| 男人舔女人的私密视频| 一本—道久久a久久精品蜜桃钙片|