Обработка результатов измерений гравитационного поля по участку

КурсоваяПомощь в написанииУзнать стоимостьмоей работы

Обработка результатов измерений гравитационного поля по участку (реферат, курсовая, диплом, контрольная)

ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ОБРАЗОВАНИЮ ГОУ ВПО УРАЛЬСКИЙ ГОСУДАРСТВЕННЫЙ ГОРНЫЙ УНИВЕРСИТЕТ Факультет геологии и геофизики Кафедра математики Направление подготовки: 657 100 ПРИКЛАДНАЯ МАТЕМАТИКА Специальность: 73 000 «Прикладная математика»

Квалификация: Инженер-математик Специализация: Применение математического моделирования и математических методов в горно-геологических науках Курсовая работа

ТЕМА: Обработка результатов измерений гравитационного поля по участку Автор: Берсенев Ярослав Александрович Учебная группа: ПМ-07

Руководитель: Серков Владимир Алексеевич Екатеринбург

2011 г.

Постановка задачи

Выделить положение объектов с различной плотностью.

1. Определить модель поля (наличие региональной и локальной составляющих).

2. Определить аналитический вид региональной составляющей поля.

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

4. Вычислить аномальное поле.

5. Определить вид и параметры оптимальных для данного аномального поля фильтров.

6. Выполнить фильтрацию аномальной составляющей и построить модель аномального поля после фильтрации.

7. Выделить объекты с различной плотностью.

8. Оценить плотность горных пород на площади исследования.

9. Написать пояснительную записку следующего содержания.

— характеристика исходных данных;

— характеристика модели поля;

— обоснование выбора модели региональной компоненты;

— характеристика аномального поля и обоснование выбора фильтров;

— характеристика аномального поля и результаты интерпретации;

— описание физико-геологической модели участка исследования.

10. Графические приложения:

— модель измеренного поля;

— модель региональной составляющей;

— модель аномального поля до и после фильтрации;

— графики автокорреляционных функций;

— модель исследуемого участка

Геофизическая характеристика

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

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

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

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

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

Техника и методика наземных гравиметрических измерений

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

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

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

Масштаб съемки определяется общем случае расстоянием между профилям. Это расстояние должно укладывать примерно в 1 см отчетной карты. При этом масштаб съемки должен соответствовать стандартным масштабам, принятым в геодезии (1:10 000, 1:25 000, 1:50 000 и др.). В зависимости от масштаба съемки разделяются на мелкомасштабные (1:100 000 и мельче), среднемасштабные (1:25 000 — 1:50 000), крупномасштабные (1:10 000 и крупнее). Съемки с мелким масштабом проводятся с целью региональных исследований, среднеи крупномасштабные съемки используются для решения вопросов геологического картирования площади исследований, при поисковых и поисково-разведочных работах.

Точность проектируемых съемок определяется величиной амплитуды минимальных аномалий, представляющих поисковый интерес. Погрешность определения аномальных значений гравитационного поля не должна быть выше 0,33 — 0,2 амплитуды минимальных искомых аномалий. Сечение изолиний отчетной карты должно превосходить среднеквадратическую погрешность съемки в 2−2,5 раз.

Затем, в зависимости от того, в какой местности, равнинной или горной, проектируется проведение гравиметрических работ, выбранные технический минимальных аномалий, Представляющих поисковый интерес. Погрешность определения аномальных значений гравитационного поля не должна быть выше 0,33 — 0,2 амплитуды минимальных искомых аномалий. Сечение изолиний отчетной карты должно превосходить среднеквадратическую погрешность съемки в 2−2,5 раз. Затем, в зависимости от того, в какой местности, равнинной или горной, проектируется проведение гравиметрических работ, выбранные технический условия съемки корректируются в соответствии с инструкцией по гравиразведке. Согласно с инструкцией по гравиразведке съемке определенного масштаба должна соответствовать требуемая точность и густота сети наблюдений, которая характеризуется числом пунктов на 1 км съемки.

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

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

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

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

Полигональная система увязки опорных точек представляет собой систему взаимно пересекающихся полигонов, образующих в конечном итоге сеть треугольников (возможно, четырехугольников). Одна из точек опорной сети должна быть привязана к опорному гравиметрическому пункту I или II класса. Согласно требованиям инструкции по гравиразведке, точность определения силы тяжести на опорных пунктах должна быть в 1,5 — 2 раза выше, чем на рядовых. Повышение точности достигается главным образом за счет многократности измерений разности значений ускорения силы тяжести между опорными точками. Кратность наблюдений (количество независимых наблюдений на пунктах опорной сети) рассчитывается на основе неравенства где — среднеквадратическая погрешность единичного измерения гравиметром,

— проектируемая точность опорной сети.

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

В процессе съемки необходимо соблюдать правила эксплуатации гравиметров.

Для оценки качества рядовой съемки проводят специальные контрольные измерения, которые предусматриваются в объеме 5 —10% от числа основных наблюдений, но не менее 50 точек по всей площади измерений. Контрольные измерения должны быть равномерно распределены по площади и во времени. Они выполняются независимыми рейсами — другим оператор, другим прибором, в другой день. Если на трех и более пунктах рейса при повторном контроле получено недопустимое расхождений основных контрольных измерений, рейс (или его звено) бракуется.

Общая точность гравиметрической съемки оценивается по формуле где и — точность рядовой и опорной сети,

- среднеквадратическая погрешность определения нормального значения силы тяжести,

- среднеквадратическая погрешность, определенная поправкой Буге,

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

Наземная гравиметрическая съемка проводится гравиметрами ГР-К1, ГР-К2 (гравиметр разведочный, класс точности первый и второй), узкодиапазонными наземными гравиметрами типа ГНУ-К1, ГНУ-К2, ГНУ-КВ, С. Гравиметрический работы в горных условиях проводят широкодиапазонными гравиметрами ПИН-2. Все эти приборы относятся к оптико-механическим астазированным кварцевым гравиметрам (ГАК). Среднеквадратическая погрешность единичного измерения разности ускорений силы тяжести таких гравиметров не должна превышать ±0,02, ±0,03 и ±0,06 мГал для приборов класса А, В и С соответственно при продолжительности рейсов не менее 4 часов.

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

модель измеренного гравитационного поля (Дg)

U рег — региональная составляющая поля (поле обусловленное влиянием объектов, размеры которых соизмеримы или превышают размеры площади исследования);

Uлок — локальная составляющая поля (поле обусловленное влиянием объектов, размеры которых значительно меньше размеров площади исследования);

Uпом — помехи, обусловленные влиянием случайных факторов (погрешности измерений, внешние шумы и т. п.).

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

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

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

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

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

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

Измерения выполняются по заранее созданной системе профилей, для опорных точек которых определены пространственные координаты.

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

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

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

Метод наименьших квадратов (МНК).

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

Матрица симметрична относительно главной диагонали.

Вектор свободных членов:

Обратная матрица системы уравнений:

Произведение обратной матрицы системы на вектор свободных членов

Уравнение регрессии будет выглядеть следующим образом:

Подставляя в левую часть уравнения значения Хi и Yi, получаем теоретические значения аппроксимирующей функции U (x).

Проверка статистической значимости регрессии

Статистическая значимость уравнения определяется на основе дисперсионного анализа и метода проверки статистических гипотез Дисперсия значений U (x, y) в области изменений независимых переменных Х и Y:

где .

Dрег=3.356

Остаточная дисперсия:

Dост=0.053

Значение статистики Фишера, которая определяет превышение дисперсии регрессии над остаточной дисперсией (т.е. дисперсией случайной составляющей):

Fexp=63.246

Пороговое значение статистики Фишера: Fпор= 1.130 161 функция EXCEL Fраспобр (0.05;727;721), если Fexp больше Fпор (что мы и наблюдаем), то принимается альтернативная гипотеза, т. е. принимается гипотеза о преобладании дисперсии регрессии над остаточной дисперсией (после учета регрессии), что свидетельствует в пользу статистической значимости уравнения регрессии в целом.

Фильтрация аномальной составляющей

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

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

Математический аппарат фильтра:

— берутся значения на каждом пикете по профилю

окно фильтра

количество значений

количество значений на выходе медианный фильтр (середина чисел, расп. в порядке возрастания) Для расчетов я выбрал m=2 для того чтобы не потерять в точности полученных данных и получить большее количество значений.

результат фильтрации по каждому профилю:

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

Графики автокорреляционных функций

Для изучения корреляционных свойств сигналов используется автокорреляционная функция:

j = 0,1,2,…, m, (2)

(3)

Процедура вычисления АКФ

По графикам АКФ видим, что прослеживается закономерность в значениях аномальной составляющей.

По графикам АКФ случайной составляющей определяем, что перед нами действительно помеха, так как все графики распределены по нормальному закону с нулевым средним

Выделение объектов с различной плотностью

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

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

Примеры интерпретации:

Профиль 134

Профиль 128

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

Оценка плотности горных пород на площади исследования

Модель представлена объектами с плотностями =2.53, =2.59, =2.69, =2.73.если это класс магматических г. п., то это кислые горные породы, например, гранит; из минералов это может быть: ортоклаз (=2.53), альбит (=2.60), кварц (=2.67), кальцит (=2.70).

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

Физико-геологическая модель исследуемого участка

регрессия автокорреляционный плотность порода

Исходные данные:

Значения ?g в каждой точке измерения (пикете) по всей системе профилей


Профиль		Профиль		Профиль		Профиль		Профиль
Пикет	?g	Пикет	?g	Пикет	?g	Пикет	?g	Пикет	?g
	— 16.66		— 16.73		— 16.80		— 16.71		— 16.48
	— 16.76		— 16.62		— 16.80		— 16.65		— 16.36
	— 16.69		— 16.49		— 16.68		— 16.56		— 16.23
	— 16.65		— 16.47		— 16.51		— 16.46		— 16.20
	— 16.62		— 16.44		— 16.41		— 16.36		— 16.13
	— 16.56		— 16.39		— 16.27		— 16.21		— 15.97
	— 16.51		— 16.32		— 16.11		— 16.05		— 15.90
	— 16.40		— 16.10		— 15.91		— 15.84		— 15.81
	— 16.26		— 15.89		— 15.86		— 15.91		— 15.61
	— 16.19		— 15.83		— 15.75		— 15.84		— 15.63
	— 16.49		— 15.65		— 15.60		— 15.79		— 15.56
	— 16.27		— 15.69		— 15.50		— 15.72		— 15.52
	— 16.15		— 15.73		— 15.46		— 15.50		— 15.66
	— 15.65		— 15.74		— 15.43		— 15.45		— 15.38
	— 15.37		— 15.82		— 15.43		— 15.32		— 15.29
	— 15.25		— 15.48		— 15.51		— 15.23		— 15.13
	— 15.17		— 15.09		— 15.61		— 15.21		— 15.03
	— 14.96		— 14.84		— 15.30		— 15.10		— 14.97
	— 14.81		— 14.65		— 15.27		— 15.09		— 14.94
	— 14.70		— 14.65		— 15.14		— 14.94		— 14.84
	— 14.55		— 14.63		— 14.94		— 14.94		— 14.68
	— 14.44		— 14.54		— 14.54		— 14.93		— 14.55
	— 14.37		— 14.47		— 14.51		— 14.82		— 14.64
	— 14.32		— 14.36		— 14.39		— 14.51		— 14.56
	— 14.18		— 14.35		— 14.27		— 14.41		— 14.23
	— 14.13		— 14.19		— 14.07		— 14.23		— 14.20
	— 13.66		— 14.14		— 13.97		— 14.09		— 14.12
	— 13.47		— 13.88		— 13.84		— 14.07		— 14.10
	— 13.36		— 13.70		— 13.77		— 13.93		— 14.08
	— 13.35		— 13.60		— 14.03		— 13.88		— 13.87
	— 13.36		— 13.68		— 13.97		— 13.97		— 13.84
	— 13.30		— 13.43		— 13.62		— 14.15		— 13.83
	— 13.24		— 13.34		— 13.43		— 13.90		— 13.86
	— 13.17		— 13.21		— 13.26		— 13.66		— 14.02
	— 13.12		— 13.12		— 13.20		— 13.51		— 14.33
	— 13.18		— 13.23		— 13.03		— 13.46		— 14.09
	— 13.08		— 13.11		— 12.98		— 13.30		— 13.77
	— 12.96		— 13.09		— 12.94		— 13.20		— 13.58
	— 12.70		— 12.90		— 13.00		— 13.02		— 13.21
	— 12.79		— 12.81		— 12.87		— 12.95		— 12.79
	— 12.61		— 12.71		— 12.73		— 12.89		— 12.68
	— 12.61		— 12.56		— 12.65		— 12.92		— 12.66
	— 12.50		— 12.56		— 12.76		— 12.69		— 12.68
	— 12.51		— 12.51		— 12.65		— 12.56		— 12.68
	— 12.40		— 12.50		— 12.55		— 12.36		— 12.71
	— 12.34		— 12.46		— 12.43		— 12.34		— 12.49
	— 12.25		— 12.35		— 12.38		— 12.06		— 12.25
	— 12.06		— 12.25		— 12.25		— 12.18		— 12.32
	— 12.01		— 12.10		— 12.17		— 12.21		— 12.20
	— 12.14		— 12.13		— 12.04		— 12.21		— 12.15
	— 12.03		— 12.12		— 12.09		— 12.12		— 12.03
	— 11.99		— 12.01		— 12.11		— 12.13		— 11.98
	— 11.93		— 11.94		— 12.16		— 12.27		— 11.93
	— 11.80		— 11.87		— 11.99		— 11.94		— 11.79
	— 11.85		— 11.80		— 11.98		— 11.84		— 11.74
	— 11.69		— 11.80		— 11.88		— 11.78		— 11.61
	— 11.63		— 11.68		— 11.83		— 11.68		— 11.54
	— 11.47		— 11.74		— 11.82		— 11.70		— 11.55
	— 11.34		— 11.52		— 11.72		— 11.60		— 11.53
	— 11.25		— 11.41		— 11.60		— 11.52		— 11.52
	— 11.09		— 11.35		— 11.52		— 11.62		— 11.46
	— 11.05		— 11.22		— 11.45		— 11.57		— 11.41
	— 11.00		— 11.17		— 11.51		— 11.44		— 11.39
	— 10.89		— 11.07		— 11.37		— 11.33		— 11.24
	— 10.81		— 10.98		— 11.18		— 11.17		— 11.17
	— 10.63		— 10.82		— 11.06		— 11.12		— 11.04
	— 10.47		— 10.71		— 10.98		— 11.04		— 11.02
	— 10.36		— 10.71		— 10.90		— 10.96		— 10.77
	— 10.13		— 10.65		— 10.79		— 10.85		— 10.64
	— 9.99		— 10.42		— 10.70		— 10.65		— 10.50
	— 9.94		— 10.00		— 10.52		— 10.48		— 10.47
	— 9.97		— 10.10		— 10.08		— 10.28		— 10.16
	— 9.97		— 9.98		— 10.01		— 10.19		— 10.00
