Файл: Богданов, В. И. Вычисление гравитационных аномалий от трехмерных тел (графические способы).pdf

Положение фиксированных размеров тел в этом случае выби­ рается, исходя из реальной геологической обстановки. Отметим также, что рабочие формулы для таких производных можно свести к комбинациям уже рассчитанных коэффициентов С31 и C3 ic но­ выми коэффициентами, определяемыми выбранными фиксирован­ ными размерами цилиндрических тел.

Процесс работы с палеткой, точность вычислений

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

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

палетке число секториальных площадок в пределах каждого полукольца или кольца концентрических окружностей отдельно. Затем для каждого кольца или полукольца снимаются с кривых поправочных коэффициентов значения Сзх и Сзг. Последние сум­ мируются в пределах каждого сектора палетки и умножаются на'число секториальных площадок в секторе и на коэффициент В — отношение разности кубов косинусов смежных углов к разности

52

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

Точность вычислений по универсальной палетке оценена рас­ четами на теоретических моделях (рис. 24) и может быть повышена при использовании палеток, построенных в более крупном мас­ штабе с дополнительным расчленением секториальных площадок на более дробные, использовании промежуточных кривых С31 и С32, а также при более детальной аппроксимации реальных тел горизонтальными или вертикальными цилиндрами.

При работе с универсальными палетками рекомендуется схема вычислений, приведенная в табл. 7.

Таблица 7

Схема вычислений при работе с универсальными палетками

С ум м а

Г л н в а VI

НЕКОТОРЫЕ ВОПРОСЫ ИНТЕРПРЕТАЦИИ ГРАВИМЕТРИЧЕСКИХ МАТЕРИАЛОВ МЕТОДОМ ПОДБОРА ПЛОТНОСТНЫХ МОДЕЛЕЙ

: . ■ ГЕОЛОГИЧЕСКОГО СТРОЕНИЯ ЗЕМНОЙ КОРЫ

' *

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

53

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

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

Вопросы методики и техники интерпретации гравиметрических материалов методом подбора, редуцирования аномалий, ослабле­ ния влияний соседних масс и регионального фона, установления контуров геологических объектов, их избыточной плотности и дру­ гие подробно изложены, например, в работах [1, 8, 19, 20], а также в соответствующих учебных пособиях и инструкциях. Однако ряд частных вопросов проблемы редуцирования, разделения аномалий и определения плотности горных пород в условиях Кольского полуострова имеют свою специфику или недостаточно полно освещены в литературе.

54

Проблема редуцирования

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

Проблема выбора редукции очень сложна и для каждой кон­ кретной области гравиметрии имеет свою специфику [36, 70—75 и др. ]. В настоящее время нашли применение редукции «в свобод­ ном воздухе», Буге, Пуанкаре—Прея, топографические, изостатические, различного рода геологические, статистические и т. д. Для целей геологических исследований Е. Н. Люстих [70] пред­

рельефа, или с поправкой Вайка, учитывающей реальное распре­ деление плотностей горных пород между наинизшей и наивысшей отметками района. Обе эти редукции весьма близки между собой. Различного рода поправки, учитывающие геологическое строение района, Е. Н. Люстих рекомендует применять только для исклю­ чения гравитационных эффектов от известных масс. Весьма близка к этой задаче постановка вопроса К. Ф. Тяпкиным [19] о разделе­

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

'дукцию за притяжение плоскопараллельного промежуточного слоя и является функцией высоты точки наблюдения [36, 76]. Однако такое редуцирование связано с большими объемами вы­ числений, и поэтому на практике пользуются картами аномалий

встандартной редукции Буге или в неполной топографической редукции. При этом влияние дальних зон совместно с гравитацион­ ными эффектами, обусловленными выбором различных нормаль­ ных формул, отличием сфероида от геоида и дефектом плотности крупных водных и осадочных бассейнов, входит в так называемый «региональный фон» и ослабляется известными способами разде­ ления аномалий. Недоучет действия этих факторов может приво­ дить к появлению зональных и региональных фиктивных аномалий интенсивностью до десятка миллигал.

