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

    基于經(jīng)驗正交函數(shù)的剖面重構(gòu)及其物理意義分析

    2022-04-25 08:37:04鐘瑋琪
    海洋技術(shù)學報 2022年1期
    關(guān)鍵詞:內(nèi)波聲速二階

    鐘瑋琪,屈 科,梁 羿

    (廣東海洋大學電子與信息工程學院,廣東 湛江 524088)

    聲速剖面是重要的海洋水聲環(huán)境參數(shù)之一,是評估海域穩(wěn)定性的基礎(chǔ)指標,也是海洋聲場分析的必要先驗信息。獲取精確的聲速剖面,對海洋環(huán)境監(jiān)測、聲源定位和軍事方面的應(yīng)用都具有重大意義。聲速剖面是有關(guān)深度和聲速的函數(shù),也可以用時間或空間隨深度變化的矩陣形式來描述,聲速剖面的矩陣越大,精度越高。然而,使用深度和聲速的數(shù)值關(guān)系模擬實際聲速剖面需要大量的參數(shù),參數(shù)增加的同時也增大了逆問題求解的難度,因此需要采用一定的降維技術(shù)盡可能使用較少的模態(tài)更為簡單、有效地重構(gòu)聲速剖面。經(jīng)驗正交函數(shù)方法(Empirical Orthogonal Function,EOF) 就被證實是十分有效的途徑之一。為了降低聲速剖面在水聲學逆問題中未知參數(shù)的維度,通常將海洋聲速剖面表示為恒定的背景剖面與擾動項的疊加,通過經(jīng)驗正交函數(shù)方法提取樣本擾動的主成分,然后可以通過若干階EOF 和投影系數(shù)描述聲速剖面。

    經(jīng)驗正交函數(shù)分析方法也稱為特征向量分析或主成分分析,是分析矩陣數(shù)據(jù)結(jié)構(gòu)特性、提取主要數(shù)據(jù)特征的一種方法,能夠?qū)㈦S時間變化的變量場分解為僅依賴時間變化的時間函數(shù)和不隨時間變化的空間函數(shù)。DAVIS R E[1]證明了EOF 是描述聲速剖面最有效的基函數(shù)。沈遠海等[2-3]、張鎮(zhèn)邁等[4]研究表明在淺海使用經(jīng)驗正交函數(shù)描述和重構(gòu)聲速剖面的可行性和有效性。在實際中聲速剖面難以使用連續(xù)函數(shù)表示,傳統(tǒng)聲速剖面重構(gòu)方法對樣本數(shù)量有一定要求,因此,通過正交分解一定數(shù)量聲速剖面樣本序列的時間向量和空間向量,可以有效地描述聲速剖面,通常情況僅需要前幾階就可以獲得足夠的精度[5-7],并且此方法可以避免描述聲速垂直結(jié)構(gòu)多參數(shù)的問題,在反演海水聲速剖面算法執(zhí)行過程中具有較高的時效性[8]。以上研究都表明了即便是在淺海聲速剖面變化劇烈的情況下,采用經(jīng)驗正交函數(shù)表示聲速剖面仍然具有較高的實用性。孫文舟等[9]解析了聲速剖面EOF 第一模態(tài)時間系數(shù)和空間函數(shù)的變化規(guī)律,分析了引起聲速剖面變化的主要因素;蘇林等[10]分析了經(jīng)驗正交函數(shù)系數(shù)分布特點,實現(xiàn)內(nèi)波環(huán)境下對聲速剖面的實時預(yù)測。從上述研究工作可以總結(jié)出,經(jīng)驗正交函數(shù)可以有效地應(yīng)用于聲速剖面時空變化分析[11]、聲速場模型構(gòu)建[12-13]和海洋活動監(jiān)測[14-15]等領(lǐng)域。

    EOF 方法作為一種數(shù)值統(tǒng)計方法,可將復(fù)雜的參考模型方程投影到少數(shù)主成分上來降低問題的維度。由于EOF 模態(tài)是樣本矩陣的主成分,并沒有物理表達式,其物理意義不直觀,難以從中獲取海洋動力學信息[16-18]。為了提升對EOF 方法反演結(jié)果的解析,本文采用經(jīng)驗正交分解法,對2001 年中美聯(lián)合實驗(The Asian Seas International Experiment,ASIAEX)在南海北部測量的溫度鏈數(shù)據(jù)進行處理,對第一、二階投影系數(shù)的物理意義進行發(fā)掘分析。結(jié)果表明前兩階主要模態(tài)與海洋動力學參數(shù)密切相關(guān),通過第一、二階投影系數(shù)的變化可以對海洋動力活動進行有效監(jiān)測。

    1 基于EOF 的聲速重構(gòu)

    本文所采用的數(shù)據(jù)選自2001 年ASIAEX 于南海北部進行聲學實驗中的溫度鏈數(shù)據(jù)。ASIAEX 是中國、美國、俄羅斯、日本、韓國、新加坡之間的綜合科學合作,其主要實驗在2001 年4—8 月進行。實驗包括兩項內(nèi)容,第一個是在東海的海底相互作用實驗,第二個是在南海的體積相互作用實驗。本次數(shù)據(jù)取自2001 年5 月3 日至2001 年5 月17 日的南海實驗,該實驗在浮標投放點(21°55.285′N,117°35.088′E)進行采樣時間為1 min 的測量。浮標的每根繩子放置11 個溫度傳感器,傳感器之間的間隔約為10 m。溫度剖面的深度范圍為21.7~135.1 m。圖1 顯示了72 h(5 月7 日00∶00 至5 月10 日00∶00)溫度剖面的變化,可以觀察到頻繁強烈的線性和非線性內(nèi)波活動,觀測到的內(nèi)孤立波振幅約為70 m。

    圖1 5 月7 日00 ∶00 至5 月10 日00 ∶00 溫度數(shù)據(jù)

    由于影響淺海聲速的主要因素是淺海海水溫度,因此聲速的計算采用Chen-Millero-Li 聲速經(jīng)驗公式[19]。由于實際采樣過程中存在實驗條件的限制,導(dǎo)致每條聲速剖面的采樣并不能完全保持在等深度處,所以利用三次樣條插值法[20]將聲速內(nèi)插到垂直標準層后得到聲速剖面。

    假設(shè)已獲得N 條聲速剖面,每條剖面有M 層垂直標準層,聲速剖面矩陣為C。計算聲速剖面矩陣C 中每層的平均聲速值,聲速剖面矩陣還可以表示成平均聲速剖面與聲速擾動X 的和。通常先將聲速剖面矩陣進行預(yù)處理,即標準化處理或轉(zhuǎn)化為聲速剖面的殘差形式。若對聲速剖面矩陣C 進行標準化處理,那么R 被稱為相關(guān)系數(shù)矩陣;若將C 進行聲速剖面殘差處理,那么R 被稱為協(xié)方差矩陣。然后計算X 與其轉(zhuǎn)置矩陣的交叉積,則X 的協(xié)方差矩陣如下。

    對協(xié)方差矩陣求解特征值向量,將R 分解如下。

    式中,E 是R 的特征值矩陣,將E 按從大到小的順序進行排列。每個非0 的特征值都對應(yīng)一個特征向量,得到F= [f1f2… fn],F(xiàn) 是特征值對應(yīng)的特征向量矩陣,即EOF 的空間函數(shù)。取前q 個EOF 系數(shù),則聲速剖面可以表示如下。

    1.1 重構(gòu)結(jié)果及誤差分析

    對獲得的所有樣本組成的聲速剖面矩陣進行經(jīng)驗正交分解,可得出各階EOF 函數(shù)和方差貢獻率。圖2 為分解得到的前五階歸一化EOF 函數(shù)。表1中第一階EOF 重構(gòu)聲速剖面的方差貢獻率為70.06%,前三階EOF 的方差累計貢獻率達到93.46%,前五階可以達到97.45%。在南海北部海內(nèi)波波動劇烈的條件下,前三階累計貢獻率也已經(jīng)達到96%以上[21],足夠較好地重構(gòu)聲速剖面。因此,從EOF 重構(gòu)聲速剖面方差貢獻率的角度考慮,本文使用前五階足夠反映聲速剖面的主要特征。

    圖2 前五階歸一化EOF

    表1 EOF 的方差貢獻率及累計方差貢獻率

    僅僅使用前幾階擬合聲速剖面和實際所測聲速剖面必然存在誤差,因此,為分析重構(gòu)剖面的準確性,將重構(gòu)的均方根誤差定義如下。

    式中,Czi為實測聲速;Ci為擬合聲速;M 為深度上的采樣點個數(shù)。

    式中,RMSEi為單個剖面均方根誤差;N 為剖面?zhèn)€數(shù)。

    圖3 為重構(gòu)聲速剖面的均方根誤差,清晰地顯示了應(yīng)用EOF 重構(gòu)聲速剖面的均方根誤差均在2.50 m/s 以下,誤差的最大值是2.19 m/s;在統(tǒng)計結(jié)果后,平均均方根誤差為0.83 m/s,42%的重構(gòu)均方根誤差在平均均方根誤差以下,即聲速剖面的重構(gòu)可以達到較好的精度。

    圖3 不同剖面序號重構(gòu)效果的均方根誤差

    海水擾動主要是由內(nèi)波、湍流等動力活動所引起的,對剖面有一定影響。如圖4 所示,選取3 種特殊情況中的3 個具有代表性的剖面進行誤差分析。當海水具有較大擾動時,海水向上運動的重構(gòu)剖面的均方根誤差為0.29 m/s,其余兩個剖面誤差略大于平均誤差。選取的3 個剖面的均方根誤差均接近平均均方根誤差,表明了在海水有較小和較大擾動的情況下,EOF 重構(gòu)聲速剖面仍有較好的效果。

    圖4 振幅最小、振幅向下最大和振幅向上最大的實測剖面(紅色線)與重構(gòu)剖面(藍色線)效果圖

    2 EOF 物理意義的分析

    前二階EOF 模態(tài)的累計方差貢獻率是86.11%,其中第一階模態(tài)占70.06%,在所有模態(tài)中所占權(quán)重最高,在重構(gòu)聲速剖面中起主導(dǎo)作用。從圖5 可見,第一階模態(tài)幅值最大值是在水深為69 m 處,在第一階模態(tài)幅值最大值時,第二階模態(tài)幅值為零,且第一階模態(tài)只有1 個極值,符號都為正值,此時表明了水質(zhì)子運動方向相同。最大振幅出現(xiàn)的位置,在海水躍層中代表著等聲速線,且等聲速線都往同一個方向振蕩。第二階模態(tài)有兩個極值,且在第一階模態(tài)幅值最大水深處附近改變符號,可以解釋為等聲速線在躍層上方和下方運動方向相反,代表了躍層內(nèi)聲速梯度的變化。

    圖5 第一、二階模態(tài)歸一化幅值

    將歸一化后的值A(chǔ) 定義如下。

    式中,an為投影系數(shù);amin為投影系數(shù)的最小值;amax為投影系數(shù)的最大值。

    在忽略第二階以上模態(tài)影響的情況下,第一階投影系數(shù)可近似描述重構(gòu)剖面中1 522 m/s 的等聲速線。根據(jù)物態(tài)方程,鹽度穩(wěn)定時可認為聲速和溫度變化成正比,等聲速線和等溫度線近似一致。如圖6 所示,將第一階投影系數(shù)和等聲速線的振幅進行歸一化,可以注意到第一階投影系數(shù)與一階模態(tài)幅值最大深度等聲速線的變化趨勢相似,并且相關(guān)系數(shù)達到0.96,這說明了第一階投影系數(shù)與該等聲速線具有很強的相關(guān)性。第一階投影系數(shù)與海水水層溫度的變化趨勢相近,因此還可以解釋為,第一階投影系數(shù)可以近似反映出海水水層的變化趨勢。

    圖6 歸一化第一階投影系數(shù)及等聲速線變化

    第二階模態(tài)具有兩個極值,通過求解得到兩個極值(39~95 m)的聲速梯度,將第二階投影系數(shù)和聲速梯度的振幅進行歸一化處理,如圖7 所示。將聲速梯度與第二階投影系數(shù)進行對比,可以明顯觀察到兩條曲線非常相似,變化規(guī)律基本一致。兩條曲線相關(guān)系數(shù)高達0.99,具有十分強的相關(guān)性。因此,第二階投影系數(shù)可以描述二階模態(tài)兩個極值間的聲速梯度,從第二階投影系數(shù)可以得到等聲速線變化的方向和聲速隨深度變化的快慢。由于聲速梯度影響聲線折射的方向,在淺海中導(dǎo)致傳播損失變大,影響水下探測,使得聲吶系統(tǒng)探測距離變短。

    圖7 歸一化第二階投影系數(shù)及聲速梯度變化

    由圖8 可以看到內(nèi)波中聲速隨時間的變化及雙振蕩結(jié)構(gòu)。如圖8(a)所示,該時段出現(xiàn)了5 個顯著連續(xù)的內(nèi)孤立波,并且在出現(xiàn)內(nèi)孤立波的情況下,該時段還產(chǎn)生了類似正余弦變化、符號相反的雙振蕩結(jié)構(gòu),如圖8(b)所示。將相鄰兩個點的振幅除以采樣時間得到投影系數(shù)的時間變化率,可以觀察到雙振蕩結(jié)構(gòu)的兩條曲線分別對應(yīng)著第一階投影系數(shù)與第二階投影系數(shù)的時間變化率,在內(nèi)波活動下第一階投影系數(shù)的時間變化率振蕩峰均略高于第二階投影系數(shù),且隨著時間變化,內(nèi)波振幅深度減小,雙振蕩結(jié)構(gòu)的振幅也隨之減小。

    圖8 內(nèi)波與雙震蕩圖

    海洋內(nèi)波是由海水溫差、鹽差造成密度穩(wěn)定分層現(xiàn)象產(chǎn)生的水下波浪,而內(nèi)波通常發(fā)生在水體的內(nèi)部,在海水表面難以觀察。形成海洋內(nèi)波需要兩個條件,第一是穩(wěn)定的密度分層,第二是擾動源。在海脊、大陸架邊緣等特殊地形條件下,海洋密度躍層中內(nèi)潮波之間的非線性相互作用會產(chǎn)生大振幅的內(nèi)孤立波。由于黑潮入侵和季風影響海洋環(huán)境,南海北部是內(nèi)孤立波活動的多發(fā)海域。內(nèi)孤立波在傳播過程中可以形成海水強烈的輻聚和突發(fā)性的波致流,嚴重威脅海洋工程活動與海洋石油工作。內(nèi)波活動還會影響聲速剖面,進而影響水下目標的探測與定位以及水聲通信。

    在海水深度較淺的情況下,溫度是聲速的主要影響因素,并且溫度和密度相互影響。EOF 模態(tài)在不同深度上的幅值代表深度擾動的大小,在有理想流體媒質(zhì)中有聲擾動的條件下,物態(tài)方程可以描述聲場中聲速與密度的變化之間的關(guān)系。

    式中,c0是聲速;P 是壓強;ρ 是密度,物態(tài)方程反映出海水受到聲擾動的壓縮特性。在理想流體媒質(zhì)中,壓強和密度的變化具有相同的方向。

    密度擾動是引起內(nèi)波的原因之一,假設(shè)海水為理想流體媒質(zhì),由于密度擾動層的密度分布是不穩(wěn)定的,擾動層中的海水密度比周圍的海水密度高,高密度海水必須向低密度海水混合,dρ <0,則dP<0,海水膨脹。此時水體朝下方運動,躍層上方幅度大于下方。水體運動速度與一階模態(tài)近似,對應(yīng)等聲速線向下運動,一階投影系數(shù)隨時間正向增大。躍層內(nèi)梯度變大,對應(yīng)二階投影系數(shù)隨時間負向增大。

    在高低密度海水混合的過程中,一階投影系數(shù)變化率向上增大到峰值則開始減小,二階投影系數(shù)變化率則是向下增大到最大后開始減小,減小到某一時刻的變化率都為0 時,表明了海水混合的速率是由快到慢地恢復(fù)到一個密度平衡的狀態(tài)。

    由于慣性的作用,平衡密度會出現(xiàn)繼續(xù)向下混合的現(xiàn)象。此時水體朝上運動,等聲速線向上運動,對應(yīng)一階投影系數(shù)隨時間負向増大。躍層內(nèi)梯度變小,二階投影系數(shù)隨時間正向增大。

    當海水停止混合時,海水密度已經(jīng)低于平衡密度,周圍高密度海水又開始向低密度海水混合,dρ>0,則dP >0,海水壓縮。與一階模態(tài)對應(yīng)的等聲速線依舊向上運動,躍層內(nèi)梯度繼續(xù)變小。一階投影系數(shù)變化率則向下增大到最大后開始減小,二階投影系數(shù)變化率向上增大到峰值后開始減小,直至某一時刻變化率再次都為0 時,海水重新恢復(fù)到平衡密度,則開始新的周期,形成內(nèi)波?;谏鲜龇治?,內(nèi)波的雙振蕩實質(zhì)為密度層擾動的一個周期。

    3 結(jié) 論

    從聲速剖面的重構(gòu)效果而言,使用EOF 方法對聲速剖面進行分解和重構(gòu)固然有效且重構(gòu)精度較高。然而,作為一種數(shù)值方法,EOF 并沒有實際的物理意義。本文應(yīng)用EOF 方法對聲速剖面進行重構(gòu),對第一、二階投影系數(shù)的物理意義進行探討和分析,得出EOF 系數(shù)擾動與海洋變化規(guī)律結(jié)果如下。

    (1) 一階EOF 模態(tài)與受溫度影響的聲速變化分布相關(guān),第一階投影系數(shù)變化趨勢與等聲速線高度相關(guān),實驗相關(guān)系數(shù)為0.96。

    (2) 二階EOF 模態(tài)主要控制二階躍層梯度,可以使用第二階投影系數(shù)描述二階模態(tài)兩個極值間深度的聲速梯段。

    (3) 內(nèi)孤立波存在時,一、二階EOF 模態(tài)具有相關(guān)性。在溫度和密度的影響下,等聲速線和梯度聯(lián)合擾動的過程可以用前兩階隨時間變化的投影系數(shù)表示成雙振蕩結(jié)構(gòu)。因此可以通過雙振蕩結(jié)構(gòu)作為內(nèi)孤立波的識別特征,進行監(jiān)測和預(yù)警。

    猜你喜歡
    內(nèi)波聲速二階
    孤立內(nèi)波對過渡海域聲場干涉結(jié)構(gòu)的影響分析
    海洋通報(2022年6期)2023-01-07 02:49:06
    內(nèi)波與死水,連潛艇都怕的海浪
    電腦報(2021年20期)2021-08-10 02:28:41
    基于MODIS 遙感影像的安達曼海內(nèi)波特征參數(shù)分布及生成周期研究
    海洋學報(2020年9期)2020-10-09 08:50:08
    一類二階迭代泛函微分方程的周期解
    一類二階中立隨機偏微分方程的吸引集和擬不變集
    二階線性微分方程的解法
    一類二階中立隨機偏微分方程的吸引集和擬不變集
    聲速是如何測定的
    跨聲速風洞全模顫振試驗技術(shù)
    航空學報(2015年4期)2015-05-07 06:43:32
    機翼跨聲速抖振研究進展
    航空學報(2015年4期)2015-05-07 06:43:28
    亚洲av日韩精品久久久久久密| 纵有疾风起免费观看全集完整版| 在线观看一区二区三区激情| 中国国产av一级| 久久精品亚洲av国产电影网| 视频区欧美日本亚洲| 少妇被粗大的猛进出69影院| 下体分泌物呈黄色| av电影中文网址| 9色porny在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩亚洲综合一区二区三区_| 久久毛片免费看一区二区三区| 国产深夜福利视频在线观看| 欧美激情久久久久久爽电影 | 日韩欧美国产一区二区入口| 飞空精品影院首页| av免费在线观看网站| 岛国在线观看网站| 欧美激情久久久久久爽电影 | 99国产精品一区二区蜜桃av | 久久精品国产综合久久久| 窝窝影院91人妻| www.熟女人妻精品国产| 成人三级做爰电影| 欧美中文综合在线视频| 国产亚洲午夜精品一区二区久久| 精品国产一区二区久久| 欧美黄色淫秽网站| 美女主播在线视频| 亚洲精品中文字幕一二三四区 | 亚洲九九香蕉| 视频区欧美日本亚洲| 美女高潮到喷水免费观看| 日本精品一区二区三区蜜桃| kizo精华| 1024香蕉在线观看| 久久毛片免费看一区二区三区| 亚洲精品一区蜜桃| 人人妻,人人澡人人爽秒播| 丝瓜视频免费看黄片| 悠悠久久av| 午夜视频精品福利| 少妇粗大呻吟视频| netflix在线观看网站| 亚洲国产av影院在线观看| www.av在线官网国产| 国产日韩欧美视频二区| 午夜久久久在线观看| 成年美女黄网站色视频大全免费| 国产一级毛片在线| 淫妇啪啪啪对白视频 | 亚洲欧美精品自产自拍| www.自偷自拍.com| 我的亚洲天堂| 欧美日韩成人在线一区二区| 日韩欧美一区二区三区在线观看 | 国产精品一区二区精品视频观看| 国产成+人综合+亚洲专区| 亚洲av成人一区二区三| 热re99久久国产66热| 丁香六月欧美| 青青草视频在线视频观看| 黄片播放在线免费| 熟女少妇亚洲综合色aaa.| 亚洲第一欧美日韩一区二区三区 | 久久久久精品人妻al黑| 国产欧美日韩精品亚洲av| 18禁国产床啪视频网站| 18禁国产床啪视频网站| kizo精华| 免费在线观看完整版高清| 精品乱码久久久久久99久播| 宅男免费午夜| kizo精华| 国产一区二区三区综合在线观看| 成人亚洲精品一区在线观看| 久久女婷五月综合色啪小说| 日韩视频一区二区在线观看| 国产有黄有色有爽视频| 精品国产乱码久久久久久男人| 一个人免费在线观看的高清视频 | 最新的欧美精品一区二区| 亚洲一码二码三码区别大吗| 亚洲av欧美aⅴ国产| 中亚洲国语对白在线视频| 亚洲欧洲精品一区二区精品久久久| 搡老乐熟女国产| 亚洲精品av麻豆狂野| 黄色a级毛片大全视频| 窝窝影院91人妻| 国产成+人综合+亚洲专区| 老司机影院成人| 久久人人爽人人片av| 另类精品久久| 亚洲av成人不卡在线观看播放网 | 国产av又大| 久久精品久久久久久噜噜老黄| 久久精品久久久久久噜噜老黄| 99热网站在线观看| 大片电影免费在线观看免费| 久久国产精品人妻蜜桃| 亚洲av成人一区二区三| 在线观看免费高清a一片| 亚洲欧美清纯卡通| 久热爱精品视频在线9| 国产一区二区 视频在线| 日本av免费视频播放| 久久女婷五月综合色啪小说| av欧美777| 精品国产超薄肉色丝袜足j| 91精品国产国语对白视频| e午夜精品久久久久久久| 亚洲国产精品一区二区三区在线| 亚洲精品一区蜜桃| 中国国产av一级| 国产欧美日韩精品亚洲av| 男女下面插进去视频免费观看| 最近中文字幕2019免费版| 久久精品aⅴ一区二区三区四区| 天天躁狠狠躁夜夜躁狠狠躁| 男女国产视频网站| 久久国产精品人妻蜜桃| 丝袜美足系列| 在线 av 中文字幕| av免费在线观看网站| 国产成人欧美| 伊人久久大香线蕉亚洲五| 精品一区二区三卡| 精品乱码久久久久久99久播| 日韩熟女老妇一区二区性免费视频| 91av网站免费观看| 男女下面插进去视频免费观看| 国产黄频视频在线观看| 如日韩欧美国产精品一区二区三区| bbb黄色大片| 国产成人啪精品午夜网站| 国产一区二区 视频在线| 50天的宝宝边吃奶边哭怎么回事| av国产精品久久久久影院| 亚洲自偷自拍图片 自拍| 97精品久久久久久久久久精品| 亚洲成人免费电影在线观看| 18禁观看日本| 一个人免费看片子| 人妻一区二区av| 肉色欧美久久久久久久蜜桃| 亚洲五月婷婷丁香| 精品一区二区三区四区五区乱码| 熟女少妇亚洲综合色aaa.| 久久久久久久精品精品| 国产在视频线精品| 亚洲专区国产一区二区| 91字幕亚洲| 精品亚洲成a人片在线观看| 久久国产精品男人的天堂亚洲| kizo精华| 欧美日本中文国产一区发布| 欧美xxⅹ黑人| 国产黄色免费在线视频| 黄色片一级片一级黄色片| 亚洲第一欧美日韩一区二区三区 | 青草久久国产| 国产熟女午夜一区二区三区| 国产精品 欧美亚洲| 黄色怎么调成土黄色| 大片免费播放器 马上看| av不卡在线播放| 免费在线观看影片大全网站| 久久久国产欧美日韩av| 亚洲视频免费观看视频| 久久久国产精品麻豆| 国产免费福利视频在线观看| 伊人亚洲综合成人网| 国产精品秋霞免费鲁丝片| 国产激情久久老熟女| 啦啦啦视频在线资源免费观看| 欧美性长视频在线观看| 99热全是精品| 成人国产一区最新在线观看| 少妇被粗大的猛进出69影院| 国产av精品麻豆| 欧美日韩一级在线毛片| 9191精品国产免费久久| 999精品在线视频| 久久中文看片网| 老鸭窝网址在线观看| 高清av免费在线| 老司机深夜福利视频在线观看 | 亚洲国产av新网站| 男男h啪啪无遮挡| 国产无遮挡羞羞视频在线观看| 久热这里只有精品99| 亚洲av日韩在线播放| 少妇人妻久久综合中文| 窝窝影院91人妻| 不卡一级毛片| 欧美在线黄色| 国产伦人伦偷精品视频| 99久久国产精品久久久| 久久热在线av| 亚洲精品成人av观看孕妇| 久久毛片免费看一区二区三区| 母亲3免费完整高清在线观看| 高清av免费在线| 色老头精品视频在线观看| 久久国产精品男人的天堂亚洲| 狠狠婷婷综合久久久久久88av| e午夜精品久久久久久久| 精品久久久精品久久久| 80岁老熟妇乱子伦牲交| 少妇人妻久久综合中文| 蜜桃国产av成人99| 麻豆av在线久日| 大片免费播放器 马上看| 久热这里只有精品99| 国产欧美日韩综合在线一区二区| 国产亚洲av高清不卡| 91成年电影在线观看| 法律面前人人平等表现在哪些方面 | av不卡在线播放| 日本五十路高清| 久久精品成人免费网站| 老司机亚洲免费影院| 在线精品无人区一区二区三| 欧美中文综合在线视频| 丝瓜视频免费看黄片| 这个男人来自地球电影免费观看| 超碰成人久久| 国产男女超爽视频在线观看| 淫妇啪啪啪对白视频 | 亚洲欧美成人综合另类久久久| 国产精品99久久99久久久不卡| 国产三级黄色录像| 制服人妻中文乱码| 美女中出高潮动态图| 国产欧美日韩综合在线一区二区| 国产亚洲av高清不卡| 叶爱在线成人免费视频播放| 天堂8中文在线网| 在线观看一区二区三区激情| 日本vs欧美在线观看视频| 亚洲精品成人av观看孕妇| 亚洲性夜色夜夜综合| av天堂久久9| 99久久人妻综合| a在线观看视频网站| 精品国产国语对白av| 精品视频人人做人人爽| 成人国语在线视频| 日韩制服丝袜自拍偷拍| 青草久久国产| 国产精品久久久久久精品电影小说| 国产精品.久久久| 国产片内射在线| 日本五十路高清| 成人手机av| 欧美精品啪啪一区二区三区 | 夫妻午夜视频| 免费少妇av软件| 国产极品粉嫩免费观看在线| 超碰97精品在线观看| 亚洲av日韩在线播放| 啦啦啦视频在线资源免费观看| 午夜福利一区二区在线看| 日韩一卡2卡3卡4卡2021年| 亚洲精品在线美女| 麻豆av在线久日| 午夜福利影视在线免费观看| 国产一区二区激情短视频 | 啦啦啦 在线观看视频| 欧美另类亚洲清纯唯美| videos熟女内射| 国产欧美日韩一区二区精品| 97精品久久久久久久久久精品| 免费观看人在逋| 精品国内亚洲2022精品成人 | 97人妻天天添夜夜摸| 日韩三级视频一区二区三区| 亚洲国产中文字幕在线视频| 国产成人系列免费观看| 精品一区二区三区四区五区乱码| 18禁观看日本| 亚洲天堂av无毛| 国产精品99久久99久久久不卡| 国产男女超爽视频在线观看| 国产99久久九九免费精品| 久久精品国产亚洲av香蕉五月 | 精品熟女少妇八av免费久了| 97人妻天天添夜夜摸| 亚洲性夜色夜夜综合| 50天的宝宝边吃奶边哭怎么回事| 久久精品国产亚洲av香蕉五月 | 精品国产乱码久久久久久小说| 侵犯人妻中文字幕一二三四区| 午夜福利在线免费观看网站| 啦啦啦视频在线资源免费观看| 少妇被粗大的猛进出69影院| 免费久久久久久久精品成人欧美视频| 国产成人系列免费观看| 欧美黄色淫秽网站| 国产日韩欧美亚洲二区| 一区二区av电影网| 亚洲精品中文字幕在线视频| 国产黄频视频在线观看| 欧美成人午夜精品| 电影成人av| 丰满少妇做爰视频| 一个人免费在线观看的高清视频 | 精品人妻在线不人妻| 亚洲av男天堂| 亚洲少妇的诱惑av| 精品国产乱码久久久久久小说| 亚洲成人免费电影在线观看| 汤姆久久久久久久影院中文字幕| 一边摸一边抽搐一进一出视频| 日韩视频一区二区在线观看| 黄色毛片三级朝国网站| 亚洲精品粉嫩美女一区| 精品一区二区三区av网在线观看 | 99re6热这里在线精品视频| 欧美日韩亚洲国产一区二区在线观看 | 桃花免费在线播放| 亚洲伊人色综图| 最黄视频免费看| 久久久久久免费高清国产稀缺| a级毛片在线看网站| 国产亚洲精品一区二区www | 精品人妻1区二区| 最黄视频免费看| 欧美日韩黄片免| 久久久久国内视频| 亚洲熟女毛片儿| 亚洲国产欧美一区二区综合| 丝袜脚勾引网站| 亚洲七黄色美女视频| 精品国产国语对白av| 国产精品.久久久| 一级毛片女人18水好多| 午夜福利,免费看| 国产片内射在线| 国产在线观看jvid| 婷婷丁香在线五月| 国产成人系列免费观看| 菩萨蛮人人尽说江南好唐韦庄| 每晚都被弄得嗷嗷叫到高潮| 久久久国产成人免费| 丰满饥渴人妻一区二区三| 一级片免费观看大全| 成年人午夜在线观看视频| 欧美激情高清一区二区三区| 99九九在线精品视频| 黄色视频不卡| 久久久久久久久免费视频了| 国产伦理片在线播放av一区| 青青草视频在线视频观看| 日韩欧美一区二区三区在线观看 | 免费观看人在逋| 男男h啪啪无遮挡| av有码第一页| 9191精品国产免费久久| 精品国产一区二区三区久久久樱花| 丁香六月欧美| 国产成人精品在线电影| 国产99久久九九免费精品| 人妻一区二区av| 亚洲国产精品一区三区| 一二三四社区在线视频社区8| 亚洲av电影在线观看一区二区三区| 欧美变态另类bdsm刘玥| 青春草亚洲视频在线观看| 欧美变态另类bdsm刘玥| 国产高清国产精品国产三级| 人妻一区二区av| 久久综合国产亚洲精品| 久久中文看片网| 中国国产av一级| 精品乱码久久久久久99久播| 黄片播放在线免费| 亚洲成av片中文字幕在线观看| 青春草亚洲视频在线观看| 一区二区日韩欧美中文字幕| 日韩中文字幕视频在线看片| 伊人久久大香线蕉亚洲五| 一区二区三区精品91| 天堂俺去俺来也www色官网| 十八禁人妻一区二区| 亚洲精品成人av观看孕妇| 日韩大码丰满熟妇| 超碰97精品在线观看| 亚洲美女黄色视频免费看| 婷婷色av中文字幕| 99久久国产精品久久久| 999久久久国产精品视频| 精品熟女少妇八av免费久了| 嫩草影视91久久| 99久久精品国产亚洲精品| 97在线人人人人妻| 爱豆传媒免费全集在线观看| 精品视频人人做人人爽| 亚洲国产中文字幕在线视频| 国产又色又爽无遮挡免| 精品一区二区三区av网在线观看 | 国产有黄有色有爽视频| 在线 av 中文字幕| 男女边摸边吃奶| 麻豆乱淫一区二区| 久9热在线精品视频| 欧美 日韩 精品 国产| 午夜精品久久久久久毛片777| 天天躁夜夜躁狠狠躁躁| 国产精品秋霞免费鲁丝片| 99久久人妻综合| 亚洲免费av在线视频| 国产精品麻豆人妻色哟哟久久| 精品久久蜜臀av无| 久久精品熟女亚洲av麻豆精品| www.熟女人妻精品国产| 国产精品二区激情视频| 脱女人内裤的视频| 亚洲性夜色夜夜综合| 亚洲va日本ⅴa欧美va伊人久久 | 我要看黄色一级片免费的| 亚洲国产精品一区二区三区在线| 亚洲精品一卡2卡三卡4卡5卡 | 成年人午夜在线观看视频| 国产精品 国内视频| 精品少妇一区二区三区视频日本电影| 中文欧美无线码| 亚洲精品国产av蜜桃| 一区二区av电影网| 天天躁狠狠躁夜夜躁狠狠躁| √禁漫天堂资源中文www| 免费在线观看完整版高清| 黄色怎么调成土黄色| 女人久久www免费人成看片| 色播在线永久视频| 一个人免费在线观看的高清视频 | av有码第一页| 亚洲欧美清纯卡通| 国产精品 国内视频| 黄网站色视频无遮挡免费观看| 久久香蕉激情| 女人高潮潮喷娇喘18禁视频| 婷婷丁香在线五月| 亚洲国产毛片av蜜桃av| 人妻久久中文字幕网| 国产av一区二区精品久久| 国产精品一区二区免费欧美 | 超碰成人久久| 久久久久国内视频| 国产欧美日韩一区二区三区在线| 午夜老司机福利片| 亚洲第一av免费看| 制服人妻中文乱码| 国产一区二区激情短视频 | 日韩大码丰满熟妇| 午夜精品国产一区二区电影| 中文字幕色久视频| 久久久精品免费免费高清| 国产高清videossex| 亚洲国产日韩一区二区| 日韩免费高清中文字幕av| 69av精品久久久久久 | 精品人妻熟女毛片av久久网站| 黄色毛片三级朝国网站| 少妇裸体淫交视频免费看高清 | 国产有黄有色有爽视频| 黄色视频在线播放观看不卡| 久9热在线精品视频| 亚洲精品一区蜜桃| 午夜久久久在线观看| 99国产精品免费福利视频| 丁香六月欧美| 国产伦人伦偷精品视频| 午夜两性在线视频| 热99re8久久精品国产| videos熟女内射| 午夜久久久在线观看| 亚洲精品美女久久久久99蜜臀| 日日摸夜夜添夜夜添小说| 国内毛片毛片毛片毛片毛片| 一二三四在线观看免费中文在| 成人手机av| 99久久人妻综合| 一级片'在线观看视频| 深夜精品福利| 人妻人人澡人人爽人人| 一级,二级,三级黄色视频| 成人av一区二区三区在线看 | 久久久精品免费免费高清| 人人妻人人添人人爽欧美一区卜| 精品国产一区二区三区久久久樱花| 色老头精品视频在线观看| 成年人黄色毛片网站| 99国产极品粉嫩在线观看| 2018国产大陆天天弄谢| 亚洲国产日韩一区二区| 不卡av一区二区三区| 18在线观看网站| 成年av动漫网址| 国产亚洲欧美精品永久| 久久av网站| 亚洲av片天天在线观看| 欧美精品人与动牲交sv欧美| 老汉色∧v一级毛片| 最新在线观看一区二区三区| 两个人免费观看高清视频| 午夜免费鲁丝| 欧美日韩亚洲综合一区二区三区_| 久久国产精品影院| av福利片在线| 巨乳人妻的诱惑在线观看| 淫妇啪啪啪对白视频 | 日本精品一区二区三区蜜桃| 国产亚洲精品一区二区www | 99国产精品一区二区蜜桃av | 亚洲精品成人av观看孕妇| 18禁黄网站禁片午夜丰满| 亚洲国产看品久久| 一区二区三区四区激情视频| 亚洲精品中文字幕一二三四区 | 日本黄色日本黄色录像| 免费在线观看视频国产中文字幕亚洲 | 免费一级毛片在线播放高清视频 | 多毛熟女@视频| 十八禁网站网址无遮挡| 无限看片的www在线观看| 一区二区三区精品91| 极品少妇高潮喷水抽搐| 久久影院123| 国产av精品麻豆| 国产视频一区二区在线看| 日本vs欧美在线观看视频| 久久久欧美国产精品| 99九九在线精品视频| 亚洲精品日韩在线中文字幕| 人人妻,人人澡人人爽秒播| www.熟女人妻精品国产| 国产精品一区二区在线观看99| 日韩 亚洲 欧美在线| 三级毛片av免费| 精品人妻在线不人妻| 精品高清国产在线一区| av网站免费在线观看视频| 黄色 视频免费看| 一级片'在线观看视频| 91精品三级在线观看| 99国产精品一区二区三区| 如日韩欧美国产精品一区二区三区| 天天影视国产精品| 国产成人av激情在线播放| 久久人妻熟女aⅴ| 国产成人欧美| 91麻豆av在线| 岛国毛片在线播放| 欧美日韩精品网址| 18禁黄网站禁片午夜丰满| 国产精品久久久人人做人人爽| 精品高清国产在线一区| 久久天躁狠狠躁夜夜2o2o| 精品亚洲乱码少妇综合久久| 精品一区在线观看国产| 精品国产一区二区三区久久久樱花| 99热全是精品| 色婷婷av一区二区三区视频| 人妻一区二区av| 热99国产精品久久久久久7| 久9热在线精品视频| 国产淫语在线视频| 亚洲 欧美一区二区三区| 欧美精品一区二区大全| 999久久久国产精品视频| av又黄又爽大尺度在线免费看| 日韩三级视频一区二区三区| 老熟妇乱子伦视频在线观看 | 一边摸一边做爽爽视频免费| 少妇人妻久久综合中文| 1024香蕉在线观看| 午夜福利影视在线免费观看| 成人黄色视频免费在线看| 少妇精品久久久久久久| 精品国内亚洲2022精品成人 | 欧美老熟妇乱子伦牲交| 人妻人人澡人人爽人人| 欧美成狂野欧美在线观看| 久久精品熟女亚洲av麻豆精品| 国产1区2区3区精品| 日韩大码丰满熟妇| 久久久国产欧美日韩av| 久久99热这里只频精品6学生| 亚洲黑人精品在线| 99久久人妻综合| 黄片播放在线免费| 99国产精品一区二区蜜桃av | 婷婷成人精品国产| 午夜激情久久久久久久| 性高湖久久久久久久久免费观看| 国产视频一区二区在线看| 成人国产av品久久久| 精品免费久久久久久久清纯 | 一级黄色大片毛片| 国产区一区二久久| 午夜福利,免费看| 久久久久网色| 日本撒尿小便嘘嘘汇集6| 亚洲成人国产一区在线观看| 老司机影院成人| 人人妻人人澡人人看| 1024香蕉在线观看| 岛国在线观看网站| 热99re8久久精品国产| 亚洲国产精品999| 国产精品久久久久久精品电影小说|