Бесплатный автореферат и диссертация по наукам о земле на тему

Алгоритмическое и программное обеспечение построения цифровых моделей магнитного поля по архивным данным аэромагнитных съемок

ВАК РФ 25.00.35, Геоинформатика

Автореферат диссертации по теме "Алгоритмическое и программное обеспечение построения цифровых моделей магнитного поля по архивным данным аэромагнитных съемок"

4852876

Середкпн Антон Борисович

АЛГОРИТМИЧЕСКОЕ И ПРОГРАММНОЕ ОБЕСПЕЧЕНИЕ ПОСТРОЕНИЯ ЦИФРОВЫХ МОДЕЛЕЙ МАГНИТНОГО ПОЛЯ ПО АРХИВНЫМ ДАННЫМ АЭРОМАГНИТНЫХ СЪЕМОК

Специальность 25.00.35 - геоинформатика

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата геолого-минералогических наук

1 5 СВн 201,

Иркутск - 2011

4852876

Работа выполнена на кафедре технологий геологической разведки» ГОУ ВПО «Иркутский государственный технический университет»

Научный руководитель:

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

Официальные оппоненты:

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

кандидат геолого-минералогических наук Турутанов Евгений Хрисанфович

Ведущая организация:

Федеральное государственное унитарное предприятие «Всероссийский научно-исследовательский геологический институт имени А.П. Карпинского» (ФГУП «ВСЕГЕИ»), г. Санкт-Петербург

Защита состоится 29 сентября 2011 г. в 15 часов на заседании диссертационного совета Д 212.073.01 при Иркутском государственном техническом университете по адресу: 664074, г.Иркутск, ул. Лермонтова, 83, ауд. Е-301. Тел./факс: 8(3952)405-112, e-mail: dis@istu.edu; Seminsky@istu.edu

С диссертацией можно ознакомиться в научной библиотеке Иркутского государственного технического университета, г. Иркутск, ул. Лермонтова, 83

Отзывы на автореферат в двух экземплярах, заверенные печатью учреждения, просим направлять по указанному адресу ученому секретарю совета Мальцевой Галине Дмитриевне, e-mail: dis@istu.edu; тел. 8(3952)405-348, 89149323049

Автореферат разослан « 0-!э » августа 2011 г.

Ученый секретарь диссертационного совета

кандидат геолого-минералогических наук ^j/Lste-eZP—Мальцева Г.Д.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность

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

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

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

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

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

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

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

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

В теоретическом отношении тема работ соответствует приоритетному направлению, утвержденному Президиумом РАН: «Разработка новых методологий, технологий, технических средств и аналитических методов исследований поверхности и недр Земли, ее гидросферы и атмосферы» (постановление № 233 от 01.07.2003г.).

Цели и задачи

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

Для достижения поставленной цели определены и решены следующие задачи:

• выявление особенностей представления основных видов архивных материалов аэромагнитных исследований;

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

• разработка методики оцифровки, преобразования и обработки аэромагнитных данных;

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

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

Фактический материал и методика исследований

Для проведения работ привлекались материалы аэромагнитных съемок территории Иркутской области, хранящиеся в архивах ФГУНПГП «Иркутскгеофизика», в том числе архивные материалы аэромагнитных съемок территории Рудногорской площади: ленты магнитограмм, полученные по результатам аэромагнитной съемки масштаба 1:50000, проведенной в 1976-77 гг. Ленской геофизической партией (Ростова Л.С., Тельнов А.К.); карта графиков аномального магнитного поля СССР, представленная листами 0-48-XXVI, О-48-ХХУП, масштаб 1:200000, 1967 г.; карта изолиний аномального магнитного поля СССР, представленная листами О-48-ХХ, 0-48-ХХ1, 0-48-XXVI, 0-48-ХХУИ масштаб 1:200000, 1967 г.

Геологическая характеристика территории Рудногорской площади производилась с привлечением геологической карты Рудногорской площади масштаба 1:200000, структурной карты поверхности кристаллического

фундамента южной части Сибирской платформы масштаба 1:1000000, проекта на выполнение работ по объекту «Комплексные региональные геофизические работы на Рудногорской площади» и отчета о результатах работ по объекту «Комплексные региональные геофизические работы на Рудногорской площади за 2007—2010 гг.».

Для решения сформулированных задач привлечены методы пространственного и регрессионного анализа данных. Разработка алгоритмического и программного обеспечения проводилась в среде ускоренной разработки программ (RAD) Borland Delphi 2006 (корпорация Borland). Оценка эффективности работы разработанного алгоритмического и программного обеспечения осуществлялась на основе сравнительного анализа результативных моделей магнитного поля с моделями, полученными по результатам работы алгоритмов, предоставляемых апробированными комплексами обработки и анализа пространственно распределенных данных Oasis montaj (Geosoft Incorporated), Surfer (Golden Software) и «Коскад 3D» (Российский Государственный Геологоразведочный Университет).

Научная новизна работы

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

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

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

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

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

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

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

Практическое значение и реализация

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

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

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

Апробация работы и публикации

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

- Региональная научно-техническая конференция, посвященная 100-летнему юбилею профессора, доктора геолого-минералогических наук М.М. Лаврова (Иркутск, Иркутский государственный технический университет, 2007г.);

- Всероссийская научно-техническая конференция «Геонауки», посвященная 145-летнему юбилею со дня рождения профессора В.А. Обручева и 120-летию геологической деятельности в Сибири (Иркутск, Иркутский государственный технический университет, 2008 г.);

- XXIII Всероссийская молодежная конференция «Строение литосферы и геодинамика» (Иркутск, Институт земной коры СО РАН, 2009 г.);

- Вторая Всероссийская научно-практическая конференция молодых ученых и специалистов «Геология, поиски и комплексная оценка твердых полезных ископаемых», посвященная 130-летию со дня рождения В.В. Аршинова (Москва, Всероссийский научно-исследовательский институт минерального сырья им. Н.М. Федоровского, 2009 г.);

- Всероссийская научно-техническая конференция «Геонауки», посвященная 80-летию факультета геологии, геоинформатики и геоэкологии Иркутского государственного технического университета (Иркутск, Иркутский государственный технический университет, 2010 г.).

По теме диссертации опубликовано 6 работ, из них 2 в изданиях, входящих в перечень ВАК.

Структура и объем работы

Диссертационная работа состоит из введения, трех глав и заключения. Текст работы изложен на 130 страницах, содержит 1 таблицу, 45 рисунков и сопровождается библиографическим списком из 95 наименований.

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

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

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

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

Благодарности

Благодарность за консультации, методическую помощь, конструктивное руководство и ценные замечания автор выражает своему научному руководителю д.ф.-м.н., проф. А.Ю. Давыденко. Особая благодарность за внимание, помощь и поддержку О.М. Середкиной. Признательность за поддержку, понимание и консультации начальнику партии комплексного анализа данных геоинформационного центра ФГУНПГП «Иркутскгеофизика» к.г.-м.н. C.B. Гаченко. Благодарность коллективу партии комплексного анализа данных геоинформационнош центра ФГУНПГП «Иркутскгеофизика» - Л.А. Кольга, Е.В. Хромых, A.A. Макарову - за помощь и поддержку при выполнении работы.

СОДЕРЖАНИЕ ДИССЕРТАЦИИ

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

Глава 1. Современные возможности цифрового моделирования магнитного поля по архивным материалам

В п. 1.1. (Современное состояние проблемы обработки архивных материалов аэромагнитных исследований) дана краткая историческая справка

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

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

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

В п. 1.2. (Существующее алгоритмическое и программное обеспечение обработки архивных материалов аэромагнитных съемок) проведена классификация программных продуктов, предлагаемых рынком, для обеспечения основных этапов обработки архивных материалов. Для каждого этапа рассмотрены ведущие представители программного обеспечения.

В общем случае процесс обработки архивных материалов аэромагнитных съемок можно разделить на три основных этапа:

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

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

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

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

- Программы - векторизаторы.

- Программы - конверторы.

- Программы математической обработки данных.

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

ArcView, Maplnfo, Панорама, AutoCAD Map и программы подготовки картографической информации Easy Trace, MapEDIT.

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

В качестве программных продуктов, применяемых для математической обработки (увязки) данных, наибольшее распространение в России приобрел» специализированные комплексы обработки и анализа пространственно распределенных данных Oasis montaj (Geosoft Incorporated, Канада) и Surfer. (Golden Software, США). Из российских разработок стоит отметить программу «Коскад 3D» (Российский государственный геологоразведочный университет).

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

Предлагаемая Oasis montaj (модуль montaj Geophysics Levelling) технология ввода поправок за систематические ошибки и увязки профильных данных на основе секущих опорных маршрутов хорошо зарекомендовала себя применительно к материалам современных аэромагнитных съемок, произведенных с использованием спутниковой навигации и достаточно высокой плотностью опорных маршрутов. Использование модуля montaj Geophysics Levelling в качестве средства для внутренней увязки архивных данных аэромагнитных исследований прошлого века практически не дает положительных результатов. Объясняется это недостаточной плотностью сети опорных съемочных маршрутов, либо частичным/полным отсутствием данных по этим маршрутам в связи с утерей части архивных материалов, а также низкой точностью плановой и высотной привязки аэромагнитных данных того времени.

Предлагаемые Oasis montaj (модуль montaj MAGMAP Filtering), Surfer и «Коскад 3D» средства внутренней увязки, основывающиеся на фильтрации данных, предполагают применение тех или иных цифровых фильтров. Фильтрация данных влечет за собой потерю части полезного сигнала, чьи частоты аналогичны частотам погрешностей, что приводит к сглаживанию, искажению, либо полному исчезновению локальных аномалий магнитного поля, при этом определенная часть погрешностей, не попавших в полосу фильтра, неизбежно пропускается.

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

Для плановой увязки результатов различных съемок наиболее эффективным является применение инструментов, предоставляемых Oasis montaj и Surfer. Эти инструменты имеют схожий принцип работы, заключающийся в обработке данных области наложения увязываемых фрагментов (область наложения данных различных полевых партий или съемок различных годов) для осуществления плавного перехода от одной регулярной сети к другой.

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

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

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

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

Глава 2. Алгоритмическое и программное обеспечение обработки архивных материалов аэромагнитных исследовании

В п. 2.1. (Обоснование методики обработки данных аэромагнитных съемок) проведен анализ основных типов архивных материалов аэромагнитных исследований. По результатам анализа обоснована универсальная методика

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

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

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

Согласно методике процесс обработки архивных материалов аэромагнитных съемок подразделяется на четыре основных этапа:

1. Оцифровка.

2. Конвертирование.

3. Внутренняя увязка.

4. Внешняя увязка.

Оцифровка производится исходя из особенностей представления того или иного типа архивных материалов.

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

В процессе внутренней увязки производится устранение погрешностей плановой и высотной привязки съемочных маршрутов.

Процесс внешней увязки подразделяется на два основных этапа:

1. Приведение разно высотных съемок к единому уровню.

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

В п. 2.2. (Техническая и программная реализация пакета прикладных программ) обосновывается выбор системы программирования, приводится математическое обеспечение обработки архивных материалов, дается описание основных модулей, разработанного автором, пакета прикладных программ.

Пакет прикладных программ для обработки архивных материалов аэромагнитных полей «Аэромаг» представляет собой самостоятельный программный продукт, разработанный в среде ускоренной разработки программ (RAD) Borland Delphi 2006 корпорации Borland.

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

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

На вход всех модулей подаются данные, записанные в файлы формата DAT. Результаты работы также сохраняются в файлы формата DAT.

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

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

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

В п. 2.3. (Выводы) по результатам второй главы сформулировано следующее защищаемое положение:

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

Глава 3. Тестирование технологии построения цифровых моделей по архивным данным на материалах Рудногорской площади

Предпосылками выбора Рудногорской площади для тестирования технологии построения цифровых моделей по архивным данным послужило следующее:

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

- в соответствии с Государственным контрактом, в рамках комплексных региональных геофизических работ, проводимых ФГУНПГП «Иркутскгеофизика» в 2007—2010 гг. на Рудногорской площади с целью изучения геологического строения и перспектив нефтегазоносности зоны сочленения Присаяно-Енисейской синеклизы и Непско-Ботуобинской антеклизы, выявления структур и зон, перспективных на нефть и газ, осуществлялась переобработка архивных материалов аэромагнитных съемок территории.

В п. 3.1. (Геологическая характеристика) приводится общая геологическая характеристика Рудногорской площади, включая стратиграфию, магматизм и тектонику территории.

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

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

В п. 3.2. (.Характеристика магнитного поля по данным аэросъемок) составлена аэромагнитная изученность Рудногорской площади по карте изолиний аномального магнитного поля СССР масштаба 1:200000 1967 г., проведен анализ особенностей магнитного поля территории.

Территория площади покрыта большим количеством разномасштабных аэромагнитных съемок, масштаб варьируется от 1:1000000 до 1:25000, проведенных в 1948—1976 гг. различными аэрогеофизическими партиями.

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

В магнитном поле территории по полосовым аномалиям картируются Каймоново-Кутская зона региона пьного глубинного разлома, Илимский, Тангуйский региональные глубинные разломы, а также глубинные разломы субширотного и северо-западного простираний, не имеющие названий. По линейным аномалиям более высокого порядка частично прослеживается Литвинцевский разлом.

Выделяется локальные аномалии от трубок взрыва с приуроченными к ним железорудными месторождениями «Рудногорское», «Читорминекое» «Татьянинское» и «Пасмурное».

Разломы, выделенные по результатам анализа карты изолиний аномального магнитного поля СССР, совпадают с разломами, вынесенными на структурную карту поверхности кристаллического фундамента южной части Сибирской платформы (Александров В.К.).

В п. 3.3. {Эффективность алгоритмов преобразования архивных данных аэромагнитных съемок) дано обоснование эффективности использования разработанного алгоритмического и программного обеспечения конвертирования архивных материалов аэромагнитных съемок, представленных картами графиков и лентами магнитограмм со схемами ориентиров в форматы, принятые в современном компьютерном картографировании.

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

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

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

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

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

распределенных данных Oasis montaj, Surfer и «Коскад 3 D». Дано обоснование эффективности использования разработанного алгоритмического и программного обеспечения внутренней и внешней увязок архивных материалов аэромагнитных съемок.

Путем расчета уровня статистической значимости критерия Манна -Уитни (U-test) и коэффициента корреляции (КК) между результативными моделями и государственной картой изолиний аномального магнитного поля СССР либо между результативными и эталонными моделями (в случае увязки разновысотных съемок) произведена численная оценка эффективности работы различных средств внутренней и внешней увязки.

Предлагаемые Oasis montaj, Surfer и «Коскад 3 D» алгоритмы внутренней увязки, как правило, основываются на фильтрации данных и предполагают применение тех или иных цифровых фильтров, недостатки применения которых изложены в гл. 1.

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

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

Для оценки эффективности работы алгоритмов внешней увязки наряду с разработанным программным обеспечением использовались программные продукты Oasis montaj и Surfer.

Процесс увязки разновысотных аэромагнитных съемок в разработанном пакете прикладных программ «Аэромаг» так же, как и в Oasis montaj, базируется на алгоритмах аналитического продолжения поля в верхнее/нижнее полупространство. Таким образом, в силу схожести принципов действия указанных алгоритмов, наблюдается практически полная идентичность результативных моделей магнитного поля, полученных с применением этих программных средств (рис.2.в, г).

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

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

Рис.2. Результаты увяжи разновысотных аэромагнитных съемок:

а - эталонная аэромагнитная съемка (по архивным данным высота - 100 м); б — обрабатываемая магнитная съемка (по архивным данным высота - 300 м); в - результат пересчета в Oasis montaj (высота пересчета - 200 м, U-test - 0.96, КК - 0.99); г - результат пересчета в «Аэромаг» (высота пересчета — 201 м, U-test — 0.97, КК — 0.99)

Рис.1. Результаты внутренней увязки:

а - полосовой фильтр (U-test - 0.088, К.К - 0.8); б - фильтр с косинусным сглаживанием (U-test - 0.645, КК - 0.97); в - направленный фильтр с косинусным сглаживанием (U-test - 0.597, КК - 0.83); г - разработанные алгоритмы (U-test - 0.863, КК - 0.99)

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

Предлагаемые Oasis montaj и Surfer алгоритмы приведения ифференцированных по уровню нормального магнитного поля съемок к бщему уровню имеют схожий принцип работы. Он заключается в обработке анных области наложения увязываемых фрагментов (область наложения энных различных полевых партий или съемок различных годов) для существлешш плавного перехода между ними. Поправка за разность уровней ормального магнитного поля увязываемых съемочных участков вносится в сходные данные путем ввода с клавиатуры.

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

Разность уровней нормального магнитного поля, восстановленная по осу дарственным каталогам, в силу значительной удаленности наблюдательных таиций от исследуемой территории несет в себе погрешность определения. !вод неверных поправок зачастую приводит к искажению магнитного поля в бласти наложения съемок (рис.3.а).

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

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

программного обеспечения, проведена оценка ее геологической ффективности.

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

О 20 40 60 80 100 120 140 160 180 200 220 240 нТ.Ч

Рис.З. Результаты увяжи съемок дифференцированных по уровню нормального

магнитного поля:

а - Oasis montaj (U-test - 0.0018. КК - 0.79); б - разработанные алгоритмы (U-test - 0.89, КК '

-0.98) |

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

Применение разработанного алгоритмического и программного обеспечения обработки архивных материалов аэромагнитных исследований позволило привлечь в качестве исходных данных полевые магнитограммы съемки масштаба 1:50000. Это привело к увеличению детальности полученной 1 в результате обработки карты магнитного поля по сравнению с уже имеющейся картой изолиний аномального магнитного поля СССР масшатаба 1:200000, 1967 г. и, как следствие, повышению эффективности геологической интерпретации.

Так же, как и карту изолиний аномального магнитного поля СССР, карту, полученную по результатам обработки архивных материалов аэромагнитных съемок, характеризуют интенсивные полосовые аномалии магнитного поля, связанные с зонами региональных глубинных разломов: Каймоново-Кутским, Илимским, Тангуйским.

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

Помимо аномалий от вулканических трубок взрыва с приуроченными к ним железорудными месторождениями «Рудногорское», «Читорминское», «Пасмурное» и «Татьянинское», имеющих место на карте изолиний аномального магнитного поля СССР, в полученной карте магнитного поля отражается аномалия от вулканической трубки взрыва железорудного месторождения «Ждановское». Кроме того, структура магнитного поля осложнена менее интенсивными локальными аномалиями, вызванными большим количеством рудогроявлений железа на территории площади.

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

В п.3.6. (Выводы) пэ результатам третьей главы сформулированы следующие защищаемые положения:

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

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

-200 -150 -100 -50

50 100 150 200 250 300 иТл

Рис. 4. Интерпретационная схема магнитного поля Рудногорской площади:

1 - оси региональных глубинных разломов; 2 - оси разломов второго порядка; 3 - зоны региональных глубинных разломов; 4 - вулканические трубки взрыва с приуроченными к ним железорудными месторождениями; 5 - контур Рудногорской площади

'¡1 \.**>* 2 3 "¡4 ГТ- "5

- "|б 7 | 9 |]0 ___|11 12 О 13 Ли14

Рис. 5. Обобщенная струкгурно-тектоническая схема Рудногорской площади:

1 - зоны глубинных разломов первого порядка; 2 - оси глубинных разломов первого порядка; 3 - оси глубинных разломов второго порядка; 4 - оси разломов третьего порядка; 5 - оси разломов четвертого порядка; б - Литвинцевский вал; 7 - оси соляных складок Непского типа; 8 - оси соляных складок Илимо-Катангского типа; 9 - линии сейсморазведочных профилей; 10 - Катангский блок фундамента; 11 - Зыряновский блок фундамента: 12 - Ленско-Непский блок фундамента; 13 - железорудные трубки взрыва; 14 -

контур Рудногорской площади

В заключении сформулированы основные результаты диссертационной работы.

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

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

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

Функциональные возможности пакета прикладных программ «Аэромаг» позволяют решать следующие задачи:

1. Преобразование архивных материалов аэромагнитных съемок в цифровые форматы, применяемые в современной компьютерной картографии.

2. Устранение погрешностей плановой и высотной привязки съемочных маршрутов (внутренняя увязка).

3. Приведение данных различных полевых партий и съемок различных годов к единому уровню (внешняя увязка).

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

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

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

Пакет прикладных программ «Аэромаг» 16 июня 2011г. успешно внедрен автором в производственный процесс геоинформационного центра ФГУНПГП «Иркутскгеофизика» в г. Иркутск.

СПИСОК ПУБЛИКАЦИЙ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ

В реферируемых журналах, рекомендованных перечнем ВАК РФ:

1. Середкин А.Б. Алгоритмическое и программное обеспечение для преобразования архивных данных аэромагнитных полей // Вестник ИрГТУ №1(Э7). - Иркутск : ИрГТУ, 2009г. - С. 38 - 42.

2. Середкчн А.Б. Алгоритмическое и программное обеспечение для увязки аэромагнитных данных // «Геоинформатика» - 2010г.- №4. - С. 49 - 53.

В научных репетируемых изданиях и сборниках трудов:

3. Середкин А.Б. Методика оцифровки и обработки аэромагнитных данных, представленных картами графиков магнитного поля // Геология, поиски и разведка полезных ископаемых и методы геологических исследований: сб. избр. тр. регион, науч.-техн. конф., посвященной 100-летнему юбилею профессора, д-ра г.-м.н. М.М. Лаврова. (Иркутск, 2007г.). - Иркутск : ИрГТУ, 2007. - Вып. №7. - С. 31 - 33.

4. Середкин А.Б. Методика обработки архивных данных аэромагнитных полей. // Геология, поиски и разведка полезных ископаемых и методы геологических исследований: материалы Всерос. науч.-техн. конф. «Геонауки», посвященной 145-летнему юбилею со дня рождения профессора В.А. Обручева и 120-летию геолопиеской деятельности в Сибири (Иркутск, 2008г.). - Иркутск : ИрГТУ, 2008. - Вып. № 8. - С. 44 - 48.

5. Середкин А.Б. Проблемы внутренней и внешней увязки аэрогеофизических данных : материалы XXIII Всерос. молодежной конф. «Строение литосферы и геодинамика» (Иркутск, 2009г.). - Иркутск : ИЗК СО РАН, 2009.-С. 295-296.

6. Середкин А.Б. Оценка эффективности работы алгоритмического и программного обеспечения обработки архивных материалов аэромагнитных исследований // Геология, поиски и разведка полезных ископаемых и методы геологических исследований: материалы Всерос. науч.-техн. конф. «Геонауки», посвященной 80-летию факультета геологии, геоинформатики и геоэкологии (Иркутск, 2010г.). - Иркутск : ИрГТУ, 2010. - Вып. № 10. - С. 46 - 52.

