WO2015094015A1 - Устройство и способ для определения положения трубопровода - Google Patents

Устройство и способ для определения положения трубопровода Download PDF

Info

Publication number
WO2015094015A1
WO2015094015A1 PCT/RU2014/000859 RU2014000859W WO2015094015A1 WO 2015094015 A1 WO2015094015 A1 WO 2015094015A1 RU 2014000859 W RU2014000859 W RU 2014000859W WO 2015094015 A1 WO2015094015 A1 WO 2015094015A1
Authority
WO
WIPO (PCT)
Prior art keywords
pipeline
determining
hardware
navigation system
inspection device
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.)
Ceased
Application number
PCT/RU2014/000859
Other languages
English (en)
French (fr)
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.)
Transneft PJSC
Transneft Diascan JSC
Original Assignee
Transneft Diascan JSC
Oil Transporting JSC Transneft
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 Transneft Diascan JSC, Oil Transporting JSC Transneft filed Critical Transneft Diascan JSC
Priority to MX2015016216A priority Critical patent/MX2015016216A/es
Priority to BR112015029372A priority patent/BR112015029372A2/pt
Priority to EA201591997A priority patent/EA033119B1/ru
Priority to EP14870724.3A priority patent/EP2985509A4/en
Publication of WO2015094015A1 publication Critical patent/WO2015094015A1/ru
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17DPIPE-LINE SYSTEMS; PIPE-LINES
    • F17D5/00Protection or supervision of installations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • G01C21/165Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments

Definitions

  • a device for measuring the spatial position of pipelines (patent No. 1809297 A1 (RU) of April 15, 1993) is known, which contains a platform with horizontal and vertical tilt angle sensors placed on it, and a distance sensor equipped with a body to increase the accuracy and performance of measurements. On the case there are elastic cuffs for basing the device in the pipeline.
  • the tilt angle sensor is designed as a vertical vertical
  • the tilt sensor in the horizontal plane is designed as a gyro semi-compass.
  • the device is equipped with a system of secondary converters, a switch, an analog-to-digital converter and a microcontroller.
  • accelerometer associated with the orthogonal X axis of the SINS, ah [m / s 2 ];
  • DKOPT consists of hardware and software.
  • AC includes accelerometers 1, 2, 3 (Fig. 1), gyroscopes 4, 5, 6 (Fig. 1) in the associated SINS and, respectively, with the VIP, rectangular coordinate system and odometer 7 (Fig. 1 ).
  • the data obtained from the ACP DKOPT is an array containing the readings of three accelerometers 1, 2, 3 (Fig. 1) and three gyros 4, 5, 6 (Fig. 1) associated with the orthogonal axes of the SINS and, accordingly, with the VIP rectangular coordinate system and odometer reading 7 (Fig. 1).
  • c is the deleted amount of ranked data in fractions of a unit.
  • the “Exhibition” algorithm that is, the determination of the position of the VIP in space 11 (Fig. 2): at the initial moment of time of the VIP pass, it remains immobile for several minutes. This time interval is tracked by the values of the speed of movement obtained in the previous stage. At this time, an exhibition BINS.
  • the initial angles of azimuth and pitch in a rectangular system of coordinates North-East-Height are determined by the rotation of the Earth and given geodesic coordinates of the VIP. In the future, when calculating the VIP passage through the pipeline, the determined values of these angles make it possible to take into account the rotation of the Earth and subtract the corresponding component of the angular velocity of the Earth from the gyroscope readings. Thus, the rotation of the Earth has practically no effect on the navigation parameters to be determined.
  • Exhibition mode is divided into two stages: rough and accurate exhibition.
  • a rough exhibition is an approximate calculation of parameters:
  • t and 2 are, respectively, the vertical coordinates of the beginning and end of the movement section of the inspected section of the pipeline at the base pass of DKOPT and Z2, respectively, the vertical coordinates of the beginning and end of the section of movement of the inspected section of the pipeline during the subsequent passage of DKOPT.
  • the direction of the calculated displacement vector 27 (Fig. 7) must be determined in the plane 26 (Fig. 7), perpendicular to the axis of the pipeline, the position “at 12 o'clock” is taken as 0 °. Inside the area of movement of the area of the inspected pipeline are calculated direction vectors. The points of the beginning and end of the vector are invited to establish the coordinates of the welds of individual sections of the pipeline.
  • the calculated displacement vector 27 (Fig. 7):
  • A (Ax p , Ay p , Az p ).
  • K - P A Ax p n ⁇ x + Ay p ⁇ ⁇ ⁇ + ⁇ ⁇ ⁇ ⁇ ⁇ . .
  • Ar ⁇ - K ⁇ I - (Ax - K ⁇ f x, Ay - K ⁇ ⁇ ⁇ , ⁇ p - K ⁇ ⁇ ⁇ ),

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Navigation (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

Изобретение относится к измерительной технике и может быть использовано для определения положения трубопровода в пространстве. Устройство состоит из аппаратной части и программной части, при этом в аппаратной части каждый из трех гироскопов и акселерометров бесплатформенной инерционной навигационной системой связан с одной из ортогональных осей внутритрубного инспекционного прибора. Программная часть состоит из алгоритмов определения навигационных параметров с использованием глобальной спутниковой навигационной системы GPS или/и GLONASS. Данные переписываются с аппаратной части, и расчет навигационных параметров производится на ЭВМ с помощью специальной программы. Способ определения относительного перемещения трубопровода состоит в том, что используются данные измерения величин и направлений линейных перемещений участка трубопровода по результатам двух и более инспекционных пропусков внутритрубного инспекционного прибора с установленным на нем аппаратной части по одному и тому же участку трубопровода. Полученные данные располагаются по таблицам и совмещаются по дистанциям. Критерием наличия перемещения трубопровода на инспектируемом участке является превышение модуля разностной кривизны заданного порогового значения.

Description

УСТРОЙСТВО И СПОСОБ ДЛЯ ОПРЕДЕЛЕНИЯ
ПОЛОЖЕНИЯ ТРУБОПРОВОДА
Описание изобретения
Изобретение относится к измерительной технике и может быть использовано для определения положения трубопровода в пространстве и в том числе в горизонтальной и вертикальной плоскостях при эксплуатации и строительстве трубопроводов.
Известно устройство для измерения пространственного положения трубопроводов (патент ΜΊ_1809297 А1 (RU) от 15.04.1993), содержащее платформу с размещенными на ней датчиками угла наклона в горизонтальной и вертикальной плоскостях, датчик пройденного расстояния, оснащенное с целью повышения точности и производительности измерений корпусом. На корпусе установлены эластичные манжеты для базирования устройства в трубопроводе. Датчик угла наклона выполнен в виде гировертикали, а датчик наклона в горизонтальной плоскости - в виде гирополукомпаса. Устройство снабжено системой вторичных преобразователей, коммутатором, аналогово-цифровым преобразователем и микроконтроллером.
Известно устройство для определения пространственного положения магистральных трубопроводов (патент N° 1404815 А1 (SU) от 23.06.1988). Устройство содержит корпус с гиросистемой, состоящей из гирокомпаса, гировертикали и датчика положения корпуса, выполненный в виде излучателя, скрепленного с гиросистемой. Корпус предназначен для установки в трубопроводе. Гиросистема снабжена фоконной матрицей, торцы которой оптически сопряжены с излучателем, со светочувтвительной лентой, обладающей возможностью перемещения в корпусе.
Известно устройство контроля внутренней поверхности, пространственного положения и напряженного состояния трубопровода (патент N°2106569 CI (RU) от 10.03.1998). Устройство имеет измерительный модуль, состоящий из контейнера, внутри которого расположена аппаратура измерения плановых и высотных меток, состоящей из цифрового вычислительного комплекса и трехосного гиростабилизатора с гироблоком и акселерометром. На наружной поверхности контейнера расположен второй датчик пути, соединенный с аппаратурой измерения плановых и высотных меток.
Недостатками выше указанных устройств являются автономность полученных данных, что влечет за собой погрешности определения положения трубопровода; а также полученные навигационные данные, не отражают всех параметров положения трубопровода, таких как углы и радиусы изгиба трубопровода в пространстве.
Известна аппаратура для измерения линейных деформаций магистрального трубопровода (патент N°2334162 С1 от 20.09.2008), содержащая струнные датчики деформации, установленные в различных сечениях магистрального трубопровода с определенным одинаковым шагом. Аппаратура дополнительно содержит приемник спутниковой радионавигационной системы, мультиплексор, микропроцессор, запоминающее устройство начального уровня, аналогово-цифровой преобразователь. Недостатком данного изобретения является стационарная установка датчиков деформации, что не позволяет получить полную информацию об изменении положения трубопровода.
Цель изобретения состоит в создании устройства диагностического комплекса для определения положения трубопровода (далее ДКОПТ) с использованием глобальной навигационной спутниковой систему GPS или/и GLONASS, которое повышает точность измерений пространственного положения трубопровода с учетом корректирующих точек, а также повышает точность измерений углов изгиба трубопровода - тангажа и азимута, и радиусов изгиба в каждой точке трубопровода с учетом вращения Земли; а также в способе определения относительного перемещения трубопровода по результатам двух и более инспекционных пропусков ДКОПТ, наглядно позволяющем определить перемещение трубопровода в пространстве путем сравнения более поздних траекторий трубопровода с базовой.
Технический результат достигается тем, что заявленное устройство ДКОПТ состоит из аппаратной части (далее АЧ) и программной части. Аппаратная часть ДКОПТ установлена на внутритрубный инспекционный прибор (далее ВИП). Аппаратная часть состоит из набора датчиков: акселерометр, связанный с осью X бесплатформенной инерционной навигационной системы (далее БИНС), акселерометр, связанный с осью Y БИНС, акселерометр, связанный с осью Z БИНС, гироскоп, связанный с осью X БИНС, гироскоп, связанный с осью Y БИНС, гироскоп, связанный с осью Z БИНС, одометр. Программная часть состоит из алгоритмов определения навигационных параметров в следующей последовательности: расчет линейной скорости при этом определение движения внутритрубного инспекционного прибора происходит по данным, полученным от аппаратной части устройства диагностического комплекса для определения положения трубопровода; выставка, то есть определение положения в пространстве ВИП производится определением начальных углов азимута и тангажа в прямоугольной системе координат Север- Восток-Высота с учетом вращения Земли и заданным геодезическим координатам ВИП в течение первых нескольких минут, когда ВИП неподвижен; расчет навигационных параметров происходит с использованием значение дискрета времени между отчетами акселерометров и гироскопов и показаниями одометра с определением соответствующих углов: крена, тангажа, азимута; расчет радиусов изгиба трубопровода как в вертикальном, так и горизонтальном направлении; коррекция траектории происходит при использовании координат корректирующих точек, полученных привязкой конкретных точек трубопровода к геодезическим координатам в системе Север-Восток-Высота с использованием глобальной спутниковой навигационной системы GPS или/и GLONASS. Данные переписываются с АЧ ДКОПТ, и расчет навигационных параметров производится на ЭВМ с помощью специальной программы. Таким образом, данные, полученные с АЧ ДКОПТ представляют собой массив, содержащий показания трех акселерометров и трех гироскопов в связанной с ортогональными осями БИНС и соответственно с ВИП прямоугольной системой координат и показания одометра. При этом ось у БИНС соответствует продольной оси ВИП. После инспекционного пропуска по трубопроводу ВИП массив данных переписывается с АЧ ДКОПТ, и расчет всех навигационных параметров производится на ЭВМ с помощью специальной программы. Используя данные ДКОПТ с использованием глобальной навигационной спутниковой системы GPS или/и GLONASS, повышается точность измерений пространственного положения трубопровода с учетом корректирующих точек. Величина и направления перемещения трубопровода определяются по результатам двух и более инспекционных пропусков по одному и тому же участку трубопровода. Способ определения относительного перемещения трубопровода по результатам двух и более инспекционных пропусков ДКОПТ состоит в том, что используются данные измерения величин и направлений линейных перемещений участка трубопровода по результатам двух и более инспекционных пропусков ВИП с установленным на нем ДКОПТ по одному и тому же участку трубопровода с одним и тем же направлением потока рабочей среды. Данные, полученные с ДКОПТ, располагаются по таблицам и совмещаются по дистанциям. Используются данные, полученные с ДКОПТ, для вычисления максимальной величины, направления и угла перемещения трубопровода.
Данные пропуска ВИП с установленной на нем АЧ ДКОПТ с более ранней датой считаются базовыми, а данные последующих пропусков сравниваются с базовыми. Критерием наличия перемещения трубопровода на инспектируемом участке является превышение модуля разностной кривизны заданного порогового значения.
Предложенный способ определения относительного перемещения трубопровода по результатам двух и более инспекционных пропусков ДКОПТ наглядно позволяет определить перемещение трубопровода в пространстве
На фиг. 1 изображена аппаратная часть ДКОПТ, включающая в себя набор датчиков.
1. акселерометр, связанный с ортогональной осью X БИНС, ах [м/ С2];
2. акселерометр, связанный с ортогональной осью Y БИНС, ау [м/с2];
3. акселерометр, связанный с ортогональной осью Z БИНС, az м/ С2];
Figure imgf000006_0001
На фиг. 2 изображена блок-схема алгоритмов определения навигационных параметров программной части ДКОПТ: 8. «Начало»;
9. «Получение данных (по акселерометрам, гироскопам и одометру) для значений сканов j=l ...N;
10. алгоритм «Расчет линейной скорости»;
11. алгоритм «Выставка», то есть определение положение ВИП в пространстве;
12. алгоритм «Расчет навигационных параметров»;
13. алгоритм «Расчет радиусов»;
14. алгоритм «Коррекция траектории»;
15. «Вывод данных»;
16. «Конец».
На фиг. 3 изображена траектория залегания трубопровода, выявленная по результатам базового пропуска ДКОПТ:
17. точка начала участка перемещения инспектируемого участка трубопровода с координатами xi,yi,zi;
18. точка конца участка перемещения инспектируемого участка трубопровода с координатами x2,y2,z2;
На фиг. 4 изображена траектория залегания трубопровода, выявленная по результатам повторного пропуска ДКОПТ:
19. точка начала участка перемещения инспектируемого участка трубопровода с координатами Х±\, Ζ^,
20. точка конца участка перемещения инспектируемого участка трубопровода с координатами 2, у2 , 22.
На фиг. 5 изображены выделенные горизонтальная и вертикальная составляющие кривизны трубопровода, выявленные по результатам базового пропуска ДКОПТ:
21. вертикальная составляющая кривизны трубопровода, EV;
22. горизонтальная составляющая кривизны трубопровода, ЕН.
На фиг. 6 изображены выделенные горизонтальная и вертикальная составляющие кривизны трубопровода, выявленные по результатам повторного пропуска ДКОПТ: б
23. вертикальная составляющая кривизны трубопровода, EV ;
24. горизонтальная составляющая кривизны трубопровода, ЕН .
На фиг. 7 изображено определение направления вектора перемещения трубопровода и угол направления перемещения трубопровода:
25. единичный вектор направления оси трубопровода, п;
26. плоскость, перпендикулярная единичному вектору П, р;
27. рассчитанный вектор перемещения оси трубопровода, Δ ;
28. проекция рассчитанного вектора перемещения оси трубопровода Δ на плоскость р, Лр;
29. единичный вектор вертикали, Z;
30. проекция вертикали z на плоскость р, ;
31. искомый угол направления перемещения трубопровода, <Х.
ДКОПТ состоит из аппаратной части и программную часть. АЧ (фиг. 1) включает в себя акселерометры 1, 2, 3 (фиг. 1), гироскопы 4, 5, 6 (фиг. 1) в связанной с БИНС и соответственно с ВИП прямоугольной системе координат и одометра 7 (фиг. 1). Данные, полученные с АЧ ДКОПТ представляют собой массив, содержащий показания трех акселерометров 1 , 2, 3 (фиг. 1) и трех гироскопов 4, 5, 6 (фиг. 1) в связанной с ортогональными осями БИНС и соответственно с ВИП прямоугольной системе координат и показания одометра 7 (фиг. 1). Программная часть состоит из алгоритмов определения навигационных параметров на основе данных с АЧ ДКОПТ для значения сканов j=l ...N 9 (фиг. 2) в следующей последовательности:
Алгоритм «Расчет линейной скорости» 10 (фиг. 2): используя показания одометра и показания таймера импульсов одометра, рассчитывается скорость движения ДКОПТ:
Sj=bjfo,
V=F{S, L,.C)/ М, где
fo - одометрический фактор;
Sj - пройденный путь; S=[Sj S2 ... Stf], V^^j Έ?2 ... νΜ] - линейная скорость;
At - дискрет времени;
F{»,LA,c} - фильтрация скользящих средним с прямоугольным окном: Pf=
среднее по выборке проранжированных по
Figure imgf000009_0001
возрастанию данных Рг, исходными для которых является массив
Figure imgf000009_0002
Для.) -отсчета;
j-L>0, j+L<N, L - половина ширины окна фильтра;
с - удаляемое количество проранжированных данных в долях единицы.
Алгоритм «Выставка», то есть определение положения ВИП в пространстве 11 (фиг. 2): в начальный момент времени пропуска ВИП остается неподвижен несколько минут. Этот промежуток времени отслеживается по значениям скорости движения, полученной на предыдущем этапе. В это время производится выставка БИНС. По вращению Земли и заданным геодезическим координатам ВИП определяются начальные углы азимута и тангажа в прямоугольной системе координат Север- Восток-Высота. В дальнейшем при расчете прохождения ВИП по трубопроводу определяемые значения этих углов позволяют учесть вращение Земли и вычесть соответствующую компоненту угловой скорости Земли из показаний гироскопов. Таким образом, вращение Земли практически не оказывает влияния на определяемые навигационные параметры. Режим выставки подразделяется на два этапа: грубая и точная выставки.
Грубая выставка- приближенное вычисление параметров:
Figure imgf000009_0003
Где arctg(») - вычисление арктангенса с учетом квадранта аргументов; α, =— 2ΐαι - > ^'У'2' - количество отсчетов, отводимых на грубую
4 y=i
выставку, во время которой ВИП неподвижен;
ус - угол крена;
0с - угол тангажа;
Точная выставка- уточнение параметров проводится в несколько этапов:
Первый этап:
cos(0c/2) sin(0c /2) 0 0 cos( c/2)'
- sin(0c/2) cos(#c/2) 0 0 0
0 0 cos(0c/2) sin(#c /2) sin(rC/2) 0 0 - sin(0c/2) cos(<9c/2) 0
Figure imgf000010_0002
Figure imgf000010_0001
KLc - начальное значение нормированного полного кватерниона, полученного из углов Эйлера;
QL - начальное значение угловой скорости Земли;
j=Lc...La, La - количество отсчетов, отводимых на точную выставку, в это время ВИП неподвижен;
ί=ω, g, >A - широта точки выставки;
Yj, θβ ψ Xj - соответственно текущие углы крена, тангажа, азимута, угол
направления на географический север - азимут навигационной системы координат;
СОЕ - угловая скорость Земли;
С\ ...с6 - константы.
Второй этап:
Figure imgf000011_0001
текущее значение нормированного полного кватерниона.
Третий этап
2k2jkAj—kljki-j)
Yj = -arctg
Четвертый этап:
Figure imgf000011_0002
Пятый этап ψ] =
Figure imgf000011_0003
Далее идет расчет текущего значения угловой скорости:
Figure imgf000011_0004
Figure imgf000011_0006
Figure imgf000011_0005
-Ale,,
1-е
1 = e'^I + Zj = arctg, {ry ./rxj ),
Figure imgf000012_0001
Bj - приращение линейной скорости на текущем кванте времени;
ltJ и г,. . - параметры петли обратной связи по оценке угловой скорости вращения прибора в горизонтальной плоскости;
- символ принадлежности к связанной системе координат.
Навигационная система координат (далее НСК) представляет собой прямоугольную систему координат, начало которой связано с точкой выставки, основная плоскость - плоскость местного горизонта, одна ось направлена вверх, При преобразовании из НСК в связанную систему и наоборот в результирующий вектор входят последние три элемента перемножения.
Алгоритм «Расчет навигационных параметров» 12 (фиг. 2) делится на несколько этапов:
Первый этап - расчет промежуточных параметров:
Figure imgf000012_0002
' cos(Ae/2) 0 0 sin(la/2)'
Figure imgf000013_0001
0 -sin(A0/2) cos( a/2) 0
-sin Aa/2) 0 0 cos(la/2)
Figure imgf000013_0002
Второй этап
По математическим выражениям первого и второго этапов «Точной выставки» алгоритма «Выставка» определяется ориентация ВИП, а по математическим выражениям третьего, четвертого и пятого этапов «Точной выставки» алгоритма «Выставка» определяются соответствующие углы:
Figure imgf000013_0004
.
Figure imgf000013_0005
Figure imgf000013_0003
=-arctglzl /z2J (Pj = arctg
Figure imgf000014_0001
I [2(W2jWA j + Wl W3y ) + [2[W jW4j - j*>2j )
Данные по скорости и углам позволяют определить скорости в системе координат Север-Восток-Высота. Интегрируя эти скорости по времени, получаются значения соответствующих координат:
vEj =vuv>i * Zj+vnuV2. sin ZJ,
VNJ = -v„„Vl. sin*, + vnuV2j coszj ,
1 Ί
Figure imgf000014_0002
- cos ^y sin Xj
COS q)j COS Zj
sin cj
Ω4, = 2у+ у,где
WL и Wj - начальное и текущее значения нормированного кватерниона, определяющего угловое положение НСК;
J=LA...N;
и Qj - начальное и текущее значения вектора угловой скорости Земли в
НСК;
Σ . - вектор транспортирующей угловой скорости;
Zj - элементы третьего столбца матрицы направляющих косинусов по компонентам нормированного кватерниона;
Vmv - вектор линейных скоростей в НСК; vE , vNj , Vy , eJ } rij , и j - соответственно скорости и координаты
Восток-Север-Высота;
ек - эксцентриситет Земли;
г - экваториальный радиус Земли.
Алгоритм «Расчет радиусов» 13 (фиг. 2): так как изгиб трубопровода может происходить как в вертикальном, так и в горизонтальном направлении, можно вьщелить горизонтальную и вертикальную составляющие кривизны, где кривизной изгиба трубопровода называется величина, обратная радиусу изгиба трубопровода. Для вычисления кривизн траектории используются значения линейной скорости прибора и углов азимута и тангажа. Значение горизонтальной кривизны определяется как отношение производной по времени от азимута к линейной скорости, вертикальной - как отношение производной по времени от тангажа к линейной скорости.
Figure imgf000015_0001
Figure imgf000015_0002
соответственно горизонтальный и вертикальный малые радиусы;
НВл и rVBt ' соответственно горизонтальный и вертикальный большие радиусы.
Figure imgf000016_0001
AS - фиксированное значение базовой дистанции для определения больших радиусов, в этом случае приращение угла определяется по значениям углов на краях AS, а значение рассчитанного радиуса ставится в соответствии точке посередине AS,
Δ^- и AQJ - соответствующие средние оценки углов на базовой дистанции;
Figure imgf000016_0002
Кнв[гнБш ГНвъа^ - нв„], RVB[rvBla rVBzab-»rVBx] и аналогично для векторов радиусов с индексом f, где
Ls, LB -соответствующие параметры окна фильтра.
Алгоритм «Коррекция траектории» 14 (фиг.2): при использовании координат корректирующих точек, полученных на основе данных спутниковой системы GPS или (и) GLONASS, происходит коррекция полученной траектории и ее перенос на координатную сетку Север-Восток-Высота:
АЕ = \{АЕ К,М},
ej=ej + Aej,
AN = 1{AN K,M],
rij = rij + An j ,
AU = \{AU K,M),
и J = и j + Auj , где
Aej, Atij, AUj - корректирующие добавки по соответствующим осям; j=La...N, AE=[AeLa AeLa+1... AeN], AN=[A.nLa AnLa+1... ΔηΜ], AV=[AuLa AuLa+t ... A N];
M} - операция линейной интерполяции К значений на M=N-La;
AE'= [Δβ,' Ae2' ... Аек <], AN'=[A«,' Аи2' ... Αηκ'],
Δί/' = [Aw, ' Аи2' ... Аик'] - соответствующие разницы по координатам корректирующих точек и точек полученной траектории для соответствующих дистанций;
Figure imgf000017_0001
- координаты корректирующей точки для ^cp^ дистанции (дистанция, для которой измерены координаты корректирующей точки);
[efc nfc щ] - координаты корректируемой траектории для sk дистанции, при этом две указанные дистанции должны быть максимально близки;
К - количество корректирующих точек.
Способ определения величины и направления перемещения трубопровода по результатам двух и более инспекционных пропусков ДКОПТ предполагает использование данных по измеренным БИНС навигационным параметрам одного и того же участка трубопровода. Предложенный способ определения относительного перемещения трубопровода по результатам двух и более инспекционных пропусков ДКОПТ наглядно позволяет определить перемещение трубопровода в пространстве. Расположенная по таблицам информация о пропусках ДКОПТ по одному и тому же участку трубопровода совмещается по дистанциям. Пропуск ДКОПТ по инспектируемому участку трубопровода с более ранней датой считается базовым (фиг. 3). Последующий пропуск ДКОПТ по инспектируемому участку трубопровода (фиг. 4) сравниваются с базовым. Прогиб участка трубопровода может происходить как в вертикальном, так и в горизонтальном направлениях. Для этого выделяются вертикальная составляющая кривизны изгиба трубопровода: базовая 21 (фиг. 5) и последующая 23 (фиг. 6), и горизонтальная составляющая кривизны изгиба трубопровода: базовая 22 (фиг. 5) и последующая 24 (фиг. 6). Кривизной изгиба участка трубопровода называется величина, обратная радиусу изгиба участка трубопровода:
E=l/R, где
R - радиус изгиба участка трубопровода.
Между повторным и базовым пропусками ДКОПТ по участку инспектируемого трубопровода вычисляется разностный вектор:
Ё=(ЕН' - EH, EV' - EV), где EH 22 (фиг. 5) и EH 24 (фиг. 6) соответственно горизонтальная составляющая кривизны, выявленная по результатам базового пропуска ДКОПТ по участку инспектируемого трубопровода, и горизонтальная составляющая кривизны, выявленная по результатам последующего пропуска ДКОПТ по участку инспектируемого трубопровода;
EV 21 (фиг. 5) и EV 23 (фиг. 6) соответственно вертикальная составляющая кривизны, выявленная по результатам базового пропуска ДКОПТ по участку инспектируемого трубопровода, и вертикальная составляющая кривизны, выявленная по результатам последующего пропуска ДКОПТ по участку инспектируемого трубопровода.
Модуль разностной кривизны:
\Ё - Ш)2 + (ЕУ - F1 2
Превышение модуля разностной кривизны выше некоторого заданного порогового значения является критерием наличия перемещения на инспектируемом участке трубопровода.
Из навигационных данных определяются координаты точки начала участка перемещения 17 (фиг. 3) и конца участка перемещения 18 (фиг. 3) инспектируемого участка трубопровода при базовом пропуске ДКОПТ. При повторном пропуске ДКОПТ по инспектируемому участку трубопровода точки начала и конца участка перемещения соответственно 19 (фиг. 4) и
20 (фиг. 4).
Для того, чтобы вычислить величину вертикального перемещения внутри участка, необходимо приравнять координаты на границах участка перемещения инспектируемого участка трубопровода, а затем вычислить разницу координат внутри участка. Для этого рассчитывается:
Figure imgf000018_0001
где t и 2 соответственно вертикальные координаты начала и конца участка перемещения инспектируемого участка трубопровода при базовом пропуске ДКОПТ; и Z2 соответственно вертикальные координаты начала и конца участка перемещения инспектируемого участка трубопровода при последующем пропуске ДКОПТ.
После этого вычисляется фактическая разность координат по вертикали в точке р инспектируемого участка трубопровода между базовым и повторном пропусках ДКОПТ на инспектируемом участке трубопровода:
Figure imgf000019_0001
L - длина участка перемещения;
L] - дистанция от начала участка перемещения до выбранной точки р.
Для того, чтобы вычислить величину горизонтального перемещения внутри участка перемещения инспектируемого участка трубопровода, необходимо развернуть траектории базового и повторного пропуска ДКОПТ в горизонтальной плоскости таким образом, чтобы координаты начала и конца участка перемещения по оси Y стали нулевыми. Для этого рассчитываются координаты развернутых траекторий:
Figure imgf000019_0002
Prol , х
Prol - ;)2 +0 - ')2 Ψ' , yP' rol = ^Ρ' - χ[ )2 +{ p - y[)2 sin ^' ,
Ух
ψ = arctg - arctg У 2 ~ Ух
x xl x2 X]
Figure imgf000019_0003
Фактическая разность коор инат в горизонтальной плоскости: ХГ
Figure imgf000019_0004
- Урго, -
Использование раздельной процедуры расчета перемещений в вертикальной и горизонтальной плоскости позволяет снизить ошибки определения перемещений, когда траектории трубопровода при базовом и последующем пропусках в горизонтальной плоскости стремятся к ортогональному положению.
Абсолютная величина перемещения:
Figure imgf000020_0001
Максимальная величина перемещения участка трубопровода может быть определена как максимальное значение абсолютного перемещения.
Направление рассчитанного вектора перемещения 27 (фиг. 7) необходимо определить в плоскости 26 (фиг. 7), перпендикулярной оси трубопровода, за 0° принимается положение «на 12 часов». Внутри участка перемещения участка инспектируемого трубопровода вычисляются векторы направления движения. Точками начала и конца вектора предлагается установить координаты сварных швов отдельных секций трубопровода.
Единичный вектор нап авления оси трубопровода 25 (фиг. 7):
Figure imgf000020_0002
Ny%, ΝΖί - координаты начала отдельной секции трубопровода; NXz , Myz l NZz - координаты конца отдельной секции трубопровода;
ПХ1 fiy,, z - координаты единичного вектора 25 (фиг. 7).
Рассчитанный вектор перемещения 27 (фиг. 7):
A = (Axp , Ayp , Azp ) .
Проекция рассчитанного вектора перемещения 28 (фиг. 7), перпендикулярная направлению движения рабочей среды, трубопровода рассчитывается по формуле:
Ар = А - (А · п) - п .
Затем рассчитывается коэффициент:
К - А- П = Ахр · пх +Аур · ηγ +Δζρ · ηζ. .
Тогда с учетом рассчитанного коэффициента проекция рассчитанного вектора перемещения 28 (фиг. 7) -
Ар = Δ - К · Я - (Ах - К · п х , Ау - К · ηγ , Αζ р - К · η ζ ) ,
\Ар\ = ~jApx 2 + ApY 2 + Ap , где Δρχ, Лру, Δρζ - координаты проекции рассчитанного вектора 28 (фиг. 7) на плоскость 26 (фиг. 7).
Проекция вертикали 30 (фиг. 7) на плоскость 28 (фиг. 7), перпендикулярную направлению движения рабочей среды в трубопроводе:
2ρ= -ηζ·ηχζ·ηγ,\-ηζ·ηζ),
Figure imgf000021_0001
ζ РХ ' ΡΥ > Δ ΡΖ - координаты проекции вертикали 30 (фиг. 7) на плоскость 28 (фиг. 7).
Косинус искомого угла направления перемещения 31 (фиг. 7) находится по формуле:
Figure imgf000021_0002
Направление перемещения:
S = % - Пх + Пуу - nY + nn · п 12 '
Пу = {ηνχ
Figure imgf000021_0003
(Ару · ζη - Αρζ · zPy , Αρχ · ζΡΖ - Αρζ · ζρχ , Αργ · ζΡΧ - Αρχ · ζΡΥ )
Если S>0, a =arccos С . Если S<0, а - 360 - arccosC .

Claims

Формула изобретения
1. Устройство диагностического комплекса для определения положения трубопровода состоящее из аппаратной части: акселерометров, гироскопов и одометра, отличающееся тем, что состоит из аппаратной части и программной части, при этом аппаратная часть установлена на внутритрубный инспекционный прибор и состоит из набора датчиков: акселерометр, связанный с осью X бесплатформенной инерционной навигационной системой, акселерометр, связанный с осью Y бесплатформенной инерционной навигационной системой, акселерометр, связанный с осью Z бесплатформенной инерционной навигационной системой, гироскоп, связанный с осью X бесплатформенной инерционной навигационной системой, гироскоп, связанный с осью Y бесплатформенной инерционной навигационной системой, гироскоп, связанный с осью Z бесплатформенной инерционной навигационной системой, одометр.
2. Устройство по п.1 , отличающееся тем, что программная часть состоит из алгоритмов определения навигационных параметров в следующей последовательности: расчет линейной скорости, при этом определение движения внутритрубного инспекционного прибора происходит по данным, полученным от аппаратной части устройства диагностического комплекса для определения положения трубопровода; выставка, то есть определение положения в пространстве внутритрубного инспекционного прибора производится определением начальных углов азимута и тангажа в прямоугольной системе координат Север-Восток-Высота с учетом вращения Земли и заданным геодезическим координатам внутритрубного инспекционного прибора в течение первых нескольких минут, когда внутритрубный инспекционный прибор неподвижен; расчет навигационных параметров происходит с использованием значения дискрета времени между отчетами акселерометров и гироскопов и показаниями одометра с определением соответствующих углов: крена, тангажа, азимута; расчет радиусов изгиба трубопровода как в вертикальном, так и горизонтальном направлении; коррекция траектории происходит при использовании координат корректирующих точек, полученных привязкой конкретных точек трубопровода к геодезическим координатам в системе Север-Восток-Высота с использованием глобальной спутниковой навигационной системы GPS или/и GLONASS.
3. Устройство по п.п. 1 и 2, отличающееся тем, что после пропуска инспекционного прибора по трубопроводу, с установленной на нем аппаратной частью диагностического комплекса для определения положения трубопровода данные, полученные с аппаратной части диагностического комплекса для определения положения трубопровода, представляют собой массив, содержащий показания трех акселерометров и трех гироскопов в связанной с ортогональными осями бесплатформенной инерционной навигационной системы и соответственно с внутритрубного инспекционного прибора прямоугольной системой координат и показания одометра.
4. Устройство по п.1 , отличающееся тем, что ось у бесплатформенной инерционной навигационной системы соответствует продольной оси внутритрубного инспекционного прибора, а массив данных переписывается с аппаратной части диагностического комплекса для определения положения трубопровода, и расчет всех навигационных параметров производится на ЭВМ с помощью специальной программы.
5. Способ определения положения трубопровода, отличающийся тем, что используются данные измерения величин и направлений линейных перемещений участка трубопровода по результатам двух и более инспекционных пропусков внутритрубного инспекционного прибора с установленным на нем диагностическим комплексом для определения положения трубопровода по одному и тому же участку трубопровода с одним и тем же направлением потока рабочей среды.
6. Способ по п. 5, отличающийся тем, что данные, полученные с диагностического комплекса для определения положения трубопровода, располагаются по таблицам и совмещаются по дистанциям.
7. Способ по п.5, отличающийся тем, что данные пропуска внутритрубного инспекционного прибора с установленной на нем аппаратной частью диагностического комплекса для определения положения трубопровода с более ранней датой считаются базовыми, а данные последующих пропусков сравниваются с базовыми.
8. Способ по п.7, отличающийся тем, что критерием наличия перемещения трубопровода на инспектируемом участке является превышение модуля разностной кривизны заданного порогового значения.
PCT/RU2014/000859 2013-12-17 2014-11-12 Устройство и способ для определения положения трубопровода Ceased WO2015094015A1 (ru)

Priority Applications (4)

Application Number Priority Date Filing Date Title
MX2015016216A MX2015016216A (es) 2013-12-17 2014-11-12 Dispositivo y metodo para determinar la posicion de una linea de tuberia.
BR112015029372A BR112015029372A2 (pt) 2013-12-17 2014-11-12 dispositivo e método para a determinação de posição da tubulação
EA201591997A EA033119B1 (ru) 2013-12-17 2014-11-12 Способ для определения положения трубопровода
EP14870724.3A EP2985509A4 (en) 2013-12-17 2014-11-12 Device and method for determining position of pipeline

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
RU2013155927 2013-12-17
RU2013155927/06A RU2558724C2 (ru) 2013-12-17 2013-12-17 Устройство диагностического комплекса для определения положения трубопровода и способ определения относительного перемещения трубопровода по результатам двух и более инспекционных пропусков диагностического комплекса для определения положения трубопровода

Publications (1)

Publication Number Publication Date
WO2015094015A1 true WO2015094015A1 (ru) 2015-06-25

Family

ID=53403198

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/RU2014/000859 Ceased WO2015094015A1 (ru) 2013-12-17 2014-11-12 Устройство и способ для определения положения трубопровода

Country Status (6)

Country Link
EP (1) EP2985509A4 (ru)
BR (1) BR112015029372A2 (ru)
EA (1) EA033119B1 (ru)
MX (1) MX2015016216A (ru)
RU (1) RU2558724C2 (ru)
WO (1) WO2015094015A1 (ru)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111624552A (zh) * 2020-05-25 2020-09-04 中国地质大学(武汉) 一种基于声波渡越时间测量的地下管道定位系统及方法
CN114166214A (zh) * 2021-11-12 2022-03-11 海南大学 基于sta/lta算法的管道连接器检测方法和系统
CN114485642A (zh) * 2022-01-25 2022-05-13 中船航海科技有限责任公司 一种基于惯性测量的油气管道故障定位方法
CN114608570A (zh) * 2022-02-25 2022-06-10 电子科技大学 一种多模式自切换的管线仪自适应精密定位方法
CN117287644A (zh) * 2023-08-28 2023-12-26 同济大学 数据采集设备以及排水管道缺陷状态诊断方法

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015147686A1 (ru) 2014-03-28 2015-10-01 Открытое акционерное общество "Акционерная компания по транспорту нефти "Транснефть" (ОАО "АК "Транснефть") Способ контроля положения трубопроводов надземной прокладки в условиях вечной мерзлоты
RU2621219C1 (ru) * 2016-05-04 2017-06-01 Федеральное государственное бюджетное образовательное учреждение высшего образования "Саратовский государственный технический университет имени Гагарина Ю.А." (СГТУ имени Гагарина Ю.А.) Способ идентификации смещений осевой линии трубопровода
RU2622619C1 (ru) * 2016-09-28 2017-06-16 Общество с ограниченной ответственностью "Институт Рациональных Технологий" Система и способ экспресс-диагностирования сетей газопотребления
RU2666387C1 (ru) * 2017-04-04 2018-09-07 Федеральное государственное бюджетное образовательное учреждение высшего образования "Уфимский государственный нефтяной технический университет" Способ выявления геодинамических зон, пересекающих магистральные трубопроводы
RU2655614C1 (ru) * 2017-05-29 2018-05-29 Публичное акционерное общество "Транснефть" (ПАО "Транснефть") Способ измерения радиусов изгиба трубопровода на основе данных диагностического комплекса для определения положения трубопровода
CN109387197B (zh) * 2017-08-03 2022-04-12 北京自动化控制设备研究所 一种螺旋前进设备导航误差补偿方法
CN109612667B (zh) * 2019-01-09 2020-10-09 上海卫星工程研究所 利用陀螺数据辨识卫星挠性附件在轨端部位移方法与系统
CN110132256B (zh) * 2019-04-29 2021-06-15 北京航空航天大学 一种基于管道内检测器的定位系统及方法
CN110793488B (zh) * 2019-11-07 2021-08-20 南昌工程学院 一种水电机组转子测圆装置及其测圆调整计算方法
CN111220113B (zh) * 2020-01-13 2021-10-19 哈尔滨工程大学 一种管道拐弯角检测方法
CN111536969B (zh) * 2020-04-16 2022-12-13 哈尔滨工程大学 一种基于初始姿态角自对准的小径管道机器人定位方法
RU2747385C1 (ru) * 2020-10-30 2021-05-04 Александр Евгеньевич Зорин Способ определения пространственного положения трубопровода
RU2750417C1 (ru) * 2020-12-03 2021-06-28 федеральное государственное бюджетное образовательное учреждение высшего образования «Санкт-Петербургский горный университет» Способ определения изгибных напряжений в стенке подземного трубопровода
CN112628524B (zh) * 2020-12-31 2023-03-24 重庆安全技术职业学院 一种基于拐弯角的小径管道机器人高精度定位方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1404815A1 (ru) 1986-07-23 1988-06-23 Башкирский государственный университет им.40-летия Октября Устройство дл определени пространственного положени магистральных трубопроводов
RU1809297C (ru) 1990-12-10 1993-04-15 Инженерный Центр По Диагностике Трубопроводов Научно-Исследовательского И Проектно-Конструкторского Института Сооружения Объектов Нефтяной И Газовой Промышленности В Прикаспийском Нефтегазовом Комплексе Устройство дл измерени пространственного положени трубопроводов
US5331578A (en) * 1990-09-14 1994-07-19 Deutsche Forschungsanstalt Fur Luft- Und Raumfahrt E.V. Procedure for measuring angles and trajectories by means of gyros and inertial systems
RU2106569C1 (ru) 1996-01-03 1998-03-10 Центральный научно-исследовательский институт "Гидроприбор" Устройство для контроля профиля внутренней поверхности, пространственного положения и напряженного состояния трубопровода
US6243657B1 (en) * 1997-12-23 2001-06-05 Pii North America, Inc. Method and apparatus for determining location of characteristics of a pipeline
US6553322B1 (en) * 1999-09-29 2003-04-22 Honeywell International Inc. Apparatus and method for accurate pipeline surveying
RU2334162C1 (ru) 2007-01-26 2008-09-20 Общество С Ограниченной Ответственностью "Газпромэнергодиагностика" Аппаратура для измерения линейных деформаций магистрального трубопровода

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6170344B1 (en) * 1999-09-29 2001-01-09 Honeywell Inc. Pipeline distortion monitoring system

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1404815A1 (ru) 1986-07-23 1988-06-23 Башкирский государственный университет им.40-летия Октября Устройство дл определени пространственного положени магистральных трубопроводов
US5331578A (en) * 1990-09-14 1994-07-19 Deutsche Forschungsanstalt Fur Luft- Und Raumfahrt E.V. Procedure for measuring angles and trajectories by means of gyros and inertial systems
RU1809297C (ru) 1990-12-10 1993-04-15 Инженерный Центр По Диагностике Трубопроводов Научно-Исследовательского И Проектно-Конструкторского Института Сооружения Объектов Нефтяной И Газовой Промышленности В Прикаспийском Нефтегазовом Комплексе Устройство дл измерени пространственного положени трубопроводов
RU2106569C1 (ru) 1996-01-03 1998-03-10 Центральный научно-исследовательский институт "Гидроприбор" Устройство для контроля профиля внутренней поверхности, пространственного положения и напряженного состояния трубопровода
US6243657B1 (en) * 1997-12-23 2001-06-05 Pii North America, Inc. Method and apparatus for determining location of characteristics of a pipeline
US6553322B1 (en) * 1999-09-29 2003-04-22 Honeywell International Inc. Apparatus and method for accurate pipeline surveying
RU2334162C1 (ru) 2007-01-26 2008-09-20 Общество С Ограниченной Ответственностью "Газпромэнергодиагностика" Аппаратура для измерения линейных деформаций магистрального трубопровода

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP2985509A4

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111624552A (zh) * 2020-05-25 2020-09-04 中国地质大学(武汉) 一种基于声波渡越时间测量的地下管道定位系统及方法
CN111624552B (zh) * 2020-05-25 2022-08-30 中国地质大学(武汉) 一种基于声波渡越时间测量的地下管道定位系统及方法
CN114166214A (zh) * 2021-11-12 2022-03-11 海南大学 基于sta/lta算法的管道连接器检测方法和系统
CN114166214B (zh) * 2021-11-12 2023-06-06 海南大学 基于sta/lta算法的管道连接器检测方法和系统
CN114485642A (zh) * 2022-01-25 2022-05-13 中船航海科技有限责任公司 一种基于惯性测量的油气管道故障定位方法
CN114485642B (zh) * 2022-01-25 2024-04-16 中船航海科技有限责任公司 一种基于惯性测量的油气管道故障定位方法
CN114608570A (zh) * 2022-02-25 2022-06-10 电子科技大学 一种多模式自切换的管线仪自适应精密定位方法
CN117287644A (zh) * 2023-08-28 2023-12-26 同济大学 数据采集设备以及排水管道缺陷状态诊断方法

Also Published As

Publication number Publication date
EA033119B1 (ru) 2019-08-30
EP2985509A4 (en) 2017-03-22
EP2985509A1 (en) 2016-02-17
BR112015029372A2 (pt) 2017-07-25
EA201591997A1 (ru) 2016-09-30
RU2558724C2 (ru) 2015-08-10
MX2015016216A (es) 2016-03-11
RU2013155927A (ru) 2015-06-27

Similar Documents

Publication Publication Date Title
WO2015094015A1 (ru) Устройство и способ для определения положения трубопровода
CN109974697B (zh) 一种基于惯性系统的高精度测绘方法
US11428823B2 (en) Methods, apparatuses, and computer programs for estimating the heading of an axis of a rigid body
JP6060642B2 (ja) 自己位置推定装置
Ladetto et al. Digital magnetic compass and gyroscope integration for pedestrian navigation
Li et al. An improved optimal method for initial alignment
CN105371844B (zh) 一种基于惯性/天文互助的惯性导航系统初始化方法
CN103900565B (zh) 一种基于差分gps的惯导系统姿态获取方法
CN109870173A (zh) 一种基于校验点的海底管道惯性导航系统的轨迹修正方法
CN106507913B (zh) 用于管道测绘的组合定位方法
CN105091907B (zh) Sins/dvl组合中dvl方位安装误差估计方法
CN106289246A (zh) 一种基于位置和姿态测量系统的柔性杆臂测量方法
CN109489661B (zh) 一种卫星初始入轨时陀螺组合常值漂移估计方法
CN103900607A (zh) 一种基于惯性系的旋转式捷联惯导系统转位方法
CN103604428A (zh) 基于高精度水平基准的星敏感器定位方法
Santana et al. Estimation of trajectories of pipeline PIGs using inertial measurements and non linear sensor fusion
US20140249750A1 (en) Navigational and location determination system
Yang et al. Application of adaptive cubature Kalman filter to in-pipe survey system for 3D small-diameter pipeline mapping
Avrutov Autonomous determination of initial latitude with an inertial measuring unit
CN103955005B (zh) 一种火箭橇轨道重力实时测量方法
CN103954288A (zh) 一种卫星姿态确定系统精度响应关系确定方法
US20170138738A1 (en) Remote location determination system
RU2539131C1 (ru) Бесплатформенная интегрированная навигационная система средней точности для мобильного наземного объекта
Zhang et al. Integrated navigation method for SINS and odometer
Jaekel et al. An inertial navigation system with acoustic obstacle detection for pedestrian applications

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 14870724

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2014870724

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 201591997

Country of ref document: EA

WWE Wipo information: entry into national phase

Ref document number: MX/A/2015/016216

Country of ref document: MX

REG Reference to national code

Ref country code: BR

Ref legal event code: B01A

Ref document number: 112015029372

Country of ref document: BR

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 112015029372

Country of ref document: BR

Kind code of ref document: A2

Effective date: 20151124