示：
[0062]
[0063] 在天然气随风速风向稳定扩散阶段，扩散范围方向基本与风速风向一致，可近似 认为天然气浓度衰减与风速风向扩散面积成正比。
[0064] 考虑到城市中实际天然气管道泄漏后，地表扩散阶段会受到楼宇，树木等障碍物 阻碍而衰减扩散速度和减小扩散极限范围。本发明经过查找大量数据并理论验证，采用阻 碍系数F = O.27修正扩散过程。天然气扩散范围随时间变化的曲线如图3所示。具体的仿真 参数如图表1所示：
[0065]
[0066] 表1
[0067] 进一步来讲，在本发明实施例中，通过如下过程来建立预测模型：
[0068] 1、建模样本数据及参数的分类处理：
[0069] 在建立预测模型之前首先对建模样本进行两类分类处理，划分为受爆管高压喷射 动能影响较大的初期高速扩散和受风力风向影响较大的中后期持续扩散模型。
[0070] 2、理化参数的确定：
[0071] 对于建模输入变量的理化参数包括:天然气管内气压、管内气体流速、天然气管泄 漏孔与管线横截面积比、地表高空风速、风向、阻力扩散系数、管内天然气浓度、泄漏系数和 时间等。
[0072]所述各理化参数数据收集的途径主要有三种：实验测定，相关手册及文献已知数 据的筛选。
[0073] 3、预测模型的建立：
[0074] 根据分类处理结果，对两类样本分别建立预测模型，具体步骤如下:将样本集随机 划分为两个部分，其中一部分用于建立模型，另一部分不参与建模，用于对所建模型进行验 证。
[0075] 针对训练集样本，以2所选用的理化参数及气体组分配比作为输入变量，样本对应 的泄漏扩散范围作为输出变量，应用气体流出量及高斯烟团多阶段动态扩散积分算法，对 两者之间的内在定量关系进行统计建模，并用matlab进行仿真，得到两者之间存在的定量 函数关系，建立相应的预测模型。
[0076] 4、模型的验证与修订：
[0077] 预测模型的稳定性的验证方式为交叉验证方法，预测模型的外推能力利用利用外 部验证方式采用建立的模型测验测试集样本的扩散范围。
[0078] 根据上述预测结果，删除误差过大样本，然后返回数据修正模型，直到模型测试结 果在误差允许范围内。
[0079] 5、预测模型的应用：
[0080] 针对样本类别，根据2选择相应的理化参数，计算出相应的数值，将其带入3所建立 的相应预测模型进行计算即可得到天然气泄漏两个阶段的扩散范围数据。天然气泄露扩散 极限距离随时间变化仿真结果如图3所示。
[0081] 对应本发明实施例中一种确定气体扩散范围的方法，本发明实施例中还提供了一 种确定气体扩散范围的系统，如图4所示为本发明实施例中一种确定气体扩散范围的系统 的结构示意图，该系统包括：
[0082] 泄露量检测模块401，用于获取设定时间内的气体总泄露量；
[0083] 范围预测模块402,用于将所述总泄漏量带入指定的高斯烟团多阶段动态扩散模 型，得到气体泄露预测范围；
[0084]处理模块403,用于在所述预测范围中删除误差超过预设阈值的预测范围，得到确 定的预测范围。
[0085]进一步，在本发明实施例中，所述泄露量检测模块401，具体用于根据地下管线参 数以及声速泄露公式，得到泄露开始至阀门关闭的时间段内的第一泄露量;根据阀门关闭 后的泄露时间以及次声速泄露公式，获取所述时间段内的第二泄露量;将所述第一泄露量 与所述第二泄露量之和作为所述总泄露量。
[0086] 阀门关闭前后天然气管道内声速和次声速泄露量可由下式计算得出，其中主要输 入参数为管内压力，管道横截面积，泄露口面积，泄露段阀门封闭区间内管线长度。
[0087] 具体来讲，根据地下管线参数以及声速泄露公式，得到泄露开始至阀门关闭的时 间段内的第一泄露量，计算公式如下：
[0089] 其中，Oms表征第一泄露量，A。!·表征管道横截面积，P2表征管道内压力，M表征管道 长度，RT 2表征泄露口面积。
[0090] 然后，根据阀门关闭后的泄露时间以及次声速泄露公式，获取所述时间段内的第 二泄露量，具体公式如下：
[0092]其中/?%表征第二泄露量，A。!·表征管道横截面积，P2表征管道内压力，M表征管道 长度，RT2表征泄露口面积。
[0094]对上式声速泄露公式可根据泄露开始至阀门关闭时间计算出一定时间段内泄露 量，对上式次声速泄露公式可根据阀门关闭后泄露时间或管内气压衰减曲线运用积分运算 法计算出一定时间内泄露量，两者之和即为阀门关闭后的总泄露量。天然气泄露流量随时 间变换曲线如图2所示。
[0095]在具有管线内压力自动监控的系统中，一般在中等泄露孔管道泄露2分钟后，压力 传感器可检查到管线内气压明显下降并自动关闭阀门或人工关闭。由上式稳态泄漏2分钟 内流量输入计算得出为1050kg/(m 2 · s )，关闭阀门后，根据质量流量随时间变化拟合出简 化快速计算公式G = 1050 XeH12t13
[0096]由于实际天然气管线泄露质量与理论值有一定偏差，可用泄露系数进一步精确结 果。在快速简化计算时泄漏系数Q取0.69左右。
[0097]进一步，在本发明实施例中，所述范围预测模块402,具体用于调取计算得出泄露 系数;将所述泄露系数与所述第一泄露量以及所述第二泄露量之和的乘积作为所述总泄露 量。
[0098]尽管已描述了本申请的优选实施例，但本领域内的普通技术人员一旦得知了基本 创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包 括优选实施例以及落入本申请范围的所有变更和修改。
[0099]显然，本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精 神和范围。这样，倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围 之内，则本申请也意图包含这些改动和变型在内。
【主权项】
1. 一种确定气体扩散范围的方法，其特征在于，所述方法包括： 获取设定时间内的气体总泄露量； 将所述总泄漏量带入指定的高斯烟团多阶段动态扩散模型，得到气体泄露预测范围； 在所述预测范围中删除误差超过预设阔值的预测范围，得到确定的预测范围。2. 如权利要求1所述的方法，其特征在于，获取设定过时间内的其他的总泄露量，包括： 根据地下管线参数W及声速泄露公式，得到泄露开始至阀口关闭的时间段内的第一泄 路里； 根据阀口关闭后的泄露时间W及次声速泄露公式，获取所述时间段内的第二泄露量； 将所述第一泄露量与所述第二泄露量之和作为所述总泄露量。3. 如权利要求2所述的方法，其特征在于，所述声速泄露公式如下：其中，钱,1。表征第一泄露量，Am表征管道横截面积，P2表征管道内压力，Μ表征管道长 度，RT2表征泄露口面积。4. 如权利要求2所述的方法，其特征在于，所述次声速泄露公式如下：5. 如权利要求2所述的方法，其特征在于，将所述第一泄露量与所述第二泄露量之和作 为所述总泄露量，包括： 调取计算得出泄露系数； 将所述泄露系数与所述第一泄露量W及所述第二泄露量之和的乘积作为所述总泄露 量。6. 如权利要求1所述的方法，其特征在于，将所述总泄露量带入指定的高斯烟团多阶段 动态扩散模型，得到气体泄露预测范围，包括： 将所述总泄露量带入如下的高斯烟团多阶段动态扩散模型得到气体泄露预测范围x，y，z，t-t'表征时间段，Q表征总泄露量。7. -种确定气体扩散范围的系统，其特征在于，包括： 泄露量检测模块，用于获取设定时间内的气体总泄露量； 范围预测模块，用于将所述总泄漏量带入指定的高斯烟团多阶段动态扩散模型，得到 气体泄露预测范围； 处理模块，用于在所述预测范围中删除误差超过预设阔值的预测范围，得到确定的预 测范围。8. 如权利要求7所述的系统，其特征在于，所述泄露量检测模块，具体用于根据地下管 线参数W及声速泄露公式，得到泄露开始至阀口关闭的时间段内的第一泄露量;根据阀口 关闭后的泄露时间W及次声速泄露公式，获取所述时间段内的第二泄露量;将所述第一泄 露量与所述第二泄露量之和作为所述总泄露量。9. 如权利要求8所述的系统，其特征在于，所述范围预测模块，具体用于调取计算得出 泄露系数;将所述泄露系数与所述第一泄露量W及所述第二泄露量之和的乘积作为所述总 泄露量。
【专利摘要】本发明公开了一种确定气体扩散范围的方法及系统，采用气体流出量和高斯烟团多阶段动态扩散积分算法并通过现有手册、数据库或实验得到的天然气管道泄漏相关的常规理化多参数的综合性数据分析来确定天然气烟团扩散范围，解决了在数据不全状态下对天然气泄露扩散范围进行较准确预测的难题，提高国家应对突发事故的应急控制能力，对有效避免重大经济损失和人员伤亡具有重要意义。
【IPC分类】G06F19/00
【公开号】CN105550505
【申请号】CN201510902073
【发明人】程揭章, 葛宝玉, 侯升飞, 公秀田, 王建勋, 逄崇雁
【申请人】哈尔滨航天恒星数据系统科技有限公司
【公开日】2016年5月4日
【申请日】2015年12月8日