Research Article  Open Access
Fuyan Guo, Yanchao Zhang, Yue Wang, Ping Wang, Peijun Ren, Rui Guo, XinYi Wang, "Fault Detection of Reciprocating Compressor Valve Based on OneDimensional Convolutional Neural Network", Mathematical Problems in Engineering, vol. 2020, Article ID 8058723, 10 pages, 2020. https://doi.org/10.1155/2020/8058723
Fault Detection of Reciprocating Compressor Valve Based on OneDimensional Convolutional Neural Network
Abstract
Reciprocating compressors are important equipment in oil and gas industries which closely relate with the healthy development of the enterprise. It is essential to detect the valve fault because valve failures account for 60% in total failures. For this field, an artificial neural network (ANN) is widely used, but a complex network is not suitable for its low accuracy and easy overfitting. This paper proposes a fault diagnosis model of a reciprocating compressor valve based on a onedimensional convolutional neural network (1DCNN). This method takes the differential pressure and differential temperature of each compressor stage as the input of 1DCNN, using the characteristics of the CNN to extract the features and finally using Softmax to classify the fault. In order to verify this method, it is compared with LMBP, RBF, and BP neural networks. The results show that the fault recognition rate of 1DCNN reaches 100%, which proves the effectiveness and feasibility of the proposed method.
1. Introduction
Reciprocating compressor is one of the most widely used compressor technologies in today’s oil and gas industries. It can compress a variety of gases and has a wide range of applications, high compression efficiency, and stable working pressure. However, due to its complicated structure and many vulnerable parts, once a failure cannot be detected and eliminated in time, it will bring huge losses to the enterprise [1–7]. The literature indicates that 60% of reciprocating compressor failures are valve failures, and the number of shutdowns caused by valve failures accounted for 36% and accounted for 50% of the total maintenance costs [8–10]. Therefore, monitoring the failure of the reciprocating compressor valve can not only reduce the overall maintenance cost but also improve the stability of the compressor operation.
Due to the characteristic of deep learning [11] which can overcome the shortcomings of traditional learning model networks such as complexity, low accuracy, and easy overfitting, it has quickly become a hot topic in the research field in the past few years. It has made breakthrough progress in image recognition, voice recognition, language processing, and other fields. In view of the achievements of deep learning in other fields, it has also received more and more attention in the fault diagnosis and intelligent identification of modern industrial machinery and equipment. Tran et al. [12] proposed a hybrid deep belief network. The hybrid algorithm was pretrained by the deep belief network, and the training results were classified by the SfAM. The single and combined faults of the suction and exhaust valves were diagnosed by acquiring the three signals of vibration, pressure, and current of the reciprocating compressor. Liu et al. [13] proposed an intelligent fault diagnosis method combining local mean decomposition (LMD) and stack noise reduction automatic encoder (SDAE) to classify faults of reciprocating compressor gas valves. The results show that the classification accuracy of this method can reach 92.72%. Zhao et al. [14] proposed a compound interpolation envelope LMD (CIELMD) algorithm to analyze the vibration signals of the reciprocating compressor in different states. The results show that this method can effectively diagnose the fault of excessive bearing clearance of reciprocating compressors.
A convolutional neural network [15], as one of the classic algorithms in deep learning, is the first learning algorithm that can truly successfully train multilayer networks. Its special convolutional layer and pooling layer can effectively reduce the complexity of the network, greatly reduce the number of free parameters during training, and have strong robustness and fault tolerance. The entire network is easy to train and optimize. These characteristics of convolutional neural networks have attracted a large number of scholars to study their theory and applications. Some experts have also tried to use convolutional neural networks to solve the problem of fault diagnosis. Yang et al. [16] used three sensors to collect vibration signals when the reciprocating compressor valve fails, and directly used as the input of the convolutional neural network to make full use of its characteristics that automatically extract characteristic signals, and then carried out fault diagnosis and obtained higher fault recognition rate. Ince et al. [17] used the onedimensional motor current signal as the input of a convolutional neural network and then diagnosed the fault of the motor. The fault recognition rate is as high as 97.8%. Zhang et al. [18] extracted vibration signals during diesel engine fire faults, and automatically extracted fault features through onedimensional convolutional neural networks, and then classified the faults. The results proved that the accuracy rate under different operating conditions can reach 100%.
At present, vibration signals are mostly used as sample data for fault diagnosis of reciprocating compressor valves. However, according to the thermal performance parameter fault diagnosis method of reciprocating compressors established by Qi [19], it is known that the most direct response of a gas valve when a leak occurs is the pressure ratio imbalance, and the temperature and pressure responses are the most obvious [20–24]. As mentioned earlier, the CNNs’ ability to learn to extract the optimal features and the system can achieve the optimal fault classification and fault detection accuracy through proper training. Therefore, this paper takes the differential pressure and differential temperature of each stage of the multistage reciprocating compressor as sample data. A fault diagnosis model of a reciprocating compressor valve based on onedimensional convolutional neural network is established, the features in the sample data are automatically extracted, and then fault classification is performed by the Softmax function. The advantages of this method are (1) directly using the temperature signal and pressure signal can well achieve online fault diagnosis and improve the quality and efficiency of detection and (2) automatically extracting the fault characteristics of the reciprocating compressor valve through the model, without the need to manually extract the characteristics and selection, which improves the accuracy of fault diagnosis.
2. Structure and Principle of CNN
Convolutional neural networks are a typical type of deep neural networks. They are supervised learning algorithms such as BP neural networks, so they need to be equipped with labeled data as training samples. The essence is to add multiple filters before the traditional fully connected network. The input data are convolved and pooled layer by layer to achieve the effect of extracting input features [25–27]. The biggest improvement of the convolutional neural network lies in the convolutional layer, which greatly reduces the number of parameters, so as to greatly reduce the calculation amount of the network.
A typical CNN topology is shown in Figure 1. The input of the network is a 2D feature map. After that, each convolution layer and pooling layer alternately appear, and the output of this layer is used as the input of the next layer, and so on. At the end of the CNN, several fully connected layers are generally used to link each neuron in the previous layer. This process can expand the feature map into a onedimensional vector for final classification.
CNN convolutional layers are used to extract local features from the input neuron data. Each type of convolution corresponds to an extracted feature. The same type of convolution kernel has the property of weight sharing, which can be learned in parallel and greatly reduces the free parameters. And the complexity of the network is reduced to a certain extent. The feature information extracted by each type of convolution kernel is limited, so multiple convolution kernels are generally used to extract feature information. The input of the convolutional neural network is X, and represents the feature of the lth layer, assuming that is the jth feature map of the convolutional layer of the lth layer, and the process of generation is as follows:where represents the weight matrix connecting the ith feature map of layers and the jth feature map of the layer l; k represents the number of feature maps in layers ; i and j are the indexes of the input and output feature maps; and is the offset corresponding to each feature map of the layer l. Finally, the jth feature map of the layer l is obtained by the excitation function f (·). The excitation function can be the ReLU function, sigmoid function, tanh (x) function, radial basis function, etc.
The CNN pooling layer is also called the downsampling layer, and its role is to scale and map the data graph of the previous layer, so as to reduce the data dimension. Suppose is the jth feature map of the lth pooling layer. The generation process is as shown in equation (2), where down (⋅) represents a pooling function, and each output feature map corresponds to its own multiplicative bias and an additive bias :
Highlevel features are extracted by sequentially stacking convolutional layers and pooling layers, and these features are flattened into onedimensional vectors. For example, the feature extracted from in Figure 1 is flattened into a onedimensional feature vector which is the output of the fully connected layer . Similar to the multilayer neural network thereafter, the output of the layer is as follows:where is an flattened onedimensional vector; k represents the number of feature maps in the layer ; and and represent the tied weights and the tied bias of the layer .
At present, convolutional neural networks are mostly a supervised deep learning model. Network training is similar to traditional artificial neural networks. In the fully connected layer, the back propagation algorithm (BP) is still used to optimize the parameters to reduce the error to the expected value. The parameters are optimized by the gradient descent method.
3. Experimental Study
3.1. Test Rig and Instrumentation
The experimental bench consists of a sixstage watercooled reciprocating compressor and a synchronous motor of the type TK280020/2150, with a maximum speed of 300r/min. The compressor consists of six cylinders, the pressure of the suction gas is 1961.33 Pa (gauge pressure), and the final exhaust pressure after compression of each cylinder is 31.4 MPa (gauge pressure). When the multistage reciprocating compressor gas valve leaks, the hot gas will blow back into the cylinder. This will not only heat the intake gas during the intake process but also continue to heat the gas at the beginning of the compression process, causing the exhaust temperature to rise. The change in temperature will cause the pressure to change, and the change in the previous stage will cause the pressure and heat to increase in the next stage of gas compression. Due to design reasons, the cooling device is fixed at the exhaust port of the cylinder, so the heat generated during the compression of the gas cannot be removed in time, which causes the volume of the gas to expand, the quality of the compressed gas to decrease, and ultimately the efficiency of the compressor to decline. In summary, when the valve leaks, the temperature and pressure of the gas entering and exiting the cylinder are directly affected. Therefore, this experiment uses 12 pressure sensors and 12 temperature sensors to collect the temperature and pressure of the gas at the inlet and outlet of each cylinder. The sensor installation diagram is shown in Figure 2. In the figure, P_{1} and T_{1} represent the sensors for collecting gas pressure and temperature at the inlet of the cylinder I, P_{2} and T_{2} represent the sensors for collecting gas pressure and temperature at the outlet of the cylinder I, and the sensors required for other cylinders are installed similarly.
3.2. Experimental Procedure
The experiment first obtains the data that are collected by all sensors when the valve of each cylinder is normal and then calculates the differential pressure and differential temperature between the inlet and outlet gases of each cylinder. For example, if the pressure and temperature of the inlet gas of the cylinder I are P_{1} and T_{1} and the pressure and temperature of the outlet gas are P_{2} and T_{2}, then the differential pressure and differential temperature between the inlet and outlet gases of the firststage cylinder are P_{d1} and T_{d1}. Because it is a 6stage reciprocating compressor, 12 numbers are obtained after all calculations. These 12 numbers form a set of sample data. Table 1 shows the symbols of the differential pressure and differential temperature between the inlet and outlet gases of each cylinder. Next, the above experiments were repeated under the same operating conditions but simulated six types of real faults, firststage valve leak, secondstage valve leak, thirdstage valve leak, fourthstage valve leak, fifthstage valve leak, and sixthstage valve leak. Table 2 is the fault code table.


