Distribution network physical form planning method, system, equipment and medium considering multi-dimensional reliability factors

Abstract

The present invention discloses a distribution network physical form planning method, system, equipment and medium taking into account multi-dimensional reliability factors, and belongs to the technical field of distribution network planning. The method includes: constructing a distribution network reliability evaluation index model taking into account operation risks; transforming the planning cost of distribution network reliability indicators; constructing a two-layer planning model of the distribution network physical form based on the reliability cost function; taking into account the KKT complementary relaxation conditions and simplifying the two-layer model based on the Big‑M method; solving the model to generate the physical form of distributed power sources and grids taking into account reliability. While improving the rationality of the planning of distributed power sources and grids in regional distribution networks, the present invention provides a more effective distribution network physical form planning method for the construction of diversified new power systems, meeting the reliability of power consumption of end users, and the safe operation of the system.

Description

Technical Field

[0001] The present invention relates to the technical field of distribution network planning, and in particular to a distribution network physical form planning method, system, equipment and medium taking into account multi-dimensional reliability factors. Background Art

[0002] Distribution network planning is to determine a scientific and reasonable distribution network construction plan by analyzing the load demand and the current status of the power grid in the planning area, so as to save construction funds while ensuring that the power demand of load users and the reliability of power grid operation are met. Most of the existing distribution network planning methods focus on a single planning approach. These models usually take economy as the main objective function, but ignore the user's demand for power supply reliability. Summary of the invention

[0003] In view of the deficiencies in the prior art, the present invention provides a distribution network physical form planning method, system, device and medium taking into account multi-dimensional reliability factors. The present invention first constructs a distribution network multi-dimensional reliability evaluation model taking into account the power supply of the distribution network system, the power consumption of users and the operation risks of nodes, converts the reliability index into one of the economic factors of distribution network planning in the form of cost, and converts it into a cost function in the distribution network double-layer planning model, then takes into account the KKT complementary relaxation condition and converts the double-layer nonlinear model into a single-layer linear model based on the Big-M method, and finally calls the Mathlab commercial solver for solving.

[0004] To achieve the above object, the present invention adopts the following technical solutions:

[0005] A distribution network physical form planning method taking into account multi-dimensional reliability factors comprises the following steps:

[0006] S1. Construct a distribution network reliability evaluation index model that considers operation risks from three dimensions: distribution network system end, distribution network user end, and distribution network node operation risk end;

[0007] S2. Based on the distribution network reliability evaluation index model, the reliability cost function of the system and user ends, the reliability cost function of the node voltage over-limit operation, and the reliability cost function of the node power over-limit operation are calculated; the reliability cost function of the system and user ends, the reliability cost function of the node voltage over-limit operation, and the reliability cost function of the node power over-limit operation are summed to obtain the comprehensive cost function of the distribution network reliability index;

[0008] S3. Use the comprehensive cost function of the distribution network reliability index to construct a two-layer planning model for the physical form of the distribution network; the upper layer of the two-layer planning model for the physical form of the distribution network is a distribution network grid expansion planning model, and the lower layer is a distribution network distributed power source site selection and capacity determination model; the upper layer takes the minimum annual comprehensive cost as the goal and provides a grid planning scheme for the lower layer; the lower layer takes the minimum annual comprehensive cost of the distribution network distributed power source planning as the goal, and selects the site and determines the capacity of the distributed power source based on the grid planning scheme provided by the upper layer, and passes it to the upper layer after determining the optimal placement address;

[0009] S4, taking into account the KKT complementary relaxation conditions, a two-level planning model of the physical form of the distribution network is simplified based on the Big-M method;

[0010] S5. Solve the simplified two-layer planning model of the physical form of the distribution network to generate the physical form of distributed generation and grid.

[0011] To optimize the above technical solutions, the specific measures taken also include:

[0012] Furthermore, S1 specifically includes:

[0013] S11. Select the reliability evaluation index of the distribution network system end, specifically:

[0014] The reliability evaluation index of the distribution network system end is taken as the first-level index, and the average power outage time and the average power supply reliability rate of the system are selected as the second-level indicators of the reliability evaluation of the distribution network system end;

[0015] System average power outage time T SAIDI The calculation formula is:

[0016]

[0017] System average power supply reliability rate f ASAI The calculation formula is:

[0018]

[0019] S12. Select the reliability evaluation index of the distribution network user end, specifically:

[0020] The reliability evaluation index of the user end of the distribution network is taken as the first-level index, and the average number of power outage users and the average power outage and power shortage of users are selected as the second-level indicators for the reliability evaluation of the user end of the distribution network;

[0021] Average number of users without power T MIC The calculation formula is as follows:

[0022]

[0023] Average power outage and power shortage for users W AENSThe calculation formula is as follows:

[0024]

[0025] S13. Construct reliability evaluation indicators for the risk end of distribution network node operation, specifically:

[0026] The node voltage over-limit rate and node power over-limit rate are selected as the reliability evaluation indicators of the risk end of the distribution network node operation. The node voltage over-limit rate f U The calculation formula is as follows:

[0027]

[0028] Where U 节点实际电压 Indicates the actual node voltage, U N Indicates the node rated voltage;

[0029] Node power excess rate f P The calculation formula is as follows:

[0030]

[0031] Where P 节点实际功率 Indicates the actual power of the node, P N Indicates the rated active power of the node.

[0032] Furthermore, S2 specifically includes:

[0033] S21. Perform planning cost conversion on the reliability evaluation index of the distribution network system end and the reliability evaluation index of the distribution network user end in the distribution network reliability evaluation index model to obtain the reliability cost function of the system end and the user end. The formula is as follows:

[0034] C1=b1p1W AENS T MIC +b2T SAIDI +b3f ASAI

[0035]

[0036]

[0037]

[0038] Among them, C1 represents the reliability cost function of the system and user sides, p1 represents the price of the unit power supply shortage at the user side under the average number of power outage users, b1 is the coefficient of the comprehensive planning cost of the user side reliability index, b2 is the cost conversion coefficient of the system average power outage time index, b3 is the cost conversion coefficient of the system average power supply possible index, r1 represents the expected value of the user side reliability evaluation index, r2 represents the expected value of the system average power outage time index, and r3 represents the expected value of the system average power supply possible index;

[0039] S22, the node voltage over-limit rate Ser (V n ) is used as the over-limit reliability index to transform into the coefficient of the node voltage operation over-limit reliability cost function, and the node voltage operation over-limit reliability cost function is obtained, and the formula is as follows:

[0040]

[0041] Where C2 represents the reliability cost function of node voltage overrun, T represents the set of all statistical periods, t represents the ordinal number of the statistical period, m represents the number of all nodes in the system, n represents the ordinal number of the node, and Ser(V n ) represents the overlimit rate of the nth node, P n represents the total active power of the nth node, c p Indicates the unit electricity price;

[0042] S23, the node power excess rate Ser (P n ) is used as the node power over-limit reliability index to transform it into the coefficient of the node power operation over-limit reliability cost function, and the node power operation over-limit reliability cost function is obtained, and the formula is as follows:

[0043]

[0044] Where C3 represents the reliability cost function of node power overrun, T represents the set of all statistical periods, t represents the ordinal number of the statistical period, m represents the number of all nodes in the system, n represents the ordinal number of the node, and Ser(P n ) represents the power excess rate of the nth node, S N Indicates the system rated apparent power, c s represents the power penalty cost coefficient;

[0045] S24. The reliability cost function of the system and user ends, the reliability cost function of the node voltage over-limit operation, and the reliability cost function of the node power over-limit operation are summed to obtain the comprehensive cost function of the reliability index of the distribution network. The formula is as follows:

[0046] C sev =C1+C2+C3

[0047] Among them C sev represents the comprehensive cost function of distribution network reliability indicators, C1 represents the reliability cost function of the system end and the user end, C2 represents the reliability cost function of the node voltage operation exceeding the limit, and C3 represents the reliability cost function of the node power operation exceeding the limit.

[0048] Further, in S3, the use of the comprehensive cost function of the distribution network reliability index to construct a two-level planning model for the physical form of the distribution network specifically includes:

[0049] S31. Construct an upper-level distribution network expansion planning model. The upper level takes the minimum annual comprehensive cost as the goal. The objective function is as follows:

[0050] min C=C line +C nl +C sev +C DG

[0051] Among them C line is the average annual cost of grid line construction, C nl is the annual network loss cost of the distribution network, C sev is the reliability index cost of distribution network, C DG Plan the annual comprehensive cost of distributed power generation for the lower distribution network;

[0052] The average annual cost of grid line construction C line The calculation formula is as follows:

[0053]

[0054] In the formula, r represents the discount rate, y represents the service life of the line, q represents the total number of new lines to be built, and x represents the total cost of the line. k represents the “0-1” decision variable of the kth line to be built, C k represents the investment cost of the kth line to be constructed;

[0055] Annual network loss cost of distribution network C nl The calculation formula is as follows:

[0056]

[0057] In the formula, T represents the set of all time periods, represents the network loss of the distribution network system during period t, represents the unit energy loss cost during period t;

[0058] S32. Construct the constraint conditions of the upper distribution network expansion planning model. The formula is as follows:

[0059]

[0060] Where P i-in The active power injected into node i by other nodes, P i-dis is the active power injected by the distributed generation at node i, P i-L is the active load at node i, Q i-in The reactive power injected into node i for other nodes, Q i-dis is the reactive power injected by the distributed generation at node i, Q i-L is the reactive load at node i, U i and U j are the voltages of the first and last nodes of any branch, G ij and B ij are the conductance and susceptance of the ij branch, respectively, and f ij is the power factor angle of the ij branch;

[0061] S33, construct the site selection and capacity determination model of the distributed power generation in the lower distribution network. The lower layer takes the minimum annual comprehensive cost of the distributed power generation planning of the distribution network as the goal, and the objective function is as follows:

[0062] min C DG =C om +C inv +C sev -F earn

[0063] In the formula, C om is the operation and maintenance cost of distributed power generation, C inv is the distributed generation investment cost, C sev is the reliability index cost of the distribution network, F earn Profit from distributed power generation;

[0064] Distributed power generation operation and maintenance cost C om The calculation formula is as follows:

[0065]

[0066] In the formula, T represents the set of all time periods, N DG represents the set of all distributed power sources, n DG represents the nth distributed power source, represents the power generation of the nth distributed generation at time t, It represents the operation and maintenance cost of the unit power generated by the nth distributed generation at time t;

[0067] Distributed power generation investment cost C inv The calculation formula is as follows:

[0068]

[0069] In the formula, r represents the discount rate, z represents the service life of the distributed power source, C DG-inv represents the fixed investment cost of distributed generation per unit capacity, S DG It is expressed as the installed capacity of distributed generation;

[0070] Distributed power generation revenue F earn The calculation formula is as follows:

[0071]

[0072] In the formula, T represents the set of all time periods, N DG represents the set of all distributed power sources, n DG represents the nth distributed power source, represents the power generation of the nth distributed generation at time t, c earn It represents the profit coefficient of distributed power generation per unit of electricity;

[0073] S34. Construct the constraint conditions of the site selection and capacity determination model of the distributed power generation in the lower distribution network, including the total installed capacity constraint of the distributed power generation in the system and the operating power constraint of the distributed power generation. The formula is as follows:

[0074] (1) Constraints on total installed capacity of distributed generation in the system

[0075]

[0076] Where: It is represented by the total installed capacity of distributed generation in the distribution network, Expressed as the maximum load of the distribution network, P PER It is the maximum penetration limit of distributed generation in distribution network.

[0077] (2) Distributed power generation operating power constraints

[0078]

[0079] Where: represents the total operating power of the distributed generation connected to node i, S i Represents the total capacity of distributed generation connected to node i.

[0080] Furthermore, the specific process of S4 is:

[0081] By constructing the lower-level Lagrangian function and based on the KKT complementary relaxation condition, the lower-level model is converted into the upper-level constraint condition, so that the physical form two-level programming model of the distribution network is transformed from a two-level nonlinear model to a single-level nonlinear model; the nonlinear terms contained in the single-level nonlinear model are converted into linear terms using the Big-M method, so that the physical form two-level programming model of the distribution network is transformed from a single-level nonlinear model to a single-level linear programming model.

[0082] Furthermore, the solution method used in S5 is: calling the solver CPLEX12.8 and YALMIP toolbox in Mathlab2016a to solve the simplified physical form two-level programming model of the distribution network.

[0083] Furthermore, in S22, the node voltage excess rate Ser (V n ) is obtained based on fuzzy membership fitting, and the expression is as follows:

[0084]

[0085] Where V n is the node voltage, and the upper and lower limits of the normal operating range are: V 1max 、V 1min , the abnormal normal operation state is: V 1min ≥V n 、V 1max ≤V n , where V 2min 、V 2max They are the upper and lower limits of the over-limit rate in abnormal operating state respectively;

