Файл: Автоматизированная система обработки и интерпретации результатов гравиметрических измерений..pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 11.04.2024
Просмотров: 200
Скачиваний: 0
Год |
работ |
Организация |
|
|
|
|
Каталог |
гравішет |
||
|
|
|
Масштаб |
ст,емкц |
|
|||||
|
Прямоугольные |
коорди |
|
Географические |
координаты |
Высота |
||||
Л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