A method and system for collaborative planning of urban distribution network and heating network

Abstract

The present invention belongs to the technical field of expansion planning of electric-thermal integrated energy systems, and provides a method and system for collaborative planning of urban distribution network-heat network, comprising the following steps: obtaining an urban distribution network-heat network model based on the same transmission network; considering the different needs of decision makers of the transmission network, the urban distribution network and the urban heat network, determining the objective function of the urban distribution network-heat network model; solving the determined objective function based on the alternating direction multiplier method of parallel splitting plus regularization term, and obtaining a collaborative planning scheme for the urban distribution network-heat network.

Description

Technical Field

[0001] The present invention belongs to the technical field of expansion planning of electric-thermal integrated energy systems, and specifically relates to a method and system for coordinated planning of urban distribution network-heat network. Background Art

[0002] The statements in this section merely provide background information related to the present disclosure and do not necessarily constitute prior art.

[0003] In recent years, with the gradual exposure of traditional energy pollution problems and the growing pressure of energy conservation and emission reduction, vigorously developing renewable energy power generation technology and replacing traditional power generation methods has become an unstoppable trend. The increase in the penetration rate of distributed renewable energy power generation in urban distribution networks has led to more complex and changeable transmission trends between distribution networks and transmission networks, and the physical coupling between the two has been enhanced. At the same time, the application of energy coupling devices such as cogeneration units and heat pumps has realized the interconnection between urban distribution networks and urban heating networks, and improved energy consumption efficiency. Based on this, in the process of urban energy network planning and design, it is of great significance to reasonably consider the coupling effect of the transmission network and realize the coordinated planning of multiple urban distribution networks and heating networks considering the influence of the transmission network.

[0004] When optimizing multiple networks, due to the expansion of system types and scales, traditional centralized algorithms face problems such as high computational pressure and the inability of centralized decision makers to obtain specific information of all networks. Distributed algorithms have gradually emerged because they can better solve these problems. As a type of distributed algorithm, the Alternative Direction Multiplier Method (ADMM) has gradually attracted the attention of researchers in the fields of distributed planning and operation of power systems in recent years due to its good convergence characteristics. However, for optimization problems with more than three operators, the direct generalization of ADMM has not yet given a convergence proof. Summary of the invention

[0005] In order to solve the above problems, the present invention proposes a method and system for collaborative planning of urban distribution network and heating network, establishes a network expansion planning model that comprehensively considers the interests of transmission network decision makers, urban distribution network decision makers and urban heating network decision makers, and designs a communication and collaboration scheme between different network decision makers based on the alternating direction multiplier method with parallel splitting plus regularization term (ADMM with Parallel SplittingAugmented Lagrangian Method, referred to as PSALM-ADMM), to ensure that each network decision maker takes into account the collaborative ability between different networks and the overall development needs of the network while planning and designing the internal network.

[0006] According to some embodiments, the first solution of the present disclosure provides a method for coordinated planning of urban distribution network-heating network, which adopts the following technical solution:

[0007] A method for collaborative planning of urban distribution network and heating network comprises the following steps:

[0008] Obtain a model of the urban distribution network-heating network based on the same transmission network;

[0009] Considering the different needs of decision makers of the transmission network, the urban distribution network and the urban heating network, an objective function of the urban distribution network-heating network model is determined;

[0010] Based on the parallel splitting plus regularization term alternating direction multiplier method, the determined objective function is solved and the coordinated planning scheme of urban distribution network and heating network is obtained.

[0011] As a further technical limitation, in the process of obtaining the urban distribution network-heating network model based on the same transmission network, the network topology structure, existing equipment status and historical operation data information of the transmission network, urban distribution network and urban heating network are collected, and based on the network node power balance equations, flow equations and equipment output upper and lower limit equations of the transmission network, urban distribution network and urban heating network, multiple urban distribution network-heating network mathematical models based on the same transmission network are constructed.

[0012] As a further technical limitation, the mathematical models of multiple urban distribution networks and heat networks based on the same transmission network include a transmission network model, a distribution network model and a heat network model.

