
книги из ГПНТБ / Боревский Б.В. Методика определения параметров водоносных горизонтов по данным откачек
.pdfскважинам |
при откачках из скважин 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 и S—t. |
В |
качестве расчетной |
выбирается |
|||||||
зависимость, |
наиболее |
близкая |
к |
прямолинейной. |
Предложенным |
|||||||
зависимостям |
могут отвечать |
следующие |
виды |
графиков: |
||||||||
|
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 мощности опробуемого гори зонта для исключения влияния несовершенства скважины. Однако,