Институт прикладной физики РАН
ДИАГНОСТИКА ПОДПОВЕРХНОСТНЫХ НЕОДНОРОДНОСТЕЙ СРЕДЫ КВАЗИСТАТИЧЕСКИМИ ЭЛЕКТРИЧЕСКИМИ ПОЛЯМИ РЕЗОНАНСНЫХ СВЧ-ЗОНДОВ
Д.В. Янин, А.Г. Галка,
А.И. Смирнов, А.В. Костров
Аннотация
Работа посвящена развитию метода диагностики подповерхностных неоднородностей диэлектрических сред квазистатическими электрическими полями резонансных СВЧ-зондов. В рамках данного метода источником информации о среде является импеданс электрически малых антенн, расположенных над поверхностью. Для модельных зондов (отрезок двухпроводной линии, краевая емкость цилиндрического конденсатора и т.п.) в рамках борновского приближения разработана схема решения обратной задачи по восстановлению пространственного распределения диэлектрической проницаемости.
Ключевые слова: квазистатическое поле, неоднородности среды, комплексная диэлектрическая проницаемость, электрически-малая антенна, обратная задача.
Abstract
The article is devoted to development of a method of diagnostics of subsurface inhomogeneity of dielectric medium by quasistatic electric fields of resonance microwave probes. Within the limits of this method the information source of mediums characteristic is the impedance of electrically-small antennas located over the mediums surface. The circuit of the decision of an inverse problem is developed for modelling probes (a cut of a two-wire line, edge capacity of the cylindrical capacitor, etc.) within the limits of Born approach on restoration of spatial distribution of a dielectric permittivity.
Key words: quasistatic electric field, medium inhomogeneity, complex dielectric permittivity, electrically-small antennas, inverse problem.
1. Введение
Дистанционная (бесконтактная) диагностика подповерхностных неоднородностей актуальна в различных областях науки и техники (например, в медицине, материаловедении, дефектоскопии, геологоразведке и т.д.). Одним из перспективных методов такой диагностики является резонансная ближнепольная СВЧ-томография, позволяющая визуализировать пространственное распределение электродинамических характеристик (диэлектрической проницаемости и проводимости) исследуемых объектов [1-4]. Основные принципы, на которых базируется этот метод диагностики, можно пояснить следующим образом.
Область среды, находящаяся в ближнем поле электрически малой антенны, оказывает влияние на ее импеданс. Если антенна включена в качестве нагрузки в резонансную систему, то по смещению резонансной частоты и изменению добротности можно судить об электродинамических параметрах среды вблизи антенного устройства. Располагая антенну рядом с поверхностью исследуемого объекта, изменяя ее размеры (глубину зондирования) и местоположение, удается получать информацию о его подповерхностной структуре. Частота зондирующего сигнала, как правило, выбирается из условия, чтобы характерный масштаб локализации квазистатического поля зонда был меньше глубины “скин-слоя” в исследуемой среде. Исходными данными для решения обратной задачи резонансной ближнепольной СВЧ-томографии являются результаты измерения резонансных характеристик диагностической системы при сканировании над поверхностью среды ее измерительной частью, представляющей собой систему электрически малых антенн с разными глубинами зондирования.
В данной статье анализируются особенности постановки обратной задачи ближнепольной СВЧ-томографии на примере резонансных датчиков с зондирующими элементами в виде цилиндрических конденсаторов с металлическим фланцем. Такие датчики успешно применялись нами для изучения изменений в структуре кожного покрова при различных патологиях [5]. томография диэлектрический проницаемость зонд
2. Резонансные измерительные системы
Измерительная система для исследования подповерхностных неоднородностей диэлектрической проницаемости представляет собой СВЧ-резонатор в виде отрезка коаксиальной линии длиной ?10см, с волновым сопротивлением =50ом, на одном конце которого находится нагрузка, включающая в себя зондирующую электрически малую антенну (измерительную ёмкость Сх), рис. 1. На противоположном конце расположена магнитная рамка диаметром 0,5 см (для этого центральный проводник коаксиальной линии замыкается на экран, образуя индуктивную нагрузку ). Возбуждение резонатора и прием его отклика осуществляется при помощи петель магнитной связи, расположенных вблизи магнитной рамки резонатора.
Рис. 1 Схема ближнепольной измерительной системы
Рис. 2 Внешний вид резонансной ближнепольной измерительной системы
Измерительная емкость Сх представляет собой краевую емкость цилиндрического конденсатора с фторопластовым заполнением. Внешняя обкладка конденсатора оканчивается плоским металлическим фланцем, рис.2 ( - емкость цилиндрического конденсатора без учета краевого эффекта). Длина конденсатора составляла 5мм, радиус внутреннего проводника - 3мм. Глубина зондирования системы определяется расстоянием между внешней и внутренней обкладками конденсатора, которое варьировалось в пределах от 0,25мм до 2мм. Для сильноконтрастных объектов такие зазоры обеспечивали глубины зондирования от 0.1мм до 1мм. Коррекция резонансной характеристики системы для конкретной исследуемой среды осуществляется при помощи индуктивности Lx, включённой параллельно с Сх. Торец цилиндрического конденсатора с плоским металлическим фланцем прикладывается к поверхности исследуемой среды. Наличие вещества в области краевого квазистатического поля цилиндрической емкости приводит к изменению резонансной частоты и добротности системы. Собственная резонансная частота измерительных систем составляла порядка 800 МГц, добротность ~ 150.
Согласно [3,4], резонансная частота и амплитуда сигнала в максимуме резонансной кривой датчика определяется величиной измерительной емкости в соответствии со следующими выражениями:
(1)
где
(2)
(3)
(4)
(5)
Таким образом, выражения (1)-(5) позволяют определить величину измерительной емкости при контакте ее с поверхностью неоднородной среды по экспериментально найденным и .
Значения для измерительных емкостей с разным масштабом локализации зондирующего поля (с разными зазорами между центральным и внешним проводниками цилиндрического конденсатора) являются исходными данными для решения обратной задачи по восстановлению пространственного распределения электродинамических характеристик среды.
3. Краевая емкость цилиндрического конденсатора
Рассмотрим цилиндрический конденсатор в виде двух проводящих коаксиальных цилиндра радиусами r1 и r2 (r2 > r1). Внешняя обкладка конденсатора оканчивается бесконечным металлическим фланцем, рис.3.
Рис. 3 Цилиндрический конденсатор с бесконечным металлическим фланцем. Распределение электрического поля и его потенциала без учета краевого эффекта
Пространство между обкладками заполнено средой с диэлектрической проницаемостью . Распределение электрического поля в такой системе без учёта краевых эффектов на торцах цилиндров хорошо известно. В плоскости поперечного сечения конденсатора радиальные линии представляют собой силовые линии электрического поля, концентрические окружности - линии равного потенциала. Электрическое поле вне конденсатора равно нулю.
Если - поверхностная плотность заряда на центральном цилиндре , разность потенциалов между обкладками и радиальное электрическое поле выражаются в следующем виде:
, (6)
(7)
Ёмкость на единицу длины конденсатора определяется выражением:
(8)
Учёт краевого эффекта приводит к изменению распределения поля в окрестности торца конденсатора. Во-первых, наличие электрического заряда на металлическом фланце и торце центрального проводника приводит к появлению электрического поля вне конденсатора. Во-вторых, внутри конденсатора наблюдается усиление поля вблизи торца с увеличением поверхностной плотности заряда на внешней и центральной обкладке и появление аксиальной компоненты электрического поля.
Для того, чтобы исключить взаимное влияние торцов конденсатора друг на друга, будем предполагать, что длина цилиндрического конденсатора неограниченна. Именно такая модель наиболее приближена к случаю, когда зазор между обкладками существенно меньше их длины.
Краевое электрическое поле цилиндрического конденсатора при контакте с поверхностью неоднородного полупространства с диэлектрической проницаемостью описывается уравнением Пуассона:
, (9)
c граничным условием в плоскости торца
(10)
Если комплексную диэлектрической проницаемость среды можно представить в виде малого отклонения , занимающего конечную область в полупространстве , от некоторого постоянного значения в остальной части этого пространства:
(11)
решение задачи (9)-(11) можно искать методом теории возмущений.
Потенциал , электрическое поле и поверхностную плотность заряда на торце центрального проводника представим в виде невозмущенных значений , и , соответствующих однородной среде и малых добавок , и , связанных с малой неоднородной компонентой комплексной диэлектрической проницаемости среды .
(12)
В нулевом приближении (среда однородная) электрическое поле вне конденсатора описывается уравнением:
(13)
с граничным условием (9). Решение этой задачи будем искать альтернативным способом. Краевые поля цилиндрического конденсатора с бесконечным металлическим фланцем могут быть вычислены как поля плоского кольцевого магнитного тока. Действительно, без учёта краевого эффекта на торце конденсатора, в плоскости z=0 имеется скачок радиального электрического поля.
(14)
В соответствии с граничным условием (14) данному разрывному полю соответствует поверхностный кольцевой магнитный ток ( - единичный орт вдоль азимутального направления), текущий в плоскости раскрыва между центральным и внешним проводниками цилиндрической емкости, рис. 4.
(15)
Стоит отметить, что таким путем, в частности, решаются задачи о возбуждении электромагнитного поля в волноводах щелями и другими отверстиями, прорезанными в их стенках.
Магнитный векторный потенциал , создаваемый током , в случае однородной среды вне конденсатора равен:
(16)
(штрихом отмечены координаты источника).
Величина , как легко установить из рис. 4, равна: . Электрическое поле выражается через магнитный векторный потенциал следующим образом:
.
При вычислении вектора рассмотрим элементарную площадку на плоском кольце. Если сначала проводить интегрирование по , а затем по , то плоское кольцо разбивается на систему магнитных кольцевых токов шириной каждый. В силу осевой симметрии задачи (все величины не зависят от азимутального угла) можно без потери общности выбрать точку наблюдения в плоскости (при этом ). При вычислении вклада в интеграл введен множитель , а интегрирование по проводится от 0 до . В результате, магнитный векторный потенциал принимает вид:
(17)
Рис. 4 Расположение источника и точки наблюдения для кольцевого магнитного тока
Составляющие поля и находятся путём дифференцирования , и выражаются следующими формулами:
(18)
(19)
Выполнение процедуры дифференцирования позволяет получить выражения для и в виде двукратных интегралов. Проведя интегрирование по переменной можно свести интегралы к однократным:
(20)
(21)
где и полные эллиптические интегралы Лежандра первого и второго рода соответственно.
(22)
и представляются в виде степенного ряда:
, (23)
В соответствии с этим, полные эллиптические интегралы можно с заданной точностью аппроксимировать полиномиальной зависимостью. Это обстоятельство позволяет сократить вычислительное время при расчете полей с помощью (20) и (21) по сравнению с численным интегрированием двукратных интегралов, соответствующих (18)-(19).
На рис. 5-8 представлен численный расчёт и компонент электрического поля вне цилиндрического конденсатора с бесконечным металлическим фланцем. Радиус внутреннего проводника составлял r1=3мм, радиус внешнего проводника - r2=3.5мм. Видно, что в плоскости раскрыва радиальная компонента поля отлична от нуля только в зазоре между внешним и внутренним проводником конденсатора. По величине совпадает с полем в конденсаторе без учёта краевого эффекта. компонента монотонно растет по координате r, достигает максимума при r=r1, затем уменьшается, меняет знак и достигает минимума при r=r2. При r=r* (r1<r*<r2) обращается в нуль. По мере удаления от плоскости раскрыва, электрическое поле спадает с характерным масштабом порядка зазора между внутренней и внешней обкладками конденсатора.
Рис. 5 Радиальное распределение (а) и (б) компонент электрического поля в плоскости раскрыва цилиндрического конденсатора с бесконечным металлическим фланцем. (r1=3мм, r2=3,5мм)
Рис. 6 Пространственная структура компоненты краевого электрического поля цилиндрического конденсатора с бесконечным металлическим фланцем (r1=3мм, r2=3,5мм)
Рис. 7 Пространственная структура компоненты краевого электрического поля цилиндрического конденсатора с бесконечным металлическим фланцем (r1=3мм, r2=3,5мм)
Рис. 8 Пространственная структура модуля краевого электрического поля цилиндрического конденсатора с бесконечным металлическим фланцем (r1=3мм, r2=3,5мм)
Поверхностная плотность заряда на торце центрального проводника определяется в соответствии с граничным условием на металле:
(24)
где - аксиальная компонента вектора индукции электрического поля. Если е - диэлектрическая проницаемость среды вне конденсатора (считается, что среда однородная), то .