СПОСОБ ОПРЕДЕЛЕНИЯ КАЛИБРОВОЧНОЙ МАТРИЦЫ ТРЕХКОМПОНЕНТНОГО МАГНИТОМЕТРА (ВАРИАНТЫ) Российский патент 2005 года по МПК G01R33/02 

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

Прецизионные измерения индукции магнитного поля основаны на выявлении и учете систематических погрешностей измерения. Одной из доминирующих составляющих основной погрешности является статическая погрешность. Эта погрешность снижется или устраняется на этапе калибровки магнитометра при его настройке и/или на последующих этапах поверки (коррекции) магнитометра. Отечественная терминология в этой области не в полной мере соответствует международной. По международному словарю основных терминов в метрологии ИСО, 1993, в результате калибровки выполняется любая из двух родственных операций:

- приписывание показаниям прибора размера единицы измеряемой величины;

- определение коррекции показаний в единицах измеряемой величины.

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

Отечественный документ "Рекомендации по межгосударственной стандартизации РМГ 29-99" принял упомянутый словарь за основу, но в отечественной измерительной практике термин "калибровка" по-прежнему понимается в узком пассивном смысле, не включающем в объем понятия термин "коррекция". Достаточно сказать, что "Толковый словарь по метрологии. Основные термины" (Корнеева Т.В., 1990 г., с.133) определяет термин "калибровка средства измерения" как "Поверка средства измерения посредством совокупных измерений". Т.е. русский термин предполагает только первую из двух упомянутых операций международного термина "калибровка". Ниже термин "калибровка" используется в широком смысле, охватывающем термин "коррекция", но иногда дополнительно используется термин "коррекция".

Известно, что статическая погрешность определяется искажающей матрицей трехкомпонентного магнитометра (см. Внучков Г.А. и др. Угловые погрешности измерения магнитной индукции. - В сб. Методы и средства исследования структуры геомагнитного поля. М., ИЗМИРАН, 1989 г., с.161-168). Там же предложен способ определения боковых элементов искажающей матрицы трехкомпонентного магнитометра с помощью трехкомпонентной меры магнитной индукции в условиях, когда магнитометр и мера неортогональны. Однако известный способ не дает возможности определить все элементы искажающей матрицы, а потому не дает возможности найти калибровочную (корректирующую) матрицу по искажающей.

Известен способ калибровки результатов измерения трехкомпонентным магнитометром (см. патент США №6130534, кл. G 01 R 35/00, G 01 R 33/02, 2000 г.), основанный на определении калибровочной (корректирующей) матрицы. Этот способ по совокупности существенных признаков наиболее близок заявляемому и принят за прототип.

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

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

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

Калибровочная матрица может зависеть от ряда условий (Подробнее см. в разделе Сведения, подтверждающие возможность осуществления изобретения). Это разнообразие условий привело к составлению заявки на способ, состоящий из четырех вариантов, обеспечивающих получение одного и того же технического результата.

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

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

Указанный технический результат при осуществлении изобретения - способа по 1-му варианту, достигается тем, что в заявляемом способе определения корректирующей матрицы трехкомпонентного магнитометра, включающем размещение датчика магнитометра в центре образцовой трехкомпонентной меры магнитной индукции, ориентацию чувствительных осей магнитометра вдоль осей меры, компенсацию внешнего магнитного поля, воспроизведение мерой точного значения магнитной индукции вдоль каждой ее оси, регистрацию всех компонент результатов измерения на выходе магнитометра, соответствующих магнитному полю, воспроизводимому мерой, построение нормированной приборной матрицы из нормированных вектор-столбцов результатов измерения и использование ее для нахождения калибровочной матрицы, в отличие от известного способа, предварительно, используя независимый модульный магнитометр, измеряют неортогональности осей образцовой меры, по которым строят матрицу осей меры b_м в собственной ортогональной системе координат меры, затем в центре меры устанавливают калибруемый трехкомпонентный магнитометр и воспроизводят магнитное поле В_мi поочередно при положительной и отрицательной полярности вдоль каждой оси меры, при этом регистрируют соответствующие результаты измерения трех компонент на выходе магнитометра и определяют нормированную приборную матрицу П_v по формуле

П_v=[(V_uм1+-V_uм1-)/2В_м1 (V_uм2+-V_uм2-)/2В_м2 (V_uм3+-V_uм3-)/2В_м3],

где V_uмi+, V_uмi- - вектор-столбцы результатов измерения напряжения на выходе магнитометра, соответствующие магнитному полю, воспроизводимому мерой при положительной и отрицательной полярности;

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

и разлагают ее на произведение

S_м=S·I_м,

где S=(S₁ S₂ S₃) - диагональная матрица чувствительности с размерностью В/Т;

- матрица единичных чувствительных осей магнитометра в виде вектор-строк I_iм в собственной системе координат меры;

с помощью чего определяют три неортогональности осей е_ij по формуле

e_ij=(l_i1l_j1+l_i2l_j2+l_i3l_j3)_м;

по которым строят матрицу осей магнитометра I₀ в его собственной ортогональной системе координат, в которой определяют матрицу чувствительности S₀=S·I₀ и находят искомую калибровочную матрицу .

Указанный технический результат при осуществлении изобретения - способа по 2-му варианту, достигается тем, что в заявляемом способе определения корректирующей матрицы трехкомпонентного магнитометра, включающем размещение датчика магнитометра в центре образцовой трехкомпонентной меры магнитной индукции, ориентацию чувствительных осей магнитометра вдоль осей меры, компенсацию внешнего магнитного поля, воспроизведение мерой точного значения магнитной индукции вдоль каждой ее оси, регистрацию всех компонент результатов измерения на выходе магнитометра, соответствующих магнитному полю, воспроизводимому мерой, построение нормированной приборной матрицы из нормированных вектор-столбцов результатов измерения и использование ее для нахождения калибровочной матрицы, в отличие от известного способа, предварительно, используя независимый модульный магнитометр, измеряют неортогональности осей образцовой меры, по которым строят матрицу осей меры b_м в собственной ортогональной системе координат меры, затем в центре меры устанавливают калибруемый трехкомпонентный магнитометр и воспроизводят магнитное поле В_мi поочередно при положительной и отрицательной полярности вдоль каждой оси меры, при этом регистрируют соответствующие результаты измерения трех компонент на выходе магнитометра и определяют нормированную приборную матрицу П_v по формуле

П_v=[(V_uм1+-V_uм1-)/2В_м1 (V_uм2+-V_uм2-)/2В_м2 (V_uм3+-V_uм3-)/2В_м3],

где V_uмi+, V_uмi- - вектор-столбцы результатов измерения напряжения на выходе магнитометра, соответствующие магнитному полю, воспроизводимому мерой при положительной и отрицательной полярности;

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

и разлагают ее на произведение

K_м=h_м·K,

где h_м=(h₁ h₂ h₃)_м - матрица единичных чувствительных осей магнитометра в виде вектор-столбцов h_iм в собственной системе координат меры;

K=(K₁ K₂ K₃) - диагональная калибровочная матрица с размерностью Т/В;

с помощью чего определяют три неортогональности осей μ_ij по формуле

μ_ij=(h_1ih_1j+h_2ih_2j+h_3ih_3j),

по которым строят матрицу осей магнитометра h₀ в его собственной ортогональной системе координат, в которой находят искомую калибровочную матрицу K_h0=h₀·K.

Указанный технический результат при осуществлении изобретения - способа по 3-му варианту, достигается тем, что в заявляемом способе определения корректирующей матрицы трехкомпонентного магнитометра, включающем размещение датчика магнитометра в центре образцовой трехкомпонентной меры магнитной индукции, ориентацию чувствительных осей магнитометра вдоль осей меры, компенсацию внешнего магнитного поля, воспроизведение мерой точного значения магнитной индукции вдоль каждой ее оси, регистрацию всех компонент результатов измерения на выходе магнитометра, соответствующих магнитному полю, воспроизводимому мерой, построение нормированной приборной матрицы из нормированных вектор-столбцов результатов измерения и использование ее для нахождения калибровочной матрицы, в отличие от известного способа, используя независимый модульный магнитометр, измеряют неортогональности осей образцовой меры, по которым строят матрицу осей меры b_м в собственной ортогональной системе координат меры, затем в центре меры устанавливают испытуемый трехкомпонентный магнитометр и воспроизводят магнитное поле В_мi поочередно при положительной и отрицательной полярности вдоль каждой оси меры, при этом регистрируют соответствующие результаты измерения трех компонент на выходе магнитометра и определяют нормированную приборную матрицу П_в по формуле

П_в=[(В_uм1+-В_uм1-)/2В_м1 (В_uм2+-В_uм2-)/2В_м2 (В_uм3+-В_uм3-)/2В_м3],

где В_uмi+, В_uмi- - вектор-столбцы результатов измерения магнитной индукции на выходе магнитометра, соответствующие магнитному полю, воспроизводимому мерой при положительной и отрицательной полярности;

затем находят искажающую матрицу магнитометра в собственной системе координат меры и_м по формуле

и разлагают ее на произведение

и_м=a·I_м,

где a=(a₁ a₂ a₃) - диагональная матрица безразмерных коэффициентов, незначительно отличающихся от единицы;

- матрица единичных чувствительных осей магнитометра в виде вектор-строк I_iм в собственной системе координат меры;

с помощью чего определяют три неортогональности осей е_ij по формуле

e_ij=(l_i1l_j1+l_i2l_j2+l_i3l_j3)_м,

по которым строят матрицу осей магнитометра I₀ в его собственной ортогональной системе координат, в которой определяют искажающую и₀=a·I₀ и находят искомую калибровочную матрицу .

Указанный технический результат при осуществлении изобретения - способа по 4-му варианту, достигается тем, что в заявляемом способе определения корректирующей матрицы трехкомпонентного магнитометра, включающем размещение датчика магнитометра в центре образцовой трехкомпонентной меры магнитной индукции, ориентацию чувствительных осей магнитометра вдоль осей меры, компенсацию внешнего магнитного поля, воспроизведение мерой точного значения магнитной индукции вдоль каждой ее оси, регистрацию всех компонент результатов измерения на выходе магнитометра, соответствующих магнитному полю, воспроизводимому мерой, построение нормированной приборной матрицы из нормированных вектор-столбцов результатов измерения и использование ее для нахождения калибровочной матрицы, в отличие от известного способа, предварительно, используя независимый модульный магнитометр, измеряют неортогональности осей образцовой меры, по которым строят матрицу осей меры b_м в собственной ортогональной системе координат меры, затем в центре меры устанавливают испытуемый трехкомпонентный магнитометр и воспроизводят магнитное поле В_мi поочередно при положительной и отрицательной полярности вдоль каждой оси меры, при этом регистрируют соответствующие результаты измерения трех компонент на выходе магнитометра и определяют нормированную приборную матрицу П_в по формуле

П_в=[(В_uм1+-В_uм1-)/2В_м1 (В_uм2+-В_uм2-)/2В_м2 (В_uм3+-В_uм3-)/2В_м3],

где В_uмi+, В_uмi- - вектор-столбцы результатов измерения магнитной индукции на выходе магнитометра, соответствующие магнитному полю, воспроизводимому мерой при положительной и отрицательной полярности;

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

и разлагают ее на произведение

k_м=h_м·k,

где h_м=(h₁ h₂ h₃)_м - матрица единичных чувствительных осей магнитометра в виде вектор-столбцов h_iм в собственной системе координат меры;

k=(k₁ k₂ k₃)_м - диагональная матрица безразмерных коэффициентов k_i;

затем определяют три неортогональности осей магнитометра μ_ij по формуле

μ_ij=(h_1ih_1j+h_2ih_2j+h_3ih_3j),

по которым строят матрицу осей магнитометра h₀ в его собственной ортогональной системе координат, в которой находят искомую калибровочную матрицу k_h0=h₀·k.

На чертеже изображена блок-схема реализации предлагаемого способа определения калибровочной матрицы трехкомпонентного магнитометра.

Блок-схема включает образцовую трехкомпонентную меру 1 магнитной индукции (МИ), в рабочей области которой вначале помещен модульный магнитометр 2 на этапе определения неортогональностей меры, затем - датчик 3 калибруемого трехкомпонентного магнитометра. Датчик 3 подключен к измерительному каналу магнитометра 4, выход которого связан с входом компьютера 5 для расчета калибровочных матриц трехкомпонентного магнитометра 3, 4 и для реализации процедур калибровки. Выход модульного магнитометра 2 также связан с входом компьютера 5.

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

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

В_+x+y=B_x+B_y В_-x+y=-B_x+B_y

В_+x-y=B_x-B_y В_-x-y=-B_x-B_y

модули которых измеряют прецизионным модульным магнитометром 2. Неортогональность меры 1 между ее осями 1 и 2 определяют (см. аналог Угловые погрешности измерения...) по формуле

Аналогично определяют остальные неортогональности Е₁₃ по комбинации полей В_х и В_z и Е₂₃ по комбинации В_у и В_z.

Затем в рабочую зону помещают датчик 3 калибруемого трехкомпонентного магнитометра 4, оси которого ориентируют вдоль осей меры, после чего вдоль каждой оси меры 1 по очереди воспроизводят мерой точно известное значение магнитного поля В_мi или (В_м1 В_м2 В_м3) при положительной и отрицательной полярности, при этом измеряют на выходе магнитометра 4 или соответствующие вектор-столбцы напряжения V_uм+=(V_uм1 V_uм2 V_uм3)₊ и V_uм-=(V_uм1 V_uм2 V_uм3)_-, или соответствующие вектор-столбцы магнитной индукции В_uм+=(В_uм1 В_uм2 В_uм3)₊ и В_uм-=(В_uм1 В_uм2 В_uм3)_-.

Результаты измерения модульным магнитометром 2 (или заранее определенные неортогональности меры 1 Е₁₂, Е₁₃, Е₂₃), значения полей, воспроизводимые мерой 1, В_мi, соответствующие результаты измерения на выходе магнитометра 4 V_uм+, V_uм- или В_uм+, В_uм- заносятся в компьютер 5. Компьютер 5 осуществляет нормирование результатов измерения, построение нормированной приборной матрицы

П_v=[(V_uм1+-V_uм1-)/2В_м1 (V_uм2+-V_uм2-)/2В_м2 (V_uм3+-V_uм3-)/2В_м3],

или

П_в=[(В_uм1+-В_uм1-)/2В_м1 (В_uм2+-В_uм2-)/2В_м2 (В_uм3+-В_uм3-)/2В_м3].

и вычисление по заданным нормированной матрице и неортогональностям меры варианта калибровочной матрицы магнитометра 3, 4:

K_0l - в выбранной собственной системе координат магнитометра (0) на магниточувствительных осях магнитометра типа (l), если результат измерения на выходе магнитометра задан в вольтах (вариант 1);

K_h0 - на магниточувствительных осях магнитометра типа (h) в выбранной собственной системе координат магнитометра (0), если результат измерения на выходе магнитометра задан в вольтах (вариант 2);

k_0l - в выбранной собственной системе координат магнитометра (0) на магниточувствительных осях магнитометра типа (l), если результат измерения на выходе магнитометра задан в теслах (вариант 3);

k_h0 - на магниточувствительных осях магнитометра типа (h) в выбранной собственной системе координат магнитометра (0), если результат измерения на выходе магнитометра задан в теслах (вариант 4).

После определения калибровочной матрицы магнитометр 3, 4 вынимают из меры 1 и последующие результаты измерения магнитометром калибруют (и корректируют) по схеме

B_0l=K_0l·(V_u-V₀) или B_h0=K_h0·(V_u-V₀)

B_0l=k_0l·(B_u-B₀) или B_h0=k_h0·(B_u-B₀),

где В_0l, В_h0 - действительные значения измеряемой величины в собственной ортогональной системе координат магнитометра соответственно (0l) или (h0);

результат измерения в вольтах; результат измерения в теслах;

V₀, В₀ - уходы нуля магнитометра соответственно в вольтах и в теслах.

Во всех вариантах вышеприведенной схемы основным корректирующим объектом являются та или другая калибровочная матрица. Очевидно, что калибровочные матрицы K_0l, K_h0 имеют размерность (Т/В), а калибровочные матрицы k_0l, k_h0 - безразмерны. Варианты 1 и 2 соответствуют случаю, когда результаты измерения на выходе магнитометра заданы в вольтах и соответственно задана нормированная приборная матрица П_v. Варианты 3 и 4 соответствуют случаю, когда результаты измерения на выходе магнитометра заданы в теслах и соответственно задана нормированная приборная матрица П_в. Различие вариантов 1 и 2 (а также 3 и 4) вызвано тем, что существуют два класса трехкомпонентных (феррозондовых) магнитометров с разнесенными и совмещенными центрами компонент датчиков. Считается, что для одного класса характерны магниточувствительные оси одного вида (I_i), для другого класса - оси другого вида (h_i). Подробнее на тему об осях I₁, и h_i см. ниже.

Процедура определения калибровочной матрицы в собственной ортогональной системе координат магнитометра по измеренным нормированной приборной матрице П_v или П_в и неортогональностям меры Е_ij неочевидна и нуждается в более подробных пояснениях.

Рассмотрим модель искажений вначале для однокомпонентного магнитометра. Результат измерения однокомпонентным магнитометром В_ul′ на некоторое единичное базовое направление l′ можно представить как проекцию измеряемого поля В на некоторую единичную ось магнитометра I_k

где коэффициент L_k учитывает отклонение от единицы в силу несовершенства настройки, нестабильности и несовершенства калибровки однокомпонентного магнитометра, - собственный уход нуля магнитометра по оси I_k. Другим источником погрешности является не строгое совпадение оси магнитометра I_k с базовой осью l′.

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

В_u=и·В+В₀,

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

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

Очевидно, что элементы каждой оси I_i зависят от выбора системы координат (СК), в которой эта ось рассматривается, соответственно матрица осей I тоже зависит от выбора СК. В отличие от I, матрица L - диагональная и не зависит от выбора системы координат. Она зависит только от качества настройки, стабильности и калибровки магнитометра. Чем лучше настройка, стабильность и калибровка i-той компоненты, тем ближе безразмерный коэффициент L_i к единице. Однако даже строгое равенство L_i=1 не гарантирует высокой точности измерения, т.к. в общем случае оси I_i неортогональны между собой и не совпадают полностью с ортами базовой СК. Угловые отклонения осей магнитометра сказываются на точности. Чем ближе матрица и к единичной матрице, тем точнее магнитометр.

Кроме того, результат измерения в теслах - это синтезированная величина, полученная из результата измерения вектор-столбца напряжения V_u на выходе магнитометра с помощью некоторой диагональной матрицы T с размерностью Т/В

В_u=T·V_u

Если в соотношении Т·V_u=и·В+В₀, измеряемое поле внешнего происхождение В заменить на поле, воспроизводимое мерой В_мi=В_мib_i вдоль ее i-той оси b_i, т.е. в положительном направлении, то получим нормированный к модулю поля меры В_мi вектор-столбец Т·V_uмi+/В_мi=и·b_i+В₀/В_мi. Если поменять полярность поля меры вдоль ее i-той оси, то ее ось меняется на противоположную и Т·V_uмi-/В_мi=-и·b_i+В₀/В_мi (модуль поля меры В_мi всегда остается со знаком плюс). Построим нормированную приборную матрицу П следующего вида П_v=[(V_uм1+-V_uм1-)/2В_м1 (V_uм2+-V_uм2-)/2В_м2 (V_uм3+-V_uм3-)/2В_м3]

Из вышеприведенных соотношений следует, что

П_v=T^-1·и·(b₁ b₂ b₃)=S·b

является произведением матрицы чувствительности магнитометра S=T^-1·и=T^-1 ·L·I=S·I на матрицу из вектор-столбцов осей меры b=(b₁ b₂ b₃). При этом уходы нулей магнитометра и остаточное внешнее поле взаимокомпенсируются. Обе матрицы S и b должны рассматриваться в одной и той же СК, каждая из них зависит от выбора СК, но их произведение П_v инвариантно к выбору СК. Итак, матрицы сомножители S и b неизвестны, но их произведение предполагается известным. Известны также неортогональности осей меры Е_ij=b_i·b_j, но сами оси b_i, b_j пока неизвестны.

Попытаемся найти матрицу осей меры b в некоторой СК. Для этого выберем такую ортогональную СК меры (ОСКМ), чтобы можно было бы упростить точное выражение всех элементов матрицы осей меры b_ij через ее неортогональности E_ij. Выберем ОСКМ с ортами ijk, например, так, чтобы орт i ОСКМ совпал с осью меры b₁, а ось меры b₂ лежала в плоскости ортов i, j

В выбранной ОСКМ, которую мы будем именовать собственной ОСКМ и обозначать через индекс "м", вся матрица осей меры b_м, следующим образом выражается через ее неортогональности Е_ij, т.е. через известные величины

где

Соотношения выше приведены в качестве примера. Собственную ОСКМ можно было бы выбрать, по крайней мере, девятью равноправными способами и получить похожие выражения матрицы осей через их неортогональности (шесть способов, когда матрица осей имеет три нулевых боковых элемента и еще три способа, когда матрица осей имеет два нулевых элемента, а два ненулевых элемента, симметричных относительно главной диагонали, равны между собой). Т.е. любую тройку неортогональных осей можно выразить точно и просто в одной из 9-ти собственных ортогональных СК через неортогональности осей.

Итак, по найденным неортогональностям меры строят матрицу осей b_м в одной из собственных ОСКМ и находят матрицу чувствительности магнитометра в той же ОСКМ по формуле

Очевидно, что найденную матрицу S_м можно однозначно разложить на сомножители: диагональную матрицу S=(S₁ S₂ S₃) и матрицу единичных вектор-строк или осей магнитометра в ОСКМ I_м, S_м=S·I_м, т.е. по известным элементам можно однозначно определить все элементы матриц S и I_м,

I_м=S^-1·S_м

Матрицу I_м, можно было бы использовать для калибровки (коррекции) в случае, если магнитометр выполнял бы последующие измерения в собственной ОСКМ, т.е. не меняя своей ориентации относительно меры. Но нас интересует возможность калибровки (коррекции) магнитометра, вынутого из меры и произвольно ориентированного, т.е. возможность коррекции в собственной ортогональной СК магнитометра (ОСКМАГ), жестко связанной с датчиком магнитометра.

Для определения ОСКМАГ находят по известной матрице

неортогональности магнитометра e_ij=(I_i·I_j)_м=(l_i1l_j1+l_i2l_j2+l_i3l_j3)_м, так как эти величины инвариантны к выбору СК. Далее по аналогии с вышерассмотренным примером выражения матрицы осей меры b_м в одной из собственных ОСКМ выбирают одну из собственных ОСКМАГ и определяют в ней матрицу осей магнитометра I₀ (индекс "0" означает принадлежность к собственной ОСКМАГ) через найденные неортогональности осей магнитометра е_ij.

Например, выберем ортонормированный базис ОСКМАГ ijk на осях магнитометра, так чтобы

Откуда получают матрицу осей магнитометра в выбранной ОСКМАГ I₀ через найденные неортогональности e_ij:

где

По полученным диагональной матрице чувствительности S и матрице осей магнитометра I₀ находят искомую калибровочную матрицу K_0l в выбранной собственной системе координат магнитометра (0) на его чувствительных осях I_i

где S^-1 - диагональная калибровочная матрица с размерностью Т/В не зависит от выбора СК.

Во втором варианте сначала определяют калибровочную матрицу K_м в выбранной ортогональной СК меры

Затем раскладывают ее, в отличие от первого варианта, на произведение матриц вида

K_м=h_м·K

где h_м=(h₁ h₂ h₃)_м - матрица единичных вектор-столбцов h_i или матрица других чувствительных осей магнитометра в выбранной собственной ОСКМ;

K=(K₁ K₂ K₃) - диагональная калибровочная матрица коэффициентов с размерностью Т/В, не зависящая от выбора СК.

По найденной матрице h_м определяют неортогональности осей

μ_ij=(h_i·h_j)_м=(h_1ih_1j+h_2ih_2j+h_3ih_3j)_м

и по найденным неортогональностям по аналогии с вышерассмотренным вариантом выбирают ту или другую собственную ОСКМАГ, в которой строят матрицу осей h₀ и находят искомую калибровочную

K_h0=h₀·K

