CN111143760A - 一种油田树枝状注水管网可靠度评价方法 - Google Patents

一种油田树枝状注水管网可靠度评价方法 Download PDF

Info

Publication number
CN111143760A
CN111143760A CN201911154177.6A CN201911154177A CN111143760A CN 111143760 A CN111143760 A CN 111143760A CN 201911154177 A CN201911154177 A CN 201911154177A CN 111143760 A CN111143760 A CN 111143760A
Authority
CN
China
Prior art keywords
node
pipe network
pipe
water injection
flow
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911154177.6A
Other languages
English (en)
Other versions
CN111143760B (zh
Inventor
于世春
杨仓海
李永长
李永清
李明江
黄天虎
李珍
闫娟
吴利利
赵晓龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201911154177.6A priority Critical patent/CN111143760B/zh
Publication of CN111143760A publication Critical patent/CN111143760A/zh
Application granted granted Critical
Publication of CN111143760B publication Critical patent/CN111143760B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/20—Administration of product repair or maintenance
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Mining

Abstract

本发明提供一种油田树枝状注水管网可靠度评价方法，其特征在于：包括以下步骤：步骤一：确定注水管网G0的连接方式；节点为单井、注水站或管网分支点，边是连接节点的管段，根据边与节点关系，绘制出整个管网系统的拓扑图；步骤二：确定管网各节点管段内流量q(i)分布；步骤三：确定各管段的故障流量的阈级T(ek)；步骤四：确定各管段基于阈级的单元重要度Iprob；步骤五：确定确定各节点i与源节点r联通的概率Pc(i)；步骤六：通过管段单元重要度Iprob与该管段节点起点与源节点联通概率乘积Pc(i)，累加计算整个树枝管网系统的可靠度R(T)；R(T)值越高，注水管网越安全可靠，本发明可快速完成树枝注水系统管网的整个可靠性计算，并根据计算结果进行管网评价，指导管网的运行与管理。

Description

一种油田树枝状注水管网可靠度评价方法
技术领域
本发明涉及油田注水开发及地面工艺管理领域，尤其涉及一种油田树枝状注 水管网可靠度评价方法。
背景技术
注水是保持油藏地层能量，维持油田高效经济开发最有效的手段。而注水管 网系统正常运行是油田注水开发的基本保障。目前各油田常用的注水管网有树枝 状管网与环状管网。长庆油田作为一个“低渗、低压、低产”的三低油气田，因 树枝状注水管网较环状管网投资少、建设周期短、管网分析简单、计算精度高等 特点，在长庆油田得到大面积应用。但树枝状注水管网系统一旦出现故障，会对 故障节点下游的注水系统造成影响，因此保障整个注水管网系统可靠性是油田注 水管理重要内容。评价注水管网系统的可靠性，一般采用管网可靠度量化计算， 即管网系统在规定的时间内，在目前条件下完成预定功能的能力的概率。
目前国内外计算可靠度的方法条件概率法、布尔法、马尔可夫方法、关联集 合分析法、故障树分析法、机理建模法、信息熵分析法、经验推理法等，各种方 法有各自的优缺点及适用条件，如概率法主要受限于不同单元失效或故障的定义 不统一，故障数据少影响失效函数建模的精度等；机理建模法多用于静态单元可 靠度评估；经验推理法在一定程度上可解决概率法需要的大样本问题，确定科学 的先验概率分布函数；马尔可夫方法在处理单元的共因失效和相关失效方面则更 具优势等。同时目前研究的可靠度计算方法在油气管道应用较多，对油田注水管 网鲜有应用。因此，建立一种适用性强、操作简单定量的油田树枝状注水管网系 统可靠度计算方法，指导现场的管网平稳安全运行管理，仍是亟待解决的技术难 题。
发明内容
为了克服现有注水管网系统一旦出现故障，会对故障节点下游的注水系统造 成影响的为题，本发明提供一种油田树枝状注水管网可靠度评价方法。该方法主 要针对油田注水树枝状管网系统，重点在于通过节点分析，遵循最大流－最小割 定理，其主要思路将注水管网系统，主要包括注水站、管道和注水井，看作一棵 以注水站为源点的有向树，以连接注水井的管道为枝，树中任意两个节点存在唯 一的路径，将树中每条枝所要输送的水量赋予该边的权，利用网络流理论与概率 论，分别求出节点最大流可靠度及基于阈级单元重要度与节点连通根源点的概 率，从而实现注水管网系统的可靠度计算与评价。其管网可靠度即注水管网所有 单元均能正常工作的概率，表示为P＝R{一定流量的水从站以规定压力到达各节 点}。
本发明采用的技术方案为：
一种油田树枝状注水管网可靠度评价方法，包括以下步骤：
步骤一：确定注水管网G0的连接方式；节点为单井、注水站或管网分支点， 边是连接节点的管段，根据边与节点关系，绘制出整个管网系统的拓扑图；
步骤二：确定管网各节点管段内流量q(i)分布；
步骤三：确定各管段的故障流量的阈级T(ek)；
步骤四：确定各管段基于阈级的单元重要度Iprob；
步骤五：确定确定各节点i与源节点r联通的概率Pc(i)；
步骤六：通过管段单元重要度Iprob与该管段节点起点与源节点联通概率乘积 Pc(i)，累加计算整个树枝管网系统的可靠度R(T)；
步骤七：根据R(T)的大小，进行管网系统设计评估、运行风险评估；R(T) 值越高，注水管网越安全可靠。
所述的步骤二中，确定管网各节点管段内流量q(i)分布，具体步骤为：
第一，根据管网拓朴图，建立管网节点与管段关联矩阵，其元素aij(i＝1～n， n表示节点；j＝1～m，m表示管段)，且规定
Figure BDA0002284353320000021
从而确定关联矩阵A为根据附图2(6个节点7个管段，且均已编号，节5 和6为已知压力节点)：
Figure BDA0002284353320000031
第二，建立节点流量方程组，其节点连续性方程，节点流出为正，流入为负；
Figure BDA0002284353320000032
式中：Qi—管段i的流量，m3/s；qi—节点i的流量，m3/s。
第三，管网的管段流量和压力降关系用矩阵G：
建立式子(2)和(3)的矩阵表达式：
AQ＝q (4)
式中：Q＝(Q1,Q2,···,Qm)T为管段流量向量；q＝(q1,q2,···,qn)T为节点载荷向量。
令p＝(p1,p2,···,pn)T为节点相对压力向量，则管段j上游压力与下游压力差Δp：
Δpj＝SjQj a＝Sj|Qj|a-1Qj (5)
式中：Sj为管段阻抗。此式(5)可以写成：
Figure BDA0002284353320000033
其中gj称为管段j的导纳。
则管网的管段流量和压力降关系用矩阵可表示为：
Q＝GΔp (6)
其中G为管段的导纳矩阵：
由式(4)、(5)和(6)可得：
AGATp＝q (7)
式(7)共n个方程，2n个水力参数，其中n个节点压力和n个节点载荷；
第四，水力学模型求解
对管网中编号为1～(n-h)的节点，假设其压力向量p1是未知量，载荷 向量q1是已知量；而对编号为(n-h+1)～n的节点，假设其压力向量p2是 已知量，载荷向量q2是未知量，令Y＝AGAT，Y是以主对角线为中心的对称 矩阵。对角元素yij，等于与节点i有关的管段导纳之和，其他非零元素yij，等 于管段(i，j)的导纳的负值。如果节点i和j之间没有管段连接，则yij＝0。
则公式(7)可表示成块矩阵的形式：
Figure BDA0002284353320000041
从而得到：
y11p1＝q1-y12p2 (9)
q2＝y21p1+y22p2 (10)
从方程组(9)和(10)可分别求出未知压力向量p1和未知载荷向量q2；然后将 整个管网节点相对压力向量p代入式(5)和(6)可求各管段压力降和流量。
所述的第四，水力学模型求解中，节点编号为1～(n-h)的节点具有已知载荷 和未知压力，而编号为(n-h+1)～n的节点具有已知压力和未知载荷。
所述的步骤三中，确定各管段的故障流量的阈级T(ek)的具体步骤为：记ψc为整个传输时间内所能传输流的总量，ψmax为起点到各终点所能传输的最大流 (端容量)，故障流量ψ(e1,e2,···,ek)为管段e1,e2,···,ek被删除后能传输到井(终点) 的最大流量；
流可靠度R(ψc)为在时间间隔ΔT内，恰有边e1,e2,···,en失效时总流的期望值<ψc>与无失效情况发生的总流ψc的比值；
Figure BDA0002284353320000042
式中1≤k≤n且k＝1,2,···,n；F(e1,···,ek)为注水管网G0中恰有管段e1,e2,···,en失效的概率；ψc(e1,e2,···,en)为总流ψc在管段e1,e2,···,en失效时所能 传输的最大流量；
令
Figure BDA0002284353320000051
当ψc＝ψmax,c，式(11)最大流可靠度可转换为：
Figure BDA0002284353320000052
当ψc＜ψmax,c时，最大流可靠度：
Figure BDA0002284353320000053
T(ek)即为管段的故障流量的阈级。
所述的步骤四中确定各管段基于阈级的单元重要度Iprob；
Figure BDA0002284353320000054
式中：T(ek)为恰有边ek失效时的阈级。
所述的步骤五中，确定各节点i与源节点r联通的概率Pc(i)；
Pc(i)＝PiPmPc(n)*100％ (15)
其中m(n，i)是树中一条链,节点n位于从节点i到根r的路径上，则 称n是i的前启，i是n的后继；Pi为节点i正常工作的概率；Pm为链m(管段) 正常工作的概率，Pc(n)为节点i的前启n与根r连通的概率。
所述的步骤六中，通过管段单元重要度Iprob与该管段节点起点与源节点联通 概率乘积Pc(i)，累加计算整个树枝管网系统的可靠度；
Figure BDA0002284353320000055
式中：ωi为注水管网中以节点i为起点基于阈级的单元重要度，
Figure BDA0002284353320000056
结果R(T)值越高，说明注水管网越安全可靠，从而根据R(T)的大小，进行 管网系统设计评估、运行风险评估。
本发明的有益效果为：
本发明利用注水管网运行数据，定期定量描述评估管网运行风险，及时掌握 管网系统的薄弱环节，采取措施消减管网运行风险，确保管网系统安全、平稳运 行，同时为下一步管网系统优化及建立管网增强方案提供技术支撑。
本发明充分利用油田物联网建设成果，将注水系统实时采集数据与静态数据 相结合，深化数据分析与应用，实现了注水管网可靠性评价分析。
本发明方法可快速完成树枝注水系统管网的整个可靠性计算，并根据计算结 果进行管网评价，指导管网的运行与管理。
本发明在管网运行时，可定期评估管网运行风险，制定维护计划、管理考核 等，掌握管网系统的薄弱环节，提出增强方案，制定应急预案，保证维护资金的 优化配置，实现管网的安全平稳运行。
以下将结合附图进行进一步的说明。
附图说明
图1为水力学计算方法管网拓扑结构说明示意图。
图2为姬二注水管网的拓扑结构图。
具体实施方式
实施例1：
为了克服现有注水管网系统一旦出现故障，会对故障节点下游的注水系统造 成影响的为题，本发明提供如图1和2所示的一种油田树枝状注水管网可靠度评 价方法。该方法主要针对油田注水树枝状管网系统，重点在于通过节点分析，遵 循最大流－最小割定理，其主要思路将注水管网系统，主要包括注水站、管道和 注水井，看作一棵以注水站为源点的有向树，以连接注水井的管道为枝，树中任 意两个节点存在唯一的路径，将树中每条枝所要输送的水量赋予该边的权，利用 网络流理论与概率论，分别求出节点最大流可靠度及基于阈级单元重要度与节点 连通根源点的概率，从而实现注水管网系统的可靠度计算与评价。其管网可靠度 即注水管网所有单元均能正常工作的概率，表示为P＝R{一定流量的水从站以规 定压力到达各节点}。
一种油田树枝状注水管网可靠度评价方法，包括以下步骤：
步骤一：确定注水管网G0的连接方式；节点为单井、注水站或管网分支点， 边是连接节点的管段，根据边与节点关系，绘制出整个管网系统的拓扑图；
步骤二：确定管网各节点管段内流量q(i)分布；
步骤三：确定各管段的故障流量的阈级T(ek)；
步骤四：确定各管段基于阈级的单元重要度Iprob；
步骤五：确定确定各节点i与源节点r联通的概率Pc(i)；
步骤六：通过管段单元重要度Iprob与该管段节点起点与源节点联通概率乘积 Pc(i)，累加计算整个树枝管网系统的可靠度R(T)；
步骤七：根据R(T)的大小，进行管网系统设计评估、运行风险评估；R(T) 值越高，注水管网越安全可靠。
本发明利用注水管网运行数据，定期定量描述评估管网运行风险，及时掌握 管网系统的薄弱环节，采取措施消减管网运行风险，确保管网系统安全、平稳运 行，同时为下一步管网系统优化及建立管网增强方案提供技术支撑。本发明充分 利用油田物联网建设成果，将注水系统实时采集数据与静态数据相结合，深化数 据分析与应用，实现了注水管网可靠性评价分析。本发明方法可快速完成树枝注 水系统管网的整个可靠性计算，并根据计算结果进行管网评价，指导管网的运行 与管理。
本发明在管网运行时，可定期评估管网运行风险，制定维护计划、管理考核 等，掌握管网系统的薄弱环节，提出增强方案，制定应急预案，保证维护资金的 优化配置，实现管网的安全平稳运行。
实施例2：
基于实施例1的基础上，本实施例中，所述的步骤二中，确定管网各节点管 段内流量q(i)分布，具体步骤为：第一，根据管网拓朴图，建立管网节点与管 段关联矩阵，其元素aij(i＝1～n，n表示节点；j＝1～m，m表示管段)，且规定
Figure BDA0002284353320000071
从而确定关联矩阵A为根据附图2(6个节点7个管段，且均已编号，节5 和6为已知压力节点)：
Figure BDA0002284353320000081
第二，建立节点流量方程组，其节点连续性方程，节点流出为正，流入为负；
Figure BDA0002284353320000082
式中：Qi—管段i的流量，m3/s；qi—节点i的流量，m3/s。
第三，管网的管段流量和压力降关系用矩阵G：
建立式子(2)和(3)的矩阵表达式：
AQ＝q (4)
式中：Q＝(Q1,Q2,···,Qm)T为管段流量向量；q＝(q1,q2,···,qn)T为节点载荷向量。
令p＝(p1,p2,···,pn)T为节点相对压力向量，则管段j上游压力与下游压力差Δp：
Δpj＝SjQj a＝Sj|Qj|a-1Qj (5)
式中：Sj为管段阻抗。此式(5)可以写成：
Figure BDA0002284353320000083
其中gj称为管段j的导纳。
则管网的管段流量和压力降关系用矩阵可表示为：
Q＝GΔp (6)
其中G为管段的导纳矩阵：
由式(4)、(5)和(6)可得：
AGATp＝q (7)
式(7)共n个方程，2n个水力参数，其中n个节点压力和n个节点载荷；
第四，水力学模型求解
对管网中编号为1～(n-h)的节点，假设其压力向量p1是未知量，载荷 向量q1是已知量；而对编号为(n-h+1)～n的节点，假设其压力向量p2是 已知量，载荷向量q2是未知量，令Y＝AGAT，Y是以主对角线为中心的对称 矩阵。对角元素yij，等于与节点i有关的管段导纳之和，其他非零元素yij，等 于管段(i，j)的导纳的负值。如果节点i和j之间没有管段连接，则yij＝0。
则公式(7)可表示成块矩阵的形式：
Figure BDA0002284353320000091
从而得到：
y11p1＝q1-y12p2 (9)
q2＝y21p1+y22p2 (10)
从方程组(9)和(10)可分别求出未知压力向量p1和未知载荷向量q2；然后将 整个管网节点相对压力向量p代入式(5)和(6)可求各管段压力降和流量。
所述的第四，水力学模型求解中，节点编号为1～(n-h)的节点具有已知载荷 和未知压力，而编号为(n-h+1)～n的节点具有已知压力和未知载荷。
所述的步骤三中，确定各管段的故障流量的阈级T(ek)的具体步骤为：记ψc为整个传输时间内所能传输流的总量，ψmax为起点到各终点所能传输的最大流 (端容量)，故障流量ψ(e1,e2,···,ek)为管段e1,e2,···,ek被删除后能传输到井(终点) 的最大流量；
流可靠度R(ψc)为在时间间隔ΔT内，恰有边e1,e2,···,en失效时总流的期望值<ψc>与无失效情况发生的总流ψc的比值；
Figure BDA0002284353320000092
式中1≤k≤n且k＝1,2,···,n；F(e1,···,ek)为注水管网G0中恰有管段e1,e2,···,en失效的概率；ψc(e1,e2,···,en)为总流ψc在管段e1,e2,···,en失效时所能 传输的最大流量；
令
Figure BDA0002284353320000101
当ψc＝ψmax,c，式(11)最大流可靠度可转换为：
Figure BDA0002284353320000102
当ψc＜ψmax,c时，最大流可靠度：
Figure BDA0002284353320000103
T(ek)即为管段的故障流量的阈级。
所述的步骤四中确定各管段基于阈级的单元重要度Iprob；
Figure BDA0002284353320000104
式中：T(ek)为恰有边ek失效时的阈级。
所述的步骤五中，确定各节点i与源节点r联通的概率Pc(i)；
Pc(i)＝PiPmPc(n)*100％ (15)
其中m(n，i)是树中一条链,节点n位于从节点i到根r的路径上，则 称n是i的前启，i是n的后继；Pi为节点i正常工作的概率；Pm为链m(管段) 正常工作的概率，Pc(n)为节点i的前启n与根r连通的概率。
所述的步骤六中，通过管段单元重要度Iprob与该管段节点起点与源节点联通 概率乘积Pc(i)，累加计算整个树枝管网系统的可靠度；
Figure BDA0002284353320000105
式中：ωi为注水管网中以节点i为起点基于阈级的单元重要度，
Figure BDA0002284353320000106
结果R(T)值越高，说明注水管网越安全可靠，从而根据R(T)的大小，进行 管网系统设计评估、运行风险评估。
本发明提供的一种树枝状注水管网系统可靠度确定方法。该方法主要针对油 田注水树枝状管网系统，重点在于通过节点分析，遵循最大流－最小割定理，其 主要思路将注水管网系统，主要包括注水站、管道和注水井，看作一棵以注水站 为源点的有向树，以连接注水井的管道为枝，树中任意两个节点存在唯一的路径， 将树中每条枝所要输送的水量赋予该边的权，利用网络流理论与概率论，分别求 出节点最大流可靠度及基于阈级单元重要度与节点连通根源点的概率，从而实现 注水管网系统的可靠度计算与评价。其管网可靠度即注水管网所有单元均能正常 工作的概率，表示为P＝R{一定流量的水从站以规定压力到达各节点}。
主要包括以下主要步骤：
步骤一：确定注水管网G0的连接方式。节点为单井、注水站或管网分支点 等，边是连接节点的管段，根据边与节点关系，绘制出整个管网系统的拓扑图。
步骤二：利用水力学计算，确定管网各节点管段内流量q(i)分布。
本发明提供了一种解节点连续方程的方式进行水力学计算，确定各节点的流 量分布，但不局限于这种方法。
1.根据管网拓朴图，建立管网节点与管段关联矩阵，其元素aij(i＝1～n，n 表示节点；j＝1～m，m表示管段)，且规定
Figure BDA0002284353320000111
从而确定关联矩阵A为根据附图2(6个节点7个管段，且均已编号，节5 和6为已知压力节点)：
Figure BDA0002284353320000112
2.建立节点流量方程组(节点连续性方程，节点流出为正，流入为负)；
Figure BDA0002284353320000121
式中：Qi—管段i的流量，m3/s；qi—节点i的流量，m3/s。
3.管网的管段流量和压力降关系用矩阵G：
建立(2)、(3)矩阵表达式：
AQ＝q (4)
式中：Q＝(Q1,Q2,···,Qm)T为管段流量向量；q＝(q1,q2,···,qn)T为节点载荷向量。
令p＝(p1,p2,···,pn)T为节点相对压力向量，则管段j上游压力与下游压力差Δp：
Δpj＝SjQj a＝Sj|Qj|a-1Qj (5)
式中：Sj为管段阻抗。此式(5)可以写成：
Figure BDA0002284353320000122
其中gj称为管段j的导纳。
则管网的管段流量和压力降关系用矩阵可表示为：
Q＝GΔp (6)
其中G为管段的导纳矩阵：
由式(4)、(5)和(6)可得：
AGATp＝q (7)
式(7)共n个方程，2n个水力参数(n个节点压力和n个节点载荷)。为求出 该方程组的唯一解，需要给定其中n个参数的值。主要做法是对每个节点给定其 一个对应参数(压力或载荷)的值，例如h个节点压力和n-h个节点载荷。通过不 断迭代修正管段流量的值，可求出各个节点的未知参数。
4.水力学模型求解(矩阵分块)。
对管网中编号为1～(n-h)的节点，假设其压力向量p1是未知量，载荷 向量q1是已知量；而对编号为(n-h+1)～n的节点，假设其压力向量p2是 已知量，载荷向量q2是未知量，令Y＝AGAT，Y是以主对角线为中心的对称 矩阵。对角元素yij，等于与节点i有关的管段导纳之和，其他非零元素yij，等 于管段(i，j)的导纳的负值。如果节点i和j之间没有管段连接，则yij＝0。
则公式(7)可表示成块矩阵的形式：
Figure BDA0002284353320000131
从而得到：
y11p1＝q1-y12p2 (9)
q2＝y21p1+y22p2 (10)
从方程组(9)和(10)可分别求出未知压力向量p1和未知载荷向量q2。进一步 将整个管网节点相对压力向量p代入式(5)和(6)可求各管段压力降和流量。
矩阵分块过程进行的前提条件是保证节点编号为1～(n-h)的节点具有已知 载荷和未知压力，而编号为(n-h+1)～n的节点具有已知压力和未知载荷。实际 上，如果与每个节点相关的2个参数(压力和载荷)中都刚好有一个已知量(压力 或载荷)，即使节点编号不符合上述规则，只需重新排列方程组(7)中方程的次序， 仍然可以采用上述分块矩阵的处理方法。
步骤三：计算各管段的故障流量的阈级T(ek)。
记ψc为整个传输时间内所能传输流的总量，ψmax为起点到各终点所能传输 的最大流(端容量)，故障流量ψ(e1,e2,···,ek)为管段e1,e2,···,ek被删除后能传输到 井(终点)的最大流量。
流可靠度R(ψc)为在时间间隔ΔT内，恰有边e1,e2,···,en失效时总流的期望值<ψc>与无失效情况发生的总流ψc的比值。
Figure BDA0002284353320000132
式中1≤k≤n且k＝1,2,···,n；F(e1,···,ek)为注水管网G0中恰有管段e1,e2,···,en失效的概率；ψc(e1,e2,···,en)为总流ψc在管段e1,e2,···,en失效时所能 传输的最大流量。
令
Figure BDA0002284353320000141
当ψc＝ψmax,c，式(11)最大流可靠度可转换为：
Figure BDA0002284353320000142
当ψc＜ψmax,c时，最大流可靠度：
Figure BDA0002284353320000143
T(ek)即为管段的故障流量的阈级。
步骤四：计算各管段基于阈级的单元重要度Iprob。
Figure BDA0002284353320000144
式中：T(ek)为恰有边ek失效时的阈级。
步骤五：计算各节点i与源节点r联通的概率Pc(i)。
Pc(i)＝PiPmPc(n)*100％ (15)
其中m(n，i)是树中一条链,节点n位于从节点i到根r的路径上，则 称n是i的前启，i是n的后继。Pi为节点i正常工作的概率；Pm为链m(管段) 正常工作的概率，Pc(n)为节点i的前启n与根r连通的概率。
步骤六：利用管段单元重要度Iprob与该管段节点起点与源节点联通概率乘积 Pc(i)，累加计算整个树枝管网系统的可靠度。
Figure BDA0002284353320000145
式中：ωi为注水管网中以节点i为起点基于阈级的单元重要度，
Figure BDA0002284353320000146
结果R(T)值越高，说明注水管网越安全可靠。从而根据R(T)的大小，进行 管网系统设计评估、运行风险评估等，进一步为树枝状注水管网系统优化及安全 运行提供指导。
本发明目的主要利用系统节点分析理论及网络流理论，提供一种树枝状注水 管网系统可靠度计算方法，用于油田注水管网系统评价管理，从而提升管网系统 管控水平及管理效率，确保管网的安全平稳运行。
实施例3：
基于实施例1和2的基础上，本实施例对长庆油田姬塬区块姬二注管网系统 进行了可靠度计算与评估。姬二注管网为典型树枝状管网系统，五条注水干线 (202、长45、长8、203配、耿63)、稳流配水阀组48个、水井129口，管线 规格包括Ф168*13、Ф114*11、Ф89*9、Ф76*8、Ф48*5。
根据现场实际情况，绘制了整个管网拓扑图，并进行了节点、管段进行了 编号；处理后整个管网拓扑结构包括普通节点41个，稳流配水阀组节点48 个，水源节点5个，管段89根。
根据本发明提供的水力学计算方法进行管网的水力学模拟计算，其中41个 普通节点为负载及压力均未知点，48个稳流配水阀组节点为负载已知节点，5 个水源节点为压力已知节点，具体求解结果见下表。
姬二注管网系统水力学模拟计算结果
Figure BDA0002284353320000151
Figure BDA0002284353320000161
Figure BDA0002284353320000171
根据水力学计算结果，利用本发明中提供的公式(11)～(14)计算各管段 的故障流量的阈级T(ek)和基于阈级的各管段单元重要度Iprob。具体值见表：
姬二注管段阈级及重要度计算结果
Figure BDA0002284353320000172
Figure BDA0002284353320000181
然后计算各节点与源节点的连通概率。一个给定大型网络系统中，节点与链 路正常工作的概率一般是一定的。为了计算简便，且考虑当节点与链路均不可靠 时，节点比链路更容易失效，这里给定Pi＝0.95，Pm＝0.99。从而利用各节点 的与源节点的连通概率公式(15)计算结果如下：
姬二注各节点与根节点的连通概率
Figure BDA0002284353320000191
最后根据公式(16)计算姬二注管网系统可靠度为0.7978。
根据计算结果可分析评价该管网的运行风险，并找到提升系统可靠度有效方 法，指导现场管理。
以上所述仅为本发明示意性的具体实施方式，并非用以限定本发明的范围。 任何本领域的技术人员，在不脱离本发明的构思和原则的前提下所作出的等同变 化与修改，均应属于本发明保护的范围。本发明中没有详细描述的方法及其步骤 均为现有技术，本发明中将不再进行一一叙述。

Claims (7)

1.一种油田树枝状注水管网可靠度评价方法，其特征在于：包括以下步骤：
步骤一：确定注水管网G0的连接方式；节点为单井、注水站或管网分支点，边是连接节点的管段，根据边与节点关系，绘制出整个管网系统的拓扑图；
步骤二：确定管网各节点管段内流量q(i)分布；
步骤三：确定各管段的故障流量的阈级T(ek)；
步骤四：确定各管段基于阈级的单元重要度Iprob；
步骤五：确定确定各节点i与源节点r联通的概率Pc(i)；
步骤六：通过管段单元重要度Iprob与该管段节点起点与源节点联通概率乘积Pc(i)，累加计算整个树枝管网系统的可靠度R(T)；
步骤七：根据R(T)的大小，进行管网系统设计评估、运行风险评估；R(T)值越高，注水管网越安全可靠。
2.根据权利要求1所述的一种油田树枝状注水管网可靠度评价方法，其特征在于：所述的步骤二中，确定管网各节点管段内流量q(i)分布，具体步骤为：
第一，根据管网拓朴图，建立管网节点与管段关联矩阵，其元素aij(i＝1～n，n表示节点；j＝1～m，m表示管段)，且规定
Figure FDA0002284353310000011
从而确定关联矩阵A为根据附图2(6个节点7个管段，且均已编号，节5和6为已知压力节点)：
Figure FDA0002284353310000012
第二，建立节点流量方程组，其节点连续性方程，节点流出为正，流入为负；
Figure FDA0002284353310000021
式中：Qi—管段i的流量，m3/s；qi—节点i的流量，m3/s。
第三，管网的管段流量和压力降关系用矩阵G：
建立式子(2)和(3)的矩阵表达式：
AQ＝q (4)
式中：Q＝(Q1,Q2,…,Qm)T为管段流量向量；q＝(q1,q2,…,qn)T为节点载荷向量。
令p＝(p1,p2,…,pn)T为节点相对压力向量，则管段j上游压力与下游压力差Δp：
Δpj＝SjQj a＝Sj|Qj|a-1Qj (5)
式中：Sj为管段阻抗。此式(5)可以写成：
Figure FDA0002284353310000022
其中gj称为管段j的导纳。
则管网的管段流量和压力降关系用矩阵可表示为：
Q＝GΔp (6)
其中G为管段的导纳矩阵：
由式(4)、(5)和(6)可得：
AGATp＝q (7)
式(7)共n个方程，2n个水力参数，其中n个节点压力和n个节点载荷；
第四，水力学模型求解
对管网中编号为1～(n-h)的节点，假设其压力向量p1是未知量，载荷向量q1是已知量；而对编号为(n-h+1)～n的节点，假设其压力向量p2是已知量，载荷向量q2是未知量，令Y＝AGAT，Y是以主对角线为中心的对称矩阵。对角元素yij，等于与节点i有关的管段导纳之和，其他非零元素yij，等于管段(i，j)的导纳的负值。如果节点i和j之间没有管段连接，则yij＝0。
则公式(7)可表示成块矩阵的形式：
Figure FDA0002284353310000031
从而得到：
y11p1＝q1-y12p2 (9)
q2＝y21p1+y22p2 (10)
从方程组(9)和(10)可分别求出未知压力向量p1和未知载荷向量q2；然后将整个管网节点相对压力向量p代入式(5)和(6)可求各管段压力降和流量。
3.根据权利要求2所述的一种油田树枝状注水管网可靠度评价方法，其特征在于：所述的第四，水力学模型求解中，节点编号为1～(n-h)的节点具有已知载荷和未知压力，而编号为(n-h+1)～n的节点具有已知压力和未知载荷。
4.根据权利要求1所述的一种油田树枝状注水管网可靠度评价方法，其特征在于：所述的步骤三中，确定各管段的故障流量的阈级T(ek)的具体步骤为：记ψc为整个传输时间内所能传输流的总量，ψmax为起点到各终点所能传输的最大流(端容量)，故障流量ψ(e1,e2,…,ek)为管段e1,e2,…,ek被删除后能传输到井(终点)的最大流量；
流可靠度R(ψc)为在时间间隔ΔT内，恰有边e1,e2,…,en失效时总流的期望值<ψc>与无失效情况发生的总流ψc的比值；
Figure FDA0002284353310000032
式中1≤k≤n且k＝1,2,…,n；F(e1,…,ek)为注水管网G0中恰有管段e1,e2,…,en失效的概率；ψc(e1,e2,…,en)为总流ψc在管段e1,e2,…,en失效时所能传输的最大流量；
令
Figure FDA0002284353310000041
当ψc＝ψmax,c，式(11)最大流可靠度可转换为：
Figure FDA0002284353310000042
当ψc＜ψmax,c时，最大流可靠度：
Figure FDA0002284353310000043
T(ek)即为管段的故障流量的阈级。
5.根据权利要求1所述的一种油田树枝状注水管网可靠度评价方法，其特征在于：所述的步骤四中确定各管段基于阈级的单元重要度Iprob；
Figure FDA0002284353310000044
式中：T(ek)为恰有边ek失效时的阈级。
6.根据权利要求1所述的一种油田树枝状注水管网可靠度评价方法，其特征在于：所述的步骤五中，确定各节点i与源节点r联通的概率Pc(i)；
Pc(i)＝PiPmPc(n)*100％ (15)
其中m(n，i)是树中一条链,节点n位于从节点i到根r的路径上，则称n是i的前启，i是n的后继；Pi为节点i正常工作的概率；Pm为链m(管段)正常工作的概率，Pc(n)为节点i的前启n与根r连通的概率。
7.根据权利要求1所述的一种油田树枝状注水管网可靠度评价方法，其特征在于：所述的步骤六中，通过管段单元重要度Iprob与该管段节点起点与源节点联通概率乘积Pc(i)，累加计算整个树枝管网系统的可靠度；
Figure FDA0002284353310000045
式中：ωi为注水管网中以节点i为起点基于阈级的单元重要度，
Figure FDA0002284353310000046
结果R(T)值越高，说明注水管网越安全可靠，从而根据R(T)的大小，进行管网系统设计评估、运行风险评估。
CN201911154177.6A 2019-11-22 2019-11-22 一种油田树枝状注水管网可靠度评价方法 Active CN111143760B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911154177.6A CN111143760B (zh) 2019-11-22 2019-11-22 一种油田树枝状注水管网可靠度评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911154177.6A CN111143760B (zh) 2019-11-22 2019-11-22 一种油田树枝状注水管网可靠度评价方法

Publications (2)

Publication Number Publication Date
CN111143760A true CN111143760A (zh) 2020-05-12
CN111143760B CN111143760B (zh) 2023-12-01

Family

ID=70517250

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911154177.6A Active CN111143760B (zh) 2019-11-22 2019-11-22 一种油田树枝状注水管网可靠度评价方法

Country Status (1)

Country Link
CN (1) CN111143760B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105721196A (zh) * 2016-01-20 2016-06-29 昆明理工大学 一种有向通信网络的链路重要性评价方法
JP2017069620A (ja) * 2015-09-28 2017-04-06 日本電信電話株式会社 信頼性評価装置、信頼性評価方法、及びプログラム

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017069620A (ja) * 2015-09-28 2017-04-06 日本電信電話株式会社 信頼性評価装置、信頼性評価方法、及びプログラム
CN105721196A (zh) * 2016-01-20 2016-06-29 昆明理工大学 一种有向通信网络的链路重要性评价方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王博弘;孟令尊;张意;张昕;张建军;王东旭;梁永图;: "某采油厂注水管网水力仿真模拟", 给水排水 *

Also Published As

Publication number Publication date
CN111143760B (zh) 2023-12-01

Similar Documents

Publication Publication Date Title
Bhave Node flow analysis distribution systems
CN107886182B (zh) 油田集输系统优化设计方法及装置
CN104466959B (zh) 电力系统关键线路辨识方法和系统
CN109344555B (zh) 一种集输管网设计方法及装置
US11326977B2 (en) Efficient method for localizing leaks in water supply pipe network based on valve operations and online water metering
CN106959608A (zh) 一种基于聚类粒子群算法的供水管网渗漏优化控制方法
CN107101087A (zh) 城市天然气管道系统管道故障传播影响评估方法及系统
CN113449993A (zh) 一种城市水源供水调度方法
CN109375050A (zh) 一种基于初始故障线路割集的电力系统关键断面辨识方法
CN111143760A (zh) 一种油田树枝状注水管网可靠度评价方法
CN107120531A (zh) 城市天然气管道系统故障传播影响评估方法及系统
Jose et al. Optimal water distribution network design accounting for valve closure
Vasilić et al. Network sectorisation through aggregation of strong connected components
CN105138753A (zh) 事故工况下高压燃气管网节点实际流量的确定方法
CN116561942B (zh) 用于城市排水系统拓扑数据的修正方法以及装置
Sabet et al. Determination of value and location of the leakage in water distribution networks using calibration (Case Study: The water distribution network of Golbahar, Khorasan Razavi)
CN115795122B (zh) 一种城市排水管网拓扑关系梳理方法
CN112597634B (zh) 一种配电网拓扑数据校验方法及系统
Yuan et al. Optimization of Expanding Gathering Pipeline Network in Gas Fields
CN114370984B (zh) 一种长距离供水管道水压试验渗漏位置的智能识别方法
KR102267558B1 (ko) 상수관로 실시간 가상수압 산정시스템 및 방법
Laucelli et al. Calibration of water distribution system using topological analysis
Cheng et al. Establishment and optimization of fluid pipe network models based on topological analysis algorithm
Liu et al. Optimal Design of District Metered Areas Based on Graph Theory and Multi-Objective Optimization:(175)
Kasyanenko et al. Selection of gas flow distribution in the gas gathering network

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant