Размещено на http://www.allbest.ru/

Реферат

Темой данной курсовой работы является тема «Обработка и интерпретация результатов геофизических работ по Неманско-Усинскому участку».

Курсовая работа содержит 5 глав, 5 рисунков, общий объем 29 листов. При ее написании использовано 22 источника.

Ключевые слова: Неманско-Усинский участок, магниторазведка, гравиразведка, трубка взрыва, геофизическая аномалия, оцифровка аномалий.

Содержание

Введение

1. Физико-географическая характеристика района работ

2. Геолого-геофизическая характеристика района работ

2.1 Геофизическая изученность

2.2 Физические свойства горных работ

2.3 Стратиграфия

3. Методика проведения геофизических работ

3.1 Методика гравиметрической съемки

3.1.1 Методы топографической привязки пунктов наблюдения

3.1.2 Полевой контроль и уравнивание наблюдений

3.1.3 Способы учета влияний рельефа

3.2 Методика проведения магнитной съемки

4. Обработка результатов полевых геофизических работ

4.1 Оцифровка геофизических аномалий

4.2 Адаптация программы по решению прямой задачи магниторазведки

4.3 Адаптация программы по решению прямой задачи гравиразведки

5. Анализ результатов обработки

5.1 Анализ результатов обработки данных магниторазведки

5.2 Анализ результатов обработки данных гравиразведки

Заключение

Список использованных источников

Приложение 1 Обзорная карта Неманско-Усинского участка

Приложение 2 Стратиграфическая колонка

Приложение 3 Геолого-геофизические разрезы

Введение

Темой данной курсовой работы является «Обработка и интерпретация результатов геофизических работ по Неманско-Усинскому участку».

Детальные геофизические работы на Неманско-Усинской группе гравимагнитных аномалий выполнялись в пределах Столбцовской равнины, относящейся к смежным с Копыльской грядой равнинам подзоны Белорусской гряды и Гомельского Полесья Полесской равнины. Административно участки работ располагались в пределах Столбцовского, Узденского районов Минской и Гомельского, Жлобинского районов Гомельской области.

Геофизические работы выполнялись в соответствии с Программой ускорения геологоразведочных работ по развитию минерально-сырьевой базы Республики Беларусь на 1996-2000 годы для решения следующих задач:

1. Сбор, систематизация и анализ материалов геолого-геофизических работ на Неманско-Усинской площади для уточнения её геологического строения.

2. Выявление объектов трубочного типа на Неманско-Усинской группе гравиметрических и магнитных аномалий.

3. Оценка формы и размеров объектов трубочного типа.

Решение поставленных задач выполнено следующим комплексом работ:

1. Разбраковка магнитных аномалий двумя ортогональными профилями магнитных наблюдений с шагом 25 м; гравиметрических - гравиметрическими наблюдениями с шагом 50 м.

2. Заверка магнитных аномалий магнитной съемки по сети 50х20 м; двумя профилями гравиметрических наблюдений шагом 50 м.

3. Заверка гравиметрических аномалий гравиметрическими профилями с шагом 50 м, магнитными наблюдениями с шагом 25 м.

В процессе выполнения курсовой работы проведена обработка и интерпретация данных грави- и магнитной съемки при помощи ЭВМ.

1. Физико-географическая характеристика района работ

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

Площадь работ расположена в пределах Столбцовского и Узденского районов Минской, Жлобинского и Гомельского районов Гомельской областей Беларуси.

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

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

Гидрографическая сеть района работ относится к бассейну Балтийского моря и представлена частью верхнего течения р. Неман и его притоков. Долина р. Неман хорошо выработана, пойма,. большей частью заболоченная, ровная с редкими старицами, русло извилистое с меандрами, ширина составляет 25-30 м в отдельных местах до 50 м, глубина 1.5-2.0 м, скорость течения 0.5-0.6 м/с. Неман принимает притоки Уссу, Сулу, Яченку, Ольховку, Красноперовку, Крыницу, Залужанку, Жатеревку, Говезнянку. Все притоки, за исключением Уссы и Сулы, представляют мелкие ручьи и в большинстве случаев канализованы. Расположение долин относительно сети геофизических наблюдений и их повсеместная заболоченность представляли значительное препятствие при производстве полевых топогеодезических и геофизических наблюдений.

Болота на участке работ в большинстве относятся к низинному типу и приурочены к поймам рек.

Ледостав наступает в середине декабря и продолжается от 100 до 130 дней. Вскрываются реки в середине марта, начале апреля.

Климатические условия района работ характеризуются как умеренно-теплые и влажные, наиболее равномерные в Беларуси, среднегодовая температура составляет +5.5⁰С. Снежный покров достигает максимальной мощности 0.5 м и более м в феврале-марте.

Через Столбцы с северо-востока на юго-запад проходит транспортная магистраль Минск-Брест, железная дорога и автомобильное шоссе. От нее отходят автодороги районного и республиканского значения с твердым покрытием. В районе работ имеется и густая сеть грунтовых дорог, которые с периодом дождей, снегопадов и весенней распутицы в большинстве становятся труднопроходимыми [см Приложение 1].

Физико-географические условия Жлобинско-Уваровичской площади определяются её отношениями к Гомельскому Полесью и южной части Чечерской равнины Предполесской провинции. Гомельское Полесье и южная часть Чечерской равнины в пределах участка работ представляет собой водно-ледниковую и озёрно-алювиальную песчаную низину с древними террасами. Абсолютные высоты 100-150 м, в наиболее пониженных местах с почти плоским рельефом находятся болота и заболоченные участки с редкими песчаными дюнами, небольшими холмами, озёрными котловинами.

Гидрографическая сеть Жлобинско-Уваровичского участка относится к бассейну Черного моря и представлена рекой Днепр, её притоками Березиной, Сожью и их многочисленными мелкими притоками: Ольса, Ола, Добысна, Друть, Уза и др. Долины рек, речек и отчасти ручьёв хорошо выработаны, поймы большей частью заболочены, русла извилистые с меандрами, старицами. Ширина поймы Днепра в районе г. Рогачёва 3-6 км, ниже по течению достигает 10 км. Поймы притоков Березины, Друти и др. достигают 5-6 км, в большей части заболочены и покрыты кустарниковой растительностью. Замерзают реки в конце ноября - начале декабря, вскрываются в конце марта - начале апреля.

Климат теплый, неустойчиво-влажный средние температуры января от -4⁰С до -36⁰С. Максимально высокие в июне - до +38⁰С. Среднегодовая температура составляет +4.5⁰С. Снежный покров достигает максимальной мощности 0.6 и более метров в феврале-марте.

Через участок работ проходит транспортная магистраль Минск-Гомель, железная дорога и автомобильное шоссе. От неё отходят дороги областного и районного значения с твердым покрытием. В районе работ имеется и густая сеть, за исключением пойменных и заболоченных частей, грунтовых дорог, которые однако, в период дождей, снегопадов и весенней распутицы труднопроходимы. Согласно существующей классификации территории Неманско-Усинской и Жлобинско-Уваровичского участков разделяются на категории III-60%; IV-40%.

Полевые геофизические и топографические работы выполнялись с полевых баз, расположенных в г. Узде и д. Тихиничи Рогачёвского района. Полевые работы были начаты в декабре 1997 года и завершены в октябре 2000 года [1].

2. Геолого-геофизическая характеристика района работ

На площади работ в 1971-1975 годах выполнены среднемасштабные геологические и гидрогеологические съемки (Стасевич А.И., Левчич В.В., Шпаков О.Н., Хоменко В.Н.).

Более детальная гидрогеологическая и инженерно-геологическая съемка выполнена в Столбцовском районе для водоснабжения г. Столбцы и др. объектов агропромышленного комплекса.

В северо-западной части Неманско-Усинского участка в пределах Околовской структурно-металлогенической зоны в 1976-1985 годах Стасевичем А.И. выполнена детальная геологическая съемка для оценки её перспектив на поиски руд черных и цветных металлов. В 1985 году завершено глубинное геологическое картирование масштаба 1:50000 Околовский структурно-металлогенической зоны (Стасевич А.И., Гришко А.И. и др.). Построен ряд геологических карт, в т.ч. геологическая кристаллического фундамента, тектоническая схема кристаллического фундамента, прогнозная карта и др. Был дан прогноз на поиски руд чёрных и цветных металлов.

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

Мощность четвертичных отложений в пределах Неманско-Усинского участка изменяется от 79 до 193 метров. На большей части территории участка работ четвертичные отложения трансгресивно залегают на архей-протерозойских образованиях. Спорадическое распространение имеют отложения мела, неогена, палеогена, мощность которых изменяется в весьма широких пределах (0-50 м).

В северо-западной и северо-восточной части участка вскрыты отложения ратайчицкой свиты венда. В северо-восточной части участка скважиной 8ж вскрыты отложения горбашевской свиты венды и пинской свиты верхнего рифея.

По результатам глубинного геологического картирования, геофизических съёмок масштабов 1:200000 и 1:50000 выполнены научно-тематические работы.

В 1974 году Ковалёвым П.А. (БГРЭ) составлен макет карты перспектив алмазоносности территории Белоруссии. В 1976 году опубликована докладная записка Махнача А.С., Веретенникова Н.В., Шкуратова В.И. «К проблеме алмазоносности территории Белоруссии в свете новых данных (о необходимости постановки научно-исследовательских работ на возможные коренные источники алмазов и рудоносные дифференцированные траппы Белоруссии)» (ИГиГ АН БССР, БелНИГРИ), 1976 год. В этом же году опубликована книга «Тектоника Белоруссии» с тектонической картой под редакцией Гарецкого Р.Г. В своих работах в 1973 году Горелик З. А. и др. доказали, что сочленении Белорусского кристаллического массива со всеми отрицательными структурами проходит по глубинным разломам, находящим отражение в нижней части осадочного чехла. Кристаллический массив (по Бондаренко Б.В.) обособлялся вплоть до кайнозоя, начиная с байкальского этапа тектогенеза, что подчёркивает важность мезокайнозойский движений.

