# Research on Application of Wax Deposition Detection in the Nonmetallic Pipeline Based on Electrical Capacitance Tomography

Authors: Nan Li; Kui Liu; Xiangdong Yang; Mingchen Cao
Journal: Journal of Sensors
 (2016)

---

## Abstract

Wax deposition detection in nonmetallic pipelines is an important requirement in the oil industry. In this paper, an ECT (electrical capacitance tomography) sensor is developed for wax deposition detection in nonmetallic pipelines. Four wax models with different concentrations were established for detection. These models were analyzed through simulations and practical experiments simultaneously and data were compared. A linear back projection algorithm is applied to reconstruct the image with both simulated and experimental data. A comparison of binary images with different concentration of stratified flow was demonstrated; this illustrates that the difference in concentration between the experimental results and profile distribution is less than 1.2%. The experimental results indicate that the ECT system is valid and feasible for detecting the degree of wax deposition in the nonmetallic pipelines.


---

## Body

## 1. Introduction

Heat and pressure are widely used in oil and gas pipeline transmission to reduce the viscosity of crude oil and the deposition of wax. However, due to difference in radial temperature, radiation loss occurs rapidly. The wax is gradually separated out and deposited on the pipe wall, which can lead to condensate tube accidents [1]. The situation is more serious in China, as the crude oil produced in China has several distinguishing characteristics, such as high waxy distillate and low temperature viscosity. Therefore, it is important and necessary to develop a technique that can detect the deposition degree of wax inside the pipe. With consideration to the operating environment of the pipeline, as well as the cost of detection, for practical application, nondestructive-based techniques are highly recommended. Many nondestructive techniques have achieved some degree of success in the pipeline detection field, for instance, the eddy current testing method [2, 3], ultrasonic-based technique [4–10], radiological image method [11], and magnetic flux leakage detection method [12–14]. However, most of these techniques are applied to detect pipeline defects, and it is difficult to detect the sediment (i.e., paraffin wax) of the pipeline by using the above technologies, due to various reasons.The detection of wax deposition in the pipeline is essentially a multiphase detection, which is composed of the wax layer, oil layer, and gas. ECT (electrical capacitance tomography) is a visual procession from measurement capacitance to object distribution. The method is a noninvasive and nondestructive technique. An ECT sensor is a typical “soft” field sensor, which means that the Electromagnetic (EM) field propagates across the entire probed volume, as shown in Figure1. The measurement at the volume surface depends on the values of the measured quantity everywhere in the volume.Figure 1
Schematic diagram of the “soft” field.The related research began in the 1980s at the University of Manchester and, from hardware, image construction algorithms to sensor design has made considerable progress [15–21]. Johansen et al. developed a dual sensor system consisting of an 8-electrode ECT sensor and γ-ray sensor. The new system was applied to reconstruct the image of oil-gas-water three-phase flow [22]. Dyakowski et al., based on this sensor structure, presented an image reconstruction algorithm, which can be used to evaluate the concentration of a measured cross-section [23]. Schmitz and Mewes applied the ECT to the chemical industry, and the system was used to detect the transient multiphase flow in bubble columns [24]. In 2010, Yang introduced an application that used an ECT sensor to analyze the drying, pellets, and coating process of pharmaceutical fluidized beds [25]. Ismail et al. applied the ECT to the oil and gas separation process in the oil industry to determine the concentration and flow rate of each stage [26]. Al Hosani et al. suggested to use a narrowband pass filter to enhance the resolution of the produced images; the experimental results indicated that the method can improve the resolution of reconstruction images [27].Sensor optimization is an important part of the ECT system. Xie et al. presented a uniformity of sensitivity distribution as an evaluation criterion to determine the performance of the ECT sensor [28]. Peng et al. analyzed the influence of electrode length on sensitivity distribution of the ECT sensor. The results showed that the electrode length of the ECT sensor should be at least equal to the width of the sensor electrode [29]. In 2012, Peng et al. discussed the influence of the number of the ECT sensor electrodes on the quality of image reconstruction and recommended that the number of electrodes be 12 [30]. Ren et al. presented a new algorithm to reconstruct the three-dimensional inclusion shapes; the method is valid and effective [31].The purpose of the present paper is to develop an ECT system for wax deposition detection of a nonmetallic pipeline. The optimal sensor structure is determined based on our previous research [32]. An ECT system is demonstrated in the experimental part, and the feasibility of the system is verified based on a comparison of the reconstructed images according to simulation and experimental results.


## 2. Methodology

### 2.1. Wax Models

In the long-range transportation process, the temperature of the crude oil decreases and wax deposition occurs due to the thermal radiation. The permittivity of wax is 1.9~2.5. The wax is deposited on the inner wall of the pipeline, and the structure could be simplified as a stratified flow pattern in a laboratory environment. A series of wax deposition models is built and the cross-sections are shown in Figure2. Note that the last model of Figure 2 is a special model, namely, heteromorphous flow. In order to further clarify the different structures of the wax, u is defined as concentration. Note that (1) is valid for 2D problems:(1)
u
=
A
w
a
x
A
0
×
100
%
,where A
w
a
x is the area of wax and A
0 refers to the area of sensing area, that is, the cross-section of the pipeline.Figure 2
Different thicknesses of wax deposition models. (a) Concentration is 11%. (b) Concentration is 18.5%. (c) Concentration is 52.2%. (d) Concentration is 62.5% (note: the last one is special stratified flow–heteromorphous flow).
(a)
(b)
(c)
(d)


### 2.2. Setup of Measurement System

The ECT sensor electrodes consist of copper, and earthed shielding plates are placed around the sensor electrodes. The measurement system and the sensor structure are shown in Figure3.Figure 3
ECT system and sensor structure. (a) Diagram of measurement system and practical ECT sensor. (b) Parameters of the ECT sensor.
(a)
(b)The shielding wire and BNC joints are used to prevent interference. In this study, the pipe consists of PVC material. The end guarding plate consists of copper, and the earthed shielding plate consists of aluminum. According to our previous research [32], the parameters of the ECT sensor are determined as listed in Table 1: N is the number of electrodes; W refers to the central angle of the electrodes; and r
1, r
2, and r
3 are the radii of the pipe’s interior and exterior and the radius of the earthed shielding plate, respectively. L
p
l
a
t
e represents the length of the sensor electrodes. The width of the end guarding plate is L
e
n
d, and the distance between the end guarding plate and sensor electrode is defined as L
g
a
p.Table 1
Parameters of the sensor.
Sensor electrodes
Pipeline
Earthed shielding plate
End guarding plate
Materials
Copper
PVC
Aluminum
Copper
Parameters
N
W
L
p
l
a
t
e
r
1
r
2
r
3
L
e
n
d
L
g
a
p
6
15°
105 mm
35 mm
40 mm
50 mm
20 mm
10 mm


## 2.1. Wax Models

In the long-range transportation process, the temperature of the crude oil decreases and wax deposition occurs due to the thermal radiation. The permittivity of wax is 1.9~2.5. The wax is deposited on the inner wall of the pipeline, and the structure could be simplified as a stratified flow pattern in a laboratory environment. A series of wax deposition models is built and the cross-sections are shown in Figure2. Note that the last model of Figure 2 is a special model, namely, heteromorphous flow. In order to further clarify the different structures of the wax, u is defined as concentration. Note that (1) is valid for 2D problems:(1)
u
=
A
w
a
x
A
0
×
100
%
,where A
w
a
x is the area of wax and A
0 refers to the area of sensing area, that is, the cross-section of the pipeline.Figure 2
Different thicknesses of wax deposition models. (a) Concentration is 11%. (b) Concentration is 18.5%. (c) Concentration is 52.2%. (d) Concentration is 62.5% (note: the last one is special stratified flow–heteromorphous flow).
(a)
(b)
(c)
(d)


## 2.2. Setup of Measurement System

