Расчет кривых вэз на эвм
Расчет кривых ВЭЗ проводим по программе JPJ, составленной на кафедре геофизики МГУ имени Ломоносова, которая содержит пакет программ для обработки, анализа, интерпретации и визуализации данных ВЭЗ и ВЭЗ-ВП.
Для решения прямой задачи ВЭЗ используется алгоритм линейной фильтрации с различными фильтрами для разных установок зондирования: Шлюмберже (Q, S), Венера (W), дипольной осевой (Д), двухэлектродной АМ и т. д. Выбор установки осуществляется по значению ключа установки (Кust) в третьей строке файла данных. Общая длина кривой зондирования не может превосходить 4 декады, что соответствует разносам питающей линии от 1 м до 10 км, а общее число разносов не может быть более 30.
Порядок расчета кривых ВЭЗ (решение прямых задач ВЭЗ) заключается в нижеследующем (все выкладки описаны согласно расположения файлов программы, установленных на ПК электроразведочной лаборатории):
Активизировать Norton Commander, для этого надо щелкнуть
дважды по вышеуказанному ярлыку.
Вывести обе панели оглавления каталогов диска D: Alt+F1 -
выводит оглавление на левую панель; Alt+F2 – на правую панель.
Найти каталог с именем JVAN VAS и войти в него, нажав ENTER.
Выделить файл dvl.dat, который структурируется следующим
образом:
В первой строке помещают два числа: число разносов, на
которых предполагается вычислить кривую ВЭЗ, KEYAB – ключ режима ввода разносов, от которого зависит способ задания значений АВ/2 во второй строке файла.
Во второй строке задается информация о сетке разносов АВ/2.
Если значения разносов предполагается задавать списком, то KEYAB = 0, и во второй строке следует набрать NR чисел, например: 0.5 1 1.5 2 2.5 3 4 6 9 12 15 20 25 30 40 60 90 (это пример оформления второй строки – 17 чисел, означающих значения АВ/2). Если положить KEYAB = 1, то во второй строке следует задать два числа: RN – начальный разнос AB/2 и КТМ, значение которого рекомендуется выбирать в интервале от 4 до 10. Величина КТМ определяет количество значений ρк на одну декаду (1 – 10, 10 – 100, 100 – 1000 и т.д.). Например, если вторую строку оформить в виде 0.5, а в первой строке задать NR = 14, то кривая ВЭЗ будет рассчитана на 14 разносах, начиная с АВ/2 = 0.5.
В третьей строке задается номер или наименование варианта,
для которого рассчитывается кривая ВЭЗ.
В четвертой строке указывается количество слоев для данного
варианта.
В пятой строке перечисляются значения удельного
электрического сопротивления слоев разреза, начиная с первого.
В шестой строке перечисляются мощности слоев в метрах,
начиная с первого.
Строки 3 – 6 повторяются по количеству вариантов.
Рекомендация. В работе программы Dvall.exe могут возникать прерывания при расчете кривой для последнего варианта. В связи с этим целесообразно последний вариант модели повторять дважды.
Ниже приведем пример исходного файла Dll.dat.
11 0 NN, KEY_AB
4.2 6 9 13 20 30 42 60 90 130 200
1
3
30 10 500
10 4
2
4
30 10 500 30
10 13 30
3
5
50 12 200 10 1000
7 18 50 40
3
5
50 12 200 10 1000
7 18 50 40
5. После набора последней модели необходимо сохранить данные,
нажав F2.
6. Выйти из редактора, нажав F10 или ESC.
7. Отредактировать файл dv_ipi.dat: выделить его, нажать F4 и
указать количество вариантов, заменив первый символ в третьей строке NVEZ на число, соответствующее количеству вариантов. Нажав F2, сохранить отредактированные данные.
8. Скопировать dv_ipi.dat в каталог UTHEB: выделив dv_ipi.dat,
перейти на соседнюю панель (Tab), найдя на ней каталог, войти в него через Enter и вернувшись к dv_ipi.dat, нажать Enter, в результате чего в этом каталоге появится файл с этим же именем.
9. Переименовать этот файл. Выделить его и нажать клавишу F6 – в
открывшемся окне перейти в конец существующего текста (D:/IPIv6/Ivan_vas/DV_QSS), стереть Dv_QSS и набрать новый текст UTHEB/<новое имя файла.dat>.
10. При нажатии Enter вместо dv_ipi.dat в нашем каталоге появится
файл с новым именем.
11. Используя строку панели две точки (. .), перейти в выше стоящий
каталог JPJ, где будет находиться файл ipi.exe.
12. Найдя файл ipi.exe, нажать на Enter и в папке UTHEB найти тот
файл, где помещен Ваш материал. Если в файле отсутствуют ошибки, то перейдя в режим интерпретации, можно проследить все теоретические кривые ВЭЗ как в цифровом, так и в графическом исполнении.
Закончив просмотр всех кривых выйти из программы IPI, нажав клавишу F10, а при сообщении «сохранять или нет» материал, нажать клавишу Enter или Esc в зависимости от необходимости этих кривых для дальнейшей работы или отсутствия ее.
Задание. Построить графическим способом с помощью двухслойных и вспомогательных палеток трехслойные кривые ВЭЗ. Вычислить ρэкв и hэкв эмпирическим способом по формулам 1, 2, 3, 4, определить погрешность определения этих параметров графическим способом. Рассчитать кривые ВЭЗ по заданным параметрам по программе IPI (табл. 3-21). Результаты представить в табличном виде (табл. 1) и в графическом на билогарифмическом бланке. Кривые ВЭЗ, полученные графическим способом и на ЭВМ, построить на одном рисунке разными цветами.
Таблица 1
№ ВЭЗ |
Тип кривой |
ρ1, Омм |
h1, м
|
ρ 2, Омм |
h1, м
|
ρ 3, Омм |
μ1=ρ2/ρ1 |
ν1=h2/h1 |
|
|
|
|
|
|
|
|
|
|
|
Продолжение таблицы 1
|
|
|
|
|
|
|
|
|
|
Ниже дается пример кривой ВЭЗ, рассчитанной по программе JPJ (табл. 2).
Таблица 2
№ ВЭЗ |
Разнос АВ/2 |
||||||||||||||
|
1.5 |
2.2 |
3 |
4.2 |
6 |
9 |
13 |
20 |
30 |
42 |
60 |
90 |
130 |
200 |
300 |
|
ρк, Омм |
||||||||||||||
1 |
150 |
150 |
135 |
110 |
90 |
52 |
27 |
16 |
17 |
23 |
29 |
29 |
25 |
15 |
9 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 3
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
150 |
5 |
20 |
25 |
340 |
2 |
500 |
1 |
60 |
18 |
250 |
3 |
20 |
2.5 |
140 |
15 |
5 |
4 |
8 |
1.5 |
60 |
20 |
10 |
5 |
20 |
3.0 |
240 |
18 |
50 |
6 |
10 |
2.0 |
50 |
12 |
180 |
7 |
24 |
1.5 |
120 |
16 |
820 |
8 |
920 |
2.4 |
210 |
9.0 |
30 |
Таблица 4
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
120 |
1.5 |
12 |
14 |
80 |
2 |
520 |
1.0 |
80 |
25 |
980 |
3 |
15 |
2.4 |
225 |
16 |
25 |
4 |
42 |
3.0 |
560 |
12 |
120 |
5 |
25 |
1.4 |
180 |
21 |
900 |
6 |
8 |
4.2 |
64 |
18 |
420 |
7 |
920 |
1.3 |
460 |
9 |
120 |
8 |
240 |
1.5 |
80 |
14 |
10 |
Таблица 5
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
950 |
3.2 |
150 |
21 |
640 |
2 |
120 |
1.5 |
25 |
15 |
270 |
3 |
12 |
2.7 |
108 |
13 |
23 |
4 |
23 |
1.4 |
241 |
18 |
52 |
5 |
8 |
0.9 |
64 |
11 |
250 |
6 |
237 |
1.9 |
81 |
20 |
340 |
7 |
1500 |
3.4 |
490 |
14 |
92 |
8 |
180 |
0.8 |
60 |
9.5 |
12 |
Таблица 6
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
12 |
3.4 |
102 |
23 |
24 |
2 |
25 |
1.2 |
250 |
14 |
12 |
3 |
12 |
0.9 |
92 |
9.5 |
340 |
4 |
23 |
1.8 |
112 |
16 |
350 |
5 |
240 |
1.5 |
60 |
18 |
120 |
6 |
1200 |
2.8 |
20 |
11 |
104 |
7 |
80 |
1.3 |
17.2 |
7.2 |
5 |
8 |
500 |
1.4 |
180 |
15.4 |
24 |
Таблица 7
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
20 |
1.5 |
120 |
28 |
11 |
2 |
5.0 |
2.6 |
48 |
18 |
10 |
3 |
9.8 |
1.2 |
49 |
11 |
206 |
4 |
20 |
3.0 |
136 |
22 |
1200 |
5 |
150 |
1.0 |
20 |
16 |
820 |
6 |
80 |
1.4 |
15 |
19 |
260 |
7 |
1500 |
2.6 |
130 |
13 |
24 |
8 |
920 |
1.4 |
306 |
22 |
54 |
Таблица 8
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
12 |
1.0 |
108 |
16 |
530 |
2 |
21 |
1.4 |
56 |
20 |
260 |
3 |
230 |
1.8 |
46 |
10 |
10 |
4 |
1120 |
1.5 |
200 |
18 |
34 |
5 |
46 |
1.4 |
10 |
21 |
120 |
6 |
560 |
2.4 |
70 |
17 |
1200 |
7 |
23 |
1.6 |
118 |
24 |
15 |
8 |
43 |
2.0 |
280 |
15 |
24 |
Таблица 9
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
2100 |
1.5 |
1400 |
12 |
250 |
2 |
180 |
1.8 |
86 |
17 |
21 |
3 |
12 |
1.4 |
96 |
21 |
24 |
4 |
30 |
2.1 |
280 |
34 |
96 |
5 |
96 |
1.8 |
24 |
22 |
180 |
6 |
210 |
2.6 |
120 |
19 |
740 |
7 |
15 |
1.6 |
90 |
26 |
360 |
8 |
48 |
3.2 |
420 |
38 |
1200 |
Таблица 10
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
640 |
2.3 |
210 |
13 |
57 |
2 |
1800 |
1.5 |
520 |
26 |
44 |
3 |
25 |
1.8 |
180 |
16 |
920 |
4 |
5.0 |
1.2 |
44 |
9.8 |
211 |
5 |
450 |
3.2 |
180 |
14.6 |
720 |
6 |
120 |
1.3 |
62 |
16.2 |
12 |
7 |
15 |
1.0 |
68 |
21 |
10 |
8 |
60 |
2.4 |
350 |
18 |
20 |
Таблица 11
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
14 |
1.0 |
260 |
12 |
680 |
2 |
34 |
2.1 |
560 |
32 |
1500 |
3 |
124 |
1.8 |
15 |
21 |
160 |
4 |
530 |
3.2 |
92 |
14 |
270 |
5 |
112 |
1.6 |
400 |
36 |
56 |
6 |
36 |
2.7 |
190 |
19 |
15 |
7 |
960 |
1.8 |
280 |
14 |
24 |
8 |
150 |
2.2 |
65 |
9.8 |
12 |
Таблица 12
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
10 |
5 |
40 |
20 |
10 |
2 |
240 |
5 |
1050 |
9 |
70 |
3 |
1000 |
2.5 |
150 |
10 |
2500 |
4 |
100 |
5 |
10 |
10 |
150 |
5 |
1000 |
7.5 |
300 |
30 |
100 |
6 |
80 |
1 |
10 |
4 |
1 |
7 |
65 |
20 |
240 |
30 |
900 |
Таблица 13
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
2 |
15 |
80 |
50 |
10 |
2 |
1 |
1 |
3 |
3 |
100 |
3 |
150 |
1.5 |
60 |
5 |
2500 |
4 |
150 |
10 |
60 |
30 |
10 |
5 |
750 |
10 |
250 |
20 |
50 |
6 |
10 |
5 |
50 |
50 |
30 |
7 |
210 |
12 |
70 |
36 |
15 |
Таблица 14
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
150 |
1 |
900 |
3 |
80 |
2 |
1 |
1 |
3 |
3 |
100 |
3 |
240 |
5 |
1050 |
9 |
70 |
4 |
250 |
3 |
1000 |
5 |
4000 |
5 |
120 |
1.5 |
60 |
7 |
5 |
6 |
1200 |
3.2 |
400 |
18 |
20 |
7 |
40 |
1.5 |
5 |
21 |
120 |
Таблица 15
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
10 |
1 |
10 |
4 |
150 |
2 |
150 |
0.8 |
60 |
7 |
2500 |
3 |
15 |
10 |
50 |
25 |
10 |
4 |
150 |
10 |
600 |
40 |
50 |
5 |
1000 |
7.5 |
300 |
30 |
100 |
6 |
15 |
7 |
45 |
21 |
2 |
7 |
65 |
20 |
240 |
30 |
900 |
Таблица 16
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
1000 |
1.5 |
250 |
12 |
10 |
2 |
250 |
3.0 |
100 |
8 |
1200 |
3 |
15 |
1.0 |
3200 |
11 |
45 |
4 |
5 |
7 |
60 |
9 |
12 |
5 |
520 |
1.5 |
120 |
12 |
480 |
6 |
15 |
0.8 |
60 |
15 |
250 |
7 |
480 |
5.4 |
48 |
18 |
5 |
Таблица 17
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
5 |
12 |
35 |
10 |
2 |
2 |
150 |
10 |
600 |
40 |
1000 |
3 |
15 |
7 |
45 |
21 |
0 |
4 |
15 |
1 |
100 |
15 |
500 |
5 |
1000 |
2.5 |
150 |
10 |
200 |
6 |
80 |
1 |
10 |
4 |
1 |
7 |
210 |
12 |
70 |
36 |
0 |
Таблица 18
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
1200 |
20 |
470 |
80 |
5 |
2 |
5.5 |
2.5 |
20 |
10 |
0 |
3 |
1000 |
10 |
100 |
50 |
3 |
4 |
100 |
15 |
50 |
30 |
10 |
5 |
100 |
1 |
10 |
2 |
150 |
6 |
70 |
3 |
40 |
10 |
200 |
7 |
12 |
110 |
120 |
100 |
5 |
Таблица 19
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
10 |
10 |
40 |
20 |
10 |
2 |
240 |
5 |
1050 |
9 |
70 |
3 |
1 |
1 |
3 |
3 |
0.1 |
4 |
100 |
1 |
10 |
2 |
150 |
5 |
1000 |
7.5 |
300 |
30 |
100 |
6 |
250 |
3 |
1000 |
5 |
4000 |
7 |
55 |
5 |
15 |
10 |
2 |
Таблица 20
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
15 |
10 |
50 |
25 |
10 |
2 |
150 |
1 |
900 |
8 |
80 |
3 |
65 |
20 |
240 |
30 |
900 |
4 |
15 |
5 |
200 |
25 |
750 |
5 |
70 |
3 |
40 |
10 |
200 |
6 |
100 |
5 |
10 |
10 |
150 |
7 |
75 |
10 |
25 |
20 |
5 |
Таблица 21
№ ВЭЗ |
ρ1, Омм |
h1, м |
ρ2, Омм |
h2, м |
ρ 3, Омм |
1 |
200 |
2.5 |
25 |
18 |
180 |
2 |
1200 |
4.8 |
300 |
24 |
800 |
3 |
20 |
5 |
720 |
42 |
50 |
4 |
75 |
2 |
450 |
16 |
40 |
5 |
800 |
6 |
400 |
12 |
150 |
6 |
720 |
1.5 |
240 |
22 |
24 |
7 |
10 |
2.4 |
60 |
14 |
360 |
8 |
26 |
5.8 |
156 |
24 |
780 |