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

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

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

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

Рис. 2.34. Положение абсцисс характерных точек графиков, используемых при количественной интерпретации

Однако чаще при интерпретации графиков наблюдаемого гравитационного поля используют так называемый метод характерных точек. Правильнее было бы называть его "метод абсцисс характерных точек графиков поля". В основе этого метода лежит связь между численными значениями какого-либо параметра объекта (например, глубины его залегания) и абсциссами наиболее характерных точек графика поля: полумаксимума D g или V _ZZ, перехода через нулевую ординату, экстремумов, точек равных значений V_ZZ и V _ZX и т.д.

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

Dg_max =k*M / h²; (2.81)

Отсюда избыточная масса будет равна M=_Dg_max*h² / k;

Объем шара V= и радиус его эпицентрального сечения R= . В итоге полная масса шара будет равна P=V*s , где _Ds - избыточная плотность объекта, s - истинная плотность вещества шара (руды), Р- запасы руды.

Характерные точки гравитационного поля простейших геометрических тел

(для эпицентрального профиля)

Таблица 2.4.

Примечание: h — глубина залегания центра масс;

k — гравитационная постоянная, равная 6,67*10^-11 м³/кг*с²;

М=V*_Ds =4/3*pR³*_Ds (избыточная масса шара);

m=S*_Ds = pR² *_Ds (избыточная масса единицы длины цилиндра);

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

(2.82)

Решение обратной задачи (т.е. определение параметров объекта) обычно выполняют в следующем порядке:

1. Построив графики V_Z,V_ZZ, V_ZX, по характерным точкам (рис.2.34) определяют глубину до центра объекта (h).

2. Формулы (2.75) или (2.78) записывают для максимального значения V_Z (V_Zmax =V_Z при x=0).

3. Сняв с графика значения V_Zmax и, зная величину h, выражают из полученной формулы избыточную массу М и вычисляют ее.

4. Учитывая, что избыточная масса для сферы определяется, как М=(s_об - s_вм)*V, а для горизонтального цилиндра, как М=(s_об - s_вм)*S, можно определить объем сферы (V) или площадь сечения (S) цилиндра.

5. Зная V и S , можно по известным формулам объема сферы и площади сечения определить радиус объекта, а затем минимальную глубину бурения (т.е. глубину до верхней точки объекта).

6. Запасы изометричной залежи руды (Р) определяются, как Р=V*s_об, а для цилиндра можно определить запасы единицы его длины, как Р =S*s_об.

7. Расчет относительных погрешностей ведется по формулам:

e_h= (çh_p - h_ист.ç/h_ист.)*100% ,

e_р = (çP_p - P_ист.ç/P_ист.)*100% ,

где e_h; e_p - относительные погрешности определения глубины и запасов объекта соответственно,

h_Р, h_ист., P_Р, P_ист. - рассчитанные и истинные значения глубин и запасов.