Рис. 31. Схемы фильтрационного течения к береговому сосредоточенному водозабору при естественном потоке, параллельном реке:
а — относительно малый расход водозабора; б — большой расход водозабора Условные обозначения см. на рис. 14
При более интенсивном водоотборе в область питания водозабора попадает и река (см. рис. 31,6). В этом случае дебит водозабора Q складывается из привлекаемых речных вод (их вклад в дебит составляет величину QP) и части естественного бытового потока Qe. Обе эти составляющие дебита водозабора можно опре-делить из выражений пQе = 0,5Q arcctgyp+хp и пQр=0,5Q(2л — — arcctgyp — хр), где Qe = Qe/ (пqxQ); QP=QP/ (пqxQ); xp и ур так же, как и для первой схемы, определяются по формулам (74) и (75).
Рис. 32. График для определения протяженности ДR ЗСО (к схеме на рис. 31):
Из этих формул следует, что при Q>2 xp2~0,5(Q+1), yр2~0,5Q2/(Q+1), при (Q>6 эти приближенные выражения еще более упрощаются: xгр=V0,5Q и yР= — V0,5Q, т. е. водораздельная точка N лежит в этом случае на прямой у= — х.
Протяженность R ЗСО вверх по потоку подземных вод можно определить для обеих схем по такой приближенной формуле:
Формула (76) получена на основе предпосылки о том, что линия тока, вдоль которой движение воды осуществляется с максимальной скоростью, совпадает с осью у. Величина R при этом несколько завышается. Значения величины АR представлены также на графиках на рис. 32.
Протяженность r ЗСО вниз по потоку можно приближенно оценить таким образом:
Tmах=yр. (77)
Ширина ЗСО для этих схем может быть приближенно установлена по формуле (59), a dmax=Qe/2q. Следует иметь в видуг что ЗСО для рассматриваемых схем асимметрична: ее ширина в сторону реки несколько больше, чем в сторону берега.
Минимальное время Tреч движения речных вод к водозабору для схемы на рис. 31,6 оценивается приближенно по формуле (68). При этом величина Tреч будет несколько занижена.
Пример расчета. Для расчета границ второго и третьего поясов ЗСО одиночного водозабора, расположенного у реки, при естественном потоке с уклоном i =0,001 (q = 0,9 м2/сут), направленном вдоль долины, используем те же исходные данные, что и в предыдущих примерах.
Найдем сначала координаты хр и ур водораздельной точки N по формулам (74), (75): Q = 2000/(3,14*0,9*100)=7,08;
В данном случае 7,08>2,01 (6,28 — arcctg 1,75), т. е. неравенство (73) не выполняется и, следовательно, схема фильтрации к водозабору соответствует схеме на рис. 31,6.
Протяженность ЗСО вниз по потоку подземных вод — величину г — приближенно можно оценить по положению водораздельной точки: r=yр=175 м.
Определим теперь протяженность R ЗСО по потоку подземных вод. Найдем сначала расстояние Rq, которое проходит частица воды, двигаясь со скоростью естественного потока: Rв = 0,9*400/(30*0,25) =48 м для Tм = 400 сут; Д,=0,9Х Х104/(30*0,25) = 1200 м для Tх=104 сут.
Далее по графику на рис. 32 находим дополнительное расстояние AR, обусловленное действием водозабора: 9=0,14; Т=0,9*400/(30*0,25-100) =0,48 для Tм = 400 сут; T=0,9*104/(30*0,25*100)=12 для Tх=104 сут.
Следовательно, ДR=1,4, а AR = 1,4*100=140 м (7М = 400 сут) и АЛ=2, а ДR=2*100=200 м (Tx=104 сут).
Далее имеем: R=484+140~190 м и R=1200+200=1400 м, а общая протяженность второго и третьего поясов составит: L= 175+190 = 365 м для Тм и L — = 175+14004-1575 м для Т*.
Найдем теперь составляющие баланса подземных вод, поступающих к водозабоу:
Qp = 2,92-6,28-100-0,9 = 1650 м3/сут.
Время поступления речных вод к водозабору определим по формуле (68): TРеч=2/3*(3,14*30*0,25*1002/2000)=79 сут.
Для расчета ширины ЗСО используем формулу 2d=Qe/q: 2d=345/0,9 = =380 м.
Рис. 33. Схемы фильтрации к линейному береговому водозабору. Естественный поток направлен к реке (а), отсутствует (б), направлен от реки (в), параллелен ей (г). Условные обозначения см. на рис. 14
Линейный береговой водозабор
Другая наиболее распространенная схема предусматривает расположение водозаборных скважин по линии вдоль уреза реки. Протяженность таких линейных береговых водозаборов обычно намного превосходит расстояние от них до реки, вследствие чего на большей части области влияния водозабора фильтрация может рассматриваться как одномерная.
На рис. 33 приведены основные кинематические схемы фильтрации подземных вод к линейным водозаборам (линейным рядам скважин, горизонтальным галереям, дренам и т. п.), работающим вблизи совершенных водотоков и водоемов.
Рис. 33. Продолжение
Естественный поток направлен к реке. В случае, когда естественный поток подземных вод направлен к реке (рис. 33,а), дебит линейного ряда скважин, как и в случае одиночного берегового водозабора, складывается из фильтрующихся речных вод и части естественного потока, т. е. Q=Qp+Qe.
Составляющие QP и Qe дебита водозабора можно найти из следующих выражений:
(78)
и QP=Q — Qe, где l=l/x0 (l — половина длины ряда скважин); УР=УР/Х0,, ур — половина ширины фронта фильтрации подземных вод к водозабору на урезе реки, определяемая по формуле
(79)
Значения составляющей расхода Qe, подсчитанные по формуле (78), графически представлены на рис. 34,а.
Протяженность R ЗСО линейного берегового водозабора вверх по естественному потоку подземных вод может быть найдена из выражения
На рис. 34,6 представлены графики изменения величины AR в зависимости от времени Т и расхода Q.
Расстояние г от водозабора до нижней границы ЗСО определяется по формуле
(81)
Очевидно, что гшах=1. Из уравнения (81), считая, что f=l, можно получить выражение для Треч, в течение которого речная вода достигает водозабора:
(82)
Рис. 34. Графики для определения расхода Qe и протяженности R ЗСО:
а — для Qe (к схеме на рис. 32, а); б — для R (к схеме на рис. 33, а) :
Формулы (81), (82) справедливы для l> (3 — 5)x0. Если длина линейного водозабора меньше указанной, то его можно рассматривать как сосредоточенный и проводить расчеты по описанной выше методике, изложенной выше.
Соотношения (81) и (82) имеют место при Q>2ql, т. е. когда дебит водозабора превышает расход естественного потока на фронте, длина которого равна длине водозабора.
Получить простые аналитические зависимости для определения ширины ЗСО при работе линейного берегового водозабора
Оказалось весьма затруднительно, поэтому для приближенной оценки величины d может быть использована формула (59), полученная на основе балансовых сопоставлений.
В общем случае для расчета ширины ЗСО линейного берегового водозабора можно использовать следующие расчетные зависимости:
d=d0+l; (83)
где d0=d1+ Ф*Дd; d1 — максимальная ширина области захвата для скважины линейного берегового водозабора в случае, если бы она была одиночной, дебит этой скважины Qo=Q/p (величина d1 определяется с использованием графика, приведенного на рис. 26,а); d2 — то же, но для дебита скважины 2Q0; Ф*=(2/п) [(сумма от р=1 до i=1) arcctg(2бj/x0)]; Ad=d2 — d1; p — число водозаборных скважин; 2а — расстояние между эксплуатационными скважинами линейного водозабора. Остальные обозначения показаны на рис. 33,а.
Коэффициент ф* в равенствах (83) физически может быть представлен как показатель взаимовлияния водозаборных скважин. Некоторые его значения в зависимости от расстояния между скважинами 2а и числа скважин приведены в табл. 15.
Полученная величина do может быть оказаться меньше расстояния 2а; это значит, что эксплуатационные скважины настолько далеко располагаются друг от друга, что зона захвата водозабора в целом разбивается на отдельные участки, примыкающие к водозаборным скважинам.
Пример расчета. Хозяйственно-питьевое водоснабжение металлургического комбината и поселка проектируется осуществить водозабором в виде линейного ряда из семи скважин на расстоянии 2а=200 м одна от другой, располагающегося вдоль реки на среднем расстоянии от нее х0=100 м. Общая длина водозаборного ряда 21=200-7= 1400м. Дебит каждой скважины Q0=2500 м3/сут; суммарный водоотбор Q=2500-7= 17 500 м3/сут. Намечаемый к эксплуатации аллювиальный водоносный горизонт представлен песками мощностью т=40 м с коэффициентом фильтрации k=50 м/сут и пористостью л = 0,2.
Гидравлическая связь подземных вод с поверхностными хорошая, уклон естественного потока к реке i= 0,0005 (q=50*40*0,0005=1 м2/сут).
Таблица 15
Значения коэффициента Ф*
2б/x0
|
Ф* при различных значениях р
|
5
|
I ?
|
10
|
1 15
|
1 20
|
30
|
1 50
|
0,25
|
2,64
|
3,44
|
4,33
|
5,37
|
6,11
|
7,15
|
8,47
|
0,5
|
1,87
|
2,32
|
2,79
|
3,33
|
3,70
|
4,24
|
4,89
|
0,75
|
1,44
|
1,74
|
2,06
|
2,42
|
2,67
|
3,02
|
3,46
|
1,0
|
1,16
|
1,39
|
1,63
|
1,90
|
2,08
|
2,35
|
2,67
|
1,5
|
0,82
|
0,98
|
1,Н
|
1,32
|
1,44
|
1,74
|
2,01
|
2,0
|
0,64
|
0,75
|
0,87
|
1,01
|
1,10
|
1,23
|
1,40
|
2,5
|
0,44
|
0,53
|
0,63
|
0,74
|
0,81
|
0,92
|
1,05
|
3,0
|
0,43
|
0,51
|
0,59
|
0,68
|
0,74 j
|
0,83
|
0,94
|
Требуется рассчитать границы третьего пояса ЗСО, чтобы обеспечить сохранение качества откачиваемой воды в течение всего срока эксплуатации (Tx =104сут).
Протяженность ЗСО вверх по потоку подземных вод находим по графику на рис. 34,6.
Параметр Q/12 определяем по формуле Q/l2=17 500*100/(3,14*1*7002) = 1,13, следовательно, безразмерная величина Тз, необходимая для определения АR по графику на рис. 34,6, составит:
По графику на рис. 34,6 находим ДR3 = 0,8, следовательно, ДR=08*700 — — [1,13*100/(l+1,13)1 =560 — 53=507.
При сохранении естественных условий частицы воды за расчетный интервал времени Т* прошли бы расстояние Rq — 1 *104/(40*0,2) = 1250 м, следовательно R=Rq+дR=1250+507=1760 м.
Расстояние до нижней границы ЗСО находим по формуле (81):
так как г>100 м, принимаем г=100 м. Общая протяженность ЗСО, таким образом, составит: 1 = 1760+10=1770 м.
Найдем теперь количество подземных вод, привлекаемых водозабором из естественного потока. По формуле (79) имеем: 2l/Q= (2*3,14*700*1)/17500= =0,25; yР=V(7002 — 1002+2*700*100ctg0,25) = 1010 м; yP = 1010/100=10,l; l= 700/100 = 7.
При полученных значениях ур и 7 по графику на рис. 34,а находим: $е = 0,2; Qe = 0,2-17 500=3500 м3/сут. Таким образом, водозабор будет в основном обес-яечиваться речными водами: Qp=17500 — 3500=14 тыс. м3/сут.
Минимальное время движения речной воды к скважинам водозабора, вычисленное по формуле (82), составит: Tреч= (2*40*0,2*700*100)/(17500 — 2-IX Х700) =70 сут.
Ширину d ЗСО находим по формуле (83), a d1 и d2 определяем по графику на рис. 26,а.
При Q0 = 2500 м3/сут T1=(2500*104)/(3,14*40*0,2*1002) = 100 и q=(3,14-100Х Х1)/2500=0,13, следовательно, d1 = 5 и d1=5*100=500 м.
При Qo=5000 м3/сут T1==200 и q=0,07, следовательно, d2=7 и d2 = 7-100= =700 м.
По табл. 15 находим, что при р=7 2б/х0 = 200 : 100=2 и ф*=0,75. Таким образом, согласно формуле (83), Лд?=700 — 500 = 200 м, d0=500+0,75* 200=650м, d=700-|-650=1350 м.
Максимальная ширина ЗСО по соотношению (59) составит: dmах = 3500/(2Х Х1)=1750 м.
Естественный поток отсутствует. При отсутствии естественного потока схема фильтрации подземных вод к линейному береговому водозабору принимает вид, изображенный на рис. 33,6. В этой схеме область питания водозабора охватывает весь водоносный пласт. Раздельные линии потока подземных вод отсутствуют. При этом, очевидно, дебит водозабора полностью компенсируется притоком речных вод: QP=Q, Qe=0.
Расстояние R до границы ЗСО в направлении берега приближенно может быть найдено по формуле (84) или по графику на рис. 34,0:
(84)
Протяженность ЗСО на участке между рекой и водозабором может быть найдена следующим образом:
r=(QT)l(2mnl); r=пQT/l, (85)
причем rmax=x0.
В соответствии с этим при г=х0 находим наименьшее время движения речной воды к скважинам водозабора:
Tpeч = 2mntx0/Q. (86)
Для приближенной оценки ширины ЗСО линейного берегового водозабора при отсутствии естественного потока подземных вод также может быть использована формула (83). Однако в данном случае величины d1 и d2 должны определяться по графикам на рис. 28, в,г.
Пример расчета. Воспользуемся исходными данными, приведенными в предыдущем примере. При этом принимаем, что естественный поток подземных вод отсутствует (i=0, q=0). Нужно рассчитать границы третьего пояса ЗСО для Tх = 104 сут.
— Протяженность ЗСО вверх по потоку подземных вод, т. е. в сторону берега, можно найти, воспользовавшись графиком на рис. 34,8.
Безразмерные параметры в этом случае следующие: T4=17 500-10V(3,14*40*0,2*7002)=14,2 и 7=700/100=7, следовательно (см. рис. 34,e), R/1=1,25, а R= 1,25-700 = 875 м.
Протяженность ЗСО на участке между рекой и водозабором находим из соотношения (85) r= 17 500*104/(2*40*0,2*700) = 15 625 м.
Полученный результат г>100 м, поэтому принимаем г=100 м. Общая протяженность ЗСО L=100+875 = 975 м.
Время продвижения речных вод к водозабору определяем по формуле (86): Tрсч = 2*40*0,2*700*100/17 500=64 сут.
Для определения ширины ЗСО по формуле (83) используем график, приведенный на рис. 28,0.
При Q0=2500 м3/сут T1=2500*104/(3,14*40*0,2*1002) = 100; Я,=5,3 и d,= = 5,3*100=530 м.
При Q0=5000 м3/сут T1=200; d2=6,7 и d2=6,7*100=670 м.
По табл. 15 при р=7 и 2б/х0=2 ф* = 0,75, следовательно, Ad — 670 — 530= = 140 м; d0=530+0,75*140=640 м и d=700-f640-1340 м.
Естественный поток подземных вод направлен от реки в сторону берега. Фильтрационная схема, соответствующая этому случаю, изображена на рис. 33, в. Область питания линейного берегового водозабора является ограниченной. Ее границами служат река и замкнутая раздельная линия тока с раздельной точкой N, имеющей координаты х=хр и y=0, причем
(87)
На урезе реки ширина области питания водозабора составляет 2yр. Величина ур определяется из следующего соотношения:
Для численного определения ширины ур можно -использовать также график на рис. 34,г. При Q>272 для этих целей удобно
применять приближенную расчетную зависимость
В данной схеме, как и в предыдущей, дебит водозабора полностью компенсируется поступлением речных вод: QP=Q и Qe=0.
Положение границы ЗСО на участке между водозабором и рекой (вверх по потоку подземных вод) устанавливается приближенно с помощью следующего равенства:
(89);
При r=xq из соотношения (89) получается минимальное вре-движения речной воды к скважинам водозабора:
(90)
Расстояние до верхней границы ЗСО, располагающейся в данном случае ниже по потоку подземных вод от водозаборного сооружения, приближенно может быть найдено из выражения
Для определения величины r по формуле (91) можно использовать график на рис. 34,д.
При выводе всех приведенных выше зависимостей по-прежнему использовалась упрощающая предпосылка об относительно большой протяженности водозабора: l> (3 — 5)x0.
Ширина ЗСО в рассматриваемой схеме может быть найдена по формуле (83). В отличие от рассмотренных ранее случаев в данной схеме величины d1 и d2 находятся по графику на рис. 30,0. Максимальная ширина ЗСО, достигаемая при весьма большом (теоретически бесконечном) времени эксплуатации водозаборного сооружения, численно будет равна ур, т. е. dmax = yp.
Пример расчета. Используем исходные данные предыдущих примеров при условии, что естественный поток (i=0,001, q=1 м2/сут) направлен в сторону берега. Требуется определить границы третьего пояса ЗСО для TХ = 104 сут.
Прежде всего найдем координату хр раздельной точки N по формуле (87). Для этого определим безразмерную величину 21/Q: 2l/Q = (2пqt) /Q=6,28 * 1* Х700/17 500=0,25, тогда хр = Х0*р = V1002 — 7002 + 2*100*700 ctg 0,25 = 260 м.
Ширину области питания водозабора на урезе реки определим по графику на оис 34 г где безразмерные параметры, необходимые для определения yf, будут следующие: Q/2l=17500/(6,28*1*700)=3,98 и l=700/100=7. Тогда по графику yP/l=1,4 и, следовательно, yр=1,4-700 = 980 м.
Протяженность ЗСО вверх по потоку (в данном случае в сторону реки) определим по формуле (89):
т. е. R>x0. В этом случае принимаем R=х0=100 м, а время движения речной воды к скважинам водозабора рассчитываем по формуле (90): 7реч = 2-40-0,2Х X 700 * 100/ (17 500+2 * 1 * 700) = 59 -сут.
Для определения протяженности ЗСО вниз по потоку воспользуемся графиком на рис. 34Д Безразмерные параметры q5 и Т$, необходимые для определения г, имеют следующие величины:
по графику л5 = 0,3, а г=0,3-700 — 100=110 м.
Общая протяженность ЗСО составит: L=100-j-l 10=210 м.
Ширину ЗСО найдем по формуле (83). Для определения величин d1 и d2 используем графики, приведенные на рис. 30,в.
При Q0=2500 м3/сут f,=2500-104/(3,14-40-0,2-1002) == 100 и ? = 3,14-100Х Х0,1/2500=0,13, следовательно, d1=2,7 и di=2,7-100 = 270 м.
При Q0=5000 м3/сут T1==200 и q = 0,07, следовательно, 52 = 3,9 и d2=3,9X X 100=390 м.
По табл. 15 при р=7 и 2б/х0=2 ф* = 0,75. Таким образом, Дс?=390 — 270= = 120 м; do=270+120*0,75=360 м и d=700+360 =1060 м.
В данном случае ширина ЗСО получилась близкой к максимальной, определяемой расстоянием yр = 980 м. Некоторое превышение ширины d над величиной dmах=yР объясняется приближенностью расчетных методик, и в частности тем, что величина yр соответствует водозаборной галерее, а величина d в формуле (83) — ряду водозаборных скважин.
Естественный поток подземных вод параллелен реке. Схема движения подземных вод к линейному береговому водозабору в условиях параллельного естественного потока приведена на рис. 33,г.
Область питания водозабора в данном случае ограничена раздельной линией и урезом реки. Раздельная линия проходит через раздельную точку N. Приближенно можно положить, что эта точка лежит ;на прямой x&xq. Тогда можно оцен-ить ординату точки N:
(92)
Параметры ЗСО в рассматриваемой схеме приближенно можно найти из следующих уравнений:
(93) (94)
Величину R по уравнению (94) можно найти также с помощью графика на рис. 34, е.
Общая ширина 2d ЗСО (см. рис. 33,г) в данном случае приближенно (с завышением) может быть определена по формуле (84) при замене в ней R на величину 2d — x0 или по графику на рис. 34,в при той же подстановке.
Пример расчета. Для расчета границ ЗСО линейного водозабора, расположенного у реки, в условиях, когда естественный поток подземных вод направлен вдоль долины и параллелен реке, воспользуемся исходными данными, изложенными в предыдущем примере.
По формуле (92) найдем координату ур раздельной точки N. Для этого сначала определим безразмерный параметр 8l/Q=8*3,14*700*1/17500=1, тогда Ур=700-|-2-100(е — 1)-°-»=850 м.
Протяженность ЗСО вниз по потоку принимаем равной г=yр=850 м, а протяженность вверх по потоку можно найти по графику на рис. 34,е.
Общая протяженность ЗСО 1=850+2376=3220 м.
Ширину ЗСО можно рассчитать по формуле (84) или по графику Ни рис. 34,0: T4=17500*104/(3,14*40*0,2*7002) = 14,2 и 7=700/100=7. По графику на рис. 34,8 находим величину Rfi. Заменяя в выражении для R/l. значение R на 2d — *о, получим: R/l=R/l=(2d — x0)/1=1,2 и 2d= 1,2-700+100=940 м.
Учет фильтрационного несовершенства речного русла при расчетах зон санитарной охраны
Несовершенные речные русла характеризуются значительными потерями напора при фильтрации поверхностных вод в подземный поток. Это может быть обусловлено малой шириной реки и ее слабым заглублением в водоносный пласт, наличием неоднородных включений в подрусловых отложениях, а, также коль-матацией и заиленностью донных отложений.
Практически гидродинамические расчеты береговых водозаборов с учетом несовершенства речных русел обычно проводятся с использованием метода дополнительного слоя. В соответствии с этим методом, понижения уровня подземных вод и структуру фильтрационного потока на берегу расположения водозабора находят по обычным зависимостям, справедливым для гидравлически совершенных рек. Однако для учета несовершенства реки урез ее при проведении расчетов сдвигается на величину.
AL=cth (2bа)/а, (95)
где 2Ь — ширина реки; а — показатель несовершенства реки, определяемый по опытным данным.
При двухслойном строении русловой зоны, когда в дне реки залегает слабопроницаемый экранирующий слой (заиленный, за-кольматированный) мощностью то с коэффициентом фильтрации &о, параметр а выражается следующим образом:
а = YkJ(kmmn), (96)
где km — водопроводимость основного водоносного горизонта.
Метод дополнительного слоя может быть использован также и при оценке понижений уровня подземных вод и исследовании фильтрационного потока на противоположном по отношению к водозабору берегу реки.
При этом действительный урез реки сдвигается от исследуемой точки в сторону водозабора на величину x0+AL*, где
AL*=l/[a sh (2ba)]. (97)
Наиболее распространенными схемами строения эксплуатируемого водоносного горизонта в разрезе являются: однородный и двухслойный пласты.
Однородный пласт. Принципиальная схема течения подземных вод к береговому водозабору вблизи гидравлически несовершенной реки и с естественным потоком, в обоих берегах направленным к реке, представлена на рис. 35.
Рис. 35. Схема водозабора у несовершенной реки:
а — план; б — разрез: 1 — слабопроницаемый слой под рекой; 2 — уровень подземны-х вод при работе водозабора. Остальные условные обозначения см. на рис. 14
Область захвата берегового водозабора у несовершенной реки может включать оба ее берега. При этом расход водозабора Q компенсируется притоком воды из следующих основных источников:
Q=Qe+QP+Q*e, (98)
где Qe — расход естественного потока, привлекаемого водозабором оо стороны берега; Q*e — то же, со стороны противоположного берега; Qp — приток к водозабору речных вод.
Величину расхода Qe в случае, когда естественный поток на участке расположения водозабора направлен к реке, приближенно можно оценить по соответствующим формулам для сосредоточенных или линейных береговых водозаборных сооружений. Для учета фильтрационного несовершенства в соответствии с методом дополнительного слоя во всех приведенных там зависимостях следует заменить XQ на XL:
xL=x0+ДL. (99)
Если бытовой поток q на участке расположения водозабора отсутствует или направлен в сторону коренного берега, то, очевидно, Qe=0.
Составляющие Q*e и Qp можно оценить следующим образом. Сначала рассчитывается величина Q**e — расход подземных вод, поступающих к урезу реки в пределах области питания водозабора со стороны противоположного берега. Величина Q**e находится по тем же зависимостям, что и Qe, однако во всех формулах следует заменить xq на x*L=kL*.
Если поток q* отсутствует или направлен от реки, то Q**e=0.
При Q**e>Q — Qe следует принимать
Q*e=Q-Qe, (100)
a QP=0.
При Q**e3 Qp=Q~Qe~Q**e.
Протяженность и ширину ЗСО на берегу расположения водозабора (параметры R, г и d) следует определять по тем же формулам, что и для сосредоточенных и линейных водозаборов, работающих вблизи совершенных водотоков и водоемов, но при указанной выше подстановке (99) с учетом направления движения потока q.
Следует иметь в виду, что при движении естественного потока к реке или q = 0 rmax=x0, а при естественном потоке, направленном от реки, Rmах=x0 (т. е. в данном случае замена xq на xl не производится).
Время Tреч находится по соответствующим зависимостям для совершенных рек при замене X0 на XL.
Важной особенностью рассматриваемой схемы фильтрации является то, что область захвата водозабора распространяется на противоположный от водозабора берег реки.
Протяженность области захвата и, следовательно, размеры ЗСО на противоположном берегу приближенно можно оценить по скорости естественного потока подземных вод q*. При этом
R*= [q*(T — Треч)/(тп)] — 2b. (101)
Ширина области захвата и ЗСО водозабора на противоположном берегу ориентировочно (с завышением) принимается равной максимальной ширине области захвата на берегу расположения водозабора d.
Расчеты и выделение ЗСО на противоположном по отношению к водозабору берегу реки следует проводить только при наличии здесь опасности загрязнения подземных вод и относительно большой доле расхода Q*e в общем расходе водоотбора Q, что может быть установлено расчетами по формуле смешения.
Рис. 36. Береговой водозабор в двухслойной толще:
1 — пески, супеси, суглинки; 2 — слабопроницаемый слон; 3 — трещиноватый известняк; 4 — водоупор; 5 — естественный уровень подземных вод; 6 — уровень подземных вод при работе водозабора; 7 — фильтр скважины
Двухслойный пласт. Отложения речных долин обычно не выдержаны по вертикали, и довольно часто наблюдается чередование сильно- и слабопроницаемых слоев. Для эрозионных долин, сформировавшихся -в результате размыва коренных пород и последующего отложения речных песчано-гравийно-галечниковых осадков, характерно двухслойное строение водоносной толщи в разрезе: нижний слой представлен цоколем поймы и террас, сложенным коренными породами, обычно более проницаемыми и водообильными; верхний — покровный слой — аллювиальными песками с прослоями глин, суглинков и супесей, в целом менее проницаемыми (рис. 36). Береговые водозаборы эксплуатируют преимущественно нижние горизонты.
Русло реки относительно нижнего эксплуатируемого горизонта может быть как совершенным, так и несовершенным. Э. М. Хох-латов установил, что наличие в области фильтрации контура питания заметно снижает роль покровной толщи как источника питания, которая имеет существенное значение в случае неограниченного по площади распространения пласта. При установившемся движении подземных вод наличие покровной толщи совсем не сказывается на эксплуатации водозабора. Питание нижнего слоя при этом полностью осуществляется за счет поступления речных вод.
Поэтому при расчетах ЗСО в таких условиях следует пользоваться теми же расчетными зависимостями и графиками, которые приведены для сосредеточенных или линейных береговых водозаборов, расположенных вблизи совершенных или несовершенных водотоков или водоемов.
Пример расчета. Используем исходные данные из примера расчета ЗСО для сосредоточенного водозабора. В дополнение к принятым там условиям положим, что река является несовершенной, показатель ее гидравлического несовершенства а = 0,01 м-1, ширина реки 2b = 20 м. Поток на противоположном берегу направлен к реке, причем q*=q=0,9 м2/сут. Требуется определить размеры ЗСО, соответствующие времени Tм = 400 сут.
Найдем величины сдвига уреза реки AL и AL*: по формуле (95) А1 = = cth (20-0,01)/0,01 =510 м и по формуле (97) AL*=1/[0,01 sh (20-0,01)] = = 500 м, следовательно, Хх. = 100+510=610 м и xL*=500 м.
Найдем составляющие расхода подземных вод, поступающих к водозабору. В данном случае Q>nqxL, так как 2000>3,14-610-0,9. Определим расход Qe. Предварительно по формуле (60) найдем параметр ур
тогда
Сопоставляя полученный результат с расчетом для случая совершенной реки, видим, что несовершенство речного русла существенно снижает долю речных вод в питании водозабора за счет дополнительного привлечения естественного потока.
С использованием тех же зависимостей найдем величину расхода Qe**. В этом случае
Принимая во внимание выражение (100), заключаем, что Qe*=2000 — 1950= =50 м3/сут и QP=0.
Таким образом, в данном случае речные воды не будут привлекаться к водозабору и расход естественного потока, поступающего к водозабору с противоположного берега, составит примерно 50 м3/сут.
Найдем теперь параметры ЗСО на берегу расположения водозабора.
Для определения положения верхней границы ЗСО — расстояние R — в соответствии с рекомендациями, данными ранее, используем график на рис 26 б Найдем сначала Rg: Rq=0,9*400/(30-0,25) =48 м.
Значения исходных параметров будут следующие: q=3,14*610*0 9/2000=0 9 и T=0,9*400/(30*0,25*610) =0,1.
По графику находим АR = 0,3, следовательно, ДR = 0,3-610= 180 м a R= = 180+48~230 м.
По графику на рис. 25,0 определим положение границы ЗСО на участке между рекой и водозабором. В данном случае q — 0,9; T1=2000-400/(3 14*30* X0,25*6102)=0,1 и r=0,25, г=0,25*610~150 м. Полученное значение r>rmax= = 100 м, поэтому принимаем r=100 м. Следовательно, L = 230+100=330 м.
Время Греч находим по формуле (62) при соответствующих указанных выше заменах:
Ширину ЗСО находим по графику на рис. 26,а. При этом q = 0,9 и 7 = 0,1, следовательно, d=0,2 и d=0,2*610= 120 м.
Приток подземных вод к водозабору с противоположного берега весьма мал (2 — 3 % дебита водозабора), поэтому расширение ЗСО на противоположный берег в данном случае, по-видимому, необязательно. Однако рассчитаем по формуле (101) возможную протяженность ЗСО и на противоположном берегу:
R*=[0,9*(400 — 140)/(30*0,25)] — 20=10 м.
ГРАФОАНАЛИТИЧЕСКИЙ МЕТОД ПОСТРОЕНИЯ
ЗОН САНИТАРНОЙ ОХРАНЫ
Расчеты ЗСО графоаналитическим способом проводятся с использованием гидродинамической сетки фильтрации подземных вод, которая представляет собой систему линий равного напора (гидроизогипс) и нормальных к ним линий тока.
Построение такой сетки начинается обычно с проведения гидроизогипс грунтового потока в условиях эксплуатации водозаборного сооружения. Для этого определяется :толожение уровня подземных вод в различных точках пласта по всей области фильтрации. Эти точки могут располагаться по некоторой заданной сетке, которая должна сгущаться на участках расположения водозаборных скважин, а также вблизи характерных точек подземного потока (раздельных точек). Уровень подземных вод в фиксированных точках определяется аналитически (в простейших случаях) или с использованием разнообразных моделирующих и цифровых устройств, широко применяемых при гидрогеологических исследованиях. Далее, путем интерполяции известных значений уровня в точках строятся изолинии уровня подземных вод при работе водозабора.
Линии тока проводятся нормально к гидроизогипсам. Для предварительного графического построения линий тока можно использовать следующий прием.
Выбирается какая-либо точка на одной из гидроизогипс с большим значением уровня, из этой точки восставляется перпендикуляр к ней (рис. 37,а). Далее на перпендикуляре находится точка, располагающаяся на середине расстояния между выбранной гидроизогипсой и соседней с ней с меньшим значением уровня; из этой точки опускается перпендикуляр на соседнюю гидроизогипсу, и его пересечение с этой гидроизогипсой дает точку, лежащую примерно на одной линии тока с выбранной вначале точкой.
Полученная таким образом система линий тока после этого уточняется с использованием следующего соотношения:
Qi - Qi+1=еi — еi+1, (102)
где Qi и Qi+1 — расходы подземного потока в соседних ячейках гидродинамической сетки (см. рис. 37,6) в пределах полосы тока; gj и e;+i — расходы воды, поступающие в рассматриваемые ячейки из внешних источников, или расходы воды, извлекаемые из этих ячеек (перетекание, суммарный расход водозаборных скважин в пределах ячейки и т. д.), в последнем случае величины ег и е;+1 должны быть отрицательными.
Из равенства (102) вытекает следующее условие, которое должно выполняться для ячеек гидродинамической сетки, располагаемых в пределах одной полосы тока:
(103)
где li и lг+1 — длины соседних ячеек сетки в пределах полосы тока; v=(km)i/(km)i+1 — отношение величин водопроводимости в этих ячейках; АHi и ДHi+1 — разность отметок гидроизогипс, ограничивающих рассматриваемые ячейки; bi и bi+1 — средняя ширина этих ячеек.
Рис. 37. Схемы к построению гидродинамической сетки:
а — построение линий тока; б — фрагмент гидродинамической сетки; 1 — линии тока; 2 — гидроизогипсы; 3 — ячейки сетки; 4 — полоса тока
Под полосой тока понимается площадь между двумя соседними линиями тока, при этом ячейки представляют собой участки, ограниченные двумя соседними гидроизогипсами и линиями тока.
В случае, когда водопроводимость пласта на исследуемом участке одинакова (v=l), а гидроизогипсы приводятся через равные интервалы (ДHi = ДHi+1) и, кроме того, еi=еi+1 соотношение (103) приобретает вид:
li/li+1=bi/bi+1. (104)
Корректировка линий тока проводится таким образом, чтобы соотношения (103) или (104) выполнялись в пределах каждой полосы тока во всей рассматриваемой области.
Анализ гидродинамической сетки начинается с отыскания особых (раздельных) точек фильтрационного потока подземных вод к водозабору. Эти точки характеризуются тем, что в них сходятся несколько линий тока и изолиний напора, а скорость движения воды равна нулю. Раздельные линии, т. е. линии тока, проходящие через раздельные точки, ограничивают участки фильтрации подземных вод к водозабору из различных источников питания (рис. 38). Соответствующие фильтрационные расходы из этих источников определяются по формулам
(105)
где Qi — расход подземных вод из данного источника в пределах i-й полосы; i=1, 2, ..., n; n — число полос тока между раздельными линиями, ограничивающими область фильтрации из данного источника; (km)i,j, ДHi,j, bi,j и li,i — параметры j-й ячейки в пределах i-й полосы; еСум i,j=еi,1+еi,2+ ... +еi,j — суммарный расход воды из дополнительных источников питания (перетекание, инфильтрация, скважины), причем суммирование проводится по ячейкам, располагающимся выше по потоку от рассматриваемой i-й ячейки.
Средняя действительная скорость v*i,,- в пределах j-й ячейки i-й полосы тока находится по зависимости
(106)
где ki,j — коэффициент фильтрации в пределах ячейки; ni,j — пористость.
Время, в течение которого частицы воды пройдут путь длиной li,j, можно найти по формуле
(107)
Полное время движения подземных вод между двумя точками, расположенными в ячейках с номерами р и т в полосе тока, получаем из выражения
(108)
Задаваясь расчетным интервалом времени Т, по формуле (108) можно найти расстояние до границ ЗСО в пределах каждой полосы тока, а по соотношениям (105) оценить расходы основных источников, питающих водозабор, и, следовательно, результирующую концентрацию Св тех или иных компонентов в воде, отбираемой водозабором (по формулам смешения). Примеры зон санитарной охраны, рассчитанных графоаналитическим методом с использованием карты гидроизогипс, построенной по данным полевых наблюдений в условиях действующего водозабора и численного моделирования фильтрации на ЭВМ, описаны в гл. 11. В моделях отражены фильтрационная неоднородность эксплуатируемого водоносного горизонта, гидравлическая связь реки с аллювиальным и эксплуатируемым водоносными горизонтами, инфильтрационное питание и другие природные факторы, одновременный учет которых в аналитических расчетах фильтрации и миграции практически невозможен. Рассчитанные контуры второго и третьего поясов ЗСО имеют довольно сложные очертания, а размеры этих зон значительны, что определяется большой скоростью движения подземных вод в трещиноватых породах с высокой водопроницаемостью и малой активной пористостью.
Рис. 38. Схема к определению расхода подземных вод по гидродинамической сетке:
1 — водозабор; 2 — раздельная линия; 3 — линии тока
Достарыңызбен бөлісу: |