Next Article in Journal
A Hybrid Motion Estimation for Video Stabilization Based on an IMU Sensor
Next Article in Special Issue
High-Accuracy Calibration Based on Linearity Adjustment for Eddy Current Displacement Sensor
Previous Article in Journal
Semi-Supervised Generative Adversarial Nets with Multiple Generators for SAR Image Recognition
Previous Article in Special Issue
Experimental Study on Residual Bending Strength of Corroded Reinforced Concrete Beam Based on Micromagnetic Sensor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Uniform Linear Multi-Coil Array-Based Borehole Transient Electromagnetic System for Non-Destructive Evaluations of Downhole Casings

1
Key Laboratory of Education Ministry for Photoelectric Logging and Detecting of Oil and Gas, Xi’an Shiyou University, Xi’an 710065, China
2
School of Marine Science and Technology, Northwestern Polytechnical University, Xi’an 710072, China
3
Petroleum Engineering Technology Research Institute, Shengli Oilfield Company, SINOPEC, Dongying 257000, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(8), 2707; https://doi.org/10.3390/s18082707
Submission received: 5 July 2018 / Revised: 15 August 2018 / Accepted: 15 August 2018 / Published: 17 August 2018
(This article belongs to the Special Issue Magnetic Sensors)

Abstract

:
Borehole transient electromagnetic (TEM) techniques have been proven to be efficient for nondestructive evaluations (NDEs) of metal casings using eddy-current properties. However, physical limitations and bad borehole conditions restrict the use of eddy-current sensors, which makes downhole casing inspections very different from those of conventional NDE systems. In this paper, we present a uniform linear multi-coil array-based borehole TEM system for NDEs of downhole casings. On the basis of the borehole TEM signal model, a numerical multi-coil array approach using the Gauss–Legendre quadrature is derived. The TEM response can be divided into two independent parts related to the transmitting-receiving distance (TRD) and the observation time and casing thickness. Using this property, the signal received by the multi-coil array is weighted to cancel the influence of the TRDs of the different array elements to obtain the optimal response according to the linearly constrained minimum variance criterion, which can be shown to be identical to that of achieving the maximum signal-to-noise ratio. The effectiveness of the proposed method was verified by applying the uniform linear multi-coil array to a borehole TEM system for NDEs of oil-well casings. Field experiments were conducted, and the results demonstrate the effectiveness of the proposed method.

1. Introduction

Transient electromagnetic (TEM) techniques have gained much attention over the past several decades owing to their use in a wide range of applications in geophysical prospecting, such as mineral and petroleum exploration [1], geotechnical and environmental investigations [2], and fundamental studies of stress and petrophysics [3]. In the field of borehole detection, TEM systems permit the rapid and accurate acquisition of broad-frequency-range data concerning the electrical and geometrical parameters of each borehole cylindrical layer due to their significant feature of the accessibility to targets [4,5]. This technique, which is also known as transient (pulsed) eddy-current testing [6,7,8], enables highly effective nondestructive evaluations (NDEs) of downhole casings [9,10]. However, unlike surface measurements that use large loop coils [1] or high performance magnetometers [11,12], the structure of borehole TEM sensors is strongly restricted by the physical limitations of the downhole conditions, including limited sensor sizes, low sampling performance due to high temperatures, and cumbrous metal tool housings against the high pressure, thereby, reducing the signal strength, as well as signal-to-noise ratio (SNR) of the borehole TEM measurements, and making NDEs of downhole casings very different.
Much research has previously been undertaken with respect to the design and use of eddy-current sensors to improve the NDE performance for downhole casings. Considering the limited space and the high temperatures, which result in a weak TEM response and low SNR, multi-turn coil approaches have been proposed to enhance the signal strength and to improve the inspection performance. In Ref. [13], the eddy-current NDE of metal pipes with an arbitrary position of the coil sensors was theoretically and experimentally investigated. To inspect different crack shapes, a combination of longitudinal and transverse sensors composed of coaxial multi-turn transmitting and receiving coils wound around a magnetic or air core was employed [14,15], and it is shown in [15] that the longitudinal sensor performs better than the transverse sensor for detection of large area. In Ref. [16], multi-pipe strings, which are required to protect the oil and gas wells from by-products for safety considerations, but which make NDEs more difficult in terms of data interpretation [15,17], were inspected using the eddy-current diffusion properties of a TEM system with a multi-turn coil-based auxiliary sensor in longitudinal direction. The above methods have been proven to be effective for improving the NDE performance in downhole casings. However, restricted by the small radial dimensions of the wellbore, the large number of turns of longitudinal multi-turn coil sensors will also increase the length of sensor, as well as the transmitting–receiving distance (TRD), which will lead to serious distortions in the receiving signal model, and strongly influence the accuracy of NDEs.
At present, there have been minimal previous attempts to eliminate the influence of TRD on downhole measurements, which is difficult to suppress using the traditional multi-turn coil sensor with single-receiver [13,14,15,16,17], since the TRD is coupled with the electrical and geometrical parameters in only one output. On the basis of Doll’s theory in induction logging system [18], the TRD was replaced by two distances from the transmitting and receiving coils to “an elementary ring”. Using this approximate theory, the multi-coil array sensors were successfully used for magnetic field focusing [19] to improve borehole detection performance [20,21]. This method is valid in an open-hole with sufficiently resistive medium by neglecting the skin effect; however, in a cased-hole with a high-conductivity metal casing, Doll’s theory may not be satisfied. In Ref. [22], a two-coil structure was proposed in the case of motion measurements, where a reference multi-turn receiver was used to cancel the effect of the magnetic background noise in the main channel by employing the correlation of the two receivers with different TRDs. Although the relationship between the two receivers was not investigated thoroughly enough to be extended to a multi-coil array sensor, the use of multiple receivers would offer a new way for solving the problem caused by TRD.
In this paper, we present a uniform linear multi-coil array-based borehole TEM system for NDEs of downhole casings. Using a Gauss–Legendre quadrature-based numerical approach to the borehole TEM signal model, we apply a weight to the signals received by the multi-coil array according to the linearly constrained minimum variance (LCMV) criterion [23]. It is shown that the influence of TRD can be eliminated using LCMV-based array signal processing, as compared to the traditional multi-turn coil sensor with the same total number of receiving coil turns, where the maximum SNR can also be achieved. The effectiveness of the proposed system was verified by applying it to an oil borehole TEM system to inspect an oil-well casing.
The rest of this paper is organized as follows. The borehole TEM signal model based on a linear uniform multi-coil array with a coaxial transmitting coil and multiple receiving coils is presented in Section 2. In Section 3, we present the Gauss–Legendre quadrature-based numerical approach to the borehole TEM signal model. The LCMV-based array signal processing method for the linear uniform multi-coil array of the borehole TEM system and its performance analyses are presented in Section 4. The experimental and simulation results are discussed in Section 5. Finally, we conclude the paper in Section 6.

