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

    渡槽水體的熱交換及水溫和含冰率的時空變化

    2022-12-05 06:28:56楊開林
    水利學(xué)報 2022年11期
    關(guān)鍵詞:渡槽外壁太陽輻射

    楊開林

    (流域水循環(huán)模擬與調(diào)控國家重點實驗室,中國水利水電科學(xué)研究院,北京 100038)

    1 研究背景

    渡槽是一種輸水的混凝土水工建筑物,槽身多為板梁結(jié)構(gòu),過水?dāng)嗝嬉话悴捎镁匦魏蚒型。渡槽又稱高架渠,通常架設(shè)于山谷、洼地、河流之上。與河渠不同,渡槽外表面裸露在大氣環(huán)境中,除水面與大氣熱交換外,水體與內(nèi)壁的熱交換也受外壁與大氣熱交換的影響。在一般情況下,渡槽過水?dāng)嗝姹冗B接渠道小得多,所以流速較大,冰期常常不封凍。這些導(dǎo)致寒潮時單位長度渡槽水溫比相同過水能力渠道下降更快,而暖潮時單位長度渡槽水溫比相同過水能力渠道上升更快。一旦受到嚴(yán)重寒潮侵襲,渡槽水溫可能很快降低到0 ℃以下,成為產(chǎn)冰場,水體中產(chǎn)生大量冰花,導(dǎo)致渡槽下游渠道發(fā)生嚴(yán)重冰塞災(zāi)害,特別是當(dāng)渡槽長度較大時。以南水北調(diào)中線輸水工程為例,沿途修建了21座大跨度渡槽,最長渡槽長11 900 m,2016年1月21日受到極端寒潮侵襲,最低氣溫達-18.6 ℃,低溫天氣持續(xù)了3天,漕河渡槽(長2300 m)水溫下降到0 ℃以下,產(chǎn)生了大量冰花,結(jié)果導(dǎo)致下游崗頭隧洞前發(fā)生了嚴(yán)重冰塞險情[1-2]。因此,當(dāng)分析長距離輸水工程的冰情時,需要考慮渡槽水體與環(huán)境的熱交換。

    渡槽水體與環(huán)境的熱交換包括水面與大氣的熱交換及水體與槽體的熱交換。在水面與大氣的熱交換計算方面,目前一般采用國際著名冰工程專家Ashton[3]和沈洪道[4]的傳統(tǒng)線性化模型,該模型只考慮氣溫和水溫的影響,并在工程計算中得以普遍應(yīng)用[5-6]。在水體與床底的熱交換方面,Jobson[7]和茅澤育等[8]采用數(shù)值求解非穩(wěn)態(tài)熱傳導(dǎo)偏微分方程的方法,研究了河床對水溫的影響。河渠與大氣的熱交換非常復(fù)雜,僅僅考慮氣溫和水溫的影響,常常不可避免會導(dǎo)致計算結(jié)果遠遠偏離實際情況。以南水北調(diào)中線工程為例,自2014年全線正式通水以來,在中線工程中進行了連續(xù)多年的冰情原型觀測[1-2,9],除2015—2016年度極端寒潮期冰情嚴(yán)重外,其它冬季以暖冬為主,主要冰情為岸冰、流冰,未形成持續(xù)冰蓋,但是采用傳統(tǒng)線性化熱交換模型的計算結(jié)果卻是形成連續(xù)冰蓋及冰蓋較厚[10-11],與實際情況不符。鑒于此,楊開林[12-15]研究了冰期水面與大氣及水體與床底的熱交換模型,考慮了太陽輻射、地面與大氣的長波輻射、水面蒸發(fā)和對流的熱交換,以及熱交換模型的線性化與水冰雪反照率的參數(shù)化,結(jié)果表明:(1)水面與大氣的傳統(tǒng)線性熱交換模型僅適用于氣候比較嚴(yán)寒地區(qū)封河期,例如我國東北,但不適用于我國華北地區(qū),包括南水北調(diào)中線,因為這些地區(qū)冬季日平均氣溫常常大于-10 ℃;(2)水體與床底的熱交換必須考慮。

    雖然迄今冰期河渠水溫模型的研究正在逐步完善,但是將其用來計算冰期渡槽水溫不可避免的會與實際產(chǎn)生很大的偏差。河渠水溫受床底地溫的影響,而地溫一般大于0 ℃。但是,由于渡槽架設(shè)空中,水溫受渡槽外壁溫度的影響,而外壁溫度取決于外壁與大氣的熱交換,包括大氣對流、外壁和大氣的長波輻射、地表的長波輻射及太陽輻射,包括地表太陽輻射的反射等。換句話說,冰期氣溫常常小于0 ℃,所以渡槽外壁溫度也常常小于0 ℃。存在的問題是目前在冰水力學(xué)中缺乏渡槽水體與槽體熱交換模型方面的研究。

    本文的主要目的是,基于太陽輻射理論、建筑物熱過程和冰水力學(xué)基礎(chǔ)理論,首先建立水體與渡槽的熱交換參數(shù)化模型,然后研究冰期渡槽的水溫模型,并分析水溫和含冰率的時空變化規(guī)律。

    2 水體與渡槽的熱交換

    下面將以圖1所示矩形渡槽為例研究水體與渡槽的熱交換,圖中:B為水面寬,H為水深,hE,1和hE,2為邊墻厚,hE,3為底板厚。

    圖1 矩形渡槽斷面示意圖

    在冰期正常輸水工況下,渡槽內(nèi)水溫一般分布均勻[9]。渡槽一般采用鋼筋混凝土結(jié)構(gòu),一般由多層材料構(gòu)成,當(dāng)假設(shè)每層厚度均勻、熱力條件是準(zhǔn)穩(wěn)態(tài)的且熱量傳遞僅在平板豎直方向進行時,根據(jù)楊開林[14]的研究,水體與渡槽的熱交換可通過多層固體平板的熱傳導(dǎo)來計算

    (1)

    (2)

    式中:φwa為渡槽外壁的凈熱通量,W/m2;kE為邊墻或底板的等效導(dǎo)熱系數(shù),W/(m·℃);hE為邊墻或底板的等效厚度,m;Ts為邊墻或底板外壁的溫度,℃;Tw為水溫,℃;hj為邊墻或底板第j層的厚度,m;kj為邊墻或底板第j層的導(dǎo)熱系數(shù),W/(m·℃);m為邊墻或底板斷面結(jié)構(gòu)的層數(shù);下標(biāo)i=1,2,3分別表示渡槽兩邊墻和底板。

    在一般情況下,兩邊墻和底板外壁的溫度存在差別,其大小取決于外壁與大氣的熱交換,包括大氣的對流、長波輻射、太陽輻射及地表的長波輻射和地表太陽輻射的反射等。

    2.1 渡槽外壁與大氣的對流熱交換當(dāng)渡槽邊墻和底板是混凝土平板時,外壁與大氣的對流熱交換可描述為

    φc,i=hc,i(Ts,i-Ta),i=1,2,3

    (3)

    式中:φc為外壁的對流熱通量,W/m2;hc為渡槽外壁對流熱交換系數(shù),W/(m2·℃);Ta為氣溫,℃。

    渡槽外壁對流熱交換系數(shù)hc可利用凱爾別克的經(jīng)驗公式計算[16]

    hc,i=3.67+3.83Va,i=1,2

    (4a)

    hc,i=2.17+3.83Va,i=3

    (4b)

    式中Va為風(fēng)速,m/s。

    2.2 渡槽外壁的長波輻射熱交換渡槽外壁上的長波輻射包括外壁的長波輻射、大氣長波逆輻射和地表的長波輻射[17]

    φabθ,i=φa,i(φa-φb)+φg,i(φg-φb),i=1,2,3

    (5)

    式中:φabθ為外壁上的長波輻射凈熱通量,W/m2;φa為大氣的長波逆輻射熱通量,W/m2;φb外壁的長波輻射熱通量,W/m2;φg為地面的長波輻射熱通量,W/m2;φa,i和φg,i分別為外壁對天空和對地面的輻射角系數(shù)。如果將渡槽外壁、地面和大氣視為一個封閉系統(tǒng),φa,i+φg,i=1.0。對于矩形渡槽邊墻外壁,φa,i=φg,i=0.5,i=1,2;對于底板外壁,φa,3=0和φg,3=1.0。

    外壁長波輻射φb和大氣長波逆輻射φa可采用下述線性方程描述[13]

    φab=φa-φb=φab0-hba(Ts,i-Ta),i=1,2

    (6)

    式中:hba為外壁長波輻射的熱交換系數(shù),W/(m2·℃);φab0為氣溫Ta的線性函數(shù),W/m2,例如對河北保定地區(qū)[13]

    φab0=-94.5-0.6Ta,hba=3.9

    (7)

    外壁長波輻射φb和地表長波輻射φg可由斯蒂弗-博爾茨曼四次方定律計算[3]

    φb=σεwb(273.15+Ts,i)4,i=1,2,3

    (8a)

    φg=σεwg(273.15+Tg)4

    (8b)

    式中:σ=5.67×10-8W/(m2·K4)為斯特凡·玻爾茲曼常數(shù);εwg為地表的黑度;Tg為地表溫度,℃;εwb為底板外壁的黑度,對混凝土可取εwb=0.94。

    當(dāng)假設(shè)εwg=εwb,根據(jù)式(8a)計算可得如圖2所示φb隨Ts變化的關(guān)系曲線(點線),采用線性回歸的方法得

    圖2 長波輻射φb的線性回歸

    φb=305.14+3.89Ts,i,i=1,2,3

    (9)

    φg=305.14+3.89Tg

    (10)

    根據(jù)式(9)(10)可得

    φg-φb=-hgb(Ts,i-Tg),i=1,2,3

    (11)

    式中hgb≈3.9。

    由于hgb=hba,把式(6)(11)代入式(5)得渡槽外壁上的長波輻射凈熱通量

    φabθ,i=φa,iφab0-φa,ihba(Ts,i-Ta)-φg,ihba(Ts,i-Tg),i=1,2,3

    (12)

    2.3 渡槽外壁與太陽輻射的熱交換渡槽外壁的太陽輻射凈熱通量等于太陽直射、散射和地表反射凈熱通量之和,即

    φcn=φDn+φdn+φRn

    (13)

    式中:φcn為外壁太陽輻射凈熱通量,W/m2;φDn為外壁太陽直射凈熱通量,W/m2;φdn為外壁太陽散射凈熱通量,W/m2;φRn為外壁接受的地表太陽輻射的反射凈熱通量,W/m2。

    如圖3所示,傾斜面上的晴天太陽直射[17-19]

    圖3 傾斜面的太陽輻射示意圖

    φDθ=I0E0Pmcos(λ)

    (14)

    式中:φDθ為傾斜面上晴天太陽直射熱通量,W/m2;I0=1367 W/m2為太陽常數(shù);E0為日地距離訂正系數(shù);P為所有波長范圍內(nèi)的大氣平均透明系數(shù);m為光學(xué)大氣質(zhì)量;λ為太陽入射角,是太陽入射線與傾斜面法線之間的夾角。

    太陽入射角λ與太陽高度角、傾斜面的方位角和太陽方位角的關(guān)系是

    cos(λ)=cos(θ)sin(α)+sin(θ)cos(α)cos(β-γ)

    (15)

    式中:α為太陽高度角,rad;γ為傾斜面的方位角,rad,是傾斜面的法線在水平面上的投影與正南方向線的夾角,規(guī)定:面向東時為負,面向西時為正;β為太陽方位角,rad,太陽至地面上某給定點的連線在水平面上的投影與正南向(當(dāng)?shù)刈游缇€)的夾角。

    對于渡槽邊墻,θ=π/2,由式(15)得

    cos(λ)=cos(α)cos(β-γ)

    (16a)

    對于渡槽水面,θ≈0,由式(15)得

    cos(λ)=sin(α)

    (16b)

    把式(16b)代入式(14)可得水平地表的晴天太陽直射計算公式。

    根據(jù)式(14)(16),考慮天氣對太陽輻射的影響及邊墻對太陽直射的反射,可得邊墻外側(cè)云天太陽直射凈熱通量[13]

    φDn=(1-ac)(1-0.65C2)I0E0Pmcos(α)cos(β-γ)

    (17)

    式中:ac為邊墻反照率,水泥砂漿抹面取0.32;C為云量,取值0~1,對晴天取0,而陰天取1。當(dāng)邊墻外側(cè)面向正東時,則γ=-π/2而cos(β-γ)=-sinβ;當(dāng)邊墻外側(cè)面向正西時,則γ=π/2而cos(β-γ)=sinβ。

    太陽方位角β與觀測地的地理緯度L、赤緯角δ和時角ω的關(guān)系是

    (18)

    當(dāng)sin(β)>1,則

    (19)

    這時,β的正負由時角ω確定。赤緯角δ和時角ω的計算可參考文獻[13]。

    傾斜面上晴天太陽散射熱通量[18-19]是

    φdθ=φsdcos2(θ/2)

    (20)

    式中:φdθ為邊墻晴天太陽散射熱通量,W/m2;φsd為水平面上晴天太陽散射熱通量,W/m2。

    對于渡槽邊墻,θ=π/2,考慮天氣對太陽輻射的影響及邊墻對太陽直射的反射[13],可得云天太陽散射時邊墻外側(cè)凈熱通量

    φdn=0.5φsd(1-ac)(1-0.65C2)

    (21)

    晴天傾斜面獲得地面反射熱通量是[17]

    φRθ=aeφssin2(θ/2)

    (22)

    式中:φRθ為邊墻接受的晴天地表反射凈熱通量,W/m2;ae為地表太陽輻射的反照率;φs為水平地表的晴天日射熱通量,W/m2。

    考慮云層對太陽輻射的影響及傾斜面對接受的地表太陽輻射的反射作用[13],可得傾斜面接受的地表太陽輻射反射的凈熱通量

    φRn=φRθ(1-ac)(1-0.65C2)

    φRn=aeφs(1-ac)(1-0.65C2)sin2(θ/2)

    (23)

    對于矩形渡槽垂直邊墻,θ=π/2,把式(17)(21)(23)代入式(13)得矩形渡槽邊墻外壁太陽輻射凈熱通量

    φcs,i=(1-ac)(1-0.65C2)(I0E0Pmcos(α)cos(β-γ)+0.5φsd+0.5aeφs)

    (24)

    式中下標(biāo)i=1和2分別表示渡槽兩邊墻。在一般情況下,渡槽兩側(cè)邊墻太陽輻射熱通量存在差異。

    對于渡槽底板,θ=π,影響外壁的太陽輻射只有地表反射熱通量

    φcs,i=ae(1-ac)φs(1-0.65C2)

    (25)

    式中下標(biāo)i=3。草地的反照率ae=0.174~0.219,混凝土面的反照率ac=0.33~0.37。

    式(24)(25)中水平地表的晴天日射和散射的熱通量可采用下述公式計算[13]

    φs=0.5(1+Pm)I0E0sin(α)

    (26)

    φsd=0.5(1-Pm)I0E0sin(α)

    (27)

    式(24)—(27)中太陽高度角α、日地距離訂正系數(shù)E0、Pm的計算可參考文獻[13],因此對于給定位置的渡槽和云量C,φcs,i可計算得到。

    下面以保定地區(qū)為例,了解水面與渡槽外壁太陽輻射的差異。假設(shè)有一矩形渡槽兩邊墻外壁分別面朝正東和正西,即上午只有面朝正東邊墻外壁接受太陽直射,而下午只有面朝正西邊墻外壁接受太陽直射。在此條件下,在上午,面朝正東邊墻外壁太陽輻射包括太陽直射、散射和地面反射,而面朝正西邊墻外壁太陽輻射只有太陽散射和地面反射;在下午,面朝正東邊墻外壁太陽輻射只有太陽散射和地面反射,而面朝正西邊墻外壁太陽輻射包括太陽直射、散射和地面反射。當(dāng)太陽輻射時間與天文輻射相同、渡槽混凝土面的反照率ac=0.33、地表反照率ae=0.2、云量C取歷史天氣資料,則可得如圖4所示保定地區(qū)2015—2016年冬季水面與渡槽外壁太陽輻射凈熱通量隨時間的變化。圖中:φsn為水面太陽輻射日平均凈熱通量,計算模型見下一節(jié);φcs,1為上午渡槽面朝正東邊墻外壁太陽輻射平均凈熱通量(等于下午渡槽面朝正西邊墻外壁太陽輻射平均凈熱通量);φcs,3為底板外壁接受的地表太陽輻射反射的日平均凈熱通量。

    從圖4可見,水面太陽輻射、邊墻外側(cè)太陽輻射和底板外側(cè)太陽輻射的關(guān)系是φsn>φcs,1>φcs,3,晴天時φcs,1=170~280 W/m2,陰天時φcs,1有時小于100 W/m2;晴天時φsn=230~400 W/m2,陰天時φsn最小接近100 W/m2;φcs,3一般在50 W/m2左右變化。

    圖4 保定地區(qū)2015年12月—2016年2月冬季水面與渡槽外側(cè)太陽輻射凈熱通量

    2.4 渡槽外壁溫度與氣溫的關(guān)系渡槽外壁與大氣的熱交換包括太陽輻射、對流和長波輻射,

    φwa,i=φcs,i-φc,i+φabθ,i,i=1,2,3

    (28)

    式中φwa為外壁的凈熱通量,W/m2。

    把式(3)(12)代入式(28)得

    φwa,i=φcs,i+φa,iφab0-(hc,i+φa,ihba)(Ts,i-Ta)-φg,ihba(Ts,i-Tg),i=1,2,3

    (29)

    聯(lián)立求解式(1)(29)得

    (Ts,i-Tw)kE,i/hE,i=φcs,i+φa,iφab0-(hc,i+φa,ihba)(Ts,i-Ta)-φg,ihba(Ts,i-Tg),i=1,2,3

    因為φa,i+φg,i=1,整理得渡槽外壁溫度

    (30)

    上式適用于渡槽架于各種地表上,包括土石地面、水面、冰面和雪面等。

    當(dāng)渡槽架于土石地面上,包括地面積冰和積雪時,為使問題簡化,一般假設(shè)地表溫度接近于氣溫[17],即Tg≈Ta,這時

    (31a)

    式中

    (31b)

    2.5 渡槽外壁的凈熱通量與氣溫和水溫的關(guān)系把式(31)代入式(1)得渡槽外壁凈熱通量

    整理得

    φwa,i=hwa,i(φcs,i+φa,iφab0)/(hc,i+hba)+hwa,i(Ta-Tw),i=1,2,3

    (32)

    hwa,i=1/[hE,i/kE,i+1/(hc,i+hba)]

    (33)

    式中hwa,i為渡槽外壁熱交換系數(shù),W/(m2·℃)。

    當(dāng)渡槽邊墻和底板沿流動方向是均值平板時,則渡槽外壁的等效凈熱通量

    (34)

    (35)

    式中:φwae為渡槽外壁的等效凈熱通量,W/m2;hwae為渡槽外壁的等效熱交換系數(shù),W/(m2·℃);χ=χ1+χ2+χ3為渡槽濕周,m;χ1和χ2分別為兩邊墻的濕周,m;χ3為底板的濕周,m。

    3 渡槽水溫的時空變化規(guī)律

    隨著渡槽水體與周圍環(huán)境的熱交換,包括水面與大氣的熱交換與水體與渡槽邊墻和底板的熱交換,熱量通過水流的運動和紊動向下游和整個過水?dāng)嗝鎮(zhèn)鬟f,在一維條件下,沿流向的對流-熱擴散方程是

    (36)

    式中:t為時間,s;x為距離,m;ρ為水的密度,kg/m3,ρ≈1000 kg/m3;Cp為水的比熱,在0 ℃時Cp=4217.7 J/(kg·℃);A為渡槽過水?dāng)嗝婷娣e,m2;V為斷面平均流速,m/s;Ex為熱擴散系數(shù);φsa為水面凈熱通量,W/m2;左邊第一項表示過水?dāng)嗝鏌崃侩S時間的變化;第二項表示過水?dāng)嗝鏌崃侩S水體運動的變化,又稱為熱量的對流傳遞;第三項表示過水?dāng)嗝鏌崃侩S水體熱擴散的變化。

    當(dāng)水面凈熱通量φsa采用參考文獻[14]線性化熱交換模型時,則水面與大氣的熱交換與水體與渡槽邊墻和底板的熱交換可描述為

    Bφsa+χφwbe=B(φsa0+hsa(Ta-Tw)-hsa2(Tw-Ta)2)+χ(φwae0+hwae(Ta-Tw))

    (37)

    φsa0=φsn-asa-bsaTa

    (38)

    式中:φsa為水面凈熱通量,W/m2;hsa為水面與大氣的熱交換系數(shù),W/(m2·℃);hsa2=0.158×10-3pa,pa為當(dāng)?shù)卮髿鈮?,hPa;φsn為水面太陽輻射凈熱通量,W/m2;系數(shù)asa?0和bsa>0,以保定地區(qū)為例,

    asa=94.6+(6.04+2.95Vz)(1-Rh)c1≥94.6,bsa=0.6+(6.04+2.95Vz)(1-Rh)c2

    (39)

    式中:Vz為水面上1.5m的風(fēng)速,m/s;Rh為大氣的相對濕度;

    在我國的北京、保定、沈陽、包頭等地區(qū),熱交換系數(shù)可采用下式計算[13]

    hsa=10.0(1.0+0.25Vz)

    (40)

    參考文獻[12-13],有

    φsn=0.5I0E0sinα(1+Pm)(1-0.65C2)(1-as)

    (41)

    式中as為水面反照率

    as=0.0564/α,α≥0.105

    (42a)

    as=0.537-4.408(α-0.105),0≤α<0.105

    (42b)

    式(37)可改寫為

    (43)

    式中:

    (44)

    在冰水力學(xué)分析中,水體熱擴散項可忽略不計[3-4]。把式(43)代入式(36),并采用特征線方法,式(36)偏微分方程可改寫為

    (45)

    dx/dt=V

    (46)

    (47)

    式中

    (48)

    式中:下標(biāo)“i”為斷面編號;m為斷面數(shù);Twp,i-1為渠段i進口時刻ti-1的水溫, ℃;Twp,i為渠段i出口時刻ti的水溫, ℃;Δti=ti-ti-1,s。由于系數(shù)r3,i>0總是成立,可得重要結(jié)論:在給定的氣象條件下,水溫Twp,i隨著時間ti=t0+Δt1+Δt2+…+Δti的增加呈指數(shù)規(guī)律變化,逐漸趨近于系數(shù)-r1,i。

    對式(46)積分并取一階近似可得

    Δti=Δxi/Vi

    (49)

    式中Δxi=xi-xi-1為渠段i的長度,m。

    把式(49)代入式(47)消去Δti得水溫的沿程分布遞推計算公式,

    (50)

    據(jù)此,可得重要結(jié)論:在給定的氣象條件下,水溫Twp,i隨離開進口距離xi=x0+Δx1+Δx2+…+Δxi的增加呈指數(shù)規(guī)律變化,逐漸趨近于系數(shù)-r1,i。

    對于均勻流,在冬季寒潮情況下,如果參數(shù)a、b、c為常數(shù),當(dāng)令Twp=0,則由式(50)可得水溫由渠道進口x0=0的溫度Tw0下降到0 ℃的條件

    (51)

    4 渡槽含冰率的時空變化規(guī)律

    當(dāng)xp>xpc時,渡槽中將開始產(chǎn)生冰花。在一維流動條件下,如果水面沒有漂浮冰塊,則沿流向的渡槽冰花對流-擴散方程可描述為

    (52)

    式中:Ci=Ai/A為含冰率,表示單位長度渡槽中冰花體積與水和冰花總體積之比,其中Ai為過水?dāng)嗝鍭中冰花所占面積;ρi為冰的密度,kg/m3,一般取ρi≈917 kg/m3;Li為冰的潛熱,在0 ℃時Li=3.33×105J/kg;Ei為冰花擴散系數(shù)。式(52)左邊第一項表示單位長度渡槽中水和冰相變的潛熱隨時間的變化;第二項表示單位長度渡槽中水和冰相變的潛熱隨水體運動的變化,又稱為熱量的對流傳遞;第三項表示單位長度渡槽中水和冰相變的潛熱隨水體熱擴散的變化;

    在冰水力學(xué)分析中,熱擴散項可忽略不計[3-4]。把式(43)代入式(52),并采用特征線方法,式(52)偏微分方程可改寫為

    (53)

    式中系數(shù)a、b、c與式(45)相同。

    根據(jù)式(45)(53)可得

    積分得

    (54)

    式中:下標(biāo)“j”為從xpc開始的斷面編號;ρCP/(ρiLi)≈0.0138。

    根據(jù)式(54)可得重要結(jié)論:渡槽斷面含冰率與負水溫成正比沿程增加,負水溫越小,水相變?yōu)楸牧吭酱螅划?dāng)水溫由負轉(zhuǎn)正時,水中冰將相變?yōu)樗?,含冰率將沿程減小。

    5 算例

    漕河渡槽是南水北調(diào)中線總干渠上的一座大型交叉建筑物,位于河北省滿城縣西北約9 km處,上游通過2000 m的總干渠與吳山隧洞銜接,下游通過750 m矩形干渠與崗頭隧洞銜接。漕河渡槽全長L=2300 m,采用3槽一聯(lián)整體多縱墻結(jié)構(gòu),簡支梁式,單孔斷面凈寬6.0 m,槽高5.4 m。槽身底板厚0.5 m,邊墻厚0.6 m。渡槽進、出口底高程為62.243~61.525 m。在2015—2016年冬季輸水流量Q=45.72 m3/s,水深H=3.76 m,平均流速0.67 m/s。在2016年1月21—23日期間發(fā)生罕見寒潮,48小時內(nèi)氣溫下降近10 ℃,實測最低氣溫-18.6 ℃,2016年1月23日,崗頭隧洞進口出現(xiàn)嚴(yán)重冰塞,進口前兩道攔冰索相繼斷裂,冰塞體整體下移,強烈撞擊隧洞進口導(dǎo)墻,威脅隧洞節(jié)制閘安全。事后分析,主要由于上游長距離明渠內(nèi)的大量流冰全部堆積在隧洞進口,形成較大體積的冰塞,順流向長度約700 m,典型剖面結(jié)構(gòu)可分為3層:上層為冰堆,厚20~40 cm,碎冰塊堆積;中層為冰蓋,厚25 cm,堅硬密實;下層為冰花,厚70~230 cm,絮狀松散,具有一定阻水作用,渠道上游水位壅高0.4 m。崗頭隧洞發(fā)生嚴(yán)重冰塞表明漕河渡槽中產(chǎn)生了大量冰花,水溫Tw≤0 ℃。下面主要計算分析:1)風(fēng)速和太陽輻射對渡槽外壁溫度的影響;2)水體與渡槽熱交換對水溫和冰花沿程分布的影響。

    原型觀測漕河渡槽3槽內(nèi)的水溫分布均勻[9],中墻對水溫分布的影響可忽略不計,只需考慮水面和渡槽邊墻和底板的熱交換。

    計算取大氣相對濕度Rh=0.9,夜間Ta=-18.6 ℃,白晝Ta=-9.0 ℃和云量C=0.14,并假設(shè)渡槽進口水溫Tw0=0.1 ℃。因水深3.76 m,邊墻濕周χ1=χ2=3.76 m,3槽底板凈寬χ3=3×6.0=18.0 m,渡槽斷面濕周χ=χ1+χ2+χ3=25.52 m。邊墻厚h1,i=0.6 m,底板厚h1,3=0.5 m,鋼筋混凝土導(dǎo)熱系數(shù)k1,i=1.74 W/(m·℃),由式(2)得

    保定地區(qū)緯度38°10′—40°00′,2016年1月23日上午,計算的渡槽面朝正東邊墻外壁平均太陽輻射φcs,1=194.9 W/m2,面朝正西邊墻外壁平均太陽輻射φcs,2=55.8 W/m2;下午,渡槽面朝正西邊墻外壁平均太陽輻射φcs,1=194.9 W/m2,面朝正東邊墻外壁平均太陽輻射φcs,2=55.8 W/m2;底板外壁的地表平均太陽輻射反射φcs,3≈44.5 W/m2。水面平均太陽輻射φsn=282.9 W/m2。

    5.1 風(fēng)速和太陽輻射對渡槽外壁溫度的影響在上述計算條件下,可得如圖5和圖6所示渡槽外壁溫度Ts與風(fēng)速Va的關(guān)系曲線,圖中:曲線1、2分別為兩邊墻外壁的Ts-Va,曲線3為底板外壁的Ts-Va。在夜間沒有太陽輻射,兩邊墻外壁溫度相同。在白晝,由于兩邊墻太陽輻射的差別,兩邊墻外壁溫度差別較大,如圖6所示。

    圖5 夜間渡槽外壁溫度隨風(fēng)速的變化

    圖6 白晝渡槽外壁溫度隨風(fēng)速的變化

    觀察圖5和圖6,可得下述結(jié)論:

    1)渡槽外壁溫度隨風(fēng)速的增加而下降,特別是在0≤Va≤2 m/s區(qū)間下降較快;

    2)在夜間,底板外壁溫度高于邊墻外壁溫度較大,而邊墻外壁溫度接近于氣溫;

    3)在白晝太陽輻射的影響下,渡槽兩邊墻外壁溫度相差較大,而底板外壁溫度介于兩邊墻外壁溫度之間,同時朝陽面邊墻外壁溫度遠遠高于氣溫,特別是風(fēng)速較小時,外壁溫度隨著風(fēng)速的下降由負轉(zhuǎn)正。

    5.2 水體與渡槽熱交換對水溫和冰花沿程分布的影響當(dāng)取風(fēng)速Vz=Va=3.0 m/s,計算的渡槽夜間寒潮氣溫Ta=-18.6 ℃時水溫沿特征線的變化如圖7所示,顯然,水溫Twp隨離開渡槽進口距離x的增加而近似線性下降,其中:曲線1考慮了水體與渡槽的熱交換,水溫由進口水溫Tw0=0.1 ℃下降到0 ℃計算的臨界長度xpc1=1954 m;曲線2沒有考慮水體與渡槽的熱交換,即取φwae0=0和hwae=0,對應(yīng)xpc2=2273 m,兩者相差較大,偏差相對值|xpc2-xpc1|/xpc1=16%。由于xpc1和xpc2小于渡槽長度L=2300 m,所以渡槽L-xpc1或者L-xpc2部分會產(chǎn)生冰花。在渡槽出口,對應(yīng)xpc1的溫度Twp1=-0.0012 ℃;對應(yīng)xpc2的溫度Twp2=-0.0018 ℃。根據(jù)式(54),可得渡槽出口含冰率的偏差相對值|Cip2-Cip1|/Cip1=|Twp2-Twp1|/|Twp1|=50%。這一結(jié)果表明,在給定氣象條件下,如果不考慮水體與渡槽熱交換的影響,則計算的渡槽出口含冰率比實際含冰率少50%。從中可得一個重要結(jié)論:考慮水體與渡槽的熱交換比不考慮的冰情嚴(yán)重得多,兩者相差較大。

    圖7 夜間寒潮時渡槽水溫沿特征線的變化

    計算的白晝氣溫Ta=-9 ℃時渡槽水溫沿特征線的變化如圖8所示。由于白晝太陽輻射的影響,當(dāng)渡槽進口水溫Tw0=0.1 ℃時,水溫沒有下降到0 ℃,渡槽內(nèi)不會產(chǎn)生冰花。此外,從圖8可見,考慮水體與渡槽的熱交換時渡槽水溫沿特征線緩慢降低,但是不考慮水體與渡槽的熱交換時渡槽水溫沿特征線反而略微增加,兩者偏差明顯。

    圖8 白晝寒潮時渡槽水溫沿特征線的變化

    綜上所述,可得一個重要結(jié)論:在寒潮期間,在渡槽中可能發(fā)生夜間產(chǎn)生冰花而白天不產(chǎn)生冰花的現(xiàn)象。需要說明的是,本算例水溫和含冰率的計算結(jié)果是在給定渡槽進口水溫的條件下得到的,不能完全反映真實的冰塞過程,只是說明考慮水體與渡槽邊墻和底板熱交換對正確模擬冰情的重要性。

    6 結(jié)論

    渡槽架設(shè)空中,水溫不僅受水面與大氣熱交換的影響,而且受水體與渡槽熱交換的影響。影響水體與渡槽熱交換的主要因素是渡槽外壁溫度。考慮渡槽外壁與大氣的熱交換,包括對流、外壁和大氣的長波輻射、地表的長波輻射、太陽輻射及地表的反射等,本文提出了渡槽外壁溫度的計算公式、水體與渡槽的熱交換參數(shù)化模型、水溫和含冰率的計算模型等。

    通過研究,得出下述主要結(jié)論:

    (1)渡槽外壁溫度隨風(fēng)速的增加而下降,特別是在0≤Vz≤2 m/s區(qū)間;在夜間,底板外壁溫度高于邊墻外壁溫度,且邊墻外壁溫度接近于氣溫;在白晝太陽輻射的影響下,渡槽兩邊墻外壁溫度相差較大,而底板外壁溫度介于兩邊墻外壁溫度之間,同時朝陽面邊墻外壁溫度遠遠高于氣溫,特別是風(fēng)速較小時,外壁溫度隨著風(fēng)速的下降由負轉(zhuǎn)正;

    (2)在給定的氣象條件下,水溫隨著時間和離開渡槽進口距離的增加呈指數(shù)規(guī)律變化;

    (3)渡槽斷面含冰率與負水溫成反比沿程增加,負水溫越小,水相變?yōu)楸牧吭酱螅划?dāng)水溫由負轉(zhuǎn)正時,水中冰將相變?yōu)樗蕦⒀爻虦p?。?/p>

    (4)考慮水體與渡槽的熱交換比不考慮的冰情嚴(yán)重得多,并且在寒潮期間,在渡槽中可能發(fā)生夜間產(chǎn)生冰花而白天不產(chǎn)生冰花的現(xiàn)象。

    猜你喜歡
    渡槽外壁太陽輻射
    淺談水利工程渡槽基礎(chǔ)承臺施工
    邯鄲太陽輻射時空分布特征
    傍山道路泥石流渡槽設(shè)計
    壁掛式鋼板立表面除銹機 在利舊鋼制儲罐外壁除銹的應(yīng)用
    袁曉峰攝影作品欣賞
    河北水利(2017年6期)2017-09-05 11:04:13
    基于PCA 的太陽輻射觀測算法研究
    燒水時燒水壺外壁為何會“出汗”
    太陽輻射作用下鋼筒倉結(jié)構(gòu)溫度場分析研究
    基于ANSYS的渡槽三維有限元分析
    非能動核電站安全殼外壁下降水膜的穩(wěn)定性分析
    国产xxxxx性猛交| 欧美日本中文国产一区发布| 中文字幕人妻丝袜制服| 亚洲伊人久久精品综合| 我的亚洲天堂| 啦啦啦啦在线视频资源| 午夜福利一区二区在线看| 99精国产麻豆久久婷婷| 亚洲伊人久久精品综合| 亚洲婷婷狠狠爱综合网| 国产欧美日韩一区二区三区在线| 亚洲久久久国产精品| 国产女主播在线喷水免费视频网站| 操美女的视频在线观看| 精品少妇内射三级| 亚洲精品国产av蜜桃| 高清视频免费观看一区二区| 日韩av不卡免费在线播放| 高清av免费在线| 免费日韩欧美在线观看| 男人添女人高潮全过程视频| 国产伦人伦偷精品视频| 各种免费的搞黄视频| 一区二区三区四区激情视频| 晚上一个人看的免费电影| 亚洲,欧美,日韩| 最新的欧美精品一区二区| 亚洲欧洲日产国产| 亚洲欧美一区二区三区国产| 亚洲人成网站在线观看播放| 十八禁人妻一区二区| 国产成人啪精品午夜网站| 国产精品熟女久久久久浪| 亚洲成人免费av在线播放| 日韩,欧美,国产一区二区三区| 青草久久国产| 91老司机精品| 如何舔出高潮| 亚洲精品美女久久av网站| 国产精品成人在线| 欧美日韩视频高清一区二区三区二| 国产日韩欧美视频二区| 日韩一区二区视频免费看| 永久免费av网站大全| 亚洲国产欧美在线一区| 美女国产高潮福利片在线看| 秋霞伦理黄片| 午夜日本视频在线| 多毛熟女@视频| 亚洲一区中文字幕在线| 欧美黑人精品巨大| 麻豆av在线久日| 欧美少妇被猛烈插入视频| 老司机在亚洲福利影院| 国产又色又爽无遮挡免| 亚洲人成电影观看| 精品一区二区三区av网在线观看 | a级片在线免费高清观看视频| a级毛片在线看网站| 熟女少妇亚洲综合色aaa.| 国产精品熟女久久久久浪| 搡老岳熟女国产| 大话2 男鬼变身卡| 亚洲国产欧美在线一区| 亚洲欧美一区二区三区国产| 18禁国产床啪视频网站| 91aial.com中文字幕在线观看| 亚洲情色 制服丝袜| 建设人人有责人人尽责人人享有的| 久久av网站| 两性夫妻黄色片| 黑人欧美特级aaaaaa片| 国产精品久久久久久精品古装| 捣出白浆h1v1| 精品一区二区三区四区五区乱码 | 男人爽女人下面视频在线观看| 国产精品一国产av| 久久久久久久精品精品| 久久人人97超碰香蕉20202| 看十八女毛片水多多多| 日韩一卡2卡3卡4卡2021年| 女人爽到高潮嗷嗷叫在线视频| 国产亚洲精品第一综合不卡| 高清黄色对白视频在线免费看| 国产97色在线日韩免费| 一区二区三区四区激情视频| 国产人伦9x9x在线观看| 精品福利永久在线观看| 国产高清不卡午夜福利| 天天躁夜夜躁狠狠久久av| 亚洲国产av新网站| 最新在线观看一区二区三区 | 美女视频免费永久观看网站| 国产又色又爽无遮挡免| a级毛片在线看网站| 日韩av免费高清视频| 男男h啪啪无遮挡| 亚洲欧美激情在线| 国产成人系列免费观看| 免费在线观看完整版高清| 中文字幕人妻丝袜制服| 欧美日韩视频精品一区| 亚洲精品乱久久久久久| 久久97久久精品| 天天躁日日躁夜夜躁夜夜| 中文字幕人妻熟女乱码| 亚洲欧洲日产国产| 中文乱码字字幕精品一区二区三区| 亚洲国产精品国产精品| 看免费成人av毛片| 91国产中文字幕| 国产亚洲午夜精品一区二区久久| 久久综合国产亚洲精品| 王馨瑶露胸无遮挡在线观看| 蜜桃国产av成人99| 精品视频人人做人人爽| 亚洲成人av在线免费| av免费观看日本| 在线观看国产h片| 一本大道久久a久久精品| 美女福利国产在线| 美女福利国产在线| 国产一级毛片在线| 人妻 亚洲 视频| 一级片免费观看大全| 欧美日韩亚洲国产一区二区在线观看 | 国产视频首页在线观看| 热99国产精品久久久久久7| 男女午夜视频在线观看| 90打野战视频偷拍视频| 少妇人妻 视频| 丝袜美足系列| 国产欧美日韩一区二区三区在线| 80岁老熟妇乱子伦牲交| 99精国产麻豆久久婷婷| av.在线天堂| 成人国产麻豆网| 日本爱情动作片www.在线观看| 国产精品免费视频内射| 人人妻人人澡人人爽人人夜夜| 电影成人av| xxx大片免费视频| 日本av手机在线免费观看| 免费在线观看视频国产中文字幕亚洲 | 自线自在国产av| 久久久久精品人妻al黑| 亚洲激情五月婷婷啪啪| 国产精品人妻久久久影院| 久久精品国产综合久久久| 在线观看www视频免费| 日韩制服丝袜自拍偷拍| 日本黄色日本黄色录像| 两性夫妻黄色片| 亚洲 欧美一区二区三区| 国产无遮挡羞羞视频在线观看| 两个人免费观看高清视频| 你懂的网址亚洲精品在线观看| 操出白浆在线播放| 日韩成人av中文字幕在线观看| 久久精品久久久久久久性| 亚洲成人国产一区在线观看 | 免费在线观看黄色视频的| 久久综合国产亚洲精品| 久久久精品区二区三区| 国产成人啪精品午夜网站| 免费观看性生交大片5| 下体分泌物呈黄色| 男女无遮挡免费网站观看| 亚洲av欧美aⅴ国产| 又粗又硬又长又爽又黄的视频| 高清欧美精品videossex| 熟妇人妻不卡中文字幕| 91精品三级在线观看| 亚洲精品,欧美精品| 国产日韩一区二区三区精品不卡| 国产亚洲av高清不卡| 国产精品久久久人人做人人爽| 亚洲av欧美aⅴ国产| www.精华液| 亚洲精品日本国产第一区| 日本欧美视频一区| 电影成人av| 黑人猛操日本美女一级片| 少妇 在线观看| 男女无遮挡免费网站观看| 亚洲精品视频女| 国产亚洲最大av| 秋霞在线观看毛片| 成年动漫av网址| 久久久国产精品麻豆| bbb黄色大片| 色视频在线一区二区三区| 自线自在国产av| 在线观看免费视频网站a站| 精品国产露脸久久av麻豆| 久久精品久久久久久噜噜老黄| 99香蕉大伊视频| 中国国产av一级| 亚洲精品第二区| 18禁国产床啪视频网站| 久久天躁狠狠躁夜夜2o2o | 大香蕉久久成人网| 精品一区在线观看国产| 国产精品嫩草影院av在线观看| 久久狼人影院| 老鸭窝网址在线观看| www.精华液| 97精品久久久久久久久久精品| 成人午夜精彩视频在线观看| www.av在线官网国产| 久久这里只有精品19| 欧美乱码精品一区二区三区| 国产亚洲精品第一综合不卡| 男女之事视频高清在线观看 | 老鸭窝网址在线观看| 亚洲一级一片aⅴ在线观看| 国产免费一区二区三区四区乱码| 男人舔女人的私密视频| 日韩熟女老妇一区二区性免费视频| 中文字幕色久视频| 美女国产高潮福利片在线看| 欧美黑人精品巨大| 国产极品天堂在线| 久久综合国产亚洲精品| 亚洲国产成人一精品久久久| 老汉色av国产亚洲站长工具| 日本色播在线视频| 在线观看国产h片| 热re99久久精品国产66热6| av视频免费观看在线观看| 欧美日韩成人在线一区二区| 国产激情久久老熟女| a级毛片黄视频| 国产日韩欧美亚洲二区| 精品国产一区二区三区久久久樱花| 亚洲精品日韩在线中文字幕| 一级毛片电影观看| 美女视频免费永久观看网站| 国产成人精品福利久久| 国产精品秋霞免费鲁丝片| 人妻人人澡人人爽人人| 巨乳人妻的诱惑在线观看| 国产精品 国内视频| 多毛熟女@视频| 国产精品一区二区精品视频观看| 亚洲欧美成人精品一区二区| 91aial.com中文字幕在线观看| 欧美日韩视频高清一区二区三区二| 国产精品一二三区在线看| 精品少妇黑人巨大在线播放| 哪个播放器可以免费观看大片| 午夜激情久久久久久久| 亚洲精品日韩在线中文字幕| 精品国产露脸久久av麻豆| 少妇人妻 视频| 久久韩国三级中文字幕| 丝瓜视频免费看黄片| 午夜福利免费观看在线| 在线观看三级黄色| 水蜜桃什么品种好| 亚洲视频免费观看视频| 亚洲美女黄色视频免费看| 国产人伦9x9x在线观看| 亚洲精品国产av蜜桃| 日韩熟女老妇一区二区性免费视频| 日韩中文字幕视频在线看片| av不卡在线播放| 尾随美女入室| 亚洲欧美成人综合另类久久久| 天天躁日日躁夜夜躁夜夜| 新久久久久国产一级毛片| 大话2 男鬼变身卡| 欧美久久黑人一区二区| 悠悠久久av| 亚洲第一青青草原| bbb黄色大片| 精品少妇内射三级| 免费少妇av软件| av网站免费在线观看视频| 丰满迷人的少妇在线观看| 精品国产超薄肉色丝袜足j| 99re6热这里在线精品视频| 一级片免费观看大全| 国产成人啪精品午夜网站| 9热在线视频观看99| 久久精品熟女亚洲av麻豆精品| 大话2 男鬼变身卡| 免费久久久久久久精品成人欧美视频| 五月天丁香电影| 久久婷婷青草| 国产黄色免费在线视频| 日日摸夜夜添夜夜爱| 在线观看免费视频网站a站| 国产黄频视频在线观看| 日韩视频在线欧美| 亚洲三区欧美一区| 欧美黑人精品巨大| 欧美黄色片欧美黄色片| 国产伦人伦偷精品视频| 亚洲欧洲日产国产| 免费黄色在线免费观看| 亚洲第一青青草原| 国产一区二区三区av在线| 精品少妇内射三级| 国产xxxxx性猛交| 大香蕉久久成人网| 中文字幕高清在线视频| 精品国产一区二区三区四区第35| 亚洲av国产av综合av卡| 国精品久久久久久国模美| 超碰97精品在线观看| 老司机深夜福利视频在线观看 | av片东京热男人的天堂| 欧美人与性动交α欧美软件| 99re6热这里在线精品视频| 爱豆传媒免费全集在线观看| 久热爱精品视频在线9| 欧美老熟妇乱子伦牲交| 国产有黄有色有爽视频| 国产成人欧美在线观看 | 永久免费av网站大全| 久久精品久久精品一区二区三区| 日韩视频在线欧美| 人人妻人人澡人人看| 水蜜桃什么品种好| 99re6热这里在线精品视频| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩国产mv在线观看视频| 午夜免费男女啪啪视频观看| 免费观看a级毛片全部| 国产成人a∨麻豆精品| 男女无遮挡免费网站观看| 中文字幕人妻丝袜制服| 欧美黑人欧美精品刺激| 中文字幕人妻熟女乱码| 2021少妇久久久久久久久久久| 久久毛片免费看一区二区三区| 亚洲美女视频黄频| 好男人视频免费观看在线| 亚洲情色 制服丝袜| 亚洲精品国产av蜜桃| 一级片免费观看大全| 亚洲av电影在线观看一区二区三区| 亚洲精品国产av蜜桃| 色网站视频免费| 激情五月婷婷亚洲| 高清在线视频一区二区三区| 汤姆久久久久久久影院中文字幕| 高清欧美精品videossex| 久久鲁丝午夜福利片| 男女边吃奶边做爰视频| 亚洲专区中文字幕在线 | 色吧在线观看| 午夜激情av网站| 日韩,欧美,国产一区二区三区| 久久久久视频综合| 久久精品久久久久久久性| 天堂中文最新版在线下载| 久久综合国产亚洲精品| 免费黄频网站在线观看国产| 久久性视频一级片| 又黄又粗又硬又大视频| 欧美激情 高清一区二区三区| 色婷婷av一区二区三区视频| 亚洲精品久久午夜乱码| 免费黄频网站在线观看国产| 成人手机av| 老司机亚洲免费影院| 中文字幕人妻丝袜一区二区 | 一级爰片在线观看| 国产一区二区三区av在线| 国产老妇伦熟女老妇高清| 九草在线视频观看| 丰满乱子伦码专区| videos熟女内射| 天天操日日干夜夜撸| 日本av免费视频播放| 亚洲成国产人片在线观看| 韩国精品一区二区三区| 亚洲欧美一区二区三区黑人| 国产日韩一区二区三区精品不卡| 在线观看一区二区三区激情| 看非洲黑人一级黄片| 一边摸一边做爽爽视频免费| 欧美日韩一区二区视频在线观看视频在线| 欧美少妇被猛烈插入视频| 久久久国产欧美日韩av| 日韩伦理黄色片| 777久久人妻少妇嫩草av网站| 老司机影院毛片| 欧美日韩视频高清一区二区三区二| 欧美av亚洲av综合av国产av | 中文字幕精品免费在线观看视频| 国产在线视频一区二区| svipshipincom国产片| 久久久国产欧美日韩av| 精品国产乱码久久久久久小说| 日韩 亚洲 欧美在线| 亚洲av日韩精品久久久久久密 | 免费在线观看黄色视频的| 一级黄片播放器| 欧美最新免费一区二区三区| 搡老乐熟女国产| 天天影视国产精品| 国产欧美日韩综合在线一区二区| 麻豆乱淫一区二区| 亚洲av日韩在线播放| 99久久综合免费| 男女床上黄色一级片免费看| 熟女少妇亚洲综合色aaa.| 亚洲av在线观看美女高潮| 国产精品无大码| 97人妻天天添夜夜摸| 在线天堂最新版资源| 久久99热这里只频精品6学生| 国产午夜精品一二区理论片| 国产精品人妻久久久影院| 伊人久久国产一区二区| 国产精品亚洲av一区麻豆 | 观看美女的网站| 久久鲁丝午夜福利片| 嫩草影院入口| 高清黄色对白视频在线免费看| 九色亚洲精品在线播放| 国产精品成人在线| 亚洲国产成人一精品久久久| 十分钟在线观看高清视频www| 欧美在线黄色| 亚洲精品国产色婷婷电影| 男女床上黄色一级片免费看| 大片免费播放器 马上看| 狠狠婷婷综合久久久久久88av| 成年人午夜在线观看视频| 日韩一本色道免费dvd| 丰满饥渴人妻一区二区三| 少妇被粗大猛烈的视频| tube8黄色片| 国产亚洲欧美精品永久| 蜜桃在线观看..| 色精品久久人妻99蜜桃| 黄片小视频在线播放| 在线观看三级黄色| 久久久国产欧美日韩av| 亚洲伊人久久精品综合| 欧美日韩亚洲高清精品| av网站免费在线观看视频| 中文字幕高清在线视频| 国产伦理片在线播放av一区| 99久久精品国产亚洲精品| 97人妻天天添夜夜摸| 成人亚洲欧美一区二区av| 成人漫画全彩无遮挡| 一个人免费看片子| 欧美日韩视频精品一区| www日本在线高清视频| 波野结衣二区三区在线| 久久热在线av| 国产日韩欧美在线精品| 国产免费现黄频在线看| 99香蕉大伊视频| 久久精品久久久久久久性| 亚洲av日韩精品久久久久久密 | tube8黄色片| 亚洲一码二码三码区别大吗| 亚洲 欧美一区二区三区| 欧美黑人精品巨大| 69精品国产乱码久久久| 国产精品久久久久久精品电影小说| 国产一区二区三区综合在线观看| 多毛熟女@视频| 在线观看人妻少妇| 国产精品一国产av| 国产精品一区二区精品视频观看| 国产精品av久久久久免费| 亚洲,一卡二卡三卡| 欧美激情 高清一区二区三区| 一级片'在线观看视频| av在线app专区| 69精品国产乱码久久久| 国产精品国产三级国产专区5o| 亚洲av日韩在线播放| 亚洲国产最新在线播放| 亚洲国产日韩一区二区| 国产精品一区二区在线不卡| 国产xxxxx性猛交| 国产国语露脸激情在线看| 久久精品久久久久久久性| 一级,二级,三级黄色视频| 亚洲精品国产av成人精品| 亚洲av日韩在线播放| 亚洲欧美精品综合一区二区三区| 国产精品一国产av| 精品亚洲乱码少妇综合久久| av网站在线播放免费| av视频免费观看在线观看| 两个人看的免费小视频| 久久久国产精品麻豆| 中文字幕精品免费在线观看视频| 久久免费观看电影| 亚洲精品在线美女| 亚洲国产av新网站| 男女无遮挡免费网站观看| 久久午夜综合久久蜜桃| 大香蕉久久网| 欧美人与性动交α欧美软件| 麻豆av在线久日| 亚洲国产欧美日韩在线播放| 少妇人妻久久综合中文| h视频一区二区三区| 亚洲国产精品一区三区| 久久久久久久久久久免费av| 在现免费观看毛片| 超色免费av| 韩国精品一区二区三区| 又粗又硬又长又爽又黄的视频| 免费看不卡的av| 啦啦啦在线免费观看视频4| 国产熟女欧美一区二区| 无限看片的www在线观看| 欧美日韩综合久久久久久| av在线播放精品| 老司机亚洲免费影院| 丝袜喷水一区| 在线观看免费午夜福利视频| 水蜜桃什么品种好| 永久免费av网站大全| 亚洲自偷自拍图片 自拍| 老汉色av国产亚洲站长工具| 性高湖久久久久久久久免费观看| 亚洲情色 制服丝袜| 国产一卡二卡三卡精品 | 亚洲情色 制服丝袜| 国产毛片在线视频| 91老司机精品| 久久鲁丝午夜福利片| 久久 成人 亚洲| 美国免费a级毛片| 人人妻人人爽人人添夜夜欢视频| 激情五月婷婷亚洲| 成年动漫av网址| 可以免费在线观看a视频的电影网站 | 建设人人有责人人尽责人人享有的| 亚洲精品自拍成人| 国产精品熟女久久久久浪| 人成视频在线观看免费观看| 一边摸一边抽搐一进一出视频| 国产精品 国内视频| 99热全是精品| 日本av免费视频播放| 欧美 日韩 精品 国产| 日韩不卡一区二区三区视频在线| 黑人巨大精品欧美一区二区蜜桃| 女人精品久久久久毛片| 国产xxxxx性猛交| 欧美另类一区| 欧美激情 高清一区二区三区| xxxhd国产人妻xxx| 国产免费现黄频在线看| 日韩免费高清中文字幕av| 最近最新中文字幕大全免费视频 | 欧美亚洲日本最大视频资源| 操出白浆在线播放| 三上悠亚av全集在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 成年人免费黄色播放视频| 不卡av一区二区三区| 中文字幕人妻熟女乱码| av视频免费观看在线观看| 亚洲欧美一区二区三区久久| 欧美日韩av久久| 久久99精品国语久久久| 少妇精品久久久久久久| 一区二区三区四区激情视频| 久久久亚洲精品成人影院| 亚洲av成人不卡在线观看播放网 | 亚洲国产欧美在线一区| 大香蕉久久网| 国产亚洲av高清不卡| 中国三级夫妇交换| 日韩,欧美,国产一区二区三区| 老司机在亚洲福利影院| 69精品国产乱码久久久| 极品少妇高潮喷水抽搐| 美女扒开内裤让男人捅视频| 久久毛片免费看一区二区三区| 丰满乱子伦码专区| 国产精品国产av在线观看| 久久久久人妻精品一区果冻| 尾随美女入室| 少妇猛男粗大的猛烈进出视频| 叶爱在线成人免费视频播放| 男人爽女人下面视频在线观看| 亚洲av中文av极速乱| 妹子高潮喷水视频| 一区二区av电影网| 精品国产一区二区三区四区第35| av天堂久久9| 制服人妻中文乱码| 高清视频免费观看一区二区| 毛片一级片免费看久久久久| 你懂的网址亚洲精品在线观看| 少妇 在线观看| 各种免费的搞黄视频| 国产一区二区激情短视频 | 男男h啪啪无遮挡| 亚洲七黄色美女视频| 国语对白做爰xxxⅹ性视频网站| 亚洲精品国产区一区二| 波多野结衣一区麻豆| 国产乱人偷精品视频| 欧美日韩亚洲国产一区二区在线观看 | 国产精品av久久久久免费| 午夜激情av网站| 亚洲一卡2卡3卡4卡5卡精品中文| 日本91视频免费播放|