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

    聲速剖面EOF重構(gòu)的實(shí)測數(shù)據(jù)采樣深度研究

    2021-08-02 06:30:58王子蘅王振杰聶志喜張遠(yuǎn)帆
    海洋科學(xué) 2021年6期
    關(guān)鍵詞:聲速波束方差

    王子蘅, 王振杰, , 聶志喜, , 張遠(yuǎn)帆

    (1. 中國石油大學(xué)(華東) 海洋與空間信息學(xué)院, 山東 青島, 266580; 2. 青島海洋科學(xué)與技術(shù)試點(diǎn)國家實(shí)驗(yàn)室海洋礦產(chǎn)資源評價與探測技術(shù)功能實(shí)驗(yàn)室, 山東 青島, 266071)

    聲波在水中的良好傳播特性使其成為獲取和傳遞水下信息最為有效的手段[1-3]。聲速是影響水下多波束系統(tǒng)作業(yè)精度的重要外部影響因素, 通過影響聲線跟蹤的精度, 最終影響到測深精度。多波束測深通過換能器實(shí)時接收其發(fā)射出的各波束經(jīng)海底反射和散射后返回的到達(dá)角和旅行時[3], 利用聲速剖面數(shù)據(jù),由公式計算得到不同波束點(diǎn)對應(yīng)的水深值。聲速剖面是聲速的垂直結(jié)構(gòu)分布, 海水的介質(zhì)特性導(dǎo)致聲波的傳播軌跡發(fā)生彎曲, 要獲取波束腳印的確切位置, 需要沿著波束的傳播路徑追蹤聲線, 計算波束腳印的水平位移和深度, 即聲線跟蹤[3]。因此, 聲速剖面的正確與否直接影響多波束測深結(jié)果的精度和可靠性[4]。

    海洋中的聲速剖面通常由以下兩種方法獲得:直接測量法和間接測量法[4]。直接測量法通過測量聲波在海水中已知的固定距離內(nèi)往返所需的時間或相位來計算聲速值。間接測量法利用溫鹽深測量系統(tǒng)(conductivity, temperature, depth, CTD) 測量海水的溫度、鹽度及深度等物理量, 再通過聲速經(jīng)驗(yàn)公式計算作業(yè)水域內(nèi)各采樣深度的聲速值, 進(jìn)而構(gòu)成聲速剖面。然而在深水大洋進(jìn)行實(shí)際作業(yè)時, 受到水深和洋流的影響, 停船投放CTD采集聲速剖面數(shù)據(jù)時間成本較高。采用走航式的拋棄式聲速剖面計雖然可以快速有效的獲得聲速剖面, 但其作業(yè)成本隨之大大增加。相較于以上方法, 聲速剖面重構(gòu)只需少量的實(shí)測聲速數(shù)據(jù)即可實(shí)現(xiàn)全海深聲速剖面的重構(gòu), 提高效率的同時還有效的降低了成本。

    由于海洋環(huán)境的復(fù)雜性, 聲速剖面難以用一個簡單的函數(shù)表示。Davis等[5]證明EOF是在最小均方差意義下反映聲速剖面最有效的基函數(shù); 孫文川等[6]論證了EOF重構(gòu)的聲速剖面具有較高的內(nèi)符合精度, 能較好地描述實(shí)際聲速剖面; 張志偉等[7]研究表明取前6階及以上的EOF表示聲速剖面可滿足多波束測深的精度要求; 張鎮(zhèn)邁等[8]采用部分深度的實(shí)測剖面, 并基于測區(qū)海域內(nèi)已有的數(shù)據(jù), 對全海深聲速剖面實(shí)現(xiàn)了重構(gòu);張孝首等[9]利用實(shí)測CTD數(shù)據(jù)通過有理擬合延伸和EOF方法實(shí)現(xiàn)了深海5 000 m聲速剖面場的重構(gòu); 李洪超等[10]對不同深度位置的部分聲速剖面的重構(gòu)效果進(jìn)行了探討; 張維等[11]在對殘缺樣本聲速合理外延的基礎(chǔ)上, 選取聲速剖面變化較劇烈深度的3個聲速值實(shí)現(xiàn)了剖面重構(gòu)。以上研究表明, 基于EOF方法利用部分實(shí)測聲速數(shù)據(jù)重構(gòu)全海深剖面是可行有效的, 但未對參與重構(gòu)的實(shí)測數(shù)據(jù)的采樣深度給出明確依據(jù)。

    本文將測區(qū)Argo浮標(biāo)的溫鹽深數(shù)據(jù)通過聲速經(jīng)驗(yàn)公式計算得到的聲速剖面作為測區(qū)歷史聲速剖面資料, 基于EOF空間函數(shù)的方差貢獻(xiàn)率選取K(EOF階次)個最大梯度對應(yīng)深度處的聲速值進(jìn)行聲速剖面重構(gòu), 在保證多波束測深精度的前提下實(shí)現(xiàn)了對全海深聲速剖面的重構(gòu)。

    1 聲速剖面EOF表示原理

    1.1 聲速剖面EOF分解

    經(jīng)驗(yàn)正交函數(shù)分析方法也稱特征向量分析或主成分分析, 是分析矩陣數(shù)據(jù)的特性、提取主要數(shù)據(jù)特征量的一種方法[12]。

    在某海域采集N個聲速剖面:C= [c(1),c(2),…,c(N)], 內(nèi)插成為M個垂直標(biāo)準(zhǔn)層, 得到聲速矩陣CM×N:

    式中, 每一列為一個聲速剖面的標(biāo)準(zhǔn)層插值, 每一行為所有剖面在同一深度的聲速。將N條聲速剖面

    擾動矩陣的協(xié)方差矩陣為:

    其中每個元素rij為:

    將RM×M特征分解, 有:

    式中,DM×M為特征值矩陣, 令特征值λi按從大到小順序排列:FM×N為特征值對應(yīng)的特征向量矩陣, 也就是EOF空間函數(shù), 可表示為:

    式中,f(k)(k= 1,2,…,N)為M×1的列向量, 即為確定的第k階EOF。特征值λi對應(yīng)的分量f(i)即為第i階EOF。特征值的大小代表了其對應(yīng)EOF空間函數(shù)對線性空間的影響權(quán)重[7], 即越大的特征值對應(yīng)的特征向量中包含的聲速重構(gòu)信息越豐富。

    1.2 聲速剖面EOF重構(gòu)

    測區(qū)內(nèi)的聲速剖面完成EOF分解后, 利用前幾階EOF即可完成測區(qū)內(nèi)任一聲速剖面的重構(gòu):

    式中,c(z)為重構(gòu)的聲速剖面;z為各層海水深度;為平均聲速;K為EOF階次;αi為重構(gòu)系數(shù);fi(z)為EOF空間函數(shù)。移項(xiàng)得:

    則重構(gòu)系數(shù)矩陣AK×N為:

    式中,ai= [ai(1),ai(2),…,ai(K)]T為每個聲速剖面對應(yīng)的重構(gòu)系數(shù)。得到重構(gòu)系數(shù)后, 即可按式(9)對測區(qū)內(nèi)任一聲速剖面進(jìn)行EOF重構(gòu)。

    2 聲速剖面EOF重構(gòu)的數(shù)據(jù)采樣深度選取方法

    本文給出了用于聲速剖面EOF重構(gòu)的數(shù)據(jù)采樣深度選取的方法步驟。算法流程圖見圖1, 詳細(xì)步驟和說明如下:

    (1) 從中國Argo實(shí)時資料中心(http://www.argo.org.cn/)獲取某海域溫鹽壓數(shù)據(jù), 預(yù)處理后利用聲速經(jīng)驗(yàn)公式計算得到采樣位置處的聲速剖面。目前國內(nèi)外較為認(rèn)可的聲速經(jīng)驗(yàn)公式有Chen-Millero,Wilson, Del Grosso及Leroy等, 各模型的溫度、鹽度、壓力的適用范圍有所差異。根據(jù)文獻(xiàn)[13]的研究, Del Grosso, Chen-Millero, C.C.Leroy、Coppens四個經(jīng)驗(yàn)聲速模型在全球海域具有較高的精度及適用性, 文獻(xiàn)[14]指出, 在350~1 000 m的深度上, Chen-Millero公式的誤差最小, 因此本文選用Chen-Millero 聲速經(jīng)驗(yàn)公式計算的聲速剖面作為樣本聲速剖面, 利用Akima插值[15]將聲速剖面數(shù)據(jù)內(nèi)插到垂直標(biāo)準(zhǔn)層。

    (2) 按照1.1節(jié)中的方法進(jìn)行EOF分解, 提取EOF函數(shù), 計算各階EOF的方差貢獻(xiàn)率和累計方差貢獻(xiàn)率。方差貢獻(xiàn)率[16]指的是每一階EOF包含聲速場信息的百分比。其計算公式為:

    式中,σi表示第i階方差貢獻(xiàn)率,λi表示協(xié)方差矩陣R經(jīng)降序排列后的第i個特征值。

    前K階EOF對區(qū)域的累計貢獻(xiàn)率可表示為:

    一般地, 如果w(K)的值大于等于95%, 則認(rèn)為利用前K階EOF能較好的表示該區(qū)域的主要信息[5]。

    (3) 計算測區(qū)內(nèi)每一聲速剖面每一深度層的聲速梯度, 然后計算平均分層聲速梯度:

    式中,ci,n和ci+1,n分別表示第n條聲速剖面第i層和第i+1層的聲速值,zi+1和zi+1,n分別表示第n條聲速剖面第i層和第i+1層的深度值。按照平均分層聲速梯度由大到小, 對聲速層進(jìn)行降序排列, 形成新的聲速層矩陣。

    (4) 從新構(gòu)建的聲速層矩陣中選取K個深度處的聲速值(即選取K個聲速變化最劇烈處的聲速值),結(jié)合K階EOF計算得到重構(gòu)系數(shù)矩陣AK×N, 重構(gòu)全海深聲速剖面。

    (5) 基于實(shí)測聲速剖面, 預(yù)設(shè)水深值H和初始入射角, 根據(jù)分層常梯度聲線跟蹤法[2], 以H作為約束迭代計算聲波的傳播時間T:

    式中,Y表示實(shí)測聲速剖面的層數(shù),θi和θi+1分別表示實(shí)測剖面第i層和第i+1層的入射角,gi表示實(shí)測剖面第i層的梯度,ci和ci+1分別表示實(shí)測剖面第i層和第i+1層的聲速, Δzi表示第i層的深度差。p為Snell系數(shù):

    (6) 用重構(gòu)的聲速剖面和T反算波束水深值[6]H′。聲波在重構(gòu)的聲速剖面第i層的傳播時間ti為:

    式中,iθ′和θi1+′分別表示重構(gòu)聲速剖面第i層和第i+1層的入射角,gi′表示重構(gòu)剖面第i層的梯度,ic′和ci1+′分別表示重構(gòu)剖面第i層和第i+1層的聲速。

    式中,J表示使等式成立的ti的個數(shù)。所以:

    (7) 計算測深誤差σ:

    定義有效波束比為滿足深度限差的波束占總波束的比率[17], 首先統(tǒng)計滿足0.25%水深限差要求的波束數(shù)數(shù), 然后判斷多波束測深的有效波束比是否達(dá)到100%, 如果是, 則輸出測深偏差, 實(shí)驗(yàn)結(jié)束,反之, 則將選取階數(shù)K加1, 重新完成重構(gòu)系數(shù)矩陣的計算以及之后的步驟, 直至符合滿足測深精度的要求。

    3 實(shí)驗(yàn)分析

    實(shí)驗(yàn)數(shù)據(jù)源于中國Argo資料實(shí)時中心發(fā)布的《全球Argo浮標(biāo)剖面觀測資料質(zhì)量再控制數(shù)據(jù)集》[18], 選取2006—2009年每年1月份的Argo浮標(biāo)數(shù)據(jù)。將153.5°W—155W°, 27.5°N—29°N的海區(qū)作為區(qū)域1;153.5°W—155.5W°, 22.5°N—24.5°N的海區(qū)作為區(qū)域2。去除不合格的剖面后, 區(qū)域1選取了18個剖面, 區(qū)域2選取了15個剖面。

    典型深海聲速剖面結(jié)構(gòu)如圖2所示[2], 主要劃分為三部分: 混合層(由表面層和季節(jié)躍變層構(gòu)成)、主躍層和深海等溫層。其中, 聲速在混合層和主躍層的變化情況復(fù)雜, 而在深海等溫層中, 聲速隨深度的增加呈趨勢穩(wěn)定的正梯度增加, 接近線性變化。因此本文重點(diǎn)對混合層和主躍層這兩部分聲速變化情況更為復(fù)雜的水層的聲速剖面進(jìn)行研究, 統(tǒng)一截取海面至海深1 000 m的深度之間作為本次實(shí)驗(yàn)的剖面范圍。

    圖2 深海典型聲速剖面圖Fig. 2 Typical deep-sea sound speed profile

    隨機(jī)選取區(qū)域內(nèi)1條剖面用作檢核重構(gòu)精度,其余剖面作為已知剖面數(shù)據(jù)。實(shí)驗(yàn)區(qū)域內(nèi)的Argo剖面浮標(biāo)分布如圖3(a)和(b)所示, 圖中“×”點(diǎn)代表已知的樣本剖面,“▲”代表用于檢核的剖面。

    圖3 實(shí)驗(yàn)海域內(nèi)浮標(biāo)剖面分布圖Fig. 3 Distribution map of the buoy profiles in the experimental sea area

    3.1 數(shù)據(jù)預(yù)處理

    (1) Argo剖面浮標(biāo)測量的是海水的溫度、鹽度和壓力數(shù)據(jù), 采用 Chen-Millero聲速經(jīng)驗(yàn)公式[19]計算得到相應(yīng)深度處的聲速值:

    式中,

    其中,T為溫度, ℃;S為鹽度;P為壓力, bar。

    (2) 將計算得到的聲速數(shù)據(jù)利用Akima插值函數(shù)內(nèi)插成1 m的等間距標(biāo)準(zhǔn)層。樣本聲速剖面結(jié)構(gòu)如圖4。

    圖4 聲速剖面結(jié)構(gòu)示意圖Fig. 4 Schematic diagrams of SSP

    3.2 確定參與重構(gòu)的數(shù)據(jù)深度及聲速值

    分別對兩個區(qū)域的樣本聲速剖面進(jìn)行EOF分解,得到測區(qū)EOF空間函數(shù)的方差貢獻(xiàn)率和累計方差貢獻(xiàn)率, 如表1所示??梢钥闯?, 在區(qū)域1中, 前5階EOF的累計方差貢獻(xiàn)率可達(dá)97.985%, 在區(qū)域2中,前5階EOF的累計方差貢獻(xiàn)率達(dá)95.813%, 均已大于95%, 因此利用前5階EOF能夠較好的表示當(dāng)前區(qū)域的主要信息。

    表1 聲速剖面不同階次EOF方差貢獻(xiàn)率Tab. 1 Variance contribution rate of the different order EOF of the SSP

    圖5 給出的是由兩個實(shí)驗(yàn)海域內(nèi)的樣本聲速剖面提取的1~5階EOF, 可以明顯看出, 第1階EOF空間函數(shù)的變化較小, 整體呈趨勢穩(wěn)定, 這表明前幾階EOF中包含了實(shí)驗(yàn)海域內(nèi)的聲速場的主要信息, 而2~5階EOF空間函數(shù)的變化明顯, 變化集中在淺層區(qū)域, 尤其是在越靠后的階次中, 這反映了聲速場相對于平均SSP的細(xì)節(jié)變化[7]。

    圖5 1~5階EOF空間函數(shù)Fig. 5 1~5 order EOFs

    計算每個樣本聲速剖面每層的梯度變化情況,得到該區(qū)域的平均分層聲速梯度及對應(yīng)的深度位置, 按照平均分層聲速梯度由大到小對聲速層降序排列, 自上而下選取5個深度處的聲速值, 選取結(jié)果如表2所示。

    表2 選取的深度及聲速Tab. 2 Selected depths and sound speed values

    3.3 聲速剖面重構(gòu)

    選取表2提取的聲速值計算重構(gòu)系數(shù)矩陣, 進(jìn)而重構(gòu)全海深聲速剖面。兩個實(shí)驗(yàn)海域內(nèi)聲速剖面的重構(gòu)結(jié)果如圖6所示, 重構(gòu)誤差如圖7所示, 區(qū)域1的均方根誤差為0.514 m/s, 區(qū)域2的均方根誤差為0.609 m/s??梢钥闯觯?重構(gòu)聲速剖面與實(shí)測聲速剖面整體符合程度較好; 在兩個區(qū)域的淺層部分, 尤其是區(qū)域1的100 m深度附近和區(qū)域2的350 m深度附近,重構(gòu)聲速剖面與實(shí)測聲速剖面有明顯出入, 這可能是由于只選取了前幾階的EOF空間函數(shù)進(jìn)行重構(gòu), 部分階次的EOF空間函數(shù)并不能包含測區(qū)聲速場的全部信息, 造成了部分信息的損失, 因而出現(xiàn)了一定的偏差。此外, 以500 m水深作為分界, 深層部分的聲速剖面重構(gòu)結(jié)果較為穩(wěn)定, 整體的重構(gòu)誤差較小; 淺層部分的聲速剖面重構(gòu)結(jié)果相比較差, 造成該現(xiàn)象的原因可能是淺層部分的海水包含混合層和主躍層, 海洋要素在其中的變化情況復(fù)雜, 而在深層海水中已趨于穩(wěn)定。

    圖6 重構(gòu)的聲速剖面Fig. 6 Reconstructed SSP

    圖7 聲速剖面重構(gòu)誤差Fig. 7 Error of reconstructed SSP

    3.4 多波束測深檢驗(yàn)

    在多波束測深中, 利用EOF重構(gòu)的聲速剖面最終都要應(yīng)用于聲線改正。因此對于利用重構(gòu)的聲速剖面進(jìn)行聲線改正能否滿足實(shí)際應(yīng)用的要求, 還需要作進(jìn)一步的檢驗(yàn)。美國國家海洋與大氣局(national oceanic and atmospheric administration, NOAA)對聲速剖面引起的水深誤差規(guī)定[20]: “若使用的聲速剖面與實(shí)際聲速剖面對測深改正造成的互差超過0.25%水深則視為超限”。本文實(shí)驗(yàn)以NOAA的測深限差作為檢驗(yàn)聲速剖面重構(gòu)精度的標(biāo)準(zhǔn), 并結(jié)合有效波束比來判斷利用EOF重構(gòu)聲速剖面的有效性。設(shè)置水深和波束的初始入射角, 基于實(shí)測聲速剖面, 分層常梯度聲線跟蹤波束的往返時間, 利用該時間與重構(gòu)的聲速剖面再通過分層常梯度聲線跟蹤法反算水深, 并與真實(shí)水深值進(jìn)行對比, 比較兩者差值是否滿足0.25%水深限差, 最后統(tǒng)計有效波束比。

    表3 為使用兩則區(qū)域重構(gòu)的聲速剖面進(jìn)行多波束測深的結(jié)果, 可以看出, 使用重構(gòu)的聲速剖面進(jìn)行聲速改正得到的各波束點(diǎn)水深均符合限差要求,有效波束比均達(dá)到了100%。

    表3 重構(gòu)聲速剖面的最大偏差和有效波束比Tab. 3 Maximum deviations and effective beam ratios of the reconstructed SSPs

    4 結(jié)論

    本文對重構(gòu)全海深聲速剖面所需實(shí)測數(shù)據(jù)的采樣深度進(jìn)行了研究, 實(shí)驗(yàn)證明在對歷史聲速剖面資料進(jìn)行分析后, 按照K階EOF函數(shù)方差累計貢獻(xiàn)率大于95%, 選取K個最大平均分層聲速梯度處的聲速值作為已知值, 重構(gòu)的聲速剖面能夠滿足NOAA對多波束測量聲速剖面造成的誤差要求, 為實(shí)際測量作業(yè)中聲剖數(shù)據(jù)的采樣深度提供了參考, 提高了工作效率。由于海洋聲速易受各種環(huán)境復(fù)雜性的影響, 因此其他海區(qū)的選取階次及深度會有所不同, 但本文介紹的方法和所得規(guī)律性結(jié)論具有一定的普適性。

    猜你喜歡
    聲速波束方差
    方差怎么算
    概率與統(tǒng)計(2)——離散型隨機(jī)變量的期望與方差
    計算方差用哪個公式
    毫米波大規(guī)模陣列天線波束掃描研究*
    圓陣多波束測角探究
    電子測試(2018年6期)2018-05-09 07:31:54
    方差生活秀
    Helix陣匹配場三維波束形成
    聲速是如何測定的
    跨聲速風(fēng)洞全模顫振試驗(yàn)技術(shù)
    機(jī)翼跨聲速抖振研究進(jìn)展
    午夜精品一区二区三区免费看| 成人性生交大片免费视频hd| 中文在线观看免费www的网站| 麻豆成人av视频| 在线观看66精品国产| 欧美日韩精品成人综合77777| 久久久精品欧美日韩精品| 成人av在线播放网站| av在线老鸭窝| 中文亚洲av片在线观看爽| 欧美精品一区二区大全| 亚洲无线在线观看| 亚洲精品影视一区二区三区av| 日本与韩国留学比较| 国产亚洲精品久久久久久毛片| 国产精品野战在线观看| 能在线免费观看的黄片| 永久网站在线| 日本五十路高清| 变态另类成人亚洲欧美熟女| 夫妻性生交免费视频一级片| 在线a可以看的网站| 一个人看的www免费观看视频| 欧美日本亚洲视频在线播放| 国产人妻一区二区三区在| 天天躁日日操中文字幕| 天美传媒精品一区二区| 精品国产三级普通话版| 久久精品久久久久久噜噜老黄 | 日韩精品青青久久久久久| 国产av麻豆久久久久久久| 国产精品1区2区在线观看.| h日本视频在线播放| 国产一区二区在线观看日韩| 国产综合懂色| 日本爱情动作片www.在线观看| 亚洲欧美日韩东京热| 免费人成视频x8x8入口观看| eeuss影院久久| 日韩大尺度精品在线看网址| 六月丁香七月| 国内揄拍国产精品人妻在线| 日本熟妇午夜| 国产精品久久视频播放| 国产亚洲5aaaaa淫片| 国产麻豆成人av免费视频| 亚洲av男天堂| 欧美激情国产日韩精品一区| 18禁在线无遮挡免费观看视频| 国产精品综合久久久久久久免费| 少妇熟女欧美另类| 亚洲熟妇中文字幕五十中出| 伦理电影大哥的女人| 久久人人爽人人片av| 日日撸夜夜添| 女人十人毛片免费观看3o分钟| 精品人妻一区二区三区麻豆| 秋霞在线观看毛片| 成人毛片a级毛片在线播放| 国产精品久久视频播放| 国产午夜福利久久久久久| 日韩高清综合在线| 麻豆久久精品国产亚洲av| 美女 人体艺术 gogo| 国内揄拍国产精品人妻在线| 日本一本二区三区精品| 中文欧美无线码| 能在线免费观看的黄片| 免费看a级黄色片| 非洲黑人性xxxx精品又粗又长| 2022亚洲国产成人精品| 国产麻豆成人av免费视频| 亚州av有码| 99热这里只有是精品在线观看| 午夜福利高清视频| 日韩视频在线欧美| 一个人免费在线观看电影| 国产极品精品免费视频能看的| 国产精品伦人一区二区| 免费看美女性在线毛片视频| 国产午夜福利久久久久久| 精品久久久久久成人av| 男女下面进入的视频免费午夜| 男女那种视频在线观看| 国产精品精品国产色婷婷| 国产乱人视频| 成年av动漫网址| 国产精品一区二区三区四区久久| 亚洲三级黄色毛片| 久久人人爽人人爽人人片va| 看黄色毛片网站| 久久精品影院6| 一本精品99久久精品77| 日本与韩国留学比较| 99久久精品国产国产毛片| 高清午夜精品一区二区三区 | 校园春色视频在线观看| 国产精品久久久久久精品电影小说 | 日韩欧美在线乱码| 如何舔出高潮| 国模一区二区三区四区视频| 91久久精品国产一区二区成人| 青春草亚洲视频在线观看| 亚洲美女视频黄频| av免费观看日本| 我的老师免费观看完整版| 天堂中文最新版在线下载 | 有码 亚洲区| 色噜噜av男人的天堂激情| 久久久久免费精品人妻一区二区| 最近手机中文字幕大全| 国产真实伦视频高清在线观看| 小说图片视频综合网站| 亚洲av不卡在线观看| 亚洲国产精品久久男人天堂| 中出人妻视频一区二区| 国产v大片淫在线免费观看| 国产单亲对白刺激| 亚洲欧美中文字幕日韩二区| 禁无遮挡网站| 精品熟女少妇av免费看| 日本在线视频免费播放| 女的被弄到高潮叫床怎么办| 国产精品永久免费网站| 国产精品久久电影中文字幕| 国产av一区在线观看免费| 99热这里只有精品一区| 欧美丝袜亚洲另类| 久久国内精品自在自线图片| 久久久久九九精品影院| 成人无遮挡网站| 一本久久中文字幕| 国产片特级美女逼逼视频| 亚洲欧美中文字幕日韩二区| 非洲黑人性xxxx精品又粗又长| 三级毛片av免费| 久久人人爽人人片av| 在线免费观看不下载黄p国产| 99热6这里只有精品| 麻豆乱淫一区二区| a级毛色黄片| 亚洲人成网站在线播| 欧美色欧美亚洲另类二区| 又爽又黄无遮挡网站| 高清毛片免费看| 亚洲国产精品国产精品| 国产又黄又爽又无遮挡在线| 精品久久久噜噜| 人人妻人人看人人澡| 亚洲,欧美,日韩| 国语自产精品视频在线第100页| 欧美最新免费一区二区三区| 欧美日韩综合久久久久久| av福利片在线观看| 一个人看视频在线观看www免费| 看免费成人av毛片| 波野结衣二区三区在线| 国产亚洲精品av在线| 亚洲无线观看免费| www.av在线官网国产| 你懂的网址亚洲精品在线观看 | 中文字幕制服av| 亚洲国产精品久久男人天堂| 国产在线精品亚洲第一网站| 国产在线男女| 日韩成人伦理影院| 日韩,欧美,国产一区二区三区 | 久久精品国产鲁丝片午夜精品| 老司机福利观看| 九九爱精品视频在线观看| 高清毛片免费看| 国产欧美日韩精品一区二区| 99精品在免费线老司机午夜| 天堂av国产一区二区熟女人妻| 成人欧美大片| 午夜免费男女啪啪视频观看| 国产伦一二天堂av在线观看| 啦啦啦韩国在线观看视频| 日本黄色片子视频| 久久久国产成人精品二区| 中文字幕免费在线视频6| 亚洲欧美中文字幕日韩二区| 亚洲婷婷狠狠爱综合网| 最近中文字幕高清免费大全6| 日韩强制内射视频| 久久综合国产亚洲精品| 岛国在线免费视频观看| 国产精品国产三级国产av玫瑰| 少妇人妻一区二区三区视频| 国产三级中文精品| 三级男女做爰猛烈吃奶摸视频| 少妇高潮的动态图| 赤兔流量卡办理| 日韩欧美国产在线观看| 精品日产1卡2卡| 国产中年淑女户外野战色| 亚洲欧美日韩高清专用| 少妇裸体淫交视频免费看高清| 国产精品国产三级国产av玫瑰| 久久99精品国语久久久| 日韩一本色道免费dvd| 天天躁夜夜躁狠狠久久av| 深夜a级毛片| 免费av毛片视频| 男女做爰动态图高潮gif福利片| 国产一区二区三区av在线 | 亚洲精品国产成人久久av| 一级毛片久久久久久久久女| 日日摸夜夜添夜夜爱| 精品欧美国产一区二区三| 99riav亚洲国产免费| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品成人久久小说 | 日本-黄色视频高清免费观看| 少妇熟女aⅴ在线视频| 插逼视频在线观看| 天天躁日日操中文字幕| 中文精品一卡2卡3卡4更新| 18禁在线播放成人免费| 日韩高清综合在线| 啦啦啦啦在线视频资源| 亚洲人成网站在线播| 97超碰精品成人国产| 日韩欧美在线乱码| 最近最新中文字幕大全电影3| 亚洲丝袜综合中文字幕| 免费看光身美女| 99视频精品全部免费 在线| 亚洲av中文av极速乱| 久久国产乱子免费精品| 日韩欧美三级三区| 狠狠狠狠99中文字幕| 久久久国产成人免费| 亚洲国产精品国产精品| 午夜免费男女啪啪视频观看| 国产真实伦视频高清在线观看| 午夜亚洲福利在线播放| 好男人视频免费观看在线| 麻豆精品久久久久久蜜桃| 国内久久婷婷六月综合欲色啪| 国产午夜精品一二区理论片| 日本黄色视频三级网站网址| 国产极品天堂在线| 丰满人妻一区二区三区视频av| 国产高清有码在线观看视频| 久久久色成人| 九九爱精品视频在线观看| 精品久久久噜噜| 日韩,欧美,国产一区二区三区 | 欧美极品一区二区三区四区| 欧美最黄视频在线播放免费| 国产探花极品一区二区| 亚洲精品乱码久久久久久按摩| 最新中文字幕久久久久| 久久久久久久午夜电影| 最近手机中文字幕大全| 亚洲人成网站在线播| 亚洲人成网站在线观看播放| 日韩欧美 国产精品| 在线免费观看的www视频| 亚洲aⅴ乱码一区二区在线播放| 国产麻豆成人av免费视频| 亚洲美女搞黄在线观看| 国产伦一二天堂av在线观看| 两个人的视频大全免费| 中文欧美无线码| 18禁在线播放成人免费| 最近手机中文字幕大全| 欧美高清性xxxxhd video| 麻豆成人av视频| 日韩精品青青久久久久久| 亚洲自偷自拍三级| 久久6这里有精品| 一本久久精品| 成人鲁丝片一二三区免费| 免费人成在线观看视频色| 国产精品久久久久久久电影| 久久精品国产自在天天线| 亚洲国产欧美人成| 非洲黑人性xxxx精品又粗又长| a级一级毛片免费在线观看| 色哟哟哟哟哟哟| 久久亚洲精品不卡| 91在线精品国自产拍蜜月| 欧美色欧美亚洲另类二区| 十八禁国产超污无遮挡网站| 日本成人三级电影网站| 亚洲欧美成人综合另类久久久 | 国产成人91sexporn| 69人妻影院| 一级av片app| 国产成人一区二区在线| 日韩人妻高清精品专区| 我的女老师完整版在线观看| 国产高清有码在线观看视频| 久99久视频精品免费| 成人一区二区视频在线观看| 变态另类成人亚洲欧美熟女| 国产不卡一卡二| 一级毛片久久久久久久久女| av女优亚洲男人天堂| 级片在线观看| 九九在线视频观看精品| 久久精品国产清高在天天线| 欧美xxxx黑人xx丫x性爽| 久久久久久久久久久丰满| 伦理电影大哥的女人| 精华霜和精华液先用哪个| 91在线精品国自产拍蜜月| 亚洲最大成人手机在线| 免费看日本二区| 欧美激情国产日韩精品一区| 国产精品无大码| 中国国产av一级| 好男人在线观看高清免费视频| 欧美丝袜亚洲另类| 一本精品99久久精品77| 网址你懂的国产日韩在线| av免费观看日本| 麻豆国产97在线/欧美| 成人性生交大片免费视频hd| 啦啦啦观看免费观看视频高清| 亚洲欧美精品自产自拍| 亚洲国产欧美在线一区| 国产亚洲精品久久久com| 午夜免费激情av| 国产精品久久电影中文字幕| 成年版毛片免费区| 久久人人精品亚洲av| 一级毛片久久久久久久久女| 尤物成人国产欧美一区二区三区| 成年女人看的毛片在线观看| a级毛片免费高清观看在线播放| 能在线免费看毛片的网站| 欧美成人精品欧美一级黄| 亚洲性久久影院| 日本五十路高清| 色哟哟·www| 久久久色成人| 看黄色毛片网站| 国产老妇女一区| 日韩制服骚丝袜av| 亚洲欧美清纯卡通| 午夜久久久久精精品| 熟女电影av网| 丰满人妻一区二区三区视频av| 六月丁香七月| 99热这里只有是精品50| ponron亚洲| 黑人高潮一二区| 国产亚洲av嫩草精品影院| 亚洲天堂国产精品一区在线| 人体艺术视频欧美日本| 免费av毛片视频| 美女内射精品一级片tv| 嘟嘟电影网在线观看| 久久久久久久久中文| 91午夜精品亚洲一区二区三区| a级毛片a级免费在线| 日韩欧美精品免费久久| 亚洲av电影不卡..在线观看| 久久久久久久久久成人| 亚洲av中文av极速乱| 国产成年人精品一区二区| 99在线视频只有这里精品首页| 亚洲欧美日韩高清在线视频| 精品不卡国产一区二区三区| 久久精品国产亚洲av香蕉五月| 日韩欧美精品免费久久| 日本一本二区三区精品| 国产精品久久久久久久电影| 中文精品一卡2卡3卡4更新| 精品无人区乱码1区二区| 精品少妇黑人巨大在线播放 | 亚洲成人久久性| 成年女人永久免费观看视频| 夜夜夜夜夜久久久久| 国产成人freesex在线| 五月伊人婷婷丁香| 久久韩国三级中文字幕| 国产 一区 欧美 日韩| 22中文网久久字幕| 99九九线精品视频在线观看视频| 免费看av在线观看网站| 99久久九九国产精品国产免费| 国产精品久久久久久精品电影| 人妻少妇偷人精品九色| 白带黄色成豆腐渣| 国产一区二区三区av在线 | 国产成人午夜福利电影在线观看| 国产一区二区三区av在线 | 亚洲在久久综合| 国产私拍福利视频在线观看| 人妻少妇偷人精品九色| 在线国产一区二区在线| 久久人人精品亚洲av| 夜夜爽天天搞| 久久人妻av系列| 久久久午夜欧美精品| 毛片一级片免费看久久久久| 一级毛片aaaaaa免费看小| 欧美区成人在线视频| 亚洲乱码一区二区免费版| 97热精品久久久久久| 国产在视频线在精品| 国语自产精品视频在线第100页| 亚洲成人久久性| 中国国产av一级| 精品一区二区三区视频在线| 老女人水多毛片| 久久欧美精品欧美久久欧美| 日韩av在线大香蕉| 国产色爽女视频免费观看| 12—13女人毛片做爰片一| 最近中文字幕高清免费大全6| 老司机福利观看| 夫妻性生交免费视频一级片| 国产私拍福利视频在线观看| 久99久视频精品免费| 岛国毛片在线播放| 国产黄a三级三级三级人| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利在线观看免费完整高清在 | 亚洲国产欧洲综合997久久,| 国产黄色视频一区二区在线观看 | 久久久久免费精品人妻一区二区| 97超视频在线观看视频| 免费在线观看成人毛片| 一级黄色大片毛片| 久久久精品94久久精品| 亚洲图色成人| 日韩一本色道免费dvd| 能在线免费看毛片的网站| 亚洲国产日韩欧美精品在线观看| 久久久久九九精品影院| 日本欧美国产在线视频| 久久亚洲国产成人精品v| 欧美日韩精品成人综合77777| 一边摸一边抽搐一进一小说| 一级毛片aaaaaa免费看小| 深夜a级毛片| 三级毛片av免费| 国产极品精品免费视频能看的| 国产一区二区三区在线臀色熟女| 极品教师在线视频| av.在线天堂| 国产精品爽爽va在线观看网站| 精品久久久久久成人av| 不卡一级毛片| 免费观看人在逋| 久久久久久久亚洲中文字幕| 日韩中字成人| 国产精品久久久久久精品电影小说 | 此物有八面人人有两片| 91在线精品国自产拍蜜月| 午夜精品一区二区三区免费看| 少妇人妻一区二区三区视频| 午夜精品在线福利| 日韩国内少妇激情av| 波多野结衣高清作品| 免费看av在线观看网站| 熟女电影av网| 在线观看美女被高潮喷水网站| 亚洲第一电影网av| av天堂在线播放| 亚洲精品久久国产高清桃花| 久久久久久久久久成人| 波多野结衣高清作品| 中出人妻视频一区二区| 国产午夜福利久久久久久| 嫩草影院新地址| 欧美潮喷喷水| 国内揄拍国产精品人妻在线| 99在线视频只有这里精品首页| 九草在线视频观看| 男女边吃奶边做爰视频| 久久久久久久久久久丰满| 蜜臀久久99精品久久宅男| 国产亚洲精品久久久com| 国产精品一区二区三区四区免费观看| 又粗又硬又长又爽又黄的视频 | 精品久久久久久久人妻蜜臀av| 男人舔女人下体高潮全视频| 最近的中文字幕免费完整| 国产精品一区二区在线观看99 | 成人二区视频| 一区二区三区高清视频在线| 嫩草影院入口| 国产伦一二天堂av在线观看| 亚洲国产欧美在线一区| 精品99又大又爽又粗少妇毛片| 亚洲精品日韩在线中文字幕 | 亚洲在久久综合| 五月玫瑰六月丁香| 国产精品久久电影中文字幕| 国产极品天堂在线| 亚洲av免费高清在线观看| 国语自产精品视频在线第100页| 乱人视频在线观看| 可以在线观看毛片的网站| 国产私拍福利视频在线观看| 国产精品蜜桃在线观看 | 亚洲av免费在线观看| 高清毛片免费看| 久久99蜜桃精品久久| 人妻少妇偷人精品九色| 日本成人三级电影网站| 亚洲国产精品成人综合色| 亚洲无线观看免费| 国产黄a三级三级三级人| 国产大屁股一区二区在线视频| 亚洲精品成人久久久久久| 免费在线观看成人毛片| 边亲边吃奶的免费视频| 一本久久精品| 三级经典国产精品| 免费观看人在逋| 成人一区二区视频在线观看| 国产蜜桃级精品一区二区三区| 非洲黑人性xxxx精品又粗又长| 精品久久久久久久人妻蜜臀av| 免费不卡的大黄色大毛片视频在线观看 | 国产亚洲精品av在线| 哪里可以看免费的av片| 国产日韩欧美在线精品| 成人三级黄色视频| 成年免费大片在线观看| 午夜福利在线在线| 午夜福利视频1000在线观看| 国产高清不卡午夜福利| 美女内射精品一级片tv| 亚洲色图av天堂| 成人漫画全彩无遮挡| 成人综合一区亚洲| 色吧在线观看| 中文在线观看免费www的网站| 村上凉子中文字幕在线| 日韩一区二区三区影片| 中文字幕av成人在线电影| 久99久视频精品免费| 91久久精品电影网| 最后的刺客免费高清国语| 国产精品一区二区在线观看99 | 中文在线观看免费www的网站| 婷婷六月久久综合丁香| 特大巨黑吊av在线直播| 最新中文字幕久久久久| 一卡2卡三卡四卡精品乱码亚洲| 日韩中字成人| 亚洲欧美精品综合久久99| 亚洲av中文av极速乱| 51国产日韩欧美| 欧美一区二区精品小视频在线| 欧美色视频一区免费| av在线蜜桃| 日韩亚洲欧美综合| 人妻夜夜爽99麻豆av| 麻豆成人午夜福利视频| 一区二区三区高清视频在线| 久久久久网色| 国产免费一级a男人的天堂| 一本久久精品| 高清午夜精品一区二区三区 | 亚洲经典国产精华液单| 最近中文字幕高清免费大全6| 精品少妇黑人巨大在线播放 | а√天堂www在线а√下载| 搞女人的毛片| 一区二区三区四区激情视频 | 欧美一区二区亚洲| 中国美白少妇内射xxxbb| 日韩强制内射视频| kizo精华| 久久综合国产亚洲精品| 伊人久久精品亚洲午夜| 一夜夜www| 国内精品一区二区在线观看| 日日摸夜夜添夜夜爱| 国产毛片a区久久久久| 一区二区三区免费毛片| 小说图片视频综合网站| 在线观看午夜福利视频| 你懂的网址亚洲精品在线观看 | 国产成人影院久久av| 国产精品久久久久久久电影| 久久九九热精品免费| 91久久精品电影网| 国产成人精品久久久久久| 亚洲精品影视一区二区三区av| 黄色日韩在线| 中国美女看黄片| 国产一级毛片七仙女欲春2| 禁无遮挡网站| 欧美三级亚洲精品| 亚洲av免费高清在线观看| 12—13女人毛片做爰片一| 国产成人aa在线观看| 国产三级中文精品| 欧美丝袜亚洲另类| 少妇人妻一区二区三区视频| 美女高潮的动态| 美女大奶头视频| 欧美极品一区二区三区四区| 男女做爰动态图高潮gif福利片| 亚洲国产精品国产精品| 亚洲18禁久久av| 日韩一区二区视频免费看| 亚洲精品国产成人久久av| 内射极品少妇av片p| 亚洲婷婷狠狠爱综合网| 美女脱内裤让男人舔精品视频 | 亚洲不卡免费看| 亚洲一级一片aⅴ在线观看| 国产成人精品一,二区 | 一级毛片aaaaaa免费看小|