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

    渤海夏季第一斜壓羅斯貝變形半徑的計算與分析

    2015-11-28 11:08:24張曉璐熊學(xué)軍
    海洋科學(xué)進(jìn)展 2015年3期
    關(guān)鍵詞:斜壓浮性渤海

    張曉璐,熊學(xué)軍*

    (1.國家海洋局 第一海洋研究所,山東 青島266061;2.海洋環(huán)境科學(xué)和數(shù)值模擬國家海洋局重點(diǎn)實(shí)驗(yàn)室,山東 青島266061)

    渤海,是一個深入中國大陸的半封閉淺海,北、西、南三面被陸地包圍,僅通過東面的渤海海峽與黃海溝通。注入渤海的沿岸河流眾多,北有遼河、大遼河、大凌河,西有海河、灤河,南有黃河、小清河,因此形成了渤海獨(dú)特的水文特征及動力結(jié)構(gòu)。相較于冬季,渤海夏季的風(fēng)速小,主要為熱鹽流。近年來,大量的數(shù)值模擬工作及調(diào)查研究表明,夏季渤海多存在一些局地范圍的環(huán)流,如渤海中部的順時針環(huán)流,長興島以西的遼東灣海域的逆時針環(huán)流,以及在渤海南部的萊州灣存在弱的順時針環(huán)流等[1]。

    在海洋環(huán)流的理論研究和數(shù)值模擬中,羅斯貝變形半徑是一個極其重要的水平尺度,它描述了地球旋轉(zhuǎn)對流體運(yùn)動的影響。這個概念,最早是由Carl-Gustav Rossby于1938年在研究地轉(zhuǎn)調(diào)整中提出的。若海面高度變化在調(diào)整過程中的尺度與羅斯貝變形半徑相比為小量,則壓強(qiáng)梯度非常大,重力支配運(yùn)動的狀態(tài);但是,當(dāng)海面高度的變化發(fā)展到具有羅斯貝變形半徑尺度的范圍時,科氏力項(xiàng)變得同壓強(qiáng)梯度項(xiàng)同等重要[2]。簡言之,即在羅斯貝變形半徑這個距離尺度上,科氏力使自由面變形的趨勢與壓強(qiáng)梯度力使自由面復(fù)原的趨勢相平衡。在海洋中,斜壓羅斯貝變形半徑是一種常與邊界現(xiàn)象(諸如邊界流、鋒面)以及渦旋有聯(lián)系的自然尺度,是設(shè)置數(shù)值模式參數(shù)的重要參考。有研究表明,在模擬中尺度運(yùn)動時,網(wǎng)格大小應(yīng)當(dāng)是羅斯貝變形半徑的一半,或者是其三分之一[3]。

    國外對斜壓羅斯貝變形半徑的研究已有不少。Emery等[4]在1984年用氣候態(tài)的溫鹽剖面數(shù)據(jù)計算了北太平洋和北大西洋的斜壓羅斯貝變形半徑;Houry等[5]在1987年用NODC的溫鹽數(shù)據(jù)計算了南大西洋的斜壓羅斯貝變形半徑;特別地,Chelton等[6]在1996年利用NODC的1°×1°的氣候態(tài)的溫鹽數(shù)據(jù)計算了全球大洋的第一斜壓羅斯貝變形半徑以及對應(yīng)的全球范圍的第一斜壓重力波相速度,對后來學(xué)者的研究工作起了重要的啟示和參考作用。除了大洋中,在內(nèi)海也有針對斜壓羅斯貝變形半徑的專門研究,如Fennel等[7]、Alenius等[8]以及Osinsik等[9]曾分別利用航測溫鹽剖面數(shù)據(jù)計算了波羅的海、芬蘭灣的斜壓羅斯貝變形半徑。但在國內(nèi),有關(guān)斜壓羅斯貝變形半徑的研究工作不多,甘子鈞和蔡樹群[9],蔡樹群等[10],曾用Levitus和NODC標(biāo)準(zhǔn)層上的氣候態(tài)的溫鹽資料計算過南海斜壓羅斯貝變形半徑的地理及季節(jié)變化以及相應(yīng)的第一斜壓重力波相速度。

    在渤海進(jìn)行第一斜壓羅斯貝變形半徑的計算和研究尚未發(fā)現(xiàn)。本文在觀測數(shù)據(jù)基礎(chǔ)上,分別利用WKB近似和數(shù)值解法,對渤海夏季第一斜壓羅斯貝變形半徑進(jìn)行了計算和分析,希望有助于渤海環(huán)流、中尺度渦旋,以及與此相關(guān)的其他海洋現(xiàn)象的研究。

    1 資料

    本文所用資料為渤海2006年夏季的溫度和鹽度大面觀測資料,觀測儀器為日本ALEC公司的AAQ1183系列的CTD,觀測時間為2006-07-18—08-12。為方便對下文計算結(jié)果分析對比,將整個渤海劃分為6個區(qū)域,分別為渤海東部海區(qū)、遼東灣海域、北戴河附近海域、天津附近海域、黃河口三角洲海域以及萊州灣海域。觀測站位分布,渤海區(qū)劃及渤海水深如圖1所示。

    2 計算方法

    羅斯貝變形半徑在數(shù)學(xué)表達(dá)上,可以理解為:在地球旋轉(zhuǎn)特征周期這一時間尺度上,長重力波傳播的特征距離[11]。在緯度大于5°的海域,羅斯貝變形半徑可以寫為:

    式中,f=2Ωsin(θ)是科氏參數(shù);Ω為地球自轉(zhuǎn)角速度;θ為地理緯度。當(dāng)n=0時,c0為海洋中正壓波的相速度,對應(yīng)的R0為正壓羅斯貝變形半徑;當(dāng)n=1時,c1為第一斜壓重力波相速度,對應(yīng)的R1為第一斜壓羅斯貝變形半徑。第一斜壓羅斯貝變形半徑在海洋環(huán)流理論方面起重要作用,是中尺度渦、沿岸射流和赤道流等的水平尺度[10],故本文主要研究第一斜壓重力波相速度以及對應(yīng)的第一斜壓羅斯貝變形半徑。通常如不特別說明,斜壓羅斯貝變形半徑均指第一斜壓羅斯貝變形半徑。

    在不考慮底摩擦和風(fēng)應(yīng)力,引入平底、海面剛蓋近似,及Boussinesq近似的假設(shè)下,穩(wěn)定層化海洋水平大尺度波動的垂直結(jié)構(gòu)滿足:

    式中,H為水深;φ(z)為垂向速度分離變量后只與z有關(guān)的垂向函數(shù);N(z)為Brunt-V?is?l?頻率,亦稱浮性頻率,忽略壓縮性的影響,其表達(dá)式為:

    式中,ρ表示密度;g表示重力加速度。

    方程(2)及邊界條件(3)組成了典型的Sturm-Liouville問題,求解羅斯貝變形半徑,即轉(zhuǎn)換成求解Sturm-Liouville本征值的問題[2]。對于連續(xù)層化的海洋,方程(2)對應(yīng)的本征值是一組有序、非負(fù)的數(shù)列:。對應(yīng)的本征函數(shù)為φ1(z),φ2(z),φ3(z)…,下標(biāo)表示模態(tài)數(shù)(mode number)。cn表征第n模態(tài)斜壓重力波相速度,由c1和科氏參數(shù)f便可確定第一斜壓羅斯貝變形半徑R1。

    2.1 WKB近似

    在求解由方程(2)及邊界條件(3)組成的Sturm-Liouville問題中,WKB近似是一種快捷有效的方法[6]。其基本思想是把φ(z)寫成下面的小參數(shù)冪級數(shù)的指數(shù)形式:

    根據(jù)本文的情況,在小參數(shù)冪級數(shù)中引入虛數(shù)單位i。將式(5)代入方程(2)中,令ε的同次冪系數(shù)相等,可得遞推方程組:

    通過遞推方程組求解χ0(z),χ1(z),代入到式(5)中,只保留冪級數(shù)顯著的前兩項(xiàng),即可得出φ(z)的WKB近似解:

    式中,A,B為任意常系數(shù),代入邊界條件(3)得:

    當(dāng)n=1時,通過式(8)和式(9)便可求解出第一斜壓重力波相速度c1和第一斜壓羅斯貝變形半徑R1。

    2.2 數(shù)值解法

    數(shù)值解法是將本征值問題離散化后求解的過程。在相鄰離散深度zk、zk+1上,利用中間密度梯度法(neutral density gradient method)[6],將現(xiàn)場密度絕熱調(diào)整到中間深度ζk=zk+1/2(圖2),然后估計浮性頻率的平方N2(z):

    圖2 數(shù)值解法中在離散深度上求解浮性頻率的示意圖Fig.2 Schematic diagram of vertical discretization in numerical method

    據(jù)Chelton等[6]的建議,為了求羅斯貝變形半徑,如果得到的N2(ζk)為負(fù)值,用相鄰較淺深度ζk-1上的N2(ζk-1)代替;若N2(ζ1)為負(fù)值,令N2(ζ1)為10-8。

    為了得到數(shù)值解,在深度ζk上離散本征值方程。首先,φ(z)在z=ζk(k=1,…,n-1)上進(jìn)行泰勒展開:

    通過式(11)求取φ(ζk-1)、φ(ζk+1),消除φ′(ζk)項(xiàng)后有:

    其中,

    把式(11)代入到方程(2)中,得到n-1個離散方程,移項(xiàng)整理得:

    式中,

    根據(jù)本文的情況,定義邊界條件為:

    至此,將求解Sturm-Liouville本征值問題轉(zhuǎn)換成求n-1階矩陣特征值的問題[6]:

    式中,

    由式(18)可解出n-1個特征值μ1,…,μn-1,對應(yīng)第一個特征值μ1,通過式(16)解出第一斜壓重力波相速度c1和第一斜壓羅斯貝變形半徑R1。

    3 結(jié)果分析

    圖3和圖4分別給出了通過WKB近似和數(shù)值解法求出的渤海夏季第一斜壓重力波相速度和第一斜壓羅斯貝變形半徑的分布。兩種算法下得到的渤海夏季第一斜壓重力波相速度和第一斜壓羅斯貝變形半徑的分布格局相似。斜壓重力波相速度與斜壓羅斯貝變形半徑依賴于浮性頻率的積分深度,即與實(shí)際水深有密切關(guān)系,故它們的分布與渤海水深分布大體一致,等值線呈耳狀分布。在北戴河和長興島附近海域,斜壓重力波相速度與斜壓羅斯貝變形半徑最大,然后向岸逐漸遞減。在渤海中部等值線發(fā)生明顯彎曲,在靠近渤海海峽內(nèi)側(cè)斜壓重力波相速度與斜壓羅斯貝變形半徑出現(xiàn)一個小值,而與之對應(yīng)在渤海中部以西出現(xiàn)一個大值,導(dǎo)致等值線呈NE-SW向傾斜的“S”狀。

    圖3 WKB近似結(jié)果Fig.3 Result of WKB method

    圖4 數(shù)值方法計算結(jié)果Fig 4 Result of numerical method

    表1 WKB近似計算得出的第一斜壓重力波相速度和第一斜壓羅斯貝變形半徑Table 1 The first baroclinic gravity wave phase speed and Rossby deformation radius calculated by WKB method

    表2 數(shù)值方法計算得出的第一斜壓重力波相速度和第一斜壓羅斯貝變形半徑Table 2 The first baroclinic gravity wave phase speed and Rossby deformation radius calculated by numerical method

    在WKB近似下,渤海夏季第一斜壓羅斯貝變形半徑的最大值出現(xiàn)在長興島附近海域,為3.13km,對應(yīng)的斜壓重力波相速度最大為0.28m/s;在萊州灣海域斜壓羅斯貝變形半徑和對應(yīng)的斜壓重力波相速度最大值僅為0.95km和0.08m/s;在遼東灣海域、北戴河附近海域、天津附近海域和黃河口三角洲海域,斜壓羅斯貝變形半徑和斜壓重力波相速度的最大值分別為1.58,2.56,2.02和1.60km,以及0.15,0.24,0.18和0.14m/s。

    在數(shù)值解法下,第一斜壓羅斯貝變形半徑與第一斜壓重力波相速度的值均比在 WKB近似下的值大。最大值出現(xiàn)在北戴河附近海域,為4.72km和0.44m/s;次最大值出現(xiàn)在長興島附近海域,為4.51km和0.42m/s;在萊州灣海域斜壓羅斯貝變形半徑和斜壓重力波相速度的最大值較小,分別為1.76km和0.16 m/s。在遼東灣海域、天津附近海域和黃河口三角洲海域的最大值分別為2.97,2.37,2.55km和0.28,0.22,0.23m/s。兩種算法得出的第一斜壓羅斯貝變形半徑的比較見圖5。

    圖5 兩種算法得出的第一斜壓羅斯貝變形半徑的比較Fig.5 Comparison between the WKB approximation and numerical estimate R1 of the first baroclinic Rossby deformation radius

    從圖5可以看出,渤海夏季第一斜壓羅斯貝變形半徑在WKB近似下的結(jié)果總體低于數(shù)值解法求出的結(jié)果,同樣的情況也出現(xiàn)在其它海區(qū)的相關(guān)研究中。Fennel等[7]和Osinski等[8]在波羅的海的斜壓羅斯貝變形半徑的計算中均表明,WKB近似的結(jié)果偏小于數(shù)值解,F(xiàn)ennel指出偏小10%~30%[7]。這主要是因?yàn)椋琖KB近似適用于浮性頻率垂向變化較弱的情況。由式(9)和式(10)可以看出,在 WKB近似中,第一斜壓重力波的相速度和第一斜壓羅斯貝變形半徑依賴于浮性頻率隨水深的積分,而在實(shí)際求解的過程中,將連續(xù)的海洋分為若干層處理,浮性頻率一般取所在層的平均值,忽略了浮性頻率的變化。故當(dāng)浮性頻率垂向變化較大,WKB近似計算誤差較大。但在深海大洋中計算羅斯貝變形半徑,WKB近似的結(jié)果偏大[6],在此不作詳細(xì)討論。

    圖3b和圖4b中標(biāo)出了斜壓羅斯貝變形半徑出現(xiàn)極值的區(qū)域,其中在北戴河和長興島附近海域的極大值區(qū)有13個站點(diǎn),在渤海中部以西海域的極大值區(qū)有7個站點(diǎn),而在渤海中部以東海域的極小值區(qū)有9個站點(diǎn)。根據(jù)這些站點(diǎn)的溫鹽數(shù)據(jù),計算了極值區(qū)域的浮性頻率的垂向分布,如圖6所示。

    圖6 渤海三極值區(qū)浮性頻率垂向分布圖Fig.6 Profile of Brunt-V?is?l?frequency in the three extreme areas

    由圖6可見,在北戴河和長興島附近海域,以及渤海中部以西海域,浮性頻率的垂向分布較大,導(dǎo)致該海域較強(qiáng)的斜壓性,因此計算得出的第一斜壓羅斯貝變形半徑在該海域出現(xiàn)極大值。而在渤海中部以東海域,浮性頻率的垂向分布整體較小,計算得出的第一斜壓羅斯貝變形半徑在該海域出現(xiàn)極小值。

    圖7繪制了整個渤海夏季浮性頻率的垂向分布??梢娫谒?dbar附近浮性頻率達(dá)到最大,達(dá)0.13s-1,這一深度也是渤海夏季躍層發(fā)生的深度。隨著深度變大,浮性頻率逐漸減小,到水深20dbar的附近,浮性頻率降至約0.04s-1,因此可以看出,渤海夏季浮性頻率垂向上變化較大,故在計算第一斜壓羅斯貝變形半徑的過程中,WKB近似的誤差較大,建議采用斜壓羅斯貝變形半徑的數(shù)值解。

    圖7 整個渤海夏季浮性頻率垂向分布圖Fig.7 Vertical distribution of Brunt-V?is?l?frequency in the whole Bohai Sea

    4 結(jié)論

    本文利用渤海2006年夏季的溫度和鹽度大面觀測資料,運(yùn)用WKB近似和數(shù)值方法計算了渤海夏季第一斜壓羅斯貝變形半徑和第一斜壓重力波相速度。結(jié)果顯示,兩種不同算法得出的斜壓羅斯貝變形半徑和斜壓重力波相速度的分布格局相似,與渤海等深線分布大體一致。數(shù)值解法下,在北戴河和長興島附近海域,斜壓羅斯貝變形半徑較大,最大為4.72km,對應(yīng)的斜壓重力波相速度最大為0.44m/s;在萊州灣海域,斜壓羅斯貝變形半徑較小,最大僅為1.76km,對應(yīng)的斜壓重力波相速度最大為0.16m/s;在渤海中部斜壓羅斯貝變形半徑和斜壓重力波相速度的等值線發(fā)生明顯彎曲,在靠近渤海海峽內(nèi)側(cè)出現(xiàn)一個小值,而與之對應(yīng)在渤海中部以西出現(xiàn)一個大值,導(dǎo)致等值線呈NE-SW向傾斜的“S”狀。在遼東灣海域、天津附近海域和黃河口三角洲海域,斜壓羅斯貝變形半徑和斜壓重力波相速度的最大值分別為2.97,2.37和2.55km,以及0.28,0.22和0.23m/s。

    渤海夏季浮性頻率的垂向變化較大,在水深8dbar附近浮性頻率達(dá)到最大。在WKB近似計算渤海夏季第一斜壓羅斯貝變形半徑的過程中,因其依賴于浮性頻率對水深的積分,得出的結(jié)果誤差較大,其值低于數(shù)值解法求出的結(jié)果。

    (References):

    [1]XIONG X J.The circulation structure and mechanism studies on the China Seas[D].Qingdao:Ocean University of China,2013.熊學(xué)軍.中國近海環(huán)流及其發(fā)生機(jī)制研究[D].青島:中國海洋大學(xué),2013.

    [2]GILL A E.Atmosphere-Ocean dynamics[M].New York:Academic Press,1982.

    [3]ALENIUS P,NEKRASOV A,MYRBERG K.Variability of the baroclinic Rossby radius in the Gulf of Finland[J].Continental Shelf Research,2003,23(6):563-573.

    [4]EMERY W J,LEE W G,MAGAARD L.Geographic and deasonal distributions of Brunt-V?is?l?frequency and Rossby radii in the North Pacific and North Atlantic[J].Journal of Physical Oceanography,1984,14(2):294-317.

    [5]HOURY S,DOMBROWSKY E,DEMEY P,et al.Brunt-V?is?l?frequency and Rossby Radii in the South Atlantic[J].Journal of Physical Oceanography,1987,17(10):1619-1626.

    [6]CHELTON D B,DESZOEKE R A,SCHLAX M G.Geographical variability of the first baroclinic Rossby radius of deformation[J].Journal of Physical Oceanography,1998,28(3):433-460.

    [7]FENNEL W,SEIFERT T,KAYSER B.Rossby radii and phase speeds in the Baltic Sea[J].Continental Shelf Research,1991,11(1):23-36.

    [8]OSINSKI R,RAK D,WALCZOWSKI W,et al.Baroclinic Rossby radius of deformation in the southern Baltic Sea[J].Oceanologia,2010,52(3):417-429.

    [9]GAN Z J,CAI S Q.Geographical and seasonal variability of Rossby radii in South China Sea[J].Journal of Tropical Oceanography,2001,20(1):1-9.甘子鈞,蔡樹群.南海羅斯貝變形半徑的地理及季節(jié)變化[J].熱帶海洋學(xué)報,2001,20(1):1-9.

    [10]CAI S Q,LONG X M,WU R H,et al.Geographical and monthly variability of the first baroclinic Rossby radius of deformation in the South China Sea[J].Journal of Marine Systems,2008,74(1-2):711-720.

    [11]YU Z H,YANG D S,HE H Y,et al.Geophysical fluid dynamics[M].Beijing:China Meteorological Press,1996.余志豪,楊大升,賀海晏,等.地球物理流體動力學(xué)[M].北京:氣象出版社,1996.

    猜你喜歡
    斜壓浮性渤海
    渤海大學(xué)作品精選
    聲屏世界(2022年17期)2022-12-18 01:42:44
    淺析船舶壓載水處理技術(shù)及處理系統(tǒng)選型
    倒扣駁船扳正過程的仿真模擬
    航海(2020年5期)2020-11-02 09:35:16
    Gemini表面活性劑體系下鈣質(zhì)磷礦中白云石的可浮性研究
    “烈焰”,在渤海邊燃燒
    斜壓渦度的變化與臺風(fēng)暴雨的關(guān)系研究
    基于Visual FoxPro的船舶靜水力計算及曲線繪制
    價值工程(2017年3期)2017-05-14 02:11:40
    遼代“斡魯朵”內(nèi)的渤海人移民
    Rossby波的線性穩(wěn)定性
    科技資訊(2016年23期)2016-05-30 20:25:49
    2005—2009年、2011年和2013年南海東北部120°E斷面秋季體積輸運(yùn)的年際變化*
    海洋與湖沼(2016年1期)2016-01-15 03:50:46
    亚洲国产高清在线一区二区三 | 777久久人妻少妇嫩草av网站| 国产高清视频在线播放一区| 亚洲第一青青草原| 一级片免费观看大全| 欧美色欧美亚洲另类二区| 久久热在线av| 亚洲国产高清在线一区二区三 | 亚洲免费av在线视频| 久久国产精品人妻蜜桃| 在线观看www视频免费| 在线永久观看黄色视频| 日韩高清综合在线| 欧美久久黑人一区二区| 999久久久精品免费观看国产| 久热这里只有精品99| 男女之事视频高清在线观看| 婷婷丁香在线五月| 女人爽到高潮嗷嗷叫在线视频| 大型黄色视频在线免费观看| 亚洲无线在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲av中文字字幕乱码综合 | xxxwww97欧美| 亚洲一区二区三区不卡视频| 午夜免费成人在线视频| 亚洲九九香蕉| 国产极品粉嫩免费观看在线| 中文字幕高清在线视频| 神马国产精品三级电影在线观看 | 在线观看免费午夜福利视频| 国产成人精品久久二区二区免费| 婷婷精品国产亚洲av在线| 久久久久久国产a免费观看| 后天国语完整版免费观看| 91大片在线观看| av免费在线观看网站| 亚洲国产看品久久| 在线观看午夜福利视频| 人人澡人人妻人| 欧美成人午夜精品| 99在线人妻在线中文字幕| 给我免费播放毛片高清在线观看| 国产在线精品亚洲第一网站| 悠悠久久av| 日本精品一区二区三区蜜桃| videosex国产| 首页视频小说图片口味搜索| 亚洲一区中文字幕在线| 久久99热这里只有精品18| 大型av网站在线播放| 欧美一级毛片孕妇| 少妇 在线观看| 久久狼人影院| 国产精品久久视频播放| 免费在线观看日本一区| 久久欧美精品欧美久久欧美| 国产伦人伦偷精品视频| 亚洲精品国产一区二区精华液| 成人国产综合亚洲| 亚洲激情在线av| 麻豆av在线久日| 丰满人妻熟妇乱又伦精品不卡| 国产精品98久久久久久宅男小说| 亚洲男人天堂网一区| 午夜福利视频1000在线观看| 国产精品电影一区二区三区| 精品电影一区二区在线| 国内毛片毛片毛片毛片毛片| 午夜福利成人在线免费观看| 欧美日韩精品网址| 中文字幕人妻丝袜一区二区| 久久久精品国产亚洲av高清涩受| 搞女人的毛片| 国产精品久久久人人做人人爽| 免费看十八禁软件| 最新美女视频免费是黄的| 国产亚洲av高清不卡| 宅男免费午夜| 日本一区二区免费在线视频| 熟妇人妻久久中文字幕3abv| 他把我摸到了高潮在线观看| 国产极品粉嫩免费观看在线| 欧美 亚洲 国产 日韩一| 午夜福利18| 精品免费久久久久久久清纯| 精品午夜福利视频在线观看一区| 美女国产高潮福利片在线看| 黄片小视频在线播放| 亚洲五月婷婷丁香| 在线国产一区二区在线| www日本黄色视频网| 欧美zozozo另类| 色av中文字幕| 久久精品91蜜桃| 免费高清视频大片| 免费在线观看亚洲国产| 可以免费在线观看a视频的电影网站| 欧美乱码精品一区二区三区| 久久亚洲精品不卡| 欧美色视频一区免费| 久久国产精品人妻蜜桃| 12—13女人毛片做爰片一| 波多野结衣高清作品| 可以在线观看的亚洲视频| 窝窝影院91人妻| ponron亚洲| or卡值多少钱| 成人国语在线视频| 女同久久另类99精品国产91| videosex国产| √禁漫天堂资源中文www| 精品不卡国产一区二区三区| 亚洲av美国av| 亚洲成a人片在线一区二区| 亚洲va日本ⅴa欧美va伊人久久| 国产男靠女视频免费网站| 搡老妇女老女人老熟妇| 91大片在线观看| 99久久久亚洲精品蜜臀av| videosex国产| 久久伊人香网站| 熟女电影av网| 一个人免费在线观看的高清视频| 亚洲国产精品久久男人天堂| 久久精品人妻少妇| 母亲3免费完整高清在线观看| www.熟女人妻精品国产| 亚洲 国产 在线| 极品教师在线免费播放| 中文字幕最新亚洲高清| av天堂在线播放| 日日干狠狠操夜夜爽| 国产精品 欧美亚洲| 黄网站色视频无遮挡免费观看| 天堂动漫精品| 精品午夜福利视频在线观看一区| 老司机靠b影院| 亚洲狠狠婷婷综合久久图片| 伦理电影免费视频| 午夜福利高清视频| 久久久久久久精品吃奶| 91大片在线观看| 大型黄色视频在线免费观看| 母亲3免费完整高清在线观看| 黄片播放在线免费| 亚洲全国av大片| 久久中文字幕一级| 最近最新中文字幕大全电影3 | 母亲3免费完整高清在线观看| 国产精品一区二区免费欧美| 国产精品影院久久| 大香蕉久久成人网| 欧美乱色亚洲激情| 久久99热这里只有精品18| 99在线视频只有这里精品首页| 亚洲国产欧美网| 欧美在线黄色| 久99久视频精品免费| 成人永久免费在线观看视频| 国产精品二区激情视频| 亚洲成a人片在线一区二区| 成年版毛片免费区| 亚洲av成人av| 99久久精品国产亚洲精品| 欧美日韩精品网址| 亚洲性夜色夜夜综合| 色av中文字幕| 波多野结衣高清无吗| 两个人视频免费观看高清| 一本精品99久久精品77| 亚洲美女黄片视频| 久久精品91无色码中文字幕| 亚洲国产日韩欧美精品在线观看 | 一级毛片女人18水好多| 久久国产乱子伦精品免费另类| 国产精品久久视频播放| 丁香六月欧美| 亚洲国产精品sss在线观看| 俺也久久电影网| 99国产精品99久久久久| 啦啦啦免费观看视频1| 亚洲成av人片免费观看| 日韩精品中文字幕看吧| 亚洲无线在线观看| 日本一本二区三区精品| 亚洲国产精品sss在线观看| a在线观看视频网站| 午夜日韩欧美国产| 亚洲第一青青草原| 国产亚洲欧美精品永久| 日韩欧美国产一区二区入口| 欧美乱妇无乱码| 亚洲午夜理论影院| 久久99热这里只有精品18| 黄色片一级片一级黄色片| 婷婷亚洲欧美| 久久精品影院6| 老司机午夜福利在线观看视频| 中文亚洲av片在线观看爽| videosex国产| 亚洲第一电影网av| 在线视频色国产色| 别揉我奶头~嗯~啊~动态视频| 熟妇人妻久久中文字幕3abv| 女性生殖器流出的白浆| 欧美国产日韩亚洲一区| 亚洲全国av大片| 久久久久久大精品| 国产区一区二久久| 99久久久亚洲精品蜜臀av| 国产欧美日韩精品亚洲av| 亚洲色图 男人天堂 中文字幕| 美国免费a级毛片| 免费电影在线观看免费观看| 亚洲欧美精品综合久久99| 色哟哟哟哟哟哟| 成人午夜高清在线视频 | 制服丝袜大香蕉在线| 久久国产精品人妻蜜桃| 国产真人三级小视频在线观看| tocl精华| 亚洲精品一卡2卡三卡4卡5卡| 琪琪午夜伦伦电影理论片6080| 欧美黄色淫秽网站| ponron亚洲| 一个人观看的视频www高清免费观看 | 男人舔奶头视频| 亚洲五月天丁香| 精品高清国产在线一区| 亚洲精品美女久久久久99蜜臀| 老司机午夜福利在线观看视频| 999精品在线视频| 亚洲一码二码三码区别大吗| 亚洲国产日韩欧美精品在线观看 | 亚洲精品av麻豆狂野| 男人舔女人下体高潮全视频| 国产97色在线日韩免费| 欧美三级亚洲精品| 日韩av在线大香蕉| 亚洲国产欧美网| 国产一卡二卡三卡精品| 校园春色视频在线观看| 国产男靠女视频免费网站| 麻豆av在线久日| 青草久久国产| 老熟妇乱子伦视频在线观看| 欧美日韩乱码在线| 欧美最黄视频在线播放免费| 日韩欧美三级三区| 午夜成年电影在线免费观看| 国产三级在线视频| 在线观看免费午夜福利视频| 亚洲中文字幕一区二区三区有码在线看 | 天天躁狠狠躁夜夜躁狠狠躁| 欧美日韩一级在线毛片| 日韩国内少妇激情av| 88av欧美| 午夜久久久久精精品| 久久人人精品亚洲av| 99在线视频只有这里精品首页| 午夜激情av网站| 久久久精品欧美日韩精品| 欧美久久黑人一区二区| 香蕉av资源在线| 国产v大片淫在线免费观看| 老汉色∧v一级毛片| 久久精品国产亚洲av高清一级| 亚洲精品av麻豆狂野| 国产精品久久久久久精品电影 | 黄色女人牲交| 99热这里只有精品一区 | 国产精品亚洲美女久久久| 长腿黑丝高跟| 国产精品一区二区免费欧美| 久久99热这里只有精品18| 性欧美人与动物交配| 国内少妇人妻偷人精品xxx网站 | 国产精品 国内视频| 亚洲av片天天在线观看| 久久久久久久久久黄片| 伦理电影免费视频| 国产av一区二区精品久久| 十分钟在线观看高清视频www| 国产成人欧美| 成人手机av| 激情在线观看视频在线高清| 日日干狠狠操夜夜爽| 久久精品国产综合久久久| 老汉色av国产亚洲站长工具| 精品一区二区三区av网在线观看| av欧美777| 亚洲免费av在线视频| 18禁观看日本| 亚洲成人久久爱视频| 精品国产亚洲在线| 桃红色精品国产亚洲av| 精品福利观看| 久久亚洲真实| 激情在线观看视频在线高清| 大型av网站在线播放| 亚洲一区二区三区不卡视频| 久久久久久久午夜电影| 精品久久久久久久人妻蜜臀av| 亚洲,欧美精品.| 久久伊人香网站| 国产又爽黄色视频| 性欧美人与动物交配| 日韩欧美国产一区二区入口| 非洲黑人性xxxx精品又粗又长| 精品一区二区三区视频在线观看免费| 亚洲五月婷婷丁香| 日本黄色视频三级网站网址| 这个男人来自地球电影免费观看| 亚洲五月色婷婷综合| 亚洲七黄色美女视频| 9191精品国产免费久久| 黑丝袜美女国产一区| 免费观看精品视频网站| 长腿黑丝高跟| 欧美久久黑人一区二区| 香蕉av资源在线| svipshipincom国产片| 中文字幕人成人乱码亚洲影| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成人精品中文字幕电影| 国产成人欧美在线观看| 亚洲成人精品中文字幕电影| 国产区一区二久久| 久久人人精品亚洲av| 久久久精品国产亚洲av高清涩受| 国产精品久久久久久亚洲av鲁大| 久99久视频精品免费| 麻豆久久精品国产亚洲av| 国产精华一区二区三区| 国产一区二区三区视频了| 一区二区三区精品91| 侵犯人妻中文字幕一二三四区| 免费搜索国产男女视频| 又黄又爽又免费观看的视频| 国产成人精品无人区| 一本大道久久a久久精品| 国产精品精品国产色婷婷| 真人做人爱边吃奶动态| 女性生殖器流出的白浆| 亚洲av成人一区二区三| 黑丝袜美女国产一区| 国产亚洲精品综合一区在线观看 | bbb黄色大片| 黄色 视频免费看| 两性夫妻黄色片| 亚洲电影在线观看av| 日韩高清综合在线| 国产午夜福利久久久久久| 日韩大码丰满熟妇| 国产99白浆流出| 给我免费播放毛片高清在线观看| 老司机靠b影院| 久久久国产成人免费| 老司机靠b影院| 麻豆国产av国片精品| 亚洲中文字幕一区二区三区有码在线看 | 男女做爰动态图高潮gif福利片| 亚洲男人的天堂狠狠| 精华霜和精华液先用哪个| 性欧美人与动物交配| 日韩成人在线观看一区二区三区| 人妻久久中文字幕网| 国产不卡一卡二| 亚洲avbb在线观看| 国产欧美日韩一区二区三| 国产久久久一区二区三区| 国产高清激情床上av| 国产久久久一区二区三区| 黄色女人牲交| 久久久久久久久免费视频了| 天堂动漫精品| 精品高清国产在线一区| 久久人妻福利社区极品人妻图片| 国产免费男女视频| 69av精品久久久久久| 欧美激情久久久久久爽电影| a级毛片在线看网站| 亚洲色图av天堂| a级毛片在线看网站| 变态另类成人亚洲欧美熟女| 禁无遮挡网站| 免费在线观看影片大全网站| 国产国语露脸激情在线看| 男人舔女人的私密视频| 桃色一区二区三区在线观看| 美女免费视频网站| 精品久久久久久成人av| 国产精品,欧美在线| 久久人妻av系列| 中文字幕人妻丝袜一区二区| 黄片播放在线免费| 国产精品日韩av在线免费观看| 热re99久久国产66热| 午夜福利免费观看在线| 又黄又粗又硬又大视频| 母亲3免费完整高清在线观看| 少妇 在线观看| 亚洲精华国产精华精| 久久香蕉精品热| 高清毛片免费观看视频网站| 黄色 视频免费看| 久久久国产成人免费| 好男人电影高清在线观看| 美女 人体艺术 gogo| 波多野结衣高清无吗| 国产成人一区二区三区免费视频网站| 村上凉子中文字幕在线| 法律面前人人平等表现在哪些方面| 他把我摸到了高潮在线观看| 怎么达到女性高潮| 日韩视频一区二区在线观看| 日本熟妇午夜| 18美女黄网站色大片免费观看| 国产精品亚洲一级av第二区| 1024手机看黄色片| 97碰自拍视频| 亚洲欧美精品综合久久99| 久久精品aⅴ一区二区三区四区| 国产精品永久免费网站| 亚洲国产欧洲综合997久久, | 两个人视频免费观看高清| 国产精品一区二区精品视频观看| 一个人观看的视频www高清免费观看 | 成熟少妇高潮喷水视频| 国产一卡二卡三卡精品| 免费一级毛片在线播放高清视频| 亚洲九九香蕉| 老熟妇乱子伦视频在线观看| 一区二区三区精品91| 国产一区二区三区视频了| 18禁黄网站禁片午夜丰满| 91成人精品电影| 99久久综合精品五月天人人| 宅男免费午夜| 精品高清国产在线一区| av视频在线观看入口| 久久国产精品影院| 中文字幕人妻熟女乱码| 国产野战对白在线观看| 国产成人啪精品午夜网站| 久久精品国产清高在天天线| 日日爽夜夜爽网站| 国产在线精品亚洲第一网站| 国产一区在线观看成人免费| 日韩欧美国产在线观看| 久久精品影院6| 精品卡一卡二卡四卡免费| 日本a在线网址| 黄色成人免费大全| 久9热在线精品视频| 91字幕亚洲| 中文字幕精品亚洲无线码一区 | 午夜福利视频1000在线观看| 黄色 视频免费看| 在线看三级毛片| 美女大奶头视频| 亚洲男人的天堂狠狠| 中文字幕人成人乱码亚洲影| 欧美国产日韩亚洲一区| 亚洲精品粉嫩美女一区| 国产成+人综合+亚洲专区| 国产精品 国内视频| 亚洲一区二区三区色噜噜| 国产麻豆成人av免费视频| 亚洲国产精品999在线| 亚洲欧洲精品一区二区精品久久久| 十八禁网站免费在线| 亚洲中文日韩欧美视频| 日日干狠狠操夜夜爽| 夜夜看夜夜爽夜夜摸| 中文字幕精品亚洲无线码一区 | 欧美日韩福利视频一区二区| 国产成人系列免费观看| 黑人欧美特级aaaaaa片| 亚洲av日韩精品久久久久久密| 国产成+人综合+亚洲专区| 免费人成视频x8x8入口观看| 亚洲 欧美 日韩 在线 免费| 无人区码免费观看不卡| 亚洲片人在线观看| 日日夜夜操网爽| 国产精品自产拍在线观看55亚洲| 欧美 亚洲 国产 日韩一| 亚洲黑人精品在线| 亚洲国产精品合色在线| 精品免费久久久久久久清纯| 国产精品乱码一区二三区的特点| 亚洲欧美一区二区三区黑人| 美国免费a级毛片| 欧美日韩亚洲综合一区二区三区_| 欧美乱色亚洲激情| 国内毛片毛片毛片毛片毛片| 又黄又粗又硬又大视频| 午夜亚洲福利在线播放| 国产免费男女视频| 成在线人永久免费视频| 国产熟女xx| 欧美日韩亚洲综合一区二区三区_| 久久久久久久久中文| a级毛片a级免费在线| 日本一本二区三区精品| 午夜免费观看网址| 黄色丝袜av网址大全| 黄色女人牲交| 国产av一区在线观看免费| 久久精品aⅴ一区二区三区四区| 午夜久久久在线观看| 国产又色又爽无遮挡免费看| 成熟少妇高潮喷水视频| 国产又色又爽无遮挡免费看| 欧美另类亚洲清纯唯美| 亚洲片人在线观看| 国产精品 欧美亚洲| 在线永久观看黄色视频| 午夜久久久久精精品| 国产精品免费一区二区三区在线| 午夜久久久久精精品| 亚洲男人天堂网一区| 欧美一级a爱片免费观看看 | 琪琪午夜伦伦电影理论片6080| 99国产精品99久久久久| 90打野战视频偷拍视频| 国产成人精品久久二区二区免费| 51午夜福利影视在线观看| 欧美日韩黄片免| 久久久久久亚洲精品国产蜜桃av| 免费在线观看完整版高清| 在线观看舔阴道视频| www.www免费av| 国产男靠女视频免费网站| 欧美大码av| 丝袜美腿诱惑在线| 亚洲精品在线美女| 首页视频小说图片口味搜索| 中文字幕人妻丝袜一区二区| 757午夜福利合集在线观看| 国产熟女午夜一区二区三区| 男人舔女人下体高潮全视频| 麻豆一二三区av精品| 男女视频在线观看网站免费 | 久久久久久国产a免费观看| 老熟妇乱子伦视频在线观看| 黄色片一级片一级黄色片| 久久婷婷人人爽人人干人人爱| 看黄色毛片网站| 国产精品 国内视频| 波多野结衣高清无吗| 午夜久久久久精精品| 欧美日韩亚洲综合一区二区三区_| 亚洲电影在线观看av| 久久草成人影院| 在线天堂中文资源库| 女性被躁到高潮视频| 亚洲成av人片免费观看| 午夜福利欧美成人| 日日摸夜夜添夜夜添小说| 国产精品 国内视频| 久久欧美精品欧美久久欧美| 最新美女视频免费是黄的| 国产亚洲精品久久久久5区| 欧美亚洲日本最大视频资源| 亚洲精品国产一区二区精华液| 亚洲av日韩精品久久久久久密| 一级黄色大片毛片| 欧美黄色淫秽网站| 国产av不卡久久| 国产精品一区二区精品视频观看| aaaaa片日本免费| 美女高潮喷水抽搐中文字幕| 亚洲国产精品合色在线| 黄色片一级片一级黄色片| 久久精品国产清高在天天线| 熟妇人妻久久中文字幕3abv| 别揉我奶头~嗯~啊~动态视频| 一区二区日韩欧美中文字幕| 精品卡一卡二卡四卡免费| 黄片大片在线免费观看| 成人三级黄色视频| 国产精品野战在线观看| 在线观看午夜福利视频| 成人三级黄色视频| 久久国产亚洲av麻豆专区| 日韩大码丰满熟妇| 国产精品久久久av美女十八| 丁香六月欧美| 18禁裸乳无遮挡免费网站照片 | 在线观看免费午夜福利视频| 在线永久观看黄色视频| 又黄又粗又硬又大视频| 欧美日韩福利视频一区二区| 亚洲激情在线av| 在线十欧美十亚洲十日本专区| 国产极品粉嫩免费观看在线| 1024香蕉在线观看| 亚洲最大成人中文| 国产成人一区二区三区免费视频网站| 亚洲第一青青草原| 岛国在线观看网站| aaaaa片日本免费| 免费高清视频大片| 少妇 在线观看| 精品久久蜜臀av无| 国产亚洲欧美在线一区二区| 黄网站色视频无遮挡免费观看| 色av中文字幕| 国产精品自产拍在线观看55亚洲| 欧美亚洲日本最大视频资源| 成人18禁在线播放| 18禁美女被吸乳视频|