Next Article in Journal
A Self-Powered Hybrid SSHI Circuit with a Wide Operation Range for Piezoelectric Energy Harvesting
Previous Article in Journal
Near-Infrared Spectroscopy and Aquaphotomics for Monitoring Mung Bean (Vigna radiata) Sprout Growth and Validation of Ascorbic Acid Content
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Marine Icing Sensor with Phase Discrimination

Faculty of Engineering and Applied Science, Memorial University of Newfoundland, St. John’s, NL A1B 3X5, Canada
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(2), 612; https://doi.org/10.3390/s21020612
Submission received: 20 November 2020 / Revised: 12 January 2021 / Accepted: 14 January 2021 / Published: 17 January 2021
(This article belongs to the Section Remote Sensors)

Abstract

:
In this paper, a novel approach is presented to the measurement of marine icing phenomena under the presence of a two-phase condition. We have developed a sensor consisting of an electrostatic array and a signal processing based on a decision tree method. A three-element electrostatic array is employed to derive signals having linearly decoupled characteristics from which two key parameters, ice and water accretion layer dimension, can be determined for the purpose of environmental monitoring. The quantified characteristics revealed a correlation with the ice layer thickness in spite of the strong influence from the top water phase layer. The decision tree model established a relationship between the signal characteristics and the two accretion thickness parameters of water and ice layer. Through experimental verification, it has been observed that our sensor array in combination with the decision tree model based signal processing provides a simple practical solution to the challenging field of a two phase composition measurement such as in the marine icing considered in this study.

1. Introduction

Marine Icing is an adverse phenomena affecting offshore vessels and other structures like wind farms [1] or oil and gas platforms [2]. Icing in general, including the atmospheric icing on airplane wings or transmission lines, poses a great deal of difficulties to many operations and can also be hazardous to personnel in the area especially the operators. The knowledge about the rate of ice growth besides the total accumulation would enable managing these hazardous conditions. One way to mitigate the challenging levels of icing is to heat trace the critical infrastructure or shutting down the operations.
Graz Technical University [3] applied a capacitive tomography to atmospheric icing on high voltage power lines. Graz team considered the two phase phenomena, ice and water, however, due to the focus on the atmospheric icing, the water phase was considered finely dispersed and embedded in the ice phase. We focus on marine icing with a significant water layer on top of the ice.
Combitech IceMonitor [4] measures the ice mass on a rotating rod by a load cell. The system requires a stationary installation which prevents this system being applied to ships and other marine structures due to the dynamic forces, vibration, wind gusts or dynamic water splashes.
The Goodrich 0871LH1 ice detectors [5] use an axially vibrating probe to detect the presence of light icing conditions. Goodrich ice detector is designed for thin ice layer applications like avionics and to the date no reports of detecting the ice under the water phase presence has been released.
HoloOptics T42 [6] employs IR signal passed through the medium and an external rain detector to eliminate the sources of false indications by water phase. No testing in marine icing conditions was conducted with the T42.
The Ice Meister Model 9734-SYSTEM industrial ice detector [7] monitors opacity and optical refraction of the ice along the contact with the probe surface. It only recognizes whether air, water or ice is present. This concept is somewhat similar to chilled mirror dew point sensor which often employs the optical reflectivity.
IDS-20 system [8] measures the complex impedance of the icing medium using capacitive plates hermetically sealed. The sensor can distinguish between water, however not in a combined multiphase state.
Jeung Sang Go with Xiang Zhi [9] employed capacitive sensors to measure ice growth in real time. The developed system has also been patented [10], however the authors conclude the water layer formation has to be prevented in order to maintain the accuracy. Charles Ezeoru [11] conducted a similar research using the same capacitive technique with interdigitated comb-style electrodes. Both approaches experienced the transformation from liquid state to solid state which has been reflected in a ramp capacitance profile in time. Unfortunately, no one has tried to quantify the transitional multiphase period.
We have developed a sensor array that allows a flexibility to be applied on curved surfaces such as wind turbine wings. Our earlier work described in [12] utilized a planar array of electrostatic sensors of different spacing that have been instrumented by a custom circuitry based on LC oscillator. Here we report on using a constant-spacing electrostatic sensor array interfaced and conditioned by a commercial four-channel capacitive pickup board by Texas Instruments.
Previously we processed the array signals by a multidimensional least squares method [12] and in another approach by artificial neural networks [13]. In this work, we applied a machine learning strategy centered around a decision tree method.
Machine learning algorithms are being used in many applications ranging from flood prediction [14], solar radiation [15], to wind generator blade monitoring [16,17]. Machine learning based on a decision tree method has been used in internet security systems [18] as well as in detecting stability of a power system voltage [19] or in classification tasks [20,21].
The ice can take different formations described in [22], however, the dielectric properties remain nearly constant which is being exploited in this research work. The problem of marine icing detection is however complicated through the presence of a water phase on top of the ice layer since water’s dielectric constant, the relative permittivity, is significantly larger than that of ice. The decision tree method has been found to meet our needs in recognizing and classifying the multiphase situations and provide a more accurate ice accretion measurement. This method translates the acquired signal data to a regression and classification model to determine the thickness of both layers of ice and water simultaneously.

