СПОСОБ НЕРАЗРУШАЮЩЕГО ОБЪЕМНОГО ИЗМЕРЕНИЯ ВЕКТОРНОЙ ФУНКЦИИ МАГНИТНОЙ ИНДУКЦИИ НЕОДНОРОДНО РАСПРЕДЕЛЕННОГО В ПРОСТРАНСТВЕ И ПЕРИОДИЧЕСКИ ИЗМЕНЯЮЩЕГОСЯ ВО ВРЕМЕНИ МАГНИТНОГО ПОЛЯ Российский патент 2013 года по МПК G01R33/00 

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

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

Известен способ сканирования магнитного поля, реализованный в устройстве [2], который основан на измерении магнитного поля поочередно в точках плоскости посредством последовательно соединенной матрицы чувствительных элементов, включающей N феррозондов, с последующей обработкой измеряемых напряжений с помощью ЭВМ. Однако известные способ и реализующее его устройство не позволяют выполнять измерения в местах недоступных для механического проникновения матрицы чувствительных элементов и позволяют получить только одну составляющую векторной функции магнитной индукции. Также отсутствует возможность выполнять измерения в произвольных точках, так как все точки измерения определены матричной структурой магниточувствительных элементов. Кроме того, способ имеет ограничения по количеству магниточувствительных элементов и, как следствие, по разрешающей способности.

Известен способ измерения и топографии магнитных полей рассеяния вблизи поверхности объекта, реализованный в устройстве [3], который основан на последовательном перемещении в соответствии с заданной траекторией с помощью измерительной штанги одного трехкомпонентного датчика Холла относительно объекта измерения посредством блока механических перемещений с поворотным столиком и подвижными каретками, приводимыми в движение шаговыми двигателями, с последующей статистической обработкой результатов измерения измерительно-вычислительным блоком. Однако известные способ и реализующее его устройство не позволяют выполнять измерения в местах недоступных для механического проникновения датчика Холла.

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

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

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

$$\{\begin{array}{l}{p}_x\left(s,\overline{n},t\right)=\left[sin\theta \cdot p\left(s,\overline{n},t\right)+cos\theta \cdot p\left(s,{\overline{n}}_{\perp \theta },t\right)\right]cos\alpha -\\ -\left[sin\theta \cdot p\left(s,{\overline{n}}_{\perp \alpha },t\right)+cos\theta \cdot p\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)\right]sin\alpha ;\\ p_y\left(s,\overline{n},t\right)=\left[sin\theta \cdot p\left(s,\overline{n},t\right)+cos\theta \cdot p\left(s,{\overline{n}}_{\perp \theta },t\right)\right]sin\alpha +\\ +\left[sin\theta \cdot p\left(s,{\overline{n}}_{\perp \alpha },t\right)+cos\theta \cdot p\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)\right]cos\alpha ;\\ p_z\left(s,\overline{n},t\right)=cos\theta \cdot p\left(s,\overline{n},t\right)-sin\theta \cdot p\left(s,{\overline{n}}_{\perp \theta },t\right),\end{array}$$

где $$p\left(s,\overline{n},t\right)$$ - исходные проекционные данные функции индукции; s - координата оси направления управляемого пространственного перемещения рабочего органа; θ и α - зенитный и азимутальный углы отклонения рабочего органа, соответственно; $${\overline{n}}_{\perp \theta }=\left(n,\theta +\pi /2,\alpha \right)$$ , $${\overline{n}}_{\perp \alpha }=\left(n,\theta ,\alpha +\pi /2\right)$$ и $${\overline{n}}_{\perp \theta ,\alpha }=\left(n,\theta +\pi /2,\alpha +\pi /2\right)$$ - векторы нормалей плоскостей, соответственно перпендикулярно повернутых по углам θ и α на величину равную π/2 относительно текущего вектора нормали рабочего органа $$\overline{n}$$ ; t - время; p_x, p_y, p_z - исходные проекционные данные x-, y-, z-компонент векторной функции магнитной индукции, которые записывают через магнитные потоки, численно равные интегралам по времени напряжений, индуцируемых в контуре:

$$\{\begin{array}{l}{p}_x\left(s,\overline{n},t\right)=-\left[sin\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt+cos\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt\right]cos\alpha +\\ +\left[sin\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \alpha },t\right)dt+cos\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)dt\right]sin\alpha ;\\ p_y\left(s,\overline{n},t\right)=-\left[sin\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt+cos\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt\right]sin\alpha -\\ -\left[sin\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \alpha },t\right)dt+cos\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)dt\right]cos\alpha ;\\ p_z\left(s,\overline{n},t\right)=-cos\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt+sin\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt,\end{array}$$

где u(t) - напряжения, индуцируемые изменением потока в контуре в момент времени t; T - период, в течение которого происходит одно полное изменение функции магнитной индукции; n=1, 2, … - номер периода, причем к исходным проекционным данным применяют обратное преобразование Радона, основанное на их свертке, осуществляющей фильтрацию с использованием свертывающей функции, являющейся обратным Фурье-преобразованием квадрата частоты пространственного спектра по формуле:

$$\{\begin{array}{l}{B}_x\left(x,y,z,t\right)={\displaystyle \int \left(p_x\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\overline{n}={\displaystyle \underset{0}{\overset{\pi /2}{\int }}\left[{\displaystyle \underset{0}{\overset{\pi }{\int }}\left(p_x\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\alpha }\right]cos\theta d\theta =}}\\ ={\displaystyle \underset{0}{\overset{\pi /2}{\int }}[{\displaystyle \underset{0}{\overset{\pi }{\int }}((-\left[sin\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt+cos\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt\right]cos\alpha +}}\\ +\left[sin\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \alpha },t\right)dt+cos\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)dt\right]sin\alpha )\ast h\left(s\right))d\alpha ]cos\theta d\theta ;\\ B_y\left(x,y,z,t\right)={\displaystyle \int \left(p_y\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\overline{n}={\displaystyle \underset{0}{\overset{\pi /2}{\int }}\left[{\displaystyle \underset{0}{\overset{\pi }{\int }}\left(p_y\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\alpha }\right]cos\theta d\theta =}}\\ ={\displaystyle \underset{0}{\overset{\pi /2}{\int }}[{\displaystyle \underset{0}{\overset{\pi }{\int }}((-\left[sin\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt+cos\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt\right]sin\alpha +}}\\ +\left[sin\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \alpha },t\right)dt+cos\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)dt\right]cos\alpha )\ast h\left(s\right))d\alpha ]cos\theta d\theta ;\\ B_z\left(x,y,z,t\right)={\displaystyle \int \left(p_z\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\overline{n}={\displaystyle \underset{0}{\overset{\pi /2}{\int }}\left[{\displaystyle \underset{0}{\overset{\pi }{\int }}\left(p_z\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\alpha }\right]cos\theta d\theta =}}\\ ={\displaystyle \underset{0}{\overset{\pi /2}{\int }}\left[{\displaystyle \underset{0}{\overset{\pi }{\int }}\left(-\left[cos\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt-sin\theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt\right]\ast h\left(s\right)\right)d\alpha }\right]cos\theta d\theta ,}\end{array}$$

где символ «*» есть оператор свертки; h(s) - свертывающая функция; $$d\overline{n}$$ - дифференциал вектора нормали; dθ и dα - дифференциалы зенитного и азимутального углов, соответственно; B_x(x,y,z,t₀), B_y(x,y,z,t₀), B_z(x,y,z,t₀) - x-, y-, z-компоненты векторной функции магнитной индукции $$\overline{B}$$ , соответственно, благодаря чему реконструируют декартовы компоненты распределения векторной функции магнитной индукции в пространстве.

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

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

Скалярное значение магнитного потока Ф_D исходной векторной функции магнитной индукции B(x,y,z,t) через поверхность D, образованную плоскостью контура, есть интеграл по плоскости, перпендикулярной вектору ее нормали $$\overline{n}$$ , и проходящей на расстоянии s от начала системы координат объема измерения:

$${Ф}_D={\displaystyle \underset{D\left(s,\overline{n}\right)}{\int }\overline{B}d\overline{\sigma },\left(1\right)}$$

где $$d\overline{\sigma }$$ - вектор нормали, выставленный к элементарной площадке плоскости интегрирования и численно равный ее площади. Переменные $$\overline{n}$$ и s задают положение контура в объеме измерения, причем в метрике сферической системы координат вектор нормали $$\overline{n}=\left(n,\theta ,\alpha \right)$$ .

Технически регистрация магнитного потока осуществляется благодаря закону Фарадея, согласно которому напряжение u(t), наводимое в контуре, определяется выражением:

$$u\left(t\right)=-\frac{d{Ф}_D\left(t\right)}{dt},\left(2\right)$$

где dФ_D(t) - дифференциал потока векторной функции индукции; dt - дифференциал времени.

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

Интегрирование по времени t выраженного из равенства (2) дифференциала dФ_D(t) приводит к записи значения магнитного потока:

$${Ф}_D\left(t\right)=-{\displaystyle \underset{nT}{\overset{nT+t}{\int }}u\left(t\right)dt,\left(3\right)}$$

где T - период, в течение которого происходит одно полное изменение магнитной индукции; n=1, 2, … - номер периода.

Измерение начинается с процедуры управляемого пространственного перемещения рабочего органа в объеме измерения с получением проекций, которая соответствует прямому преобразованию Радона для трехмерного пространства (фиг.1). Управляемое пространственное перемещение реализует способ параллельного формирования исходных проекционных данных, для которого рабочий орган совершает поступательно-поворотное движение, предполагающее чередование дискретных параллельных перемещений в направлении вектора нормали органа $$\overline{n}$$ по оси OS и поворотов направления этих перемещений, задаваемых зенитным θ и азимутальным α углами сферической системы координат. Таким образом, дискретные параллельные перемещения многократно повторяются под разными углами, причем для зенитного в интервале от 0 до ½π, а для азимутального в интервале от 0 до π. В ходе управляемого пространственного перемещения рабочего органа после каждого его дискретного передвижения вдоль оси OS для текущего значения переменной положения s под углами, определенными вектором нормали $$\overline{n}$$ , регистрируется плоскостная проекция p - значение интеграла исходной функции индукции В по плоскости D, численно равное с учетом (1) значению пронизывающего магнитного потока Ф_D(t):

$$p\left(s,\overline{n},t\right)={\displaystyle \underset{D\left(s,\overline{n}\right)}{\int }\overline{B}d\overline{\sigma }={Ф}_D\left(t\right).\left(4\right)}$$

Принимая во внимание выражение (3), на основании (4) проекция значений интеграла исходной функции индукции $$p\left(s,\overline{n},t\right)$$ с учетом введенной зависимости от переменной времени записывается посредством интегралов по времени напряжений, индуцируемых изменением магнитного потока в соответствии с законом электромагнитной индукции Фарадея в контуре:

$$p\left(s,\overline{n},t\right)={\displaystyle \underset{D\left(s,\overline{n}\right)}{\int }Bd\overline{\sigma }=-{\displaystyle \underset{nT}{\overset{nT+t}{\int }}u\left(s,\overline{n},t\right)dt.\left(5\right)}}$$

Уравнение (5) приведено для сферической системы координат, при этом в декартовой системе координат (X,Y,Z) выражение компонент данного равенства имеет вид:

$$\{\begin{array}{l}{p}_x\left(s,\overline{n},t\right)=\left[\sin \theta \cdot p\left(s,\overline{n},t\right)+\cos \theta \cdot p\left(s,{\overline{n}}_{\perp \theta },t\right)\right]\cos \alpha -\\ -\left[\sin \theta \cdot p\left(s,{\overline{n}}_{\perp \alpha },t\right)+\cos \theta \cdot p\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)\right]\sin \alpha ;\\ p_y\left(s,\overline{n},t\right)=\left[\sin \theta \cdot p\left(s,\overline{n},t\right)+\cos \theta \cdot p\left(s,{\overline{n}}_{\perp \theta },t\right)\right]\sin \alpha +\\ +\left[\sin \theta \cdot p\left(s,{\overline{n}}_{\perp \alpha },t\right)+\cos \theta \cdot p\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)\right]\cos \alpha ;\\ p_z\left(s,\overline{n},t\right)=\cos \theta \cdot p\left(s,\overline{n},t\right)-\sin \theta \cdot p\left(s,{\overline{n}}_{\perp \theta },t\right),\end{array}\left(6\right)$$

где $${\overline{n}}_{\perp \theta }=\left(n,\theta +\pi /2,\alpha \right)$$ , $${\overline{n}}_{\perp \alpha }=\left(n,\theta ,\alpha +\pi /2\right)$$ и $${\overline{n}}_{\perp \theta ,\alpha }=\left(n,\theta +\pi /2,\alpha +\pi /2\right)$$ - векторы нормалей плоскостей, соответственно перпендикулярно повернутых по углам θ и α на величину равную π/2 относительно текущего вектора нормали $$\overline{n}$$ .

С учетом выражения (5), система уравнений (6) записывается в виде:

$$\{\begin{array}{l}{p}_x\left(s,\overline{n},t\right)=-\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt\right]\cos \alpha +\\ +\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \alpha },t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)dt\right]\sin \alpha ;\\ p_y\left(s,\overline{n},t\right)=-\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt\right]\sin \alpha -\\ -\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \alpha },t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)dt\right]\cos \alpha ;\\ p_z\left(s,\overline{n},t\right)=-\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt+\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt.\end{array}\left(7\right)$$

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

