RU2818011C1 - Способ подавления сигнала прямого поля в электромагнитном индукционном зондировании и устройство для его реализации - Google Patents
Способ подавления сигнала прямого поля в электромагнитном индукционном зондировании и устройство для его реализации Download PDFInfo
- Publication number
- RU2818011C1 RU2818011C1 RU2023121274A RU2023121274A RU2818011C1 RU 2818011 C1 RU2818011 C1 RU 2818011C1 RU 2023121274 A RU2023121274 A RU 2023121274A RU 2023121274 A RU2023121274 A RU 2023121274A RU 2818011 C1 RU2818011 C1 RU 2818011C1
- Authority
- RU
- Russia
- Prior art keywords
- signal
- direct field
- generator
- coil
- amplitude
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000005674 electromagnetic induction Effects 0.000 title claims abstract description 5
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 16
- 238000001914 filtration Methods 0.000 claims abstract description 16
- 230000003044 adaptive effect Effects 0.000 claims abstract description 14
- 230000001629 suppression Effects 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000013461 design Methods 0.000 abstract description 7
- 238000005259 measurement Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 abstract 1
- 230000005611 electricity Effects 0.000 abstract 1
- 239000000126 substance Substances 0.000 abstract 1
- 230000006870 function Effects 0.000 description 20
- 230000006698 induction Effects 0.000 description 11
- 230000014509 gene expression Effects 0.000 description 6
- 239000000523 sample Substances 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000013178 mathematical model Methods 0.000 description 4
- 230000005672 electromagnetic field Effects 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 244000309464 bull Species 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010230 functional analysis Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Abstract
Изобретение относится к электромагнитным индукционным зондированиям, в частности к подавлению сигнала прямого поля в приемном сигнале. Сущность: подавление сигнала прямого поля выполняют в блоке фильтрации с использованием многоканального адаптивного RLS-алгоритма путем определения весовых коэффициентов при минимизации невязки целевой функции, характеризующей разность между измеренным и моделируемым сигналами с приемной катушки, формированием по полученным значениям коэффициентов симплекса в соответствии с методом Нелдера-Мида для нахождения безусловного минимума целевой функции. Проводят процентную оценку наличия сигнала прямого поля в регистрируемых данных и вычисляют амплитуды и фазы полезного сигнала на основе весовых коэффициентов. Устройство состоит из генераторной катушки, генератора гармонических сигналов, датчика тока, установленного внутри генераторной катушки и дающего информацию о прямом поле, а также одной или набора приемных катушек. Сигналы от катушек через соответствующие цепи, имеющие полосовой фильтр и аналого-цифровой преобразователь, соединены с блоком фильтрации, где амплитуда и фаза полезного сигнала вычисляются с использованием многоканального адаптивного RLS-алгоритма. Технический результат: повышение надежности выделения полезного сигнала и упрощение конструкции геофизической аппаратуры электромагнитного зондирования. 2 н.п. ф-лы, 3 ил., 2 табл.
Description
Область техники
Изобретение относится к геофизической аппаратуре электромагнитных индукционных зондирований, в частности, к способам подавления сигнала прямого поля в приемном сигнале приборов.
Уровень техники
В геофизических устройствах для электромагнитного индукционного зондирования локальные неоднородности удельного сопротивления породы выявляются при создании электромагнитного поля генераторной катушкой и измерении сигнала от токов, индуцированных в земле [1, 2]. Переменный ток в генераторной катушке индуцирует э.д.с. прямого поля в приемной катушке и наводит ток в породе или проводящих объектах, находящихся вблизи устройства. Переменное магнитное поле от индуцированного тока называется вторичным полем, и оно создает в приемной катушке полезный сигнал, величина которого обратно пропорциональна удельному сопротивлению окружающей горной породы. Соответственно, на приемную катушку действует переменное магнитное поле нескольких источников. Сложность выделения полезного сигнала с приемной катушки заключается в превосходстве величины сигнала от прямого поля, на фоне которого полезный сигнал от вторичного поля сопоставим с помехой. Схемы подавления первичного поля имеют различное строение, настройка и калибровка, которых требует применения сложных электротехнических и конструкционных решений, например, таких как, в патенте RU2461850 [3], где применяют замкнутое проволочное кольцо, над которым выполняют измерение сигнала от индуцированных токов в кольце и, сравнивая с расчетом, подбирают моменты катушек и расстояния от центра генераторного диполя до центров приемных диполей.
В многочастотных приборах, патент RU2152058 [1], переменное магнитное поле создается последовательно на многих частотах. На каждой частоте выполняют измерение вторичного магнитного поля двумя приемными датчиками, расположенными вдоль одной прямой с генераторной петлей. Сигнал от прямого поля предварительно подавляют во всем частотном диапазоне (от 2,5 кГц до 250 кГц) методом вычитания сигналов двух катушек, соблюдая условие равенства отношения момента первого датчика к кубу расстояния от центра генераторной петли до центра первого датчика и отношения этих же величин для второго датчика. Недостатком данного метода подавления сигнала прямого поля является уменьшение полезного сигнала. Расчеты показывают, что при идеальной фазировке реальной компоненты, разностный сигнал трехкатушечного зонда над грунтом в 1.7 раза меньше, чем теоретическая реальная компонента полезного сигнала двухкатушечного зонда [4].
Устройство, описанное в патенте RU2502092 [2], предназначено для измерения электропроводности почвы и содержит плоский излучающий контур - источник переменного тока, непроводящую штангу, расположенную под углом к плоскости излучающего контура, набор приемных катушек, закрепленных на непроводящей штанге. Близкий к нулю вертикальный магнитный поток прямого поля через приемные катушки достигается с помощью наклона непроводящей штанги под углом близким к 35 градусам. На этой образующей вертикальная компонента магнитной индукции отсутствует. Недостатком аппаратуры является изменение первоначальной настройки и уровня подавление прямого поля при транспортировке и эксплуатации, ввиду слабой устойчивости к механическим повреждениям [5].
Наиболее близким аналогом к изобретению является устройство со способом обнаружения присутствия слабого магнитного поля, описанным в патенте US9389281 [6], заключающимся в адаптивном формировании компенсационного сигнала от первичного магнитного поля. При такой схеме подавления, недостатком аппаратуры является многократное цифро-аналоговое преобразование (ЦАП) данных от магнитного датчика, в ходе которого может происходить «смещение нуля» согласно описанию характеристик ЦАП, а также согласно теории адаптивной фильтрации, в установившемся режиме значения весовых коэффициентов адаптивных фильтров могут находиться в небольшой окрестности истинного минимума [7], что в совокупности и по отдельности влияет на точность выделения полезного сигнала.
Раскрытие сущности изобретения
Задачей изобретения является повышение надежности измеряемых данных и упрощение конструкции геофизической аппаратуры.
Повышение надежности и упрощение конструкции достигается за счет того, что подавление сигнала прямого поля в индукционной аппаратуре, содержащей генераторную катушку, генератор гармонических сигналов, датчик тока, установленный внутри генераторной катушки и дающий информацию о прямом поле, а также одну или набор приемных катушек (в случае геометрического зондирования - изменение расстояния между генератором и приемником), сигналы от которых через соответствующие цепи, имеющие АЦП (аналого-цифровой преобразователь), преобразовываются в цифровой вид, отличается тем, что подавление сигнала прямого поля выполняется в блоке фильтрации в процессе выделения основных параметров полезного сигнала (сигнал от среды) путем нахождения весовых коэффициентов математической модели, соответствующей сигналу с приемной катушки, что позволяет, например, использовать одну приемную катушку, если сравнивать с патентом RU2152058 [1], причем генератор и приемник жестко закреплены горизонтально в одной плоскости, поэтому настройки аппаратуры не будут изменяться в ходе транспортировки и эксплуатации как в патенте RU2502092 [2].
Ввиду постоянного изменения регистрируемых данных и их зависимости от частоты, целесообразно в аппаратуре индукционного зондирования использовать адаптивные системы фильтрации [7], позволяющие подстраивать изменяемые параметры (весовые коэффициенты) передаточной функции в процессе работы алгоритма. Так как критерием работы адаптивного фильтра часто является достижение минимума целевой функции (функция входного, выходного и требуемого сигналов [7]), определяющего набор параметров фильтра [7], то для повышения надежности и упрощения конструкции критерий работы фильтра должен учитывать преобладание помехи над полезным сигналом и его смещение по фазе, обусловленное свойством исследуемого объекта [8].
В реализованном алгоритме целевая функция содержит математическую модель сигнала с приемной катушки с неизвестными весовыми коэффициентами, при изменении которых невязка между измеренным сигналом с приемной катушки и теоретическим сигналом математического моделирования в целевой функции стремилась к минимуму, - такой подход широко известен и применяется для решения обратных задач [9].
Поскольку при обработке используются стационарные сигналы, то, согласно [7], наиболее подходящей целевой функцией является наименьшие квадраты ошибок:
где Uнаблюденное - сигнал с приемной катушки, Uтеоретическое - математическая модель сигнала, которая представляет реальный процесс взаимодействия электромагнитных полей в генераторной петле, датчике тока и приемной катушке.
Согласно теории, индукция электромагнитного поля внутри генераторной петли равна [10, 11]:
где μ0 - магнитная постоянная, Гн/м; Sg - площадь генераторной петли, м2; Ng - количество витков в генераторной петле; I0 - амплитуда тока, протекающего в цепи, A; Rg - радиус генераторной катушки, м; ƒ - частота, Гц; t - время, сек.
Э.д.с, в датчике тока, который изготовлен в виде небольшой катушки, установленной внутри генераторной петли, равен производной от потока по времени со знаком минус и определяется следующим выражением [11]:
где St - площадь катушки датчика тока; Nt - количество витков в датчике тока.
Сигнал в приемной катушке равен сумме полезного сигнала и сигнала помехи. Э.д.с. прямого поля в приемной катушке отличается по амплитуде и находится в противофазе относительно э.д.с. от датчика тока, по причине их расположения друг относительно друга. Датчик тока находится внутри генераторной петли, а приемная катушка находится за ее пределами, причем обе катушки находятся в одной плоскости. Полезный сигнал, возникающий от индуцированных токов в среде за пределами генераторной петли, отстает на π/2+ϕ (фаза) от первичного поля, а также имеет свою амплитуду [8].
Математическое выражение модели сигнала с приемной катушки можно представить в виде:
где А - амплитуда сигнала прямого поля в приемной катушке; С - амплитуда полезного сигнала.
Последнее выражение можно преобразовать, применяя формулу косинуса разности двух углов и приведя подобные слагаемые. Таким образом математическая модель примет вид:
где:
Для оценки величины э.д.с. прямого поля в приемной катушке применяются данные с токовой катушки. Сигнал с датчика тока описывается выражением (3), с помощью которого, зная амплитуду образцового сигнала, можно вычислить амплитуду тока в генераторной петле:
где εt - амплитудное значение э.д.с. с датчика тока.
При достаточном удалении генераторной петли от центра приемной катушки [10], равном больше трех диаметров генератора, последний можно считать диполем и тогда амплитуда сигнала прямого поля в приемной катушке будет равна:
L - расстояние от центра генераторной петли до центра приемной катушки; Sp - площадь генераторной петли; Np - количество витков в генераторной петле;.
Также стоит отметить, что постоянный контроль за током в генераторе необходим для корректной работы аппаратуры, поскольку любое электронное устройство так или иначе подвергается внешним воздействиям, которые влияют на общий функционал. Таким образом конечный вид целевой функции примет вид:
где k - это индекс элемента дискретной записи сигнала.
После нахождения неизвестных весовых коэффициентов (A1, C1) целевой функции, амплитуду и фазу полезного сигнала рассчитываются по следующим формулам:
Отличительной особенностью предлагаемого способа подавления прямого поля, является выделение параметров полезного сигнала, в частности фазы сигнала, которая определяется не через отношение мнимой компоненты сигнала к реальной, как это принято в большинстве приборов индукционного зондирования.
Нахождение истинной амплитуды и фазы полезного сигнала возможно при условии, что целевая функция будет иметь единственный минимум в области допустимых значений. Достаточность доказывается вычислением производных Фреше [12]:
Для удобства вычисления знак суммы в целевой функции был заменен интегралом с границами интегрирования от -∞ до +∞. После вычислений был получен следующий результат:
где ω=2⋅π⋅ƒ - циклическая частота.
Дальнейшее преобразование для вычисления производных Фреше не обязательно, поскольку на данном этапе вычислений видно, что выражения похожи на уравнение прямой, имеющих единственную точку минимума в области допустимых значений для измеряемых данных.
Ввиду того, что в реальных условиях данные с приемной катушки могут содержать собственные помехи аппаратуры или внешние помехи, которые не будут учитываться датчиком тока, после нахождения безусловного минимума функции (9), перед расчетом характеристик полезного сигнала по формулам (10), выполняется проверка наличия полезного сигнала с приемной катушки, основанная на том, что коэффициент А1 определяется характеристиками полезного сигнала и прямого поля согласно выражению (6). Поэтому нахождение относительной погрешности между рассчитанным амплитудным значением сигнала прямого поля по формуле (8) и найденным коэффициентом A1, позволяет оценить наличие полезного сигнала в регистрируемых данных, а если вычесть это значение из 100%, то будет получена процентная оценка сигнала прямого поля, которая позволяет после нахождения коэффициентов А1 и С1 принять решение о вычислении амплитуды и фазы полезного сигнала в соответствии с формулами (10) или пропустить этот этап, полагая что полезный сигнал в измеренных данных отсутствует.
Техническим результатом, при выполнении совокупности приведенных признаков, является повышение надежности выделения полезного сигнала и упрощение конструкции геофизической аппаратуры электромагнитного зондирования, за счет применения цифровой фильтрации, обеспечивающей глубокое подавление сигнала прямого поля на всех частотах и позволяющей расширить динамический диапазон аппаратуры, путем нахождения минимального значения невязки в целевой функции между теоретическим сигналом математического моделирования и измеренным сигналом с приемной катушки, для определения параметров полезного сигнала, в совокупности дающих информацию о среде и объектах подземного пространства, а также простота конструкции аппаратуры по предлагаемому способу уменьшает влияние внешних механических воздействий, не сбивающие настроек и калибровок аппаратуры.
Краткое описание графических материалов
На фиг. 1 представлена структурная схема предлагаемого устройства -аппаратуры индукционного зондирования подземного пространства с блоком подавления прямого поля (фильтрации).
Предлагаемое устройство аппаратуры индукционного зондирования включает в себя генераторную катушку 1, которая подключена к электроразведочному генератору гармонических сигналов 2, который в свою очередь управляется блоком управления и регистрации 3, питающимся от общего блока питания аппаратуры 4. Приемная катушка 5 преобразовывает переменное магнитное поле в э.д.с, соединена с полосовым фильтром 6 и с АЦП 7. Цифровой входной сигнал от приемной катушки 5, представляет собой аддитивную смесь полезного сигнала и сигналов помех, наибольший из которых происходит от первичного поля генератора. Образцовый сигнал получается от комплекса - датчик тока 8, полосовой фильтр 6, АЦП 7, поступает на вход блока фильтрации 9 и регистратор тока 3. Стоит отметить, что полосовые фильтры 6 и АЦП 7 для приемной катушки 5 и датчика тока 8 имеют одинаковые параметры.
Алгоритм работы блока фильтрации 9, в заявленном изобретении, поясняется чертежом на фиг. 2, где вычисление амплитуды и фазы полезного сигнала происходит после нахождения весовых коэффициентов целевой функции поэтапно:
1. Подготовка регистрируемых (измеряемых) данных 1;
2. Поиск начального приближения для метода Нелдера-Мида, используя многоканальный адаптивный фильтр RLS 2;
3. Нахождение безусловного минимума целевой функции методом Нелдера-Мида 3;
4. Процентная оценка наличия прямого поля в регистрируемых данных 4;
5. Вычисление амплитуды и фазы полезного сигнала 5.
Пример осуществления изобретения
Для апробации предлагаемого технического решения был изготовлен макет, моделирующий работу аппаратуры частотного электромагнитного зондирования, схема макета представлена на Фиг. 3. Макет состоит из генератора гармонических сигналов Г3-112/1; усилителя мощности ЕР4000; усилителя сигнала от приемной катушки, разработанный на основе микросхемы INA163 с коэффициентом усиления 119; а также быстродействующий модуль АЦП Е20-10 (разрядность АЦП 14 бит, частота работы АЦП до 10 МГц) с пределом допускаемой относительной основной погрешности частоты преобразований АЦП - ±0.005%, допускаемой относительной основной погрешности измерений напряжения переменного тока от ±0.2% и с возможностью подключения к компьютеру для управления и сбора данных.
Для проверки способа подавления прямого поля выполнены эксперименты, моделирующие ситуацию расположения прибора в воздухе, когда полезный сигнал отсутствует в регистрируемых данных, и с замкнутым кольцом, лежащим параллельно плоскости макета на глубине 0.22 м, создающим полезный сигнал в приемной катушке. Кольцо сделано из медного провода диаметром 1.5 мм со средним радиусом 0.247 м. Расстояние между центром генераторной катушки и центром приемной катушки 0.605 м. Параметры датчика тока, приемной катушки и генераторной катушки представлены в таблице 1. Далее описываются результаты экспериментов при частоте тока 3 и 10 кГц в генераторной петле и при различной амплитуде тока.
Для поиска весовых коэффициентов заданной целевой функции подготовка измеренных данных и вычислить следующие значения:
1. Амплитуду тока в генераторной петле по формуле 7.
2. Время существования сигнала и шаг дискретизации:
где N - количество периодов,
где m - количество точек данных.
В качестве алгоритма фильтрации используется многоканальный адаптивный фильтр на основе RLS-алгоритма [7, 13]. Выбор алгоритма рекурсивных наименьших квадратов обуславливается тем, что входные сигналы считаются детерминированными, а также данный алгоритм демонстрирует быструю сходимость по сравнению с большинством своих конкурентов [7]. Многоканальность достигается за счет того, что теоретический сигнал математического моделирования представляется в виде двух сигналов, а именно: cos(ω⋅t) и - sin(ω⋅t).
Стоит отметить, что нахождение минимального значения невязки в целевой функции (9) может быть выполнено разными адаптивными алгоритмами, однако, как отмечалось ранее, в установившемся режиме значения весовых коэффициентов адаптивных фильтров могут находиться в окрестности истинного минимума [7]. Поэтому для проверки и уточнения ВК используется метод Нелдера-Мида, также известный как метод деформируемого многогранника и симплекс-метод. Алгоритм позволяет находить безусловный минимум функции нескольких переменных, не используя производной (точнее - градиента) функции, а поэтому может применяться к негладким и/или зашумленным функциям [14]. Последовательность шагов для реализации алгоритма метода Нелдера-Мида подробно описана в [15]. Суть метода заключается в последовательном перемещении и деформировании симплекса вокруг точки экстремума и требует начального приближения. Поэтому, после реализации адаптивного фильтра по полученным значениям коэффициентов, формируется симплекс, в соответствии с методом Нелдера-Мида, и далее используется для нахождения безусловного минимума функции (9). После нахождения весовых коэффициентов рассчитывалась процентная оценка сигнала прямого поля в регистрируемых данных и определялись амплитуда и фаза полезного сигнала.
При размещении прибора в воздухе в отсутствии замкнутого кольца, приемная катушка регистрировала только сигнал прямого поля, соответственно, согласно логике алгоритма, при заданной оценке наличия прямого поля (помехи) в регистрируемых данных больше 95%, характеристики полезного сигнала приравнивались к нулю. Так для частоты 3000.4 Гц коэффициент A1=8.662⋅10-3, а амплитудное значение прямого поля по формуле (8) равняется 8.66⋅10-3, соответственно процентная оценка наличия прямого поля равна 99.98%, и для частоты 10003.2 Гц, амплитудное значение прямого поля равно 9.101⋅10-3, коэффициент A1=9.152⋅10-3, и тогда, на этой частоте, прямое поле удалось подавить на 99.44%.
Для сравнения, при замкнутой петле на частоте 3000.4 Гц, коэффициент A1=8.75⋅10-3, а амплитудное значение прямого поля равно 8.28⋅10-3, тогда процентная оценка наличия прямого поля равна 94.63%, что указывает на наличие полезного сигнала в данных с приемной катушки. Также для частоты 10003.2 Гц коэффициент A1=9.22⋅10-3 и амплитудное значение прямого поля равно 8.67⋅10-3, соответственно оценка наличия прямого поля равна 93.93%.
Данные регистрировались с частотой дискретизации 500 кГц. На каждой рабочей частоте записывалось десять периодов с синхронизацией по каналу сигнала от датчика тока. Результаты фильтрации сравнивались с теоретическим расчетом магнитных полей с дипольными источниками и приемниками [3, 8, 10, 11, 16, 17] (таблица 2).
Литература
1. Способ индукционного частотного зондирования: пат.RU №2152058 С1 / А.К. Манштейн, М.И. Эпов, В.В. Воевода, К.В. Сухорукова; заявитель и патентообладатель Институт геофизики СО РАН; опубл. 27.06.2000, Бюл. №18. 4 с.
2. Способ и устройство для индукционного частотного зондирования: пат.RU №2502092 С2 / А.К. Манштейн, Е.В. Балков; заявитель и патентообладатель Институт нефтегазовой геологии и геофизики СО РАН; опубл. 20.12.2013, Бюл. №35.
3. Способ калибровки устройства для наземного электромагнитного частотного зондирования: пат.RU №2461850 С2/ А.К. Манштейн, М.И. Эпов, Е.В. Балков, К.В. Сухорукова; заявитель и патентообладатель Институт нефтегазовой геологии и геофизики СО РАН; опубл. 20.09.2012 Бюл. №26. 4 с.
4. Балков Е.В., Манштейн А.К. Сравнение характеристик двух- и трехкатушечной реализации индукционных зондов для малоглубинного частотного зондирования // Геофизический вестник. - 2006. - №1. - С. 12-17.
5. Фадеев Д.И., Балков Е.В., Панин Ю.Г., Карин Ю.Г. Калибровка и экспресс-трансформация сигналов компактной аппаратуры малоглубинного электромагнитного профилирования // Геофизика. - 2018. - №2. - С. 52-57.
6. Magnetic compensation circuit and method for compensating the output of a magnetic sensor, responding to changes in a first magnetic field: pat. US Pub. No.: 9389281 B2 / G. F. West, P. W. Walker, B. D. Polzer; Application Vale S. A.; pub. Date 12.07.2016. 12 p.
7. Джиган В.И. Адаптивная фильтрация сигналов: теория и алгоритмы. - М.: техносфера, 2013. - 528 с.
8. Параснис Д.С. Принципы прикладной геофизики; перевод с английского С.Д. Виноградова; под редакцией В.Н. Жаркова. - М.: Мир, 1965. - 200 с.
9. Митрофанов Г.М. Обратные задачи геофизики. - Новосибирск: РИЦ НГУ, 2015. - 101 с.
10. Кауфман А.А. Введение в теорию геофизических методов. Часть 2. Электромагнитные поля. - М.: Недра, 2000. - 483 с.
11. Яворский Б.М., Детлаф А.А. Справочник по физике: 3-е изд., испр. - М.: Наука, 1990. - 624 с.
12. Корпусов М.О., Панин А.А. Лекции по линейному и нелинейному функциональному анализу. Том 3. Нелинейный анализ. - М.: Физический факультет МГУ, 2016. - 235 с.
13. Джиган В.И. Параллельные регуляризированные RLS - алгоритмы многоканальной адаптивной фильтрации // Цифровая Обработка Сигналов. - 2004. -№2. - С. 7-13.
14. Гладилин П.Е., Боченина К.О. Технологии машинного обучения. - СПб.: Университет ИТМО, 2020. - 75 с.
15. Банди Б. Методы оптимизации. Вводный курс: перевод с английского. -М.: Радио и связь, 1988. - 128 с.
16. Тамм И.Е. Основы теории электричества. - М.: Наука, 1966. - 624 с.
17. Балков Е.В., Манштейн А.К. Трехкатушечный индукционный зонд в частотном зондировании // Геофизический вестник. - 2001. - №12. - с. 17-20.
Claims (2)
1. Способ подавления сигнала прямого поля при электромагнитном индукционном зондировании, включающий измерение суммы сигналов прямого поля и полезного сигнала от среды в приемной катушке и сигнала датчика тока, установленного внутри генераторной катушки и дающего информацию о прямом поле, полосовую фильтрацию и оцифровку измеренных данных, отличающийся тем, что подавление сигнала прямого поля выполняют в блоке фильтрации с использованием многоканального адаптивного RLS-алгоритма путем определения весовых коэффициентов при минимизации невязки целевой функции, характеризующей разность между измеренным и моделируемым сигналами с приемной катушки, формированием по полученным значениям коэффициентов симплекса в соответствии с методом Нелдера-Мида для нахождения безусловного минимума целевой функции, процентной оценкой наличия сигнала прямого поля в регистрируемых данных и вычислением амплитуды и фазы полезного сигнала на основе весовых коэффициентов.
2. Устройство подавления сигнала прямого поля при электромагнитном зондировании, состоящее из генераторной катушки, генератора гармонических сигналов, датчика тока, установленного внутри генераторной катушки и дающего информацию о прямом поле, а также одной или набора приемных катушек, сигналы от которых через соответствующие цепи, имеющие полосовой фильтр и аналого-цифровой преобразователь, соединены с блоком фильтрации, где амплитуда и фаза полезного сигнала вычисляются с использованием способа подавления сигнала прямого поля по п. 1.
Publications (1)
Publication Number | Publication Date |
---|---|
RU2818011C1 true RU2818011C1 (ru) | 2024-04-23 |
Family
ID=
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU702334A1 (ru) * | 1976-01-22 | 1979-12-05 | Государственный Геофизический Трест "Татнефтегеофизика" | Измерительное устройство дл геоэлектроразведки |
RU2152058C1 (ru) * | 1998-06-24 | 2000-06-27 | Институт геофизики СО РАН | Способ индукционного частотного зондирования |
UA32855U (ru) * | 2007-06-27 | 2008-06-10 | Общество С Ограниченной Ответственностью "Юг - Нефтегазгеология" | Устройство для удаления сильных узкополосных помех в приемниках естественного электромагнитного излучения Земли |
US20090295391A1 (en) * | 2008-06-02 | 2009-12-03 | Geonics Limited | Combined Electromagnetic Sensor and Magnetometer |
US20100052685A1 (en) * | 2008-08-29 | 2010-03-04 | Geotech Airborne Limited | Bucking coil and b-field measurement system and apparatus for time domain electromagnetic measurements |
CN103424780B (zh) * | 2013-08-27 | 2016-05-18 | 中国航空无线电电子研究所 | 基于感应线圈的载机磁环境补偿方法 |
US9389281B2 (en) * | 2013-03-21 | 2016-07-12 | Vale S.A. | Magnetic compensation circuit and method for compensating the output of a magnetic sensor, responding to changes in a first magnetic field |
RU2629705C1 (ru) * | 2016-04-19 | 2017-08-31 | Закрытое акционерное общество "Аэрогеофизическая разведка" | Способ высокоточных электромагнитных зондирований и устройство для его осуществления |
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU702334A1 (ru) * | 1976-01-22 | 1979-12-05 | Государственный Геофизический Трест "Татнефтегеофизика" | Измерительное устройство дл геоэлектроразведки |
RU2152058C1 (ru) * | 1998-06-24 | 2000-06-27 | Институт геофизики СО РАН | Способ индукционного частотного зондирования |
UA32855U (ru) * | 2007-06-27 | 2008-06-10 | Общество С Ограниченной Ответственностью "Юг - Нефтегазгеология" | Устройство для удаления сильных узкополосных помех в приемниках естественного электромагнитного излучения Земли |
US20090295391A1 (en) * | 2008-06-02 | 2009-12-03 | Geonics Limited | Combined Electromagnetic Sensor and Magnetometer |
US20100052685A1 (en) * | 2008-08-29 | 2010-03-04 | Geotech Airborne Limited | Bucking coil and b-field measurement system and apparatus for time domain electromagnetic measurements |
US9389281B2 (en) * | 2013-03-21 | 2016-07-12 | Vale S.A. | Magnetic compensation circuit and method for compensating the output of a magnetic sensor, responding to changes in a first magnetic field |
CN103424780B (zh) * | 2013-08-27 | 2016-05-18 | 中国航空无线电电子研究所 | 基于感应线圈的载机磁环境补偿方法 |
RU2629705C1 (ru) * | 2016-04-19 | 2017-08-31 | Закрытое акционерное общество "Аэрогеофизическая разведка" | Способ высокоточных электромагнитных зондирований и устройство для его осуществления |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5999884A (en) | Method for determining the axial position of formation layer boundaries using measurements made by a transverse electromagnetic induction logging instrument | |
US9715034B2 (en) | Method for multi-tubular evaluation using induction measurements | |
US6044325A (en) | Conductivity anisotropy estimation method for inversion processing of measurements made by a transverse electromagnetic induction logging instrument | |
US5065099A (en) | Coil array for a high resolution induction logging tool and method of logging an earth formation | |
US5666057A (en) | Method of skin effect correction and data quality verification for a multi-frequency induction well logging instrument | |
EA012792B1 (ru) | Способ фазовой и амплитудной коррекции данных электромагнитного исследования с управляемым источником | |
CN109085651B (zh) | 一种探测井下套管损伤的方法和系统 | |
Poliakov et al. | The range of induction-coil magnetic field sensors for geophysical explorations | |
Tasič et al. | Seismometer self-noise estimation using a single reference instrument | |
Chen et al. | Weak coupling technology with noncoplanar bucking coil in a small-loop transient electromagnetic system | |
CN109814163A (zh) | 一种基于扩展补偿模型的航磁张量数据抑噪方法及系统 | |
Davis et al. | Quantifying AEM system characteristics using a ground loop | |
Mester et al. | Development and drift-analysis of a modular electromagnetic induction system for shallow ground conductivity measurements | |
Ambruš et al. | Automatic compensation of primary field coupling for a frequency-domain electromagnetic induction sensor | |
RU2818011C1 (ru) | Способ подавления сигнала прямого поля в электромагнитном индукционном зондировании и устройство для его реализации | |
CN110273675B (zh) | 一种瞬变电磁差分测井方法及其系统 | |
Kolobov et al. | The KVVN-7 multifunction digital measuring station for electromagnetic monitoring of seismoactive zones | |
US5736860A (en) | Inductive logging system having multiple coils from which are selected a receiver and a counter-receiver | |
US5500597A (en) | Multiple transmit frequency induction logging system with closed loop conversion circuitry for phase and gain variation elimination | |
USRE47622E1 (en) | High-sensitivity subsurface sensing system | |
Bryakin et al. | Electromagnetic Acceleration Transducer with Eddy Current Transformation | |
RU2152058C1 (ru) | Способ индукционного частотного зондирования | |
RU2461850C2 (ru) | Способ калибровки устройства для наземного электромагнитного индукционного частотного зондирования | |
Larchenko et al. | Three-component measurements of the structure of the electromagnetic field in the extremely-low-frequency and ultralow-frequency bands | |
Zhou et al. | A calibration method of ground-airborne frequency domain electromagnetic receving system based on magnetic source excitation |