Подписано в печать 16.08.2011. Формат 60 х 90 / 16. Бумага офсетная. Печать трафаретная. Усл. печ. л. 1,75. Тираж 100 экз. Зак. 171. Поз. плана 34н.

Лицензия ИД № 06506 от 26.12.2001 Иркутский государственный технический университет 664074, г. Иркутск, ул. Лермонтова, 83

Содержание диссертации, кандидата геолого-минералогических наук, Середкин, Антон Борисович

ВВЕДЕНИЕ.

Глава

СОВРЕМЕННЫЕ ВОЗМОЖНОСТИ ЦИФРОВОГО МОДЕЛИРОВАНИЯ МАГНИТНОГО ПОЛЯ ПО АРХИВНЫМ МАТЕРИАЛАМ.

1.1. Современное состояние проблемы обработки архивных материалов аэромагнитных исследований.

1.2. Существующее алгоритмическое и программное обеспечение обработки архивных материалов аэромагнитных съемок.

1.2.1. Программы - векторизаторы.

1.2.2. Программы - конверторы.

1.2.3. Программы математической обработки данных.

1.3. Выводы.

Глава

АЛГОРИТМИЧЕСКОЕ И ПРОГРАММНОЕ ОБЕСПЕЧЕНИЕ ОБРАБОТКИ АРХИВНЫХ МАТЕРИАЛОВ АЭРОМАГНИТНЫХ ИССЛЕДОВАНИЙ.

2.1. Обоснование методики обработки данных аэромагнитных съемок.

2.1.1. Типы архивных материалов аэромагнитных съемок.

2.1.2. Методика обработки данных аэромагнитных съемок.

2.2. Техническая и программная реализация пакета прикладных программ.

2.2.1. Система программирования. 2.2.2. Математическое обеспечение обработки архивных материалов.

2.2.3. Преобразование архивных данных аэромагнитных съемок в цифровой формат.

2.2.4. Автоматизация плановой и высотной привязки съемочных маршрутов (внутренняя увязка).

2.2.5. Особенности интеграции аэромагнитных съемок (внешняя увязка).

2.3. Выводы.

Глава З

ТЕСТИРОВАНИЕ ТЕХНОЛОГИИ ПОСТРОЕНИЯ ЦИФРОВЫХ МОДЕЛЕЙ

ПО АРХИВНЫМ ДАННЫМ НА МАТЕРИАЛАХ РУДНОГОРСКОИ ПЛОЩАДИ.

3.1. Геологическая характеристика.

3.1.1. Стратиграфия.

3.1.2. Магматизм.

3.1.3. Тектоника.

3.2. Характеристика магнитного поля по данным аэросъемок.

3.2.1. Аэромагнитная изученность.

3.2.2. Характеристика магнитного поля.

3.3. Эффективность алгоритмов преобразования архивных данных аэромагнитных съемок.

3.4. Математическая эффективность алгоритмов увязки аэромагнитных данных.

3.4.1. Внутренняя увязка аэромагнитных данных.

3.4.2. Внешняя увязка аэромагнитных данных.

3.5. Геологическая эффективность результативных моделей магнитного поля.

3.5.1. Характеристика использованных материалов аэромагнитных съемок территории Рудногорской площади.

3.5.2. Оценка геологической эффективности цифровой модели магнитного поля Рудногорской площади.

3.6. Выводы.

Введение Диссертация по наукам о земле, на тему "Алгоритмическое и программное обеспечение построения цифровых моделей магнитного поля по архивным данным аэромагнитных съемок"

Актуальность

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

В*- связи с этим' большое значение имеет возможность* привлечения архивных материалов, ранее проведенных на изучаемых территориях аэромагнитных, исследований. Подобные материалы имеют значительный' потенциал для решения^ поисковых и картировочных задач благодаря« тому, что этот сравнительно недорогостоящий способ получения информации о магнитном поле применялся в массовом порядке в качестве основы для геологического картирования* масштабов 1:200000, 1: 50000* и поисков полезных ископаемых.

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

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

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

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

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

В: теоретическом: отношении тема работ соответствует приоритетному направлению, утвержденному Президиумом РАН: «Разработка новых методологий; технологий,, технических средств и аналитических методов исследований, поверхности и недр Земли, ее гидросферы, и. атмосферы» (постановление № 233 от 01.07.2003 г.).

Цели и задачи

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

Для достижения поставленной цели определены и решены следующие задачи:

• выявление особенностей представления основных видов архивных материалов аэромагнитных исследований;

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

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

• оценка эффективности применения' методики и разработанного алгоритмического и: программного обеспечения1 на фактических архивных материалах аэромагнитных съемок территории Рудногорской площади., . '

Фактический материал и методика исследований I

Для проведения работ привлекались материалы аэромагнитных съемок территории^ Иркутской« области, - хранящиеся в архивах ФГУНГТГП «Иркутскгеофизика», в том; числе, архивные материалы; аэромагнитных съемок территории Рудногорской - площади: ленты магнитограмм; полученные по . результатам аэромагнитной? съемки; масштаба 1:50000, проведенной в 1976-77 гг. Ленской геофизической партией (Ростова Л.С., Тельнов Л.К.); карта графиков аномального, магнитного) поля СССР, представленная^ листами 0-48-ХХУ1; О-48-ХХУП, масштаб 1:200000, 1967г. [79]; карта изолинию аномального магнитного поля СССР, представленная : листами О-48-ХХ, 0-48-ХХ1, 0-48-ХХУ1, 0-48-ХХУИ масштаб 1:200000, 1967г. [80].

Реологическая характеристика территории Рудногорской площади производилась с привлечением геологической: карты; Рудногорской площади масштаба 1:200000 [72], структурной карты поверхности кристаллического фундамента южной части Сибирской платформы масштаба 1:1000000 [2], проекта на выполнение работ по объекту «Комплексные региональные геофизические работы на Рудногорской площади» [71] и отчета о результатах работ по объекту «Комплексные региональные геофизические работы на Рудногорской площади за 2007 - 2010 гг.» [34].

