Next Article in Journal
Magnetohydrodynamics Williamson Nanofluid Flow over an Exponentially Stretching Surface with a Chemical Reaction and Thermal Radiation
Previous Article in Journal
Non-Ferrous Metal Price Point and Interval Prediction Based on Variational Mode Decomposition and Optimized LSTM Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parameter Tuning for WECC Generic Photovoltaic System Models Using Latin Hypercube Sampling and Pareto Optimality

Department of Electrical Engineering, National Taiwan University of Science and Technology, Taipei 10607, Taiwan
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(12), 2739; https://doi.org/10.3390/math11122739
Submission received: 9 May 2023 / Revised: 7 June 2023 / Accepted: 14 June 2023 / Published: 16 June 2023

Abstract

:
With attention increasing towards worldwide ecological conservation, the energy policies of various countries are moving towards renewable energy systems, primarily with solar and wind as the primary power sources. The Western Electricity Coordinating Council (WECC) has developed a generic model as a simulation framework for a typical solar wind generation system to analyze the response of the actual equipment. In this paper, we propose a parameter tuning process for a generic model controller that combines Latin hypercube sampling and Pareto optimization to enable a generic model to effectively match the output of an actual inverter device. Different operating scenarios were considered to match the dynamic characteristics of the inverter device using the multi-objective particle swarm algorithm (MOPSO) along with Pareto optimization. The effectiveness of the proposed parameter tuning process was verified on a user-defined test system and, later, the IEEE 39 bus system with a large solar power facility.

1. Introduction

The rise of major global issues, such as climate change, environmental protection, the depletion of fossil fuels and global warming, has resulted in the topic of renewable energy receiving considerable attention from experts [1]. Accordingly, renewable power generation systems have opened a vast scope of fields for researchers [2]. Currently, solar and wind power have the capability to replace traditional fossil fuels in future power systems [3]. In the latest International Energy Agency (IEA) report, globally, the renewable energy sector was predicted to increase almost 60% from 2020 to 2026. Solar and wind systems account for more than 50% of power generation [4]. Taiwan, for example, has abundant wind energy and solar energy resources owing to its geographical location [5]. In view of this, the Taiwan government plans to expand renewable energy generation to a total capacity of 27,423 MW by 2025, with renewable energy accounting for 20% of the total electricity generation. However, installations of large-scale renewable sources into the existing power systems may pose several challenges to the system operations [6], owing to uncertainties in solar and wind power generation caused by climatic conditions [7,8]. Moreover, a risk of unstable system frequency exists when there is an imbalance between power supply and load demand [9]. In addition, the low inertia of renewable energy generation leads to poor frequency responses during abnormal events in power systems, affecting their stability [10,11]. Therefore, system operators must conduct stability analysis and dynamic assessment in accordance with the relevant grid codes and regulations before large-scale renewable energy sources are interconnected to the power grid [12,13]. In recent years, many studies have been performed on renewable systems and their advanced control theories. To assess the stability of inverter-based devices, the authors of [14] modeled a hybrid solar–wind generation system and presented its impact on power grid integration. In [15], different mathematical models of photovoltaic (PV) systems were developed based on power flow, harmonic and electromagnetic transient responses. In [16], a solar power generation model with a phase-locked loop (PLL) and the LCL filter was proposed to achieve low voltage ride-through (LVRT) capability. For examining the effects of solar power on distribution systems, Refs. [17,18] developed a tool to model solar power generation systems via the electromagnetic transient (EMTP) software and analyzed the model response in the Open Distribution System Simulator (OpenDSS) software. However, using different models from different equipment vendors may complicate the fast-screening process for power system operators. The Western Electricity Coordinating Council (WECC) and the International Electrotechnical Commission (IEC) have focused on developing generic models, as in [19,20], to enable power system operators and researchers to analyze the renewable power generation systems on a common simulation platform. Similarly, in [21], the Electric Power Research Institute (EPRI) has produced a technical guide to develop generic models based on the WECC. The generic models proposed by the WECC are already used in commercial software packages, such as PLSF, PSS/E, PowerWorld and DIgSILENT. In [22], wind turbine models based on the IEC and the WECC were proposed, and the compatibility of these two models was verified by adjusting the parameters between them. In [23], the development of solar dynamic models was introduced, and the dynamic response of the WECC model was verified with real-time test data. In [24], the output response of wind turbines was recorded via on-field measurements, and the actual measured data were verified using the WECC model. In [25], the WECC model was used to simulate a practical large-scale solar power system, and the frequency response of the power grid was successfully predicted under a contingency event. In [26], the WECC model was used to test the LVRT ability of large-scale solar power generation systems. Dynamic indices were used to evaluate system responses under different models to improve voltage stability.
Dynamic models are effective in assessing the response and stability of practical generation systems. Regarding the dynamic model configuration, the controller gains, time constants and other control parameters must be carefully tuned to ensure the stable operation of the generation systems and the responses of the actual equipment. However, the parameter tuning process is often performed by trial and error, and considerable time is required to tune the controller parameters. In electrical engineering, controller gain tuning is an optimization problem for determining the most suitable combination of controller gains. In related problems, optimization algorithms such as the genetic algorithm (GA), particle swarm optimization (PSO), ant colony optimization [27], Rao optimization [28] and the bat algorithm (BA) are used to search for solutions that are compatible with the operational situations. In [29], a lighting attachment procedure optimization (LAPO) algorithm was proposed to optimize the system voltage and current output responses for tuning PID controller gains in hybrid renewable energy systems. In [30,31], the design of the controller parameters of a static var compensator (SVC) and power stabilizer (PSS) was proposed using a GA and PSO. In addition, the rotor stability and dynamic stability of the power systems were considered. In [32], the Ziegler–Nichols PSO (ZN-PSO) algorithm was proposed for designing the dynamic parameters of the PI controller in the WECC generic model. The design results were validated against the actual measurement data.
However, traditional algorithms, which convert multiple objectives into a single-objective problem using a weighting approach, may not satisfy all objective criteria. In [33], NSGA-II was used to design the controller parameters of the static synchronous compensator (STATCOM) based on voltage stability, transient stability and dynamic robustness by considering the uncertainties of wind power generation. Therefore, to fit the expected response under different operating conditions, this paper proposes a parameter tuning process that combines multi-objective particle swarm optimization (MOPSO) with Pareto optimality. Generally, optimization techniques such as PSO and the GA often become stuck as the local optima or local best. In the proposed method, MOPSO is also equipped with the process of crowd distance calculation and the use of an external archive, which helps the optimization method to jump out of the local best and approach the global best. According to the characteristics of different response zones, a combination of controller parameters was selected using the minimum Manhattan distance method. The contributions of this study are as follows:
(1)
The parameter tuning of the WECC generic model controller was performed using MOPSO along with Pareto optimality, and the final solutions from all Pareto fronts were determined using the minimum Manhattan distance method.
(2)
Latin hypercube sampling (LHS) was used instead of the conventional random sampling method to increase the diversity of the variable initialization.
(3)
The proposed method was verified in an equivalent test system and the IEEE 39 bus system, considering six cases under two operating scenarios. Notably, the tuned parameters of any case could simultaneously satisfy different scenarios.
The remainder of this paper is organized as follows. Section 2 describes the multi-objective optimization algorithm, including Pareto optimality and hypercube sampling. Section 3 illustrates the proposed parameter tuning process and presents a WECC generic model. Section 4 describes the test system and presents the simulation results. Section 5 shows the application of the proposed method to the standard IEEE 39 bus system. Finally, conclusions are presented in Section 6.

