CN101982734A - Calculation method for underground magnetic navigation - Google Patents

Calculation method for underground magnetic navigation Download PDF

Info

Publication number
CN101982734A
CN101982734A CN 201010523432 CN201010523432A CN101982734A CN 101982734 A CN101982734 A CN 101982734A CN 201010523432 CN201010523432 CN 201010523432 CN 201010523432 A CN201010523432 A CN 201010523432A CN 101982734 A CN101982734 A CN 101982734A
Authority
CN
China
Prior art keywords
magnetic
signals
space
point
generated
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN 201010523432
Other languages
Chinese (zh)
Inventor
程强强
胡博
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanchang Hangkong University
Original Assignee
Nanchang Hangkong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanchang Hangkong University filed Critical Nanchang Hangkong University
Priority to CN 201010523432 priority Critical patent/CN101982734A/en
Publication of CN101982734A publication Critical patent/CN101982734A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)
  • Measuring Magnetic Variables (AREA)

Abstract

The invention discloses a calculation method for underground magnetic navigation, comprising the following steps: step 1, preprocessing acquired magnetic signals by an underground magnetic navigation method; step 2, establishing a model for calculating the strength of induced magnetic fields generated at any point in space by a magnetic source; step 3, starting from a magnetic dipole model; step 4, solving the undetermined coefficient in a relational expression shown in the specification by adopting the least square linear fitting method; step 5, adopting the model established in step 3 to solve the relation among an azimuth, a hole deviation angle and the strengths of magnetic inductions generated by the magnetic source at any point in space in the directions of x, y and z through back calculation; step 6, adopting three pairs of magnetic field strengths in the directions of x, y and z acquired by three fluxgate type sensors; and step 7, judging the position of the magnetic source relative to a magnetic vector sensor. The method has the advantages of adopted artificial magnetic beacon, long navigation distance and high angle measurement precision. The calculation method is simple, can be programmed conveniently by a computer, and has high operation efficiency.

Description

Calculation method for underground magnetic navigation
Technical Field
The invention relates to a navigation calculation method, in particular to a calculation method for underground magnetic navigation.
Background
The coal bed gas reserves in China are rich, the social demand and the growth are great, and the utilization prospect is wide. At present, the yield of the coal bed gas single well is low and the number of the wells is small. Current development attempts show that: in the future, a mode with a multi-branch horizontal well as a main well and a vertical well as an auxiliary well is formed. The multi-branch horizontal well technology in China basically depends on foreign companies, the domestic drilling technology is not matched, special tools and equipment are deficient, the scale development is severely restricted, the key problem is how to realize the communication between the main branch of the horizontal well and the vertical well, the underground magnetic navigation system is developed to solve the problem, the development yield of the oil and gas field and the coal bed gas can be greatly improved, and the exploitation efficiency is improved. The calculation method for underground magnetic navigation is that whether an underground navigation system can successfully navigate a drill collar or not, and the communication of two wells is important.
Disclosure of Invention
The invention aims to provide a calculation method for underground magnetic navigation, which can provide an effective high and new technical means for realizing the communication between a main branch of a horizontal well and a vertical well.
The invention is realized in the following way, which is characterized by comprising the following steps:
step 1, preprocessing the collected magnetic signals by an underground magnetic navigation method: the method comprises the steps of performing FIR digital filtering on original magnetic signals acquired by a magnetic vector sensor to filter geomagnetic signals and noise contained in the signals, wherein the noise signals are generally high-frequency signals, so that the signals are firstly passed through a low-pass filter to filter high-frequency noise, then performing FIR digital filtering on the signals, and filtering the geomagnetic signals to obtain magnetic signals generated by a permanent magnet pup joint;
step 2, establishing a magnetic dipole model on the basis of two basic laws, namely an ampere molecular circulation hypothesis and a biot-savart law, and establishing an induced magnetic field intensity calculation model generated by a magnetic source at any point in space;
step 3, deriving a relational expression of the distance and the magnetic induction intensity of the magnetic source in the x, y and z directions generated at any point in space from the magnetic dipole model;
step 4, solving a coefficient to be determined in the relational expression by adopting a least square linear fitting method;
step 5, inverting a relational expression of the azimuth angle, the inclination angle and the magnetic induction intensity of the magnetic source in the x, y and z directions generated at any point of the space by using the calculation model of the magnetic induction intensity generated by the magnetic source at any point of the space, which is established in the step 3, by adopting an extreme value taking method;
step 6, utilizing three pairs of magnetic field intensity of x, y and z acquired by the three fluxgate sensors, and applying a multi-sensor data fusion technology to further improve the angle measurement accuracy;
and 7, judging the position of the magnetic source relative to the magnetic vector sensor by using the three pairs of magnetic field strengths of x, y and z acquired by the three fluxgate sensors.
The invention has the advantages that: and the artificial magnetic beacon is adopted, so that the navigation distance is long, and the angle measurement precision is high. The calculation method is simple, the computer programming is convenient to realize, and the calculation efficiency is high.
Drawings
FIG. 1 is a graph of the induced magnetic field strength generated at any point in space by a magnetic source of the present invention.
FIG. 2 is a schematic diagram illustrating the determination of the direction of deviation of the magnetic source according to the present invention.
FIG. 3 is a basic schematic diagram of the underground magnetic navigation system of the present invention.
Detailed Description
The invention is realized by the following steps:
step 1, preprocessing the collected magnetic signals by an underground magnetic navigation method: the method comprises the steps of performing FIR digital filtering on original magnetic signals acquired by a magnetic vector sensor to filter geomagnetic signals and noise contained in the signals, wherein the noise signals are generally high-frequency signals, so that the signals are firstly passed through a low-pass filter to filter high-frequency noise, then performing FIR digital filtering on the signals, and filtering the geomagnetic signals to obtain magnetic signals generated by a permanent magnet pup joint;
step 2, establishing a magnetic dipole model on the basis of two basic laws, namely an ampere molecular circulation hypothesis and a biot-savart law, and establishing an induced magnetic field intensity calculation model generated by a magnetic source at any point in space; as shown in figure 1, when the induced magnetic field intensity generated by a magnetic source at any point in space is calculated, the permanent magnet pup joint is regarded as a pair of magnetic dipoles, and three directional components of the magnetic field intensity x, y and z generated by the pair of magnetic dipoles at any point in space are calculated
Figure DEST_PATH_IMAGE001
Figure 695446DEST_PATH_IMAGE002
Figure DEST_PATH_IMAGE003
Figure 674904DEST_PATH_IMAGE004
Wherein,
Figure 923482DEST_PATH_IMAGE006
Figure 912167DEST_PATH_IMAGE008
Figure DEST_PATH_IMAGE010
Figure DEST_PATH_IMAGE012
the angular frequency of the rotation of the permanent magnet short section is shown, and x, y and z are coordinates from the center of the permanent magnet short section to any point in space in three directions.
Step 3, fromStarting from a magnetic dipole model, deducing a relational expression of the distance and the magnetic induction intensity of the magnetic source in the x, y and z directions generated at any point in space; the distance between the center of the magnetic source and the magnetic vector sensor and the one-third root of the magnetic field intensity in the y direction have a linear relation, namely
Figure DEST_PATH_IMAGE014
Figure 352638DEST_PATH_IMAGE014
Step 4, solving a coefficient to be determined in the relational expression by adopting a least square linear fitting method; the coefficients to be determined, a =38.28 and B =0.2192, in the relation can be obtained by fitting a test data curve by a least square method, so that the distance relation between the magnetic source center and the magnetic vector sensor is:
Figure DEST_PATH_IMAGE016
step 5, inverting a relational expression of the azimuth angle, the inclination angle and the magnetic induction intensity of the magnetic source in the x, y and z directions generated at any point of the space by using the calculation model of the magnetic induction intensity generated by the magnetic source at any point of the space, which is established in the step 3, by adopting an extreme value taking method; the relation between the azimuth angle and the magnetic induction intensity of the magnetic source in the x, y and z directions generated at any point in the space is inverted by adopting an extreme value taking method, and the relation is as follows:
Figure DEST_PATH_IMAGE018
wherein
Figure DEST_PATH_IMAGE020
Figure DEST_PATH_IMAGE022
The x-direction magnetic induction and the z-direction magnetic induction generated by the magnetic source collected by the magnetic vector sensor of claim 3.
The relation between the inclination angle and the magnetic induction intensity of the magnetic source in the x, y and z directions generated at any point in space can be known from the symmetry as follows:
Figure DEST_PATH_IMAGE024
wherein
Figure 768576DEST_PATH_IMAGE020
The x-direction magnetic induction and the y-direction magnetic induction generated by the magnetic source collected by the magnetic vector sensor of claim 3.
Step 6, utilizing three pairs of magnetic field intensity of x, y and z acquired by the three fluxgate sensors, and applying a multi-sensor data fusion technology to further improve the angle measurement accuracy;
and 7, judging the position of the magnetic source relative to the magnetic vector sensor by utilizing the magnetic field intensity of three pairs of x, y and z acquired by the three fluxgate sensors, wherein as shown in fig. 2, the three fluxgate sensors in the magnetic vector sensor are arranged on the same axis at an angle different from 7.5 degrees, when the sensors are aligned with the magnetic source, one fluxgate sensor is aligned with the magnetic source, the other fluxgate sensor is shifted by 7.5 degrees leftwards, and the other fluxgate sensor is shifted by 7.5 degrees rightwards. The magnetic field intensity generated by the arbitrary point p in the y direction of the rectangular coordinate system where the three sensors are located is
Figure DEST_PATH_IMAGE028
Therefore, in a period of data collected by the three sensors, which sensor has the minimum magnetic induction intensity indicates that the drill bit is deviated to which sensor, so that the direction of deviation of the drill bit from the normal track can be judged.
The distance and the relative position between the drill bit and the magnetic vector sensor can be calculated through the steps, and then the distance and the relative position are displayed through a computer, so that people on the ground can know the position of the drill bit in the horizontal well at any time.
As shown in fig. 3, the drill bit 1 is fixedly connected to the permanent magnet sub 2, the magnetic vector sensor 4 receives the magnetic signal of the permanent magnet sub 2, and the magnetic signal acquired by the magnetic vector sensor 4 = the geomagnetic signal + the signal generated by the magnetic source + the noise, so that the acquired magnetic signal is digitally filtered to remove the geomagnetic signal and the noise to obtain the sub signal generated by the magnetic source. And (3) obtaining the distance and the angle of the magnetic source relative to the magnetic vector sensor by adopting the distance and angle calculation formula derived in the step (5) and the step (6) for the preprocessed signals, obtaining a more accurate angle value by using the multi-sensor data fusion technology in the step (7), and finally transmitting the obtained distance and angle to a computer on the ground, wherein the computer analyzes the navigation track 3 for directional drilling of the navigation drill collar.

