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

    曹妃甸二期圍海造地工程取沙物理模型設(shè)計及驗證

    2017-11-07 12:13:31嘯,毛寧,張
    海洋工程 2017年5期
    關(guān)鍵詞:岸灘波高曹妃甸

    徐 嘯,毛 寧,張 磊

    (南京水利科學(xué)研究院,江蘇 南京 210024)

    曹妃甸二期圍海造地工程取沙物理模型設(shè)計及驗證

    徐 嘯,毛 寧,張 磊

    (南京水利科學(xué)研究院,江蘇 南京 210024)

    為解決曹妃甸首鋼圍海造地工程用沙問題,計劃在曹妃甸甸頭東側(cè)灘涂取沙,因關(guān)系到曹妃甸岸灘的穩(wěn)定性,擬通過波浪泥沙物理模型進(jìn)行試驗研究,探討取沙對岸灘向/離岸運(yùn)動的影響?;诤涌诤0段锢砟P拖嗨评碚摚接懥四P捅瘸叩拇_定方法及模型沙的選??;在缺乏合適驗證資料的情況下,通過三種途徑對物理模型的合理性和可靠性進(jìn)行驗證試驗,驗證試驗表明,本模型可將二維沙灘剖面特性研究成果應(yīng)用于三維條件。

    曹妃甸;取沙;物理模型設(shè)計;驗證試驗;圍海造地;模型沙

    圖1 曹妃甸首鋼圍海造地擬選取沙區(qū)及物理模型布置Fig.1 The model layout and the selected sand pits

    首鋼京唐鋼鐵聯(lián)合有限責(zé)任公司圍海造地二期工程造地面積約10 km2,圍堤長度約15 km,吹填砂約4 000萬m3。設(shè)計單位擬在曹妃甸甸頭東側(cè)岸坡取沙,初步選定2#、4#及8#三個沙源區(qū),位于底高程+1~-12 m左右范圍,底質(zhì)d50在0.10~0.20 mm,即取沙區(qū)為沙質(zhì)岸灘,如圖1所示。波浪是該海域主要動力因素。為了研究不同取沙方案實施后對曹妃甸岸灘穩(wěn)定性影響程度(即侵蝕程度),在物理模型中需要模擬波浪條件下的泥沙向岸-離岸的輸移運(yùn)動??紤]到曹妃甸灘、槽緊鄰的特點(diǎn),物理模型還需滿足流場相似要求。

    1 取沙工程物理模型簡介

    如上所述,波浪作用是直接影響曹妃甸近岸岸灘穩(wěn)定性的主要動力條件,本研究重點(diǎn)是在曹妃甸特定自然條件下,各取沙方案的實施是否會對近岸岸灘造成較強(qiáng)的岸灘侵蝕,進(jìn)而影響曹妃甸岸灘的穩(wěn)定性;其次還需考慮是否會影響附近深槽即是否會發(fā)生較強(qiáng)的淤積。為此,在進(jìn)行物理模型設(shè)計時,主要考慮波浪相似(包括波浪近岸傳播變形相似、破碎波相似和波浪輸沙相似),同時也要滿足潮汐水流相似。

    1.1潮汐水流、泥沙運(yùn)動基本比尺關(guān)系

    1.2波浪模型相似要求

    波浪是導(dǎo)致沙質(zhì)岸灘近岸區(qū)(特別是破波區(qū))泥沙運(yùn)動的主要動力因素,根據(jù)本課題研究特點(diǎn),在模型中應(yīng)保證波浪折射相似、破波形態(tài)相似、波浪掀沙、橫向輸沙相似等要求[2]。

    1.2.1 波浪傳播運(yùn)動相似

    1) 折射相似[3]

    2) 破波形態(tài)相似[4]

    3) 波動水質(zhì)點(diǎn)運(yùn)動速度相似

    式中:λH為波高比尺,λL為波長比尺,λC為波速比尺,λT為波周期比尺,Um為床面最大波動速度。

    1.2.2 波浪條件下泥沙運(yùn)動相似

    1) 破波掀沙相似[5]

    式中:λρ為流體密度比尺,λρs為泥沙顆粒密實容重比尺。

    2) 碎波區(qū)內(nèi)岸灘剖面沖淤趨勢相似

    根據(jù)Hattori & Kawamata公式[6]:

    可導(dǎo)得波動條件下泥沙沉降速度比尺:

    3) 波浪條件下泥沙運(yùn)動方式相似

    1.3模型比尺及模型沙的選擇

    1.3.1 水平比尺λl

    根據(jù)場地條件及模型相似要求,經(jīng)過各方面因素綜合考慮,初步考慮模型水平比尺λl范圍為600~800。

    1.3.2 垂直比尺λh

    圖2 曹妃甸取沙區(qū)及測量斷面位置Fig.2 The selected sand pits and the measuring sections

    從水流角度考慮,為滿足模型與原型流態(tài)相似,需要做成變態(tài)模型,由張友齡公式可以算得:當(dāng)λl=600時,λh≤116;當(dāng)λl=800時,λh≤143。

    為滿足波浪傳播相似,需考慮表面張力相似要求,據(jù)研究,模型中波高不小于2 cm,周期不小于0.35 s,即可避免表面張力引起波浪衰減。

    首鋼二期圍海工程計劃取沙區(qū)附近近岸灘坡度為1/250~1/500(見圖2及表1),在實驗室條件下進(jìn)行如此平緩的坡度正態(tài)或小變率模型的試驗幾乎不可能,因為過分平緩的坡度會使波浪迅速衰減,也就無法保證波浪形態(tài)相似及破碎位置相似。結(jié)合現(xiàn)場波浪條件及以往試驗工作的經(jīng)驗,初步確定水平比尺λl=720,垂直比尺λh=120,即變率為6.0。模型布置情況見圖1。

    表1 取沙區(qū)附近部分?jǐn)嗝?~15 m等深線范圍內(nèi)岸灘平均坡度Tab.1 The mean slope of some sections in 0~15 m bathymetric range near the sand pit

    1.3.3 波高及波長比尺

    1.3.4 模型沙的選擇

    考慮到本項目重點(diǎn)研究近岸區(qū),特別是破波區(qū)岸灘沖淤變化,根據(jù)這一特點(diǎn)和要求,在選擇模型沙時需分別考慮以下相似要求:

    如前所述,主要研究區(qū)域位于+1~-10 m水深范圍;d50=0.10~0.20 mm,下面按dp=0.15 mm考慮,應(yīng)用以上相似關(guān)系式(6)、(7)來估算模型沙粒徑范圍,計算結(jié)果見表2。

    表中按沉降相似要求λω=λuλH/λl算得λω,然后應(yīng)用武漢水院統(tǒng)一沉速公式計算不同容重泥沙的對應(yīng)粒徑[1]。

    參考以上計算結(jié)果和以往波浪動床模型經(jīng)驗,模型中可采用顆粒密實容重γs=1.33 g/cm3,中徑d50=0.27 mm的煤粉作模型沙,這種規(guī)格煤粉干容重γo=0.7 g/cm3,沉速ω=0.86 cm/s。需要說明,個別組次曾采用天然沙(d50=0.12 mm)進(jìn)行比較試驗,后發(fā)現(xiàn)本次研究區(qū)域水深較大,大部分在破波區(qū)外,采用天然沙難以滿足泥沙起動相似要求,為此本文未予介紹。

    表2 不同容重模型沙計算結(jié)果(dp=0.15 mm)Tab.2 Calculation results of different bulk density model sand(dp=0.15 mm)

    1.3.5 含沙濃度比尺λs和沖淤時間比尺λt2

    根據(jù)以上選沙結(jié)果及有關(guān)比尺,可得:

    正如文獻(xiàn)[1]所指出,由于實際輸沙涉及因素十分復(fù)雜,不宜單純依據(jù)計算值來確定λt2,需要根據(jù)地形演變相似來確定λt2。物理模型主要比尺列于表3、表4。

    表3 主要比尺情況Tab.3 Main scales

    表4 泥沙比尺情況Tab.4 Sediment transport scales

    2 關(guān)于模型變率的討論[2]

    模型幾何變率不等于波浪變率,兩者可采用不同的變率。

    在模擬寬淺的海域時,整體模型需按幾何變態(tài)設(shè)計和制作,否則無法進(jìn)行縮尺模型模擬。而一般沙質(zhì)岸灘坡度為1/15~1/50,由于床面坡度較陡,床面阻力損耗與波浪紊動損耗相比是次一級的,因此常采用幾何正態(tài)模型或變率較小的變態(tài)模型。

    圖3為幾何正態(tài)模型中波浪自-15 m等深線向岸傳播時波浪衰減情況,計算初始波高H=1.2 m,波周期T=5.0 s,岸灘坡度按1/300考慮;可以看出,波浪沿程衰減十分迅速,很難保證破碎相似。圖4為幾何變態(tài)模型中波高沿程衰減情況,計算時水平比尺λI=720,垂直比尺分別考慮為720、360、240、120及80??梢钥闯觯怪北瘸擀薶=120時,基本可以保證近岸區(qū)波高衰減不到10%,而且試驗波周期也基本可以滿足生波機(jī)要求。

    要解決上述問題只有二種選擇,一是采用小比尺正態(tài)模型,要求模型做得很大,可能喪失縮尺模型的基本優(yōu)點(diǎn)。另一選擇就是采用幾何變態(tài)模型,但模型中波浪是正態(tài)或接近正態(tài)。這時根據(jù)需要滿足折射相似(取λL=λH)或繞射相似(λL=λl)。在動床模型情況下,如研究區(qū)域不位于波影區(qū),折射相似一般更重要。

    圖3 幾何正態(tài)模型中波高沿程衰減情況(Hp=1.2 m,TP=5 s,床面坡度1∶300,初始水深15 m)Fig.3 The wave height attenuation along the wave direction in an undistorted model

    圖4 幾何變態(tài)模型中波高沿程衰減情況(Hp=1.2 m,TP=5 s,床面坡度1∶300,初始水深15 m)Fig.4 The wave height attenuation along the wave direction in a distorted model

    3 驗證試驗

    為了掌握曹妃甸甸頭附近岸坡取沙對岸灘穩(wěn)定性的影響,需要在整體模型中模擬波動條件下泥沙向岸-離岸運(yùn)動,在這方面可以借鑒的經(jīng)驗和成果很少,同時工程現(xiàn)場適合物理模型驗證的資料也不多,而這時選擇和進(jìn)行合適可靠的驗證試驗直接關(guān)系到物理模型的可靠性。經(jīng)過對曹妃甸現(xiàn)場資料和泥沙運(yùn)動機(jī)理的反復(fù)研究,決定在物理模型中進(jìn)行以下三方面的驗證試驗:

    1)擬選沙源區(qū)試挖槽回淤試驗與理論結(jié)果的比較驗證以確定時間比尺。

    2)擬選沙源區(qū)平衡剖面特性的驗證試驗。

    3)擬選沙源區(qū)岸坡2003~2005年實際沖淤變化驗證試驗。

    3.1物理模型中試挖槽回淤驗證

    根據(jù)劉家駒研究,只要海岸泥沙運(yùn)動屬于懸移質(zhì)運(yùn)動,用于淤泥質(zhì)海岸航道的淤積計算式也可用于粉沙淤泥質(zhì)海岸,計算航道破波區(qū)淤積問題[5]。

    為此,在模型中布置試挖槽,測量大浪過程后挖槽內(nèi)的淤積,再與理論公式計算結(jié)果進(jìn)行比較,以檢驗和修正模型中泥沙沖淤時間比尺。

    劉家駒泥沙回淤計算關(guān)系式為:

    其次,在得到移動IoT設(shè)備分簇列表后,則需要計算K個UAV的最優(yōu)布署坐標(biāo),最小化大規(guī)模移動異構(gòu)物聯(lián)網(wǎng)的通信功率.由式(2)、(6)可得,

    模型中試挖槽橫穿2#、4#、8#沙源區(qū),從-15 m到+1 m垂直等深線布置,挖槽寬150 m,深5 m(見圖5)。

    “大浪”作用兩天后的理論關(guān)系式計算回淤厚度,與模型試驗?zāi)嗌郴赜俸穸鹊谋容^結(jié)果見圖6。由圖可以看出,在破波區(qū)范圍試驗結(jié)果與計算結(jié)果相當(dāng)一致;深水區(qū)計算值明顯偏小。分析認(rèn)為,深水區(qū)計算值偏小的原因是在沙質(zhì)床面條件下,深水區(qū)挖槽中有一定的底沙回淤量,而理論關(guān)系式只考慮懸沙回淤。

    以上試驗表明,模型可以較好地模擬現(xiàn)場波浪作用下的泥沙運(yùn)動,最后確定以煤粉為模型沙的泥沙沖淤時間比尺為λt2煤=400。

    圖5 試挖槽布置情況Fig.5 Trial tunnel layout

    圖6 大浪作用后試挖槽泥沙回淤分布與理論計算值比較Fig.6 The siltation distribution in a trial tunnel (the test results compared with theoretical value after high wave impact)

    3.2擬選沙源區(qū)剖面特征的驗證試驗

    由Hattori & Kawamata[2],即:

    根據(jù)擬選沙源區(qū)各處岸灘地形特點(diǎn)岸灘坡度(tgβ)和泥沙條件(ω),代入Hattori & Kawamata判數(shù)關(guān)系式,可以算得擬選沙源區(qū)對應(yīng)于不同剖面形態(tài)的“特征波高(H0)”,將之換算為模型試驗波要素。在物理模型中按照這些波要素進(jìn)行試驗,驗證沙源區(qū)岸灘剖面特征是否符合上述規(guī)律,以便從另一個角度證明模型的合理性。

    將擬選沙源區(qū)岸灘參數(shù),代入Hattori等判別式,可得沙源區(qū)西(D)、中(E)、東(F)三個斷面處對應(yīng)于典型剖面的特征波高,列于表5。

    由表5可以看出,位于沙源區(qū)西側(cè)的1#(D)斷面處,當(dāng)波高大于1.2 m(模型中為1.4 cm),即可能發(fā)生侵蝕。從泥沙起動條件看,物理模型中波高小于1.4 cm后,深水區(qū)泥沙即不能起動,因此,更小的波浪沒有實際意義。

    位于沙源區(qū)中間的3#(E)斷面處,當(dāng)深水波高為2.0 m時,將呈現(xiàn)“過渡型”剖面特征;當(dāng)深水波高大于2.7 m時,將呈現(xiàn)“侵蝕型”剖面特征。

    表5 沙源區(qū)各斷面處剖面特征波高Tab.5 The characteristic wave height in each cross-section in sand-taking area

    而位于沙源區(qū)東側(cè)的4#(F)斷面處,當(dāng)深水波高小于2.7 m時,將呈“淤積型”剖面特征;當(dāng)深水波高為4.3 m,呈“過渡型”剖面特征;當(dāng)深水波高大于6.2 m時,才呈現(xiàn)“侵蝕型”剖面特征。

    在模型中分別進(jìn)行了深水波高1.2、2.2和4.3 m的波浪泥沙動床試驗。

    試驗表明,在深水波高1.2 m的小浪條件下,1#(D)斷面近岸區(qū)依然發(fā)生侵蝕,符合“侵蝕型”岸灘剖面特征;3#(E)斷面近岸區(qū)發(fā)生淤積,也完全符合“淤積型”剖面特征;4#(F)斷面床面坡度平緩, 小浪在向岸行進(jìn)過程中沿程減弱,對近岸區(qū)泥沙的作用已十分微弱,床面已看不出明顯變化。

    圖7為深水波高2.2 m的中浪條件下岸灘沖淤趨勢。可以看出,此時1#(D)斷面近岸區(qū)發(fā)生較強(qiáng)侵蝕,為典型“侵蝕型”剖面;3#(E)斷面近岸區(qū)有淤有沖,基本屬于“微淤型”(Ⅲ-1型)剖面特征;在2.2 m波浪條件下,4#(F)斷面呈現(xiàn)為典型的“淤積型”剖面特征。

    圖7 深水波高2.2 m條件下沙源區(qū)3個主要斷面處沖淤趨勢Fig.7 The erosion trends of 3 main sections in sand-taking region under deep water wave height of 2.2 m

    深水波高4.3 m的大浪條件下,1#(D)斷面近岸2 500 m范圍均發(fā)生較強(qiáng)侵蝕,為“侵蝕型”剖面;3#(E)斷面上在-6 m水深處形成較大“沙壩(Bar)”,為典型的“Bar型剖面”,即“侵蝕型”剖面;最后,4#(F)斷面近岸區(qū)有淤有沖,基本呈“過渡型”剖面特征。

    以上驗證試驗結(jié)果分析表明,本波浪泥沙動床模型,較好地符合岸灘特征剖面理論,可以保證不同波浪條件下,岸灘沖淤趨勢的可靠性和合理性。

    3.3擬選沙源區(qū)岸坡2003~2005年地形沖淤變化驗證

    對已有地形測圖整理比較后發(fā)現(xiàn),2003年5月和2005年7月測圖《沙源調(diào)查水深》有部分范圍重疊,為動床模型驗證提供了難得的擬選沙源區(qū)近期實測地形變化資料。

    兩次測圖的間隔時間約2年,由于缺乏2003年~2005年曹妃甸海域波浪資料,根據(jù)1996年、1997年和1999年曹妃甸波浪資料統(tǒng)計出工程海域海向來浪年內(nèi)分布情況。具體試驗波浪要素見表6。現(xiàn)場0.6 m以下的波浪對泥沙運(yùn)動作用很小,模型中0.4 cm的波浪對模型沙運(yùn)動作用也很小,而且波高0.4 cm的波浪在模型試驗中也很難實現(xiàn),為此只考慮現(xiàn)場波高0.6 m以上波浪作用。

    表6 模型試驗波要素Tab.6 The wave elements in model experiments

    由于兩次測圖范圍的局限性,模型中只能驗證3條斷面(見圖8),驗證部分結(jié)果見圖9。模型沖淤趨勢及量級和現(xiàn)場實測資料基本一致,個別區(qū)域沖淤差別稍大,總體上,驗證結(jié)果良好。

    圖8 水深對比斷面Fig.8 Water depth comparison

    圖9 03-05A斷面地形驗證結(jié)果Fig.9 The results of 03-05A section terrain verification

    4 結(jié) 語

    1)首鋼造陸取沙區(qū)域位于曹妃甸甸頭東側(cè),由于取沙范圍大,可以預(yù)測在擬定區(qū)域取沙后將改變近岸波浪動力條件,由此可能影響附近岸灘和堤坊的穩(wěn)定性,因此研究是必要的。

    2)取沙區(qū)底質(zhì)為d50=0.10~0.20 mm沙質(zhì),有2大特點(diǎn):①取沙區(qū)岸灘等深線基本與主控浪向垂直,因此,泥沙運(yùn)動以向、離岸運(yùn)動為主;②擬取沙區(qū)坡度變化大,約1∶50~1∶500,三類波浪作用下泥沙向、離岸運(yùn)動均會發(fā)生。在兼顧潮流、波浪動力及相應(yīng)動力條件下泥沙相似的變態(tài)模型模擬難度非常大。

    3)從以上技術(shù)特點(diǎn)和難點(diǎn)出發(fā),研究首先基于正、變態(tài)物理模擬波浪折射的摩阻損耗所導(dǎo)致的相似和不相似性,進(jìn)行了不同坡度岸灘的正態(tài)、變態(tài)模型波浪傳播情況的比較,得到在一定比尺條件下,變態(tài)模型更易滿足折射相似的結(jié)論。其次,使用三重驗證:①劉家駒航槽回淤公式與試驗比較;②根據(jù)Hattori判別式對擬取沙區(qū)進(jìn)行不同波要素驗證試驗;③用2003~2005年實測地形變化進(jìn)行地形驗證試驗。證明了在λl=720、λh=120比尺條件下完全可以將二維沙灘剖面特性研究成果應(yīng)用于三維條件。

    [1] 武漢水利電力學(xué)院.河流泥沙工程學(xué)(下)[M].北京: 水利出版社,1982.(Wuhan University of Hydraulic and Electrical Engineering.River sedimentation engineering[M].Beijing: China Water & Power Press,1982.(in Chinese))

    [2] 徐嘯.波流共同作用下渾水動床整體模型的比尺設(shè)計及模型沙選擇[J].泥沙研究,1998(2):17-25.(XU Xiao.Scale design and model material selection for a turbid movable bed model under simultaneous actions of wave and current[J].Journal of Sediment Research,1998(2): 17-25.(in Chinese))

    [3] 河海大學(xué).海岸動力學(xué)[M].北京:交通出版社,1980.(Hohai University.Coastal dynamics[M].Beijing: China Communications Press,1980.(in Chinese))

    [4] BATTJES J A.Surf similarity[C]//Proceedings of the 14th Conf.CEC..1978:218-222.

    [5] 劉家駒.海岸泥沙運(yùn)動研究及應(yīng)用[M].北京:海洋出版社,2009.(LIU Jiaju.The study and application of coastal sediment movement[M].Beijing: China Ocean Press,2009.(in Chinese))

    [6] HATTORI M,KAWAMATA R.Onshore-offshore trandport and beach profile changes[C]//Proceedings of the 17th CEC..1980.

    [7] ENGELUND F.Turbulent energy and suspended load[R].Coastal Eng.Lab.,Tech.Univ.of Denmark,1965.

    Hydraulic model design and reliability analysis for sand excavation in Caofeidian

    XU Xiao,MAO Ning,ZHANG Lei

    (Nanjing Hydraulic Research Institute,Nanjing 210024,China)

    Sand excavation will be conducted on the east of Caofeidian east intertidal zone for Shougang Group,related to the stability of Caofeidian shore.A hydraulic model experiment has been carried out on the sand movement under waves off the beach,and the impact of the sand excavation analyzed,especially,the offshore transport of sand.Three methods are used for reliability analysis.The use of sediment 2-D motion of beach for 3-D is creative.

    Caofeidian; sand excavation; hydraulic model design; verification test; coastal reclamation; model sand

    TV148

    A

    10.16483/j.issn.1005-9865.2017.05.008

    1005-9865(2017)05-0071-08

    2017-02-28

    徐 嘯(1943-),男,江蘇淮安人,教授級高級工程師,主要從事海洋工程泥沙的研究。E-mail:xxu1943@163.com

    毛 寧。E-mail:nnao@nhri.cn

    猜你喜歡
    岸灘波高曹妃甸
    基于FHDI-GNWM 數(shù)據(jù)的全球超越概率波高宏觀分布特征分析
    風(fēng)暴浪作用下沙質(zhì)岸灘穩(wěn)定機(jī)制物理模型試驗研究*
    海洋與湖沼(2022年4期)2022-07-28 01:07:32
    中國銀行股份有限公司曹妃甸自貿(mào)區(qū)分行
    基于漂流浮標(biāo)的南大洋衛(wèi)星高度計有效波高研究
    海洋通報(2021年3期)2021-08-14 02:20:46
    非平整港池的多向不規(guī)則波試驗研究
    中國銀行股份有限公司曹妃甸分行
    崛起的曹妃甸
    岸灘溢油監(jiān)測評價指標(biāo)體系研究*
    飽和秋色
    女報seaside(2014年10期)2014-04-29 21:02:51
    潛堤影響下的沙質(zhì)岸灘剖面變化實驗研究
    海岸工程(2014年2期)2014-02-27 12:51:02
    国产老妇伦熟女老妇高清| 性高湖久久久久久久久免费观看| 日韩免费av在线播放| 国产成人精品久久二区二区免费| 欧美一级毛片孕妇| 国产精品久久久久久人妻精品电影 | 国产片内射在线| 亚洲一区二区三区欧美精品| 1024香蕉在线观看| av网站在线播放免费| 高清黄色对白视频在线免费看| 精品亚洲成国产av| 亚洲精品一卡2卡三卡4卡5卡| 国产高清国产精品国产三级| 99香蕉大伊视频| 日韩欧美免费精品| 久热这里只有精品99| 黄色 视频免费看| 成人永久免费在线观看视频 | 18禁观看日本| 麻豆乱淫一区二区| 精品国产国语对白av| av一本久久久久| 日本一区二区免费在线视频| 在线看a的网站| 国产精品九九99| 亚洲精品久久成人aⅴ小说| 三上悠亚av全集在线观看| 亚洲精品成人av观看孕妇| 午夜福利影视在线免费观看| 91大片在线观看| 久久久国产成人免费| 国产精品一区二区在线观看99| 亚洲国产看品久久| 国产无遮挡羞羞视频在线观看| 美女国产高潮福利片在线看| av国产精品久久久久影院| 美女扒开内裤让男人捅视频| 国产成人av教育| 中文字幕人妻熟女乱码| 亚洲av第一区精品v没综合| 欧美日韩亚洲综合一区二区三区_| 免费在线观看完整版高清| avwww免费| 亚洲欧美色中文字幕在线| 国产欧美日韩一区二区三区在线| 国产无遮挡羞羞视频在线观看| 飞空精品影院首页| 少妇粗大呻吟视频| 不卡一级毛片| 在线观看66精品国产| 国产精品亚洲一级av第二区| 一边摸一边抽搐一进一出视频| 丝袜美足系列| 黄色怎么调成土黄色| 亚洲欧洲日产国产| 老司机午夜福利在线观看视频 | 久久人妻av系列| 欧美大码av| 精品国内亚洲2022精品成人 | 国产精品久久久av美女十八| 悠悠久久av| 每晚都被弄得嗷嗷叫到高潮| 国产日韩欧美视频二区| 激情在线观看视频在线高清 | 国产成人av激情在线播放| 一级毛片电影观看| 国产精品 欧美亚洲| 成年女人毛片免费观看观看9 | 另类亚洲欧美激情| 国产又爽黄色视频| 狠狠精品人妻久久久久久综合| 蜜桃国产av成人99| 成年版毛片免费区| av电影中文网址| 热re99久久国产66热| 日日夜夜操网爽| 两个人看的免费小视频| 法律面前人人平等表现在哪些方面| 国产成人免费观看mmmm| 成人手机av| 午夜福利在线观看吧| 国产免费福利视频在线观看| 丝袜美腿诱惑在线| 性高湖久久久久久久久免费观看| 亚洲中文av在线| 热99久久久久精品小说推荐| 国产精品国产av在线观看| 一本一本久久a久久精品综合妖精| 黄片小视频在线播放| 国产不卡av网站在线观看| 精品一品国产午夜福利视频| 久久狼人影院| 国产欧美亚洲国产| 性高湖久久久久久久久免费观看| 国产精品自产拍在线观看55亚洲 | 他把我摸到了高潮在线观看 | 日本撒尿小便嘘嘘汇集6| 一本色道久久久久久精品综合| 高清av免费在线| 黄色a级毛片大全视频| 高清在线国产一区| 亚洲情色 制服丝袜| 啦啦啦视频在线资源免费观看| 欧美日韩国产mv在线观看视频| 精品少妇内射三级| 一本久久精品| 9热在线视频观看99| 国产日韩欧美视频二区| 国产精品久久电影中文字幕 | 日韩熟女老妇一区二区性免费视频| 国产在线免费精品| 18禁国产床啪视频网站| 国产在线免费精品| 老汉色av国产亚洲站长工具| 如日韩欧美国产精品一区二区三区| 视频在线观看一区二区三区| 成人手机av| 亚洲色图 男人天堂 中文字幕| 日韩有码中文字幕| 日本撒尿小便嘘嘘汇集6| 桃花免费在线播放| 精品午夜福利视频在线观看一区 | 正在播放国产对白刺激| 大片电影免费在线观看免费| 黑人猛操日本美女一级片| 搡老乐熟女国产| 国产免费视频播放在线视频| 亚洲专区字幕在线| 一级a爱视频在线免费观看| 亚洲精品在线美女| 精品亚洲乱码少妇综合久久| 男女无遮挡免费网站观看| 欧美激情极品国产一区二区三区| 亚洲精品中文字幕一二三四区 | 黄色视频在线播放观看不卡| 亚洲精品国产区一区二| 涩涩av久久男人的天堂| 少妇精品久久久久久久| 精品人妻熟女毛片av久久网站| 国内毛片毛片毛片毛片毛片| 夜夜爽天天搞| 波多野结衣一区麻豆| 首页视频小说图片口味搜索| 国产在线精品亚洲第一网站| 久久青草综合色| 高清欧美精品videossex| 老汉色∧v一级毛片| 亚洲国产看品久久| 天堂动漫精品| 久久久国产精品麻豆| 久久99一区二区三区| 一边摸一边抽搐一进一小说 | 亚洲性夜色夜夜综合| 国产福利在线免费观看视频| 久久精品91无色码中文字幕| 大香蕉久久成人网| 两性夫妻黄色片| 久久久精品国产亚洲av高清涩受| 国产成人免费观看mmmm| 亚洲精品国产色婷婷电影| 亚洲精品久久午夜乱码| 国产老妇伦熟女老妇高清| 人人妻,人人澡人人爽秒播| 亚洲欧美色中文字幕在线| 巨乳人妻的诱惑在线观看| 精品乱码久久久久久99久播| 国产一区二区三区在线臀色熟女 | 欧美精品啪啪一区二区三区| 又紧又爽又黄一区二区| 另类精品久久| 久久精品国产a三级三级三级| 色在线成人网| 久热这里只有精品99| 亚洲精品粉嫩美女一区| 大型av网站在线播放| 国产精品偷伦视频观看了| 极品人妻少妇av视频| 制服人妻中文乱码| 黄片播放在线免费| 热re99久久国产66热| 九色亚洲精品在线播放| 亚洲美女黄片视频| 深夜精品福利| 国内毛片毛片毛片毛片毛片| 亚洲天堂av无毛| av片东京热男人的天堂| 久久久精品94久久精品| 国产1区2区3区精品| 欧美日韩精品网址| 黄色 视频免费看| 日韩视频在线欧美| 人人妻人人爽人人添夜夜欢视频| 午夜福利在线免费观看网站| 狠狠狠狠99中文字幕| 国产成人影院久久av| a级片在线免费高清观看视频| 久久久久久人人人人人| 欧美av亚洲av综合av国产av| 亚洲成人国产一区在线观看| 亚洲中文av在线| 亚洲一码二码三码区别大吗| 日本av免费视频播放| 亚洲国产欧美在线一区| 欧美久久黑人一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 高清毛片免费观看视频网站 | 手机成人av网站| 亚洲国产欧美日韩在线播放| 涩涩av久久男人的天堂| 国产一卡二卡三卡精品| 在线天堂中文资源库| 黑人巨大精品欧美一区二区蜜桃| 亚洲一区中文字幕在线| 每晚都被弄得嗷嗷叫到高潮| 波多野结衣av一区二区av| 久久中文看片网| 国产有黄有色有爽视频| 在线观看免费午夜福利视频| 免费一级毛片在线播放高清视频 | 婷婷成人精品国产| www.999成人在线观看| 日韩大码丰满熟妇| 99久久精品国产亚洲精品| 久久青草综合色| 侵犯人妻中文字幕一二三四区| 999久久久国产精品视频| 国产真人三级小视频在线观看| 91成年电影在线观看| 两人在一起打扑克的视频| 99久久99久久久精品蜜桃| 亚洲全国av大片| 欧美黄色片欧美黄色片| a级毛片在线看网站| 一级片'在线观看视频| 亚洲午夜理论影院| 日韩成人在线观看一区二区三区| 亚洲七黄色美女视频| 极品人妻少妇av视频| 五月天丁香电影| 丁香六月天网| 亚洲,欧美精品.| 在线看a的网站| 国产一区二区在线观看av| 国产高清国产精品国产三级| 国产老妇伦熟女老妇高清| 女人被躁到高潮嗷嗷叫费观| 午夜福利一区二区在线看| 狠狠婷婷综合久久久久久88av| 桃花免费在线播放| 午夜91福利影院| 丁香六月天网| 18禁美女被吸乳视频| 亚洲国产欧美在线一区| 国产精品久久久久久精品电影小说| www.自偷自拍.com| 黄色视频不卡| 制服诱惑二区| 欧美中文综合在线视频| 一区二区三区国产精品乱码| 少妇 在线观看| 国产男靠女视频免费网站| 欧美黑人欧美精品刺激| 免费女性裸体啪啪无遮挡网站| 欧美精品一区二区免费开放| 自拍欧美九色日韩亚洲蝌蚪91| 窝窝影院91人妻| 一级a爱视频在线免费观看| 97在线人人人人妻| a级片在线免费高清观看视频| 免费女性裸体啪啪无遮挡网站| 男人操女人黄网站| 亚洲欧洲日产国产| 咕卡用的链子| 动漫黄色视频在线观看| 国产精品香港三级国产av潘金莲| 久久精品人人爽人人爽视色| 18在线观看网站| 亚洲免费av在线视频| 色视频在线一区二区三区| 免费日韩欧美在线观看| 91字幕亚洲| 考比视频在线观看| 叶爱在线成人免费视频播放| 免费观看a级毛片全部| 国产精品 欧美亚洲| 黑人操中国人逼视频| 激情视频va一区二区三区| 国产淫语在线视频| 日韩熟女老妇一区二区性免费视频| 亚洲精品自拍成人| 久久久久久亚洲精品国产蜜桃av| 丰满迷人的少妇在线观看| 亚洲专区字幕在线| 欧美日韩亚洲高清精品| 国产欧美日韩一区二区三| 51午夜福利影视在线观看| 黄色片一级片一级黄色片| 制服人妻中文乱码| 曰老女人黄片| 涩涩av久久男人的天堂| 成年人免费黄色播放视频| 亚洲一区二区三区欧美精品| 日韩大片免费观看网站| 黄色a级毛片大全视频| 欧美日韩av久久| 新久久久久国产一级毛片| 91精品三级在线观看| 啪啪无遮挡十八禁网站| 国产一区二区在线观看av| 黄色怎么调成土黄色| 女同久久另类99精品国产91| 1024香蕉在线观看| 99riav亚洲国产免费| 9色porny在线观看| 亚洲成人免费电影在线观看| 18禁国产床啪视频网站| a在线观看视频网站| 欧美大码av| 建设人人有责人人尽责人人享有的| 在线永久观看黄色视频| 麻豆av在线久日| 久久午夜综合久久蜜桃| 国产老妇伦熟女老妇高清| 99国产极品粉嫩在线观看| 成年人免费黄色播放视频| 别揉我奶头~嗯~啊~动态视频| 丰满迷人的少妇在线观看| 久热爱精品视频在线9| 一区福利在线观看| 亚洲av第一区精品v没综合| 极品人妻少妇av视频| 久久av网站| 亚洲av国产av综合av卡| 国产精品欧美亚洲77777| 久久精品亚洲熟妇少妇任你| 99久久人妻综合| 亚洲中文av在线| 国产熟女午夜一区二区三区| 两个人看的免费小视频| 亚洲专区国产一区二区| 日韩精品免费视频一区二区三区| 中文字幕人妻丝袜制服| aaaaa片日本免费| 免费一级毛片在线播放高清视频 | 王馨瑶露胸无遮挡在线观看| 老熟妇仑乱视频hdxx| 一二三四社区在线视频社区8| 最近最新免费中文字幕在线| 午夜日韩欧美国产| 国产91精品成人一区二区三区 | 精品人妻1区二区| 久久精品91无色码中文字幕| 大片免费播放器 马上看| 国产免费av片在线观看野外av| 国产一区二区在线观看av| 亚洲欧美激情在线| 俄罗斯特黄特色一大片| 欧美精品av麻豆av| 国产日韩欧美亚洲二区| cao死你这个sao货| 黑人欧美特级aaaaaa片| 亚洲伊人色综图| 国产精品久久久久久人妻精品电影 | 又紧又爽又黄一区二区| 亚洲七黄色美女视频| 免费日韩欧美在线观看| 电影成人av| av免费在线观看网站| a级毛片在线看网站| 国产精品免费视频内射| 大香蕉久久成人网| 12—13女人毛片做爰片一| 色精品久久人妻99蜜桃| av不卡在线播放| 久久精品国产亚洲av高清一级| av不卡在线播放| 久久精品国产亚洲av高清一级| 超碰成人久久| 99国产精品一区二区蜜桃av | 精品人妻在线不人妻| 多毛熟女@视频| 啦啦啦免费观看视频1| 狂野欧美激情性xxxx| 国产精品久久久av美女十八| 亚洲国产中文字幕在线视频| 欧美午夜高清在线| 日本av手机在线免费观看| 久久久水蜜桃国产精品网| 日日爽夜夜爽网站| 国产av精品麻豆| 欧美人与性动交α欧美软件| 精品福利永久在线观看| 侵犯人妻中文字幕一二三四区| 99riav亚洲国产免费| 精品国产亚洲在线| 麻豆国产av国片精品| 欧美 日韩 精品 国产| 国产极品粉嫩免费观看在线| 人人妻,人人澡人人爽秒播| 桃花免费在线播放| 啦啦啦免费观看视频1| 国产激情久久老熟女| 另类亚洲欧美激情| xxxhd国产人妻xxx| 欧美激情 高清一区二区三区| 久久国产精品人妻蜜桃| 久久中文看片网| 日本av手机在线免费观看| 国产在线免费精品| 一本久久精品| 老鸭窝网址在线观看| 亚洲欧美激情在线| 在线观看免费午夜福利视频| videos熟女内射| 91麻豆av在线| 91麻豆精品激情在线观看国产 | av网站在线播放免费| avwww免费| 久久久精品94久久精品| 欧美黄色片欧美黄色片| 真人做人爱边吃奶动态| 中文字幕另类日韩欧美亚洲嫩草| 国产精品一区二区精品视频观看| 超色免费av| 不卡一级毛片| 极品人妻少妇av视频| 亚洲欧美色中文字幕在线| 亚洲欧洲日产国产| 精品少妇黑人巨大在线播放| 午夜91福利影院| 满18在线观看网站| 亚洲欧美精品综合一区二区三区| 男人操女人黄网站| 新久久久久国产一级毛片| 自拍欧美九色日韩亚洲蝌蚪91| 欧美变态另类bdsm刘玥| 亚洲色图综合在线观看| 日韩欧美免费精品| 久热爱精品视频在线9| 超色免费av| 在线亚洲精品国产二区图片欧美| 亚洲国产av新网站| 极品少妇高潮喷水抽搐| 高清av免费在线| 国产一卡二卡三卡精品| 女人爽到高潮嗷嗷叫在线视频| 国产精品av久久久久免费| 国产老妇伦熟女老妇高清| av超薄肉色丝袜交足视频| 精品国产亚洲在线| 极品人妻少妇av视频| 国产xxxxx性猛交| 国产人伦9x9x在线观看| 美女高潮到喷水免费观看| 国产在视频线精品| 69av精品久久久久久 | 热99re8久久精品国产| 色在线成人网| 久久久久精品国产欧美久久久| 久久久久精品人妻al黑| 成人永久免费在线观看视频 | 久久久久久久国产电影| av又黄又爽大尺度在线免费看| 久热这里只有精品99| 人人澡人人妻人| 成人免费观看视频高清| 亚洲一区中文字幕在线| 如日韩欧美国产精品一区二区三区| 韩国精品一区二区三区| 国产精品久久电影中文字幕 | 国内毛片毛片毛片毛片毛片| 一级黄色大片毛片| 免费在线观看黄色视频的| 欧美av亚洲av综合av国产av| 久久影院123| 成人免费观看视频高清| 夜夜骑夜夜射夜夜干| 久久亚洲真实| 中文字幕制服av| 老司机午夜福利在线观看视频 | 高清av免费在线| 欧美黑人欧美精品刺激| 一本久久精品| 欧美亚洲 丝袜 人妻 在线| 亚洲国产成人一精品久久久| 色综合欧美亚洲国产小说| 亚洲,欧美精品.| 午夜日韩欧美国产| 亚洲精华国产精华精| 久久久久久人人人人人| 成人三级做爰电影| 91精品三级在线观看| av电影中文网址| 中国美女看黄片| 天天躁夜夜躁狠狠躁躁| www.精华液| 午夜久久久在线观看| 人妻一区二区av| 色综合婷婷激情| 久久亚洲真实| 亚洲全国av大片| 欧美日韩亚洲国产一区二区在线观看 | 亚洲少妇的诱惑av| av一本久久久久| 欧美激情极品国产一区二区三区| 成年人午夜在线观看视频| 欧美日韩视频精品一区| 免费少妇av软件| 黄色a级毛片大全视频| 精品国产乱码久久久久久男人| 国产视频一区二区在线看| 国产麻豆69| 一区二区日韩欧美中文字幕| 日日摸夜夜添夜夜添小说| 男女下面插进去视频免费观看| 美女视频免费永久观看网站| 狂野欧美激情性xxxx| 人妻一区二区av| av电影中文网址| 高清毛片免费观看视频网站 | 老司机深夜福利视频在线观看| 高清毛片免费观看视频网站 | 一区二区三区国产精品乱码| 老鸭窝网址在线观看| 999久久久国产精品视频| 国产精品免费大片| 欧美日韩中文字幕国产精品一区二区三区 | 欧美av亚洲av综合av国产av| 国产亚洲午夜精品一区二区久久| 欧美激情极品国产一区二区三区| 下体分泌物呈黄色| 18禁黄网站禁片午夜丰满| 亚洲国产看品久久| 久久 成人 亚洲| 久久亚洲真实| 乱人伦中国视频| 久久ye,这里只有精品| 天天添夜夜摸| 建设人人有责人人尽责人人享有的| 桃花免费在线播放| 人人妻人人爽人人添夜夜欢视频| 不卡av一区二区三区| 亚洲色图 男人天堂 中文字幕| 丁香欧美五月| 日日夜夜操网爽| 夜夜爽天天搞| 久久精品国产a三级三级三级| 正在播放国产对白刺激| 亚洲九九香蕉| 久久久久久免费高清国产稀缺| 成年女人毛片免费观看观看9 | 中文字幕人妻丝袜制服| 国产精品久久久久久精品古装| 欧美成人免费av一区二区三区 | 最近最新中文字幕大全免费视频| 欧美 日韩 精品 国产| 免费久久久久久久精品成人欧美视频| 欧美日韩中文字幕国产精品一区二区三区 | 人妻 亚洲 视频| 中文字幕高清在线视频| 汤姆久久久久久久影院中文字幕| 一个人免费看片子| 人人妻人人添人人爽欧美一区卜| 大型av网站在线播放| 91字幕亚洲| 亚洲欧洲精品一区二区精品久久久| 久久中文字幕人妻熟女| 婷婷成人精品国产| 美女扒开内裤让男人捅视频| 色94色欧美一区二区| 在线观看66精品国产| 亚洲免费av在线视频| 精品少妇黑人巨大在线播放| 成年版毛片免费区| 老司机午夜十八禁免费视频| 欧美黑人欧美精品刺激| a级毛片黄视频| 国产精品秋霞免费鲁丝片| 1024香蕉在线观看| 不卡av一区二区三区| 99热国产这里只有精品6| 久久精品aⅴ一区二区三区四区| 9191精品国产免费久久| 成年人黄色毛片网站| 精品视频人人做人人爽| 亚洲国产欧美网| 大香蕉久久成人网| 美女午夜性视频免费| 黄色毛片三级朝国网站| 视频在线观看一区二区三区| 亚洲国产毛片av蜜桃av| 国产成人免费无遮挡视频| 最新的欧美精品一区二区| 丁香欧美五月| 一区二区三区国产精品乱码| 嫁个100分男人电影在线观看| 日本欧美视频一区| 精品高清国产在线一区| 91国产中文字幕| avwww免费| 丰满少妇做爰视频| 久9热在线精品视频| 国产精品国产av在线观看| 亚洲自偷自拍图片 自拍| 亚洲人成电影免费在线| 亚洲精品一二三| 欧美精品啪啪一区二区三区| 日本wwww免费看| 国产精品二区激情视频| 欧美精品亚洲一区二区| av片东京热男人的天堂| 99在线人妻在线中文字幕 | 免费在线观看完整版高清| 亚洲午夜精品一区,二区,三区|