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

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 197

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

1

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

В дополнение к рис. 25, на котором изображена структура си­ стемы обработки и хранения данных, на рис. 29 показаны некоторые связи и элементы этой структуры, т. е. нечто вроде технологического процесса обработки на базе ЭВМ. При его построении автор исполь­

зовал не только

собственный опыт и практику применения ЭВМ

в гравиразведке,

но и опыт работы с АСО в некоторых экспедициях

{Н. И. Бакланов) и в трестах Киевгеология и Укргеофизразведка <Г. С. Бутаков, В. И. Старостенко, В. И. Дядюра).

Ч А С Т Ь В Т О Р А Я

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

Программы, составляющие автоматизированную систему обра­ ботки, делятся на основные (ОсП) и вспомогательные (ВсП). Основ­ ная программа решает одну из основных гравиметрических задач. Каждая ОсП имеет ряд модификаций, решающих на основе того же численного метода более узкие задачи. Функциями ВсП являются формирование и организация числовых данных в виде, необходи­ мом для ОсП, а также настройка ОсП в соответствии с информацией. По основным и вспомогательным программам подвергаются обра­ ботке числовые массивы ограниченных размеров (МИД массивы исходных данных).

Помимо числовых данных, каждой из программ необходимо за­ давать некоторую информацию, характеризующую массивы исход­ ных и результативных данных, и режим работы. Эта информация задается в карте информации о счете и в карте информации о маг­ нитной ленте. В первой из них всегда указываются параметры счета Ii управляющие признаки, которые делятся на три группы: первая — характеризует вид исходных данных (например, исходные данные

заданы по

неравномерной

сетке),

вторая — модификацию

основной

программы,

третья — вид

выдачи

результативных

данных (узкая

печать, выдача на АЦПУ и т. д.).

 

 

 

Возможны д в а р е ж и м а и с п о л ь з о в а н и я

основных

программ:

 

 

 

 

 

 

1. О д и н о ч н ы й

р е ж и м,

когда работает одна программа,

выполняющая один этап обработки данных.

 

 

2. Г р у п п о в о й

(или пакетный)- . р е ж и м,

когда

работает

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

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

126


Если основные программы используются во втором режиме, то при этом результаты предыдущей основной программы являются

исходными данными для

следующей программы. Передача дан­

ных от одной программы

к другой ведется с помощью магнитной

ленты (МЛ).

 

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

решаемых

ими

конкретных

задач.

 

 

 

 

 

ОсШ — первоначальная

обработка

рейсов.

 

 

 

ОсП2 — уравнивание

опорных сетей.

 

 

 

 

Модификации:

ОсП2.1 — передача

абсолютных

значений

силы

тяжести в

рядовые пункты;

Ос1І2.2 — уравнивание

опорных

сетей

с весом Рц = 1.

 

 

 

 

 

 

 

 

 

ОсПЗ — вычисление аномальных значений

силы тяжести и

соста­

вление каталога

 

пунктов.

 

 

 

 

 

 

Модификации:

ОсПЗ.1 — перевод координат { х, у } т±

{

ф, А);

ОсПЗ. 2 — вычисление нормальных значений

силы тяжести

по

фор­

мулам Гельмерта

 

и Красовского; ОсПЗ.З — вычисление аномальных

значений

силы

тяжести

в

редукции

Фая;

ОсПЗ.4 — вычисление

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

постоянных, заданных пользователем,

плотностях.

 

 

ОсП4 — редуцирование

гравиметрических измерений

со

слож­

ного рельефа на плоскость ОТНОСИМОСТИ.

 

 

 

ОсП5 — интерполяция

(восстановление

исходной функции

в уз­

лах квадратной сети).

 

 

 

 

 

Модификации: ОсП5.1 — вычисление

горизонтальных

градиен­

тов функции в узлах квадратной сети;

ОсП5.2 — вычисление

векто­

ров исходной функции в узлах квадратной сети; ОсП5.3 — вычисле­

ние осредненных значений исходной функции в узлах

квадратной

сети; ОсП5.4 — вычисление вторых вертикальных

производных

исходных функций в узлах квадратной сети.

 

ОсП7 — трансформации потенциальных

полей.

 

Модификации: ОсП7.1 — вычисление

потенциальных полей

в верхнем полупространстве и остаточных аномалий (с розами про­ стирания изолиний); ОсП7.2 — вычисление вертикальной производ­ ной исходной функции (с розами простирания изолиний).

ОсП8 — устойчивое продолжение исходной функции на плоскости нижнего полупространства.

Модификации: ОсП8.1 — вычисление вертикальных производных на плоскостях нижнего полупространства; ОсП8.2 — сглаживание исходной функции; ОсПЗ.З — пересчет исходной функции в верхнее

полупространство;

ОсП8.4 — вычисление вертикальных

производ­

ных на плоскостях верхнего полупространства;

 

ОсП9 — решение

прямой задачи — вычисление Ѵг от

аномаль­

ных тел типа контактной поверхности на горизонтальную плоскость. Модификации: ОсГО.І — вычисление поправки за рельеф мест­ ности в средней ИЛИ дальней зоне по равномерной сети; ОсП9.2 — вычисление поправки за рельеф местности по неравномерной сети;

127



ОеПО.о — вычисление F 2 2 1 5 от тел на горизонтальную плоскость; ОсШ.4 — вычисление поправки за рельеф для Уг 2 „ по равномерной

сети; ОсІ19.5 — вычисление поправки за рельеф для Ѵгг9

по неравно­

мерной сети;

ОсП9.6 — вычисление Ѵг9 и Ѵгг9