Для решения сформулированных задач привлечены методы пространственного и регрессионного анализа' данных. Разработка алгоритмического и программного обеспечения проводилась^ в» среде ускоренной разработки программ (RAD) Borland-Delphi- 2006 (корпорация Borland). Оценка эффективности работы разработанного алгоритмического и программного обеспечения осуществлялась на основе сравнительного анализа результативных моделей магнитного поля с моделями, полученными по результатам работы алгоритмов, предоставляемых апробированными комплексами обработки и анализа пространственно распределенных данных Oasis montaj,(Geosoft Incorporated), Surfer (Golden Software) и «Коскад - 3D» (Российский Государственный Геологоразведочный Университет).

Научная новизна работы

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

В' основу технологии* положены новые алгоритмические решения для компьютерного картографирования, при построении цифровых моделей магнитного поля по архивным материалам аэромагнитных съемок. К таким решениям относятся:

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

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

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

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

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

Практическое значение и реализация

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

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

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

Апробация работы и публикации

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

- Региональная научно-техническая конференция, посвященная 100-летнему юбилею профессора, доктора геолого-минералогических наук М.М*. Лаврова (Иркутск, Иркутский Государственный Технический Университет, 2007г.);

- Всероссийская научно-техническая; конференция «Геонауки», посвященная 145-летнему юбилею со дня; рождения' профессора В.А. Обручева и 120-летию геологической; деятельности в Сибири (Иркутск, Иркутский Государственный.Технический Университет, 2008г.);

- XXIII Всероссийская молодежная конференция «Строение литосферы и геодинамика» (Иркутск, Институт земной коры.СО РАН, 2009г.);

- Вторая Всероссийская1 научно-практическая; конференция; молодых ученых и специалистов «Геология; поиски» и комплексная оценка- твердых полезных ископаемых», посвященная 130-летию; со: дня рождениям В.В. Аршинова (Москва, Всероссийский! научно-исследовательский: институт минерального сырья им. Н;М; Федоровского, 2009г.);

- Всероссийская научногтехническая конференция «Геонауки», посвященная 80-летию факультета геологии, геоинформатики и- геоэкологии Иркутского: Государственного Технического Университета (Иркутск, Иркутский Государственный:Технический Университет, 2010г.).

По теме-диссертации? опубликовано 6 работ, из них 2 в; изданиях, входящих в перечень ВАК.

Структура и объем работы

Диссертационная работа состоит из введения, трех, глав и заключения. Текст работы изложен на 130 страницах, содержит 1 таблицу, 45 рисунков и сопровождается библиографическим списком из 95 наименований;

Заключение Диссертация по теме "Геоинформатика", Середкин, Антон Борисович

3.6. Выводы

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

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

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

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

Математическая эффективность применения разработанного алгоритмического и программного обеспечения доказана посредством получения наибольших, относительно Oasis montaj, Surfer и «Коскад - 3 D», значений уровня статистической значимости критерия Манна-Уитни и коэффициента корреляции между результативными моделями магнитного поля Рудногорской площади и государственной картой изолиний аномального магнитного поля СССР.

Геологическая эффективность- доказана на примере Рудногорской площади, где на основании анализа данных магнитного поля, полученных посредством обработки архивных материалов аэромагнитных съемок в пакете прикладных программ, «Аэромаг», была составлена интерпретационная схема тектонических нарушений территории. Полученная схема подтверждает результаты совместной интерпретации данных сейсморазведки, гравиметрии, цифровой модели рельефа и дистанционного зондирования Земли [34].

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

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

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

ЗАКЛЮЧЕНИЕ

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

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

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

Функциональные возможности пакета прикладных программ «Аэромаг» позволяют решать следующие задачи:

1. преобразование архивных материалов аэромагнитных съемок в цифровые форматы, применяемые в современной компьютерной картографии; "

2. устранение погрешностей плановой и высотной привязки съемочных маршрутов (внутренняя увязка);

3. приведение данных различных полевых партий и съемок различных годов к единому уровню (внешняя увязка).

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

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

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

Пакет прикладных программ «Аэромаг» 16 июня 2011г. успешно внедрен в производственный процесс геоинформационного центра ФГУНПГП «Иркутскгеофизика» в г. Иркутск.

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

1. Айвазян С.А., Бухштабер В.М., Енюков И.С., Мешалкин Л.Д. Прикладная статистика.: Справочное издание. М.: Финансы и статистика, 1989г. - 607с.

2. Александров. В.К. Структурная карта поверхности кристаллического фундамента южной части Сибирской платформы. Масштаб 1:1000000. 1982г.

3. Альберт А. Регрессия, псевдоинверсия и рекуррентное оценивание. — М.: Наука, 1977г.-224с.

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

5. Берлянд Н.Г., Цирель B.C. Анализ и использование материалов опорной картографической аэромагнитной сети. — М.: Недра, 1972г. -66с.

6. Булина JI.B. О выделении разрывных нарушений на Сибирскойплатформе по данным аэромагнитной съемки // Геология и геофизика

7. Новосибирск: СО РАН 1964г. - С. 105 - 114.t '

8. Введение в Arc View GIS. Рязань: РИНФО, 1997г. - 608с.

9. Вентцель Е.С. Теория вероятностей. -М.: Наука, 1969г. 576с.

10. Воеводин В.В., Кузнецов Ю.А. Матрицы и вычисления. М.: Наука, 1984г.-320с.

11. Временная инструкция по составлению и подготовке к изданию карт аномального магнитного поля масштаба 1:1000000 и 1:200000. М.:V

12. Госгеолтехиздат, 1959г. 28с. П.Выгодский М.Я. Справочник по элементарной математике. - М.:

13. Наука, 1966г. 424с. 12.Выгодский М.Я. Справочник по высшей математике. - М.: Наука, 1977г.-872с.

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

15. Н.Гафаров P.A. Тектоника фундамента и типы магнитных полей Сибирской платформы. // Известия АН СССР. Серия Геология. №7. -М.: АН СССР, 1965г. - С. 95 - 108.

16. Гдадский К.В. Гравиразведка и магниторазведка. М.: Недра, 1967г. -319с.

17. Геоинформационная система «Карта 2008».: Руководство пользователя.^ КБ Панорама, 2009г. - 138с.

18. Геология и полезные ископаемые России в 6-ти томах. Том III. Восточная Сибирь. / Под ред. Малич Н.С. С-Пб.: ВСЕГЕИ, 2002г. -396с.

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

20. Гордин В.М. Очерки по истории геомагнитных измерений. М.: ОИЗФ РАН, 2004г. - С. 63 - 65.

21. Демиденко Е.З. Линейная и нелинейная регрессия. М.: Финансы и статистика, 1981г. - 302с.

22. Демиденко Е.З. Нелинейная регрессия. М.: ИМЭМО, 1984г. - 72с.

23. Дрейпер Н., Смит Г. Прикладной регрессионный анализ. 3- изд. М.: Диалектика-Вильяме, 2007г. - 912с.

24. Дэвис Дж. Статистика и анализ геологических данных. М.: Мир, 1977г.-572с.

25. Зайцев Н.С. О тектоники южной части Сибирской платформы. // Вопросы геологии Азии. Том I. М.: АН СССР, 1954г. - С. 399 - 443.

26. Изаров В.Т. Статистический анализ магнитного и гравитационного полей траппов. // Вопросы рудной геофизики Сибири. Выпуск 92. -Новосибирск: СНИИГГиМС, 1969г.

27. Инструкция по магниторазведке: наземная магнитная съемка, аэромагнитная съемка, гидромагнитная съемка. / Под ред. Глебовского Ю.С., Никитского В.Е. Л.: Недра, 1981г. - 263с.

28. Инструкция по составлению и подготовке к изданию карт аномального магнитного поля (АТа) и масштаба 1:500000 1:25000. - Л.: ВСЕГЕИ, 1977г. - 24с.

29. Клаербоут Дж. Ф. Теоретические основы обработки геофизической информации с приложением к разведке нефти. М.: Недра, 1981г. -316с.

30. Комплекс спектрально-корреляционного анализа'данных «КОСКАД -ЗБ» V. 2008.1.: Руководство пользователя. М.: РГГУ, 2008г. - 119с.

31. Косыгин Ю.А. Тектоника. М.: Недра, 1988г. - 462с.

32. Краснов А.И. Вопросы геологического использования материалов крупномасштабных аэромагнитных съемок. // Аэросъемка и ее применение. Л.: Наука, 1967г.

33. Кремер Н.Ш. Теория вероятностей и математическая статистика.: Учебник. 2-е изд. М.: Юнити, 2004г. - 576с.

34. Кривощеков А.Л., Пастухов Н.П., Руденко Г.Е. Отчет о результатах работ по объекту «Комплексные региональные геофизические работы на Рудногорской площади за 2007 2010 гг.». Том I. - Иркутск: ФГУНПГП «Иркутскгеофизика», 2010г. - 284с.

35. Кулханек О. Введение в цифровую фильтрацию в геофизике. М.: Недра, 1981г.-200с.'

36. Логачев A.A., Захаров В.П. Магниторазведка. Л.: Недра, 1973г. — 352с.

37. Мазмишвили А.И. Способ наименьших квадратов. М.: Недра, 1968г.-437с.

38. Макаров A.A., Гаченко C.B. Опыт анализа потенциальных полей на Сибирской платформе. // Материалы XXIII Всероссийской молодежной конференции «Строение литосферы и геодинамика». Иркутск: ИЗК СО РАН, 2009г. - С. 288 - 289.

39. Макаров A.A., Количественная оценка эффектов, обусловленных увязкой данных аэромагнитных съёмок. // Вестник ИрГТУ №1(41). — Иркутск: ИрГТУ, 2010г. С. 144 - 147.

40. Моисеенко Ф.С. Основы глубинной геологии.: Учебное пособие для вузов. Л.: Недра, 1981г. - 279с.

41. Мостеллер Ф., Тьюки Дж. Анализ данных и регрессия в 2-х выпусках. Выпуск I. М.: Финансы и статистика, 1982г. - 317с.

42. Мэйндоналд Дж. Вычислительные алгоритмы в прикладной статистике. М.: Финансы и статистика, 1988г. - 350с.

43. Наймарк Б.М. , Погребинский Г.А., Резников Е.Л. Практические методы, вычисления преобразования Фурье. // Теоретическая и вычислительная геофизика. М.: Наука, 1974г. - С. 55 - 69.

44. Никитин A.A. Основные принципы обработки геофизических полей с целью выделения аномалий на фоне помех. // Известия вузов, геология и разведка №5. М.: РГГУ, 1976г. - С. 115 - 121.

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

46. Никитин A.A., Трофимова Т.А. Теоретические основы обработки геофизических данных.: Учебное пособие. МГГА (МГРИ). М.: МГГА, 1987г.-99с.

47. Никитский В;Е. Аэромагнитные методы в геофизике. М.: Недра, 1966г.-224с.

48. Нейве A.B. Тектоника и магнетизм; // Известия: АН СССР. Серия -Геология №3v-М: АН СССР, 1961г. С. 36 - 55.

49. Петров A.B., Трусов А.Н. Компьютерная технология статистического и спектрально-корреляционного; анализам трехмерной, информации — КОСКАД 3D.// Геофизика №4 - М.: ЕАГО, 2000г. - С. 29 - 33';