[0013] As a further technical limitation, the objective function of the urban distribution network-heating network model is related to the transmission network operation cost, distribution network construction cost, distribution network operation cost, heating network construction cost and heating network operation cost.

[0014] As a further technical limitation, a regularization term is added to the prediction-correction method of the parallel split augmented Lagrangian to restrict it, thereby obtaining an alternating direction multiplier method of the parallel split plus regularization term.

[0015] As a further technical limitation, in the process of solving the determined objective function based on the alternating direction multiplier method with parallel splitting and regularization term, the coupling variables and Lagrange multipliers are initialized at the beginning of the iteration, the coupling variables and Lagrange multiplier information are updated during the iteration, the deviation in the iteration process is calculated, and the iteration ends when the maximum deviation is less than the allowable deviation, and the coordinated planning scheme of distribution network-heat network in multiple cities is output.

[0016] As a further technical limitation, during the distribution network-heating network collaborative planning process, decision makers of the transmission network, urban distribution network and urban heating network exchange network coupling node status information during the iterative update process.

[0017] According to some embodiments, the second solution of the present disclosure provides a city distribution network-heat network collaborative planning system, which adopts the following technical solutions:

[0018] A city distribution network-heating network collaborative planning system, comprising:

[0019] An acquisition module is configured to acquire a city distribution network-heat network model based on the same transmission network;

[0020] A building module is configured to consider different needs of decision makers of the transmission network, the urban distribution network and the urban heating network, and determine the objective function of the urban distribution network-heating network model;

[0021] The planning module is configured to solve the determined objective function based on the alternating direction multiplier method of parallel splitting plus regularization term to obtain a coordinated planning scheme of the urban distribution network-heating network.

[0022] According to some embodiments, a third solution of the present disclosure provides a computer-readable storage medium, which adopts the following technical solution:

[0023] A computer-readable storage medium stores a program, which, when executed by a processor, implements the steps in the urban distribution network-heating network collaborative planning method as described in the first aspect of the present disclosure.

[0024] According to some embodiments, a fourth solution of the present disclosure provides an electronic device, which adopts the following technical solution:

[0025] An electronic device comprises a memory, a processor and a program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps in the method for collaborative planning of an urban distribution network and a heating network as described in the first aspect of the present disclosure are implemented.

[0026] Compared with the prior art, the present invention has the following beneficial effects:

[0027] The PSALM-ADMM-based collaborative planning method for urban electricity and heat networks proposed in the present disclosure realizes the collaborative planning of multiple urban distribution networks and heat networks with coupling relationships. This method comprehensively considers the existing resources within multiple networks such as the transmission network, urban distribution network, and urban heat network, as well as the development needs of each network. Compared with the conventional independent planning of multiple networks, the proposed method can more effectively utilize the internal resources of multiple networks, reduce the redundant configuration of equipment capacity within multiple networks, and propose a more economical and resource-efficient collaborative planning solution. In the multi-network collaborative planning and design process, each network decision maker only needs to exchange the network coupling node status information during the iterative update process, which not only reduces the communication burden between network decision makers, but also ensures the privacy and security of the internal information of each network. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] The accompanying drawings constituting a part of the present disclosure are used to provide a further understanding of the present disclosure. The illustrative embodiments of the present disclosure and their descriptions are used to explain the present disclosure and do not constitute an improper limitation on the present disclosure.

[0029] Figure 1 is a flow chart of the urban distribution network-heating network collaborative planning method in the first embodiment of the present disclosure;

[0030] Figure 2 It is a diagram of the urban distribution network-heating network structure in the second embodiment of the present disclosure;

[0031] Figure 3 is a flowchart of the PSALM-ADMM algorithm in the first embodiment of the present disclosure;

[0032] Figure 4 It is a typical application system in the first embodiment of the present disclosure;

[0033] Figure 5 It is a structural block diagram of the urban distribution network-heating network collaborative planning system in the second embodiment of the present disclosure. DETAILED DESCRIPTION

[0034] The present disclosure is further described below in conjunction with the accompanying drawings and embodiments.

[0035] It should be noted that the following detailed descriptions are all illustrative and are intended to provide further explanation of the present disclosure. Unless otherwise specified, all technical and scientific terms used herein have the same meanings as those commonly understood by those skilled in the art to which the present disclosure belongs.

[0036] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present disclosure. As used herein, unless the context clearly indicates otherwise, the singular form is also intended to include the plural form. In addition, it should be understood that when the terms "comprising" and / or "including" are used in this specification, it indicates the presence of features, steps, operations, devices, components and / or combinations thereof.

[0037] In the absence of conflict, the embodiments in the present disclosure and the features in the embodiments may be combined with each other.

[0038] Embodiment 1

[0039] Embodiment 1 of the present disclosure introduces a method for collaborative planning of urban distribution network and heating network.

[0040] like Figure 1 A method for collaborative planning of urban distribution network and heat network is shown, comprising the following steps:

[0041] Obtain a model of the urban distribution network-heating network based on the same transmission network;

[0042] Considering the different needs of decision makers of the transmission network, the urban distribution network and the urban heating network, an objective function of the urban distribution network-heating network model is determined;

[0043] Based on the parallel splitting plus regularization term alternating direction multiplier method, the determined objective function is solved and the coordinated planning scheme of urban distribution network and heating network is obtained.

[0044] like Figure 2 As shown in the figure, a mathematical model of multiple urban distribution networks and heat networks connected to the same transmission network is constructed. Secondly, the overall objective function of the system is constructed by comprehensively considering the needs of each network decision maker. Then, the parallel splitting plus regularization term is introduced to construct the PSALM-ADMM algorithm. Finally, according to the proposed collaborative planning process based on PSALM-ADMM, each network decision maker transmits coupling variables and Lagrange multiplier information to each other, iteratively updates its own planning and operation plan until convergence, and obtains the collaborative planning plan of distribution networks and heat networks in multiple cities.

[0045] Mathematical model

[0046] The mathematical model of multiple urban distribution networks and heat networks connected to the same transmission network includes several parts, such as the transmission network model, the distribution network model, and the heat network model.

[0047] The transmission network model is:

[0048]

[0049]

[0050]

[0051]

[0052]

[0053] Where P L,i,s is the active power of the load at node i in scenario s; P ij,s is the active power flowing through line ij in scenario s; θ i,s is the phase angle of the voltage at node i in scenario s; b ij is the susceptance of line ij; is the maximum value of the active output of the generator at node i; is the maximum active power allowed to flow through line ij; is the maximum value of active power exchanged between the transmission network and the distribution network q; T b , T l , T g are the transmission network node set, line set, and generator set respectively; D11 It is the set of distribution network serial numbers connected to the transmission network.

[0054] Formula (1) is the node power balance equation, formula (2) is the line DC power flow equation; formula (3), formula (4) and formula (5) represent the upper and lower limit equations of the power exchange between generators, lines and transmission and distribution networks respectively.

[0055] The distribution network model is:

[0056]

[0057]

[0058]

[0059]

[0060]

[0061]

[0062]

[0063]

[0064] V min ≤V i,s ≤V max (14)

[0065]

[0066] In the formula, Db is the set of distribution network nodes; s is the substation node of the distribution network, that is, the set of nodes connected to the transmission network; V i,s represents the voltage amplitude of node i in typical scenario s; V ref is the reference node voltage; They represent the active and reactive power flowing through line ij in a typical scenario s respectively; They represent the active power generated by the distributed wind turbine and photovoltaic at node i in scenario s, and the active power obtained by the substation from the transmission grid. is the corresponding reactive power; represents the active power of the coupling node interacting with the heat desire at node i in scenario s; Respectively represent the active and reactive power base values ​​of the load at node i; They represent the normalized electric load, the maximum output of photovoltaic power, and the maximum output of wind power in the typical scenario s respectively; Z f,ijIt is a 0-1 variable, indicating whether there is a type f feeder at the distribution network line ij after planning; They represent the maximum capacity of f-type feeder and s-type substation respectively; is the maximum power of the cth coupling node between the distribution network and the heating network; b ij,s is an auxiliary slack variable. When the line at line ij is not built, b ij,s For [(V min -V max ),(V max -V min )], the auxiliary variable changes arbitrarily within the range, at this time, formula (18) is relaxed; r f 、x f They represent the resistance and reactance per unit length of type f feeder respectively; L ij Represents the length of line ij; V max 、V min Respectively represent the maximum and minimum values ​​of node voltage; They represent whether there is a forward or reverse virtual power flow on line ij.

[0067] Formula (6) is the linear Distflow power flow equation, which characterizes the relationship between the voltage at both ends of the line and the line power flow. Formula (7) and Formula (8) are the node active and reactive power balance equations respectively. Formula (9) is the auxiliary relaxation variable b ij,s The value range equation; Formula (10) ensures that there is only one type of feeder at line ij and only one type of substation at the substation node. Formula (11), Formula (12), Formula (13) and Formula (14) are the line capacity equation, substation capacity equation, new energy station output equation and node voltage upper and lower limit equation respectively. Formula (15) is a tree topology equation, which ensures that the planned distribution network is a tree distribution network.

[0068] The heat network model is:

[0069]

[0070]

[0071]

[0072]

[0073]

[0074]

[0075]

[0076] In the formula, Gb is the set of heating network nodes; Gb g A collection of nodes for heating network distribution stations; They represent the heat input from the heat distribution station and the electric-thermal coupling device at the node m in the scenario s respectively; represents the heat flowing through the m end of the line mn in the typical scenario s; δ is the heat loss rate per unit length of the pipeline; Represents the base value of heat load at node m; represents the normalized heat load in a typical scenario s; Z p,mn It is a 0-1 variable, indicating whether there is a p-type pipeline at the mn of the heating network line after planning; They represent the maximum flow of the p-type pipeline, the g-type heat distribution station, and the maximum power of the cogeneration device and heat pump in the c-th coupling node respectively; are the volume of natural gas consumed and the active power and heat generated by the CHP unit in the cth coupling node in the typical scenario s, respectively; are respectively the active power consumed and the heat generated by the heat pump in the cth coupling node in the typical scenario s; They are the electricity and heat production efficiencies of the cogeneration unit, respectively; is the heat generation efficiency of the heat pump; They respectively represent whether there is forward and reverse virtual heat in line mn.

[0077] Formula (16) is the node heat balance equation. Formula (17) is the pipeline heat loss equation. Formula (18) ensures that there is only one type of pipeline at line mn and only one type of gas distribution station at the gas distribution station node; Formula (19) is the upper and lower limit equation of equipment capacity; Formula (20) and Formula (21) are the operating equations of the cogeneration equipment and heat pump in the coupling node; Formula (22) is the tree topology equation.

[0078] Objective Function

[0079] The overall objective function of the system includes the transmission network operation cost Distribution network construction cost Distribution network operating costs Heat network construction cost Heating network operating costs As shown in formula (23):

[0080]

[0081] Transmission grid operating costs It can be expressed as:

[0082]

[0083] In the formula, H sis the number of hours that the typical scenario s occurs in a year; C G,i (g) is the power generation cost function of the generator at node i; P G,i,s is the active power output of the generator at node i in scenario s; C e is the unit electricity purchase cost; a G,i 、b G,i 、c G,i are the quadratic term, linear term coefficient and constant term of the generator cost function at node i respectively; Z on,i,s is a 0-1 variable representing the start / stop status of the generator at node i in scenario s.

[0084] Distribution network construction cost Distribution network operating costs It can be expressed as:

[0085]

[0086]

[0087] In the formula, Zn f,ij , Zc f,ij is a 0-1 variable, indicating whether there is a new or expanded f-type feeder at the distribution network line ij; Z s It is a 0-1 variable, indicating whether the distribution network substation is expanded to type s; Cn f , Cc f , C s They represent the construction and expansion costs of type f feeders per unit length and the expansion cost of type s substations respectively; L ij is the length of line ij; κ represents the annual investment cost conversion coefficient of the equipment; r represents the annual interest rate; T represents the planning period; Dl, Df, Ds, and Dc represent the line set, alternative feeder type set, alternative substation type set, and distribution-heating network coupling node set of the qth distribution network, respectively.

[0088] Heat network construction cost and heating network operating costs It can be expressed as:

[0089]

[0090] In the formula, Zn p,mn , Zc p,mn is a 0-1 variable, indicating whether there is a new or expanded p-type pipeline at the heating network line mn; Z g It is a 0-1 variable, indicating whether the heat distribution station of the heat network is expanded to type g; Cn p , Cc p , C gThey represent the construction and expansion costs of a unit length p-type pipeline and the expansion cost of a g-type heat distribution station; C ng , C n They are the unit gas purchase cost and heat purchase cost respectively; is the gas volume consumed by the CHP unit in the cth coupling node in the typical scenario s; F g,s The amount of heat purchased by the heat distribution station in the typical scenario s; L mn is the length of line mn; Gl, Gp, Gg represent the line set, alternative pipeline type set, and alternative gas distribution station type set of the qth heating network respectively; Dc represents the set of coupling nodes between the qth distribution network and the heating network.

[0091] PSALM-ADMM algorithm

[0092] For the single operator problem, x is the primal variable and λ is the dual variable in the variational inequality. In the algorithm flow, the k-th iteration of ALM starts from the given λ k Start with formula (28) to provide a new λ for the next iteration k+1 , x k+1 According to λ k The result obtained by calculation, so x is called the intermediate variable of the algorithm and λ is the core variable.

[0093]

[0094] For the three-operator problem, when yz separation is not considered, the method shown in formula (29) is convergent. When it is directly separated (called the direct generalization form of ADMM three operators), because the core variables include y, z, λ. y and z are the original variables, and λ is the dual variable. For the two original variables, the algorithm form is unfair. When solving the y subproblem, only the information of xk+1, zk, λk can be used, because zk has not been updated. When solving the z problem, y has been updated. There are examples where the calculation is successful, but there are also examples where it does not converge.

[0095]

[0096] For the (x, λ) problem, the inequality condition for the shrinkage of the core variable sequence λk is shown in equation (30). This equation guarantees the convergence of ALM.

[0097] ||λ k+1 -λ * || 2 ≤||λ k -λ * || 2 -||λ k -λ k+1 || 2 (30)

[0098] For the two-operator problem, this condition holds for both y and λ, so the two-operator optimization converges.

[0099] The cross-direction multiplication method, which is directly generalized, also has good computational effects for problems with more than three operators, but no convergence proof has been given so far. Therefore, some modified algorithms are proposed for the three-operator problem.

[0100] The prediction-correction algorithm of the partially parallel split ALM is as follows:

[0101]

[0102]

[0103] In this algorithm, the y and z subproblems are parallelized, but too freely, so they are corrected. In this embodiment, α=0.55 is taken.

[0104] The direct generalization relaxes the solution of the original part of the core variable, so when making up the overall, an appropriate reduction is added to ensure convergence. Based on the prediction and correction algorithm of the partial parallel split ALM, a regularization term is added to restrict it, and the ADMM algorithm of parallel split plus regularization term (PSALM) is obtained, as shown in the following formula:

[0105]

[0106] In the formula, τ is greater than 1.

[0107] The collaborative planning scheme process based on the PSALM-ADMM algorithm iteratively updates its own planning and operation scheme until convergence, and the specific process of obtaining the collaborative planning scheme of distribution networks and heating networks in multiple cities includes: first, initializing coupling variables, Lagrange multipliers and other information; second, solving the transmission network sub-problem, distribution network sub-problem and heating network sub-problem respectively; then, updating the Lagrange multipliers; finally, making a convergence judgment. If converged, the distribution network and heating network planning schemes are returned, otherwise the next iterative solution is performed.

[0108] like Figure 3 The planning method based on the PSALM-ADMM algorithm shown in FIG. 1 is as follows:

[0109] Step 1: Initialize variable P t,k-1 =P t,0 , P d,k =P d,0 , P h,k =P h,0 , Lagrange multiplier λ 1,0 =λ0,λ 2,0 =λ0,λ 3,0 =λ0, λ0 is the initial value, and the number of iterations k is set to 1.