В связи с выявлением на Жлобинской седловине трубок взрыва активизировались научные исследования ИГН АН РБ. БелНИГРИ и производственных организаций по темам, связанным с региональным изучением глубинного строения земной коры по всей западной части Восточно-Европейской платформы [1].

2.1 Геофизическая изученность.

До 60-х годов геофизические исследования в пределах Неманско-Усинского участка носили эпизодический характер, выполнялись в маршрутном или мелкомасштабном площадном варианте. С 1970 года Геофизической экспедицией в описываемом районе проводятся планомерные крупномасштабные съёмки масштаба 1:50000 с целью выбора и подготовки площадей для глубинного геологического картирования применительно к масштабу 1:50000. По результатам этих работ получена качественно новая информация о структуре поверхности и внутреннем строении кристаллического фундамента участка работ. Приведённые в таблице съёмки сопровождались буровыми работами. Керновый материал подвергался всестороннему анализу, в т.ч. выполнен большой объём определений физических свойств горных пород. Партией № 40 Геофизической экспедиции выполнена систематизация и обобщение физических параметров горных пород.

Все вышеперечисленные работы проводились до выявления на территории Белоруссии трубок взрыва. Детальному изучению подверглись, в основном интенсивные гравимагнитные аномалии. На слабоинтенсивные локальные аномалии особого внимания не обращалось, детальными геофизическими работами они не изучались и буровые работы с целью обнаружения трубок взрыва и алмазов в них не производились [1].

2.2 Физические свойства горных пород

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

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

а) существенно основные породы - габбро, нориты, габбро-диабазы, незначительно измененные; средняя плотность - 3.0-3.1 г/см³, магнитная восприимчивость - 10·10^-3 ед. СГС;

б) породы условно среднего состава (или гранитизированные породы класса «а») - амфиболовые, биотит-амфиболовые, биотит-пироксеновые плагиогнейсы, чарнокиты, диориты; средняя плотность - 2.75-2.85 г/см³, магнитная восприимчивость 1500-2200·10^-6 ед. СГС;

в) кислые породы - граниты, гранито-гнейсы биотитовые, мигматиты, аплиты; плотность преобладает 2.65-2.70 г/см³, магнитная восприимчивость 10-20 до 500·10^-6 ед. СГС;

г) амфиболиты - выделены по аномальному соотношению плотности 2.95-3.0 г/см³ и магнитной восприимчивости 50-250·10^-6 ед. СГС;

д) среднего состава - диориты, гранодиориты, обогащенные магнетитом, разуплотненные; невысокая плотность - 2.7-2.75 г/см³, значительная магнитная восприимчивость 2000-2500·10^-6 ед. СГС.

Характеристика физических полей определяется горными породами указанных классов:

1) породы класса «в» характеризуются минимальными значениями физических полей, «мозаичным» рисунком аномалий;

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

3) породы среднего состава картируются полями промежуточного слоя.

Породы осадочного чехла имеют магнитную восприимчивость от 0 до 10-20·10^-6 ед. СГС и характеризуются неоднородностью магнитной восприимчивости, плотность близка к 2.0 г/см³, доломиты, известняки характеризуются плотностью 2.6 - 2.7 г/см³. Подробнее физические свойства осадочного чехла и кристаллического фундамента освещены в таблице 1.2.

Кимберлиты со средней плотностью (2.4 г/см³) и магнитной восприимчивостью 1000-1500·10^-6 ед. СГС дадут разные типы аномалий:

а) кимберлиты в осадочном чехле - положительные аномалии магнитного поля и силы тяжести при интенсивности в десятки и первые сотни нТл, 0.3-0.8 мГал;

б) кимберлиты в кристаллическом фундаменте - слабоинтенсивные положительные магнитные либо локальные магнитные минимумы интенсивностью в сотни и более нТл (трубка среди основных пород); в любом случае - локальный тип силы тяжести.

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

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

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

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

Выделенные линейно-вытянутые структуры отличаются вещественным составом и формационной принадлежностью слагающих их комплексов пород. Буровыми работами доказано, что в максимумах поля G развиты комплексы пород преимущественно метабазитового ряда со средней плотностью от 2.72 до 2.86 г/см³ и интенсивностью намагниченности от 0 до 300010^-6 ед. СГС. Здесь вскрыты амфиболиты, пироксеновые, амфиболовые и биотитовые плагиогнейсы, кристаллические сланцы, часто с гранатом, и переходные разности между ними. Повсеместно с гнейсами основного состава ассоциируют согласные тела габбро, габбро-норитов, пироксенитов, габбро-амфиболитов и диабазов.

Зона регионального гравитационного минимума и ослабленного по напряженности слабо дифференцированного магнитного поля обусловлена породами гнейсово-мигматитового комплекса с преобладанием интенсивно гранитизированных разностей гнейсов кислого и среднего состава и мигматитов. Однако, и в этой зоне встречаются отдельные тела измененных габбро и амфиболитов, выделяемые по локальным гравитационным и магнитным аномалиям наибольшей интенсивности, что свидетельствует о происходивших магматических интрузивных процессах и в этой области. Для комплекса характерны средние значения плотности 2.65 - 2.67 г/см³ и магнитной восприимчивости от 0 до 60010^-6 ед. СГС.

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

Юго-восточную большую часть Неманско-Усинского участка занимает Бобовнянский гранито-гнейсовый блок, для которого характерны отрицательные спокойные магнитные и гравитационные поля [1].

2.3 Стратиграфия

Архей - нижний протерозой.

Околовская серия (AR₂ок). Слагающие ее породы представлены амфиболитами, амфибол-биотитовыми, роговообманково-биотитовыми, биотитово-роговообманковыми, куммингтонитовыми, пироксено-амфиболовыми, биотитово-куммингтонитовыми, гранатово-биотитовыми, гранатово-силлимонитово-биотитовыми и гранатово-биотитово-актинометовыми плагиогнейсами; гранатово-биотитовыми, двуслюдеными, биотитовыми, амфиболовыми, хлорито-биотитовыми сланцами и железистыми кварцитами. В связи со столь широким разнообразием отложений и неравномерностью их распространения Околовская эффузивно-осадочная серия разделяется на три свиты (толщи): нижнюю (Гуменовщинскую) - амфиболитовую, среднюю (Шашковскую) - железорудную, верхнюю (Яченскую) - гнейсовую (11). В генетическом отношении они отнесены (Махнач А.С. и др. 1974 г.) к образованиям соответственно: метабазитовой, железисто-кремнисто-метабазитовой и сланцевой формации.

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

В пределах Бобовнянского (AR₁bb) выступа на участке работ развиты биотитовые, амфибол-биотитовые, микроклин-плагиоклазовые гнейсограниты, гранитогнейсы, граниты. Гнейсы роговообманковые занимают северную и восточную часть площади. Гнейсо-граниты амфибол-биотитовые, биотитовые архейского возраста распространены на юго-востоке участка, центральную часть занимают амфибол-биотитовые, биотитовые гранитогнейсы, в центральной и юго-западной части выступа прорванные интрузиями биотитовых гранитов Жуховичского комлекса верхнего архея. (Рис. 1.5).

На границе Бобовнянского блока и Околовской зоны имеют место интрузивные тела габброидов. Юго-восточнее тел габброидов в зоне Минского разлома, в южной части участка, в северо-восточном направлении располагаются тела Ужанской дайки нагорновского комплекса верхнего протерозоя диабазов (Приложение 2 ).

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

Осадочный чехол

Верхний протерозой, рифей

Белорусская серия. Представлена пинской свитой (R₂pn), вскрытой скважиной №8ж, находящейся в северо-восточном углу участка работ. Отложения пинской свиты мощностью 12.8 м залегают с угловым несогласием и стратиграфическим перерывом на породах кристаллического фундамента. Ближайшими скважинами на участке не обнаружены. Распространены в северо-восточном и восточном направлении (). Отложения свиты, вскрытые скважиной № выполнены алевритами, кварцевыми песчаниками.

Волынская серия. Горбашевская свита (V₂gr). Отложения мощностью 41.5 м вскрыты скважиной № 8, представлены песчаниками, развиты в северо-восточном и восточном направлении. Породы залегают со стратиграфическим перерывом на отложениях пинской свиты, перекрыты несогласно отложениями ратайчицкой свиты.

Ратайчицкая свита (V₂rt) представлена вулканогенно-осадочными породами: туффитами, туфогенными алевролитами, туфопесчаниками, алевролитами, глинами с прослоями туфогенных пород. Породы ратайчицкой свиты вскрыты и пройдены на полную мощность скважиной № (4.9 м), 49п (6.0 м) в северо-восточном углу участка и скважинами №№ 152, 167, 207, 160, 185, 155, 184, 206 (Рис. № 1.5). Мощность отложений в северо-восточном углу площади, вскрытых скважиной № 49п, составляет 6.0 м. В юго-восточном направлении в скважине № на листе мощность достигает 40 м. В северо-западной части площади работ мощность ратайчицких отложений, представленных, в основном, туффитами, достигает 56.6 м. Породы ратайчицкой свиты залегают трансгрессивно на породах горбашевской свиты в северо-восточной части участка и на архей-протерозойских образованиях на северо-западе площади работ. Перекрыты отложения венда, в основном, меловыми отложениями турона или сеномана. В отдельных случаях четвертичными отложениями.

Палеозой.

