1. Introduction
There are numerous variables of uncertainty in power systems [
1,
2], which are typically handled as planning and forecasting problems. Power systems involve more complex analyses, such as load flow [
3,
4,
5], and using deterministic methods for analysis cannot provide accurate results. Borkowska [
6] introduced the concept of probabilistic flow, which is applied to the DC model. The injection power of uncertain nodes is modeled using probability functions. Subsequently, the probability functions are applied to the AC model, initiating numerical analysis [
7]. When analyzing probabilistic flow, boundary flow analysis [
8] is conducted to evaluate the maximum and minimum values that flow state variables can assume.
However, there are still many uncertain variables in power systems that can be considered as non-probabilistic in certain operating conditions. The magnitude of loads is constantly changing at any given moment due to economic growth, changes in social conditions, and human behavior. Therefore, in some cases, using past experience to evaluate future events is an incomplete approach. On the contrary, uncertainty in actual power systems is often associated with imprecise qualitative information and is also linked to certain common expressions, such as “rated power of the motor is between 15 and 25 MW” or “negative power is approximately 5 MW”, which are typically derived from the experience of system operators. Uncertainty analysis of existing models is necessary when non-probabilistic uncertainties arise in the power system [
9]. Fuzzy sets are commonly used for analysis as they can represent and analyze variables with fuzzy and imprecise mathematical precision. Miranda and Mato proposed fuzzy power flow calculation [
10,
11,
12,
13] to study fuzzy variables in power system flows, aiming to compute the uncertainties in power flows by utilizing fuzzy membership functions for node injection powers and addressing potential distribution issues. The possibility analysis model introduced in power flow calculations has been widely accepted. In [
14], the fuzzy interval algorithm [
15,
16,
17,
18,
19] is used to solve power flow calculation problems. Although it is a fuzzy model in specific cases, it represents an alternative way to handle uncertain boundaries. It paves the way for new probabilistic power flow solutions, where this method requires imprecise and fuzzy information as input variables instead of probability density functions.
Fuzzy power flow has been widely applied in the operation and planning of power systems since its proposal, including power transmission systems [
20,
21], distribution networks, and distributed generation systems [
22,
23,
24,
25]. This has led to new requirements for fuzzy power flow algorithms. Consequently, over time, traditional fuzzy power flow algorithms have undergone changes. In [
26], improvements to the power flow calculation method were made based on numerical values such as branch losses, voltage, and current between nodes. Although there are many traditional applications, their use in large-scale systems is limited. This paper aims to improve the traditional fuzzy power flow algorithm to make it suitable for large-scale power transmission systems. Research on power flow calculation has a long history, and there is already a certain foundation in studying the convergence of power flow calculation, especially the convergence issue of the Newton method for power flow calculation, which has been favored by researchers. In order to address the sensitivity of the Newton method to the selection of initial values, in [
27], the authors propose using the PO decoupling method for 1–2 iterations and using the updated solution as the initial value to start the Newton iteration, which has achieved good results. The power flow module in the PSD-BPA (2.9) software of the Electric Power Research Institute still retains this option. In [
28], the authors replace the node admittance with the injected power as the elements of the Jacobian matrix and improve the convergence of power flow when calculating nodes with large loads by enlarging the smallest eigenvalue of the Jacobian matrix. Another group of algorithms takes a different approach. In [
29], the elimination tree theory for power flow calculation is proposed. This method does not involve inverting the admittance matrix, thus avoiding convergence issues. However, current power system networks are generally complex with a large number of nodes, which increases the computational burden and slows down computer performance when applying this method to complex power networks, making it difficult to improve convergence. In [
30], a power flow calculation method based on adjoint theory is proposed, which uses numerical integration to solve the power flow equations and constructs a least squares formulation. This method has achieved certain results in challenging-to-converge small systems.
The Levenberg–Marquardt (LM) method is an iterative algorithm for solving non-linear least squares problems. It combines the features of the most rapid descent method with the Gauss–Newton method and aims to find the parameter that minimizes the sum of squares of the residuals between the predicted and observed values of the model. In [
31], they first applied the LM method to power flow calculations, and through numerical simulations, it was demonstrated that the LM method expanded the range of initial value selection for power flow convergence. In the Levenberg–Marquardt method, it is first necessary to define an objective function and transform the fitting problem into a problem of minimizing the objective function. The algorithm then progressively updates the parameters by iterating to minimize the objective function. During each iteration, the gradient of the objective function under the current parameters, i.e., the partial derivative of the loss function with respect to each parameter, is first calculated. Then, based on the current gradient and previous iteration information, a new parameter update direction is calculated. Then, based on this update direction, a new parameter value is computed and used to update the model. The main feature of the Levenberg–Marquardt method is that in each iteration, the Hessian matrix of the objective function is multiplied by a positive definite matrix to obtain a more accurate parameter update direction. The parameters of this positive definite matrix are adjusted in each iteration to control the step size and direction of the model update.
The above references analyze the improvement of convergence for traditional Newton power flow calculation methods in conventional power grids. This paper will be based on the Levenberg–Marquardt method, focusing on the energy interconnection system, to study methods to improve the convergence of fuzzy Newton method power flow algorithms.
In this paper, by analyzing the convergence theorem and convergence tendency of the fuzzy Newton tidal current algorithm and combining it with existing methods to improve the convergence of the tidal current calculation, we propose an analytical method to improve the convergence of the fuzzy Newton tidal current calculation based on the Levenberg–Marquardt method, which reduces the number of iterations and the running time of the tidal current calculation and expands the convergence region. This method reduces the number of iterations and the running time of the tidal current calculation and enlarges the convergence region. In conclusion, the analytical method to improve the convergence of the fuzzy Newton method based on the Levenberg–Marquardt algorithm proposed in this paper has important theoretical significance and practical value. It provides a new and effective tidal current calculation method in the field of power system analysis, which helps to improve the stability and operation efficiency of the power system.
The rest of the paper is organized as follows: In
Section 2, the fuzzy set theory is applied to construct a fuzzy tidal current model for the generating units of the power system. Subsequently, in
Section 3, the convergence problem of fuzzy current calculation is investigated, and it is shown that the Jacobian matrix is an important factor affecting the convergence of the current. And we analyze the convergence of the fuzzy tidal current calculation method and propose corresponding methods to improve the convergence. In
Section 4, we use the Levenberg–Marquardt method applied to the fuzzy tidal current calculation. Given the content of the previous analysis, it is known that the influence of the Jacobian matrix on the convergence of the tidal current calculation is huge. In order to improve the convergence of the tidal currents, we reconstruct the Jacobian matrix based on this method. Finally, the simulation results show that our method reduces the number of calculation iterations and expands the convergence region, which verifies that our proposed method is effective.
2. Fuzzy Modeling of the Uncertainty in the Output and Availability of Power Generation Device
In this paper, the fuzzy modeling of the power generation unit is achieved using the trapezoidal fuzzy number. The estimates of the power generation unit appear between
, but most likely between
, and the uncertainty of the power generation unit is described through the affiliation function by using the following trapezoidal fuzzy number
.
The fuzzy number center value is The mean of the intercept is .
The fuzzy availability of the kth generator at the ith node is denoted by
, the active fuzzy output and fuzzy load are denoted by
,
, and the reactive fuzzy output and fuzzy load are denoted by an
, respectively. In this case, the fuzzy active and fuzzy reactive powers injected by the node
i can be represented by the following equation:
In Equation (2), is the number of power-generating devices on node . By performing fuzzy number operations on , , it can be determined that they are trapezoidal fuzzy numbers.
Distributed generation device wind turbine generators are driven by the wind, and the wind blades drive the horizontal or longitudinal axis to rotate. This converts mechanical torque into mechanical energy, and the generator uses the transmission device to obtain mechanical power. The stator induces an electric current via the rotation of the rotor, and the power generated is subsequently transferred to the grid on the stator side.
The wind speed, which conforms to the Weibull distribution, is given by the following equation:
where
denotes the shape parameter,
denotes the scale parameter, and
is the random variable of the Weibull distribution.
In fuzzy tidal current calculation, the wind speed needs to be fuzzified to derive the fuzzy power emitted by the wind turbine. In order to obtain the fuzzy distribution function of wind power, a probability distribution function is used, which is transformed into a fuzzy distributed function by fuzzifying the wind speed. That is to say, the left and right sides of the maximum value
are fuzzified, as are the left side of any value
and the right side of any value
:
where
is the fuzzy affiliation function.
Nodes with loads in energy interconnection systems usually have transformers connected to them. Generally, the size of the load is related to the capacity of the connected transformer, since the size of the load of the node will not be lower than 20% or higher than 120% of the capacity of the transformer, but normally varies from 60% to 80% of the capacity of the transformer. So, the magnitude of the node load is described by a set of fuzzy interval numbers corresponding to the determination intervals and . , where where the rated values of active power P and reactive power Q (denoted by x) at node of ab are denoted by .
The affiliation function is shown below:
where
stands for P or Q.
For the electric vehicle charging pile load in the energy interconnection system, assuming that it can run continuously from entering the operation state to exiting, the number of fuzzy intervals that determines the boundary [0, ] or with uncertainty can describe the operation process, where denotes the maximum capacity of the charger.
For the industrial loads of an electric arc furnace and aluminum electrolysis in the energy interconnection system, its fuzzy interval model can be obtained using the same method, where denotes the maximal capacity of an electric arc furnace and aluminum electrolysis. When the variables of the energy interconnection system are in the fuzzy interval , , , and the subordination value of the variables is 1, the fuzzy trend calculation is changed to the deterministic trend calculation; when , the fuzzy trend calculation is changed to the triangular fuzzy number; when , it becomes the triangular fuzzy number. When , it becomes the traditional Newtonian tidal current calculation.
3. Convergence Analysis of the Fuzzy Newton Method
Since the convergence of Newton’s method is mainly related to the Jacobian matrix, to analyze the convergence of the fuzzy Newton’s method, the Jacobian matrix has to satisfy the following conditions.
Assumption 1. The Jacobian matrix satisfies the Lipschitz condition after the kth Newton’s method of selecting generations:where is a real number. Lemma 1. Let .
is a non-singular matrix; ,
and , then Theorem 1. The initial value of voltage and the Jacobian matrix in the fuzzy trend calculation satisfy the following conditions:where , are real numbers. In the neighborhood of the voltage magnitude , if the fuzzy tidal current calculation convergence operator satisfies the following conditions: The solution sequence of the fuzzy Newton method for power flow calculation converges at , indicating the existence of a solution within the range of .
Proof of Theorem 1. From the quadratic merit function, we get:
The two analytical solutions can be obtained as
,
, which is carried into the Equation (10) and the Equation (11) to get:
The iterative sequence , , is denoted by , which confirms that . Next, we analyze .
Let
,
,
, when
,
,
. Derivation of Equation (11), we get:
we substitute the
and
into Equations (13) and (14), and we get:
Then, we construct a function of the following form:
From Equation (15), we have:
Substituting the variables
,
into Equations (16) and (17) yields:
Noting that
,
, there are:
From (6), (16), (21) and
we find that:
Due to
and
. Based on Lemma 1, it follows that
exists and that
Since the magnitude of the voltage amplitude is changing as the iteration proceeds, it is important to note what
represents, as shown by Equations (10), (11), and (21).
The amplitude of the voltage is expressed as .
Similarly to Equations (25) and (26), it follows identically that
,
k = 0, 1, 2 … and we get:
For any positive integers
,
, we have:
It follows from the convergence of
tk that {
V(k)} is a Cauchy convergent sequence, and thus {
V(k)} has a limit, which is set to be
V∗, so that
in the above equation leads to
Observing the convergence of
and through Equation (27), it can be seen that , the solution of the fuzzy tidal current computation is , which can end the whole proof of inference process. □
Using this theorem, the difficulty of choosing the initial value is reduced, and fewer computational steps are required to calculate the trend operator ρ before the trend calculation to determine whether the initial value can find a convergent solution for the trend calculation. If ρ is less than 1/2, the next trend computation can continue; otherwise, it is necessary to repeatedly train other initial values.
4. Levenberg–Marquardt Method Applied to Fuzzy Tidal Current Calculation
4.1. Convergence Improvement Proof
In this paper, the Levenberg–Marquardt method is used to improve the convergence of the fuzzy Newton method of tidal current computation, which belongs to one of the non-linear least squares computational methods. A first-order
Taylor expansion of the basic fuzzy tidal current equation
is performed at the iteration point
x(
k):
where
denotes the iterative correction and
,
, the least squares of the trending equations are notated:
We substitute the Equation (30) into Equation (31). Taking a step size of 1/2 for constraints, the computational model for the Levenberg–Marquardt method can be obtained as:
where
δ is a dynamic variable that can be changed according to the demand; to reduce the sum of squares of iterative errors, the magnitude of
has to be determined first. Using the Levenberg–Marquardt method to limit the magnitude of iterative corrections by adding a damping factor
μ as an additional constraint
can be expressed as:
The introduction of the damping factor μ effectively alleviates the divergence problem in the convergence process, and at the same time gives the Levenberg–Marquardt method the following features:
- (1)
When λ ≥ 0, the number of second-paradigm conditions of the matrix decreases monotonically as the damping factor μ increases.
Proof. Assuming that the
matrix is an
order matrix with eigenvalues
then for any x,
, so the matrix
is a semi-positive definite matrix, and thus all of its eigenvalues can satisfy
. Therefore,
. This is because the singular values of a positive definite matrix and the absolute values of its eigenvalues are equal. Since the condition number of the matrix
is satisfied:
Thus, the second-fan condition number of matrix A decreases monotonically as μ increases. □
- (2)
The iterative correction for the Levenberg–Marquardt method decreases monotonically with the increase in damping factor μ.
Proof . Due to
. The result obtained by taking the derivation of Equation (27) with respect to
can be obtained:
where
is a positive–definite matrix, then
, and
decreases monotonically as the damping factor
μ increases. □
4.2. Convergence Improvement—Reconstructing the Jacobian Matrix
Based on the Levenberg–Marquardt method of improving the fuzzy Newton method of tidal current calculation, the iterative corrections are expressed in the form of (34), which is rewritten in the form of an equation as follows:
From the above equation, the matrix
is rearranged to obtain the same sparse matrix as the conductor matrix, and the sparsification process is carried out according to the following steps:
where
y is the intermediate variable,
.
Substituting Equation (39) into Equation (38) yields:
Joining Equations (39) and (40) yields the matrix equation shown below:
where
, and the coefficient matrix of the matrix equation is denoted by
A.
The coefficient matrix A of Equation (12) is sparsely reordered, and the sparse coefficient matrix
A is as follows:
where the chunking matrix has the same structure as the conductivity matrix.
Figure 1 shows the coefficient matrix sparse rearrangement undirected graph. The corresponding undirected graph of the system’s conductivity matrix is
M1,
M2 has the same structure as
M1, and the associated nodes of
M2 and
M1 are connected to obtain the topology graph denoted as
M3, which is described by
M3, which also describes the coefficient matrix
A.
M2 has the same structure as M1, and the nodes of M3 can be regarded as the generalized nodes of the associated nodes M1 and M2. Therefore, the structure of M3 is the same as that of M1, so the node ordering of M1 applies to M3 as well.
5. Simulation Studies
In this paper, an algorithm for improving the convergence of the fuzzy Newton method of tidal current computation based on the Levenberg–Marquardt method is validated by means of the IEEE 118 standard node system.
Figure 2 illustrates the IEEE 118 node structure. The generating unit nodes in this system are 1, 6, 10, 25, 26, 38, 69, 72, 89, 91, 105, 107, 112, 116, etc., and load unit nodes are 16, 45, 47, 58, 70, 85, etc. Where
type nodes are 72, 91;
type nodes are 25, 116;
type nodes are 38;
type nodes are 1, 16, 47, 58;
type nodes are 10, 26, 89, 105;
type nodes are 70, 85;
type nodes are 69; and
type nodes are 6, 45, 107, 112.
This paper analyzes the convergence of the tidal current calculation based on the initial values of voltage and phase angle of the nodes and the convergence factor
K given in the following formula. When
K = 1, it means that the convergence has been achieved and there is no need for additional iterations, and when
K = 0, it is equivalent to the initialization of all the unknown voltages to 1. In addition, it should be noted that the fuzzy tidal current calculation cannot be converged normally outside the range of a certain convergence factor
K.
There are many methods in the literature for determining the damping factor
μ. The Levenberg–Marquardt method for determining the damping factor has a wide range of usability, and it is generally accepted that when the initial value is poorly chosen, it is necessary to choose a larger damping factor
μ. In this paper, the damping factor chosen for the simulation varies with the number of iterative corrections, which can be obtained as:
where
,
.
Figure 3 and
Figure 4 depict three different algorithms, namely, the fuzzy Newton method (FNR), the fuzzy Newton method based on the Levenberg–Marquardt method to improve convergence of tidal currents (FNRM), and the fuzzy Newton method based on the Levenberg–Marquardt method to improve convergence of tidal currents after Jacobian matrix reconstruction (FNRMR), with the change in the number of iterations and change in running time graphs, respectively. Algorithms, the number of iterations, and the running time are plotted as the convergence factor changes.
From the comparative analysis of the above figures, it can be seen that the FNRM and FNRMR methods contain the convergence situation when all unknown voltages are initialized to 1 when = 0. Moreover, the FNRM and FNRMR methods do not need longer iteration times and running times compared with the traditional FNR method, but the convergence region is three times that of the traditional FNR method, and the FNRMR method is able to converge very well even in the case of poorly selected initial values. In addition, the FNRMR method reconstructs the Jacobian matrix through the sparse matrix. Although the process of the sparse matrix increases the running time, the FNRMR significantly reduces the number of iterations, which leads to the conclusion that the fuzzy Newton tidal current algorithm based on the Levenberg–Marquardt algorithm can solve the problem of the deviation of the initial value very well, and expand the convergence region, which improves the tidal current. For the simulation results obtained above, we summarize our advantages from the following two aspects: (1) Number of iterations: Since the fuzzy Newton method with adjusted Jacobian matrix is characterized by dynamically adjusted convergence factors, it can usually approximate the optimal solution faster, and therefore fewer iterations will be required. Meanwhile, the fuzzy Newton method may require more iterations to converge. (2) Running Time: Since the fuzzy Newton method, after adjusting the Jacobian matrix, requires additional computations to update the convergence factors, it may require more running time. However, if the method is able to converge faster, the total running time may be shorter than the fuzzy Newton method. The Levenberg–Marquardt fuzzy Newton tidal current algorithm can solve the problem of initial value deviation well and expand the convergence region to improve the convergence of the tidal current calculation.
6. Conclusions
By analyzing the convergence theorem and convergence tendency of the fuzzy Newton tidal current algorithm and combining them with the existing main methods to improve the convergence of tidal current calculation, an analytical method to improve the convergence of fuzzy Newton tidal current calculation based on the Levenberg–Marquardt method is proposed and verified on the IEEE standard 118 points. The convergence of fuzzy Newton tidal current calculation, fuzzy Newton tidal current calculation based on Levenberg–Marquardt method, and fuzzy Newton tidal current calculation based on Jacobian matrix reconstruction after the Levenberg–Marquardt method are analyzed by comparison with the simulation. The number of iterations and running time of the fuzzy Newton method, the fuzzy Newton method based on the Levenberg–Marquardt method, and the fuzzy Newton method based on the Levenberg–Marquardt method after the Jacobian matrix reconstruction are compared and analyzed by simulation, and the improved method reduces the number of iterations and enlarges the area of convergence, which verifies that the method can improve the convergence of the fuzzy Newton method.