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

в виду, что значения функций в точках, расположенных по пери­ метру всей (шириной я » 2s-f-3s) результативной карты, вычисляются с меньшей точностью, чем в других частях поля. Это обусловлено тем, что при вычислении значений в этих пунктах точки счета на­ ходятся на краю окрестности выборки.

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

Общее время счета для теста составляет 20 мин.

164

При работе ОсП7 используется СП с номерами: 0000; 0010; 0012; 0027; 0042; 0154.

1.ПОДГОТОВКА ИСХОДНОЙ ИНФОРМАЦИИ К СЧЕТУ

Параметрами счета к ОсП7 являются: шаг задания исходной функции s, высоты пересчета z, сечение трансформированных карт С. Шаг задания исходной функции s выбирают равным 1 см в масштабе исходной карты. Оптимальные в смысле точности высоты пересчета и сечение трансформированных карт выбирают по табл. 13.

1°. Первый массив (Ml) информации к задаче трансформации содержит выбранные параметры счета. Он состоит из трех перфо­ карт.

Первая перфокарта массива M l содержит следующие двенадцать

десятичных чисел:

1-ая строка — пг — число уровней пересчета исходного поля в верхнее полупространство.

165

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

заданных в узлах квадратной сети. Этот массив готовится вручную, если карта задана в изолиниях, или с помощью ЭВМ, если задан каталог гравиметрических пунктов (см. гл. X V , описание ОсП.5).

166

Рассмотрим ручную подготовку массива М2. На исходной карте разбивается квадратная сеть точек с шагом s, равным 1 см в масштабе карты. Это осуществляется либо копированием исходной карты на прозрачную миллиметровку, либо нанесением на исходную карту квадратной сети точек. Карта интерполируется и выписываются значения Ѵг в узлах равномерной сети. Точность интерполяции является функцией от точности и масштаба исходной карты. В за­ висимости от изменения градиента Ѵхг и Ѵуг точность интерполя­ ции будет переменной по площади. Если горизонтальные градиенты невелики (0,2—0,4 мгл/км), то при вычислении Ѵг в узлах равномер­ ной сети сечение карты для облегчения интерполяции искусственно уменьшается рисовкой дополнительных изоаномал, форма которых повторяет форму двух соседних изолиний. Если на исходной карте имеются крупные замкнутые аномалии, то проводятся профили вкрест их простираний, строятся графики по этим профилям в про­ извольно выбранном вертикальном масштабе и оценивается зна­ чение V. в их центре по максимальному значению на графике. Ин­ терполяция внутри замкнутой аномалии проводится с учетом этого значения. Проинтерполированная карта разбивается на квадраты размером 16sX'16s, где s — равно 1 см в масштабе карты. Квадратам присваиваются порядковые номера, начиная с левого верхнего ква­ драта, от нуля до п — 1 , где іі — общее число квадратов (см. рис. 32). Значения Ѵг в узлах квадратной сети переписываются на бланки

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

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

2.Аналогично выписываются значения Ѵг в следующих по по­ рядку квадратах.

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

Для каждого квадрата получается своя контрольная сумма, которая пробивается на перфокарте, помещаемой в конце массива. Таким образом, массив М2 будет состоять из перфокарт, содержа­ щих значения исходной функции по квадратам. Каждый квадрат содержит максимум 22 перфокарты (если запись значений по одному числу в ячейку); затем следует перфокарта с KS квадрата, а после

167