
книги из ГПНТБ / Шарапанов, Н. Н. Методика геофизических исследований при гидрогеологических съемках с целью мелиорации земель
.pdfТ а б л и ц а 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
При использовании описанного способа смещение компенса ционного нуля не сказывается на результате определения постоян ной времени спада. Однако точность расчета т по формуле (36) сильно зависит от ошибок измерения отдельных величин Д£/вп (£)> уменьшаясь с увеличением интервала времени Дt снятия вели чии Д£/Вп (t) с фотоленты. Поэтому для получения значений т с точностью, удовлетворяющей требованиям выполнения работ, необходимо, чтобы амплитуда записи Д£/вп (0 была не менее
100мм.
2.Более трудоемкий, но зато более точный, значительно
меньше зависящий от ошибок в определении отдельных амплитуд Д£/вп (<) — способ расчета т по площади. В этом случае кривую записи спада Д£/вп (£)> как и выше, разбивают на два участка,
заключенных, |
например, |
между tx = 1 с, tm = |
6 с и между t,n = |
|||||
— б е и |
tn = |
11 с; |
для |
каждого участка находят площади S г |
||||
и S |
заключенные между осью абсцисс и соответствующим отрез |
|||||||
ком кривой спада Д£7ВП. |
Среднюю постоянную времени спада т |
|||||||
определяют по формуле |
|
|
|
|||||
|
|
|
|
|
т — |
• |
(37) |
|
Площади S 1 и |
S 2 находят |
любым из известных способов, |
||||||
например, как сумму площадей отдельных трапеций: |
||||||||
|
|
|
Si = |
|
+ |
+ • • • + У т - 1 + |
|
|
|
|
So = |
|
|
■• •+ У п - 1 + - | 5_) |
|
||
где у г, у 2 , ■■., уп |
— амплитуда записи ДЕ/ВП |
для времен спада |
||||||
1,2, |
. . >, /I, a |
At |
~ |
12 |
t ^ — tg |
12 . . • |
|
|
По результатам обработки и вычислений строят кривые ВЭВ |
||||||||
ВП (т1к, |
рк, А к, А*, тк) в зависимости от ABI2 либо соответству |
|||||||
ющие графики электропрофилирования. |
иного параметра |
|||||||
Целесообразность |
использования того или |
в различных гидрогеологических условиях будет рассмотрена ниже.
С е й с м о р а з в е д к а . Первичную обработку сейсмограмм (заполнение паспорта, маркировка лент, проверка правильности отметки момента) выполняют в соответствии с действующей ин струкцией. Аналогично проводят и корреляцию волн. Построение годографов рекомендуется вести в масштабах: по оси рассто яний — от 1 : 200 до 1 : 1000, по оси времен — от 1 см = 0,005 с до 1 см = 0,02 с. При этом более детальные масштабы выбирают при использовании более детальной расстановки сейсмоприемни ков (т. е. меньшей глубины до уровня грунтовых вод) [9].
Р а д и о в о л н о в о е п р о ф и л и р о в а н и е . Основным первичным полевым материалом при работах этим методом яв ляется диаграмма РВП, записанная на бумажной ленте. На диаг-
87
грамме должны быть отражены пометками номер профиля, дата, чувствительность, все пикеты (точки ВЭЗ, ВЭЗ ВП, сейсмозонди рований, скважин, шурфов), расположенные вдоль профиля, границы смен растительности и почв, а также микро- и мезо рельефа (по визуальным наблюдениям геоботаника), контрольные калибровки.
Обработку диаграмм (лент) РВП начинают со сверки подшток с записями в полевых журналах геоботаника. Затем проверяют расхождения в уровнях контрольных калибровочных отметок и проводят линию отсчетного нуля с учетом колебаний в уровне записи по калибровочным меткам. Дальнейшая обработка лент заключается в нанесении на них данпых тарировок по почвогрунтам района и по величинам рк. На этом, по сути дела, первич ная обработка диаграмм РВП заканчивается.
КАЧЕСТВЕННАЯ ИНТЕРПРЕТАЦИЯ РЕЗУЛЬТАТОВ ЭЛЕКТРОРАЗВЕДОЧНЫХ РАБОТ
Остановимся лишь на двух, по нашему мнению, наиболее эффек тивных приемах качественной интерпретации результатов электроразведочных работ, а именно, типизации кривых ВЭЗ ВП и ВЭЗ н построении карт кажущихся сопротивлений и комплекс ного параметра для отдельных разносов АВ.
В различных геолого-гидрогеологических и геоэлектрпческпх условиях тип кривой рк является носителем информации о раз резе [25]. Примером могут служить наши работы в существенно отличных по гидрогеологическим условиям районах — Саратов ском Заволжье и Чуйской впадине.
Анализ кривых рк на площади съемки в Саратовском Заволжье как по их типам (Н, К, A, Q), так и по соотношению кажущихся сопротивлений, а также сопоставление кривых рк с разрезами скважин и результатами их каротажа позволили выделить семь типов кривых, соответствующих площадям распространения раз личных литологофацпальных комплексов пород (рис. 19).
Первый тип кривых рк (KQ) характеризуется достаточно высо кими значениями рк в левой ветви кривой (более 30 Ом-м) и по стоянством рк на малых разносах (до ABI2 = 15 25 м). Такой тип характерен для участков сравнительно неглубокого залегания песков (10—15 м), перекрытых легкими суглинками и супесями. Пески средне- и крупнозернистые, их мощность 40—50 м. Ниже залегает песчано-глинистая толща пород п глины. В геоморфоло гическом отношении рассматриваемый тип кривых соответствует второй надпойменной террасе р. Волга.
Второй тип кривых рк (HKQ) является частным случаем пер вого. Отличительная черта разреза, соответствующего данному типу — наличие в верхней части средних суглинков и меньшая глубина залегания кровли песков (3—4 м).
88
Третий тип кривых (KQ) характеризуется резкой дифферен цированностью по величине рк (14—15 Ом-м в начальной и ко нечной частях кривой, 100—150 Ом-и в средней). Кривые харак терны для участков, где пески выходят на поверхность земли; мощность их составляет 50—70 м; ниже залегает песчано-глини стая толща пород. Пески крупнозернистые, иногда с включениями.
Рис. 19. Карта типов кривых рк.
1—7 — площади распространения I—VII типов кривых р|{; S —
точки ВЭЗ; 0 — опорные скважины; 10 — супеси; П — легкие и
средние суглпнкп; 12 —■пески; 13 — тяжелые суглинки; 14 — глины; 15 — мергели; 1в — известняки
гальки. Рассматриваемый тип кривых соответствует первой над-^
пойменной террасе р. Волга.
Четвертый, пятый и шестой типы кривых рк достаточно близки между собой. Они представлены кривыми вида НК и KHQ (ше стой тип). Кажущееся электрическое сопротивление по кривой меняется от 20—30 до 10—15 Ом-м. Все они характерны для площадей глубокого залегания кровли песчаного горизонта (30—
Sft