Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Боревский Б.В. Методика определения параметров водоносных горизонтов по данным откачек

.pdf
Скачиваний:
154
Добавлен:
25.10.2023
Размер:
15.76 Mб
Скачать

скважинам

при откачках из скважин 400, 426 с расходом

148 л/сек

в ноябре,

когда среднемесячный расход реки составлял

130 л/сек,

и из скважин 400, 428 с расходом 237 л/сек в феврале, когда расход реки был равен 370 л/сек. На графике рис. 65, б видно, что фев­ ральские определения систематически выше ноябрьских, хотя при большей величине возмущения они должны были быть меньше.

Следовательно, при уве-

окт,м1кутки

зооо-

25002000' 15001000500 -

0 20 40 60 ВО 100 120 140 160 Q, л/сен

кщ м'/суткй

|

^

300

 

 

 

 

^

^

200200

400

600

000

1000 кт.мУсутки

 

 

Ноябрьские определения (среднемесяч­

 

 

ный

расход

реки

130 л

/сек)

Р и с . 65. Графики зависимости коэффициен­ тов водопроводимости, рассчитанных по вре­ менному прослеживанию понижения, от сте­ пени возмущения (а) и режима реки (б) (усло­ вия см. на рис. 64)

.

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

Таким образом, дан­ ный пример показывает, что определение действи­ тельных параметров пла­ ста при неизвестной вели­ чине дополнительного пи­ тания по временным графикам становится не­ определенным.

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

Поэтому данные проведенной откачки могут

быть использованы лишь при прогнозном расчете понижения при расходе близком к опытному расходу в условиях межени.

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

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

2. В прогнозных расчетах это приводит к занижению расчетных понижений уровня при увеличении проектного расхода по сравнению с опытным.

3. Определение дополнительного питания и расчет параметров с учетом этого питания могут быть выполнены по данным опытных откачек, проведенных при двух разных ступенях возмущения.

4 . У Ч А С Т К И В П Л А С Т А Х СО СЛОЖНОЙ К О Н Ф И Г У Р А Ц И Е Й Р А З Н О Р О Д Н Ы Х Г Р А Н И Ц

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

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

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

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

Таким образом, при обработке результатов опытно-эксплуата­ ционных откачек в рассматриваемых условиях основными вопросами являются:

1)установление эмпирического закона изменения уровня во времени;

2)установление, насколько полно в достигнутой закономерности изменения уровня отражена реальная природная обстановка;

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

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

Вопросы установления закономерности изменения уровня во

времени рассматривались в работах

[48, 106].

В работе Н. М. Елмановой и Л.

С. Язвина [48], посвященной

оценке эксплуатационных запасов термальных радоновых вод Белокурихинского месторождения на Алтае, предлагается для выбора эмпирического закона изменения уровня строить графики в коор­

динатах S—lg

t, SУ

t и St.

В

качестве расчетной

выбирается

зависимость,

наиболее

близкая

к

прямолинейной.

Предложенным

зависимостям

могут отвечать

следующие

виды

графиков:

 

S = A +

C]gt,

S =

A + CVI

S =

A +

Ct.

 

В

частном

случае

может

быть

либо

А = 0

график проходит

через

начало

координат, либо С = 0

стабилизация

уровня.

Такой прием вполне

справедлив, однако он

не

учитывает воз­

можности достижения промежуточных закономерностей. В общем случае если график изменения уровня не аппроксимируется прямой в полулогарифмическом масштабе, должна подбираться степенная зависимость вида S — f (tc) (где С <С 1).

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

подбором

следует сначала последовательно построить

графики

S

--- f (I)

в полулогарифмическом, коренном и линейном масштабах,

а

затем в

случае необходимости строить промежуточные

степенные

закономерности, выбирая показатель степени в зависимости от характера отклонения графика от прямой в координатах S ~\ft. Если график выпуклый вниз, выбирается ряд 1 > С > 0,5, если график выпуклый вверх, С < 0,5.

Этот способ в принципе позволяет решить поставленную задачу. Его неудобство связано лишь с неопределенностью возможного количества вариантов построения графиков *.

После установления вида зависимости S = f (t), аппроксимиру­ ющегося прямой, угловой коэффициент графика может служить своеобразным расчетным гидравлическим параметром.

Наиболее существенной погрешностью в условиях, когда изме­ нение понижения во времени определяется степенной зависимостью,

* Вид уравнения может быть также установлен по кривой изменения уровня способом наименьших квадратов.

является

стремление привести криволинейные графики

S

— l g t

к серии

прямолинейных отрезков с более пологими первыми н кру­

тыми конечными

участками

[106]. При прогнозах это

приводит

к существенному занижению расчетных понижений уровня.

 

Рассмотренные

примеры

Бердичевского, Щучинского

и

Попе­

речного участков иллюстрируют установление характера их гра­ ничных условий на основе установления закономерностей изменения уровня.

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

затель

степени времени. Для этого в работе

[106] предложено исполь­

зовать

построение графиков в координатах

lg S—lg t. В этом случае

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

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

S =

S0+Atc,

(8.67)

где SQ — постоянное во времени

понижение при работе

скважины,

определяемое конфигурацией границ и положением относительно этих границ опытной и расчетных наблюдательных скважин. Это хорошо видно из структуры формул (8.9) и (8.19) для полосообразного и замкнутого кругового пласта.

При построении графиков в координатах l g S-—lg

t игнорируется

величина S0. Действительно, логарифмируя выражение (8.67), имеем:

lg(S-S0)

= lgA + Clgt.

(8.68)

Таким образом, построение билогарифмических графиков следует

проводить

в координатах

l g (S—S0)-—lg

t.

Однако

величина S0 в

рассматриваемых условиях обычно неиз­

вестна, что не позволяет использовать этот вид графиков. О величине 5 0 можно судить лишь косвенно. Если начальное понижение уровня в расчетной наблюдательной скважине мало по сравнению с пони­

жением во времени St,

т. е. выполняется условие S0

<JC St,

вполне

допустимо

построение

графиков в координатах l g S—lg t.

Именно

к

таким случаям относится большинство примеров,

рассмотренных

в

работе

[106].

 

 

 

 

Иногда

<S0 может быть определено приближенным расчетом и уч­

тено при построении билогарифмических графиков. В противных случаях использование изложенного способа неправомерно.

При построении билогарифмических графиков следует также учитывать, что в координатах lg S—lg t ввиду сильного сжатия шкалы понижений криволинейные графики часто имеют видимость прямолинейных. Поэтому выявленная закономерность всегда должна лроверяться при построении соответствующего графика S -f- t°.

Одной из форм выражения гидравлических показателей, ха­ рактеризующих закономерности изменения понижения уровня при его линейной зависимости от логарифма времени, является предло­ женный В . А. Грабовниковым и Б . М. Зильберштейном метод обоб­ щенных параметров [38, 39].

При использовании этого метода реальные неоднородные пласты со сложной конфигурацией границ заменяются для прогноза неко­ торым условным неограниченным однородным пластом с парамет­ рами, определяемыми по конечным участкам графиков S—lg t.

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

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

показатель —f-, как это и сделано в работе [39]. Противопоставление

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

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

При ограниченном дополнительном питании экстраполяция рас­

хода должна

производиться с

учетом изменения

темпа

снижения

уровня при

выявленном законе его

изменения,

как это

показано

в главе 8. Примерами такого

типа

могут служить данные откачек

на Щучинском, Поперечном

и Бердичевском участках.

 

Вопрос о граничных условиях решается на основе достаточно детального изучения гидрогеологических условий участка и анализа развития депрессионной воронки по площади. Если какие-либо границы не нашли отражения в достигнутых закономерностях

изменения уровня, то при использовании

обобщенных параметров

их предлагается учитывать аналитически

[38].

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

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

В целом вопрос о возможности распространения полученных опытных закономерностей изменения уровня на другие условия водоотбора решается приближенно.

Г л а в а 9

ОСОБЕННОСТИ МЕТОДИКИ ОБРАБОТКИ И ИНТЕРПРЕТАЦИИ РЕЗУЛЬТАТОВ О П Ы Т Н Ы Х РАБОТ НА П Р И Б Р Е Ж Н Ы Х УЧАСТКАХ

Как уже указывалось в главе 1, одной из основных задач опытных работ, проводимых при опробовании водоносных горизонтов, свя­ занных с поверхностными водотоками и водоемами, является уста­ новление степени взаимосвязи подземных и поверхностных вод. Степень гидравлической взаимосвязи определяется суммарным со­ противлением русловых отложений, обусловленным несовершен­ ством вреза реки в водоносный горизонт, наличием слабопрони­ цаемого заиленного слоя и неоднородностью русловых отложений. В настоящее время благодаря работам Ф. М. Бочевера [27], В . М. Шестакова [135], В . М. Григорьева [40], М. С. Хантуша [129] и других исследователей получены решения, позволяющие охарактеризовать закономерности движения подземных вод в речных долинах с учетом закольматированности русловых отложений, ширины реки и несо­ вершенства ее вреза при установившемся движении, а также для некоторых случаев неустановившейся фильтрации. Случай отрыва уровня грунтовых вод от русла реки не рассматривается.

Гидравлическое сопротивление, связанное с наличием заиленного (закольматированного) слоя, В . М. Григорьев [40] предложил учи­ тывать с помощью двух параметров: коэффициента сопротивления заиленного слоя 0) и коэффициента перетекания (А). Эти пара­ метры определяются по следующим зависимостям:

А 0 = ^ ;

(9.1)

где К0

и т0

— соответственно, коэффициент фильтрации и мощность

 

km

заиленного слоя;

 

— водопроводимость водоносного горизонта

В

работе

В. М. Григорьева [40] приводятся крайне ориентиро­

вочные значения А 0 , которые он рекомендует принимать для грубых

приближенных расчетов. Эти значения следующие: для горных рек 1—5 суток, для предгорных участков рек 5—10, для равнинных рек 10—15 и для зарегулированных небольших рек 20—30 суток.

Ф. М. Бочевером было получено решение, в котором суммарное сопротивление русловых отложений учитывается путем введения в формулу Форхгеймера дополнительного гидравлического сопро­ тивления:

 

 

 

с _

Q

 

(9.3)

 

 

 

 

2пкт •{Ri + Rt),

 

где

Л 2

— гидравлическое

сопротивление в

условиях совершенной

 

Я2

связи подземных и поверхностных вод;

 

 

— дополнительное

гидравлическое

сопротивление, обусло­

 

 

вленное несовершенством и заиленностью русловых отло­

 

 

жений.

 

 

 

Величина этого дополнительного сопротивления зависит от ко­

эффициента

перетекания

(по В . М. Григорьеву) и ширины

реки.

Как

показал

анализ решения, полученного

Ф. М. Бочевером,

сопро­

тивление русловых отложений оказывает наибольшее влияние на понижение уровня при приближении скважины к реке. Так, если принять сопротивление русловых отложений при расстоянии сква­

жины до реки равном 100 м за 100%, то при уменьшении

расстояния

до 50 и 25 м, сопротивление русловых отложений будет

соответст­

венно составлять 165% и 250%. С приближением к реке значительно увеличивается и относительная роль сопротивления русловых отло­ жений в общей величине сопротивления. Так, при расстоянии до

реки

25

м, сопротивление

русловых отложений

будет

составлять

32%

от

общего, при расстоянии

50 м — 20%

и при расстоянии

100 м —

12%. (Все сопоставления

приведены для широких

рек при

коэффициенте"

перетекания

Ю - 4 1/м, что при

водопроводимости

500 м2 /сутки,

соответствует

коэффициенту сопротивления

А0

20 су­

ток. Значительное влияние на величину сопротивления русловых от­ ложений может оказывать также ширина реки. Как показали ис­

следования

Ф. М. Бочевера,

ширину

реки

можно не

учитывать

только при выполнении следующего

соотношения:

 

 

 

Ь^^ТТ'

 

 

<9-4>

где b — половина

ширины реки.

 

 

 

Для нешироких

рек неучет

ее ширины может занизить сопроти­

вление русловых отложений в 1,5—2,0 раза.

 

Таким

образом,

для учета

сопротивления

русловых

отложений

при оценке эксплуатационных запасов подземных вод, необходимо определить водопроводимость основного водоносного горизонта, коэффициент сопротивления заиленного слоя и ширину реки. Если определение ширины реки и водопроводимости водоносного гори­ зонта во многих случаях больших затруднений не вызывает, то ме­ тодика установления параметра А0 ъ полевых условиях разработана

•очень слабо. Имеющиеся на этот счет

предложения ряда авторов

(Ф. М. Бочевер, Н. Н. Лапшин, Э. М.

Хохлатов [28], Е . Л . Мин-

кин [96]) основаны на предпосылке фильтрационной однородности водоносного горизонта на обоих берегах реки и, кроме того, не учи­ тывают реального строения заиленного слоя. Точность расчетов по формулам, предлагаемым указанными исследователями, сущест­ венно зависит от правильности определения коэффициента .водопро­ водимости.

Поэтому практически в расчетах для условий установившейся фильтрации наибольшее распространение получил метод учета об­ щего сопротивления русловых отложений, предложенный В. М. Шестаковым [135] и основанный на увеличении истинного расстояния от водозабора до реки на величину AL (при этом принимается, что дополнительное сопротивление, связанное с увеличением длины по­ тока на величину AL, будет эквивалентно сопротивлению русловых отложений).

Величина дополнительного сопротивления A L , которая исполь­ зуется при этом методе, связана с параметрами пласта и заиленного •слоя следующей теоретической зависимостью:

AL =

Vkm~A0cth

/1— .

(9.5)

 

 

У ктА0

 

Методика определения

величины

AL по данным опытных

работ

и наблюдений за естественным режимом подземных вод была разра­ ботана В . М. Шестаковым, Ю. О. Зеегофером и Е . Л . Минкиным [135, 53, 96]. Этими же авторами было показано, что при большой

ширине реки этот

метод

по сравнению со строгим

решением

Л . М. Альтшулера

[1]

дает

хорошую сходимость

при

AL < 7L

{ошибка в этом случае

не превышает 20%). Массовые

практические

определения показали, что в реальных гидрогеологических усло­ виях величина AL, определяемая по предлагаемой методике, полу­ чается существенно зависимой от положения опытной и наблюдатель­ ной скважин, а также от метода определения (наблюдения за естественным режимом или опытные работы). Это связано с тем, что теоретические зависимости, принимаемые для расчетов, полностью не отражают особенностей геологического строения речных долин: фильт­ рационной неоднородности водоносных горизонтов, различий в мощ­ ности и проницаемости заиленного слоя вдоль русла, извилистости русла, несоответствия двухслойной схемы строения русловых отло­ жений, принятой при выводе основных зависимостей, реальным природным условиям.

Особое значение из этих факторов, очевидно, имеет фильтрацион­ ная неоднородность водоносного горизонта. В то же время методика •определения AL, предложенная в работах [53, 96], основана на использовании соотношения понижений уровня в точках, располо­ женных между рекой и опытной скважиной, т. е. как раз по тому направлению, в котором особенно сильно проявляется фильтра­ ционная неоднородность.

В

связи

этим в последние годы при оценке запасов на ряде объ*

ектов стал применяться метод определения AL путем решения обрат­

ной

задачи

с использованием

формулы:

 

 

 

 

 

 

 

 

 

с

_

Q _ ы 2(L+AL)-ri

 

 

 

 

 

( 9 6 )

 

 

 

 

 

2пкт

 

 

 

 

 

 

 

При проведении откачки из нескольких скважин с равными де-

битами вместо формулы

(9.6)

используется

(9.7):

 

 

 

 

 

 

In

2(Ь + АЬ)-п

l n V[2

(L+A^-r^

+ r^

,

 

 

2пкт

 

 

 

 

 

 

 

 

 

 

 

.

V[2(L

+ AL)-r^

+

rl_3

+

-J l n ^ [ 2 ( ^

+ A L ) - r 1 1 2 +

r 2 _ T i

,

(9-7)

 

 

 

 

 

 

 

 

 

 

 

 

где

r x — расстояние между

наблюдательной

скважиной,

по

которой

определяется AL,

и

опытной,

расположенной с ней в створе, пер­

пендикулярном реке,

и имеющей номер

1 ;

 

 

 

 

 

г і - 2 ' г і-з» • • •» ri-n

— соответственно

расстояния

между

опыт­

ной скважиной, расположенной в створе с наблюдательной, и опыт­ ной скважиной, имеющей номер 2, 3, . . ., п.

Определение величины AL по формуле (9.7) проводится неслож­ ным подбором.

Величина AL, рассчитанная по формулам (9.6) или (9.7), в отли­ чие от AL, определяемой по теоретической зависимости (9.5), яв­ ляется интегральным параметром, в совокупности учитывающим фильтрационную неоднородность, заиленность русловых отложе­ ний, несовершенства вреза русла.

Для оценки изменчивости этой величины в зависимости от места расположения наблюдательной скважины и количества опытных скважин был проведен численный анализ с использованием строгого решения Ф. М. Бочевера [27]. При расчетах принималось, что на­ блюдательная скважина расположена между центральной скважи­ ной и рекой, по направлению, перпендикулярному руслу, а цен­ тральная скважина расположена в 100 м от уреза реки. Расчеты про­ ведены для значений А0, равных 20 и 80 суткам, водопроводимости водоносного горизонта 500 м2 /сутки и ширине реки 55 м.

Как следует из проведенных расчетов, наибольшее значение величина AL принимает на урезе реки (наименьшее в ближайшей наблюдательной скважине). Расхождение в величинах AL, опреде­ ленных по разным наблюдательным скважинам, незначительное (не превышает 10—20%). Однако эти расчеты были выполнены" для однородного в фильтрационном отношении водоносного горизонта. Поэтому для того чтобы величина AL была действительно гидравли­ ческим параметром, ее следует определять по наблюдательной сква­ жине, расположенной в непосредственной близости от центральной, но на расстоянии, превышающем 0,7 мощности опробуемого гори­ зонта для исключения влияния несовершенства скважины. Однако,

Соседние файлы в папке книги из ГПНТБ