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

    冰下爆破預(yù)防冰壩的理論探索及實踐

    2020-05-14 02:05:26楊開林郭新蕾郭永鑫彭旭明吳煜楠
    水利學(xué)報 2020年2期
    關(guān)鍵詞:冰蓋水深融化

    楊開林,郭新蕾,王 濤,郭永鑫,彭旭明,吳煜楠

    (1.流域水循環(huán)模擬與調(diào)控國家重點實驗室,中國水利水電科學(xué)研究院,北京 100038;2.黑龍江省防汛抗旱保障中心,黑龍江 哈爾濱 150001)

    1 研究背景

    我國高寒地區(qū)冬季江河冰情嚴重,開河時容易發(fā)生冰壩洪水災(zāi)害,如我國的黃河、黑龍江、松花江[1-2]。冰壩是開河期由大量的流冰在河道中局部堆積而成的[3],一旦冰壩形成,上游水位短期內(nèi)往往接近或超過歷史上最高洪水位[4-5],可導(dǎo)致水位壅高漫出河堤,淹沒農(nóng)田房屋,使沿岸水工建筑物和構(gòu)筑物遭到破壞。采用爆破作業(yè)消除冰壩洪水危險是行之有效的非工程措施。

    爆破作業(yè)清理冰的方法可以追溯到200多年前[6-8]。Fonstad等[8]通過冰蓋下爆破試驗提出了經(jīng)驗公式:

    式中:d為炸藥放置在冰下的最優(yōu)水深,m;Wopt為最佳裝藥量,kg;h為冰蓋厚度,m;D為爆炸后爆破坑直徑,m。Mellor[7]對以前資料進行了整編,給出的最佳裝藥量和爆破坑直徑與冰厚的經(jīng)驗公式:

    他們的研究表明爆破的效果與炸藥的類型無關(guān),炸藥筒可以放置在冰上的鉆孔中,然后慢慢順進冰下去,之后把它綁到木制的橫擋上,橫跨在鉆孔上,如圖1所示。

    圖1 炸藥筒放置示意圖

    在我國黃河開河期,一直采用飛機投彈破除危險性冰壩,以防止嚴重冰壩洪水發(fā)生。黑龍江也曾經(jīng)采用爆破破除危險性冰壩。近年來,已經(jīng)開展了大量的黃河冰蓋下的爆炸破冰試驗。梁向前等[9-10]試驗研究了冰面接觸聚能爆炸和冰下水中爆炸兩種破冰方式,歸納提出了冰蓋下單藥包和多藥包爆破水中沖擊波壓力經(jīng)驗公式。

    爆破消除冰壩洪水風(fēng)險作業(yè)可劃分為兩種:爆破破除冰壩和爆破預(yù)防冰壩,前者是在開河時冰壩已經(jīng)形成條件下被動采取的應(yīng)急措施,后者是在開河前主動采取的預(yù)防措施。

    爆破破除冰壩需要考慮4方面問題[7]:時間、爆破效果、破冰順序和冰的輸運。一旦冰壩形成,河/江道過水?dāng)嗝婢蜁乐囟氯瑢?dǎo)致上游水位迅速上升,嚴重時將發(fā)生漫堤,甚至決堤事故,因此,在冰壩剛剛形成的幾小時之內(nèi)實施爆破破除冰壩是最好的時間。為了獲得最大的爆破效果,炸藥應(yīng)該安裝在冰壩下,但是如果不能夠把人員放在冰壩上,這將很難實現(xiàn)。如果炸藥不能安裝在冰壩的下方,那就應(yīng)該將炸藥安裝在有空洞或裂縫的地方,炸藥一旦安裝完畢,同時引爆炸藥可獲得最大的效果。破冰順序應(yīng)該從冰壩趾部開始向上游延伸,此時水流可以將破碎的冰塊向下游輸運。問題在于,在冰壩剛剛形成的幾小時之內(nèi),水位還會不斷的上漲,冰壩的穩(wěn)定性一般還不能滿足人員和設(shè)備的要求,實施冰下爆破存在很大的風(fēng)險。

    由于冰壩形成后采取冰下爆破作業(yè)存在很大的風(fēng)險,因此,在水利部公益性行業(yè)科研專項“黑龍江冰情預(yù)報及災(zāi)害防治研究”資助下,把冰、水、熱動力學(xué)理論與爆破效果相結(jié)合,開展了冰蓋下爆破預(yù)防冰壩的理論探索和實踐。

    2 基礎(chǔ)理論

    為了使得冰蓋下爆破預(yù)防冰壩的效果最大化,首先需要解決5個關(guān)鍵科學(xué)技術(shù)問題:爆破時間、爆破點、爆破效果、爆破孔的布置和溫度場的作用。

    2.1 爆破時間冰下爆破后,爆破孔周圍冰蓋破碎,形成爆破坑,爆破后的碎冰除部分像雪花一樣散落在爆破坑周圍外,大部分堆積在爆破坑內(nèi),形成如圖2所示水浸冰。當(dāng)爆破后氣溫在0.0℃以下,則爆破坑內(nèi)水會重新結(jié)冰形成冰蓋,并把碎冰塊凍結(jié)在一起,使得爆破效果大打折扣。因此,實施冰下爆破的最佳時間是氣溫由負溫穩(wěn)定轉(zhuǎn)為正溫的初期,以使爆破河段及下游冰蓋有足夠長的時間融化,達到爆破間接效果最大化。例如在開河前幾個星期鉆孔作業(yè),這時冰蓋承載能力較大,作業(yè)人員掉落冰窟窿的風(fēng)險較小。

    2.2 爆破點冰壩易發(fā)點一般是河道變窄、彎道、淺灘、多分叉河段、跨河建筑物及水深由淺到深過渡區(qū),當(dāng)這些河段下游具有完整冰蓋時,則上游來冰就可能在冰蓋前堆積形成冰壩[2,11]。

    彎道的水力學(xué)特點是:在水流的離心力作用下,凹入的河岸受到侵蝕,侵蝕下來的泥沙,通過彎遭環(huán)流,搬運到對岸沉積,逐漸形成圓弧狀坡緩的凸岸。凹岸通常水深、流速較大;凸岸坡度緩和,近岸水流緩慢。在寒冷的冬季,凸岸一側(cè)封河期易于形成岸冰,并向河心發(fā)展,縮窄河道,當(dāng)上游大量浮冰來臨,則會形成冰橋,成為冰蓋的起點并向上游發(fā)展。當(dāng)“武開河”時,即開河順序是由上游到下游,由于彎道冰蓋較厚,不易破損,大量水面流冰易于在這里上擠下潛形成冰壩。

    淺灘河道的水力學(xué)特點是:上游河道底坡較陡,水深較淺,平均流速較大,而淺灘河寬較大,水深淺,平均流速小。在封河期,在淺灘流速<0.3 m/s的水面形成靜態(tài)冰[12],在上游水面漂浮大量來冰的情況下形成封河冰蓋。如果上游河道水流急,短時間內(nèi)水面不能冰封,則在寒冷空氣的作用下,通過空氣與水面的熱交換,水中將產(chǎn)生大量的冰花和水內(nèi)冰流到下游淺灘,通過下潛堆積在冰蓋前緣下面形成冰塞,從而導(dǎo)致上游水位升高,流速減小,冰蓋向上游發(fā)展。

    圖2 黑龍江爆破坑照片

    如果在開河前在冰壩易發(fā)點實施冰蓋下爆破,迫使這些河段冰蓋或者冰塞在上游開河前融化消除,則可能防止冰壩洪水的發(fā)生。

    2.3 爆破效果冰下爆破的效果可劃分為直接效果和間接效果。直接效果的常用評估指標是用藥量與形成爆破坑的面積大小。間接效果是指爆破坑形成后,由于出現(xiàn)水與空氣的直接接觸,增強了爆破坑處水與大氣、冰與水的熱交換,不僅使得爆破坑內(nèi)冰融化,而且使得爆破坑周圍冰蓋加速融化,進而使得自由水面不斷增加,同時,由于冰蓋的融化,原來的冰域也轉(zhuǎn)變?yōu)樗颍钤黾?,槽蓄量增加,有利于開江時冰排的順利通過。在黑龍江進行的冰蓋下爆破結(jié)果證明,在水深較淺的地方開展爆破不能很好發(fā)揮炸藥的能力,在水深較大的地方可以提高爆破的直接效果[13]。下面將通過分析冰蓋厚度融化速率與氣溫、水溫、流速、水深的關(guān)系,確定爆破孔的最佳位置,以提高冰下爆破的間接效果。

    2.3.1 冰蓋厚度融化速率與氣溫、水溫、流速、水深的關(guān)系 在沒有雪蓋的條件下,冰蓋厚度融化速率與氣溫和水溫的關(guān)系是[14]

    式中:h為冰蓋厚度,m;t為時間,s;TW為水溫,℃;Ta為空氣的日平均氣溫,℃;Ts為冰蓋頂?shù)臏囟龋?;Tm為冰點溫度,℃,Tm≈ 0.0℃;ρi為冰的密度,ρi=915~917 kg/m3;Li為單位質(zhì)量冰的潛熱,當(dāng)水結(jié)冰時,Li=3.33×105J/kg的潛熱被釋放;hai為空氣與冰的熱交換系數(shù),W/(m2·℃),是一個綜合考慮太陽輻射、空氣的溫度和濕度、風(fēng)速、氣壓、云及降雨和降雪等因素的參數(shù);hwi為水與冰的熱交換系數(shù),W/(m2·℃)。觀察式(3),右邊第一項是冰表面與大氣的熱交換產(chǎn)生的冰蓋厚度減小的速率,右邊第二項是冰底面水與冰的熱交換產(chǎn)生的冰蓋厚度減小的速率,顯然當(dāng)氣溫和水溫均為正時,冰蓋厚度上下都會發(fā)生融化減小。

    熱交換系數(shù)hwi與流速和水力半徑的關(guān)系可描述為[3]:

    式中:U為沿流向的水深平均流速,m/s;H為過水?dāng)嗝鏅M向坐標y處的水深,m;R為水力半徑,m。假設(shè)河道單位長度上冰蓋與河底濕周相同,則R=0.5H,這時式(4)可改寫為:

    丹麥DHI的流速橫向分布公式[15]:

    式中:Um為整個河道過水?dāng)嗝娴钠骄魉?,m/s;Hm為過水?dāng)嗝嫫骄?,m;α為常數(shù)。顯然,在河道過水?dāng)嗝嫔?,水深H越大,流速U越大,當(dāng)然,一些局部深坑除外。

    把式(6)代入式(5)消去H得:

    觀察式(7)、式(6)、式(3)可得重要結(jié)論1:當(dāng)水溫為正時,在河道同一過水?dāng)嗝嫔?,流速與水深成正比,冰蓋厚度減小的速率隨流速的增加而增加。

    2.3.2 爆破坑側(cè)壁冰融化的速率與水溫和流速的關(guān)系 類似地,當(dāng)采用柱坐標系時,爆破坑側(cè)壁冰融化的速率與水溫的關(guān)系可描述為:

    式中:r為水深z處側(cè)壁到坐標原點的距離,m;hwi,z為爆破坑側(cè)壁的水與冰的熱交換系數(shù),W/(m2·℃),假設(shè)hwi,z與流速具有式(7)關(guān)系,則:

    式中:v為爆破坑內(nèi)側(cè)壁附近的流速,m/s;c和β為待定系數(shù)。

    對于寬深比較大的河道,冰蓋下流速的縱向分布特點是:冰蓋底和河床上的流速為零,最大流速在冰蓋和河床之間,最大流速與冰蓋及最大流速與河床之間近似呈對數(shù)規(guī)律分布[16],因此,在爆破坑尺寸相對河寬較小的情況下,冰蓋下的較高流速將拖拽爆破坑內(nèi)水流運動,增強了質(zhì)量和熱量的對流交換,水溫趨于均勻,在坑內(nèi)形成上下和左右環(huán)流,且迎水側(cè)v大于背水側(cè),下部v大于上部。假設(shè)爆破坑內(nèi)水溫摻混均勻,則從式(8)和式(9)可得重要結(jié)論2:當(dāng)水溫為正時,爆破坑側(cè)壁范圍內(nèi)冰融化的速率隨水深的增加而增加,迎水側(cè)融化速率大于背水側(cè)。

    圖3是美國陸軍寒冷地區(qū)研究與工程實驗所(CRREL)在冷凍水槽中觀測到的直徑2.54 cm鉆孔邊緣的融化現(xiàn)象,迎水側(cè)融化速率明顯大于背水側(cè),下部大于上部。根據(jù)結(jié)論2,這種鉆孔邊緣的冰融化形狀應(yīng)是由于鉆孔側(cè)壁流速分布不同產(chǎn)生的。類似地,可以得出,爆破坑融化的形狀類似吊掛的頂部開口的不對稱的喇叭型。

    2.3.3 爆破孔的最佳位置 由上面結(jié)論1和2,當(dāng)爆破孔位于流速較大的地方時,爆破坑周圍冰蓋融化的速率較大,即冰下爆破間接效果較大,同時,對同一河道過水?dāng)嗝妫钤酱蟮牡胤?,水深平均流速也大,從中可得重要結(jié)論3:爆破孔的最佳位置是河道冰壩易發(fā)點的主槽深泓線,那里是水深流速大的地方。

    2.4 爆破孔的布置冰壩易發(fā)點的河道一般為復(fù)式河道,可劃分為主槽和漫灘,主槽水深和流速較大,漫灘水深和流速較小。在灘槽交界面附近,水深和流速發(fā)生急劇變化,水流紊動強度大,形成復(fù)雜的橫向二次流和螺旋流,使灘槽水體發(fā)生大量的質(zhì)量交換和對流熱量交換[17-18]。

    在黑龍江現(xiàn)場的觀察發(fā)現(xiàn),在冰下爆破實施江段,開河前會發(fā)生下游相當(dāng)長江段的冰蓋突然崩解的現(xiàn)象,這表明爆破提高的水溫對下游冰蓋融化有很大作用。為了使得冰下爆破預(yù)防冰壩的成功性和投入的經(jīng)濟性,理想的目標是通過爆破形成具有一定規(guī)模的溫度場,不僅消除爆破河段的冰蓋,而且消除下游相當(dāng)長河道的冰蓋。換句話說,冰下爆破預(yù)防冰壩不僅要考慮爆破的直接效果,更要考慮冰、水、熱動力學(xué)作用的間接效果。為此,基于復(fù)式河道的流動和熱交換特點及爆破坑周圍冰蓋融化的特點,可以采用下述爆破孔布置方式:

    (1)根據(jù)河寬,只在深泓線附近沿流向布設(shè)幾列爆破孔,例如1~3 排,以形成一定規(guī)模的溫度場,然后利用主槽和漫灘的質(zhì)量和熱量對流以促使兩側(cè)漫灘和下游冰蓋的解凍;

    (2)考慮爆破間接效果的作用,爆破孔之間的間距以大于相鄰爆破坑的半徑之和為宜,沿縱向的相鄰爆破孔的間距宜大于橫向相鄰爆破孔的間距,其原因一是縱向流速大,二是上游爆破坑吸收的熱量將傳遞到下游,增強了下游爆破坑周圍冰蓋的融化速率。

    2.5 溫度場作用在按照上述爆破時間和爆破孔布置實施冰下爆破后,爆破坑內(nèi)碎冰和周圍冰蓋將首先融化,然后爆破坑自由水面相互連通形成清溝,形成局部溫度場,隨后,隨著清溝兩側(cè)冰蓋的融化,形成全斷面自由水面的溫度場。

    下面分析溫度場形成后,場水溫與長度的關(guān)系,場水溫對下游冰蓋下水溫的影響長度,下游冰蓋厚度與場水溫和時間的關(guān)系等。為使問題簡化,假設(shè):(1)過水?dāng)嗝鏋榫匦?;?)溫度場明流和下游冰蓋下流動為均勻流。

    2.5.1 溫度場水溫與氣溫和溫度場長度的關(guān)系 冰下爆破河段形成溫度場后,忽略溫度擴散的影響,則一維熱對流方程可描述[14]為:

    圖3 冰蓋上鉆孔周圍的融化狀態(tài)

    式中:x為沿流向的距離,m;U1為溫度場的流速,m/s;H1為溫度場的水深,m;ρ為水的密度,一般取ρ=1000 kg/m3;Cp為水的比熱,J/kg℃,在0℃時,Cp=4217.7 J/kg℃;hwa為水向空氣的熱交換系數(shù),是一個綜合考慮太陽輻射、空氣的溫度和濕度、風(fēng)速、氣壓、云及降雨和降雪等因素的參數(shù),對北美地區(qū),大氣和水交界面的熱交換系數(shù)hwa=20 W/(m2·℃)[3]。

    采用特征線方法,式(10)可改寫為:

    式中x為水流質(zhì)點離溫度場進口的距離,m。

    在一天時間內(nèi),日平均氣溫為常數(shù),求解式(11)常微分得:

    式中TW0為t=0時刻溫度場進口水溫,℃。

    對于靜水,x=0,式(13)表示水溫隨時間的增加成自然指數(shù)增加。在爆破初期的爆破坑內(nèi)的水體可以認為是靜水,這時水溫隨氣溫的增加而快速增加,但隨著碎冰的減少,冰蓋下流動將使得爆破坑內(nèi)水體流動,把熱傳遞到下部,導(dǎo)致爆破坑內(nèi)的水溫下降。

    當(dāng)流速不為零時,把式(14)代入式(13)可得溫度場水溫與距離x的關(guān)系:

    當(dāng)溫度場長度為x1,則對應(yīng)水溫:

    式中:TW1為溫度場出口水溫或下游冰蓋前緣流入水溫,℃。

    在一般情況下,x1<104m,即成立,式(16)可改寫為:

    在自然環(huán)境下,冰蓋下水溫TW0非常接近Tm,即TW0≈ 0℃[19],因此,可得重要結(jié)論4:溫度場出口水溫TW1與溫度場長度、氣溫及水與空氣熱交換系數(shù)的積TaHWax1成正比,與單寬流量成反比。

    2.5.2 溫度場水溫對下游冰蓋下水溫的影響長度 對于溫度場下游冰蓋下的流動,一維熱對流方程可描述[14]為:

    式中:U2為溫度場下游冰蓋下流速,m/s;H2為冰蓋下水深,m,即冰蓋底高程與江底高程之差。

    與節(jié)2.5.1類似可得:

    假設(shè)x=x2時,由于水與冰的熱交換,冰蓋下水溫從TW1下降到溫度場進口水溫TW0,則稱L=x2-x1為溫度場水溫對下游冰蓋下水溫的影響長度,這時:

    當(dāng)L<104m時,,上式改寫為:

    把式(5)代入式(21):

    由上式可得重要結(jié)論5:場水溫、冰蓋下水深和流速越大,場水溫對下游冰蓋下水溫的影響長度L越大。

    當(dāng)取ρ=1000 kg/m3和Cp=4217.7 J/kg℃,則可得:

    2.5.3 溫度場下游冰蓋厚度與溫度場水溫、氣溫和時間的關(guān)系 式(3)也適合于溫度場下游冰蓋厚度的變化,在沒有雪蓋的條件下,冰蓋厚度的變化與氣溫和水溫的關(guān)系為:

    當(dāng)取Ts=Tm=0℃,將式(19)代入式(23)可得溫度場下游x處冰蓋厚度隨時間變化的計算公式:

    3 黑龍江冰下爆破預(yù)防冰壩的實踐

    3.1 基本情況黑龍江位于我國最北端,為中俄界河,全長4363 km。黑龍江氣候寒冷,最低氣溫可達-50℃以下,封江日期長,從每年的11初到次年的4月末,冰蓋厚度超過1 m。開江時,通常是上游先開,下游后開,形成“武開江”。自1950年代以來,黑龍江局部江段卡塞幾乎年年發(fā)生,平均3年左右形成一次具有一定規(guī)模的冰壩洪水,嚴重威脅到沿岸人民群眾生命財產(chǎn)安全。

    “黑龍江冰情預(yù)報及災(zāi)害防治研究”課題組與黑龍江防汛抗旱保障中心和黑龍江省水文局合作,從2015—2019年在黑龍江開展了冰下爆破預(yù)防冰壩的理論探索及實踐,采用的是巖石乳化炸藥。每年,在黑龍江上游漠河縣、呼瑪縣和黑河市容易發(fā)生卡冰的22 處冰壩易發(fā)點實施冰下爆破。初期,參考國內(nèi)外的實踐經(jīng)驗,破冰寬度一般為天然河道冬季河寬的60%左右,因為兩岸附近40%左右的冰蓋,其有效水深已經(jīng)很小,甚至形成連底冰,在這些地方實施爆破的效果一般不好[13]。后期,隨著把冰、水、熱動力學(xué)理論與爆破效果相結(jié)合的研究,采用了2.4節(jié)爆破孔的布置,減小河道橫向爆破寬度,增加了縱向爆破的長度,以便形成一定規(guī)模的溫度場。下面分析2018年4月在呼瑪縣金山鄉(xiāng)新街基王八大漢江段淺灘交界處的冰下爆破及其效果。

    為了提高冰下爆破預(yù)防冰壩的效果,首先采用雙頻雷達[20]測量了爆破實施江段的冰蓋厚度和水深的分布,以確定深泓線大致位置。如圖4所示,所研究淺灘上游小島將江道分成兩部分,順流向右側(cè)部分屬于我國,小島及左側(cè)部分屬于俄國。交匯口上游水淺流速較快,主槽(有水區(qū)域)過水?dāng)嗝嫫骄罴s1 m,平均冰蓋厚度1.1 m,漫灘(無水區(qū))鋪滿冰蓋。交匯口下游淺灘主槽過水?dāng)嗝嫫骄罴s1.5 m,流速較小,平均冰蓋厚度1.5 m,如圖5所示。由于黑龍江冰厚較大,融化后會使得槽蓄量大大增加,對開江洪水的影響不容忽視。

    圖4 淺灘上游典型斷面實測冰面、冰底、江床高程

    圖5 淺灘中部縱向和橫向典型斷面實測冰面、冰底、江床高程

    3.2 炸藥量及爆破孔布置為了研究需要,把長約2300 m的江段按照從上游到下游順序劃分為6個區(qū),爆破孔布置如下:

    1區(qū)(交匯口上游):長500 m×寬100 m,布置2列,孔距50 m,共20孔,每列每孔炸藥量分別是24、27 kg;

    2區(qū):長400 m×寬80 m,布置2列,孔距40 m,共20孔,每列每孔炸藥量分別是24、27 kg。

    3區(qū):長400 m×寬80 m,布置2列,孔距40 m,共20孔,每列每孔炸藥量分別是21、27 kg。

    4區(qū):長400 m×寬80 m,布置2列,孔距40 m,共20孔,每列每孔炸藥量分別是21、27 kg;

    5區(qū):長300 m×寬90 m,布置3列,共40孔,每列每孔炸藥量分別是3、9、18 kg,分別對應(yīng)孔距15、30、30 m;

    6區(qū):長度300 m×寬90 m,布置3列,共40孔,每列每孔炸藥量分別是6、9、18 kg,分別對應(yīng)孔距15、30、30 m。

    呼瑪合計炸藥用量2660 kg,爆破孔160個。在下面統(tǒng)計分析的過程中,把冰蓋厚度變化在0.1 m范圍的相同炸藥量的爆破坑歸為同一類,例如:爆破孔裝藥量為3 kg,冰蓋厚度在1.0~1.1 m之間,則歸為同一類;爆破孔裝藥量為3 kg,冰蓋厚度在1.1~1.2 m之間,則歸為另一類;以此類推,爆破孔裝藥量為6 kg,冰蓋厚度在1.0~1.1 m之間,也歸為同一類等。為了減小觀測誤差和冰厚不均的影響,對同一類爆破坑的冰厚、水深、爆破坑直徑分別取平均值,則可得表1爆破區(qū)平均冰厚h、平均水深H、炸藥量W、爆破坑平均直徑d的一覽表,其中炸藥量24 kg爆破坑冰厚變化從0.3 m到1.2 m,爆破坑直徑變化較大,表中只列出冰厚超過1 m的統(tǒng)計結(jié)果。

    觀察表1可得下述結(jié)論:

    (1)爆破坑平均直徑隨炸藥量的增加而增加,但在炸藥量超過18 kg時,爆破坑平均直徑變化不大,甚至減小。例如,在炸藥量18 kg、冰厚1.16~1.44 m 時,爆破坑平均直徑在12.6~13.32 m 之間,但是,在炸藥量分別為21、24、27 kg、冰厚超過1 m 時,爆破坑最大平均直徑分別為12.5、11.7、11.86 m。一般來說,炸藥量增大到一定程度,爆破坑直徑可能增大微小,但爆破坑直徑隨炸藥量增加減小屬于異常,產(chǎn)生這一原因可能是21~27 kg炸藥量爆破孔下面過水?dāng)嗝嫠钇?,不能充分發(fā)揮炸藥的爆破的能力,今后冰下爆破應(yīng)盡可能在水深處實施。

    表1 爆破直接效果一覽表

    (2)在炸藥量小于等于18 kg且冰厚在1~1.5 m時,同一炸藥量產(chǎn)生的爆破坑平均直徑隨冰厚的變化不大。例如,當(dāng)炸藥量為6 kg 時,在冰厚1.02~1.34 m 范圍內(nèi),最小與最大爆破坑直徑分別為6.82 m 和7.32 m,僅相差0.5 m。利用空氣、冰與水物理特性差異實現(xiàn)水情全天候自動化監(jiān)測[21-22]表明,在實施冰下爆破前,氣溫表現(xiàn)為白天正溫,夜晚負溫,冰蓋內(nèi)部出現(xiàn)液態(tài)水,冰蓋層疏松[23],這可能是產(chǎn)生前述現(xiàn)象的主要原因。

    從國外經(jīng)驗式(1)和(2)來看,評定冰下爆破的直接效果是炸藥量和爆破坑直徑與冰厚的最優(yōu)關(guān)系。當(dāng)冰厚1.1 m 和1.45 m 時,如果采用式(2)確定炸藥量Wopt和爆破坑直徑D,當(dāng)冰厚1.1 m 時,Wopt≈28 kg 和D=16.5 m;當(dāng)冰厚1.45 m 時Wopt=63.6 kg 和D=21.7 m,但是,在呼瑪縣江段的爆破實測,在冰厚1.07 m和炸藥量為27 kg時,D=11.9 m,比式(2)計算的D=16.5 m小得多,其爆破直接效果不如炸藥量18 kg的。其主要原因可能是式(1)和(2)來源的爆破試驗是氣溫長期在負溫下進行的,那時冰蓋堅實,結(jié)構(gòu)強度高,而冰下爆破預(yù)防冰壩是在冰蓋疏松條件下進行的。可以得出重要結(jié)論6:國外經(jīng)驗公式(1)和(2)炸藥量和爆破坑直徑與冰厚的最優(yōu)關(guān)系不適用于黑龍江開江期冰下爆破實際,當(dāng)冰厚超過1 m時,將造成炸藥量的很大浪費。

    根據(jù)表1,在冰厚超過1 m 時可得圖6 爆破坑直徑與炸藥量關(guān)系圖,其中曲線為多項式趨勢線,爆破坑直徑與炸藥量的函數(shù)關(guān)系近似為:

    圖6 爆破坑直徑與炸藥量關(guān)系

    為了評估爆破的直接效果,表1也列出了單位炸藥量的破冰量V/W,其中:為炸藥量W破碎的冰蓋體積。V/W大,表示1kg 炸藥量破碎的冰蓋體積大,爆破的直接效果好。觀察表1 可知,當(dāng)W=3 kg時,V/W=(12.1~19.1)m3/kg,爆破的直接效果最好;當(dāng)W=9 kg時,V/W=(11.6~14.1)m3/kg,爆破的直接效果次之;當(dāng)W=18kg時,V/W=(8.1~11.1)m3/kg;當(dāng)W>18kg時,V/W<7.0 m3/kg,爆破的直接效果較差。

    綜上所述,考慮實施冰下爆破前冰蓋層疏松的特點,建議黑龍江冰下爆破預(yù)防冰壩的單孔炸藥量不超過9 kg。

    3.3 爆破后冰情的變化呼瑪縣天氣預(yù)報日平均溫度將在2018年4月7日轉(zhuǎn)正,所以冰下爆破于這天實施。表2 列出了呼瑪縣2018年4月實際歷史天氣,從中可知呼瑪縣日平均氣溫直到2018年4月14日才穩(wěn)定轉(zhuǎn)正。需要說明的是,黑龍江爆破江段的氣溫通常比呼瑪縣氣溫低,另外,由于缺乏日平均氣溫資料,采用了以最高氣溫和最低氣溫平均的方法估計。

    2018年4月爆破江段的冰情變化可劃分為四個階段:

    第一階段:7日至13日期間,日平均氣溫在-2.5℃與5.5℃之間交替變化,當(dāng)夜晚氣溫在零度以下時,爆破坑水面會結(jié)冰;當(dāng)白天氣溫轉(zhuǎn)正后,冰蓋表面融化形成水膜。

    第二階段:14日至17日期間,日平均氣溫從2.5℃增加到11℃,自17日開始爆破江段1 區(qū)和2區(qū)爆破區(qū)域冰蓋融化形成明水清溝,如圖7所示;但是,其他爆破區(qū)域爆破坑沒有出現(xiàn)明顯水面。產(chǎn)生這一現(xiàn)象的主要原因是:(1)1區(qū)和2區(qū)的冰蓋較薄,爆破區(qū)域冰厚在0.3~1.1 m之間,雖然爆破孔間距較大,達40~50 m,但是由于冰下過水?dāng)嗝嫫骄罴s為1 m,流速較高,爆破坑內(nèi)碎冰同時受到水溫融化和底部流速的沖蝕,另外,由2.3節(jié)結(jié)論1和2:當(dāng)水溫為正時,冰蓋厚度和爆破坑側(cè)壁冰融化的速率隨流速的增加而增加,所以1、2區(qū)爆破坑融化速度快;(2)其他爆破區(qū)域冰蓋較厚,流速較小,爆破坑融化較慢。

    第三階段:18日至23日期間,19日爆破江段下雪,延緩了爆破區(qū)域的解凍,經(jīng)歷20—23日的持續(xù)正溫,特別是23日平均氣溫達12.5℃,到24日,除江中暗礁外,爆破江段完全解凍,形成明流,見圖8,但爆破江段上下游仍然處于冰封狀態(tài)。

    第四階段:25日至28開江期間,凌情總體平穩(wěn),僅黑龍江上游漠河縣洛古河江段和呼瑪縣鷗浦江段由于河道水位偏低,河道中出現(xiàn)淺灘,導(dǎo)致水流動力不足,流凌不暢出現(xiàn)卡塞并形成冰壩。洛古河卡塞長度30多公里,歷時5天后逐漸解體;鷗浦江段卡塞當(dāng)日內(nèi)解除,未造成險情和災(zāi)情。

    綜上所述,在冰蓋較厚而流速較小時,考慮到初期爆破坑直徑一般小于15 m:(1)應(yīng)適當(dāng)減小爆破孔間距;(2)在江面較寬且江中存在暗礁的情況下,應(yīng)適當(dāng)增加爆破坑縱向列數(shù);(3)爆破孔應(yīng)布置在水深的地方,以使得爆破江段盡早解凍,形成溫度場相應(yīng),促使下游更大范圍提前解凍,減小開江過程中發(fā)生冰凌卡塞的風(fēng)險。

    3.4 溫度場的效果下面應(yīng)用式(17)、式(22)、式(24)計算分析爆破區(qū)形成明流后溫度場出口水溫TW1、溫度場對下游冰蓋下水溫的影響長度L及冰蓋厚度隨氣溫、熱交換系數(shù)、溫度場縱向長度的變化。

    表2 黑龍江呼瑪縣2018年4月天氣

    圖7 2018年4月17日1區(qū)和2區(qū)形成清溝

    圖8 2018年4月24日爆破江段冰蓋融化

    計算時假設(shè)江道為矩形,大氣和水的熱交換系數(shù)hWa=20 W/(m2·℃),大氣與冰蓋的熱交換系數(shù)hai=20W/(m2·℃);冰面溫度Ts=0.0℃,冰蓋底溫度Tm=0.0℃,溫度場進口水溫TW0=0.0℃;水的密度ρ=1000 kg/m3,冰的密度ρi=917 kg/m3,水的比熱Cp=4217.7 J/kg℃,冰的潛熱Li=3.33×105J/kg,溫度場(爆破區(qū))沿流向的長度x1=2300 m。

    由于缺乏實測的流速數(shù)據(jù),采用兩個典型情況確定溫度場流速與水深的關(guān)系:

    3.4.1hWi和L與流速和水深的關(guān)系 爆破前,溫度場進口冰下平均水深約為1 m 和平均冰蓋厚1.1 m,溫度場出口平均水深約為1.5 m和平均冰蓋厚1.5 m。爆破區(qū)冰蓋融化形成溫度場后,原來的冰域轉(zhuǎn)變?yōu)樗?,槽蓄量增加,由于下游冰蓋尚未解凍,在上下游冰蓋下流量不變的情況下,溫度場出口水深將由原來冰下1.5 m增加到,溫度場進口水深小于,由于溫度場水深增加,流速將下降。

    假設(shè)溫度場下游冰蓋下單寬流量與溫度場相同,即U2H2=U1H1時,可得表3列出的冰和水熱交換系數(shù)hWi與溫度場對下游冰蓋下水溫的影響長度L的值。

    觀察表3可得下述結(jié)論:

    (1)hWi隨流速U2的增加而增加。例如:對情況1,當(dāng)U2由0.13 m/s增加到0.65 m/s時,hWi由326 W/(m2·℃)增加到1218 W/(m2·℃);對情況2,當(dāng)U2由0.19 m/s 增加到0.98 m/s 時,hWi由451 W/(m2·℃)增加到1685 W/(m2·℃)。

    (2)因為hWi>>hWa和hai,所以微小的水溫升高都會加快冰蓋的融化。因hWa和hWa約為20W/(m2·℃),對情況1,hWi>(16~60)hWa;對情況2,hWi>(22~80)hWa。

    (3)溫度場對下游冰蓋下水溫的影響長度L隨水深和流速的增加而增加,其影響長度接近或者大于溫度場長度。

    表3 hWi和L與溫度場及下游冰蓋下流速和水深的對應(yīng)關(guān)系

    3.4.2 溫度場出口水溫和下游冰蓋前緣冰厚與氣溫和水深的關(guān)系 在溫度場長度x1=2300m 條件下,對于情況1和2,在1天時間T=24×3600=84400 s時,可得圖9和10溫度場出口水溫TW1、下游冰蓋前緣冰厚改變量隨氣溫和水深的變化,圖中Dh1為下游冰蓋前緣底部在水溫TW1時受冰和水熱交換影響減小的厚度(沒有考慮冰蓋迎水面融化的影響),Dh2為下游冰蓋前緣頂部受氣溫和水熱交換影響減小的厚度,H=H1。

    觀察圖9和10可得下述結(jié)論:(1)溫度場出口水溫隨水深的增加而減小,隨氣溫的增加而增大,隨流速的增加而減?。唬?)溫度場融化下游冰蓋底部的能力Dh1與氣溫融化冰蓋頂部的能力Dh2處于同一數(shù)量級,例如:在TW1=0.1~0.5℃、Ta=10℃和H=H2=1~3 m時,在1天時間內(nèi)下游冰蓋前緣底部融化量Dh1=0.035~0.052 m,而冰蓋前緣頂部氣溫融化量Dh2=0.055 m;(3)在1 天之內(nèi)影響長度L融化的冰量=(0.5Dh1+Dh2)L。

    圖9 情況1時溫度場出口水溫、下游冰蓋前緣冰厚隨氣溫和水深的變化

    圖10 情況1時溫度場出口水溫、下游冰蓋前緣冰厚隨氣溫和水深的變化

    4 結(jié)論

    把冰、水、熱動力學(xué)理論與爆破效果相結(jié)合,進行了冰蓋下爆破預(yù)防冰壩的理論探索。建立了描述爆破坑融化擴大與氣溫、水溫、流速關(guān)系的數(shù)學(xué)模型,證明了冰蓋底部和側(cè)壁冰融化的速率隨流速的增加而增加,以及爆破孔布置的最佳位置是河道冰壩易發(fā)點的主槽深泓線附件。冰下爆破預(yù)防冰壩不僅要考慮爆破的直接效果,更要考慮冰、水、熱動力學(xué)作用的間接效果,提出了在深泓線附近沿流向布設(shè)1~3列爆破孔形成一定規(guī)模溫度場的方法,以便利用主槽和漫灘的質(zhì)量和熱量對流使漫灘和下游冰蓋加速融化。然后,理論上證明溫度場出口水溫與溫度場長度、氣溫及水與空氣熱交換系數(shù)的積成正比,與單寬流量成反比;冰蓋下水深和流速越大,場水溫對下游冰蓋下水溫的影響長度越大;并提出了冰厚隨時間變化的解析式。最后,應(yīng)用理論成果指導(dǎo)在黑龍江冰下爆破預(yù)防冰壩的實踐,結(jié)果證明了理論的正確性,并得到了爆破坑直徑與炸藥量關(guān)系的經(jīng)驗公式。實踐結(jié)果也表明,相同炸藥量產(chǎn)生的爆破坑直徑隨冰厚的變化不大,每個爆破孔炸藥量不宜超過9 kg;冰水熱交換系數(shù)hWi的是氣溫和水或者冰面的熱交換系數(shù)hWa的數(shù)十倍,微小的水溫升高都會加快下游冰蓋的融化,并且溫度場融化下游冰蓋底部的能力與氣溫融化冰蓋頂部的能力處于同一數(shù)量級。

    猜你喜歡
    冰蓋水深融化
    軍事文摘(2024年6期)2024-04-30 03:13:59
    書法靜水深流
    河北水利(2022年10期)2022-12-29 11:48:12
    基于水深分段選擇因子的多光譜影像反演水深
    海洋通報(2022年4期)2022-10-10 07:40:32
    格陵蘭島的冰蓋悄悄融化
    參花(下)(2022年1期)2022-01-15 00:45:01
    長距離輸水工程的冰期冰蓋數(shù)值模擬研究
    一起融化
    民族音樂(2016年1期)2016-08-28 20:02:52
    融化的Ice Crean
    健康女性(2014年8期)2015-05-12 03:10:47
    GPS RTK技術(shù)在水深測量中的應(yīng)用
    浸入式水深監(jiān)測儀器的設(shè)計
    冰如何開始融化
    国产淫片久久久久久久久| av视频在线观看入口| 国产亚洲91精品色在线| 天堂√8在线中文| 国产黄a三级三级三级人| 亚洲av成人av| 亚洲熟妇中文字幕五十中出| 床上黄色一级片| 亚洲精品国产成人久久av| 97在线视频观看| 亚洲性夜色夜夜综合| 午夜老司机福利剧场| 亚洲欧美精品综合久久99| 亚洲av中文字字幕乱码综合| 亚洲精品亚洲一区二区| 中文字幕熟女人妻在线| 麻豆av噜噜一区二区三区| 午夜视频国产福利| 午夜福利在线在线| 91久久精品国产一区二区三区| 亚洲国产色片| 在线国产一区二区在线| 12—13女人毛片做爰片一| 亚洲av美国av| 简卡轻食公司| 欧美高清成人免费视频www| 国产高清视频在线播放一区| 91久久精品国产一区二区三区| 国产成人aa在线观看| 日韩大尺度精品在线看网址| 一区二区三区高清视频在线| 亚洲内射少妇av| 久久久色成人| 午夜视频国产福利| 深夜a级毛片| 偷拍熟女少妇极品色| 麻豆乱淫一区二区| 午夜免费男女啪啪视频观看 | 狠狠狠狠99中文字幕| 桃色一区二区三区在线观看| 日韩,欧美,国产一区二区三区 | 成人av一区二区三区在线看| 国产精品日韩av在线免费观看| 最近在线观看免费完整版| 国产又黄又爽又无遮挡在线| 久久久久久九九精品二区国产| 免费高清视频大片| 久久精品国产亚洲网站| 国内精品一区二区在线观看| 香蕉av资源在线| 嫩草影院入口| 亚洲人成网站高清观看| 22中文网久久字幕| 噜噜噜噜噜久久久久久91| 美女 人体艺术 gogo| 午夜视频国产福利| 五月伊人婷婷丁香| 国产精品久久视频播放| 久久久精品大字幕| 男人舔奶头视频| 91久久精品国产一区二区三区| 一级黄片播放器| 亚洲中文日韩欧美视频| av国产免费在线观看| 国产高清有码在线观看视频| 久久久久精品国产欧美久久久| 亚洲人成网站在线播放欧美日韩| 亚洲国产欧美人成| 一级av片app| 一个人看视频在线观看www免费| 亚洲欧美日韩无卡精品| 91久久精品国产一区二区成人| 中文字幕免费在线视频6| 黄色配什么色好看| 99久久中文字幕三级久久日本| 久久久久九九精品影院| 久久久久久久亚洲中文字幕| 国产极品精品免费视频能看的| 日韩三级伦理在线观看| 国产成人aa在线观看| 亚洲精华国产精华液的使用体验 | 我要看日韩黄色一级片| 欧美成人精品欧美一级黄| 欧美成人精品欧美一级黄| 可以在线观看毛片的网站| 国内少妇人妻偷人精品xxx网站| 亚洲内射少妇av| 狂野欧美白嫩少妇大欣赏| 赤兔流量卡办理| 永久网站在线| 丰满乱子伦码专区| 亚洲精品456在线播放app| 国产伦在线观看视频一区| 国产真实伦视频高清在线观看| 两个人视频免费观看高清| 免费看日本二区| 欧美绝顶高潮抽搐喷水| 日日撸夜夜添| 可以在线观看的亚洲视频| 真人做人爱边吃奶动态| 97人妻精品一区二区三区麻豆| av在线播放精品| 欧美国产日韩亚洲一区| av在线蜜桃| 国产成人福利小说| 免费看av在线观看网站| 小说图片视频综合网站| 九九热线精品视视频播放| 97超级碰碰碰精品色视频在线观看| 国产伦在线观看视频一区| 少妇被粗大猛烈的视频| 亚洲最大成人手机在线| 国产在线男女| 3wmmmm亚洲av在线观看| 国产亚洲精品久久久com| 国产欧美日韩精品一区二区| 久久国产乱子免费精品| 久久午夜亚洲精品久久| 亚洲无线观看免费| 日韩欧美一区二区三区在线观看| 一进一出抽搐动态| 亚洲性久久影院| 婷婷色综合大香蕉| 亚洲av免费高清在线观看| 麻豆一二三区av精品| av在线观看视频网站免费| 天美传媒精品一区二区| 十八禁国产超污无遮挡网站| 自拍偷自拍亚洲精品老妇| 搡老妇女老女人老熟妇| 真实男女啪啪啪动态图| 日本与韩国留学比较| av在线播放精品| 亚洲天堂国产精品一区在线| 亚洲真实伦在线观看| 亚洲七黄色美女视频| 午夜精品国产一区二区电影 | 国产三级中文精品| 午夜福利在线观看吧| 免费观看精品视频网站| 亚洲成人中文字幕在线播放| 国产精品嫩草影院av在线观看| 亚洲中文日韩欧美视频| 丰满的人妻完整版| 日本撒尿小便嘘嘘汇集6| 免费看光身美女| 久久99热这里只有精品18| 1024手机看黄色片| 少妇的逼水好多| 不卡一级毛片| 日韩精品有码人妻一区| 我的女老师完整版在线观看| 久久久久久久久大av| 亚洲内射少妇av| 午夜亚洲福利在线播放| 男插女下体视频免费在线播放| 桃色一区二区三区在线观看| 国产乱人偷精品视频| 成人永久免费在线观看视频| 身体一侧抽搐| 国产精品av视频在线免费观看| 嫩草影院精品99| 日本欧美国产在线视频| 午夜福利成人在线免费观看| av在线老鸭窝| 中文资源天堂在线| 99国产极品粉嫩在线观看| 国产一区二区亚洲精品在线观看| 美女高潮的动态| 人人妻人人澡人人爽人人夜夜 | 国内久久婷婷六月综合欲色啪| 永久网站在线| 男女之事视频高清在线观看| 日韩大尺度精品在线看网址| 五月伊人婷婷丁香| 51国产日韩欧美| 欧美绝顶高潮抽搐喷水| 亚洲精品粉嫩美女一区| 尤物成人国产欧美一区二区三区| 午夜影院日韩av| 国产精品久久视频播放| 亚洲,欧美,日韩| 日韩欧美一区二区三区在线观看| 美女高潮的动态| 九九在线视频观看精品| 国产精品女同一区二区软件| 人妻丰满熟妇av一区二区三区| 亚洲美女黄片视频| 好男人在线观看高清免费视频| 搞女人的毛片| 亚洲精品粉嫩美女一区| 午夜久久久久精精品| 日韩一本色道免费dvd| av在线播放精品| 青春草视频在线免费观看| 无遮挡黄片免费观看| 一级毛片我不卡| 看黄色毛片网站| 菩萨蛮人人尽说江南好唐韦庄 | 中文字幕熟女人妻在线| 欧美一级a爱片免费观看看| 久久人人精品亚洲av| 午夜老司机福利剧场| 亚洲精华国产精华液的使用体验 | 亚洲精品国产av成人精品 | 午夜老司机福利剧场| 亚洲精华国产精华液的使用体验 | 男人的好看免费观看在线视频| 小蜜桃在线观看免费完整版高清| 国产精品国产三级国产av玫瑰| 韩国av在线不卡| 国产老妇女一区| 久久人人爽人人爽人人片va| 国产一区二区三区在线臀色熟女| 免费观看在线日韩| 波野结衣二区三区在线| 成人国产麻豆网| 国产成人a∨麻豆精品| 成人鲁丝片一二三区免费| 男女那种视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区在线观看日韩| 我要搜黄色片| 免费看日本二区| 亚洲欧美日韩东京热| 亚洲精品成人久久久久久| 日日干狠狠操夜夜爽| 国产精品久久久久久久久免| 欧美+日韩+精品| 国产 一区精品| 亚洲无线在线观看| 日韩 亚洲 欧美在线| 成人永久免费在线观看视频| 久久久精品94久久精品| 午夜精品一区二区三区免费看| 五月伊人婷婷丁香| 成年免费大片在线观看| 成年女人看的毛片在线观看| 亚洲欧美日韩卡通动漫| 国产伦一二天堂av在线观看| 国产黄a三级三级三级人| 亚洲国产精品久久男人天堂| 国产午夜精品久久久久久一区二区三区 | 两个人视频免费观看高清| 成人欧美大片| 日韩欧美三级三区| 精品久久久久久成人av| 亚洲精品日韩av片在线观看| 99热6这里只有精品| 国产亚洲精品久久久久久毛片| 超碰av人人做人人爽久久| 天堂网av新在线| 日韩精品有码人妻一区| 老司机午夜福利在线观看视频| 久久精品91蜜桃| 99国产精品一区二区蜜桃av| 国产伦在线观看视频一区| 国产国拍精品亚洲av在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产精品福利在线免费观看| 热99re8久久精品国产| av福利片在线观看| 免费搜索国产男女视频| 久久久久国产网址| 欧美一区二区亚洲| 18+在线观看网站| 国产精品一区二区三区四区久久| 乱码一卡2卡4卡精品| 五月伊人婷婷丁香| 亚洲精品456在线播放app| 日韩一本色道免费dvd| 春色校园在线视频观看| 亚洲专区国产一区二区| 欧美成人a在线观看| 欧美日韩国产亚洲二区| 舔av片在线| 精品一区二区三区av网在线观看| 丰满的人妻完整版| 成熟少妇高潮喷水视频| 人人妻,人人澡人人爽秒播| 免费观看的影片在线观看| 亚洲av五月六月丁香网| 日本a在线网址| 成人精品一区二区免费| 久久久久久久久久黄片| 又爽又黄a免费视频| 精品久久久久久久久亚洲| 国产中年淑女户外野战色| 日本-黄色视频高清免费观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产日韩欧美精品在线观看| 国产精品,欧美在线| 亚洲av电影不卡..在线观看| 国产 一区精品| 精品无人区乱码1区二区| 欧美三级亚洲精品| 午夜老司机福利剧场| 又黄又爽又免费观看的视频| 一个人免费在线观看电影| 人人妻人人看人人澡| 亚洲18禁久久av| 日本在线视频免费播放| 午夜a级毛片| 亚洲精品456在线播放app| 波多野结衣高清无吗| 好男人在线观看高清免费视频| 精品不卡国产一区二区三区| 18禁在线无遮挡免费观看视频 | 插阴视频在线观看视频| 一级黄色大片毛片| 午夜老司机福利剧场| 看十八女毛片水多多多| 丝袜美腿在线中文| 久久久久国产精品人妻aⅴ院| 免费看美女性在线毛片视频| 搡老熟女国产l中国老女人| 97在线视频观看| 午夜精品一区二区三区免费看| 99热网站在线观看| 国产乱人偷精品视频| 一区福利在线观看| 成年av动漫网址| 国产真实伦视频高清在线观看| 久久精品国产亚洲网站| 自拍偷自拍亚洲精品老妇| 在线免费观看不下载黄p国产| 我的老师免费观看完整版| 国产精华一区二区三区| 欧美日韩综合久久久久久| 色播亚洲综合网| 精品少妇黑人巨大在线播放 | 色综合站精品国产| 真人做人爱边吃奶动态| 露出奶头的视频| 免费一级毛片在线播放高清视频| 波多野结衣高清作品| 又黄又爽又刺激的免费视频.| 九九久久精品国产亚洲av麻豆| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品日韩av片在线观看| 老司机午夜福利在线观看视频| 亚洲美女黄片视频| 91精品国产九色| 亚洲乱码一区二区免费版| 中文字幕精品亚洲无线码一区| av在线蜜桃| 久久午夜亚洲精品久久| 日韩欧美 国产精品| 亚洲欧美日韩高清专用| 91久久精品电影网| 狂野欧美白嫩少妇大欣赏| 日本免费一区二区三区高清不卡| 欧美xxxx性猛交bbbb| 日本-黄色视频高清免费观看| 久久久久九九精品影院| 在线天堂最新版资源| 一区二区三区免费毛片| 亚洲电影在线观看av| 日本一二三区视频观看| 永久网站在线| 熟女人妻精品中文字幕| 日日啪夜夜撸| 男人狂女人下面高潮的视频| 一a级毛片在线观看| 晚上一个人看的免费电影| 精华霜和精华液先用哪个| 如何舔出高潮| 久久婷婷人人爽人人干人人爱| 亚洲精品乱码久久久v下载方式| 亚洲欧美成人综合另类久久久 | 国产v大片淫在线免费观看| 亚洲欧美日韩无卡精品| 精品久久久久久成人av| 高清日韩中文字幕在线| 成人综合一区亚洲| 亚洲精品国产成人久久av| 18+在线观看网站| 欧美性猛交╳xxx乱大交人| 亚洲精华国产精华液的使用体验 | 成人鲁丝片一二三区免费| 日韩欧美精品v在线| 在现免费观看毛片| 久久鲁丝午夜福利片| 国产精品一区二区三区四区久久| 国产精品伦人一区二区| 一个人看的www免费观看视频| 黑人高潮一二区| 免费看av在线观看网站| 国产免费一级a男人的天堂| 一区二区三区四区激情视频 | 热99在线观看视频| 国产人妻一区二区三区在| 伦理电影大哥的女人| 国产精品日韩av在线免费观看| 乱码一卡2卡4卡精品| 色在线成人网| 久久午夜亚洲精品久久| 国内精品一区二区在线观看| АⅤ资源中文在线天堂| 久久久久久久久久久丰满| 久久久久久伊人网av| 精品人妻一区二区三区麻豆 | 中出人妻视频一区二区| 一区福利在线观看| 国产色爽女视频免费观看| 天堂动漫精品| 18+在线观看网站| 免费大片18禁| h日本视频在线播放| 性色avwww在线观看| 亚洲中文字幕一区二区三区有码在线看| 亚洲四区av| 1000部很黄的大片| 麻豆乱淫一区二区| 国产一区二区在线观看日韩| 中文在线观看免费www的网站| 国产精品免费一区二区三区在线| 亚洲在线观看片| 哪里可以看免费的av片| 亚洲国产日韩欧美精品在线观看| 欧美zozozo另类| 国产黄片美女视频| 久久午夜亚洲精品久久| 嫩草影视91久久| 女人被狂操c到高潮| 人人妻人人澡欧美一区二区| 特大巨黑吊av在线直播| 色吧在线观看| 99久久精品一区二区三区| 欧美性猛交╳xxx乱大交人| 日本熟妇午夜| 国产成人影院久久av| 国产成人精品久久久久久| 简卡轻食公司| 深夜a级毛片| 亚洲精品乱码久久久v下载方式| 99久久中文字幕三级久久日本| 亚洲三级黄色毛片| 偷拍熟女少妇极品色| 国产精品国产高清国产av| 一边摸一边抽搐一进一小说| 老熟妇仑乱视频hdxx| 99热这里只有是精品在线观看| eeuss影院久久| 午夜久久久久精精品| 亚洲熟妇熟女久久| 欧美日韩国产亚洲二区| 亚洲高清免费不卡视频| 色哟哟哟哟哟哟| 免费在线观看成人毛片| 久久精品久久久久久噜噜老黄 | 成熟少妇高潮喷水视频| 亚洲一级一片aⅴ在线观看| av专区在线播放| 99九九线精品视频在线观看视频| 美女黄网站色视频| 免费人成在线观看视频色| 听说在线观看完整版免费高清| 国内精品一区二区在线观看| 毛片女人毛片| 91午夜精品亚洲一区二区三区| 日本黄色片子视频| 丰满的人妻完整版| 亚洲av中文字字幕乱码综合| 国产高清有码在线观看视频| 亚洲国产高清在线一区二区三| 亚洲人与动物交配视频| eeuss影院久久| 男女那种视频在线观看| 舔av片在线| 白带黄色成豆腐渣| 男人和女人高潮做爰伦理| 婷婷精品国产亚洲av在线| 成年女人毛片免费观看观看9| 一进一出好大好爽视频| 赤兔流量卡办理| 岛国在线免费视频观看| av.在线天堂| 99热6这里只有精品| 少妇丰满av| 久久久精品大字幕| 免费观看人在逋| 波野结衣二区三区在线| 国产成人aa在线观看| 亚洲精品成人久久久久久| 国产又黄又爽又无遮挡在线| 噜噜噜噜噜久久久久久91| 成人高潮视频无遮挡免费网站| 亚洲欧美清纯卡通| 久久欧美精品欧美久久欧美| 亚洲图色成人| 女生性感内裤真人,穿戴方法视频| 国产成人freesex在线 | 亚洲美女黄片视频| 老师上课跳d突然被开到最大视频| 一级毛片电影观看 | 精品午夜福利视频在线观看一区| 久久久久久久久久黄片| av在线亚洲专区| 美女内射精品一级片tv| 又粗又爽又猛毛片免费看| 国产v大片淫在线免费观看| 日韩欧美三级三区| 国产伦在线观看视频一区| 久久久久国内视频| 日产精品乱码卡一卡2卡三| 别揉我奶头 嗯啊视频| 男女边吃奶边做爰视频| 国产精品一区二区三区四区免费观看 | 久久亚洲国产成人精品v| 亚洲高清免费不卡视频| 一进一出抽搐gif免费好疼| 精品午夜福利视频在线观看一区| 亚洲av成人av| 久久久久国产精品人妻aⅴ院| 在现免费观看毛片| 亚洲成人中文字幕在线播放| 亚洲欧美日韩高清在线视频| 午夜亚洲福利在线播放| 亚洲中文字幕一区二区三区有码在线看| 日本黄色片子视频| 日韩欧美三级三区| 搞女人的毛片| 国产精品嫩草影院av在线观看| 深夜a级毛片| 成人美女网站在线观看视频| 久久综合国产亚洲精品| 日本黄色视频三级网站网址| 桃色一区二区三区在线观看| 在线免费观看的www视频| 97超碰精品成人国产| 亚洲精品久久国产高清桃花| 国产精品久久久久久亚洲av鲁大| 欧美国产日韩亚洲一区| 噜噜噜噜噜久久久久久91| 亚洲成人av在线免费| 免费在线观看成人毛片| 免费电影在线观看免费观看| av在线亚洲专区| 免费观看的影片在线观看| 欧美日韩在线观看h| 精品一区二区三区av网在线观看| 日本爱情动作片www.在线观看 | 国产片特级美女逼逼视频| 国产午夜精品久久久久久一区二区三区 | 国产视频一区二区在线看| 中文字幕精品亚洲无线码一区| 久久精品夜夜夜夜夜久久蜜豆| 麻豆一二三区av精品| 久久天躁狠狠躁夜夜2o2o| 亚洲激情五月婷婷啪啪| 国产在视频线在精品| 偷拍熟女少妇极品色| 午夜福利在线观看免费完整高清在 | 九九爱精品视频在线观看| 欧美一区二区国产精品久久精品| 成年版毛片免费区| 精品无人区乱码1区二区| 国产私拍福利视频在线观看| 国产精品亚洲美女久久久| 久久午夜亚洲精品久久| 国产不卡一卡二| 成人综合一区亚洲| 久久韩国三级中文字幕| 国产乱人视频| 成熟少妇高潮喷水视频| 成人特级av手机在线观看| 免费av不卡在线播放| 少妇熟女欧美另类| 我的女老师完整版在线观看| 欧美绝顶高潮抽搐喷水| 亚洲欧美中文字幕日韩二区| 日本 av在线| 久久亚洲国产成人精品v| 99久久精品国产国产毛片| 成人三级黄色视频| 午夜精品一区二区三区免费看| 亚洲av第一区精品v没综合| 人妻夜夜爽99麻豆av| 成熟少妇高潮喷水视频| 91久久精品国产一区二区三区| 国产高清不卡午夜福利| 又爽又黄无遮挡网站| 精品久久久久久久久av| 在线a可以看的网站| 欧美色视频一区免费| 18禁黄网站禁片免费观看直播| 联通29元200g的流量卡| 毛片一级片免费看久久久久| 亚洲精品一区av在线观看| 日本一本二区三区精品| 欧美在线一区亚洲| 日产精品乱码卡一卡2卡三| 在线播放国产精品三级| 哪里可以看免费的av片| 欧美高清成人免费视频www| 日本a在线网址| 婷婷精品国产亚洲av| 免费观看的影片在线观看| 久久久久久九九精品二区国产| 国产亚洲欧美98| 草草在线视频免费看| 成人综合一区亚洲| 一级毛片久久久久久久久女| .国产精品久久| 午夜爱爱视频在线播放| 日韩欧美精品v在线| 精品人妻一区二区三区麻豆 | 两个人视频免费观看高清| 亚洲人成网站在线播放欧美日韩| 欧美xxxx性猛交bbbb| 色噜噜av男人的天堂激情| 亚洲欧美日韩高清在线视频| 久久综合国产亚洲精品| 精品不卡国产一区二区三区|