Файл: Автоматизированная система обработки и интерпретации результатов гравиметрических измерений..pdf

и с 484 до 124 As. С увеличением размеров s эффективность алгоритма возрастает [65].

Г Л А В А V

ЗАДАЧИ ИНТЕРПОЛЯЦИИ

В гл. I было отмечено, что основным результатом процесса обра­ ботки является функция U (п, Сп), т. е. значения И З О Л И Н И Й И их координаты. Функция U (п, Сп) изображается в виде карты изо­ линий аномальных значений силы тяжести в определенной редукции.

Можно представить два пути получения U (п, Сп). Первый состоит в том, чтобы работали два оператора:

ная во всех последующих задачах: различного вида трансформа­ циях функции U (s, s) и решении обратных задач. В этом случае существенно упрощаются алгоритм u организация программ решения этих задач.

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

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

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

1.МЕТОД ПАРАБОЛИЧЕСКОЙ ИНТЕРПОЛЯЦИИ ФУНКЦИИ

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

46

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

Будем считать, что на небольшой области изменения Іх, у] функ­ ция U (х, у) может быть достаточно хорошо приближена параболой

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

S Pk[Uk(x, y)-F£>(x, у)]* = тш; (Ѵ.З)

47

ветственно их горизонтальными производными. В (Ѵ.З) суммиро­ вание ведется по всем точкам Ak, лежащим в окрестности радиуса R вокруг точки В (см. рис. 6).

Дифференцируя выражения (Ѵ.З) по неизвестными, Ь, с, d, е, / и приравнивая нулю частные производные, получим систему уравне­

Решая систему (Ѵ.5) одним из известных методов, например методом исключения Гаусса по стандартной подпрограмме, найдем искомые коэффициенты. Задав шаг s между точками счета и каждый раз помещая начало координат в точку счета (центр окрестности радиуса R), получим систему элементарных перекрывающихся параболоидов, аппроксимирующих значения U (х, у) во всей области

еезадания.

Ввыражении (Ѵ.З) через Рк обозначена некоторая весовая функция. Как известно, за вес можно принять функцию, обратно пропорциональную квадрату среднеквадратических погрешностей каждой точки [109]. При гравиметрических наблюдениях значения

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

48

авторов показали, что более предпочтительна введенная Лапортом [119] весовая функция вида

где R — радиус окрестности элементарного параболоида; dk— рас­ стояние от точки В (0,0) до точек Ak (xk, yk); л , ѵ — некоторые коэффициенты.

2. АЛГОРИТМ ОБРАБОТКИ ПЛОЩАДНЫХ ГЕОФИЗИЧЕСКИХ НАБЛЮДЕНИЙ

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

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

порядка 2R должна выбираться в соответствии с текущими коорди­ натами рассчитываемой точки.

4. Изменять размеры элементарной матрицы порядка 2R в зави­

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