Next Article in Journal
Recent Advances in Intelligent Data Analysis and Its Applications
Next Article in Special Issue
A Modern Approach to Securing Critical Infrastructure in Energy Transmission Networks: Integration of Cryptographic Mechanisms and Biometric Data
Previous Article in Journal
Grid-Based Non-Uniform Probabilistic Roadmap-Based AGV Path Planning in Narrow Passages and Complex Environments
Previous Article in Special Issue
Fuzzy Controller in the Products Collecting System of the Jig for Minerals Beneficiation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Signal Processing from the Radiation Detector of the Radiometric Density Meter Using the Low-Pass Infinite Impulse Response Filter in the Measurement Path in the Coal Enrichment Process Control System

by
Jarosław Joostberens
1,
Aurelia Rybak
1,*,
Aleksandra Rybak
2 and
Paulina Gwoździk
3
1
Department of Electrical Engineering and Automation in Industry, Faculty of Mining, Safety Engineering and Industrial Automation, Silesian University of Technology, 44-100 Gliwice, Poland
2
Department of Physical Chemistry and Technology of Polymers, Faculty of Chemistry, Silesian University of Technology, 44-100 Gliwice, Poland
3
Silesian University of Technology, 44-100 Gliwice, Poland
*
Author to whom correspondence should be addressed.
Electronics 2024, 13(1), 227; https://doi.org/10.3390/electronics13010227
Submission received: 4 November 2023 / Revised: 21 December 2023 / Accepted: 2 January 2024 / Published: 4 January 2024

Abstract

:
One of the most common coal preparations is enrichment in a jig using a float regulation system. The latest solutions propose to compensate for significant measurement errors of the float by introducing a radiometric density meter operating on the principle of gamma radiation absorption into the bottom product discharge zone of the jig. The signal from the radiometric density meter detector is in the form of a sequence of pulses with a Poisson time distribution, which are counted by a counter, as a form of digital low-pass filter. The requirement to maintain accuracy at an appropriate level forces the measurement time to be extended, which worsens the dynamic properties. Changes in the density of the coal–water medium have an unsteady, cyclical course, resulting from the principle of operation of the jig. The research goal was to develop an algorithm for processing the signal from the radiation detector using an IIR filter in the measurement path in a way that ensures optimization of the dynamic properties of the radiometric density meter operating in the control system of the coal enrichment process in the jig. For this purpose, a low-pass IIR filter was introduced into the measurement path to process the signal from the pulse counter. The identified course of the medium density for one cycle (the first) served as a reference signal. A first-order IIR filter was proposed, with a constant parameter selected on the basis of the reference signal and a parameter depending on the time derivative of the identified density of the medium. The mean squared error MSE was adopted as an indicator for assessing the dynamic properties of the radiometric density meter. The results of simulation tests showed that introducing an IIR filter into the measurement path gives better results in terms of the adopted criterion than using a counter with a constant measurement time. The best results (MSE = 2.05 × 10−4) were obtained using an IIR filter with a parameter that is a linear function of the derivative of the medium density over time, determined for one air pulsation cycle and applied in four subsequent cycles. These results were obtained for an adaptive first-order filter with a variable parameter a from the designated range from 0.833 to 0.999, for a measurement time of 2 ms.

1. Introduction

Radiometric meters have been widely used in the processing of mineral raw materials. In the field of hard coal production, both energy and coke, they are used in the control systems for enriching this raw material, as well as to monitor the quality parameters of coal as a commercial product. In these systems, radiometric meters are used for measuring the quantitative and qualitative parameters of coal, such as density, ash content, concentration of solid parts in transmission for flotation, as well as the density of heavy liquid [1,2,3]. In connection with recent events taking place in Ukraine and Russia, the energy security of the European Union member states, and therefore also of Poland [4,5], was being endangered. Poland is in the privileged situation that it has natural coal resources of 5% of world resources. Coal from Polish deposits is the only one hundred percent guarantee of the country’s energy security. The main assumption of the European Green Deal [6] is the de-carbonization of energy EU systems. However, as of today, this is directly related to the liquidation of the heating and coal-powered power plants [7]. It should be noted that there is a need to ensure a mild transformation of the energy mix, before a source of energy is created which is ready to take over the role of coal in the structure of production and energy consumption in Poland. The solution in this matter is clean coal technologies that will allow the supply of high-quality coal to the combustion process [8,9]. One of the most widespread technologies is enriching coal in a pulsating jig. In this process, the grains are distributed according to their density as a result of cyclical water movement up and down, causing adequate loosening of the material, and then stratification of grains due to differences in the speeds of their falling depending on the density and size. A properly conducted process ensures stratification of the material in such a way that the layers consisting of greater density grains are in the lower part, and the grains with lower density form the upper layers [10,11]. In the automatic control system of the coal enrichment process usually the basic sensor indicating the position of material layer with a selected density is a metal float. In this system, the location of the separation layer is identified by the position of the float and should be stabilized at the level of the overflow threshold using the jig flow rate control. Correct setting of the distribution layer ensures that its upper half together with the material of layers above the distribution layer moves to the upper product above the overflow threshold. In turn, the lower distribution layer along with layers below should go to the lower product [12].
Due to the occurrence of significant measurement errors of the float in automatic control systems, stabilizing the position of the separation layer at the overflow threshold level by the flow rate of the lower product of the jig is a difficult goal to achieve. In order to improve the operation of the system, the latest solutions propose replacing the float or compensating its measurement errors with a radiometric density meter placed at an appropriate height [13,14]. Thanks to this, the signal from the density meter can be used to achieve two main control purposes: stabilizing the course of changes in the density of the medium and stabilizing the separation density. This action will be all the more effective the smaller the deviation of the signal at the density meter output from the actual density value at each moment. In such a situation, a cascade control system can be used with a P controller in the loop with feedback from the float and PI in the main loop from the radiometric density meter, in which the set value is the separation density [14].
Due to the Poisson noise present in the signal from the radiometric density meter, the value of the measured quantity determined at a given moment is subject to error. Therefore, the issue of separating the useful component from the radiometric density meter signal is important when using a radiometric meter to compensate for float measurement errors in automatic control systems of the coal enrichment process in jigs. This shows that achieving the two main goals presented in the automatic control system of the coal enrichment process in the jig using a radiometric density meter requires optimization of the dynamic properties of this meter. In a radiometric density meter, the signal at the output of a gamma ray detector (often in the form of a scintillation counter) is a series of pulses. This sequence is stochastic with a Poisson time distribution, the average value of which is a measure of the analysed density.
One of the basic elements of the measurement path is a counter that counts pulses from the gamma radiation detector in a cyclical manner in a defined unit of time, called the measurement time τs. The value of the counted pulses is converted into the value of the measured quantity using the calibration characteristic. In the measurement path, the counter plays the role of a digital low-pass filter with a parameter in the form of measurement time. The traditional method of filtering the signal from the detector using a counter with a constant but correctly selected measurement time will always be a compromise solution between the measurement accuracy (in the steady state) and the dynamic properties of the radiometric meter. Another solution may be digital filtering of the signal from the gamma radiation detector of a radiometric density meter by supplementing the measurement path with an additional low-pass digital filter. In such a situation, it is possible to shape the dynamic properties of the meter using an appropriately selected low-pass filter connected to the measurement path with a pulse counter. The article describes the method of designing a low-pass IIR filter for use in processing the signal from the detector of a radiometric density meter operating in the control system of the coal enrichment process in a jig. An adaptive form of the IIR filter was also proposed, with the filter parameter selected on the basis of a model of changes in the density of the coal–water medium during the pulsation cycle in the jig, which is a new approach to signal processing from a gamma radiation detector. The results of simulation tests using an adaptive IIR filter in the measurement path for digital signal processing from a radiometric density meter installed at the overflow threshold in the jig are presented.

