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

Дементьев, Юрий Викторович

Бесплатный автореферат и диссертация по наукам о земле на тему
Развитие теории и разработка технологии определения аномалий силы тяжести в полной топографической редукции
ВАК РФ 25.00.32, Геодезия

Автореферат диссертации по теме "Развитие теории и разработка технологии определения аномалий силы тяжести в полной топографической редукции"

На правах рукописи

Дементьев Юрий Викторович

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

25.00.32 - «Геодезия»

Автореферат диссертации на соискание ученой степени доктора технических наук

1 9 ДПР 2012

Новосибирск - 2012

005018262

Работа выполнена в ФГБОУ ВПО «Сибирская государственная геодезическая академия».

Научный консультант -

Официальные оппоненты:

доктор технических наук, старший научный сотрудник Капеницкий Анатолий Иванович.

Мазурова Елена Михайловна, доктор технических наук, профессор, Московский государственный университет геодезии и картографии, заведующая кафедрой геодезии;

Пимшин Юрий Иванович, доктор технических наук, профессор, Волгодон-ский инженерно-технический институт Национального исследовательского ядерного университета, заместитель руководителя по научной работе и инновации;

Тимофеев Владимир Юрьевич, доктор физико-математических наук, старший научный сотрудник, Институт нефтегазовой геологии и геофизики им. A.A. Трофимука СО РАН, заведующий лабораторией физических проблем геофизики.

ФГУП «Центральный ордена "Знак Почета" научно-исследовательский институт геодезии, аэросъемки и картографии им. Ф.Н. Красовского» (г. Москва).

Защита состоится 29 мая 2012 г. в 13-00 часов на заседании диссертационного совета Д 212.251.02 при при ФГБОУ ВПО «Сибирская государственная геодезическая академия» (СГТА) по адресу: 630108, Новосибирск, ул. Плахот-ного, д. 10, ауд. 403.

С диссертацией можно ознакомиться в библиотеке ФГБОУ ВПО «CITA». Автореферат разослан 10 апреля 2012 г.

Ведущая организация -

Ученый секретарь диссертационного совета

Середович В. А.

Изд. лиц. ЛР № 020461 от 04.03.1997. Подписано в печать 28.02.2012. Формат 60 х 84 1/16. Пет. л. 1,75. Тираж 100 экз. Заказ 18 Редакционно-издательский отдел СГГА 630108, Новосибирск, Плахотного, 10. Отпечатано в картопечатной лаборатории СГГА 630108, Новосибирск, Плахотного, 8.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы исследования. Гравиметрия является наукой, изучающей поле силы тяжести Земли и методы её измерения. Измеренные значения силы тяжести в точках земной поверхности используются для определения фигуры и внутреннего строения Земли, выявления причин и выполнения оценок последствий геодинамических процессов природного и техногенного характера, выполнения геологического картирования и тектонического районирования территорий, поиска и разведки месторождений полезных ископаемых. Раздел гравиметрии, посвященный приложению гравиметрических данных к решению геодезических задач по определению фигуры и гравитационного поля Земли, приведению геодезических измерений на поверхность эллипсоида, называется геодезической гравиметрией. Высока роль в развитии геодезической гравиметрии Центрального ордена «Знак Почета» научно-исследовательского института геодезии, аэросъемки и картографии им. Ф.Н. Красовского (ЦНИИГАиК). Весомый вклад в решение задач геодезической гравиметрии внесли советские и российские ученые: Бровар В. В., Бровар Б. В., Бузук В. В., Вовк И. Г., Грушин-ский Н. П., Демьянов Г. В., Еремеев В. Ф., Каленицкий А. И., Кузьмин Ю.О., Магницкий В. А., Мазурова Е. М., Машимов М. М., Михайлов А. А., Молоден-ский М. С., Огородова Л. В., Остач О. М., Пеллинен Л. П., Шимбирев Б. П., Юркина М. И. и др., а также зарубежные ученые Гофман-Велленгоф Б., Мориц Г.

В настоящее время гравиметрическая аппаратура позволяет измерять силу тяжести с погрешностью (0,01 - 0,03) мГал, или (1 - 3)х10"8 относительно абсолютного значения. Методы пространственного координирования позволяют определить положение гравиметрического пункта и значение нормального гравитационного поля с такой же точностью. Однако, аномалии силы тяжести, особенно в районах со значительным перепадом высот, определяются в 50 - 100 раз грубее. Поэтому, проблема повышения точности определения аномального гравитационного поля Земли, как на её поверхности, так и вне её, в том числе . во внешнем пространстве является актуальной, особенно в последнее время \

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

знании.

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

иной области знаний.

При вычислении аномалий высот и уклонений отвесной линии требуется

знание величины аномалии силы тяжести в редукции за свободный воздух в точке, где измеряются геодезическая высота и сила тяжести. В условиях дискретной гравиметрической съемки для определения аномалий в промежуточных пунктах, где наблюдения силы тяжести не производились, появляется необходимость в интерполировании результатов измерений. Однако аномалии в свободном воздухе существенно зависят от изменения высоты рельефа местности. Как показано в книге Шимбирева Б. П. «Теория фигуры Земли», при расстояниях между пунктами от 20 до 30 км линейная интерполяция аномалий в свободном воздухе в горных районах приводит к погрешностям порядка (20-25) мГал.

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

ям в свободном воздухе. При этом погрешность интерполяции аномалий Буге в горных районах уже на порядок меньше и составляет от 2 до 3 мГал.

При использовании в косвенной интерполяции аномалий силы тяжести в неполной топографической редукции погрешности интерполяции ещё уменьшаются. Как показали исследования, при выполнении неполной топографической редукции с учетом влияния топографических масс постоянной плотности, расположенных в промежуточном слое, ограниченном снизу поверхностью отно-симости нормального поля, сверху - физической поверхностью Земли, для горной области в радиусе 200 км от точки наблюдений погрешность интерполяции составляет от 0,7 до 0,8 мГал. Вместе с тем, установлено, что значение топографической поправки постоянно возрастает с увеличением размеров учитываемой области, и нет предела, при котором можно «отбросить» влияние масс дальних зон.

В работе Пеллинена Л.П. и Остача О.М. (ЦНИИГАиК) показано, что определение высот квазигеоида с погрешностью ±1 м приводит к погрешности вычисления аномалий силы тяжести порядка ±0,25 мГал. Очевидно, что для расчета высот квазигеоида с субметровой точностью требуется знание аномалий силы тяжести с погрешностью не ниже 0,02 мГал. Такую точность может обеспечить использование аномалий силы тяжести в полной топографической редукции.

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

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

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

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

При выполнении полной топографической редукции силы тяжести предполагается, что весь учитываемый промежуточный слой разделяется на внутреннюю, внешнюю и дальнюю области. Внутренней областью считается участок местности, расположенный в непосредственной близости от точки наблюдений. Размер её колеблется от нескольких метров до нескольких километров в зависимости от изрезанное™ рельефа. Граница внешней области рекомендована Инструкцией по гравиразведке (1980 г.). Её радиус в среднем составляет 200 км. Дальняя область представляет всю поверхность Земли за исключением внутренней и внешней областей.

В настоящее время накоплен значительный опыт учета влияния рельефа местности во внутренней и внешней областях на значение силы тяжести. Систематизация этого опыта и новая автоматизированная методика выбора оптимального шага цифровой модели рельефа (ЦМР) местности во внешней учитываемой области в зависимости от масштаба и точности гравиметрических работ изложена в Методических рекомендациях по учету влияния рельефа местности, утвержденных МинГео СССР как практическое руководство в 1981 г. Эта методика позволяет определять аномалии силы тяжести как в редукции Буге, так и в неполной топографической редукции. Однако влияние топографических

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

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

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

Для достижения поставленной цели решались следующие основные задачи.

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

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

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

- разработка технологии реализации полной топографической редукции.

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

масс Земли.

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

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

Основные положения диссертации, выносимые на защиту:

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

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

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

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

Научная новизна исследований:

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

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

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

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

Научная значимость работы. Расширена тематика и область применения аномалий силы тяжести: в геодезии - при косвенной интерполяции и осредне-

НИИ аномалий силы тяжести в редукции за свободный воздух; в разведочной геофизике - при геологическом картировании и тектоническом районировании территорий, поисках и разведке полезных ископаемых. Получены необходимые формулы для вычисления топографической редукции в дальней области и разработана технология её реализации.

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

Соответствие диссертации паспорту научной специальности. Тематика диссертации соответствует пунктам: 1 - «Определение параметров земного эллипсоида и квазигеоида» и 7 - «Математическая обработка результатов измерений и информационное обеспечение топографо-геодезических работ» паспорта научной специальности 25.00.32 - «Геодезия».

Апробация и реализация результатов исследований. Основные положения диссертационной работы докладывались, обсуждались и получили одобрения на научно-практическом семинаре «Геометрия, физика, динамика и гравитационное поле Земли» (г. Новосибирск, 2009), на Международных научных конгрессах «ГЕО-Сибирь» (г. Новосибирск, 2008, 2009, 2010 и 2011 гг.).

Результаты выполненных исследований внедрены в учебный процесс СГГА (акт о внедрении от 15 декабря 2011 г.).

Публикации по теме диссертации. По теме диссертации опубликовано 18 статей (в том числе 10 статей в реферируемых изданиях, утвержденных ВАК для защиты докторских диссертаций).

Структура и объем работы. Диссертационная работа состоит из введения, пяти разделов, заключения, списка литературы из 57 наименований. Основной текст диссертации изложен на 169 страницах машинописного текста, содержит 30 рисунков и 29 таблиц.

Диссертация и автореферат диссертации оформлены в соответствии с СТО СГГА 012-2011.

ОСНОВНОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИИ

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

Первый раздел «Аномалии силы тяжести и их применение в геодезии и разведочной геофизике (аналитический обзор)» посвящен аналитическому обзору методов учета поправок за рельеф при вычислении аномалий силы тяжести. Приведена зависимость поправок за плоский и сферический слои от их толщины и радиуса учитываемой зоны. Описано применение аномалий силы тяжести в геодезии и разведочной геофизике. Дополнительно рассмотрены некоторые физические параметры, прямо или косвенно осложняющие определение и интерпретацию аномалий силы тяжести. К ним отнесены приливные деформации поверхности Земли и атмосферные нагрузки (снежный покров, уровень грунтовых вод и атмосферное давление). Приведены алгоритмы учета влияния этих факторов.

Во втором разделе «Методика учета влияния промежуточного слоя при выполнении полной топографической редукции» на теоретической основе реализована новая методика учета влияния всех топографических масс, заключенных в промежуточном слое, на вертикальную составляющую силы тяжести. При этом под промежуточным слоем подразумевается слой постоянной плотности, ограниченный снизу поверхностью принятого земного эллипсоида, сверху - реальной поверхностью Земли. Показано, что для реализации полной топографической редукции требуется построение общеземной цифровой модели рельефа по сфероидическим трапециям определенного размера. В узлах этих трапеций должны быть известны эллипсоидальные (геодезические) высоты. Разделение поверхности эллипсоида на сфероидические блоки позволяет один раз определить высотные характеристики местности и затем использовать их

многократно при выполнении редукции в любых точках наблюдений. 10

В третьем разделе «Построение ПЦМР и определение интегральных параметров топографической поправки дальней области» изложена технология реализации разработанной методики для учета влияния масс промежуточного слоя дальней области на вертикальную составляющую ускорения силы тяжести. Построена планетарная цифровая модель рельефа с регулярной сеткой задания высот по широте и долготе. Показано, что поправка за дальнюю область в точке наблюдений на физической поверхности Земли имеет значимую величину, носит нерегулярный характер, а интерполяция её значений в плане приводит к неудовлетворительным результатам. Предложено эту поправку разделить на две составляющие - топографическую поправку на поверхности отсчетного эллипсоида и её вертикальный градиент. Эти параметры интегральны по сути и в плане изменяются плавно, что дает возможность их линейной интерполяции.

В четвертом разделе «Технология реализации полной топографической редукции» изложена реализация технологии учета влияния масс промежуточного слоя всей Земли на вертикальную составляющую ускорения силы тяжести.

В пятом разделе «Результаты реализации топографической редукции для дальней области на территорию Сибири» получены значения интегральных параметров в узлах регулярной сетки номенклатурных трапеций с шагом, равным 10 угловым минутам для территории Сибири (от 60 до 120 градусов восточной долготы и от 44 до 76 градусов северной широты). Составлены таблицы этих значений, которые могут быть использованы для линейной интерполяции с погрешностью, не превосходящей погрешности измерения силы тяжести. Приведены обзорные карты изоаномал.

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

СОДЕРЖАНИЕ ЗАЩИЩАЕМЫХ ПОЛОЖЕНИЙ .

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

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

Измерения силы тяжести выполняются на физической поверхности Земли, во внешнем пространстве, под водой и под земной поверхностью, например, в скважинах и шахтах. Выполненные в столь разных условиях наблюдения несопоставимы между собой и с нормальным гравитационным полем, построенным для фигуры относимости (уровенный эллипсоид). Поэтому из измеренных значений силы тяжести вычитают некоторые нормальные величины, представляющие гравитационный эффект теоретической Земли. Такие разности ^g = g-"1 получили название аномалий силы тяжести. Здесь £ - измеренное значение силы тяжести в точке на земной поверхности; у- вычисленное (нормальное) значение силы тяжести в этой же точке. При этом если в расчетах у используются нормальные высоты, то аномалии называют «смешанными», если геодезические (высоты точек наблюдения относительно эллипсоида) - то «чистыми».

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

В общем случае, редукции силы тяжести сводятся к введению ряда поправок в нормальную часть силы тяжести:

где

gA - измеренное значение силы тяжести в точке наблюдений «А»; у0 - нормальное значение силы тяжести на поверхности уровенного эллипсоида (нормальной Земли);

Аул - поправка к нормальному полю на этом же эллипсоиде за высоту его пересчета в точку наблюдений;

8gA - поправка за «нормальный» промежуточный слой с постоянным значением плотности ст0, ограниченный сверху физической поверхностью Земли, снизу - уровнем относимости 7о (в редукции за свободный воздух она равна нулю);

- поправка за плоский (в редукции Буге), сферический или сферои-дический (в топографической редукции) параллельный слой;

5Й - поправка за рельеф.

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

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

аЕб = - (Уо - + ), (3)

называют аномалиями Буге.

Здесь поправка

5гс=271/о0Я,-5&' (4)

отражает притяжение материальной точки «А» плоской пластиной бесконечного простирания (пластиной Буге) толщиной НА, уточненное поправкой за рельеФ %Г ■ в выражении (4) величина/ есть гравитационная постоянная.

Аномалии Буге применяют преимущественно для геологической интерпретации результатов гравиметрических съемок. В геодезии они широко ис-

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

Если обозначить через р* радиус области учета гравитационного влияния промежуточного слоя вокруг точки измерения силы тяжести, то при условии НА « р^ соотношение (3) представляется в виде

(5)

где

=2к/а0НА(Рк +НА -л/р* +Н2Л). (6)

Здесь поправку за рельеф следует рассчитывать в «плоском» варианте. Отличие формулы (5) от выражения (3) заключается в том, что здесь гравитационный эффект плоскопараллельного слоя бесконечного простирания заменен на эффект «плоского» промежуточного слоя, когда уровенная поверхность относимости нормального поля считается плоскостью.

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

где

5«г =5гг[1 + яп|^-0,00012Я^-8я?. (8)

Теперь уже поправка за рельеф 8gf должна вычисляться с учетом сферичности Земли.

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

Принято мнение, что при выполнении редукции гравитационного поля Земли в области, ограниченной радиусом р* круговой зоны до 200 км сферичность Земли можно не учитывать. Однако, как показали вычисления, если выполнять редукцию с точностью, соизмеримой с точностью измерения силы тяжести примерно от 0,01 до 0,03 мГал, то сферичностью Земли можно пренебречь для На — 5 км только в зоне с р* < 100 м, а при НА = 2 км - в зоне с радиусом p/t < 1 000 м. Для pt = 200 км разность плоской и сферической поправок на два-три порядка превосходит погрешность измерений. Очевидно, что более «качественной» является топографическая редукция.

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

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

Внутренней областью считается участок местности расположенный в непосредственной близости от точки вычисления поправки за рельеф. Внутренняя область может быть ограничена площадью круга. При этом её центр (результативная точка «А») совмещается с конкретным пунктом гравиметрических измерений. Размер области в плане определяется или величиной радиуса р„.

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

Дальняя область представляет всю поверхность Земли за исключением внутренней и внешней. Для вычисления поправок за рельеф в дальней области требуется построение планетарной цифровой модели рельефа (ПЦМР). При этом использование эллипсоидальной системы координат здесь является наиболее предпочтительным.

Таким образом, полная топографическая поправка за гравитационное влияние масс промежуточного слоя будет складываться из поправок за внутреннюю внешнюю 5§ц и дальнюю §£Ш области, то есть

= ^ + + ■ (9)

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

При определении поправок за рельеф центральной зоны рекомендуется использовать алгоритм Смирнова В. П., позволяющий решать задачу практически с любой точностью. Здесь используется цилиндрическая система координат. В этой системе учитываемая область делится концентрическими окружностями, проведенными вокруг точки наблюдений «А», на зоны, которые, в свою 16

Внутренняя область; \ Внешний область; И Дальняя область.

Рисунок 1 - Внутренняя, внешняя и дальняя области при вычислении полной топографической редукции

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

Поверхность рельефа каждого сектора в промежутке между двумя соседними точками задания высот местности по «лучу» является конической. В рассматриваемом случае рабочей формулой служит выражение

№

(10)

где к] -значение превышения на расстоянии р; отточки наблюдений;

