• <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è)計
    冰如何開始融化
    www.自偷自拍.com| 又黄又粗又硬又大视频| 久久精品国产亚洲av香蕉五月| 久久精品国产亚洲av香蕉五月| 级片在线观看| 午夜福利,免费看| 丝袜美足系列| 桃红色精品国产亚洲av| 亚洲av成人不卡在线观看播放网| 国产欧美日韩综合在线一区二区| 欧美午夜高清在线| 少妇 在线观看| 成年女人毛片免费观看观看9| 亚洲久久久国产精品| 欧美国产精品va在线观看不卡| 日韩高清综合在线| 天天躁夜夜躁狠狠躁躁| 色哟哟哟哟哟哟| 不卡av一区二区三区| 一级毛片女人18水好多| 咕卡用的链子| 在线观看一区二区三区激情| 色婷婷av一区二区三区视频| 精品国产美女av久久久久小说| 可以在线观看毛片的网站| 黄片大片在线免费观看| 亚洲中文av在线| 9热在线视频观看99| 多毛熟女@视频| 男女之事视频高清在线观看| avwww免费| 国产一卡二卡三卡精品| 免费高清视频大片| 精品一区二区三卡| 国产亚洲精品一区二区www| 中文字幕色久视频| 少妇裸体淫交视频免费看高清 | 女性生殖器流出的白浆| 麻豆成人av在线观看| 亚洲欧美一区二区三区久久| 欧美一级毛片孕妇| 后天国语完整版免费观看| 亚洲欧美一区二区三区黑人| 国产一区在线观看成人免费| 欧美日本亚洲视频在线播放| 色哟哟哟哟哟哟| 超碰97精品在线观看| 亚洲专区国产一区二区| 国产亚洲精品一区二区www| 欧美精品亚洲一区二区| 真人一进一出gif抽搐免费| 桃色一区二区三区在线观看| 无人区码免费观看不卡| 亚洲男人的天堂狠狠| 午夜成年电影在线免费观看| 9191精品国产免费久久| 99久久久亚洲精品蜜臀av| 亚洲精品成人av观看孕妇| 欧美性长视频在线观看| 美女福利国产在线| 国产国语露脸激情在线看| 在线观看日韩欧美| 欧美中文日本在线观看视频| 国产成人欧美在线观看| 中文字幕色久视频| 女人被狂操c到高潮| 精品第一国产精品| www国产在线视频色| 国产熟女午夜一区二区三区| 国产精华一区二区三区| 露出奶头的视频| 久久久久久久午夜电影 | 午夜福利一区二区在线看| 日本免费一区二区三区高清不卡 | 精品久久蜜臀av无| 母亲3免费完整高清在线观看| 日本精品一区二区三区蜜桃| 少妇被粗大的猛进出69影院| 午夜福利欧美成人| 亚洲精品在线观看二区| 国产无遮挡羞羞视频在线观看| 黄色毛片三级朝国网站| 久久草成人影院| 99久久精品国产亚洲精品| 久久香蕉激情| 国产一区二区激情短视频| 成人免费观看视频高清| 久久欧美精品欧美久久欧美| 欧美一区二区精品小视频在线| 纯流量卡能插随身wifi吗| 国产精品影院久久| 国产xxxxx性猛交| 国产成年人精品一区二区 | 精品日产1卡2卡| 我的亚洲天堂| 一本大道久久a久久精品| 热re99久久国产66热| 欧美一级毛片孕妇| 一级片'在线观看视频| 久久草成人影院| 久久久国产成人精品二区 | 亚洲精品国产区一区二| av欧美777| 黄色视频不卡| 成人免费观看视频高清| 亚洲片人在线观看| 人人澡人人妻人| 日韩欧美免费精品| 午夜亚洲福利在线播放| 亚洲一码二码三码区别大吗| 日韩高清综合在线| 交换朋友夫妻互换小说| 90打野战视频偷拍视频| 免费av毛片视频| 日本 av在线| 久久精品国产清高在天天线| www.999成人在线观看| 欧美激情 高清一区二区三区| 国产精品久久久人人做人人爽| 亚洲国产精品999在线| 国产精品影院久久| 他把我摸到了高潮在线观看| 亚洲精品中文字幕一二三四区| 国产精华一区二区三区| 最新美女视频免费是黄的| 免费日韩欧美在线观看| 每晚都被弄得嗷嗷叫到高潮| 日韩欧美一区二区三区在线观看| 国产一区二区激情短视频| 色精品久久人妻99蜜桃| 首页视频小说图片口味搜索| 久久青草综合色| 欧美精品啪啪一区二区三区| 999精品在线视频| 新久久久久国产一级毛片| 亚洲伊人色综图| 亚洲精品美女久久av网站| 亚洲精华国产精华精| 中出人妻视频一区二区| netflix在线观看网站| 亚洲专区字幕在线| 首页视频小说图片口味搜索| 精品国内亚洲2022精品成人| 国产精品99久久99久久久不卡| 夜夜看夜夜爽夜夜摸 | 在线观看66精品国产| 精品一区二区三卡| 极品教师在线免费播放| 国产极品粉嫩免费观看在线| 人人妻人人添人人爽欧美一区卜| 免费日韩欧美在线观看| 久久精品亚洲熟妇少妇任你| 日韩成人在线观看一区二区三区| 国产高清激情床上av| 欧美日韩亚洲综合一区二区三区_| 交换朋友夫妻互换小说| 美女 人体艺术 gogo| 国产aⅴ精品一区二区三区波| 久久精品国产综合久久久| 黑人巨大精品欧美一区二区mp4| 老熟妇仑乱视频hdxx| 国产精品二区激情视频| 久久九九热精品免费| 国产精品一区二区在线不卡| 在线十欧美十亚洲十日本专区| 日韩欧美国产一区二区入口| 久久精品国产清高在天天线| 精品国产美女av久久久久小说| av在线天堂中文字幕 | 高清黄色对白视频在线免费看| 视频区图区小说| 可以在线观看毛片的网站| 激情视频va一区二区三区| 亚洲三区欧美一区| 精品人妻在线不人妻| www.999成人在线观看| 久久精品91无色码中文字幕| 亚洲第一欧美日韩一区二区三区| 欧美性长视频在线观看| 别揉我奶头~嗯~啊~动态视频| 好男人电影高清在线观看| 搡老熟女国产l中国老女人| 久久国产精品人妻蜜桃| 男女高潮啪啪啪动态图| 青草久久国产| 久久精品91无色码中文字幕| 一边摸一边抽搐一进一小说| 男女下面插进去视频免费观看| 在线观看免费视频网站a站| 亚洲三区欧美一区| av网站在线播放免费| a级片在线免费高清观看视频| 一进一出抽搐gif免费好疼 | 五月开心婷婷网| 1024香蕉在线观看| 五月开心婷婷网| 美女高潮到喷水免费观看| 久久天堂一区二区三区四区| 久久国产精品影院| 亚洲精品av麻豆狂野| 久久香蕉精品热| 国产精品久久久久久人妻精品电影| 国产av在哪里看| 欧美不卡视频在线免费观看 | 色综合站精品国产| 午夜福利欧美成人| 色在线成人网| 精品国产一区二区三区四区第35| 久久香蕉精品热| 一级黄色大片毛片| 国产伦一二天堂av在线观看| 人妻丰满熟妇av一区二区三区| а√天堂www在线а√下载| 免费女性裸体啪啪无遮挡网站| 亚洲自拍偷在线| 丝袜在线中文字幕| 女同久久另类99精品国产91| 在线观看66精品国产| 又黄又爽又免费观看的视频| 一边摸一边抽搐一进一出视频| 脱女人内裤的视频| 久久青草综合色| 欧美黄色淫秽网站| bbb黄色大片| 久久久久久亚洲精品国产蜜桃av| e午夜精品久久久久久久| 免费不卡黄色视频| 亚洲avbb在线观看| 成人18禁在线播放| 亚洲中文日韩欧美视频| 九色亚洲精品在线播放| 亚洲色图av天堂| 纯流量卡能插随身wifi吗| 91字幕亚洲| 亚洲精品在线美女| 999久久久国产精品视频| 久久国产乱子伦精品免费另类| 亚洲国产欧美日韩在线播放| 热re99久久国产66热| 国产高清videossex| 国产欧美日韩一区二区精品| 香蕉国产在线看| 高清毛片免费观看视频网站 | 国产91精品成人一区二区三区| 国产视频一区二区在线看| 好看av亚洲va欧美ⅴa在| 色婷婷久久久亚洲欧美| 精品久久久精品久久久| 美女福利国产在线| 女人被狂操c到高潮| 一级作爱视频免费观看| 午夜精品久久久久久毛片777| 一级片免费观看大全| 久久精品91无色码中文字幕| 精品熟女少妇八av免费久了| 欧美中文日本在线观看视频| 国产在线精品亚洲第一网站| 91九色精品人成在线观看| 欧美大码av| 夜夜夜夜夜久久久久| 超色免费av| 18禁美女被吸乳视频| 亚洲一区二区三区欧美精品| 国产成人精品无人区| 国产成人影院久久av| 香蕉国产在线看| 男人操女人黄网站| 久久人妻福利社区极品人妻图片| 老汉色av国产亚洲站长工具| 黄色女人牲交| 巨乳人妻的诱惑在线观看| 午夜精品国产一区二区电影| 亚洲专区国产一区二区| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品国产一区二区精华液| 视频在线观看一区二区三区| 久久精品人人爽人人爽视色| 国产成年人精品一区二区 | 搡老岳熟女国产| 99久久久亚洲精品蜜臀av| 每晚都被弄得嗷嗷叫到高潮| 乱人伦中国视频| 日韩中文字幕欧美一区二区| 老汉色∧v一级毛片| 新久久久久国产一级毛片| 无遮挡黄片免费观看| 99国产精品免费福利视频| 亚洲久久久国产精品| 老汉色av国产亚洲站长工具| 久久精品国产综合久久久| 91老司机精品| 丝袜美足系列| a级毛片黄视频| 成人三级黄色视频| 亚洲五月天丁香| 久久久久久久久免费视频了| 欧美+亚洲+日韩+国产| 久久中文字幕人妻熟女| 操美女的视频在线观看| 成人精品一区二区免费| 亚洲一区二区三区色噜噜 | 国产欧美日韩精品亚洲av| 国产精品成人在线| 女人被狂操c到高潮| 两性午夜刺激爽爽歪歪视频在线观看 | 极品人妻少妇av视频| 精品午夜福利视频在线观看一区| 日韩有码中文字幕| 国产日韩一区二区三区精品不卡| 久久久久久久久免费视频了| 久久精品91蜜桃| 丰满的人妻完整版| 电影成人av| 欧洲精品卡2卡3卡4卡5卡区| 日韩精品中文字幕看吧| 亚洲国产看品久久| 少妇裸体淫交视频免费看高清 | 男男h啪啪无遮挡| 久久久久国产精品人妻aⅴ院| 999久久久精品免费观看国产| 一级a爱片免费观看的视频| 51午夜福利影视在线观看| 最新在线观看一区二区三区| 免费日韩欧美在线观看| 久久久久久人人人人人| 国产亚洲欧美在线一区二区| 精品久久久久久电影网| 欧美乱码精品一区二区三区| 悠悠久久av| 亚洲熟妇熟女久久| 一区二区三区精品91| 大码成人一级视频| 伦理电影免费视频| 久久精品亚洲精品国产色婷小说| 777久久人妻少妇嫩草av网站| 国产精品久久久av美女十八| 久久99一区二区三区| 很黄的视频免费| 色播在线永久视频| 巨乳人妻的诱惑在线观看| 久久亚洲真实| 国产人伦9x9x在线观看| 一进一出抽搐gif免费好疼 | 露出奶头的视频| 久久天堂一区二区三区四区| 午夜免费激情av| 丁香六月欧美| 国产真人三级小视频在线观看| 天堂动漫精品| 亚洲一卡2卡3卡4卡5卡精品中文| 黄片小视频在线播放| 叶爱在线成人免费视频播放| 亚洲免费av在线视频| 高潮久久久久久久久久久不卡| 国产亚洲av高清不卡| 久久精品人人爽人人爽视色| 亚洲成国产人片在线观看| 国产精品乱码一区二三区的特点 | 伊人久久大香线蕉亚洲五| 国产成人精品久久二区二区91| 多毛熟女@视频| 999精品在线视频| 女性生殖器流出的白浆| 国产精品亚洲av一区麻豆| 国产精品国产高清国产av| 91麻豆精品激情在线观看国产 | 水蜜桃什么品种好| 亚洲激情在线av| 99热国产这里只有精品6| 国产成人免费无遮挡视频| 日韩欧美国产一区二区入口| 久久久精品欧美日韩精品| 在线观看日韩欧美| 99精品在免费线老司机午夜| 日韩欧美国产一区二区入口| 日韩av在线大香蕉| 脱女人内裤的视频| 黄色 视频免费看| 亚洲国产欧美日韩在线播放| 老熟妇仑乱视频hdxx| 久久人妻熟女aⅴ| 亚洲精品国产一区二区精华液| 久久精品91蜜桃| 亚洲av第一区精品v没综合| 午夜免费鲁丝| 老司机亚洲免费影院| 岛国在线观看网站| 老司机深夜福利视频在线观看| 午夜成年电影在线免费观看| 男男h啪啪无遮挡| 国产精品影院久久| 在线观看免费日韩欧美大片| 国产黄色免费在线视频| 黄色毛片三级朝国网站| 亚洲伊人色综图| cao死你这个sao货| 亚洲 欧美一区二区三区| 99国产极品粉嫩在线观看| 午夜福利在线观看吧| 国产99久久九九免费精品| 激情视频va一区二区三区| 色综合婷婷激情| 亚洲激情在线av| 很黄的视频免费| 好看av亚洲va欧美ⅴa在| 精品福利永久在线观看| 免费一级毛片在线播放高清视频 | 亚洲 欧美 日韩 在线 免费| 一夜夜www| 丝袜美足系列| 亚洲国产看品久久| 99久久精品国产亚洲精品| 97人妻天天添夜夜摸| 亚洲av日韩精品久久久久久密| 精品一区二区三区四区五区乱码| 在线观看舔阴道视频| 岛国视频午夜一区免费看| av网站在线播放免费| 国产三级在线视频| 国产国语露脸激情在线看| 好看av亚洲va欧美ⅴa在| 黄色视频不卡| 正在播放国产对白刺激| 国产真人三级小视频在线观看| 久久久国产成人免费| 成人特级黄色片久久久久久久| 狠狠狠狠99中文字幕| 成人精品一区二区免费| 亚洲av成人av| 丰满迷人的少妇在线观看| 精品久久久久久,| 色婷婷av一区二区三区视频| 少妇 在线观看| 精品国产乱码久久久久久男人| 久久国产精品人妻蜜桃| 80岁老熟妇乱子伦牲交| www国产在线视频色| 黑人巨大精品欧美一区二区蜜桃| 久久亚洲真实| 色综合婷婷激情| 久久久水蜜桃国产精品网| 黄色丝袜av网址大全| 亚洲国产精品999在线| 精品一区二区三区四区五区乱码| 国产精品免费视频内射| 天堂中文最新版在线下载| 国产精品永久免费网站| 国产一区在线观看成人免费| 久久久精品国产亚洲av高清涩受| 国产精品一区二区精品视频观看| 亚洲狠狠婷婷综合久久图片| xxxhd国产人妻xxx| 丰满的人妻完整版| av中文乱码字幕在线| 亚洲人成77777在线视频| 女人被狂操c到高潮| 免费看a级黄色片| 99香蕉大伊视频| 中文亚洲av片在线观看爽| 亚洲av日韩精品久久久久久密| 国产精品永久免费网站| 国产xxxxx性猛交| 色婷婷av一区二区三区视频| 侵犯人妻中文字幕一二三四区| 黄片小视频在线播放| 国内久久婷婷六月综合欲色啪| 国产乱人伦免费视频| 日本免费a在线| 国产在线观看jvid| 国产成人精品在线电影| 欧洲精品卡2卡3卡4卡5卡区| a级毛片黄视频| 亚洲,欧美精品.| 日本三级黄在线观看| 18禁观看日本| 亚洲人成网站在线播放欧美日韩| 嫩草影视91久久| 天堂√8在线中文| 18禁裸乳无遮挡免费网站照片 | 99久久国产精品久久久| 精品国产国语对白av| 777久久人妻少妇嫩草av网站| 亚洲欧美日韩无卡精品| 国产av精品麻豆| 欧美大码av| 亚洲精品在线美女| 亚洲国产欧美一区二区综合| 一本综合久久免费| 国产成人系列免费观看| 亚洲欧美激情综合另类| ponron亚洲| 午夜福利免费观看在线| 午夜免费鲁丝| 手机成人av网站| 国产成人精品久久二区二区91| 国产av又大| 人人澡人人妻人| 老司机深夜福利视频在线观看| 超碰97精品在线观看| 一区二区三区精品91| 多毛熟女@视频| 在线视频色国产色| 亚洲性夜色夜夜综合| 可以在线观看毛片的网站| 97碰自拍视频| 国产精品美女特级片免费视频播放器 | 日本三级黄在线观看| 久久久国产欧美日韩av| 99精品欧美一区二区三区四区| 在线观看免费日韩欧美大片| 9191精品国产免费久久| 在线观看一区二区三区激情| 国产一区二区激情短视频| 亚洲精华国产精华精| 超碰成人久久| videosex国产| 久久香蕉国产精品| 国产无遮挡羞羞视频在线观看| 99香蕉大伊视频| 亚洲专区中文字幕在线| 久久中文字幕一级| 免费在线观看视频国产中文字幕亚洲| 欧美午夜高清在线| 欧美乱妇无乱码| svipshipincom国产片| 久久精品国产清高在天天线| 国产av一区在线观看免费| 亚洲精品一卡2卡三卡4卡5卡| 级片在线观看| 亚洲精品成人av观看孕妇| 欧美激情极品国产一区二区三区| 久久久国产成人精品二区 | 精品少妇一区二区三区视频日本电影| 欧美黑人欧美精品刺激| 成人18禁在线播放| 日本五十路高清| 最新在线观看一区二区三区| 亚洲 欧美 日韩 在线 免费| 午夜免费激情av| 国产成人啪精品午夜网站| 国产不卡一卡二| 亚洲人成电影观看| 国产伦一二天堂av在线观看| 国产成人精品在线电影| 少妇的丰满在线观看| 亚洲美女黄片视频| 国产成人影院久久av| 777久久人妻少妇嫩草av网站| 男女床上黄色一级片免费看| 亚洲在线自拍视频| 老熟妇乱子伦视频在线观看| 中文字幕人妻丝袜制服| 在线观看免费午夜福利视频| 激情视频va一区二区三区| 久久午夜综合久久蜜桃| 中文字幕av电影在线播放| 看片在线看免费视频| 久久人人精品亚洲av| 一二三四社区在线视频社区8| 欧美精品一区二区免费开放| 一级黄色大片毛片| cao死你这个sao货| 999久久久精品免费观看国产| 90打野战视频偷拍视频| 中文字幕精品免费在线观看视频| 亚洲专区中文字幕在线| 精品午夜福利视频在线观看一区| 久久国产精品男人的天堂亚洲| 久久天躁狠狠躁夜夜2o2o| 在线观看免费视频日本深夜| 成在线人永久免费视频| 欧美 亚洲 国产 日韩一| 亚洲片人在线观看| 国产激情久久老熟女| 国产一区二区三区视频了| 三级毛片av免费| 老司机亚洲免费影院| 在线看a的网站| 麻豆av在线久日| 99国产精品一区二区三区| 免费看十八禁软件| 搡老乐熟女国产| 在线免费观看的www视频| av福利片在线| 亚洲精品成人av观看孕妇| 男女下面插进去视频免费观看| 国产熟女xx| 成人亚洲精品一区在线观看| 国产精品一区二区免费欧美| 亚洲av美国av| 午夜老司机福利片| 久久欧美精品欧美久久欧美| 国产色视频综合| 日韩欧美国产一区二区入口| 女人被狂操c到高潮| 欧美一区二区精品小视频在线| 久久久久久久午夜电影 | 精品一区二区三区四区五区乱码| 老熟妇乱子伦视频在线观看| 国产精品爽爽va在线观看网站 | 欧美一区二区精品小视频在线| 视频在线观看一区二区三区| 九色亚洲精品在线播放| 日韩欧美在线二视频| 亚洲精品美女久久av网站| www.www免费av| 黄色毛片三级朝国网站| 国产成年人精品一区二区 | 日韩一卡2卡3卡4卡2021年| 一本综合久久免费| 一个人观看的视频www高清免费观看 | 亚洲专区国产一区二区| 精品一区二区三卡| 一级,二级,三级黄色视频| 精品日产1卡2卡|