2. Borehole TEM Signal Model

Consider a uniform linear multi-coil array-based borehole TEM system equipped with a single transmitter and M receivers with an inter-element spacing of Δz and TRDs ranging from z1 to zM, which are comprised of coaxial coils wound around a soft magnetic core in a cylindrically layered medium. The multi-cylindrically layered structures of the borehole TEM system including the electrical and geometrical parameters of the jth layer (μj, εj, σj, and rj with j = 1,2,…,J) are illustrated in Figure 1. We consider the magnetic core to be the innermost layer. The transmitter and the M receivers are located in the second layer, with their number of turns given by NT and NR, respectively, where each receiver has the same number of turns. In this paper, M receivers are utilized instead of the traditional single receiver [13,14,15,16,17] with the same total number of receiving coil turns of MNR. Thereby, with a large value of M, NR, as well as NT, can be small enough to ensure that the change in the TRD for each single receiver can be ignored and treated approximately as a single value. It is assumed that the diameters of all the coils are sufficiently small so that the source region contains only the second layer, and the induced electromagnetic force (EMF) of the receiving coils is only related to the vertical component of the magnetic field of the first layer (the magnetic core). Moreover, all of the other layers, such as the liquid mud, casing, and formation, are regarded as source-free regions. In addition, to measure the thickness of the metal pipes, we assume that the electrical parameters of all the layers and the inner radius of the metal pipes are fixed.
The response of a TEM system in such a multi-cylindrically layered geometry consists of reflection and transmission components with both standing and outgoing waves. Using the vector potential A, the homogeneous and inhomogeneous Helmholtz equations are given by [18]:
2 A 2 + k 2 2 A 2 = J e ,
2 A j + k j 2 A j = 0 ,    j 2 ,
where kj2 = μjεjω2jσjω and Je denotes the electrical source. With the introduction of the variables xj and λ, which satisfy xj2 = λ2kj2, the vector potential A can be calculated by solving the Helmholtz equations. Note that directional measurements cannot be achieved using the proposed borehole TEM system due to the cylindrical symmetry of the model in Figure 1. For the primary field [18], the solution of the inhomogeneous Helmholtz equation in the source region has the form:
A 2 ( ω , z , r ) = { I T r 0 π 0 K 1 ( x 2 r ) I 1 ( x 2 r 0 ) cos ( λ z ) d λ , r > r 0 I T r 0 π 0 K 1 ( x 2 r 0 ) I 1 ( x 2 r ) cos ( λ z ) d λ , r < r 0 ,
where IT denotes the transmitting current, z denotes the distance between the excitation (transmitter) and observation (receiver) positions along the borehole axis, r0 is the radius of the transmitter, and I1(·) and K1(·) are modified Bessel functions of the first and second kind of order one, respectively. For the secondary field [18], the solution of the homogeneous Helmholtz equation has the form:
A j ( ω , z , r ) = I T r 0 π 0 [ C j I 1 ( x j r j ) + D j K 1 ( x j r j ) ] cos ( λ z ) d λ ,
where Cj and Dj denote the reflection and transmission coefficients, respectively, which are related to the geometrical and electrical parameters of all the layers. Note that CJ and D1 are zero due to the absence of transmission and reflection in the innermost and outermost layers, respectively. Then, Cj and Dj can be calculated using the boundary conditions of the multilayered cylindrical structures [18]. In the source region, the response contains both primary and secondary fields, while in the source-free region, only the secondary field is involved. Therefore, the vertical component of the magnetic field [18] in the innermost layer with radius r (0 < r < r1) can be calculated by combining the solutions of the Helmholtz equations, and can be written as:
H z 1 ( ω , z , d ) = N T r 1 I T π 0 x 1 C 1 I 0 ( x 1 r ) cos ( λ z ) d λ ,
where I0(·) denotes a modified Bessel function of the first kind of order zero, and d denotes the thickness of the metal casings, and can be expressed by r5r4 if only one casing is used; all of these geometrical parameters are included in C1. On the basis of the above models [18], the induced EMF in the frequency domain in the mth receiver of multi-coil array with a TRD of zm can then be expressed as:
U m ( ω , z m , d m ) = i ω ξ 0 r 1 0 2 π r f ( λ , r , ω , d m ) cos ( λ z m ) d λ d r ,
where ξ = μ1NRNTIT/π, f(λ,r,ω,dm) = x1C1I0(x1r), and dm denotes the thickness of the metal casing with respect to the observation position of the mth receiver. Given a ramp signal with a turn-off time of tof, the induced EMF Um(t,zm,dm) can be obtained by converting Equation (6) into the time domain. Using the Gaver–Stehfest inverse Laplace transform [24] with S stages as an example, we obtain:
U m ( t , z m , d m ) = ln 2 t s = 1 S D s e ( s   ln 2 t of ) / t 1 t of ( s   ln 2 / t ) 2 U m ( s   ln 2 / i t , z m , d m ) ,
where = sln2/t and t and Ds denote the observation time and the integral coefficient of the Gaver–Stehfest inverse Laplace transform, respectively. We can see that the induced EMF of the borehole TEM system is related not only to the observation time and the thickness of the downhole casing, but also to the TRD of the mth receiver. In other words, the three variables in Equation (7) are coupled together in the TEM response. Of these three variables, it is the coupling between t and d, also known as the eddy-current diffusion or the time-domain property, that offers several benefits to the TEM system to achieve better performance [16,25]. Conversely, when a single receiver with a large number of turns or a multi-coil array is used, the TRD will strongly influence the interpretation of the NDE, where the change in the TRD cannot be ignored and must be compensated for to avoid model distortions in the borehole TEM system.

3. Numerical Approximation of the Borehole TEM Signal Model

In previous studies, multi-turn coils and multi-coil arrays have been used to improve the inspection performance of borehole TEM systems. However, as shown in Section 2, the effect of the TRD on the TEM response makes the interpretation of NDEs more difficult. In Ref. [22], even though the correlation between two channels with different TRDs is employed for noise cancellation with a registration matrix, this correlation was not thoroughly investigated because the two parts with respect to the TRD and geometrical–electrical parameters are still coupled to each other. In this section, we decouple the TRD and the other two variables in Equation (7) using a Gauss–Legendre quadrature-based [26] numerical approximation for the borehole TEM signal model. Considering the monotonically decreasing characteristic of the integrand in Equation (6) with respect to the modified Bessel functions, the upper limit of the infinite integral can be reduced to an approximately limited value of λ0 (λ0 = 6000 in this paper). Then, by adjusting the integration interval of the two integrations, we can rewrite Equation (6) as:
U m ( ω , z m , d m ) = i ω ξ χ 1 1 1 1 r + 1 2 f ( λ 0 λ + 1 2 , r 1 r + 1 2 , ω , d m ) cos ( λ 0 z m λ + 1 2 ) d λ d r ,
where χ = πr12λ0/2. On the basis of the Gauss–Legendre quadrature equation, the above two integrations can both be expanded as multi-stage Legendre polynomials such that:
U m ( ω , z m , d m ) = i ω ξ χ q = 1 Q A q B q + 1 2 p = 1 P A p f ( λ 0 B p + 1 2 , r 1 B q + 1 2 , ω , d m ) cos ( λ 0 z m B p + 1 2 ) ,
where P and Q represent the number of stages of the two Legendre polynomials, with their quadrature coefficient and zero-point denoted by A and B, respectively. Substituting Equation (9) into Equation (7), we have:
U m ( t , z m , d m ) = ξ χ t of s = 1 S q = 1 Q p = 1 P g s , q , p ( t , d m ) v p ( z m ) ,
where
g s , q , p ( t , d m ) = D s e ( s ln 2 t of ) / t 1 s A q B q + 1 2 A p f ( λ 0 B p + 1 2 , r 1 B q + 1 2 , s   ln 2 / i t , d m ) ,
v p ( z m ) = cos ( λ 0 z m B p + 1 2 ) ,
In matrix form, Equation (10) can be expressed as:
U m ( t , z m , d m ) = ξ χ t of x ( z m ) g T ( t , d m ) ,
where
x ( z m ) = [ v ( z m ) , , v ( z m ) S Q ] 1 × S Q P ,
v ( z m ) = [ v 1 ( z m ) , , v P ( z m ) ] 1 × P ,
g ( t , d m ) = [ g 1 , 1 , 1 ( t , d m ) , , g S , Q , P ( t , d m ) ] 1 × S Q P .
In Ref. [13], x(zm) and g(t,dm) represent the effect of the TRD and the geometrical–electrical parameters of the multi-cylindrical structure, respectively. Using the Gauss–Legendre quadrature-based numerical approach, the two parts can be approximately separated and are independent from each other, which means that the influence of the TRD on the induced EMFs can be eliminated by cancelling x(zm). To formulate a clear expression, we assume that the induced EMFs of each coil are discretely sampled with a sampling length L and a sampling interval ∆t by a 16-bit analog-to-digital converter; then, the induced EMF of the mth receiver with sampling time tl and TRD zm can be rewritten as:
U m , l ( t l , z m , d m ) = ξ χ t of x ( z m ) g T ( t 1 + ( l 1 ) Δ t , d m ) .
The vector form of the induced EMF of the mth receiver can be expressed as:
U m , 1 L = [ U m , 1 ( t 1 , z m , d m ) U m , 2 ( t 2 , z m , d m ) U m , L ( t L , z m , d m ) ] = ξ χ t of x ( z m ) [ g T ( t 1 , d m ) g T ( t 1 + 1 Δ t , d m ) g T ( t 1 + ( L 1 ) Δ t , d m ) ] 1 × L ,
and the matrix form of the received signal of the uniform linear multi-coil array can be expressed by stacking the induced EMF of each receiver, such that:
U 1 M , 1 L = [ U 1 , 1 L T ( t , z 1 , d 1 ) U 2 , 1 L T ( t , z 2 , d 2 ) U M , 1 L T ( t , z M , d M ) ] T = ξ χ t of [ x ( z 1 ) g T ( t 1 , d 1 ) x ( z 1 ) g T ( t 2 , d 1 ) x ( z 1 ) g T ( t L , d 1 ) x ( z 2 ) g T ( t 1 , d 2 ) x ( z 2 ) g T ( t 2 , d 2 ) x ( z 2 ) g T ( t L , d 2 ) x ( z M ) g T ( t 1 , d M ) x ( z M ) g T ( t 2 , d M ) x ( z M ) g T ( t L , d M ) ] M × L
In this paper, we assume that the length of the multi-coil array is much smaller than the distance between two neighboring thickness changes, which means that the thicknesses of the metal casing with respect to the observation position of each receiving coil are the same along the borehole axis with d1 = d2 =…= dM = d0. Then, the received signal of the multi-coil array can be expressed as:
U 1 M , 1 L = ξ χ t of [ x ( z 1 ) x ( z 2 ) x ( z M ) ] [ g ( t 1 , d 0 ) g ( t 2 , d 0 ) g ( t L , d 0 ) ] T .
Considering the system noise, we rewrite the received signals of the lth sampling time, as well as the lth column of Equation (20) as follows:
U 1 M , l = ξ χ t of X ( z ) g T ( t l , d 0 ) + N ,
where
X ( z ) = [ x ( z 1 ) T x ( z 2 ) T x ( z M ) T ] M × S Q P T ,
N = [ n 1 n 2 n M ] M × 1 T .
N is assumed to follow a Gaussian distribution, and each element of N is independently identically distributed. Obviously, for each sampling time, the difference between each receiver is the TRD. Notably, in Equation (21), the matrix X(z) plays a similar role as the array manifold of a phased array in radar signal processing, and will contaminate the matrix g(tl,d0), which includes the thicknesses of the metal casings to be estimated. Moreover, because the array manifold X(z) is independent of the observation time and the casing thickness, the cancellation of the TRD effect will not impact NDEs of downhole casings, and therefore, offers a new way to process the TEM data collected from the uniform linear multi-coil array with each coil corresponding to different TRDs. In the next section, we use the array signal processing principle to solve this problem on the basis of the proposed numerical approximation.

4. LCMV-Based Multi-Coil Array Weighting for a Borehole TEM System

4.1. LCMV-Based Multi-Coil Array Weighting Model

In Section 3, it was illustrated that X(z) is only related to the TRD for a uniform linear multi-coil array, and acts similarly to an array manifold in a phased-array radar. Therefore, we can apply a weight to the received signal to eliminate the influence of the TRD by cancelling out the array manifold X(z). Under this configuration, we employ the principle of the LCMV criterion to weight the array signals in Equation (21), where a receiving weight vector WR1 is applied to optimize the array data.
y l = ξ χ t of W T X ( z ) g T ( t l , d 0 ) + W T N ,
where yl is the weighted array output of a sampling time tl. Note that, when z = 0, we have x(0) = F, and therefore, FgT(tl, d0) is just the desired part of the received signal of the receiver collocated with the transmitter, where FR1×SQP is a row vector with all SQP elements equal to 1. Then, according to the LCMV criterion-based adaptive digital beamforming method [23], we let the desired array output of Equation (13) be FgT(tl, d0). The constraint should be subject to WTX(z) = F. Moreover, because X(z) is only related to the TRD, the array weight will not change regardless of the sampling time and casing thickness as long as the multi-coil array structure is fixed. Using this property, the array weighting can be extended by applying the weight to the received signal of the uniform linear multi-coil array for all L sampling times such that:
Y = [ y 1 y 2 y L ] 1 × L = W T U 1 M , 1 L ,
where the weighting structure of the proposed method is described in Figure 2.
In this paper, the array weighting can be applied to the output of each receiver by an operation amplifier, where the weight coefficient should be a real value as the cosine but not complex function of X(z). Then, we can directly optimize the received signal to obtain the optimal array output by minimizing the variance of Equation (24), and the optimal weight can be calculated by using the following optimization problem with the LCMV criterion as
{ min W T R U W s . t W T X ( z ) = F ,
where
R U = E { U 1 M , l U 1 M , l T } M × M = X ( z ) E { g ( t l , d 0 ) g T ( t l , d 0 ) } X T ( z ) + R N = δ s 2 X ( z ) X T ( z ) + R N ,
denotes the auto-correlation matrix of U1-M,l and RN is the noise correlation matrix (NCM). Moreover, because each channel has the same number of turns of receiving coils, the NCMs of all the channels are also the same.

4.2. Approximate Optimal Solution for the LCMV Criterion Problem in Borehole TEM System

Using the Lagrange multiplier method [23], the upper optimization problem can be transformed by introducing the Lagrange operator vector β∈C1×SQP, such that:
Ξ = W T R U W + β ( X T ( z ) W F T ) ,
where each of the elements of β satisfy 0 < βsqp < 1 and the optimal weight is:
W = R U 1 X ( z ) ( X T ( z ) R U 1 X ( z ) ) 1 F T .
In Equation (29), the rank of XT(z)RU−1X(z) should be the smaller of SQP and M. In Section 3, considering the numerical approach that converts the finite integral to a matrix multiplication using the Gauss–Legendre quadrature, larger P and Q stage numbers will result in smaller approximation errors. Conversely, the manufacture of the multi-coil array will restrict M. As a result, SQP is usually much larger than M, and XT(z)RU−1X(z) is not invertible when solving the optimization problem. In this paper, the characteristics of the numerical approximation-based borehole TEM model are used to approximately solve this problem. According to Equation (13), the constraint in Equation (26) can be expanded as:
W T X ( z ) = [ w 1 w 2 w M ] T [ v ( z 1 ) v ( z 1 ) v ( z 1 ) v ( z 2 ) v ( z 2 ) v ( z 2 ) v ( z M ) v ( z M ) v ( z M ) S Q ] M × S Q P .
Equation (30) indicates that the column of X(z) is repeated with an interval of P. Substituting Equation (15) into Equation (30), the constraint of the optimization problem in Equation (30) can be converted to:
W T Λ ( z ) = [ w 1 w 2 w M ] T [ v 1 ( z 1 ) , , v P ( z 1 ) v 1 ( z 2 ) , , v P ( z 2 ) v 1 ( z M ) , , v P ( z M ) ] M × P = Φ ,
where Λ∈RM×P and Φ∈R1×P is a row vector with all P elements equal to 1. Now, as long as MP, which can be achieved in the fabrication of the multi-coil array sensor, XT(z)RU−1X(z) can be inverted. Substituting Equation (31) into Equation (29), the optimal weight can then be calculated as:
W = R U 1 Λ ( z ) ( Λ T ( z ) R U 1 Λ ( z ) ) 1 Φ T .

4.3. Performance Analyses

On the basis of the LCMV-based multi-coil array for borehole TEM system, we demonstrated that the influence of the TRD can be effectively eliminated, where the received signals of the uniform linear multi-coil array can be weighted to obtain the optimal response that would be received by a receiver with z = 0. Now, we show performance analyses of the proposed LCMV-based multi-coil array method. According to Equation (26), the SNR of the received signal of the multi-coil array follows:
SNR = W T R U W W T R N W .
Assuming that the NCM can be decomposed as RN = ΓTΓ, where Γ is invertible, Equation (33) can then be rewritten as:
SNR = ( Γ W ) T ( Γ 1 ) T X ( z ) X T ( z ) Γ 1 ( Γ W ) ( Γ W ) T ( Γ W ) .
According to the Schwartz inequality, we have:
( Γ W ) T ( Γ 1 ) T X ( z ) X T ( z ) Γ 1 ( Γ W ) ( Γ W ) T ( Γ W ) ( X T ( z ) Γ 1 ) ( X T ( z ) Γ 1 ) T .
Therefore, the max SNR can only be achieved when (XT(z)ΓΓ−1)Tμ = ΓW and the optimal weight is
W SNR = R N 1 X ( z ) μ ,
where μ∈CSQP×1. Conversely, assuming that the noise and the signal are not correlated, the inversion of RU can be rewritten as:
R U 1 = R N 1 δ s 2 R N 1 X ( z ) X T ( z ) R N 1 1 + δ s 2 X T ( z ) R N 1 X ( z ) .
The optimal weight in Equation (29) can be rewritten as:
W = ( R N 1 δ s 2 R N 1 X ( z ) X T ( z ) R N 1 1 + δ s 2 X T ( z ) R N 1 X ( z ) ) X ( z ) ( X T ( z ) ( R N 1 δ s 2 R N 1 X ( z ) X T ( z ) R N 1 1 + δ s 2 X T ( z ) R N 1 X ( z ) ) X ( z ) ) 1 F = R N 1 X ( z ) ( X T ( z ) R N 1 X ( z ) ) 1 F = R N 1 X ( z ) μ ,
where μ’∈CSQP×1. Comparing Equations (36) and (38), we find that, when the optimization problem is subjected to WTX(z) = F to maximize the SNR, we have μ = μ’, which means that the proposed LCMV-based multi-coil array method is equal to the maximum SNR of the borehole TEM system; therefore, the SNR, as well as the inspection performance of the downhole casings, can be effectively improved.

5. Field Experiments

5.1. Experimental Results

The validity of the uniform linear multi-coil array-based borehole TEM system for downhole casing inspections was confirmed using field experiments. The experiments were conducted in standardized 51/2-inch (with thicknesses of 7.72 mm and outer diameters of 139.7 mm) metal casings. In our experiments, two metal casings with different structures with thickness changes, including increases in the thickness from 1 mm to 4 mm corresponding to lengths from 5 cm to 15 cm in the longitudinal direction, respectively, as shown in Figure 3, were designed and analyzed to illustrate the effectiveness of the proposed method.
The parameters of the multi-coil array sensor and the experiment are shown in Table 1, where the actual outer radius r5 can be calculated by adding the thickness change in Figure 3 to the standardized casing outer radius.
Figure 4 and Figure 5 show the experimental results of the eight receivers of the uniform linear multi-coil array at an early time of 20 ms and a late time of 40 ms in the inspection of the experimental casing structures in Figure 3a,b, respectively. In Figure 4 and Figure 5, the received signals of each receiver are normalized to make them clearer and easier to distinguish and compare. Regardless of the decrease in the SNR with increasing TRD, we can still see that the induced EMFs of the eight receivers at the same sampling time have similar shapes, excepting a longitudinal (borehole axis) shift, which is primarily caused by the different TRDs of the different receivers, and performs just like a “phase shift” in a phased-array radar. Obviously, if the receiving data of the eight receivers are directly summed (the arithmetic mean), the resolution of the casing inspection will be substantially decreased due to the shift in the receiving data along the borehole axis with respect to the different TRDs. However, comparing the early and late time data in Figure 4 and Figure 5, we find that, even though the early time data of the TEM system have better resolution than the late time data [16], the shifting also appears more seriously, which means that the size of this type of shifting is related not only to the TRDs and the inter-element spacing but also to the eddy-current diffusion range. Therefore, it will be difficult to compensate for the shifting via an ordinary curve alignment for all the diffusion times (sampling times). Using the numerical approach in Section 3, the proposed LCMV-based array processing is only related to the TRDs and will not be affected by the sampling time and casing thickness, therefore making the compensation more effective and efficient.

5.2. Analysis and Discussion

In order to demonstrate the effectiveness of the proposed method, the inspection performance of a traditional borehole TEM sensor with a large number of turns in the receiving coils was also analyzed to compare to the proposed uniform linear multi-coil array sensor, where the transmitters of the two sensors were the same, and the receiver of the traditional borehole TEM sensor was approximately substituted using the arithmetic mean of the eight receivers of the linear multi-coil array, so that the number of receiving coil turns of two sensors were also the same. Note that, even though the use of a receiver with a TRD of zero (z = 0, collocated with the transmitter) in the linear multi-coil array can substantially improve the NDE performance, the effectiveness of the proposed method is doubtful in the case of a high SNR measurement, because the response of the collocated receiver will then be nearly the same as the desired array output. Therefore, without loss of generality, the collocated receiver is only used as a comparison to the proposed method, but is not included in the eight receivers of the uniform linear multi-coil array. Nevertheless, note that the collocated receiver, as well as the symmetric array with respect to the transmitter, will be of great importance in improving the NDE performance in a real system.
Figure 6 and Figure 7 compare five cases of array output for the experimental results. The first case shows the original data of the receiver collocated with the transmitter with z = 0, where no shifting exists, regardless of the sampling time and motion velocity. In the second case, the induced EMFs of the eight receivers are directly summed to obtain the arithmetic mean. For comparison, we show the weighted array output using the proposed LCMV-based multi-coil array method in the last three cases, where eight receivers (Receivers 1–8) with inter-element spacings of 2 cm, four receivers (Receivers 1–4) with inter-element spacings of 2 cm, and four receivers (Receivers 2, 4, 6, and 8) with inter-element spacings of 4 cm, are employed. In this paper, we present a uniform linear multi-coil array-based borehole TEM system for the NDE of downhole casings, where the main contribution is the elimination of the influence of TRD; this paper does not focus on the interpretation algorithms. Therefore, the processed array outputs for different casing structures and different sampling times are analyzed instead of the interpreted thicknesses, because different interpretation methods using data even at different sampling times may have quite different results. Moreover, by analyzing the results produced by the proposed method with respect to the experimental casing structures, the improvement of array outputs could substantially correspond to a better NDE performance, even though the interpreted thickness changes are not illustrated.
From Figure 6 and Figure 7, we find that a serious distortion exists when the induced EMFs of the eight receivers are directly summed compared to the receiver with z = 0. Using the LCMV-based uniform linear multi-coil array for the borehole TEM system, the weighted array outputs in the last three cases perform more smoothly than in the first case, which reveals that a higher SNR can be achieved due to the accumulation of the coherent received signals. Notably, if the signals obtained from the first case were averaged, although there will be no phase shift and the averaging can increase the SNR substantially, repeated measurements are needed to achieve the same SNR as the proposed method, which will be quite a time-consuming procedure since the motion velocity should be slowed down compared to that of the multi-coil array-based borehole TEM system. It has also to be mentioned that since the two experimental casing structures have only one pipe per case, the amplitude of the late time data usually monotonically increases with the thickness of the metal pipe [22], and the shape of late time data will be more similar to that of the real casing thickness, as shown in Figure 3. Furthermore, comparing the experimental results for the two types of casing structures at 40 ms, the TRD will have a stronger influence on the NDE performance of Casing B than that of Casing A, which reveals that the TRD will also influence the longitudinal resolution of the NDE of downhole casings. As a result, the proposed method seems to have more effect on improving the NDE accuracy of the early time data (20 ms) and Casing A with small spacings than that of the late time data (40 ms) and Casing B with large spacings. Nevertheless, the proposed method is effective to achieve better performance for the experimental results of all cases, regardless of the casing structures and sampling times.
In addition, the weighted array outputs of the third and fourth cases with Δz = 2 cm have nearly the same performance as the first case, where the longitudinal shifting caused by the influence of the TRD can be effectively cancelled. Furthermore, comparing the last three cases, it is obvious that the inter-element spacing Δz and the number of receivers M influence the performance of the NDE. Because the signal strength is decreased by increasing the TRD, the SNR of each single receiver is inversely proportional to the inter-element spacing, where receivers with too large TRDs become useless. Conversely, in Section 4, it is shown that the number of receivers M is associated with the number of the Legendre polynomial stage P. Therefore, with a large value of M, in addition to the improved SNR due to the accumulation of the coherent signals, a large value of P can be achieved so that the numerical approach can approximate the borehole TEM model more closely. However, a larger number of receivers also requires a high level of accuracy in the fabrication of the uniform linear multi-coil array. Without a loss of generality, we employ the normalized root-mean-square error (RMSE) of the different inter-element spacings to further illustrate the effect of the number of receivers, where the RMSE is defined as:
RMSE = 1 S Q P s q p = 1 S Q P ( F s q p [ W T X ( z ) ] s q p ) 2 ,
and Fsqp and [WTX(z)]sqp denote the sqpth element of the vectors F and WTX(z), respectively. The simulation results of the RMSE are shown in Figure 8.
With increasing M, we find that the width of the null of the RMSE curves becomes much narrower. Taking M = 8 as an example, even though the neighboring receivers with small inter-element spacing can be weighted to obtain a better performance, the tolerance of the proposed method for the array geometric error will also be reduced so that a high accuracy array fabrication is required, where even a 1-mm error of the inter-element spacing in the fabrication may introduce a much larger error than that of M = 6 or M = 4. Therefore, a reasonable number of receivers should consider not only the SNR and the numerical approximation, but also the array geometric error to achieve a better NDE performance. The optimization of the receiver number and the inter-element spacing of the multi-coil array need to be investigated further in future studies.

6. Conclusions

A uniform linear multi-coil array-based borehole TEM system was proposed to improve the NDE performance for downhole casings. We presented a Gauss–Legendre quadrature-based numerical approach for the borehole TEM signal model. It was shown that the TEM response is coupled to the TRD, which greatly influences the NDE performance. Moreover, on the basis of the proposed numerical approximation, the received signals of the uniform linear multi-coil array were weighted according to the LCMV criterion to cancel the influence of the TRD, where the optimization and performance analyses of the proposed array weighting were also investigated. Simulations and experiments for a standardized oil-well casing inspection demonstrated the effectiveness of the proposed system.

Author Contributions

B.D., L.Y., and C.L. conceived and designed the experiments; Y.Z., H.L., and B.S. performed the experiments; R.D. and L.Y. analyzed the data; and B.D. wrote the paper.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 51504194), the Natural Science Basic Research Plan in Shaanxi Province of China (Grant No. 2018JQ5133), and the Natural Science Foundation of Shaanxi Provincial Department of Education (Grant No. 17JS106).

Acknowledgments

The authors would like to thank the reviewers for their helpful suggestions, which have considerably improved the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, C.; Liu, F.; Lin, J.; Zhu, K.G.; Wang, Y.Z. An Optimized Air-Core Coil Sensor with a Magnetic Flux Compensation Structure Suitable to the Helicopter TEM System. Sensors 2016, 16, 508. [Google Scholar] [CrossRef] [PubMed]
  2. Ezersky, M.G.; Frumkin, A. Evaluation and mapping of Dead Sea coastal aquifers salinity using Transient Electromagnetic (TEM) resistivity measurements. C. R. Geosci. 2017, 349, 1–11. [Google Scholar] [CrossRef]
  3. Danielsen, J.E.; Auken, E.; Jørgensen, F.; Jørgensen, F.; Søndergaard, V.; Sørensen, K.I. The application of the transient electromagnetic method in hydrogeophysical surveys. J. Appl. Geophys. 2003, 53, 181–198. [Google Scholar] [CrossRef]
  4. Dutta, S.M.; Reiderman, A.; Schoonover, L.G. New Borehole Transient Electromagnetic System for Reservoir Monitoring. Petrophysics 2012, 53, 222–232. [Google Scholar]
  5. Spies, B.R. Electrical and electromagnetic borehole measurements: A review. Surv. Geophys. 1996, 17, 517–556. [Google Scholar] [CrossRef] [Green Version]
  6. Huang, S.; Wang, S. The Pulsed Eddy Current Testing. In New Technologies in Electromagnetic Non-Destructive Testing; Springer: Singapore, 2016; pp. 41–80. [Google Scholar]
  7. Garcíamartín, J.; Gómezgil, J.; Vázquezsánchez, E. Non-Destructive Techniques Based on Eddy Current Testing. Sensors 2011, 11, 2525–2565. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Mao, X.F.; Lei, Y.Z. Thickness measurement of metal pipe using swept-frequency eddy current testing. Ndt E Int. 2016, 78, 10–19. [Google Scholar] [CrossRef]
  9. Brill, T.M.; Demichel, C.; Nichols, E.A.; Bermudez, F.Z. Electromagnetic Casing Inspection Tool for Corrosion Evaluation. In Proceedings of the International Petroleum Technology Conference, Bangkok, Thailand, 7–9 February 2012. [Google Scholar]
  10. Bateman, R.M. Casing Inspection Cased-Hole Log Analysis and Reservoir Performance Monitoring; Springer: New York, NY, USA, 2015; pp. 227–243. [Google Scholar]
  11. Park, D.G.; Kishore, M.B.; Kim, J.Y.; Jacobs, L.J.; Lee, D.H. Detection of Corrosion and Wall Thinning in Carbon Steel Pipe Covered With Insulation Using Pulsed Eddy Current. J. Magn. 2016, 21, 57–60. [Google Scholar] [CrossRef] [Green Version]
  12. Rifai, D.; Abdallam, A.N.; Razali, R.; Ali, K.; Faraj, M.A. An eddy current testing platform system for pipe defect inspection based on an optimized eddy current technique probe design. Sensors 2017, 17, 579. [Google Scholar] [CrossRef] [PubMed]
  13. Marinov, S.G. Theoretical and experimental investigation of eddy current inspection of pipes with arbitrary position of sensor coils. In Review of Progress in Quantitative Nondestructive Evaluation; Thompson, D.O., Chimenti, D.E., Eds.; Springer: Boston, MA, USA, 1986; Volume 5A, pp. 225–232. [Google Scholar]
  14. Wu, T.; Bowler, J.R.; Theodoulidis, T.P. Eddy-Current induction by a coil whose axis is perpendicular to that of a Tube. IEEE Transact. Magn. 2017, 53, 1–9. [Google Scholar] [CrossRef]
  15. Fu, Y.; Yu, R.; Peng, X.; Ren, S. Investigation of casing inspection through tubing with pulsed eddy current. NDT E Int. 2012, 27, 353–374. [Google Scholar] [CrossRef]
  16. Dang, B.; Yang, L.; Du, N.; Liu, C.; Dang, R.; Wang, B.; Xie, Y. Auxiliary sensor-based borehole transient electromagnetic system for the nondestructive inspection of multipipe strings. Sensors 2017, 17, 1836. [Google Scholar] [CrossRef] [PubMed]
  17. Marinov, S.G. Improved interpretation of the downhole casing inspection logs for two strings of pipes. In Review of Progress in Quantitative Nondestructive Evaluation; Thompson, D.O., Chimenti, D.E., Eds.; Springer: Boston, MA, USA, 1987; Volume 6A, pp. 1673–1679. [Google Scholar]
  18. Dashevsky, A.; Yu, A. Principles of Induction Logging; Elsevier: Amsterdam, The Netherlands, 2003. [Google Scholar]
  19. Bo, H.C.; Ji, H.K.; Cheon, J.P.; Rim, C.T. Synthesized magnetic field focusing using a current-controlled coil array. IEEE Magn. Lett. 2016, 7, 1–4. [Google Scholar] [CrossRef]
  20. Tian, X.; Liu, C.; Shen, L.C. A stable algorithm for simulation of array induction and measurement-while-drilling logging tools. Radio Sci. 2016, 33, 949–956. [Google Scholar] [CrossRef]
  21. Onegova, E.V. Effect of multicoil electromagnetic tool eccentricity on measured signals. Russian Geol. Geophys. 2010, 51, 423–427. [Google Scholar] [CrossRef]
  22. Dang, B.; Yang, L.; Dang, R.; Xie, Y. Borehole electromagnetic induction system with noise cancelation for casing inspection. IEICE Electron. Express 2016, 13, 20160714. [Google Scholar] [CrossRef]
  23. Nemani, A.; Eidiani, M. Linearly Constrained Minimum Variance Beamforming. In Time-Domain Beamforming and Blind Source Separation; Bourgeois, J., Minker, W., Eds.; Springer: New York, NY, USA, 2009; pp. 27–38. [Google Scholar]
  24. Li, J.; Farquharson, C.G.; Hu, X. Three effective inverse Laplace transform algorithms for computing time-domain electromagnetic responses. Geophysics 2016, 81, 75–90. [Google Scholar] [CrossRef]
  25. Desjardins, D.R.; Vallières, G.; Whalen, P.P.; Krause, T.W. Advances in transient (pulsed) eddy current for inspection of multi-layer aluminum structures in the presence of ferrous fasteners. Am. Inst. Phys. 2012, 1, 400–407. [Google Scholar] [CrossRef]
  26. Swarztrauber, P.N. On computing the points and weights for Gauss-Legendre Quadratur. Soc. Ind. Appl. Math. 2002, 24, 945–954. [Google Scholar] [CrossRef]
Figure 1. Linear multi-coil array-based borehole transient electromagnetic (TEM) system.
Figure 1. Linear multi-coil array-based borehole transient electromagnetic (TEM) system.
Sensors 18 02707 g001
Figure 2. The weighting structure of the linear multi-coil array output.
Figure 2. The weighting structure of the linear multi-coil array output.
Sensors 18 02707 g002
Figure 3. Experimental metal casing structures with different types of thickness changes: (a) three types of thickness changes with spacings of 20 cm (marked with casing A) and (b) four types of thickness changes with spacings of 10 cm (marked with casing B).
Figure 3. Experimental metal casing structures with different types of thickness changes: (a) three types of thickness changes with spacings of 20 cm (marked with casing A) and (b) four types of thickness changes with spacings of 10 cm (marked with casing B).
Sensors 18 02707 g003
Figure 4. Normalized induced electromagnetic forces (EMFs) of the eight receivers of the linear multi-coil array for casing A at (a) 20 ms and (b) 40 ms.
Figure 4. Normalized induced electromagnetic forces (EMFs) of the eight receivers of the linear multi-coil array for casing A at (a) 20 ms and (b) 40 ms.
Sensors 18 02707 g004
Figure 5. Normalized induced electromagnetic forces (EMFs) of the eight receivers of the linear multi-coil array for casing B at (a) 20 ms and (b) 40 ms.
Figure 5. Normalized induced electromagnetic forces (EMFs) of the eight receivers of the linear multi-coil array for casing B at (a) 20 ms and (b) 40 ms.
Sensors 18 02707 g005
Figure 6. Comparison of the normalized induced electromagnetic forces (EMFs) between the traditional receiver and the linear multi-coil array for casing A at (a) 20 ms and (b) 40 ms.
Figure 6. Comparison of the normalized induced electromagnetic forces (EMFs) between the traditional receiver and the linear multi-coil array for casing A at (a) 20 ms and (b) 40 ms.
Sensors 18 02707 g006
Figure 7. Comparison of the normalized induced electromagnetic forces (EMFs) between the traditional receiver and the linear multi-coil array for casing B at (a) 20 ms and (b) 40 ms.
Figure 7. Comparison of the normalized induced electromagnetic forces (EMFs) between the traditional receiver and the linear multi-coil array for casing B at (a) 20 ms and (b) 40 ms.
Sensors 18 02707 g007
Figure 8. Normalized root-mean-square error (RMSE) of different inter-element spacings.
Figure 8. Normalized root-mean-square error (RMSE) of different inter-element spacings.
Sensors 18 02707 g008
Table 1. Parameters of the multi-coil array sensor and the field experiment.
Table 1. Parameters of the multi-coil array sensor and the field experiment.
ParameterSymbolValue
Radius of the multi-coil array sensor r112 mm
Number of receiving coilsM8
Inter-element spacingΔz20 mm
Transmitting–receiving distancesz1–zM20–160 mm
Number of transmitting coil turnsNT19
Number of receiving coil turnsNR62
Tool housing inner radiusr218.5 mm
Tool housing outer radiusr321.5 mm
Standardized casing inner radiusr462.13 mm
Standardized casing outer radiusr569.85 mm
Cement ring outer radiusr688.9 mm

Share and Cite

MDPI and ACS Style

Dang, B.; Yang, L.; Liu, C.; Zheng, Y.; Li, H.; Dang, R.; Sun, B. A Uniform Linear Multi-Coil Array-Based Borehole Transient Electromagnetic System for Non-Destructive Evaluations of Downhole Casings. Sensors 2018, 18, 2707. https://doi.org/10.3390/s18082707

AMA Style

Dang B, Yang L, Liu C, Zheng Y, Li H, Dang R, Sun B. A Uniform Linear Multi-Coil Array-Based Borehole Transient Electromagnetic System for Non-Destructive Evaluations of Downhole Casings. Sensors. 2018; 18(8):2707. https://doi.org/10.3390/s18082707

Chicago/Turabian Style

Dang, Bo, Ling Yang, Changzan Liu, Yahong Zheng, Hui Li, Ruirong Dang, and Baoquan Sun. 2018. "A Uniform Linear Multi-Coil Array-Based Borehole Transient Electromagnetic System for Non-Destructive Evaluations of Downhole Casings" Sensors 18, no. 8: 2707. https://doi.org/10.3390/s18082707

APA Style

Dang, B., Yang, L., Liu, C., Zheng, Y., Li, H., Dang, R., & Sun, B. (2018). A Uniform Linear Multi-Coil Array-Based Borehole Transient Electromagnetic System for Non-Destructive Evaluations of Downhole Casings. Sensors, 18(8), 2707. https://doi.org/10.3390/s18082707

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop