Capacitive Water-Cut Meter with Robust Near-Linear Transfer Function

: The water content in fuel–water emulsions can vary from 10% to 30%, and is under control during the process of emulsiﬁcation. The main task of this study was to obtain near-linear static function for a water-cut meter with capacitive sensors, and to provide it with effective type-uncertainty compensation during the process of water–fuel emulsion moisture control. To fulﬁll the capacitive measurements, two capacitive sensors in the measuring channel and two capacitive sensors in the reference channel were used. The method of least squares and general linear regression instruments were used to obtain robust and near-linear transfer function of the capacitive water-cut meter. The prototype product of the water-cut meter was developed with the purpose of fulﬁlling multiple moisture measurements and checking the workability of the new transfer function. Values of moisture content for the new transfer function and the closest analog were compared with the help of dispersion analysis. The new transfer function provided minimal dispersions of repeatability and adequacy, and minimal F -test values, proving its better capability for type-uncertainty compensation and better adequacy for the nominal linear transfer function of the water-cut meter.


Introduction
Emulsified fuels have a clear effect on the combustion process, and are one of the modern combustible systems that offer prospects for new technologies. This has been exploited recently in diesel engines, which have been used for a long time as the main source of driving power due to their reliable design and high fuel-saving economy [1,2]. Significant results in the toxicity reduction of diesel engines' exhaust fumes were achieved by the application of fuels with water and alcohol additives [3]. Hence, emulsified fuels are mixtures of a combustible liquid, either oil or fuel as a continuous phase, and a smaller amount of water with or without a surfactant as a dispersed phase. Emulsion of petroleum fuel in water is the most common form of such emulsion [4]. In emulsified diesel fuel, the heat absorption by water vaporization causes a decrease in the local adiabatic flame temperature, and this reduces the chemical reaction in the gas phase to produce thermal NO x [5]. Moreover, this helps to reduce the formation of soot, PM, CO, and HC [6,7]. Water-fuel emulsion is a system that includes water (dispersed phase with 0.1 . . . 10 µm diameter of droplets) and diesel fuel as a dispersion medium [8][9][10], and the volumetric water content in the fuel-water emulsions can vary from 5% to 30%, and is an object of strict control during the process of emulsification [4,[11][12][13].
Methods of moisture measurement can be divided into direct-based on dividing the substance under research into free water and dry material-and indirect, where moisture content can be defined based on the values of other measurements functionally connected with moisture [14][15][16]. To identify popular principles of moisture measurement, a detailed analysis of the world moisture meter market, produced in a recent period, was carried out for bulk substances and oil products as objects of measurement. In total, 358 moisture meters for bulk materials and 64 moisture meters for oil products were detected. The dielcometer (capacitive) principle of moisture measurement remains forward-looking among all indirect methods [17].
A large number of factors influence moisture measurement accuracy. Among them, we traditionally emphasize the physicochemical composition of the combustible liquid (oil fuel as a continuous phase of the emulsion) [16,18,19]. It is not only the type of fuel that has an effect on its physicochemical composition, but the conditions of extraction and processing, etc., and obtaining a full list of factors for analytic forecasting is usually complicated.
A list of analytical methods has been developed to determine the moisture content of oils, based on the fact that water has appreciably different bulk physical characteristics than the oil matrix, e.g., density, dielectric permittivity, or refractive index. These methods are usually only suitable for the analysis of substances in which the composition of the oil matrix does not change significantly, but where the water-to-oil ratio of the matrix changes. For example, the water content of an oil-in-water emulsion can be determined by measuring the density or dielectric permittivity, because the density and dielectric permittivity of water are significantly higher than those of oil. If the composition of the oil matrix changes as well as the water content, then it may not be possible to accurately determine the moisture content of the oil, because more than one oil composition may give the same value for the physical property being measured. Unfortunately, it is impossible to obtain an analytic forecast of the chemical composition and different features of all of the emulsified fuels under research. This is why most indirect methods of moisture measurement have local effectiveness, and are not versatile.
At present, approximately 33% of moisture meters are represented by dielcometer measuring instruments with capacitive sensors. Moisture meters of this type have method error (type uncertainty), the value of which directly depends on the dielectric permittivity of the dry oil as a continuous phase of the emulsion under research (oil type, for example) [16,18]. The influence of type uncertainty on the result of moisture measurement is usually significant, and traditional methods of compensation can be effective when we know the composition of a substance. Such compensation, in general, can be fulfilled using a secondary measuring transducer by using complementary reference capacitors, dedicated analytic calculations, reference calibration curves stored in the secondary measuring transducers' memory, etc. It can all be used for the limited number of materials under research, and does not allow the moisture meter to be versatile for a wide range of fuel substances. In other cases, the effectiveness of traditional compensation methods essentially decreases, because of which further improvement of the existing methods of moisture measurement with a view to solving the problem of type-uncertainty compensation is a relevant and perspective mission [20]. In-depth analysis of modern moisture control approaches with the purpose of further development and improvement is still a relevant task, and the process of searching for new methods of water-cut sensing optimization is under constant continuation [21][22][23].
"Successful" modifications of capacitive moisture meters with capacitive sensors appeared years ago [24][25][26]. An idea to obtain capacitive measurements using two or more positions for one of the capacitor plates, with further direct comparison, seemed to be of special interest [27,28]. This approach helps to eliminate parasitic capacitances almost completely, compensate for leakage currents, and reduce the influence of fringe electric fields. The main task of this research is in achieving robust, near-linear transfer function for the comparison of methods of moisture measurement, as described in [17], and to provide effective type-uncertainty compensation during the process of water-fuel emulsion moisture control.

Materials and Methods
In capacitive sensors, where one of the capacitor plates should be placed in two or more positions [27], we first take an empty sensor, where a change in electric capacitance ∆C 0 should be defined before and after the distance between two plates is changed. Then, we define the change in electric capacitance ∆C 1 when a capacitive sensor is filled with a tested substance for the same positions of the capacitor plates. After that, the relation ∆C 1 /∆C 0 should be calculated. In this method, only the accuracy of the capacitor plates' positioning would influence the uncertainty of the measurements. Thus, the transfer function of the moisture meter should be described with the following formula: where K is a normalizing coefficient [17] equal to K = 28.599 when we take the values of electric capacitances, equal to 15 pF for C 1 and C 4 , and 50 pF for C 2 and C 3 when the capacitive sensors are empty. The suggested design of the capacitive instrument measuring transducer is described below ( Figure 1). As we can see, both of its sensors consist of a system of flat plates (1), where two pairs of flat plates belong to measuring capacitors C 1 and C 4 , and the rest of the flat plates create another pair of measuring capacitors C 2 and C 3 . All flat plates of equal length l are assembled inside two fluoroplastic rings (2) at an equal gap, designated as Z.

Materials and Methods
In capacitive sensors, where one of the capacitor plates should be placed in two or more positions [27], we first take an empty sensor, where a change in electric capacitance ΔC0 should be defined before and after the distance between two plates is changed. Then, we define the change in electric capacitance ΔC1 when a capacitive sensor is filled with a tested substance for the same positions of the capacitor plates. After that, the relation ΔC1/ΔC0 should be calculated. In this method, only the accuracy of the capacitor plates' positioning would influence the uncertainty of the measurements. Thus, the transfer function of the moisture meter should be described with the following formula: where K is a normalizing coefficient [17] equal to K = 28.599 when we take the values of electric capacitances, equal to 15 pF for C1 and C4, and 50 pF for C2 and C3 when the capacitive sensors are empty. The suggested design of the capacitive instrument measuring transducer is described below ( Figure 1). As we can see, both of its sensors consist of a system of flat plates (1), where two pairs of flat plates belong to measuring capacitors C1 and C4, and the rest of the flat plates create another pair of measuring capacitors C2 and C3. All flat plates of equal length l are assembled inside two fluoroplastic rings (2) at an equal gap, designated as Z. The measuring channel should be filled with a probe of liquid fuel under moisture control, and the reference channel should be filled with the same substance, but previously dehydrated. Using the method mentioned above, it is necessary to measure the values of electric capacitances C1, C2, C3, and C4, and then to calculate the differences C2 − C1 and C3 − C4. When this is done, we can calculate the relationship of these differences (C3 − C4)/(C2 − C1), and this would be an informative parameter for further processing, directly connected with moisture content.
To implement this idea in practice, moisture content measuring instruments should have a structure, as given in Figure 2. As we can see, four capacitive sensors of the measuring transducer are connected to the inputs of corresponding monostable multivibrators MMV1…MMV4, designed to convert values of electric capacitance into a duration of rectangular pulses. The microcontroller CPU controls the monostable multivibrators' work by generating rectangular pulses of a stable frequency on one of its output pins. The measuring channel should be filled with a probe of liquid fuel under moisture control, and the reference channel should be filled with the same substance, but previously dehydrated. Using the method mentioned above, it is necessary to measure the values of electric capacitances C 1 , C 2 , C 3 , and C 4 , and then to calculate the differences C 2 − C 1 and C 3 − C 4 . When this is done, we can calculate the relationship of these differences (C 3 − C 4 )/(C 2 − C 1 ), and this would be an informative parameter for further processing, directly connected with moisture content.
To implement this idea in practice, moisture content measuring instruments should have a structure, as given in Figure 2. As we can see, four capacitive sensors of the measuring transducer are connected to the inputs of corresponding monostable multivibrators MMV1 . . . MMV4, designed to convert values of electric capacitance into a duration of rectangular pulses. The microcontroller CPU controls the monostable multivibrators' work by generating rectangular pulses of a stable frequency on one of its output pins.
Each monostable multivibrator generates a sequence of rectangular pulses, the duration of which is in direct proportion with the value of electric capacitance, connected to its input (diagrams in Figure 3).  Each monostable multivibrator generates a sequence of rectangular pulses, the duration of which is in direct proportion with the value of electric capacitance, connected to its input (diagrams in Figure 3). Low-pass filters convert pulse duration into DC voltage in proportion to the pulse duration amplitude. The multiplexer MUX, controlled by the CPU, commutes the outputs of the low-pass filters LPF1…LPF4 with the input of the normalizing amplifier NA   Each monostable multivibrator generates a sequence of rectangular pulses, the duration of which is in direct proportion with the value of electric capacitance, connected to its input (diagrams in Figure 3). Low-pass filters convert pulse duration into DC voltage in proportion to the pulse duration amplitude. The multiplexer MUX, controlled by the CPU, commutes the outputs of the low-pass filters LPF1…LPF4 with the input of the normalizing amplifier NA with different types of moist substances, we took four values of dielectric permittivity for imaginary dehydrated fuels (ε n = 2.0; ε n = 2.5; ε n = 3.0; ε n = 3.5) and four values of moisture content (W = 0%; W = 10%; W = 20%; W = 30%). The values of dielectric permittivity for moist substances, necessary to calculate the values of the four sensors' capacitances, were estimated with the help of the universal Wiener equation [29][30][31][32]. The values of capacitances C 1 , C 2 , C 3 , and C 4 , and the calculated values of moisture content defined with the help of the static function (1), are given in Table 1. Graphs of received transfer functions for different moist substances, and their comparison with an ideal linear transfer function, are given in Figure 4.
As we can see, all transfer functions received with the help of Equation (1) are nonlinear, and the calculated values of moisture content are significantly different from nominal in all points, except for W = 0%. Hence, at first, Equation (1) should be linearized.
To achieve near-linear dependence between moisture content W and the (C 3 − C 4 )/ (C 2 − C 1 ) relationship, the method of least squares and instruments of general linear regression were used. The modified transfer function for the water-cut meter is given in Equation (2), where a and b are coefficients preliminarily defined in [33] (a = 0.9622, b = 0.0371): (2)  To achieve near-linear dependence between moisture content W and the (C3 − C4)/(C2 − C1) relationship, the method of least squares and instruments of general linear regression were used. The modified transfer function for the water-cut meter is given in Equation (2) 1.019 1.269 0.00527 0.000112 .
New values of moisture content, calculated using Formula (2), are given in Table 2. The graphs of received transfer functions for different moist substances can be seen in Figure 5. If we compare Figures 4 and 5, it would be possible to say that the modified transfer function (2) is far more effective than the initial transfer function (1). New values of moisture content, calculated using Formula (2), are given in Table 2. The graphs of received transfer functions for different moist substances can be seen in Figure 5. If we compare Figures 4 and 5, it would be possible to say that the modified transfer function (2) is far more effective than the initial transfer function (1).
To prove this analytically, it seemed to be rational to calculate one of the goodness-of-fit parameters, e.g., root mean estimator (Equation (3)): For the data in Table 1, we can write: For the data in Table 2, the result would be S = 0.4158%. To prove this analytically, it seemed to be rational to calculate one of the goodnessof-fit parameters, e.g., root mean estimator (Equation (3)): For the data in Table 1, we can write: For the data in Table 2, the result would be S = 0.4158%. To carry out experimental research with the purpose of checking the workability of a transfer function (2), and to estimate the provided level of uncertainty in the moisture content measurements, it was necessary to obtain stable reference samples of oil products with different dielectric permittivity values εn in a dehydrated state, and values of moisture content equal to W = 0%, W = 10%, W = 20%, and W = 30%. We used diesel fuel with εn = 2.01 and mazut with εn = 2.67 for this purpose. Stable reference samples of diesel fuel and mazut with 10 μm diameter of water particles were obtained by adding appropriate volumes of water, with the help of a second-class precision pipette [34] and 30 min of mechanical mixing at n = 3200 rpm [35,36].
Sensors of a measuring transducer were assembled from four identical fluoroplastic rings with slots to insert flat stainless steel plates, which were glued to the internal sur- To carry out experimental research with the purpose of checking the workability of a transfer function (2), and to estimate the provided level of uncertainty in the moisture content measurements, it was necessary to obtain stable reference samples of oil products with different dielectric permittivity values ε n in a dehydrated state, and values of moisture content equal to W = 0%, W = 10%, W = 20%, and W = 30%. We used diesel fuel with ε n = 2.01 and mazut with ε n = 2.67 for this purpose. Stable reference samples of diesel fuel and mazut with 10 µm diameter of water particles were obtained by adding appropriate volumes of water, with the help of a second-class precision pipette [34] and 30 min of mechanical mixing at n = 3200 rpm [35,36].
Sensors of a measuring transducer were assembled from four identical fluoroplastic rings with slots to insert flat stainless steel plates, which were glued to the internal surface of the appropriate pair of fluoroplastic rings with the help of superglue gel, and soldered with wires to create four capacitive sensors ( Figure 6).
The process of moisture measurement was performed following the substitution method, where the value of the quantity being measured is not found directly from a reading of the measuring instrument, but rather from the magnitude of the standard, which is regulated in such a way that the reading of the measuring instrument remains the same when the quantity being measured is replaced by the standard. Commonly known, this method eliminates systematic errors and provides high accuracy, and the substitution method is extensively used in measuring electrical quantities, such as resistance, capacitance, and inductance. The connection scheme of the devices involved in the process of moisture measurement is given in Figure 7.
The setup consists of four capacitive sensors, designated as C 1 , C 2 , C 3 , and C 4 , one channeled capacitance into a DC-voltage transducer that performs a sequence of transformations "capacitance-pulse duration-dc voltage", a variable air capacitor (reference capacitor), an oscilloscope, a digital voltmeter, and an RLC meter. face of the appropriate pair of fluoroplastic rings with the help of superglue gel, and soldered with wires to create four capacitive sensors ( Figure 6). The process of moisture measurement was performed following the substitution method, where the value of the quantity being measured is not found directly from a reading of the measuring instrument, but rather from the magnitude of the standard, which is regulated in such a way that the reading of the measuring instrument remains the same when the quantity being measured is replaced by the standard. Commonly known, this method eliminates systematic errors and provides high accuracy, and the substitution method is extensively used in measuring electrical quantities, such as resistance, capacitance, and inductance. The connection scheme of the devices involved in the process of moisture measurement is given in Figure 7. The setup consists of four capacitive sensors, designated as C1, C2, C3, and C4, one channeled capacitance into a DC-voltage transducer that performs a sequence of transformations "capacitance-pulse duration-dc voltage", a variable air capacitor (reference capacitor), an oscilloscope, a digital voltmeter, and an RLC meter.
From the very beginning, a part of the measuring transducer with two capacitive sensors C1 and C2 (reference channel) was filled with the original (with traces of water) sample of oil product (Figure 8a). Both capacitive sensors of the reference channel were connected one by one to the input of capacitance into the DC-voltage transducer, and appropriate DC voltage values were taken from the screen of a digital voltmeter and fixed by the operator. Then, the variable air capacitor (reference capacitor or standard) was connected instead of the capacitive sensors. Its capacitance was slowly changed in correspondence with the method of substitution until the moment when the DC voltage value on the voltmeter's screen became equal to the values detected previously ( Figure  8b).  The process of moisture measurement was performed following the substitution method, where the value of the quantity being measured is not found directly from a reading of the measuring instrument, but rather from the magnitude of the standard, which is regulated in such a way that the reading of the measuring instrument remains the same when the quantity being measured is replaced by the standard. Commonly known, this method eliminates systematic errors and provides high accuracy, and the substitution method is extensively used in measuring electrical quantities, such as resistance, capacitance, and inductance. The connection scheme of the devices involved in the process of moisture measurement is given in Figure 7. The setup consists of four capacitive sensors, designated as C1, C2, C3, and C4, one channeled capacitance into a DC-voltage transducer that performs a sequence of transformations "capacitance-pulse duration-dc voltage", a variable air capacitor (reference capacitor), an oscilloscope, a digital voltmeter, and an RLC meter.
From the very beginning, a part of the measuring transducer with two capacitive sensors C1 and C2 (reference channel) was filled with the original (with traces of water) sample of oil product (Figure 8a). Both capacitive sensors of the reference channel were connected one by one to the input of capacitance into the DC-voltage transducer, and appropriate DC voltage values were taken from the screen of a digital voltmeter and fixed by the operator. Then, the variable air capacitor (reference capacitor or standard) was connected instead of the capacitive sensors. Its capacitance was slowly changed in correspondence with the method of substitution until the moment when the DC voltage value on the voltmeter's screen became equal to the values detected previously ( Figure  8b). From the very beginning, a part of the measuring transducer with two capacitive sensors C 1 and C 2 (reference channel) was filled with the original (with traces of water) sample of oil product (Figure 8a). Both capacitive sensors of the reference channel were connected one by one to the input of capacitance into the DC-voltage transducer, and appropriate DC voltage values were taken from the screen of a digital voltmeter and fixed by the operator. Then, the variable air capacitor (reference capacitor or standard) was connected instead of the capacitive sensors. Its capacitance was slowly changed in correspondence with the method of substitution until the moment when the DC voltage value on the voltmeter's screen became equal to the values detected previously (Figure 8b). After that, the standard capacitor was disconnected from the transducer, and its capacitance was measured with the help of an accurate RLC meter at a 10 kHz frequency. The same process took place for the capacitive sensors in the measuring channel, filled with moist samples of oil products (Figure 8d). As a result, 10 measurements of electric capacitance were received from each capacitive sensor. The oscilloscope helped to control the correct operation of the capacitance into the pulse duration transducer. At first, this helped to choose the correct operating mode of the secondary transducer, when the duration of rectangular pulses was in direct proportion with the C1, C2, C3, and C4 capacitance values in the whole range of their variation. Furthermore, the presence of the oscilloscope was necessary to detect possible oscillation stops when measuring high moisture contents (i.e., over 20%). After that, the standard capacitor was disconnected from the transducer, and its capacitance was measured with the help of an accurate RLC meter at a 10 kHz frequency. The same process took place for the capacitive sensors in the measuring channel, filled with moist samples of oil products (Figure 8d). As a result, 10 measurements of electric capacitance were received from each capacitive sensor. The oscilloscope helped to control the correct operation of the capacitance into the pulse duration transducer. At first, this helped to choose the correct operating mode of the secondary transducer, when the duration of rectangular pulses was in direct proportion with the C 1 , C 2 , C 3 , and C 4 capacitance values in the whole range of their variation. Furthermore, the presence of the oscilloscope was necessary to detect possible oscillation stops when measuring high moisture contents (i.e., over 20%).

Theory/Calculation
The next step was to check the ability of the new transfer function (2) to retain stability when working with experimental results of measurements with natural random variation. The values of electric capacitances taken from capacitive sensors were not the ideal deterministic values used in Table 1, but would always have some random variation. The presence of random variation (i.e., capacitance measurement uncertainty) influences the uncertainty of the calculated moisture content values; it can decrease the robustness of the transfer function (2) up to the moment when it becomes irrelevant [37][38][39]. To check the robustness of the transfer function (2) theoretically, we introduced symmetric variation into the values of electric capacitance. The range of variation was set at a level of 0.01% from the values of electric capacitance after the analysis of metrological performance for modern capacitance and RLC meters. The calculated values of electric capacitance for capacitive sensors placed in the reference channel are shown in Table 3, while those for the sensors of the measuring channel are shown in Table 4 (initial values of electric capacitance calculated for different substances are marked with bold print). Table 3. Capacitance values for the sensors in the reference channel, with 0.01% variation. W ε n = 2.0 ε n = 2.5 ε n = 3.0 ε n = 3.5 0%   The robustness of the transfer function (2) was checked for the conditions when each value of sensors' capacitance was the maximum or minimum within the symmetric variation range ±0.01%. The signs "+" and "−" in later designations represent the maximal or minimal values of electric capacitance. As we had four sensors, 16 different combinations were possible (Table 5). For example, the designation "+-++" means that the electric capacitances C 1 , C 3 , and C 4 take maximal values, while C 2 is minimal [40]. As we can see from Table 5, it would be nearly impossible to use separate random results of capacitance measurements together with Formula (2), because the relative fullscale error of moisture content would be rather large. Here, it would be rational to use such robust estimators as the mean (average), median, and truncated mean-in other words, to increase the accuracy of the moisture measurements, we would have to use some kind of estimator, the uncertainty of which would be many times smaller compared with the uncertainty of separate measurements.
Measured values of capacitances for the sensors C 1 , C 2 , C 3 , and C 4 are given in Table 6. Calculated values of moisture content with mean and median values of the sensors' capacitances for the transfer functions (1) and (2) can be found in Table 7.
To compare the robustness of the results in Table 7, it would be rational to use dispersion analysis. At first, four values of mean square error (dispersions of repeatability) should be calculated (4): Formula (4) requires the mean values of moisture content to be calculated for the results in Table 7. Appropriate calculations are given in Table 8.
After the application of Formula (4), we receive:  (1) with mean values of C 1 , C 2 , C 3 , and C 4 ; S 2 rep2 is a dispersion of repeatability for the transfer function (1) with median values; S 2 rep3 is a dispersion of repeatability for the transfer function (2) with mean values; and S 2 rep4 is a dispersion of repeatability for the transfer function (2) with median values.
Then four dispersions of adequacy (for the data in Table 8) should be calculated as follows (5): where N is the total number of moisture content values W calculated ; and m is the number of nominal points of moisture content W nom . is a dispersion of adequacy for the transfer function (2) with mean values; and S 2 ad4 is a dispersion of adequacy for the transfer function (2) with median values.
To compare the dispersions of adequacy with the dispersions of repeatability, we used the following F-test: F = S 2 ad /S 2 rep . All calculations for the results given in Tables 7 and 8 can be found in Table 9. Table 6. Measured values of capacitance for the sensors C 1 , C 2 , C 3 , and C 4 with diesel fuel and mazut.

Results/Discussion
As mentioned in the "Materials and Methods" section, the first attempt to provide effective "type uncertainty" compensation for capacitive moisture meters was described in [8], and transfer function (1) is the result.
The calculated values of the sensors' capacitances and moisture content for the transfer function (1) were received with a help of the universal Wiener equation as one of the most popular dielectric mixture models (6): where ε mix is the dielectric permittivity of a binary material-water mixture; ε n is the dielectric permittivity of dry material; ε W = 80, and is the dielectric permittivity of water; and W is the volumetric moisture content, taken in absolute values. If we apply the universal Wiener equation to the values of dielectric permittivity given in Table 1, and use it as the water-cut meter's transfer function, we obtain a family of curves ( Figure 9). Figure 9 illustrates the situation when the oil matrix composition changes as well as the water content (oil type uncertainty). If no type-uncertainty compensation happens, we have a family of transfer functions for the water-cut meter instead of one, taken as a nomi-nal.
The transfer function (1) helps to solve this problem, and provides the possibility of measurements in emulsified fuels, where the dielectric permittivity of a dry oil ε n varies in a wide range from 2 to 3.5, without any preliminary calibration. However, Figure 1 illustrates that a transfer function received with the help of Equation (1) is still nonlinear, the influence of fuel type uncertainty remains significant, and the calculated values of moisture content do not correspond to the nominal linear model.
To determine the linear dependence between moisture content W and the (C 3 − C 4 )/(C 2 − C 1 ) relationship, we used the method of least squares: where a and b are unknown coefficients of a first-order polynomial (7). To build a system of conditional equations, we had to choose average values of the (C 3 − C 4 )/(C 2 − C 1 ) relationship (Table 10).
A system of linear conditional equations, based on the first-order polynomial (7), is given below.   Calculated values of moisture W calc , received from a first-order polynomial (7): are given in Table 11.  Figure 9. Influence of the changes in the composition of a dry substance εn on the value of the dielectric permittivity of a binary mixture εmix.
However, Figure 1 illustrates that a transfer function received with the help of Equation (1) is still nonlinear, the influence of fuel type uncertainty remains significant, and the calculated values of moisture content do not correspond to the nominal linear model.
To determine the linear dependence between moisture content W and the (C3 − C4)/(C2 − C1) relationship, we used the method of least squares: where a and b are unknown coefficients of a first-order polynomial (7). To build a system of conditional equations, we had to choose average values of the (C3 − C4)/(C2 − C1) relationship (Table 10). A system of linear conditional equations, based on the first-order polynomial (7), is given below.  After obtaining the results of Table 11, it was possible to conclude that the least squares method is not effective for our case, as it was expected to obtain the linear transfer function of the water-cut meter ( Figure 10). The next step of the linearization was to calculate the discrepancies ∆W between nominal points of moisture content and values from Table 11, and to approximate them by applying the instruments of general linear regression. Discrepancies between nominal points of moisture content and values from Table 11, taken for ε n = 3.0, are given in Table 12. To approximate the values of discrepancies ∆W from Table 11, a sum of four functions (1, W, W 2 , and W 3 ) was taken with appropriate coefficients, defined with the help of Mathcad software (linfit function (x, y, Y)): where W is a volumetric moisture content taken from the first-order polynomial (7). Formula (8) provides an ideal approximation of the discrepancies ∆W, as can be seen in Table 13.
squares method is not effective for our case, as it was expected to obtain the linear transfer function of the water-cut meter ( Figure 10). The next step of the linearization was to calculate the discrepancies ΔW between nominal points of moisture content and values from Table 11, and to approximate them by applying the instruments of general linear regression. Discrepancies between nominal points of moisture content and values from Table 11, taken for εn = 3.0, are given in Table 12. Figure 10. Graphs of a first-order polynomial (7), placed together with an ideal transfer function of the water-cut meter. Taking into account Equations (7) and (8), we received a modified transfer function of a water-cut meter (function 2), as displayed in Figure 5.
During the comparison of the water-cut meter's transfer functions (1) and (2), it was necessary to estimate how both of them compensated for the type uncertainty, and how close both of them were to the nominal linear transfer function. Dispersions of repeatability, calculated for transfer function (1) and transfer function (2) with mean and median values, describe a possible variation in the measured moisture values for two chosen substances (diesel fuel and mazut). We can see that the values of S 2 rep1 and S 2 rep2 are approximately 1.35 times greater in comparison with S 2 rep3 and S 2 rep4 . This means that transfer function (2) provides a lower level of variability of the measured values, depending on the type of oil fuel under research. In other words, transfer function (2) is not as sensitive to fuel type as transfer function (1).
Dispersions of adequacy reflect the correspondence of experimental transfer functions (1) and (2) (taken for diesel fuel and mazut as dielectric substances with a significant difference in relative permittivity) with the nominal linear transfer function of a water-cut meter (the ideal transfer function given in Figures 4 and 5). Again, it is possible to see that dispersions of adequacy calculated for the transfer function (2) give us values that are approximately 4.5 times smaller in comparison with the dispersions of transfer function (1). This means that Equation (2) provides 4.5 times better adequacy to the ideal than Equation (1).
As a general result, we can conclude that it would be rational to choose transfer function (2) with median values of moisture content as a nominal transfer function of the water-cut meter, because it has the smallest F-test value (F = 6.19).
Moisture measurement uncertainty can be estimated by using the method of mathematical programing [41]. If the function f (x) is continuous, it is possible to find maximal and minimal values of f (x) within the limit value of the function's argument error.

Conclusions
At first, the workability of a transfer function suggested in [17] with different types of oil fuels was checked. Four values of dielectric permittivity for dehydrated substances (ε n = 2.0; ε n = 2.5; ε n = 3.0; ε n = 3.5) and four values of moisture content (W = 0%; W = 10%; W = 20%; W = 30%) were taken for that purpose. The obtained graphs of transfer functions for the mentioned substances happened to be nonlinear, and the calculated values of moisture content were significantly different from the nominal at all points, except for W = 0%. To achieve near-linear transfer function of the capacitive water-cut meter, the method of least squares was used. The next step of linearization was in calculating discrepancies between the nominal points of moisture content and values received after the application of the LS method, with their approximation using general linear regression instruments.
Analyzing the graphs of the modified transfer functions for different moist substances, it was possible to say that they happened to be far more effective than the initial transfer function suggested in [17]. The root mean estimator was calculated for the initial transfer function, and the static function-received after general linear regression-as an integral difference between the nominal and calculated values of moisture content. The corresponding root mean estimator values were 1.3062% and 0.4158%, proving the effectiveness of transfer function (2), modified with the help of general linear regression instruments.
To check the workability of the modified transfer function in real conditions, two different fuels were taken: diesel fuel and mazut. The values of moisture content for two transfer functions, calculated for the mean and median values of C 1 , C 2 , C 3 , and C 4 , were compared with the help of dispersion analysis. After the dispersions of repeatability and adequacy had been calculated and compared, it was possible to conclude that the modified transfer function of a water-cut meter, calculated for median values of C 1 , C 2 , C 3 , and C 4 , provided minimal dispersions of repeatability and adequacy, and a minimal F-test value. Values of both dispersions for the modified transfer function were significantly smaller, as were the F-test values, proving its more robust properties (i.e., capacity for type-uncertainty compensation) and better adequacy to the nominal linear transfer function of the water-cut meter.