4. Established 1DCNN Fault Diagnosis Model
Unlike image data, temperature and pressure data are one dimensional, and it is more reasonable to use the onedimensional convolutional CNN to process these data. In this paper, the feature extraction and fault diagnosis of reciprocating compressor valves are combined, and a onedimensional convolutional neural network (1DCNN) model is proposed. As shown in Figure 3, the model uses the original temperature and pressure signals without the need for manual feature extraction. The classic VGG model replaces the convolution kernel size with 3 3 and the pooling kernel size with 2 2, which can greatly reduce the number of parameters and improve the accuracy of the model, while maintaining the spatial structure of the input of the previous layer. Because the input data of the model established in this paper are 12 1 (the 6stage reciprocating compressor and the differential pressure and differential temperature between the gas inlet and outlet of each stage), the size of the convolution kernel in this model is set to 3 1 and the size of the pooling kernel is set to 2 1. The Keras framework is used to construct a training model in Python. After several verifications and adjustments, the 1DCNN model for reciprocating compressor valve fault diagnosis uses the network structure shown in Table 3, including 2 convolutional layers, 2 maximum pooling layers, 1 flatten layer, and 2 fully connected layers whose length is 128 and 64, respectively. Finally, the Softmax function is used to classify the 7 faults in Table 2. In order to overcome the disappearance of gradients and speed up network training, the activation functions of each layer use the ReLU function.

The 1DCNN model in this paper uses a multiclass crossentropy function to measure the Softmax classification results. represents the target class probability distribution, represents the predicted probability distribution, and and crossentropy function formula is as follows:
In order to minimize the loss function, this article uses the Adam algorithm proposed by Kingma et al. [28] to optimize the gradient descent method. Its advantage lies in the design of independent adaptive learning rates for different parameters when calculating the gradient, and basically only very few parameters need to be adjusted.
1DCNN diagnosis flowchart is shown in Figure 4. The differential pressure and differential temperature data of each stage of gas when a certain level of valve failure is collected are divided into a training set and a test set. The training set data are trained by the 1DCNN to perform forward and back propagation to minimize errors. The 1DCNN training model that meets the requirements uses the test set data for evaluation and verification.
5. Test and Verification
5.1. Description of the Data Set
Through the experimental process described in Section 3.2, a total of 2100 sets of sample data are obtained. Each set of data is 12 numbers representing the differential pressure and differential temperature between the inlet and outlet gases of each cylinder. Table 4 shows some sample data. According to the holdout method [29], the data are decomposed and recombined: 70% are used as training data, including 1470 data samples; 30% are used as test samples, including 630 data samples. The training samples are used for model establishment, and the test samples are used for model evaluation.

5.2. Experimental Results and Discussion
70% of the sample data are randomly selected as the input for network training, and the parameters are adjusted and updated after forward and back propagation. The remaining 30% of the sample data are used to test the fault recognition performance of the trained model, and the accuracy curve and loss curve are shown in Figure 5. The red dashed line in the loss curve represents that the set error is . In order to prevent the random factors from affecting the training results, 10 experiments are repeated. The results are shown in Figure 6. The accuracy of each diagnosis is stable at 100%.
(a)
(b)
The comparison methods include the Levenberg–Marquardt neural network (LMBP), radical basis function neural network (RBF), and conventional BP neural network [30]. Statistics and analysis of the training results are shown in Table 5, and Figure 7 is a comparison of the results of the fault diagnosis test. When using the 1DCNN model for training, the required error is reached when the number of iterations reaches 23, the average test error is 0.0001, and the test accuracy rate reaches the highest level of 100%. LMBP achieves the required error after 87 iterations. The test error and test accuracy are slightly lower than the 1DCNN model. The average test accuracy is 96.78%, and the average test error is 0.001. Although the RBF method and the BP method do not reach the required error, the test accuracy also reaches a high level. The average test accuracy of the RBF method is 92.25%, and the average test error is 0.005. The average test accuracy of the BP method is 89.72%, and the average test error is 0.064. The above results show that the method proposed in this paper can not only have high diagnostic accuracy but also the diagnostic results are very stable.

