Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Вольфганг Торге - Гравиметрия - 1999.pdf
Скачиваний:
244
Добавлен:
06.06.2015
Размер:
30.51 Mб
Скачать

Гравиметрические съемки

297

Предложение Леваллуа в 1971 г. [410] о создании мировой абсолютной гравиметриче­ ской сети было развито в нескольких резолюциях МАГ. Учитывая геодинамические явле­

ния в разных регионах земного шара, специальная исследовательская группа 3.87 МАГ

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

дополнительные пункты [60]. В настояшее время начинаются работы по созданию такой

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

блюдений (разд. 5.2.2).

9.2. Гравиметрические сети

9.2.1.Назначение, выбор nунктов

Гравиметрические сети создаются для образования систем опорных пунктов; эти пункты должны быть закреплены центрами. Значения силы тяжести на опорных пунктах определяют из измерений (разд. 9.2.2) с последующей их обработкой (разд. 9.2.4). Координаты пунктов, значения силы тяжести и другие необходимые

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

гравиметрических данных (разд. 9.4.1, 9.4.2).

Опорные пункты служат основой для детальных гравиметрических съемок,

как профильных, так и площадных; повторные измерения на этих пунктах ис­

пользуются для мониторинга временных вариаций силы тяжести. Различают сле­

дующие типы сетей:

-мировые гравиметрические сети с расстояниями между пунктами от несколь­

ких сотен до нескольких тысяч километров. Они являются главными элемен­

тами референцных гравиметрических систем и создаются по международным программам (разд. 9.1);

-региональные сети с расстояниями между пунктами от нескольких километров

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

(разд. 9.2.6);

-локальные гравиметрические сети с расстояниями от нескольких сотен метров

до десятков километров. Их создают преимущественно для целей геофизики и геодинамики (разд. 9.3.1, 10.1.1).

По возможности должна быть обеспечена долговременная сохранность опор­

ных пунктов мировых, региональных и геодинамических сетей (несколько десяти­

летий и более), а при повторных измерениях условия наблюдений должны

обеспечивать достижение точности не хуже 0,01 - 0,1 мкм·с- 2

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

-равномерное, насколько возможно, расположение опорных пунктов на терри­

тории работ, исключением могут быть локальные геодинамические сети;

-мировые и региональные сети должны содержать определенное число пунктов

для непрерывных и повторных наблюдений;

-геологическая, сейсмическая и гидрологическая стабильность в окрестностях

пунктов;

298

Глава 9

-стабильное основание для установки прибора (наземные основания в зданиях,

пилоны, скальные основания, бетонный пол), закрепление пункта; выгодно

использовать существующие центры плановых или высоких геодезических

сетей;

-определение 2 -3 лунктов-спутников (±0,01 - 0,1 мкм·с- 2) для контроля

неизменности основного пункта относительно фундаментальной гравиметриче­

ской сети (расстояние до спутников от нескольких сотен метров до нескольких

километров, дg < 100 мкм·с- 2);

-определение плановых координат и высот относительно пунктов националь­

ных геодезических сетей;

-плановая ( ± 1 см) и высотная ( ± 1 мм) привязка к ближайшим геодезическим

пунктам.

Для пунктов, на которых предстоят наблЮдения с абсолютными гравиметра­

ми, должны соблюдаться дополнительные условия (внутреннее устройство зда­

ния, подвод электропитания и др.). В сетях СГУШения и при пониженной точности работ требования к выбору опорных пунктов могут быть ослаблены. Однако и в таких сетях необходимо свести к минимуму влияние близких притягивающих масс, чтобы полученные аномалии силы тяжести были представительными для большой территории (разд. 9.3.1).

9.2.2.Измерения

Гравиметрические сети создают как абсолютными, так и относительными грави­

метрами.

По затратам средств и времени на определение одного пункта с одинаковой точнос­

тью (имея в виду необходимое число повторных измерений) в настоящее время (1987 r.) абсолютные и относительные измерения соотносятся как l О: l. В работах по абсолютным определениям используют около lO транспортабельных приборов, тогда как число приме­

няемых пружинных гравиметров составляет несколько тысяч.

Гравиметры до начала измерений в сети должны пройти исследования и кали­

бровку. К ним относятся:

-лабораторные исследования для определения зависимости показаний прибора

от атмосферного давления, температуры, магнитного поля, сотрясений и дру­

гих факторов;

- изучение характера смещения нуль-пункта пружинных гравиметров

(разд. 6.6.3);

- калибровка стандартов длины и времени у абсолютных гравиметров

(разд. 5.1.4);

-калибровка пружинных гравиметров на калибровочной системе, диапазон

которой превышает изменение силы тяжести на территории съемки

(разд. 6.4.5);

-сравнение различных приборов из измерений на одних и тех же пунктах или

из измерений одних и тех же приращений силы тяжести.

Такие исследования позволяют:

-оценить точность, достижимую с данным прибором;

Гравиметрические съемки

299

-выявить источники систематических ошибок, которые необходимо учитывать

поправками или инструментальными мерами;

-сделать заключение о частоте повторных измерений в сети.

Мировые и региональные фундаментальные гравиметрические сети, а также

геодинамические сети должны иметь высокую точность. Для таких сетей нужны

избыточные данные из повторных наблюдений при различных условиях (разные

гравиметры, различный характер внешних влияний). При этом

-выявляются грубые и систематические ошибки;

-влияние остаточных систематических ошибок приобретает случайный

характер;

-оценка точности дает величины ошибок, близкие к реальным.

Теория оптимизации помогает составить рациональный график построения се­ тей (разд. 9.2.5).

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

пример, гравиметры Ла Коста - Ромберга), при этом рекомендуется немедленно nосле измерения на определяемом пункте возвращаться на исходный пункт (раз­

ностный способ); при транспортировке по воздуху удобен метод звезды (разд. 6.6.3). Таким образом создают сети, состоящие из треугольников и четы­ рехугольников. Абсолютные измерения и связи с перекрытиями стабилизируют сеть, особенно на краях.

Сочетание абсолютных и относительных измерений позволяет калибровать гравимет­

ры и осуществлять взаимный контроль приборов. Более того, такое сочетание позволит

уточнить абсолютный масштаб (уровень) силы тяжести в сети. Например, если 4 пункта

абсолютных определений ( ± 0,1 мкм·с- 2 , систематические ошибки отсутствуют) связаны

«безошибочно» измеренными приращениями силы тяжести, то ошибка масштаба силы

тяжести уменьшается и становится равной ±0,1/д = ±0,05 мкм·с- 2

Измерения в сетях более низких классов выполняют относительными грави­

метрами, а nункты этих сетей часто располагают между пунктами фундамен­

тальной сети в виде ходов, измерения выполltяют методом профилей или

ступенчатым методом (разд. 6.6.3). Контроль смещения нуль-пункта следует осу­ ществлять через несколько часов. Примеры фундаментальной сети и сетей сгуще­ ния можно найти в разд. 9.2.6 и 9.3.6.

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

позволят контролировать изменения силы тяжести за длительные промежутки

времени (разд. 10.1.1).

Возможная точность измерений сильно зависит от внешних условий. Неблагоприят­

ные условия (очень высокая или низкая температура, быстрые ее изменения, ветер, силь­

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

300

Глава 9

Измерения силы тяжести частично автоматизируют благодаря цифровой ре­

гистрации и предварительной обработке данных. Обработав полученные на

пункте результаты (вычисление среднего и т.д.), можно оценить полученную точ­

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

кие регистрирующие устройства лишь в редких случаях.

В Национальном географическом институте, Брюссель [539], была разработана и ис­

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

Коста - Ромберга. Она состоит из цифрового блока регистрации и обработки, установ­ ленного на общем основании с гравиметром. Кроме того, имеются блок регистрации вре­

мени, источник питания и цифровой кассетный накопитель. Блок регистрации и обработки

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

туру, номер пункта) и контролирует источник питания. В полевых условиях при

осреднении данных, полученных в !О-секундных интервалах, ошибка составила

± 5 нм· с- 2 Модификацию такой системы можно использовать для регистрации земных

нриливов (разд. 10.1.4). Представителем кварцевых пружинных гравиметров с микропро­

цессором является гравиметр Сиитреке CG-3 (разд. 6.5.3).

9.2.3.Приведения отсчетов и поправки

Прежде чем приступить к обработке сети, в отсчеты по гравиметру следует вве­

сти различные поправки. Таким образом осуществляют по крайней мере прибли­

женный перевод отсчетов в единицы ускорения силы тяжести, учитывают

систематические инструментальные влияния, а также изменения силы тяжести во

времени и, наконец, относят данные к исходному пункту.

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

учитывает данные калибровки. Отсчеты по относительному гравиметру преобра­

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

ные измерения учитывают в настоящее время конечность скорости света

(разд. 5.1.6). Для гравиметров с металлической пружиной можно построить зави­

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

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

как на относительные измерения влияют лишь короткопериодические вариации,

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

Поправки за приливные изменения силы тяжести можно определить в зависи­ мости от нужной точности различными способами (разд. 3.4.2, 3.4.3):

- вычисление по разложению в ряд приливного потенциала для жесткой Земли по Картрайту-Тейлеру-Эддену (рекомендация МАГ N2 11, 1971 г.), причем гло-

Гравиметрические съемки

301

бальный гравиметрический фактор о = 1, 16, а фазовая задержка равна нулю. Ря­

ды с ограниченным числом членов (например, ежегодно публикуемые таблицы

в «Geophysical Prospecting» с о = 1,2) могут дать ошибки до О,1 мкм·с- 2 ;

- вычислецие по разложению в ряд потенциала жесткой ЗемЛИ! с амплитудными факторами для упругой Земли (модель 1066 Гильберта и Дзевонского [219]) по

Вару [745] и с учетом влияния океанских приливов на амплитуды и фазы главных

волн [619] в соответствии с резолюцией N2 9 МАГ, 1983 г. [552];

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

или интерполяция параметров земных приливов по измерениям на пунктах, рас­

положенных вблизи района работ (разд. 10.3.3).

Два последних метода дают поправки за земные приливы с ошибкой не хуже

нескольких сотых долей мкм ·с- 2, время измерений необходимо фиксировать с

точностью :2:: 1 мин.

Гравиметрическая поправка за движения полюса имеет долгопериодический характер (разд. 3.4.1). Необходимые мгновенные координаты полюсах и у (со­ ставляющие отклонения мгновенного полюса от МУН в направлении меридианов

Л = оо и Л = 270°) можно экстраполировать в пределаХ одного года по данным,

