一种基于遗传算法的贝叶斯网络天然气管道泄漏概率计算方法与流程

技术特征：
1.一种基于遗传算法的贝叶斯网络天然气管道泄漏概率计算方法，其特征在于，包括以下步骤：s1、基于统计数据库分析天然气管道失效原因因素,并分析影响天然气管道失效的实际状态因素，确定管道失效的实际状态因素与历史失效原因因素间的因果关联关系；s2、根据因果关联关系，建立天然气管道失效故障树模型，利用故障树与贝叶斯网络的转换关系确定贝叶斯网络结构，完成天然气管道失效贝叶斯网络模型结构的构建；s3、引入遗传算法完成贝叶斯网络参数学习以此获取最优网络条件概率；s4、采用genie贝叶斯网络仿真软件建立模型，输入贝叶斯网络根节点概率与条件概率，计算天然气管道在实际状态下的失效概率。2.根据权利要求1所述的一种基于遗传算法的贝叶斯网络天然气管道泄漏概率计算方法，其特征在于，还包括s5、利用贝叶斯网络的反向推理功能，对每个风险因素进行重要度分析并得到影响管道发生失效的关键致因链。3.根据权利要求1所述的一种基于遗传算法的贝叶斯网络天然气管道泄漏概率计算方法，其特征在于，所述s1中，采用思维导图分别分析基于统计数据库的管道失效原因因素与管道失效的实际状态因素，并建立关联关系图。4.根据权利要求1所述的一种基于遗传算法的贝叶斯网络天然气管道泄漏概率计算方法，其特征在于，所述s2中，故障树与贝叶斯网络的转变步骤如下：s2.1、故障树中的事故与贝叶斯网络的节点进行一一对应；s2.2、按照故障树中事故的逻辑因果关系用有向边进行贝叶斯网络节点间的连接；s2.3、故障树中的逻辑门与节点条件概率分布表进行对应。5.根据权利要求1所述的一种基于遗传算法的贝叶斯网络天然气管道泄漏概率计算方法，其特征在于，所述s3中，基于遗传算法的条件概率θ
i
的计算步骤如下：s3.1、染色体编码，集合θ中每个参数均在[0,1]区间内，精度为二进制表示小数点后8位；s3.2、种群初始化，在计算前以随机方式生成初始种群p(t)，个体为a,设置染色体长度为k,其值取决于条件概率表中参数的个数n,k＝8n；设置种群规模m,即随机产生m个k位随机数形成初始种群,种群中每个个体对应一组条件概率分布；s3.3、适应度函数，选取最大似然函数l(θ)作为适应度函数f(x)，公式如下：式中：i为网络拓扑结构g的节点个数；j为拓扑结构下节点δ(x
i
)的组合数目；k为当前节点的取值、为一个常数；n
ijk
为管道风险样本的特征函数；θ
ijk
为非根节点对应的一组条件概率，对于初始种群p(t)中的每一个个体p
i
(t)，构建节点x
i
与其父节点p
g
(x
i
)的适应度函数f(x)，计算每个个体的适应度值；s3.4、选择操作，利用轮盘赌方式从p(t)中选择最优个体，得到选择种群p1(t 1)；s3.5、交叉操作，确定交叉概率p
c
,一般为0.65～0.90，采用离散重组式多切点交叉法进行交叉运算，得到交叉种群p2(t 1)，s3.6、变异操作，在交叉种群p2(t 1)中以变异概率p
m
选择个体进行变异操作,变异过程采用反转位值法，即对二进制位进行反转操作，得到变异种群p3(t 1)、变异种群p3(t 1)作
为新种群p(t 1)代；s3.7、终止条件，对每一代群体进行选择、交叉和变异运算，设最大遗传代数为t,当进行t轮选择、交叉和变异运算后，得到新种群p(t t)代、则满足终止条件，输出最优网络参数。6.根据权利要求5所述的一种基于遗传算法的贝叶斯网络天然气管道泄漏概率计算方法，其特征在于，所述s3.4中的轮盘赌方式，具体操作如下：s3.4.1、根据适应度函数计算初始种群p(t)中每个个体的适应度值f(a
i
)，i＝1,2,...,m,m为种群大小；s3.4.2、计算每个个体被遗传到下一代的入选概率f(a
i
)：s3.4.3、计算每个个体的累积概率f
i
：s3.4.4、在[0，1]区间内产生一个均匀分布的随机数γ；s3.4.5、若γ<f1，则选择个体a1，否则，选择个体q，使得：f
q-1
<γ<f
q
成立；s3.4.6、重复步骤s3.4.4和s3.4.5多次，被选中次数最多的个体为初始种群p(t)代的最优个体，被选中次数最少且入选概率最小的个体为初始种群p(t)代最劣个体；淘汰最劣个体，让最优个体替代最劣个体得到选择种群p1(t 1)。7.根据权利要求5所述的一种基于遗传算法的贝叶斯网络天然气管道泄漏概率计算方法，其特征在于，所述s3.5中的交叉概率p
c
，具体计算公式如下：式中：p
c
为交叉概率；f0为当前个体的适应度值；δ为设定的阈值；f
max
为当前种群中最高适应度值；f
av
为当前种群的平均适应度值；f
std
为当前种群的适应度值标准差。8.根据权利要求5或7所述的一种基于遗传算法的贝叶斯网络天然气管道泄漏概率计算方法，其特征在于，所述的离散重组式多切点交叉法进行交叉运算如下：s3.5.1、根据标准差f
std
自动选择交叉方式，s3.5.2、选择种群p1(t 1)中的每两个个体p1和p2随机配对,以两个体中较高概率p
c
判断是否进行交叉操作；在[0,1]区间产生一个随机数e，若e<p
c
，则进行交叉操作，否则继续产生随机数判断之后两个个体；s3.5.3、对进行交叉操作的两个个体随机选取k/8个切点,切点位置相同，等位交换切点之间的子串，若f
std
≥1.2δ，交换子串数为双子串,若f
std
≤0.8δ，交换子串数为单子串,交叉后得到交叉种群p2(t 1)。9.根据权利要求5所述的一种基于遗传算法的贝叶斯网络天然气管道泄漏概率计算方
法，其特征在于，所述s3.6变异操作，具体步骤如下所示：s3.6.1、假设变异概率p
m
为0.01～0.3，式中：p
m
为变异概率；s3.6.2、变异方式的设定：若交叉种群的适应度值标准差f
std
≤0.8δ时，说明交叉种群适应度集中，应增大变异概率，对交叉种群p2(t 1)中被选中个体进行两点变异；当f
std
≥1.2δ时，对被选中个体进行单点变异；s3.6.3、在[0,1]区间产生一个随机数e，若e<p
m
，则进行变异操作，否则继续产生随机数判断之后的个体；通过变异操作得到变异种群p3(t 1)，变异种群p3(t 1)作为新一代的新种群p(t 1)代。

技术总结
本发明公开了一种基于遗传算法的贝叶斯网络天然气管道泄漏概率计算方法，涉及天然气管道安全技术领域，为了解决天然气管道泄漏概率计算脱离管道实际状态的问题，包括以下步骤：基于统计数据库分析天然气管道失效原因因素,并分析影响天然气管道失效的实际状态因素，根据因果关联关系，建立天然气管道失效故障树模型，利用故障树与贝叶斯网络的转换关系确定贝叶斯网络结构，完成天然气管道失效贝叶斯网络模型结构的构建；引入遗传算法完成贝叶斯网络参数学习以此获取最优网络条件概率；采用GeNIe贝叶斯网络仿真软件建立模型，输入贝叶斯网络根节点概率与条件概率，计算天然气管道在实际状态下的失效概率。道在实际状态下的失效概率。道在实际状态下的失效概率。


技术研发人员：韩林序 唐春凌 张林霞 古冉 赵松 孙明炬
受保护的技术使用者：中国石油天然气股份有限公司
技术研发日：2021.12.20
技术公布日：2022/3/22