Файл: Автоматизированная система обработки и интерпретации результатов гравиметрических измерений..pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 11.04.2024
Просмотров: 203
Скачиваний: 0
процесса. Если значение Д х велико, то будет продолжен итерацион ный процесс по вычислению ф; . Если значение Аг меньше заданной точности, то счет перейдет к вычислению Vz (z0 ) на плоскости от носимости.
По желанию пользователя, после достижения необходимой точ ности, ЭВМ перейдет к вычислению Vz или в нужный момент — в при нудительном порядке. В последнем случае ЭВМ останавливают и вводят перфокарты: «вызов ИС-2»; чистая; 157; чистая; M l .
В конце счета выдается Vz на плоскости относимости в узлах квадратной сети и печатаются шгСП-0154 значения и координаты изолиний. Останов происходит в ячейке 1112 по команде 77 после
записи результата. |
|
|
Время вычисления Vz |
определяется как t = 0,6m-п-10~3 |
мин, |
где m — число исходных |
точек, а п — число результативных |
точек. |
При сбое машины в момент вычисления ср; или Ѵг вводят перфо карту 158. Она перепишет с МБ-0 содержимое нулевого и первого МОЗУ и счет продолжится.
При нехватке времени останавливают ЭВМ и вводят пакет пере рыва счета — перфокарты: 159; 160; чистая; M l . Содержимое МБ-0 переписывается на МЛ.В этом случае для продолжения счета вводят пакет перфокарт: 160; 162; чистая; M l . Содержимое нулевого и пер вого кубов перепишется с МЛ на МБ, и счет продолжится.
При расчетах по тесту для Ос114 будут напечатаны данные, при веденные нами в табл. 50 и 52. Время счета по тесту составит 3 ч 48 мин.
|
|
|
|
|
|
|
Т а б л и ц а 52 |
|
' |
_І_ |
01 |
548 |
926 |
000 |
|
|
|
~г~і— |
01 |
620 |
174 |
001 |
|
То |
||
t |
1 |
|
|
|
|
|
||
01 |
530 |
963 |
030 |
|
|
|||
1 i |
|
|
||||||
- t - + - |
01 |
151 |
399 |
000 |
|
|
||
++- |
01 |
335 |
229 |
001 |
|
Фі |
||
++- |
01 |
138 |
S74 |
030 |
|
|||
Д |
заданное |
|||||||
+ + - |
•Ol |
749 |
999 |
999 |
||||
+ - + |
01 |
122 |
447 |
156 |
Ді |
|||
+ + - |
00 |
254 |
230 |
897 |
Д |
2 |
||
+ |
— |
02 |
289 |
911 |
000 |
Фа |
||
-++ |
01 |
415 |
089 |
323 |
|
|
||
-++ |
01 |
406 |
410 |
379 |
|
|
||
- + + |
01 |
800 |
00 |
Значения |
и координаты изолішии Vг |
|||
- + + |
00 |
974 |
|
|||||
18 |
019 |
|
|
|||||
154
Г Л А В А X V
ИНТЕРПОЛЯЦИЯ (ВОССТАНОВЛЕНИЕ ИСХОДНОЙ ФУНКЦИИ В УЗЛАХ КВАДРАТНОЙ СЕТИ)
При работе программы используются СП с номерами: 0000; 0002; 0007; 0010; 0027; 0075; 0154.
1.ПОДГОТОВКА ИСХОДНОЙ ИНФОРМАЦИИ
1°. Исходным материалом для задачи служат схема расположения пунктов наблюдений и каталог гравиметрических пунктов. В послед нем представлены значения аномалии силы тяжести и прямоуголь ные координаты точек наблюдений &.gK, хк, ук. Обрабатываемая площадь должна состоять из целого числа (по столбцу и по строке) квадратов (рис. 32) размером 16sXl6s, где s — шаг результатив ных значений, равный 1 см карты. Квадратам присваивают поряд ковые номера слева направо и сверху вниз, начиная с нуля. В каждом квадрате должно быть не менее 40 точек. Задают прямо угольные координаты х0, у0 левого верхнего угла всей площади.
Следует помнить, что применяется прямоугольная система коор
динат, |
принятая в геодезии, |
а именно: положительное направление |
||||||||
оси X — на север, положительное .направление оси у — на |
восток. |
|||||||||
2°. Подготавливаются три информационных массива: информация |
||||||||||
о магнитной |
ленте — M l ; информация |
о |
счете |
— М2; числовой |
||||||
массив |
— |
МЗ. |
|
|
|
|
|
|
|
|
Карта |
информации о МЛ (Ml) содержит |
следующие данные, за |
||||||||
писываемые на отдельной перфокарте |
в |
восьмеричном виде: |
||||||||
1-ая |
строка — номер магнитофона и номер зоны, где расположена |
|||||||||
|
|
|
программа |
в |
нулевом |
состоянии. |
|
|||
2-ая |
строка — номер магнитофона и номер первой зоны, |
начиная |
||||||||
|
|
|
с которой будут последовательно записаны неупо |
|||||||
|
|
|
рядоченные |
триады |
исходных |
данных. |
|
|||
3-я |
строка — номер магнитофона и номер первой зоны, |
начиная |
||||||||
|
|
|
с которой будут последовательно записаны |
упоря |
||||||
|
|
|
доченные триады. |
|
|
|
|
|
||
4-ая |
строка — номер магнитофона |
и номер зоны, начиная с ко |
||||||||
|
|
|
торой будут записаны массивы результативных |
|||||||
|
|
|
данных. |
|
|
|
|
|
|
|
5-ая |
строка — номер магнитофона и номер первой зоны, на иная |
|||||||||
|
|
|
с которой будут храниться ОсП и исходные данные |
|||||||
|
|
|
при перерыве |
счета. |
|
|
|
|
||
Поэтому |
перед счетом необходимо |
разметить минимум две МЛ |
||||||||
в соответствии с объемом исходных данных. Рекомендуется следу ющая разметка: МЛ-0 — 1 зона — 7200 кодов для ОсП в нулевом виде, 2, 3 и 4 зоны — по 7777 кодов для хранения ОсП и числовых
массивов при перерыве счета, от 5 до 5 + |
К зоны — по 2000 |
кодов |
для .хранения неупорядоченных массивов |
исходных данных, |
от 5 |
155
до К + |
1 и от 5 т Я т MN зоны — по 410 кодов для записи масси |
||||||||||||||||
вов результативных |
данных; МЛ-1 от 1 до MN зоны — по 1410 ко |
||||||||||||||||
|
|
|
|
|
|
|
|
дов |
для хранения |
упорядо |
|||||||
|
|
|
Т а б л и ц а 53 |
|
ченных |
матриц. |
|
|
о МЛ |
||||||||
ООО |
0000 |
0001 |
|
0000 |
|
|
|
|
Вид |
информации |
|||||||
|
|
|
|
с учетом табл. 54 при указан |
|||||||||||||
ООО |
0000 |
0005 |
|
0000 |
|
|
|
||||||||||
|
|
|
|
ном |
выше |
расположении |
|||||||||||
001 |
0000 |
0001 |
|
0000 |
|
|
|
||||||||||
ООО |
0000 |
0016 |
' 0000 |
|
|
|
на МЛ приведен |
в табл. 53, • |
|||||||||
ООО |
0000 |
0002 |
|
0000 |
K S |
|
|
где |
помер |
магнитофона |
ука |
||||||
001 |
оооо • |
0027 . |
0000 |
|
|
зан |
в |
КОП, а номер |
|
зоны |
|||||||
|
|
|
|
|
|
|
|
в А2. |
|
|
|
|
|
|
|
||
3°. Выбираются параметры счета из табл. |
11 и составляется ин |
||||||||||||||||
формация о счете (М2). |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
В М2 записываются |
следующие |
данные: |
|
|
|
|
|
|
|
||||||||
1-ая |
строка |
|
|
шаг |
результативных |
значений |
функции |
||||||||||
|
|
|
|
в км. |
|
|
|
|
|
|
|
|
|
|
|
||
2-ая |
строка - |
|
M |
- число квадратов исходного поля в строке, |
|||||||||||||
|
|
|
|
причем |
38 |
|
M s* 1. |
|
|
|
|
|
|
||||
3-я |
строка |
|
N |
число квадратов исходного поля в столб |
|||||||||||||
|
|
|
|
це. Следует иметь в виду, что алгоритмом |
|||||||||||||
|
|
|
|
наложено ограничение на форму исход |
|||||||||||||
|
|
|
|
ного |
поля: |
если |
M — 1, |
то |
N = 1, |
||||||||
|
|
|
|
если M > - 1 , то N ^ |
1, т. е. |
алгоритм |
|||||||||||
|
|
|
|
не рассчитан на обработку |
материалов, |
||||||||||||
|
|
|
|
когда исходные квадраты имеют вид |
|||||||||||||
|
|
|
|
вертикального столбца. В этом случае |
|||||||||||||
|
|
|
|
следует поменять координатные оси и |
|||||||||||||
|
|
|
|
столбец |
квадратов расположить |
|
гори |
||||||||||
|
|
|
|
зонтально. |
|
|
|
|
|
|
|
|
|
|
|||
4-ая |
строка |
|
|
точность определения |
аномалий для дан |
||||||||||||
|
|
|
|
ного |
масштаба |
съемки |
в мгл. |
|
|
|
|||||||
о-ая |
строка |
е2 /е3 |
отношение |
|
величин |
погрешностей |
на |
||||||||||
|
|
|
|
блюдений рядовых и опорных точек. |
|||||||||||||
|
|
|
|
Если |
точек |
опорной |
сети нет, то в 5-й |
||||||||||
|
|
|
|
строке |
записывается |
единица, |
|
|
|
||||||||
6-ая |
строка — |
Сг |
сечение |
изолиний |
результативной |
карты |
|||||||||||
|
|
|
|
в мгл или в этвешах, или 1 - Ю - 5 ед. СГС |
|||||||||||||
|
|
|
|
(при вычислении' Ѵ^) |
в зависимости от |
||||||||||||
|
|
|
|
решаемой |
задачи. |
|
|
|
|
|
|
|
|||||
7-ая |
строка — х0 |
|
координаты верхнего левого угла обра |
||||||||||||||
8-ая |
строка — у0 |
||||||||||||||||
батываемой |
площади |
в км. |
|
|
|
|
|||||||||||
9-ая |
строка — |
R |
безразмерный |
параметр, |
выбираемый |
||||||||||||
|
|
|
|
программой автоматически на основании |
|||||||||||||
|
|
|
|
средней плотности съемки. В строке на |
|||||||||||||
|
|
|
|
бланке записывается нуль. Если необхо |
|||||||||||||
|
|
|
|
димо провести вычисления с заранее |
|||||||||||||
|
|
|
|
выбранным |
радиусом, |
например |
при |
||||||||||
156
осреднении, то на бланке в 9-й строке записывается величина радиуса в см карты.
10-ая строка — т] |
| |
безмерные |
параметры, выбираемые на |
v |
' |
основании |
масштаба съемки и точности |
исходных данных (табл. 11). 12-ая строка — число массивов исходных данных.
Выбранные параметры счета составляют М2 и записываются в де сятичном виде по форме, показанной в табл. 54.
+++
+++
+++
+++
+++
+++
+++
+++
-++
+++
+++
|
|
|
|
|
Т а б л и ц а 54 |
00 |
500 |
000 |
|
000 |
5 |
01 |
300 |
000 |
|
000 |
M |
01 |
300 |
000 |
|
000 |
N |
00 |
400 |
000 |
• |
000 |
s |
01 |
100 |
000 |
|
000 |
s2 /e3 |
01 |
100 |
000 |
|
000 |
Ci |
02 |
240 |
000 |
|
000 |
x0 |
00 |
000 |
000 |
|
000 |
Уо |
01 |
300 |
000 |
|
000 |
R |
01 |
100 |
000 |
|
000 |
4 |
01 |
200 |
000 |
|
000 |
V |
01 |
900 |
000 |
|
000 |
MN |
Следует иметь в виду, что величина е участвует только в отбра ковке исходных ошибочных точек (V.8) и только посредством этого влияет на результат. Очень часто в соотношении (V.8) величина Л чуть больше Зе и точка попадает в ошибочные. При расчетах по сильно дифференцированным полям рекомендуется указывать в ин формации о счете величину s значительно большую, чем рекомендо вана в [97]. Это иногда улучшает результат, так как избегают жесткой браковки многие точки, лежащие на пределе точности съемки.
4°. При подготовке массива исходных данных (МЗ) из |
каталога' |
||
на бланки выписываются триады чисел хк, |
ук, AgK для |
опорных |
|
точек, а затем хк, |
і/к , AgK для рядовых точек. Триады чисел выписы |
||
ваются в любой |
последовательности. |
|
|
Обычно значения хк, ук имеют в каталоге |
7—8 значащих цифр. |
||
Часть старших разрядов координат точек наблюдений, постоянную для всей площади, можно не вписывать, так как вычисления прово дятся в условной системе координат. Например, если в каталоге записано
^ = 35 678,975
хг^ 35 677,384 х к = 35 687,735
157
где три старших разряда постоянны, то координаты можно записать так:
^ = 78,975; жа = 77,384; . . .; агкв = 87,735.
Координаты и значения функции выписываются па бланки в де сятичном виде по одному числу в ячейку (табл. 55). Координаты должны быть выражены в км, Ag — в мгл.
|
|
|
|
|
Т а б л и ц а 55 |
||
+ + + |
04 |
757 |
690 |
960 |
|
|
|
04 |
801 |
370 |
450 |
г/к |
K |
|
|
_ |
02 |
980 |
410 |
000 |
Ag |
опорная |
|
+ +! +і |
|
||||||
+ + І |
04 |
757 |
691 |
370 |
|
|
|
- + + |
04 |
SOI |
369 |
430 |
Ук+1 |
|
|
+ + - |
02 |
979 |
460 |
000 |
Agic+i опорная |
||
+ + + |
04 |
757 |
693 |
130 |
*п |
|
|
+ + + |
04 |
SOI |
380 |
100 |
Ѵп |
|
|
+ + + |
01 |
321 |
500 |
000 . |
|
|
|
Значениям AgK опорных точек присваивается признак (табл. 55). |
|||||||
Весь числовой массив исходных данных для счета |
подразделяется |
||||||
на несколько массивов по 75 перфокарт в колоде (300 пунктов на блюдений).
На отдельных перфокартах пробивается количество пунктов на блюдений каждого массива. Должно быть 300 пунктов в каждом и менее 300 в последнем массиве. В первом обязаны быть все опорные точки. Образованные массивы исходных данных перфорируют и
суммируют. |
|
|
|
|
|
|
|
Порядок сбора перфокарт одного массива исходных данных (МЗ) |
|||||
следующий: первая перфокарта с количеством |
пунктов |
наблюдений |
||||
в |
массиве, |
вторая и последующие — трпады |
чисел |
(хк, ук, àgK) |
||
в |
любой последовательности; |
последняя |
перфокарта — контроль |
|||
ная сумма |
массива. |
|
|
|
|
|
|
5°. С учетом возможностей ОсП, решающей как основную задачу |
|||||
(вычисление |
функции в узлах |
квадратной |
сетп), так и ряд ее моди |
|||
фикаций, была выбрана система признаков, с помощью которых управление автоматически передается на блок программы, настраи
вающий ОсП на вычисления по одной из |
задач. |
|
|
||
В общем случае, |
когда н и |
о д н о м у |
и з |
ч и с е л |
в п е р |
ф о к а р т е и н ф о р м а ц и и |
о с ч е т е |
(М2) н е |
п р и |
||
с в а и в а е т с я |
п р и з н а к , |
вычисляются |
значения |
функции |
|
в узлах квадратной сети (ОсП5). При решении этой задачи исходная информация задается в виде координат функции и значений функции в узлах неравномерной сети. В результате работы программы получаются значения функции в узлах квадратной сети и выдаются на печать значения изолиний функции и координаты изолиний.
Ниже приводятся описание модификаций этой задачи и номера строк в М2, в которых ставится признак (минус в 45 разряде) для
158