The ECT sensor electrodes consist of copper, and earthed shielding plates are placed around the sensor electrodes. The measurement system and the sensor structure are shown in Figure3.Figure 3
ECT system and sensor structure. (a) Diagram of measurement system and practical ECT sensor. (b) Parameters of the ECT sensor.
(a)
(b)The shielding wire and BNC joints are used to prevent interference. In this study, the pipe consists of PVC material. The end guarding plate consists of copper, and the earthed shielding plate consists of aluminum. According to our previous research [32], the parameters of the ECT sensor are determined as listed in Table 1: N is the number of electrodes; W refers to the central angle of the electrodes; and r
1, r
2, and r
3 are the radii of the pipe’s interior and exterior and the radius of the earthed shielding plate, respectively. L
p
l
a
t
e represents the length of the sensor electrodes. The width of the end guarding plate is L
e
n
d, and the distance between the end guarding plate and sensor electrode is defined as L
g
a
p.Table 1
Parameters of the sensor.
Sensor electrodes
Pipeline
Earthed shielding plate
End guarding plate
Materials
Copper
PVC
Aluminum
Copper
Parameters
N
W
L
p
l
a
t
e
r
1
r
2
r
3
L
e
n
d
L
g
a
p
6
15°
105 mm
35 mm
40 mm
50 mm
20 mm
10 mm


## 3. Results and Discussion

