Dynamic Optimal Power Flow of Active Distribution Network Based on LSOCR and Its Application Scenarios

: Optimal power ﬂow (OPF) is a crucial aspect of distribution network planning and operation. Conventional heuristic algorithms fail to meet the system requirements for speed and accuracy, while linearized OPF approaches are inadequate for distribution networks with high R/X ratios. To address these issues and cater to multi-period scenarios, this study proposes a dynamic linearized second-order cone programming-based (SOCP) OPF model. The model is built by ﬁrst establishing a dynamic OPF model based on linearized second-order conic relaxation (LSOCR-DOPF). The components of the active distribution network, such as renewable energy power generation units, energy storage units, on-load-tap-changers, static var compensators, and capacitor banks, are then separately modeled. The model is implemented in MATLAB and solved by YALMIP and GUROBI. Finally, three representative scenarios are used to evaluate the model accuracy and effectiveness. The results show that the proposed LSOCR-DOPF model can ensure calculation time within 3 min, voltage stability, and error control within 10 − 6 for all three applications. This method has strong practical value in the ﬁelds of active distribution network day-ahead dispatch, accurate modeling of ZIP load, and real-time operation.


Introduction
Recently, researchers have shown an increased interest in the active distribution network.The integration of various distributed generations, energy storage units, and active management devices has presented new challenges to the planning and operation of distribution networks [1], especially in the field of active management (AM) of distribution networks [2].It is particularly urgent to develop optimization algorithms and high-performance computing tools applicable to various fields of active distribution networks.Ref. [3] analyzed three kinds of optimization problems of the smart grid: optimal power flow (OPF), unit commitment, and operation planning.Their essence is distribution network optimization, while having different optimization scales.The OPF is of great significance in the development process of the distribution network, and is the most common and fundamental optimization problem in power systems [4].Research on distribution network OPF has mainly focused on the alternating current power flow (AC-OPF).Exploring a solution method to enhance the solution speed of distribution network AC-OPF while ensuring its optimal operation and fulfilling the requirements of active distribution network planning and operation has been a major concern in the field of power system research.As a non-convex optimization problem, the OPF is difficult to solve.It is easy to fall into local optimum in the process of solving and has been proven to be an NP-hard problem [5,6].The power flow constraints are characterized by nonlinearity, and hence the essence of OPF lies in nonlinear programming.The methods for solving OPF problems can be broadly categorized into the following three categories: (1) Pursuit of local optimal solutions, including early classical methods [7] (such as Simplified Gradient Method, Newton Method, Sequential Quadratic Programming, Interior Point Method) and recent rapidly-developing heuristic algorithms [3,8]; (2) Approximation of power flow equality constraints.For example, the AC-OPF constraints can be approximately linearized as direct current power flow constraints, and the resultant direct current optimal power flow (DC-OPF) problem can be accordingly solved [9]; (3) Relaxation of the power flow equality constraint using convex relaxation techniques [10].
Methods of seeking local optimal solutions, due to their advantages such as simplicity and ease in simulating complex constraints, have been widely applied in solving nonlinear programming models [11].However, due to the non-convex nature of the OPF problem, these methods cannot guarantee the quality of the solution, and it is impossible to measure the gap between the locally optimal solution and the globally optimal solution.Methods that approximate the power flow equation constraints, such as the DC-OPF, present several obvious disadvantages.Firstly, they can be challenging to apply to research areas related to voltage and reactive power, as well as distribution networks with high R/X ratios.Secondly, the optimal solution of the DC-OPF problem may not be a feasible solution of the original OPF problem, leading to the need for constant adjustment of the tightness of DC-OPF constraints and the need to solve again during the actual optimization process [12].With the numerous issues associated with seeking local optima and approximating the current flow equality constraints, it requires new methods to deal with the current flow equality constraints.
The convex relaxation technique has gained significant attention in recent years due to its advantages and potential in solving OPF problems in the field of power system optimization.The use of convex relaxation techniques, particularly second-order cone programming (SOCP) relaxation, has become increasingly prevalent.The SOCP mainly converts the original model into a convex programming form, thus obtaining the globally optimal solution and a good computational speed.Refs.[13,14] systematically established a branch flow model (BFM) based on Distflow [15] to solve the OPF model framework, and presented two relaxation steps: (1) Elimination of voltage and current phase angles; (2) Second Order Conic Relaxation (SOCR).The authors also demonstrated the relaxation accuracy of SOCR.As a typical representative of SOCP, SOCR can be summarized in the mathematical field as the classic "dimensionality relaxation-return mapping" process: A new variable is introduced to elevate the dimensionality of the original problem; then, in the elevated problem, non-convex constraints are relaxed and a solution is obtained; finally, the solution of the original problem is recovered through a return mapping.The main challenge in using SOCR to solve the optimal power flow lies in accurately satisfying the conditions of the relaxed model.Refs.[16,17] provided a comprehensive summary of the sufficient conditions for the accuracy of SOCR under radiated network conditions, which are divided into three categories: power injection constraints, voltage amplitude constraints, and node voltage phase angle deviations, and corresponding explanations are given for the optimal power flow objective function.After the examination of Ref. [17], further research on the exact relaxation sufficient condition has been extensively conducted.
Ref. [18] expanded upon the sufficient conditions for the accuracy of SOCR in the presence of the high penetration of distributed generation.Refs.[19,20] highlighted the shortcomings of the second sufficient condition in Ref. [17], as it neglected the impact of the grounding branch power flow on line capacity, and proposed an improved sufficient condition.
Additionally, although the OPF model of the active distribution network based on MISOCP has achieved a relatively high solution efficiency [21], strictly speaking, it is still a nonlinear model, and the efficiency of solving nonlinear second-order cone constraints will decrease with the increase in the number of distribution network nodes.In order to achieve efficient solving of the OPF problem in future large-scale distribution networks, it is necessary to simplify the second-order cone programming constraints while ensuring the accuracy and efficiency of distribution network solving.Ref. [22] has proven that the second-order cone in dimension N, can be outer-approximated to an arbitrary accuracy by a polyhedral cone in an extended space.Polyhedral approximations are very powerful for solving MISOCP as the benefits of warm-starting for LPs can be utilized throughout the branch-and-bound algorithm [23].This is particularly important for problems that can be formulated as MISOCP, such as the OPF problem in this paper.Meanwhile, Refs.[24,25] extended the constant power load model to a ZIP load model under a rectangular coordinate system and applied it to the solution of the OPF.Meanwhile, most of the existing distribution network models only include power generation units and energy storage devices, and rarely consider reactive power compensation devices and other active management devices at the same time [26].In addition, in the active distribution network planning and operation optimization model, discrete variables will inevitably appear with the active management devices considered, which turns the original problem into a mixed integer linear programming (MILP) [27].With continuous development and maturity, commercial optimization software (GUROBI, CPLEX, MOSEK, etc.) has been widely used in distribution network reconstruction [28], reactive power optimization [29], distribution network planning [11], etc.
It is noteworthy that the above studies were mostly limited to the traditional singleperiod static OPF category, while the actual optimization requires the overall coordination of multi-periods, which is actually a dynamic optimal power flow (DOPF).Additionally, compared to the rectangular coordinate system, the current form in polar coordinates is more common.Therefore, in order to meet the operational requirements of multi-period scenarios, various complex constraints, and fast large-scale solutions for active distribution networks, this paper proposes a linearized second-order cone relaxation dynamic optimal power flow (LSOCR-DOPF) model of the distribution network and explores the linear modeling method of key constraints of active distribution network participating elements (such as on-load-tap-changer (OLTC), static var compensator (SVC), capacitor banks (CB), the energy storage system (ESS), etc.).The nonlinear OPF problem of the distribution network is transformed into a computationally efficient solution utilizing MILP, with a comprehensive explanation of the network's constraints on radiation and connectivity.On this basis, this paper further presents an innovative approximation of the ZIP load model under polar coordinates and validates the effectiveness of the OPF framework through three scenarios: power coordination optimization, network reconfiguration, and ZIP load application.The main contributions of this study can be summarized as follows:

•
Present the LSOCR-DOPF model, which is based on branch power flow analysis, for the active distribution network.

•
Explore the linear modeling method for the constraints of various active management units, including OLTC, SVC, CB, ESS, etc.

•
Validate the LSOCR-DOPF model through simulation experiments in three typical scenarios: power coordination optimization, network reconfiguration, and ZIP load application.
The article is organized as follows: The Section 2 presents the LSOCR-DOPF model for active distribution networks and the design of various active management units based on branch power flow analysis.The Section 3 discusses the results of three simulation experiments and provides a comprehensive analysis and discussion of the findings.Finally, the Section 4 presents the conclusion and future perspectives.In most distribution networks, the steady-state power operation mode is radial and its structure is depicted in Figure 1.For a radial topology network, the node directed graph can be used for equivalent analysis.Furthermore, S ij and S i represent complex power, S ij = P ij + Q ij i and S i = p i + q i i. Branch complex impedance Z ij = r ij + x ij i. Set B represents the set of all nodes in the network.In the traditional OPF, the voltage remains constant.If OLTC is installed, the voltage will change with the OLTC transformation ratio; E represents the collection of all branches in the network.There are N sub substations, N bus nodes, and N Line line branches in the network.
based on branch power flow analysis.The third section discusses the results of thr ulation experiments and provides a comprehensive analysis and discussion of th ings.Finally, the fourth section presents the conclusion and future perspectives.

Basic Structure of Distribution Networks
In most distribution networks, the steady-state power operation mode is rad its structure is depicted in Figure 1.For a radial topology network, the node directed can be used for equivalent analysis.Furthermore, and represent complex = + i and = + i. Branch complex impedance = + i. Set resents the set of all nodes in the network.In the traditional OPF, the voltage r constant.If OLTC is installed, the voltage will change with the OLTC transformatio E represents the collection of all branches in the network.There are sub subs bus nodes, and Line line branches in the network.

Basic OPF Model Based on BFM
Generally, the basic model of optimal power flow based on branch power flow OPF, Figure 1) is expressed as follows [14]: , , , ,

Basic OPF Model Based on BFM
Generally, the basic model of optimal power flow based on branch power flow (BFM-OPF, Figure 1) is expressed as follows [14]: min f (p, q, P, Q, V, I) where, p j and q j represent the active and reactive power injections, respectively, at each node; the branch ij represents the positive direction of flow direction from node i to node j; δ(j) is the collection of branch end nodes with j as the head node, and π(j) is the collection of branch end nodes with j as the end node; P jk and Q jk denote the active and reactive power at the head node of the branch ij, respectively; P ij and Q ij denote the active and reactive current flow in each branch; r ij and x ij correspond to the individual resistance and reactance of each branch; g j and b j are the separate ground conductance and ground susceptance of node j; V j and − V j are the upper and lower limits of node voltage, respectively; I ij and − I ij are the upper and lower limits of branch ij current, respectively.
From Equations ( 1)-( 7), it can be inferred that: (1) The optimization variables for the OPF consist of node injection power (p, q), branch power flow (P, Q), node voltage (V), and branch current (I), with the substation node voltage not considered as an optimization variable.(2) Equation ( 1) represents the objective function, which can be the minimization of network losses and substation node power purchases.(3) The well-known branch flow equation [30] is expressed by Equations ( 2) and (3), while Equation (4) represents the power calculation equation.Equations ( 5) and ( 6) denote the safety constraint equations for branch currents and nodal voltages, respectively.Equation ( 7) introduces node-dependent constraints that are subject to change based on the employed model.

Multi-Period LSOCR-OPF Model
In order to convert power flow constraints into quadratic cone constraints, additional variables of cone optimization need to be set: Substitute Equation (8) into Equations ( 2)-( 6) to change the power flow constraint to the second-order cone constraint [31] as follows: min f (p, q, P, Q, V, I) (9) where, || || 2 represents the L2 norm.With this, the SOCR-OPF (a typical Mixed-Integer Second Order Cone Programming, MISOCP) model is fully modeled.Among them, the second-order cone programming represented by Equation ( 12) can be expanded and simplified as: Equation (15) can be uniformly described as: Ben-Tal and Nemirovski [22] showed that Equation ( 16) could be approximated by a system of linear homogeneous equalities and inequalities in terms of x 1 , x 2 , x 3 , and 2(k + 1) variables a k , b k for k = 0, . . ., K. K is a parameter of the polyhedral-relaxed approximation.The approximate expression of the polyhedron of the three-dimensional SOCR constraint (Equation ( 16)) is: where, K represents the number of facets in a polyhedron.a 0 and b k are auxiliary variables.The polyhedral approximation given by ( 17) can be reduced by using the linear equality constraints a k = a k−1 cos (π/2 k+1 +b k−1 sin (π/2 k+1 to solve for a k (k = 0, . . ., K) in terms of a 0 and a k (k = 0, . . ., K) and then substitute a k out of Equation ( 17).The resulting system will only have linear inequality constraints in terms of x 1 , x 2 , x 3 , a 0 , and the (k + 1) variables b k for k = 0, . . ., K. The error (ε(K)) [22] of the polyhedral approximations of the three-dimensional second-order cone constraint (Equation ( 16)) is: According to Equation (18), when K = 11, the error is about 3 × 10 −7 .Therefore, the MISOCP model for active distribution network reconfiguration is approximately equivalent to the MILP model.At this point, the improved version of SOCR-OPF, namely the LSOCR-OPF model, is fully modeled.
So far, the power flow constraint has been transformed from a nonlinear programming model to a MISOCP model composed of Equations ( 9)-( 14).Then, through the polyhedron approximation method, which is used to linearize Equation ( 16), the MISOCP model can be approximately converted into a MILP model for solution.
The previously described model represents the conventional single-period OPF model.However, since most practical applications involve multi-period optimization, this study converts the single-period (static) power flow model into a multi-period (dynamic) OPF model.For clarity of expression, the following dynamic OPF model, based on LSOCR, can be expressed in vector form: where, t is the period identification, T is the total number of time periods.x t ∈ X t , ∀t represents the constraint relationship in the traditional single-period OPF model, such as upper and lower limit constraints, power flow equation constraints, etc. ∑ t∈T B t x t ≤ c is the second-order cone constraint relationship under each branch at each time.
Equation (19) adds the linear coupling relationship between multi-period periods to the objective function and constraint conditions.Some elements will be described in detail in the next section, such as OLTC, CB, ESS, etc.Hence, the further improved version of LSOCR-OPF, namely the LSOCR-DOPF model, is fully modeled.
In order to facilitate understanding, additional explanation is required for the entire model transformation process (from AC-OPF to LSOCR-DOPF): Equations ( 8) and ( 16) respectively embody the phase angle relaxation and second-order cone relaxation of LSOCR-DOPF, and Figure 2 depicts the schematic diagram of the two-step relaxation process.The non-convex feasible region C AC−OPF of the original AC-OPF problem will be relaxed into a convex second-order cone feasible region C OPF−cr after phase angle relaxation and secondorder cone relaxation.Then, the convex feasible region of the second-order cone is further linearized into the convex feasible region C OPF−linear of the integer programming by the polyhedral approximations.At this time, the optimal power flow problem in the original formulation has already been transformed into a convex optimization problem.Numerous studies, as demonstrated in Refs.[13,14,16,17], have substantiated the strict accuracy of the second-order cone relaxation (SOCR) approach for most distribution network structures, when the objective function is both a convex and strictly increasing function.
in the next section, such as OLTC, CB, ESS, etc.Hence, the further improved version o LSOCR-OPF, namely the LSOCR-DOPF model, is fully modeled.
In order to facilitate understanding, additional explanation is required for the entire model transformation process (from AC-OPF to LSOCR-DOPF): Equations ( 8) and (16 respectively embody the phase angle relaxation and second-order cone relaxation o LSOCR-DOPF, and Figure 2 depicts the schematic diagram of the two-step relaxation pro cess.The non-convex feasible region  of the original AC-OPF problem will be re laxed into a convex second-order cone feasible region  after phase angle relaxa tion and second-order cone relaxation.Then, the convex feasible region of the second order cone is further linearized into the convex feasible region  of the intege programming by the polyhedral approximations.At this time, the optimal power flow problem in the original formulation has already been transformed into a convex optimi zation problem.Numerous studies, as demonstrated in Refs.[13,14,16,17], have substan tiated the strict accuracy of the second-order cone relaxation (SOCR) approach for mos distribution network structures, when the objective function is both a convex and strictly increasing function.

Active Distribution Network Modeling
In this section, the active management components of the active distribution network are considered, including: ① reactive compensation device (SVC and CB); ② active

Active Distribution Network Modeling
In this section, the active management components of the active distribution network are considered, including: 1 reactive compensation device (SVC and CB); 2 active power regulation device (energy storage system and electric vehicle mobile energy storage system); 3 on-load tap changer (OLTC); 4 distributed generation power regulation.Given that the objective function in the OPF is either linear or quadratic, the quadratic form can be effectively addressed by means of piecewise linearization.This paper linearizes the constraints related to the active management equipment.Given that the network reconfiguration is of great significance for active distribution network planning and optimal operation [32], this paper discusses the constraints related to grid reconfiguration, such as radial constraints.In addition, in order to ensure the applicability of the model, this paper further considers the ZIP load.
where, B CB is the set of CB nodes; y CB j,t is the number of groups put into operation and the discrete variable value; Y CB,max j is the upper limit of the number of CB groups connected by node j; Q CB, step j is the compensation power of each group of CB.Considering factors such as equipment life or economy, discrete reactive compensation is mostly limited by the number of adjustments, so it generally includes the total number of operations in multiple periods; N CB,max j is the upper limit of operation times: In addition, for the absolute value constraint in the above equation, add an auxiliary variable δ CB j,t = y CB j,t − y CB j,t−1 that represents the change in CB compensation capacity between adjacent periods, corresponding to: (2) Modeling of continuous reactive power regulation device (SVC).
where, B SVC is the node set containing SVC; Q SVC, min j and Q SVC,max j are the lower limit and upper limit of SVC compensation power, respectively.Considering that, in the process of active distribution network operation, with the increasing penetration of distributed generation (DG) such as photovoltaic power generation, the system power flow may be reversed and overvoltage problems may occur, the lower limit Q of SVC compensation in this paper is:

OLTC model
The OLTC is used to adjust the voltage value at the low-voltage side of the bus node.Therefore, the substation bus node V 0 is further converted to the adjustable variable: where, B OLTC refers to the node set of substation containing OLTC; V Base j,t is the voltage value at the high voltage side of the transformer, which is a constant value; r max j and r min j are the square of the upper and lower limit of the OLTC adjustable transformation ratio; r j,t is the square of the OLTC transformation ratio, defined as the ratio of the secondary side to the primary side, which is actually a discrete value variable, and can be further treated as the following relationship including 0-1 variables: where, r j,s represents the difference between OLTC gear s and the square of gear s − 1, which is the adjacent adjustment increment.σ OLTC j,s,t is a 0-1 identification variable.If it is considered to be constrained by the limit of adjustment times in practice, it can be further constrained as: where, δ is the maximum allowable adjustment times of the OLTC gear at time T.

ESS model
During this part, the modeling of the energy storage system takes into account multiple-period constraints, including restrictions on its charging and discharging status, charging and discharging power, as well as capacity limitations.
(3) Capacity constraints. where As a novel active management technique, electric vehicles can be considered as mobile active power energy storage systems [33].Their basic model is largely similar to that of Energy Storage Systems (ESS).The equivalent injection power at each bus node in the distribution network can be represented as the clustering outcomes of individual electric vehicles.

4.
Distributed generation model Respectively modeling DG with or without reactive power: (1) DG modeling without considering reactive power.
Currently, the modeling form of active management for DG primarily takes into account the possibility of DG allowing power shedding under certain conditions, and assumes that DG is only related to active power output [34], that is: , ∀t, ∀j ∈ B DG (30) where, B DG is the node set containing DG; P DG, PRE j,t is the predicted active power output of node j at time t.
With the maturity of active power regulation as a primary method in the study of DG and the advent of new technology, some DGs can have a certain impact on reactive power in the power grid, including outputting and absorbing reactive power.As a result, active management research focusing on the reactive power of DG has emerged, mainly divided into constant power factor control and variable power factor control.Variable power factor control where, F DG, max j and F DG min j are the upper and lower limits of DG transformation ratio adjustment; F DG j,t represents the ratio of reactive power to active power, and its optimization range can be converted from the power factor control range.At present, most of the literature will limit the reactive power Q DG,min j and Q DG,max j is set as a constant [27].

Radial Constraints
Network reconfiguration is of great significance for the active distribution network with high renewable energy penetration, and radiation and islanding constraints are more important constraints in the process of the optimal operation of the distribution network.Since the BFM model assumes the positive direction of the power flow, that is, P ij > 0 indicates that the current flows from i to j, if P ij < 0, the current direction is opposite.In addition, suppose σ ij (0-1 variable) represents the state of branch ij, if σ ij,t = 0, then the branch ij switch is open, and vice versa.It can be obtained in the same way in the distribution network planning, if σ ij,t = 1, then the candidate branch ij is connected, and vice versa.Based on this method, the LSOCR-DOPF model in this paper needs to add the following constraints: Equation ( 33) imposes a stringent restriction on the branch power flow.It serves as both necessary and sufficient conditions, which guarantee a certain level of load power balance.In a purely load-based grid, these conditions can ensure stability.However, in an active distribution network integrating DERs, there may arise the possibility of ring networks or isolated islands.The introduction of connectivity constraints can effectively eliminate such isolated island operations in the power grid.This can be achieved by setting the injection power of all non-substation nodes in the network to a small normal value ε and incorporating it into the power flow constraint, ensuring connectivity between each node and the substation node.The corresponding auxiliary equation is as follows: So far, Equations ( 33) and ( 34) ensure the radiation and connectivity in the network.At the same time, in the process of power grid reconstruction, the number of switch adjustments is also limited.Similar to the constraint on the number of OLTC adjustments in Section 2.2.1, the switch adjustment limit is modeled as follows: where, δ is the maximum allowable adjustment times of the switch at time T.

ZIP Load Application
It is evident that the traditional approach of only considering constant power loads in demand is no longer adequate in the era of advanced technology and the growing demand for refined simulation.Thus, it is imperative to develop a ZIP load modeling that accounts for voltage static characteristics.As noted in [25], the load can be modeled using the ZIP model, which is comprised of constant power, constant current, and constant impedance models.For ease of expression, the time marker t has been omitted.
where, α A j , α B j , and α c j are the proportion of constant impedance load, constant current load, and constant power load; P load j,0 , Q load j,0 are the power demand under the rated voltage.According to the second-order cone relaxation principle, the model is equivalent to: where, for constant current load, − V j , since the unit value of voltage is near the rated voltage of 1.0, the square term is also near the rated value of 1.0, set By expanding it in Taylor series: It should be noted that: The calculated equivalent error order is 10 −3 , which fully meets the optimization requirements.Based on the above analysis,

Proposed LSOCR-DOPF Optimization Flowchart
The comprehensive optimization flowchart of LSOCR-DOPF is depicted in Figure 3.The comprehensive LSOCR-DOPF model described in this paper is comprised of the preliminary LSOCR-DOPF model (Equation ( 19)) introduced in Section 2.1.3and all of the linearized models outlined in Section 2.2.The LSOCR-DOPF is a classic Mixed Integer Linear Programming (MILP) problem, which can be modeled using the YALMIP toolkit in MATLAB and can be efficiently solved by invoking advanced solvers.The network loss, main network power purchase cost and power abandonment Eqs. ( 10),( 11), (16) Eq. ( 13) Eq. ( 14) Eqs. ( 20 Eq. ( 17

Validation and Discussions
This paper mainly analyzes the application of three scenarios based on LSOCR-DOPF: (1) power coordination optimization; (2) network reconstruction; (3) ZIP load application.The numerical simulation analysis was carried out using an active distribution

Validation and Discussions
This paper mainly analyzes the application of three scenarios based on LSOCR-DOPF: (1) power coordination optimization; (2) network reconstruction; (3) ZIP load application.The numerical simulation analysis was carried out using an active distribution system based on the IEEE33 and IEEE69 bus systems data.The simulation parameters of the basic simulation system are presented in Appendix A. The simulation was performed on a 64-bit Windows 11 operating system with an Intel (R) Core (TM) i7-6700HQ CPU @ 2.64 GHz and 24 GB RAM, using MATLAB R2022b and calling YALMIP and GUROBI 10.0.The dynamic simulation period was set to 24 h.The voltage safety range (p.u.) was set to [0.94, 1.06].

Scenario Description
This scenario (Figure 4) adds elements such as OLTC, ESS, CB, SVC, Wind, and photovoltaic (PV) to the basic IEEE33 network to comprehensively verify the application effect of each active management object in LSOCR-DOPF; assuming the OLTC transformation ratio range is 1-6%.
where,  Loss ,  TR ,  WTG , and  PVG are the network loss price, the main network power purchase price, the wind power abandonment price, and the photovoltaic power abandonment penalty price, respectively;  , PVG and  , WTG are the predicted output of wind power and photovoltaic, respectively.

Analysis of Optimization Results
The multi-period operational results of the power coordination and optimization of each active management unit are obtained based on the original data of the scenario, as depicted in Figures 5-10.Figures 5-9 show the active load demand, the active output of each wind turbine unit, the PV active output, the active output of the main grid, and the node injection power, respectively, as well as the active power of energy storage charging and discharging.Figure 10 presents the total 24-h CB reactive power output and SVC reactive power output.The validity of the similar relaxation model has been demonstrated in Refs.[14,17].To further verify its universality, this scenario incorporates the cost of power purchase from the main network, as well as penalty costs for the abandonment of wind and photovoltaic power.The resulting total operational cost (C OPE ) is as follows: where, c Loss , c TR , c WTG , and c PVG are the network loss price, the main network power purchase price, the wind power abandonment price, and the photovoltaic power abandonment penalty price, respectively; P PVG j,t and P WTG j,t are the predicted output of wind power and photovoltaic, respectively.

Analysis of Optimization Results
The multi-period operational results of the power coordination and optimization of each active management unit are obtained based on the original data of the scenario, as depicted in Figures 5-10.Figures 5-9 show the active load demand, the active output of each wind turbine unit, the PV active output, the active output of the main grid, and the node injection power, respectively, as well as the active power of energy storage charging and discharging.Figure 10 presents the total 24-h CB reactive power output and SVC reactive power output.
each active management unit are obtained based on the original data of the scenario, as depicted in Figures 5-10.Figures 5-9 show the active load demand, the active output of each wind turbine unit, the PV active output, the active output of the main grid, and the node injection power, respectively, as well as the active power of energy storage charging and discharging.Figure 10 presents the total 24-h CB reactive power output and SVC reactive power output.It can be seen from Figures 5-10 that: (1) During the peak load periods (08:00-15:00, 19:00-21:00), the active of the system can be met while the Energy Storage System (ESS) is unable surplus clean energy due to its own charging limitations, resulting in pow ment.The ESS discharges during both the peak load period and the low p renewable energy generation, effectively reducing the peak-valley difference alent load; (2) During periods when the proportion of renewable energy output to tively high (6:00-15:00, 18:00-22:00), reactive power compensation device Compensators (SVC) and Capacitor Banks (CB)) absorb the excess reactive system, avoiding overvoltage.
Hence, the LSOCR-DOPF model proposed in this paper demonstrates s fects on the optimization of both active and reactive power.

Model Validity Analysis
The LSOCR-DOPF is analyzed under three situations: "relaxation accu tion efficiency, and comparison of different optimization cases".
(1) Relaxation accuracy The relaxation accuracy of the objective function, such as network lo proved and verified in the pure load network environment.However, consid ther analysis of the model accuracy after increasing the main network output It can be seen from Figures 5-10 that: (1) During the peak load periods (08:00-15:00, 19:00-21:00), the active load demand of the system can be met while the Energy Storage System (ESS) is unable to absorb the surplus clean energy due to its own charging limitations, resulting in power abandonment.The ESS discharges during both the peak load period and the low peak period of renewable energy generation, effectively reducing the peak-valley difference of the equivalent load; (2) During periods when the proportion of renewable energy output to load is relatively high (6:00-15:00, 18:00-22:00), reactive power compensation devices (Static Var Compensators (SVC) and Capacitor Banks (CB)) absorb the excess reactive power of the system, avoiding overvoltage.
Hence, the LSOCR-DOPF model proposed in this paper demonstrates significant effects on the optimization of both active and reactive power.

Model Validity Analysis
The LSOCR-DOPF is analyzed under three situations: "relaxation accuracy, calculation efficiency, and comparison of different optimization cases". (

1) Relaxation accuracy
The relaxation accuracy of the objective function, such as network loss, has been proved and verified in the pure load network environment.However, considering the further analysis of the model accuracy after increasing the main network output, power abandonment, and load loss penalty costs, it is necessary to define the error index: 11 shows the error scatter diagram of each branch in one day.Obviously, the deviation after relaxation meets the requirements of accurate operation, which is 10 -6 .
ther analysis of the model accuracy after increasing the main network output, power donment, and load loss penalty costs, it is necessary to define the error index: , + , − , , .
Figure 11 shows the error scatter diagram of each branch in one day.Obvious deviation after relaxation meets the requirements of accurate operation, which is 1 (2) Calculation timeliness (2) Calculation timeliness Table 1 presents the solution speed and respective simulation results under various objective functions.Specifically, the objective functions of each simulated operating conditions are as follows: Case 1: minimized network loss; Case 2: minimized power purchase cost of the main grid; Case 3: minimum power loss; Case 4: minimized network loss and main grid power purchase cost; Case 5: minimized network loss, main grid power purchase cost, and power abandonment.Table 1 shows the solution speed and corresponding simulation results under different objective functions.Among them, the objective function Case 1 of each simulated operation condition is the minimum network loss; Case 2: the minimum power purchase cost of the main network; Case 3: minimal power loss; Case 4: network loss and main network power purchase cost are the least; Case 5: the network loss, main network power purchase cost, and power abandonment are the minimum.Based on Table 1, with the increase in the objective function, the calculation time continues to increase, but the calculation speed is still acceptable, the maximum is 200 s, which meets the time requirements for day-ahead dispatching and real-time optimization of the active operation in the active distribution network.
Furthermore, based on Case 1, this part compares the solving efficiency of the MINLP (Interior Point-DOPF, Mixed-Integer Nonlinear Programming) model, the MISOCP (SOCR-DOPF) model, and the MILP (LSOCR-DOPF) model of the DOPF problem, and the results are summarized in Table 2.It can be seen from Table 2 that the accuracy and efficiency of LSOCR and SOCR are higher than those of the traditional interior point method.The distribution network loss optimization result of the MILP model linearized by SOCR is equal to the result of the MISOCP model, but the solution speed has been improved to a certain extent.This fully shows that the optimization efficiency can be improved by the LSOCR, which can reach 25~30% (compared with SOCR).It is shown that the LSOCR method is slightly better than the SOCR method in terms of comprehensive solution efficiency, but its solution timeliness is entirely better than the traditional interior point method.(3) Comparison of different optimization cases According to Tables 1 and 2, it can be seen that: 1 In Case 1, only the optimization of network loss is considered, leading to substantial waste of clean energy power.On the other hand, Case 2 prioritizes the minimization of main grid power purchase costs, which results in increased consumption of clean energy, reducing the main grid power purchase and, in turn, reducing power waste.Conversely, Case 3 focuses on minimizing power loss.In comparison to Case 1, it requires more utilization of clean energy to reduce main grid power purchases, but results in an increase in network loss; 2 In Case 4, the objective is to simultaneously minimize both the network loss and the main network's power purchase.The results show that while the main network's power purchase remains unchanged, the network loss decreases, but the amount of power abandonment increases.The paper assumes that the cost of network loss and the penalty cost of power abandonment are both equal to 5000 $/MWh, which suggests that excessive access to clean energy may increase the network loss of the system; 3 Case 5 aims to minimize the network loss, main network power purchase, and power abandonment.The main network power purchase remains the same, while the power abandonment correspondingly decreases, and the network loss increases, which means that the network loss caused by increasing the access to clean energy is less than the penalty of clean energy power abandonment, so the power abandonment is further reduced. 4To assess the impact of the active management equipment (SVC, CB, ESS) introduced in this paper on the performance of the distribution network, three tests were designed and included: (a) No addition of energy storage and reactive power compensation equipment; (b) Addition of an energy storage device in place of a reactive compensation device; (c) Simultaneous addition of both an energy storage device and a reactive compensation device.It can be seen from Figure 12 that adding the active management unit can improve the uniform distribution of voltage in the distribution network to some extent, which shows some functions of the active management unit.It is worth mentioning that the relatively elevated node voltage depicted in Figure 12 can be attributed to the network loss specified by the objective function.
Electronics 2023, 12, x FOR PEER REVIEW 18 of 27 more utilization of clean energy to reduce main grid power purchases, but results in an increase in network loss; ② In Case 4, the objective is to simultaneously minimize both the network loss and the main network's power purchase.The results show that while the main network's power purchase remains unchanged, the network loss decreases, but the amount of power abandonment increases.The paper assumes that the cost of network loss and the penalty cost of power abandonment are both equal to 5000 $/MWh, which suggests that excessive access to clean energy may increase the network loss of the system; ③ Case 5 aims to minimize the network loss, main network power purchase, and power abandonment.The main network power purchase remains the same, while the power abandonment correspondingly decreases, and the network loss increases, which means that the network loss caused by increasing the access to clean energy is less than the penalty of clean energy power abandonment, so the power abandonment is further reduced.
④ To assess the impact of the active management equipment (SVC, CB, ESS) introduced in this paper on the performance of the distribution network, three tests were designed and included: (a) No addition of energy storage and reactive power compensation equipment; (b) Addition of an energy storage device in place of a reactive compensation device; (c) Simultaneous addition of both an energy storage device and a reactive compensation device.It can be seen from Figure 12 that adding the active management unit can improve the uniform distribution of voltage in the distribution network to some extent, which shows some functions of the active management unit.It is worth mentioning that the relatively elevated node voltage depicted in Figure 12 can be attributed to the network loss specified by the objective function.Based on the simulation results of the scenario in this section, an analysis of three important parameters, calculation time, solution error, and node voltage distribution, is conducted to evaluate the effectiveness of the proposed method in this paper.
Firstly, in terms of power grid stability, as shown in Figure 12c, the voltage distribution of each node in the IEEE33 node network can meet the set range requirement of [0.94,Based on the simulation results of the scenario in this section, an analysis of three important parameters, calculation time, solution error, and node voltage distribution, is conducted to evaluate the effectiveness of the proposed method in this paper.
Firstly, in terms of power grid stability, as shown in Figure 12c, the voltage distribution of each node in the IEEE33 node network can meet the set range requirement of [0.94, 1.06], and the distribution is relatively uniform.
Secondly, in terms of power grid scheduling frequency, as shown in Tables 1 and 2, the six sets of comparison cases (IEEE33 + IEEE69) set in this scenario can run within 3 min (minimum of 8.555 s), which can meet the real-time power grid scheduling frequency requirements.
Thirdly, in terms of solving accuracy, as shown in Figure 11, the magnitude of the error in the simulation results of this scenario is strictly controlled within 10-6, which can ensure the requirement of power grid solving accuracy.

Network Reconfiguration
The distribution network reconfiguration scenario is analyzed based on the standard IEEE 33 systems (Appendix A), and the load value is the original system data.In order to verify the effectiveness of second-order cone relaxation in distribution network reconstruction, this example uses static single-period reconstruction for analysis.In the two examples, node 1 is a substation node with a voltage amplitude of 1.06 p.u..It is assumed that each branch is equipped with a section switch, and the optimization calculation is carried out with the goal of minimizing the network loss.The results are shown in Table 3.This also solves the error scatter diagram of each branch, as shown in Figure 13.The distribution network reconfiguration scenario is analyzed based on the standard IEEE 33 systems (Appendix A), and the load value is the original system data.In order to verify the effectiveness of second-order cone relaxation in distribution network reconstruction, this example uses static single-period reconstruction for analysis.In the two examples, node 1 is a substation node with a voltage amplitude of 1.06 p.u..It is assumed that each branch is equipped with a section switch, and the optimization calculation is carried out with the goal of minimizing the network loss.The results are shown in Table 3.This also solves the error scatter diagram of each branch, as shown in Figure 13.As shown in Figure 14, it is not difficult to find that the error is 10 -orders of magnitude, meeting the reconstruction requirements.It is worth noting that, in this example, no active management equipment is added.For IEEE 33 nodes, all nodes have load values, so only Equation ( 33) can be added to the radiation constraint.For some nodes in the IEEE69 system without load, Equation (34) must be added to ensure that all nodes are connected and operate without islands and radiation.
Based on the simulation results of this section, the three important parameters of calculation time, solution error, and node voltage distribution are evaluated.Firstly, in terms of power grid stability, as shown in Figure 15, the voltage distribution of each node in the IEEE33 node network can meet the set range requirement of [0.94, 1.06], and the voltage distribution in multi-period is relatively uniform.Secondly, regarding power grid scheduling frequency, as shown in Table 3, the two sets of comparative cases (single-period and multi-period) ensure that the operating time is within 3 min (the shortest being 1.355 s), meeting the requirements of real-time power grid scheduling frequency.Thirdly, in terms of solving accuracy, as shown in Figure 14, the error level of the simulation results in this scenario is strictly controlled within 10 −9 , which can well guarantee the accuracy require- As shown in Figure 14, it is not difficult to find that the error is 10 -9 orders of magnitude, meeting the reconstruction requirements.It is worth noting that, in this example, no active management equipment is added.For IEEE 33 nodes, all nodes have load values, so only Equation (33) can be added to the radiation constraint.For some nodes in the IEEE69 system without load, Equation (34) must be added to ensure that all nodes are connected and operate without islands and radiation.

ZIP Load Application
This scenario adopts the PG69 node system test (Figure 16).In order test th voltage characteristics of the load, this part does not consider the active manageme and only analyzes the calculation results after citing the static voltage characteri the load.Three types of loads, namely constant power, constant current, and const pedance, are added to each node.The voltage distribution of each node at each calculated as shown in Figure 17. Figure 18 shows the comparison of load active d considering voltage static characteristics or not.Based on the simulation results of this section, the three important parameters of calculation time, solution error, and node voltage distribution are evaluated.Firstly, in terms of power grid stability, as shown in Figure 15, the voltage distribution of each node in the IEEE33 node network can meet the set range requirement of [0.94, 1.06], and the voltage distribution in multi-period is relatively uniform.Secondly, regarding power grid scheduling frequency, as shown in Table 3, the two sets of comparative cases (single-period and multi-period) ensure that the operating time is within 3 min (the shortest being 1.355 s), meeting the requirements of real-time power grid scheduling frequency.Thirdly, in terms of solving accuracy, as shown in Figure 14, the error level of the simulation results in this scenario is strictly controlled within 10 −9 , which can well guarantee the accuracy requirements of power grid solving.

ZIP Load Application
This scenario adopts the PG69 node system test (Figure 16).In order test the static voltage characteristics of the load, this part does not consider the active management unit, and only analyzes the calculation results after citing the static voltage characteristics of the load.Three types of loads, namely constant power, constant current, and constant impedance, are added to each node.The voltage distribution of each node at each time is calculated as shown in Figure 17. Figure 18 shows the comparison of load active demand considering voltage static characteristics or not.

ZIP Load Application
This scenario adopts the PG69 node system test (Figure 16).In order test the static voltage characteristics of the load, this part does not consider the active management unit, and only analyzes the calculation results after citing the static voltage characteristics of the load.Three types of loads, namely constant power, constant current, and constant impedance, are added to each node.The voltage distribution of each node at each time is calculated as shown in Figure 17. Figure 18 shows the comparison of load active demand considering voltage static characteristics or not.As can be seen from Figures 17 and 18, since most of the node voltage is lower tha 1.0 p.u., the load demand will be reduced after considering the static voltage characteris tics of the load.It is important to consider the static voltage characteristics of active distr bution network fine simulation.In addition, the second-order cone relaxation error of eac period is shown in Figure 19.Obviously, the relaxation effect is also highly satisfactory.As can be seen from Figures 17 and 18, since most of the node volt 1.0 p.u., the load demand will be reduced after considering the static vo tics of the load.It is important to consider the static voltage characteristic    As can be seen from Figures 17 and 18, since most of the node vo 1.0 p.u., the load demand will be reduced after considering the static v tics of the load.It is important to consider the static voltage characterist bution network fine simulation.In addition, the second-order cone relax period is shown in Figure 19.Obviously, the relaxation effect is also hig As can be seen from Figures 17 and 18, since most of the node voltage is lower than 1.0 p.u., the load demand will be reduced after considering the static voltage characteristics of the load.It is important to consider the static voltage characteristics of active distribution network fine simulation.In addition, the second-order cone relaxation error of each period is shown in Figure 19.Obviously, the relaxation effect is also highly satisfactory.Based on the simulation results of this scenario, the calculation ti and node voltage distribution are evaluated.Firstly, in terms of powe shown in Figure 17, the voltage distribution of each node in the IEEE33 meet the requirements of the [0.94, 1.06] range.Secondly, in terms of po frequency, the running time of this case is 5.394 s, which can meet th real-time power grid dispatch frequency.Thirdly, in terms of solution a in Figure 19, the error level of the simulation result in this scenario is within 10 -, which can guarantee the requirement of power grid sol nally, regarding the accurate modeling of distribution networks, as sh the introduction of ZIP loads can change the load demand, demonstrat of the active and detailed modeling of distribution networks.

Conclusions
In this paper, a polyhedral linear approximation method of the relaxation (LSOCR-DOPF) for the dynamic optimal power flow prob branch power flow model of the distribution network is proposed and v linear constraint modeling of active management units is performed, an of the proposed solution method is verified through comparison with three major application scenarios.The results show that the strategy p per can meet the engineering standards for typical application scenar networks: on one hand, the calculation time is controlled within 3 min, requirement of the general power grid scheduling refresh rate and en of OPF problems in future large-scale active distribution networks; on calculation result error is controlled within 10 −6 , effectively avoiding th sequences caused by calculation errors in previous work; moreover, th tion of each node in the distribution network meets the condition con Based on the simulation results of this scenario, the calculation time, solution error, and node voltage distribution are evaluated.Firstly, in terms of power grid stability, as shown in Figure 17, the voltage distribution of each node in the IEEE33 node network can meet the requirements of the [0.94, 1.06] range.Secondly, in terms of power grid dispatch frequency, the running time of this case is 5.394 s, which can meet the requirements of real-time power grid dispatch frequency.Thirdly, in terms of solution accuracy, as shown in Figure 19, the error level of the simulation result in this scenario is strictly controlled within 10 -8 , which can guarantee the requirement of power grid solution accuracy.Finally, regarding the accurate modeling of distribution networks, as shown in Figure 18, the introduction of ZIP loads can change the load demand, demonstrating the importance of the active and detailed modeling of distribution networks.

Conclusions
In this paper, a polyhedral linear approximation method of the second-order cone relaxation (LSOCR-DOPF) for the dynamic optimal power flow problem based on the branch power flow model of the distribution network is proposed and validated.Detailed linear constraint modeling of active management units is performed, and the effectiveness of the proposed solution method is verified through comparison with other solutions in three major application scenarios.The results show that the strategy proposed in this paper can meet the engineering standards for typical application scenarios of distribution networks: on one hand, the calculation time is controlled within 3 min, which satisfies the requirement of the general power grid scheduling refresh rate and ensures fast solution of OPF problems in future large-scale active distribution networks; on the other hand, the calculation result error is controlled within 10 −6 , effectively avoiding the catastrophic consequences caused by calculation errors in previous work; moreover, the voltage distribution of each node in the distribution network meets the condition constraints [0.94, 1.06] set by authors, which can achieve the optimization goal of reducing network loss.It is noteworthy that the LSOCR-DOPF model satisfies the requirements of calculation time and accuracy in the fields of daily scheduling, real-time operation, and the control of active distribution networks, demonstrating strong practical application value, and the case study suggests that considering the ZIP model for the fine simulation of active distribution networks is also of great significance.Currently, the inclusion of discrete variables, which are non-convex sources, has an impact on the accuracy of the relaxation model, and the universality condi-tion of the second-order cone relaxation model requires further theoretical investigation in future work.
Furthermore, it is difficult to ensure the accuracy of convex relaxation using a relatively universal method due to the various boundary conditions that need to be set in the model for practical applications.The technical challenges of applying convex relaxation techniques in engineering can be summarized as follows: further exploring the influence of the objective function and feasible region of the OPF problem on convex relaxation and seeking sufficient conditions to guarantee accurate convex relaxation in theory; constructing tighter and more precise relaxations based on SOCP relaxation techniques and exploring the possibility of combining convex relaxation techniques with other OPF solution methods.
LSOCR-DOPF Model of the Distribution Network 2.1.1.Basic Structure of Distribution Networks

Figure 1 .
Figure 1.Structure of radial distribution network.

Figure 2 .
Figure 2. The proposed modeling process of LSOCR for solving AC-OPF.

Figure 8 .
Figure 8. Active power output of main power grid and active power diagram of node injection in 24-h.

Figure 8 .
Figure 8. Active power output of main power grid and active power diagram of node injection in 24-h.

Figure 8 .
Figure 8. Active power output of main power grid and active power diagram of node injection in 24-h.

Figure 8 .
Figure 8. Active power output of main power grid and active power diagram of node injection in 24-h.

Figure 8 .
Figure 8. Active power output of main power grid and active power diagram of n 24-h.

( a )Figure 9 .
Figure 9. SOC and charging and discharging of ESS.

6 Figure 11 .
Figure 11.Scatter chart of error under each branch in each period.

Figure 11 .
Figure 11.Scatter chart of error under each branch in each period.

Figure 13 .
Figure 13.Optimized and reconstructed IEEE33 network (node serial number is given by black numbers, branch serial number is shown by red numbers, and red cross and red dotted line indicate disconnected branch).

Figure 14 .
Figure 14.Scatter diagram of distribution network reconfiguration error.

Figure 14 .
Figure 14.Scatter diagram of distribution network reconfiguration error.

Figure 17 .
Figure 17.Voltage distribution diagram of each node.

Figure 18 .
Figure 18.Comparison of load active demand considering voltage static characteristics or not.

Figure 17 .
Figure 17.Voltage distribution diagram of each node.

Figure 18 .
Figure 18.Comparison of load active demand considering voltage static charact

Figure 17 .
Figure 17.Voltage distribution diagram of each node.

Figure 17 .
Figure 17.Voltage distribution diagram of each node.

Figure 18 .
Figure 18.Comparison of load active demand considering voltage static chara

Figure 18 .
Figure 18.Comparison of load active demand considering voltage static characteristics or not.

8 Figure 19 .
Figure 19.Error scatter chart of ZIP load model.

Figure 19 .
Figure 19.Error scatter chart of ZIP load model.

Table 1 .
Optimization results under four types of objective functions.

Table 2 .
Comparison of optimization results of DOPF.

Table 3 .
Reconfiguration scheme and optimization results.

Table 3 .
Reconfiguration scheme and optimization results.

Table A2 .
System data for 33-bus radial distribution network.