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

    平原河網(wǎng)十字型汊口出入流斷*面污染物濃度響應(yīng)關(guān)系研究

    2011-10-24 10:59:12韓龍喜
    關(guān)鍵詞:河網(wǎng)水槽流速

    韓龍喜,陸 東,計(jì) 紅

    (1.河海大學(xué)水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇南京 210098;2.河海大學(xué)淺水湖泊綜合治理與資源開發(fā)教育部重點(diǎn)實(shí)驗(yàn)室,江蘇南京 210098;3.河海大學(xué)環(huán)境科學(xué)與工程學(xué)院,江蘇南京 210098)

    平原河網(wǎng)十字型汊口出入流斷*面污染物濃度響應(yīng)關(guān)系研究

    韓龍喜1,2,3,陸 東3,計(jì) 紅3

    (1.河海大學(xué)水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇南京 210098;2.河海大學(xué)淺水湖泊綜合治理與資源開發(fā)教育部重點(diǎn)實(shí)驗(yàn)室,江蘇南京 210098;3.河海大學(xué)環(huán)境科學(xué)與工程學(xué)院,江蘇南京 210098)

    在平原河網(wǎng)地區(qū),十字型汊口是不同河道間實(shí)行水量與污染物交換的重要水域,而出、入流斷面污染物濃度間的響應(yīng)關(guān)系,是模擬平原河網(wǎng)水質(zhì)時(shí)空變化的一個(gè)十分重要的技術(shù)參數(shù)。本文通過十字型水槽物理模型,對(duì)“兩進(jìn)兩出”十字型汊口二維水動(dòng)力、水質(zhì)數(shù)學(xué)模型進(jìn)行驗(yàn)證。結(jié)果表明,在“兩進(jìn)兩出”十字型汊口不同出流斷面的污染物平均濃度存在較大差異。在此基礎(chǔ)上,提出了新的該類型汊口出入流斷面污染物平均濃度間的響應(yīng)關(guān)系,給出了響應(yīng)系數(shù)計(jì)算公式。將利用本文方法得出的響應(yīng)系數(shù)與常用的三級(jí)聯(lián)解法基于“汊口污染物均勻混合假設(shè)”得出的響應(yīng)系數(shù)進(jìn)行了對(duì)比,突顯了三級(jí)聯(lián)解法的誤差。本文的研究成果可為平原河網(wǎng)水質(zhì)模擬方法的改進(jìn)提供依據(jù)。

    河網(wǎng);十字型汊口;污染物濃度;響應(yīng)關(guān)系

    平原河網(wǎng)中河流交匯現(xiàn)象十分普遍,許多單一河道由汊口 (交匯點(diǎn))相互連接形成河網(wǎng)系統(tǒng)。目前,數(shù)學(xué)模型仍然是研究平原河網(wǎng)地區(qū)污染物濃度的時(shí)間、空間分布規(guī)律的一個(gè)不可或缺的重要工具。一般來說,汊口污染物遷移特征十分復(fù)雜,而汊口出流斷面與入流斷面污染物濃度間的響應(yīng)關(guān)系,是河網(wǎng)水質(zhì)數(shù)學(xué)模擬的一個(gè)關(guān)鍵技術(shù)問題。在現(xiàn)有的平原河網(wǎng)水質(zhì)模擬方法中,三級(jí)聯(lián)解法是最為常用的方法,而該方法在處理汊口處污染物混合擴(kuò)散這一關(guān)鍵問題時(shí),主要采用“均勻混合假設(shè)”——流入汊口的污染物在汊口充分混合達(dá)到均勻程度,故汊口所有流出斷面的污染物平均濃度均等于入流斷面污染物充分混合后的濃度[1-2]。顯然,此假設(shè)有一定的適用范圍——連接汊口河道為“一進(jìn)多出”或 “多進(jìn)一出”。而實(shí)際上,平原河網(wǎng)地區(qū)汊口河道以 “兩進(jìn)兩出”最為常見,由于達(dá)到均勻混合所需縱向距離較長(zhǎng)[3-4],或流入汊口的河道污染物濃度差異較大時(shí),污染物根本不能達(dá)到均勻混合,假設(shè)必然帶來系統(tǒng)誤差。因此,應(yīng)充分考慮汊口出流斷面濃度的不均勻性,并研制與之相匹配的平原河網(wǎng)水質(zhì)模型[5]。

    本文采用物理模型和數(shù)學(xué)模型相結(jié)合的方法研究“兩進(jìn)兩出”十字型汊口流體的水動(dòng)力特性和污染物平面二維分布特征。在此基礎(chǔ)上,提出了新的該類型汊口出、入流斷面的污染物濃度響應(yīng)關(guān)系,給出了響應(yīng)系數(shù)計(jì)算公式。構(gòu)造典型算例,將所得結(jié)果與常用的 “汊口污染物均勻混合假設(shè)”的相應(yīng)結(jié)果進(jìn)行對(duì)比,突顯了 “均勻混合假設(shè)”的系統(tǒng)誤差。本文研究成果可為平原河網(wǎng)水質(zhì)模擬方法的改進(jìn)提供技術(shù)依據(jù)。流流量為 0.61 L/s,水槽 A出流流量為 0.82 L/s,水槽B入流流量為 0.61 L/s,水槽 A出流水深為6.6 cm;工況二:水槽A入流流量為 0.49 L/s,水槽A出流流量為 0.85 L/s,水槽 B入流流量為0.73 L/s,水槽 A出流水深為 6.6 cm。水槽模型尺寸及測(cè)點(diǎn)布置見圖 1。

    圖 1 水槽模型尺寸和測(cè)點(diǎn)布置圖 (單位:cm)Fig.1 Sizes of flume model and distribution ofmeasurement points

    1 汊口水動(dòng)力學(xué)特性研究

    1.1 汊口水動(dòng)力特性試驗(yàn)設(shè)置

    1.1.1 物理模型設(shè)計(jì) 采用室內(nèi)水槽模擬河道的水流、污染物輸運(yùn)特征。物理模型試驗(yàn)成果按一定的比尺換算后,可用于水動(dòng)力數(shù)值模擬結(jié)果的率定和驗(yàn)證。根據(jù)試驗(yàn)場(chǎng)地及供水條件,采用變態(tài)比尺設(shè)計(jì)水槽模型,水平比尺為 100∶1;垂向比尺為30∶1[6-7]。具體模型尺寸及流速測(cè)點(diǎn)布置見圖 1。共布置了 10個(gè)監(jiān)測(cè)斷面 35個(gè)測(cè)點(diǎn),其中,在摻混劇烈的汊口處測(cè)點(diǎn)較密集,共布置了 4個(gè)監(jiān)測(cè)斷面,19個(gè)測(cè)點(diǎn)。

    1.1.2 物理模型試驗(yàn)工況 對(duì)汊口入流、出流水動(dòng)力條件設(shè)置了 2種試驗(yàn)工況,研究分析汊口水動(dòng)力特征。共有 2種試驗(yàn)工況,工況一:水槽 A入

    1.2 汊口水動(dòng)力數(shù)值模擬

    1.2.1 控制方程 平原河網(wǎng)地區(qū)河流一般為寬淺型,可采用水深平均的平面二維水動(dòng)力數(shù)學(xué)模型數(shù)值模擬垂線平均流速的平面分布。通過求解黎曼問題計(jì)算汊口流場(chǎng)[8]。

    1.2.2 邊界條件 根據(jù)物理模型幾何相似率,得到流速比尺為 5.48,流量比尺為 16 640,據(jù)此對(duì)物理模型兩種工況下的入流、出流邊界條件進(jìn)行換算,得到數(shù)學(xué)模型模擬的原型河道的水動(dòng)力邊界條件。

    1.2.3 網(wǎng)格剖分 采用均勻矩形網(wǎng)格,網(wǎng)格大小設(shè)定為 4 m,計(jì)算域?yàn)閮蓷l各為長(zhǎng) 920 m,寬 40 m的十字交叉形河道,其中汊口上游 80 m,下游 800 m。共 4 500個(gè)計(jì)算單元,4 591個(gè)計(jì)算節(jié)點(diǎn)。

    1.2.4 水動(dòng)力模擬結(jié)果分析 依據(jù)流速比尺,將水槽物理試驗(yàn)結(jié)果換算成相應(yīng)的原型河道水流流速,采用工況 1物理試驗(yàn)數(shù)據(jù)對(duì)汊口二維水動(dòng)力模型進(jìn)行率定,采用工況 2進(jìn)行驗(yàn)證。流場(chǎng)率定結(jié)果見圖 2(a),驗(yàn)證結(jié)果見圖 2(b)。由模型率定、驗(yàn)證結(jié)果對(duì)比可見,流速相對(duì)誤差小于 30%的占80%以上,計(jì)算流場(chǎng)與實(shí)測(cè)流場(chǎng)基本一致。

    圖 2 模型率定及驗(yàn)證時(shí)水槽試驗(yàn)及數(shù)值模擬流速對(duì)比Fig.2 Comparison of velocity between physical and mathematicalmodels in confirmation and validation

    1.3 試驗(yàn)結(jié)果及計(jì)算結(jié)果對(duì)比分析

    本文僅給出工況 2數(shù)值模擬及物理試驗(yàn)流場(chǎng)(見圖 3)。由圖 3可知,對(duì)于汊口入流河道,流線基本順直,流場(chǎng)受汊口影響相對(duì)較小;對(duì)于汊口,流速的流向、大小具有很大的空間分布不均勻性,表現(xiàn)為流向急劇變化,最大變幅近 90°,流速差異較大,具體表現(xiàn)為近岸流速較小,汊口中心流速較大;而對(duì)于汊口出流河道,由于受到汊口水流流態(tài)的影響,靠近汊口的出流斷面水流流速呈現(xiàn)出較大的不均勻性,具體表現(xiàn)為水流流向沿橫斷面方向存在較大的切變,流速大小沿橫斷面分布也很不均勻:臨近入流斷面一側(cè)流速較小,遠(yuǎn)離入流斷面一側(cè)由于兩入流斷面水流的擠壓及水流的慣性作用,流速相對(duì)較大。隨著水流向下游流動(dòng),水流在東面分布上的差異特征逐步消失。

    圖 3 工況 2數(shù)值模擬及物理試驗(yàn)流場(chǎng)Fig.3 Flow fields bymathematicalmodel and by physicalmodel

    2 汊口擴(kuò)散質(zhì)輸運(yùn)規(guī)律研究

    2.1 汊口擴(kuò)散質(zhì)輸運(yùn)規(guī)律試驗(yàn)研究

    進(jìn)行污染物在水體中的對(duì)流擴(kuò)散的模型試驗(yàn)時(shí),需考慮所研究的擴(kuò)散物質(zhì)在不同比尺水體中的對(duì)流擴(kuò)散相似規(guī)律,這樣才能較為準(zhǔn)確的將模型試驗(yàn)數(shù)據(jù)推廣到原型水體中。王美敬等[9-10]在水動(dòng)力學(xué)條件相似的基礎(chǔ)上,推導(dǎo)出了污染物擴(kuò)散-降解模型所遵循的相似準(zhǔn)則并進(jìn)行了驗(yàn)證,結(jié)果表明:在具有明顯二維特征的原型與模型試驗(yàn)中,模型所獲得成果可較好地應(yīng)用到實(shí)際河流中。本文借鑒前人的研究成果采用適當(dāng)?shù)南嗨茥l件進(jìn)行汊口污染物擴(kuò)散的物理模型試驗(yàn)。

    2.1.1 試驗(yàn)工況 水流工況分 2種:工況 1和工況2。

    2.1.2 示蹤劑及濃度設(shè)定 試驗(yàn)中利用碳素墨水作為示蹤劑,通過濁度儀測(cè)量透光度作為流體中污染物的當(dāng)量濃度值。在A水槽中投放示蹤劑,B水槽進(jìn)水為清水。設(shè)定清水中示蹤劑的濃度為 0,投放的示蹤劑碳素墨水濃度設(shè)為 6。

    2.2 汊口擴(kuò)散質(zhì)輸運(yùn)數(shù)值模擬

    應(yīng)用平面二維對(duì)流擴(kuò)散方程描述污染物的輸運(yùn)擴(kuò)散,通過求解黎曼問題得到污染物垂線平均濃度在平面上的分布特征。

    2.2.1 控制方程

    式中,h為水深;u、ν分別為x、y方向垂線平均水平流速分量;Ci為污染物的垂線平均濃度;Eix、Eiy分別為x、y方向各污染物的擴(kuò)散系數(shù);KCi為污染物綜合降解系數(shù)。

    2.2.2 定解條件 給定初始條件為c(x,y,0)=c0(x,y),對(duì)岸邊界采用邊界不可入條件,入流邊界給出各節(jié)點(diǎn)的污染物濃度,出流邊界采用第二類邊界條件,即=0。網(wǎng)格的布置同流場(chǎng)計(jì)算。計(jì)算濃度時(shí)兩個(gè)上游入口邊界給定濃度同物模試驗(yàn),A河道上游邊界濃度為 6,B河道上游邊界濃度為0??v向擴(kuò)散系數(shù)Ex、Ey由經(jīng)驗(yàn)公式確定[11]。污染物降解系數(shù)KCi取為 0(不考慮降解)。

    2.3 試驗(yàn)、計(jì)算結(jié)果對(duì)比及汊口濃度分布分析

    本文僅給出工況 1物理試驗(yàn)、數(shù)值模擬的示蹤劑濃度場(chǎng) (見圖 4)。根據(jù)所有觀測(cè)點(diǎn)上物理試驗(yàn)數(shù)據(jù)、數(shù)值模擬結(jié)果的統(tǒng)計(jì)分析可知,觀測(cè)濃度與數(shù)值模擬結(jié)果差值小于 1mg/L的點(diǎn)據(jù)占 83%,計(jì)算濃度場(chǎng)與實(shí)測(cè)濃度場(chǎng)基本一致。

    根據(jù)工況 1、2示蹤劑濃度平面分布的物理試驗(yàn)、數(shù)值模擬結(jié)果可知,在汊口區(qū)域,示蹤劑濃度沿水流運(yùn)動(dòng)的法向,存在較大的濃度梯度,靠近低濃度入流河道一側(cè)汊口水域濃度相對(duì)較低,反之則較高。且由于流程較短,示蹤劑在汊口混合程度較低,因此,相應(yīng)地導(dǎo)致了不同出流河道斷面平均濃度的巨大差異,即緊鄰低入流濃度河道的出流河道,由于其來水組成大部分為低濃度來水,因此該出流斷面示蹤劑濃度相對(duì)較低,而緊鄰高入流濃度河道的出流河道,由于其來水組成大部分為高濃度來水,因此出流濃度相對(duì)較高。

    3 出入流斷面濃度響應(yīng)關(guān)系分析

    汊口出流斷面的平均濃度與汊口流入斷面的污染物平均濃度間的響應(yīng)關(guān)系,是平原河網(wǎng)水質(zhì)模擬的關(guān)鍵技術(shù)。本文提出基于汊口污染物非均勻混合的出入流斷面濃度響應(yīng)系數(shù)的求解方法,并構(gòu)造典型算例采用本文的計(jì)算方法計(jì)算響應(yīng)系數(shù),然后與常用的基于均勻混合假設(shè)的三級(jí)聯(lián)合解法的響應(yīng)關(guān)系進(jìn)行比較分析。

    圖 4 工況 1物理模型及數(shù)值模擬示蹤劑濃度場(chǎng)Fig.4 Concentration field bymathematicalmodel and by physicalmodel

    3.1 均勻混合假設(shè)概述

    假設(shè)有如圖 5所示的 “兩進(jìn)兩出”十字型汊口 (圖中箭頭指向表示水流方向),Q1、Q2分別為 A河道、B河道的入流流量,C1、C2為相應(yīng)的入流污染物質(zhì)量濃度;Q3、Q4為 A河道、B河道的出流流量,C3、C4為相應(yīng)的出流污染物質(zhì)量濃度。

    均勻混合假設(shè)認(rèn)為,所有汊口出流斷面的污染物平均濃度,均等于汊口節(jié)點(diǎn)濃度——入流斷面污染物在交叉點(diǎn)均勻混合后的濃度。數(shù)學(xué)表示為[1]:

    可將出流斷面平均濃度簡(jiǎn)單表示為入流斷面平均濃度的線性函數(shù),有:

    圖 5 汊口污染物混合輸運(yùn)示意圖Fig.5 Schematic illustration of pollutantmixing and transferring at joint nodes

    3.2 非均勻混合假設(shè)響應(yīng)關(guān)系數(shù)學(xué)表達(dá)

    十字形汊口污染物非均勻混合條件下,出流斷面濃度可表示為入流斷面平均濃度的線性組合,數(shù)學(xué)表示為:

    式中,α′1,α′2分別為河道 A出流污染物濃度對(duì)河道 A、B入流污染物濃度的響應(yīng)系數(shù),β′1,β′2分別為河道 B出口污染物濃度對(duì)河道 A、B入流污染物濃度的響應(yīng)系數(shù)。

    3.3 響應(yīng)系數(shù)計(jì)算方法

    天然河道進(jìn)入汊口的污染物濃度通常不等,不妨假設(shè)C2<C1,令C2=Cb,則C1可表示為C1=Cb+△C,代入式 (5)、(6)得:

    A河道入流斷面污染物濃度與B河道相比存在變幅△C,不妨設(shè)由此在A河道出流斷面產(chǎn)生的相應(yīng)濃度增量為△C3,在 B河道出流斷面產(chǎn)生的相應(yīng)濃度增量響應(yīng)為△C4。則有:C3=Cb+△C3,C4=Cb+△C4。將此關(guān)系帶入到方程 (7),得到:

    同理可得,

    由式 (8)、(9)得:

    式中,B3為A河道出流斷面的寬度,B4表示B河道出流斷面的寬度。u(b)、h(b)、c(b)分別表示出流河道斷面的流速、水深及垂線平均濃度沿橫向的分布,可由上述汊口平面二維數(shù)值模擬計(jì)算得到。

    3.4 典型算例響應(yīng)系數(shù)對(duì)比

    計(jì)算了 3種典型算例的出流斷面污染物濃度響應(yīng)系數(shù),并與基于均勻混合假設(shè)的響應(yīng)系數(shù)進(jìn)行了對(duì)比,結(jié)果見表 1。

    表 1 響應(yīng)系數(shù)對(duì)照表Table1 Comparison of response coefficients be tween different hypotheses

    由表 1可知,出流河道的濃度響應(yīng)系數(shù)與入、出流河道的流量分配有關(guān)?;诰鶆蚧旌霞僭O(shè)的響應(yīng)系數(shù)與基于非均勻混合理論的響應(yīng)系數(shù)相比,在數(shù)值上存在較大的偏差。

    4 小 結(jié)

    1)水槽試驗(yàn)及數(shù)學(xué)模型模擬結(jié)果顯示,十字型汊口水域,由于受到進(jìn)入汊口水體污染物濃度差異的影響,擴(kuò)散質(zhì)摻混很不均勻,污染物濃度在汊口存在很大的濃度梯度,出流斷面污染物濃度沿橫向的分布存在顯著差異。

    2)提出了汊口出流斷面平均濃度對(duì)入流斷面污染物平均濃度的響應(yīng)系數(shù)的計(jì)算方法,研究結(jié)果表明,不同出流斷面的污染物濃度響應(yīng)系數(shù)不等。

    3)研究結(jié)果對(duì)比表明,基于汊口污染物均勻混合假設(shè)的出流斷面濃度響應(yīng)系數(shù)存在較大偏差。在進(jìn)行平原河網(wǎng)水質(zhì)模擬時(shí),應(yīng)考慮入流斷面水質(zhì)不均勻性對(duì)流出汊口水質(zhì)的影響,研制更為合理的平原河網(wǎng)水質(zhì)模型。

    [1]韓龍喜,朱黨生 .河網(wǎng)地區(qū)水環(huán)境規(guī)劃中的污染源控制方法[J].水利學(xué)報(bào),2001(10):28-31.

    [2]吳作平.河網(wǎng)水流數(shù)值模擬方法研究 [J].水科學(xué)進(jìn)展,2003,14(3):350-353.

    [3]張書農(nóng).環(huán)境水力學(xué) [M].南京:河海大學(xué)出版社,1988:134-139.

    [4]CHO I G W,MOL INAS A.Simultaneous solution algorithm for channel network modeling[J].Water Resource,1993,29(2):321-328.

    [5]韓龍喜,陸冬.平原河網(wǎng)水流水質(zhì)數(shù)值模擬研究進(jìn)展[J].河海大學(xué)學(xué)報(bào):自然科學(xué)版,2004,32(2):127-130.

    [6]汪德勝 .水力變態(tài)模型相似律問題的探討[J].泥沙研究,1987(4):72—79.

    [7]鄭國(guó)棟,鄭邦民.變態(tài)模型水流相似的精度及誤差分析[J].水動(dòng)力研究與進(jìn)展,2005,20(1):33-37.

    [8]胡四一,譚維炎.一維不恒定明流計(jì)算的三種高性能差分格式[J].水利學(xué)報(bào),1991,2(1):11-21.

    [9]王美敬,羅麟.水中污染物擴(kuò)散模型實(shí)驗(yàn)中的相似理論[J].四川大學(xué)學(xué)報(bào):工程科學(xué)版,2004,36(2):25-28.

    [10]王美敬,羅麟,黃小雪 .污染物橫向擴(kuò)散在變態(tài)模型與河流中的相似性[J].四川大學(xué)學(xué)報(bào):工程科學(xué)版,2006,38(4):15-18.

    [11]李錦秀,黃真理,呂平毓 .三峽庫(kù)區(qū)江段縱向離散系數(shù)研究[J].水利學(xué)報(bào),2000(8):89-94.

    Response Relation of Pollutant Concentration Between Outlet and Inlet Sections at Joint Nodes of Pla in River Network

    HAN Longxi1,2,3,LU Dong3,JI Hong3
    (1.State KeyLaboratory of Hydrology-Water Resource and Hydraulic Engineering,Nanjing 210098,China;2.KeyLaboratory of Integrated Regulation and Resource Development on Shallow Lakes,MOE,Nanjing 210098,China;3.College of Environmental Science and Engineering,HohaiUniversity,Nanjing 210098,China)

    In plain river ne twork,joint nodes can not be ignored for exchange of water flow and pollutants from differentwatercourses.Response relation between outlet sections and inlet sections is an important technical parameter forwater quality simulation.Aiming at the response relation research of crossed joint nodeswith two inlet sections and two outlet sections,thispaper builds2-D hydraulic andwater qualitymathematicalmodels and validates them by physical models.The results show that pollutant concentration differs greatly for different outlet sections.Based on the work above,a new expression of response relation and calculating equations of response coefficients are put for ward.By comparing the response coefficientswith those inferred from even mixing hypothesis mainly used in three-step method,obvious errors caused by even mixing hypothesis are disclosed.This study helps to improve the simulation method of plain river network water quality.

    plain river network;crossed joint nodes;pollutant concentration;response relation

    TU991

    B

    0529-6579(2011)01-0123-06

    2009-11-16

    水利部公益性行業(yè)科研專項(xiàng)經(jīng)費(fèi)資助項(xiàng)目 (201001028)

    韓龍喜 (1964年生),男,教授;E-mail:hanlongxi@sina.com

    猜你喜歡
    河網(wǎng)水槽流速
    基于小世界網(wǎng)絡(luò)的海河流域河網(wǎng)結(jié)構(gòu)及功能響應(yīng)
    “流體壓強(qiáng)與流速的關(guān)系”知識(shí)鞏固
    『流體壓強(qiáng)與流速的關(guān)系』知識(shí)鞏固
    可升降折疊的飲水機(jī)水槽
    山雨欲來風(fēng)滿樓之流體壓強(qiáng)與流速
    基于PSR模型的上海地區(qū)河網(wǎng)脆弱性探討
    為什么水槽管要做成彎曲狀
    愛虛張聲勢(shì)的水
    不同引水水源對(duì)平原河網(wǎng)影響分析
    要挑好水槽,就看這里了!
    幸福(2016年6期)2016-12-01 03:08:13
    成人特级av手机在线观看| 国产极品天堂在线| av国产免费在线观看| 五月伊人婷婷丁香| 国产精品三级大全| 18禁在线播放成人免费| 视频中文字幕在线观看| 精品一区二区免费观看| 我的女老师完整版在线观看| 国产精品国产三级国产av玫瑰| 综合色丁香网| 亚洲高清免费不卡视频| 国产精品美女特级片免费视频播放器| 日韩高清综合在线| 99热这里只有是精品在线观看| 老师上课跳d突然被开到最大视频| 国产精品野战在线观看| 看十八女毛片水多多多| kizo精华| 久久亚洲国产成人精品v| 真实男女啪啪啪动态图| 久久99热6这里只有精品| 午夜免费激情av| 超碰av人人做人人爽久久| 精品久久久久久久末码| 久久精品久久精品一区二区三区| 精品人妻视频免费看| 只有这里有精品99| 亚洲第一区二区三区不卡| 一级黄色大片毛片| 亚洲人成网站在线观看播放| 国产精品电影一区二区三区| 天堂av国产一区二区熟女人妻| 高清视频免费观看一区二区 | 男女那种视频在线观看| 国产在线男女| 欧美不卡视频在线免费观看| 边亲边吃奶的免费视频| 日本午夜av视频| 国产欧美另类精品又又久久亚洲欧美| 一级二级三级毛片免费看| 天天躁日日操中文字幕| 国产片特级美女逼逼视频| 日本黄色片子视频| 热99re8久久精品国产| 亚洲精品456在线播放app| 乱码一卡2卡4卡精品| 国产探花极品一区二区| 久久久久精品久久久久真实原创| 午夜老司机福利剧场| 亚洲图色成人| 亚洲av二区三区四区| 久久久午夜欧美精品| 哪个播放器可以免费观看大片| 看非洲黑人一级黄片| av国产免费在线观看| 亚洲精品久久久久久婷婷小说 | 一区二区三区免费毛片| 十八禁国产超污无遮挡网站| av在线蜜桃| 国产男人的电影天堂91| 久久精品夜色国产| 久久亚洲国产成人精品v| av在线观看视频网站免费| 一级二级三级毛片免费看| 啦啦啦啦在线视频资源| 99久久九九国产精品国产免费| 欧美日韩精品成人综合77777| 日本爱情动作片www.在线观看| 又爽又黄无遮挡网站| 国产伦在线观看视频一区| 亚洲乱码一区二区免费版| 欧美性猛交╳xxx乱大交人| 国产视频内射| 岛国毛片在线播放| 国产精品一及| 人妻少妇偷人精品九色| 特级一级黄色大片| 热99在线观看视频| 99久国产av精品国产电影| 天堂网av新在线| 日本黄大片高清| 天堂√8在线中文| 精品人妻熟女av久视频| 国产免费男女视频| 久久99热这里只有精品18| 国模一区二区三区四区视频| 大又大粗又爽又黄少妇毛片口| 日韩欧美 国产精品| 精品久久国产蜜桃| 女人被狂操c到高潮| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲美女视频黄频| 日韩亚洲欧美综合| 亚洲精华国产精华液的使用体验| av免费观看日本| 久久韩国三级中文字幕| 中文字幕久久专区| av播播在线观看一区| 免费观看人在逋| 91精品国产九色| www.av在线官网国产| 国内揄拍国产精品人妻在线| 一级毛片我不卡| av专区在线播放| 永久免费av网站大全| 九色成人免费人妻av| 欧美xxxx性猛交bbbb| 免费一级毛片在线播放高清视频| 欧美日韩精品成人综合77777| 国产成人a∨麻豆精品| 高清在线视频一区二区三区 | 久久6这里有精品| 精品久久久久久成人av| 国内精品美女久久久久久| 亚洲欧洲日产国产| 99久国产av精品| 网址你懂的国产日韩在线| 国内揄拍国产精品人妻在线| 国产精华一区二区三区| 久久精品熟女亚洲av麻豆精品 | 男人舔女人下体高潮全视频| 久久久色成人| 大又大粗又爽又黄少妇毛片口| 91午夜精品亚洲一区二区三区| 激情 狠狠 欧美| 干丝袜人妻中文字幕| 亚洲在久久综合| 亚洲美女视频黄频| 国产精品久久久久久精品电影小说 | 日本色播在线视频| 少妇猛男粗大的猛烈进出视频 | 内地一区二区视频在线| 国产免费一级a男人的天堂| 久久久久久大精品| 国产一级毛片在线| 亚洲av日韩在线播放| 狂野欧美激情性xxxx在线观看| 超碰97精品在线观看| 亚洲av二区三区四区| 国产成年人精品一区二区| 乱系列少妇在线播放| 91狼人影院| 午夜a级毛片| 国产精华一区二区三区| 熟妇人妻久久中文字幕3abv| 爱豆传媒免费全集在线观看| 国产精品国产三级国产专区5o | 99久国产av精品国产电影| 中文字幕av在线有码专区| 国产一区二区三区av在线| 嫩草影院新地址| 国产亚洲一区二区精品| 在线播放国产精品三级| 91精品一卡2卡3卡4卡| 夫妻性生交免费视频一级片| 日韩欧美精品v在线| 一级毛片电影观看 | 精品久久久久久久人妻蜜臀av| 国产成人精品久久久久久| 久久久久久久久大av| 国产高清有码在线观看视频| 麻豆乱淫一区二区| av专区在线播放| 精品一区二区三区人妻视频| 亚洲精品乱码久久久久久按摩| 久久久久精品久久久久真实原创| 亚洲内射少妇av| 18禁在线播放成人免费| 国产麻豆成人av免费视频| 欧美激情国产日韩精品一区| 国产精品人妻久久久影院| 午夜激情福利司机影院| 小说图片视频综合网站| 18禁在线播放成人免费| 最近最新中文字幕大全电影3| 国产爱豆传媒在线观看| 午夜福利在线观看吧| 成人一区二区视频在线观看| 草草在线视频免费看| 国产老妇伦熟女老妇高清| 久久久久免费精品人妻一区二区| 久久久久久久久大av| 少妇人妻精品综合一区二区| 你懂的网址亚洲精品在线观看 | 亚洲色图av天堂| 亚洲国产欧美在线一区| 青春草国产在线视频| 网址你懂的国产日韩在线| 九色成人免费人妻av| 91精品国产九色| 永久网站在线| 欧美成人a在线观看| 免费人成在线观看视频色| 婷婷色综合大香蕉| 99热全是精品| 国产精品一区二区性色av| 午夜亚洲福利在线播放| 国产精品电影一区二区三区| 美女国产视频在线观看| 免费看av在线观看网站| 全区人妻精品视频| 卡戴珊不雅视频在线播放| 大又大粗又爽又黄少妇毛片口| 亚洲国产精品sss在线观看| 日韩欧美精品免费久久| 亚洲国产精品合色在线| 国产国拍精品亚洲av在线观看| 干丝袜人妻中文字幕| 精品熟女少妇av免费看| 成人午夜高清在线视频| 人人妻人人看人人澡| 久久精品夜夜夜夜夜久久蜜豆| 国产人妻一区二区三区在| 少妇熟女aⅴ在线视频| 一边亲一边摸免费视频| 在线播放国产精品三级| 国产 一区精品| 午夜福利视频1000在线观看| 国产成人freesex在线| 久久精品影院6| 99热这里只有精品一区| 亚洲精品乱码久久久久久按摩| 亚洲五月天丁香| 国产午夜精品一二区理论片| 亚州av有码| 免费看日本二区| 永久免费av网站大全| 99国产精品一区二区蜜桃av| 联通29元200g的流量卡| 国产精品一及| 内地一区二区视频在线| 美女大奶头视频| 亚洲欧美日韩高清专用| av天堂中文字幕网| 最近中文字幕高清免费大全6| 免费人成在线观看视频色| 亚洲欧洲国产日韩| 美女高潮的动态| 成年女人看的毛片在线观看| 一个人观看的视频www高清免费观看| 国产在线男女| 国产国拍精品亚洲av在线观看| 免费黄网站久久成人精品| 久久久亚洲精品成人影院| 亚洲精品国产av成人精品| 国产伦一二天堂av在线观看| 国产亚洲最大av| 久久久久久久久久久丰满| 国产精品福利在线免费观看| 午夜激情欧美在线| 我的老师免费观看完整版| 一级二级三级毛片免费看| 成人美女网站在线观看视频| 国产精品av视频在线免费观看| 热99re8久久精品国产| 精品国产露脸久久av麻豆 | 亚洲欧美清纯卡通| 网址你懂的国产日韩在线| 听说在线观看完整版免费高清| 久久久精品大字幕| 日韩大片免费观看网站 | 男女啪啪激烈高潮av片| 日本欧美国产在线视频| 中文亚洲av片在线观看爽| 亚洲av一区综合| 国产免费又黄又爽又色| 免费不卡的大黄色大毛片视频在线观看 | 九九爱精品视频在线观看| 日本五十路高清| 久久久色成人| 永久网站在线| 免费观看的影片在线观看| 亚洲最大成人中文| 天美传媒精品一区二区| 麻豆一二三区av精品| 亚洲一级一片aⅴ在线观看| 午夜精品国产一区二区电影 | 两性午夜刺激爽爽歪歪视频在线观看| 亚洲av电影在线观看一区二区三区 | 午夜老司机福利剧场| 亚洲在久久综合| 国产精品日韩av在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 国产精品国产三级专区第一集| 国产精品熟女久久久久浪| 不卡视频在线观看欧美| 免费看美女性在线毛片视频| 国产激情偷乱视频一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美zozozo另类| 看黄色毛片网站| 亚洲精品色激情综合| 别揉我奶头 嗯啊视频| 2021天堂中文幕一二区在线观| 欧美色视频一区免费| 黄色日韩在线| 亚洲精品影视一区二区三区av| 国产午夜精品论理片| 97在线视频观看| 国产日韩欧美在线精品| 国产私拍福利视频在线观看| 日本三级黄在线观看| 日本色播在线视频| 国产精品久久久久久久久免| av播播在线观看一区| 最近2019中文字幕mv第一页| 91久久精品电影网| 好男人在线观看高清免费视频| 国语自产精品视频在线第100页| 免费大片18禁| 亚洲无线观看免费| 亚洲av电影在线观看一区二区三区 | 国产精品一区www在线观看| 人人妻人人澡人人爽人人夜夜 | 国内精品宾馆在线| 欧美三级亚洲精品| 午夜福利视频1000在线观看| 久久精品91蜜桃| 91精品伊人久久大香线蕉| 久久久久久久久久成人| 久久综合国产亚洲精品| 国产伦理片在线播放av一区| 91久久精品国产一区二区三区| 免费观看精品视频网站| 变态另类丝袜制服| 观看免费一级毛片| 久久久久久久久久黄片| 又黄又爽又刺激的免费视频.| 亚洲国产精品专区欧美| 亚洲精品一区蜜桃| 色综合亚洲欧美另类图片| 午夜福利成人在线免费观看| 成人三级黄色视频| 亚洲精品456在线播放app| 亚洲真实伦在线观看| 午夜激情福利司机影院| 亚洲真实伦在线观看| 精品无人区乱码1区二区| 欧美一区二区精品小视频在线| 亚洲欧美日韩高清专用| 亚洲国产精品合色在线| 亚洲欧洲日产国产| 最近最新中文字幕大全电影3| 精品久久久久久电影网 | 精品久久国产蜜桃| 中文字幕av在线有码专区| 亚洲国产欧美人成| 赤兔流量卡办理| 国产亚洲午夜精品一区二区久久 | 高清在线视频一区二区三区 | 国产一级毛片在线| 3wmmmm亚洲av在线观看| 欧美极品一区二区三区四区| 色噜噜av男人的天堂激情| 18禁动态无遮挡网站| 日本色播在线视频| 特大巨黑吊av在线直播| 国产不卡一卡二| 欧美成人免费av一区二区三区| 91av网一区二区| 免费看av在线观看网站| 黄片wwwwww| 亚洲精品456在线播放app| 深夜a级毛片| 国产一级毛片在线| 久久久久精品久久久久真实原创| 亚洲精品一区蜜桃| 久久久久久久久久黄片| 亚洲丝袜综合中文字幕| 国产高清国产精品国产三级 | 丰满少妇做爰视频| 老司机影院毛片| 欧美xxxx黑人xx丫x性爽| 网址你懂的国产日韩在线| 嫩草影院入口| 欧美区成人在线视频| 亚洲精品国产av成人精品| av.在线天堂| 男人的好看免费观看在线视频| 国产老妇女一区| 丰满少妇做爰视频| 日本黄大片高清| av女优亚洲男人天堂| 永久免费av网站大全| 成人性生交大片免费视频hd| 久久精品人妻少妇| 成人鲁丝片一二三区免费| 亚洲乱码一区二区免费版| 国产美女午夜福利| 夜夜看夜夜爽夜夜摸| 国产成人a∨麻豆精品| 亚洲成色77777| 中文字幕熟女人妻在线| 国产日韩欧美在线精品| 国产一区二区亚洲精品在线观看| 18+在线观看网站| 18禁在线播放成人免费| www.av在线官网国产| 国产成人a∨麻豆精品| 三级男女做爰猛烈吃奶摸视频| 午夜福利在线观看吧| 国产av码专区亚洲av| 亚洲精品国产av成人精品| 永久网站在线| 在线免费观看不下载黄p国产| 蜜臀久久99精品久久宅男| 韩国高清视频一区二区三区| 男人舔女人下体高潮全视频| 国产亚洲av片在线观看秒播厂 | 亚洲av一区综合| 久久鲁丝午夜福利片| 国内精品宾馆在线| 日韩av不卡免费在线播放| 久久久久性生活片| 日本色播在线视频| av黄色大香蕉| 精品久久久久久久久av| 日韩在线高清观看一区二区三区| 岛国在线免费视频观看| 插逼视频在线观看| 乱系列少妇在线播放| 麻豆久久精品国产亚洲av| 狂野欧美白嫩少妇大欣赏| 神马国产精品三级电影在线观看| 成人二区视频| 婷婷色综合大香蕉| 老师上课跳d突然被开到最大视频| 亚洲国产欧美在线一区| 国产精品人妻久久久久久| 日韩成人av中文字幕在线观看| 国产精品一区二区三区四区免费观看| 一本久久精品| 你懂的网址亚洲精品在线观看 | 97人妻精品一区二区三区麻豆| 精品久久久久久成人av| 国产成年人精品一区二区| 免费一级毛片在线播放高清视频| 综合色丁香网| 色网站视频免费| 毛片女人毛片| 看十八女毛片水多多多| 观看美女的网站| 色5月婷婷丁香| 国产淫语在线视频| 精品人妻视频免费看| 亚洲国产色片| 国产黄片美女视频| 久久韩国三级中文字幕| 在线播放无遮挡| 麻豆国产97在线/欧美| 成人美女网站在线观看视频| 黑人高潮一二区| 精品久久久久久久人妻蜜臀av| 少妇裸体淫交视频免费看高清| 久久久久久久久久久免费av| 1024手机看黄色片| 亚洲成人精品中文字幕电影| 国产乱人偷精品视频| 欧美精品国产亚洲| 亚洲国产日韩欧美精品在线观看| 丝袜美腿在线中文| av.在线天堂| 久久精品夜色国产| 久久精品久久久久久噜噜老黄 | 亚洲欧美成人综合另类久久久 | 中国国产av一级| 少妇高潮的动态图| 精品欧美国产一区二区三| 亚洲激情五月婷婷啪啪| 欧美极品一区二区三区四区| av又黄又爽大尺度在线免费看 | 在线免费观看的www视频| 成人午夜精彩视频在线观看| 日韩成人av中文字幕在线观看| 欧美最新免费一区二区三区| 在现免费观看毛片| 免费看a级黄色片| 欧美成人午夜免费资源| 国产成人福利小说| 看免费成人av毛片| 国产亚洲av嫩草精品影院| 久久欧美精品欧美久久欧美| 亚洲精品乱码久久久v下载方式| 亚洲精品乱久久久久久| 亚洲精华国产精华液的使用体验| av卡一久久| 成人毛片60女人毛片免费| 欧美丝袜亚洲另类| 日韩精品有码人妻一区| 免费人成在线观看视频色| 综合色丁香网| 久久精品91蜜桃| 日韩一区二区三区影片| 色综合亚洲欧美另类图片| 99久久精品一区二区三区| 国产高清有码在线观看视频| 大话2 男鬼变身卡| 婷婷色av中文字幕| 黑人高潮一二区| 日日摸夜夜添夜夜添av毛片| 成人无遮挡网站| 日本黄大片高清| 免费看av在线观看网站| 日韩欧美三级三区| 麻豆成人av视频| 国语对白做爰xxxⅹ性视频网站| 亚洲激情五月婷婷啪啪| 亚洲久久久久久中文字幕| 少妇的逼好多水| 亚洲最大成人中文| 精品一区二区三区视频在线| 久久国产乱子免费精品| 久久久精品大字幕| 国产黄片视频在线免费观看| 寂寞人妻少妇视频99o| 国产单亲对白刺激| 国产精品一区二区三区四区久久| 国产成人一区二区在线| 春色校园在线视频观看| av在线观看视频网站免费| 国产乱人偷精品视频| 综合色av麻豆| av福利片在线观看| 欧美xxxx性猛交bbbb| 国产精品av视频在线免费观看| 免费一级毛片在线播放高清视频| 一个人观看的视频www高清免费观看| 91久久精品电影网| 成年av动漫网址| 国产一区二区亚洲精品在线观看| h日本视频在线播放| 能在线免费观看的黄片| 亚洲欧美成人综合另类久久久 | 一区二区三区四区激情视频| 亚洲精品aⅴ在线观看| 国产精品伦人一区二区| 天堂√8在线中文| 亚洲国产精品合色在线| 精品人妻熟女av久视频| 九九久久精品国产亚洲av麻豆| 欧美变态另类bdsm刘玥| 男女视频在线观看网站免费| 22中文网久久字幕| 直男gayav资源| 丰满乱子伦码专区| 精品不卡国产一区二区三区| 99热这里只有是精品在线观看| 免费电影在线观看免费观看| 国产免费男女视频| 久久久国产成人免费| 久久久久精品久久久久真实原创| 精品久久久久久久人妻蜜臀av| 嫩草影院入口| av在线观看视频网站免费| 丝袜喷水一区| 最近中文字幕高清免费大全6| 寂寞人妻少妇视频99o| 女的被弄到高潮叫床怎么办| 又粗又硬又长又爽又黄的视频| 国内精品一区二区在线观看| 熟女人妻精品中文字幕| 国产伦精品一区二区三区视频9| 亚洲av中文字字幕乱码综合| 一边摸一边抽搐一进一小说| 黄片无遮挡物在线观看| 日本一本二区三区精品| 天美传媒精品一区二区| av在线蜜桃| 成人av在线播放网站| 日本免费一区二区三区高清不卡| 国产亚洲一区二区精品| 免费观看精品视频网站| 日韩大片免费观看网站 | 精品不卡国产一区二区三区| 熟女人妻精品中文字幕| 看非洲黑人一级黄片| 级片在线观看| 丰满乱子伦码专区| 岛国毛片在线播放| 男人和女人高潮做爰伦理| 中文欧美无线码| 国产乱人偷精品视频| 淫秽高清视频在线观看| 亚洲成色77777| 狠狠狠狠99中文字幕| .国产精品久久| 春色校园在线视频观看| 久久久久久九九精品二区国产| 亚洲一区高清亚洲精品| 最新中文字幕久久久久| 欧美一级a爱片免费观看看| 乱人视频在线观看| 三级男女做爰猛烈吃奶摸视频| 国产麻豆成人av免费视频| 亚洲欧洲国产日韩| 99久久精品一区二区三区| 亚洲av熟女| 看免费成人av毛片| 久久久国产成人免费| 亚洲国产精品sss在线观看| 中文在线观看免费www的网站| 色哟哟·www| 日韩,欧美,国产一区二区三区 | 中文资源天堂在线| 青春草亚洲视频在线观看| 国产探花极品一区二区| 国产亚洲av嫩草精品影院| 中文天堂在线官网| 久99久视频精品免费| 国产乱人视频| 99久久人妻综合| 亚洲欧洲日产国产| 亚洲成色77777| 久久人人爽人人片av|