Верхний мел. Сеноманский ярус. (K₂s). Отложения сеноманского яруса вскрыты скважинами в южной, северо-западной и северо-восточной части площади работ. Представлены отложения песками кварцевыми, глауконитовыми, в отдельных местах конгломератом (скв. № 45), песчаником карбонатным, мергелем, запесоченным мелом, переходящим в верхней части в мел глинистый, мел. Глубины залегания сеноманских отложений изменяются от 102.8 до 174.6м, мощность отложений достигает 76.8 м. Залегание пород несогласное на породах кристаллического фундамента, пинской свиты белорусской серии рифея, горбашевской, ратайчицкой свит волынской серии венда. Перекрываются согласно - туронскими отложениями, либо несогласно - четвертичными.

Туронский ярус (K₂t) представлен двумя отделами. Нижний отдел сложен глинистым мелом, мелоподобными мергелями. Верхний отдел - чистым мелом, в низах разреза, при отсутствии нижнего отдела - глинистым мелом. Туронский ярус имеет более широкое распространение, чем сеноман. Его отложения встречены даже в центре площади работ в скв. №36. Наибольшие мощности отложений турона (до 60 м) пройдены скважинами в северо-западной части участка работ. Преобладают мощности: в северо-западной, западной части - 10 - 30 м, в северо-восточной, северной, восточной - 3 - 10 м. Глубина залегания кровли отложений турона изменяется от 90 до 190 м. Залегают отложения турона согласно на сеноманских отложениях, несогласно - на горных породах протерозоя и архея. Перекрываются несогласно в южной, юго-восточной части участка отложениями киевской, харьковской свит палеогена, белицкой серии неогена на юге и северо-западе центральной части участка. На большей части территории перекрываются несогласно отложениями антропогена.

Кайнозой.

Палеоген.Эоцен. Киевская свита (Р₂kv). Отложения представлены кварцевыми песками с глауконитом. Встречен в южной части участка в двух скважинах №№, 9.Мощность соответственно 49.8 м и 2 м. В скв. №72 пески Р₂kv непосредственно несогласно залегают на породах кристаллического фундамента. В скв. № 9 на отложениях мела К₂s₂-t. Перекрываются в скв. № 27 глинами харьковской свиты, в скв. № 9 - четвертичными отложениями. Кровля поверхности отложений киевской свиты находится на глубине 97.2 м.

Олигоцен. Харьковская свита. (Р₂hr). Отложения глин харьковской свиты мощностью 11.4 м встречены на глубине 85.8 м на юге площади работ в скв. №. Залегают согласно на отложениях киевской свиты, перекрываются несогласно четвертичными отложениями.

Неоген. Миоцен. (N₁). Отложения не расчленены. Встречены на юге и западе центральной части участка в скважинах №№34, 70. Представлены глинами. Залегают несогласно на горных породах кристаллического фундамента. Кровля отложений в скважине №67 находится на глубине 123.6 м., в скважине №56 - на глубине 103.4 м. Перекрываются в обоих случаях четвертичными отложениями.

Четвертичные отложения (Q) развиты повсеместно. Мощность отложений изменяется от 78 м до 190 м, преобладают мощности 100 - 140 м. Отложения представлены ледниковыми, водно-ледниковыми, озерно-ледниковыми, озерно-болотными, флювиогляциальными, аллювиальными и др. фациями березинской, днепровской и московской морен. Залегают несогласно преимущественно на горных породах кристаллического фундамента или меловых отложениях [см Приложение 2].

3. Методика проведения геофизических работ

3.1 Методика гравиметрической съемки

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

Исследований с целью получить картину распределения аномалий силы тяжести называется гравиметрической съемкой [2].

Изучаемая площадь покрыта высокоточной гравиметрической съемкой масштаба 1:10000 [1].

Гравиметрические наблюдения осуществляют отдельными рейсами. Каждый рейс представляет собой непрерывную последовательность наблюдений, объединенную общим учетом смещения нуль-пункта гравиметра [2].

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

Пункты опорной сети необходимо располагать на площади исследований равномерно. Расстояния между пунктами опорной сети выбирают исходя из методики съемки и способа учета нуль-пункта. В зависимости от густоты рядовой сети расстояние между опорными пунктами колеблется от 1 до 20 км [2].

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

При этом каждый последующий рейс должен включать несколько пунктов, наблюденных в последующих рейсах.

Опорные рейсы планируются таким образом, чтобы они образовывали систему замкнутых полигонов [2].

Методика построения рейсов может быть двух типов: методика однократных наблюдений и методика наблюдений с повторением.

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

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

при резко нелинейном смещении нуль-пункта, недостаточной густоте опорной сети и т.д [2].

3.1.1 Методы топографической привязки пунктов наблюдения

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

- из триангуляционных построений;

- засечками (угловыми, линейными, линейно-угловыми и графическими);

- полярным способом;

- проложением теодолитных ходов и профильных линий.

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

Триангуляционные построения должны опираться не менее чем на 2 исходных пункта, количество треугольников в системе не должно превышать 20.

Определение координат прямой засечкой выполняется не менее чем с трех исходных пунктов, углы между направлениями при определенной точке 30-150 градусов. Определение координат обратной засечкой выполняется не менее чем по 4 исходным пунктам при условии, что определяемая точка не находится около окружности, проходящей через 3 исходных пункта.Комбинированная засечка выполняется не менее чем по 3 исходным пунктам.

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