т - число точек определения высот на луче; I/- количество лучей;

А =А + (о -о ^ 1 .

"т г га-1) 1

Рт+1 — Рт-1

уГ» А о о при

И^'°'Ру'Ро) + Ч0'Лу+1'Ро'Ру+1) при (АД+1)<0;

у+1 н] 2 К/О0

К=2л/о0

+ —7=1П

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

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

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

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

(р., р*) ■= а (2НЛЩ - Я2)+Яр,, рк )Н, (11)

где а = 7С/о0 ———; ^(Р„,Р*) = 271/о0 ;

Р»Р*

м '=1 1_.

Р*Р,

п т

№

рл ,рл - величины внутреннего и внешнего радиусову'-й зоны;

Я,- j - высота рельефа в г'-м секторе;

Я - среднее значение радиуса Земли;

0О - постоянное значение плотности промежуточного слоя.

Аналогично записывается приближенное выражение для поправки за рельеф:

5^-(Р„Рь) = а(нгА -2НАН0 +Й-!) + ^(р„,р,)(Я,-Я) (12) и поправка за параллельный слой

= (13)

Видно, что правые части соотношений (11)—(12) содержат две составляющие - плоскую и сферическую.

Сферическая составляющая при определении поправки за промежуточный слой (формула (11)) не зависит от высоты ЯА результативной точки. Значения Я являются, по существу, в каждой точке земной поверхности средневзвешенной величиной высоты рельефа местности. Их величины очень плавно изменяются в плане, не коррелируя с высотами точек гравиметрических измерений. Эффект, связанный с этим фактором, носит сугубо региональный характер, поэтому при разделении аномалий он исключается с региональным фоном.

Плоская составляющая для поправки за рельеф (выражение (12)) может быть представлена в виде суммы двух слагаемых:

а (я2 - 2Я,Я0 + Я2) = а{НА - Я0 )2 + а(Я2 - Я2) = + . (14)

Первое слагаемое (5#') равняется 0, когда НА-Н0 и возрастает по параболическому закону, если разница (НА-Н0), независимо от её знака, увеличивается. Это указывает на то, что значения Н0 соответствуют интегральной для учитываемой области поверхности, на которой поправки за рельеф минимальны и соответствуют плавно изменяющимся в плане значениям 8§0. Вместе с тем, очевидно, что значение 8g0 дисперсионно отражает степень изрезанности рельефа в учитываемой области вокруг каждой результативной точки. Это говорит о недопустимости весьма трудоемкой процедуры осреднения высот рельефа местности в пределах элементарных площадок вокруг узловых точек задания исходной цифровой модели рельефа.

Значения #0, Н2 и как и значения Я, интегральны (изменяются линейно в плане на заданном интервале) по своей сути. Имеется возможность их определения на конкретную территорию заранее и один раз (если в рельефе не произошло существенных изменений) с построением в графическом или электронном виде соответствующих карт. Величины бя'также заранее могут быть представлены аналитически в виде таблиц или электронном. Тогда по мере проведения гравиметрических съемок с получением плановых координат и высот точек измерений можно сразу же получать поправки за рельеф и промежуточный слой в целом, как в редукции Буге, так и в неполной топографической редукции.

Для вычисления отмеченных выше интегральных параметров используют различного вида ЦМР. Пусть одна из них представляется вертикальными параллелепипедами, которые имеют горизонтальные основания в виде квадрата (рисунок 2).

Высоты ЦМР здесь заданы в узлах квадратной сетки с шагом Д, отражающим размер параллелепипеда в плане. Каждая зона здесь состоит из ряда поясов, минимальный г и максимальный и, номера которых определяют ее положение относительно результативной точки.

Тогда величины а, Н0 и Н определяются из соотношений

Результати и пая точка

Внешняя граница учитываемой -зоны

гп ¡=п )-п

^ Ы-п }--п

1=П J-'^

ЕЕ^

нп =

_ і=-п У=-п

Г=П у-'»

ЕЕ^-

І—П ]=П

_ ЕЕ^

ці _ *=-п у=-я_

|=В у-"

ЕЕ^-

ГЗ N Л.. ' 5

'' і

і'й:ч!.

Внутренняя граница учитываемой зоны

Номера поясов зоны

Рисунок 2 - Фрагмент задания ЦМР в узлах квадратной сетки

где Р,= —

Vа? +Р2 У^ГЙІГ ,>/«,*+ Р? д/^Тр «іРі «.Р. а2р! а2|32

Ну - отметки высот в г',у-м узле зоны.

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

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

Если задана отдельная (і-я) трапеция, ограниченная двумя меридианами с долготами Ь\, Ь2, и двумя параллелями с геодезическими широтами Ви В2, то вертикальная составляющая притяжения масс соответствующей призмы высо-

той Н, в предположении, что направление отвесной линии совпадает с направлением вектора , определяется следующим соотношением:

(R-RjRjAfiVcosfl

i А о

dBdLdh,

(15)

где R =

(N + h)cosBcosL (N + h)cosBsinL [(l-e2)W + /i]sinB

; R,=

(Na+Ha)cosBacosLa (NA+HA)cosBAsmLA [(l-e^+H^sinS,

N, Na - радиусы кривизны первого вертикала в точках тиА соответст-

венно;

М- радиус кривизны меридиана в точке т;

е2 - квадрат первого эксцентриситета земного эллипсоида.

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

RcRJ-FI

(i6)

где V/ - приближенный объем притягивающих масс сфероидической призмы, определяемый соотношением V, =МС (В2 -5,)--соз5с ■(Ь1-Ь1)-Н;

R с =

(Nc +H/2)cosBc cosLc (Nc + Я/2) cos Bc sin Lc [(l-e2)jVc +H¡2]sm Bc

Nc, Mc -радиусы кривизны первого вертикала и меридиана в точке С.

Пусть 2 - общее количество призм с высотами Д (/ = 1, 2, 3, ..., 0 в дальней области. Тогда суммарное влияние топографических масс промежуточного слоя на значение силы тяжести в точке наблюдений «А» определится соотношением

8*ш=26&(# = я;)- (17)

На основе формулы (17) были рассчитаны значения поправок 5я/;; в трех точках земной поверхности с широтами ВА = 40, 50 и 60°, лежащих на меридиане с долготой, равной 80°. Массы промежуточного слоя учитывались, начиная с расстояния р* от результативной точки, принимающего значения от 200 до 20 000 километров. Результаты вычислений по каждому ряду трех результативных точек приведены в таблице 1.

Таблица 1 - Значения поправки в зависимости от радиуса р* «неучтенной» внешней области промежуточного слоя

Ьл (градусы) 80

ВА (градусы) 40 50 60 Р*

Нл (км) 1,060 0,330 0,060 (км)

48,385 39,698 32,895 200

43,522 38,270 32,617 500

33,331 34,860 31,745 1 000

20,959 24,800 27,368 2 000

(мГал) 16,339 16,039 21,602 3 000

13,694 13,620 14,507 4 000

12,119 12,053 12,193 5 000

6,435 6,076 5,306 10 000

1,463 3,146 3,242 15 000

0 0 0 20 000

Как видно из таблицы 1, массы, расположенные за пределами зоны радиуса даже в 10 000 километров, оказывают ощутимое переменное по величине гравитационное влияние. По данным таблицы 1 построен график изменения ££/// для исследуемых точек (рисунок 3).

Рисунок 3 - Изменение значения поправки с изменением радиуса р* неучтенной области промежуточного слоя

Из рисунка 3 видно, что построенные кривые имеют нерегулярный характер.

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

меньше (таблица 2).

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

Таблица 2 - Значения поправки (в зависимости от радиуса р4) с шагом по широте, равным 1 градусу

LA (градусы) 80 Р» (км)

ВА (градусы) 41 42 43 44 45 46

НЛ (км) 1,420 5,780 1,965 0,530 3,120 0,500

ъ3и, (мГал) 48,606 51,315 47,038 44,815 45,933 43,252 200

43,801 44,543 42,611 41,488 41,503 40,486 500

33,821 34,691 34,778 35,139 35,768 35,555 1 000

20,872 20,827 20,780 20,943 21,426 21,948 2 000

16,277 16,217 16,147 16,082 16,014 15,950 3 000

13,655 13,633 13,609 13,595 13,590 13,594 4 000

12,114 12,101 12,083 12,067 12,048 12,028 5 000

6,401 6,363 6,333 6,298 6,264 6,234 10 000

1,515 1,637 1,821 | 2,023 2,219 2,434 15 000

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

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

На основе анализа множества цифровых моделей рельефа для построения ПЦМР в качестве основы была выбрана модель радарной топографической съемки Shuttle Radar Topography Mission (SRTM), покрывающая около 80 % всей поверхности суши земного шара. В этой модели значения высот заданы в узлах регулярной сетки с шагом в 3 угловые секунды. Декларируемая погрешность высот модели соответствует величине ±16 м. Данные модели SRTM имеют свободный доступ в Интернете.

На основании этой модели были сформированы блоки (трапеции) размером 15°х15° по широте (В) и долготе (L) с шагом значений высот в 5 угловых минут в интервале 45° <В < 60°; 30° <L < 180°. Для остальной территории шаг составил Г (рисунок 4).

Данные SRTM в нашем варианте занимают 192 блока размером 15°х15°.

Высоты для двадцати трапеций, лежащих в интервале 60° <Б < 90°; 30° <L < 180°, получены с топографических карт масштаба 1 : 1 ООО ООО. Остальная информация - с физических карт более мелкого масштаба.

Таким образом, сформированная из 288 блоков планетарная цифровая модель рельефа (ПЦМР-СГГА), покрывает всю поверхность Земли. Следует отметить, что высоты в узлах сетки, попадающие на водные пространства, приняты равными нулю. Общий объем масс промежуточного слоя, согласно ПЦМР-

СГГА, составляет 1,175 115 ■ 1017 м3.

Исследование влияния систематических и случайных погрешностей высот разработанной ПЦМР на значение топографической поправки за дальнюю область (при р4 = 200 км) привело к следующим результатам:

- изменение всех высот модели на 1 м приводит к изменению топографической поправки примерно на 0,047 мГал, при этом максимальное отличие поправок от их среднего значения составляет 0,002 мГал;

- случайная погрешность высот модели в 1 м приводит к погрешности топографической поправки, равной 0,0005 мГал.

Следовательно, декларированная погрешность (16 м) высот модели БЯТМ, которая является основой модели ПЦМР-СГГА, вполне приемлема.

Результаты расчета значений топографической поправки для всей дальней области по формуле (17) с использованием построенной планетарной ЦМР показали следующее:

- процесс вычисления на современном персональном компьютере значений поправок весьма трудоемкий;

- линейная интерполяция значений поправок между результативными точками, расположенными через 10 угловых минут, отличается от непосредственно определяемой величины в пределах ±0,7 мГал, когда высоты результативных точек изменяются от первых сотен метров до двух километров;

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

Выяснилось, что изменение величины топографической поправки 5&п в зависимости от высоты Нл связано со значением вертикального градиента Сш. Определение значений величин вт в равномерно расположенных узловых точках ПЦМР по широте и долготе производилось с учетом разницы топографических поправок на земной поверхности 8§1П (высота точки Ни) и на поверхности принятого эллипсоида 5^®, (Ни = 0)

(ни )--77-1-• (18)

п>,1

Здесь индекс «г» определяет номер узла по широте, <</» - по

долготе.

G (мГал/км)

25 1 —-----------

20 •

15 •

10 •

5 • 0 - 1

100 150 200 250

?к <КМ)

Рисунок 5 - График изменения величины Gui в зависимости от значения р*

На рисунке 5 приведен график изменения параметров вщ в зависимости от значения р* в результативной точке с геодезическими координатами В = 35°,

Ь = 90°.

Из рисунка 5 видно, что с увеличением внутреннего радиуса дальней области значение градиента уменьшается по экспоненте.

Из рассмотрения графика видно, что изменение значений С(р, ,В,НА)

весьма дифференцировано по величине при р* < 150 км, особенно в зависимости от высоты ни. Максимум значений, соответствующий величине Н0, как уже отмечалось, по мере увеличения величины р* смещаются к нулю, переходя, по-видимому, к отрицательному значению.

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

рисунке 6.

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

20 40 60

Знамения широты в градусах

Рисунок 6 - Обобщенный график изменения величины G с широтой

ально существующих высот при величине р* = 200 км показали, что значение па-

раметра <7(р,, В, Нл) изменяется практически линейно по вертикали (всего 0,004 мГал/км) и весьма плавно в пределах от 0 до 2,030 мГал/км в плане на территории России и в ближайшем её окружении. Это свидетельствует о том, что параметр <?(р,,В,НА) при р, > 200 км является интегральной величиной.

Величина топографической поправки на поверхности эллипсоида зависит как от радиуса р* внутренней границы дальней области, так и от мощности промежуточного слоя в ней вокруг района исследования. Возьмем в качестве примера точку с координатами ВА = 45°, ЬА = 80° и для различных значений величины р* рассчитаем величины поправок (таблица 3).

Таблица 3 - Изменение значения топографической поправки на поверхности эллипсоида в зависимости от радиуса рк

р*> 8*4,

км мГал

25 35,213

50 43,033

100 44,192

150 44,042

200 43,794

500 40,766

1 000 35,442

40 •

0 200 400 600 800 1000 Радиус неучтенной области, в км

Как видно из результатов таблицы 3, наибольшее значение поправки соответствует величине р* равной 100 км. С уменьшением радиуса рь поправка резко убывает и, как показали дополнительные вычисления, при р* < 10 км переходит в отрицательные значения. Рост радиуса также приводит к уменьшению топографической поправки, которая при р4 = пКср, где Яср - средний радиус Земли, обращается в ноль. При этом, начиная с радиуса рь равного 200 км, функция б^щ(р^) почти линейна.

Предварительный анализ показал, что величины и §£ш линейны в плане (с погрешностью, не превосходящей погрешность измерения силы тяжести) на интервале 10 угловых минут. Для подтверждения этого были выполнены расчеты параметров <?ш и по меридианам с долготами ЬА, равными 70, 80, 90 и 100 градусов (профили 1, 2, 3 и 4 соответственно) в широтном интервале [40°, 60°] с шагом в 5 угловых минут. При этом р* полагался равным 200 км.

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

и определялись разности

V, = (Я,)---

От(Вм) + Ош(В„)

Очевидно, что величины V, и в, характеризуют погрешности линейной интерполяции рассчитанных значений Ъ8°Ш{В,) и СШ(В{) соответственно в интервале 10 угловых минут. В таблице 4 приведены средние квадратические погрешности (1^0 и )1С, вычисленные по значениям V, и .

Таблица 4 - Средние квадратические погрешности интерполяции значений и °т(Ю по заданным профилям

№ М-6г«>

профиля мГал мГал/км

1 0,0008 0,0007

2 0,0014 0,0010

3 0,0010 0,0006

4 0,0011 0,0006

Из таблицы 4 видно, что средние квадратические погрешности интерполяции по всем профилям не превосходят 0,0015 мГал для и 0,001 мГал/км -для параметра (?ш.

Таким образом, линейная интерполяция величин бо-Ц, и С,,, допустима для узлов регулярной сетки с шагом в 10 угловых минут при рк = 200 км. При этом топографическая поправка 5gш за влияние промежуточного слоя дальней области может быть рассчитана по формуле

^ш^СЛНС^Я,. (19)

Здесь ^щ(А) и (7Ш(А) - интерполированные к результативной точке значения интегральных параметров, определенных в узлах регулярной сетки номенклатурных трапеций; НА - отметка высоты результативной точки.

Как уже отмечалось, значение полной топографической поправки 8^ за гравитационное влияние масс всего промежуточного слоя складывается из поправок за внутреннюю (5^,), внешнюю (<%,/) и дальнюю (<%/я ) области (формула (9)).

Обобщенная технологическая схема определения аномалий силы тяжести в полной топографической редукции представлена на рисунке 7.

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

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

Рисунок 7 - Обобщенная технологическая схема определения аномалий силы тяжести в полной топографической редукции

Учет гравитационного влияния промежуточного слоя в центральной зоне во многих случаях требует индивидуального подхода, который сложно автоматизировать. Наиболее рациональным, малозатратным является корреляционный способ, когда имеется возможность использовать отметки высот соседних точек гравиметрических измерений. Топографическая поправка 5*. за промежуточный слой внутренней области может быть представлена суммой двух поправок, описанных выражениями (6) и (10).

Основой для реализации топографической редукции ео внешней области служат цифровые модели рельефа на заданную область и соотношения (11)-<14). Здесь значения Я0, 5g0, как и величина Я интегральны по своей сути. Эти параметры можно определить заранее и один раз (если в рельефе не произошло существенных изменений) с построением в графическом или электронном виде соответствующих карт и таблицы значений 5g' на конкретную территорию. Тогда по мере проведения гравиметрических съемок с получением плановых координат и высот точек измерений можно сразу же получать поправки за рельеф и промежуточный слой в целом, как в редукции Буге, так и в неполной топографической редукции. При расчете поправки 5gn внешняя область должна быть ограничена радиусом рь равным 200 км.

Учет гравитационного влияния промежуточного слоя * дальней области осуществляется с использованием ПЦМР. Вычисление поправки SgDI осуществляется по формуле (19). Суммарная топографическая поправка в каждой результативной точке определяется из соотношения (9), на основании которого рассчитываются аномалии силы тяжести в полной топографической редукции (выражение (1)).

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

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

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

координат точек наблюдений.

К настоящему время выполнен огромный объем производственных работ по определению аномалий силы тяжести, в которых учитываются поправки за рельеф в зонах внутренней и внешней области (в редукции Буге и в неполной топографической редукции). Согласно Инструкции по гравиразведке, радиус р* внешней области при этом должен был приниматься равным 200 км. Поэтому будем считать, что поправки 5g, и 5g„ в выражении (9) определены на всю

территорию России.

Ранее было показано, что параметры Sg°, и СП1, составляющие общую

поправку 5gni за дальнюю область при р* - 200 км, линейны в пределах 10 угловых минут по широте и долготе. Это дает возможность составить электронные таблицы значений отмеченных величин в узлах регулярной сетки номенклатурных трапеций с целью их дальнейшей интерполяции к результативным

точкам.

Для вычисления поправки Sg,n на территорию Западной и Восточной Сибири согласно формуле (19), рассчитаны таблицы значений параметров 5g°m и в номенклатурной разграфке карт масштаба! : 1 000 000, представленных

в виде документа Microsoft Word.

До 60-й параллели таблицы интегральных параметров рассчитаны с шагом

в 10 угловых минут по широте и долготе. Далее (при широте В >60°) шаг по

долготе составляет 20 угловых минут.

Для примера в таблицах 5 и 6 представлены фрагменты номенклатурной

трапеции L-44.

Таблица 5 - Значения параметра (мГал) в узлах регулярной сетки фрагмента номенклатурной трапеции Ь-44

B\L 78 10 20 30 40 50 79 10

48 40,587 40,644 40,699 40,754 40,808 40,861 40,913 40,965

50 40,743 40,799 40,854 40,909 40,962 41,015 41,067 41,118

40 40,902 40,958 41,013 41,067 41,120 41,172 41,224 41,274

30 41,063 41,119 41,173 41,226 41,279 41,330 41,380 41,431

20 41,226 41,281 41,335 41,387 41,438 41,488 41,538 41,587

10 41,389 41,443 41,496 41,547 41,597 41,645 41,693 41,739

47 41,548 41,602 41,654 41,704 41,751 41,798 41,843 41,888

50 41,702 41,753 41,803 41,851 41,898 41,942 41,986 42,029

Таблицы 5, 6 документа Word предназначены для непосредственной (табличной) интерполяции значений интегральных величин 5 gf„ и Gm.

По результатам вычислений составлены обобщенные карты изолиний интегральных параметров на территорию Сибири (рисунки 8 и 9).

Таблица 6 - Значения параметра <7га (мГал/км) в узлах регулярной сетки фрагмента номенклатурной трапеции Ь-44

В \ L 78 10 20 30 40 50 79 10

48 0,398 0,402 0,406 0,409 0,413 0,416 0420 0,423

50 0,408 0,412 0,416 0,420 0,424 0,428 0,431 0,435

40 0,419 0,424 0,428 0,432 0,436 0,441 0,444 0,448

30 0,432 0,437 0,441 0,446 0,450 0,454 0,458 0,462

20 0,446 0,451 0,456 0,460 0,465 0,469 0,474 0,478

10 0,461 0,466 0,471 0,476 0,481 0,485 0,490 0,494

47 0,476 0,481 0,487 0,493 0,497 0,502 0,506 0,510

50 0,490 0,495 0,501 0,506 0,511 0,516 0,519 0,523

Рисунок 8 - Обобщенная карта изолиний параметра 8g0ш

Рисунок 9 - Обобщенная карта изолиний параметра вш

'Єкатеринбург? ♦Челябинск

'Омск

Из-за ограничений размеров страницы, изолинии значений и Gm проведены с сечением 1 мГал и 0,05 мГал/км соответственно. Фактически же на конкретную номенклатурную трапецию они могут быть представлены с любым сечением, позволяющим выполнить интерполяцию с заданной точностью.

Из рисунков видно, что на величину вертикального градиента поправки 5Яш (рисунок 9) существенное влияние оказывают массы горных массивов Средней и Южной Азии (Тянь-Шань, Памир, Куньлунь и др.). В меньшей степени - Уральские горы и горная система Алтае-Саянской складчатой области.

Поправка 5^t°„ (рисунок 8) более регулярна и менее чувствительна к «мелким» горным хребтам. Основное влияние на неё оказывают массы, расположенные в южной части Азии (Куньлунь, Тибет).

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

Следует отметить, что разработанный комплекс вычислительных программ в среде Delphi на языке Object Pascal позволяет построить аналогичные таблицы и карты на всю территорию России, а при необходимости - и на всю поверхность Земли.

Таким образом, на основании полученных результатов можно сделать следующие выводы и рекомендации:

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

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

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

• для учета влияния рельефа области, ограниченной радиусом в 200 км, рекомендуется использовать методику, разработанную в СНИИГГиМСе, позволяющую, с одной стороны, существенно упорядочить процесс полевых измерений и проводить заранее значительный объем камеральных работ, с другой, - выбрать с требуемой точностью оптимальные условия учета влияния топографических масс.

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

ЗАКЛЮЧЕНИЕ

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

значение в развитии наук о Земле.

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

Основные научные и практические результаты, полученные в диссертации,

заключаются в следующем.

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

правки за дальнюю область, которая является необходимой составляющей полной топографической поправки, постоянно возрастает с увеличением размеров учитываемой области и носит нерегулярный характер. Так, например, поправка только за гравитационное влияние масс горных массивов Мексики и западной части США на точку, расположенную на меридиане с долготой 80°, на широте 45° (Восточный Казахстан) составляет 2,22 мГал, а на широте 46° этого же меридиана - уже 2,43 мГал.

2. Создана цифровая модель рельефа поверхности всей Земли, не имеющая аналогов, в которой значения геодезических высот представлены в узлах регулярной сетки сфероидических трапеций с погрешностью высот, не превосходящих 20-50 м на территории России. На основе статистического анализа установлено, что случайная погрешность высот построенной планетарной цифровой модели рельефа в 1 м приводит к погрешности расчета топографической поправки за дальнюю область, составляющую 0,0005 мГал, что подтверждает достаточную точность цифровой модели при расчете соответствующей поправки с погрешностью, соизмеримой с погрешностью измерений (0,01-0,03 мГал).

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

4. Обосновано, что топографическая поправка за влияние промежуточного слоя дальней области может быть разделена на две интегральные составляющие: 1 - поправку на поверхности отсчетного эллипсоида 5g", и 2 - поправку за высоту точки наблюдений с учетом её вертикального градиента GUI.

5. Экспериментально доказано, что параметры 5^°, и Gm линейны в интервале 10 угловых минут по широте и долготе, что дает возможность получения значений поправок за дальнюю область с погрешностью, не превосходящей погрешность гравиметрических измерений.

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

7. На основе разработанной технологии созданы цифровые модели величин и Gm для территории Западной и Восточной Сибири, позволяющие достоверно получать значения топографических поправок за дальнюю область, которые являются составной частью полной топографической редукции.

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

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

1 Дементьев, Ю. В. Подтверждение к алгоритму 1336 [Текст] / Ю. В. Дементьев, Н.К. Шендрик // Сб. Библиотека алгоритмов 1516 - 2006. Под ред. М. И. Агеева. - М.: Радио и связь, 1981. - С. 176.

2 Дементьев, Ю. В. Изменение ускорения силы тяжести под действием приливных деформаций Земли и атмосферных нагрузок [Текст] / Ю. В. Дементьев И Сб. материалов Ш Международного научного конгресса «ГЕО-Сибирь-2007». - Том 1, ч. 2. - Новосибирск: СГГА, 2007. - С. 266-272.

3 Дементьев, Ю. В. Математическое моделирование движения точек земной поверхности [Текст] / Ю. В. Дементьев // Известия вузов. Горный журнал. -2007.-№5.-С. 81-85.

4 Дементьев, Ю. В. Учет сжатия Земли в топографических редукциях силы тяжести [Текст] / Ю. В. Дементьев // Геодезия и картография. - 2008. - № 2. -С. 3-6.

5 Дементьев, Ю. В. Расчет топографических редукций силы тяжести по съемочным трапециям земного эллипсоида [Текст] / Ю. В. Дементьев // Геодезия и картография. - 2008. - № 7. - С. 14-16.

6 Дементьев, Ю. В. Изменение ускорения силы тяжести под действием атмосферного давления [Текст] / Ю. В. Дементьев // Сборник материалов V Международного научного конгресса «ГЕО-Сибирь-2009». - Новосибирск, 2009. -С. 80-81.

7 Дементьев, Ю. В. Численно-аналитический метод расчета топографической редукции в центральной и ближней зонах [Текст] / Ю. В. Дементьев // Геодезия и картография. - 2009. - № 4. - С. 15-17.

8 Дементьев, Ю. В. Об уточнении формул для вычисления вертикальной

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

[Текст] / Ю. В. Дементьев, А. И. Каленицкий // Геодезия и картография. - 2010. -№ 3. - С. 7-9.

9 Дементьев, Ю. В. Связь плоских, сферических и сфероидических поправок за гравитационное влияние промежуточного слоя постоянной плотности и толщины [Текст] / Ю. В. Дементьев // Геодезия и картография. - 2010. - № 4. -С. 19-21.

10 Дементьев, Ю. В. Учет влияния топографических масс при вычислении возмущающего потенциала [Текст] / Ю. В. Дементьев// Вестник СГГА, вып. 1(12). - Новосибирск: СГГА, 2010. - С. 44-48.

11 Дементьев, Ю. В. О редукциях силы тяжести [Текст] / Ю. В. Дементьев // Сборник материалов VI Международного научного конгресса «ГЕО-Сибирь-2010». - Новосибирск, 2010. - С. 158-163.

12 Дементьев, Ю. В. Построение планетарной цифровой модели рельефа и её приложения [Текст] / Ю. В. Дементьев, Е. Н. Кулик, Е. В. Дергачева // Сбор-

ник материалов VI Международного научного конгресса «ГЕО-Сибирь-2010». -

Новосибирск, 2010. - С. 170-173.

13 Дементьев, Ю. В. Зависимость поправок за плоский и сферический слои в неполной топографической редукции от толщины слоя и радиуса учитываемой зоны [Текст] / Ю. В. Дементьев // Вестник СГГА, вып. 2(13). - Новосибирск: СГГА, 2010. - С. 13-17.

14 Дементьев, Ю. В. Построение планетарной цифровой модели рельефа Земли для выполнения полной топографической редукции гравитационного поля [Текст] / Ю. В. Дементьев, А. И. Каленицкий, Е. Н. Кулик, А. В. Чере-мушкин // Геодезия и картография. - 2010. - № 12. - С. 17-19.

15 Дементьев, Ю. В. Оптимизация планетарной цифровой модели рельефа Земли для выполнения полной топографической редукции гравитационного поля [Текст] / Ю. В. Дементьев, Е. Н. Кулик, Е. С. Акулич // Сборник материалов VII Международного научного конгресса «ГЕО-Сибирь-2011». - Новосибирск,

2011.-С. 202-204.

16 Дементьев, Ю. В. Рекомендации по выбору оптимальной модели рельефа для гравиметрии на территории Российской Федерации [Текст] / Ю. В. Дементьев, А.И. Каленицкий, E.H. Кулик // Геодезия и картография. - 2011. -Na 7.-С. 9-13.

17 Дементьев, Ю. В. О величине и значимости поправки за гравитационное влияние масс промежуточного слоя внешней области [Текст] / Ю. В. Дементьев // Геодезия и картография. - 2011. - № 9. - С. 8-10.

18 Дементьев, Ю. В. Возможность интерполяции поправок, обусловленных влиянием масс промежуточного слоя внешней области при топографической редукции силы тяжести [Текст] / Ю. В. Дементьев // Геодезия и картография. -2011.-№10.-С. 2-6.

Содержание диссертации, доктора технических наук, Дементьев, Юрий Викторович

ВВЕДЕНИЕ.

1 АНОМАЛИИ СИЛЫ ТЯЖЕСТИ И ИХ ПРИМЕНЕНИЕ В ГЕОДЕЗИИ

И РАЗВЕДОЧНОЙ ГЕОФИЗИКЕ (АНАЛИТИЧЕСКИЙ ОБЗОР).

1.1 Понятие об аномалиях силы тяжести и их редукциях.

1.2 Аномалии в свободном воздухе. Поправка за высоту.

1.3 Поправка за плоский параллельный слой. Аномалия Буге.

1.4 Поправки за промежуточный слой. Топографические редукции.

1.5 Поправки за рельеф.

1.6 Уточнение формул для вычисления вертикальной составляющей величины притяжения масс сферического или плоского слоя.

1.7 Зависимость поправок за плоский и сферический слои от их толщины и радиуса учитываемой зоны.

1.8 Применение аномалий силы тяжести в геодезии.

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

1.10 Учет изменения ускорения силы тяжести со временем.

1.10.1 Постановка задачи.

1.10.2 Учет приливных деформаций Земли.

1.10.3 Поправки за атмосферные нагрузки.

2 МЕТОДИКА УЧЕТА ВЛИЯНИЯ ПРОМЕЖУТОЧНОГО СЛОЯ ПРИ ВЫПОЛНЕНИИ ПОЛНОЙ ТОПОГРАФИЧЕСКОЙ РЕДУКЦИИ.

2.1 Общие положения.

2.2 Поправка за рельеф во внутренней области.

2.3 Топографическая поправка во внешней области.

2.4 Поправка за топографию в дальней области.

2.5 Связь плоских, сферических и сфероидических поправок за гравитационное влияние параллельного слоя.

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

3.1 Общие положения.

3.2 Построение планетарной цифровой модели рельефа Земли.

3.3 Обоснование оптимальной ЦМР для территории Российской Федерации.

3.4 Влияние случайных и систематических погрешностей высот ПЦМР на значение топографической поправки.

3.5 Предельная погрешность задания внутренней границы дальней области.

3.6 Влияние масс промежуточного слоя дальней области на значение топографической поправки.

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

3.7.1 Постановка задачи.

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

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

3.8 Линейная интерполяция в сфероидической трапеции.

4 ТЕХНОЛОГИЯ РЕАЛИЗАЦИИ ПОЛНОЙ ТОПОГРАФИЧЕСКОЙ РЕДУКЦИИ.

5 РЕЗУЛЬТАТЫ РЕАЛИЗАЦИИ ТОПОГРАФИЧЕСКОЙ РЕДУКЦИИ ДЛЯ ДАЛЬНЕЙ ОБЛАСТИ НА ТЕРРИТОРИЮ СИБИРИ.

Введение Диссертация по наукам о земле, на тему "Развитие теории и разработка технологии определения аномалий силы тяжести в полной топографической редукции"

Актуальность темы исследования. Гравиметрия является наукой, изучающей поле силы тяжести Земли и методы её измерения. Измеренные значения силы тяжести в точках земной поверхности используются для определения фигуры и внутреннего строения Земли, выявления причин и выполнения оценок последствий геодинамических процессов природного и техногенного характера, выполнения геологического картирования и тектонического районирования территорий, при поисках и разведке месторождений полезных ископаемых. Раздел гравиметрии, посвященный приложению гравиметрических данных к решению геодезических задач по определению фигуры и гравитационного поля Земли, приведению геодезических измерений на поверхность эллипсоида, называется геодезической гравиметрией. Высока роль в развитии геодезической гравиметрии Центрального ордена «Знак Почета» научно-исследовательского института геодезии, аэросъемки и картографии им. Ф.Н. Красовского (ЦНИИГАиК), что нашло отражение в работах [3,4,27,28,29,30,32,41,44,45,46,47, 52,53,54]. Весомый вклад в решение задач геодезической гравиметрии внесли советские и российские ученые: Бровар В. В., Бровар Б. В., Бузук В. В., Вовк И. Г., Грушинский Н. П., Демьянов Г. В., Еремеев В. Ф., Каленицкий А. И., Кузьмин Ю.О., Магницкий В. А., Мазурова Е. М., Машимов М. М., Михайлов А. А., Молоденский М. С., Огородова Л. В., Остач О. М., Пеллинен Л. П., Шимбирев Б. П., Юркина М. И. и др., а также зарубежные ученые Гофман-Велленгоф Б., Мориц Г.

В настоящее время гравиметрическая аппаратура позволяет измерять о силу тяжести с погрешностью (0,01 - 0,03) мГал, или (1 - 3)х10 относительно абсолютного значения. Методы пространственного координирования позволяют определить положение гравиметрического пункта и значение нормального гравитационного поля с такой же точностью. Однако, аномалии силы тяжести, особенно в районах со значительным перепадом высот, определяются в 50 - 100 раз грубее. Поэтому, проблема повышения точности определения аномального гравитационного поля Земли, как на её поверхности, так и вне её, в том числе во внешнем пространстве является актуальной, особенно в последнее время в связи с необходимостью обеспечения более точного решения задач физической геодезии, общей и прикладной геофизики, изучения процессов глобальной, локальной, в том числе техногенной геодинамики, определения начальных параметров запуска и корректировки траектории космических аппаратов, а также целого ряда других, имеющих принципиальное значение во многих отраслях знаний.

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

При вычислении аномалий высот и уклонений отвесной линии требуется знание величины аномалии силы тяжести в редукции за свободный воздух в точке, где измеряются геодезическая высота и сила тяжести. В условиях дискретной гравиметрической съемки для определения аномалий в промежуточных пунктах, где наблюдения силы тяжести не производились, появляется необходимость в интерполировании результатов измерений. Однако аномалии в свободном воздухе существенно зависят от изменения высоты рельефа местности. Как показано в книге [50], при расстояниях между пунктами от 20 до 30 км линейная интерполяция аномалий в свободном воздухе в горных районах приводит к погрешностям порядка (2025) мГал.

При использовании в косвенной интерполяции аномалий силы тяжести в неполной топографической редукции погрешности интерполяции ещё уменьшается. Как показали исследования, при выполнении неполной топографической редукции с учетом влияния топографических масс постоянной плотности, расположенных в промежуточном слое, ограниченном снизу поверхностью относимости нормального поля, сверху -физической поверхностью Земли, для горной области в радиусе 200 км от точки наблюдений погрешность интерполяции составляет от 0,7 до 0,8 мГал. Вместе с тем, установлено [19,25], что значение топографической поправки постоянно возрастает с увеличением размеров учитываемой области, и нет предела, при котором можно «отбросить» влияние масс дальних зон.

В работе [47] показано, что определение высот квазигеоида с погрешностью ± 1 м приводит к погрешности вычисления аномалий силы тяжести порядка ± 0,25 мГал. Очевидно, что для расчета высот квазигеоида с субметровой точностью требуется знание аномалий силы тяжести с погрешностью не ниже 0,02 мГал. Такую точность может обеспечить использование аномалий силы тяжести в полной топографической редукции.

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

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

При выполнении полной топографической редукции силы тяжести предполагается, что весь учитываемый промежуточный слой разделяется на внутреннюю, внешнюю и дальнюю области. Внутренней областью считается участок местности, расположенный в непосредственной близости от точки наблюдений. Размер её колеблется от нескольких метров до нескольких километров в зависимости от изрезанности рельефа. Граница внешней области рекомендована Инструкцией по гравиразведке [33]. Её радиус в среднем составляет 200 км. Дальняя область представляет всю поверхность Земли за исключением внутренней и внешней областей.

В настоящее время накоплен значительный опыт учета влияния рельефа местности во внутренней и внешней областях на значение силы тяжести. Систематизация этого опыта и новая автоматизированная методика выбора оптимального шага цифровой модели рельефа (ЦМР) местности во внешней учитываемой области в зависимости от масштаба и точности гравиметрических работ изложена в Методических рекомендациях по учету влияния рельефа местности, утвержденных МинГео СССР как практическое руководство [39]. Эта методика позволяет определять аномалии силы тяжести как в редукции Буге, так и в неполной топографической редукции. Однако влияние топографических масс, расположенных в дальней области, остается мало изученным. Поэтому отмеченные выше Методические рекомендации требуют дополнения методикой определения аномалий силы тяжести в полной топографической редукции.

Для достижения поставленной цели решались следующие основные задачи.

- разработка технологии реализации полной топографической редукции.

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

Основные положения диссертации, выносимые на защиту:

- обоснованная возможность определения «чистых» аномалий силы тяжести в полной топографической редукции существенно повышает как качество геодезических исследований в изучении фигуры и гравитационного поля Земли, так и эффективность использования гравиметрии в прикладных целях;

Научная новизна исследований:

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

Результаты выполненных исследований внедрены в учебный процесс СГГА (акт о внедрении от 15 декабря 2011 г.).

Заключение Диссертация по теме "Геодезия", Дементьев, Юрий Викторович

ЗАКЛЮЧЕНИЕ

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

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

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

- создана цифровая модель рельефа поверхности всей Земли, не имеющая аналогов, в которой значения геодезических высот представлены в узлах регулярной сетки сфероидических трапеций с погрешностью высот, не превосходящих (20 - 50) м на территории России. На основе статистического анализа установлено, что случайная погрешность высот построенной планетарной цифровой модели рельефа в 1 м приводит к погрешности расчета топографической поправки за дальнюю область, составляющую 0,0005 мГал, что подтверждает достаточную точность цифровой модели при расчете соответствующей поправки с погрешностью, соизмеримой с погрешностью измерений (0,01 - 0,03) мГал\

- эмпирически обосновано, что линейная интерполяция значений топографических поправок за дальнюю область в предгорных и горных районах в интервале 10' приводит к погрешности ±0,7 мГал, что является характеристикой точности при интерполировании аномалий силы тяжести в неполной топографической редукции;

- обосновано, что топографическая поправка за влияние промежуточного слоя дальней области может быть разделена на две интегральные составляющие: поправку на поверхности отсчетного эллипсоида ( 8g(jп) и поправку за высоту точки наблюдений с учетом её вертикального градиента (£///);

- экспериментально доказано, что параметры 8g(|/I и линейны в интервале 10' по широте и долготе, что дает возможность получения значений поправок за дальнюю область с погрешностью, не превосходящую погрешность гравиметрических измерений;

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

- на основе разработанной технологии созданы цифровые модели значений поправок на поверхности отсчетного эллипсоида ( бg^//) и вертикального градиента (С?///) для территории Западной и Восточной Сибири, позволяющие достоверно получать значения топографических поправок за дальнюю область, которые являются составной частью полной топографической редукции;

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

Библиография Диссертация по наукам о земле, доктора технических наук, Дементьев, Юрий Викторович, Новосибирск

1. Антонович, K.M. Использование спутниковых радионавигационных систем в геодезии Текст.: в 2т. Т.1: Монография /K.M. Антонович; ГОУ ВПО «Сибирская государственная геодезическая академия». М.: Картгеоцентр, 2005.-334 е.: ил.

2. Бахвалов, Н.С. Численные методы Текст. / Н.С. Бахвалов. М.: Наука, 1975. - 632 е.: ил.

3. Бровар, B.B. К решению краевой задачи Молоденского сотносительной погрешностью 5-Ю-5 Текст. / В.В. Бровар // Научно-технический сборник по геодезии, аэрокосмическим съемкам и картографии. Физическая геодезия. Кн. 2. М.: ЦНИИГАиК, 1996. - С. 156-164.

4. Голенко, Д-И. Моделирование и статистический анализ псевдослучайных чисел на электронных вычислительных машинах Текст. / Д.И. Голенко М.: Наука, 1965. - 227 с.

5. Гофман-Велленгоф, Б. Физическая геодезия Текст.: учебник / Б. Гофман-Велленгоф, Г. Мориц; пер. с англ. Ю.М. Неймана, Л.С. Сугаиповой. М.: МИИГАиК, 2007, - 426 с.: илл.

6. Грушинский, Н.П. Теория фигуры Земли Текст. / Н.П. Грушинский -М.: Физматгиз, 1963. 448 е.: илл.

7. Двулит, П.Д. Сезонные деформации уровенной поверхности, обусловленные снежным покровом Текст. / П.Д. Двулит, Б.Л. Скуин // Тезисы докладов. Симпозиум Комиссии Академии Планетарной Геофизики по изучению движения земной коры. Воронеж, 1988. - С. 95-96.

8. Дементьев, Ю.В. Подтверждение к алгоритму 1336 Текст. / Ю.В. Дементьев, Н.К. Шендрик // Библиотека алгоритмов: сб. 1516 2006. / под ред. М.И. Агеева-М.: Радио и связь, 1981,-С. 176.

9. Дементьев, Ю.В. Изменение ускорения силы тяжести под действием приливных деформаций Земли и атмосферных нагрузок Текст. / Ю.В. Дементьев // Материалы III международного конгресса «ГЕО-Сибирь-2007». Новосибирск: СГГА, 2007. - Т. 1, ч. 2. - С. 266 - 272.

10. Дементьев, Ю.В. Математическое моделирование движения точек земной поверхности Текст. /Ю.В. Дементьев // Изв. вузов. Горный журнал. -2007. -№ 5. С. 81 -85.

11. Дементьев, Ю.В. Учет сжатия Земли в топографических редукциях силы тяжести Текст. /Ю.В. Дементьев // Геодезия и картография. 2008. -№ 2. - С. 3 - 6.

12. Дементьев, Ю.В. Расчет топографических редукций силы тяжести по съемочным трапециям земного эллипсоида Текст. / Ю.В. Дементьев // Геодезия и картография. 2008. -№ 7. - С. 14 - 16.

13. Дементьев, Ю.В. Изменение ускорения силы тяжести под действием атмосферного давления Текст. /Ю.В. Дементьев // Сборник материалов V Международного научного конгресса «ГЕО-Сибирь-2009». Новосибирск, 2009. -Т.1, 4.1 - С. 80-81.

14. Дементьев, Ю.В. Численно-аналитический метод расчета топографической редукции в центральной и ближней зонах Текст. /Ю.В. Дементьев // Геодезия и картография. 2009. -№ 4. - С. 15-17.

15. Дементьев, Ю.В. Об уточнении формул для вычисления вертикальной составляющей величины притяжения масс сферического или плоского слоя Текст. / Ю.В. Дементьев, А.И. Каленицкий // Геодезия и картография. 2010. -№ 3. - С.7 - 9.

16. Дементьев, Ю.В. Связь плоских, сферических и сфероидических поправок за гравитационное влияние промежуточного слоя постояннойплотности и толщины Текст. / Ю.В. Дементьев // Геодезия и картография. 2010. -№ 4. С. 19 - 21.

17. Дементьев, Ю.В. Учет влияния топографических масс при вычислении возмущающего потенциала Текст. /Ю.В. Дементьев// Вестник СГГА. Новосибирск: СГГА, 2010. Вып. 1(12). - С.44 - 48.

18. Дементьев, Ю.В. О редукциях силы тяжести Текст. / Ю.В. Дементьев // Сборник материалов VI Международного научного конгресса «ГЕО-Сибирь-2010». Новосибирск, 2010. - С. 158 - 163.

19. Дементьев, Ю.В. Построение планетарной цифровой модели рельефа и её приложения Текст. /Ю.В. Дементьев, E.H. Кулик, Е.В. Дергачева // Сборник материалов VI Международного научного конгресса «ГЕО-Сибирь-2010».-Новосибирск, 2010.-Т. 1,ч.1.-С. 170- 173.

20. Дементьев, Ю.В. Зависимость поправок за плоский и сферический слои в неполной топографической редукции от толщины слоя и радиуса учитываемой зоны Текст. / Ю.В. Дементьев// Вестник СГГА. вып. 2(13). Новосибирск: СГГА, 2010. Вып. 2(13). - С. 13-17.

21. Дементьев, Ю.В. Рекомендации по выбору оптимальной модели рельефа для гравиметрии на территории Российской Федерации Текст. / Ю.В. Дементьев, А.И. Каленицкий, E.H. Кулик // Геодезия и картография.2011. -№.7.-С. 9-13.

22. Дементьев, Ю.В. О величине и значимости поправки за гравитационное влияние масс промежуточного слоя внешней области Текст. / Ю.В. Дементьев // Геодезия и картография. 2011. -№ 9. - С. 8-10.

23. Дементьев, Ю.В. Возможность интерполяции поправок, обусловленных влиянием масс промежуточного слоя внешней области при топографической редукции силы тяжести Текст. / Ю.В. Дементьев // Геодезия и картография. 2011. -№ 10. - С. 2-6.

24. Демьянов, Г.В. Построение общеземной системы нормальных высот Текст. / Г.В. Демьянов, А.Н. Майоров, М.И. Юркина // Геодезия и картография. 2009. -№.1. - С. 12 - 16.

25. Еремеев, В.Ф. Применение способа моделей к исследованию формул, определяющих форму Земли Текст. / В.Ф. Еремеев // Тр. ЦНИИГАиК. М.: Геодезиздат, 1950. - Вып. 71.-С. 11-77.

26. Еремеев, В.Ф. Теория ортометрических, динамических и нормальных высот Текст. /В.Ф. Еремеев // Тр. ЦНИИГАиК. М.: Геодезиздат, 1951. -Вып. 86.-С. 11-51.

27. Еремеев, В.Ф. Теория высот в гравитационном поле Земли Текст. / В.Ф. Еремеев, М.И. Юркина. М.: Недра, 1972. - 145 с.

28. Жарков, В.Н. Внутреннее строение Земли и планет Текст. / В.Н. Жарков. М.: Наука, 1983. - 416 с.

29. Игнатова, М.П. Оценка точности формул Нумерова для вертикального градиента силы тяжести на моделях Земли Текст. / М.П. Игнатова // Тр. ЦНИИГАиК М.: Геодезиздат, 1965. - Вып. 157. - С. 125134.

30. Инструкция по гравиразведке Текст. / отв. ред. К.Е. Веселов. М., 1980.-80 с.

31. Каленицкий, А.И. К проблеме повышения качества редуцирования гравитационного поля для геологической интерпретации Текст. / А.И. Каленицкий; -СНИИГГиМС. Новосибирск: 1987. - 28 с. - Деп. в ВИНИТИ, № 6404 - В87.

32. Каленицкий, А.И. Еще раз о редукциях в гравиметрии Текст. / А.И. Каленицкий// Вестник СГГА. Новосибирск: СГГА, 2006. - Вып.11. С. 98110.

33. Крылов, В.И. Международные системы отсчета: Современное состояние Текст. / В.И. Крылов, Буй Йен Тинь // Труды международной научно-технической конференции, посвященной 225-летию МИИГАиК. Геодезия. -М.: МИИГАиК, 2004. С. 187-195.

34. Машимов, М.М. Теоретическая геодезия Текст. /. М.М. Машимов -М.: Недра, 1991. -268 е.: ил.

35. Машимов, М.М. Высшая геодезия. Методы изучения Земли и создания общеземной системы геодезических координат Текст.: учебник / М.М. Машимов. -М.: ВИА, 1991. -552 с.

36. Методические рекомендации по учету влияния рельефа местности в гравиразведкеТекст.: практ. руководство МинГео СССР. Новосибирск: СНИИГГиМС, 1981. - 171 с.

37. Михайлов, A.A. Курс гравиметрии и теории фигуры Земли Текст.: учебник / A.A. Михайлов. М.: РЕДБЮРО ГУГК при СНК СССР, 1939. -432 с.

38. Молоденский, М.С. Методы изучения внешнего гравитационного поля и фигуры Земли Текст. /М.С. Молоденский, В.Ф. Еремеев, М.И. Юркина // Тр. ЦНИИГАиК. М.: Геодезиздат, 1960. - Вып. 131. - 264 с.

39. Мудрецова, Е.А. Гравиразведка: Справочник геофизика Текст. / Е.А. Мудрецова, К.Е. Веселов и др. -М.: Недра, 1990. 607 е.: ил.

40. Нумеров, Б.В. Зависимость между местными аномалиями силы тяжести и производными от потенциала Текст. / Б.В. Нумеров // ДАН СССР, 1929. №24,-С. 101-105.

41. Пеллинен, Л.П. Исследование уклонений отвеса и вывод фигуры квазигеоида на Кавказе Текст. / Л.П. Пеллинен // Тр. ЦНИИГАиК. М.: Геодезиздат, 1951. Вып. 86. - С. 62-96.

42. Пеллинен, Л. П. Влияние топографических масс на вывод характеристик гравитационного поля Земли Текст. / Л.П. Пеллинен // Тр. ЦНИИГАиК. М.: Геодезиздат, 1962. - Вып. 145. - С. 23-42.

43. Пеллинен, Л.П. Об учете влияния топографических масс при вычислении уклонений отвеса и высот квазигеоида Текст. / Л.П. Пеллинен, О.М. Остач // ЭШсНа geophysika & geodaetica, 1974. № 4. - С. 319-328.

44. Смирнов, В.П. Формула для вычисления поправки за гравитационное влияние дневного рельефа в центральной зоне Текст. / В.П. Смирнов // Материалы по геологии, геофизике и полезным ископаемым Сибири. -Новосибирск, 1974. С. 62-64.

45. Шимбирев, Б.П. Теория фигуры Земли Текст.: учебник / Б.П. Шимбирев. М.: Недра, 1975. - 432 с.

46. Эскобал, П. Методы астродинамики Текст. / П. Эскобал. М.: Мир, 1971.-342 с.

47. Юркина, М.И. Вычисление первой и второй вертикальных производных силы тяжести по картам её аномалий Текст. /М.И. Юркина// Тр. ЦНИИГАиК. -М.: Геодезиздат, 1965. Вып. 157. - С. 116-124.

48. Юркина, М.И. Оценка на модели Земли отдельных членов поправки к стоксову приближению коэффициентов разложения возмущающегопотенциала Текст. / М.И. Юркина, А.Б. Старостина// Тр. ЦДИИГАиК. М.: Геодезиздат, 1965.-Вып. 157-С. 58-65.

49. Юркина, М.И. К вопросу об учете влияния центральной зоны на величину поправок Gx Молоденского Текст. / М.И. Юркина, JI.E. Кропотова// Тр. ЦНИИГАиК. М.: Геодезиздат, 1969. - Вып. 176. - С. 89-97.

50. IERS Technical Note 32. IERS Conventions (2003) Electronic resource. / D.D.Mc/ Carthy and G. Petit (eds.) Англ. - Режим доступа: ftp://maia.usno.navy.mil/conv2000/. - Загл. с экрана.

Информация о работе

Дементьев, Юрий Викторович
доктора технических наук
Новосибирск, 2012
ВАК 25.00.32

Диссертация

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

Автореферат

Похожие работы