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

4

АНАЛИЗ И СИНТЕЗ СЛОЖНЫХ СИСТЕМ

УДК 531.5

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

О.М. Яшникова

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

Ключевые слова: аномальное гравитационное поле, точечные массы, карты изолиний.

Введение

При решении ряда прикладных задач требуются карты различных характеристик гравитационного поля Земли (ГПЗ), их трехмерные изображения и соответствующие им изолинии, а также профили вдоль различных направлений. В частности, такая потребность возникает при исследовании плотностных не-однородностей земной коры, разведке полезных ископаемых, оценке сейсмологической обстановки и т.п. Измерение величин, характеризующих ГПЗ, может проводиться гравиметрическими методами с использованием данных гравиметров, гравитационных вариометров и градиентометров. По данным съемки строятся гравиметрические карты, характеризующие аномалии ГПЗ [1, 2]. Поскольку гравиметрические съемки проведены далеко не на всей поверхности Земли, особенно в труднодоступных районах с аномальной изменчивостью поля, возникают задачи описания параметров гравитационного поля в условиях отсутствия данных гравиметрической съемки. Для этого создаются различного рода математические модели поля исследуемого района.

Один из вариантов построения модели аномального ГПЗ может быть основан на так называемом методе точечных масс. Теоретические основы этого метода известны и довольно подробно освещены в работах [3-5] применительно к описанию возмущающего потенциала [6].

В настоящей работе ставится задача построения карт изолиний аномального гравитационного поля на основе метода точечных масс и определение параметров ГПЗ, характеризующих изменчивость аномальной составляющей поля как в пространственной, так и во временной областях. Для этих целей в прикладном пакете МаИаЪ разработана компьютерная программа, реализующая описываемую ниже математическую модель гравитационного поля и позволяющая на основе заданных характеристик точечных масс строить карты изолиний интересующих параметров. Приводятся результаты апробации методики построения карт для района озера Байкал.

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

Согласно закону тяготения Ньютона проекции вектора силы притяжения Ж, направленного от притягиваемой точки Р с единичной массой и координатами (хР, уР, хР) к притягивающей точке М с массой т и координатами (хм, ум, хм), могут быть представлены [7] в векторной форме в виде

т

= О —СО8(г,] ), ] = х, у, х , г

или в скалярной форме в виде

Рш = Ст( УР  ]М ) , у = х, у, х ,

] г3 ^ ,

где О - гравитационная постоянная; г = ^(хр - хм)2 + (ур - ум)2 + (хр - хм)2 .

Помимо силы притяжения, на каждую точку на Земле действует центробежная сила, однако она не зависит от наличия притягивающих масс, поэтому можно считать, что сила тяжести для единичной массы определяется только силой притяжения, т.е. gx = Ех, gy = ¥у, gх = Ех. Если на точку Р действуют не одна, а п притягивающих точечных масс, то выражения для компонент вектора действительной силы тяжести g = gy, gх )Т (здесь индекс «Т» означает транспонирование) могут быть представлены в виде

<ти

?] = О 7

= О X" " , ] = х, у, х .

*=1 гк

Проекции вектора силы тяжести на плоскости географического сопровождающего трехгранника (ГСТ) позволяют определить направление вектора силы тяжести. Примем точку Р за начало прямоуголь-

ной системы координат, ось г совместим с направлением вектора нормальной силы тяжести у = [0,0,у] ось х направим по касательной к меридиану на север, а ось у - на восток (рис. 1).

Р

8у V / -