публикуемым МБВ (±0,05") [633]. Преобразовав (3.18), получим поправку в

виде

ognoл = ОПОЛU'2R sin 2'Р(Х cos л - у sin Л).

(9.2)

Для амплитудного фактора опол можно принять значения 1,О или 1,2 (прилив­

ный фактор) [746]; имеются результаты первых определений этой величины

(разд. 10.3.2).

Прямое и косвенное влияние вариаций атмосферного давления вызывает из­

менение силы тяжести на - 3 - -4 нм·с- 2/ГПа (разд. 10.2.1). Если параметры

регрессионной зависимости для данной местности неизвестны, то поправку сле­

дует вычислять с коэффициентом - 3 нм ·с- 2/ГПа (резолюция МАГ N2 9,

1983 г.). Такая поправка необходима, чтобы учесть отличие реального давления воздуха р (измерения барометром с ощибкой ± 1 ГПа) от нормального давления

Рн на пункте; ее вычисляют по формуле

(9.3а)

Величину Рн находят по высоте пункта Н для стандартной атмосферы (модель DIN 5450, которая хорошо согласуется со стандартной атмосферой США 1962 г.)

по формуле [61]-

Рн

= 1013 25 ( 1 -

00065 н(м) )

5,2559 гпа_

(9-3б)

'

288,15

 

 

Атмосферное давление обычно меняется не более чем на

±50 ГПа, и влияние

его изменений на абсолютную величину g достигает О,15 мкм ·с- 2 _ При измере­

нии приращений силы тяжести вариации давления в основном остаются меньше

10 ГПа, что соответствует ,: 0,03 мкм ·с-2

Изменения уровня грунтовых вод и влажности почвы могут вызвать измене­

ния силы тяжести до нескольких десятых долей мкм ·с- 2 При создании фунда­

ментальных гравиметрических сетей и геодинамических сетей рекомендуется

302

Глава 9

приводить измеренную силу тяжести к нормальным для данного пункта услови­

ям. Разработка способов приведения пока еще в начальной стадии (разд. 10.2.2). Для приведения измерений к единой точке (марке пункта) необходима поправ­

ка за смещение прибора в плане и по высоте; точность поправки должна состав­

лять ::t:O,Ol МКМ·с- 2

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

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

ниях относительными гравиметрами прибор нужно всегда устанавливать на одну и ту же

высоту (3 мм ~ 0,01 мкм·с- 2). Если считать, что вертикальный градиент силы тяжести

на пунктах одинаков, редуцирование за высоту выполнять не нужно. При неодинаковых

высотах гравиметра, как правило, достаточна редукция с нормальным градиентом

- 3086 им·с- 2/м. При больших различиях высоты прибора и/или значительной аномалии

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

Если на пункте одновременно работают абсолютный и относительный гравиметры либо

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

ние силы тяжести в плане [567] либо определять его по локальной модели градиента [173].

9.2.4.Обработка результатов

Величины силы тяжести и другие неизвестные находят из уравнивания сети с из­ быточными измерениями (число измерений больше числа неизвестных). Уравни­

вание выполняют по методу наименьших квадратов [52, 305]. При этом

реальный физический объект аппроксимируется математической моделью, со­

держащей функциональную и случайную части. В большинстве случаев при урав­ нивании применяют параметрический способ [530], который будет рассмотрен

ниже.

Функциональная модель описывает связь между результатами измерений и не­

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

1 + v =Ах.

(9.4)

Здесь 1 - вектор n измерений, а х - вектор и неизвестных, которые необходимо

оценить (п > и). Матрица А размерностьЮ n х и, состоящая из коэффициентов

модели, связывает результаты измерений с неизвестными. Вектор поправок v

определяется влиянием неучтенных ошибок измерений. Как обычно, при обра­

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

нейной функцией. Если зависимость нелинейная, ее линеаризуют разложением в ряд Тейлора, используя при этом приближенные величины.

Стохастическая модель описывает характер ошибок измерений через ковари­

ационную матрицу

 

ат

QI2a1a2

Q1na1an

 

 

 

r;/1 =

Qz1aza1

а~

 

(9.5)

 

 

 

 

 

Qп1апа1

 

а~

 

 

 

 

 

 

 

Гравиметрические съемки

303

где и2 - дисперсия, а е- коэффициенты корреляции. В общем случае нормиро­

вание матрицы Eu выполняют после введения весов измерений

 

 

 

 

об

(9.6)

 

 

 

р; = of'

где обдисперсия единицы веса (р =

1). Получим матрицу обратных весов (ко­

факторную

матрицу)

 

 

 

 

 

 

Qu =

-1

Eu

(9.7а)

 

 

 

 

об

 

 

и весовую

матрицу

 

 

 

 

 

 

 

Pu = Qii 1

(9.7б)

Из основного условия уравнивания

 

 

 

 

 

 

Pv =>

min

(9.8)

следует алгоритм уравнивания. Из системы нормальных уравнений

 

 

 

 

Nx- Ат Pl =О,

N =Ат РА

(9.9)

получим

величины

неизвестных

 

 

 

 

 

 

х = N-IAT Pl,

(9.10)

матрицу обратных

весов уравненных

значений неизвестных

 

 

 

 

Qxx = N- 1 = (АтРА)- 1

(9.11)

и вектор

поправок

 

 

 

 

 

 

 

v =Ах- 1.

(9.12)

Оценкой величины u0 служит эмпирическое среднеквадратическое отклонение еди­ ницы веса; оно определяется формулой

So = J-~TpV.

(9.13)

п-и

 

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

 

Ехх = .roQxx.

(9.14)

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

можно, используя закон переноса ошибок.

Для линейных функций

'Р = fтх

(9.15)

дисперсия вычисляется по формуле

 

s; = fтExxf.

(9.16)

Уравнивание по методу нанменьших квадратов устраняет невязки в сети и дает на­

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

304

Глава 9

В гравиметрических сетях распределение ошибок близко к нормальному, если измерения

выполняются выверенными гравиметрами и в максимально . различных условиях

(разд. 9.2.2).

Установим математическую модель измерений силы тяжести [698).

При абсолютных определениях время является непосредственно измеряемой

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

ние для каждого пункта используется как квазиизмерение /. Уравнение связи име-

ет ВИД

+ v = g.

(9.17)

1

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

эти корреляции, то полагают, что ковариационная матрица ошибок (9.5) абсо­ лютных измерений диагональная, и на ее главной диагонали располагаются дис­

персии, значения которых определены априори.

При относительных измерениях силы тяжести непосредственно измеренными

величинами являются «деления шкалы счетчика» z. После предварительной кали­

бровки, редуцирования (разд. 9.2.3), определения калибровочной функции (6.52) и смещения нуль-пункта (6.75) уравнение связи для исправленного отсчета имеет

ВИД

 

1 + v = g - No - д.F(z) + D(t).

(9.18)

Помимо значения силы тяжести в (9.18) входят дополнительные параметры, та­ кие, как неизвестный исходный уровень N 0 , калибровочные параметры (корректи­

рующая функция калибровки д.F(z)) и параметры смешения нуль-пункта (функция

дрейфа D(t)). Если поправки за калибровку были определены с помощью кали­

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

точным (разд. 6.4.5).

Смещение нуль-пункта в большинстве случаев учитывают введением неизвестного об­ щего уровня для каждого прибора и на каждый день измерений (разд. 6.6.3). Если измере­ ния в сети выполняются длительное время, модель (9.18) будет содержать большое число неизвестных, определяющих общий уровень и смещение нуль-пункта; это приводит к уве­ личению объема уравнительных вычислений. Объем вычислений можно уменьшить, если образовать квазиизмеренные величины как разности последовательных отсчетов.

Уравнение связи для разности отсчетов д.l;,J = 11 -

l;

имеет

вид

bl;,1 + v;, 1 = g1 - g; - (д.F(z1) - д.F(z;)) + (D(t1)

-

D(t;)).

(9.19)

Неизвестные, определяющие общий начальный уровень, исключаются. Посколь­ ку разность моментов измерений l; и /1 мала, достаточно описывать дрейф поли­

номами низких степеней.

Существуют многочисленные подходы к исследованию точности относитель­

ных измерений по результатам полевых и лабораторных наблюдений. Неполный

Гравиметрические съемки

305

учет смещения нуль-пункта и влияние других остаточных систематических оши­

бок приводят к возникновению физической корреляции результатов наблюдений.

При увеличении времени транспортировки среднеквадратическая ошибка разно­ стей отсчетов увеличивается (разд. 6.6.4). Для сетей с большим числом избыточ­

ных измерений эти корреляции можно оценить по поправкам, полученным после

предварительного уравнивания; последние используют для вычисления эмпириче­

ских коэффициентов корреляции, их описывают зависящей от времени ковариаци­ онной функцией [33].

Для приращений силы тяжести, полученных в Северной Исландии из следующих друг за другом измерений с гравиметрами Ла Коста - Ромберга по последовательной схеме (время транспортировки 5 - 15 мин), определен средний коэффициент корреляции, рав­ ный г = -0,30 ± 0,04 [349]; для D.g, полученных из наблюдений по непоследовательным схемам, коэффициент уменьшается до г< 0,1. Необходимо отметить, что·при использова­

нии формулы (9.19) смежные D.g с идентичным отсчетом на общем пункте алгебраически

коrрелированы с коэффициентом г = - 0,5. Игнорирование этого обстоятельства приво­

дит .к искажению значений силы тяжести на величину порядка среднеквадратической

ошибки; оцененные точности оказываются завыщенными [153]. Зачастую не имеет смысла

пр~1нимать во внимание эту алгебраическую корреляцию, поскольку на нее накладывается более сильная физическая корреляция, которую обычно определить точно невозможно.

В общем случае в уравнивание сети входят абсолютные и относительные из­ мерения силы тяжести. Роль измеренных абсолютных значений g могут играть величины, полученные из уравнивания сети более высокого класса; их включают в уравнивание как квазиизмерения со своей кофакторной матрицей (9 .11) или по крайней мере с их среднеквадратическими ошибками. Это случай иерархического пос-троения национальных гравиметрических сетей (разд. 9.2.6). Если имеются лишь абсолютные измерения, то в уравнивании нет необходимости. Выполнить

уравнивание сети, в которой сделаны лишь относительные измерения силы тяже­

сти, можно, зафиксировав одно значение g (например, принимая значение g на одном из пунктов как безошибочное); можно 1 акже уравнять такую сеть как сво­ бодную [461]. В последнем случае исходный уровень определяется по аппрокси­

мации значений силы тяжести на всех пунктах сети.

На рис. 9.2 приведены измеренные приращения силы тяжести (средние значе­

ния), невязки полигонов и уравненные приращения дg в локальной опорной сети,

измеренной с гравиметрами Ла Коста - Ромберга. Исходный уровень силы тя­

жести в этой сети фиксируется ее величиной на пункте Акюрейри, связанном с

МГСС-71.

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

ших и мини-ЭВМ. Приведем две разработки Института геодезии Ганнаверского

университета.

Пакет программ GNLSA предназначен для большого компьютера CDC Cyber 76/18 [764]. Программа предварительной обработки GNPREN выполняет пересчет отсчетов по

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

(рис. 9.3). Исправленные отсчеты по гравиметрам записывают на магнитный диск и вво­

дят в главную программу GNLSA, выполняющую совместное уравнивание абсолютных

и относительных измерений в сети (рис. 9.4). Процедура уравнивания основана на пара­ метрическом способе, в котором в качестве измеренных величин принимают разности от-

306

Глава 9

Рис. 9.2. Местная гравиметрическая опорная сеть (масштаб не сохранен) в Северной Исландии, в скоб·

ках указано число связей, верхнее число - среднее значение измеренной величины, нижнее

число - ее уравненное значение (мкм·с- 2 ), невязки даны в единицах 0,01 мкм·с- 2 (из базы

данных Института геодезии, Ганновер).

счетов (9.19). Система нормальных уравнений (9.9) строится последовательно, по мере поступления данных. После обращения матрицы коэффициентов нормальных уравнений

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

также поправки из уравнивания (9.10), (9.12); кофакторную матрицу (9.11) можно записать

впамять. Эти поправки можно использовать для статистического анализа (разд. 9.2.5).

Число относительных измерений практически неограниченно, количество абсолютных зна­ чений не должно превышать 31, а максимальное число относительных гравиметров - 40.

Нанбольшее число неизвестных равно 950.

Пакет Программ GRAVIТY был разработан для микрокомпьютера типа НР 9816.

Программы составлены на языке HP-Basic 3.0 и обеспечивают обработку и оценку данных

вквазиреальном масштабе времени. Поэтому качество измерений в сети можно оценивать прямо в ходе полевых работ (рис. 9.5). Входными данными служат файлы результатов измерений (MDAT), параметров калибровки (САЦ и параметров пункта (STAT). Програм­ ма DATAIN служит для формирования и преобразования файлов данных. Программа

Калибровочные таблиць1

изготовителя

Разложение приnивного

nотенциала

Координаты nункта,

параметрылриливов

Исправленные отсчеты

по гравиметру

Рис. 9.3. Блок-схема программы вычисления и учета поправок в грав••метричссю•х наш1ы.х [764].

Гравиметрические съемки

307

Иэвестн~18 кеnиброео~Н~18

.._....__,._ _., параметр~•. немэвестн~•е

nараметры каnиОровки м ~рейфа

Аliсоnютн~•е оnредеnени"

силь1 тяжести. значения и ковармацмоннаR матрица

Уравненн~•е 3Ha'I8HMR м

ковармацмоннаяматрица

Рис. 9.4. Блок-схема программы уравнивания гравиметрической сети [764].

 

PREPRO

 

ADJUST

 

 

 

 

 

Рис. 9.5. Блок-схема программы GRAVIТY для микрокомпьютера НР 9816 (Институт геоде11111,

Ганновер).

 

 

PREPRO выполняет предварительную обработку

(разд. 9.2.3), а программа ADJUST -

уравнивание результатов, при этом блок GEN содержит дополнительные данные (абсо­

лютные значении, веса). При стандартной оперативной памяти компьютера (512 Кбайт)

в обработке могут участвовать 1400 измерений (абсолютные величины g, приращения) и

может быть до 75 неизвестных. Для предварительной обработки одного измерения с вве­

дением поправки за земные приливы с ошибкой ± 1 (или ± 1О) нм· с- 2 необходимо соот­

ветственно 18 с (или 6 с), а дли уравнивания гравиметрической сети самых больших

размеров - 10 мин.

9.2.5.Оценка качества сети м оптимизация

Методики, разработанные для оценки качества и оптимизации геодезических се­ тей [511, 638], можно использовать и для гравиметрических сетей.

308

Глава 9

О качестве сети судят, с одной стороны, по ее точности и надежности, а с

другойпо денежным затратам. Для геодинамических сетей особое значение

приобретает их чувствительность при обнаружении временных изменений. Точность сети описывается кофакторной матрицей (матрицей обратных ве­

сов) Qxx (9.11), позволяющей получить дисперсии отдельных значений силы тяже­

сти g; по формуле

(9.20а)

С использованием (9 .15) и (9 .16) получим выражение для дисперсии приращений

силы тяжести gj - g;:

 

 

s5 =

ilJ(q;; + Qjj - 2qu).

(9.20б)

где q;;, Qjj , Qu- элементы

матрицы Qxx·

 

Надежность сети характеризует степень ее устойчивости по отношению к

ошибкам математической модели (разд. 9.2.4). При геодезических работах высо­ кую надежность обеспечивают, выполняя избыточные измерения при различных условиях наблюдений (разд. 9.2.2). «Локальную» надежность обеспечивают до­ статочным числом (более двух) контрольных измерений на каждом пункте.

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

мерений это обстоятельство иллюстрируется табл. 5.1 в разд. 5.1.6. Таблица 9.1 содержит

данные по региональной опорной сети, созданной с относительными гравиметрами (Ла Коста - Ромберга) при неблагоприятных условиях. Эти данные указывают на большие

расхождения между результатами, полученными одним и тем же прибором и разными

приборами. Они достигают 0,5 мкм·с- 2, что гораздо больше величин, получаюшихся при

нормальном распределении. В микрогравиметрической сети (благоприятные условия,

ошибка средних значений по внутренней сходиМости ± 0,01 мкм·с- 2 ) встречаются расхож­

дения, превышающие 0,1 мкм·с- 2 (табл. 9.2).

Таблица 9. 1. Изменение измеренного

приращения силы тяжести в гравиметрической

сети в Северной Исландии [598]

Приращения силы тяжести, мкм·с- 2

Ла Коста-Ромберr G-79

Ла Коста-Ромберr G-85

 

 

51,71

51,90

,72

67

,92

,91

,57

,72

,83

,63

,87

,65

,87

,67

,92

52,05

Таблица 9.2. Изменения средних прирашений

силы тяжести в микросети (МБМВ, Севр) [34]

Среднее приращение силы тяжести, мкс ·с- 1

Прибор

1 - 3

4 - 5

Ла Коста

- Ромберга

 

 

 

 

D

14

-0,80

+0,02

D

21

,82

,04

D

38

,71

,08

Gl31

,85

-0,02

G

253

,75

-0,03

G

258

,79

+0,02

Среднее 51,80

51,78

 

Среднее

-0,79

+0,02

 

 

 

 

 

 

Гравиметрические съемки

309

Точность и надежность сети можно оценить статистическими тестами,

основанными на вычислении точностных характеристик и определении поправок

из уравнивания. С помощью таких тестов оцениваются: -однородность сети (одинаковая точность по всей сети);

-значимость параметров моделей, которые невозможно определить физически

(члены калибровочной функции и функции дрейфа более высоких порядков);

-наличие грубых ошибок. Для их исключения эффективен метод контроля гру­

бых ошибок, называемый «data snooping» [19];

-параметры стохастической модели (9.5) по оценкам дисперсионно­

ковариационных компонент [339]. Эта модель особенно эффективна при оцен­ ке отношений дисперсий между достаточно большими и сильно связанными группами измерений (например, выполненными разными приборами); реже поддается оценке физическая корреляция.

Врезультате тестирования выявляют слабые участки сети. Их можно уси­

лить, выполнив дополнительные измерения и уточниi3 параметры моделей. Луч­

ше, если при создании сети выполняется постоянный качественный контроль с

дополнительными измерениями на слабых участках.

При создании высокоточных гравиметрических сетей (фундаментальные, гео­

динамические сети) рекомендуется использовать методы математической опти­ мизации сети; соответствующие методики для геодезических сетей уже

разработаны [238, 600]. Поскольку в общем случае исходный уровень силы тяже­ сти устанавливается из абсолютных измерений на пункте, условия размещения которого определены в разд. 9.2.1, оптимизации подлежит лишь проект распре­ деления измерений в сети (создаваемой или обновляемой) [236]. При этом целе­ вая функция определяет точностные требования для данной сети (кофакторная матрица Qxx), которые должны удовлетворять заданным условиям (максималь­

ная стоимость). После определения сетевой матрицы А находят весовую матрицу

Pu. Если точность отдельных измерений известна, задача сводится к определе­

нию числа повторных измерений для каждой связи.

В качестве целевых функций можно использовать обобщенную точностную

характеристикуслед матрицы Qxx(tr Qxx) или критериальную матрицу Qxx [20].

Последняя характеризует заданную точность сети; помимо среднеквадратических

ошибок значений g в эту матрицу могут быть также введены' например, точност­

ные характеристики дg, зависящие от расстояния между пунктами. Оптимальное

согласование между Qxx и Qxx можно получить по методу наименьших квадратов

или с помощью линейного программирования. Затем выполняется оптимизация

сети по имеющимся программам уравнивания, в которых точность, получаемая

в результате каждого дополнительного наблюдения, максимизируется с каждым шагом. В работе [590] исследована процедура формирования и аппроксимации

критериальных матриц для гравиметрических сетей.

Результаты оптимизации сильно зависят от вида целевой функции, припятой модели ошибок измерений (например, зависимости измеренных приращений си­

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

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

310

Глава 9

чески и при необходимости корректировать; желательно решать эту задачу по­ следовательными приближениями [457].

При создании фундаменrrшльной гравиметрической сети ФРГ 1976 г. оптимизация была выполнена с помощью линейного пршраммирования [61] последовательным мето­ дом [760]; ее результаты в значительной мере определили окончательную структуру сети (разд. 9.2.6). Планирование основывалось на использовании k = 21 пункта с приближен­

ными значениями g и на 4 пунктах абсолютных измерений (±0,15 мкм·с- 2). Пункты

предстояло соединить 40 связями с относительными гравиметрами (ошибка

среднего из

16 измерений ±0,05 мкм·с- 2, ошибка единичного измерения ±0,2 мкм·с- 2).

При после­

довательной оптимизации точность сети характеризовалась среднеквадратической ошиб­

кой силы тяжести:

 

 

-

~rQ,,

(9.21)

Sg = So

-- .

 

k

 

Целевая функция задавалась следующим образом:

 

tr Qotr Q;

 

 

=> max.

(9.22а)

С;

Здесь Qo и Q; - кофакторю>Iе матрицы величин сиЛы тяжести, полученных для начально­

го решения и после наблюдения i стоимостью с;. В качестве стоимостного параметра ис­ пользовалось время, необходимое для одной гравиметрической связи (расстояние s;):

С; = LJ.l;(h) = 1

s.,~..,

 

+ ~-

(9.22б)

Систематический анализ всех возможных k(k - 1) : 2

связей позволяет определить опти­

мальное место i-го измерения в сети, а последовательная обработка всех 40 запланирован­

ных измерений в конце концов определяет структуру сети. Пренебрежение с; приводит к

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

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

Ромберга позволяет определить общий хаvактер распределения точности; он ил­

люстрирован рис. 9.6. В будущем благодаря созданию пунктов абсолютных из­

мерений можно . будет существенно снизить рост ошибок при увеличении

величины дg (или соответственно при увеличении времени транспортировки).

(мкм·с-2)

Ошибка6g

1,00

Рис. 9.6.

Точность уравненных прнраwеннй си­

лы тяжести в высокоточной сети, по­

строенной с несколькими калиброван­ ными гравиметрами Ла Коста - Ромберrа.

Гравиметрические съемки

311

9.2.6.Национальные гравиметрические сети

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

можно условно разделить на три этапа, каждый из которых имеет свои

особенности.

До 1930-х rr. национальные сети строились в Потедамской системе относи­

тельными маятниковыми измерениями (разд. 9.1.1). В дальнейшем стали все

больше применять пружинные гравиметры. При создании фундаментальных се­

тей вначале использовали и маятниковые измерения, однако в конце концов их роль свелась к определению гравиметрических калибровочных базисов. Начиная с 1970-х гг. начали создавать новые национальные сети, привязывая относитель­ ные измерения к сети МГСС-71 либо используя совместно относительные и абсо­ лютные измерения; при этом число абсолютных измерений росло (разд. 9.1.2). По­ строение новых гравиметрических сетей стало особенно необходимым вс.;1едствие:

-появления системы МГСС-71 ;

-утери заметного количества пунктов имеющихся сетей;

-повышения требований к точности и надежности сетей;

-возрастающего интереса к мониторингу временных измерений силы тяжести;

-возможности использовать новые абсолютные и относительные приборы и

новые методики.

Национальные гравиметрические сети, как правило, делят на несколько клас­ сов; сети разных классов создаются последовательно, а при их обработке значе­

ния силы тяжести на пунктах высших классов считают твердыми для пунктов

более низких классов (разд. 9.2.4). Преимущества такой иерархической структуры

втом, что гравиметрические данные могут использоваться до завершения по­

следнего этапа сгущения сети. Среднее расстояние между пунктами опорных се­

тей в значительной степени определяется территорией страны (для

фундаментальной гравиметрической сети это 100300 км, в сети 1-го класса­

несколько десятков километров, в сети 2-го классаоколо 10 км, в сети 3-го

классанесколько километров). Приведем примеры.

В Федеративной Республике Германии гравиметрическая опорная сеть 1976 г.

(DSGN 76) была создана Немецким геодезическим исследовательским институ­

том под эгидой Немецкой геодезической комиссии [639]. Сеть состоит из 21 пунк­ та, каждый из них имеет три пункта-спутника; на 4 пунктах выполнены

абсолютные определения Институтом метрологии (Турин) (рис. 9.7). Согласно

оптимизации сети, были сделаны 44 гравиметравые связи. Все они измерены

дважды инезависимо в последовательности 1-2, 2-1 двумя бригадами, каждая из которых располагала двумя гравиметрами Ла Коста - Ромберга; всего полу­

чено 656 приращений силы тяжести. Пункты закреплены металлическими марка­

ми в пилонах или фундаментах зданий (рис. 9.8); как правило, пункты располагались в капитальных общественных зданиях; пример описания пункта приведем в разд. 9.4.1. Общее время измерений составило 200 дней. В результате

уравнивания по методу наименьших квадратов величины силы тяжести получены

с ошибкой ±0,06- 0,11 мкм·с- 2; максимальное приращение силы тяжести

(8000 мкм ·с- 2) определено с ошибкой ± 0,15 мкм ·с- 2Последующие измерения

предполагается проводить примерно каждые 10 лет.

312

 

 

Глава 9

 

 

 

 

 

 

 

Северное

 

 

 

 

 

море

Баnтийское море

 

 

54°N

 

 

 

 

 

 

 

100 200 км

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 9.7 (левый). Гравиметрическая опорная сеть 1976 (DSGN76) Федеративной Республики Германии

[692].

Рис. 9.8 (правый). Марка пункта гравиметрической опорной сети 1976 (ФРГ).

Геодезические ведомства выполиили сгущение DSGN 76 сетями трех классов.

В табл. 9.3 приведеиы информация о nлотности пунктов и другие данные для

земли Нижняя Саксония (47 000 км2) [285]. Гравиметрические пункты совмещены

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

и 2-го классов имеются пункты-спутники. На рис. 9.9 приведена схема гравимет­ рической сети 1-го класса в Нижней Саксонии. Сеть ФРГ 1-го класса DHSN 82 была полностью измерена и уравнена в системе DSGN 76 с применением соответ­

ствующей ковариационной матрицы; сеть содержит 289 пунктов, которые опре­

делены с ошибкой ±0,08- 0,14 мкм·с- 2 [132].

Сети DSGN 76 предшествовали маятниковая сеть, которая создавалась начиная с

1934 г. Потедамским геодезическим институтом для обеспечения гравиметрической съемки Германии (около 100 пунктов с шагом 50100 км) [596] (разд. 9.3.6), а также Немецкая гравиметрическая опорная сеть DSGN 62 [374]. Последняя содержала 43 пункта (9 маятни­ ковых пунктов, измерения с гравиметрами Аскания, Северная Америка и Уордена), она

была построена с 1950 по 1957 г. в виде петель большой протяженности и приведена в

Потсдамскую систему (± 1 мкм·с- 2).

В США гравиметрическая опорная сеть, состоящая из 59 пунктов (у каждого не менее трех пунктов-спутников), была создана Береговой и геодезической съем-

Гравиметрические съемки

313

Рис. 9.9. Сеть, 1-ro класса земли Нижняя Саксония, Федеративная Ресnублика Германии (385].

кой США в 1966-1967 rr. и обновлена в 1975-1976 rr. с 4 гравиметрами Ла

КостаРомберга; измерения выполняли по схеме петли. После уравнивания по­

лучена точность ±0,2 мкм·с- 2• Абсолютные измерения силы тяжести образова­

ли Среднеконтинентальный калибровочный базис от Форт-Дейвис, Техас, до

Грейт-Фолс, Монтана (диапазон силы тяжести 16 000 мкм·с- 2). Национальная

геодезическая съемка Национального управления США по изучению океана и ат­ мосферы, а также ГИдрографический и топографический центр Боенно­ картографического агентства МО США планируют создание фундаментальной абсолютной гравиметрической сети; она будет служить отсчетной для всех от­ носительных измерений [531] (рис. 9.10). Сеть состоит из 11 основных и 40 допо­ лнительных пунктов, среднее расстояние между пунктами 300 км. Измерения выnолняются абсолютным гравиметром JILA (разд. 5.2.3). В Канаде националь­ ная гравиметрическая сеть также укрепляется пунктами абсолютных определений

[403].

Китайская гравиметрическая опорная сеть 1985 г., содержащая 57 пунктов,

была создана в 19821985 rr. Государственным геодезическим и картографиче-

Таблица 9.3. Гравиметрические сети сгущения земли Нижняя Саксония [385]

Класс

Плотность

Число

Число

Метод

Продолжительность

 

пунктов

пунктов

приборов

измерений

наблюдений

1-й

111000 км2

56

~]

Разности 1-2, 2-1

-1 год

2-й

1/100 км2

425

В треугольниках

-5лет

3-й

115 км2

-9000

1 - 2

Измерения no схеме

-13 лет

 

 

 

 

nетли no nрофилям

 

в nрямом и обратном

наnравлениях

314

Глава 9

Рис. 9.10. Проект rравиметрйческой оnорной сети абсолютных оnределений США:

t. - nервокласс­

ные,

о -

nервоклассные доnолнительные, о- второклассные nункты

[531J.

ским бюро

[650]. Она включает 6 пунктов абсолютных

определений

(±О,1 мкм ·с- 2 ).

Относительные измерения сделаны с 9 гравиметрами Ла

Коста - Ромберга по схеме петли ( ± 0,2 мкм ·с- 2); выполнены связи с внешними

пунктами сети МГСС-71 , которая не распространялась на территорию Китая. Ре­

зультаты уравнивания имеют ошибку ± 0,08 мкм ·с- 2 (для максимального прира­ шения силы тяжести ошибка равна ±О,13 мкм ·с- 2 ).

9.3.Региональные и локальные гравиметрические с"Ъемки

9.3.1.Цели работ, распределение пунктов

Региональные и локальные гравиметрические съемки выполняют для решения за­

дач геодезии, геофизики и геодинамики (разд. 4.2- 4.4). Для трансформации и интерпретации данных о гравитационном поле обычно по· измеренным величи­

нам силы тяжести находят гравитационные аномалии, на основе которых с по­

мошью интерполяции получают единообразное представление поля. Сами значения силы тяжести на земной поверхности необходимы для изучения времен­ ных вариаций. Такие исследования рассмотрены в разд. 10.1.1.

Плотность пунктов и точность гравиметрической съемки зависят от разме­

ров территории, вида и величины исследуемых особенностей поля силы тяжести или распределения аномальных масс, порождающих это поле. При региональных

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

те (например, в горной местности пункты не выбирают только в долинах).

Ошибка точечных аномалий определяется ошибками измерений силы тяжести и

гравиметрической редукции (разд. 9.3.3, 9.3.4).

Вгеодезии для вычисления высот гравиметрического геоида или квазигеоида,

атакже уклонений отвеса необходимо знать аномалии в свободном воздухе для

Гравиметрические съемки

315

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

кальных задач проводится лишь для исследуемой территории и ее окрестностей

(разд. 4.2.2, 4.2.3). На линиях нивелирования значения силы тяжести требуются для определения ортаметрических или нормальных высот (разд. 4.2.4). Для гео­

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

±10 мкм·с- 2 при среднем расстоянии между пунктами в несколько километров.

В таком случае можно определять превышения регионального геоида или квазигеои­

па, а также высоты в гравитационном поле на сантиметровом уровне точности, а уклоне­

ния отвеса -

с ошибкой не более ± l н ; остаточные ошибки из-за влияния

высокочастотных составляющих гравитационного поля не превысят соответственно

0,001 м и 0,4" (разд. 3.3.2). В области со сложным рельефом для точного определения

уклонения отвеса на пункте необходима дополнительная съемка повышенной плотности

вокруг этого пункта в ближней зоне радиусом в несколько километров [518]. Если притя­ жение топографических масс можно учесть по моделям рельефа с высоким разрешением,

то расстояние между пунктами съемки может быть до 5 -10 км.

При национальных гравиметрических съемках в целях геофизики желательны расстояния между пунктами в 1 - 5 км; точность аномалий Буге должна быть

± 5 - 1О мкм ·с- 2 На больших территориях в начальной стадии съемки плот­

ность пунктов меньше (шаг 1020 км) и располагаются они в основном вдоль дорог; так же поступают при рекогносцировочных съемках в прикладной геофи­

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

(определение глубины кристаллического фундамента, поиск локальных осадочных

бассейнов, соляных структур, тектонических нарушений) и при разведке полезных ископаемых (разд. 4.3.5). В зависимости от конкретного объекта могут предъяв­

ляться специфические требования к шагу съемки, который меняется от 50 до

500 м, и точности измерений от ± 1 до ± 5 мкм ·с- 2 В угольных районах и при

рудной разведке, как правило, расстояние между пунктами составляет несколько

десятков метров, а ошибка измерений не превышает ± 1 мкм ·с- 2При микрогра­

виметрических исследованиях, создании инженерных сооружений и поиске пустот

шаг съемки равен нескольким метрам, а ошибка составляет

±О,1 мкм ·с- 2

Если национальная гравиметрическая сеть имеет

высокую

плотность пунктов

(разд. 9.2.6), она может полностью или частично служить

региональной гравиметриче­

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

Поэтому чаше всего задачу

решают поэтапно путем сбора, предварительной обработки

и комбинирования данных,

предоставляемых различными ведомствами (разд. 9.4.2).

9.3.2.Методика с"Ьемки и обработка результатов

Региональные и локальные гравиметрические съемки выполняют исключительно с пружинными гравиметрами. Удобно, когда к гравиметрическому пункту мож­ но подъехать на автомобиле (рассчитанном на переезды по .шоссе, в некоторых случаяхна автомобиле, пригодном и для бездорожья). В труднодоступных

районах возможна транспортировка на вертолете. При локальных и специальных

съемках (например, подземных съемках в шахтах) с гравиметром выполняют и

316

Глава 9

Рис. 9.11. Местная гравиметрическая опорная сеть для региональной съемки; невя-зки в мк\l·с ' (535].

пешие переходы. Для снижения затрат на определение плановых координат и вы­

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

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

Если национальная гравиметрическая сеть не очень плотная или не обеспечи­

вает точность, нужную для специальных задач, создают опорные сети для де­

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

кратные измерения), преимущественно по схеме петли; пример приведен на

рис. 9.11 [535]. Значения силы тяжести находят по методу наименьших квадратов

(разд. 9.2.4).

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

каждые два часа, с тем чтобы контролировать смещение нуль-пункта и обнару­

живать его скачки. Невязки распределяют пропорционально времени, а прилив­

ными поправками часто пренебрегают. Такая методика обеспечивает точность по

внутренней сходимости ±0,5 мкм ·с- 2 , при повышении точностных требований

соответственно возрастает объем дополнительных работ.

Вре.мя измерений на пункте составляет всего несколько минут (отыскание пункта, установка прибора, нивелирование, дезарретирование, приведение индекса

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

мер для повышения точности (разд. 6.6. - 6.6.3) время работы на пункте не бу­ дет более 10 мин. В зависимости от условий транспортировкИ за день можно

определить (при однократном наблюдении) 50 (или 30) пунктов при расстояниях между ними в несколько сотен метров (или несколько км). Для быстрого выпол­

нения детальной съемки можно использовать неисчерпанные возможности верто­

летной гравиметрии (разд. 7.4.2), а также наземной и аэроградиентометрии

(разд. 8.3.2).

3.2.1, 4.3.3).
в масштабе
При локальных съемках может потребоватьси точность

Гравиметрические съемки

317

9.3.3.Определение планового положения и высот пунктов

Плановые и высотные координаты гравиметрических пунктов устанавливают их

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

тограмметрии. Выбранный метод должен соответствовать требуемой точности

с учетом экономических аспектов и рассчитан на имеющиеся материалы (катало­ ги пунктов геодезических сетей, топографические карты, аэрофотоснимки).

В отличие от пунктов геодезических сетей положение гравиметрических пунктов, как правило, не требуется устанавливать с максимальной точностью; исключение составляют

сложные опорные сети, предназначенные для мониторинга временных измерений силы тя­

жести (разд. 10.'1.1). Лучше выбирать такие методы, при которых временные затраты со­

измеримы с затратами на измерения силы тяжести и обеспечивается требуемая точность.

Опыт работ показывает, что геодезическое обеспечение съемки часто недооценивают, осо­

бенно в неизученных районах.

Для определения положения гравиметрических пунктов нужно определить их

плановые координаты, с тем чтобы эти пуикты можно было опознавать на кар­

тах, а последующие оценки выполнять в численном виде. При региональных съемках, как правило, Достаточна точность планового положения ± 1 О - 100 м

(что соответствует ошибкам от ± 0,2 мм в масштабе карты 1 : 50 000 до ± 0,4 мм

1 : 250 000).

координат ±0,1 - 1 м.

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

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

разд. 9.3.4.

Дифференцируя формулу (3.7б) для редукции в свободном воздухе оgсв.в' по­ лучим зависимость ошибки этой редукции от ошибки определения высоты

(9.23)

Если d(оgсв.в.> = О,1 мкм ·с- 2 (точность измерения силы тяжести), то высоту

нужно находить с ошибкой ±0,03 м. Для вычисления топографической редукции

(4.36) требуется знать отметку пункта и высоты рельефа. Разделив на части

(4.38), объединим поправки за притяжение пластины ognл (4.37) и редукцию в

свободном воздухе в общее выражение, зависящее от высоты. Дифференцируя

это выражение, найдем

dН<м> = (3 - 0,0004е)- 1d(оgсв.в - оgпл><мкм·с-'>·

(9.24)

Чтобы получить редукцию с ошибкой О,1 мкм ·с- 2 = 2700 кгм- 3 ),

высоту не­

обходимо знать с ошибкой ±0,05 м. Для вычисления поправки за рельеф с ошиб­

кой 1 мкм·с- 2 (равнинный район) или 10 мкм·с- 2 (горы) разрешающая

способность и точность топографических моделей (цифровых моделей рельефа)

должны быть различными в зависимости от сложности рельефа и расстояния

до притягиваемой точки.

318 Глава 9

В ближней зоне (0- 100 м) либо не должно быть сильных нарушений рельефа, либо

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

с ошибкой

±5-

10 м и разрешением до 0,1 км х 0,1 км; это возможно,

если имеются

топографические карты масштабов 1 : 5000 -

1 : 25 000. Влияние зоны с радиусами 2 -

50 км можно учесть по модели местности

с разрешением 1 км х 1 км

(ошибки отме­

ток ± 10100 м), такая модель может

основываться на картах масштаба 1 : 25 000-

1 : 250 000.

Для

области за пределами

50

км необходимо разрешение

10 км х 10 км

( ± 100 м). Предпочтительно для разных зон использовать разные оценочные алгоритмы

[55].

При вычислении аномалий в свободном воздухе и аномалий Буге необходимо

находить нормальную силу тяжести -у0, которая является функцией геодезиче­

ской широты ~. На основании (3.1в) имеем

 

 

 

 

1

 

19

d-yo

_2

(9.25)

 

 

 

d<P(pa ) == --~~~~ d-yo

== - . --

 

 

 

д

)'еfЗ Sin 2'{)

 

2'{)

(мс

)

 

Для

d-yo ==

0,1

мкм·с~ 2

наибольшая

требуемая

точность (<Р == 45°)

составит

d<P==

±0,4"

==

12м.

 

 

 

 

 

 

Геодезические методы определения координат и высот можно разделить сле­ дующим образом: методь1 определений по имеющимся картам и аэрофотосним­

кам, методы космической геодезии и традиционные наземные геодезические

методы. Для обеспечения региональных и локальных гравиметрических съемок,

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

[64, 243, 258, 377, 626, 691].

Координаты и высоты точек можно взять с топографической карты, если положение

гравиметрических пунктов опознается на ней однозначно. По картам масштабов 1 : 25 ()(Х)

(1 : 50 000) плановые координаты можно получать с ошибкой примерно ± 10 м ( ± 20 - 50 м), а высоты- с ошибкой ± 1 - 5 м (для точек с подписанными отметками ±0,3- 1 м). Для определения координат можно также использовать трансформированные аэро­ фотоснимки масштаба 1 : 10 000 - 1 : 50 000, по стереопарам можно получить высоты пунктов [325]. Фотографии из космоса (станция Space1ab 1) уже сейчас обеспечивают пла­

новые и высотные координаты с ошибкой ± lV м [376]. Пункты обоснования, необходи­

мые для пользования аэрофотоснимками, определяют методами, которые описаны ниже. Среди методов, основанных на наблюдениях космических объектов, для определения

координат гравиметрических пунктов методы геодезической астрономии, требующие

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

широты и долготы этими методами ± 1" соответствуют в линейной мере ± 30 м; в дей­

ствительности эта ошибка на порядок больше из-за неизвестных уклонений отвеса, кото­

рые определены с точностью ± 1" - 2" лишь для некоторых территорий. Большое значение приобрели методы космической геодезии [192, 623]. Можно приобрести приемо­

индикаторы, позволяющие оперативно определять пространствеиные декартовы коорди­

наты пунктов в общеземной геоцентрической системе (разд. 2.1.1 ). Эти координаты

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

высоту) (разд. 2.4.2). Навигационная спутниковая система ВМС США (система Транзит)

начала работать в штатном режиме с 1970 г. и будет функционировать примерно до

1995 г.; при измерении доплеровской частоты с использованием бортовых эфемерид она

обеспечивает определение координат с ошибкой ± 1О - 30 м по одному прохождению

Гравиметрические съемки

319

ИСЗ; накапливая измерения в течение нескольких суток,

можно уменьшить ошибку до

± 3 - 5 м. Из одновременных наблюдений на двух или нескольких станциях, расположен­

ных не далее 500 км друг от друга, разности их координат получают с ошибками ± 0,2 - 0,5 м [378]. Вводится в строй (завершение около 1990 г.) Глобальная система определения местоположения Новстар (GPS); в любой точке на Земле можно будет одн~временно на­

блюдать не менее 4 спутников этой системы.

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

масштабе времени [365, 622]. При использовании кода С/А (coarse acquisition- код моду­

ляции несушей частоты L1), который предназначен для рядовых гражданских пользовате­

лей, местоположение в реальном масштабе времени определяется с ошибкой ± 30 м. При

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

Наземными .методами измерений в~шолняют по существу интерполяцию между

опорными пунктами, координаты которых определены национальными геодезическими

службами или получены спутниковыми методами. Для определения пространствеиных ко­ ординат на больших расстояниях с ошибками ±0,1 - 0,5 м на 1050 км можно исполь­ зовать инерциальные геодезические системы [613, 621]. Рекомендуется такие измерения сочетать с GРS-измерениями. При локальных съемках можно использовать метод тахео­ метрии (дециметровая точность). Для определения плановых координат можно приме­ нять классические методы, основанные на измерении расстояний и углов (электронный дальномер и теодолит); измерения можно выполнять по различным схемам (например, ходами), ошибки определения координат исчисляются сантиметрами или дециметрами.

Особо важно определение высот. Высокую точность (±0,01 м/1 км- ±0,1 м/100 км и точнее) обеспечивает геометрическое нивелирование, однако оно требует больших затрат времени. Более экономично, особенно в горных областях, тригонометрическое нивелиро­

вание (или тахеометрия); оно почти на порядок грубее геометрического нивелирования. Барометрическое нивелирование позволяет определять превышения с ошибкой примерно

±дН(м)/100 при условии, что надежно учитываются изменения атмосферного давления и дрейф прибора (применение барографа, привязка к реперам, использование нескольких ба­

ронивелиров). Характерно, что этот метод не зависит от наличия прямой видимости меж­

ду пунктами; в этом отношении он аналогичен гравиметрическому методу и потому

весьма экономичен. Даже при неблагоприятных условиях (транспортировка на вертолете

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

±1- 10 м [201].

При обработке геодезических измерений необходимо помнить, что наряду с общезем­

ной системой координат существует множество национальных координатных систем. Они,

как правило, не являются геоцентрическими и в качестве поверхностей относимости испо­

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

да сочетают методы космической геодезии с наземными методами и когда пользуются

имеющимися опорными пунктами и картами [64]. Параметры перехода (от локальной гео­ дезической системы к геоцентрической системе WGS 84) даны в работе [147].

320

Глава 9

При определении плановых координат и высот для национальных гравимет­

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

9.3.4.Определение плотности

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

чины плотности необходимо и для вычисления аномалий в свободном воздухе

с поправкой за рельеф (разд. 4.2.3, 4.3.3).

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

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

Для вычисления поправок нужна трехмерная модель реального распределения

плотности горных пород. Такая обьемная плотность определяется на основе

плотности минерального скелета ем горных пород с учетом содержания в их

порах воды р (ОТо) по формуле

(9.26)

где е.- плотность воды; предполагается, что в естественных условиях почти

всегда наступает насыщение. Объемные плотности основных горных пород и ми­

нералов приведены в табл. 4.1 (разд. 4.3.5). Реально удается моделировать лишь

наиболее крупные комплексы горных пород (порядка км3). В этом случае необхо­

димо знать обьемную плотность формации, определяемую объемной плотное-

тью горных пород с учетом трещин и степени выветривания.

Рассмотрим точность определения плотности, необходимую для вычисления

аномалий Буге. Топографическую редукцию можно представить состоящей из по­

правки за притяжение плоского промежуточного слоя ognл и поправки за рельеф

ogP (разд. 4.3.3):

 

 

ДgБ = g - -дg

Н -

{jg + {jg

-

"'О ,

 

 

 

дН

 

р

1

 

или

дgБ = g-

3,086Н- 0,000419еН + ogp-

'УО

МКМ·С- 2 .

(9.27)

 

Из этого соотношения следует связь ошибки определения плотности и ошибки

поправки за промежуточный слой:

(9.28)

 

 

Гравиметрические съемки

 

321

Если поправка

необходима с

точностью ± 1 мкм ·с - 2, то

при Н = 100 м

(или

1000 м) допустима

ошибка

определения

плотности

± 24 кг · м- 3

(или

± 2,4 кг· м- 3),

или в

относительной мере 1OJo

(или 0,1% ).

 

 

Поправка за рельеф линейно зависит от плотности. Если нужно определить

эту поправку с точностью 1 мкм ·с - 2, то допустимая относительная ошибка в

определении плотности составит 1OOJo для равнинных территорий, а в горных

районах уменьшится до 1%.

Для вычисленШI топографической редуiЩИИ используются различные подходы. Как пра­ вило, для больших геологических формаций или даже для всего района считают, что распре­ деление плотности с глубиной однородно; это приводит к понятию однородной плотности

Буге. Часто используют значение плотности 2670 кг· м- 3 (кристаллические горные поро­

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

иногда также выполняют специальные определения плотности (см. далее). Так, плотность

консолидированных осадочных пород примерно равна 2300 kг ·м- 3, а плотность магмати­ ческих и метаморфических пород - 2700 кг· м- 3 Из-за смешанного состава горных пород

и сильных локальных изменений плотности при региональной съемке не имеет смысла

подразделение по плотности < 50 кг· м - 3 [ 656]. ДетализацШI плотностной модели по

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

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

тальные и наклонные плоскости). Если вычислять топографическую редукцию, отсчиты­

вая высоту от выбранной поверхности (например, проходяшей через самую низкую точку),

то можно уменьшить ошибки аномалий силы тяжести из-за погрешности в принятой

плотности глубоких масс. ·

Для вычисления редукции в свободном воздухе необходимо знать вертикаль­ ный градиент силы тяжести дg/дН. Из (9.27) имеем

 

1

(9.29)

 

d(дg/дН)с.z = - -- d(оgсв.в>(мс·'>·

 

н(м)

 

При d(оgсв.в.> =

1 мкм·с- 2 и Н= 100 м (или 1000 м) допустимая ошибка гради­

ента ±10 (или

±1) нс- 2, или в относительной мере 0,3% (или 0,03%). Если по­

льзоваться стандартным нормальным значением градиента, ошибка может быть

существенно больше (разд. 8.1.2) 0 .

В региональных и локальных исследованиях особую роль играют разности аномалий

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

гипотезы о постоянных величинах плотности и вертикального градиента сильно коррели­

рованы. При расстоянии между пунктами в несколько километров коэффициент корреля­

ции г вейнчины плотности может составлять 0,5 - 1 [152]. Из выражения (9.20)

следует,

что в таком случае ошибка разности аномалий существенно уменьшается в (1 -

r) 112 раз.

Методы определения плотности подробно описаны в геофизической литера­ туре [462, 505]. Дадим лишь их обзор.

11 ) Аномалию в свободном воздухе находят на физической поверхности Земли, пользуясь нор­

мальным значением д-у!дН вертикального градиента силы тяжести, поэтому указанная ошибка не

возникает. - Прим. ред.

322

Глава 9

Объемную плотность можно определить непосредственно, измерив массу и

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

пружинных весах дважды: сначала в воздухе, а затем погрузив его в воду. Проис­

ходит потеря веса, равная весу воды, вытесненной образцом. Уменьшение веса определяется формулой

(9.30а)

где т- масса образца, v - его объем, твмасса вытесненной воды. По­

скольку

(9.306)

то плотность образца будет равна

(9.31)

где евплотность воды.

Взвешивание nозволяет оnределить nлотность с высокой точностью. Однако из-за слабой nредставительности образца небольших размеров (вес его менее нескольких кило­ граммов) результаты сильно меняются как на ограниченном участке, так и по региону. К систематическим ошибкам nриводят изменения nри извлечении образца из горной фор­ мации (уменьшение давления, nотеря каnиллярной воды). Для восстановления насышения

образца водой его nриходится смачивать. Если число взятых образцов достаточно пред­

ставительно для территории, то средняя nлотность горной формации может быть оnреде­

лена с ошибкой ±50- 100 кг·м- 3 (табл. 9.4).

Таблица 9.4. Плотность горных nород, nолученная взвешиванием их образцов, Северная Исландия [598]

Горная порода

Количество

Плотность, кг· м- 3

СредНЯЯ ПЛОТНОСТЬ,

пунктов

образuов

 

кг·м -J

Третичные базальты

5

27

2

500

-

3 ()()()

2

830

±

70

Третичные туфы

4

19

2

100

 

2 600

2 330

±

110

Кварцевые базальты

14

89

2

100

 

3 ()()()

2

650

±

60

ГИалокластические

6

40

1 700

 

2 500

2

100

±

80

лороды (туфы

и брекчия)

Объемную плотность формации для структур ограниченной протяженности

можно определять по измерениям силы тяжести с использованием соотноше­

ния, связывающего зависимость силы тяжести от высоты с плотностью пород

(разд. 3.2.2). Этот метод при~еним при измерениях как на земной поверхности,

так и на вертикальных профилях.

Как показано в работе Неттлтона [503], определение плотности топографи­

ческих структур гравиметрическим методом основывается на предположении,

Гравиметрические съемки

323

280

 

2400~

 

 

 

240~~=

 

200

 

--=-~--2700~

 

 

 

 

- 2800

- •

1~

~

 

Н(мl

~~~о 2 4 6 -8 '10 12 км

Рис. 9.12. Построенный по методу Неттлтона профиль Вадлахейди, Северная Исландия: слои .1авы 11 J

третичных базальтов с туфовыми включениями, nлотность 2700 кг-м-·'; см. табл. 9.3, ..:о..:тав­

ленную по данным [599].

что аномалии Буге, вычисленные с истинным значением плотности пород, ведут

себя как случайные величины, не коррелированные с высотами пунктов. Возника­

ющую же корреляцию можно устранить, изменяя припятое значение плотности.

Этот метод применим для гравиметрических профилей (рис. 9.12) и для площад­ ных съемок; оценки можно получить графически или из вычислений [342, 418].

Приняв nриближенное значение плотности Qo по (9.27), можно вычислить nриближенные аномалии Буге дg~. Сравнение их с «истинными» аномалиями,

полученными с плотностью Буге

 

 

(9.32)

приводит

к выражению

 

 

дg6 = дg~- 0,000 419·Н·дg- og: де МКМ·С-2,

(9.33)

 

е

 

где og~ -

поправка за рельеф, вычисленная с плотностью ео. После исключения

региональной составляющей (например, среднего значения аномалии) остаточные

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

равен

(9.34)

-

где дg6 и Н- соответственно средние значения аномалии и высоты. Подстав-

ляя выражение (9.33) и считая, что для поправки за рельеф достаточна плотность

ео, получим

l:((дg~- дti)(H- Н))

 

 

 

.

(9.35)

 

де = ----- ·-·-----~

 

0,000419·l:(H- Н)2

'

 

могут потребоваться несколько приближений.

324

Глава 9

Одновременно с определением плотности из уравнивания по методу наимень­ ших квадратов можно определить и региональную составляющую поля. Если представить функцию LlgБ поверхностным полиномом, уравнение поправок в со­

ответствии с (9.27) имеет вид

n

т

V = ~

~ aij ~yi- (3,086- 0,000419е) Н- g- ogp +')'О МКМ·С- 2(9.36)

i=O

}=О

В результате уравнивания находят коэффициенты полинома aij и плотность е.

Гравиметрический метод Неттлтона - очень эффективный метод определения объем­ ной плотности топографических масс в области с гравиметрической съемкой; ошибка

определения плотности ± 2050 кг·м- 3 (т.е. относительная ошибка 1 - 20Jo). При ис­

пользовании метода предполагается следуюшее:

-распределение плотности не зависит от рельефа. Если же рельеф определяется строени­

ем нижележащих геологических структур, то метод не работает;

-однородное распределение плотности в исследуемой (ограниченной) области;

-достаточно большие (более 100 м) перепады высот н представнтельное распределение

гравиметрических пунктов как в плане, так н по высоте;

--'- удовлетворительная модель регионального поля, влияние которого необходимо исклю­

чать нз используемых данных.

Вертикальные гравиметрические профили можно построить в шахтах и сква­

жинах (прнменение скважинных гравиметров (разд. 7.1.2)) [398]. Из измерений силы тяжести в двух точках Р1 и Pz, лежащих на одной вертикали, можно опре­ делить объемную плотность для слоя (толщина Llz = zz - ZI) между этими точка­ ми (рис. 9.13). Вычисления выполняются в три этапа. Сначала введеннем

редукции og~~n исключают влияние топографических масс (притягиваемая точка

Р1), затем для разности высот между Р1 н Р2 вводят редукцию в свободном воз­

духе, и после этого редукцией og~~n восстанавливается влияние масс (притягивае­

мая точка Р2). На основании (9.27) получим

Рис. 9.13.

Вертикальный гравиметрический nрофиль (шахтный nрофиль) llЛЯ

оnределения nлотности nород.

Гравиметрические съемки

325

или

(9.37)

Переходя к дифференциалам, выражение (9.37) можно преобразовать в дифферен­

циальное уравнение Пуассона (2.22).

Обьемная плотность пород вычисляется по формуле:

(дglдz)дz -

дg + og< 0

+ og<2>

(9.38)

е =

4-л·Gдz Р

Р '

причем вместо дglдz обычно используется нормальное значение 3086 не- 2

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

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

включить также параметры изменения плотности по высоте и аномалию вертикального

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

Относительная ошибка определения плотности пород не превышает 1OJo.

В скважинах можно также определять плотность близлежащих пород (в пределах от нескольких сантиметров до нескольких дециметров), используя метод гамма-гамма каро­

тажа, основанного на эффекте Комптона. Гамма-излучение радиоактивного источника те­

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

пород.

Результаты разных методов определения плотности могут различаться до

50Jo. Расхождения вызваны локальными изменениями плотности и ее отличием

от плотности формации, изменением плотности извлеченного образца от nлотности породы in situ и модельными погрешностями гравиметрического метода (допушение об однородности слоев, использование нормального градиен­

та силы тяжести). Рекомендуется применять разные методы и сочетать их ре­

зультаты.

9.3.5.Морская и аэрогравиметрическая с-ьемки

В морской и воздушной гравиметрии выполняют квазинепрерывные измерения на подвижном основании (разд. 7.3). На съемочных профилях накапливается

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

на отрезке маршрута около l км. Расстояние между профилями выбирается для

каждого проекта и лежит в пределах от 10 до 50 км при работах на море и от 1 до 5 км - на вертолете. Для определения положения профиля, вычисления по­

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

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

шинстве случаев достаточен лишь линейный масштабный коэффициент

(разд. 6.4.5). Абсолютный уровень силь1 тяжести определяют привязкой к опор­ ным пунктам. Для морской гравиметрии эти пункты создают в портах, в аэро-

326

Глава 9

гравиметрии исходные значения получают путем аналитического продолжения

силы тяжести на высоту полета над данным пунктом. Повторные привязки в

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

съемке в виде сетки пересекающихся профилей (рис. 9.14). В большинстве случаев

профили проектируются в направлениях с севера на юг и с запада на восток;

иногда их ориентируют в зависимости от исследуемой структуры (перпендику­

лярно ее простиранию). Невязки в точках пересечения профилей зависят от точ­

ности измерений силы тяжести, поправки Этвеша и определения координат. Эти

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

янную поправку. При измерениях на море расхождения можно уменьшить до

±5- 20 МКМ·С- 2 (разд. 7.4.1, 7.4.2).

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

Для определения положения пунктов, как правило, достаточно знать коорди­

наты с ошибкой

±50100 м. Для вычисления поправки Этвеша в средних ши­

ротах с

ошибкой

± 5 мкм ·с- 2

необходимо соблюдать

следующие требования:

скорость судна должна определяться с ошибкой ±0,05

м/с, а курс- с ошибкой

±0,5°,

скорость

самолета необходима с ошибкой ±0,04 м/с, а его курс -

с

ошибкой ±0,03° (разд. 7.2.4). Глубины должны измеряться с ошибкой ± 1 -

5

м

(при этом ошибка поправки за промежуточный слой не будет превышать

1 -

4 мкм ·с- 2 ), а высота полета -

с ошибкой не более

± 3 м (что соответствует

ошибке редукции за высоту ± 10 мкм ·с- 2 ). Более жесткие требования предъявля­

ются к точности определения изменений высоты полета (до ±О, 1 м), если их ис­

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

N

~ог----------------------,--~

Испыта­

тельный

540

9 Морские 1?0

мили

Рис. 9.14. Морские гравиметрические профили в голландской части Северного моря (Институт геодезии Делфтского технологического университета [660).)

Гравиметрические съемки

327

Измерение

sz

силы тАжес:ти

/J,

.:;;;:г·~-

~р>О

Рис. 9.15.

Вычисление аномалии Буге на море: е н- nлотность воды,

и - nлотность земной коры, bg" - nоnравка за рельеф:

Аномалии в свободном воздухе и аномалии Буге на акваториях используют для гео­

физической интерпретации. При вычислении аномалий Буге плотность масс воды допол­

н·яют до стандартной плотности земной коры или плотности, характерной для

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

димо помнить, что последняя может быть как положительной, так и отрицательной

(рис. 9.15). Считая,

что плотность

земной

коры

е= 2670 кг·м- 3,

а плотность воды

е.=

1030 кг·м- 3, из

соотношений

(4.36) и

(4.37)

получим

 

 

 

 

АgБ = g + 0,687Т + ogp- -уо мкм·с-

2,

(9.39)

где

Т- глубина, ogP- поправка

за рельеф.

 

 

 

Известно много методов навигации морских и воздушных судов [168, 321, 624]. Их можно разделить на определения положения (непрерывного или дискретного)

и определения скорости и курса. Все они могут базироваться или на внешних

навигационных системах (глобальные, региональные, временные локальные), или

на автономных. Особую роль играют спутниковые, радионавигационные и инер­

циальные методы, которые часто объединяют в интегральные навигационные

системы.

Для спутниковой навигации по всей поверхности земного шара создана Спутниковая навигационная система ВМФ США NNSS (система ТРАНЗИТ). Наблюдения спутников, которые можно выполнять через интервалы в несколько часов (в зависимости от широты места) дают координаты с ошибкой ±50100 м. Следуюшая новая системаГлобаль­ ная система определения местоположения (GPS). Система GPS обеспечивает определение

местоположения в реальном масштабе времени (±50 м), скорости

(±О, 1 м/с) и привязку

ко времени (±50 не) [162]. Зная высоты геоида, полученные по

глобальным моделям

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

ров), установи.в на суше (на исходном пункте) опорный приемоиндикатор, а на судне -

бортовой приемоиндикатор, методом дифференциальной навигации можно определять их

взаимное положение с ошибкой ± 5 м и менее. После развертывания системы GPS (при­

мерно к 1990 г.) она обеспечит непрерывную навигацию на всем земном шаре. Если будет открыт доступ к точному коду системы для любого пользователя, то появится возмож-

328 Глава 9

ность обойти многие трудности при гравиметрических съемках•>. Астрономические опре­

деления координат из наблюдений звезд и Солнца ( ± 2 - 4 км) в значительной мере

утратили свою роль (разд. 9.3.3).

В радионавигационных методах применяют наземные базисные станции, излучающие

электромагнитные волны; сигналы примимаются приемоиндикатором подвижной станции

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

они пересекаются) и от состояния атмосферы. Существуют различные системы: система Омега (10 кГц, дальность 8000 км, ошибка ± 1 - 3 км) является глобальной, большая часть поверхности Мирового океана входит в зону действия системы Лоран-С (100 кГц,

дальность до 2000 км, ошибка ± 100 м- 1 км), базисные станции цепочек системы Декка (100 кГц, дальность до 450 км, ±20100 м) установлены по основным судовым фарвате­ рам. Рабочие зоны высокоточных систем ( ± несколько метров), например Силедне (450 МГц) или Хай-Фикс (2 МГц), охватывают прибрежные акватории (до 150 -200 км). Для локальных съемок можно также применить подвижные системы морской и воздуш­ ной навигации. К их числу, например, относятся система Мини-Рэйнджер (из измерений

времени прохождения микроволновых радиоимпульсов определяют расстояния от подвиж­

ной станции до базисных, дальность 180 км, ошибка ±3 м), а также дальномерная систе­

ма ANA (фирма PraklaSeismos, Ганновер) с односторонним распространением радиоволн. В этой системе, разработанной для морских и аэрогеофизических исследова­

ний, на передающей и приемной станциях установлены атомные стандарты частоты;

дальность действия системы 900 км, ошибка ± (10 м + 8 м/час). С использованием судо­

вого радара и антенн, установленных на суше или буях, местоположение судна при удале­

нии на несколько десятков километров можно определять засечкой с ошибкой ±50 м.

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

применить доплеровский гидролокатор, использующий отражение звуковых волн от мор­

ского дна (при глубинах до 200 м скорость судна определяется с ошибкой ±0,05 м/с). При аэронавигации можно использовать метод доплеравекого радара. Он обеспечивает опреде­ ления курса и скорости относительно земной поверхности на высотах О, 1 - 1 км, а в соче­

тании с гиросистемой и при интегрировании по времени позволяет определять положение

относительно исходных пунктов (ошибка не более 1'Jfo от длины маршрута) [63]. Инерци­

альные навигационные системы представляют собой универсальное средство навигации

(ошибка ±0,05 м/с), однако при увеличении времени наблюдений их ошибки быстро воз­ растают (разд. 4.5.2, 7.5).

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

ния курса и скорости, для которых характерна кратковременная стабильность. Объединив

их с автопилотом, можно исключить быстрые изменения траектории на съемочном про-

l) Завершен этаn развертывания космического сетмента GPS. По состоянию на февраль 1994 г. он содержит 24 рабочих и 3 резервных ИСЗ. В отечественной практике широко используется система

ГЛОНАСС, во многом схожая с системой GPS; точность мгновенного определения плановых коорди·

нат -100 м; в феврале 1994 г. система содержала 13 ИСЗ. После завершения развертывания (1995 r.) она будет состоять из 24 ИСЗ (включая 3 резервных), размещенных в трех плоскостях на высоте 19100 км с nериодом обращения 11,25 ч и наклонением 64,8°.- Прим. ред.

 

 

 

 

 

 

 

Гравиметрические съемки

 

 

 

 

 

329

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Морской гравиметр

 

 

 

 

 

Источник 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

KssЭO

 

 

 

 

1 ПМТIНМА

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Несrnаженн~оtе

г---....1...'---,2 МIМ8МАПА

 

 

 

 

 

 

 

~Г;::оо=о=мме=том=-,см=~'::-rм

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

==~=~·г-Г.-f»фо----,

 

 

 

 

 

"ММСММЙ AIT"'MI(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rмен~о1й процессор

 

 

 

nocrpourwл•

 

 

 

 

 

Gss ~~бnо ·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~•!!!nе~к~тро~нм!!!к~м_(;кi<iон;;;;;троЮ"n;;;..;;;0-;:.i

 

 

 

с0"tJ:':.:..~~"~мм

В~о•бор

 

 

 

~=':,

 

 

 

 

 

 

Гмроnмтформзмермтеn~t,Н~А

 

 

 

 

 

мнтерфейсt.•

оnерациИ

 

 

 

 

 

 

 

 

 

М1 с бnоком

 

 

информация

 

 

 

 

 

 

 

 

Д.ннwем

 

 

'-.._д.:;;,;,;]).

 

 

 

 

 

 

:мектронмкм

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

КОМ8НА1о:~

 

 

 

 

 

 

 

 

 

 

 

 

СГМ*8ННМ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Гр88М...Т~КМ8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

М АР'/ГМ8 До8ННW8

 

 

 

 

 

 

 

 

 

 

 

 

====...:..~=...,.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

на...гечмоннан ммформацмfl

 

 

г-...

 

AII.,"коррекциигироскоnа

 

 

 

L..--.1--

 

 

 

 

iНавм..==~rацмоннlоiЙ информация

 

 

 

 

 

Интерфейс

 

в

 

 

 

 

 

 

 

 

 

 

 

Скорост~о

 

 

.----------

~~~-~

 

 

 

 

 

 

 

 

 

 

 

 

 

CVAМI

 

 

 

 

 

 

 

 

СVАНОМ М ГIОА8314'18СКОГО

 

 

1

 

 

IM.-niГ

 

 

 

 

 

 

 

 

 

 

oбecnr.ettмA

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Гnавнwй

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

комnwотер An"

 

 

 

 

 

 

 

 

 

 

 

 

Эxol18'f

 

 

Гnу6мн1

 

 

 

 

 

8W....СМНМА

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

KOOPAtfН8T,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

текушей

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~бoTKtfM

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ч1а..1

 

 

Д.та/ерем11

 

 

 

 

 

tмеТ8МIТм.JIЦМИ

 

_/rМа-----rнмтнаlll''i

 

 

 

1

 

 

 

 

 

 

 

 

 

 

A8HHioiX

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

НР &845

 

 

 

 

 

 

 

 

 

 

 

 

 

Курс CVAНI

 

 

 

 

 

Ре~мкм

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n.тот..

 

 

 

 

 

 

 

 

 

 

1 Гирокомп«

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~~

 

--------~

 

 

 

 

мта

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

СмС:ТеМI НIВИГIЦММ

 

 

 

 

 

 

 

 

 

 

 

 

т

--

 

 

 

 

 

 

 

М обреботкм AIHH ..IX

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~~====~---------

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ПocneAYIOW8"

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

обработка

Рис. 9.16. Схема nолуче1шя ннформашш в морско~• rраuнмстрс 1...:"30

(uключая <.:IKTC\1/,1 11<11!111 ;1111111 11

обработк11 данных) [540]

 

 

(материал ф11рмы Prakla

- Sci,mo,. Гашюuср).

филе. Особенно удачным оказалось сочетание спутниковых навигационных систем (ТРАН­ ЗИТ, GPS) и/или радионавигационных систем с компасными блоками и инерциальными навигационными системами [775]. Интегральные системы обеспечивают в любом месте

земного шара точность ± 100200 м, в прибрежных областях достижима точность не

хуже ± 10 м.

Измерения глубины выполняют акустическими эхолотами с точностью примерно ±O,IOJo (на мелководье около ±0,1 м). Высоту полета измеряют, главным образом ис­ пользуя одновременно баравысотомер и радиовысотомер или лазерный высотомер (при

полете над водной поверхностью); на сравнительнонебольших высотах ошибка составляет

± 1 м. В качестве вспомогательных средств могут использоваться система GPS и инерци­ альные геодезические системы; в будущем эти системы позволят регистрировать и учиты­

вать короткопериодические вертикальные возмущающие ускорения по высокоточным·

измерениям вариаций высоты полета.

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

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

8.2.1).

330

Глава 9

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

гравиметрических редукций. Невязки в точках пересечения профилей и аналоговая

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

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

вузлы координатной сетки и построением карты изоаномал (разд. 9.4.3).

9.3.6. Примеры региональных и локальных

гравиметрических СDемок

Рассмотрим некоторые примеры региональных и локальных съемок.

В Германии гравиметрическая съемка была осуществлена в 1935-1945 гг. в рамках национальной геофизической съемки под руководством Геологической службы [125]. Связями с пунктом абсолютных определений в Потедаме была со­

здана маятниковая сеть (±5- 10 мкм·с- 2). После включения в нее опорных се­

тей, созданных с гравиметрами Тиссена-Шлезенера и Аскания, плотность этой

сети была повышена и расстояния между пунктами составили от 2 до 5 км

(±1- 3 мкм·с- 2) (разд. 9.2.6, 9.3.2). В обширных равнинных районах Германии

была также проведена съемка с вариометрами (съемка по профилям, расстояние между профилями несколько километров, шаг по профилю несколько сотен мет­ ров) (разд. По результатам 45 тыс. измерений силы тяжести были созда­ ны карты масштаба 1:200 000 (аномалии Буге). На территории ФРГ съемка была продолжена в 1950-х гг. Государственной геологической службой земли Нижняя

Саксония (Ганновер). Вслед за построением опорных сетей была выполнена съем­

ка территорий, ранее недостаточно обеспеченных гравиметрическими пунктами

либо совсем не изученных. Применялись приборы фирмы «Аскания» и кварцевые

пружинные гравиметры ( ± 0,2 - 0,3 мкм ·с- 2 ). По этим данным, включающим

также съемки с вариометрами и протяженные локальные съемки, выполненные

частными организациями (поиски нефти), создается карта аномалий Буге масшта­

ба 1:500 000 (3 листа) с сечением изоаномал 10 мкм ·с- 2 [538] (рис. 9.17). Значения

силы тяжести отнесены к системе МГСС-71, а нормальная сила тяжести опреде­

лена в Геодезической референц-системе 1967 г. (разд. 3.1.2). Притяжение проме­

жуточного слоя вычисляли для сферического сегмента (радиус 166,7 км) с

плотностью 2670 кг· м- 3 , влияние топографических масс при вычислении поправ­

ки за рельеф учитывалось на расстояниях до 23 км. В ГДР создание гравиметри­

ческой сети (с гравиметрами Аскания) было завершено в 1961 г. [592]. Она

разделена на 3 класса (расстояния между пунктами примерно 50, 15 и 3 - 4 км),

в последующей региональной съемке среднее расстояние между пунктами соста­

вило 1,5 км.

Во Франции национальная гравиметрическая сеть была создана Управлением

научно-исследовательских работ по геологии и шахтам (Орлеан). Сеть содержит

Гравиметрические съемки

331

Рис. 9.17.Фраrмент карты аномалий Буrе ФРГ, М. 1:500000, сечение изоаномал 10 мкм·с- 2 [538].

6 пунктов абсолютных определений; относительные измерения выполнены че­

тырьмя гравиметрами Ла Коста - Ромберга на 51 пункте первого класса и 324 пунктах второго. После сгущения сетями 3-го и 4-го классов съемка будет содер­

жать примерно 370 тыс. пунктов, на ее основе будут созданы гравиметрические

карты Франции масштабов 1:320 000

и 1: 1 000 000 (плотность осадочных пород

2300 кг·м- 3 , плотность фундамента

2700 кг· м - 3 , сечение изоаномал 1О -

100мкм·с- 2).

ВСША в 1975 г. были начаты совместные работы Общества геофизиков­

разведчиков (SEG), Национальной геологической съемки, Картографического

управления МО и Национального управления по исследованиям океана и атмос­ феры. Собраны данные, имевшиеся в различных организациях; после обра­

ботки составлены карты аномалий масштаба 1:2 500 000 на территорию США

[301]. Карта сечением 50 мкм·с- 2 основана примерно на 1 млн. гравиметриче­

ских пунктов на суше (аномалии Буге, плотность пород 2670 кг·м- 3 ) и примерно

на 0,8 млн. прибрежных пунктов (аномалии в свободном воздухе) [515]; в число

этих пунктов вошли измерения силы тяжести по линиям нивелирования, выпо­

лненные Национальной геодезической съемкой. Аномалии определены в системе МГСС-71 и Геодезической референц-системе 1967 г. В горных районах поправки

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

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

картографировании положения и строения фундамента, кристаллических пород,

332

Глава 9

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

пределения экструзивных и интрузивных магматических пород.

Примерам поэтапного создания региональной гравиметрической карты для

большой территории может служить карта Индии [3]. Национальная программа

гравиметрических работ предусматривает площадную съемку со средним шагом 15 км, опирающуюся на фундаментальную сеть (с шагом -100 км). К 1986 г. ре­

гиональная съемка покрывала примерно 55% территории страны. Она дополнена

специальными геофизическими съемками и измерениями по нивелирным линиям, в результате к 1986 г. имелось 26 тыс. пунктов; распределение пунктов иллюстри­ рует рис. 9.18. Гравиметрические данные представлены аномалиями в свободном воздухе, Буге и изостатическими аномалиями.

Во многих областях Земли были выполнены .региональные гравиметрические съемки с целью изучения процессов тектоники литосферных плит, эти работы, как правило, сочетались с другими геофизическими методами (разд. 4.4.1). Отме­ тим региональные исследования в области Центральных Анд (Кордильер), пред­ ставляющих собой классический пример горной формации, образовавшейся в

результате субдукции океанической коры под континентальную [225] (рис. 9.19).

С 1982 по 1986 г. гравиметрами Ла Каста - Ромберга было определено около

2600 пунктов на площади 850 х 500 км. Обычно измерения выполнялись вблизи

транспортных магистралей (с шагом около 5 км), высоты пунктов определяли,

как правило, барометрическим нивелированием с привязкой к нивелирной сети.

Плановое положение находили по топографическим картам .масштаба 1:50 000-

1:500 000. Вслед за созданием базы гравиметрических данных, содержащей все

имеющиеся значения силы тяжести и цифровую модель местности (для трапеций

3' х 3'), была создана карта аномалий Буге с сечением 100 мкм·с- 2 (МГСС-71,

~10-20км

0>2Окм

90' Е

Рис. 9.18 (левый). Схема региональных гравиметрических съемок Индии по состо11нию на 1986 г. [3].

Рис. 9.19 (правый).Гравиметрические профили н аномалии Буrе региональной съемки в Центральных

Андах (система МГСС-71, плотность 2670 кr·м- 3), сечение сглаженных изоаномал

500 мкм·с- 2 [228].