CC BY
24
УДК 550.382.4
ВЛИЯНИЕ СФЕРИЧНОСТИ ЗЕМЛИ НА РЕЗУЛЬТАТЫ АЭРОМАГНИТНЫХ ИЗМЕРЕНИЙ НАД ПЛАТОБАЗАЛЬТАМИ (НА ПРИМЕРЕ НОРИЛЬСКОГО РАЙОНА)
Долгаль Александр Сергеевич,
dolgal@mi-perm.ru
Учреждение Российской академии наук «Горный институт» Уральского отделения РАН, Россия, 620007, г. Пермь, ул. Сибирская, 78-А.
Актуальность исследования обусловлена целесообразностью использования модельных представлений, адекватных реальным физико-геологическим условиям, при редуцировании результатов магниторазведки в областях развития платобазальтов. Цель: оценка влияния сферичности Земли на аномальное магнитное поле в условиях резко расчлененного рельефа дневной поверхности, сложенного магматическими горными породами, обладающими повышенной намагниченностью (на практических материалах).
Методы: вычисление поправок за влияние магнитного рельефа местности путем решения прямой задачи магниторазведки для области с радиусом 20 км в двух вариантах: для плоской и шарообразной моделей Земли. Выполнено сопоставление полученных результатов и пространственно-статистический анализ выявленных различий. Предполагалось, что вектор намагничивания слагающих рельеф горных пород направлен по современному геомагнитному полю региона. Переменная по латерали эффективная намагниченность верхней части разреза определялась путем решения линейной обратной задачи магниторазведки. Результаты. Кратко охарактеризована проблема учета влияния рельефа местности при измерениях магнитного поля над вулканогенными образованиями. На примере материалов крупномасштабной высокоточной аэромагнитной съемки, выполненной в центральной части Норильского рудного района, показаны нелинейные искажения сигнала, связанные с криволинейностью земной поверхности. Амплитуда этих искажений (по модулю) достигает 100-200 нТл и более при сравнительно небольших различиях высотных отметок для разных моделей. Это объясняется высокими значениями вертикального градиента магнитного поля. Выводы. С целью подавления интенсивных аномалий-помех, обусловленных магнитными породами, слагающими верхнюю часть геологического разреза, предложено вычисление топографических поправок с использованием сферической модели Земли. Представляется также перспективным использование статистического подхода, в неявном виде учитывающего форму земной поверхности, для выделения аномалий «рельефной» природы.
Ключевые слова:
Платобазальты, рельеф местности, аэромагнитная съемка, аномалия, топографическая поправка, редуцирование, модель Земли.
Введение
По сравнению с гравиразведкой, в магниторазведке гораздо слабее проработаны вопросы, касающиеся учета влияния рельефа земной поверхности, сложенного интенсивно намагниченными горными породами, в частности - платобазальтами. Плато-базальты, образовавшиеся при линейных излияниях лавы в периоды тектоно-магматической активации платформ, развиты в Сибири, Индостане, в Южной Америке, Южной Африке, Тасмании, Антарктиде и занимают около 5 млн км2 площади континентов [1]. В частности, 250000 км2 территории России занимает трапповое плато Путорана, второе (после плато Декан) по величине в мире, абсолютные отметки высот в пределах которого достигают 1500-1700 м. С процессами платформенного магматизма связаны месторождения цветных и редких металлов [2], железных руд, камнецвет-ного сырья и др. полезных ископаемых, для поиска которых используются геофизические методы, в т. ч. аэромагнитная съемка (АМС).
Характерной особенностью платобазальтов является сложный характер аномального магнитного поля ДТ, обусловленный их высокими магнитными свойствами (эффективная намагниченность J8(1)=3-5 А/м и более) и резко расчлененным рельефом местности, практически полностью исключающий «возможность интерпретации аномалий,
связанных с более крупными и глубинными объектами» [3]. Моделирование магнитных аномалий-помех «рельефной» природы осуществляется путем решения прямой задачи магниторазведки [4] для цифровой модели местности (ЦММ) при априорно заданной намагниченности верхней части геологического разреза, как правило, зависящей от координат: J#=J#(x,y) [5]. Полученные таким образом значения магнитного поля являются поправками за влияние рельефа местности 8Тр (топографическими поправками) [5-7]. Физический смысл топографической редукции в гравиразведке и магниторазведке одинаков - приближенное определение аномального эффекта от приповерхностного слоя горных пород с целью его последующего исключения из наблюденных значений поля [8]. Данная редукция может применяться совместно с пересчетом аномального магнитного поля на горизонтальную плоскость [5] или с его аналитическим продолжением на поверхность, конкордант-ную рельефу земной поверхности («draping» в зарубежной терминологии) [9, 10].
При решении прикладных задач в качестве моделей Земли используются: плоскость, шар, эллипсоид вращения, квазигеоид, геоид. Вычисление гравитационных аномалий в редукции Буге может базироваться как на плоском [11], так и на сферическом [12] представлениях о слое однородных по
26
DOI 10.18799/24131830/2019/11/2345
плотности горных пород о=со^ в окрестностях пункта наблюдений. Установлено, что значения соответствующих поправок имеют замет-
ные различия [13]. При вычислении поправок за влияние магнитного рельефа STр традиционно используется модель «плоской Земли», геометрические параметры которой характеризует прямоугольная система координат Zxyz, а нижней границей ЦММ является горизонтальная поверхность. Такая идеализация физико-геологических реалий на первый взгляд кажется вполне оправданной, т. к. область учитываемого влияния рельефа 5 в магниторазведке, как правило, меньше, чем в гра-виразведке (для средне- и мелкомасштабных гравиметрических съемок радиус области 5 равен 200 км). Однако влияние формы Земли вносит существенный вклад в амплитуду магнитных аномалий, отвечающих влиянию выходящих на дневную поверхность вулканогенных пород, и его требуется принимать во внимание при разработке методов редуцирования результатов аэромагнитных измерений.
Методика и результаты обработки данных
аэромагнитной съемки
В данной статье будет рассмотрен вопрос о целесообразности использования «сферической» модели Земли при расчете топографических поправок 8Тр на примере материалов крупномасштабной АМС нового поколения, выполненной НФ ВСЕГЕИ в 2012-2013 гг. в центральной части Норильского рудного района на площади 10000 км2. Комплексная аэрогеофизическая съемка выполнялась на самолете Ан-3 по системе прямолинейных широтных маршрутов с обтеканием рельефа, расстояние между маршрутами 250 м, средняя высота полета Япол - 122 м. Для измерений магнитного поля использовался аэромагнитометр ГТ-МАГ, укомплектованный квантовым цезиевым датчиком С8-2 фирмы SCINTREX (Канада). Среднеквадратиче-ская погрешность аэромагнитной съемки для поля с горизонтальным градиентом менее 50 нТл/км составила ±0,89 нТл. Точность определения координат с использованием GPS-технологий составила ±1,55 м.
Рис. 1. Карта рельефа дневной поверхности центральной части Норильского района: 1 - границы зон развития вулканогенных образований; 2 - рудоносные интрузии; 3 - контур площади аэромагнитной съемки; 4 - профиль АБ; 5 - «палетка» для расчета топографических поправок STp
Fig. 1. Map of the surface relief of the Norilsk region central part: 1 are the boundaries of the zones of development of volcanic formations; 2 are the ore-bearing intrusions; 3 is the contour of aeromagnetic shooting area; 4 is the AB profile; 5 is the «palette» for calculation of topographical corrections STp
В пределах площади съемки находятся Хан-тайско-Рыбнинское поднятие, сложенное преимущественно терригенно-карбонатными, сульфато-ностными и терригенными угленосными отложениями палеозоя, а также брахиосинклинали (мульды), выполненные туфобазальтовыми формациями поздней перми и раннего триаса - Норильская, Вологчанская, Хараелахская, Иманг-динская, Иконская [14]. Зоны развития вулканогенных пород, занимающие около 48 % рассматриваемой территории, характеризуются преимущественно расчлененным рельефом дневной поверхности с абсолютными отметками высот рельефа до 700-1000 м и более. Широко развиты интрузивные образования триаса, в т. ч. дифференцированные тела габброидов, с которыми связаны крупнейшие в мире платино-медно-никелевые месторождения Норильск 1, Талнахское, Октябрьское (рис. 1). Амплитуда аномального магнитного поля в районе исследований изменяется от -1500 до 2150 нТл, отмечается пространственная корреляция интенсивных положительных аномалий ДТ с наиболее приподнятыми участками в пределах трапповых мульд. Следует отметить, что интенсивность локальных магнитных аномалий ДТ, обусловленных дифференцированными интрузиями, как правило, не превышает 50-100 нТл. Выделение аномалий «рудной» природы в пределах трапповых мульд крайне затруднительно и традиционные методы трансформации магнитного поля для этой цели не дают положительных результатов.
При вычислении топографических поправок 8ТР в качестве элементарных тел, аппроксимирующих рельеф земной поверхности, использовались вертикальные прямоугольные призмы. Размер оснований призм 100x100 м отвечал шагу сети ЦММ. Вектор намагниченности каждой призмы был направлен по современному магнитному полю Т0 региона исследований: наклонение !=83,5°, склонение 5=17,6°. Апликаты верхних оснований призм отвечали абсолютным отметкам высот Н земной поверхности, нижние основания располагались на уровне 2=0. Путем векторизации крупномасштабных топографических карт была сформирована ЦММ, имеющая площадь 30500 кв. км, представляющая собой прямоугольную матрицу высот Н, размером 1581 строка, 1931 столбец. Полученное в процессе АМС аномальное магнитное поле было представлено в виде матрицы ДТ размером 911 столбцов, 1229 строк. С использованием данных радиовысотомера была сформирована матрица абсолютных отметок точек измерений Нпол.
С целью оценки размеров области й учитываемого влияния рельефа («палетки»), имеющей форму круга с радиусом Я, применялось вычисление топографических поправок 8ТР при последовательно увеличивающихся значениях Я: 2<Я<25, км. Определялась энергия Е2 магнитных аномалий для каждой из результативных матриц 8Тр. Выход функции/=Е2(Я) на асимптоту послужил критери-
ем для выбора скользящего окна размером Яопт=20 км для дальнейших расчетов.
Для приближенной оценки петромагнитных характеристик верхней части геологического разреза использовались материалы ранее выполненной аэромагнитной съемки масштаба 1:100000 на постоянной барометрической высоте 2400 м. Величина модуля намагниченности Jдф=Jф(x,y) слагающих рельеф горных пород определялась путем решения линейной обратной задачи магниторазведки. Изучаемый объем геологической среды для территории 30500 кв. км (площадь АМС-25 и зона обрамления, необходимая для расчета топографических поправок) был представлен в виде прямоугольных блоков с размерами в плане 1000x1000 м, достигнутая невязка наблюденного и модельного полей в евклидовой метрике составила ~0,7 нТл. Расчетные значения намагниченности вулканогенных образований изменяются от 0,07 до 4,22 А/м, что согласуется с имеющимися петрофизическими данными.
Амплитуда поправок за влияние рельефа изменяется в диапазоне от -620 до 822 нТл, наибольшие положительные значения 8ТР отвечают водораздельным зонам развития молодых базальтов с высокими содержаниями ферромагнитных минералов, над долинами рек отмечаются линейные зоны отрицательных значений 8ТР. В редуцированном (полученном после ввода топографических поправок) магнитном поле ДТред и его трансформантах нашли свое отражение почти все известные медно-никелевые месторождения Норильского района. В результате подавления аномалий-помех, созданных верхней частью разреза, была получена принципиально новая геологическая информация, которую можно использовать для решения прогнозно-поисковых задач.
Рис. 2. Изолинии разности высот «плоской» и «сферической» моделей Земли в пределах области S («палетки»), м
Fig. 2. Isolines of the height difference between «flat» and «spherical» models of the Earth within the area S («palettes»), m
Весь процесс обработки данных, как было сказано выше, выполнялся в рамках модели «плоской» Земли - нижней границей ЦММ являлась горизонтальная плоскость (уровенная поверхность z=0, относительно которой осуществлялся отсчет высот рельефа Я). Более точным приближением к форме нашей планеты является шар с радиусом г=6371 км - «сферическая» модель Земли. Рассмотрим случай, когда поверхности «плоской» и «сферической» моделей Земли соприкасаются в центре О области й «палетки». Разность высот Дz для этих моделей в произвольно выбранной точке с
координатами (х,у) в проекции Гаусса-Крюгера определяется по формуле
Д = гг2 - х2 - у2
и на внешней границе области й составляет ~31,4 м (рис. 2). Можно осуществить изгиб уровенной поверхности Я=0, в вычислительном плане сводящийся к корректировке высот Я в пределах «палетки», в соответствии с параметром Дг, без изменения планового положения точек ЦММ. Таким образом, центры нижних оснований моделирующих рельеф призм теперь будут находиться на поверхности сфе-
Рис. 3. Графики геолого-геофизических параметров по профилю АБ: а) аномальное магнитное поле АТ; б) топографические поправки STp для «плоской» (сплошная линия) и для сферической (точки) Земли; в) разность А топографических поправок для «плоской» и «сферической» Земли; г) рельеф дневной поверхности: 1 - осадочные отложения; 2 - вулканогенные образования
Fig. 3. Graphs of geological and geophysical parameters along the AB profile: а) anomalous magnetic field АТ; б) topographic corrections STp for the «flat» (solid line) and for the spherical (point) of the Earth; в) difference А of topographic corrections for the «flat2 and «spherical» Earth; г) surface topography: 1 - sediment; 2 - volcanogenic formations
рического сегмента высотой ~31,4 м. Этот подход был реализован при втором цикле вычислений поправок за влияние магнитного рельефа, отвечающих «сферической» модели Земли.
Рассмотрим различия в значениях 8Тр для разных моделей Земли на примере профиля Аб, пересекающего южную часть Хараелахской мульды (рис. 1), с перепадом высот рельефа 740 м (рис. 3, г). Аномальное магнитное поле ДТ по профилю изменяется от - 465 до 1296 нТл (рис. 3, а); размах значений топографических поправок составляет ~1400 нТл; графики 8Тр для двух алгоритмов расчета, в представленном на рис. 3, б вертикальном масштабе, очень близки между собой. Однако разность Д топографических поправок для «плоской» и «сферической» Земли изменяется в диапазоне [-115, 115] нТл, при среднем значении т=0,8 нТл и среднеквадратическом отклонении 8=±25,6 нТл (рис. 3, б).
-120 -80 -40 О 40 80 120 160
Разность Л топопоправок, нТл
Рис. 4. Гистограмма разности А поправок 8Tt за влияние рельефа местности в пределах центральной части Норильского района
Fig. 4. Histogram of the difference A of corrections STp for the influence of the terrain in the central part of the Norilsk region
Для всей площади АМС этот диапазон еще шире и составляет [-252, 251] нТл, при m=2,6 нТл и 5=±29,7 нТл (рис. 4). Максимальная величина A может достигать ~16,5 % от наибольшего (по модулю) значения поправки 5ТРСФ для сферической Земли и примерно в 2 раза превышает аномальный эффект от рудоносных интрузий. Среднее значение квадрата амплитуды (энергия) для разности A топографических поправок составляет ~2,3 % от средней энергии 5ТРСФ.
Следовательно, пренебрежение сферичностью Земли в данных физико-геологических условиях приводит к весьма заметным высокочастотным искажениям магнитного поля ДТред. Эти искажения обусловлены высокими значениями вертикального градиента й(8Тр)/йг над платобазальтами, достигающими 5-8 нТл/м и более. Градиент магнитного поля может быть как положительным, так и отрицательным по знаку.
Заключение
Представленные результаты приводят к выводу о необходимости использования представлений о «сферической» (а в дальнейшем - «эллипсоидальной») Земле при разработке компьютерных технологий учета влияния магнитного рельефа местности, базирующихся на решении прямой задачи магниторазведки (в т. ч. с использованием сферических аппроксимирующих тел). Это полностью согласуется с современными тенденциями в редуцировании данных гравиметрических съемок [15] и обеспечивается вычислительными возможностями современных компьютеров [16]. Еще одним перспективным направлением исследований является развитие статистических приемов определения топографических поправок 8Тр, в частности - использующих метод EMD (эмпирической модовой декомпозиции) [17, 18]. Применение статистических методов целесообразно при обработке фондовых материалов аэромагнитных съемок, характеризующихся низкой точностью плановой привязки съемочных маршрутов и отсутствием данных о фактических высотах полетов. В этом случае учет криволинейной формы земной поверхности будет происходить в неявном виде. При выполнении АМС на больших территориях [19, 19] «сферическая» модель Земли может повысить достоверность результатов количественной интерпретации полученных материалов.
Автор выражает глубокую благодарность директору НФ ВСЕГЕИ канд. геол.-минерал. наук ФД. Лазареву и главному геофизику НФ ВСЕГЕИ П.В. Кирплюку, благодаря которым был апробирован в производственных объёмах ввод поправок за влияние магнитного рельефа местности в результаты высокоточных аэромагнитных наблюдений. Полученные результаты включены в отчет по объекту «Опережающие геофизические исследования в пределах Хантайско-Рыбнинского поднятия с целью оценки перспектив медно-никелевого оруденения».
Работа выполнена при поддержке РФФИ (проекты № 19-05-00654А, 18-05-70094 «РесурсыАрктики»).
СПИСОК ЛИТЕРАТУРЫ
1. Арискин А.А. Формации платобазальтов Земли. Океанические плато и континентальные траппы. Лекция № 17. URL: http://wiki.web.rU/images/b/b0/Magmatic_Petrology_Lec11_ 2013.pdf (дата обращения 17.07.2018).
2. Krivolutskaya N.A. Siberian Traps and Pt-Cu-Ni deposits in the Noril'sk area. - Springer, 2016. - 361 p.
3. Методические рекомендации по средне- и крупномасштабной магнитной картографии. - Л.: НПО «Рудгеофизика», 1990. -84 с.
4. Автеньев Г.К., Ерофеев Л.Я., Сысоев А.Т. Аналитическая связь между магнитными аномалиями и формой рельефа местности // Известия Томского политехнического института. -1976. - Т. 260. - С. 63-67.
5. Долгаль А.С. Магниторазведка: компьютерные технологии учета влияния рельефа местности. - Пермь: Перм. гос. нац. ис-след. ун-т, 2014. - 92 с.
6. Иванов В.К. Учет влияния рельефа при магниторазведке // Разведка и охрана недр. - 1977. - № 5. - С. 42-47.
7. Ахметов Е.М., Нусипов Е.Н., Сологуб Т.В. Учет влияния рельефа в аэромагниторазведке // Развитие методов обработки и интерпретации геофизической информации. - Алма-Ата: Казахский политех-кий ин-тут, 1991. - С. 60-70.
8. Ремпель Г.Г. Актуальные вопросы методики введения поправок, связанных с рельефом местности, в данные гравиразведки и магниторазведки // Изв. АН СССР. Сер. Физика Земли. -1980. - № 12. - С. 75-89.
9. Pilkington M., Thurston J.B. Draping corrections for aeromag-netic data: line-versus grid-based approaches // Exploration Geophysics. - 2001. - V. 32. - № 2. - P. 95-101.
10. Cowan D., Gordon C. Drape-related problems in aeromagnetic surveys: the need for tight-drape surveys // Exploration Geophysics. - 2003. - V. 34. - № 1. - P. 87-92.
11. Гравиразведка: справочник геофизика / под ред. Е.А. Мудре-цовой, К.Е. Веселова. - М.: Недра, 1990. - 607 с.
12. LaFehr, T.R., Nabighian, M.N. Fundamentals of Gravity Exploration. - Tulsa: SEG, 2012. 218 p.
13. Бычков С.Г., Долгаль А.С., Симанов А.А. Вычисление аномалий силы тяжести при высокоточных гравиметрических съемках. - Пермь: УрО РАН, 2015. - 142 с.
14. Геология и рудоносность Норильского района / О.А. Дюжи-ков, В.В. Дистлер, Б.М. Струнин и др. - М.: Наука, 1988. -498 с.
15. New standards for reducing gravity data: the North American gravity database / W.J. Hinze, C. Aiken, J. Brozena, B. Coakley, D. Dater, G. Flanagan, R. Forsberg, T.G. Hildenbrand, R. Keller, J. Kellogg, R. Kucks, X. Li, A. Mainville, R. Morin, M. Pilkington, D. Plouff, D. Ravat, D. Roman, J. Urrutia-Fucugauchi, M. Veronneau, M. Webring, D. Winester // Geophysics. -2005.- V. 70. - №4. - P. 25-32.
16. Старостенко В.И., Пятаков Ю.В., Исаев В.И. Решение прямых задач гравиметрии для сферических аппроксимирующих тел. Тестирование алгоритмов // Известия Томского политехнического университета. - 2013. - Т. 322. - № 1. - С. 35-39.
17. Долгаль А.С., Христенко Л.А. Применение эмпирической мо-довой декомпозиции при обработке геофизических данных // Известия Томского политехнического университета. Инжиниринг георесурсов. - 2017. - Т. 328. - № 1. - С. 100-108.
18. Dolgal A.S., Muravina O.M., Hristenko L.A. The reduction of the magnetic field within development areas of the plateaubasalts // Geoinformatics 2017. - Kiev, Ukraine, 15-17 May 2017. -11143_ENG. URL: http://geoinformatics.org.ua/eng/conferen-ces/pages-and-navigation/geoinformatics-2017 (дата обращения 26.12.2018).
19. Rasmussen T.M., Thorning L. Airborne geophysical surveys in Greenland in 1998 // Geology of Greenland Survey Bulletin. -1999. - № 183. - P. 34-38.
20. Cuma M., Zhadanov M.S. Continental-scale joint inversion of Alaska and Yukon gravity and magnetic data // First Break. -September 2017. - V. 35. - № 9. - P. 55-62.
Поступила 09.01.2019 г.
Информация об авторах
Долгаль А.С., доктор физико-математических наук, главный научный сотрудник лаборатории геопотенциальных полей Учреждения Российской академии наук «Горный институт» Уральского отделения РАН.
UDC 550.382.4
EFFECT OF EARTH SPHERICITY ON THE RESULTS OF AEROMAGNETIC MEASUREMENTS OVER PLATOBASALTS (ON THE EXAMPLE OF NORILSK REGION)
Alexander S. Dolgal,
dolgal@mi-perm.ru
Mining Institute of the Ural Branch Russian Academy of Sciences, 78-a, Sibirskaya street, Perm, 614007, Russia.
The research relevance is caused by the expediency of using model concepts adequate to real physical and geological conditions, while reducing the results of magnetic prospecting in the areas of development of platobasalts.
The aim of the research is to evaluate the influence of the Earth's sphericity on the anomalous magnetic field under conditions of a sharply dissected relief of the day surface which is composed of magmatic mountain rocks that have an increased magnetization (on practical materials).
Methods: calculation of corrections for the influence of magnetic terrain by solving the direct problem of magnetic exploration for a region with a radius of 20 km in two versions: for flat and spherical models of the Earth. Comparison of the results and spatial-statistical analysis of the identified differences were carried out. It was assumed that the magnetization vector of the rocks forming the relief is directed along the modern geomagnetic field of the region. Variable laterally effective magnetization of the upper part of the section was determined by solving the linear inverse problem of magnetic exploration.
Results. The paper describes briefly the problem of taking into account the influence of the terrain relief in measuring the magnetic field over volcanogenic formations. On the example of materials of a large-scale high-precision aeromagnetic survey, performed in the central part of the Norilsk ore region, the nonlinear signal distortions associated with the curvilinearity of the earth's surface were shown. The amplitude of these distortions (absolute value) reaches 100-200 nT and more, with relatively small differences in elevations for different models. This is due to the high values of the vertical gradient of the magnetic field.
Findings. In order to suppress the intense noise anomalies caused by magnetic rocks composing the upper part of the geological section, the author proposed the calculation of topographic corrections using the spherical model of the Earth. It also seems promising to use a statistical approach, implicitly taking into account the shape of the earth's surface, to isolate anomalies of «relief» nature. A comparison of the results obtained and a spatial-statistical analysis of the revealed differences are performed.
Key words:
Platobasalts, terrain, aeromagnetic survey, anomaly, topographic correction, reduction, model of the Earth.
The author expresses deep gratitude to the director of the NF VSEGEI Cand. Sc. F.D. Lazarev and the chief geophysicist of the NF VSEGEI P.V. Kirpluk, who tested in production volumes the introduction of corrections on the influence of the magnetic terrain as a result of high-precision aeromagnetic observations. The results obtained are included in the report on the object «Advanced geophysical surveys within the Khantaysk-Rybninsky uplift in order to assess the prospects of copper-nickel mineralization». The research was supported by the RFBR (projects no. 19-05-00654A, 18-05-70094 «Arctic Resources»).
REFERENCES
1. Ariskin A.A. Formatsii platobazaltov Zemli. Okeanicheskie plato i kontinentalnye trappy. Lektsiya №17 [Formations of the Earth plateau-basalts. Ocean plateaus and continental trapps. Lecture no. 17]. Available at: http://wiki.web.ru/images/b/b0/Magma-tic_Petrology_Lec11_2013.pdf_(acessed 17 July 2018).
2. Krivolutskaya N.A. Siberian Traps and Pt-Cu-Ni deposits in the Noril'sk area. Springer, 2016. 361 p.
3. Metodicheskie rekomendatsii po sredne- i krupnomasshtabnoy magnitnoy kartografii [Guideline in average- and largescale magnetic mapping]. Leningrad, Rudgeofizika Publ., 1990. 84 p.
4. Avtenev G.K., Erofeev L.Ya., Sysoev A.T. Analiticheskaya svyaz mezhdu magnitnymi anomaliyami i formoy relyefa mestnosti [Analytical relation between magnetic anomalies and local topography]. Izvestiya Tomskogo politekhnicheskogo instituta, 1976, vol. 260, pp. 63-67.
5. Dolgal A.S. Magnitorazvedka: kompyuternye tekhnologii ucheta vliyaniya relyefa mestnosti [Magnetic survey: computer technologies of considering local topography effect]. Perm, Perm State National Research University Press, 2014. 92 p.
6. Ivanov V.K. Uchet vliyaniya relyefa pri magnitorazvedke [Considering relief effect at magnetic survey]. Razvedka i ohrana nedr, 1977, no 5, pp. 42-47.
7. Akhmetov E.M., Nusipov E.N., Sologub T.V. Uchet vliyaniya relyefa v aeromagnitorazvedke [Considering relief effect in aero-
magnetics]. Razvitie metodov obrabotki i interpretatsii geofi-zicheskoy informatsii [Development of methods of processing and interpreting geophysical information]. Alma-Ata, Kazakhstan Polytechnic Institute Publ., 1991. pp. 60-70.
8. Rempel G.G. Aktualnye voprosy metodiki vvedeniya popravok, svyazannykh s relyefom mestnosti, v dannye gravirazvedki i magnitorazvedki [Relevant issues in technique of corrections, related to location topography, in the data of gravitation prospecting and magnetic survey]. Izv. AN SSSR. Ser. Fizika Zemli, 1980, no. 12, pp. 75-89.
9. Pilkington M., Thurston J.B. Draping corrections for aeromag-netic data: line-versus grid-based approaches. Exploration Geophysics, 2001, vol. 32, no 2, pp. 95-101.
10. Cowan D., Gordon C. Drape-related problems in aeromagnetic surveys: the need for tight-drape surveys. Exploration Geophysics, 2003, vol. 34, no 1, pp. 87-92.
11. Gravirazvedka: spravochnik geofizika [Gravitation prospecting: the reference book of a physicist]. Ed. by E.A. Mudretsova, K.E. Veselov. Moscow, Nedra Publ., 1990. 607 p.
12. LaFehr T.R., Nabighian M.N. Fundamentals of Gravity Exploration. Tulsa, SEG, 2012. 218 p.
13. Bychkov S.G., Dolgal A.S., Simanov A.A. Vychislenie anomaliy sily tyazhesti pri vysokotochnykh gravimetricheskih semkah [Calculations of gravity anomaly at high-accuracy gravimetrical survey]. Perm, UrO RAN Publ., 2015. 142 p.
14. Dyuzhikov O.A., Distler V.V., Strunin B.M. Geologiya i rudonos-nost Norilskogo rayona [Geology and ore-bearing capacity of Norilsk region]. Moscow, Nauka Publ., 1988. 498 p.
15. Hinze W.J., Aiken C., Brozena J., Coakley B., Dater D., Flanagan G., Forsberg R., Hildenbrand T.G., Keller R., Kellogg J., Kucks R., Li X., Mainville A., Morin R., Pilkington M., Plouff D., Ravat D., Roman D., Urrutia-Fucugauchi J., Veron-neau M., Webring M., Winester D. New standards for reducing gravity data: The North American gravity database. Geophysics, 2005, vol. 70, no. 4, pp. 25-32.
16. Starostenko V.I., Pyatakov Yu.V., Isaev V.I. Solving direct gravimetric problems for spherical approximating bodies. Algorithm Testing. Bulletin of the Tomsk Polytechnic University. Geo Assets Engineering, 2013, vol. 322, no. 1, pp. 35-39. In Rus.
17. Dolgal A.S., Hristenko L.A. The application of method of empirical mode decomposition (EMD) in processing of geophysical data. Bulletin of the Tomsk Polytechnic University. Geo Assets Engineering, 2017, vol. 328, no. 1, pp. 100-108. In Rus.
18. Dolgal A.S., Muravina O.M., Hristenko L.A. The reduction of the magnetic field with-in development areas of the plateaubasalts. 16th International Conference Geoinformatics - Theoretical and Applied Aspect. 2017 Available at: https://www.scopus.com/ inward/record.uri?eid=2-s2.0-85040014521&partne-rID=40&md5=38ad3b5bd6b4d434eb6f4422dcfa41e5 (accessed 26 December 2018).
19. Rasmussen T.M., Thorning L. Airborne geophysical surveys in Greenland in 1998. Geology of Greenland Survey Bulletin, 1999, no. 183, pp. 34-38.
20. Cuma M., Zhadanov M.S. Continental-scale joint inversion of Alaska and Yukon gravity and magnetic data. First Break, September 2017, vol. 35, no. 9, pp. 55-62.
Received: 9 January 2019.
Information about the authors
Alexander S. Dolgal, Dr. Sc., senior researcher, Mining Institute of the Ural Branch Russian Academy of Sciences.
