Design and Optimization of Annular Flow Electromagnetic Measurement System for Drilling Engineering

Using the downhole annular ﬂ ow measurement system to get real-time information of downhole annular ﬂ ow is the core and foundation of downhole micro ﬂ ux control drilling technology. The research work of electromagnetic ﬂ owmeter in recent years creates a challenge to the design of downhole annular ﬂ ow measurement. This paper proposes a design and optimization of annular ﬂ ow electromagnetic measurement system for drilling engineering based on the ﬁ nite element method. Firstly, the annular ﬂ ow measuring and optimization principle are described. Secondly, a simulation model of an annular ﬂ ow electromagnetic measurement system with two pairs of coil is built based on the fundamental equation of electromagnetic ﬂ owmeter by COMSOL. Thirdly, simulations of the structure of excitation system of the measurement system are carried out, and simulations of the size of the electrode ’ s radius are also carried out based on the optimized structure, and then all the simulation results are analyzed to evaluate the optimization e ﬀ ect based on the evaluation indexes. The simulation results show that optimized shapes of the excitation system and electrode size can yield a better performance in the annular ﬂ ow measurement.


Introduction
In recent decades, oil and gas exploration is being carried out in some extremely harsh and challenging environmental conditions [1].Drilling safety issues become increasingly prominent when exploring complex and deep formations.Downhole microflux control drilling technology can effectively solve drilling accident such as kick and lose in narrow density window drilling scenarios, so using the measurement system to get the real-time information of downhole annular flow is the core and foundation of downhole microflux control drilling technology.
A huge array of flow technology options is on offer which provides options in selecting the correct annular flow measurement for the application of drilling engineering.A broad range of factors regarding the special environment of downhole drilling, such as downhole space, velocity profile, temperature, and fluid properties, should be considered.Electromagnetic measurement has the advantages of simple structure, no moving parts, and no obstruction of fluid flow throttle parts.Also, the flow path does not cause any additional pressure loss, and it does not cause wear or blockage, in particular when measuring slurry with solid particles, sewage and other liquid-solid two-phase bodies, or a variety of viscous slurry, and so on.In addition, because the structure has no moving parts, so any corrosion will be attached to the insulation lining.After selecting corrosion-resistant electrode material with a very good corrosion resistance, it can be used for a variety of corrosive media measurements.In 2017, Liang et al. propose a new method for an annular flow measurement system based on the electromagnetic induction principle [2].However, this paper does not describe how to design and optimize the downhole annular flow electromagnetic measurement system.Based on the above reasons, this study mainly focuses on the design and optimization of downhole annular flow electromagnetic measurement system for drilling engineering.
For the traditional electromagnetic flowmeter used in a round pipe, the signal voltage is dependent on the average flow velocity, the magnetic flux density, and the pipe diameter.The signal voltage is expected to be linearly related to the average flow velocity if the magnetic field is a uniform magnetic field.In this ideal case, the flow rate can be considered to be immune to the velocity profile of the pipe flow, especially when the flow has been fully developed.However, it is difficult for the annular flow electromagnetic measurement system to yield a uniform magnetic field for an annular flow path, and also, the velocity profile of the annular flow also cannot be considered axisymmetrically distributed under this special drilling environment.These effects affect the accuracy of the annular flow electromagnetic measurement system.According to Bevie's vector weight function theory [3], if the result of the magnetic flux density cross-product density of virtual current is constant, the annular flow electromagnetic measurement system can be considered to be immune to the velocity profile of the annular flow.When the electrode and the structure of the flow path are fixed, the density of the virtual current is fixed.So the shape of the excitation system can be derived based on this constant condition.Similarly, when the shape of excitation system and the structure of flow path are fixed, the density of virtual current only depends on the electrodes.In this case, the size of the electrode was selected to be optimized to reduce the distortion caused by the velocity profile.
In the past few decades, great efforts about the structure of excitation systems and electrodes have been made to reduce the effects of the velocity profile.Horner B improved the measurement accuracy of flow measurement by increasing the number of electrodes [4,5].Michalski et al. investigated the design of the coils of an electromagnetic flowmeter in 1998 [6].Wang et al. analyzed the relationship between velocity profile and distribution of induced potential for an electromagnetic flow meter in 2007 [7].An optimum excitation coil for an open channel electromagnetic flowmeter was reported by Michalski et al. in 2001 [8].Vieira et al. developed an enhanced ellipsoid method for electromagnetic device optimization and design in 2010 [9].At the same time, some simulation methods were designed with the development of finite element software for multifields.Lim and Choong analyzed the relative errors in evaluating the electromagnetic flowmeter signal using the weight function method and the finite volume method [10].Michalski et al. applied 3D approach to designing the excitation coil of an electromagnetic flowmeter in 2012 [11].Yin et al. investigated the theoretical and numerical approaches to the sensitivity calculation of a novel contactless inductive flow tomography [12].Although great efforts about the structure of the excitation system and electrodes have been made to reduce the effects of velocity profile, most of the design and optimization was focused on the round pipe.The main objective of this paper is to design and optimize the downhole annular flow electromagnetic measurement system with two pairs of electrodes using the method of numerical finite element analysis.

