Evaluation of the Influence of Water Content in Oil on the Metrological Performance of Oil Flow Measurement Systems

: According to current Brazilian regulations, if the volumes of oil produced used as a reference for the payment of government shares and third parties contain a water content greater than 2% v / v , these volumes must be arbitrarily increased between 1.44% and 10.89% due exclusively to water content, which has caused operational problems for oil companies such as differences between volumes produced and volumes sold, and additional payments from government shares and third parties. This study aimed to evaluate the metrological performance of oil measurement systems with ultrasonic, Coriolis and positive displacement flow meters when subjected to varied water content, fluid temperature and flow rate conditions using the Design of Experiments and the Response Surface Methodology. The analysis of variance showed that the models presented good fits for the ultrasonic meter (coefficient of determination R 2 of 97.96%, p -value of 0.001, and a standard deviation of 5.89 × 10 − 5 ); Coriolis meter (R 2 of 90.91%, p -value of 0.037, and a standard deviation of 5.88 × 10 − 5 ); and positive displacement meter (R 2 of 99.07%, p -value of 0.000, and a standard deviation of 4.85 × 10 − 5 ). The results of the experiments carried out indicate that the contribution of each parameter analyzed to the metrological performance of the measurement system varies depending on the measurement technology used by the flow meter. However, the fluid temperature proved to be a relevant parameter common to all flow measurement technologies evaluated. All measuring technologies evaluated were influenced by water content in the range of 0% to 10% v / v , with the measurement error being less than 0.2% when compared to a standard positive displacement type meter in almost all experimental conditions. The Coriolis-type flow meter was the one that presented the smallest error among the measuring technologies evaluated.


Introduction
Globally, the oil and gas industry presents several challenges, especially in multiphase flow, where measurement uncertainties are due to different flow rates and BSW [1]; on the other hand, different operating conditions (flow and temperature) can be correlated with a maximum BSW content to ensure reliable measurements [2].Furthermore, in offshore oil production, fluctuations in flow rates and BSW content can destabilize the process, resulting in operational losses and environmental fines [3].
Government participation generated by the exploration of oil and natural gas in Brazil is an important source of income for the development of Brazilian society, with these resources being allocated especially to education and health.In 2021, the amount collected by the Union in the form of government participation was more than USD 15.5 billion [4], which makes the quality of the measurement systems used to control the production of oil and natural gas and the existence of a robust measurement management system to guarantee the reliability of the results obtained.
In Brazil, the measurement systems used to control the production of oil and natural gas in an oil field are governed by the Technical Regulation for the Measurement of Oil and Natural Gas (RTM), approved by Joint Resolution ANP (Brazilian National Agency of Petroleum, Natural Gas and Biofuels) and INMETRO (Brazilian Institute for Metrology, Quality and Technology) 1/2013, which establishes the minimum technical, constructive and metrological conditions and requirements that these measuring systems must comply with, with a view to guaranteeing the credibility of the measurement results [5].
Online petroleum fiscal measurement systems must consist of at least the following equipment: a flow meter; an instrument for measuring temperature; an instrument for measuring static pressure; and a flow computer for correction of volumes based on flow, temperature and static pressure information.
As determined by this Brazilian regulation, fiscal petroleum measurement systems must be located immediately downstream of petroleum treatment systems to ensure that the measured petroleum is stabilized and does not contain a water content (BSW) greater than 1% v/v.Additionally, for events where the average BSW is greater than 2% v/v, the volumes of oil produced must be increased by a correction factor pre-defined by the ANP that can vary between 1.44% and 10.89% depending on exclusively from the average daily BSW [6].
The increase in the volumes of oil produced because of oil non-compliance events due to arbitrary factors can cause operational problems for oil companies.A one-day oil decommissioning event can result in an additional payment of USD 144,000 in government participation in Brazil.
The increase in produced volumes by arbitrary correction factors based exclusively on average BSW may not be sufficient to characterize the metrological performance of online fiscal petroleum measurement systems operating in high BSW conditions.Hallanger et al. [7] and Skea and Hall [8] indicate that constructive aspects of the measurement system, such as flow meter technology and characteristics of the measured oil, such as temperature, kinematic viscosity and density, can have an influence on the metrological performance of measurement systems oil in line operating under high BSW conditions.
Water content in crude oil is the subject of several studies in the most diverse areas, such as high-performance biodegradable [9], the development of new microwave sensors [10,11] and quality control [12].
Economically, it is important for the oil industry that the amount paid as government shares and to third parties reflects the volumes of oil produced and that there is consistency between these volumes and the volumes of oil sold.Thus, a gap was identified in the literature regarding the performance of a sensitivity study for the simultaneous assessment of factors influencing the metrological performance of online oil measurement systems.
This study aimed to evaluate, by means of the Design of Experiments and Response Surface Methodology, the influence of BSW, fluid temperature and flow rate on the metrological performance of oil measurement systems with ultrasonic, Coriolis and positive displacement flow meters.

