Research article Topical Sections

Optimization and cost-benefit analysis of a grid-connected solar photovoltaic system

  • Growing energy demand has exacerbated the issue of energy security and caused us to necessitate the utilization of renewable resources. The best alternative for promoting generation in Bangladesh from renewable energy is solar photovoltaic technology. Grid-connected solar photovoltaic (PV) systems are becoming increasingly popular, considering solar potential and the recent cost of PV modules. This study proposes a grid-connected solar PV system with a net metering strategy using the Hybrid Optimization of Multiple Electric Renewables model. The HOMER model is used to evaluate raw data, to create a demand cycle using data from load surveys, and to find the best cost-effective configuration. A sensitivity analysis was also conducted to assess the impact of differences in radiation from the solar (4, 4.59, 4.65, 5 kWh/m2/day), PV capacity (0 kW, 100 kW, 200 kW, 300 kW, 350 kW, 400 kW, 420 kW), and grid prices ($0.107, $0.118, $0.14 per kWh) upon that optimum configuration. Outcomes reveal that combining 420 kW of PV with a 405-kW converter and connecting to the utility grid is the least expensive and ecologically healthy configuration of the system. The electricity generation cost is estimated to be 0.0725 dollars per kilowatt-hour, and the net present value is 1.83 million dollars with a payback period of 6.4 years based on the system's 20-year lifespan. Also, compared to the existing grid and diesel-generator system, the optimized system, with a renewable fraction of 31.10%, provides a reduction in carbon dioxide emissions of 191 tons and 1,028 tons, respectively, each year.

    Citation: Md. Mehadi Hasan Shamim, Sidratul Montaha Silmee, Md. Mamun Sikder. Optimization and cost-benefit analysis of a grid-connected solar photovoltaic system[J]. AIMS Energy, 2022, 10(3): 434-457. doi: 10.3934/energy.2022022

    Related Papers:

    [1] Bei Liu, Wenbin Tan, Xian Zhang, Ziqi Peng, Jing Cao . Recognition study of denatured biological tissues based on multi-scale rescaled range permutation entropy. Mathematical Biosciences and Engineering, 2022, 19(1): 102-114. doi: 10.3934/mbe.2022005
    [2] David Cuesta–Frau . Permutation entropy: Influence of amplitude information on time series classification performance. Mathematical Biosciences and Engineering, 2019, 16(6): 6842-6857. doi: 10.3934/mbe.2019342
    [3] J. Murillo-Escobar, Y. E. Jaramillo-Munera, D. A. Orrego-Metaute, E. Delgado-Trejos, D. Cuesta-Frau . Muscle fatigue analysis during dynamic contractions based on biomechanical features and Permutation Entropy. Mathematical Biosciences and Engineering, 2020, 17(3): 2592-2615. doi: 10.3934/mbe.2020142
    [4] Xuyang Xie, Zichun Yang, Lei Zhang, Guoqing Zeng, Xuefeng Wang, Peng Zhang, Guobing Chen . An improved Autogram and MOMEDA method to detect weak compound fault in rolling bearings. Mathematical Biosciences and Engineering, 2022, 19(10): 10424-10444. doi: 10.3934/mbe.2022488
    [5] Enas Abdulhay, Maha Alafeef, Hikmat Hadoush, V. Venkataraman, N. Arunkumar . EMD-based analysis of complexity with dissociated EEG amplitude and frequency information: a data-driven robust tool -for Autism diagnosis- compared to multi-scale entropy approach. Mathematical Biosciences and Engineering, 2022, 19(5): 5031-5054. doi: 10.3934/mbe.2022235
    [6] Fu Tan, Bing Wang, Daijun Wei . A new structural entropy measurement of networks based on the nonextensive statistical mechanics and hub repulsion. Mathematical Biosciences and Engineering, 2021, 18(6): 9253-9263. doi: 10.3934/mbe.2021455
    [7] Shaojun Zhu, Jinhui Zhao, Yating Wu, Qingshan She . Intermuscular coupling network analysis of upper limbs based on R-vine copula transfer entropy. Mathematical Biosciences and Engineering, 2022, 19(9): 9437-9456. doi: 10.3934/mbe.2022439
    [8] Bei Liu, Hongzi Bai, Wei Chen, Huaquan Chen, Zhen Zhang . Automatic detection method of epileptic seizures based on IRCMDE and PSO-SVM. Mathematical Biosciences and Engineering, 2023, 20(5): 9349-9363. doi: 10.3934/mbe.2023410
    [9] Ziqi Peng, Seiroh Okaneya, Hongzi Bai, Chuangxing Wu, Bei Liu, Tatsuo Shiina . Proposal of dental demineralization diagnosis with OCT echo based on multiscale entropy analysis. Mathematical Biosciences and Engineering, 2024, 21(3): 4421-4439. doi: 10.3934/mbe.2024195
    [10] Erik M. Bollt, Joseph D. Skufca, Stephen J . McGregor . Control entropy: A complexity measure for nonstationary signals. Mathematical Biosciences and Engineering, 2009, 6(1): 1-25. doi: 10.3934/mbe.2009.6.1
  • Growing energy demand has exacerbated the issue of energy security and caused us to necessitate the utilization of renewable resources. The best alternative for promoting generation in Bangladesh from renewable energy is solar photovoltaic technology. Grid-connected solar photovoltaic (PV) systems are becoming increasingly popular, considering solar potential and the recent cost of PV modules. This study proposes a grid-connected solar PV system with a net metering strategy using the Hybrid Optimization of Multiple Electric Renewables model. The HOMER model is used to evaluate raw data, to create a demand cycle using data from load surveys, and to find the best cost-effective configuration. A sensitivity analysis was also conducted to assess the impact of differences in radiation from the solar (4, 4.59, 4.65, 5 kWh/m2/day), PV capacity (0 kW, 100 kW, 200 kW, 300 kW, 350 kW, 400 kW, 420 kW), and grid prices ($0.107, $0.118, $0.14 per kWh) upon that optimum configuration. Outcomes reveal that combining 420 kW of PV with a 405-kW converter and connecting to the utility grid is the least expensive and ecologically healthy configuration of the system. The electricity generation cost is estimated to be 0.0725 dollars per kilowatt-hour, and the net present value is 1.83 million dollars with a payback period of 6.4 years based on the system's 20-year lifespan. Also, compared to the existing grid and diesel-generator system, the optimized system, with a renewable fraction of 31.10%, provides a reduction in carbon dioxide emissions of 191 tons and 1,028 tons, respectively, each year.



    Regardless of its relatively young age in comparison with other entropy statistics, Permutation Entropy (PE) has soon become one of the most utilised time series entropy–related measures. It was proposed in the well known paper by Bandt and Pompe [1] in 2002, and since then, it has given rise to a number of applications and further algorithm developments. This number is growing exponentially [2], as a clear evidence of the utility of the PE approach.

    Regarding PE applications, it has been used mainly in physiological records classification. There are many scientific papers that illustrate this point. For example, we have used PE to successfully classify body temperature records from febrile and healthy patients [3] and to classify glucose records from patients at diabetes risk [4]. More frequent applications are based on electroencephalogram [5,6,7,8] and heart rate variability [9,10,11] analysis.

    Other PE fields of application are also receiving a great deal of attention. In econometrics, there are quite many interesting examples already available. In [12], the authors used PE to try to unearth the dynamical properties of time series featuring Dow Jones Industrial Average data from 1901 to 2016, obtaining a very high degree of randomness related to the market efficiency. They also used PE on temporal windows to identify market events. Along this line, [13] applied PE to stock market time series of several countries to estimate their investment attractiveness. Their findings confirmed the validity of this approach, and they also found a clear correlation between market crisis and PE value decline. The paper [14] studied the evolution of stock market efficiency during the last financial crisis. It used some specific stock exchange indices, with its main focus on the differences between pre and post crisis data processed by PE. In the same context, [15] used the complexity-entropy causality plane and a permutation statistical complexity to analyse financial time series. Using permutation Shannon entropy and permutation Fisher information measure, [16] analysed the possible deterministic components in Libor rates time series induced by manipulation. The study in [17] applied a dynamic approach to detect structural changes in time series using different entropy measures related to PE, including Gaussian, Rényi, Tsallis, and Shannon entropies. They first developed a theoretical study over synthetic time series with abrupt changes in order to assess the transition detection power. Then they applied the methods to real time series, seismic data, and economic data (exchange rates between US dollar and gold, and Nasdaq time series).

    In mechanical engineering, PE has also clearly demonstrated its usefulness, mainly in the framework of fault diagnosis. This is the case in [18]. This study uses vibration signals for bearing fault diagnosis based on multiscale PE and multinomial logit model. The classification accuracy achieved was close to 100%. The researchers in [19] also used multiscale PE for automatic recognition of weak faults in hydraulic systems from vibration signals as well. All the diagnosis methods tested achieved at least a recognition rate of 89%, and were able to discern among normal state, slight, moderate, and severe leakage.

    The basic PE algorithm has also been improved since its initial version. Two possible algorithm weaknesses were almost immediately detected by researchers regarding PE. First, PE is based on relative frequency of ordinal patterns, but it does not take into account the possible influence of amplitude differences. For example, the subsequences (0.25,1.7,0.33) and (100,99,99.5) correspond to the same ordinal pattern (1,2,0), but from an amplitude perspective they are very different. Some PE algorithm improvements have been proposed in the last years to have the amplitude patterns somehow covered in PE. For example, the Weighted Permutation Entropy method [20,21,22] applies a correction factor before computing the relative frequencies, based on the subsequence variance. The Fine Grained Permutation method [23] adds a new symbol to the ordinal motif of the subsequence that accounts for the amplitude differences. The amplitude information overlooked by the standard PE method has been demonstrated to be significant in many classification tests [24].

    The other claimed PE flaw is the impossibility to assign a single motif to subsequences that contain equal values or ties. For example, the subsequence (2,1,2) could be assigned the ordinal pattern or motif (1,0,2), but also the (1,2,0). In the seminal paper [1], Bandt and Pompe already acknowledged this possible issue, and proposed to add a small random noise to break ties in the unlikely event that equal values fall within the same subsequence. However, such ties are not that unlikely, and some improvements have also been suggested since then to minimise the possible histogram bias due to ordinal ambiguities. Thus, the modified PE method [11] generates more motifs corresponding to the possible ties. For example, the motif (0,1,2,3) generates also ordinal patterns (0,1,22), (0,11,3), (0,11,2), and (0,111). The method in [25] uses a Bayesian missing data imputation to learn from the unambiguous subsequences what are the most likely motifs in case of ties. The presence of such ties can lead to the incorrect interpretation of the chaotic nature of the time series, but they exert a minor impact on classification performance [26].

    Another line of research related to PE is the use of complex networks to map time series [27]. Specifically, the use of transition networks with nodes featuring ordinal patterns as for PE, connected by edges based on temporal succession information from a time series, is a very promising research topic for the future in the field of non–linear dynamics analysis [28]. We also have results in this regard already, using Hidden Markov Models to synthetically generate time series based on the transition probabilities between consecutive ordinal patterns [29].

    Other methods have been proposed to improve the performance and robustness of PE. Special attention deserves the Amplitude Aware Permutation Entropy method [30], a method that addresses simultaneously the two PE problems described above. This is also the case for the recently published Improved Permutation Entropy method [31], which adds an amplitude quantization stage to account for amplitude differences, ties, and noise. A different approach is the Bubble Entropy (BE) method [32], devised to reduce the dependence of PE on input parameters, such as the data length or the embedding dimension, by counting the number of sample swaps necessary to achieve the ordered subsequences instead of counting order patterns. It seems that BE exhibits more stability and discriminating power than the standard PE method [32].

    Since the BE approach could be a game–changer, based on sorting relations instead of order relations, we hypothesized that PE and BE could exhibit some kind of synergy, since they are not looking exactly at the same aspects of the time series dynamics. For example, the patterns {0,2,1} and {1,0,2} are different for PE, but obtained with the same number of swaps (1) from the original ordinal pattern {0,1,2}, therefore, they are considered the same from the BE perspective.

    The present study is based on the original PE algorithm described in [1], for a single time scale. This method computes a normalised histogram of ordinal patterns found in the subsequences drawn from a time series, when sorted in ascending order, from which the Shannon Entropy is calculated. The length of these subsequences is defined by an input parameter, the embedded dimension m.

    Formally, the input time series under analysis is defined as a vector of N components x={x0,x1,,xN1}. A generic subsequence extracted commencing at sample xj of x is defined as a vector of m components xmj={xj,xj+1,,xj+m1}. In its original state, the samples in xmj can be assigned a default growing set of indices given by πm={0,1,,m1}. The subsequence xmj undergoes then an ascending sorting process, and the sample order changes in it, are mirrored in the vector of indices πm. The resulting new version of this vector, πmj={π0,π1,,πm1}, with xj+π0xj+π1xj+π2xj+πm1, is compared, in principle, with all the possible m! ordinal patterns of length m. When a coincidence is found, a specific associated counter to that pattern, cic, is increased. This process is repeated with all the possible N(m1) subsequences (0j<Nm+1) until the complete histogram is obtained. Each bin of the histogram is finally normalised by N(m1) in order to obtain an estimation of the probability of each ordinal pattern: p={p0,p1,,pm!1}|pi=ciN(m1). This vector of probabilities is used to calculate PE as (assuming log 0=0):

    PE(x,m,N)=m!1k=0pklog pk (1)

    For example, let x={0.45,1.9,0.87,0.91,2.3,1.1,0.75,1.3,1.6,0.47,0.15,0.65,0.55,1.1,0.3} be a time series of length 15 whose PE has to be calculated using m=3. The procedure of subsequence extraction and sorting is illustrated in Table 1. Many other numerical examples of PE computation can be found in the literature. For example, in works such as [33,34].

    Table 1.  Motifs found for PE computation for the example input series x={0.45,1.9,0.87,0.91,2.3,1.1,0.75,1.3,1.6,0.47,0.15,0.65,0.55,1.1,0.3}.
    Subsequence Swap 1 Swap 2 Swap 3 Motif
    x30={0.45,1.9,0.87} {0.45,0.87,1.9} {0,2,1}
    x31={1.9,0.87,0.91} {0.87,1.9,0.91} {0.87,0.91,1.9} {0.91,0.87,1.9} {2,1,0}
    x32={0.87,0.91,2.3} {0.91,0.87,2.3} {1,0,2}
    x33={0.91,2.3,1.1} {0.91,1.1,2.3} {0,2,1}
    x34={2.3,1.1,0.75} {1.1,2.3,0.75} {1.1,0.75,2.3} {0.75,1.1,2.3} {2,1,0}
    x35={1.1,0.75,1.3} {0.75,1.1,1.3} {1,0,2}
    x36={0.75,1.3,1.6} {0.75,1.6,1.3} {1.6,0.75,1.3} {2,0,1}
    x37={1.3,1.6,0.47} {1.6,1.3,0.47} {1.6,0.47,1.3} {1,2,0}
    x38={1.6,0.47,0.15} {1.6,0.15,0.47} {0,2,1}
    x39={0.47,0.15,0.65} {0.15,0.47,0.65} {1,0,2}
    x310={0.15,0.65,0.55} {0.15,0.55,0.65} {0,2,1}
    x311={0.65,0.55,1.1} {0.55,0.65,1.1} {0.55,1.1,0.65} {1.1,0.55,0.65} {2,1,0}
    x312={0.55,1.1,0.3} {1.1,0.55,0.3} {1.1,0.3,0.55} {1,2,0}

     | Show Table
    DownLoad: CSV

    BE is a very recently proposed entropy measure [32] that has not received the attention it deserves yet, but it will surely become an indispensable tool in the field of non–linear dynamics analysis due to the possible improvements over PE it introduces.

    The main objective of BE was to minimise the dependence of entropy measures on input parameters. In general, many of the most utilised measures require at least two for their computation, an embedded dimension m, and a similarity threshold r. Depending on the values of these parameters, the performance can vary significantly in terms of discriminating power, robustness to artifacts, or any other desirable feature. Along this line, the Rank–Entropy method (RankEn) described in [35] tried to reduce the dependence on r by calculating the amount of shuffling that the distances between two subsequences under comparison had to undergo to sort such distances in ascending order. In this case, the parameter r determines the maximum rank of the set of distances that contribute to the entropy measure, and therefore, the r value is less critical.

    In [32], researchers went one step further in order to remove completely the r parameter from the entropy computation. Taking advantage of the PE method, already independent of r, and the shuffling used for RankEn, the authors proposed first a new method, called Conditional Rényi Permutation Entropy, which combined Conditional Entropy (CE) [36], and Rényi Permutation Entropy (RPE) [37]. CE can be computed as CE(m)=PE(m+1)PE(m). RPE achieved the best results when using a quadratic approach, as for Sample Entropy, and outperformed other entropy methods in terms of discriminating power [38]. Then, BE was finally defined by using two consecutive m values as CE, with a quadratic Rényi definition of entropy, and processing ordinal patterns, as PE. This way, BE does not need parameter r, is less dependent on m, or even independent for large m values, it gives more emphasis to peaks without neglecting lower values, the number of unique possible states is reduced (ties and amplitude influence are less critical since more matches to compute relative frequencies can be found), and converges faster than PE, with a higher discriminating power [32].

    The core of the BE algorithm is that of PE, but instead of computing the relative frequency of the ordinal patterns, it computes the relative frequency of the necessary swaps to get an ordered subsequence. First, each subsequence xmj is also sorted in ascending order using a bubble sort algorithm. The counter vector c stores instead all the swaps necessary in each case, with a maximum given by [0,m(m1)2]. Each bin is normalised by Nm+1 too. From all the resulting relative frequencies pi, accounting for how likely a number of swaps is [32], the Rényi entropy of order 2 is computed as:

    Hm2(x)=logm(m1)2k=0p2k (2)

    The embedding dimension is then increased by 1, mm+1, and the procedure is repeated again to obtain a new entropy value from Eq. 2, Hm+12. Finally, BE is obtained in a similar way as for Approximate Entropy, ApEn [39]:

    BE(x,m,N)=(Hm+12Hm2)logm+1m1 (3)

    The procedure of subsequence extraction, sorting, and swap computation, for the same input time series as for PE, is illustrated in Table 2. The maximum number of swaps in this case is 3.

    Table 2.  Necessary swaps to sort in ascending order all the subsequences that can be extracted from the example data to compute BE.
    Subsequence Swap 1 Swap 2 Swap 3 Swaps
    x30={0.45,1.9,0.87} {0.45,0.87,1.9} 1
    x31={1.9,0.87,0.91} {0.87,1.9,0.91} {0.87,0.91,1.9} {0.91,0.87,1.9} 3
    x32={0.87,0.91,2.3} {0.91,0.87,2.3} 1
    x33={0.91,2.3,1.1} {0.91,1.1,2.3} 1
    x34={2.3,1.1,0.75} {1.1,2.3,0.75} {1.1,0.75,2.3} {0.75,1.1,2.3} 3
    x35={1.1,0.75,1.3} {0.75,1.1,1.3} 1
    x36={0.75,1.3,1.6} {0.75,1.6,1.3} {1.6,0.75,1.3} 2
    x37={1.3,1.6,0.47} {1.6,1.3,0.47} {1.6,0.47,1.3} 2
    x38={1.6,0.47,0.15} {1.6,0.15,0.47} 1
    x39={0.47,0.15,0.65} {0.15,0.47,0.65} 1
    x310={0.15,0.65,0.55} {0.15,0.55,0.65} 1
    x311={0.65,0.55,1.1} {0.55,0.65,1.1} {0.55,1.1,0.65} {1.1,0.55,0.65} 3
    x312={0.55,1.1,0.3} {1.1,0.55,0.3} {1.1,0.3,0.55} 2

     | Show Table
    DownLoad: CSV

    From the results in Table 2, the counter vector components for the first iteration of the method are obtained as: c0=0, c1=7, c2=3, and c3=3. These values, once normalised, will become the probability vector p from which entropy Hm can be computed. The process is then repeated for mm+1=4 to obtain Hm+1, with results: c0=0, c1=1, c2=2, c3=2, c4=6, c5=1, and c6=0. Finally, applying Eq. 3, the value for BE obtained is 0.3115, with H4=1.1411, and H3=0.9252. An implementation of this method is detailed in Algorithm 1.

    Algorithm 1 Bubble Entropy (BE) Algorithm
    Input: x, m>2, N>m+1
    Initialisation: BE 0, c{0}, p{0}
    for Mm,,m+1 do
      for j0,,NM do
        xMj{xj,xj+1,,xj+M1}
        bSorted false
        nSwaps 0
        while (bSorted=false) do      ▷Bubble sort
          bSorted true
          for ij,,j+M2 do
            if (xi>xi+1) then
              swap(xi,xi+1)
              nSwaps=nSwaps+1
              bSorted false
            end if
          end for
        end while
        cnSwapscnSwaps+1       ▷Update counter
    end for
    HM20
    for k0,,m(m1)2 do
        pkckNM+1
        HM2HM2(p2k)
     end for
      HM2logHM2
    end for
    BE=Hm+12Hm2
    BE=BE/logm+1m1
    Output: BE(x,m,N)

    The performance analysis of the methods under analysis will be assessed both in quantitative and qualitative terms. The classification accuracy will be quantified as the ratio of correctly assigned time series to classes. For example, a 0.80 classification accuracy will refer to the fact that 80% of the time series were labelled as their true class after the classification analysis has been completed. The remaining 20% has been incorrectly included in another class distinct to the true class to which they belong. The classification accuracy was quantified using the classical sensitivity, specificity and accuracy performance indicators. They were obtained using a threshold from the closest point in the ROC curve [40] to the (0, 1) point [41]. We first proceeded by computing the number of True Positives (TP), False Negatives (FN), True Negatives (TN), and False Positives (FP). Then, the classical performance metrics Sensitivity = TP / (TP+FN), Specificity = TN / (TN+FP), and Accuracy = (TN + TP) / (TN + TP + FN + FP) [42], were obtained.

    However, the numerical result of the classification accuracy does not suffice to provide a complete picture of the discriminating power of the method. For example, a relatively high accuracy of 0.75 can be achieved for the classification of two size–balanced classes, but with a sensitivity close to 0.5, pure guess, and a specificity close to 1, a useless performance since there must be a trade–off between both parameters to avoid statistical uncertainty.

    Therefore, in order to assess the statistical significance of the results quantified in terms of classification accuracy, we used a bootstrap version [43] of the equality of means (there is no difference between the classes under analysis according to the values of PE or BE obtained for classification) hypothesis testing. By using this bootstrap method, it is not necessary to assume any specific distribution of the data, and it can be applied even to small sizes of data.

    Given two samples Y={Y1,,YnY} and Z={Z1,,ZnZ}, the steps to carry out this test are:

    1. Calculate the T statistic as in the standard t–Student test [44].

    2. Resample the input data randomly with replacement to create bootstrapped versions of the data, Y and Z.

    3. Calculate a bootstrap statistic T as:

    T=(ˆμYˆμZ)(ˆμYˆμZ)ˆσ2YnY1+ˆσ2ZnZ1,

    where μ is the mean, and σ2 the variance of the data.

    4. Repeat steps 2 and 3 several times (k=200 in this case) to obtain the corresponding statistics T1,T2,,Tk.

    5. Sort the previous statistics in increasing order, T(1),T(2),,T(k).

    6. Reject the equal means hypothesis, H:μY=μZ, if T<T(q1) or T>T(q2), where q1=kα/2 and q2=kq1+1.

    The null hypothesis of equal means of the classes in terms of PE or BE was assessed using this test. The acceptance threshold was set at α=0.05.

    The possible synergy between PE and BE can be exploited in many ways. They can become independent variables of the same classification function, as in logistic regression applications [45,3], among many other similar methods. The input dataset can be split into training and test sets to obtain more complex polynomial classification functions, as is the case when using neural networks [46]. They can also conform the input feature vector for an unsupervised classification approach featured by a clustering algorithm [47].

    We have chosen the clustering approach for its simplicity and good results in previous similar studies [48]. The specific method employed is the Max–Min algorithm [49]. The input number of clusters was set to 2, and therefore 2 centroids are selected, one from each experimental dataset. In order to further simplify the clustering method, since the goal was not to design a classification method but rather explore the possible synergy between PE and BE, the centroids are chosen from those time series with the maximum and minimum PE value, the furthest points in the PE feature space. Finally, each time series is assigned to the class of the closest centroid in a single iteration, quantified using a Euclidean distance. If the goal was to maximise the classification performance, a more evolved clustering algorithm could be chosen, such as a genetic method that converges to a global minimum cost [50], or a density–based method, that enables the use of highly complex partition regions [51]. Other distance metrics, or more iterations with centroid updating could also be used, such as in the k–Means method or its variants [52,53].

    The experimental dataset was composed of 7 biomedical records previously used in other studies. In addition, 4 out of the 7 datasets are publicly available at Physionet [54]. The information regarding the specific datasets employed is summarised in Table 3, for our proprietary data, and in Table 4, for the publicly available data at Physionet. Each row includes the dataset assigned name, a short description, and references to further studies or information about the specific dataset in each case. An example plot of one record is also included.

    Table 3.  Summary of the own experimental datasets employed in the present study.
    Name Description Plot
    TEMPERATURE 8h temperature records of 16 healthy patients and 14 patients that developed a fever 24h before the recording. Sampling rate of 1 sample per minute [55,3].
    SURVIVAL Several days of temperature records of critically ill patients. A total of 36 records, 18 for survivors, and 18 for non-survivors. Sampling rate of 10 minutes [56].
    GLUCOSE 24h records of 206 patients at risk of developing a diabetes. After 3 years, 18 of these patients became diabetic. Sampling rate of 5 minutes [4,57].

     | Show Table
    DownLoad: CSV
    Table 4.  Summary of the publicly available experimental datasets employed in the present study.
    Name Description Plot
    BONN 4097 samples long electroencephalograph records. Two classes, with and without epileptic seizures, 100 records each. Sampling rate of 173.61 Hz [58,59].
    EMG Three classes, healthy (10 records), myopathy (22 records), and neuropathy (29 records) electromyogram records of 5000 samples. Sampling rate of 4KHz [54,59].
    FANT Two classes of RR--interval records from 20 young and 20 elderly healthy subjects monitored during 120 minutes. The length of the records varies, but they are longer than 5000 samples in all cases [54,24,60].
    PAF Two classes of RR--interval records, 25 paroxysmal atrial fibrillation episode free, and 25 with. The duration of the records is 5 minutes, with lengths mostly around 400 samples [54,59].

     | Show Table
    DownLoad: CSV

    The experiments were devised to assess the performance of BE and PE using a different and more varied dataset than in the seminal paper of BE [32], and also to explore the possible synergies between both measures. These experiments also included a variation of the parameter m, from 3 to 8, since input parameter influence is another topic of intense debate and research in the scientific literature. All the experiments and algorithms were implemented in C++ programming language, using MinGW 4.9.2 32 bit compiler (www.mingw.org).

    In this regard, Table 5 shows the classification performance (if statistically significant) of the results achieved by PE and BE using the TEMPERATURE dataset. In this case, PE exhibits a high classification accuracy, around 0.8 in some m cases. BE performance is not significant for any m value.

    Table 5.  Classification results achieved using PE and BE, applied to the TEMPERATURE database.
    m PE BE
    Sensitivity Specificity Accuracy Sensitivity Specificity Accuracy
    3
    4
    5
    6 0.81 0.57 0.70
    7 0.93 0.64 0.80
    8 0.81 0.85 0.83

     | Show Table
    DownLoad: CSV

    For the SURVIVAL dataset, the classification results achieved are shown in Table 6. The results for BE are not statistically significant for m=8, being PE more stable with m. The maximum performance achieved by both methods was 0.77.

    Table 6.  Classification results achieved using PE and BE, applied to the SURVIVAL database.
    m PE BE
    Sensitivity Specificity Accuracy Sensitivity Specificity Accuracy
    3 0.72 0.72 0.72 0.72 0.77 0.75
    4 0.72 0.77 0.75 0.72 0.77 0.75
    5 0.72 0.77 0.75 0.77 0.77 0.77
    6 0.72 0.77 0.75 0.72 0.77 0.75
    7 0.72 0.83 0.77 0.61 0.61 0.61
    8 0.72 0.83 0.77

     | Show Table
    DownLoad: CSV

    For the last dataset in Table 3, GLUCOSE, the results are listed in Table 7. Both statistics yielded significant classification results in all m cases tested, but the performance of BE was slightly superior to that of PE, 0.80 vs. 0.79. This small difference is more significant if sensitivity is taken into account, since with PE it was 0.61, but with BE it was 0.66, a clear difference for a small class size of 18 records.

    Table 7.  Classification results achieved using PE and BE, applied to the GLUCOSE database.
    m PE BE
    Sensitivity Specificity Accuracy Sensitivity Specificity Accuracy
    3 0.77 0.56 0.58 0.66 0.71 0.70
    4 0.72 0.67 0.67 0.72 0.69 0.69
    5 0.61 0.80 0.79 0.72 0.77 0.77
    6 0.77 0.63 0.65 0.77 0.66 0.67
    7 0.72 0.73 0.73 0.66 0.74 0.73
    8 0.72 0.72 0.72 0.66 0.81 0.80

     | Show Table
    DownLoad: CSV

    Table 8 shows the classification performance achieved by PE and BE using the BONN dataset. In this case, PE exhibits a very high classification accuracy, close to 0.9 in all m cases. BE performance is lower, 0.785 at most, and even with no significant results for m=5,6,7.

    Table 8.  Classification results achieved using PE and BE, applied to the BONN database.
    m PE BE
    Sensitivity Specificity Accuracy Sensitivity Specificity Accuracy
    3 0.93 0.90 0.91 0.74 0.83 0.78
    4 0.93 0.89 0.91 0.70 0.67 0.68
    5 0.92 0.89 0.90
    6 0.91 0.89 0.90
    7 0.93 0.85 0.89
    8 0.90 0.85 0.87 0.75 0.52 0.63

     | Show Table
    DownLoad: CSV

    In Table 9, the classification results for the three classes in the EMG database are shown (healthy, myopathy, neuropathy). PE was not able to achieve significant results for the second case (Healthy vs. Myopathy), whereas BE found differences in all cases for m=3 and 5. Therefore, BE can be considered to outperform PE in this experiment.

    Table 9.  Classification results achieved using PE and BE, applied to the EMG database. There are three rows for each case, first row for results comparing healthy and myopathy records, second row for healthy and neuropathy, and the third row for myopathy and neuropathy comparison.
    m PE BE
    Sensitivity Specificity Accuracy Sensitivity Specificity Accuracy
    3 1 1 1 1 1 1
    0.9 0.58 0.66
    0.55 0.77 0.67 0.96 1 0.98
    4 1 1 1 1 1 1
    1 1 1 1 1 1
    5 1 1 1 0.86 0.90 0.88
    0.75 0.80 0.78
    1 1 1 0.51 0.77 0.66
    6 1 1 1 0.70 0.77 0.75
    1 1 1 0.75 0.77 0.76
    7 1 1 1 1 1 1
    1 1 1 1 1 1
    8 1 1 1 1 1 1
    1 1 1 1 1 1

     | Show Table
    DownLoad: CSV

    Table 10 shows the results for the FANT database. In this case, BE was clearly better than PE, with significant results for m=7,8, and a maximum classification accuracy of 0.875. PE was not capable of achieving the significance threshold, probably because in many experiments either Sensitivity or Specificity were close to 0.5, despite a good accuracy in the other measure.

    Table 10.  Classification results achieved using PE and BE, applied to the FANT database.
    m PE BE
    Sensitivity Specificity Accuracy Sensitivity Specificity Accuracy
    3
    4
    5
    6
    7 0.75 1 0.875
    8 0.80 0.85 0.82

     | Show Table
    DownLoad: CSV

    Finally, the results of the last dataset are shown in Table 11. The classification accuracy of PE was quite stable around 0.8, but BE clearly failed in this case.

    Table 11.  Classification results achieved using PE and BE, applied to the PAF database.
    m PE BE
    Sensitivity Specificity Accuracy Sensitivity Specificity Accuracy
    3 0.76 0.88 0.82
    4 0.80 0.84 0.82
    5 0.80 0.80 0.80
    6 0.92 0.72 0.82
    7 0.96 0.68 0.82
    8 0.92 0.68 0.80

     | Show Table
    DownLoad: CSV

    A summary of the best statistically significant accuracy results achieved using the two measures assessed (PE, BE) is shown in Table 12. Of the 7 datasets, PE yielded a better accuracy in 3 cases, and BE in other 3. In 2 datasets each, PE and BE were not capable of finding significant differences. There was no dataset where none of the measures studied failed to find significant differences, but in 4 of the experiments, one of the measures achieved a good performance, whereas the other was unable to achieve statistically significant results. It is not possible to say what measure was better in this case. However, what became more apparent was the fact that there exists a clear complementarity between PE and BE, each one seems to look at different properties of the time series dynamics, and this should be exploited.

    Table 12.  Summary of the performance achieved by PE and BE. Only the best significant case is reported in terms of overall accuracy.
    Dataset PE BE
    TEMPERATURE 0.80
    SURVIVAL 0.77 0.77
    GLUCOSE 0.79 0.80
    BONN 0.91 0.78
    EMG 1, 0.66, 0.98
    FANT 0.87
    PAF 0.82

     | Show Table
    DownLoad: CSV

    The second group of experiments were focused on exploring the possible synergy between PE and BE, used jointly for classification. Without any specific customisation, the clustering algorithm described in Section 2.3 was applied to some datasets in a completely unsupervised way. The quantitative results achieved are shown in Tables 13, 14 and 15 to provide an example of beneficial synergy or no synergy.

    Table 13.  Clustering results using the GLUCOSE database and features PE and BE.
    m Clustering PE BE
    3 0.91 0.58 0.70
    4 0.91 0.67 0.69
    5 0.91 0.79 0.77
    6 0.91 0.65 0.67
    7 0.91 0.73 0.73
    8 0.91 0.72 0.80

     | Show Table
    DownLoad: CSV
    Table 14.  Clustering results using the TEMPERATURE database and features PE and BE.
    m Clustering PE BE
    3 0.53
    4 0.60
    5 0.53
    6 0.53 0.70
    7 0.56 0.80
    8 0.70 0.83

     | Show Table
    DownLoad: CSV
    Table 15.  Clustering results using the SURVIVAL database and features PE and BE.
    m Clustering PE BE
    3 0.69 0.72 0.75
    4 0.75 0.75 0.75
    5 0.72 0.75 0.77
    6 0.58 0.75 0.75
    7 0.58 0.77 0.61
    8 0.50 0.77

     | Show Table
    DownLoad: CSV

    The results for the GLUCOSE dataset in Table 13 probably best illustrate the optimal synergy between BE and PE. The classification performance achieved combining these two measures is far higher than that achieved individually, which confirms the hypothesis that there exits at least some PE–BE complementarity. Moreover, it is independent of the m value, an almost ideal case.

    However, one size does not fit all, and the method should be customised for each dataset. For example, what worked in Table 13, for the GLUCOSE dataset, did not work in other cases. Table 14 shows the classification accuracy achieved using the TEMPERATURE dataset, where the individual performance using PE is higher than using both measures. Since BE did not achieve any significant result, it seems it has a blurring effect on the differences between classes. The discerning method should be customised to automatically detect the measure that better features the differences, and therefore increase not only the accuracy, but also the robustness of the classification by reducing the number of datasets where no significance is achieved.

    The results in Table 15 also correspond to a case where the possible synergy is not seamlessly exploited. The clustering method achieves a lower performance than the methods individually, but with a small difference for m<6. As in the previous case, with clustering method optimisation and customisation, or using other more sophisticated pattern recognition methods, it can be hypothesised that there is room for performance improvement.

    The experiments in Section 3 were first aimed at confirming the performance of the BE approach from a broader point of view than in [32], using a set of records with markedly different properties and behaviour: slow varying signals such as TEMPERATURE, SURVIVAL, and GLUCOSE datasets, faster varying signals such as the BONN dataset, spiky records in the EMG dataset, and a combination of probably all the previous, in the RR datasets.

    According to the results, the performance of BE was, in general, similar but complementary to that of PE, in terms of classification accuracy. Namely, both statistics achieved similar performance, but on different datasets. Taking into account all the m values tested, PE seems to be more robust and stable with slow varying signals, and BE for spiky records (as theoretically stated in Section 2.2), with the rest of records in between. The dependence of both PE and BE on parameter m remains an unsolved problem, since there was a great accuracy variability depending on the specific m value chosen in many cases.

    PE yielded a higher accuracy than BE for datasets TEMPERATURE (Table 5), BONN (Table 8), and PAF (Table 11), whereas BE outperformed PE when applied to GLUCOSE (Table 7), EMG (Table 9), and FANT (Table 10) time series. The performance was the same for the SURVIVAL dataset (Table 6). Obviously, a generalisation can not be made however representative the experimental dataset is, but it can arguably be stated that the combination of PE and BE is more likely to provide significant results than each measure isolatedly, mainly based on the fact that non–significance on one measure is often accompanied by high performance on the other (Table 12).

    Along this line, we explored a potential synergy using a very simple approach based on a clustering algorithm, and PE and BE as the two extracted features. If the object distribution in the PE–BE plane matches the centroid computation and the clustering algorithm geometric properties, the possible synergy can be exploited successfully, as demonstrated in Table 13. However, this is not always the case, and the combination of PE and BE can also exert a blurring effect and have a detrimental impact on classification performance (Table 14 and Table 15). Therefore, whatever the discriminating function chosen is (linear, polynomial, logistic, etc.), or the classification algorithm (clustering, neural networks, random forest, etc.), it should be customised in each case to make the most of the information provided by PE and BE. Obviously, there will be cases where no synergy will be found, whatever the method is, but in some cases as in Table 13, combining more than a single method can make a significant difference.

    Fig. 1 depicts an example of the interaction between BE and PE for the dataset GLUCOSE, with m=8. With a suitable computed class segmentation function, it is possible to increase the accuracy achieved by each method individually. As stated above, this is not always the case, and the possible synergy has to be assessed on a case–by–case basis. It will also depend on the complexity of the segmentation function employed, and on the match between the shape of the clusters and the clustering algorithm applied: spherical, linear, or other non–linear shapes. In a more general way, when one of the methods fails, either PE or BE, it is more likely the other yields a significantly better performance, since complementarity seems to be even more frequent than synergy (Tables 5, 8, 10 11).

    Figure 1.  Example of BE–PE plane for GLUCOSE dataset, with m=8. Big black squares represent patients that became diabetic. Small red squares correspond to patients that did not develop a diabetes during the study.

    Every year quite many new tools to quantify the dynamical features of time series are described in the scientific literature. These new tools are claimed to be more efficient in algorithmic terms, more sensitive, more robust, or less dependent on input parameters, among many other possible benefits. In this regard, BE is a recently proposed measure that exploits the effort to sort subsequences instead of the ordinal patterns obtained [32].

    BE was presented as an improvement over PE, less influenced by the time series length, the specific value of the input parameter m, and with a better discriminating power [32]. The present study assessed this discriminant power using several different datasets, and we would conclude that the BE discriminating capability was not clearly better, but complementary, BE succeeded where PE failed, and the other way round. This led us to try to combine both measures in a single method to take advantage of BE and PE strengths and simultaneously minimise their weaknesses. The scheme used was based on a clustering algorithm, equivalent to a linear discriminant function.

    According to the results obtained, the combination of methods can be considered a new strategy worth examining in the field of time series classification. It should not be claimed to be a cure–all method, but the classification performance confirmed the hypothesis in some cases, and it seems to be able to exploit the synergy between PE and BE provided there is some customisation to the problem at hand.

    Moreover, classification performance should not be the only item to be maximised in quantitative terms. Since there are many datasets where only one measure is able to find significant differences, other customisations should be devised to increase the robustness of the classification. In other words, new algorithms able to automatically exploit the differences provided by a measure, and minimise the influence of the confounding measure, should be proposed. This way, it would be possible to have a method able to work with a wider input range of time series or properties, with, for example, statistically significant results in more datasets used in the present study.

    There are many clustering algorithms described in the scientific literature, most of them more robust than the example used in this paper. In addition to implementation issues to save memory requirements or computational cost, the algorithms could be improved in terms of accuracy, unknown number of classes, optimality or convergence.

    This approach will need further studies using other databases and other non–linear features. The direct combination of ordinal patterns and sorting relations in a single statistic should also be investigated, or the introduction of the BE information into PE computation as a some kind of histogram weighting. The influence of time series length [59], equal values in the subsequences [26], and time delay τ, could also be characterised. Further integration with other PE improvements could be worth exploring [30,23,20,4,11].

    No funding was received to support this research work.

    The authors declare that they have no conflict of interest.



    [1] Power system master plan-2016 (2021) Available from: https://powerdivision.gov.bd/site/page/f68eb32d-cc0b-483e-b047-13eb81da6820/Power-System-Master-Plan-2016.
    [2] Nallapaneni MK, Chopra S, Chand A, et al. (2020) Hybrid renewable energy microgrid for a residential community: A techno-economic and environmental perspective in the context of the SDG7. Sustainability 12: 3944. https://doi.org/10.3390/su12103944 doi: 10.3390/su12103944
    [3] Rofiqul Islam, Rabiul Islam, Beg RA (2008) Renewable energy resources and technologies practice in Bangladesh. Renewable Sustainable Energy Rev 12: 299-343. Available from: https://ideas.repec.org/a/eee/rensus/v12y2008i2p299-343.html.
    [4] Agyekum E, Mehmood U, Kamel S, et al. (2022) Technical performance prediction and employment potential of solar PV systems in cold countries. Sustainability 14: 3546. https://doi.org/10.3390/su14063546 doi: 10.3390/su14063546
    [5] Falih H, Hamed AJ, Khalifa AHN (2022) Techno-economic assessment of a hybrid connected PV solar system. Int J Air-Cond Ref 3: 30. https://doi.org/10.1007/s44189-022-00003-7 doi: 10.1007/s44189-022-00003-7
    [6] Rafique MM, Rehman S (2017) National energy scenario of Pakistan: current status, future alternatives, and institutional infrastructure: An overview. Renewable Sustainable Energy Rev 69: 156-167. https://doi.org/10.1016/j.rser.2016.11.057 doi: 10.1016/j.rser.2016.11.057
    [7] Kumar NM, Elavarasan RM, Shafiullah GM, et al. (2020) Hybrid renewable energy microgrid for a residential community: A techno-economic and environmental perspective in the context of the SDG7. Sustainability 12: 3944. https://doi.org/10.3390/su12103944 doi: 10.3390/su12103944
    [8] Poudyal R, Loskot P, Parajuli R (2021) Techno-economic feasibility analysis of a 3-kW PV system installation in Nepal. Renewables 5: 8. https://doi.org/10.1186/s40807-021-00068-9 doi: 10.1186/s40807-021-00068-9
    [9] Adetokun BB, Ojo JO Muriithi CM (2021) Application of large‑scale grid‑connected solar photovoltaic system for voltage stability improvement of weak national grids. Sci Rep 11: 24526. https://doi.org/10.1038/s41598-021-04300-w doi: 10.1038/s41598-021-04300-w
    [10] Uddin MM, Faysal A, Raihan, et al. (2018) Present energy scenario, necessity and future prospect of renewable energy in Bangladesh. American J Eng Res 7: 45-51. Available from: https://www.researchgate.net/publication/327103112_Present_Energy_Scenario_Necessity_and_Future_Prospect_of_Renewable_Energy_in_Bangladesh.
    [11] Ministry of Power, Energy and Mineral Resources (2021) Available from: https://mpemr.gov.bd/.
    [12] Agyekum EB (2021) Techno-economic comparative analysis of solar photovoltaic power systems with and without storage systems in three different climatic regions, Ghana. Sustainable Energy Technol Assess 43: 100906. https://doi.org/10.1016/j.seta.2020.100906 doi: 10.1016/j.seta.2020.100906
    [13] PV magazine-Photovoltaics Markets and Technology (2021) Available from: https://www.pv-magazine.com/2020/04/06/world-now-has-583-5-gw-of-operational-pv/.
    [14] Uddin MN, Rahman MA, Mofijur M, et al. (2019) Renewable energy in Bangladesh: Status and prospects. Energy Procedia 160: 655-661. https://doi.org/10.1016/j.egypro.2019.02.218 doi: 10.1016/j.egypro.2019.02.218
    [15] Sustainable and renewable energy development authority (2021) Available from: http://www.sreda.gov.bd/.
    [16] Country profile data by world bank (2021) Available from: https://data.worldbank.org/country/Bangladesh?view=chart.
    [17] Central Intelligence Agency (CIA) (2022) South Asia: Bangladesh-The World Factbook. Available from: https://www.cia.gov/the-world-factbook/countries/bangladesh/.
    [18] Md. Belayet Hossain (2019-2020 Annual Report) Bangladesh Power Development Board. Available from: https://bdcom.bpdb.gov.bd/bpdb_new/resourcefile/annualreports/annualreport_1605772936_AnnualReport2019-20.pdf.
    [19] An J, Mikhaylov A, Jung SU (2020) The strategy of South Korea in the global oil market. Energies 24: 2491. https://doi.org/10.3390/en13102491 doi: 10.3390/en13102491
    [20] Mutalimov, Verdi, Kovaleva, et al. (2021) Assessing regional growth of small business in Russia. Entre Business Econos Re 9: 119-133. https://doi.org/10.15678/EBER.2021.090308 doi: 10.15678/EBER.2021.090308
    [21] Noor Abir MA, Rahman MT (2018) Energy scenario of Bangladesh & future challenges. IJSER 9: 1001-1005. Available from: https://www.ijser.org/researchpaper/Energy-Scenario-of-Bangladesh-Future-Challenges.pdf.
    [22] GM Shafiullah, Tjedza Masola, Remember Samu, et al. (2021) Prospects of hybrid renewable energy-based power system: A case study, post analysis of Chipendeke Micro-Hydro, Zimbabwe. IEEE Access 9: 73433-73452. https://doi.org/10.1109/ACCESS.2021.3078713 doi: 10.1109/ACCESS.2021.3078713
    [23] Hossain CA, Chowdhury N, Michela Lo, et al. (2019) System and cost analysis of stand-alone solar home systems applied to a developing country. Sustainability 11: 1403. Available from: https://ideas.repec.org/a/gam/jsusta/v11y2019i5p1403-d211588.html.
    [24] Global Solar Atlas (2021) Available from: https://globalsolaratlas.info/download/Bangladesh.
    [25] Photovoltaic power potential of Bangladesh (2021) Available from: https://datacatalog.worldbank.org/dataset/bangladesh-solar-irradiation-and-pv-power-potential-map.
    [26] Hossain S, Rahman M (2021) Solar energy prospects in Bangladesh: Target and current status. Energy Power Eng. 13: 322-332. https://doi.org/10.4236/epe.2021.138022 doi: 10.4236/epe.2021.138022
    [27] PV magazine-Photovoltaics markets and technology (2021) Available from: https://www.pv-magazine.com/2021/04/30/Bangladesh's-largest-pv-plant-comes-online/.
    [28] HOMER Pro-Microgrid Software for Designing Optimized Hybrid Microgrids (2021) Available from: https://www.homerenergy.com/products/pro/index.html.
    [29] Das, Nipu, Chakrabartty, et al. (2020) Present energy scenario and future energy mix of Bangladesh. Energy St Rev 32: 1-11. https://doi.org/10.1016/j.esr.2020.100576 doi: 10.1016/j.esr.2020.100576
    [30] Shuvho B, Chowdhury M, Ahmed S, et al. (2019) Prediction of solar irradiation and performance evaluation of grid-connected solar 80KWp PV plant in Bangladesh. Energy Reports 5: 714-722. https://doi.org/10.1016/j.egyr.2019.06.011 doi: 10.1016/j.egyr.2019.06.011
    [31] Kumar Das B, Hasan M, Fazlur Rashid (2021) Optimal sizing of a grid-independent PV/diesel/pump-hydro hybrid system: A case study in Bangladesh. Sustainable Energy Technol Asses 44: 100997. http://dx.doi.org/10.1016/j.seta.2021.100997 doi: 10.1016/j.seta.2021.100997
    [32] Islam A, Shima FA, Khanam A (2013) Analysis of grid-connected solar PV systems in the southeastern part of Bangladesh. Appl Sol Energy 49: 116-123. https://doi.org/10.3103/S0003701X13020035 doi: 10.3103/S0003701X13020035
    [33] Mohammad Shuhrawardy, Kazi Tanvir Ahmmed (2014) The feasibility study of a grid-connected PV system to meet the power demand in Bangladesh-A case study. AJEE 2: 59-64. Available from: https://www.sciencepublishinggroup.com/journal/paperinfo.aspx?journalid=168&doi=10.11648/j.ajee.20140202.12.
    [34] Mahmud, Nasif (2013) Modeling and economic analysis of grid-connected solar photovoltaic systems in Bangladesh. IJAET 6: 1452-1463. Available from: https://www.researchgate.net/profile/nasif-mahmud/publication/312218967_modeling_and_economic_analysis_of_grid_connected_solar_photo_voltaic_system_in_bangladesh/links/594b18cc458515225a831d7c/modeling-and-economic-analysis-of-grid-connected-solar-photo-voltaic-system-in-bangladesh.pdf.
    [35] Nurunnabi Md, Roy N (2015) Grid connected hybrid power system design using HOMER in 3rd ICAEE. https://doi.org/10.1109/ICAEE.2015.7506786
    [36] Iqra P, Abdul Razaque S, Sukru D (2017) Designing off-grid and on-grid renewable energy systems using HOMER Pro software. J Int Envirn Appl Sci 12: 270-276. Available from: https://www.acarindex.com/pdfs/433383.
    [37] Soumya Mandal, Das Barun K, Najmul Hoque (2018) Optimum sizing of a stand-alone hybrid energy system for rural electrification in Bangladesh. J Cleaner Pro 200: 12-27. https://doi.org/10.1016/j.jclepro.2018.07.257 doi: 10.1016/j.jclepro.2018.07.257
    [38] Nandi SK, Ghosh HR (2010) Prospect of wind-PV-battery hybrid power system as an alternative to grid extension in Bangladesh. Energy 35: 3040-3047. https://doi.org/10.1016/j.energy.2010.03.044 doi: 10.1016/j.energy.2010.03.044
    [39] Lipu MSH, Uddin MS, Miah MAR (2013) A feasibility study of solar-wind-diesel hybrid system in rural and remote areas of Bangladesh. IJRER 3: 4. https://doi.org/10.20508/ijrer.v3i4.898.g6220 doi: 10.20508/ijrer.v3i4.898.g6220
    [40] Pradhan SR, Amit Kumar, Ashutosh Sahoo, et al. (2017) Optimization of grid-connected hybrid energy (solar and biomass) system using HOMER Pro software. Int J Inno Sci Res Tech 2: 50-57. Available from: https://ijisrt.com/wp-content/uploads/2017/04/Optimization-of-Grid-connected-Hybrid-Energy-solar.pdf.
    [41] Ali F, Muahnmmad A, Jiang, et al. (2021) A techno-economic assessment of hybrid energy systems in rural Pakistan. Energy 215: 119103. https://doi.org/10.1016/j.energy.2020.119103 doi: 10.1016/j.energy.2020.119103
    [42] El-Tous Y (2012) A study of a grid-connected PV household system in Amman and the effect of the incentive tariff on the economic feasibility. Int J Appl Sci Tech 2: 100-105. Available from: http://www.ijastnet.com/update/journals/Vol_2_No_2_February_2012/14.pdf.
    [43] Li D, Cheung KL, Lam T, et al. (2012) A study of grid-connected photovoltaic (PV) system in Hong Kong. Appl Energy 90: 122-127. https://doi.org/10.1016/j.apenergy.2011.01.054 doi: 10.1016/j.apenergy.2011.01.054
    [44] Dawoud S, Lin XN, Sun JW, et al. (2015) Feasibility study of isolated PV-Wind hybrid system in Egypt. Ad Materials Res 1092-1093: 145-151. https://doi.org/10.4028/www.scientific.net/AMR.1092-1093.145 doi: 10.4028/www.scientific.net/AMR.1092-1093.145
    [45] Chouki Ghenai, Maamar Bettayeb (2019) Grid-tied solar PV/Fuel cell hybrid power system for university building. Energy Procedia 159: 96-103. https://doi.org/10.1016/j.egypro.2018.12.025 doi: 10.1016/j.egypro.2018.12.025
    [46] Shafiullah GM (2016) Hybrid renewable energy integration (HREI) system for subtropical climate in central Queensland, Australia. Renewable Energy 96: 1034-1053. https://doi.org/10.1016/j.renene.2016.04.101 doi: 10.1016/j.renene.2016.04.101
    [47] Gm Shafiullah, Sanjeev Kumar P, Kumar NM, et al. (2020) A comprehensive review on renewable energy development, challenges, and policies of leading Indian states with an international perspective. IEEE Access 8: 74432-74457. http://dx.doi.org/10.1109/ACCESS.2020.2988011 doi: 10.1109/ACCESS.2020.2988011
    [48] Ali I, Shafiullah GM, Urmee T, et al. (2018) A preliminary feasibility of roof-mounted solar PV systems in the Maldives. Renewable Sustainable Energy Rev 83: 18-32. Available from: https://EconPapers.repec.org/RePEc:eee:rensus:v:83:y:2018:i:c:p:18-32.
    [49] Shaahid SM, El-Amin I (2009) Techno-economic evaluation of off-grid hybrid photovoltaic-diesel-battery power systems for rural electrification in Saudi Arabia-A way forward for sustainable development. Renewable Sustainable Energy Rev 13: 625-633. https://doi.org/10.1016/j.rser.2007.11.017 doi: 10.1016/j.rser.2007.11.017
    [50] Al-Ghussain L, Samu R, Taylan O, et al. (2020) Techno-economic comparative analysis of renewable energy systems: Case study in Zimbabwe. Inventions 5: 27. https://doi.org/10.3390/inventions5030027 doi: 10.3390/inventions5030027
    [51] Tareq Salameh, Chaouki Ghenai, Adel Merabet, et al. (2020) Techno-economical optimization of an integrated stand-alone hybrid solar PV tracking and diesel generator power system in Khorfakkan, United Arab Emirates. Energy 190: 116475. https://doi.org/10.1016/j.energy.2019.116475 doi: 10.1016/j.energy.2019.116475
    [52] Sina Makhdoomi, Alireza Askarzadeh (2020) Optimizing operation of a photovoltaic/diesel generator hybrid energy system with pumped hydro storage by a modified crow search algorithm. J Energy Storage 27: 101040. https://doi.org/10.1016/j.est.2019.101040 doi: 10.1016/j.est.2019.101040
    [53] Javed MS, Ma T, Jurasz J, et al. (2020) Solar and wind power generation systems with pumped hydro storage: Review and future perspectives. Renewable Energy 148: 176-192. Available from: https://ideas.repec.org/a/eee/renene/v148y2020icp176-192.html.
    [54] Swaminathan Ganesan, Umashankar Subramaniam, Ghodke AA, et al. (2020) Investigation on sizing of voltage source for a battery energy storage system in microgrid with renewable energy sources. IEEE Access 8: 188861-188874. https://doi.org/10.1109/ACCESS.2020.3030729 doi: 10.1109/ACCESS.2020.3030729
    [55] Yashwant Sawle, Siddharth Jain, Sanjana Babu, et al. (2021) Prefeasibility economic and sensitivity assessment of hybrid renewable energy system. IEEE Access 9: 28260-28271. https://doi.org/10.1109/ACCESS.2021.3058517 doi: 10.1109/ACCESS.2021.3058517
    [56] Oladigbolu JO, Al-Turki YA, Olatomiwa L (2021) Comparative study and sensitivity analysis of a standalone hybrid energy system for electrification of rural healthcare facility in Nigeria. Alexandria Eng. J 60: 5547-5565. https://doi.org/10.1016/j.aej.2021.04.042 doi: 10.1016/j.aej.2021.04.042
    [57] Niyonteze JDD, Zou FM, Osarumwense Asemota GN, et al. (2020) Key technology development needs and applicability analysis of renewable energy hybrid technologies in off-grid areas for the rwanda power sector. Heliyon 6: e03300. https://doi.org/10.1016/j.heliyon.2020.e03300 doi: 10.1016/j.heliyon.2020.e03300
    [58] Faiz FUH, Shakoor R, Raheem A, et al. (2021) Modeling and analysis of 3 MW solar photovoltaic plant using PVSyst at Islamia University of Bahawalpur, Pakistan. Int J Photoenergy 2021: 1-14. http://dx.doi.org/10.1155/2021/6673448
    [59] Ajith Gopi, Sudhakar K, Ngui WK, et al. (2021) Performance modeling of the weather impact on a utility-scale PV power plant in a tropical region. Int J Photoenergy 2021: 1-10. https://doi.org/10.1155/2021/5551014 doi: 10.1155/2021/5551014
    [60] Imasiku K (2021) A solar photovoltaic performance and financial modeling solution for grid-connected homes in Zambia. Int J Photoenergy 2021: 1-13. https://doi.org/10.1155/2021/8870109 doi: 10.1155/2021/8870109
    [61] Mahmoud FE, Elkadeem MR, Kotb KM, et al. (2021) Optimal design and energy management of an isolated fully renewable energy system integrating batteries and supercapacitors. Energy Con Manage 245: 114584. https://doi.org/10.1016/j.enconman.2021.114584 doi: 10.1016/j.enconman.2021.114584
    [62] Ahmed Al-Sarraj, Kareem KM (2020) Simulation design of hybrid system (grid/PV/wind turbine/battery/diesel) with applying HOMER: A case study in Baghdad, Iraq. Int J Electron Commun Eng. 7: 10-18. https://doi.org/10.14445/23488549/IJECE-V7I5P103 doi: 10.14445/23488549/IJECE-V7I5P103
    [63] Natei Ermias Benti, Yedilfana Setarge Mekonnen, Ashenafi Abebe Asfaw, et al. (2022) Techno-economic analysis of solar energy system for electrification of rural school in Southern Ethiopia. Cogent Engineering 9: 1-21. https://doi.org/10.1080/23311916.2021.2021838 doi: 10.1080/23311916.2021.2021838
    [64] Md. Shafiuzzaman KK (2007) Solar and wind energy resource assessment (SWERA)-Bangladesh, UNEP/GEF. Available from: https://www.researchgate.net/publication/282665355_Solar_and_Wind_Energy_Resource_Assessment_SWERA_-_Bangladesh.
    [65] Rezzouk H, Mellit A (2015) Feasibility study and sensitivity analysis of a stand-alone photovoltaic-diesel-battery hybrid energy system in the north of Algeria. Renewable Sustainable Energy Rev 43: 1134-1150. https://doi.org/10.1016/j.rser.2014.11.103 doi: 10.1016/j.rser.2014.11.103
    [66] Abdul Momin BM (2019) Summary of rainfall in Bangladesh for the year 2017 & 2018. Surface Water Processing Branch. Available from: http://www.hydrology.bwdb.gov.bd/img_upload/ongoing_project/756.pdf.
    [67] Gebrehiwot K, Hossain Mondal MA, Claudia Ringler, et al. (2019) Optimization and cost-benefit assessment of hybrid power systems for off-grid rural electrification in Ethiopia. Energy 177: 234-246. https://doi.org/10.1016/j.energy.2019.04.095 doi: 10.1016/j.energy.2019.04.095
    [68] Division P (2018) "Net Metering Guidelines-2018", Sustainable and Renewable Energy Development Authority (SREDA), December 2018. Available from: https://www.bd.undp.org/content/dam/bangladesh/docs/Projects/srepgen/2018.11.28%20-%20Net%20Metering%20Guidelien%202018%20(English).pdf.
    [69] Mondal AH, Denich M (2010) Hybrid systems for decentralized power generation in Bangladesh. Energy Sustainable Develp 14: 48-55. https://doi.org/10.1016/j.esd.2010.01.001 doi: 10.1016/j.esd.2010.01.001
    [70] Jäger-Waldau A (2021) Snapshot of photovoltaics. EPJ Photovoltaics 12: 1-7. https://doi.org/10.1051/epjpv/2021002 doi: 10.1051/epjpv/2021002
    [71] Retail electricity tariff rate (2020) Bangladesh Energy Regulatory Commission. Available from: http://www.berc.org.bd/.
    [72] Inflation rate in Bangladesh (2021) Available from: https://www.statista.com/statistics/438363/inflation-rate-in-bangladesh/.
  • This article has been cited by:

    1. Zilin Xia, Jinan Gu, Wenbo Wang, Zedong Huang, Research on a lightweight electronic component detection method based on knowledge distillation, 2023, 20, 1551-0018, 20971, 10.3934/mbe.2023928
    2. Hanyi Shi, Ningzhi Wang, Xinyao Xu, Yue Qian, Lingbin Zeng, Yi Zhu, HeMoDU: High-Efficiency Multi-Object Detection Algorithm for Unmanned Aerial Vehicles on Urban Roads, 2024, 24, 1424-8220, 4045, 10.3390/s24134045
    3. Bao Liu, Wenqiang Jiang, RS-YOLOx: target feature enhancement and bounding box auxiliary regression based object detection approach for remote sensing, 2024, 18, 1931-3195, 10.1117/1.JRS.18.016514
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4761) PDF downloads(367) Cited by(12)

Figures and Tables

Figures(16)  /  Tables(3)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog