Компьютерные технологии интерпретации гравитационного и магнитного полей в условиях горной местности тема диссертации и автореферата по ВАК РФ 25.00.10, доктор физико-математических наук Долгаль, Александр Сергеевич

• Долгаль, Александр Сергеевич
• доктор физико-математических наук
• 2002

Оглавление диссертации доктор физико-математических наук Долгаль, Александр Сергеевич

ВВЕДЕНИЕ

1. ОПРЕДЕЛЕНИЕ ПОПРАВОК ЗА ВЛИЯНИЕ РЕЛЬЕФА МЕСТНОСТИ ПРИ ГРАВИМЕТРИЧЕСКОЙ И МАГНИТНОЙ СЪЕМКАХ

1.1. Моделирование погрешностей учета влияния рельефа при гравиметрической съемке

1.2. Компьютерная технология определения поправок за влияние рельефа местности при гравиметрической съемке

1.3. Вычисление топопоправок при магнитной съемке

1.4. Моделирование погрешностей учета влияния рельефа при магнитной съемке

1.5. Определение намагниченности горных пород рельефа по наблюденному магнитному полю

1.6. Приближенный способ вычисления топопоправки 5Тр при аэромагнитной съемке

2. АППРОКСИМАЦИИ ГЕОПОТЕЬЩИАЛЬНЬГХ ПОЛЕЙ ЭКВИВАЛЕНТНЫМИ ИСТОЧНЖАМИ

2.1. Краткая характеристика проблемы истокообразной аппроксимации полей

2.2. Теоретические основы аппроксимации полей сеточными эквивалентными моделями

2.3. Истокообразная аппроксимация гравитационного и магнитного полей, заданных в узлах равномерной сети

2.4. Истокообразная аппроксимация гравитационного и магнитного полей, заданных в узлах неравномерной сети

2.5. Истокообразная аппроксимация аномалий ЕЭП

3. КОМПЬЮТЕРНЫЕ ТЕХНОЛОГИИ ОБРАБОТКИ И ИНТЕРПРЕТАЦИИ ДАННЫХ ГРАВИМЕТРИЧЕСКИХ И МАГНИТОМЕТРИЧЕСКИХ ИЗМЕРЕНИЙ В ГОРНОЙ МЕСТНОСТИ

3.1. Краткая характеристика комплексов программ RELMAG и

ЯБЬО^

3.2. Технология учета влияния резкорасчлененного рельефа местности при магнитометрии

3.2.1. Имитационное моделирование: обнаружение и локализация рудоносных интрузий в условиях Норильского района

3.2.2. Талнахский рудный узел (Норильский район)

3.2.3. Участок Огинраваям (п-ов Камчатка)

3.2.4. Дербинская площадь (Восточный Саян)

3.3. Технология учета влияния резкорасчлененного рельефа местности при гравиметрии

4. ПРИМЕРЫ РЕШЕНИЯ ПРАКТИЧЕСКИХ ЗАДАЧ С ПОМОЩЬЮ КОМПЬЮТЕРНЫХ ИНТЕРПРЕТАЦИОННЫХ ТЕХНОЛОГИЙ

4.1. Общая характеристика процесса интерпретации геофизических данных

4.2. Комплексная интерпретация геопотенциальных полей при поисках медно-никелево-платинового оруденения в Норильском районе

4.2.1. Краткая характеристика проблемы интерпретации гравитационного и магнитного полей исследуемой территории

4.2.2. Краткая физико-геологическая характеристика территории

4.2.3. Методика интерпретации геофизических полей

4.2.4. Некоторые геологические результаты

4.3. Локальное прогнозирование золотого оруденения в Ольховско-Чибижекском рудном районе 257 4.3.1. Геолого-геофизические особенности площади, исходные материалы и последовательность их интерпретации

4.3.2. Методика интерпретации геофизических материалов

4.3.3. Основные геологические результаты

4.4. Геологическое редуцирование поля силы тяжести Талнахского рудного узла

4.5. Комплексирование методов решения обратной задачи для пространственной локализации глубокозалегающих интрузий в Норильском районе

4.6. Решение обратной задачи гравиметрии монтажным методом при поисковых работах на Боотанкагском участке

Введение диссертации (часть автореферата) на тему «Компьютерные технологии интерпретации гравитационного и магнитного полей в условиях горной местности»

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

Геологическая эффективность геофизических исследований во многом определяется методологией и технологией извлечения информации из данных полевых геофизических наблюдений. Методам интерпретации аномалий геопотенциальных полей посвящено огромное количество работ, значение которых трудно переоценить, как в теоретическом, так и в практическом плане. По оценкам В.Н. Страхова, вопросами теории и практики интерпретации данных гравиметрии и магнитометрии в XX веке в СССР и России занималось около 550 исследователей [151], в число которых входят: М.А. Алексидзе, Е.М. Ананьева, Б.А. Андреев, СВ. Аплонов, В.И. Аронов, М. Е. Артемьев, Б.М. Афанасьев, Н.Л. Афанасьев, Б.К. Балавадзе, П.И. Балк, Т.В. Балк, В.И. Бауман, А.Н. Бахвалов, И.М. Бахурин, Ю.И. Блох, H.A. Беляевский, В. А. Болдырева,

A. A. Борисов, В.В. Бродовой, М.А. Бродский, Ю.П. Булашевич, Е.Г. Булах, H.A. Булгаков, Г.М. Валяшко, Ю.Я. Ващилов, Г.С Вахромеев, К.Е. Веселов, О.В. Витвицкий, Г.М. Воскобойников, А.Г. Гайнанов, Г.А. Гамбурцев, И.С Гельфанд, О.Б. Гинтов, К.В. Гладкий, В.Н. Глазнев, В.В. Глазунов, В.Б. Гла-ско, Ю.Н. Годин, Г.Я. Голиздра, В.И. Гольдшмидт, Ф.М. Гольцман, В.М. Гордин, A.M. Городницкий, Г.И. Гринкевич, A.A. Грознова, А.Ю. Давыденко,

B. Л. Данилов, В.М. Девицын, В.А. Дядюра, М.С Жданов, H.A. Журавлев, А.И.

Заборовский, В.Н. Завойский, A.A. Заморев, Ю.А. Зорин, В.К. Иванов, H.A. Иванов, СИ. Иванов, Н.И. Идельсон, В.М. Исаков, М.К. Кабан, Т.Б. Калинина, A.M. Карасик, Г.И. Каратаев, K.M. Картвелишвили, Н.Г. Кассин, И.Г. Клушин, А.И. Кобрунов, А.Б. Коган, Л.А. Коган, В.Г. Козленко, Г.Н. Константинов, И.Н. Корчагин, Г.Г. Кравцов, С.С. Красовский, З.А. Крутиховская, П.Н. Кропоткин, A.B. Кудря, Г.К. Кужелов, М.М. Лаврентьев, П.П. Лазарев, М.И. Лапина, O.K. Литвиненко, A.A. Логачев, В.В. Ломтадзе, П.И. Лукавченко, А.И. Лучицкий, A.A. Любимов, E.H. Люстих, В.А. Магницкий, А.К. Маловичко, Н.Р. Малкин, A.C. Маргулис, М.Н. Маркова, П.С. Мартышко, В.Р. Мелихов, Д.С. Миков, В.О. Михайлов, Т.Л. Михеева, Л.Н. Морозов, Е.А. Мудрецова, Л.И. Надежка, Н.В. Неволин, Л.Д. Немцов, A.A. Непомнящих, A.A. Никитин, n.M. Никифоров, Ф.И. Никонова, П.С. Новиков, В.М. Новоселицкий, Б.В. Нумеров, E.H. Нусипов, A.B. Овчаренко, СМ. Оганесян, Н.И. Павленкова, В.И. Павловский, И.К. Пашкевич, Л.Г. Перфильев, A.A. Попов, А.И. Прилепко, И.Л. Пруткин, И.М. Пудовкин, Е.П. Пучков, А.Б. Раевский, И.М. Рапопорт, Г.Г. Ремпель, E.H. Розе, В.М. Рыбалко, Т.В. Романюк, И.Д. Савинский, К.А. Савинский, М.У. Сагитов, М.А. Садовский, Н.И. Свитальский, М.Г. Сербулен-ко, CA. Серкеров, Т.Н. Симоненко-Розе, З.М. Слепак, A.C. Случановский, П.К. Соболевский, К.П. Соколов, К.И. Соколовский, O.A. Соловьев, Л.В. Сорокин, Л.Н. Сретенский, В.И. Старостенко, В.И. Страхов, В.Н. Страхов, Н.К. Ступак, СИ. Субботин, B.C. Сурков, Б.Б. Тальвирский, О. Л. Тарунина, М.А. Телепин, А.Н. Тихонов, Г.А. Трошков, К.Ф. Тяпкин, Д.Г. Успенский, CA. Ушаков, В.В. Федынский, В.А. Филатов, В.Г. Филатов, Л.Е. Фильштинский, Э.Э. Фотиади, СВ. Шалаев, O.A. Шванк, Э.Л. Шен, Г.Ш. Шенгелая, Н.Г. Шмидт, В.И. Шрайбман, И.И. Шульман, Б.Э. Хесин, И.О. Цимельзон, A.B. Цирульский, В.Г. Чередниченко, A.A. Чернов, А.Б. Черный, H.A. Чуйкова, Ю.В. Юнаковская, A.A. Юньков, H.A. Якимчук, Б.М. Яновский, А.Я. Ярош и другие.

Разработаны различные методы (алгоритмы) интерпретации аномальных потенциальных полей, сформулированы методологические основы создания интерпретационных технологий и разработаны использующиеся в производственной практике интерпретационно - обрабатывающие комплексы программ (системы), ориентированные на решение различных геологических задач. Однако, объективно существующий определенный разрыв между теорией и практикой интерпретации данных гравиразведки и магниторазведки влечет за собой постоянное развитие новых компьютерных технологий, необходимых для решения конкретных геологических задач. Большинство широко использующихся в настоящее время рабочих методов и компьютерных интерпретационных технологий создано в период формирования "парадигмы ранней компьютерной эпохи (условно 1960 - 1985 гг.)" [182] и не являются полностью адекватными геофизической практике. С одной стороны, это связано с идеализированными теоретическими представлениями о геофизических полях и аномалиеобразующих объектах, положенными в основу ряда методов интерпретации [189]; с другой - ограниченными возможностями вычислительной техники указанного периода времени («больших ЭВМ»), не позволяющими реализовать решение многих интерпретационных задач без такого рода идеализации. Вычислительные возможности современных компьютеров позволяют разрабатывать необходимые для производства технологии интерпретации потенциальных полей, измеренных в пределах ограниченной по размерам площади, на криволинейной границе раздела Земля-воздух, в точках нерегулярной сети и осложненных помехами различного происхождения; при количестве наблюденных значений, достигающем десятков - сотен тысяч, а иногда и более.

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

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

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

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

• разработка усовершенствованной технологии вычисления топопоправок при гравиметрической съемке, позволяющей использовать цифровые модели местности большой размерности (до 10Л высотных отметок) и полностью автоматизировать весь процесс вычислений (включая определение поправок для центральной зоны);

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

• разработка приближенного способа расчета топопоправок при аэромагнит

НОЙ съемке, не требующего априорных сведений о петромагнитных характеристиках слагающих рельеф горных пород и высоте съемочных полетов;

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

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

• исследование возможностей разработанных алгоритмов истокообразной аппроксимации для различных преобразований геофизических полей, заданных на криволинейной земной поверхности;

• разработка комплексов программ для зЛета влияния магнитного и гравити-рующего рельефа, объединяющих в себе расчет топопоправок и аналитические аппроксимации полей, а также методических рекомендаций по их практическому применению;

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

• разработка типовой схемы процесса интерпретации данных гравиметрической и аэромагнитной съемок при поисках медно-никелево-платиновых руд в Норильском районе и ее практическая реализация;

• практическая реализация компьютерных технологий интерпретации геопотенциальных полей, зафиксированных в условиях горной местности, при локальном прогнозировании золотого оруденения в Ольховско-Чибижекском районе (Восточный Саян), изучении флангов Талнахского рудного узла (Норильский район), поисках медно-никелевых руд на Дербинской площади (Восточный Саян) и Боотанкагской площади (Центральный Таймыр).

Научная новизна работы заключается в разработке следующих положений:

• сформулирована и решена задача определения точности учета поправок за влияние рельефа местности (топографических масс) при гравиметрии и магнитометрии путем стохастического моделирования, позволяющая осуществлять учет влияния (в т.ч. и совокупного) различных мешающих факторов (латеральной изменчивости петрофизических характеристик верхней части геологического разреза, отклонений в плановой и высотной привязке точек наблюдений, колебаний высот съемочных полетов при аэромагнитной съемке и т.д.);

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

• экспериментально доказана возможность аппроксимации квадратичной регрессионной зависимостью взаимосвязи между высотными отметками рельефа местности и аномальным магнитным полем в условиях субвертикально намагниченной верхней части геологического разреза, что положено в основу алгоритма приближенного определения поправок за влияние магнитного рельефа для аэромагнитных съемок при отсутствии априорной информации о петромагнитных параметрах и высотах съемочных полетов;

• разработаны высокотехнологичные адаптивные алгоритмы построения эквивалентных моделей геологической среды, позволяющие осуществлять аппроксимацию гравитационного и магнитного полей, заданных в узлах равномерной и неравномерной сети, позволяющие осуществлять решение практических задач, размерность которых достигает пхЮ"* значений поля (1 < п < 10) при сравнительно небольших затратах машинного времени;

• на модельных и практических примерах доказаны преимущества методов интерполяции и трансформации потенциальных полей, основанных на построении эквивапентных сеточных моделей геологической среды;

• сформулировано понятие «эффекта разновысотности» - произведения вертикального градиента поля на высотные отметки точек измерений и обоснована целесообразность приведения аномальных геофизических полей к горизонтальной плоскости, что устраняет ложные экстремумы аномалий, расположенные вне проекции искомых объектов на дневную поверхность;

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

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

• создана типовая схема процесса интерпретации геопотенциальных полей Норильского района, основанная на комплексном многовариантном формализованном анализе геопотенциальных полей, в значительной степени очищенных от влияния известных геологических объектов, регионального фона, помех; приведенных к единой горизонтальной плоскости («поверхности наблюдений») и преобразованных в тождественные функции;

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

Практическая значимость работы определяется ее направленностью на решение важных прикладных задач. Разработанное программно-алгоритмическое обеспечение (комплексы программ КЕЬОКУ и КЕЬМАО) характеризуется высокими технико-экономическими показателями, прошло широкую практическую апробацию и используется в ряде производственных и научных организаций. Применение созданных компьютерных технологий интерпретации потенциальных полей позволяет: с высокой точностью осуществлять вычисление поправок за влияние рельефа местности при гравиметрической и магнитной съемках; проводить интерполяцию и трансформацию полей в ЗВ-варианте (с учетом высот поверхности наблюдений) с одновременным подавлением аномалий помех, нарушающих гармонический характер поля; извлекать принципиально новую геологическую информацию из данных полевых наблюдений за счет сочетания в единой технологической цепочке расчета топопоправок и аналитической аппроксимации полей. Достоверность интерпретации данных аэромагнитных съемок, проводящихся в пределах складчатых областей с обтеканием рельефа местности, повышается за счет приведения аномалий магнитного поля к горизонтальной плоскости. Повышение геологической информативности материалов гравиразведки и магниторазведки при решении геокартировочных и прогнозно-поисковых задач обеспечивает сочетание в процессе комплексной интерпретации алгоритмов истокообразной аппроксимации полей и распознавания образов (в случае необходимости дополняющееся решением прямых и обратных задач). Практическая реализация указанного сочетания методов в типовой схеме процесса интерпретации гравитационного и магнитного полей в Норильском районе доказывает принципиальную возможность выявления глубокозалегающих крупных месторождений медно-никелево-платиновых руд по геофизическим данным, что позволяет резко сократить объем дорогостоящего поискового бурения.

Использование представленных в работе компьютерных технологий в практике геофизических исследований, проводящихся в различных регионах России, позволило подготовить геофизическую основу для последующих стадий геологоразведочных работ; локализовать участки, перспективные на обнаружение рудных полезных ископаемых; определить точки заложения буровых скважин и в ряде случаев - выявить при поисковых работах конкретные рудные объекты. В частности, в Норильском районе скважинами Б-12, -13, заданными в пределах выявленной аномальной зоны, вскрыто промышленное никелевое оруденение.

Апробация и публикации. Основные положения и результаты работы докладывались на семинаре "Геологическая интерпретация магнитных и гравитационных аномалий" (Москва, 1992 г.); на семинаре "Компьютерные технологии прогнозирования при геологосъемочных работах масштабов 1:200 CCD -1:50 ООО" (Красноярск, 1992 г.); на научных чтениях, посвященных 100-летию со дня рождения H.H. Урванцева (Норильск, 1993 г.); на семинаре «Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей» (Екатеринбург, 1999); на научно-технической конференции «Уральской высшей школе разведочной геофизики 70 лет» (Екатеринбург, 1999 г.); на региональном совещании «Направление геологоразведочных работ на никель, медь и платиноиды в Таймырском автономном округе в 2000 - 2005 гг.» (Норильск, 1999 г); на научно-практической конференции, посвященной 80-летию геологической службы Западной Сибири «Геологическое строение и полезные ископаемые западной части Алтае-Саянской складчатой области» (Кемерово, 1999 г.); на Таймырской региональной конференции «Минерально-сырьевая база и экология Таймырского автономного округа: состояние и перспективы» (Норильск, 2000 г.); на региональной конференции геологов Сибири, Дальнего Востока и Северо-Востока России «III века горногеологической службе России»; на заседании Научно-методического Совета по геолого-геофизическим технологиям поисков и разведки твердых полезных ископаемых МНР России (Санкт-Петербург, 29.09.2000 г.); на международной Геофизической конференции, посвященной 300-летию горно-геологической службы России (Санкт-Петербург, 2000 г.); на международной школе-семинаре «Вопросы теории и практики комплексной геологической интерпретации гравитационных, магнитных и электрических полей» (Ухта, 2000 г); на международной научно-технической конференции «Горно-геологическое образование в Сибири. 100 лет на службе науки и производства» (Томск, 2001 г.); на общероссийском семинаре «Платина в геологических формациях Сибири» (Красноярск, 2001 г.); на Второй Всероссийской конференции «Геофизика и математика» (Пермь, 2001 г.), на заседаниях геолого-геофизических секций научно-технических советов Норильской КГРЭ, Центрально-Арктической

ГРЭ, ЕГО "Красноярскгеология", Норильского горно-металлургического комбината; на заседаниях научно-технических советов Южной геофизической экспедиции и ОАО «Красноярской горно-геологической компании» и Ученом Совете Красноярского научно-исследовательского института геологии и минерального сырья (КНИИГиМС).

По теме диссертации опубликовано 35 печатных работ, написано 5 отчетов по опытно-методическим и тематическим работам и целый ряд геологических отчетов по производственным работам.

Защищаемые положения:

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

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

3. Использование интерпретационных компьютерных технологий, объединяющих в себе различные методы формализованного анализа геопотенциальных полей, приведенных к единой горизонтальной плоскости («поверхности наблюдений») и в значительной степени очищенных от влияния изученной части геологического разреза, регионального фона, помех негармонического характера повышает геологическую эффективность геофизических исследова

НИИ в сложных физико-геологических условиях, характерных для большинства рудных районов.

Фактический материал и личный вклад автора. Основу диссертационной работы составляют результаты исследований автора, выполненных в период с 1987 по 2001 гг. в Норильской КГРЭ, Таймыргеолкоме и Южной геофизической экспедиции. В том числе по темам: «Разработка и внедрение автоматизированной прогнозно-поисковой системы "Норильск-Поиск"» (гос. per. № 16-88-5/23, 1990 г.); «Адаптация автоматизированных систем обработки и интерпретации геолого-геофизических данных применительно к условиям Норильского района» (гос. per. № 16-90-14/11, 1993 г.); «Разработка программно - алгоритмического обеспечения для учета влияния магнитного рельефа в Норильском промышленном районе» (гос. per. № 56-97-12/7, 1998 г.); «Разработка программно - алгоритмического обеспечения для обработки и интерпретации данных гравиразведки в горной местности» (гос. per. № 56-98-18/6, 1999 г.); «Оптимизация комплекса наземных геофизических наблюдений для целей поисков рудного золота на территории республики Хакасия» (гос. per. № 58-98-2/3, 2000 г.).

Под руководством и при непосредственном участии автора на вычислительном центре Норильской комплексной геологоразведочной экспедиции внедрялось и адаптировалось к условиям региона программно-алгоритмическое обеспечение для обработки и интерпретации потенциальных геофизических полей, разработанное в различных научных организациях бывшего СССР. В результате этой работы была выполнена интерпретация геолого-геофизических материалов, обеспечившая создание надежной крупномасштабной геофизической основы для геологического картирования и поисков медно-никелево-платиновых руд на наиболее перспективных плош;адях Норильского района - Дьянгинской, Тальми-Кумгинской, Хараелахско-Далдыканской и Вологочанской, а также проведено решение разнообразных геологических задач в пределах Северо-Сибирской никеленосной провинции. Фактическим материалом для интерпретационных построений являлись результаты геологоразведочных работ, выполненных в различные годы Норильской КГРЭ, Полярной ГРЭ, Центрально-Арктической ГРЭ, Северо-Енисейской ГФЭ, АО «Корякгеолдобыча» и Южной ГФЭ. Результативная информация представлена в главах и разделах производственных отчетов по геофизическим, геологосъемочным и поисковым работам.

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

Большинство методологических, теоретических и технических идей, представленных в диссертации, базируется на научных работах В.И. Аронова, П.И. Балка, Е.Г. Булаха, Г.С. Вахромеева, В.В. Ломтадзе, A.A. Никитина, Г.Г. Ремпеля, В.Н. Страхова. В процессе проведения исследований, а также при изложении полученных результатов в работе автор опирался на известные элементы и положения из теории потенциала, математического анализа, математической статистики, вычислительной математики, методологии интерпретации геофизических полей, петрофизики, структурной геологии, геологии рудных месторождений, математического программирования и программирования на ЭВМ.

Структура и объем диссертации. Диссертация состоит из введения, четырех глав, заключения и списка литературы из 223 наименований. Она содержит 351 странице основного текста, в том числе 93 рисунка и 32 таблицы.

Заключение диссертации по теме «Геофизика, геофизические методы поисков полезных ископаемых», Долгаль, Александр Сергеевич

ЗАКЛЮЧЕНИЕ

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

Основные научные и практические результаты выполненных исследований сводятся к следующему:

1. Создана методика оценки точности учета влияния рельефа местности, базирующаяся на стохастическом моделировании процесса определения топопо-правок, использующая метод Монте-Карло. Получены реальные оценки погрешностей вычисления поправок за рельеф, возникающих за счет латеральной изменчивости физических характеристик верхней части геологического разреза и колебаний в плановом и высотном положении точек измерений поля. В частности, для Норильского района установлено, что величины этих погрешностей могут быть сопоставимы или даже превышать по амплитуде аномальные эффекты от искомых геологических объектов.

2. Усовершенствована технология вычисления поправок за влияние рельефа местности при гравиметрической съемке, ориентированная на использование ЦММ большой размерности (» 10'' и более высотных отметок), формируемых путем векторизации скан-образов топографических карт. Впервые автоматизирован учет влияния рельефа в пределах т.н. центральной зоны - сравнительно малой области, охватывающей пункт фавиметрических наблюдений. Повышенная точность определения топопоправок в пунктах гравиметрических наблюдений достигается за счет применения ЗВ-интерполяции. Созданное программное обеспечение характеризуется высокими технико-экономическими показателями и используется в ряде производственных организаций.

3. Разработаны алгоритмы и программы определения поправок за влияние рельефа местности при наземной и аэромагнитной съемках. Сведения о пространственном распределении намагниченности горных пород, слагающих рельеф, могут задаваться на основе априорной информации, либо могут быть получены путем пошагового построения регрессионной зависимости между наблюденным магнитным полем и полем, обусловленным ЦММ. Модельные и практические примеры свидетельствую о принципиальной возможности достоверного определения аномальных эффектов от резкорасчлененного магнитного рельефа дневной поверхности (топопоправок). Во многих случаях исключение топопоправок из наблюденного магнитного поля позволяет выделять слабоконтрастные аномалии от геологических тел, которые крайне сложно выявить в исходном поле. Предложен приближенный способ определения топо-поправок при аэромагнитной съемке, основанный на построении квадратичной рефессионной зависимости между высотными отметками рельефа местности и логарифмами амплитуды аномального магнитного поля.