2. Methodology

Our research into sensing and detecting the marine icing as a two phase phenomena has been initially founded on the principles of linearly decoupled array of electrostatic sensors of variable electrode gap spacing. The experimental verification confirmed our hypothesis of dissimilar gap spacing being capable to encode the stray electrostatic field above the sensor plane in a unique and linearly independent way. This phenomena was utilized in discriminating each individual phase, ice and water. Our team developed a multi-dimensional least-squares-method [12] to map the measured data to the measurands. The least-squares method was later substituted by a more robust and more accurate technique based on artificial neural networks method [13].
In this paper, we modified our earlier methods in four ways, (1) at the transducer array level (constant gap spacing, variable insulator height), (2) at the signal pick up circuitry (commercial 4ch converter), (3) at the signal processing algorithm level (decision tree method) and (4) at the spacial adaptation to curved surfaces (structural flexibility). Our earlier setup resulted in a large variation of capacitance range which required a custom build circuitry for capacitance to frequency conversion. In this work, we adapted an off-the-shelf signal detection board FDC2214EVM from Texas Instruments [23] which features four independent channels simultaneously scanned and a simple interfacing (Figure 1). FDC2214EVM module features 4-ch 28 bit capacitance-to-digital signal converters (FDC). The manufacturer also provides a software enabling the data streaming into a PC. In order to keep all array elements within a consistent accuracy range, we have built a new three element array of uniform elements (uniform spacing) that was attached to a curved surface of 80 mm radius (Figure 2).
In order to linearly decouple the measured capacitance across the array, we introduced three different dielectric layer heights. The capacitance of two elements was analyzed by a finite element method (FEM) using Ansys Maxwell software (Figure 3). The solid line characteristics represent the capacitance for an element having 0.25 mm coating of PET, plotted as constant capacitance contours across a range of ice and water accretion levels. The dotted line corresponds to 0.35 mm dielectric coating layer. The graph demonstrates a linear independence in a similar way to the gap variation analysis in [12] though not as profoundly. The FEM parameters are listed in Table 1.
Any traditional insulation material that is not permeable to water will provide the function of eliminating the conductance (real part) from the impedance measurement. PET sheets were readily available to us and their dielectric constant is similar to the epoxy resin used in the past. The critical part in obtaining the linearly independent characteristics across the array is in a different number of stacked layers of PET sheets above each array element.

3. Experimental Validation

We have conducted a series of ice layer deposition experiments in a similar way described in [12] using a water saturated roller. The system was placed in a deep freezer at −20 ° C and was initially tempered for 1 h. Then, over a period of two hours, we have periodically rolled a water layer onto the sensor surface which is illustrated in Figure 4. Each peak corresponds to a new wetting cycle during which the capacitance rises sharply. If we were to consider only the ice layer as a number of previously developed sensors were based upon, the reported ice layer accretion would follow the same sharp increase causing a large discrepancy from the real situation. A low pass filter could be applied to remedy these time limited deviations however the water phase could also be present at all times in situations of icy waters constantly battering the place of interest.
Figure 4 shows six transitional cycles where each cycle starts with a fresh water layer that gradually solidifies into a new ice layer. Each cycle took approximately fifteen minutes to settle. Our experimental approach is based on an assumption that each deposited layer is consistent with all the other layers without the need for knowing its exact dimensions which constitutes a parametric approach.

