• <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)的分期分河段混合算法研究
    免费电影在线观看免费观看| 大香蕉久久网| 午夜视频国产福利| 国产精品电影一区二区三区| 国产激情偷乱视频一区二区| 观看免费一级毛片| 亚洲国产日韩欧美精品在线观看| 人妻制服诱惑在线中文字幕| 国产伦精品一区二区三区四那| 国产aⅴ精品一区二区三区波| 看片在线看免费视频| 伊人久久精品亚洲午夜| 变态另类成人亚洲欧美熟女| 欧美日韩综合久久久久久| 熟女电影av网| 日本免费一区二区三区高清不卡| 中文亚洲av片在线观看爽| 女同久久另类99精品国产91| 精品福利观看| 搡老熟女国产l中国老女人| 麻豆成人午夜福利视频| 床上黄色一级片| 国产精品一及| 99热网站在线观看| 免费av毛片视频| 久久人人爽人人爽人人片va| 亚洲性夜色夜夜综合| 欧美激情在线99| 久久午夜福利片| 免费看美女性在线毛片视频| 日本爱情动作片www.在线观看 | 国产精品永久免费网站| 中文资源天堂在线| 成人特级黄色片久久久久久久| .国产精品久久| 女同久久另类99精品国产91| 午夜影院日韩av| 国产精品人妻久久久影院| 色综合站精品国产| 有码 亚洲区| 国产一区二区亚洲精品在线观看| 日韩欧美免费精品| 性色avwww在线观看| 日韩欧美精品免费久久| 成人av在线播放网站| 亚洲美女搞黄在线观看 | 亚洲自偷自拍三级| 长腿黑丝高跟| 搡老妇女老女人老熟妇| 免费搜索国产男女视频| 禁无遮挡网站| 久久人人爽人人片av| 久久久久久久久久久丰满| 亚洲无线观看免费| 在线天堂最新版资源| 自拍偷自拍亚洲精品老妇| 国产精品人妻久久久影院| 女人十人毛片免费观看3o分钟| 亚洲成人久久爱视频| 夜夜夜夜夜久久久久| 麻豆av噜噜一区二区三区| 99热精品在线国产| 欧美日本视频| 你懂的网址亚洲精品在线观看 | 最近的中文字幕免费完整| 好男人在线观看高清免费视频| 国产精品一区www在线观看| 色视频www国产| 别揉我奶头 嗯啊视频| 老女人水多毛片| 免费黄网站久久成人精品| 亚洲美女搞黄在线观看 | 九九热线精品视视频播放| 少妇裸体淫交视频免费看高清| 99国产精品一区二区蜜桃av| 一本一本综合久久| 99热6这里只有精品| aaaaa片日本免费| 十八禁国产超污无遮挡网站| 97超碰精品成人国产| 亚洲成人久久性| 国产精品,欧美在线| 美女xxoo啪啪120秒动态图| 一本久久中文字幕| 老司机福利观看| 午夜视频国产福利| 国产熟女欧美一区二区| 美女大奶头视频| 亚洲av美国av| 岛国在线免费视频观看| 亚洲国产色片| 美女高潮的动态| 欧美一级a爱片免费观看看| 在线国产一区二区在线| 看十八女毛片水多多多| av福利片在线观看| 亚洲精品日韩av片在线观看| ponron亚洲| 我的女老师完整版在线观看| 99热网站在线观看| ponron亚洲| 国产亚洲精品久久久久久毛片| 亚洲精品日韩av片在线观看| or卡值多少钱| 午夜老司机福利剧场| a级毛色黄片| 级片在线观看| 国产伦精品一区二区三区四那| 精品熟女少妇av免费看| 人妻夜夜爽99麻豆av| 欧美一区二区精品小视频在线| 免费av毛片视频| 舔av片在线| 亚洲欧美日韩卡通动漫| 成人二区视频| 国产 一区 欧美 日韩| 一个人看视频在线观看www免费| 国产午夜福利久久久久久| 五月伊人婷婷丁香| 国产男人的电影天堂91| 天天躁日日操中文字幕| 麻豆av噜噜一区二区三区| av黄色大香蕉| 麻豆一二三区av精品| 亚洲中文字幕一区二区三区有码在线看| 村上凉子中文字幕在线| 长腿黑丝高跟| 少妇裸体淫交视频免费看高清| 麻豆国产av国片精品| 亚洲av不卡在线观看| 老女人水多毛片| 熟女电影av网| 日韩精品有码人妻一区| 天堂√8在线中文| 狠狠狠狠99中文字幕| 久久韩国三级中文字幕| 国产视频内射| 中文在线观看免费www的网站| 亚洲性久久影院| 婷婷精品国产亚洲av| 日韩中字成人| 精品久久久久久久久av| 18禁裸乳无遮挡免费网站照片| 麻豆精品久久久久久蜜桃| 精品久久久久久久久久免费视频| 亚洲精品一卡2卡三卡4卡5卡| 蜜臀久久99精品久久宅男| 日韩欧美三级三区| 成人无遮挡网站| 老司机福利观看| 97超碰精品成人国产| 一进一出抽搐gif免费好疼| 国产高清激情床上av| 国产高清有码在线观看视频| 国产精华一区二区三区| 插阴视频在线观看视频| 美女cb高潮喷水在线观看| 精品99又大又爽又粗少妇毛片| 亚洲av中文字字幕乱码综合| 国产精品一区二区性色av| 91av网一区二区| 国产黄a三级三级三级人| 简卡轻食公司| 一本精品99久久精品77| 日日摸夜夜添夜夜添av毛片| 又粗又爽又猛毛片免费看| 丰满人妻一区二区三区视频av| 国产av不卡久久| 国产亚洲精品久久久com| 亚洲天堂国产精品一区在线| 伊人久久精品亚洲午夜| 日本五十路高清| 亚洲内射少妇av| 午夜爱爱视频在线播放| 久久久精品大字幕| 亚洲无线观看免费| 一区二区三区高清视频在线| 亚洲精品亚洲一区二区| 国产一区二区亚洲精品在线观看| 春色校园在线视频观看| 特大巨黑吊av在线直播| 亚洲欧美日韩高清在线视频| 国产精品1区2区在线观看.| 少妇被粗大猛烈的视频| 亚洲最大成人av| 男插女下体视频免费在线播放| 欧美国产日韩亚洲一区| 色综合色国产| 少妇人妻一区二区三区视频| 久久久久久久午夜电影| 男人和女人高潮做爰伦理| 国产久久久一区二区三区| 97超碰精品成人国产| 麻豆久久精品国产亚洲av| 国产真实伦视频高清在线观看| 成人特级av手机在线观看| 国内精品宾馆在线| 天堂影院成人在线观看| 精品福利观看| 亚洲美女视频黄频| 六月丁香七月| 久久久久精品国产欧美久久久| 久久精品国产亚洲av涩爱 | 夜夜看夜夜爽夜夜摸| 精品乱码久久久久久99久播| 99国产精品一区二区蜜桃av| 国产v大片淫在线免费观看| 国内精品一区二区在线观看| 午夜福利18| 两个人视频免费观看高清| 高清毛片免费看| 亚洲色图av天堂| 久久久久久伊人网av| 国产亚洲精品久久久com| 精品一区二区三区av网在线观看| 全区人妻精品视频| 国产精品一区二区三区四区久久| 亚洲va在线va天堂va国产| 在线播放国产精品三级| 99精品在免费线老司机午夜| .国产精品久久| 亚洲av成人av| 麻豆av噜噜一区二区三区| 国产精品人妻久久久久久| 久久欧美精品欧美久久欧美| 一级毛片久久久久久久久女| 免费在线观看成人毛片| 日韩大尺度精品在线看网址| 亚洲国产高清在线一区二区三| 成人国产麻豆网| 美女内射精品一级片tv| 一本精品99久久精品77| 三级经典国产精品| 日产精品乱码卡一卡2卡三| 欧美性猛交╳xxx乱大交人| 99热这里只有是精品在线观看| 黄色一级大片看看| 一a级毛片在线观看| 亚洲av第一区精品v没综合| 高清毛片免费看| av.在线天堂| 国产亚洲91精品色在线| 少妇人妻精品综合一区二区 | 女人被狂操c到高潮| 国产av麻豆久久久久久久| 1000部很黄的大片| 91久久精品国产一区二区三区| 日本黄色片子视频| 亚洲专区国产一区二区| 国产男人的电影天堂91| 久久久色成人| 91在线观看av| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲无线观看免费| 美女大奶头视频| 亚洲av第一区精品v没综合| 久久久国产成人精品二区| 欧美人与善性xxx| 国产精品无大码| 欧美日韩乱码在线| 色av中文字幕| 成人亚洲欧美一区二区av| 国产精品国产三级国产av玫瑰| 99久国产av精品国产电影| 精品久久久久久久人妻蜜臀av| 久久精品人妻少妇| 久久韩国三级中文字幕| 级片在线观看| 高清午夜精品一区二区三区 | 久久精品综合一区二区三区| 国产黄色视频一区二区在线观看 | 毛片女人毛片| 精品久久久噜噜| 欧美日韩一区二区视频在线观看视频在线 | 久久久色成人| 级片在线观看| 日本与韩国留学比较| 丝袜美腿在线中文| 中文在线观看免费www的网站| 久久欧美精品欧美久久欧美| 99riav亚洲国产免费| 在线播放无遮挡| 亚洲av美国av| 女的被弄到高潮叫床怎么办| 亚洲成人久久性| 夜夜爽天天搞| 亚洲av中文字字幕乱码综合| 国产伦精品一区二区三区视频9| 精品无人区乱码1区二区| 国产黄片美女视频| a级毛片免费高清观看在线播放| 亚洲精品乱码久久久v下载方式| 此物有八面人人有两片| 亚洲va在线va天堂va国产| 国产探花极品一区二区| 五月玫瑰六月丁香| 精品欧美国产一区二区三| 天堂网av新在线| 大型黄色视频在线免费观看| 超碰av人人做人人爽久久| 欧美成人免费av一区二区三区| 亚洲av五月六月丁香网| 91久久精品国产一区二区三区| 国产单亲对白刺激| 熟妇人妻久久中文字幕3abv| 丰满的人妻完整版| 日韩亚洲欧美综合| 欧美人与善性xxx| 亚洲熟妇熟女久久| 久久精品国产鲁丝片午夜精品| 国产三级在线视频| 精品一区二区三区视频在线| 婷婷精品国产亚洲av| 99九九线精品视频在线观看视频| 丝袜喷水一区| 在线观看一区二区三区| 国产麻豆成人av免费视频| 中文亚洲av片在线观看爽| 日日干狠狠操夜夜爽| av专区在线播放| 亚洲综合色惰| 99久国产av精品| 亚洲av中文av极速乱| 国内精品宾馆在线| 亚洲色图av天堂| 免费看美女性在线毛片视频| 99精品在免费线老司机午夜| 久久亚洲精品不卡| 日本欧美国产在线视频| 亚洲,欧美,日韩| 床上黄色一级片| 中文在线观看免费www的网站| 欧美日韩精品成人综合77777| 99热6这里只有精品| 日本-黄色视频高清免费观看| 夜夜爽天天搞| 伦精品一区二区三区| 有码 亚洲区| av黄色大香蕉| 欧美日韩精品成人综合77777| 免费在线观看成人毛片| 我的老师免费观看完整版| 色播亚洲综合网| 成人永久免费在线观看视频| 亚洲内射少妇av| 国产 一区 欧美 日韩| 国产不卡一卡二| 秋霞在线观看毛片| 色综合色国产| 国产精品人妻久久久影院| 久久精品国产亚洲网站| 日韩国内少妇激情av| 色吧在线观看| a级毛色黄片| 亚洲av五月六月丁香网| 国产亚洲91精品色在线| 亚洲国产精品久久男人天堂| 免费在线观看影片大全网站| 搡老妇女老女人老熟妇| 国产熟女欧美一区二区| av国产免费在线观看| 一级毛片我不卡| 亚洲中文日韩欧美视频| 国产精品av视频在线免费观看| 熟女人妻精品中文字幕| 欧美激情国产日韩精品一区| 亚洲国产欧美人成| 色尼玛亚洲综合影院| 国产亚洲精品久久久久久毛片| 亚洲无线在线观看| 国产乱人偷精品视频| 搡老岳熟女国产| 99热全是精品| 欧美区成人在线视频| 精品久久国产蜜桃| 中出人妻视频一区二区| 国产高清不卡午夜福利| 亚洲性夜色夜夜综合| 国产色爽女视频免费观看| 1000部很黄的大片| 国产精品久久久久久av不卡| 搞女人的毛片| 国产一区二区激情短视频| 欧美高清性xxxxhd video| 免费看a级黄色片| 在线观看午夜福利视频| 午夜精品一区二区三区免费看| 久久精品夜夜夜夜夜久久蜜豆| 狂野欧美白嫩少妇大欣赏| 一个人看视频在线观看www免费| 国产精品一区二区三区四区久久| 中出人妻视频一区二区| 国产一级毛片七仙女欲春2| 噜噜噜噜噜久久久久久91| 中文字幕免费在线视频6| 成人av一区二区三区在线看| 69人妻影院| 观看免费一级毛片| 亚洲美女黄片视频| 日韩欧美精品免费久久| 免费av观看视频| 少妇的逼好多水| 久久久精品大字幕| 97超视频在线观看视频| 99久国产av精品| av.在线天堂| 毛片一级片免费看久久久久| 日韩成人伦理影院| 少妇人妻一区二区三区视频| 男人和女人高潮做爰伦理| h日本视频在线播放| 国产一区二区激情短视频| 欧美极品一区二区三区四区| 国产高潮美女av| videossex国产| 天堂av国产一区二区熟女人妻| 97超视频在线观看视频| 国产在视频线在精品| 波多野结衣高清无吗| 久久国内精品自在自线图片| 亚洲中文日韩欧美视频| 国产精品国产三级国产av玫瑰| 成年女人毛片免费观看观看9| 又粗又爽又猛毛片免费看| 久久人人爽人人片av| 欧美日韩在线观看h| 一卡2卡三卡四卡精品乱码亚洲| 搡女人真爽免费视频火全软件 | 免费在线观看影片大全网站| 国产视频内射| 日本 av在线| 天堂影院成人在线观看| 精品久久久久久久久av| 欧美人与善性xxx| 亚洲aⅴ乱码一区二区在线播放| 国语自产精品视频在线第100页| 一进一出抽搐gif免费好疼| 最近的中文字幕免费完整| 亚洲自拍偷在线| 国产精品久久电影中文字幕| 国产精品人妻久久久影院| 亚洲av中文字字幕乱码综合| 男女之事视频高清在线观看| 亚洲成人av在线免费| 久久午夜亚洲精品久久| 日本免费一区二区三区高清不卡| 久久精品国产99精品国产亚洲性色| 男人舔女人下体高潮全视频| 久久久久久久亚洲中文字幕| 国产精品日韩av在线免费观看| 国产伦一二天堂av在线观看| 夜夜爽天天搞| 免费在线观看影片大全网站| 亚洲一区高清亚洲精品| www.色视频.com| 久久久a久久爽久久v久久| av免费在线看不卡| 久久精品国产鲁丝片午夜精品| 18禁黄网站禁片免费观看直播| 成人二区视频| 日本与韩国留学比较| 桃色一区二区三区在线观看| 久久精品人妻少妇| 国产中年淑女户外野战色| 午夜福利在线在线| 久久人人爽人人片av| 自拍偷自拍亚洲精品老妇| 欧美丝袜亚洲另类| 亚洲精品日韩在线中文字幕 | 啦啦啦啦在线视频资源| 99久久精品热视频| 大又大粗又爽又黄少妇毛片口| 男人舔奶头视频| 成人鲁丝片一二三区免费| 国产伦精品一区二区三区视频9| 国内精品久久久久精免费| 听说在线观看完整版免费高清| 日韩精品中文字幕看吧| 在现免费观看毛片| 亚洲乱码一区二区免费版| 色尼玛亚洲综合影院| 九色成人免费人妻av| 国产伦一二天堂av在线观看| 最后的刺客免费高清国语| 精品人妻偷拍中文字幕| 女的被弄到高潮叫床怎么办| 国产精品一区二区三区四区久久| 伦精品一区二区三区| 国模一区二区三区四区视频| 又粗又爽又猛毛片免费看| 欧美极品一区二区三区四区| 波多野结衣高清作品| 国产精品不卡视频一区二区| 一区二区三区免费毛片| 日韩av不卡免费在线播放| 成熟少妇高潮喷水视频| 免费高清视频大片| 亚洲欧美日韩东京热| 久久国内精品自在自线图片| 国产男靠女视频免费网站| 成年av动漫网址| 国产精品一区二区三区四区久久| 欧美人与善性xxx| 18禁裸乳无遮挡免费网站照片| 嫩草影视91久久| 中文字幕免费在线视频6| 国产亚洲91精品色在线| 内射极品少妇av片p| 午夜福利视频1000在线观看| 国产精品福利在线免费观看| 国产又黄又爽又无遮挡在线| 成人漫画全彩无遮挡| 黄色日韩在线| 国产精品美女特级片免费视频播放器| www日本黄色视频网| 久久久精品大字幕| 最近中文字幕高清免费大全6| 精品熟女少妇av免费看| 又爽又黄无遮挡网站| 日日撸夜夜添| 少妇的逼水好多| 天天躁日日操中文字幕| 久久精品人妻少妇| 日本在线视频免费播放| a级毛片a级免费在线| 国产精品1区2区在线观看.| 一进一出抽搐动态| 久久婷婷人人爽人人干人人爱| 真人做人爱边吃奶动态| 亚洲高清免费不卡视频| 欧美又色又爽又黄视频| 久久久久精品国产欧美久久久| 淫秽高清视频在线观看| 国产精品1区2区在线观看.| 欧美xxxx性猛交bbbb| 婷婷六月久久综合丁香| 18+在线观看网站| 国模一区二区三区四区视频| av中文乱码字幕在线| 中文字幕人妻熟人妻熟丝袜美| 嫩草影院精品99| 国产亚洲精品久久久久久毛片| 99热这里只有是精品在线观看| 我的女老师完整版在线观看| 国产成人一区二区在线| 男人的好看免费观看在线视频| 午夜福利18| 亚洲在线观看片| 亚洲七黄色美女视频| 日本黄大片高清| 女人被狂操c到高潮| 国产精品一区二区三区四区久久| 久久久久国产网址| 亚洲国产精品成人久久小说 | 久久天躁狠狠躁夜夜2o2o| 日韩精品青青久久久久久| 尤物成人国产欧美一区二区三区| 国产精品爽爽va在线观看网站| 欧美区成人在线视频| 美女高潮的动态| 亚洲国产欧洲综合997久久,| 国产毛片a区久久久久| 一进一出抽搐动态| 香蕉av资源在线| 国产伦一二天堂av在线观看| 久久久色成人| 日韩精品有码人妻一区| 国产精品人妻久久久久久| 午夜免费激情av| 欧美又色又爽又黄视频| 精品午夜福利视频在线观看一区| 日产精品乱码卡一卡2卡三| 在线免费观看的www视频| 十八禁国产超污无遮挡网站| 亚洲av成人av| 久久久精品大字幕| 国产一区二区在线av高清观看| 中文字幕av在线有码专区| 欧美日韩乱码在线| 国产高清不卡午夜福利| 夜夜夜夜夜久久久久| 十八禁网站免费在线| 网址你懂的国产日韩在线| 国产中年淑女户外野战色| 99热只有精品国产| 日本熟妇午夜| 美女 人体艺术 gogo| 丰满人妻一区二区三区视频av| 亚洲最大成人中文| 午夜激情福利司机影院| 亚洲av美国av| 亚洲av五月六月丁香网| 麻豆精品久久久久久蜜桃| 波多野结衣高清作品| 一个人看的www免费观看视频| 秋霞在线观看毛片| 国产高清视频在线观看网站| 国产白丝娇喘喷水9色精品| 亚洲成a人片在线一区二区| 国产亚洲91精品色在线| 黄色日韩在线| 又爽又黄无遮挡网站| 亚洲国产欧美人成| 国产精品久久久久久久电影| 欧美极品一区二区三区四区| 狂野欧美激情性xxxx在线观看| 亚洲第一区二区三区不卡| 老司机福利观看| 亚洲性夜色夜夜综合| 两个人的视频大全免费| 国产精品永久免费网站| 搡老妇女老女人老熟妇| 波多野结衣高清无吗| 国产色婷婷99| 一边摸一边抽搐一进一小说|