CN111027184A - Ship power grid fault reconstruction convex optimization model considering reliability constraint - Google Patents

Ship power grid fault reconstruction convex optimization model considering reliability constraint Download PDF

Info

Publication number
CN111027184A
CN111027184A CN201911148323.4A CN201911148323A CN111027184A CN 111027184 A CN111027184 A CN 111027184A CN 201911148323 A CN201911148323 A CN 201911148323A CN 111027184 A CN111027184 A CN 111027184A
Authority
CN
China
Prior art keywords
load
ship
constraint
line
power grid
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
CN201911148323.4A
Other languages
Chinese (zh)
Other versions
CN111027184B (en
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.)
719th Research Institute of CSIC
Original Assignee
719th Research Institute of CSIC
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 719th Research Institute of CSIC filed Critical 719th Research Institute of CSIC
Priority to CN201911148323.4A priority Critical patent/CN111027184B/en
Publication of CN111027184A publication Critical patent/CN111027184A/en
Application granted granted Critical
Publication of CN111027184B publication Critical patent/CN111027184B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] 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
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

A ship power grid fault reconstruction convex optimization model considering reliability constraint comprises the following steps: establishing a reliability index calculation model based on the line-distribution board incidence matrix; aiming at reducing the load loss and the network loss of the ship power grid, considering the voltage constraint of the ship power grid, the reliability constraint of a system and an important load, and establishing a ship power grid fault reconstruction model; and converting the original model into a ship power grid fault reconstruction convex optimization model through second-order cone relaxation so as to solve and obtain a reconstruction strategy with the maximum load recovery quantity. The invention can minimize the power loss after the important load is recovered, ensure the reliable power supply of the system and the important load and improve the voltage quality. The method has high precision and can quickly obtain the global optimal solution.

Description

Ship power grid fault reconstruction convex optimization model considering reliability constraint
Technical Field
The invention relates to the technical field of power system analysis, in particular to a ship power grid fault reconstruction convex optimization model considering reliability constraint.
Background
When a ship executes a navigation task at sea, the electric power system of the ship is completely separated from the land electric power system, and the whole ship electric equipment can be powered only by the power supply system of the ship. In actual operation, due to problems of navigation damage, improper operation or equipment, various faults or abnormal operation states may occur in a ship power system, which affects safe and reliable operation of the whole network, and may even cause equipment damage or power interruption of the whole power grid, which affects fighting performance and navigation safety of the ship. The fault reconstruction strategy is a common power restoration method, but is different from a conventional power system, if a ship power system has a fault, the reliability of a reconstructed power grid needs to be considered when the power grid is reconstructed and the load is restored. In addition, the ship is supplied with power by a power station, so that the load capacity generally has no too much redundancy, and if major faults occur, the problems of line transmission capacity and voltage safety constraint are also a problem which cannot be ignored.
In conclusion, the method has very important significance in rapidly and maximally recovering the power supply of important loads and ensuring the safety of the ship operation mode aiming at the fault state of the ship power grid. The existing method for reconstructing the fault of the ship power system does not consider the reliability and voltage safety of the power grid after recovery, and the existing fault reconstruction model of the ship power system is usually a non-convex model and is usually solved by an intelligent algorithm based on population evolution, but has the problems of low convergence speed, easy falling into local optimization, low search efficiency and the like.
Therefore, it is necessary to provide a ship power grid reconstruction model capable of rapidly obtaining a global optimal solution and considering reliability and voltage safety operation constraints.
Disclosure of Invention
In order to solve the technical problems, the technical scheme adopted by the invention is to provide a ship power grid fault reconstruction convex optimization model considering reliability constraint, and the model is characterized by comprising the following steps of:
dividing the load of the ship power grid into different important levels based on a reliability index concept and a ship power grid fault reconstruction requirement, and defining the important degrees of the different load levels;
selecting I-grade loads according to the load classification result obtained in the step 1, and establishing a reliability index calculation model through a line-distribution board incidence matrix;
aiming at reducing the load loss and the network loss of the ship power grid, considering the voltage constraint of the ship power grid, the reliability constraint of a system and an important load, and establishing a ship power grid fault reconstruction model;
converting the ship power grid fault reconstruction model into a convex optimization model through second-order cone relaxation;
and solving to obtain the action switch which meets the constraint and minimizes the loss load quantity.
In the above method, the importance levels of the ship power grid load include:
the ship power grid is divided into a class I load, a class II load and a class III load according to the importance of the load.
I-level load: this load controls the vitality of the ship and the safety of the passengers belonging to the un-unloadable load; in any case a first order load must be provided. The primary loads of the ship system comprise a boiler, a generator, a weapon system, an electronic countermeasure system, a medical operating room and the like, and at least two power supply paths are required for connecting important loads. If a certain load is certified as a primary load (e.g. a ship propulsion system) in any one of the ship's operational tasks, it has to be connected to the ship's electrical system by Automatic Bus Transfer (ABT). Automatic bus transfer is a device that detects power loss from a normal power supply. When the normal power supply fails, the ABT can automatically and quickly disconnect the load from the normal power supply and connect the load power supply with the standby power supply. Non-critical loads have only one power supply path.
And II-stage load: this load is important to the proper operation of the vessel, but allows it to be offloaded or transferred to other platforms for power if necessary, in case of substantial loss of load. These loads include: cargo lifts, sea water pumps, some radar loads, etc.
Stage III load: the ship can be unloaded immediately when necessary, and the service life and the safety of the ship cannot be adversely affected; the method mainly comprises the following steps: residential heating systems, kitchen electrical loads, and the like.
In the above method, the ship grid reliability constraint includes:
first, the average outage frequency (SAIFI) and the expected low battery (EDNS) of the system are defined by equation (1):
Figure BDA0002282840570000031
in the formula, NIRepresenting a set of important distribution boards; lambda [ alpha ]jRepresenting the outage probability of panel j; cjRepresents the number of loads to which panel j is connected; pLjRepresenting the real load of panel j.
The reliability constraint of the ship power grid is embodied in that SAIFI and EDNS are smaller than given standard values, as shown in formula (2):
Figure BDA0002282840570000032
in the formula, SAIFI*Representing a given average outage frequency standard value; EDNS*A standard value representing a desired value of a given charge deficit.
And defining a line-panel association matrix by equation (3):
Figure BDA0002282840570000033
equation (4) represents a reliability index calculation model based on the line-distribution board correlation matrix:
Figure BDA0002282840570000034
in the formula (I), the compound is shown in the specification,
Figure BDA0002282840570000041
representing the failure rate of the mth power supply path of the power distribution board j; lambda [ alpha ]ijRepresents the failure rate of line ij; n is a radical ofm(j) A set of power supply paths representing a panel j;
Figure BDA0002282840570000042
the binary variable represents whether the distribution board j selects the mth power supply path or not; xjA binary variable representing whether the distribution board j is put into use; y isijThe input condition of the line ij is represented by a binary variable; n is a radical ofLRepresents a collection of distribution boards (meaning distribution boards) in a marine power grid; l denotes the line set in the vessel's power grid.
In the above method, the mathematical model for reconstructing the ship power grid fault comprises:
the goal of ship grid fault reconstruction is to recover as important loads as possible, and also to reduce the operating losses of the grid, since the ship grid needs to operate in a reconstructed state for a longer time.
The constraint conditions of the ship power grid fault reconstruction comprise active power balance constraint, reactive power balance constraint, voltage drop equality constraint, apparent power equality constraint, voltage safety operation constraint, line transmission capacity constraint, radiation network topology constraint and reliability constraint.
Maximizing the recovered load:
Figure BDA0002282840570000043
where S is a set of load levels including primary, secondary and tertiary loads, ωsIs the weight of each type of load.
And (3) minimizing the network loss:
Figure BDA0002282840570000044
in the formula iijIndicating lineCurrent at way ij; r isijRepresenting the resistance of line ij.
Overall objective:
Figure BDA0002282840570000045
the active power balance constraint is:
Figure BDA0002282840570000046
where w (j) is the parent panel set of panel j; pijIs the active power of line ij; pkjIs the active power of line kj; y iskjThe input condition of the line kj is represented by a binary variable; v (j) is a sub-panel set of panel j.
The reactive power balance constraint is:
Figure BDA0002282840570000051
in the formula, QijIs the reactive power of line ij; qkjIs the reactive power of line kj; x is the number ofijIs the impedance of line ij; qLjIs the reactive load of panel j.
The voltage drop equality constraints are:
Figure BDA0002282840570000052
in the formula ujIs the voltage of panel j.
The apparent power equation constrains:
Figure BDA0002282840570000053
voltage safety operation constraint:
Figure BDA0002282840570000054
in the formula uminAnd umaxIs a safe voltage boundary of a ship power system.
Line transmission capacity constraint:
Figure BDA0002282840570000055
in the formula (I), the compound is shown in the specification,
Figure BDA0002282840570000056
is the maximum transmission capacity of line ij.
And (3) topological constraint of the radiation network:
Figure BDA0002282840570000057
in the method, the second-order cone relaxation strategy of the ship power grid fault reconstruction model comprises the following steps:
converting the power flow model of the power distribution network into a second-order conical form:
Figure BDA0002282840570000058
Figure BDA0002282840570000059
wherein N is a distribution board set of the ship power grid; i isij、UjIs the square of the line current and panel voltage in the power distribution network flow model;
to consider whether a distribution line is switched on, the voltage is associated with the branch to which it is connected by the equation (17):
Figure BDA00022828405700000510
in the formula (I), the compound is shown in the specification,
Figure BDA00022828405700000511
is the square of the link voltage of panel j to line ij;
the objective function is converted into:
Figure BDA0002282840570000061
constraint (8) -formula (12) is transformed into formula (19) -formula (23):
Figure BDA0002282840570000062
Figure BDA0002282840570000063
Figure BDA0002282840570000064
in the formula (I), the compound is shown in the specification,
Figure BDA0002282840570000065
is the square of the link voltage of panel i to line ij;
Figure BDA0002282840570000066
in the formula (I), the compound is shown in the specification,
Figure BDA0002282840570000067
is the maximum current allowed for line ij.
Figure BDA0002282840570000068
To avoid islanding, the addition (24) serves as a constraint:
Figure BDA0002282840570000069
in the formula (I), the compound is shown in the specification,
Figure BDA00022828405700000610
line ij virtual active power, ξ is the virtual active load of panel j.
And (3) projecting the solution space onto the cone by the equation (25) to realize the relaxation of the power flow equation:
Figure BDA00022828405700000611
in the method, the ship power grid fault reconstruction convex optimization model comprises the following steps:
after the second-order cone is relaxed, the ship power grid fault reconstruction model is converted into a form of a formula (26) -a formula (27), an objective function of the ship power grid fault reconstruction model is a linear expression, constraint conditions comprise the linear expression and the second-order cone expression, a mixed integer second-order cone plan is formed, and the ship power grid fault reconstruction model is a convex optimization model.
Figure BDA0002282840570000071
Figure BDA0002282840570000072
Drawings
FIG. 1 is a flow chart provided by the present invention.
FIG. 2 is a schematic diagram of a statistical marine power system configuration according to the present invention;
Detailed Description
The invention is described in detail below with reference to specific embodiments and the accompanying drawings.
As shown in FIG. 1, the invention provides a ship power grid fault reconstruction convex optimization model considering reliability constraint, which comprises the following steps:
s1, dividing the load of the ship power grid into different important levels based on the reliability index concept and the ship power grid fault reconstruction requirement, and defining the important degrees of the different load levels;
and S11, dividing the ship power grid into a class I load, a class II load and a class III load according to the importance of the load.
I-level load: this load controls the vitality of the ship and the safety of the passengers belonging to the un-unloadable load; in any case a first order load must be provided. The primary loads of the ship system comprise a boiler, a generator, a weapon system, an electronic countermeasure system, a medical operating room and the like, and at least two power supply paths are required for connecting important loads. If a certain load is certified as a primary load (e.g. a ship propulsion system) in any one of the ship's operational tasks, it has to be connected to the ship's electrical system by Automatic Bus Transfer (ABT). Automatic bus transfer is a device that detects power loss from a normal power supply. When the normal power supply fails, the ABT can automatically and quickly disconnect the load from the normal power supply and connect the load power supply with the standby power supply. Non-critical loads have only one power supply path.
And II-stage load: this load is important to the proper operation of the vessel, but allows it to be offloaded or transferred to other platforms for power if necessary, in case of substantial loss of load. These loads include: cargo lifts, sea water pumps, some radar loads, etc.
Stage III load: the ship can be unloaded immediately when necessary, and the service life and the safety of the ship cannot be adversely affected; the method mainly comprises the following steps: residential heating systems, kitchen electrical loads, and the like.
S12, setting the importance degree of the loads with different importance degrees.
In the present invention, the importance level of the class I load is defined as 20, the importance level of the class II load is defined as 5, and the importance level of the class III load is defined as 1.
S2, establishing a reliability index calculation model through the line-distribution board incidence matrix;
and S21, defining the reliability constraint of the ship power grid.
First, the average outage frequency (SAIFI) and the expected low battery (EDNS) of the system are defined by equation (1):
Figure BDA0002282840570000091
in the formula, NIRepresenting a set of important distribution boards; lambda [ alpha ]jRepresenting the outage probability of panel j; cjRepresents the number of loads to which panel j is connected; pLjRepresenting the real load of panel j.
The reliability constraint of the ship power grid is embodied in that SAIFI and EDNS are smaller than given standard values, as shown in formula (2):
Figure BDA0002282840570000092
in the formula, SAIFI*Representing a given average outage frequency standard value; EDNS*A standard value representing a desired value of a given charge deficit.
And S22, obtaining a reliability index calculation model based on the line-distribution board incidence matrix.
First, a line-panel association matrix is defined by equation (3):
Figure BDA0002282840570000093
equation (4) represents a reliability index calculation model based on the line-distribution board correlation matrix:
Figure BDA0002282840570000101
in the formula (I), the compound is shown in the specification,
Figure BDA0002282840570000102
representing the failure rate of the mth power supply path of the power distribution board j; lambda [ alpha ]ijRepresents the failure rate of line ij; n is a radical ofm(j) A set of power supply paths representing a panel j;
Figure BDA0002282840570000103
the binary variable represents whether the distribution board j selects the mth power supply path or not; xjA binary variable representing whether the distribution board j is put into use; y isijThe input condition of the line ij is represented by a binary variable; n is a radical ofLRepresents a collection of distribution boards (meaning distribution boards) in a marine power grid; l denotes the line set in the vessel's power grid.
S3, establishing a ship power grid fault reconstruction model by taking the voltage constraint of a ship power grid, the system and the reliability constraint of important loads into consideration with the aim of reducing the load loss and the network loss of the ship power grid as targets;
and S31, establishing an objective function of the ship power grid fault reconstruction model.
The goal of ship grid fault reconstruction is to recover as important loads as possible, and also to reduce the operating losses of the grid, since the ship grid needs to operate in a reconstructed state for a longer time.
Maximizing the recovered load:
Figure BDA0002282840570000104
where S is a set of load levels including primary, secondary and tertiary loads, ωsIs the weight of each type of load.
And (3) minimizing the network loss:
Figure BDA0002282840570000111
in the formula iijRepresents the current of line ij; r isijRepresenting the resistance of line ij.
Overall objective:
Figure BDA0002282840570000112
and S32, establishing a constraint condition of the ship power grid fault reconstruction model.
The constraint conditions of the ship power grid fault reconstruction comprise active power balance constraint, reactive power balance constraint, voltage drop equality constraint, apparent power equality constraint, voltage safety operation constraint, line transmission capacity constraint, radiation network topology constraint and reliability constraint.
The active power balance constraint is:
Figure BDA0002282840570000113
where w (j) is the parent panel set of panel j;Pijis the active power of line ij; pkjIs the active power of line kj; y iskjThe input condition of the line kj is represented by a binary variable; v (j) is a sub-panel set of panel j.
The reactive power balance constraint is:
Figure BDA0002282840570000114
in the formula, QijIs the reactive power of line ij; qkjIs the reactive power of line kj; x is the number ofijIs the impedance of line ij; qijIs the reactive load of the panel.
The voltage drop equality constraints are:
Figure BDA0002282840570000115
in the formula ujIs the voltage of panel j.
The apparent power equation constrains:
Figure BDA0002282840570000116
voltage safety operation constraint:
Figure BDA0002282840570000117
in the formula uminAnd umaxIs a safe voltage boundary of a ship power system.
Line transmission capacity constraint:
Figure BDA0002282840570000118
in the formula (I), the compound is shown in the specification,
Figure BDA0002282840570000119
is the maximum transmission capacity of line ij.
And (3) topological constraint of the radiation network:
Figure BDA0002282840570000121
s4, converting the ship power grid fault reconstruction model into a convex optimization model through second-order cone relaxation;
s41, second-order cone relaxation strategy.
Converting the power flow model of the power distribution network into a second-order conical form:
Figure BDA0002282840570000122
Figure BDA0002282840570000123
wherein N is a distribution board set of the ship power grid; i isij、UjIs the square of the line current and panel voltage in the power distribution network flow model;
to consider whether a distribution line is switched on, the voltage is associated with the branch to which it is connected by the equation (17):
Figure BDA0002282840570000124
in the formula (I), the compound is shown in the specification,
Figure BDA0002282840570000125
is the square of the contact voltage of panel j to line ij.
The objective function is converted into:
Figure BDA0002282840570000126
constraint (8) -formula (12) is transformed into formula (19) -formula (23):
Figure BDA0002282840570000127
Figure BDA0002282840570000128
Figure BDA0002282840570000129
in the formula (I), the compound is shown in the specification,
Figure BDA00022828405700001210
is the square of the link voltage of panel i to line ij;
Figure BDA00022828405700001211
in the formula (I), the compound is shown in the specification,
Figure BDA00022828405700001212
is the maximum current allowed for line ij.
Figure BDA00022828405700001213
To avoid islanding, the addition (24) serves as a constraint:
Figure BDA0002282840570000131
in the formula (I), the compound is shown in the specification,
Figure BDA0002282840570000132
line ij virtual active power, ξ is the virtual active load of panel j.
And (3) projecting the solution space onto the cone by the equation (25) to realize the relaxation of the power flow equation:
Figure BDA0002282840570000133
and S42, reconstructing the convex optimization model of the ship power grid fault.
After the second-order cone is relaxed, the ship power grid fault reconstruction model is converted into a form of a formula (26) -a formula (27), an objective function of the ship power grid fault reconstruction model is a linear expression, constraint conditions comprise the linear expression and the second-order cone expression, a mixed integer second-order cone plan is formed, and the ship power grid fault reconstruction model is a convex optimization model.
Figure BDA0002282840570000134
Figure BDA0002282840570000141
Figure BDA0002282840570000142
And S5, solving to obtain the action switch which satisfies the constraint and minimizes the load loss.
The present invention is not limited to the above-mentioned preferred embodiments, and any structural changes made under the teaching of the present invention shall fall within the protection scope of the present invention, which has the same or similar technical solutions as the present invention.

