Improvement of the Iraqi Super Grid Performance Using HVDC/HVAC Links by the Integration of Large-Scale Renewable Energy Sources

: The increasing use of renewable energy sources (RES) and their integration into transmission systems requires extensive studies to get more beneﬁts from these sources. Stability analysis is one of the most crucial issues in these systems, and therefore needs to be studied in depth. In this paper, an analytical study was conducted for the stability assessment of the Iraqi super grid after adding a solar power plant with a capacity of 1000 MW using both HVAC and HVDC links for comparison purposes. Thus, the objective of this study was to determine the best way to integrate RES into the power grid while improving its performance. The system’s stability was evaluated for two types of faults: line and bus disturbances. The strength of the grid and its ability to absorb the newly added energy was also studied. The simulations were performed with the PSS/E software using the Newton-Raphson method for a 500 km transmission line length. The results showed that the HVDC link outperforms the HVAC link in improving the overall performance of the power system.


Introduction
Currently, the electric power system faces new challenges, such as combating air pollution and reducing dependence on fossil fuels. Renewable energy sources (RES), although unconventional, can provide promising solutions to these challenges, but the integration of RES into the regular power system planning and operation is still limited [1,2].
Renewable energy sources are critical components of contemporary systems and one of the answers for meeting rising demand. Renewable energy is included in the generating mix, making it a suitable solution to the issue of electricity scarcity [3]. The ever-increasing need for electrical energy has been satisfied mainly through burning fossil fuels and converting their heat to electrical power. Excessive exploitation of natural resources has raised worries about greenhouse gas emissions, pollution, and harmful environmental consequences. Additionally, it emphasizes the sustainability of the generation necessary to fulfil the world's energy demand in the following decades using non-renewable fossil fuels. Renewable energy sources such as the sun, wind, and tides answer these problems [4].
Recent expansion in the electric power market and increased integration of RES into the power system have created new network design and management problems. As a result, it is necessary to conduct thorough research and analysis of the influence of integrated renewable energy sources on the current power system [5][6][7].
As electricity demand increases, the power system has extended [8]. As a result, the continuous load increase and the power system expansion may affect voltage stability and make the reactive power unstable [9][10][11]. The HVAC power transmission line is one major components of the electric power system, whose primary function is to transport electric power from the production centre to a load centre (usually separated by long distances) with minimal losses [12][13][14]. According to Kirchhoff's law, the relationship between the injected current I into the network nodes and the bus voltage V is described by the admittance matrix Y as follows [28,29]: where Y ij is the line's admittance between the ith and jth buses, V i is the voltage on the ith bus, and I i is the current on the ith bus. The injected current at node i can be expressed as follows: The active and reactive powers injected at node i are rewritten as follows [29]: where δ i and θ ik are, respectively, the voltage and the admittance matrix phase angles.
The power flow problem consists of solving a system of non-linear equations to determine the voltage magnitude and phase angle at each bus i of the network. Thus, the power problem is formulated as follows: where P gi and Q gi are the generated active and reactive powers, respectively, and P ci and Q ci are the active and reactive loads, respectively.
The voltage magnitude and phase angle of each bus are computed as follows: where r is the number of iterations.

Short-Circuit Analysis
Short-circuit level (SCL) is an electrical fault that occurs when a large current flows down an alternate channel, resulting in severe equipment damage, power failure, physical injury, or death [10]. Short circuit analysis determines the system's ability to create a short circuit current, the electrodynamic resistance capability of the wire and equipment, and the magnitude of the short circuit current relative to the interrupting power of overcurrent protection devices. The effects of a short-circuit often depend on the type, duration, and location (near sparking or arcing locations) of the fault. Phase-to-ground faults account for 80% of faults; three-phase faults account for less than 5% of initial faults, and 15% are phase-to-phase faults [30].

•
Phase-to-ground fault This form of failure activates the zero-sequence impedance, Z 0 . Unless spinning equipment is involved, short-circuit current, I SC0 is less than three-phase fault current (lower zero-sequence impedance). Calculations for I SC0 may be necessary depending on the neutral system (earthing system configuration), assuming the established thresholds for zero-sequence (high voltage) or earth-fault (low voltage) safety devices [31].

•
Three-phase fault This fault manifests itself in all three phases; short-circuit current can be expressed as: where U refers to the phase-to-phase voltage that corresponds to the transformer's no-load voltage. It is approximately 3% to 5% higher than the on-load voltage across the terminals. Computing the short-circuit current involves the calculation of Z SC , which is equal to the sum of all impedances through which I SC flows from the generator to the point of failure. The impedances of the power sources and the lines are defined as the impedance per phase that is "Positive-sequence": where ∑ R and ∑ X are the total resistance and reactance of a sequence, respectively.

HVDC System: Modelling and Control
The prospects of a future HVDC grid are anticipated to meet the rising power demand and incorporate renewable energy sources into the grid [32][33][34][35]. An HVDC grid can provide the extra required power capacity, and its economy for long-distance transmission is ideal for integrating remote renewable sources into the grid.
Popular HVDC converter topologies are the line commutated converter (LCC), the voltage source converter (VSC), and the newly developed modular multilevel converter (VSC-MMC) [36][37][38][39]. Historically, LCCs have been employed as rectifiers and inverters in an HVDC system. MMCs, however, offer many advantages over LCCs, including a smaller footprint, low harmonic output, control of both real and reactive power, black start capability, and no chance of commutation failure [32,35,36,40]. These characteristics make MMCs excellent for use as inverters in urban load centres where compactness, fault recovery, and power quality are critical [34,35]. However, LCCs offer a higher current capacity and are more economical and technologically mature compared to MMCs [41]. The drawbacks of an LCC are typically not a concern at the rectifier due to the general availability of land at a remote location and proximity to generation [34,35]. The hybrid LCC-MMC system can be a practical solution for riding through DC side faults and multiterminal connections, which are foreseen in the future HVDC grid. This paper pays attention to the LCC-HVDC topology applied to the Iraqi super grid.

LCC-HVDC System
In Figure 2, a typical LCC-HVDC link between two AC systems is illustrated. It consists of two terminals linked by a direct current connection. Each terminal has converters, transformers, filters, reactive power equipment, a control station, and other components. This figure shows two direct current lines with an electrode in each converter station to supply the ground return path [42]. An LCC converter's fundamental structure is a three-phase full-wave bridge referred to as a "Graetz Bridge" or a six-pulse bridge. In practice, two of these Graetz bridges are connected in series to make a 12-pulse bridge, eliminating the need for harmonic filtering and producing twice the DC voltage. Harmonic reduction of DC and AC voltages is accomplished by introducing a 30-degree phase shift between the AC voltages provided at the upper and lower 6-pulse bridges. The 30-degree phase shift is achieved by connecting a Y-connected transformer to the top 6-pulse bridge's secondary and a delta-connected transformer to the lower 6-pulse bridge's secondary. The switching from one thyristor switch to another on the same DC pole is referred to as commutation. The converter is called 'line commutated' because it depends on the line voltage for commutation to occur. Due to the use of unidirectional thyristor valves, the LCC cannot operate with reverse currents [35].
The DC terminal voltage, V dc , is the difference between the positive and negative pole voltages, V dcp and V dcn , respectively. For rectifier operation, V dc may be mathematically described as follows: (8), V dc0 is the voltage at zero ignition angle α, R c is the commutation resistance, B denotes the number of six-pulse bridges, I dc is the DC line current, V LL is the line-to-line voltage of the AC grid, and T is the transformer ratio. The ignition angle is used to postpone valve commutation. The commutation resistance R c = 3 π × X c is equal to the voltage drop across the source side line inductance, L c .
For inverter operation, a comparable equation is defined as: where γ is the extinction angle.