Задача реконструкции декартовых компонент в способе неразрушающего объемного измерения сводится к решению основного интегрального уравнения (7) с нахождением распределения компонент B_x, B_y, B_z по измеренным значениям плоскостных проекций p_x, p_y, p_z. Решение основного интегрального уравнения (7) предполагает использование алгоритма реконструкции, основанного на обратном преобразовании Радона, посредством обратной проекции с Фурье-фильтрацией. Фурье-фильтрация осуществляется путем свертки проекций компонент непосредственно в пространстве оригинала Фурье-преобразования с соответствующей реализующей фильтрацию пространстве оригинала Фурье-преобразования с соответствующей реализующей фильтрацию свертывающей функцией h(s), являющейся обратным Фурье-преобразованием F^-1[] квадрата частоты K пространственного спектра [6]. Так из (7) получают:

$$\{\begin{array}{l}{B}_x\left(x,y,z,t\right)={\displaystyle \int \left(p_x\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\overline{n}={\displaystyle \underset{0}{\overset{\pi /2}{\int }}\left[{\displaystyle \underset{0}{\overset{\pi }{\int }}\left(p_x\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\alpha }\right]\cos \theta d\theta =}}\\ ={\displaystyle \underset{0}{\overset{\pi /2}{\int }}[{\displaystyle \underset{0}{\overset{\pi }{\int }}((-\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt\right]\cos \alpha +}}\\ +\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \alpha },t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)dt\right]\sin \alpha )\ast h\left(s\right))d\alpha ]\cos \theta d\theta ;\\ B_y\left(x,y,z,t\right)={\displaystyle \int \left(p_y\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\overline{n}={\displaystyle \underset{0}{\overset{\pi /2}{\int }}\left[{\displaystyle \underset{0}{\overset{\pi }{\int }}\left(p_y\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\alpha }\right]\cos \theta d\theta =}}\\ ={\displaystyle \underset{0}{\overset{\pi /2}{\int }}[{\displaystyle \underset{0}{\overset{\pi }{\int }}((-\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt\right]\sin \alpha +}}\\ +\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \alpha },t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)dt\right]\cos \alpha )\ast h\left(s\right))d\alpha ]\cos \theta d\theta ;\\ B_z\left(x,y,z,t\right)={\displaystyle \int \left(p_z\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\overline{n}={\displaystyle \underset{0}{\overset{\pi /2}{\int }}\left[{\displaystyle \underset{0}{\overset{\pi }{\int }}\left(p_z\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\alpha }\right]\cos \theta d\theta =}}\\ ={\displaystyle \underset{0}{\overset{\pi /2}{\int }}\left[{\displaystyle \underset{0}{\overset{\pi }{\int }}\left(-\left[\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt-\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt\right]\ast h\left(s\right)\right)d\alpha }\right]\cos \theta d\theta ,}\end{array}\left(8\right)$$

где символ «*» есть оператор свертки; $$h(s)=F^{-1}\lfloor K^2\rfloor ={\displaystyle \int K^2\cdot e^{2\pi Ks}dK}$$ , $${\displaystyle \int d\overline{n}={\displaystyle \int {\displaystyle \int d\alpha cos\theta d\theta }}}$$ - дифференциал вектора нормали; dθ и dα - дифференциал зенитного и азимутального углов, соответственно; B_x(x,y,z,t), B_y(x,y,z,t), B_z(x,y,z,t) - x-,y-, z-компоненты векторной функции магнитной индукции $$\overline{B}$$ , соответственно.

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

Литература

1. Авторское свидетельство СССР №1652951, кл. G01R 33/02, опубл. 30.05.1991.

2. Авторское свидетельство СССР №1762282, кл. G01R 33/02, опубл. 15.09.1992.

3. Авторское свидетельство СССР №1684761, кл. G01R 33/06, опубл. 15.10.1991.

4. Патент РФ №2179323, кл. G01R 33/02, опубл. 10.02.2002.

5. Патент РФ №2174235, кл. G01R 33/02, опубл. 27.09.2001.

6. J. Radon. Uber die Bestimmung von Funktionen durch ihre Integralwerte langs gewisser Mannigfaltigkeiten // Berichte Sachsische Akademie der Wissenschaften, Bande 29, s.262-277, Leipzig, 1917.

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

Способ неразрушающего объемного измерения векторной функции магнитной индукции неоднородно распределенного в пространстве и периодически изменяющегося во времени магнитного поля, заключающийся в измерениях распределения векторной функции магнитной индукции периодически изменяющегося во времени поля в определенных точках исследуемого пространства для произвольно выбранных моментов времени на периоде и основанный на последовательно-поступательных перемещениях и поворотах на углы рабочего магнитоизмерительного органа и регистрации индуцируемых в нем напряжений, отличающийся тем, что измерения мгновенных объемных состояний распределения неоднородного в пространстве магнитного поля осуществляют в местах, недоступных для механического проникновения, а рабочий магнитоизмерительный орган выполняют в виде ориентируемого в пространстве одного контура, привязанного к сферической системе координат, причем в ходе его управляемого пространственного перемещения осуществляют способ параллельного формирования исходных проекционных данных функции индукции, для чего рабочим органом совершают движение, предполагающее чередование дискретных параллельных перемещений в направлении оси вектора нормали органа и поворотов направления этих перемещений, задаваемых зенитным и азимутальным углами сферической системы координат, таким образом, что дискретные параллельные перемещения многократно повторяют под разными углами, причем для зенитного в интервале от 0 до $$\frac{1}{2}\pi$$ , а для азимутального в интервале от 0 до π, а необходимые для алгоритма реконструкции исходные проекционные данные декартовых компонент распределения векторной функции магнитной индукции в объеме измерения в декартовой системе координат получают посредством тригонометрических преобразований:
$$\{\begin{array}{l}{p}_x\left(s,\overline{n},t\right)=\left[\sin \theta \cdot p\left(s,\overline{n},t\right)+\cos \theta \cdot p\left(s,{\overline{n}}_{\perp \theta },t\right)\right]\cos \alpha -\\ -\left[\sin \theta \cdot p\left(s,{\overline{n}}_{\perp \alpha },t\right)+\cos \theta \cdot p\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)\right]\sin \alpha ;\\ p_y\left(s,\overline{n},t\right)=\left[\sin \theta \cdot p\left(s,\overline{n},t\right)+\cos \theta \cdot p\left(s,{\overline{n}}_{\perp \theta },t\right)\right]\sin \alpha +\\ +\left[\sin \theta \cdot p\left(s,{\overline{n}}_{\perp \alpha },t\right)+\cos \theta \cdot p\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)\right]\cos \alpha ;\\ p_z\left(s,\overline{n},t\right)=\cos \theta \cdot p\left(s,\overline{n},t\right)-\sin \theta \cdot p\left(s,{\overline{n}}_{\perp \theta },t\right),\end{array}$$
где $$p\left(s,\overline{n},t\right)$$ - исходные проекционные данные функции индукции;
s - координата оси направления управляемого пространственного перемещения рабочего органа;
θ и α - зенитный и азимутальный углы отклонения рабочего органа, соответственно;
$${\overline{n}}_{\perp \theta }=\left(n,\theta +\pi /2,\alpha \right)$$ , $${\overline{n}}_{\perp \alpha }=\left(n,\theta ,\alpha +\pi /2\right)$$ и $${\overline{n}}_{\perp \theta ,\alpha }=\left(n,\theta +\pi /2,\alpha +\pi /2\right)$$ - векторы нормалей плоскостей, соответственно перпендикулярно повернутых по углам θ и α на величину, равную π/2 относительно текущего вектора нормали рабочего органа $$\overline{n}$$ ;
t - время;
p_x, p_y, p_z - исходные проекционные данные x-, y-, z-компонент векторной функции магнитной индукции, которые записывают через магнитные потоки, численно равные интегралам по времени напряжений, индуцируемых в контуре:
$$\{\begin{array}{l}{p}_x\left(s,\overline{n},t\right)=-\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt\right]\cos \alpha +\\ +\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \alpha },t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)dt\right]\sin \alpha ;\\ p_y\left(s,\overline{n},t\right)=-\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt\right]\sin \alpha -\\ -\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \alpha },t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)dt\right]\cos \alpha ;\\ p_z\left(s,\overline{n},t\right)=-\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt+\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt,\end{array}$$
где u(t) - напряжения, индуцируемые изменением потока в контуре в момент времени t;
T - период, в течение которого происходит одно полное изменение функции магнитной индукции;
n=1, 2, … - номер периода, причем к исходным проекционным данным применяют обратное преобразование Радона, основанное на их свертке, осуществляющей фильтрацию с использованием свертывающей функции, являющейся обратным Фурье-преобразованием квадрата частоты пространственного спектра по формуле
$$\{\begin{array}{l}{B}_x\left(x,y,z,t\right)={\displaystyle \int \left(p_x\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\overline{n}={\displaystyle \underset{0}{\overset{\pi /2}{\int }}\left[{\displaystyle \underset{0}{\overset{\pi }{\int }}\left(p_x\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\alpha }\right]\cos \theta d\theta =}}\\ ={\displaystyle \underset{0}{\overset{\pi /2}{\int }}[{\displaystyle \underset{0}{\overset{\pi }{\int }}((-\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt\right]\cos \alpha +}}\\ +\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \alpha },t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)dt\right]\sin \alpha )\ast h\left(s\right))d\alpha ]\cos \theta d\theta ;\\ B_y\left(x,y,z,t\right)={\displaystyle \int \left(p_y\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\overline{n}={\displaystyle \underset{0}{\overset{\pi /2}{\int }}\left[{\displaystyle \underset{0}{\overset{\pi }{\int }}\left(p_y\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\alpha }\right]\cos \theta d\theta =}}\\ ={\displaystyle \underset{0}{\overset{\pi /2}{\int }}[{\displaystyle \underset{0}{\overset{\pi }{\int }}((-\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt\right]\sin \alpha +}}\\ +\left[\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \alpha },t\right)dt+\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta ,\alpha },t\right)dt\right]\cos \alpha )\ast h\left(s\right))d\alpha ]\cos \theta d\theta ;\\ B_z\left(x,y,z,t\right)={\displaystyle \int \left(p_z\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\overline{n}={\displaystyle \underset{0}{\overset{\pi /2}{\int }}\left[{\displaystyle \underset{0}{\overset{\pi }{\int }}\left(p_z\left(s,\overline{n},t\right)\ast h\left(s\right)\right)d\alpha }\right]\cos \theta d\theta =}}\\ ={\displaystyle \underset{0}{\overset{\pi /2}{\int }}\left[{\displaystyle \underset{0}{\overset{\pi }{\int }}\left(-\left[\cos \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,\overline{n},t\right)dt-\sin \theta \cdot {\displaystyle \underset{nT}{\overset{nT+t}{\int }}u}\left(s,{\overline{n}}_{\perp \theta },t\right)dt\right]\ast h\left(s\right)\right)d\alpha }\right]\cos \theta d\theta ,}\end{array}$$
где символ «*» есть оператор свертки;
h(s) - свертывающая функция;
$$d\overline{n}$$ - дифференциал вектора нормали;
dθ и dα - дифференциалы зенитного и азимутального углов соответственно;
B_x(x,y,z,t₀), B_y(x,y,z,t₀), B_z(x,y,z,t₀), где x, y, z - компоненты векторной функции магнитной индукции В соответственно, благодаря чему реконструируют декартовы компоненты распределения векторной функции магнитной индукции в пространстве.