2. Changes in the Density of the Medium in the Jig during the Air Pulsation Cycle

In a pneumatic jig, medium pulsations are caused by the supplied air using appropriately controlled air valves—inlet and outlet. Automatic control in an open system is then implemented. In the first phase of the air pulsation cycle, when the air inlet valve is opened, the entire bed is lifted with no loosening of grains. A clear increase in the density of the medium is then observed as a result of rising layers of increasing density. After closing the inlet valve, there is a phase of relaxation, as a result of which the grains fall (the density decreases). As a result, they gradually reach a state of constriction, which is accelerated by opening the outlet valve, which makes the density increase visible. The exhaust valve then closes to ensure the same conditions for the next cycle. The cycle duration is from 1 to 1.3 s. The stratified material is separated into two products in the final part of the jig compartment, i.e., in the product discharge zone, which is closed by an overflow threshold.
The cyclical nature of changes in the density of the coal–water medium during a properly conducted coal enrichment process in the jig can be described by a mathematical model. Industrial research carried out at one of the Polish mines has shown that the time course of changes in the density of the medium during the air pulsation cycle can be described sufficiently precisely by a mathematical model [15,16]:
ρ ( t ) A j · e β j · t · s i n ω j · t ψ j + ρ u j
where Aj, βj, ωj, ψj, ρuj—equation parameters (1) for the pulsation cycle.
Values of model parameters (1), determined on the basis of industrial data recorded during five cycles of air pulsation, and corresponding to the times of each cycle separately, standardized from the value of 0 s (start of a given cycle) to the end of the pulsation period (ψj = 0), are listed in Table 1. The identification results are shown in Figure 1.
As can be seen, changes in the density of the coal–water medium are characterized by an unsteady cyclical pattern (except for very short moments of confinement at the end of each cycle), resulting from the principle of operation of the jig.
In the analysis presented below, a simulation model is presented of the measurement method (source, absorber, and scintillation counter) in the form of discrete Poisson noise and controlled average pulse intensity described in [17]. Due to the exponential distribution of random variables in the Poisson process, the time Ti, after which the next impulse appears is determined from the relationship:
T i = l n ( R ) λ
where:
R—random number from 0 to 1;
λ—pulse intensity corresponding to the measured density.
Ti values calculated according to Formula (2) are transformed into subsequent, narrow component time windows and are created on the timeline pulse series with Poisson distribution. This model, verified in laboratory conditions with high accuracy, represents the actual radiometric method of density measurement [15].

3. Signal Processing from the Gamma Radiation Detector of a Radiometric Density Meter

3.1. Radiometric Density Meter

The radiometric densimeter used to measure the density of the medium in the bottom product discharge zone of the jig uses the phenomenon of gamma radiation absorption. The construction of such a density meter is described in [18]. Radiation passing through a specific material (absorber) is weakened, which can be expressed by the formula:
I = I 0 · e μ · x · ρ
where:
I0—radiation intensity without absorber;
μ—mass attenuation coefficient of gamma radiation;
x—absorbent thickness;
ρ —density of the absorbing layer.
The basic elements of an industrial radiometric density meter are a measuring head containing a gamma radiation source, a detector, and an electronics block, the main element of which is a pulse counter used to process the signal from the detector (Figure 2). In the radiometric meter, the radiation source is the 137Cs isotope, and a scintillation counter is used as a gamma radiation detector [19,20,21]. The signal at the output of the radiation detector is a series of pulses s(t). These pulses are counted using a counter (block 2, Figure 2) in a cyclical manner during a time interval called measurement time τs, which can be written:
k i = l = 1 τ s Δ t S τ s Δ t l
where:
S—counted impulse occurring in an elementary time interval;
Δt—elementary pulse counting time determined by the dead time of the radiation detector Δt < τs.
Figure 2. Schematic diagram of a density meter with a pulse counter and a constant measurement time, with the possibility of including an additional filter in the measurement path. Key: 1—measuring head; 2—pulse counter; 3—block for pulse counting to measured value recalculation; 4—additional low-pass digital filter.
Figure 2. Schematic diagram of a density meter with a pulse counter and a constant measurement time, with the possibility of including an additional filter in the measurement path. Key: 1—measuring head; 2—pulse counter; 3—block for pulse counting to measured value recalculation; 4—additional low-pass digital filter.
Electronics 13 00227 g002
The value of the counted pulses is converted into the value of the quantity measured using the calibration characteristic (block 3, Figure 2), which is the dependence of the measured density on the number of counted pulses during the measurement time, usually standardized to one second. The measurement path can be supplemented with an additional low-pass filter (block 4, Figure 2).