Background Theory of the Annular Four-Electrode Electromagnetic Flowmeter
Following Faraday's law, the flow of a conductive liquid through a magnetic field will cause a voltage signal to be sensed by electrodes located on the flow pipe walls.Faraday's formula can be expressed as where E is the signal voltage in a conductor, v is the average flow velocity, B is the magnetic flux density, and D is the pipe diameter.
Equation (1) indicates that the signal voltage is dependent on the average flow velocity, the magnetic flux density, and the pipe diameter.
The signal voltage is expected to be linearly related to the average flow velocity if the magnetic field is a uniform magnetic field.So (1) can be rewritten as Here, C represents a constant.However, the annular flow electromagnetic measurement system is difficult to yield a uniform magnetic field for an annular flow path, and the velocity profile of the annular flow also cannot be considered axisymmetrically distributed under this special drilling environment.This affects the accuracy of the annular flow electromagnetic measurement system.So the design and optimization of the downhole annular flow electromagnetic measurement system with two pairs of electrodes cannot be investigated based on the traditional Faraday theory.
According to Bevir's theory [3], the theoretical expression of the induced potential in annular flow electromagnetic measurement system can be given as a volume integral of the weight function vector W and the annular flow velocity ν as follows: where W is given by W = B × j .
Here, U is the induced potential, W is the weight function vector, ν is the velocity of the annular flow, B is the magnetic flux vector, j is the virtual current vector, and τ is the integration of annular volume.The velocity of the annular flow weight function vector W is dependent on the magnetic flux density vector B and virtual current vector j .
According to Bevie's vector weight function theory, if the result of the magnetic flux density cross-product density of virtual current is constant, the annular flow electromagnetic measurement system can be considered to be immune to the velocity profile of the annular flow.When the electrode and the structure of the flow path are fixed, the density of virtual current is fixed.Thus, the shape of excitation system can be derived based on this constant condition.Similarly, when the shape of the excitation system and the structure of the 2 Journal of Sensors flow path are fixed, the density of the virtual current only depends on the electrodes.In this case, the radius of the electrode was selected to be optimized to reduce the distortion caused by the velocity profile.
Figure 1 shows a schematic diagram of the downhole annular four-electrode electromagnetic flowmeter and the flow path.A1, A2, B1, and B2 represent four electrodes placed on the outer wall at intervals of 2γ inside the annular flow path (γ = 45 °while in simulation), a and b are the distance of the inner and outer surfaces to the center, and ∑1 and ∑2 represent the inner and outer surfaces of the annular flow path.
To simplify the boundary conditions, we suppose that the electronic conductivity of the outer surface of the annular flow path is much smaller than the flow.The partial differential equation and the boundary conditions of the outer surface of the annular flow path and electrodes can be written as follows: By using the segregation variable method, (4) can be expressed as When K = n 2 , by solving differential equations, the results can be written as 6) can be transformed into the following equation: If A1 and B1 are the positive electrodes and A2 and B2 are the negative electrodes, the boundary conditions of the inner surface of the annular flow path and electrodes can be written as follows: Based on boundary condition (8), we have Here, we suppose that R = r/b and τ = a/b, so G can be written as 3

Journal of Sensors
The density of the virtual current can be derived using the gradient method based on (11).Theoretically, the shape and the parameters of the exciting coil can be solved based on the constant condition of the weight function vector.However, it is very complicated and difficult to derive the shape and the parameters of the exciting coil by the analytical method.Even if the result is solved, it is also complicated and difficult for the coil to wind.In recent years, numerical simulation has become an important process in a new system design.These difficulties can be worked out by the finite element software for multifields.

Optimization of the Annular Four-Electrode Electromagnetic Flowmeter
3.1.Optimization of Excitation Structure.To avoid winding the complicated coil, the iron core was introduced to the excitation part of the annular flow electromagnetic measurement system.Considering the special downhole environment, the schematic diagram of the excitation part of annular flow electromagnetic measurement system is shown in Figure 2. Here, the inner radius of the system is 4 cm, the outer radius is 9 cm, and the inner and outer wall thickness of the system is 1 cm.In the annular flow electromagnetic measurement system, four excitation coils were positioned to generate a magnetic field.Four pick-up electrodes were located around the diameter of the external system wall, so that the flow-induced voltage could be detected.The 3D model of the annular four-electrode flow electromagnetic measurement system and the flow path by COMSOL are shown in Figure 3.
In this paper, the structure of the iron core and the size of the coil among the annular domain between the inner and outer surfaces of the system were investigated, and an optimized excitation structure was designed.Therefore, the width of the core and the height of core protrusion could be used as variables.In the simulations, the round electrodes were chosen with a radius of 0.7 cm, and the length of the iron core was 10 cm.The excitation coils were modified with the width of the core changing from 2 cm to 5 cm and the height of core protrusion from 5.5 cm to 8.5 cm.The width of the core was incremented every 1 cm, while the height of core protrusion every 1 cm.
To validate the effectiveness of the optimization, some performances of the selected optimum annular flow electromagnetic measurement system were compared.Compared simulation results of the magnetic flux vector and weight function vector are shown from Figures 4-7.
Figures 4-7 show that the distribution of the magnetic flux vector and weight function vector changes with the width of the core and height of core protrusion changing.Considering the important role of weight function vector in optimization, the weight function vector distributions from Figures 4-7 are analyzed closely.The analysis results show that the weight function vector gets the smallest distribution range (from −0.03 to 0.32) when core_X is equal to 5.5 cm and Core_w is equal to 2 cm; the weight function vector gets the biggest distribution range (from −0.05 to 0.66) when core_X is equal to 8.5 cm and Core_w is equal to 5 cm; the distribution range of the weight function vector increases while the width of the core and height of core protrusion increase steadily, and increasing the width of the coil is much more pronounced.

The Change of Electrode Radius.
After the structure of the excitation system of the measurement system was optimized and in order to improve the accuracy of the system, some simulations about the radius of the electrodes were used to improve the system performance.The electrodes were modified, with the radius changing from 1 cm to 5 cm, and the radius of the electrodes was incremented every 0.5 cm.Compared simulation results of the virtual current vector and weight function vector are shown in Figures 8 and 9.
Figures 8 and 9 show that the virtual current density and weight function vector changed while the electrode radius changed.Based on Figures 8 and 9, we can find that the weight function vector gets the smallest distribution range (from 0 to 0.03) when the electrode radius is 5 cm; the weight function vector gets the biggest distribution range (from −0.01 to 0.16) when the electrode radius is 1 cm; the distribution range of weight function vector decreases while the electrode radius increases steadily.The ideal excitation system design is to make the weight function vector to be a constant, and this is helpful to improve the measurement accuracy of the system.However, the weight function vector is not easily designed to be a constant.In order to design a magnetic field that makes the weight function vector as constant as possible, it is necessary to define a design quantity that measures the degree of nonuniformity of the weight function vector over the cross-sectional area of the annular domain.The definition of the design quantity was discussed by Dennis and Wyatt in 1972 [13], but owing to integration difficulties, this quantity is not used here.Instead, a similar quantity W σ named the weight function vector standard deviation is used which is defined by the following: where W is the average weight function vector value of all points in the XY plane of the annular while z = 0 and N is the number of the points.

Homogeneity Range Ratio of Weight Function Vector.
To evaluate the homogeneity range in the annular area, the ratio of the homogeneity range is used as an evaluation criterion [14].The equation of this evaluation criterion is as follows:  where W x,y is the weight function vector of the central cross section of the electrode in the annular area, W is the average weight function vector value of all points in the XY plane of the annular while z = 0, and P is a uniform percentage value that depends on actual distribution and which is 30% in this research.If a finite element meets (13), we may regard the area of this finite element to be homogeneous; otherwise, we regard that the area of this finite element is inhomogeneous.If the number of inhomogeneous finite elements is N 1 and the number of homogeneous finite elements is N 2 , the homogeneity range ratio can be defined as follows: 4.1.3.The Coefficient of Variation of Weight Function Vector.
The coefficient of variation is a useful statistic for comparing the degree of variation from one data series to another, even if the means are drastically different from one another.A coefficient of variation of the weight function vector is a statistical measure of the dispersion of the weight function vector data points in a data series around the mean.It is calculated as follows: where W σ is the standard deviation of the weight function vector and W is the average weight function vector value of all points.
4.1.4.The Sensitivity of Output Voltage.The fourth performance criterion is the sensitivity of the output voltage which gives the strength of the system output.In this research, we assure that the turns and current of the coil are constant, and we get the sensitivity of the output voltage while the average flow velocity is 1 meter per second.
During the optimal design and quantitative evaluation process, smaller standard deviation of the weight function vector is better, as is a bigger homogeneity range ratio of the weight function vector is better, smaller coefficient of variation of the weight function vector is better, and bigger sensitivity of the output voltage is better.In these four evaluation indexes, the standard deviation of the vector weight function is the most important evaluation index, and the other 3 indexes are auxiliary reference indexes.In Figure 10, Core_w represents the width of the core and core_X represents the height of core protrusion.When the width of the core is fixed, the weight function vector standard deviation increases while the height of core protrusion increases steadily.When the height of core protrusion is fixed, the weight function vector standard deviation increases while the position of the width of the core increases steadily.
Using Figure 10, we can find the minimum standard deviation of the weight function vector while core_X is 5.5 cm and Core_w is 2 cm.Similarly, based on the simulation data and mathematical operations, the variation curve of the homogeneity range ratio, coefficient of variation of the weight function vector, and the sensitivity of the output voltage are shown from Figures 11-13.
It can be seen from Figures 11-13 that when the width of the core and height of core protrusion change, the evaluation indexes change.When core_X is equal to 5.5 cm and Core_w     Journal of Sensors is equal to 2 cm, the homogeneity range ratio is relatively large, but the coefficient of variation of the weight function is large and the sensitivity of output voltage is small.Considering that the standard deviation of the vector weight function is the most important evaluation index, optimum results of annular flow electromagnetic measurement system under the double pairs of coil excitation structure and current finite structure space can be obtained, with core_X = 5.5 cm and Core_w = 2 cm being the best coil excitation structures.However, from the range of numerical magnitude of evaluation indexes, we can find that changing the width of the core and height of core protrusion has a limited effect on improving the system measurement performance.

Analysis of Electrode Radius Effect
. In order to intuitively analyze and compare the effect of excitation system optimization, the value of the weight function vector standard deviation can be analyzed using the definition of the design quantity.To validate the effectiveness of the optimization, the results of the weight function vector standard deviation are compared in Figure 14.The data result in Figure 14 shows that the weight function vector standard deviation decreases while the electrode radius increases steadily, so the big radius   Similarly, based on simulation data and mathematical operations, the variation curve of the homogeneity range ratio, the coefficient of variation of the weight function vector, and the sensitivity of output voltage are shown from Figures 15-17.It is clear that with the increases of the electrode radius from 0.5 cm to 6 cm, the homogeneity range ratio keeps increasing from 6.43% to 79.49%, the coefficient of variation of the weight function vector keeps decreasing from 2410 to 498, and the sensitivity of the output voltage increases first from 3.07E − 05 V to 3.22E − 05 V and then decreases from 3.22E − 05 V to 2.81E − 05 V. Thus, the electrode radius need to be chosen according to the requirement of the system performance, and the electrode radius of 4.5 cm has a better optimum comprehensive effect during the annular flow measuring process.

Conclusions
This paper described the design and optimization of annular flow electromagnetic measurement system for drilling engineering.The following conclusions can be drawn according to the above-mentioned analysis: (1) The theory of annular flow electromagnetic measuring and optimization principles were described, and an annular flow electromagnetic measurement system simulation model was built based on this theory by COMSOL.Journal of Sensors (2) Simulations on the structure of excitation system of measurement system were carried out, as well as some simulations on the size of the electrode's radius were also carried based on the optimized structure, and then all the simulation results were analyzed to evaluate the optimization effects based on the evaluation indexes.The simulation results showed that the optimized structure of excitation system and electrode radius size can yield a better performance during the annular flow measurement process.

Figure 1 :
Figure 1: The schematic diagram of the annular four-electrode flow electromagnetic measurement system and the flow path.

Figure 2 :
Figure 2: The schematic diagram of the excitation part of the annular flow electromagnetic measurement system.

Figure 3 : 1 .
Figure 3: The 3D model of the annular four-electrode flow electromagnetic measurement system and the flow path.

Figure 4 :
Figure 4: The simulation diagram of the magnetic flux vector and weight function vector with core width equal to 2 cm.

( b )Figure 5 :
Figure 5: The simulation diagram of the magnetic flux vector and weight function vector with core width equal to 3 cm.

Figure 6 :
Figure 6: The simulation diagram of the magnetic flux vector and weight function vector with core width equal to 4 cm.

Figure 7 :
Figure 7: The simulation diagram of the magnetic flux vector and weight function vector with core width equal to 5 cm.
(f) Radius equal to 5 cm in the YZ plane

Figure 8 :
Figure 8: The simulation diagram of virtual current density with electrode radius equal to 1 cm, 3 cm, and 5 cm separately.

8
Radius equal to 5 cm in the YZ plane

Figure 9 :
Figure 9: The simulation diagram of weight function vector with electrode radius equal to 1 cm, 3 cm, and 5 cm separately.

Figure 10 :
Figure 10: The curve of weight function vector standard deviation and structure of excitation system.

Figure 11 :Figure 12 :Figure 13 :
Figure 11: The curve of the homogeneity range ratio and structure of the excitation system.

Figure 14 :Figure 15 :Figure 16 :
Figure 14: The curve of the weight function vector and electrode radius.

Figure 17 :
Figure 17: The curve of the sensitivity of the output voltage and electrode radius.
(b) Radius equal to 1 cm in the YZ plane Radius equal to 3 cm in the YZ plane