# Numerical Investigation on ELF Electromagnetic Field Distribution of Pipeline Robot Tracking and Positioning System Using UAV

Authors: Bao Tian; Zhaolin Yang; Shouhong Ji; Xiang Chen; Kun Zhao
Journal: Journal of Sensors
 (2022)

---

## Abstract

Pipeline robot, as a new type of equipment for pipeline operations such as pigging and detection, will play an increasingly important role in the operation and maintenance of oil and gas pipeline networks. The tracking and positioning technology during its operation process is one of the essential topics to improve the operating performance of pipeline robots and eliminate potential pipeline accidents. This paper presents the overall design of the wireless tracking and positioning system for pipeline robots based on extremely low frequency (ELF) electromagnetic method using the unmanned aerial vehicle (UAV). Starting from the classical electromagnetic theory, a mathematical model for the distribution of ELF electromagnetic field in buried metal pipeline environment is established. According to the characteristics of ELF electromagnetic wave transceiver and the model based on classical theory, the equivalent magnetic dipole transmission model is deduced. Based on the equivalent model, the attenuation characteristics of the ELF electromagnetic wave in the external space of the pipeline are analyzed by numerical simulation. The influence of the geometrical dimensions, environmental and working parameters, and other factors on the electromagnetic field intensity outside the pipeline is given at the same time. Finally, the distribution of ELF electromagnetic field in pipeline environment is discussed for horizontal-laying pipeline and inclined-laying pipeline, and the tracking method of pipeline robot is proposed on this basis. The proposed scheme is practical and effective, and it is suitable for real-time tracking and positioning of robot working in the pipelines whose slope is no more than 60 degrees with known or unknown distribution.


---

## Body

## 1. Introduction

As the primary way of oil and natural gas transportation, pipeline plays a vital role in industrial production and civil use. During the laying period and long-term operation, pipelines will inevitably produce dirt sedimentation, deformation, corrosion, and other hidden troubles, which not only reduce the transportation capacity and increase the production cost but also may cause oil or gas leakage and even lead to safety accidents in severe cases. Consequently, strengthening the regular cleaning, testing, and maintenance of oil and gas pipelines is of great significance to ensure the safe operation and prolong the service life of pipelines. Pipeline robot, represented by intelligent pig and internal detector, has been more and more used in the daily maintenance of pipelines due to its unique advantages and is increasingly becoming an essential application carrier in pipeline safety assurance system. However, when the pipeline robot moves in the pipeline, it may get stuck in the pipe due to the pipeline deformation, sediment accumulation, the failure of its system, etc., which may lead to pipeline blockage, affect the normal transportation of oil and gas, and even endanger the safety of pipe network [1]. Therefore, when pipeline robot performs pigging or inspection operations, technicians need to grasp its movement status and accurate position in real-time through effective tracking and positioning technology, so as to take corresponding measures to quickly and accurately handle the accident and resume production.With the development of pipeline robot technology, its tracking and positioning technology has experienced a process from “through instruction”, “accurate positioning” to “tracking” [2, 3]. “Through instruction” means that the pipeline robot sends a specific prompt when passing through a particular position of the pipeline [4]. The typical implementation method is mechanical method, and the equipment used is passage indicator. This method has a low-operating cost, but the false alarm rate is high. It can only obtain the information of whether the robot passes through a specific position, but cannot obtain the accurate position information at its runtime. “Accurate positioning” is to determine the exact position of the robot by using equipment such as transmitters and receivers based on “through instruction” in response to accidents such as the jamming of the robot. Typical implementation methods include magnetic method, acoustic method, pressure method, and radioisotope method. The principle of the magnetic method is to install a transmitter that emits an ELF electromagnetic pulse signal on the robot, detect the pulse signal outside the tube through the receiving coil, and indicate the location of the pipeline robot after processing [5–7]. The advantage of this method is that the signal has strong penetrating power and can penetrate the soil of the pipe wall. It is suitable for many occasions, but it has higher requirements for the sensitivity, response speed, and portability of the receiving device. Acoustic-based positioning methods include noise method, acoustic wave feature discrimination method, and ultrasonic method. Its principle is to record and analyze the acoustic wave characteristics generated by the pipeline robot when it is running or to receive the pulse signal emitted by the ultrasonic transmitter installed on the robot for positioning. The disadvantage is that the signal attenuation is fast, the noise signal interference is significant, and the application range is limited [8, 9]. The pressure method is to realize positioning through different pressure values detected by two pressure transmitters when the robot passes. The pressure method can use the pressure and flow data obtained in the field to judge the operation of the robot, but the monitoring distance is relatively short, which requires high work experience for operators, and its practicability needs further research. The principle of the radioactive isotope method is to install a radioactive source on the robot and indicate the position of the robot by detecting the rays emitted by the radioactive isotope tracer [10]. The method has high-positioning accuracy and strong anti-interference ability, but it will pose a threat to human health and the ecological environment.The technology of “accurate positioning” is generally used only when the robot has a jamming accident. The complexity and diversity of pipeline engineering promote the development of pipeline robot technology towards “autonomous movement,” “internal and external communication,” and “measurable position,” so as to expand its application scope and reduce the risk of pipeline operation. In the process of cleaning, flaw detection, and other procedures of the pipeline robot, when it detects the pipe wall defect or its system failure, the staff outside the pipe needs to know its accurate position in the pipe in time, so as to quickly and accurately troubleshoot and resume operations. Therefore, the realization of wireless, reliable, and real-time “tracking” of pipeline robots has become one of the essential topics to improve the operating performance and practical value of pipeline robots. In theory, both the magnetic method and the pressure method can achieve the “tracking” function. In recent years, the tracking and positioning technology based on optical fiber vibration has also attracted the attention of researchers [11]. The technology needs to lay an optical fiber simultaneously as the pipeline, which will significantly increase the cost and limit the engineering test. Signal processing methods and vibration interference problems still need further research. In summary, the current various real-time tracking and positioning technologies have issues that need to be solved, such as reducing errors, operability, and cost issues. Further research is required at the basic theoretical and application levels.To meet the requirement of real-time tracking and positioning, this paper presents a design scheme of the wireless tracking and positioning system for pipeline robots based on the ELF electromagnetic method using UAV, establishes the distribution model of ELF electromagnetic field in a complex pipeline environment, studies the propagation and attenuation characteristics of ELF electromagnetic wave through numerical and simulation analysis methods, and discusses the tracing algorithm of robot under the two conditions of horizontal-laying pipeline and inclined-laying pipeline, which lays a theoretical foundation for future practical application of the system.


## 2. Pipeline Robot Tracing and Positioning System Using UAV