Claims (1)

1. A calculation method for underground magnetic navigation is characterized by comprising the following steps:
step 1, preprocessing the collected magnetic signals by an underground magnetic navigation method: the method comprises the steps of performing FIR digital filtering on original magnetic signals acquired by a magnetic vector sensor to filter geomagnetic signals and noise contained in the signals, wherein the noise signals are generally high-frequency signals, so that the signals are firstly passed through a low-pass filter to filter high-frequency noise, then performing FIR digital filtering on the signals, and filtering the geomagnetic signals to obtain magnetic signals generated by a permanent magnet pup joint;
step 2, establishing a magnetic dipole model on the basis of two basic laws, namely an ampere molecular circulation hypothesis and a biot-savart law, and establishing an induced magnetic field intensity calculation model generated by a magnetic source at any point in space;
step 3, deriving a relational expression of the distance and the magnetic induction intensity of the magnetic source in the x, y and z directions generated at any point in space from the magnetic dipole model;
step 4, solving a coefficient to be determined in the relational expression by adopting a least square linear fitting method;
step 5, inverting a relational expression of the azimuth angle, the inclination angle and the magnetic induction intensity of the magnetic source in the x, y and z directions generated at any point of the space by using the calculation model of the magnetic induction intensity generated by the magnetic source at any point of the space, which is established in the step 3, by adopting an extreme value taking method;
step 6, utilizing three pairs of magnetic field intensity of x, y and z acquired by the three fluxgate sensors, and applying a multi-sensor data fusion technology to further improve the angle measurement accuracy;
and 7, judging the position of the magnetic source relative to the magnetic vector sensor by using the three pairs of magnetic field strengths of x, y and z acquired by the three fluxgate sensors.
CN 201010523432 2010-10-29 2010-10-29 Calculation method for underground magnetic navigation Pending CN101982734A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010523432 CN101982734A (en) 2010-10-29 2010-10-29 Calculation method for underground magnetic navigation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010523432 CN101982734A (en) 2010-10-29 2010-10-29 Calculation method for underground magnetic navigation

Publications (1)

Publication Number Publication Date
CN101982734A true CN101982734A (en) 2011-03-02

Family

ID=43619636

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010523432 Pending CN101982734A (en) 2010-10-29 2010-10-29 Calculation method for underground magnetic navigation

Country Status (1)

Country Link
CN (1) CN101982734A (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102998708A (en) * 2011-05-25 2013-03-27 利宇Tec株式会社 Apparatus and method for detecting location of underground facility
CN104806166A (en) * 2015-03-25 2015-07-29 宁夏煤炭勘察工程公司 Butt joint method for ground multi-branch horizontal well and underground horizontal drilling far end
CN105928511A (en) * 2016-04-18 2016-09-07 哈尔滨工业大学 Positioning and orientation method based on double magnetic beacons
CZ307060B6 (en) * 2016-10-31 2017-12-20 České vysoké učení technické v Praze - Fakulta elektrotechnická A vehicle presence detector
CN113124882A (en) * 2021-06-17 2021-07-16 天津大学 Multi-dipole magnetic source inversion positioning method under condition of unknown background magnetic field
US20220120169A1 (en) * 2020-10-16 2022-04-21 Halliburton Energy Services, Inc. Use of residual gravitational signal to perform anomaly detection
CN114749700A (en) * 2021-12-10 2022-07-15 盐城工学院 Magnetic auxiliary control method for straightness error of non-magnetic material gun drill processing
CN115963038A (en) * 2022-12-14 2023-04-14 中国科学院空间应用工程与技术中心 Magnetic particle motion track measuring system and method based on space microgravity condition

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001046554A1 (en) * 1999-12-21 2001-06-28 Utilx Corporation Long range electronic guidance system for locating a discrete in-ground boring device
US20030111268A1 (en) * 1999-09-24 2003-06-19 Vermeer Manufacturing Company Underground boring machine employing navigation sensor and adjustable steering
US20050140373A1 (en) * 2003-05-22 2005-06-30 Schlumberger Technology Corporation Directional electromagnetic wave resistivity apparatus and method
CN201095981Y (en) * 2007-10-10 2008-08-06 新疆石油管理局钻井工艺研究院 Permanent magnetism short joint for well drilling
CN101361660A (en) * 2008-05-16 2009-02-11 深圳先进技术研究院 Multi-magnetic target positioning method and system
US20100234685A1 (en) * 2009-03-16 2010-09-16 Aleksandar Juloski Coil system for the contact-free magnetic navigation of a magnetic body in a working space

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030111268A1 (en) * 1999-09-24 2003-06-19 Vermeer Manufacturing Company Underground boring machine employing navigation sensor and adjustable steering
WO2001046554A1 (en) * 1999-12-21 2001-06-28 Utilx Corporation Long range electronic guidance system for locating a discrete in-ground boring device
US20050140373A1 (en) * 2003-05-22 2005-06-30 Schlumberger Technology Corporation Directional electromagnetic wave resistivity apparatus and method
CN201095981Y (en) * 2007-10-10 2008-08-06 新疆石油管理局钻井工艺研究院 Permanent magnetism short joint for well drilling
CN101361660A (en) * 2008-05-16 2009-02-11 深圳先进技术研究院 Multi-magnetic target positioning method and system
US20100234685A1 (en) * 2009-03-16 2010-09-16 Aleksandar Juloski Coil system for the contact-free magnetic navigation of a magnetic body in a working space

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《中国优秀硕士学位论文全文数据库(电子期刊)》 20061231 杨翡 定向钻井导航系统 全文 1 , 第3期 2 *
《中国博士学位论文全文数据库(电子期刊)》 20071231 徐涛 水平定向钻进随钻测量方法及定位技术研究 全文 1 , 第5期 2 *
《仪器仪表学报》 20090930 徐涛等 水平定向钻进随钻测量系统研究与设计* 第1976-1980页 1 第30卷, 第9期 2 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102998708A (en) * 2011-05-25 2013-03-27 利宇Tec株式会社 Apparatus and method for detecting location of underground facility
CN102998708B (en) * 2011-05-25 2016-01-20 利宇Tec株式会社 For detecting equipment and the method for location of underground facility
CN104806166A (en) * 2015-03-25 2015-07-29 宁夏煤炭勘察工程公司 Butt joint method for ground multi-branch horizontal well and underground horizontal drilling far end
CN105928511A (en) * 2016-04-18 2016-09-07 哈尔滨工业大学 Positioning and orientation method based on double magnetic beacons
CN105928511B (en) * 2016-04-18 2018-11-23 哈尔滨工业大学 A kind of positioning and orienting method based on double magnetic beacons
CZ307060B6 (en) * 2016-10-31 2017-12-20 České vysoké učení technické v Praze - Fakulta elektrotechnická A vehicle presence detector
US20220120169A1 (en) * 2020-10-16 2022-04-21 Halliburton Energy Services, Inc. Use of residual gravitational signal to perform anomaly detection
CN113124882A (en) * 2021-06-17 2021-07-16 天津大学 Multi-dipole magnetic source inversion positioning method under condition of unknown background magnetic field
CN113124882B (en) * 2021-06-17 2021-09-28 天津大学 Multi-magnetic dipole magnetic source inversion positioning method under condition of unknown background magnetic field
CN114749700A (en) * 2021-12-10 2022-07-15 盐城工学院 Magnetic auxiliary control method for straightness error of non-magnetic material gun drill processing
CN115963038A (en) * 2022-12-14 2023-04-14 中国科学院空间应用工程与技术中心 Magnetic particle motion track measuring system and method based on space microgravity condition
CN115963038B (en) * 2022-12-14 2023-07-28 中国科学院空间应用工程与技术中心 Magnetic particle motion trail measurement system and method based on space microgravity condition

Similar Documents

Publication Publication Date Title
CN101982734A (en) Calculation method for underground magnetic navigation
US10161246B1 (en) Rotatable sensors for measuring characteristics of subterranean formation
CN101798918B (en) Method for determining relative spatial position of adjacent well parallel segment
CN104131808B (en) A kind of device based on transient electromagnetic method coordinate detection accident well
CN101852078B (en) Electromagnetic distance measurement guide system for double solenoid set during drilling
CN104343438B (en) Measure the rotating excitation field rangefinder and its measurement method of drilling well relative distance
CN101713285B (en) Calculation method for measuring distance between adjacent wells by electromagnetic detection while drilling
CN101799558B (en) Electromagnetic surveying system while drilling of adjacent-well parallel intervals
AU2008248145B2 (en) Distance determination from a magnetically patterned target well
CN108166971B (en) A kind of Multifunctional inclination measuring apparatus and the method based on inclinometer progress offset well casing relative position judgement
CN101806211B (en) Calculation method using solenoid groups to achieve electromagnetic guiding distance measurement while drilling
CN102052069A (en) Near-bit measurement while drilling (MWD) system and method
CN104594881B (en) Method for determining relative special positions of adjacent well parallel sections
CN105928511B (en) A kind of positioning and orienting method based on double magnetic beacons
CN212083694U (en) Downhole time domain electromagnetic multi-component far detection instrument and detection system based on same
CN108873073B (en) A kind of across hole resistivity tomography method of three-dimensional based on electrical method of network concurrency
CN113687428B (en) Precise electromagnetic calibration method for ultra-deep underground pipeline azimuth
CN106291745B (en) Granite type U-ore target area preferred method based on Geophysical-chemical exception
CN103995286B (en) A kind of acoustic receiver transducer and orientation, down-hole noise leakage detection apparatus and leak hunting method
CN104234697A (en) Intelligent drilling angular difference measurement device and angular different measurement method
CN106014385B (en) A kind of guidance method of no-dig technique guide instrument
CN214170529U (en) Multi-sensor fusion attitude angle resolving circuit
Tian et al. A high-precision magnetic induction through-the-earth positioning scheme
CN106197237B (en) A kind of judgment method of the measurement point of no-dig technique guide instrument
CN117090558B (en) Rescue well track adjusting method and device

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent of invention or patent application
CB03 Change of inventor or designer information

Inventor after: Yu Runqiao

Inventor after: Cheng Qiangqiang

Inventor after: Hu Bo

Inventor after: Ying Zheng

Inventor after: Dai Pupu

Inventor before: Cheng Qiangqiang

Inventor before: Hu Bo

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: CHENG QIANGQIANG HU BO TO: YU RUNQIAO CHENG QIANGQIANG HU BO YING ZHENG DAI PUPU

C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20110302