一种天然气管网供气可靠性的评价方法及装置的制造方法

一种天然气管网供气可靠性的评价方法及装置的制造方法
【专利摘要】本发明提供了一种天然气管网供气可靠性的评价方法及装置，涉及天然气传输技术领域。方法包括：建立天然气管网的容量网络；确定各管道单元在各预设失效模式下的第一失效概率，并确定各管道单元在一预设时间内的第一状态分布；确定各天然气站场单元的第二失效概率，并确定各天然气站场单元在一预设时间内的第二状态分布；实时更新所述容量网络；向天然气管网中的各用户供气，确定各用户的供气情况信息；根据所述各用户的供气情况信息确定向各用户供气的可靠度以及各用户的平均缺气量。通过本发明可以准确评价在不同干扰因素下天然气管网的供气可靠性，解决了当前天然气管网的供气优劣难以确认的问题。
【专利说明】
一种天然气管网供气可靠性的评价方法及装置
技术领域
[0001] 本发明涉及天然气传输技术领域，尤其涉及一种天然气管网供气可靠性的评价方 法及装置。
【背景技术】
[0002] 目前，天然气管网是天然气输送的主要手段，也是连接资源和市场的纽带。近年 来，中国天然气管网规模不断扩大，对天然气这种清洁能源的需求不断提高，确保天然气管 网可靠供气的重要性日益凸显。天然气管网的供气可靠性是指在一段时间内，在各种不确 定性因素的干扰下，天然气管网能够满足用户用气需求的能力。目前天然气管网的可靠性 评价处于起步阶段。其中，对天然气管网的运行可靠性关注较多，然而针对供气可靠性评价 方法的研究还不够成熟。此处的运行可靠性和供气可靠性是两个不同的概念。运行可靠性 的评价对象是天然气管网的运行状态是否达到期望或要求;供气可靠性更多关注管网在天 然气供、需之间发挥的作用是否达标。运行可靠时，供气不一定可靠，这是由于管网的拓扑 结构复杂性和物理输送能力的限制作用；而供气可靠时运行也并不一定可靠，这是由于管 网的物理传输能力可能存在冗余。
[0003] 当前，有很多研究者试图采用运行可靠性领域的评价方法来分析供气可靠性。考 虑到运行环境与评价对象存在巨大差异，研究者一般都会认识到经典可靠性理论不足以解 决天然气管网可靠性评价过程中面临的诸多问题，如环境因素、人为活动引起的不确定性 等。当前，还没有准确的天然气管网供气可靠性评价的方法，使得天然气管网的供气优劣难 以确认。

【发明内容】

