Способ геологической разведки минералов

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

 

Изобретение относится к области дистанционного мониторинга подстилающей поверхности аэрокосмическими средствами.

Известны геопатогенные зоны Земли, так называемые места силы [см. Internet, Wikipedia, Аномальные зоны], в которых происходят необъяснимые явления.

Существует множество методов и средств мониторинга объектов путем их дистанционного зондирования для приема собственного излучения объектов либо отраженного от них электромагнитного поля (солнечного потока).

Информационные признаки собственного излучения подстилающей поверхности, такие как мощность восходящего излучения, его поляризация, диапазон частот, связаны с механическими напряжениями в земной коре и другими аномалиями: изменение удельного сопротивления грунта, плотности теллурических токов, магнитной напряженности. Известно [см. Физический энциклопедический словарь, под ред. A.M. Прохорова, изд. «Сов. Энциклопедия», М, 1983 г., стр. 827], что при механических напряжениях среда становится анизотропной, величина анизотропии пропорциональна механическим напряжениям.

Известен «Способ обнаружения очагов землетрясений», Патент RU №2181495, 2002 г. - аналог. В способе-аналоге осуществляют регистрацию собственного излучения подстилающей поверхности, преобразуют зарегистрированную функцию электрического сигнала в цифровые матрицы |m×n| отсчетов зависимости амплитуды I(х, у) от пространственных координат, выделяют, методами пространственного дифференцирования, контуры на изображении, отличающиеся тем, что собственное излучение регистрируют в двух, взаимно ортогональных по поляризации, плоскостях, формируют результирующую матрицу изображения из попиксельных отношений амплитуд сигнала в двух, взаимно ортогональных по поляризации, каналах приема, выделяют контуры на результирующем изображении, вычисляют функцию фрактальной размерности внутри выделенных контуров, фиксируют очаг землетрясения при совпадении фрактальной размерности участка текущего контура с эталонным, либо при отклонениях на величину не более пороговой.

Недостатками аналога являются:

- глубина проникновения СВЧ излучения (аналога) в диапазоне удельных сопротивлений почвогрунтов не превышает 1 м, в то время как для разведки минералов, глубина проникновения должна быть несколько километров;

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

Ближайшим аналогом к заявляемому техническому решению является «Коррелятор сигналов - предвестников землетрясений», Патент RU №2272306, 2006 г., G01V 9/00.

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

Недостатком ближайшего аналога являются:

- невозможность непосредственного использования из-за физического различия измеряемых величин;

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

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

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

Изобретение поясняется чертежами, где:

фиг. 1 - дискретные отсчеты восходящего низкочастотного излучения вида пород а) жильный кварц, б) агрегаты из зернистого кварца;

фиг. 2 - интервал канального поддиапазона (ΔF) соответствующий полосе захвата автогенератора;

фиг. 3 - функции взаимной корреляции разведываемого сигнала с эталонными сигналами вида пород в) жильный кварц, г) агрегаты из зернистого кварца, д) мелкозернистый кварц сплошной массы;

фиг. 4 - контурная карта разведанного месторождения;

фиг. 5 - функциональная схема устройства, реализующая способ.

Техническая сущность заявленного технического решения состоит в следующем.

Одним из самых распространенных породообразующих минералов является кварц (SiO₃) и его модификация - гранит (в переводе с итальянского - зернистый) [см., например, Советский энциклопедический словарь, под ред. A.M. Прохорова, стр. 570 Кварц, стр. 338 Гранит, изд. Сов. Энциклопедия, 1989 г.] Кварц существует в нескольких структурных формах: зерна, зернистые кристаллы, агрегаты, жильный кварц, прожилки и сплошные массы. Жильный кварц находится в наибольшем геологическом родстве почти со всеми минералами, особенно с золотом (кварцевые пески), кремнезем содержит оливин, пироп, нередко в россыпях встречаются платина, иридий, титановые руды, а также, алмазы (кимберлитовые трубки), рубины и др.

С электрической точки зрения, кварц является пьезоэлектриком, при механических напряжениях на его противоположных гранях появляются заряды противоположного знака, величина которых пропорциональна механическим напряжениям [см., например, Справочник по радиоэлектронике, том 2, под ред. А.А. Куликовского, изд. Энергия, М., 1968 г., стр. 43, Кварцевые резонаторы]

В аномальных зонах тектонических напряжений из-за волн сейсмического фона и приливных волн изменяется плотность электрических зарядов. Кварцевые жилы становятся излучателями низкочастотного восходящего излучения.

Чем больше механические напряжения, тем больше анизотропия пород [см., A.M. Прохоров, Физический энциклопедический словарь, изд. Сов. Энциклопедия, 1983 г., стр. 227]. Анизотропия пород пропорциональна мощности и поляризации восходящего низкочастотного излучения. Частота излучения зависит от толщины пласта и структуры кварцевых зерен. Аномальные изменения электрических параметров литосферы, как то плотности теллурических токов, удельной электрической проводимости пород сказываются на параметрах собственного восходящего излучения. В соответствии с приближенными граничными условиями Леонтовича, глубина проникновения (δ) электромагнитного поля в вещество определяется из соотношения:

ƒ- частота электромагнитного поля, Гц;

μ - магнитная проницаемость породы;

g - удельная проводимость породы Ом⋅м³

На основе принципа взаимности [см. Драбкин А.Л., Зузенко В.Л. «Антенно-фидерные устройства», М., Сов. радио, 1964 г., стр. 136…137, Принцип взаимности] характеристики электрической цепи при излучении (передаче) идентичны характеристикам цепи при приеме. Из принципа взаимности можно рассчитать глубину излучающего пласта породы. Глубина залегания пород (кимберлитовых трубок) - порядка нескольких км, удельная проводимость (мрамор, гранит, кварц g ≈5⋅10¹² Ом⋅м³). Расчетная частота собственного излучения, в зависимости от структурных конгломератов кварца, находится в диапазоне от сотен Гц до нескольких кГц.

Радиоразведку сигналов осуществляют посредством многоканального приемника [см., например, С.А. Вакин, Л.Н. Шустов «Основы радиопротиводействия и радиотехнической разведки, изд. Сов. радио, М., 1968 г., стр. 410, рис. 10.21]. Весь диапазон разведываемых частот разбит на несколько поддиапазонов (каналов приема) с полосой слежения ΔF в каждом, равной полосе захвата i-го автогенератора, как это иллюстрируется рисунком фиг. 2. Вместо гетеродина в аналоге, в реализуемом приемнике используют автогенератор, синхронизацию частоты которого осуществляют внешней ЭДС, наводимой в антенном контуре. В результате такого режима, сигнал автогенератора, модулируемый спектром восходящего излучения, оказывается частотно-модулированным. Модулирующую функцию восходящего излучения выделяют частотным детектором. Следующей проблемой является создание антенного контура для приема низкочастотного восходящего излучения. Известно [см., например, А.Л. Драбкин, В.Л. Зузенко «Антенно-фидерные устройства», Сов. радио, М., 1964 г., стр. 142-143, 256-258, 268], что для эффективного приема (излучения), геометрическая длина антенны должна быть соизмерима с длиной волны (порядка λ2…λ/6). Собственная резонансная частота антенного контура определяется соотношением где L, C - индуктивность и емкость контура. При килогерцовой частоте восходящего излучения, параметры контура должны соответствовать Максимальная емкость конденсатора на основе сегнетоэлектриков (BaTiO₃) достигает ≈10^-4 (Ф). Индуктивность соленоида L(Гн) рассчитывают из соотношения где μ - магнитная проницаемость сердечника (типа пермаллой) может достигать ~10⁴, - длина сердечника, S - площадь сечения сердечника, n - количество витков. При количестве витков несколько сотен, длина - несколько см и сечении 0,1 м², требуемые параметры антенного контура могут быть реализованы на существующей технической базе.

Кроме разведки спектра сигналов, измеритель позволяет идентифицировать вид породы. Основным методом идентификации объектов является сравнение с эталоном. Математической процедурой сравнения подобия двух функций служит вычисление взаимной корреляции их сигналов [см., например, Заездный В.М., «Основы расчетов по статистической радиотехнике», Связь-издат, М., 1969 г., стр. 92-95].

При современном развитии вычислительной техники, методы корреляционного анализа легко реализуются специальными математическими программами для ПЭВМ. Предварительно осуществляют преобразование зарегистрированных сигналов из аналогового вида в цифровой посредством АЦП. Текст программы вычисления функции взаимной корреляции.

