CC BY
46
УДК 528:[551.2/.3+550.38+55]
А.И. Каленицкий, А.П. Карпик, В.С. Середович СГГ А, Новосибирск
НЕОБХОДИМОСТЬ ИНТЕГРАЦИИ НАУК В РЕШЕНИИ ЗАДАЧ ГЕОДИНАМИКИ ДЛЯ ОБЕСПЕЧЕНИЯ УСТОЙЧИВОГО РАЗВИТИЯ ТЕРРИТОРИЙ
Человеческое общество обеспокоено происходящими на Земле природными катаклизмами (землетрясениями, извержением вулканов, цунами, ураганами). Они наносят огромный моральный и материальный урон, приводят к гибели людей, разрушениям поселений и различных сооружений. Все в большей степени на стихийные бедствия налагаются негативные проявления техногенного характера, вызванные внедрением интенсивных технологий освоения и переработки природных ресурсов Земли, часто без объективной оценки последствий для живой природы, растительности и населения. Уже становится ясным, что антропогенные факторы становятся соизмеримыми с природными, а иногда и превышают их по вредному воздействию на систему: человек - природа.
Все острее ощущается необходимость в определенной степени оперативном, но объективно-эффективном комплексном изучении и анализе ситуации с целью принятия соответствующих упреждающих решений на различных уровнях управления, хозяйствования и ориентации населения по предотвращению последствий негативных процессов, нарушающих или осложняющих экологическую обстановку, устойчивость развития территорий и сложных инженерно-технических систем.
Речь, таким образом, идет о необходимости формирования новой модели развития цивилизации, гарантирующей устойчивость развития территорий, природно-технических систем на основе регулируемого управления процессами, происходящими в ландшафтной оболочке Земли, объединяющей литосферу, атмосферу, гидросферу и биосферу [1, 2, 3].
Наблюдения за состоянием ландшафтной оболочки Земли осуществляются в ходе проведения мониторинга окружающей среды. При этом изучением механизма проявления сейсмотектонических процессов, их отражения в изменении пространственного положения поверхности Земли и объектов на ней, а также (гео) физических полей призвана заниматься геодинамика.
Известно, что геодинамика возникла как наука междисциплинарная. Ее развитие, по крайней мере, в СССР, а затем - в России осуществлялось и финансировалось в основном, по отраслевому принципу. Координация усилий, обобщение результатов исследований в определенной степени проводились в институтах АН СССР (в настоящее время - РАН). Но при этом ведомственный принцип, по сути, оставался превалирующим. Геология и геофизика, геодезия, геомеханика, геоморфология, геоэкология и другие науки о Земле интерпретировали и интерпретируют получаемые результаты со своих позиций с “вялым” привлечением результатов и выводов смежников. Это в определенной степени сказывается на снижении эффективности прогнозов проявления и оценки последствий стихийных и техногенных
воздействий. Вместе с тем, комплексирование методов исследований всегда эффективнее, чем применение одного, признаваемого ведущим. Например, сейсмология надежно регистрирует положение очага напряженной деформации в земной коре, когда в нем начинается “треск”, проявляемый в виде сейсмических волн. Но это происходит относительно незадолго до разгрузки напряжения, выражаемой разломом слоев земной коры, или их сдвигом. Задолго до этого начинает происходить перемещение огромных масс в подкоровом пространстве, нарушается изостатическое равновесие земной коры. Все это должно находить отражение в различных геофизических полях, в том числе в гравитационном (возможно - в магнитном, тепловом, радиационном) и может повлиять на изменение положения земной поверхности. Комплексное применение геолого-геофезических, геодезических и других методов, совместная обработка и интерпретация результатов может и должна, в конечном результате, обеспечить выявление признаков приближения гидравлического удара, извержения вулкана, деформации земной поверхности и других угрожающих явлений. Очевидно, что при этом должна быть обеспечена общая координация работ.
Геодинамика как научный комплекс претерпевает в последнее время существенные изменения в методологическом обосновании и технологическом обеспечении в связи с интенсивным внедрением спутниковых методов пространственно-координатных определений, дистанционного зондирования Земли, ГИС-технологий, 3D-компьютерного моделирования среды изучения по результатам пространственного координирования поверхности и геофизического зондирования земных недр и так далее. Однако при этом возникает ряд проблем, особенно в части, касающейся интерпретации результатов с привлечением полученных ранее производных материалов и данных, когда первичная исходная информация не может быть восстановлена в связи с утилизацией или из-за применения устаревших технологий обработки с потерей точности на выходе. По нашему мнению, накопился огромный объем устаревшей информации, которая не соответствует современным требованиям как по точности первичных данных, так и, особенно, по производным результатом. Эта информация содержится в различных базах данных и знаний. Рано или поздно ее придется или переобрабатывать или утилизировать.
Особого внимания заслуживает вопрос о необходимости пересмотра процедуры подготовки содержания и порядка доведения до исполнителей нормативных документов. Они должны оперативно обновляться и содержать в основном ключевые положения конечных результатов, давая исполнителям в тесном взаимодействии с наукой возможность творческой реализации программ исследований.
Рассмотрим в связи с этим некоторые примеры, касающиеся применения гравиметрии как одного из основных методов в физической, динамической геодезии и в прикладной геофизике применительно к решению задач глобальной, региональной и локальной геодинамики.
Геодезия как научный и производственный комплекс призвана обеспечивать все виды творческой и хозяйственной деятельности человеческого общества надежной пространственно-координатной основой, соответствующей по точности и полноте содержания запросам любых потребителей. Кроме того, основной научной задачей геодезии изначально являлась и является определение (с последующим уточнением) параметров фигуры Земли и ее гравитационного поля. При этом важность ее решения постоянно возрастает в связи с необходимостью оценки степени стабильности взаимного пространственного положения материковых плит (глобальная проблема), массивов суши и островов в океане (глобальнорегиональная проблема), поверхности платформенных и складчатых образований, равнин, отдельных сейсмоактивных территорий (региональная проблема), районов и участков, интенсивно осваиваемых в связи с добычей и переработкой природных ресурсов, строительством крупных промышленно -энергетических комплексов и инженерных коммуникаций (локальная проблема).
Повсеместное внедрение спутниковых методов координатных определений, введение новой референцной системы координат СК-95 открыло широкие возможности для совершенствования геодезических технологий, повышения точности определения пространственного положения поверхности Земли и различных объектов на ней применительно к решению любых задач глобального, регионального и локального характера [4, 5, 6].
Спутниковые координатные определения позволяют получать теперь пространственное положение любой определяемой точки (М) непосредственно относительно центра масс Земли и в геодезической системе координат (В, Ь, НМ). НМ - геодезическая высота определяемой точки М над принятым за нормальную Землю эллипсоидом (референцным или общим), отсчитываемая по нормали к нему [7].
Проблема состоит в том, что большинство потребителей не испытывает потребности в геодезических высотах. Необходимо, строго говоря, обеспечивать их ортометрическими высотами местности, отсчитываемыми от поверхности геоида (поверхности гидростатического равновесия) по силовым линиям силы тяжести (отвесным линиям),
где g - текущее значение силы тяжести вдоль силовой линии, начиная от отсчетного геоида (точка О) до определяемой точки М;
gm - среднее интегральное значение силы тяжести на отрезке ОМ.
До настоящего времени из-за невозможности определения положения геоида внутри топографических масс, слагающих рельеф местности, по результатам высокоточного нивелирования определяются и используются в основном нормальные высоты, отсчитываемые от близкой к геоиду поверхности квазигеоида. При этом в (1) величина gm заменяется величиной ут - среднего значения нормальной силы тяжести на отрезке О'М, где О'-
(1)
положение точки на эллипсоиде (референцном или общем), образуемой пересечением с ним нормали, проходящей через определяемую точку М.
Разница между геодезической и нормальной высотой определяет аномалию высот квазигеоида над эллипсоидом относимости нормальной силы тяжести (у0)
См = Нм - нМ. (2)
Квазигеоид, в отличии от геоида, не является, строго говоря, уровенной поверхностью, характеризуемой возмущающим гравитационным потенциалом как интегральной (вторичной) функцией силы тяжести. Его отличие (А£) от геоида оценивается в идеальном случае для равнинных районов в несколько сантиметров, в горных 2 - 3 м.. Фактически эти величины на практике могут отличаться в большую сторону в несколько раз.
До последнего времени аномалия высоты определялась через возмущающий потенциал
См = тмум , (3)
где ум - значение нормальной силы тяжести в какой-то точке М,
отличающейся по высоте от точки М на величину АС;
ТМ - разность потенциала силы тяжести в точке М и пересчитанного в нее (с ошибкой в высоте, равной АС) потенциала нормальной Земли в виде или референц-эллипсоида, или общеземного.
Заметим, что для решения самых различных хозяйственных задач нормальные высоты отвечают предъявленным требованиям, так как изменение величины АС на ограниченных территориях незначительно. Все топографические карты отражают рельеф земной поверхности в нормальных высотах.
При мониторинге геодинамических процессов, особенно в глобальном, а во многих случаях - и в региональном плане сейчас, несомненно, выгоднее использовать геодезические высоты, определяемые не по гравиметрическим данным (что весьма сложно и не всегда корректно), а непосредственно спутниковым методом. Сложность и некорректность связаны с тем, что значения ТМ определяются интегрально вокруг точки М по смешанным аномалиям силы тяжести, когда используются нормальные высоты (нужно было бы геодезические), а пересчет в них нормального гравитационного поля должен был бы производиться с поверхности эллипсоида относимости, а не с квазигеоида. При этом интегрирование аномалий, как правило, осуществляется уже вообще без учета высот результативных точек и пунктов, в которых производилось измерение силы тяжести. В конечном итоге возникает очевидная неопределенность: к какому же уровню относить значение Т?
Положение усугубляется еще и тем, что при расчетах использовались и используются аномалии силы тяжести в различных редукциях. При этом априори полагается что, например, аномалии силы тяжести в редукции Буге по точности соответствуют сечению изоаномал отчетных и изданных карт. Однако такое мнение во многих случаях не соответствует реальному
положению. Это неоднократно доказывалось и показывалось [8, 9, 10]. Внедрение новой высокопроизводительной технологии редуцирования [8] задержалось в силу ряда причин объективного и субъективного характера. В частности, завершалось издание гравиметрической карты страны в масштабе 1:200000. Номенклатурные листы карты, изданные ранее, не “сводились” бы с новыми, причем существенно; требовала пересмотра нормативно-техническая база гравиметрических работ. В период перестройки вопросы гравиметрии на определенный срок отошли на второй план. Далее сказалась инерция: нормативная база осталась “старой”, учебные программы в ВУЗах и техникумах геологоразведочного и геодезического профиля - тоже. Однако реалии современной действительности все настойчивее требуют приведения результатов гравиметрии в соответствие с современными требованиями. Остановимся, в связи с этим, на основных, принципиальных позициях, которые необходимо учитывать при обработке и интерпретации результатов гравиметрии, как в геодезии, так и в приложении к разведочной геофизике и геодинамике.
1. Поправки за гравитационное влияние рельефа местности которые учитываются практически во любых редукциях исходного гравитационного поля, изменяются как в плане, так и по высоте нелинейно, а значит - не могут интерполироваться. В операторном выражении их величина может быть представлена в виде
Г. Гк Г. Гк )+А&рф, (4)
где гн, гк - соответственно внутренний и внешний радиусы внешней учитываемой области;
дР - оператор вычисления поправок за рельеф в плоском варианте при
исправлении плоско-параллельного слоя Буге на реальный объем топографических масс;
Ад^р? - оператор приведения величины “плоской” поправки к
сферическому виду.
В свою очередь, значения дР и Адрф выражаются в виде
дРГ,Гк,Нц) = + н2 - 2Нц ■ Но| _
Адр = е{нц - Н) ]
где а и F - соответственно постоянные коэффициенты, характеризующие геометрию учитываемой области в плоском и сферическом представлении;
Нц - высота точки вычисления поправок (результативной точки),
Н2 и Н, Н0 - соответственно среднее весовое квадратическое значение и средние весовые значения высот рельефа местности вокруг результативной точки в пределах учитываемой области;
Ну - отметка рельефа в центре элементарной площадки, на которые
разделяется учитываемая область в радиальном и угловом направлениях (і = 1, 2, ..., пі,у = 1, 2, ..., ті). =
В качестве весовых функций для определения Н0 и Н используются соотношения плановых размеров зон, на которые разделяется учитываемая область, а для определения Н - соотношения коэффициента F для каждой зоны.
Оператор для поправки за влияние промежуточного слоя топографических масс имеет вид
4гсф пр сл(г,,Гк,Щ)= оГ2НцНо -Н21 + ^• Н. (6)
Вместе с тем выражение (5) для плоской поправки за рельеф можно представить в виде двух слагаемых
'=2 '
!Т(г«’ГкН)=%о +%' = о( Н - Но2| + а(Нц - Но)2. (7)
Первое слагаемое не зависит от величины Нц, но зависит от степени
изрезанности окружающего рельефа, второе меняется по параболическому закону, равняясь нулю при Нц = Н0 и возрастает по мере увеличения
абсолютной разности Нц - Н0 . Другими словами поверхность Н0 - эта
поверхность минимальных (но не нулевых, если рельеф - не плоскость) значений поправок.
На основе приведенных выражений можно сделать ряд важных выводов.
- Сферическая составляющая в (6) общей поправки за влияния топографических масс промежуточного слоя не зависит от высоты результативной точки. Вычисленные один раз значения Н однозначно позволяют для любой территории перейти от редукции Буге к неполной топографической, особенно важной для геодезии, и, наоборот - от неполной топографической редукции к редукции Буге, используемой в прикладной гравиметрии.
- Имеется возможность определения значений Н0 (гн, Гк ), Н(гн, Гк ), и $80 (гн, гк) на конкретную территорию всего один раз, если рельеф на ней не претерпел существенных изменений, для гравиметровых работ самого разного масштаба. Величины Н0 и 0 интегральны по сути и могут интерполироваться в плане, то есть изображаться в виде планов и карт в электронном и, при необходимости, в графическом представлении. Значения <8' могут быть также заранее представлены в виде электронной таблицы или аналитически. В результате, по мере получения пространственных координат точек гравиметровых измерений, для них сразу же можно определять поправки за рельеф во всей внешней учитывая области согласно (7).
Приведенные операторные соотношения для круговой учитываемой области полностью соответствуют квадратной, когда цифровая модель рельефа
(ЦМР) определяется массивом отметок высот в узлах квадратной сетки или тригональной.
- Отпадает как ненужная, и даже вредная трудоемкая операция осреднения высот рельефа в элементарных ячейках палеточного поля учитываемой области. Другими словами, генерализация снижает величину до и в целом всей поправки на значение, переменное в каждой
результативной точке. Вместе с тем дискретность шага (Аг или Ах, Ау) задания ЦМР также обусловливает генерализацию и искажение поправок за рельеф и, как следствие - за промежуточный слой топографических масс. В связи с этим были выявлены информативные характеристики изрезанности рельефа, разработан математический аппарат автоматизированного определения оптимальных условий задания ЦМР и редуцирования гравитационного поля с условием определения поправок за рельеф, влияние промежуточного слоя и аномалий силы тяжести с необходимой точностью [8].
2. В гравитационных аномалиях Буге, представляемых на изданных картах различного масштаба, как правило, отмечается сложная картина корреляции поля силы тяжести с высотами рельефа местности. Детальный анализ показал [9], что физико-геологических причин наблюдаемой картины может быть, в основном, две: 1 - унаследовательность земной поверхностью формы границ раздела геоструктур с различной плотностью (складок, тектонических блоков, изостатической поверхности основания земной коры и других), 2 - глобальные плотностные неоднородности за пределами территории гравиметровых измерений (водные массы крупных озер, водохранилищ, морей и океанов, горных блоковых массивов), отличающихся по плотности от горных пород в районе работ. В отличие от бытующего еще мнения такие массы обусловливают не только региональный гравитационный фон, но и сугубо локальную составляющую, прямо зависящую от высоты точек гравимеировых измерений. Следует отметить, что учет и исключение гравитационного влияния таких масс должен тщательно выполняться и учитываться только после редуцирования измеряемого гравитационного поля - на этапе интерпретации аномалий силы тяжести, но до выделения локальных аномалий. В связи с этим особо подчеркнем, что в наблюденные значения силы тяжести на поверхности моря нельзя вводить никакие редукции поправки за влияния водных масс (или рельефа морского дна), за исключением поправок, связанных с волнением водной поверхности. Измеренные значения силы тяжести здесь после прямого вычитания из них величины у(#0)являются одновременно и аномалиями в полной топографической редукции, и в свободном воздухе, так как морская поверхность - уровенная, то есть геоид.
К искажающим факторам, обусловленным несовершенством методики редуцирования измеряемого гравитационного поля, относятся: 1 -
использование переменной величины плотности промежуточного слоя топографических масс или поправок за рельеф, тогда как она должна
считаться постоянной (а0 = const); 2 - учет сферичности Земли при определении поправок за рельеф для “плоской” редукции Буге; 3 - различие размера учитываемой области в редукции Буге, когда плоско-параллельный слой полагается бесконечным, а учет поправок за рельеф осуществляется в области, ограниченной радиусом rK, - создается искусственно ложный эффект “боковых” масс; 4 - необоснованное произвольное или директивное (нормативное) установление величины rK; 5 - замена геодезических высот точек гравиметровых измерений на нормальные.
Без исключения влияния первых четырех факторов в процессе редуцирования гравитационного поля (выделения аномалий) не может быть и речи о точной гравиметрии при разведочных работ и мониторинге геодинамических процессов на локальных территориях. При проведении геодинамических исследований в региональном и, особенно - в глобальном масштабе, надо стремиться исключать влияние пятого фактора.
3. Важнейшее теоретическое и методологическое значение на современном этапе имеет необходимость учета разновысотности точек гравиметровых измерений и соответствующих им аномальных значений силы тяжести в процессе математического и геоинформационного моделирования объектов исследований в аномальном гравитационном поле. Должно учитываться и время измерений силы тяжести. Пространственная фиксация аномального поля и его вариаций во времени не только не осложняют, а наоборот, обеспечивают повышение качества прогнозных оценок в геодинамике и при разведке недр, изучения фигуры Земли и формы ее поверхности.
4. Принципиальное методологическое, методическое и практическое значение для геодезической и прикладной гравиметрии имеет решение двух вопросов.
Первый относится к процедурам: а - редуцирования исходного гравитационного поля при его разделении на нормальную и аномальную составляющую, б - редуцирования гравитационных аномалий (разделяя их на региональную и локальные составляющие). Требования к выполнению процедуры “а” сформулированы ранее в п.2. Подчеркнем только еще раз, что на этапе выполнения этой процедуры плотность промежуточного слоя (с0) должна полагаться постоянной и равной даже при раздельном учете влияния параллельного слоя и поправки за рельеф.
Второй связан с определением последовательности учета и исключения гравитационного влияния известных аномалеобразующих факторов - при геологическом редуцировании (термин впервые введен Э.Э. Фотиади). При этом сразу же следует подчеркнуть, как также отмечалось в п.2, что геологическое редуцирование, включая учет и исключение гравитационного эффекта геоструктур, в том числе “боковых” масс, должно выполнятся после процедуры “а”, но до процедуры “б”. Это позволяет получить остаточное поле гравитационных аномалий, которое при его разделении (редуцировании) на региональную и локальную составляющие обеспечивает количественную (а не качественную, как при разделении “полного”
аномального поля силы тяжести) сопоставимость локальных аномалий с модельными гравитационными эффектом объекта исследований. Степень их соответствия является критерием качества моделирования исследуемого объекта [10].
Требование к постоянству величины оь не относится к случаю учета поправок (А£ и А^) в уклонения отвеса за влияние топографических масс в центральной зоне (до радиуса 48.5 км) [7]. Более того, она должна быть переменной, так как на изменение положения отвеса влияют не только превышения рельефа относительно уровня результативной точки, но и дифференциация топографических масс по плотности, что и отражает аномальное гравитационное поле во внешней учитываемой области (от 48.5 до 1000 км), используемое при определении возмущающего потенциала (Т), уклонений отвеса и аномалий высот в геодезии. И это соответствует отмеченным в п.2 требованиям, так как процедура определения значений Т является процессом моделирования гравитационного поля. Следует в связи с этим отметить, что в геодезии, к сожалению, часто выполняется математическое моделирование оценочное (приближенное), как предварительное в прикладной гравиметрии. Поэтому на территориях сопряжения осадочных бассейнов с образованиями платформенного или складчато-орогенного типа ошибки вычисленных поправок в уклонения отвеса и высоту квазигеоида над эллипсоидом возрастают, искажая конечные результаты (в триангуляции, нивелировании).
Вместе с тем предварительное моделирование допустимо в случае оценки требований к точности топографо-геодезического обеспечения высокоточной гравиметрии при изучении геодинамических процессов, особенно в локальном плане, когда ведется мониторинг техногенной геодинамики. В этом случае нормативные требования, основанные на точности определения нормальной силы тяжести и ее вертикального и горизонтального градиентов, становятся неприемлемыми. Требуется оценка по аномальным градиентам [11].
5. Вышеизложенное показывает, что методология и технология обработки и интерпретации результатов гравиметрии требуют коренной корректировки. По большому счету, карты гравитационных аномалий для территорий с рельефом средней и высокой степенью изрезанности, перепадом высот от 200 м и более должны быть исправлены с исключением проявления методических факторов, указанных в п.2.
По-видимому настало время пересмотреть систему определения нормальных высот на основе переориентации некорректной в некоторой степени технологии определения положения уровня квазигеоида внутри топографических масс на однозначную и корректную концепцию объемного представления поля ортометрических (нормальных) и геодезических высот во внешнем относительно земной поверхности пространстве. Вычисление высот ниже этой поверхности не требуется. Предпосылки для перехода к такой концепции имеются: 1 - прямое определение геодезических высот точек местности спутниковым методом, 2 - совмещение точек спутниковых
определений с реперами нивелирования, 3 - возможность совмещения нивелирных ходов с гравиметрическими и спутниковыми, что при имеющихся сейчас высокопроизводительных, портативных гравиметрах несложно, в том числе по затратам, а математическое обеспечение возможности определения ортометрических (нормальных) высот достигается (при оптимизации дискретности шага измерений) численным интегрированием согласно выражению (1); 4 - возможность определения возмущающего потенциала на земной поверхности, его аппроксимации и пересчета вверх. Аппроксимация может производиться методом подбора подсобной модели, например, точечных масс или тонкого слоя. Гармоническая аппроксимация на неровной поверхности может дать загрубление результатов.
Базовыми опорными площадками (площадными реперами) для дальнейших построений полей высот могут служить геодинамические, в том числе локальные техногенные полигоны, где регулярно проводится комплекс геодезическо-гравиметрических измерений и построений картосхем локальных (квази) геоидов и в пределах которых затем может осуществляться спутниковое нивелирование для нужд потребителей. Вместе с комплексом ФАГС и ВГС они, в конечном итоге, позволят обеспечить решение задач глобалиного , регионального и локального значения.
В заключение отметим, что один из основных этапов сложного процесса геодинамических исследований является интерпретация полученных результатов. Ключевым процессом при этом является математическое объемно-геоинформационное моделирование. Важнейшее место при этом отводится сбору всех имеющихся данных о физических и механических свойствах горных пород, слагающих геологический разрез района работ, а также данные о его геологическом строении и тектонике. Вся эта информация, будучи систематизированной и обработанной с выявлением различных особенностей изменения физических и механических свойств в плане и по глубине, их связи с отражением в физических полях, в том числе гравитационном, используется при пространственных модельных построениях. На основе установленных связей выполняется комплексная интерпретация результатов геодезических и геофизических определений. По мере проведения повторных циклов геодинамических исследований, позволяющих перейти к построению пространственно-временных вариаций, в том числе земной поверхности и гравитационных аномалий, выявляются на основе комплексирования с данными других методов закономерности развития геодинамики, в том числе осложнения в сейсмотектонической обстановке. При региональных комплексных исследованиях устанавливается идентификация движений и напряженно-деформированного состояния геодинамических систем [12].
Таким образом, практика работ показывает, что комплексирование геодезических и геолого-физических исследований геодинамических процессов, комплексная интерпретация результатов обеспечивает повышение качества прогнозных построений, способствуя своевременному принятию
мер по предотвращению негативных последствий воздействия природных и техногенных факторов. Это подтверждает необходимость координации и интеграции наук о Земле при решении задач геодинамики с целью обеспечения устойчивого развития территорий.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Конференция ООН по окружающей среде и развитию (Рио-де Жанейро, июнь-1992) [Текст]: Иформ.обз. / В.А. Коптюг. - Новосибирск: СО РАН, 1992. - 69 с.
2. Устойчивое развитие: ресурсы России [Текст]/Под.общ.ред. акад. РАН Н.П. Лаверова. - М.: Издат. центр РХТУ им. Д.И. Менделеева, 2004. - 212 с.
3. Карпик А.П.Метологические и технологические основы геоинформационного обеспечения территорий [Текст]: Монография. - Новосибирск: СГГА, 2004. - 250 с.
4. Бородко, А.В. Развитие системы геодезического обеспечения в современных условиях [Текст] / А.В. Бородко, Н.Л. Макаренко, Г.В. Демьянов // Геодезия и картография. - 2003. - № 10. - с. 7 - 13.
5. Хабаров В.Ф. О развитии ФАГС и ВГС и постоянно действующих спутниковых станций [Текст] / В.Ф. Хабаров // Геодезия и картография. - 2004. - № 10. - С. 15 - 19.
6. Единая государственная система геодезических координат 1995 года (СК-95) [Текст] / Н.Л. Макаренко и др.; под ред. А.А. Дражнюка. - М.: Федеральная служба геодезии и картографии России, 2000. - 34 с.
7. Шимбирев Б.П., Теория фигуры Земли [Текст]. - М.: Недра, 1975. - 432 с.
8. Методические рекомендации по учету влияния рельефа местности в гравиразведке [Текст]: Практическое руководство МинГео СССР / А.И. Каленицкий,
В.П.Смирнов; Отв. за выпуск А.И. Каленицкий. - Новосибирск: СНИИГГиМС, 1981. - 171 с.
9. Каленицкий А.И. К проблеме повышения качества редуцирования гравитационного поля для геологической интерпретации [Текст] / А.И. Каленицкий. - Новосибирск: СНИИГГиМС, 1987. - 28 с. - Деп. в ВИНИТИ, № 6404-В.87.
10. Каленицкий А.И. Геодезическо-гравиметрический мониторинг техногенной геодинамики инженерных сооружений [Текст] / А.И. Каленицкий // Геодезия и картография. - № 8. - 2000. - С. 24 - 27.
11. Каленицкий А.И. О требованиях к топографо-геодезическому обеспечению инженерно-гравиметрических исследований [Текст]/А.И.Каленицкий//Геодезия и картография.-№ 7.-1995.-с.12-14.
12. Середович В.А. Идентификация движений и напряженно-деформированного состояния самоорганизующихся геодинамических систем по комплексным геодезическим и геофизическим наблюдениям [Текст]: Монография / В.А. Середович, В.К. Панкрушин, Ю.И. Кузнецов. Б.Т. Мазуров, В.Ф. Ловягин; под общей редакцией В.К. Панкрушина.-Новосибирск: СГГА, 2004. - 356 с.
© А.И. Каленицкий, А.П. Карпик, В.С. Середович, 2005