от аномальных тел

на

заданную

негоризонтальную

поверхность;

ОсП9.7 — вычисле­

ние Ѵг9 и Ѵгг9

от ограниченных аномальных тел произвольной формы;

ОсШ.8 — вычисление

Ѵг9, Vzz9 от тел типа контактной

поверхности

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

местности

при известном законе изменения

плотности, т. е. вычисление функции при переменной

плотности.

 

В целях удобства освоения как подготовки информации для ОсП,

так

и их эксплуатации

каждая программа снабжена тестовым при­

мером.

 

 

 

 

 

Кроме того, в целях автоматической и быстрой подготовки число­ вых массивов для тестов, когда объемы их оказались слишком боль­ шими, были составлены три программы подготовки числовых мас­

сивов для тестов. С помощью

п е р в о й

т е с т о в о й

п р о -

г р а м м ы

считается и выдается на перфорацию числовой

массив

для ОсП.4,

содержащий

{х, у, z, Vz}9, где Ѵг9 (х, у, z) вычисляется

в узлах неравномерной

сети.

В т о р а я

т е с т о в а я

 

п р о ­

г р а м м а

готовит числовой

массив для ОсП5, который

 

содер­

жит {.г, у, У2 }, где Ѵг9 (х, у) вычислены в узлах неравномерной

сети.

Т р е т ь я

т е с т о в а я п р о г р а м м а

готовит числовой

мас­

сив для ОсП7 и ОсП8, содержащий г (s, s)}, где Ѵг9 (s, s) вычис­ лена в узлах равномерной сети.

При разработке программ был создан ряд стандартных подпро­ грамм, используемых в нескольких ОсП. Им присвоены номера от 0154 до 0161 (под такими номерами в ОсП идет к ним обращение):

СП-0154 — вычисление значений и координат изолиний. СП-0155 — вычисление поправок за лунно-солнечное притяжение. СП-0156 — вычисление экваториальных координат Луны. СП-0157 — вычисление экваториальных координат Солнца. СП-0160 — формирование адресов команд по параметрам (автор

А.А. Корнейчук).

СП-0161 — печать матрицы чисел с переводом из десятичной

системы в двоичную и нумерацией

по строкам матрицы (автор

А. А. Корнейчук).

 

Программы оформлены в системе

ИС-2 и приводятся в данной

монографии.

 

Кроме того, используются следующие программы из стандартной библиотеки системы ИС-2:

СП-0000 — программа обмена; СП-0001 — вычисление функ­ ции z = хУ; СП-0002 — перевод числа из десятичной системы в двоич­

ную; СП-0003 — вычисление функции ех;

СП-0004 — вычисление

логарифма; СП-0005 — вычисление функции

sin х; СП-0006 — вы­

числение функции aresin х; СП-0007 — решение системы линейных уравнений; СП-0010 — перевод числа из двоичной системы в деся­ тичную; СП-0011 — вычисление функции tg х; СП-0012 — вычисле­ ние функции arctg х; СП-0015 — печать восьмеричного и десятпч-

•128


ного материала с любым числом кодов; СП-0016 —- перфорация восьмеричного материала с любым числом кодов и контроль; СП-0027 — нечать материала с переводом из двоичной системы в деся­ тичную и сохранением двоичного материала в МОЗУ; СП-0040 — перевод материала из двоичной системы в десятичную и перфо­ рация чисел; СП-0042 — групповой перевод чисел и десятичной

системы

в

двоичную;

СП-0071 — вычисление функции

У={1^1>

СП-00.74 — перевод

градусов

в

радианы;

СП-0075 — перевод ра­

дианов в

градусы;

СП-0433 — генератор

выдач

на

АЦПУ

(автор

Б.

Пухлинский).

 

 

 

 

 

 

 

 

 

 

В связи

с тем,

что в

отдельных ВЦ

обращение

к СП-0007 не

стандартно,

ниже

приводим

обращение

к СП-0007,

используемое

в

наших

программах. Оно имеет вид

 

 

 

 

 

 

 

х —1

0

0

0

16

х

7501

7610

 

 

 

 

 

х

я і

0

я 3

52

п-\-т

0007

п

 

 

 

 

 

x-f-1

я1

0

я 3

52

а и

0000

ct l

 

 

где il — порядок системы; m — число правых частей; п (п + 2т) ^

^ 3700;

0007 —• номер данной СП;

а— адрес первого

элемента

расширенной матрицы системы, записанной по строкам; с х 1

— адрес

первого

элемента рещений, которые

будут записываться построчно.

 

Г Л А В А

X I

 

 

ПЕРВОНАЧАЛЬНАЯ ОБРАБОТКА РЕЙСОВ

 

При работе программы использѵются СП с номерами: 0000 0002; 0006; 0007; 0010; ООН; 0012;"0015; 0027; 0042; 0071; 0075 0155; 0156; 0157.

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

Первоначально должны обрабатываться опорные рейсы. После

обработки рейсов опорной сети проводится уравнивание

опорной

сети (см. гл. X I I ) по

ОсІ12. Результаты

счета

по ОсП2

(массив

{ПР ПК, gyB}) составят

массив М2 при обработке рядовых

рейсов.

Исходная информация для ОсПІ состоит из

четырех массивов.

1°. Первый массив информации (Ml) для каждого прибора со­

стоит из одной перфокарты:

 

 

 

1-ая строка — каа — коэффициент сползания

нуль-пункта гра­

 

виметра в мгл/ч (паспортный или получен­

 

ный в результате

псследоваппй прибора).

9 Заказ 76

129