Next Article in Journal
Mapping Maize Cultivated Area Combining MODIS EVI Time Series and the Spatial Variations of Phenology over Huanghuaihai Plain
Next Article in Special Issue
Stochastic Predictive Energy Management of Multi-Microgrid Systems
Previous Article in Journal
Multi-Objective Design Optimization of the Reinforced Composite Roof in a Solar Vehicle
Previous Article in Special Issue
An Online Reactive Power-Optimization Strategy Based on Load-Curve Prediction and Segmentation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Location and Capacity Selection Method for Electric Thermal Storage Heating Equipment Connected to Distribution Network Considering Load Characteristics and Power Quality Management

1
Department of Electrical Engineering, College of Information Science and Engineering, Northeastern University, Shenyang 110819, China
2
State Grid Liaoning Electric Power Supply Co. Ltd., Shenyang 110003, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(8), 2666; https://doi.org/10.3390/app10082666
Submission received: 25 February 2020 / Revised: 2 April 2020 / Accepted: 9 April 2020 / Published: 12 April 2020
(This article belongs to the Special Issue Control, Optimization and Planning of Power Distribution Systems)

Abstract

:
High-permeability distributed wind power and photovoltaic systems are connected to the distribution network, which exacerbates the volatility and uncertainty of the distribution network. Furthermore, with the increasing demand of heating in winter and environmental protection, the wide use of electric thermal storage heating equipment (ETSHE) can promote distributed renewable energy utilization. However, an unplanned ETSHE connection to the distribution network may cause serious power quality problems. A new method of equipment location and capacity is proposed, which considered the improvement of power quality and load demand characteristics of the distribution network. First, based on heat load portrait technology, the node’s thermal load classification prediction was carried out to provide the data basis for the model solution. Second, the multi-objective optimal location and capacity programming model including harmonic distortion rate, voltage deviation, voltage fluctuation, and ETSHE cost was established. Then, the system nodes were preprocessed based on the sensitivity analysis method to reduce the number of installation nodes to be selected, and a feasible alternative set of installation nodes for the optimal configuration model of ETSHE could be obtained. Finally, the improved multi-objective particle swarm optimization algorithm was used to solve the model, and the data envelope analysis method was used to evaluate the power quality of each access scheme. The analysis of the numerical example shows that it can not only satisfy the user′s heat demand, but also effectively improve the power quality by rationally planning the location and capacity of ETSHE, which achieves the safe and efficient utilization of energy.

1. Introduction

In recent years, large-scale use of fossil fuels has led to exacerbating the environmental pollution problems, thus energy reform is imperative. New energy power generation is widely developed worldwide, among them, distributed photovoltaic and wind power generation are widely used because they have the characteristics of abundant available resources and reliable devices [1]. However, from another perspective, distributed photovoltaic and wind power also have the characteristics of fluctuation and stochastic. When the output of a distributed photovoltaic or wind power connected to a node in the power distribution system exceeds the node′s self-consumption, the node voltage will exceed the limit. At the same time, a large number of distributed photovoltaics use power electronic equipment, which will bring a lot of harmonics to the distribution network, and can cause serious power quality problems [2,3]. Due to the intensification of environmental problems and the oversupply of electrical energy, electric thermal storage heating equipment, which holds the function of time shifting and storage of electrical energy, has been widely applied. In fact, the location and capacity of ETSHE is closely related to the characteristics of heat load. The characteristics of heat load mainly include the duration of heat load and the distribution of heat load within the day. There are various factors that affect load such as population, weather, and holidays, etc. The relationship between these factors and load is non-linear, which is even more difficult to describe and predict. The difference of energy use characteristics of thermal load brings new challenges to distribution network planning. For different thermal loads, the optimal configuration planning of ETSHE will also change accordingly, therefore, load portrait is the premise and foundation of optimal operation [4].
In light of the improvement of power quality, scholars at home and abroad have found that the goal of improving power quality can be achieved through the reasonable allocation of energy storage devices at corresponding nodes. However, traditional storage battery energy storage equipment, with poor peak-regulating capacity and high cost, cannot meet the needs of heat users. However, the storage of energy through the ETSHE can make good use of its time-shifting and storage characteristics to meet the multi-energy needs of users. However, the different access locations and capacity of the energy storage systems have a great influence on the power quality governance, so the choice of access location and the determination of capacity are problems that need to be solved urgently [5]. In [6], a constant capacity optimization model for energy storage equipment location was presented and solved by genetic algorithm. Meanwhile, the total cost was saved by adding the capacity correction link. In [7], a two-layer optimization method for multi-point distributed energy storage system location and fixed capacity in high photovoltaic penetration distribution networks was proposed. The authors in [6,7] only considered traditional electric energy storage. In addition to the basic energy storage planning presented above, some have also studied the energy storage planning method considering the demand characteristics. Such as [8], the proposed polygeneration microgrids system includes trigeneration plants, renewables, district heating, thermal and electrical storage systems and flexible loads, building demand-response model and storage systems model. Thermal storage proved to potentially have a very relevant impact on the microgrid operation, and the analysis also showed the potential reduction of the operating costs that could be achieved by demand-response. Furthermore, the authors in [9,10] realized the development of an energy management system (EMS), which was implemented and run at the smart polygeneration microgrid (SPM), that modified the existing EMS model and predictive control to take into account the existence of a thermal storage system with a suitable capacity. It is worth noting that the heat storage has been represented by resorting to a new battery-type model that is able to predict the transient behavior of the grid, according to the historical series of heat and electricity demand data during 2017. In such a way, it has been possible to calculate the optimal size of the thermal storage that maximizes the global saving. The study in [11] investigated the influence of energy storage location and capacity on the overall economy of the grid and wind energy utilization, where the premise is a third party owns the energy storage. However, this did not consider the improvement of energy quality caused by energy storage. In [12,13], according to the data collected by the power quality detection device, an appropriate evaluation method to calculate the power quality level of the access point after the energy storage system was connected. However, it did not combine the planning of the energy storage system with a power quality assessment.
In summary, considering the power quality and the impact of load characteristics, there are still few works concerning solving the location and capacity selection problem for ETSHE (not basic electric energy storage). Based on the above discussions, this paper proposed a method for selecting the location and capacity of the ETSHE connected to the distribution network, which avoids the disadvantages of high costs and the existence of pollution. The adaptive inertial weight improved particle swarm optimization (IMOPSO) method was proposed, which is suitable for distribution network systems with larger nodes. The influence of load characteristics on the planning was considered, thus the research results were closer to the actual situation. The model uses the power quality as the optimization target, which effectively solves the power quality problem caused by the new energy generation with high permeability.
The main contributions are introduced from four aspects, as follows:
  • First, the minimum intra-class distance is used as the objective function to perform fuzzy clustering on the heat load. Second, on the basis of conventional forecasts, the influence weight of the meteorological factors on the outdoor temperature of the day is calculated to correct the outdoor temperature, improve the linear fitting goodness with heating heat load, and provide initial classified heat load data for the optimal configuration of ETSHE.
  • The sensitivity coefficient is calculated based on the sensitivity radar chart and sensitivity calculation formula. The sensitivity discrimination of each node in the distribution network can effectively eliminate the nodes that do not meet the sensitivity and prevents the solution from falling into the wrong node, which can improve the efficiency and accuracy of the solution.
  • Establish a location and capacity planning model for ETSHE and multi-objective functions including harmonic distortion, voltage deviation, voltage fluctuation, and energy storage equipment costs.
  • In the planning stage, data envelopment analysis (DEA) is used to evaluate the power quality of access schemes with multiple indicators and single indicators, and an effective and unified method for distribution network planning and power quality assessment is proposed.

2. Load Characteristic Portrait

Traditional load forecasting methods generally only consider the typical daily curve of heat load, but not consider the differences in the heating characteristics of classified heat loads. Under the background of the rapid development of big data technology, portrait technology has been widely used in the fields of e-commerce and the internet [14]. This paper applied portrait technology to the extraction of heating user characteristics for the first time. In this paper, the thermal user load portrait was defined as follows. Different thermal loads were predicted by a high-precision prediction method so that the thermal loads were treated differently according to the different heating characteristics. To comprehensively consider the user heating characteristics and influencing factors, we established a user heating behavior label library. Based on the fuzzy C-means (FCM) algorithm for label clustering, different types of heating users used warming mode portraits, which makes the results of user heating behavior analysis more intuitive [15].

2.1. Clustering with Warm Behavior Patterns

The goal of users label clustering is to quickly grasp the heating characteristics of different thermal users to achieve different treatments of different heating groups. After combing and screening many times, the following six labels were obtained, as shown in Table 1.
In this paper, FCM clustering algorithm was used to calculate each thermal user’s label data to be clustered, and four memberships of thermal user clustering centers to realize the clustering division of the sample thermal user’s label data. During the iterative process, FCM takes the minimum intra-cluster distance and the maximum distance between clusters as the optimization goals, and constantly adjusts the categories of various thermal user data, so that it optimizes the clustering results continuously. The object to be clustered is the D -dimensional dataset of N samples:
X = [ x 11 x 1 j x 1 D x i 1 x i j x i D x N 1 x N j x N D ]
The data collection class is divided into K subsets, where the value of K is 4, the subset is represented as C = { c 1 , , c i , , c 4 } ; and c 1 ,   c 2 ,   c 3 , and c 4 represent residents, education, business and industry, respectively. Here, we set D = 7, which represents six items of user label data and one category data, respectively. Assuming that the probability of each thermal user j belonging to the i th clustering subset is u i j , then FCM can be defined as an optimization problem with equality constraints [15]:
{ m i n   J = i = 1 K j = 1 N u i j m x j c i 2 s t . i = 1 K u i j = 1
where J is the value function of FCM. m is the membership factor, which is used to adjust the degree of influence of the membership indicator in the optimization target. In most cases, we chose m = 2. The iterative formulas for obtaining membership coefficients and cluster centers are shown in Equations (2) and (3) [15]:
c i = j = 1 N ( x j u i j m ) j = 1 N u i j m
u i j = 1 k = 1 K ( | | x j c i | | | | x j c k | | ) 2 m 1
Clustering steps of the warm behavior patterns for thermal users based on the FCM algorithm:
  • Initialization: Set the number of clustering K and membership factor m, randomly initialize the membership matrix U , U is a matrix made up of u i j , the center of the cluster.
  • Parameter update: According to Equation (2), calculate the clustering center C of each label of thermal users; according to Equation (3), calculate the membership matrix U of all thermal user labels and each clustering center C , classify thermal users into the clustering center with the highest membership.
  • Convergence judgment: If the change of the objective function is less than the preset threshold, then output the clustering result, otherwise return to Step 2.
  • End of clustering: Save and display the clustering results.

2.2. Load Portrait

Load portrait refers to the method of high-precision prediction of the heating heat load by considering the main factors affecting heating load such as outdoor temperature, solar radiation, natural wind, etc. Correct outdoor temperature can be performed by using solar radiation and wind speed to improve the fit between outdoor temperature and thermal load [16]. Wind speed and solar radiation are the main factors affecting the heating heat load. The correction of wind speed for outdoor temperature can be converted into an equivalent cooling temperature [17]. The equivalent cooling temperature formula considering the influence of wind speed is shown in Equation (4):
Δ T w i n d = α 1 ( l g ( β W w i n d ) ) 3 α 2 ( l g ( β W w i n d ) ) 2 + α 3 l g ( β ( W w i n d ) )
where W w i n d is the external wind speed value, m / s ; Δ T w i n d is the equivalent temperature of wind speed cooling; °C. α 1 ,   α 2 ,   α 3 , and β are the wind speed correction coefficients, where the values are 0.0246, 0.4525, 3.2398, and 7.23, respectively.
The correction of outdoor temperature by solar radiation can be converted into an equivalent heating temperature as shown in Equation (5):
Δ T s o l a r = S s o l a r / S k
where S s o l a r is the external light value; W / m 2 , S k is the light conversion coefficient; W / ( m 2 ° C ) ; according to [17], we set the unit of S k as 100 W / ( m 2   ° C ) .
T o u t is the uncorrected outdoor temperature. When considering both the outdoor wind speed and solar radiation of the day, the calculation formula for outdoor temperature T w , s is given as follows:
T w , s = T o u t + Δ T s o l a r Δ T w i n d
In order to improve the accuracy of forecasting heating load, this study not only considered the influence of wind speed and solar radiation on outdoor temperature of the day, but also considered the influence of wind speed and solar radiation on outdoor temperature of the past day. In this paper, we calculated the influence weight of previous meteorological factors on the outdoor temperature of the day to modify the outdoor temperature and linearly fit it with the heating heat load, so that their fitting goodness was further improved [18,19,20].
According to heat transfer knowledge, the influence of meteorological factors on the outdoor temperature of the day decreases with the increase in days. In order to facilitate the prediction of heating heat load, the influence of meteorological factors of all past days on the outdoor temperature of that day cannot be considered. Statistical results showed that heating heat load was related to meteorological factors of the previous 4–5 days. We used MATLAB programming to find the optimal days and the reasonable weight of past weather factors to correct the outdoor temperature of the day. When the outdoor temperature and the heating heat load were linearly fitted, the goodness of fit value was the maximum. This means that when the number of days considered is n days, the outdoor temperature of the day is given as follows:
T w , i , e = A i T w , i + A ( i 1 ) T w , ( i 1 ) + A ( i 2 ) T w , ( i 2 ) + + A ( i n ) T w , ( i n )
where A i , A ( i 1 ) , A ( i 2 ) , represent the weight value of day i , day i − 1, day i − 2, and day i n, respectively, and A i > A ( i 1 ) > A ( i 2 ) > …… > A ( i n ) , and the sum of weight values is 1. T w , ( i n ) is the outdoor temperature after solar radiation and wind speed correction of day i − n, and T w , i , e is the outdoor temperature of day i after correction on the current day and the previous day.
The “trial” method was adopted to determine the index weight of the influence of past meteorological factors on the outdoor temperature of the day to correct the outdoor temperature and obtain the relationship between the heat load and the temperature of the day, then obtained the relationship between the heat load and time. For example, if n = 2 in Equation (7), we can set A i = m 0 / ( m 0 + m 1 + m 2 ) , A ( i 1 ) = m 1 / ( m 0 + m 1 + m 2 ) , and A ( i 2 ) = m 2 / ( m 0 + m 1 + m 2 ) to simplify the calculation, and m 0 ,   m 1 ,   m 2 is an integer between 1 and 10.

3. Optimal Configuration Model of Electric Thermal Storage Heating Equipment

High-permeability distributed photovoltaic and wind power connected to the distribution network cause a series of power quality problems, which can be improved by installing ETSHE at appropriate nodes of the distribution network. However, the power quality of the power grid involves many aspects such as voltage and current. The location and capacity of the energy storage system needs to consider the effect of improving the power quality, the consumption of new energy, and the economics of the energy storage system. It is a multi-objective optimization plan model.

3.1. Objective Function

3.1.1. Total Current Harmonic Distortion

The fundamental current flowing through the power load is equal to the sum of the fundamental current of the power grid and the fundamental current of the inverter. When the power load is unchanged, the current harmonics are unchanged, but the fundamental current of the distribution network is reduced, and the total current harmonic distortion rate at the access point is increased. With the increase of access capacity, the total current harmonic distortion rate of the distribution network may exceed the range of grid access energy. Therefore, the total distortion rate of current harmonics was selected as one of the objective functions of the energy storage system location and capacity. The mathematical expression is [21]:
f T H D I a = h = 2 M ( I h ) 2 I 1 × 100 %
where I 1 is the root mean square value of the fundamental current and I h is the root mean square value of the h -th harmonic current. Here, we set M = 7, where only the 7th and lower harmonics were considered.

3.1.2. Voltage Deviation

Voltage is one of the most important indicators of power quality, and voltage deviation is a main indicator to measure the normal operation of the power supply system [22]. Distributed photovoltaic access to the distribution network can increase the voltage of the nodes, and excessive photovoltaic output will cause some nodes to exceed the limit. In this paper, the node voltage deviation was selected as one of the objective functions of the energy storage system location and capacity determination, its mathematical expression is [13]:
f δ U = i = 1 N j = 1 T | U i j U N |
where U i j is the voltage of node i at time j ; U N is the system standard voltage; and N is the number of system nodes. T is the number of detection moments.

3.1.3. Voltage Fluctuation

As the output power of the device is fluctuating, the voltage fluctuation caused by the connection to the distribution network is Δ d , so the voltage fluctuation was selected as one of the objective functions of this article. The mathematical expression is [23]:
f Δ d = Δ I Z U N × 100 %
where U N is the rated voltage; Δ I is the current change caused by the output power change at the access point; and Z is the equivalent impedance of the two-port network distribution network viewed at the access point.

3.1.4. Total Cost of Energy Storage System

The ETSHE must also consider the cost issue while ensuring the quality of the power. The total cost of the energy storage system was selected as one of the objective functions.

Equipment Cost of Electric Boiler

The equipment cost of an electric boiler is proportional to its rated heating capacity, which is equal to the design daily heat supply rating, thus Equation (11) is:
C S ( Q S e ( Q S ) ) = ρ S Q S e ( Q S ) = ρ S max t = 1 , , 24 { Q S ( t ) }
where ρ S is the unit price corresponding to the unit heat supply of the electric boiler. Q S e ( Q S ) is the rated heat supply of the electric boiler. Q S = [ Q S ( 1 )    Q S ( 2 )    Q S ( 24 ) ] T is a vector composed of the hourly heating capacity of the designed electric power boiler.

Equipment Cost of Heat Storage Body

The cost of the heat storage equipment is proportional to the volume of the equipment, then C X ( V X ) = ρ X V X , and ρ X ,   V X are the unit price and volume of the heat storage body, respectively.

Equipment Cost of Heat Storage Pump

The key parameters of heat pump selection are head and flow. The head is only related to the building structure and was selected according to typical values, and the heat pump investment was proportional to the rated flow at a certain head [24,25]. Therefore, the heat pump equipment cost can be expressed by Equation (12):
C B ( F B e ( Q S , θ ) ) = ρ B F B e ( Q S , θ ) = ρ B γ C 1 max t = 1 , , 24 { Q S ( t ) ( 1 θ ( t ) ) } C W ρ W Δ τ T
where ρ B is the increased equipment investment for each additional unit of the rated flow of the heat pump. γ is the safety factor when calculating the rated flow, in this paper, γ = 1.05 .   C 1 is the conversion constant of kilowatt hours and joules, here C 1 = 3.6 × 10 6 . θ = [ θ ( 1 )    θ ( 2 )    θ ( 24 ) ] T is a vector formed by designing the hourly direct supply ratio of the electric power boiler, and C W ρ W are the specific heat and density of the heat medium. The difference in heat storage temperature is usually Δ τ T = 40   ° C [26].

Equipment Cost of Control Device

The equipment cost of the control device is related to the complexity of the electric boiler system. It was found that the equipment cost of the control device is basically proportional to the rated heating capacity of the heating device of the electric boiler, which is shown as follows:
C K ( Q S e ( Q S ) ) = ρ C Q S e ( Q S )
where ρ C is the investment cost of the control device for each additional 1 kW of the rated heating capacity of the electric boiler.
This article does not consider the operation and maintenance costs of energy storage equipment, only the investment and installation costs of energy storage equipment were considered, so the total cost of ETSHE is given as Equation (14):
f C T ( Q S , V X , θ ) = k i n s [ C S ( Q S e ( Q S ) ) + C X ( V X ) + 2 C B ( F B e ( Q S , θ ) ) ] + C K ( Q S e ( Q S ) )
where k i n s is the equipment installation cost coefficient. Usually, the equipment installation cost is about 20% of the total equipment cost, so we set k i n s = 1.2.
Therefore, the location and capacity multi-objective optimization model of the energy storage system comprehensively considered the voltage deviation of the distribution network nodes, the total distortion rate of the current harmonics, the voltage fluctuation, and the total cost of the energy storage system. The objective function is shown as Equation (15):
{ min   F = ϵ f T H D I a + ϑ f δ U + τ f Δ d + ς f C T ( Q S , V X , θ ) ϵ + ϑ + τ + ς = 1
where ϵ ,   ϑ ,   τ ,   ς is the normalized weight coefficient.

3.2. Restrictions

3.2.1. Equality Restrictions

1. Power flow constraint equation
The initial power flow calculation of the planned distribution network is performed, the voltage value and active power of each node are calculated. The power flow constraint is as follows [27]:
P i = U i j i U j ( G i j c o s θ i j + B i j s i n θ i j ) Q i = U i j i U j ( G i j s i n θ i j B i j c o s θ i j )
where P i and Q i are the active and reactive power injections at node i . U i and U j are the voltage amplitudes of nodes i and j . G i j and B i j are the conductance and susceptance of branch i j . θ i j is the voltage phase angle difference between the i -th and j -th node.
2. The hourly balance of heat supply and heat demand is:
θ ( t ) Q S ( t ) + Q q ( t ) = Q l l , m a x
where θ = [ θ ( 1 )    θ ( 2 )   θ ( 24 ) ] T is a vector formed by the design of the hourly direct supply ratio of the electric power boiler. Q S = [ Q S ( 1 )    Q S ( 2 )    Q S ( 24 ) ] T is the vector for designing the hourly heating capacity of the electric power boiler. Q q ( t ) is the heat provided by the heat storage body to the user during period t . Q l , m a x is the maximum design daily heat load. r l ( t ) is the ratio of the heat load to the maximum value for the t days on the design day.
3. The thermal equilibrium equation of the heat storage body is:
Q N ( t + 1 ) = Q N ( t ) + Q S ( t ) [ 1 θ ( t ) ] Q q ( t ) t = 1 , 2 , , 23 Q N ( 1 ) = Q N ( 24 ) + Q S ( 24 ) [ 1 θ ( 24 ) ] Q q ( 24 )
where Q N ( t ) is the amount of heat stored in the heat storage body at the beginning of time t.

3.2.2. Inequality Constraint

1. The maximum heat storage capacity is limited by the volume of the heat storage body:
0     Q N ( t )     V X Δ τ T ρ W C W υ X C 2 t = 1 , 2 , , 24
where V X is the volume of the heat storage body, and C W and ρ W are the specific heat and density of the heat medium, respectively. C 2 is the conversion constant of kilowatt-hours and joules, here C 2 = 3.6 × 10 6 . Usually, the difference in heat storage temperature is Δ τ T = 40 °C. υ X is the effective volume ratio of the heat storage body, which refers to the ratio of the volume that can actually contain the heat medium to the total external volume.
2. The upper limit of the heat supply of the heat storage body is:
Q q ( t )     Q N ( t ) t = 1 , 2 , , 24
3. The range of the value of the direct heat supply ratio of electric boilers is constrained as:
0     θ ( t )     1 t = 1 , 2 , , 24
4. Constraints of voltage levels before and after the ETSHE is connected to the distribution network:
L D G     ε 1 × L d i s
where L D G and L d i s are the voltage level indicators before and after the energy storage is connected to the distribution network, respectively; ε 1 is between 0 and 1.
5. Constraints on active loss before and after the ETSHE is connected to the distribution network:
P l o s s , D G     ε 2 × P l o s s , d i s
where P l o s s , D G and P l o s s , d i s are the active power loss before and after the ETSHE is connected to the distribution network, respectively. Similarly, ε 2 is between 0 and 1.
The objective functions and constraints of the optimization model in this paper were grouped as shown in Figure 1:

4. Node Sensitivity Analysis and Calculation Method

For power networks with a small number of nodes, the optimal configuration scheme for the location and capacity of ETSHE can be easily calculated by exhaustive methods. When the power network is complicated, there are many installation nodes available for installing energy storage at this time. If the optimal node is still selected according to the exhaustive method, the calculation of the optimal allocation scheme is large, which cannot be solved or the result falls into a local optimum. For this reason, the sensitivity analysis method was proposed to calculate the sensitivity coefficient of electric regenerative heating equipment installed at the node, and the node with higher sensitivity was selected as the alternative set of the installation node of ETHSE, which will greatly improve the optimal allocation efficiency of electrical thermal storage.
Sensitivity coefficient refers to the degree of sensitivity of the system, represented by state variables to changes in the control variables and disturbance variables [28]. For the electric storage thermal heating equipment location and capacity optimization configuration model, it refers to the sensitivity of the objective function of the optimal configuration to changes in energy storage installation capacity. The optimization objective of this paper is multi-objective function and the fundamental purpose of installing energy storage is to effectively control the power quality problem of the distribution system. When a large number of single-phase photovoltaic power sources are connected to the system, it will cause voltage three-phase imbalance problems. However, as single-phase photovoltaic power generation has not been widely used, it has little influence on the three-phase unbalance degree of power grid voltage. Therefore, only the three-phase imbalance is considered when defining sensitivity.
In a three-phase three-wire system, there is no zero-sequence component. For a three-phase system without a zero-sequence component, the simplified calculation method of the three-phase imbalance recommended by the national standard is as follows [29]:
f ε 2 = 1 3 6 L 1 + 3 6 L × 100 %
L = a 4 + b 4 + c 4 ( a 2 + b 2 + c 2 ) 2
where f ε 2 represents the three-phase negative sequence imbalance, and a , b , and c represent the effective value of three-phase voltage, respectively.
In this paper, the radar chart analysis method was used to visually and comprehensively analyze and evaluate the sensitivity of the nodes. Finally, the nodes with higher sensitivity were selected as candidate nodes. The schematic diagram of the sensitivity radar is shown in Figure 2. The solid arrows in the figure indicate the power quality indicators; the value increases according to the arrow direction. The outer circle indicates the advanced level value of each indicator, which can be expressed by 1.5 times of the indicator standard value. The middle circle indicates the standard value of each indicator. The inner circle indicates the unqualified value. The dashed line indicates the actual value of each indicator. If the actual value is less than the standard value, this type of indicator is indicated by “−”. If the actual value is greater than the standard value, this type of indicator is indicated by “+”.
Therefore, in this paper, when the four power quality indicators of any node were defined as “−”, it was considered that the sensitivity of this node to power quality reached the standard and could be used as an alternative node. The specific calculation method is:
α s e n = { 1 ,   Δ f Y H D I a , Δ f δ U , Δ f Δ d , Δ f ε 2 s t a n d a r d   v a l u e 0 , o t h e r
When the sensitivity coefficient α s e n = 1, the node can be used as a candidate node. If α s e n = 0, the node is eliminated.
The sensitivity coefficient α s e n reflects the degree of sensitivity of any node to the four aspects of harmonic content, voltage deviation, voltage fluctuation, and three-phase imbalance. If the energy storage system is installed at this node, the higher the value, the better the power quality will be. Therefore, by calculating the sensitivity coefficients of all nodes in the system, it is possible to find in advance a set of alternative nodes that are better at improving the power quality of the system, which can effectively reduce the number of computing nodes, greatly reduce the amount of calculations, and improve the optimization of the economy and efficiency of configuration.
Nodes are ranked according to the sensitivity coefficient, and higher-sensitivity nodes are preferentially selected to form an initial optimized candidate set. Subsequently, the heuristic optimization algorithm can be applied to solve the model quickly. This paper used the adaptive inertial weight IMOPSO algorithm to solve the model.
The main steps of the adaptive inertia weight IMOPSO algorithm are as follows:
  • The network topology, switch state, node voltage range, distributed photovoltaic, wind parameters, and parameters of ETSHE in the distribution network system are initialized. The sensitivity coefficients of each node are calculated and sorted separately. The top 40% of the sorted nodes are selected to form a candidate set of installation nodes.
  • The parameters such as the size of the particle swarm, the maximum number of iterations of the algorithm, the acceleration factor, the value range of the adaptive inertia weight, and the speed range of the particles are set in the improved particle swarm optimization.
  • The particles in the candidate set are evaluated, and the fitness function value of each particle is calculated. After ranking, the individual optimal value and the global optimal value of the particles are determined and updated.
  • Update the particle′s velocity and displacement, calculate the value of the adaptive inertia parameter.
  • Check if the maximum number of iterations is exceeded. If it is exceeded, stop iteration and output the calculation result, otherwise skip to Step 3 to continue the iterative calculation.
In this algorithm program, the total current harmonic distortion rate, the sum of voltage deviations, the voltage fluctuations, the investment, and installation cost of the energy storage system are optimized objective functions. The number of iterations is the cycle control variable, and the location and capacity of the ETSHE are taken as decision variables. By improving the solution of the particle swarm algorithm, the optimal position and optimal capacity for the energy storage installation can be obtained.
The solution flow chart of IMOPSO algorithm in this paper is shown in Figure 3.

5. Simulation Analysis

This paper used IEEE-33 and IEEE-118 power distribution systems for verification, respectively, which were used to analyze the impact of high-permeability wind power and photovoltaic access on the power quality of the power distribution network. Through the rational optimization of flexible resources such as ETSHE, the power quality of high-permeability wind power and photovoltaic power distribution networks is improved.

5.1. Load Portrait Result

Using FCM clustering algorithm for label clustering, where the value of N was 200 in this paper, we obtained four typical thermal users, as shown in Table 2. Figure 4 shows the distribution diagram of the FCM clustering results. According to the membership matrix in Figure 4, there were some differences in the characteristics of the four types of heating customers, but there were also some overlaps, which will affect the planning of the ETSHE. Equations (27) and (28) were used to calculate the fitting relationship between T w , i , e and Q , and the four types of typical users’ heat consumption over time were further calculated, as shown in Figure 5. The results of the thermal load portrait were input into the PSO as the initial data.
It can be seen from the results that when considering the influence of meteorological factors on the outdoor temperature of the day in the past two days, the goodness-of-fit R 2 of the linear fitting with the heating heat load was a maximum value after the correction of the meteorological factors on that day. It shows that the curve at this time can best reflect the relationship between Q and T , and the value of Q can be predicted from the modified T . For this set, the heating heat load can be predicted by the following formula:
T w , i , e = 0.7 T w , i + 0.2 T w , ( i 1 ) + 0.07 T w , ( i 2 )
Q = 187.6 10.13 T w , i , e
The predicted values of different heating heat loads were obtained from Equations (27) and (28). Within a day, the temperature changes with time, and the heating heat load changes with air temperature, so the curve of the heating load with time can be obtained equivalently. The energy consumption characteristics of four typical users are shown in Figure 5, where the typical users′ energy consumption characteristic curves show different trends.

5.2. Simulation Scenario Setting

We used two scenarios to analyze the effects of the optimal scheduling of electrical thermal storage heating equipment and load demand characteristics on the improvement in the power quality of the distribution network.
Scenario 1: In the absence of flexible resources, the indicators of power quality include high-permeability wind power and photovoltaic power distribution networks.
Scenario 2: In the case of ETSHE and load demand characteristics affecting the optimal dispatching scenario, various indicators of power quality include high-permeability wind power and photovoltaic distribution networks.

5.3. IEEE-33 System Simulation Analysis

1. Basic data
The distribution network load, photovoltaic, and wind power data are shown in Table 3. The installation locations of wind power and photovoltaic are shown in Figure 6. The distribution network contains twelve distributed photovoltaics with a capacity of 0.5 MW, which are connected to nodes 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, and 24. At the same time, the distribution network also contained four distributed wind power generation with a capacity of 0.5 MW, which were connected to nodes 26, 28, 30, and 32. In order to ensure the maximum consumption of renewable energy, each distributed power supply is operated with actual output.
2. Result analysis
(1) Optimal scheduling results of ETSHE considering load characteristics
In principle, the nodes that can be connected to ETSHE are 1–32, but harmonics, voltage deviations, and voltage fluctuations caused by high-permeability distributed wind power and photovoltaic will cause some nodes to cancel each other out. Considering the objective function, the node sensitivity analysis method was applied, and the set of nodes for energy storage access included five nodes: nodes 4, 8, 13, 25, and 31. In order to highlight the impact of load characteristics on optimal planning settings, we considered that nodes 4 and 13 were residential loads, node 8 was a commercial load, node 25 was an industrial load, and node 31 was an educational load.
According to the optimization model constructed in this paper, the IMOPSO algorithm was used to solve the calculation example. The population size of the IMOPSO algorithm was 80, and μ m i n = 0.5, μ m a x = 0.8, σ = 0.2 [30,31]. The installation position and capacity of the electric regenerative heating equipment are shown in Table 4.
(2) Node sensitivity analysis before and after comparison
As can be seen from Table 5, the range of optimized nodes is effectively reduced and the effectiveness of the solution is greatly improved by the use of the node sensitivity analysis method.
(3) Power quality improvement analysis
Through the optimization of two flexible resources, the power quality of distributed wind power with high penetration and photovoltaic power distribution networks has been significantly improved. Figure 7, Figure 8 and Figure 9 are the comparison curves of the total current harmonic distortion rate, voltage deviation, voltage fluctuation of the distribution network with high penetration wind power, and photovoltaic access in two scenarios. Line charts show the comparison of different scenarios of node power quality indicators, and the bar charts show the percentage improvement of the two scenarios relative to the national standard value.
According to the GB/T14549 standard, the total current distortion rate of the electric heat storage heating equipment after being connected to the distribution network should be less than or equal to 2.2% [21,32]. Figure 7a,b are the comparison curves of the total distortion rate of the node current harmonics in the two scenarios, and Figure 7c is the percentage increase of the harmonic distortion rate of the node current relative to the national standard value. As can be seen from Figure 7a,b, in Scenario 1, nodes 4 and 31 have a high light intensity between 8:00 am and 7:00 pm. Due to the introduction of a large number of power electronic equipment by photovoltaics, the current harmonic distortion rate exceeded the standard. To ensure the stability of the system, wind abandonment, light abandonment, and load shedding must occur. In Scenario 2, the current harmonic distortion rate was successfully reduced below the standard value, and the improvement rates were 45% and 47%, respectively. Although nodes 8 and 25 did not exceed the standard, the percentage of improvement through energy storage optimization can reach up to 50% and 26%. As can be seen from Figure 7c, the percentage of improvement of power quality relative to the standard value at each time was a positive value, which means that the power quality had improved. From 8 am to 7 pm, ETSHE mainly used the abandoned light resources, and other times mainly used abandoned wind resources.
According to GB/T12325, the absolute value of the voltage deviation should be less than or equal to 7% [33]. As can be seen from Figure 8, taking into account the load characteristics and the reasonable planning of the ETSHE, the voltage deviation of each node was significantly reduced and the voltage deviation of the access point at each moment was within the allowed range of the national standard. As can be seen from Figure 8a,b, at 12:00 am, the light intensity was the highest, but the corresponding voltage deviation was the lowest. The voltage deviation was inversely proportional to both wind speed and light intensity. Furthermore, as can be seen from Figure 8c, the percentage of each node′s improvement relative to the standard value was positive, with an average improvement rate of 38%, and the voltage deviation of each node was improved.
According to the GB/T12326-2008 standard, the voltage fluctuation of a 10 kV distribution network should be less than or equal to 4% [23]. As can be seen from Figure 9a,b, through the optimized configuration of energy storage, the node voltage fluctuation was improved. At 12:00 am every day, the light intensity was the highest, but the voltage fluctuation is the lowest. It can be seen from Figure 9c that at the time when there was more wind power and photovoltaic power generation, the voltage fluctuation was greatly improved relative to the standard value. Among them, node 4 increased by 23% on average, nodes 8 and 25 increased by about 25%, and node 31 increased by 11.2% on average.

5.4. IEEE-118 System Simulation Analysis

1. Basic data
The standard IEEE-118 distribution system was adopted for simulation analysis [34]. The distribution network contained 19 distributed photovoltaics with a capacity of 0.5 MW, which were connected to nodes 7, 9, 14, 31, 35, 36, 52, 54, 59, 60, 63, 78, 85, 87, 95, 96, 113, 114, and 117. The distribution network also contained 10 distributed wind power generation with a capacity of 0.5 MW, which were connected to nodes 4, 18, 25, 26, 28, 30, 32, 40, 41, and 42. The rest of the data is consistent with the IEEE-33 system, which is not repeated here.
2. Results analysis
(1) Optimal scheduling results of ETSHE considering the load characteristics
The same method was used for node sensitivity analysis, and the set of nodes that were selected to obtain the energy storage access included {3, 6, 10, 35, 38, 44, 47, 48, 61, 66, 77, 78, 79, 91, 92, 93, 114, 115, 117}, a total of 19 nodes. In order to highlight the impact of load characteristics on the planning settings, it was considered that nodes 3, 6, 10, 61, 66, 77, 91, 92, and 93 were residential loads; nodes 35, 38, 78, and 79 were commercial loads; nodes 44, 47, 48 were industrial loads; and nodes 114, 117, and 119 were educational loads.
According to the optimization model constructed in this paper, the IMOPSO algorithm was used to solve the calculation example. The installation position and capacity of the electric regenerative heating equipment are shown in Table 6.
(2) Node sensitivity analysis before and after comparison
This system is more complex than IEEE-33, so it takes a longer time to calculate. The comparative analysis before and after the node sensitivity analysis is shown in Table 7.
(3) Power quality improvement analysis
From the above simulation results, it can be seen that the power quality of residential and educational loads was similarly affected by new energy generation, while commercial and industrial loads were similar. Therefore, the residential and industrial loads were selected to give the results for node 3 and node 44, respectively.
It can be seen from the simulation results that although the IEEE-118 system had a long and complex calculation time, the node sensitivity calculation could obtain a set of candidate nodes in advance, which improved the efficiency and accuracy of the solution. The same results as those of the IEEE-33 nodes’ power distribution system can be seen in Figure 10, Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15. Through the simulation analysis of the above IEEE-33 and IEEE-118 systems, by configuring a certain capacity of ETHSE at specific nodes of the distribution network, it can not only solve the voltage deviations and voltage fluctuations caused by the photovoltaic and wind power, but can also restrain harmonics, which can improve the power quality of distribution network effectively. Finally, it can be seen that the evaluation levels were all above five levels by using the data envelope analysis method, which indicates that the access scheme is feasible and that the indicators of the access point can reach the prescribed standards. For specific evaluation methods and evaluation levels, see Appendix A.

6. Conclusions

Aiming at the power quality problems caused by high-permeability distributed photovoltaic and wind power connected to the grid, this paper proposed a method for selecting and determining the capacity of the ETSHE connected to the distribution network, which considered the load characteristics and power quality management. A load forecasting method based on load profile technology was proposed, and the node sensitivity analysis method was used to accurately and quickly solve the multi-objective optimization model, which realized the management of power quality problems such as node voltage overrun, voltage fluctuation, and harmonic pollution. Moreover, the method satisfied the high demand for penetration of photovoltaic and wind power into the distribution network. The main innovations of the work are: (a) comprehensively considering the impact of demand-side response and power quality issues on distribution network planning, and establishing an optimization model that includes power quality indicators and economic indicators; (b) based on the node sensitivity analysis method, so that it can be applied to complex power distribution systems; and (c) the heat load portrait technology is defined to provide a data basis for subsequent optimization. It can provide references for the promotion and application of distributed photovoltaic and wind power, and provide solutions for power quality management methods caused by high-power density distributed energy grid-connection. The next work will consider multiple types of energy storage equipment to further optimize the economy and safety.

Author Contributions

X.-R.L. conceived the idea for the manuscript and wrote the manuscript with input from F.-J.Z., Q.-Y.S., and P.J. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by National Key R&D Program of China under grant (2018YFA0702200), the Fundamental Research Funds for the Central Universities (N2004013).

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

ETSHEElectric Thermal Storage Heating Equipment
DEAData Envelopment Analysis
FCMFuzzy C-means
JValue function of FCM
mMembership factor
KNumber of clustering
UMembership matrix
CCluster center
W w i n d External wind speed value [m/s]
Δ T w i n d Equivalent temperature of wind speed cooling [°C]
α , β Wind speed correction coefficients
S s o l a r External light value [ W / m 2 ]
S k Light conversion coefficient [ W / ( m 2   ° C ) ]
f T H D I a Total current harmonic distortion
I 1 Root mean square value of the fundamental current [A]
I h Root mean square value of the h -th harmonic current [A]
U i j Voltage of node i at time j [V]
U N System standard voltage [V]
Δ I Current change caused by the output power change [A]
Z Equivalent impedance of the two-port network distribution network [Ω]
ρ S Unit price corresponding to the unit heat supply of the electric boiler [¥/kWh]
Q S e ( Q S ) Rated heat supply of the electric boiler [kWh]
ρ X Unit price of the heat regenerator [¥]
V X Volume of the heat accumulator [m3]
ρ B The increased equipment investment for each additional unit of the rated flow [¥]
γ Safety factor when calculating the rated flow
C 1 Conversion constant of kilowatt hours and joules
C W Specific heat of the heat medium [J/(kg · °C)]
ρ W Density of the heat medium [kg/m3]
Δ τ T Difference in heat storage temperature [°C]
ρ C The investment cost of the control device for each additional 1kW of the rated heating capacity of the electric boiler [¥]
k i n s Equipment installation cost coefficient
ϵ , ϑ , τ , ς Normalized weight coefficient
P i Active power of node i [kW]
Q i Reactive power of node i [KVar]
U i Voltage amplitudes of node i [V]
U j Voltage amplitudes of node j [V]
G i j Conductance of branch i j [S]
B i j Susceptance of branch i j [S]
θ i j Voltage phase angle difference between i -th and j -th node
Q q ( t ) Heat provided by the heat storage body to the user during time t [kWh]
Q l , m a x Maximum design daily heat load [kW]
r l ( t ) The ratio of heat load to maximum value in period t of design day
Q N ( t ) Heat stored in the heat storage body at the beginning of time t [kWh]
α s e n Sensitivity coefficient
w j Relative production efficiency
u r Weight of the r -th power quality indicator
v i Weight of the i -th quality level
x r j r -th power quality input indicator in the j -th DMU
y i j i -th quality level output in the j -th DMU indicator

Appendix A

  • Comprehensive evaluation method of power quality based on data envelope analysis
DEA is a comprehensive evaluation method for evaluating the relative efficiency of the same type of multi-index inputs and multi-index outputs [35]. In this method, the weights of the input and output variables are determined from the actual data of each decision unit, and the evaluation is performed from the perspective of the most favorable for the decision unit, which can avoid the influence of subjective factors. Additionally, it is not necessary to determine the display expressions of input and output, which simplifies the internalized relationship between input indicators and output indicators. In the case where the relationship between input indicators and output indicators is complex and difficult to determine, it can reflect its obvious advantages. The data envelopment analysis method was used to evaluate the power quality of each access point, and the access plan was evaluated for power quality at the planning stage, and finally, an ETSHE access plan with a good power quality level was obtained. The basic model and simplified transformation of DEA are as follows:
(1) C 2 R Model
The C 2 R model is the most basic DEA model. Suppose there are n evaluation objects, each Decision Making Units (DMU) has m types of input indicators and q types of output indicators. Construct the following C 2 R model [35]:
w j = i = 1 q v i y i j r = 1 m u r x r j
{ i = 1 q v i y i j / r = 1 m u r x r j     1 u r     0 i = 1 , 2 , , n ; v i     0    r = 1 , 2 , , m
where w j is the relative production efficiency; u r is the weight of the r -th power quality indicator ( r = 1,2, …, m ); v i indicates the weight of the i -th quality level ( i   = 1,2, …, q ); x r j indicates the r -th power quality input indicator in the j -th DMU; and y i j indicates the i -th quality level output in the j -th DMU indicator ( j = 1,2, …, n ). Choose a reasonable input weight vector v i and output weight vector u r so that w j     1 . When w j = 1, the j -th DMU is considered to be the most effective.
(2) Charnes–Cooper transformation
Due to the difficulty of solving the C 2 R model, the fractional programming is transformed into an equivalent linear programming problem using the Charnes–Cooper transformation, and the dual model is obtained using dual theory.
Carry out the Charnes–Cooper transformation on Equation (A2), and obtain a model in the form of equivalent linear programming as in Equation (A3):
{ min [ θ ε ( e ^ T s + e ^ T s + ) ] j = 1 n λ j x j + s = θ x 0 ;    j = 1 n λ j y j + s + = y 0 λ j 0 , j = 1 , 2 , , m ; s 0 , s + 0
where x j and y j represent the “input” and “output” of the j -th DMU, respectively. In order to solve in a larger feasible region, s = ( s 1 , s 2 , s m ) T and s + = ( s 1 + , s 2 + , s q + ) T are the relaxation variables. λ j represents the weight of D M U j 0 . θ is the efficiency value of the j -th DMU. Generally take ε = 10−7 in practice. e ^ T = (1,1,…,1) R s is used to adjust the relaxation variable.
  • Power quality evaluation of the ETSHE planning scheme
The power quality evaluation is carried out for the planning scheme of the ETSHE to be connected to the distribution network. First, the classification of grades should be carried out. There are many ways of classification, among which domestic literatures are mainly divided into grades 3, 5, 7, 9, and 10. This paper divided the power quality indicators into five grades. The assessment level was Q , which decreased from Q 1– Q 5 in order: high quality, good, average, qualified, and unqualified.
Considering the influence of the ETSHE connected to the distribution network, the evaluation indicators selected here were the total current harmonic distortion rate, voltage deviation, and voltage fluctuation. X 1, X 2, and X 3 represent the evaluation indicators, respectively. According to the GB/T14549 standard, the total current distortion rate of ETSHE, after it is connected to the distribution network, should be less than or equal to 2.2% [21,32]. According to the standard of GB/T12326-2008, the voltage fluctuation of the 10 kV distribution network should be less than or equal to 4% [23]. According to GB/T12325, the absolute value of the voltage deviation is less than or equal to 7% [33]. According to the national standard, the upper limit is equally divided into grade intervals, and the specific grade standards are shown in Table A1.
Table A1. Classification of power quality indicators.
Table A1. Classification of power quality indicators.
Evaluation LevelEvaluation Indicators
X1X2X3
Q1[0,0.55%][0,1.75%][0,1%]
Q2[0.55%,1.1%][1.75%,3.5%][1%,2%]
Q3[1.1%,1.65%][3.5%,5.25%][2%,3%]
Q4[1.65%,2.2%][5.25%,7%][3%,4%]
Q5[2.2%,30%][7%,30%][4%,30%]
In this paper, the data envelopment analysis method was used to evaluate the power quality, with X 1, X 2, and X 3 as the inputs and the evaluation level Q as the output to establish a DMU. The DMU format was ( X 1, X 2, X 3, Q ). The median and evaluation level of each level interval in Table 6 were taken as a single decision unit. The data of the standard decision unit of the classification are shown in Table A2. The power quality index grade standard in Table 6 and the DMU of the hierarchical standard in Table A2 were used as the total decision unit to perform the data envelope validity analysis.
Table A2. Hierarchical decision-making units.
Table A2. Hierarchical decision-making units.
DMUInput IndicatorOutput Indicators
X1X2X3Q
DMU10.275%0.875%0.5%1
DMU20.825%2.625%1.5%2
DMU31.375%4.375%2.5%3
DMU41.825%6.125%3.5%4
DMU516.1%18.5%4.5%5
Calculate the quality indicators of the access node according to the formula and the average of each indicator. According to the calculation formulas of the total current harmonic distortion rate, voltage deviation, voltage fluctuation, the quality indicators of the power quality of the access node on that day were calculated, respectively, and the average values of the indicators were calculated, as shown in Table A3.
Table A3. Data of various indicators of power quality.
Table A3. Data of various indicators of power quality.
Access PointPower Quality Index
X1X2X3
43.340%2.075%0.118%
84.008%4.753%0.198%
134.359%5.191%0.215%
255.200%2.054%0.208%
315.391%2.258%0.247%
As the efficiency value of the power quality indicators of each access point relative to each level was different from the decision unit, it was recorded as δ Q , and the level with the smallest absolute value of δ Q was used as the evaluation level. Table A4 shows the power quality evaluation results of each access point calculated by the data envelope analysis method.
Table A4. Proximity of comprehensive efficiency of each decision unit and standard unit.
Table A4. Proximity of comprehensive efficiency of each decision unit and standard unit.
Access Point δ Q 1 δ Q 2 δ Q 3 δ Q 4 δ Q 5 Evaluation Level
4−0.581−0.248−0.181−0.1520.1834
8−0.462−0.129−0.062−0.0330.3024
13−0.2580.0750.1420.1710.5062
25−0.050.2830.350.3790.7141
3100.3330.40.4290.7941
As can be seen from the result of the power quality evaluation of each access point by the data envelope analysis method, only the power quality evaluation level of access points 4 and 8 was a level 4 deviation. The level of access point 13 was level 2, and the power quality evaluation level of the access points at 25 and 31 was a quality of level 1. The evaluation level without access points was a level 5 failure, indicating that the access scheme is feasible and that all indicators of the access point can meet the requirements of the standard.

References

  1. Bracco, S.; Delfino, F.; Pampararo, F.; Robba, M.; Rossi, M. A dynamic optimization-based architecture for polygeneration microgrids with trigeneration, renewables, storage systems and electrical vehicles. Energy Convers. Manag. 2015, 96, 511–520. [Google Scholar] [CrossRef]
  2. Chathurangi, D.; Jayatunga, U.; Rathnayake, M. Potential power quality impacts on LV distribution networks with high penetration levels of solar PV. In Proceedings of the International Conference on Harmonics and Quality of Power (ICHQP), Ljubljana, Slovenia, 13–16 May 2018; Volume 8, pp. 1–6. [Google Scholar]
  3. Chai, Y.; Li, G.; Wang, C. Network partition and voltage coordination control for distribution networks with high penetration of distributed PV units. IEEE Trans. Power Syst. 2018, 33, 3396–3407. [Google Scholar] [CrossRef]
  4. Jin, M.; Feng, W.; Liu, P.; Marnay, C.; Spanos, C. MOD-DR: Microgrid optimal dispatch with demand response. Appl. Energy 2017, 187, 758–776. [Google Scholar] [CrossRef] [Green Version]
  5. Mokgonyana, L.; Zhang, J.; Li, H. Optimal location and capacity planning for distributed generation with independent power production and self-generation. Appl Energy 2017, 188, 140–150. [Google Scholar] [CrossRef] [Green Version]
  6. Xiao, J.; Zhang, Z.; Liang, H. Optimization method for location and capacity of public energy storage in distribution network. Autom. Electr. Power Syst. 2015, 19, 60–66 + 73. [Google Scholar]
  7. Cui, Y.; Liu, W.; Tang, Y. Optimal capacity allocation method for energy storage batteries to improve the voltage level of photovoltaic high-permeability distribution networks. J. Solar Energy 2017, 38, 1157–1165. [Google Scholar]
  8. Bracco, S. An optimization model for polygeneration microgrids with renewables, electrical and thermal storage: Application to the savona campus. In Proceedings of the 2018 IEEE International Conference on Environment and Electrical Engineering and 2018 IEEE Industrial and Commercial Power Systems Europe (EEEIC/I&CPS Europe), Palermo, Italy, 12–15 June 2018; Volume 5, pp. 1–6. [Google Scholar]
  9. Brignone, M.; Federico, D.; Francesco, D.; Fossa, M.; Pampararo, F. Modelling and simulating a thermal storage system for the Savona campus smart polygeneration micro grid. Model. Meas. Control C 2018, 79, 83–89. [Google Scholar] [CrossRef]
  10. Delfino, F. Identification and optimal control of an electrical storage system for microgrids with renewables. Sustain. Energy Grids Netw. 2019, 17, 100183. [Google Scholar] [CrossRef]
  11. Ghofrani, M.; Arabali, A.; Etezadi-Amoli, M.; Fadali, M.S. A Framework for Optimal Placement of Energy Storage Units within a Power System with High Wind Penetration. IEEE Trans. Sustain. Energy 2013, 4, 434–442. [Google Scholar] [CrossRef]
  12. Xie, X.; Wang, K.; Liu, Q. A method for location and capacity planning of distributed power sources considering harmonics. J. Wuhan Univ. 2019, 52, 1671–8844. [Google Scholar]
  13. Chen, C.; Jin, Y.; Cai, Z. Node sensitivity analysis method for power quality management of high-density distributed photovoltaic access distribution network. Power Capacit. React. Compens. 2019, 3, 137–142. [Google Scholar]
  14. Rahman, S.; Hazim, O. A generalized knowledge-based short-term load-forecasting technique. IEEE Trans. Power Syst. 1993, 8, 508–514. [Google Scholar] [CrossRef]
  15. Krinidis, S.; Chatzis, V. A robust fuzzy local information c-means clustering algorithm. IEEE Trans. Image Process. 2010, 19, 1328–1337. [Google Scholar] [CrossRef] [PubMed]
  16. Qin, L.; Zhao, N.; Sun, D.F. Analysis of measured results of heating energy consumption of residential buildings in severe cold areas. J. Northeast Dianli Univ. 2016, 36, 36–40. [Google Scholar]
  17. Cai, Q. Study on the relationship between meteorological factors and heating load. Dist. Heat. 2016, 4, 27–32. [Google Scholar]
  18. Seyed, M.S. Electric load forecasting using an artificial neural network. IEEE Trans. Power Syst. 2013, 6, 442–449. [Google Scholar]
  19. Amjady, N. Short-term hourly load forecasting using time-series modeling with peak load estimation capability. IEEE Trans. Power Syst. 2016, 4, 798–805. [Google Scholar]
  20. Xu, X.; Zhang, Q.; Mou, Y. Server load prediction algorithm based on CM-MC for cloud systems. J. Syst. Eng. Electron. 2018, 29, 1069–1078. [Google Scholar]
  21. National Standard of the People’s Republic of China. GBT 24337-2009, Power Harmonics among Public Power Grids; Jointly Issued by General Administration of Quality Supervision, Inspection and Quarantine (GAQSIQ) and Standardization Administration (SAC) of the People’s Republic of China: Beijing, China, 2009.
  22. Wang, Y.; Wen, W.; Wang, C. Adaptive Voltage Droop Method of Multiterminal VSC-HVDC Systems for DC Voltage Deviation and Power Sharing. IEEE Trans. Power Deliv. 2019, 34, 169–176. [Google Scholar]
  23. National Standard of the People’s Republic of China. GBT 12326-2008, Voltage Fluctuations and Flickers; Jointly Issued by General Administration of Quality Supervision, Inspection and Quarantine (GAQSIQ) and Standardization Administration (SAC) of the People’s Republic of China: Beijing, China, 2008.
  24. Cui, Q.; Bai, X.; Dong, W. Collaborative planning of distributed wind power generation and distribution network with large-scale heat pumps. CSEE J. Power Energy Syst. 2019, 5, 335–347. [Google Scholar]
  25. Barrett, E.; Eustis, C.; Bass, R.B. A Dual-Heat-Pump Residential Heating System for Shaping Electric Utility Load. IEEE Power Energy Technol. Syst. J. 2018, 5, 56–64. [Google Scholar] [CrossRef]
  26. Zhang, W. Application examples of water heat storage technology for electric boilers. HVAC 2008, 94–97. [Google Scholar]
  27. Dong, X. Power Flow Analysis Considering Automatic Generation Control for Multi-Area Interconnection Power Networks. IEEE Trans. Ind. Appl. 2017, 53, 5200–5208. [Google Scholar] [CrossRef]
  28. Christakou, K.; LeBoudec, J.; Paolone, M.; Tomozei, D. Efficient Computation of Sensitivity Coefficients of Node Voltages and Line Currents in Unbalanced Radial Electrical Distribution Networks. IEEE Trans. Smart Grid 2013, 4, 741–750. [Google Scholar] [CrossRef] [Green Version]
  29. National Standard of the People’s Republic of China. GBT 15543-2008, Power Quality Three-Phase Voltage Unbalance; Jointly Issued by General Administration of Quality Supervision, Inspection and Quarantine (GAQSIQ) and Standardization Administration (SAC) of the People’s Republic of China: Beijing, China, 2008.
  30. Sato, M.; Fukuyama, Y.; Iizaka, T.; Matsui, T. Total Optimization of Energy Networks in a Smart City by Multi-Swarm Differential Evolutionary Particle Swarm Optimization. IEEE Trans. Sustain. Energy 2019, 10, 2186–2200. [Google Scholar] [CrossRef]
  31. Duan, H.; Li, P.; Yu, Y. A predator-prey particle swarm optimization approach to multiple UCAV air combat modeled by dynamic game theory. IEEE/CAA J. Autom. Sin. 2015, 2, 11–18. [Google Scholar]
  32. National Standard of the People’s Republic of China. GBT 29319-2012, Technical Regulations for Access to Distribution Networks of Photovoltaic Power Generation Systems; Jointly Issued by General Administration of Quality Supervision, Inspection and Quarantine (GAQSIQ) and Standardization Administration (SAC) of the People’s Republic of China: Beijing, China, 2012.
  33. National Standard of the People’s Republic of China. GBT 12325-2008, Power Quality Voltage Deviation; Jointly Issued by General Administration of Quality Supervision, Inspection and Quarantine (GAQSIQ) and Standardization Administration (SAC) of the People’s Republic of China: Beijing, China, 2008.
  34. Nourizadeh, S.; Karimi, M.J.; Ranjbar, A.M. Power system stability assessment during restoration based on a wide area measurement system. Gener. Trans. Distrib. 2012, 6, 1171–1179. [Google Scholar] [CrossRef]
  35. Matin, R.K.; Ghahfarokhi, M.I. A two-phase modified slack-based measure approach for efficiency measurement and target setting in data envelopment analysis with negative data. IMA J. Manag. Math. 2015, 26, 83–98. [Google Scholar] [CrossRef]
Figure 1. Optimization configuration model.
Figure 1. Optimization configuration model.
Applsci 10 02666 g001
Figure 2. Schematic of the sensitivity radar.
Figure 2. Schematic of the sensitivity radar.
Applsci 10 02666 g002
Figure 3. Improved particle swarm optimization algorithm flowchart.
Figure 3. Improved particle swarm optimization algorithm flowchart.
Applsci 10 02666 g003
Figure 4. Fuzzy c-means clustering results of four types of users.
Figure 4. Fuzzy c-means clustering results of four types of users.
Applsci 10 02666 g004
Figure 5. Energy consumption characteristics of the four typical users.
Figure 5. Energy consumption characteristics of the four typical users.
Applsci 10 02666 g005
Figure 6. IEEE 33-node power distribution system wiring diagram.
Figure 6. IEEE 33-node power distribution system wiring diagram.
Applsci 10 02666 g006
Figure 7. Current harmonic distortion rate curve. (a) Comparison of sub-scenarios of nodes 4 and 8. (b) Comparison of sub-scenarios of nodes 25 and 31. (c) Percentage improvement of current harmonic distortion rate relative to the standard value.
Figure 7. Current harmonic distortion rate curve. (a) Comparison of sub-scenarios of nodes 4 and 8. (b) Comparison of sub-scenarios of nodes 25 and 31. (c) Percentage improvement of current harmonic distortion rate relative to the standard value.
Applsci 10 02666 g007
Figure 8. Voltage deviation curve. (a) Comparison of sub-scenarios of nodes 4 and 8. (b) Comparison of sub-scenarios of nodes 25 and 31. (c) Percentage improvement of voltage deviation relative to standard value.
Figure 8. Voltage deviation curve. (a) Comparison of sub-scenarios of nodes 4 and 8. (b) Comparison of sub-scenarios of nodes 25 and 31. (c) Percentage improvement of voltage deviation relative to standard value.
Applsci 10 02666 g008
Figure 9. Voltage fluctuation curve. (a) Comparison of sub-scenarios of nodes 4 and 8. (b) Comparison of sub-scenarios of nodes 25 and 31. (c) Percentage improvement of voltage fluctuation relative to standard value.
Figure 9. Voltage fluctuation curve. (a) Comparison of sub-scenarios of nodes 4 and 8. (b) Comparison of sub-scenarios of nodes 25 and 31. (c) Percentage improvement of voltage fluctuation relative to standard value.
Applsci 10 02666 g009aApplsci 10 02666 g009b
Figure 10. Current harmonic distortion rate curve.
Figure 10. Current harmonic distortion rate curve.
Applsci 10 02666 g010
Figure 11. Voltage deviation curve.
Figure 11. Voltage deviation curve.
Applsci 10 02666 g011
Figure 12. Voltage fluctuation curve.
Figure 12. Voltage fluctuation curve.
Applsci 10 02666 g012
Figure 13. Percentage improvement of current harmonic distortion rate relative to the standard value.
Figure 13. Percentage improvement of current harmonic distortion rate relative to the standard value.
Applsci 10 02666 g013
Figure 14. Percentage improvement of voltage deviation relative to the standard value.
Figure 14. Percentage improvement of voltage deviation relative to the standard value.
Applsci 10 02666 g014
Figure 15. Percentage improvement of voltage fluctuation relative to the standard value.
Figure 15. Percentage improvement of voltage fluctuation relative to the standard value.
Applsci 10 02666 g015
Table 1. Users label description.
Table 1. Users label description.
No.LabelsDescription
1Type of heatingAccording to the characteristics of load timing sequence, it is divided into residents, education, commerce, and industry.
2Seasonal characteristicsAccording to users sensitivity to different seasons, it can be divided into summer sensitivity, winter sensitivity, summer–winter sensitivity and non-sensitivity.
3Temperature sensitivity:According to the load changing trend of users at different temperatures, it is divided into low temperature sensitivity, high temperature sensitivity, high–low temperature sensitivity and non-sensitivity.
4Growth rate of heatingAccording to the heating growth of users, it can be divided into rising, falling, and flat.
5Peak and valley characteristicsAccording to the users sensitivity to the peak–valley electricity price system, it can be divided into peak–valley equilibrium, valley preference, and peak preference.
6Weekly rest characteristicsAccording to the characteristics of the users week, can be divided into five days, six days, seven days.
Table 2. Four types of customer label characteristics.
Table 2. Four types of customer label characteristics.
SeasonTemperature SensitivityHeating Growth RatePeak and Valley PreferencesWeekly Rest Characteristics
ResidentsSensitive in summer and winterLow temperature sensitivityRisePeak-valley equilibriumSeven-day work
EducationWinter sensitiveLow temperature sensitivityFlatValley preferenceSeven-day work
BusinessWinter sensitiveHigh and low temperature sensitiveFlatPeak-valley equilibriumSix-day work
IndustrySummer and winter sensitiveHigh and low temperature sensitiveRiseValley preferenceSeven-day work
Table 3. Load, wind power and photovoltaic power of a typical day (MW).
Table 3. Load, wind power and photovoltaic power of a typical day (MW).
PeriodTotal LoadWind PowerPhotovoltaic PowerPeriodTotal LoadWind PowerPhotovoltaic Power
14.80.420.00135.30.380.50
24.80.400.00145.40.300.44
34.80.480.00155.40.340.45
44.80.490.00165.40.380.28
54.80.320.00175.60.330.24
64.80.280.04186.20.330.02
75.30.200.21196.30.130.00
85.30.360.31206.70.360.00
95.50.420.41215.80.450.00
105.60.380.35225.30.410.00
115.60.330.43235.30.430.00
125.50.340.45244.80.400.00
Table 4. Optimized Configuration Results.
Table 4. Optimized Configuration Results.
Node Number P D G / M W
4 (Resident load)0.3
8 (Business load)1.0
13 (Resident load)0.3
25 (Industrial load)1.1
31 (Education load)0.16
Table 5. Comparison analysis before and after node sensitivity analysis.
Table 5. Comparison analysis before and after node sensitivity analysis.
Optimal Planning Model of ETSHEBefore AnalysisAfter AnalysisEffectiveness Enhancement
Calculation time/s153.6121.560%
Table 6. Optimized configuration results.
Table 6. Optimized configuration results.
Node Number P D G / M W
3 (Resident load)0.32
6 (Resident load)0.35
10 (Resident load)0.33
35 (Business load)1.2
38 (Business load)1.3
44 (Industrial load)1.5
47 (Industrial load)1.6
48 (Industrial load)2.1
61 (Resident load)0.41
66 (Resident load)0.35
77 (Resident load)0.43
78 (Business load)1.1
79 (Business load)0.9
91 (Resident load)0.45
92 (Resident load)0.32
93 (Resident load)0.5
114 (Education load)0.41
115 (Education load)0.35
117 (Education load)0.28
Table 7. Comparison analysis before and after node sensitivity analysis.
Table 7. Comparison analysis before and after node sensitivity analysis.
Optimal Planning Model of ETSHEBefore AnalysisAfter AnalysisEffectiveness Enhancement
Calculation time/s267.3165.973.7%

Share and Cite

MDPI and ACS Style

Liu, X.-R.; Zhang, F.-J.; Sun, Q.-Y.; Jin, P. Location and Capacity Selection Method for Electric Thermal Storage Heating Equipment Connected to Distribution Network Considering Load Characteristics and Power Quality Management. Appl. Sci. 2020, 10, 2666. https://doi.org/10.3390/app10082666

AMA Style

Liu X-R, Zhang F-J, Sun Q-Y, Jin P. Location and Capacity Selection Method for Electric Thermal Storage Heating Equipment Connected to Distribution Network Considering Load Characteristics and Power Quality Management. Applied Sciences. 2020; 10(8):2666. https://doi.org/10.3390/app10082666

Chicago/Turabian Style

Liu, Xin-Rui, Fu-Jia Zhang, Qiu-Ye Sun, and Peng Jin. 2020. "Location and Capacity Selection Method for Electric Thermal Storage Heating Equipment Connected to Distribution Network Considering Load Characteristics and Power Quality Management" Applied Sciences 10, no. 8: 2666. https://doi.org/10.3390/app10082666

APA Style

Liu, X. -R., Zhang, F. -J., Sun, Q. -Y., & Jin, P. (2020). Location and Capacity Selection Method for Electric Thermal Storage Heating Equipment Connected to Distribution Network Considering Load Characteristics and Power Quality Management. Applied Sciences, 10(8), 2666. https://doi.org/10.3390/app10082666

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

Article Metrics

Back to TopTop