Файл: Автоматизированная система обработки и интерпретации результатов гравиметрических измерений..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

- У

 

 

Л і '

 

 

Д/г. ошибка построения

структурной

карты; 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