Профиль		Профиль		Профиль		Профиль		Профиль
Пикет	?g	Пикет	?g	Пикет	?g	Пикет	?g	Пикет	?g
	— 16.14		— 15.81		— 15.61		— 15.39		— 15.18
	— 16.19		— 15.83		— 15.75		— 15.53		— 15.05
	— 16.18		— 15.81		— 15.81		— 15.58		— 15.11
	— 15.92		— 15.65		— 15.86		— 15.77		— 15.18
	— 15.80		— 15.65		— 15.85		— 16.02		— 15.27
	— 15.70		— 15.64		— 15.88		— 15.99		— 15.36
	— 15.59		— 15.67		— 15.78		— 15.89		— 15.13
	— 15.55		— 15.71		— 15.54		— 15.55		— 14.73
	— 15.61		— 15.46		— 15.38		— 15.36		— 14.83
	— 15.61		— 15.32		— 15.30		— 15.34		— 15.03
	— 15.49		— 15.33		— 15.02		— 15.24		— 15.03
	— 15.38		— 15.24		— 14.96		— 15.21		— 14.98
	— 15.33		— 15.17		— 14.86		— 15.04		— 14.92
	— 15.26		— 15.03		— 14.86		— 14.86		— 14.91
	— 15.12		— 14.67		— 14.80		— 14.75		— 14.89
	— 14.87		— 14.65		— 14.73		— 14.51		— 14.78
	— 14.73		— 14.65		— 14.68		— 14.44		— 14.73
	— 14.73		— 14.60		— 14.53		— 14.25		— 14.48
	— 14.71		— 14.60		— 14.39		— 14.09		— 14.01
	— 14.62		— 14.66		— 14.08		— 13.99		— 14.05
	— 14.54		— 14.44		— 13.92		— 13.87		— 13.91
	— 14.47		— 14.02		— 13.79		— 13.85		— 13.79
	— 14.30		— 14.11		— 13.80		— 13.75		— 13.68
	— 14.31		— 14.02		— 13.61		— 13.59		— 13.45
	— 14.28		— 13.91		— 13.54		— 13.61		— 13.29
	— 14.22		— 13.95		— 13.44		— 13.43		— 13.18
	— 13.92		— 13.93		— 13.43		— 13.21		— 13.11
	— 13.84		— 13.58		— 13.37		— 13.11		— 13.08
	— 13.77		— 13.58		— 13.35		— 13.04		— 12.82
	— 13.56		— 13.44		— 13.49		— 13.07		— 12.73
	— 13.59		— 13.58		— 13.21		— 12.94		— 12.66
	— 13.49		— 13.28		— 12.92		— 12.86		— 12.71
	— 13.38		— 12.90		— 12.73		— 12.81		— 12.63
	— 13.44		— 12.88		— 12.62		— 12.69		— 12.57
	— 13.51		— 12.76		— 12.60		— 12.61		— 12.55
	— 13.48		— 12.74		— 12.54		— 12.64		— 12.40
	— 13.57		— 12.83		— 12.38		— 12.78		— 12.33
	— 13.22		— 12.60		— 12.56		— 12.65		— 12.01
	— 12.95		— 12.60		— 12.52		— 12.89		— 12.01
	— 12.62		— 12.65		— 12.17		— 12.12		— 11.87
	— 12.40		— 12.78		— 12.06		— 12.00		— 11.92
	— 12.25		— 12.74		— 11.98		— 11.87		— 11.62
	— 12.20		— 12.74		— 11.99		— 11.80		— 11.60
	— 13.13		— 12.39		— 12.05		— 11.79		— 11.42
	— 12.12		— 12.29		— 12.14		— 11.73		— 11.74
	— 12.13		— 12.21		— 12.13		— 11.59		— 11.74
	— 12.09		— 12.11		— 12.02		— 11.64		— 12.09
	— 11.96		— 12.03		— 11.85		— 11.48		— 12.16
	— 11.94		— 11.99		— 11.78		— 11.42		— 12.23
	— 11.80		— 11.95		— 11.78		— 11.52		— 12.30
	— 11.85		— 11.83		— 11.73		— 11.68		— 12.23
	— 11.71		— 11.62		— 11.60		— 11.60		— 12.15
	— 11.70		— 11.71		— 11.41		— 11.51		— 11.98
	— 11.61		— 11.50		— 11.50		— 11.29		— 11.99
	— 11.52		— 11.60		— 11.45		— 11.28		— 11.75
	— 11.46		— 11.58		— 11.35		— 11.11		— 11.59
	— 11.36		— 11.45		— 11.21		— 11.13		— 11.56
	— 11.41		— 11.15		— 11.10		— 11.15		— 11.40
	— 11.32		— 11.23		— 11.11		— 10.96		— 11.20
	— 11.35		— 11.48		— 11.00		— 10.97		— 11.24
	— 11.35		— 11.28		— 10.94		— 10.74		— 10.49
	— 11.31		— 11.21		— 10.72		— 10.51		— 10.26
	— 11.25		— 11.14		— 10.75		— 10.58		— 10.16
	— 11.09		— 11.11		— 10.63		— 10.51		— 10.25
	— 10.97		— 10.97		— 10.65		— 10.61		— 10.25
	— 10.87		— 10.75		— 10.59		— 10.60		— 10.29
	— 10.78		— 10.69		— 10.51		— 10.49		— 10.16
	— 10.61		— 10.63		— 10.40		— 10.30		— 10.16
	— 10.51		— 10.60		— 10.28		— 10.31		— 10.15
	— 10.36		— 10.48		— 10.16		— 10.31		— 10.08
	— 10.21		— 10.34		— 10.12		— 10.07		— 10.04
	— 10.10		— 9.96		— 9.89		— 10.02		— 10.02
	— 10.02				— 10.12		— 10.03		— 9.98

Пространственные координаты опорных точек системы профилей


Проф.	ПК	X	Y	Z

Значения региональной составляющей поля


пр 120	пр 122	пр 124	пр 126	пр 128	пр 130	пр 132	пр 134	пр 136	пр 138	пикет
— 15.637	— 15.918	— 16.107	— 16.267	— 16.57	— 16.729	— 16.848	— 16.91	— 16.889	— 16.774
— 15.535	— 15.811	— 15.999	— 16.159	— 16.451	— 16.615	— 16.736	— 16.796	— 16.772	— 16.657
— 15.434	— 15.706	— 15.891	— 16.053	— 16.333	— 16.501	— 16.625	— 16.683	— 16.656	— 16.54
— 15.334	— 15.601	— 15.785	— 15.948	— 16.216	— 16.388	— 16.514	— 16.571	— 16.541	— 16.425
— 15.234	— 15.498	— 15.679	— 15.843	— 16.101	— 16.277	— 16.405	— 16.46	— 16.426	— 16.311
— 15.136	— 15.395	— 15.575	— 15.74	— 15.986	— 16.166	— 16.296	— 16.349	— 16.313	— 16.197
— 15.039	— 15.294	— 15.472	— 15.637	— 15.873	— 16.057	— 16.188	— 16.24	— 16.201	— 16.085
— 14.943	— 15.193	— 15.369	— 15.535	— 15.761	— 15.948	— 16.082	— 16.132	— 16.09	— 15.973
— 14.847	— 15.093	— 15.268	— 15.435	— 15.649	— 15.841	— 15.976	— 16.024	— 15.98	— 15.863
— 14.753	— 14.995	— 15.167	— 15.335	— 15.539	— 15.734	— 15.871	— 15.918	— 15.871	— 15.753
— 14.659	— 14.897	— 15.068	— 15.236	— 15.43	— 15.629	— 15.767	— 15.812	— 15.762	— 15.644
— 14.567	— 14.8	— 14.969	— 15.138	— 15.322	— 15.524	— 15.664	— 15.708	— 15.655	— 15.536
— 14.475	— 14.704	— 14.871	— 15.041	— 15.215	— 15.42	— 15.562	— 15.604	— 15.548	— 15.429
— 14.384	— 14.609	— 14.774	— 14.945	— 15.109	— 15.317	— 15.461	— 15.501	— 15.443	— 15.323
— 14.294	— 14.515	— 14.678	— 14.849	— 15.004	— 15.215	— 15.36	— 15.399	— 15.338	— 15.218
— 14.204	— 14.422	— 14.583	— 14.754	— 14.9	— 15.114	— 15.261	— 15.297	— 15.234	— 15.113
— 14.116	— 14.329	— 14.489	— 14.66	— 14.797	— 15.014	— 15.162	— 15.197	— 15.133	— 15.009
— 14.028	— 14.238	— 14.395	— 14.567	— 14.694	— 14.914	— 15.064	— 15.097	— 15.032	— 14.906
— 13.941	— 14.147	— 14.302	— 14.475	— 14.593	— 14.816	— 14.967	— 14.998	— 14.933	— 14.804
— 13.855	— 14.057	— 14.21	— 14.383	— 14.492	— 14.718	— 14.87	— 14.9	— 14.834	— 14.702
— 13.769	— 13.968	— 14.119	— 14.293	— 14.393	— 14.621	— 14.774	— 14.802	— 14.736	— 14.602
— 13.685	— 13.879	— 14.028	— 14.202	— 14.294	— 14.525	— 14.679	— 14.706	— 14.638	— 14.503
— 13.601	— 13.791	— 13.939	— 14.113	— 14.196	— 14.429	— 14.585	— 14.61	— 14.542	— 14.405
— 13.518	— 13.704	— 13.85	— 14.024	— 14.098	— 14.334	— 14.491	— 14.514	— 14.446	— 14.307
— 13.436	— 13.618	— 13.761	— 13.936	— 14.002	— 14.24	— 14.398	— 14.42	— 14.35	— 14.21
— 13.355	— 13.532	— 13.674	— 13.849	— 13.906	— 14.147	— 14.306	— 14.326	— 14.256	— 14.114
— 13.274	— 13.447	— 13.587	— 13.762	— 13.811	— 14.054	— 14.215	— 14.232	— 14.161	— 14.019
— 13.194	— 13.363	— 13.501	— 13.675	— 13.717	— 13.962	— 14.124	— 14.14	— 14.068	— 13.924
— 13.114	— 13.279	— 13.415	— 13.589	— 13.624	— 13.871	— 14.034	— 14.048	— 13.975	— 13.83
— 13.035	— 13.196	— 13.33	— 13.503	— 13.531	— 13.78	— 13.944	— 13.956	— 13.883	— 13.737
— 12.956	— 13.114	— 13.246	— 13.418	— 13.439	— 13.69	— 13.855	— 13.865	— 13.791	— 13.644
— 12.879	— 13.032	— 13.162	— 13.334	— 13.347	— 13.601	— 13.767	— 13.775	— 13.7	— 13.551
— 12.801	— 12.951	— 13.079	— 13.25	— 13.256	— 13.512	— 13.679	— 13.685	— 13.61	— 13.46
— 12.725	— 12.87	— 12.997	— 13.167	— 13.166	— 13.424	— 13.592	— 13.596	— 13.52	— 13.368
— 12.648	— 12.79	— 12.915	— 13.085	— 13.076	— 13.336	— 13.505	— 13.508	— 13.43	— 13.278
— 12.573	— 12.71	— 12.833	— 13.003	— 12.987	— 13.249	— 13.419	— 13.42	— 13.341	— 13.188
— 12.496	— 12.632	— 12.754	— 12.921	— 12.911	— 13.162	— 13.333	— 13.332	— 13.253	— 13.098
— 12.419	— 12.554	— 12.675	— 12.84	— 12.835	— 13.076	— 13.248	— 13.245	— 13.164	— 13.009
— 12.343	— 12.477	— 12.597	— 12.759	— 12.759	— 12.99	— 13.163	— 13.159	— 13.077	— 12.921
— 12.268	— 12.4	— 12.52	— 12.679	— 12.683	— 12.905	— 13.079	— 13.073	— 12.99	— 12.832
— 12.193	— 12.324	— 12.442	— 12.6	— 12.608	— 12.82	— 12.995	— 12.987	— 12.903	— 12.745
— 12.119	— 12.248	— 12.365	— 12.52	— 12.534	— 12.736	— 12.912	— 12.902	— 12.817	— 12.658
— 12.045	— 12.173	— 12.289	— 12.441	— 12.459	— 12.652	— 12.829	— 12.817	— 12.731	— 12.571
— 11.971	— 12.098	— 12.213	— 12.363	— 12.385	— 12.569	— 12.746	— 12.733	— 12.645	— 12.484
— 11.898	— 12.023	— 12.137	— 12.285	— 12.312	— 12.486	— 12.664	— 12.649	— 12.56	— 12.398
— 11.825	— 11.949	— 12.062	— 12.207	— 12.238	— 12.403	— 12.582	— 12.565	— 12.475	— 12.313
— 11.753	— 11.875	— 11.987	— 12.13	— 12.165	— 12.321	— 12.501	— 12.482	— 12.391	— 12.228
— 11.68	— 11.802	— 11.912	— 12.053	— 12.093	— 12.239	— 12.42	— 12.399	— 12.307	— 12.143
— 11.609	— 11.729	— 11.838	— 11.977	— 12.02	— 12.157	— 12.339	— 12.317	— 12.223	— 12.058
— 11.537	— 11.656	— 11.764	— 11.9	— 11.948	— 12.076	— 12.258	— 12.235	— 12.14	— 11.974
— 11.466	— 11.584	— 11.69	— 11.824	— 11.876	— 11.995	— 12.178	— 12.153	— 12.057	— 11.89
— 11.396	— 11.512	— 11.617	— 11.749	— 11.804	— 11.914	— 12.098	— 12.071	— 11.974	— 11.807
— 11.325	— 11.44	— 11.544	— 11.673	— 11.733	— 11.834	— 12.019	— 11.99	— 11.891	— 11.723
— 11.255	— 11.368	— 11.471	— 11.598	— 11.661	— 11.754	— 11.939	— 11.909	— 11.809	— 11.64
— 11.185	— 11.297	— 11.398	— 11.523	— 11.59	— 11.674	— 11.86	— 11.828	— 11.727	— 11.557
— 11.116	— 11.226	— 11.326	— 11.449	— 11.519	— 11.594	— 11.781	— 11.747	— 11.645	— 11.475
— 11.045	— 11.156	— 11.253	— 11.374	— 11.443	— 11.514	— 11.702	— 11.667	— 11.563	— 11.392
— 10.975	— 11.085	— 11.181	— 11.3	— 11.367	— 11.435	— 11.623	— 11.588	— 11.482	— 11.31
— 10.904	— 11.015	— 11.11	— 11.226	— 11.291	— 11.356	— 11.545	— 11.508	— 11.401	— 11.228
— 10.834	— 10.945	— 11.038	— 11.152	— 11.216	— 11.277	— 11.466	— 11.429	— 11.319	— 11.146
— 10.765	— 10.875	— 10.966	— 11.078	— 11.14	— 11.198	— 11.388	— 11.349	— 11.238	— 11.065
— 10.695	— 10.805	— 10.895	— 11.004	— 11.065	— 11.119	— 11.31	— 11.27	— 11.157	— 10.983
— 10.625	— 10.735	— 10.824	— 10.931	— 10.989	— 11.04	— 11.232	— 11.191	— 11.077	— 10.902
— 10.556	— 10.666	— 10.752	— 10.858	— 10.914	— 10.962	— 11.154	— 11.112	— 10.996	— 10.821
— 10.487	— 10.597	— 10.681	— 10.784	— 10.839	— 10.883	— 11.076	— 11.033	— 10.915	— 10.739
— 10.417	— 10.528	— 10.61	— 10.711	— 10.764	— 10.805	— 10.998	— 10.955	— 10.835	— 10.658
— 10.348	— 10.458	— 10.539	— 10.638	— 10.689	— 10.726	— 10.92	— 10.876	— 10.754	— 10.577
— 10.279	— 10.389	— 10.468	— 10.565	— 10.614	— 10.648	— 10.842	— 10.797	— 10.673	— 10.496
— 10.21	— 10.32	— 10.397	— 10.492	— 10.539	— 10.57	— 10.764	— 10.718	— 10.593	— 10.415
— 10.141	— 10.252	— 10.326	— 10.419	— 10.464	— 10.491	— 10.687	— 10.64	— 10.512	— 10.334
— 10.072	— 10.183	— 10.255	— 10.346	— 10.389	— 10.413	— 10.609	— 10.561	— 10.432	— 10.254
— 10.003	— 10.114	— 10.184	— 10.273	— 10.314	— 10.334	— 10.531	— 10.482	— 10.351	— 10.173
— 9.934	— 10.045	— 10.114		— 10.239	— 10.256	— 10.453	— 10.403	— 10.271	— 10.092

Значения аномальной составляющей поля и помехи:


пр 120	пр 122	пр 124	пр 126	пр 128	пр 130	пр 132	пр 134	пр 136	пр 138	пикет
0.457	0.528	0.497	0.457	0.43	0.249	0.138	0.11	0.159	0.114
0.485	0.281	0.249	0.329	0.261	0.255	0.086	— 0.004	0.152	— 0.103
0.324	0.126	0.081	0.243	0.153	0.271	0.065	0.003	0.166	— 0.15
0.154	— 0.169	— 0.075	0.298	0.296	0.188	0.054	0.061	0.071	— 0.225
— 0.036	— 0.522	— 0.171	0.193	0.301	0.147	0.045	0.05	— 0.014	— 0.309
— 0.224	— 0.595	— 0.305	0.1	0.286	0.196	0.086	0.079	— 0.077	— 0.363
— 0.091	— 0.596	— 0.308	— 0.033	0.283	0.157	0.138	0.13	— 0.119	— 0.425
0.213	— 0.357	— 0.171	— 0.175	0.211	0.138	0.242	0.222	— 0.01	— 0.427
0.017	— 0.267	— 0.112	— 0.025	0.039	0.231	0.066	0.164	0.09	— 0.397
— 0.277	— 0.345	— 0.133	0.015	— 0.071	0.104	0.031	0.168	0.041	— 0.437
— 0.371	— 0.343	0.048	— 0.094	— 0.06	0.069	— 0.023	0.212	0.112	— 0.846
— 0.413	— 0.41	0.009	— 0.102	— 0.058	0.004	— 0.056	0.208	— 0.035	— 0.734
— 0.445	— 0.336	0.011	— 0.129	— 0.115	— 0.24	0.062	0.144	— 0.182	— 0.721
— 0.526	— 0.251	— 0.086	— 0.085	— 0.151	— 0.063	0.011	0.071	— 0.297	— 0.327
— 0.596	— 0.235	— 0.122	0.179	— 0.116	— 0.075	0.04	— 0.031	— 0.482	— 0.152
— 0.576	— 0.088	— 0.147	0.104	0.03	— 0.016	0.031	— 0.213	— 0.246	— 0.137
— 0.614	— 0.111	— 0.191	0.01	0.067	— 0.016	— 0.048	— 0.413	0.043	— 0.161
— 0.452	— 0.012	— 0.135	— 0.033	— 0.036	— 0.056	— 0.036	— 0.203	0.192	— 0.054
— 0.069	0.057	— 0.088	— 0.125	— 0.117	— 0.124	— 0.123	— 0.272	0.283	— 0.006
— 0.195	0.067	0.13	— 0.277	— 0.128	— 0.122	— 0.07	— 0.24	0.184	0.002
— 0.141	0.098	0.199	— 0.147	— 0.147	— 0.059	— 0.166	— 0.138	0.106	0.052
— 0.105	0.029	0.238	0.182	— 0.176	— 0.025	— 0.251	0.166	0.098	0.063
— 0.079	0.041	0.139	0.003	— 0.104	— 0.211	— 0.235	0.1	0.072	0.035
0.068	0.114	0.24	0.004	— 0.212	— 0.226	— 0.019	0.124	0.086	— 0.013
0.146	0.008	0.221	0.026	— 0.278	0.01	— 0.012	0.15		0.03
0.175	0.102	0.234	— 0.101	— 0.314	— 0.053	0.076	0.256	0.066	— 0.016
0.164	0.237	0.157	— 0.168	— 0.109	— 0.066	0.125	0.262	0.021	0.359
0.114	0.253	0.131	0.095	— 0.123	— 0.138	0.054	0.3	0.188	0.454
0.294	0.239	0.065	0.009	— 0.146	— 0.209	0.104	0.278	0.275	0.47
0.305	0.126	— 0.16	0.063	— 0.029	— 0.09	0.064	— 0.074	0.283	0.387
0.296	0.174	0.036	— 0.162	— 0.151	— 0.15	— 0.115	— 0.105	0.111	0.284
0.169	0.172	0.242	0.054	— 0.143	— 0.229	— 0.383	0.155	0.27	0.251
0.171	0.141	0.349	0.35	— 0.124	— 0.348	— 0.221	0.255	0.27	0.22
0.155	0.18	0.377	0.287	— 0.274	— 0.596	— 0.068	0.336	0.31	0.198
0.098	0.18	0.315	0.325	— 0.434	— 0.994	— 0.005	0.308	0.31	0.158
0.173	0.07	0.293	0.263	— 0.493	— 0.841	— 0.041	0.39	0.111	0.008
0.166	— 0.148	0.374	0.091	— 0.659	— 0.608	0.033	0.352	0.143	0.018
0.409	— 0.096	0.115	0.24	— 0.385	— 0.504	0.048	0.305	0.074	0.049
0.333	— 0.413	0.077	0.159	— 0.191	— 0.22	0.143	0.159	0.177	0.221
0.398	0.28	0.35	0.029	0.063	0.115	0.129	0.203	0.18	0.042
0.273	0.324	0.382	— 0.18	0.208	0.14	0.105	0.257	0.193	0.135
0.499	0.378	0.385	— 0.22	0.284	0.076	— 0.008	0.252	0.257	0.048
0.445	0.373	0.299	— 0.299	0.259	— 0.028	0.139	0.057	0.171	0.071
0.551	0.308	0.163	— 0.027	— 0.745	— 0.111	0.186	0.083	0.135	— 0.026
0.158	0.293	— 0.003	— 0.005	0.192	— 0.224	0.304	0.099	0.06	— 0.002
0.085	0.359	— 0.068	— 0.003	0.108	— 0.087	0.242	0.135	0.015	— 0.027
— 0.337	0.235	— 0.033	0.02	0.075	0.071	0.441	0.102	0.041	— 0.022
— 0.48	0.322	0.062	0.023	0.133	— 0.081	0.24	0.149	0.057	0.083
— 0.621	0.309	0.058	— 0.013	0.08	— 0.043	0.129	0.147	0.123	0.048
— 0.763	0.136	— 0.016	— 0.05	0.148	— 0.074	0.048	0.195	0.01	— 0.166
— 0.764	— 0.096	— 0.04	— 0.006	0.026	— 0.035	0.058	0.063	— 0.063	— 0.14
— 0.754	— 0.088	0.017	0.129	0.094	— 0.066	— 0.032	— 0.039	— 0.036	— 0.183
— 0.655	— 0.07	0.134	— 0.037	0.033	— 0.096	— 0.251	— 0.17	— 0.049	— 0.207
— 0.735	0.078	— 0.029	0.098	0.051	— 0.036	— 0.001	— 0.081	— 0.061	— 0.16
— 0.565	0.017	— 0.052	— 0.077	0.07	— 0.066	0.02	— 0.152	— 0.073	— 0.293
— 0.474	0.116	— 0.024	— 0.131	0.059	— 0.016	0.001	— 0.133	— 0.155	— 0.215
— 0.515	0.026	0.043	— 0.076	0.083	— 0.026	0.022	— 0.163	— 0.117	— 0.238
— 0.425	— 0.065	0.081	0.15	— 0.043	— 0.115	— 0.077	— 0.232	— 0.258	— 0.16
— 0.296	0.055		— 0.004	— 0.029	— 0.174	— 0.055	— 0.212	— 0.119	— 0.112
— 0.406	— 0.025	0.038	— 0.328	— 0.134	— 0.243	— 0.054	— 0.171	— 0.091	— 0.104
0.275	0.135	0.026	— 0.202	— 0.21	— 0.262	— 0.232	— 0.171	— 0.112	— 0.025
0.435	0.295	0.175	— 0.206	— 0.245	— 0.291	— 0.26	— 0.18	— 0.063	— 0.067
0.465	0.155	0.074	— 0.209	— 0.261	— 0.35	— 0.208	— 0.319	— 0.093	— 0.098
0.306	0.156	0.122	— 0.252	— 0.176	— 0.278	— 0.176	— 0.258	— 0.074	— 0.069
0.237	— 0.013	0.031	— 0.186	— 0.131	— 0.287	— 0.094	— 0.147	— 0.065	— 0.071
0.127	— 0.072	0.02	— 0.039	— 0.106	— 0.235	— 0.122	— 0.105	0.015	0.028
0.188	— 0.032	0.029	— 0.052	— 0.091	— 0.294	— 0.12	— 0.104	0.044	0.107
0.119	0.089	0.068	— 0.065	0.004	— 0.122	— 0.118	— 0.103	— 0.037	0.136
0.06	0.01	0.117	— 0.108	0.029	— 0.07	— 0.086	— 0.072	— 0.057	0.285
0.061	— 0.058	0.166	— 0.061	0.104	— 0.009	0.037	— 0.06	0.092	0.344
0.032	0.113	0.135	0.006	0.179	— 0.057	0.129	0.041	0.432	0.314
— 0.017	0.094	0.294	0.313	0.214	0.174	0.251	0.402	0.251	0.203
— 0.046	0.015	— 0.006		0.219	0.256	0.263	0.393	0.291	0.122

Полученные значения показывают, что 98% от общего значения поля занимает его региональная составляющая.

Заключение

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

Список используемой литературы

1. Серков.В.А., Курс лекций по дисциплине «Математическая обработка результатов геофизического эксперимента» — УГГУ, 2010;2011 гг.

2. Кошелев И. Н., Гравитационная и магнитная разведка. Практикум.-К.: Вища шк. Головное изд-во, 1984.

Показать весь текст

Заполнить форму текущей работой