[0086] In S23, the node power excess rate Ser (P n ) is obtained based on fuzzy membership fitting, and the expression is as follows:

[0087]

[0088] Where P n represents the node power, and Ser(Pn) represents the node power excess rate.

[0089] The present invention also proposes a distribution network physical form planning system taking into account multi-dimensional reliability factors, comprising:

[0090] Evaluation index construction module, used to construct a distribution network reliability evaluation index model that takes into account operation risks from three dimensions: distribution network system end, distribution network user end, and distribution network node operation risk end;

[0091] The planning cost conversion module is used to calculate the reliability cost function of the system and user ends, the reliability cost function of the node voltage operation exceeding the limit, and the reliability cost function of the node power operation exceeding the limit based on the distribution network reliability evaluation index model; the system and user end reliability cost function, the node voltage operation exceeding the limit reliability cost function, and the node power operation exceeding the limit reliability cost function are summed to obtain the comprehensive cost function of the distribution network reliability index;

[0092] A power grid physical form two-layer planning model construction module is used to construct a distribution network physical form two-layer planning model using a distribution network reliability index comprehensive cost function; the upper layer of the distribution network physical form two-layer planning model is a distribution network grid expansion planning model, and the lower layer is a distribution network distributed power source site selection and capacity determination model; the upper layer takes the minimum annual comprehensive cost as the goal and provides a grid planning scheme for the lower layer; the lower layer takes the minimum annual comprehensive cost of the distribution network distributed power source planning as the goal, and performs site selection and capacity determination of distributed power sources on the basis of the grid planning scheme provided by the upper layer, and passes the optimal placement address to the upper layer after determining it;

[0093] The power grid physical form two-level programming model simplification module is used to take into account the KKT complementary relaxation conditions and simplify the distribution network physical form two-level programming model based on the Big-M method;

[0094] The model solving module is used to solve the simplified double-layer planning model of the distribution network physical form to generate the physical form of distributed power sources and grid.

[0095] The present invention also proposes an electronic device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the distribution network physical form planning method taking into account multi-dimensional reliability factors as described above is implemented.

[0096] A computer-readable storage medium stores a computer program, wherein the computer program enables a computer to execute the distribution network physical form planning method taking into account multi-dimensional reliability factors as described above.

[0097] The beneficial effects of the present invention are:

[0098] The distribution network planning method proposed in the present invention not only takes into account the economic efficiency of planning and construction, but also improves the rationality of the distributed power source and grid planning of the regional distribution network, and takes into account the reliability of the system, users, and operation in an all-round and multi-level manner. On the basis of building a reliability evaluation model, the reliability index is converted into one of the economic factors of distribution network planning in the form of cost, and the economy and reliability are integrated into one, providing a more effective distribution network physical planning method to support the construction of diversified new power systems and meet the needs of end users for reliable power supply and safe operation of the system. BRIEF DESCRIPTION OF THE DRAWINGS

[0099] Figure 1 It is the overall flow chart of the method of the present invention;

[0100] Figure 2 It is a flow chart of constructing a distribution network reliability evaluation index model considering operation risks;

[0101] Figure 3 It is the cost conversion flow chart of distribution network reliability indicator planning;

[0102] Figure 4 It is a flow chart of constructing a two-level planning model of the physical form of the distribution network based on the reliability cost function;

[0103] Figure 5 It is a simplified flow chart of a two-layer model based on the Big-M method taking into account the KKT complementary relaxation conditions;

[0104] Figure 6 It is a regional grid structure diagram in a typical scenario of a new type of power system distribution network to which the present invention is applied;

[0105] Figure 7 It is a schematic diagram of the planning results of distributed power sources and grid physical form. DETAILED DESCRIPTION

[0106] The present invention will now be described in further detail with reference to the accompanying drawings.

[0107] In one embodiment, the present invention proposes a distribution network physical form planning method taking into account multi-dimensional reliability factors. The overall process is as follows: Figure 1 As shown, the following steps are included:

[0108] S1. Construct a distribution network reliability evaluation index model that considers operation risks from three dimensions: distribution network system end, distribution network user end, and distribution network node operation risk end; Figure 2 As shown, specifically including:

[0109] S11. Select the reliability evaluation index of the distribution network system end, specifically:

[0110] The reliability evaluation index of the distribution network system end in the "Power Supply Reliability Evaluation Regulations for Power Supply System" is used as the first-level index, and the average power outage time and the average power supply reliability rate of the system are selected as the second-level indicators for the reliability evaluation of the distribution network system end;

[0111] System average power outage time T SAIDI The calculation formula is:

[0112]

[0113] System average power supply reliability rate f ASAI The calculation formula is:

[0114]

[0115] S12. Select the reliability evaluation index of the distribution network user end, specifically:

[0116] The reliability evaluation index of the user end of the distribution network in the "Power Supply Reliability Evaluation Regulations for Power Supply System" is used as the first-level index, and the average number of power outage users and the average power outage and power shortage of users are selected as the second-level indicators for the reliability evaluation of the user end of the distribution network;

[0117] Average number of users without power T MIC The calculation formula is as follows:

[0118]

[0119] Average power outage and power shortage for users W AENS The calculation formula is as follows:

[0120]

[0121] S13. Construct reliability evaluation indicators for the risk end of distribution network node operation, specifically:

[0122] The node voltage over-limit rate and node power over-limit rate are selected as reliability evaluation indicators for the risk end of the node operation of the distribution network. In the distribution network, the size of the node voltage more intuitively represents the safety and reliability of the current distribution network operation. The over-limit of the node power will cause the distribution network to lose load, and in severe cases, it may cause the regional distribution network to be disconnected. Therefore, the node voltage over-limit rate and the node power over-limit rate are selected as the distribution network over-limit risk indicators.

[0123] Node voltage overlimit rate f U The calculation formula is as follows:

[0124]

[0125] Where U 节点实际电压 Indicates the actual node voltage, U N Indicates the node rated voltage;

[0126] Node power excess rate f P The calculation formula is as follows:

[0127]

[0128] Where P 节点实际功率 Indicates the actual power of the node, P N Indicates the rated active power of the node.

[0129] S2. Based on the distribution network reliability evaluation index model, the reliability cost function of the system and user ends, the reliability cost function of the node voltage over-limit operation, and the reliability cost function of the node power over-limit operation are calculated; the system and user end reliability cost function, the node voltage over-limit operation reliability cost function, and the node power over-limit operation reliability cost function are summed to obtain the comprehensive cost function of the distribution network reliability index; specifically including:

[0130] S21. According to the user's expectations for different reliability indicators and the opinions of relevant experts, the reliability evaluation indicators of the distribution network system end and the reliability evaluation indicators of the distribution network user end in the distribution network reliability evaluation indicator model are transformed into planning costs to obtain the reliability cost functions of the system end and the user end. The formula is as follows:

[0131] C1=b1p1W AENS T MIC +b2T SAIDI +b3f ASAI

[0132]

[0133]

[0134]

[0135] Among them, C1 represents the reliability cost function of the system and user sides, p1 represents the price of the unit power supply shortage at the user side under the average number of power outage users, b1 is the coefficient of the comprehensive planning cost of the user side reliability index, b2 is the cost conversion coefficient of the system average power outage time index, b3 is the cost conversion coefficient of the system average power supply possible index, r1 represents the expected value of the user side reliability evaluation index, r2 represents the expected value of the system average power outage time index, and r3 represents the expected value of the system average power supply possible index;

[0136] S22, based on fuzzy membership fitting, the relationship between node voltage and over-limit rate is obtained, such as Figure 3 As shown; get the node voltage and the excess rate Ser (V n ) is expressed as follows:

[0137]

[0138] Under the premise of ensuring the quality of power consumption for users, the node voltage is allowed to fluctuate within a small range and is considered to be in normal system operation. n is the node voltage, and the upper and lower limits of the normal operating range are: V 1max 、V 1min , the abnormal normal operation state is: V 1min ≥V n 、V 1max ≤V n , where V 2min 、V 2max They are the upper and lower limits of the over-limit rate in abnormal normal operation state. It can be seen that the node voltage has different over-limit rates in different ranges.

[0139] The node voltage excess rate Ser (Vn ) is used as the over-limit reliability index to transform into the coefficient of the node voltage operation over-limit reliability cost function, and the node voltage operation over-limit reliability cost function is obtained, and the formula is as follows:

[0140]

[0141] Where C2 represents the reliability cost function of node voltage overrun, T represents the set of all statistical periods, t represents the ordinal number of the statistical period, m represents the number of all nodes in the system, n represents the ordinal number of the node, and Ser(V n ) represents the overlimit rate of the nth node, P n represents the total active power of the nth node, c p Indicates the unit electricity price;

[0142] S23, about the relationship between node power size and over-limit rate, such as Figure 4 As shown;

[0143] Based on fuzzy membership fitting, the node power and over-limit rate Ser(P n ) is expressed as follows:

[0144]

[0145] The size of the node power reflects the quality of the system operation status in the area. A small range of fluctuations is considered to be a normal operating state under the premise of meeting the user's power quality requirements. In the formula: P n Expressed as node power, Ser(P n ) is expressed as the node power excess rate. It can be seen from the formula that the node power has different excess rates in different ranges.

[0146] The node power excess rate Ser(P n ) is used as the node power over-limit reliability index to transform it into the coefficient of the node power operation over-limit reliability cost function, and the node power operation over-limit reliability cost function is obtained, and the formula is as follows:

[0147]

[0148] Where C3 represents the reliability cost function of node power overrun, T represents the set of all statistical periods, t represents the ordinal number of the statistical period, m represents the number of all nodes in the system, n represents the ordinal number of the node, and Ser(P n ) represents the power excess rate of the nth node, S N Indicates the system rated apparent power, c s represents the power penalty cost coefficient;

[0149] S24. The reliability cost function of the system and user ends, the reliability cost function of the node voltage over-limit operation, and the reliability cost function of the node power over-limit operation are summed to obtain the comprehensive cost function of the reliability index of the distribution network. The formula is as follows:

[0150] C sev =C1+C2+C3

[0151] Among them C sev represents the comprehensive cost function of distribution network reliability indicators, C1 represents the reliability cost function of the system end and the user end, C2 represents the reliability cost function of the node voltage operation exceeding the limit, and C3 represents the reliability cost function of the node power operation exceeding the limit.

[0152] S3. Use the comprehensive cost function of the distribution network reliability index to construct a two-layer planning model for the physical form of the distribution network; the upper layer of the two-layer planning model for the physical form of the distribution network is the distribution network grid expansion planning model, and the lower layer is the distribution network distributed power site selection and capacity determination model; the upper layer takes the minimum annual comprehensive cost as the goal and provides a grid planning scheme for the lower layer; the lower layer takes the minimum annual comprehensive cost of the distribution network distributed power planning as the goal, and selects the site and determines the capacity of the distributed power based on the grid planning scheme provided by the upper layer, and passes it to the upper layer after determining the optimal placement address. Then the upper model optimizes the grid planning scheme again, taking into account the influencing factors of the distributed power access.

[0153] The construction of a two-level planning model for the physical form of a distribution network using the comprehensive cost function of distribution network reliability indicators specifically includes:

[0154] S31. Construct an upper-level distribution network expansion planning model. The upper level takes the minimum annual comprehensive cost as the goal. The objective function is as follows:

[0155] min C=C line +C nl +C sev +C DG

[0156] Among them C line is the average annual cost of grid line construction, C nl is the annual network loss cost of the distribution network, C sev is the reliability index cost of distribution network, C DG Plan the annual comprehensive cost of distributed power generation for the lower distribution network;

[0157] The average annual cost of grid line construction C line The calculation formula is as follows:

[0158]

[0159] In the formula, r represents the discount rate, y represents the service life of the line, q represents the total number of new lines to be built, and x represents the total cost of the line. k represents the “0-1” decision variable of the kth line to be built, C k represents the investment cost of the kth line to be constructed;

[0160] Annual network loss cost of distribution network C nl The calculation formula is as follows:

[0161]

[0162] In the formula, T represents the set of all time periods, represents the network loss of the distribution network system during period t, represents the unit energy loss cost during period t;

[0163] S32. Construct the constraint conditions of the upper distribution network expansion planning model. The formula is as follows:

[0164]

[0165] Where P i-in The active power injected into node i by other nodes, P i-dis is the active power injected by the distributed generation at node i, P i-L is the active load at node i, Q i-in The reactive power injected into node i for other nodes, Q i-dis is the reactive power injected by the distributed generation at node i, Q i-L is the reactive load at node i, U i and U j are the voltages of the first and last nodes of any branch, G ij and B ij are the conductance and susceptance of the ij branch, respectively, and f ij is the power factor angle of the ij branch;

[0166] S33, construct the site selection and capacity determination model of the distributed generation (DG) of the lower distribution network. The lower layer takes the minimum annual comprehensive cost of the distributed generation planning of the distribution network as the goal, and the objective function is as follows:

[0167] min C DG =C om +C inv +C sev -F earn