[0110] Step 2: Place P d,k-1 , P h,k-1 Substitute it into the transmission network sub-problem, that is, solve formula (26) to get P t,k+1 .

[0111]

[0112] Step 3: P t,k+1 , P h,k-1 , bring it into the distribution network sub-problem, that is, solve formula (27) to get P d,k+1 .

[0113]

[0114] Step 4: P t,k+1 , P d,k-1 , bring it into the heat network problem, that is, solve formula (28) to get P h,k+1 .

[0115]

[0116] Step 5: Update the Lagrange multiplier λ according to formula (29) 1,k+1 ,λ 2,k+1 and λ 3,k+1 .

[0117] λ 1,k+1 =λ 1,k -β(P t,k+1 -η e P d,k+1 )

[0118] λ 2,k+1 =λ 2,k -β(P d,k+1 -η h P h,k+1 ) (37)

[0119] λ 3,k+1 =λ 3,k -β(P h,k+1 -η q P t,k+1 )

[0120] Step 6: Calculate the maximum deviation in k iterations using formula (30) and verify whether the maximum deviation is less than the allowable deviation. If the condition is met, terminate the iteration; otherwise, update the number of iterations k and return to step 2 to continue the iterative calculation.

[0121]

[0122] The typical system used in this embodiment is as follows Figure 4As shown in the figure, the transmission network is a dual-machine six-node system, connected to three urban distribution networks and a constant load. The topological structures of the three distribution-heating networks to be planned are the same, only the distributed power sources in the distribution networks are different. Only the structure of the distribution-heating network No. 1 to be planned is shown in detail in the figure. The distribution network to be planned is a 23-node system, with 4 existing lines, represented by black solid lines; 31 alternative lines, represented by gray dotted lines. The heat network to be planned is a 30-node system, with 10 existing pipelines, represented by black solid lines; 27 alternative lines, represented by gray dotted lines.

[0123] In the context of increasingly close coupling of energy networks in various cities, this embodiment proposes a system and method for collaborative planning of urban electricity-heat networks based on PSALM-ADMM. Based on PSALM-ADMM, a communication and collaboration scheme between different network decision makers is designed to ensure that each network decision maker takes into account the collaborative capabilities between different networks and the overall development needs of the network while planning and designing the internal network.

[0124] Embodiment 2

[0125] The second embodiment of the present disclosure introduces a city distribution network-heating network collaborative planning system.

[0126] like Figure 5 The urban distribution network-heating network collaborative planning system shown includes:

[0127] An acquisition module is configured to acquire a city distribution network-heat network model based on the same transmission network;

[0128] A building module is configured to consider different needs of decision makers of the transmission network, the urban distribution network and the urban heating network, and determine the objective function of the urban distribution network-heating network model;

[0129] The planning module is configured to solve the determined objective function based on the alternating direction multiplier method with parallel splitting plus regularization term to obtain the coordinated planning scheme of the urban distribution network-heating network

[0130] The detailed steps are the same as those of the urban distribution network-heating network collaborative planning method provided in Example 1, and will not be repeated here.

[0131] Embodiment 3

[0132] A third embodiment of the present disclosure provides a computer-readable storage medium.

[0133] A computer-readable storage medium stores a program, which, when executed by a processor, implements the steps in the urban distribution network-heating network collaborative planning method as described in the first embodiment of the present disclosure.

[0134] The detailed steps are the same as those of the urban distribution network-heating network collaborative planning method provided in Example 1, and will not be repeated here.

[0135] Embodiment 4

[0136] A fourth embodiment of the present disclosure provides an electronic device.

[0137] An electronic device comprises a memory, a processor and a program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps in the urban distribution network-heating network collaborative planning method as described in the first embodiment of the present disclosure are implemented.

[0138] The detailed steps are the same as those of the urban distribution network-heating network collaborative planning method provided in Example 1, and will not be repeated here.

