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

    那蘇瓦卡里水電站沉沙池沉沙及沖沙效果研究

    2022-04-16 10:48:06高東紅鄒浩
    人民長江 2022年3期
    關(guān)鍵詞:沉沙池泥沙粒徑

    高東紅 鄒浩

    摘要: 在水利水電工程中一般都采用沉沙池來降低泥沙含量,減小泥沙對引水發(fā)電系統(tǒng)產(chǎn)生的不利影響,沉沙池中含沙水流流態(tài)及沉沙沖沙效果的研究就顯得極為重要?;谀岵礌柲翘K瓦卡里水電工程的沉沙池,采用理論公式、RNG ?k~ε 紊流模型及體積分辨率函數(shù)VOF模擬計(jì)算了沉沙池內(nèi)水沙兩相流的水流形態(tài),并詳細(xì)展示了多組分泥沙顆粒組合中不同顆粒直徑泥沙的沉積效果。結(jié)果表明:不小于0.2 mm粒徑的泥沙沉沙率達(dá)到95%以上,數(shù)值計(jì)算結(jié)果與采用理論公式計(jì)算得出的沉沙率保持較高的吻合性,滿足工程設(shè)計(jì)要求。

    關(guān) 鍵 詞: 沉沙池; 沉沙率; 數(shù)值計(jì)算; 那蘇瓦卡里水電站; 尼泊爾

    中圖法分類號(hào): ?TV68

    文獻(xiàn)標(biāo)志碼: ?A

    DOI: 10.16232/j.cnki.1001-4179.2022.03.019

    ?0 引 言

    自然界河流中多含有泥沙,水流中的泥沙通常會(huì)對水利水電工程帶來諸多不利影響。近年來部分學(xué)者陸續(xù)對三峽水庫、小浪底水庫的泥沙淤積? [1-2] 、水工建筑物及水輪機(jī)組磨損? [3] 等問題開展了研究。

    19世紀(jì)開始,含沙水流的理論研究就取得了飛速發(fā)展。Duboys? [4] 提出了推移質(zhì)運(yùn)動(dòng)的拖拽力理論,Rouse? [5] 則提出了著名的懸移質(zhì)泥沙濃度分布公式。中國學(xué)者竇國仁? [6] 則提出了表示縱向脈動(dòng)流速沿豎直線分布的公式;蔡樹棠? [7] 從流體力學(xué)基本理論出發(fā),討論了泥沙在湍流中對水相對速度、脈動(dòng)速度的改變,泥沙在二維渠道均勻定常流動(dòng)中含沙量的垂直分布和通常的擴(kuò)散理論修正等3個(gè)問題。

    從本質(zhì)上講,含沙水流屬于固液兩相流理論研究范疇,泥沙數(shù)值模擬主要基于兩相流理論。隨著計(jì)算機(jī)技術(shù)的發(fā)展,數(shù)值計(jì)算越來越廣泛地用于含沙水流計(jì)算中。中國西部地區(qū)及在世界眾多國家的河流中均含有大量泥沙,泥沙問題成為眾多水利水電工程的制約因素。解決引水系統(tǒng)中的泥沙問題,通常是依靠修建沉沙池及其沖沙系統(tǒng),而沉沙池中水沙流動(dòng)分析常借助于數(shù)值計(jì)算。在沉沙池?cái)?shù)值計(jì)算方面,楊紅等? [8] 進(jìn)行了沉沙池立面流場模擬計(jì)算,胡斌等? [9] 進(jìn)行了沉沙池的三維流場模擬,周曉泉、張軍等? [10-11] 分別對復(fù)雜邊界條件下及圓中環(huán)沉沙池的三維紊流場進(jìn)行了數(shù)值模擬研究。

    尼泊爾地區(qū)多數(shù)為多泥沙河流,修建水電站工程通常需考慮沉沙池。岳竹等? [12] 等研究了尼泊爾內(nèi)有壓隧洞內(nèi)沉沙池泥沙沉降特性。本文結(jié)合那蘇瓦卡里水電站工程,對沖沙設(shè)施的結(jié)構(gòu)設(shè)計(jì)理論及布置進(jìn)行了詳細(xì)論述,并運(yùn)用大型流體計(jì)算軟件對沉沙池中多組分泥沙顆粒的流態(tài)及沉沙效果進(jìn)行數(shù)值模擬驗(yàn)算,選用較為成熟的湍流模型 RNG ?k~ε 模型,應(yīng)用流體積函數(shù)VOF追蹤自由液面研究了沉沙池中泥沙沉積效果? [13] ,為工程設(shè)計(jì)和運(yùn)行提供了依據(jù)。

    1 工程概況

    那蘇瓦卡里水電站位于尼泊爾東北部地區(qū),背靠喜馬拉雅山脈,地勢北高南低,河谷狹窄、岸坡陡立、階地較為發(fā)育,且河床落差較大。此外由于山體植被覆蓋較少,小型崩塌、滑坡及古滑坡體均較為發(fā)育,水土流失較嚴(yán)重,尤其是在雨季,山體滑坡頻發(fā),河水泥沙含量較大。

    工程從Bhotekoshi河引水,通過4.194 km長的引水隧洞后進(jìn)入地下發(fā)電廠房。發(fā)電水頭達(dá)168.0 m,設(shè)計(jì)發(fā)電流量為80 m? 3 /s,總裝機(jī)容量110 MW。

    整個(gè)樞紐河段河谷狹窄,受限于現(xiàn)場地形地質(zhì)條件,河流兩側(cè)沒有足夠緩坡和灘地布置沉沙池,因此采用地下式沉沙池。

    初擬地下沉沙池為3個(gè)獨(dú)立的洞室,位于河流右岸的山體里面,其上游連接段與取水隧洞相相接。在正常情況下,3個(gè)獨(dú)立的沉沙池連續(xù)運(yùn)行。沉沙池主要尺寸為:上游連接段長度 32.0 m,工作段長度125.0 m,工作寬度15.0 m,下游連接段 24.5 m ;沉沙池工作水深13.0 m,底坡坡率0.02。

    下游連接段布置有出口閘及閘室、排沙閘及閘室、沉沙池操作洞、排沙道等。下游連接段經(jīng)斷面收縮后與沉沙池出水隧洞相接,沉沙池總體平面布置如圖1所示。

    2 沉沙池中沉沙率理論分析

    2.1 泥沙基本資料

    實(shí)測首部樞紐所在河段汛期泥沙濃度為5.624 kg/m? 3 ,泥沙粒徑級(jí)配曲線見圖2。

    2.2 泥沙沉降理論分析

    2.2.1 沉沙池設(shè)計(jì)要求

    根據(jù)工程特點(diǎn)及業(yè)主要求,單個(gè)沉沙池的設(shè)計(jì)應(yīng)滿足以下要求:

    ① 粒徑大于等于0.2 mm的泥沙沉降率應(yīng)大于95%;

    ② 工作段平均最大流速為0.2 m/s;

    ③ 沿程沖洗段的平均流速應(yīng)大于2.5 m/s。

    2.2.2 沉沙率理論分析

    (1) 泥沙沉速計(jì)算? [14] 。

    當(dāng)泥沙粒徑小于等于 0.062 ?mm時(shí),采用斯托克斯公式:

    ω= g 1800 ( ρ? s -ρ? w? ρ? w? ) d 2 γ? (1)

    γ= 0.01775 1+0.0337t+0.000221t 2? (2)

    當(dāng)粒徑為0.062~2.000 mm時(shí),采用沙玉清天然沙沉速公式:

    ( lg S a+3.790)? 2+ ( lg ψ-5.777)? 2=39.0 (3)

    沉速判數(shù):

    S a= ω g? 1/3? ( ρ? s? ρ? w? -1)?? 1/3 γ? 1/3? ?(4)

    粒徑判數(shù):

    ψ= g? 1/3? ( ρ? s? ρ? w? -1)?? 1/3 d 10γ? 2/3? ?(5)

    當(dāng)粒徑大于2.0 mm時(shí),采用沙玉清紊流沉速計(jì)算公式:

    ω=4.58 10d? (6)

    上述式中: γ 為水的運(yùn)動(dòng)黏性系數(shù),cm? 2 /s,20 ℃時(shí)為 0.010 07 ?cm? 2 /s; t 為水溫,℃; d 為泥沙粒徑,mm; ρ?? s 為泥沙密度,取2.65 g/cm? 3 ; ρ?? w 為水密度,取1.0 g/cm? 3 ; ω 為泥沙沉速,cm/s; g 為重力加速度,取9.81 m/s? 2 ; S? a? 為沉速判數(shù); ψ 為粒徑判數(shù)。

    粒徑組平均沉速:

    ω i? = ω iω? i+1? ?(7)

    式中: ω? i? 為粒徑組平均沉速,cm/s; ω? i? 為粒徑組下限粒徑沉速,cm/s; ω? i+1? 為粒徑組上限粒徑沉速, cm/s 。

    (2) 泥沙率計(jì)算。

    各粒徑組懸移質(zhì)沉降率為

    η? ik =1- e?? -α i ?ω — ??iL k q k ??(8)

    q k= Q k B? (9)

    式中: η? ik 為k池段粒徑組沉降率;α? i 為恢復(fù)飽和系數(shù);q? k? 為池段單寬流量,m? 2 /s; Q? k? 為池段平均流量,m? 3 /s,水電站定期沖洗式沉沙池為引用流量,連續(xù)沖洗式沉沙池為引用流量和沖沙流量之和; L? k? 為池段長,m; B 為工作寬度,m,對于矩形斷面沉沙池,為工作段閘室的凈寬度,對于梯形斷面沉沙池,可近似取正常水位以下的平均寬度。

    恢復(fù)飽和系數(shù):

    α i=K (? ω — ??i u * )?? 0.25? (10)

    u *= nv g? R? 1/6? ?(11)

    式中: u?? * 為摩阻流速,m/s; ω — ??i? 為粒徑組平均沉速, cm/s ; K 為綜合經(jīng)驗(yàn)系數(shù),一般取 K =1.2; v 為斷面平均流速,m/s; R 為池段平均過水?dāng)嗝嫠Π霃?,m; n 為池壁糙率; g 為重力加速度,m/s? 2 。

    不小于 i 粒徑級(jí)的沉降率為

    η l= ??l ?i=m? η i Δ P i ??l ?i=m?? Δ P i? (12)

    式中:l為第i粒徑組下限粒徑編號(hào),按粒徑由大到小排列; Δ P? i 為第i粒徑組沙重百分?jǐn)?shù) ,%。經(jīng)計(jì)算,在沉沙池計(jì)算長度為125.0 m時(shí),沉沙池內(nèi)粒徑不小于 0.2 mm 泥沙沉降率為97.8%,滿足設(shè)計(jì)要求。

    (3) 計(jì)算結(jié)果。

    沉沙池沉降率計(jì)算結(jié)果見表1。

    3 數(shù)值模擬

    3.1 控制方程

    水沙兩相流模型基于連續(xù)介質(zhì)假設(shè),將顆粒相作為連續(xù)相處理。兩相均由各自的連續(xù)性方程和動(dòng)量方程描述,兩相相互作用通過相間作用力耦合,本文采用質(zhì)量平均來處理。流經(jīng)沉沙池內(nèi)的流體假設(shè)為不可壓縮的牛頓流體,本次研究的連續(xù)方程和動(dòng)量方程分別為

    ρ ?t + ?ρ U i? ??x i =0 (13)

    ρ U i? ??t + ?ρ U iU j? ??x j =- ?P -? ?x i +

    x j?? μ+ρυ T?? ??U i? ??x j + ??U j? ??x i?? +ρg i (14)

    υ T=C μ k 2 ε? (15)

    式中:t表示時(shí)間, U i? 表示 U? i方向的平均速度i=(1,2,3,…,n),P - 為平均壓力,ρ為流體密度,g? i 為不同方向的重力加速度,k為紊流動(dòng)能,ε為紊流動(dòng)能消散率,μ為水分子黏性系數(shù),υ? T 為紊流運(yùn)動(dòng)黏性系數(shù),C? μ = 0.0085 ,為半經(jīng)驗(yàn)系數(shù)。

    3.2 紊流模型

    本文基于大型流體計(jì)算軟件CFX對沉沙池的含沙水流進(jìn)行數(shù)值模擬。目前連續(xù)相湍流模型中最具代表性的兩方程模型是標(biāo)準(zhǔn) k~ε 模型,然而由于水沙兩相流的不均勻性,水相與泥沙相之間的相互作用又難以精確模擬,本次紊流模型選用以控制方程的時(shí)間平均形式和RNG ?k~ε 湍流模型來實(shí)現(xiàn)雷諾應(yīng)力的封閉。RNG ?k~ε 模型是由一種嚴(yán)謹(jǐn)?shù)慕y(tǒng)計(jì)方法即重整化群理論推導(dǎo)出的,它與Standard ?k~ε 湍流模型保持形式上的一致并有所改進(jìn),在其 ε 方程中增加了額外項(xiàng),這樣顯著提高了快速應(yīng)變流動(dòng)的準(zhǔn)確性。

    3.3 數(shù)學(xué)模型

    沉沙池幾何模型極其復(fù)雜,難點(diǎn)之一是幾何尺度跨度非常大且不規(guī)則,需進(jìn)行必要的簡化。本次研究采用MS來建立沉沙池三維模型。模型簡化如圖3所示。

    模擬高濃度泥沙沉積過程,需要同時(shí)求解泥沙和水的連續(xù)性方程和動(dòng)量方程,因此計(jì)算量很大。網(wǎng)格劃分原則是首先保證網(wǎng)格尺寸均勻過渡,不會(huì)使計(jì)算因網(wǎng)格縱橫比太大而出現(xiàn)流量不守恒。其次,為保證泥沙沉降計(jì)算的穩(wěn)定性,需選擇合適的計(jì)算時(shí)間步長來確定網(wǎng)格尺寸。經(jīng)過多次嘗試,在保證Courant Number最大值在40上下波動(dòng)的情況下,選擇全局網(wǎng)格最大尺寸為1 200 mm,進(jìn)口段最大為650 mm,進(jìn)口截面最大為400 mm,沉沙池底部最大為450 mm,沉沙池尾部最大為650 mm,排沙管最大為350 mm。網(wǎng)格采用正四面體網(wǎng)格,網(wǎng)格單元為 2 035 047 ,網(wǎng)格節(jié)點(diǎn)數(shù)為396 739,最大縱橫比14,最小角度16°,最小行列式值0.8。模型劃分網(wǎng)格節(jié)點(diǎn)如圖4所示。

    3.4 邊界條件

    3.4.1 進(jìn)口邊界條件

    (1) ?入口流速。

    水相和泥沙相均采用相同的入口流速條件。平行于進(jìn)口邊界的速度為零,垂直于進(jìn)口邊界的流速由竇國仁明渠流速度分布公式給定:

    U U *? ????inlet? =2.5 ln? 1+y * +7.05? ?y * 1+y *? ??2+

    2.5? y * 1+y *? -B *+0.55 1- cos? 2 π? y H??? (16)

    式中:y為垂向高度;U? * 為摩阻流速;B? * 為壁面粗糙度對流速分布的影響,y *=ρU *y/sμ,本次計(jì)算分析中壁面糙率在粗糙區(qū),B? * 取值為2.3,H為明渠流水深。

    (2) 入口體積分?jǐn)?shù)。

    進(jìn)口邊界水相體積分?jǐn)?shù)為

    r? w????? inlet? =1-? r? s????? inlet?? (17)

    r? s = S? w? ρ? s? ?(18)

    式中: r?? w? 為水相體積分?jǐn)?shù),r?? s? 為泥沙相體積分?jǐn)?shù),S?? w? 為泥沙濃度,ρ?? s? 為泥沙密度。本次計(jì)算分析中取進(jìn)口流量Q =36.67 m? 3 /s,因此泥沙相體積分?jǐn)?shù) r?? s =5.6/2650×100%=0.2113%

    本次研究主要選取0.20,0.25,0.50,1.00 mm 4種顆粒級(jí)配泥沙粒徑,不同粒徑泥沙的體積分?jǐn)?shù)如表2所列。

    3.4.2 出口邊界條件

    (1) 出口壓力。

    水沙兩相出口平均靜壓均為零。

    (2) 出口湍流。

    因流動(dòng)為充分發(fā)展的湍流,出口給定湍動(dòng)能和耗散率零梯度邊界條件,即湍動(dòng)能和耗散率沿出口方向一階導(dǎo)數(shù)為零。

    (3) 出口濃度。

    沉沙池應(yīng)滿足平衡輸沙條件,出口濃度邊界條件給定垂直于出口邊界的兩相體積分?jǐn)?shù)方向?qū)?shù)為零。

    3.4.3 自由表面邊界

    因數(shù)學(xué)模型考慮重力,本次研究采用VOF體積分辨率函數(shù)來追蹤自由表面。

    3.4.4 壁面邊界條件

    壁面采用無滑移邊界條件,即壁面上兩相速度均為零。

    3.5 結(jié)果分析

    3.5.1 泥沙沉積率

    本次研究采用多組分泥沙顆粒模型對沉沙池沉沙效果進(jìn)行分析計(jì)算,由于很難單純直接對多組分泥沙顆粒沉沙率進(jìn)行查看,故采用含沙水流中泥沙體積分?jǐn)?shù)變化率來反映沉沙池中泥沙沉降效果。某顆粒粒徑的泥沙沉降率用 η 表示:

    η= 1- V?? outlet?? V?? inlet?? ?×100 %? (19)

    式中:V?? outlet? 表示某粒徑出口泥沙體積分?jǐn)?shù);V?? inlet? 表示某粒徑進(jìn)口泥沙體積分?jǐn)?shù)。

    3.5.2 沉沙池不同顆粒直徑泥沙沉積效果

    為了達(dá)到更好的收斂效果,本次計(jì)算時(shí)間步長取為0.5 s,在650~800 s時(shí)基本達(dá)到動(dòng)態(tài)平衡穩(wěn)定狀態(tài)。為更精細(xì)地了解泥沙的沉降情況,以沉沙池出水口截面的泥沙平均濃度的體積分?jǐn)?shù)為監(jiān)測點(diǎn),分析不同粒徑泥沙體積分?jǐn)?shù)隨時(shí)間的變化,其變化趨勢如圖5所示。

    根據(jù)出水口截面泥沙體積分?jǐn)?shù)在不同時(shí)刻變化情況計(jì)算得出泥沙的沉降率,不同顆粒直徑泥沙隨時(shí)間變化的沉降率如表3所列。

    計(jì)算成果表明:沉沙池內(nèi)水流含沙量沿垂線分布是上小下大,即水流含沙濃度接近池底最高,在水面含沙量最低,符合泥沙的運(yùn)動(dòng)擴(kuò)散規(guī)律;且通過以上計(jì)算對比可知,不同粒徑泥沙的沉沙率均滿足設(shè)計(jì)要求。相關(guān)研究表明:混合泥沙中粗顆粒對細(xì)顆粒有一定的遮蔽效應(yīng),即與等粒徑的均勻泥沙顆粒相比,混合泥沙中同等粒徑的細(xì)顆粒由于受到粗顆粒的制約而沉降更快。沉沙池進(jìn)口區(qū)域受水流紊動(dòng)影響,有輕微的絮凝現(xiàn)象發(fā)生,符合泥沙沉淀一般規(guī)律。

    4 結(jié) 論

    (1) 基于傳統(tǒng)理論分析,結(jié)合國內(nèi)外研究學(xué)者總結(jié)得出的沉沙池中含沙水流泥沙顆粒沉降理論,推導(dǎo) 出該水電站沉沙池不同顆粒直徑泥砂的沉沙率,其中不小于0.2 mm顆粒直徑泥沙沉降率為97.8%,滿足工程設(shè)計(jì)要求。

    (2) 數(shù)值模擬計(jì)算結(jié)果表明:該工程不小于0.2 mm粒徑泥沙沉積率為95.59%,與理論分析結(jié)果較為吻合。而且數(shù)值方法還能模擬沉沙池沖砂效果,比理論分析方法更直觀。

    參考文獻(xiàn):

    [1]? 袁晶,許全喜,童輝.三峽水庫蓄水運(yùn)用以來庫區(qū)泥沙淤積特性研究[J].水力發(fā)電學(xué)報(bào),2013,32(2):139-145.

    [2] 王普慶.小浪底水庫運(yùn)用初期庫區(qū)泥沙淤積分布特征[J].人民黃河,2012,34(10):26-28.

    [3] 廖庭庭,陳和春,高甜,等.三峽水電站過機(jī)泥沙粒徑對水輪機(jī)葉片空化空蝕的影響[J].中國農(nóng)村水利水電,2012(2):121-124.

    [4] DUBOYS M P.One dimensional bed wave movement in lowland rivers[J].Water Resources Research,1983,19(3):627-631.

    [5] ROUSE H.Experiments on the mechanics of sediment suspension[C]∥Proceedings of the 5th International Congress for Applied Mechanics,Cambridge,1939:550-554.

    [6] 竇國仁.論河流紊動(dòng)與流速分布[J].水利學(xué)報(bào),1959,5(4):49-50.

    [7] 蔡樹棠.泥沙對水的相對速度和含沙濃度對脈動(dòng)速度的影響[J].應(yīng)用數(shù)學(xué)和力學(xué),1982,3(1):1-10.

    [8] 楊紅,潘光在.沉沙池立面流場數(shù)值模擬研究[J].泥沙研究,2000,18(2):55-59.

    [9] 胡斌,桂林.斜板式沉沙池水流流態(tài)數(shù)值模擬[J].西南民族大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,38(3):425-430.

    [10]? 周曉泉.復(fù)雜邊界條件下的三維紊流數(shù)值模擬研宄[D].成都:四川大學(xué),2003.

    [11] 張軍,侍克斌,高亞平,等.圓中環(huán)沉沙排沙池流場三維數(shù)值模擬研究[J].水力發(fā)電學(xué)報(bào),2015,34(3):103-110.

    [12] 岳竹,肖鈺.有壓隧洞內(nèi)沉沙池泥沙沉降特性研究:以尼泊爾尼亞迪水電站為例[J].人民長江,2020,51(增1):200-203.

    [13] CHEN Q,DAI G,LIU H.Volume of fluid model for turbulence numerical simulation of stepped spillway overflow[J].Journal of Hydraulic Engineering,2002,128:683-8.

    [14] 朱鑒遠(yuǎn).水利水電工程泥沙設(shè)計(jì)[M].北京:中國水利水電出版社,2010.

    (編輯:胡旭東)

    Study on sand depositing and sand flushing effect in grit chamber of ?Rasuwagadhi Hydropower Station

    GAO Donghong,ZOU Hao

    ( Shanghai Investigation,Design & Research Institute Co.,ltd.,Shanghai 200434,China )

    Abstract:

    In water conservancy and hydropower projects,a grit chamber is generally used to reduce the sediment content,so as to reduce the adverse impact of the sediment on the water diversion and power generation system.It is very important to study the flow pattern of sand-containing water flow,sand depositing and sand flushing effect in the sand grit chamber.Based on a grit chamber of Rasuwagadhi Hydropower Station,using the theoretical formula,RNG ?k~ε ?turbulence model and volume resolution function VOF,we simulated the flow pattern of the two-phase flow of sediment in the grit chamber,and showed the deposition effect of different sand particle in the multi-diameter particle composition in detail.The results showed that the sand depositing efficiency of no less than 0.2 mm was more than 95%,and the numerical results agreed with the theoretical calculation and met the engineering design requirements.

    Key words:

    grit chamber;sand depositing efficiency;numerical simulation;Rasuwagadhi Hydropower Station;Nepal

    猜你喜歡
    沉沙池泥沙粒徑
    泥沙做的父親
    木屑粒徑對黑木耳栽培的影響試驗(yàn)*
    新疆多泥沙河流水庫泥沙處理措施
    土壤團(tuán)聚體對泥沙沉降速度的影響
    滲透破壞沖淤式沉沙池技術(shù)研究
    基于近場散射的顆粒粒徑分布測量
    引黃濟(jì)青沉沙池運(yùn)行存在問題及對策
    山東水利(2017年5期)2017-04-07 06:29:51
    自排式沉沙池上游渠道水力學(xué)特性研究
    泥沙滅火
    兒童繪本(2015年2期)2015-05-25 18:10:15
    Oslo結(jié)晶器晶體粒徑分布特征的CFD模擬
    尾随美女入室| 免费观看精品视频网站| 国产欧美日韩精品一区二区| 久久久久久久久中文| 日韩国内少妇激情av| 中文在线观看免费www的网站| 精品一区二区三区视频在线观看免费| 成人精品一区二区免费| 欧美日本视频| 免费无遮挡裸体视频| 成人特级黄色片久久久久久久| 亚洲国产高清在线一区二区三| 日本免费一区二区三区高清不卡| 美女免费视频网站| 在线免费观看的www视频| 色哟哟·www| 我要看日韩黄色一级片| 岛国在线免费视频观看| 国模一区二区三区四区视频| 国产淫片久久久久久久久| 精品午夜福利视频在线观看一区| 日本一本二区三区精品| 中文资源天堂在线| 在线a可以看的网站| 久久欧美精品欧美久久欧美| 色吧在线观看| 男人和女人高潮做爰伦理| 级片在线观看| 亚洲一区二区三区色噜噜| 亚洲一级一片aⅴ在线观看| 六月丁香七月| 免费黄网站久久成人精品| 日韩国内少妇激情av| 最好的美女福利视频网| 日日干狠狠操夜夜爽| 男人狂女人下面高潮的视频| 白带黄色成豆腐渣| 一a级毛片在线观看| 综合色丁香网| 精品久久久久久成人av| 日本成人三级电影网站| 国产精品国产三级国产av玫瑰| videossex国产| 高清毛片免费看| 黄色一级大片看看| 色综合亚洲欧美另类图片| 一个人观看的视频www高清免费观看| 色视频www国产| 中国美白少妇内射xxxbb| av中文乱码字幕在线| 嫩草影院精品99| 99热全是精品| 日本熟妇午夜| 日韩欧美一区二区三区在线观看| 搡老妇女老女人老熟妇| 禁无遮挡网站| 成人av在线播放网站| 免费看日本二区| 国产av一区在线观看免费| 久久鲁丝午夜福利片| 伊人久久精品亚洲午夜| 国产熟女欧美一区二区| 欧美又色又爽又黄视频| 国产私拍福利视频在线观看| 欧美一区二区亚洲| 亚洲欧美日韩高清专用| 成年女人毛片免费观看观看9| 日韩大尺度精品在线看网址| 日韩成人伦理影院| 国内少妇人妻偷人精品xxx网站| 少妇的逼好多水| 乱系列少妇在线播放| 一级毛片aaaaaa免费看小| 国产精品综合久久久久久久免费| 狂野欧美激情性xxxx在线观看| 一级黄片播放器| 日韩欧美在线乱码| 搡老岳熟女国产| 一区二区三区四区激情视频 | 久久精品国产99精品国产亚洲性色| 尾随美女入室| 亚洲精品日韩av片在线观看| 午夜亚洲福利在线播放| 欧美丝袜亚洲另类| 日本免费a在线| 亚洲欧美成人综合另类久久久 | 简卡轻食公司| 亚洲国产精品国产精品| 国产色婷婷99| 免费观看在线日韩| 久久久久久久久大av| 99九九线精品视频在线观看视频| 最后的刺客免费高清国语| 麻豆精品久久久久久蜜桃| 亚洲国产精品国产精品| 少妇的逼水好多| 亚洲人成网站在线播放欧美日韩| 午夜福利18| 久久久久久久久久黄片| 超碰av人人做人人爽久久| 日韩亚洲欧美综合| 亚洲精华国产精华液的使用体验 | 噜噜噜噜噜久久久久久91| 国产亚洲欧美98| 久久精品国产亚洲av涩爱 | 69人妻影院| 老司机福利观看| 成年版毛片免费区| 国产成年人精品一区二区| 国产欧美日韩一区二区精品| av在线老鸭窝| 国产精品久久久久久久电影| 国产av不卡久久| 久99久视频精品免费| 成人精品一区二区免费| 日本 av在线| 俺也久久电影网| 中文字幕久久专区| 2021天堂中文幕一二区在线观| 日本a在线网址| 国产91av在线免费观看| 成年女人看的毛片在线观看| 九色成人免费人妻av| 少妇人妻精品综合一区二区 | 欧美日韩综合久久久久久| av免费在线看不卡| 久久精品国产亚洲网站| 在线观看66精品国产| 长腿黑丝高跟| 国产精品伦人一区二区| 在线天堂最新版资源| 美女黄网站色视频| 可以在线观看的亚洲视频| 欧美日本视频| 久久中文看片网| 国产黄色视频一区二区在线观看 | 日韩欧美一区二区三区在线观看| 我要搜黄色片| 精品无人区乱码1区二区| 国产在视频线在精品| 日韩成人伦理影院| 一级黄色大片毛片| 免费看a级黄色片| 久久精品综合一区二区三区| 久久九九热精品免费| 内地一区二区视频在线| 亚洲性久久影院| 在线观看av片永久免费下载| 国产精品国产三级国产av玫瑰| 夜夜爽天天搞| 欧美成人免费av一区二区三区| 精品久久久久久久人妻蜜臀av| 亚洲国产高清在线一区二区三| 草草在线视频免费看| 欧美日韩乱码在线| 99热这里只有是精品在线观看| 中文字幕av在线有码专区| 身体一侧抽搐| 成人国产麻豆网| 日本免费a在线| 99热全是精品| 国产伦精品一区二区三区四那| 国产成人精品久久久久久| 亚洲欧美成人综合另类久久久 | 午夜福利18| 色综合站精品国产| 亚洲欧美清纯卡通| 日本三级黄在线观看| 成人高潮视频无遮挡免费网站| av中文乱码字幕在线| 女生性感内裤真人,穿戴方法视频| 欧美色视频一区免费| 国产大屁股一区二区在线视频| 亚洲av五月六月丁香网| 亚洲欧美日韩无卡精品| 国产精品精品国产色婷婷| 人妻少妇偷人精品九色| 亚洲av成人精品一区久久| .国产精品久久| 婷婷六月久久综合丁香| 在线观看美女被高潮喷水网站| 中文字幕免费在线视频6| 免费观看在线日韩| 一区二区三区四区激情视频 | 成人av一区二区三区在线看| 嫩草影院新地址| 日韩欧美国产在线观看| 91av网一区二区| 久久九九热精品免费| 欧美日韩国产亚洲二区| 欧美精品国产亚洲| 99热6这里只有精品| 中文字幕av在线有码专区| 最新中文字幕久久久久| 非洲黑人性xxxx精品又粗又长| 网址你懂的国产日韩在线| 国产亚洲91精品色在线| 99热这里只有精品一区| 久久久午夜欧美精品| 久久韩国三级中文字幕| 国产精品电影一区二区三区| 亚洲成人中文字幕在线播放| 好男人在线观看高清免费视频| 国产精品不卡视频一区二区| www.色视频.com| 欧美xxxx性猛交bbbb| 老熟妇乱子伦视频在线观看| 深爱激情五月婷婷| 欧美潮喷喷水| 国产成人精品久久久久久| 日韩欧美精品v在线| 亚洲av免费在线观看| 久久精品国产清高在天天线| 国产成人a∨麻豆精品| 99久国产av精品| 国产精品一区二区三区四区免费观看 | 精品福利观看| 日韩欧美一区二区三区在线观看| 国产精品久久视频播放| 99久久成人亚洲精品观看| 小说图片视频综合网站| 男女做爰动态图高潮gif福利片| 日韩av在线大香蕉| 午夜精品国产一区二区电影 | 97超级碰碰碰精品色视频在线观看| 久久精品国产亚洲av天美| 男女啪啪激烈高潮av片| 给我免费播放毛片高清在线观看| 天天一区二区日本电影三级| 天堂动漫精品| 久久久久久久久久久丰满| 欧美又色又爽又黄视频| 在线免费十八禁| 99久久九九国产精品国产免费| 国产一区二区激情短视频| 日本在线视频免费播放| 欧美日韩精品成人综合77777| 国产真实乱freesex| 亚洲色图av天堂| 中文字幕熟女人妻在线| 99热精品在线国产| 麻豆久久精品国产亚洲av| 狠狠狠狠99中文字幕| 亚洲三级黄色毛片| 久久婷婷人人爽人人干人人爱| 免费在线观看成人毛片| 成人av一区二区三区在线看| 亚洲婷婷狠狠爱综合网| 国产精品乱码一区二三区的特点| 有码 亚洲区| 在线播放国产精品三级| 婷婷色综合大香蕉| 夜夜夜夜夜久久久久| 日本-黄色视频高清免费观看| 久久久精品欧美日韩精品| 男女那种视频在线观看| 草草在线视频免费看| 亚洲成a人片在线一区二区| 精品熟女少妇av免费看| 91久久精品国产一区二区成人| av在线老鸭窝| 欧美又色又爽又黄视频| 一区二区三区四区激情视频 | 亚洲天堂国产精品一区在线| 午夜福利在线观看免费完整高清在 | 男女啪啪激烈高潮av片| 乱人视频在线观看| 国产精品久久久久久久久免| 国产黄片美女视频| www.色视频.com| 精品欧美国产一区二区三| 在线国产一区二区在线| 午夜影院日韩av| 日日摸夜夜添夜夜爱| 精品久久久久久久久久久久久| ponron亚洲| 一级毛片久久久久久久久女| 精品久久久噜噜| 亚洲美女视频黄频| 免费看日本二区| 日韩一本色道免费dvd| 一区二区三区免费毛片| 国产精品亚洲美女久久久| 变态另类丝袜制服| 人妻丰满熟妇av一区二区三区| 男女下面进入的视频免费午夜| 99热全是精品| 一进一出抽搐gif免费好疼| 亚洲婷婷狠狠爱综合网| 国产白丝娇喘喷水9色精品| 国产精品爽爽va在线观看网站| 国产色爽女视频免费观看| 最新中文字幕久久久久| 成人亚洲精品av一区二区| av在线观看视频网站免费| 国产欧美日韩一区二区精品| 99视频精品全部免费 在线| 久久6这里有精品| 搡老熟女国产l中国老女人| 麻豆国产97在线/欧美| 国产大屁股一区二区在线视频| 五月伊人婷婷丁香| 免费看日本二区| 日产精品乱码卡一卡2卡三| 欧美色视频一区免费| 国产精品亚洲一级av第二区| 91在线精品国自产拍蜜月| 亚洲人成网站在线观看播放| 亚洲国产欧洲综合997久久,| 一本精品99久久精品77| 亚洲综合色惰| 天堂影院成人在线观看| 一区福利在线观看| 亚洲人成网站在线播| 亚洲欧美日韩无卡精品| 男人狂女人下面高潮的视频| 深夜a级毛片| 麻豆乱淫一区二区| 亚洲真实伦在线观看| 午夜爱爱视频在线播放| 国产爱豆传媒在线观看| 国产美女午夜福利| 欧美+亚洲+日韩+国产| 日本欧美国产在线视频| 日本黄色视频三级网站网址| 一区二区三区免费毛片| 国产老妇女一区| 嫩草影院入口| h日本视频在线播放| 天美传媒精品一区二区| 久久精品国产亚洲av涩爱 | 亚洲一区高清亚洲精品| 国产三级在线视频| 最新在线观看一区二区三区| 亚洲av中文av极速乱| 乱系列少妇在线播放| 人妻夜夜爽99麻豆av| 国产精华一区二区三区| 国产午夜福利久久久久久| 亚洲va在线va天堂va国产| av福利片在线观看| 18禁裸乳无遮挡免费网站照片| 色综合色国产| av卡一久久| 老熟妇乱子伦视频在线观看| 久久6这里有精品| 国产精品乱码一区二三区的特点| 欧美日韩乱码在线| 久久鲁丝午夜福利片| 男女之事视频高清在线观看| 在线国产一区二区在线| av免费在线看不卡| 99热全是精品| 看十八女毛片水多多多| 免费高清视频大片| 欧美另类亚洲清纯唯美| 欧美+亚洲+日韩+国产| 午夜免费激情av| 欧美人与善性xxx| 亚洲av成人av| 桃色一区二区三区在线观看| 小蜜桃在线观看免费完整版高清| 18禁黄网站禁片免费观看直播| 成人av在线播放网站| 亚洲av五月六月丁香网| 中文字幕熟女人妻在线| 国产成人a区在线观看| 国产乱人视频| 99久国产av精品国产电影| 黄色欧美视频在线观看| 国产一区二区在线观看日韩| 国产美女午夜福利| 亚洲欧美日韩东京热| 波多野结衣高清作品| 老司机影院成人| 国产探花极品一区二区| 国产国拍精品亚洲av在线观看| 久久精品国产99精品国产亚洲性色| 熟女电影av网| 美女xxoo啪啪120秒动态图| 精品久久久久久久久久久久久| av在线观看视频网站免费| 日韩欧美精品免费久久| 在线播放无遮挡| 搡女人真爽免费视频火全软件 | 日韩一区二区视频免费看| avwww免费| 午夜激情欧美在线| 亚洲成人av在线免费| 免费看美女性在线毛片视频| 国产大屁股一区二区在线视频| 亚洲美女黄片视频| 亚洲精品国产av成人精品 | 69人妻影院| 99久久精品一区二区三区| 99久久精品热视频| 久久精品国产亚洲网站| 1000部很黄的大片| 国产精品久久视频播放| 国产伦在线观看视频一区| 亚洲av免费在线观看| 在线播放国产精品三级| 3wmmmm亚洲av在线观看| 中出人妻视频一区二区| 夜夜爽天天搞| 亚洲av五月六月丁香网| 亚洲欧美清纯卡通| 日本免费一区二区三区高清不卡| 成年免费大片在线观看| 男女下面进入的视频免费午夜| 日本在线视频免费播放| 国产成年人精品一区二区| 99热精品在线国产| 亚洲性夜色夜夜综合| 精品少妇黑人巨大在线播放 | 亚洲国产欧美人成| 亚洲无线观看免费| 亚洲av熟女| 人人妻人人澡人人爽人人夜夜 | 国产白丝娇喘喷水9色精品| 欧美成人一区二区免费高清观看| 你懂的网址亚洲精品在线观看 | 久久久久国内视频| 12—13女人毛片做爰片一| 少妇人妻精品综合一区二区 | 日韩欧美三级三区| 真人做人爱边吃奶动态| 俄罗斯特黄特色一大片| 伦理电影大哥的女人| 成人无遮挡网站| 寂寞人妻少妇视频99o| 成人av在线播放网站| 人人妻,人人澡人人爽秒播| 美女高潮的动态| 在线观看一区二区三区| 真实男女啪啪啪动态图| 男女边吃奶边做爰视频| 在线播放无遮挡| 狠狠狠狠99中文字幕| 变态另类成人亚洲欧美熟女| 国产精品嫩草影院av在线观看| 精品人妻偷拍中文字幕| 女人被狂操c到高潮| 精品国产三级普通话版| 亚洲欧美日韩高清在线视频| 丰满乱子伦码专区| 久久久久国产精品人妻aⅴ院| 97热精品久久久久久| 国产老妇女一区| 亚洲国产精品成人久久小说 | 午夜亚洲福利在线播放| 日本成人三级电影网站| 一进一出抽搐动态| 一级a爱片免费观看的视频| 别揉我奶头 嗯啊视频| 美女黄网站色视频| 在线国产一区二区在线| 日韩av在线大香蕉| 身体一侧抽搐| 成人鲁丝片一二三区免费| 国产精品av视频在线免费观看| 日韩欧美国产在线观看| 成人特级av手机在线观看| 国产亚洲精品av在线| 日韩av在线大香蕉| 欧美日韩一区二区视频在线观看视频在线 | 免费搜索国产男女视频| av在线蜜桃| 国产精品一二三区在线看| 少妇被粗大猛烈的视频| 精品久久久噜噜| 黄色视频,在线免费观看| 少妇熟女aⅴ在线视频| 午夜老司机福利剧场| 悠悠久久av| 干丝袜人妻中文字幕| www.色视频.com| 成人欧美大片| 最近中文字幕高清免费大全6| 九九在线视频观看精品| 精品无人区乱码1区二区| 久久精品久久久久久噜噜老黄 | 日本在线视频免费播放| 国产一区二区在线av高清观看| 国产精品亚洲美女久久久| 中国美女看黄片| 成人毛片a级毛片在线播放| 99热全是精品| 国产蜜桃级精品一区二区三区| 啦啦啦韩国在线观看视频| 观看美女的网站| av天堂中文字幕网| 三级毛片av免费| 欧美bdsm另类| 欧美最黄视频在线播放免费| 91在线精品国自产拍蜜月| 成人特级黄色片久久久久久久| 99热6这里只有精品| 日本 av在线| 啦啦啦韩国在线观看视频| 日本爱情动作片www.在线观看 | 久久99热6这里只有精品| 国产伦精品一区二区三区视频9| 最近最新中文字幕大全电影3| 狂野欧美激情性xxxx在线观看| 免费无遮挡裸体视频| 亚洲国产色片| 日韩中字成人| 联通29元200g的流量卡| 校园人妻丝袜中文字幕| 免费观看在线日韩| 啦啦啦啦在线视频资源| 99热只有精品国产| 亚洲无线观看免费| 美女cb高潮喷水在线观看| 亚洲精品一区av在线观看| 尤物成人国产欧美一区二区三区| a级毛片a级免费在线| 免费无遮挡裸体视频| 能在线免费观看的黄片| 成人国产麻豆网| 国产精品野战在线观看| 国产成人a∨麻豆精品| 精品人妻偷拍中文字幕| 高清毛片免费看| 97在线视频观看| 国语自产精品视频在线第100页| 亚洲国产色片| 午夜a级毛片| 国产精品野战在线观看| av卡一久久| 亚洲精品亚洲一区二区| 国产精品永久免费网站| 欧美一区二区国产精品久久精品| av福利片在线观看| 欧美色欧美亚洲另类二区| 在线观看一区二区三区| 亚洲图色成人| 国产精品免费一区二区三区在线| 直男gayav资源| 高清日韩中文字幕在线| 国产av一区在线观看免费| 最新中文字幕久久久久| 免费黄网站久久成人精品| www.色视频.com| 99热这里只有是精品50| 天堂√8在线中文| 国产在线男女| 综合色av麻豆| 搞女人的毛片| 一个人观看的视频www高清免费观看| 精品福利观看| 一本久久中文字幕| 我的老师免费观看完整版| 一进一出抽搐动态| 亚洲真实伦在线观看| 极品教师在线视频| 久久精品夜夜夜夜夜久久蜜豆| 国产精品永久免费网站| 97超视频在线观看视频| 国产精品久久久久久久电影| 久久九九热精品免费| 波多野结衣高清无吗| 久久久久精品国产欧美久久久| 亚洲最大成人av| 全区人妻精品视频| 国产v大片淫在线免费观看| 变态另类丝袜制服| 波野结衣二区三区在线| 欧美色欧美亚洲另类二区| 亚洲av美国av| 看免费成人av毛片| 狂野欧美白嫩少妇大欣赏| 久久久久久伊人网av| 欧美日本视频| 男人和女人高潮做爰伦理| 在线观看一区二区三区| 18+在线观看网站| 美女cb高潮喷水在线观看| 18禁黄网站禁片免费观看直播| 国产精品野战在线观看| 亚洲精品久久国产高清桃花| 亚洲欧美日韩无卡精品| 美女被艹到高潮喷水动态| 欧美xxxx黑人xx丫x性爽| 村上凉子中文字幕在线| 精品一区二区三区av网在线观看| 日韩av在线大香蕉| 日本五十路高清| 亚洲在线观看片| 毛片一级片免费看久久久久| 国产男靠女视频免费网站| 亚洲av成人精品一区久久| 91久久精品国产一区二区成人| 久久久午夜欧美精品| 亚洲在线自拍视频| 欧美zozozo另类| 国产一区二区亚洲精品在线观看| 欧美最新免费一区二区三区| 麻豆av噜噜一区二区三区| 少妇的逼好多水| 小说图片视频综合网站| 一级av片app| 男女边吃奶边做爰视频| 插逼视频在线观看| 91精品国产九色| 91久久精品国产一区二区三区| 白带黄色成豆腐渣| 久久99热6这里只有精品| АⅤ资源中文在线天堂| 国产高清不卡午夜福利| 国产精品伦人一区二区| 丰满乱子伦码专区| 国产一级毛片七仙女欲春2| 亚洲精品乱码久久久v下载方式|