4. Signal Processing

The Decision tree model is constructed from the main root branches consisting of the interior nodes and the final nodes also called the decision nodes [24]. The tree represents a mapping characteristic between the signal data and the target parameter, the ice accretion besides the water phase information as a byproduct.
We have split the acquired signal data and the corresponding ice and water accretion estimate into the training data set and the test data set. The training data is used to derive the decision tree parameters during the training phase. Figure 5 shows both the training and the testing processes. The training data set consists of 50% of data points from all the acquired data (odd rows in a column vector). Once the relationship between input data and the measurand was established, the testing phase has been conducted on the test data (even rows). Figure 6 shows the ’pairplot’ data set illustrating the relationship and its weight between the sensor data and the output classes.
Our decision tree implementation is based on the Entropy rule described in [25]. In comparison to other models like the Gini rule, the entropy model performed better in our specific application. Both rules are characterized in Equation (1) with p i being the entropy measure.
Gini = 1 i = 1 c ( p i ) 2 Entropy = i = 1 c p i l o g 2 ( p i )
Figure 7 shows our assumption in the time profile of our anticipated data of ice accretion (dotted line) and the water level height (solid line) [12]. At the end of each cycle, the water layer is completely converted into the ice.

5. Results

We implemented a twelve level deep decision tree processes the measured data based on entropy value where a low entropy value leads to the best leaf node in terms of purity. Applying the trained decision tree model to the test data results in the ice accretion estimate depicted in Figure 8 and the water layer estimate depicted in Figure 9. Both profiles show a significant improvement over the earlier results reported in [12] and a corresponding results to the ANN approach in terms of the ice accretion [13].
The Precision, Recall, F-Measure and Accuracy are commonly used for evaluating classification methods in Machine learning [26]. Precision is the ratio between the correct positive predicted classes for a layer and the total number that includes the correct and incorrect positive classes predicted for the same layer. Recall is the ratio between the correct positive predicted classes for a layer and the total number of the correct and incorrect negative classes predicted for the same layer. F-Measure makes a harmonic relation between Precision and Recall. These parameters were calculated as follows:
Precision = T P T P + F P Recall = T P T P + F N F - Measure = 2 × Precision × Recall Precision + Recall
True Positive (TP) is a classification method outcome where the decision tree method generates a correct positive class layer. Likewise, True Negative (TN) is a classification method outcome where the decision tree method generates a correct negative class layer. False Positive (FP) and False Negative (FP) are a classification method outcome where the decision tree method generates incorrect positive and negative class, respectively.
The Confusion-matrix yields the most ideal suite of metrics for evaluating the performance of a classification algorithm. Gupta [27] provides a detailed description on the confusion matrix as a measure of determining the accuracy of the decision tree classification model. The confusion matrix for ice and water classifier is a 6 × 6 matrix illustrated in Figure 10 with the elements along the diagonal representing ’True Predicted’ data. The overall accuracy for the prediction model of Ice and Water Classifier recorded accuracy of 92.8%.
The conventional method of using the least squares method to our experimental data as reported in [12] is presented in Figure 11 and Figure 12. Comparing these results with the decision tree results clearly indicate the level of improvement by using the new signal processing method. A comparison with the ANN method applied to the same data (Figure 13 and Figure 14) reveals a similar error characteristic in the ice accretion. The comparison analysis is quantitatively summarized by a root mean squares error analysis presented in Table 2. Both the decision tree and the ANN methods are equivalent in terms of the ice accretion measurement estimation, however, the current tree approach offers more lightweight implementation suitable for a microcontroller based system.

6. Conclusions

A novel sensor array has been developed to measure the marine icing accretion and to allow flexibility in deployment over curved surfaces. A new signal processing method based on a decision tree method has been applied to the array data. A theoretical analysis confirmed linearly decoupled signal characteristics that are critical in discriminating the water phase influence on the ice accretion measurement. An experimental analysis validated the high relevance of employing the decision tree method for signal processing. The ice accretion estimate obtained through the decision tree method demonstrates a significant improvement over the conventional least square method and a similar characteristic with the Neural Network method, however of a lesser computational overhead and a smaller footprint. Our three element capacitive array is easy to fabricate with off-the-shelf components and the 4 ch 28 bit signal processing interface circuit from TI is also aligned with a very low cost characteristics. The nature of our non-contact based design without any moving parts and the low cost characteristics promises a robust solution towards the marine icing measurement under the harsh environmental conditions.

Author Contributions

A.E. and V.M. conceptualized the initial hypothesis and designed the experiments. The formal analysis and theoretical concept validation was performed by A.E. under V.M. and S.B. supervision. A.E. constructed the experimental setup and conducted the experiments. A.E. proposed the data processing algorithm based on the decision tree method and implemented it in MATLAB and Python. The original draft was prepared by A.E. while V.M. and S.B. contributed to the manuscript revision. All authors read and agreed to the published version of the manuscript.

Funding

This research received funding from Mitacs “Accelerate Lab2Market” grant IT20713 (Sponsor: Springboard Atlantic) held by V.M., and from the research grant “ICE WEAR AND ABRASION EFFECTS ON CONCRETE SURFACES” held by S.B (Sponsors: Kvaerner Canada Ltd., NL Innovation Council, NSERC and Memorial University).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The MATLAB and Python code introduced in the paper can be found at https://sciprofiles.com/profile/1438257.

Acknowledgments

