Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Подземка

.pdf
Скачиваний:
29
Добавлен:
22.05.2015
Размер:
2.22 Mб
Скачать
o x0 y0 z0 .

Для придания невозмущаемости

погрешностей системы линейными

 

 

 

ускорениями составляющие вектора

h

в (15) должны вычисляться по

значениям линейной скорости (VN ,VE ) и значениям широты места,

выработанным в БИНС на основе информации, полученной от акселерометров СП.

 

 

E

VN

 

 

 

 

 

 

VE

 

 

 

 

 

 

 

, N cos cos

 

,

 

R

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

VE

 

 

 

 

 

 

 

h sin sin

 

tg ,

 

 

(16)

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

7,292116 10 51/ c

- угловая скорость суточного вращения Земли, а

 

 

a(1 e2 )

 

 

 

 

 

a

 

 

 

 

 

 

R

 

 

h,

R

 

 

 

 

 

h, -

радиусы кривизны

 

(1 e2 sin 2 )3 / 2

 

 

(1 e2 sin 2

)1/ 2

эллипсоида Красовского:

e2

0,0066934 ,

a 6378245

м.

 

 

 

Задача преобразования

 

кажущихся ускорений

в

бесплатформенном

гироинклинометре в режиме БИНС формулируется как задача нахождения

проекций nN , nE , nh вектора

 

на оси географического трехгранника по

n

информации

об измеренных

акселерометрами

составляющих nх0, ny 0, nz 0

 

 

 

 

 

 

o x0 y0 z0 и по

вектора n в осях связанного с

корпусом СП трехгранника

вычисленной

матрице

C 0

,

определяющей

взаимную

ориентацию

 

 

h

 

 

 

 

трехгранников oNEh и

В блоке выработки параметров поступательного движения осуществляется вычисление составляющих линейной скорости и координат

места путем интегрирования составляющих nN , nE , nh .

Исходя из выражения (15) могут быть получены кинематические уравнения углового движения бесплатформенного инклинометра в углах

Эйлера:

N sin A E cos A x0 sin y0 cos ,

Asin h sin cos N cos A E sin A x0 cos y0 sin ,

(17)

sin N cos A E sin A z 0 sin cos y0 sin x0 cos .

71

ориентацию
72
o x0 y0 z0

Из (15) можно получить соотношения, связывающие угловые скорости

вращения трехгранника o x0 y0 z0 с угловыми скоростями вращения географического трехгранника oNEh и скоростями изменения углов азимута

A , зенитного угла и угла поворота :

 

 

 

x0

sin h A sin cos N sin A sin cos A cos cos

 

 

 

y0

E cos A sin sin A cos cos ,

cos h A sin sin N sin A cos cos A cos sin (18)

 

 

 

E cos A cos sin A cos sin ,

z0 h A cos sin N cos A E sin A .

Для автокомпенсации инструментальных погрешностей чувствительных элементов в гироинклинометре могут быть предусмотрены модуляционные развороты измерительного блока (ИБ), включающего блок акселерометров и блок гироскопов, вокруг продольной оси СП в диапазоне 180 . Тогда

положение связанной

с ИБ системы координат

o xb yb zb

относительно

 

 

 

 

 

 

 

 

 

cos

sin

0

 

 

трехгранника

o x

 

y

 

z

 

характеризуется матрицей C b sin

cos

0

, где

 

 

 

0

 

0

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- текущее значение угла разворота вокруг продольной оси СП, которое вычисляется с помощью датчика угла (ДУ), расположенного на этой оси.

Матрица направляющих косинусов Coh , определяющая ориентацию географического трехгранника oNEh (h) относительно трехгранника

(o), вычисляется как:

 

c11

c12

c13

 

 

 

 

 

h

 

 

 

 

 

 

 

h

b

 

 

C0

c21

c22

c23

 

Cb

Co

,

(19)

 

c

31

c

32

c

33

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Cob - матрица направляющих косинусов, определяющая системы координат o xb yb zb (b) относительно o x0 y0 z0 (o).

Вычисление азимутального, зенитного углов A , и угла может быть

осуществлено по вычисленным элементам матрицы направляющих

косинусов Ch .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A arctg

 

c

23

 

 

1 (c )2

arctg

 

d

32

 

 

 

 

 

, arctg

 

 

33

 

,

 

 

.

(20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c13

 

 

c33

 

 

 

 

 

 

d 31

 

 

Модель погрешностей определения азимута в трехосной схеме

 

 

 

 

 

 

 

гироинклинометра в режиме БИНС

 

 

Матрицы Ch , Cb , C h определяются в БИНС с погрешностями, т. е. могут

 

 

 

 

 

b o

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

быть определены лишь некоторые оценки C h , C b ,

C h этих матриц, причем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

o

 

0

 

 

C h C h

,

 

 

 

 

 

Ch Ch ,

 

 

 

 

 

(21)

C h

Cb Cb Cb ,

Ch

 

 

 

 

 

o

o

o

 

 

 

o

o

 

o

b

 

 

 

b

b

 

 

 

 

 

 

где

C h , C b

, C h

-

матрицы,

элементы

 

которых

определяются

 

 

b

o

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

варьированием соответствующих матриц Ch

, Cb

, C h .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

o

 

 

 

0

 

 

Подставляя соотношения (21) в алгоритм (19) и, при этом не учитывая

величин второго порядка малости относительно рассматриваемых переменных, получим следующее уравнение:

C h C h C b C h

C b

(22)

 

0

 

 

b o

 

 

b

 

o

 

Умножая левую и правую части равенства (21) на матрицу C 0

, получим

 

 

 

 

 

 

 

 

 

h

 

следующее выражение для кососимметрической матрицы погрешностей:

Ch

Ch

Ch Cb

Cb ,

(23)

 

0

 

b

b

0

 

 

h

 

и соответствующему ей вектору погрешностей:

 

 

 

 

 

C h

 

 

 

 

 

(24)

 

0h

 

bh

 

0b

,

 

 

 

b

 

 

 

 

которые характеризуют суммарные погрешности гироинклинометра в построении географической системы координат относительно осей корпуса СП.

Составляющие веторов в выражении (24) имеют вид:

0h

sin A cos Asin

cos A sin Asin ,A cos

 

 

 

,

 

0

 

(25)

 

 

 

0

,

 

 

 

 

 

 

 

 

 

bh

 

 

0b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

73

Если считать, что ДУ измеряет угол без ошибок (т.е. =0) и

представить векторное уравнение (24) в скалярной форме, то получим уже известные кинематические соотношения (12) и (13) для погрешностей в выработке параметров ориентации.

Воспользовавшись выводом модели погрешнстей БИНС на ДУС,

приведенным в , можно получить уже известное кинематическое уравнение,

в котором «дрейф восточного гироскопа» E в бесплатформенном гироинклинометре определяется выражением (11).

74

Лекция № 8

Двухосные гироинклинометры. Режим точечного компасирования в

продольной схеме гироинклинометра: алгоритмы, погрешности

Применение бесплатформенной инклинометрической системы,

построенной на одном двухосном ДУС c расположением вектора кинетического момента по оси скважины (продольная схема) - наиболее естественное решение в условиях малого диаметра скважин.

Вычисление азимута скважины при использовании продольной схемы с

одним двухосным ДУС, измеряющим только две составляющие x0 , y 0

угловой скорости вращения Земли, традиционно производится так же, как и в трехосной схеме, т.е. с помощью функцию arctg :

A arctg

y cos

(1)

 

 

x sin sin

где: y x0 sin y0 cos , x x0

cos y0 sin

в основе которой лежат следующие выражения, вытекающие из первых двух соотношений в (1) лекции №9:

sin A

y

(2)

 

 

 

cos

 

 

cos A

x

sin sin

(3)

cos cos

 

 

Для погрешности определения азимута нетрудно получить:

A

x sin A y cos cos A

 

 

 

 

 

(4)

 

cos cos

 

 

 

(cos tg sin cos A sin 2 A sin tg ) sin A(tg cos A tg ),

 

где: y

x 0 sin y 0

cos , x x 0 cos y 0 sin

 

Первое слагаемое уравнения (4) определяет «гироскопическую» составляющую ошибки определения азимута, выражение для которой подтверждает известное положение о неработоспособности продольной схемы гироинклинометра при приближении зенитного угла к 90 .

Второе и третье слагаемые уравнения (4) определяют составляющую

75

ошибки азимута, обусловленную погрешностями акселерометров. «Акселерометрическая» составляющая A в (4) содержит члены, которые пропорциональны tg и при значении 90 , также как и «гироскопическая» составляющая (4), стремятся к бесконечности. Таким образом, даже при идеальном гироскопе погрешность вычисления азимута по алгоритму (1) в

продольной схеме будет расти с увеличением угла из-за погрешностей акселерометров.

Сегодняшнее применение продольной схемы в режиме точечного компасирования ограничено зенитными углами - 60 -75 .

76

Лекция № 9

Продольная схема гироинклинометра: измерение угловых параметров наклонно-направленных и горизонтальных скважин,

особенности построения траекторий скважин.

Если в качестве источника информации об азимуте A рассматривать только выражение (2) (см. лекцию №8), то можно получить следующий

алгоритм его вычисления:

A arcsin

y

 

 

 

(5)

 

 

 

 

 

 

cos

 

 

 

Варьируя (5), получим:

 

 

 

y

1

 

 

A

 

 

(cos tg sin

 

)

(6)

cos cos A

cos A

Как видим, в отличие от (4) «гироскопическая» составляющая ошибки вычисления азимута (5) вообще не зависит от значений зенитного угла, а ее

«акселерометрическая» составляющая свободна от зависимости от tg .

Таким образом, ошибки определения азимута (1) и (5) обладают разными свойствами. Пользуясь этими соотношениями и при этом, учитывая только составляющие, обусловленные погрешностями ДУС, можно определить условие, при котором для вычисления азимута предпочтительней

пользоваться алгоритмом (5).

Значения ошибки азимута A, полученные из выражения (4) и (6)

обозначим соответственно At

и As .

 

Проверяя условия выполнения неравенства:

 

 

 

 

A

 

 

 

 

A

 

 

At

 

1,

 

 

 

 

 

 

<

 

 

или

 

(7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

s

 

 

 

 

As

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

получим:

 

cos2

A

 

sin Acos A

 

1, и, далее, используя

свойство

 

 

 

 

 

 

 

 

 

 

 

x

cos

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

транзитивности и известное соотношение модуля суммы и суммы модулей,

77

получим: cos2 A

 

y

 

sin Acos A

 

 

1. Теперь

условие, при котором

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

cos

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

cos A

 

 

 

 

 

выполняется A A , будет иметь вид:

 

 

 

 

 

 

 

sin A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

s

 

 

 

 

 

 

x

 

 

 

cos

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И, наконец, рассматривая (7) применительно к максимально допустимой

ошибке

 

и

принимая

 

условие

равноточности погрешностей ДУС

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

x

 

max

 

y

 

 

max ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получим:

 

 

tgA

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

(8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, при выполнении условия (8) следует пользоваться для вычисления A алгоритмом (1), в противном случае ─ алгоритмом (5).

Области предпочтительного использования обоих алгоритмов вычисления угла A наглядно иллюстрируются на картинной плоскости

(рис.1). Очевидно, что суммарно эти области гораздо шире, чем области использования традиционного алгоритма, что делает перспективным применение комбинированного алгоритма, предполагающего переключение с одного алгоритма на другой.

, град

 

 

 

 

 

 

 

 

 

 

 

180

 

 

 

 

 

 

 

 

 

 

 

 

160

 

 

 

 

Область

 

 

 

 

 

 

 

 

 

 

использования

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

140

 

 

 

 

"полного"

 

 

 

 

 

 

120

 

 

 

 

алгоритма

 

 

 

 

 

 

 

Область

 

 

 

 

 

 

 

Область

 

 

 

 

 

 

 

 

 

 

 

 

 

100

 

использования

 

 

 

 

 

 

использования

 

80

 

"усеченного"

 

 

 

 

 

 

"усеченного"

 

 

алгоритма

 

 

 

 

 

 

алгоритма

 

 

 

 

 

 

 

 

 

 

60

 

 

 

 

 

Область

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

использования

 

 

 

 

 

40

 

 

 

 

 

"полного "

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

 

 

 

 

 

алгоритма

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0

15

30

45

60

75

90

105

120

135

150

165

180

 

 

 

A, град

 

 

 

 

Рис.1

 

 

 

 

 

 

 

В то же время при значении зенитного угла =

90 азимут можно

 

 

вычислить только в соответствии с алгоритмом (5).

 

 

 

 

 

 

 

Анализ показал, что при пересечении траекторией

скважины

зоны,

граничные значения азимутов которой равны 90 (270 )

 

 

, 90 (270 )

 

 

 

 

2

y

 

2

y

 

 

 

 

 

 

78

(где

y

y

 

), приборное

значение

Ac

перестает меняться.

cos

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Гироинклинометр

оказывается в

«мертвой

зоне»

протяженностью 2 2 y

(рис.2, рис.3). Такую же величину имеет диапазон возможных значений,

принимаемых Ac (рис.2, рис.3). При попадании в этот «квадрат» и нахождении в нем сверх той протяженности проходки, которая обеспечивает выполнение критерия гладкости:

 

dA

(l l

 

) 2

 

 

 

 

 

 

 

0

2

y

,

 

 

(9)

 

dl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где:

dA

─ допустимое искривление скважины (

 

);

(l l0 ) ─ протяженность

 

dl

м

 

 

 

 

 

 

 

 

 

проходки на рассматриваемом интервале,

теряется возможность измерения угловых параметров скважины. При выполнении условия (9) выбор между значениями Ac и (180 - Ac )

(«определение четверти») неоднозначен - выход азимута из «мертвой зоны» становится неконтролируемым, а дальнейшая траектория скважины – ненаблюдаемой и приходится констатировать невозможность решения задач проводки или контроля скважины гироинклинометром в режиме точечного компасирования.

A (град)

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

1

при значении

 

 

 

 

при значении 0

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

A,град

0

15

30

45

60

75

90

105

120

135

150

165

180 град

-1

при значении 0

 

 

 

 

при значении 0

 

 

 

 

 

 

 

 

 

 

 

 

-2

 

 

 

 

 

 

 

 

 

 

 

 

-3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 2

 

 

 

 

 

 

79

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

A

 

 

 

 

 

 

 

 

 

(град)

-5

-4

-3

-2

-1

1

2

3

4

5

 

 

 

 

 

-1

 

 

 

 

 

 

 

 

 

 

-2

 

 

 

 

 

 

 

 

 

 

-3

 

 

 

 

 

Рис.3

Таким образом, подводя итог всему вышесказанному, можно сделать

следующие выводы:

1. Для измерения угловых параметров скважин произвольной ориентации с помощью продольной схемы в режиме точечного компасирования следует пользоваться комбинированным, т. е. меняющимся в зависимости от углов ориентации ( A , ) алгоритмом ( arcsin или arctg ), при

этом условием «перехода» с

одного алгоритма на другой является

 

cos tgA

 

1. Следует отметить,

что область применения алгоритма arcsin

 

 

 

 

 

 

 

гораздо шире, чем arctg , а для горизонтальной скважины с зенитным углом

= 90 этот алгоритм является единственным, позволяющим вычислить ее азимут.

2.Анализ погрешностей компасирования участков скважин, лежащих в плоскости горизонта и приближающихся к линии «восток - запад» показывает, что и в этом случае сопровождение бурения скважины одногироскопным скважинным прибором возможно, а его погрешности ограничены. Однако, при достаточной протяженности такого участка может произойти потеря наблюдаемости его траектории и, в силу этого, дальнейшая невозможность проводки любого другого горизонтального участка.

80