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

    山區(qū)寬窄相間河段輸水能力變化規(guī)律

    2022-11-28 11:53:46閆旭峰姚致東旦增平措
    工程科學(xué)與技術(shù) 2022年6期
    關(guān)鍵詞:水躍相間沿程

    閆旭峰,姚致東,徐 輝,旦增平措

    (1.四川大學(xué) 水力學(xué)與山區(qū)河流開(kāi)發(fā)保護(hù)國(guó)家重點(diǎn)實(shí)驗(yàn)室,四川 成都 610065;2.四川大學(xué) 水利水電學(xué)院,四川 成都 610065;3.國(guó)家氣象中心,北京 100000;4.西藏自治區(qū)水文水資源勘測(cè)局,西藏 拉薩 850000)

    山區(qū)河流寬窄相間河段分布廣泛,相較于順直河道,河寬變化引起水流運(yùn)動(dòng)重新調(diào)整[1–2]。寬窄相間河段水流特性復(fù)雜,正確理解寬窄相間河道水流特性與輸水能力,是促進(jìn)山區(qū)河道水沙災(zāi)害防治的一項(xiàng)重要內(nèi)容[3–5]。國(guó)內(nèi)外眾多學(xué)者對(duì)河寬變化下的水沙運(yùn)動(dòng)特性進(jìn)行了研究,例如:閆旭峰等[2]通過(guò)寬窄相間水槽水流運(yùn)動(dòng)試驗(yàn),發(fā)現(xiàn)水位沿程分布呈現(xiàn)出明顯的展寬段壅水、窄段跌水特征,由窄到寬的展寬段易發(fā)生水躍現(xiàn)象。Valiani和Caleffi[6]采用線性角動(dòng)量守恒理論,得到了線性漸變水槽水躍空間結(jié)構(gòu)的解析解。Gandhi[7]基于突擴(kuò)型水槽水流運(yùn)動(dòng)試驗(yàn),分析了不同窄寬比的水躍尺度特征。王淑英[8]、Wang[9]、王文娥[10]等以水槽試驗(yàn)為基礎(chǔ),研究了寬窄相間河道橫向流速、紊動(dòng)特性及斷面環(huán)流分布。高永勝等[11]采用數(shù)值模擬方法,分析了河寬變化下的行洪過(guò)程及其河床演變特征。閆旭峰等[12]基于2維淺水方程模擬研究了寬窄相間河道因推移質(zhì)輸沙誘發(fā)床面調(diào)整的水面線時(shí)空變化特征,指出寬段和窄段存在不同的洪水過(guò)程模式。Wu[13]、Duró[14]、Nelson[15]等探究了寬窄相間河段局部河床演變特征。

    對(duì)于寬窄相間河道,窄段–寬段河寬比(窄寬比)是控制河段平面形態(tài)的主要參數(shù)。窄寬比越小,兩個(gè)斷面的過(guò)水面積相差越大,基于流量恒定條件,窄–寬段的流速差也就越大。水頭損失是影響河道輸水能力的主控因素,水頭損失越高,河道輸水能力就越低[16]。對(duì)于順直河道,水頭損失來(lái)源于邊界阻力;對(duì)于非順直河道,水頭損失不僅受邊界阻力影響,河道平面邊界形態(tài)的作用也非常大[17–19]。因此,研究寬窄相間河道水頭損失特性不僅有利于提高1維淺水方程的模擬精度[20],而且有助于理解該類型河道的輸水輸沙特性和機(jī)制。

    天然河流自然條件下,下游水位和寬–窄段河寬變化引起的河流形態(tài)變化是影響寬窄相間河段輸水特性的重要因素。本文基于平面2維水流數(shù)值模擬方法,探討下游水位與河寬變化如何影響寬窄相間河道水流特性,從而揭示寬窄相間河道輸水能力、空間形態(tài)非均勻阻力特征及洪水過(guò)程機(jī)制,為寬窄相間型河段的水沙災(zāi)害防治提供科學(xué)依據(jù)。

    1 平面2維水流數(shù)學(xué)模型及其驗(yàn)證

    采用2維淺水力學(xué)模型對(duì)寬窄相間河道水流運(yùn)動(dòng)進(jìn)行計(jì)算,2維淺水方程可表示為:

    連續(xù)方程:

    動(dòng)量方程:

    其中,

    式(1)~(5)中:h為 水深;t為 時(shí)間; ρ 為水的密度;U、V為水深平均流速;g為重力加速度;Zs為水面高程;τbx、 τby為床面切應(yīng)力;Cf為床面阻力系數(shù),由曼寧阻力系數(shù)n(與床面材料粗糙程度有關(guān),采用n=0.02)和水深表示,即Cf=gn2/h1/3; νt為渦流系數(shù),采用k–ε 模型進(jìn)行計(jì)算:

    式中,Cμ=0.09,k和 ε通過(guò)水深平均紊動(dòng)能運(yùn)輸方程及其耗散方程進(jìn)行描述:

    其中,

    式(7)~(11)中:C1ε=1.44,C2ε=1.92, σk=1.0,σε=1.3,Cεβ=3.6,Cε=1.3 ,均為經(jīng)驗(yàn)常數(shù)[21–22];U?=為摩擦流速。

    數(shù)值計(jì)算利用淺水模型,采用貼體網(wǎng)格對(duì)計(jì)算區(qū)進(jìn)行劃分,并對(duì)控制方程進(jìn)行相應(yīng)的轉(zhuǎn)換[23];采用有限差分法對(duì)控制方程進(jìn)行離散求解;采用3階顯性CIP算法進(jìn)行插值,精度高,穩(wěn)定性好[24–25]。模型驗(yàn)證利用河段室內(nèi)試驗(yàn)對(duì)沿程水位進(jìn)行比較,水槽試驗(yàn)在四川大學(xué)水力學(xué)與山區(qū)河流開(kāi)發(fā)保護(hù)國(guó)家重點(diǎn)實(shí)驗(yàn)室進(jìn)行,具體可參考文獻(xiàn)[2]。平面幾何形態(tài)及網(wǎng)格劃分如圖1(a)所示。水槽坡度S=0.002,最寬段寬度B= 1.4 m,最窄段寬度b= 0.6 m,寬窄段間距為3 m,并采用正弦曲線銜接。入口邊界條件采用流量邊界Q=0.1 m3/s,出口邊界條件采用水位邊界Hd= 0.148 m,水槽壁面為光滑水泥面,曼寧系數(shù)設(shè)為0.02,模型驗(yàn)證結(jié)果如圖1(b)所示。數(shù)值計(jì)算的水位沿程變化與試驗(yàn)值較為一致,誤差率±2.5%,表明該2維淺水模型及網(wǎng)格系統(tǒng)對(duì)于寬窄相間河道水流運(yùn)動(dòng)模擬具有良好的效果。此外,經(jīng)過(guò)網(wǎng)格收斂性分析,計(jì)算結(jié)果已經(jīng)收斂。

    圖1 寬窄相間河道數(shù)值模型設(shè)置及驗(yàn)證Fig.1 Sketch of diverging-converging channel and simulation result validation

    2 寬窄相間河段水流運(yùn)動(dòng)數(shù)值模擬結(jié)果分析

    利用驗(yàn)證后的模型研究窄寬比和下游出口水位對(duì)寬窄相間河道水流特性和輸水能力的影響,設(shè)計(jì)不同數(shù)值試驗(yàn)工況進(jìn)行敏感性分析。以驗(yàn)證工況為基準(zhǔn),保持流量和最寬段寬度不變,即Q=0.1 m3/s,B=1.4 m。通過(guò)改變窄段最小寬度b和下游出口水位Hd,研究窄寬比(b/B)變化與下游水位(Hd/B)效應(yīng)對(duì)寬窄相間河道水流特性的影響。設(shè)計(jì)7個(gè)窄寬比和8個(gè)不同水位,共56組試驗(yàn)工況。

    2.1 寬窄相間河段平均水位沿程變化

    圖2為窄寬比變化條件下斷面平均水位沿程H的分布與出口水位的關(guān)系。

    圖2 不同工況下寬窄相間河道斷面平均水位沿程分布Fig.2 Cross-sectional averaged water level profile along the channel under different scenarios

    河段呈現(xiàn)出典型的窄段跌水、寬段壅水;且當(dāng)下游水位效應(yīng)低時(shí),展寬段出現(xiàn)淹沒(méi)式水躍;跌水–壅水–水躍效應(yīng)隨著河道窄寬比增大而減小。當(dāng)b/B=5/7時(shí),水躍在給定淹沒(méi)條件下消失。當(dāng)出現(xiàn)水躍時(shí),水位沿程變化受下游出口水位變化影響較小,特別是上游第1個(gè)河段(窄—寬—窄)水面線基本重合,表明當(dāng)水躍出現(xiàn)時(shí),沿程水位分布主要受邊界平面形態(tài)影響,水躍位置出現(xiàn)在最窄段下游;而未產(chǎn)生水躍的跌水最低點(diǎn)則剛好出現(xiàn)在最窄斷面。相較而言,窄寬比越小,寬段的水位越高,表明壅水效應(yīng)越嚴(yán)重,洪水淹沒(méi)風(fēng)險(xiǎn)嚴(yán)重。隨著出口水位持續(xù)增大,上下游寬段水位在不同窄寬比條件下的差異變小,表明大水深情況下邊界平面形態(tài)對(duì)其影響較小,下游水位效應(yīng)顯著。而窄寬比越大(越接近于b/B=1),寬段水位隨下游出口水位變化幅度越大。

    2.2 寬窄相間河段平均流速沿程變化

    由于計(jì)算過(guò)程保持流量不變,因此,斷面平均流速主要由水位決定,水位越大,流速越小。圖3為不同工況下的斷面平均流速Um沿程分布。可知:斷面平均流速沿程分布基本與水位分布規(guī)律相反。全局高流速始終出現(xiàn)在窄斷面附近,即:如果產(chǎn)生水躍,最高流速則出現(xiàn)在最窄段下游,對(duì)應(yīng)水躍位置;如果未產(chǎn)生水躍,最高流速則出現(xiàn)在最窄斷面。窄寬比越小,寬、窄段流速差異越大(特別是寬段),且對(duì)下游出口水位變化不敏感,表明平面邊界形態(tài)效應(yīng)突出;窄寬比越大,下游水位效應(yīng)的主導(dǎo)作用越顯著。河寬縮窄整體上表現(xiàn)為輸水能力明顯降低,即寬段斷面流速在下游低水位時(shí)大幅降低,則相同流量對(duì)應(yīng)的水位更高,而這一輸水降低效應(yīng)隨著下游水位的增大而差異變小。此外,窄段河寬越小,最高流速越大,特別是產(chǎn)生水躍現(xiàn)象之后,這表明窄段河床沖刷越大;此時(shí)寬段流速較低,窄段沖刷出的泥沙較易淤積在寬段,形成深潭–淺灘結(jié)構(gòu)。在大水深條件下,窄段流速明顯降低,上游泥沙可能會(huì)重新淤積于窄段。

    圖3 不同工況下寬窄相間河道斷面平均流速沿程分布Fig.3 Cross-sectional averaged flow velocity profile along the channel under different scenarios

    2.3 寬窄相間河段水流流態(tài)沿程變化

    寬窄相間河道的跌水–壅水–水躍效應(yīng)是該平面邊界形態(tài)河道的主要特征。通常地,水躍的出現(xiàn)表明流態(tài)從急流(弗勞德數(shù)Fr>1)向緩流(Fr<1)轉(zhuǎn)變,因此,可利用Fr解釋寬窄相間河道窄段下游水躍現(xiàn)象。圖4為不同工況下的Fr沿程分布規(guī)律??芍篎r分布規(guī)律與斷面流速沿程分布規(guī)律相似,即窄段Fr較大,寬段較小;低窄寬比情況下,寬段Fr對(duì)于下游出口水位變化不敏感,而高窄寬比下Fr隨下游水位變化較大;最高Fr出現(xiàn)在最窄斷面或下游。計(jì)算結(jié)果還表明:當(dāng)Fr>1,即急流時(shí),水躍現(xiàn)象出現(xiàn);對(duì)于Fr<1或接近于1工況,水躍不會(huì)出現(xiàn),嚴(yán)格對(duì)應(yīng)水躍形成的臨界條件,即上游急流向下游緩流轉(zhuǎn)變。一般情況下,因急流向緩流轉(zhuǎn)變形成的水躍會(huì)出現(xiàn)于陡變緩的河道及因水工建筑物(閘門、溢流壩等)約束的流態(tài)下,而寬窄相間河道基于平面邊界形態(tài)的變化及下游出口水位效應(yīng)的影響產(chǎn)生了急流向緩流轉(zhuǎn)變而出現(xiàn)水躍的流態(tài)。

    圖4 不同工況下寬窄相間河道斷面平均弗勞德數(shù)沿程分布Fig.4 Cross-sectional averaged Froude number profile along the channel under different scenarios

    3 寬窄相間河道水流運(yùn)動(dòng)特性討論

    綜上可知,寬窄相間河道水流特性(水位、流速、弗勞德數(shù)等)受平面邊界形態(tài)和下游出口水位共同影響,出現(xiàn)窄段跌水、寬段壅水的水力現(xiàn)象,當(dāng)展寬段從急流(Fr>1)向緩流(Fr<1)轉(zhuǎn)變時(shí),出現(xiàn)淹沒(méi)式水躍現(xiàn)象。對(duì)于低窄寬比河道(即窄段河寬變?。?,水躍較易出現(xiàn),表明寬段易出現(xiàn)壅水。因此,當(dāng)河道由于滑坡等因素導(dǎo)致河道局部束窄,即使下游水位較低,亦會(huì)出現(xiàn)寬段壅水現(xiàn)象,導(dǎo)致洪水淹沒(méi)風(fēng)險(xiǎn)。而對(duì)于高窄寬比河道(即窄段河寬變大),寬段水位主要受下游出口水位影響,具有較為一致的同步性。

    寬窄相間河道存在寬段高水位、窄段低水位等特征水位,其水位分布特征受窄寬比和下游水位共同影響。通過(guò)分析其相關(guān)性,有助于理解寬窄相間河段上下游水位沿程分布規(guī)律。利用下游寬段水位(H3)對(duì)上游寬、窄段特征水位(H1、H2)進(jìn)行無(wú)量綱化,可以得到上游寬、窄段特征水位與下游寬段特征水位比(H1/H3、H2/H3)。圖5(a)為在高窄寬比條件下,上下游寬段特征水位比(H1/H3)隨下游出口水位的變化曲線,可知:H1/H3隨下游出口水位(Hd/B)的增大而降低,表明上下游河道因水位上漲形成了壅水效應(yīng);在低窄寬比條件下,H1/H3整體上接近于1,表明在出口低水位條件下,因河寬縮窄出現(xiàn)壅水效應(yīng),但隨下游出口水位的變化幅度不大。此外,低窄寬比條件下,H1/H3隨下游出口水位的增大呈現(xiàn)先增大后減小的趨勢(shì),表明低窄寬比寬窄相間河道可能存在最優(yōu)輸水效率。圖5(b)為在低窄寬比條件下,上游窄段與下游寬段特征水位比(H2/H3)隨下游出口水位的變化,可知:H2/H3先隨著下游出口水位的增大而幾乎保持不變甚至減??;超過(guò)一定水位之后,H2/H3隨出口水位的繼續(xù)增大而顯著增大。前者表現(xiàn)為水躍出現(xiàn)的水力條件,后者則對(duì)應(yīng)水躍消失的水力條件。

    圖5 上游寬(窄)段和下游寬段特征水位比與下游水位和河寬變化的關(guān)系Fig.5 Water level ratios of upstream wide (narrow) section to downstream wide section under different scenarios

    河道的輸水能力主要由河道阻力直接反映,當(dāng)河道阻力較大時(shí),河道輸水能力較低,而河道阻力可以由單位水頭損失表征。對(duì)于順直河道,水頭損失即能量損失主要由床面阻力、紊動(dòng)交換和斷面二次流造成[17,19]。但因河道平面邊界不變,河道沿程的水深變化較小,甚至形成均勻流,因此,順直河道的水頭損失與寬窄相間河道水頭損失研究機(jī)制不同。將縮窄–展寬段作為一個(gè)整體(圖1(a)中的x=3~9 m),進(jìn)行平均水頭損失分析,有助于從本質(zhì)上厘清寬窄相間河道的輸水機(jī)理。單位水頭損失計(jì)算方法如下。選擇在縮窄段中首先計(jì)算寬段和窄段的總水頭Ht,其計(jì)算公式為:

    式中,Htw1、Htn1(Htn2)、Htw2依次為x=3、6、9 m控制斷面總水頭,Hw1、Hn1、Hw2依次為x=3、6、9 m控制斷面水位,Uw1、Un1、Uw2依次為x=3、6、9 m控制斷面平均流速,g為重力加速度。縮窄段與展寬段的平均水頭損失為:

    式中,Scon、Sdiv依次為縮窄段、展寬段平均水頭損失,L1、L2依次為縮窄段、展寬段長(zhǎng)度。

    圖6為不同工況下的縮窄段與展寬段平均水頭損失變化規(guī)律。總體上看,展寬段的平均水頭損失明顯大于縮窄段,寬窄相間河道的水頭損失主要發(fā)生于展寬段,并與窄寬比及下游出口水位直接相關(guān)。由圖6(a)可知:在窄寬比變化條件下,水頭損失呈現(xiàn)出相似變化趨勢(shì),即隨下游出口水位增大先保持不變,然后逐漸減??;整體上,河道窄寬比越小,水頭損失越小,并最終由于下游出口水位的持續(xù)增大(Hd/B>0.14)使水頭損失收斂至同一個(gè)范圍。該變化趨勢(shì)表明,縮窄段由勢(shì)能轉(zhuǎn)變?yōu)閯?dòng)能,窄寬比越低,能量轉(zhuǎn)換效率越高,水頭損失較小。由圖6(b)可知:下游出口水位相同時(shí),窄寬比越小,水頭損失越高,這一趨勢(shì)與縮窄段相反。這是因?yàn)楫?dāng)動(dòng)能轉(zhuǎn)換為勢(shì)能時(shí),較小窄寬比可能會(huì)引起水躍這一高耗能水力現(xiàn)象。因此,可以發(fā)現(xiàn):當(dāng)下游水位(Hd/B)增大、較小窄寬比時(shí),展寬段的水躍效應(yīng)會(huì)不斷弱化甚至消失,相應(yīng)水頭損失會(huì)不斷降低并趨近于一致;進(jìn)一步比較發(fā)現(xiàn),對(duì)于低窄寬比河道,展寬段的水頭損失遠(yuǎn)大于收縮段,其比值最高約為6(圖6(c)),展寬段水頭損失占整個(gè)河段的近85%(圖6(d))。隨著水位的持續(xù)上升,展寬段水頭損失會(huì)大幅降低,對(duì)應(yīng)水躍現(xiàn)象的消失,此時(shí)展寬段水頭損失占全河段水頭損失比例大幅降低。相對(duì)而言,展寬段水頭損失維持在60%左右,始終高于縮窄段,但兩者差異在高水位條件下不斷縮小。

    圖6 縮窄段與展寬段平均水頭損失與下游水位和河寬變化的關(guān)系Fig.6 Head loss characteristics in diverging-converging channels under different scenarios

    4 結(jié) 論

    寬窄相間屬于天然河道平面形態(tài)的一種,其水流特性主要受河寬變化和下游水位影響,而當(dāng)前研究尚未對(duì)其輸水特性進(jìn)行系統(tǒng)探討。本文利用2維數(shù)值模型對(duì)不同條件下的寬窄相間型河道水流運(yùn)動(dòng)進(jìn)行了模擬,分析了不同窄寬比和下游水位效應(yīng)對(duì)河道輸水特性的影響。研究發(fā)現(xiàn):寬窄相間河道整體表現(xiàn)為寬段壅水低流速和窄段跌水高流速的水流特性;對(duì)于低窄寬比河道容易出現(xiàn)擴(kuò)段的淹沒(méi)式水躍,而下游出口水位上漲會(huì)抑制水躍的形成。河道窄寬比越小,寬段的壅水效應(yīng)越顯著,從而也提高了洪水淹沒(méi)的風(fēng)險(xiǎn)。通過(guò)分析收縮段–展寬段的水頭損失,發(fā)現(xiàn)展寬段的水頭損失始終大于收縮段,且其占全河段總水頭損失的比例隨著水躍的出現(xiàn)大幅提高,嚴(yán)重影響了河道的輸水輸沙能力,并對(duì)寬窄相間河段床面進(jìn)行塑造。本文從機(jī)理上分析了寬窄相間河道的輸水機(jī)制和洪水演進(jìn)及分布規(guī)律,可為山區(qū)河流水沙災(zāi)害防治提供科學(xué)依據(jù)。

    猜你喜歡
    水躍相間沿程
    那吉航運(yùn)樞紐大壩消力池?fù)p毀成因和搶修方案分析
    不同微納米曝氣滴灌入口壓力下迷宮流道沿程微氣泡行為特征
    典型生活垃圾爐排焚燒鍋爐沿程受熱面飛灰理化特性分析
    基于井下長(zhǎng)管線沿程阻力損失的計(jì)算研究
    突然擴(kuò)散水躍方程的改進(jìn)與比較
    新型分支相間導(dǎo)流排
    漸擴(kuò)式消力池水躍長(zhǎng)度的分析與計(jì)算
    變壓器相間過(guò)流保護(hù)不正確動(dòng)作原因的探討
    電子制作(2017年19期)2017-02-02 07:08:58
    矩形平底明渠水躍長(zhǎng)度公式的分析與應(yīng)用
    輸電線路相間距離保護(hù)應(yīng)用探討
    河南科技(2014年5期)2014-02-27 14:05:44
    亚洲性夜色夜夜综合| 午夜福利一区二区在线看| 亚洲少妇的诱惑av| 国产三级黄色录像| 一区在线观看完整版| 老司机午夜十八禁免费视频| 啦啦啦观看免费观看视频高清 | 亚洲成av片中文字幕在线观看| 狠狠狠狠99中文字幕| 久99久视频精品免费| 国产人伦9x9x在线观看| 国产精品一区二区精品视频观看| 精品国产美女av久久久久小说| 国内毛片毛片毛片毛片毛片| 亚洲av日韩精品久久久久久密| 99国产综合亚洲精品| 两个人视频免费观看高清| 国内毛片毛片毛片毛片毛片| 91av网站免费观看| 淫秽高清视频在线观看| 女性生殖器流出的白浆| 国产精品日韩av在线免费观看 | 无限看片的www在线观看| 亚洲国产欧美网| 国产成人精品久久二区二区免费| 亚洲国产欧美网| 超碰成人久久| 亚洲自偷自拍图片 自拍| 亚洲 欧美一区二区三区| 91九色精品人成在线观看| 亚洲人成77777在线视频| 精品高清国产在线一区| 国产在线精品亚洲第一网站| 国产在线精品亚洲第一网站| 欧美一级a爱片免费观看看 | 午夜福利高清视频| 中文字幕人妻熟女乱码| 欧美日韩亚洲综合一区二区三区_| 中文字幕久久专区| 午夜久久久在线观看| 两个人看的免费小视频| 国产精品电影一区二区三区| 好看av亚洲va欧美ⅴa在| 国产亚洲av高清不卡| 嫩草影视91久久| 天堂√8在线中文| 亚洲五月婷婷丁香| 脱女人内裤的视频| 人人妻,人人澡人人爽秒播| 色综合站精品国产| 很黄的视频免费| 91在线观看av| 午夜福利欧美成人| 国产成人影院久久av| 国内久久婷婷六月综合欲色啪| 欧美日韩精品网址| 国产精品免费一区二区三区在线| 久9热在线精品视频| 啦啦啦观看免费观看视频高清 | 久久性视频一级片| 日日夜夜操网爽| 少妇粗大呻吟视频| 琪琪午夜伦伦电影理论片6080| 99在线视频只有这里精品首页| 怎么达到女性高潮| 每晚都被弄得嗷嗷叫到高潮| 久久久精品欧美日韩精品| 美女 人体艺术 gogo| 人妻丰满熟妇av一区二区三区| 成人av一区二区三区在线看| 欧美激情极品国产一区二区三区| 国产单亲对白刺激| 精品少妇一区二区三区视频日本电影| 欧美黄色片欧美黄色片| 99精品久久久久人妻精品| 日日摸夜夜添夜夜添小说| 99久久久亚洲精品蜜臀av| 变态另类丝袜制服| 久久这里只有精品19| 亚洲av片天天在线观看| 欧美av亚洲av综合av国产av| 久久人妻熟女aⅴ| 自拍欧美九色日韩亚洲蝌蚪91| 欧美精品啪啪一区二区三区| 欧美在线黄色| 两性夫妻黄色片| 50天的宝宝边吃奶边哭怎么回事| 成人国产一区最新在线观看| 国产高清有码在线观看视频 | 最近最新免费中文字幕在线| 日韩大尺度精品在线看网址 | 亚洲人成伊人成综合网2020| 97超级碰碰碰精品色视频在线观看| 一进一出好大好爽视频| 欧美激情极品国产一区二区三区| 久久香蕉国产精品| 免费观看人在逋| 免费人成视频x8x8入口观看| 午夜福利视频1000在线观看 | 欧美 亚洲 国产 日韩一| 在线观看www视频免费| 亚洲自拍偷在线| 岛国视频午夜一区免费看| 久久天躁狠狠躁夜夜2o2o| 国产激情久久老熟女| 国产精品1区2区在线观看.| 黄色女人牲交| 制服诱惑二区| 咕卡用的链子| 亚洲国产看品久久| 麻豆国产av国片精品| 麻豆一二三区av精品| 波多野结衣高清无吗| av在线天堂中文字幕| 狂野欧美激情性xxxx| 激情在线观看视频在线高清| 人人妻人人爽人人添夜夜欢视频| 欧美色欧美亚洲另类二区 | 日韩大码丰满熟妇| 热re99久久国产66热| 男人舔女人的私密视频| 亚洲精品中文字幕一二三四区| 国产激情久久老熟女| 亚洲熟女毛片儿| 亚洲av成人av| 欧美中文日本在线观看视频| 一区二区三区激情视频| 亚洲性夜色夜夜综合| 中出人妻视频一区二区| 亚洲中文字幕日韩| 亚洲狠狠婷婷综合久久图片| 多毛熟女@视频| 在线观看www视频免费| 琪琪午夜伦伦电影理论片6080| 欧美丝袜亚洲另类 | 好男人在线观看高清免费视频 | av中文乱码字幕在线| 两个人视频免费观看高清| 亚洲欧洲精品一区二区精品久久久| 精品国产超薄肉色丝袜足j| 美女高潮到喷水免费观看| 最新在线观看一区二区三区| 成年版毛片免费区| 亚洲精品粉嫩美女一区| 精品久久久久久成人av| 一进一出抽搐gif免费好疼| 黑丝袜美女国产一区| 亚洲欧美精品综合久久99| 国产精品99久久99久久久不卡| 美女高潮喷水抽搐中文字幕| 久久草成人影院| 国语自产精品视频在线第100页| 欧美黄色片欧美黄色片| 免费在线观看亚洲国产| 国产亚洲精品久久久久久毛片| 欧美日本中文国产一区发布| 成人三级做爰电影| 黑人巨大精品欧美一区二区mp4| 一二三四社区在线视频社区8| 国产成人av教育| 久久天堂一区二区三区四区| 精品国产乱码久久久久久男人| 一级作爱视频免费观看| 国产aⅴ精品一区二区三区波| 欧美国产精品va在线观看不卡| 此物有八面人人有两片| 亚洲精品美女久久av网站| 国产精华一区二区三区| 在线播放国产精品三级| 天天躁夜夜躁狠狠躁躁| 亚洲精品中文字幕一二三四区| 国产日韩一区二区三区精品不卡| 波多野结衣高清无吗| 亚洲国产精品成人综合色| 亚洲情色 制服丝袜| 大型黄色视频在线免费观看| 人成视频在线观看免费观看| 久久青草综合色| 久久久国产欧美日韩av| 久9热在线精品视频| 国产免费男女视频| 亚洲av成人av| 法律面前人人平等表现在哪些方面| 国产一区二区三区视频了| 国产精品av久久久久免费| 亚洲精品国产色婷婷电影| 黑人巨大精品欧美一区二区蜜桃| 亚洲五月婷婷丁香| 国产亚洲精品av在线| 久久精品国产亚洲av香蕉五月| 一区在线观看完整版| 美国免费a级毛片| 一进一出好大好爽视频| 久久草成人影院| 久久久久久免费高清国产稀缺| 亚洲狠狠婷婷综合久久图片| 精品国产乱子伦一区二区三区| 久久人人精品亚洲av| 这个男人来自地球电影免费观看| 亚洲国产欧美网| 国产免费av片在线观看野外av| 99精品久久久久人妻精品| 又黄又爽又免费观看的视频| 欧美激情极品国产一区二区三区| 成熟少妇高潮喷水视频| 亚洲av成人不卡在线观看播放网| 久久精品成人免费网站| 日韩国内少妇激情av| 免费在线观看完整版高清| 日本撒尿小便嘘嘘汇集6| 九色亚洲精品在线播放| 久久中文字幕人妻熟女| 人人澡人人妻人| 涩涩av久久男人的天堂| 久久九九热精品免费| 国产伦一二天堂av在线观看| 美女高潮到喷水免费观看| 欧美中文综合在线视频| 亚洲成人精品中文字幕电影| 制服诱惑二区| 一级a爱视频在线免费观看| 国产亚洲欧美在线一区二区| 国产日韩一区二区三区精品不卡| 男女下面进入的视频免费午夜 | 成年人黄色毛片网站| 咕卡用的链子| 国产欧美日韩精品亚洲av| 中出人妻视频一区二区| 亚洲伊人色综图| a在线观看视频网站| 美国免费a级毛片| 深夜精品福利| www.精华液| 人人妻,人人澡人人爽秒播| 99riav亚洲国产免费| 国产成+人综合+亚洲专区| 亚洲av美国av| 国产精品永久免费网站| 国产成人系列免费观看| 国产国语露脸激情在线看| 久久久久国产精品人妻aⅴ院| 亚洲精品美女久久av网站| www.精华液| 国产av在哪里看| 桃红色精品国产亚洲av| 黄色成人免费大全| а√天堂www在线а√下载| 亚洲精品av麻豆狂野| 午夜a级毛片| x7x7x7水蜜桃| 亚洲国产欧美日韩在线播放| a在线观看视频网站| 免费在线观看视频国产中文字幕亚洲| 久久久久久久午夜电影| 亚洲在线自拍视频| 久久欧美精品欧美久久欧美| 国产精品自产拍在线观看55亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲在线自拍视频| 亚洲成人免费电影在线观看| 久久天堂一区二区三区四区| 亚洲色图 男人天堂 中文字幕| 亚洲 欧美一区二区三区| 乱人伦中国视频| 少妇裸体淫交视频免费看高清 | 欧美在线黄色| 亚洲精品久久国产高清桃花| 最新美女视频免费是黄的| 丝袜在线中文字幕| 女性生殖器流出的白浆| 久久久久久亚洲精品国产蜜桃av| 少妇熟女aⅴ在线视频| 国产高清激情床上av| 久久性视频一级片| 麻豆一二三区av精品| 村上凉子中文字幕在线| 满18在线观看网站| 成人18禁高潮啪啪吃奶动态图| 成人特级黄色片久久久久久久| 男女床上黄色一级片免费看| 高清毛片免费观看视频网站| av天堂久久9| 99久久久亚洲精品蜜臀av| 国产精品九九99| 免费不卡黄色视频| 又黄又爽又免费观看的视频| 午夜免费观看网址| 波多野结衣一区麻豆| 一级毛片精品| 一边摸一边做爽爽视频免费| 啦啦啦 在线观看视频| 两人在一起打扑克的视频| 亚洲午夜精品一区,二区,三区| 欧美激情 高清一区二区三区| 精品熟女少妇八av免费久了| 国产乱人伦免费视频| 国内毛片毛片毛片毛片毛片| 一个人免费在线观看的高清视频| 亚洲欧美日韩高清在线视频| 岛国视频午夜一区免费看| 久久久国产成人精品二区| 欧美不卡视频在线免费观看 | 在线播放国产精品三级| 在线观看www视频免费| www.999成人在线观看| 黄色片一级片一级黄色片| 91九色精品人成在线观看| 久久亚洲精品不卡| 青草久久国产| 9191精品国产免费久久| 亚洲国产看品久久| 国产色视频综合| 一区二区三区国产精品乱码| 日本欧美视频一区| 国产色视频综合| 国产亚洲精品综合一区在线观看 | 一区二区三区国产精品乱码| 啦啦啦 在线观看视频| 久久人人爽av亚洲精品天堂| 色综合亚洲欧美另类图片| 国产亚洲欧美精品永久| 97人妻精品一区二区三区麻豆 | 天堂影院成人在线观看| 黄色女人牲交| 亚洲精品国产色婷婷电影| 在线观看一区二区三区| 欧美另类亚洲清纯唯美| 欧美在线黄色| 美女扒开内裤让男人捅视频| 日韩 欧美 亚洲 中文字幕| 亚洲成国产人片在线观看| 国产av在哪里看| 无人区码免费观看不卡| 自线自在国产av| 久99久视频精品免费| 国产成人欧美| 国产熟女午夜一区二区三区| 亚洲成人久久性| 给我免费播放毛片高清在线观看| 亚洲中文字幕日韩| 久久九九热精品免费| 日本黄色视频三级网站网址| 久久性视频一级片| 亚洲精品中文字幕在线视频| 国产欧美日韩一区二区精品| ponron亚洲| 如日韩欧美国产精品一区二区三区| 日韩中文字幕欧美一区二区| 久久九九热精品免费| 一区二区三区激情视频| 后天国语完整版免费观看| 日韩中文字幕欧美一区二区| 99riav亚洲国产免费| 在线观看66精品国产| 母亲3免费完整高清在线观看| 国产精品乱码一区二三区的特点 | 成人特级黄色片久久久久久久| 视频在线观看一区二区三区| 日本黄色视频三级网站网址| 国产av又大| 不卡一级毛片| 免费在线观看黄色视频的| 亚洲成人久久性| 麻豆国产av国片精品| 黄色 视频免费看| 香蕉国产在线看| 久久精品91蜜桃| 国产亚洲av嫩草精品影院| 夜夜看夜夜爽夜夜摸| 18禁裸乳无遮挡免费网站照片 | 日韩精品青青久久久久久| 可以在线观看的亚洲视频| 精品电影一区二区在线| 精品国产一区二区久久| 国产av在哪里看| av片东京热男人的天堂| 搡老岳熟女国产| 色综合欧美亚洲国产小说| 97人妻精品一区二区三区麻豆 | 成人18禁在线播放| 涩涩av久久男人的天堂| 桃红色精品国产亚洲av| 中文字幕av电影在线播放| 每晚都被弄得嗷嗷叫到高潮| 一级片免费观看大全| 国产亚洲精品av在线| 精品乱码久久久久久99久播| 午夜老司机福利片| 欧美精品亚洲一区二区| 90打野战视频偷拍视频| 久久伊人香网站| 天天躁夜夜躁狠狠躁躁| 久久久久国内视频| 91九色精品人成在线观看| 欧美黑人欧美精品刺激| 校园春色视频在线观看| 亚洲专区中文字幕在线| 麻豆av在线久日| 老司机深夜福利视频在线观看| 日韩三级视频一区二区三区| 国产亚洲欧美精品永久| 亚洲av成人不卡在线观看播放网| 久久精品91无色码中文字幕| 黄色丝袜av网址大全| 亚洲一区二区三区色噜噜| 狠狠狠狠99中文字幕| 亚洲av片天天在线观看| 99国产精品免费福利视频| 久久中文字幕人妻熟女| 国产三级在线视频| 亚洲熟妇中文字幕五十中出| 后天国语完整版免费观看| 亚洲自偷自拍图片 自拍| 美女免费视频网站| 乱人伦中国视频| 久久久久久久久久久久大奶| 久久青草综合色| av天堂在线播放| 满18在线观看网站| 91精品国产国语对白视频| 欧美激情 高清一区二区三区| 欧美色欧美亚洲另类二区 | 岛国视频午夜一区免费看| 成人三级做爰电影| 精品人妻1区二区| 给我免费播放毛片高清在线观看| 一级,二级,三级黄色视频| 淫妇啪啪啪对白视频| 欧美激情极品国产一区二区三区| 9191精品国产免费久久| av视频在线观看入口| 9热在线视频观看99| 一级,二级,三级黄色视频| 久久精品91蜜桃| 精品熟女少妇八av免费久了| 久久人人精品亚洲av| 欧美成人免费av一区二区三区| 久热这里只有精品99| 午夜免费成人在线视频| 成人三级做爰电影| 久久中文看片网| 怎么达到女性高潮| 国产精品二区激情视频| 国产成人啪精品午夜网站| 午夜福利高清视频| 日韩成人在线观看一区二区三区| 欧美成狂野欧美在线观看| 91在线观看av| 亚洲熟妇中文字幕五十中出| 精品久久久久久成人av| 波多野结衣高清无吗| 久久久国产成人免费| 在线观看免费视频日本深夜| 亚洲aⅴ乱码一区二区在线播放 | 久久久久久大精品| 免费在线观看视频国产中文字幕亚洲| 欧美成人一区二区免费高清观看 | 久久精品国产99精品国产亚洲性色 | 巨乳人妻的诱惑在线观看| 多毛熟女@视频| www.999成人在线观看| 亚洲av成人一区二区三| 啪啪无遮挡十八禁网站| 欧美日韩黄片免| 人人妻人人澡人人看| 久久久国产欧美日韩av| av天堂久久9| 国产精品九九99| 欧美绝顶高潮抽搐喷水| 高清毛片免费观看视频网站| www.自偷自拍.com| 精品久久久久久成人av| 天堂√8在线中文| 亚洲午夜理论影院| 一个人观看的视频www高清免费观看 | 亚洲视频免费观看视频| 精品久久久久久久人妻蜜臀av | 日日摸夜夜添夜夜添小说| 日韩欧美免费精品| 亚洲专区国产一区二区| 18禁黄网站禁片午夜丰满| 一级a爱片免费观看的视频| 日韩欧美国产在线观看| 久久久久久国产a免费观看| 国产精品乱码一区二三区的特点 | 亚洲伊人色综图| 美女扒开内裤让男人捅视频| АⅤ资源中文在线天堂| 两性午夜刺激爽爽歪歪视频在线观看 | 国产主播在线观看一区二区| 欧美最黄视频在线播放免费| 18禁美女被吸乳视频| 亚洲三区欧美一区| 日韩一卡2卡3卡4卡2021年| 国产又色又爽无遮挡免费看| 50天的宝宝边吃奶边哭怎么回事| 香蕉丝袜av| 久热爱精品视频在线9| 久9热在线精品视频| 精品久久久久久,| 国产精品一区二区在线不卡| 精品久久久久久久毛片微露脸| 老熟妇乱子伦视频在线观看| 国产熟女午夜一区二区三区| 搡老岳熟女国产| 国产在线精品亚洲第一网站| 国内毛片毛片毛片毛片毛片| 亚洲中文字幕日韩| 国产区一区二久久| 亚洲第一青青草原| netflix在线观看网站| 母亲3免费完整高清在线观看| 12—13女人毛片做爰片一| 久久久久国产一级毛片高清牌| 国产一卡二卡三卡精品| 黄色毛片三级朝国网站| 久久精品人人爽人人爽视色| 午夜免费鲁丝| 热re99久久国产66热| 国产欧美日韩精品亚洲av| 淫妇啪啪啪对白视频| 99在线人妻在线中文字幕| 国产在线精品亚洲第一网站| 久久中文看片网| 高清在线国产一区| 国产一卡二卡三卡精品| 欧美黄色片欧美黄色片| 亚洲男人的天堂狠狠| 国产精品久久久人人做人人爽| 久久午夜综合久久蜜桃| 国产精品免费一区二区三区在线| 亚洲第一电影网av| 欧美黄色片欧美黄色片| 一级毛片精品| 极品教师在线免费播放| 亚洲一区二区三区不卡视频| 中文字幕另类日韩欧美亚洲嫩草| 搡老熟女国产l中国老女人| 久久久水蜜桃国产精品网| 黄片播放在线免费| 国产91精品成人一区二区三区| 久久久久精品国产欧美久久久| 国产私拍福利视频在线观看| 久久午夜综合久久蜜桃| 在线免费观看的www视频| 精品久久久久久久久久免费视频| 亚洲第一青青草原| 高潮久久久久久久久久久不卡| www.999成人在线观看| 啦啦啦韩国在线观看视频| 女生性感内裤真人,穿戴方法视频| 一级毛片高清免费大全| 大陆偷拍与自拍| 欧美国产精品va在线观看不卡| 两个人视频免费观看高清| 九色亚洲精品在线播放| 男人的好看免费观看在线视频 | 桃红色精品国产亚洲av| 亚洲成人免费电影在线观看| 精品一区二区三区四区五区乱码| av视频免费观看在线观看| netflix在线观看网站| 午夜免费激情av| 欧美日本亚洲视频在线播放| 精品一区二区三区av网在线观看| 最近最新中文字幕大全免费视频| 999久久久国产精品视频| 午夜久久久在线观看| 国产成人精品在线电影| 国产亚洲精品久久久久久毛片| 在线av久久热| 中文字幕人妻丝袜一区二区| 久久精品亚洲熟妇少妇任你| 久99久视频精品免费| 国产蜜桃级精品一区二区三区| 午夜日韩欧美国产| 在线观看免费日韩欧美大片| 欧美在线黄色| 亚洲中文日韩欧美视频| 久久久久国内视频| 欧美av亚洲av综合av国产av| 可以在线观看毛片的网站| 亚洲黑人精品在线| 精品人妻1区二区| 亚洲在线自拍视频| 免费人成视频x8x8入口观看| 亚洲精品在线观看二区| 黄片播放在线免费| 日韩av在线大香蕉| 日韩精品青青久久久久久| 最新美女视频免费是黄的| 亚洲精品中文字幕一二三四区| 日本免费一区二区三区高清不卡 | 又紧又爽又黄一区二区| 久久这里只有精品19| 欧美黄色片欧美黄色片| 成人国产综合亚洲| 老司机午夜十八禁免费视频| 精品乱码久久久久久99久播| 亚洲国产精品合色在线| x7x7x7水蜜桃| 在线免费观看的www视频| 国产精品国产高清国产av| 免费在线观看亚洲国产| 国产精品二区激情视频| 国产精品99久久99久久久不卡| 亚洲免费av在线视频| АⅤ资源中文在线天堂| 侵犯人妻中文字幕一二三四区| 狂野欧美激情性xxxx| 琪琪午夜伦伦电影理论片6080| 久久香蕉精品热| 成人亚洲精品av一区二区| 国产精品自产拍在线观看55亚洲| 国产成人欧美|