The practical stratified models are shown in Figure4, and the concentrationsu of the deposited wax are 11%, 18.5%, 52.2%, and 62.5%, respectively. In order to calibrate the measurement results, two more situations should be considered, namely, those for an empty pipe and a pipe full of wax. Moreover, simulation and practical experiments are processed and analyzed simultaneously, presented, respectively, as C
s
i
m and C
m. The results are listed in Table 2. The simulation model is 3D model; FEM (finite element method) is used to analyze the sensitivity distribution. The conditions of simulation models are the same as practical models. However, because the electrostatic module is applied, the frequency of excitation is not considered during simulation.Table 2
Comparison of simulation results with experimental results.
Flow patterns
Empty[fF]
Full wax[fF]
Stratified flow (a)[fF]
Stratified flow (b)[fF]
Stratified flow (c)[fF]
Stratified flow (d)[fF]
u
0%
100%
11%
18.5%
52.2%
62.5%
C
i
j
C
s
i
m
C
m
C
s
i
m
C
m
C
s
i
m
C
m
C
s
i
m
C
m
C
s
i
m
C
m
C
s
i
m
C
m
1
2
300.9
288.3
419.8
382.3
380.0
335.6
427.6
350.5
448.9
383.0
444.8
385.3
1
3
45.9
47.7
100.9
100.7
56.0
53.2
58.7
56.8
72.1
71.2
70.2
71.4
1
4
30.8
30.5
70.3
69.3
37.2
35.8
39.0
37.1
42.7
42.7
45.8
50.4
1
5
45.3
44.7
99.7
85.3
55.1
44.7
58.4
47.6
71.4
60.8
96.4
90.5
1
6
298.0
264.1
416.5
351.4
376.0
317.7
430.4
335.0
445.3
358.0
432.0
346.6
2
3
304.4
307.5
424.4
405.8
300.6
301.1
300.3
297.0
343.7
331.1
340.3
338.3
2
4
45.2
48.7
99.4
100.8
43.1
47.6
44.7
48.1
56.8
59.2
57.0
66.1
2
5
30.9
29.1
70.7
68.7
29.0
29.1
31.0
32.4
45.8
44.2
62.9
72.4
2
6
45.1
46.6
99.3
99.1
48.3
60.2
68.0
77.3
116.0
110.5
111.9
104.3
3
4
298.1
297.0
416.4
397.5
296.8
294.5
295.6
295.6
280.2
285.3
281.3
274.8
3
5
45.3
44.3
99.5
96.6
43.9
43.2
42.5
43.5
38.3
38.4
47.9
62.4
3
6
30.7
30.4
70.2
67.7
29.0
30.7
31.1
33.0
45.8
47.0
47.1
45.7
4
5
300.4
295.9
419.6
390.3
299.6
296.7
298.7
284.9
288.1
285.9
291.7
334.8
4
6
45.2
41.5
99.5
90.4
43.6
41.0
45.2
42.7
57.9
54.4
69.9
72.1
5
6
298.6
171.5
417.2
268.7
296.9
168.4
298.7
173.2
340.3
209.1
397.1
273.4Figure 4
Practical models of stratified flow wax. (a) Concentration is 11%. (b) Concentration is 18.5%. (c) Concentration is 52.2%. (d) Concentration is 62.5%.
(a)
(b)
(c)
(d)In order to illustrate the comparison results more intuitively, the normalized capacitances are calculated according to (2), and the results are shown in Figure 5. Consider (2)
n
o
r
C
=
C
m
-
C
l
C
h
-
C
l
,wherenorC is the normalized value and C
l and C
h are measured capacitance filled with low dielectric constant and high permittivity materials, respectively. For example, in this study C
l is the capacitance when the pipe is empty, and C
h is the capacitance when the pipe is full of wax.Figure 5
Normalized capacitance trend analysis. (a) Concentration is 11%. (b) Concentration is 18.5%. (c) Concentration is 52.2%. (d) Concentration is 62.5%.
(a)
(b)
(c)
(d)According to Figure5, the following can be derived. The same tendencies for the measurement results from the simulation and experiments are observed. For regular stratified flow patterns, the difference of the two curves decreases with the increase of concentration u, as shown in Figures 5(a)–5(c). The reason for this comes from two aspects: (i) the sensitivity distribution of the ECT sensor which was nonuniform; the sensitivity close to the sensor electrodes was much higher than that in the center of the sensing area; (ii) the inaccurate forward modeling or the noise in measurements. Note that the difference of the two curves increases in Figure 5(d) compared with Figure 5(c), the reason for which is that a heteromorphous model affects the measurement results.In order to examine the effects of the design factors, it was necessary to perform the image reconstructions. Large numbers of algorithms have been developed to derive reconstructed images over the past several decades, such as linear back projection (LBP), Landweber iteration, and the sparsity-inspired image reconstruction method. In [19], Yang and Peng have analyzed the characterizations of five image reconstruction algorithms and three criteria were used to evaluate the performance of reconstructed images:(1)
relative image error,
(2)
relative capacitance residual,
(3)
correlation coefficient between the test object and the reconstruction.Consider(3)
I
m
a
g
e
e
r
r
o
r
=
g
^
-
g
g
,
C
a
p
a
c
i
t
a
n
c
e
r
e
s
i
d
u
a
l
=
λ
-
S
·
g
^
λ
,
C
o
r
r
e
l
a
t
i
o
n
c
o
e
f
f
i
c
i
e
n
t
=
∑
i
=
1
N
g
^
i
-
g
^
-
g
i
-
g
-
∑
i
=
1
N
g
^
i
-
g
^
-
2
∑
i
=
1
N
g
i
-
g
-
2
,where g is the true permittivity distribution of the test object, g
^ is the reconstructed permittivity distribution, and g
- and g
^
- are the mean values of g and g
^, respectively. λ is the normalized capacitance vector, S represents the normalized sensitivity matrix, and N refers to the number of pixels. According to [19], the image error and correlation coefficient of LBP are better than SVD and Tikhonov and worse than Iterative Tikhonov method. The capacitance residual of LBP is similar to Iterative Tikhonov and Projected Landweber, but worse than SVD and Tikhonov methods. Due to the fact that noniterative algorithms are typically used in practical applications and that the time consumption is less than in iterative algorithms, in order to simplify the comparison, the LBP algorithm is used to reconstruct the image [32]. Moreover, the LBP method is more suitable than iterative algorithms for dealing with stratified flow patterns. The reconstructed images are shown in Figure 6. Note that the images are based on the normalized data of the simulation and experiment data, separately. Figure 6(a) is a series of reconstructed images based on the simulation results, and Figure 6(b) is a series of reconstructed images based on the experimental results.Figure 6
Different concentrations of stratified flow image results. (a) Simulation results. (b) Experimental results.
(a)
(b)The imaging results of the experiments are consistent with the reconstructed images of the simulations, which signifies that the designed ECT system is valid and feasible. For different concentrations of the wax, the reconstructed images are varied. Especially for the heteromorphous flow pattern, the shape of the semicircular convex is difficult to identify. The boundary of the heteromorphous flow pattern is difficult to determine, and the stratified flow resembles an inclined plane. According to [19, 21], this issue could be solved in terms of two aspects, namely, algorithm improvement and ECT sensor optimization. In order to evaluate the validation of using reconstructed images to identify the concentration of the wax, threshold processing was completed, and binary images of different concentrations of stratified flow are presented and compared with the standard models. The threshold processing is based on adaptive threshold filtering method which was presented by Xie et al. [33] as shown in Figure 7.Figure 7
Flowchart of the linear back projection threshold processing method.In Figure7, β
m is the initial concentration, which can be determined via a nonthreshold imaging results and experience. In this study, we define β
m
=
N
(
G
s
(
e
k
)
m
a
x
)
/
N
(
G
s
(
e
k
)
), where N is the number of pixels, G
s represents standard binary image, N
(
G
s
(
e
k
)
m
a
x
) refers to the number of pixels whose value is 1 in standard binary image, and N
(
G
s
(
e
k
)
) represents the total number of pixels in standard binary image. The concentration increment Δ
β
=
|
β
t
-
β
m
| and the initial value is set to 0, where β
t
=
N
(
G
n
e
w
(
e
k
)
m
a
x
)
/
N
(
G
n
e
w
(
e
k
)
); N
(
G
n
e
w
(
e
k
)
m
a
x
) represents number of pixels whose value is 1 in processed image and N
G
n
e
w
e
k represents the total number of pixels in processed image. T
r refers to threshold, and Δ
T is the threshold increment whose initial value is set to 0. Taking the flow chart as an example, each iteration step increment is 0.0005; steps can be corrected according to the actual operation. Threshold is determined by the following equation:(4)
T
r
=
G
m
i
n
+
Δ
T
G
m
a
x
-
G
m
i
n
.The threshold value takes the minimum pixel value as a reference; the increment is the product of pixel value difference and threshold increment. In this paper, the threshold values are 0.705, 0.78, 0.38, and 0.395, respectively, for the four reconstructed images.The binary images of different concentrations of stratified flows are shown in Figure8, and the concentrations and correlation coefficients are listed in Table 3. The difference of concentrations between the standard models and reconstructed images is less than 1.2%, which indicates the validation of the ECT system. However, the definition of concentration is the ratio of wax deposition to the total sensing cross-section. It may be improper to evaluate the quality of reconstructed image. Therefore, the correlation coefficients of different models are listed in Table 3.Table 3
Contrast between standard images and reconstructed images.
Model 1
Model 2
Model 3
Model 4
Concentration of standard images
11%
18.5%
52.5%
62.5%
Concentration of reconstructed images
12.2%
19.2%
52.97%
62.9%
Correlation coefficient
93.79%
92.78%
98.01%
87.04%Figure 8
Binary images of different concentrations of stratified flow. (a) Profile distribution. (b) Threshold processing image based on experimental data.
(a)
(b)As we can see from Table3, when concentration is 52.5% (model 3), the correlation coefficient is the best (98.01%). When the concentration is 11% and 18.5% (model 1 and model 2), the correlation coefficients are similar to each other and the quality of images is acceptable. However, when the concentration is 62.5% (model 4), the correlation coefficient is the worst (87.04%). The main reason is the limitation of the LBP algorithm which is not good at imaging of heteromorphous flow, especially for edge imaging.


## 4. Conclusions

For wax deposition detection in nonmetallic pipelines, the use of an ECT sensor was demonstrated. LBP algorithm was applied to reconstruct images. Four models were tested, and the reconstructed images showed that the sensor is feasible. However, due to the limitations of the algorithm, it was difficult to identify the special stratified model. The boundary of the heteromorphous flow pattern was recognized as an incline in the reconstructed images. For the concentrations of 11%, 18.5%, 52.5%, and 62.5%, the differences between the profile and reconstructed image were 1.2%, 0.7%, 0.47%, and 0.4%, respectively. The correlation coefficients between standard images and reconstructed images were 93.79%, 92.78%, 98.01%, and 87.04%, separately. This proves that the system designed in this paper can be applied to wax deposition detection of nonmetallic pipelines.



---

*Source: 7390470-2016-08-07.xml*