Rationale
The Constitution of the Federative Republic of Brazil of 1988 determines that mineral resources, which include oil and natural gas deposits, constitute assets of the Union.In this sense, in the case of the exploitation of these non-renewable resources by authorized companies, financial compensation to the Union is provided for, which varies according to the exploration and production contract signed and which, generically, is called Government Participation (PGov).
Based on Brazilian laws, the following government and third-party participations are considered: signature bonus, royalties, special participation, payment for occupation or retention of area, payment to landowners, and surplus oil.

Technical Regulations for the Oil and Natural Gas Area
The Brazilian RTM established the minimum conditions and requirements for oil and natural gas measurement systems with a view to ensuring accurate and complete results in line with the criticality of the process of determining oil and natural gas production.
The Technical Regulation for Notification of Failures of Oil and Natural Gas Measurement Systems and Oil Framing Failures (RTNF), approved by ANP Resolution No. 18 of 27 March 2014, defined deadlines and procedures to be adopted with relation to the occurrence of failure events associated with the measurement systems provided for in the RTM [13].
To guarantee the Brazilian interests, this regulation established standardized methodologies for estimating the volumes affected in the event of a measurement failure event.
For oil classification failure events, i.e., where volumes produced with average daily BSW greater than 1% v/v are measured, aiming to mitigate a potential increase in measurement uncertainty in a fiscal oil measurement system, it is mandatory to apply correction factors to the volume produced.Almeida et al. [14] concluded that BSW uncertainty and flow meter calibration factor uncertainty are the main sources of uncertainty for in-line petroleum measurement systems.
Table 1 presents the correction factors for the volume of oil produced according to the average daily BSW predicted in Circular Letter No. 2/2021/NFP/ANP-e [6].
The Metrological Technical Regulation for Dynamic Measurement Systems equipped with meters for quantities of petroleum and its liquid derivatives, approved by INMETRO Ordinance No. 291 of 7 July 2021, established the metrological and technical requirements, such as accuracy class and maximum error permissible, which flow measurement systems must comply with for their use in determining the volumes of oil provided for in the RTM [15].
Fines are the main penalty applied to companies that produce oil and natural gas as part of the inspection process of measurement systems conducted by the ANP's Production Inspection Center (NFP).Fines applicable to infractions associated with oil and natural gas measurement systems range from 5 thousand to 2 million reais.The value of the fine is defined by the dosimetry provided for in Decision No. 2/2021/NFP [6] and depends, among other factors, on the severity of the infraction, the offender's background, and the volume of oil and natural gas moved in the installation as detailed in Equation (1).
where VP is the value of the penalty in reais; V min is the minimum value of the penalty in reais according to Law no.9.847/99 in which the infraction was framed; V max is the maximum value of the penalty in reais according to Law no.9.847/99 in which the infraction was framed; FGI is the severity factor of the offense committed; FM is the production movement factor for the installation in which the infraction was committed; FCEI is the factor of the economic condition of the company responsible for operating the field where the infraction occurred; FAI is the background factor of the company responsible for operating the field where the infraction occurred; FVA is the economic advantage factor obtained by the company responsible for operating the field where the infraction occurred; A conf is the mitigating factor due to spontaneous communication or confession of the infraction to the ANP; A rec is the mitigating factor in the case of recalculation of the volumes affected by the occurrence of the infraction.API MPMS [16] defines requirements for the application of liquid flow measurement systems.According to this standard, the volumetric flow of oil in the standard measurement condition, that is, at a temperature of 20 • C and a static pressure of 101.325 kPa, is given by Equation (2).
where Q b is volumetric oil flow in standard measurement conditions in m 3 /h, Q is the gross volumetric flow of oil in operating conditions in m 3 /h, MF is the meter factor, CTL is the conversion factor from the volume measured at the operating temperature to the volume at the standard temperature condition, CPL is the conversion factor from the volume measured at operating pressure to the volume at standard pressure condition and BSW is the water content in the oil in the operating condition in % v/v.

General Design of Experiments and Response Surface Methodology
An experiment or trial can be defined as a set of actions to analyze a process or system with the aim of obtaining new information, evaluating hypotheses, or reaffirming previously obtained results.The basic concept of an experiment is, based on the selection of control factors (process input), observing the response variable (process output) to analyze the relationship between them.
Depending on the process analyzed, an experiment can be complex and costly due to the high number of possible combinations between control factors and response variables.In this context, the use of a statistical tool such as the Design of Experiments (DoE) becomes critical for obtaining the most information at the lowest cost.Montgomery [17] cites the advantages of using the DoE tool, such as reducing cost and product development time and improving process performance.The DoE tool has been widely used in various branches of science and industry: Guerra et al. [18] used the Central Composite Design to optimize acceptance criteria for instrument calibration in oil and natural gas measurement systems; Sengel-Turk et al. [19] applied the DoE tool to evaluate new techniques for cancer treatment; Sibiya et al. [20] used DoE concepts to optimize water treatment; a central composite design was used for optimizing turbidity removal in a reactor system [21]; a Box-Behnken design was used to analyze the effects of pollutant interactions in pollutants oilfield [22]; a Box-Behnken design was used to identify optimum conditions in drilling and completion operations [23]; and a 3 2 factorial design was used for assessing the relevance of changing in test method for kinematic viscosity of opaque oils [24].
Montgomery [17] proposed a basic protocol for using the DoE approach as follows: (i) definition of the objectives of the experiment; (ii) parameter mapping; (iii) selection of control factors and response variables; (iv) definition of the experimental matrix; (v) carrying out the experiment; (vi) data analysis; and (vii) interpretation of results.

Box-Behnken Design
The Box-Behnken Design (BBD), which is connected to Response Surface Methodology (RSM), is a very powerful type of experimental design when considering a model with second-order adjustment.To use the BBD, there must be at least three control factors, and each factor must always be observed at three levels (minimum, medium and maximum).The experiment must be conducted so that each of the levels is positioned equidistant from the center of the model.In a BBD with three control factors, this is achieved by always fixing one of the levels at the midpoint and obtaining the response variables for each combination of the other factors at the minimum and maximum levels.Additionally, in BBD, observations are made at the central point, with all control factors at the average level, to determine the experimental error.
Manohar et al. [25] described the BBD as spherical since all points of the model are inscribed on a sphere of radius √ 2 and no points are located at the vertices of the area delimited by the minimum and maximum levels of each control factor.
Table 2 presents the generic experimental matrix for a BBD with three control factors and three central points.The number of treatments for a BBD is defined by Equation (3) [25]: where N is the number of treatments for BBD; k is the number of control factors in the experiment; and C 0 is the number of observations of the response variable at the center point.One advantage of BBD is that it requires fewer observations than other types of DoE, such as full factorial design.Additionally, in a BBD, extreme conditions-that is, when all control factors are at their minimum or maximum levels-are not necessary.On the other hand, this design is not suitable for experiments where the observation of control factors at their extreme levels may be relevant.

Experimental Procedure
To evaluate the metrological performance of oil flow measurement systems, a test bench was built especially for the experiment on the Technological Research Institute (IPT) premises in São Paulo, Brazil.This test bench consisted of a closed circuit with all flow meters to be evaluated arranged in series in a horizontal line.
Table 3 presents the oil flow measurement technologies evaluated and the details of each meter model used in the experiment.All measurement systems were calibrated by accredited laboratories [26].The schematic diagram of the closed circuit used as a test bench is shown in Figure 1.

