Расчет параметров локальных первичных преобразователей для измерения вариаций электропроводности Земли Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

НАУЧНОЕ ИЗДАНИЕ МЕТУ ИМ. Н. Э. БАУМАНА

НАУКА и ОБРАЗОВАНИЕ

Эл № ФС77 - 48211. Государственная регистрация №0421200025. КБМ 1994-0408

электронный научно-технический журнал

Расчет параметров локальных первичных преобразователей для измерения вариаций электропроводности Земли

# 09, сентябрь 2012

Б01: 10.7463/0912.0476222 Мисеюк О. И.

УДК 550.380.001

Россия, МГТУ им. Н.Э. Баумана omiseyuk@mail.ru

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

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

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

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

наблюдении за вариациями электропроводности. В работе [3] приведено обоснование метода измерения вариаций электрической проводимости и приведена расчетная модель формирования электрического поля излучателя при наличии вертикального контакта, имитирующего деформацию среды, приводящую к скачку электропроводности. Рассматривается ситуацию, когда первоначальная структура объекта уже неоднородна. Тогда при заданных и фиксированных параметрах внешнего источника поле определяется удельным сопротивлением однородных областей и зависит от свойств и расположения этих областей. Пусть Е1- поле неоднородного объекта до появления новой неоднородности, а Е2- поле объекта после появления новой неоднородности. Для обнаружения вариаций электрического состояния среды представляет интерес разностное поле А Е= Е 2 - Е1. Разностное поле А Е можно представить как поле некоторого источника, обусловленного и определяемого возникшей неоднородностью. Именно это поле необходимо измерить. При этом исходное поле не должно восприниматься прибором.

Таким образом, измерительная установка должна удовлетворять следующим требованиям:

— максимальная чувствительность к разностной составляющей поля;

— подавление сигнала от исходного поля, которое рассматривается как помеха;

— необходимость оптимизации по отношению А Е/Е 1, т.е. по отношению сигнал/помеха. Расчетная модель исследуемого поля (рис.1) , предложенная в [3], содержала два задающих поле

электрода А и В, через которые протекает ток I от источника V, расположенные вблизи границы раздела земля-воздух на расстоянии 1АВ = 21 друг от друга. На расстояниях г ы = п и гМ = т от

каждого из них расположены точки наблюдения N и М возле электрода А и точки N1, М1 около электрода В. На расстоянии d от центра диполя расположен плоский вертикальных контакт, имитирующий деформацию среды. В точках наблюдения для однородного полупространства при наличии плоского вертикального контакта были получены выражения потенциала и напряженности электрического поля [3]. Измерение величин напряженности электрического поля в точках наблюдения осуществляется с помощью локальных первичных преобразователей (ЛИП).

21 М< < [В М N ... А А < 21 А М <4

т п \ ^ 0 |д] £ , Л7 е! . ¥

1 X - 1.

<зг

к

21

е

л/

21 А и

в!

ЛКп 21Кп.

£

Рис. 1. Расчетная модель исследуемого поля при наличии вертикального контакта

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

Локальный первичный преобразователь

Пусть исходное состояние среды характеризуется наличием в области измерений неоднородности в виде вертикального контакта (рис.1), разделяющего среду I с проводимостью сх и среду II, где проводимость равна с2 = а1 ±АО. Индексом «О» обозначены значения величин в исходном состоянии среды.

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

Напряженность исследуемого поля в точках М и N при наличии плоского контакта равна [3]

Е = у 4пох

{(у -1)Р + (у -1)2]-3/2 - (у +1)Р + (у +1)Г2 -

-3/2

-K12{(2d - у +1)[р2 + (2^ - у +1)2]-3/2 - (2й - у -1)[р2 + (2^ - у -1)2]-3/2}} (1)

2л -3/2'

о = \ х: - г:, коэффициент отражения

К

12

_(о,1-<т2)_

удельная электрическая проводимость

Земли;

(Т2 =(7! +Л - удельная проводимость среды после появления новой неоднородности.

Обозначив относительное изменение электропроводности —^ = 5, и с учетом того, что

2с1

5

5 <<1, получим К 12 =-« -5 . Распределение напряженности электрического поля

1 + 5

при наличии вертикального контакта представлено на рис.2.

Рис. 2 Распределение напряженности электрического поля двух точечных источников

при наличии вертикального контакта

21

4пст1

(—-------------1

т2 (21 + т)

2 /0 7 , л2_) - ^12

1

1

(2$ — т) (2$ — 21 — т)

(2)

В точке наблюдения N

Е = -

J^YN

21

4пст1

■ + ■

1

и

(2/ — п)2

12

1

+

1

(2$ + п) (2$ — 2/ + п)

(3)

Тогда, обозначив коэффициент преобразования локального ПП- 1Э ( эквивалентная база); иХ-напряжение холостого хода на выходе первичного преобразователя; ЕМ, Ем -напряженность исследуемого поля, получим выражения для сигналов на выходе ЛПП, расположенных в точках М и К, при наличии плоского контакта при отсутствии возмущающего воздействия

2

1

2

Vмо — Ем’ *эМ=2ЫзМ([_---------------------------

ІТТ(Ті 11т2 (21+т)2

- К

12

— Еы • 1ъм —

N LэN

гг

4^(7]^ (. I

4 п(г1 11_п2 (21—п)2.

- К

(2 с1-т)2

1

12

(2(1-21-т)2 1

(2<1+п)2 (2гі-2!+п)2

]}

і

(4)

(5)

Вследствие различных знаков напряжений ИМо и ИКо условие компенсации исходного поля в этом случае запишется ихм?1 = = 0 .

С учетом поля помехи сигнал ЛПП в точке М будет равен

а в точке N и'а=и?}0 + 17ш = (—Ен + Еп) - 1ЭК.

Напряжение холостого хода на выходе двух ЛПП, расположенных в точках М и К, при выполнении условия компенсации составит их = ^'М0+и^0 = Еп- 1эм+^пт ^эя ■

С целью компенсации сигнала помехи на выходе измерителя при отсутствии возмущающего воздействия следует дополнительно к основным использовать еще два ЛПП и разместить их вблизи излучающего электрода В в точках М1 и N1, расположенных на расстоянии гм1=ш и гК1=п соответственно от электрода В на линии, соединяющей электроды А и В. Размещение измерительных электродов в точках М, N и М1,К1, где

напряженности равны по величине, но противоположны по знаку, компенсирует поле излучающих электродов на выходе каждой пары ЛПП, но удваивает сигнал помехи при условии 1ЭМ=1ЭМ. Для компенсации сигнала помехи сигналы с основного и дополнительного датчиков должны вычитаться при дальнейшей обработке сигнала. При появлении возмущающего скачка проводимости, возникающего на расстоянии ё от центра излучающего физического диполя на границе вертикального контакта, в точках наблюдения изменяются величины напряженности ЕМ и Е№ Как следует из рис 2 величина ЕМ увеличивается, а Ем уменьшается, что приводит к появлению сигнала на выходе датчика, образованного ЛППМ и ЛППМ. Предположив 1эм=1эм=1э, получим

и

XMN

. (6)

При этом на выходе датчика, образованного ЛПП, расположенными в точках М1 и N1, напряжение холостого хода будет равно

и

XM1N1

Для определения соотношений между т,п и 1э, которые обеспечили бы выполнение условий компенсации исходного поля излучателей, были рассчитаны зависимости = ДДШ11) при фиксированных значениях параметров с1 и Д0. Причем в случае использования локальных ПП соотношение Ме и N определяется формулой

Ме= Ne - Д.

(8)

л . тп п гча тт

где Ме=у; Ne—у; ^~у Для компенсации исходного поля при измерении напряженности с

помощью локальных ПП расстояние Меот датчика, расположенного в точке М, до излучающего электрода А должно быть меньше расстояния Ne, отделяющего датчик в точке N от того же электрода, в отличие от измерителя на основе ИПП, где с той же целью должно выполняться соотношение Mp = Ny + Атп. В таблице 1 приведены значения Ne для фиксированных

Дтт1, D иД0 и при отсутствии возмущающего воздействия ДБ0,М= 0. Из зависимости Ыс(Дт„.) видно, что при Ne< 0,2 Дтп < 0,005 т.е. равенство Ме= Ne выполняется с погрешностью менее 5 % и не зависит от значений D иД0.При Ne>0,2 Атп возрастает, и при выборе места расположения ЛПП следует пользоваться соотношением (8) и данными таблицы 1.

Значения N6 для локального ПП

А 0 Б \ А тп 0,01 0,02 0,04

0,0 1,5 0,3370986 0,4214993 0,5253700

3,0 0,3370988 0,4214990 0,5253707

0,01 1,5 0,3395692 0,4249140 0,5302380

3,0 0,3371667 0,4215783 0,5254575

0,1 1,5 0,3659744 0,4645886 0,5262067

3,0 0,3377527 0,4222564 0,5062067

При появлении возмущающего воздействия ДВ0:5М^ 0 на выходе ЛПП в точках М и N и Мь N1 появляются сигналы, пропорциональные соответственно

г _ ^хМ! N!

и 21Ц/4ШГ!'

Сравнительный анализ данных расчета Ае(^тп) и ^зеС^тл) Для О = 1,5 и для Б= 3 при Д0= 0,0; 0,01; 0,1 позволяет сделать вывод, что изменение Д0 в указанных пределах слабо влияет на величину /1е. Удаление вертикального контакта от центра физического диполя приводит к уменьшению сигнала на выходе датчика. Сигнал с дополнительных ЛПП, расположенных в точках М] и N1, пропорциональный величине /2е, на порядок меньше сигнала с основных ЛПП в точках М и N при прочих равных условиях. Величина /1е монотонно возрастает с увеличением Дтп и ДВ03М, /2е также возрастает с увеличением Дтп и ДВ03М при До.=0,0 и 0=1,5 ...3, а при Д0.= 0,01 и 0=1,5 /зе(Дтп) уменьшается и меняет знак.

Сопоставительный анализ датчиков на основе интегрального и локального первичных преобразователей

Функциональная блок-схема датчиков на основе ЛПП представлена на рис. 3. Датчик на базе ЛПП содержит четыре ЛПП, размещенных в точках М и N (ЛПП1 и ЛИ! 12) и Мь^ (ЛПП3 и

ЛПП4). С целью компенсации поля излучающих электродов при отсутствии возмущающего воздействия первичные преобразователи должны быть установлены в точках, где напряженность исходного поля имеет одинаковые значения. Как видно из рис. 2 в точках М и N как и в М1 и напряженность поля имеет противоположные по знаку и равные по величине значения. Это позволяет, подав сигнал с ЛПП1 и ЛПП2 через усилители У1 и У2 на сумматор С1, скомпенсировать поле излучающих электродов при отсутствии возмущающего воздействия. Однако сигнал однородной помехи при суммировании им и ^ удваивается. Для компенсации сигнала помехи используются два дополнительных ПП - ЛПП3 и ЛПП4, расположенные вблизи излучающего электрода В идентично ЛПП1 и ЛПП2. После суммирования их сигналов на сумматоре С2 на выходе этого сумматора будет сформирован удвоенный сигнал помехи, равный по величине и одинакового знака с сигналом на выходе сумматора С1. Для компенсации сигналов помехи в устройстве используется вычитатель В. На выходе датчика при этом формируется сигнал, пропорциональный разности ([1в — [2ё)- При отсутствии возмущающего воздействия этот сигнал равен 0. При появлении скачка проводимости, вызванного деформацией среды, изменяется картина поля в месте расположения измерительных электродов и на выходе датчика появляется полезный сигнал, пропорциональный возмущающему воздействию.

Рис.3. Функциональная блок-схема датчика на основе ЛПП

С целью проведения сравнительного анализа датчиков на основе ИПП и ЛПП были построены зависимости

и

Ґе Хв/2І13І4па1 &■**-&* £?((Ао + ЛвоЗМ)

при фиксированных Э, Д0 и Атп. Графики этих зависимостей приведены на рис. 4, 5, где Э принималось равным 1,5 и 3; Д0 = 0,0; 0,01; ОД. На каждом графике приведено семейство/^, (Д0 + ДВ03М) либо /е({Д0 + ДВозм).- параметром семейства является величина Атп,

однозначно определяющая расположение электродов в точках М и N при котором выполняется условие компенсации поля излучателей, когда возмущающее воздействие отсутствует. Дтп задавалось равным 0,01; 0,02 и 0,04. Приведенные зависимости характеризуют изменение коэффициента преобразования измерителей при возрастании возмущающего воздействия ДБ03М. На графиках наблюдается линейное возрастание сигнала на выходе датчика при увеличении ДБ03М. В таблице 2 приведены значения коэффициентов преобразования измерителей на основе ИПП и ЛПП соответственно К<р=[<р/(А0 + ДБозм) и Ке = Те/№а + ^возм)? их относительного изменения

0,02; 0,04. Расчеты проведены при Двозм=0,1 и при Д(ДВ03М)= 0,04. Приведенные численные данные подтверждают преимущество использования ЛПП в измерителях изменений электропроводности среды. Наряду с отмеченными в [2] достоинствами ЛПП, такими как меньший, чем у ИПП уровень собственных шумов и большая предельная чувствительность, в ЛПП появляется возможность компенсации исходного поля излучателей не только за счет геометрии установки, как у интегральных ПП, но и за счет возможности регулировки коэффициентов усиления в каналах ЛПП. Этим же путем устраняется не идентичность каналов измерения. Следует отметить превышение коэффициента преобразования Ке над К^ более чем в 4

раза при расположении вертикального контакта на расстоянии сІ=1,5і, и в 3,6 2,3 раза при сІ=Зі.

то

0.05

аю

0/5

020

-ом

-0№

-оде-

-Ц06

Рис. 4. Зависимость /р(ДЕ03б) ПРИ ¿о=0,01;

Рис. 5. Зависимость/*, (Двозб) при До=0,01; Б=1,5

На основании полученных данных и с учетом данных, приведенных в [2],можно оценить порог разрешения измерительной системы по относительному изменению удельной электрической проводимости при использовании локальных ПП

Если принять 10= ОД... 0,2 А, = 10 2См/м; и^= 8- 10 15В2[4],то Д

л/2-810-15-2-ЗД410

1-2

= 0,585 ■ Ю“4.

Лвозмтт

1,160,2 0.2340,25

Таким образом, измеритель способен реагировать на скачки удельной электропроводности порядка 10-4%.

Коэффициенты преобразования измерителей на базе И1111 (К^) и J11111 (Ке)

Кр 0,235 0, 3 1 5 0,445

0,00 Щр 0,2 0,3125 0,4375

ке 1,16 1,393 1,84

ьке 1,03 1,25 1,65

Ке/Кф 4,936 4,42 4,13

1,5 Ку 0,23 0,31 0,445

0,01 Щр 0,22 0,3 0,43

Ке 1,125 1,38 1,83

д ке 1,04 1,36 1,76

Ке/Кф 4,89 4,45 4,11

Ку 0,238 0,334 0,49

0,1 Щр 0,188 0,281 0,434

ке 1,17 1,5 1,82

ьке 0,95 1,3625 1,77

Ке/Кф 4,93 4,49 4,10

К<р 0,006 0,0076 0,0097

0,00 Щр 0,005 0,007 0,00875

Ке 0,02182 0,02227 0,02318

д ке 0,02045 0,02103 0,02046

Ке/Кф 3,637 2,93 2,3896

3,0 Ку 0,0006 0,0076 0,0098

0,01 Щр 0,0054 0,0068 0,0084

Ке 0.022 0,0225 0,023

лке 0,021 0,021 0,022

Ке/% 3,67 2,96 2,347

К<р 0,0056 0,00726 0,00917

0,1 Щр 0,00446 0,00565 0,00804

Ке 0,02125 0,0215 0,022

лке 0,01875 0,01875 0,018

Ке/Кр 3,795 2,96 2,4

Выводы

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

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

Показано, что с целью компенсации поля однородной помехи следует использовать дополнительные ЛПП, расположенные около второго излучающего электрода. Основные и дополнительные ЛПП следует располагать на прямой, проходящей через излучающие электроды. Если измерительные ПП расположены на расстоянии Гм (м) < 0,1/ от

излучающих электродов, то условие компенсации исходного поля выполняется при гм = гк, т.е. измерительные ПП расположены на одинаковом расстоянии от излучающих электродов. При увеличении этого расстояния т.е. при Гм(м) > 0,2/ Гм должно быть меньше на Атп для того,

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

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

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

Получены зависимости выходных сигналов основных и дополнительных ЛПП при различных возмущающих воздействиях и различной геометрии расположения ЛПП.

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

Список литературы

1. Зимин Е.Ф., Мисеюк О.И., И.И. Плаксин И.И. Способ измерения временных вариаций кажущейся удельной электрической проводимости Земли и устройства для его осуществления: пат. 94025990/25 Российская Федерация.1996. Бюл. № 18. 4 с.

2. Зимин Е.Ф., Кочанов Э.С. Измерение параметров электрических и магнитных полей в проводящих средах. М.: Энергоатомиздат, 1985. 255 с.

3. Мисеюк О.И.О способе измерения временных вариаций удельной электрической проводимости Земли // Наука и образование. МГТУ им. Н.Э.Баумана. Электрон. журн. 2012. № 3. Режим доступа: http://technomag^и.ru/doc/338840.html (дата обращения 18.03.2012).

4. Гохберг М.Б., Гуфельд И.Л., Гершензон Н.И., Пилипенко В.А. Электромагнитные эффекты при разрушении земной коры // Известия АНСССР. Физика Земли. 1985. № 1. С. 72 - 87.

SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU

SCIENCE and EDUCATION

EL № FS77 - 48211. №0421200025. ISSN 1994-0408

electronic scientific and technical journal

Calculation of parameters of local transducers for measuring electrical conductivity variations of Earth

# 09, September 2012 DOI: 10.7463/0912.0476222 Miseyk O.I.

Russia, Bauman Moscow State Technical University

omiseyuk@mail.ru

The author considers problems of designing local transducers of specific electric conduction which are mounted in industrial impact zones. The author obtained the ratio for determining configuration of electrodes in local transducers which provide signal interference compensation. The author proposes a flowchart for transducers measuring electrical conductivity of the medium, on the basis of local transducers which have compensation of signal interference fields and fields of unperturbed medium. A comparative analysis of sensors based on local and integral transducers is given in the article.

Publications with keywords:potential, intensity of electric field, primary converter, the specific electrical conductivity, the disturbance, the measuring electrode, electromagnetic field, the boundary-fence of the media, signal interference

Publications with words:potential, intensity of electric field, primary converter, the specific electrical conductivity, the disturbance, the measuring electrode, electromagnetic field, the boundary-fence of the media, signal interference

References

1. Zimin E.F., Miseiuk O.I., Plaksin I.I. Sposob izmereniia vremennykh variatsii kazhushcheisia udel'noi elektricheskoi provodimosti Zemli i ustroistva dlia ego osushchestvleniia [Method for measuring the temporal variations of the apparent electrical conductivity of the Earth and devices for its implementation]. Patent RF, no. 94025990/25/, 1996.

2. Zimin E.F., Kochanov E.S. Izmerenieparametrov elektricheskikh i magnitnykhpolei v provodiashchikh sredakh [Measurement of electric and magnetic fields parameters in conducting media]. Moscow, Energoatomizdat, 1985. 255 p.

3. Miseiuk O.I.O sposobe izmereniia vremennykh variatsii udel'noi elektricheskoi provodimosti Zemli [About the way of measurement of temporal variations of Earth’s electric conductivity]. Nauka i obrazovanieMGTUim. N.E. Baumana [Science and Education of the Bauman MSTU], 2012, no. 3. Available at: http://technomag.edu.ru/doc/338840.html , accessed 18.03.2012.

4. Gokhberg M.B., Gufel'd I.L., Gershenzon N.I., Pilipenko V.A. Elektromagnitnye effekty pri razrushenii zemnoi kory [Electromagnetic effects in the destruction of the Earth's crust]. Izvestiia ANSSSR. Fizika Zemli [Journal of the USSR Academy of Sciences. Physics of the Earth], 1985, no. 1, pp. 72-87.