The tracking and positioning technology studied in this paper includes two aspects: monitoring real-time motion state of the robot and indicating its precise position. To achieve the goal, ELF electromagnetic method with strong penetrating power for metal pipe wall and soil is a better choice. Still, traditional technology requires technicians to follow the pipeline trend for inspection. When the surface of the pipeline section is a harsh geographical environment such as rivers, lakes, alpine jungles, marshes, and gorges, it will bring great difficulties to the tracking and positioning. In order to design a system that does not require personnel tracking and can accurately grasp the real-time position information of the robot, this paper introduces the UAV technology, which has recently become a research hotspot. By exploiting its advantages such as flexibility, wide inspection range, and timely information transmission to the full, the pipeline robot can be tracked and positioned in real time, accurately and quickly.The wireless tracking and positioning system of the pipeline robot using UAV is shown in Figure1. Its working principle is as follows: an ELF electromagnetic wave transmitter is installed on the pipeline robot, and the periodic emission of ELF electromagnetic pulse is realized by transmitting antenna (solenoid coil); the time-varying ELF electromagnetic field is established in the buried pipeline-ground environment by the strong penetrability of ELF wave; the receiving antenna (solenoid coil) carried by the tracking UAV converts the electromagnetic signal penetrating the metal pipe wall, overlying soil and air into voltage signal. After amplification and filtering, the electromagnetic field intensity is measured, and the location of ELF electromagnetic wave emitter is determined, and then the robot’s position is determined.Figure 1
Wireless tracing and positioning system of pipeline robot using UAV.When the pipeline robot gets stuck, the UAV can quickly find and approach the accident point according to the ELF electromagnetic field signal detected and report the location and environment of the accident to the operators at the upstream and downstream stations through the airborne wireless communication system. The operators can rapidly formulate the accident treatment plan based on the location, environmental photos of the incident site, and other information received and take timely actions, which significantly saves staffing and time. If the conditions permit, the UAV can work in tracking mode. The airborne-receiving device detects the ELF electromagnetic field signal in real-time, sets the parameters of flight control system according to the solution result to maintain the synchronization of flight tracking path and the robot, and simultaneously reports the information such as path and speed to operators through wireless communication, so that operators can keep track of the movement status of drone and pipeline robot. In the tracking process, if the abnormal motion status of the robot is found, such as moving more slowly, etc., the alarm can be sent out to remind the operator to make an emergency action plan in time.The hardware part of the system can be realized by installing ELF signal transmitting on UAV and receiving modules on existing robots. The UAV equipped with a receiver to track and locate the robot can be flexibly selected by carrying a receiving module according to the working environment and conditions such as monitoring distance and flight environment, without developing a unique UAV system, which will significantly reduce the initial cost of the system. At the same time, the system can select the accident location positioning or real-time tracking according to the operation purpose. Based on experience, the robot does not need to work every time for real-time tracking, more importantly, in the event of a robot failure or found abnormal pipeline can accurately obtain its position information. Compared with the current precise positioning method, the drone can find the target more quickly and accurately. The operation range is incomparable to other methods, and the total operation cost is lower. Currently, the latest optical fiber early warning system also has high-positioning accuracy, extensive monitoring range, and low-operating cost, but the initial cost is enormous. It needs to lay a fiber with the pipeline at the same time and needs to build a distributed optical fiber-sensing system [12]. Especially for buried pipelines that have already been laid, reexcavation to lay the fiber is not the best option at extra cost. Therefore, considering the positioning effect, application difficulty, and economic factors, the pipeline robot tracking and positioning system based on UAV has unique advantages.According to the working principle of the system, it can be seen that establishing the mathematical model of ELF electromagnetic field in pipeline environment and analyzing its distribution and attenuation characteristics are the theoretical basis for studying the position relation between the signal received by UAV (receiving coil) and pipeline robot (transmitting coil) to realize tracking and positioning, as well as the basis for the development algorithm of the application software of the tracking and positioning system.


## 3. Distribution Model of the ELF Magnetic Field in the Buried Pipe-Ground Environment

### 3.1. Establishment of Electromagnetic Field Distribution Model

The environment geometry of the tracking and positioning system is shown in Figure2. M1 in the pipeline is the transmission medium (oil or natural gas), pipe wall M2 is the ferromagnetic material, M3 is the overlying soil above the pipeline, M4 is the air layer, and the anticorrosion insulation layer is ignored. The electromagnetic parameters of each layer are εi,μi,σi (i=1, 2, 3, and 4), εi is the permittivity, μi is the permeability, σi is the conductivity, and assuming that each layer is homogeneous and isotropic. The inner diameter of the buried pipeline is r1, the external diameter is r2, the distance between the surface and pipeline axis is r3, and the thickness of the buried medium (overlying soil) is h=r3−r2. The transmitting solenoid coil S on the robot is placed on the central axis of the pipeline by default, the distance from the surface is r3; its radius a is far less than the distance between the transmitting coil and the receiving coil and far less than the length of the transmitting coil 2 l; the number of coil turns is N, and the number of turns in unit length is n; the exciting current is I=I0e−jωt, where ω is the current angular frequency. The cylindrical-coordinate system is established with the coil center as the origin, and the Z axis is on the central axis of the pipeline. The coordinate of the receiving coil on the UAV is Pρ,θ,z, the space vector from the transmitting coil center to the receiving coil is R=rer, and the distance is r=ρ2+z21/2. The magnetic field distribution generated by the transmitting coil is axisymmetric.Figure 2
The environment geometry of the tracking and positioning system.The electric field intensity in the four-layer medium isE1, E2, E3, and E4, and the magnetic field intensity is H1, H2, H3, and H4, respectively. Due to the different properties of the two media, the field quantities varies at the interface. Knowing from the electromagnetic theory, the tangential and normal components of electric field E, D and H, B of magnetic field are continuous at the interface of transmission medium-pipe wall, pipe wall-overlying soil, and overlying soil-air in the cylindrical-coordinate system. When the transmitting solenoid coil is placed horizontally along the axis of the pipe, for the axisymmetric medium environment, the electric field intensity has only Eθ component, while the magnetic field intensity has Hρ and Hz components. According to the Helmholtz equation, the expression of electric field intensity in each layer can be obtained by using the method of separation of variables as follows:
(1)E1θρ,z=−jωμ1M2π2∫0∞A1λJ1ρx1+x1Y1ρx1cosλzdλ,E2θρ,z=−jωμ2M2π2∫0∞A2λJ1ρx2+B2λY1ρx2cosλzdλ,E3θρ,z=−jωμ3M2π2∫0∞A3λJ1ρx3+B3λY1ρx3cosλzdλ,E4θρ,z=−jωμ4M2π2∫0∞B4λY1ρx4cosλzdλ.FromH=1/jωμ∇×E, the expression of magnetic field intensity in each layer can be obtained as follows:
(2)H=eρ−1jωμ∂Eθ∂z+ez1jωμEθρ+∂Eθ∂ρ.The mathematical model of ELF electromagnetic field distribution under buried pipeline environment is obtained by expanding the expression as follows:(3)H1ρρ,z=−M2π2∫0∞λA1λJ1ρx1+x1Y1ρx1sinλzdλ,H1zρ,z=−M2π2∫0∞x1−A1λJ0ρx1+x1Y0ρx1cosλzdλ,H2ρρ,z=−M2π2∫0∞λA2λJ1ρx2+B2λY1ρx2sinλzdλ,H2zρ,z=−M2π2∫0∞x2−A2λJ0ρx2+B2λY0ρx2cosλzdλ,H3ρρ,z=−M2π2∫0∞λAλJ1ρx3+B3λY1ρx3sinλzdλ,H3zρ,z=−M2π2∫0∞x3−A3λJ0ρx3+B3λY0ρx3cosλzdλ,H4ρρ,z=−M2π2∫0∞λB4λY1ρx4sinλzdλ,H4zρ,z=−M2π2∫0∞x4B4λY0ρx4cosλzdλ,where M=πa2NI is the magnetic moment of the transmitting solenoid coil, λ is an arbitrary real variable, J0 and J1 are the first kinds of modified Bessel functions, Y0 and Y1 are the second kinds of modified Bessel functions, and Ai and Bi are the undetermined coefficients which are derived from initial excitation conditions and boundary conditions, xi=λ2−ω2μiεi+jωμiσi21/2.Therefore, when the receiving coilPρ,θ,z in the UAV is placed vertically (the transmitting coil is parallel to the pipe axis and the receiving coil is perpendicular to the transmitting coil, as shown in Figure 3(a)), the magnetic field component of the ρ direction is mainly received, and the induced electromotive force obtained is as follows:
(4)εvρ,θ,z=−∂ΦV∂t=−jωμ4MAcosθv2π2∫0∞x4B4λY0ρx4cosλzdλ.Figure 3
Coil relative position diagram.
(a)(b)When the receiving coilPρ,θ,z in the UAV is placed horizontally (the transmitting coil is parallel to the pipe axis, and the receiving coil is parallel to the transmitting coil, as shown in Figure 3(b)), the magnetic field component of the z direction is mainly received, and the induced electromotive force obtained is as follows:
(5)εhρ,θ,z=−∂ΦH∂t=−jωμ4MAcosθh2π2∫0∞λB4λY1ρx4sinλzdλ,where A is the cross-sectional area of the receiving coil, θv is the angle between the cross-sectional normal and the pipeline axial direction when the receiving coil is placed vertically on the ground, and θh is the angle between the cross-sectional normal and the pipeline radial direction when the receiving coil is placed horizontally.


### 3.2. Dipole Equivalent Magnetic Model

It can be seen from the above derivation process that the inversion method of transmitting source position by using the precise model to calculate the electromagnetic distribution at a certain point provides a theoretical basis for the subsequent research. Still, the computational process is too complex and not conducive to real-time analysis. Therefore, it is necessary to seek an equivalent model for analysis and calculation.When the geometrical size of the transmitting coil (solenoid) is much smaller than the detection range, the length of the solenoid is much larger than the radius of the solenoid, and the excitation magnetic field is ELF field (or quasistatic field); the electromagnetic field generated by the excitation source can be equivalent to that generated by a magnetic dipole composed of a pair of point magnetic charges±qm with equal quantity and opposite sign [13].Bessel function in the precise model, Formula (3) is written in power series first, and the expression is expanded by Bessel integral formula. Under the far-field condition, i.e., ρ,z≫a,l, the series expression converges rapidly, so only the n=0 term is needed in the series expressions of radial and axial magnetic field intensity, which is the precondition of using the equivalent magnetic dipole model. When n=0, the radial and axial components of the magnetic induction intensity at the receiving coil Pρ,θ,z can be deduced as follows:
(6)Bρρ,z=nμ0Ia2ρ41r2−2lz3/2−1r2+2lz3/2=μ0m4πr5∙3ρz,Bzρ,z=nμ0Ia24z−lr2−2lz3/2−z+lr2+2lz3/2=nμ0Ia24r32lr22z2−ρ2=μ0m4πr5∙2z2−ρ2,i.e.,(7)Bs=Ber=Bρ2+Bz21/2er=μ04πr53m∙RR−r5m,where m=2lnπa2I is the magnetic moment of the transmitting solenoid coil.The Formula (7) is precisely an expression of the magnetic field generated by a magnetic dipole with the magnetic moment m=ezm in cylindrical-coordinate system. It shows that it is theoretically feasible to use the equivalent magnetic dipole model when calculating the magnetic field in the far-field region.According toHρρ,z, Hzρ,z obtained above, numerical simulation of spatial magnetic field distribution is carried out in the X-o-Z plane. The selected parameters are as follows: n=10turns/mm, l=500mm, a=50mm, I=0.3sin46πtA, taking the free space permeability μ0=4π×10−7. When the receiving coil is ρ=8m away from the pipeline axis, the variation curve of the amplitude H of the magnetic field intensity component with z is shown in Figure 4.Figure 4
Variation of magnetic field intensityH with z.
(a)
Radial direction(b)
Axial directionIt can be seen from Figure4(a) that the component Hρ of the spatial magnetic induction intensity in the ρ direction generated by the equivalent magnetic dipole model presents a double-peak symmetric distribution in the axial direction of the pipeline. Near the center of the transmitting coil, the component of magnetic induction intensity in the ρ direction is the global minimum value; as the distance from the central point increases, the element of the magnetic induction intensity in this direction also increases, and two maximum signals can be detected at the front and rear positions of the transmitting coil. That is, in the case of single-coil reception, when the receiving coil moves directly above the transmitting coil, the detected signal is the weakest; when the receiving coil moves away from the transmitting coil, the detected signal is strengthened accordingly, which is the theoretical basis of positioning the pipeline robot with ELF electromagnetic wave. It can be seen from Figure 4(b) that the component Hz of magnetic induction intensity in the z direction presents a main-peak symmetric distribution and has a global maximum value near the center of the transmitting coil; that is, a global maximum signal can be measured in the z direction when the receiving coil is directly above the transmitting antenna, and there are two local maximum signals at the front and rear positions of the transmitting coil. The depth information of the transmitting coil can be obtained by measuring the amplitude width of the global maximum signal.


## 3.1. Establishment of Electromagnetic Field Distribution Model

The environment geometry of the tracking and positioning system is shown in Figure2. M1 in the pipeline is the transmission medium (oil or natural gas), pipe wall M2 is the ferromagnetic material, M3 is the overlying soil above the pipeline, M4 is the air layer, and the anticorrosion insulation layer is ignored. The electromagnetic parameters of each layer are εi,μi,σi (i=1, 2, 3, and 4), εi is the permittivity, μi is the permeability, σi is the conductivity, and assuming that each layer is homogeneous and isotropic. The inner diameter of the buried pipeline is r1, the external diameter is r2, the distance between the surface and pipeline axis is r3, and the thickness of the buried medium (overlying soil) is h=r3−r2. The transmitting solenoid coil S on the robot is placed on the central axis of the pipeline by default, the distance from the surface is r3; its radius a is far less than the distance between the transmitting coil and the receiving coil and far less than the length of the transmitting coil 2 l; the number of coil turns is N, and the number of turns in unit length is n; the exciting current is I=I0e−jωt, where ω is the current angular frequency. The cylindrical-coordinate system is established with the coil center as the origin, and the Z axis is on the central axis of the pipeline. The coordinate of the receiving coil on the UAV is Pρ,θ,z, the space vector from the transmitting coil center to the receiving coil is R=rer, and the distance is r=ρ2+z21/2. The magnetic field distribution generated by the transmitting coil is axisymmetric.Figure 2
The environment geometry of the tracking and positioning system.The electric field intensity in the four-layer medium isE1, E2, E3, and E4, and the magnetic field intensity is H1, H2, H3, and H4, respectively. Due to the different properties of the two media, the field quantities varies at the interface. Knowing from the electromagnetic theory, the tangential and normal components of electric field E, D and H, B of magnetic field are continuous at the interface of transmission medium-pipe wall, pipe wall-overlying soil, and overlying soil-air in the cylindrical-coordinate system. When the transmitting solenoid coil is placed horizontally along the axis of the pipe, for the axisymmetric medium environment, the electric field intensity has only Eθ component, while the magnetic field intensity has Hρ and Hz components. According to the Helmholtz equation, the expression of electric field intensity in each layer can be obtained by using the method of separation of variables as follows:
(1)E1θρ,z=−jωμ1M2π2∫0∞A1λJ1ρx1+x1Y1ρx1cosλzdλ,E2θρ,z=−jωμ2M2π2∫0∞A2λJ1ρx2+B2λY1ρx2cosλzdλ,E3θρ,z=−jωμ3M2π2∫0∞A3λJ1ρx3+B3λY1ρx3cosλzdλ,E4θρ,z=−jωμ4M2π2∫0∞B4λY1ρx4cosλzdλ.FromH=1/jωμ∇×E, the expression of magnetic field intensity in each layer can be obtained as follows:
(2)H=eρ−1jωμ∂Eθ∂z+ez1jωμEθρ+∂Eθ∂ρ.The mathematical model of ELF electromagnetic field distribution under buried pipeline environment is obtained by expanding the expression as follows:(3)H1ρρ,z=−M2π2∫0∞λA1λJ1ρx1+x1Y1ρx1sinλzdλ,H1zρ,z=−M2π2∫0∞x1−A1λJ0ρx1+x1Y0ρx1cosλzdλ,H2ρρ,z=−M2π2∫0∞λA2λJ1ρx2+B2λY1ρx2sinλzdλ,H2zρ,z=−M2π2∫0∞x2−A2λJ0ρx2+B2λY0ρx2cosλzdλ,H3ρρ,z=−M2π2∫0∞λAλJ1ρx3+B3λY1ρx3sinλzdλ,H3zρ,z=−M2π2∫0∞x3−A3λJ0ρx3+B3λY0ρx3cosλzdλ,H4ρρ,z=−M2π2∫0∞λB4λY1ρx4sinλzdλ,H4zρ,z=−M2π2∫0∞x4B4λY0ρx4cosλzdλ,where M=πa2NI is the magnetic moment of the transmitting solenoid coil, λ is an arbitrary real variable, J0 and J1 are the first kinds of modified Bessel functions, Y0 and Y1 are the second kinds of modified Bessel functions, and Ai and Bi are the undetermined coefficients which are derived from initial excitation conditions and boundary conditions, xi=λ2−ω2μiεi+jωμiσi21/2.Therefore, when the receiving coilPρ,θ,z in the UAV is placed vertically (the transmitting coil is parallel to the pipe axis and the receiving coil is perpendicular to the transmitting coil, as shown in Figure 3(a)), the magnetic field component of the ρ direction is mainly received, and the induced electromotive force obtained is as follows:
(4)εvρ,θ,z=−∂ΦV∂t=−jωμ4MAcosθv2π2∫0∞x4B4λY0ρx4cosλzdλ.Figure 3
Coil relative position diagram.
(a)(b)When the receiving coilPρ,θ,z in the UAV is placed horizontally (the transmitting coil is parallel to the pipe axis, and the receiving coil is parallel to the transmitting coil, as shown in Figure 3(b)), the magnetic field component of the z direction is mainly received, and the induced electromotive force obtained is as follows:
(5)εhρ,θ,z=−∂ΦH∂t=−jωμ4MAcosθh2π2∫0∞λB4λY1ρx4sinλzdλ,where A is the cross-sectional area of the receiving coil, θv is the angle between the cross-sectional normal and the pipeline axial direction when the receiving coil is placed vertically on the ground, and θh is the angle between the cross-sectional normal and the pipeline radial direction when the receiving coil is placed horizontally.


## 3.2. Dipole Equivalent Magnetic Model

It can be seen from the above derivation process that the inversion method of transmitting source position by using the precise model to calculate the electromagnetic distribution at a certain point provides a theoretical basis for the subsequent research. Still, the computational process is too complex and not conducive to real-time analysis. Therefore, it is necessary to seek an equivalent model for analysis and calculation.When the geometrical size of the transmitting coil (solenoid) is much smaller than the detection range, the length of the solenoid is much larger than the radius of the solenoid, and the excitation magnetic field is ELF field (or quasistatic field); the electromagnetic field generated by the excitation source can be equivalent to that generated by a magnetic dipole composed of a pair of point magnetic charges±qm with equal quantity and opposite sign [13].Bessel function in the precise model, Formula (3) is written in power series first, and the expression is expanded by Bessel integral formula. Under the far-field condition, i.e., ρ,z≫a,l, the series expression converges rapidly, so only the n=0 term is needed in the series expressions of radial and axial magnetic field intensity, which is the precondition of using the equivalent magnetic dipole model. When n=0, the radial and axial components of the magnetic induction intensity at the receiving coil Pρ,θ,z can be deduced as follows:
(6)Bρρ,z=nμ0Ia2ρ41r2−2lz3/2−1r2+2lz3/2=μ0m4πr5∙3ρz,Bzρ,z=nμ0Ia24z−lr2−2lz3/2−z+lr2+2lz3/2=nμ0Ia24r32lr22z2−ρ2=μ0m4πr5∙2z2−ρ2,i.e.,(7)Bs=Ber=Bρ2+Bz21/2er=μ04πr53m∙RR−r5m,where m=2lnπa2I is the magnetic moment of the transmitting solenoid coil.The Formula (7) is precisely an expression of the magnetic field generated by a magnetic dipole with the magnetic moment m=ezm in cylindrical-coordinate system. It shows that it is theoretically feasible to use the equivalent magnetic dipole model when calculating the magnetic field in the far-field region.According toHρρ,z, Hzρ,z obtained above, numerical simulation of spatial magnetic field distribution is carried out in the X-o-Z plane. The selected parameters are as follows: n=10turns/mm, l=500mm, a=50mm, I=0.3sin46πtA, taking the free space permeability μ0=4π×10−7. When the receiving coil is ρ=8m away from the pipeline axis, the variation curve of the amplitude H of the magnetic field intensity component with z is shown in Figure 4.Figure 4
Variation of magnetic field intensityH with z.
(a)
Radial direction(b)
Axial directionIt can be seen from Figure4(a) that the component Hρ of the spatial magnetic induction intensity in the ρ direction generated by the equivalent magnetic dipole model presents a double-peak symmetric distribution in the axial direction of the pipeline. Near the center of the transmitting coil, the component of magnetic induction intensity in the ρ direction is the global minimum value; as the distance from the central point increases, the element of the magnetic induction intensity in this direction also increases, and two maximum signals can be detected at the front and rear positions of the transmitting coil. That is, in the case of single-coil reception, when the receiving coil moves directly above the transmitting coil, the detected signal is the weakest; when the receiving coil moves away from the transmitting coil, the detected signal is strengthened accordingly, which is the theoretical basis of positioning the pipeline robot with ELF electromagnetic wave. It can be seen from Figure 4(b) that the component Hz of magnetic induction intensity in the z direction presents a main-peak symmetric distribution and has a global maximum value near the center of the transmitting coil; that is, a global maximum signal can be measured in the z direction when the receiving coil is directly above the transmitting antenna, and there are two local maximum signals at the front and rear positions of the transmitting coil. The depth information of the transmitting coil can be obtained by measuring the amplitude width of the global maximum signal.


## 4. ELF Electromagnetic Wave Attenuation Model

ELF electromagnetic field propagates outward in the form of wave. When ELF electromagnetic wave passes through conducting medium, part of it is absorbed by the medium, which causes electromagnetic wave attenuation, and part of it is consumed in the form of heat. In this section, the attenuation model of ELF electromagnetic wave in the buried pipe environment is derived using the equivalent magnetic dipole model.The expression of the ELF magnetic field intensity in lossy medium can be obtained from the equivalent dipole model according to the wave equation(8)Hr=Hse−αde−jβd=Bsμ0e−αde−jβd,where α is the attenuation constant, indicating the attenuation of the amplitude of the wave per unit length in the direction of propagation, and β is the phase-shift constant, meaning the change of the wave’s phase per unit length in the direction of propagation.
(9)α=ωμε′21+ε′′ε′2−11/2,β=ωμε′21+ε′′ε′2+11/2,whereε’=ε, ε”=σ/ω, and ω=2πf are angular frequency, μ is magnetic permeability, ε is the dielectric constant of the medium, and σ is the electrical conductivity of the medium. The phase lag effect in low frequency electromagnetic wave transmission can be neglected, and the attenuation of amplitude is only needed care. The expression of α can be simplified for different mediums passing through the ELF electromagnetic wave propagation path.In practice, ifε”/ε’>102, the medium can be considered as a good conductor, such as pipe wall and overlying soil (especially for moist soil), and α can be approximated as
(10)α≈ωμε′′2=ωμσ2=πfμσ.It can be seen from the Formula (10) that α increases with the increase of σ and transmitting frequency f. The ability of electromagnetic wave to penetrate the conductor is described by skin depth δ=1/α; that is, the thickness of electromagnetic wave propagating when the amplitude of wave decreases to 0.37 times of the initial value. In order to penetrate the metal pipe wall, the thickness of the pipe wall must be guaranteed to be less than the skin depth. In order to penetrate the metal pipe wall, it must ensure that the wall thickness is less than the skin depth.
(11)d<δ=1α=1πfμσ.Therefore,f<1/d2πμσ indicates that increasing the skin depth by reducing the emission frequency can improve the penetration ability of electromagnetic waves.Ifε”/ε’<10−2, the medium can be considered as low-loss dielectrics, such as transmission fluid and air, and α can be approximated as
(12)α≈ωε′′2με′=σ2με.In low-loss dielectric,β is the same as in the case of lossless, and α is independent of frequency. Then, in the buried pipeline environment, the pipe wall is regarded as a good conductor, and the transmission fluid, soil, and air are considered as low-loss dielectric. The expression of ELF electromagnetic wave from the transmitting coil to the receiving coil obtained by using the equivalent magnetic dipole model is as follows:
(13)Hr=Hse−αde−jβd=14πr53m∙RR−r2me−r1σ1/2μ1/ε1e−πfμ2σ2de−πfμ3σ3he−ρ−r3σ3/2μ3/ε3.


## 5. Analysis of the Influencing Factors of ELF Electromagnetic Field

Based on the above analysis, this section studies the effects of environmental parameters and working parameters on the ELF electromagnetic field intensity through numerical simulation. The research results will provide a theoretical basis for the selection of transmitter frequency, the optimization design of transmitting coil and receiving coil, the subsequent positioning algorithm research, and field experiments in the future.In the buried pipeline environment, the electromagnetic constitutive parameters of each layer through which the electromagnetic wave passes are shown in Table1. For both crude oil pipelines and gas pipelines, e−r1σ1/2μ1/ε1≈1 can be obtained by substituting the parameters in the table; that is, the medium in the pipeline will not cause the attenuation of electromagnetic wave. Similarly, for the air layer, e−ρ−r3σ3/2μ3/ε3≈1; that is, the air layer will not cause the rapid attenuation too. Therefore, the following analysis mainly considers the influence of pipe wall and overlying soil layer on electromagnetic wave transmission. The selection of transmitting antenna and environment parameters is shown in Tables 2 and 3.Table 1
Electromagnetic constitutive parameters (low frequency and room temperature).
MediumRelative permittivity (εr)Conductivity (σ)Relative permeability (μr)Oil M12.13×10−131Gas M133.11×10−161Metal pipe (carbon steel) M2∞1×10710-300Overlying soil (moist) M3161×10−21Air M41.00061×10−161Notes:ε=εrε0, ε0=8.854×10−12F/m, μ=μrμ0, and μ0=4π×10−7H/m.Table 2
Transmitting antenna parameters.
Length (mm)Coil radius (mm)Coil turns (mm)Excitation current (A)Frequency (Hz)50050200.3sin(50πt)A10 ~ 100Table 3
Environment parameters.
Internal diameter of the pipeline (mm)Pipe thickness (mm)Overlying soil thickness (m)Elevation (m)8000 ~ 600 ~ 65Without losing generality, the simulation calculation is carried out with the receiving antenna located directly above the transmitting antenna. The influence of different pipe thicknesses and overlying soil thicknesses on the received magnetic field intensity is analyzed, as shown in Figures5 and 6.Figure 5
Effect of pipe thickness.Figure 6
Effect of overlying soil thickness.Figure5 shows the exponential decay curve of the received magnetic field intensity with the increase of the pipe thickness when the transmitting frequency is 25 Hz and the overlying soil thickness is 2 m. Figure 6 shows the exponential decay curve of the received magnetic field intensity with the increase of overlying soil thickness when the transmitting frequency is 25 Hz and pipe thickness is 12 mm. It can be seen from the gradient of the two decay curves that the pipe thickness has a more significant attenuation on the magnetic field intensity. When designing the transmitting and receiving antenna, the influence of the pipe size and burial depth on the robot positioning should be fully considered.Figure7 shows the decay curve of the received magnetic field intensity with the increase of transmitting frequency when pipe thickness is 12 mm and overlying soil thickness is 3 m. The attenuation of the magnetic field intensity is more linear and smoother. However, the lower the frequency is not the better. For a transmitting antenna in a moving state, too low frequency will reduce the resolution of the receiving antenna signal processing. Figure 8 shows the exponential decay curves of the received magnetic field intensity with the increase of relative permeability and conductivity when the transmitting frequency is 25 Hz, pipe thickness is 12 mm, and overlying soil thickness is 3 m. The electromagnetic parameters, such as relative permeability and conductivity, are determined by the material properties of pipelines, including the carbon content and the doping of alloys. In addition, the variation characteristics of the surrounding conductivity and permeability of pipes with cathodic protection or corrosion need special attention.Figure 7
Effect of transmitting frequency.Figure 8
Effect of electromagnetic parameter.


## 6. Tracking and Positioning Method under Different Pipeline-Laying Conditions

The deduction and analysis are carried out when the transmitting antenna and the receiving antenna are perpendicular or parallel to each other and in the same plane. Still, in more general, the laying of buried or submarine pipelines is more complex and changeable. When the robot runs in the pipe, the axis direction of the transmitting coil can be approximately parallel to the axis direction of the pipeline, but the direction of the pipeline itself is constantly changing, so the relative position between the receiving antenna outside the pipeline (horizontal and vertical) and the transmitting coil is also changing, and the distribution law of the spatial magnetic field is also changing. Then, the positioning method of the pipeline robot is discussed in more general cases aiming at horizontal laying of pipelines and inclined laying of pipelines.


### 6.1. Horizontal Pipeline Environment

The distribution of the magnetic field in horizontal laying of pipeline is analyzed first. The axis of the transmitting solenoid coil is always parallel to the axis of pipeline by default; that is, the transmitting antenna is in the horizontal position. The receiving antenna is still in the vertical or horizontal direction to measure the horizontal and vertical components of magnetic field intensity, but it is not necessarily in the same plane as the transmitting antenna. The combination of horizontal transmission and vertical reception is shown in Figure9. For the convenience of intuitive analysis, the output voltage of the receiver is directly used to represent the received magnetic field intensity, and the influence of the different relative positions on the voltage change is mainly discussed.Figure 9
Horizontal transmission-vertical reception.Based on the analysis of the ELF electromagnetic field distribution model, the components received in three directions (regardless of the position of the receiving antenna) can be obtained(14)Hx=3M4πl2rsin3αcosαcos3β,Hy=−3M4πl2rsin2αcosαsinβcos2β,HzM4πl2rsin2αcos2β3cos3α−1.The relationship between the output voltage of the sensor and the relative position of the transmitting antenna and the receiving antenna (represented bya, β, l, and r) is as follows:
(15)Uv=CHx=3CM4πl2rsin3αcosαcos3β,where C is the signal gain coefficient of the receiving sensor, and in this case, C=150V/T; M is the equivalent magnetic moment considering the attenuation caused by the pipe thickness and overlying soil, taking M=30A·m2; l=5m, α∈[0, 180°], and β∈−90°,90°.The distribution and plane isogram ofUv with the change of relative position under horizontal transmission-vertical reception is shown in Figure 10. It can be seen from the figure that the distribution of Uv shows a double-peak characteristic, which is the same as the distribution feature on the X-o-Z plane obtained in the preceding. The UAV receives the electromagnetic field signal during the positioning process. According to the distribution regulation shown in the figure, UAV always travels along the gradient direction (the direction of the fastest descent or ascent of signal intensity) and explores the peak and trough positions shown in the figure. The midpoint of the peak and trough is the position directly above the robot. If the voltage is taken as an absolute value, the trough becomes a peak, and the robot is at the midpoint of the two peaks.Figure 10
Distribution and plane isogram ofUv under horizontal transmission-vertical reception.When the horizontal-receiving mode is adopted, the receiving antenna is parallel to the ground during the flight, but the course is not necessarily the same as the direction of the pipeline; that is, the axis of the receiving antenna and the axis of the transmitting antenna may produce a certain angleγ on the same horizontal plane. γ=0°, the receiving antenna axis, is parallel to the transmitting antenna axis, and γ=90°, the receiving antenna axis, is vertical to the transmitting antenna axis, as shown in Figure 11.Figure 11
Horizontal transmission-horizontal reception.The output voltage of the sensor under this working condition is(16)Uh=CHh=CHzcosγ−Hysinγ=3Cm4πl2rsin2αcos2β3cos3α−1cosγ−cosαsinβsinγ.The distribution and plane isogram ofUh with the change of relative position under horizontal transmission-horizontal reception is shown in Figure 12. When γ=0°, it can be seen that the distribution of Uh has single peak and double troughs and the absolute output voltage has a global maximum value with two local maximum values on both sides of it, which is the same as the analysis conclusion on the two-dimensional plane. When tracking and positioning, the UAV searches for the global maximum value along the steepest ascent curve of the output voltage amplitude to locate the pipeline robot. When γ increases gradually, the contour deflects and the positions of local maximum values rotate. The troughs representing the global maximum gradually change from single trough to connected trough, and finally forms double troughs. When γ=90°; that is, the axis of receiving antenna is perpendicular to the direction of the pipe; two global maximum values and two local maximum values are distributed at an equal interval of 45°. To sum up, the vertical reception should be used primarily in the positioning process by finding the midpoint of the peak and trough. Then, the horizontal reception is used as the supplement by searching for the maximum value of the sensor output signal to adjust the course of the UAV and assist the position verification.Figure 12
Distribution and plane isogram ofUh under horizontal transmission-horizontal reception.
(a)
γ=0°(b)
γ=25°(c)
γ=75°(d)
γ=90°


### 6.2. Inclined Pipeline Environment