! 1 1 /*-. Л / £ Я; Г \ \

I / 1 / :/ У 8

у г

Рис. 1. Проекции вектора силы тяжести

Из (1) можем получить выражения для интересующих нас составляющих аномального ГПЗ:

£ = -?х, п = -?у,

8 г

а? = & - у,

где £ и л - так называемые продольная и поперечная составляющие ГПЗ, а А? - аномалия силы тяжести. В выражении (2) ввиду малости углов тангенсы заменены значениями самих углов, выраженными в радианах.

Уравнения (1)-(3) и определяют используемую далее математическую модель аномального ГПЗ. Для ее конкретизации необходимо выбирать точечные массы и соответствующие им координаты. Поскольку процедура подбора этих параметров носит итерационный характер и требует многократных вычислений, была разработана программа, позволяющая автоматизировать процесс построения таких моделей и соответствующих им карт изолиний в зависимости от числа задаваемых точечных масс и их характеристик.

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

Процедура построения карт изолиний аномального ГПЗ с использованием модели (1)-(3) на основе точечных масс заключается в следующем. Географическая карта интересующего района разбивается на равные участки, в центре каждого из них размещается точечная масса. С использованием географических карт определяются координаты хм., ум. основных притягивающих масс Mi (рис. 2, а). Для выбранного количества точечных масс строятся изолинии аномалий. Необходимое количество точечных масс определяется из условия непересечения изолиний и их сглаживания с критерием волнистости изолинии не более 5%. Данный критерий характеризует степень приближенности расчетной изолинии к реальной и представляет собой соотношение длин этих изолиний [8, 9]. При добавлении точечных масс разработанная программа автоматически рассчитывает величину достигнутого критерия и отображает новые изолинии.

Описанная процедура была апробирована при построении карт изолиний сильно аномального гравитационного поля в районе озера Байкал (рис. 2). Байкал лежит в глубоком тектоническом понижении, и местность вокруг озера имеет сложный, глубоко расчлененный рельеф. Горные хребты тянутся здесь параллельно один другому в направлении с северо-запада на юго-восток и разделены котловинообраз-ными понижениями. Высота большинства хребтов Забайкалья составляет около 1300-1800 м, а наиболее высокие хребты достигают значений около 3000 м.

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

В этом случае величина массы определялась выражением р-£ • к

™ы, =

3

где р - плотность вмещающей среды; к - площадь основания и высота пирамиды соответственно. Согласно работам [10, 11], за среднюю плотность вмещающей среды, составляющей Байкальскую рифто-вую зону, принимается р = 3-103 кг/м3. В качестве высоты залегания точечной массы выбиралась высота расположения центра масс пирамиды.

X

<

о •

/о

, • «Притягивающая» масса

о «Отталкивающая» масса

^г^даялгавмш. в?^--¿г /у//. а б

Рис. 2. Схема расположения точечных масс

Аналогичная аппроксимация проводилась и для глубин озера с заменой в формуле (4) средней плотности среды на плотность воды р = 1-103 кг/м3. Масса в этом случае бралась с обратным знаком, так как является «отталкивающей», т.е. уменьшающей значение полной силы тяжести.

Далее с использованием значений масс т, и соответствующих им координат хм., ум,, хм, в

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

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

В ряде случаев требуется знание аномалий ГПЗ вдоль заданной траектории, например, когда точка Р прямолинейно перемещается в исследуемом районе. Удобно, чтобы зависимость параметров ГПЗ представляла собой функцию от времени. Для этого необходимо задаться траекторией движения внутри поля. Предположим, что объект перемещается прямолинейно из точки в точку Е (рис. 3).

Рис. 3. Траектория движения

В этом случае координаты вычисляются как Хрг = X0 + V • (/ - /о) • яп(К), Ур, = У + V • (/ - /о) • СО8(К), где Хо,Уо - координаты точки старта (точка 5 на рис. 3); V - скорость движения объекта; / - текущее время; /о - начальное время; К - курс объекта.

В программе предусмотрена процедура реализации траектории в виде координат Хр,,Ур, для каждого момента времени и вычисления параметров ГПЗ л, Дg для каждой точки траектории.

Полученные параметры ГПЗ в проекциях на плоскости ГСТ приведены на рис. 4 и представляют собой карты изолиний, описывающие характер аномального гравитационного поля в данном районе.

х105

6 5 4

^ 3

2 1

х105

х105

м

х105

б

Рис. 4. Параметры ГПЗ в проекциях на меридиональную

Аномалия силы тяжести, а также продольная и поперечная составляющие аномального ГПЗ относительно касательной к траектории движения были получены для курса К =135° (рис. 3) и скорости V =8 узлов и приведены на рис. 5.

20 10

о 0

н

и

-10

-20 -30

1

1

2

12 10 8

1 6 4 2

0 4000 8000 12000 16000

Время, с

4000 8000 12000 16000 Время, с

б

Рис. 5. Продольная £

Из рис. 5, б, видно, что аномалия силы тяжести для приведенной траектории (рис. 3) имеет вид параболы, значения которой возрастают до 12 мГал по мере приближения к горному массиву и опускаются фактически до нуля в середине траектории. Подобный характер кривой, описывающей аномалии ГПЗ в данном районе, полностью согласуется с данными, полученными на основе гравиметрической съемки [11, 12].

Заключение

Описаны процедура и реализующая ее в пакете МаИаЪ программа построения карт аномального гравитационного поля Земли на основе метода точечных масс.

Практическая апробация процедуры выполнена для акватории озера Байкал, характеризуемой сложным, глубоко расчлененным рельефом. Характер изолиний параметров ГПЗ, полученных с помощью используемой модели, с достаточной степенью точности соответствует аналогичным картам, полученным гравиметрическим методом [11, 12].

Литература

1. Непоклонов В. Б. О создании цифровых карт уклонений отвесной линии // Геодезия и картография.

- 1996. - № 9. - С. 2-6.

м

а

0

а

Ю.М. Гордеев, Т.В. Падерина, Д.А. Соколов

2. Хмелевской В.К. Геофизические методы исследований: Учебное пособие для геофизических специальностей вузов. - Петропавловск-Камчатский: изд-во КГПУ, 2оо4. - 232 с.

3. Фомин В.Н. О представлении гравитационного поля простейших тел притяжением точечных масс // Астрономия и геодезия. - Томск: Изд-во ТГУ, 198о. - Вып. 8. - С. Ш2-1Ю.

4. Антонов В. А. Представление гравитационного поля планеты потенциалом системы точечных масс // Труды астрономической обсерватории ЛГУ. - 1978. - Т. 34. - С. 145-155.

5. Полещиков С.М. О построении вещественной системы точечных масс, представляющих гравитационное поле планеты // Вестник ЛГУ. - 1984. - Т. 14. - С. 95-99.

6. Гравиметрия и геодезия / Под ред. Б.В. Бровара. - М.: Научный мир, 2о1о. - 57о с.

7. Огородова Л.В., Шимбирев Б.П., Юзефович А.П. Гравиметрия. - М: Недра, 1978. - 325 с.

8. Бурик О.В. Исследование методов автоматического построения изолиний // Сборник трудов магистрантов Донецкого национального технического университета. Вып. 2. - Донецк: ДонНТУ. Министерство образования и науки Украины, 2оо3. - С. 331-342.

9. Чернова Л.И. Обработка геопространственной информации при цифровом моделировании топографических задач // Автореферат дис. к.т.н. - Иркутск, 2ооб. - 2о с.

10. Савиных В.П., Ященко В.Р., Ямбаев Х.К. Геодезические исследования геодинамики рифтовой зоны озера Байкал // Изв. вузов. Геодезия и аэрофотосъемка. - 2ооб. - № 5. - С. 3-1о.

11. Кабан М.К. Гравитационная модель коры и верхней мантии Северной Евразии // Российский журнал наук о Земле. - 2оо1. - Т. 3. - № 2. - С. 143-163.

12. Милановский Е.Е. Геология России и Ближнего Зарубежья (Северной Евразии): Учебник. - М.: Изд-во МГУ. - 1996. - 448 с.

Яшникова Ольга Михайловна - ОАО «Концерн «ЦНИИ «Электроприбор», аспирант, olga_evstifeeva@mail.ru УДК 531.746

ПРЕЦИЗИОННАЯ ГИРОСКОПИЧЕСКАЯ СИСТЕМА ПОДЗЕМНОЙ НАВИГАЦИИ ДЛЯ РАБОТЫ В ВЫСОКИХ ШИРОТАХ. ОСНОВНЫЕ СХЕМНО-КОНСТРУКТИВНЫЕ РЕШЕНИЯ Ю.М. Гордеев, Т.В. Падерина, Д.А. Соколов

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

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

Введение

Исследования российских геологов доказывают, что акватория Северного Ледовитого океана может стать крупнейшей нефтегазоносной провинцией на Земле. Однако насколько богат этот «арктический Клондайк», настолько и труден в освоении. Уникальность этих месторождений обусловлена, прежде всего, необходимостью строительства скважин в шельфовой зоне морей, имеющей сотни метров глубины, находящейся под воздействием чрезвычайно суровых природно-климатических условий. На сегодняшний день в мировой практике не существует аналогов работ, связанных с проектированием и разработкой подобных месторождений.

Строительство скважин арктического шельфа, как правило, с большими отклонениями ствола от вертикали, практически сложно осуществить без применения гироскопической системы подземной навигации, работающей в процессе бурения (так называемой gMWD-системы). Специфика режима точечного компасирования в такой системе обусловливает снижение ее точности в высоких широтах. Это означает, что на широтах 74-78° с.ш. (местоположение наиболее крупных месторождений углеводородов) погрешность вычисления азимута скважины gMWD-системой при условии применения гироскопа с одним и тем же уровнем ухода вырастет в 1,8-2,4 раза по сравнению с аналогичной погрешностью на широте 6о° с.ш. [1].

Для обеспечения подземной навигации в упомянутых широтах требуется разработка прецизионной гироскопической системы нового типа. И решение для создания такой системы есть. Это применение так называемой диаметральной схемы, предложенной некоторое время назад учеными ЦНИИ «Электроприбор» в процессе разработки серии универсальных гироскопических инклинометров (ГИ) для работы в скважинах с малым диаметром и произвольной ориентацией. Диаметральная схема представляется наиболее совершенной в ряду технических решений для бесплатформенных ГИ. Она предусматривает разворот корпуса гироскопа вокруг поперечной оси скважинного прибора (СП), при котором вектор кинетического момента остается перпендикулярным оси разворота и находится в диаметральной плоскости