CC BY
104
ПРОБЛЕМЫ ЭКОЛОГИИ И РАЦИОНАЛЬНОГО ПРИРОДОПОЛЬЗОВАНИЯ
УДК 630.378, 539.3
А. Г. Поздеев, Ю. А. Кузнецова, Ю. В. Лоскутов, А. Е. Трухан
АНАЛИЗ ПОСЛЕДСТВИЙ РАЗРУШЕНИЯ ПОДВОДНЫХ ПЕРЕХОДОВ НЕФТЕПРОВОДОВ И ИХ ВЛИЯНИЯ НА ЭКОЛОГИЧЕСКОЕ СОСТОЯНИЕ ВОДОТОКОВ
Рассмотрен процесс разбавления загрязненных вод и составлена модель кислородного баланса водоемов при загрязнении нефтяными выбросами на основе системного подхода. Произведена оценка распространения нефтяного загрязнения в водном потоке. Разработана системная модель оценки ущерба окружающей природной среде от загрязнения нефтепродуктами для сбора нефти с поверхности реки при аварийных разливах. Предложена конструкция нефтесборщика.
Ключевые слова: экологическая безопасность, системная модель, нефтепроводы, подводные переходы, нефтесборщик.
Введение. Эксплуатация магистральных нефтепроводов является одним из наиболее опасных источников нарушения экологического равновесия природной среды. В процессе эксплуатации линейные части нефтепроводов могут разрушаться, приводя к загрязнению окружающей природной среды. Ликвидация последствий разрушения нефтепроводов требует больших экономических затрат и приводит к необратимым экологическим последствиям.
Анализ причин возникновения чрезвычайных ситуаций при транспорте нефти и нефтепродуктов [1] показал, что особую опасность представляют участки через водные преграды: реки, ручьи и болота. Нарушение гидрологического режима водотоков, связанное с прокладкой в их створах подводных переходов, вызывает появление силовых нагрузок, не учитываемых при проектировании линейных частей трубопроводов [2, 3]. Это существенно снижает надежность и экологическую безопасность нефтепроводов в целом. Поэтому разработка мероприятий для обеспечения бесперебойной работы подводных переходов нефтепроводов и ликвидации последствий нарушений их герметичности с целью соблюдения экологического равновесия водотоков является актуальной задачей.
© Поздеев А. Г., Кузнецова Ю. А., Лоскутов Ю. В., Трухан А. Е., 2011.
В имеющихся литературных источниках содержатся разрозненные данные, имеющие отношение к анализу последствий разрушения подводных переходов нефтепроводов и их влияния на экологическое состояние водотоков. В частности, ряд работ посвящен определению причин разрыва стенок трубопроводов и оценке объемов залповых выбросов нефтепродуктов [2]. В ряде работ делаются попытки количественной оценки содержания нефтепродуктов в речной воде при аварийных ситуациях [4]. Наиболее полная оценка влияния разбавления загрязнений водами природных объектов выполнена в источниках, посвященных расчету разбавления сточных вод промышленных предприятий [5,6]. В отдельных работах предлагаются инженерные методы прогнозирования качества воды водоемов при загрязнении сточными водами [6], при этом в работах, относящихся к защите водоемов от нефтяных загрязнений, описываются технические устройства преимущественно мобильного типа, не привязанные к конкретным гидрологическим параметрам водоемов [7]. По нашему мнению, исследования процессов распространения загрязнений в водных потоках должны осуществляться на основе системного подхода [8], при котором параметры оцениваются во взаимосвязи друг с другом с учетом их изменения во времени [9]. В связи с этим должны создаваться реализуемые в прикладной программной среде математические модели [10], которые основываются на применении эффективных численных методов решения дифференциальных уравнений [11]. При этом численные значения исходных для расчета параметров должны соответствовать нормативным и экспериментальным данным [1], а также принимать во внимание правила охраны окружающей природной среды при транспорте нефти [3].
Целью работы является получение количественных оценок экологических последствий разрушения подводных переходов при эксплуатации магистральных нефтепроводов на основе системного анализа.
Объектом исследования является распространение нефтяного загрязнения при повреждении подводного перехода. Предметом исследования являются модели, алгоритмы и программы решения задач распространения нефтяного загрязнения при разрушении подводных переходов нефтепроводов и их влияния на экологическое состояние водотоков.
Решаемые задачи:
1) проанализировать характер распространения нефтяного загрязнения в водном потоке при разрушении подводных переходов нефтепроводов;
2) произвести решение обыкновенного дифференциального уравнения одномерной турбулентной диффузии загрязнения методами информационных технологий для оценки эффективности применения прикладных программных пакетов по предмету исследования;
3) выявить методы численного дифференцирования уравнений нестационарной турбулентной диффузии, свободные от неустойчивости решений;
4) разработать балансовую системную модель распространения загрязнения в водном потоке и на основании известных параметров залпового выброса нефти и гидрологических характеристик водотока произвести расчет зависимостей дефицита кислорода от времени и расстояния от места выброса нефти до створа максимального дефицита кислорода;
5) разработать системную модель ущерба окружающей природной среде от загрязнения нефтепродуктами и на основе программы расчета ее показателей в имитационном режиме оценить динамику массы углеводородов;
6) наметить способы сбора нефтепродуктов с поверхности водотоков с учетом анализа характера распространения нефтяных загрязнений при аварийных разливах.
Методы исследования. Для решения сформулированных задач исследования использованы методы математического и системного анализа, элементы теории дифференциальных уравнений и их численного решения, а также приемы программирования в прикладной программной среде MathCad.
Характер распространения нефтяного загрязнения в водном потоке. Практика показывает [1], что основным негативным последствием аварий на магистральных нефтепроводах является воздействие на окружающую природную среду загрязнения водных объектов, почвы и атмосферы.
Нефтяное пятно при залповом выбросе из подводного трубопровода представляет собой полосу, в средней области более широкую, чем по краям [4]. Из трещин нефть уходит каплями и на поверхности образует отдельные пятна [2]. У берегов нефтяное пятно прижимается к линии уреза.
На реках со слабым течением скорость движения нефтяного пятна составляет 3...4 % скорости ветра. Полное перемешивание нефтяных продуктов с пресной водой происходит за 4.. .7 суток. Изменение концентрации нефтяного загрязнения в водотоке связано с процессом перемешивания и разбавления нефтяного пятна водами основного потока, с одной стороны, а также разложением нефти в результате биохимических процессов, с другой стороны.
Модель процесса разбавления загрязненных вод в проточных водных объектах. Залповый выброс нефти приводит к загрязнению водоема, которое распространяется под действием течения вниз по потоку. В процессе движения нефтяное загрязнение перемешивается с водой.
В случае одномерного приближения распространение загрязнения вдоль потока под влиянием турбулентного перемешивания описывается уравнением турбулентной диффузии [11]:
дС-у, дС- = 0, (1)
дt ' д' ' д'2
где С - концентрация загрязнения, кг/м3; V - скорость водного потока, м/с; А - коэффициент турбулентной диффузии, м2/с; ' - продольная координата потока, м; t - время, с.
Стационарное уравнение турбулентной диффузии имеет вид:
V' — + Ц ^ = 0, (2)
' д' ' д'2
где концентрация С является функцией продольной координаты потока, поэтому
дС = о.
дг
Уравнение (2) может быть преобразовано к виду
SC д 2C — + —
д1 д1
Р— + — = 0. (3)
ГГ ^ "
При постоянном параметре p =- получим линеиное однородное уравнение
Di
второго порядка, которое допускает аналитическое решение. Будем искать решение в виде
C = ekl, где k = const , (4)
dC , kl тогда — = kekl; , dl dl2
kl. d!C = k2ekl.
Подставляя полученные выражения производных в уравнение (3), находим:
ekl (к2 + pk) = 0.
Так как e Ф 0, то
k2 + pk = 0.
Из характеристического уравнения (7) находим два корня
к = 0; k 2 =- p = - *.
Корни характеристического уравнения действительны и различны: к ф к2. В этом случае частными решениями будут функции
- ^
С1 = ekll = 1; C2 = ek2l = e * .
(5)
(6)
(7)
(8)
Эти решения линейно независимы, так как
С2 --1 ■ = e
С
2 - " Dl ф const.
Следовательно, общий интеграл имеет вид
C = Cp + Ccme Dl
(9) (10)
(11)
В данном случае величина Ср является фоновой концентрацией загрязнения. Примем фоновую концентрацию равной предельно допустимой концентрации (ПДК) нефтяного загрязнения в чистой воде, которая, согласно руководящему документу [3], составляет 0,1...0,2 мг/л=1,0..,2,0-10"4 кг/м3.
Величина Сст определяет концентрацию загрязнения в створе выброса. Характер изменения концентрации загрязнения в зависимости от расстояния от створа выброса вниз по потоку вычислен в среде MathCad (рис.1). При необходимости расчет может быть выполнен для конкретных значений параметров, полученных в результате натурных измерений.
Рис. 1. Вычислительная процедура и график функции концентрации загрязнения в зависимости от расстояния до створа выброса вниз по потоку
Аналогичный результат решения обыкновенного дифференциального уравнения одномерной турбулентной диффузии загрязнения (2) может быть получен в системе MathCad с помощью встроенной функции Odesolve (рис. 2).
Рис. 2. Решение дифференциального уравнения турбулентной диффузии загрязнения
с помощью функции Odeso've
Весьма странным является то, что в литературе, посвященной исследованиям распределения концентрации загрязнения водных потоков при разливах нефти и сбросах сточных вод [2, 4-6], рассматривается только стационарное уравнение турбулентной диффузии. Одномерное уравнение математической физики параболического типа (1) может быть решено численно, например, методом конечных разностей (рис. 3). Обозначения исходных данных решения уравнения соответствуют использованным выше (рис. 1 и 2).
Отметим, что неустойчивость операции численного дифференцирования методом конечных разностей вызывает необходимость применения вычислительной схемы Кранка-Николсона [11]. Впрочем, алгоритмы указанной вычислительной схемы достаточно широко известны [10].
Рис. 3. Решение в программной среде MathCad дифференциального уравнения нестационарной турбулентной диффузии загрязнения
Системный подход к анализу процессов на подводных переходах. Выбор адекватной модели исследуемого процесса - один из важных шагов в моделировании, поскольку процесс упрощения может привести к потере исходной информации, в то время как излишняя детализация приводит к созданию сложного, не имеющего решения алгоритма.
В системной балансовой модели вводятся понятия уровня накопления субстанций и темпа потока, представляющего расход этой субстанции [8, 9]. Уровни графически изображаются в виде прямоугольников, а темпы потока в виде вентилей (рис. 4).
Рис. 5. Системная диаграмма процесса распространения загрязнения в водном потоке
Рис. 4. Обозначения, принятые в модели Форрестера
Построим балансовую модель изучаемого процесса на основе системной диаграммы распространения нефтяного загрязнения в водном потоке, использовав обозначения: PDK - предельно допустимая концентрация загрязнения, кг/м3; С - концентрация нефтяного загрязнения в расчетном створе, кг/м3, КЯ, КБ, ко - коэффициенты разбавления, разложения и генерации загрязнения объекта, 1/м; ЯС, БС, GC - темпы разбавления, разложения и генерации загрязнения, кг/(м3-м) (рис. 5).
На основе системной диаграммы процесса распро-
странения загрязнения в водном потоке запишем уравнение темпов и уровней:
С = Ср + Сст.
?
С = С + Б1 (ОС - ЯС - БС); ОС = КО ■ С;
ЯС = КЯ ■ С;
БС = КБ ■ С
(12)
(13)
(14)
(15)
(16)
где Сст, Ср, С - концентрация загрязнения в створе выброса нефти, фоновая концентрация водного объекта и концентрация в расчетном створе ниже выброса нефти соответственно, кг/м3, Ср = РБК; РБК - предельно допустимая концентрация (ПДК) нефтяного загрязнения, кг/м3; ОС, ЯС, БС - темпы изменения концентрации загрязнения за счет залпового выброса нефти, разбавления и разложения загрязнения водотока соответственно, (кг/м3)/м; КО, КЯ, КБ - параметры генерации, разбавления и разложения загрязнения водного объекта соответственно, 1/м.
Программа и результаты расчета показателей распространения загрязнения в водном потоке в графической форме реализованы в среде MathCad (рис. 6).
Исходные параметры модели подобраны так, чтобы результаты расчета концентрации загрязнения по балансовой модели можно было сравнить с ее аналитическим выражением, носят иллюстративный характер. С
1:= 0.. 100 PDK =2-10 Ср = PDK Сет = 10 GC = 0 RC = 0 DC = 0 KG, KR, KD
- концентрация з агрязнения б расчетном створе, кгЛтГ1
- расстояние от створа выброса нефти до расчетного!, кг/м' _4 - пр ед ельно д опустимая к онцентрация
загрязнения в одного объекта (ПДК), кгЛ^
- фоновая концентрация загрязнения в одного объекта, кг/м^
- концентрация загрязнения в створе выброса нефти, кг/м5
- исходные темпы изменения концентрации загрязнение разбавления и разложения загрязнения водотока, (кг/ьг)/м;
С(1)
С Сет + Ср KG 0.01 KR 0.001 KD 0.011 for Die 0..1 GC KG-RC KR DC KD
параметры генерации, разбавления и разложения з агрязнения водотока,
10
С(1)
5
С i- С + D1 ■ (GC - RC - DC)
Рис. 6. Программа и результаты расчета зависимости концентрации нефтяного загрязнения
от расстояния до створа выброса нефти
Системная модель кислородного баланса в реке при наличии органических загрязнений. При попадании нефти в водоток начинается процесс ее окисления, на что затрачивается растворенный в воде кислород. В результате концентрация кислорода в воде может значительно понизиться, что приводит к эвтрофированию водного объекта.
В одномерном приближении процессы, формирующие кислородный баланс в реке при наличии нефтяных загрязнений, описываются дифференциальными уравнениями [5]: - биохимической потребности в кислороде:
А (<L) + —(QL) = - LdD, — I- k<L;
дГ ' J öl { l 8l 1 1
растворенного кислорода:
8 8 8 l 8c I / \
— (<ac)h--(Qc) = — I aDl — I - k1<L + k2ю (cp - c)
8
8
8t
8l
8l
8c
8l
(17)
(18)
где Ю - площадь поперечного сечения русла, м2; Q - расход воды, м3/с; kl - коэффициент скорости биохимического окисления нефти,1/с; k2 - коэффициент скорости реа-эрации, 1/с; Dl - коэффициент продольной диффузии, м2/с; L - биохимическое потребление растворенного кислорода, кг/м3; c - концентрация растворенного кислорода, кг/м3; cp - концентрация растворенного кислорода выше створа выброса нефти,
кг/м3. Для ряда практических задач можно предположить, что формирование качества воды не зависит от конвективного переноса и диффузии вещества, тогда процесс кислородного баланса в реке при наличии нефти описывается дифференциальными уравнениями Стритера - Фелпса [6]:
dt - klL
— = klL -k2D dt 1 2
(19)
где Б, L - дефицит и концентрация кислорода, кг/м ; £ - коэффициент скорости биохимического окисления нефти, 1/с; k2 - коэффициент скорости реаэрации, 1/с.
Величины коэффициентов £ и k2 можно определить по данным натурных исследований или по эмпирическим зависимостям.
Например, из уравнения Г. В. Стритера [5, 6] следует:
2,3, Сп
ki lg^0, х С '
(20)
где С0 - начальная концентрация загрязняющего вещества, кг/м3; Сх - концентрация того же вещества через время X (сут), кг/м3.
Коэффициент скорости реаэрации £2 характеризует скорость поступления кислорода из атмосферы при постоянной температуре. Этот процесс протекает до полного насыщения воды кислородом в пределах его растворимости при данных значениях температуры и давления.
Шедшилл, Элмор, Букингем [6] предложили определять коэффициент реаэрации £2 по зависимости:
k 2 = 5,75 -10~5Vl0 969 H
-1,673
(21)
где Vl - средняя скорость потока, м/с; Н - средняя глубина русла, м.
Для модели Стритера - Фелпса может быть построена системная диаграмма (рис. 7).
Рис. 7. Системная диаграмма модели Стритера - Фелпса
По системной диаграмме могут быть записаны конечно-разностные уравнения:
L = L - БТ ■ LD;
ЬБ = £1 Ь; Б = Б + БТ ■ (БО - ББ); БО=£1 Ь;
ББ = £ 2 ■ Б .
При БТ ^ 0 получим систему дифференциальных уравнений (19). Систему дифференциальных уравнений (19) можно также представить в виде:
(22)
(23)
(24)
(25)
(26)
— ,-k< L dt 1
dc dt
dL dt
+ k 2 (Cp - C)
(27)
Решение конечно-разностных уравнений в среде MathCad (рис. 8) позволяет построить зависимость дефицита кислорода от времени.
- текущее время, с
- период определения Б ПК, сут
- продольный коэффициент турбулентной диффузии, м/с
- средняя глубина потока, м
- продольная скорость водного потока, м^с
- фоновая концентрацил загрязнения в одного объекта, кг/м^
- концентрация загрязнения в створе выброса нефти, кг/м^
0..3600 г = 5
D1 = 10 Н = 2 VI = 0.5 Ср = 2■10 Сет = 10
D(t)
Li- 15 ■ 10 LD i- 0 , D 60 ■ 10 DG 0 DD ^ 0 kl^O k2 10 Сет Сет Cp Cp VI VI for DT e 0..t ti- DT 1 VI■t
L, LD - величина и скорость биохимического потребления растворенного кислорода, кг/М3 0,00, ОО - величина и скорости роста деградации дефицита кислорода, кгЛт!3 к1 - коэффициент скорости биохимического окисления нефти, 1/с к2 - коэффициент скорости реаэрации, 1/с 1 - продольная к о ордината потока, м
D
D1
Ct Ср + Сет ■ е
23 , (Сст ,п
kl i--log -,10
т 24 -3600 V Ct
0.02
n-5
yjO.969 ц-1.673
k2 5.75 ■ 10 L L - DT ■ LD LO kl ■ L
D D + DT ■ (DG - DD) DG kl ■ L DD ^ k2 ■ D
4000
Рис. 8. Программа и результаты расчета зависимости дефицита кислорода от времени
Приняв начальные условия при t = 0 (в створе выброса нефти) Ь(г,=0) = Ь0, Б(г,=0) = Б0, с (1=0) = с0, можно получить следующие решения систем уравнений (19) и (27)
L = Lo e k1t:
(28)
D =
k 2 - k1
Lo (
o ve-k1 - e-k2)+ Doe
t)+ Do
-k,t ■
с = cp -
k 2 - k1
Lo -k1t - e\-(cp - Co) e
' Y(cp- co^
(29)
(30)
где Ь), Б0, с0 - соответственно биохимическое потребление, дефицит и концентрация кислорода в начальный момент времени (в створе выброса нефти), кг/м3; ср - фоновая концентрация кислорода в створе выброса нефти, кг/м3.
k
Дефицит кислорода в створе реки определяется подстановкой величины /, соответствующей времени добегания загрязнения от створа выброса до контрольного, в зависимость (29).
Между конечным Dt и начальным D0 дефицитом кислорода существует следующая зависимость [3]:
Д = D010"*2'.
Рис. 9. Расчет расстояния от створа выброса нефти до максимального дефицита кислорода
Продифференцировав уравнение (29) и приравняв его к нулю, найдем время через которое наступит максимальный дефицит кислорода:
1
t =-
max 7 7
- k
ln
k
- D0 k2 - k1
L0 k1 J
(32)
1
Считая движение воды в реке равномерным с постоянной скоростью течения Vl,
можно найти расстояние lmax от створа выброса нефти до створа максимального дефицита кислорода (рис. 9) в виде:
l = Vt (33)
max l max
Системная динамическая модель оценки ущерба окружающей природной среде от загрязнения нефтепродуктами. В таблице идентификаторов приведены переменные, описывающие ущерб, наносимый природной среде от загрязнения почв и акваторий.
Обозначения элементов системы ущерба окружающей природной среде от загрязнения нефтепродуктами
Идентификатор Вербальное описание Размерность Численное значение
SP Масса углеводородов, испарившихся с поверхности почвы кг/м2 1,732 (t ° =20°С)
SPG Темп испарения углеводородов с поверхности почвы (кг/м2)/сут 0,072
SW Масса углеводородов, испарившихся с поверхности водоемов кг/м2 0,318 (t°=20°С)
SWG Темп испарения углеводородов с поверхности водоемов (кг/м2)/сут 0,003
FP Площадь загрязнения почв нефтепродуктами м2 1
FPG Темп увеличения площади загрязнения почв нефтепродуктами м2/сут 1,5
FW Площадь загрязнения акваторий нефтепродуктами м2 1
FWG Темп увеличения площади загрязнения водоемов нефтепродуктами м2/сут 100
TV Температура воздуха оС 20
TVG Темп изменения температуры воздуха оС/сут 0,1
TP Температура почвы оС 20
TPG Темп изменения температуры почвы оС/сут 0,1
TW Температура воды оС 20
TWG Темп изменения температуры воды оС/сут 0,1
CPW Ущерб природной среде от загрязнения атмосферы р 1
HB Базовый норматив за выброс загрязняющих веществ в атмосферу р/кг 1575
КИ Коэффициент индексации ущерба окружающей среде - 2
КЕ Коэффициент экологической ситуации и экологической значимости атмосферного воздуха - 1,1
РК Повышающий коэффициент - 5
SPT Множитель влияния температуры воздуха и почвы на массу нефти, испарившейся с почвы кг/(оСсут) 0,0036
SWT Множитель влияния температуры воздуха и воды на массу нефти, испарившейся с поверхности воды кг/(оСсут) 0,00015
В состав переменных включены семь уровней состояния, характеризующих физические процессы загрязнения преимущественно воздуха. Пять величин определяют значение экономической характеристики ущерба в зависимости от происходящих физических процессов загрязнения. Семь темпов изменения физических показателей характеризуют прирост (динамику) происходящих процессов. Два множителя влияния связывают подсистемы комплекса в единое целое. Таким образом, система оперирует с двадцатью одной переменной.
Переменные и постоянные параметры системы могут быть объединены в структурную диаграмму, на основании которой составляется система дифференциальных уравнений.
Элементы системы ущерба окружающей среде от загрязнения атмосферы определяются семью уровнями (SP, SW, FP, FW, TV, TP, TW) и семью темпами (SPG, SWG, FPG, FWG, TVG, TPG, TWG).
Ущерб природной среде (CPW) связан с нормативом за выброс загрязняющих веществ в атмосферу (НВ) и коэффициентами индексации ущерба окружающей среде (КИ) и экологической значимости атмосферного воздуха (КЕ), а также с повышающим коэффициентом (РК). Два множителя влияния (SPT и SWT) температуры связывают подсистемы между собой (рис. 10).
Рис. 10. Системная диаграмма оценки ущерба природной среде от загрязнения нефтепродуктами
Системная диаграмма позволяет записать уравнения темпов и уровней:
т=т+БТTWО; ТЖ=0.. .30; FP=FP+DT•FPО; FP=0.. ,102; Ш=Ш+БТШО;
FW=0... 104; ТР=ТР+БТТРО; ТР=0.. .30; ТУ=ТУ+БТТУО; ТУ=0.. .30; SP=SP+DTSPО;
SPО=0,5 SPT•(TУ+TP); Ш=Ш+БТШО; SWG=0,5 SWT•(TW+TУ);
CPW=PK•КЕ•КИ•НВ(SP+SW).
Подстановка численных значений переменных и постоянных величин системного комплекса (табл. 1) произведена на основании данных «Методики определения ущерба природной среде от разрушения нефтепроводов» и «Методических рекомендаций по определению платы за выбросы, сбросы, размещение загрязняющих веществ в природную среду» [1, 12].
Программа расчета показателей системного комплекса выполнена в среде MathCad (рис. 11).
CPW(t) :-
SPi- 1.732 SPGi- 0J072 SWi- 0318 SWG<-0J003 FPi- 1 FPG<- 15 Pi(- 1 FWGi- 100 TV<h 20 TVGi- 0.1 TP i-20 TP Cr*- 0.1 TW 20
тиге«- o.i
CPGi- 1 HB t- 1575 KUi-2
KE<- 1.1 PK<- 5 SPTi- 0Я036 SWT <-0.00015 for DT « 0..t
TW i-TW +DT TWG FPi- FP+DT FPG FWf- FW+ DT-FWG TPi- TP + DT-TP Cr TVi- TV+DT-TVG SPi- SP+DT-SPG SPGi- 05SPT-(TV+TP) SW<h SW + DT-SWG-SWS «- 0 5 SWT- (TW + TW) CPW PEKEKUHD (SP + S W) CPW
Рис. 11. Программа расчета показателей системного комплекса В имитационном режиме выполнена оценка динамики показателей системного комплекса (рис. 12). В качестве временного интервала прогноза выбран десятилетний период. Анализ динамики показателей указывает на их нелинейный рост. Например, масса углеводородов, испарившихся с поверхности почвы и акваторий, возрастает в несколько раз (рис. 12, а, б). Изменение площади загрязнений почв и водоемов также имеет тенденцию к возрастанию (рис. 12, в, г). При этом температуры воздуха, почвы и воды возрастают (рис. 12, д, е, ж).
в)
FP(t) 50
Жj
TW(t)25 -
Рис. 12. Динамика массы углеводородов, кг/м (а - испарившихся с поверхности почвы, б - испарившихся с поверхности водоемов, кг/м2), площади загрязнения нефтепродуктами, м2 (в - почв, г - акваторий), температуры, оС (д - воздуха, е - почвы; ж - воды), ущерба природной среде, р (з - от загрязнения атмосферы)
Указанные физические процессы осуществляются за счет положительных темпов роста их показателей.
Процессы загрязнения окружающей среды приводят к экономическому ущербу, определяемому затратами на ликвидацию последствий аварий нефтепроводов и разлива нефти. Экономический ущерб также является возрастающей во времени величиной (рис. 12, з).
Сбор нефти с поверхности реки при аварийных разливах. ВНИИСПТнефть разработан опытный образец нефтесборщика НВ-1 ТУ 39-1208-87, предназначенного для
сбора нефти с поверхности воды при аварийных разливах на подводных переходах через судоходные реки [7].
При всех достоинствах нефтесборщика ВНИИСПТнефть для условий малых рек необходимо предусмотреть стационарные системы сбора нефти при аварийных разливах. Например, может быть предложено нефтесборное устройство на основе гибкой бо-новой системы с потокообразующими лопастями (рис. 13).
При образовании нефтяного пятна 1 приводные барабаны 7 приводятся во вращение. При этом рабочая ветвь гибкого бона 2 приходит в движение и с помощью системы потокообразующих лопастей 4, соединенных с гибким боном при помощи шарниров 5, возбуждает поток, увлекающий нефть к нефтезаборнику 9. Холостая ветвь гибкого бона 3 движется встречно рабочей ветви 2. На холостой ветви бона потокообра-зующие лопасти 4 примыкают к бону, обеспечивая снижение гидравлического сопротивления. Гибкий экран 6, закрепленный на гибком боне, обеспечивает предотвращение ухода нефти вниз по течению. Для удержания на поверхности воды ветвей гибкого бона 2, 3 использованы сферические емкости с запасом плавучести 10.
1. В математической постановке задачи одномерного приближения распространения загрязнения вдоль потока составлено уравнение математической физики гиперболического типа, которое решено методом конечных разностей и с применением методов прогонки типа Кранка-Николсона для устранения неустойчивости решения при изменении коэффициента турбулентной диффузии в широком диапазоне.
2. В работе последовательно решены задачи стационарной турбулентной диффузии на основе прямого аналитического решения и реализована процедура решения стационарного уравнения средствами прикладного программного пакета MathCad с помощью встроенной функции Odesolve.
Рис. 13. Нефтесборное устройство: 1 - нефтяное пятно; 2 - рабочая ветвь гибкого бона; 3 - холостая ветвь гибкого бона; 4 - потокообразующая лопасть; 5 - шарнир раскрытия лопасти; 6 - гибкий экран; 7 - приводные барабаны; 8 - береговые опоры; 9 - нефтезаборник, 10 - плавучие буи
Выводы.
3. В связи с отсутствием в литературе прикладного характера эффективных методов решения уравнений нестационарной турбулентной диффузии как загрязнений вызванных сбросами сточных вод, так и нефтепродуктами, в работе приведен пример вычислительной процедуры Кранка-Николсона, доведенный до численных результатов.
4. В связи с тем, что описанные физические модели не определяют характер изменения концентрации нефтяного загрязнения в зависимости от расстояния до створа выброса на природном фоне, в работу включена также программа решения задачи в указанных условиях.
5. Процесс окисления нефти в водотоке вызывает процессы интенсивного потребления растворенного в воде кислорода, поэтому решение задачи формирования кислородного баланса, описанное в работе дифференциальными уравнениями биохимической потребности в кислороде и концентрации растворенного кислорода, в одномерном приближении показывает динамику этих показателей во времени.
6. С практической точки зрения весьма полезным является решение задачи по вычислению расстояния от створа выброса нефти до створа максимального дефицита кислорода. Отмеченная задача также решена в работе в приближении равномерной и постоянной скорости течения.
7. Непосредственное применение параметров, определенных вышеперечисленными средствами, является затруднительным, поскольку ущерб, наносимый природной среде от загрязнения акватории почв, необходимо учитывать в комплексе взаимосвязанных влияющих факторов, поэтому в работе составлена, формализована и программно реализована системная динамическая модель оценки ущерба окружающей природной среде от загрязнения нефтепродуктами. В программе расчета показателей системного комплекса, реализованной в среде MathCad, определены экономические характеристики ущерба природной среде. Анализ динамики массы углеводородов в различных природных объектах показал относительно простой характер функций связи параметров модели со временем. Результаты тестирования модели свидетельствуют о негативном влиянии процессов загрязнения водотоков на их экологическое состояние.
8. В связи с вероятностным характером нарушения целостности трубопровода подводного перехода расположение вниз по течению относительно створа нефтесборщиков известных конструкций является малоцелесообразным. По этой причине на несудоходных реках незначительной ширины целесообразно устанавливать либо стационарные, либо нестационарные средства задержания и сбора нефти с поверхности воды при аварийных разливах. В работе предложена схема нефтесборного устройства динамического типа, представляющего потокообразователь с распределенными рабочими элементами.
По итогам выполненной работы можно сформулировать перспективные направления анализа процессов по предмету исследования. Прежде всего, следует разработать методы информационно-технологического обеспечения решения задач трехмерной турбулентной диффузии с краевыми условиями не только первого рода, но также с условиями второго и третьего рода. Существующие в настоящее время параметры оценки водотоков должны быть экспериментально изучены на предмет их применимости к характеру динамики загрязнения водотоков нефтепродуктами.
Список литературы
1. Методика определения ущерба, причиненного окружающей природной среде при авариях на магистральных нефтепроводах. - М.: Минтопэнерго РФ, 1996. - 68 с.
2. Гросс, С.А. Определение расхода и времени вытекания жидкости из щели при разрыве стенки трубопровода / С.А. Гросс, Б.Г. Янов // Транспорт и хранение нефти и нефтепродуктов. - 1982. - № 11. - С. 9-10.
3. РД 39-0147098-005-88 «Правила охраны окружающей среды при сборе, подготовке и транспорте
нефти». - Уфа: ВостНИИТБ, 1988. - 25 с.
4. Вознесенский, Г.Ф. К оценке содержания нефтепродуктов в речной воде при разливах на подводных переходах/ Г.Ф. Вознесенский, И.А. Колосков // Транспорт и хранение нефти и нефтепродуктов. -1979. - № 6. - С. 32 - 35.
5. Караушев, А.В. Практические рекомендации по расчету разбавления сточных вод в реках, озерах, водохранилищах / А.В. Караушев, А.Я. Шварцман, М.А. Бесценная. - Л.: ГГИ, 1973. - 101 с.
6. Родзиллер, И.Д. Научные и инженерные основы прогнозирования качества воды водоемов и их защита от загрязнения сточными водами / И.Д. Родзиллер. - М.: Наука, 1976. - 247 с.
7. Спиридонов, В.П. Технические средства предотвращения загрязнения водоемов нефтью / В.П. Спиридонов, И.Д. Черкасов. - М.: ЦБНИТИ Минречфлота, 1983. - 48 с.
8. Дмитриев, Ю.А. Математическое моделирование экологических систем: Учебное пособие / Ю.А. Дмитриев, А.Г. Поздеев. - Йошкар-Ола: МарГТУ, 1997. - 206 с.
9. Форрестер, Дж. Мировая динамика / Дж.Форрестер. - М.: Наука, 1978. - 167 с.
10. Семененко, М.Г. Математическое моделирование в MathCad / М.Г. Семененко. - М: Альтекс-А, 2003. - 208 с.
11. Мэтьюс, Д.Г. Численные методы. Использование MATLAB / Д.Г. Мэтьюс, Д.К. Финк. - Изд. 3. Пер. с англ. - М.: Вильямс, 2001. - 720 с.
12. Методические рекомендации по определению платы за выбросы, сбросы, размещение загрязняющих веществ в природную среду. - М.: Минприроды РФ, 1993. - 56 с.
Статья поступила в редацию 19.01.11.
A. G. Pozdeev, Ju. A. Kuznetsova, Yu. V. Loskutov, A. E. Trukhan
ANALYSIS OF OIL-TRUNK PIPELINE UNDERWATER PASSAGES DESTRUCTION CONSEQUENCES AND THE INFLUENCE ON WATERWAY ECOLOGICAL STATE
The process of contaminated water dilution is considered and a model of oxygen balance of the basins, contaminated with oil emissions on the basis of system approach is composed. Oil pollution spreading in the waterway is estimated. A system model of assessment of damage to environment from oil products pollution is developed. In order to gather oil from a river surface in case of accidental spills, a skimmer construction design is offered.
Key words: ecological safety, system model, oil-trunk pipelines, underwater passages, skimmer.
ПОЗДЕЕВ Анатолий Геннадиевич - доктор технических наук, профессор, заведующий кафедрой водных ресурсов МарГТУ. Область научных интересов - проблемы водного транспорта и комплексного освоения водных ресурсов, математическое моделирование в гидродинамике и экологии. Автор более 60 публикаций.
E-mail: PozdeevA.G.@marstu.net
КУЗНЕЦОВА Юлия Анатольевна - кандидат технических наук, доцент кафедры водных ресурсов МарГТУ. Область научных интересов - исследование и моделирование русловых процессов в нижних бьефах гидроузлов, математическое моделирование в гидродинамике и экологии. Автор 14 публикаций, в том числе одной монографии и одного патента РФ.
E-mail: 2103Julia@mail.ru
ЛОСКУТОВ Юрий Васильевич - кандидат технических наук, доцент кафедры сопротивления материалов и прикладной механики. Область научных интересов - механика конструкций и материалов, методы математического моделирования, системы автоматизированного проектирования (в машиностроении). Автор 25 публикаций, в том числе одной монографии.
E-mail: LoskutovYV@marstu.net
ТРУХАН Артем Евгеньевич - аспирант кафедры водных ресурсов МарГТУ. Область научных интересов - экологическое состояние водотоков и водоемов, технические средства мелиорации сплавных путей. Автор одной публикации.
E-mail: kvr@marstu.net