Test Bench Validation
The Technological Research Institute has a bench accredited by ISO/IEC 1702 for calibrating oil flow meters, which uses a compact tester-type reference standa calibrations in the range of 0.8 to 800 m 3 /h with mineral oil with a kinematic vis between 7 and 350 cSt under standard measurement conditions (20 °C and 101.325 k The bench validation was carried out by comparing the results obtained for t trasonic and Coriolis type meters under evaluation and the standard positive disp ment type meter with the results obtained on the accredited bench.
For comparison between benches, the normalized error En according to Equatio By this criterion, if | | ≤ 1, the results can be considered compatible.

𝐸 = 𝑋 − 𝑋 𝑈 + 𝑈
where Xlab is the experiment's result, Xref is the reference value, Ulab is the expanded u tainty associated with the experiment's outcome, and Uref is the expanded uncertain sociated with the reference value.Tables 4-6 present the comparison results for each meter.

Test Bench Validation
The Technological Research Institute has a bench accredited by ISO/IEC 17025 [26] for calibrating oil flow meters, which uses a compact tester-type reference standard for calibrations in the range of 0.8 to 800 m 3 /h with mineral oil with a kinematic viscosity between 7 and 350 cSt under standard measurement conditions (20 • C and 101.325 kPa).
The bench validation was carried out by comparing the results obtained for the ultrasonic and Coriolis type meters under evaluation and the standard positive displacement type meter with the results obtained on the accredited bench.
For comparison between benches, the normalized error E n according to Equation (4).By this criterion, if |E n | ≤ 1, the results can be considered compatible.
where X lab is the experiment's result, X ref is the reference value, U lab is the expanded uncertainty associated with the experiment's outcome, and U ref is the expanded uncertainty associated with the reference value.Tables 4-6 present the comparison results for each meter.For the Coriolis meter, it was only possible to evaluate three gradients due to operational limitations.
In Tables 4-6, Q test is the flow rate of the meter on the test bench in m 3 /h; E test is the average flow error of the meter on the test bench in %; U test is the expanded uncertainty associated with the result of the flow measurement on the test bench in %; Q IPT is the flow rate of the meter on the IPT accredited bench in m 3 /h; E IPT is the average error of the meter flow on the IPT accredited bench in %; U IPT is the expanded uncertainty associated with the result of the flow measurement on the IPT accredited bench in %; E n is the normalized error.
Based on the results presented, given that the absolute normalized error |E n | was less than unity for all experimental conditions, the absence of systematic errors, the test bench built for the experiment, was considered validated.

Experimental Matrix
Using BBD, an experimental matrix was proposed for each flow meter to be evaluated, as presented in Table 2.The control factors chosen to evaluate the metrological performance of the flow meters were the BSW, the fluid temperature and the volumetric flow rate (Table 7).The levels for the BSW were defined considering the most frequent range of values from oil classification failure events in fiscal measurement systems communicated to the ANP.They were the same for all flow meters evaluated.
The fluid temperature levels considered for the experiment were chosen based on the average temperature conditions found in online petroleum fiscal measurement systems, and, like BSW, they were the same for all oil flow meters evaluated.
Finally, the flow rate levels are different for each meter evaluated and were chosen considering the test bench's flow rate limit (200 m 3 /h) and the flow limits of each meter, as detailed in Table 3. Due to the logistical difficulties of carrying out the experiment with oil in the laboratory, mineral oil with a viscosity like that produced in the Brazilian pre-salt fields (20 cSt under standard measurement conditions, 20 • C and 101.325 kPa) was decided upon.
Initially, the closed circuit was filled with pure mineral oil, that is, with BSW equal to 0% v/v.The volume of mineral oil inserted into the closed circuit was calculated based on the fluid's density, and the fluid's mass was determined using a scale.

Fluid Change
The change in the BSW level was evaluated by inserting water into the closed circuit.Figure 2 shows the monitoring of the results of measuring the density of the fluid using the Coriolis flow meter.After the insertion of water into the circuit, until the fluid is homogenized, the density tends to vary until it reaches a region of stability where it is possible to continue with the other stages of the experiment.
The fluid temperature levels considered for the experiment were chosen based on the average temperature conditions found in online petroleum fiscal measurement systems, and, like BSW, they were the same for all oil flow meters evaluated.
Finally, the flow rate levels are different for each meter evaluated and were chosen considering the test bench's flow rate limit (200 m 3 /h) and the flow limits of each meter, as detailed in Table 3.
Due to the logistical difficulties of carrying out the experiment with oil in the laboratory, mineral oil with a viscosity like that produced in the Brazilian pre-salt fields (20 cSt under standard measurement conditions, 20 °C and 101.325 kPa) was decided upon.
Initially, the closed circuit was filled with pure mineral oil, that is, with BSW equal to 0% v/v.The volume of mineral oil inserted into the closed circuit was calculated based on the fluid's density, and the fluid's mass was determined using a scale.

Fluid Change
The change in the BSW level was evaluated by inserting water into the closed circuit.Figure 2 shows the monitoring of the results of measuring the density of the fluid using the Coriolis flow meter.After the insertion of water into the circuit, until the fluid is homogenized, the density tends to vary until it reaches a region of stability where it is possible to continue with the other stages of the experiment.

Results and Discussion
The objective of the experiment was to evaluate the impact of different measurement technologies on the metrological performance of oil flow measurement systems calibrated under different conditions of BSW, fluid temperature and flow rate (control factors).

Results and Discussion
The objective of the experiment was to evaluate the impact of different measurement technologies on the metrological performance of oil flow measurement systems calibrated under different conditions of BSW, fluid temperature and flow rate (control factors).
For this analysis, the response variable of the experiment was the absolute error between the difference in the volumetric flow rate indicated by meters with different technologies and the volumetric flow rate of the standard meter as defined in Equation ( 5): where Q meter is the corrected flow indicated by the flow meter being evaluated in m 3 /h, and Q standard is the corrected flow indicated by the standard flow meter in m 3 /h.The corrected flow rates Q meter and Q standard are expressed in the standard measurement condition provided for in the RTM, that is, at a temperature of 20 • C and a static pressure of 101.325 kPa.The gross flow rates recorded by the flow meters were corrected using the temperature correction factors (CTL) and pressure correction factors (CPL) provided for in the API MPMS Chapter 11.1 standard [16].
Here, DoE is considered a 0.05 significance level.Next, the results for each oil flow meter technology evaluated from the experimental matrix for BBD using Minitab 19 software and Microsoft Excel 365 are presented, and the analysis of variance (ANOVA) method and the surface methodology response (RSM) are applied.
The control factors and their interactions are detailed below: x 1 is the linear term of the fluid temperature; x 2 is the linear term of the flow rate; x 3 is the linear term of BSW; x 2 1 is the quadratic term of the fluid temperature; x 2 2 is the quadratic term of the flow rate; x 2 3 is the quadratic term of BSW; x 1 × x 2 is the linear interaction between fluid temperature Energies 2024, 17, 2355 9 of 19 and flow rate; x 1 × x 3 is the linear interaction between fluid temperature and BSW; and x 2 × x 3 is the linear interaction between the flow rate and the BSW.

Ultrasonic Meter
BBD was used with the control factors, and the response variable analyzed for the ultrasonic flow meter is presented in Table 8 below: The effects of control factors and interactions between these factors on the response variable using ANOVA are presented in Table 9: The Pareto chart for the estimated standardized effects is presented in Figure 3: The Pareto chart for the estimated standardized effects is presented in Figure 3: From Tables 8 and 9 and Figure 3, the influence of the fluid temperature ( ) in the response variable is statistically significant.This statement is also valid for the quadratic contribution ( ) of this factor of control.The mathematical model, Equation ( 6) that describes the relationship between control factors and the variable in response obtained was as follows:  From Tables 8 and 9 and Figure 3, the influence of the fluid temperature (x 1 ) in the response variable is statistically significant.This statement is also valid for the quadratic contribution (x 2 1 ) of this factor of control.The mathematical model, Equation ( 6) that describes the relationship between control factors and the variable in response obtained was as follows: This model presented a good fit with a coefficient of determination R 2 of 97.96% and a standard deviation of 5.89 × 10 −5 .
Using the RSM, it is possible to observe the response variable.Tables S1-S3 present the model's results mathematically detailed in Equation ( 6) considering the variation from the fluid temperature and the flow rate for the BSW equals 0% v/v, 5% v/v and 10% v/v.
Below are presented the response surfaces and contour plots of the ultrasonic flow meter error for each BSW condition evaluated (Figures 4-6).
Energies 2024, 17,2355 11 of 20 This model presented a good fit with a coefficient of determination R 2 of 97.96% and a standard deviation of 5.89 × 10 −5 .
Using the RSM, it is possible to observe the response variable.Tables S1-S3 present the model's results mathematically detailed in Equation ( 6) considering the variation from the fluid temperature and the flow rate for the BSW equals 0% v/v, 5% v/v and 10% v/v.
Below are presented the response surfaces and contour plots of the ultrasonic flow meter error for each BSW condition evaluated (Figures 4-6).The metrological performance of ultrasonic flow meters based on the transit time method, such as the flow meter evaluated in the experiment, is associated with the occurrence of a fully developed flow in a turbulent regime (Re > 4000).Since the lower the fluid temperature, the higher its viscosity and the higher the fluid viscosity, the lower the Reynolds number of the flow, it can be seen that the highest errors were observed close to the lowest levels for the fluid temperature (25 • C) and the flow rate (125 m 3 /h), which may be an indication that the meter is operating in non-ideal conditions (laminar regime).
Additionally, it is worth highlighting that these non-contact mechanical flow meters are recommended for clean and homogeneous fluids [31].The presence of a two-phase flow can lead to the attenuation of the ultrasonic wave and, consequently, an increase in measurement error.This scenario can be evidenced by the increase in measurement error from the increase in BSW, as shown in Figures 4-6 for a fully developed flow in a turbulent regime.
Table 10 presents a comparison of the error for the same flow rate condition (175 m 3 /h) considering two different fluid temperature conditions (34 • C and 43 • C).The error increases as the BSW increases, which may be associated with the occurrence of a two-phase flow and an eventual attenuation of the ultrasonic wave in this system.

Coriolis Meter
The experimental matrix for BBD used with the control factors and the response variable analyzed for the Coriolis flow meter is presented in Table 11: The effects of control factors and interactions between these factors on the response variable using ANOVA are presented in Table 12: The Pareto chart for the estimated standardized effects is presented in Figure 7: The Pareto chart for the estimated standardized effects is presented in Figure 7: From Tables 11 and 12 and Figure 7, the influence of fluid temperature ( ) and its interaction with BSW ( ×  ) on the response variable are statistically significant.The mathematical model that described the relationship between control factors and variables in response obtained was as follows: From Tables 11 and 12 and Figure 7, the influence of fluid temperature (x 1 ) and its interaction with BSW (x 1 × x 3 ) on the response variable are statistically significant.The mathematical model that described the relationship between control factors and variables in response obtained was as follows: This model presented a good fit with a coefficient of determination R 2 of 90.91% and a standard deviation of 5.88 × 10 −5 .
Using the RSM, it was possible to observe the response variable.Tables S4-S6 presented the results of the mathematical model detailed in Equation ( 7) considering the variation in fluid temperature and flow rate for the BSW equal to 0% v/v, 5% v/v and 10% v/v.
Below are presented the response surfaces and the contour plots of the error of the Coriolis meter for each of the BSW conditions evaluated (Figures 8-10).
This model presented a good fit with a coefficient of determination R 2 of 90.91% and a standard deviation of 5.88 × 10 −5 .
Using the RSM, it was possible to observe the response variable.Tables S4-S6 presented the results of the mathematical model detailed in Equation ( 7) considering the variation in fluid temperature and flow rate for the BSW equal to 0% v/v, 5% v/v and 10% v/v.
Below are presented the response surfaces and the contour plots of the error of the Coriolis meter for each of the BSW conditions evaluated (Figures 8-10).For higher BSW content, the error increases with extreme values of flow and temperature.Within the same BSW contents, the error increases for the smallest and largest values of flow and temperature.
Coriolis flow meters are mass meters; that is, they measure the mass flow rate of the flow.The volumetric flow rate of the flow is obtained by determining the density of the fluid using the flow meter itself.In this sense, this type of flow meter must be calibrated with a fluid like the one that was measured by it to minimize the error associated with the volumetric flow.
Based on Figure 2, changing the fluid's BSW level causes a change in its density.Additionally, changes in fluid temperature lead to changes in its density.Such premises confirm the mathematical model described in Equation ( 7) and the results presented in Tables 11 and 12 and Figure 7, which indicate the fluid temperature and its interaction with the BSW as the main sources of influence on the measurement error.

Positive Displacement Meter
The experimental matrix for BBD used with the control factors and the response variable analyzed for the positive displacement type flow meter is presented in Table 13:  For higher BSW content, the error increases with extreme values of flow and temperature.Within the same BSW contents, the error increases for the smallest and largest values of flow and temperature.
Coriolis flow meters are mass meters; that is, they measure the mass flow rate of the flow.The volumetric flow rate of the flow is obtained by determining the density of the fluid using the flow meter itself.In this sense, this type of flow meter must be calibrated with a fluid like the one that was measured by it to minimize the error associated with the volumetric flow.
Based on Figure 2, changing the fluid's BSW level causes a change in its density.Additionally, changes in fluid temperature lead to changes in its density.Such premises confirm the mathematical model described in Equation ( 7) and the results presented in Tables 11 and 12 and Figure 7, which indicate the fluid temperature and its interaction with the BSW as the main sources of influence on the measurement error.

Positive Displacement Meter
The experimental matrix for BBD used with the control factors and the response variable analyzed for the positive displacement type flow meter is presented in Table 13: The effects of control factors and interactions between these factors on the response variable using ANOVA are presented in Table 14: The Pareto chart for the estimated standardized effects is presented in Figure 11: Energies 2024, 17, 2355 16 of 20 The effects of control factors and interactions between these factors on the response variable using ANOVA are presented in Table 14: The Pareto chart for the estimated standardized effects is presented in Figure 11: From Tables 13 and 14 and Figure 11, the influence of fluid temperature ( ) and flow rate ( ) on the response variable is statistically significant.Additionally, the quadratic portion of the fluid temperature ( ) and the interaction between the fluid temperature and the flow rate ( ×  ) influence the second-order model estimated for the response variable.
The mathematical model that describes the relationship between the control factors From Tables 13 and 14 and Figure 11, the influence of fluid temperature (x 1 ) and flow rate (x 2 ) on the response variable is statistically significant.Additionally, the quadratic portion of the fluid temperature (x 2  2 ) and the interaction between the fluid temperature and the flow rate (x 1 × x 2 ) influence the second-order model estimated for the response variable.
Using the RSM, the response variable can be observed.Tables S7-S9 present the results of the mathematical model detailed in Equation ( 8), considering the variation in fluid temperature and flow rate for the BSW, which is equal to 0% v/v, 5% v/v and 10% v/v.
The response surfaces and the contour plots of the error of the positive displacement flow meter for each of the BSW conditions evaluated are presented in Figures 12-14.
This model presented a good fit with a coefficient of determination R 2 of 99.07% and a standard deviation of 4.58 × 10 −5 .
Using the RSM, the response variable can be observed.Tables S7-S9 present the results of the mathematical model detailed in Equation ( 8), considering the variation in fluid temperature and flow rate for the BSW, which is equal to 0% v/v, 5% v/v and 10% v/v.
The response surfaces and the contour plots of the error of the positive displacement flow meter for each of the BSW conditions evaluated are presented in Figures 12-14 This model presented a good fit with a coefficient of determination R 2 of 99.07% and a standard deviation of 4.58 × 10 −5 .
Using the RSM, the response variable can be observed.Tables S7-S9 present the results of the mathematical model detailed in Equation ( 8), considering the variation in fluid temperature and flow rate for the BSW, which is equal to 0% v/v, 5% v/v and 10% v/v.
The response surfaces and the contour plots of the error of the positive displacement flow meter for each of the BSW conditions evaluated are presented in Figures 12-14     Positive displacement flow meters that use oval gears work by discretizing the flow into portions of known volume through the action of the flow on a rotating element.For the rotating element to fulfill its role, that is, to transport the discretized volume from the flow meter inlet to its outlet, it is necessary that the rotating element is not in contact with the meter shell, which can cause a small amount of volume (slide).
The greater the slip, the greater the measurement error.Slippage is directly proportional to the viscosity of the fluid and inversely proportional to the flow rate, especially when the meter operates close to the minimum flow limit for which it was designed.
Additionally, the measurement error can be impacted by physical changes in the flow meter due to the thermal expansion of the materials used, which leads to changes in the discretized volume with each rotation of the measuring element.
This premise is validated by the results obtained, as the mathematical model detailed in Equation (8) indicates that the main influencing factors are the fluid temperature (and its quadratic portion), the flow rate and the linear interaction between the fluid temperature and the flow rate, with the largest error (0.179%) being observed for the largest fluid temperature level (43 °C) and the lowest flow rate level (35 m 3 /h).

Conclusions
This study aimed to evaluate, using DoE and RSM, the metrological performance of oil measurement systems with ultrasonic flow meters, Coriolis and positive displacement when subjected to different BSW operating conditions, fluid temperature and flow rate.
From the execution of the experiment detailed in this study, it was found that all measuring technologies evaluated suffer low influence from BSW in the range of 0% to 10% v/v, with the measurement error being less than 0.2% in almost all operational conditions evaluated, when compared to a standard positive displacement type meter in almost all experimental conditions.
Additionally, it can be proven that the influence of the control factors evaluated on the response variable varies depending on the measurement technology used by the flow meter.However, it was found that fluid temperature is a relevant parameter in evaluating the metrological performance of oil flow measurement systems for all measurement technologies evaluated.
Finally, it is concluded that the Coriolis flow meter presented the lowest error among the measuring technologies evaluated.Additionally, ultrasonic flow meters should be avoided with flows that are not completely developed or in laminar flows.
To deepen discussions on the metrological performance of oil measurement systems, it is suggested to evaluate the issues by undertaking the following: (i) carry out other experiments using petroleum instead of mineral oil; (ii) evaluate the metrological performance of petroleum measurement systems with other measurement technologies, such as Positive displacement flow meters that use oval gears work by discretizing the flow into portions of known volume through the action of the flow on a rotating element.For the rotating element to fulfill its role, that is, to transport the discretized volume from the flow meter inlet to its outlet, it is necessary that the rotating element is not in contact with the meter shell, which can cause a small amount of volume (slide).
The greater the slip, the greater the measurement error.Slippage is directly proportional to the viscosity of the fluid and inversely proportional to the flow rate, especially when the meter operates close to the minimum flow limit for which it was designed.
Additionally, the measurement error can be impacted by physical changes in the flow meter due to the thermal expansion of the materials used, which leads to changes in the discretized volume with each rotation of the measuring element.
This premise is validated by the results obtained, as the mathematical model detailed in Equation (8) indicates that the main influencing factors are the fluid temperature (and its quadratic portion), the flow rate and the linear interaction between the fluid temperature and the flow rate, with the largest error (0.179%) being observed for the largest fluid temperature level (43 • C) and the lowest flow rate level (35 m 3 /h).

Conclusions
This study aimed to evaluate, using DoE and RSM, the metrological performance of oil measurement systems with ultrasonic flow meters, Coriolis and positive displacement when subjected to different BSW operating conditions, fluid temperature and flow rate.
From the execution of the experiment detailed in this study, it was found that all measuring technologies evaluated suffer low influence from BSW in the range of 0% to 10% v/v, with the measurement error being less than 0.2% in almost all operational conditions evaluated, when compared to a standard positive displacement type meter in almost all experimental conditions.
Additionally, it can be proven that the influence of the control factors evaluated on the response variable varies depending on the measurement technology used by the flow meter.However, it was found that fluid temperature is a relevant parameter in evaluating the metrological performance of oil flow measurement systems for all measurement technologies evaluated.
Finally, it is concluded that the Coriolis flow meter presented the lowest error among the measuring technologies evaluated.Additionally, ultrasonic flow meters should be avoided with flows that are not completely developed or in laminar flows.
To deepen discussions on the metrological performance of oil measurement systems, it is suggested to evaluate the issues by undertaking the following: (i) carry out other experiments using petroleum instead of mineral oil; (ii) evaluate the metrological performance of petroleum measurement systems with other measurement technologies, such as turbine-type flow meters; (iii) evaluate the metrological performance of oil measurement systems under more severe BSW conditions; (iv) carry out other experiments with different models of the same measurement technology; and (v) evaluate the influence of the Reynolds number on the metrology of oil measurement systems with ultrasonic flow meters.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/en17102355/s1.Table S1.Experiment results for the ultrasonic flow meter under BSW conditions were equal to 0% v/v.Table S2.Experiment results for the ultrasonic flow meter under BSW conditions were equal to 5% v/v.Table S3.Experiment results for the ultrasonic flow meter under BSW conditions were equal to 10% v/v.Table S4.Experiment results for Coriolis meter under BSW condition equal to 0% v/v.Table S5.Experiment results for Coriolis meter under BSW condition equal to 5% v/v.Table S6.Experiment results for Coriolis meter under BSW condition equal to 10% v/v.Table S7.Experiment results for positive displacement type flow meter under BSW condition equal to 0% v/v.Table S8.Experiment results for positive displacement type flow meter under BSW condition equal to 5% v/v.Table S9

Funding:
The authors thank the financial support provided by the scholarship from the Brazilian agency CNPq (305479/2021-0).This study was financed in part by the Coordination for the Improvement of Higher Education Personnel, Brazil (CAPES), Finance Code 001.
of the closed circuit used as a test bench is shown in F 1.

Figure 1 .
Figure 1.Schematic diagram of the test bench.

Figure 1 .
Figure 1.Schematic diagram of the test bench.

Figure 2 .
Figure 2. Monitoring the fluid density using the Coriolis flow meter after changing the BSW level.

Figure 2 .
Figure 2. Monitoring the fluid density using the Coriolis flow meter after changing the BSW level.

Figure 3 .
Figure 3. Pareto chart for the estimated standardized effects of the ultrasonic flow meter.

Figure 3 .
Figure 3. Pareto chart for the estimated standardized effects of the ultrasonic flow meter.

Figure 4 .Figure 4 .
Figure 4. (A) Response surface and (B) contour plot of the error of the ultrasonic flow meter in the BSW condition equal to 0% v/v.

Figure 4 .
Figure 4. (A) Response surface and (B) contour plot of the error of the ultrasonic flow meter in the BSW condition equal to 0% v/v.

Figure 5 .Figure 5 .
Figure 5. (A) Response surface and (B) contour plot of the error of the ultrasonic flow meter in the BSW condition equal to 5% v/v.

Figure 4 .
Figure 4. (A) Response surface and (B) contour plot of the error of the ultrasonic flow meter in the BSW condition equal to 0% v/v.

Figure 5 .Figure 6 .
Figure 5. (A) Response surface and (B) contour plot of the error of the ultrasonic flow meter in the BSW condition equal to 5% v/v.

Figure 7 .
Figure 7. Pareto chart for the estimated standardized effects for meter in flow rate of Coriolis type.

Figure 7 .
Figure 7. Pareto chart for the estimated standardized effects for meter in flow rate of Coriolis type.

Figure 8 .Figure 8 .
Figure 8. (A) Response surface and (B) contour plot of the error of the Coriolis meter in the BSW condition equal to 0% v/v.

Figure 8 .
Figure 8. (A) Response surface and (B) contour plot of the error of the Coriolis meter in the BSW condition equal to 0% v/v.

Figure 9 .Figure 9 . 20 Figure 10 .
Figure 9. (A) Response surface and (B) contour plot of the error of the Coriolis meter in the BSW condition equal to 5% v/v.

Figure 10 .
Figure 10.(A) Response surface and (B) contour plot of the error of the Coriolis meter in the BSW condition equal to 10% v/v.

Figure 11 .
Figure 11.Pareto chart for the estimated standardized effects for positive displacement meter.

Figure 11 .
Figure 11.Pareto chart for the estimated standardized effects for positive displacement meter. .

Figure 12 .
Figure 12. (A) Response surface and (B) contour plot of the error of the positive displacement meter in the BSW condition equal to 0% v/v.

Figure 13 .
Figure 13.(A) Response surface and (B) contour plot of the error of the positive displacement meter in the BSW condition equal to 5% v/v.

Figure 12 .
Figure 12. (A) Response surface and (B) contour plot of the error of the positive displacement meter in the BSW condition equal to 0% v/v. .

Figure 12 .
Figure 12. (A) Response surface and (B) contour plot of the error of the positive displacement meter in the BSW condition equal to 0% v/v.

Figure 13 .
Figure 13.(A) Response surface and (B) contour plot of the error of the positive displacement meter in the BSW condition equal to 5% v/v.

Figure 13 .
Figure 13.(A) Response surface and (B) contour plot of the error of the positive displacement meter in the BSW condition equal to 5% v/v.

Figure 14 .
Figure 14.(A) Response surface and (B) contour plot of the error of the positive displacement meter in the BSW condition equal to 10% v/v.

Figure 14 .
Figure 14.(A) Response surface and (B) contour plot of the error of the positive displacement meter in the BSW condition equal to 10% v/v.

Table 1 .
Factors correction for events in failure in the framework of petroleum.

Table 2 .
The headquarters experimental generic is for one BBD with three control factors and three center points (C).

Table 3 .
Evaluation of the flow meter rate in petroleum.

Table 4 .
Validation from the bench in test with an ultrasonic meter.

Table 5 .
Validation from the bench in test with a Coriolis meter.

Table 6 .
Validation from the bench in test with a positive displacement meter (meter standard).

Table 7 .
Experimental conditions for flow meters.

Table 8 .
Experimental matrix and results obtained for the ultrasonic flow meter.

Table 9 .
ANOVA for the ultrasonic flow meter.

Table 10 .
Ultrasonic flowmeter error due to increased BSW.

Table 11 .
Experimental matrix and results obtained for the Coriolis flow meter.

Table 12 .
ANOVA for Coriolis flow meter.

Table 13 .
Experimental matrix and results obtained for the positive displacement type flow meter.

Table 13 .
Experimental matrix and results obtained for the positive displacement type flow meter.

Table 14 .
ANOVA for positive displacement type flow meter.

Table 14 .
ANOVA for positive displacement type flow meter.