3.2. Radiometric Meter with a Pulse Counter and Constant Measurement Time

The traditional approach to processing the signal from the detector is to use a counter with a constant measurement time in the measurement path. In such a situation, the radiometric density meter can be rearranged schematically as in Figure 2, excluding block 4, i.e., the additional digital filter. With respect to a radiometric meter, the Poisson distribution describes the probability of k pulses appearing at the output of a radiation detector in the form of a scintillation counter in a time interval equal to the measurement time τs. The series of pulses from the radiation detector s(t) is always a stochastic signal with a Poisson time distribution, regardless of the measured quantity and the nature of its changes. In the radiometric density meter system, the counter performs cyclical counting of pulses from the detector in a specific time interval, called the measurement time τs, and the obtained values are converted according to the calibration characteristics. In practice, this can be expressed as:
y c i = a 0 a 1 · ln 1 τ s · k [ i ]
where: a0, a1—calibration characteristic coefficients determining the relationship with the parameters of Equation (3), where a0 is proportional to (μx)−1 and a1 = (μx)−1; i—sampling step.
The value of the counted pulses k using the counter (block 3, Figure 2) is a measure of the measured density. The signal from the radiation detector is always burdened with stochastic noise. In the measurement path, the counter functions as a digital low-pass filter with the parameter τs, enabling a reduction in the share of noise in the measurement signal. In the case of a radiometric density meter, the selection of the measurement time comes down to determining a value that will allow obtaining good dynamic properties of the meter without significantly impairing its accuracy. The measurement time should have an optimal value that simultaneously minimizes the dynamic error and maximizes the suppression of fluctuations resulting from the stochastic nature of the signal from the density meter detector.
Analysing the results of identifying the medium density model during the pulsation cycle—presented in Figure 1 and Table 1—it can be seen that the model coefficients change slightly from cycle to cycle in quantitative terms. It can therefore be assumed that the measurement time selected using the identified coal–water density model for the first cycle can be used to filter the signal in subsequent cycles without significant deterioration of the dynamic properties and accuracy of the meter.

3.3. IIR Filter in the Measurement Path of a Radiometric Density Meter

In the case of signal processing only using a pulse counter, the measured values (yc) are maintained for the next measurement time, during which subsequent pulses are counted. Therefore, a step reconstruction of the signal takes place, i.e., the sample value (converted to density) is maintained for a time corresponding to the measurement time τs. This is an analogy to the use of a zero-order extrapolator with a time to maintain the value obtained from the counter equal to the measurement time τs, whose frequency response includes side lobes next to the main lobe. In order to reduce their negative impact, an additional low-pass filter should be included in the measurement path with the pulse counter. One of the possible variants is the use of a low-pass IIR filter due to its satisfactory frequency response and, in particular, its speed.
To design an additional IIR (Infinite Impulse Response) filter included in the measurement path of a radiometric density meter, a method based on the prototype of an analogue filter can be used:
K s = α n s + α n
where:
α—filter parameter;
n—filter order.
The prototype analogue filter with structure (6) was chosen due to its advantageous form, which requires only the selection of the α and the filter order.
A characteristic feature of these filters is that for the first-order filter limit pulsation ω = 1/T, i.e., equal to the filter parameter α. The attenuation increases for this pulsation for higher-order filters. The attenuation in the stopband in the case of a first-order filter is 20 dB/decade. As the order of the filter in Equation (6) increases, the phase delay and attenuation in the stopband increase for the same value of α. This shows that for the same value of α the attenuation in the stopband will be the smallest, with the smallest phase delay for the first-order filter.
The discrete transmittance equation of the n-th order filter, corresponding to Equation (6) of the analogue filter, was determined using the impulse response invariance method. The transmittance function equation in such a situation takes the following form:
K z = α · τ s n 1 · e α · τ s n 1 D n 1 · z z e α · τ s n
In expression (7), Dn—1 is a determinant that is a function of the expression (z e α · τ s ) and with dimensions depending on the filter order n [22]:
D n 1 = 1 1 z e α τ s 0 0 1 2 ! 1 1 z e α τ s 0 1 3 ! 1 2 ! 1 0 1 n ! 1 ( n 1 ) ! 1 ( n 2 ) ! 1
The filter parameters (7) obtained by applying the impulse response invariance method should be scaled and standardized. When using an additional low-pass digital filter to process the signal from the radiometric density meter (Figure 2), the output signal from the pulse counter (after conversion to the value of the measured density) becomes the input signal of the IIR filter, and then the sampling period is equal to the measurement time τs. To ensure good dynamic properties, the measurement time should be short, but its value should also be determined. In such a situation, the order and parameters of the IIR filter must be selected, as well as the value of the measurement time. Also, in this case, the selection of the filter order and parameters based on the medium density model for the first cycle of air pulsation can be used in practice to filter the signal in subsequent cycles.
In the case of a first-order filter, it is preferable to use an equation of the form:
y i = a · y i 1 + 1 a · u [ i ]
where a is the filter parameter.
Parameter a is selected, and in relation to the signal processing from the gamma radiation detector, the value of the measurement time should also be determined.

3.4. A First-Order Filter with a Parameter Depending on the Modulus of the Time Derivative of the Medium’s Density

Analysing the shape of the medium density curve (Figure 1), it can be seen that the cut-off pulsation of the IIR filter should have a large value in moments of sudden, significant density changes, so that the dynamics of changes are not dampened by the filter, and a small value in time intervals characterized by small changes ρ, in order to increase the degree of noise reduction from the measurement signal. This shows that when a filter with Equation (9) is used, the parameter a can be related to the rate of change in the density of the medium. Due to the fact that the empirical model of the density of the medium is known from Equation (1), the derivative of the density over time for a given pulsation cycle is expressed by the equation:
d ρ t d t = A j · e β j · t · ω j · c o s ω j · t ψ j β j · s i n ω j · t ψ j
The use of Equation (10) allows one to avoid difficulties with numerically determining the derivative of a signal burdened with noise. This makes it possible to improve the dynamic properties of the radiometric density meter while maintaining its accuracy at the required level by making the parameter a of the filter Equation (9) dependent on the time derivative of the medium density in such a way that for a high rate of change of the medium density the coefficient a will have a small value, and for a small value of /dt a large value. Knowing the parameters of the model (1) for one cycle, the derivative of density with respect to time (10) in this cycle can be determined. In order to use the derivative of the density of the medium, determined for one cycle, to filter the signal from the density meter in subsequent cycles, it is advantageous to standardize it in the range [0, 1].
d ρ d t N = d ρ ( t ) d t d ρ d t m i n · d ρ ( t ) d t m a x d ρ d t m i n 1
Due to the fact that the filter parameter (9) should depend on the value and not the sign of the medium density derivative with respect to time, this function should be expressed as follows:
a = f d ρ d t N
where aN is the normalized parameter filter value.
With reference to the empirical model (1), the standardized course of the density derivative over time determined for the first cycle is shown in Figure 3.
By identifying the model of the density course of the medium (1) for the first cycle, it is possible to determine the derivative of this quantity with respect to time in this cycle and use its standardized module to filter the signal in subsequent cycles. The proposed method of processing the signal from the detector using a first-order filter with Equation (9) and a parameter depending on the rate of change of medium density during the air pulsation cycle in the jig is shown in Figure 4.
In a properly operating coal enrichment system in the jig, the modulus of the standardized derivative of the medium density determined for the first cycle will be qualitatively similar to the rate of changes in subsequent cycles. As the starting time of each subsequent air pulsation cycle (t0) is known and results from the operation of the control system, the process of changing the filter parameter a is carried out cyclically in the same way. In the proposed algorithm, function (12) was adopted as a linear dependence of the normalized value of the filter parameter with Equation (9) on the normalized modulus of the model of the derivative of the medium density over time, which can be expressed by the formula:
a N = 1 d ρ d t N
The normalized aN value, determined on the basis of Equation (13), is converted into an absolute value (so-called denormalization) using the equation:
a = a m a x a m i n · a N + a m i n
When using a filter (9) with a parameter changing according to the relationship (12), the measurement time τs should be as short as possible, and assuming that amax ≈ 1 (e.g., 0.999), the value of amin should be selected.

4. Simulation Studies of the Use of the IIR Filter to Process the Signal from a Radiometric Density Meter

4.1. Test Conditions

The research task included a comparative analysis of the results of using the classic form of signal processing from the detector, i.e., using a counter with a constant measurement time τs, with the results obtained when the measurement path was supplemented with a low-pass IIR filter.
The research was carried out on the basis of simulated waveforms representing the course of the density of the medium in the jig during five successive cycles of air pulsation. A series of pulses, representing the signal from the gamma radiation detector with changing density of the medium in the jig (Table 1), was generated using the simulation model described in Section 2. Taking into account the elementary pulse counting time Δt equal to 1 μs, a data sequence with a length of 6,350,000 points was obtained, corresponding to the duration of six pulsation cycles, equal to 6.35 s. The lengths of the data series for cycles I to V were, respectively: 1,300,000; 1,250,000; 1,200,000; 1,300,000; 1,300,000. The output signals from the counter yc, the filter y, and the density waveform ρ were observed discretely with a sampling period of Ts = 1 ms. This ensured the possibility of conducting a comparative analysis of the obtained results, regardless of the method used to process the signal from the detector. To convert the number of pulses counted in time τs, Equation (5) was used, with the parameter values used in the industrial device being a0 = 5.25010 and a1 = 0.40206, respectively.
With regard to the processing of the signal from the detector using a counter with a constant measurement time, first, the optimal values of τs in terms of criterion (9) were determined for each pulsation cycle. This has practical justification due to the fact that, in industrial conditions, the starting time of each subsequent air pulsation cycle is known and taken into account in the control algorithm. Then, the value of the measurement time was determined, for which the minimum was reached by expression (9) for the entire data sequence. With respect to cycles I to V, for the shortest and longest measurement times, calculations were carried out to determine the value of the expression (9).
In the second stage of the research, the process of selecting the measurement time and parameters of IIR filters with n orders from 1 to 5 for the first air pulsation cycle was carried out. Also, in this case, a criterion value was used (9). For the case for which the minimum value of expression (9) was achieved, the order and values of the filter parameters were accepted as correct and used to process the signal from the detector in subsequent air pulsation cycles. In the last variant of the conducted simulation tests of signal processing from the gamma radiation detector into the measurement path with a pulse counter, a first-order filter was introduced with a parameter changed depending on the value |/dt|N. In the proposed algorithm, the modulus of the standardized density time derivative identified for the first pulsation cycle was used in subsequent cycles. In such a case, the measurement time should be as short as possible so that the pulse counter does not introduce an additional phase shift. Therefore, τs = 2 ms was assumed. The largest value of the filter parameter (15) amax close to 1, i.e., 0.999, was assumed, and the parameter amin was determined by simulation, changing its value from 0.750 to 0.925. For the correct value of amin, the one for which the value of the indicator (15) was the smallest was adopted.

4.2. Criterion for Selecting the Measurement Time and IIR Filter Parameters

