一种基于数学模型的气体管道泄漏检测定位的方法与流程

本发明涉及管道检漏技术领域，尤其涉及一种基于数学模型的气体管道泄漏检测定位的方法。

背景技术：

在类似于托卡马克这种大型装置中，常常配有各种复杂形状的输运管道，管道的良好密封性才能保证整个装置的正常运行，现实中常用真空检漏定位的方法来检查管道的密封性，现有的真空检漏定位手段主要是正压冒泡法、正压吸枪法及负压喷氦气方法，但对于长距离且置于复杂结构中的大型管道的检测效果不佳且操作困难，并且这些检漏定位方法必须要破除托卡马克真空室的真空环境才能够进行，这就导致了在管道检漏定位完成后，需要重新清洁真空室，并还原真空状态，整个检漏定位过程流程繁琐且耗时过长。

技术实现要素：

本发明目的就是为了弥补已有技术的缺陷，提供一种基于数学模型的气体管道泄漏检测定位的方法。

本发明是通过以下技术方案实现的：

一种基于数学模型的气体管道泄漏检测定位的方法，根据不同时刻管道压强的值并结合初始的物理参数，通过数学模型得出管道泄漏点的位置；所述的测试管道初始内部为正压状态，测试过程中管道两端封闭的；具体步骤如下：

步骤s01：记录下管道的原始参数；

步骤s02：进行测试并收集管道两端和中部的压力随时间变化的三组数据；

步骤s03：将记录下的管道原始参数作为数学模型的输入量，测试得到的任一组数据作为算法中数学模型的输入变量，管道泄漏位置作为数学模型的输出量；

步骤s04：选择另一组测试得到的数据作为算法中数学模型的输入变量，重复步骤s03，得出新的输出量；当将步骤s02测得的数据都代入数学模型中后，通过智能算法软件matlab处理数据，得出管道泄漏位置。

步骤s01中，所述的管道的原始参数包括有：测试管道的直径、长度、待测点位置和管道内外的初始压力。

步骤s02中，收集的三组数据是在同一次实验中的同一时刻收集的；压力的测量采用压力传感器测量。

步骤s03中采用的数学模型的推导是依据流体力学运动方程连续性方程状态方程所述数学模型的推导用到的数学公式包括偏微分方程求解和拉普拉斯变换，最终得出的数学模型为

当0≤x≤l时

当l≤x≤l时

上述式中p为压力，pl为漏孔左侧的压力，pr为漏孔右侧的压力，pa；ρ为气体密度，kg/m³；v为气体流速，m/s；λ为水力摩阻系数；d为管道直径，m；g为重力加速度，m/s²；α为管道水平倾角；c为等温条件下声音在气体中传播速度，m/s；z为气体的压缩系数；r为气体常数；t为绝对温度，k；a²为线性化系数，m²/s，ph为管道内气体的初始压力，pa；l为假设的泄露点位置，km；为泄露函数；l为管道的长度，km；t为时间变量，s；x为气体在管道中流动的位置坐标，km。

本发明的优点是：本发明通过将测量的管道内压力作为数学模型的输入变量，模型的输出量为泄漏点的位置，操作过程简单，定位精度高，特别的是该方法是在管道内部进行的，不会破坏管道外部的真空腔的真空度，特别是该方法可以用于管径很小的管道，耗时少，并能达到高效的检漏定位水平。

附图说明

图1为本发明的流程图。

图2是本发明的结构示意图。

图3是本发明实例测试管道的重要参数示意图。

具体实施方式

一种基于数学模型的气体管道泄漏检测定位的方法，根据不同时刻管道压强的值并结合初始的物理参数，通过数学模型得出管道泄漏点的位置；所述的测试管道初始内部为正压状态，测试过程中管道两端封闭的；具体步骤如下：

步骤s01：记录下管道的原始参数；

步骤s02：进行测试并收集管道两端和中部的压力随时间变化的三组数据；

步骤s03：将记录下的管道原始参数作为数学模型的输入量，测试得到的任一组数据作为算法中数学模型的输入变量，管道泄漏位置作为数学模型的输出量；

步骤s04：选择另一组测试得到的数据作为算法中数学模型的输入变量，重复步骤s03，得出新的输出量；当将步骤s02测得的数据都代入数学模型中后，通过智能算法软件matlab处理数据，得出管道泄漏位置。

如图1所示步骤s01，记录下管道的原始参数；管道的原始数据包括:如图2所示，测试管道2的直径d，长度l，待测点位置l1、l2、l3，管道内外的初始压力pi、po；如图2将测试管道放在一段大管道1中，开始会对大管道1抽气至po，从而模拟真空腔室内的冷却水管的正压状态。

步骤s02，进行测试并收集管道两端和中部的压力随时间变化的三组数据；测试中会将压力传感器3放置在如图2、3所示的待测点位置l1、l2、l3处，测得的压力为p1、p2、p3；压力传感器3会记录不同时刻下的待测点压力值，从而得到很多组数据例如(t1、p1t1、p2t1、p3t1)；(t2、p1t2、p2t2、p3t2)；(t3、p1t3、p2t3、p3t3)等等。

步骤s03,将步骤s01记录下的管道原始参数(d、l、l1、l2、l3、pi、po)作为数学模型的输入量，步骤s02测试得到的t1时刻数据(t1、p1t1、p2t1、p3t1)作为数学模型的输入变量，数学模型的输出量为管道泄漏位置lt1。

步骤s03中采用的数学模型的推导是依据流体力学运动方程连续性方程状态方程所述数学模型的推导用到的数学公式包括偏微分方程求解和拉普拉斯变换。最终得出的数学模型为

当0≤x≤l时

当l≤x≤l时

上述式中p为压力，pl为漏孔左侧的压力，pr为漏孔右侧的压力，pa；ρ为气体密度，kg/m³；v为气体流速，m/s；λ为水力摩阻系数；d为管道直径，m；g为重力加速度，m/s²；α为管道水平倾角；c为等温条件下声音在气体中传播速度，m/s；z为气体的压缩系数；r为气体常数；t为绝对温度，k；α²为线性化系数，m²/s，ph为管道内气体的初始压力，pa；l为假设的泄露点位置，km；为泄露函数；l为管道的长度，km；t为时间变量，s；x为气体在管道中流动的位置坐标，km。

步骤s04,选择另一时刻下数据作为数学模型的输入变量，重复步骤s03得出新的输出量；当将步骤s02测得的所有数据都代入数学模型中后，通过智能算法软件matlab处理所有输出量(lt1，lt2，lt3，lt4...)并得出其分布情况，最终得出出现的频率最大的l值或者范围，这即为所需要的精确泄漏位置。