Курсовая работа по гравиразведке и магниторазведке

«Геологическая задача и геологическая обстановка»
В курсовой работе условно ставится следующая геологическая задача: требуется комплексом методов магнитометрии и гравиметрии для целей детальной геологической разведки определить форму, размеры, глубину залегания и физические свойства тел геологических объектов, сложенных сильномагнитными плотными рудами. Тела расположены на достаточном удалении друг от друга так, что взаимное влияние одной аномалии на другую не наблюдается. Рудные залежи могут быть аппроксимированы телами правильной геометрической формы.
1. Вычисление гравитационных аномалий от тел правильной формы (решение прямой задачи для цилиндра)

Вычисление аномалий от тел с заданными параметрами необходимо для выбора точности наблюдений и густоты съемочной сети. При решении прямой задачи по гравиразведке избыточная плотность не должна превышать 1 г/см3, а размеры, глубина залегания, плотность должны быть такими, чтобы гравитационная аномалия на максимуме составляла не менее 0,5 мГал. Вычисления производятся по известной из теоретического курса формуле с точностью 0,5%, приведенной ниже. Решение прямой задачи осуществляется в среде Mathcad.

Формула для решения прямой задачи для цилиндра:

Рассчитанное аномальное поле изображается в виде графика вдоль профиля наблюдения (рис.1).

Рис.1. Гравитационное поле цилиндра


2. Определение точности наблюдений и густоты съёмочной сети


Точность наблюдений
Точность наблюдений обычно характеризуется величиной средней квадратической ошибки ε. Соотношение между средней квадратической и средней арифметической σср ошибками имеет вид: σср=0,8٠ε.
По величине средней квадратической ошибки гравиметрические съемки соответственно подразделяются на виды:
– съемки пониженной точности (ε>0,5 мГал);
– съемки средней точности (0,5<ε<15 мГал);
– съемки высокой точности (ε<0,1 мГал).
Ошибки наблюдений сказываются на результатах интерпретации аномалий, так как в расчетные формулы для определения элементов зале-гания и других характеристик геологических объектов в качестве исходных данных входят абсциссы и ординаты характерных точек аномальных графиков. Чем больше погрешности съемок, тем сильнее искажены ре- зультаты интерпретации.
Из выше изложенного следует, что при выборе точности съемки необходимо исходить из требуемой точности определения тех или иных параметров возмущающих объектов. На практике ошибки интерпретации в пределах 5% считаются приемлемыми.
Для определения точности съемки построим часть кривой Δg, включающую точку (Δg 1/2; х1/2), где Δg 1/2=0,5Δg max, в более крупном масштабе (рис. 2). Загрубим значение Δg 1/2 на 5%: Δg '1/2= Δg 1/2+0,05*Δg 1/2. Отложим точку (Δg '1/2;х1/2) и проведем через нее прямую параллельную части теоретической кривой. Найдем абсциссу точки пересечения загрубленной части графика с прямой у= Δg 1/2, обозначим ее х'1/2. Определим приращение Δабс=х'1/2 – х1/2 и найдем:







Рис.2. Определение точности наблюдений


Определение шага наблюдения
В основу выбора расстояния меду точками наблюдений положен принцип, согласно которому все части графика аномального поля, как в точках наблюдения, так и между ними, должны быть равноточными. Поле между точками наблюдений содержит ошибки за нелинейность. Эти ошибки называются ошибками интерполяции и обозначаются δинт. Расстояние между точками наблюдений должно быть таким, чтобы δинт<ε.
Ошибка за нелинейность поля в линейной части аномалии стремится к нулю, а в зоне экстремумов достигает максимального значения. В экстремальной части аномалии по теоретическому графику необходимо изучить характер изменения δинт в зависимости от шага наблюдений Δхi. Суть анализа показана на рисунке 3 и 4.



Рис.3 Определение шага наблюдений

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

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

Рис.4. Определение Δхрабочее