УДК 556.53:626.34
МЕХАНИЗМ ОБРАЗОВАНИЯ ПРОРАНА ПРИ ГИДРОДИНАМИЧЕСКИХ АВАРИЯХ НА НАМЫВНЫХ ГИДРОТЕХНИЧЕСКИХ СООРУЖЕНИЯХ ПРОМЫШЛЕННЫХ ПРЕДПРИЯТИЙ
© 2011 г. М.Г. Морозов, А.В. Комин, Н.Б. Шакирова
ФГУП «Российский научно-исследовательский институт комплексного использования и охраны водных ресурсов», г. Екатеринбург
Ключевые слова: безопасность гидротехнических сооружений, Mike 11, намывные гидротехнические сооружения, гидродинамическая авария, проран, процесс образования прорана.
В статье показан механизм образования прорана в результате гидродинамической аварии на намывных гидротехнических сооружениях промышленных предприятий, основываясь на расчетах, выполненных авторами при разработке «Расчетов вероятного вреда, который может быть причинен жизни, здоровью физических лиц, имуществу физических и юридических лиц в результате аварии ГТС» по двум объектам: хвостохрани-лище Качканарского ГОКа и шламохранилище № 2 Богословского Алюминиевого завода филиала «БАЗ-СУАЛ». Сопоставляя расчеты образования прорана по этим объектам было установлено, что формирование прорана протекает в 3 стадии, в связи с чем расчет развития прорана производится по определенному алгоритму. А также были выявлены четыре позиции, влияющие на формирование прорана на намывных гидротехнических сооружениях промышленных предприятий, которые необходимо подтвердить в дальнейшем изучении механизма развития гидродинамической аварии с образованием прорана на ГТС намывного типа.
При выполнении комплексов работ по оценке безопасности гидротехнических сооружений (ГТС) промышленности (хвостохранилища, шламонакопители, золоотвалы) как на стадии предпроектных проработок и проектирования, так и на стадии разработки исполнительной документации (декларация безопасности, расчет вероятного вреда, критерии безопасности и т. д.) возникает необходимость определения наиболее вероятных и наиболее опасных сценариев возникновения гидродинамической аварии на сооружениях [1] и, как следствие, параметров волны прорыва и границ зоны возможного затопления.
На сегодняшний день используется множество программных продуктов, с помощью которых можно наиболее точно и быстро определять параметры волны прорыва и границу зоны возможного затопления, на-
Водное хозяйство России
пример, разработанная С.А. Школьниковым программа «SV-1», «БОР» и «River» B.B. Беликова, а также зарубежная — «Mike 11» Датского гидравлического института [2].
Однако использование указанных программных продуктов применительно к намывным гидротехническим сооружениям промышленности имеет ряд существенных ограничений, поскольку механизм образования прорана и дальнейшее развитие гидродинамической аварии на них существенно отличаются от данных параметров на грунтовых плотинах или дамбах наливных накопителей. Указанные отличия возникновения и протекания гидродинамической аварии заключаются в следующем:
1. На наливных накопителях (грунтовых плотинах водохранилищ) при аварии с разрушением сооружений напорного фронта в нижний бьеф изливается только жидкая фаза (вода), в то время как на намывных сооружениях происходит истечение взвеси (пульпы), т. е. наблюдается селевой поток.
2. При расчетах в программных комплексах и использовании существующих методик на начальной стадии расчета задается размер прорана (сразу на всю высоту сооружения или определенное соотношение длины и ширины), в то время как на намывных сооружениях предполагается постепенное по времени развитие прорана и изменение параметров волны прорыва.
B связи с вышеуказанным, для того чтобы оценить достоверность результатов расчетов, полученных с помощью программного продукта, необходимо детально изучить весь цикл гидродинамической аварии, начальной фазой которой является процесс образования прорана и истечения через него неуправляемого потока жидких отходов в нижний бьеф накопителя.
К намывным гидротехническим сооружениям относятся ограждающие дамбы, которые [3] сооружаются намывом из переработанных отходов (хвостов или шламов) после строительства пионерных дамб, возводимых из местных грунтов, в основном из вскрышных пород. Пионерная дамба замывается отходами 1-го яруса намыва таким образом, чтобы крупные частицы отходов откладывались ближе к упорной призме, а мелкие частицы далее к пруду-отстойнику. Наращивание хранилища ограждающими намывными дамбами производится после окончания намыва нижнего яруса и образования соответствующего пляжа — зоны между ограждающими намывными дамбами и прудком-отстойником (рис. 1).
Механизм образования прорана в результате гидродинамической аварии на намывных ГТС промышленных предприятий показан на основании расчетов, выполненных авторами при разработке «Расчетов веро-
Водное хозяйство России № 3, 2011
Водное хозяйство России
Рис. 1. Схема возведения намывных сооружений: 1 — пионерная дамба (первичная дамба обвалования) из вскрышных пород с обратным фильтром; 2 — граница упорной призмы; 3 — намывные дамбы (вторичные дамбы обвалования); 4 — пляж; 5 — прудок-отстойник; 6 — канава; 7 — заскладированные отходы.
ятного вреда, который может быть причинен жизни, здоровью физических лиц, имуществу физических и юридических лиц в результате аварии ГТС» по двум объектам: хвостохранилище Качканарского ГОКа [4] и шламохранилище № 2 Богословского Алюминиевого завода филиала «БАЗ-СУАЛ» [5].
Известно, что при намыве ограждающих дамб накопителей жидких промышленных отходов применяются два вида намыва пульпы: рассредоточенный (через отдельные выпуски) и сосредоточенный (через торец трубы). При этом на хвостохранилище Качканарского ГОКа ограждающие дамбы намывают торцевым способом намыва хвостов, а на шламохранилище № 2 Богословского Алюминиевого завода («БАЗ-СУАЛ») — рассредоточенным намывом шламов.
Согласно указанным документам [4, 5], на хвостохранилище Качканарского ГОКа (рис. 2) сценарием наиболее тяжелой и наиболее вероятной аварии является [4] возникновение суффозионных процессов на участках интенсивного наращивания дамб. При этом образование прорана возможно лишь за счет суффозионного разрушения низового откоса дамбы обвалования, если угол его наклона будет превышать критический (ф/2), что, как правило, имеет место при торцевом намыве отходов.
Что касается шламохранилища № 2 «БАЗ-СУАЛ» (рис. 3), то по результатам расчетов коэффициентов надежности и рисков разрушения ограждающей дамбы [5], было выявлено, что наиболее вероятным сценарием возникновения гидродинамической аварии, является оползень низового откоса ограждающей дамбы, а наиболее опасным сценарием — размыв ограждающей дамбы водным потоком, переливающимся через ее гребень.
Водное хозяйство России № 3, 2011
Водное хозяйство России
На основании сопоставления расчетов образования прорана по вышеуказанным двум объектам, приведенных в документах «Расчетов вероятного вреда, который может быть причинен жизни, здоровью физических лиц, имуществу физических и юридических лиц в результате аварии ГТС», установлено, что формирование прорана (рис. 4) протекает в 3 стадии:
1. Начальная стадия. Образование прорана происходит в результате создания аварийных ситуаций на намывных ГТС промышленных предприятий. Например, на хвостохранилище Качканарского ГОКа [4] возникновение аварийных ситуаций связано с возможностью размыва гребней дамб при разрыве разводящего пульпопровода или сверхрасчетного
Водное хозяйство России № 3, 2011
Водное хозяйство России
Рис. 3. ГТС шламохранилища № 2 Богословского алюминиевого завода.
_. У.В.
-.- - ::-■■"■'" I стадия \_ho- ■ V-, " V II стадия ■..:;., III стадия 1 С;:, ф ;> ан
Рис. 4. Схема формирования прорана в намывной дамбе накопителей: Н0 — напор на подходе к прорану; Н0 — глубина воды в проране; кн — толщина слоя стока в проране на низовом откосе дамбы; ф — угол внутреннего трения; ап — угол наклона низового откоса.
подъема уровня воды в приоткосной зоне, где ведется наращивание высоты дамб торцевым намывом. А на шламохранилище № 2 «БАЗ-СУАЛ» возникновение аварийных ситуаций связано с грубым нарушением Правил безопасности ГТС накопителей жидких промышленных отходов (ПБ 03-438-02) [6], в частности несоблюдением пунктов 8.15 и 8.17, когда на участке производимого намыва не соблюдается технология намыва (практически отсутствует пляж и имеет место обратный уклон,
Водное хозяйство России
превышение гребня ограждающего валика на участке намыва над пляжем составляет 0,2 м и менее) [5].
2. Размыв гребня дамбы и пляжа (горизонтальный участок прорана) на уровне уреза воды [4], сопровождающийся медленным углублением дна прорана из-за небольшой глубины воды в проране, относительно низкой скорости стока и значительной протяженностью прорана, предопределяющей большой объем подлежащих размыву и транспортировке хвостов.
3. После размыва низового откоса дамбы происходит быстрое понижение дна прорана, вследствие чего имеет место лавинообразное увеличение расхода воды через него [4].
Расчет развития прорана
Расчет развития прорана производится по следующему алгоритму:
1. В связи с тем, что дамбы обвалования формируются [4] с помощью бульдозера без уплотнения крупных фракций хвостов (шламов), они практически не препятствуют движению воды через них. При оценке динамики стока воды с участка намыва в сторону низового откоса можно использовать уравнение расхода через водослив с широким порогом, рассчитываемый по формуле, взятой из справочника по гидравлическим расчетам под редакцией П.Г. Киселева [7]:
Qо =
3
О '
(1)
где Q0 — расход пульпы через проран, м3/с; т — коэффициент водослива [7];
Ь0 — начальная ширина прорана, определяется в зависимости от вида сценария аварии на ГТС; в дальнейшем ширина прорана определяется для каждого /-го слоя по формуле [4, 5], м:
Ь, = Ь/ -1 + 2 АН/ • с^ ф,
(2)
где АН, — величина углубки прорана;
ф — угол внутреннего трения грунтов, слагающих борта прорана. Из уравнения (1) вычисляется Н0 (м) — напор на подходе к прорану по формуле:
Но =
Qo
.2/3
тЬод/^?
(3)
Водное хозяйство России № 3, 2011
ссии
2
Водное хозяйств
Глубина воды в проране К0 (м), в свою очередь, рассчитывается по формуле [8]:
2
Ко = 3 но.
(4)
2. Средняя скорость V« (м/с) течения воды по прорану рассчитывается по формуле [8]:
V« =
6о_
Ьо\
(5)
Для определения будет ли иметь место дальнейшая углубка прорана и его развитие в сторону отстойного пруда, по формуле В.Н. Гончарова [9], рассчитывается неразмывающая скорость турбулентного потока vn (м/с):
Vn = Ш
8,8 Н0
Л95 "У
2£(Р, -Ро) ( 3,5ро
50
(6)
где р^ и р0 — плотность минеральной части хвостов (шламов) и воды, соответственно, кг/м3;
Л50 и Л95 — диаметр частиц хвостов (шламов) соответствующей накопленной частости 50 % и 95 %, м.
Если V« > V,,, то дальнейшая углубка прорана и его развитие в сторону отстойного пруда неизбежны.
3. Интенсивность расширения и углубки прорана на любой стадии его развития определяется массой хвостов, которые стекающая по нему в сторону низового откоса дамбы вода будет в состоянии транспортировать. Транспортирующая способность потока рт (кг/м3) рассчитывается по формуле [9]:
Рт = 3(1 + ф)
&
К
95
( 3 \ г \
1 щ -1 |1 - -1
3
1 V 1 VоJ
(7)
где ф — параметр турбулентности.
4. После начала перелива воды через верхнюю бровку низового откоса дамбы (наклонный участок прорана) [4] на участке ее пересечения прораном вследствие резкого увеличения уклона параметры селевого потока существенно изменятся за счет:
— увеличения скорости течения;
— уменьшения толщины стока;
— увеличения транспортирующей способности потока.
Водное хозяйство России № 3, 2011
яйство России
Скорость движения воды по низовому откосу дамбы в общем случае определяется из формулы Шези [10], которая описывает связь средней скорости V (м/с) со средней глубиной и уклоном потока:
= С
'И, 1И,
1И =-
Ь21Ъ •
пИ, 1И,
(8)
где ^ — толщина слоя стока в проране на низовом откосе дамбы, м; ¡И. — уклон низового откоса, %е;
пс — коэффициент шероховатости низового откоса, определялся из справочника [7].
Значение ^ находится из равенства расходов [4] при стоке воды /-го слоя через горизонтальную часть прорана (порог водослива) и его наклонную часть (низовой откос):
а = 0,35 • Ь• И?'2= Ь ■
• VИ =-
(9)
откуда =
г п Л3/5
а • пс
(10)
Ч ' V у
Значения турбулентных неразмывающих скоростей и транспортирующей способности потока вычисляются по формулам (6) и (7), соответственно.
При расчете развития гидродинамической аварии с образованием прорана принимаются следующие характеристики:
— значение расхода пульпы;
— физико-механические свойства и гранулометрический состав грунтов, слагающих пляжную зону и дамбы обвалования;
— объем каждого слоя сработки пруда и площадь зеркала отстойного пруда на отметке, соответствующей /-му слою;
— значения уклона пляжной зоны и заложений откосов дамб обвалования.
Таким образом, для того чтобы результаты расчетов, полученных с помощью программных продуктов, были наиболее достоверны, необходимо учитывать следующие позиции, влияющие на процесс формирования прорана, и которые необходимо подтвердить при дальнейшем изучении механизма развития гидродинамической аварии с образованием прорана на намывных ГТС промышленности:
1. Секционность накопителей жидких отходов промышленных предприятий. Различают накопители односекционные — состоящие из од-
Водное хозяйство России № 3, 2011
ссии
V
И
с
с
Водное хозяйств
ной секции; двух- и более секционные — состоящие из двух и более секций (отсеков).
Оценка динамики разрушения [4] ограждающих дамб односекци-онных накопителей при снижении уровня воды в пруду после начала ее стока через проран описывается дифференциальным уравнением, справедливым при условии, что накопители имеют дамбы кольцевого типа [5] и у них нет дополнительного притока с водосборной площади:
дг = - (г, - хг )3/2
дг ^ (х)
(11)
где г, — отметка зеркала воды в пруду на момент времени г; хг — отметка гребня дамбы до начала его понижения; Р(х) — площадь зеркала воды в пруде при его отметке г; т — коэффициент водослива (прорана); Ъ, — ширина прорана на момент времени г.
У накопителей, состоящих из системы отсеков, в верхнем отсеке при гидродинамической аварии оценка динамики разрушения ограждающих дамб будет описываться дифференциальным уравнением (11), а в нижерасположенных отсеках, когда гидродинамическая авария происходит в результате повышения уровня воды в их отстойных прудах за счет сброса воды из выше расположенных отсеков, т. е. в накопитель будет поступать дополнительный приток воды, оценка динамики разрушения ограждающих дамб будет описываться следующим дифференциальным уравнением [4]:
дг = йо - (г, - гг)
д, ^ (г)
3/2
(12)
где й0 — приток воды в пруд накопителя, м3/с.
Дифференциальное уравнение (11) справедливо для шламохранили-ща № 2 Богословского Алюминиевого завода и Рогалевского отсека хвостохранилища Качканарского ГОКа, а уравнение (12) применяется для Промежуточного и Выйского отсеков хвостохранилища Качканар-ского ГОКа.
2. Начальная ширина прорана и расход зависят от выбранного сценария гидродинамической аварии. Например, при суффозионном размыве начальная ширина прорана приравнивается к участку, на котором скорость развития суффозионного размыва будет наибольшей из-за максимальной интенсивности намыва. При оползневых процессах начальная ширина прорана рассчитывается, исходя из параметров оползня (высо-
Водное хозяйство России № 3, 2011
яйство России
ты и угла наклона откоса на оползневом участке дамбы), а при размыве ограждающей дамбы водным потоком, переливающимся через ее гребень, за начальную ширину прорана принимают протяженность участка, на котором уровень воды совпадает или превышает отметку основания валика.
3. Способ возведения дамбы. При намывном способе ограждающие дамбы обвалования возводятся поочередно из крупнозернистых отходов без уплотнения крупных фракций. Скорость формирования прорана будет значительно выше, чем при наливном способе возведения дамб, при котором ограждающие дамбы строятся сразу на полную высоту из связных грунтов путем их послойного уплотнения.
4. Тип грунтов, слагающих тело дамбы. На второй стадии формирования прорана [4], для несвязных грунтов характерно, что размыв низового откоса происходит с выполаживанием или увеличением его угла наклона до достижения им угла внутреннего трения, после чего смещение при размыве низового откоса происходит с сохранением угла внутреннего трения.
СПИСОК ЛИТЕРАТУРЫ
1. Методика определения размера вреда, который может быть причинен жизни, здоро-
вью физических лиц, имуществу физических и юридических лиц в результате аварии гидротехнического сооружения (РД 03-626-03): утв. Приказом Ростехнадзора от 01.08.2006. < 738.
2. Каганов Г.М., Волков В.И., Секисова И.А. Сравнение параметров волны прорыва, оп-
ределенных различными методами // Гидротехническое строительство. 2010. < 1. С. 31—37.
3. Рекомендации по проектированию и строительству шламонакопителей и хвостохра-
нилищ металлургической промышленности / ВНИИ ВОДГЕО. М.: Стройиздат, 1986. 128 с.
4. Расчет вероятного вреда, который может быть причинен жизни, здоровью физичес-
ких лиц, имуществу физических и юридических лиц в результате аварии ГТС хвос-тохранилища Качканарского ГОКа / ЭЦ УГГГА. Екатеринбург, 2005. 58 с.
5. Расчет вероятного вреда, который может быть причинен жизни, здоровью физичес-
ких лиц, имуществу физических и юридических лиц в результате аварии ГТС шла-мохранилища < 2 Богословского Алюминиевого завода филиала ОАО «Сибирско-Уральская алюминиевая компания» / ЭЦ УГГУ. Екатеринбург, 2010. 66 с.
6. Правила безопасности гидротехнических сооружений накопителей жидких промыш-
ленных отходов (ПБ 03-438-02): утв. Постановлением Госгортехнадзора России от 28.01.2002. < 6 в Минюсте РФ 16.04.02, регистрационный < 3372.
7. Справочник по гидравлическим расчетам / под ред. П.Г. Киселева. Изд. 5-е. М.: Энер-
гия, 1974. 313 с.
8. Методические рекомендации по расчету развития гидродинамических аварий на нако-
пителях жидких промышленных отходов (РД 03-607-03): утв. Постановлением Госгортехнадзора России от 05.06.03 < 51. Введены в действие с 01.08.03.
Водное хозяйство России № 3, 2011
Водное хозяйство России
9. Гончаров В Н. Динамика русловых потоков. Л.: Гидрометеоиздат, 1962. 374 с.
10. Бебенина Т.П. Гидравлика. Техническая гидромеханика. Конспект лекций. Екатеринбург: изд-во УГГУ, 2006. 180 с.
Сведения об авторах:
Морозов Михаил Григорьевич, к. т. н., заведующий сектором гидротехнических сооружений, ФГУП РосНИИВХ, г. Екатеринбург; e-mail: mg_frost@mail.ru
Комин Алексей Владимирович, к. т. н., научный сотрудник сектора гидротехнических сооружений, ФГУП РосНИИВХ, г. Екатеринбург
Шакирова Наталья Борисовна, ведущий инженер сектора гидротехнических сооружений, ФГУП «Российский научно-исследовательский институт комплексного использования и охраны водных ресурсов» (ФГУП РосНИИВХ), 620049, г. Екатеринбург, ул. Мира, 23; e-mail: nb61@yandex.ru
Водное хозяйство России