HVDC System Control with LCC Converters
The control's purpose is to regulate the flow of electricity over the HVDC circuit. Adjusting the DC internal voltages V dc0 cos α for the rectifier and −V dc0 cos γ, the inverter allows for the adjustment of the voltage and current. The rectifier's ignition angle α, and the inverter extinction angle γ, can be quickly adjusted through controls to respond within 1 to 10 ms [35]. The tap changers may also be adjusted after the initial quick response from the ignition angle to ensure that α is brought back within a desirable range. The tap change response ranges from 5 s to 6 s per step [35]. Key factors to be satisfied by the controls are [35]: prevention of large fluctuations in DC, I dc under AC disturbances; operation with high power factors; operation near rated voltages; and prevention of commutation failure, which is a major concern in LCC inverters. The DC voltage is a function of the power factor: V dc ≈ −V dc0 cos ∅, where cos ∅ is the power factor. Operation with a high power factor is desirable: to increase rated power; minimize valve stress; reduce losses and current ratings of AC system equipment connected to the converter; reduce the reduction in AC voltage due to higher loading, and minimize reactive power supply costs [35]. A high power factor is achieved with low ignition/extinction angles, α for the rectifier and γ for the inverter. Typically, α is maintained within 10 • ≤ α ≤ 20 • to ensure a high power factor without limiting the control range [42]. Conventional HVDC characteristic control curves are depicted in Figure 3. The rectifier and inverter controls are visually shown in these curves under normal and AC fault conditions. The characteristic curves for the rectifier and inverter are taken concerning the DC terminal of the rectifier, i.e., the rectifier DC terminal voltage on the y-axis and the x-axis, V dc and current, I dc , respectively. The inverter voltage with respect to the rectifier DC terminal considers the voltage drop across the HVDC line: Consequently, the operating point of the system is the intersection of the rectifier and inverter characteristic curves. As seen in Figure 3, the rectifier runs under constant current control, whereas the inverter functions under constant extinction angle control (i.e., constant γ control). The abbreviations CIA, CEA, CC, CM, and VDCOL represent the constant ignition angle, constant extinction angle, constant current, current margin, and voltage-dependent current order limit, respectively.

Fault Handling Strategy
This section summarizes an LCC's technique for dealing with AC and DC problems.

AC Faults
The DC voltage is a function of the AC line-to-line voltage; hence, a decrease in the AC voltage caused by a fault results in a drop in the DC voltage. Figure 4 depicts the HVDC LCC characteristic curves under AC fault conditions [43]. As shown in this figure, the current is regulated by a voltage-dependent current order limit (VDCOL), which is used to limit the reactive power demand from the AC system during the fault and hence minimizes the danger of commutation failure. The current margin (CM), which is commonly set from 0.1 pu to 0.15 pu, is established to guarantee that the two characteristic curves overlap at a single point during decreased voltage operation. If the AC voltage remains low, tap changers will be changed to bring the system voltages and currents up to parity with the ignition angle, α, back to their normal operating range. Figure 4a,b shows the operating point during rectifier and inverter side faults. For rectifier side faults close to the converter station, there is a significant voltage drop, which would require the inverter station to operate with a large extinction angle, γ, resulting in large reactive power demand. As a result, the DC system is shut down under extremely low voltages. A voltage drop from 10% up to 15% may cause commutation failure [35,42]. If the voltage drop is significant, the converter valves are blocked until the fault is cleared and the voltage recovers to its normal value.

DC Pole-to-Ground Fault
At a DC pole-to-ground fault, the converter switches from a constant angle extinction command to constant current control. The converter operates at a reduced current order for the rectifier of the CM quantity. The reduced current restricts the reactive power demand of the AC system. The rectifier switches its operating mode from constant current control to constant ignition angle control, α. During the DC fault condition, a fault current equal to the CM is sustained. Otherwise, fast-acting line protection can be employed during a poleto-ground DC fault. The rectifier and inverter both operate to invert their DC terminal voltage. The fault current is seeking to exit the fault trail and enter the converters. However, due to the unidirectional diodes, the current rapidly decreases to zero. When it comes to incorporating renewable energy sources, system stability is a prevalent issue. The system strength, which indicates the system variables' sensitivity to different disturbances, affects the performance of various components in a power system. Short circuit ratio (SCR) is a metric that indicates the strength of a network bus compared to the rated power of a device. Frequently, the SCR is used to indicate the system's strength. When disturbances are introduced, a more robust system generates few fluctuations in voltage and other variables. The equipment's SCR is defined as its short circuit capacity ratio to its MW rating on the bus [44,45]. A robust system has an SCR of more than 3, whereas a poor or severely weak system has an SCR between 3 and 2 or less. The short circuit ratio is a critical parameter in power electronic applications such as the power converters that connect renewable energy sources (wind and solar energy) to conventional power grids. When used with insufficient AC systems, HVDC/FACTs devices based on current source converters exhibit problems. When HVDC/FACTs devices are coupled to weak AC systems, particular technologies are frequently needed (short circuit ratio of less than 3).

