Next Article in Journal
RM-ADR: Resource Management Adaptive Data Rate for Mobile Application in LoRaWAN
Next Article in Special Issue
Design of a New Seismoelectric Logging Instrument
Previous Article in Journal
Orthogonality-Constrained CNMF-Based Noise Reduction with Reduced Degradation of Biological Sound
Previous Article in Special Issue
The Minimum AC Signal Model of Bipolar Transistor in Amplification Region for Weak Signal Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Ultrasonic-Capacitive System for Online Characterization of Fuel Oils in Thermal Power Plants

by
Mateus Mendes Campos
1,2,
Luiz Eduardo Borges-da-Silva
2,
Daniel de Almeida Arantes
2,
Carlos Eduardo Teixeira
1,
Erik Leandro Bonaldi
1,
Germano Lambert-Torres
1,
Ronny Francis Ribeiro Junior
1,
Gabriel Pedro Krupa
1,
Wilson Cesar Sant’Ana
1,*,
Levy Ely Lacerda Oliveira
1 and
Renato Guth de Paiva
3
1
Gnarus Institute, R&D Department, Itajuba 37500-052, MG, Brazil
2
Institute of Engineering Systems and Information Technology, Itajuba Federal University, Pro-Reitoria de Pesquisa e Pos-Graduacao (PRPPG), Itajuba 37500-903, MG, Brazil
3
Norte Energia S.A., Vitoria do Xingu 68383-000, PA, Brazil
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(23), 7979; https://doi.org/10.3390/s21237979
Submission received: 4 October 2021 / Revised: 1 November 2021 / Accepted: 24 November 2021 / Published: 29 November 2021
(This article belongs to the Special Issue Sensors in Electronic Measurement Systems)

Abstract

:
This paper presents a ultrasonic-capacitive system for online analysis of the quality of fuel oils (FO), which are widely used to produce electric energy in Thermal Power Plants (TPP) due to their elevated heating value. The heating value, in turn, is linked to the quality of the fuel (i.e., the density and the amount of contaminants, such as water). Therefore, the analysis of the quality is of great importance for TPPs, either in order to avoid a decrease in generated power or in order to avoid damage to the TPP equipment. The proposed system is composed of two main strategies: a capacitive system (in order to estimate the water content in the fuel) and an ultrasonic system (in order to estimate the density). The conjunction of the two strategies is used in order to estimate the heating value of the fuel, online, as it passes through the pipeline and is an important tool for the TPP in order to detect counterfeit fuel. In addition, the ultrasonic system allows the estimation of the flow rate through the pipeline, hence estimating the amount of oil transferred and obtaining the total mass transferred as a feature of the system. Experimental results are provided for both sensors installed in a TPP in Brazil.

1. Introduction

The Brazilian electrical energy mix, since the 1990s, has shown a strong growth in Thermal Power Plants (TPP) based on fossil fuels [1]. According with the Brazilian regulatory agency for energy (ANEEL) [2], the granted electrical power to fossil fuels is approximately 35 GW (representing about 16% of the Brazilian energy mix). Among the fossil fuels, fuel oils (FOs, which include heavy and light fuel oils) represent a share of 23% (around 8 GW).
The efficiency of power generation in a TPP is linked with the quality of the FOs, besides the fact that oils with contaminants may induce to long-term damage to the equipment in the TPP. The most common contaminant is water. The process of contamination might be unintentional (due the processes involved in extraction and storage)—however, it might also be due to counterfeit. The counterfeit/adulteration of a high-quality fuel with some other substance of less quality and cheaper economical value is a common (and deceptive) practice around the world, as presented in scientific reports from Bangladesh [3], Ghana [4], India [5], and Tanzania [6]. In these aforementioned references, gasoline or diesel is mixed with subsidized kerosene or lubricants and affect automotive engines. In the case of the Brazilian TPPs, the fuel oils are usually loaded from tank trucks. These trucks are filled with an agreed amount of oil at a distribution center or refinery. During the transportation, it is reported [7] that some truck drivers remove a small quantity of oil (in order to sell it in the black market) and replace this amount with water.
The literature presents some methods to detect damage on the internal combustion engines that convert the chemical energy of the fuels into rotational mechanical energy at their crankshafts (which are tied to synchronous generators in order to produce electrical energy). Such methods may rely on the direct measurement of internal parameters of the engine (such as pressure and temperature at different parts of the engine) [8,9] or even the use of voltage and current signals of the synchronous generator [10] in order to indirectly detect mechanical issues within the engine. However, these methods cannot test the fuel quality, which may cause a problem that will only be manifested in the long term.
In order to check the quality of their fuels, the Brazilian TPPs perform laboratory analyses (which follow the standards ASTM D95 [11] for water content, ASTM D1298 [12] for density, and ASTM D4868 [13] for heat value) on samples. However, the samples taken are only a small fraction of the total amount of fuel that is loaded to the TPP pipelines. Moreover, in the case of TPPs located far from the large centers (which is a common scenario in the North and Northeast regions of Brazil), the laboratory tests are not even performed (due to lack of certified laboratories nearby). Hence, the use of online methods for the estimation of fuel quality (especially heavy fuel oils—HFO) is a general demand from the TPP shareholders.

1.1. Literature Review

1.1.1. Online Measurement of Water Content in Fuel Oils

The work of Han et al. [14] presents a method based on differential pressure. Its advantage is that it this is a simple (and low cost) method, requiring only a differential pressure transducer. However, the accuracy of this method for estimating water content is strongly dependent on the difference in the specific mass of fuel oil and water.
The literature also presents methods based on absorption of electromagnetic waves, such as infrared (as presented by Al-Saiyed et al. [15]) and gamma ray (as presented by Scheers [16] and Johansen and Tjugum [17]). The main advantage of methods based on electromagnetic waves is the frequency range that can be used, whose wavelengths reach molecular levels, detecting the fraction of water in molecules as well as other contaminants, sulfur, vanadium, etc. However, these techniques are more common for laboratory applications that use small samples of fuel oils to be analyzed and obtain their characteristics. Considering online models, these techniques are limited to small ducts (from mm to few cm) [18], which differs from the applications proposed in the work (where the ducts are larger than 3 ). Another drawback of these techniques is the cost of processing.
The works of Gregory and Clarke [19] and Makeev et al. [20] present a method based on microwaves. The main advantage is the possibility to measure the whole range of water fractions (from 0% to 100%, with an extra advantage at the precision at lower fractions). However, this technology has a high initial cost in relation to the others, and it is also sensitive to variations in salinity in the range of higher fractions of water.
The work of Teixeira et al. [21] presents a method based on ultrasound. The advantage of this method is that it allows the gathering of information online, adaptable to conduits with larger diameters. However, the disadvantages are a relatively high cost and also the need for a cumbersome calibration process in order to classify water percentages.
Another technology for the estimation of the water content is the measurement of the capacitance of the emulsion. The fluid being measured acts as the dielectric of the sensor, making it possible to measure the electrical permittivity of the fluid and, consequently, the capacitance of the sensor–fluid system. The percentage of water in oil can be calculated based on a predictable relationship in the electrical properties of the emulsions. Hammer et al. [22] proposed a non-intrusive capacitive sensor, capable of measuring water contents in the range between 0% and 80%. Shuller et al. [23] proposed a measurement principle of the dielectric properties of the emulsion based on a single electrode excited by an oscillator (in the range of Mega Hertz). Libert et al. [24,25] presented a system based on capacitive sensors installed at the internal wall of the pipes.
Systems such as EIT (Electrical Impedance Tomography) or ECT (Electrical Capacitance Tomography) are techniques widely researched and disseminated in the literature [26,27,28] with great advances in the field of flow measurement and physical– characterization of multiphase fluids. The advantages found in these types of technique are mainly focused on the potential to be applied in multiphase fluids, to present the flow image in the penstock cross section and the possibility to be used on flow conditions unfavorable for other methods [26,27,28]. The disadvantages, however, are linked to the development of more consolidated techniques for commercial application, in addition to the complexity of electronics [26].
The capacitive sensor used in this current work was proposed by the authors of [29] and is based on paralleled disks, which are immersed in the emulsion inside the pipeline. The advantage of the system proposed in [29] is the electronic topology based on a Howland current source, capable of eliminating parasitic effects caused by influences from cables, connectors, and constructive characteristics of the capacitive sensor that affect the capacitance measurement. In addition, as will be seen in Section 2.2, the sensor is based on several disks in order to perform a paralleled association of the capacitances, which increases the total capacitance of the sensor. The fluid to be inspected fills the spaces between each pair of disks, forming small associated capacitors. Higher capacitance values allow it to have greater precision in measurements, as their proportion to the parasitic elements (such as those from cables, traces of plates and electronic parts) is greater, thus reducing the magnitude of errors caused by them.

1.1.2. Online Measurement of Density in Fuel Oils

The literature presents several methods for measuring the density: many methods are based on direct physical relationships and others are based on indirect parameter determinations [30]. Most of the established methods have limitations inherent to the system, which often result in application restrictions in the sensing implementation (physical process limitations such as tube diameter, limitations in process deviations, limitation in the flow velocity of the fluid to be monitored, and others) [31]. In other cases, the presence of bubbles, particles, or incrustations may cause problems in specific mass measurements. Methods involving ultrasonic waves appear as alternatives to replace the standard methods, where these cannot be applied [32]. The most widespread ultrasonic techniques for measurement of density are based on three different principles: speed of sound in the propagation medium, acoustic impedance, and propagation on waveguides.
Waveguide approaches generally use time variations in the propagation of ultrasonic waves on a transmission line immersed in the measurement medium [33].
The acoustic impedance method approach is based on the determination of the reflection and transmission coefficient and the transmission coefficients of flat ultrasonic waves that propagate through different media [32].
Considering the determination of the density of any medium just by measuring the speed of sound, Swoboda et al. [34] used measurements of speed of sound and temperature in a sulfuric acid solution and found that the assessment of density was possible over a relative range greater than 1.10–1.30 and over a temperature range greater than 10 C to 50 C. The density and speed of sound are quantities that vary with temperature and pressure [35]. For determination of density in controlled systems, where there is little variation in temperature, pressure, and knowledge of the compressibility of the medium to be measured, the application of this concept is quite interesting, as it reduces the complexity of electronics in measuring the speed of sound and, by consequence, in the calculation of the density.

1.2. Aim and Objectives

The proposition of this work is to present a methodology for measurement of the quality of fuel oils based on two different strategies (capacitive and ultrasonic), in order to estimate the water content and heating values. The methodology is applied in an ultrasonic-capacitive hybrid on-line measurement system, which analyzes the fuel oil constantly throughout the entire fueling process. The advantage of this system is that it estimates the fuel oil quality parameters in real time, unlike the methods currently applied in Brazilian TPPs, which are based on laboratory analysis of small samples of fuel oil. However, laboratory results are not instantaneous, and some TPPs in Brazil do not even have certified laboratories nearby. Hence, the proposition of this paper allows for an immediate estimation of the fuel quality, while the fuel oil is flowing through the pipelines.
The originality of this paper lies precisely in the application of the proposed system at the loading bays of thermoelectric power plants for online monitoring of fuel oil used for energy generation.

1.3. Outline of the Paper

The paper is organized as follows: Section 2 presents the fundamental concepts about the two main strategies (capacitive and ultrasonic) and its correlation in order to infer about the fuel quality. Section 3 presents the experimental results obtained with the prototype installed in a Brazilian power plant. Finally, Section 4 presents the main conclusions of the work and some opportunities for future research.

2. Materials and Methods

Figure 1 presents the conceptual drawing of the proposed system. This system must be installed in a section of the pipeline between the TPP fuel storage and the point where the tank truck discharges the oil fuel. Section 2.1 briefly discusses the main characteristics of an emulsion of water-oil. As the emulsion passes through the pipe, the proposed capacitive sensor (described in detail in Section 2.2) estimates the water content of the emulsion, and the proposed ultrasonic sensors (described in detail in Section 2.3) estimate the density and the flow rate. The box indicated as “control boards” contains an FPGA board that controls and receives data sent from the capacitive system control board, ultrasonic system control board, and, also, from the temperature sensor. Through ethernet, the FPGA makes the data available to a remote analysis software, where the information of the aforementioned sensors is used to estimate the water content, the density and the heating value (hence, the quality) of the oil fuel. Section 2.4 presents the proposed procedure to estimate the heating value.

2.1. Characteristics of Fuel Oils

Fuel oils (FO) can be categorized in relation to viscosity as Light Fuel Oils (LFO) and Heavy Fuel Oils (HFO). FOs also have several other parameters that are used in their characterization—three of them, particularly in this work, are of special interest:
  • Water content is probably the main contaminant of fuel oils, since contamination can occur in all parts of the production process, from the moment of extraction of crude oil to the moment of sending it to the TPPs [21]. The lower the water-in-oil content, the better the oil quality for power generation. The Resolution number 3 of the Brazilian National Agency of Petroleum, Natural Gas and Biofuels (ANP) stipulates a maximum of 2% water-in-oil content [36].
  • Density is the ratio of the mass of oil to its volume. The density of the FO is an important property, which directly influences its combustion power and in the energy production, as well as in the removal of contaminants in centrifuges [37]. Density can be used to define other intrinsic characteristics of FO, such as the heat of combustion.
  • Heating Value shows the amount of energy that could be obtained from fuel combustion. It is, usually defined in terms of a Higher Heating Value (HHV, or Gross Heat of Combustion (GHC)) and a Lower Heating Value (LHV, or Net Heat of Combustion (NHC)). The difference between the two is that HHV includes the latent heat of evaporation from the water vapor formed during combustion, while the LHV considers the loss of the latent heat needed to evaporate the water formed during combustion [38,39]. The LHV is the indicator used by Brazilian TPPs in order to estimate the generation capacity [36] and is the one that is going to be estimated in Section 2.4.

2.2. Estimation of Water Content through Capacitance

The estimation of water content in fuel oils is a welcomed measure in the power generation industry, as the presence of this contaminant reduces the efficiency of energy production. Moreover, the presence of water in fuel oil tends to be corrosive to system equipment, also causing environmental problems. Still, from an economic point of view, water has no economic value for energy generation.
The authors proposed in [29] a capacitive sensor based on paralleled disk, as presented in Figure 2. As the sensor is immersed into the emulsion, the fluid act as the dielectric of an array of capacitors. The importance of having an array of paralleled capacitors is that higher values of capacitance increase the precision of the measurements, as the parasitic elements (from the cables, circuit boards, etc.) have less influence on the process. In addition, the proposed circuit uses a Howland current source topology, with enhancements that allow its continuous oscillating operation and easy adjustment of output current and frequency.
The sensor must be installed in a transversal section of the pipeline, as presented in Figure 3. Its sensitive region is limited to the metallic disks that will be filled with the fluid. The length of this sensitive region is 8 cm.
The permittivity of water is approximately 40 times greater than the permittivity of the fuel oil [40,41]—hence, the greater the water content in the emulsion, the greater its measured capacitance will be.
In order to estimate the water content from the measured capacitance, a calibration procedure is performed, based on the capacitances of four emulsion samples with known water percentages (measured with a Karl Fischer laboratory equipment). (The calibration procedure must be performed for each different sensor, as there might occur some mechanical differences (such as the diameter of the metallic disks and the spaces between them) from one sensor to another.) Table 1 presents the data obtained in the calibration of the prototype of the capacitive sensor (used in the experimental results of Section 3). The left column contains the water concentrations (obtained in a Karl Fischer laboratory test) of the four fuel samples (prepared with pure HFO mixed with different amounts of water), and the right column contains the respective capacitances measured with the sensor.
Through a linear interpolation of the data of Table 1, Equation (1) is obtained. This equation relates a given measurement of capacitance (in pF) to the water content (in %).
f w a t e r % = 0.0087 · C p F 2.0712 ,
where C p F is the measured capacitance value (in pF) and f w a t e r % is the estimated water content (in%).

2.3. Estimation of Density and Flow through Ultrasonic System

2.3.1. Density

According to [42,43], ultrasonic methods can determine the density ρ u of a fluid through the use of the Newton–Laplace Equation (2) and measurement of the speed of sound in that medium.
ρ u = 1 κ s · c 2 = B c 2 ,
where κ s is the isentropic (adiabatic) compressibility of the fluid and c is the speed of sound at this fluid. Usually, the compressibility is also expressed in terms of its reciprocal, i.e., the bulk modulus B = 1 / κ s [32,44].
The speed of sound can be obtained either through a pulse-echo system [45] or through a system of two transducers (emitter–receiver) [46,47]. Both principles are based on the measurement of the wave propagation time in the fluid, differing only in the number of ultrasonic transducers used.
In case of a pulse-echo system, the basic principle involves the use of an ultrasonic transducer that outputs a very short pulsed signal of the proper ultrasonic frequency. The ultrasonic wave travels through the propagation medium to the other side, where a reflector material is positioned. Hence, there is a reflection of this ultrasonic wave, causing a large part of the energy of the ultrasonic wave to return through the medium to the transducer that emitted it.
In case of an emitter–receiver system, there are two transducers positioned on separate sides forming an acoustic trajectory between them. While one ultrasonic transducer emits, the other receives the ultrasonic signal that has propagated through the measurement medium.
In the proposed system, multiple acoustic trajectories are used, as illustrated in Figure 4. In the figure, two transducers are opposed at the conduit and form an acoustic trajectory of length L and angle φ with the conduit wall. During transmission, the ultrasonic pulse that travels in favor of the fluid flow through the conduit (in downstream direction) travels the distance L in a time period t d o w n smaller than the pulse that travels in the opposite direction (against the fluid flow, in upstream direction) in a time t u p .
Considering that, in Figure 4, v a x is the axial velocity of the fluid flowing through the conduit (whose amplitude has a value of v ¯ a x ), its projection on the trajectory L (at an angle φ with the flow direction) is given by v ¯ a x · cos φ . Hence, whenever there is a fluid flowing through the conduit, the speed of the ultrasonic pulses downstream ( v d o w n ) and upstream ( v u p ) can be calculated as (3) and (4), respectively.
v d o w n = c + v ¯ a x · cos φ ,
v u p = c v ¯ a x · cos φ .
In addition, the propagation times of the ultrasonic pulses downstream ( t d o w n ) and upstream ( t u p ) can be calculated as (5) and (6), respectively.
t d o w n = L c + v ¯ a x · cos φ ,
t u p = L c v ¯ a x · cos φ .
Combining Equations (5) and (6), it can be shown that the speed of ultrasonic wave c can be obtained as a function of the downstream and upstream travel times as per Equation (7).
c = L 2 · t u p + t d o w n t u p · t d o w n .
As presented in Figure 1, the ultrasonic system has three layers of planes equivalent to Figure 4, and each layer has two crossed ultrasonic trajectories. Hence, Equation (7) can be modified to (8), in order to obtain an average velocity c ¯ considering each of the six trajectories L i and each of the six pairs of downstream and upstream times ( t d o w n , i and t u p , i , respectively).
c ¯ = 1 6 i = 1 6 L i 2 · t u p , i + t d o w n , i t u p , i · t d o w n , i ,
where i is an index corresponding to one of the six trajectories.
The technique used to measure the ultrasonic wave propagation times is based on the zero-crossing detection [48] of the pulses at the receivers. Figure 5 presents the detection of the zero-crossing of an ultrasonic wave. After the transmission of the ultrasonic pulse, the FPGA starts counting the time of several zero-crossings at the receiver. In order to avoid any initial transient, the first cycles are skipped. This is performed by detecting a pre-defined threshold (indicated by (1) in the figure). After the threshold, more three zero-crossings are skipped (until the first valid zero-crossing, indicated by (2) at the figure). Then the time from the starting of the pulse until next M number of zero-crossings ( t T o F 1 until t T o F M ) are computed. The average time of the valid zero crossings is calculated (recursively) using Equation (9).
t ¯ T o F = 1 M · t T o F 1 + j = 2 M t T o F j ( j 1 ) · t ¯ c / 2 ,
where t ¯ T o F is the average of the zero crossing times. It is important to note that t ¯ T o F can be any of the t u p , i or t d o w n , i times in Equation (8). Furthermore, t ¯ c / 2 is the average time of the half cycle of oscillation of the signal, and is calculated with (10).
t ¯ c / 2 = 1 M 1 · j = 1 M 1 t T o F j + 1 t T o F j .
Finally, the average speed of the ultrasound waves (calculated as (8)) at the fluid is used in Equation (2) in order to obtain the density ρ u of the fluid.

2.3.2. Flow Rate

According with [49,50], the flow rate Q v using multiple acoustic trajectories arranged at different heights in the cross section of the conduit can be calculated as (11). This calculation performs a numerical integration of N trajectories arranged in distinct heights in the cross section of the penstock. The most used integration methods [51] are: Gauss–Legendre (for rectangular sections) and Gauss–Jacobi (for circular sections). The OWICS (Optimal Weighted Integration, for circular sections) and the OWIRS (Optimal Weighted Integration for rectangular sections) [48] can also be highlighted. What differentiates the methods are the speed profiles [49,52], whereas the Gauss–Legendre and the Gauss–Jacobi propose uniform speed profiles [51].
Q v = D 2 · i = 1 N w i · v ¯ a x , i ( z i ) · L i · sin ( φ i ) ,
where z i is the position of the trajectory in relation to the cross section of the conduit. v ¯ a x , i ( z i ) is the average axial speed along trajectory i, calculated according to Equation (12). w i are weighting coefficients (which depend on the number of trajectories and on the integration technique in use, according to Equation (13)). D is the dimension of the parallel conduit for the intersection of two acoustic planes (i.e., the diameter of a circular section conduit). N is the total number of acoustic trajectories in a measurement plane. L i is the distance (length) of the acoustic trajectory i, and φ i is the angle between the acoustic trajectories.
v ¯ a x , i ( z i ) = L i 2 · cos ( φ i ) · t u p t d o w n t u p · t d o w n ,
where t u p is the measured propagation time upstream and t d o w n is the measured propagation time downstream. Their difference is known as Δ t = t u p t d o w n .
According to [49,51], the weighting coefficients can be calculated as (13).
w i = 1 1 4 · z i 2 D 2 κ · 2 D · D / 2 D / 2 1 4 · z i 2 D 2 κ · L i . d z ,
where the parameter κ is dependent on the type of integration used ( κ = 0.5 for the Gauss–Legendre method and κ = 0.6 for the OWICS method, considering circular cross-section conduits in both cases). The variable L i is the Lagrange polynomial integration, calculated as (14).
L i z = j = 1 j i N z z j z i z j .
The developed system was based on a measurement configuration using six acoustic trajectories introduced in a circular spool, as presented in Figure 6. The acoustic trajectories are distributed in two crossed measurement planes (A and B), containing three trajectories in each plane. Therefore, the flow rate for a measurement plan (A or B) is determined from Equation (11) and results in (15). Furthermore, each sensor forms an angle of φ i = 45 with the conduit wall.
Q A | B = D 2 · i = 1 3 w i · v ¯ a x , i ( z i ) · L i · sin ( φ i ) .
Using (15) for both planes A and B, the resultant flow rate can be obtained as the average (16).
Q v = Q A + Q B 2 .
Each of the two planes has a trajectory positioned at the center of the circular spool and two other trajectories equidistant from the center, forming the positions z 1 = 0.707107 · D / 2 , z 2 = 0 and z 3 = + 0.707107 · D / 2 . For each trajectory, based on Equations (13) and (14), the respective weights w i are presented in Table 2, either considering Gauss–Jacobi or OWICS integration methods.
Table 3 presents the measured lengths ( L i ) of the trajectories for either planes A and B. Ideally, the lengths L 1 and L 3 of both planes should have the same value, as well as the lengths L 2 of both planes. However, due to a not so precise drilling, there are these 1 mm discrepancies.

2.3.3. Evaluation Setup

In order to evaluate the ultrasonic system, a reduced scale flow setup has been assembled, as presented in Figure 7. This setup is formed by a pump driven by an electric motor (whose speed is controlled by an inverter) and a closed loop hydraulic system.
With the setup of Figure 7, it was possible to vary the flow of the hydraulic circuit by controlling the rotation of the motor that drives the pump. In total, five motor speeds were used (0 RPM, 450 RPM, 500 RPM, 600 RPM, and 700 RPM) in order to achieve a more stable flow in the system. The motor speeds are reduced to the pump using a gear train (of ratio 2.54).
Water has been used as the circulating fluid, as it is it is cheaper than oil and its physical–chemical parameters are well known in the literature. For each motor speed, 100 complete transit time measurements were acquired for each of the 12 trajectories. Means and standard deviations were calculated from the measurements taken, as well as the uncertainties associated with the measurement.
Evaluation of Flow Measurements
Table 4 presents the obtained flows (measured with the ultrasonic system) for each motor rotation (as well as the respective pump rotations).
Ideally, the flow measurements obtained with the proposed sensors should be validated using a certified flow measurement system. As this was not possible, the solution found was to compare the measured flows against the data sheet of the pump [53], which relates flow with pump rotation. From the data of Table 4, a linear regression is performed, resulting in the Equation (17).
Q p r o t o t y p e = 6.91275 × 10 2 × R P M 2.43041 × 10 1 ,
where R P M is the motor rotation (in RPM) and Q p r o t o t y p e is the interpolated flow (in m 3 /h) measured with the prototype.
The manufacturer curve of the pump can be obtained in [53] and results in Equation (18).
Q p u m p = 6.88505 × 10 2 × R P M 2.32743 × 10 1 ,
where Q p u m p is the reference flow (in m 3 /h) for a motor rotation R P M .
Based on Equations (17) and (18), for some values of motor rotation, the flows of the proposed system can be compared against the theoretical flows at the pump. This comparison is presented in Table 5.
Evaluation of Speed-of-Sound and Density Measurements
As the fluid used in the evaluation setup was water, and the fluid temperature at the moment of the tests was 19.40 C, and the value of speed of sound can be obtained as 1480.36 m/s [54] and the value of density can be obtained as 998.34 kg/m 3 [55]. These theoretical values can be used in order to validate the measurements of speed of sound and density performed with the proposed prototype. These comparisons are presented in Table 6 and Table 7.
As observed in Table 6 and Table 7, the deviations are below 0.2%.

2.4. Estimation of Heating Value

The heating values of any organic compound are associated with the bonding energies between the atoms that form the chemical structure of the compound and, therefore, the character of the bonds [56]. However, the possibility of calculation of the heating value for petroleum-derived fuels with reasonable accuracy from elementary composition alone took many researchers to establish empirical correlations from their commonly measured characteristics [57,58]. These correlations are often expressed in the form of linear combinations of the percentages by weight of the elements of the atoms of carbon (C), hydrogen (H), and oxygen (O), and sometimes expanded to sulfur (S) and nitrogen (N). The reason is that the main elements in the chemical composition of fossil fuels are, in fact, limited to C, H, O, N, and S in their organic part [58]. Furthermore, the contaminants present in the oils have a great influence on the efficiency of their combustion. Water is the contaminant that has the greatest influence on HHV and LHV levels, where the greater the amount of water, the lower the values of the two factors will be. Furthermore, sulfur and ashes are two other important contaminants.
The standard ASTM D4868 [13] presents an empirical method to estimate the HHV and the LHV based on the fuel density and the amounts of water, sulfur, and ashes. It is described by Equations (19) and (20).
H H V | A S T M D 4868 = 51.916 8.792 · 10 6 · ρ 2 · 1 f w a t e r + f s u l f u r + f a s h e s + 9.420 · f s u l f u r ,
where ρ is the density of the fuel at 15 C (in kg/m 3 ), obtained from a laboratory test. f w a t e r is the water content (in V / V ), f s u l f u r is the sulfur content (in V / V ), and f a s h is the ash content (in V / V ).
L H V | A S T M D 4868 = 46.423 8.792 · 10 6 · ρ 2 + 3.170 · 10 3 · ρ · 1 f w a t e r + f s u l f u r + f a s h e s + 9.420 · f s u l f u r 2.449 · f w a t e r .
The amount of sulfur can be obtained through laboratory analysis using technologies such as ultraviolet fluorescence, non-dispersive infrared, or X-ray fluorescence spectroscopy [59]. These technologies measure the sulfur content at the molecular level and are very specific and difficult to use online. The same can be said for the ash content. Hence, this work considers only the water content.
The water content is obtained online (with the sensor and method presented in Section 2.2). Furthermore, the density is obtained online (with the sensor and method presented in Section 2.3). These two are applied in simplified versions of Equations (19) and (20), without considering the sulfur and ashes contents, resulting in Equations (21) and (22).
H H V = 51.916 8.792 · 10 6 · ρ 0 2 · 1 f w a t e r ,
L H V = 46.423 8.792 · 10 6 · ρ 0 2 + 3.170 · 10 3 · ρ 0 · 1 f w a t e r 2.449 · f w a t e r ,
where ρ 0 is the density of the fuel compensated in relation to the amount of water using Equation (23).
ρ 0 = ρ u ρ w · f w a t e r 1 f w a t e r ,
where ρ u is the density of the fuel, as obtained with the ultrasonic sensor of Section 2.3. ρ w is the density of water.

3. Results and Discussion

The proposed system has been installed in a Brazilian TPP. Figure 8 presents a photograph of the setup, installed in a section of the pipeline between the TPP fuel storage and the point where the tank truck unloads the fuel oil. The system is composed of one capacitive sensor and six pairs of ultrasonic sensors, as described in Section 2.
Figure 9 presents a photograph of the proposed capacitive sensor for determination of water content, as described in Section 2.2.
Figure 10 presents a photograph of the proposed ultrasonic sensor for determination of density and flow, as described in Section 2.3.
The following results were obtained at three distinct operations of loading, with three different tank trucks. These tests were carried out at different periods due to the requirement of the plant itself for energy generation. During these periods, the measurement system was turned on and followed the supply process carried out by TPP. General loading information is presented in Table 8.
The prototype was configured to generate information in a 20-s update cycle at the time of loading. In these 20 s, acquisitions of signals and data processing were performed. In case of the ultrasonic system (which is responsible for measurement of flow measurement, sound velocity and, consequently, density), 600 UP and DOWN transit time samples were acquired for each acoustic trajectory, and from these measurements, the corresponding calculations were performed. Therefore, each 20 s update presents data for an average of 600 samples of the ultrasonic system. In case of the capacitive system, 5000 capacitance measurement samples were acquired per update cycle that are taken into account to calculate the water content of the emulsion.
The uncertainties presented in the graphs that follows are defined according to the guidelines given in reference [60]. In the case of the ultrasonic system process, the uncertainties are correlated to the measurement of transit times and take into account their mean, standard deviation, and number of samples for each update cycle. The other uncertainties (path length, angle, etc.) related to uncertainty propagation are static, as they involve measurements based on calibrated equipment. The same principle is used in the capacitive system, where the capacitance uncertainty is estimated by the mean, standard deviation, and number of samples of measurements performed in each update cycle. Uncertainty propagation for the other quantities, density, and LHV use the concept of uncertainty propagation [61].

3.1. Flow

Figure 11 and Figure 12 present the results of the flow measurements performed for the loading operations of Table 8. Figure 11 presents the general perspective of the three loading operations. It can be noted that the flow measurements of the three operations are very close over time.
Figure 12 presents a zoom in the region of Figure 11 where the loading is actually happening, in order to visualize more clearly the calculated uncertainties. It is possible to observe that some samples present high levels of uncertainty. For the first loading (blue color) it is observed that the uncertainty levels at the beginning of the loading (from the start of the procedure until near the 15 min mark) had some samples with large levels of uncertainty, while the remaining samples had lower levels. These variations in the uncertainties are related to the dispersion of transit time differences (Equation (12)) with high standard deviations and, consequently, higher levels of uncertainty in the flow.
Table 9 presents the summary of the flow measurements of the three loading operations of Table 8. As noted, the average values of the means are close to each other and have similar standard deviations from the mean. Concerning the uncertainties, it can be seen that the maximum uncertainty observed between loading operations was 6.74%. Variations in uncertainties are related to the standard deviation of the difference in transit time ( Δ t = t u p t d o w n ) of the trajectories.
The average fluid temperatures during the loading operations were 17.75 C, 19.58 C, and 21.55 C, for loading operations 1, 2, and 3, respectively.

3.2. Capacitance and Water Content

Figure 13 presents the results of the flow measurements performed for the loading operations of Table 8. Table 10 presents the summary of the capacitance measurements of the three loading operations.
Based on the measured capacitances, applying Equation (1), the water content for the loading operations of Table 8 can be obtained. Figure 14 presents the estimated water contents. Table 11 presents the summary of the water content measurements of the three loading operations. It can be noted that the average water contents were 0.456%, 0.503%, and 0.514%, for loading operations 1, 2, and 3, respectively.
The results obtained using the capacitive sensor should have, ideally, been compared against a laboratory test. However, the only laboratory results available had been obtained almost three months before the installation of the prototype, but they are used here, in order to give an idea of how close the proposed system compares to a laboratory analysis. The complete laboratory results are presented in the Appendix A. The results obtained with the prototype indicate around 0.5% of water content. The laboratory results indicated 0.05%. This a (somehow) large discrepancy—although both values are well below the limit of 2% determined by the Brazilian National Agency of Petroleum, Natural Gas and Biofuels (ANP) [36].

3.3. Speed of Sound and Density

Figure 15 presents the results of the online measurements of the speed of sound at the emulsion, performed for the loading operations of Table 8. Table 12 presents the summary of the speed of sound measurements of the three loading operations.
Based on the measured speeds of sound, applying Equation (23), the densities for the loading operations of Table 8 can be obtained. Figure 16 presents the estimated densities. Table 13 presents the summary of the density measurements of the three loading operations.
The results obtained with the prototype indicate average density values from 928.2 to 949.2 kg/m 3 . The laboratory results of the Appendix A indicated 927.7 kg/m 3 . This indicates a maximum discrepancy of 2.3%, although the laboratory tests have been performed three months before the installation of the prototype.

3.4. Heating Value

Using the measured values of the water content (data from Figure 14) and the density (data from Figure 16), the lower heating value can be estimated online, as presented in Figure 17.
Table 14 presents the summary of the LHV measurements of the three loading operations. It can be observed that the calculated uncertainties had values below 0.25%.
The results obtained with prototype indicate average LHV values from 41.31 to 41.57 MJ/kg. The laboratory results of the Appendix A indicate 41.59 MJ/kg. This indicates a discrepancy below 1%, although the laboratory tests have been performed three months before the installation of the prototype.

4. Conclusions

This paper presented a hybrid online monitoring system for the determination of fuel oil quality. It is proposed the association of two different measurement techniques (i.e., capacitive and ultrasonic techniques) and the correlation of their information in order to assess the quality parameters. The capacitive technique aims to measure the water content in the fuel oil, while the ultrasonic technique directly measures the fuel’s density and flow. The combination of information from both techniques serves to infer about other quality indicators and carry out calibration processes.
The estimation of water content through the capacitive technique is based on the fact that the dielectric properties of the fuel oil changes according to its water content. Hence, a capacitive sensor (composed of multiple metallic disks) is proposed in order to measure the capacitance of the emulsion (fuel oil with water contamination) that fills the spaces between the disks.
The estimation of density and flow through the ultrasonic technique is based on the propagation times of the ultrasonic pulses on the fuel oil. A system composed of an array of sensors is proposed in order to estimate the velocity of the pulses, even with the fuel in motion through the pipeline where the sensors are installed.
With the information of both the water content and the density, the heating value of the fuel oil can be estimated. The heating value is an important parameter that shows the quality of the fuel, which implies in the efficiency of the combustion cycle of the power generators. Usually, combustion-based power plants have to perform laboratory tests in order to determine the quality of their fuels. In these cases, samples are collected and sent to a laboratory for the analysis. With the proposed system, the analysis can be performed instantaneously, while the tank truck is loading the plant.
The proposed system, at the current stage of research, is an invasive procedure. As a future work, it is being studied a non invasive procedure, for both the capacitive and ultrasonic systems. Furthermore, as future work, it is planned to analyze the influence of flow conditions on the measurements of capacitance and the effects of the temperature on the density of the fuel oil. To date, due to intermittent loading operations at the TPP (few operations per month due to low demand), some long-term measurements, such as the evaluation of drift in the sensors, cannot be performed.

Author Contributions

Conceptualization, M.M.C., L.E.B.-d.-S., D.d.A.A., C.E.T. and E.L.B.; data curation, L.E.B.-d.-S., C.E.T., E.L.B., G.L.-T. and L.E.L.O.; formal analysis, M.M.C., D.d.A.A., C.E.T., R.F.R.J. and G.P.K.; investigation, M.M.C., D.d.A.A., C.E.T., L.E.B.-d.-S., R.F.R.J. and G.P.K.; methodology and software, M.M.C., D.d.A.A., C.E.T., G.P.K., L.E.L.O. and E.L.B.; project administration, L.E.B.-d.-S., C.E.T., E.L.B., G.L.-T., L.E.L.O. and R.G.d.P.; validation, M.M.C., D.d.A.A., C.E.T. and R.F.R.J.; supervision, L.E.B.-d.-S., C.E.T., E.L.B., G.L.-T., L.E.L.O. and R.G.d.P.; writing, M.M.C. and W.C.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Norte Energia S.A. and the Brazilian Research Agencies CNPq, CAPES, and ANEEL R&D for their support of this project.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
FOFuel Oil
GHCGross Heat of Combustion
HFOHeavy Fuel Oil
HHVHigher Heating Value
LFOLight Fuel Oil
LHVLower Heating Value
NHCNet Heat of Combustion
TPPThermal Power Plant

Appendix A. Laboratory Tests Performed on a Sample of HFO

Table A1 presents an excerpt of the results sheet of a laboratory analysis performed on a sample of HFO of the same TPP where the prototype has been installed. The test, however, has been performed some months before the installation of the system.
Table A1. Results of laboratory tests performed on a sample of HFO (from the results sheet provided by the Brazilian laboratory Intertek).
Table A1. Results of laboratory tests performed on a sample of HFO (from the results sheet provided by the Brazilian laboratory Intertek).
MethodsTestsResultsUnits
ASTM D1298density (20/4 C)927.7kg/m 3
ASTM D95water by distillation0.05% of volume
ASTM D4868lower heating value41.59kJ/kg

References

  1. Santos, A.Q.O.; da Silva, A.R.; Ledesma, J.J.G.; de Almeida, A.B.; Cavallari, M.R.; Junior, O.H.A. Electricity Market in Brazil: A Critical Review on the Ongoing Reform. Energies 2021, 14, 2873. [Google Scholar] [CrossRef]
  2. ANEEL. Sistemas de Geracao da ANEEL—SIGA. [S.1.] 2021. Available online: http://www.aneel.gov.br/siga (accessed on 17 September 2021).
  3. Ehsan, M.; Rahman, M.M.; Saadi, H. Effect of fuel adulteration on engine crankcase dilution. J. Mech. Eng. 2010, 41, 114–120. [Google Scholar] [CrossRef] [Green Version]
  4. Boadu, K. Effects of adulteration on diesel oil with kerosene fuel in Ghana. J. Appl. Sci. Environ. Manag. 2019, 23, 1195–1200. [Google Scholar] [CrossRef]
  5. Gawande, A.P.; Kaware, J.P. Fuel adulteration consequences in India: A review. Sci. Rev. Chem. Commun. 2013, 3, 161–171. [Google Scholar]
  6. Massawe, E.A.; Kilavo, H.; Sam, A. Fuel Adulteration in Tanzania and its Consequencies: An Overview. Res. J. Eng. Appl. Sci. 2013, 4, 8. [Google Scholar]
  7. Mazivila, S.J.; Gontijo, L.C.; Santana, F.B.D.; Mitsutake, H.; Santos, D.Q.; Borges Neto, W. Fast Detection of Adulterants/Contaminants in Biodiesel/Diesel Blend (B5) Employing Mid-Infrared Spectroscopy and PLS-DA. Energy Fuels 2015, 29, 227–232. [Google Scholar] [CrossRef]
  8. Vilas Boas, F.M.; Borges-da Silva, L.E.; Villa-Nova, H.F.; Bonaldi, E.L.; Oliveira, L.E.L.; Lambert-Torres, G.; Assuncao, F.D.O.; Costa, C.I.D.A.; Campos, M.M.; Sant’Ana, W.C.; et al. Condition Monitoring of Internal Combustion Engines in Thermal Power Plants Based on Control Charts and Adapted Nelson Rules. Energies 2021, 14, 4924. [Google Scholar] [CrossRef]
  9. Assuncao, F.D.O.; Borges-da Silva, L.E.; Villa-Nova, H.F.; Bonaldi, E.L.; Oliveira, L.E.L.; Lambert-Torres, G.; Teixeira, C.E.; Sant’Ana, W.C.; Lacerda, J.; Silva Junior, J.L.M.D.; et al. Reduced Scale Laboratory for Training and Research in Condition-Based Maintenance Strategies for Combustion Engine Power Plants and a Novel Method for Monitoring of Inlet and Exhaust Valves. Energies 2021, 14, 6298. [Google Scholar] [CrossRef]
  10. Mendonca, P.; Bonaldi, E.; de Oliveira, L.; Lambert-Torres, G.; da Silva, J.B.; da Silva, L.B.; Salomon, C.; Santana, W.; Shinohara, A. Detection and modelling of incipient failures in internal combustion engine driven generators using Electrical Signature Analysis. Electr. Power Syst. Res. 2017, 149, 30–45. [Google Scholar] [CrossRef]
  11. American Society for Testing and Materials. ASTM D95-13 (2018)—Standard Test Method for Water in Petroleum Products and Bituminous Materials by Distillation; ASTM International: West Conshohocken, PA, USA, 2018. [Google Scholar]
  12. American Society for Testing and Materials. ASTM D1298-12b (2017)—8 ASTM. ASTM D1298-12b(2017)—Standard Test Method for Density, Relative Density, or API Gravity of Crude Petroleum and Liquid Petroleum Products by Hydrometer Method; ASTM International: West Conshohocken, PA, USA, 2017. [Google Scholar]
  13. American Society for Testing and Materials. ASTM D4868-00 (2010)—Standard Test Method for Estimation of Net and Gross Heat of Combustion of Burner and Diesel Fuels; ASTM International: West Conshohocken, PA, USA, 2010. [Google Scholar]
  14. Han, Y.; Zhao, A.; Zhang, H.; Ren, Y.; Liu, W.; Jin, N. Differential pressure method for measuring water holdup of oil–water two-phase flow with low velocity and high water-cut. Exp. Therm. Fluid Sci. 2016, 72, 197–209. [Google Scholar] [CrossRef]
  15. Al-Saiyed, M.R.; Warsi, S.; Phillips, J.; Gilani, S.; Zareef, M.; Lievois, J. Infrared Water Cut Meter for Cost Effective Individual Monitoring of Complex Wells. In Proceedings of the 26th International North Sea Flow Measurement Workshop, St. Andrews, Scotland, 21–24 October 2008. [Google Scholar]
  16. Scheers, A. An Oil/Water/Gas Composition Meter Based on Multiple Energy Gamma Ray Absorption (MEGRA) Measurement; International Atomic Energy Agency (IAEA): Vienna, Austria, 2000; pp. 105–122. [Google Scholar]
  17. Johansen, G.A.; Tjugum, S.A. Fluid composition analysis by multiple gamma-ray beam and modality measurements. In Proceedings of the 25th International North Sea Flow Measurement Workshop, NFSOGM—Norwegian Society for Oil and Gas Measurement, Oslo, Norway, 16–18 October 2007. [Google Scholar]
  18. Rigaku Applied Technologies, Inc. NEX OL Process Elemental Analyzer for Liquid Streams—On-line, Real-Time Analysis by Energy Dispersive X-ray Fluorescence (EDXRF)—Online Catalog. 2021. Available online: https://rigakuedxrf.com/nex-ol.php?gclid=Cj0KCQjw7pKFBhDUARIsAFUoMDZ9rZZZnKoF9lLcjd8FPYUlaNiSO7OXMWHKJ57xc1WCqM1RNT3W4RgaAopxEALw_wcB (accessed on 17 September 2021).
  19. Gregory, A.; Clarke, R. A review of RF and microwave techniques for dielectric measurements on polar liquids. IEEE Trans. Dielectr. Electr. Insul. 2006, 13, 727–743. [Google Scholar] [CrossRef]
  20. Makeev, Y.V.; Lifanov, A.; Sovlukov, A.S. Microwave measurement of water content in flowing crude oil. Autom. Remote Control 2013, 74, 157–169. [Google Scholar] [CrossRef]
  21. Teixeira, C.E.; da Silva, L.E.B.; Veloso, G.F.; Lambert-Torres, G.; Campos, M.M.; Noronha, I.; Bonaldi, E.L.; de Oliveira, L.E.L. An ultrasound-based water-cut meter for heavy fuel oil. Measurement 2019, 148, 106907. [Google Scholar] [CrossRef]
  22. Hammer, E.; Tollefsen, J.; Olsvik, K. Capacitance transducers for non-intrusive measurement of water in crude oil. Flow Meas. Instrum. 1989, 1, 51–58. [Google Scholar] [CrossRef]
  23. Schuller, R.; Gundersen, T.; Halleraker, M.; Engebretsen, B. Measurement of water concentration in oil/water dispersions with a circular single-electrode capacitance probe. IEEE Trans. Instrum. Meas. 2004, 53, 1378–1383. [Google Scholar] [CrossRef]
  24. Libert, N.; Morales, R.E.; da Silva, M.J. Capacitive measuring system for two-phase flow monitoring. Part 1: Hardware design and evaluation. Flow Meas. Instrum. 2016, 47, 90–99. [Google Scholar] [CrossRef]
  25. Libert, N.; Morales, R.E.; da Silva, M.J. Capacitive measuring system for two-phase flow monitoring. Part 2: Simulation-based calibration. Flow Meas. Instrum. 2016, 50, 102–111. [Google Scholar] [CrossRef]
  26. Hansen, L.S.; Pedersen, S.; Durdevic, P. Multi-Phase Flow Metering in Offshore Oil and Gas Transportation Pipelines: Trends and Perspectives. Sensors 2019, 19, 2184. [Google Scholar] [CrossRef] [Green Version]
  27. Wang, F.; Marashdeh, Q.; Fan, L.S.; Warsito, W. Electrical Capacitance Volume Tomography: Design and Applications. Sensors 2010, 10, 1890–1917. [Google Scholar] [CrossRef]
  28. Mosorov, V.; Zych, M.; Hanus, R.; Sankowski, D.; Saoud, A. Improvement of Flow Velocity Measurement Algorithms Based on Correlation Function and Twin Plane Electrical Capacitance Tomography. Sensors 2020, 20, 306. [Google Scholar] [CrossRef] [Green Version]
  29. De Almeida Arantes, D.; Borges da Silva, L.E.; Teixeira, C.E.; Campos, M.M.; Lambert-Torres, G.; Bonaldi, E.L.; de Lacerda de Oliveira, L.E.; da Costa, G.A. Relative Permittivity Meter Using a Capacitive Sensor and an Oscillating Current Source. In Proceedings of the IECON 2019—45th Annual Conference of the IEEE Industrial Electronics Society, Lisbon, Portugal, 14–17 October 2019; Volume 1, pp. 806–811. [Google Scholar] [CrossRef]
  30. Bjørndal, E. Acoustic Measurement of Liquid Density with Applications for Mass Measurement of Oil. Ph.D. Thesis, University of Bergen, Bergen, Norway, 2007. [Google Scholar]
  31. Bjorndal, E.; Froysa, K.E. Acoustic methods for obtaining the pressure reflection coefficient from a buffer rod based measurement cell. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2008, 55, 1781–1793. [Google Scholar] [CrossRef] [Green Version]
  32. Hoche, S.; Hussein, M.; Becker, T. Ultrasound-based density determination via buffer rod techniques: A review. J. Sens. Sens. Syst. 2013, 2, 103–125. [Google Scholar] [CrossRef] [Green Version]
  33. Kim, J.O.; Bau, H.H. Instrument for simultaneous measurement of density and viscosity. Rev. Sci. Instrum. 1989, 60, 1111–1115. [Google Scholar] [CrossRef] [Green Version]
  34. Swoboda, C.A.; Cannon, P.H.; Fredrickson, D.R.; Gabelnick, S.D.; Hornstra, F.; Phan, K.A.; Yao, N.P. Development of an ultrasonic technique to measure specfic gravity in lead-acid battery electrolyte. IEEE Trans. Sonics Ultrason. (USA) 1983, 30, 69–77. [Google Scholar] [CrossRef]
  35. Davis, L.A.; Gordon, R.B. Compression of mercury at high pressure. J. Chem. Phys. 1967, 46, 2650–2660. [Google Scholar] [CrossRef]
  36. Agencia Nacional do Petroleo, Gas Natural e Biocombustiveis. Resolucao ANP nº 3 DE 27/01/2016. 2016. Available online: https://www.gov.br/anp/pt-br/assuntos/producao-de-derivados-de-petroleo-e-processamento-de-gas-natural/producao-de-derivados-de-petroleo-e-processamento-de-gas-natural/oleo-combustivel (accessed on 28 September 2021).
  37. Speight, J.G. The Chemistry and Technology of Petroleum; CRC Press: Boca Raton, FL, USA, 2014. [Google Scholar]
  38. Al-Ghouti, M.; Al-Degs, Y.; Mustafa, F. Determination of hydrogen content, gross heat of combustion, and net heat of combustion of diesel fuel using FTIR spectroscopy and multivariate calibration. Fuel 2010, 89, 193–201. [Google Scholar] [CrossRef]
  39. Demirbas, A.; Ak, N.; Aslan, A.; Sen, N. Calculation of higher heating values of hydrocarbon compounds and fatty acids. Pet. Sci. Technol. 2018, 36, 712–717. [Google Scholar] [CrossRef]
  40. Malmberg, C.; Maryott, A. Dielectric constant of water from 0 to 100 C. J. Res. Natl. Bur. Stand. 1956, 56, 1–8. [Google Scholar] [CrossRef]
  41. Warren, W. BS and W Measurement Principles and Practices. J. Pet. Technol. 1962, 14, 1207–1212. [Google Scholar] [CrossRef]
  42. Kang, S.; You, Y.; Feng, M. Study of high-pressure physical properties of marine heavy oil. In Key Engineering Materials; Trans Tech Publ.: Bach, Switzerland, 2016; Volume 693, pp. 411–418. [Google Scholar]
  43. Akahama, Y.; Kawaguchi, S.; Hirao, N.; Ohishi, Y. Elastic softening of bulk modulus of monoclinic HfO2 under high pressure. Appl. Phys. Lett. 2020, 117, 182903. [Google Scholar] [CrossRef]
  44. Boehman, A.L.; Morris, D.; Szybist, J.; Esen, E. The impact of the bulk modulus of diesel fuels on fuel injection timing. Energy Fuels 2004, 18, 1877–1882. [Google Scholar] [CrossRef]
  45. McClements, D.; Fairley, P. Ultrasonic pulse echo reflectometer. Ultrasonics 1991, 29, 58–62. [Google Scholar] [CrossRef]
  46. Hale, J. Ultrasonic density measurement for process control. Ultrasonics 1988, 26, 356–357. [Google Scholar] [CrossRef]
  47. Henning, B.; Daur, P.C.; Prange, S.; Dierks, K.; Hauptmann, P. In-line concentration measurement in complex liquids using ultrasonic sensors. Ultrasonics 2000, 38, 799–803. [Google Scholar] [CrossRef]
  48. Campos, M.M. Medidor de Vazao de Multi-Trajetorias (in Portuguese). Master’s Thesis, Universidade Federal de Itajuba, Itajuba, Brazil, 2017. [Google Scholar]
  49. Gruber, P.; Wermelinger, F.; Hug, S. Experience with the ATT method for discharge measurement with adaptive weighting of the path velocities. In Proceedings of the HYDRO 2013 Conference, Innsbruck, Austria, 7–9 October 2013. [Google Scholar]
  50. Gruber, P.; Marushchenko, S. Comparative study of 4(8)-path and 5(10)-path configurations for ATT flow measurements in circular conduits. In Proceedings of the IGHEM 2014 Conference, Itajuba, Brazil, 16–19 September 2014. [Google Scholar]
  51. Tresch, T.; Lüscher, B.; Staubli, T.; Gruber, P. Presentation of optimized integration methods and weighting corrections for the acoustic discharge measurement. In Proceedings of the IGHEM 2008 Conference, Milano, Italy, 3–6 September 2008. [Google Scholar]
  52. Hug, S.; Staubli, T.; Gruber, P. Comparison of measured path velocities with numerical simulations for heavily disturbed velocity distributions. In Proceedings of the IGHEM 2012 Conference, Trondheim, Norway, 27–30 June 2012. [Google Scholar]
  53. RZR Bombas Positivas. Bombas de Engrenagens Internas Modelo RZR-30000. 2014. Available online: https://rzrbombas.com.br/wp-content/themes/rzr/catalogos/RZR-30.000-SVA.pdf (accessed on 31 October 2021).
  54. The Engineering ToolBox. Water—Speed of Sound. 2021. Available online: https://www.engineeringtoolbox.com/sound-speed-water-d_598.html (accessed on 31 October 2021).
  55. Jones, F.E.; Harris, G.L. ITS-90 density of water formulation for volumetric standards calibration. J. Res. Natl. Inst. Stand. Technol. 1992, 97, 335. [Google Scholar] [CrossRef]
  56. Buck, W.; Rudtsch, S. Thermal Properties. In Springer Handbook of Materials Measurement Methods; Czichos, H., Saito, T., Smith, L., Eds.; Springer: Berlin/Heidelberg, Germany, 2006; pp. 399–429. [Google Scholar] [CrossRef]
  57. Marlair, G.; Cwiklinski, C.; Tewarson, A. An analysis of some practical methods for estimating heats of combustion in fire safety studies. In Proceedings of the Interflam 99 Conference, Edinburgh, UK, 29 June–1 July 1999. [Google Scholar]
  58. Ponomarev, D.A.; Mettee, H.D.; Miller, J. Empirically estimated heats of combustion of oxygenated hydrocarbon bio-type oils. BioResources 2015, 10, 3638–3656. [Google Scholar] [CrossRef] [Green Version]
  59. Corilo, Y.E.; Rowland, S.M.; Rodgers, R.P. Calculation of the total sulfur content in crude oils by positive-ion atmospheric pressure photoionization Fourier transform ion cyclotron resonance mass spectrometry. Energy Fuels 2016, 30, 3962–3966. [Google Scholar] [CrossRef]
  60. BIPM; IEC; IFCC; ILAC; ISO; IUPAC; IUPAP; OIML. Evaluation of Measurement Data—Guide to the Expression of Uncertainty in Measurement. 2008. Available online: https://www.bipm.org/utils/common/documents/jcgm/JCGM_100_2008_E.pdf (accessed on 30 October 2021).
  61. BRITISH STANDARD. BS ISO 5168:2005—Measurement of Fluid Flow—Procedures for the Evaluation of Uncertainties; International Organization for Standardization (ISO): Geneva, Switzerland, 2005. [Google Scholar]
