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

    基于地基GPS的武漢地區(qū)電離層TEC建模和變化分析

    2016-09-24 08:03:33孫偉喬煒周凱魏致富姚鵬
    城市勘測 2016年4期
    關(guān)鍵詞:區(qū)域模型

    孫偉,喬煒,周凱,魏致富,姚鵬

    (武漢市測繪研究院,湖北武漢 430022)

    基于地基GPS的武漢地區(qū)電離層TEC建模和變化分析

    孫偉*,喬煒,周凱,魏致富,姚鵬

    (武漢市測繪研究院,湖北武漢 430022)

    利用武漢地區(qū)的GPS觀測數(shù)據(jù),建立和比較了多項式函數(shù)模型、廣義三角級數(shù)、低階球諧函數(shù)模型等三種常用區(qū)域電離層模型在武漢地區(qū)電離層TEC建模的適用性,結(jié)果顯示三者擬合精度相當(dāng)。在此基礎(chǔ)上,選用多項式函數(shù)模型建立了武漢地區(qū)近一個太陽周期的電離層TEC模型,對武漢地區(qū)電離層TEC的年度變化、半年度異常、周日變化等進行了分析,得出了一些有益的結(jié)論。

    地基GPS;武漢地區(qū);電離層TEC

    1 引 言

    相比電離層測高儀、非相干散射雷達等傳統(tǒng)的電離層探測方法,利用GPS研究電離層具有探測范圍更廣、產(chǎn)品更豐富、更易形成空間一體化監(jiān)測網(wǎng)絡(luò)等優(yōu)勢,特別是利用地基GPS觀測值建立電離層模型的方法已成為電離層空間物理研究的主要手段[1]。常用的區(qū)域電離層函數(shù)擬合模型有多項式模型、三角級數(shù)函數(shù)模型和低階球諧函數(shù)模型等[1,2]。柳景斌[3]基于中國區(qū)域IGS站GPS觀測數(shù)據(jù),比較分析了三種區(qū)域電離層TEC模型的擬合效果,認(rèn)為在一定條件下三種模型可以認(rèn)為是等價的;楊凱等[4]利用中國地殼運動監(jiān)測網(wǎng)絡(luò)的GPS觀測數(shù)據(jù),建立和評估了三種電離層模型,認(rèn)為所建的中國區(qū)域上空TEC模型符合較好且差別較小。目前,我國很多大型城市均建有較為密集的、常年運行的GPS基準(zhǔn)站,極大地豐富了區(qū)域上空的GPS觀測資料,為建立區(qū)域范圍長時間序列的電離層模型,分析電離層的精細(xì)變化提供了有利條件。

    為了分析、比較各種模型在武漢區(qū)域范圍的適用性,本文利用武漢地區(qū)2005年~2014年的GPS實測數(shù)據(jù),實現(xiàn)了上述三種函數(shù)方法的電離層TEC建模,在此基礎(chǔ)上對武漢地區(qū)電離層TEC的季節(jié)變化、周日變化、半年異常等進行了分析。

    2 基于地基GPS的區(qū)域電離層建模方法

    為簡化計算,在實際研究和計算中通常采用單層電離層模型,即假定所有的自由電子集中在一個無限薄層的地球表面上空,用這個薄層替代整個電離層,用H表示該單層高度。為使H更準(zhǔn)確地描述電離層特征,一般設(shè)定為電離層最大電子密度層所在的高度,約為300 km~400 km[2],本文所述三種函數(shù)的單層模型H取值均為350 km。

    2.1多項式函數(shù)模型

    多項式擬合模型是一種結(jié)構(gòu)簡單,顧及電離層隨緯度、太陽時角變化特征的一種模型,能夠在小區(qū)域范圍內(nèi)取得較好的擬合效果。該模型是以穿刺點與擬合位置地理中心的緯度差、太陽時角作為參數(shù),按階對電離層TEC進行擬合,其模型公式如下[1,2]:

    式中,Eij表示模型擬合系數(shù),φ、φ0分別表示穿刺點處和擬合區(qū)域中心的地理緯度,S、S0分別表示穿刺點處和擬合區(qū)域中心的太陽時角,n、m表示階數(shù),本文取值為4、5。

    2.2廣義三角級數(shù)模型

    多項式擬合模型的局限性在于需要利用數(shù)個小時的數(shù)據(jù)才能達到一個理想的精度水平。建立一種基于三角級函數(shù)的電離層擬合方法(Trigonometric Series Function,TSF)可以有效改善區(qū)域建模效果[4],但是模型中部分參數(shù)的固定卻限制了區(qū)域電離層特征變化的擬合能力。有學(xué)者在此基礎(chǔ)上對模型進行了擴展,提出在地磁坐標(biāo)系下建立參數(shù)可變的廣義三角級數(shù)模型(GTSF),其模型公式如下[4]:

    式中,Ai表示函數(shù)模型系數(shù)(本文建立模型時計算系數(shù)為15個),φm表示穿刺點的地磁緯度,h表示時間相關(guān)變量。h=2π(t-14)/T,T=14h,t表示穿刺點的地方時;φm=φ+0.064cos(λ-1.617),φ和λ分別表示地理緯度和地理經(jīng)度。

    2.3低階球諧函數(shù)模型

    局部區(qū)域的函數(shù)模型一般不能用來描述全球電離層TEC特征,國外學(xué)者上世紀(jì)90年代提出用球諧函數(shù)(SH)來建立全球電離層TEC[5]。球諧函數(shù)是利用球諧展開級數(shù)函數(shù)來描述單層電離層模型的,可以用來建立全球區(qū)域的電離層模型,其數(shù)學(xué)表達式為[6]:

    這里需要指出的是,球諧函數(shù)的本意是應(yīng)用于全球范圍映射,當(dāng)進行區(qū)域范圍研究時雖然能得到正確的擬合結(jié)果,但是其模型的參數(shù)已經(jīng)不能夠用來描述區(qū)域范圍的電離層特性,不再具有特定的物理意義[7]。

    3 武漢地區(qū)電離層TEC建模與分析

    3.1電離層TEC建模

    本文建模數(shù)據(jù)來自武漢市連續(xù)運行衛(wèi)星定位服務(wù)系統(tǒng)(WHCORS)6個基準(zhǔn)站的GPS觀測數(shù)據(jù),6個基準(zhǔn)站均勻分布于武漢市域范圍,建模區(qū)域為112°~115°E、29°~31°N。電離層TEC建模時,采用IGS事后精密星歷計算衛(wèi)星位置;采用基于卡爾曼濾波的載波相位平滑P碼偽距方法提高偽距觀測值精度。采用每2 h的觀測數(shù)據(jù)與單日固定值的接收機硬件延遲和衛(wèi)星硬件延遲,共同建立法方程,再用最小二乘求解;解算出每2 h的模型參數(shù)和單日硬件延遲后,便可以得出測站上空高精度的電離層TEC。

    電離層TEC的幅值變化和異常發(fā)生主要受太陽輻射、地磁活動以及其他天氣現(xiàn)象的各種擾動的影響,其中太陽輻射是最主要的因素。這里給出了2005年~2014年的太陽活動指數(shù),包括太陽黑子數(shù)和F10.7指數(shù),如圖1所示,從圖可以看出2007年~2009年為太陽活動低年,2010開始太陽黑子數(shù)逐漸增多,2011年~2014年為太陽活動高年,太陽輻射量較大。由于電離層函數(shù)模型的擬合精度隨著太陽輻射的增強而降低[8],所以選擇在太陽活動高年比較和分析電離層模型在武漢地區(qū)的適用性更為適當(dāng)。

    本文選取2013年WHCORS基準(zhǔn)站的GPS實測數(shù)據(jù),分別采用多項式擬合、廣義三角級數(shù)和低階球諧函數(shù)模型進行計算,模型參數(shù)個數(shù)按上述建模方法設(shè)置。求得2013年一年時間的模型系數(shù)(每日12組)之后利用模型系數(shù)反算回各個穿刺點處的TEC值,統(tǒng)計得到各個模型的每日殘差均值和殘差標(biāo)準(zhǔn)差均值,每日殘差均值接近于0。表1給出的是殘差均值和殘差標(biāo)準(zhǔn)差的年均值,從標(biāo)準(zhǔn)差精度來看三個模型在武漢地區(qū)范圍的擬合精度水平相當(dāng),廣義三角級數(shù)模型稍差。究其原因,本文統(tǒng)一采用每2 h的觀測數(shù)據(jù)進行建模分析,而廣義三角級數(shù)模型由于其模型定義的物理意義,一般更適用于用一天的數(shù)據(jù)進行單日擬合[9]。

    圖1 2005年~2014年太陽黑子數(shù)和F10.7 cm指數(shù)

    三種常用單層電離層模型的擬合精度比較(單位/TECU) 表1

    為了進一步驗證利用模型擬合電離層TEC的準(zhǔn)確性,這里選取武漢市地理位置中心(30°32′,114° 20′)上空的擬合結(jié)果與CODE GIM發(fā)布值進行比較。三種模型是利用解算的模型參數(shù),通過輸入坐標(biāo)和時間獲取該位置的2013年的電離層TEC日均值,而CODE GIM值是采用IGS推薦的四點格網(wǎng)內(nèi)插方法,內(nèi)插的該位置2013年每日12個時刻的TEC然后取平均所得,如圖2所示。圖中綠色實線代表CODE GIM,可以看出在武漢地區(qū)范圍內(nèi)三種模型的擬合效果差異不大,比CODE GIM值相比全年整體變化趨勢保持一致,相關(guān)系數(shù)分別達到了0.990 5、0.973 4、0.988 7,可見幾種模型的擬合結(jié)果是可靠、穩(wěn)定、無偏的。由于CODE是全球擬合結(jié)果,采用的地基GPS數(shù)據(jù)在中國區(qū)域范圍內(nèi)僅僅使用了不到10個跟蹤站的數(shù)據(jù),相比其他幾種模型其在局部小范圍的擬合精度并不高,由此導(dǎo)致與實測數(shù)據(jù)的差異。

    圖2 三種電離層函數(shù)模型電離層每日均值比較

    3.2電離層TEC變化分析

    根據(jù)上述分析結(jié)論,在區(qū)域范圍三種電離層函數(shù)模型擬合精度相當(dāng)?shù)幕A(chǔ)上,本節(jié)選用多項式函數(shù)模型,利用2005年~2014年WHCORS基準(zhǔn)站的GPS實測數(shù)據(jù),建立了高精度、高時空分辨率的武漢地區(qū)電離層TEC模型。下面利用電離層TEC實測數(shù)據(jù)對武漢地區(qū)電離層TEC季節(jié)變化和周日變化進行了分析,如圖3所示。

    圖3 2005年~2014年武漢地區(qū)電離層TEC每日均值變化

    由于繞太陽公轉(zhuǎn)和黃赤交角的存在,太陽直射和斜射對電離層具有不同程度的影響,從而使得電離層呈現(xiàn)出年周期性和季節(jié)性變化[10]。圖3繪制了武漢地區(qū)近一個太陽周期年的電離層TEC每日均值(縱坐標(biāo)單位為TECU、橫坐標(biāo)單位為日)。電離層TEC值域表現(xiàn)和太陽黑子數(shù)量大小高度相關(guān),太陽活動低年(2005年~2010年)電離層TEC值明顯小于太陽活動高年(2011年~2014年),且全年電離層TEC變化較為平緩。以2008年為例,全年電離層TEC每日均值在0 TECU~10 TECU之間緩慢變化,可見影響電離層TEC值大小的主要因素是太陽輻射強度。整體上來看,不管是太陽活動低年還是高年都呈現(xiàn)出周年變化特性:每年4月份~5月份為電離層活躍期,太陽活動高年2012年電離層TEC峰值超過50 TECU;每年1月份、12月份為電離層平靜期,太陽活動低年2008年這2個月的電離層TEC每日均值僅為1 TECU~3 TECU。

    圖4 2013年武漢地區(qū)電離層TEC每日均值變化

    為了進一步分析武漢地區(qū)電離層TEC的周日變化規(guī)律,圖4給出了武漢地區(qū)2013年20 m時間分辨率的電離層TEC值(縱坐標(biāo)單位為小時、橫坐標(biāo)單位為日)。圖中全年TEC時序變化與圖3中2013年日均值變化相互吻合,2013年4月~5月電離層TEC達到全年峰值。從每日時間分布來看,2013年2月~8月期間白天TEC峰值出現(xiàn)時間為當(dāng)?shù)貢r間14:00~15:00(UTC06:00~07:00),而在秋分過后,白天TEC峰值時間變?yōu)楫?dāng)?shù)貢r間13:00~14:00(UTC05:00~06:00)。全年夜間電離層TEC谷值3月~8月出現(xiàn)時間為當(dāng)?shù)貢r間凌晨04:00 (UTC20:00)左右,1月~2月、9月~12月的當(dāng)?shù)貢r間20:00~次日08:00(UTC12:00~24:00)電離層TEC均接近于0,上述變化分析符合中低緯地區(qū)“正午過后達到峰值”的電離層周日變化規(guī)律[10]。結(jié)合圖3和圖4,可以看到武漢地區(qū)全年電離層TEC值變化存在半年規(guī)律,即每年1月~2月TEC較低,4月~5月達到一個波峰值,然后在7月~8月又顯著降低,來到11月左右時又達到一個波峰,隨著冬季來臨又逐漸降低,而且前一個波峰值總是大于后一個波峰值。這種半年的周期變化規(guī)律在太陽活動高年比太陽活動低年更為顯著,在電離層物理研究中,這是典型的中低緯地區(qū)電離層半年異?,F(xiàn)象,即全年電離層TEC峰值出現(xiàn)在春秋分時期,而非夏季,在北半球通??梢杂^測到該現(xiàn)象。有學(xué)者從該現(xiàn)象產(chǎn)生的物理機制出發(fā),認(rèn)為可能受到電場相關(guān)的電動力學(xué)過程控制,導(dǎo)致中低緯和赤道地區(qū)電離層呈現(xiàn)出相對于周年變化幅度更大的半年變化幅度及其明顯的地磁對稱特征[11,12]。

    4 結(jié) 語

    (1)在武漢地區(qū)范圍內(nèi)利用常用的三種電離層模型擬合電離層TEC時,擬合精度相當(dāng),其中廣義三角級數(shù)模型擬合效果稍差。相比CODE GIM值相比全年整體變化趨勢保持一致,可見三種模型的擬合結(jié)果是可靠、穩(wěn)定、無偏的。由于CODE是全球擬合結(jié)果,相比其他幾種模型其在局部小范圍的擬合精度并不高。

    (2)通過與太陽指數(shù)數(shù)據(jù)的對比可見,武漢地區(qū)電離層變化和太陽輻射密切相關(guān)。在不同的時間周期內(nèi):武漢地區(qū)電離層TEC變化表現(xiàn)出顯著的半年異常,即全年峰值出現(xiàn)在兩分季節(jié);周日變化屬于典型的中低緯地區(qū)電離層變化規(guī)律,即每日14:00左右電離層TEC達到日峰值。

    [1]李征航,黃勁松.GPS測量與數(shù)據(jù)處理[M].武漢:武漢大學(xué)出版社,2005.

    [2]李征航,張小紅.衛(wèi)星導(dǎo)航定位新技術(shù)及高精度數(shù)據(jù)處理方法[M].武漢:武漢大學(xué)出版社,2009.

    [3]柳景斌.基于地基GPS的區(qū)域電離層TEC球冠諧分析及預(yù)報[D].武漢:武漢大學(xué),2008.

    [4]楊凱,薛軍琛,宋淑麗等.基于CMONOC建立和評估中國區(qū)域電離層模型[J].大地測量與地球動力學(xué),2016,36(2):138~142.

    [5]Georgiadiou,Y.1994.Modeling the ionosphere for an active control network of GPS station.LGR-Series,(7):1.

    [6]Schaer,S.,W.Gurtner,J.Feltens.1997.IONEX:The IONosphere Map Exchange Format Version 1,Proceedings of the 1998 IGS Analysis Centers Workshop,ESOC,Darmstadt,Germany,9~11 February 1998.

    [7]Schaer,S.1999.Mapping and Predicting the Earth′s Ionosphere Using the Global Positioning System.Bern,Switzerland:University of Bern.

    [8]章紅平.基于地基GPS的中國區(qū)域電離層監(jiān)測與延遲改正研究[D].上海:中國科學(xué)院上海天文臺.

    [9]陳鵬,陳家君.全球電離層格網(wǎng)模型精度影響因素分析[J].大地測量與地球動力學(xué),2014,34(2):105~110.

    [10]袁運斌.基于GPS的電離層監(jiān)測及延遲改正理論與方法的研究[D].武漢:中國科學(xué)院測量與地球物理研究所,2002.

    [11]熊年祿,唐存琛,李行健.電離層物理概論[M].武漢:武漢大學(xué)出版社,1999.

    [12]余濤,萬衛(wèi)星,劉立波等.利用IGS數(shù)據(jù)分析全球TEC的周年和半年變化特性[J].地球物理學(xué)報,2006,49(4):943~949.

    [13]霍星亮,袁運斌,歐吉坤.基于GPS資料研究中國區(qū)域電離層TEC的周日變化、半年度及冬季異?,F(xiàn)象[J].自然科學(xué)進展,2005,15(5):626~630.

    Modeling and Analysis of Ionospheric TEC in Wuhan Area Based on the Ground-based GPS

    Sun Wei,Qiao Wei,Zhou Kai,Wei Zhifu,Yao Peng
    (Wuhan Geomatics Institute,Wuhan 430022,China)

    Using GPS observation data in Wuhan area,the polynomial function model,generalized trigonometric series and low order spherical harmonic function model of three commonly used regional ionospheric model applicability in Wuhan Ionospheric TEC modeling are established and compared.The results showe that 3 fitting accuracy is based on GPS data in Wuhan.On this basis,the paper choose polynomial function model was established in Wuhan area near a solar cycle of the Ionospheric TEC model,annual change of TEC in Wuhan area,abnormal semi-annual and diurnal changes were analyzed,and some useful conclusions are drawn.

    Ground-based GPS;wuhan area;ionospheric TEC

    1672-8262(2016)04-114-05

    P228.4

    A

    2016—03—17

    孫偉(1983—),男,博士,高級工程師,主要從事GNSS應(yīng)用研究、城市基礎(chǔ)控制測量等工作。

    國家自然科學(xué)基金青年項目(41504023)

    猜你喜歡
    區(qū)域模型
    一半模型
    永久基本農(nóng)田集中區(qū)域“禁廢”
    分割區(qū)域
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    3D打印中的模型分割與打包
    關(guān)于四色猜想
    分區(qū)域
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    基于嚴(yán)重區(qū)域的多PCC點暫降頻次估計
    電測與儀表(2015年5期)2015-04-09 11:30:52
    ponron亚洲| 桃色一区二区三区在线观看| 亚洲熟妇中文字幕五十中出| 国产欧美日韩一区二区三| 人人妻,人人澡人人爽秒播| 国产成+人综合+亚洲专区| 赤兔流量卡办理| 在线看三级毛片| 亚洲黑人精品在线| 99久久无色码亚洲精品果冻| 一区二区三区免费毛片| 69av精品久久久久久| 国产精品日韩av在线免费观看| av欧美777| 婷婷六月久久综合丁香| 久久精品国产清高在天天线| 99在线人妻在线中文字幕| 一个人免费在线观看的高清视频| 午夜福利视频1000在线观看| 男女那种视频在线观看| 99久久久亚洲精品蜜臀av| 我的女老师完整版在线观看| 90打野战视频偷拍视频| 国产亚洲av嫩草精品影院| 国内久久婷婷六月综合欲色啪| 免费看光身美女| 自拍偷自拍亚洲精品老妇| 成人特级av手机在线观看| 国产欧美日韩一区二区精品| 看免费av毛片| 国内精品美女久久久久久| 99久久精品热视频| 夜夜爽天天搞| 黄色丝袜av网址大全| 少妇高潮的动态图| 久久精品国产清高在天天线| 男人和女人高潮做爰伦理| 国产在视频线在精品| 国产黄片美女视频| 91麻豆精品激情在线观看国产| 国产伦人伦偷精品视频| 亚洲欧美精品综合久久99| 色5月婷婷丁香| 人人妻,人人澡人人爽秒播| 淫秽高清视频在线观看| x7x7x7水蜜桃| 在线观看舔阴道视频| 好男人电影高清在线观看| 国产成人福利小说| 亚洲国产色片| 亚洲中文日韩欧美视频| 国产男靠女视频免费网站| a级毛片免费高清观看在线播放| 欧美成狂野欧美在线观看| 国产 一区 欧美 日韩| 免费av不卡在线播放| 五月伊人婷婷丁香| 中国美女看黄片| 亚洲第一电影网av| 99视频精品全部免费 在线| 色综合亚洲欧美另类图片| 夜夜爽天天搞| 午夜精品一区二区三区免费看| 色哟哟哟哟哟哟| 欧美黑人欧美精品刺激| 久9热在线精品视频| 午夜精品一区二区三区免费看| 身体一侧抽搐| 国产精品日韩av在线免费观看| 亚洲人成电影免费在线| 免费av观看视频| 91午夜精品亚洲一区二区三区 | 亚洲国产精品成人综合色| 亚洲一区高清亚洲精品| 蜜桃久久精品国产亚洲av| 欧美又色又爽又黄视频| 午夜福利高清视频| 久久亚洲精品不卡| 最近视频中文字幕2019在线8| 黄色视频,在线免费观看| 中国美女看黄片| aaaaa片日本免费| 国产aⅴ精品一区二区三区波| 波多野结衣巨乳人妻| 日日摸夜夜添夜夜添小说| 中文字幕av在线有码专区| 波野结衣二区三区在线| 亚洲最大成人中文| 国产精品亚洲一级av第二区| 成人国产一区最新在线观看| 精品国产三级普通话版| 在线免费观看的www视频| 极品教师在线免费播放| 亚洲精华国产精华精| 三级毛片av免费| 特大巨黑吊av在线直播| 中文字幕久久专区| 最好的美女福利视频网| 老鸭窝网址在线观看| 亚洲精品乱码久久久v下载方式| 一级黄色大片毛片| 成年女人永久免费观看视频| 性欧美人与动物交配| 亚洲最大成人中文| 一二三四社区在线视频社区8| 欧美日韩福利视频一区二区| 免费av观看视频| 少妇人妻一区二区三区视频| 高清毛片免费观看视频网站| 赤兔流量卡办理| 此物有八面人人有两片| 99久久无色码亚洲精品果冻| 欧美潮喷喷水| 一级a爱片免费观看的视频| 色播亚洲综合网| 一本精品99久久精品77| 神马国产精品三级电影在线观看| 啪啪无遮挡十八禁网站| 国产在视频线在精品| 色哟哟哟哟哟哟| 淫秽高清视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 一区二区三区激情视频| 99国产综合亚洲精品| 国产精品伦人一区二区| 久久香蕉精品热| 变态另类丝袜制服| 色av中文字幕| 国产免费一级a男人的天堂| 欧美日韩中文字幕国产精品一区二区三区| 色综合欧美亚洲国产小说| av天堂中文字幕网| 床上黄色一级片| 免费在线观看成人毛片| 小说图片视频综合网站| 国产一区二区在线av高清观看| 国产日本99.免费观看| 黄色配什么色好看| 天堂√8在线中文| 精品无人区乱码1区二区| 国产免费男女视频| 亚洲自偷自拍三级| 51国产日韩欧美| 特级一级黄色大片| 欧美精品国产亚洲| 色综合欧美亚洲国产小说| 国产亚洲精品久久久久久毛片| 变态另类丝袜制服| 国产成人啪精品午夜网站| 亚洲中文字幕一区二区三区有码在线看| 一二三四社区在线视频社区8| 一二三四社区在线视频社区8| 国产熟女xx| 国产黄片美女视频| 成人三级黄色视频| 婷婷六月久久综合丁香| 国产伦精品一区二区三区视频9| 色在线成人网| 他把我摸到了高潮在线观看| 久久精品影院6| 久久久久久久久久黄片| 美女xxoo啪啪120秒动态图 | 校园春色视频在线观看| av专区在线播放| 九九在线视频观看精品| 少妇丰满av| 搡老妇女老女人老熟妇| 日本撒尿小便嘘嘘汇集6| 免费一级毛片在线播放高清视频| 小蜜桃在线观看免费完整版高清| 精品久久久久久久久亚洲 | 少妇高潮的动态图| 午夜福利在线观看免费完整高清在 | 草草在线视频免费看| 成年人黄色毛片网站| 日日摸夜夜添夜夜添av毛片 | 老鸭窝网址在线观看| 啦啦啦观看免费观看视频高清| 最好的美女福利视频网| 麻豆国产av国片精品| 老司机深夜福利视频在线观看| 亚洲av.av天堂| 少妇人妻精品综合一区二区 | 欧美日韩乱码在线| 亚洲人成网站高清观看| 午夜福利18| 俺也久久电影网| 午夜日韩欧美国产| www.999成人在线观看| 麻豆av噜噜一区二区三区| 女同久久另类99精品国产91| 嫩草影院新地址| 淫妇啪啪啪对白视频| 亚洲美女黄片视频| 男插女下体视频免费在线播放| 看黄色毛片网站| 日本 av在线| 国产成人aa在线观看| 国内精品一区二区在线观看| 欧美性猛交黑人性爽| 亚洲国产精品久久男人天堂| 日韩成人在线观看一区二区三区| 亚洲国产色片| 亚洲成人免费电影在线观看| 最近中文字幕高清免费大全6 | 国产伦在线观看视频一区| 色综合亚洲欧美另类图片| 少妇人妻精品综合一区二区 | 免费搜索国产男女视频| 91在线观看av| 日本精品一区二区三区蜜桃| a在线观看视频网站| 欧美乱妇无乱码| 免费看a级黄色片| 欧美精品国产亚洲| 十八禁人妻一区二区| 高清在线国产一区| 国产精品久久久久久久电影| 精品久久久久久久久av| 18禁黄网站禁片午夜丰满| 真人做人爱边吃奶动态| 一区二区三区高清视频在线| 国产av不卡久久| 国产精华一区二区三区| 啦啦啦观看免费观看视频高清| 黄色日韩在线| 久久这里只有精品中国| aaaaa片日本免费| 国产亚洲精品久久久久久毛片| 美女xxoo啪啪120秒动态图 | 国产欧美日韩一区二区三| 亚洲成a人片在线一区二区| 亚洲国产高清在线一区二区三| 欧美日韩乱码在线| 18+在线观看网站| 免费av毛片视频| 97超级碰碰碰精品色视频在线观看| 欧美日韩综合久久久久久 | 欧美xxxx黑人xx丫x性爽| 在线国产一区二区在线| 欧美精品国产亚洲| 天堂影院成人在线观看| 午夜福利欧美成人| 午夜激情福利司机影院| 免费在线观看影片大全网站| 美女cb高潮喷水在线观看| 国产伦在线观看视频一区| 亚洲av成人av| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 麻豆国产av国片精品| 欧美不卡视频在线免费观看| 观看美女的网站| 看免费av毛片| 亚洲自偷自拍三级| 国产久久久一区二区三区| 日本五十路高清| 真人一进一出gif抽搐免费| 国产探花极品一区二区| 俄罗斯特黄特色一大片| 麻豆一二三区av精品| 日韩欧美一区二区三区在线观看| 老司机福利观看| 日韩欧美国产一区二区入口| av女优亚洲男人天堂| 三级男女做爰猛烈吃奶摸视频| 乱人视频在线观看| 国产色婷婷99| 国产高清视频在线观看网站| 真实男女啪啪啪动态图| 久久99热这里只有精品18| 又爽又黄无遮挡网站| 日韩欧美 国产精品| 97热精品久久久久久| 国产亚洲精品综合一区在线观看| 国产成人福利小说| 免费观看精品视频网站| 99热只有精品国产| 精品久久久久久久久久久久久| 亚洲成人精品中文字幕电影| 99热这里只有是精品在线观看 | 欧美丝袜亚洲另类 | 精品久久久久久久久av| 精品一区二区三区av网在线观看| 久久中文看片网| 欧美日韩亚洲国产一区二区在线观看| 搡老岳熟女国产| x7x7x7水蜜桃| 日韩精品青青久久久久久| 国产精品嫩草影院av在线观看 | 国产乱人视频| 麻豆国产av国片精品| 在线a可以看的网站| 69人妻影院| 中文字幕av在线有码专区| 亚洲熟妇中文字幕五十中出| 亚洲av成人精品一区久久| 最新中文字幕久久久久| 如何舔出高潮| 18禁黄网站禁片免费观看直播| 热99re8久久精品国产| 成年女人看的毛片在线观看| 俺也久久电影网| 又黄又爽又免费观看的视频| 黄色视频,在线免费观看| 日日摸夜夜添夜夜添小说| 97超视频在线观看视频| 日韩欧美国产在线观看| 国产欧美日韩一区二区三| 欧美日韩亚洲国产一区二区在线观看| eeuss影院久久| xxxwww97欧美| av天堂在线播放| 别揉我奶头~嗯~啊~动态视频| 欧美日韩综合久久久久久 | 久久久国产成人免费| 99视频精品全部免费 在线| 亚洲黑人精品在线| 国产黄片美女视频| 国产精品久久视频播放| 三级男女做爰猛烈吃奶摸视频| 欧美zozozo另类| 欧美日韩国产亚洲二区| 欧美性感艳星| 韩国av一区二区三区四区| 国产三级黄色录像| 丁香六月欧美| 女生性感内裤真人,穿戴方法视频| 真实男女啪啪啪动态图| 亚洲自偷自拍三级| 国产色婷婷99| 精品久久久久久久久av| 国产伦精品一区二区三区四那| 97超视频在线观看视频| 成人美女网站在线观看视频| 男女床上黄色一级片免费看| 免费观看精品视频网站| 色哟哟·www| 国内精品久久久久精免费| 免费电影在线观看免费观看| 观看免费一级毛片| 在线免费观看的www视频| 精品久久久久久久久久久久久| 精品国产三级普通话版| 亚洲第一电影网av| 丰满人妻一区二区三区视频av| 久久伊人香网站| 一a级毛片在线观看| 婷婷丁香在线五月| 国内揄拍国产精品人妻在线| 蜜桃久久精品国产亚洲av| 色噜噜av男人的天堂激情| 男人狂女人下面高潮的视频| 国产一区二区在线观看日韩| 99久国产av精品| 悠悠久久av| 亚洲国产欧洲综合997久久,| 99久久精品热视频| 国模一区二区三区四区视频| 国产欧美日韩精品一区二区| 午夜影院日韩av| 最后的刺客免费高清国语| 中文字幕高清在线视频| 亚洲国产欧洲综合997久久,| 色尼玛亚洲综合影院| 亚洲美女黄片视频| 午夜激情欧美在线| 亚洲成人久久性| 色播亚洲综合网| 国产aⅴ精品一区二区三区波| 国产三级在线视频| 老鸭窝网址在线观看| 搡老熟女国产l中国老女人| 久久精品国产亚洲av香蕉五月| 两人在一起打扑克的视频| 在线观看66精品国产| 精品久久国产蜜桃| 老司机午夜福利在线观看视频| 黄色日韩在线| 性色av乱码一区二区三区2| 搞女人的毛片| 夜夜爽天天搞| 夜夜看夜夜爽夜夜摸| 国产亚洲av嫩草精品影院| 永久网站在线| bbb黄色大片| 中文字幕av成人在线电影| 无人区码免费观看不卡| av天堂在线播放| 最新在线观看一区二区三区| 3wmmmm亚洲av在线观看| 欧美国产日韩亚洲一区| 亚洲,欧美精品.| 国内精品久久久久精免费| 国产av不卡久久| 99久久精品一区二区三区| 亚洲精品粉嫩美女一区| 18+在线观看网站| 99国产综合亚洲精品| 老司机深夜福利视频在线观看| 狠狠狠狠99中文字幕| 国内精品一区二区在线观看| 18美女黄网站色大片免费观看| 一级黄色大片毛片| 1024手机看黄色片| 国产乱人视频| 又爽又黄无遮挡网站| 精品久久久久久久久亚洲 | 免费人成在线观看视频色| 精品久久国产蜜桃| 欧美潮喷喷水| 国产精品女同一区二区软件 | 久久精品久久久久久噜噜老黄 | av福利片在线观看| 欧美成狂野欧美在线观看| 老熟妇仑乱视频hdxx| bbb黄色大片| 人人妻人人澡欧美一区二区| 中文字幕免费在线视频6| 亚洲片人在线观看| 国产v大片淫在线免费观看| 一区二区三区高清视频在线| xxxwww97欧美| 成人三级黄色视频| 国产亚洲欧美98| 人妻夜夜爽99麻豆av| 色播亚洲综合网| 搞女人的毛片| 国产精品野战在线观看| 一级a爱片免费观看的视频| 国产成人av教育| 国产老妇女一区| 我要搜黄色片| 国产亚洲精品av在线| 国产高清视频在线观看网站| 亚洲中文日韩欧美视频| 狠狠狠狠99中文字幕| 三级男女做爰猛烈吃奶摸视频| 给我免费播放毛片高清在线观看| 亚洲av成人不卡在线观看播放网| 婷婷丁香在线五月| 99热只有精品国产| 免费av观看视频| 欧美成人性av电影在线观看| av天堂中文字幕网| 成人性生交大片免费视频hd| 无遮挡黄片免费观看| 亚洲av电影不卡..在线观看| 久久久久久久久久黄片| 99久久99久久久精品蜜桃| 国产高清有码在线观看视频| 日韩中文字幕欧美一区二区| 一进一出抽搐动态| 欧美在线黄色| 亚洲av成人精品一区久久| 九九热线精品视视频播放| 欧美+日韩+精品| 99久久无色码亚洲精品果冻| 麻豆久久精品国产亚洲av| 午夜激情福利司机影院| 男女下面进入的视频免费午夜| 九色国产91popny在线| 亚洲熟妇熟女久久| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲美女视频黄频| 午夜久久久久精精品| www.999成人在线观看| netflix在线观看网站| 欧美激情在线99| 亚洲国产高清在线一区二区三| 一进一出抽搐动态| 国产一区二区在线观看日韩| 精品一区二区三区av网在线观看| av在线天堂中文字幕| 黄色一级大片看看| 好男人在线观看高清免费视频| 日本五十路高清| 亚洲美女搞黄在线观看 | 亚洲专区中文字幕在线| 国产精品伦人一区二区| 在线免费观看的www视频| 小蜜桃在线观看免费完整版高清| 69av精品久久久久久| 永久网站在线| 久久久国产成人精品二区| 久久婷婷人人爽人人干人人爱| 看免费av毛片| 91在线精品国自产拍蜜月| .国产精品久久| 99在线视频只有这里精品首页| 午夜影院日韩av| av黄色大香蕉| 欧美一区二区国产精品久久精品| 黄色一级大片看看| 别揉我奶头~嗯~啊~动态视频| 亚洲国产精品sss在线观看| 亚洲天堂国产精品一区在线| 一个人免费在线观看的高清视频| 一个人看的www免费观看视频| 欧美乱妇无乱码| 国产野战对白在线观看| 久久久久精品国产欧美久久久| 成年女人永久免费观看视频| 伦理电影大哥的女人| 国产欧美日韩一区二区精品| 一级a爱片免费观看的视频| av在线老鸭窝| 亚洲av一区综合| 欧美日韩综合久久久久久 | 日韩欧美在线二视频| 国产av不卡久久| 日本熟妇午夜| 亚洲成av人片免费观看| 亚洲乱码一区二区免费版| 午夜久久久久精精品| 国产一区二区三区在线臀色熟女| 三级男女做爰猛烈吃奶摸视频| 亚洲人与动物交配视频| 草草在线视频免费看| 国产激情偷乱视频一区二区| 91久久精品国产一区二区成人| 国产高潮美女av| 成年女人毛片免费观看观看9| 国产91精品成人一区二区三区| 日韩欧美在线乱码| 色精品久久人妻99蜜桃| 国产精品久久久久久亚洲av鲁大| 日本在线视频免费播放| a级毛片a级免费在线| 精品福利观看| 久99久视频精品免费| 国产一区二区在线观看日韩| 可以在线观看毛片的网站| 日韩av在线大香蕉| 亚洲 国产 在线| 日本在线视频免费播放| 嫩草影院精品99| 熟女电影av网| 久久久国产成人精品二区| 88av欧美| 757午夜福利合集在线观看| 麻豆成人午夜福利视频| 一a级毛片在线观看| 99久久精品一区二区三区| 老鸭窝网址在线观看| 久久精品人妻少妇| www.熟女人妻精品国产| 99国产综合亚洲精品| 日本黄色视频三级网站网址| 成人鲁丝片一二三区免费| 精品午夜福利在线看| 人妻制服诱惑在线中文字幕| 国产欧美日韩一区二区精品| 男女下面进入的视频免费午夜| 真人做人爱边吃奶动态| 一边摸一边抽搐一进一小说| 精品人妻1区二区| 国产在视频线在精品| 偷拍熟女少妇极品色| 亚洲中文字幕一区二区三区有码在线看| 成年女人毛片免费观看观看9| 国产欧美日韩精品亚洲av| 国内精品一区二区在线观看| 国产精品美女特级片免费视频播放器| 一区福利在线观看| 国产在线男女| 在线观看午夜福利视频| 亚洲乱码一区二区免费版| 亚洲欧美日韩高清在线视频| 亚洲av熟女| 91午夜精品亚洲一区二区三区 | 永久网站在线| 色av中文字幕| 午夜福利免费观看在线| 国产精品久久久久久久电影| 99久国产av精品| 国产亚洲欧美在线一区二区| 女同久久另类99精品国产91| 国产单亲对白刺激| 国产一区二区在线av高清观看| av天堂中文字幕网| 国产精品亚洲av一区麻豆| 日本精品一区二区三区蜜桃| 欧美不卡视频在线免费观看| 免费在线观看成人毛片| 欧美黑人欧美精品刺激| 99久久九九国产精品国产免费| 桃色一区二区三区在线观看| 天堂动漫精品| 国语自产精品视频在线第100页| 午夜激情福利司机影院| 久久久久久久久久黄片| 国产精品一区二区三区四区免费观看 | 日韩欧美三级三区| 日日夜夜操网爽| 超碰av人人做人人爽久久| 国产毛片a区久久久久| 18美女黄网站色大片免费观看| 全区人妻精品视频| 成人性生交大片免费视频hd| 久久久久亚洲av毛片大全| 老女人水多毛片| 成人午夜高清在线视频| 日韩精品青青久久久久久| 亚洲人成网站高清观看| 色播亚洲综合网| 2021天堂中文幕一二区在线观| 脱女人内裤的视频| 欧美成人a在线观看| 亚洲自拍偷在线| a级一级毛片免费在线观看| 午夜免费激情av| 欧美在线黄色| 在线播放无遮挡| 午夜精品久久久久久毛片777| 99久久精品热视频|