Latin Hypercube Sampling Method for Location Selection of Multi-Infeed HVDC System Terminal

: Owing to the stochastic states of power systems with large-scale renewable generation, the impact of high-voltage direct current (HVDC) systems on the stability of the power system should be examined in a probabilistic manner. A probabilistic small signal stability assessment methodology to select the best locations for multi-infeed high-voltage direct current systems in alternating current (AC) grids is proposed in this paper. The Latin hypercube sampling-based Monte Carlo simulation approach is taken to generate the stochastic operation scenarios of power systems with the consideration of several stochastic factors, i.e., load demand and power generation. The damping ratio of the critical oscillation modes and the controllability of power injection to oscillation modes are analyzed by the probabilistic small signal stability. A probabilistic index is proposed to select the best locations of high-voltage direct current systems for improving the damping of the oscillation modes. The proposed methodology is applied to an IEEE 39 bus system considering the stochastic load demand and power generation. The results of probabilistic small signal stability assessment and a time-domain simulation show that the installation of a high-voltage direct current system on the selected locations can e ﬀ ectively improve the system damping.


Introduction
With the progress of power systems, especially the extension of large-scale interconnected power systems, the inter-area oscillation issue is becoming a major problem. The weak connection between subsystems may cause oscillation issues of the interconnected power system, particularly during the initial period of the interconnected power system. High-voltage direct current (HVDC) transmission systems are a promising method to enhance the connection of the interconnected power systems. The integration of HVDC systems makes the structure of the receiving-side power systems more complex. Moreover, the large-scale renewable generation can be integrated by HVDC systems, which increases the uncertainties of the system operation states. Thus, the location of HVDC, especially multi-infeed HVDC systems, on the system security and stability should be examined carefully [1,2].
Different methods have been proposed to investigate the best placement of HVDC converters. Most of the methods focus on the optimal power flow with the integration of HVDC power systems. In [3], an optimization algorithm is proposed to identify the HVDC line placement for the relief of congestion of transmission lines and saving the operational cost of the whole system. The trade-off between the minimization of operational cost and the controllability of the power flow is discussed in [4]. A comprehensive algorithm based on the flexibility matrix is proposed to determine the location of HVDC converters for balancing the cost and flexibility. In [5], the N-1 criterion is considered in the location selection of HVDC systems in meshed power systems to reduce the congestion. The problem is transferred to a mixed integer linear programming optimization problem and an adaptive approach is employed to improve the flexibility of the proposed method in large-scale power systems. With the penetration of renewable energy in power systems, the uncertainty of the system operational point is increased, which should be taken into account in the placement of HVDC. In [6], a double stage stochastic optimal algorithm is proposed to deal with uncertainties of load demand and power generation in a HVDC placement problem. The proposed method can find the best HVDC locations for minimizing the load of transmission lines without requiring many parameters. In [7,8], the multi-infeed short-circuit ratio (MISCR) defined by International Council on Large Electric systems (CIGRE) is used to identify the appropriate location of HVDC systems to improve the security of hybrid alternating current/direct current (AC/DC) power systems. However, the aforementioned papers did not consider the impact of HVDC locations on system stability, e.g., the oscillation issues.
Small signal stability analysis (SSSA) is a promising method to investigate the system stability. The conventional SSSA is based on the linearized model and the eigenvalue analysis method of the system is performed at a certain operating point, which belongs to deterministic stability analysis. To examine the stability of a system under different operating conditions, a great number of deterministic states need to be analyzed [9]. The main problem of the deterministic method is that it cannot objectively reflect the inherent uncertainty nature of the system's network structure, power-generation mode, load level and element parameters. Therefore, it is difficult to comprehensively and accurately evaluate the overall system stability [10].
In 1978, Burchett applied the probability method to power system probabilistic small signal stability analysis (PSSSA) for the first time and analyzed the influence of uncertainty of system parameters following normal distribution on the stability of power systems in the case of small disturbance [11]. PSSSA refers to considering the uncertainty of input random variables, conducting SSSA on the power system, and using the probabilistic analysis method to obtain the probabilistic characteristics of the eigenvalues of the electromechanical mode. A quasi-Monte Carlo-based PSSSA method to evaluate the dynamic performance of electric vehicles and energy conversion systems is presented in [12]. A PSSSA method based on Monte Carlo (MC) approach for iterative assessment is proposed in [13] and the PSSSA considers the uncertainties of the operational states by modal analysis of small signal stability. Cumulative probability density function of damping ratios of oscillatory modes (OMs) is used as the evaluation index.
MC simulation is a statistical test method of parameter estimation when the sample size is large. MC has been widely used in power systems such as reliability assessment, and the number of tests has nothing to do with the scale of the system. MC is used to produce a great deal of random computational scenarios in accordance with the distribution density of the random sources [14]. Recent literature has used MC simulation for SSSA of power systems with small sets of random input variables [13,15]. Latin hypercube sampling (LHS) was first proposed by McKay et al. [16] in 1979. It is an efficient stratified sampling technique and it can obtain an accurate result with a much smaller sample size than simple random sampling (SRS) [17].
Owing to the multiple stochastic factors of interconnected power systems, e.g., load demand and power generation, the damping ratio of the OMs may vary with the change of operating states. The Latin hypercube sampling-based Monte Carlo simulation (LHS-MCS) can be used to identify the critical OMs and the locations for the DC system terminal. This paper proposes a PSSSA methodology to evaluate the best location for multi-infeed HVDC systems. The LHS-MCS approach is adopted to generate the stochastic operation scenarios of power systems with the consideration of several stochastic factors, i.e., load demand and power generation. The damping ratio of the critical OMs and the controllability of power injection to OMs are analyzed by the probabilistic small signal stability. A probabilistic index is proposed to select the best locations of HVDC systems for improving the damping of the OMs. The proposed methodology is tested on an IEEE 39 bus system by selecting the stochastic load demand and power generation. The results of PSSSA and the time-domain simulation show that the installation of the HVDC system on the selected locations can effectively improve the system damping.
The organization of this paper is as follows. Section 2 describes the theoretical background of uncertainty analysis, MC simulation, and the LHS method. Section 3 presents the proposed methodology and discusses the main issues associated with the conditions for application of LHS method in PSSSA and the calculation of probabilistic indices for selecting the best locations of multi-infeed HVDC systems. In Section 4, the selection results of the multi-infeed HVDC systems obtained from a New England 39 bus system are shown. In Section 5, comparisons and discussions with the results of other papers are given. Finally, Section 6 concludes the paper.

Theoretical Background
An HVDC system is a promising technique for renewable energy integration. With the increasing penetration of renewable energy integrated by HVDC systems, the operational point of whole power systems is time varying. The stability of hybrid AC/DC power systems could be influenced by the uncertainty originating from the load demand and power generation. In the planning stage of an HVDC system, the influence of HVDC station locations on the system stability must be investigated in a probabilistic way. To analyse the probabilistic stability of hybrid AC/DC power systems and evaluate the candidate locations of HVDC stations, the Monte Carlo method is employed to obtain the random operational scenarios and Latin Hypercube sampling method is used to improve the calculation efficiency.

Uncertainty Analysis Method
The general physical process of a system can be represented by establishing one physical model. This physical model can usually use an abstract function to represent input and output, and there is a certain mapping relationship, the function is described as: where the vector h is a set of functions describing the behavior of the model, such as the dynamic behavior of power system; q = [q 1 , q 2 , . . . ] T is a vector of input variables; and p = [p 1 , p 2 , . . . ] T is a vector of output variables. From the perspective of probability theory, uncertainty analysis includes determination of the probability density function (PDF) of the output variable p generated by the function h and the PDFs defining the probability space of the input variable q. In addition, the PDF of the output variable p can be represented as a cumulative probability density function (CPDF). Both of the functions summarize the uncertainty of output variables, and yet CPDF provides more illustrative abstracts in sample-based studies. CPDF is formally defined by integration [13] prob where prob(p i ≤ P) is the probability that an output p i under P will happen, PDF(q) is the PDF of input variables, the differential dq is the integration domain, and: Energies 2020, 13, 1646 4 of 15 In fact, the integral (2) is obtained by some approximate procedure. Although there are some new methods have been proposed in other literature to Equation (2), the effectiveness of MC simulation in uncertainty evaluation has been proved in many fields [18,19].