[0168] In the formula, C om is the operation and maintenance cost of distributed power generation, C inv is the distributed generation investment cost, C sev is the reliability index cost of the distribution network, F earnProfit from distributed power generation;

[0169] Distributed power generation operation and maintenance cost C om The calculation formula is as follows:

[0170]

[0171] In the formula, T represents the set of all time periods, N DG represents the set of all distributed power sources, n DG represents the nth distributed power source, represents the power generation of the nth distributed generation at time t, It represents the operation and maintenance cost of the unit power generated by the nth distributed generation at time t;

[0172] Distributed power generation investment cost C inv The calculation formula is as follows:

[0173]

[0174] In the formula, r represents the discount rate, z represents the service life of the distributed power source, C DG-inv represents the fixed investment cost of distributed generation per unit capacity, S DG It is expressed as the installed capacity of distributed generation;

[0175] Distributed power generation revenue F earn The calculation formula is as follows:

[0176]

[0177] In the formula, T represents the set of all time periods, N DG represents the set of all distributed power sources, n DG represents the nth distributed power source, represents the power generation of the nth distributed generation at time t, c earn It represents the profit coefficient of distributed power generation per unit of electricity;

[0178] S34. Construct the constraint conditions of the site selection and capacity determination model of the distributed power generation in the lower distribution network, including the total installed capacity constraint of the distributed power generation in the system and the operating power constraint of the distributed power generation. The formula is as follows:

[0179] (1) Constraints on total installed capacity of distributed generation in the system

[0180]

[0181] Where: It is represented by the total installed capacity of distributed generation in the distribution network, Expressed as the maximum load of the distribution network, P PERIt is the maximum penetration limit of distributed generation in distribution network.

[0182] (2) Distributed power generation operating power constraints

[0183]

[0184] Where: represents the total operating power of the distributed generation connected to node i, S i Represents the total capacity of distributed generation connected to node i.

[0185] S4. In the two-level programming model of the physical form of the distribution network constructed by the present invention based on the reliability cost function, the upper model contains two types of variables, continuous and integer, and there are nonlinear constraints in the lower layer. There is a coupling relationship between the upper and lower models, which is difficult to solve directly using traditional methods. The method proposed by the present invention is: taking into account the KKT complementary relaxation conditions, simplifying the two-level programming model of the physical form of the distribution network based on the Big-M method; the solution process is as follows Figure 5 The specific process is as follows:

[0186] By constructing the lower-level Lagrangian function and based on the KKT complementary relaxation condition, the lower-level model is converted into the upper-level constraint condition, so that the physical form two-level programming model of the distribution network is transformed from a two-level nonlinear model to a single-level nonlinear model; the nonlinear terms contained in the single-level nonlinear model are converted into linear terms using the Big-M method, so that the physical form two-level programming model of the distribution network is transformed from a single-level nonlinear model to a single-level linear programming model.

[0187] S5. Call the CPLEX12.8 solver and the YALMIP toolbox in Mathlab 2016a to solve the simplified two-level programming model of the physical form of the distribution network to generate the physical form of the distributed generation and the grid.

[0188] In another embodiment, the present invention provides a distribution network physical form planning system taking into account multi-dimensional reliability factors corresponding to the method of embodiment 1, including:

[0189] Evaluation index construction module, used to construct a distribution network reliability evaluation index model that takes into account operation risks from three dimensions: distribution network system end, distribution network user end, and distribution network node operation risk end;

[0190] The planning cost conversion module is used to calculate the reliability cost function of the system and user ends, the reliability cost function of the node voltage operation exceeding the limit, and the reliability cost function of the node power operation exceeding the limit based on the distribution network reliability evaluation index model; the system and user end reliability cost function, the node voltage operation exceeding the limit reliability cost function, and the node power operation exceeding the limit reliability cost function are summed to obtain the comprehensive cost function of the distribution network reliability index;

[0191] A power grid physical form two-layer planning model construction module is used to construct a distribution network physical form two-layer planning model using a distribution network reliability index comprehensive cost function; the upper layer of the distribution network physical form two-layer planning model is a distribution network grid expansion planning model, and the lower layer is a distribution network distributed power source site selection and capacity determination model; the upper layer takes the minimum annual comprehensive cost as the goal and provides a grid planning scheme for the lower layer; the lower layer takes the minimum annual comprehensive cost of the distribution network distributed power source planning as the goal, and performs site selection and capacity determination of distributed power sources on the basis of the grid planning scheme provided by the upper layer, and passes the optimal placement address to the upper layer after determining it;

[0192] The power grid physical form two-level programming model simplification module is used to take into account the KKT complementary relaxation conditions and simplify the distribution network physical form two-level programming model based on the Big-M method;

[0193] The model solving module is used to solve the simplified double-layer planning model of the distribution network physical form to generate the physical form of distributed power sources and grid.

[0194] The implementation methods of each module and the functions of the modules in the system are completely consistent with the steps of the method in Example 1, so they will not be repeated here.

[0195] In another embodiment, the present invention proposes an electronic device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the method for physical form planning of a distribution network taking into account multi-dimensional reliability factors as described in Embodiment 1 is implemented.

[0196] In another embodiment, the present invention provides a computer-readable storage medium storing a computer program, wherein the computer program enables a computer to execute the distribution network physical form planning method taking into account multi-dimensional reliability factors as described in the first embodiment.

[0197] A specific embodiment of a novel method for constructing a typical scenario physical form of a power system distribution network using the present invention is as follows:

[0198] Assume that the grid structure of the region is Figure 6 As shown: There are three power nodes (also called substation nodes), marked as nodes 1, 2 and 3, and their capacities are all 2*100MVA. The system currently contains 19 load nodes, and 5 new load nodes have been added, which are marked as nodes 12 to 16. There are 19 lines represented by solid lines in the system, and 12 new lines represented by dotted lines can be selected for construction. All lines use the same model, and their impedance parameters are 0.08+0.444jΩ / km. The thin solid lines in the system represent existing lines, the bold solid lines represent newly built extension lines, and the dotted lines represent interconnection lines. The dotted lines represent lines that have not yet been built. In the system, the gray nodes represent the locations where node generation (DG) is installed.

[0199] Enter the user's expected values ​​for all reliability indicators as shown in Tables 1 to 4:

[0200] Table 1

[0201]

[0202] Table 2

[0203]

[0204] Table 3

[0205]

[0206] Table 4

[0207] System Parameters Numeric discount rate 4 Line construction cost (10,000 yuan / km) 110 Photovoltaic cell investment cost (yuan / kW) 4500 Photovoltaic cell operation and maintenance cost (yuan / (kW·year)) 126 Unit network loss cost (yuan / kWh) 0.5 Unit power outage cost (yuan / kWh) 1.4 Unit electricity purchase cost (yuan / kWh) 0.34 Electricity price (yuan / kWh) 0.5 Return on investment 0.14

[0208] Finally, we get Figure 7 The planning results are shown.

[0209] In the embodiments disclosed in the present application, the computer storage medium can be a tangible medium that can contain or store a program for use by or in conjunction with an instruction execution system, device or equipment. The computer storage medium can include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices or equipment, or any suitable combination of the above. More specific examples of computer storage media can include electrical connections based on one or more lines, portable computer disks, hard disks, random access memories (RAM), read-only memories (ROM), erasable programmable read-only memories (EPROM or flash memory), optical fibers, portable compact disk read-only memories (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination of the above.

[0210] Those of ordinary skill in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed in this application can be implemented in electronic hardware or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this application.

[0211] The above are only preferred embodiments of the present invention, and the protection scope of the present invention is not limited to the above embodiments. All technical solutions under the concept of the present invention belong to the protection scope of the present invention. It should be pointed out that for ordinary technicians in this technical field, some improvements and modifications without departing from the principle of the present invention should be regarded as the protection scope of the present invention.

Claims

1. A distribution network physical form planning method taking into account multi-dimensional reliability factors, characterized in that: The steps include: S1. Construct a distribution network reliability evaluation index model that considers operation risks from three dimensions: distribution network system end, distribution network user end, and distribution network node operation risk end; S2. Based on the distribution network reliability evaluation index model, the reliability cost function of the system and user ends, the reliability cost function of the node voltage over-limit operation, and the reliability cost function of the node power over-limit operation are calculated; the reliability cost function of the system and user ends, the reliability cost function of the node voltage over-limit operation, and the reliability cost function of the node power over-limit operation are summed to obtain the comprehensive cost function of the distribution network reliability index; S2 specifically includes: S21. Perform planning cost conversion on the reliability evaluation index of the distribution network system end and the reliability evaluation index of the distribution network user end in the distribution network reliability evaluation index model to obtain the reliability cost function of the system end and the user end. The formula is as follows: C1=β1p1W AENS T MIC +β2T SAIDI +β3f ASAI Among them, C1 represents the reliability cost function of the system and user side, p1 represents the price of unit power supply shortage at the user side under the average number of power outage users, β1 is the coefficient of the comprehensive planning cost of the user side reliability index, and W AENS is the average power outage and power shortage of users, T MIC is the average number of users without power, β2 is the cost conversion coefficient of the average power outage time indicator of the system, T SAIDI is the average power outage time of the system, β3 is the cost conversion coefficient of the average power supply reliability index of the system, f ASAI is the average power supply reliability rate of the system, r1 is the expected value of the reliability evaluation index of the user end, r2 is the expected value of the average power outage time index of the system, and r3 is the expected value of the average power supply reliability rate index of the system; S22, the node voltage over-limit rate Ser (V n ) is used as the over-limit reliability index to transform into the coefficient of the node voltage operation over-limit reliability cost function, and the node voltage operation over-limit reliability cost function is obtained, and the formula is as follows: Where C2 represents the reliability cost function of node voltage overrun, T represents the set of all statistical periods, t represents the ordinal number of the statistical period, m represents the number of all nodes in the system, n represents the ordinal number of the node, and Ser(V n ) represents the voltage over-limit rate of the nth node, P n represents the total active power of the nth node, c p Indicates the unit electricity price; S23, the node power excess rate Ser (P n ) is used as the node power over-limit reliability index to transform it into the coefficient of the node power operation over-limit reliability cost function, and the node power operation over-limit reliability cost function is obtained, and the formula is as follows: Where C3 represents the reliability cost function of node power overrun, T represents the set of all statistical periods, t represents the ordinal number of the statistical period, m represents the number of all nodes in the system, n represents the ordinal number of the node, and Ser(P n ) represents the power excess rate of the nth node, S N Indicates the system rated apparent power, c s represents the power penalty cost coefficient; S24. The reliability cost function of the system and user ends, the reliability cost function of the node voltage over-limit operation and the reliability cost function of the node power over-limit operation are summed to obtain the comprehensive cost function of the reliability index of the distribution network. The formula is as follows: C sev =C1+C2+C3 Among them C sev represents the comprehensive cost function of distribution network reliability indicators, C1 represents the reliability cost function of the system end and the user end, C2 represents the reliability cost function of node voltage over-limit operation, and C3 represents the reliability cost function of node power over-limit operation; S3, using the comprehensive cost function of the distribution network reliability index to construct a two-layer planning model for the physical form of the distribution network; the upper layer of the two-layer planning model for the physical form of the distribution network is a distribution network grid expansion planning model, and the lower layer is a distribution network distributed power site selection and capacity determination model; the upper layer takes the minimum annual comprehensive cost as the goal and provides a grid planning scheme for the lower layer; the lower layer takes the minimum annual comprehensive cost of the distribution network distributed power planning as the goal, and selects the site and determines the capacity of the distributed power on the basis of the grid planning scheme provided by the upper layer, and passes it to the upper layer after determining the optimal placement address; in S3, the construction of the two-layer planning model for the physical form of the distribution network using the comprehensive cost function of the distribution network reliability index specifically includes: S31. Construct an upper-level distribution network expansion planning model. The upper level takes the minimum annual comprehensive cost as the goal. The objective function is as follows: minC=C line +C nl +C sev +C DG Among them C line is the average annual cost of grid line construction, C nl is the annual network loss cost of the distribution network, C sev is the reliability index cost of distribution network, C DG Plan the annual comprehensive cost of distributed power generation for the lower distribution network; The average annual cost of grid line construction C line The calculation formula is as follows: In the formula, r represents the discount rate, y represents the service life of the line, q represents the total number of new lines to be built, and x represents the total cost of the line. k represents the "0-1" decision variable of the kth line to be built, C k represents the investment cost of the kth line to be constructed; Annual network loss cost of distribution network C nl The calculation formula is as follows: In the formula, T represents the set of all time periods, W t nl represents the network loss of the distribution network system during period t, represents the unit energy loss cost during period t; S32. Construct the constraint conditions of the upper distribution network expansion planning model. The formula is as follows: Where P i-in The active power injected into node i by other nodes, P i-dis is the active power injected by the distributed generation at node i, P i-L is the active load at node i, Q i-in The reactive power injected into node i for other nodes, Q i-dis is the reactive power injected by the distributed generation at node i, Q i-L is the reactive load at node i, U i and U j are the voltages of the first and last nodes of any branch, G ij and B ij are the conductance and susceptance of the ij branch, respectively, ij is the power factor angle of the ij branch; S33, construct the site selection and capacity determination model of the distributed power generation in the lower distribution network. The lower layer takes the minimum annual comprehensive cost of the distributed power generation planning of the distribution network as the goal. The objective function is as follows: minC DG =C om +C inv +C sev -F earn In the formula, C om is the operation and maintenance cost of distributed power generation, C inv is the distributed generation investment cost, C sev is the reliability index cost of the distribution network, F earn Profit from distributed power generation; Distributed power generation operation and maintenance cost C om The calculation formula is as follows: In the formula, T represents the set of all time periods, N DG represents the set of all distributed power sources, n DG represents the nth distributed power source, represents the power generation of the nth distributed generation at time t, It represents the operation and maintenance cost of the unit power generated by the nth distributed generation at time t; Distributed power generation investment cost C inv The calculation formula is as follows: In the formula, r represents the discount rate, z represents the useful life of the distributed power source, C DG-inv represents the fixed investment cost of distributed generation per unit capacity, S DG It is expressed as the installed capacity of distributed generation; Distributed power generation revenue F earn The calculation formula is as follows: In the formula, T represents the set of all time periods, N DG represents the set of all distributed power sources, n DG represents the nth distributed power source, represents the power generation of the nth distributed generation at time t, c earn It represents the profit coefficient of distributed power generation per unit of electricity; S34. Construct the constraint conditions of the site selection and capacity determination model of the distributed power generation in the lower distribution network, including the total installed capacity constraint of the distributed power generation in the system and the operating power constraint of the distributed power generation. The formula is as follows: The total installed capacity constraint of distributed power generation in the system is: Where: It is represented by the total installed capacity of distributed generation in the distribution network, Expressed as the maximum load of the distribution network, P PER is the maximum penetration limit of distributed generation in the distribution network; The operating power constraint of distributed generation is: Where: represents the total operating power of the distributed power connected to node i, S i represents the total capacity of distributed power sources connected to node i; S4, taking into account the KKT complementary relaxation conditions, a two-level planning model of the physical form of the distribution network is simplified based on the Big-M method; S5. Solve the simplified two-layer planning model of the physical form of the distribution network to generate the physical form of distributed power sources and grid.

2. The method for physical form planning of a distribution network taking into account multi-dimensional reliability factors according to claim 1, characterized in that: S1 specifically includes: S11. Select the reliability evaluation index of the distribution network system end, specifically: The reliability evaluation index of the distribution network system end is taken as the first-level index, and the average power outage time and the average power supply reliability rate of the system are selected as the second-level indexes of the reliability evaluation of the distribution network system end; System average power outage time T SAIDI The calculation formula is: System average power supply reliability rate f ASAI The calculation formula is: S12. Select the reliability evaluation index of the distribution network user end, specifically: The reliability evaluation index of the user end of the distribution network is taken as the first-level index, and the average number of power outage users and the average power outage and power shortage of users are selected as the second-level indicators for the reliability evaluation of the user end of the distribution network; Average number of users without power T MIC The calculation formula is as follows: Average power outage and power shortage for users W AENS The calculation formula is as follows: S13. Construct reliability evaluation indicators for the risk end of distribution network node operation, specifically: The node voltage excess rate and node power excess rate are selected as reliability evaluation indicators for the risk end of distribution network node operation.

3. The method for physical form planning of a distribution network taking into account multi-dimensional reliability factors according to claim 1, characterized in that: The specific process of S4 is: By constructing the lower-level Lagrangian function and based on the KKT complementary relaxation condition, the lower-level model is converted into the upper-level constraint condition, so that the physical form two-level planning model of the distribution network is transformed from a two-level nonlinear model to a single-level nonlinear model. The Big-M method is used to convert the nonlinear terms contained in the single-layer nonlinear model into linear ones, so that the physical form two-layer programming model of the distribution network is transformed from a single-layer nonlinear model to a single-layer linear programming model.

4. The method for physical form planning of a distribution network taking into account multi-dimensional reliability factors according to claim 1, characterized in that: The solution method used in S5 is to call the solver CPLEX12.8 and YALMIP toolbox in Mathlab 2016a to solve the simplified physical form two-level programming model of the distribution network.

5. The method for physical form planning of distribution network taking into account multi-dimensional reliability factors according to claim 1, characterized in that: In S22, the node voltage excess rate Ser (V n ) is obtained based on fuzzy membership fitting, and the expression is as follows: Where V n is the node voltage, and the upper and lower limits of the normal operating range are: V 1max 、V 1min , the abnormal normal operation state is: V 1min ≥V n 、V 1max ≤V n , where V 2min 、V 2max They are the upper and lower limits of the over-limit rate in abnormal operating state respectively; In S23, the node power excess rate Ser (P n ) is obtained based on fuzzy membership fitting, and the expression is as follows: Where P n represents the node power, and Ser(Pn) represents the node power excess rate.

6. A distribution network physical form planning system taking into account multi-dimensional reliability factors, characterized in that: include: Evaluation index construction module, used to construct a distribution network reliability evaluation index model that takes into account operation risks from three dimensions: distribution network system end, distribution network user end, and distribution network node operation risk end; The planning cost conversion module is used to calculate the reliability cost function of the system and user ends, the reliability cost function of the node voltage operation exceeding the limit, and the reliability cost function of the node power operation exceeding the limit based on the distribution network reliability evaluation index model; the system and user end reliability cost function, the node voltage operation exceeding the limit reliability cost function, and the node power operation exceeding the limit reliability cost function are summed to obtain the comprehensive cost function of the distribution network reliability index; the formula of the system and user end reliability cost function is as follows: C1=β1p1W AENS T MIC +β2T SAIDI +β3f ASAI Among them, C1 represents the reliability cost function of the system and user side, p1 represents the price of unit power supply shortage at the user side under the average number of power outage users, β1 is the coefficient of the comprehensive planning cost of the user side reliability index, and W AENS is the average power outage and power shortage of users, T MIC is the average number of users without power, β2 is the cost conversion coefficient of the average power outage time indicator of the system, T SAIDI is the average power outage time of the system, β3 is the cost conversion coefficient of the average power supply reliability index of the system, f ASAI is the average power supply reliability rate of the system, r1 is the expected value of the reliability evaluation index of the user end, r2 is the expected value of the average power outage time index of the system, and r3 is the expected value of the average power supply reliability rate index of the system; The node voltage over-limit reliability cost function is specifically: the node voltage over-limit rate Ser (V n ) is used as the over-limit reliability index to transform into the coefficient of the node voltage operation over-limit reliability cost function, and the node voltage operation over-limit reliability cost function is obtained, and the formula is as follows: Where C2 represents the reliability cost function of node voltage overrun, T represents the set of all statistical periods, t represents the ordinal number of the statistical period, m represents the number of all nodes in the system, n represents the ordinal number of the node, and Ser(V n ) represents the voltage over-limit rate of the nth node, P n represents the total active power of the nth node, c p Indicates the unit electricity price; The node power operation over-limit reliability cost function is specifically: the node power over-limit rate Ser (P n ) is used as the node power over-limit reliability index to transform it into the coefficient of the node power operation over-limit reliability cost function, and the node power operation over-limit reliability cost function is obtained, and the formula is as follows: Where C3 represents the reliability cost function of node power overrun, T represents the set of all statistical periods, t represents the ordinal number of the statistical period, m represents the number of all nodes in the system, n represents the ordinal number of the node, and Ser(P n ) represents the power excess rate of the nth node, S N Indicates the system rated apparent power, c s represents the power penalty cost coefficient; The comprehensive cost function of distribution network reliability index is obtained by summing the reliability cost functions of the system and user ends, the reliability cost function of the node voltage over-limit operation, and the reliability cost function of the node power over-limit operation. The formula is as follows: C sev =C1+C2+C3 Among them C sev represents the comprehensive cost function of the reliability index of the distribution network, C1 represents the reliability cost function of the system end and the user end, C2 represents the reliability cost function of the node voltage operation exceeding the limit, and C3 represents the reliability cost function of the node power operation exceeding the limit; a power grid physical form double-layer planning model construction module is used to construct a distribution network physical form double-layer planning model using the comprehensive cost function of the distribution network reliability index; the upper layer of the distribution network physical form double-layer planning model is the distribution network grid expansion planning model, and the lower layer is the distribution network distributed power site selection and capacity determination model; the upper layer takes the minimum annual comprehensive cost as the goal and provides a grid planning scheme for the lower layer; the lower layer takes the minimum annual comprehensive cost of the distribution network distributed power planning as the goal, and selects the site and determines the capacity of the distributed power on the basis of the grid planning scheme provided by the upper layer, and passes it to the upper layer after determining the optimal placement address; the construction of the distribution network physical form double-layer planning model using the comprehensive cost function of the distribution network reliability index specifically includes: S31. Construct an upper-level distribution network expansion planning model. The upper level takes the minimum annual comprehensive cost as the goal. The objective function is as follows: minC=C line +C nl +C sev +C DG Among them C line is the average annual cost of grid line construction, C nl is the annual network loss cost of the distribution network, C sev is the reliability index cost of distribution network, C DG Plan the annual comprehensive cost of distributed power generation for the lower distribution network; The average annual cost of grid line construction C line The calculation formula is as follows: In the formula, r represents the discount rate, y represents the service life of the line, q represents the total number of new lines to be built, and x represents the total cost of the line. k represents the "0-1" decision variable of the kth line to be built, C k represents the investment cost of the kth line to be constructed; Annual network loss cost of distribution network C nl The calculation formula is as follows: In the formula, T represents the set of all time periods, W t nl represents the network loss of the distribution network system during period t, represents the unit energy loss cost during period t; S32. Construct the constraint conditions of the upper distribution network expansion planning model. The formula is as follows: Where P i-in The active power injected into node i by other nodes, P i-dis is the active power injected by the distributed generation at node i, P i-L is the active load at node i, Q i-in The reactive power injected into node i for other nodes, Q i-dis is the reactive power injected by the distributed generation at node i, Q i-L is the reactive load at node i, U i and U j are the voltages of the first and last nodes of any branch, G ij and B ij are the conductance and susceptance of the ij branch, respectively, ij is the power factor angle of the ij branch; S33, construct the site selection and capacity determination model of the distributed power generation in the lower distribution network. The lower layer takes the minimum annual comprehensive cost of the distributed power generation planning of the distribution network as the goal. The objective function is as follows: minC DG =C om +C inv +C sev -F earn In the formula, C om is the operation and maintenance cost of distributed power generation, C inv is the distributed generation investment cost, C sev is the reliability index cost of the distribution network, F earn Profit from distributed power generation; Distributed power generation operation and maintenance cost C om The calculation formula is as follows: In the formula, T represents the set of all time periods, N DG represents the set of all distributed power sources, n DG represents the nth distributed power source, represents the power generation of the nth distributed generation at time t, It represents the operation and maintenance cost of the unit power generated by the nth distributed generation at time t; Distributed power generation investment cost C inv The calculation formula is as follows: In the formula, r represents the discount rate, z represents the useful life of the distributed power source, C DG-inv represents the fixed investment cost of distributed generation per unit capacity, S DG It is expressed as the installed capacity of distributed generation; Distributed power generation revenue F earn The calculation formula is as follows: In the formula, T represents the set of all time periods, N DG represents the set of all distributed power sources, n DG represents the nth distributed power source, represents the power generation of the nth distributed generation at time t, c earn It represents the profit coefficient of distributed power generation per unit of electricity; S34. Construct the constraint conditions of the site selection and capacity determination model of the distributed power generation in the lower distribution network, including the total installed capacity constraint of the distributed power generation in the system and the operating power constraint of the distributed power generation. The formula is as follows: The total installed capacity constraint of distributed power generation in the system is: Where: It is represented by the total installed capacity of distributed generation in the distribution network, Expressed as the maximum load of the distribution network, P PER is the maximum penetration limit of distributed generation in the distribution network; The operating power constraint of distributed generation is: Where: represents the total operating power of the distributed power connected to node i, S i represents the total capacity of distributed power sources connected to node i; The power grid physical form two-level programming model simplification module is used to take into account the KKT complementary relaxation conditions and simplify the distribution network physical form two-level programming model based on the Big-M method; The model solving module is used to solve the simplified double-layer planning model of the distribution network physical form to generate the physical form of distributed power sources and grid.

7. An electronic device, characterized in that: include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the method for physical form planning of a distribution network taking into account multi-dimensional reliability factors as described in any one of claims 1 to 5 is implemented.

8. A computer-readable storage medium, characterized in that: A computer program is stored, and the computer program enables a computer to execute the distribution network physical form planning method taking into account multi-dimensional reliability factors as described in any one of claims 1 to 5.

Images