55

С другой стороны, поправка Буге не учитывает всех изменений силы тяжести, связанных с плотностными неоднородностями зем­ ной коры. Г. С. Бутаков [77 ], представляя вертикальный градиент силы тяжести в виде

Amax

Amin

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

СССР: И'г“ рм =(3086 + 2) Е — нормальный градиент силы тя­ жести; Т'У"р= —0.0419о=(—712-:— 1131)2? — вертикальный гради­

и локальных неоднородностей плотности земной коры. При особонеблагоприятных условиях редукция Буге может вносить погреш­ ность до 0.05 мгл/м, создавая тем самым ложные аномалии ампли­ тудой в несколько миллигал.

Приблизительно на такую же величину могут отличаться гра­ витационные аномалии, определенные на физической поверх­ ности Земли и редуцированные на плоскость. Вопросы приведе­ ния аномалий к одному уровню очень сложны и в настоящее время слабо разработаны. Некоторые способы приближенного решения их приведены в работах [8, 20, 78, 79]. При подборе плотностных разрезов больших территорий приходится, дополнительно к выше­ изложенному, проводить учет гравитационных эффектов, возни­ кающих вследствие замены поверхности уровня плоскостью. Приемы такого учета изложены в работах [43, 75 и др.].

Важно также подчеркнуть, что специфической особенностьюметода подбора является возможность выполнения интерпретации при наличии гравитационных аномалий в любой редукции, как в процессе региональных, так и детальных исследований. Необ­ ходимо только ясно представлять, каким образом осуществлялось редуцирование поля, и в зависимости от этого оперировать при подборе плотностных разрезов или избыточными плотностями (аномалии Буге), или истинными плотностями для промежуточногослоя, а для остальной части разреза — избыточными плотностями (аномалии «в свободном воздухе»). Точно так же в процессе под­ бора разрезов могут быть исключены гравитационные эффекты от рельефа дневной поверхности.

Плотность горных пород

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

56

границ геологических разрезов. В настоящее время в гравираз­ ведке хорошо разработаны методы лабораторного измерения плот­ ности и пористости сравнительно небольших по размерам (порядка

•1.0 дм3) образцов и методы статистической обработки результатов измерений [8, 38, 80]. Однако при подборе плотностных моделей геологического строения какого-либо района исследователей ин­ тересует плотность больших объемов горных пород, или, как ее называют горняки, плотность пород в массиве [81]. Последняя зависит от гетерогенности петрографического состава породы, от

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

Сведения о плотности основных структурно-фациальных ком­ плексов кристаллических пород Кольского полуострова полу­ чены в результате выполнения большого объема работ Западным геофизическим трестом, Северо-Западным территориальным гео­ логическим управлением, Ленинградским горным институтом, Кольским филиалом АН СССР и другими научными и производ­ ственными организациями. Для определения плотности отбира­ лись небольшие по размерам образцы горных пород из обнажений или керна скважин. Общее число и схема отбора образцов для какой-либо геологической разновидности пород обеспечивали представительность статистической выборки. Измерения прово­ дились путем гидростатического взвешивания на технических

•весах или денситометрах с абсолютной ошибкой определения плотности +0.02 г/см3 [82, 83].

Сведения о пористости образцов очень немногочисленны. Коэффициенты открытой пористости редко превышают 1.0% [82, 83]. Это позволяет считать, что плотность кристаллических пород Кольского полуострова вполне определяется их минерало­ гическим составом [82—84 и др. ]. Однако полевыми исследова­ ниями на Кольском геофизическом полигоне установлена повсе­ местная раздробленность и трещиноватость кристаллических пород вблизи дневной поверхности [85—88 и др.]. Это — расколы, со­ путствующие региональным и глубинным разломам в земной коре, трещиноватость средних размеров, фиксируемая при дешифриро­ вании аэрофотоснимков, трещины отдельности и тектонической трещиноватости, расчленяющие породы на все более и более мел­ кие ячейки.

57