Knowing the mathematical model of changes in the density of the medium for one cycle, it is possible to select the measurement time by minimizing the value of the criterion in the form of a square error expressed by the formula:
M S E j = 1 N j i = 1 N j y j i ρ j [ i ] 2
where:
yc(j)—signal at the output of a density meter with a constant measurement time counter or with a constant measurement time counter and an additional IIR filter connected to the measurement path;
Nj—number of data during the j-th pulsation cycle;
j—cycle number (from I to V).
The indicator (15) can also be used to select the measurement time or assess its correctness for all pulsation cycles. In that case, j stands for cycles I to V.
To select the order and parameters of the IIR filter with Equation (4), a reference signal in the form of model (1) can also be used. In such a situation, the correct order and parameters of the filter can be considered those for which the value of the criterion in the form of the mean square error is minimized (15).

4.3. Results

The results of the first stage of simulation tests, in the form of measurement time values (τs) determined by minimizing criterion (15) in the situation of using a pulse counter with a constant measurement time to process the signal from the detector, are summarized in Table 2. The table also gives the optimal measurement time determined for all pulsation cycles, and the results of using the shortest and longest of the optimal (in the sense of MSE) measurement times. Treating the measurement time determined for the first cycle as the value used in subsequent air pulsation cycles, the values of the MSE index were determined for cycles II-V and for the entire course using this parameter, and the results are also summarized in Table 2.
The course of the signal from the density meter using a constant measurement time of 45 ms is shown in Figure 5.
The assessment of the order of the IIR model with Equation (7) was carried out in relation to density changes during the first air pulsation cycle, and the results are summarized in Table 3. The values of filter parameter α are given, and measurement times for which the minimum value of the mean square error was obtained for a given filter order.
Based on the results in Table 2 and Table 3, the first-order filter was adopted as the most advantageous from the point of view of obtaining good dynamic properties of the radiometric meter while maintaining accuracy at an appropriate level. The values of the parameter a of the filter using Equation (9) determined to be optimal, in the sense of MSE for individual air pulsation cycles, are included in Table 4. Due to the proposed filter parameter selection algorithm, which involves determining the parameter a (as well as the measurement time τs) based on the identified course of the medium density during the first cycle and using this setting for signal processing during subsequent cycles, calculations were carried out aiming at determining the MSE value for signal filtering in pulsation cycles II to V for a = 0.884 and τs = 5 ms. The results are also summarized in Table 4.
The calculation results from Table 4 for the first air pulsation cycle are shown in Figure 6. The value of the mean square error for cycles II to V, using the IIR filter with the parameter calculated for the first air pulsation cycle, is 2.25 × 10−4. The waveform of the signal from the first-order filter for all pulsation cycles is shown in Figure 7.
In the last stage of the research, calculations were carried out to determine the range of variability of parameter a in accordance with the algorithm described in Section 3.4. The results obtained are summarized in Table 5 and presented graphically (Figure 8). However, the waveforms of the signal from the filter with parameter a depending on |/dt| are shown in Figure 9. The value of criterion (15) determined for the course from cycle II to V is 2.05 × 10−4.
The mean square error values obtained for all applied methods of signal processing from the gamma radiation detector of the radiometric density meter are presented graphically in Figure 10.

5. Discussion

The comparative analysis of the results was carried out based on the quantitative assessment of the criterion in the form of the mean square error (Equation (15)). This was possible thanks to the use of the sampling period Ts with an identical value in the observation of the resulting signals (after processing) for each of the considered cases, regardless of the value of the measurement time τs.
The traditional method of processing the signal from a gamma radiation detector, which involves using a counter with a selected measurement time, requires determining its value. In a comparative analysis of the results, they can be treated as reference values. The results in the form of MSE = f(τs) waveforms presented in Table 2 show that the optimal values of measurement time in the sense of MSE for five air pulsation cycles range from 38 ms to 61 ms. It should be noted that as many as three values of τs obtained from the calculation results for five cycles amount to approximately 40 ms. This affects the measurement time determined on the basis of the density curves of the medium of all considered pulsation cycles. It is 45 ms with an MSE value of 3.92 × 10−4. It is true that it is equal to the result obtained for the first cycle, but this similarity should be treated as accidental. As the results in Table 2 show, introducing a measurement time of 45 ms for processing the signal from the detector for all cycles obviously leads to the deterioration of the MSE results, increasing their values by 1.2 to 1.4 times.
Based on the results obtained in relation to the first pulsation cycle (Table 3), it can be concluded that supplementing the measurement path with an IIR filter with Equation (7) is justified only for first and second order filters, because for these algorithms the mean square error is, respectively, 2.47 × 10−4 and 2.85 × 10−4, i.e., their values are lower than when using only a pulse counter with a constant measurement time. Due to the lowest MSE value obtained among the considered IIR filters with Equation (7), a first-order filter was adopted in the further part of the research, using its more convenient form—Equation (9). Comparing the results in Table 3 and Table 4, it should be concluded that supplementing the measurement path with a first-order IIR filter with a parameter selected for each cycle separately leads to significantly lower MSE values than in the case of using a constant measurement time counter. The MSE values in this case range from 60% to 88% of the corresponding criterion values obtained for classical signal processing from the detector. The parameter a of the filter (9) has values ranging from 0.780 to 0.902.
Here again, the determined values of the measurement time (Table 4), as the second parameter in the situation of processing the signal from the detector using a first-order IIR filter, range from 4 to 10 ms and, unfortunately, have different values in each pulsation cycle.
In the case of the proposed algorithm for selecting the filter parameter (and possibly measurement time) only for the first cycle and applying it in subsequent air pulsation cycles, differences in the determined optimal values of parameters a and τs in terms of the mean square error have an unfavourable impact on the value of this criterion (15).
Based on the results from Table 4, it can be concluded that the filter with Equation (9) and parameter a = 0.884 and τs = 5 ms, selected for the first cycle while minimizing MSE and used in subsequent pulsation cycles, gives criterion values from 1.46 × 10−4 (cycle III) to 3.41 × 10−4 (cycle V), which is smaller than when using the pulse counter alone. In this case, the MSE value determined for pulsation cycles II to V is 2.25 × 10−4. This leads to the conclusion that supplementing the measurement path with an additional first-order IIR filter is justified. A further reduction in the MSE criterion was achieved when the value of the parameter a of the filter (9) was dependent on the rate of change in the density of the medium in the first pulsation cycle. By selecting the minimum value of the filter parameter 0.833 from the range of its changes [amin, amax] using the mean square error criterion (as shown in Figure 8b), assuming amax = 0.999, the MSE index values were lower than in the case of classic signal processing with a detector (counter). In this case, the indicator values range from 31% to 74% of the MSE of the reference value of each cycle (Table 2, row 3). Comparing the results of using the IIR filter with Equation (9) with a constant value of parameter a with the algorithm in which a = f|/dt|, it should be stated that lower MSE values were obtained in four out of five cycles when a changed form of parameter a was used. Only a slightly more favourable result was obtained with a constant value of parameter a for cycle IV. Calculated for cycles II to V, i.e., during the period of use of the range of variability of the filter parameter a (9) determined in cycle I, the value of index (15) is 2.05 × 10−4 and is the smallest value. From this point of view, supplementing the measurement path of the radiometric density meter with a first-order IIR filter with parameter a, the value of which depends on the time derivative of the density of the medium, is the best solution.

6. Conclusions

In traditional industrial signal processing systems from a gamma ray detector with a pulse counter (with a constant measurement time), due to the fast-changing component of Poisson noise in the measurement signal of a radiometric density meter, the density values determined in subsequent moments corresponding to the measurement time are burdened with a measurement error, as the greater the value, the shorter the measurement time. Therefore, the issue of separating the useful component from the radiometric density meter signal in the coal enrichment process control system in the jig is important, where changes in the density of the coal–water medium are unsteady and cyclical (except for very short moments of confinement at the end of each cycle), resulting from the principle of operation of this coal enrichment process. In such a case, the radiometric density meter should have good dynamic properties and sufficient accuracy. Based on the research carried out, it can be concluded that the use of the measurement time selected based on the density course of the medium of the first cycle to process the signal from the detector in subsequent cycles is possible, but it only provides intermediate features of a radiometric density meter between the appropriate measurement accuracy and its dynamic properties. The advantage of this solution is its simplicity. The disadvantages include fluctuations in the signal and a phase shift in relation to the density of the medium.
The introduction of an IIR filter designed based on a prototype analogue filter (Equation (6)) into the measurement path of a radiometric density meter using the impulse response invariance method is justified only in relation to the first-order filter. Filters with higher orders (n > 2), despite better attenuation in the stopband, introduce a larger phase shift, resulting in an increase in the dynamic error. This is especially important in the first phase of the pulsation cycle, when a significant increase in the density of the medium is observed. A practical algorithm consisting in identifying the course of the density of the coal–water medium during the first cycle of air pulsation and using this model to select the parameter a of the filter with Equation (9) used in the subsequent four cycles is a better solution than using only a counter with a constant measurement time. Further improvement of the dynamic properties (without deteriorating the accuracy) of the radiometric density meter was achieved by changing the parameter a of the IIR filter from the model of the derivative of the medium density over time in the form of Equation (10), considering the identified model parameters (Table 1) for cycle I. Using the relationship aN = f(|dρ/dt|N) standardized in the range [0, 1] in the form of a linear function, the filter parameter a is related to the rate of change of the medium density at a given moment of the pulsation cycle in the jig. The method of using the model of the derivative of the medium density over time, determined for cycle I and used in the subsequent four cycles to change the parameter a of the IIR filter, turned out to be the best method among those considered.

Author Contributions

Conceptualization, J.J. and A.R. (Aurelia Rybak); methodology, J.J.; software, J.J.; validation J.J.; formal analysis J.J.; investigation, J.J., A.R. (Aurelia Rybak) and A.R. (Aleksandra Rybak); resources A.R. (Aurelia Rybak) and P.G.; data curation J.J.; writing—original draft preparation, J.J., A.R. (Aurelia Rybak), A.R. (Aleksandra Rybak) and P.G.; writing—review and editing J.J., A.R. (Aurelia Rybak) and A.R. (Aleksandra Rybak); visualization, J.J.; supervision, J.J.; project administration J.J. and A.R. (Aurelia Rybak); funding acquisition A.R. (Aleksandra Rybak). All authors have read and agreed to the published version of the manuscript.

Funding

