CN113704936A - 考虑通道容量和转弯半径约束的电缆路径多目标优化方法 - Google Patents
考虑通道容量和转弯半径约束的电缆路径多目标优化方法 Download PDFInfo
- Publication number
- CN113704936A CN113704936A CN202110856624.3A CN202110856624A CN113704936A CN 113704936 A CN113704936 A CN 113704936A CN 202110856624 A CN202110856624 A CN 202110856624A CN 113704936 A CN113704936 A CN 113704936A
- Authority
- CN
- China
- Prior art keywords
- cable
- channel
- constraint
- optimization
- follows
- 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
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 95
- 238000000034 method Methods 0.000 title claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims description 24
- 238000005094 computer simulation Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000009413 insulation Methods 0.000 claims description 8
- 238000012067 mathematical method Methods 0.000 abstract description 2
- 238000010187 selection method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 244000050510 Cunninghamia lanceolata Species 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 241000764238 Isis Species 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- -1 j ∈ [1 Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/16—Cables, cable trees or wire harnesses
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Description
技术领域
本发明涉及电缆路径规划领域,更具体地说它是一种考虑通道容量和转弯半径约束的电缆路径多目标优化方法。
背景技术
无论对于大型水电站、变电站或者小型厂房,电缆路径规划一直是一项复杂而高难度的工作;当前,电缆路径规划主要依靠人工手动选择电缆路径,然后生成电缆敷设表;该做法不仅工作量巨大,耗时耗力,还对设计人员的经验和素质有较高的要求,尤其对于大型厂房而言,各类电缆总量可能高达上万根,设计人员需要对整个厂房的电缆走向和数量以及电缆通道的走向和容量有一个清晰全面的了解,往往需要花费数周甚至更长的时间才能完成整个厂房的电缆路径规划工作;若规划不合理,则很容易出现电缆堆积、通道容量无法满足要求等问题;此外,单纯依赖人工进行规划无法使电缆总长度或总成本达到最优,还可能出现电缆总量估算不准确,导致后期电缆成本增加或者浪费。
对于以上问题,目前相关研究主要集中在对单条电缆的最短路径搜索方案上,较少从全局出发对所有电缆的路径进行统一规划;最短路径搜索方法虽然可以找到所有电缆的最短路径,但可能无法满足电缆通道的容量约束、电缆转弯半径约束等,即可能出现大量电缆通向某一通道,超过该通道的最大容量,而其他部分通道无电缆通过的情况,或出现电缆转弯半径过小的情况,后续还需人工进行大量调整以满足实际敷设需求,电缆路径规划质量和效率均无法达到最优。
因此,亟需从全局角度对所有电缆进行统一规划,提出一种能满足通道容量和电缆转弯半径等约束且电缆总长度或总成本和换层数量最优的多电缆路径规划算法很有必要。
发明内容
本发明的目的是为了克服上述背景技术的不足之处,而提供一种基于图论的多电缆路径规划方法。
为了实现上述目的,本发明的技术方案为:考虑通道容量和转弯半径约束的电缆路径多目标优化方法,其特征在于,包括以下步骤:
步骤1:对所有的电缆通道进行编号,其中,不含分支的顺直通道视作一个通道段;
步骤3:以电缆总长度为第一优化目标,目标函数如下:
其中,fd表示以电缆总长度最小为优化目标,di表示通道i的长度,单位为米;
步骤4:以电缆换层数最小为第二优化目标,目标函数如下:
其中,fl表示以电缆换层数最小为优化目标,li表示每根通道所在层数,un表示计算矩阵内不同元素的数量;
步骤5:对电缆路径设置约束条件
约束1,起点通道和终点通道约束:根据电缆清册确定每条电缆的起点通道和终点通道,对于起点通道和终点通道而言,有且仅有一条通道与其相连,约束条件如下:
其中,sk表示电缆k的起点通道,ek表示电缆k的终点通道;Mi,j表示通道i与通道j之间的连通性矩阵,当通道i与通道j之间存在连接,则Mi,j=1时,否则Mi,j=0,当i=j时,设Mi,j=1;
约束2,连通性约束:对于第k条电缆经过的非起点和非终点的通道,必有两通道与其相连;即如果电缆通过某条非起点和非终点的通道,那么在电缆经过的通道中,必定有两条通道与该通道相连,约束条件如下:
约束3,单一性约束:对于每条电缆而言,每个通道只经过一次或者不经过,约束条件如下:
约束4,通道容量约束:每个通道所敷设的电缆数量有一定的限制,需添加通道容量约束,约束条件如下:
其中,rk表示电缆k的截面参数,考虑敷设空隙,该参数可为按矩形计算的电缆截面积,Si表示通道i所能铺设的最大电缆数量,可根据实际情况给出可摆放的最大截面积;
约束5,电缆转弯半径约束:对于某些电缆而言,为保障其绝缘性能,对其转弯半径有一定的要求;设定约束条件如下:
其中,Ai,j表示通道i与通道j之间的转弯半径,当Mi,j=0时,Ai,j=0;当i=j时,设表示电缆k的最小转弯半径要求;以上公式即表示,对于某一条电缆而言,若其经过通道i和通道j,且通道i和j之间有连接,则需要满足最小转弯半径约束;
步骤6:针对上述电缆路径多目标规划模型,利用求解器分别针对第一优化目标和第二优化目标进行求解,得到最优结果f1min和f2min,采用加权法将多目标优化转化为单目标f3优化,如下所示:
f3=w1f1/f1min+w2f2/f2min
w1+w2=1
0≤w1,w2≤1
按一定步长对w1进行赋值,循环调用求解器对以f3为目标函数的优化模型进行求解,即可得到所有电缆的Pareto解集。
考虑通道容量和转弯半径约束的电缆路径多目标优化方法,其特征在于,包括以下步骤:
步骤1:对所有的电缆通道进行编号,其中,不含分支的顺直通道视作一个通道段;
步骤3,以电缆总成本最小为第一优化目标,构建目标函数如下:
其中,fp表示以电缆总成本最小为优化目标,di表示通道i的长度,单位为米;pk表示第k条电缆的单位长度价格,单位元/米;
步骤4:以电缆换层数最小为第二优化目标,目标函数如下:
其中,fl表示以电缆换层数最小为优化目标,li表示每根通道所在层数,un表示计算矩阵内不同元素的数量;
步骤5:对电缆路径设置约束条件
约束1,起点通道和终点通道约束:根据电缆清册确定每条电缆的起点通道和终点通道,对于起点通道和终点通道而言,有且仅有一条通道与其相连,约束条件如下:
其中,sk表示电缆k的起点通道,ek表示电缆k的终点通道;Mi,j表示通道i与通道j之间的连通性矩阵,当通道i与通道j之间存在连接,则Mi,j=1时,否则Mi,j=0,当i=j时,设Mi,j=1;
约束2,连通性约束:对于第k条电缆经过的非起点和非终点的通道,必有两通道与其相连;即如果电缆通过某条非起点和非终点的通道,那么在电缆经过的通道中,必定有两条通道与该通道相连,约束条件如下:
约束3,单一性约束:对于每条电缆而言,每个通道只经过一次或者不经过,约束条件如下:
约束4,通道容量约束:每个通道所敷设的电缆数量有一定的限制,需添加通道容量约束,约束条件如下:
其中,rk表示电缆k的截面参数,,考虑敷设空隙,该参数可为按矩形计算的电缆截面积,Si表示通道i所能铺设的最大电缆数量,可根据实际情况给出可摆放的最大截面积;
约束5,电缆转弯半径约束:对于某些电缆而言,为保障其绝缘性能,对其转弯半径有一定的要求;设定约束条件如下:
其中,Ai,j表示通道i与通道j之间的转弯半径,当Mi,j=0时,Ai,j=0;当i=j时,设Ai,j=0;表示电缆k的最小转弯半径要求;以上公式即表示,对于某一条电缆而言,若其经过通道i和通道j,且通道i和j之间有连接,则需要满足最小转弯半径约束;
步骤6:针对上述电缆路径多目标规划模型,利用求解器分别针对第一优化目标和第二优化目标进行求解,得到最优结果f1min和f2min,采用加权法将多目标优化转化为单目标f3优化,如下所示:
f3=w1f1/f1min+w2f2/f2min
w1+w2=1
0≤w1,w2≤1
按一定步长对w1进行赋值,循环调用求解器对以f3为目标函数的优化模型进行求解,即可得到所有电缆的Pareto解集。
在上述技术方案中,在步骤2中,针对敷设的不同类型的电缆路径规划,需将不同类型的电缆分别建模,分别设置不同的决策变量。
在上述技术方案中,步骤4中,为后续仿真建模需要,需对fl进行线性化,目标函数可转化为以下形式:
在上述技术方案中,步骤5的约束2中,为后续仿真建模需要,需对约束2进行线性化,约束条件可转化为以下形式:
在上述技术方案中,步骤5的约束5中,为后续仿真建模需要,需对约束5进行线性化,约束条件可转化为以下形式:
与现有技术相比,本发明具有以下优点:
1)本发明的电缆路径规划方法适用于水电站、变电站以及其他厂房的动力站等不同场合的电缆路径规划,适用各种类型、规格的电缆,同时适用桥架、电缆沟、穿管等各种形式的电缆通道。
2)本发明将电缆路径规划问题转换为混合整数线性组合优化问题,采用数学方法替代传统的依赖人工的电缆路径规划方法,提高了电缆路径规划的效率和质量,减轻了设计人员的工作量。
3)本发明采用的优化算法,可在满足通道容量、电缆转弯半径等约束条件下,获得以电缆总长度或总成本和换层数最小为优化目标的电缆路径最优解集,相较于当前的解决方法,具有规划效率和质量高的优点,相较于仅以电缆总长度或总成本最小或仅以换层数最小为目标的单目标优化方案,本发明可为设计人员提供更多可选的路径规划方案。
附图说明
图1为本发明实施例的应用场景示意图;图中1、2、3、4、5、6、7、8、9、10分别表示不同的通道。
图2为本发明实施例的所有电缆的Pareto解集;图中横坐标表示电缆总长度或总价格,纵坐标表示电缆总换层数。
具体实施方式
下面结合附图详细说明本发明的实施情况,但它们并不构成对本发明的限定,仅作举例而已。同时通过说明使本发明的优点变得更加清楚和容易理解。
参阅附图可知:考虑通道容量和转弯半径约束的电缆路径多目标优化方法,包括以下步骤:
步骤1:对所有的电缆通道进行编号,其中,不含分支的顺直通道视作一个通道段;
针对敷设在不同通道的不同类型的电缆路径规划,如高压电缆与低压电缆的路径规划,需看作两个问题分别建模,分别设置不同的决策变量;
步骤3:以电缆总长度为第一优化目标,目标函数如下:
其中,fd表示以电缆总长度最小为优化目标,di表示通道i的长度,单位为米;
步骤4:以电缆换层数最小为第二优化目标,目标函数如下:
其中,fl表示以电缆换层数最小为优化目标,li表示每根通道所在层数,un表示计算矩阵内不同元素的数量;
步骤5:对电缆路径设置约束条件
约束1,起点通道和终点通道约束:根据电缆清册确定每条电缆的起点通道和终点通道,对于起点通道和终点通道而言,有且仅有一条通道与其相连,约束条件如下:
其中,sk表示电缆k的起点通道,ek表示电缆k的终点通道;Mi,j表示通道i与通道j之间的连通性矩阵,当通道i与通道j之间存在连接,则Mi,j=1时,否则Mi,j=0,当i=j时,设Mi,j=1;
约束2,连通性约束:对于第k条电缆经过的非起点和非终点的通道,必有两通道与其相连;即如果电缆通过某条非起点和非终点的通道,那么在电缆经过的通道中,必定有两条通道与该通道相连,约束条件如下:
约束3,单一性约束:对于每条电缆而言,每个通道只经过一次或者不经过,约束条件如下:
约束4,通道容量约束:每个通道所敷设的电缆数量有一定的限制,需添加通道容量约束,约束条件如下:
其中,rk表示电缆k的截面参数,考虑敷设空隙,该参数可为按矩形计算的电缆截面积;Si表示通道i所能铺设的最大电缆数量(根据实际情况给出可摆放的最大截面积);
约束5,电缆转弯半径约束:对于某些电缆而言,为保障其绝缘性能,对其转弯半径有一定的要求;设定约束条件如下:
其中,Ai,j表示通道i与通道j之间的转弯半径,当Mi,j=0时,Ai,j=0;当i=j时,设Ai,j=0;表示电缆k的最小转弯半径要求;以上公式即表示,对于某一条电缆而言,若其经过通道i和通道j,且通道i和j之间有连接,则需要满足最小转弯半径约束;
步骤6:针对上述电缆路径多目标规划模型,利用求解器分别针对第一优化目标和第二优化目标进行求解,得到最优结果f1min和f2min,采用加权法将多目标优化转化为单目标f3优化,如下所示:
f3=w1f1/f1min+w2f2/f2min
w1+w2=1
0≤w1,w2≤1
按一定步长对w1进行赋值,循环调用求解器对以f3为目标函数的优化模型进行求解,即可得到所有电缆的Pareto解集(可选路径集)。
考虑通道容量和转弯半径约束的电缆路径多目标优化方法,包括以下步骤:
步骤1:对所有的电缆通道进行编号,其中,不含分支的顺直通道视作一个通道段;
针对敷设在不同通道的不同类型的电缆路径规划,如高压电缆与低压电缆的路径规划,需看作两个问题分别建模,分别设置不同的决策变量;
步骤3,以电缆总成本最小为第一优化目标,构建目标函数如下:
其中,fp表示以电缆总成本最小为优化目标,di表示通道i的长度,单位为米;pk表示第k条电缆的单位长度价格,单位元/米。
步骤4:以电缆换层数最小为第二优化目标,目标函数如下:
其中,fl表示以电缆换层数最小为优化目标,li表示每根通道所在层数,un表示计算矩阵内不同元素的数量;
步骤5:对电缆路径设置约束条件
约束1,起点通道和终点通道约束:根据电缆清册确定每条电缆的起点通道和终点通道,对于起点通道和终点通道而言,有且仅有一条通道与其相连,约束条件如下:
其中,sk表示电缆k的起点通道,ek表示电缆k的终点通道;Mi,j表示通道i与通道j之间的连通性矩阵,当通道i与通道j之间存在连接,则Mi,j=1时,否则Mi,j=0,当i=j时,设Mi,j=1;
约束2,连通性约束:对于第k条电缆经过的非起点和非终点的通道,必有两通道与其相连;即如果电缆通过某条非起点和非终点的通道,那么在电缆经过的通道中,必定有两条通道与该通道相连,约束条件如下:
约束3,单一性约束:对于每条电缆而言,每个通道只经过一次或者不经过,约束条件如下:
约束4,通道容量约束:每个通道所敷设的电缆数量有一定的限制,需添加通道容量约束,约束条件如下:
其中,rk表示电缆k的截面参数,考虑敷设空隙,该参数可为按矩形计算的电缆截面积;Si表示通道i所能铺设的最大电缆数量(根据实际情况给出可摆放的最大截面积);
约束5,电缆转弯半径约束:对于某些电缆而言,为保障其绝缘性能,对其转弯半径有一定的要求;设定约束条件如下:
其中,Ai,j表示通道i与通道j之间的转弯半径,当Mi,j=0时,Ai,j=0;当i=j时,设Ai,j=0;表示电缆k的最小转弯半径要求;以上公式即表示,对于某一条电缆而言,若其经过通道i和通道j,且通道i和j之间有连接,则需要满足最小转弯半径约束;
步骤6:针对上述电缆路径多目标规划模型,利用求解器分别针对第一优化目标和第二优化目标进行求解,得到最优结果f1min和f2min,采用加权法将多目标优化转化为单目标f3优化,如下所示:
f3=w1f1/f1min+w2f2/f2min
w1+w2=1
0≤w1,w2≤1
按一定步长对w1进行赋值,循环调用求解器对以f3为目标函数的优化模型进行求解,即可得到所有电缆的Pareto解集(可选路径集)。
步骤2中,针对不同类型的电缆路径规划,需将不同类型的电缆分别建模,分别设置不同的决策变量。
步骤4中,为后续仿真建模需要,需对fl进行线性化,目标函数可转化为以下形式:
步骤5的约束2中,为后续仿真建模需要,需对约束2进行线性化,约束条件可转化为以下形式:
步骤5的约束5中,为后续仿真建模需要,需对约束5进行线性化,约束条件可转化为以下形式:
实施例1
考虑通道容量和转弯半径约束的电缆路径多目标优化规划方法,具体步骤如下:
步骤1:对所有的电缆通道进行编号,其中,不含分支的顺直通道视作一个通道段;如图1所示,图中数字所对应的黑色连线表示1-10号电缆通道,其中,1号电缆通道长度为2米,2号电缆通道长度为3米,3号电缆通道长度为4米,4号电缆通道长度为5米,5号电缆通道长度为10米,6号电缆通道长度为4米,7号电缆通道长度为7米,8号电缆通道长度为3米,9号电缆通道长度为1米,10号电缆通道长度为1米;
步骤2:设置决策变量为为布尔型变量,表示第k条电缆所经过的通道标志位,当时,表示第k条电缆经过通道i,否则不经过;如图1所示,共有10条通道,则对每条电缆而言,为1×10的矩阵,i∈[1,10];本实施例仅考虑可以敷设在同一通道的电缆,针对敷设在不同通道的不同类型的电缆路径规划,如高压电缆与低压电缆的路径规划,需看作两个问题分别建模,分别设置不同的决策变量;
步骤3:本实施例以电缆总长度最小为第一优化目标,目标函数如下:
其中,fd表示以电缆总长度最小为优化目标,di表示通道i的长度,单位为米;
根据附图1,可得di如下所示:
步骤4:以电缆换层数最小为第二优化目标,目标函数如下:
其中,fl表示以电缆换层数最小为优化目标,li表示每段通道所在层数,un表示计算矩阵内不同元素的数量;本实施例设各段通道的层数如下所示:
为后续仿真建模需要,需对fl进行线性化,上述目标函数可转化为以下形式:
步骤5:对电缆路径设置约束条件
约束1,起点通道和终点通道约束:根据电缆清册确定每条电缆的起点通道和终点通道,对于起点通道和终点通道而言,有且仅有一条通道与其相连,约束条件如下:
其中,sk和ek分别表示电缆k的起点通道和终点通道;Mi,j表示通道i与通道j之间的连通性,当通道i与通道j之间存在连接,则Mi,j=1时,否则Mi,j=0;当i=j时,设Mi,j=1;
根据附图1,可得通道连通性矩阵如下所示:
M=[[1,1,0,0,1,1,1,0,1,0],[1,1,1,0,0,1,1,1,0,0],[0,1,1,1,0,1,0,1,0,1],[0,0,1,1,1,0,1,1,0,1],[1,0,0,1,1,1,1,1,1,0],[1,1,1,0,1,1,0,1,1,0],[1,1,0,1,1,0,1,1,0,0],[0,1,1,1,1,1,1,1,0,0],[1,0,0,0,1,1,0,0,1,0],[0,0,1,1,0,0,0,0,0,1]];
本实施例中设置电缆清册如下表所示:
电缆编号k | 起点通道s<sub>k</sub> | 终点通道e<sub>k</sub> |
1 | 9 | 10 |
2 | 1 | 3 |
3 | 9 | 10 |
约束2,连通性约束:对于非起点和非终点的通道,必有两通道与其相连;即如果电缆通过某条非起点和非终点的通道,那么在电缆经过的通道中,必定有两条通道与该通道相连,约束条件如下:
为后续仿真建模需要,需对上述约束进行线性化,上述约束条件可转化为以下形式:
约束3,单一性约束:对于每条电缆而言,每个通道只经过一次或者不经过,约束条件如下:
约束4,通道容量约束:每个通道所敷设的电缆数量有一定的限制,需添加通道容量约束,约束条件如下:
其中,rk表示电缆k的截面参数,考虑敷设空隙,该参数可为按矩形计算的电缆截面积,Si表示通道i所能铺设的最大电缆面积(根据实际情况给出可摆放的最大截面积);
本实施例中设所有电缆规格一致,设rk=1,每段电缆通道的容量为3,则Si如下:
约束5,电缆转弯半径约束:对于某些电缆而言,为保障其绝缘性能,对其转弯半径有一定的要求。设定约束条件如下:
为后续仿真建模需要,需对上述约束进行线性化,上述约束条件可转化为以下形式:
Amin=[0.5 0.5 0.5]
A=[[1,1,0,0,1,1,1,0,1,0],[1,0,1,0,0,1,1,1,0,0],[0,1,0,1,0,1,0,1,0,1],[0,0,1,0,1,0,1,1,0,1],[1,0,0,1,0,1,1,1,1,0],[1,1,1,0,1,0,0,1,1,0],[1,1,0,1,1,0,0,1,0,0],[0,1,1,1,1,1,1,0,0,0],[1,0,0,0,1,1,0,0,0,0],[0,0,1,1,0,0,0,0,0,0]];
步骤5:针对上述电缆路径多目标规划模型,利用现有求解器(如中国杉树科技开发的基数求解器,COPT)分别针对第一优化目标和第二优化目标进行求解,得到最优结果f1min和f2min,分别为25米和4层,采用加权法将多目标优化转化为单目标f3优化,如下所示:
f3=w1f1/f1min+w2f2/f2min
w1+w2=1
0≤w1,w2≤1
按0.05的步长对w1进行赋值,循环调用求解器对以f3为目标函数的优化模型进行求解,即可得到所有电缆的Pareto解集(可选路径集),如图2所示;
由图2可知,对电缆路径多目标规划模型求解后可得到3个可选路径解;
若追求电缆总路径最短,则可选择第二个路径解,即电缆总长度为29米,电缆的总层数为6层;此时,各条电缆的路径如下:
第1条电缆的最优路径为[9,6,3,10]
第2条电缆的最优路径为[1,2,3]
第3条电缆的最优路径为[9,6,3,10]
若追求电缆总层数最少,则可选择第三个路径解,即电缆总长度为43米,电缆的总层数为4层。此时,各条电缆的路径如下:
第1条电缆的最优路径为[9,5,4,10]
第2条电缆的最优路径为[1,2,3]
第3条电缆的最优路径为[9,5,4,10]
根据不同的需求有不同的电缆路径选择方法,具体选取方法还需根据现场情况和工程技术人员的具体需求确定。
实施例2
考虑通道容量和转弯半径约束的电缆路径多目标优化规划方法,具体步骤如下:
步骤1:对所有的电缆通道进行编号,其中,不含分支的顺直通道视作一个通道段;如图1所示,图中数字所对应的黑色连线表示1-10号电缆通道,其中,1号电缆通道长度为2米,2号电缆通道长度为3米,3号电缆通道长度为4米,4号电缆通道长度为5米,5号电缆通道长度为10米,6号电缆通道长度为4米,7号电缆通道长度为7米,8号电缆通道长度为3米,9号电缆通道长度为1米,10号电缆通道长度为1米;
步骤2:设置决策变量为为布尔型变量,表示第k条电缆所经过的通道标志位,当时,表示第k条电缆经过通道i,否则不经过;如图1所示,共有10条通道,则对每条电缆而言,为1×10的矩阵,i∈[1,10];本实施例仅考虑可以敷设在同一通道的电缆,针对敷设在不同通道的不同类型的电缆路径规划,如高压电缆与低压电缆的路径规划,需看作两个问题分别建模,分别设置不同的决策变量;
步骤3,本实施例以电缆总成本最小为第一优化目标,构建目标函数如下:
其中,fp表示以电缆总成本最小为优化目标,di表示通道i的长度,单位为米;pk表示第k条电缆的单位长度价格,单位元/米;
根据附图1,可得di如下所示:
步骤4:以电缆换层数最小为第二优化目标,目标函数如下:
其中,fl表示以电缆换层数最小为优化目标,li表示每段通道所在层数,un表示计算矩阵内不同元素的数量;本实施例设各段通道的层数如下所示:
为后续仿真建模需要,需对fl进行线性化,上述目标函数可转化为以下形式:
步骤5:对电缆路径设置约束条件
约束1,起点通道和终点通道约束:根据电缆清册确定每条电缆的起点通道和终点通道,对于起点通道和终点通道而言,有且仅有一条通道与其相连,约束条件如下:
其中,sk和ek分别表示电缆k的起点通道和终点通道;Mi,j表示通道i与通道j之间的连通性,当通道i与通道j之间存在连接,则Mi,j=1时,否则Mi,j=0;当i=j时,设Mi,j=1;
根据附图1,可得通道连通性矩阵如下所示:
M=[[1,1,0,0,1,1,1,0,1,0],[1,1,1,0,0,1,1,1,0,0],[0,1,1,1,0,1,0,1,0,1],[0,0,1,1,1,0,1,1,0,1],[1,0,0,1,1,1,1,1,1,0],[1,1,1,0,1,1,0,1,1,0],[1,1,0,1,1,0,1,1,0,0],[0,1,1,1,1,1,1,1,0,0],[1,0,0,0,1,1,0,0,1,0],[0,0,1,1,0,0,0,0,0,1]];
本实施例中设置电缆清册如下表所示:
电缆编号k | 起点通道s<sub>k</sub> | 终点通道e<sub>k</sub> | 电缆价格 |
1 | 9 | 10 | 1 |
2 | 1 | 3 | 1 |
3 | 9 | 10 | 1 |
可得电缆价格矩阵pk=[1,1,1]。
约束2,连通性约束:对于非起点和非终点的通道,必有两通道与其相连;即如果电缆通过某条非起点和非终点的通道,那么在电缆经过的通道中,必定有两条通道与该通道相连,约束条件如下:
为后续仿真建模需要,需对上述约束进行线性化,上述约束条件可转化为以下形式:
约束3,单一性约束:对于每条电缆而言,每个通道只经过一次或者不经过,约束条件如下:
约束4,通道容量约束:每个通道所敷设的电缆数量有一定的限制,需添加通道容量约束,约束条件如下:
其中,rk表示电缆k的截面参数,考虑敷设空隙,该参数可为按矩形计算的电缆截面积,Si表示通道i所能铺设的最大电缆面积(根据实际情况给出可摆放的最大截面积);
本实施例中设所有电缆规格一致,设rk=1,每段电缆通道的容量为3,则Si如下:
约束5,电缆转弯半径约束:对于某些电缆而言,为保障其绝缘性能,对其转弯半径有一定的要求。设定约束条件如下:
为后续仿真建模需要,需对上述约束进行线性化,上述约束条件可转化为以下形式:
Amin=[0.5 0.5 0.5]
A=[[1,1,0,0,1,1,1,0,1,0],[1,0,1,0,0,1,1,1,0,0],[0,1,0,1,0,1,0,1,0,1],[0,0,1,0,1,0,1,1,0,1],[1,0,0,1,0,1,1,1,1,0],[1,1,1,0,1,0,0,1,1,0],[1,1,0,1,1,0,0,1,0,0],[0,1,1,1,1,1,1,0,0,0],[1,0,0,0,1,1,0,0,0,0],[0,0,1,1,0,0,0,0,0,0]];步骤5:针对上述电缆路径多目标规划模型,利用现有求解器(如中国杉树科技开发的基数求解器,COPT)分别针对第一优化目标和第二优化目标进行求解,得到最优结果f1min和f2min,分别为25元和4层,采用加权法将多目标优化转化为单目标f3优化,如下所示:
f3=w1f1/f1min+w2f2/f2min
w1+w2=1
0≤w1,w2≤1
按0.05的步长对w1进行赋值,循环调用求解器对以f3为目标函数的优化模型进行求解,即可得到所有电缆的Pareto解集(可选路径集),如图2所示;
由图2可知,对电缆路径多目标规划模型求解后可得到3个可选路径解;
若追求电缆总成本最小,则可选择第二个路径解,即电缆总成本为29米,电缆的总层数为6层;此时,各条电缆的路径如下:
第1条电缆的最优路径为[9,6,3,10]
第2条电缆的最优路径为[1,2,3]
第3条电缆的最优路径为[9,6,3,10]
若追求电缆总层数最少,则可选择第三个路径解,即电缆总成本为43元,电缆的总层数为4层。此时,各条电缆的路径如下:
第1条电缆的最优路径为[9,5,4,10]
第2条电缆的最优路径为[1,2,3]
第3条电缆的最优路径为[9,5,4,10]
根据不同的需求有不同的电缆路径选择方法,具体选取方法还需根据现场情况和工程技术人员的具体需求确定。
其它未说明的部分均属于现有技术。
Claims (6)
1.考虑通道容量和转弯半径约束的电缆路径多目标优化方法,其特征在于,包括以下步骤:
步骤1:对所有的电缆通道进行编号,其中,不含分支的顺直通道视作一个通道段;
其中,fd表示以电缆总长度最小为第一优化目标,di表示通道i的长度,单位为米;
步骤4:以电缆换层数最小为第二优化目标,目标函数如下:
其中,fl表示以电缆换层数最小为第二优化目标,li表示每根通道所在层数,un表示计算矩阵内不同元素的数量;步骤5:对电缆路径设置约束条件
约束1,起点通道和终点通道约束:根据电缆清册确定每条电缆的起点通道和终点通道,对于起点通道和终点通道而言,有且仅有一条通道与其相连,约束条件如下:
其中,sk表示电缆k的起点通道,ek表示电缆k的终点通道;Mi,j表示通道i与通道j之间的连通性矩阵,当通道i与通道j之间存在连接,则Mi,j=1时,否则Mi,j=0,当i=j时,设Mi,j=1;
约束2,连通性约束:对于第k条电缆经过的非起点和非终点的通道,必有两通道与其相连;即如果电缆通过某条非起点和非终点的通道,那么在电缆经过的通道中,必定有两条通道与该通道相连,约束条件如下:
约束3,单一性约束:对于每条电缆而言,每个通道只经过一次或者不经过,约束条件如下:
约束4,通道容量约束:每个通道所敷设的电缆数量有一定的限制,需添加通道容量约束,约束条件如下:
其中,rk表示电缆k的截面参数,考虑敷设空隙,该参数可为按矩形计算的电缆截面积,Si表示通道i所能铺设的最大电缆数量,可根据实际情况给出可摆放的最大截面积;
约束5,电缆转弯半径约束:对于某些电缆而言,为保障其绝缘性能,对其转弯半径有一定的要求;设定约束条件如下:
其中,Ai,j表示通道i与通道j之间的转弯半径,当Mi,j=0时,Ai,j=0;当i=j时,设Ai,j=0;表示电缆k的最小转弯半径要求;以上公式即表示,对于某一条电缆而言,若其经过通道i和通道j,且通道i和j之间有连接,则需要满足最小转弯半径约束;
步骤6:针对上述电缆路径多目标规划模型,利用求解器分别针对第一优化目标和第二优化目标进行求解,得到最优结果f1min和f2min,采用加权法将多目标优化转化为单目标f3优化,如下所示:
f3=w1f1/f1min+w2f2/f2min
w1+w2=1
0≤w1,w2≤1
按一定步长对w1进行赋值,循环调用求解器对以f3为目标函数的优化模型进行求解,即可得到所有电缆的Pareto解集。
2.考虑通道容量和转弯半径约束的电缆路径多目标优化方法,其特征在于,包括以下步骤:
步骤1:对所有的电缆通道进行编号,其中,不含分支的顺直通道视作一个通道段;
步骤3,以电缆总成本最小为第一优化目标,构建目标函数如下:
其中,fp表示以电缆总成本最小为优化目标,di表示通道i的长度,单位为米;pk表示第k条电缆的单位长度价格,单位元/米;
步骤4:以电缆换层数最小为第二优化目标,目标函数如下:
其中,fl表示以电缆换层数最小为优化目标,li表示每根通道所在层数,un表示计算矩阵内不同元素的数量;
步骤5:对电缆路径设置约束条件
约束1,起点通道和终点通道约束:根据电缆清册确定每条电缆的起点通道和终点通道,对于起点通道和终点通道而言,有且仅有一条通道与其相连,约束条件如下:
其中,sk表示电缆k的起点通道,ek表示电缆k的终点通道;Mi,j表示通道i与通道j之间的连通性矩阵,当通道i与通道j之间存在连接,则Mi,j=1时,否则Mi,j=0,当i=j时,设Mi,j=1;
约束2,连通性约束:对于第k条电缆经过的非起点和非终点的通道,必有两通道与其相连;即如果电缆通过某条非起点和非终点的通道,那么在电缆经过的通道中,必定有两条通道与该通道相连,约束条件如下:
约束3,单一性约束:对于每条电缆而言,每个通道只经过一次或者不经过,约束条件如下:
约束4,通道容量约束:每个通道所敷设的电缆数量有一定的限制,需添加通道容量约束,约束条件如下:
其中,rk表示电缆k的截面参数,考虑敷设空隙,该参数可为按矩形计算的电缆截面积,Si表示通道i所能铺设的最大电缆数量,可根据实际情况给出可摆放的最大截面积;
约束5,电缆转弯半径约束:对于某些电缆而言,为保障其绝缘性能,对其转弯半径有一定的要求;设定约束条件如下:
其中,Ai,j表示通道i与通道j之间的转弯半径,当Mi,j=0时,Ai,j=0;当i=j时,设Ai,j=0;表示电缆k的最小转弯半径要求;以上公式即表示,对于某一条电缆而言,若其经过通道i和通道j,且通道i和j之间有连接,则需要满足最小转弯半径约束;
步骤6:针对上述电缆路径多目标规划模型,利用求解器分别针对第一优化目标和第二优化目标进行求解,得到最优结果f1min和f2min,采用加权法将多目标优化转化为单目标f3优化,如下所示:
f3=w1f1/f1min+w2f2/f2min
w1+w2=1
0≤w1,w2≤1
按一定步长对w1进行赋值,循环调用求解器对以f3为目标函数的优化模型进行求解,即可得到所有电缆的Pareto解集。
3.根据权利要求1或2所述的考虑通道容量和转弯半径约束的电缆路径多目标优化方法,其特征在于,在步骤2中,针对不同类型的电缆路径规划,需将不同类型的电缆分别建模,分别设置不同的决策变量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110856624.3A CN113704936B (zh) | 2021-07-28 | 2021-07-28 | 考虑通道容量和转弯半径约束的电缆路径多目标优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110856624.3A CN113704936B (zh) | 2021-07-28 | 2021-07-28 | 考虑通道容量和转弯半径约束的电缆路径多目标优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113704936A true CN113704936A (zh) | 2021-11-26 |
CN113704936B CN113704936B (zh) | 2023-06-27 |
Family
ID=78650715
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110856624.3A Active CN113704936B (zh) | 2021-07-28 | 2021-07-28 | 考虑通道容量和转弯半径约束的电缆路径多目标优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113704936B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116708447A (zh) * | 2023-08-08 | 2023-09-05 | 天津市奥讯通电缆科技发展有限公司 | 一种综合分叉电缆的数据耦合方法 |
CN118134217A (zh) * | 2024-05-08 | 2024-06-04 | 成都建工雅安建设有限责任公司 | 一种基于bim技术的电缆桥架施工管理方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106096079A (zh) * | 2016-05-26 | 2016-11-09 | 中国能源建设集团广东省电力设计研究院有限公司 | 提高电缆敷设长度精度的设计方法及系统 |
CN107766667A (zh) * | 2017-10-27 | 2018-03-06 | 中国空间技术研究院 | 一种空间运动电缆走向设计方法 |
CN111400899A (zh) * | 2020-03-13 | 2020-07-10 | 广州轨道交通建设监理有限公司 | 一种用于电缆敷设建模的处理方法、系统和存储介质 |
WO2021102623A1 (zh) * | 2019-11-25 | 2021-06-03 | 上海电气风电集团股份有限公司 | 风电场的电缆路径的规划方法、系统、介质及电子设备 |
US11036899B1 (en) * | 2020-02-11 | 2021-06-15 | CyberSecure IPS, LLC | Method and system for determining optimal pathways for installing cables in an infrastructure |
-
2021
- 2021-07-28 CN CN202110856624.3A patent/CN113704936B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106096079A (zh) * | 2016-05-26 | 2016-11-09 | 中国能源建设集团广东省电力设计研究院有限公司 | 提高电缆敷设长度精度的设计方法及系统 |
CN107766667A (zh) * | 2017-10-27 | 2018-03-06 | 中国空间技术研究院 | 一种空间运动电缆走向设计方法 |
WO2021102623A1 (zh) * | 2019-11-25 | 2021-06-03 | 上海电气风电集团股份有限公司 | 风电场的电缆路径的规划方法、系统、介质及电子设备 |
US11036899B1 (en) * | 2020-02-11 | 2021-06-15 | CyberSecure IPS, LLC | Method and system for determining optimal pathways for installing cables in an infrastructure |
CN111400899A (zh) * | 2020-03-13 | 2020-07-10 | 广州轨道交通建设监理有限公司 | 一种用于电缆敷设建模的处理方法、系统和存储介质 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116708447A (zh) * | 2023-08-08 | 2023-09-05 | 天津市奥讯通电缆科技发展有限公司 | 一种综合分叉电缆的数据耦合方法 |
CN116708447B (zh) * | 2023-08-08 | 2023-10-20 | 天津市奥讯通电缆科技发展有限公司 | 一种综合分叉电缆的数据耦合方法 |
CN118134217A (zh) * | 2024-05-08 | 2024-06-04 | 成都建工雅安建设有限责任公司 | 一种基于bim技术的电缆桥架施工管理方法及系统 |
CN118134217B (zh) * | 2024-05-08 | 2024-07-12 | 成都建工雅安建设有限责任公司 | 一种基于bim技术的电缆桥架施工管理方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN113704936B (zh) | 2023-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113704936A (zh) | 考虑通道容量和转弯半径约束的电缆路径多目标优化方法 | |
CN104268682B (zh) | 一种有源配电网的规划方法及装置 | |
CN106503839B (zh) | 一种海上风电场环形集电网络分层规划方法 | |
CN111900732B (zh) | 考虑风电不确定性及无功潮流与电压约束的风电集群接入方式规划方法 | |
CN106655227B (zh) | 一种基于智能软开关的有源配电网馈线负载平衡方法 | |
CN108023364B (zh) | 基于凸差规划的配电网分布式电源最大接入能力计算方法 | |
CN102545258A (zh) | 一种大规模风电并网的电网优化规划方法 | |
CN109586303A (zh) | 一种配电网区域电压分布式协调控制方法及系统 | |
CN102522752A (zh) | 面向供电能力提高的配电网联络结构优化方法 | |
CN109728603A (zh) | 有源配电网分布式电源分区就地电压控制策略整定方法 | |
CN109829560B (zh) | 一种配电网可再生能源发电集群接入规划方法 | |
CN111082466B (zh) | 考虑风电不确定性的新能源接入与网架扩建优化方法 | |
CN110851932B (zh) | 一种低压配电网三相单线图初始布局自动生成方法 | |
CN103824124B (zh) | 一种针对电网公司的节能潜力评估方法 | |
CN113343426A (zh) | 考虑现实需求导向的城市配电网线路规方法及系统 | |
US8738332B2 (en) | System and process for client driven automated circuiting and branch circuit wiring | |
Hardy et al. | A greedy algorithm for optimizing offshore wind transmission topologies | |
CN113704934B (zh) | 基于图论的多电缆路径规划方法 | |
Bresesti et al. | An OPF-based procedure for fast TTC analyses | |
CN112418732B (zh) | 一种含多能量枢纽综合能源系统规划方法及系统 | |
CN111105089B (zh) | 一种考虑主配网协调的城市电网规划方法 | |
Hennig et al. | Ancillary services analysis of an offshore wind farm cluster–technical integration steps of a simulation tool | |
Çakmak et al. | An optimization-based approach for automated generation of residential low-voltage grid models using open data and open source software | |
Vrana et al. | Improved investment cost model and overhead cost consideration for high voltage direct current infrastructure | |
CN113704935A (zh) | 计及三维建模随机性的电缆通道属性及关联关系提取方法 |
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 |