UDC 621.396.41
E.I. Trenkal, A.G. Loshchilov
New approach for increasing the precision of TDR analysis of multilayer environments
The authors propose the solution allowing to determine the structure and the layers parameters of the multiphase liquid in real-time mode and thereby to increase the measurement precision. The method is based on using the periodic controllable loads in the measuring probe construction, which can change their own parameters at the required time moment. In the paper a description of the mathematic model of the proposed solution is presented; an algorithm for determination of the layers parameters of the multiphase liquid is proposed; the calculation example of the model multiphase liquid parameters using the proposed algorithm is shown. Also the comparison of the precision of the proposed and standard approaches is presented.
Keywords: level measurement, multiphase liquid, TDR, method, measuring probe, controllable loads, parameters measurement.
doi: 10.21293/1818-0442-2016-19-4-5-9
One of the standard processes of modern industry is measurement of liquid levels in tanks. In spite of the availability of a variety of methods and techniques for implementation of this procedure, their application efficiency and measurement precision are often related to the parameters of measuring liquids and their change under the influence of external factors. One of the actual tasks is the analysis of a multilayer medium, which is the combination of several liquids with different properties. An example is the monitoring task of layer-by-layer composition of liquids at the stage of oil-product clearing requiring the determination of the following levels: the oil, the bottom water and the emulsion layer.
One of the promising methods for solving this task is the time domain reflectometry (TDR) method, which was initially applied for the diagnostics of cable lines [1]. The method is as follows: the investigated object is affected by a video pulse or step function. Analyzing the object response to the test signal, its properties are determined.
In the 1960s the TDR method was first applied to determine the level of the single liquid [2]. In [3] the application of the method for measuring of the multiphase liquid's levels (fuel oil and seawater) is described. The article [4] shows the detailed description of physical principles of the liquid level measurement using the TDR method. In addition, in this article the comparison with the alternative methods is presented. Numerous papers [5-8] describing the practical implementation of the TDR method for liquids' level measurement are known.
The restriction for the existing solutions is that priori information about the parameters of the multiphase liquid for the accurate determination of its levels is used. Consideration of these parameters is performed either on the basis of reference data or as a result of the calibration procedure during the mounting of the measurement system. In the latter case the highest accuracy is achieved.
However, in practice the liquids' properties are not stationary and depend on external conditions and can be changed during the measurement. The interface between
adjacent layers is often an emulsion, whose parameters are quite difficult to predict. A high temperature of the measured liquids leads to intense evaporation of the upper layers and the change of the air layer properties. This is rarely taken into account when the method is used.
Problem definition
It may be concluded that the values of the predefined parameters of a liquid may be true only in a limited time interval. To improve the accuracy of the TDR analysis of multiphase liquids, it is necessary to dynamically determine the parameters of each layer during the measurement.
Thus, a relevant task is to find a method for determining the multiphase liquid's parameters. The method must allow getting the properties of the inner layers of this liquid in real time.
Theory
One of the ways to improve the accuracy of level sensors based on the TDR is shown in the patents [9, 10]. The authors offer to use the periodical discontinuities in the probe design for the calibration of signal propagation speed in the investigated environments using the delays of the signals reflected from these discontinuities. The disadvantage of such solution is the attenuation of the energy of the useful signal on this discontinuities and complication of the reflectogram due to the presence of a lot of reflections between them.
Based on this calibration method, we offer the following solution of the problem: to use the controllable loads with the known parameters in the probe design. It allows to realize two operating modes of the measuring probe: the calibration mode in which the condition of the controllable loads is successively changed and the parameters of the investigated environments are determined, and the measurement mode in which the controllable loads have no effect on the signal propagation in the probe, or distortion caused by their presence is not significant.
1. Mathematical model of the proposed solution
In Fig. 1 an example of the measuring probe with the controllable loads is shown.
For simplicity the case in which each layer of the measuring liquid contains at least one controllable load in the measuring probe will be considered. The equivalent circuit diagram of such measuring probe represents the cascading segments of the transmission lines Lj and the controllable loads, which represent the lumped obstacle Zk.
controllable load
1
( 1=5-
\ir
1=5-
1=5-
s Oil
1=5
1=5- Water
Fig.1. Example of the measuring probe with controllable loads
The model of such probe with the controllable loads can be described by the method of classic transmission matrix. The full transmission matrix of the measuring probe AP is equal to the product of the matrices AL, each of which complies with the segment of the measuring probe in a layer of the modeled liquid:
N
AP=H ALi , 1=1
where i - the index of the layer; N - the number of the modeling layers.
The equation of the transmission matrix of the measuring probe segment ALi can be written as:
M
ALj = (n4,j • AZj,j) • AjMi +i, j=i
where j - the index of the controllable load in i-th layer; Mi - the number of the controllable loads in i-th layer; Aj - the transmission matrix of the measuring probe segment; AZ,j - the transmission matrix of the controllable load.
The transmission matrix of the measuring probe segments L j can be written as:
4 j =
ch(y, ■ lij) Pi ■ sh(y, ■ ltj) — ■ li,j) ch(y, ■ li,j)
Pi
where l, - the length of the segment Lj pi - the impedance of the transmission line segment; yi> - the propagation coefficient equal to:
Y = a+j'P.
In the simplest case, the controllable load can be written as transmission matrix of a quadripole:
AZk =
1 0
"_V Zpk i
where Zpk - the impedance of the controllable load Zk.
Voltage on the input of the measuring probe immersed in a multiphase liquid is determined as
i r, , ~ \
Uin (t ) = F
-1
E ( f )
Ztn ( f ) Zn ( f ) + Zg
where E(f) - the spectrum of the test signal; Zg - the output resistance of the generator; F- inverse Fourier transform operator; Zin(f) - the frequency response of complex input resistance of the probe.
The frequency response of complex input resistance of the probe can be define by formula Z f = APo,o( f) • Zi + APo,i( f) APi,o( f) • Zi + APi,i( f), where Zl - the load resistance on the measuring probe output.
The main function of the controllable loads Zk is to perform the mode of full or partial reflection of the test signal for determination of the layer parameters of the liquid under investigation, which are calculated based
on the known distance i* from the measuring probe entrance to the controllable load.
2. Algorithm of determining the structure and the parameters of a multiphase liquid
The algorithm of determining the layer parameters of a multiphase liquid contains the following steps:
1. Formation of the vector of reflectograms.
The sequence of the calculation of the measuring probe reflectograms for the definite set of the controllable loads states is conducted. The reflectogram vector Uk is formed, where k - the number of a controllable load providing the mode of the full or partial reflection of the test signal. Index k = 0 corresponds to the mode in which all controllable loads have off-state.
2. Definition of the response array.
Based on the received reflectogram, the two arrays are formed: R - the array containing the information about the responses reflected from the phase boundaries and the end of the measuring probe. R* - the array containing the information about the responses reflected from the controllable loads.
The elements of the vectors R<i>, R*<1> are response delays relative to the test signal, and the elements of the vectors R<2>, R<2> are the amplitudes of the responses.
During formation of the array R, to avoid the parasitic responses caused by signals multiple reflections from the phase boundaries, the signal selection algorithm is applied [ii].
The value of the element R0 complies with test signal response.
To form the array R*, differences of reflectograms AU(/) determined in the step i are calculated:
AUk (t) = Uk (t) - Uo(t), where k e[i, K]; K - the number of the controllable loads in the measuring probe.
ffomadbi TYCYPa, № 1 (39), Mapm 2016
The values of the elements R* comply with delays and amplitudes of the first responses on each difference reflectograms AU(t)k.
3. Definition of the controllable loads belonging. Belonging of the controllable loads to the analyzed
liquid layers are defined. For that, element-by-element comparison of the delay vectors R*<1> u R<:> is provided.
For each element R*, the appropriate layer index i, in which the controllable load is located and the number of the controllable load j in the i-th layer are provided. As a result, each element of the array R* and the vector l* may be characterized by indexes i and j, writing them in
the form of R* j and l* j .
4. Definition of the layers parameters of the multiphase liquid.
The layers' parameters of the measured liquid and ai are determined for each layer in the following sequence:
a) based on the parameter of the responses reflected from the controllable loads, the average value of the propagation velocity u,- in i-th layer is calculated:
2 • (l* y - I ln)
=-!
=1
P*<1> R<1> Ri, j — Ri-1
(1)
b) the dielectric constant and the length of the layer are determined:
\2
(2)
li =-
e, =1 —
v U
•(R,<1> - R<1> ) 2
(3)
where c - the speed of light.
c) the average value of the attenuation coefficient ai [N] in i-th layer is calculated:
ln
=-!
R<2> R0
ff (
«=1
— r« ) • e-21« a"
R
!<2>
2 • (l,
i—1
i, j—! ln) n=1
(4)
where r - reflection coefficient from i-th layer determined by the formula:
3<2>
-. (5)
t, =-
R-
R0
<2>3—1((1
n=1
— r2)^ e-2 •ln
—2 • lj • ai
Modeling
The developed algorithm based on using the periodic controllable loads in the measuring probe construction has been tested on the mathematic model. The parameters of the model liquid are shown in Table 1.
The model of the measuring probe contained three controllable loads (K = 3), which are located in the different layers of the multiphase liquid and provide the
mode of full reflection in the on-state. The distances from the measuring probe entrance to the controllable loads are the following: l* = 0,375 m; l2 = 0,75 m; l3*= 1,125 m.
Table 1
Multiphase liquid parameters
Layer Dielectric Losses a, Layer length,
number constant e dB/m l, m
1 1 0,5 0,5
2 3 1 0,5
3 80 2 0,5
To determine the layer parameters of the multiphase liquid, the described algorithm is applied. At the first stage the reflectograms Uk(t) in the modes k = 0 (Fig. 2) and k = 1, 2, 3 (Fig. 3) are calculated and the response arrays R and R* are formed
R =
- t0 A ' " 0 1
"1 =t1 -10 A 3,33-10-9 —0,25
"2 = t2 —10 A 9,M0-9 —0,53
"3 = t3 —10 A3 _ 3,89 40-8 —0,34
Uo(t)
ei, ai, l1 e2, a2, l2 e3, a3, l3
t1
Fig. 2. Forming the reflectogram U0(t) in the mode k = 0
R =
* * * *
"1 = f[ — to A
* * * *
"2 = t2 —10 A2
* * * *
"3 = t3 —10 A3
2,5 T0-9 —0,96
6,22 T0-9 —0,83
1,66 40-8 —0,4
After forming the arrays R and R*, belonging of the controllable loads to the analyzed liquid layers is defined by comparison of the vectors elements R*<1> and R<*>:
Z*: d* <d* - belongs to the first layer, i = 1, j = 1;
Z2: d1 < d* < d2 - belongs to the second layer, i = 2, j = 1;
Z3: d2 < d3 < d3 - belongs to the third layer, i = 3, j = 1.
In Table 2 the parameters A*, d*, l* corresponding to the controllable loads are represented in the form of a table.
Finally, the layer parameters of the multiphase liquid are defined in sequence by formulas (1)-(5). Calculation results are shown in Table 3.
n
U
t
2
U (t)
Si, Œi, /1 S2, ^2, /2 S3, «3, /3
03
• 1 1 1 1 1 1 1 1 1 1 1
U (t)
U3 (t)
Table 4
t
r,
Fig. 3. Forming the reflectograms Ui(t) - U3(t), where a - the reflectogram in the mode k = 1 (Zi = 0); b - the reflectogram in the mode k = 2 (Z2 = 0); c - the reflectogram in the mode k = 3 (Z3 = 0).
Table 2
i J R j Distance to the controllable * load j , m
Response delay dij, s Response amplitude Ay, V
1 1 2,5-10-9 -0,96 0,375
2 1 6,22-10-9 -0,83 0,75
3 1 1,6610-8 -0,4 1,125
Table 3
layer numbei Layer length /i, m Losses ai, dB/m Impedance pi, Ohm Dielectric constant s Propagation velocity u, m/s
1 0,5 0,5 50 1 3-108
2 0,5 1 28,86 3 1,73-108
3 0,5 2 5,59 80 3,36-10'
Comparison of the standard approach and the new solution
For the evaluation of efficiency of the new approach it was performed the modeling of the measurement of the multiphase liquid parameters using the standard and new measurement approaches. The key difference between the approaches is method for evaluation of the environment dielectric constant: in the standard approach the reference data is used; in the proposed solution the measurement of the environment parameters is performed in calibration mode.
The parameters of the model multiphase environment and reference data about the dielectric constant of the layers are presented in the Table 4.
Layer number Dielectric constant S Reference dielectric constant s Losses a, dB/m Layer length /, m
1 1,05 1 0,5 0,3
2 2,3 2,1 1 0,5
3 70 80 2 0,5
The measuring probe which is realized the new approach contained six controllable loads (K = 6). The positions of the loads in the probe are shown in the Table 5. The probe length is equal to height of the tank.
Table 5
Controllable load number The distance to the controllable load
1 0,186
2 0,371
3 0,557
4 0,743
5 0,929
6 1,114
In the modeling the noise of generator and receiver are taken into consideration. The noises parameters are following: normal law of distribution, RMS-noise of the generator - 1 mV, RMS-noise of the receiver - 600 ^V. Also in the model the quantization procedure with resolution 16 bit was added.
The calculation result of the length of the layers for both approaches is presented in the Table 6. Also in the table the relative measurement errors are shown.
Table 6
Calculated length of miilti| phase liquid layers
Layer number Standard approach Proposed approach
Layer length /, m Relative error 5, % Layer length /, m Relative error 5, %
1 0,307 2,333 0,299 0,045
2 0,523 4,6 0,5003 0,062
3 0,468 6,4 0,499 0,102
The presented calculation results show that the proposed method has higher accuracy compared to standard method. The higher accuracy of the proposed solution pre-conditioned by the ability to measure of the parameters of investigated environment during the measurement procedure.
Conclusion
To determine the levels of multiphase liquids using TDR method with high precision, it is necessary to know the structure of the multiphase liquid and parameters of its layers. Taking into consideration the fact that the parameters of the measured liquid are often not stationary and depend on the environment condition, the procedure of parameters determination, at the best case, must be provided at the moment of the levels measurement.
In the paper, we propose a solution allowing to determine the structure and the layers parameters of the multiphase liquid in real-time mode and thereby to in-
a
t
b
t
t
c
t
2
t
*
fiomadbi TYCYPa, № 1 (39), Mapm 2016
crease the measurement precision. The method is based on using the periodic controllable loads in the measuring probe construction, which can change their own parameters at the required time moment. In the paper a description of the mathematic model of the proposed solution is presented; the algorithm for determination of the layers parameters of the multiphase liquid is proposed; the calculation example of the model multiphase liquid parameters using the proposed algorithm is shown; the comparison with standard approach is presented. The presented calculation result shows high efficiency of the proposed method, so using the controllable loads in the measuring probe construction is a perspective solution for the task of extracting the multiphase liquids parameters.
References
1. Glebovich G.V. Issledovanie obektov s pomoshchyu pikosekundnyh impulsov [Research of objects using picosecond pulses]. - M.: Radio i svyaz, 1984. - 256 p.
2. Pat. US3474337. Petrick J.R. System for sensing levels and electrical characteristics of fluent materials. Patent, No. US3474337, 1969.
3. Harney W.J., Nemarich C.P. Electromagnetic level indicating (EMLI) system using time domain reflectometry. OCEANS '83. - Proceedings, 1983. - PR 233-236.
4. Nemarich C.P Time domain reflectomerty liquid level sensors // IEEE Instrumentation & Measurement Magazine. -2001. - Vol. 4. - PP. 40-44.
5. Gerding M., Musch T., Schiek B. Precision level measurement based on time-domain reflection (TDR) measurements // Advances in Radio Science. - 2002. - PP. 27-31.
6. Cataldo A. Lay-Ekuakille A., De Carlo C. Remote sensing of liquid characteristics using Time Domain Reflec-tometry // Proceedings of the SPIE. - 2002. - Vol. 4814. -PP. 465-473.
7. Yahya S.I. A multi-level gauging system for crude oil settling mass tanks using a one-port time-domain technique. WIT Transactions on Engineering Sciences. - 2012. - Vol. 81. -PP. 255-266.
8. Cataldo A. Catarinucci L. Tarricone L., Attivissimo F., Trotta A. A frequency-domain method for extending TDR performance in quality determination of fluids // Measurement science and technology. - 2007. - Vol. 18, No. 3. - PP. 675-688.
9. Dozer B.E. Self-calibrating liquid-level measuring device. Patent, No. US3398578, 1968.
9
10. Vanyashin V.A., Zakonov M.A., Zagidullin M.Z., Mezikov A.K., Habibullin R.A. Urovnemer [Level transmitter]. Patent RF, No. 2491519, 2013.
11. Trenkal E.I., Loshchilov A.G. Algoritm selektsii poleznyh signalov v zadache TDR-izmereniya urovney mnogofaznyh zhidkostey [The algorithm of desired signal selection for TDR analysis of multilayer liquids] // Elektronika i mikroehlektronika SVCh. - 2016. - No. 1. - PP. 222-226.
Trenkal Evgeniy Igorevich
Postgraduate student of the chair of design of radio electronic components in Tomsk state university of control system and radioelectronics (TUSUR) Phone: +7-913-816-57-16 E-mail: trenkal@mail.ru
Loshchilov Anton Gennadevich
Ph.D., the head of the chair of design of radio electronic components in Tomsk state university of control system and radioelectronics (TUSUR), the head of special design bureau «Smena» Phone: +7-906-947-69-20 E-mail: yogl@mail.ru
Тренкаль Е.И., Лощилов А.Г.
Новый метод повышения точности TDR-анализа многослойных сред
Предложено решение, позволяющее в режиме реального времени определять структуру и параметры слоев многофазной жидкости, тем самым повысить точность измерения. Метод основан на использовании в конструкции измерительного зонда управляемых нагрузок, способных изменять свои параметры в необходимый момент времени. В работе представлено описание математической модели предложенного решения; предложен алгоритм определения параметров слоев многофазной жидкости; показан пример расчета параметров модельной многофазной жидкости по описанному алгоритму. Также представлено сравнение точности предлагаемого решения и классического подхода измерения.
Ключевые слова: измерение уровня, многофазная жидкость, TDR, способ, измерительный зонд, управляемые нагрузки, измерение параметров.
Доклады ТУСУРа, том 19, № 4, 2016