Файл: Шарапанов, Н. Н. Методика геофизических исследований при гидрогеологических съемках с целью мелиорации земель.pdf

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

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

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

Добавлен: 21.10.2024

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

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

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

Т а б л и ц а 16

Показатели гранулометрического состава (число пластичности Я, глинистости Г) и величина параметра А * для пород Саратовского Заволжья

пУравнение регрессии

25

In Я =

1,87+

0,96 In Л*

35

In Г, =

3 ,1 +

0,65

In А*

46

In Г 3 =

3,9 — 0,84 In А*

38

1пГ3=

2 ,5 +

0,68

In А*

35

In Г4 =

3 ,1 +

0,6 In А*

45

In Г6 =

9,76 — 2,3

In Л*

 

П р и м е ч а н и е .

Процентное

содержание

 

Досто­

Довери­

 

верный

Г

коэффи­

тельныii

циент

интервал

 

корреля­

для г,

 

ции

% от Г

0,80

0,45

+ 8

0,59

0,39

+ 12

-0,75

0,32

+ 8,5

+0,78

0,38

±10

0,89

0,38

6

-0,44

-0,32

± 9

частиц размером: Г , — 0,005 —

—0,05 мм, Г- —0,05 -2, мм, Г а—менее 0,05, Га —менее 0,01 мм, Г 5 —менее 0,005 мм.

в табл. 16 (по материалам исследований на площадях Волковской оросительной системы в Заволжье).

В результате составления этих зависимостей оказалось воз­ можным по данным интерпретации кривых ВП выделять в районе следующие литологические разности пород (табл. 17).

Т а б л и ц а 17

Интерпретационная таблица для литологического расчленения разреза

 

 

Фракция

 

 

 

 

Порода

глина

пыль

песок

п

л * • 10-*,

Р,

См/м

Ом • м

 

(<0,005 мм),

(0,005-

(0,05-2,0 мм)

 

 

 

 

%

0,05 мм)

 

 

 

 

 

 

 

 

Глина

30—60

Содержится:

17-25

> 4

< 9

 

 

Меньше,

Меньше,

 

 

 

Суглпн-

20-30

чем глины

чем глины

13—17

2 ,5 -4

9 -1 5

Больше,

Меньше,

КИ тя-

 

чем глины

чем пыли

 

 

 

желые

15-20

Меньше,

Больше,

10-13

1,5-2,5

12—22

Суглин-

ки сред-

 

чем песка

чем пыли

 

 

 

ние

10-15

Больше,

Меньше,

7 -10

0,3-1,5

17-35

Суглин-

ки легкие

6 -10

нем песка

чем пыли

4 - 7

0,3-0,9

25-40

Супесь

Меньше,

Больше,

 

 

чем песка

чем пыли

 

 

 

Следовательно, в тех случаях, когда результаты опробования позволяют составить зависимость геофизических параметров от процентного содержания отдельных фракций в породе, точность

80


п однозначность дальнейшей геолого-гидрогеологической интер­ претации по результатам геофизических работ существенно повы­

шается.

В л и я н и е ф и л ь т р а ц и о н н ы х с в о й с т в в о д о ­ н а с ы щ е н н ы х п о р о д и п о р о д з о н ы а э р а ц и и н а г е о ф и з и ч е с к и е п а р а м е т р ы . Изучение зависи­ мости между геофизическими параметрами и фильтрационными свойствами горных пород является важной задачей при проведении съемок с целью мелиорации.

Коэффициент фильтрации как динамическая характеристика пласта не может быть непосредственно определен по какому-либо геофизическому параметру. Однако определив влияние таких показателей пород, как пористость, степень глинистости и удельная поверхность, можно получить уверенные корреля­ ционные зависимости между параметрами ВП и коэффициентом фильтрации пород. Соответствующие корреляционные зависи­ мости, представленные на рис. 18, получены в результате стати­ стической обработки данных параметрических полевых исследо­ ваний в различных районах страны. Водоносный горизонт повсеместно представлен песками, глубина залегания его кровли колеблется в различных районах от 20 до 40 м.

Как следует из графиков на рис. 18, параметр А* уменьшается с увеличением коэффициента фильтрации пород, что хорошо согласуется с результатами лабораторных экспериментов. Зави­ симости в первом приближении располагаются параллельно друг другу; смещение их вверх вызвано уменьшением удельного элек­ трического сопротивления пород.

Результаты статистической обработки этих данных приведены в табл. 18.

Т а б л и ц а 18

Коэффициент фильтрации Л-ф и параметр А * для водонасыщенных пород

 

 

 

 

 

 

 

Достовер­

Довери­

 

 

 

71

 

 

Г

ный

тельный

Район работ

Уравпение регрессии

коэффи­

интервал

 

 

 

 

 

 

 

циент

для г,

 

 

 

 

 

 

 

корреля­

% от Ьф

 

 

 

 

 

 

 

ции

 

Сарпинская

низ­

26

1п£ф =

5,98 — 2,3 In А*

0,72

0,80

±12

менность

-

 

 

 

 

 

 

= 2 ж 6 Ом -м)

20

 

2,99 — 2,3 In Л*

0,84

0,62

± 9

Днепропетровская

1пкф =

область

=

 

 

 

 

 

 

= 20-5-30 Ом-м)

22

In =

1,94 — 2,3 In А*

0,78

0,62

±11

Саратовское

За­

волжье

=

 

 

 

 

 

 

= 20-5-60 Ом -м)

43

In Л-ф =

1,70 — 2,45 In А*

0,84

0,41

± 16

Чуйская

долина

(Р =

 

 

 

 

 

 

 

 

= 25-5-45 Ом -м)

6 Заказ 1073

81


Как и следовало ожидать, наименьший коэффициент корреля­ ции наблюдается при малом удельном сопротивлении пород.

Помимо данных откачек из скважин в некоторых случаях (например, в Чуйской впадине) в качестве опорных использова-

э ° Z * Ч

Рис. 18. Зависимость параметра А* от коэффициента фильтрации Агф:

а — водонасыщеиные породы; 1 р = 2 -f- 6 Ом- м, 2

р = 20 -Ь 30 Ом-м, 3 — р = 20 -f- 60 Ом-м, 4 — р = 25 -г

45 Ом-м

шурфах); б — породы

зоны

аэрации: 1

суглинки (верхний горизонт 0—1 м),

2 — суглинки (го­

р и зо н т

1—5 м), 3 — пески (горизонт

1—5 м)

лись данные откачек из шурфов, когда глубина до воды составляла в среднем 2—3 м.

Зависимость относительной величины комплексного параметра от коэффициента фильтрации пород зоны аэрации в настоящее время изучена недостаточно. На рис. 18, б приведены графики зависимости А* = / (&ф), полученные в результате сопоставления данных количественной интерпретации кривых ВЭЗ ВП с дан­ ными наливов в шурфы для Саратовского Заволжья. Наливы

82

были проведены в интервалах глубин 0,5—1,0 и 2 - 3 м. Разрез здесь представлен супесями, различными суглинками и песками.

Из приведенных графиков следует, что зависимости четко разделяются для суглинков и песков, а также для суглинков, насыщенных и не насыщенных водой. Последнее, по-видимому, объясняется различием в макро- н микроструктуре пород и нали­ чием значительного числа известковых включений в самой верхней части разреза. Внутри выделенных интервалов глубин п литоло­ гических разностей связь достаточно определенная.

 

 

 

 

 

Т а б л и ц а 19

Уравнения множественной корреляции между параметрами ВП

 

 

п водно-физическими свойствами пород

 

 

п

 

Г

а

Закон распреде­

Т и п с в я з и

У р а в н е н и е р е г р е с с и и

ления

*■ II

М= /(р ,

А*)

M= f ( p , D

оо

:

г= ! ( р , а ,

А*)

Г= Пр, А )

Г= 1 ( А ,

А* )

*ф = /(Р.

А , А *)

*Ф = /(Л * .

р)

*Ф= /(Л *,

А)

Ц = 1{А *)

136

In (ЛйГ+1,7) = 4,46 —

0,971

0,0072 •

 

—0,04р—0.066Л *—

 

 

 

—0,01Г + 0,001р2 +

 

 

 

+ 0,001Л *2+

 

 

 

+ 0,00004Г2 —0,951пр +

 

 

136

+ 0,32 In (Л * + 1,5)

0,908

0,021

In (JW+1,7) = 3,76—

 

—0,23р—0,013Л * +

 

 

70

+ 0,004р2

0,904

0,015

In (А/+ 1,7) = 3,63 —

 

—0,14р—0,0137’ +

 

 

136

+ 0,016р2

0,836

0,025

In ( М + 1,7) = 2,94—

201

—0,12р+0,0016р2

0,848

0,020

In J1= —21,95 —0,03р —

 

—14,1Л —0,024Л * +

 

 

 

+ 1,38Л2 + 0,78 In р +

 

 

 

+ 33,1 In (Л f 1,93) +

 

 

201

+ 0,69 In (Л *+ 1,85)

0,800

0,025

1пГ = 2,34+0,04р +

201

+ 2,1Л —0,6Л2

0,780

0,027

1п Г = —6,5 —3.7Л +

 

+ 0.009Л *+

 

 

40

+ 13,01 In (Л+

1,93)

0,801

0,06

In (А ф + 1,98) =

4,09 +

 

+ 0712Л *-0,78Л

 

 

 

-1 ,8 3 In (Л *+1,5) +

 

 

 

+ 2,03 In (Л +

1,8)

 

 

40

0,381пр

0,79

0,06

In (А-ф + 1,98) = 4,63 +

 

+ 0.029Л * —

 

 

 

— 1,24 In (Л * +1,5) —

 

 

40

—0,32 In р

0,76

0,07

1п(А-ф + 1,98) = 3,б9 +

 

+ 0,1Л * —1,32Л —

 

 

40

—1,31 In (Л * +1,5)

0,72

0,08

In (Лф+ 1,98) = 2,93—

 

—0,73 In (Л * + 1,5)

 

 

М -нормальный, р-пормальный, А ‘-нормаль­ ный, Г-нормальный

Т-нормальный, р-нормальный. Л-нормальный, А ‘-нормаль­ ный

/сф-логнор-

мальный, А ‘ -логнор­ мальный,

Л-нормальиый, р-логнормаль- ный

6*

83


Таковы эмпирические корреляционные зависимости, которые используются при геолого-гидрогеологической интерпретации результатов электроразведочных работ. Вместе с тем необходимо отметить следующее. Несмотря на то, что в каждом отдельном случае мы имеем возможность использовать парную корреляцию, взаимосвязь между геофизическими параметрами и литологиче­ ским составом и водно-физическими свойствами пород носит многомерный характер и требует применения множественной корреляции [37, 38]. В табл. 19 приведен ряд уравнений такого рода, характеризующих взаимосвязь электрических параметров с минерализацией грунтовых вод, степенью глинистости пород и с коэффициентами фильтрации водонасыщенных песков. Урав­ нения составлены по данным статистической обработки материалов полевых исследований, выполненных в различных по геолого­ гидрогеологическим условиям районах.

Применение множественной корреляции способствует повы­ шению объективности и надежности изучаемых эмпирических связей, а следовательно, и выполняемой при их помощи геолого­ гидрогеологической интерпретации. Анализ результатов свиде­ тельствует о том, что величина р является крайне устойчивым параметром, характеризующим минерализацию подземных вод п степень засоленности пород зоны аэрации; параметры же А и А* определяются в основном степенью глинистости пород. Данное обстоятельство позволяет применять корреляционные зависи­ мости при обработке материалов геофизических исследований, выполненных в различных районах с песчано-глинистым типом разрезов.

4. ГЕОЛОГО-ГИДРОГЕОЛОГИЧЕСКАЯ ИНТЕРПРЕТАЦИЯ РЕЗУЛЬТАТОВ КОМПЛЕКСНЫХ ИССЛЕДОВАНИЙ

В данном разделе изложены некоторые приемы геолого-гидро­ геологической интерпретации результатов геофизических исследо­ ваний, так как они являются одним из существенных элементов предлагаемой методики работ. При этом рассмотрены приемы как качественной, так и количественной интерпретации с использова­ нием эмпирических корреляционных зависимостей между геофизическими параметрами и водно-физическими свойствами пород.

ПЕРВИЧНАЯ ОБРАБОТКА РЕЗУЛЬТАТОВ НАБЛЮДЕНИЙ

М е т о д в ы з в а н н о й п о л я р и з а ц и и . Основным пер­ вичным (полевым) материалом в методе ВП является осцилло­ грамма (фотолента). На фотоленте должно быть отражено поло­ жение компенсационных нулей каналов ВП и ПР перед началом замера, величина Д£/пр и кривая спада разности потенциалов Д^впДополнительно в процессе производства измерений должна быть получена либо кривая рк (ВЭЗ ВП), либо значение рк на определенном разносе (электропрофилирование). Обработка фото­ лент заключается в вычислении ДС/ап для различных моментов спада потенциалов ВП и уточнении значения Д£/пр. По этим дан­ ным определяют различные параметры вызванной поляризации.

Обработка фотоленты начинается с временной разбивки интер­ вала записи спада Д£/вп. За нулевое время принимают момент выключения поляризующего тока, который определяют по ка­ налу ПР. Начиная с этого момента, на фотоленте отмечаются каждые 0,5 с до конца записи спада Д£/ВпПосле выполнения этой операции приступают к вычислению величин ДС/Вп на Раз_ личные моменты спада. Нулевой линией отсчета при этом является положение компенсационного нуля канала ВП перед записью; произвольная корректировка линии недопустима.

Прежде всего находят коэффициент поляризуемости т]к. Коэф­ фициент поляризуемости вычисляют по формуле (11) на момент времени t = 1 с после выключения тока в цепи АВ. Точность расчета коэффициента поляризуемости определяется точностью

вычисления величин Д£7Вп и Д£/пр.

Применяемая аппаратура позволяет определять величину Д£7пр с ошибкой, не превышающей 5%. Ошибка вычисления ДС/Вп

8 5


зависит от амплитуды ДС/ВП и постоянства компенсационного нуля линии MN. Амплитуда записи Л£/Вп должна быть такой, чтобы ошибки ее графического определения в минимальной мере сказались на результатах вычислений. Если при графическом определении амплитуды записи Af/Bn допускается абсолютная ошибка ±0,5 мм, то амплитуда записи Д£/вп не должна быть меньше 30 мм. В этом случае относительная ошибка графических вычислений не превысит ±2% .

Сползание компенсационного нуля линии M N в процессе измерений — самый существенный источник погрешностей при вычислении АС/ВпДля того чтобы относительная ошибка опре­ деления ДЕ/Вп за счет этого не превысила ±5% , сдвиг нуля по абсолютной величине не должен превышать ±5% измеряемой разности потенциалов вызванной поляризации.

При соблюдении указанных условий измерений возможно вычисление коэффициента поляризуемости с погрешностью не более ±10%, что соответствует требованиям действующей инструкции.

Следующим параметром ВП, подлежащим вычислению,

является комплексный

параметр

А к, который находят

по фор­

муле (13). Ошибка вычисления

А к зависит лишь

от

точности

определения А£/вп (t) при t = 1с, А£7ВП

(t) при t = 11

с и Д£/пр

и при амплитуде записи А[/вп

(/) при

t = 1с более 100 мм не

превышает 5%.

 

для каждого разноса питающих

Используя значения рк и А к,

электродов вычисляют

относительную

величину

комплексного

параметра на каждом разносе /1* = Ак/рк -100 См/м, где вели­ чина А к берется в процентах.

В настоящее время неясен окончательно вопрос о целесо­ образности использования постоянной времени т спада потен­ циала ВП. Но так как этот параметр в отдельных случаях может иметь информационное значение, целесообразно описать наиболее точные способы его определения [35, 45].

1. На кривой спада Д£/Вп (£) выбирают три равноотстоящих

друг

от

друга во времени значения, например

Д£/вп (П) при

11 =

1 с,

Д£7ВП (t2) при t2 = 6 с и Д£/вп (t3)

при t3 = 11 с.

В соответствии с принятым нами экспоненциальным законом спада эти значения должны удовлетворять равенствам

AUBn(t2) = AUBU(t1)e-M/\

АНвп(<з) = АНвп (У е -^/^

где At = t 2 t x = t3 t2 по условию.

Логарифмируя разность этих уравнений, найдем постоянную времени

(36)

£6