Результат расчета иллюстрируется графиками (в, г, д) фиг 3. Селектируемым параметром функции взаимной корреляции B(ƒ₁, ƒ₂) служит интервал корреляции (ширина функции) на уровне 0,1 от ее максимального значения. Для быстрофлюктуирующего процесса функция имеет острую форму. Поэтому график фиг. 3 в) может быть идентифицирован как сигнал жильного кварца, интервал корреляции Δτ порядка единиц μс. График фиг. 3 г) имеет интервал корреляции порядка десятков μс и может быть идентифицирован как участок из агрегатов зернистого кварца. График фиг. 3 д) с интервалом корреляции порядка долей секунды, может быть идентифицирован как участок мелкозернистого кварца сплошной массы.

Пример реализации способа.

Заявленный способ может быть реализован по схеме фиг. 5. Функциональная схема содержит аэрокоптер 1 (беспилотный летательный аппарат, с установленным на его борту измерителем 2 восходящего низкочастотного излучения). Трассовая покадровая съемка запланированного участка разведки осуществляется от бортового комплекса управления (БКУ) 3 по радиолинии командного управления 4 наземного комплекса управления 5 с одновременной привязкой регистрограмм измерений по координатам от бортовой аппаратуры потребителей 6 глобальной космической позиционирующей системы ГЛОНАСС.

Измеритель 2 выполнен по схеме многоканального разведывательного радиоприемника, каждый из N каналов которого разнесен по диапазону на полосу ΔF захвата сигнала соответствующим автогенератором и представляет собой последовательно соединенные элементы: специализированный антенный контур 7, автогенератор 8, операционный усилитель 9, электронный ключ 10, частотный детектор 11, аналогово-цифровой преобразователь 12, буферное запоминающее устройство 13. Синхронизацию работы элементов измерителя 2 осуществляют посредством программируемой выборки измерений 14, в которую закладывают телекоммуникационную программу от БКУ 3. Покадровую дискретизацию измерений во времени в каждом канале реализуют посредством электронного ключа 10, выполненного по схеме [см. «Справочник по радиоэлектронным устройствам» под ред. А.А. Куликовского, изд. Энергия, М., 1978 г., стр. 33, Электронные ключи].

После посадки носителя 1, массивы измерений перекачивают на сервер хранения данных 15 и, через устройство ввода 16, передают в Центр тематической обработки 17. Обработку массивов зарегистрированной информации осуществляют на средствах автоматизированного центра 17.

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

Затем осуществляют визуализацию разведанных территорий с привязкой областей к географическим координатам. Предварительно, в оперативное запоминающее устройство процессора центра записывают программы специализированного программного обеспечения MATH CAD. Затем формируют кадры синтезированных матриц из попиксельных значений дискретных отсчетов в ЗУ каждого канала для соответствующих точек координат [см., например, Векторизация элементов матрицы. «Специализированное программное обеспечение MATH CAD 7.0» PLVS изд. 4-е стереотипное, М,, Информ. Изд. «Филинъ» 2000 г., стр. 50-68]. Программной обработкой выделяют контуры на поле синтезированных матриц, как это иллюстрируется рисунком фиг. 4.

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

Способ геологической разведки минералов, включающий измерение восходящего низкочастотного излучения кварцевых пород, содержащих сопутствующие минералы, измерителем, установленном на аэрокоптере по запланированным трассам полета с привязкой регистрограмм измерений по координатам от бортовой аппаратуры потребителей системы ГЛОНАСС; определение частоты и спектра принимаемого излучения многоканальным радиоприемником с несколькими каналами, разнесенными по поддиапазонам на интервал ΔF захвата частоты автогенератора в режиме затягивания, каждый из которых содержит специализированный антенный контур, автогенератор, операционный усилитель, электронный ключ дискретизации отсчетов измеряемой величины во времени, частотный детектор, аналогово-цифровой преобразователь, запоминающее устройство; идентификацию типа сигнала на соответствующий минерал путем расчета функции взаимной корреляции зарегистрированного сигнала с соответствующим эталонным сигналом от типа пород; формирование результирующей матрицы измерений разведанного сигнала из дискретных отсчетов; выделение контуров границ залегания вида пород на синтезированном изображении.