Claims (9)

1. A ship power grid fault reconstruction method based on a convex optimization model is characterized in that the ship power grid fault reconstruction method considers reliability constraint and comprises the following steps:
step 1, dividing the load of a ship power grid into different important levels based on reliability indexes and ship power grid fault reconstruction requirements, and defining the important degrees of the different load levels;
step 2, according to the load classification result obtained in the step 1, establishing a reliability index calculation model through a line-distribution board incidence matrix;
step 3, taking the reliability index calculation model in the step 2 as a constraint condition, taking reduction of the loss load and the network loss of the ship power grid as a target, considering voltage constraint of the ship power grid, reliability constraint of a system and important loads, and establishing a ship power grid fault reconstruction model;
step 4, converting the ship power grid fault reconstruction model into a convex optimization model through second-order cone relaxation;
step 5, solving the convex optimization model to obtain an action switch which meets the constraint and minimizes the load loss of the ship power grid;
and 6, executing the action switch.
2. The method of claim 1, wherein defining a load level of a vessel's power grid comprises:
dividing a ship power grid into a class I load, a class II load and a class III load according to the importance of the load;
the I-level load is a first-level load of a ship system, each first-level load is connected with at least two power supply paths in a distributed mode, and each power supply path is connected to a ship electrical system through automatic bus transmission; the automatic bus transmission can automatically and quickly disconnect the load from the normal power supply and connect the load power supply with the standby power supply;
the II-level load is a secondary load of the ship system, and is allowed to be unloaded or transferred to other platforms for power supply when necessary, so that the load is prevented from being greatly lost;
the class III load is a three-level load of the ship system which can be unloaded immediately when necessary.
3. The method of claim 1, wherein considering the ship grid reliability constraints comprises:
first, the average outage frequency SAIFI and the expected low battery value EDNS of the system are defined by equation (1):
Figure FDA0002282840560000021
in the formula, NIRepresenting a set of important distribution boards; lambda [ alpha ]jRepresenting the outage probability of panel j; cjRepresents the number of loads to which panel j is connected; pLjRepresents the active load of panel j;
the reliability constraint of the ship power grid is embodied in that SAIFI and EDNS are less than or equal to given standard values, as shown in formula (2):
Figure FDA0002282840560000022
in the formula, SAIFI*Representing a given average outage frequency standard value; EDNS*A standard value representing a desired value of a given charge deficit.
4. The method according to claim 3, wherein the establishing of the reliability indicator calculation model specifically comprises the steps of:
first, a line-panel association matrix is defined by equation (3):
Figure FDA0002282840560000023
equation (4) represents a reliability index calculation model based on the line-distribution board correlation matrix:
Figure FDA0002282840560000031
in the formula (I), the compound is shown in the specification,
Figure FDA0002282840560000032
representing the failure rate of the mth power supply path of the power distribution board j; lambda [ alpha ]ijRepresents the failure rate of line ij; n is a radical ofm(j) A set of power supply paths representing a panel j;
Figure FDA0002282840560000033
the binary variable represents whether the distribution board j selects the mth power supply path or not; xjA binary variable representing whether the distribution board j is put into use; y isijThe input condition of the line ij is represented by a binary variable; n is a radical ofLRepresenting a set of distribution boards in a marine power grid; l denotes the line set in the vessel's power grid.
5. The method of claim 1, wherein the goal of reconstructing the ship grid fault reconstruction model further comprises reducing operational losses of the grid;
the reliability constraint conditions of the ship power grid fault reconstruction model comprise active power balance constraint, reactive power balance constraint, voltage drop equality constraint, apparent power equality constraint, voltage safety operation constraint, line transmission capacity constraint, radiation network topology constraint and reliability constraint.
6. The method according to claim 5, wherein the objective function of the ship grid fault reconstruction model is specifically:
maximizing the recovered load:
Figure FDA0002282840560000034
where S is a set of load levels including primary, secondary and tertiary loads, ωsIs the weight of each type of load.
And (3) minimizing the network loss:
Figure FDA0002282840560000041
in the formula iijRepresents the current of line ij; r isijRepresenting the resistance of line ij.
Overall objective:
Figure FDA0002282840560000042
7. the method of claim 5, wherein the vessel fault reconfiguration constraints are as follows:
the active power balance constraint is:
Figure FDA0002282840560000043
where w (j) is the parent panel set of panel j; pijIs the active power of line ij; pkjIs the active power of line kj; y iskjIs a binary systemA variable representing the input condition of the line kj; v (j) is a sub-panel set of panel j;
the reactive power balance constraint is:
Figure FDA0002282840560000044
in the formula, QijIs the reactive power of line ij; qkjIs the reactive power of line kj; x is the number ofijIs the impedance of line ij; qLjIs a reactive load of the panel;
the voltage drop equality constraints are:
Figure FDA0002282840560000045
in the formula ujIs the voltage of the panel j, uiIs the voltage of panel i;
the apparent power equation constrains:
Figure FDA0002282840560000046
voltage safety operation constraint:
Figure FDA0002282840560000047
in the formula uminAnd umaxIs a ship power system safe voltage boundary;
line transmission capacity constraint:
Figure FDA0002282840560000048
in the formula (I), the compound is shown in the specification,
Figure FDA0002282840560000049
is the maximum transmission capacity of line ij;
and (3) topological constraint of the radiation network:
Figure FDA0002282840560000051
8. the method of claim 1, wherein converting the power distribution network reconstruction model shown in the formulas (1) to (14) into a second-order cone planning model through second-order cone relaxation specifically comprises the steps of:
converting the power flow model of the power distribution network into a second-order conical form:
Figure FDA0002282840560000052
Figure FDA0002282840560000053
wherein N is a distribution board set of the ship power grid; i isij、UjIs the square of the line current and panel voltage in the power distribution network flow model;
to consider whether a distribution line is switched on, the voltage is associated with the branch to which it is connected by the equation (17):
Figure FDA0002282840560000054
in the formula (I), the compound is shown in the specification,
Figure FDA0002282840560000055
is the square of the link voltage of panel j to line ij;
the objective function is converted into:
Figure FDA0002282840560000056
constraint (8) -formula (12) is transformed into formula (19) -formula (23):
Figure FDA0002282840560000057
Figure FDA0002282840560000058
Figure FDA0002282840560000059
in the formula (I), the compound is shown in the specification,
Figure FDA00022828405600000510
is the square of the link voltage of panel i to line ij;
Figure FDA00022828405600000511
Figure FDA00022828405600000512
in the formula (I), the compound is shown in the specification,
Figure FDA00022828405600000513
is the maximum current allowed for line ij.
Figure FDA00022828405600000514
To avoid islanding, the addition (24) serves as a constraint:
Figure FDA00022828405600000515
in the formula (I), the compound is shown in the specification,
Figure FDA0002282840560000061
line ij virtual active power, ξ is the virtual active load of panel j.
And (3) projecting the solution space onto the cone by the equation (25) to realize the relaxation of the power flow equation:
Figure FDA0002282840560000062
9. the method of claim 8, wherein the ship power grid fault reconstruction convex optimization model after the second-order cone relaxation is constructed by the following steps:
after the second-order cone is relaxed, the ship power grid fault reconstruction model is converted into a form of a formula (26) -a formula (27), an objective function of the ship power grid fault reconstruction model is a linear expression, constraint conditions comprise the linear expression and the second-order cone expression, a mixed integer second-order cone plan is formed, and the ship power grid fault reconstruction model is a convex optimization model;
Figure FDA0002282840560000063
Figure FDA0002282840560000071
CN201911148323.4A 2019-11-21 2019-11-21 Ship power grid fault reconstruction convex optimization model considering reliability constraint Active CN111027184B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911148323.4A CN111027184B (en) 2019-11-21 2019-11-21 Ship power grid fault reconstruction convex optimization model considering reliability constraint

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911148323.4A CN111027184B (en) 2019-11-21 2019-11-21 Ship power grid fault reconstruction convex optimization model considering reliability constraint

Publications (2)

Publication Number Publication Date
CN111027184A true CN111027184A (en) 2020-04-17
CN111027184B CN111027184B (en) 2023-02-28

Family

ID=70206070

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911148323.4A Active CN111027184B (en) 2019-11-21 2019-11-21 Ship power grid fault reconstruction convex optimization model considering reliability constraint

Country Status (1)

Country Link
CN (1) CN111027184B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112271726A (en) * 2020-10-15 2021-01-26 北京交通大学 Power distribution system fault recovery method considering electricity-water-gas coupling relation
CN112653240A (en) * 2020-12-15 2021-04-13 贵州电网有限责任公司 Power distribution network protection configuration method considering network performance constraint
CN116381406A (en) * 2023-03-16 2023-07-04 武汉船舶职业技术学院 Ship power grid fault positioning method, device, equipment and readable storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130159309A1 (en) * 2009-06-25 2013-06-20 University Of Tennessee Research Foundation Method and apparatus for predicting object properties and events using similarity-based information retrieval and modeling
CN105243427A (en) * 2015-09-02 2016-01-13 南京航空航天大学 Aircraft power supply network dynamic programming management method
CN107844628A (en) * 2017-09-26 2018-03-27 上海电力学院 A kind of Large Scale Offshore Wind Farm collector system redundancy optimization method
CN109667677A (en) * 2018-12-26 2019-04-23 中国船舶重工集团公司第七〇九研究所 A kind of Poewr control method when Ship Electrical Power System diesel engine lubricating oil abnormal parameters
CN110350508A (en) * 2019-05-16 2019-10-18 东南大学 Method that is a kind of while considering the active distribution network fault recovery unified model that reconstruct is divided with isolated island

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130159309A1 (en) * 2009-06-25 2013-06-20 University Of Tennessee Research Foundation Method and apparatus for predicting object properties and events using similarity-based information retrieval and modeling
CN105243427A (en) * 2015-09-02 2016-01-13 南京航空航天大学 Aircraft power supply network dynamic programming management method
CN107844628A (en) * 2017-09-26 2018-03-27 上海电力学院 A kind of Large Scale Offshore Wind Farm collector system redundancy optimization method
CN109667677A (en) * 2018-12-26 2019-04-23 中国船舶重工集团公司第七〇九研究所 A kind of Poewr control method when Ship Electrical Power System diesel engine lubricating oil abnormal parameters
CN110350508A (en) * 2019-05-16 2019-10-18 东南大学 Method that is a kind of while considering the active distribution network fault recovery unified model that reconstruct is divided with isolated island

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王云帆，等: "基于逆变电源组网的电力系统短路电流改进算法", 《舰船科学技术》, 30 June 2018 (2018-06-30), pages 89 - 94 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112271726A (en) * 2020-10-15 2021-01-26 北京交通大学 Power distribution system fault recovery method considering electricity-water-gas coupling relation
CN112653240A (en) * 2020-12-15 2021-04-13 贵州电网有限责任公司 Power distribution network protection configuration method considering network performance constraint
CN112653240B (en) * 2020-12-15 2023-10-03 贵州电网有限责任公司 Power distribution network protection configuration method considering network performance constraint
CN116381406A (en) * 2023-03-16 2023-07-04 武汉船舶职业技术学院 Ship power grid fault positioning method, device, equipment and readable storage medium
CN116381406B (en) * 2023-03-16 2024-06-04 武汉船舶职业技术学院 Ship power grid fault positioning method, device, equipment and readable storage medium

Also Published As

Publication number Publication date
CN111027184B (en) 2023-02-28

Similar Documents

Publication Publication Date Title
CN111027184B (en) Ship power grid fault reconstruction convex optimization model considering reliability constraint
CN107862405B (en) Power system grid reconstruction optimization method taking microgrid as black-start power supply
CN108183545B (en) Multi-energy micro-grid power supply system for data center
CN109713704B (en) Communication interruption safety control method, system and medium for power grid side battery energy storage power station
WO2019096048A1 (en) Method and system for controlling energy consumption means of converter
CN111725829A (en) Hierarchical composite energy storage system for ship
CN110943525A (en) DC power supply system with optimal topology and intrinsic safety
CN117040037A (en) LNG (liquefied Natural gas) ship engine room monitoring power management system
CN113657619A (en) Key elastic lifting element identification and fault recovery method considering fault linkage
CN112865261A (en) Energy storage battery, application device thereof and charge-discharge control method
US20220166218A1 (en) Method of controlling a power distribution system including a microgrid
Qi et al. Post-disaster distribution system restoration considering uav-based communication recovery based on multi-agent reinforcement learning
CN114552579B (en) Power distribution network maximum power supply capacity calculation method considering low-voltage distribution area flexible interconnection
CN110867946A (en) Alternating current-direct current hybrid power supply integrated power supply
CN113572189B (en) Bipolar flexible direct current system for offshore wind power and transformer fault switching method thereof
CN113270871B (en) Flexible interconnection device capacity configuration optimization method based on power distribution network N-1 safety assessment
CN115513935A (en) Module fault bypass control method and device for high-voltage cascade energy storage system
CN110737872A (en) comprehensive evaluation method, device and readable storage medium for non-tree backbone net rack of system recovery after power failure
CN107453466B (en) Direct-current power supply system and control method thereof
Sennewald et al. Curative Actions by embedded bipolar HVDC-interconnections
CN215071655U (en) Energy storage system
TW202115987A (en) Method of transfer supply containing green energy for distribution feeder
Edström et al. Impact of remote control failure on power system restoration time
Wu et al. Service restoration of active distribution network considering the islanded operation of distributed generation and micro-grid
CN116073523B (en) Power supply system, low-voltage distribution line monitoring method and monitoring device thereof

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