Figure 1.
The diagram of test system.
Citation: Kazushige Yokoyama, Christa D. Catalfamo, Minxuan Yuan. Reversible peptide oligomerization over nanoscale gold surfaces[J]. AIMS Biophysics, 2015, 2(4): 649-665. doi: 10.3934/biophy.2015.4.649
| [1] | Li Bin, Muhammad Shahzad, Qi Bing, Muhammad Ahsan, Muhammad U Shoukat, Hafiz MA Khan, Nabeel AM Fahal . The probabilistic load flow analysis by considering uncertainty with correlated loads and photovoltaic generation using Copula theory. AIMS Energy, 2018, 6(3): 414-435. doi: 10.3934/energy.2018.3.414 |
| [2] | Sonia Longo, Maurizio Cellura, Francesco Guarino, Vincenzo La Rocca, Giuseppe Maniscalco, Massimo Morale . Embodied energy and environmental impacts of a biomass boiler: a life cycle approach. AIMS Energy, 2015, 3(2): 214-226. doi: 10.3934/energy.2015.2.214 |
| [3] | Abdulrahman Th. Mohammad, Amer Ibrahim Ismael . An equivalent photovoltaic solar system to solve the problems of electricity in Iraqi houses. AIMS Energy, 2019, 7(5): 660-670. doi: 10.3934/energy.2019.5.660 |
| [4] | Ryuto Shigenobu, Oludamilare Bode Adewuyi, Atsushi Yona, Tomonobu Senjyu . Demand response strategy management with active and reactive power incentive in the smart grid: a two-level optimization approach. AIMS Energy, 2017, 5(3): 482-505. doi: 10.3934/energy.2017.3.482 |
| [5] | Mowffaq M. Oreijah, Mohammed Yunus . A parametric analysis to evaluate the performance metrics of power generation system involving Trilateral Flash Cycle using three different working fluids for low grade waste heat. AIMS Energy, 2019, 7(4): 483-492. doi: 10.3934/energy.2019.4.483 |
| [6] | Albert K. Awopone, Ahmed F. Zobaa . Analyses of optimum generation scenarios for sustainable power generation in Ghana. AIMS Energy, 2017, 5(2): 193-208. doi: 10.3934/energy.2017.2.193 |
| [7] | Abshir Ashour, Taib Iskandar Mohamad, Kamaruzzaman Sopian, Norasikin Ahmad Ludin, Khaled Alzahrani, Adnan Ibrahim . Performance optimization of a photovoltaic-diesel hybrid power system for Yanbu, Saudi Arabia. AIMS Energy, 2021, 9(6): 1260-1273. doi: 10.3934/energy.2021058 |
| [8] | Jinhua Tian, Yueyuan Zhang, Yanan Gao, Yu Qin, Bihan Fan, Cheng Zhang, Qiqi Zang . Distributed optimization and scheduling strategy for source load storage distribution grid based on alliance chain. AIMS Energy, 2024, 12(5): 946-967. doi: 10.3934/energy.2024044 |
| [9] | Khuthadzo Kgopana, Olawale Popoola . Improved utilization of hybrid energy for low-income houses based on energy consumption pattern. AIMS Energy, 2023, 11(1): 79-109. doi: 10.3934/energy.2023005 |
| [10] | Matteo Moncecchi, Davide Falabretti, Marco Merlo . Regional energy planning based on distribution grid hosting capacity. AIMS Energy, 2019, 7(3): 264-284. doi: 10.3934/energy.2019.3.264 |
The application of renewable sources (RS) contributes to the economic benefits and technical advances in the distribution system (DS). The RS are conducive when reducing the losses of power and voltage, improving flexibility and reliability, extending equipment life of the system as well as lessening the costs of purchasing energy from the market [1,2,3,4,5,6]. Additionally, the RS integrated in energy systems can reduce the pollutant emissions including CO, CO2, SO2 and NOx from traditional energy sources [7,8,9]. Therefore, the many researches have considered the application of RS in the planning and operating of DS with different approaches used for the optimization of RS generation and integration subproblems [10]. The long-term planning model is proposed to maximize the benefit of a project integrating the RS [11]. Similarly, the model planning the DS considering the distributed generators is also introduced with the objective is minimizing the investment and operation cost of distributed generators, cost to upgrade feeder and substations during the planning period, energy expenses purchased from markets and emission costs [12]. The planning models to simultaneously select optimal allocation of RS and capacitor banks in DS are also introduced in researches [13,14]. The objective function is minimizing the total cost of investment and operation of DS, pollutants emission cost of resources. The results showed that the RS have great effect in DS planning by reducing for power and energy loss, energy cost and emissions.
A new methodology is proposed to optimize the hybrid renewable energy system with various system configurations based on reliability and economic constraints [15]. The objective function is utilized to minimize the net present cost of project computed by HOMER software. Similarly, HOMER is also applied to optimize the configuration for grid-connected hybrid power systems considering RS such as PV and wind turbine [16]. The results is for computing the optimum system architecture and size of equipment, at the same time, the hybrid system reduces the significant amount of greenhouse gas emissions leading to greater environmental impact compared to the existing system. However, the above studies only compute with average load, electrical price, solar radiation and wind speed while they always change and are stochastic values. In addition, the lifetime and uptime of equipment on the planning period are not the same, the planning time is usually shorter than the equipment life cycle and the impact of reactive power is ignored in optimal problems.
Hence, life cycle cost (LCC) during computation time is introduced to analyze planning and operating of RS in energy systems [17,18,19,20]. Beside, the output power of wind turbine (WT) and photovoltaic (PV) as well as the load and electrical price are stochastic parameters that are usually expressed by the probability density function (pdf) in which the normal probability density function is chosen for modeling the uncertainties of the load and electrical price [21]. Similarly, the wind speed is modeled by the weibull probability density function and the beta probability density function is selected to present the irradiance data [22,23,24,25,26,27]. The output power of WT and PV depends on wind speeds and radiations, respectively, and thus those are stochastic parameters. The stochastic parameters are divided into levels (states) by clustering technique with specific values and different probabilities. In this context, a combined model are utilized to integrate the output power of WT and PV, the load and energy price [25,28,29]. The planning of the DS considering uncertainties is presented at several researches with stochastic parameters including the load, electrical, output power of the RS. The siting and sizing of the RS, as well as upgrading capacity and time of feeders, are concurrently chosen with the cost of the object at the present value [22,23,24,25,26,27,28,29,30,31] or power loss [32]. Similarly, the stochastic optimization approaches are also utilized to examine uncertainties in optimizing the configuration of energy systems. The objective is minimizing the cost of energy of over a particular time horizon [33] or the LCC which constitutes the capital cost, the replacement cost, the resale cost, and the operations and maintenance costs of the system [34]. However, the lifetime and uptime of the equipment have not been considered and so the error of the planning problem is large. Additionally, the optimal site selection of electric vehicle charging stations based on a bayesian network model is proposed, which simultaneously consider the uncertainty and the economic, environmental, and social criteria [35]. The optimal integration of electric vehicles and RS in a smart residential district is also introduced. The results showed that the amounts of power from RS reduce the annual energy cost of district [36].
Therefore, this research proposes a planning framework to analyze the optimal siting, sizing and technology of the RS with the objective function being minimum LCC during the planning period and considers uncertainties of the DS. The lifetime and uptime of equipment examined as well as the LCC objective function increase the accuracy of planning results. The LCC includes the energy cost, emission taxes, the investment and operation cost as well as the residual value of RS at the end of the calculation period. The technical parameters and economic effectiveness of the DS are determined and the normal operation of the DS is guaranteed by constraints. A combination model is utilized to integrate the uncertain parameters as demand, electrical price and output power of RS. The loss and operation parameters of the DS are analyzed at the all of states with the alternating current (AC) constraints contributing to the reduction of errors. The planning problem with uncertainties is very complex and thus GAMS/SCENRED is employed to quickly reduce the scenario number of the problem while guaranteeing the small error of calculation results as shown in researches [37,38,39,40].
Next section of the paper represents the modeling of uncertain parameters, the modeling problem with objective function and constraints presents in section 3. Section 4 shows results analyzed for the test system and the conclusion demonstrates in section 5.
The output power of RS changes random depending on natural intermittent of primary sources. Similarly, the energy price and electrical demand also are stochastic parameters. Hence, the pdf is often utilized to model uncertainties of parameters in DS as introduced in researchers [31,41]. Because the pdf is the continuous function, the clustering technique is utilized to divide stochastic parameters into different states. In each state, there is a specific value with the related probability.
The load at buses randomly change according to demand of consumers and is modeled by the normal probability density function as present at Eq 1 in which the µ is mean of the distribution, σ is the standard deviation and σ2 is variance [28,41,22,26]. The electricity price is also a stochastic parameter at the competitive market and expressed by the normal probability density function as Eq 1 [22,26].
| $ P(X = x|\mu , {\sigma ^2}) = f(x) = \frac{1}{{\sqrt {2\pi \sigma } }}\exp \left[ { - \frac{{{{(x - \mu )}^2}}}{{2{\sigma ^2}}}} \right] $ | (1) |
A value related with probability at each state of solar irradiance is expressed by beta probability density function as Eq 2 [29,44] in which the Iir is the solar irradiance, µ is mean and σ is the standard deviation of the stochastic variable.
|
$
\begin{array}{l} {f_b}({I_{ir}}) = \left\{ \begin{array}{l} \frac{{\Gamma (\alpha + \beta )}}{{\Gamma (\alpha )\Gamma (\beta )}}.{I_{ir}}^{(\alpha - 1)}.{(1 - {I_{ir}})^{(\beta - 1)}}\;{\mkern 1mu} if{\mkern 1mu}\; 0 \le {I_{ir}} \le 1,\alpha ,\beta \\ 0{\mkern 1mu} \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;else \end{array} \begin{array}{l} \right.\\
\beta = (1 - \mu ).\left( {\frac{{\mu .(1 + \mu )}}{{{\sigma ^2}}} - 1} \right);\;\;\;\;\;\;\alpha = \frac{{\mu .\beta }}{{1 - \mu }}
\end{array}
$
{f_b}({I_{ir}}) = \left\{ \begin{array}{l} \frac{{\Gamma (\alpha + \beta )}}{{\Gamma (\alpha )\Gamma (\beta )}}.{I_{ir}}^{(\alpha - 1)}.{(1 - {I_{ir}})^{(\beta - 1)}}\;{\mkern 1mu} if{\mkern 1mu}\; 0 \le {I_{ir}} \le 1,\alpha ,\beta \\ 0{\mkern 1mu} \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;else \end{array} |
(2) |
In each state, the output power is analyzed as expression (3) in which the kF.s is the fill factor and the n is the number of modules. The Is, Us are the generated current and voltage of the PV module at state s. Besides, the ISC and UOC are short circuit current and open circuit voltage, the UMPP and IMPP are current and the voltage at the maximum power point, respectively. The kU, kI are the coefficients depending on temperature, voltage and current of PV. The θA, θN are the ambient and nominal operating temperature of module PV.
Similarly, the wind speed also changes random represented under rayleigh probability density function which is a special case of weibull probability density function as expression (4) with the shape index k = 2 and the scale index c = 1.128vm [28,29,23].
| $ P_s^{PV}({I_{ir}}) = n.{k_{F.s}}.{U_s}({I_{ir}}).{I_s}({I_{ir}})\, \\ {k_{F.s}} = \frac{{{U_{MPP.s}}.{I_{MPP.s}}}}{{{U_{OC}}.{I_{SC}}}} \\ {U_s}({I_{ir}}) = {U_{OC}} - {k_U}.\left( {{\theta _A} + {I_{ir}}_{.s}.\frac{{{\theta _N} - 20}}{{0.8}}} \right) \\ {I_s}({I_{ir}}) = {I_{ir}}_{.s}.\left[ {{I_{SC}} + {k_I}.\left( {{\theta _A} + {I_{ir}}_{.s}.\frac{{{\theta _N} - 20}}{{0.8}} - 25} \right)} \right] \\ $ | (3) |
| $ {f_r}(v) = \left( {\frac{{2v}}{{{c^2}}}} \right)\exp \left[ { - {{\left( {\frac{v}{c}} \right)}^2}} \right] $ | (4) |
The value and probability of wind speed in each state are determined when the clustering technique is utilized to divide the pdf into the states. In this context, the output power of the WT is calculated as Eq 5.
|
$
P_{s}^{WT}(v) = \left\{ 0 vs≤vci or vco≤vsPWTr.vs−vcivr−vci vci≤vs≤vrPWTr vr≤vs≤vco \right.\,
$
|
(5) |
where, the $P_s^{WT}$ and $P_r^{WT}$ are output and rated power of WT. The vci, vcr, and vco are the cut-in speed, rated speed and cut-off speed of WT, respectively.
The multi-state parameters of the DS are examined by utilizing a combinable model to integrate all possible states of parameters [28,29]. Hence, the matrix M is represented as Eq 6 to enumerate all combination of states of the load, electricity price and output power of RS. In which, the Cs is the matrix that enumerates the possible values of parameters, the $\lambda _s^{}{\text{\{ }}{{\text{C}}_s}{\text{\} }}$ is one-column matrix that represents probability corresponding with matrix Cs and the Ns is the sum of individual states of the matrix M which is computed by product the possible states of the parameters.
|
$
M= { Cs, λs { Cs } } λs { Cs } = λρs { ρs } .λLs { PLs } .λWTs { PWTs } .λPVs { PPVs } Ns=nρs.nLs.nWTs.nPVs∀s∈Ns $
|
(6) |
where, the $\rho _s^{}$, $P_s^L$ are the electricity price and load corresponding with probabilities $\lambda _s^\rho $ and $\lambda _s^L$. The $P_s^{WT}$, $P_s^{PV}$ are the output power of WT and PV corresponding with probabilities $\, \, \lambda _s^{{\text{W}}T}$ and $\, \lambda _s^{PV}$, respectively. The $n_s^\lambda $, $n_s^L$ are the numbers of states assumed of the electricity price and load while the$n_s^{WT}, n_s^{PV}$are the numbers of states of the output power of RS, respectively.
The numbers of scenarios in matrix M are very big depending on choosing state number in each pdf. Therefore, the GAMS/SCENRED tool is employed to quickly reduce the scenario number of problem lead to decrease the computable volume and time while the results are guaranteed with the small error [42]. Beside, the choosing state numbers are very important because of its affect to the accuracy of results and the ability to calculate. The higher numbers of states selected, the higher accuracy is and increasing complexity of computation and vice versa.
When the renewable sources participate in the DS, the economic benefits and the electrical quality of this system are considerably enhanced. However, the invested cost of the RS is often high and siting, sizing as well as technique of the RS significant effect to efficiency of planning problem and thus a planning framework should be created to optimize the invested parameters of the RS in the DS planning. In this study, the mathematic model with the LCC objective function can examine the different lifetime and uptime of the RS in any planning period. Moreover, the binary variables and uncertainty of parameters in the DS added in objective and constrains greatly enhance the accuracy and appropriateness of investment decisions in practice. The mathematic model of the problem, which determines siting, sizing and locating of RS in DS, includes the objective and constrains detailed as follows.
The objective function of the model is the LCC of the project during the planning period as presented in Eq 7.
|
$
LCC=T∑t=11(1+r)t(NRS∑i=1NK∑k=1CRSP.k.PRSk,i,t+8760NK∑k=1NRS∑i=1NS∑s=1λs.ρRSOM.k.PRSk,i,t,s+8760NT∑i=1NS∑s=1kρ.λs(ρPs.PTi,t,s+ρQs.QTi,t,s)+8760NRS∑i=1NK∑k=1NS∑s=1β.(ξT−ξRS).λs.PRSk,i,t,s−TRSC.k−tRSinv.kTRSC.kNRS∑i=1NK∑k=1CRSP.k.PRSk,i,t)
$
|
(7) |
In Eq 7, the first and second parts are the investment and operation cost of RS, the third part is cost to purchase energy from the grid, the fourth part is emission taxes cost reduced due to low emissions of RS and the last part is residual value of equipment at the end of planning period. All of them are calculated at the base year by discount rate r.
Where, $C_{P.k}^{RS}$ is the invested cost of the RS with technology k; $P_{k, i, t, s}^{RS}$ is the out power variables of the RS; $\rho _{OM.k}^{RS}$ is operation and management (O & M) of the RS; $P_{i, t, s}^T, \, Q_{i, t, s}^T$ are the active and reactive powers received from main grid; $\rho _s^P, \, \rho _s^Q$ are electrical energies price with rise factor in each year ${k_\rho }$; ${\xi _T}, \, {\xi _{RS}}$ are CO2 emission coefficient of the traditional sources and RS with emission tax of CO2 is β; NT is total buses connected to main grid; NRS is total buses for installing the RS; k, NK are technology type and technologies set of the RS considered for selection; t, T are year and planning period; $T_{C.k}^{RS}, \, t_{inv.k}^{RS}$ are lifetime and invested time of the RS with technology k.
The nonlinear power flow model AC is selected to calculate both active and reactive power loss and thus the calculating operation cost of the DS is more accurate. They are reformed to integrate the RS as represented in Eq 8 with state index s of uncertainty parameters. The load at buses and the output power depending on technology k of RS are computed according to factors and rated powers present in Eqs 9 and 10. Similarly, the electricity prices are computed as Eq 11 with the peak price of active power is $\rho _{}^P$ and reactive power is$\rho _{}^Q$.
|
$
PTi,t,s−PLi,t,s+NK∑k=1PRSk,i,t,s=N∑j=1|Yij,t|.|Ui,t,s|.|Uj,t,s|.cos(θij,t−δj,t,s−δi,t,s)QTi,t,s−QLi,t,s+NK∑k=1QRSk,i,t,s=−N∑j=1|Yij,t|.|Ui,t,s|.|Uj,t,s|.sin(θij,t−δj,t,s−δi,t,s)∀i∈N,t∈T,s∈NS
$
|
(8) |
| $ P_{i, t, s}^{L} = k_{s}^{L}.P_{i, t}^{L}Q_{i, t, s}^{L} = k_{s}^{L}.Q_{i, t}^{L}\forall i\in N, \, \, \, t\in T, \, s\in {{N}_{S}}\, $ | (9) |
| $ P_{k, i, t, s}^{RS} = k_{k, s}^{RS}.P_{k, i, t}^{RS};\, \, \, \, \, \, Q_{k, i, t, s}^{RS} = P_{k, i, t, s}^{RS}.\tan {{\varphi }_{k}}\ \ \ \ \forall k\in {{N}_{K}}, \, i\in N, \, t\in T, \, s\in {{N}_{S}}\, $ | (10) |
| $ \rho _{s}^{P} = k_{s}^{\rho }.\rho _{{}}^{P} \;\;\; \rho _{s}^{Q} = k_{s}^{\rho }.\rho _{{}}^{Q} \;\;\; \forall s\in {{N}_{S}}\, $ | (11) |
where $P_{i, t, s}^L, \, Q_{i, t, s}^L$ are the active and reactive powers of the loads with demand factor at state s $k_s^L$; $\left| {{U_{i, t, s}}} \right|, \, {\delta _{i, t, s}}$are module and angle of bus voltages; $\left| {{Y_{ij, t}}} \right|, \, {\theta _{ij, t}}$ are magnitude and angle of admittance matrix element formulated by the feeder impedances in year t; The out power of the RS at state s is reflected by coefficient $k_{k, h}^{RS}$; $\rho _{}^P, \, \rho _{}^Q$ are maximum price of electrical energies at hour h and $k_s^\rho $ is coefficient analyzing at state s.
The binary variable ${\gamma _{{\text{k, i, t}}}}$ is utilized to limit sizing and power of RS that those can be chosen as shown in Eq 12. Each type of RS should only be invested one time at each bus during the planning period.
| $ P_{k, i, t}^{RS}\ge 0;\, \, \, P_{k, i, t}^{RS}-P_{k.\max }^{RS}.{{\gamma }_{\text{k, i, t}}}\le 0;\, \, \, \sum\limits_{t = 1}^{T}{{{\gamma }_{\text{i, k, t}}}\le 1} \;\;\; \forall k\in {{N}_{K}}, \, \, i\in N, \, t\in T\, $ | (12) |
Beside, the each feeder or connectable substation has a capacity limit for total power flow through it in all states. Because of that, the constraints (13) are proposed to guarantee that the equipment of the DS is not overloaded with the limited capacity of feeders $S_{ij.\max }^F$and connectable substation$S_{i.\max }^T$. Similarly, the voltage profile at buses must be guaranteed in a limit allowing the ordinary operation of the system as shown in Eq 14. At the substation bus connected to the grid, the voltage is often stability so it is assumed the constant. The voltage at load buses is usually varied under change of load. Hence, the voltage is limited as constrains (14) with minimum voltage profile ${U_{\min }}$ and maximum voltage profile ${U_{\max }}$.
| $ S_{ij, t, s}^{F}\le S_{ij.\max }^{F};\ \ S_{i, t, s}^{T}\le S_{i.\max }^{T}\ \ \ \forall \, ij\in N, \, t\in T, \, s\in S\, $ | (13) |
|
$
|Ui,t,s|=constant ∀i∈NT,t∈T,s∈NSUmin≤|Ui,t,s|≤Umax ∀i∈N,t∈T,s∈NS
$
|
(14) |
The proposed model is Mixed-Integer Nonlinear Programming (MINLP) and thus the GAMS/BONMIN solver is utilized to optimize the planning problem [40].
The radial diagram 9-bus with a connectable substation introduced by researches [13,15] as represented in Figure 1 is selected to investigate the feasibility and efficiency of the proposed model. The data of the system is modified to match the problem in which the normal voltage is 22 kV, the load at buses presents in the appendix and grows 5% annually. The peak demand total at the base year is 17.95 MW and 14.08 MVAr, respectively. A substation is connected to the main grid at bus 1 to feed energy for loads with maximum capacity 25 MVA.
Because of quick installation and small space, the candidate site chose to install the RS are the load buses. The power limit at each bus assumes 5 MW and the discrete step selected to calculate is 0.1 MW [28,29]. The cost of investment, average O & M and the analyzed parameter of RS are presented on table 1 [43,44]. Besides, the emission coefficient and tax of resources are also introduced on table 2 [12,31].
| Type of RS | Parameters | Value | Parameters | Value |
| PV | Peak power, Pr (kW) | 100 | kU (mV/℃); kI (mA/℃) | 14.4; 1.22 |
| UOC; UMPP (V) | 22; 17.32 | ISC; IMPP (A) | 5.3; 4.76 | |
| Life time (year) | 25 | θN; θA (℃) | 25; 43 | |
| Invested cost (M＄/MW) | 1.5 | O & M cost ($/MWh) | 20 | |
| WT | Rate power, Pr (kW) | 100 | vr (m/s) | 14 |
| vci (m/s) | 4 | vco (m/s) | 24 | |
| Invested cost (M＄/MW) | 1.8 | O & M cost ($/MWh) | 40 | |
| Life time (year) | 25 |
DownLoad:
CSV
| Parameters | Value |
| Emission tax of CO2, β | 10 ＄/ton |
| Emission coefficient of traditional sources, ξT | 0.68 kg/kWh |
| Emission coefficient of RS, ξRS | 0.02 kg/kWh |
| Discount rate, r | 10% |
DownLoad:
CSV
The voltage at the substation bus is assumed the constant equaling 1.05 pu while the voltages at load buses allow change from 0.9 pu to 1.1 pu.
The stochastic parameters, which are conceptualized under the pdf, are divided by clustering technique into states under the annual data as shown in Table 3. The 15 levels (statuses) with the specific value and related probability to guarantee a reasonable balance between accuracy and computed speed are selected to analyze in this research [28,29,22,23,24,25,26].
| State | PV | WT | Load | Electrical price | ||||
| $k_s^{RS}$ | $\lambda _s^{PV}$ | $k_s^{RS}$ | $\lambda _s^{WT}$ | $k_s^L$ | $\lambda _s^L$ | $k_s^\rho $ | $\lambda _s^\rho $ | |
| 1 | 0 | 0.450 | 0 | 0.199 | 0.30 | 0.001 | 0.30 | 0.012 |
| 2 | 0.06 | 0.002 | 0.1 | 0.112 | 0.35 | 0.003 | 0.35 | 0.023 |
| 3 | 0.12 | 0.005 | 0.2 | 0.113 | 0.40 | 0.010 | 0.40 | 0.041 |
| 4 | 0.18 | 0.012 | 0.3 | 0.106 | 0.45 | 0.028 | 0.45 | 0.065 |
| 5 | 0.23 | 0.023 | 0.4 | 0.096 | 0.50 | 0.062 | 0.50 | 0.092 |
| 6 | 0.29 | 0.039 | 0.5 | 0.085 | 0.55 | 0.111 | 0.55 | 0.116 |
| 7 | 0.35 | 0.057 | 0.6 | 0.072 | 0.60 | 0.159 | 0.60 | 0.131 |
| 8 | 0.40 | 0.073 | 0.7 | 0.059 | 0.65 | 0.184 | 0.65 | 0.132 |
| 9 | 0.45 | 0.083 | 0.8 | 0.047 | 0.70 | 0.172 | 0.70 | 0.119 |
| 10 | 0.50 | 0.084 | 0.9 | 0.036 | 0.75 | 0.130 | 0.75 | 0.096 |
| 11 | 0.55 | 0.072 | 1.0 | 0.027 | 0.80 | 0.079 | 0.80 | 0.070 |
| 12 | 0.60 | 0.053 | 1.0 | 0.020 | 0.85 | 0.039 | 0.85 | 0.045 |
| 13 | 0.65 | 0.031 | 1.0 | 0.014 | 0.90 | 0.015 | 0.90 | 0.026 |
| 14 | 0.70 | 0.013 | 1.0 | 0.009 | 0.95 | 0.005 | 0.95 | 0.014 |
| 15 | 0.75 | 0.003 | 1.0 | 0.006 | 1.00 | 0.001 | 1.00 | 0.006 |
DownLoad:
CSV
The agammaegation of uncertainty parameters is made by matrix M with 154 (50625) scenarios, the GAMS/SCENRED tool utilized to reduce number of scenarios to 15 [42].
The GAMS/BONMIN solver is applied to optimize the planning problem. The simulation result determines the siting, sizing and time invested of each RS as table 4. The WT and PV are chosen to install at buses, which are far connectable substation as buses 7, 8 and 9 because they decrease the loss and energy cost. The total power of PV chose is 7.0MW in the first year to salvage the low operation cost of this source. Similarly, the 8.1MW is the rated power of WT installed in DS at the 6th year.
| Type of RS | Siting (bus) | Rated power (MW) | Planning year |
| PV | 7 | 2.3 | 1 |
| 8 | 2.3 | 1 | |
| 9 | 2.4 | 1 | |
| WT | 7 | 2.7 | 6 |
| 8 | 2.6 | 6 | |
| 9 | 2.8 | 6 |
DownLoad:
CSV
The voltage profile at buses is improved because of loss in the reduction of DS as shown in Figure 2. The lowest voltage profile when without RS is 0.88 pu at bus 7 of the 10th planning year. Therefore, the electrical quality does not guarantee supply to consumers. Conversely, the voltage profiles at all buses are guaranteed during planning period with value being greater than 0.92 pu in case with RS.
The power and energy loss also decreased due to a decrease in power flow in the DS as presented in Figure 3. The PV installed at first year resulted in the decrease of maximal power loss at beginning of planning time by 2%, and increase by 3.97% at the 10th year when the WT was applied at the 6th year. Besides, the RS reduces energy loss of the DS as well as energy cost from the market. The electrical energy obtained from the market decreased from 1,098,100.0 MWh to 858,310.0 MWh in correspondence with 239,790.0 MWh or 21.84%. The CO2 emission of the RS equals only 2.9% of the traditional sources or the reduction is about 0.66 tons/MWh and thus the CO2 emission drops 158,261.4 tons at the case with RS as shown in table 5. In this case, the emission from the traditional sources is 583,650.8 tons while from the RS is only 4,795.8 tons.
| Case | Total life cycle cost, M＄ | Emission taxes cost, M＄ | CO2 emission, ton |
| Without RS | 82.92 | 5.34 | 746,708.0 |
| With RS | 76.58 | 4.22 | 588,446.6 |
| Comparison of two cases | 6.34 | 1.12 | 158,261.4 |
DownLoad:
CSV
As presented in table 5, the effect of investment project considering the RS is verified along with the reduction of LCC. The LCC decreased by 7.65% (6.34 M＄) when the RS is involved. Similarly, the emission taxes in the case of RS is only 4.22 M＄, 1.12 M＄ less compared to the situation without RS. The simulation result shows the feasibility and effectiveness of the proposed model as well as the benefit brought by RS to the DS when the stochastic parameters are included in the problem.
This research introduced a model to optimize the invested siting, sizing, setting and time of the RS in the DS considering the uncertain parameters. The developed model features a number of important aspects:
● The planning framework proposed with the objective function that minimizes the LCC of investment project including the investment and operation cost of the RS, the cost for purchasing the energy from market, the emission taxes cost and the residual value of equipment at the end of the analysis period. Therefore, the different lifetime and uptime of the RS are examined which improves the accuracy and suitability for the practical planning problem.
● A combinable model with the clustering technique is utilized to integrate the multi-state uncertain parameters of planning problem as the electrical price, load, out power of the PV and WT.
● The binary variables are integrated in MINLP model to perform the invested decision of the RS with limited power and selected power, which agree with the actual parameters.
● The simulation results demonstrate how the application of the RS in DS could promote the reduction of LCC, and minimize power losses and emission cost during the planning stage. The voltage profile at buses and the flexibility of system operation could also be improved when the RS is considered when investing in DS planning problem.
The cases study has illustrated the feasibility and effectiveness of the proposed model with uncertainty parameters of the DS, which indicates that the optimization method could be brought into practice. However, upgrading the feeders and substations during planning period can be deferred when RS are selected to invest because it reduces the power transmitted from the utility grid. The simultaneous selection the siting, sizing of the RS and upgrading the feeders and substations of the DS would also be an interesting work for future research.
| Bus i-j | $S_{ij.\max }^F$(MVA) | Lij (km) | Rij (Ω) | Xij (Ω) |
| 1-2 | 6.67 | 5.3 | 4.10 | 2.27 |
| 1-3 | 6.67 | 9.2 | 7.11 | 3.95 |
| 1-4 | 6.67 | 10.3 | 7.96 | 4.42 |
| 1-5 | 6.67 | 10.5 | 8.12 | 4.50 |
| 2-6 | 6.67 | 11.7 | 9.04 | 5.02 |
| 3-7 | 6.67 | 18.2 | 14.07 | 7.81 |
| 4-8 | 6.67 | 17.2 | 13.30 | 7.38 |
| 5-9 | 6.67 | 15.2 | 11.75 | 6.52 |
DownLoad:
CSV
| Bus | $P_i^L$(MW) | $Q_i^L$(MVAr) | Bus | $P_i^L$(MW) | $Q_i^L$(MVAr) |
| 1 | - | - | 6 | 1.58 | 1.37 |
| 2 | 2.83 | 2.29 | 7 | 1.31 | 1.14 |
| 3 | 2.78 | 2.38 | 8 | 1.16 | 0.98 |
| 4 | 2.95 | 1.99 | 9 | 1.22 | 1.05 |
| 5 | 2.86 | 1.90 | Total | 16.69 | 13.10 |
DownLoad:
CSV
The research is supported by Thai Nguyen University of Technology (TNUT), Thai Nguyen, Viet Nam.
The author declares no conflict of interest in this paper.
| [1] | Zahs KR, Ashe KH (2013) β-Amyloid oligomers in aging and Alzheimer's disease. Front Aging Neurosci 5: 28. |
| [2] | Ramamoorthy A, Lim MH (2013) Structural Characterization and Inhibition of Toxic Amyloid-β Oligomeric Intermediates. Biophys J 105: 287-288. |
| [3] | Pham JD, Chim N, Goulding CW, et al. (2013) Structures of Oligomers of a Peptide from β-Amyloid. J Am Chem Soc 135: 12460-12467. |
| [4] |
Hendrik H, Farmer BL, Pandey RB, et al. (2009) Nature of molecular interactions of peptides with gold, palladium, and Pd-Au bimetal surfaces in aqueous solution. J Am Chem Soc 131: 9704-9714. doi: 10.1021/ja900531f
|
| [5] | Nooren IM, Thornton JM (2003) Structural characterisation and functional significance of transient protein-protein interactions. J Mol Biol 325: 991-1018. |
| [6] | Nooren IM, Thornton JM (2003) Diversity of protein-protein interactions. EMBO J 14: 3486-3492. |
| [7] |
Goodsell DS, Olson AJ (1993) Soluble Proteins: Size, Shape and Function. Trends Biochem Sci 18: 65-68. doi: 10.1016/0968-0004(93)90153-E
|
| [8] | Long M, Betrán E, Thornton K, et al. (2003) The origin of new genes: glimpses from the young and old. Nat Rev Genet 4: 865-875. |
| [9] | Chothia C, Gough J, Vogel C, et al. (2003) Evolution of the protein repertoire. Science 300: 1701-1703. |
| [10] | Liu Y, Eisenberg D (2002) 3D domain swapping: as domains continue to swap. Protein Sci 11: 1285-1299. |
| [11] |
Gill S, Beck K (1965) Differential heat capacity calorimeter for polymer transition studies. Rev Sci Instrum 36: 274. doi: 10.1063/1.1719552
|
| [12] | Goodsell DS, Olson AJ (2000) Structural Symmetry and Protein Function. Annu Rev Biophys Biomol Struct 29: 105-153. |
| [13] |
Jones S, Thornton JM (1996) Principles of protein-protein interactions. Proc Natl Acad Sci U S A 93: 13-20. doi: 10.1073/pnas.93.1.13
|
| [14] |
Searle MS, Ciani B (2004) Design of β-sheet systems for understanding the thermodynamics and kinetics of protein folding. Curr Opin Struct Biol 14: 458-464. doi: 10.1016/j.sbi.2004.06.001
|
| [15] |
Snow CD, Nguyen H, Pande VS, et al. (2002) Absolute comparison of simulated and experimental protein-folding dynamics. Nature 420: 102-106. doi: 10.1038/nature01160
|
| [16] | Zagrovic B, Pande VS (2003) Structural correspondence between the alpha-helix... proteins can have native-like properties. Nat Struct Biol 10: 955-961. |
| [17] | Qiu L, Pabit SA, Roitberg AE, et al. (2002) Smaller and Faster: The 20-Residue Trp-Cage Protein Folds in 4 μs. J Am Chem Soc 124: 12952-12953. |
| [18] | Snow CD, Zagrovic B, Pande VS (2002) Folding kinetics and unfolded state topology via molecular dynamics simulations. J Am Chem Soc 124: 14548-14549. |
| [19] | Hartl FU, Hayer-Hartl M (2002) Molecular chaperones in the cytosol: from nascent chain to folded protein. Science 295: 1852-1858. |
| [20] | Dobson CM (2003) Protein folding and misfolding. Nature 426: 884-890. |
| [21] | Dobson CM ( 2004) Principles of protein folding, misfolding and aggregation. Semin Cell Dev Biol 15: 3-16. |
| [22] | Daggett V, Fersht A (2003) The present view of the mechanism of protein folding. Nat Rev Mol Cell Biol 4: 497-502. |
| [23] |
Shea J-E, Brooks CL, III (2001) From folding theories to folding proteins: a review and assessment of simulation studies of protein folding and unfolding. Annu Rev Phys Chem 52: 499-535. doi: 10.1146/annurev.physchem.52.1.499
|
| [24] | Saunders JA, Scheraga HA (2003) Ab initio structure prediction of two alpha-helical oligomers with a multiple-chain-united residue force field and global search. Biopolymers 68: 300-317. |
| [25] |
Zitzewitz JA, Bilsel O, Luo J, et al. (1995) Probing the folding mechanism of a leucine zipper peptide by stopped-flow circular dichroism spectroscopy. Biochemistry 34: 12812-12819. doi: 10.1021/bi00039a042
|
| [26] |
Dürr E, Jelesarov I, Bosshard HR (1999) Extremely fast folding of a very stable leucine zipper with a strengthened hydrophobic core and lacking electrostatic interactions between helices. Biochemistry 38: 870-880. doi: 10.1021/bi981891e
|
| [27] | Ibarra-Molero B, Makhatadze GI, Matthews CR (2001) Mapping the energy surface for the folding reaction of the coiled-coil peptide GCN4-p1 Biochemistry 40: 719-731. |
| [28] |
Zhu H, Celinski SA, Scholtz JM, et al. (2001) An engineered leucine zipper a position mutant with an unusual three-state unfolding pathway. Protein Sci 10: 24-33. doi: 10.1110/ps.30901
|
| [29] |
Wendt H, Leder L, Härmä H, et al. (1997) Very rapid, ionic strength-dependent association and folding of a heterodimeric leucine zipper. Biochemistry 36: 204-213. doi: 10.1021/bi961672y
|
| [30] |
Jang H, Hall CK, Zhou Y (2004) Thermodynamics and stability of a beta-sheet complex: molecular dynamics simulations on simplified off-lattice protein models. Protein Sci 13: 40-53. doi: 10.1110/ps.03162804
|
| [31] |
Dürr E, Jelesarov I (2000) Thermodynamic Analysis of Cavity Creating Mutations in an Engineered Leucine Zipper and Energetics of Glycerol-Induced Coiled Coil Stabilization. Biochemistry 39: 4472-4482. doi: 10.1021/bi992948f
|
| [32] |
Norde W (2008) My voyage of discovery to proteins in flatland and beyond. Colloids and Surfaces B: Biointerfaces 61: 1-9. doi: 10.1016/j.colsurfb.2007.09.029
|
| [33] | Sarkar B, Das AK, Maiti S (2013 ) Thermodynamically stable amyloid-β monomers have much lower membrane affinity than the small oligomers. Front Physiol 4: 84. |
| [34] | Zhu X, Yan D, Fang Y (2001) In situ FTIR spectroscopic study of the conformational change of isotactic polypropylene during the crystallization process. J Phys Chem B 105: 12461-12463. |
| [35] | Harper SM, Neil LC, Gardner KH (2003) Structural basis of a phototropin light switch. Science 301: 1541-1544. |
| [36] | Ohba S, Hosomi H, Ito Y (2001) In situ X-ray observation of pedal-like conformational change and dimerization of trans-cinnamamide in cocrystals with phthalic acid. J Am Chem Soc 123: 6349-6352. |
| [37] | Gupta R, Ahmad F (1999) Protein stability: functional dependence of denaturational gibbs energy on urea concentration. Biochemistry 38: 2471-2479. |
| [38] | Aisenbrey C, Borowik T, Bystrom R, et al. (2008) How is protein aggregation in amyloidogenic diseases modulated by biological membranes? Eur Biophys J 37: 247-255. |
| [39] | Lepère M, Chevallard C, Hernandez JF, et al. (2007) Multiscale surface self-assembly of an amyloid-like peptide. Langmuir 23: 8150-8155. |
| [40] | Kusumoto Y, Lomakin A, Teplow DB, et al. (1998) Temperature dependence of amyloid beta-protein fibrillization. Proc Natl Acad Sci U S A 95: 12277-12282. |
| [41] | Coles M, Bicknell W, Watson AA, et al. (1998) Solution structure of amyloid beta-peptide (1-40) in a water-micelle environment. Is the membrane spanning domain where we think it is?. Biochemistry 37: 11064-11077. |
| [42] | Shao HY, Jao SC, Ma K, et al. (1999) Solution structures of micelle-bound amylois beta-(1-40) and beta-(1-42) peptides of Alzhimer's disease. J Mol Biol 285: 755-773. |
| [43] | Giacomelli CE, Norde W (2005) Conformational changes of the amyloid beta-peptide (1-40) adsorbed on solid surfaces. Macromolecular Bioscience 5: 401-407. |
| [44] | Rocha S, Krastev R, Thunemann AF, et al. (2005) Adsorption of amyloid beta-peptide at polymer surfaces: a neutron reflectivity study. Chem Phys Chem 6: 2527-2534. |
| [45] | Kowalewski T, Holtzman DM (1999) In situ atomic force microscopy study of Alzheimer's beta-amyloid peptide on different substrates: new insights into mechanism of beta-sheet formation. Proc Natl Acad Sci U S A 96: 3688-3693. |
| [46] | Schladitz C, Vieira EP, Hermel H, et al. (1999) Amyloid-beta-sheet formation at the air-water interface. Biophys J 77: 3305-3310. |
| [47] | Rocha S, Thünemann AF, Pereira MC, et al. (2008) Influence of fluorinated and hydrogenated nanoparticles on the structure and fibrillogenesis of amyloid beta-peptide. Biophys Chem 137: 35-42. |
| [48] | Yokoyama K (2010) Nanoscale Surface Size Dependence in Protein Conjugation. In: Chen EJ, Peng N, editors. Advances in Nanotechnology: Nova Publisher, 65-104. |
| [49] | Yokoyama K (2011) Modeling of Reversible Protein Conjugation on Nanoscale Surface. In: Musa SM, editor. Computational Nanotechnology: Modeling and Applications with MATLAB: CRC Press-Taylor and Francis Group, LLC. |
| [50] | Yokoyama K (2012) Nano Size Dependent Properties of Colloidal Surfaces. In: Ray PC, editor. Colloids: Classification, Properties and Applications: Nova Science Publishing. |
| [51] | Yokoyama K (2010) Nanoscale Surface Size Dependence in Protein Conjugation In: Chen EJ, Peng N, editors. Advances in Nanotechnology: Nova Publisher, 65-104. |
| [52] | Yokoyama K (2011) Modeling of Reversible Protein Conjugation on Nanoscale Surface In: Musa SM, editor. Computational Nanotechnology: Modeling and Applications with MATLAB: CRC Press-Taylor and Francis Group, LLC. |
| [53] | Yokoyama K (2012) Nano Size Dependent Properties of Colloidal Surfaces In: Ray PC, editor. Colloids: Classification, Properties and Applications Nova Science Publishing. |
| [54] | Yokoyama K, Briglio NM, Sri Hartati D, et al. (2008) Nanoscale size dependence in the conjugation of amyloid beta and ovalbumin proteins on the surface of gold colloidal particles. Nanotechnology 19: 375101-375108. |
| [55] | Yokoyama K, Welchons DR (2007) The conjugation of amyloid beta protein on the gold colloidal nanoparticles' surfaces. Nanotechnology 18: 105101-105107. |
| [56] | Andrieux K, Couvreur P (2013) Nanomedicine as a promising approach for the treatment and diagnosis of brain diseases: the example of Alzheimer's disease. Ann Pharm Fr 71: 225-233. |
| [57] |
Politi J, Spadavecchia J, Iodice M, et al. (2015) Oligopeptide-heavy metal interaction monitoring by hybrid gold nanoparticle based assay. Analyst 140: 149-155. doi: 10.1039/C4AN01491J
|
| [58] | Yokoyama K (2011) Modeling of Reversible Protein Conjugation on Nanoscale Surface In: Musa SM, editor. Computational Nanotechnology: Modeling and Applications with MATLAB: CRC Press, 381-409. |
| [59] | Yokoyama K (2012) Nano Size Dependent Properties of Colloidal Surfaces. In: Ray PC, editor. Colloids: Classification, Properties and Applications: Nova Publishers, 25-58 |
| [60] | Yokoyama K (2013) Controlling Reversible Self Assembly Path of Amyloid Beta Peptide over Gold Colloidal nanoparticle’s Surface. In: Musa SM, editor. Nanoscale Spectroscopy with Applications CRC Press-Taylor and Francis Group, LLC, 279-304. |
| [61] | Link S, El-Sayed M (1999) Spectral properties and relaxation dynamics of surface plasmon electronic oscillations in gold and silver nanodots and nanorods J Phys Chem B 103: 8410-8426. |
| [62] | Kelly KL, Coronado E, Zhao LL, et al. (2003) The optical properties of metal nanoparticles: the influence of size, shape, and dielectric environment J Phys Chem B 10:7 668-677. |
| [63] | Jensen TR, Schatz GC, Van Duyne RP (1999) Nanosphere lithography: surface plasmon resonance spectrum of a periodic array of silver nanoparticles by ultraviolet-visible extinction spectroscopy and electrodynamic modeling. J Phys Chem B 103: 2394-2401. |
| [64] | Ackermann T (1969) The calorimeters: adiabatic calorimeters; Brown HD, Ed., editor. New York: Academic Press. |
| [65] | Barrow CJ, Yasuda A, Kenny PT, et al. (1992) Solution conformations and aggregational properties of synthetic amyloid beta-peptides of Alzheimer's disease. Analysis of circular dichroism spectra. J Mol Biol 225: 1075-1093. |
| [66] | Wood SJ, Maleeff B, Hart T, et al. (1996) Physical, morphological and functional differences between pH 5.8 and 7.4 aggregates of the Alzheimer's amyloid peptide Abeta. J Mol Biol 256: 870-877. |
| [67] | Jang S, Shin S (2008) Computational study on the structural diversity of amyloid beta-peptide (Aβ10-35) oligomers. J Phys Chem B 112: 3479-3484. |
| [68] | Nichols MR, Moss MA, Reed DK, et al. (2005) Amyloid-beta aggregates formed at polar-nonpolar interfaces differ from amyloid-beta protofibrils produced in aqueous buffers. Microsc Res Tech 67: 164-174. |
| [69] |
Bucciantini M, Giannoni E, Chiti F, et al. (2002) Inherent toxicity of aggregates implies a common mechanism for protein misfolding diseases. Nature 416: 507-511. doi: 10.1038/416507a
|
| [70] |
Brovchenko I, Burri RR, Krukau A, et al. (2009) Thermal expansivity of amyloid β16-22 peptides and their aggregates in water. Phys Chem Chem Phys 11 11: 5035-5040. doi: 10.1039/b820340g
|
| [71] | Melquiond A, Dong X, Mousseau N, et al. (2008) Role of the region 23-28 in Abeta fibril formation. Curr Alzheimer Res 5: 244-250. |
| [72] | Meier M, Seelig J (2008) Length dependence of the coil - β-sheet transition in a membrane environement. J Am Chem Soc 130: 1017-1024. |
| [73] |
Yokoyama K, Gaulin NB, Cho H, et al. (2010) Temperature Dependence of Conjugation of Amyloid Beta Peptide on the Gold Colloidal Nanoparticles. J Phys Chem A 114: 1521-1528. doi: 10.1021/jp907880f
|
| [74] |
Sabate R, Estelrich J (2003) Pinacyanol as effective probe of fibrillar beta-amyloid peptide; Comparative study with Congo Red. Biopolymers 72: 455-463. doi: 10.1002/bip.10485
|
| [75] | Carrotta R, Manno M, Bulone D, et al. (2005) Protofibril formation of amyloid beta-protein at low pH via a non-cooperative elongation mechanism. J Biol Chem 280: 30001-30008. |
| 1. | Nooriya A. Mohammed, Ammar Al-Bazi, Management of renewable energy production and distribution planning using agent-based modelling, 2021, 164, 09601481, 509, 10.1016/j.renene.2020.08.159 | |
| 2. | Abdul Matin Ibrahimi, Mohammad Masih Sediqi, Harun Or Rashid Howlader, Mir Sayed Shah Danish, Shantanu Chakraborty, Tomonobu Senjyu, Generation expansion planning considering renewable energy integration and optimal unit commitment: A case study of Afghanistan, 2019, 7, 2333-8334, 441, 10.3934/energy.2019.4.441 | |
| 3. | Merlin Sajini M.L., Suja S., Merlin Gilbert Raj S., Impact analysis of time‐varying voltage‐dependent load models on hybrid DG planning in a radial distribution system using analytical approach, 2021, 15, 1752-1416, 153, 10.1049/rpg2.12013 | |
| 4. | Angelo Raciti, Santi Agatino Rizzo, Giovanni Susinni, Circuit model of LED light bulb suitable for typical voltage THD on LV distribution networks, 2020, 212, 03787788, 109665, 10.1016/j.enbuild.2019.109665 | |
| 5. | 2020, Bi-Level Distribution Network Planning Integrated with Energy Storage to PV-Connected Network, 978-1-7281-5635-4, 1325, 10.1109/ISIE45063.2020.9152508 | |
| 6. | Vu Van Thang, Nguyen Hien Trung, 2020, Chapter 61, 978-3-030-37496-9, 534, 10.1007/978-3-030-37497-6_61 | |
| 7. | V. V. Thang, N. H. Trung, 2021, Chapter 45, 978-3-030-64718-6, 399, 10.1007/978-3-030-64719-3_45 | |
| 8. | Mahmoud M. Gamil, Mohammed Elsayed Lotfy, Ashraf M. Hemeida, Paras Mandal, Hiroshi Takahashi, Tomonobu Senjyu, Optimal sizing of a residential microgrid in Egypt under deterministic and stochastic conditions with PV/WG/Biomass Energy integration, 2021, 9, 2333-8334, 483, 10.3934/energy.2021024 | |
| 9. | V. V. Thang, N. H. Trung, Probabilistic optimization of planning and operation of networked microgrids with renewable energy resources considering demand response programs, 2022, 1868-3967, 10.1007/s12667-022-00559-8 | |
| 10. | Abdolmajid Erfani, Mehdi Tavakolan, Ali Hassandokht Mashhadi, Pouria Mohammadi, Heterogeneous or homogeneous? A modified decision-making approach in renewable energy investment projects, 2021, 9, 2333-8334, 558, 10.3934/energy.2021027 | |
| 11. | Zahir Barahmand, Marianne S. Eikeland, Techno-Economic and Life Cycle Cost Analysis through the Lens of Uncertainty: A Scoping Review, 2022, 14, 2071-1050, 12191, 10.3390/su141912191 | |
| 12. | Aprajay Verma, K Shanti Swarup, An analysis of distribution planning under a regulatory regime: An integrated framework, 2023, 4, 2634-1581, 179, 10.1049/enc2.12088 | |
| 13. | Mohammad Ahmad A. Al-Ja’Afreh, Bilal Amjad, Kirkland Rowe, Geev Mokryani, Jorge L. Angarita Marquez, Optimal planning and forecasting of active distribution networks using a multi-stage deep learning based technique, 2023, 10, 23524847, 686, 10.1016/j.egyr.2023.07.014 | |
| 14. | Alejandro Valencia-Díaz, Sebastián García H., Ricardo A. Hincapie I., Stochastic optimization for siting and sizing of renewable distributed generation and D-STATCOMs, 2025, 27726711, 101026, 10.1016/j.prime.2025.101026 |
Figures(8) / Tables(1)
Kazushige Yokoyama, Christa D. Catalfamo, Minxuan Yuan. Reversible peptide oligomerization over nanoscale gold surfaces[J]. AIMS Biophysics, 2015, 2(4): 649-665. doi: 10.3934/biophy.2015.4.649
| Type of RS | Parameters | Value | Parameters | Value |
| PV | Peak power, Pr (kW) | 100 | kU (mV/℃); kI (mA/℃) | 14.4; 1.22 |
| UOC; UMPP (V) | 22; 17.32 | ISC; IMPP (A) | 5.3; 4.76 | |
| Life time (year) | 25 | θN; θA (℃) | 25; 43 | |
| Invested cost (M＄/MW) | 1.5 | O & M cost ($/MWh) | 20 | |
| WT | Rate power, Pr (kW) | 100 | vr (m/s) | 14 |
| vci (m/s) | 4 | vco (m/s) | 24 | |
| Invested cost (M＄/MW) | 1.8 | O & M cost ($/MWh) | 40 | |
| Life time (year) | 25 |
DownLoad:
CSV
| Parameters | Value |
| Emission tax of CO2, β | 10 ＄/ton |
| Emission coefficient of traditional sources, ξT | 0.68 kg/kWh |
| Emission coefficient of RS, ξRS | 0.02 kg/kWh |
| Discount rate, r | 10% |
DownLoad:
CSV
| State | PV | WT | Load | Electrical price | ||||
| $k_s^{RS}$ | $\lambda _s^{PV}$ | $k_s^{RS}$ | $\lambda _s^{WT}$ | $k_s^L$ | $\lambda _s^L$ | $k_s^\rho $ | $\lambda _s^\rho $ | |
| 1 | 0 | 0.450 | 0 | 0.199 | 0.30 | 0.001 | 0.30 | 0.012 |
| 2 | 0.06 | 0.002 | 0.1 | 0.112 | 0.35 | 0.003 | 0.35 | 0.023 |
| 3 | 0.12 | 0.005 | 0.2 | 0.113 | 0.40 | 0.010 | 0.40 | 0.041 |
| 4 | 0.18 | 0.012 | 0.3 | 0.106 | 0.45 | 0.028 | 0.45 | 0.065 |
| 5 | 0.23 | 0.023 | 0.4 | 0.096 | 0.50 | 0.062 | 0.50 | 0.092 |
| 6 | 0.29 | 0.039 | 0.5 | 0.085 | 0.55 | 0.111 | 0.55 | 0.116 |
| 7 | 0.35 | 0.057 | 0.6 | 0.072 | 0.60 | 0.159 | 0.60 | 0.131 |
| 8 | 0.40 | 0.073 | 0.7 | 0.059 | 0.65 | 0.184 | 0.65 | 0.132 |
| 9 | 0.45 | 0.083 | 0.8 | 0.047 | 0.70 | 0.172 | 0.70 | 0.119 |
| 10 | 0.50 | 0.084 | 0.9 | 0.036 | 0.75 | 0.130 | 0.75 | 0.096 |
| 11 | 0.55 | 0.072 | 1.0 | 0.027 | 0.80 | 0.079 | 0.80 | 0.070 |
| 12 | 0.60 | 0.053 | 1.0 | 0.020 | 0.85 | 0.039 | 0.85 | 0.045 |
| 13 | 0.65 | 0.031 | 1.0 | 0.014 | 0.90 | 0.015 | 0.90 | 0.026 |
| 14 | 0.70 | 0.013 | 1.0 | 0.009 | 0.95 | 0.005 | 0.95 | 0.014 |
| 15 | 0.75 | 0.003 | 1.0 | 0.006 | 1.00 | 0.001 | 1.00 | 0.006 |
DownLoad:
CSV
| Type of RS | Siting (bus) | Rated power (MW) | Planning year |
| PV | 7 | 2.3 | 1 |
| 8 | 2.3 | 1 | |
| 9 | 2.4 | 1 | |
| WT | 7 | 2.7 | 6 |
| 8 | 2.6 | 6 | |
| 9 | 2.8 | 6 |
DownLoad:
CSV
| Case | Total life cycle cost, M＄ | Emission taxes cost, M＄ | CO2 emission, ton |
| Without RS | 82.92 | 5.34 | 746,708.0 |
| With RS | 76.58 | 4.22 | 588,446.6 |
| Comparison of two cases | 6.34 | 1.12 | 158,261.4 |
DownLoad:
CSV
| Bus i-j | $S_{ij.\max }^F$(MVA) | Lij (km) | Rij (Ω) | Xij (Ω) |
| 1-2 | 6.67 | 5.3 | 4.10 | 2.27 |
| 1-3 | 6.67 | 9.2 | 7.11 | 3.95 |
| 1-4 | 6.67 | 10.3 | 7.96 | 4.42 |
| 1-5 | 6.67 | 10.5 | 8.12 | 4.50 |
| 2-6 | 6.67 | 11.7 | 9.04 | 5.02 |
| 3-7 | 6.67 | 18.2 | 14.07 | 7.81 |
| 4-8 | 6.67 | 17.2 | 13.30 | 7.38 |
| 5-9 | 6.67 | 15.2 | 11.75 | 6.52 |
DownLoad:
CSV
| Bus | $P_i^L$(MW) | $Q_i^L$(MVAr) | Bus | $P_i^L$(MW) | $Q_i^L$(MVAr) |
| 1 | - | - | 6 | 1.58 | 1.37 |
| 2 | 2.83 | 2.29 | 7 | 1.31 | 1.14 |
| 3 | 2.78 | 2.38 | 8 | 1.16 | 0.98 |
| 4 | 2.95 | 1.99 | 9 | 1.22 | 1.05 |
| 5 | 2.86 | 1.90 | Total | 16.69 | 13.10 |
DownLoad:
CSV
| Type of RS | Parameters | Value | Parameters | Value |
| PV | Peak power, Pr (kW) | 100 | kU (mV/℃); kI (mA/℃) | 14.4; 1.22 |
| UOC; UMPP (V) | 22; 17.32 | ISC; IMPP (A) | 5.3; 4.76 | |
| Life time (year) | 25 | θN; θA (℃) | 25; 43 | |
| Invested cost (M＄/MW) | 1.5 | O & M cost ($/MWh) | 20 | |
| WT | Rate power, Pr (kW) | 100 | vr (m/s) | 14 |
| vci (m/s) | 4 | vco (m/s) | 24 | |
| Invested cost (M＄/MW) | 1.8 | O & M cost ($/MWh) | 40 | |
| Life time (year) | 25 |
| Parameters | Value |
| Emission tax of CO2, β | 10 ＄/ton |
| Emission coefficient of traditional sources, ξT | 0.68 kg/kWh |
| Emission coefficient of RS, ξRS | 0.02 kg/kWh |
| Discount rate, r | 10% |
| State | PV | WT | Load | Electrical price | ||||
| $k_s^{RS}$ | $\lambda _s^{PV}$ | $k_s^{RS}$ | $\lambda _s^{WT}$ | $k_s^L$ | $\lambda _s^L$ | $k_s^\rho $ | $\lambda _s^\rho $ | |
| 1 | 0 | 0.450 | 0 | 0.199 | 0.30 | 0.001 | 0.30 | 0.012 |
| 2 | 0.06 | 0.002 | 0.1 | 0.112 | 0.35 | 0.003 | 0.35 | 0.023 |
| 3 | 0.12 | 0.005 | 0.2 | 0.113 | 0.40 | 0.010 | 0.40 | 0.041 |
| 4 | 0.18 | 0.012 | 0.3 | 0.106 | 0.45 | 0.028 | 0.45 | 0.065 |
| 5 | 0.23 | 0.023 | 0.4 | 0.096 | 0.50 | 0.062 | 0.50 | 0.092 |
| 6 | 0.29 | 0.039 | 0.5 | 0.085 | 0.55 | 0.111 | 0.55 | 0.116 |
| 7 | 0.35 | 0.057 | 0.6 | 0.072 | 0.60 | 0.159 | 0.60 | 0.131 |
| 8 | 0.40 | 0.073 | 0.7 | 0.059 | 0.65 | 0.184 | 0.65 | 0.132 |
| 9 | 0.45 | 0.083 | 0.8 | 0.047 | 0.70 | 0.172 | 0.70 | 0.119 |
| 10 | 0.50 | 0.084 | 0.9 | 0.036 | 0.75 | 0.130 | 0.75 | 0.096 |
| 11 | 0.55 | 0.072 | 1.0 | 0.027 | 0.80 | 0.079 | 0.80 | 0.070 |
| 12 | 0.60 | 0.053 | 1.0 | 0.020 | 0.85 | 0.039 | 0.85 | 0.045 |
| 13 | 0.65 | 0.031 | 1.0 | 0.014 | 0.90 | 0.015 | 0.90 | 0.026 |
| 14 | 0.70 | 0.013 | 1.0 | 0.009 | 0.95 | 0.005 | 0.95 | 0.014 |
| 15 | 0.75 | 0.003 | 1.0 | 0.006 | 1.00 | 0.001 | 1.00 | 0.006 |
| Type of RS | Siting (bus) | Rated power (MW) | Planning year |
| PV | 7 | 2.3 | 1 |
| 8 | 2.3 | 1 | |
| 9 | 2.4 | 1 | |
| WT | 7 | 2.7 | 6 |
| 8 | 2.6 | 6 | |
| 9 | 2.8 | 6 |
| Case | Total life cycle cost, M＄ | Emission taxes cost, M＄ | CO2 emission, ton |
| Without RS | 82.92 | 5.34 | 746,708.0 |
| With RS | 76.58 | 4.22 | 588,446.6 |
| Comparison of two cases | 6.34 | 1.12 | 158,261.4 |
| Bus i-j | $S_{ij.\max }^F$(MVA) | Lij (km) | Rij (Ω) | Xij (Ω) |
| 1-2 | 6.67 | 5.3 | 4.10 | 2.27 |
| 1-3 | 6.67 | 9.2 | 7.11 | 3.95 |
| 1-4 | 6.67 | 10.3 | 7.96 | 4.42 |
| 1-5 | 6.67 | 10.5 | 8.12 | 4.50 |
| 2-6 | 6.67 | 11.7 | 9.04 | 5.02 |
| 3-7 | 6.67 | 18.2 | 14.07 | 7.81 |
| 4-8 | 6.67 | 17.2 | 13.30 | 7.38 |
| 5-9 | 6.67 | 15.2 | 11.75 | 6.52 |
| Bus | $P_i^L$(MW) | $Q_i^L$(MVAr) | Bus | $P_i^L$(MW) | $Q_i^L$(MVAr) |
| 1 | - | - | 6 | 1.58 | 1.37 |
| 2 | 2.83 | 2.29 | 7 | 1.31 | 1.14 |
| 3 | 2.78 | 2.38 | 8 | 1.16 | 0.98 |
| 4 | 2.95 | 1.99 | 9 | 1.22 | 1.05 |
| 5 | 2.86 | 1.90 | Total | 16.69 | 13.10 |