Теодолитные ходы прокладываются между исходными пунктами преимущественно по заданному направлению, а также по дорогам, тропам и другим местам, удобным для линейных и угловых измерений.Длины линий измеряют радио- и светодальномерами, стальными лентами и рулетками [4]

3.1.2 Полевой контроль и уравнивание наблюдений

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

Контроль сьемки считается представительным, если он выполнен в достаточном объеме и все звенья рейсов проконтролированы(при измерениях одним гравиметром). Количество независимых наблюдений, выполненных для контроля качества работ, устанавливается техническим проектом в соответствии с требованиями по данному виду съемки.Контрольные наблюдения проводятся в объеме 5-10%.При наблюдениях на пунктах двумя или более гравиметрами одновременно количество повторных наблюдений устанавливается от 3% (при сьемках с большим количеством пунктов) до 5%, но их число не менее 50.

Существуют различные способы уравнивания опорных сетей.Они используются в геодезии при уравнивании триангуляционных и нивелирных ходов.Наиболее часто используют 2 способа:способ полигонов и способ угловых точек.

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

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

Eоп=м(M-m) (3.1.1);

Где М-средняя удаленность опорных пунктов от ближайших исходных пунктов, выраженная в числе приращений; m-среднее число наблюдений силы тяжести между двумя смежными пунктами.

Коэффициент пропорциональности м вычисляется по формуле:

М=Уд/(N-S) (3.1.2);

Где S-число сторон, образующих все полигоны опорной сети [3].

3.1.3 Способы учета влияний рельефа

Аномалия силы тяжести Дg-разность между наблюденным значением силы тяжести и ее нормальным значением г в одной и той же точке физической поверхности Земли в одно и то же время.Однако нормальные значения силы тяжести рассчитаны только для поверхности сфероида (г0), а измерения проводят в точках, не соответствующих им.Поэтому вычисляют г в точках наблюдения или определяют полное значение силы тяжести г0 в точке с теми же горизонтальными координатами, но на поверхности сфероида.Таким образом,

Дg=g-г=go-гп (3.1.3);

В пределах небольших перепадов высот считают, что г0 постоянно изменяется с высотой.С увеличением высоты на 1 м значение г0 уменьшается на 0,03086м/с.Аномалия в свободном воздухе будет:

Дg=g-г0+0,003086Н (3.1.4);

Где Н-высота точки наблюдения в метрах, (0,003086Н)-поправка за высоту.

Редукция Буге учитывает не только поправку за высоту, но и притяжение промежуточного слоя (между поверхностью Земли и поверхностью геоида). Гравитационное притяжение такого слоя 0,000000419уH, где у-средняя плотность пород слоя,Н-мощность слоя, равная высоте точки наблюдения.

Аномалия в редукции Буге будет:

Дg=g-г0+(0,3086-0,0419)З (3.1.5);

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

Обе поправки вводятся с одним и тем же значением плотности.При обработке съемок масштаба 1:50000 и мельче вычисляют аномалии силы тяжести при двух значениях плотности пород промежуточного слоя -2,3 и 2,67 г/см3, получая 2 карты аномалии силы тяжести.

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

(-0,3086+2*0,01419уВ)0,0001р (3.1.6);

где уВ-плотность морской воды,принимают равной 1,03 г/см3;р-глубина точки измерения.

Учет вариаций силы тяжести.

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

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

О сущности вековых вариаций силы тяжести говорят на основании теоретических предположений и некоторых косвенных наблюдений.Неотектонические движения вызывают значительные изменения высоты отдельных точек земной поверхности.Изменение высоты точки наблюдения на 1 м приводит к изменению силы тяжести на 0,2 мГл.Внутри Земли перемещения масс также обуславливают значительные вариации силы тяжести.

Периодические изменения силы тяжести вызываются изменением положения Земли относительно Луны и, менее, Солнца.Считая, что центр тяжести Земли в течение суток испытывает постоянное притяжение, периодические изменения силы тяжести, вызванные Луной, будем определять как разность значений силы тяжести на поверхности и в центре Земли.Приняв Мл за массу Луны, r1 и r за расстояние ее центра тяжести от пункта наблюдений В и от центра тяжести Земли, z1 и z -зенитные расстояния центра луны для пункта наблюдений и для центра Земли, R-радиус-вектор из центра Земли О в пункт наблюдений В, можно записать:

Дg=kMrR(3cos z-1)/r (3.1.7);

Аналогичная формула может быть записана и для Солнца. Определим численное значение коэффициента М(R/r) для Луны и Солнца. Известно что масса Луны Мл=Мз.81,6;

масса Солнца Мс=332000Мз.Для определения численных значений силы тяжести необходимо вычислить для данного момента времени зенитные расстояния z Луны и Солнца по формулам:

cos z=cosцcosт+sinцcosт(t-a) (3.1.8);

где ц-географическая широта пункта наблюдения;а-прямое восхождение;у-склонение Солнца (Луны) в момент наблюдения;t-момент наблюдения. Величины а и у берутся из астрономического ежегодника.Для учета лунно-суточных вариаций составлены специальные таблицы, по которым получают величины поправок Дg по известному в данный момент времени зенитному расстоянию z.

Проверка гравитационных аномалий также осуществлялась в два этапа - разбраковки и заверки. Разбраковка осуществлялась ортогональными профилями с шагом наблюдений равным 50 м. Для заверки подтверждённых и перспективных аномалий создавалась топооснова по сети 100х50 м или 50х50 м (для аномалии № БВ-72) если размер аномалии был соизмерим с сетью наблюдений 100х50 м. На профилях выполнялось техническое нивелирование с погрешностью не превышающей ±10 см и гравиметрические наблюдения шагом 50 м. Гравиметрические наблюдения и нивелирование выполнялись в условном уровне от центральной точки.

По результатам площадных гравимагнитных наблюдений намечались расчётные профили, наблюдения на которых выполнялись по двукратной методике шагом 50 м с опорой на центральную точку, выбранную при площадной съёмке. Точность площадной съёмки и гравиметрических наблюдений на расчетных профилях определялась контрольными наблюдениями, объём которых составил 4.6%. Погрешность площадной съёмки и наблюдений на расчетных профилях не превышает ±0.04 мГл при проектной ±0.04 мГл [1].

Плановая привязка осуществлялась по топокартам масштаба 1:10000, погрешность которой не превышает ±10 мГл. Концы расчётных профилей и центр аномалий закреплялись столбами [см прил 3].

3.2 Методика магнитной съемки.

Наземные магнитные съемки проводятся в масштабе 1:10000 и крупнее [3].

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

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

Сеть измерений и допустимая погрешность съемки обеспечивается соответствующей методикой полевых работ. Прежде всего решается вопрос о необходимом типе магнитометра, так как точность измерений с оптико-механическими и квантовыми приборами существенно разная. Решается вопрос и об опорной сети. Все современные магнитометры обладают изменением нуль-пункта во времени, причем у оптико-механических приборов оно может быть существенным и отличным от линейного. Опыт показывает, что близким к линейному можно считать смещение нуль-пункта у магнитометра М-27 в течение примерно 2-2,5ч.Поэтому при съемке надо контролировать смещение нуль-пункта через такие интервалы времени. Для этого вблизи или на самом участке съемки выбирают контрольный пункт и выполняют на нем измерения ежедневно перед началом работы и после ее окончания. Возвращение на этот пункт через каждые 2 часа нецелеобразно, поэтому, если точность съемки этого требует, на площади работ разбивают и взаимно увязывают сеть опорных пунктов [3].

В зависимости от размеров площади съемки, допустимой погрешности и характера рельефа местности опорная сеть может быть создана различными методами. Если площадь сьемки, выполняемой с оптико-механическими магнитометрами, значительна, то последовательно разбивают сначала каркасную. Затем заполняющую опорную сеть. Пункты каркасной сети, чаще всего расположенные на магистралях, являются исходными для последующего разбития заполняющей сети. Погрешность определения приращений поля для каркасной сети должна быть приблизительнов 4 раза, а для заполняющей сети -в 2 раза меньше, чем для рядовой съемки. Если площадь работ невелика, или если проводится съемка средней точности, ограничиваются одноступенчатой разбивкой сети. В данном случае погрешность должна быть в 2 раза ниже [3].

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

Опорная сеть должна быть увязана с исходным контрольным пунктом. Густота опорной сети должна быть такой, чтобы оператор при работе оказался вблизи какого-либо из опорных пунктов. Измерение и учет вариаций при разбивке опорной сети обязательны [3].

Проектом предусматривалось проверка выявленных при анализе магнитометрических материалов съемок предыдущих лет на Неманско-Усинской площади и выявленных по «Проекту на переинтерпретацию », 1977 г.

Проверка выявленных локальных аномалий осуществлялась в два этапа. На первом этапе выполнялась разбраковка магнитных аномалий, заключавшаяся в отыскании эпицентра аномалий, отработки через найденный центр аномалий двух ортогонально ориентированных профилей шагом 25 м с выходом в нормальное поле. При этом направление одного из них соответствовало азимуту съёмного маршрута предшествующих работ. Полевые материалы обрабатывались с введением поправки за вариации магнитного поля, полученных с помощью вариационной станции М-33 или магнитометра ММП-203 с регистрацией магнитного поля через 3-5 мин.

При подтверждении наличия перспективной аномалии на ней выполнялись дальнейшие заверочные работы. В зависимости от размера аномалии, определённого крестообразными разбраковочными профилями подготавливалась площадка наблюдений по сети 50х20 м, размером от 700х700 м, до 1500х1500 м. В отдельных случаях (аномалии: БВ-68, 69, 67), где аномалии располагались на близком расстоянии, для них создавалась общая площадка размером, обеспечивающим выход наблюдений магнитного поля до нормального уровня. Кроме введения поправки за вариации для достижения требуемой точности съёмки на крайних магистралях или на центральной, разбивалась двойными наблюдениями опорная сеть. [1]

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

После отработки площадки, обработки полевых материалов и построения магнитного поля намечались расчетные профиля длиной, обеспечивающей выход в нормальное поле. Наблюдения магнитного поля на расчетных профилях выполнялись шагом 25 м с введением поправок за вариации магнитного поля. Среднеквадратичная погрешность магнитометрических наблюдений составила ±2.125 нТл при допустимой проектной ±2.0 нТл. Проектная точность магнитометрических наблюдений не достигнута по причине неустойчивой работы приборов из-за их естественного износа. Объём контрольных наблюдений составил 4.6% при проектном 5.0% [1].

4. Обработка результатов полевых геофизических работ

4.1 Оцифровка геофизических аномалий

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

Оцифровка геофизических аномалий велась для правой части РП-1с шагом дискретизации 100 м. Значения Za и Ga были получены для 41 пикета при расстоянии между пикетами 50м. Результаты оцифровки представлены в таблице 3.

Таблица 3

Результаты оцифровки геофизических аномалий

ПК	Za, нТл	Ga, мГал	ПК	Za, нТл	Ga, мГал
80	-410	-5,4	122	-200	-1,13
82	-416	-4,95	124	-135	-0,8
84	-420	-4,4	126	-40	0,2
86	-411	-4,51	128	-65	1,0
88	-408	-4,65	130	-135	1,04
90	-410	-4,82	132	-210	1,1
92	-411	-5,1	134	-90	0,7
94	-412	-5,1	136	-265	0,01
96	-402	-5,2	138	-365	0,65
98	-404	-4,95	140	-380	-1,25
100	-400	-4,65	142	-382	-1,40
102	-411	-4,63	144	-385	-1,55
104	-415	-4,6	146	-390	-2,30
106	-410	-4,44	148	-395	-3,1
108	-395	-4,25	150	-385	-3,25
110	-410	-4,17	152	-375	-3,4
112	-405	-4,1	154	-371	-3,7
114	-381	-3,62	156	-365	-3,9
116	-340	-3,0	158	-362	-4,31
118	-297	-2,33	160	-360	-4,63
120	-255	-1,65

4.2 Адаптация программы по решению прямой задачи магниторазведки

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

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

В учении о земном магнетизме напряженность магнитного поля принято рассматривать по стставляющим прямоугольной системы координат (рисунок 1), в которой ось Х горизонтальна и направлена на географический север, ось У горизонтальна и направлена вниз, ось Z вертикальна и направлена вниз. Составляющие вектора напряженности Т по осям соответственно называются северной Х, восточной У и вертикальной Z; проекция вектора т на горизонтальную плоскость называется горизонтальной составляющей

H= х2+у2 (4.2.1)

Угол, образуемый Н и осью Х, называется склонением и обозеачается D.Угол между плоскостью ХОУ и вектором Т называется наклонением и обозначается I [3].

Аномалии магитного поля на Неманско-Усинском участке обусловлены, предположительно, трубками взрыва. Для определения элементов их залегания используем решение прямой задачи магниторазведки для вертикального стержня.

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

Zа= (4.2.2)

Полагая сечение стержня столь малым, что все магнитные массы сосредоточены в одной точке, можно рассматривать координаты его вершины и как координаты фиксированной точки; это позволяет вынести величину (z-)/r3 за знак интеграла. Тогда, обозначая через М=I*S магнитный момент стержня единичной длины, можем записать:

Ха=М*[-(X-)]/r3; Ya=M*[-Y-)]/r3; Za=M*[-(Z-)]/r3 (4.2.3)

Перенося начало координат в точку, являющуюся проекцией оси стержня на плоскость ХОУ и учитывая, что глубина залегания =h, получаем:

Ха=М*(-Х)/r3; Ya=M*(-Y)/r3; Za=T=M*h/r3; Ta= M/r3 (4.2.4);

где r=h2+x2+y2.

При произвольном направлении намагниченности составляющие напряженности магнитного поля Ха и У а имеют достаточно сложный вид и здесь не приводятся. Составляющая

Za=M*[h*sin I-(x cos A+ y sinA)*cos I]/r3 ( 4.2.5)

Где I-наклонение вектора намагниченности, т е угол между вектором I и его проекцией не плоскость ХОУ; А-угол между проекциями вектора I на ось х и на плоскость ХОУ.

По профилю, проходящему через экстремумы аномалий над осью вертикального стержня, Za и Ha упрощаются:

На=M*[(h* cos I+x* sin I)/ r3-(r-h)/(x*2*r)*cos I] (4.2.6);

Za=M*(h*sin I- x* cos I)/r3 (4.2.7).

Выражение (4.2.5) было запрограммировано на языке Фортран. При помощи программы, написанной на основе выражения, решалась прямая задача магниторазведки для вертикального стержня.

...