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

    基于JASON-1高度計(jì)的海況偏差校正參數(shù)模型

    2013-07-07 15:39:16李書光王云海苗洪利周曉光任浩然王桂忠
    關(guān)鍵詞:高度計(jì)交叉點(diǎn)海況

    李書光,王云海,,苗洪利,周曉光,任浩然,王桂忠,張 杰

    (1.中國石油大學(xué)理學(xué)院,山東青島 266580;2.國家海洋局第一海洋研究所海洋物理與遙感研究室,山東青島 266061; 3.中國海洋大學(xué)信息科學(xué)與工程學(xué)院,山東青島 266100)

    基于JASON-1高度計(jì)的海況偏差校正參數(shù)模型

    李書光1,王云海1,2,苗洪利3,周曉光3,任浩然3,王桂忠3,張 杰2

    (1.中國石油大學(xué)理學(xué)院,山東青島 266580;2.國家海洋局第一海洋研究所海洋物理與遙感研究室,山東青島 266061; 3.中國海洋大學(xué)信息科學(xué)與工程學(xué)院,山東青島 266100)

    基于JASON-1高度計(jì)研究海況偏差校正參數(shù)模型。采用交叉點(diǎn)處測(cè)量值差分的方法消除高度計(jì)中除海況偏差之外的噪聲信號(hào)。基于泰勒展開構(gòu)造以有效波高、風(fēng)速為變量的32種海況偏差參數(shù)模型,利用線性回歸的方法得到各模型的參數(shù)估計(jì)值,并對(duì)其進(jìn)行評(píng)價(jià)、篩選得到最優(yōu)模型。通過將該模型的海況偏差校正值與JASON-1的地球物理數(shù)據(jù)集(GDR)數(shù)據(jù)相對(duì)比,進(jìn)一步驗(yàn)證該參數(shù)模型的有效性。結(jié)果表明:該參數(shù)模型有效,可以用于JASON-1高度計(jì)的海況偏差校正。

    高度計(jì);海況偏差;參數(shù)模型;線性回歸;交叉點(diǎn)差分

    海況偏差由電磁偏差、斜偏差和追蹤偏差三部分組成[1],是衛(wèi)星高度計(jì)測(cè)量海表面高度主要的誤差源之一[2],主要有理論和經(jīng)驗(yàn)兩種研究方法。理論方法得到的海況偏差模型[3-12]通常只能用來做定性分析,并且其需要的參數(shù)(如偏斜系數(shù)、統(tǒng)計(jì)動(dòng)差和波譜)并不能從高度計(jì)數(shù)據(jù)中得到,理論模型并不實(shí)用[13]。經(jīng)驗(yàn)方法認(rèn)為海況偏差是高度計(jì)可測(cè)得的與海況相關(guān)的變量(如風(fēng)速v、有效波高Hsw、后向散射系數(shù)σ0以及它們的組合等)的函數(shù),其最大的優(yōu)勢(shì)就是可以利用高度計(jì)所測(cè)得的覆蓋全球、數(shù)量極大的數(shù)據(jù)集。主要有參數(shù)模型[1,14-18]和非參數(shù)模型[2,19-23]等。在JASON-1高度計(jì)中,海況偏差校正采用的是Gaspar等提出的非參數(shù)模型[20]。非參數(shù)模型無法計(jì)算超出建模所用海況數(shù)據(jù)范圍的海況下的海況偏差值,并且在建模數(shù)據(jù)密度稀少的海況區(qū)域,海況偏差估計(jì)精度較差[21]。參數(shù)模型相對(duì)非參數(shù)模型具有表達(dá)式外延性好、簡(jiǎn)單直觀、分析容易以及使用方便的優(yōu)點(diǎn)。筆者基于JASON-1高度計(jì)數(shù)據(jù)開展海況偏差校正參數(shù)模型的研究。

    1 原 理

    1.1 信號(hào)噪聲的處理

    研究海況偏差校正參數(shù)模型時(shí)所用的初始數(shù)據(jù)是高度計(jì)相對(duì)于參考橢球面的海表面高度(Hss)。未校正的海表面高度測(cè)量值H′ss可以表示為

    式中,B為海況偏差;hg為大地水準(zhǔn)面高度;hd為海表面動(dòng)力高度;ε'為除海況偏差之外的所有測(cè)高誤差之和,包括軌道誤差、儀器誤差和其他地球物理校正誤差。

    海況偏差只是高度計(jì)數(shù)據(jù)中一個(gè)小信號(hào),大地水準(zhǔn)面高度是數(shù)據(jù)中最大的信號(hào),海表面動(dòng)力高度也是高度計(jì)數(shù)據(jù)中一個(gè)很大的信號(hào)。通常采用對(duì)衛(wèi)星軌道地面軌跡交叉點(diǎn)處測(cè)量數(shù)據(jù)進(jìn)行差分的方法來消除大地水準(zhǔn)面信號(hào)和海表面動(dòng)力高度信號(hào)中的定常部分。對(duì)交叉點(diǎn)之間不同時(shí)間的測(cè)量值進(jìn)行差分,則公式(1)可表示為

    式中,ΔH′ss稱為交叉點(diǎn)處海表面高度的不符值,它消除了數(shù)據(jù)中的所有不隨時(shí)間改變的信號(hào)部分,如大地水準(zhǔn)面,動(dòng)力地形的定常部分等;Δhd為海表面動(dòng)力高度的變化量。在JASON-1數(shù)據(jù)集中,升軌和降軌在交叉點(diǎn)處平均時(shí)間間隔為3.5 d,hd在此短時(shí)間間隔內(nèi)的變化一般是很有限的。

    1.2 建模依據(jù)

    參數(shù)模型基于海況偏差B與有效波高Hsw線性相關(guān)這一基本假設(shè),則海況偏差參數(shù)模型可寫為

    式中,b為一個(gè)無量綱的負(fù)值,稱為海況偏差系數(shù);X為與海況相關(guān)的變量組成的向量;θ為參數(shù)向量。在以往的經(jīng)驗(yàn)研究中,X一般直接從高度計(jì)測(cè)得的與海況相關(guān)的變量中選取,即有效波高Hsw、風(fēng)速v或后向散射系數(shù)σ0以及它們的組合。由于v和σ0具有高度相關(guān)性(通常用σ0反演v),因此,兩者只能選其一。本研究選擇使用v。將b對(duì)Hsw和v進(jìn)行泰勒展開,則可得到海況偏差校正參數(shù)模型為

    公式(4)所表達(dá)的各模型均保留常數(shù)項(xiàng)a1,依次增加變量的個(gè)數(shù),則可以得到數(shù)量眾多的海況偏差校正參數(shù)模型。按照選取的變量名對(duì)模型進(jìn)行命名,例如,選取a1、a3和a6項(xiàng)的模型命名為M136。

    2 模型的建立

    2.1 數(shù)據(jù)集

    所用的數(shù)據(jù)集是從JASON-1的Ku波段cycle40~cycle140的地球物理數(shù)據(jù)集(GDR)中提取的交叉點(diǎn)數(shù)據(jù),提取交叉點(diǎn)采用的是分段擬合衛(wèi)星地面軌道軌跡的方法。對(duì)交叉點(diǎn)處的海表面高度測(cè)量值添加了儀器、干對(duì)流層、濕對(duì)流層、電離層、大氣逆壓、高頻振蕩、海洋潮汐、極潮、固體地球潮以及負(fù)荷潮等校正。并對(duì)數(shù)據(jù)集采用以下質(zhì)量控制:

    (1)去除指向角大于0.3°的數(shù)據(jù);

    (2)去除后向散射系數(shù)σ0<7 dB、σ0>20 dB的所有數(shù)據(jù);

    (3)去除有效波高Hsw>11 m的全部數(shù)據(jù)。

    經(jīng)過質(zhì)量控制后,共得到253152組未經(jīng)海況偏差校正的海表面高度、風(fēng)速、有效波高測(cè)量值的交叉點(diǎn)數(shù)據(jù)集(H′ss,Hsw,v)。根據(jù)模型形式對(duì)Hsw、v進(jìn)行組合并對(duì)其進(jìn)行差值,最終獲得確定模型系數(shù)所需要的(ΔH′ss,ΔHsw,Δv,Δ,Δv2,Δ(Hswv),Δ(v),…)數(shù)據(jù)集。

    2.2 模型參數(shù)的確定

    高度計(jì)數(shù)據(jù)中海況偏差信號(hào)可以表示為

    式中,εB為參數(shù)模型所不能解釋的海況偏差信號(hào)部分;Xi為公式(4)中的變量;ai為變量Xi對(duì)應(yīng)的參數(shù)。則公式(2)變?yōu)?/p>

    將所有誤差合并為偏差a0和平均值為零的噪聲ε的和,則式(6)可重新表示為

    假定噪聲ε滿足高斯-馬爾柯夫條件和正態(tài)分布條件,那么此問題成為一個(gè)典型的多元線性回歸問題。給定n組(ΔH′ss,ΔX)觀測(cè)值,則參數(shù)的標(biāo)準(zhǔn)線性最小二乘估計(jì)為

    如果ΔX和ε不相關(guān),則估計(jì)量無偏。

    2.3 模型評(píng)價(jià)及優(yōu)選

    在給定0.05的顯著性水平下,對(duì)所有模型進(jìn)行F、Fi和決定系數(shù)(R2)檢驗(yàn),若F、Fi和決定系數(shù)(R2)其中之一的值小于其臨界值,則淘汰該模型。其中

    式中,m為模型中變量的數(shù)量;n為數(shù)據(jù)數(shù)量;Rss為回歸平方和;Ess為殘差平方和;Pi為偏回歸平方和; Tss為總變差。

    按照上述方法發(fā)現(xiàn),當(dāng)海況偏差系數(shù)b展開到三階時(shí),其模型性能相對(duì)二階模型沒有得到提升,并且海況偏差計(jì)算值出現(xiàn)了正值,這與海況偏差為負(fù)值相悖。因此,將b的泰勒展開式限制在二階。則公式(4)可以簡(jiǎn)化為

    公式(9)所表達(dá)的模型均保留常數(shù)項(xiàng)a1,依次增加變量的個(gè)數(shù),則可以得到1個(gè)常數(shù)模型,5個(gè)二參數(shù)模型,10個(gè)三參數(shù)模型,10個(gè)四參數(shù)模型,5個(gè)五參數(shù)模型,1個(gè)六參數(shù)模型,共計(jì)32種形式。各模型的F、R2以及解釋方差值D見表1。解釋方差D定義為未經(jīng)海況偏差校正前交叉點(diǎn)海表面高度H′ss不符值的方差減去經(jīng)過海況偏差校正后海表面高度Hss不符值的方差,可以理解為交叉點(diǎn)海表面高度不符值的方差中海況偏差能夠解釋的部分,用其來評(píng)價(jià)海況偏差模型的有效性,其值越大說明模型越有效。

    表1 模型性能評(píng)價(jià)Table 1 Evaluation of model functions

    按照上述評(píng)價(jià)方法,經(jīng)檢驗(yàn),所有的模型對(duì)應(yīng)的F值和R2均大于其臨界值,因此這32個(gè)模型總體自變量對(duì)因變量的線性影響和線性回歸的效果都是顯著的;Fi值在表中沒有列出。除模型M1256的第3項(xiàng)(F3值為0.04)、M12356的第5項(xiàng)(F5值為2.03)對(duì)因變量作用不顯著外,其余模型的各項(xiàng)作用都顯著。而從表1中也可以看出,這兩個(gè)模型與去掉不顯著變量對(duì)應(yīng)的模型相比其性能基本無任何改觀,因此淘汰這兩個(gè)模型。

    在以上模型參數(shù)檢驗(yàn)的基礎(chǔ)上,對(duì)剩余的6組30個(gè)模型進(jìn)行優(yōu)選,優(yōu)選的標(biāo)準(zhǔn)是:(1)決定系數(shù)值最大;(2)解釋方差最大。按此標(biāo)準(zhǔn),從表1中可以看出,隨著變量數(shù)量的增加,模型的性能整體上有所提升。因而,研究得到的最優(yōu)模型是M123456,其R2

    2.4 最優(yōu)模型特征分析

    圖1和圖2分別給出了模型M123456的平均回值及D在30個(gè)模型中都是最大?;貧w得到它的參數(shù)為:歸殘差()與有效波高不符值ΔHsw、風(fēng)速不符值Δv關(guān)系的散點(diǎn)圖。隨ΔHsw和Δv的變化越小,說明模型中參數(shù)估計(jì)偏差越小。

    圖1 平均回歸殘差ˉε隨ΔHsw變化Fig.1 Variation of mean regression residuals with ΔHsw

    圖2 平均回歸殘差ˉε隨Δv變化Fig.2 Variation of mean regression residuals with Δv

    由圖1可以看到,絕大部分ε分布在-0.1~0.5 cm之間;ΔHsw在-4~4 m范圍內(nèi)波動(dòng)較小,都在-0.5~0.5 cm之間。由圖2可以看出,全部分布在0線附近。因此,模型的參數(shù)估計(jì)值的偏差較小。

    圖3 b隨Hsw和v的分布特征Fig.3 b as a function of Hswand v

    圖3是由模型M123456計(jì)算得到的海況偏差系數(shù)b乘以100后的值隨Hsw和v的分布特征圖。從中可以看到,b隨Hsw變化明顯。Hsw小于4 m時(shí),隨著Hsw的增大,b的絕對(duì)值從3.6增大到4. 7;Hsw大于4 m時(shí),隨Hsw的增大b的絕對(duì)值從4.7遞減到1.9。這與Witter和Chelton[24]基于GEOSAT數(shù)據(jù)研究得出的結(jié)論相近。b隨v的變化基本以13 m/s為分界點(diǎn),v小于13 m/s時(shí),隨v的增加,b的絕對(duì)值增大,v大于13 m/s時(shí),隨v的增加,b的絕對(duì)值減小,這與Rodriguez等[6]通過數(shù)值模擬試驗(yàn)得出的結(jié)論相近。

    3 模型驗(yàn)證

    式中,n為選取數(shù)據(jù)的數(shù)量;Bi為參數(shù)模型校正值; Byi為GDR數(shù)據(jù)中對(duì)應(yīng)的海況偏差校正值。

    經(jīng)計(jì)算,全局?jǐn)?shù)據(jù)的標(biāo)準(zhǔn)偏差S為2.39 cm,相對(duì)誤差為11.18%,說明參數(shù)模型計(jì)算值和GDR數(shù)據(jù)整體吻合度較好。

    (2)計(jì)算B和By兩者之差ΔB,繪制其概率密度分布圖(圖4)及隨有效波高和風(fēng)速的變化圖(圖5)。

    為了檢驗(yàn)參數(shù)模型的適用性,本文中將M123456模型海況偏差校正值與JASON-1的GDR第41~76cycle覆蓋一年的數(shù)據(jù)進(jìn)行比對(duì)研究。

    (1)計(jì)算全局?jǐn)M合標(biāo)準(zhǔn)偏差S,用來整體上判斷兩者的符合度。其中

    圖4 ΔB概率密度分布Fig.4 Probability density distribution of ΔB

    由圖4可以看出,ΔB主要集中在-4~1 cm之間,此范圍內(nèi)的數(shù)據(jù)占總數(shù)據(jù)的87.1%,說明B和By在絕大多數(shù)的數(shù)據(jù)內(nèi)符合度良好。

    圖5 ΔB隨Hsw和v的變化Fig.5 Variation of ΔB with Hswand v

    由圖5可看出,B和By在絕大多數(shù)的Hsw、v范圍內(nèi)符合度良好,尤其是Hsw在0~5 m范圍內(nèi),此海況條件下ΔB主要集中在-4~1 cm之間。

    4 結(jié)束語

    基于JASON-1數(shù)據(jù)研究了海況偏差校正參數(shù)模型,利用線性回歸的方法確定了各模型的參數(shù)估計(jì)值,經(jīng)評(píng)價(jià)、篩選得到最優(yōu)模型是M123456,該模型解釋方差為39.61 cm2。通過與JASON-1 GDR數(shù)據(jù)中采用非參數(shù)模型所得的海況偏差值比對(duì),得到整體標(biāo)準(zhǔn)偏差S為2.39 cm,在絕大多數(shù)海況條件下兩者符合度良好。因此,該模型可以用于JASON-1高度計(jì)的海況偏差校正。

    [1] GASPAR P,OGOR F,Le TRAON P Y,et al.Estimating the sea state bias of the TOPEX and POSEIDON altimeters from crossover differences[J].Journal of Geophysical Research,1994,99:24981-24994.

    [2] LABROUE S,GASPAR P,DORANDEU J,et al.Non parametric estimates of the sea state bias for the Jason-1 radar altimeter[J].Marine Geodesy,2004,27(3/4): 453-481.

    [3] JACKSON F C.The reflection of impulses from a nonlinear random sea[J].Journal of Geophysical Research, 1979,84:4939-4943.

    [4] BARRICK D E,LIPA B J.Analysis and interpretation of altimeter sea echo[J].Adv Geophys,1985,27:60-99.

    [5] SROKOSZ M A.On the joint distribution of surface elevation and slopes for a nonlinear random sea,with an application to radar altimetry[J].Journal of Geophysical Research,1986,91:995-1006.

    [6] RODRIGUEZ E,KIM Y,MARTIN J M.The effect of small-wave modulation on the electromagnetic bias[J]. Journal of Geophysical Research,1992,97:2379-2389.

    [7] GLAZMAN R,FABRIKANT A,SROKOSZ M.Numerical analysis of the sea state bias for satellite altimetry [J].Journal of Geophysical Research,1996,101:3789-3799.

    [8] ELFOUHAILY T,CHAPRON B,KATSAROS K,et al. Weakly nonlinear theory and sea state bias estimations [J].Journal of Geophysical Research,1999,104:7641-7647.

    [9] ELFOUHAILY T,THOMPSON D R,CHAPRON B,et al.Improved electromagnetic bias theory[J].Journal of Geophysical Research,2000,105:1299-1310.

    [10] WARNICK K F,MILLET F W,ARNOLD D V.Physical and geometrical optics for 2-D surface with power-law height spectra[J].IEEE Trans Antennas Propag, 2005,53(3):922-932.

    [11] MILLET F W,WARNIC K F,NAGEL J R,et al. Physical optics-based electromagnetic bias theory with surface height-slope cross-correlation and hydrodynamic modulation[J].IEEE Transactions on Geoscience and Remote Sensing,2006,44(6):1470-1483.

    [12] NAENNA Praphun,JOHNSON Joel T.A Monte Carlo study of altimeter pulse returns and the electromagnetic bias[J].IEEE Transactions on Geoscience and Remote Sensing,2010,48(8):3218-3224.

    [13] FU L L,CAZENAVE Anny.Satellite altimetry and earth sciences:a handbook of techniques and applications[M].Academic Press,2001:57-64.

    [14] WITTER D L,CHELTON D B.A Geosat altimeter wind speed algorithm and a method for altimeter wind speed algorithm development[J].Journal of Geophysical Research,1991a,96:8853-8860.

    [15] RAY Richard D,KOBLINSKY C J,CHESTER J.On the sea-state bias of the Geosat altimeter[J].Journal of Atmospheric and Oceanic Technology,1991,8:397-408.

    [16] FU L L,GLAZMAN R.The effect of the degree of wave development on the sea state bias in radar altimetry measurement[J].Journal of Geophysical Research, 1991,96:829-834.

    [17] RAY R D,KOBLINSKY C J,BECKLEY B D.On the effectiveness of Geosat altimeter corrections[J].Int J Remote Sens,1991,12:1979-1984.

    [18] CHELTON D B.The sea state bias in altimeter estimates of sea level from collinear analysis of TOPEX data [J].Journal of Geophysical Research,1994,99: 24995-25008.

    [19] GASPAR P,FLORENS J P.Estimation of the sea state bias in radar altimeter measurements of sea level:results from a new non parametric method[J].Journal of Geophysical Research,1998,103:15803-15814.

    [20] GASPAR P,LABROUE S,OGOR F,et al.Improving nonparametric estimates of the sea state bias in radar altimeter measurements of sea level[J].JAOT,2002, 19:1690-1707.

    [21] SCHARROO R,LILLIBRIDGE J.Non-parametric seastate bias models and their relevance to sea level change studies.Proceedings of the 2004 Envisat&ERS Symposium,Salzburg,Austria,September 6-10,2004[C/ OL].[2011-05-10]http://earth.esa.int/workshops/ salzburg04/papers_posters/2P05_1_scharroo_198r.pdf.

    [22] FENG Hui,YAO Shan,LI Lin-yuan,et al.Splinebased nonparametric estimation of the altimeter sea-state bias correction[J].IEEE Geoscience and Remote Sensing Letters,2010,7(3):577-581.

    [23] TRAN N,LABROUE S,PHILIPPS S,et al.Overview and update of the sea state bias corrections for the Jason-2,Jason-1 and TOPEX missions[J].Marine Geodesy, 2010,33(S1):348-362.

    [24] WITTER D L,CHELTON D B.An apparent wave height dependence in the sea-state bias in GEOSAT altimeter range measurements[J].Journal of Geophysical Research,1991b,96:8861-8867.

    (編輯 修榮榮)

    A parametric model of estimating sea state bias based on JASON-1 altimetry

    LI Shu-guang1,WANG Yun-hai1,2,MIAO Hong-li3,ZHOU Xiao-guang3, REN Hao-ran3,WANG Gui-zhong3,ZHANG Jie2

    (1.College of Science in China University of Petroleum,Qingdao 266580,China;
    2.The First Institute of Oceanography,SOA,Qingdao 266061,China;
    3.College of Information Science and Engineering in Ocean University of China,Qingdao 266100,China)

    Based on the data of JASON-1 altimetry,the parametric model of estimating the sea state bias(SSB)was studied. Non-SSB signals within the altimeter data were eliminated by means of forming differences between measurements taken at crossover points.According to Taylor expansion,32 parametric models of SSB were developed as the function of both the significant wave height and wave speed.The estimation values of each parametric model were derived from the linear regression, and then the optimal model was obtained via the process of evaluation and selection.Finally,the effectiveness of the parametric model was validated through comparing the SSB estimation of the model with the geophysical data records(GDR)of JASON-1.The results show that the parametric model is effective,and can be used for JASON-1 SSB correction.

    JASON-1 altimetry;sea state bias;parametric model;linear regression;difference at crossover points

    TP 722. 6;P 228.3

    A

    1673-5005(2013)02-0181-06

    10.3969/j.issn.1673-5005.2013.02.030

    2011-07-22

    海洋公益性科研專項(xiàng)(201105032-1);國家自然科學(xué)基金項(xiàng)目(41176157)

    李書光(1964-),男,教授,碩士,主要研究方向?yàn)槁晫W(xué)與彈性波理論和應(yīng)用。E-mail:lshguang@upc.edu.cn。

    猜你喜歡
    高度計(jì)交叉點(diǎn)海況
    基于漂流浮標(biāo)的南大洋衛(wèi)星高度計(jì)有效波高研究
    典型海況下艦載發(fā)射箱結(jié)構(gòu)強(qiáng)度仿真分析
    MIMU/GNSS/ODO/高度計(jì)/航姿儀組合導(dǎo)航微系統(tǒng)硬件設(shè)計(jì)
    航天控制(2020年5期)2020-03-29 02:10:36
    同化衛(wèi)星高度計(jì)觀測(cè)對(duì)CAS-ESM-C上層海洋溫度模擬的改進(jìn)
    圍棋棋盤的交叉點(diǎn)
    合成孔徑雷達(dá)高度計(jì)與傳統(tǒng)高度計(jì)精度比對(duì)分析與機(jī)載試驗(yàn)驗(yàn)證
    基于高中生命科學(xué)知識(shí)交叉點(diǎn)的教學(xué)方法研究
    惡劣海況下海洋石油116內(nèi)轉(zhuǎn)塔式FPSO裝配載優(yōu)化
    船海工程(2015年4期)2016-01-05 15:53:38
    極端海況下軟質(zhì)海崖大規(guī)模蝕退計(jì)算模型研究
    極限海況下單點(diǎn)系泊系統(tǒng)纜索動(dòng)張力研究
    国产一区二区三区av在线| 人妻少妇偷人精品九色| 男女国产视频网站| 国产伦精品一区二区三区视频9| 看非洲黑人一级黄片| 麻豆乱淫一区二区| 久久久久久久精品精品| 我的女老师完整版在线观看| 全区人妻精品视频| 日韩,欧美,国产一区二区三区| 免费av观看视频| 成人亚洲精品av一区二区| 少妇人妻精品综合一区二区| 青春草亚洲视频在线观看| eeuss影院久久| 日韩欧美精品v在线| 国产精品久久久久久久电影| 夜夜爽夜夜爽视频| 国产精品一二三区在线看| 国产成人a∨麻豆精品| 又粗又硬又长又爽又黄的视频| freevideosex欧美| 在线观看av片永久免费下载| 久久久午夜欧美精品| 久久久午夜欧美精品| 日韩中字成人| 亚洲国产精品专区欧美| 成人国产av品久久久| 国产精品成人在线| 久久久久久久久久久免费av| 日韩强制内射视频| 国产成人精品婷婷| 卡戴珊不雅视频在线播放| 毛片女人毛片| 亚洲成人一二三区av| 国产精品蜜桃在线观看| 日韩 亚洲 欧美在线| 天堂俺去俺来也www色官网| 大香蕉97超碰在线| 91精品一卡2卡3卡4卡| 免费黄网站久久成人精品| 国产精品99久久99久久久不卡 | 亚洲最大成人手机在线| 日韩成人av中文字幕在线观看| 偷拍熟女少妇极品色| 男女边吃奶边做爰视频| 五月玫瑰六月丁香| 草草在线视频免费看| 国产精品成人在线| 国产 一区精品| 国产真实伦视频高清在线观看| 日本与韩国留学比较| 亚洲精品久久久久久婷婷小说| 国产高清不卡午夜福利| 大香蕉久久网| 日韩,欧美,国产一区二区三区| 伊人久久精品亚洲午夜| 国产伦理片在线播放av一区| 涩涩av久久男人的天堂| 中文字幕亚洲精品专区| 久久久久国产精品人妻一区二区| 久久国内精品自在自线图片| 成人特级av手机在线观看| 国产精品一及| av国产免费在线观看| 三级男女做爰猛烈吃奶摸视频| 99久久九九国产精品国产免费| 成人亚洲精品一区在线观看 | 日日摸夜夜添夜夜爱| 中国三级夫妇交换| 精品人妻视频免费看| 肉色欧美久久久久久久蜜桃 | 国产男人的电影天堂91| 国产精品无大码| 网址你懂的国产日韩在线| 国产伦精品一区二区三区四那| 中文字幕制服av| 亚洲精品一区蜜桃| 国产中年淑女户外野战色| 久久久久国产精品人妻一区二区| 麻豆精品久久久久久蜜桃| 插阴视频在线观看视频| 色播亚洲综合网| 久久韩国三级中文字幕| 插阴视频在线观看视频| 久久国内精品自在自线图片| 亚洲av免费高清在线观看| 欧美成人a在线观看| 亚洲av一区综合| 精品午夜福利在线看| 亚洲av.av天堂| 亚洲美女视频黄频| 久久久久久久久大av| 日韩大片免费观看网站| av国产免费在线观看| 三级国产精品欧美在线观看| 岛国毛片在线播放| 中文资源天堂在线| 日日摸夜夜添夜夜添av毛片| 中国美白少妇内射xxxbb| 国产白丝娇喘喷水9色精品| 女人久久www免费人成看片| 中国美白少妇内射xxxbb| 看黄色毛片网站| 欧美高清成人免费视频www| 亚洲人成网站在线观看播放| 国产综合懂色| 国产综合懂色| 97在线视频观看| 啦啦啦在线观看免费高清www| 激情五月婷婷亚洲| 国产欧美另类精品又又久久亚洲欧美| 亚洲av福利一区| 国产精品爽爽va在线观看网站| 日韩欧美精品v在线| 九九在线视频观看精品| 国产精品.久久久| 网址你懂的国产日韩在线| av福利片在线观看| 欧美人与善性xxx| av在线亚洲专区| 精品视频人人做人人爽| 欧美成人午夜免费资源| 97精品久久久久久久久久精品| 男插女下体视频免费在线播放| 色播亚洲综合网| 日韩成人av中文字幕在线观看| 免费大片18禁| 麻豆国产97在线/欧美| 激情五月婷婷亚洲| 午夜亚洲福利在线播放| 卡戴珊不雅视频在线播放| 一区二区三区免费毛片| 好男人在线观看高清免费视频| 亚洲欧美一区二区三区国产| 成年人午夜在线观看视频| 国模一区二区三区四区视频| 深爱激情五月婷婷| 嫩草影院精品99| 91久久精品国产一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 一级毛片电影观看| 久久久久久伊人网av| 少妇人妻一区二区三区视频| 日韩av在线免费看完整版不卡| 男人舔奶头视频| 街头女战士在线观看网站| 亚洲精品成人久久久久久| 神马国产精品三级电影在线观看| 国产精品福利在线免费观看| 日韩中字成人| eeuss影院久久| 成人国产av品久久久| 成年版毛片免费区| 日韩国内少妇激情av| 黄色日韩在线| av在线蜜桃| 日本免费在线观看一区| 国产探花极品一区二区| 十八禁网站网址无遮挡 | 国产精品国产三级国产av玫瑰| 黄片无遮挡物在线观看| 国产 一区精品| 日韩不卡一区二区三区视频在线| 丰满人妻一区二区三区视频av| 综合色av麻豆| 男女边吃奶边做爰视频| 国产成人免费无遮挡视频| 久久久久精品性色| 国产乱来视频区| av黄色大香蕉| 午夜精品国产一区二区电影 | 一个人看视频在线观看www免费| 伊人久久精品亚洲午夜| 大话2 男鬼变身卡| 五月伊人婷婷丁香| 亚洲av福利一区| 精品久久久久久久末码| 精品久久久久久久末码| 午夜老司机福利剧场| 午夜爱爱视频在线播放| 国产淫片久久久久久久久| 久久久久久久国产电影| 国产免费一区二区三区四区乱码| 国产国拍精品亚洲av在线观看| 男女边吃奶边做爰视频| 亚洲自偷自拍三级| 亚洲国产欧美在线一区| 舔av片在线| 综合色丁香网| 丝袜美腿在线中文| 麻豆成人午夜福利视频| 婷婷色av中文字幕| 联通29元200g的流量卡| 成人午夜精彩视频在线观看| 亚洲精华国产精华液的使用体验| 亚洲无线观看免费| 人妻制服诱惑在线中文字幕| 国产永久视频网站| 日本色播在线视频| 久久久久久伊人网av| 久久久精品欧美日韩精品| 亚洲av福利一区| 超碰av人人做人人爽久久| 嫩草影院新地址| 亚洲高清免费不卡视频| 欧美成人一区二区免费高清观看| 国产老妇伦熟女老妇高清| 久久精品国产a三级三级三级| 亚洲丝袜综合中文字幕| 免费黄网站久久成人精品| 精品久久久久久久人妻蜜臀av| 国产av码专区亚洲av| 97人妻精品一区二区三区麻豆| 99久久精品热视频| 成人亚洲欧美一区二区av| 国内精品美女久久久久久| 日本黄大片高清| 日韩国内少妇激情av| 成人午夜精彩视频在线观看| tube8黄色片| 国产成人午夜福利电影在线观看| 少妇高潮的动态图| 夜夜看夜夜爽夜夜摸| 男女啪啪激烈高潮av片| 亚洲天堂av无毛| 久久久久久久久大av| 一级a做视频免费观看| 久久综合国产亚洲精品| 大片免费播放器 马上看| 日韩免费高清中文字幕av| 乱系列少妇在线播放| 中文字幕亚洲精品专区| 看黄色毛片网站| 欧美区成人在线视频| 精品人妻偷拍中文字幕| 日韩不卡一区二区三区视频在线| 亚洲国产欧美在线一区| 亚洲精品久久午夜乱码| 五月天丁香电影| 一本一本综合久久| 国产乱来视频区| 亚洲av在线观看美女高潮| 日韩成人伦理影院| 国产极品天堂在线| 免费av毛片视频| 国产一区二区亚洲精品在线观看| 久久国内精品自在自线图片| 秋霞在线观看毛片| 伊人久久国产一区二区| 精品久久久久久久人妻蜜臀av| 91午夜精品亚洲一区二区三区| 久久久久久久久大av| 超碰av人人做人人爽久久| 99久久精品国产国产毛片| 久久精品国产鲁丝片午夜精品| 久久鲁丝午夜福利片| 亚洲av日韩在线播放| 久久久成人免费电影| 少妇裸体淫交视频免费看高清| 久久久久九九精品影院| 国产精品久久久久久久久免| 夜夜爽夜夜爽视频| av.在线天堂| 97人妻精品一区二区三区麻豆| 亚洲自偷自拍三级| 黄色配什么色好看| 嘟嘟电影网在线观看| 精品久久国产蜜桃| 国产成人精品婷婷| 久久人人爽人人爽人人片va| 最新中文字幕久久久久| 91午夜精品亚洲一区二区三区| 在线天堂最新版资源| 18+在线观看网站| 国产高清三级在线| 女的被弄到高潮叫床怎么办| 婷婷色综合www| 国产精品女同一区二区软件| 日韩三级伦理在线观看| 美女主播在线视频| 日韩av在线免费看完整版不卡| 欧美最新免费一区二区三区| 成人综合一区亚洲| 国产久久久一区二区三区| 建设人人有责人人尽责人人享有的 | 国产国拍精品亚洲av在线观看| 欧美一区二区亚洲| 午夜福利视频1000在线观看| 麻豆乱淫一区二区| 国产精品麻豆人妻色哟哟久久| 国产黄片美女视频| 亚洲国产精品专区欧美| 一区二区三区免费毛片| 成人亚洲精品av一区二区| 国产精品99久久久久久久久| 精品熟女少妇av免费看| 日本wwww免费看| 午夜免费鲁丝| 狂野欧美白嫩少妇大欣赏| 久久这里有精品视频免费| 国产女主播在线喷水免费视频网站| 美女主播在线视频| 欧美另类一区| 一级毛片我不卡| 国产精品久久久久久精品古装| 99热这里只有精品一区| 成人综合一区亚洲| 三级男女做爰猛烈吃奶摸视频| 精品人妻视频免费看| 激情 狠狠 欧美| 国产人妻一区二区三区在| 日日啪夜夜爽| 日本熟妇午夜| 国产午夜精品久久久久久一区二区三区| 日韩av免费高清视频| 婷婷色综合www| 熟妇人妻不卡中文字幕| 80岁老熟妇乱子伦牲交| 欧美日韩一区二区视频在线观看视频在线 | 18+在线观看网站| 成人毛片a级毛片在线播放| 亚洲精品乱久久久久久| 最后的刺客免费高清国语| 亚洲av不卡在线观看| 国产精品99久久久久久久久| 午夜日本视频在线| 少妇的逼好多水| 亚洲av免费高清在线观看| 国产日韩欧美亚洲二区| 尾随美女入室| 天天躁夜夜躁狠狠久久av| 国产成年人精品一区二区| 精品久久久精品久久久| 国产一级毛片在线| 日韩国内少妇激情av| 久久精品国产自在天天线| 国产成人免费无遮挡视频| 久久精品国产鲁丝片午夜精品| 成人免费观看视频高清| 久久久久精品久久久久真实原创| 青青草视频在线视频观看| 久久人人爽人人片av| 亚洲av成人精品一区久久| 肉色欧美久久久久久久蜜桃 | 美女视频免费永久观看网站| 久久精品熟女亚洲av麻豆精品| 免费少妇av软件| 真实男女啪啪啪动态图| 免费av观看视频| 亚洲色图av天堂| 肉色欧美久久久久久久蜜桃 | 久久久精品欧美日韩精品| 在线观看国产h片| 欧美成人a在线观看| 精品一区二区三区视频在线| 久久精品国产亚洲网站| 蜜桃久久精品国产亚洲av| 丝袜喷水一区| 国产欧美日韩精品一区二区| 看免费成人av毛片| 国产黄频视频在线观看| 婷婷色麻豆天堂久久| 国国产精品蜜臀av免费| 好男人在线观看高清免费视频| 一级爰片在线观看| 国产精品一区二区性色av| 2018国产大陆天天弄谢| 国产淫片久久久久久久久| 日韩欧美一区视频在线观看 | 色播亚洲综合网| 亚洲国产精品国产精品| 看十八女毛片水多多多| eeuss影院久久| 在线观看免费高清a一片| 久久精品国产鲁丝片午夜精品| 色网站视频免费| 亚洲av.av天堂| 高清欧美精品videossex| tube8黄色片| 亚洲精品aⅴ在线观看| 中国国产av一级| 一级av片app| 天天一区二区日本电影三级| 少妇人妻一区二区三区视频| 欧美激情在线99| 亚洲精品中文字幕在线视频 | 欧美日韩一区二区视频在线观看视频在线 | 日韩av不卡免费在线播放| 观看美女的网站| 国产色爽女视频免费观看| 自拍偷自拍亚洲精品老妇| 国产高清有码在线观看视频| 欧美97在线视频| 久久精品国产亚洲av天美| 成人免费观看视频高清| 亚洲av一区综合| av在线app专区| 麻豆久久精品国产亚洲av| 校园人妻丝袜中文字幕| 人妻少妇偷人精品九色| 成人免费观看视频高清| 午夜亚洲福利在线播放| 赤兔流量卡办理| 日本wwww免费看| 国产亚洲一区二区精品| 国产老妇女一区| 欧美zozozo另类| 看黄色毛片网站| 午夜福利高清视频| 大香蕉97超碰在线| 欧美最新免费一区二区三区| 一级毛片久久久久久久久女| a级一级毛片免费在线观看| 99久久中文字幕三级久久日本| 亚洲精品国产av蜜桃| 国产精品成人在线| 黄色视频在线播放观看不卡| 色网站视频免费| av在线播放精品| 欧美日本视频| 777米奇影视久久| 99久久人妻综合| 成人黄色视频免费在线看| 最近2019中文字幕mv第一页| 夜夜看夜夜爽夜夜摸| 久久97久久精品| av在线亚洲专区| 麻豆国产97在线/欧美| 在线观看人妻少妇| 黄片wwwwww| 免费看av在线观看网站| av网站免费在线观看视频| 国产亚洲av嫩草精品影院| h日本视频在线播放| 观看美女的网站| 天堂网av新在线| 在线观看一区二区三区| 亚洲欧美日韩另类电影网站 | 在线观看人妻少妇| 最新中文字幕久久久久| 啦啦啦啦在线视频资源| 久久精品人妻少妇| 亚洲av免费高清在线观看| 国产精品一区二区在线观看99| 校园人妻丝袜中文字幕| 亚洲精华国产精华液的使用体验| 国产av不卡久久| 91精品国产九色| 亚洲av中文av极速乱| 在线免费十八禁| 日本-黄色视频高清免费观看| 91精品伊人久久大香线蕉| 十八禁网站网址无遮挡 | 国产成人精品久久久久久| 91在线精品国自产拍蜜月| 欧美老熟妇乱子伦牲交| 97在线视频观看| 欧美97在线视频| 亚州av有码| 熟妇人妻不卡中文字幕| 亚洲精品日韩av片在线观看| 免费在线观看成人毛片| 成人毛片a级毛片在线播放| 99久久精品国产国产毛片| 国产男女超爽视频在线观看| 女人十人毛片免费观看3o分钟| kizo精华| 成人特级av手机在线观看| 大又大粗又爽又黄少妇毛片口| 欧美性感艳星| 国产成人午夜福利电影在线观看| 91aial.com中文字幕在线观看| 色婷婷久久久亚洲欧美| 午夜免费观看性视频| 国产真实伦视频高清在线观看| 国产老妇伦熟女老妇高清| 国产毛片a区久久久久| 国产av不卡久久| av在线app专区| 亚洲自拍偷在线| 免费观看的影片在线观看| 看免费成人av毛片| 中文精品一卡2卡3卡4更新| 亚洲精品乱码久久久v下载方式| 日本-黄色视频高清免费观看| 国产黄片视频在线免费观看| 老师上课跳d突然被开到最大视频| 日本与韩国留学比较| 欧美性感艳星| 久久久久久久大尺度免费视频| 亚洲激情五月婷婷啪啪| 最后的刺客免费高清国语| 岛国毛片在线播放| 色哟哟·www| 亚洲国产精品国产精品| 国产真实伦视频高清在线观看| 九色成人免费人妻av| 成人鲁丝片一二三区免费| 91久久精品电影网| 亚洲va在线va天堂va国产| 人人妻人人看人人澡| 日韩av免费高清视频| 免费黄频网站在线观看国产| 国产亚洲5aaaaa淫片| 婷婷色综合www| 伦理电影大哥的女人| 午夜免费鲁丝| 一级毛片aaaaaa免费看小| 99久久中文字幕三级久久日本| 美女视频免费永久观看网站| 蜜臀久久99精品久久宅男| 一区二区三区四区激情视频| 国产成人午夜福利电影在线观看| 国产成年人精品一区二区| 高清在线视频一区二区三区| 免费观看a级毛片全部| 免费在线观看成人毛片| 大话2 男鬼变身卡| av卡一久久| 亚洲国产日韩一区二区| 中国国产av一级| 男女无遮挡免费网站观看| 国产精品一及| 久久精品国产亚洲av天美| 看非洲黑人一级黄片| 97精品久久久久久久久久精品| 日本-黄色视频高清免费观看| 69av精品久久久久久| 国产成人freesex在线| 又粗又硬又长又爽又黄的视频| 直男gayav资源| 日本午夜av视频| 久热久热在线精品观看| 2022亚洲国产成人精品| 国产成人aa在线观看| 夜夜看夜夜爽夜夜摸| 王馨瑶露胸无遮挡在线观看| 99久久九九国产精品国产免费| 久久综合国产亚洲精品| kizo精华| 亚洲av福利一区| 成人鲁丝片一二三区免费| 久久精品人妻少妇| 一本一本综合久久| 高清午夜精品一区二区三区| 最近手机中文字幕大全| 女人久久www免费人成看片| 嫩草影院新地址| 人妻制服诱惑在线中文字幕| 校园人妻丝袜中文字幕| 免费人成在线观看视频色| 九九久久精品国产亚洲av麻豆| 亚洲av在线观看美女高潮| 九九久久精品国产亚洲av麻豆| 中文字幕免费在线视频6| 在线观看一区二区三区| 晚上一个人看的免费电影| av免费观看日本| 欧美丝袜亚洲另类| 午夜福利在线在线| 三级国产精品片| 亚洲av欧美aⅴ国产| 黄色视频在线播放观看不卡| 成人毛片a级毛片在线播放| 丝袜脚勾引网站| 中文天堂在线官网| 黄色视频在线播放观看不卡| 22中文网久久字幕| 51国产日韩欧美| 亚洲欧美精品专区久久| 超碰97精品在线观看| 亚洲av中文av极速乱| 热re99久久精品国产66热6| 五月伊人婷婷丁香| 欧美极品一区二区三区四区| 国产精品无大码| 亚洲精品456在线播放app| 国产 一区精品| 国产91av在线免费观看| av国产免费在线观看| 韩国av在线不卡| 日韩人妻高清精品专区| 亚洲va在线va天堂va国产| 在线a可以看的网站| 91久久精品国产一区二区成人| 国产欧美亚洲国产| 91在线精品国自产拍蜜月| 日韩国内少妇激情av| 日韩成人av中文字幕在线观看| 国产精品99久久久久久久久| 久久人人爽av亚洲精品天堂 | 国产 精品1| 欧美精品人与动牲交sv欧美| 国产亚洲av嫩草精品影院| 国产精品久久久久久久电影| 嫩草影院精品99| 日日摸夜夜添夜夜添av毛片| 美女xxoo啪啪120秒动态图| 中文资源天堂在线| 国产精品精品国产色婷婷| 成人欧美大片| 搡女人真爽免费视频火全软件| 女的被弄到高潮叫床怎么办| eeuss影院久久| 内地一区二区视频在线| 国产精品av视频在线免费观看| 中文字幕免费在线视频6| 最新中文字幕久久久久| 亚洲不卡免费看| 一本色道久久久久久精品综合| 国产一区有黄有色的免费视频| 国产国拍精品亚洲av在线观看| 午夜免费观看性视频| 欧美一区二区亚洲|