50. Принципы магнитной картографии и методика составления-карт. Л.: ВСЕГЕИ, 1985г. - 327с.

51. Программа автоматизированной векторизации картографических материалов MapEDIT.: Руководство пользователя. — М.: Резидент, 2005г. 100с.,

52. Савинский И.Д. Программные системы обработки и интерпретации гравитационных и магнитных данных. // Геофизика №1. М.: ЕАГО, 1995г.-С. 24-31.

53. Середкин А.Б. Алгоритмическое, и программное обеспечение для преобразования: архивных; данных; аэромагнитных полей. // Вестник ИрГТУ №1(37). Иркутск: ИрГ'ГУ, 2009г. - С. 38 - 42. .

54. Ссредкин А.Б. Проблемы внутренней и внешней; увязки аэрогеофизических данных. // Материалы XXIII Всероссийской молодежной конференции «Строение литосферы и геодинамика». — Иркутск: ИЗК СО РАН, 2009г. С. 295 - 296.

55. Середкин А.Б. Алгоритмическое и программное обеспечение для увязки аэромагнитных данных. // «Геоинформатика» № 4 М: ВНИИГЕОСИСТЕМ, 2010г. - С. 49 - 53.

56. Серкеров С.А. Гравиразведка и магниторазведка. М.: Недра, 2006г. — 437с.

57. Смирнов Н.В., Дунин Барковский И.'В. Курс теории вероятностей и математической статистики для технических приложений. - М.: Наука, 1969г.-511с.

58. Соколов' К.П. Геологическое истолкование магниторазведочных данных. М.: Госгеолтехиздат, 1956г. - 127с.

59. Тафеев Г.П., Соколов К.П. Геологическая интерпретация магнитных аномалий. Л.: Недра, 1981г. - 327с.

60. Техническая инструкция по магнитной разведке. М.: Госгеолтехиздат, 1963г.-247с.

61. Трондин В.И., Пастухов Н.П. Проект на выполнение работ по объекту «Комплексные региональные геофизические работы на Рудногорской площади». Иркутск: ФГУН1И П Иркутскгеофизика, 2007г. - 123с.

62. Трондин В.И., Пещерова Н.С. Геологическая карта Рудногорской площади. Масштаб 1:200000. Листы: 0-48-ХХ, O-48-XXI, 0-48-XXVI, 0-48-XXVII, O-48-XXXII, 0-48-XXXIIL Иркутск: ФГУНПГП Иркутскгеофизика, 2007г.

63. Фаронов В.В. Delphi. Программирование на языке высокого уровня. -С-Пб.: ПИТЕР, 2004г. 640с.

64. Федынский В.В. Разведочная геофизика. М.: Недра, 1967г. - 672с.

65. Фленов М.Е. Библия Delphi. С-Пб.: БХВ-Петербург, 2004г. - 880с.

66. Хемминг Р.В. Цифровые фильтры. М.: Недра, 1987г. - 221с.

67. Цирульский А.В. Вопросы теории и методы интерпретации потенциальных геофизических полей.: Учебное пособие. Л.: Горный институт им. Плеханова, 1989г. - 95с.

68. Черемисина Е.Н., Липилин А.В., Никитин А.А. Компьютерные технологии комплексной интерпретации геолого-геофизических данных: современное состояние и перспективы развития. // Геоинформатика №3. М.: ВНИИГЕОСИСТЕМ, 2000г. - С. 98 - 105.

69. Щупак Л.М., Давыдов В.Ф. Карта графиков аномального магнитного поля СССР. Масштаб 1:200000. Листы: O-48-XXVI, 0-48-XXVII. 1967г.

70. Щупак Л.М., Давыдов В.Ф. Карта изолиний аномального магнитного поля СССР. Масштаб 1:200000. Листы: 0-48-ХХ, 0-48-XXI, 0-48-XXVI, 0-48-XXVII. 1967г.

71. Яновский Б.М. Земной магнетизм в 2-х томах. Том II. Теоретические основы магнитометрического метода исследования Земной коры и геомагнитные измерения. Л.: ЛОЛГУ, 1963г. — 462с.

72. Blackman R.B. Data Smoothing and Processing. Reading. Mass.: Addison-Wesley, 1965.

73. Brigham E.O. The Fast Fourier Transform. Englewood Cliffs, N.J.: Prentice-Hall, 1974.

74. Digital Signal Processing. / Edited by Rabiner L.R., Rader C.M. -New-York: IEEE Press, 1975.

75. Golden Software Surfer.: Руководство пользователя. Часть 1. Golden Software Inc., 1997r. - 198c.

76. Kendall M.G., Stuart A. The Advanced Theory of Statistics v.3. -New-York: Hafher, 1968.

77. Lighthill M.J. Fourier Analysis and Generalized Functions. New-York: Cambrige University Press, 1960.

78. Lilley F.E.M., Hitchman A.P., Lie Jun Wang. Time-varying effects in magnetic mapping: Amphidromes, doldrums and induction hazard. // Geophysics vol. 64 №6, 1999. P. 1720 - 1729.

79. MapInfo Professional v. 9.0.: User Guide. Maplnfo Corporation, 2007. -560p.

80. Montaj Geophysics v. 7.2.: Tutorial and User Guide. Geosoft Inc., 2010. -42p.

81. Montaj MagMap Filtering v. 7.1.: Tutorial. Geosoft Inc., 2009. - 83p.

82. Oasis montaj v. 7.2. Mapping and Processing System.: Quick Start Tutorials.- Geosoft Inc., 2010. 325p. 93.Papoulis A. The Fourier Integral and Its Applications. - New-York:

83. McGroaw-Hill, 1962. 94.Spector A., Grant F.S. Statistical models for interpreting aeromagnetic data.

84. Geophysics №2, 1970. P. 293 - 302. 95.Surfer.: Getting Started Guide. - Golden Software Inc., 2009. - 56p.