Monte Carlo Simulation
MC simulation aims at estimating the probability of an event by a lot of experiments. The basic thought is firstly to establish a probabilistic model for the problem. The next step is to calculate the statistical characteristics of the problem through a large number of sampling experiments. The final step is to get all kinds of estimators to find the approximate answer.
The application of MC in uncertainty analysis is justified by two statistical theories: (1) the Law of Large Numbers, and (2) the Central Limit Theorem.
Let X 1 , X 2 , . . . X k be a sequence of independent and identically distributed random variables with expectation E(X k ) = µ (k = 1, 2, . . . n). The arithmetic mean of n variables is calculated by: For any small value ε > 0, the Law of Large Numbers states that: Equation (4) indicates that the sampling average of X k approaches converge in probability towards the expected value µ when n is large enough.
The distribution of the sample average can be described by the Central Limit Theorem. Assume that the variance of X k is σ 2 (0 < σ 2 < ∞). For any real number x, the distribution function F n (x) of sample average is calculated by: where φ(x) is the standard normal distribution. When n is large enough, the average of the samples is similar to the normal distribution whose expectation is µ and the variance is σ 2 .

Latin Hypercube Sampling
LHS is an effective probabilistic sampling method which can be used to generate the necessary scenarios for PSSSA. X 1 , X 2 , . . . , X M is the M input random variables. Each random variable obeys a certain probability distribution: where Y m is the probability of X m and is scaled from 0 to 1, and F m is the probability cumulative function of X m . The cumulative curve is divided into N equal intervals and only one sampling value should be selected randomly in each interval, which is illustrated in Figure 1.
Then the nth sample value of the mth random variable x mn is calculated by the inverse function of F m: where a is an internal variable which meets the 0-1 distribution. The sample values of X m can be written in a column vector as [x m1 x m2 . . . x mN ] T . All the sample values of input variables construct a N*M sampling matrix. Then each column of the sampling matrix should be shuffled randomly, and each row constitute a random scenario [20]. LHS can effectively improve the computation efficiency of MC as the sampling number are reduced significantly without the loss of accuracy.

Probabilistic Small Signal Stability Analysis
PSSSA can evaluate the dynamic performance of power systems considering uncertainty operation states. The stochastic operation scenario of power systems are generated firstly, and then the SSSA is used for each of the random scenarios. The probabilistic dynamic characteristics of power systems can be acquired by the probability statistics method, such as distribution of eigenvalues, expectation and variance of damping ratio and oscillation frequency, etc.
The linearized model and the random variables of power systems are the bases of PSSSA. The state matrix A of large-scale complex power systems is obtained by DIgSILENT/PowerFactory, and the input matrix B and output matrix C is obtained by a numerical calculation method introduced in [21].
Several stochastic variables can affect the operation states of power systems, which leads to different distributions of eigenvalues. Random scenarios are composed of different combinations of stochastic variables. In this paper, the following random variables are considered: (1) the output power of synchronous generators; (2) the load demand. Both stochastic variables obey the normal distribution. Then the oscillation characteristics of power systems are analyzed by the statistical method and the location which has the best controllable performance will be distinguished.

Modal Analysis
The dynamic behavior of one power system is often represented by a series of non-linear ordinary differential equations [22]: ,..., , where x is the state variables; xi is a element of the state vector; ui is the system input variables, which are external signals that can influence the behavior of the system;  x is the derivative of the state variable with respect to time; n is the system order, r is the number of inputs to the system; and t is time. Equation (9) can be rewritten as follows:  LHS can effectively improve the computation efficiency of MC as the sampling number are reduced significantly without the loss of accuracy.

Probabilistic Small Signal Stability Analysis
PSSSA can evaluate the dynamic performance of power systems considering uncertainty operation states. The stochastic operation scenario of power systems are generated firstly, and then the SSSA is used for each of the random scenarios. The probabilistic dynamic characteristics of power systems can be acquired by the probability statistics method, such as distribution of eigenvalues, expectation and variance of damping ratio and oscillation frequency, etc.
The linearized model and the random variables of power systems are the bases of PSSSA. The state matrix A of large-scale complex power systems is obtained by DIgSILENT/PowerFactory, and the input matrix B and output matrix C is obtained by a numerical calculation method introduced in [21].
Several stochastic variables can affect the operation states of power systems, which leads to different distributions of eigenvalues. Random scenarios are composed of different combinations of stochastic variables. In this paper, the following random variables are considered: (1) the output power of synchronous generators; (2) the load demand. Both stochastic variables obey the normal distribution. Then the oscillation characteristics of power systems are analyzed by the statistical method and the location which has the best controllable performance will be distinguished.

Modal Analysis
The dynamic behavior of one power system is often represented by a series of non-linear ordinary differential equations [22]: where x is the state variables; x i is a element of the state vector; u i is the system input variables, which are external signals that can influence the behavior of the system; .
x is the derivative of the state variable with respect to time; n is the system order, r is the number of inputs to the system; and t is time. Equation (9) can be rewritten as follows: . If the derivative of the state variable is not an explicit function of time, the system is autonomous [18]. At this point, Equation (10) can be simplified as: .
Energies 2020, 13, 1646 6 of 15 Output variables are the ones which can be directly measured in power systems. They can be described by functions of state and input variables: where m is the number of the output variables; y = [y 1 , y 2 , . . . y m ] T is the vector of output variables; and g = [g 1 , g 2 , . . . g m ] T is the vector of non-linear functions which define the outputs.
In small signal analysis, a linear technique can give useful information for the intrinsic dynamic characteristics of power systems. The linearized model of Equations (11) and (12) is [23]: where A is an n*n state matrix; B is an n*r input matrix; C is an m*n output matrix; D is a feedforward matrix. The oscillation modes can be identified by the eigenvalues of state matrix A [24]. Lyapunov's first law indicates that if all eigenvalues have negative real parts, the system is asymptotically stable at equilibrium point. The real eigenvalue λ i = α i corresponds to a non-oscillating mode. Negative real eigenvalues (α i < 0) indicate a decay mode. The larger the amplitude, the faster the decay. Positive eigenvalues (α i > 0) indicate non-periodic instable modes. The complex eigenvalues appear as conjugation pairs λ i = α i + jβ i , and each pair corresponding to one OM. The real part is related to the damping ratio ς i , which is defined as [25]: The negative real part means damped oscillation, and the positive real part means instable OM [24]. The decay rate of the oscillation amplitude is determined by the damping ratio ς i . The time constant of amplitude decay is 1/|α i |. Because ς i embodies the changes of the real and imaginary parts of the eigenvalues, it is considered a more appropriate measurement of the degree of the oscillation damping than α i which only considers the changes of the real parts. The mode corresponding to each critical eigenvalue which satisfies the requirements can represent different types of oscillation. The imaginary part β i is related to the oscillation frequency of the OMs.

The Index of Best Buses Selection
After performing modal analysis for all the random scenarios, the oscillation characteristics of OMs can be obtained for further analysis. The probabilistic critical OMs (PCOM) should be identified first. PCOM is defined as the OMs whose probability of being poorly damped exceeds 50%: where S pcom represents the set of PCOMs, N om is the number of OMs, P(·) means probability, which can be calculated by: where R im represents the damping state of ith OM in mth scenario. Then the probabilistic controllability indices of all the potential buses can be calculated. The controllability index reflects the influence of input variables on OMs. Since the main aim of the AC/DC interconnection system systems is power transmission, the element of input vector u is the injection power at different buses. The controllability index is defined as: where CI ob is the magnitude of the products of v o and B b , v o represents the oth left eigenvector, and B b is the bth column of input matrix B. Superscript s represents the s-th operation scenario. CI ob evaluates the controllability of the power injected at bus b for oth mode. CI o evaluates the controllability of the power injected at all buses for oth mode. CI is updated as [26]: where max CI o s is the maximum CI ob out of all locations on the s-th operation scenario.
Then average value µ ob and variance σ 2 ob can be calculated by: The best locations can be selected as the ones which has the highest µ ob and σ 2 ob < σ 2 max .

The Process of Probabilistic Small Signal Stability Analysis (PSSSA) and Location Selection
The flow chart of proposed method is shown in Figure 2. The detailed processes are as follows: (1) Set the total sample number N; (2) Generate random load and generator's output by the LHS method; (3) The DC system is equivalent to a −100 MW load. The −100 MW load is regarded as a mobile load. Put the mobile load and a breaker to one bus in once simulation;

Case Study
The proposed method is applied in the IEEE 39-bus system. The result of modal analysis is obtained by the Modal Analysis Package of DIGSILENT/PowerFactory. The random load of the system is assumed to follow the normal distribution. The expected value and the variance of loads are Pl0,i and 0.1Pl0,i. The output of generator changes according to the load level; 1000 random operating scenarios are generated by the LHS method. For the test system, there are 9 poorly damped low frequency oscillation (LFO) modes. The probability of these modes is shown in Table 1. A box plot of their damping ratio is shown in Figure  3.