2. Multi-Objective Optimization Algorithm

The optimization method follows a set of constraints to reach a desired single objective or goal. With more complex optimization problems, the number of objectives has increased. This shift from a single to multilevel problem is termed a multi-objective optimization problem. However, dealing with a multi-objective problem is always more complicated than dealing with a single-objective problem. This problem is often solved by converting multiple objectives into a single-objective function. A commonly used conversion method is the weighting sum method in which each objective is assigned a weighting factor in advance. The objective functions are weighted and combined into a single-objective function. However, the weighting sum method often provides unsatisfactory results for individual objectives, mainly because of the uneven and unknown selection of the weighting factors. The proposed parameter tuning algorithm can be divided into two steps: (1) parameter optimization and (2) parameter selection. MOPSO is used to generate a Pareto front set for parameter optimization. For parameter selection, the minimum Manhattan distance method is used to select the final solution from the feasible solution set. The framework of the proposed decision-making process is shown in Figure 1.

2.1. Latin Hypercube Sampling (LHS)

Generally, PSO can generate the first population using general random sampling (GRS). However, simple random sampling can lead to uneven sampling results. As shown in Figure 2 and Figure 3, LHS [34] is a stratified sampling concept. Each variable is divided into n non-overlapping intervals, and the parameter values are randomly selected in each interval, according to the probability density. Compared with GRS, LHS can increase the diversity of variable initialization in the sampling space and help avoid the over-concentration of sampling data in a specific area.
Based on LHS, the vector set of the initial random variables X can be expressed in (1) as
  X = x 1 , x 2 , , x m
Equation (2) is the cumulative distribution function of the initial set of random variables X:
F x = P r X x
For any discrete random variable, the probability mass function is expressed as (3). If there are m variables in the set of variables X, then the probability condition of the discrete random variables in (4) must be satisfied:
p x = P r X = x
i = 1 p x i = 1 , i   1 , 2 , m
Using (5), a sample of each variable at different intervals can be determined:
x i , j = F 1 δ j i U i , j n ,   j   1 , 2 , n
where x i , j is the ith variable in the jth sampling interval, m is the total number of variables, n is the number of sampling intervals, F 1 x is the inverse function of the cumulative distribution, δ j i is a random permutation of 1 , 2 , n , and U i , j are random variables generated by the standard uniform distribution.

2.2. Pareto Optimality

In multi-objective problems, conflicting solutions may occur when the objective functions are interdependent or proportional. In the Pareto optimality process, all feasible solutions are recorded using an external archive [35,36]. For multi-objective minimization, the individual solutions are compared for different objective functions using (6). If a set of feasible solutions outperforms other sets of feasible solutions in all objectives, it can be considered as a set of non-dominated solutions:
f n x < f n x ,             n 1 , 2 , N
where x X and x X are feasible solutions. If the fitness of x is better than that of x for all objectives, x is called the dominant solution, and x is called the dominated solution. Conversely, x is called a non-dominated solution if none of the existing solutions is better than the dominant solution. All non-dominated solutions regarded as Pareto fronts in the solution set space can be represented as (7):
P S A = x p 1 , x p 2 , x p i ,             i 1 , 2 , I
where x p i is the ith non-dominated solution, and I is the total number of non-dominated solutions.

2.3. MOPSO Algorithm

The traditional PSO is designed with single-objective optimization and cannot be directly applied to multi-objective problems. In this study, MOPSO was adopted for multi-objective problems to avoid conflicting solutions owing to its dependency on individual objectives [37]. In MOPSO, the search capability of the swarm-based algorithm is combined with the Pareto optimality properties. The nondominated solutions are recorded in the external archive during the iterative process. Any nondominated solution in the external archive can be selected as the global best (Gbest). The velocity and position of the particles are updated using (8) and (9). For single-objective PSO, there is only one global best solution; however, for MOPSO, the global best solution is a set of solutions. Particle swarms gradually move towards the main goal, moving through the local and global best solutions. During the iterative process, the external archive accumulates all nondominated solutions. The final result is the Pareto front solution. In addition, the crowding distance procedure is used to select the appropriate global best solution and to evaluate the insertion and deletion of non-dominated solutions, especially when the external archive is full. Figure 4 shows the flowchart of the proposed MOPSO algorithm.
v i k + 1 = w × v i k + c 1 × r a n d 1 × P b e s t i x i k + c 2 × r a n d 2 × G b e s t x i k
x i k + 1 = x i k + v i k + 1
where P b e s t i is the best solution of the ith particle itself, Gbest is the best solution of the group, w is the inertia weight, k is the kth iteration, c 1 and c 2 are acceleration constants, r a n d 1 and r a n d 2 are random variables between 0 and 1, x i k is the original position of the ith particle, x i k + 1 is the new position of the ith particle, v i k is the original velocity of the ith particle, and v i k + 1 is the new velocity of the ith particle.

2.4. Manhattan Distance Method

The Manhattan distance is the sum of the absolute axis distances of two points on a standard coordinate system. For multi-objective optimization problems, selecting the final solutions from the Pareto front set is called multiple-criteria decision-making. The minimum Manhattan distance (MMD) method selects a solution with the minimum distance from the normalized ideal vector as the final solution [38]. The Pareto front set P S A is obtained by MOPSO. The minimum of each objective function is calculated using (10), where x p i P S A is the ith solution in the Pareto front set:
n = min x p i P S A f n x p i
Similarly, the difference between the maximum and minimum of each objective function is calculated using (11):
L n = max x p i P S A f n x p i min x p i P S A f n x p i
The ideal vector y o p t is defined as follows:
y o p t = l 1 L 1       l 2 L 2             l N L N     T
The normalized fitness vector y x p i for each objective is calculated using (13):
y x p i = f 1 x p i L 1       f 2 x p i L 2             f N x p i L N     T
Following (12) and (13), the minimum Manhattan distance is defined as follows:
M N 1 = y x p i y o p t 1
M M D = min x p i P S A M N 1
where · 1 is the 1-norm, also known as the Manhattan norm.

3. Proposed Parameter Tuning Methodology

To achieve the desired response of PV systems, power grid engineers must simulate the equipment behaviors and tune their controller parameters before the PV systems are connected to the existing power system. However, to satisfy both the actual responses of PV systems and the safety of power grid operation, operators must adjust the controller parameters to meet these two conditions continuously. To improve parameter tuning efficiency, this study proposes a parameter tuning method based on MOPSO and Pareto optimality.

3.1. Test System

The test system used in this study is shown in Figure 5. Section A is the equivalent circuit of a large-scale PV power system connected to the Section B grid via an equivalent transmission line from the point of interconnection (POI). During the parameter tuning process, contingency events result in various voltage disturbances in power systems. A contingency study was conducted to confirm the security of the PV connected to the power grid.

3.2. WECC Generic Model

During the design phase of a PV generation system, dynamic models are used to ensure that the grid-connected system satisfies the grid codes. Several inverters have been used in practical power systems. User-defined models are required for various inverters on the corresponding simulation platform. As many PV systems are integrated into the power grid, it is difficult to manage different types of user-defined models and parameters. Therefore, generic models are used to simulate the actual response of user devices in the power grid.
The generic renewable energy generation models developed by WECC can effectively simulate the actual response in dynamic simulations. The WECC generic model, shown in Figure 6, consists of a generator model (REGCA), an area control model (REECB) and a central control model (REPCA), which can simulate the response of various inverters in practical power systems by adjusting the model parameters. Each model can be described as follows:
  • REGCA model: The REGCA model represents the generator model of the PV generation systems for generating Iq and Ip based on Iqcmd and Ipcmd, as shown in Figure 6. Furthermore, it maintains the voltage profile of the power systems within a predefined range based on Vt.
  • REECB model: As shown in Figure 7, the REECB model can simulate inverter behavior under different system states. The inverter controller strategies can be divided into real and reactive power control. Real power control is determined by the control signal at the previous level and remains between the maximum and minimum references. Reactive power control includes regulating reactive powers by iqinj during high- and low-voltage disturbance events and stabilizing the output using the PI controller under normal scenarios. The model output can be used to approximate the system responses by adjusting the controller parameters.
  • REPCA model: As shown in Figure 8, the REPCA model primarily receives the feedback response from the POI. The real P/F and Q/V control modes of all inverters in the PV systems can be realized using the REPCA model. The bus voltage and reactive power flow at the POI can be regulated in the reactive power control mode, in which PI controllers provide a primary frequency response according to the system frequency variations.

3.3. Problem Description

Under different operating scenarios, the generic model may provide real and reactive power responses to power grids according to voltage deviations. It can also be used as a fault compensator or a stability controller. The dynamic response of the test system with a three-phase short-circuit fault for 5 s at the POI is illustrated in Figure 9. Massive fluctuations were observed in both the real and reactive powers. This study proposes tuning the controller parameters to meet the preset reference powers by the system operator in both the shaded zones, as shown in Figure 9.
To tune the control parameters, the objective functions, equation constraints, inequality constraints and solution sets can be expressed as follows:
min     f n x , u ,       n 1 , 2 , N
subject   to   g i x , u = 0 ,       i 1 , 2 , , I h j x , u 0 ,       j 1 , 2 , , J  
x X
where   f n x , u is the nth objective function; g i x , u is the ith equation constraint; h j x , u is the jth inequality constraint; x and u are the control and dependent variables, respectively.

3.3.1. Control Variables

The large-scale PV model developed by the WECC performs real and reactive power–control calculations based on the specified bus voltage. To obtain a stable output response, engineers may adjust the proportional and integral gains of the PI controllers at each stage to achieve the desired response. In this study, the parameters of the PI controllers in REECB and REPCA for Q/V control were designed to satisfy the performance of the reference responses under different operating voltage scenarios, as shown in Figure 7 and Figure 8 within the red dashed blocks. The control variable x and dependent variable u based on MOPSO can be expressed as (19) and (20), respectively:
x = K v p , K v i , K q p , K q i , K p , K i T
u = P p v , 1 , P p v , i , Q p v , 1 , Q p v , i T
where K v p , K v i , K q p , K q i are the parameters of the Q/V control PI controller in REECB; K p , K i are the parameters of the Q/V control PI controller in REPCA; and   P p v , j and Q p v , j are the real and reactive power responses for the ith voltage operation scenario, respectively.

3.3.2. Objective Functions

To achieve the objective of tuning the PI controller parameters under different operating voltage scenarios, in this study, the real and reactive power dynamic responses were used as the objective functions. As shown in Figure 9, the dynamic features of the real power responses ( P p v , i o c c u r t , P p v , i c l r t ) and reactive power responses ( Q p v , i o c c u r t , Q p v , i c l r t ) under fault disturbances during the occurrence period ( t s t r t , 1 t e n d , 1 ) and recovery period ( t s t r t , 2 t e n d , 2 ) were used to minimize the average absolute error between the power responses of each shaded zone (see Figure 9) and the reference powers ( P p v b a s e , i o c c u r t , P p v b a s e , i c l r t , Q p v b a s e , i o c c u r t and Q p v b a s e , i c l r t ). The objective functions used in this study are shown in (22)–(25):
min     f p p v , i o c c u r ,   f Q p v , i o c c u r , f p p v , i c l r , f Q p v , i   c l r
f p p v , i o c c u r = 1 N k t = t s t r t , k t e n d , k P p v , i o c c u r t P p v b a s e , i o c c u r t
f Q p v , i o c c u r = 1 N k t = t s t r t , k t e n d , k Q p v , i o c c u r t Q p v b a s e , i o c c u r t
f p p v , i c l r = 1 N k t = t s t r t , k t e n d , k P p v , i c l r t P p v b a s e , i c l r t
f Q p v , i c l r = 1 N k t = t s t r t , k t e n d , k Q p v , i c l r t Q p v b a s e , i c l r t
where i denotes different operating voltage scenarios; N k is the total number of samples considered in the simulation time, which represents the number of samples taken at different times in 1 , 2 , K .
The sampling time can be represented by t s t r t , k , t s t r t , k + t s t e p , t s t r t , k + 2 t s t e p , , t s t r t , k + N k t s t e p . The variable t s t r t , k is the sampling start time, t e n d , k is the sampling end time, and t s t e p is the sampling interval of the simulation process; i.e., t e n d , k = t s t r t , k + N k t s t e p . In Figure 9, the first and second sampling intervals are represented by   N 1 and N 2 . From N 1 and N 2 , the average absolute errors between the power responses obtained by the WECC model and the desired references for different operating voltage scenarios can be calculated to obtain a feasible solution.

3.3.3. Equality Constraints

For any power system, there must be a balance between power generation and load demand, as follows:
P G i = P L i + V i j = 1 , j i n V j ( G i j cos θ i j + B i j sin θ i j )
Q G i = Q L i + V i j = 1 , j i n V j ( G i j sin θ i j B i j cos θ i j )
where P G i and Q G i are the active and reactive powers at bus i , V i is the voltage at bus i , V j is the voltage at bus j , and n is the number of buses; P L i and Q L i are the active and reactive power demands at bus i , respectively, G i j and B i j are the conductance and inductance between bus i and bus j , respectively, and θ i j is the impedance angle between bus i and bus j .

3.3.4. Inequality Constraints

Considering the planning and design, the network components in power systems must operate within a specified range, including the output limits of each unit, bus voltage limits, and line capacity limits, all of which must satisfy the maximum and minimum limits, as follows:
P G i m i n P G i P G i m a x         i 1 , 2 , , N G e n
Q G i m i n Q G i Q G i m a x       i 1 , 2 , , N G e n
V B u s i m i n V B u s i V B u s i m a x       i 1 , 2 , , N b u s
L i m i n L i L i m a x       i 1 , 2 , , N l i n e
where N G is the number of generating equipment, N b u s is the total number of buses, N l i n e is the total number of lines, P G i m i n and P G i m a x are the upper and lower limits of active power for the i th generating unit, Q G i m i n and Q G i m a x are the upper and lower limits of reactive power for the i th generating unit, V B u s i m i n and V B u s i m a x are the upper and lower limits of the allowable operating voltage of bus i , and L i m i n and L i m a x are the upper and lower limits of the thermal capacity for the i th transmission line, respectively.
Additionally, to ensure that the active and reactive power responses can effectively match the response references, fitness should be less than the average tolerance threshold, as follows:
0 f ε ,   f   f p p v , i o c c u r ,   f Q p v , i o c c u r , f p p v , i c l r , f Q p v , i   c l r
where ε is the average tolerance threshold.

3.4. Application of MOPSO to Parameter Tuning

The parameters of a dynamic model often require repeated tuning. In this paper, a tuning process is proposed to search for the Pareto front solutions of dynamic parameters using MOPSO. The most suitable solution is selected using the minimum Manhattan distance method. A flowchart of the proposed parameter tuning method is shown in Figure 10 with the following steps:
PART (A) Initialization
Step (A-1)
System initialization: The system architecture, line parameters, scenarios, objectives and constraints are initialized.
Step (A-2)
Control variable initialization: The control variables in (19) are initialized based on Latin hypercube sampling.
Step (A-3)
Dynamic simulation: Control variables are used for the dynamic simulations by considering different operating scenarios. The dependent variables in (20) are obtained from the simulation results.
Step (A-4)
Fitness calculation: The fitness values are calculated using (22)–(25) using the dependent variables obtained in Step A-4.
Step (A-5)
Initial external archive: According to Step A-4, the set of Pareto front solutions is obtained using (6) in the external archive.
PART (B) MOPSO algorithm
Step (B-1)
Gbest selection: The best solution is selected from the external archive using the crowding distance procedure.
Step (B-2)
PSO algorithm: As shown in Figure 4, the control variable x is updated according to (8) and (9), and the new fitness values are evaluated.
Step (B-3)
Pareto optimality: The fitness generated in Step B-2 is compared with the non-dominated solutions in the external archive, and the external archive is updated if required.
Step (B-4)
Termination condition: The iterative process is stopped if a solution with all fitness values less than the preset threshold exists in the external archive.
PART (C) MMD decision-making
Step (C-1)
MMD solution selection: The final result is selected from the final external archive using (10)–(15).

4. Results and Discussion

To verify the performance of the proposed parameter tuning method, the responses of the controller parameters were simulated for different operating scenarios.

4.1. Scenarios

Figure 11 and Figure 12 are the reference responses under different operating scenarios set in this study, which can also be regarded as field measurements. Table 1 lists the operating scenarios of the test system, and Table 2 lists the controller parameters under the reference response.
In this study, based on the reference responses (field measurements), several different design scenarios were considered to validate the proposed parameter tuning process. The different controller parameters and cases are listed in Table 3.
The descriptions of the cases are as follows:
Case (1)
The controller parameters for the V/Q control in the REECB model are tuned. In other words, some parameters are used to adjust the responses of the PV systems.
Case (2)
The controller parameters for V/Q control in REECB and REPCA are tuned. In other words, all parameters are used to adjust the responses of the PV systems.
Case (3)
Some parameter settings for the inverter are already known from the device manual for applications. In this case, the controller parameters of the V/Q controls in the REECB and REPCA models are tuned using the known parameters of the PI controllers.
Case (4)
The limitations between the WECC generic model and the practical device are considered. The differences between the control schemes of the WECC generic model and the practical device are examined. In this case, the controller parameters of the V/Q control in the REECB and REPCA models are tuned with a false parameter of the PI controller.
Case (5)
The controller parameters for V/Q control in REECB and REPCA are tuned using only two dynamic features, f p p v , i o c c u r and f Q p v , i o c c u r , considering the occurrence period of the fault disturbances.
Case (6)
Similarly to Case 5, the controller parameters for V/Q control in REECB and REPCA are tuned using f p p v , i c l r and f Q p v , i c l r , considering the recovery period of fault disturbances.

4.2. Test System and Parameter Setting

The different operating scenarios in the test system shown in Figure 5 were used to verify the performance of the proposed optimization strategy to match the reference response based on the dynamic features of the active and reactive powers in different regions. Table 4 lists the relevant parameters and operating scenarios for the MOPSO algorithm.
In this study, the proposed parameter tuning strategy was modeled on a power system simulator for an engineering (PSS®E) dynamic simulation platform. Controller parameter optimization was implemented using MATLAB. A large-scale PV power system of 100 MW was modeled as shown in Figure 5. The WECC generic PV system was modeled under the V/Q control mode to tune the local control model (REECB) and central control model (REPCA) parameters.
Considering the stability of the output response, the PI controller parameters mentioned in the previous section were tuned, and the range of the controller variables is listed in Table 5. The other parameters in the general model were obtained directly from the WECC configuration [19,20,21]. The parameter tuning process is shown in Figure 10, and the two operating scenarios are listed in Table 1. The fault disturbance occurred at the POI at exactly 1 s and cleared at 1.5 s, and the dynamic simulation continued for 3 s.

4.3. Test Results

The well-tuned controller parameters according to the proposed method for different cases are listed in Table 6. The fitness values of the objective functions are listed in Table 7 according to the active and reactive power responses in different characteristic zones. Based on the responses of Figure 13, Figure 14, Figure 15 and Figure 16 in different cases, the following conclusions can be drawn:
(1)
The response results and fitness values in Cases 1 and 2 matched the reference response in the base case. The tuned controller parameters approximate the control parameters of the base response and thus can be used as a benchmark for inverter design.
(2)
From the results of Case 3, it can be observed that in practical applications, if some parameters of the inverter are known, the reference values can be effectively tracked, and the remaining control parameters can be tuned.
(3)
Case 4 shows that if the equipment conditions limit the tuning parameters, the difficulty of the parameter tuning may increase.
(4)
Cases 5 and 6 were tuned using only two dynamic features under the occurrence and recovery periods of fault disturbances. The difficulty of parameter tuning may have increased in Case 5.
Generally, in multi-objective problems, PSO evaluates the fitness value using the weighting sum method. The weighting sum method involves assigning weights to each single-objective function. In this case, each objective function value is assigned an equal weight for equal representation. MOPSO continuously updates the Pareto front solutions in the external archives according to Pareto optimality, and the best solutions are selected. To compare the searchability of the multi-objective and single-objective algorithms, the MOPSO and PSO algorithms were compared by normalizing the results of the original objective functions with equal weighting to all objectives for Cases 1 and 2. Table 8 lists the controller parameters tuned using MOPSO and PSO, and Table 9 lists the fitness values calculated for MOPSO and PSO with different objective functions. Figure 17 and Figure 18 show the performance comparison for Cases 1 and 2, respectively. This shows that, owing to the continuous updating and replacement of the Pareto front solutions, MOPSO has better search performance and convergence than PSO and can be effectively applied for tuning controller parameters.

5. Applications

To test the parameter tuning method proposed in this paper for practical application, the IEEE 39 bus power system shown in Figure 19 was used as the sample system, and a 100 MW PV power plant was installed at Bus 2. The effects of the reference cases in Table 1 on the IEEE 39 bus system are shown in Figure 20 and Figure 21. The controller parameters were adjusted based on the base scenarios and cases, as described in Table 3.
The tuned controller parameters in the IEEE 39 bus system are listed in Table 10 according to the different cases. The fitness values of the objective functions are tabulated in Table 11 according to the active and reactive power responses in the different feature regions.
According to the response and tuned parameters for Figure 22, Figure 23, Figure 24 and Figure 25 in different cases, the observations are as follows:
(1)
Table 10 shows that the Case 1–Case 3 tuning results can approximate the reference cases. The dynamic response of the actual inverter in the power system can be simulated effectively.
(2)
In Case 4, if the equipment conditions limit the inverter responses, the difficulty in parameter tuning may increase.
(3)
Cases 5 and 6 are tuned by only two dynamic features under the occurrence and recovery periods of fault disturbances. The difficulty of parameter tuning in Cases 5 and 6 is more serious than that in Cases 1–3.
Table 12 shows the controller parameters tuned by MOPSO and PSO in Cases 1 and 2 for the IEEE 39 bus system. The fitness values obtained by MOPSO and PSO under different objectives are listed in Table 13. Figure 26 and Figure 27 show the performance comparisons for Cases 1 and 2, respectively. MOPSO has better search performance and convergence than PSO. Table 14 presents the GD and IGD values of the proposed method. It is observed that for both the systems and cases the performance values are quite low.

6. Conclusions

With a large number of renewable energy generation systems connected to power systems, a regressive stability analysis is imperative. During the simulation stage, the WECC generic model was used to evaluate system stability. To ensure that the model response can meet the expected response of the equipment, a generic model parameter tuning method was proposed using MOPSO combined with Pareto optimality. The effectiveness of the proposed method was tested on a user-defined test system and an IEEE 39 bus system. The major conclusions drawn are as follows:
(1)
MOPSO can effectively improve the ability to search for the best solution in the solution space by continuously updating and replacing an external archive.
(2)
The performance of each objective feature in MOPSO is better than that of the single-objective weighting sum method. The lowest fitness values for IEEE 39 using MOPSO and PSO were 6.82 × 10−10 and 7.34 × 10−6, respectively, whereas the highest value using MOPSO was 1.76 × 10−7, and the highest value using PSO was 4.88 × 10−4. In other words, the MOPSO method is superior to the standard PSO method.
(3)
For different cases with different parameter selections and objective features, the proposed method can effectively tune the parameters to approximate the expected response. For the IEEE 39 system, the level of error between the MOPSO-tuned parameter values and the reference values was under 1%.
(4)
The proposed method can operate efficiently under various voltage scenarios. In other words, the proposed method is flexible under different field measurement conditions.
As an extended work, we plan to design a new objective function in accordance with the control system. In addition, more response characteristics, operation scenarios and different control modes will be considered. Moreover, the acquisition of real-time inverter data is planned, which will improve the tracking and performance of the method even further.

Author Contributions

Conceptualization, N.-C.Y.; funding acquisition, N.-C.Y.; investigation, N.-C.Y. and C.-W.H.; methodology, N.-C.Y. and C.-W.H.; software, N.-C.Y. and C.-W.H.; supervision, N.-C.Y.; validation, N.-C.Y. and C.-W.H.; writing—original draft preparation, N.-C.Y., C.-W.H. and A.S.; writing—review and editing, N.-C.Y. and Abhilash Sen. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by the National Science and Technology Council (NSTC) in Taiwan (NSTC 111-2622-8-011-014-SB) and the DELTA-NTUST Joint Research Center.

Data Availability Statement

Not applicable.

Acknowledgments

This work was supported by the National Science and Technology Council (NSTC) in Taiwan. The authors are thankful for the funding provided by NSTC (NSTC 111-2622-8-011-014-SB) and Delta Electronics (DELTA)—National Taiwan University of Science of Technology (NTUST) Joint Research Center.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Qazi, A.; Hussain, F.; Rahim, N.A.; Hardaker, G.; Alghazzawi, D.; Shaban, K.; Haruna, K. Towards sustainable energy: A systematic review of renewable energy sources, technologies, and public opinions. IEEE Access 2019, 7, 63837–63851. [Google Scholar] [CrossRef]
  2. Rizzi, F.; van Eck, N.J.; Frey, M. The production of scientific knowledge on renewable energies: Worldwide trends, dynamics and challenges and implications for management. Renew. Energy 2014, 62, 657–671. [Google Scholar] [CrossRef]
  3. Sadeghi, H.; Rashidinejad, M.; Moeini-Aghtaie, M.; Abdollahi, A. The energy hub: An extensive survey on the state-of-the-art. Appl. Therm. Eng. 2019, 161, 114071. [Google Scholar] [CrossRef]
  4. Renewables 2021; IEA: Paris, France, 2021.
  5. Chen, H.H.; Lee, A.H.J.R.; Reviews, S.E. Comprehensive overview of renewable energy development in taiwan. Renew. Sustain. Energy Rev. 2014, 37, 215–228. [Google Scholar] [CrossRef]
  6. Alam, M.S.; Al-Ismail, F.S.; Salem, A.; Abido, M.A.J.I.A. High-level penetration of renewable energy sources into grid utility: Challenges and solutions. IEEE Access 2020, 8, 190277–190299. [Google Scholar] [CrossRef]
  7. Kumar, K.P.; Saravanan, B. Recent techniques to model uncertainties in power generation from renewable energy sources and loads in microgrids—A review. Renew. Sustain. Energy Rev. 2017, 71, 348–358. [Google Scholar] [CrossRef]
  8. Athari, M.H.; Wang, Z. Modeling the uncertainties in renewable generation and smart grid loads for the study of the grid vulnerability. In Proceedings of the 2016 IEEE Power & Energy Society Innovative Smart Grid Technologies Conference (ISGT), Minneapolis, MN, USA, 6–9 September 2016; IEEE: Minneapolis, MN, USA, 2016; pp. 1–5. [Google Scholar]
  9. Dreidy, M.; Mokhlis, H.; Mekhilef, S. Inertia response and frequency control techniques for renewable energy sources: A review. Renew. Sustain. Energy Rev. 2017, 69, 144–155. [Google Scholar] [CrossRef]
  10. Yu, S.; Zhang, L.; Lu, H.H.-C.; Fernando, T.; Wong, K.P. A dse-based power system frequency restoration strategy for pv-integrated power systems considering solar irradiance variations. IEEE Trans. Ind. Inform. 2017, 13, 2511–2518. [Google Scholar] [CrossRef]
  11. Johnson, S.C.; Papageorgiou, D.J.; Mallapragada, D.S.; Deetjen, T.A.; Rhodes, J.D.; Webber, M.E. Evaluating rotational inertia as a component of grid reliability with high penetrations of variable renewable energy. Energy 2019, 180, 258–271. [Google Scholar] [CrossRef]
  12. Etxegarai, A.; Eguia, P.; Torres, E.; Iturregi, A.; Valverde, V. Review of grid connection requirements for generation assets in weak power grids. Renew. Sustain. Energy Rev. 2015, 41, 1501–1514. [Google Scholar] [CrossRef]
  13. Lammert, G.; Heß, T.; Schmidt, M.; Schegner, P.; Braun, M. Dynamic grid support in low voltage grids—Fault ride-through and reactive power/voltage support during grid disturbances. In Proceedings of the 2014 Power Systems Computation Conference, Wroclaw, Poland, 18–22 August 2014; IEEE: Wroclaw, Poland, 2014; pp. 1–7. [Google Scholar]
  14. Zhu, S.; Piper, D.; Ramasubramanian, D.; Quint, R.; Isaacs, A.; Bauer, R. Modeling inverter-based resources in stability studies. In Proceedings of the 2018 IEEE Power & Energy Society General Meeting (PESGM), Portland, OR, USA, 5–10 August 2018; IEEE: Portland, OR, USA, 2018; pp. 1–5. [Google Scholar]
  15. Ćuk, V.; Ribeiro, P.F.; Cobben, J.F.; Kling, W.L.; Isleifsson, F.R.; Bindner, H.W.; Martensen, N.; Samadi, A.; Söder, L. Considerations on the Modeling of Photovoltaic Systems for Grid Impact Studies. In Proceedings of the 1st International Workshop on the Integration of Solar Power into Power Systems, Aarhus, Denmark, 24 October 2011; pp. 35–42. [Google Scholar]
  16. Nanou, S.I.; Papathanassiou, S.A. Modeling of a pv system with grid code compatibility. Electr. Power Syst. Res. 2014, 116, 301–310. [Google Scholar] [CrossRef]
  17. Mahseredjian, J.; Dinavahi, V.; Martinez, J.A. Simulation tools for electromagnetic transients in power systems: Overview and challenges. IEEE Trans. Power Deliv. 2009, 24, 1657–1669. [Google Scholar] [CrossRef]
  18. Hariri, A.; Faruque, M.O. A hybrid simulation tool for the study of pv integration impacts on distribution networks. IEEE Trans. Sustain. Energy 2016, 8, 648–657. [Google Scholar] [CrossRef]
  19. Ellis, A.; Behnke, M.; Elliott, R. Generic Solar Photovoltaic System Dynamicsimulation Model Specification; Sandia Report; Sandia National Laboratories: Albuquerque, NM, USA, 2013. [Google Scholar]
  20. Sørensen, P.; Andresen, B.; Fortmann, J.; Pourbeik, P. Modular structure of wind turbine models in iec 61400-27-1. In Proceedings of the 2013 IEEE Power & Energy Society General Meeting, Vancouver, BC, Canada, 21–25 July 2013; IEEE: Vancouver, BC, Canada, 2013; pp. 1–5. [Google Scholar]
  21. Farantatos, E. Model User Guide for Generic Renewable Energy Systemmodels; Tech. Rep. 3002014083; Electric Power Research Institute: Washington, DC, USA, 2018. [Google Scholar]
  22. Göksu, Ö.; Sørensen, P.; Morales, A.; Weigel, S.; Fortmann, J.; Pourbeik, P. Compatibility of iec 61400-27-1 and wecc 2nd generation wind turbine models. In Proceedings of the 15th International Workshop on Large-Scale Integration of Wind Power into Power Systems as well as on Transmission Networks for Offshore Wind Power Plants, Vienna, Austria, 15–17 November 2016; pp. 15–17. [Google Scholar]
  23. Eguia, P.; Etxegarai, A.; Torres, E.; San Martin, J.; Albizu, I. Modeling and validation of photovoltaic plants using generic dynamic models. In Proceedings of the 2015 International Conference on Clean Electrical Power (ICCEP), Taormina, Italy, 16–18 June 2015; IEEE: Taormina, Italy, 2015; pp. 78–84. [Google Scholar]
  24. Pourbeik, P.; Etzel, N.; Wang, S. Model validation of large wind power plants through field testing. IEEE Trans. Sustain. Energy 2017, 9, 1212–1219. [Google Scholar] [CrossRef]
  25. Machlev, R.; Batushansky, Z.; Soni, S.; Chadliev, V.; Belikov, J.; Levron, Y. Verification of utility-scale solar photovoltaic plant models for dynamic studies of transmission networks. Energies 2020, 13, 3191. [Google Scholar] [CrossRef]
  26. Lammert, G.; Premm, D.; Ospina, L.D.P.; Boemer, J.C.; Braun, M.; Van Cutsem, T. Control of photovoltaic systems for enhanced short-term voltage stability and recovery. IEEE Trans. Energy Convers. 2018, 34, 243–254. [Google Scholar] [CrossRef] [Green Version]
  27. Abubakar, H.; Abdullahi, M.; Smaiala, B. Ants colony optimization algorithm in the hopfield neural network for agricultural soil fertility reverse analysis. Iraqi J. Comput. Sci. Math. 2022, 3, 32–42. [Google Scholar] [CrossRef]
  28. Shamis, N.A.; Mohammad, A.; Hind Raad, I. Rao-svm machine learning algorithm for intrusion detection System. Iraqi J. Comput. Sci. Math. 2020, 1, 23–27. [Google Scholar] [CrossRef]
  29. Saleh, B.; Yousef, A.M.; Ebeed, M.; Abo-Elyousr, F.K.; Elnozahy, A.; Mohamed, M.; Abdelwahab, S.A.M. Design of pid controller with grid connected hybrid renewable energy system using optimization algorithms. J. Electr. Eng. Technol. 2021, 16, 3219–3233. [Google Scholar] [CrossRef]
  30. Fetouh, T.; Zaky, M.S. New approach to design svc-based stabiliser using genetic algorithm and rough set theory. IET Gener. Transm. Distrib. 2017, 11, 372–382. [Google Scholar] [CrossRef]
  31. Shahgholian, G.; Movahedi, A. Power system stabiliser and flexible alternating current transmission systems controller coordinated design using adaptive velocity update relaxation particle swarm optimisation algorithm in multi-machine power system. IET Gener. Transm. Distrib. 2016, 10, 1860–1868. [Google Scholar] [CrossRef]
  32. Soni, S. Solar pv Plant Model Validation for Grid Integration Studies; Arizona State University: Tempe, AZ, USA, 2014. [Google Scholar]
  33. Chi, Y.; Xu, Y. Multi-objective robust tuning of statcom controller parameters for stability enhancement of stochastic wind-penetrated power systems. IET Gener. Transm. Distrib. 2020, 14, 4805–4814. [Google Scholar] [CrossRef]
  34. Helton, J.C.; Davis, F.J. Latin hypercube sampling and the propagation of uncertainty in analyses of complex systems. Reliab. Eng. 2003, 81, 23–69. [Google Scholar] [CrossRef] [Green Version]
  35. Yang, N.C.; Le, M.D. Multi-objective bat algorithm with time-varying inertia weights for optimal design of passive power filters set. IET Gener. Transm. Distrib. 2015, 9, 644–654. [Google Scholar] [CrossRef]
  36. Yang, N.-C.; Le, M.-D. Optimal design of passive power filters based on multi-objective bat algorithm and pareto front. Applied Soft Computing 2015, 35, 257–266. [Google Scholar] [CrossRef]
  37. Reyes-Sierra, M.; Coello, C.C. Multi-objective particle swarm optimizers: A survey of the state-of-the-art. Int. J. Comput. Intell. Res. 2006, 2, 287–308. [Google Scholar]
  38. Chiu, W.-Y.; Yen, G.G.; Juan, T.-K. Minimum manhattan distance approach to multiple criteria decision making in multiobjective optimization problems. IEEE Trans. Evol. Comput. 2016, 20, 972–985. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Decision-making process framework.
Figure 1. Decision-making process framework.
Mathematics 11 02739 g001
Figure 2. General random sampling distribution.
Figure 2. General random sampling distribution.
Mathematics 11 02739 g002
Figure 3. Latin hypercube sampling distribution.
Figure 3. Latin hypercube sampling distribution.
Mathematics 11 02739 g003
Figure 4. Flow chart of the proposed MOPSO algorithm.
Figure 4. Flow chart of the proposed MOPSO algorithm.
Mathematics 11 02739 g004
Figure 5. Single line diagram for test system.
Figure 5. Single line diagram for test system.
Mathematics 11 02739 g005
Figure 6. WECC generic model block diagram.
Figure 6. WECC generic model block diagram.
Mathematics 11 02739 g006
Figure 7. REECB control block diagram.
Figure 7. REECB control block diagram.
Mathematics 11 02739 g007
Figure 8. REPCA control block diagram.
Figure 8. REPCA control block diagram.
Mathematics 11 02739 g008
Figure 9. Dynamic response of WECC model during fault disturbances.
Figure 9. Dynamic response of WECC model during fault disturbances.
Mathematics 11 02739 g009
Figure 10. Flow chart of the proposed parameter tuning method.
Figure 10. Flow chart of the proposed parameter tuning method.
Mathematics 11 02739 g010
Figure 11. Reference response of PV power generation under operating scenario 1.
Figure 11. Reference response of PV power generation under operating scenario 1.
Mathematics 11 02739 g011
Figure 12. Reference response of PV power generation under operating scenario 2.
Figure 12. Reference response of PV power generation under operating scenario 2.
Mathematics 11 02739 g012
Figure 13. Active power dynamic response in each case under operating scenario 1.
Figure 13. Active power dynamic response in each case under operating scenario 1.
Mathematics 11 02739 g013
Figure 14. Reactive power dynamic response in each case under operating scenario 1.
Figure 14. Reactive power dynamic response in each case under operating scenario 1.
Mathematics 11 02739 g014
Figure 15. Active power dynamic response in each case under operating scenario 2.
Figure 15. Active power dynamic response in each case under operating scenario 2.
Mathematics 11 02739 g015
Figure 16. Reactive power dynamic response in each case under operating scenario 2.
Figure 16. Reactive power dynamic response in each case under operating scenario 2.
Mathematics 11 02739 g016
Figure 17. Performance comparison in Case 1 using different algorithms.
Figure 17. Performance comparison in Case 1 using different algorithms.
Mathematics 11 02739 g017
Figure 18. Performance comparison in Case 2 using different algorithms.
Figure 18. Performance comparison in Case 2 using different algorithms.
Mathematics 11 02739 g018
Figure 19. IEEE 39 bus power system.
Figure 19. IEEE 39 bus power system.
Mathematics 11 02739 g019
Figure 20. Reference response of PV system installed in the IEEE 39 bus system under operating scenario 1.
Figure 20. Reference response of PV system installed in the IEEE 39 bus system under operating scenario 1.
Mathematics 11 02739 g020
Figure 21. Reference response of PV system installed in the IEEE 39 bus system under operating scenario 2.
Figure 21. Reference response of PV system installed in the IEEE 39 bus system under operating scenario 2.
Mathematics 11 02739 g021
Figure 22. Active power dynamic response in the IEEE 39 bus system under operating scenario 1.
Figure 22. Active power dynamic response in the IEEE 39 bus system under operating scenario 1.
Mathematics 11 02739 g022
Figure 23. Reactive power dynamic response in the IEEE 39 bus system under operating scenario 1.
Figure 23. Reactive power dynamic response in the IEEE 39 bus system under operating scenario 1.
Mathematics 11 02739 g023
Figure 24. Active power dynamic response of the IEEE 39 bus system under operating scenario 2.
Figure 24. Active power dynamic response of the IEEE 39 bus system under operating scenario 2.
Mathematics 11 02739 g024
Figure 25. Reactive power dynamic response of the IEEE 39 bus system under operating scenario 2.
Figure 25. Reactive power dynamic response of the IEEE 39 bus system under operating scenario 2.
Mathematics 11 02739 g025
Figure 26. Performance comparison of the IEEE 39 bus system in Case 1 using different algorithms.
Figure 26. Performance comparison of the IEEE 39 bus system in Case 1 using different algorithms.
Mathematics 11 02739 g026
Figure 27. Performance comparison of the IEEE 39 bus system in Case 2 using different algorithms.
Figure 27. Performance comparison of the IEEE 39 bus system in Case 2 using different algorithms.
Mathematics 11 02739 g027
Table 1. Description of the operating scenarios.
Table 1. Description of the operating scenarios.
Operating Scenarios
Scenario 1POI voltage dip to 0.9 pu
Scenario 2POI voltage dip to 0.35 pu
Table 2. Controller parameters under reference response.
Table 2. Controller parameters under reference response.
ParameterValue
K v p 0.353
K v i 155.388
K q p 47.552
K q i 3.076
K p 47.476
K i 319.045
Table 3. Tuned controller parameters under each case.
Table 3. Tuned controller parameters under each case.
ParameterREECBREPCA
K v p K v i K q p K q i K p K i
Case 1
Case 2
Case 3
(known)

