Data of the 28bus grid.
Abstract
This chapter presents two probabilistic planning tools developed for the longterm analysis of distribution networks. The first one focuses on the lowvoltage (LV) level and the second one addresses the issues occurring in the mediumvoltage (MV) grid. Both tools use Monte Carlo algorithms in order to simulate the distribution network, taking into account the stochastic nature of the loading parameters at its nodes. Section 1 introduces the probabilistic framework that focuses on the analysis of LV feeders with distributed photovoltaic (PV) generation using quarterhourly smart metering data of load and generation at each node of a feeder. This probabilistic framework is evaluated by simulating a real LV feeder in Belgium considering its actual loading parameters and components. In order to demonstrate the interest of the presented framework for the distribution system operators (DSOs), the same feeder is then simulated considering future scenarios of higher PV integration as well as the application of mitigation solutions (reactive power control, P/V droop control thanks to a local management of PV inverters, etc.) to actual LV network operational issues arising from the integration of distributed PV generation. Section 2 introduces the second planning tool designed to help the DSO, making the best investment for alleviating the MVnetwork stressed conditions. Practically, this tool aims at finding the optimal positioning and sizing of the devices designed to improve the operation of the distribution grid. Then a centralized control of these facilities is implemented in order to assess the effectiveness of the proposed approach. The simulation is carried out under various load and generation profiles, while the evaluation criteria of the methodology are the probabilities of voltage violation, the presence of congestions and the total line losses.
Keywords
 Design of experiments
 Dispersed generation
 Distribution systems
 Monte Carlo simulation
 Smart metering
1. Introduction
Following the Kyoto Protocol and in the context of a liberalized energy market, decentralized generation, often based on renewable energy sources, is emerging in medium voltage (MV) and lowvoltage (LV) distribution grids. Those networks are becoming more and more active with power flows and voltage profiles influenced by both generation and consumption. In the future, given the ‘20–20–20’ objectives of the European Union and the more ambitious objectives in Wallonia (Belgian region), the expected penetration of decentralized generation in distribution grids could lead to a critical behaviour of the system that has not been initially designed and sized to face power injections coming from dispersed units. This major paradigm shift in distribution grids is even strengthened by the numerous financial incentives granted to the renewable generation for ensuring their profitability. Indeed, this financial support has progressively driven down the electricity prices on power markets to the point of making the conventional generation less profitable, which already led to the closing or the mothballing of more than 50 GW of gas plants in Europe during the last few years [1].
Consequently, in the MV grids, the main problems due to a massive integration of distributed generation (mainly wind parks) are due to local congestions of power lines and/or to overvoltages in the neighbourhood of the dispersed generation units. Concerning the LV grids, the major problem that can be met by the distribution system operators (DSOs) comes from overvoltage in the neighbourhood of photovoltaic (PV) installations. Such a situation could happen during periods of low consumption and high generation of those decentralized units, especially on long power lines. Furthermore, reverse power flows going up the line towards the MV/LV transformer can then degrade the network stability.
The DSO is responsible for the security of its system and power quality. It is therefore important to develop planning tools that ensure secure and reliable operating conditions within the whole grid. Moreover, it is also essential for DSOs to avoid early damages of the network equipment as well as to optimally manage the investment decisions and the maintenance plans. However, as the intermittent generation is often difficult to predict due to its stochastic behaviour, the planning studies become more complex to carry out. Indeed, it is now indispensable to use probabilistic approaches to model the different network components, since the traditional deterministic worstcase scenarios yield much more restrictive conclusions that can lead to overestimated reinforcements [11].
In that way, two probabilistic tools are described in this work. Both are based on Monte Carlo (MC) simulation whose general principle is to generate a large number of system states in order to provide consumption and generation patterns that are representative of the actual behaviour of all customers. However, it is important to mention that longterm analyses are focused on modelling realistic scenarios of the system evolution rather than making an accurate prediction. It is indeed utopian to expect a precise hourly model of the wind farm generation for a whole year. Practically, two different types of Monte Carlo simulation are implemented. The first one is a nonsequential simulation that generates a set of system states independent of each other. This method is combined with the design of experiments (DOE) methodology in order to proceed to an optimal positioning as well as to a presizing of voltage compensation devices (capacitor banks, storage units, static VAR compensators (SVCs), etc.) at the MV level. The purpose is to help the DSO in optimally managing its systems with adequate investments. The first tool requires aggregated load profiles at the MV/LV interface. Those statistical profiles are typically obtained, thanks to a second tool developed for the safe integration of PV units in LV grids. The latter is implemented in a pseudochronological Monte Carlo environment and is based on quarterhourly energy flows recorded by smart meters (SMs). Practically, thanks to statistical profiles established by the use of these measurements, the tool is able to compute different parameters throughout a specified period of time such as the voltage profile along the studied feeder or the distribution of the power profile at the MV/LV substation. In this chapter, in order to demonstrate its interest for DSOs and to introduce the potential of some technical solutions (reactive power control, P/V droop control thanks to a local management of PV inverters, etc.) to alleviate the influence of PV generation in LV grids, this tool will be applied to a practical case study.
This chapter is thus divided into two main parts, each one focusing on one of both the abovementioned probabilistic tools. Section 1 concerns the description of the planning at the LV level. More particularly, the extraction of the energy flow data coming from SM installed in the LVBelgian distribution grid as well as their use in the Monte Carlo procedure is presented. Then, the potential of several power or voltagecontrol strategies is evaluated for an existing LV area subjected to abnormal operating conditions. Section 2 relates to the planning tool studying investment plans that may have to be taken in critical areas of the MV grid.
2. Probabilistic analysis of lowvoltage networks
The technoeconomic analysis of LV feeders is typically governed by a set of basic rules that define voltage, voltage unbalance and current limits. According to these rules, which are most often precisely determined by regional or national standards, voltage along the feeder should fall within a defined band during a minimum percentage of time. As far as the voltage unbalance is concerned, an upper limit should also be respected during a minimum percentage of time. The limits that concern the value of the current do not often apply in LV feeders apart from urban networks with high load density.
When adding distributed photovoltaic generation to an LV feeder, voltage profile is the parameter that is mostly affected. Indeed, during periods of high PV injection and low consumption, reverse power flows towards the head of the feeder can occur. These reverse power flows result in a voltage rise at the end of the feeder. According to the European EN 50160 standard [2] and its national implementations (i.e. the voltage cannot vary more than 10% around its nominal value), if the upper limit of the root mean square (rms) voltage is exceeded at a certain PV node, the unit must be temporarily cut off [2], [3]. This causes a loss of generated PV power, which means a loss of income for the PV owner. Moreover, these frequent cutoffs, combined with the randomly changing loading profile along LV feeders, affect the distributed power quality and accelerate the degradation of network components. The operational cost of the network, currently assigned to the DSO, is therefore increased due to the accelerated ageing of its components.
The development of technical solutions to such problems has, from the DSO point of view, a threefold objective, namely the security level of the network, the increase of PVhosting capacity of the feeder and, naturally, the investment cost. For a given PVhosting capacity, the optimal design of a technical solution could be roughly determined by a graph similar to Figure 1. On the other hand, for a selected security level, the optimal hosting capacity in a feeder could be based on a graph similar to Figure 2.
For designing adequate technical solutions to current problems while increasing the PVhosting capacity, the DSO needs to locate the critical points and quantify the voltage magnitude and unbalance violation risk in every LV feeder. Currently, the DSO evaluates the LV network operation using deterministic analysis tools: the energy exchanges at the nodes of a network with distributed PV generation are determined based on worstcase scenario. For example, as far as the evaluation of overvoltage risk is concerned, the DSO makes its calculations considering the scenario of the highest expected PV injection and the lowest consumption load along each specific feeder. Although worstcase approach is safe for the protection (100% security level) of the electrical system, it focuses on rare extreme operation states of the network instead of the most possible ones. As a result, the worstcase approach often leads to oversized and costly technical solutions that do not consider the time fluctuation of PV injection and the randomness of consumption loads. Taking this time dependence into account will allow a refined optimized design of mitigation techniques in terms of efficiency and investment cost, closer to the optimal level that is roughly represented in Figure 1. For the above reasons, a doubt concerning the effectiveness of the deterministic worstcase models has recently arisen in literature [4], [5].
In order to take into account the time variation of the PV injection and the random behaviour of the loads and PV installations along an LV feeder, probabilistic analyses were introduced in late research works [6, 9]. These methodologies are based on analytical or numerical approaches, such as the Monte Carlo simulation, aiming at simulating the uncertainty of PV injection and residential loads. Nevertheless, most of these studies propose advanced probabilistic models based on solar irradiation data for the PV panels and meteorological data for the consumption loads. Such data are rarely directly accessible to the DSO and most of them do not consider the efficiency of the PV cells. Currently, DSOs in several countries are installing smart metering devices that record voltage, current as well as active and reactive powers with a high time resolution (10 or 15 min), even for domestic customers. According to European Photovoltaic Industry Association (EPIA), large amounts of such data will be available at the horizon of 2020 [10]. The processing of these data is, nevertheless, still a challenging task that will have to be addressed.
2.1. Statistical modelling of the energy flows at LV nodes
A procedure for building statistical profiles of energy flows at LV nodes and of the voltage at the MV/LV substations is presented in references [11, 13] and will be summarized here. The data are recorded by SM devices that have been installed by the DSO in the Belgian grid but the procedure can be extended to any LV grid with such metering devices. The SM configuration of every PV customer is presented in Figure 3.
The SM1 records the total quarterhourly (15 min) injected energy towards the network (
For the purpose of evaluating the overvoltage risk along a feeder, the SM1 data, namely the energy flows at the coupling point of the PV customer to the network, are most preferably used. These 15min values are therefore used for the creation of typical day profiles (TDPs) for every individual quarter of an hour in a day (96 quarters of an hour). Created independently for each node (customer), the TDPs reflect the variation of energy flow, separately for injection and consumption, for every 15min interval in a day, during the studied period. Similar TDPs are also created to reflect the variation of the voltage at the head of the feeder, based as well on 15min SM voltage data recorded at the LV side of the MV/LV transformer. It is worth noting that such profiles can be very easily created for various averaging time steps of SM devices. In most cases, SM devices in LV networks in Western Europe record energy flows on a 15 or 10min time resolution.
The TDPs are then statistically transformed into 96 cumulative distribution functions (CDFs), that is, one for each quarter of an hour of the typical day, and they will be used in the framework of an MC probabilistic load flow.
2.2. Probabilistic load flow using a Monte Carlo method
The Monte Carlo simulation is a mathematical technique that allows accounting for the random behaviour of some variables within the studied system. In our case, the stochastic variables are the consumption and the generation of each customer of the LV grid during the day. The principle relies on the repeated random sampling of these variables in order to accurately simulate the system and to obtain numerical outcomes (i.e. probability of overvoltage at each node of the LV grid).
The abovementioned CDFs are sampled under an MC framework in order to generate randomly 15min network states. A balanced powerflow algorithm then computes the voltage profile along the feeder considering the sampled energy flow values at each node and the sampled voltage at the head of the feeder. The procedure is repeated numerous times in order to test a large amount of possible combinations. For obtaining a good convergence threshold (<0.1%) on the results, it is shown that 10,000 simulations is a good tradeoff to keep acceptable computation times [11].
However, for performing the powerflow computation, the total 15min energyflow values must be transformed into instantaneous power values at each node. The injected or consumed power is rarely stable during the 15min interval, which means that the total recorded energy is very often injected or consumed in a very short (<<15 min) period of time. Indeed, it has been shown in reference [12], thanks to advanced grid analysers installed by the DSO at critical network points that the power flow at the coupling point of a grid customer is highly fluctuating within the 15min intervals. Therefore, the peak values of injected or consumed power might be much higher than the average ones that consider a stable power flow.
In the proposed probabilistic framework, for the sake of security, these instantaneous power values are not considered equal to the average values, which would be computed by dividing the total injected or consumed energy by the respective period of time. As explained in reference [12], TDPs are also created for the time repartition of the energy flow within the 15min intervals based on 3min recordings of advanced grid analysers at limited critical points of the network. These TDPs are also integrated in the MC algorithm, once they have been transformed into CDFs, in order to provide the sampling frame for the time repartition factor
The computed peak power values are therefore the ones considered in the powerflow algorithm. The rms voltage in this way is calculated node by node for each 15min network state, which is therefore characterized by the following values:
These 15min values are compared to the voltage band that is imposed by national, regional or local values. In the case of Europe, such limits are imposed by the EN 50160. In this way, the overvoltage probabilities can be determined, node by node, for the studied period.
2.3. Illustration on an LVBelgian network
The LV feeder, which is illustrated in Figure 4, makes part of the LVBelgian distribution network. This feeder (technically described in reference [12]) is simulated with the presented probabilistic framework for every month of a typical year, considering the 15min SM data recorded over a period of two years. In total, 16 clients are connected to the feeder, among which three are equipped with a PV installation (5 kVA at node 4, 2.65 kVA at node 12 and 5 kVA at node 13). In order to clearly demonstrate the influence of the PV generation on the voltage variation, six more PV clients have been considered at nodes 5, 6, 8, 10, 11 and 14 for this simulation. For these additional units, the same PV injection SM data as for the PV unit at node 13 were considered.
Time variation for a given MC simulation of the voltage profile at node 14 is shown in Figure 5. Calculated with a deterministic analysis, the voltage would have been a stable worstcase value that would characterize the whole period of a day. However, based on the current analysis and the use of SM data, it is clearly shown that the voltage rises only during certain hours of the day, usually between 11:00 A.M. and 16:00 P.M.
Thanks to the probabilistic simulation, the overvoltage risk can be defined for each node during the studied period. This risk is defined as the probability of exceeding the limit (>1.10
A significant voltage drop takes place between 18:00 and 21:00 hrs, due to peak electricity demand, which increases undervoltage risk (<0.90
In Figure 6, it can be seen that nodes adjacent to PV units are also affected by rapid voltage variations and oscillations (caused by the on–off control of the units), which tend to accelerate the degradation of many network components. On the other hand, it has been shown in reference [14] that the number of line congestions and the transformer load are significantly impacted with the emergence of PV generation and electric vehicle (EV)charging devices. In this way, Figure 8 illustrates the decrease of line congestion when PV generation is included into the grid. Indeed, the power produced by these units is consumed locally by the neighbouring loads, which reduces the line power flows. Conversely, the introduction of EVcharging devices has an adverse impact on the grid. Indeed, these are powerful appliances (typically 230 V and 16 A) that significantly increase the current in different lines. Figure 9 shows the statistical distribution of power flows at the MV/LV substation. One can see that the power at the MV/LV interface decreases when generation is installed into the grid and substantially rises with the integration of EVcharging appliances. Therefore, three different scenarios are investigated in this study. Scenario 1 considers the feeder represented in Figure 4 but with no PV generation. Scenario 2 considers the same feeder with PV units at nodes 4, 5, 6, 13, 14 (rated power 5 kVA) and at node 12 (rated power 2.65 kVA), whereas scenario 3 considers also the EVs connected to the same households than those with PV panels. It should be noted that the EVs are simulated based on typical EVcharging profiles and on the statistical behaviours of current EV owners.
The previous results highlight the necessity to adopt strategies aiming at limiting these important voltage and current variations. Otherwise, actions can be too restrictive for potential PV hosting capacity of the network or much less optimized in terms of efficiency and cost as roughly represented in Figures 1 and 2. Moreover, the probabilistic framework can also simulate future scenarios of PV integration increase, or technoeconomically evaluate multiple solutions to network operational problems. It can therefore provide answers for an efficient and secure operation of LV networks, taking into account the optimal level that is represented in Figures 1 and 2.
2.4. Evaluation of the strategies for mitigating the overvoltage risk
Currently, certain DSOs evaluate the potential benefit of loadshifting strategies on mitigating overvoltage risk and increasing the PV energy capture. Incentives that promote selfconsumption of the generated PV power by its producer should be designed based on the time variation of the network states. In line with this assumption, the selfconsumption rate of 20 residential PV customers has been studied by using the 15min recordings of the SM2 devices. The selfconsumption ratio of each PV customer can be calculated for each quarter of an hour in a day by using both the SM1 and SM2 data and applying the following relation:
The average selfconsumption ratio of three residential PV customers is shown in Figure 10 for each quarter of an hour, based on SM1 and SM2 data for the month of July. According to this diagram, the selfconsumption ratio decreases during hours of peak PV injection, differently for each one of the clients. One could therefore assume a quite high potential for strategies that could shift consumption loads towards high PV injection hours. Such strategies could achieve an important mitigation of the overvoltage risk that usually increases during these hours as shown in Figure 5. The design of incentives could be therefore based on a thorough study of the selfconsumption behaviour of customers connected to similar LV feeders, by using SM recordings and suitable clustering methodologies.
Apart from load management strategies, certain power or voltage control strategies are currently discussed by the DSOs for facilitating the integration of renewables in the LV network. Among them, reactive power control is a strategy that has been recently implemented in the LV networks of some European countries. In these cases, the distributed PV units have to be equipped with Q/V droop functions in order to provide voltage support as the conventional large power plants. The reactive power of the PV inverter is then adjusted, based on the local voltage, according to a
In this study, a simpler application of this control scheme, that does not consider the actual local voltage, has been simulated in order to evaluate its benefit on the mitigation of the overvoltage risk. The purpose of this example is to illustrate the potential of this probabilistic tool in order to constitute a benchmark for evaluating and comparing voltage control solutions. In the previously presented results, the reactive power generated by PV inverters has been set equal to zero in order to minimize power losses. Turytsin et al. demonstrate that the optimal adjustment of PV inverter reactive power could reduce voltage variations. The maximum available reactive power at each PV node
where
In European countries where such a control scheme has been implemented in the LV network, the value of the power factor (
Presently, the PV power injection towards the network is solely subject to the maximum power point (MPP) of the PV units. In case of voltage violation, the PV inverters passively undergo the on–off control that is required by the European or national standards, without considering the current network state. In order to address this operational issue, T. L. Vandoorn and colleagues developed a fastacting primary control scheme based on voltage droops [22]. The droops are applied by P/V controllers requiring neither interunit communication nor voltage tracking for synchronization. In this study, P/V droop control is evaluated with the proposed probabilistic framework in the same LV feeder in order to analyse its benefit on the voltage level, the captured renewable energy and the on–off oscillations of the PV inverters. The evaluation considers the time variation of nodal injection and consumption and is therefore based on the analysis of multiple network states and not only the worstcase ones.
The P/V control scheme modifies the injected active power of the PV unit in function of the local voltage in order to prevent and eliminate overvoltage situations. In this way, the total cutoffs of PV units, applied in the conventional on–off control, and their subsequent effects (voltage and current transients as well as significant PV power loss) are avoided. The P/V droop controller adapts
For the purpose of evaluating the benefit of the P/V droop control, all PV inverters were considered equipped with a P/V controller. The application of this control in the studied LV feeder, considering SM data for the month of July, demonstrated its efficiency for eliminating the overvoltage risk along the feeder [23]. In Figure 13, the voltage profile at node 14 for the droop control and for the on–off control scenario is indicatively presented for a period of two typical days of July. Generally, the voltage profile of the P/V control scenario coincides with one of the on–off scenarios apart from the three 15min intervals of the first day and two 15min intervals of the second day during which the P/V control was enabled. Indeed, during these periods, the rms voltage at node 14 exceeded the reference value of the droop control
Obviously, this gradual smoothing of the voltage profile induces the curtailment of an amount of active PV power. As demonstrated in reference [23], the instantaneous curtailed PV energy can often be much more direct and higher for the on–off control. This argument can foster the acceptance of a certain control strategy, such as the P/V droop control, by the concerned PV customers. The longterm impact of such mitigation solutions on the income of the PV producer and on the operation and maintenance cost of the network, which is currently assigned to the DSO, should be evaluated keeping in mind the principle of Figures 1 and 2.
The presented evaluation that takes into account, in a fast optimized manner, multiple network states and not only the worstcase ones can lead to technical solutions that should be tailored to the special needs and necessary security level of each individual network. Moreover, the analysis and design of the LV distribution network with the presented probabilistic framework contribute to a more rigorous analysis of the MV distribution network. Indeed, when studying the MV network, downstream LV feeders are reduced to P–Q buses defined by the power flow at their root nodes. Most often, the time variation of this power flow, due to PV injection and random loads, is not taken into account. Section 2 introduces a probabilistic methodology that aims at analysing the MV distribution network and uses aggregated load and generation profiles at the MV/LV interface obtained thanks to the presented probabilistic framework. The impact of time variation of load and generation profiles in the LV network in this way is considered in the analysis of the MV network as well.
3. Investment planning at the MV level
Originally, distribution networks were oversized in order to easily ensure the connection of new clients without requiring the upgrade or the replacement of the equipment. However, the constant increase of load and generation, although the network infrastructure has not increased accordingly, induces a global system operation closer to its limits. In particular, this may lead to excessive line power flows that can be destructive not only for the network equipment but also for the enduser material. Indeed, important current flows may infer unintended line congestions as well as an extensive increase of the transformer load that can potentially reduce its lifetime [24]. Moreover, as the voltage drop depends not only on the impedance of the line but also on the line power flows, the voltage profile can locally fluctuate outside its acceptable limits. Although the overvoltages are more dangerous for electrical devices, both low and highvoltage issues have to be taken into account in the regulation of distribution networks. However, this objective of ensuring a good voltage profile can be in conflict with the minimization of line losses that are very costly for the DSO, in the range of 50 euros/MWh, concerning the Belgian situation. Indeed, for a same demand of power, the reduction of line losses is obtained by a voltage rise. Consequently, the DSO will look for a tradeoff between the quality of the delivered power and the economic losses.
Currently, several strategies can be used to improve the realtime operation of the grid. The most popular method is the use of onload tap changer (OLTC) of transformer. The principle is to mechanically adjust the turn ratio of the transformer winding, usually with regard to a local measurement of the current. It should be noted that the dynamic of an OLTC is relatively fast with a total operation time between 3 and 10 s, depending on its design [25].
Then, the reactive power compensation devices are also very useful to get a better voltage profile in a power system with an X/R ratio higher than one. Devices such as the static VAR compensator or static synchronous compensator (STATCOM) have the interesting ability to work in both inductive and capacitive modes. Furthermore, the technology is well developed and financially affordable compared to other techniques [26].
In recent years, the possibility of investing in storage is gaining credibility [27]. Indeed, although the most mature technology (i.e. the pumped hydro) is not appropriate for the needs of distribution systems, other technologies are expected to become increasingly important in the near future. In this way, stationary batteries seem really appealing due to their good flexibility regarding the input and output power. Nonetheless, the cost is still very high and their success will largely depend on the cost reduction. Then, the compressedair energy storage (CAES) could be a solution in the coming years, but its integration is hampered by its low efficiency and limited operational flexibility. Finally, the hydrogen storage technology is showing promise but is not expected to be sufficiently developed to be competitive before 2020.
Then, the possibility to use the reactive power produced by distributed generation (DG) units is also an investigated alternative [28]. However, this method has a limited margin due to the current legislation restricting the reactive power flows at the connection node of each client (generators as well as industrial customers). Moreover, as there is currently no market for ancillary services at the distribution level, the customers have no incentive to help the DSO in its task of maintaining optimal operating conditions.
Finally, if the power flows far exceed the network capacity, the system reinforcement may constitute the most viable solution.
Note that the DG units connected to the distribution grid are automatically disconnected when the network frequency or the local voltage exceeds the limits defined by the legislation. Therefore, another major asset of a local control of the network operating conditions is the avoidance of such disconnections, which may have to be financially compensated, depending on the contract between the producer and the DSO.
In this study, a new procedure aiming at ensuring the respect of the operational constraints of the distribution grid at any time is implemented. This amounts to avoiding any congestion or voltage violation. The crucial role of maintaining the frequency balance is indeed assigned to the transmission system operator (TSO). The general structure of the proposed method is represented in Figure 14 and is divided into three parts. First, a twostep planning procedure is conducted for identifying the optimal strategies that the DSO has to implement in order to optimally improve the longterm operation of the MV network. As most of the Belgian HV/MV transformers are already equipped with OLTC devices, and given the current context not fostering the reactive control of DG units, the study is mainly focusing on the installation of VAR compensators and storage utilities. In this way, the first part of the planning process evaluates the best positioning of these network regulation devices, whereas the second part consists in sizing the devices by taking cost constraints into account. Obviously, the final choice of the positioning will depend on the geographic and administrative considerations, which will not be considered here, as they have no repercussions on the philosophy and the principle of the method. The second part consists in simulating the realtime centralized control of the system, which is performed by optimizing the command of the network regulation devices in function of the network state. This operation is thus carried out for a large number of simulated states. Finally, the last part is postprocessing whose purpose is to compare the probability of constraint violations as well as the total line losses before and after the installation of the network regulation devices.
The developed methodologies for the planning procedure and the realtime simulation of the network states are both developed in a nonsequential Monte Carlo procedure due to the lack of realtime measurements, especially at the MV/LV interface. This issue concerns less the industrial area, as quarterhourly energy flow recordings are mandatory for the Belgian consumers with an installed power higher than 100 kVA. Indeed, these clients are not only invoiced for their actual energy consumption but also for their peak power, regardless of its duration.
When working in a Monte Carlo environment, it is essential that the underlying simulations are fast and efficient in order to avoid an excessive simulation time. Therefore, because of its excellent tradeoff between simulation time and quality of the results, the experimental design method was privileged. Indeed, as explained hereafter, this method, which is also referred to as design of experiments, is perfectly suited for conducting the desired studies (i.e. planning procedure and simulation of the centralized control of network regulation devices). Then, as both proposed studies require performing loadflow studies in order to determine the voltage profile and the power losses within the system, a fast, simplified computation, based on the analytical approach presented in reference [29], was implemented. However, in its initial form, this algorithm was designed for a single radial feeder without ramifications; it had to be thus adapted in order to be used on the traditional treeshape structure of MV distribution networks.
Since the first issues originating from the progressive introduction of renewable energy generation over the past years, the improvement of MV network operation has become a real matter of concern. In this way, in recent studies [30], [31], new voltagecontrol methods combining the advantages of OLTC action and DSTATCOM response were proposed. However, in practice, this procedure can hardly be implemented as the transformers are not necessarily directly managed by the distribution system operator (i.e. most of the Belgian transformers, even those installed at the distribution level, are the property of the TSO).
This study is intended to be more complete and takes into account the reality of the economic and technical constraints by evaluating the impact of a local centralized control downstream each of the HV/MV transformers. The control is thus at the interface between a decentralized control performed at each connection node of the network regulation devices and a centralized control applied on a national scale.
Practically, this work is structured as follows. First, the principle and the basic concepts of the DOE method are introduced. Then, the proposed methodology is fully presented and applied to an existing MVBelgian network. Finally, the collected results are exposed and discussed.
3.1. General principle of the DOE method
Design of experiments constitutes a powerful tool to establish and study the effects of multiple inputs on a desired output. The method thus involves two types of variables, the response of interest
where
There are different experimental designs, which are commonly divided into two main categories, given the desired objective. The first one constitutes a screening study, which is used to determine the factors that really influence the outcome. This procedure allows therefore keeping only the most influential factors in the context of a selection process. The second one, which is called response surface methodology (RSM), is performed when the objective is to find the optimal settings of the factors in order to obtain the desired system optimization. Both methods can therefore be used in a complementary way.
Each design is characterized by its own mathematical approximation of the response
The parameters
where,
Finally, the fitted regression model is:
where
3.2. The screening of significant parameters
The screening study is used to quantify the significance of each factor on the response and is therefore very useful in the context of selection procedure or reduction of the problem size. Here, this method is applied in order to find the best area for installing the network regulation devices. Note that the goal is not to precisely determine an exact location, since the probability that this node can really be used in practice is highly insignificant. Furthermore, this approach allows to substantially decrease the number of factors and to generate important time savings. The underlying aim of the study is thus to avoid the upgrading of the feeders, as such a solution is not economically viable for DSOs, especially considering the current trend to install underground power cables. Practically, two different screening studies are carried out, one for each of the two envisaged control means (i.e. VAR compensators and storage devices). Recall that, due to the relatively low extra costs compared with normal transformers as well as their high efficiency, OLTC devices are already installed in HV/MV substations. Moreover, the control of the reactive power of DG units does not seem feasible at the present time. At the end of the screening study, we will therefore obtain the most favourable areas for installing VAR compensators and storage facilities.
The response
where
where
where
Beyond its simplicity of implementation, the DOE approach presents another highly interesting asset. The coefficients of the mathematical model have indeed a physical meaning, and their values can be interpreted as a quantitative reflection of the effect of the factors on the response. This interesting and useful property is illustrated in Figure 15 for
The number
First, the network topology is used to define all the possible positioning configurations of the network regulation devices. In this study, these are either reactive power compensators or storage devices. This leads to the definition of the experimental matrix
In order to take the nondeterministic effects of the load and the generation into account, the screening process is implemented in a nonchronological Monte Carlo environment. Such a study allows simulating a large number of operating states of the network. For each of them, the influence of the network regulation devices is evaluated. The process thus simulates the random behaviour of the different MV clients as well as the MV/LV substation by the means of a random sampling of the consumed and/or injected power at each bus. More details about the procedure are given in Section 3.5.1.
The importance of the location of each area for installing network regulation devices is then computed by averaging the effects provided by the screening procedure conducted for each of the
3.3. Sizing procedure of the selected devices
After the determination of the most favourable areas for respectively installing a reactive power compensator and a storage utility, the sizing of these devices is carried out. As one of the objectives of the method is to illustrate the assets of implementing a centralized control at the MV level, the sizing procedure is applied for more than one network regulation device. For simplicity as well as economic reasons, only two devices are considered in this work. Moreover, as the storage technology is still too expensive, the possibility to install two storage utilities is not considered. As a result, two different scenarios are investigated. The first one focuses on the combined installation of a VAR compensator and a storage device, whereas the second one evaluates the possibility to invest in two reactive power compensators.
The sizing of the storage utility concerns only its output power. Indeed, the energy sizing is less critical for the DSO than it would be for an operator eager to maximize its return on investment, since the DSO is not to take the market considerations into account. Its only focus is to use the storage device when it is necessary for alleviating stressed conditions within its grid and to take advantage of the noncritical situations to adjust the amount of stored energy. However, such a control strategy is outside the framework of this study.
As in the screening process, the sizing procedure makes use of Monte Carlo simulations for modelling the different network states. For each simulation, an optimal control of the selected network regulation devices is carried out. Thus, this requires the use of a fast optimization tool, which led to the choice of the response surface methodology. The response is the same as the one used in the screening process (i.e. a weighted sum of the constraints violations); and the response surface thus models the variation of this response with regard to the output power of the two selected regulation devices. The mathematical model used for the study is of the form:
where
The purpose of optimization is to find the configuration of the control devices that eliminates all constraint violations while minimizing their interaction (i.e. power exchanges) with the grid. Mathematically, it amounts to find the intersection between the response surface and the plane
The optimal values given at the end of each Monte Carlo simulation are logged. In this way, at the end of the simulations, the DSO disposes of useful information. Indeed, the output power of devices can be directly translated in installation costs and, given the different values of the regulation devices, the DSO is then able to adapt its investment decision as regards the technical and economic considerations.
3.4. Centralized realtime management of the distribution grid
The method proposes to use a centralized control of each of the network regulation devices of the considered distribution grid. Indeed, the current SVC/STATCOM devices compare the voltage at their connection node to a fixed reference voltage in order to bring back the voltage at the defined reference level. But sometimes, it could be interesting to change that reference depending on the location of the device in the grid. For example, to lower it when there is an increased power generation coming from DG units combined with a low demand. In this context, the output power of the storage devices would be adapted by a remote control of power electronics.
The strategy of the control depends directly on the presence of operational constraints violations. On the one hand, if there are such violations, the aim is to overcome the issue with an optimal use of the network regulation devices. However, with several devices, several configurations can potentially solve the problem. In such a case, the DSO is interested to opt for the configuration that minimizes the line losses. On the other hand, if there is no problem, different strategies can be envisaged. The most straightforward option is to do nothing, which will minimize the number of cycles of the devices and, as a result, increase their lifetime. For the sake of simplicity, this strategy was implemented in this study. The second one is to use the control devices for reducing the line losses. Such an action has to be carried out only if it generates worthwhile financial profit. Finally, if a storage facility is available, its energy level can be adapted in anticipation of the future needs.
The optimization is carried out by using the same response surface methodology than the one used in the sizing procedure (i.e. same response and same variables).
3.5. Application of the DOE method to an existing radial MV grid
3.5.1. Modelling of the MV network
The simulation of the statistical behaviour of the different industrial loads comes from energy flows directly collected at the connection node by the DSO. The power injections at MV/LV substation is provided by the aggregated power profiles determined in the LV planning tool in Section 1. As the PV installations are the property of industrial companies, these are also equipped with SM. This is not the case with most of the wind farms that often belong to independent operators that directly participate in the electricity market and do not depend on the DSO. Consequently, the power generated by wind farms is inferred from openly accessible wind speeds. As data for a large number of years are available, those are perfectly adapted to a Monte Carlo procedure. The wind speed values can then be converted into a value of active power according to the following power curve [34]:
where
Therefore, considering the theoretical number of the data available, the same pseudochronological Monte Carlo procedure as the aforementioned LV planning tool could be performed. However, several impediments such as the difficulty to gather the data of all clients combined with the noticeable loss of data and measurement errors led to the implementation of a nonsequential Monte Carlo.
Concerning the reactive power of DG installations, it must comply with the operating point of DG as well as the national legislation. This work is based on Belgian rules [3], but could easily be adapted in order to suit other legislation.
The reactive power compensators are modelled by their reactive power consumed (inductive mode) or injected (capacitive mode) into the grid. Then, given the price and the maturity of other storage technologies, the battery is currently the most appropriate solution for the desired application of providing ancillary services to the distribution grid. The pricing of these devices is determined by manufacturers, which will not be revealed here for confidentiality purposes.
As previously mentioned, the permitted range of voltage fluctuations is defined by the EN 50160 standard that allows a variation of ±10% around the nominal value. The thermal limits of power cables are defined as the maximum current that can flow in the conductors without inducing excessive heating. Indeed, the current flowing in the conductors generates heat that needs to be eliminated in the surrounding environment. This diffusion is made difficult by the important thermal resistances of the insulating materials in cables and by the surrounding environment, especially in the case of underground cables. In order to preserve the network security, the DSO must ensure that such situations do not occur.
3.5.2. Definition of the network and of the factors
The DOE method is used on a 28bus radial network, which is highly critical regarding the risks of congestion and voltage violation, because of the high penetration rate of DG units such as wind farms or large PV plants. These are installed at nodes 6, 10, 15, 23 and 28. This situation is indeed stimulated by the current financial attractiveness of investing in renewable energies. Consequently, in order to avoid an expensive reinforcement of the grid, the introduction of network regulation devices constitutes an interesting solution to investigate. The singlephase diagram of this 20kV system is shown in Figure 17. The representation of the grid is purposefully divided into eight different subareas, which are arbitrarily chosen and correspond to the areas whose influence is evaluated in the screening process. The resistance and reactance per kilometre of each line section are, respectively, 0.098 and 0.106 Ω/km, whereas the distance between nodes varies from 2 to 10 km. The definition of the rated power of each load and DG unit is shown in Table 1.






01  0  0  15  0  3900 
02  240  0  16  190  0 
03  190  0  17  160  0 
04  220  0  18  150  0 
05  100  0  19  230  0 
06  0  2700  20  100  0 
07  160  0  21  160  0 
08  100  0  22  140  0 
09  180  0  23  170  2700 
10  260  3900  24  110  0 
11  190  0  25  260  0 
12  130  0  26  170  0 
13  90  0  27  160  0 
14  180  0  28  140  3900 
3.5.3. The screening step
A preliminary analysis of the network behaviour showed that the congestion issues are exclusively observed in areas I and II of the grid, whereas the voltage issues tend to occur at the most distant nodes of the HV/MV substation. This gives useful information even before starting the screening procedure. However, the intuitive solution that consists in positioning one control device at the beginning of the grid and the other one at the extremity is likely to be inefficient. Indeed, the power flows originating from the transformer are representative of the power demand of the downstream loads. Therefore, the benefit of a control device installed in the neighbourhood of the substation is highly localized (i.e. area between the transformer and the control device) and this solution has a good probability to be suboptimal. The determination of the best solution is thus not trivial and requires a more thorough study.
Two screening studies are therefore carried out, both aiming at finding the optimal location to install respectively a reactive power compensator and a storage device. The influence of the eight areas represented in Figure 17 is investigated. For consistency, the same sizing is considered for all devices of the same type. In this way, the maximal sizing of the VAR compensators is adapted to the maximum reactive load of the case study and the maximal output power of the battery is 1 MW, which corresponds to the current upper limit proposed by the manufacturers.
The results of both screening studies are represented in Figure 18. The bar chart representation of the effects is privileged for its simplicity of interpretation. The height of a bar associated to a factor represents indeed the effect of this factor on the response. Therefore, the highest bars correspond to the most influential (and thus efficient) factors.
As expected, the results are quite similar for both scenarios. Indeed, the areas located near the end of the grid are more appropriate than those near the HV/MV transformers in the context of investing in network regulation devices for improving the MV grid operation. This can be explained by the importance of the line impedance. However, while the best place for the storage device is the end of the main feeder in area III, the VAR compensator will have more effect if it is installed near node 28 in area VIII.
3.5.4. The sizing step
The purpose of the sizing step is to use the network regulation devices in order to improve the operation of the grid for a large number of simulated network states. All the optimal values of these devices are then collected as they will constitute the basis for the sizing decision. Indeed, thanks to manufacturer offers, the output power can be associated to a price. The objective is to allow the DSO to make the best decision with regard to the investment costs and the resulting effects on the network operation.
For each network state, the centralized control strategy is carried out by using the aforementioned response surface methodology. The goodness of fit of the model (13) is quantified through the coefficient of determination
As previously mentioned, two scenarios are considered here. In the first one, the combined control of a VAR compensator with a storage device is studied while the second one focuses on the optimal regulation of two reactive power compensators. It should be noted that the range of variation of the reactive power compensator is symmetrical (
The general model of the polynomial approximation of the response for the first scenario becomes:
The optimal values of network storage devices for both defined scenarios are represented in Figure 19.
Before interpreting the results, it is important to focus on their accuracy by evaluating the quality of the model used to determine the optimal values of the control devices for each of the Monte Carlo simulation. To that end, the parameter
In the first scenario, as the current maturity of the battery technology does not allow having an output power higher than 1 MW, a significant part of the optimal value is close or equal to the maximum value. However, if the sizing of the VAR compensator is limited to 1.6 MVAR, more than 85% of the constraint violations are eliminated and this solution can therefore constitute a good compromise for the DSO. In the second scenario, it can be seen that a large part of the optimal point corresponds to opposite operating modes of the two VAR compensators (i.e. one is injecting reactive power while the other one is in inductive mode). This can be explained by the conflicting nature of the concurrent elimination of voltage and congestion issues. Indeed, as the line currents are fixed by the power demand, the removal of line congestions requires increasing the voltage. However, if overvoltages are already observed in other parts of the grid, the issue can hardly be solved and necessitates opposite actions of control devices. The nature of their actions then depends on their positioning with regard to the operational constraints violations.
3.5.5. The simulation of the realtime management of the grid
In this method, network regulation devices are used in order to ensure an optimal operation of the network. This objective is translated in the removal of all constraint violations of network operational conditions. The process consists therefore in a statebystate optimization based on the response surface methodology carried out for a large number of typical years. The effectiveness of the proposed methodology is then tested by taking as evaluation criteria the probabilities of voltage violation and the presence of congestions. The total line losses with and without the centralized control are also compared.
The simulation of the realtime management of the grid is performed for both defined scenarios with their predetermined sizing. In scenario 1, the battery has a maximum output power of 1 MW and the VAR compensator can vary between +1.6 and −1.6 MVAR. The second scenario simulates the centralized control of two reactive power compensators of 2 MVAR. Figures 20 and 21 show the comparison of constraint violations with and without the implemented control for the first and second scenario, respectively.
In both investigated variants, the undervoltage issues have been completely removed. However, the number of overvoltages has significantly decreased, but they could not be totally eradicated. The probability of congestion has even increased for some nodes, which highlights the difficulty to combine the objectives of simultaneously alleviating overvoltage and congestion issues. In such a case, the DSO would have no other choice than to have recourse to a curtailment of DG units, which may be very expensive. In this context, the research concerning the implementation of loadshifting strategies (e.g. with financial incentives in order to drive customers consuming when there is a lot of generation within the system) is of major importance.
The impact of control strategy on the line losses is represented in Table 2. The results are averaged to represent the total losses of a representative year. Recall that their minimization was not included in the objective function of the study since the reduction of the number of operating cycles of the control devices was privileged.



Without control  With control  Without control  With control  
Total line losses [MW]  632.89  822.73  665.30  843.07 
It can be seen that the amount of power losses increases with the network regulation devices. This can be explained by the presence of several DG units near the end of the grid. Indeed, in order to avoid the overvoltages inferred by the generation of these installations, the storage and the VAR compensators consumed locally active or reactive power. This energy was used by some loads that consequently need to be supplied by the HV/MV transformer. The power flows at the MV output of the substation are thus increased and, as the line losses are proportional to the square of the electrical current, these losses, together with the probability of congestion, rise significantly.
4. Conclusions and perspectives
The natural evolution of distribution networks with high integration of distributed renewable energy sources and random nature of consumptions loads, both uncertain and time variable parameters, requires the development of new accurate optimized methodologies for their analysis. To this end, probabilistic approaches are highly recommended since they can simulate distribution networks taking into account the time variation of renewable generation and consumption loads. This chapter introduced two probabilistic planning tools, both developed in a Monte Carlo algorithm environment, for the analysis and planning of LV and MV distribution networks. Section 1 presented a probabilistic framework that analyses in a fast, optimized way the voltage profile along LV feeders, taking into account the uncertainty of their loading parameters node by node, based on real SM data. It can therefore be deployed for the technoeconomic evaluation and refinement of solutions to operational problems in LV feeders with high PV penetration. The only prerequisite for its deployment is the availability of SM measurements in the studied LV feeder. This prerequisite goes along with the perspective of a wide rolling out of SM devices in the European LV networks. Section 2 introduced a new method intended to constitute an efficient and realistic planning tool for the DSO. In this way, the study takes three of the main objectives of the DSO into account. Indeed, the screening study aims at finding the best area for installing network regulation devices for minimizing the probability of violations of the operational constraints. Then, the sizing procedure is more complete by accounting for the investment costs. Finally, the centralized control of the number of operating cycles could also easily be implemented for reducing the line losses if such an action turns out to be financially worthwhile.
References
 1.
P. Lemaire, “Electricity market changes, impact on thermal generation”, SRBE/KBVE seminar, 24 Sept. 2014.  2.
EN50160, Voltage characteristics of electricity supplied by public electricity networks, 2012.  3.
SYNERGRID, Prescriptions techniques spécifiques de raccordement d’installations de production décentralisée fonctionnant en parallèle sur le réseau de distribution, C10/11, 2012.  4.
M.H.J. Bollen, F. Hassan, Integration of Distributed Generation Sources in the Power System, 1st ed. IEEE Press Series on Power Engineering, Hoboken, New Jersey: Wiley, 2011.  5.
EPIA, Global Market Outlook for Photovoltaics 20142018, 2014.  6.
S. Conti, S. Raiti, ”Probabilistic load flow using Monte Carlo techniques for distribution networks with photovoltaic generators”, Solar Energy, vol. 81, no. 12, 1473–1481, 2007.  7.
F.J. RuizRodriguez, J.C. Hernandez, F. Jurado, “Probabilistic load flow for radial distribution network with photovoltaic generators”, IET Renewable Power Generation, vol. 6, no. 2, 110–121, 2011.  8.
R. Billinton, R. Karki, “Reliability/cost implications of utilizing photovoltaics in small isolated power systems”, Reliability Engineering & System Safety, vol. 79, no. 1, 11–16, 2003.  9.
R.M. Moharil, P.S. Kulkarni, ”Reliability analysis of solar photovoltaic system using hourly mean solar radiation data”, Solar Energy, vol. 84, no. 4, 691–702, 2010.  10.
EPIA, Connecting the Sun: Solar Photovoltaic on the Road to LargeScale Grid Integration, 2012.  11.
F. Vallee, V. Klonari, T. Lisiecki, O. Durieux, F. Moiny, J. Lobry, “Development of a probabilistic tool using Monte Carlo simulation and smart meters measurements for the long term analysis of low voltage distribution grids with photovoltaic generation”, International Journal of Electrical Power & Energy Systems, vol. 53, 468–477, 2013.  12.
V. Klonari, F. Vallee, O. Durieux, Z. De Greve, J. Lobry, “Probabilistic modeling of short term fluctuations of photovoltaic power injection for the evaluation of overvoltage risk in low voltage grids”, EnergyCon, Croatia, 2014.  13.
V. Klonari, J.F. Toubeau, O. Durieux, Z. De Gréve, J. Lobry, F. Vallée, “Probabilistic analysis tool of the voltage profile in low voltage grids”, In Cired, France, 2015.  14.
J.F. Toubeau, V. Klonari, Z. De Gréve, J. Lobry, F. Vallée, “Probabilistic study of the impact on the network equipment of changing load profiles in modern low voltage grids”, accepted for presentation in ICREPQ’15, Spain, 2015.  15.
CEI 021, Regola tecnica di riferimento per la connessione di Utenti attivi e passivi alle reti BT delle imprese distributrici di energia elettrica, 2012.  16.
K. Turytsin, P. Sul, S. Backhaus, M. Chertkov, “Options for control of reactive power by distributed photovoltaic generator”, Proceedings of IEEE, vol. 99, no. 6, 1063–1073, 2011.  17.
B. Blazic, B. Uljanic, B. Bletterie, K. De Brabandere, C. Dierckxsens, T. Fawzy, W. Deprez, “Active and autonomous operation of networks with high PV penetration”, MetaPV FP7 project, 2011.  18.
C. Dierckxsens, K. De Brabandere, W. Deprez, “Economic evaluation of grid support with photovoltaics”, MetaPV FP7 project, 2012.  19.
J. AuYeung, G.M.A. Vanalme, J.M.A. Myrzik, P. Karaliolios, M. Bongaerts, J. Bozelie, W.L. Kling, “Development of a voltage and frequency control strategy for an autonomous LV network with distributed generators”, UPEC 2009, Scotland, 2009.  20.
H. Laaksonen, P. Saari, R. Komulainen, “Voltage and frequency control of inverter based weak LV network microgrid”, in FPS, The Netherlands, 2005.  21.
A. Engler, O. Osika, M. Barnes, N. Hatziargyriou, DB2 Evaluation of the local controller strategies, www.microgrids.eu/micro2000, 2005, accessed in June 2014.  22.
T.L. Vandoorn, J. De Kooning, B. Meersman, L. Vandevelde, “Voltagebased droop control of renewables to avoid onoff oscillations caused by overvoltages”, IEEE Transactions on Power Delivery, vol. 28, no. 2, 845–854, 2013.  23.
V. Klonari, J.F. Toubeau, T. L. Vandoorn, B. Meersman, De Grève, Zacharie, J. Lobry, and F. Vallée, “Probabilistic Framework for Evaluating Droop Control of Photovoltaic Inverters,” Electr. Power Syst. Res. , vol. 129, pp. 1–9.  24.
P.T. Staats, W.M. Grady, A. Arapostathis, R.S. Thallam, “A procedure for derating a substation transformer in the presence of widespread electric vehicle battery charging”, IEEE Transactions on Power Delivery, vol. 12, no. 4, Oct. 1997.  25.
D. Dohnal, “Onload tapchangers for power transformers”, Maschinenfabrik Reinhausen GmbH, 2013.  26.
J. Dixon, L. Moran, E. Rodriguez, R. Domke, “Reactive power compensation technologies: stateoftheart review”, IEEE Transactions on Industry Applications, vol. 93, no. 12, 2144–2164, Dec. 2005.  27.
C. Pieper, H. Rubel, “Revisiting energy storage. There is a business case”, Report of the Boston Consulting Group, Feb. 2011.  28.
V. Calderaro, G. Conio, V. Galdi, G. Massa, A. Piccolo, "Optimal decentralized voltage control for distribution systems with inverterbased distributed generators", IEEE Transactions on Power Systems, vol. 29, no. 1, 230, 241, Jan. 2014.  29.
E. Tzimas, I. Papaioannou, A. Purvins, "Demand shifting analysis at high penetration of distributed generation in low voltage grids", Electrical Power and Energy Systems, 2012.  30.
B. Bakhshideh Zad, J. Lobry, F. Vallée, “Coordinated control of onload tap changer and DSTATCOM for voltage regulation of radial distribution systems with DG units,” 3rd International Conference on Electric Power and Energy Conversion Systems, Istanbul, Turkey, October 2013.  31.
J.F. Toubeau, F. Vallée, Z. De Grève, J. Lobry, “Optimal allocation process of voltage control devices and operational management of the voltage in distribution systems using the experimental design method”, Proceedings of IEEE International Energy Conference (ENERGYCON), 1083–1090, 13–16 May 2014.  32.
D.C. Montgomery, Design and Analysis of Experiments, 5th ed. Hoboken, New York: Wiley, 2001.  33.
A.J. Wood, B.F. Wollenberg, Power Generation, Operation, and Control, New York: John Wiley & Sons, Inc., pp. 430–432, 1996.  34.
LB. Fang, JD. Cai, “Reliability assessment of microgrid using sequential Monte Carlo simulation”, Journal of Electronic Science & Technology, vol. 9, no. 1, 31–34, March 2011.