In the practical long-distance pipeline-laying process, affected by rugged terrain, obstacles, and other factors, it is difficult for the pipeline to keep parallel with the horizontal plane for a long distance, and there are a lot of slopes and corners. Therefore, it is necessary to study the electromagnetic field distribution characteristics in the case of the inclined pipeline. As shown in Figure13, the angle between the pipeline and the horizontal plane (inclination angle) is θ, and θ∈0°~90° has the following relationship:
(17)r=lsinθcosα+cosθsinαcosβ.Figure 13
Inclined transmission-vertical reception.The components received in three directions can be obtained(18)Hx=3M4πl3sinαcosαcosβsinθcosα+cosθsinαcosβ3,Hy=3M4πl3sinαcosαsinβsinθcosα+cosθsinαcosβ3,HzM4πl33cos2α+1sinθcosα+cosθsinαcosβ3.Considering the vertical-receiving mode shown in Figure13, the vertical magnetic field intensity is as follows:
(19)Hv=Hxcosθ+Hzsinθ.The relationship between the output voltage of the sensor and the relative position of the pipeline-receiving antenna is as follows:(20)Uv=CHv=CM4πl3sinθ2cos2α−sin2α+3cosθsinαcosαcosβsinθcosα+cosθsinαcosβ3.The distribution and plane isogram ofUv with the change of relative position under inclined transmission-vertical reception is shown in Figure 14. It can be seen from the figure that when the gradient of pipeline is small, such as θ=10 degrees (see Figures 14(a) and 14(b)), the peak-trough characteristic is still apparent, and the pipeline robot can be located by finding the midpoint of the wave peak and trough; when the pipeline gradient increases gradually (see Figure 14(c)), the wave peak topples and falls to the side of the trough, and the trough rises gradually (corresponding amplitude gets gradually smaller); when θ>60° (see Figure 14(d)), the trough disappears and the peak appears only at the position where the receiver is far from the robot; and when θ>80° (Figure 14(e)), the wave peak and trough have disappeared, and it is not easy to locate the position of the robot. However, considering that the gradient of long-distance pipeline is generally not more than 45° in practice, ELF tracking and positioning method is still available.Figure 14
Distribution and plane isogram ofUv under inclined transmission-vertical reception.
(a)
θ=10°(b)
θ=20°(c)
θ=60°(d)
θ=80°When the receiving antenna is parallel to the ground, as shown in Figure15, the angle between the axis of the receiving antenna and the axis of the transmitting is still γ. The magnetic field intensity in the horizontal direction is
(21)Hh=Hzcosθcosγ−Hysinγ−Hxsinθcosγ.Figure 15
Inclined transmission-horizontal reception.The output voltage of the sensor under this working condition by substituting Formula (18) is
(22)Uh=CHh=CM4πl3cosθ2cos2α−sin2αcosγ−3cosθsinαcosαcosβcosγ+3sinαcosαcosβsinγsinθcosα+cosθsinαcosβ3.There are four variablesα, β, γ, and θ in the above formula, so the distribution of Uh is complicated. Figure 16 shows the distribution and plane isogram of Uh when θ=20°, 40°, and 60° under γ=0°, 30°, 60°, and 90°, respectively. It can be seen from the figure that when θ=20° (see Figure 16(a)), the distribution of the global maximum (trough) and local maximum (peak) is pronounced. With the increase of γ, the positions of the two peaks are deflected, and the trough is gradually split into two troughs, and finally evolve into double peaks and double troughs. The robot is at the midpoint of the peak and trough, which can be used for positioning. When θ=40° (see Figure 16(b)), there is only one global maximum value (single trough) under γ=0°. The position of the trough deflects to the rear of the robot, and a local maximum (peak) is generated, which increases with γ. The axial position of the pipe can be judged by the midpoint of the wave and trough. Although the precise position of the robot cannot be located, it can be considered whether it is in front or rear of the UAV. When θ=60° (see Figure 16(c)) and if γ is small, only the global maximum (single trough) can no longer locate the robot and the direction of the pipeline cannot be determined; if γ is further increased, a local maximum value (peak) appears, and the boundary line of the peak and valley is the pipe axis. Therefore, when the pipeline gradient is large, the horizontal-receiving antenna can hardly locate the pipeline robot. Still, it can help the UAV to return to the pipeline direction when it is yawing. At the same time, as mentioned above, the practical pipeline does not have such a large gradient.Figure 16
Distribution and plane isogram ofUh under inclined transmission-horizontal reception.
(a)
θ=20°, γ=0°, 30°, 60°, and 90°(b)
θ=40°, γ=0°, 30°, 60°, and 90°(c)
θ=60°, γ=0°, 30°, 60°, and 90°


## 6.1. Horizontal Pipeline Environment

The distribution of the magnetic field in horizontal laying of pipeline is analyzed first. The axis of the transmitting solenoid coil is always parallel to the axis of pipeline by default; that is, the transmitting antenna is in the horizontal position. The receiving antenna is still in the vertical or horizontal direction to measure the horizontal and vertical components of magnetic field intensity, but it is not necessarily in the same plane as the transmitting antenna. The combination of horizontal transmission and vertical reception is shown in Figure9. For the convenience of intuitive analysis, the output voltage of the receiver is directly used to represent the received magnetic field intensity, and the influence of the different relative positions on the voltage change is mainly discussed.Figure 9
Horizontal transmission-vertical reception.Based on the analysis of the ELF electromagnetic field distribution model, the components received in three directions (regardless of the position of the receiving antenna) can be obtained(14)Hx=3M4πl2rsin3αcosαcos3β,Hy=−3M4πl2rsin2αcosαsinβcos2β,HzM4πl2rsin2αcos2β3cos3α−1.The relationship between the output voltage of the sensor and the relative position of the transmitting antenna and the receiving antenna (represented bya, β, l, and r) is as follows:
(15)Uv=CHx=3CM4πl2rsin3αcosαcos3β,where C is the signal gain coefficient of the receiving sensor, and in this case, C=150V/T; M is the equivalent magnetic moment considering the attenuation caused by the pipe thickness and overlying soil, taking M=30A·m2; l=5m, α∈[0, 180°], and β∈−90°,90°.The distribution and plane isogram ofUv with the change of relative position under horizontal transmission-vertical reception is shown in Figure 10. It can be seen from the figure that the distribution of Uv shows a double-peak characteristic, which is the same as the distribution feature on the X-o-Z plane obtained in the preceding. The UAV receives the electromagnetic field signal during the positioning process. According to the distribution regulation shown in the figure, UAV always travels along the gradient direction (the direction of the fastest descent or ascent of signal intensity) and explores the peak and trough positions shown in the figure. The midpoint of the peak and trough is the position directly above the robot. If the voltage is taken as an absolute value, the trough becomes a peak, and the robot is at the midpoint of the two peaks.Figure 10
Distribution and plane isogram ofUv under horizontal transmission-vertical reception.When the horizontal-receiving mode is adopted, the receiving antenna is parallel to the ground during the flight, but the course is not necessarily the same as the direction of the pipeline; that is, the axis of the receiving antenna and the axis of the transmitting antenna may produce a certain angleγ on the same horizontal plane. γ=0°, the receiving antenna axis, is parallel to the transmitting antenna axis, and γ=90°, the receiving antenna axis, is vertical to the transmitting antenna axis, as shown in Figure 11.Figure 11
Horizontal transmission-horizontal reception.The output voltage of the sensor under this working condition is(16)Uh=CHh=CHzcosγ−Hysinγ=3Cm4πl2rsin2αcos2β3cos3α−1cosγ−cosαsinβsinγ.The distribution and plane isogram ofUh with the change of relative position under horizontal transmission-horizontal reception is shown in Figure 12. When γ=0°, it can be seen that the distribution of Uh has single peak and double troughs and the absolute output voltage has a global maximum value with two local maximum values on both sides of it, which is the same as the analysis conclusion on the two-dimensional plane. When tracking and positioning, the UAV searches for the global maximum value along the steepest ascent curve of the output voltage amplitude to locate the pipeline robot. When γ increases gradually, the contour deflects and the positions of local maximum values rotate. The troughs representing the global maximum gradually change from single trough to connected trough, and finally forms double troughs. When γ=90°; that is, the axis of receiving antenna is perpendicular to the direction of the pipe; two global maximum values and two local maximum values are distributed at an equal interval of 45°. To sum up, the vertical reception should be used primarily in the positioning process by finding the midpoint of the peak and trough. Then, the horizontal reception is used as the supplement by searching for the maximum value of the sensor output signal to adjust the course of the UAV and assist the position verification.Figure 12
Distribution and plane isogram ofUh under horizontal transmission-horizontal reception.
(a)
γ=0°(b)
γ=25°(c)
γ=75°(d)
γ=90°


## 6.2. Inclined Pipeline Environment

In the practical long-distance pipeline-laying process, affected by rugged terrain, obstacles, and other factors, it is difficult for the pipeline to keep parallel with the horizontal plane for a long distance, and there are a lot of slopes and corners. Therefore, it is necessary to study the electromagnetic field distribution characteristics in the case of the inclined pipeline. As shown in Figure13, the angle between the pipeline and the horizontal plane (inclination angle) is θ, and θ∈0°~90° has the following relationship:
(17)r=lsinθcosα+cosθsinαcosβ.Figure 13
Inclined transmission-vertical reception.The components received in three directions can be obtained(18)Hx=3M4πl3sinαcosαcosβsinθcosα+cosθsinαcosβ3,Hy=3M4πl3sinαcosαsinβsinθcosα+cosθsinαcosβ3,HzM4πl33cos2α+1sinθcosα+cosθsinαcosβ3.Considering the vertical-receiving mode shown in Figure13, the vertical magnetic field intensity is as follows:
(19)Hv=Hxcosθ+Hzsinθ.The relationship between the output voltage of the sensor and the relative position of the pipeline-receiving antenna is as follows:(20)Uv=CHv=CM4πl3sinθ2cos2α−sin2α+3cosθsinαcosαcosβsinθcosα+cosθsinαcosβ3.The distribution and plane isogram ofUv with the change of relative position under inclined transmission-vertical reception is shown in Figure 14. It can be seen from the figure that when the gradient of pipeline is small, such as θ=10 degrees (see Figures 14(a) and 14(b)), the peak-trough characteristic is still apparent, and the pipeline robot can be located by finding the midpoint of the wave peak and trough; when the pipeline gradient increases gradually (see Figure 14(c)), the wave peak topples and falls to the side of the trough, and the trough rises gradually (corresponding amplitude gets gradually smaller); when θ>60° (see Figure 14(d)), the trough disappears and the peak appears only at the position where the receiver is far from the robot; and when θ>80° (Figure 14(e)), the wave peak and trough have disappeared, and it is not easy to locate the position of the robot. However, considering that the gradient of long-distance pipeline is generally not more than 45° in practice, ELF tracking and positioning method is still available.Figure 14
Distribution and plane isogram ofUv under inclined transmission-vertical reception.
(a)
θ=10°(b)
θ=20°(c)
θ=60°(d)
θ=80°When the receiving antenna is parallel to the ground, as shown in Figure15, the angle between the axis of the receiving antenna and the axis of the transmitting is still γ. The magnetic field intensity in the horizontal direction is
(21)Hh=Hzcosθcosγ−Hysinγ−Hxsinθcosγ.Figure 15
Inclined transmission-horizontal reception.The output voltage of the sensor under this working condition by substituting Formula (18) is
(22)Uh=CHh=CM4πl3cosθ2cos2α−sin2αcosγ−3cosθsinαcosαcosβcosγ+3sinαcosαcosβsinγsinθcosα+cosθsinαcosβ3.There are four variablesα, β, γ, and θ in the above formula, so the distribution of Uh is complicated. Figure 16 shows the distribution and plane isogram of Uh when θ=20°, 40°, and 60° under γ=0°, 30°, 60°, and 90°, respectively. It can be seen from the figure that when θ=20° (see Figure 16(a)), the distribution of the global maximum (trough) and local maximum (peak) is pronounced. With the increase of γ, the positions of the two peaks are deflected, and the trough is gradually split into two troughs, and finally evolve into double peaks and double troughs. The robot is at the midpoint of the peak and trough, which can be used for positioning. When θ=40° (see Figure 16(b)), there is only one global maximum value (single trough) under γ=0°. The position of the trough deflects to the rear of the robot, and a local maximum (peak) is generated, which increases with γ. The axial position of the pipe can be judged by the midpoint of the wave and trough. Although the precise position of the robot cannot be located, it can be considered whether it is in front or rear of the UAV. When θ=60° (see Figure 16(c)) and if γ is small, only the global maximum (single trough) can no longer locate the robot and the direction of the pipeline cannot be determined; if γ is further increased, a local maximum value (peak) appears, and the boundary line of the peak and valley is the pipe axis. Therefore, when the pipeline gradient is large, the horizontal-receiving antenna can hardly locate the pipeline robot. Still, it can help the UAV to return to the pipeline direction when it is yawing. At the same time, as mentioned above, the practical pipeline does not have such a large gradient.Figure 16
Distribution and plane isogram ofUh under inclined transmission-horizontal reception.
(a)
θ=20°, γ=0°, 30°, 60°, and 90°(b)
θ=40°, γ=0°, 30°, 60°, and 90°(c)
θ=60°, γ=0°, 30°, 60°, and 90°


## 7. Conclusion

The tracking and positioning of pipeline robot are closely related to the safe and reliable operation of pipeline cleaning and inspection equipment loaded by the robot. The tracking and positioning system using UAV presented in this paper is to improve the automation level and engineering applicability of pipeline robot and meet the urgent needs of trenchless detection and risk-elimination in pipeline engineering. The mathematical model of the ELF electromagnetic field distribution in the metal pipe environment established from classical electromagnetic theory provides a theoretical basis for the equivalent magnetic dipole model to describe the real system. The equivalent model is used to analyze the influence of electromagnetic parameters, geometrical dimensions, and transmitting frequency on the ELF electromagnetic field distribution characteristics in engineering pipeline environment. The investigation of ELF electromagnetic field distribution law under two conditions of laying pipeline horizontally and laying pipeline inclined is the prerequisite of realizing the accurate positioning of pipeline robot. The tracking method of moving robot inside the pipeline based on this research has the advantages of convenience, efficiency, and high accuracy, which not only ensures the safe operation of robot but also provides convenient conditions for pipeline construction operation such as excavation, maintenance, and robot’s removal when blocking accidents occur. The tracking and positioning method is mainly studied by numerical analysis in this paper, and the engineering application of the technology will be discussed in-depth in the future.



---

*Source: 4247277-2022-12-19.xml*