Количественная интерпретация

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

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

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

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

Косвенные методы, применяемые для обработки как простых, так и сложных аномалий, включают методы, основанные на применении ЭВМ, а также палеточные, графические и др.

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

На наблюденной или трансформированной карте выделяется отдельная (региональная или локальная) аномалия. Строго через ее центр перпендикулярно простиранию изолиний строится график ∆g_Б. Иногда вкрест полученных аномалий проводятся более точные и детальные полевые работы для получения интерпретационных графиков.

Если на карте имеются изометрические аномалии (длина и ширина отличаются не более чем в 2-3 раза), то, исходя из общего геологического строения района, их аппроксимируют шаром или вертикальным столбом.

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

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

Если на карте ∆g наблюдаются вытянутые аномалии (длина более чем в 3-6 раз превышает ширину), то создающие их геологические объекты могут быть аппроксимированы бесконечно длинным цилиндром или пластом. Это могут быть вытянутые, горизонтально (или полого) расположенные геологические объекты, поперечное сечение которых мало меняется (антиклинальные и синклинальные складки, линзообразные залежи полезных ископаемых, пласты и т.п.).

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

Для перечисленных моделей простой геометрической формы количественная интерпретация методом характерных точек дает следующие результаты. Центр тела, создавшего аномалию, располагается под максимумом ∆gmax за исключением уступа, который располагается там, где аномалия равна половине максимума. Абсцисса максимума (или 0,5∆g_max для уступа) принимается за начало координат, а слева и справа от нее находятся абсциссы x_1/2 точек, для которых ∆g равно 0,5∆g_max (для уступа 0,25∆g_max). Зная величины x_1/2 можно определить глубину залегания или центра возмущающего тела (для шара, цилиндра), или верхней кромки (для столба), или середины высоты уступа (для сброса) с помощью таблицы 1.

Таблица1

Аппроксимируемый объект	шар	цилиндр	столб	уступ
Глубина залегания

Определив h, можно найти избыточную массу M=V(σ-σ₀). Зная из посторонних источников (например, по данным измерений плотности образцов пород, взятых из керна скважин или обнажений) избыточную плотность, можно по избыточной массе рассчитать поперечные размеры разведываемых тел. Большинство отдельных аномалий, наблюдаемых на гравиметрических картах, с той или иной степенью приближения может быть отнесено к рассмотренным выше типичным аномалиям.

Результаты количественной интерпретации этим методом будут близки к истинным (с точностью до 10 - 20%), если имеется площадное распределение аномалий (получены карты ∆g), возмущающие массы близки по форме к телам перечисленной простой геометрической формы, расстояния между отдельными геологическими неоднородностями превышают их размеры и известна избыточная плотность. Однако если эти условия не выполняются, то методы характерных точек позволяют оценить глубины и полную аномальную массу лишь приближенно (с погрешностью до 30% и больше).

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

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

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

где a- коэффициент, меняющийся от 0,7 для вытянутых до 0,9 для изометрических тел; ∆g_max и ∆g’_max - максимальная аномалия силы тяжести и максимальный горизонтальный градиент (∆g’_max=∆g_max/∆x) снятые по графику ∆g, проходящему вкрест простирания структур.

Дата добавления: 2015-08-31