Abstract
With the extensive integration of high-penetration renewable energy resources, more fast-response frequency regulation (FR) providers are required to eliminate the impact of uncertainties from loads and distributed generators (DGs) on system security and stability. As a high-quality FR resource, community integrated energy station (CIES) can effectively respond to frequency deviation caused by renewable energy generation, helping to solve the frequency problem of power system. This paper proposes an optimal planning model of CIES considering FR service. First, the model of FR service is established to unify the time scale of FR service and economic operation. Then, an optimal planning model of CIES considering FR service is proposed, with which the revenue of participating in the FR service is obtained under market mechanism. The flexible electricity pricing model is introduced to flatten the peak tie-line power of CIES. Case studies are conducted to analyze the annual cost and the revenue of CIES participating in FR service, which suggest that providing ancillary services can bring potential revenue.
WITH the increasing penetration of distributed generators (DGs) at the user side, large-scale intermittent and random power generation of DGs aggravate frequency deviation and bring stability problem to power system [
The application of combined heat and power (CHP) technology brings heating, natural gas, and electric power systems together [
A number of studies have investigated CIES planning. In [
With the development of the ancillary service market, CIES can participate in the FR market competition as a fast-response load. Previous works have been proposed to consider FR demand in the planning stage. In [
Thus, an optimal planning model of CIES considering FR service is proposed in this paper. First, the mechanism and simplified transaction procedure of the FR market are presented. An optimal planning model of CIES considering FR service is proposed with a minimal annual cost. Conversion devices, energy storage and distributed generators are considered in the model. In addition, flexible electricity pricing model is introduced to mitigate the fluctuation of tie-line power. Taking a real industrial park as the example, the optimal configuration of equipment and annual comprehensive cost are compared.
The main contributions are summarized as follows.
1) An optimal planning model of CIES considering FR is established. First, the model of FR service is established to unify the time scale of FR service and economic operation of CIES. Then, considering FR service and operation constraints, the planning model of CIES is proposed, with the objective function of minimizing annual cost. CIES can provide FR service, while efficiently satisfying its internal cooling, heating and electric demands. Both equipment capacity and FR bidding strategy can be optimized. Case study on real load data is conducted to analyze the benefit of the CIES participating in the FR market.
2) Given that the participation of CIES in the FR market and concentrated power purchasing may lead to power fluctuation in distribution network, a flexible electricity pricing model is introduced to flatten the peak power demand of the CIES. The simulation results indicate that the pricing model can effectively mitigate the tie-line power fluctuation and over-limit problem.
The remainder of this paper is organized as follows. In Section II, the mechanism and simplified transaction procedure of the FR market are introduced. In Section III, a planning model of CIES participating in the FR market is established. In Section IV, taking a practical demand of an industrial park as the example, numerical results and comparisons demonstrate the effectiveness of the proposed method. Finally, conclusions are drawn in Section V.
FR service is an important part of the ancillary service market. The main function of FR service is to track the load change, maintain the real-time balance of electric power, and keep the system frequency stable [
In the bidding process, each FR service provider reports the information such as available capacity, response speed, callable time and price [

Fig. 1 PJM FR signal.
The bidding interval is 5 min in PJM market [
(1) |
(2) |
where is the FR power at time ; is the value of the FR signal at time ; is the FR capacity at period ; and is the upper limit of the FR capacity.
The other variable to be optimized in FR service is the corresponding electric power adjustment value. To quantify the FR power adjustment, the RegUp efficiency and RegDown efficiency are defined, which represent the fractions of the assigned capacity actually deployed for RegUp and RegDown, respectively, and are calculated as follows:
(3) |
(4) |
where and are the RegUp and RegDown efficiencies at period , respectively; and and are the RegUp and RegDown signals, respectively.
In the planning problem, we concentrate on the power consumption related to FR service in each hour, without consideration on the dynamic FR power change in time scale of several seconds. The electric power adjustment at period corresponding to the FR service can be calculated as follows:
(5) |
In PJM ancillary service market, FR revenue is settled according to the actual contributions of each participant [
(6) |
(7) |
(8) |
where and are the capacity clearing price and the performance clearing price at period , respectively; is the performance score [
RegD is a fast and dynamic signal that requires the resources to respond nearly instantaneously, such as electric energy storage equipment. RegA is a slower signal for the resources such as a gas turbine. It is used to compensate larger and longer fluctuations in system conditions. For a RegD system, the mileage ratio is defined as [
(9) |
where and are the mileages of RegD and RegA signals, respectively.
Mileage is defined as the accumulated movement requested by the regulation control signal in a period. For example, the RegD mileage over a one-hour period is defined as [
(10) |
where is the value of RegD signal at time .
Similarly, the RegA mileage is calculated as follows:
(11) |
where is the value of RegA signal at time .
In this section, a planning model of a CIES is established considering its participation in the FR market.
The structure of CIES is shown in

Fig. 2 Structure of CIES with candidate devices.
The loads consist of electric load, cooling load and heating load (includes hot water load and steam load), which cover the general needs of users. For a CIES to be constructed, the main energy resources are electricity, gas and renewable energy. Typical electric energy conversion equipments such as the electric boiler (EB) and electric chiller (EC) are considered to supply the cooling and heating loads. CHP is chosen as the coupling equipment between power and natural gas. Ground source heat pump (HP) transfers the heat from the ground for space heating or cooling, which can effectively reduce annual maintenance costs and carbon emissions. Moreover, DGs such as PV and wind turbine (WT), electrical storage (ES), and heat storage (HS) are considered to achieve better economy and sustainability [
Time-of-use electricity price is widely adopted in distribution system. It provides an incentive to charge the CIES during off-peak periods in order to achieve the lowest operation cost [
A real-time electricity pricing model is introduced to mitigate the power fluctuation and over-limit at the tie-line [
(12) |
where is the price of electricity affected by both fixed and flexible loads at period ; is the basic electricity price at period ; is the price sensitivity coefficient at period ; and is the flexible load at period . This approximation allows for using quadratic programming to solve the CIES planning problem, while it includes the feedback of purchased power on the electricity prices. By considering this feedback mechanism, the power demand of CIES will not be supplied in a short time interval at the lowest prices. This paper introduces the flexible electricity pricing model to mitigate the power fluctuation at tie-line. CIES can regulate its power purchase and respond to the distribution network based on market pricing mechanism. When the purchased power does not exceed the limit, the electricity purchase price is equal to the system basic electricity price. When the purchased power exceeds the limit, the electricity purchase price will be exactly relevant to the power demand, which is higher than the system basic electricity price. For the CIES, the electricity price and purchased power can be expressed as follows:
(13) |
where is the purchased power at period ; and is the upper limit of purchased power.
The minimum annual cost is taken as the objective function, which consists of investment cost , maintenance cost , operation cost , and FR revenue .
(14) |
(15) |
where is the discount rate; is the service life; is the set of candidate devices; is the investment cost of per unit capacity of facility ; is the minimum configuration unit of facility; and is the amount of which is an integer variable.
(16) |
where is the maintenance cost of per-unit power of facility; is the set of typical day data; and is the output power of facility at period of day in month .
The operation cost consists of the electricity purchase cost and gas purchase cost . The electricity purchase cost related to FR service is included in (18).
(17) |
(18) |
(19) |
where is the purchased power of the CIES at period of day in month ; and are the RegUp and RegUp efficiencies, respectively; is the FR capacity; is the electricity purchase cost of FR service at period ; is the gas price at period t; is the input gas volume of CHP; and is the calorific value of natural gas.
1) Operation Constraints of EB
(20) |
(21) |
where and are the electric and heating power of the EB at period , respectively; is the efficiency of the EB; is the minimum planning unit of the EB; and is the amount of .
2) Operation Constraints of EC
(22) |
(23) |
where and are the electric and cooling power of the EC at period , respectively; is the coefficient of performance (COP) of the EC; is the minimum planning unit of the EC; and is the amount of .
3) Operation Constraints of CHP
In this planning model, we consider CHP supplies heating and electric power at a fixed ratio.
(24) |
(25) |
(26) |
where and are the heating and electric power of CHP at period , respectively; is the input gas volume of CHP at period ; and are the efficiencies of heating and electricity of CHP, respectively; is the minimum planning unit of CHP; and is the amount of .
4) Operation Constraints of HP
Unlike CHP, the HP can only produce either heating or cooling power in a given time interval. This paper assumes that the HP supplies cooling power in summer and heating power in winter.
(27) |
(28) |
(29) |
where and are the heating and cooling power of the HP at period , respectively; is the electric power of the HP at period ; and are the COPs of heating and cooling of HP, respectively; is the minimum planning unit of HP; and is the amount of .
5) Operation Constraints of HS
(30) |
(31) |
(32) |
where is the heat energy stored in HS at period ; T is the number of periods in one planning cycle; is the heat loss coefficient of HS; and are the charging and discharging power of HS at period , respectively; and are the charging and discharging efficiencies of HS, respectively; is the time interval; is the minimum planning unit of HS; and is the amount of .
6) Operation Constraints of ES
To ensure the long life-time of ES, the upper and lower limits of state of charge (SOC) are used to constrain the ES from being fully charged or discharged [
(33) |
(34) |
(35) |
(36) |
(37) |
(38) |
(39) |
where is the electricity stored in ES at period ; is the heat loss coefficient of ES; and are the charging and discharging power at period , respectively; and are the charging and discharging efficiencies of ES, respectively; is the minimum planning unit of ES; is the amount of ; is the SOC of ES at period ; and are the lower and upper limit of SOC, respectively; and are the upper limits of charging and discharging power, respectively; and and are Boolean variables of charging and discharging modes, respectively.
7) Operation Constraints of Converter
When the ES is charged or discharged, the energy conversion is realized through the converter. The capacity of the converter should be larger than the charging/discharging power and the FR capacity.
(40) |
(41) |
where is the minimum planning unit of converter; and is the amount of .
8) Operation Constraints of PV
The electric power of PV depends on many factors such as solar irradiance and temperature. However, to reduce the complexity, the electric power of PV is regarded to be directly determined by the capacity and the solar irradiance [
(42) |
where is the solar irradiance at period ; is the standard solar irradiance; is the minimum planning unit of PV; and is the amount of .
9) Operation Constraints of WT
The electric power of WT is directly determined by the capacity and the wind speed.
(43) |
where is the wind speed at period ; , and are the cut-in, rated and cut-off wind speeds, respectively; is the minimum planning unit of WT; and is the amount of .
10) Operation Constraints of FR
(44) |
(45) |
where is a decision variable. When , the purchased power exceeds the limit, and CIES cannot provide FR service. In the optimization model, (45) indicates whether the purchased power exceeds the limit.
11) Operation Constraints of Tie-line
For the planning of the CIES, it is generally considered that the electric, heating, and cooling loads remain constant in one-hour period. Thus, the tie-line power is also considered to be a constant value. However, while participating in the FR service, the CIES is required to provide power that can track the rapid change of the FR signal. Therefore, the tie-line power can be divided into two parts: the constant purchased power and the dynamic FR power.
(46) |
where is the tie-line power at period ; and is the FR power at time .
The FR power is not larger than the FR capacity , so the maximum power of the tie-line is as follows:
(47) |
(48) |
where is the upper limit of the tie-line power; and is the maximum power transfer limit of the tie-line.
12) Power Balance Constraints
1) Electric load balance
(49) |
where is the electric load at period .
2) Heating load balance
(50) |
where is the heating load at period .
3) Cooling load balance
(51) |
where is the cooling load at period .
The decision variables mainly include amounts of equipment planning unit, FR capacity, purchased power, gas volume purchased at period , output power of each equipment, SOC, and charging and discharging modes of ES. The planning model of the CIES can be described as follows:
(52) |
Formulae (1)-(5) and (44)-(45) are the FR constraints; (20)-(43) are the operational constraints of equipment; (46)-(48) are the tie-line power constraints; and (49)-(51) are the power balance constraints.
In this paper, a practical demand of an industrial park [

Fig. 3 Typical daily load in June.

Fig. 4 Typical daily load in November.
The proposed planning model belongs to mixed-integer quadratic programming (MIQP) model. The model is implemented in the OPTI optimization toolbox using MATLAB R2016a and solved by IBM ILOG CPLEX 12.6 [
To analyze the influence of FR service and the flexible electricity pricing model on the planning results, three scenarios are selected and shown as follows.
Scenario 1: the CIES does not participate in the FR market, and the flexible electricity pricing model and tie-line constraints are not considered.
Scenario 2: the CIES participates in the FR market, but the flexible electricity pricing model and tie-line constraints are not considered.
Scenario 3: the CIES participates in the FR market, and the flexible electricity pricing model and tie-line constraints are both considered.
The optimal configurations in the three scenarios are shown in
Comparing Scenarios 1 and 2, since the CIES in Scenario 2 provides FR service, the power exchanged between the ES and the electric bus increases, so the capacity of the converter increases from 4400 kW to 5900 kW. However, the electric power adjustment values of charging and discharging corresponding to the FR service are almost balanced in each time period, so nearly no additional ES capacity is needed.
It can be observed that ES capacity in Scenario 3 is reduced by 700 kWh compared with Scenario 2, which proves that the power purchase congestion problem has been mitigated due to the introduction of flexible electricity pricing model. Meanwhile, the capacity of converter is reduced by 1400 kW, indicating that the maximum value of charging/discharging decreases and the FR capacity of the CIES is reduced compared with Scenario 2. Because of the input electric power limitation of EB and HP, HS is required to supply more heating power, so its capacity increases by 2200 kWh. In addition, the capacities of PV and WT increase by 400 kW and 200 kW, respectively, indicating that CIES can accommodate higher DG penetration.
In
It is obvious that CIES participation in the FR market can effectively reduce the annual cost. Through CIES participation in the ancillary service market, additional revenue is obtained while other costs change very little.
The annual cost increases by 605300 CNY in Scenario 3 with an increase of 9.50% compared with the optimization result of Scenario 2. Due to the influence of the flexible electricity pricing model, additional capacities of DGs and HS are required, so the investment cost increases by 273000 CNY. Considering that less electricity is purchased when the loads are at the valley periods, the electricity consumption is reduced by 325300 kWh. In addition, as the tie-line constraints are considered, when the purchased power reaches the limit, the CIES cannot provide FR service, so the FR revenue decreases by 295500 CNY.
In Scenario 1, due to the time-of-use electricity price, the CIES concentrates on purchasing electricity when the prices are low through ES and stops purchasing electricity when the loads are at the peak periods.
In Scenario 2, because the CIES provides the FR service, the tie-line power is a combination of the purchased power and FR power. The power fluctuation and over-limit conditions are not alleviated. In contrast, the maximum tie-line power further increases. Thus, the flexible electricity pricing model is introduced in Scenario 3 to mitigate the tie-line power fluctuation and prevent the purchased power from exceeding the limit. A comparison of the optimal operation strategies between Scenarios 2 and 3 is presented in

Fig. 5 Annual purchased power in Scenarios 2 and 3. (a) Scenario 2. (b) Scenario 3.
The introduction of DGs results in reduced electricity purchase. The power of the DGs can cover the electric load when the demand is low. However, when the load demands are high or the electricity prices are low, the purchased power in Scenario 2 exceeds the limit, while in Scenario 3 CIES can adjust the power purchasing strategy and energy storage operation mode so that the purchased power reaches but does not exceed the limit. The comparison suggests that the flexible electricity pricing model adopted in this paper can effectively mitigate the tie-line power fluctuation and prevent the power from exceeding the limit.
The rates of the purchased power exceeding the limit in Scenarios 2 and 3 are shown in
The power purchasing strategies of the CIES on the typical day in June are compared in

Fig. 6 Comparison of purchased power on typical days in June and November. (a) June. (b) November.
In
Based on the previous power purchase plan in

Fig. 7 Comparison of electricity prices in Scenarios 2 and 3.
Since the CIES also provides FR service, the instantaneous tie-line power is equal to the sum of the purchased power and the real-time FR power. When the CIES provides the RegDown service (absorbing the excess power), the tie-line power may exceed the limit while the purchased power does not exceed the limit. These short-duration peaks will affect the secure operation of the distribution network.

Fig. 8 Maximum tie-line power in Scenarios 2 and 3 in June. (a) Scenario 2. (b) Scenario 3.
The maximum tie-line power of the typical day in November is compared in

Fig. 9 Maximum tie-line power in Scenarios 2 and 3 in November. (a) Scenario 2. (b) Scenario 3.
Due to the large heating and electricity loads of CIES, the purchased power is at the peak in periods of the 1
This paper proposes an optimal planning model of CIES considering FR service. After modeling the FR service, an optimal planning model of CIES considering FR service is established to optimize the type and capacity of candidate devices. With the minimum annual comprehensive cost as the objective function, it is considered that the CIES can participate in the FR market to obtain revenue. The flexible electricity pricing model is introduced to mitigate the power fluctuation problem of tie-line. Using practical demand of an industrial park, case study is conducted to verify the effectiveness of the proposed model. The result shows that the participation in FR service can provide fast-response regulation capacity as well as effectively reduce the annual cost of CIES. Additionally, the feedback mechanism of electricity price can also mitigate the power fluctuation of tie-line and prevent power over-limit. In the further work, several notable issues are worth studying. The siting of CIES can be considered in the planning model, which will affect the configuration of device and trading strategy. In FR market, it can be investigated that slow-response equipment in CIES participates in the traditional regulation market.
References
W. Yao, L. Jiang, J. Wen et al., “Wide-area damping controller for power system interarea oscillations: a networked predictive control approach,” IEEE Transactions on Control Systems Technology, vol. 23, no. 1, pp. 27-36, Jan. 2015. [百度学术]
H. T. Ma and B. H. Chowdhury, “Working towards frequency regulation with wind plants: combined control approaches,” IET Renewable Power Generation, vol. 4, no. 4, pp. 308-316, Aug. 2010. [百度学术]
Z. Tang, Y. Lim, S. Morris et al., “A comprehensive work package for energy storage systems as a means of frequency regulation with increased penetration of photovoltaic systems,” International Journal of Electrical Power & Energy Systems, vol. 110, pp. 197-207, Sept. 2019. [百度学术]
H. Bevrani, A. Ghosh, and G. Ledwich, “Renewable energy sources and frequency regulation: survey and new perspectives,” IET Renewable Power Generation, vol. 4, no. 5, pp. 438-457, Oct. 2010. [百度学术]
F. Brahman, M. Honarmand, and S. Jadid, “Optimal electrical and thermal energy management of a residential energy hub, integrating demand response and energy storage system,” Energy and Building, vol. 90, pp. 65-75, Mar. 2015. [百度学术]
T. Niknam, M. Bomapour, and A. Gheisari, “Combined heat, power and hydrogen production optimal planning of fuel cell power plants in distribution networks,” Energy Conversion and Management, vol. 66, pp. 11-25, Feb. 2013. [百度学术]
C. Wang, C. Lv, P. Li et al., “Modeling and optimal operation of community integrated energy systems: a case study from China,” Applied Energy, vol. 230, pp. 1242-1254, Nov. 2018. [百度学术]
X. Hu, C. Zou, C. Zhang et al., “Technological developments in batteries: a survey of principal roles, types, and management needs,” IEEE Power and Energy Magazine, vol. 15, no. 5, pp. 20-31, Aug. 2017. [百度学术]
H. Ren and W. Gao, “A MILP model for integrated plan and evaluation of distributed energy systems,” Applied Energy, vol. 87, no. 3, pp. 1001-1014, Mar. 2010. [百度学术]
X. Wu, X. Hu, X. Yin et al., “Optimal battery sizing of smart home via convex programming,” Energy, vol. 140, pp. 444-453, Dec. 2017. [百度学术]
W. Huang, N. Zhang, J. Yang et al., “Optimal configuration planning of multi-energy systems considering distributed renewable energy,” IEEE Transactions on Smart Grid, vol. 10, no. 2, pp. 1452-1464, Mar. 2019. [百度学术]
M. Sharafi and T. Y. ELMekkawy, “Multi-objective optimal design of hybrid renewable energy systems using PSO-simulation based approach,” Renewable Energy, vol. 68, pp. 67-79, Aug. 2014. [百度学术]
D. Zhang, S. Evangelisti, P. Lettieri et al., “Optimal design of CHP-based microgrids: multiobjective optimisation and life cycle assessment,” Energy, vol. 85, pp. 181-193, Jun. 2015. [百度学术]
P. Gaur, D. Bhowmik, and N. Soren, “Utilisation of plug-in electric vehicles for frequency regulation of multi-area thermal interconnected power system,” IET Energy Systems Integration, vol. 1, no. 2, pp. 88-96, Mar. 2019. [百度学术]
X. Pan, H. Xu, J. Song et al., “Capacity optimization of battery energy storage systems for frequency regulation,” in Proceedings of 2015 IEEE International Conference on Automation Science and Engineering (CASE), Gothenburg, Sweden, Oct. 2015, pp. 1139-1144. [百度学术]
X. Zhang, K. W. Chan, H. Wang et al., “Game-theoretic planning for integrated energy system with independent participants considering ancillary services of power-to-gas stations,” Energy, vol. 176, pp.249-264, Jun. 2019. [百度学术]
B. Vatandoust, A. Ahmadian, M. A. Golkar et al., “Risk-averse optimal bidding of electric vehicles and energy storage aggregator in day-ahead frequency regulation market,” IEEE Transactions on Power Systems, vol. 3, no. 3, pp. 2036-2047, May 2019. [百度学术]
P. Manjarres and O. Malik, “Frequency regulation by fuzzy and binary control in a hybrid islanded microgrid,” Journal of Modern Power Systems and Clean Energy, vol. 3, no. 3, pp. 429-439, Aug. 2015. [百度学术]
R. H. Byrne, R. J. Concepcion, and C. Silva-Monroy, “Estimating potential revenue from electrical energy storage in PJM,” in Proceedings of 2016 IEEE PES General Meeting (PESGM), Boston, USA, Jul. 2016, pp. 1-5. [百度学术]
H. Ma, B. Wang, W. Gao et al., “Optimal scheduling of a regional integrated energy system with energy storage systems for service regulation,” Energies, vol. 11, no. 1, pp. 195-213, Jan. 2018. [百度学术]
PJM. (2019. Aug.). Regulation requirement definition. [Online]. Available: https://www.pjm.com/-/media/markets-ops/ancillary/regulation-requirement-definition.ashx?la=en [百度学术]
J. Shen, C. Jiang, Y. Liu et al., “Microgrid operation optimization with regulation of grid tie-line power fluctuations and risk management,” International Transactions on Electrical Energy Systems, vol. 26, no. 11, pp. 2308-2321, Mar. 2016. [百度学术]
PJM. (2019, Dec.). PJM manual 28: operating agreement accounting. [Online]. Available: https://www.pjm.com/directory/manuals/m28/index.html#Sections/42%20Regulation%20Credits.html [百度学术]
E. Yao, V. W. S. Wong, and R. Schober, “Optimization of aggregate capacity of PEVs for frequency regulation service in day-ahead market,” IEEE Transactions on Smart Grid, vol. 9, no. 4, pp. 3519-3529, Dec. 2016. [百度学术]
PJM. (2019, Oct.). Performance score template. [Online]. Available: https://www.pjm.com/-/media/markets-ops/ancillary/mkt-based-regulation/performance-score-template.ashx?la=en [百度学术]
X. Jin, T. Jiang, Y. Mu et al., “Scheduling distributed energy resources and smart buildings of a microgrid via multi-time scale and model predictive control method,” IET Renewable Power Generation, vol. 13, no. 6, pp. 816-833, Apr. 2019. [百度学术]
H. K. Nunna and S. Doolla. “Multiagent based distributed energy resource management for intelligent microgrids,” IEEE Transactions on Industrial Electronics, vol. 60, no. 4, pp. 1678-1687, Apr. 2013. [百度学术]
J. Hu, S. You, M. Lind et al., “Coordinated charging of electric vehicles for congestion prevention in the distribution grid,” IEEE Transactions on Smart Grid, vol. 5, no. 2, pp. 703-711, Mar. 2014. [百度学术]
R. A. Verzijlbergh, L. J. De Vries, and Z. Lukszo, “Do we need congestion management in the distribution grid?” IEEE Transactions on Power Systems, vol. 29, no. 5, pp. 2119-2128, Mar. 2014. [百度学术]
S. Huang, Q. Wu, S. S. Oren et al., “Distribution locational marginal pricing through quadratic programming for congestion management in distribution networks,” IEEE Transactions on Power Systems, vol. 30, no. 4, pp. 2170-2178, Jul. 2015. [百度学术]
T. K. Kristoffersen, K. Capion, and P. Meibom, “Optimal charging of electric drive vehicles in a market environment,” Applied Energy, vol. 88, no. 5, pp. 1940-1948, May 2011. [百度学术]
Y. Zou, X. Hu, H. Ma et al., “Combined state of charge and state of health estimation over lithium-ion battery cell cycle lifespan for electric vehicles,” Journal of Power Sources, vol. 273, pp. 793-803, Jan. 2015. [百度学术]
L. Luo, W. Gu, X. Zhang et al., “Optimal siting and sizing of distributed generation in distribution systems with PV solar farm utilized as STATCOM (PV-STATCOM),” Applied Energy, vol. 210, pp. 1092-1100, Aug. 2017. [百度学术]
IBM. (2015, Jul.). IBM ILOG CPLEX optimization studio. [Online]. Available: https://www.ibm.com/products/ilog-cplex-optimization-studio [百度学术]