Лабораторная работа № 8 решение прямой и обратной задач гравиразведки

Цель работы: знакомство с принципами решения прямых и обратных задач гравиразведки при количественной интерпретации гравиметрических данных для тел правильной геометрической формы.

1. Основные теоретические сведения.

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

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

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

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

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

1.1. Способы решения прямой задачи.

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

Неправильные и сложные объекты могут аппроксимироваться набором правильных тел. Это удобно тем, что для правильных тел аналитические выражения поля силы тяжести и его градиентов известны. Например, выражения поля силы тяжести _g и его градиентов V_ZX и V_ZZ для сферы имеют вид (при у=0):

V_Z =_g=; (1)

V_ZZ=; (2)

V_ZX=; (3)

Для горизонтального кругового цилиндра:

V_Z =_g=; (4)

V_ZZ=; (5)

V_ZX=. (6)

где k - гравитационная постоянная, равная 6,673 10^-8 см³/г с²;

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

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

x - текущая координата; h - глубина центра сферы;

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