способ локализации магнитного диполя

Формула изобретения

Способ локализации магнитного диполя, включающий предварительные измерения расстояний между датчиками трехкомпонентных магнитометров, затем измерения с помощью этих датчиков градиента индукции магнитного поля, определение собственных чисел градиента _{2 1 3}, которыми определяют скалярную интенсивность его , а также определение соответствующих собственных векторов u₁, u₂, u₃, по которым находят четыре единичных направления на диполь n_i и направления его магнитного момента m_i, где i=1, 2, 3, 4 - индекс пары решений (одно истинное, три ложных), отличающийся тем, что с помощью датчиков трехкомпонентных магнитометров, размещенных у одного из концов переносной штанги на полуосях ее системы координат, соосно и на одинаковом удалении от ее центра, измеряют пространственные приращения индукции магнитного поля относительно опорного датчика на другом конце штанги, по которым вычисляют среднее приращение по всем полуосям, а также разности приращений для каждой одноименной пары полуосей, строят симметричную матрицу разностей, при этом из каждого диагонального элемента матрицы вычитают треть ее следа, а единичные направления находят по формулам, зависящим от размера числа = ₁/N:

при -0,95 <0,95:



при -1 <-0,95:



при 0,95< 1:



после чего, нормируя вектор среднего приращения, получают векторное уравнение с одним неизвестным k=а/r:



где В₀₇ - вектор среднего приращения;

а - расстояние между центром опорного датчика и центром системы координат;

=(1+k²-2kn_a·n)^1/2 ;

n_a - единичный вектор направления штанги,

в свою очередь, векторное нормированное уравнение сводят к модульному |В₀₇|/N_a=f(k), решая которое для каждой из четырех пар найденных векторов n_i и m_i получают соответствующее значение k_i , затем каждую пару величин n_i, m_i, если ей соответствует физически приемлемое значение k_i, подставляют в векторное нормированное уравнение, рассчитывают его и отбирают тройку истинных решений n, m и k по минимуму модуля разности рассчитанной и измеренной версий этого уравнения, после чего определяют искомое расстояние до источника r=а/k.

Описание изобретения к патенту

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

Многие сосредоточенные источники магнитного поля можно аппроксимировать (точечным) диполем с магнитным моментом М=Mm, где М - модуль магнитного момента, m - единичный вектор направления магнитного момента. Различают задачи обнаружения, пеленгации и локализации (неподвижного или движущегося) диполя по параметрам его магнитного поля, измеренным в одной или нескольких точках наблюдения. При локализации определяют радиус-вектор X=Xn, где X - модуль расстояния, n - единичный вектор направления. Другими словами, при локализации определяют неизвестные координаты вектора X в системе координат наблюдения. При пеленгации определяют единичный вектор направления n или направляющие косинусы углов его в системе координат наблюдения. При обнаружении определяют параметры неоднородного магнитного поля на фоне поля Земли, по которым судят о том, что «где-то что-то есть». Некоторые специалисты под обнаружением понимают и пеленгацию и локализацию. Однако мы для определенности будем различать эти задачи.

Задачу пеленгации магнитного диполя, по-видимому, впервые поставил О.П.Хвостов [Магнитный пеленгатор. // Геофизическое приборостроение. - Л.: Гостехиздат. 1963. С.3-11], который также предложил приближенное решение ее по градиенту модуля дипольной индукции.

Первое аналитическое (явное) решение задачи пеленгации диполя по градиенту вектора дипольной индукции получил C.Frahm [W.M.Wynn а.о. Advanced Superconducting Gradiometer. / Magnetometer Arrays and a Novel Signal Processing Technique. // IEEE Transactions on Magnetics, vol. MAG-11, no.2, March 1975 Р.701-707]. Решение Frahm a было неоднозначным. Оно включало 4 направления, одно из которых было истинным, остальные три - ложными. Вместе с тем решение Frahm a осталось нераскрытым, поскольку не было найдено удобной формы представления этого решения.

Затем T.Kolbenhayer [Prispevok krieseniu obratenej magnetometrickey ulohy predipol. // Banickey listy. 1976. #4 31-35 (Реферативный Журнал Геология. 1980. № 12. C.20)] нашел явное решение задачи локализации по заданным вектору и градиенту вектора дипольной индукции, опираясь на теорему Эйлера для однородных функций. Однако решение Kolbenhayer a оказалось неустойчивым к ориентации искомого диполя.

Позднее оба упомянутых решения были радикально усовершенствованы в работе [В.Г.Семенов. Решение обратной задачи для источника физического поля дипольной или квадрупольной модели. // Методы и средства точных магнитных измерений. / Сб. научных трудов НПО «ВНИИМ им. Д.И.Менделеева». 1980. С.3-19]. В результате задача локализации обрела устойчивость, а для задачи пеленгации было найдена простая и удачная форма выражения этого решения. В последующие годы работы Frahm a [1975] и Семенова [1980] подверглись несправедливой критике в части неоднозначности решения в статье [В.В.Тихонов и Б.Н.Шалаев Особенности решения обратных задач при электромагнитном исследовании подводных сред. Известия ЛЭТИ им.Ульянова (Ленина). Выпуск 327 Л-д 1983 г. Стр.28] и монографии [под редакцией лауреата Гос. премии СССР Е.Г.Пащенко «Оптико-физические средства исследования океана». Л-д. Судостроение 1984. Стр.207]. Авторы этих публикаций предложили якобы однозначное решение задачи вместо решения [Семенов 1980]. Ответная статья [В.Г.Семенов и А.Я.Зайончковский. Сравнительный анализ градиентных методов локализации источников магнитного поля дипольной модели. Сб. научных трудов «Создание средств измерений для метрологического обеспечения прецизионных нанотесламетров». ВНИИМ им. Менделеева Л-д. 1988 г. С.56-67] показала несостоятельность претензий на однозначность решения Тихонова и Шалаева. Вместе с тем, Семенов и Зайончковский [1988] получили аналитическое выражение погрешности способа Frahm a [1975]/Семенова [1980], выявившее некоторую неустойчивость этого способа, который ранее считался полностью устойчивым. Много лет спустя это аналитическое выражение было перепроверено компьютерным моделированием [Ю.М.Иванов и В.Г.Семенов. Оценки погрешностей пеленгации и локализации источника магнитного поля дипольной модели. // Измерительная техника № 5. 2008. С.30-35 и № 7. 2008. С.72]. Результаты компьютерного моделирования подтвердили существование наличие области неустойчивости.

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

а) способ локализации [Семенов, 1980]. Этот способ устойчив и однозначен, но он не позволяет осуществлять локализацию в движении носителя средств измерения;

б) способ пеленгации [Семенов и Зайончковский, 1988]. Этот способ позволяет выполнять измерения при движении средств измерения, но он неоднозначен, так как дает четыре решения (какое из них истинное неизвестно), к тому же в некоторой узкой области этот способ неустойчив. Но самое главное, это способ пеленгации, а не локализации.

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

- Патент РФ № 2256930 Способ минимизации погрешностей в переносном дифференциальном магнитометре. Приоритет изобретения от 07.04.2003 г.;

- Патент РФ № 2257593 Способ минимизации погрешностей в дифференциальном магнитометре. Приоритет изобретения от 07.04.2003 г.;

- Патент РФ № 2290657 Способ определения местоположения и ориентации трехкомпонентного датчика магнитометра относительно меры магнитного момента. Приоритет изобретения от 09.08.2005 г.;

- Патент РФ № 2313801 Способ измерения приращения индукции магнитного поля трехкомпонентным дифференциальным магнитометром. Приоритет изобретения от 04.05.2006 г.

В качестве прототипа заявляемого способа выбран способ пеленгации источника магнитного поля дипольной модели [Семенов и Зайончковский, 1988] (далее, просто [1988]) по совокупности существенных признаков, наиболее близких к заявляемому способу.

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

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

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

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

Указанный технический результат достигается тем, что заявляемый способ локализации магнитного диполя, включающий предварительные измерения расстояний между датчиками трехкомпонентных магнитометров, затем измерения с помощью этих датчиков градиента индукции магнитного поля, определение собственных чисел градиента _{2 1 3}, которыми определяют скалярную интенсивность его , а также определение соответствующих собственных векторов u₁, u₂, u₃, по которым находят четыре единичных направления на диполь n_i и направления его магнитного момента m_i по формулам

где i=1, 2, 3, 4 - индекс пары решений (одно истинное, три ложных), отличается тем, что с помощью датчиков трехкомпонентных магнитометров, размещенных у одного из концов переносной штанги на полуосях ее системы координат, соосно и на одинаковом удалении от ее центра, измеряют пространственные приращения индукции магнитного поля относительно опорного датчика на другом конце штанги, по которым вычисляют среднее приращение по всем полуосям, а также разности приращений для каждой одноименной пары полуосей, строят симметричную матрицу разностей, при этом из каждого диагонального элемента матрицы вычитают треть ее следа, а единичные направления находят по формулам, зависящим от размера числа = ₁/N:

при -0.95 <0.95:

при -1 <-0,95:

при 0,95< 1:

после чего, нормируя вектор среднего приращения, получают векторное уравнение с одним неизвестным k=a/r:

где В₀₇ - вектор среднего приращения;

a - расстояние между центром опорного датчика и центром системы координат;

=(1+k²-2kn_a·n)^1/2 ;

n_a - единичный вектор направления штанги,

в свою очередь, векторное нормированное уравнение сводят к модульному |B₀₇|/Na=f(k), решая которое для каждой из четырех пар найденных векторов n _i и m_i, получают соответствующее значение k _i, затем каждую пару величин n_i, m_i , если ей соответствует физически приемлемое значение k_i , подставляют в векторное нормированное уравнение, рассчитывают его и отбирают тройку истинных решений n, m и k по минимуму модуля разности рассчитанной и измеренной версий этого уравнения, после чего определяют искомое расстояние до источника r=a/k.

На Фиг.1 изображена схема реализации заявляемого способа.

Схема, изображенная на Фиг.1, содержит датчики трехкомпонентных магнитометров, установленные на переносной штанге. Датчики 1-6 расположены у одного конца штанги, опорный датчик 7 - у другого конца ее. Датчики 1-2 размещены на полуосях х₊ и х_- системы координат (СК) штанги, на одинаковом расстоянии от центра 0 СК, датчики 3-4 - на полуосях y₊ и y_-, датчики 5-6 - на полуосях z ₊ и z_-. Все датчики 1-7 сориентированы своими осями вдоль осей СК штанги.

Заявляемый способ действует следующим образом. Датчики 1-6 включены по схеме измерения приращений магнитной индукции относительно опорного датчика 7. Это означает, что при измерении вектора приращения магнитной индукции, например, парой 1-7, регистрирующие устройства отсчитывают три проекции разности В₁₇=В₁-В₇ , но не сами векторы В₁ и В₇. Предварительно измеряют расстояния между датчиками, а также их взаимную ориентацию с помощью внешней трехкомпонентной меры магнитного момента (ММ). Для этого используют способ определения местоположения и ориентации трехкомпонентного датчика магнитометра относительно меры ММ [Патент РФ № 2290657]. При работе оперируют следующими результатами измерения: приращениями магнитной индукции В₁₇, В ₂₇, В₃₇, В₄₇, В₅₇, В ₆₇ и заранее определенными расстояниями между соответствующими датчиками a₁₂, a₃₄, a ₅₆, а также расстоянием между центрами СК 0 и датчика 7 - a.

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

Сначала с помощью измеренных приращений строят матрицу 3×3:

с помощью которой находят матрицу градиента

где А - транспонированная матрица А,

A₁₁ +A₂₂+A₃₃ - след матрицы А,

I - единичная матрица 3×3.

Операции (2) придают матрице градиента симметрию и обнуляют его след, поскольку известно, что градиент вектора магнитной индукции симметричен и бесследен. Затем стандартными приемами по матрице градиента находят его собственные числа _{2 1 3} и векторы u₁, u₂, u ₃, которыми определяют скалярную интенсивность градиента N, число , которое является скалярным произведением единичных векторов направлений на источник из центра СК n и магнитного момента m, а также сами направления n и m по формулам, зависящим от размера :

где µ₀=4 ×10^-7 Гн/м - магнитная постоянная;

М - размер магнитного момента источника;

r - расстояние между источником и центром СК;

(Обратим внимание на то, что сами величины М и r по градиенту не определяются);

при -0.95 <0.95:

при -1 <-0,95:

при 0.95< 1:

где i=1, 2, 3, 4.

Формулы (5)-(7) дают неоднозначное решение, т.к они определяют 4 пары решений, одна пара - истинное решение, три - ложные. Какая пара истинная, неизвестно, так как все решения удовлетворяют градиенту. Проявляется это в том, что собственные векторы определяются с точностью до знака.

Обратим внимание, что формулы (5)-(7) применяют в зависимости от размера . Анализ погрешностей показал, что первые две строчки (5), каждая из двух слагаемых, методически безупречны, но при наличии погрешностей одно слагаемое расходится при -1, а другое при +1. Поэтому выгодно, при | | 1÷0,95 использовать одно из этих слагаемых, допуская незначительную методическую погрешность, но повышая устойчивость к погрешностям входных данных. Таким образом, разбивка (5) на подвыражения, зависящие от размера , устраняет неустойчивость, имеющую место в известном способе [1988].

Далее определяют вектор среднего приращения

Равенство (8) основано на общеизвестной теореме Гаусса о среднем. Для источника дипольной модели вектор среднего приращения можно записать как

В правой части (9) помимо неизвестного k=a/r фигурирует еще два неизвестных: М и r. Для того чтобы избавиться от этих двух неизвестных, используют отношение

Выражение (10) представляет собой векторное нелинейное уравнение относительно только одного скалярного неизвестного k, учитывая, что величины , n_i, m_i известны из вышеприведенных решений, а вектор n_a известен как орт СК вдоль продольной оси штанги. Далее из (10) составляют модульное уравнение вида

где _a=m·n_a,

решая которое для каждой пары n_i, m_i, получают соответствующее k_i. Иногда полученное k_i не является физически приемлемым решением, например слишком большим, (NaN - термин МатЛаб а), такие решения отбрасывают. Каждую тройку решений n _i, m_i, k_i, где k_i физически приемлемое решение, подставляют в уравнение (10) и рассчитывают его. Истинную тройку решений определяют по наименьшему модулю разности _i измеренной и рассчитанной версий уравнения (10).

Таким образом в заявляемом способе устраняется неоднозначность решений, присущая известному способу [1988].

Затем по найденному истинному k определяют искомое расстояние до источника r=a/k, что завершает переход от (устойчивой и однозначной) пеленгации к (устойчивой и однозначной) локализации источника.

Заявляемый способ проверен с помощью компьютерного моделирования в пакете МатЛаб. По МатЛаб у рассчитывают компоненты градиента, к ним добавляют случайные погрешности с нулевым средним и заданным уровнем СКО компоненты , и осуществляют пеленгацию по известному способу [1988] 500 раз для каждого из 21 значения =-1, -0.9, -0.8 , 1. Всю совокупность обрабатывают для определения СКО пеленгации _n.

На Фиг.2 кривая 1 соответствует СКО пеленгации по уравнению (5), кривая 2 - по (6), 3 - по (7). Таким образом, при только одно из слагаемых (5) содержит почти всю полезную информацию, а другое слагаемое (5) - только «шумит», поэтому его отбрасывают. Отбрасывание второго слагаемого (5) создает методическую погрешность пеленгации, которая в диапазоне 0.95 | | 1 еще не велика и ею можно пренебречь в этой части диапазона. Кривая 2 представляет СКО первого члена (5), кривая 3 - СКО второго члена (5), кривая 4 - сумму методической погрешности плюс удвоенное СКО первого члена (5), кривая 5 - сумму методической погрешности плюс удвоенное СКО второго члена (5). Как показано на Фиг.2, при -1 -0.95 пеленгация по методически неточному уравнению (6) выгоднее, чем по точному уравнению (5), а при 0.95 1 по уравнению (7) выгоднее, чем по (5). Значение определяют в начале пеленгации по уравнению (4).

Пример, представленный на Фиг.2, соответствует погрешностям пеленгации диполя с магнитным моментом М=0.8 Ам² на расстоянии r=2.5 м от центра 0 СК штанги при расстоянии между датчиками 1-6 1=0.75 м. Уровень СКО шумов компоненты датчика _B=10^-10 Тл.

Аналитическое выражение СКО пеленгации (определения направления на искомый источник) по уравнению (5) имеет следующий вид:

Как видно из кривой 6 (гладкая линия) на Фиг.2, аналитическое выражение (12) практически полностью совпадает с результатами компьютерного моделирования. Это совпадение доказывает, что выражение (12) выведено правильно, а также что моделирование проведено правильно.

Таким образом, разбивка по заявляемому способу уравнения (5) на два дополнительных (6) и (7) снижает неустойчивость пеленгации, присущую известному способу [1988] при 0.95 | | 1.

Далее в Таблице 1 представлены примеры устранения по заявляемому способу неоднозначности пеленгации, присущей известному способу [1988]. Эти примеры подготовлены также компьютерным моделированием пеленгации/локализации источника, удаленного на расстояние r=3.5 м от центра 0 СК штанги, у которой опорный датчик 7 расположен на расстоянии а=1.5 м от центра 0.

Таблица 1
Определение параметра ki (11) и соответствующего ему - модуля разности измеренной и рассчитанной частей уравнения (10) для четырех пар решений (5), определяемых знаками собственных векторов u2, u3
	+u2, +u3	-u2, +u3	-u2, -u3	+u2, -u3
k1	1, нТл	k2	2, нТл	k3	3, нТл	k4	4, нТл
-0.9	0.4286	0.000	0.2573	21.4	0.2619	24.8	0.4574	15.2
-0.7	0.4286	0.000	0.2328	11.5	0.2499	21.1	0.0411	13.4
-0.5	0.4286	0.000	0.1927	7.0	0.2286	17.8	0.0678	17.4
-0.3	0.4286	0.000	0.1285	12.4	0.1885	15.3	0.0751	20.4
-0.1	0.4286	0.000	0.0415	18.5	0.1066	13.8	0.0403	21.2
0	0.4286	0.000	NaN*	-	NaN	-	NaN	-
+0.1	0.4286	0.000	0.0403	21.2	0.1066	13.8	0.0415	18.5
+0.3	0.4286	0.000	0.0751	20.4	0.1885	15.3	0.1285	12.4
+0.5	0.4286	0.000	0.0678	17.4	0.2286	17.8	0.1927	7.0
+0.7	0.4286	0.000	0.0411	13.4	0.2499	21.1	0.2328	11.5
+0.9	0.4286	0.000	0.4574	15.2	0.2619	24.8	0.2573	21.4
*NaN - not a number (термин МатЛаб а).

Как видно из каждой строки Таблицы 1, значение модуля разности ₁ - наименьшее, поэтому тройку решений n ₁, m₁, k₁ в данном примере принимают за истинные направления n, m и истинное отношение k=a/r. Остальные тройки отбрасывают как ложные.

Таким образом, заявляемый способ устраняет неоднозначность решений, присущую известному способу [1988].

Наконец, по найденному и отобранному истинному значению отношения k=a /r определяют искомое расстояние до источника r=a/k.

Для примера в Таблице 1 r=1.5/0.4286=3.4998 м. Отклонение от расчетного расстояния 3.5 м произошло из-за округления числа 0.4285714286 до 0.4286.

Таким образом, заявляемый способ осуществляет не только однозначную (и устойчивую) пеленгацию, но однозначную (и устойчивую) локализацию источника, в то время как известный способ [1988] выполняет только (неоднозначную и неустойчивую) пеленгацию источника.

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