Figure 1. Conceptual 3D drawing of the proposed system.
Figure 1. Conceptual 3D drawing of the proposed system.
Sensors 21 07979 g001
Figure 2. Proposed capacitive sensor and its equivalent circuit.
Figure 2. Proposed capacitive sensor and its equivalent circuit.
Sensors 21 07979 g002
Figure 3. Proposed capacitive sensor installed in a transversal section of the pipeline.
Figure 3. Proposed capacitive sensor installed in a transversal section of the pipeline.
Sensors 21 07979 g003
Figure 4. Measurement of time-of-flight in a conduit using multiple acoustic trajectories of ultrasonic transducers.
Figure 4. Measurement of time-of-flight in a conduit using multiple acoustic trajectories of ultrasonic transducers.
Sensors 21 07979 g004
Figure 5. Measurement of the transit time through multiple zero crossings.
Figure 5. Measurement of the transit time through multiple zero crossings.
Sensors 21 07979 g005
Figure 6. Ultrasonic sensors forming two planes of three trajectories.
Figure 6. Ultrasonic sensors forming two planes of three trajectories.
Sensors 21 07979 g006
Figure 7. Evaluation setup for the ultrasonic sensors.
Figure 7. Evaluation setup for the ultrasonic sensors.
Sensors 21 07979 g007
Figure 8. Proposed system installed at the fuel input of a Brazilian TPP.
Figure 8. Proposed system installed at the fuel input of a Brazilian TPP.
Sensors 21 07979 g008
Figure 9. Proposed capacitive sensor for determination of water content.
Figure 9. Proposed capacitive sensor for determination of water content.
Sensors 21 07979 g009
Figure 10. Proposed ultrasonic sensor for determination of density and flow.
Figure 10. Proposed ultrasonic sensor for determination of density and flow.
Sensors 21 07979 g010
Figure 11. Measurement of flow using the proposed ultrasonic sensor of Section 2.3 for the three loading operations of Table 8.
Figure 11. Measurement of flow using the proposed ultrasonic sensor of Section 2.3 for the three loading operations of Table 8.
Sensors 21 07979 g011
Figure 12. Zoomed view on the measurements of flow, in order to better visualize the uncertainties.
Figure 12. Zoomed view on the measurements of flow, in order to better visualize the uncertainties.
Sensors 21 07979 g012
Figure 13. Measurement of capacitance using the proposed capacitive system of Section 2.2 for the three loading operations of Table 8.
Figure 13. Measurement of capacitance using the proposed capacitive system of Section 2.2 for the three loading operations of Table 8.
Sensors 21 07979 g013
Figure 14. Measurement of water content using the proposed capacitive system of Section 2.2 for the three loading operations of Table 8.
Figure 14. Measurement of water content using the proposed capacitive system of Section 2.2 for the three loading operations of Table 8.
Sensors 21 07979 g014
Figure 15. Measurement of speed of sound using the proposed ultrasonic system of Section 2.3 for the three loading operations of Table 8.
Figure 15. Measurement of speed of sound using the proposed ultrasonic system of Section 2.3 for the three loading operations of Table 8.
Sensors 21 07979 g015
Figure 16. Measurement of density using the proposed ultrasonic system of Section 2.3 for the three loading operations of Table 8.
Figure 16. Measurement of density using the proposed ultrasonic system of Section 2.3 for the three loading operations of Table 8.
Sensors 21 07979 g016
Figure 17. Measurement of Lower Heating Value using both the proposed capacitive system of Section 2.2 and the proposed ultrasonic system of Section 2.3 for the three loading operations of Table 8.
Figure 17. Measurement of Lower Heating Value using both the proposed capacitive system of Section 2.2 and the proposed ultrasonic system of Section 2.3 for the three loading operations of Table 8.
Sensors 21 07979 g017
Table 1. Calibration example of the capacitive sensor.
Table 1. Calibration example of the capacitive sensor.
Karl Fischer Water ContentMeasured Capacitance
0.19%261.8 pF
0.40%286.1 pF
1.01%353.2 pF
2.20%482.9 pF
Table 2. Weighting coefficients for a circular cross section conduit.
Table 2. Weighting coefficients for a circular cross section conduit.
Trajectory iPosition z i / ( D / 2 ) w i for Gauss–Jacobi w i for OWICS
1−0.7071070.5553600.546150
200.7853980.792715
3+0.7071070.5553600.546150
Table 3. Lengths of the acoustic trajectories.
Table 3. Lengths of the acoustic trajectories.
Trajectory iLength L i for Plane ALength L i for Plane B
1162.26 mm163.54 mm
2203.00 mm202.83 mm
3163.05 mm162.93 mm
Table 4. Relation between flow and motor/pump angular speed.
Table 4. Relation between flow and motor/pump angular speed.
Motor RotationPump RotationFlow
0 RPM0.00 RPM0.00 m 3 /h
450 RPM117.17 RPM11.47 m 3 /h
500 RPM196.85 RPM13.19 m 3 /h
600 RPM236.22 RPM16.07 m 3 /h
700 RPM275.59 RPM19.29 m 3 /h
Table 5. Comparison between the theoretical flows against the interpolated flows of the prototype.
Table 5. Comparison between the theoretical flows against the interpolated flows of the prototype.
Motor Rotation Q pump Q prototype Deviation
400 RPM27.41 m 3 /h27.31 m 3 /h0.37%
500 RPM34.32 m 3 /h34.19 m 3 /h0.37%
600 RPM41.23 m 3 /h41.08 m 3 /h0.38%
700 RPM48.15 m 3 /h47.96 m 3 /h0.38%
Table 6. Validation of measurements of speed of sound of water circulating in the evaluation setup at 19.40 C, whose theoretical value is 1480.36 m/s [54].
Table 6. Validation of measurements of speed of sound of water circulating in the evaluation setup at 19.40 C, whose theoretical value is 1480.36 m/s [54].
Motor RotationSpeed of SoundDeviation
0 RPM1479.01 m 3 /h0.091%
450 RPM1480.52 m 3 /h0.011%
500 RPM1480.31 m 3 /h0.004%
600 RPM1480.73 m 3 /h0.025%
700 RPM1481.10 m 3 /h0.034%
Table 7. Validation of measurements of density of water circulating in the evaluation setup at 19.40 C, whose theoretical value is 998.34 kg/m 3 [55].
Table 7. Validation of measurements of density of water circulating in the evaluation setup at 19.40 C, whose theoretical value is 998.34 kg/m 3 [55].
Motor RotationDensityDeviation
0 RPM1000.16 kg/m 3 0.183%
450 RPM998.13 kg/m 3 0.021%
500 RPM998.41 kg/m 3 0.007%
600 RPM997.84 kg/m 3 0.050%
700 RPM997.35 kg/m 3 0.069%
Table 8. Recorded loading operations at the TPP.
Table 8. Recorded loading operations at the TPP.
LoadingDateTimeDuration
19 June 202110:50–11:3545 20
230 June 202110:52–11:3644 00
330 June 202112:05–12:5045 20
Table 9. Summary of the flow measurements of the three loading operations of Table 8.
Table 9. Summary of the flow measurements of the three loading operations of Table 8.
LoadingAverage of the MeansStandard Deviation of the MeansMinimum UncertaintyMaximum Uncertainty
161.3 m 3 /h0.3 m 3 /h0.23%6.74%
260.8 m 3 /h0.3 m 3 /h0.23%5.93%
360.7 m 3 /h0.3 m 3 /h0.23%1.31%
Table 10. Summary of the capacitance measurements of the three loading operations of Table 8.
Table 10. Summary of the capacitance measurements of the three loading operations of Table 8.
LoadingAverage of the MeansStandard Deviation of the MeansMinimum UncertaintyMaximum Uncertainty
1290.5 pF0.3 pF0.032%0.054%
2297.4 pF0.5 pF0.028%0.152%
3295.9 pF0.3 pF0.024%0.048%
Table 11. Summary of the water content measurements of the three loading operations of Table 8.
Table 11. Summary of the water content measurements of the three loading operations of Table 8.
LoadingAverage of the MeansStandard Deviation of the MeansMinimum UncertaintyMaximum Uncertainty
10.456%0.002%0.18%0.30%
20.503%0.003%0.14%0.77%
30.514%0.003%0.12%0.25%
Table 12. Summary of the speed of sound measurements of the three loading operations of Table 8.
Table 12. Summary of the speed of sound measurements of the three loading operations of Table 8.
LoadingAverage of the MeansStandard deviation of the MeansMinimum UncertaintyMaximum Uncertainty
11352.4 m/s0.5 m/s0.122%0.122%
21367.6 m/s1.2 m/s0.121%0.142%
31363.5 m/s1.3 m/s0.121%0.122%
Table 13. Summary of the density measurements of the three loading operations of Table 8.
Table 13. Summary of the density measurements of the three loading operations of Table 8.
LoadingAverage of the MeansStandard Deviation of the MeansMinimum UncertaintyMaximum Uncertainty
1949.2 kg/m 3 0.6 kg/m 3 0.556%0.556%
2928.2 kg/m 3 1.6 kg/m 3 0.556%0.575%
3933.8 kg/m 3 1.7 kg/m 3 0.556%0.566%
Table 14. Summary of the Lower Heating Values measurements of the three loading operations of Table 8.
Table 14. Summary of the Lower Heating Values measurements of the three loading operations of Table 8.
LoadingAverage of the MeansStandard Deviation of the MeansMinimum UncertaintyMaximum Uncertainty
141.31 MJ/kg0.01 MJ/kg0.192%0.222%
241.57 MJ/kg0.02 MJ/kg0.181%0.210%
341.49 MJ/kg0.02 MJ/kg0.180%0.211%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Campos, M.M.; Borges-da-Silva, L.E.; Arantes, D.d.A.; Teixeira, C.E.; Bonaldi, E.L.; Lambert-Torres, G.; Ribeiro Junior, R.F.; Krupa, G.P.; Sant’Ana, W.C.; Oliveira, L.E.L.; et al. An Ultrasonic-Capacitive System for Online Characterization of Fuel Oils in Thermal Power Plants. Sensors 2021, 21, 7979. https://doi.org/10.3390/s21237979

AMA Style

Campos MM, Borges-da-Silva LE, Arantes DdA, Teixeira CE, Bonaldi EL, Lambert-Torres G, Ribeiro Junior RF, Krupa GP, Sant’Ana WC, Oliveira LEL, et al. An Ultrasonic-Capacitive System for Online Characterization of Fuel Oils in Thermal Power Plants. Sensors. 2021; 21(23):7979. https://doi.org/10.3390/s21237979

Chicago/Turabian Style

Campos, Mateus Mendes, Luiz Eduardo Borges-da-Silva, Daniel de Almeida Arantes, Carlos Eduardo Teixeira, Erik Leandro Bonaldi, Germano Lambert-Torres, Ronny Francis Ribeiro Junior, Gabriel Pedro Krupa, Wilson Cesar Sant’Ana, Levy Ely Lacerda Oliveira, and et al. 2021. "An Ultrasonic-Capacitive System for Online Characterization of Fuel Oils in Thermal Power Plants" Sensors 21, no. 23: 7979. https://doi.org/10.3390/s21237979

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