全尺寸干线输气管道爆炸地振动的空间分布和时频分布特性*

2019-06-05马华原李兴华钟明寿吴建源

爆炸与冲击 2019年4期

马华原，龙源，李兴华，钟明寿，周辉，吴建源

（1.陆军工程大学，江苏南京 210007；2.白城兵器试验中心，吉林白城 137000；3.中国人民解放军75833部队，广州广东 510000）

随着世界能源趋向低碳化、无碳化，天然气作为一种高效、清洁的能源已备受青睐。目前，天然气在世界能源结构中已跃居第二位，预计将会超过石油而居第一[1]。尽管高压天然气管道具有高效、经济等特点，但由于近域施工行为、环境腐蚀、机械或材料失效、自然灾害及其他未知因素的影响，天然气管道爆裂事故时有发生（见图1）[2-4]。

图1 天然气管道爆炸事故Fig.1 Pipeline gas explosion accidents

目前，关于天然气管道爆炸危害效应的研究主要以热辐射和火焰危害为主[5-7]。关于振动的研究大多停留在数值模拟阶段：Su[8]通过建立土层与管道耦合效应的非线性数学模型，对地下管道在复合振动加载模式下的动态响应进行研究，分析了振动对管道的影响，得到评估埋地管道抗震性能的方法；Mahdavi等[9]研究了地下管线在黏性土中的荷载传递特性，分析了管道回填土和沟道宽度的几何效应，得到了土壤强度分布和响应因子。目前尚缺乏对全尺寸干线输气管道爆炸地震波传播规律和时频分布特性的研究，本文中针对高压天然气管道爆炸瞬间物理现象及其产生的振动危害效应展开研究，分析地震波的频谱特性、衰减规律及时频特性。研究成果不仅对并行管道安全设计具有理论指导意义和工程应用价值，同时对管线附近建（构）筑物和重要设施的抗震与防护设计也具有参考作用。

1 全尺寸富气管道爆炸试验

试验选用第三代大输量天然气管道，该型管道大量应用于西气东输二线工程。管身材料采用X90级别管线钢，弹性模量206 GPa，泊松比0.3，屈服强度715 MPa，直径1 219 mm，壁厚16.3 mm，内压12 MPa，埋深1.2 m。管道总长约430 m，中间部分是长130 m的试验段，为尽可能地模拟无限长输气管道爆炸时的泄压状态，在试验段两端各设置长150 m的储气段，如图2所示。管内填充天然气，压力12 MPa，气体摩尔分数分别为：x(CH4)=94.91%，x(C2H6)=2.55%，x(CO2)=0.94%，x(N2)=1.4%，x(其他)=0.2%。在试验段中点沿轴线方向，设置长0.5 m的线型聚能切割器，以引入初始裂纹使管道发生爆裂。

试验测线布设，以起爆点为原点，在与管道轴线分别呈0°、30°、60°和90°的4个方向上设置测线。其中，30°、60°和90°测线上设置6个测点，0°测线上设置5个测点。每个测点处安装振动速度测试仪和振动加速度测试仪各1台。具体布设方案如图3所示，其中60°、30°以及0°测线的测点间距同90°测线。

图2 试验场Fig.2 Test field

图3 测点布设方案Fig.3 Layout scheme of measuring points

试验选用TC-4850和Blast-UM型振动测试仪。两型仪器均配备集成式三通道传感器，可同时记录水平径向（X轴）、水平切向（Y轴）和铅垂方向（Z轴）3个方向的振动数据。试验前，所有仪器均在江苏省计量科学研究院进行了校准，以确保数据的真实性。

起爆瞬间以及爆后管道裂口如图4所示。试验起爆过程为：发射燃烧弹，聚能切割器点火。射流刀穿过管壁形成初始裂缝，管壁在初始裂纹和巨大内压共同作用下裂纹产生扩展。管内高压天然气迅速喷出，气团接触到燃烧弹后被瞬间点燃。起爆后约3 s，气团发生完全爆燃，现场形成巨大蘑菇云，直径达150 m左右。爆炸后无人机传回画面，可看到有残余气体持续燃烧约5 min后熄灭。现场留下长约20 m、宽约5 m的漏斗坑，爆后测得裂口长度为17.43 m。

图4 爆炸瞬间及爆后管道Fig.4 Explosions and pipeline affected

2 天然气管道爆炸地振动的空间分布特性

试验得到了各测点数据详细信息，见表1～4（其中6# 和12# 测点设备未正常触发）。表中，l为距离，VX、VY、VZ为X、Y、Z方向最大振动速度，fX、fY、fZ为X、Y、Z方向主频，V为矢量和。为研究振动能量的衰减规律以及空间分布特性，选取各测点振动信号的矢量和数据进行深入研究。

表1 90° 测线振动速度Table 1 Data of vibration velocity of measuring line 90°

表2 60° 测线振动速度Table 2 Data of vibration velocity of measuring line 60°

表3 30° 测线振动速度Table 3 Data of vibration velocity of measuring line 30°

表4 0° 测线振动速度Table 4 Data of vibration velocity of measuring line 0°

由图5可以看出，振动峰值的衰减呈现先快后慢的模式，为了进一步研究其衰减规律，分别用指数拟合和幂率拟合对各个测线上的数据进行拟合并对比效果，如图6所示。

图5 振动峰值分布图Fig.5 Peak vibration distribution diagram

在文献[10-11]中，爆破地震波的衰减规律遵循幂率分布。然而，由图6的拟合效果对比可以看出，4条测线上数据的指数拟合结果均优于其幂率拟合结果，指数拟合有着更小的拟合误差值（SSE），并且其确定系数（R-square）更接近1。因此，天然气管道爆炸地振动衰减更符合指数分布，本次试验中90°、60°、30°、0°测线上的衰减拟合结果分别为：

式中：V为振动速度，单位为cm/s；r为距爆心距离，单位为m。

从图5～6以及拟合结果可以看出，振动强度和测线角度显著相关。其中，90o测线振动强度最大，之后依次为30°测线、60°测线和0°测线。为了更直观的反应振动强度在地面的分布，首先对采集到的数据在空间内进行两次镜像延拓，然后采用重调和方程对所得数据进行了三维拟合，得到了振动强度在空间内的分布场，效果如图7所示。

图7中，X轴（Y=0）与管道重合。可以看出，地震波在空间内不呈圆对称分布，这一点与普通集团装药爆炸地震波特性不同。管道爆炸地震波在90°方向和30°方向较强，60°方向和0°方向较弱。这样的分布特性意味着，在讨论天然气管道爆炸地震波损害范围时，不能只考虑距离因素，还需考虑角度因素。

图6 拟合效果对比Fig.6 Comparison of fitted effects

3 天然气管道爆炸地振动的时频分布特性

选用周辉等[12]优化的双参数匹配追踪算法联合WVD分布的时频分解算法。其中WVD分布是由Ville引入信号分析领域的。根据特征函数方法，推导信号X(t)的WVD分布为：

图7 振动场三维拟合效果图Fig.7 3D fitting diagram of vibration field

WVD拥有极高的时间-频率分辨率，但由于其自身为二次型时频解析，因此不具备线型可加性，严重影响了它在非平稳随机信号分析领域的运用。为解决此问题，将匹配追踪算法引入，其原理是将多成分混杂的原始信号分解为基础子波，分别对子波进行WVD分解，再将所有子波的WVD结果相加，便可得到原始信号的真实WVD分布。另外，选用的双参数匹配追踪算法是结合了Hilbert变换，对原始信号的优势频率和瞬时相位先行确定，从而大大降低了匹配追踪的运算量，为大数据量处理提供了基础。

将1#、7#、13#和19#测点的振动信号采用改进型MP-WVD算法进行处理，得到其时频分布如图8～11所示。

由图8～11可看出，各个测点振动数据的时频分布呈现一定的规律性：振动能量频率中心都分布在10～20 Hz范围内，振动持续时长约为0.1～0.2 s。同时，X轴振动在时频域内相对集中，而Y轴和Z轴振动在时频域的分布相对比较分散。振动加强区与削弱区的时频特性并没有显著区别。将时频谱对时间轴做积分，得到频率边际谱[13-14]，相比于通过傅里叶变换得到的幅频谱，频率边际谱可以更直观地看到各个频率成分在全局上的累加，如图12所示。

图8 1#测点振动信号时频分布Fig.8 Time frequency distribution of measuring point 1#

图9 7#测点振动信号时频分布Fig.9 Time frequency distribution of measuring point 7#

图10 13#测点振动信号时频分布Fig.10 Time frequency distribution of measuring point 13#

图11 19#测点振动信号时频分布Fig.11 Time frequency distribution of measuring point 19#

图12 频率边际谱Fig.12 Marginal spectrum

由图12可以看出，天然气管道爆炸地震波绝大部分能量分布在0～40 Hz频带，40 Hz以上频带所占比重极小。并且代表瑞利波的X轴、Z轴分量明显强于代表勒夫波的Y轴分量，因此可以判断管道爆炸地震波主要以瑞利波形式传播。通过对振动速度信号的滤波和积分得到其质点运动轨迹如图13所示，可以看出明显的瑞利波“椭圆形”轨迹的特征，并且水平方向位移较大。