The research leading to these results has received funding from the Norway Grants 2014–2021 via the National Centre for Research and Development. Grant number NOR/SGS/MOHMARER/0284/2020-00.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Cierpisz, S. Computer-based monitoring and control systems in coal preparation plants. In Mineral Processing on the Verge of the 21st Century; Routledge: London, UK, 2017; pp. 385–389. [Google Scholar] [CrossRef]
  2. Clarkson, C.; Hornsby, D.; Walker, D. Automatic Flotation Control using On-Stream Ash Analysis. Coal Prep. 1993, 12, 41–64. [Google Scholar] [CrossRef]
  3. Cierpisz, S.; Kryca, M.; Sobierajski, W. Control of coal separation in a jig using a radiometric meter. Miner. Eng. 2016, 95, 59–65. [Google Scholar] [CrossRef]
  4. Benton, T.G.; Froggatt, A.; Wellesley, L.; Grafham, O.; King, R.; Morisetti, N.; Schröder, P. The Ukraine war and threats to food and energy security. Chatham House Int. Aff. Think Tank 2022, 4, 1–47. [Google Scholar]
  5. Żuk, P. National energy security or acceleration of transition? Energy policy after the war in Ukraine. Joule 2022, 6, 709–712. [Google Scholar] [CrossRef]
  6. Fetting, C. The European green deal. ESDN Rep. 2020, 53. [Google Scholar] [CrossRef]
  7. Burke, A.; Fishel, S. A coal elimination treaty 2030: Fast tracking climate change mitigation, global health and security. Earth Syst. Gov. 2020, 3, 100046. [Google Scholar] [CrossRef]
  8. Yang, Z.; Xing, Y.; Wang, D.; Yangchao, X.; Gui, X. A new process based on a combination of gravity and flotation for the recovery of clean coal from flotation tailings. Energy Sources Part A Recovery Util. Environ. Eff. 2018, 40, 420–426. [Google Scholar] [CrossRef]
  9. Brożek, M.; Surowiak, A. Effect of particle shape on jig separation efficiency. Physicochem. Probl. Miner. Process. 2007, 41, 397–413. [Google Scholar]
  10. Surowiak, A.; Niedoba, T.; Wahman, M.; Hassanzadeh, A. Optimization of Coal Production Based on the Modeling of the Jig Operation. Energies 2023, 16, 1939. [Google Scholar] [CrossRef]
  11. Xie, G.Y.; Wu, L.; Ou, Z.S.; Yu, H.S. Research on fine coal classified flotation process and key technology. Procedia Earth Planet. Sci. 2009, 1, 701–705. [Google Scholar] [CrossRef]
  12. Joostberens, J.; Cierpisz, S. Nonlinear control of coal separation in a jig. IFAC-PapersOnLine 2021, 54, 31–36. [Google Scholar] [CrossRef]
  13. Lyman, G.J. Review of jigging principles and Control. Coal Prep. 1992, 11, 41–72. [Google Scholar] [CrossRef]
  14. Cierpisz, S. Control of a Coal Separation Process in Jigs; Wydawnistwo Politechniki Śląskiej: Gliwice, Poland, 2012. [Google Scholar]
  15. Joostberens, J. Dynamical Properties Optimization of Radiometric Density Meter for Monitoring of a Coal Separation Process in a Jig; Monograph; Silesian University of Technology Publishing House: Gliwice, Poland, 2019. (In Polish) [Google Scholar]
  16. Joostberens, J.; Cierpisz, S. An adaptive radiometric meter with variable measurement time for monitoring of coal jigs operation. IFAC-PapersOnLine 2020, 53, 488–493. [Google Scholar] [CrossRef]
  17. Knuth, E.D. The Art of Computer Programming; Addison-Wesley: Boston, MA, USA, 1998. [Google Scholar]
  18. Cierpisz, S.; Kryca, M.; Gola, M.; Sobierajski, W. OS-C Radiometric Density Meter for a Jig—Principle of Operation. In Proceedings of the Przeróbka Kopalin, 2014 Conference, Szczyrk, Poland, 15–16 May 2014; pp. 41–52. (In Polish). [Google Scholar]
  19. Cierpisz, S.; Kryca, M.; Gola, M.; Sobierajski, W. Radiometric Methods of Monitoring of a Coal Separation Process in A Pulsating Jig. Miner. Resour. Manag. 2016, 2, 125–134. [Google Scholar] [CrossRef]
  20. Knoll, G.F. Radiation Detection and Measurement; John Wiley & Sons: New York, NY, USA, 2000. [Google Scholar]
  21. Bonczyk, M.; Wysocka, M.; Bukowska, M.; Michali, B.; Prusek, S.; Wadas, M. Radiometric Method for the Research on Geomechanical Parameters of Rocks. Pure Appl. Geophys. 2017, 174, 1025–1031. [Google Scholar] [CrossRef]
  22. Kuzin, L.T. Analysis and Synthesis of Discrete Control Systems; Maszgiz: Moscow, Russia, 1962. [Google Scholar]
Figure 1. Changes in the medium density during five cycles of air pulsation in the jig. Key: 1—model; 2—measurement data; I–V—air pulsation cycles.
Figure 1. Changes in the medium density during five cycles of air pulsation in the jig. Key: 1—model; 2—measurement data; I–V—air pulsation cycles.
Electronics 13 00227 g001
Figure 3. Derivative of the medium density over time determined on the basis of model (1) for the first pulsation cycle, standardized in the range [0, 1].
Figure 3. Derivative of the medium density over time determined on the basis of model (1) for the first pulsation cycle, standardized in the range [0, 1].
Electronics 13 00227 g003
Figure 4. Processing the signal from the gamma radiation detector using a filter with a parameter depending on the modulus of the time derivative of the medium density identified for the first cycle and used in subsequent cycles.
Figure 4. Processing the signal from the gamma radiation detector using a filter with a parameter depending on the modulus of the time derivative of the medium density identified for the first cycle and used in subsequent cycles.
Electronics 13 00227 g004
Figure 5. Processing of the signal from the gamma radiation detector using a filter with a parameter depending on the modulus of the time derivative of the medium density identified for the first cycle and used in subsequent cycles.
Figure 5. Processing of the signal from the gamma radiation detector using a filter with a parameter depending on the modulus of the time derivative of the medium density identified for the first cycle and used in subsequent cycles.
Electronics 13 00227 g005
Figure 6. Results of determining the optimal value of the amin parameter of a first-order filter: (a) waveform of the signal at the filter output against the background of the medium density during the first cycle of air pulsation; (b) relationship between MSE and the a parameter and the τs measurement time.
Figure 6. Results of determining the optimal value of the amin parameter of a first-order filter: (a) waveform of the signal at the filter output against the background of the medium density during the first cycle of air pulsation; (b) relationship between MSE and the a parameter and the τs measurement time.
Electronics 13 00227 g006
Figure 7. Signal waveforms from the IIR filter with a = 0.884 and a measurement time of 5 ms against the background of medium density waveforms for five air pulsation cycles.
Figure 7. Signal waveforms from the IIR filter with a = 0.884 and a measurement time of 5 ms against the background of medium density waveforms for five air pulsation cycles.
Electronics 13 00227 g007
Figure 8. Results of determining the optimal value of the amin parameter of a first-order filter with a parameter changed as a function of the derivative of the medium density over time: (a) waveform of the signal at the filter output against the background of the medium density during the first cycle of air pulsation; (b) dependence of MSE on the amin parameter.
Figure 8. Results of determining the optimal value of the amin parameter of a first-order filter with a parameter changed as a function of the derivative of the medium density over time: (a) waveform of the signal at the filter output against the background of the medium density during the first cycle of air pulsation; (b) dependence of MSE on the amin parameter.
Electronics 13 00227 g008
Figure 9. Signal processing from the gamma radiation detector using a filter with a parameter depending on the modulus of the time derivative of the medium density identified for the first cycle and used in subsequent cycles.
Figure 9. Signal processing from the gamma radiation detector using a filter with a parameter depending on the modulus of the time derivative of the medium density identified for the first cycle and used in subsequent cycles.
Electronics 13 00227 g009
Figure 10. The results of signal processing from the gamma radiation detector of the radiometric density meter used in the enrichment process control system in the jig, obtained for various filtration methods. Key: 1—constant measurement time determined for each cycle of air pulsation separately—reference values; 2—constant measurement time determined on the basis of the first cycle; 3—first-order IIR filter with equation (9) with a parameter selected on the basis of ρ(t) of the first cycle; 4—first-order IIR filter with a parameter value resulting from the relationship: a = f(|/dt|) determined for the first air pulsation cycle.
Figure 10. The results of signal processing from the gamma radiation detector of the radiometric density meter used in the enrichment process control system in the jig, obtained for various filtration methods. Key: 1—constant measurement time determined for each cycle of air pulsation separately—reference values; 2—constant measurement time determined on the basis of the first cycle; 3—first-order IIR filter with equation (9) with a parameter selected on the basis of ρ(t) of the first cycle; 4—first-order IIR filter with a parameter value resulting from the relationship: a = f(|/dt|) determined for the first air pulsation cycle.
Electronics 13 00227 g010
Table 1. Values of the model of the medium density changes parameters in the jig, determined on the basis of industrial data recorded during ten cycles of air pulsation [15].
Table 1. Values of the model of the medium density changes parameters in the jig, determined on the basis of industrial data recorded during ten cycles of air pulsation [15].
CycleAβωρut
I0.17112.525.721.3830from 0 to 1.30 s
II0.20883.746.831.3978from 0 to 1.25 s
III0.16813.877.001.4004from 0 to 1.20 s
IV0.14762.055.641.4034from 0 to 1.30 s
V0.20164.216.661.4018from 0 to 1.30 s
Mean
value
0.17943.286.371.3973-
Table 2. Values of optimal measurement times in terms of criterion (15) for individual pulsation cycles and the entire course of the medium density and MSE values, including the shortest and longest measurement times and the measurement time determined for cycle I.
Table 2. Values of optimal measurement times in terms of criterion (15) for individual pulsation cycles and the entire course of the medium density and MSE values, including the shortest and longest measurement times and the measurement time determined for cycle I.
CyclesIIIIIIIVVfrom I to V
τs, (ms)454161384145
MSE 10−43.183.332.242.753.493.92
τs, (ms)45
MSE 10−43.184.533.223.314.813.92
τs, (ms)38
MSE 10−44.044.613.152.756.024.36
τs, (ms)61
MSE 10−45.304.402.243.565.364.89
Table 3. IIR filter parameters of Equation (7) of orders one to five, determined by minimizing criterion (15) determined for the first air pulsation cycle.
Table 3. IIR filter parameters of Equation (7) of orders one to five, determined by minimizing criterion (15) determined for the first air pulsation cycle.
n12345
τs, (ms)555510
1/α (ms)9.003.602.401.801.22
MSE × 10−42.472.853.643.864.54
Table 4. The values of the filter parameter with Equation (9) determined by minimizing the MSE criterion for all air pulsation cycles and the value of the indicator (15) using the parameter a determined for the first cycle and used to process the signal for the remaining pulsation cycles.
Table 4. The values of the filter parameter with Equation (9) determined by minimizing the MSE criterion for all air pulsation cycles and the value of the indicator (15) using the parameter a determined for the first cycle and used to process the signal for the remaining pulsation cycles.
CycleIIIIIIIVV
τs, (ms)589410
a0.8840.7980.8280.9020.780
MSE × 10−42.472.701.391.653.07
τs, (ms)5
a0.884
MSE × 10−42.472.821.461.673.41
Table 5. Values of the mean square error obtained using the changed parameter a of the filter in the range calculated for the first pulsation cycle [amin, amax] = [0.833, 0.999] and used in processing the signal from the detector for the remaining cycles (τs = 2 ms).
Table 5. Values of the mean square error obtained using the changed parameter a of the filter in the range calculated for the first pulsation cycle [amin, amax] = [0.833, 0.999] and used in processing the signal from the detector for the remaining cycles (τs = 2 ms).
CycleIIIIIIIVV
MSE × 10−41.941.370.701.862.57
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Joostberens, J.; Rybak, A.; Rybak, A.; Gwoździk, P. Signal Processing from the Radiation Detector of the Radiometric Density Meter Using the Low-Pass Infinite Impulse Response Filter in the Measurement Path in the Coal Enrichment Process Control System. Electronics 2024, 13, 227. https://doi.org/10.3390/electronics13010227

AMA Style

Joostberens J, Rybak A, Rybak A, Gwoździk P. Signal Processing from the Radiation Detector of the Radiometric Density Meter Using the Low-Pass Infinite Impulse Response Filter in the Measurement Path in the Coal Enrichment Process Control System. Electronics. 2024; 13(1):227. https://doi.org/10.3390/electronics13010227

Chicago/Turabian Style

Joostberens, Jarosław, Aurelia Rybak, Aleksandra Rybak, and Paulina Gwoździk. 2024. "Signal Processing from the Radiation Detector of the Radiometric Density Meter Using the Low-Pass Infinite Impulse Response Filter in the Measurement Path in the Coal Enrichment Process Control System" Electronics 13, no. 1: 227. https://doi.org/10.3390/electronics13010227

APA Style

Joostberens, J., Rybak, A., Rybak, A., & Gwoździk, P. (2024). Signal Processing from the Radiation Detector of the Radiometric Density Meter Using the Low-Pass Infinite Impulse Response Filter in the Measurement Path in the Coal Enrichment Process Control System. Electronics, 13(1), 227. https://doi.org/10.3390/electronics13010227

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

Article Metrics

Back to TopTop