• <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模擬
    精品福利观看| 一个人免费在线观看的高清视频| 一级a爱视频在线免费观看| 日韩精品中文字幕看吧| 亚洲少妇的诱惑av| 深夜精品福利| 99久久人妻综合| 一进一出好大好爽视频| 日本免费a在线| а√天堂www在线а√下载| 亚洲av第一区精品v没综合| 国产精品av久久久久免费| 久久精品影院6| 他把我摸到了高潮在线观看| 亚洲国产精品999在线| 亚洲伊人色综图| 久久人人97超碰香蕉20202| 在线观看66精品国产| 国产精品国产高清国产av| 在线十欧美十亚洲十日本专区| 一级毛片女人18水好多| a在线观看视频网站| 欧洲精品卡2卡3卡4卡5卡区| 国产成人精品无人区| 色综合站精品国产| 欧美黑人精品巨大| 亚洲一区中文字幕在线| 久久午夜综合久久蜜桃| 高清在线国产一区| 中文字幕最新亚洲高清| 中文字幕最新亚洲高清| 超碰成人久久| 国产av在哪里看| 亚洲精品粉嫩美女一区| 久久久国产成人免费| 9色porny在线观看| 国产无遮挡羞羞视频在线观看| 久热这里只有精品99| 日韩欧美免费精品| 精品久久蜜臀av无| 欧美国产精品va在线观看不卡| 国产深夜福利视频在线观看| 99精国产麻豆久久婷婷| 女警被强在线播放| 韩国av一区二区三区四区| 国产欧美日韩一区二区三区在线| 国产一区二区在线av高清观看| 久久九九热精品免费| 国产精品永久免费网站| 亚洲欧美日韩另类电影网站| 亚洲精品久久成人aⅴ小说| 大陆偷拍与自拍| 欧美色视频一区免费| 级片在线观看| 91大片在线观看| 91大片在线观看| 可以在线观看毛片的网站| 午夜影院日韩av| 咕卡用的链子| 韩国av一区二区三区四区| 韩国av一区二区三区四区| 国产熟女xx| 91麻豆av在线| 超碰97精品在线观看| 久久 成人 亚洲| xxx96com| 日韩一卡2卡3卡4卡2021年| 午夜免费激情av| 叶爱在线成人免费视频播放| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲性夜色夜夜综合| 免费高清视频大片| 如日韩欧美国产精品一区二区三区| 国产精品永久免费网站| 亚洲精品中文字幕在线视频| 成人免费观看视频高清| 免费在线观看亚洲国产| 99国产精品一区二区三区| 男人舔女人的私密视频| 丰满的人妻完整版| 91精品国产国语对白视频| 国产一区二区在线av高清观看| 国产在线观看jvid| 波多野结衣av一区二区av| av欧美777| 欧美老熟妇乱子伦牲交| 国产亚洲欧美在线一区二区| 国产国语露脸激情在线看| 久久 成人 亚洲| 精品电影一区二区在线| 嫩草影院精品99| 日本vs欧美在线观看视频| 男男h啪啪无遮挡| 一进一出抽搐gif免费好疼 | 亚洲精品国产色婷婷电影| 亚洲色图av天堂| 国产不卡一卡二| 国产亚洲精品久久久久5区| 久9热在线精品视频| 交换朋友夫妻互换小说| 色在线成人网| 啦啦啦免费观看视频1| 好男人电影高清在线观看| 99久久久亚洲精品蜜臀av| av视频免费观看在线观看| 亚洲狠狠婷婷综合久久图片| 欧美av亚洲av综合av国产av| 免费看a级黄色片| 一区福利在线观看| 精品国产乱码久久久久久男人| 亚洲五月婷婷丁香| 精品一区二区三区四区五区乱码| 精品一区二区三区av网在线观看| 老司机在亚洲福利影院| 久久久久亚洲av毛片大全| 免费av毛片视频| 欧美激情高清一区二区三区| 欧美色视频一区免费| 一区二区日韩欧美中文字幕| 18禁观看日本| 国产成人系列免费观看| 日韩欧美一区视频在线观看| 国产成人系列免费观看| 亚洲av美国av| 欧美激情高清一区二区三区| 视频在线观看一区二区三区| 精品熟女少妇八av免费久了| 香蕉久久夜色| 午夜福利一区二区在线看| 中国美女看黄片| 久久国产精品人妻蜜桃| 久久人人97超碰香蕉20202| 日韩视频一区二区在线观看| 日韩欧美在线二视频| 在线观看免费视频日本深夜| 长腿黑丝高跟| 国产精品亚洲av一区麻豆| 老司机亚洲免费影院| 精品一区二区三卡| 香蕉国产在线看| 免费在线观看日本一区| 成人三级黄色视频| 日韩av在线大香蕉| 午夜成年电影在线免费观看| 老熟妇仑乱视频hdxx| 久久精品人人爽人人爽视色| 欧美中文综合在线视频| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利一区二区在线看| 桃红色精品国产亚洲av| 国产精品一区二区免费欧美| 久久久国产欧美日韩av| 久久久久国内视频| 丝袜美足系列| 丰满的人妻完整版| 国产成人av教育| 亚洲欧美一区二区三区久久| 亚洲第一av免费看| 男人操女人黄网站| 国产免费男女视频| 午夜福利,免费看| 操出白浆在线播放| 岛国在线观看网站| 88av欧美| 在线观看免费午夜福利视频| 黄色a级毛片大全视频| 热re99久久精品国产66热6| 欧美激情极品国产一区二区三区| 91麻豆av在线| 一本大道久久a久久精品| 久久久久亚洲av毛片大全| 老汉色av国产亚洲站长工具| 日韩人妻精品一区2区三区| 久久人人爽av亚洲精品天堂| 亚洲精品一区av在线观看| 久久人妻熟女aⅴ| 黑人猛操日本美女一级片| 不卡一级毛片| 天堂中文最新版在线下载| 日本三级黄在线观看| 亚洲一区二区三区色噜噜 | 日本一区二区免费在线视频| 激情在线观看视频在线高清| 亚洲av成人av| 首页视频小说图片口味搜索| 999久久久精品免费观看国产| 视频区欧美日本亚洲| 大码成人一级视频| 亚洲国产欧美日韩在线播放| 最近最新中文字幕大全免费视频| 人人澡人人妻人| 12—13女人毛片做爰片一| 嫩草影院精品99| 美女扒开内裤让男人捅视频| 日韩精品青青久久久久久| 亚洲七黄色美女视频| 国产高清国产精品国产三级| 女人高潮潮喷娇喘18禁视频| 麻豆av在线久日| 国产99久久九九免费精品| 男女午夜视频在线观看| 两性夫妻黄色片| 亚洲欧美一区二区三区久久| 国产高清激情床上av| 一夜夜www| 亚洲人成电影观看| 日本精品一区二区三区蜜桃| 久久香蕉国产精品| 免费女性裸体啪啪无遮挡网站| 久久精品成人免费网站| 国产成人一区二区三区免费视频网站| 久久伊人香网站| 麻豆一二三区av精品| 中文亚洲av片在线观看爽| 国产一区二区激情短视频| 91麻豆av在线| 嫩草影视91久久| 亚洲黑人精品在线| 色综合婷婷激情| 国产xxxxx性猛交| 国产精华一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 巨乳人妻的诱惑在线观看| 中文欧美无线码| 亚洲一区二区三区不卡视频| 久久精品国产清高在天天线| 成年女人毛片免费观看观看9| 侵犯人妻中文字幕一二三四区| 国产午夜精品久久久久久| 欧美成人免费av一区二区三区| 一区二区三区激情视频| 高清黄色对白视频在线免费看| 怎么达到女性高潮| av视频免费观看在线观看| 人人妻人人澡人人看| 老司机亚洲免费影院| 亚洲精品av麻豆狂野| 亚洲成人精品中文字幕电影 | 一本综合久久免费| 国产欧美日韩精品亚洲av| 可以免费在线观看a视频的电影网站| 久久香蕉精品热| av在线天堂中文字幕 | 麻豆久久精品国产亚洲av | 又紧又爽又黄一区二区| 欧美日韩中文字幕国产精品一区二区三区 | 久久精品成人免费网站| 亚洲午夜理论影院| 亚洲av成人av| av天堂久久9| 亚洲少妇的诱惑av| 欧美日韩一级在线毛片| 日韩中文字幕欧美一区二区| 在线看a的网站| 免费观看精品视频网站| 1024香蕉在线观看| www国产在线视频色| 999久久久精品免费观看国产| 天堂中文最新版在线下载| xxxhd国产人妻xxx| 国产欧美日韩一区二区三| 色婷婷av一区二区三区视频| 国产高清视频在线播放一区| 男人的好看免费观看在线视频 | 久久99一区二区三区| 高潮久久久久久久久久久不卡| 欧美在线黄色| 国产亚洲欧美98| 激情视频va一区二区三区| 久久人妻熟女aⅴ| 精品国内亚洲2022精品成人| a级片在线免费高清观看视频| 久久久久久亚洲精品国产蜜桃av| 亚洲专区中文字幕在线| 久久精品亚洲av国产电影网| 亚洲欧美精品综合一区二区三区| 久久精品国产亚洲av高清一级| 日日夜夜操网爽| 久热这里只有精品99| 日韩一卡2卡3卡4卡2021年| 亚洲 欧美 日韩 在线 免费| a级毛片在线看网站| 韩国av一区二区三区四区| 欧美激情 高清一区二区三区| 成年女人毛片免费观看观看9| 每晚都被弄得嗷嗷叫到高潮| 亚洲五月婷婷丁香| 侵犯人妻中文字幕一二三四区| 免费观看人在逋| 免费观看精品视频网站| 国产精品免费一区二区三区在线| 午夜久久久在线观看| 亚洲少妇的诱惑av| 午夜老司机福利片| ponron亚洲| 一边摸一边抽搐一进一小说| 老熟妇乱子伦视频在线观看| 欧美日韩亚洲综合一区二区三区_| 怎么达到女性高潮| 亚洲av美国av| 国产精品一区二区在线不卡| 亚洲国产毛片av蜜桃av| 极品人妻少妇av视频| 一区二区三区激情视频| 1024视频免费在线观看| 9热在线视频观看99| 90打野战视频偷拍视频| 可以在线观看毛片的网站| 欧美黑人精品巨大| 国产亚洲欧美精品永久| av中文乱码字幕在线| 波多野结衣一区麻豆| 亚洲av五月六月丁香网| 中文字幕另类日韩欧美亚洲嫩草| 色播在线永久视频| 久久精品人人爽人人爽视色| 99精品在免费线老司机午夜| 国产成年人精品一区二区 | 在线天堂中文资源库| 夫妻午夜视频| 亚洲片人在线观看| 欧美日韩亚洲国产一区二区在线观看| 久久人人精品亚洲av| 欧美色视频一区免费| 亚洲五月色婷婷综合| 18禁美女被吸乳视频| 精品高清国产在线一区| 19禁男女啪啪无遮挡网站| 国产日韩一区二区三区精品不卡| 免费av毛片视频| 精品日产1卡2卡| 久久久久亚洲av毛片大全| 如日韩欧美国产精品一区二区三区| 欧美不卡视频在线免费观看 | 亚洲av五月六月丁香网| 高清黄色对白视频在线免费看| 咕卡用的链子| 俄罗斯特黄特色一大片| 91在线观看av| 女生性感内裤真人,穿戴方法视频| 91大片在线观看| 精品人妻在线不人妻| 性欧美人与动物交配| 天天躁夜夜躁狠狠躁躁| 操出白浆在线播放| 日韩免费av在线播放| 一区二区三区精品91| 一二三四在线观看免费中文在| 日本免费a在线| 在线国产一区二区在线| 琪琪午夜伦伦电影理论片6080| 成人18禁高潮啪啪吃奶动态图| 男女做爰动态图高潮gif福利片 | 91大片在线观看| 久久亚洲精品不卡| 国产麻豆69| 国产精品国产高清国产av| 国产一卡二卡三卡精品| 桃红色精品国产亚洲av| 在线观看免费高清a一片| 国产野战对白在线观看| 一个人观看的视频www高清免费观看 | 9热在线视频观看99| 深夜精品福利| av片东京热男人的天堂| 日韩免费高清中文字幕av| 久热这里只有精品99| 如日韩欧美国产精品一区二区三区| 久久久久久久久中文| 女人高潮潮喷娇喘18禁视频| 亚洲一码二码三码区别大吗| 国产av一区在线观看免费| 国产高清视频在线播放一区| 日日干狠狠操夜夜爽| 亚洲精品av麻豆狂野| 国产欧美日韩精品亚洲av| 久久久久久久久免费视频了| 人人妻人人爽人人添夜夜欢视频| 精品一区二区三区av网在线观看| 国产黄a三级三级三级人| 国产在线精品亚洲第一网站| 88av欧美| 啦啦啦韩国在线观看视频| 深夜精品福利| 成人一区二区视频在线观看| 丰满乱子伦码专区| 国产熟女xx| 中文字幕熟女人妻在线| 91在线观看av| 国产黄色小视频在线观看| 国产免费男女视频| 99国产精品一区二区蜜桃av| 看十八女毛片水多多多| 日韩欧美精品免费久久 | 99在线视频只有这里精品首页| 国产午夜福利久久久久久| 中亚洲国语对白在线视频| 网址你懂的国产日韩在线| 丝袜美腿在线中文| 久久国产精品影院| 亚洲最大成人手机在线| 日韩欧美 国产精品| 精品欧美国产一区二区三| 日本黄大片高清| 1024手机看黄色片| 成人国产一区最新在线观看| 老司机午夜十八禁免费视频| 日韩欧美 国产精品| 成人av在线播放网站| 欧美一区二区精品小视频在线| 精品人妻熟女av久视频| 91字幕亚洲| 国产中年淑女户外野战色| 99国产极品粉嫩在线观看| 国内精品久久久久久久电影| 色在线成人网| 精品99又大又爽又粗少妇毛片 | 免费人成视频x8x8入口观看| 欧美日韩综合久久久久久 | 国产淫片久久久久久久久 | 俺也久久电影网| 成年女人看的毛片在线观看| 亚洲人成电影免费在线| av中文乱码字幕在线| 欧美3d第一页| 国产不卡一卡二| 久久精品国产亚洲av天美| 欧美午夜高清在线| 国产色婷婷99| 亚洲无线在线观看| 亚洲国产日韩欧美精品在线观看| 男女之事视频高清在线观看| 国产精品98久久久久久宅男小说| 国产亚洲av嫩草精品影院| 999久久久精品免费观看国产| 国产淫片久久久久久久久 | 久久性视频一级片| 精品人妻偷拍中文字幕| aaaaa片日本免费| 最好的美女福利视频网| 色av中文字幕| 久久中文看片网| 亚洲在线观看片| 精品久久国产蜜桃| 搡老熟女国产l中国老女人| 制服丝袜大香蕉在线| 国内精品一区二区在线观看| www.www免费av| 永久网站在线| 精品无人区乱码1区二区| 九色成人免费人妻av| 成熟少妇高潮喷水视频| 色吧在线观看| 午夜精品一区二区三区免费看| 精品久久久久久久久av| 亚洲一区二区三区不卡视频| 丝袜美腿在线中文| 精品国产亚洲在线| 亚洲人与动物交配视频| 成人三级黄色视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产欧洲综合997久久,| 在线国产一区二区在线| 制服丝袜大香蕉在线| 少妇高潮的动态图| 日本一二三区视频观看| 久久热精品热| 欧美中文日本在线观看视频| 男人舔女人下体高潮全视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 99久久精品国产亚洲精品| 久久国产乱子伦精品免费另类| 精品午夜福利视频在线观看一区| 国内揄拍国产精品人妻在线| 久久久精品大字幕| 淫秽高清视频在线观看| 熟妇人妻久久中文字幕3abv| 毛片一级片免费看久久久久 | 日本成人三级电影网站| 国产亚洲av嫩草精品影院| 91字幕亚洲| 蜜桃亚洲精品一区二区三区| 国产高清激情床上av| 久久久精品大字幕| 淫秽高清视频在线观看| 99久久九九国产精品国产免费| 欧美xxxx黑人xx丫x性爽| 身体一侧抽搐| 免费在线观看亚洲国产| 国产高清激情床上av| 日本一本二区三区精品| 90打野战视频偷拍视频| 成人国产综合亚洲| 精品欧美国产一区二区三| 国产av在哪里看| 亚洲美女搞黄在线观看 | 久久精品久久久久久噜噜老黄 | 久久国产精品影院| 亚洲国产色片| 国产野战对白在线观看| 国产在视频线在精品| 床上黄色一级片| 亚洲无线观看免费| 免费搜索国产男女视频| 色哟哟哟哟哟哟| xxxwww97欧美| 一进一出好大好爽视频| 中文字幕高清在线视频| 女人十人毛片免费观看3o分钟| 女同久久另类99精品国产91| 简卡轻食公司| 亚洲成人中文字幕在线播放| 观看免费一级毛片| 黄色视频,在线免费观看| 亚洲色图av天堂| 真人一进一出gif抽搐免费| 国产三级中文精品| 精品久久久久久成人av| www.999成人在线观看| 国产精品亚洲美女久久久| 色精品久久人妻99蜜桃| 国产爱豆传媒在线观看| 如何舔出高潮| 一二三四社区在线视频社区8| 91在线观看av| av在线天堂中文字幕| 看黄色毛片网站| 色综合站精品国产| 成人美女网站在线观看视频| 国产蜜桃级精品一区二区三区| 欧美另类亚洲清纯唯美| 精品人妻视频免费看| 夜夜夜夜夜久久久久| av中文乱码字幕在线| 在线免费观看不下载黄p国产 | 久久精品91蜜桃| a在线观看视频网站| 在线观看免费视频日本深夜| 俺也久久电影网| 99视频精品全部免费 在线| a在线观看视频网站| 在线观看免费视频日本深夜| 欧美xxxx黑人xx丫x性爽| 日本成人三级电影网站| 午夜福利在线观看吧| 亚洲av五月六月丁香网| 日韩av在线大香蕉| 伦理电影大哥的女人| 欧美成人性av电影在线观看| 黄色一级大片看看| 免费搜索国产男女视频| h日本视频在线播放| a在线观看视频网站| 亚洲真实伦在线观看| 日本三级黄在线观看| 欧美性猛交黑人性爽| 成人国产综合亚洲| 国产精品美女特级片免费视频播放器| 搡老岳熟女国产| 亚洲三级黄色毛片| 搡老妇女老女人老熟妇| 国产午夜福利久久久久久| 亚洲精品456在线播放app | 欧美xxxx黑人xx丫x性爽| 91字幕亚洲| 欧美高清性xxxxhd video| АⅤ资源中文在线天堂| 搞女人的毛片| 1024手机看黄色片| 搡老岳熟女国产| 免费人成在线观看视频色| 欧美成人性av电影在线观看| 91久久精品电影网| 99久久精品国产亚洲精品| 亚洲国产精品sss在线观看| 亚洲久久久久久中文字幕| 我的老师免费观看完整版| 在线天堂最新版资源| 亚洲成av人片免费观看| 又粗又爽又猛毛片免费看| 亚洲第一欧美日韩一区二区三区| 久久精品国产自在天天线| 欧洲精品卡2卡3卡4卡5卡区| 婷婷色综合大香蕉| 成熟少妇高潮喷水视频| 国产不卡一卡二| 日韩欧美在线二视频| 久久中文看片网| 久久久久久久亚洲中文字幕 | 国产单亲对白刺激| 色噜噜av男人的天堂激情| 亚洲不卡免费看| 亚洲av.av天堂| 国产欧美日韩精品一区二区| 听说在线观看完整版免费高清| 日韩精品中文字幕看吧| 日韩av在线大香蕉| 黄色日韩在线| 一本久久中文字幕| 成人av一区二区三区在线看| 精品人妻熟女av久视频| 国产野战对白在线观看| 精品午夜福利视频在线观看一区| 九九久久精品国产亚洲av麻豆| 中文字幕人妻熟人妻熟丝袜美| 久久人人爽人人爽人人片va | 欧美一区二区精品小视频在线| 国产成人av教育| 无遮挡黄片免费观看| 国产v大片淫在线免费观看| 亚洲美女视频黄频| 中文字幕人妻熟人妻熟丝袜美| 少妇的逼好多水| 国产亚洲精品av在线| 亚洲黑人精品在线|