(known)
Case 4
(false setting)
Case 5
Case 6
Table 4. MOPSO-related parameters and other variables of the proposed method.
Table 4. MOPSO-related parameters and other variables of the proposed method.
ParameterValueParameterValue
Population size20 t s t e p 1.66 × 10 3
LHS level20 t s t r t , 1 1
External archive size100 t e n d , 1 1.25
Maximum iteration500 t s t r t , 2 1.5
Feature tolerance ( ε )0.1 t e n d , 2 1.75
Simulation time3 s
Table 5. Feasible ranges of controller parameters.
Table 5. Feasible ranges of controller parameters.
ParameterValue
K v p 0–0.8
K v i 0–200
K q p 0–100
K q i 0–10
K p 0–100
K i 0–400
Table 6. Tuned controller parameters in different cases.
Table 6. Tuned controller parameters in different cases.
ParameterREECBREPCA
K v p K v i K q p K q i K p K i
Reference0.353155.38847.5523.07647.476319.045
Case 10.353155.38847.5523.078--
Case 20.353155.38847.5493.11547.476319.014
Case 30.353155.39147.5523.07647.482314.935
Case 40.500127.09177.2340.08549.56964.626
Case 50.353155.38847.0193.43917.407301.658
Case 60.353155.34947.4253.46517.9965.004
Table 7. Fitness values in different cases.
Table 7. Fitness values in different cases.
Fitness f p p v , 1 o c c u r f p p v , 1 c l r f Q p v , 1 o c c u r f Q p v , 1 c l r f p p v , 2 o c c u r f p p v , 2 c l r f Q p v , 2 o c c u r f Q p v , 2 c l r
Case 10.00 × 10+07.84 × 10−100.00 × 10+09.80 × 10−100.00 × 10+01.73 × 10−80.00 × 10+01.38 × 10−8
Case 20.00 × 10+04.55 × 10−81.56 × 10−97.08 × 10−80.00 × 10+04.25 × 10−80.00 × 10+06.15 × 10−8
Case 31.21 × 10−76.65 × 10−74.86 × 10−74.22 × 10−61.22 × 10−96.46 × 10−75.70 × 10−93.39 × 10−6
Case 48.72 × 10−41.47 × 10−35.49 × 10−38.96 × 10−38.46 × 10−52.75 × 10−32.45 × 10−41.77 × 10−2
Case 58.77 × 10−8-6.69 × 10−7-1.14 × 10−8-4.62 × 10−8-
Case 6-5.41 × 10−7-7.29 × 10−6-8.70 × 10−7-8.90 × 10−6
Table 8. Tuned controller parameters using multi-objective algorithms and single-objective algorithms.
Table 8. Tuned controller parameters using multi-objective algorithms and single-objective algorithms.
ParameterREECBREPCA
K v p K v i K q p K q i K p K i
Reference0.353155.38847.5523.07647.476319.045
Case 1MOPSO0.353155.38847.5523.078--
PSO0.351155.20630.5168.507--
Case 2MOPSO0.353155.38847.5493.11547.476319.014
PSO0.366156.03735.8869.89870.913399.223
Table 9. Fitness values in the test system using multi-objective algorithms and single-objective algorithms.
Table 9. Fitness values in the test system using multi-objective algorithms and single-objective algorithms.
Reference f p p v , 1 o c c u r f p p v , 1 c l r f Q p v , 1 o c c u r f Q p v , 1 c l r f p p v , 2 o c c u r f p p v , 2 c l r f Q p v , 2 o c c u r f Q p v , 2 c l r
Case 1MOPSO0.00 × 10+07.84 × 10−100.00 × 10+09.80 × 10−100.00 × 10+01.73 × 10−80.00 × 10+01.38 × 10−8
PSO5.78 × 10−66.61 × 10−63.53 × 10−54.44 × 10−51.43 × 10−69.80 × 10−64.05 × 10−66.65 × 10−5
Case 2MOPSO0.00 × 10+04.55 × 10−81.56 × 10−97.08 × 10−80.00 × 10+04.25 × 10−80.00 × 10+06.15 × 10−8
PSO2.99 × 10−56.00 × 10−52.01 × 10−43.74 × 10−43.87 × 10−51.67 × 10−41.09 × 10−41.14 × 10−3
Table 10. Tuned controller parameters of the IEEE 39 bus system under different scenarios.
Table 10. Tuned controller parameters of the IEEE 39 bus system under different scenarios.
ParameterREECBREPCA
K v p K v i K q p K q i K p K i
Reference0.353155.38847.5523.07647.476319.045
Case 10.353155.38847.5563.085--
Case 20.353155.38847.5533.07547.477319.075
Case 30.353155.38847.5523.07647.476319.046
Case 40.500143.54349.6953.39647.227326.763
Case 50.353155.40247.5373.72547.472318.500
Case 60.353155.38747.3506.74947.470318.340
Table 11. Fitness values of the IEEE 39 bus system in different design scenarios.
Table 11. Fitness values of the IEEE 39 bus system in different design scenarios.
Fitness f p p v , 1 o c c u r f p p v , 1 c l r f Q p v , 1 o c c u r f Q p v , 1 c l r f p p v , 2 o c c u r f p p v , 2 c l r f Q p v , 2 o c c u r f Q p v , 2 c l r
Case 19.43 × 10−81.13 × 10−71.86 × 10−73.85 × 10−76.82 × 10−107.49 × 10−89.74 × 10−111.76 × 10−7
Case 28.88 × 10−89.55 × 10−81.16 × 10−71.47 × 10−76.82 × 10−106.65 × 10−89.74 × 10−118.47 × 10−8
Case 30.00 × 10+08.27 × 10−80.00 × 10+01.91 × 10−80.00 × 10+03.92 × 10−100.00 × 10+04.90 × 10−10
Case 41.19 × 10−51.57 × 10−41.83 × 10−49.11 × 10−43.53 × 10−46.25 × 10−46.76 × 10−44.66 × 10−3
Case 51.13 × 10−7-2.59 × 10−7-1.07 × 10−7-1.71 × 10−7-
Case 6-4.14 × 10−7-3.02 × 10−6-4.41 × 10−7-2.47 × 10−6
Table 12. Tuned controller parameters in the IEEE 39 bus system using multi-objective algorithms and single-objective algorithms.
Table 12. Tuned controller parameters in the IEEE 39 bus system using multi-objective algorithms and single-objective algorithms.
ParameterREECBREPCA
K v p K v i K q p K q i K p K i
Reference0.353155.38847.5523.07647.476319.045
Case 1MOPSO0.353155.38847.5563.085--
PSO0.334157.85670.52710.00--
Case 2MOPSO0.353155.38847.5533.07547.477319.075
PSO0.290161.48646.0354.55048.333374.811
Table 13. Fitness values in the IEEE 39 bus system using multi-objective algorithms and single-objective algorithms.
Table 13. Fitness values in the IEEE 39 bus system using multi-objective algorithms and single-objective algorithms.
Fitness f p p v , 1 o c c u r f p p v , 1 c l r f Q p v , 1 o c c u r f Q p v , 1 c l r f p p v , 2 o c c u r f p p v , 2 c l r f Q p v , 2 o c c u r f Q p v , 2 c l r
Case 1MOPSO9.43 × 10−81.13 × 10−71.86 × 10−73.85 × 10−76.82 × 10−107.49 × 10−89.74 × 10−111.76 × 10−7
PSO1.30 × 10−52.51 × 10−58.64 × 10−51.46 × 10−42.70 × 10−58.77 × 10−54.80 × 10−57.36 × 10−4
Case 2MOPSO8.88 × 10−89.55 × 10−81.16 × 10−71.47 × 10−76.82 × 10−106.65 × 10−89.74 × 10−118.47 × 10−8
PSO7.34 × 10−68.39 × 10−51.00 × 10−44.88 × 10−44.64 × 10−53.23 × 10−41.05 × 10−42.98 × 10−3
Table 14. Generational distance (GD) and inverted generational distance (IGD) using MOPSO for Test system IEEE 39 system.
Table 14. Generational distance (GD) and inverted generational distance (IGD) using MOPSO for Test system IEEE 39 system.
MOPSOGDIGD
Test SystemCase 14.63 × 10−73.67 × 10−6
Case 26.39 × 10−94.86 × 10−6
IEEE 39Case 11.26 × 10−82.19 × 10−8
Case 23.82 × 10−87.16 × 10−9
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, N.-C.; Hsu, C.-W.; Sen, A. Parameter Tuning for WECC Generic Photovoltaic System Models Using Latin Hypercube Sampling and Pareto Optimality. Mathematics 2023, 11, 2739. https://doi.org/10.3390/math11122739

AMA Style

Yang N-C, Hsu C-W, Sen A. Parameter Tuning for WECC Generic Photovoltaic System Models Using Latin Hypercube Sampling and Pareto Optimality. Mathematics. 2023; 11(12):2739. https://doi.org/10.3390/math11122739

Chicago/Turabian Style

Yang, Nien-Che, Chun-Wei Hsu, and Abhilash Sen. 2023. "Parameter Tuning for WECC Generic Photovoltaic System Models Using Latin Hypercube Sampling and Pareto Optimality" Mathematics 11, no. 12: 2739. https://doi.org/10.3390/math11122739

APA Style

Yang, N. -C., Hsu, C. -W., & Sen, A. (2023). Parameter Tuning for WECC Generic Photovoltaic System Models Using Latin Hypercube Sampling and Pareto Optimality. Mathematics, 11(12), 2739. https://doi.org/10.3390/math11122739

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop