Analysis of a Dataset for Modeling a Transport Conveyor Oleh Pihnastyi 1 and Anna Burduk 2 1 National Technical University "Kharkiv Polytechnic Institute", 2 Kyrpychova, Kharkiv, 61002, Ukraine 2 Wroclaw University of Science and Technology, 27 W. Wyspianskiego, Wrocław, 50370, Poland Abstract The analysis of the works, which considered the use of neural networks for modeling a multi- section transport conveyor, was carried out. The prospects for the use of neural networks for the design of highly efficient control systems for the flow parameters of a multi-section transport conveyor are studied. The problem that limits the use of neural networks for building control systems for the flow parameters of a multi-section transport conveyor is considered. The possibility of constructing generators for generating a data set for the process of training a neural network is being studied. A method for generating a data set based on experimentally obtained measurements of the instantaneous values of the input material flow as a result of the operation of industrial transport systems is proposed. Using dimensionless variables, a statistical analysis of a stochastic flow of material entering the input of the transport system was performed. An estimate of the correlation time of a stochastic process characterizing the input flow of material is given. The recommendations on choosing the type of correlation function for the model of the input material flow were confirmed. It is demonstrated that the input flow of material is a non-stationary stochastic process. Approximations for modeling the input flow of materials of the operating transport system are considered. Keywords 1 transport conveyor, neural network, non-stationary stochastic process, dataset generator 1. Introduction Conveyor-type transport systems are widely used in the mining industry due to the low cost of material transportation [1]. The unit cost of transportation is significantly affected by the length of the transport route and the material load factor of the transport system. Within each section, to reduce the specific energy consumption, belt speed control systems [2, 3, 4], material flow control systems [5, 6] coming from the bunker to the section input, as well as control systems based on the energy management methodology [7, 8]. One of the types of models that are used to design control systems for the flow parameters of the transport conveyor and diagnose the state of the transport system are models based on neural networks. To diagnose the wear level of a conveyor belt, a neural network with the 13-5-1 architecture (13 nodes in the input layer, 4 nodes in the hidden layer, and one node in the output layer) was proposed [9]. The conveyor belt speed control algorithm is represented by a model based on a neural network with 3- 4-3 [10] and 3-10-1 [11] architectures. Transport system state models using a neural network are presented in [12] (architecture 4-9-14) and [13] (architecture 4-20-1). The control system of the flow parameters of a multi-section transport contains a model of a conveyor section using a neural network with the 9-3-2 architecture [14]. It is shown that for the case of the non-stochastic input flow of material for the transport conveyor, the control algorithm of the input flow value has satisfactory accuracy. ITTAP’2022: 2nd International Workshop on Information Technologies: Theoretical and Applied Problems, November 22–24, 2022, Ternopil, Ukraine EMAIL: pihnastyi@gmail.com (A. 1); anna.burduk@pwr.edu.pl (A. 2) ORCID: 0000-0002-5424-9843 (A. 1); 0000-0003-2181-4380 (A. 2) ©️ 2021 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). CEUR Workshop Proceedings (CEUR-WS.org) A comparative analysis of the computational performance of the models is carried out, followed by the conclusion that for transport systems with a large number of individual sections [15, 16, 17, 18], it is advisable to use a neural network when building models of the transport conveyor. It is multi-section transport conveyors with a large number of sections (several tens or hundreds of sections) that open up prospects for the widespread use of models based on neural networks. Additionally, it should be pointed out that models using a neural network scale well. 2. Problem statement One of the problems that limit the use of neural networks for building control systems for the flow parameters of a multi-section transport pipeline is the lack of a sufficient data set for training a neural network. This circumstance is due, on the one hand, to the fact that under production conditions the transport conveyor operates in a narrow range of flow parameters, which is determined by the technological features of the material extraction process, on the other hand, the different structure of the technological route, and, accordingly, the location and connection of the conveyor sections. It is assumed that the flow of material (t ) , incoming into the transport system is subject to stochasticity and can be described as a random process ( t )   d ( t )   s ( t ) , (1) where  d (t ) is a deterministic function of time t ;  s (t ) is stationary centered ergodic process, which is determined by the one-dimensional distribution density [18, 19, 20] 1 2  f  ( )  exp  , 2  (2)  2   2  with mathematical average m  0 , standard deviation  , and correlation function [20] K ()   exp  . 2 (3) The correlation time kor 1  , (4) kor   (ti  t j ) . depends on the method of organizing the production process, determines the characteristic time of the process, at which the correlation between the sections of the random process ti and t j can be neglected. Having performed experimental measurements of the values of the input (t ) and output 1 (t , Sd ) flows of material from sections of the transport conveyor with a length of Sd , under the conditions of production activity, it is possible to form a data set for training a neural network, , on the basis of that a model of the transport system is built. However, as already emphasized above, due to technological limitations on the values of the belt speed a (t ) of the conveyor section and the input flow (t ) , as well as significant differences between the length and branching of transport routes, the data set, constructed in this way for training the neural network, will be very useful for optimizing the control system of the existing transport conveyor. This method is not suitable for designing a control system for the flow parameters of a new transport conveyor with a completely different material transportation route, both in length and in branching. Another way to build an exact correspondence between the flow parameters (t ) and 1 (t, Sd ) is to use the analytical model of the conveyor section [21, 22], which allows you to calculate the value of the output flow 1 (t, Sd ) from the known value of the input flow and the law of belt speed change. In [23] ], a method for generating a data set for training a neural network on the example of a branched eight- section conveyor with a deterministic material flow for the input sections of the transport system (t )   d (t ) ~  d 0   d 1 sin t    , (5) and conveyor belt speed a(t ) ~ a0  a1 sin a t  a  , (6) with arbitrarily given values of coefficients  d 0 ,  d 1 , a 0 , a1 , angular velocities  , a and angles   ,  a . Further development of the research presented in [23] is to build a random value generator (t ) , based on the experimental data of the instantaneous values of the input material flow. This generator can be used to form the required number of sets of input flows with specified properties (2), (3) with subsequent calculations of the output flow values using an analytical model of the conveyor section. This work is devoted to the first step in this study, namely, the use of a set of experimental data to determine the form of a deterministic function of time  d (t ) for a given stationary centered ergodic process  s (t ) , determined by a one-dimensional distribution density (2) and a correlation function (3). 3. Main material Experimental data characterizing the dynamics of the flow of material incoming the transport system input are presented in [24–30] and can be used to construct random process generators (t ) (1). To determine the form of the deterministic function of time  d (t ) let's use the results of [29, 31]. After scanning a graphic image, a set of points  i , ti  for time values is obtained ti  tmin  it  , i  0..N , (7) (t  t ) t  max min , ( N  1) where tmin , tmax are initial and final values of time; N  11591is number of sections  i , ti  used when scanning a graphic image. The tabular data set  i , ti  obtained as a result of scanning is shown in Fig.1. Figure 1: Actual volume capacity for the SRs 2000 bucket-wheel excavator recorded in the Belchatow surface lignite mine during the excavation of one terrace [29] Let us introduce dimensionless variables  (t )  ( )  , (8)   s (t )  s ( )  ,   (t )  d ( )  d ,  where t  t min  , (9) t max  t min taking into account which the flow of material (t ) can be represented in a dimensionless form ( )   d ( )   s ( ) , (10) where  s () is a stationary centered ergodic process, which is determined by the one-dimensional distribution density 1  2  f  ()  exp   , (11) 2  2 with mathematical average m  0 , standard deviation   1 , and correlation function    K  ()  exp   ,   kor  (12) kor  kor  , t max  t min  ( ti  t j )    ( i   j ) . (13) (tmax  tmin ) (t max  tmin ) Taking into account the introduced dimensionless variables, the implementation  () takes the form shown in Fig. 2 Figure 2: Implementation of a stochastic process For a stochastic process  () , known on the time interval [0,1] let's look for a deterministic time function  d () in the form of a Fourier series expansion of a function with period T=1 a0   2    2   d ( )    an cos n    bn sin  n  , (14) 2 n 1 T  n 1 T  1  a cos(2n)d  0 , 0 n (15) 1  b sin(2n)d  0 . 0 n The expansion coefficient a0 is determined from the condition of ergodicity of a stationary centered process  s ()  1 k k   m  lim  s ( )d  0 , (16) k 0 whence, taking into account (15), it follows 1 1 a0 a0 0  (  )d  0 2 d  2 , (17) Let's write the expression for the correlation function  1 k K d (, an )  lim  k    ( )   d ()(   )   d (   )d  k 0 1  1 1  ( )  (   )d    ( ) d (   )   d (   ) d    d ( )  d (   )d , 1 1   0  (18) 0 0 Let's make preliminary simplifications   d (   )   d (   )  a0   bn sin 2n(   )   sin 2n(   )  a0  n 1   2 cos2nan cos2n  bn sin 2n . (19) n 1 1 2  1  () d (  )   d (  )d   2 an cos2n  ( ) cos2nd  a0 0 2 n 1 0  1  2 bn cos2n  ( ) sin 2nd . (20) n 1 0 a0   d (   )    an cos2ncos2n  sin 2nsin 2n  2 n 1    bn sin 2ncos2n  cos2nsin 2n . (21) n 1 By virtue of the fulfillment of the identity (15) 1 2 a0 a0 0 2 d (   ) d  4 , (22) 1 2  an cos2n d (  )d  cos2n  n n sin 2n , an ab (23) 0 2 2 1 2 sin 2n  n cos2n . an bn b 0  bn sin(2n) d (  )d   2 2 (24) Let's substitute the prepared expressions (20), (22)–(24) into (18) and obtain the form of the correlation function   K d (, an )  A0   An cos2n   Bn sin 2n , (25) n 1 n 1 where 1  2 1 a A0   1  0  ( )  (   )d  0 , 4  1 1 a  bn  2 2 An  2 an   ( ) cos2nd  bn   ( ) sin 2nd  n ,   0 0 4  Bn  an bn . The coefficients a n are determined from the equality equating the coefficients for the same expansion harmonics K  ()  K d (, an ) . (26) where K  () is the given function (12) represented by the Fourier series. 4. Analysis of results Consider the form of the correlation function in the zeroth approximation 1  2 1 a0 1   0 K d (, a0 )  A0   (  )  (   ) d  . (27) 4 The distribution density f s ( ) of a random variable  s next to the distribution density f  ( ) (11) is shown in Fig.3. To construct the diagram, 50 intervals were selected in accordance with the recommendations [32, 33] with the number of elements in the histogram N  11591. Figure 3: The density of distribution f s ( ) of a random variable For the process (t ) let's define the values t max 1  (t )dt  ms  4663 t max  t min tmin (28) tmax  1743.76 , tmin  1685.825 t max 1  (t )  ms 2 dt  s 2  14702  (29) t max  t min tmin whence the form for a random process  () in the zeroth approximation can be represented by the equation a 4663  d ( )  0   3.2  (30) 2 1470 ( )  3.2   s ( )  An analysis of the distribution histogram presented in Fig. 3 shows that the zero approximation is not enough to simulate a random input material flow [29]. To describe the non-stationarity of a random process [29], additional terms of expansion (25) are needed. For the zero approximation, the correlation function K d (, a0 ) (27) of the random process can be calculated (30). The chart of the function K d (, a0 ) is shown in Fig.4. Figure 4: Correlation functions K d (, a0 ) and K  () For comparative analysis, a monotonically decreasing function K  () has been added to the Fig.4    K  ()  exp    (31)   kor  with correlation time  kor =1/20. Analysis of the results in Fig. 4 shows that the random process, as expected, has a fairly short correlation time. Adding non-stationary terms to series (14) will make it possible to smooth the correlation functions K d (, a0 ) . 5. Conclusion The quality of the neural network training process, and accordingly the accuracy of the multi-section pipeline model based on the neural network, directly depends on the data set used for training. In this regard, the formation of a data set for training a neural network is an important practical and theoretical task, the solution of which is directly related to the development of the theory of control systems for the flow parameters of a multi-section conveyor. As a method for solving this problem, this paper considers the construction of generators of a random process that simulates the input flow at the incoming to the transport system. The solution to this problem is complicated by the fact that for the existing transport systems the incoming material flow is characterized not only by stochasticity but is also a non-stationary flow. Thus, the method of constructing an input material flow generator based on a set of experimental data will allow taking into account these features in close relationship with the technological features of production. Generators of this kind will make it possible to form the required number of data sets for high-quality training of neural networks in multi-sectional conveyor models. In addition, one and the same input data set can be converted into different output data sets depending on the structure of the transport route and the law of interaction between sections of the transport conveyor. In this paper, the implementation of a random process is analyzed, which is supposed to be used to construct a generator of the input material flow. The stochastic process is represented by the canonical expansion of the non-stationary part in a Fourier series. A statistical analysis of the implementation of the input material flow is carried out and a model of the input material flow in the zero approximation is presented. The recommendation to use the correlation function in exponential form for modeling the input material flow has been confirmed. The prospect of further research is to improve the accuracy of approximation of the statistical characteristics of the generator to the input material flow for the available sets of experimental data. An additional task is to determine the law of distribution of the value of the input cargo flow of material entering the input of the transport system per unit time. 6. References [1] M. Alspaugh, “Latest developments in belt conveyor technology,” MINExpo 2004, New York, USA, 2004. [2] M. Alspauch, “The evolution of intermediate driven belt conveyor,” Bulk Solids Handling, vol. 23(3), pp. 168–173, 2003. http://www.overlandconveyor.com/pdf/bsh_AlspaughM_3_2003.pdf [3] I. Halepoto, M Shaikh, and B. Chowdhry, “Design and Implementation of Intelligent Energy Efficient Conveyor System Model Based on Variable Speed Drive Control and Physical Modeling, “ Journal of Control and Automation, vol. 9(6), pp. 379–388, 2016. [4] O. Pihnastyi, V. Khodusov, “Optimal Control Problem for a Conveyor-Type Production Line,” Cybern. Syst. Anal., vol. 54(5), pp. 744–753, 2018. https://doi.org/10.1007/s10559-018-0076-2 [5] O. Pihnastyi, V. Khodusov, “The optimal control problem for output material flow on conveyor belt with input accumulating bunker,” Bulletin of the South Ural State University. Ser.Mathematical Modelling, Programming & Computer Software, vol. 12(2), pp. 67–81, 2019. http://dx.doi.org/10.14529/mmp190206 [6] P. Bardzinski, P. Walker, and W. Kawalec, “Simulation of random tagged ore flow through the bunker in a belt conveying system,” International Journal of Simulation Modelling, vol. 4, pp. 597– 608 (2018). https://doi.org/10.2507/IJSIMM17(4)445 [7] T. Mathaba, X. Xia, and J. Zhang, “Optimal scheduling of conveyor belt systems under Critical Peak Pricing,” 10th International Power & Energy Conference (IPEC), pp. 315–320, 2012. https://doi.org/10.1109/ASSCC.2012.6523285 [8] T. Cousins, “Using Time of Use (TOU) Tariffs in Industrial, Commercial and Residential Applications Effectively,” TLC Engineering Solutions, pp. 1–15, 2010. [9] A. Kirjanow, “The possibility for adopting an artificial neural network model in the diagnostics of conveyor belt splices,” Interdisciplinary issues in mining and geology, vol. 6, , pp. 1–11, 2016. [10] Xi. Pingyuan, S. Yandong, “Application Research on BP Neural Network PID Control of the Belt Conveyor,” JDIM, vol. 9(6), pp. 266–270, 2011. [11] I. Shareef, H. Hussein, “Implementation of artificial neural network to achieve speed Control and Power Saving of a Belt Conveyor System,” Eastern-European Journal of enterprise technologies, vol. 2(110), pp. 44–53, 2021. https://doi.org/10.15587/1729-4061.2021.224137 [12] Y. Yuan, W. Meng, and X. Sun, “Research of fault diagnosis of belt conveyor based on fuzzy neural network,” The Open Mechanical Engineering Journal, vol. 8, pp. 916–921, 2014. https://doi.org/10.2174/1874155X01408010916 [13] S. Abdollahpor, A. Mahmoudi, and A. Mirzazadeh, “Artificial neural network prediction model for material threshing in combine harvester,” Elixir Agriculture, vol. 52, pp. 11621–1162, 2012. https://www.elixirpublishers.com/articles/1353477517_52%20(2012)%2011621-11626.pdf [14] O. Pihnastyi, G. Kozhevnikov, “Control of a Conveyor Based on a Neural Network,” International Conference on Problems of Infocommunications. Science and Technology, pp. 295– 300, 2020. https://doi.org/10.1109/PICST51311.2020.9467934 [15] M. Koman, Z. Laska, “The constructional solution of conveyor system for reverse and bifurcation of the ore flow, Rudna mine KGHM Polska Miedź SA,” CUPRUM, vol. 3(72), pp. 69– 82, 2014. [16] M. Bajda, R. Błażej, and L. Jurdziak, “Analysis of changes in the length of belt sections and the number of splices in the belt loops on conveyors in an underground mine,” Engineering Failure Analysis, vol. 101, pp. 439–446, 2019. https://doi.org/10.1016/j.engfailanal.2019.04.003 [17] R. Król, W. Kawalec, L. Gładysiewicz, “An effective belt conveyor for underground ore transportation systems,” Earth and Environmental Science, vol. 95(4), pp. 1–9, 2017. https://doi.org/10.1088%2F1755-1315%2F95%2F4%2F042047 [18] V. Prokuda, ”Research and assessment of cargo flows on the main conveyor transport “PSP Pavlogradskaya Mine, DTEK Pavlogradugol,” Mining electromechanics, vol. 288, pp. 107–111, 2012. [19] V. Zaika, Yu. Razumny, V. Prokuda, ”Regulated drives influence on coal flow and energy efficiency of mine conveyor transport system,” Naukovyi Visnyk Natsionalnoho Hirnychoho Universytetu, vol. 3, pp. 82–88, 2015. [20] V. Kondrakhin, N. Stadnik, P. Belitsky, ”Statistical Analysis of Mine Belt Conveyor Operating Parameters,” Naukovi pratsi DonNTU, vol. 2(26), pp. 140–150, 2013. [21] O. Pihnastyi, V. Khodusov, ”Model of a composite magistral conveyor line,” IEEE International Conference on System analysis & Intelligent computing, Ukraine, pp. 68–72, 2018. https://doi.org/10.1109/saic.2018.8516739 [22] O. Pihnastyi, V. Khodusov, ”Calculation of the parameters of the composite conveyor line with a constant speed of movement of subjects of labour,” Scientific bulletin of National Mining University, vol. 4(166), pp. 138–146, 2018 –https://doi.org/10.29202/nvngu/2018-4/18 [23] O. Pihnastyi, G. Kozhevnikov and T. Bondarenko, “An Analytical Method for Generating a Data Set for a Neural Model of a Conveyor Line,” IEEE 11th International Conference on Dependable Systems, Services and Technologies, Ukraine, pp. 202–206, 2020. https://doi.org/10.1109/DESSERT50317.2020.9125041 [24] B. Jeftenić, L. Ristić, M. Bebić, S. Štatkić, D. Jevtić, I. Mihailović, and N. Rašić, “Realization of system of belt conveyors operation with remote control,” Structural integrity and life, vol. 10(1), pp. 21–30, 2010. http://divk.inovacionicentar.rs/ivk/ivk10/021-030-IVK1-2010-BJ-LR-MB-SS- DJ-IM-NR.pdf [25] B. Doroszuk, R. Król, and J. Wajs, “Simple design solution for harsh operating conditions: redesign of conveyor transfer station with reverse engineering and DEM simulations,” Energies, vol. 14, pp. 1–13, 2021. https://doi.org/10.3390/en14134008 [26] F. Zeng, C. Yan, Q. Wu, T. Wang, “Dynamic behaviour of a conveyor belt considering non- uniform bulk material distribution for speed control,” Applied sciences, vol. 10, pp. 1–13, 2020. https://doi.org/10.3390/app10134436 [27] N. Stadnik, “Frequency-Controlled Electric Drive of Band Conveyors Based on Self- Ventilating Engines,” Scientific Bulletin of the Donetsk National Technical University, vol. 2, pp. 226–232, 2012. [28] A. Semenchenko, M. Stadnik, P. Belitsky, D. Semenchenko, and O. Stepanenko, “The impact of an uneven loading of a belt conveyor on the loading of drive motors and energy consumption in transportation,” Eastern-European Journal of Enterprise Technologies, vol. 82, pp. 42–51, 2016. https://doi.org/10.15587/17294061.2016.75936 [29] W. Kawalec, R. Król, “Generating of Electric Energy by a Declined Overburden Conveyor in a Continuous Surface Mine,” Energies, vol. 14, pp. 1–13, 2021. https://doi.org/10.3390/en14134030 [30] V. Kondrakhin, N. Stadnik, P. Belitsky, “Traffic measurement on the belt conveyor through removable tensor device with belt tension accounting,” DonNTU scientific works, vol. 1, pp. 79– 87, 2013. [31] L. Gladysiewicz, W. Kawalec, “Optimised Selection Of A Belt Conveyor Loaded By A BWE,” 8th InternationalSymposium of Continuous Surface Mining ISCSM, Aachen, Germany, pp. 24–27, September 2006. [32] A. Tynynyka, “Application of the entropic coefficient for interval number optimization during interval assessment,” Technology and design in electronic equipment, vol. 3, pp. 49–54, 2017. https://www.readcube.com/articles/10.15222%2Ftkea2017.3.49 [33] I. Heinhold, K. Gaede, Ingeniur statistic.Springler Verlag, 1964.Wang, Xin, Tapani Ahonen, and Jari Nurmi. "Applying CDMA technique to network-on-chip." IEEE transactions on very large scale integration (VLSI) systems 15.10 (2007): 1091-1100.