Файл: Булах Е.Г. Автоматизированная система интерпретации гравитационных аномалий (метод минимизации).pdf

Скачать файл
Заказать решение

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 28.07.2024

Просмотров: 114

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

 

 

 

 

 

 

 

Т а б л и ц а 21

Параметры

 

 

 

 

 

 

 

 

геологиче­

 

1

2

3

4

 

5

б

ской схемы

 

 

 

 

 

 

 

 

 

а

 

0,25

—0,25

0,25

—0,25

 

0,25

- 0,2 5

 

h

 

0,25

0,25

0,5

0,5

 

0,875

0,875

 

H

 

0,5

0,5

0,875

0,875

 

1,375

1,375

 

/,

 

—3,0

—3,0

—3,0

—3,0

—3,0

—3,0

 

к

 

3,0

3,0

3,0

3,0

 

3,0

3,0

^задаіі

 

- 1,7 5

2,0

— 1,75

2,0

—

1,75

2,0

^вычнсл

 

 

0,65

— 1,05

0,44

—1,32

0,68

I

ЦИКЛ

 

—0,93

I I

цикл

 

—0,94

0,81

— 1,03

0,79

—1,29

0,83

Параметры

 

 

 

 

10

 

п

12

геологиче­

 

7

8

9

 

ское! схемы

 

 

 

 

 

 

 

 

 

а

 

0,25

—0,25

0,25

—0,25

 

0,25

—0,25

 

h

 

1,375

1,375

2,0

2,0

 

2,75

2,75.

 

H

 

2,0

2,0

2,75

2,75

 

4,0

4,0

 

/,

 

—3,0

—3,0

—3,0

—3,0

—3,0

—3,0

 

k

 

3,0

3,0

3,0

3,0

 

3,0

3,0

^задан

 

—1,75

2,0

— 1,75

2,0

— 1,75

2,0

^вычнсл

 

— 1,54

1,70

—1,70

2,28

—1,87

2,65

I

ЦИКЛ

 

I I

ЦИКЛ

 

— 1,46

1,65

—1,62

2,24

— 1,77

2,60

 

 

 

 

 

 

Т а б л и ц а

22

 

Координаты

 

 

Значение Д g

 

 

 

 

 

 

 

 

n начальном

подобранное

подобранное

 

У

X

наблюденное

приближении

 

D I цикле

 

во II цикле

 

 

 

 

(I цикл)

 

 

 

 

 

 

 

 

 

 

—8,0

0

 

0,82

0,62

0,76

 

0,78

—6,0

0

 

1,85

1,40

1,73

 

1,76

—4,5

0

 

3,61

3,06

3,46

 

3,53

—3,0

0

 

6,92

8,75

6,88

 

6,96

—2,0

0

10,52

13,26

10,49

 

10,45

- 1 , 5

0

12,90

14,27

12,72

 

12,88

—1,0

0

14,27

14,87

14,26

 

14,33

—0,5

0

14,79

15,20

15,10

 

14,89

 

0

0

14,94

15,30

15,35

 

15,04

 

0,5

0

14,79

15,20

15,14

 

14,87

 

1,0

0

14,27

14,87

14,15

 

14,19

 

2,0

0

10,52

13,26

10,38

 

10,47

 

3,0

0

 

6,92

8,75

7,00

 

7,05

 

4,0

0

 

4,50

4,16

4,37

 

4,43

 

5,5

0

 

2,30

1,79

2,12

 

2,15

 

8,0

0

 

0,82

0,62

0,74

 

0,75

—1,83 и +1,13 км. Чтобы выяснить, как это отразится на резуль­ татах решения задачи, нужно повторить расчет, сделать своеобраз­

ный второй цикл вычислений. При этом в табл.

21 параметры ^

и 12 взяты из результатов вычислений первого

цикла по оси Oy

(табл. 23).

цикле по оси Ох,

Значения параметров, вычисленных во втором

записаны в табл. 21.

 

Теперь могут быть уточнены параметры и в плоскости уОг. Язменив параметры Іг и 1.2, проводим второй цикл расчетов по оси

Параметры

геологиче­ ской схемы

а

h H

'x

и

''задан

^вычнсл I ЦИКЛ 11 ЦИКЛ

Параметры геологиче­ ской схемы

а

h H

' I

^задан

^вычнсл

I

ЦИКЛ

I I

ЦИКЛ

Т а б л и ц а 23

1 2 3 4 5 6

0,25

—0,25

0,25

—0,25

0,25

—0,25

0,25

0,25

0,5

0,5

0,875

0,875

0,5

0,5

0,875

0,875

1,375

1,375

—0,93

—0,93

—1,05

—1,05

—1,32

—1,32

0,65

0,65

0,44

0,44

0,68

0,68

—3,0

3,0

—3,0

3,0

—3,0

3,0

—1,83

1,13

—1,79

1,66

—2,37

2,48

—1,69

1,32

—1,64

1,81

—2,45

2,53

7

8

9

10

u

12

0,25

—0,25

0,25

—0,25

0,25

—0,25

1,375

1,375

2,0

2,0

2,75

2,75

2,0

2,0

2,75

2,75

4,0

4,0

—1,54

—1,54

—1,7

—1,7

—1,87

—1,87

1,70

1,70

2,28

2,28

2,65

2,65

—3,0

3,0

—3,0

3,0

—3,0

3,0

—2,83

2,99

—3,44

3,42

—3,82

3,64

—2,91

3,06

—3,52

3,50

—3,93

3,76

-Oy. Результаты их приведены в табл. 23. При анализе этих данных установлено, что теперь расхождения в искомых параметрах не­ большие (составляют менее 200 м).

Был произведен и третий цикл. Его результаты мало отлича­ ются от второго и в таблицах не приведены.

Результаты вычислений показаны на рис. 15 и 16.

64

Аномалия горизонтального градиента силы тяжести

На рис. 17 показаны аномалии горизонтального градиента силы тяжести для того же геологического объекта. Методика подбора остается прежней.

В табл. 24—27 приводятся данные для вычисления обратных задач по аномалиям ѴХ2 и Ѵуг. За начальное приближение принят параллелепипед, грани которого совпадают с экстремальными точ­ ками. В самом начале, используя приближенные данные о размерах

Номер точки

X

Значение аномалии V XZ

наблюден­ ное

в началь­ ном приб­ лижении (I цикл)

подобран­ ное в 1 цикле

подобран­ ное во 11 цикле

Номер точки

X

 

Т а б л и ц а

24

Значение аномалии

Vxz

наблюден­ ное

в началь­ ном приб­ лижении (I цикл)

подобран­ ное в 1 цикле

подобран­ ное во 11 цикле

ï

- 4 , 0

 

45

18

33

40

9

0,5

—160

— 114

—149

—158

2

—3,0

79

38

67

75

10

1,0

—176

—242

—164

—172

3

—2,5

107

57

9.6

104

11

1.5

—160

—153

— 152

—158

4

—2,0

146

90

137

143

12

2,0

— 134

—90

— 136

—134

5

—1,5

202

153

189

195

13

2,5

—109

—57

— 106

- 1 1 1

6

—1,0

255

242

237

250

14

3,0

- 8 8

—38

—77

—87

7

—0,5

124

114

96

112

15

4,0

—55

—18

—39

—49

8

0

—22

0

—11

—14

16

5,0

—34

—10

—20

—27

тела

вдоль оси 0у, определяем контур в сечении xOz (табл. 24 и 25).

Затем переходим к аномалии

вдоль оси 0у. Размеры тела по про­

стиранию

теперь уже уточнены

предыдущим

вычислением, а в

табл. 25 и 27 приводятся результаты

расчетов. На этом

заканчи­

вается

первый

цикл

вычислений.

 

 

 

 

 

 

Второй

цикл

вновь

начинается с подбора контура

в

плоскости

хОг. Используя

результаты

расчетов вдоль оси 0у, уточняем

разме-

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

25

 

Параметр

 

 

 

 

 

 

 

 

 

 

 

геологиче­

1

2

 

3

 

4

5

 

 

6

 

ской схемы

 

 

 

 

 

 

 

 

 

 

 

 

а

 

1

—1

 

1

 

—1

1

 

— 1

 

 

h

 

0,25

0,25

0,5

 

0,5

0,875

 

0,875

 

 

H

 

0,5

0,5

 

0,875

 

0,875

1,375

 

1,375

 

 

 

- 1,7 5

—1,75

—1,75

 

—1,75

—1,75

 

—1,75

 

 

 

 

1,75

1,75

1,75

 

1,75

1,75

 

1,75

 

^задан

—1,0

1,0

—1,0

 

1,0

- 1 , 0

 

1.0

 

 

 

 

 

 

 

 

 

 

 

 

 

^вычисл

 

 

0,54

—1,10

 

0,98

.-^1,51

 

 

1,77

 

I

ЦИКЛ

—1,0

 

 

 

 

I I

ЦИКЛ

- 0,9 3

0,47

—1,08

 

0,96

тг-1,37

 

 

1,36

Б 2—1443

 

 

 

 

 

 

 

 

 

 

 

65

 

 

 

 

 

Продолжение

табл. 25

Параметр

7

s

9

10

п

12

геилогпче

скоп схемы

 

 

 

 

 

 

 

0

1

—1

1

—1