В третьем варианте сначала находят искажающую матрицу магнитометра и_м в выбранной ОСКМ и разлагают ее на произведение диагональной безразмерной матрицы L и матрицы осей I_м

Затем по найденной матрице I_м, определяют неортогональности е_ij между осями по аналогии с первым вариантом, после чего выбирают собственную ОСКМАГ, в которой по найденным неортогональностям е_ij строят матрицу осей I₀ и определяют искажающую и₀=L·I₀, по которой находят искомую калибровочную

В четвертом варианте сначала находят калибровочную в ОСКМ k_м и разлагают ее на произведение матрицы осей h_м и диагональной безразмерной матрицы k

Затем по найденной матрице h_м определяют неортогональности μ_ij между осями по аналогии со вторым вариантом, после чего выбирают собственную ОСКМАГ, в которой по найденным неортогональностям μ_ij строят матрицу осей h₀ и определяют искомую калибровочную k_h0=h₀·k.

Калибровочные матрицы K_0l и K_h0, а также k_0l и k_h0 могут не совпадать друг с другом в зависимости от выбора собственной ОСКМАГ на единичных осях I_i или h_i.

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

В некоторых случаях дополнительно определяют и матрицу связи σ или ортогональную матрицу преобразования координат между выбранными собственными ОСКМ и ОСКМАГ

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

где σ′ - транспонированная σ, поскольку матрица σ ортогональна, то σ·σ′=I.

Найденная матрица σ характеризует тщательность угловой ориентации испытуемого неортогонального трехкомпонентного магнитометра в неортогональной трехкомпонентной мере.

Для проверки возможностей заявляемого способа в сравнении со способом прототипом с помощью системы автоматизации расчетов MatLab вычислены точностные показатели коррекции в виде максимального элемента матриц

для способа прототипа,

где - корректирующая по прототипу, записанная в наших обозначениях;

для заявляемого способа

в условиях, характеризуемых максимальными элементами матриц: осей магнитометра , осей меры и углового рассогласования между мерой и магнитометром . Результаты расчета сведены в таблицу 1.

Таблица 1Хар-ка магнитометраХар-ка мерыХар-ка тщательности ориентации магнитометра в мереРасчетная остаточная погрешность коррекции по прототипуУлучшение по прототипуРасчетная остаточная погрешность коррекции для заявляемого способа0.0150.000290.0580.05840.26^*)<2.2*10^-16--0.00580.00612.5-«»---0.000580.000916.7-«»---0.0000580.000350-«»-0.00150.00290.0580.06100.025^*)-«»---0.00580.00870.17^*)-«»---0.000580.00350.43^*)-«»---0.0000580.00300.5^*) 0.00150.000290.0580.05840.026^*)-«»---0.00580.00610.25^*) --0.000580.00091.7-«»---0.0000580.00035 *) улучшения нет**) расчеты с моделированием всех операций заявляемого способа.

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

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

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

1. Способ определения калибровочной матрицы трехкомпонентного магнитометра, включающий размещение датчика магнитометра в центре образцовой трехкомпонентной меры магнитной индукции, ориентацию чувствительных осей магнитометра вдоль осей меры, компенсацию внешнего магнитного поля, воспроизведение мерой точного значения магнитной индукции вдоль каждой ее оси, регистрацию всех компонент результатов измерения на выходе магнитометра, соответствующих магнитному полю, воспроизводимому мерой, построение нормированной приборной матрицы из нормированных вектор-столбцов результатов измерения и использование ее для нахождения калибровочной матрицы, отличающийся тем, что предварительно, используя независимый модульный магнитометр, измеряют неортогональности осей образцовой меры, по которым строят матрицу осей меры b_м, в собственной ортогональной системе координат меры, затем в центре меры устанавливают калибруемый трехкомпонентный магнитометр и воспроизводят магнитное поле B_мi поочередно при положительной и отрицательной полярностях вдоль каждой оси меры, при этом регистрируют соответствующие результаты измерения трех компонент на выходе магнитометра и определяют нормированную приборную матрицу П_v по формуле:

П_v=[(V_uм1+-V_uм1-)/2В_м1(V_uм2+-V_uм2-)/2В_м2(V_uм3+-V_uм3-)/2В_м3],

где V_uмi+, V_uмi- - вектор-столбцы результатов измерения напряжения на выходе магнитометра, соответствующие магнитному полю, воспроизводимому мерой при положительной и отрицательной полярностей;

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

и разлагают ее на произведение

S_м=S·I_м,

где S=(S₁ S₂ S₃) - диагональная матрица чувствительности с размерностью В/Т;

- матрица единичных чувствительных осей магнитометра в виде вектор-строк I_iм в собственной системе координат меры;

с помощью чего определяют три неортогональности осей е_ij по формуле:

e_ij=(l_i1l_j1+l_i2l_j2+l_i3l_j3)_м;

по которым строят матрицу осей магнитометра I₀ в его собственной ортогональной системе координат, в которой определяют матрицу чувствительности S₀=S·I₀ и находят искомую калибровочную матрицу .

П_v=[(V_uм1+-V_uм1-)/2В_м1(V_uм2+-V_uм2-)/2В_м2(V_uм3+-V_uм3-)/2В_м3],

где V_uмi+, V_uмi- - вектор-столбцы результатов измерения напряжения на выходе магнитометра, соответствующие магнитному полю, воспроизводимому мерой при положительной и отрицательной полярностях;

затем находят калибровочную матрицу магнитометра К_м в собственной системе координат меры по формуле

и разлагают ее на произведение

К_м=h_м·К,

где h_м=(h₁ h₂ h₃)_м - матрица единичных чувствительных осей магнитометра в виде вектор-столбцов h_iм в собственной системе координат меры;

K=(K₁ К₂ К₃) - диагональная калибровочная матрица с размерностью Т/В;

с помощью чего определяют три неортогональности осей μ_ij по формуле:

μ_ij=(h_1ih_1j+h_2ih_2j+h_3ih_3j),

по которым строят матрицу осей магнитометра h₀ в его собственной ортогональной системе координат, в которой находят искомую калибровочную матрицу К_h0=h₀·К.

П_в=[(В_uм1+-В_uм1-)/2В_м1 (В_uм2+-В_uм2-)/2В_м2 (В_uм3+-В_uм3-)/2В_м3],

где B_uмi+, B_uмi- - вектор-столбцы результатов измерения магнитной индукции на выходе магнитометра, соответствующие магнитному полю, воспроизводимому мерой при положительной и отрицательной полярностях;

затем находят искажающую матрицу магнитометра в собственной системе координат меры и_м по формуле:

и разлагают ее на произведение

и_м=L·I_м,

где L=(L₁ L₂ L₃) - диагональная матрица безразмерных коэффициентов, незначительно отличающихся от единицы;

- матрица единичных чувствительных осей магнитометра в виде вектор-строк I_iм в собственной системе координат меры;

с помощью чего определяют три неортогональности осей е_ij по формуле:

e_ij=(l_i1l_j1+l_i2l_j2+l_i3l_j3)_м,

по которым строят матрицу осей магнитометра I₀ в его собственной ортогональной системе координат, в которой определяют искажающую и₀=L·I₀ и находят искомую калибровочную матрицу .

П_в=[(В_uм1+-В_uм1-)/2В_м1(В_uм2+-В_uм2-)/2В_м2(В_uм3+-В_uм3-)/2В_м3],

где В_uмi+, B_uмi- - вектор-столбцы результатов измерения магнитной индукции на выходе магнитометра, соответствующие магнитному полю, воспроизводимому мерой при положительной и отрицательной полярностях;

затем находят калибровочную матрицу магнитометра в собственной системе координат меры k_м по формуле:

и разлагают ее на произведение

k_м=h_м·k,

где h_м=(h₁ h₂ h₃)_м - матрица единичных чувствительных осей магнитометра в виде вектор-столбцов h_iм в собственной системе координат меры;

k=(k₁ k₂ k₃)_м - диагональная матрица безразмерных коэффициентов k_i;

затем определяют три неортогональности осей магнитометра μ_ij по формуле:

μ_ij=(h_1ih_1j+h_2ih_2j+h_3ih_3j),

по которым строят матрицу осей магнитометра h₀ в его собственной ортогональной системе координат, в которой находят искомую калибровочную матрицу k_h0=h₀·k.