Case Study
The proposed method is applied in the IEEE 39-bus system. The result of modal analysis is obtained by the Modal Analysis Package of DIGSILENT/PowerFactory. The random load of the system is assumed to follow the normal distribution. The expected value and the variance of loads are P l0,i and 0.1P l0,i . The output of generator changes according to the load level; 1000 random operating scenarios are generated by the LHS method. For the test system, there are 9 poorly damped low frequency oscillation (LFO) modes. The probability of these modes is shown in Table 1. A box plot of their damping ratio is shown in Figure 3.  According to Table 1, only OM 6 is considered to be the probabilistic critical OM. Other OMs whose damping ratio is most likely lower than 0.1 can be included in Spcom. Thus, Spcom = {OM 3, OM 4, OM 5, OM 6, OM 7, OM8}. The mean values of the mode controllability for Spcom are shown in Table 2. It can be observed from Table 2 that the top three buses are bus 39, bus 1 and bus 9 which are sorted by the mean value of controllability index and the probability of 6 modes. The variance of these three buses is not the maximum among all buses. Thus, the multi-infeed HVDC system is considered to be connected to the bus 39, bus 1 and bus 9. The box plots of controllability index results from modal analysis on critical modes are shown in Figure 4. A total of 39 buses were observed. According to Table 1, only OM 6 is considered to be the probabilistic critical OM. Other OMs whose damping ratio is most likely lower than 0.1 can be included in S pcom . Thus, S pcom = {OM 3, OM 4, OM 5, OM 6, OM 7, OM 8}. The mean values of the mode controllability for S pcom are shown in Table 2. It can be observed from Table 2 that the top three buses are bus39, bus1 and bus9 which are sorted by the mean value of controllability index and the probability of 6 modes. The variance of these three buses is not the maximum among all buses. Thus, the multi-infeed HVDC system is considered to be connected to the bus39, bus1 and bus9. The box plots of controllability index results from modal analysis on critical modes are shown in Figure 4. A total of 39 buses were observed.  From Figure 4, the box plot area of bus 39, bus 1 and bus 9 is smaller than other buses, which indicates that the controllability index distribution of the three buses is relatively concentrated. The value of these three buses' controllability index are closer to 1 than other buses. It can be seen the From Figure 4, the box plot area of bus39, bus1 and bus9 is smaller than other buses, which indicates that the controllability index distribution of the three buses is relatively concentrated. The value of these three buses' controllability index are closer to 1 than other buses. It can be seen the result of box plots and mean value table are same. Thus, Bus39, bus1 and bus9 are the selected location.
To demonstrate the damping performance of the HVDC system, a three-phase fault is caused on the transmission line 14-15 of the 39-bus system at t = 1 s and lasts for 0.1 s. The simulation results are shown in Figure 5. Case 1: one HVDC system connects to bus39; Case 2: two HVDC systems connect to bus39 and bus1; Case 3: three HVDC systems connect to bus39, bus1 and bus39. To demonstrate the damping performance of the HVDC system, a three-phase fault is caused on the transmission line 14-15 of the 39-bus system at t = 1 s and lasts for 0.1 s. The simulation results are shown in Figure 5. Case 1: one HVDC system connects to bus 39; Case 2: two HVDC systems connect to bus 39 and bus 1; Case 3: three HVDC systems connect to bus 39, bus 1 and bus 39. It can be observed from Figure 5 that the HVDC system has the effective damping performance. When there is no HVDC to be connected, the oscillations last for long time. When bus 39 connects with one HVDC system, the oscillations time decreases significantly. When the candidate three buses all connect with HVDC systems, it takes less time to return to the equilibrium point. The dynamic response of the rotor angular velocity can be shown by Figure 6. It can be observed from Figure 5 that the HVDC system has the effective damping performance. When there is no HVDC to be connected, the oscillations last for long time. When bus39 connects with one HVDC system, the oscillations time decreases significantly. When the candidate three buses all connect with HVDC systems, it takes less time to return to the equilibrium point. The dynamic response of the rotor angular velocity can be shown by Figure 6. Figure 6 shows that the system connected with HVDC system can have better damping effect than there is no HVDC to be connected. When the three buses all connect with HVDC systems, the system has the best damping effect.  Figure 6 shows that the system connected with HVDC system can have better damping effect than there is no HVDC to be connected. When the three buses all connect with HVDC systems, the system has the best damping effect.

Discussion
By contrast with the previous location selection method of HVDC systems, the proposed method mainly focusses on the impact of the locations on system stability. In [5,6], the best locations of HVDC links in 39-bus systems are identified as {bus3, bus22, bus35} and {bus8, bus20, bus27}, respectively. To compare the system dynamic responses with different HVDC locations, three cases are set: 1) Case A: the locations selected by this paper {bus1, bus9, bus39}; (2) Case B: the locations found by [5], i.e., {bus3, bus22, bus35}; (3) Case C: the locations found by [6], i.e., {bus8, bus20, bus27}. A three-phase fault is applied on the transmission line L14-15 at t = 1 s and lasts for 0.1 s. The system dynamics responses are shown in Figure 7. Since the power flow with different HVDC locations is different, only the power deviation of the transmission lines after the disturbance are shown. It can be observed that the locations identified by this paper, i.e., Case A, has the best damping performance. The locations in [5,6] has little influence on the power oscillations as the aim of these locations are cost reduction and power flow controllability.
The grid losses are summarized in Table 3. One can observe that the grid losses with the HVDC locations identified in [5,6] are lower than with the locations selected in this paper. Although the grid losses are increased, the stability of the hybrid AC/DC power system is improved.

Discussion
By contrast with the previous location selection method of HVDC systems, the proposed method mainly focusses on the impact of the locations on system stability. In [5,6], the best locations of HVDC links in 39-bus systems are identified as {bus3, bus22, bus35} and {bus8, bus20, bus27}, respectively. To compare the system dynamic responses with different HVDC locations, three cases are set: (1) Case A: the locations selected by this paper {bus1, bus9, bus39}; (2) Case B: the locations found by [5], i.e., {bus3, bus22, bus35}; (3) Case C: the locations found by [6], i.e., {bus8, bus20, bus27}. A three-phase fault is applied on the transmission line L14-15 at t = 1 s and lasts for 0.1 s. The system dynamics responses are shown in Figure 7. Since the power flow with different HVDC locations is different, only the power deviation of the transmission lines after the disturbance are shown. It can be observed that the locations identified by this paper, i.e., Case A, has the best damping performance. The locations in [5,6] has little influence on the power oscillations as the aim of these locations are cost reduction and power flow controllability.

Conclusions
In this paper, a probabilistic small signal stability assessment methodology based on the Latin hypercube sampling-based Monte Carlo simulation approach has been proposed. The heavy computation burden of the traditional MC method is significantly decreased by the LHS technique. The influence of potential locations of HVDC systems on the system stability is studied in a probabilistic way. The probabilistic critical oscillatory modes are identified and the candidates of the potential locations of multi-infeed HVDC systems can be distinguished by the probabilistic controllability index to improve the system damping performance. This method can be used for the initial selection of the scheme of arbitrary locations of the HVDC station in the planning stage. The effectiveness of the selected candidate locations is validated by the non-linear time-domain simulations in DIgSILENT/PowerFactory. Compared with other methods, the locations identified by the proposed method can effectively improve the stability of the hybrid AC/DC power system. The results show that the system damping to the probabilistic critical modes are enhanced and the power oscillations on the transmission lines are suppressed effectively. However, the grid losses are slightly magnified. Thus, in our future work, both the economic factors and the stability factors will be considered in the HVDC placement tasks.   The grid losses are summarized in Table 3. One can observe that the grid losses with the HVDC locations identified in [5,6] are lower than with the locations selected in this paper. Although the grid losses are increased, the stability of the hybrid AC/DC power system is improved.

Conclusions
In this paper, a probabilistic small signal stability assessment methodology based on the Latin hypercube sampling-based Monte Carlo simulation approach has been proposed. The heavy computation burden of the traditional MC method is significantly decreased by the LHS technique. The influence of potential locations of HVDC systems on the system stability is studied in a probabilistic way. The probabilistic critical oscillatory modes are identified and the candidates of the potential locations of multi-infeed HVDC systems can be distinguished by the probabilistic controllability index to improve the system damping performance. This method can be used for the initial selection of the scheme of arbitrary locations of the HVDC station in the planning stage. The effectiveness of the selected candidate locations is validated by the non-linear time-domain simulations in DIgSILENT/PowerFactory. Compared with other methods, the locations identified by the proposed method can effectively improve the stability of the hybrid AC/DC power system. The results show that the system damping to the probabilistic critical modes are enhanced and the power oscillations on the transmission lines are suppressed effectively. However, the grid losses are slightly magnified. Thus, in our future work, both the economic factors and the stability factors will be considered in the HVDC placement tasks.