Научная статья на тему 'ИДЕНТИФИКАЦИЯ ПАРАМЕТРОВ ИНКЛИНОМЕТРИЧЕСКИХ ДАТЧИКОВ'

ИДЕНТИФИКАЦИЯ ПАРАМЕТРОВ ИНКЛИНОМЕТРИЧЕСКИХ ДАТЧИКОВ Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
11
5
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ИНКЛИНОМЕТР / ВЕКТОР / ДАТЧИК КООРДИНАТ ВЕКТОРА ГЕОФИЗИЧЕСКОГО ПОЛЯ / ТЕНЗОР / МАТРИЦА

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Рогатых Н.П.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

IDENTIFICATION OF THE PARAMETERS OF INCLINOMETER SENSORS

The article deals with simple ways to identify the parameters of inclinometer sensors using rotary units with inaccurate spatial orientation and arbitrary installation of sensors in inclinometers.

Текст научной работы на тему «ИДЕНТИФИКАЦИЯ ПАРАМЕТРОВ ИНКЛИНОМЕТРИЧЕСКИХ ДАТЧИКОВ»

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

Blockchain в свою очередь способствует децентрализации, а при помощи сети компьютеров происходит перевод денег, все переводы чисты и видны сообществу системы. Решается проблема отмывки денег, решается проблема воровства в банковской системе. Так как каждый компьютер представляет собой записную книжку, в которую будет записан перевод наших денег. Из выше сказанного может показать, что система blockchain плохо защищена, но это не так. Blockchain использует криптографические ключи, причём, например в BitCoin используется два ключа. Открытый(puЫic) и закрытый (Private). Что такое Public key и для чего он нужен? Публичный ключ поможет расшифровать хеш транзакции, который имеет у всех пользователей системы. Можно сказать, что это и есть номер нашей карты при переводе для системы то есть программно. Что такое Private key и для чего он нужен? Приватный ключ имеется только у хозяина транзакции, собственно при помощи его и отправляются деньги, если грубо, он формирует начальную хеш строку, при помощи публичного ключа, можно подтвердить (расшифровать) транзакцию. Секретный ключ не может быть раскрыт. [1]

Так где же BitCoin применят Blockchain?

Так где же BitCoin применят Blockchain? Наши транзакции с хеш функциями и есть цепь блоков, которая идёт друг за другом и если один из блоков будет изменён, мошенник попытается изменить в нашей книге сумму перевода, то предыдущей и последующий блок отвергнет изменённый блок, так как запись нового хеша будет принципиально другой.

Использованные источники: 1. Блокчейн, Википедия [Электронный ресурс] — режим доступа: https://ru.wikipedia.org/wiki/Блокчейн (дата обращения: 17.03.2017);

УДК 53.089.6

Рогатых Н.П., к. техн. н. доцент, ведущий инженер-конструктор Катав-Ивановский приборостроительный завод

Россия, г. Катав-Ивановск ИДЕНТИФИКАЦИЯ ПАРАМЕТРОВ ИНКЛИНОМЕТРИЧЕСКИХ ДАТЧИКОВ

Аннотация:

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

Ключевые слова:

Инклинометр, вектор, датчик координат вектора геофизического поля, тензор, матрица.

Rogatykh N. P., candidate of technical sciences, associate Professor

Leading design engineer Katav-Ivanovsk instrument engineering plant

Russia, Katav-Ivanovsk IDENTIFICATION OF THE PARAMETERS OF INCLINOMETER SENSORS

Annotation:

The article deals with simple ways to identify the parameters of inclinometer sensors using rotary units with inaccurate spatial orientation and arbitrary installation of sensors in inclinometers.

Keywords:

Inclinometer, vector, coordinate sensor of geophysical field vector, tensor, matrix.

Современные инклинометры, а также магнитометры, компасы и т.п. строятся на основе трёхкомпонентных датчиков координат векторов геомагнитного и гравитационного полей. Такие датчики включают в себя три элементарных однокомпонентных датчика, сигнал каждого из которых U линейно связан с проекцией измеряемого вектора физического поля а на его ось чувствительности с и имеет вид

U = кс ■ а + U0 ^

где к и Ц0 - коэффициент преобразования и смещение нулевого уровня сигнала датчика. Соответственно вектор сигналов трёхкомпонентного датчика, объединяющий сигналы трёх элементарных датчиков, представляется в виде уравнения

О=кса + и0 ф

где К - тензор коэффициентов преобразования трёхкомпонентного датчика с диагональной матрицей, С - тензор ориентации осей чувствительности трёхкомпонентного датчика с матрицей направляющих косинусов этих осей в системе координат, связанной с корпусом

инклинометра, и °° - вектор смещений нулевых уровней сигналов датчика, а измеряемый вектор при работе инклинометра определяется по алгоритму

а = с-К~\и - ио). (3)

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

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

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

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

и = (к + Ак )(с + Ас) • а + (и0 + Аио) ^

где к, с, и0 - известные паспортные или заданные значения параметров датчиков, а Дк, Ас, ДЦ0 - погрешности соответствующих параметров датчиков или поправки, приводящие заданные значения параметров к реальным. Здесь

а = (А + АЛ)(а0 + Аа0) ^^

где А - известный тензор, описывающий задаваемые с помощью поворотной установки пространственные повороты инклинометра, ДА -тензор, учитывающий инструментальные погрешности поворотной

установки, а° - исходный измеряемый вектор с известными координатами, Аа

0 - вектор, учитывающий погрешности задания координат исходного вектора.

Погрешности параметров датчиков определяются из (4) в зависимости от задаваемых значений вектора а и соответствующих значений сигнала датчика и. При этом для упрощения математических выражений и вычислений часто делается ряд допущений, искажающих реальные характеристики датчиков. В частности, в [1, 2] считается, что коэффициент преобразования датчика имеет точное паспортное значение, т.е. Дк = 0, единичный направляющий вектор оси чувствительности датчика непременно совпадает с одной из осей связанного с цилиндрическим корпусом инклинометра ортогонального базиса, а неточность установки осей чувствительности рассматривается в виде малых приращений направляющих векторов.

Такой подход имеет широкое распространение, как в отечественной, так и в зарубежной практике, и лежит в основе многих инструкций по калибровке инклинометров. Он даёт вполне удовлетворительные результаты, если погрешности датчиков малы и находятся в пределах технологических допусков, но при существенных отклонениях параметров датчиков от заданных значений может привести к погрешностям, которые не компенсируются при калибровке.

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

с=(1,0,°), (6)

то его малое приращение будет иметь вид

Ас = (0,8Х, 32), (7)

где 51 и 52 - величины, характеризующие погрешности установки оси чувствительности. При этом длина суммы векторов (6,7) составит

|с + М| = ,112 + £2 + ^2 >1, (8)

а коэффициент преобразования изменится на величину

Ак = к(^1 + 3? + -1). (9)

Это приращение коэффициента преобразования не будет учтено при калибровке датчика.

Отметим, что рассмотренные в [1, 2] приёмы калибровки нельзя применить к тем датчикам, у которых оси чувствительности принципиально не совпадают с указанным выше ортогональным базисом, как, например, в [5], а и к датчикам с произвольной ориентацией осей чувствительности, т.к. в этом случае требуются другие представления математических моделей.

Задание координат исходного вектора а0 при использовании известных способов калибровки инклинометров осуществляется за счёт предварительной ориентации поворотных установок относительно сторон света и местной вертикали. При этом используются специальные приборы, например, оптические квадранты и магнитные буссоли, собственные погрешности которых и погрешности, связанные с креплением этих приборов в поворотных установках, обуславливают погрешность

Аа

определения координат исходного вектора 0 , которая также снижает точность калибровки датчиков.

Таким образом, существенными недостатками известных способов калибровки датчиков являются применение несовершенных математических моделей, которые не соответствуют реальным характеристикам датчиков, и необходимость точной пространственной ориентации поворотных установок, связанная с использованием дополнительных приборов и затратами времени. В результате этих недостатков снижаются эффективность калибровки датчиков и точность инклинометров.

Исключить указанные недостатки позволяет рассматриваемый ниже способ идентификации параметров датчиков, частный случай которого был изложен в [5]. В основе способа лежит более общая и в то же время простая по сравнению с известными математическая модель элементарного датчика, а также ряд математических приёмов, позволяющих получить результаты, не

зависящие от координат исходного вектора геофизического поля а°.

Эта модель представляется в виде

и = п ■а + и°, (10)

где п = кс - обобщённый параметрический вектор, характеризующий одновременно коэффициент преобразования датчика, расположение осей чувствительности датчика и длину измеряемого вектора, т.е. напряжённость измеряемого физического поля. Введение такого вектора оправдано тем, что целью идентификации параметров инклинометрических датчиков, в конечном счёте, являются не значения коэффициентов преобразования и координат осей чувствительности датчиков, а значения абстрактных коэффициентов, устанавливающих зависимость между сигналами датчиков и направлением измеряемого вектора.

В соответствии с (10) вектор сигналов трёхкомпонентного датчика будет иметь вид

и = Ш + Оо, (11)

где N - параметрический тензор, а единичный направляющий вектор измеряемого поля будет определяться как

а = N-\и- ио). (12)

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

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

п = (п, п, П )

а = Аср Ав Ааа0 а0 _ (а01' а02, а03 )

(13)

(14)

(15)

А Ав Аа

где п1, п2, п3 -координаты параметрического вектора, тензоры, описывающие пространственные повороты датчика с инклинометром, ф, в, а - углы, задаваемые с помощью стандартной поворотной установки, а01, а02, а03 - координаты исходного вектора.

Если ввести неподвижный ортогональный базис Я0(0Х0У070), направив его ось 070 по гравитационной вертикали места вниз, а ось 0У0 -по касательной к магнитному меридиану на Север, и при этом совместить эти оси соответственно с осью вращения рамки а поворотной установки и нулевой отметкой шкалы отсчёта угла а (Рис.1), а также ввести подвижный ортогональный базис Я1(0Х1У171), совместив его ось 071 с осью вращения рамки

Рис.1

Ф и одновременно с продольной осью корпуса инклинометра, а 0X1 - с нулевой отметкой шкалы ф, то указанные углы будут представлять собой известные в инклинометрии визирный угол ф (Tool Face), зенитный угол в и магнитный азимут а. При этом матрицы соответствующих тензоров будут иметь вид [3]

cos^ sin^

\ =

sm^ 0

cos^ 0

A =

cos в 0

- sin в

0 1 0

sin в 0

cos в

А =

cos« - sin« 0 sin« - cos« 0

0

0

1

(16)

В дальнейшем мы будем использовать поворотную установку с произвольной пространственной ориентацией относительно поверхности Земли и в связи с этим рассматривать углы ф, в, а только в качестве углов относительных поворотов рамок установки.

Заметим, что оси вращения рамок в базисе R0 в любом случае задаются направляющими векторами = (cos a sin в, sin a sin в, cos в), = (sin a, cos а, 0) и = (0,0,1) (Рис.2).

Рис.2

При вращении датчика вокруг некоторой оси его сигнал в зависимости от угла поворота у может быть представлен в виде

U Acos^ + B sin^ + C (17)

где A, B, C - постоянные коэффициенты, зависящие от направления оси вращения датчика. Разбиение диапазона изменения угла у (00...3600) на p равных частей позволяет вычислить эти коэффициенты с помощью известных формул, применяемых для анализа гармоник рядов Фурье, -

2 p 2 p 1 p a = -E U cosB = -X U sin^. C = -X U

p 1=i , p 1=i , p 1=i , (18) =—(i -1)

где p и Ui - значения угла и сигнала датчика,

соответствующие текущему номеру их отсчёта i,

Используя (10, 13, 14, 15, 16), представим сигнал датчика аналогичным образом в виде функции угла ф (у=ф)

U = A (e, a) cos ( + B (e, a) sin( + C (e, a)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(19)

где

A = ai [n (- sinecosa) + n2 (sina)]

+ a02 [n (- sinesina) + n2 (- cosa)]+ a03 [n sine] B( = ai [n (sina) + n2 (cosecosa)]

+ ao2 [n (- cosa) + n (cosesina)]+ a03 [n (- sine)] C = n [am (sinecosa) + a02 (sinesina) + a03 (cose)] + U0

(20)

Теперь эти коэффициенты определим в виде функций угла а, как Av = A^ cosa + B^ sina + С

Bp = A cosa + B sina + С ^^

С = Ac cosa + Bc sina + Сс

Аа = -а{)Хпх - я02я2 ВА = а01п2 - а02п бш^ С = а^щ

где а "01' А = атп2 соб6

А = атпъ

а02П1

Вв = атпх + а^2п2 собО Сд = -а03п2

В =

^ 03 2 1

сс = аоъпъ соб^ + и0

Полученные коэффициенты (22) запишем для двух значений в, равных 00 и 900, в результате чего для в =00 будем иметь

Аа (0) = -ао2П2, Ва (0) = ^, Са (0) = 0,

Ав (0) = а0Лп2 - а0Ппх Вв (0) = а0П + а0)2п2 Св (0) = 0

Вс (0) = 0

Сс (0)

а03п3 + Ц0

Ас (0) = 0 , а для в =900 -

Аа(90) = - а0)2п2, ВА(90) = а0)Хп2 - а0)2п1, Са(90) = а0зп,

ав (90) = -а02п1 вв (90) = а01п1 св (90) = -а03п2

(23)

вс (90) = а02п3

Сс (90) = Ц

(24)

(25)

ас (90) = а01пз, Отсюда следует, что п = д/А2В (90) + ВВ(90) + СА (90),

п2 ^л/[аВ(90)-ВА(90)]2Т[АА(90)+ВВ(90)]ГТСВ(90)

пз =у1 ас (90) + ВС (90) + [СС (0) - СС (90)]2 , Ц = Сс (90).

Очевидно, что вычисленные по формулам (25) абсолютные значения координат параметрического вектора и значение смещения нулевого сигнала

датчика не зависят от абсолютных значений координат исходного вектора а0 . Однако знаки координат параметрического вектора всё-таки определяются

знаками координат исходного вектора, Например, если а(п>0, "°2>0, "°3>0 (Рис.1), то

п1 > 0 при ВВ(90) > 0, п1 < 0 при ВВ(90) < 0,

п2 > 0 при СВ(90) < 0, п2 < 0 при СВ(90) > 0, (26)

п3 > 0 при АС(90) > 0, п3 < 0 при АС(90) < 0.

Поэтому для полной идентификации параметров датчика следует

знать, в каком из квадрантов базиса Я0 изначально находится вектор а0. Это условие является единственным для исходной ориентации поворотной установки.

Таким образом, для идентификации параметров датчика измеряются значения сигнала датчика при вращении инклинометра вокруг оси 070 (в=00) и нескольких осей в плоскости Х00У0 (в=900), направления которых в этой плоскости задаются рядом значений угла а. При этом количество значений а выбирается достаточным для определения необходимых коэффициентов (23, 24).

Интересно отметить, что, в случае использования прецизионных

ап

ап

датчиков и достаточно точной идентификации их параметров инклинометры в дальнейшем могут служить в качестве эталонных приборов для задания пространственной ориентации поворотных установок.

Для идентификации параметров современных датчиков инклинометров, обладающих высокой линейностью характеристик, нет необходимости в большом количестве измерений сигнала датчика. В частности, для определения коэффициентов А, В, С в уравнении (17) достаточно измерить всего три значения сигнала, соответствующие трём равноотстоящим друг от друга значениям угла например, 00, 1200 и 2400, и вычислить коэффициенты с помощью простых формул

А=и (о0)-С В[и (12°0)" и(240°)],

С =1 [и (00) + и (1200) + и (2400)]

3 . (27)

Применяя эти формулы для идентификации параметров датчика, достаточно будет измерить значения сигнала датчика в двенадцати указанных в таблице 1 положениях инклинометра и с помощью полученных Таблица 1

в а

Ф

00 900

00 00 1200 2400

00 U1 U4 U7 U10

1200 U2 U5 U8 U11

2400 U3 U6 U9 U12

значений U1... U12 вычислить коэффициенты (23,24) по формулам Cc (0) = ±(Ц + U 2 + Uз)

Aa (90) = -(U7 + Uw) + 1(U7 + U8 + U9 + Uw + Un + Uj2)

(90) ^

(U7 - UJ - 1(U7 + U8 + U9 - Uw - Un - U12) 1

Ca (90) = (U4 + U7 + U10)--(U4 + U5 + U6 + U7 + U8 + U9 + Uw + Un + U12)

As(90) = -1(U5 - U6) - ^[(U5 + U8 + Un) - (U6 + U9 + U12)]

BB(90) = 1 [(U8 - U9) - (Un - U12)]

^ (90) ^ + U8 + Un) " {U6 + 09 + U12^

ac (90) = 1(U4 + u5 + U6) - 1(U4 + u5 + U6 + U7 + U8 + U9 + U10 + U11 + U12)

ßc(90) = j-j= (U7 + U8 + U - U10 - Un - U12)

Cc (90) = 1(U4 + U5 + U6 + U7 + U8 + U9 + U10 + Un + U12)

Использование этих значений в (25) позволяет полностью идентифицировать все параметры датчика. Очевидно, что для одновременной идентификации параметров всех датчиков инклинометра достаточно будет также задать ему только двенадцать пространственных положений.

Рассмотренный способ идентификации параметров

инклинометрических датчиков свободен от недостатков, связанных с применением некорректных математических моделей датчиков и необходимостью точного определения координат исходного вектора, т.е. точной пространственной ориентации поворотных установок. Погрешности данного способа в основном зависят от уровня квантования сигналов датчиков при аналого-цифровом преобразовании, степени нелинейности статических характеристик датчиков и точности используемых поворотных установок. В частности, использование 12-разрядных АЦП и средней по точности поворотной установки для идентификации трёхкомпонентных гравиметрических датчиков на основе акселерометров позволяет получить погрешность измерения зенитного угла в среднем Ад =(0,1...0,15)0, а визирного угла - Аф =(0,15.0,2)0 при д =1,00. Идентификация магнитометрических датчиков для тех же условий обеспечивает погрешность измерения азимута Аа =(1,5.2,0)0 при использовании феррозондов и Аа =(0,2.0,5)0 при использовании магниторезистивных преобразователей.

Таким образом, рассмотренный способ идентификации параметров датчиков допускает произвольное размещение датчиков в корпусе инклинометра и использование поворотных установок без строгой пространственной ориентации. Это позволяет более гибко решать задачи построения конструкций инклинометров и осуществлять оперативные настройки инклинометров в полевых условиях, что имеет немаловажное значение для развития отечественного геофизического приборостроения.

Использованные источники:

1. Ковшов Г.Н. Инклинометры (Основы теории и проектирования) / Ковшов Г.Н., Алимбеков Р.И., Жибер А.В. - Уфа: Гилем, 1998. - 380 с.;

2. Миловзоров Г.В. Анализ инструментальных погрешностей

инклинометрических устройств / Г.В. Миловзоров. - Уфа: Гилем, 1997. -184с.;

3. Рогатых Н.П. Векторные математические модели гравимагнитных инклинометров / Н.П. Рогатых // Изв. Вузов / Нефть и Газ.-1990.-№11.- с. 8185;

4. Рогатых Н.П. Определение параметров трёхкомпонентных инклинометрических датчиков // Методы и средства измерения в системах контроля и управления: Труды международной научно-технической конференции под. ред. М.А.Щербакова. - Пенза: Информационно-издательский центр Пенз. гос. ун-та, 2002 - с.167;

5. Рогатых Н.П. Гравиметричекий датчик // Патент РФ №37523, 2004, Бюл. №12.

i Надоели баннеры? Вы всегда можете отключить рекламу.