1. Introduction
The modern power network has changed in many aspects; the power reliability, stability and control has improved dramatically. However, the supervision and energy management systems have became more demanding and complicated [
1]. The power grid is continuously expanding with the deployment of large-scale low-carbon technologies, super-grids and continental transmission networks [
1,
2,
3,
4,
5]. As a result, a definition of hybrid Alternating Current (AC)/Direct Current (DC) networks has become realistic with cost-effective and high power capability advantages. More specifically, the intensive integration of HVDC grids has provided a controllable network with an overall better performance over the conventional high voltage networks [
1,
3,
6]. Besides, it has allowed asynchronous power networks to interconnect, and to achieve higher power transfer capability [
2,
6,
7]. As a result, the demand for an hybrid HVDC/AC network SCADA/Energy Management Systems (EMS) became essential with new and high priority functions for controlling HVDC alongside AC [
1,
8,
9].
Back in 1970s, a centralized (monolithic) SCADA with remote data equipment started to be used to control the different power network sectors (generation, transmission and distribution) [
2,
10,
11]. However, this has changed when each section became larger, and dedicated but connected SCADAs were provided for each energy level, later called hierarchical control and supervision (stations and sub-stations) [
2,
12]. The same architecture is used with the HVDC stations and transmission grids, they are operated by separate SCADA systems as shown in
Figure 1 [
1,
2,
3,
13,
14,
15]. As a result, the research on SCADA for hybrid HVDC/AC transmission (mainly VSC based) is divided into two main approaches, a distributed or decentralized approach such as in [
16,
17,
18,
19], and centralized or hierarchical structure such as in [
1,
2,
20,
21,
22]. In the centralized approach, both AC and DC sides are controlled by one unified SCADA, leading to several challenges and modifications have to be addressed such as:
New timescales requirements in data transmission [
23];
Advanced Remote Terminal Unit (RTU) and Intelligent Electronic Devices (IED) characteristics (e.g., hybrid HVDC/AC sensors/transducers) [
23];
Unified Human Machine Interface (HMI) and power system applications for AC and DC network (i.e., modified toolboxes such as the state estimator and fault detection blocks) [
1];
Cyber-security and communication overload traffic [
1].
A state estimator toolbox is part of the power system applications available in SCADA/EMS; it is used for better observation of the network states by processing noisy and redundant measurements [
24,
25,
26]. The tool requires accurate modeling of the network, including the network circuit breakers states and settings [
9].
A hybrid HVDC/AC system state estimator requires modification over the traditional AC method. First, the DC side components and measurements have to be included. Second, the coupling between the AC and DC sides has to be modeled (hence, the converter impact). Third, a robust estimation has to be guaranteed using an AC and DC bad data detection block [
19,
27,
28] or using an estimation algorithm claimed to be robust (e.g., Least Absolute Value (LAV) [
29,
30]).
In this paper, formulation and implementation of a state estimator are provided based on WLS for hybrid VSC-HVDC/AC transmission networks. The additional modifications on the AC and DC sides coupling are derived and explained. The converters at the Point of Common Coupling (PCC) are modeled as VSC based. The algorithm estimates the AC and DC system states and the converter’s sides coupling in a unified form.
The structure of the paper is as follows:
Section 2 goes through a literature review on VSC converters and the state estimation toolbox.
Section 3 reviews the WLS state estimator and the mathematical formulations.
Section 4 introduces the unified multi-system hybrid VSC-HVDC/AC state estimator.
Section 5 presents the results and the validation of the modified method.
Section 6 concludes the work of this paper.
4. The Unified WLS for HV-(AC/VSC-DC/AC) Transmission Network
Considering the presented AC and DC side coupling approaches, a unified state estimator can be constructed. It is a single state estimator that solves one algorithm for the whole VSC-HVDC/AC system together, using a single unified Jacobian matrix. The modified measurements input for the WLS is formulated below:
where
are the AC side measurements and can be in the form of ;
are the DC side measurements and can be in the form of ;
are zero measurements and represent the right side of converter power coupling constraints;
are the to ratios (measurements), calculated and transmitted by the converter;
n, m and k are the number of AC systems, DC systems and converters respectively.
The measurement functions in the unified WLS have the same format of the z vector.
The new Jacobian matrix is formatted in the following manner:
Which can be reformatted differently to simplify the structure as shown below:
The expanded unified Jacobian matrix in (
29) represents all hybrid VSC-HVDC/AC system components in one matrix. It includes the AC, DC and converter power and voltage coupling elements, and it is structured for
AC systems,
DC systems and
Converters.
where
is the partial derivative of to and
is the partial derivative of to and ↔ 0-matrix
is the partial derivative of to ↔ 0-matrix
is the partial derivative of to
is the partial derivative of the power coupling constraint to and
is the partial derivative of the power coupling constraint to
is the partial derivative of to and
is the partial derivative of to
The following subsection presents the partial derivatives of the additional components to the unified Jacobian matrix
. The traditional AC and DC derivative components are available in
Appendix A.1 and
Appendix A.2 respectively.
4.1. Converter Components: Power Coupling
presents the partial derivative of the power coupling constraint at both sides of the converter. For the DC side: Taking the partial derivatives to the DC voltage magnitude as shown in Equation (
30).
Taking the partial derivatives to the phase angles as shown in Equation (
32), which can be fully expressed in Equation (
33).
Taking the partial derivatives to the AC voltage magnitude as shown in Equation (
34), which can be fully expressed in Equation (
35).
4.2. Converter Components: Voltage Coupling
presents the partial derivative of the voltage coupling at both sides of the converter. For the DC side: Taking the partial derivatives to the DC voltage magnitude as shown in Equation (
36).
Taking the partial derivatives to the AC voltage magnitude, it can be expressed as shown in Equation (
38).
Since the
constraint is angle free, then taking the partial derivatives to the phase angles is zero as shown in Equation (
39).
It is worthy to note the can be included in the system states instead of system measurements, which makes the system states vector contains: AC phase angles and voltages, DC voltages and .
5. Decentralized vs. Power Coupling vs. Unified: State Estimation Simulations
In this section, a validation of the unified WLS is provided along with a comparison with the decentralized approach. The unified WLS uses the power and voltage coupling of the converter to connect the decentralized AC and DC grids.
The unified approach of the WLS state estimator was implemented on Julia optimization programming language [
55]. It was structured in a flexible way that allows easy switching between the different test cases and scenarios (Decentralized, converter power/voltage coupling and unified). The WLS algorithm was implemented in recursion loops with dictionaries and data structures for better memory optimization and arithmetical processing. The true measurements of the hybrid VSC-HVDC/AC systems were obtained by the power-flow solver from PowerModelsACDC.jl (Julia library) [
54]. The noisy measurements are calculated by adding Gaussian noise to the true measurements based on Equation (
40).
The Gaussian noise
is generated with
equals to the true measurement and
is calculated from Equation (
41) [
56]. Further details on the equation derivation are available in [
25].
where
is assumed to cover
of the Gaussian distribution curve (
), and
is the additional error percentage. In this work, all noisy measurements were corrupted with 3% percentage error, except
measurements with 1%.
The simulation study compares between the different state estimation scenarios as follow:
Decentralized: The systems are assumed to be separated with no data/measurements exchange (no coupling). The Jacobian matrix is reduced to contain only AC and DC components, and multi-thread WLSs are run for each AC and DC system.
Converter P-Coupling: A single WLS is run centralized with the converter power coupling constraints added to the Jacobian matrix.
Converter PV-Coupling (unified): The power constraints and voltage coupling measurements are taken into the centralized WLS estimation, forming the unified estimation approach.
The WLS algorithm is configured to have a maximum number of iteration equal to 20 and the tolerance stop condition, related to the change of correction (
), equal to
. The weights of the measurements are are shown in
Table 1.
5.1. Hybrid VSC-HVDC/AC Networks Test Cases
Two hybrid VSC-HVDC/AC networks were modeled to test the algorithm (
Section 5.1.1 and
Section 5.1.2). Both networks are assumed to have bipolar DC links (hence
). The following subsections present each network WLS simulation results for three scenarios: Decentralized AC and DC systems, AC/DC power coupling and unified.
5.1.1. Four(4)-AC/Four(4)-DC/Four(4)-AC Network
The first network is comprised of two 4-bus AC systems, numbered from 1 to 4 and 5 to 8, respectively, connected by a DC grid, as shown in
Figure 7. The AC systems have two AC generators at bus 1 and 6; both considered slack buses (references). The AC load is presented only at bus 8. The DC grid has 4 buses, with bus 1 as a slack bus. Bus 3 has a virtual converter is used to force the power flow direction.
In
Table 2 and
Table 3, the true, noisy and estimated measurements of the DC and AC systems are presented respectively, where T1 is: Decentralized, T2: With converter P-coupling, and T3: Unified.
Figure 8 visualizes the accuracy performance for the three scenarios in some of the AC and DC power injection measurements. It shows the error between the true measurements against the noisy and estimated measurements. The bars represent the relative error in percentages, lower means closer to the true measurement. In
Figure 8a, the estimation of the DC power injection at bus 3 has shown high relative error in the three estimation approaches because this measurement represents a virtual converter.
The converter power coupling measurement is presented as a power constraint equal to zero (Equation (
20)), and it aims to correct the power measurements related to the converter. While the voltage coupling measurement is defined as a ratio (
) that relates the voltage magnitudes of both sides of the converter. The unified WLS was able to accurately estimates the
measurements and the power coupling constraints as shown in
Table 4.
Figure 9 shows the Mean Absolute Error (MAE) of the converter measurements. In (a) the converter power coupling has minimized the error in the DC grid measurements compared to the decentralized method, it has reduced from −38.2754 dB to −46.7865 dB (14.09% less in linear form). This case concludes that on the DC side a dramatical improvement can be achieved whenever the noisy measurements are close to the converter. In (b) it shows a slightly similar impact on the AC side, the MAE has decreased by −0.6857 dB.
The DC system states estimation are shown in
Table 5. For the three different scenarios, the estimations are almost identical, and the error is lower than 5 digits (
). That is due to the low number of the DC buses and the linear feature of the DC components.
The AC side voltages and phase angles (in radian) estimations are available in
Appendix B.1 for the three scenarios. The unified approach has shown better estimates at bus 7 and 8 compared to the other scenarios. However, the overall error is expected to be low since the AC systems do not have many buses and branches.
Figure 10 shows the relative error in the phase angles for each scenario. In (a), the P-coupling and unified approaches were able to provide 5-digits accuracy in the phase angles compared to the decentralized approach. The unified approach was able to preserve the power coupling corrections on the measurements level while improving the estimates of the voltage states using the voltage coupling (e.g., Bus 7 in
Table A3).
5.1.2. Cigre B4 AC/DC/AC Network
The second network case is a modified version of the Cigre B4 DC grids test system [
57]. It has 6 AC systems—22 buses-, 2 DC systems—15 buses-, 8 AC generators (onshore and offshore), 11 converters (rectifiers and inverters) and 5 demand nodes/buses as shown in
Figure 11. A total of 109 measurements (
z) were used in the unified WLS, divided into 66 AC, 21 DC and 22 converter measurements. Note that the decentralized approach used AC and DC measurements only. Gaussian noise was added to all measurements except the slack buses voltages. Measurements types and counts are shown in
Table 6.
Table 7 shows the estimates of the converter measurements and constraints. A similar conclusion to the previous test case can be found, the power coupling minimizes the error in the power measurements near the converter (injected
and
). However, in some cases that correction can lead to errors in the voltage magnitude at both sides of the converter, which the
voltage coupling can correct.
In
Figure 12, the unified WLS has estimated the AC and DC measurements better than the decentralized approach, the MAE has reduced by −0.7577 dB and −3.5393 dB respectively. This pattern of accuracy performance in the Cigre B4 model is similar to the previous test case in
Figure 9.
In addition, the unified WLS has shown relatively low errors in the estimates of the power couplings constraints and the voltage coupling measurements as shown in
Figure 13a,b respectively. In (a), the estimated constraint values are compared with zero, while in (b) the estimated
are compared with the true measurements.
Table 8 shows the DC systems states estimation results for the three scenarios respectively, where T1 is: Decentralized, T2: With converter P-coupling, and T3: Unified. It can be noticed that the error in the decentralized method is higher than the other two scenarios (e.g., bus 6, 11, 12 and 13). In addition, whenever the system has many buses but less trusted redundant measurements, the possibility of having errors is more elevated.
Appendix B.2 contains the numerical values of the AC side estimations for the three scenarios in tables. In this test case, the 2nd AC system has the largest errors, mainly on the phase angle states. That is because the only trusted measurement is the voltage magnitude of the slack bus (#5), while the rest of the measurements are noisy. However, the unified state estimation was able to reduce that error compared to the other two scenarios, due to the converter coupling -redundant- measurements and constraints.
Figure 14 shows the relative error for each scenario in two AC systems with the largest errors in phase angles (radian) up to 5 digits.
The time performance of the three scenarios was calculated for 100 simulations on Cigre B4 test case. It was concluded that the unified WLS time performance was a downside for this approach.
Figure 15 shows that the decentralized approach was 62% faster than the unified method. In addition, the decentralized approach showed fewer memory allocations and storage consumption, as shown in
Table 9. These outcomes were expected since the unified approach has 65.1% more measurements to be processed. However, a pre-processing algorithm can build a trade-off between the estimation accuracy and the computational performance by selectively pick the most dominant measurements. In other words, the unified WLS computational drawback can be minimized by having sufficient and the most effective measurements.
6. Main Conclusions
The work in this paper presents a unified state estimator toolbox based on WLS algorithm that is suitable for hybrid VSC based HVDC/AC transmission networks. The estimator solves a centralized Newton-Raphson that uses AC, DC and converter measurements to provide better system states estimates for the system operators. The modified method has extended the WLS measurements function and the Jacobian matrix to include VSC-HVDC/AC coupling components, such as power constraints and voltage ratios. The algorithm demands that the AC, DC and converter measurements to be available in synchronized time. The approach was simulated in two different network cases with multiple AC and DC systems. A comparison between the unified WLS and the decentralized method was investigated. It was shown that the accuracy of the system states estimations has been improved and benefited from the converter coupling -redundant- measurements. Moreover, the overall errors in the measurements have been reduced compared to the decentralized WLS. However, higher computational power was required for the unified WLS and the power coupling method. Besides, it is worthy to note that the current stage of the algorithm is still missing the robustness feature against bad data.