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

    露天礦富水裂隙巖體臺階爆破的殉爆機(jī)理和防殉爆研究

    2025-02-19 00:00:00費(fèi)鴻祿王天恒荊廣杰
    爆炸與沖擊 2025年1期
    關(guān)鍵詞:臨界壓力

    摘要: 殉爆現(xiàn)象會影響露天礦臺階爆破作業(yè)安全、邊坡穩(wěn)定性和爆破效果。在炸藥沖擊起爆機(jī)理基礎(chǔ)上,并結(jié)合露天礦實(shí)際富水裂隙巖體臺階爆破振動監(jiān)測結(jié)果,通過對比爆破振動信號波動差異來判別殉爆現(xiàn)象。為研究殉爆產(chǎn)生的機(jī)理和防殉爆方法,采用數(shù)值模擬和現(xiàn)場試驗(yàn)分析主發(fā)藥量、裂隙寬度及藥包之間的距離等參數(shù)對被發(fā)藥包孔壁壓力的影響。結(jié)果表明:孔壁沖擊壓力隨著裝藥耦合系數(shù)的減小、炮孔間裂隙寬度(0.25~1.00 cm)的增大以及炮孔間距離的減小而提高。在裂隙位置裝藥使用阻波管、充填巖粉或設(shè)置空氣間隔器,能顯著降低通過富水裂隙傳遞到被發(fā)炮孔的沖擊壓力,并使其低于乳化炸藥的起爆壓力臨界值。當(dāng)炮孔內(nèi)只有單條裂隙時,選擇填充巖粉是便捷且有效防殉爆方法;當(dāng)炮孔內(nèi)有多條裂隙時,該試驗(yàn)條件下,炮孔內(nèi)放置厚度為2.6 mm 的阻波管是最佳防殉爆方法,并能保證爆破效果。

    關(guān)鍵詞: 殉爆;臺階爆破;富水裂隙巖體;臨界壓力

    中圖分類號: O383 國標(biāo)學(xué)科代碼: 13505 文獻(xiàn)標(biāo)志碼: A

    影響被發(fā)藥包殉爆的因素主要有爆炸產(chǎn)生的爆轟波、沖擊波、沖擊波的傳播介質(zhì)、溫度、作用時間以及被發(fā)炸藥的感度。炸藥的起爆一方面受到爆炸時釋放熱量的影響,另一方面受到?jīng)_擊波壓力的作用[1-2]。李錚等[3] 通過試驗(yàn)發(fā)現(xiàn),對于不敏感炸藥,爆轟產(chǎn)物的沖擊是導(dǎo)致被發(fā)炸藥產(chǎn)生熱點(diǎn)并引發(fā)爆炸或爆轟的主要原因。爆轟產(chǎn)物是一種混合物,包含固體顆粒(如炭顆粒)和多種氣體(如N2、CO2、CO、H2O 和NOx 等);而敏感炸藥殉爆主要是由于沖擊波超壓和正壓作用時間引發(fā)的爆轟。炸藥的殉爆過程通常經(jīng)歷燃燒、加速燃燒和爆炸3 個階段。當(dāng)被發(fā)裝藥受到主發(fā)裝藥爆炸產(chǎn)生的爆轟產(chǎn)物作用時,其表面溫度升高,導(dǎo)致局部分解。分解產(chǎn)生的熱量引發(fā)炸藥燃燒,燃燒釋放的熱量進(jìn)一步提高溫度,加速燃燒,并通過孔隙傳遞至炸藥內(nèi)部;同時,沖擊波壓縮孔隙內(nèi)的空氣,形成熱點(diǎn),進(jìn)一步發(fā)生急劇的化學(xué)反應(yīng),發(fā)生加速燃燒直至整個炸藥爆炸[4]?;诒_擊波壓力殉爆理論[5] 可知,當(dāng)主發(fā)藥包爆炸產(chǎn)生的沖擊波壓力大于等于炸藥臨界起爆壓力時,被發(fā)藥包將會被引爆。因此,判定被發(fā)炮孔內(nèi)炸藥是否殉爆需確定炸藥起爆的臨界壓力。

    關(guān)于露天巖體爆破殉爆事件的研究很少,爆破中即使出現(xiàn)幾個孔殉爆也沒有引起注意,且值得借鑒研究方法和防止殉爆的技術(shù)也鮮有報道。余德運(yùn)等[6] 為了檢驗(yàn)空氣間隔中,混裝銨油炸藥(ammoniumnitrate fuel oil, ANFO)主發(fā)藥包殉爆被發(fā)藥包的可能性,設(shè)計殉爆試驗(yàn)并通過實(shí)驗(yàn)結(jié)果得出孔徑為311 mm 的炮孔比孔徑為165 mm 的炮孔殉爆距離更大。Zhang 等[7] 和Yang 等[8] 通過試驗(yàn)研究TNT 炸藥在水下爆破時,采用壓力演化、壓力峰值以及脈沖周期這3 種方式來評估炸藥殉爆,得到炸藥殉爆的臨界距離和安全距離分別為99.55 和128.92 mm。姜穎資等[9] 使用點(diǎn)火增長模型對3 種不同運(yùn)動速度的TNT 和PBX9404 主發(fā)炸藥進(jìn)行數(shù)值模擬研究,發(fā)現(xiàn)主發(fā)炸藥的運(yùn)動速度越快,破壞能力越強(qiáng),被發(fā)炸藥的臨界殉爆距離越大。Shin 等[10] 從撞擊產(chǎn)生的沖擊波誘發(fā)炸藥爆炸的角度,研究了不同約束材料對沖擊波誘發(fā)炸藥爆炸的影響, 得出增大約束材料的彈性阻抗會降低由于沖擊波產(chǎn)生殉爆的概率。Starkenberg 等[11] 對彈丸撞擊、殉爆和其他危險場景的進(jìn)行了模擬研究,發(fā)現(xiàn)爆炸部件上的沖擊載荷與用于校準(zhǔn)許多沖擊起爆模型的楔形試驗(yàn)中產(chǎn)生的沖擊載荷有顯著差異。李凱等[12] 通過數(shù)值模擬以及試驗(yàn)驗(yàn)證,發(fā)現(xiàn)爆炸沖擊波經(jīng)隔板衰減后起爆炸藥的能力隨著隔板中沖擊波傳播速度的提高而提高。李順波等[13] 利用ANSYS/LS-DYNA 軟件研究了爆炸沖擊波在混凝土、水、土中的傳播衰減規(guī)律,發(fā)現(xiàn)介質(zhì)的波阻抗對沖擊波初始壓力峰值有很大影響,土中沖擊波的持續(xù)時間最長,混凝土中沖擊波衰減最快,水中沖擊波能量損耗最低。趙根等[14] 通過理論推導(dǎo)以及水中炸藥爆炸的沖擊波測試,得到了較符合實(shí)際情況的水中沖擊波傳播規(guī)律?;▽毩岬萚15] 首次將拉氏量計用于研究乳化炸藥的沖擊起爆過程,同時針對乳化炸藥的沖擊起爆特點(diǎn)建立了反應(yīng)速率方程,并且通過試驗(yàn)得到乳化炸藥的臨界起爆壓力為2 GPa。李建軍等[ 1 6 ] 通過開展乳化炸藥的沖擊起爆試驗(yàn)得到乳化炸藥的沖擊起爆臨界能量值在(13~35)×1012 Pa2·s 之間。馬晶晶等[17] 和陳慶凱等[18] 通過殉爆測試,得到了乳化炸藥在不同約束條件下的殉爆距離以及2 種不耦合系數(shù)下的殉爆距離。陳朗等[19] 和張所碩等[20] 開展了殼裝固黑鋁炸藥殉爆試驗(yàn),通過殘留炸藥、見證板和被發(fā)裝藥殼體破壞情況,判斷裝藥的殉爆反應(yīng)等級,從而確定炸藥的臨界殉爆距離;并通過數(shù)值模擬分析了殼體厚度對裝藥殉爆距離的影響。胡宏偉等[21] 和劉曉文等[22] 開展了殼裝炸藥水下殉爆試驗(yàn),獲得了水下裝藥的殉爆距離和殉爆安全距離,并基于氣泡能分析了被發(fā)裝藥的反應(yīng)率。Ko 等[23] 研究了水下裝藥結(jié)構(gòu)的殉爆特點(diǎn)和機(jī)理,得到了兩藥包間的臨界距離,解決了水下聚能爆破的殉爆問題。

    目前,解決殉爆的研究主要在軍事、炸藥生產(chǎn)、炸藥儲存和運(yùn)輸?shù)确矫妫瑢τ诘V山方面的殉爆問題及其危害處理未見報道,也沒有對露天臺階爆破產(chǎn)生的殉爆現(xiàn)象提出具體科學(xué)、合理、系統(tǒng)的分析和研究方法及用于防止殉爆的具體措施。因此,本文中通過數(shù)值模擬和現(xiàn)場試驗(yàn),探討富水裂隙巖體的爆破殉爆機(jī)理,給出現(xiàn)場易殉爆區(qū)域降低(防止)殉爆的具體技術(shù)和施工措施。

    1 基于爆破振動監(jiān)測的現(xiàn)場工況與殉爆因素分析

    1.1 工況

    元寶山露天礦走向長度為4.68 km;傾向?qū)挾绕骄鶠?.56 km,總面積約為12 km2。地理地形來源于燕山山脈向北部延伸的部分,地貌形態(tài)以低山丘陵為主,整體地勢呈西南高、東北低,西北部有固定和半固定沙丘間隔分布。英金河和老哈河是元寶山礦區(qū)表面上的2 條主干支流,英金河貫穿整個礦區(qū)的中間地帶,地下水系分布豐富,如圖1 所示。

    2020 年以來,內(nèi)蒙古平莊煤業(yè)元寶山露天礦爆破開采時連續(xù)出現(xiàn)多次殉爆情況,殉爆產(chǎn)生導(dǎo)致爆破區(qū)域周圍振動過大。通過比較所收集的殉爆區(qū)域地質(zhì)資料,得知殉爆出現(xiàn)的位置集中在采場作業(yè)區(qū)間2 200 線以西、344~404 平盤,此區(qū)域地質(zhì)條件為細(xì)砂巖、粉砂巖,巖層裂隙較發(fā)育,殉爆位置周圍坡面均有水點(diǎn)出現(xiàn)。

    1.2 殉爆因素

    元寶山露天礦每次殉爆情況都出現(xiàn)于巖石臺階,此區(qū)域巖石巖性為細(xì)砂巖、粉砂巖,巖層內(nèi)含水量大。巖石臺階中節(jié)理裂隙發(fā)育復(fù)雜,巖層中存在裂縫,裂縫可能貫通相鄰炮孔。

    (1) 元寶山露天礦易殉爆區(qū)巖體存在一定寬度的貫通富水裂隙

    元寶山露天礦位于老哈河與英金河的交匯處,地層裂隙水豐富。對易殉爆區(qū)進(jìn)行現(xiàn)場勘察,發(fā)現(xiàn)易殉爆區(qū)巖體完整性較差,裂隙發(fā)育程度更高(圖2~3),為沖擊波傳遞至被發(fā)炮孔提供了有利條件。

    (2) 小藥量爆破除水

    現(xiàn)場使用的炸藥為2#巖石乳化炸藥,其密度為0.95~1.30 g/cm3,與炮孔內(nèi)水的密度接近,導(dǎo)致炸藥不易放置孔底,因此進(jìn)行裝藥工序時需要將炮孔內(nèi)的積水排出?,F(xiàn)場采用爆破除水的方式進(jìn)行排水,由于爆區(qū)巖質(zhì)較軟且裂隙較為發(fā)育,爆破除水時可能進(jìn)一步引起相鄰炮孔裂隙擴(kuò)大甚至貫通,易引起現(xiàn)場爆破作業(yè)時殉爆的發(fā)生。

    (3) 寶馬礦采空區(qū)瓦斯爆炸

    2016 年12 月3 日,赤峰寶馬礦業(yè)越界開采元寶山露天煤礦資源時發(fā)生特別重大瓦斯爆炸事故。經(jīng)調(diào)查發(fā)現(xiàn),該礦業(yè)采用房柱式炮采、長壁炮采、掘進(jìn)機(jī)巷采等方式越界盜采煤炭資源,實(shí)際形成的采空區(qū)(見圖4)總面積約為1.49 km2。經(jīng)現(xiàn)場實(shí)地勘察,發(fā)現(xiàn)元寶山露天礦易殉爆區(qū)在“寶馬礦”采空區(qū)上方區(qū)域,推測瓦斯爆炸進(jìn)一步使采空區(qū)頂板上覆巖體裂隙擴(kuò)張。

    1.3 基于爆破振動監(jiān)測的殉爆判別

    GB 6722—2014《爆破安全規(guī)程》[24] 規(guī)定爆破質(zhì)點(diǎn)振動速度計算采用Sadovski 公式[25]:

    v = K(Q1/3/R)α (1)

    式中:v 為臺階質(zhì)點(diǎn)峰值振速;R 為爆破振動安全允許的爆心距;Q 為單段起爆最大藥量;K 和α 分別為與地質(zhì)條件以及高程相關(guān)的系數(shù)和衰減指數(shù),結(jié)合現(xiàn)場地質(zhì)條件,K=150~250,α=1.5~1.8。

    在逐孔起爆時,爆破條件、爆破參數(shù)和爆心距R 不變的條件下:如果不發(fā)生殉爆,根據(jù)式(1) 計算得到的藥量僅為單孔藥量Q;如果n 個炮孔發(fā)生殉爆,則會導(dǎo)致一次起爆藥量增大至nQ,質(zhì)點(diǎn)峰值振速為單孔質(zhì)點(diǎn)峰值振速的nα/3 倍,質(zhì)點(diǎn)峰值振速呈指數(shù)增大。因此,通過對監(jiān)測得到的爆破振動波形中的峰值振速與計算得到的峰值振速進(jìn)行比較,如果監(jiān)測得到的質(zhì)點(diǎn)峰值振速顯著增大,可認(rèn)定發(fā)生了殉爆。

    基于1.2 節(jié)中的殉爆因素分析,對含裂隙水巖體和無裂隙水巖體區(qū)域分別進(jìn)行爆破試驗(yàn),具體爆破參數(shù)如表1 所示。

    為判別現(xiàn)場爆破作業(yè)時是否發(fā)生殉爆,在距爆心68.00、80.89、121.01 和171.69 m 處均放置TC-4 850測振儀,測得的質(zhì)點(diǎn)峰值振動速度如圖5 所示。

    由圖5 可知,第1 組試驗(yàn)測得的質(zhì)點(diǎn)峰值振動速度均大于第2 組相同爆心距處的質(zhì)點(diǎn)峰值振動速度。由于篇幅有限,以距爆心80.89 m的測點(diǎn)為例進(jìn)行分析。將2 組試驗(yàn)中距爆心均為80.89 m 的測點(diǎn)處的Z 方向質(zhì)點(diǎn)振動速度進(jìn)行比較,如圖6 所示。由圖6 可以看出,第1 組試驗(yàn)中距爆心8 0 . 8 9 m 的質(zhì)點(diǎn)峰值振速為12.42 cm/s,第2 組試驗(yàn)中距爆心80.89 m 的質(zhì)點(diǎn)峰值振速為3.25 cm/s。由式(1) 計算得知距爆心80.89 m 處的正常質(zhì)點(diǎn)峰值振速為0.97~3.76 cm/s,遠(yuǎn)小于第1 組試驗(yàn)中該爆心距處的質(zhì)點(diǎn)峰值振速12.42 cm/s,與第2 組試驗(yàn)中該爆心距處的質(zhì)點(diǎn)峰值振速3.25 cm/s 接近??梢?,第1 組試驗(yàn)區(qū)域發(fā)生了殉爆。從而驗(yàn)證了1.2 節(jié)中所述:當(dāng)巖體存在貫通富水裂隙時,易殉爆區(qū)會產(chǎn)生殉爆現(xiàn)象。

    2 富水裂隙巖體爆破殉爆的數(shù)值模擬

    2.1 材料參數(shù)

    (1)2#乳化炸藥的狀態(tài)方程及材料參數(shù)

    在數(shù)值模擬中選用與現(xiàn)場爆破一致的2#乳化炸藥,采用材料模型*MAT_HIGH_EXPLOSIVE_BURN 和Jones-Wilkins-Lee (JWL) 狀態(tài)方程[26] 表征2#乳化炸藥。針對2#乳化炸藥,JWL 狀態(tài)方程的具體形式為:

    式中: 為炸藥爆轟產(chǎn)物的壓力,V 為爆轟產(chǎn)物的相對比容,A、B、R1、R2 和ω 為待定參數(shù),Ee0 為炸藥的初始比內(nèi)能。2#乳化炸藥材料及狀態(tài)方程的參數(shù)[26] 如表2 所示, 為炸藥初始密度,D 為炸藥爆速。

    (2) 巖石材料模型

    選用Riedel-Hiermaier-Thoma (RHT) 模型[27] 表征巖石,巖石材料參數(shù)[26] 如表3 所示。

    (3) 水介質(zhì)狀態(tài)方程及材料參數(shù)

    選用材料模型*MAT_MULL 和Grüneisen 狀態(tài)方程[28] 表征水介質(zhì),相關(guān)表達(dá)式為:

    式中:c 為粒子速度曲線的截距,Ew 為水的初始體積內(nèi)能,ρw為水的初始密度,γ0為Grüneisen 系數(shù),a 為γ0的一階修正系數(shù),μ為壓縮比,S1、S2、S3 為粒子速度曲線的斜率參數(shù)。具體參數(shù)[28] 如表4 所示。

    2.2 模擬模型建立

    為驗(yàn)證巖體爆破數(shù)值模擬材料模型及參數(shù)的準(zhǔn)確性,利用有限元分析軟件,按照所選的材料本構(gòu)模型及參數(shù),建立一個尺寸為3 000 cm×1 800 cm×1 cm 的巖樣爆破準(zhǔn)三維數(shù)值模型(見圖7),在巖樣模型中心線兩側(cè)布置2 個孔距為600 cm 的炮孔,炮孔長度均為1 350 cm,炮孔直徑均為20 cm,兩炮孔間富水裂隙寬度為1 cm。裝藥長度和巖粉填塞長度分別為450 和400 cm。有限元模型共劃分單元430 362 個,節(jié)點(diǎn)865 956 個。為了模擬無限介質(zhì),將模型上邊界定義為自用面,左右兩側(cè)及底部邊界設(shè)置為無反射性邊界。模型采用流固耦合算法,其中炸藥、空氣及水介質(zhì)采用ALE (aribitrary Lagrange-Euler) 算法,巖石與填塞采用Lagrange 算法。

    3 殉爆機(jī)理分析與判定

    3.1 被發(fā)藥包殉爆位置節(jié)點(diǎn)確定

    由殉爆能量閾值[16] 可知,被發(fā)藥包是否發(fā)生殉爆與藥包內(nèi)任意有限元點(diǎn)在單位時間內(nèi)的壓力值相關(guān)。為進(jìn)一步判斷被發(fā)藥包是否達(dá)到殉爆臨界能量值,選擇貫穿裂隙單元與被發(fā)藥包接觸面的壓力時程曲線進(jìn)行分析。根據(jù)被發(fā)裝藥受沖擊方向,提取被發(fā)藥包內(nèi)部不同位置單元的壓力曲線。被發(fā)藥包壓力提取點(diǎn)分布如圖8 中A~I(xiàn) 所示。

    圖9 為無防護(hù)措施時貫穿富水裂隙條件下巖體爆破的應(yīng)力云圖。主發(fā)炮孔起爆后,在t=0.4 ms 時炸藥起爆后產(chǎn)生的沖擊波迅速向外擴(kuò)散,其形狀呈“橢圓形”,這說明爆炸產(chǎn)生的應(yīng)力波在巖體中的衰減與傳播是均勻的。在0.8 和1.0 ms 時,炮孔間裂隙處壓力值在沖擊波在傳播過程中明顯高于周圍巖體壓力值。當(dāng)炸藥起爆后,由于水介質(zhì)的不可壓縮性,爆炸沖擊波用于壓縮水介質(zhì)所消耗的能量占比小,導(dǎo)致沖擊波傳播至相鄰炮孔被發(fā)藥能量大,進(jìn)而使鄰近炮孔內(nèi)炸藥發(fā)生殉爆現(xiàn)象。

    根據(jù)圖8 中監(jiān)測點(diǎn)位置,提取被發(fā)藥包內(nèi)部單元的壓力時程曲線,如圖10 所示:在0~2 ms 時,各壓力曲線緩慢上升,2~3 ms 時,壓力曲線發(fā)生突躍變化,3 ms 以后迅速衰減,并在4 ms 后趨于平緩。其中在裂隙處最易殉爆單元F 壓力峰值最大,達(dá)到2.3 GPa,大于乳化炸藥的臨界起爆壓力[15]。即在貫穿裂隙處被發(fā)炸藥表面最有可能首先引起殉爆,進(jìn)而導(dǎo)致整體殉爆。

    3.2 不同參數(shù)對殉爆的影響

    根據(jù)3.1 節(jié)對殉爆過程的數(shù)值模擬分析,可以確定導(dǎo)致被發(fā)炮孔炸藥發(fā)生殉爆的關(guān)鍵因素是沖擊波對被發(fā)藥包施加的壓力超過了其臨界起爆壓力,或者沖擊波向炸藥傳遞的能量超過了被發(fā)炸藥的臨界起爆能量。在爆破過程中,被發(fā)炮孔中沖擊波壓力與主發(fā)炮孔中炸藥質(zhì)量、炮孔間的距離以及裂隙寬度有關(guān)。因此,進(jìn)行不同參數(shù)條件下被發(fā)藥包沖擊起爆的研究,具體的數(shù)值模擬方案如表5 所示,模擬結(jié)果如圖11~16 所示。

    方案Ⅰ:圖11 為炮孔間不同炮孔間距條件下巖體爆破的應(yīng)力云圖,繪制被發(fā)藥包內(nèi)部F 單元的峰值壓力隨炮孔間距變化的曲線,如圖12 所示。由圖11~12 可知,隨著炮孔間距的增大,被發(fā)藥包所受到的沖擊壓力降低。當(dāng)炮孔間距為6 m 時,被發(fā)藥包所受的沖擊壓力為2.27 GPa,大于乳化炸藥的臨界起爆壓力;因此,當(dāng)孔距小于6 m 時,被發(fā)藥包所受的沖擊壓力達(dá)到炸藥沖擊起爆條件,從而導(dǎo)致殉爆。

    方案Ⅱ:圖13 為2 ms 時炮孔間不同裂隙寬度條件下巖體爆破的應(yīng)力云圖,繪制被發(fā)藥包內(nèi)部F 單元的峰值壓力隨炮孔間裂隙寬度變化的曲線,如圖14 所示。由圖13~14 可知:當(dāng)裂縫寬度在0.25~1.00 cm 之間時,被發(fā)藥包的峰值壓力曲線呈現(xiàn)線性增長,并在裂縫寬度為1.00 cm時達(dá)到最大值,峰值為2.17 GPa;當(dāng)裂縫寬度為1.00~2.00 cm 之間時,隨著裂縫寬度的增大,被發(fā)藥包的峰值壓力曲線的下降趨勢逐漸變緩。在上述模擬條件下,炮孔之間裂隙寬度為1.00 cm時,被發(fā)藥包所受的沖擊壓力達(dá)到被發(fā)藥包的沖擊起爆條件。

    方案Ⅲ:圖15 為2 ms 時不同主發(fā)藥量條件下巖體爆破的應(yīng)力云圖,繪制被發(fā)藥包內(nèi)部F 單元的峰值壓力隨主發(fā)炮孔內(nèi)藥量變化的曲線,如圖16 所示。由圖15~16 可知,被發(fā)裝藥所受的沖擊壓力隨著主發(fā)藥量的增加而不斷升高,呈現(xiàn)一個先急速上升后增速變緩的趨勢。主發(fā)炮孔的藥量從20 kg 增加到120 kg,被發(fā)藥包所受到的峰值壓力增大了約1.3 倍。在上述模擬條件下,主發(fā)炮孔的藥量大于72 kg 時,被發(fā)藥包所受的沖擊壓力達(dá)到被發(fā)藥包的沖擊起爆條件。

    3.3 阻波管和巖粉防殉爆設(shè)計

    從安全性防護(hù)設(shè)計的角度出發(fā),難以控制現(xiàn)場爆破試驗(yàn)過程中主發(fā)炮孔炸藥爆炸產(chǎn)生的沖擊波壓力沿著富水裂隙傳播至被發(fā)炮孔的大小,但可通過在主發(fā)藥包外放置約束材料對被發(fā)藥包進(jìn)行防護(hù)。此外,由Shin 等[10] 的研究可知:不同約束材料對沖擊波誘發(fā)炸藥爆炸的影響,隨著約束材料的彈性阻抗增大,會減少由于沖擊波產(chǎn)生殉爆的機(jī)會。

    考慮到實(shí)際裝藥過程中,要求約束材料需具備質(zhì)量輕、成本低、材料來源廣泛及乳化炸藥取放方便等特點(diǎn),本文中主發(fā)炮孔內(nèi)的約束材料選擇2.0、2.6 和3.2 mm 等3 種不同厚度的阻波管。主發(fā)炮孔防殉爆設(shè)計如圖17~18 所示。

    針對元寶山露天礦臺階爆破出現(xiàn)的殉爆現(xiàn)象,提出了2 類防止殉爆的方案,即主發(fā)藥包內(nèi)放置阻波管和裂隙處填塞巖粉/空氣間隔器方案。前者針對主發(fā)炮孔內(nèi)多條裂隙情況下殉爆災(zāi)害的控制,而后者是針對主發(fā)炮孔內(nèi)單條裂隙情況下殉爆災(zāi)害的控制。上述2 種方案分別為:

    (1) 主發(fā)藥包內(nèi)放置阻波管防殉爆模擬方案:根據(jù)現(xiàn)場實(shí)際情況,模擬貫穿裂隙2 條,寬度均為1 cm,炮孔間距為6 m,底部裝藥量為96 kg,炸藥裝填長度為4.5 m,巖石碎屑填充長度為4 m。

    (2) 主發(fā)藥包內(nèi)裂隙處填巖粉防殉爆模擬方案:根據(jù)現(xiàn)場實(shí)際情況,模擬貫穿裂隙寬度為1 cm,炮孔間距為6 m,底部裝藥量為96 kg,炸藥裝填長度為4.5 m,巖石碎屑填充長度為4 m。

    3.3.1 不同厚度阻波管防殉爆結(jié)果分析

    為了分析不同厚度阻波管的防殉爆效果,在原模型基礎(chǔ)上(圖7)主發(fā)炮孔裂隙處添加2.0、2.6 和3.2 mm 等3 種不同厚度的阻波管(圖19),進(jìn)行數(shù)值計算,計算結(jié)果如圖20~21 所示。

    根據(jù)文獻(xiàn)[17-18] 可知,阻波管對炸藥爆炸沖擊波的傳播有約束作用,可以改變炸藥爆炸沖擊波的傳遞方向。由圖20 可知,隨著阻波管厚度的增大,其約束作用也越來越顯著,炮孔間沖擊波在傳播過程中相同時刻裂隙處的壓力值逐漸降低。這表明,阻波管厚度的增大,在一定程度上降低了爆破能量的傳遞,降低殉爆發(fā)生的概率。

    由圖21 可知,隨著阻波管厚度的增大,被發(fā)藥包內(nèi)部的壓力峰值不斷降低。當(dāng)阻波管厚度為3.2 mm 時,峰值壓力為835 MPa,與未加阻波管條件相比,被發(fā)藥包內(nèi)部的峰值壓力降低了73.71%。

    3.3.2 巖粉防殉爆結(jié)果分析

    為了分析在主發(fā)炮孔內(nèi)裂隙處填塞巖粉的防殉爆效果,在原模型基礎(chǔ)上(圖7)將主發(fā)炮孔裂隙處主發(fā)藥包替換為巖粉(圖22),進(jìn)行數(shù)值計算,結(jié)果如圖23~24 所示。

    在空氣間隔和填塞巖粉兩種不同填充條件下,炸藥起爆開始階段壓力值區(qū)別不明顯。由圖23(a) 可以發(fā)現(xiàn),t=1.4 ms 時,空氣間隔條件下炮孔間裂隙處出現(xiàn)了壓力集中現(xiàn)象。這是因?yàn)?,被空氣間隔的上、下兩部分炸藥爆炸后,沖擊波向中間運(yùn)動,當(dāng)兩者相遇時形成高壓區(qū),即表現(xiàn)為壓力集中現(xiàn)象。比較圖23(a) 和(b) 可以發(fā)現(xiàn),t=2.0 ms 時,在空氣間隔條件下裂隙處的壓力比巖粉填充條件下高。這表明空氣間隔中儲存的能量釋放使得裂隙處壓力升高。

    由圖24 可以看出,填塞巖粉模型的壓力時程曲線呈現(xiàn)一個波峰,而空氣間隔模型的壓力時程曲線出現(xiàn)2 個波峰。這是因?yàn)?,炸藥爆炸后,空氣因其可壓縮性暫時儲存能量,炸藥爆炸沖擊波傳播結(jié)束后,壓縮空氣儲存的能量釋放,釋放的能量通過裂隙水傳遞作用于被發(fā)藥包。在填塞巖粉的情況下,最易殉爆單元的峰值壓力為224 MPa,使用空氣間隔器時為276 MPa,均未達(dá)到乳化炸藥的起爆壓力2 GPa。

    綜合被發(fā)裝藥內(nèi)部壓力和應(yīng)力云圖判斷,在有防護(hù)措施的條件下,被發(fā)裝藥未發(fā)生殉爆。這表明,針對主發(fā)炮孔的防殉爆設(shè)計可有效降低殉爆發(fā)生的概率。

    4 現(xiàn)場試驗(yàn)驗(yàn)證

    4.1 殉爆檢測方法

    現(xiàn)場主發(fā)炮孔和被發(fā)炮孔深度均為13.5 m,主發(fā)炮孔裝藥量為96 kg,被發(fā)炮孔裝藥量為24 kg。以主發(fā)炮孔為中心,設(shè)置1 號被發(fā)炮孔與其距離為2 m,2~6 號被發(fā)炮孔與主發(fā)炮孔的距離依次遞增1 m,并且隨著距離的增大被發(fā)炮孔與主發(fā)炮孔的連線逆時針旋轉(zhuǎn)依次增大60°,如圖25 所示。為了避免試驗(yàn)結(jié)果的偶然性,在距離主發(fā)炮孔30、35 和40 m 處分別設(shè)置1 臺測振儀,對炸藥爆炸產(chǎn)生的振動進(jìn)行監(jiān)測。設(shè)置孔間起爆延遲時間為1 s,中心炮孔為瞬發(fā),按照逆時針對被發(fā)藥包依次起爆,當(dāng)被發(fā)炮孔內(nèi)炸藥未按照既定時間爆炸時,則判定該孔內(nèi)的炸藥發(fā)生了殉爆。

    4.2 主發(fā)藥包內(nèi)放置阻波管防殉爆試驗(yàn)

    根據(jù)殉爆檢測現(xiàn)場布置,分別在主發(fā)炮孔內(nèi)放置厚度為2.0、2.6、3.2 mm 的阻波管,并進(jìn)行現(xiàn)場殉爆試驗(yàn)和爆破振動監(jiān)測。提取監(jiān)測點(diǎn)①的振動數(shù)據(jù),繪制不同厚度阻波管的防殉爆試驗(yàn)振速時程曲線,如圖26 所示。由圖26(a) 可以看出,在3 s 時振動波形缺失。這表明,3 號被發(fā)炮孔內(nèi)發(fā)生了殉爆??梢姡褂煤穸?.0 mm 的阻波管仍存在殉爆風(fēng)險。而由圖26(b) 和(c) 可以看出,1~6 s 均發(fā)生了振動。這表明,6 個被發(fā)炮孔內(nèi)均未發(fā)生殉爆。試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果基本一致。

    4.3 主發(fā)炮孔裂隙處填塞巖粉/空氣間隔器防殉爆試驗(yàn)

    根據(jù)殉爆檢測現(xiàn)場布置,分別在主發(fā)炮孔內(nèi)裂隙處填塞巖粉或者使用空氣間隔器,并進(jìn)行現(xiàn)場殉爆試驗(yàn)和爆破振動監(jiān)測。提取監(jiān)測點(diǎn)①的振動數(shù)據(jù),繪制填塞巖粉和使用空氣間隔的防殉爆試驗(yàn)振速時程曲線圖如圖27 所示。由圖27 可以看出,1~6 s 時間點(diǎn)均發(fā)生了振動。這表明,當(dāng)主發(fā)炮孔內(nèi)部裂隙處填塞巖粉或者使用空氣間隔器后,即使被發(fā)藥包處于裂隙位置也未發(fā)生殉爆現(xiàn)象。可見,在炮孔內(nèi)的裂隙區(qū)域填塞巖粉或者使用空氣間隔器,避免了主發(fā)藥包產(chǎn)生的沖擊波直接作用于富水裂隙,降低了被發(fā)藥包受到的沖擊壓力,有效地避免了被發(fā)藥包殉爆的發(fā)生。

    試驗(yàn)結(jié)果表明:在主發(fā)炮孔裂隙處添加阻波管和巖粉/空氣間隔器可防止殉爆。

    5 結(jié) 論

    以元寶山露天礦殉爆問題為研究背景,開展了數(shù)值模擬和現(xiàn)場殉爆試驗(yàn),主要工作及結(jié)論如下。

    (1) 該露天礦臺階爆破殉爆產(chǎn)生的原因是炮孔間存在一定寬度的貫通富水裂隙,炸藥爆破產(chǎn)生的水沖擊波傳播至相鄰炮孔時,最大壓力值超過了乳化炸藥的臨界起爆壓力,引發(fā)殉爆現(xiàn)象。通過研究不同參數(shù)(炸藥質(zhì)量、炮孔間的距離以及裂隙寬度)對被發(fā)藥包沖擊起爆條件的影響,發(fā)現(xiàn)了孔壁沖擊壓力隨裝藥耦合系數(shù)的減小、裂隙寬度(0.25~1.00 cm)的增大以及炮孔間距的減小而提高。

    (2) 通過對3 種防止殉爆措施進(jìn)行模擬研究,發(fā)現(xiàn)被發(fā)藥包在3 種防止殉爆措施下被施加的壓力均低于臨界起爆壓力。結(jié)合現(xiàn)場實(shí)際情況,當(dāng)炮孔內(nèi)炸藥放置處存在單條裂隙時,采用填塞巖粉的方法來防止殉爆的發(fā)生;當(dāng)存在多條裂隙時,炮孔內(nèi)放置厚度為2.6 mm 阻波管的方法防止殉爆的發(fā)生。

    (3) 基于數(shù)值模擬和現(xiàn)場殉爆試驗(yàn),提出了現(xiàn)場判定殉爆的檢測方法;在易殉爆區(qū)主發(fā)炮孔裝藥使用阻波管和巖粉,可降低被發(fā)藥包所受的沖擊壓力,防止殉爆的發(fā)生。

    在調(diào)研露天礦臺階爆破殉爆現(xiàn)象和開展現(xiàn)場殉爆試驗(yàn)過程中,得到了內(nèi)蒙古平莊煤業(yè)(集團(tuán))有限責(zé)任公司元寶山露天煤礦的大力支持。礦方工程技術(shù)人員王晨光、王勇、張秋園、劉競超、敖然松和武浩不僅多次參與了線上、線下的討論會,還在試驗(yàn)現(xiàn)場進(jìn)行了精心的組織協(xié)調(diào)。他們的參與確保了整個研究過程的順利進(jìn)行,在此謹(jǐn)表謝意。

    參考文獻(xiàn):

    [1]王玉杰. 爆破工程 [M]. 武漢: 武漢理工大學(xué)出版社, 2007: 20–22.

    [2]DERIBAS A A, MEDVEDEV E A, RESHETNYAK Y A, et al. Detonation of emulsion explosives containing hollow microspheres [J]. Doklady Physics, 2003, 48(4): 163–165. DOI: 10.1134/1.1574370.

    [3]李錚, 項(xiàng)續(xù)章, 郭梓熙. 各種炸藥的殉爆安全距離 [J]. 爆炸與沖擊, 1994, 14(3): 231–241. DOI: 10.11883/1001-1455(1994)03-0231-11.

    LI Z, XIANG X Z, GUO Z X. Various explosives of safety distance of unsympathetic detonation [J]. Explosion and Shock Waves, 1994, 14(3): 231–241. DOI: 10.11883/1001-1455(1994)03-0231-11.

    [4] 費(fèi)鴻祿. 爆破理論及其應(yīng)用 [M]. 2 版. 北京: 煤炭工業(yè)出版社, 2018: 89–92.

    [5]汪成運(yùn), 魏志豐, 何鵬鵬. 炸藥殉爆的研究進(jìn)展與展望 [J]. 爆破器材, 2022, 51(6): 1–8. DOI: 10.3969/j.issn.1001-8352.2022.06.001.

    WANG C Y, WEI Z F, HE P P. Research progress of sympathetic detonation of explosives [J]. Explosive Materials, 2022,51(6): 1–8. DOI: 10.3969/j.issn.1001-8352.2022.06.001.

    [6]余德運(yùn), 謝烽, 王旭耀. ANFO 在炮孔中的殉爆起爆試驗(yàn)研究 [J]. 爆破器材, 2020, 49(5): 59–64. DOI: 10.3969/j.issn.1001-8352.2020.05.011.

    YU D Y, XIE F, WANG X Y. Experimental study on sympathetic detonation of ANFO in hole [J]. Explosive Materials, 2020,49(5): 59–64. DOI: 10.3969/j.issn.1001-8352.2020.05.011.

    [7]ZHANG Z F, WANG C, HU H L, et al. Investigation of underwater sympathetic detonation [J]. Propellants, Explosives,Pyrotechnics, 2020, 45(11): 1736–1744. DOI: 10.1002/prep.202000099.

    [8]YANG J X, SHI C, YANG W K, et al. Numerical simulation of column charge explosive in rock masses with particle flow code [J]. Granular Matter, 2019, 21(4): 96. DOI: 10.1007/s10035-019-0950-2.

    [9]姜穎資, 王偉力, 黃雪峰, 等. 帶殼炸藥在高速運(yùn)動炸藥作用下殉爆效應(yīng)研究 [J]. 工程爆破, 2014, 20(3): 1–4. DOI:10.3969/j.issn.1006-7051.2014.03.001.

    JIANG Y Z, WANG W L, HUANG X F, et al. Research on the sympathetic detonation effect of shelled explosive by highspeed movement explosive [J]. Engineering Blasting, 2014, 20(3): 1–4. DOI: 10.3969/j.issn.1006-7051.2014.03.001.

    [10]SHIN H, LEE W. Material design guidelines for explosive confinements to control impact shock-induced detonations based on shock transmission/reflection analysis [J]. International Journal of Impact Engineering, 2003, 28(5): 465–478. DOI: 10.1016/S0734-743X(2)00075-1.

    [11]STARKENBERG J, HUANG Y, ARBUCKLE A. Numerical modeling of projectile impact shock initiation of bare and covered composition-B [J]. Journal of Energetic Materials, 1984, 2(1/2): 1–41. DOI: 10.1080/07370658408012327.

    [12]李凱, 詹勇, 程波, 等. 爆炸沖擊波經(jīng)隔板衰減后的起爆能力數(shù)值研究 [C]//2014’(第六屆) 含能材料與鈍感彈藥技術(shù)學(xué)術(shù)研討會論文集. 北京: 中國兵工學(xué)會爆炸與安全技術(shù)專業(yè)委員會, 2014: 474–478.

    [13]李順波, 東兆星, 齊燕軍, 等. 爆炸沖擊波在不同介質(zhì)中傳播衰減規(guī)律的數(shù)值模擬 [J]. 振動與沖擊, 2009, 28(7): 115–117.DOI: 10.13465/j.cnki.jvs.2009.07.001.

    LI S B, DONG Z X, QI Y J, et al. Numerical simulation for spread decay of blasting shock wave in different media [J]. Journal of Vibration and Shock, 2009, 28(7): 115–117. DOI: 10.13465/j.cnki.jvs.2009.07.001.

    [14]趙根, 季榮, 鄭曉寧, 等. 乳化炸藥水中爆炸沖擊波傳播規(guī)律試驗(yàn)研究 [J]. 爆破, 2011, 28(2): 1–4. DOI: 10.3963/j.issn.1001-487X.2011.02.001.

    ZHAO G, JI R, ZHENG X N, et al. Experimental investigation on propagation rule of shock wave by emulsion explosives underwater blasting [J]. Blasting, 2011, 28(2): 1–4. DOI: 10.3963/j.issn.1001-487X.2011.02.001.

    [15]花寶玲, 李建軍, 丁淳彤. 乳化炸藥沖擊起爆過程的研究 [J]. 工程爆破, 1998, 4(1): 30–33.

    HUA B L, LI J J, DING C T. Study on shock initiation process for emulsion explosives [J]. Engineering Blasting, 1998, 4(1):30–33.

    [16]李建軍, 汪旭光, 歐育湘, 等. 乳化炸藥沖擊起爆的實(shí)驗(yàn)研究 [J]. 工程爆破, 1995, 1(1): 14–19.

    LI J J, WANG X G, OU Y X, et al. Experimental study on shock initiations of emulsion explosives [J]. Engineering Blasting,1995, 1(1): 14–19.

    [17]馬晶晶, 龍運(yùn)杰, 唐虹靖, 等. 炮孔約束下炸藥殉爆距離試驗(yàn)研究 [J]. 采礦技術(shù), 2023, 23(6): 165–169. DOI: 10.13828/j.cnki.ckjs.2023.06.037.

    MA J J, LONG Y J, TANG H J, et al. Experimental study on explosive detonation distance under blast hole constraints [J].Mining Technology, 2023, 23(6): 165–169. DOI: 10.13828/j.cnki.ckjs.2023.06.037.

    [18]陳慶凱, 夏亞偉, 劉占富, 等. 約束條件對乳化炸藥殉爆距離影響的研究 [J]. 礦業(yè)研究與開發(fā), 2017, 37(1): 45–49. DOI:10.13827/j.cnki.kyyk.2017.01.011.

    CHEN Q K, XIA Y W, LIU Z F, et al. The effect of constraint conditions on the gap distance of emulsion explosive [J].Mining Research and Development, 2017, 37(1): 45–49. DOI: 10.13827/j.cnki.kyyk.2017.01.011.

    [19]陳朗, 王晨, 魯建英, 等. 炸藥殉爆實(shí)驗(yàn)和數(shù)值模擬 [J]. 北京理工大學(xué)學(xué)報, 2009, 29(6): 497–500,524. DOI: 10.15918/j.tbit1001-0645.2009.06.004.

    CHEN L, WANG C, LU J Y, et al. Experiment simulation of sympathetic detonation tests [J]. Transactions of Beijing Institute of Technology, 2009, 29(6): 497–500,524. DOI: 10.15918/j.tbit1001-0645.2009.06.004.

    [20]張所碩, 聶建新, 張劍, 等. 約束空間內(nèi)殼裝炸藥殉爆及防護(hù) [J]. 爆炸與沖擊, 2023, 43(8): 085101. DOI: 10.11883/bzycj-2022-0456.

    ZHANG S S, NIE J X, ZHANG J, et al. Sympathetic detonation of explosive charge in confined space and its protection [J].Explosion and Shock Waves, 2023, 43(8): 106–119. DOI: 10.11883/bzycj-2022-0456.

    [21]胡宏偉, 王健, 卞云龍, 等. 帶殼裝藥水中殉爆特性分析 [J]. 水下無人系統(tǒng)學(xué)報, 2022, 30(3): 308–313. DOI: 10.11993/j.issn.2096-3920.2022.03.005.

    HU H W, WANG J, BIAN Y L, et al. Experiments of sympathetic detonation performance of explosives with shell in water [J].Journal of Unmanned Undersea Systems, 2022, 30(3): 308–313. DOI: 10.11993/j.issn.2096-3920.2022.03.005.

    [22]劉曉文, 高玉剛. 炸藥在水介質(zhì)中殉爆特性分析 [J]. 工程爆破, 2022, 28(4): 102–107. DOI: 10.19931/j.EB.20210241.

    LIU X W, GAO Y G. Analysis of explosive martyrdom in water medium [J]. Engineering Blasting, 2022, 28(4): 102–107.DOI: 10.19931/j.EB.20210241.

    [23]KO Y H, KIM S J, YANG H S. Assessment for the sympathetic detonation characteristics of underwater shaped charge [J].Geosystem Engineering, 2017, 20(5): 286–293. DOI: 10.1080/12269328.2017.1323679.

    [24]"中華人民共和國國家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局. 爆破安全規(guī)程: GB 6722–2014 [S]. 北京: 冶金工業(yè)出版社, 2014.

    [25]張忠偉, 任舸, 李洪濤. 錦屏二級水電站導(dǎo)流隧洞進(jìn)口圍堰拆除爆破 [J]. 爆破, 2011, 28(4): 77–80. DOI: 10.3963/j.issn.1001-487X.2011.04.021.

    ZHANG Z W, REN G, LI H T. Explosive demolition of intake cofferdam of diversion tunnel on Jingping Ⅱ Hydropower Station [J]. Blasting, 2011, 28(4): 77–80. DOI: 10.3963/j.issn.1001-487X.2011.04.021.

    [26]李文焱. 元寶山露天礦富水裂隙臺階爆破的殉爆機(jī)理及防止殉爆技術(shù)研究 [D]. 遼寧工程技術(shù)大學(xué), 2023.DOI:10.27210/d.cnki.glnju.2023.000910.

    [27]凌天龍, 王宇濤, 劉殿書 等. 修正RHT 模型在巖體爆破響應(yīng)數(shù)值模擬中的應(yīng)用 [J]. 煤炭學(xué)報, 2018, 43(S2): 434–442.DOI: 10.13225/j.cnki.jccs.2017.1698.

    LING T L, WANG Y T, LIU D S, et al. Modified RHT model for numerical simulation of dynamic response of rock mass under blasting load [J]. Journal of China Coal Society, 2018, 43(S2): 434–442. DOI: 10.13225/j.cnki.jccs.2017.1698.

    [28]SHINY S, LEE M, LAM K Y, et al. Modeling mitigation effects of watershield on shock waves [J]. Shock and Vibration,1998, 5(4): 225–234. DOI: 10.1155/1998/782032.

    (責(zé)任編輯 張凌云)

    基金項(xiàng)目: 遼寧省教育廳基本科研項(xiàng)目(青年項(xiàng)目)基金(JYTQN2023206)

    猜你喜歡
    臨界壓力
    加勁環(huán)式地下埋管非線性有限元屈曲分析
    易拉罐容器外壓失穩(wěn)實(shí)驗(yàn)裝置設(shè)計與應(yīng)用
    僅考慮自重的細(xì)長受彎構(gòu)件是否需滿足長細(xì)比要求的研究
    基于燒結(jié)物塌陷壓力檢測預(yù)測成品關(guān)鍵指標(biāo)①
    聚乙烯絕緣CO2物理發(fā)泡淺析
    科技資訊(2019年8期)2019-06-18 01:06:18
    封閉采空區(qū)瓦斯涌出影響因素及防治措施
    水下環(huán)肋功能梯度材料圓柱殼穩(wěn)定性研究
    縱骨對環(huán)肋圓柱殼肋間殼板穩(wěn)定性的影響
    隧洞岔管高壓壓水試驗(yàn)及成果分析
    煤粉壓力燃燒特性及動力學(xué)分析
    国产极品粉嫩免费观看在线 | 国产视频首页在线观看| 美女cb高潮喷水在线观看| 亚洲精品中文字幕在线视频 | 亚州av有码| av国产久精品久网站免费入址| 国产亚洲一区二区精品| 黑人巨大精品欧美一区二区蜜桃 | 丝袜脚勾引网站| 国模一区二区三区四区视频| 九色成人免费人妻av| 伊人亚洲综合成人网| 日韩av不卡免费在线播放| 寂寞人妻少妇视频99o| 日韩av不卡免费在线播放| 亚洲久久久国产精品| 欧美日本中文国产一区发布| 国产成人精品一,二区| a级片在线免费高清观看视频| 高清午夜精品一区二区三区| 蜜臀久久99精品久久宅男| 午夜日本视频在线| 看非洲黑人一级黄片| 色吧在线观看| 国产精品欧美亚洲77777| av国产精品久久久久影院| 久久久久网色| 亚洲精品国产av蜜桃| 丝袜喷水一区| 国产精品99久久久久久久久| 亚洲欧美一区二区三区黑人 | 欧美人与善性xxx| 人妻制服诱惑在线中文字幕| 天天操日日干夜夜撸| 搡女人真爽免费视频火全软件| a级毛片免费高清观看在线播放| 国产成人免费无遮挡视频| 国产老妇伦熟女老妇高清| 婷婷色av中文字幕| av在线app专区| 久久鲁丝午夜福利片| 国产精品国产三级专区第一集| 精品午夜福利在线看| 嫩草影院新地址| 水蜜桃什么品种好| 夜夜爽夜夜爽视频| 夫妻性生交免费视频一级片| av在线观看视频网站免费| 在线观看三级黄色| 日日爽夜夜爽网站| 久久人人爽人人片av| 一级爰片在线观看| 九草在线视频观看| 黄片无遮挡物在线观看| 亚洲国产精品专区欧美| 天堂8中文在线网| 欧美日韩av久久| 又黄又爽又刺激的免费视频.| 99国产精品免费福利视频| 亚洲欧美一区二区三区黑人 | 午夜免费鲁丝| 国产美女午夜福利| 中文字幕免费在线视频6| 免费看不卡的av| 你懂的网址亚洲精品在线观看| 3wmmmm亚洲av在线观看| 777米奇影视久久| 91久久精品国产一区二区成人| 性色av一级| 久久婷婷青草| 亚洲中文av在线| 久久免费观看电影| 亚洲精品色激情综合| 国模一区二区三区四区视频| 亚洲av不卡在线观看| 在线播放无遮挡| 永久免费av网站大全| 日本av免费视频播放| 五月玫瑰六月丁香| 在线观看美女被高潮喷水网站| 亚洲综合精品二区| 丝瓜视频免费看黄片| 免费看光身美女| 亚洲av国产av综合av卡| 久久久国产一区二区| 日本欧美国产在线视频| 一级爰片在线观看| 色网站视频免费| 国产一区二区三区av在线| 一本—道久久a久久精品蜜桃钙片| 人人妻人人添人人爽欧美一区卜| 午夜激情福利司机影院| 国产精品蜜桃在线观看| kizo精华| 男女边摸边吃奶| 桃花免费在线播放| 久久午夜福利片| 插阴视频在线观看视频| 最近2019中文字幕mv第一页| 伊人久久国产一区二区| 亚洲精品一二三| 中文资源天堂在线| 男女啪啪激烈高潮av片| 男女边摸边吃奶| 久久久久视频综合| 91久久精品国产一区二区成人| 国产精品一二三区在线看| 午夜免费观看性视频| 欧美激情极品国产一区二区三区 | 香蕉精品网在线| 日韩一区二区视频免费看| 夫妻性生交免费视频一级片| 日本av免费视频播放| 男女国产视频网站| 亚洲国产精品成人久久小说| 亚洲av成人精品一区久久| 老司机影院毛片| 五月伊人婷婷丁香| 黄色毛片三级朝国网站 | av免费观看日本| 伦理电影免费视频| 日韩亚洲欧美综合| 少妇精品久久久久久久| 亚洲经典国产精华液单| 欧美亚洲 丝袜 人妻 在线| 中国三级夫妇交换| 日韩强制内射视频| 毛片一级片免费看久久久久| 国产 精品1| av线在线观看网站| 日日摸夜夜添夜夜添av毛片| 我的老师免费观看完整版| 日韩精品有码人妻一区| 国产高清国产精品国产三级| 亚洲欧美清纯卡通| 欧美成人午夜免费资源| 亚洲国产精品999| 我的老师免费观看完整版| 久久久久久伊人网av| av国产精品久久久久影院| 在现免费观看毛片| 赤兔流量卡办理| 亚洲熟女精品中文字幕| 一级毛片我不卡| 99久久精品一区二区三区| av在线播放精品| 久久久久久人妻| 制服丝袜香蕉在线| 亚洲情色 制服丝袜| 高清毛片免费看| 免费观看av网站的网址| 欧美日韩视频高清一区二区三区二| 99热这里只有精品一区| 男女边摸边吃奶| 亚洲欧美日韩卡通动漫| 涩涩av久久男人的天堂| 国产一区二区三区av在线| 久久午夜福利片| 亚洲丝袜综合中文字幕| 成人国产麻豆网| 大香蕉97超碰在线| 精品人妻熟女毛片av久久网站| 国产精品秋霞免费鲁丝片| 一级毛片我不卡| 赤兔流量卡办理| 免费观看a级毛片全部| 性色avwww在线观看| 成人综合一区亚洲| 搡女人真爽免费视频火全软件| 在线 av 中文字幕| 伦理电影大哥的女人| 亚洲自偷自拍三级| 久久久久久人妻| 国产成人精品无人区| 人人妻人人添人人爽欧美一区卜| 涩涩av久久男人的天堂| 免费观看a级毛片全部| av不卡在线播放| 99久久精品一区二区三区| 美女主播在线视频| 免费黄色在线免费观看| 久久久国产精品麻豆| 精品久久久久久久久av| av在线app专区| 日韩av不卡免费在线播放| 一二三四中文在线观看免费高清| 人妻 亚洲 视频| 中文欧美无线码| 精品少妇久久久久久888优播| 黑人猛操日本美女一级片| 国产黄片美女视频| 国内少妇人妻偷人精品xxx网站| 精品熟女少妇av免费看| 久久国产精品大桥未久av | 久久国产亚洲av麻豆专区| 黑人猛操日本美女一级片| 国产精品一区二区在线不卡| 五月玫瑰六月丁香| 99热网站在线观看| 久久99蜜桃精品久久| 亚洲久久久国产精品| 80岁老熟妇乱子伦牲交| 国产探花极品一区二区| 亚洲经典国产精华液单| 看十八女毛片水多多多| 成年人免费黄色播放视频 | 伦理电影大哥的女人| 国产极品天堂在线| 日韩中文字幕视频在线看片| av国产久精品久网站免费入址| 久久精品国产a三级三级三级| 国产欧美亚洲国产| 久久精品国产自在天天线| 午夜免费鲁丝| 女的被弄到高潮叫床怎么办| 亚洲,欧美,日韩| 免费观看a级毛片全部| 亚洲色图综合在线观看| 极品教师在线视频| 黄色配什么色好看| 国产伦精品一区二区三区视频9| 日日爽夜夜爽网站| 成人二区视频| 久久亚洲国产成人精品v| 一区二区三区精品91| 国产综合精华液| av在线老鸭窝| 国产极品粉嫩免费观看在线 | 丝袜脚勾引网站| 在线播放无遮挡| 18禁裸乳无遮挡动漫免费视频| 亚洲av综合色区一区| 天堂中文最新版在线下载| 免费观看性生交大片5| 男人添女人高潮全过程视频| 最近最新中文字幕免费大全7| 亚洲成人一二三区av| 亚洲欧美精品自产自拍| 欧美日韩国产mv在线观看视频| 国产伦理片在线播放av一区| 观看av在线不卡| 国产成人a∨麻豆精品| 国产精品麻豆人妻色哟哟久久| 交换朋友夫妻互换小说| 久久ye,这里只有精品| 国产精品久久久久久精品古装| 看免费成人av毛片| 国产真实伦视频高清在线观看| 国产精品三级大全| 最近2019中文字幕mv第一页| 国产又色又爽无遮挡免| 国产黄频视频在线观看| 亚洲av电影在线观看一区二区三区| 色视频www国产| 国产在线视频一区二区| 日韩伦理黄色片| av卡一久久| 亚洲av不卡在线观看| 在线天堂最新版资源| 大码成人一级视频| 中文乱码字字幕精品一区二区三区| 精品一区二区三卡| 欧美精品一区二区大全| √禁漫天堂资源中文www| 一级黄片播放器| 亚洲国产精品999| 观看av在线不卡| 日日爽夜夜爽网站| 又爽又黄a免费视频| 日本-黄色视频高清免费观看| 国产国拍精品亚洲av在线观看| 99热6这里只有精品| 伊人久久精品亚洲午夜| 搡女人真爽免费视频火全软件| 能在线免费看毛片的网站| 欧美精品高潮呻吟av久久| 国产精品一区二区性色av| 欧美国产精品一级二级三级 | 国产在线一区二区三区精| 2022亚洲国产成人精品| 人妻少妇偷人精品九色| 中文天堂在线官网| 国产精品嫩草影院av在线观看| 女人久久www免费人成看片| 亚洲国产精品一区三区| 国产成人午夜福利电影在线观看| 一级a做视频免费观看| 久久精品国产鲁丝片午夜精品| 成人黄色视频免费在线看| 热99国产精品久久久久久7| 久热这里只有精品99| 国产黄色免费在线视频| 嫩草影院入口| av视频免费观看在线观看| 欧美bdsm另类| 麻豆成人午夜福利视频| a级片在线免费高清观看视频| 日韩中文字幕视频在线看片| 亚洲国产毛片av蜜桃av| 美女国产视频在线观看| 97在线人人人人妻| 极品人妻少妇av视频| 久热久热在线精品观看| 国产精品无大码| 亚洲精品第二区| 两个人的视频大全免费| 在现免费观看毛片| 国产熟女欧美一区二区| 最近手机中文字幕大全| 人体艺术视频欧美日本| 丰满少妇做爰视频| 老司机影院毛片| 91久久精品国产一区二区成人| 在线精品无人区一区二区三| 9色porny在线观看| 亚洲国产精品999| 丰满乱子伦码专区| 成人午夜精彩视频在线观看| xxx大片免费视频| 老女人水多毛片| 少妇人妻 视频| 亚洲av不卡在线观看| 国产一区二区三区av在线| 国产日韩欧美视频二区| 22中文网久久字幕| 97精品久久久久久久久久精品| 中文字幕人妻丝袜制服| 我要看黄色一级片免费的| 亚洲经典国产精华液单| 精品卡一卡二卡四卡免费| 免费大片黄手机在线观看| 国产精品女同一区二区软件| 亚洲久久久国产精品| 成人漫画全彩无遮挡| 成人午夜精彩视频在线观看| 91久久精品国产一区二区三区| 久久婷婷青草| 久久久精品94久久精品| 国产午夜精品一二区理论片| 欧美精品国产亚洲| 国产一级毛片在线| freevideosex欧美| 性色av一级| 日韩一本色道免费dvd| 下体分泌物呈黄色| 国产成人精品无人区| 国产精品女同一区二区软件| 久热这里只有精品99| 少妇的逼好多水| 美女视频免费永久观看网站| 中文天堂在线官网| 丝袜脚勾引网站| 午夜福利网站1000一区二区三区| 欧美日韩av久久| 在线看a的网站| 亚洲四区av| 免费观看av网站的网址| 国产老妇伦熟女老妇高清| 高清午夜精品一区二区三区| 一本一本综合久久| 久久久久久久久久久久大奶| 亚洲欧美一区二区三区国产| 蜜臀久久99精品久久宅男| 亚洲精品第二区| 91午夜精品亚洲一区二区三区| 久久久久精品性色| 午夜日本视频在线| 久久精品久久久久久久性| 自拍偷自拍亚洲精品老妇| 免费黄色在线免费观看| 国产精品福利在线免费观看| 国产精品久久久久久精品古装| 国产一区二区三区av在线| 亚洲av综合色区一区| 如日韩欧美国产精品一区二区三区 | 亚洲怡红院男人天堂| 久久精品国产亚洲av天美| 日韩精品有码人妻一区| 涩涩av久久男人的天堂| 777米奇影视久久| 国产日韩欧美在线精品| 男女边吃奶边做爰视频| 18+在线观看网站| 91久久精品国产一区二区三区| 另类亚洲欧美激情| 丰满迷人的少妇在线观看| 久久av网站| 青青草视频在线视频观看| 在线精品无人区一区二区三| 午夜福利视频精品| 九九在线视频观看精品| 精品一区二区三区视频在线| 亚洲精品亚洲一区二区| 亚洲成人一二三区av| 免费观看a级毛片全部| 女的被弄到高潮叫床怎么办| 菩萨蛮人人尽说江南好唐韦庄| 亚洲av在线观看美女高潮| 亚洲av.av天堂| 国产亚洲av片在线观看秒播厂| 免费不卡的大黄色大毛片视频在线观看| 国产亚洲av片在线观看秒播厂| 爱豆传媒免费全集在线观看| 国产欧美日韩精品一区二区| 久久ye,这里只有精品| 国产伦在线观看视频一区| 少妇裸体淫交视频免费看高清| 99国产精品免费福利视频| av在线播放精品| 成人亚洲精品一区在线观看| 久久国内精品自在自线图片| 日韩欧美 国产精品| 简卡轻食公司| 日日啪夜夜撸| 中文字幕亚洲精品专区| 色婷婷av一区二区三区视频| 在线观看一区二区三区激情| 一级a做视频免费观看| 久久久欧美国产精品| 久久人人爽av亚洲精品天堂| 性高湖久久久久久久久免费观看| 一本大道久久a久久精品| 我的老师免费观看完整版| 最近2019中文字幕mv第一页| 少妇人妻一区二区三区视频| 99久久人妻综合| 亚洲精品乱久久久久久| 午夜视频国产福利| 亚洲精品久久午夜乱码| 人人妻人人澡人人爽人人夜夜| 午夜老司机福利剧场| 久久人妻熟女aⅴ| 国产老妇伦熟女老妇高清| 亚洲精品视频女| 免费看av在线观看网站| 九九久久精品国产亚洲av麻豆| 在线观看人妻少妇| 黑人猛操日本美女一级片| 国产毛片在线视频| 91精品伊人久久大香线蕉| 国产熟女欧美一区二区| 最近2019中文字幕mv第一页| av免费观看日本| 伊人久久国产一区二区| 国产精品免费大片| 亚洲久久久国产精品| av卡一久久| 永久网站在线| 国产乱来视频区| 欧美bdsm另类| 国产日韩欧美在线精品| 亚洲av电影在线观看一区二区三区| 亚洲欧美日韩东京热| 美女内射精品一级片tv| 亚洲av日韩在线播放| 麻豆乱淫一区二区| 成人国产麻豆网| 亚洲av成人精品一区久久| a级毛片免费高清观看在线播放| 高清视频免费观看一区二区| 男男h啪啪无遮挡| 亚洲欧美一区二区三区黑人 | 亚洲欧美成人综合另类久久久| 久久久国产欧美日韩av| 91久久精品国产一区二区三区| 国产在线视频一区二区| 春色校园在线视频观看| 一级,二级,三级黄色视频| 亚洲精品456在线播放app| 美女国产视频在线观看| 乱系列少妇在线播放| 男的添女的下面高潮视频| 18禁在线无遮挡免费观看视频| 精品少妇黑人巨大在线播放| 高清av免费在线| 亚洲丝袜综合中文字幕| 自拍欧美九色日韩亚洲蝌蚪91 | 国产精品国产三级专区第一集| 成年人午夜在线观看视频| 天堂中文最新版在线下载| 国产综合精华液| 国产精品国产av在线观看| 午夜久久久在线观看| 成人综合一区亚洲| 亚洲精品乱久久久久久| tube8黄色片| 丰满人妻一区二区三区视频av| 国产精品蜜桃在线观看| 一区在线观看完整版| 妹子高潮喷水视频| av播播在线观看一区| 久久精品久久久久久久性| 免费高清在线观看视频在线观看| 麻豆乱淫一区二区| freevideosex欧美| 日韩一区二区三区影片| 免费高清在线观看视频在线观看| 免费黄色在线免费观看| 最新的欧美精品一区二区| a级毛色黄片| 国产精品无大码| 久久久久久久大尺度免费视频| 精华霜和精华液先用哪个| 18禁动态无遮挡网站| av免费观看日本| 亚洲精品国产av蜜桃| 视频区图区小说| 亚洲精品国产av蜜桃| 免费播放大片免费观看视频在线观看| 成年av动漫网址| av福利片在线| 熟妇人妻不卡中文字幕| 亚洲欧洲国产日韩| 丝瓜视频免费看黄片| 亚洲国产精品专区欧美| 91成人精品电影| 亚洲自偷自拍三级| 日韩精品免费视频一区二区三区 | 97在线人人人人妻| 一级毛片 在线播放| 久久久国产欧美日韩av| 777米奇影视久久| 欧美激情国产日韩精品一区| 最近手机中文字幕大全| 只有这里有精品99| 国产乱人偷精品视频| 精品一区在线观看国产| 看十八女毛片水多多多| 欧美精品亚洲一区二区| 精品国产乱码久久久久久小说| 日本欧美视频一区| 国产综合精华液| 天堂中文最新版在线下载| 超碰97精品在线观看| 黑人高潮一二区| 99国产精品免费福利视频| 2021少妇久久久久久久久久久| 中文欧美无线码| 国产一区有黄有色的免费视频| 青青草视频在线视频观看| 免费观看av网站的网址| 亚洲av综合色区一区| 国产一区亚洲一区在线观看| 人妻一区二区av| 亚洲精品国产av成人精品| 街头女战士在线观看网站| 国产精品麻豆人妻色哟哟久久| 视频中文字幕在线观看| 久久人妻熟女aⅴ| 国产一级毛片在线| 又粗又硬又长又爽又黄的视频| 国产成人freesex在线| 老司机影院成人| 最黄视频免费看| 男人添女人高潮全过程视频| 日日摸夜夜添夜夜爱| a级毛片在线看网站| 99久久人妻综合| 我要看日韩黄色一级片| 国产有黄有色有爽视频| 婷婷色av中文字幕| 国产成人精品婷婷| 亚洲中文av在线| 欧美精品亚洲一区二区| 亚洲伊人久久精品综合| 国产伦理片在线播放av一区| 欧美 亚洲 国产 日韩一| 五月伊人婷婷丁香| 国产伦精品一区二区三区四那| 18禁在线无遮挡免费观看视频| 日产精品乱码卡一卡2卡三| 国产色爽女视频免费观看| 视频中文字幕在线观看| 日韩欧美 国产精品| 久久精品国产亚洲av天美| 在线观看国产h片| 精品熟女少妇av免费看| 特大巨黑吊av在线直播| 97在线视频观看| 99视频精品全部免费 在线| 男人狂女人下面高潮的视频| 一级片'在线观看视频| 成人免费观看视频高清| 久久精品国产亚洲网站| 国产一区亚洲一区在线观看| 黄色一级大片看看| 亚洲精品成人av观看孕妇| 色吧在线观看| 日韩,欧美,国产一区二区三区| 亚洲国产色片| 日韩不卡一区二区三区视频在线| 亚洲精华国产精华液的使用体验| av国产久精品久网站免费入址| 欧美精品一区二区免费开放| 久久精品熟女亚洲av麻豆精品| 国产免费一区二区三区四区乱码| 亚洲精品,欧美精品| 亚洲三级黄色毛片| 国产免费一区二区三区四区乱码| 欧美日韩在线观看h| 精品一品国产午夜福利视频| 美女内射精品一级片tv| 91aial.com中文字幕在线观看| 精品视频人人做人人爽| 午夜91福利影院| 如日韩欧美国产精品一区二区三区 | 女的被弄到高潮叫床怎么办| 丰满乱子伦码专区| 日本-黄色视频高清免费观看| 熟妇人妻不卡中文字幕| 9色porny在线观看| 美女大奶头黄色视频| 乱码一卡2卡4卡精品| 精品少妇久久久久久888优播| 男女边摸边吃奶|