[0139] Although the above describes the specific implementation methods of the present disclosure in conjunction with the accompanying drawings, it is not intended to limit the scope of protection of the present disclosure. Technical personnel in the relevant field should understand that on the basis of the technical solution of the present disclosure, various modifications or variations that can be made by those skilled in the art without creative work are still within the scope of protection of the present disclosure.

Claims

1. A method for collaborative planning of urban distribution network and heating network, characterized in that: The following steps are involved: Obtain a model of the urban distribution network-heating network based on the same transmission network; Considering the different needs of decision makers of the transmission network, the urban distribution network and the urban heating network, an objective function of the urban distribution network-heating network model is determined; Based on the parallel splitting plus regularization term alternating direction multiplier method, the determined objective function is solved to obtain the coordinated planning scheme of urban distribution network and heating network. On the basis of the prediction-correction method of parallel split augmented Lagrangian, regularization terms are added to restrict it, and the alternating direction multiplier method of parallel split plus regularization terms is obtained. In the process of solving the determined objective function based on the alternating direction multiplier method with parallel splitting and regularization term, the coupling variables and Lagrange multipliers are initialized at the beginning of the iteration, the coupling variables and Lagrange multiplier information are updated during the iteration, and the deviation in the iteration process is calculated until the maximum deviation is less than the allowable deviation. The iteration ends and the coordinated planning scheme of distribution network-heat network in multiple cities is output.

2. A method for coordinated planning of urban distribution network and heating network as claimed in claim 1, characterized in that: In the process of obtaining the urban distribution network-heating network model based on the same transmission network, the network topology structure, existing equipment status and historical operation data information of the transmission network, urban distribution network and urban heating network are collected, and according to the network node power balance equations, flow equations and equipment output upper and lower limit equations of the transmission network, urban distribution network and urban heating network, multiple urban distribution network-heating network mathematical models based on the same transmission network are constructed.

3. A method for coordinated planning of urban distribution network and heating network as claimed in claim 1, characterized in that: The mathematical model of multiple urban distribution networks-heating networks based on the same transmission network includes a transmission network model, a distribution network model and a heating network model.

4. A method for collaborative planning of urban distribution network and heating network as claimed in claim 1, characterized in that: The objective function of the urban distribution network-heating network model is related to the transmission network operation cost, distribution network construction cost, distribution network operation cost, heating network construction cost and heating network operation cost.

5. A method for coordinated planning of urban distribution network and heating network as claimed in claim 1, characterized in that: In the distribution network-heating network collaborative planning process, decision makers of the transmission network, urban distribution network and urban heating network exchange network coupling node status information during the iterative update process.

6. A city distribution network-heating network collaborative planning system, characterized in that: include: An acquisition module is configured to acquire a city distribution network-heat network model based on the same transmission network; A building module is configured to consider different needs of decision makers of the transmission network, the urban distribution network and the urban heating network, and determine the objective function of the urban distribution network-heating network model; A planning module is configured to solve the determined objective function based on the alternating direction multiplier method with parallel splitting plus regularization term to obtain a coordinated planning scheme of the urban distribution network-heating network; On the basis of the prediction-correction method of parallel split augmented Lagrangian, regularization terms are added to restrict it, and the alternating direction multiplier method of parallel split plus regularization terms is obtained. In the process of solving the determined objective function based on the alternating direction multiplier method with parallel splitting and regularization term, the coupling variables and Lagrange multipliers are initialized at the beginning of the iteration, the coupling variables and Lagrange multiplier information are updated during the iteration, and the deviation in the iteration process is calculated until the maximum deviation is less than the allowable deviation. The iteration ends and the coordinated planning scheme of distribution network-heat network in multiple cities is output.

7. A computer-readable storage medium having a program stored thereon, characterized in that: When the program is executed by a processor, the steps in the urban distribution network-heating network collaborative planning method as described in any one of claims 1 to 5 are implemented.

8. An electronic device comprising a memory, a processor, and a program stored in the memory and executable on the processor, characterized in that: When the processor executes the program, the steps in the urban distribution network-heating network collaborative planning method as described in any one of claims 1-5 are implemented.

Images