Stability Analysis of the HVDC/HVAC Power System in the Presence of Large-Scale
When the SCR is close to one, capacitor-commutated or voltage-source-based converters are used for HVDC applications. If such technologies are not accessible, specialized studies will be necessary to assess the impact and provide solutions for avoiding or mitigating harmful consequences. Low levels of SCR may result in difficulties such as high overvoltages, low-frequency resonances, and control system instability [46,47].
SCR is often described as the ratio of the short circuit capacity of the bus on which the device is placed to the MW rating of the device. According to this definition, SCR is expressed as follows: where S i denotes the short circuit capacity of the bus before the device's connection, and P dc denotes the device's MW rating. Table 1 summarizes the system strength assessment categorization. Table 1. Short circuit ratio and AC system strength.

Classification Short Circuit Ratio (SCR)
Strong system SCR > 3 Weak system 2 < SCR ≤ 3 Very weak system SCR < 2 4.1.2. Weighted Short Circuit Ratio (W SCR ) The weighted short circuit ratio (W SCR ) is a more suitable indicator, as it considers the influence of system interactions and provides a more accurate estimate of system strength. W SCR is given by: where S i denotes the short circuit capacity of bus i before connecting solar plant i, P dci denotes the MW rating of the solar plant i to be connected, N is the number of solar plants that are completely interconnected, and i denotes the solar plant index [45].

Stability Analysis
The term "power system stability" refers to studying the dynamics of the power system in the presence of failures. After being exposed to a disturbance, the stability of a power system refers to its capacity to resume regular or steady operation. When a specific disturbance occurs in the system, power system instability may be defined as the loss of synchronism (i.e., specific synchronous machines moving out of phase). Three distinct kinds of stability are affected: steady-state, transient, and dynamic stability [27,29].
System frequency and voltage levels are two critical parameters in alternating current power systems. Furthermore, most electric power produced in power systems is a threephase alternating current where generators operate synchronously. During the normal operation of a power system, the mechanical and electromagnetic torques acting on each generator must balance, and all synchronous generators must revolve at the same electrical speed in phase with the system frequency [48]. If the power system is subjected to severe disturbances such as short-circuit defects in transmission lines, it is possible to preserve generator synchronism. Transient stability is described in [48,49] as the system's capacity to sustain generator synchronization in the presence of substantial failures. Otherwise, a lack of transient stability might continually induce rotor-angle variations and self-exciting power oscillations, further damaging power system components such as generators and transmission lines or even resulting in power system blackout [48].
As a result, the physical nature and methodology for transient stability should be thoroughly investigated to ensure power system security, particularly when new technologies such as high voltage direct current (HVDC) and renewable energy are integrated into power systems and have their own dynamic characteristics for transient stability. Consequently, the dynamic implications of such new technologies on the transient stability of power systems should be thoroughly researched for the intent of power system security [48].
Within this context, this paper assesses the power system's stability and evaluates its strengths versus disturbances after integrating large-scale renewable energy. Small signal rotor angle stability, which refers to the power system's ability to retain its synchronism under disturbances, is the considered one. The instability appears in rotor oscillations with increased amplitude owing to the lack of damping torque and an aperiodic mode stemming from the lack of synchronizing torque. The study is conducted on the Iraqi super grid, where the dynamic system behaviour is investigated before and after integrating the solar power plant using HVDC and HVAC links. The impact of line and bus outages on the system stability is examined for different cases using the PSS/E software. The latter utilizes statistics channels called spread angle, largest angle, smallest angle, and average angle. The smallest angle can be defined as the trace for smaller rotor angles recorded when applying disturbance. The largest rotor angle represents the trace of a larger angle recorded through disturbance. In contrast, the average rotor angle is the average value of the rotor angle recorded when applying disturbance. Finally, the spread rotor angle is the deviation of the angle when applying disturbance with respect to the steady-state rotor angle. If any generator goes out of step in the simulation, the angle spread will increase. Angle spread is, therefore, a good indicator of rotor angle instability.
The reference of stability limits, and all issues noted above, is the steady state analysis of the system in load flow equilibrium, which means that the smallest deviation in spread angle, largest angle, smallest angle, and average angle led to a good stability response.

