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

    伊丹河輸水河道封凍期冰情演變數(shù)值模擬

    2017-11-10 07:55:12茅澤育齊文彪劉樹峰
    水利水電科技進(jìn)展 2017年6期
    關(guān)鍵詞:冰情冰蓋河段

    樊 霖,茅澤育,齊文彪,劉樹峰

    (1.清華大學(xué)水利水電工程系,北京 100084; 2.吉林省水利水電勘測(cè)設(shè)計(jì)研究院,吉林 長(zhǎng)春 130021)

    伊丹河輸水河道封凍期冰情演變數(shù)值模擬

    樊 霖1,茅澤育1,齊文彪2,劉樹峰2

    (1.清華大學(xué)水利水電工程系,北京 100084; 2.吉林省水利水電勘測(cè)設(shè)計(jì)研究院,吉林 長(zhǎng)春 130021)

    根據(jù)水動(dòng)力學(xué)、熱力學(xué)、河冰水力學(xué)及固體力學(xué)等基本理論,針對(duì)吉林省伊丹河輸水河段具體特征,建立了伊丹河的河冰數(shù)值模型,并應(yīng)用1958—1988共30年的水文、氣象資料對(duì)該輸水河段進(jìn)行了數(shù)值模擬及分析。結(jié)果表明,伊丹河輸水河段冰期封凍形式為平封,可形成穩(wěn)定冰蓋從而實(shí)現(xiàn)冰期冰蓋下輸水,且滿足冰期正常輸水的要求。

    冰期輸水;調(diào)水工程;冰情演變;數(shù)值模擬; 伊丹河

    河冰是寒冷地區(qū)尤其是高緯度區(qū)域冬季河道中普遍存在的自然現(xiàn)象。北方河道冰期輸水時(shí),嚴(yán)寒的氣候往往出現(xiàn)冰害,導(dǎo)致運(yùn)行和管理上的許多困難。因此研究冰期輸水所涉及的相關(guān)冰害問題對(duì)于河渠道冬季安全運(yùn)行具有重要意義。

    河冰生消演變?yōu)樗畡?dòng)力、機(jī)械力和熱力等相互耦合作用且極為復(fù)雜的物理過程。多年來,隨著河冰水力學(xué)理論的發(fā)展,國(guó)內(nèi)外學(xué)者們相繼提出了一些模擬河道冰情的數(shù)值模型。例如Lal等[1]根據(jù)熱力和兩層輸冰模型的理論分析框架,建立了冰情演變河冰數(shù)值模型—RICE模型。Shen等[2]提出了二維與一維結(jié)合的河冰動(dòng)力學(xué)模型DynaRICE,主要用于模擬冰壩形成的動(dòng)力學(xué)過程等。Zufelt等[3]建立了一維冰水耦合的冰塞動(dòng)力學(xué)模型,模擬封河期和開河期水位變化。茅澤育等[4-5]針對(duì)天然河道彎曲復(fù)雜的特點(diǎn),首次建立了適體坐標(biāo)下的二維河冰數(shù)值模型,對(duì)黃河河曲段、黃河萬(wàn)家寨水庫(kù)庫(kù)區(qū)河段、松花江上游白山河段、新疆額爾齊斯河段等進(jìn)行了數(shù)值模擬研究,取得了較好的效果。Fu等[6]在RICE和RICEN模型的基礎(chǔ)上,開發(fā)了用于模擬冰期洪水的YRIDM模型。靳國(guó)厚等[7]采用一維非恒定流水力學(xué)及熱力學(xué)相耦合的模型,預(yù)報(bào)了京密引水渠、引黃濟(jì)青等輸水工程的冰期輸水過程。卞雪軍等[8]以多元線性回歸理論與人工神經(jīng)網(wǎng)絡(luò)理論為基礎(chǔ),設(shè)計(jì)開發(fā)了黃河內(nèi)蒙古河段冰情預(yù)報(bào)系統(tǒng)。

    如上所述,河流冰情演變受水力、氣象、熱力、河勢(shì)及河道特征等眾多因素的綜合影響,這些影響因素的不同組合必將導(dǎo)致發(fā)生不同的冰情。所以,對(duì)于任一具體河道,需根據(jù)具體的氣候、地形、水力等條件建立模型、率定參數(shù)并進(jìn)行數(shù)值計(jì)算。

    伊丹河全流域面積為483 km2,其河源至新立城水庫(kù)入庫(kù)斷面全長(zhǎng)47.71 km(圖1)。伊丹河為伊通河的最大支流,沿途流經(jīng)伊通縣二道鎮(zhèn)、伊通縣伊丹鎮(zhèn)、長(zhǎng)春市南關(guān)區(qū)新湖鎮(zhèn)等,在長(zhǎng)春市南關(guān)區(qū)新湖鎮(zhèn)榆樹屯注入伊通河。

    圖1 伊丹河輸水河段示意圖

    伊丹河冰期輸水工程為吉林省中部城市群引松供水長(zhǎng)春干線工程的重要組成部分。長(zhǎng)春干線全長(zhǎng)14.1 km,設(shè)計(jì)流量13.2 m3/s,由豐滿水庫(kù)經(jīng)110 km有壓隧洞(隧洞取水口位于水庫(kù)正常蓄水位以下21 m),經(jīng)馮家?guī)X分水樞紐進(jìn)入長(zhǎng)春干線,再在C09斷面處泄入伊丹河,并自流入新立城水庫(kù)。匯入斷面C09至新立城水庫(kù)入庫(kù)斷面河段長(zhǎng)度為7.3 m。由于工程地處嚴(yán)寒地區(qū),河道輸水一旦形成冰害,將抬高長(zhǎng)春干線末端即出口斷面C09的水位,減小長(zhǎng)春干線的輸水作用水頭,影響長(zhǎng)春干線設(shè)計(jì)輸水量。由于長(zhǎng)春干線與馮家?guī)X分水樞紐連通,還將影響四平泵站、遼源泵站設(shè)計(jì)參數(shù),并對(duì)馮家?guī)X調(diào)壓井結(jié)構(gòu)產(chǎn)生不利影響,進(jìn)而影響整個(gè)總干線的設(shè)計(jì)參數(shù)。此外,長(zhǎng)春干線伊丹河下游末端為水庫(kù),極易形成冰害,一旦形成冰害,對(duì)河岸、堤防產(chǎn)生影響,影響輸水穩(wěn)定性。因此,實(shí)現(xiàn)長(zhǎng)春干線伊丹河段冰期安全輸水對(duì)于長(zhǎng)春干線、馮家?guī)X分水樞紐乃至整個(gè)中部供水工程正常運(yùn)行至關(guān)重要。本文根據(jù)伊丹河河道的地形、氣象及研究河段實(shí)際情況,建立河冰數(shù)值模型,并應(yīng)用伊丹河1958—1988年的水文、氣象資料對(duì)伊丹河冰期輸水進(jìn)行了冰情數(shù)值模擬。

    1 河冰數(shù)學(xué)模型及計(jì)算方法

    封凍期河冰演變的模型主要包括河道水力學(xué)模型、輸冰模型、熱力學(xué)模型及冰凍模型等。

    封凍河道的水流一維流動(dòng)近似由以下連續(xù)性方程和動(dòng)量方程描述[4]:

    (1)

    式中:Q為流量;A為過流斷面面積;ql為單寬側(cè)向入流量;x為距離;t為時(shí)間;g為重力加速度;z為水位;zb為河床高程;ti為冰蓋厚度;ρ、ρi分別為水和冰的密度;h為冰蓋下水深;nc為綜合糙率系數(shù);P為濕周。

    根據(jù)熱量守恒原理,沿河水溫η的時(shí)空分布近似由以下對(duì)流-擴(kuò)散方程表述[4]:

    (3)

    其中η=ρCpATw

    式中:v為斷面平均流速;B0為水面寬;∑S為水體與周圍環(huán)境的單位面積熱交換量(增熱或失熱);Cp為水的比熱容,取4 185 J/(kg·℃);Tw為斷面平均水溫。

    根據(jù)分層輸冰理論,河流中的流冰由面冰和懸浮冰兩部分組成,其質(zhì)量守恒方程可參見文獻(xiàn)[9]。根據(jù)冰蓋前緣斷面的水動(dòng)力條件,以及上游來冰情況,冰蓋的發(fā)展一般有并置推進(jìn)、水力增厚推進(jìn)、機(jī)械增厚推進(jìn)[1,10]。冰蓋下冰的輸移和積聚采用文獻(xiàn)[11]中的輸冰能力公式計(jì)算。當(dāng)積聚冰蓋體向上游發(fā)展過程中,一方面,冰蓋體中的孔隙水凍結(jié)而形成固狀冰體外殼,增大冰體強(qiáng)度;另一方面,由于冰蓋表面和底面發(fā)生熱交換,冰體厚度將發(fā)生熱力增厚(或消融)。冰蓋體的熱力增厚或消融過程由熱力平衡原理得到,即:

    (4)

    式中:ti為水溫度;Ts為冰層表面溫度;Ta為氣溫;Tm為冰層熱融度,即0℃;hia為大氣與冰層表面的熱交換系數(shù),取19.71;hwi為水體與冰層底面的熱交換系數(shù)。

    模型中式(1)(2)(3)采用顯式有限差分MacCormack步進(jìn)格式求解[12],其在時(shí)間和空間上具有二階精度o(Δx2,Δt2)。MacCormack方法屬兩步預(yù)測(cè)-修正型,其預(yù)測(cè)和修正步驟,可以交替地向前和向后差分。各子模型組成的聯(lián)合方程組采用同步求解。為準(zhǔn)確模擬伊丹河輸水河段的冰情,筆者曾應(yīng)用白山河段完整詳盡的氣象、冰情原型觀測(cè)資料,采用工程類比驗(yàn)證方法對(duì)上述數(shù)學(xué)模型進(jìn)行了模型參數(shù)率定及驗(yàn)證[13]。

    自豐滿水庫(kù)取水口至長(zhǎng)春干線出水口斷面C09之間的埋地輸水管道水溫,采用以下一維水溫對(duì)流-擴(kuò)散方程計(jì)算[14],得到長(zhǎng)春干線出水口斷面水溫隨時(shí)間的變化過程:

    (5)

    2 伊丹河輸水河段數(shù)值模擬

    2.1 河段基本概況

    對(duì)伊丹河(全長(zhǎng)47.7 km)上游河源、下游新立城水庫(kù)及豐滿水庫(kù)有壓隧洞取水口作全系統(tǒng)研究。長(zhǎng)春干線全長(zhǎng)14.1 km(PCCP管道,平均埋深6 m),匯入斷面C09;豐滿水庫(kù)至長(zhǎng)春干線的有壓隧洞全長(zhǎng)110 km。計(jì)算域共布置斷面154個(gè),平均斷面間距約300 m。斷面1布置在伊丹河河源斷面(C09斷面),編號(hào)為84,距河源40.40 km;新立城水庫(kù)入庫(kù)斷面編號(hào)為154,距河源47.7 km;斷面127位于河源下游44.93 km。河段平均坡降3.2×10-4。

    2.2 氣象水文資料選取與整理

    采用新立城壩上水文站、伊通氣象站1958—1988年水文、氣象多年平均觀測(cè)數(shù)據(jù)。伊通氣象站冬季(11月至次年3月)多年日均氣溫如圖2所示(橫坐標(biāo)時(shí)間以11月1日為起點(diǎn)(下同))。計(jì)算總時(shí)長(zhǎng)151 d,時(shí)間步長(zhǎng)60 s,對(duì)于實(shí)測(cè)斷面間距過大河段進(jìn)行內(nèi)插,最大計(jì)算斷面間距取50 m。各實(shí)測(cè)橫斷面的過水面積與濕周通過微斷面法[15]進(jìn)行計(jì)算。

    圖2 伊通氣象站冬季多年平均日均氣溫

    2.3 初始條件及邊界條件確定

    初始條件為計(jì)算開始時(shí)刻計(jì)算域內(nèi)流量、水位、水溫、浮冰濃度、水內(nèi)冰濃度和冰蓋體厚度。伊丹河河道初始流量及水位采用上游邊界及下游邊界對(duì)應(yīng)的流量和水位值,通過長(zhǎng)時(shí)間迭代運(yùn)算直至穩(wěn)定,最終得到沿程流量和水位值;初始水體溫度采用6℃。根據(jù)水文、氣象資料,初始時(shí)刻沿程浮冰濃度、水內(nèi)冰濃度分布及冰蓋體厚度均設(shè)定為零。流場(chǎng)計(jì)算邊界條件:上游邊界為伊丹河河源斷面處的流量過程以及中途匯入的長(zhǎng)春干線設(shè)計(jì)引水流量(13.2 m3/s),下游邊界為入庫(kù)斷面(斷面154)流量水位關(guān)系曲線(圖3)。水流溫度場(chǎng)計(jì)算的邊界條件為伊丹河河源斷面處以及豐滿水庫(kù)取水口處實(shí)測(cè)水溫值(圖4),初始條件通過模型迭代計(jì)算,直至初始溫度場(chǎng)基本達(dá)到熱交換的準(zhǔn)平衡狀態(tài)得到。

    圖3 新立城入庫(kù)斷面冬季流量水位關(guān)系線

    圖4 豐滿水庫(kù)有壓隧洞取水口處實(shí)測(cè)水溫隨時(shí)間變化

    2.4 計(jì)算結(jié)果及分析

    圖5為封河過程中不同時(shí)刻輸水河段沿程水位計(jì)算結(jié)果,圖中橫坐標(biāo)以C09斷面為起點(diǎn)(下同)。隨著氣溫降低,冰蓋前緣逐漸向上游推進(jìn),前緣斷面附近水位隨之壅高,同時(shí)冰體厚度增長(zhǎng)。1月29日水位壅高達(dá)到最大;3月下旬隨著氣溫回暖、水溫升高,冰體開始消融,沿程水位降低。圖6給出的是C09斷面封凍期(11月至次年3月)水位變化過程。由圖6可知,12月末1月初,由于下游河段冰蓋的生成,C09斷面水位開始壅高,最高壅水達(dá)0.16 m。隨著冰蓋進(jìn)一步固結(jié),冰蓋糙率系數(shù)降低,計(jì)算壅水位逐漸降低。入庫(kù)斷面封凍期流量變化過程如圖7所示,12月初,入庫(kù)斷面流量有所減小,這是由于河道槽蓄水量增加所致,隨后入庫(kù)斷面流量穩(wěn)定在13.4 m3/s。入庫(kù)斷面水溫隨時(shí)間的變化過程如圖8所示。由于受長(zhǎng)春干線輸水水體水溫的影響,伊丹河輸水河段的水溫變化滯后于日均氣溫變化。從圖8可知,12月15日受大氣溫度下降及河道水體與大氣熱交換的影響,水溫持續(xù)降至0℃以下,水體進(jìn)入過冷卻狀態(tài),隨著冰層生成,水溫又逐漸穩(wěn)定在0℃附近。

    圖5 不同時(shí)間沿程水位變化

    圖6 C09斷面冬季水位過程

    圖7 入庫(kù)斷面冬季流量過程

    圖8 入庫(kù)斷面水溫變化過程

    圖9 1月29日河道弗勞德數(shù)沿程分布

    圖10 冰蓋前緣位置隨時(shí)間變化

    圖11 入庫(kù)斷面冰厚隨時(shí)間變化

    圖12 1月29日冰厚沿程分布

    1月29日輸水河段沿程弗勞德數(shù)如圖9所示。由圖9可知,輸水河段的弗勞德數(shù)均在第二臨界弗勞德數(shù)(0.09)以下,表明在氣象及水溫等條件下,河道大部分河段都可形成冰蓋。圖10給出的是冰蓋推進(jìn)前緣位置隨時(shí)間的變化,可以看出,12月30日下游入庫(kù)斷面首先形成冰蓋,6 d后冰蓋前緣到達(dá)127斷面,由于水溫的影響,冰蓋前緣到達(dá)127斷面之后便不在向前發(fā)展。圖11給出的是封凍期入庫(kù)斷面冰蓋平均厚度隨時(shí)間的變化,可見入庫(kù)斷面的初封時(shí)間為12月30日。隨后,冰蓋體在冰面與大氣及冰底面與水體熱交換的共同作用下,開始熱力增厚。2月20日入庫(kù)斷面冰厚達(dá)最大,為0.86 m,此后受氣溫回暖影響,冰厚開始減小。圖12為伊丹河輸水河段1月29日冰蓋體厚度沿程分布。

    3 結(jié) 論

    a. 伊丹河輸水河段冰期沿程弗勞德數(shù)均在第二臨界弗勞德數(shù)(0.09)以下,在氣象以及水溫條件等適宜條件下,河道大部分河段都可形成冰蓋,進(jìn)而實(shí)現(xiàn)河道冰期冰蓋下輸水,輸水河段沿程弗勞德數(shù)大多小于0.06,河段封凍形式為平封形式。

    b. 河段封凍始于入庫(kù)斷面,然后逐漸上溯,6 d之后冰蓋前緣到達(dá)127斷面,冰期最大冰厚為0.86 m。

    c. 在封凍過程中,冰蓋前緣推進(jìn)速度及冰體厚度分布受諸多因素影響,且隨冰蓋前緣逐步向上游推進(jìn),沿程水位也逐漸壅高,最高壅水達(dá)0.16 m。冰期入庫(kù)斷面流量值在冰蓋形成期間略有下降(河道槽蓄水量增加所致),隨后一直維持在13.4 m3/s附近,伊丹河輸水河段能夠滿足冰期設(shè)計(jì)流量下的正常輸水要求。

    [1] LAL A M,SHEN H T.A mathematical model for river ice processes[J].Journal of Hydraulic Engineering,1993,117(7): 851-867.

    [2] SHEN H T,SU J,LIU L.SPH simulation of river ice dynamics[J].Journal of Computational Physics,2000,165(2): 752-770.

    [3] ZUFELT J E,ETTEMA R.Fully coupled model of ice-jam dynamics[J].Journal of Cold Regions Engineering,2000,14(1): 24-41.

    [4] 茅澤育,吳劍疆,張磊,等.天然河道冰塞演變發(fā)展的數(shù)值模擬[J].水科學(xué)進(jìn)展,2003,14(6): 700-705.(MAO Zeyu,WU Jianjiang,ZHANG Lei,et al.Numerical simulation of river ice jam[J].Advance in Water Science,2003,14(6): 700-705.(in Chinese))

    [5] 茅澤育,許昕,王愛民,等.基于適體坐標(biāo)變換的二維河冰模型[J].水科學(xué)進(jìn)展,2008,19(2):214-223.(MAO Zeyu,XU Xin,WANG Aimin,et al.2-D numerical model for river-ice processes based upon body-fitted coordinate[J].Advance in Water Science,2008,19(2):700-705.(in Chinese))

    [6] FU C,POPESCU I,WANG C,et al.Challenges in modeling ice floods on the Ningxia-Inner Mongolia reach of the Yellow River,China[J].Hydrology and Earth System Sciences Discussions,2013,10(10): 12293-12329.

    [7] 靳國(guó)厚,高霈生,呂斌秀.明渠冰情預(yù)報(bào)的數(shù)學(xué)模型[J].水利學(xué)報(bào),1997,28(10):1-9.(JIN Guohou,GAO Peisheng,LU Binxiu.Ice regime forecast model in open channel[J].Journal of Hydraulic Engineering,1997,28(10):1-9.(in Chinese))

    [8] 卞雪軍,冀鴻蘭,姜新華,等.黃河內(nèi)蒙古段冰情預(yù)報(bào)系統(tǒng)的開發(fā)[J].水利水電科技進(jìn)展,2014,34(4):62-65.(BIAN Xuejun,JI Honglan,JIANG Xinhua,et al.Development and application analyze of forecasting system of ice conditions in Inner Mongolia reach of Yellow River[J].Advances in Science and Technology of Water Resources,2014,34(4):62-65.(in Chinese))

    [9] 王永填.河冰數(shù)值模型及河渠冬季輸水冰情特性研究[D].北京:清華大學(xué),1999.

    [10] PARISET E,HAUSSER R.Formation and evolution of ice covers on rivers[J].Transaction of Engineering Institute of Canada,1961,5(1): 41-49.

    [11] SHEN H T,WANG D S,LAL A M.Numerical simulation of river ice processes[J].Journal of Cold Regions Engineering,1995,9(3): 107-118.

    [12] MACCORMACK R W.A numerical method for solving the equations of compressible viscous flow[J].AIAA Journal,1982,20(9): 1275-1281.

    [13] 樊霖,茅澤育,吳劍疆,等.松花江白山河段冰情數(shù)值計(jì)算及分析[J].水科學(xué)進(jìn)展,2016,27(6):890-897.(FAN Lin,MAO Zeyu,WU Jianjiang,et al.Study on the numerical model of river ice for the Baishan reach of Songhua River,China[J].Advance in Water Science,2016,27(6):890-897.(in Chinese))

    [14] 崔慧,吳長(zhǎng)春.熱油管道非穩(wěn)態(tài)工況傳熱與流動(dòng)的耦合計(jì)算模型[J].中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2005,29(3):101-105.(CUI Hui,WU Changchun.Heat transfer and flow coupling calculation model of transient scenario for hot oil pipeline[J].Journal of China University of Petroleum,2005,29(3):101-105.(in Chinese))

    [15] FAN L,MAO Z Y,BAO J,et al.A river-ice model for Wanjiazhai reach of the Yellow River[C]//Proceedings of the 22nd IAHR International Symposium on Ice.Singapore:NUS,2014:869-876.

    Numericalsimulationofriver-iceprocessesinawater-transferriverofYidanRiverduringfreezingperiod//

    FAN Lin1, MAO Zeyu1, QI Wenbiao2, LIU Shufeng2

    (1.DepartmentofHydraulicEngineering,TsinghuaUniversity,Beijing100084,China; 2.JilinInvestigationandDesignInstituteofWaterResourcesandHydropower,Changchun130021,China)

    Based upon the theory of hydrodynamics, river mechanics, thermodynamics, river-ice hydraulics and mechanics of solids etc., a mathematic model of the river-ice processes for the Yidan reach is constructed. 30-year long-term series comprehensive field observation data, including geomorphy, meteorology, hydrology, thermo-dynamics and ice regime, are used to simulate the water transfer process of the Yidan River during freezing period. The results show that steady ice cover can be formed by juxtaposition on the Yidan River, which achieves the goal of water transfer beneath the ice cover. In addition, the design discharge of the Yidan River water transfer project can be guaranteed during freezing period.

    water transfer during freezing period; water diversion project;river ice process; numerical simulation; Yidan River

    國(guó)家重點(diǎn)研發(fā)計(jì)劃(2016YFC0402504);清華大學(xué)博士生短期出國(guó)訪學(xué)基金(2013118)

    樊霖(1990—),男,博士研究生,主要從事河冰水力學(xué)研究。E-mail:l-fan11@mails.tsinghua.edu.cn

    茅澤育(1962—),男,教授,博士,主要從事水力學(xué)及河流動(dòng)力學(xué)研究。E-mail:maozeyu@tsinghua.edu.cn

    10.3880/j.issn.1006-7647.2017.06.003

    TV143

    A

    1006-7647(2017)06-0014-05

    2017-02-21 編輯:鄭孝宇)

    猜你喜歡
    冰情冰蓋河段
    軍事文摘(2024年6期)2024-04-30 03:13:59
    格陵蘭島的冰蓋悄悄融化
    參花(下)(2022年1期)2022-01-15 00:45:01
    松花江干流哈爾濱江段封凍
    SL流量計(jì)在特殊河段的應(yīng)用——以河源水文站為例
    石瀘高速公路(紅河段)正式通車
    石油瀝青(2020年1期)2020-05-25 06:54:10
    詠菊致友人
    長(zhǎng)距離輸水工程的冰期冰蓋數(shù)值模擬研究
    海冰基礎(chǔ)知識(shí)及船舶冰區(qū)航行的注意事項(xiàng)
    淺析黃河河段維修養(yǎng)護(hù)施工工藝
    河南科技(2014年14期)2014-02-27 14:11:44
    防凌期寧蒙河段流量演進(jìn)的分期分河段混合算法研究
    青春草亚洲视频在线观看| 边亲边吃奶的免费视频| 欧美日韩一区二区视频在线观看视频在线 | 国产高清三级在线| 日日撸夜夜添| 久久精品国产自在天天线| 你懂的网址亚洲精品在线观看 | 91精品一卡2卡3卡4卡| 一个人免费在线观看电影| 国产真实伦视频高清在线观看| 观看免费一级毛片| 一级毛片电影观看 | 日韩欧美在线乱码| 有码 亚洲区| 欧美极品一区二区三区四区| 久久久精品94久久精品| 麻豆精品久久久久久蜜桃| 91av网一区二区| 内地一区二区视频在线| 最近中文字幕高清免费大全6| 99久久九九国产精品国产免费| 亚洲三级黄色毛片| 久久久精品大字幕| 久久亚洲国产成人精品v| 亚洲人成网站在线观看播放| 中文字幕久久专区| 午夜视频国产福利| 久久鲁丝午夜福利片| 99久久成人亚洲精品观看| 秋霞在线观看毛片| 国产精品永久免费网站| 国产午夜精品论理片| 免费看光身美女| 大又大粗又爽又黄少妇毛片口| 日韩在线高清观看一区二区三区| 少妇高潮的动态图| av.在线天堂| 精品熟女少妇av免费看| 婷婷色麻豆天堂久久 | 免费观看性生交大片5| 日本免费a在线| 一级二级三级毛片免费看| 我要搜黄色片| 午夜福利视频1000在线观看| 国产黄a三级三级三级人| 亚洲国产精品sss在线观看| av播播在线观看一区| 午夜福利在线在线| 美女大奶头视频| 七月丁香在线播放| 激情 狠狠 欧美| 最近中文字幕高清免费大全6| 最后的刺客免费高清国语| 亚洲欧美一区二区三区国产| 成年av动漫网址| 亚洲美女搞黄在线观看| 亚洲丝袜综合中文字幕| 亚洲,欧美,日韩| 少妇熟女aⅴ在线视频| 久久久精品94久久精品| 亚洲怡红院男人天堂| 免费不卡的大黄色大毛片视频在线观看 | 日韩强制内射视频| 国产高清有码在线观看视频| 好男人在线观看高清免费视频| 男人狂女人下面高潮的视频| 亚洲人成网站高清观看| 好男人视频免费观看在线| 麻豆成人午夜福利视频| 欧美另类亚洲清纯唯美| 欧美日韩国产亚洲二区| 色尼玛亚洲综合影院| 国产伦一二天堂av在线观看| 久久久午夜欧美精品| 国产一区二区在线av高清观看| 最近手机中文字幕大全| 久久人妻av系列| 国产白丝娇喘喷水9色精品| 日本色播在线视频| 亚洲欧美日韩卡通动漫| 内地一区二区视频在线| 成人亚洲欧美一区二区av| 亚洲成人中文字幕在线播放| av在线播放精品| 韩国高清视频一区二区三区| 久久久久久久久久成人| 深夜a级毛片| 激情 狠狠 欧美| 午夜激情欧美在线| 中国美白少妇内射xxxbb| 国产成人免费观看mmmm| 听说在线观看完整版免费高清| 国产老妇伦熟女老妇高清| 国产综合懂色| 亚洲无线观看免费| 日本黄色视频三级网站网址| 热99在线观看视频| 最近视频中文字幕2019在线8| 国产成人a∨麻豆精品| 插阴视频在线观看视频| 晚上一个人看的免费电影| 看片在线看免费视频| 岛国毛片在线播放| 高清日韩中文字幕在线| 国产精品嫩草影院av在线观看| 日本一本二区三区精品| 久久热精品热| 午夜亚洲福利在线播放| 午夜福利成人在线免费观看| 日韩国内少妇激情av| 中文欧美无线码| 乱人视频在线观看| 如何舔出高潮| 亚洲色图av天堂| 久久久久久九九精品二区国产| 少妇丰满av| 伦精品一区二区三区| 午夜福利高清视频| 国产成人91sexporn| www日本黄色视频网| 99在线人妻在线中文字幕| 舔av片在线| 日日干狠狠操夜夜爽| 国产成人a∨麻豆精品| 欧美xxxx性猛交bbbb| 一本一本综合久久| 免费搜索国产男女视频| 国产老妇伦熟女老妇高清| 久久久午夜欧美精品| 精品午夜福利在线看| 久久6这里有精品| 国产欧美日韩精品一区二区| 精品久久久久久久久久久久久| 亚洲精华国产精华液的使用体验| 亚洲精品日韩在线中文字幕| 国产真实伦视频高清在线观看| 精品不卡国产一区二区三区| 91aial.com中文字幕在线观看| 免费不卡的大黄色大毛片视频在线观看 | 精品酒店卫生间| 日韩人妻高清精品专区| 日韩三级伦理在线观看| 99国产精品一区二区蜜桃av| 成人亚洲精品av一区二区| 亚洲人成网站高清观看| 亚洲精华国产精华液的使用体验| 婷婷色麻豆天堂久久 | 久久精品夜色国产| 在线天堂最新版资源| 免费观看人在逋| 最后的刺客免费高清国语| 久久国内精品自在自线图片| 九九爱精品视频在线观看| 99热这里只有是精品50| 小说图片视频综合网站| 最近最新中文字幕免费大全7| 欧美性感艳星| 非洲黑人性xxxx精品又粗又长| 国产精品美女特级片免费视频播放器| 直男gayav资源| 国产大屁股一区二区在线视频| 亚洲在线观看片| 国产精品1区2区在线观看.| 亚洲真实伦在线观看| 少妇的逼好多水| 能在线免费看毛片的网站| 蜜桃亚洲精品一区二区三区| 一个人看视频在线观看www免费| 国产乱人视频| 美女黄网站色视频| 综合色av麻豆| 非洲黑人性xxxx精品又粗又长| 啦啦啦观看免费观看视频高清| 色网站视频免费| 日韩欧美精品v在线| 国产精品久久久久久久久免| 女的被弄到高潮叫床怎么办| 夜夜看夜夜爽夜夜摸| 久久精品影院6| 最后的刺客免费高清国语| 99热这里只有是精品50| 亚洲欧美精品专区久久| or卡值多少钱| 三级毛片av免费| 日韩一区二区视频免费看| 久久精品国产亚洲av涩爱| 少妇的逼水好多| 国产av码专区亚洲av| 精品一区二区免费观看| 国产精品日韩av在线免费观看| 亚洲综合精品二区| 国产真实乱freesex| 亚洲av成人精品一区久久| 三级经典国产精品| 久久精品国产自在天天线| 视频中文字幕在线观看| 九九爱精品视频在线观看| 亚洲最大成人中文| 只有这里有精品99| 97在线视频观看| 成人性生交大片免费视频hd| 18禁动态无遮挡网站| 国产欧美日韩精品一区二区| 日日啪夜夜撸| 美女国产视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 精品无人区乱码1区二区| 最近中文字幕2019免费版| 欧美日韩在线观看h| 性插视频无遮挡在线免费观看| 欧美潮喷喷水| 亚洲欧美精品自产自拍| 欧美性感艳星| 成人综合一区亚洲| 亚洲精品成人久久久久久| 一级毛片我不卡| 天堂av国产一区二区熟女人妻| 日韩欧美精品v在线| 日韩欧美精品免费久久| 色5月婷婷丁香| 99热网站在线观看| 熟女人妻精品中文字幕| 三级经典国产精品| 男女边吃奶边做爰视频| 久久久久久久久久久丰满| 一区二区三区高清视频在线| 观看美女的网站| 亚洲天堂国产精品一区在线| 亚洲精品乱码久久久久久按摩| 嫩草影院新地址| 亚洲在线观看片| 少妇人妻一区二区三区视频| 国产精品美女特级片免费视频播放器| 精品国产一区二区三区久久久樱花 | 搡女人真爽免费视频火全软件| 午夜精品在线福利| 国产伦一二天堂av在线观看| 一边亲一边摸免费视频| 级片在线观看| 精华霜和精华液先用哪个| 久久久久久九九精品二区国产| 狠狠狠狠99中文字幕| videos熟女内射| 亚洲精品亚洲一区二区| 亚洲va在线va天堂va国产| 国产淫语在线视频| 永久网站在线| 亚洲av男天堂| 久久久久九九精品影院| 特大巨黑吊av在线直播| 国产亚洲精品av在线| 国产精品一区二区三区四区免费观看| 亚洲国产高清在线一区二区三| 亚洲精品日韩av片在线观看| 亚洲国产欧洲综合997久久,| 亚洲国产最新在线播放| 别揉我奶头 嗯啊视频| 看黄色毛片网站| 99视频精品全部免费 在线| 免费看a级黄色片| av免费在线看不卡| 校园人妻丝袜中文字幕| 日韩亚洲欧美综合| 人妻系列 视频| 韩国av在线不卡| 少妇熟女aⅴ在线视频| 长腿黑丝高跟| 国产av一区在线观看免费| 久久精品人妻少妇| 日本一二三区视频观看| 国产精品一区二区三区四区免费观看| 美女xxoo啪啪120秒动态图| 中文在线观看免费www的网站| 久久久午夜欧美精品| 一区二区三区乱码不卡18| 国产精品1区2区在线观看.| 97超视频在线观看视频| 国产91av在线免费观看| 亚洲欧美日韩无卡精品| 国产精品美女特级片免费视频播放器| 亚洲第一区二区三区不卡| 我要搜黄色片| 日韩在线高清观看一区二区三区| АⅤ资源中文在线天堂| 国产乱人偷精品视频| 久久99蜜桃精品久久| 99久久精品国产国产毛片| 看片在线看免费视频| 国产成人免费观看mmmm| 日本三级黄在线观看| 精品人妻视频免费看| 国产亚洲精品久久久com| 看免费成人av毛片| 激情 狠狠 欧美| 久久久精品大字幕| 国产伦精品一区二区三区视频9| 在线免费观看不下载黄p国产| 波野结衣二区三区在线| 欧美性感艳星| 免费av观看视频| 免费电影在线观看免费观看| 超碰97精品在线观看| 永久网站在线| 男人的好看免费观看在线视频| 天天一区二区日本电影三级| 超碰av人人做人人爽久久| 美女被艹到高潮喷水动态| 日韩精品有码人妻一区| 成人鲁丝片一二三区免费| 亚洲av日韩在线播放| 亚洲精品影视一区二区三区av| 高清视频免费观看一区二区 | av视频在线观看入口| 午夜免费男女啪啪视频观看| 国产精品一区二区三区四区久久| 别揉我奶头 嗯啊视频| 91aial.com中文字幕在线观看| 成人性生交大片免费视频hd| 国产免费一级a男人的天堂| 六月丁香七月| 国产黄片美女视频| 狠狠狠狠99中文字幕| 亚洲欧美成人精品一区二区| 成年女人永久免费观看视频| 成人欧美大片| 亚洲av福利一区| 蜜桃亚洲精品一区二区三区| 国产真实伦视频高清在线观看| 国产精品久久久久久精品电影小说 | 久久久久久久久久黄片| 建设人人有责人人尽责人人享有的 | 日韩大片免费观看网站 | 日韩人妻高清精品专区| 国产精品一区www在线观看| 高清午夜精品一区二区三区| 国语对白做爰xxxⅹ性视频网站| 日本熟妇午夜| 国产亚洲最大av| 久久久精品大字幕| 亚洲熟妇中文字幕五十中出| 天堂网av新在线| 高清日韩中文字幕在线| 乱码一卡2卡4卡精品| 欧美日韩精品成人综合77777| 国产视频内射| 欧美又色又爽又黄视频| 麻豆乱淫一区二区| 国产免费一级a男人的天堂| 一本久久精品| 国产老妇女一区| 女人十人毛片免费观看3o分钟| 欧美日本视频| 狂野欧美激情性xxxx在线观看| 丝袜喷水一区| 1000部很黄的大片| 国产在视频线精品| 夜夜看夜夜爽夜夜摸| 嫩草影院入口| 特级一级黄色大片| 99九九线精品视频在线观看视频| 亚洲av中文字字幕乱码综合| 久久久精品94久久精品| 国产又色又爽无遮挡免| 国产一区二区亚洲精品在线观看| 成人一区二区视频在线观看| 国产精品一区二区三区四区久久| 国产在视频线在精品| 亚洲av成人精品一二三区| 久久精品夜色国产| 男人舔奶头视频| 国产淫片久久久久久久久| 最近中文字幕2019免费版| 超碰av人人做人人爽久久| 欧美日本视频| 一级毛片久久久久久久久女| 中国美白少妇内射xxxbb| 一级毛片电影观看 | 亚洲怡红院男人天堂| 国产91av在线免费观看| 精品国产三级普通话版| 一卡2卡三卡四卡精品乱码亚洲| 2021天堂中文幕一二区在线观| 国产片特级美女逼逼视频| 中文亚洲av片在线观看爽| 六月丁香七月| 精品不卡国产一区二区三区| 亚洲婷婷狠狠爱综合网| 亚洲,欧美,日韩| 偷拍熟女少妇极品色| 国产成人午夜福利电影在线观看| 色哟哟·www| 啦啦啦观看免费观看视频高清| 中文亚洲av片在线观看爽| 亚洲aⅴ乱码一区二区在线播放| 日日摸夜夜添夜夜添av毛片| 国产一区二区亚洲精品在线观看| 国产激情偷乱视频一区二区| 欧美一区二区国产精品久久精品| 99热这里只有是精品50| 又粗又硬又长又爽又黄的视频| 丝袜喷水一区| 中文精品一卡2卡3卡4更新| 国产免费又黄又爽又色| 国产黄a三级三级三级人| 联通29元200g的流量卡| 国产精品三级大全| 成人午夜精彩视频在线观看| 国国产精品蜜臀av免费| 人人妻人人看人人澡| 亚洲国产欧美人成| 日韩成人伦理影院| 国产精品女同一区二区软件| 欧美+日韩+精品| 免费大片18禁| 秋霞伦理黄片| 少妇被粗大猛烈的视频| 精品一区二区三区视频在线| 免费av观看视频| eeuss影院久久| 又爽又黄a免费视频| 亚洲精品乱码久久久v下载方式| 亚洲中文字幕一区二区三区有码在线看| 亚洲欧美日韩卡通动漫| 赤兔流量卡办理| 亚洲国产精品专区欧美| 少妇熟女欧美另类| 99久久九九国产精品国产免费| 国产精品一区二区三区四区久久| 亚洲一区高清亚洲精品| 日本免费在线观看一区| 国产精品三级大全| 成人亚洲精品av一区二区| av在线老鸭窝| 天堂中文最新版在线下载 | 精品久久久久久久末码| 成人av在线播放网站| 欧美另类亚洲清纯唯美| 国产极品精品免费视频能看的| 亚洲aⅴ乱码一区二区在线播放| 国产精品人妻久久久久久| 久久久国产成人免费| 国产亚洲av嫩草精品影院| 日本一本二区三区精品| 男女啪啪激烈高潮av片| 秋霞在线观看毛片| 免费观看a级毛片全部| 精品99又大又爽又粗少妇毛片| 26uuu在线亚洲综合色| 特大巨黑吊av在线直播| 精品久久久久久电影网 | 日本免费一区二区三区高清不卡| 国产视频内射| 午夜a级毛片| 国产大屁股一区二区在线视频| 亚洲18禁久久av| 青春草视频在线免费观看| 欧美精品国产亚洲| 亚洲欧美成人精品一区二区| 欧美一区二区亚洲| 久久亚洲国产成人精品v| 久久人人爽人人爽人人片va| 精品久久久噜噜| av又黄又爽大尺度在线免费看 | 亚洲在久久综合| 91精品伊人久久大香线蕉| 三级经典国产精品| 男人的好看免费观看在线视频| 午夜福利网站1000一区二区三区| 午夜激情欧美在线| 国产真实伦视频高清在线观看| 亚洲不卡免费看| 美女黄网站色视频| 看黄色毛片网站| 免费观看性生交大片5| 三级男女做爰猛烈吃奶摸视频| 在现免费观看毛片| 国产精品综合久久久久久久免费| 久久精品久久久久久久性| 色噜噜av男人的天堂激情| 国产又黄又爽又无遮挡在线| 97在线视频观看| 伊人久久精品亚洲午夜| 最后的刺客免费高清国语| 1000部很黄的大片| 日本黄色片子视频| 亚洲欧洲国产日韩| 十八禁国产超污无遮挡网站| 一级毛片电影观看 | 高清av免费在线| 亚洲欧美成人精品一区二区| 91久久精品国产一区二区成人| 精品国产一区二区三区久久久樱花 | 啦啦啦啦在线视频资源| 婷婷六月久久综合丁香| 免费人成在线观看视频色| 偷拍熟女少妇极品色| 一个人看的www免费观看视频| 禁无遮挡网站| 日日啪夜夜撸| 高清毛片免费看| 久久人人爽人人爽人人片va| 麻豆一二三区av精品| 一边亲一边摸免费视频| 不卡视频在线观看欧美| 国模一区二区三区四区视频| 永久网站在线| 超碰97精品在线观看| 色综合亚洲欧美另类图片| 欧美潮喷喷水| av在线播放精品| 亚洲第一区二区三区不卡| 少妇熟女欧美另类| 国产精品av视频在线免费观看| 亚洲人成网站高清观看| 99九九线精品视频在线观看视频| 在线a可以看的网站| 成人av在线播放网站| 91av网一区二区| 最近最新中文字幕大全电影3| 国产成人福利小说| 一个人看视频在线观看www免费| 一区二区三区乱码不卡18| 亚洲精品自拍成人| 蜜桃久久精品国产亚洲av| 亚洲精品色激情综合| 尾随美女入室| 亚洲高清免费不卡视频| 日韩av不卡免费在线播放| 18禁动态无遮挡网站| 亚洲av福利一区| 国产精品不卡视频一区二区| 国产免费视频播放在线视频 | 国产成人免费观看mmmm| 午夜免费男女啪啪视频观看| 嫩草影院入口| 午夜日本视频在线| 午夜福利成人在线免费观看| 九色成人免费人妻av| 亚洲最大成人中文| 美女大奶头视频| 亚洲天堂国产精品一区在线| 舔av片在线| 大香蕉久久网| 国产精品嫩草影院av在线观看| 日韩欧美精品v在线| 老女人水多毛片| 国产成年人精品一区二区| 男女那种视频在线观看| 天堂√8在线中文| 亚洲三级黄色毛片| 丰满乱子伦码专区| 日韩成人伦理影院| 久久99热这里只频精品6学生 | 久久精品国产自在天天线| 免费黄网站久久成人精品| 亚洲伊人久久精品综合 | 日本免费在线观看一区| 国产精品熟女久久久久浪| 身体一侧抽搐| 最近的中文字幕免费完整| 91午夜精品亚洲一区二区三区| 精品无人区乱码1区二区| 日本五十路高清| 欧美丝袜亚洲另类| 午夜福利高清视频| 久久人人爽人人片av| 亚洲欧美日韩无卡精品| 午夜福利成人在线免费观看| eeuss影院久久| 国产乱人偷精品视频| 免费不卡的大黄色大毛片视频在线观看 | 久久国内精品自在自线图片| 午夜福利网站1000一区二区三区| 午夜激情福利司机影院| 日韩欧美精品v在线| 亚洲人成网站在线播| 久久久久久久久大av| 夜夜爽夜夜爽视频| 色哟哟·www| av在线亚洲专区| 丝袜美腿在线中文| 97超碰精品成人国产| 五月玫瑰六月丁香| 亚洲av电影在线观看一区二区三区 | 日韩欧美在线乱码| 国产精品国产三级国产av玫瑰| 又爽又黄a免费视频| 一级毛片我不卡| av福利片在线观看| 18禁裸乳无遮挡免费网站照片| 身体一侧抽搐| 午夜激情欧美在线| 国产三级在线视频| 中文字幕免费在线视频6| 亚洲伊人久久精品综合 | 伊人久久精品亚洲午夜| 少妇人妻一区二区三区视频| 亚洲婷婷狠狠爱综合网| 成年女人看的毛片在线观看| 国产黄片视频在线免费观看| 九草在线视频观看| 精品久久久久久久末码| 日韩,欧美,国产一区二区三区 | 久久久精品大字幕| 噜噜噜噜噜久久久久久91| 日日摸夜夜添夜夜添av毛片| 国国产精品蜜臀av免费| 欧美日韩国产亚洲二区| 色综合站精品国产| 精品一区二区三区视频在线| 波多野结衣高清无吗| 99热这里只有精品一区| 有码 亚洲区| 国产午夜精品一二区理论片| 日日干狠狠操夜夜爽|