Решение прямой задачи Аэромагниторазведки для набора тел

Лабораторная работа №1

Вариант 5

Решение прямой задачи Аэромагниторазведки для набора тел

В 1967-1968 годах аэрогеофизическая группа ВИМС под руководством Э.Я. Островского начала первые работы по изучению радиоактивности месторождений нерадиоактивных руд. Первые же полеты над золото-серебряными и молибденовыми месторождениями Охотского района и Магаданской области показали, что рудные зоны неизменно выделяются с воздуха мощными аномалиями калия. Исследования оловянных, других месторождений выявили "зеркальную" картину в поведении калия и тория. В рудных зонах происходила заметная концентрация тория на фоне низкого содержания калия. Накопление тория при низком уровне калия оказалось признаком таких важных полезных ископаемых, как кимберлиты и алюминиевые руды - бокситы.

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

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

Jx=ӕX0 Jy=ӕY0 Jz=ӕZ

где X0, Y0 и Z0 – горизонтальные и вертикальная компоненты нормального магнитного поля Земли, определенные для района работ, ӕ-магнитная восприимчивость. Используя выражение для компонент магнитного поля шара

[pic 1]

Рисунок 1. Аномалия магнитных полей

Решение

Расчеты произведем с помощью он-лайн калькулятора https://geomag.bgs.ac.uk/data_service/models_compass/wmm_calc.html.

Точка с указанными в задании координатами представлена на карте ниже.

[pic 2]

Рисунок 2. Расчет магнитных полей

Полученные данные представлены ниже.

По рассчитанным компонентам магнитного поля вычисляют компоненты намагниченности.

Xa=[pic 6]

Ya=[pic 7]

Za=[pic 8]

Шар №1.

Jx=ӕX0=310* 24722= 7663820

Jy=ӕY0=310*848= 262880

Jz=ӕZ=310*38898= 12058380

Xa=[pic 9]

Ya=[pic 10]

Za=[pic 11]

Шар №2.

Jx=ӕX0=180* 24718=4449240

Jy=ӕY0=180*848=152640

Jz=ӕZ=180*38992= 7000560

Xa=[pic 12]

Ya=[pic 13]

Za=[pic 14]

Шар №3.

Jx=ӕX0=650* 24725= 16071250

Jy=ӕY0=650*848=551200

Jz=ӕZ=650*38904= 25287600

Xa=[pic 15]

Ya=26194013899[pic 16]