The authors would like to express a sincere gratitude to the Libyan Ministry of Education for providing a postgraduate fellowship to the first author. Equal gratitude goes to the Memorial University for supporting and enabling this research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gregow, E.; Bernstein, B.; Wittmeyer, I.; Hirvonen, J. LAPS–LOWICE: A Real-Time System for the Assessment of Low-Level Icing Conditions and Their Effect on Wind Power. J. Atmos. Ocean. Technol. 2015, 32, 1447–1463. [Google Scholar] [CrossRef]
  2. Dehghani-Sanij, A.R.; Dehghani, S.R.; Naterer, G.F.; Muzychka, Y.S. Marine icing phenomena on vessels and offshore structures: Prediction and analysis. Ocean. Eng. 2017, 143, 1–23. [Google Scholar] [CrossRef]
  3. Neumayer, M.; Bretterklieber, T.; Flatscher, M. Signal Processing for Capacitive Ice Sensing: Electrode Topology and Algorithm Design. IEEE Trans Instrum. Meas. 2019, 68, 1458–1466. [Google Scholar] [CrossRef]
  4. Combitech, The Ice Load Surveillance Sensor IceMonitor, Product Sheet, 2013. Available online: http://www.rwis.net/res/pdffiles/IceMonitor_Product_Sheet.pdf (accessed on 23 December 2020).
  5. Goodrich Ice Detector Models 0871LH1; UTC Aerospace Systems: Burnsville, MN, USA, 2013.
  6. T40 Series of Icing Rate Sensors. HoloOptics, User Guide. Katarinavägen 22, S116 45 Stockholm Sweden/Suede, 2010. Available online: http://holooptics.utrymmet.com/Products_T20series.htm (accessed on 23 December 2020).
  7. Ice Meister Model 9734—Industrial Ice Detecting Sensor System. New Avionics Corporation, Technical Data Sheet. 2501 East Commercial Blvd, Fort Lauderdale FL 33308 USA, 2016. Available online: http://newavionics.com/9734-open-sky/ (accessed on 23 December 2020).
  8. IDS-20 Product Information. SOMMER Messtechnik. Sommer GmbH, Strassenhaeuser 27, A-6842 Koblach, V-11/2016. Available online: https://www.elscolab.com/en/products/non-contact-ice-detection-sensor-ids-20 (accessed on 23 December 2020).
  9. Zhi, X.; Cho, H.C.; Wang, B.; Ahn, C.H.; Moon, H.S.; Go, J.S. Development of a Capacitive Ice Sensor to Measure Ice Growth in Real Time. Sensors 2015, 15, 6688–6698. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Go, J.S. Ice Thickness Measurement Sensor. U.S. Patent 10,066,923 B2, 4 September 2018. Available online: http://patentimages.storage.googleapis.com/bd/75/23/9016fcbfb5cbc2/US10066923.pdf (accessed on 23 December 2020).
  11. Ezeoru, C. Marine Icing Sensor Design Using Capacitive Techniques. Master’s Thesis, Memorial University of Newfoundland, St. John’s, NL, Canada, August 2016. [Google Scholar]
  12. Elzaidi, A.; Masek, V.; Muzychka, Y. Phase Discrimination in Marine Icing Using a Coplanar Capacitive Array. IEEE Sens. J. 2019, 19, 11301–11308. [Google Scholar] [CrossRef]
  13. Elzaidi, A.; Masek, V.; Bruneau, S. Water and Ice Detection in Marine Icing by Capacitive Sensor Array and the Artificial Neural Network Model. In Proceedings of the 2019 IEEE 10th Annual Information Technology, Electronics and Mobile Communication Conference (IEMCON), Vancouver, BC, Canada, 17–19 October 2019. [Google Scholar]
  14. Mosavi, A.; Ozturk, P.; Chau, K.-w. Flood prediction using machine learning models: Literature review. Water 2018, 10, 1536. [Google Scholar] [CrossRef] [Green Version]
  15. Torabi, M.; Mosavi, A.; Ozturk, P.; Varkonyi-Koczy, A.; Istvan, V. A hybrid machine learning approach for daily prediction of solar radiation. In Recent Advances in Technology Research and Education; Springer: Cham, Switzerland, 2018; pp. 266–274. [Google Scholar]
  16. Jiménez, A.A.; Márquez, F.P.G.; Moraleda, V.B.; Muñoz, C.Q.G. Linear and nonlinear features and machine learning for wind turbine blade ice detection and diagnosis. Renew. Energy 2019, 132, 1034–1048. [Google Scholar] [CrossRef] [Green Version]
  17. Mahmoud, T.; Dong, Z.Y.; Ma, J. An advanced approach for optimal wind power generation prediction intervals by using self-adaptive evolutionary extreme learning machine. Renew. Energy 2018, 126, 254–269. [Google Scholar] [CrossRef]
  18. Alharbi, S.; Rodriguez, P.; Maharaja, R.; Iyer, P.; Subaschandrabose, N.; Ye, Z. Secure the internet of things with challenge response authentication in fog computing. In Proceedings of the 2017 IEEE 36th International Performance Computing and Communications Conference (IPCCC), San Diego, CA, USA, 10–12 December 2017. [Google Scholar]
  19. Vanfretti, L.; Narasimham Arava, V.S. Decision tree-based classification of multiple operating conditions for power system voltage stability assessment. Int. J. Electr. Power Energy Syst. 2020, 123, 106251. [Google Scholar] [CrossRef]
  20. Phiri, D.; Simwanda, M.; Nyirenda, V.; Murayama, Y.; Ranagalage, M. Decision Tree Algorithms for Developing Rulesets for Object-Based Land Cover Classification. ISPRS Int. J. Geo-Inf. 2020, 9, 329. [Google Scholar] [CrossRef]
  21. Polat, K.; Güneş, S. Classification of epileptiform EEG using a hybrid system based on decision tree classifier and fast Fourier transform. Appl. Math. Comput. 2007, 187, 1017–1026. [Google Scholar] [CrossRef]
  22. Knight, C.A. Curved growth of ice on surfaces. J. Appl. Phys. 1962, 33, 1808–1815. [Google Scholar] [CrossRef]
  23. Texas Instruments, FDC2214 with Two Capacitive Sensors Evaluation Module. Available online: http://www.ti.com/tool/FDC2214EVM (accessed on 23 December 2020).
  24. Song, Y.-Y.; Ying, L.U. Decision tree methods: Applications for classification and prediction. Shanghai Arch. Psychiatry 2015, 27, 130. [Google Scholar] [PubMed]
  25. Tan, P.-N.; Steinbach, M.; Kumar, V. Introduction to Data Mining; Pearson New International Edition PDF eBook; Pearson Schweiz Ag: London, UK, 2013; Volume 4, pp. 145–206. [Google Scholar]
  26. Han, J.; Pei, J.; Kamber, M. Data Mining: Concepts and Techniques; Elsevier: Amsterdam, The Netherlands, 2011. [Google Scholar]
  27. Gupta, D.L.; Malviya, A.K.; Singh, S. Performance analysis of classification tree learning algorithms. Int. J. Comput. Appl. 2012, 55, 39–44. [Google Scholar]
Figure 1. Interfacing Curved Capacitive Sensors Array with Computer by using FDC2214EVM.
Figure 1. Interfacing Curved Capacitive Sensors Array with Computer by using FDC2214EVM.
Sensors 21 00612 g001
Figure 2. Spiral Capacitive Sensor Element.
Figure 2. Spiral Capacitive Sensor Element.
Sensors 21 00612 g002
Figure 3. Contours of Constant Capacitance for Sensor Geometries 0.25 mm and 0.35 mm.
Figure 3. Contours of Constant Capacitance for Sensor Geometries 0.25 mm and 0.35 mm.
Sensors 21 00612 g003
Figure 4. Experimental Data from Curved Capacitive Sensors Array.
Figure 4. Experimental Data from Curved Capacitive Sensors Array.
Sensors 21 00612 g004
Figure 5. The Decision Tree Processes.
Figure 5. The Decision Tree Processes.
Sensors 21 00612 g005
Figure 6. Pairplot data set for curved capacitive sensors array.
Figure 6. Pairplot data set for curved capacitive sensors array.
Sensors 21 00612 g006
Figure 7. Ice and Water accretion in time during the experiment.
Figure 7. Ice and Water accretion in time during the experiment.
Sensors 21 00612 g007
Figure 8. The ice accretion estimate in comparison to the assumed data characteristic.
Figure 8. The ice accretion estimate in comparison to the assumed data characteristic.
Sensors 21 00612 g008
Figure 9. The water layer estimate in comparison to the assumed data characteristic.
Figure 9. The water layer estimate in comparison to the assumed data characteristic.
Sensors 21 00612 g009
Figure 10. Confusion Matrix for Decision Tree Ice Classifier.
Figure 10. Confusion Matrix for Decision Tree Ice Classifier.
Sensors 21 00612 g010
Figure 11. Detected Ice vs. Predictions using Least Squares.
Figure 11. Detected Ice vs. Predictions using Least Squares.
Sensors 21 00612 g011
Figure 12. Detected Water vs. Predictions using Least Squares.
Figure 12. Detected Water vs. Predictions using Least Squares.
Sensors 21 00612 g012
Figure 13. Detected Ice vs. Predictions using Neural Network.
Figure 13. Detected Ice vs. Predictions using Neural Network.
Sensors 21 00612 g013
Figure 14. Detected Water vs. Predictions using Neural Network.
Figure 14. Detected Water vs. Predictions using Neural Network.
Sensors 21 00612 g014
Table 1. FEM Model Parameters.
Table 1. FEM Model Parameters.
Material ε r
air1.0006
water81
ice4.2
PET3.6
Table 2. Error Analysis.
Table 2. Error Analysis.
MethodMediumRoot Mean Square ErrorMean Absolute Error
Decision Treewater0.02310.0082
ice0.02910.01704
Least Squares [12]water0.04530.0336
ice0.05560.0424
Neural Network [13]water0.01440.0072
ice0.02140.01702
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Elzaidi, A.; Masek, V.; Bruneau, S. Marine Icing Sensor with Phase Discrimination. Sensors 2021, 21, 612. https://doi.org/10.3390/s21020612

AMA Style

Elzaidi A, Masek V, Bruneau S. Marine Icing Sensor with Phase Discrimination. Sensors. 2021; 21(2):612. https://doi.org/10.3390/s21020612

Chicago/Turabian Style

Elzaidi, Abdulrazak, Vlastimil Masek, and Stephen Bruneau. 2021. "Marine Icing Sensor with Phase Discrimination" Sensors 21, no. 2: 612. https://doi.org/10.3390/s21020612

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