[发明专利]基于逆瞬态模型的城市非金属管道泄漏定位方法

权利要求书

1.一种基于逆瞬态模型的城市非金属管道泄漏定位方法，其特征在于：包括如下步骤：

S1：用AnsysIcem对管道进行建模，使用Fluent软件模拟非金属管道泄漏和无泄漏状态下的压力、流量和流速参数分布，分别找出非金属管道泄漏和无泄漏状态下的压力、流量和流速参数随时间和空间变化的规律；

S2：在实验室进行非金属埋地管道泄漏模拟试验，分别在泄漏和未泄漏两种状态下测量管道流量、温度及多处的压力参数数据，与步骤S1中获得的模拟结果对比和验证，突出泄漏状态下的规律，获得实验管道测量到的测压头；

S3：根据步骤S1得到的压力、流量和流速的变化规律，建立非金属管道泄漏状态的逆瞬态模型；

为实现计算出的测压头和实验测量到的测压头之间的加权均方误差最小化，逆瞬态模型具体如下：

其中，OF为目标函数；DT为观测到的水力瞬变事件的时间步数；H为通过瞬态求解器计算出的测压头；H^*为实验管道测量到的测压头；n^H为系统中压力观测点的数量；H_t,j为不同测量节点和泄漏时间下的通过瞬态求解器计算出的测压头，H^*_t,j不同测量节点和泄漏时间下的实验设施中测量到的测压头；j为测量节点；t为泄漏后的时间；

S4：应用全局和局部算法，得出Kelvin-Voigt元素的数目以及收敛性最好的目标函数值J_k和τ_k，确定逆瞬态模型中通过瞬态求解器计算出的测压头H，具体步骤包括：

S4.1：在matlab环境下，使用有效的全局搜索方法PSO粒子群算法，适应度线性定标，随机余数采样SRS作为选择方案，每个PSO评估使用九十代，找到Kelvin-Voigt元素的数目；

S4.2：分别投放1，10，20，50，100，200，500个随机种子，观察目标函数OF的函数值收敛情况，得出全局搜索下收敛性最佳时的J_k和τ_k值，其中Z＝{J_k，τ_k}，k＝1到N_KV，其中，N_KV为Kelvin-Voigt元素的数目；

S4.3：使用局部搜索LM算法来估计和优化系数J_k，Z＝{J_k，τ_k}；k＝1到N_KV，并考虑由PSO确定的区域中τ_k的值，计算出三种常用的K值和使用LM方法的目标函数值，并绘折线图，挑选收敛性最好的目标函数值J_k和τ_k；

S4.4：将Kelvin-Voigt元素的数目和收敛性最好的目标函数值J_k和τ_k值代入蠕变函数方程及MOC特征方程，从而得出H值；

S5：将实验测量到的测压头与计算出的测压头共同代入逆瞬态模型，得到目标函数值OF；

S6：在实验所得多处的OF值中，挑选最小值的位置节点，并在最小值位置节点前后布置等距的且距离小于前一次，个数与其相同的多个测量节点，重复步骤计算各节点H值，并得到OF，当各节点距离小于管道检测长度的4％时终止重复步骤。

2.如权利要求1所述的基于逆瞬态模型的城市非金属管道泄漏定位方法，其特征在于：步骤4.4所述的蠕变函数方程具体为：

J(t)＝e_r(t)/σ₀

所述蠕变函数方程用以描述材料的延迟应变e_r(t)与载荷σ₀之间的比率，为解决特征方程中无法直接得到的延迟应变时间导数而引入。