1

— 1

 

/(

1,375

1,375

2,0

2,0

2,75

2.75

 

H

2,0

2.0

2,75

2,75

4,0

4,0

 

/,

- 1,7 5

— 1.75

- 1,7 5

—1,75

—1,75

— 1,75

 

k

1,75

1,75

1,75

1,75

1,75

1,75

^задан

—1,0

1,0

—1,0

1,0

—1,0

1,0

 

 

 

 

 

 

^вычисл

 

1,96

— 1,68

1,90

—1,67

 

I

ЦИКЛ

— 1.69

1,92

I l

цикл

— 1,63

2,13

—1,73

2,13

- 1 , 7 8

2,20

J Номер точки

1

2

3

4

5

6

7

8

9

10

11

 

Значение

аномалии

 

 

I

 

Ѵуг

 

Номерточки

 

наблюден­ ное

началь­в номпрнближенпи

цикл)(1

CL—

подобран ниево II цикле

 

У

 

 

 

 

 

 

У

 

 

 

 

° 5

 

 

 

 

 

 

 

о ä X

 

 

 

 

і

і

 

= г ГІ

 

 

 

—5,0

51

19

 

26

48

12

0

—4,0

78

37

 

77

86

13

0,5

—3,0

115

80

 

112

118

14

1,0

—2,5

140

127

 

139

148

15

1,25

—2,0

179

219

 

177

177

16

1,5

— 1,75

200

260

 

176

207

17

1,75

- 1 . 5

173

213

 

150

172

18

2,0

— 1,25

104

153

 

114

109

19

2,5

—1,0

67

109

 

83

69

20

3,0

—0,75

40

75

 

57

44

21

5,0

—0,5

25

47

 

35

26

 

 

 

Т а б л и ц а

26

Значение аномалии

 

наблюден­ ное

о началь­ ном приб­ лижении (I цикл)

подобран­ ное в 1 цикле

подобран­ ное во И цикле

0

0

—2

0

—25

—47

—39

—27

—67

— 109

—86

—70

—104 —153 —117 — 111

—173 —213 — 151 — 174

—200

—260

- 1 7 8

—204

—179

- 2 1 9

— 179

—172

—

140

— 127

— 140

— 144

- 1

1 5

—80

—116

—119

—52

— 19

—25

—48

Т а б л и ц а 27

Параметр

 

 

 

 

 

 

геологической

1

2

3

4

5

6

 

схемы

 

 

 

 

 

 

 

О"

1

—1

1

— 1

1

—1

 

h

0,25

0,25

0,5

0,5

0,875

0,875

 

H

0,5

0,5

0,875

0,875

1,375

1,375

 

к

—1

— 1

—1,1

—1,1

—1,5

—1,5

 

к

0,54

0,54

1,0

1,0

1,77

1,77

^задан

—1,75

1,75

—1,75

1,75

—1,75

1,75

 

 

 

 

 

 

''вы чнсл

—1,60

 

 

 

 

2,00

I

цикл

1,50

—1,77

1,79

- 2,0 0

I I

цикл

—1,67

1,66

—1,74

1,72

—2,82

2,96

 

 

J

 

 

 

 

 

66

5*

 

 

 

 

 

Продолжение

табл.

27

Параметр

7

8

ü

10

i l

 

 

12

геологическом

 

 

схемы

 

 

 

 

 

 

 

 

 

а

1

— 1

I

—1

1

 

—1

 

Л

1,375

1,375

2,0

2,0

2,75

 

2,75

 

H

2,0

2,0

2,75

2,75

4,0

 

 

4,0

 

Ii

—1,7

— 1,7

— 1,68

—1,68

—1,67

— 1,67

 

k

1,96

1,96

1,90

1,90

1,92

 

1,92

^задан

—1,75

1,75

- 1 , 7 5

1,75

—1,75

 

1,75

 

 

 

 

 

 

 

 

^вы ч пел

 

 

 

 

 

 

 

 

1

цикл

—2,17

2,09

—2,08

1,97

—2,00

 

1,90

I I

цикл

—2,71

2,69

—3,42

3,39

—3,20

 

3,15

ры тела по простиранию. Результаты

расчетов второго цикла приве­

дены в тех же таблицах.

 

 

 

 

 

 

Вычисленная аномальная функция достаточно хорошо

совпада­

ет с наблюденной, и

нет необходимости повторять

итерационный

процесс.

 

 

 

 

 

 

 

 

Результаты подбора приведены на рис.

17. Здесь

совмещено

два графика, поэтому не приведены значения

аномалий

начального

приближения. Их можно найти по таблицам. На рисунке

приведены

контуры теоретической модели, для которой производились все рас­ четы.

Смотрите также файлы