Файл: Автоматизированная система обработки и интерпретации результатов гравиметрических измерений..pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 11.04.2024
Просмотров: 209
Скачиваний: 0
+• 1 |
|
01 |
499999999 |
' |
|
'г — |
|||||
+ J—L |
Ol |
499999999 |
|
||
00 |
• 659745359 ' |
||||
1 |
1 |
||||
-+— |
Ol |
999999999 |
|
||
+•-j |
|
Ol |
999999999 |
|
|
•-] |
|
00 |
495613927 |
|
|
|
Ol |
|
|||
++-f - |
|
|
624999999 |
|
|
+ -Т ~ |
|
Ol |
624999999 |
|
|
+-+• + |
Ol |
129315396 |
|
||
1 |
1 |
02 |
24000 |
|
|
+ 1 |
1 |
|
|||
+- 1 _ |
j |
Ol |
65789 OOS |
+02 10681 ООО
+ -+•+ |
02 |
12012 ООО |
|
+ |
|
02 |
24000 |
~ |
г |
Ol |
71557 007 |
1 |
|
00 |
95597 015 |
+ J |
|
||
|
Ol |
499999999 |
|
-1 |
|
Ol |
499999999 |
+• |
|
00 |
175000000 |
1 |
|||
+•1 |
— |
|
|
- Ь |
Ol |
599823807 |
|
01 |
499999999 |
|
01 |
499999999 |
|
00 |
659745359 |
|
01 |
999999999 |
|
01 |
999999999 |
|
00 |
495613927 |
++- |
01 |
624999999 |
01 |
624999999 |
|
|
01 |
129315396 |
|
03 |
10000 |
-++ |
01 |
66676 008 |
|
|
|
-++ |
01 |
96411 ООО |
03 |
10000 |
|
|
01 |
72147 007 |
|
01 |
85048 007 |
-+ - 01 499999999
01499999999 ++ 00 17500000000 47945 009
Т а б л и ц а 71
s, s - i
Изолиния |
Интерполяция |
|
|
в мгл |
|
||
по X |
Регулярпзи- |
||
Координата х |
|||
|
рованныѳ |
||
|
|
||
|
|
приближения |
|
Изолиппя |
|
на 1-ом |
|
Интерполяция |
уровне |
||
в мгл |
|||
по у |
|
||
Координата у |
|
||
|
|
Т а б л и ц а 72
as, ßs . es
Изолиния |
|
Интерполяция |
|
|
в этвешах |
|
|
||
|
по |
X |
Регуляризи- |
|
Координата х |
|
|
рованное |
|
|
|
|
|
J приближение |
Изолиния |
|
Интерполяция |
на 1-ом |
|
в этвешах |
|
уровне |
||
|
по |
у |
|
|
Координата |
у |
|
||
|
|
|
— Н - |
01 825904176 |
Г Л А В А X V I I I
РЕШЕНИЕ ПРЯМОЙ ЗАДАЧИ - ВЫЧИСЛЕНИЕ Ѵ2
НА ГОРИЗОНТАЛЬНУЮ ПЛОСКОСТЬ ОТ АНОМАЛЬНЫХ ТЕЛ
ТИПА КОНТАКТНОЙ ПОВЕРХНОСТИ
В программе используются СП с номерами: 0004; 0012; 0027; 0040; 0042; 0154.
|
|
1. ПОДГОТОВКА |
ИСХОДНОЙ |
ИНФОРМАЦИИ |
|
||||
Параметрами счета по ОсП9 являются: шаг цифровой модели |
|||||||||
структуры |
(s); шаг таблицы |
высот |
(к); шаг таблицы |
плотностей |
|||||
(Іга); внешний и внутренний |
радиусы при вычислении |
поправок за |
|||||||
рельеф |
местности; размеры |
краев. |
|
|
|
|
|||
Исходная величина погрешности для выбора параметра е счи |
|||||||||
тается |
равной Ѵ4 величины допустимой |
ошибки результата. |
|||||||
Ш а г |
р а з б и е н и я |
с т р у к т у р ы |
(s) выбирается в пре |
||||||
делах sm i n |
=^ s < smax |
таким, |
чтобы |
по возможности выполнялись |
|||||
следующие |
условия: |
1) минимальное |
время счета; 2) получение ре |
||||||
зультата с заданной |
точностью; 3) минимальное время |
подготовки |
|||||||
исходных |
данных. |
|
|
|
|
|
|
|
|
Минимальный шаг (sm i n ) желательно |
брать равным 1 см в мас |
||||||||
штабе |
карты. Максимальная |
величина |
шага |
sm a x оценивается как |
|||||
smax = |
2>-, где 1- — рассчитывается по формуле |
|
|||||||
|
|
_ |
if |
2hi Ah + Ш-Фг |
~ |
|
|
||
|
|
7 |
- У |
|
4Ф |
|
Л і ' |
|
|
Д/г. — ошибка построения |
структурной |
карты; h ± — средняя глу |
|||||||
бина до верхней поверхности структуры; |
|
|
|||||||
|
|
|
|
Ф = О 4 |
АЛ. |
|
|
||
|
|
|
|
|
2п/а |
|
|
|
|
Для задачи вычисления поправки за рельеф местности шаг раз биения структуры выбирается минимальным, в частности, он может быть равным 1 см в масштабе карты.
М и н и м а л ь н ы й ш а г т а б л и ц ы в ы с о т определяется по формуле
г, - ^ г п а х ^ ш і п
где h w a x и h m i n — соответственно максимальная и минимальная от метки структурной карты. Если по формуле получается величина,
183
меньшая 1 / 3 сечения карты, то шаг принимается приблизительно равным Ѵз сечения карты, так как иначе цифровая модель структуры' будет задана с завышенной точностью.
М и н и м а л ь н ы й |
ш а г |
т а б л и ц ы |
п л о т н о с т е й |
||
(при расчетах с |
переменной плотностью) оценивается по формуле |
||||
|
|
|
64 |
» |
|
где стпшх, сттім — |
соответственно максимальная и минимальная] от |
||||
метки плотностей |
карты. |
|
|
|
|
Если по формуле получается, что величина шага меньше г / 3 сече ния карты, то шаг таблицы плотностей рекомендуется брать равным Ѵз сечения карты изоплотностей.
Р а з м е р ы |
к р а е в |
определяются по формуле |
|
|
||||
|
|
|
•hi) Ах |
, Да* |
|
|
|
|
где X — искомое |
расстояние |
от края результативного |
квадрата |
до |
||||
исходного (размер краев); |
Ах |
— расстояние |
от |
центра |
до края |
ре |
||
зультативного квадрата; hx |
— средняя глубина |
на |
рассматриваемом |
|||||
крае до поверхности аномалпеобразующего |
тела; |
h2 — глубина |
до |
|||||
нижней кромки |
тела. |
|
|
|
|
|
|
|
Прп решении задачи вычисления поправки за рельеф местности вместо оценки размеров краев проводится оценка внутреннего и внеш него радиусов. Внутренний радиус следует брать равным 1 см в мас штабе карты. Внешний радиус, ограничивающий область, в которой необходимо учитывать влияние рельефа на показания прибора, сле дует выбирать, исходя из характера местности и точности определе ния аномалий. Для расчетов внешнего радиуса авторы рекомендуют
использовать формулу, предложенную |
Н. И. Дергачевым: |
|
|
2ah'M' |
|
Л = |
2dôg |
» |
где а — средняя плотность пород дальней зоны; h' — среднее пре вышение высот дальней зоны над высотами локального поля; Ah' — приращение высоты от пункта к пункту в пределах площади учета поправки за рельеф; dbg — допустимое искажение аномального поля, вызванное недоучетом влияния рельефа дальней зоны.
|
1°. Первый массив (Ml) |
исходной информации к 0сП9 составляется |
из |
выбранных параметров |
счета и состоит из четырех перфокарт. |
На |
первых двух перфокартах (22 строки) записываются параметры |
счета в |
десятичном виде: |
|
1- ая |
строка — е — Ѵ4 величины допустимой ошибки результата |
|
|
в |
мгл. |
2-ая |
строка — к — шаг таблицы высот в км. |
|
3-ая |
строка — х0 |
\ — координаты первой точки счета относи |
мая |
строка — у0 |
J тельно выбранного начала координат в км. |
184
5-ая |
строка |
|
— глубина от плоскости результата до плоско |
|||||||||||
|
|
|
|
|
сти, принятой за нуль (ось z направлена |
|||||||||
6-ая |
строка |
|
|
|
вниз) в км. |
|
|
|
|
|
|
|
||
|
|
|
— шаг |
результативных |
точек по г |
в км. • |
||||||||
7-ая |
строка |
|
|
|
— шаг результативных точек по у в км. |
|||||||||
8- ая |
строка |
|
|
|
— шаг |
результативных точек по z в км. |
||||||||
9-ая |
строка |
|
go |
|
— координаты левого нижнего угла |
струк |
||||||||
10-ая |
строка |
|
40 |
|
туры |
относительно |
выбранного |
начала |
||||||
|
|
|
|
|
координат в км. |
|
|
|
|
|
|
|
||
11-ая |
строка |
— S| |
|
— шаг |
по £ |
шаг разбиения структуры в км. |
||||||||
12-ая |
строка |
— sn |
|
— шаг |
по т) |
|
|
|
|
по х. |
|
|
||
13-ая |
строка |
— |
nx |
|
— число результативных |
точек |
|
|
||||||
14-ая |
строка — nu |
|
— число результативных точек по |
у. |
|
|
||||||||
15-ая |
строка — n, |
|
— число плоскостей, на которых |
вычисляется |
||||||||||
|
|
|
|
|
гравитационный эффект от |
структуры, |
||||||||
16-ая |
строка — |
п0 |
|
число постоянных избыточных |
плотностей |
|||||||||
|
|
|
|
|
(единица, если ведутся расчеты |
с перемен |
||||||||
|
|
|
|
|
ной плотностью). |
|
|
|
|
|
|
|
||
17-ая |
строка |
— а, |
значения |
первой |
избыточной |
плотности |
||||||||
|
|
(или шаг таблицы плотностей) в г/см3 , |
||||||||||||
|
|
— |
|
|
|
|
|
|
|
|
|
|||
18-ая |
строка |
|
с» |
|
значение |
второй |
избыточной |
плотности |
||||||
|
|
|
|
|
(или минимальная плотность |
сттіп) в |
г/см3 , |
|||||||
19ая |
строка — а 3 |
|
— значение |
третьей |
избыточной |
плотности |
||||||||
|
строка — С |
|
в г/см3 . |
|
|
|
|
|
|
|
|
|||
20-ая |
|
— сечение |
изолиний |
результативной |
функ |
|||||||||
|
|
|
|
|
ции, |
если результат выдается в изолиниях |
||||||||
|
|
|
|
|
(при выдачи результата |
в узлах |
равномер |
|||||||
|
|
|
|
|
ной |
сети |
записывается |
нуль), |
в мгл |
или |
||||
|
|
|
|
|
в этвешах. |
|
|
|
|
|
|
|
||
21-ая |
строка — г |
|
— внутренний радиус |
для расчета поправок |
||||||||||
|
|
|
|
|
за рельеф местности (для всех |
остальных |
||||||||
|
строка — R |
|
задач — нуль) в км. |
|
|
|
|
|
|
|||||
22-ая |
— внешний радиус в км. |
|
|
|
|
|
|
|||||||
Третья перфокарта массива M l содержит К 2 |
первых |
перфокарт, |
||||||||||||
последняя перфокарта |
чистая. |
|
|
|
|
|
|
|
|
|
||||
Две первые перфокарты M l в зависимости от модификации |
ОсП9 |
снабжаются признаками («минус» в 45 разряде в определенной строке): в первой строке признак ставится при вычислении поправки за рельеф местности в средней или дальней зоне по равномерной сети для Vz (ОсПЭ.І), во второй — при вычислений поправки за рельеф в сред
ней или дальней зоне для |
Ѵг |
по неравномерной сети (ОсП9.2), |
в третьей — при вычислении |
Ѵгг |
на заданную горизонтальную пло |
скость (ОсПЭ.З), в первой и третьей строках (два признака) — при вычислении поправки за рельеф для Ѵгг по равномерной сети (ОсП9.4),
во второй и третьей — |
при вычислении поправки за рельеф для |
Ѵгг |
по неравномерной сети |
(ОсП9.5), в четвертой — при вычислении |
Ѵг |
иѴгг на заданную негоризонтальную поверхность (ОсП9.6), в пятой
ишестой — при вычислении Ѵг и Ѵгг от ограниченных аномальных
185
тел |
произвольной |
формы |
(ОсП9.7), в |
седьмой — при |
вычислении |
|
Ѵг, |
V'гг и поправки за рельеф от тел с переменной плотностью (ОсП9.8); |
|||||
в*восьмой — при перфорации результата (вид выдачи), в |
девятой — |
|||||
при накапливании |
промежуточного результата на МБ, в десятой — |
|||||
прп вычислении Ѵг |
и V'а с повышенной точностью с использованием |
|||||
стандартных |
программ, в одиннадцатой — при выдаче результата на |
|||||
печать в пзолпниях (вид выдачи). |
|
|
||||
|
В том случае, |
когда |
о т с у т с т в у ю т п р и з н а к и , р е |
|||
ш а е т с я о с н о в н а я |
з а д а ч а |
вычисления Ѵг от тел типа |
||||
контактной |
поверхности |
иа заданную |
горизонтальную |
плоскость |
||
(ОсП9). |
|
|
|
|
|
|
Vz |
Приводим примеры массивов M l для двух задач. Для вычисления |
|||||
на горизонтальной плоскости от заданного тела (ОсП9) две первые |
перфокарты M l показаны в табл. 73. Используемые при вычислении поправки за рельеф местности для Vz по неравномерной сети (ОсП9.2) две первые перфокарты M l приведены в табл. 74.
|
|
|
|
|
|
|
|
Т а б л и ц а 73 |
|
|
— |
Ol |
200 |
000 |
000 |
e = |
0,02 M I л |
||
- j |
|
|
02 |
500 |
000 |
000 |
к — —0,005 км |
||
Л _ |
J |
! _ |
00 |
|
|
|
2,0 = 0 |
||
; |
_ j |
t_ |
00 |
|
|
|
Уо = 0 |
||
|
|
|
|
|
|
||||
--+ |
00 |
250 |
000 |
000 |
го = —0,25 км |
||||
|
! L |
00 |
250 |
000 |
000 |
s^ = |
0,25 км |
||
_ t |
ji |
ij |
00 |
250 |
000 |
000 |
s,, = |
0,25 км |
|
|
I |
i |
00 |
125 |
000 |
000 |
sz = |
0 |
|
- — t - |
01 |
îo = -1,25 км |
|||||||
~ г |
" i — Г |
|
|
|
|
|
|
||
_ j |
|
i _ |
01 |
125 |
000 |
000 |
•По = —1.25 км |
||
|
! |
p |
00 |
250 |
000 |
000 |
sç = 0,25 км |
||
+++ |
00 |
250 |
000 |
000 |
s„ = 0,25 км |
||||
J |
! |
l_ |
02 |
100 |
000 |
000 |
nj = 10 |
||
t |
i |
I |
|
|
|
|
пу = 9 |
||
+'++I 1 1 |
01 |
900 |
. 000 |
000 |
|
|
|||
1 |
1 |
1 |
01 |
100 |
000 |
000 |
пв=і |
||
t + t |
01 |
100 |
000 |
000 |
|||||
01 |
100 |
000 |
000 |
aL — i т/си3 |
|||||
T — 1 — Г |
|||||||||
++ + |
00 |
|
|
|
|
|
|||
+++ |
00 |
|
|
|
|
|
|||
01 |
200 |
000 |
000 |
С = 2 мгл |
|||||
4—1—г |
|||||||||
+++ |
01 |
350 |
• 000 |
000 |
Д = 3,5 км |
||||
+++ |
00 |
|
|
|
|
|
|||
|
|
|
|
|
|
||||
|
|
П р и м е ч а н и е . |
Третья карта |
массива Ml содержит 110 16 3452 0000 |
|||||
0002 К 2 . |
|
|
|
|
|
|
2°. Массив М2 содержит записанные в восьмеричном виде высоты элементарных параллелепипедов. Число перфокарт массива зависит от размеров структуры. Размеры структуры лимитируются объемом оперативной памяти машины: число элементарных параллелепипедов не должно превышать 34651 0 чисел. Массив М2 используется во всех модификациях ОсП9 без исключения. Подготовка М2 включает сле дующие операции:
186