• <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孫偉喬煒周凱魏致富姚鵬
    城市勘測(cè) 2016年4期
    關(guān)鍵詞:區(qū)域模型

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

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

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

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

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

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

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

    1 引 言

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    3.1電離層TEC建模

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

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

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

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

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

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

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

    3.2電離層TEC變化分析

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

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

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

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

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

    4 結(jié) 語

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

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

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

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

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

    [4]楊凱,薛軍琛,宋淑麗等.基于CMONOC建立和評(píng)估中國區(qū)域電離層模型[J].大地測(cè)量與地球動(dòng)力學(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)測(cè)與延遲改正研究[D].上海:中國科學(xué)院上海天文臺(tái).

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

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

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

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

    [13]霍星亮,袁運(yùn)斌,歐吉坤.基于GPS資料研究中國區(qū)域電離層TEC的周日變化、半年度及冬季異?,F(xiàn)象[J].自然科學(xué)進(jìn)展,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—),男,博士,高級(jí)工程師,主要從事GNSS應(yīng)用研究、城市基礎(chǔ)控制測(cè)量等工作。

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

    猜你喜歡
    區(qū)域模型
    一半模型
    永久基本農(nóng)田集中區(qū)域“禁廢”
    分割區(qū)域
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    關(guān)于四色猜想
    分區(qū)域
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    基于嚴(yán)重區(qū)域的多PCC點(diǎn)暫降頻次估計(jì)
    免费看十八禁软件| 自线自在国产av| 777久久人妻少妇嫩草av网站| 91九色精品人成在线观看| 免费一级毛片在线播放高清视频 | 日日摸夜夜添夜夜添小说| 久久热在线av| 久久草成人影院| 无人区码免费观看不卡| 精品熟女少妇八av免费久了| 午夜福利在线观看吧| 国产精品av久久久久免费| av有码第一页| 一二三四社区在线视频社区8| 午夜a级毛片| 一级a爱片免费观看的视频| 欧美av亚洲av综合av国产av| 日韩精品青青久久久久久| 男人舔女人的私密视频| 精品一区二区三区视频在线观看免费| 美女国产高潮福利片在线看| 国内精品久久久久精免费| 最新美女视频免费是黄的| 久久影院123| 亚洲情色 制服丝袜| 亚洲情色 制服丝袜| 成人18禁高潮啪啪吃奶动态图| 高清在线国产一区| www.熟女人妻精品国产| 色综合欧美亚洲国产小说| 精品久久蜜臀av无| 成人三级做爰电影| 成熟少妇高潮喷水视频| 成人三级做爰电影| www.熟女人妻精品国产| 中文字幕最新亚洲高清| 天堂影院成人在线观看| xxx96com| 精品高清国产在线一区| 欧美国产精品va在线观看不卡| 看免费av毛片| 欧美国产日韩亚洲一区| 久久午夜综合久久蜜桃| 多毛熟女@视频| 亚洲国产精品久久男人天堂| 午夜久久久久精精品| 精品第一国产精品| 国产成人av激情在线播放| 黄片播放在线免费| 正在播放国产对白刺激| 精品第一国产精品| 国产激情欧美一区二区| 亚洲av成人一区二区三| 男男h啪啪无遮挡| 美国免费a级毛片| 女性被躁到高潮视频| 一级a爱视频在线免费观看| 精品国产超薄肉色丝袜足j| 欧美一区二区精品小视频在线| 纯流量卡能插随身wifi吗| 日韩欧美国产一区二区入口| 久久久久亚洲av毛片大全| 人人妻,人人澡人人爽秒播| 女人高潮潮喷娇喘18禁视频| 久久国产精品男人的天堂亚洲| 他把我摸到了高潮在线观看| 91精品三级在线观看| 久久久久久国产a免费观看| 久久中文字幕人妻熟女| 日本免费a在线| 亚洲精品美女久久久久99蜜臀| 日日摸夜夜添夜夜添小说| 国产日韩一区二区三区精品不卡| 国产精品二区激情视频| 女人高潮潮喷娇喘18禁视频| 麻豆成人av在线观看| 天天躁夜夜躁狠狠躁躁| av中文乱码字幕在线| 巨乳人妻的诱惑在线观看| 国产精品免费一区二区三区在线| 久久久久久大精品| 欧美在线一区亚洲| 久久久国产成人免费| 99久久精品国产亚洲精品| 日本三级黄在线观看| 看免费av毛片| 亚洲专区字幕在线| 黑人巨大精品欧美一区二区mp4| 亚洲成人久久性| 亚洲国产看品久久| 色综合亚洲欧美另类图片| 中文字幕人妻丝袜一区二区| 精品久久久久久久毛片微露脸| 18禁美女被吸乳视频| 久久亚洲真实| 禁无遮挡网站| 97超级碰碰碰精品色视频在线观看| 成人亚洲精品av一区二区| 淫妇啪啪啪对白视频| 黑人操中国人逼视频| 国产成+人综合+亚洲专区| 一级a爱视频在线免费观看| 午夜免费成人在线视频| 久久国产精品男人的天堂亚洲| 精品乱码久久久久久99久播| 国产精品精品国产色婷婷| 亚洲成av人片免费观看| 成人三级做爰电影| 脱女人内裤的视频| av电影中文网址| 夜夜看夜夜爽夜夜摸| 欧美+亚洲+日韩+国产| 国产又色又爽无遮挡免费看| 日本 av在线| 国产三级黄色录像| 亚洲国产精品成人综合色| 欧美成狂野欧美在线观看| 性少妇av在线| 免费久久久久久久精品成人欧美视频| 精品一区二区三区av网在线观看| 国产精品国产高清国产av| 精品国产乱子伦一区二区三区| 中文字幕av电影在线播放| 欧美 亚洲 国产 日韩一| 久久国产精品影院| 久久欧美精品欧美久久欧美| 久久狼人影院| 一二三四社区在线视频社区8| 黄频高清免费视频| www.www免费av| 国产成人av教育| 日本三级黄在线观看| 一级毛片精品| 一二三四在线观看免费中文在| 亚洲精品在线观看二区| 在线永久观看黄色视频| 午夜福利免费观看在线| 亚洲人成电影观看| 欧美丝袜亚洲另类 | 男人舔女人的私密视频| 色在线成人网| 欧美激情 高清一区二区三区| 国产精品自产拍在线观看55亚洲| 亚洲精品中文字幕一二三四区| 日本 欧美在线| 日韩一卡2卡3卡4卡2021年| 国产av一区在线观看免费| 欧美最黄视频在线播放免费| 国产乱人伦免费视频| 男人舔女人的私密视频| 18禁国产床啪视频网站| 一夜夜www| 国产成人一区二区三区免费视频网站| 国产精品永久免费网站| 亚洲最大成人中文| 大型黄色视频在线免费观看| 搡老岳熟女国产| 两个人免费观看高清视频| 国产精品九九99| 一级,二级,三级黄色视频| 国产区一区二久久| 免费在线观看日本一区| 精品不卡国产一区二区三区| 91成人精品电影| 久久精品影院6| www.www免费av| 成人亚洲精品av一区二区| 久久精品国产综合久久久| 夜夜看夜夜爽夜夜摸| 日韩欧美国产一区二区入口| 成人特级黄色片久久久久久久| 国产又爽黄色视频| 视频区欧美日本亚洲| 亚洲熟妇熟女久久| 一夜夜www| 亚洲欧美一区二区三区黑人| 这个男人来自地球电影免费观看| 免费不卡黄色视频| 国产99久久九九免费精品| 久久久久国产精品人妻aⅴ院| 大型黄色视频在线免费观看| 精品乱码久久久久久99久播| 黄频高清免费视频| 国产欧美日韩综合在线一区二区| 久久久久久久久免费视频了| 9热在线视频观看99| 村上凉子中文字幕在线| 一二三四社区在线视频社区8| 在线观看午夜福利视频| 精品久久蜜臀av无| 一级毛片高清免费大全| 亚洲精品国产区一区二| 国产av一区二区精品久久| 搡老妇女老女人老熟妇| 在线观看日韩欧美| 成人18禁高潮啪啪吃奶动态图| 国产精品国产高清国产av| 精品国产一区二区三区四区第35| 91av网站免费观看| 女警被强在线播放| 99精品在免费线老司机午夜| 久久精品91蜜桃| 国产成年人精品一区二区| 色播亚洲综合网| 色av中文字幕| 国产伦一二天堂av在线观看| www.自偷自拍.com| 欧美激情极品国产一区二区三区| 中国美女看黄片| 免费久久久久久久精品成人欧美视频| 国产伦一二天堂av在线观看| 日韩免费av在线播放| 国产成年人精品一区二区| 久久久久九九精品影院| 日本在线视频免费播放| 不卡av一区二区三区| 日本五十路高清| 亚洲无线在线观看| 如日韩欧美国产精品一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 不卡av一区二区三区| av中文乱码字幕在线| 国产一区二区在线av高清观看| 在线观看免费日韩欧美大片| 午夜久久久久精精品| 一级黄色大片毛片| 色在线成人网| 亚洲人成77777在线视频| 国产成人av教育| 久久国产亚洲av麻豆专区| 国产一区二区三区视频了| 精品人妻在线不人妻| av片东京热男人的天堂| 欧美一级a爱片免费观看看 | 在线观看免费日韩欧美大片| 亚洲欧美激情综合另类| 国产精品一区二区在线不卡| 久久久久国产精品人妻aⅴ院| 久久狼人影院| av网站免费在线观看视频| 久久草成人影院| 久久 成人 亚洲| 好男人电影高清在线观看| 亚洲精品美女久久av网站| 欧美成狂野欧美在线观看| 日本在线视频免费播放| 手机成人av网站| 亚洲国产精品999在线| 男女之事视频高清在线观看| aaaaa片日本免费| 看免费av毛片| 午夜精品国产一区二区电影| 婷婷六月久久综合丁香| 国产av精品麻豆| 国产麻豆成人av免费视频| 国产成人免费无遮挡视频| 不卡一级毛片| 欧美激情 高清一区二区三区| 91国产中文字幕| 免费在线观看日本一区| 12—13女人毛片做爰片一| 黄色女人牲交| 欧美日韩乱码在线| 国产成人欧美在线观看| 老司机深夜福利视频在线观看| 国产成人系列免费观看| 性色av乱码一区二区三区2| 性欧美人与动物交配| 在线观看舔阴道视频| 人人妻人人澡欧美一区二区 | 亚洲第一电影网av| 亚洲专区国产一区二区| 黄色女人牲交| 嫁个100分男人电影在线观看| 欧美国产日韩亚洲一区| 久久国产精品影院| 一进一出抽搐gif免费好疼| 99国产精品免费福利视频| 中文字幕人妻熟女乱码| 在线国产一区二区在线| 精品人妻1区二区| а√天堂www在线а√下载| 国产激情久久老熟女| 国产亚洲av高清不卡| 黄色视频不卡| 操美女的视频在线观看| www日本在线高清视频| 国产主播在线观看一区二区| 18禁裸乳无遮挡免费网站照片 | av片东京热男人的天堂| 超碰成人久久| 国产91精品成人一区二区三区| 日本撒尿小便嘘嘘汇集6| 欧美精品亚洲一区二区| 天堂√8在线中文| 成熟少妇高潮喷水视频| 国产一区二区三区视频了| 精品久久久久久久毛片微露脸| 精品国产一区二区久久| 无限看片的www在线观看| 久久婷婷人人爽人人干人人爱 | 搞女人的毛片| 一区福利在线观看| www.www免费av| 人妻丰满熟妇av一区二区三区| 国产欧美日韩一区二区三区在线| 精品久久久久久久毛片微露脸| 首页视频小说图片口味搜索| 欧美激情久久久久久爽电影 | 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人久久性| 极品教师在线免费播放| 亚洲国产精品久久男人天堂| 久久国产精品影院| 午夜福利,免费看| 一进一出抽搐动态| 又黄又粗又硬又大视频| 国产三级在线视频| 夜夜爽天天搞| 高潮久久久久久久久久久不卡| 国产三级在线视频| 久久精品91无色码中文字幕| 一级毛片女人18水好多| ponron亚洲| 18禁美女被吸乳视频| 国产av一区在线观看免费| 日韩免费av在线播放| 久久久久久久精品吃奶| 人人澡人人妻人| 人人妻人人澡欧美一区二区 | 亚洲专区国产一区二区| 久久精品影院6| 亚洲成a人片在线一区二区| 欧美精品亚洲一区二区| 老汉色av国产亚洲站长工具| 99久久精品国产亚洲精品| 国产精品99久久99久久久不卡| 变态另类丝袜制服| 欧美黄色淫秽网站| 亚洲欧美激情在线| 久久精品人人爽人人爽视色| 欧美中文综合在线视频| 精品国产超薄肉色丝袜足j| 制服丝袜大香蕉在线| 伦理电影免费视频| 大型黄色视频在线免费观看| 久久精品91无色码中文字幕| 男男h啪啪无遮挡| 丰满的人妻完整版| 别揉我奶头~嗯~啊~动态视频| 在线十欧美十亚洲十日本专区| 免费高清视频大片| 性欧美人与动物交配| 777久久人妻少妇嫩草av网站| 女性生殖器流出的白浆| 欧美中文综合在线视频| 999久久久国产精品视频| 制服丝袜大香蕉在线| 淫妇啪啪啪对白视频| 亚洲专区中文字幕在线| 大型av网站在线播放| 精品一区二区三区av网在线观看| 亚洲少妇的诱惑av| 日韩大码丰满熟妇| 国产精品永久免费网站| ponron亚洲| 亚洲精品中文字幕一二三四区| 国产精品电影一区二区三区| 人人妻人人澡欧美一区二区 | 亚洲欧美精品综合一区二区三区| 十八禁网站免费在线| 久久久国产成人精品二区| 国产精品亚洲一级av第二区| 国产亚洲欧美在线一区二区| 亚洲av电影不卡..在线观看| 12—13女人毛片做爰片一| 欧美日韩乱码在线| 黄片播放在线免费| 97人妻精品一区二区三区麻豆 | 黄色视频,在线免费观看| 丰满人妻熟妇乱又伦精品不卡| 婷婷精品国产亚洲av在线| 中亚洲国语对白在线视频| 91在线观看av| 宅男免费午夜| 最新在线观看一区二区三区| 亚洲 欧美 日韩 在线 免费| 国内毛片毛片毛片毛片毛片| 日韩有码中文字幕| 搞女人的毛片| 变态另类丝袜制服| 国产av一区在线观看免费| 国产三级在线视频| 1024视频免费在线观看| 免费一级毛片在线播放高清视频 | 成人亚洲精品av一区二区| 纯流量卡能插随身wifi吗| 国产精品久久久久久人妻精品电影| 精品国内亚洲2022精品成人| 亚洲片人在线观看| 99久久国产精品久久久| 亚洲欧美日韩另类电影网站| 久久久精品欧美日韩精品| 亚洲欧洲精品一区二区精品久久久| 人人妻,人人澡人人爽秒播| 1024视频免费在线观看| 国产熟女xx| 老汉色av国产亚洲站长工具| √禁漫天堂资源中文www| 日本免费一区二区三区高清不卡 | 一级a爱视频在线免费观看| 极品教师在线免费播放| 精品一区二区三区av网在线观看| 午夜福利18| 国产熟女午夜一区二区三区| 亚洲自偷自拍图片 自拍| 女生性感内裤真人,穿戴方法视频| 多毛熟女@视频| 久久久国产欧美日韩av| 国产欧美日韩一区二区精品| 又黄又爽又免费观看的视频| 久久久久久国产a免费观看| 中文字幕最新亚洲高清| 国产成人系列免费观看| 女警被强在线播放| 黄网站色视频无遮挡免费观看| 香蕉丝袜av| 丝袜美足系列| 深夜精品福利| 自拍欧美九色日韩亚洲蝌蚪91| ponron亚洲| 1024视频免费在线观看| 女性被躁到高潮视频| 亚洲狠狠婷婷综合久久图片| 久久人妻av系列| 国产又爽黄色视频| 欧美丝袜亚洲另类 | 国产乱人伦免费视频| 午夜亚洲福利在线播放| 麻豆一二三区av精品| 深夜精品福利| 在线播放国产精品三级| 久久狼人影院| 无人区码免费观看不卡| 成人特级黄色片久久久久久久| 久久久久久久精品吃奶| 亚洲专区字幕在线| 长腿黑丝高跟| 国产精品 欧美亚洲| 欧美av亚洲av综合av国产av| 国产成人av教育| 丝袜在线中文字幕| 桃红色精品国产亚洲av| 777久久人妻少妇嫩草av网站| 国产精品av久久久久免费| 午夜精品在线福利| 成熟少妇高潮喷水视频| av天堂在线播放| 免费在线观看黄色视频的| 成人特级黄色片久久久久久久| 亚洲国产欧美一区二区综合| 波多野结衣巨乳人妻| 岛国视频午夜一区免费看| 在线播放国产精品三级| 亚洲国产看品久久| 少妇的丰满在线观看| 亚洲 国产 在线| 69精品国产乱码久久久| 亚洲三区欧美一区| 国产精品久久久久久人妻精品电影| 精品久久久精品久久久| 国产一区二区在线av高清观看| 黄片播放在线免费| 欧美老熟妇乱子伦牲交| 欧美国产日韩亚洲一区| avwww免费| 色综合欧美亚洲国产小说| 亚洲欧美日韩无卡精品| 国产精品秋霞免费鲁丝片| 每晚都被弄得嗷嗷叫到高潮| 91国产中文字幕| 一级a爱视频在线免费观看| 满18在线观看网站| 女人爽到高潮嗷嗷叫在线视频| 国产欧美日韩精品亚洲av| av天堂久久9| 久久精品国产清高在天天线| 国内久久婷婷六月综合欲色啪| 欧美黑人精品巨大| 亚洲精品国产一区二区精华液| 久久国产精品男人的天堂亚洲| 久久久久国产一级毛片高清牌| 亚洲成人精品中文字幕电影| 桃色一区二区三区在线观看| 黑人巨大精品欧美一区二区蜜桃| 精品熟女少妇八av免费久了| 又紧又爽又黄一区二区| 老鸭窝网址在线观看| 久久久久国产精品人妻aⅴ院| 久久精品国产亚洲av高清一级| 国产伦一二天堂av在线观看| 国产精品av久久久久免费| 久久久久久久久中文| 19禁男女啪啪无遮挡网站| 亚洲精品av麻豆狂野| 亚洲伊人色综图| 欧美色视频一区免费| 国产一区二区三区视频了| 欧美乱妇无乱码| 欧美亚洲日本最大视频资源| 国产色视频综合| 欧美日韩黄片免| 好男人在线观看高清免费视频 | 久久久久久人人人人人| 如日韩欧美国产精品一区二区三区| 午夜两性在线视频| 禁无遮挡网站| 亚洲中文字幕一区二区三区有码在线看 | 少妇粗大呻吟视频| 一级作爱视频免费观看| 悠悠久久av| 两性夫妻黄色片| 精品国产超薄肉色丝袜足j| 黄片大片在线免费观看| 在线观看免费视频日本深夜| 欧美成人一区二区免费高清观看 | 免费看a级黄色片| 日本vs欧美在线观看视频| 精品人妻1区二区| 国产极品粉嫩免费观看在线| 丝袜在线中文字幕| 午夜福利在线观看吧| 九色亚洲精品在线播放| 一级,二级,三级黄色视频| 中文亚洲av片在线观看爽| 制服丝袜大香蕉在线| 成人永久免费在线观看视频| 国产一区二区三区在线臀色熟女| 一区二区三区国产精品乱码| 琪琪午夜伦伦电影理论片6080| 一级毛片高清免费大全| 十八禁网站免费在线| 啦啦啦观看免费观看视频高清 | 国产aⅴ精品一区二区三区波| 久久精品国产清高在天天线| 制服人妻中文乱码| 亚洲,欧美精品.| 国产精品一区二区在线不卡| 日韩欧美三级三区| 黄色a级毛片大全视频| 免费高清视频大片| 黄色毛片三级朝国网站| 亚洲av电影在线进入| 精品久久久精品久久久| 18禁美女被吸乳视频| 欧美日韩瑟瑟在线播放| 一个人观看的视频www高清免费观看 | 一区在线观看完整版| 国产高清视频在线播放一区| 中文字幕人妻熟女乱码| 成人av一区二区三区在线看| 国产欧美日韩一区二区精品| ponron亚洲| 97超级碰碰碰精品色视频在线观看| 男女午夜视频在线观看| 国产97色在线日韩免费| 国产又爽黄色视频| 日本vs欧美在线观看视频| 一级a爱视频在线免费观看| 91av网站免费观看| 精品久久久久久久人妻蜜臀av | 在线天堂中文资源库| 亚洲欧洲精品一区二区精品久久久| 涩涩av久久男人的天堂| 午夜精品久久久久久毛片777| 美女国产高潮福利片在线看| 露出奶头的视频| 久久热在线av| 国产不卡一卡二| 精品一区二区三区四区五区乱码| 精品国产超薄肉色丝袜足j| 人成视频在线观看免费观看| 91九色精品人成在线观看| 国产精品一区二区三区四区久久 | 婷婷六月久久综合丁香| 国产区一区二久久| 日本免费一区二区三区高清不卡 | 国产成人免费无遮挡视频| 国产91精品成人一区二区三区| 久久久久国产一级毛片高清牌| 黑人巨大精品欧美一区二区蜜桃| 黄片大片在线免费观看| 99国产精品一区二区三区| 天天添夜夜摸| 成人欧美大片| 免费在线观看完整版高清| 少妇的丰满在线观看| 久久久久亚洲av毛片大全| 精品熟女少妇八av免费久了| 变态另类成人亚洲欧美熟女 | 成人18禁高潮啪啪吃奶动态图| 一a级毛片在线观看| av片东京热男人的天堂| 大型黄色视频在线免费观看| 欧美色欧美亚洲另类二区 | 久久这里只有精品19| 亚洲国产毛片av蜜桃av| 中文字幕色久视频| 婷婷六月久久综合丁香| 50天的宝宝边吃奶边哭怎么回事| 成人免费观看视频高清| 两性午夜刺激爽爽歪歪视频在线观看 | 国产成人啪精品午夜网站| 国产91精品成人一区二区三区| 成人欧美大片|