Optimal Connection of Offshore Wind Farm with Maximization of Wind Capacity to Power Systems considering Losses and Security Constraints

)e technical, economic, and environmental constraints related to the construction of new transmission lines are complex issues related to the definition of points for connecting new offshore wind farms (OWFs) to the grid. In this context, it has become an important research topic to choose the best OWF connection point to a power system, among some geographically close to each other within a given region, aiming at ensuring maximum generation capacity of the wind farm and safe use of existing transmission network.)e objective of this work is to present a methodology to determine the optimal OWF connection point in a power system, with maximum penetration of firm wind power and minimum loss, considering security constraints related to the “N− 1” contingency criterion, exchange limits between areas, and a strategy to reduce the number of constraints in the optimization problem. )e algorithm is modeled using a Mixed Integer Nonlinear Programming (MINLP), and it is evaluated in a tutorial system and three well-known other networks from literature: IEEE 14-Bus, IEEE RTS-79, and Southern Brazilian System.


Introduction
Renewable energy is receiving more attention due to growing concerns about the environment and the depletion of conventional energy sources based primarily on fossil fuels.In this sense, wind energy is one of the most prominent options in renewable sources investment scenarios.Moreover, many countries are planning ambitious scenarios with this energy for the future [1,2].
Onshore wind energy is one of the cheapest types of renewable options.It is a good choice for investors in many countries due to its relatively short installation time, government incentives, and easy operation procedures.However, in most cases, regions with high wind generation potential are far from the load centers.In some cases, it justifies the installation of wind farms in a marine environment, which is known as offshore wind farms (OWFs).On the other hand, this option needs a suitable transmission system to inject the wind power into the existing network.
e problem is to determine the best network point for connecting these farms to the main onshore grid, considering that they are installed in places with high energy potentials.Building a new structure for OWF requires not only investment but also regulatory issues and environmental barriers.Also, for remote wind farms, the produced wind energy may not be fully usable due to constraints of transmission lines as, for instance, their thermal limits.
Optimal strategies for connecting wind farms to the grid should be studied to overcome all the mentioned barriers.However, a comprehensive analysis is required over an extended period to find the ideal connection points to the main grid.e aforementioned analysis should consider the load pro le, existing energy sources, transmission lines' infrastructure, and wind speed scenario of all potential regions under study.
is paper proposes a novel methodology to determine the optimal connection point of a transmission system for injecting the power from OWF with the purpose of maximizing the wind energy penetration with minimum network loss, among some candidate buses "k" geographically close to each other within a given region, connected or not by a "n" possibilities of switches "sw k ," as illustrated by Figure 1. e proposed approach considers important aspects and constraints related to the problem and uses Mixed Integer Nonlinear Programming (MINLP) for the decisions about the connection points within a given region.e approach is assessed through studies using three well-known networks besides a tutorial case.

Contributions of is Paper.
e main contributions of the present work can be summarized as follows: (i) A novel and e cient approach based on MINLP to determine the best network connection point for injecting OWF power that maximizes the wind energy penetration and minimizes the transmission losses (ii) e proposal of a strategy for representing the transmission line constraints in the optimization problem, which allows the e cient incorporation of the "N − 1" contingency criterion to the problem for security analysis

Paper Organization.
e research work is structured as follows.Section 2 shows a brief overview of the literature about the themes related to this paper.Section 3 presents the proposed approach explaining in details the strategies used in each step of the optimization process.Section 4 gives a tutorial case for explaining the steps of the proposed algorithm, whereas other results for practical systems are shown in Section 5. Finally, Section 6 gives the main conclusions.

Background and Related Works
In this section, a brief overview is made on the maximum penetration of wind power in electric power systems (EPSs) (Section 1.1) and optimal allocation and connection of wind power generation (Section 1.2).Besides, it is highlighted how this paper relates to these issues.

Maximum Penetration of Wind Power in Electric Power Systems (EPS).
An optimization cost function is de ned in [3] to minimize an economic cost criterion for wind turbine farm, comprising the choice of the best wind resource regions.A speci c target related to wind energy penetration is initially assumed to be feasible for a rm wind capacity connection.However, it would not be realistic in most of the energy systems with current transmission limitations.As an improvement of [3], the problem described in [4] attempts to maximize the wind power potential of an existing system in the short term, before the transmission expansion in the long run.Incremental penetration targets of rm wind energy are applied from a lower initial level to a threshold when the optimization model becomes unfeasible.
Reference [5] analyzes the possibility of considering the wind as a non rm power source, where historical wind data are used to model wind energy variations and spatial interdependence.
e Benders decomposition method was applied to explore the structure of the large-scale constraint diagonal matrix.e wind and load demand pro les, as well as fuel price and contingency management sensitivities, are also evaluated for demonstrating the importance of an ideal wind power capacity allocation solution.
In [3][4][5], the "N − 1" contingency criterion is considered and the thermal limits of the existing transmission lines and transformers are observed.ese aspects are also addressed in the approach proposed in the present paper.Moreover, this paper proposes the inclusion of transmission losses in the economic dispatch model and an e cient strategy to include only the active constraints into the optimization problem.

Optimal Allocation and Connection of Wind Power
Generation.As previously discussed, wind power is site dependent, and in many cases, good wind locations are far from the grid.Additionally, the connection of wind power to a power system can be hampered by transmission capacity limitations.In this scope, reference [6] presents a methodology for nding optimal positions in an existing transmission system to connect rm wind capacity for achieving desired renewable and secure penetration targets.
Following this line, paper [7] addresses the problem of optimal allocation of wind capacity and determines the capacity of the transmission line to connect the wind farm to the grid.In [8], a network reinforcement model incorporating the optimal integration of wind farms (WFs) into power systems is proposed.Based on the reinforcement Considering the diverse interests of investors and power systems' operators in wind farms, Zeńczak [9] presents issues related to the choice of the best location for connecting wind farms.In [2], it is presented a robust planning model to determine the optimal allocation of wind farms in a multiarea power system, so that the expected energy not supplied (EENS) is minimized under uncertain wind power and unit forced outages.
In this scenario, a MINLP is introduced in the present paper to handle the nonlinear nature of the problem and the integer binary variables that are proposed to represent switches' statuses, which are associated to the options for connecting OWF to a candidate bus of a power system.

Proposed Methodology
It is required a feasibility technical study to verify if the wind energy connection does not congest the transmission network after possible network contingencies.Note that a rm wind connection does not observe wind reduction due to system limitations.However, when a connection is non rm, the system operator can reduce the wind in case of contigencies [7].As mentioned before, the present work considers a rm wind power connection.
In this sense, it is assumed initially that the transmission system is decongested (not violating the established limits) and that the addition of new wind generating capacity must preserve this situation.us, a low target for wind energy penetration is de ned to assess the system feasibility.As the optimization problem gives a feasible solution for this condition, the wind energy is increased through small discrete steps until the connection of additional wind power to any of the candidate buses is not possible without congestion in the system.e ow chart of the proposed algorithm, based on [4], is presented in Figure 2.
In Step 1, the initial wind generation time series are generated, which are considered as being from historical data for the candidate regions.ese series are used in Step 2 for the economic dispatch of conventional and dispatchable generation in the system, aiming at minimum losses.In Step 3, a direct current power ow (DCPF) model [10] uses these wind power series to obtain a simpli ed power ow solution that allows identifying the lines that have their active power capacity violated.en, the violated lines de ne linear inequality constraints related to transmission limits for a MINLP problem in Step 4. e purpose of Step 4 is to determine the connection point in the power system for the o shore wind farm, aiming at maximizing the wind power utilization through a feasible connection solution for the current wind power penetration level.In this step, the decision variables of the optimization problem are the wind power capacity (C k ) for each candidate bus k.If the MINLP nds a feasible solution, the wind power time series are updated in Step 5 as described in the next section.Notice that the updated series cannot t the corresponding series that were introduced in the economic dispatch.In this case, Steps 2 to 5 are carried out again until the power wind time series of Step 5 are close to the input series of Step 2. e di erence between these series is further described in the next section.When it occurs, the convergence criterion is achieved.
After obtaining the convergence of the process that comprises Steps 2 to 5, if the solution is feasible, an evaluation is made to check the presence of violated transmission lines in Step 6.If a violation is identi ed, the aforementioned iterative process is repeated until no violation occurs.Otherwise, the wind energy penetration is increased in Step 7, and the whole process is repeated until a no feasible wind penetration level is achieved.At this infeasible point, the previous point is identified as that presents the optimal connection of the OWF, associated with maximum bearable wind power capacity for the system under analysis.
In the following topics, each step of the methodology proposed in this paper is discussed in detail.

Approximate System Wind Power Time Series (Step 1).
In order to carry out the economic dispatch of the conventional generations, it is necessary to have information about the time series of wind generation of the system.Each potential location candidate for the wind farm installation has its capacity factor (λ) and individual variations in its wind power generation series.erefore, the system wind power time series that are introduced in the Economic Dispatch Step cannot be precisely defined before completing the Optimization Step (Step 4, MINLP), where the individual wind power capacities are defined.However, it can initially be assumed that the system wind power series converges to a geographically smoothed or approximate time series, denoted by t APR h .e calculation of the value of t APR h is presented in detail in [4].t APR h satisfies the goal of wind energy penetration δ respecting hourly wind variations and, therefore, this is suitable for Step 2.

Economic Dispatch with Minimal Losses (Step 2).
e connection of a wind power to a transmission system not only affects the power flows due to the corresponding power injections, but also impacts on the entire operational dispatch of the system.us, the interdependence between load, existing and new, wind generation, and present conventional generation must be determined with the task of finding optimal locations for penetration of wind power.e increase in wind energy penetration changes the load demand requirements, which needs a dispatch response from conventional plants.
e geographically smoothed time series of wind power (Step 1) is subtracted from the system load, and the resulting net load time series serve as input data for the Economic Dispatch Step (Step 2).In this step, the economic dispatch for minimum losses corresponds to the lowest cost for conventional generation powers G i , considering the current wind energy penetration target.
A relevant technical issue of system operators regarding energy efficiency consists of trying to impact as little as possible the transmission losses due to generation redispatches, with no violation of lines limits, which justifies the loss minimization problem proposed in this paper.In this scope, the transmission losses are included in DC Economic Dispatch as shown by equations ( 1) and ( 2), according to methodology introduced by [11]: where P loss j is the loss of the branch j (MW), x j is the reactance of the branch j, g j is the conductance of the branch j, and θ j is the angular difference between buses of the branch j.
us, the inclusion of the loss of line j in each iteration t of the DCPF is done as in [11]: where θ (t−1) j indicates the value of θ j obtained in the previous iteration.is procedure is needed to avoid the inclusion of the nonlinearities related to losses in the DCPF.e proposal of the economic dispatch, instead of the unit commitment, given in [4], is justified by the focus of the present paper on the medium-or long-term planning, as well as by the adopted context of the Brazilian system, whose profile is mostly hydraulic.

Strategy for Formulation of the Transmission Line Constraints (
Step 3).After the DCPF, the obtained power injections P i in each bus i, except at the reference bus r, are used to calculate the power flows f j in each line j, as a linear combination of a set of linear coefficients α j,i ("power transfer distribution factors" or "distribution displacement generation factors") [12], equation ( 4).An important advantage of the DCPF is that linear constraints can be formulated to represent network power flow safety criteria in the optimization (MINLP) stage; this advantage stands out by reducing the complexity of the optimization problem.
e factor α j,i can be calculated according to the reactance parameters of the system [12], as shown in the following equation: where X � B ′ −1 (B ′ is the nodal admittance matrix, X ki is the element referring to the row k and column i of the matrix X, X mi represents the element in the position mi of the matrix X, and x j is the reactance of line j) [13].
With the information of the conventional generations coming from the economic dispatch of Step 2, the partial flows c j,h,s are calculated (for branch j, hour h, contingency scenario s) considering no injection of wind power.In this context, the inequality constraints of equation ( 6) are defined for every hour to ensure that the wind power penetration does not overwhelm any transmission line by violating its thermal capacities.Notice that these constraints are defined both in the direct and reverse flow directions.
where for every hour h, every branch j, every contingency scenario s, and considering k as the bus index of the wind farm in the network, L j (MW) is the thermal limit of line j,  6) can be manipulated algebraically and represented by (6a) and (6b).Equation (6c) ensures that the turbine capacity optimization variables are positive.
e computational requirements of the optimization algorithm are sensitive to both the number of variables and the number of constraints in the applied mathematical model [14].In this sense, the present paper proposes to include only the active constraints, i.e., the constraints for lines whose thermal limits are violated in the economic dispatch, to the MINLP of Step 5. is procedure seeks to obtain computational efficiency in the optimization process avoiding a large and unnecessary number of constraints.e strategy proposed in the present paper significantly reduces the percentage of the number of constraints that are active in the optimization problem, as will be seen in the case studies in items 3 and 4.

Mixed Integer Nonlinear Programming (MINLP) Formulation (Step 4).
In the proposed approach, the MINLP model is applied to determine the connection point for offshore wind farm that makes the wind power penetration level feasible with its highest utilization as possible.For this task, the wind power generation capacity C k for the current incremental wind energy penetration target δ is defined as the objective function (OBF), equation (7a), which is minimized.It can be highlighted that the minimization of C k implies in the choice of the best place for wind farm, i.e., with the best capacity factor λ, which means the best utilization of the wind resource.Considering the possibilities for connection of OWF, the proposed optimal connection model for maximum penetration of wind power capacity is formulated by the following equations: where n defines the total number of possible switches to connect the wind farm to the grid.e output of the MINLP model (7a)-(7g) is the optimal connection of OWF considering the maximum capacity of wind power that guarantees the feasibility of the wind power penetration target, i.e., so that the wind power goal is achieved and no line is overloaded for all contingency scenarios ("N − 1" criterion).
In theory, the feasibility of each firm wind penetration target δ is defined only by the constraints and should not be sensitive to the applied cost function.However, the choice of the cost function can affect the convergence of the general methodology and the allocation of wind generation capacity.
e energy contribution of a wind farm is defined by its capacity factor value. e optimal connection of OWF with maximum wind power capacity in the MINLP should fit the total wind energy penetration δ. is is ensured by including equation (7b) where each switching devices sw k and wind capacity C k are the optimization variables.
e equality constraints (7c) and (7d) are inserted into the mathematical model to ensure the optimal connection of one of the OWF connection options.Equation (7d) has two at a time combination of n switches that results in a set of n • (n − 1)/2 constraints, where a and b represent distinct OWF connections.

Update Wind Power Time Series (Step 5).
e formulation of the economic dispatch problem, in Step 2, uses as input a time series of geographically smoothed wind power.However, this may not necessarily correspond precisely to the time series of wind power resulting from the MINLP.us, to preserve the energy balance of the power system and correct dispatch of the conventional units at each operating time, Steps 2 to 5 are performed until the input series of Step 2 is equivalent to the output of Step 5.

Evaluation of Transmission Lines Limits (Step 6).
e purpose of this step is to determine whether or not the thermal limits of the network lines are violated after the solution found in the optimization process.If any line has been violated, Steps 2 to 6 are reiterated aiming at adding to the optimization model constraints for such line, to ensure that the solution found is feasible at the current penetration target.If no line is overloaded at all periods and contingency scenarios, the process moves to Step 7.

Increment the Wind Energy Penetration Target (Step 7).
After the wind power series have been updated and the power flows have met the thermal lines' limits, the wind energy penetration δ is increased in small discrete steps Δδ � 1%.us, as in Figure 1, Steps 1 to 7 of the methodology are reiterated until a feasible solution is obtained for the optimum connection of OWF with maximum firm wind power capacity.In this point, the algorithm ends and the optimal results are from the last feasible connection obtained.

Tutorial Example
4.1.General System Information.Figure 3 shows the unifilar diagram of the tutorial system used to illustrate the application of the proposed approach for maximizing the wind power generation.
is system is inspired by the 11-Bus described in [13].To consider wind generation, buses "12" and "13" are inserted into the system, where the reactances and resistances of the lines connecting these buses to the original system are equal to 0.01667 p.u. and 0.00100 p.u., respectively.e values of the active power demanded PL as well as the behavior of the wind power time series t k h , normalized in 1 MW, for a period of three hours under study are shown in Figure 4. e conventional generation capacity of the generators connected in buses "1," "2," "3," and "4" are 2100, 2100, 2154, and 2100 MW, respectively.e peak load demanded by the system is 3007.4MW. e average load over the period under analysis is 2734 MW. e thermal limits of the lines are shown in Table 1.
e analysis of optimal connection of OWF with maximum wind energy penetration is performed considering two options to connect the wind farm: at bus "12" or bus "13." e capacity factor λ k of the offshore wind farm is 28.37%.

Optimal Connection of Offshore Wind Farm with Maximum Wind Generation Capacity.
e MINLP model applied to determine the connection point among buses "12" and "13" is written for this tutorial system in (8a)-(8d).e inequality constraints follow the formulation presented in (7e)-(7g): As described in Section 1, it is initially considered the system operating with a low wind power penetration target, equal to 1%.In this case, the active power flows in the lines are within the respective thermal limits with no wind energy penetration.
en, this penetration is increased in small steps of 1% until the target of 11% makes the system operation not feasible.us, penetration increments of 0.1% are performed from the 10% level to find the maximum penetration value of wind power that respects the thermal line limits.Table 2 shows the values of the allocated wind power capacity C k for each penetration level and each candidate connection point.
As it can be observed, the OWF connected to bus "12" is the optimal option to obtain the maximum penetration (δ max � 10.3%) and maximum wind generation capacity (C 12 max � 992.67 MW), with the switches being arranged as follows: sw 12 � 1 and sw 13 � 0.
e infeasible level identified by the proposed optimization model for this case is given by δ � 10.4%.At this point, the power flow in Lines 7-12 violates the corresponding thermal limit (1000 MW) in all contingency scenarios, making it unfeasible for an additional injection of wind energy into the system.Any increase in penetration is also not supported by bus "13." e described results for this tutorial case are the same found by using the linear programming optimization presented in [4], considering only bus "12" as option for the OWF, which can validate the proposed approach in the present paper.

MINLP Inequality Constraints Used in
Step 5.For each penetration target, the insertion of line constraints, equations (7e) and (7f ), is considered to ensure that their thermal limits are not violated.In this sense, for the base case (without contingency), there are "14" limit constraints for the forward direction, equation (7e), added by "14" constraints for the inverse direction, equation (7f ), by hour of the period under analysis.It should also be taken into account that a number of "12" "N − 1" contingencies are analyzed (contingencies for all lines except lines 7-12 and 9-13 because they would generate conflicting constraints for the optimization process).us, for the period of three hours considered, there is a total number of constraints by each penetration level given in the following equations: where N constrains is the total number of constraints, N lines base is the number of transmission lines existing in the base case, N lines "N−1" is the number of transmission lines in a "N − 1" contingency scenario, and N contingencies is the number of contingencies under analysis.However, with the proposed strategy of selecting only the constraints related to lines that have their thermal limits violated, only "52" inequality constraints are added in this case throughout the optimization process, in which only the constraints concerning the lines 7-12 and 9-13 are violated, in the first hour, in both directions of the power flow.en, obtaining a percentage reduction of 95.67% thus improves the computational efficiency in the process.

Results and Discussion
is section presents the results obtained by applying the proposed methodology for optimal connection of OWF with maximum penetration of wind energy in well-known systems in the literature: IEEE 14-Bus [15], IEEE RTS-79 (a 24-Bus system), [16] and a Southern Brazilian system [17,18].
e IEEE 14-Bus system has originally 14 Buses and 20 lines, the total load demand of 259 MW and the maximum generation capacity as 2021.6 MW. e IEEE RTS-79 consists of a 24 buses system and 38 candidate lines for expansion, widely used in studies of transmission systems expansion, in which each branch can receive a maximum of three reinforcements; the total load demand is 2850 MW, and the maximum generation capacity is 3405 MW. e Southern Brazilian system has a basic network topology with  Journal of Electrical and Computer Engineering 46 buses and 79 branch candidates for expansion, where each branch can receive a maximum of three reinforcements; the total demand is 6880 MW, and the maximum generation capacity is 10545 MW.
For all the systems, it was considered that the behavior of the wind power series (1 MW) for the offshore wind generator in the 3-hour study period is similar to that shown in Figure 4; and that the load behavior followed the guideline of light, medium, and heavy load scenarios in the three hours under analysis (it can be extended to a larger situation of hours, without impairing the efficiency and robustness of the proposed model).that, the capacity factor (λ k ) of the offshore wind farm is 28.37%.
It is important to note that the proposed approach and analyses focus on the optimal connection of OWF, among some candidate buses geographically close to each other within a given region.e regions and respective buses are chosen at random and assuming that such regions are potential for wind energy, a strategy that is similarly used in [19], to obtain maximum penetration of wind power on a network whose structure is already consolidated.e network reinforcements required to achieve a higher penetration rate of wind power is not addressed in the present paper.
e proposed methodology presents robustness, both for smaller systems, e.g., IEEE 14-Bus and IEEE RTS-79, and for large systems, e.g., the Southern Brazilian system that is a more complex network.
e algorithm used to implement the methodology was developed in MATLAB [20], version R2016a, and utilized the FMINCON solver to solve the problems of Nonlinear Programming, using the "Sequential Quadratic Programming (SQP)" simulation model in the optimization algorithm configuration.e simulations were performed on a computer with an Intel Core i7-4500U, 2.40 GHz processor.

Case 1: System IEEE 14-Bus.
e first test system used to evaluate the methodology proposed in this work is a modified version of the IEEE 14-Bus system, similar to [3], presented in Figure 5.In this case, the candidate connection points for an OWF are the original buses "12," "'13," and "14," which are considered in a region with high wind potential, and the respective terminal buses of the OWF are buses "15," "'16," and "17." e "N − 1" contingency criterion is adopted.
Initially, it is considered the system operating with a low target of penetration of wind energy, δ � 1%.
e active power flows in the lines are within the thermal limits with no wind energy penetration.e values of the generation capacity of the conventional generators and the thermal limits of the lines are presented in Tables 3 and 4, respectively.
e wind energy penetration is increased in steps of 1% until a target whose the system operation is unfeasible, given by 22% in this case.us, penetration increments of 0.1% are performed from the target of 21% to find the maximum feasible penetration value.Table 5 shows the allocated wind generation capacities C k for each penetration level and candidate connection point.As it can be observed, the OWF connected to bus "16" is the best option, since this leads to the maximum penetration (δ max � 21.3%) and maximum wind generation capacity (C 16 max � 194.47 MW).For this optimal solution, the switching devices are sw 15 � 0, sw 16 � 1, and sw 17 � 0.
e maximum penetration level obtained in this case is given by δ � 21.3%.Beyond this point, the power flow in Lines 13-14 violates its limit (90 MW), in the first hour, for the contingency of Lines 9-14.Any increase in penetration is also not supported by buses "15" and "17." By adding all inequality constraints in Step 5 of the proposed optimization algorithm for both power flow directions, the number of hours in analysis, the number of lines (N lines base � 23 and N lines "N−1" � 22) and the number of contingencies (N contingencies � 20), the resulting number of inequality constraints is given by On the other hand, by using the proposed constraint selection strategy for transmission lines whose thermal limits are violated, the number of constraints is significantly reduced to 44, obtaining a percentage reduction of 98.42%, thus demonstrating the computational efficiency improvement through the proposed methodology.In order to validate the solution found about the choice of the optimal connection point for the OWF, the solutions for each of the possible connection paths were realized individually, wherein the maximum penetration targets found are equal to δ 15 max 17.0%, δ 16 max 21.3%, and δ 17 max 11.3%, associated with the OWF connected through switches sw 15 , sw 16 , and sw 17 , respectively, validating the optimal solution found through the proposed methodology.
e active power ows in the transmission lines are within the thermal limits with no wind energy penetration.
e values of resistances and reactances of the circuits of this system, as well as the generating capacity of the conventional generators and the lines active power ow, can be found in [16].For Lines 10-25, 6-26, and 8-27, the values of resistance and reactance are equal to r j 0.0023 p.u. and x j 0.0839 p.u., respectively, and the transmission limit is equal to 400 MW.
In this case, the wind power penetration is not feasible for a target of 4%.us, penetration increments of 0.1% are performed from the 3% to nd the maximum wind power penetration that meets the line limits.Table 6 shows the values of the wind capacity allocated C k for each penetration level and candidate connection point.
As it can be observed, the bus "25" is the best option with the maximum penetration (δ max 3.9%) and maximum wind generation capacity (C 25 max 391.81 MW).
e corresponding switches decision is given by sw 25 1, sw 26 0, and sw 27 0. Beyond the wind penetration limit (δ ≥ 4.0%), the power ow in Lines 10-25 overcome its limit (400 MW) in all "N − 1" contingencies, thus making additional wind  Journal of Electrical and Computer Engineering power unfeasible.Any increase in penetration is also not supported by buses "26" and "27."For this system, the number of transmission lines are given by N lines base � 41 and N lines "N−1" � 40, whereas the number of contingencies is N contingencies � 38, resulting in a number of inequality constraints given by equation (11).On the other hand, the proposed strategy that consists of considering only the active constraints leads to the significantly reduced number of 1530, obtaining a percentage reduction of 83.66%, as expected for the computational efficiency.
In order to validate the solution found about the choice of the optimal connection point for the OWF, the solutions for each of the possible connection paths sw 25 , sw 26 , and sw 27 were realized individually, wherein the maximum penetration targets found are equal to δ 25 max � 3.9%, δ 26 max � 2.9%, and δ 27 max � 2.3%, respectively.

Case 3: Southern Brazilian System.
e third test system is the modified version of the Southern Brazilian System [17,18], presented in Figure 7. is system is widely used for transmission expansion planning studies, so the thermal limits of the lines are very close to the values of the active power flows found in the base case.In this sense, it is not possible to penetrate wind energy in the base topology without violating the limits of some transmission lines, even at shallow levels.Due to such characteristic, it is considered, in this case, that the system is expanded by adding to it some candidate transmission lines, as shown in Figure 7 (dashed lines).It is important to note that this measure does not interfere in the efficiency of the methodology proposed in this paper since the proposal is to evaluate the maximum penetration of wind energy into a system whose structure is previously defined.e values of resistances and reactances of the circuits, as well as the generating capacity of the conventional generators and line limits can be found in [17,18].For Lines 44-47, 43-48, and 34-49, the resistance and reactance values are equal to r j � 0.01250 p.u. and x j � 0.08205 p.u., respectively, whereas their limits are equal to 600 MW.Buses "47," "48," and "49" are the candidate points to connect an OWF and the "N − 1" contingency criterion is also considered.For this practical system, an increase of 3% in wind energy penetration is already infeasible.us, the penetration increments of 0.1% are performed from the level 2% to find the maximum feasible penetration target.Table 7 shows the values of the wind capacity allocated C k for each penetration level and candidate connection point.
As it can be observed, the offshore wind park connected to the bus 48 was chosen as the optimal connection site to obtain maximum penetration (δ max � 2.4%) and maximum wind generation capacity (C 48 max � 581.94 MW), with the switches being arranged in the optimum configuration: sw 47 � 0, sw 48 � 1, and sw 49 � 0. Any increase in penetration is also not supported by buses "47" and "49." Above the maximum feasible penetration (δ ≥ 2.5%), the power flow in Lines 43-48 overcomes the respective limit (600 MW) for all "N − 1" contingency cases, thus making an additional injection of wind power into the system unfeasible.
In Case 3, there is N lines base 123, N lines "N−1" 122, N contingencies 120, and when the strategy for selecting only the active line constraints is not applied: On the other hand, when only the active lines are selected as the proposed approach, the resulted number of constraints is reduced to 3232 for increasing the computational e ciency.
e solutions found for each OWF connection possibility, when they are realized individually through the sw 47 , sw 48 , and sw 49 switches, are δ 47 max 2.1%, δ 48 max 2.4%, and δ 49 max 1.5%, respectively.With these results, the solution found through the proposed approach is validated.

Runtimes.
e typical runtimes for the evaluated test systems by using the proposed methodology are presented in Table 8.

Conclusions
In this paper was proposed a novel approach for determining the optimal connection of an o shore wind farm to the  onshore network, among some regionally close candidate points in a given region, aiming at obtaining the maximum wind power penetration, observing the network transmission limits and the "N − 1" contingency criterion.In addition, the formulation allows representing the stochastic nature of winds and the di erent behavior of the load, by using hourly time series.Network security criteria were analyzed for each hour of the time series, thus generating a large dimensionality of power ow constraints.For handling this feature, an efcient strategy is proposed to select only active line limit

Figure 2 :
Figure 2: Flow chart of the proposed algorithm.

Figure 5 :
Figure 5: Uni lar diagram of the modi ed IEEE 14-Bus system.

Figure 6 :
Figure 6: Uni lar diagram of the modi ed IEEE RTS-79 system.

Figure 7 :
Figure 7: Uni lar diagram of the modi ed Southern Brazilian system.

4
(1)rnal of Electrical and Computer Engineeringt k h (MW) is the wind power time series (normalized in 1 MW) in hour h, c j,h,s (MW) is the predefined partial flow, sw k represents the switch that connects(1)or not (0) the wind farm to the grid, and C k (MW) is the optimization variable of wind capacity.e inequalities on both sides of equation (

Table 2 :
Allocated capacity (MW) for each penetration level and candidate connection, tutorial case.

Table 3 :
Capacity of the conventional generators, case 1.

Table 5 :
Allocated capacity (MW) for each penetration level and candidate connection, case 1.

Table 6 :
Allocated capacity (MW) for each penetration and candidate connection, case 2.

Table 7 :
Allocated capacity (MW) for each penetration level and candidate connection, case 3.