Loading [MathJax]/jax/element/mml/optable/GreekAndCoptic.js
Research article

Geostatistical modelling of soil contamination with arsenic, cadmium, lead, and nickel: the Silesian voivodeship, Poland case study

  • Received: 22 January 2020 Accepted: 23 March 2020 Published: 01 April 2020
  • The increasing number of spatial data sets permits their application for minimising the duration and cost of such research. An example of such application is geostatistical modelling. Data on the quality of atmospheric air can be used for the assessment of the quality of soil in a given area. The objective of the study was an attempt to apply geostatistical methods in the estimation of the degree of soil contamination with selected heavy metals resulting from deposition from atmospheric air. This paper uses data obtained from the State Environmental Monitoring (data on the quality of atmospheric air) and other collections of publicly available spatial data for the purpose of analysis of the state of quality of soils in the province of Silesia in Poland. Conducted analyses revealing that contamination with lead in atmospheric air in the Silesian voivodeship considerably exceeds acceptable values, and the load of lead deposition is largely transferred to the soil. The paper also presents geochemical maps necessary to understand sources of soil contamination, including their natural content.

    Citation: Jolanta Kwiatkowska-Malina, Andrzej Szymon Borkowski. Geostatistical modelling of soil contamination with arsenic, cadmium, lead, and nickel: the Silesian voivodeship, Poland case study[J]. AIMS Geosciences, 2020, 6(2): 135-148. doi: 10.3934/geosci.2020009

    Related Papers:

    [1] Rasool M. Imran, Kadhim Hamzah Chalok, Siraj A. M. Nasrallah . Innovative two-stage thermal control of DC-DC converter for hybrid PV-battery system. AIMS Electronics and Electrical Engineering, 2025, 9(1): 26-45. doi: 10.3934/electreng.2025002
    [2] Said Oucheriah, Abul Azad . Current-sensorless robust sliding mode control for the DC-DC boost converter. AIMS Electronics and Electrical Engineering, 2025, 9(1): 46-59. doi: 10.3934/electreng.2025003
    [3] José M. Campos-Salazar, Roya Rafiezadeh, Felipe Santander, Juan L. Aguayo-Lazcano, Nicolás Kunakov . Comprehensive GSSA and D-Q frame dynamic modeling of dual-active-bridge DC-DC converters. AIMS Electronics and Electrical Engineering, 2025, 9(3): 288-313. doi: 10.3934/electreng.2025014
    [4] Tarun Naruka, Debasis Tripathy, Prangya Mohanty . Power quality enhancement by a solar photovoltaic-based distribution static compensator. AIMS Electronics and Electrical Engineering, 2025, 9(2): 192-214. doi: 10.3934/electreng.2025010
    [5] Muamer M. Shebani, M. Tariq Iqbal, John E. Quaicoe . Comparison between alternative droop control strategy, modified droop method and control algorithm technique for parallel-connected converters. AIMS Electronics and Electrical Engineering, 2021, 5(1): 1-23. doi: 10.3934/electreng.2021001
    [6] Suliana Ab-Ghani, Hamdan Daniyal, Abu Zaharin Ahmad, Norazila Jaalam, Norhafidzah Mohd Saad, Nur Huda Ramlan, Norhazilina Bahari . Adaptive online auto-tuning using Particle Swarm optimized PI controller with time-variant approach for high accuracy and speed in Dual Active Bridge converter. AIMS Electronics and Electrical Engineering, 2023, 7(2): 156-170. doi: 10.3934/electreng.2023009
    [7] I.E.S. Naidu, S. Srikanth, A. Siva sarapakara Rao, Adabala Venkatanarayana . A novel mine blast optimization algorithm (MBOA) based MPPT controlling for grid-PV systems. AIMS Electronics and Electrical Engineering, 2023, 7(2): 135-155. doi: 10.3934/electreng.2023008
    [8] Mohamed M. Albarghot, Mohamed T. Iqbal, Kevin Pope, Luc Rolland . Dynamic modeling and simulation of the MUN Explorer autonomous underwater vehicle with a fuel cell system. AIMS Electronics and Electrical Engineering, 2020, 4(1): 114-131. doi: 10.3934/ElectrEng.2020.1.114
    [9] José M. Campos-Salazar, Juan L. Aguayo-Lazcano, Roya Rafiezadeh . Non-ideal two-level battery charger—modeling and simulation. AIMS Electronics and Electrical Engineering, 2025, 9(1): 60-80. doi: 10.3934/electreng.2025004
    [10] José M. Campos-Salazar, Roya Rafiezadeh, Juan L. Aguayo-Lazcano, Constanza Márquez . Reduction of harmonic distortion in electromagnetic torque of a single-phase reluctance motor using a multilevel neutral-point-clamped DC-AC converter. AIMS Electronics and Electrical Engineering, 2025, 9(2): 215-242. doi: 10.3934/electreng.2025011
  • The increasing number of spatial data sets permits their application for minimising the duration and cost of such research. An example of such application is geostatistical modelling. Data on the quality of atmospheric air can be used for the assessment of the quality of soil in a given area. The objective of the study was an attempt to apply geostatistical methods in the estimation of the degree of soil contamination with selected heavy metals resulting from deposition from atmospheric air. This paper uses data obtained from the State Environmental Monitoring (data on the quality of atmospheric air) and other collections of publicly available spatial data for the purpose of analysis of the state of quality of soils in the province of Silesia in Poland. Conducted analyses revealing that contamination with lead in atmospheric air in the Silesian voivodeship considerably exceeds acceptable values, and the load of lead deposition is largely transferred to the soil. The paper also presents geochemical maps necessary to understand sources of soil contamination, including their natural content.


    Power converters are increasingly employed in numerous applications, such as industrial electronics, portable devices, energy storage systems, and consumer electronics. In recent applications, such as artificial intelligence, networks related to the Internet of Things, and wireless sensors, power converters are critical components. Switch-mode DC-DC converters transfer the DC voltage from one level to another. High-level DC voltage is reduced to a low-level DC voltage by a DC-DC buck converter, which is widely used in modern-era applications such as photovoltaic systems, electric vehicles etc.

    Regulating the output voltage of DC-DC converters in the presence of external noise and system parameter variations has proven to be challenging. Rigorous investigations are still ongoing in this subject, and different control algorithms have been proposed for controlling the output voltage of the converter. Chebyshev's neural network-based adaptive controller was designed to control the output voltage of a buck converter [1]. In another work, a new control algorithm was proposed to track the output voltage of a buck converter with overcurrent protection to hardware elements [2]. In Miao (2019) and Wang (2019), current-constrained controllers were proposed for the buck converter [3,4]. A linear matrix inequality-based H-infinity controller has been demonstrated to regulate the output voltage of buck and boost converters [5,6]. An artificial neural network (ANN)-based control algorithm was proposed for a DC-DC buck converter, and an optimal controller was designed using approximate dynamic programming [7]. Several observer-based control algorithms have been proposed, both for buck and boost converters [8,9,10,11]. A data-driven control algorithm was suggested for the interleaved DC-DC boost converter based on online integral reinforcement learning [12]. In 2022, the non-minimum-phase model of a boost converter was modified to a minimum-phase model for the purpose of controller design [13].

    The proportional-integral-derivative (PID) controller and its modified versions are widely used in different applications. A sigmoid function-based data-driven proportional-integral (PI) controller was proposed to control the angular speed of a DC motor that is energized by a DC-DC converter [14]. Here, a sigmoid function represents the limits of the parameters of the PI controller. In [15], a neuroendocrine PID controller was designed for a multivariable crane system where the hormone discharge rate of neuroendocrine-PID can be regulated as per the error change. The output voltage of the DC-DC boost converter has been tracked using a brain emotional learning-based intelligent controller (BELBIC), where parameters are determined using the particle swarm optimization technique [16]. The performance of this controller outperformed the classical PID controller. Fractional-order proportional-integral (FOPI) and fractional-order proportional-integral-derivative (FOPID) controllers are very popular in industrial applications [17,18,19]. The FOPID controller is designed to stabilize the output power of the buck-boost converter for proton exchange membrane fuel cells [20]; an optimal FOPI controller was proposed for a solar photovoltaic system with a DC-DC boost converter [21]. The FOPID controller was also implemented in a real-time system, like a DC-DC boost converter and a magnetic levitation system [22,23]. A modified version of the FOPID controller, namely the tilt integral derivative (TID) and hybrid proportional integral derivative (PID) controllers, was proposed for load frequency control in multi-area interconnected power systems and second-order time delay systems, respectively [24,25].

    Though different control techniques have been proposed for tracking the output voltage of the DC-DC converter, the Kalman filter-based control algorithm has not been addressed before. The measurement of the output voltage of a DC-DC converter is essential for control purposes. However, measured outputs may be affected by the sensor noise. To deal with model uncertainty, noise must be considered. As a result, this research work designed a Kalman filter (KF) along with a controller to reduce the effects of sensor and process noises on the converter output. Modified versions of PID controllers have been employed for DC-DC converters, but some complications limit their widespread use. Instead, fractional-order PID controllers and variants have shown superior control abilities. Fractional-order proportional-integral (FOPID) controllers suffer from tuning complexity, and TID controllers are sensitive to parameter variations. To avoid these problems, a KF-based new controller, named blended proportional-integral (BPI) controller, is proposed to track the output voltage of the DC-DC buck converter. The BPI controller is a combination of FOPI and integer-order integral (I) controllers.

    For parameter determination of the PID and FOPID controllers, optimization techniques are essential; for that purpose, metaheuristic algorithms like particle swarm optimization (PSO), artificial bee colony (ABC), genetic algorithm (GA), or firefly algorithm (FA) play a critical role [26]. The PID, FOPID, and FOPI parameters have been determined via the PSO technique [27], and a modified version of the sine-cosine algorithm was employed to estimate the PID controller parameters [28]. In [29], six parameters of a complex controller, which was a combination of the PID controller with an N-filter and a PD controller, were predicted using the golden eagle optimization technique. A new metaheuristic optimization technique, which is a combination of pattern search with opposition-based snake optimizer, was introduced and utilized to calculate the parameters of two cascaded first-order low-pass filters with the PID (PID-FF) controller [30]. Due to fast convergence rate, ease of implementation, and simplicity, the firefly algorithm (FA) was applied to estimate the parameters of the proposed BPI controller, which was introduced by Xin-She Yangin, 2008 [31]. This FA was then modified by different researchers and has been used in different applications for optimization purposes [32,33,34,35].

    The contributions of the proposed work are presented below.

    ●  A new Kalman filter-based controller is designed for the DC-DC buck converter. To control the output voltage of the converter, a blended proportional-integral (BPI) controller is proposed, which is the fractional order proportional-integral (FOPI) plus integer integral (I) controller.

    ●  The parameters of the suggested controller are estimated using four metaheuristic optimization techniques: firefly algorithm (FA), artificial bee colony (ABC), particle swarm optimization (PSO), and Harris Hawks optimization (HHO). The FA-based controller outperforms the other three metaheuristic algorithms.

    ●  The Kalman filter minimizes the effect of sensor and process noises on the output of the system, while the Kalman filter-based BPI controller tracks the output voltage of the converter with excellent transient response specifications in the presence of sensor and process noises. When the performance of the suggested BPI controller is compared with the PI, PID, FOPI, and FOPID controllers, BPI outperforms the others. The output of the closed-loop system is also compared with the artificial neural network-based controller presented by Dong et al., 2021 [7]; the proposed controller outperforms this controller.

    Figure 1 represents the circuit diagram of the buck converter. The circuit elements include the controlled switch MOSFET (SBu), diode (DBu), inductor (LBu) with internal resistance RL, and capacitor (CBu) with internal resistance RC, load resistance R, and DC input voltage Vin. The converter circuits with controlled switch (SBu) during on and off conditions are also shown. The parameter values for the buck converter circuit are given in Table 1 and are based on [7].

    Figure 1.  Buck converter circuit (a); on-mode circuit (b); off-mode circuit (c).
    Table 1.  Values of different parameters of the buck converter [7].
    Parameter Value
    Input voltage, Vin (𝑉) 42
    Load resistance, R (Ω) 10
    Internal resistance of inductor, RL (Ω) 0.3
    Internal resistance of capacitor, RC (Ω) 0.02
    Inductance, LBu(m𝐻) 5.63
    Capacitance, CBu (μF) 5

     | Show Table
    DownLoad: CSV

    The voltage current of the converter during the on state is represented by the following equations:

    Vin=RLiL+LBudiLdt+v0 (1a)
    CBudvcdt=iLvoR (1b)
    vo=RC(iLvoR)+vc (1c)

    where iL and vc are the inductor current and capacitor voltage, respectively, and vo is the output voltage.

    The equations below characterize the converter during the off state:

    RLiL+LBudiLdt+v0=0 (2a)
    CBudvcdt=iLvoR (2b)
    vo=RC(iLvoR)+vc (2c)

    Thus, the average model of the converter is given by the following equations:

    DVin=RLiL+LBudiLdt+v0(t) (3a)
    CBudvcdt=iLvoR (3b)
    vo=RC(iLvoR)+vc (3c)

    Here, Drepresents the duty ratio (control input), which is provided by the pulse width modulator (PWM), andD[0,1].

    The state-space model of the converter is illustrated below:

    [diLdtdvcdt]=[1LBu(RRCR+Rc+RL)RLBu(R+Rc)1CBu(1+RCR+Rc)1(R+Rc).CBu][iLvc]+[VinLBu0]D (4)
    vo=[RRcR+RcRR+Rc][iLvc] (5)

    The Kalman filter-based controller is used to regulate the output voltage of the DC-DC buck converter. The parameters of the proposed BPI controller are determined by applying the metaheuristic firefly optimization technique. The block diagram of the buck converter with the Kalman filter-based controller is shown in Figure 2. The Kalman filter eliminates the effects of sensor and process noise on the output voltage, which is used for comparison with the reference voltage. The control input u(t) is applied to the buck converter through the PWM.

    Figure 2.  Buck converter with the Kalman filter-based BPI controller.

    The Kalman filter identifies the state variables of the system from the input and output data. This estimated output is used for control purposes. R. E. Kalman proposed the idea of the Kalman filter in 1960 [36]; since then, the Kalman filter and its modified versions have been used extensively in real-time systems for state estimation [37,38,39,40,41].

    In discrete time domain, the state-space model of a linear time invariant is represented as

    xk+1=Axk+Buk+wk (6)
    yk=Cxk+nk=yt+nk (7)

    Here, k represents the time instant; xRnx, uRnu, and yRny are the state, input, and measurement vectors, respectively; yt is the output without sensor noise; wk and nk are the process and sensor noise, respectively.

    It is assumed that process (wk) and sensor (nk) noises are Gaussian white noises with zero-mean and are mutually uncorrelated signals. The process and sensor noise covariance matrices are E(wkwTk)=Q and E(nknTk)=R, respectively. The algorithm of the Kalman filter is given in Table 2. The Kalman filter estimates the state variables of the system from the information of input and noisy output of the system. It estimates the state variables by minimizing the mean squared error E[(xkˆxk)(xkˆxk)T]=PK where ˆxk is the prior estimated state vector and PK is the error covariance matrix at k. These estimated states are free from sensor and process noise and are used for control purposes. Accordingly, the Kalman filter minimizes the effects of noise on the output of the system. The algorithm of the Kalman filter is given in Table 2 [37].

    Table 2.  Kalman filter (KF) algorithm.
    Step 1:
    Initialization (k = 0)
    Mean value of states: ˆx0=E(x0)=xin
    Error covariance matrix: P0=E[(x0ˆx0)(x0ˆx0)T]=Pin
    Step 2:
    Prediction
    (time update)
    Predicted state vector: ˆxk+1k=Aˆxkk+Buk
    A priori covariance matrix: Pk+1k=APkkAT+Qw
    Step 3:
    Correction
    (measurement update)
    The updated state vector: ˆxk+1k+1=ˆxk+1k+Kk+1[yk+1kCˆxk+1k]
    Posteriori covariance matrix: Pk+1k+1=(IKk+1C)Pk+1k
    The Kalman filter gain: Kk+1=Pk+1kCT(Qv+CPk+1kCT)1
    The algorithm returns to step 2.

     | Show Table
    DownLoad: CSV

    It is challenging to select the process noise covariance matrix Q as it depends on model uncertainty. For this, Q is chosen by the trial-and-error method as diag(0.00001,0.00001). The sensor noise covariance (R) is chosen as 0.1, which depends on the accuracy of the voltage sensor. The Kalman filter is implemented using MATLAB software, and the state estimation error covariance matrix becomes [2.728×1072.725××1082.725×1082.975×109]. The resulting Kalman filter estimates the state variables iL and vc accurately, and accuracy is computed in terms of root mean square error (RMSE) given by RMSE=(Tsk=1(x1ˆx1)2Ts)1/2, where x1 and ˆx1 are the actual and estimated state variables, respectively. The RMSE values of state variables iL and vc are 0.1261 and 0.2167, respectively. The covariance of the measured error (error between yt and measured output, y) and the estimated error (error between yt and estimated output, ˆy) are 0.0994 and 0.0465, respectively. From these covariances, it can be concluded that KF minimizes the effect of noise on the output of the system.

    Linear controllers, as the integer and fractional-order proportional-integral (PI) and proportional-integral-derivative (PID) controllers, are used extensively for the fulfillment of the transient as well as steady-state response specifications of the system. The transfer function of PID and fractional-order PID (FOPID) controllers is as follows [17,18].

    GPID(s)=kp+kis+kds;GFOPID(s)=kp+kisλ+kdsμ

    where kp, kd, and ki are the gain of the proportional, derivative, and integral controllers, respectively; λ and μ are real fractional numbers (0(λ,μ)1).

    The FOPID has more control on the system than PID, since the former has five control parameters, as opposed to the three from PID. Still, the issue lies in optimizing these five parameters in FOPID instead of the three in PID.

    For this study, a new fractional-order controller is proposed, which is a combination of integer and fractional-order PI controllers and is named the blended proportional-integral (BPI) controller. The structure of the controller is shown in Figure 3, and the transfer function of the controller is given by

    GBPI(s)=kp+kis+kfsλ

    where kp is the gain of the proportional controller, ki is the gain of the integer-order integral controller, kf is the gain of the fractional-order integral controller, and the exponent λ is a real number (0λ1).

    Figure 3.  Structure of BPI controller.

    The number of parameters that need to be optimized for the BPI controller is four, which is less than the FOPID (five parameters). The merit of the BPI controller over other fractional-order controllers like FOPID and TID is that it has no derivative term. Actually, the derivative term is required to damp out oscillations, but it can strengthen noise if parameter estimations are incorrect. The results section proves that BPI presents robust performance when facing parameter variations and noise.

    The firefly algorithm (FA) is an extremely popular metaheuristic optimization technique that takes inspiration from nature and is extensively used for solving engineering problems. In this algorithm, the objective function changes proportionally with the brightness of the firefly; this means that when brightness is reduced, the firefly will travel or search for a brighter firefly. If there is no brighter firefly, then the firefly will move aimlessly.

    Consider there are n fireflies where f(d) represents the fitness function or objective function. The brightness of a firefly (Ii) determines the current position di. The distance between two fireflies is given [31,32] by

    rij=didj=k=nk=1(dikdjk)2 (8)

    where di and dj represent the position of ith and jth fireflies, respectively.

    The relationship between the attractiveness function β and distance r is β=β0eγr2, where β0 signifies the attractiveness at r=0 and the light absorption coefficient is denoted by γ.

    Suppose that in a t-th iteration, the firefly i at position di is moved to a brighter firefly j. Now, the new position of firefly i in (t+1)-th iteration is expressed by the following equation [33,34]:

    {d}_{i}(t+1) = {d}_{i}\left(t\right)+{\beta }_{0}{e}^{-\gamma {r}_{ij}^{2}}\left({d}_{j}\left(t\right)-{d}_{i}\left(t\right)\right)+\alpha {ϵ}_{i} (9)

    where {d}_{i}\left(t\right) represents the old position of firefly i; {ϵ}_{i} is a random number that can be determined from Gaussian, uniform, or any other distribution; \alpha is called the randomization parameter \left(0\le \alpha \le 1\right) .

    The pseudocode of FA is illustrated in Table 3. Optimization problems like extremely nonlinear and multi-modal can be handled efficiently using FA; its speed of convergence is also extraordinary.

    Table 3.  Algorithm of the firefly optimization technique.
    Pseudocode of firefly algorithm (FA)
    Input: Number of fireflies (n), no. of maximum iterations (MaxIter), attractiveness factor ( {\beta }_{0}) , randomization parameter ( \alpha ), light absorption coefficient ( \gamma ).
    Output: Best result originated by the FA.
    1.       Initialization of n, MaxIter, {\beta }_{0} , \alpha , \gamma .
    2.       Determine initial population of fireflies {d}_{i} , where i=\mathrm{1, 2}, \dots, n .
    3.       Define objective function f\left(d\right) where d={[{d}_{1}, \dots, {d}_{n}]}^{T} .
    4.       Calculate light intensity {I}_{i} for firefly at {d}_{i} using f\left({d}_{i}\right) for all fireflies
    Main Loop
    5. While ( t < MaxIter )
                For i=1\;\mathrm{t}\mathrm{o}\;n
                For j=1\;\mathrm{t}\mathrm{o}\;n
                if ({I}_{j} > {I}_{i})
                Travel firefly i towards firefly j according to equation (9)
                                    end if
                                    Evaluate new f\left({d}_{i}\right) and light intensity {I}_{i}
                                end for j
                end for i
    Rank the fireflies and determine the best global solution
                    end while

     | Show Table
    DownLoad: CSV

    For the optimization, the fitness function or cost function is considered as integral time absolute error (ITAE), which computes the integral of the product of the absolute value of difference between reference input and actual output and time. The ITAE is given by

    f\left({k}_{p}, {k}_{i}, {k}_{f, }\lambda \right) = ITAE = {\int }_{0}^{t}t\left|{v}_{r}-{\widehat{v}}_{o}\right|dt

    where {v}_{r} and {\widehat{v}}_{o} are the reference voltage and estimated output voltage by Kalman filter respectively.

    One of the key objectives of this research work is to minimize the time response specifications, like maximum peak overshoot, settling time, and steady-state error. The cost function ITAE fulfils these requirements better than other cost functions, like integral square error (ISE) and integral absolute error (IAE). The design procedure of the BPI controller using the FA algorithm is presented in Figure 4.

    Figure 4.  Design procedure of the BPI controller using the FA algorithm.

    The values of different parameters used for the buck converter are presented in Section 2.1 [7]. The transfer function of the considered DC-DC buck converter is G\left(s\right) = \frac{(148.9s+1.489\times {10}^{9})}{({S}^{2}+2.002\times {10}^{4}s+3.652\times {10}^{7})} .

    For the BPI controller, four parameters, {k}_{p} , {k}_{i} , {k}_{f} , and \lambda need to be optimized using the optimization technique FA. The lower limit (LL) and upper limit (UL) of these four parameters are chosen from experience as follows:

    {k}_{p} = \left[\begin{array}{cc}0.001\left(LL\right)& 0.09\left(UL\right)\end{array}\right] ; {k}_{i} = \left[\begin{array}{cc}0.001\left(LL\right)& 10\left(UL\right)\end{array}\right] ; {k}_{f} = \left[\begin{array}{cc}0.001\left(LL\right)& 6\left(UL\right)\end{array}\right] ; \lambda = \left[\begin{array}{cc}0.001\left(LL\right)& 0.999\left(UL\right)\end{array}\right] .

    Applying FA, the gains of the controller BPI are optimized considering the cost function ITAE, ISE, and IAE. The best cost function versus the number of iteration curves is shown in Figure 5a for 25 iterations.

    Figure 5.  Fitness function vs. number of iteration curve (a) for different cost functions and (b) for different controllers.

    In this case, PI, PID, FOPI, and FOPID controllers are also designed for the buck converter using FA to compare the potential of BPI. Table 4 represents the values of the parameters of the different controllers. The objective function ITAE vs. the number of iterations (25 iterations) of different controllers is shown in Figure 5b; it shows that, among all controllers, BPI provides the minimum value of ITAE. The fractional-order controller blocks are obtained using the FOMCON toolbox. Figure 6 represents the closed-loop output voltage and inductor current responses of the buck converter with Kalman filter-based controllers for an 18 V reference voltage. Performance measures of different controllers are given in Table 5.

    Table 4.  Parameters of different controllers obtained using different metaheuristic algorithms.
    Optimization technique Control algorithm Value of the parameters
    {k}_{p} {k}_{i} {k}_{d} {k}_{f} \lambda \mu
    PI 0.001 9.8729 - - - -
    PID 0.0998 10 0.0362 - - -
    FOPI 0.0012 9.9803 - - 0.999 -
    FA FOPID 0.001 9.9317 0.0157 0.999 0.001
    BPI 0.0017 9.99 - 5.9858 0.999 -
    ABC BPI 0.0112 10 - 6 0.999 -
    PSO BPI 0.001 10 - 6 0.999 -
    HHO BPI 0.09 10 - 6 0.999 -

     | Show Table
    DownLoad: CSV
    Figure 6.  Output voltage (a) and inductor current (b) of the buck converter with different controllers.
    Table 5.  Performance measures of different controllers using different metaheuristic algorithms.
    Time response specifications Settling-time ( \mathit{s}) Maximum overshoot Under shoot Integral time absolute error (ITAE) Integral square error (ISE) Integral absolute error (IAE)
    FA-based PI 0.0095 0 0 9.108×10-5 0.4736 0.04471
    PID 0.001 4.54% 0 0.0007329 0.05703 0.03065
    FOPI 0.0095 0 0 9.125×10-5 0.4636 0.04412
    FOPID 0.02 0 0 0.0001613 0.3343 0.04433
    BPI 0.0055 0.06% 0 3.101×10-5 0.3177 0.02766
    ABC-based BPI 0.01 0% 0 4.672×10-5 0.2337 0.02756
    PSO-based BPI 0.0051 0.236% 0 3.044×10-5 0.3261 0.02774
    HHO-based BPI 0.028 0% 0 0.0001776 0.08423 0.02754

     | Show Table
    DownLoad: CSV

    To test the performance of the Kalman filter, outputs are collected for a reference voltage of 18V with a sensor noise power of 0.1 and process noise power of 0.000001. The output voltage and inductor current of the buck converter with a KF-based BPI controller and without KF (only with BPI controller) are shown in Figure 7.

    Figure 7.  Output voltage of the closed-loop system with and without the Kalman filter.

    Besides the firefly algorithm (FA), three other metaheuristic algorithms [artificial bee colony (ABC), particle swarm optimization (PSO), and Harris Hawks optimization (HHO)] were used to estimate parameters ( {k}_{p} , {k}_{i} , {k}_{f} , and \lambda ) of the BPI controller. The algorithms of these optimization techniques are given in the literature [42]. The estimated parameters of the BPI controller using these four optimization techniques are presented in Table 4. Performance measures of different optimization techniques for the KF-based BPI controllers are given in Table 5 and shown in Figure 8.

    Figure 8.  Performance measures of different algorithm-based BPI controllers with KF.

    The statistical analysis of the fitness function (ITAE) with different optimization algorithms is presented in Table 6. As the metaheuristic algorithms FA, PSO, ABC, and HHO are stochastic in nature, the best value, standard deviation, and mean value of ITAE are calculated for 15 individual trials. The p-value of Wilcoxon's rank-sum test is also calculated and provided in Table 6. Wilcoxon's test, a nonparametric test, is used to compare two paired groups. A p-value lower than 0.05 means there is no variance between the two compared groups.

    Table 6.  Statistical analysis of fitness function (ITAE) for different optimization algorithms.
    Algorithm Best Standard deviation Mean Wilcoxon's rank-sum test p-value FA (BPI) vs
    FA (BPI) 1.75× 10-5 1.81548× 10-6 1.88× 10-5
    ABC (BPI) 1.77× 10-5 6.36946× 10-6 2.54× 10-5 6.0111× 10-4
    HHO (BPI) 1.78× 10-5 4.55625× 10-6 2.34× 10-5 0.0012
    PSO (BPI) 1.77× 10-5 5.10655× 10-6 2.30× 10-5 0.0047

     | Show Table
    DownLoad: CSV

    For comparison with the proposed KF-based BPI controller, an artificial neural network (ANN) -based controller, suggested by Dong (2021), was considered. The same system is considered, and simulation results are compared in both cases. The KF-based BPI controller is tested with the same load and parameter variations as mentioned in [7]. Responses under various conditions are illustrated in Figure 9, and performance measures are given in Table 7.

    Figure 9.  Output voltage and inductor current responses when (a) load changes from 7.33 to 11 Ω; (b) reference voltage varies from 18 to 24 V; (c) input voltage {V}_{in} changes from 42 to 47 V; (d) {L}_{Bu} decreases by 50% with the same load change as in (a); (e) {L}_{Bu} decreases by 85% with the same load change as in (a).
    Table 7.  Performance comparison with the ANN-based controller.
    Changed condition % Overshoot Settling time (s)
    ANN controller Proposed controller ANN controller Proposed controller
    Load resistance R changes from 7.33 to 11 Ω 38.89% 1.39% 0.002 0.002
    Reference voltage {v}_{r} varies from 18 to 24 V 0% 0% 0.007 0.005
    Input voltage {V}_{in} is changed from 42 to 47 V 4.44% 17.2% 0.006 0.0035
    Inductance {L}_{Bu} decreases by 50% and load resistance R changes from 7.33 to 11 Ω 38.89% 1.39% 0.005 0.003
    Inductance {L}_{Bu} decreases by 85% and load resistance R changes from 7.33 to 11 Ω 44.44% 1.66% 0.0017 0.0018

     | Show Table
    DownLoad: CSV

    Table 6 shows that the proposed BPI controller performs better than the ANN-based controller in terms of overshoot and settling time. Some ringing (oscillations) can be seen in both voltage and current waveforms for the ANN-based controller, which is not present at all in the proposed BPI controller.

    The best cost function values for IAE, ISE, and ITAE were 0.011533, 0.010004, and 1.75 × 10-5, respectively. From the convergence curve (Figure 5a) and these cost values, it is evident that ITAE provides better results than the other cost functions.

    From Table 5, it is clear that the proposed BPI controller demonstrates better performance in terms of settling time, ITAE, and IAE than the other designed controllers. Figure 6 shows that the output voltage with the KF-based BPI controller accurately tracks the reference voltage of 18 V, with a negligible amount of overshoot in both the output voltage and inductor current waveforms. Hence, the BPI controller is considered for the regulation of the output voltage of the buck converter.

    From Figure 7, it is clear that the converter outputs with the KF-based BPI controller are not influenced by the sensor and process noise; Additionally, the controller provides overcurrent protection to the hardware components of the buck converter. This indicates that the KF eliminates the influence of sensor and process noise on the system output.

    Table 5 shows that the ABC and HHO-based BPI controllers are capable of producing zero overshoot, but the settling time are higher than with the FA and PSO-based controllers. Also, the maximum peak overshoot of the PSO-based controller is greater than that of the FA-based controller. By comparing these performance measures, it is clear that the FA-based BPI controller outperforms the other three metaheuristic algorithm-based controllers. From Table 6, it is clear that the FA-based BPI controller provides the lowest best value, standard deviation, and mean value compared to the other techniques.

    Figure 9 shows the responses of the KF-based BPI controller under various parameter variations. A performance comparison between the proposed controller and the ANN-based controller [7] is presented in Table 7. The BPI controller improves maximum overshoot by 96.42% compared to the ANN controller under load resistance variations.

    This study describes the design of a new Kalman filter-based blended proportional-integral (BPI) controller and demonstrates the regulation of the output voltage of a DC-DC buck converter. The responses with and without the Kalman filter reveal that it eliminates the effect of sensor and process noises on the output voltage. The proposed fractional-order KF-based BPI controller gives better time response with respect to settling time, maximum overshoot, and ITAE than other KF-based controllers like PI, PID, FOPI, and FOPID. Statistical analysis shows that the performance of the proposed FA-based controller is superior to that of other metaheuristic algorithms. A comparative study also shows that the KF-based BPI controller outperforms other control techniques like ANN: the peak value of the voltage and transient time is much lower for the BPI controller than the ANN controller under different load conditions and parameter variations. The BPI controller decreases the maximum overshoot by 96.26% when the inductance of the buck converter is reduced by 85% under a load change from 7.33 Ω to 11 Ω. Hence, the proposed controller is capable to provide excellent transient responses and robust performance. In the future, the proposed control algorithm may be experimentally validated with a real-time system.

    Debarchita Mishra: Conceptualization, simulation, writing – original draft; Sharmistha Mandal: Analysis, drafting and editing of the manuscript.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    On behalf of all authors, the corresponding author declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.



    [1] Modrzewska B, Wyszkowski M (2014) Trace metals content in soils along the state road 51 (north-eastern Poland). Environ Monit Assess 186: 2589-2597. doi: 10.1007/s10661-013-3562-z
    [2] Osakwe SA (2013) Chemical partitioning of iron, cadmium, nickel and chromium in contaminated soils of south-eastern Nigeria. Chem Speciation Bioavailability 25: 71-78. doi: 10.3184/095422913X13581872822530
    [3] Bień JD, Meer J, Rulkens WH, et al. (2005) A GIS-based approach for the long-term prediction of human health risks at contaminated sites. Environ Model Assess 9: 221-226. doi: 10.1007/s10666-005-0909-z
    [4] Jalali M, Karami S, Marj AF (2016) Geostatistical Evaluation of Spatial Variation Related to Groundwater Quality Database: Case Study for Arak Plain Aquifer, Iran. Environ Model Assess 21: 707-719. doi: 10.1007/s10666-016-9506-6
    [5] Krami LK, Amiri F, Sefiyanian A, et al. (2013) Spatial patterns of heavy metals in soil under different geological structures and land uses for assessing metal enrichments. Environ Monit Assess 185: 9871-9888. doi: 10.1007/s10661-013-3298-9
    [6] Goovaerts P (1997) Geostatistics for natural resources evaluation. Oxford Uni. Press, New York, 483.
    [7] Lin WC, Lin YP, Wang YC (2016) A decision-making approach for delineating sites which are potentially contaminated by heavy metals via joint simulation. Environ Pollut 211: 98-110. doi: 10.1016/j.envpol.2015.12.030
    [8] Shen F, Liao R, Ali A, et al. (2017) Spatial distribution and risk assessment of heavy metals in soil near a Pb/Zn smelter in Feng County, China. Ecotoxicol Environ Saf 139: 254-262. doi: 10.1016/j.ecoenv.2017.01.044
    [9] Xie Y, Chen TB, Lei M, et al. (2011) Spatial distribution of soil heavy metal pollution estimated by different interpolation methods: Accuracy and uncertainty analysis. Chemosphere 82: 468-476. doi: 10.1016/j.chemosphere.2010.09.053
    [10] Zawadzki J, Fabjańczyk P (2012) Geostatistical evaluation of lead and zinc concentration in soils of an old mining area with complex land management. Int J Environ Sci Technol 10: 729-742 doi: 10.1007/s13762-012-0132-9
    [11] Namysłowska-Wilczyńska B (2019) Application of Geostatistical Techniques for the Determining of an Anomalous Zone of Copper Ore Deposit in the Area of Polkowice Mine (Region of Lubin-Sieroszowice, SW Part of Poland). Geoinfor Geostat Overv 7: 1-22
    [12] Namysłowska-Wilczyńska B, Wilczyński A (2002) Multivariate estimation and simulationfor environmental data modelling: processing of heavy metals concentration data in soil. Data Sci J 1: 27-44 doi: 10.2481/dsj.1.27
    [13] Ripin SNM, Hasan S, Kamal M (2014) Environmental Geochemical Mapping on Distribution of Metal Contamination in Topsoils Perlis, Malaysia. J Med Bioeng 3: 277-281.
    [14] Cai L, Xu Z, Bao P, et al. (2015). Multivariate and geostatistical analyses of the spatial distribution and source of arsenic and heavy metals in the agricultural soils in Shunde, Southeast China. J Geochem Explor 148: 189-195. doi: 10.1016/j.gexplo.2014.09.010
    [15] Salminen R, Tarvainen T, Demetriades A, et al. (1998) FOREGS Geochemical Mapping Field Manual. Geological Survey of Finland.
    [16] De Vivo B, Lima A, Bove AM, et al. (2006) Environmental Geochemical Atlas of Italy. European Congress on Regional Geoscientific Cartography and Earth Information and Systems Water, 126-127.
    [17] Guillén MT, Delgado J, Albanese S, et al. (2011) Environmental geochemical mapping of Huelva municipality soils (SW Spain) as a tool to determine background and baseline values. J Geochem Explor 109: 59-69. doi: 10.1016/j.gexplo.2011.03.003
    [18] Frentiu T, Ponta M, Levei E, et al. (2008) Preliminary study on heavy metals contamination of soil using solid phase speciation and the influence on groundwater in Bozanta-Baia Mare Area, Romania. Chem Speciation Bioavailability 20: 99-109. doi: 10.3184/095422908X324471
    [19] Malina G (2004) Ecotoxicological and environmental problems associated with the former chemical plant in Tarnowskie Gory, Poland. Toxicology 205: 157-172. doi: 10.1016/j.tox.2004.06.064
    [20] Saito H, Goovaerts P (2000) Geostatistical Interpolation of Positively Skewed and Censored Data in a Dioxin-Contaminated Site. Environ Sci Technol 34: 4228-4235. doi: 10.1021/es991450y
    [21] Fowler D, Skiba U, Nemitz E, et al. (2004). Measuring Aerosol and Heavy Metal Deposition on Urban Woodland and Grass Using Inventories of 210Pb and Metal Concentrations in Soil. Water Air Soil Pollut Focus 4: 483-499. doi: 10.1023/B:WAFO.0000028373.02470.ba
    [22] Brunner PH, Rechberger H (2016) Practical Handbook of Material Flow Analysis. CRC Press, 173-178.
    [23] Pleim J, Venkatram A, Yamartino R (1984) ADOM/TADAM model development program. ERT P-b980-520, prepared for OME AES of Canada and the Umweltbundesamt, West Germany, 111.
    [24] Sportisse B (2007) A review of parameterizations for modelling dry deposition and scavenging of radionuclides, Atmos Environ 41: 2683-2698.
    [25] Hopenhayn C (2006) Arsenic in Drinking Water: Impact on Human Health. Elements 2: 103-107. doi: 10.2113/gselements.2.2.103
    [26] Kapaj S, Peterson H, Liber K, et al. (2006) Human Health Effects From Chronic Arsenic Poisoning-A Review. J Environ Sci Health Part A Toxic/Hazard Subst Environ Eng 41: 2399-2428. doi: 10.1080/10934520600873571
    [27] Zawierucha I, Kozłowski C, Malina G (2016) Immobilized materials for removal of toxic metal ions from surface/groundwaters and aqueous waste streams. Environ Sci Processes Impacts 18: 429. doi: 10.1039/C5EM00670H
    [28] Kucharczak-Moryl E, Moryl A, Żmuda R (2014) Environment influence arsenic content in croplands in the zgorzelecki-bogatyński region. Ecological Engineering, No 37: 107-116. Available from: http://www.ineko.net.pl/Wplyw-srodowiska-na-zawartosc-arsenu-w-glebach-uprawnych-rejonu-zgorzelecko-bogatynskiego,283,0,1.html. Accessed 5 May 2016.
    [29] Mikulski S (2010) Gold, arsenic and dopeore. Polish Geological Institute-National Research Institute.
    [30] Kuys K, Ralston J, Sobieraj S, et al. (1989) Methods for the enrichment of bauxite ultrafines. Physicochem Metall Probl 21: 23-26.
    [31] Bernard A (2008) Cadmium & its adverse effects on human health. Indian J Med Res 128: 557-564.
    [32] Kabata-Pendias A, Piotrkowska M (1997) Transfer of cadmium, zinc and lead from soils to plants. The 5th international conference proceedings: transport, fate and effects of silver in the environment, 171-176.
    [33] Kwiatkowska-Malina J, Maciejewska A (2013) Uptake of heavy metals by darnel multifloral (Loliummultiflorum Lam.) at diverse soil reaction and organic matter content. Soil Sci Annu 64: 19-23.
    [34] Maciejewska A, Kwiatkowska J (2001) Contamination of soil with heavy metals along highways as well as in the city of Warsaw, In Obieg pierwiastków w przyrodzie (in Polish), Instytut Ochrony Środowiska: 49-54.
    [35] Maciejewska A, Ociepa E, Kwiatkowska J (2003) Changes in the availability of heavy metals and content of organic carbon in polluted soil after the application of an organic-mineral fertilizer obtained from brown coal. Pol J Soil Sci 36:129-136.
    [36] Duda-Chodak A, Błaszczyk U (2008) The impact of nickel on human health. J Elem 13: 685-696.
    [37] Cempel M, Nikel G (2006) Nickel: A Review of Its Sources and Environmental Toxicology. Pol J Environ Stud 15: 375-382.
    [38] Järup L, Åkesson A (2009) Current status of cadmium as an environmental health problem. Toxicol Appl Pharmacol 238: 201-208. doi: 10.1016/j.taap.2009.04.020
    [39] Mohammed AS, Kapri A, Goel R (2011) Heavy Metal Pollution: Source, Impact, and Remedies: Biomanagement of Metal-Contaminated Soils. Environ Pollut 20: 1-28. doi: 10.1007/978-94-007-1914-9_1
    [40] Pistocchi A, Groenwold J, Lahr J, et al. (2011) Mapping Cumulative Environmental Risks: Examples from the EU NoMiracle Project. Environ Model Assess 16: 119-133. doi: 10.1007/s10666-010-9230-6
    [41] Strzelecki R (2004) Directions and problems of abiotic environment research in Poland conducted by the State Geological Survey. Pol Geol Rev 52: 716-719.
    [42] Institute of Soil Science and Plant Cultivation-State Research Institute (2010) The assessment of soil contamination in Poland. Available from http://www.iung.pulawy.pl/index.php?option=com_content&view=article&id=147:zanieczyszczenia&catid=45:oferta&Itemid=115. Accessed 11 June 2016.
  • Reader Comments
  • © 2020 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(6013) PDF downloads(553) Cited by(7)

Figures and Tables

Figures(8)  /  Tables(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog