Научная статья на тему 'Моделирование процесса нагнетания шахтных вод в водоносные пласты'

Моделирование процесса нагнетания шахтных вод в водоносные пласты Текст научной статьи по специальности «Энергетика и рациональное природопользование»

CC BY
123
31
i Надоели баннеры? Вы всегда можете отключить рекламу.
i Надоели баннеры? Вы всегда можете отключить рекламу.

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

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Моделирование процесса нагнетания шахтных вод в водоносные пласты»

УДК 622.08 В.В. Сенкус

МОДЕЛИРОВАНИЕ ПРОЦЕССА НАГНЕТАНИЯ ШАХТНЫХ ВОД В ВОДОНОСНЫЕ ПЛАСТЫ

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

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

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

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

Н 2 — И2 Q = 1,36к Н2 >СКВ

где Q - дебит скважины, м /с; к - коэффициент фильтрации водоносного горизонта, м/с; Н2 - мощность водоносного горизонта, м; ИСКВ - глубина воды в скважине, м; гСкВ - радиус скважины, м; ЯСКВ - радиус влияния скважины, м.

Радиус влияния скважины ЯСКВ определяется из выражения [1]

Ясш = 3000 Б,4к,

где Б0 = ИСКВ — Н0, м; Но - естественный напор, м; к - коэффициент фильтрации, м/с.

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

т = —-.

Ш

Пористость различных грунтов приведена в табл. 1, а вмещающих пород угольные пласты приведены в табл. 2 [3].

ИсКв = 3000 Б^у/к ; [м]; Я [м], Бо [м].

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

Л = 3000 50 V* ; М; /? [л<], 50 [м]. Рис. 1. Поглощающая скважина

различными авторами предложены эмпирические формулы для коэффициента фильтрации [2]:

m

(1 - 2) т

т

(1 1 ) 2

1

т 3

31 - т

т

формула Конепп k = а1--------

М

формула Цункера k = а2 —

М

d2

формула Терчаги k = а3 —

М

d2

формула Шмиделя к = а —

4 м (1 - т)1Ь '

где а а2, а3, а4 - опытные числа; d -диаметр частиц грунта; м - статическая вязкость;

по Слихерту коэффициент фильтрации следует определять по формуле [2]

^2

к = 10,21-

и

где V - кинематическая вязкость воды; и - скорость фильтрации.

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

и = к],

где ] - гидравлический уклон, т.е. отношение разности давления к расстоянию между рассматриваемыми точками . dH ] = •

Проницаемость грунтов характеризуется коэффициентом кп, которой в общем виде есть функция [1]

кп = д{(т) ,

где d - эффективный диаметр гранул грунта; т - пористость грунта.

V

Грунты Пористость т

Гравий ^ = 2 ^20 мм) 0,3^0,4

Песок ^ = 0,05 +2 мм) 0,3^0,45

Супесь 0,35^0,45

Суглинок 0,35^0,5

Глина 0,4^0,55

Торф 0,6^0,95

Бурый уголь 0,17

Уголь 0,23^0,25

Таблица 2 Пористость пород

Породы Пористость т Примечание

Песчаник почвы пласта Песчаник кровли пласта Сланец песчано-глинистый (кровля) Сланец песчано- глинистый (почва) Известняк (кровля) Алевролиты Аргелиты 0,04-0,16 0,01-0,06 0,04-0,07 0,06-0,12 0,06 0,05 0,09 Восточный Донбасс

Алевролиты карбонатные Песчаники Суглинки 0,13-0,24 0,04-0,14 0,45 Канско-Аческий бассейн

Алевролиты Песчаники мелко- и среднезернистые Суглинки Глины 0,05-0,07 0,05-0,09 0,44 0,13-0,37 Северный Кузбасс

Коэффициент фильтрации к и пористости кп связаны зависимостью [1]

к=,

V

где д - ускорение свободного падения, м/с2; V - кинематическая вязкость воды.

Закон фильтрации Дарси, представленный через коэффициент пористости будет иметь вид

и = 1

Для определения скорости фильтрации необходимо знать коэффициент фильтрации к или коэффициент пористости кп.

В табл. 3 приведены значения коэффициента фильтрации к для различных грунтов [3].

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

Радиус влияния скважины, в первую очередь, зависит от гидравлического напора (линейно) и во вторую очередь, от коэффициента фильтрации нелинейно (степенная функция 1/2).

В табл. 5 приведены примерные величины радиусов влияния скважины.

Грунт к, м/сут.

Глины 0,001

Суглинок тяжелый 0,05

Суглинок легкий 0,05-0,1

Супесь 0,1-0,5

Лесс 0,25-0,5

Песок пылеватый 0,5-1,0

Песок мелкозернистый 1-5

Песок среднезернистый 5-20

Песок крупнозернистый 20-50

Гравий 20-150

Галечник 100-500

Крупный галечник без песочного заполнителя 500-100 и более

Таблица 4

Коэффициенты фильтрации грунтов к [2]

Грунт к, м/сут

Песчаные грунты с примесью глины 5,2-6,1

Песчано-глинистые грунты 2,6-3,5

Грунты проницаемые глинистые 0,86

Таблица 5

Радиусы влияния скважины [1]

Грунт Радиус влияния скважины, м

Мелкозернистые пески 100-200

Среднезернистые пески 250-500

Крупнозернистые пески 500-600

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Мелкий гравий 700-1000

Крупный гравий 1500-3000

где Игор - глубина горизонта от уровня поверхности, м3 уп - плотность горных пород, т/м ; ув - плотности воды т/м3.

Среднее отношение плотности пород к плотности воды составляет 2,6, что в дальнейшем используется в расчетах.

Например. Глубина водоносного горизонта Игор = 200 м; грунт водоносного горизонта - среднезернистые пески с коэффициентом фильтрации к = 5-20 м/сут (0,058 • 10~3 - 0,232 •

10-3 м/с).

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

Дополнительное давление (м. вод.ст.), при котором сохраняется ламинарный режим фильтрации определяем из формулы [1]

Е>

С = ^скв

С " 3000 4к ’

где к - коэффициент фильтрации, м/с.

Полное давление насоса для нагнетания воды в водоносный горизонт представляется выражением с учетом,

что 5о = Инас + Игор

=^Иоор +- Нскв

ув гор 3000 4к

Таблица 6

Скорость фильтрации и м/сут в зависимости от гидравлического уклона ( для грунтов со средним диаметром гранул d мм [2].

Гидравлический

Средний диаметр гранул грунта d, мм

уклон і 1,0 1,8 3 5 10 20 30 40

0,0005 0,5 1,5 3,5 10 30 80 150 250

0,001 1,0 2,9 7 21 58 145 275 450

0,002 1,9 6 14 40 110 275 480 710

0,006 5,2 18 40 113 275 620 930 1240

0,01 8,5 30 67 175 385 830 1220 -

0,02 16 56 125 300 580 1080 - -

0,03 30 100 183 400 860 1800 - -

0,06 54 200 280 660 1060 - - -

0,1 85 320 500 930 1550 - - -

0,2 160 550 970 1530 - - - -

0,3 290 930 1600 - - - - -

0,6 460 1600 - - - - - -

0,8 580 - - - - - - -

Инас1 = 200 • 2,6 + 250

3000 ^0,232 • 10-Инж2 = 200 • 2,6 + 500

553

м. вод. ст.

532

м. вод. ст.

3000 ^0,232 • 10-

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

Дополнительное давление Б0 составляет от 12,3 до 33,0 м. вод. ст., т.е. 0,12-0,33 МПа.

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

В табл. 6 сведены данные, полученные разными авторами для грунтов со средним значением гранул С в диапазоне от 1 до 40 мм; диапазон

і = ■

гидравлического уклона 1 охватывает величины от 0,0005 до 0,8.

Если представить данные в логарифмической шкале (1, и) (рис. 2), то можно выделить линейные зависимости для гранул различной величины, расположенные под углом, уменьшающимся с увеличением диаметров гранул.

Аналитически эту зависимость можно записать уравнением

І9 и - 1д и0 'д 1 - ід 10

= а (6),

где

и = 6 [46 - 0,5) ;

1д 10 = 0,0005; а(6) =

100 - 6 105

При преобразовании получают ~^Ц = а{С),

* [ Г

1д| и: )=1д (г

а( 6)

0 У V 0

Следовательно

и_

и0

а( 6)

Гидравличсский уклон /

(100-С

I--------- ( 1 Ч 105

и = С и С - 0,5) I----------------I

у \ 0,0005 )

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

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

Лордом Релем [4], а затем Стефа-нюком Б. М. [5] при решении капиллярного уравнения получены радиусы

Таблица 7

Дальнодействие капиллярных сил воды [4]

Рис. 2. Скорость фильтрации в зависимости от гидравлического уклона

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

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

Радиус дальнодействия капиллярных сил меняется в диапазоне 6,65-7,68 мм, а высота дальнодействия (- 6,65 - + 6,65 мм).

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

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

Краевой угол смачивания воды и материала в11 Радиус дальнодействия капиллярных сил гд, мм Высота дальнодействия капиллярных сил А/гд, мм

0 6,65 6,65

15 6,85 5,33

25 6,95 4,46

30 7,04 4,07

45 7,30 3,02

60 7,52 2,02

90 7,68 0,00

120 7,52 -2,02

150 7,04 -4,07

180 6,65 -6,65

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

В линейном законе Дарси, первом приближении процесса фильтрации, отмечается, что кажущееся невлияние величины гранул на пористость т (табл. 1) не подтверждается экспериментально, о чем свидетельствует табл. 3, дающая значение коэффициентов фильтрации к. Однако величины этих коэффициентов меняются для схожих грунтов от 2 до 7 раз. Разброс значений, по видимому, вносят силы адгезии, проявляющиеся через краевой угол смачивания водой грунта водоносного горизонта.

С учетом вышесказанного скорость фильтрация равна

, ^ (100-СЧ

и = С (С - 0,5)2 • 1031) 105

или, умножив и разделив на величину гидравлического уклона 1, можно получить

и = 6 (6 - 0,5)(2 • 1031)1

Сравниваем выражение с формулой Дарси и = к I, если считать, что линейный закон Дарси справедлив, то можно получить значение к в следующем виде

к =

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

6 (6 - 5)2 • 103)1

'100-61 , 105

к =

6 (6 - 5)(2 • 103)1

'100-61 , 105

Например, для грунта с величиной гранул 6 = 5 мм (табл. 6), согласно закону Дарси

к = и.

Тогда при і = 0,001 к = ■

21

при і = 0,01 к =

0,001

175

0,01

= 21000 м/сут;

= 17500 м/сут;

930

Следовательно, коэффициент филь-трации к зависит от размеров гранул водопроницаемого грунта и гидравлического уклона.

Так как всегда d >0, то формула будет иметь вид

при 1 = 0,1 к = "01 = 9300 м/сут;

при 1 = 0,2 к = 10>3° = 7650 м/сут.

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

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

При сопоставлении столбцов табл. 6 для гранул с диаметрами 1, 3, 5 и 10 мм, можно заметить, что «инерционность» изменения коэффициента фильтрации к при изменении гидравлического уклона 1 от 0,001 до 0,1 обладает линейной зависимостью. Результат представлен в табл. 8. График изменения кратности уменьшения коэффициента фильтрации представлен на рис. 3.

Выводы

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

Таблица 8

Кратность изменения коэффициента фильтрации к в зависимости от диаметра гранул грунта d

Диаметр гранул грунта d, мм Коэффициент фильтрации при гидравлическом уклоне і, м/сут Кратность изменения коэффициента фильтрации к

0,001 0,1

1 1000 850 1,17

3 7000 5000 1,40

5 21000 9300 2,25

10 58000 15500 3,75

2. Скорость фильтрации пропорциональна коэффициенту фильтрации и гидравлическому уклону, но коэффициент фильтрации имеет сложную нелинейную зависимость от

Рис. 3. Кратность уменьшения коэффициента фильтрации к в зависимости от диаметра гранул грунта d при увеличении гидравлического уклона г от 0,001 до 0,1

величины гранул грунта и гидравли-

ческого уклона.

3. Коэффициент фильтрации заданного грунта уменьшается с увеличением гидравлического уклона.

4. Кратность уменьшения коэффициента фильтрации с увеличением гидравлического уклона пропорциональна диаметру гранул грунта.

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

------------- СПИСОК ЛИТЕРАТУРЫ

1. Штеренлихт, Д.В. Гидравлика. Книга 2 [Текст] / Штеренлихт, Д.В. - М.: Энер-гоатомиздат, 1991. - 367 с.

2. Гидравлика [Текст] / Лейбензон Л.С., Валькер Д. С., Шумилов П.П., Яблонский В. С. - М.- Л.- Новосибирск: Гос. науч.-техн. горно-геолог-нефтяное изд-во, 1934. - 370 с.

3. Справочник (кадастр) физических свойств горных пород [Текст]. - М.: Недра, 1975. - 279 с.

4. Адамсон, А. Физическая химия поверхностей [Текст]. - М.: Мир, 1979. -630 с.

5. Стефанюк, Б. М. Снижение энергозатрат гидравлической технологии добычи угля [Текст]: диссертация на соискание уч. степени доктора техн. наук, глава 6. - Новокузнецк: ВНИИгидроуголь, 1998. - С.215 -263. ГТТТте

— Коротко об авторах--------------------------------------------------------

Сенкус В.В. - Новокузнецкий филиал-институт Кемеровского государственного университета, г. Новокузнецк.

i Надоели баннеры? Вы всегда можете отключить рекламу.