CC BY
93
© Ю.Б. Башкуев, Д.Г. Буянова, 2011
УДК 550.837.7
Ю.Б. Башкуев, Д.Г. Буянова
ЭЛЕКТРИЧЕСКИЕ СВОЙСТВА ГОРНЫХ ПОРОД АЛДАНО-СТАНОВОГО НАГОРЬЯ ПО ДАННЫМ СВЕРХДЛИННОВОЛНОВОГО РАДИОИМПЕДАНСНОГО ЗОНДИРОВАНИЯ С БОРТА САМОЛЕТА
Рассмотрены статистические характеристики эффективного сопротивления р_ различных комплексов горных пород Алдано-Станового нагорья в сверхдлинноволновом диапазоне радиоволн, полученные методом СДВ радиоимпедансного зондирования с борта самолета в районах развития мерзлоты и кристаллических массивов. Исследованы 6участков площадью 30 тыс. км2.
Ключевые слова: Алдано-Становое нагорье, радиоимпедансное зондирование, сверх-длинные (СДВ) радиоволны, аэросъемка, электрические свойства, статистические характеристики горных пород.
ЖЭ связи со строительством нефте-
-Ш-} провода «Восточная Сибирь-Тихий Океан», началом эксплуатации Амуро-Якутской железной дороги на участке «Беркакит-Томмот» и развитием горно-рудной отрасли Южной Якутии актуальной стала задача изучения электрических свойств горных пород Алда-но-Станового нагорья. Статистические характеристики эффективного сопротивления р_ различных комплексов горных пород необходимы для решения практических задач транспортного строительства, инженерной геофизики, геологии, электроэнергетики и телекоммуникаций.
В литературе имеется очень мало данных по электрическим свойствам горных пород Алдано-Станового нагорья [1]. Цель статьи - дать обобщенное статистическое описание электрических свойств основных комплексов горных пород Алдано-Станового нагорья по данным радиоимпедансных зондирований с борта самолета в сверхдлинновол-
новом (СДВ) диапазоне электромагнитных волн.
В бездорожной горно-лесистой, пустынной и тундровой областях возникает необходимость использования методов определения электрических свойств верхней части земной коры с применением летательных аппаратов. Достоинство этих методов состоит в том, что они дают пространственное распределение эффективного сопротивления р~ и поверхностного импеданса 5 исследованной территории в виде крупномасштабных геоэлектрических карт. Исследования эффективного сопротивления подстилающей среды в СДВ-СВ диапазонах электромагнитных волн с борта самолета выполнены в США, Канаде и России [2, 3].
Аппаратура и методика измерений
Аппаратура представляет трехканальный селективный микровольтметр-фазометр, измеряющий амплитуднофазовую структуру СДВ электромагнитного поля вблизи границы раздела «земля - воздух» [3, 4]. Антенная систе-
ма состоит из двух несимметричных антенн для приема вертикальной Е2 и горизонтальной Ех составляющих электрического поля и рамочной антенны для приема горизонтальной составляющей магнитного поля Ну. В аппаратуре использован метод фазоквадратурной компенсации помехи от вертикальной составляющей Е2 в канале для приема Ех [3]. Однако при этом регистрируется не полный вектор Ех, а только его мнимая часть 1тЕх, и, следовательно, измеряется мнимая часть приведенного поверхностного импеданса 5 = Re5 - Пт5, где 1т8=(1тЕх)/(Ну^0), Zo - характеристический импеданс вакуума. Над однородным проводящим полупространством 1т8 = ReS, ф5 = arctg Im5/Re5 = -45° и эффективное сопротивление р~ вычисляется по формуле:
р~ = 120Я^( 1тEx) / (Hy - 2о )]2
р~ = 120Я^( 1тEx)/(Hy - 2о)]2, (1)
где р~ - эффективное электрическое сопротивление среды, X - длина волны в воздухе, м. Формула для определения р_ справедливы с погрешностью не более 10% при условии 60Х/р_е >>10, е - относительная диэлектрическая проницаемость однородной среды. В частности на частотах 10-30 кГц при е<10 этому условию соответствуют р_ <
18000^6000 Омм. Источниками электромагнитного поля служили СДВ радиостанции, работающие на частотах 17,4 и 22,3 кГц. При съемке использовалась равномерная сеть прямолинейных маршрутов длиной до 150-250 км и с расстоянием между ними от 0,5 до 10,0 км. Полеты выполнялись на высотах от 50 до 400 м со скоростью 140 - 180 км/ч. Калибровка измерительной установки проводилась над однородной в пределах скин-слоя водной средой оз. Байкал при сопротивлении воды р = 150 Ом-м.
Методические погрешности, обусловленные вычислением р_ по ImS при ф5= -45°, а также погрешности при калибровке установки в полете и обработке СДВ сигнала дают общую погрешность определения р_ не хуже ±(10-15) %. Эта оценка подтверждается сравнением данных на контрольных маршрутах и результатами наземных параметрических измерений. Важное достоинство радиоимпедансного зондирования с борта самолета в СДВ диапазоне состоит в его высокой пространственной разрешающей способности при определении геоэлектрических границ между проводящими и слабопроводящими типами горных пород, значительной глубинности и в высокой производительности аэросъемки для целей геоэлектриче-ского картирования и инженерной геофизики в труднодоступных местностях.
Методика обработки результатов
Особенностью геоэлектрического строения Алдано-Станового нагорья является сложное пространственное распределение эффективного сопротивления верхней части земной коры, обусловленное его геолого-тектони-ческим строением. Исходные данные для статистических расчетов получены в результате аэроработ на самолете АН-2 в Южной Якутии (рис. 1) общей площадью свыше 30 тыс. км2.
Под эффективным сопротивлением р~ проводящей неоднородной среды (геоэлектрического разреза) на фиксированной частоте СДВ диапазона мы понимаем сопротивление некоторого фиктивного однородного изотропного полупространства, по
Рис. 1. Схема аэросъемочных работ на самолете АН-2. Республика Саха - Якутия (участки: 1 - Центральный Алдан, 1а - Нимгер-канский, 1б - Канкунский, 2 - Восточный Алдан, 2а - Эльконкский, 3 - Становой)
значению равное его удельному электрическому сопротивлению (УЭС) на данной частоте. Эффективное сопротивление определяется по величине мнимой части импеданса (1). Слой скин-эффекта
нс - 503^^77 (2)
изменяется в широких пределах и на частоте f = 17,4 кГц при = 100-20000 Ом-м составляет 58-540 м, охватывая приповерхностную часть земной коры. Определение р~ проводилось в дискретных последовательностях точек, расположенных на маршруте аэросъемки через 100-400 м.
С целью получения статистических оценок р^ однородной в геологическом отношении области (комплекса горных пород) значения р^ объединялись в вариационные ряды, содержащие от 25 до 1544 величин р^, по которым строились гистограммы распределения р_ Форма гистограммы служила основным качественным критерием для выяснения правильности выделенной петрофизической группы. После уточнения выделенных групп вычислялись обобщенные показатели распределения р^ [5]:
среднее арифметическое значение
__ 1 N
р~ - N § р~';
среднее значение логарифма
_____ 1 N
^ Р~ - Т7 § ^ Р~у;
1У у-1
среднее геометрическое
рс -105*
' с.г.
Сведения о законе распределения эффективного сопротивления р^ в литературе отсутствуют. Для проверки соответствия эмпирического распределения р~ нормальному или логнормальному законам использован приближенный метод, основанный на вычислении показателя асимметрии у3 и эксцесса у4 выборочного распределения [5].
Результаты дистанционного зондирования и их обсуждение
В Южной Якутии исследованы 6 участков общей площадью 30 тыс. км2. Участки сложены метаморфическими и изверженными породами архея, кар-бонатно-терригенными породами кембрия и терригенными отложениями юрского возраста. Толщина покровных рыхлых отложений колеблется от единиц до нескольких десятков метров, удельное сопротивление - от десятков до первых тысяч Ом-м. На всей площади развиты тектонические нарушения различной ориентировки, зоны распространения которых отличаются узкими протяженными областями пониженных значений эффективного сопротивления р~. Участки, сложенные метаморфическими и изверженными породами архея, такие как Эльконк-ский и Становой, имеют высокие значения р~. Области развития кембрийских и юрских отложений, расположенные, в основном, в северной и южной частях исследуемого района, обладают более низкими относительно архейских пород значениями р~
^п Горные породы Рабочая частота, кГц N Р~ /ре.™ Ом-м рмакс/^мин, Ом-м Участок, маршрут, район
1 2 3 4 5 6 7
Кристаллические комплексы
1. Гранитоиды архейские 22,3 1052 6200 4700 4600 0,35 36000 250 Эльконкский
2. Гранитоиды архейские 22,3 50 4500 3300 3500 0,32 15000 990 Становой, маршрут М3
3. Гранитоиды архейские с непрерывным чехлом рыхлых отложений 17,4 1544 3600 3700 0,42 31700 Центральный Алдан
2500 160
Осадочные комплексы
4. Терригенные угленосные отложения юрского возраста 17,4 602 750 500 730 0,45 5400 50 Чульманская впадина
5. Карбонатно-терригенные отложения кембрия 22,3 549 1160 840 940 0,38 7400 80 Северный склон Алданского нагорья
В геоэлектрике обычно принято исследовать электрические свойства горных пород по регионам с выделением типичных комплексов горных пород. Существенным затруднением при систематизации горных пород по электрическим свойствам является то обстоятельство, что их петрографическая классификация составлена по признакам (генезис, химизм, структура и др.), которые не находятся в прямой связи с их электрическими свойствами. В связи с этим описание горных пород (таблица) дано без излишней детализации. Группирование по комплексам пород выполнено по геологическим картам. Результаты определения р_, полученные на переходных участках от одного комплекса пород к другому, не рассматривались. Участки аэроработ расположены на Ал-дано-Становом нагорье (рис. 1). В целом на всех участках наблюдается очень
широкий диапазон изменения р~: от 20 до 30 000 Ом-м и более. Установлена значительная дифференциация р~ между осадочными и кристаллическими комплексами горных пород (таблица). Для горных районов, сложенных кристаллическими породами, средние значения р_ составляют 4000-6000 Ом-м и выше, тогда как для межгор-ных впадин, заполненных осадочными породами, средние значения р_ обычно ниже 1000-2000 Ом-м. Электрические границы различных комплексов горных пород довольно точно совпадают с геологическими границами.
Распределения р_ основных комплексов горных пород региона положительно асимметричные, чаще всего крутые слева и плоские справа (рис. 2). Распределение р~ ограничено условием р_ >
0, а изменения р_ происходят в широких пределах и достигают нескольких по-
рядков (от 0,25 до 105 кОм-м). Распределение УЭС в однотипных комплексах горных пород обычно подчиняется логнормальному закону.
Полученные статистические характеристики позволяют оценить степень геоэлектрической неоднородности однотипных в геологическом отношении комплексов, а также степень отличия между различными комплексами горных пород. Для большинства рассмотренных в таблице комплексов кристаллических пород закон распределения р^ соответствует логнормальному, например, для карбонатно- терригенных пород кембрия северного склона Алданского щита.
Величина рсг. в случае логнормального распределения случайной величины р~ является медианой этого распределения, оценки среднего значения р~ и моды определяются по формулам [5]:
— = 10(^ +и5134
-)
то— = 10('8— -130264—■).
(3)
Р,%
20-
10-
N=597
й
Следует отметить, что среднее геометрическое не является оценкой неизвестного математического ожидания, так как при увеличении числа наблюдений оно стремится к величине меньшей, чем математическое ожидание. Среднее же арифметическое при увеличении числа наблюдений в пределе приближается к неизвестному значению математического ожидания, независимо от вида функции распределения случайной величины. На этом основании в таблице приводятся значения р~ и рсг. В частности, для рассмотренных в табл. 1 карбо-натно-терригенных пород кембрия расчетные по формулам (3) величины совпадают с табличными.
Анализ таблицы и гистограмм на рис. 2 - 3 показывает, что кристаллические горные породы характеризуются
'по ‘о оо Оо
З' ^ о сЧ»-н' г\Г ГП КГ
О О О ’—н
р^кОм-м
Рис 2. Гистограмма распределения р~ кристаллических пород кембрия северного склона Алданского щита на частоте 22,3 кГц
высокими фоновыми значениями р~, превышающими 1-2 тыс. Ом-м. Низкие значения р~ обусловлены влиянием вы-сокопроводящих осадочных пород верхней части разреза и связаны с наличием чехла рыхлых покровных отложений и коры выветривания толщиной от единиц до нескольких десятков метров (рис. 4). Сравнительно низкие сопротивления р~ некоторых кристаллических горных пород обусловлены также наличием проводящих включений (графити-стые сланцы, магнетитовые породы и др.). Обнаженные участки кристаллических пород, соответствующие обычно возвышенным частям рельефа, имеют, как правило, высокие значения р~, характерные для кристаллического фундамента. В целом можно отметить рост средних значений р~ архейских грани-тоидов от 2500 до 4700 Ом-м с увеличением широты места участка работ и степени расчлененности рельефа. Причиной может быть появление развитой криолитозоны (климатический фактор) и понижение минерализации поверхностных и грунтовых пород.
Р’°/о ґ = 17,4 кГц И= 1544
р,%
30-
20-
10-
0
^ *0 Г\| ^ оо Ю N
О гп
р_,кОм*м
а
N=408
15 -
3<
10 -
2(
5 -
1( 0
(
ГМ ІЛ-НІЛГІГО^ЮОО^
Гч|
р^, кОм • М
Рис. 5. Гистограмма р~ маршрута 31 на ЦентральноМАЛдОНй * М
б
Рис. 3. Гистограммы эффективного сопротивления р~: а - гранитоиды архейские с непрерывным чехлом рыхлых отложений (участок Центральный Алдан); б - гранитоиды архейские (участок Эльконкский)
На некоторых участках гранитоид-ных массивов, перекрытых чехлом рыхлых покровных отложений, получены бимодальные гистограммы распределения р^, на которых частота меньшего модального интервала соизмерима с частотой большего модального интервала (рис. 5). Бимодальные гистограммы, полученные в результате определения р^ на длинных маршрутах, проходящих над однотипным комплексом пород, позво-
ляют дать оценки параметров среднестатистического двухслойного геоэлек-трического разреза. Для гистограммы р^ рис 5 значения р1 и р2, определенные по средней величине модального интервала, будут равны 750 и 3500 Ом-м. Толщина первого слоя ^ может быть определена по величине скин-слоя на рабочей частоте 22,3 кГц по формуле (2). Она будет изменяться случайным образом от 0 до 92 м. В качестве значения ^ условно можно принять четверть величины скин-слоя ^=23 м. Такая приближенная оценка параметров геоэлектрического разреза может быть полезна при проверке прогнозных карт параметров слоистого геоэлектрического разреза и расчетах поверхностного импеданса слоистой среды.
Заключение
1. Проведена классификация и получены стати-
Рис.4. Гистограмма распределения р~ осадочных пород Алданского нагорья
стические характеристики эффективного сопротивления для различных комплексов горных пород Алдано-Станового нагорья в СДВ диапазоне радиоволн. Установлена значительная дифференциация эффективного сопротивления между осадочными и кристаллическими породами, выявлено хорошее совпадение электрических границ различных комплексов горных пород (гранитоидов, кембрийских и юрских отложений и др.) с геологическими границами. Для однотипных кристаллических и осадочных пород статистическое распределение р~ соответствует логнормальному закону.
2. Для горных районов, сложенных кристаллическими породами, средние значения — составляют 2500-4700 Ом-м, тогда как для впадин, образованных осадочными породами, средние значения — обычно ниже 1000 Ом-м. Наиболее высокими значениями — (>8000 Ом-м) характеризуются выходы
кристаллических пород архея, не перекрытые слоем наносов, а также узкие речные долины в северной и северовосточной Алданского нагорья. Резкое увеличение р~ на узких речных долинах можно объяснить, если предположить, что в их пределах распространены многолетнемерзлые породы.
3. Использованный в экспериментах метод дистанционного зондирования позволяет быстро исследовать большие территории в труднодоступной местности, дает возможность оперативно составлять карты эффективного сопротивления крупных регионов, которые могут использоваться при инженерно-геологических изысканиях. Перспективен данный метод и в мерзлотоведении. Результаты аэроработ использованы для решения практических задач инженерной геофизики, геологии и транспортного строительства.
Статья подготовлена при финансовой поддержке грантов РФФИ № 08-02-98007, 0905-98611 и интеграционного проекта СО РАН №56.
1. Сикорский В.А. Измерение сопротивления (р) и диэлектрической проницаемости (s) среды радиоимпедансным методом (Алданский щит) // Геология и геофизика. - 1977. -№ 10. -P. 95-101.
2. King R. J. On airborne wave tilt measurements // Radio Science. - 1977. - Vol. 12, №
3. - P. 405-414.
3. Вешев А. В., Ивочкин В. Г., Яковлев А. В. и др. Аэроэлектроразведочная аппаратура
--------------- СПИСОК ЛИТЕРАТУРЫ
РЭМП и результаты ее полевых испытаний // Геофизическая аппаратура. - Л., 1981. - Вып. 73. - С. 77-87.
4. Башкуев Ю.Б., Дембелов М.Г., Буянова Д.Г. Высотное изменение импеданса подстилающей среды.// Сб. "Солнечно-земная физика", 2004, вып. 5 (118), с. 147-149.
5. Каждан А. Б., Гуськов О. И. Математические методы в геологии. - М.: Недра, 1990. - 252 с. шгЛ
— Коротко об авторах ----------------------------------------------------------------
Башкуев Ю.Б. - профессор, доктор технических наук, заведующий лабораторией геоэлектромагнетизма Отдела физических проблем Бурятского научного центра Сибирского отделения РАН, buddich@mail.ru
Буянова Д.Г. - кандидат физико-математических наук, доцент, ведущий научный сотрудник лаборатории геоэлектромагнетизма Отдела физических проблем Бурятского научного центра Сибирского отделения РАН, dbuy@mail.ru
------------------------------------------------------ РУКОПИСИ,
ДЕПОНИРОВАННЫЕ В ИЗДАТЕЛЬСТВЕ
МОСКОВСКОГО ГОСУДАРСТВЕННОГО ГОРНОГО УНИВЕРСИТЕТА
Голиков А.В. - инженер-программист, НФККИПКРО, e-mail: ap4u@ap4u.ru
МЕТОДИКА АВТОМАТИЧЕСКОГО ПРЕОБРАЗОВАНИЯ УЗЛОВЫХ МОДЕЛЕЙ В ЗАДАЧАХ КОНЕЧНОЭЛЕМЕНТНЫХ МОДЕЛЕЙ ПРИ КОНЕЧНОЭЛЕМЕНТНОЙ АППРОКСИМАЦИИ В ЗАДАЧАХ ПОИСКА РУД (784/01-11 от 20.10.2010) 16 с.
В настоящее время методы электроразведки широко применяются в исследовании геологического строения земной коры и поисках месторождений полезных ископаемых. Для расчета электромагнитного поля в областях, имеющих сложную структуру и трехмерные неоднородности, наиболее эффективным является. Метод конечных элементов (МКЭ). Актуальной проблемой» является создание быстрых алгоритмов построения конечноэлементных сеток и удобных для реализации вычислительных схем при работе на постоянном токе и по методу ВП для решения задачи поисков и разведки залежей углеводородов и рудных объектов в сложных геологических условиях.
Ключевые слова: метод конечных элементов, разведка, залежей углеводородов.
Golikov A. V. THE METHOD OF DATA CONVERSION OF THE KEY MODELS FOR SOLVING PROBLEMS OF APPROXIMATION DURING ORE SURVEY
Currently the methods of the electrical survey are widely used in the studies of the geological structure of the Earth crust andfor the exploration of mineral deposits. The finite element method (FEM) is the most effective one for calculation of parameters of the electromagnetic field in areas with the complex structure and three-dimensional heterogeneity. The topical problem today is the creation of fast algorithms for constructing the finite-element nets and easy to implement numerical schemes while working at DC, and using the overvoltage method to solve the problem of prospecting and exploration for hydrocarbon deposits and ore objects in complex geological conditions.
Key words: The finite element method, exploration of hydrocarbon deposits.