(a)
(b)
In order to further verify the diagnostic performance of the 1DCNN model, the proportion of training samples and test samples in the data is adjusted, and the proportions of training samples are set to 70%, 50%, and 30%. The best three sets of results are selected from the training results of these three sample sets, and the loss curve and accuracy curve of the three sets of training results are compared, as shown in Figure 8. The results show that, whether it is a loss curve or an accuracy curve, 70% of the training data converge the fastest, followed by 50% of the training data, and 30% of the training data are the slowest.
(a)
(b)
5.3. Visualization of Learning Features
The advantage of 1DCNN is that it can automatically extract the features of the sample data for training. In order to verify the feature extraction capability of this model, the tdistributed stochastic neighbor embedding (tSNE) method is used to reduce the highdimensional data to twodimensional data and visualize. Because the output layer of the 1DCNN model is Softmax classification, the fully connected layer of the previous layer is used as the sample feature extracted by the model, and the 64dimensional feature is reduced to twodimensional feature and visualized by the tSNE method. Figure 9 shows a visualization of the features extracted using the 1DCNN model. The colored numbers in the figure correspond to the fault numbers in Table 2. It is obvious that each fault condition is well aggregated, and the different characteristics are effectively separated. The above results show that the 1DCNN model has excellent feature extraction capabilities and can adaptively extract features required for fault diagnosis.
6. Conclusions
This paper proposes a onedimensional convolutional neural networkbased fault diagnosis method for reciprocating a compressor gas valve. Two types of data such as the differential pressure and differential temperature between each stage of a multistage reciprocating compressor are tested. And compared with the fault classification results of BP neural network, RBF neural network, and LMBP neural network, the experimental results show that the fault recognition rate of the BP neural network reaches 89.72%, and the test error is 0.064; the fault recognition rate of the RBF neural network is 92.25%, and the test error is 0.005; the fault recognition rate of the LMBP neural network is 96.78%, and the test error is stable at 0.001; and the 1DCNN fault recognition rate reaches 100%, and the test error is stable at 0.0001. By comparison, it is concluded that the accuracy of the 1DCNN is higher, which can effectively monitor the operating conditions of the reciprocating compressor valve, thereby reducing or even avoiding accidents. At present, when establishing a model in this paper, it is necessary to continuously adjust the parameters manually. In the future, I hope that a universal model can appear.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was funded by the Tianjin Natural Science Foundation. The project name is “Development of Energy Management and Fault PreDiagnosis System for Industrial System Equipment” (project number 17JCTPJC52200).
References
 M. Elhaj, F. Gu, A. D. Ball, A. Albarbar, M. AlQattan, and A. Naid, “Numerical simulation and experimental study of a twostage reciprocating compressor for condition monitoring,” Mechanical Systems and Signal Processing, vol. 22, no. 2, pp. 374–389, 2008. View at: Publisher Site  Google Scholar
 Y. Tang, Q. Liu, J. Jing, Y. Yang, and Z. Zou, “A framework for identification of maintenance significant items in reliability centered maintenance,” Energy, vol. 118, pp. 1295–1303, 2017. View at: Publisher Site  Google Scholar
 H. Y. Zhao, M. Q. Xu, J. D. Wang, and Y. B. Li, “A Parameters optimization method for planar joint clearance model and its application for dynamics simulation of reciprocating compressor,” Journal of Sound and Vibration, vol. 334, no. 5, pp. 416–433, 2015. View at: Google Scholar
 H. Y. Zhao, M. Q. Xu, and J. D. Wang, “Local mean decomposition based on rational hermite interpolation and its application for fault diagnosis of reciprocating compressor,” Journal of Mechanical Engineering, vol. 51, no. 1, pp. 83–89, 2015. View at: Publisher Site  Google Scholar
 S. M. Ali, K. H. Hui, L. M. Hee et al., “Automated valve fault detection based on acoustic emission parameters and support vector machine,” Alexandria Engineering Journal, vol. 57, no. 1, pp. 491–498, 2017. View at: Google Scholar
 H. Cui, L. Zhang, R. Kang, and X. Lan, “Research on fault diagnosis for reciprocating compressor valve using information entropy and SVM method,” Journal of Loss Prevention in the Process Industries, vol. 22, no. 6, pp. 864–867, 2009. View at: Publisher Site  Google Scholar
 J. Ma, Z. N. Jiang, and J. J. Gao, “Feature extraction method based on chaotic fractal theory and its application in fault diagnosis of gas valves,” Journal of Vibration and Shock, vol. 31, no. 19, pp. 26–30, 2012. View at: Google Scholar
 P. Loukopoulos, G. Zolkiewski, I. Bennett et al., “Abrupt fault remaining useful life estimation using measurements from a reciprocating compressor valve failure,” Mechanical Systems and Signal Processing, vol. 121, pp. 359–372, 2019. View at: Publisher Site  Google Scholar
 H. P. Bloch, “A Practical Guide to Compressor Technology,” Chemical Engineering Progress, vol. 114, no. 1, pp. 89, 2007. View at: Google Scholar
 J. R. Kolodziej and J. N. Trout, “An imagebased pattern recognition approach to condition monitoring of reciprocating compressor valves,” Journal of Vibration and Control, vol. 24, no. 19, pp. 4433–4448, 2018. View at: Publisher Site  Google Scholar
 G. E. Hinton, S. Osindero, and Y.W. Teh, “A fast learning algorithm for deep belief nets,” Neural Computation, vol. 18, no. 7, pp. 1527–1554, 2006. View at: Publisher Site  Google Scholar
 V. T. Tran, F. AlThobiani, T. Tinga, A. Ball, and G. Niu, “Single and combined fault diagnosis of reciprocating compressor valves using a hybrid deep belief network,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 232, no. 20, pp. 3767–3780, 2018. View at: Publisher Site  Google Scholar
 Y. Liu, L. X. Duan, Z. Yuan et al., “An intelligent fault diagnosis method for reciprocating compressors based on LMD and SDAE,” Sensors, vol. 19, no. 5, p. 1041, 2019. View at: Google Scholar
 H. Y. Zhao, J. D. Wang, L. Jay et al., “A compound interpolation envelope local mean decomposition and its application for fault diagnosis of reciprocating compressors,” Mechanical Systems and Signal Processing, vol. 110, pp. 273–295, 2018. View at: Google Scholar
 Y. Lecun, L. Bottou, Y. Bengio, and P. Haffner, “Gradientbased learning applied to document recognition,” Proceedings of the IEEE, vol. 86, no. 11, pp. 2278–2324, 1998. View at: Publisher Site  Google Scholar
 H.b. Yang, J.a. Zhang, L.l. Chen, H.l. Zhang, and S.l. Liu, “Fault diagnosis of reciprocating compressor based on convolutional neural networks with multisource raw vibration signals,” Mathematical Problems in Engineering, vol. 2019, Article ID 6921975, 7 pages, 2019. View at: Publisher Site  Google Scholar
 T. Ince, S. Kiranyaz, L. Eren, M. Askar, and M. Gabbouj, “Realtime motor fault detection by 1D convolutional neural networks,” IEEE Transactions on Industrial Electronics, vol. 63, no. 11, pp. 7067–7075, 2016. View at: Publisher Site  Google Scholar
 K. Zhang, J. F. Tao, C. J. Qin et al., “Fault diagnosis of diesel engine misfire based on random discard and batch standardized deep convolutional neural network,” Journal of Xi'an Jiaotong University, vol. 53, no. 8, pp. 159–166, 2019. View at: Google Scholar
 W. M. Qi, “A method for judging thermal failure of reciprocating compressors,” Mechanical and Electrical Equipment, vol. 5, pp. 32–34, 2002. View at: Google Scholar
 P. Loukopoulos, G. Zolkiewski, I. Bennett et al., “Reciprocating compressor prognostics of an instantaneous failure mode utilizing temperature only measurements,” Applied Acoustics, vol. 147, pp. 77–86, 2019. View at: Google Scholar
 G. D. Dang, Fault Diagnosis and Failure Analysis of Reciprocating Compressor Gas Valve, Xi’an Shiyou University, Xi’an, China, 2018.
 M. FarzanehGord, A. Niazmand, M. DeymiDashtebayaz, and H. R. Rahbari, “Effects of natural gas compositions on CNG (compressed natural gas) reciprocating compressors performance,” Energy, vol. 90, pp. 1152–1162, 2015. View at: Publisher Site  Google Scholar
 M. FarzanehGord and H. Khoshnazar, “Valve fault detection for singlestage reciprocating compressors,” Journal of Natural Gas Science and Engineering, vol. 35, pp. 1239–1248, 2016. View at: Publisher Site  Google Scholar
 J. Townsend, M. A. Badar, and J. Szekerces, “Updating temperature monitoring on reciprocating compressor connecting rods to improve reliability,” Engineering Science and Technology, an International Journal, vol. 19, no. 1, pp. 566–573, 2016. View at: Publisher Site  Google Scholar
 A. Krizhevsky, I. Sutskever, and G. Hinton, “ImageNet classification with deep convolutional neural networks,” Communications of the ACM, vol. 25, pp. 1097–1105, 2012. View at: Google Scholar
 O. Abdeljaber, O. Avci, S. Kiranyaz, M. Gabbouj, and D. J. Inman, “Realtime vibrationbased structural damage detection using onedimensional convolutional neural networks,” Journal of Sound and Vibration, vol. 388, pp. 154–170, 2017. View at: Publisher Site  Google Scholar
 J. An, P. Ai, S. Xu et al., “An intelligent fault diagnosis method for rotating machinery based on onedimensional convolutional neural network,” Journal of Nanjing University (Natural Science), vol. 55, no. 1, pp. 133–142, 2019. View at: Google Scholar
 D. P. Kingma and J. L. Ba, “Adam: A Method for Stochastic Optimization,” Computer Science, vol. 8, pp. 2333–9721, 2015. View at: Google Scholar
 C. Dwork, V. Feldman, M. Hardt, T. Pitassi, O. Reingold, and A. Roth, “The reusable holdout: preserving validity in adaptive data analysis,” Science, vol. 349, no. 6248, pp. 636–638, 2015. View at: Publisher Site  Google Scholar
 H. Zhang, Research on Fault Diagnosis of MultiStage Reciprocating Compressor Air Valve Based on LMBP Neural Network, Donghua University, Shanghai, China, 2015.
Copyright
Copyright © 2020 Fuyan Guo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.