一种综合管廊燃气泄漏浓度场预测矫正及泄漏率估计方法技术

本发明专利技术公布一种综合管廊燃气泄漏浓度场预测矫正及泄漏率估计方法。该方法结合综合管廊舱内燃气传感器、风速风向传感器及温湿度传感器实时监测数据，通过构建综合管廊舱内燃气泄漏扩散的数值计算模型，运用数据同化算法进行燃气泄漏扩散过程及燃气管道泄漏率的预测估计。所述数据同化算法为改进的集合卡尔曼滤波算法，该算法的状态向量由燃气浓度值和泄漏率组成。本发明专利技术方法结合风速风向传感器和温湿度传感器的监测数据，对风流速度与扩散系数加入符合实际情况的噪声后，可抑制风流速度和扩散系数不确定性对结果的影响，在获得管廊内部燃气浓度分布情况的同时，能准确对燃气管道泄漏率进行反演估计，更加符合实际状况。

Prediction correction and leakage rate estimation method for gas leakage concentration field in integrated Gallery

The invention discloses a method for predicting correction and leakage rate of gas leakage concentration field in comprehensive pipe gallery. The method combines the real-time monitoring data of gas sensor, wind speed and wind sensor and temperature and humidity sensor in the comprehensive Gallery cabin. By constructing the numerical model of gas leakage diffusion in the comprehensive Gallery cabin, the data assimilation algorithm is used to predict the gas leakage and diffusion process and the prediction of gas pipeline leakage rate. The data assimilation algorithm is an improved ensemble Calman filtering algorithm. The state vector of the algorithm is composed of gas concentration and leakage rate. The method combined with the monitoring data of wind speed and wind direction sensor and temperature and humidity sensor, and can suppress the influence of wind velocity and diffusivity uncertainty on the result after adding the noise of wind velocity and diffusion coefficient to the actual situation. It can be accurate to gas pipe while obtaining the distribution of gas concentration in the tube corridor. The estimation of leakage rate is more consistent with the actual situation.

【技术实现步骤摘要】
一种综合管廊燃气泄漏浓度场预测矫正及泄漏率估计方法
本专利技术属于城市公共安全领域，具体涉及一种综合管廊舱内燃气管道泄漏扩散浓度场预测矫正及泄漏率估计方法。
技术介绍
城市地下综合管廊将设置在地面、地下或架空的各类公用类管线集中容纳于一体，充分利用了城市地下空间，避免了路面的反复开挖、确保了道路交通功能的充分发挥，创造了良好的市民生活环境。综合管廊的建设对充分地利用地下空间，节省投资，对拉动经济发展、改变城市面貌、保障城市安全都具有重要作用。但同时，在综合管廊内一种灾害的发生可能会招致其他类型灾害的发生，造成更大的危害，尤其是管廊内燃气管道的泄漏可能导致人员窒息伤亡和爆炸等事故，并可能引发一系列次生灾害，对城市的安全运行、公民的生命安全造成巨大威胁。目前，综合管廊内燃气管道发生泄漏后的主要应急措施为启动事故段及相邻段的通风设备，或紧急切断燃气管道分段阀门。因此，实现对燃气管道泄漏扩散浓度场的精准预测及泄漏率的估计对燃气泄漏后的应急处理十分关键。
技术实现思路
本专利技术的目的是提供一种综合管廊燃气泄漏浓度场预测矫正及泄漏率估计的技术方法，为城市地下综合管廊内燃气泄漏的应急处理与决策提供依据。为了快速、精确实现浓度场的预测矫正及泄漏率的估计，本专利技术制定了以下技术方案：本专利技术通过数值模拟方法结合现场实际监测数据实现了对燃气在管廊内部的浓度分布及泄漏率的计算预测。使用的计算方法为燃气泄漏扩散模型和改进的集合卡尔曼滤波算法。现场实际监测的数据包括：1.燃气探测器所监测的在探测器有效探测范围内的燃气浓度实时数据；2.风速风向传感器所监测的管廊内部风速、风向信息；3.温湿度传感器所监测的管廊内部空气的温度和湿度。由于在本专利技术的数值模拟部分中使用改进的集合卡尔曼滤波算法将燃气探测器的实时监测数据作为观测值对燃气扩散浓度场的预测值进行更新，并对燃气管道泄漏率进行反演，因此数值模拟的预测效果将受到传感器精度的影响，并且会受到燃气探测器的布置位置和间距的影响。依据对燃气扩散浓度场和泄漏率的计算效果，合理的设计燃气探测器的布置位置与间隔，可达到本专利技术的最佳预测效果。燃气在管廊内的浓度分布经过泄漏扩散模型初步计算后，借助改进的集合卡尔曼滤波算法利用实时燃气探测器监测数据进行数据同化更新后将得到燃气在管廊内部浓度场分布的准确预测值。由于各类传感器均存在一定误差，并且管廊内部空气流动为湍流，在风流速度和扩散系数上均存在很大的不确定性，本专利技术考虑了管廊内部风流速度和扩散系数的不确定性对燃气浓度场分布及泄漏率估计结果的影响，结合风速风向传感器和温湿度传感器的监测数据，在使用数值模拟部分时对风流速度与扩散系数加入符合实际情况的噪声，有效抑制了风流速度和扩散系数的不确定性对预测估计结果的影响。本专利技术可快速，精确的预测出燃气在管廊内的浓度分布并且可以实现对泄漏率的估计，这将大大提高管廊内燃气泄漏应急处理的效率和精确度。附图说明图1为本专利技术提出的实现管廊内部燃气管道泄漏扩散浓度场预测矫正及泄漏率估计的技术方案框架图。图2为本专利技术实施例提供的综合管廊燃气管道泄漏扩散浓度场预测矫正及泄漏率估计方法传感器布置示意图。图3为本专利技术提出的燃气泄露扩散浓度场预测矫正及泄漏率估计数值计算模型的程序流程图。图4为本专利技术实施例通过本专利技术方法预测矫正得到的450s时的燃气浓度场与真实燃气浓度场的对比情况。具体实施方式如图1本专利技术提出的管廊内部燃气管道泄漏扩散浓度场预测矫正及泄漏率估计的结构框图所示，本专利技术结合了数值模型计算与实时观测数据，将管廊内部燃气传感器实时监测数据(包括燃气探测器有效探测范围内的燃气浓度、管廊内部的风速、风向信息、管廊内部空气温度和湿度信息)输入数值计算模型(包括燃气泄漏扩散模型和改进的集合卡尔曼滤波算法)，进行数据同化计算后可以产生管廊内的燃气浓度分布情况以及泄漏率的预测估计结果。并且每进行一次数据同化将输出一次结果，该模型可连续进行模拟计算，直至使用者结束模拟计算。图2具体展示了本专利技术的一个实施例，列出了综合管廊燃气管道泄漏扩散浓度场预测矫正及泄漏率估计方法所使用的传感器布置。燃气探测器(2～10)从第20m处布置，间隔为20m。13和14分别为燃气管道和管座(墩)，15为综合管廊燃气仓的通风口，温湿度传感器11和风速风向传感器12均布置在管廊燃气仓内部。假定燃气管道45m处有泄漏源1，泄漏率为0.147m3/s，下面详细阐述本专利技术的实施方式：一维的对流扩散方程(ADE)作为该实施例燃气泄漏扩散模拟模型的控制方程；式中，c表示管廊燃气仓断面的燃气平均浓度，u表示风速，Dx表示扩散系数，q表示源项。利用MATLAB编程，使用有限差分方法(FDM)对上述控制方程离散后显式求解，为了保证计算的稳定，时间步长将会自动调整，使柯郎数始终小于1。由该一维对流扩散模型得到的计算结果将作为改进的集合卡尔曼滤波算法的预测值。改进的集合卡尔曼滤波算法的预测过程如下：式中，为k时刻状态集合中第i个状态预测值，M为模型算子，即k-1时刻到k时刻状态的变化关系，为k-1时刻状态集合中第i个状态分析值，ω为模型误差。当k时刻有观测值时，将利用观测值对集合中所有的状态进行更新：Pe＝X′X′T式中，是k时刻状态集合中第i个状态分析值；Ke是增益矩阵；yi(tk)是k时刻的观测数据；H表示观测算子；ν是观测误差，Re是其协方差矩阵；Pe为状态集合的误差方差矩阵；x表示状态向量，下标i(i＝1,2,…,N)表示每个状态向量在集合中的编号，N表示集合的大小，在本实施例中N取60。本专利技术改进的集合卡尔曼滤波算法中状态向量x由两部分组成：式中，在该模型中为管廊各处的燃气浓度平均值，y在该模型中为泄漏率，是不可直接观测的变量。泄漏率y通过状态向量集合的误差方差矩阵Pe更新，实现对泄漏率的反演。如图3本专利技术所构建的数值模拟模型的程序流程图所示，在本专利技术的方法中，第一次数据同化时的泄漏率由使用者完成初始化，在本实施例中泄漏率初始化为1m3/s，其标准差为0.1m3/s，而后每个数据同化过程的泄漏率将由改进的集合卡尔曼滤波算法通过误差方差矩阵完成更新后自动产生。本专利技术的数值模拟模型的计算流程为：1.使用者完成燃气泄漏扩散模型参数初始化和状态向量集合初始化，管廊断面的燃气平均浓度c初始化为0，风速u初始化为0.33m/s，扩散系数Dx初始化为0.0047m2/s，各状态向量初始化为：管廊各处的燃气浓度平均值均设为0，第一步数据同化的泄漏率y设为平均值为1m3/s，标准差为0.1m3/s的随机数。2.燃气泄漏扩散模型向前计算，将依据温湿度传感器和风速风向传感器的监测数据所得出的添加有噪声的风流速度和扩散系数带入模型计算，在本实施例中风流速度为0.33m/s，添加符合N(1,0.1)分布的白噪声，扩散系数为0.0047m2/s，同样添加符合N(1,0.1)分布的白噪声，得出的计算结果将作为改进的集合卡尔曼滤波算法的预测值保存。此时若有燃气浓度的观测数据，便使用改进的集合卡尔曼滤波算法进行数据同化，经过算法更新状态向量后，便得到这一阶段数据同化的结果。若没有燃气浓度观测数据，则跳过改进的集合卡尔曼滤波算法部分，直接判断模拟是否结束。在本实施例中每30s取一组燃气探测器监测的燃气浓度值作为本文档来自技高网...

【技术保护点】
1.一种综合管廊燃气泄漏浓度场预测矫正及泄漏率估计方法，其特征在于，结合综合管廊舱内实时监测数据，通过构建综合管廊舱内燃气泄漏扩散的数值计算模型，运用数据同化算法进行燃气泄漏扩散过程的计算预测，该方法能有效抑制管廊舱内风流速度和扩散系数不确定性对计算预测结果的影响，在获得管廊内部燃气浓度分布情况的同时，能准确对燃气管道泄漏率进行估计。

【技术特征摘要】
1.一种综合管廊燃气泄漏浓度场预测矫正及泄漏率估计方法，其特征在于，结合综合管廊舱内实时监测数据，通过构建综合管廊舱内燃气泄漏扩散的数值计算模型，运用数据同化算法进行燃气泄漏扩散过程的计算预测，该方法能有效抑制管廊舱内风流速度和扩散系数不确定性对计算预测结果的影响，在获得管廊内部燃气浓度分布情况的同时，能准确对燃气管道泄漏率进行估计。2.根据权利要求1所述的综合管廊燃气泄漏浓度场预测矫正及泄漏率估计方法，其特征在于，所述综合管廊舱内监测数据包括燃气传感器实时采集的燃气浓度信息和风速风向传感器所监测到的风向、风速信息以及管廊内的温度、湿度信息，燃气传感器实时采集的燃气浓度作为数据同化算法的观测值矫正燃气扩散浓度场并反演泄漏率，风速风向传感器所监测到的风向、风速信息以及管廊内的温度、湿度信息，用于准确地确定管廊舱内风速和管廊舱内燃气扩散系数。3.根据权利要求1所述的综合管廊燃气泄漏浓度场预测矫正及泄漏率估计方法，其特征在于，所述综合管廊内燃气泄漏扩散数值计算模型采用有限差分方法对燃气对流扩散控制方程进行离散求解，得到管廊内各处的燃气浓度预测值，而后运用数据...

【专利技术属性】
技术研发人员：吴建松，原帅琪，张灿，
申请(专利权)人：中国矿业大学北京，
类型：发明
国别省市：北京,11