Iraqi Power System
The electrical power system of Iraq is divided into three operational subsystems: from the control point of view, the North, Middle, and South regions, in addition to the Kurdistan region network, which is separate from the Iraqi network. The HVAC Iraqi transmission system operates with a frequency of 50 Hz and a transmission voltage of 132 kV up to 400 kV. The power system under study is the Iraqi power grid (400 kV), which contains 41 buses, 25 of which are generating buses, and 72 transmission lines, as illustrated in Figure 5. The peak generation is 17,513 MW, and the peak loading is 24,679 MW, so the shortage is 7166 MW based on the Iraqi centre load on 5 July 2020 [50].
Considering the information provided by the Iraqi Ministry of Electricity database, this work suggests adding a new transmission line with 1000 MW supplied from a new solar plant (SP) connected to the secondary station named West Baghdad (BGW4). It has been assumed that this new solar plant, named SHMS, has a capacity of 1000 MW and is used to feed the loads on bus BGW4 through a 500 km transmission line. To evaluate the strength of the Iraqi super grid and its stability versus disturbances after the integration of large-scale renewable energy, the dynamic system behaviour was assessed before and after integrating a solar power plant using bipolar HVDC and double line circuit HVAC links using the PSSE simulation tool. The test system parameters used in this study are summarized in Table 2.

Power Flow and Short Circuit Levels Simulation Results
Our previous work [51] included the study of the impact of the renewable energy penetration on the power load flow and short circuit levels using double circuit lines for HVAC and bipolar HVDC links. Simulation results demonstrated that the HVDC system surpasses the HVAC in power losses and SCL minimization. Additionally, results have shown that, for the HVAC link case, the short-circuit levels increase after adding a solar power plant and decrease progressively as the distance between this generation station and the load centre increases. Moreover, it was noted that the double circuit line provides acceptable power loss and SCL results. When HVDC links are used, the results showed that SCL considerably decreased whether the distance between the power plant and the load centre increase or decreases. Furthermore, it was noted that the power losses and the SCLs for bipolar links are significantly less than those provided by HVAC links.

System Strengths Evaluation for Large-Scale Solar Plant Integration
Numerical simulations were carried out using PSS/E software for the base system configuration (without any newly added power), the short-circuit levels for three-phase faults at the BGW4 bus, and its connections and are presented in Table 3. Subsequently, a new large-scale solar power plant with various capacities was added to the base system at BGW4 (bus no. 16407) using an HVDC link. The short circuit ratio is evaluated as indicated in Table 4. The above results show that the system is strong until adding 3900 MW and operates without problems after the integration of a solar power plant using an HVDC link; beyond this value, the system will be weak. Therefore, it can be concluded that the integration of solar power plant is bounded by the system strengths, which are further affected and exhibit weak and very weak conditions when the solar power generation capacity reaches 4000 MW.
With the integration of a 3900 MW solar power capacity into the system at various points, the short circuit ratio is as shown in Table 5. Considering full interaction, by Equation (11) the actual weight of the short circuit ratio is: W SCR = 2.65.
The corresponding point of interconnection (POI) W SCR is just 2.65, although all SCRs at each point of interconnection exceed 3, which would be an acceptable characteristic of connection to a strong system. This implies that when these solar plants are combined, the system strengths are significantly reduced compared to when each solar plant oscillates individually.

Dynamic Simulation Results and Discussion
A two-terminal HVDC line was simulated using the CDC4T-DC line control model, as shown in Figure 6. The CDC4T model assumes that an instantaneous response of the DC system to disturbances is coming from adjacent grids. It uses the same algebra as the power flow model to update converter parameters (α, γ) at every integration step. It assumes that converters can maintain the steady-state control algorithm using the same setpoints as in power flow. There is a special algorithm for simulating blocking and by-passing HVDC that uses some threshold voltages on AC and DC sides. Setting up these thresholds requires experience because they are very sensitive to how strong the interconnection to the system is.
The parameters of HVDC system control are shown in Table 6. In this study, the impact of line and bus outages on the system stability has been simulated for three cases: 1.
First case: the base system; 2.
Second case: after the integration of the solar plant using an overhead AC doublecircuit line; 3.
Third case: after integrating the solar plant using a bipolar DC link for each pole.
• Line disturbance First, stability analysis of the system was carried out by applying a line fault at the point node named 19436HMD4-19437ANBP loaded 300MW; the fault is applied after the simulation initiation (at 1 s) and is cleared after 0.06s. Figure 7 shows the system response during the onset of the short circuit: the spread angle, the largest angle, the smallest angle, and the average angle of the system, respectively. It can be noted in Figure 7a that the overall spread angle of the system had much better results for case 3 than for the other cases. HVDC proved its ability to improve the total spread angle of the system and make it smaller than its value even before adding a solar power plant. As for Figure 7b, which represents the average angle of the system, it can be noted that there is a clear difference that case 3 achieved better results than case 2. Figure 7c represents the smallest angle in the system, where the HVDC also kept the angle value from falling compared to the other cases. Although Figure 7d represents the largest angle of the system, it was proved that case 3 behaves close to the system before adding a solar power plant, thus keeping the value of the largest angle almost unchanged.
The voltage response during the short circuit occurrence at different points in the power system was also studied for all test cases. Figure 8 illustrates the results of the test points at buses no. 16840-GQD11 and 16227-YSP13. These plots prove that the voltage at bus no. 16840-GQD11 is estimated to be 1.010 pu, 1.035 pu, and 1.020 pu for cases 1, 2, and 3 respectively (Figure 8a); these values falls to 0.8 pu during the fault. The voltage at bus no. 16227-YSP13 is estimated to be 1.02 pu for cases 1 and 3 (Figure 8b) during the pre-fault. This value falls to 1.01 pu during the fault. After that, the voltage shows some oscillations but always lies within the stability limits. However, for case 2, the voltage gap is slightly larger compared to other cases. The fail is not far from the other cases, and it is caused by the decreased system inertia due to HVDC characteristics (HVDC does not have inertia).

•
Bus disturbance For this fault scenario, the power system stability has been evaluated for a bus disturbance at buses 18806-ZBP6, 16225-YSP11, and 16226-YSP12 loaded 558, 286, 286 MW, respectively. The same fault conditions, time occurrence, and duration are maintained. Figure 9 illustrates the transient response of all systems during the short circuit occurrence, angle spread, largest angle, smallest angle, and average angle of the system, respectively. The simulation results show a similar response to the first fault scenario with a slight improvement in the HVDC that shows better oscillations damping of the generator angles. Figure 10 illustrates the transient voltage response during the bus disturbance scenario at system points no. 16840-GQD11, and no. 16227-YSP13. Figure 10b shows that the voltage drop for this fault scenario has reached 0.18 pu. However, the power system has regained its stability after removing the fault.

•
Line and bus disturbance In this scenario, both line and bus disturbances have been applied. The Iraqi super grid stability has been examined for a line outage at 19436HMD4-19437ANBP and bus disconnection at buses no. 18806-ZBP6, no. 16225-YSP11, and no. 16226-YSP12. For this fault scenario, a total generation of 1730 MW is disconnected. Figures 11 and 12 illustrate the transient response of the power system angles and voltages to this fault scenario.
The simulation results described in the figures above show that when the line and bus faults are applied simultaneously, the angle deviation of the power system is greater than for the case where each fault is applied separately. Nevertheless, the Iraqi super grid has proven its resilience to such a severe fault. Moreover, as for previous fault scenarios, the HVDC system has outperformed the HVAC system, showing a better resilience to faults, resulting in a small deviation of the angle and better damping of the oscillations. Furthermore, the transient voltage response depicted in Figure 12 shows a similar result as for the bus disturbance case, which demonstrates that this type of fault affects the system voltage more than the line outage. However, after the fault clears, the system voltage shows some oscillations but lies within the stability limits.

Conclusions
This paper discussed the impact of large-scale renewable energy integration with both HVAC and HVDC systems, considering stability analysis. Simulation results prove that despite the severity of the applied faults, the Iraqi super grid has demonstrated its resilience to such a fault. Additionally, the HVDC system has outperformed the HVAC system, showing a better resilience to faults, resulting in a small deviation of the angle and better damping of the oscillations. The dynamic system behaviour in the presence of HVDC lines has shown better performance and a more stable state for permanent faults than pure HVAC networks. Moreover, the HVDC lines have a better performance when comparing transient with permanent faults. It provides additional voltage support to congested transmission lines, leading to system stability and reliability improvement. However, the most significant benefit of adding an HVDC line in accidents is the system's ability to withstand a considerable loss of power without significantly impacting the network's stability, which was observed when studying the line and bus disturbance.