大口徑高壓管道物理爆炸沖擊波傳播規律的試驗研究

程良玉， 龍 源， 毛益明， 徐全軍， 紀 沖， 吳建源

(陸軍工程大學 野戰工程學院， 南京 210007)

大口徑高壓管道物理爆炸沖擊波傳播規律的試驗研究

程良玉， 龍 源， 毛益明， 徐全軍， 紀 沖， 吳建源

(陸軍工程大學 野戰工程學院， 南京 210007)

根據動態測試理論，構建基于DH5960的網絡動態信號測試分析系統，開展直徑為1 422 mm的大口徑全尺寸高壓氣體管道爆炸試驗。獲取各測點空氣沖擊波的P(t)曲線，經小波變換濾波后采用最小二乘法對沖擊波在空間上的傳播特性進行數值擬合，得到了衰減曲線和衰減公式。對比試驗組和對照組的沖擊波曲線，發現高壓管道物理爆炸空氣沖擊波在近地面自由場的傳播具有非圓形波陣面的特征。

物理爆炸；小波變換去噪；超壓衰減曲線；非圓形波陣面

天然氣作為一種高效、清潔的能源為各國所青睞。未來10余年內，我國經濟發展和居民對天然氣的需求量仍將呈爆炸式增長。截至2015年底，我國天然氣主干管道總里程達8.5萬 km，已經形成相對完善的天然氣管網。由于天然氣管道受到近域施工行為、環境腐蝕、機械或材料失效、自然災害及其他未知因素的影響，導致開裂爆炸事故時有發生。

天然氣管道爆炸事故有兩種：①物理爆炸，管道破裂后內部高壓氣體解除約束逸出并向自由場迅速膨脹[1]，瞬間壓縮空氣形成沖擊波并向外傳播；②化學爆炸，天然氣與空氣混合后，在一定外界條件作用下發生爆燃，形成化學爆炸事故。其中后者需要可燃氣體在空氣中擴散并受到外力作用，發生爆燃的可能性受引燃條件(如溫度，混合比)的影響。而物理爆炸不需要可燃氣體與氧化性氣體混合，也不需要外界條件引燃。相較于天然氣管道泄漏后發生的化學爆炸，物理爆炸發生的幾率更高。

天然氣管道泄漏后發生物理爆炸的問題一直為國內外一些專家學者所關注。Bariha等[2]理論分析了石油天然氣管道發生泄漏時的氣體沖擊波壓力場的狀態；李偉等[3]采用CFD(Computer Fluid Dynamics)數值模擬計算了高壓氫氣輸運裝置物理爆炸事故后的空氣沖擊波狀態場；Skacel等[4]通過激波管縮比試驗研究了物理爆炸沖擊波在管內傳播規律，并驗證了TNT當量法對于物理爆炸沖擊波一維傳播估算的可行性。這類研究通過不同的方法近似分析了物理爆炸沖擊波的狀態，但對于物理爆炸形成空氣沖擊波在自由場中傳播的天然氣管道全尺寸試驗研究卻鮮見報道。

本文基于大口徑高壓氣體管道爆炸的全尺寸試驗，獲得了物理爆炸產生空氣沖擊波測試的數據。采用小波變換去噪方法，對物理爆炸空氣沖擊波在自由場的超壓衰減規律及其壓力場狀態進行了研究。

1 試驗設計與實施

1.1 高壓氣體管道

為研究高壓天然氣管道爆炸產生的空氣沖擊波的傳播規律，開展了大口徑高壓氣體管道爆炸的全尺寸試驗。試驗的管道采用X80鋼材，管道長20 m，直徑為1 422 mm，如圖1所示。管道淺埋地下，并向管內充入一定量空氣使壓力升高至12 MPa，試驗現場如圖2所示。采用聚能切割器在管道的中心位置瞬間切割開一個50 cm的長條形切縫，管內高壓氣體解除約束瞬間逸出，形成物理爆炸。逸出氣體到達地面上空自由場后向周圍膨脹形成空氣沖擊波向外傳播。

圖1 高壓氣體管道中間段Fig.1 The middle section of high-pressure gas pipe

圖2 高壓氣體管道爆炸試驗現場Fig.2 The site of high-pressure gas pipe burst measuring

1.2 測點布設

為對比研究高壓氣體管道爆炸空氣沖擊波的傳播規律，設置試驗組和對照組兩條沖擊波壓力測試線，分別為垂直于管道方向的垂直測試線和與管道成45°傾斜角的傾斜測試線。兩條測試線各布設10個測點，共20個自由場壓力測點。垂直測試線過起爆點并垂直于管道軸線，各測點距起爆點距離如表1所示；斜測試線過起爆點并沿管道東南方向，與管道軸線成45°夾角，各測點距起爆點距離如表2所示。各測點的位置示意圖如圖3所示。

為將壓力傳感器固定于近地面自由場中，避免受沖擊波壓力產生運動而影響測試，設計了自由場壓力傳感器支架。將桿式自由場傳感器固定于各個對應的測點，使其距離地面1.9 m[5]。

1.3 測試系統

考慮到試驗環境的特殊性，試驗構建了基于DH5960的網絡型動態信號測試分析系統，如圖4所示。該系統采用集成的CYG411型壓阻式壓力傳感器提高物理瞬態響應的靈敏度，并通過計算機對采集的數據進行實時存儲，具有瞬態采樣頻率高，頻率響應范圍大等特點。為操作安全，設置兩個采集點進行數據存儲與收集，一個位于爆炸場內的測試間對數據進行直接采集，另一個位于安全范圍的控制端，二者通過光纖網絡聯接。試驗人員通過控制端計算機操控測試間計算機采集與傳輸試驗數據，具有安全性高的特點。

表1 垂直測試線上各測點距起爆點距離

表2 斜測試線上各測點距起爆點距離

圖3 測點布設示意圖(單位：m)Fig.3 Schematic diagram layout of measuring point (unit: m)

圖4 DH5960動態信號測試分析系統Fig.4 DH5960 dynamic signal measurement and analysis system

2 試驗結果及數據處理分析

試驗共采集到20組沖擊波信號。研究發現原始P(t)曲線有很多“毛刺”(噪聲信號)，嚴重干擾有效信號的讀取。通過分析發現這些噪聲信號主要有兩個來源：①大幅度擾動，頻率范圍在傳感器支架固有頻率附近，根據速度規律，推測為管道爆炸在土壤介質中形成的應力波擾動[6-7]。該擾動通過大地傳至傳感器支架底座，再通過傳感器支架對傳感器形成噪聲干擾。這部分噪聲可以通過測點布設規避，使之不出現在沖擊波信號到達的時間段內，也可以通過濾波濾去。②小幅度擾動，頻率范圍在傳感器固有頻率附近，以聲速在各個測點間傳播，推測為沖擊波形成與傳播過程中對周圍空氣引起的擾動。該擾動通過空氣介質直接對傳感器形成干擾。這部分噪聲無法規避，但是可以通過濾波除去。

采用小波變換對試驗數據進行濾波處理[8-9]，以除去“毛刺”，獲得清晰的沖擊波信號，增強數據的可讀性。

傳感器對沖擊波壓力的物理響應為連續信號，但數字采集儀只能采集離散的數字信號，因此采用離散小波變換(Discrete Wavelet Transform，DWT)對原始信號進行處理。離散小波變換可通過離散化連續小波變換(Continuous Wavelet Transform，CWT)中的伸縮因子和平移因子得到[10]。離散小波為

(1)

相應的離散小波變換為

(2)

式中：ψ(t)∈L2(R);a為離散化的伸縮因子，b為離散化的平移因子，a,b∈R；m,n為離散因子，m,n∈R。

通過檢驗殘差的收斂性，驗證去噪信號與原始沖擊波信號的相關性。原始P(t)曲線，去噪后的P(t)曲線和殘差曲線對比如圖5和圖6所示。圖中第一層為原始信號，第二層為去噪信號，第三層為殘差曲線。

圖5 垂直測線部分測點的原始曲線、濾波曲線和殘差曲線Fig.5 Original signal、de-noised signal and residuals of points in vertical measuring line

圖6 傾斜測線部分測點的原始曲線、濾波曲線和殘差曲線Fig.6 Original signal、de-noised signal and residuals of points in bevel measuring line

從P(t)曲線可以看出，當沖擊波陣面到達測點時，壓力突躍上升形成強間斷面。隨著氣體向外膨脹壓力迅速衰減，衰減至大氣壓力時，由于慣性作用繼續下降形成一個壓力稍小于大氣壓力的負壓區。負壓區被周圍高壓區壓縮，壓力升高至稍高于大氣壓力。接著又開始第二次膨脹和壓縮的脈動過程。通過多次有阻尼脈動，壓力衰減至大氣壓力；并且同一條測線上距離越遠的測點所測到的超壓峰值越小，衰減速度越慢。這與常規爆炸沖擊波的傳播規律類似。

從殘差曲線圖中可以看出，時域范圍內，殘差一直在零點附近波動，其主值基本收斂于直線P=0，說明小波變換去噪沒有改變原始信號的沖擊突躍部分，去噪結果清晰可靠。將讀取的沖擊波數據列于表3和表4(其中B05#測點未測到有效信息)。

表3 垂直測試線上各測點超壓峰值

表4 斜測試線上各測點超壓峰值

3 沖擊波傳播特性分析

3.1 超壓衰減曲線擬合

空氣沖擊波在自由場中傳播具有氣體壓力大，傳播速度快等特點，此過程中氣體的重力和黏性可忽略不計，氣體可以假定為完全氣體[11]，根據相似理論有

Δpm=f(E,r,P0,ρ0,k)

(3)

式中：被定量ΔPm為沖擊波超壓；主定量E為爆炸總能量；r為沖擊波波陣面距爆心距離；P0為空氣的初始壓力；ρ0為空氣的初始密度；k為絕熱指數。

由π定理寫出三個π項

(4)

根據量綱和諧條件求解得到

(5)

k一般取常數，P0一定的條件下，式(5)可寫成

(6)

將式(6)展開

(7)

爆炸能量由管內氣體唯一確定，沖擊波傳播過程中E不發生改變，可當作常量，則式(7)變為

(8)

已知邊界條件：當r→∞時，ΔPm=0，故A0=0，則

(9)

結合試驗數據，采用最小二乘法分別對兩條測線上的超壓進行曲線擬合，獲得擬合曲線和擬合函數，擬合曲線如圖7和圖8所示。

圖7 Vertical測線的擬合曲線Fig.7 Fitting curve of vertical measuring line

圖8 斜測線的擬合曲線Fig.8 Fitting curve of bevel measuring line

垂直方向的超壓衰減擬合函數

(10)

45°角傾斜方向的超壓衰減函數

(11)

3.2 沖擊波陣面的非圓形特征

研究試驗組與對照組的沖擊波超壓發現，相同距離時垂直測線上測點超壓明顯高于傾斜測線上對應測點的超壓，即截取相同超壓時，垂直測線上的距離比斜測線上的距離更遠。說明在近地面自由場中空氣沖擊波超壓的等壓線為非圓形，如圖9所示。沖擊波陣面為壓力最大的等壓線，在近點面自由場中的分布也具有非圓形的特征。

圖9 沖擊波等壓線Fig.9 Shockwave isobars schematic diagram

分析發現，非圓形波陣面的形成與管道破裂氣體逸出及沖擊波傳播過程有關。如圖10所示，切割器沿管道軸線切開一個長條形切縫。切縫軸向尺寸大于徑向尺寸，即徑向自由面大于軸向自由面，致使其徑向強度小于軸向強度。在內部高壓作用下，切縫先從中間位置沿徑向發生塑性變形，撕裂成裂口如圖11所示。由于徑向塑性變形程度更大，氣體從裂口逸出時，徑向泄漏更快，更早形成沖擊波，導致初始波陣面為非圓形波陣面。

各方面非均勻的空氣沖擊波在傳播過程中有一個各方面漸近均勻化作用[12]，將非圓形波陣面均勻化為近似圓形的波陣面。化學爆炸形成的空氣沖擊波超壓一般很大，因此各方面漸近均勻化作用很明顯，很快形成近似圓形的沖擊波陣面；而管道爆炸空氣沖擊波超壓較小，各方面漸近均勻化作用不明顯，因此在近地面自由場中具有明顯的非圓形波陣面特征。

圖10 爆炸切縫Fig.10 Blasting incision

圖11 裂口變形及氣體噴流示意圖Fig.11 Crack deformation and gas jet schematic diagram

4 結 論

(1)沖擊波測試過程中干擾信號主要有兩部分：①管道爆炸在土壤介質中形成的應力波擾動；②沖擊波在周圍空氣中形成的不規則擾動，這兩部分噪聲都與傳感器和支架的固有頻率有關。

(2)根據試驗數據，擬合得到本試驗物理爆炸空氣沖擊波在與管道垂直方向和與管道成45°角方向上的超壓衰減曲線及其函數。可以為實際應用提供有效參考。

(3)高壓氣體管道物理爆炸形成空氣沖擊波，其波陣面形狀受管道破裂變形過程和沖擊波傳播過程雙重影響而具有非圓形特征，如本次試驗所得的沖擊波傳播呈現與管道垂直方向最強，其他方向稍弱的非圓形波陣面特征。說明高壓管道爆炸形成空氣沖擊波的安全防護不同方向具有不同的防護等級。

[ 1 ] STAWCZYK J. Experimental evaluation of LPG tank explosion hazards[J]. Journal of Hazardous Materials，2003，96 (2/3)：189-200.

[ 2 ] BARIHA N，MISHRA I M，SRIVASTAVA V C. Hazard analysis of failure of natural gas and petroleum gas pipelines[J]. Journal of Loss Prevention in the Process Industries ，2016，40： 217-226.

[ 3 ] 李偉，張奇. 高壓氫氣輸運裝置物理爆炸狀態場特征及災害效應研究[J]. 高壓物理學報，2009，23(3)：203-208.LI Wei，ZHANG Qi. Study on the physical explosion character and ejection effect of the high-pressurized hydrogen transport device[J]. Chinese Journal of High Pressure Physics，2009，23(3)：203-208.

[ 4 ] SKACEL R，JANOVSKY B，DOSTAL L，et al. Small-scale physical explosions in shock tubes in comparison with condensed high explosive detonations[J]. Journal of Loss Prevention in the Process Industries，2013，26(6)：1590-1596.

[ 5 ] 杜紅棉，祖靜，馬鐵華，等. 自由場傳感器外形結構對沖擊波測試的影響研究[J]. 振動與沖擊，2011,30(11)：85-89.

DU Hongmian，ZU Jing，MA Tiehua，et al. Effect of mount configuration of free-field transducers on shock wave measurement[J]. Journal of Vibration and Shock，2011,30(11)：85-89.

[ 6 ] 張志鎮，高峰，林斌，等. 巖石沖擊傾向與其波速變化的相關性研究[J]. 巖石力學與工程學報，2012，31(2)：3527-3532.

ZHANG Zhizhen，GAO Feng，LIN Bin，et al. Correlation study of rock burst tendency and ultrasonic velocity[J]. Chinese Journal of Rock Mechanics and Engineering，2012，31(2): 3527-3532.

[ 7 ] 邱志剛，薄景山，羅奇峰. 土壤剪切波速與埋深關系的統計分析[J]. 世界地震工程，2011，27(3): 81-88.

QIU Zhigang，BO Jingshan，LUO Qifeng. Statistical analysis of relationship between shear wave velocity and depth of soil[J]. World Earthquake Engineering，2011，27(3): 81-88.

[ 8 ] 賈云飛,張春棋,康金. 爆炸沖擊波信號測試與小波分析[J]. 彈箭與制導學報，2014，34(5): 86-89.

JIA Yunfei，ZHANG Chunqi，KANG Jin. Shock wave signal measuring and wavelet analysis [J]. Journal of Projectiles，Rockets，Missiles and Guidance, 2014, 34(5): 86-89.

[ 9 ] 張衍芳. 沖擊波信號處理方法的研究[D].太原：中北大學,2011.

[10] 張德豐. MATLAB小波分析[M].北京：機械工業出版社，2009.

[11] 鄭哲敏，談慶明，王補宣. 相似理論與模化[C]∥ 鄭哲敏文集. 北京：科學出版社，2004: 777-808.

[12] 奧爾連科. 爆炸物理學[M]. 孫承緯, 譯. 北京：科學出版社，2011.

Anexperimentalstudyonlarge-diameterhigh-pressurepipelinephysicalexplosionsshockwavepropagation

CHENG Liangyu, LONG Yuan，MAO Yiming，XU Quanjun，JI Chong，WU Jianyuan

(College of Field Engineering，Army Engineering University, Nanjing 210007, China)

According to the theory of dynamic test, this paper presented a DH5960-based dynamic signal analysis system, in which a diameter of 1 422 mm full-size high-pressure gas pipeline explosion experiments was performed. The air blast curve of each measuring point was obtained. After the wavelet de-noising, the least squares method was adopted to carry out numerical fitting of the peak overpressure shock wave propagation characteristics in space. Decay curves and attenuation formula were obtained. Shock curve comparison between experimental group and the control group were performed and discovered that high-pressure pipeline physics explosion air shock has the feature of a noncircular wave front propagation characteristics on the free-field near the ground.

physical explosions; wavelet de-noising; overpressure decay curve; noncircular wave front

國家自然科學基金(11672331)；中國博士后科學基金(2015M58279)

2016-06-24 修改稿收到日期： 2016-07-22

程良玉 男，碩士生，1992年生

龍源 男，博士，教授，博士生導師，1958年生

O383； X937

A

10.13465/j.cnki.jvs.2017.22.007