книги из ГПНТБ / Автоматизированная система обработки и интерпретации результатов гравиметрических измерений
..pdfГод |
работ |
Организация |
|
|
|
|
Каталог |
гравішет |
||
|
|
|
Масштаб |
ст,емкц |
|
|||||
|
Прямоугольные |
коорди |
|
Географические |
координаты |
Высота |
||||
Л1. |
наты в системе |
1942 г. |
|
|||||||
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
в Бал |
|
n//u |
|
|
У, |
Широта В |
|
Долгота |
L |
тийск. |
||
|
метры |
|
|
системе, |
||||||
|
|
метры |
град |
млн |
сек |
град |
міш |
сек |
метры |
|
|
|
|
|
|||||||
1 |
38000000 |
12400000.0 |
034 |
19 |
19.8 |
067 |
54 |
48.5 |
50.00 |
|
2 |
3S000500 |
12400000.0 |
034 |
19 |
21.4 |
067 |
54 |
48.5 |
50.80 |
|
3 |
38001000 |
12400000.0 |
034 |
19 |
23.0 |
067 |
54 |
4S.5 |
51.30 |
|
4 |
38001500' |
12400000.0 |
• 034 |
19 |
24,6 |
067 |
54 |
48.5 |
52.00 |
|
5 |
3800200. |
12400000.0 |
034 |
19 |
26.2 |
067 |
54 |
48.4 |
52.30 |
|
6 |
3800250.0 |
12400000.0 |
034 |
19 |
27.9 |
067 |
54 |
48.4 |
53.00 |
|
7 |
3800000.0 |
12400200.0 |
034 |
19 |
19.8 |
067 |
54 |
56.3 |
49.00 |
|
20" |
3800080.0 |
12400700.0 |
034 |
19 |
22,6 |
067 |
55 |
15.9 |
51.30 |
|
П р и м е ч а н и е . Вид таблицы соответствует выдаче данных на АЦПУ.
|
|
|
|
|
|
|
|
|
Т а б л и ц а 48 |
|
-Г-Г + |
03 |
0000 0010 0000 |
|
|
|
|
||||
+++ |
03 |
|
0000 |
0015 |
0001 |
Информация о магнитной лейте (МЗ) |
||||
+++ |
03 |
0000 0023 0000 |
|
|
|
|
||||
+++ |
01 |
|
200 |
000 |
001 |
|
|
|
|
|
+ |
+ -Г |
01 |
|
230 |
000 |
000 |
Информация о счете (Ml) |
|
||
_1 |
і !_ |
01 |
|
267 |
000 |
000 |
|
|||
|
|
01 |
|
100 |
000 |
000 |
|
|
|
|
+++ |
00 |
|
000 |
000 |
000 |
|
|
|
|
|
01 |
|
100 |
000 |
000 |
п[ = |
\—количество |
отбракованных |
точек |
||
1 1 1 |
|
|
|
|
|
|||||
+++ |
07 |
|
380 |
017 |
000 |
хіъ |
I |
|
|
|
+ + + |
08 |
|
104 |
005 |
000 |
jj16 |
\ Информация |
о браковаппоіі |
точке |
|
H |
1 г |
02 |
|
520 |
000 |
000 |
||||
1 |
1 т |
06 |
|
979 |
662 |
100 |
gu J |
|
|
|
ч—і—г |
|
|
|
|||||||
+++ |
02 |
|
230 |
000 |
ООО |
|
ni = 23 — количество чисел в массиве |
|||
01 |
|
141 |
999 |
531 |
|
|||||
|
|
|
|
|
|
|
Д с ? Б ( С Т І )І = — 1 , 4 1 9 |
Ы Г Л |
|
|
+ |
— |
00 |
|
233 |
017 |
621 |
Д { ? Б ( с г і ) і 5 = — 0 , 2 3 3 мел |
|
||
+++ |
01 |
|
216 |
857 |
728 |
Д # Б ( f f ib4= 2,168 мгл |
|
|||
+-+ |
00 ' |
834 |
293 |
721 |
Д - £ Б |
(о"і)і7= -0,834 мгл |
|
|||
+++ 01 |
|
153 |
909 |
002 |
Д ^ Б ( 0 і Ь з = 1,539 мгл |
|
||||
++- |
01 |
|
677 |
542 |
594 |
Д{ТБ |
(°"З)23 = 0,067 мгл |
|
||
точек (точек из другой зоны) и соответствующие этим точкам тетрады десятичных чисел x, у, H, gH или ср, %, Я , gH.
Время счета (t в секундах) по ОсПЗ при решении основной за дачи определяется по эмпирической формуле t = пх, где пх — ко личество тетрад массива М2. Время счета теста массива М2 (табл. 45)
|
|
|
|
Т а б л и ц а 47 |
||
рпчеекпх пунктов |
|
|
Лист |
|
||
|
Автор отчета |
|
|
|
||
|
Нормальное |
|
Аномалия |
Аномалия |
Аномалия |
|
Ускорение |
значение |
|
||||
Анома |
Буге |
Буге |
Буге |
|||
силы тяжести |
силы тяжести |
|||||
в ГІотсдамск. |
по формуле |
лия |
Плотность |
Плотность |
Плотность |
|
системе, мгл |
Гельмерта |
Фая, мгл |
2.00 г/см куб, |
2,30 г/см чуб, 2,07 г/см куб, |
||
|
19011909, мгл |
|
мгл |
мгл |
мгл |
|
|
|
|
|
|
||
979660.000 |
979672.650 |
2.771 |
—1.420 |
-2.049 |
—2.824 |
|
979660.200 |
979672.688 |
3.180 |
—1.078 |
—1.717 |
—2.505 |
|
979660.401 |
979672.727 |
3.496 |
—0.804 |
—1.449 |
—2.244 |
|
979660.401 |
979672.765 |
3.674 |
—0.684 |
—1.338 |
-2.145 |
|
979660.301 |
979672.802 |
3.629 |
-0.755 |
—1.413 |
-2.224 |
|
0.000 |
979672.839 |
|
|
|
|
|
979663.502 |
979672.652 |
5.962 |
•1.855 |
1.239 |
0.479 |
|
979660.700 |
979672.716 |
3.805 |
—0.494 |
—1.139 |
—1.935 |
|
приблизительно равно 25 с. При решении модификаций основной задачи время счета меньше приведенного. При «авосте» рекомендуется делать повторный просчет, введя заново предварительно проверен ный пакет 3. Пакет 1 вводится в машину один раз, независимо от информации о счете и количества обрабатываемых массивов М2. Пакеты 2 и 3 собираются и вводятся в машину для каждой шестерки массивов М2 (если все шесть массивов М2 обрабатываются по единой информации о счете).
П р о г р а м м н ы е о с т а н о в ы |
в процессе счета: 1) в ячейке |
0403 по команде 77 0000 7777 0000; |
закончен счет по задачам: |
вычисление координат, вычисление нормальных значений силы тя жести, вычисление аномальных значений силы тяжести в редукции
Фая и Буге, |
составление массивов триад десятичных |
чисел п {х, |
у, Д#ф} или п |
{x, у, àgB (о-,)}; 2) в ячейке ИС-2 при выдаче массива |
|
результативных данных на перфорацию по СП-0016 |
(перфорация |
|
с контролем). В последнем случае необходимо поставить в ЧУ (чи тающее устройство) отперфорированный массив без первой карты (число кодов в результативном массиве) и нажать кнопку «Пуск». При совпадении контрольных сумм продолжается счет в соответ ствии с решаемой задачей. Если контрольные суммы при чтении не совпали, нажать кнопку «Пуск», и результаты (без карты с числом кодов) будут выданы на перфорацию еще раз. Процесс повторить до совпадения контрольных сумм. При двух-трехкратном повторении сообщить оператору ВЦ о неустойчивой работе перфоратора или читающего устройства.
По окончании счета массив результативных данных при выдаче в виде триад п {х, у, Ago} или п {х, у, А^Б(О",-)} записывается на МЛ. При решении других модификаций задачи результаты выдаются на АЦПУ, или на узкую бумажную ленту, или на перфорацию и на уз кую бумажную ленту.
148
149
Иа АЦПУ массив результативных данных выдается в виде стан дартных листов каталога гравиметрических пунктов. На перфорацию и узкую бумажную ленту массив результативных данных выдается при решении модификаций основной задачи в десятичном виде. Тестовые примеры выдачи массива результативных данных при ре шении основной задачи и одной из ее модификаций (табл. 40) пока заньг соответственно в табл. 47 и 48.
Порядок печати результатов по остальным тестовым примерам массива M l аналогичен печати результатов, приведенной в табл. 48.
Г Л А В А X I V
РЕДУЦИРОВАНИЕ ГРАВИМЕТРИЧЕСКИХ ИЗМЕРЕНИЙ СО СЛОЖНОГО РЕЛЬЕФА НА ПЛОСКОСТЬ ОТИОСИМОСТИ
В программе используются СП с номерами: 0000; 0002; 0007; 0027; 0042; 0154; 0160; 0161.
1.ПОДГОТОВКА ИСХОДНОЙ ИНФОРМАЦИИ
На площади гравиметрической съемки выбирается участок пря моугольной формы для редуцирования. Максимальные размеры участка 45sX 45s, где s — шаг съемки, а число шагов по осям х и у должно быть нечетным. Подготавливаются три информационных массива.
1°. |
На отдельной |
перфокарте |
подготавливается информация |
|
о МЛ (Ml), где в восьмеричном виде |
записывается следующее: |
|||
1-ая |
строка — номер магнитофона и номер первой зоны, начиная |
|||
|
с которой расположены в пулевом виде четыре |
|||
|
блока |
программы. |
|
|
2-ая |
строка — номер магнитофона и номер первой зоны, начиная |
|||
|
с которой расположены в настроенном виде четыре |
|||
|
блока |
программы. |
|
|
3-я |
строка — номер |
магнитофона и номер зоны, где записаны |
||
|
исходные |
данные. |
|
|
4-ая |
строка — номер магнитофона и номер первой зоны, начиная |
|||
|
с которой будут записаны программа и числовые |
|||
|
массивы |
при перерыве счета. |
||
5-ая строка — номер магнитофона и номер зоны для записи результата.
Перед началом счета нужно разметить МЛ на зоны следующих размеров: 0700, 7010, 0170 и 0470 кодов, для записи программы в ну левом виде; 0700, 7010, 0170, 7010, 4000, 4000, 7730 кодов для записи
|
|
|
|
программы в настроенном |
||
|
|
Т а б л и ц а |
49 |
виде |
и числовых |
массивов; |
000 |
0000 |
0023 |
0000 |
7730 кодов для записи исход |
||
002 |
0000 |
0001 |
0000 |
ных данных; 7740,7740 и 7740 |
||
001 |
0000 |
0001 |
0000 |
кодов для записи |
программы |
|
002 |
0000 |
0010 |
0000 |
при |
перерыве счета; 4000 ко |
|
002 |
0000 |
0014 |
0000 |
|||
007 |
0000 |
0051 |
0000 |
дов |
для записи |
результата. |
150
Вид M l приведен в табл. 49, где в КОП записывается номер магни-
тофона, |
а в |
А2 — номер |
зоны. |
|
|
|
|
|
|
|
|
|
|
|
|
|||
2°. Второй массив (М2) состоит из параметров счета, и в нем за- |
||||||||||||||||||
даются следующие |
данные: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
1-ая |
строка — |
s |
|
шаг результативной функции в |
км. |
|
||||||||||||
2-ая |
строка — |
M |
- размер |
исходного |
поля |
в |
сантиметрах |
|||||||||||
|
строка — |
N |
|
карты |
по |
строке. |
|
|
|
|
|
|
|
|
|
|||
3-я |
- размер |
исходного |
поля |
в |
сантиметрах |
|||||||||||||
|
|
|
àx |
|
карты |
по |
столбцу. |
|
|
|
|
|
|
х, |
у |
|||
4-ая |
строка — |
- точность |
определения |
|
координат |
|||||||||||||
|
|
|
|
|
в |
км. |
|
|
|
|
|
|
|
|
|
|
|
|
5-ая строка — |
àz |
- точность определения координаты z в км. |
||||||||||||||||
6-ая |
строка — |
|
|
координаты левого |
верхнего |
угла |
пло |
|||||||||||
7-ая строка — |
У о I |
|
||||||||||||||||
|
щади редуцирования |
в |
км. |
|
|
|
|
|
||||||||||
8-ая |
строка — |
гтіп |
- минимальное |
значение |
высоты |
|
на |
пло |
||||||||||
|
|
|
|
|
щади |
в |
км. |
|
|
|
|
|
|
|
|
|
|
|
9-ая |
строка — |
|
- минимальное |
значение |
|
Aga |
на |
нлощади |
||||||||||
|
|
|
|
|
редуцирования в |
мгл. |
|
|
|
|
|
|
Aga |
|||||
10-ая |
строка — |
8 |
- точность |
определения |
аномалий |
|
||||||||||||
11-ая |
строка — |
хг \ |
|
в |
мгл. |
|
|
|
|
|
|
|
|
|
|
|
||
|
координаты |
первой |
|
результативной |
||||||||||||||
12-ая строка — |
|
|
|
|||||||||||||||
|
|
точки |
(левый |
верхний |
угол) |
в |
|
км. |
|
|
||||||||
13-ая строка — |
|
- высота |
плоскости |
редуцирования. |
|
|
||||||||||||
14-ая строка — |
m |
- число |
результативных |
точек |
по |
строке. |
||||||||||||
15-ая строка — |
il |
- число результативных точек по столбцу. |
||||||||||||||||
16-ая |
строка — |
Сг |
- сечение результативной карты Ѵг в мгл. |
|||||||||||||||
17-ая |
строка — |
t |
= 1 при |
высокоточных |
съемках |
|
|
|
|
|||||||||
|
: 2 при съемках средней и низкой точности. |
|||||||||||||||||
|
|
|
|
|||||||||||||||
Пример записи для теста исходных данных М2 приведен в табл. 50. |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
Т а б л и ц а |
50 |
|||||
+++ |
00 |
500 |
|
000 |
|
000 |
|
|
|
|
|
|
s = |
0,5 |
км |
|||
02 |
310 |
|
000 |
|
000 |
|
|
|
|
|
|
M = |
31 |
см |
||||
4+++- + + |
02 |
110 |
|
000 |
|
000 |
|
|
|
|
|
|
N = |
11 |
см |
|||
++- |
01 |
100 |
|
000 |
|
000 |
|
|
|
|
|
|
A* = |
0,01 |
|
|||
03 |
100 |
|
000 |
|
000 |
|
|
|
|
|
|
Дг |
= |
0,0001 |
||||
-ь++ |
01 |
250 |
|
000 |
|
000 |
|
|
|
|
|
|
г 0 |
= |
2,5 |
км |
||
++- |
|
|
|
|
|
|
|
|
||||||||||
+-+ |
01 |
750 |
|
000 |
|
ООО |
|
|
|
|
г |
|
|
= |
—7,5 |
км |
||
4-4-4- |
00 |
000 |
|
000 |
|
000 |
|
|
|
|
|
т і п = |
0 |
|
|
|||
+++ |
00 |
000 |
|
000 |
|
000 |
|
|
|
|
Д£тіп = |
0 |
|
|
||||
00 |
600 |
|
000 |
|
000 |
|
|
|
|
|
|
e = |
0,6 |
мгл |
||||
+++ |
01 |
150 |
|
000 |
|
000 |
|
|
|
|
|
|
Zi = |
1,5 |
км |
|||
+ - Н - |
01 |
500 |
|
000 |
|
000 |
|
|
|
|
|
|
іГі = |
—0,5 |
|
|||
+-+ |
|
|
|
|
|
|
|
|
км |
|||||||||
00 |
750 |
|
000 |
|
000 |
|
|
|
|
|
|
го = |
0,75 |
|||||
+++ |
02 |
210 |
|
000 |
|
000 |
|
|
|
|
|
|
m = |
21 |
|
|
||
+++ |
01 |
700 |
|
000 |
|
000 |
|
|
|
|
|
|
ге = 7 |
|
|
|||
+++ |
01 |
200 |
|
000 |
|
000 |
|
|
|
|
|
|
|
= |
0,2 |
|
|
|
4 - + + |
01 |
20Ü |
|
000 |
|
000 |
|
|
|
|
|
|
Т) = |
2 |
|
|
||
4-4-+ |
|
|
|
|
|
|
|
|
|
|
||||||||
151
3°. Числовой массив исходных данных готовится но каталогу гравиметрических пунктов. На бланки переписываются х, у, z, Aga в десятичном виде, как указано в табл. 51.
Т а б л и ц а 51
+++
+++
+++
+-+
+++
+++
++-ь
02 |
165 |
600 |
хі |
= |
16,56 км |
02 |
366 |
300 |
г/і = 36,63 км |
||
00 |
804 ' |
700 |
=1 = 0,8047 км |
||
03 |
203 |
380 |
Д ^ = —203,38 мгл |
||
02 |
351 |
000 |
я 2 |
= |
35,1 км |
01 |
S4Ç |
000 |
.!/2 |
= |
8,46 км |
00 |
379 |
700 |
=г = 0,3797 км |
||
03 |
110 |
580 |
Ag„ = —110 58 мгл |
||
Подготовленные числовые |
данные |
у, |
z, |
Ag„} делятся на |
несколько массивов, в каждом из которых |
должно быть ио более |
|||
500 пунктов наблюдений. В |
каждом массиве |
на |
отдельной перфо |
|
карте указывается число пунктов наблюдений в массиве. В послед нем массиве числу пунктов присваивается признак.
Рис. 31. Изолинии в рельефа в км (1) и исходного поля в мгл (2) для теста ОсП4.
Все массивы должны иметь свои контрольные суммы.
4°. В предыдущих главах для опробования ОсП приводились тестовые примеры (информация.о МЛ, информация о счете и числовые массивы) в том виде, в котором они могут использоваться в про грамме. Для данной ОсП4 числовой массив (МЗ) оказался слишком обширным, и для его составления используется небольшая вспомога тельная программа, описанная в гл. X I X , разделе 3 (рис. 31). С по-ѵ мощью ее выдается массив {х, і/, z, àg) на перфокарты, который используется в качестве МЗ для ОсП4. Информация о МЛ (Ml) подготавливается пользователем; М2 приведен в табл. 50.
152
2. СЧЕТ НА ЭВМ
Перед началом счета необходимо разметить МЛ-0 на 16 зон по числу кодов, указанному в начале данного раздел_а, либо выделить 16 зон, указанных или больших размеров в ранее размеченной ленте.
Необходимо проверить, чтобы |
были записаны СП-0007, СП-0154, |
|
СП-ОКО, СП-0161. |
|
|
1°. Собираются три пакета |
перфокарт: |
|
Пакет 1 — для записи на МЛ. Порядок |
сбора перфокарт: 1; 2; |
|
чистая; M l ; чистая; 3—39; чистая; 40—118; |
чистая; 119—128; чи |
|
стая; 129-154. |
|
|
Пакет 2 — для настройки программ к счету. Порядок сбора пер фокарт: «вызов ИС-2», чистая; 155; чистая; M l ; чистая; М2.
Пакет 3 — для ввода числового массива и начала счета. Порядок сбора перфокарт: «вызов ИС-2»; чистая; 156; чистая; M l ; чистая; первый числовой массив МЗ; чистая; второй числовой массив МЗ; чистая; . . .; последний числовой массив МЗ.
Пакет 1 ставится в ЧУ и вводится в машину. Все подпрограммы последовательно запишутся в нулевом виде на МЛ. Останов по
команде |
77 в ячейке |
0020. |
|
|
|
Пакет |
2 ставится |
в ЧУ и вводится в машину. Печатается М2, |
|
все подпрограммы в течение 2 мин последовательно |
настраиваются |
|||
по |
параметрам и записываются на МЛ. Останов |
по команде 77 |
||
в |
ячейке 0151. |
|
|
|
|
Пакет 3 вводится в машину. Последовательно вводятся числовые |
|||
массивы, обрабатываются и начинается счет. В течение 3—4 мин происходит вычисление ср0 и результат записывается на МЛ в не сколько зон. Печатается результативная матрица по СП-0161. Сле
дует просмотреть |
для контроля запись значений ф 0 . |
|
||||
Результативная матрица представляет собой массив чисел, со |
||||||
стоящий из + |
+ |
+ |
11 343 597 000 (код из ячейки |
7752) и значе |
||
ний ф 0 |
. Условный код находится в узлах, в которых нет измеренных |
|||||
(х> |
У)- Значения ф 0 должны плавно изменяться от точки к точке. |
|||||
В ином случае следует еще раз проверить числовой массив. |
||||||
2°. После печати матрицы ф 0 вызывается подпрограмма вычисле |
||||||
ния ф г и начинается |
процесс расчета ф , , который |
обычно бывает |
||||
очень |
длителен. |
Время одного приближения оценивается как t = |
||||
= 0,57??іг -10~3 |
мин, где m = MN — число |
исходных точек. |
||||
В процессе |
работы ЭВМ после расчета |
каждой |
результативной |
|||
строки нулевой и первый кубы записываются на МБ. Счетчик точек по строке находится в ячейке 0075, счетчик по столбцу в ячейке 0074.
Когда будут вычислены ф , во всех точках, печатается результативпая матрица, ЭВМ останавливается в ячейке 0222 по команде 77.
После нажатия кнопки «Пуск» |
результат записывается на МЛ. п |
счет перейдет на подпрограмму |
Ф = 2 Фі- |
В этой подпрограмме за 1—2 мин будут вычислены и напечатаны
А заданное, А Г = max | ф ; — ф і _ і | и А2 , равная среднему | фг- — ф і _ і |, по которым можно следить за скоростью сходимости итерационного
153
процесса. Если значение Д х велико, то будет продолжен итерацион ный процесс по вычислению ф; . Если значение Аг меньше заданной точности, то счет перейдет к вычислению 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