4. Развит подход В.И. Аронова к преобразованиям аномалий потенциальных полей, зафиксированным на криволинейной поверхности наблюдений. Разработаны алгоритмы и профаммы истокообразной аппроксимации геопотенциальных полей, заданных в узлах равномерной и неравномерной сети ( размерность задач пхЮ'А значений поля, 1 < и < 10 и более). Посфоение аппрок-симационной конструкции сводится к решению обратной задачи в линейной постановке, при этом устойчивость решения СЛАУ большой размерности обеспечивается за счет геометрии эквивалентных источников, а повышенная скорость вычислений достигается за счет специальных адаптивных процедур. На модельных и практических примерах показана возможность решения на основе истокообразной аппроксимации широкого круга важных практических задач, например:

• Восстановления значений потенциального поля в узлах регулярной сети с учетом различий в высотных отметках исходных и результативных точек (интерполяция в ЗВ-варианте).

• Исключения искажающего влияния аномального вертикального градиента («эффекта разновысотности») при пересчете наблюденного поля на горизонтальную плоскость. Следует отметить, что при этом устраняется часто отмечающееся смещение эпицентров аномалий от аномалиеобразующих тел, которое может привести к пропуску искомых объектов при проведении заверочных горно-буровых работ.

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

• Фильтрации помех, нарушающих гармонический характер наблюденного поля.

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

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

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

6. Разработаны комплексы программ КЕЬОКУ и КЕЬМАО, предназначенные, соответственно, для интерпретации данных гравиметрической и магнитной съемок в условиях горной местности. Основными методами выделения полезного сигнала при наличии аномалий-помех «рельефной» природы, случайных помех и регионального фона являются вычисление топопоправок (которому при магниторазведке обычно предшествует определение петромагнит-ных характеристик среды) и аппроксимационные преобразования полей с использованием эквивалентных источников. Реализованы различные графы обработки практических геофизических материалов с помош;ью созданного программного обеспечения и проведено решение геологических задач картиро-вочного и поискового характера.

7. Практически реализованы различные схемы проведения комплексной и методной интерпретации геопотенциальных полей, зафиксированных в горной местности:

• Подготовлена геофизическая основа для геологического картирования и поисков медно-никелево-платинового оруденения в сложных физико-геологических условиях Норильского района путем комплексного многовариантного формализованного анализа геопотенциальных полей, приведенных к единой горизонтальной плоскости («поверхности наблюдений»), в значительной степени очиш;енных от влияния известных геологических объектов, регионального фона и помех. Полученные материалы по центральной части территории свидетельствуют о принципиальной возможности выявления крупных месторождений медно-никелево-платиновых руд, залегающих под мощным покровом туфолавовых образований, по данным гравиметрической и аэромагнитной съемок при использовании указанной методики интерпретации. Скважинами Б-12, Б-13, пробуренными в пределах геофизической аномальной зоны, вскрыто промышленное никелевое оруденение.

• Проведено локальное прогнозирование золотого оруденения в Ольховско-Чибижекском районе (Восточный Саян) путем применения различных методов распознавания образов к пересчитанным на горизонтальную плоскость гравитационному и магнитному полям и их трансформантам, также полученным с помощью истокообразной аппроксимации. Определены формальные геофизические критерии прогнозирования (ГФКП) рудного золо

327 та, выделены перспективные участки для проведения дальнейших поисковых работ.

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

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

• Выполнено решение обратной задачи гравиметрии монтажным методом по данным профильных наблюдений на Боотанкагской площади (Центральный Таймыр); в результате четыре из пяти рекомендованных скважин вскрыли дифференцированная интрузия основного состава в заданном интервале глубин.

Список литературы диссертационного исследования доктор физико-математических наук Долгаль, Александр Сергеевич, 2002 год

1. Аведисян В.И. Пути повышения геолого-экономической эффективности гравиметрических съемок. //Разведка и охрана недр. 1991. № 5. С. 29-34.

2. Алексидзе М.А. Приближенные методы решения прямых и обратных задач гравиметрии. М.: Наука. Гл. ред. физ. мат. лит., 1987. 336 с.

3. Анискович М.Ф., Бойко А.З. О возможности автоматической безэталонной классификации гравитационных и магнитных полей в условиях УЩ. // Геофизический журнал. 1989. № 6. С.61-70.

4. Аронов В.И., Бородатый И.И., Фильштинский Л.Е. Опыт вычисления поправок за рельеф местности в горной области при помощи электронных счетных машин //Геофизическая разведка. Вып. 15. М.: Недра, 1964. С. 104-110.

5. Аронов В.И. Обработка на ЭВМ значений аномалий силы тяжести при произвольном рельефе поверхности наблюдений. М.: Недра, 1976. 131 с.

6. Аронов В.И. Методы математической обработки геологических данных на ЭВМ. М.: Недра, 1977. 168 с.

7. Аронов В.И. Методы построения карт геолого-геофизических признаков и геометризации залежей нефти и газа на ЭВМ. М.: Недра, 1990. 300 с.

8. Аронов В.И. Трехмерная аппроксимация как проблема обработки, моделирования и интерпретации геофизических и геологических данных. // Геофизика. 2000 г. №4.0.21-25 .

9. Афанасьев Б.М., Грищенко В. А., Долгаль A.C., Мухаметшин A.M., Панов В.В., Садур О.Г. Геофизический мониторинг при эксплуатации железорудных месторождений. / Геофизика. 2000 г. № 2. С. 55 61.

10. Ю.Балк П.И., Балк Т.В. Детерминистские модели интерпретации для анализа разрешающей способности гравиметрического метода в линейной постановке обратной задачи. // Изв. АН СССР. Сер. Физика Земли. 1987. № 5 С. 41-54.

11. Балк П.И., Долгаль A.C., Балк T.B. Сеточные методы решения обратных задач и опыт их применения при прослеживании дифференцированных интрузий по данным гравиразведки.//Геология и геофизика. 1993. № 5 . С. 127-134.

12. Балк П.И. Использование априорной информации о топологических особенностях источников поля при решении обратной задачи гравиметрии. // Докл. АН СССР. 1989. Т.309. № 5 . С.1082-1084.

13. Балк П.И., Долгаль A.C., Балк Т.В. Сеточные модели плотностной среды и опыт их применения при прослеживании дифференцированных интрузий по данным гравиразведки. // Геология и геофизика. 1993 . № 5. С. 127-134.

14. Балк П.И. «Столкновение геофизических и математических интересов -главный источник противоречий в современной теории интерпретации потенциальных полей».// Геофизический журнал. 2000. Т 22. № 4. С .3 20.

15. Бережная Т.Л., Любимов A.A., Телепин М.А. Пути решения проблемы автоматизации хранения, обработки и интерпретации гравиметрических данных. // Разведочная геофизика. Вып. 104. М.: Недра. 1986. С. 97-104.

16. Бережная Т. Л., Телепин М.А. Определение плотности промежуточного слоя по значениям поля силы тяжести и отметкам рельефа. // Прикладная геофизика. Вып. 104. М.: Недра. 1982. С. 127-134.

17. Блох Ю.И. Комплексирование методов интерпретации гравитационных аномалий при определении формы геологических объектов.//Разведочная геофизика. -М.: Недра, 1984. -Вып. 97- С. 50-54.

18. Блох Ю.И. Решение прямых задач гравиразведки и магниторазведки. М.: Моск. гос. геол.-развед. акад. 1993. 79 с.

19. Блох Ю.И. Решение прямой задачи магниторазведки для трехмерных аниза-тропных геологических объектов с учетом размагничивания // Изв. АН СССР. Физика Земли. 1987. № 12. С. 49-55.

20. Блох Ю.И. Возможности интерпретации магнитных аномалий с учетом размагничивания.// Изв. АН СССР. Физика Земли. 1987. № 4. С. 56-62.

21. Блох Ю.И. Комплексное моделирование сильномагнитных объектов. // Геофизика. 1997. №2. С. 60-63.

22. Богданов Л.А., Забелин В.Г., Петрова A.A., Яновская Ю.А. Объемное моделирование для локального прогноза хромитов. // Разведка и охрана недр. 1993. №5. С. 6-8.

23. Боровко H.H. Количественный анализ поисковых критериев крупных эндогенных рудных месторождений. / Обзор. Сер. геол. методы поисков и разведки месторождений метал, полезных ископаемых. М.: изд. ВИЭМС. 1973. 53 с.

24. Боровко H.H. Оптимизация геофизических исследований при поисках рудных месторождений. Л: Недра. 1979. 230 с.

25. Бродовой В.В. Геофизические исследования в рудных провинциях. М.: Недра. 1984. 269 с.

26. Бугров Н.С., Никольский СМ. Высшая математика. Дифференциальное и интегральное исчисление. М.: Наука. Гл. ред. физ.-мат. лит. 1984. 432 с.

27. Булага В.Х., Ксенофонтов В.А. Интерпретация гравитационных аномалий Припятской впадины методом геологического редуцирования. // Разведочная геофизика. Вып. 100. М.: Недра. 1985. С. 85-89.

28. Булах Е.Г., Ржаницын В.А., Маркова М.Н. Применение метода минимизации для решения задач структурной геологии по данным гравиразведки. Киев: Наук, думка. 1976. 220 с.

29. Булах Е.Г., Зейгельман М.С., Корчагин И.Н. Автоматизированный подбор гравитационных и магнитных аномалий: программно-алгоритмическое обеспечение и методические рекомендации. Деп. в ВИНИТИ № 8363-В86. 1986. 235 с.

30. Булах Е.Г., Маркова М.Н., Бойко П.Д. Математическое обеспечение автоматизированной системы интерпретации гравитационных аномалий. Киев: Наук, думка. 1984. 112 с.

31. Булах Е.Г., Маркова М.Н. Обратные задачи гравиметрии в классе тел, аппроксимируемых прямыми уступами. Программное обеспечение и методические рекомендации. Деп. в УкрИНТЭИ 08.07.92. 1992. 110 с.

32. Булах Е.Г., Маркова М.Н. Подбор геологических моделей по гравитационному полю, осложненному региональным влиянием. // Геофизический журнал. 1990.T.12.N1. С. 9-17.

33. Булах Е.Г., Маркова М.Н. Решение обратных задач гравиметрии методом подбора. //Геофизический журнал. 1992. Т. 14. N4. С 9-19.

34. Булах Е. Г., Маркова М. Н. Решение обратных задач гравиметрии при неизвестном значении избыточной плотности // Геофизический журнал. 1994. Т. 16. №6. С. 3-11.

35. Булах Е.Г., Михеева Т.Л. Прямые задачи гравиметрии в классе звездных тел. Алгоритмическое и программное обеспечение. Методические рекомендации // Деп. в ГНТБ Украины 04.04.96. №892 -Ук96. 1996. 114 с.

36. Булах Е.Г., Шуман В.Н. Основы векторного анализа и теория поля. Киев: Наук, думка. 1998. 359 с.

37. Булах Е.Г., Шиншин И.В. Прямые и обратные задачи гравиметрии для совокупности локальных объектов и построение аналитической модели исходного поля. // Доклады HAH Украины. 1999. № 1. С. 112 -115.

38. Булах Е.Г. Об одном алгоритме решения обратной задачи гравиметрии по аномальному полю, осложненному фоновым влиянием. // Доклады HAH Украины. 1999. № 2. С. 122 126.

39. Булах Е.Г., Шиншин И.В. Об одном аппроксимационном подходе к решению задач структурной геологиии по данным гравиразведки // Теоретичн! та прикладш проблемы нафтогазово! геолог!, Киев: «Карбон-лтд». 2000 г, Т,1, С, 78 -84.

40. Булах Е.Г. Обратная задача магнитометрии в связи с построением аналитической модели исходного поля. // Доклады HAH Украины. 2000. № 9. с. 115 -119.

41. Булах Е.Г., Маркова М.Н., В.В. Прилуков. Об одной задаче гравиметрии для совокупности локальных возмущающих тел. // Геофизический журнал. 2000. T.22.N6. С. 137-147.

42. Бурсиан В.Р. Теория электромагнитных полей, применяемых в электроразведке. Л: Недра. 1972. 368 с.

43. Варламов A.C. Автолокализация аномалий силы тяжести.// Разведочная геофизика. Вып.ЮЗ. М.: Недра. 1986. С. 104-106

44. Вахромеев Г.С. Основы методологии комплексирования геофизических исследований при поисках рудных месторождений, М.: Недра. 1973, 152с.

45. Вахромеев Г.С, Давыденко А.Ю. Моделирование в разведочной геофизике. М.: Недра. 1987. 192 с.

46. Веселов К.Е. Гравиметрическая съемка. М.: Недра. 1986. 312 с.

47. Временные методические указания по комплексированию геолого-геофизических и геохимических методов при поисках медно-никелевых месторождений в Норильском промышленном районе. Л.: НПО "Рудгеофизика". 1983.88 с.

48. Вычислительная математика и техника в разведочной геофизике: Справочник геофизика / Под ред. В.М. Дмитриева. -2-е изд., перераб. и доп. М.: Недра. 1990. 498 с.

49. Голиздра Г.Я. Комплексная интерпретация геофизических полей при изучении глубинного строения земной коры, М.: Недра. 1988. 212 с.

50. Голиздра Г.Я. Основные методы решения прямой задачи на ЭВМ.//Региональная, разведочная и промысловая геофизика. М.: ВИЭМС. 1977. 98 с.

51. ГОЛОМОЛЗИН В.Е. О связи параметров аномалий композиции информативных геофизических признаков с прогнозными ресурсами полезных ископаемых (на примере урановых месторождений) // Российский геофизический журнал. 2000. №17-18.0.19-25.

52. Гольдшмидт В.И. Оптимизация процесса количественной интерпретации данных гравиразведки. М.: Недра. 1984. 184 с.

53. Гольцман Ф.М., Калинина Т.Б. Априорная неопределенность выбора геолого-геофизических моделей и эффективность интерпре- тации гравитационных и магнитных аномалий. // Геология и геофизика. 1991. № 6. С. 127-134.

54. Гольцман Ф.М., Калинина Т.Б. Статистическая интерпретация магнитных и гравитационных аномалий. Л.: Недра. 1983. 248 с.

55. Гольцман Ф.М., Калинин Д.Ф., Калинина Т.Б. Компьютерная технология М иЪТАЬТ многоальтернативной классификации и прогноза по комплексу признаков // Российский геофизический журнал. 2000. № 17-18. с. 64-70.

56. Гордин В.М. Способы учета влияния рельефа дневной поверхности при высокоточных гравитационных измерениях. Обзор ОНТИ ВИЭМС. сер. IX. М.: ВИЭИС. 1974. 89 с.

57. Демидович В.И., Чехович K.M. Подготовка геофизической основы для поисков и ГСР-50 с общими поисками на Вологочанской и Хараелахско-Далдыканской площадях. Т. 1. Норильск: 1995. 439 с.

58. Долгаль A.C., Маркова М.Н., Якимчук H.A. Решение прямых задач грави-разведки с помощью программируемых микрокалькуляторов. /Деп. ВИНИТИ №5425-В29. 1989. 12 с.

59. Долгаль A.C., Маркова М.Н., Якимчук H.A. Решение обратных задач гра-виразведки с помощью программируемых микрокалькуляторов. / Деп. ВИНИТИ № 5426 -В29. 1989. 12 с.

60. Долгаль A.C., Маркова М.Н., Якимчук H.A. Расчет компонент ÄZ и АГ от двухмерных тел произвольной формы с помощью ПМК «Электроника МК-52». / Деп. ВИНИТИ № 7432 -В29. 1989. 9 с.

61. Долгаль A.C., Кашкаров A.A. Петрофизические исследования при поисках золотого оруденения. Методическая разработка. Свердловск. Изд-во СГИ. 1990. 15 с.

62. Долгаль A.C., Кашкаров A.A. Решение прямой и обратной задачи гравираз-ведки с помощью программируемых микрокалькуляторов. Методическая разработка. Свердловск. Изд-во СГИ. 1990. 19 с.

63. Долгаль A.C. Решение обратной задачи гравиразведки при поисках медно-никелевых руд. // Геофизический журнал. 1993. № 6. С. 83-88.

64. Долгаль A.C. О системном подходе к интерпретации гравитационных и магнитных полей (на примере Хараелахской трапповой мульды) // Геофизический журнал. № 5. 1994. С. 58-64.

65. Долгаль A.C., Булах Е.Г., Русанов Э.Б. Объемное моделирование геологического строения Талнахского рудного узла (Норильский район) по гравитационному полю. // Геофизический журнал.1995. № 5. С. 62-66.

66. Долгаль A.C. Состояние и перспективы региональных геофизических исследований в Таймырском автономном округе./ В сб. «Недра Таймыра». Вып. 1. 1995. С. 164-173.

67. Долгаль A.C. Моделирование погрешностей учета влияния рельефа при гравиметрической съемке.// Известия РАН. Сер. Физика Земли. 1997. № 8. С. 88-93.

68. Долгаль A.C. Оценка точности учета влияния рельефа местности при гравиметрической и магнитной съемках.// Доклады академии наук. 1997. Т. 354. № 3. С. 389-391.

69. Долгаль A.C., Христенко Л.А. Учет влияния рельефа при обработке магнито-разведочных данных. //Геофизика. 1997. № 1. С. 51-57.

70. Долгаль A.C. Аналитические аппроксимации геопотенциальных полей и их практическое применение. / В сб. «Геофизические исследования в Средней Сибири». Красноярск: 1997. С. 214 222.

71. Долгаль A.C. Магнитная съемка в условиях горного рельефа, сложенного базальтовой формацией. / В сб. «Недра Таймыра». Вып. 2. 1997. С. 123 -135.

72. Долгаль A.C. Усовершенствование технологии учета влияния рельефа местности при гравиметрической съемке. // Геофизический журнал. 1998. №2. С. 51-57.

73. Долгаль A.C., Чехович K.M. Комплексная интерпретация геопотенциальных полей при поисках медно-никелево-платинового оруденения (Норильский район). // Геология и геофизика. 1998. Т. 39. № 11. С. 1615 1625.

74. Долгаль A.C. Новые возможности учета гравитационного влияния рельефа центральной зоны. / В сб. « Геология и полезные ископаемые Красноярского края». Красноярск. 1998. С. 291-296.

75. Долгаль A.C. Учет влияния рельефа местности на магнитное поле Земли.//В сб. «Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических аномалий». Екатеринбург. 1999. С. 199-200.

76. Долгаль A.C. Аппроксимация геопотенциальных полей эквивалентными источниками при решении практических задач. // Геофизический журнал. 1999. Т. 21. №4. С. 71-80.

77. Долгаль A.C. Ностроение аналитической модели естественного электрического поля. // В сб. « Геология и полезные ископаемые Красноярского края». Красноярск. 1999. С. 199-205.

78. Долгаль A.C. Истокообразная аппроксимация потенциальных геофизических полей, заданных в узлах нерегулярной сети // Геология и минеральные ресурсы Центральной Сибири. Красноярск, 2000 г. С. 193-198.

79. Долгаль A.C. Новые компьютерные технологии учета влияния рельефа местности при гравиметрической съемке. // Материалы региональной конференции геологов Сибири, Дальнего Востока и Северо-Востока России. Томск. 2000. С. 489-491.

80. Долгаль A.C. Практическое применение аппроксимационных преобразований потенциальных геофизических полей. // 300 лет горно-геологической службе России. Международная геофизическая конференция. Санкт-Петербург. 2000. С. 515-516.

81. Долгаль A.C. Аппроксимация аномалий естественного электрического поля совокупностью эквивалентных источников //Геофизический журнал. 2001.Т 23. №1.С. 66-76.

82. Долгаль A.C. Компьютерные технологии интерпретации геопотенциальных полей при поисках медно-никелево-платинового оруденения // Геофизический журнал. 2001. Т23 .№2. С. 106-112.

83. Долгаль A.C. Комплексная интерпретация материалов гравиметрической и магнитной съемок при поисках медно-никелево-платинового оруденения в Норильском районе.// Тезисы докладов общероссийского семинара 20-21 сентября 2001 г. Красноярск. 2001. С. 51-53.

84. Долгаль A.C. Истокообразная аппроксимация аномалий естественного электрического поля. // Материалы Второй Всероссийской конференции «Геофизика и математика». Пермь, Горный институт УрО РАН. 2001. С. 84-94.

85. Долгаль A.C. Разработка и внедрение автоматизированной прогнозно-поисковой системы "Норильск-Поиск". Отчет по опытно-методическим работам за 1988-1990 гг. Т. 1. Норильск: 1990. 209 с.

86. Долгаль A.C. Адаптация автоматизированных систем обработки и интерпретации геолого-геофизических данных применительно к условиям Норильского района. Отчет по тематическим работам за 1990-1993 гг. Т.1. Норильск: 1993.242 с.

87. Долгаль A.C. Разработка программно алгоритмического обеспечения для учета влияния магнитного рельефа в Норильском промышленном районе. Т1. Красноярск: 1998. 87 с.

88. Долгаль A.C. Разработка программно алгоритмического обеспечения для обработки и интерпретации данных гравиразведки в горной местности. Т1: Красноярск. 1999. 82 с.

89. Долгаль A.C., Михеева М.А. Оптимизация комплекса наземных геофизических наблюдений для целей поисков рудного золота на территории республики Хакасия. Т1: Абакан. 2000. 286 с.

90. Дэвис Дж. С. Статистический анализ данных в геологии. Пер. с англ. В 2 кн. /Пер. В. А. Голубевой. Под ред. Д. А. Родионова. Кн. 2. М.: Недра. 1990. 427 с.

91. Зайцев В.Е. Палетки для учета гравитационного влияния высокогорного рельефа местности, аппроксимируемого наклонной плоскостью. // Разведочная геофизика. Вып.68. М.: Недра. 1975. С. 87-92.

92. Иванов В.К. Учет влияния рельефа при магниторазведке. //Разведка и охрана недр. 1977. № 5. С. 42 47.

93. Инструкция по гравиметрической разведке. М.: Недра. 1975. 88 с.

94. Казане В.И. Проблема физических неоднородностей разреза и подходы к ее решению не севере Средней Сибири. //В сб. Геофизические исследования в Средней Сибири. Красноярск. 1997. С. 119-134.

95. Каленицкий А.И., Смирнов В.П. Методические рекомендации по учету влияния рельефа местности в гравиразведке. Новосибирск, СНИИГГиМС. 1981. 160 с.

96. Керимов И.А. Метод Р-аппроксимаций при решении задач гравиметрии и магнитометрии. // Материалы Второй Всероссийской конференции «Геофизика и математика». Пермь, Горный институт УрО РАН. 2001. С. 133-147.

97. Кобрунов А. И. Теория интерпретации данных гравиметрии для сложно-построенных сред. Киев: УМК ВО. 1989. 100 с.

98. Кобрунов А. И. Теоретические основы решения обратных задач геофизики. Ухта: Ухтин. индустр. ин-т. 1995. 228 с.

99. Коваль Л.А. Вычисление поправок за влияние рельефа в гравиметрии с помощью электронных цифровых машин //Изв. АН Каз. ССР, сер. геолог. 1963. №4(55). С. 37-41.

100. Койфман Л.И., Кореневич К.А. Объемное моделирование при решении задач разведочной геофизики. // Геологическое истолкование потенциальных полей. Киев: Наук, думка. 1983. С. 121-130.

101. Комплексирование методов разведочной геофизики: Справочник геофизика. / Под ред. В.В. Бродового, A.A. Никитина. М.: Недра. 1984. 384 с.

102. Комплексная территориальная научно-производственная программа «Но-рильск-95»: «Укрепить и расширить минерально-сырьевую базу Норильского горно-металлургического комбината им. А.П. Завенягина». Кн.1: Норильск, Ленинград, Москва. 1989. 331 с.

103. Константинов Г.Н., Константинова Л. С. Моделирование в рудной магниторазведке (методическое руководство).Новосибирск: СНИИГиМС. 1973. 113 с.

104. Коцур O.e., Марушко Г.В. Подготовка геофизической основы для поисков медно-никелевых руд на Дербинской площади. Отчет Южной геофизической экспедиции за 1992-1999 гг. Кн. 1. Абакан: 1999. 236 с.

105. Кочнев В.А. Адаптивные методы решения обратных задач геофизики. Учебное пособие. Красноярск. Красноярский гос. университет. 1993. 126 с.

106. Культин Н. Delphi 3. Программирование на Object Pascal. СПб.: BHV -Санкт-Петербург. 1998. 304 с.

107. Ломтадзе В.В., Большедворский Г.Г. Вычисление гравиметрических поправок за рельеф местности с применением полиномиальной аппроксимации

108. Гравитационного влияния параллелепипеда. //Разведочная геофизика. Вып. 94. М.: Недра. 1982. С. 119-127.

109. Ломакин А.Б. Петрофизическое картирование слабоконтрастных сред и прогноз месторождений полезных ископаемых. СПб: Изд-во С-Петер. ун-та. 1998. 144 с.

110. Львовский Л.Н. Статистические методы построения эмпирических формул: Учеб. пособие для вузов. 2-е изд., перераб. и доп. М.: Высш. шк. 1988. 239 с.

111. Любимов Г.А., Любимов A.A. Методика гравимагнитных исследований с использованием ЭВМ. М.: Недра. 1988. 303 с.

112. Магниторазведка: Справочник геофизика. / Под ред. Б.Е. Никитского, Ю.С. Глебовского. 2-е изд., перераб. и доп. М.: Недра. 1990. 470 с.

113. Макаров СВ. Отчет о поисках сульфидных медно-никелевых руд в пределах Боотанкагского массива на Центральном Таймыре. Т.1. Норильск: 1991. 167 с.

114. Маловичко А.К., Костицин В.И., Тарунина О.Л. Детальная гравиразвед-ка на нефть и газ. -2-е изд., перераб. и доп. М.: Недра. 1989. 224 с.

115. Маргулис A.C. Вопросы эквивалентности и единственности в обратных задачах гравиметрии.// Свердловск. Тр. Института геофизики УНЦ АН СССР. 1987. 156 с.

116. Мартышко П.С О решении прямой и обратной задачи магниторазведки.// Геофизический журнал. 1982. Т. 4. № 8. С. 39-49.

117. Марченко В.В. Человеко-машинные методы геологического прогнозирования. М.: Недра. 1988. 292 с.

118. Матусевич A.B. Объемное моделирование геологических структур на ЭВМ. М.: Недра. 1988. 184 с.

119. Марченко А.И., Марченко Л.М. Программирование в среде Turbo Pascal 7.0. Киев, Век+, М.: Бином универсал. 1998. 496 с.

120. Методика и результаты крупномасштабной гравиразведки в горных районах. // Методы разведочной геофизики. Л., НПО «Рудгеофизика». 1980. 90 с.

121. Методические рекомендации по средне- и крупномасштабной магнитной картофафии. Л.: НПО "Рудгеофизика". 1990. 84 с.

122. Методические рекомендации по применению комплекса методов интерпретации фавимагнитных данных с использованием компьютерных технологий. /Под ред. И.Д. Савинского. М.: ТОО «МЦАИ». 1995. 93 с.

123. Методические рекомендации по геофизическому обеспечению геологосъемочных работ масштаба 1:200 ООО. СПб: Министерство природных ресурсов РФ, ВИРГ-Рудгеофизика, 2000. 240 с.

124. Мудров А.Е. Численные методы для ПЭВМ на языках Бейсик, Фортран и Паскаль. Томск: МП «РАСКО». 1991. 272 с.

125. Наугольников В.Б. О введении поправки за рельеф местностив фавимет-ровые наблюдения. //Разведочная геофизика. Вып. 8. М.: Недра. 1965. С. 99-104.

126. Никитин A.A. Статистические методы выделения геофизических аномалий. М.: Недра. 1979. 280 с.

127. Никитин A.A. Теоретические основы обработки геофизической информации. М: Недра. 1986. 342 с.

128. Новоселицкий В.М. К теории определения изменения плотности в горизонтальном пласте по аномалиям силы тяжести. // Изв. АН ССР. Физика Земли. 1965. №5. С. 25-32.

129. Новоселицкий В.М., Губайдулин М.Г., Койфман Л.И. Изучение строения осадочного чехла севера Урало-Поволжья на основе гравитационного моделирования. // Геофизический журнал. 1979. № 2. С. 99-104.

130. Нусипов E.H., Ахметов Е.М. Учет влияния рельефа в магниторазведке. // В сб. Развитие методов обработки и интерпретации геофизической информации. / Казахский политехнический институт. Алма-Ата. 1991. С. 60 70.

131. Орлов В.К., Рокотян Е.В. Выделение геологически содержательного регионального фона // Вопросы теории и результаты применения методов интерпретации и моделирования геофизических полей. Свердловск: УрО АН СССР. 1989. С. 110-117.

132. Пигулевский П.И., Тяпкин O.K. Объемное моделирование Володарского щелочного массива (Приазовский блок Украинского щита) по гравитационному полю// Геофизический журнал. 2001. Т. 23. № 1. С. 102 107.

133. Петрищевский А. М. Опыт аппроксимации сложных геологических сред массивом материальных точек // Геология и геофизика. 1981. № 5. С. 105 115.

134. Принципы и методы прогноза скрытых месторождений меди, никеля и кобальта / М.Б. Бородаевская, A.M. Кривцов, А.П. Лихачев и др. М.: Недра. 1987. 246 с.

135. Программы по математическому обеспечению обработки и интерпретации геолого-геофизических данных на ЭВМ. / В.Н. Яковлев, С.Ф. Безруков, А.К. Вайнберг. Л.: НПО "Рудгеофизика". 1982. 354 с.

136. Развитие гравиметрии и магнитометрии в XX веке: Труды конференции. Москва, 23 25 сентября 1996 г. М.: ОИФЗ РАН. 1997. 234 с.

137. Ревякин П.С., Бродовой В.В., Ревякина Э.А. Высокоточная магниторазведка. М.: Недра, 1986. 172 с.

138. Ремпель Г.Г. О введении поправки за рельеф при интерпретации данных аэромагнитной съемки // Труды СНИИГиМС, вып. 30. Новосибирск: СНИИ-ГиМС. 1964. С 36-47.

139. Ремпель F.F. Актуальные вопросы методики введения поправок, связанных с рельефом местности в данные гравиразведки и магниторазведки. // Изв. АН СССР. Сер. Физика Земли. 1980. № 12. С. 75-89.

140. Ремпель F.F., Паршуков Н.П., Вайвод Е.А. Объемное моделирование ту-фогенно-эффузивной толщи Норильского района по данным аэромагнитной съемки и прогноз медно-никелевого оруденения.// Реология и геофизика. 1990. № 10. С. 87 98.

141. Результаты применения моделирования в рудной геофизике в различных районах Сибири. / Под ред. B.C. Моисеева, F.F. Ремпеля. М.: Недра. 1988. 219 с.

142. Савинский И.Д. Программные системы обработки и интерпретации гравитационных и магнитных данных.//Геофизика. 1995. №1. С. 24 31.

143. Салов В.М., Сухов Л.Р. К проблеме интерпретации аэромагнитных данных для прогноза никеленосности в условиях развития трапповой формации. // Реология и геофизика. 1982. № 10. С. 138-142.

144. Семенов A.C. Электроразведка методом естественного электрического поля. Л: Недра. 1980. 446 с.

145. Старостенко В.И., Козленко В.Г., Оганесян СМ. Трехмерное распределение плотности в коре Днепровского грабена. // Геофизический журнал. 1986. Т.8. №6. С. 3-19.

146. Старостенко В.И. Устойчивые численные методы в задачах гравиметрии. Киев: Наук, думка. 1978. 227 с.

147. Степанова Н.Э. Метод S-аппроксимаций и его использование при решении задач гравиметрии. // Материалы Второй Всероссийской конференции «Геофизика и математика». Пермь, Горный институт УрО РАН. 2001. С. 266267.

148. Стохастические модели в морфоструктурном анализе. М.: Недра. 1985. 152 с.

149. Страхов В.Н., Лапина М.И. Определение интегральных характеристик возмущающих масс аппроксимационным методом в задачах гравиметрии и маг-нитометрии//Изв. АН СССР. Физика Земли. 1975. №4. С. 35 -58.

150. Страхов В.Н., Лапина М.И. Монтажный метод решения обратной задачи гравиметрии. // Докл. АН СССР. 1976. Т. 227. № 2. С. 344-347.

151. Страхов В.Н. Вариационные методы в теории линейных трансформаций гравитационных и магнитных аномалий. // Докл. АН СССР. 1990. № 1. С. 63-67.

152. Страхов В.Н. К теории метода подбора. // Изв. АН СССР, серия геофизическая. 1964. С. 494-509.

153. Страхов В.Н. О методе подбора при решении обратных задач гравиметрии и магнитометрии. // Докл. АН УССР. Серия Б. 1975. № 13. С. 990-995.

154. Страхов В.Н, Основные идеи и методы извлечения информации из данных гравитационных и магнитных наблюдений, // Теория и методика интерпретации гравитационных и магнитных аномалий, М,: Изд, ИФЗ АН СССР, 1979, С. 146-269,

155. Страхов В.Н. Алгоритмы редуцирования и трансформации аномалий силы тяжести, заданных на физической поверхности Земли. //Интерпретация гравитационных и магнитных аномалий. Киев: Наук, думка. 1992. С. 4 81.

156. Страхов В.Н. Основные направления теории и методологии интерпретации геофизических данных на рубеже XXI столетия. I. // Геофизика. 1995. № 3. С. 9-18.

157. Страхов В.Н. Основные направления теории и методологии интерпретации геофизических данных на рубеже XXI столетия. П. // Геофизика. 1995. № 4. С. 10-20.

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

159. Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей. М.: ОИФЗ РАН, 1997. С.38-42.

160. Страхов В.Н. Третья парадигма в теории и практике интерпретации потенциальных полей (гравитационных и магнитных аномалий). Ч. III // Электр, науч.-инф. журн. "Вестник ОГГГГН РАН". М: ОИФЗ РАН. 1998. № 1(3). С. 100-152.

161. Страхов В.Н., Страхов A.B. О регуляризации метода наименьших квадра-тов//Геофизический журнал. 1998. Т.20. №6. С. 18-38.

162. Страхов В.Н. Современное состояние и перспективы развития теории интерпретации гравитационных и магнитных аномалий. //Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей. Воронеж. 1998. С. 4 35.

163. Страхов В.Н. Что делать ? (о развитии гравиметрии и магнитометрии в России в начале XXI века) М.: ОИФЗ РАН, 1998. 24 с.

164. Страхов В.Н. Критический анализ классической теории линейных некорректных задач.// Геофизика, 1999. № 3. С. 3 9.

165. Страхов В.Н., Керимов И. А., Страхов A.B. Линейные аналитические аппроксимации рельефа поверхности Земли //Геофизика и математика: Материалы 1-й Всероссийской конференции. М.: ОИФЗ РАН. 1999. С. 198-212.

166. Страхов В.Н. Страхов A.B. О регуляризации систем линейных алгебраических уравнений, возникающих в линейных задачах гравиметрии и магнитометрии // О некоторых вопросах теории интерпретации потенциальных полей. М.: ОИФЗ РАН. 1999. С. 212-218.

167. Страхов В.Н. Геофизика и математика. Методологические основы математической геофизики. // Геофизика. 2000. № 1. С. 3-18.

168. Страхов В.Н. В чем причины различий в понимании взаимоотношений геофизики и математики ? // Геофизика. 2000. № 3. С. 39 47.

169. Страхов В.Н. Главнейшая задача в развитии теории и практики интерпретации потенциальных полей в начале XXI века разрушение господствующего стереотипа мышления. // Геофизика, 2001. № 1. С. 3 - 18.

170. Страхов В.Н. Смена парадигмы в теории линейных некорректных задач. М.: ОИФЗ РАН. 2001. 48 с.

171. Страхов В.Н., Керимов И. А., Степанова И.Э., Страхов A.B., Гричук Л.В. Новый информационный базис фавиметрии и магнитометрии.// Материалы Второй Всероссийской конференции «Геофизика и математика». Пермь, Горный институт УрО РАН. 2001. С. 274 277.

172. Сухов Л.Г., Городнянский A.A. Взаимосвязь формализованного и содержательного подходов к оценке факторов никеленосности трапповых формаций. // Поиски и прогнозирование сульфидных медно-никелевых месторождений. Л.: ПГО Тевморгеология". 1982. С. 60-78.

173. Тарасов Г.А. Электрическое поле над комплексом вертикально поляризованных роводящих сфер. / Вопр. рудн. геофизики. Вып. 2. 1961. С. 61 -67.

174. Тарунина О.Л. Структурно-картировочные возможности фавиразведки в комплексе геолого-геофизических исследований. Пермь: Изд-во Пермского университета. 1993. 200 с.

175. Тархов А.Г., Бондаренко В.М., Никитин A.A. Комплексирование геофизических методов. Учебник для вузов. М.: Недра. 1982. 295 с.

176. Телло Э.Р. Объектно-ориентированное профаммирование в среде Windows: пер. с англ. М.: Наука-Уайли. 1993. 347 с.

177. Тихонов А.Н., Костомаров Д.П. Вводные лекции по прикладной математике. М.: Наука. 1984. 192 с.

178. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. 3-е изд. М.: Наука. 1986. 288 с.

179. Фаронов В.В. Турбо Паскаль (в 3-х книгах). Практика профам-мирования. Книга 2. М., УИЦ МВТУ "ФЕСТО-ДИДАКТИК", 1993. 256 с.

180. Фаронов В.В. Турбо Паскаль (в 3-х книгах). Практика профам-мирования. М., Книга 3. УИЦ МВТУ "ФЕСТО-ДИДАКТЖ". 1993. 304 с

181. Фаронов В.В. Delphi 4. Учебный курс. М.: Нолидж. 1998. 448 с.

182. Филатов В.Г., Захаров СВ., Жбанков Ю.В. Способы простран-ственной обработки и интерпретации гравитационных и магнитных полей. Обзор ВИ-ЭМС Сер. разведочная геофизика. М.: ВИЭМС. 1991. 84 с.

183. Филатова В.Т. Объемная модель Мончегорского рудного района на основе гравимагнитных данных.//Отечественная геология. 1995, №10. С. 65-72.

184. Хесин Б. Э. Рудная геофизика в горных областях. М.: Недра. 1969. 200 с.

185. Цирульский A.B. О редукции потенциальных геофизических полей на внешнюю плоскость //Изв. АН СССР. Физика Земли. 1975. № 7. С. 43 47.

186. Цирульский A.B., Никонова Ф.И., Федорова Н.В. Метод интерпретации гравитационных и магнитных аномалий с построением эквивалентных семейств решений. Свердловск: Изд. Ин-та геофизики АН СССР. 1980. 135 с.

187. Шадрин Л.М. Отчет о результатах поисковых работ по оценке промышленных перспектив на богатые руды северных флангов Талнахского рудного узла. Т.1. Норильск: 1990. 500 с.

188. Шапиро В.Б., Киршин A.B., Мелькановицкий И.М. Усовершенствование методики вычислений поправок за гравитационное влияние рельефа местности на электронно-вычислительных машинах. // Разведочная геофизика. Вып. 59. М.: Недра. 1970. С. 125-129.

189. Шефер У., Балк Т.В. Монтажный метод решения совмещенной обратной задачи грави- и магнитометрии. // Докл. АН. 1992. Т. 327. № 1. С. 79-83.

190. Шрайбман В.И., Жданов М.С, Витвицкий О.В. Корреляционные методы преобразования и интерпретации геофизических аномалий. М.: Недра. 1977. 137 с.

191. Юровских И.Н., Попова Е.Г. Опытно-методические работы по внедрению программ комплексной интерпретации на ЭВМ материалов объемных аэромагнитных съемок в полях развития эффузивных образований Норильского района. Т.1. Лесосибирск: 1992. 74 с.351

192. Юровских В.Н. Вычисление вектора эффективной намагниченности рельефа траппов по аэромагнитным данным. // В сб. Геофизические исследования в Средней Сибири. Красноярск. 1997. С. 277-281.

193. Bjerhammer. On graviti. // Stocholm: Royal Inst, of Techn., 1968. 130 p.

194. Graber Brunner V., Klinger E., Marson I. An improved solution for the problem of upward continuation of gravity field data in rudder topography. // Boll, geofis. teor. et appl. 1991. 33. № 130 - 131. p. 135 - 144.

195. Зидаров Д.П., Обратна гравиметрична задача в геопроучването и геоде-зията. Изд. На Българската Академия на науките. София, 1984. 287 с.

196. Naidu P.S., Mathew M.P.//Fast reduction oofpotential fields measured over an uneven surface to a plane surface.// IEEE Trans. Geosci. and Remote Sens. 1994. 32. №3 . p. 508-512.

197. Pilkington Mark, Urquhart W. E. S. Reduction of potential field data to a horizontal plane.//Geofizics. 1990. 55. № 5 . p. 549 555.

198. Zhou X., Zhong В., Li X. Gravimetric terrain correction by triangular -element method. // Geophisics. 1990. № 2. C. 232 238.