[0004] 本发明的实施例提供一种天然气管网供气可靠性的评价方法及装置，以解决当前 还没有准确的天然气管网供气可靠性评价的方法，天然气管网的供气优劣难以确认的问 题。
[0005] 为达到上述目的，本发明采用如下技术方案：
[0006] -种天然气管网供气可靠性的评价方法，包括：
[0007] 获取天然气管网的拓扑信息以及天然气管网拓扑中各管道单元信息和各天然气 站场单元信息；
[0008] 根据天然气管网的拓扑信息建立天然气管网的容量网络;所述容量网络包括多个 节点和各节点之间的弧;所述容量网络中的节点为各天然气站场单元;各节点之间的弧为 各管道单元；
[0009] 根据所述各管道单元信息确定各管道单元在各预设失效模式下的第一失效概率；
[0010] 根据所述第一失效概率确定各管道单元在一预设时间内的第一状态分布;所述第 一状态分布包括各管道单元处于正常状态或各预设失效模式状态；
[0011] 根据所述各天然气站场单元信息确定各天然气站场单元的第二失效概率；
[0012] 根据所述第二失效概率确定各天然气站场单元在一预设时间内的第二状态分布； 所述第二状态分布包括各天然气站场单元处于正常状态或失效状态；
[0013] 根据各管道单元信息以及所述第一状态分布确定各管道单元在正常状态和各预 设失效模式状态下的输气能力，并将所述输气能力确定为所述容量网络的各弧容量；
[0014] 根据所述各弧容量和所述第二状态分布，实时更新所述容量网络；
[0015] 根据实时更新的容量网络向天然气管网中的各用户供气，确定各用户的供气情况 信息；
[0016] 根据所述各用户的供气情况信息确定向各用户供气的可靠度以及各用户的平均 缺气量。
[0017] 具体的，所述预设失效模式包括管道腐蚀；
[0018] 根据所述各管道单元信息确定各管道单元在各预设失效模式下的第一失效概率， 包括：
[0019] 根据极限状态方程建立管道失效概率计算模型，并根据蒙特卡洛模拟方法计算所 述管道失效概率计算模型，获取缺陷运行年限为t时的管道上每个缺陷的失效概率Ff(t);
[0020] 其中，？￡(〇=?&￡(〇-?0<0]^￡(〇为缺陷运行年限为七时的管道缺陷的失 效压力；Po为管道的实际运行压力；
-;其中，M为管道膨胀系数，
;〇f为管道的流变应力；
[0021] D为管道公称直径;Wt为管道壁厚;L(t)为t时刻管道上缺陷的长度;L(t)=L(0) + rut; L (0)为0时刻管道上缺陷的长度;n为管道轴向腐蚀速率；
[0022] 根据齐次泊松模型，确定每千米管道上的缺陷数目k;
[0023] 其中，该齐次泊松模型为
1，2~k;其中，Pr(n)为每千米管道 上出现n个缺陷的概率;k为每千米管道上的缺陷数目；x为根据内检测结果统计得到的每千 米管道上的缺陷数目的平均值；
[0024] 根据公式:A(t)=k ? Ff(t)确定一管道上在缺陷运行年限为t时的失效概率A(t)。
[0025] 具体的，根据所述第一失效概率确定各管道单元在一预设时间内的第一状态分 布，包括：
[0026] 根据公式：
Y =0，1，2…，确定各管道单元在一预设时 间〇至t内发生故障次数为Y的概率;其中，PKWOir/]为各管道单元在一预设时间0至t 内发生故障次数为Y的概率；为管道单元长度;A(t)为管道上在缺陷运 行年限为t时的失效概率；
[0027]根据公式：>确定各管道单元从一to时刻开始下一个 故障出现的时间大于t时间的概率;Pr[Y(to)>t]为各管道单元从一 to时刻开始下一个故障 出现的时间大于t时间的概率;Y(to)为从一 to时刻开始下一个故障出现的时间；
[0028]根据PrtWOir/]和Pr[Y(t0)>t]确定各管道单元在所述预设时间内的第一状 态分布。
[0029] 此外，所述各天然气站场单元信息包括各天然气站场单元的历史失效情况数据；
[0030] 根据所述各天然气站场单元信息确定各天然气站场单元的第二失效概率，包括： [0031 ]根据各天然气站场单元的历史失效情况数据，确定各天然气站场单元的第二失效 概率。
[0032]具体的，根据所述第二失效概率确定各天然气站场单元在一预设时间内的第二状 态分布，包括：
[0033] 根据公式
Y =0，1，2…，确定各天然气站场单元在一 预设时间0至t内发生故障次数为Y的概率;其中，PV [Y (0=!/ ]为各天然气站场单元在 一预设时间〇至t内发生故障次数为Y的概率；#(/)=￡./+(/)dt ;f (t)为各天然气站场单元的 第二失效概率；
[0034] 根据公式：丨丨> ^ = 确定各天然气站场单元从一to时刻开始下 一个故障出现的时间大于t时间的概率;PV [V (to)>t]为各天然气站场单元从一to时刻开 始下一个故障出现的时间大于t时间的概率;V (to)为从一to时刻开始下一个故障出现的时 间；
[0035]根据Pr' [Y (0=!/ ]和PV [Y' (to)>t]确定各天然气站场单元在所述预设时间 内的第二状态分布。
[0036]具体的，所述各管道单元信息包括各管道的状态与输气能力的对应关系；
[0037]根据各管道单元信息以及所述第一状态分布确定各管道单元在正常状态和各预 设失效模式状态下的输气能力，包括：
[0038]根据第一状态分布确定一预设时间点各管道单元的状态；
[0039]根据预设时间点各管道单元的状态，在所述对应关系中确定各管道单元在正常状 态和各预设失效模式状态下的输气能力。
[0040] 具体的，所述各用户的供气情况信息包括：向用户模拟的供气次数N;向用户模拟 的供气次数中供气不足的次数K;向用户模拟的供气年份Y;每次向用户模拟的供气不足造 成的缺气量Z;
[0041] 根据所述各用户的供气情况信息确定向各用户供气的可靠度以及各用户的平均 缺气量，包括：
[0042] 根据公式： 确定在年份Y向各用户供气的可靠度R;
[0043] 根据公式：
，确定向一用户模拟的供气平均缺气量Q。
[0044] -种天然气管网供气可靠性的评价装置，包括：
[0045] 信息获取单元，用于获取天然气管网的拓扑信息以及天然气管网拓扑中各管道单 元信息和各天然气站场单元信息；
[0046] 容量网络建立单元，用于根据天然气管网的拓扑信息建立天然气管网的容量网 络;所述容量网络包括多个节点和各节点之间的弧;所述容量网络中的节点为各天然气站 场单元;各节点之间的弧为各管道单元；
[0047] 第一失效概率确定单元，用于根据所述各管道单元信息确定各管道单元在各预设 失效模式下的第一失效概率；
[0048] 第一状态分布确定单元，用于根据所述第一失效概率确定各管道单元在一预设时 间内的第一状态分布;所述第一状态分布包括各管道单元处于正常状态或各预设失效模式 状态；
[0049] 第二失效概率确定单元，用于根据所述各天然气站场单元信息确定各天然气站场 单元的第二失效概率；
[0050] 第二状态分布确定单元，用于根据所述第二失效概率确定各天然气站场单元在一 预设时间内的第二状态分布;所述第二状态分布包括各天然气站场单元处于正常状态或失 效状态；
[0051] 管道单元输气能力确定单元，用于根据各管道单元信息以及所述第一状态分布确 定各管道单元在正常状态和各预设失效模式状态下的输气能力，并将所述输气能力确定为 所述容量网络的各弧容量；
[0052] 容量网络更新单元，用于根据所述各弧容量和所述第二状态分布，实时更新所述 容量网络；
[0053]模拟供气单元，用于根据实时更新的容量网络向天然气管网中的各用户供气，确 定各用户的供气情况信息；
[0054]可靠性评价单元，用于根据所述各用户的供气情况信息确定向各用户供气的可靠 度以及各用户的平均缺气量。
[0055] 此外，所述预设失效模式包括管道腐蚀；
[0056] 该第一失效概率确定单元，具体用于：
[0057]根据极限状态方程建立管道失效概率计算模型，并根据蒙特卡洛模拟方法计算所 述管道失效概率计算模型，获取缺陷运行年限为t时的管道上每个缺陷的失效概率Ff(t); [0058] 其中，Ff(t)=P[pf(t)-Po彡0] ;pf(t)为缺陷运行年限为t时的管道缺陷 的失效压力；Po为管道的实际运行压力；
；其中，M为管道膨胀系数，
;〇f为管道的流变应力；
[0059] D为管道公称直径;Wt为管道壁厚;L(t)为t时刻管道上缺陷的长度;L(t)=L(0) + nt; L (0)为0时刻管道上缺陷的长度;n为管道轴向腐蚀速率；
[0060] 根据齐次泊松模型，确定每千米管道上的缺陷数目k;
[0061] 其中，该齐次泊松模型为：
，n = 0，l，2…k;其中，Pr(n)为每千米管道 上出现n个缺陷的概率;k为每千米管道上的缺陷数目；x为根据内检测结果统计得到的每千 米管道上的缺陷数目的平均值；
[0062] 根据公式:A(t)=k ? Ff(t)确定一管道上在缺陷运行年限为t时的失效概率A(t)。
[0063] 此外，所述第一状态分布确定单元，具体用于：
[0064] 根据公式：
r/=0,1，2…，确定各管道单元在一预设时 间〇至t内发生故障次数为Y的概率;其中，PKWOir/]为各管道单元在一预设时间0至t 内发生故障次数为Y的概率；#(〇=|^.7(/迚丄为管道单元长度4(〇为管道上在缺陷运 行年限为t时的失效概率；
[0065] 根据公式：^ ~)dt >确定各管道单元从一to时刻开始下一个 故障出现的时间大于t时间的概率;Pr[Y(to)>t]为各管道单元从一 to时刻开始下一个故障 出现的时间大于t时间的概率;Y(to)为从一 to时刻开始下一个故障出现的时间；
[0066]根据PrtWOir/]和Pr[Y(t0)>t]确定各管道单元在所述预设时间内的第一状 态分布。
[0067] 此外，所述各天然气站场单元信息包括各天然气站场单元的历史失效情况数据；
[0068] 所述第二失效概率确定单元，具体用于：
[0069] 根据各天然气站场单元的历史失效情况数据，确定各天然气站场单元的第二失效 概率。
[0070] 此外，所述第二状态分布确定单元，具体用于：
[0071] 根据公式：
=0,1，2…，确定各天然气站场单元在 一预设时间〇至t内发生故障次数为Y的概率;其中，PV [Y (0=!/ ]为各天然气站场单元 在一预设时间〇至t内发生故障次数为Y的概率，(/)〇(/)* ;f⑴为各天然气站场单元 的第二失效概率；
[0072] 根据公式：q = e_Ir/(i>kA确定各天然气站场单元从一to时刻开始下 一个故障出现的时间大于t时间的概率;PV [疒（to)>t]为各天然气站场单元从一to时刻开 始下一个故障出现的时间大于t时间的概率;V (to)为从一to时刻开始下一个故障出现的时 间；
[0073]根据Pr' [Y (0=!/ ]和PV [Y7 (to)>t]确定各天然气站场单元在所述预设时间 内的第二状态分布。
[0074]另外，所述各管道单元信息包括各管道的状态与输气能力的对应关系；
[0075] 该管道单元输气能力确定单元，具体用于：
[0076] 根据第一状态分布确定一预设时间点各管道单元的状态；
[0077]根据预设时间点各管道单元的状态，在所述对应关系中确定各管道单元在正常状 态和各预设失效模式状态下的输气能力。
[0078] 具体的，所述各用户的供气情况信息包括：向用户模拟的供气次数N;向用户模拟 的供气次数中供气不足的次数K;向用户模拟的供气年份Y;每次向用户模拟的供气不足造 成的缺气量Z;
[0079] 该可靠性评价单元，具体用于根据公式： 1确定在年份Y向各用户 ，
供气的可靠度R;
[0080] 根据公式 确定向一用户模拟的供气平均缺气量Q。
[0081] 本发明实施例提供的天然气管网供气可靠性的评价方法及装置，确定各管道单元 和天然气站场单元的失效概率，以及实时更新天然气管网的容量网络，并向天然气管网中 的各用户供气，确定各用户的供气情况信息；从而根据各用户的供气情况信息确定向各用 户供气的可靠度以及各用户的平均缺气量。通过本发明可以准确评价在不同干扰因素下天 然气管网的供气可靠性，解决了当前天然气管网的供气优劣难以确认的问题。
【附图说明】
[0082] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可 以根据这些附图获得其他的附图。
[0083] 图1为本发明实施例提供的一种天然气管网供气可靠性的评价方法的流程图一；
[0084] 图2为本发明实施例提供的一种天然气管网供气可靠性的评价方法的流程图二；
[0085] 图3为本发明实施例中的天然气管网的抽象图；
[0086] 图4为本发明实施例中的天然气管网中的一用户获取天然气量的累积概率密度函 数示意图；
[0087] 图5为本发明实施例提供的一种天然气管网供气可靠性的评价装置的结构示意 图。
【具体实施方式】
[0088] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0089] 如图1所示，本发明实施例提供一种天然气管网供气可靠性的评价方法，包括：
[0090] 步骤101、获取天然气管网的拓扑信息以及天然气管网拓扑中各管道单元信息和 各天然气站场单元信息。
[0091]步骤102、根据天然气管网的拓扑信息建立天然气管网的容量网络。
[0092 ]容量网络包括多个节点和各节点之间的弧;容量网络中的节点为各天然气站场单 元;各节点之间的弧为各管道单元。
[0093] 步骤103、根据各管道单元信息确定各管道单元在各预设失效模式下的第一失效 概率。
[0094] 步骤104、根据第一失效概率确定各管道单元在一预设时间内的第一状态分布。 [0095]第一状态分布包括各管道单元处于正常状态或各预设失效模式状态。
[0096] 步骤105、根据各天然气站场单元信息确定各天然气站场单元的第二失效概率。
[0097] 步骤106、根据第二失效概率确定各天然气站场单元在一预设时间内的第二状态 分布。
[0098]其中，第二状态分布包括各天然气站场单元处于正常状态或失效状态。
[0099]步骤107、根据各管道单元信息以及第一状态分布确定各管道单元在正常状态和 各预设失效模式状态下的输气能力，并将输气能力确定为容量网络的各弧容量。
[0100] 步骤108、根据各弧容量和第二状态分布，实时更新容量网络。
[0101] 步骤109、根据实时更新的容量网络向天然气管网中的各用户供气，确定各用户的 供气情况信息。
[0102] 步骤110、根据各用户的供气情况信息确定向各用户供气的可靠度以及各用户的 平均缺气量。
[0103] 本发明实施例提供的一种天然气管网供气可靠性的评价方法，确定各管道单元和 天然气站场单元的失效概率，以及实时更新天然气管网的容量网络，并向天然气管网中的 各用户供气，确定各用户的供气情况信息；从而根据各用户的供气情况信息确定向各用户 供气的可靠度以及各用户的平均缺气量。通过本发明实施例提供的天然气管网供气可靠性 的评价方法，可以准确评价在不同干扰因素下天然气管网的供气可靠性，解决了当前天然 气管网的供气优劣难以确认的问题。
[0104] 为了使本领域的技术人员更好的了解本发明，下面列举一个更为详细的实施例， 如图2所示，本发明实施例提供一种天然气管网供气可靠性的评价方法，包括：
[0105] 步骤201、获取天然气管网的拓扑信息以及天然气管网拓扑中各管道单元信息和 各天然气站场单元信息。
[0106] 该各天然气站场单元信息可以包括各天然气站场单元的历史失效情况数据、各天 然气站场单元的配置信等。该各管道单元信息包括各管道的状态与输气能力的对应关系， 还可以包括各管道的参数(长度，管壁厚度、管壁周长等）、管道内的检测结果（压力、温度 等)。
[0107] 步骤202、根据天然气管网的拓扑信息建立天然气管网的容量网络。
[0108] 此处，步骤202采用图论中的网络流表征方法建立天然气管网的容量网络。其中， 每个弧的容量为管道的输气能力，可以通过水力、热力分析模型计算。水力、热力的模型一 般通过商用软件建立，方法中采用TGNET软件为建模基础。将管道的管径、管长、压缩机配 置、天然气组成输入软件，即可以获取管道内的压力P和输气能力c。不同管道的输气能力不 同；同一管道在不同运行状态(正常或失效）的输气能力也不相同。
[0109] 基于图论可以对天然气管网进行抽象化处理，例如将其转化为容量网络[Clj]。具 体方法是根据天然气管网的实际拓扑结构，建立一个抽象的网络图。该网络图中，容量网络 包括多个节点和各节点之间的弧;容量网络中的节点为各天然气站场单元;各节点之间的 弧为各管道单元。然后，给图中的每个弧赋一个权值 Clj，代表节点i到节点j的管道输气能 力。
[0110] 此处，各管道单元或者天然气站场单元的失效会导致整个天然气管网的拓扑结构 改变，从而影响供气能力。
[0111] 步骤203、根据各管道单元信息确定各管道单元在各预设失效模式下的第一失效 概率。
[0112] 此处，由于天然气管道地域跨度广、影响因素复杂多变，对于不同区域的管道不能 等同确定其失效概率，因此需基于失效机理，分析不同管段的失效模式与失效概率。不同的 失效模式可以有管道腐蚀等情况。
[0113] 此处，当预设失效模式为管道腐蚀时，该步骤203可以通过如下方式实现：
[0114] 由于管道腐蚀导致的老化是管道失效的主要原因。腐蚀可能导致管体出现缺陷， 提高了失效风险。然而，管体的老化并不等同于管道一定会失效。管道失效是由多方面因素 决定的，如管壁厚度、管材强度、运行压力等。这些因素之间的相互作用是非线性的，运用解 析方法很难求解。因此，此处可以以极限状态方程构建失效率计算模型，同时采用蒙特卡洛 模拟求解该模型。如下所述：
[0115]根据极限状态方程建立管道失效概率计算模型，并根据蒙特卡洛模拟方法计算所 述管道失效概率计算模型，获取缺陷运行年限为t时的管道上每个缺陷的失效概率Ff(t);
[0116]其中，Ff⑴=P[pf⑴K0] ;pf⑴为缺陷运行年限为t时的管道缺陷的失效压 力，单位为Mpa;Po为管道的实际运行压力，单位为MPa，其可以通过管道的水里、热力模拟获 取，但不仅局限于此；
；其中，M为管道膨胀系数，也称为Folias系数，
为管道的流变应力（即与屈服 强度相应的材料属性为管道公称直径，单位为mm;Wt为管道壁厚，单位为mm;L(t)为t时 刻管道上缺陷的长度，单位为1111]1;以1:)=以0)+1^1:;以0)为0时刻管道上缺陷的长度，单位为 mm; ri为管道轴向腐蚀速率，单位为mm/a。
[0117] 根据齐次泊松模型，确定每千米管道上的缺陷数目k。
[0118] 其中，该齐次泊松模型为：
％n = 0，l，2…k;其中，Pr(n)为每千米管道 上出现n个缺陷的概率;k为每千米管道上的缺陷数目；x为根据内检测结果统计得到的每千 米管道上的缺陷数目的平均值。
[0119] 根据公式:A(t)=k ? Ff(t)确定一管道上在缺陷运行年限为t时的失效概率A(t)。
[0120] 步骤204、根据第一失效概率确定各管道单元在一预设时间内的第一状态分布。
[0121 ]第一状态分布包括各管道单元处于正常状态或各预设失效模式状态。
[0122] 此处，该步骤204可以通过如下方式实现：
[0123] 根据公式
，1/ = 0，1，2…，确定各管道单元在一预设时 间0至t内发生故障次数为Y的概率;其中，PKWOir/]为各管道单元在一预设时间0至t 内发生故障次数为Y的概率；;L为管道单元长度;A(t)为管道上在缺陷运 行年限为t时的失效概率。
[0124] 根据公式：p丨-丨j ^ 确定各管道单元从一to时刻开始下一个 故障出现的时间大于t时间的概率;Pr[Y(to)>t]为各管道单元从一 tQ时刻开始下一个故障 出现的时间大于t时间的概率;Y(to)为从一 to时刻开始下一个故障出现的时间。
[0125] 根据PHWOir/]和Pr[Y(t0)>t]确定各管道单元在所述预设时间内的第一状 态分布。
[0126] 步骤205、根据各天然气站场单元的历史失效情况数据，确定各天然气站场单元的 第二失效概率。
[0127] 此处，由于天然气站场单元的运行环境较为单一，因此可以采用历史失效情况数 据来确定其失效概率。
[0128] 步骤206、根据第二失效概率确定各天然气站场单元在一预设时间内的第二状态 分布。
[0129] 该步骤206可以通过如下方式实现：
[0130] 根据公式：
，：=0,1，2…，确定各天然气站场单元在 一预设时间〇至t内发生故障次数为Y的概率;其中，PV [Y (0=!/ ]为各天然气站场单元 在一预设时间〇至t内发生故障次数为Y的概率；W(/)=J^/'(/)dt;f(t)为各天然气站场单元 的第二失效概率。
[0131 ]根据公式：pj^yd 确定各天然气站场单元从一to时刻开始下 一个故障出现的时间大于t时间的概率;PV [疒（to)>t]为各天然气站场单元从一to时刻开 始下一个故障出现的时间大于t时间的概率;V (to)为从一to时刻开始下一个故障出现的时 间。
[0132] 根据PV [Y (0=!/ ]和PV [Y' (to)>t]确定各天然气站场单元在所述预设时间 内的第二状态分布。
[0133] 其中，第二状态分布可以包括各天然气站场单元处于正常状态或失效状态。
[0134] 步骤207、根据第一状态分布确定一预设时间点各管道单元的状态。
[0135] 步骤208、根据预设时间点各管道单元的状态，在所述对应关系中确定各管道单元 在正常状态和各预设失效模式状态下的输气能力，并将输气能力确定为容量网络的各弧容 量。
[0136] 步骤209、根据各弧容量和第二状态分布，实时更新容量网络。
[0137] 步骤210、根据实时更新的容量网络向天然气管网中的各用户供气，确定各用户的 供气情况信息。
[0138] 该各用户的供气情况信息包括：向用户模拟的供气次数N;向用户模拟的供气次数 中供气不足的次数K;向用户模拟的供气年份Y;每次向用户模拟的供气不足造成的缺气量 Z〇
[0139] 在向天然气管网中的各用户供气时可以采用就近供气的原则，即采用最短路径算 法来确定每一个随机生成的天然气管网容量网络的流向，即确定天然气在管网中的流动方 向。步骤可以如下：（1)建立天然气管网的管道长度矩阵。构建方法是在邻接矩阵的基础上， 将每条管道所对应的元素赋值为管道长度。（2)采用Bellman-Ford算法计算出每个气源到 所有用户的最近距离，并对这些距离进行升序排序；（3)根据排序结果，由低到高以依次用 户供气。经过该算法的调整，原容量矩阵[ Clj]中的元素排列方式发生改变，命名为[Clj]ne3W。 在容量矩阵中，元素的排列方式代表了天然气的流动方向。
[0140]之后，对于每个确定供气方向后的容量网络[cu]ne5W，都存在一个最大供气能力 Fmax。方法采用图论中的Ford-Fulkerson最大流算法计算每个容量网络[cij]ne5W最大供气能 力Fmax。将容量矩阵带入Ford-Fulkerson最大流算法就能够计算出该容量网络中，气源能够 向用户提供的最大天然气量总量以及每个用户能够获取的最大天然气量。这样就可以得到 一个充分的供气能力样本空间。该样本空间包括管网总供气能力[F max]与每个用户可获取 气量。之后可以采用统计方法，计算出管网总供气能力样本空间[Fmax]的累计概率分布函 数、期望、方差、中位数等统计特征。从而确定各用户的供气情况信息。
[0141] 步骤211、根据各用户的供气情况信息确定向各用户供气的可靠度以及各用户的 平均缺气量。
[0142] 此处，该步骤211可以通过如下方式实现：
[0143] 根据公式： 确定在年份Y向各用户供气的可靠度R。 >：
[0144] 根据公式：
，确定向一用户模拟的供气平均缺气量Q。
[0145] 为了使本领域的技术人员更好的了解本发明，下面列举一个应用上述步骤201至 步骤211的具体实例，如图3所示，其为某天然气管网的抽象图。该管网由管道、压气站场、管 道气源、地下储库及LNG站构成。其中，节点11和节点19为上游管道来气输入点，供气能力分 别为31mcm/d和25mcm/d;节点51为地下储气库，供气能力为7.1mcm/d;节点10为LNG站，正常 供气能力为 4mcm/d，经改造可提高至 10mcm/d。节点5,6,7，10，13，17，18,21，25,26,28,30， 33,34,35,36,37,38,39,41，42,43,45,47,48,49,51，52为用气需求点，对天然气的需气量 分别为 1 ? 43，1 ? 57，1 ? 66，1 ? 46，4 ? 4，1 ? 54，0 ? 5，1 ? 5，1 ? 6，1 ? 8，2 ? 5，2，0 ? 8，0 ? 8，1，1，1 ? 74， 1.3，1，2，1.4,0.5，1.06，1.82,0.68，1.17,2,0.98(mcm/d)。用户总体需求为41 ? 21mcm/d。
[0146] 下表1为某天然气管网的基本信息，其中输气能力可以通过TGNET软件模拟算得。 其中压缩机信息与天然气组成选择软件的默认设置。每个站场均假设为两台压缩机并联。 From表示输气的起始节点，To表示输气的接收点，Capacity表示输气量，Length表示输气长 度。
[0147] 表1:
[0149] 目前，管道内检测技术已比较成熟，通过分析检测数据，可以对管道缺陷的数量、 位置及其几何形状做出估计。表2是根据该管网的管道内检测结果获取的管体缺陷几何参 数统计特征。
[0150] 表2:
[0151]
[0152] 天然气管网由站场、管道、气源(上游管道、LNG站，地下储气库)构成以及天然气用 户构成。受环境因素、第三方活动、维护检修等影响，站场与管道为运行状态不确定的单元。 根据一般性认识，单元运行不稳定可能导致其存在多种运行状态，进而可以认为其输送天 然气的能力存在多种情况，其具体程度应按实际情况确定。因此，假设天然气管道存在三种 运行状态：正常运行、局部失效(输气能力下降至原输气能力的0.75,由TGNET算得）与完全 失效;采用基于失效机理的失效分析方法(蒙特卡罗模拟一百万次），得出在第5年局部失效 发生的概率是〇.〇〇〇l/(km ? a);完全失效的概率为3 X l(T5/(km ? a)。这里需要强调，在本模 型中，局部失效可以根据实际需要设定为多种状态；即只要有需要，任何单元可以存在任意 多个输气状态。根据历史统计数据，压气站的失效率为0.025。根据现场工作经验，发生失效 后，与其连接的管道输气能力下降至原输气能力的0.80。
[0153] 分别分析了原始情况、地下储气库失效情况、LNG站失效情况、地下储气库与LNG站 同时失效情况、提升LNG站供气能力这五种情况下的天然气管网供气能力统计特征与供气 可靠性。采用泊松过程模拟50万次。
[0154] 下表3表示供气能力统计分析与供气可靠性(mcm/d)。其中min表示天然气管网的 最小供气量；max表示天然气管网的最大供气量;mean表示天然气管网的平均供气量； median表示天然气管网的中间值;sta表示一致性参数;Supply reliabi 1 ity表示供其可靠 性。
[0155] 表3:

[0158] 进一步的，可以分析出整个管网在不确定性运行状态下的供气能力与供气可靠 性。如果给出总需求量，则可进一步计算其能够满足需求的概率。此外，可以看出，该模型可 以分析不同气源对管网供气能力的贡献;通过调整数据还可得知不同优化方案(增加气源、 提高起源供气能力、提高管输能力、降低管道失效率等)对提高供气安全性的作用大小。
[0159] 此外，还能够分析所有用户的需求满足程度。图4是管网中某用户(节点49)在不同 情况下能够获取的最大天然气量。可以看出，虽然地下储气库的失效对管网总体供气可靠 性行影响不大，但是对该用户影响比较明显。根据模型采用供气可靠度和平均缺气量两个 指标评价该天然气管网对于每个用户的供气可靠性，部分结果如下表4所示：
[0160] 表4:
[0162] 本发明实施例提供的一种天然气管网供气可靠性的评价方法，确定各管道单元和 天然气站场单元的失效概率，以及实时更新天然气管网的容量网络，并向天然气管网中的 各用户供气，确定各用户的供气情况信息；从而根据各用户的供气情况信息确定向各用户 供气的可靠度以及各用户的平均缺气量。通过本发明实施例提供的天然气管网供气可靠性 的评价方法，可以准确评价在不同干扰因素下天然气管网的供气可靠性，解决了当前天然 气管网的供气优劣难以确认的问题。
[0163] 对应于上述图1和图2所述的方法实施例，如图5所示，本发明实施例还提供一种天 然气管网供气可靠性的评价装置，包括：
[0164] 信息获取单元301，可以获取天然气管网的拓扑信息以及天然气管网拓扑中各管 道单元信息和各天然气站场单元信息。
[0165] 容量网络建立单元302,可以根据天然气管网的拓扑信息建立天然气管网的容量 网络;所述容量网络包括多个节点和各节点之间的弧;所述容量网络中的节点为各天然气 站场单元;各节点之间的弧为各管道单元。
[0166] 第一失效概率确定单元303,可以根据所述各管道单元信息确定各管道单元在各 预设失效模式下的第一失效概率。
[0167] 第一状态分布确定单元304,可以根据所述第一失效概率确定各管道单元在一预 设时间内的第一状态分布;所述第一状态分布包括各管道单元处于正常状态或各预设失效 模式状态。
[0168] 第二失效概率确定单元305,可以根据所述各天然气站场单元信息确定各天然气 站场单元的第二失效概率。
[0169] 第二状态分布确定单元306,可以根据所述第二失效概率确定各天然气站场单元 在一预设时间内的第二状态分布;所述第二状态分布包括各天然气站场单元处于正常状态 或失效状态。
[0170] 管道单元输气能力确定单元307,可以根据各管道单元信息以及所述第一状态分 布确定各管道单元在正常状态和各预设失效模式状态下的输气能力，并将所述输气能力确 定为所述容量网络的各弧容量。
[0171] 容量网络更新单元308,可以根据所述各弧容量和所述第二状态分布，实时更新所 述容量网络。
[0172] 模拟供气单元309,可以根据实时更新的容量网络向天然气管网中的各用户供气， 确定各用户的供气情况信息。
[0173] 可靠性评价单元310,可以根据所述各用户的供气情况信息确定向各用户供气的 可靠度以及各用户的平均缺气量。
[0174]此外，所述预设失效模式包括管道腐蚀。
[0175]该第一失效概率确定单元303,具体可以：
[0176]根据极限状态方程建立管道失效概率计算模型，并根据蒙特卡洛模拟方法计算所 述管道失效概率计算模型，获取缺陷运行年限为t时的管道上每个缺陷的失效概率Ff(t)。
[0177] 其中，Ff(t)=P[pf(t)-Po彡0] ;pf(t)为缺陷运行年限为t时的管道缺陷 的失效压力；Po为管道的实际运行压力；
其中，M为管道膨胀系数，
;of为管道的流变应力。D为管道 公称直径;Wt为管道壁厚;L(t)为t时刻管道上缺陷的长度;L(t)=L(0)+rLt;L(0)为0时刻 管道上缺陷的长度;^为管道轴向腐蚀速率。
[0178] 根据齐次泊松模型，确定每千米管道上的缺陷数目k。
[0179] 其中，该齐次泊松模型为
，n = 0，l，2~k;其中，Pr(n)为每千米管道 上出现n个缺陷的概率;k为每千米管道上的缺陷数目；x为根据内检测结果统计得到的每千 米管道上的缺陷数目的平均值。
[0180] 根据公式:A(t)=k ? Ff(t)确定一管道上在缺陷运行年限为t时的失效概率Mt)。 [0181 ]此外，该第一状态分布确定单元304，具体可以：
[0182]根据公式
…，确定各管道单元在一预设时 间〇至t内发生故障次数为Y的概率;其中，PKWOir/]为各管道单元在一预设时间0至t 内发生故障次数为Y的概率；妒(0=f'/-A/)dt :L为管道单元长度;A(t)为管道上在缺陷运 J〇 行年限为t时的失效概率。
[0183] 根据公式：G 确定各管道单元从一to时刻开始下一个 故障出现的时间大于t时间的概率;Pr[Y(to)>t]为各管道单元从一 to时刻开始下一个故障 出现的时间大于t时间的概率;Y(to)为从一 to时刻开始下一个故障出现的时间。
[0184] 根据PHWOir/]和Pr[Y(t0)>t]确定各管道单元在所述预设时间内的第一状 态分布。
[0185] 此外，该各天然气站场单元信息包括各天然气站场单元的历史失效情况数据。
[0186] 该第二失效概率确定单元305,具体可以：
[0187] 根据各天然气站场单元的历史失效情况数据，确定各天然气站场单元的第二失效 概率。
[0188] 此外，改第二状态分布确定单元306,具体可以：
[0189] 根据公式
，：=0,1，2…，确定各天然气站场单元在 一预设时间〇至t内发生故障次数为Y的概率;其中，PV [Y (0=!/ ]为各天然气站场单元 在一预设时间〇至t内发生故障次数为Y的概率；故⑴=f'/'(〇dt;f(t)为各天然气站场单元 Jo 的第二失效概率；
[0190] 根据公式：q 确定各天然气站场单元从一 to时刻开始下 一个故障出现的时间大于t时间的概率;PV [V (to)>t]为各天然气站场单元从一to时刻开 始下一个故障出现的时间大于t时间的概率;V (to)为从一to时刻开始下一个故障出现的时 间。
[0191]根据PV [Y (0=!/ ]和PV [Y' (to)>t]确定各天然气站场单元在所述预设时间 内的第二状态分布。
[0192] 另外，该各管道单元信息可以包括各管道的状态与输气能力的对应关系。
[0193] 该管道单元输气能力确定单元307，具体可以：
[0194] 根据第一状态分布确定一预设时间点各管道单元的状态。
[0195] 根据预设时间点各管道单元的状态，在所述对应关系中确定各管道单元在正常状 态和各预设失效模式状态下的输气能力。
[0196] 具体的，该各用户的供气情况信息包括：向用户模拟的供气次数N;向用户模拟的 供气次数中供气不足的次数K;向用户模拟的供气年份Y;每次向用户模拟的供气不足造成 的缺气量Z。
[0197] 该可靠性评价单元310,具体可以根据公式： 确定在年份Y向各用 ，+ 户供气的可靠度R;
[0198] 根据公式：
，确定向一用户模拟的供气平均缺气量Q。
[0199] 值得说明的是，本发明实施例提供的天然气管网供气可靠性的评价装置的具体实 现方式可以参见上述的方法实施例，此处不再赘述。
[0200] 本发明实施例提供的天然气管网供气可靠性的评价装置，确定各管道单元和天然 气站场单元的失效概率，以及实时更新天然气管网的容量网络，并向天然气管网中的各用 户供气，确定各用户的供气情况信息；从而根据各用户的供气情况信息确定向各用户供气 的可靠度以及各用户的平均缺气量。通过本发明实施例提供的天然气管网供气可靠性的评 价装置，可以准确评价在不同干扰因素下天然气管网的供气可靠性，解决了当前天然气管 网的供气优劣难以确认的问题。
[0201] 本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序 产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实 施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机 可用存储介质（包括但不限于磁盘存储器、CD-ROM、光学存储器等）上实施的计算机程序产 品的形式。
[0202] 本发明是参照根据本发明实施例的方法、设备（系统）、和计算机程序产品的流程 图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流 程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序 指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产 生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实 现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0203] 这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特 定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指 令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或 多个方框中指定的功能。
[0204] 这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计 算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或 其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一 个方框或多个方框中指定的功能的步骤。
[0205] 本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例 的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员， 依据本发明的思想，在【具体实施方式】及应用范围上均会有改变之处，综上所述，本说明书内 容不应理解为对本发明的限制。
【主权项】
1. 一种天然气管网供气可靠性的评价方法，其特征在于，包括： 获取天然气管网的拓扑信息W及天然气管网拓扑中各管道单元信息和各天然气站场 单元信息； 根据天然气管网的拓扑信息建立天然气管网的容量网络;所述容量网络包括多个节点 和各节点之间的弧;所述容量网络中的节点为各天然气站场单元;各节点之间的弧为各管 道单元； 根据所述各管道单元信息确定各管道单元在各预设失效模式下的第一失效概率； 根据所述第一失效概率确定各管道单元在一预设时间内的第一状态分布;所述第一状 态分布包括各管道单元处于正常状态或各预设失效模式状态； 根据所述各天然气站场单元信息确定各天然气站场单元的第二失效概率； 根据所述第二失效概率确定各天然气站场单元在一预设时间内的第二状态分布;所述 第二状态分布包括各天然气站场单元处于正常状态或失效状态； 根据各管道单元信息W及所述第一状态分布确定各管道单元在正常状态和各预设失 效模式状态下的输气能力，并将所述输气能力确定为所述容量网络的各弧容量； 根据所述各弧容量和所述第二状态分布，实时更新所述容量网络； 根据实时更新的容量网络向天然气管网中的各用户供气，确定各用户的供气情况信 息； 根据所述各用户的供气情况信息确定向各用户供气的可靠度W及各用户的平均缺气 量。2. 根据权利要求1所述的天然气管网供气可靠性的评价方法，其特征在于，所述预设失 效模式包括管道腐蚀； 根据所述各管道单元信息确定各管道单元在各预设失效模式下的第一失效概率，包 括： 根据极限状态方程建立管道失效概率计算模型，并根据蒙特卡洛模拟方法计算所述管 道失效概率计算模型，获取缺陷运行年限为t时的管道上每个缺陷的失效概率Ff (t); 其中，。^*)=?[9^*)斗。《〇1.^(*)^缺陷运行年限为*时的管道缺陷的失效压 力；Po为管道的实际运行压力；其中，M为管道膨胀系数，;0f为管道的流变应力； D为管道公称直径;Wt为管道壁厚;L(t)为t时刻管道上缺陷的长度;L(t)=L(0)+rLt;L (0)为加寸刻管道上缺陷的长度;a为管道轴向腐蚀速率； 根据齐次泊松模型，确定每千米管道上的缺陷数目k; 其中，该齐次泊松模型为：n = 0，l，2…k;其中，Pr(n)为每千米管道上出 现n个缺陷的概率;k为每千米管道上的缺陷数目；X为根据内检测结果统计得到的每千米管 道上的缺陷数目的平均值； 根据公式:A(t)=k ? Ff(t)确定一管道上在缺陷运行年限为t时的失效概率A(t)。3. 根据权利要求2所述的天然气管网供气可靠性的评价方法，其特征在于，根据所述第 一失效概率确定各管道单元在一预设时间内的第一状态分布，包括： 根据公式确定各管道单元在一预设时间O至 t内发生故障次数为n/的概率;其中，Pr阳(t)=n/]为各管道单元在一预设时间O至t内发生 故障次数为n/的概率L为管道单元长度;Mt)为管道上在缺陷运行年限 为t时的失效概率；根据公式确定各管道单元从一 to时刻开始下一个故障 出现的时间大于t时间的概率;Pr[Y(to)>t]为各管道单元从一 to时刻开始下一个故障出现 的时间大于t时间的概率;Y(to)为从一 to时刻开始下一个故障出现的时间； 根据Pr[N(t)=n/]和Pr[Y(t〇)>t]确定各管道单元在所述预设时间内的第一状态分 布。4. 根据权利要求3所述的天然气管网供气可靠性的评价方法，其特征在于，所述各天然 气站场单元信息包括各天然气站场单元的历史失效情况数据； 根据所述各天然气站场单元信息确定各天然气站场单元的第二失效概率，包括： 根据各天然气站场单元的历史失效情况数据，确定各天然气站场单元的第二失效概 率。5. 根据权利要求4所述的天然气管网供气可靠性的评价方法，其特征在于，根据所述第 二失效概率确定各天然气站场单元在一预设时间内的第二状态分布，包括： 根据公式：确定各天然气站场单元在一预设 时间0至t内发生故障次数为n/的概率;其中，Pr/ [沪（t) =n/]为各天然气站场单元在一预 设时间0至t内发生故障次数为n/的概率：f(t)为各天然气站场单元的第二 失效概率； 根据公式确定各天然气站场单元从一 to时刻开始下一个 故障出现的时间大于t时间的概率;Pr/ [Y^ (to)>t]为各天然气站场单元从一to时刻开始下 一个故障出现的时间大于t时间的概率;(to)为从一to时刻开始下一个故障出现的时间； 根据Pr/阳/ (t)=n/]和Pr/ [r (to) >t]确定各天然气站场单元在所述预设时间内的第 二状态分布。6. 根据权利要求5所述的天然气管网供气可靠性的评价方法，其特征在于，所述各管道 单元信息包括各管道的状态与输气能力的对应关系； 根据各管道单元信息W及所述第一状态分布确定各管道单元在正常状态和各预设失 效模式状态下的输气能力，包括： 根据第一状态分布确定一预设时间点各管道单元的状态； 根据预设时间点各管道单元的状态，在所述对应关系中确定各管道单元在正常状态和 各预设失效模式状态下的输气能力。7. 根据权利要求6所述的天然气管网供气可靠性的评价方法，其特征在于，所述各用户 的供气情况信息包括：向用户模拟的供气次数N;向用户模拟的供气次数中供气不足的次数 K;向用户模拟的供气年份Y;每次向用户模拟的供气不足造成的缺气量Z;根据所述各用户的供气情况信息确定向各用户供气的可靠度W及各用户的平均缺气 量，包括： 根据公式 确定在年份Y向各用户供气的可靠度R; 根据公式 I 一用户模拟的供气平均缺气量Q。8. -种天然气管网供气可靠性的评价装置，其特征在于，包括： 信息获取单元，用于获取天然气管网的拓扑信息W及天然气管网拓扑中各管道单元信 息和各天然气站场单元信息； 容量网络建立单元，用于根据天然气管网的拓扑信息建立天然气管网的容量网络；所 述容量网络包括多个节点和各节点之间的弧;所述容量网络中的节点为各天然气站场单 元;各节点之间的弧为各管道单元； 第一失效概率确定单元，用于根据所述各管道单元信息确定各管道单元在各预设失效 模式下的第一失效概率； 第一状态分布确定单元，用于根据所述第一失效概率确定各管道单元在一预设时间内 的第一状态分布;所述第一状态分布包括各管道单元处于正常状态或各预设失效模式状 态； 第二失效概率确定单元，用于根据所述各天然气站场单元信息确定各天然气站场单元 的第二失效概率； 第二状态分布确定单元，用于根据所述第二失效概率确定各天然气站场单元在一预设 时间内的第二状态分布;所述第二状态分布包括各天然气站场单元处于正常状态或失效状 态； 管道单元输气能力确定单元，用于根据各管道单元信息W及所述第一状态分布确定各 管道单元在正常状态和各预设失效模式状态下的输气能力，并将所述输气能力确定为所述 容量网络的各弧容量； 容量网络更新单元，用于根据所述各弧容量和所述第二状态分布，实时更新所述容量 网络； 模拟供气单元，用于根据实时更新的容量网络向天然气管网中的各用户供气，确定各 用户的供气情况信息； 可靠性评价单元，用于根据所述各用户的供气情况信息确定向各用户供气的可靠度W 及各用户的平均缺气量。9. 根据权利要求8所述的天然气管网供气可靠性的评价装置，其特征在于，所述预设失 效模式包括管道腐蚀； 该第一失效概率确定单元，具体用于： 根据极限状态方程建立管道失效概率计算模型，并根据蒙特卡洛模拟方法计算所述管 道失效概率计算模型，获取缺陷运行年限为t时的管道上每个缺陷的失效概率Ff (t);其中，Ff(t)=P[pf(t)-Pp《n1?n*ft^?缺陷运行年限为t时的管道缺陷的失效压 力；Po为管道的实际运行压力 ；其中，M为管道膨胀系数，;〇f为管道的流变应力； D为管道公称直径;Wt为管道壁厚;L(t)为t时刻管道上缺陷的长度;L(t)=L(0)+rLt;L (0)为加寸刻管道上缺陷的长度;为管道轴向腐蚀速率； 根据齐次泊松模型，确定每千米管道上的缺陷数目k; 其中，该齐次泊松模型为n = 0，l，2…k;其中，Pr(n)为每千米管道上出 现n个缺陷的概率;k为每千米管道上的缺陷数目；X为根据内检测结果统计得到的每千米管 道上的缺陷数目的平均值； 根据公式:A(t)=k ? Ff(t)确定一管道上在缺陷运行年限为t时的失效概率A(t)。10. 根据权利要求9所述的天然气管网供气可靠性的评价装置，其特征在于，所述第一 状态分布确定单元，具体用于： 根据公式?，确定各管道单元在一预设时间0至 t内发生故障次数为n/的概率:巧中，PWN(t)=n/]为各管道单元在一预设时间0至t内发生 故障次数为n/的概率L为管道单元长度;A(t)为管道上在缺陷运行年限 为t时的失效概率； 根据公式插定各管道单元从一 to时刻开始下一个故障 出现的时间大于t时间的概率;Pr[Y(to)>t]为各管道单元从一 to时刻开始下一个故障出现 的时间大于t时间的概率;Y(to)为从一 to时刻开始下一个故障出现的时间； 根据Pr[N(t)=n/]和Pr[Y(t〇)>t]确定各管道单元在所述预设时间内的第一状态分 布。11. 根据权利要求10所述的天然气管网供气可靠性的评价装置，其特征在于，所述各天 然气站场单元信息包括各天然气站场单元的历史失效情况数据； 所述第二失效概率确定单元，具体用于： 根据各天然气站场单元的历史失效情况数据，确定各天然气站场单元的第二失效概 率。12. 根据权利要求11所述的天然气管网供气可靠性的评价装置，其特征在于，所述第二 状态分布确巧单元，具化巧干： 根据公式/= O，1，2…，确定各天然气站场单元在一预设 时间O至t内发生故障次数为n/的概率;其中，Pr/ [沪（t) =n/]为各天然气站场单元在一预 设时间O至t内发生故障次数为n/的概率；（t)为各天然气站场单元的第二 失效概率；根据公式确定各天然气站场单元从一 to时刻开始下一个 故障出现的时间大于t时间的概率;Pr/ [Y^ (to)>t]为各天然气站场单元从一to时刻开始下 一个故障出现的时间大于t时间的概率;(to)为从一to时刻开始下一个故障出现的时间； 根据Pr/阳/ (t)=n/]和Pr/ [r (to) >t]确定各天然气站场单元在所述预设时间内的第 二状态分布。13. 根据权利要求12所述的天然气管网供气可靠性的评价装置，其特征在于，所述各管 道单元信息包括各管道的状态与输气能力的对应关系； 该管道单元输气能力确定单元，具体用于： 根据第一状态分布确定一预设时间点各管道单元的状态； 根据预设时间点各管道单元的状态，在所述对应关系中确定各管道单元在正常状态和 各预设失效模式状态下的输气能力。14. 根据权利要求12所述的天然气管网供气可靠性的评价装置，其特征在于，所述各用 户的供气情况信息包括：向用户模拟的供气次数N;向用户模拟的供气次数中供气不足的次 数K;向用户模拟的供气年份Y;每次向用户模拟的供气不足造成的缺气量Z; 该可靠性评价单元，具体用于根据公式确定在年份Y向各用户供气的 可靠度R; 根据公：确定向一用户模拟的供气平均缺气量Q。
【文档编号】G06F19/00GK105912844SQ201610213385
【公开日】2016年8月31日
【申请日】2016年4月7日
【发明人】张劲军, 苏怀, 杨楠, 王馨艺, 张宗杰
【申请人】中国石油大学(北京)