автореферат диссертации по безопасности жизнедеятельности человека, 05.26.03, диссертация на тему:Прогнозирование последствий аварийных истечений нефтепродуктов через дефектные отверстия из подземных продуктопроводов

005047173

На правах рукописи

(А'

КИСЕЛЕВ СЕРГЕЙ ЮРЬЕВИЧ

ПРОГНОЗИРОВАНИЕ ПОСЛЕДСТВИЙ АВАРИЙНЫХ ИСТЕЧЕНИЙ НЕФТЕПРОДУКТОВ ЧЕРЕЗ ДЕФЕКТНЫЕ ОТВЕРСТИЯ ИЗ ПОДЗЕМНЫХ ПРОДУКТОПРОВОДОВ

Специальность 05.26.03 - Пожарная и промышленная безопасность

(нефтегазовый комплекс)

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

1 з СЕН 2012

Уфа 2012

005047173

Работа выполнена в ООО «Волго-Урапьский научно-исследовательский и проектный институт нефти и газа» (ООО «ВолгоУралНИПИгаз»)

Научный руководитель - доктор технических наук

Клейменов Андрей Владимирович

Официальные оппоненты: - Идрисов Роберт Хабибович,

доктор технических наук, профессор, ГУП «Институт проблем транспорта энергоресурсов», заведующий отделом «Безопасность сложных технических систем»

- Тляшева Резеда Рафисовна,

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

Ведущее предприятие - РГУ нефти и газа имени И.М. Губкина

Защита состоится 20 сентября 2012 г. в 1400 часов на заседании диссертационного совета Д 222.002.01 при Государственном унитарном предприятии «Институт проблем транспорта энергоресурсов» (ГУП «ИПТЭР») по адресу: 450055, г. Уфа, проспект Октября, 144/3.

С диссертацией можно ознакомиться в библиотеке ГУП «ИПТЭР».

Автореферат разослан 20 августа 2012 г.

Ученый секретарь

диссертационного совета

доктор технических наук, профессор

— Худякова Лариса Петровна

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы

При эксплуатации промысловых и магистральных продуктопроводов одной из актуальных проблем является обеспечение безопасной эксплуатации линейной части, которая решается, главным образом, за счет своевременного проведения восстановительных и ремонтных работ на основе данных диагностики. Однако проводимые мероприятия не могут полностью исключить возможность возникновения аварийных ситуаций. Тяжесть последствий аварий зависит от оперативности принятия мер по локализации и ликвидации разлившихся нефтепродуктов, а также готовности эксплуатирующей организации к таким авариям. Для этих целей осуществляется прогноз масштабов возможных аварийных ситуаций, их последствий для социальной сферы и окружающей среды.

Существующие модели и методики позволяют оценить последствия от уже произошедших аварий или, применительно к проливам нефтепродуктов, дать примерные характеристики геометрии, площади разлива и объемов загрязненных грунтов. В настоящее время отсутствует доступный инженерный метод моделирования разлива жидких углеводородов при истечении нефтепродуктов из подземных трубопроводов через дефектные отверстия максимальными линейными размерами до 0,75 от диаметра трубопровода. При нарушении целостности заглубленного трубопровода существенно затруднена количественная оценка ареала распространения нефтепродукта в грунте, поскольку данный вид аварий характеризуется продолжительностью во времени от момента начала утечки до момента ее обнаружения и остановки транспортировки.

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

Цель работы - обеспечение пожарной и промышленной безопасности трубопроводного транспорта путем разработки метода прогнозирования объемного загрязнения почвогрунтов при аварийном истечении жидких

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

Основные задачи работы:

- сбор и систематизация характеристик аварий, анализ существующих методов и подходов к оценке объемов загрязненных почвогрунтов и методов расчета массового расхода жидких углеводородов из типичных дефектных отверстий в теле трубопроводов при аварийном истечении;

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

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

- проверка адекватности разработанного метода на примере реально произошедших аварий.

Методы решения поставленных задач

Применялись методы аналитической оценки движения жидких углеводородов в пористой среде, скорости их фильтрации и действительной скорости движения жидкости в грунте; линейный закон Дарси; методы определения фильтрационных параметров пористых сред - коэффициентов фильтрации и проницаемости; методы построения объективных зависимостей между емкостными и фильтрационными параметрами пористых сред (между коэффициентами пористости и проницаемости); математические методы моделирования процессов фильтрации и процессов формирования опасных факторов техногенных аварий, а также численное решение комплексов математических уравнений и их практическая проверка на основе материалов реализованных аварий.

Научная новизна результатов работы:

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

- определены условия перехода подземного истечения на стадию поверхностного истечения;

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

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

На защиту выносятся:

- усовершенствованная теоретическая модель подземного фильтрационного движения жидких нефтепродуктов в месте прокладки продуктопроводов;

- модель подземного истечения жидких углеводородов из продуктопроводов, разработанная с учетом закономерностей строения и особенностей грунтов в районе их прокладки;

- метод оценки последствий аварийных истечений нефтепродуктов через дефектные отверстия максимальными линейными размерами до 0,75 от диаметра трубы;

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

Практическая ценность и реализация результатов работы

Метод может быть использован при оценке рисков объектов трубопроводного транспорта нефтепродуктов и адекватном прогнозировании потенциального ущерба, нанесенного социальной сфере и окружающей природной среде, при проектировании объектов нефтегазового комплекса. Метод опробован при разработке Деклараций промышленной безопасности к проектам обустройства Филипповской и Среднекаменноугольной залежей Оренбургского нефтегазоконденсатного месторождения ООО «Газпром добыча Оренбург».

Результаты оценки последствий распространения нефтепродуктов в толще грунта позволили уточнить масштабы возможных загрязнений почвогрунтов, разработать эффективные Планы локализации и ликвидации аварийных разливов нефти (ПЛАРН), обосновать состав сил и средств для

проведения работ по локализации и ликвидации аварийных разливов нефти и нефтепродуктов для нефтегазодобывающих предприятий.

Апробация результатов работы

Основные научные положения и практические результаты работы неоднократно доложены, обсуждены, одобрены и рекомендованы к применению на международных и российских научно-технических конференциях и симпозиумах, включая VI Всероссийскую научно-техническую конференцию «Актуальные проблемы развития нефтегазового комплекса России» (Москва, 2008); VIII Всероссийскую научно-техническую конференцию «Актуальные проблемы состояния и развития нефтегазового комплекса России», посвященную 80-летию Российского государственного университета нефти и газа имени И.М. Губкина (Москва, 2010); Ш и IV научно-технические конференции с международным участием «Основные проблемы освоения и обустройства нефтегазовых месторождений и пути их решения» (Оренбург, 2009,2010).

Публикации

Основные результаты диссертационной работы опубликованы в 14 научных трудах, в том числе 6 статей в изданиях, рекомендованных ВАК Министерства образования и науки РФ.

Структура и объем работы

Диссертационная работа состоит из введения, четырех глав, основных выводов, библиографического списка использованной литературы, включающего 104 наименования, и 2 приложений. Работа изложена на 140 страницах машинописного текста, содержит 11 таблиц, 41 рисунок.

Благодарности

Автор выражает искреннюю благодарность сотрудникам отдела промышленной и экологической безопасности ООО «ВолгоУралНИПИгаз»: зав. сектором Даньшиной O.A., зав. сектором Донецковой A.A.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе представлен краткий анализ аварийности на трубопроводном транспорте нефти и нефтепродуктов за период с 2003 по 2009 годы на территории Российской Федерации. Данные получены из отчетов о деятельности Федеральной службы по экологическому, технологическому и атомному надзору.

Выявлено, что более 60 % аварий (коррозия, брак строительно-монтажных работ, заводской брак) от общего их количества приходится на аварии, связанные с выходом нефтепродуктов на (дневную) поверхность через дефектные отверстия искусственного или естественного происхождения максимальными линейными размерами до 0,75 от диаметра трубопровода.

На рисунке 1 представлена диаграмма удельной аварийности на нефтепродуктопроводах за период с 2003 по 2009 годы. За рассматриваемый временной период удельная аварийность на нефтепродуктопроводах составила около 0,2 аварий на 1000 км в год.

Рисунок 1 - Удельная аварийность на нефтепродуктопроводах за период с 2003 по 2009 годы

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

Известно, что основными характеристиками пористой среды, определяющими процесс движения жидкости в ней, являются открытая

пористость и проницаемость грунта. Представление о структуре порового пространства дает рисунок 2.

Рисунок 2 - Срез пористой среды

Впервые закон фильтрации жидкости в пористом грунте экспериментально установил французский инженер Анри Дарси. Впоследствии этот закон назван в его честь. Дарси установил, что объемный расход жидкости пропорционален разности уровней (напору) АН и площади сечения образца 5 и обратно пропорционален длине образца Ь. Закон Дарси описывает линейную фильтрацию.

Существенный вклад в развитие теории фильтрации в разное время внесли Буссинеск Ж., Форхгеймер Ф. Слихтер Ч. впервые предложил модели идеального и фиктивного грунтов и показал, что пористость и просветность фиктивного грунта зависят не от диаметра частиц, а лишь от плотности их укладки.

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

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

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

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

Изучению негативного воздействия на окружающую природную среду опасных производственных объектов, разработке методов оценки риска и снижения отрицательных последствий аварий, в том числе определению объемов загрязненных грунтов при авариях на нефтепродуктопроводах, посвящены работы отечественных ученых Ларионова В.И., Котляревского В.А., Брушлинского H.H., Сущева С.П., Кутукова С.Е., Павлова C.B., Гумерова P.C., Гумерова А.Г., Елохина А.Н., Идрисова Р.Х., Швыряева A.A., Прусенко Б.Е., Козлова М.А. и др. В разное время разливами нефти занимались в УНИ, ВНИИСПТнефть, РГУ нефти и газа имени И.М. Губкина, БашНИИ, СибНИИНП, УфНИИ и пр. Указанные работы в большинстве своем направлены на оценку площади и глубины инфильтрации нефтепродуктов в грунт при поверхностном истечении.

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

При вычислениях производительности аварийного источника истечения через дефектное отверстие учитываются технические характеристики трубопровода, размещение запорной арматуры и свойства транспортируемой нефти. Согласно нормативно-методической литературе процесс аварийного истечения нефтепродуктов делится на три характерных периода: период аварийного истечения под рабочим давлением до отключения насосов (Vj\ период между отключением насосов и закрытием запорно-регулирующей арматуры (V2), период самотека после закрытия задвижек (Vj). Полный объем выброса равен сумме объемов за все три периода.

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

На практике наиболее распространенным способом определения скорости выхода нефтепродуктов через дефектное отверстие является использование зависимости Бернулли-Эйлера. Максимальный объемный расход прямо пропорционален площади дефектного отверстия, скорости выхода жидкости через дефектное отверстие и коэффициенту расхода, который зависит от режима течения жидкости и формы дефектного отверстия.

При адаптации методов теории фильтрации к решению поставленных задач исходим из следующих предположений:

- все дефекты линейной части трубопровода приведены к типовым аварийным отверстиям: свищ, трещина, гильотинный разрыв;

- в начальный момент времени при «внезапном» появлении подземного источника истечения фильтрационное движение будет подчиняться нелинейному закону фильтрации - закону Форхгеймера, но протяженность во времени данного периода будет незначительной. Исходя из этого, считаем, что фильтрация будет линейной и описывается законом Дарси;

- плотность грунта в траншее увеличивается ко дну траншеи, соответственно уменьшаются пористость и проницаемость;

- основное распространение жидких углеводородов в толще грунта будет происходить по профилю траншеи;

- движение грунтового потока будет вызываться перепадом давления, действием гравитационных и капиллярных сил;

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

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

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

не только от координат дг^иг.нои от времени г.

Модифицируя уравнение Буссинеска для описания фильтрации в ненасыщенном грунте с известными расходом и давлением в начальный момент времени, получим:

дх2 ду2 дг2

гДе Чх. Чу. Чг - проекции объемного расхода в направлении осей х, у иг соответственно, м3/с; ипот - потеря скорости, м/с; 0 - производительность аварийного источника, м3/с; 5 - площадь фильтрации, м2. Представленное уравнение описывает распределение скоростей фильтрации в объеме расчетной области. В качестве физического закона, выражающего скорость движения жидкости в пористом грунте под действием градиента давления, принимаем закон Дарси.

Третья глава диссертации посвящена разработке метода оценки последствий загрязнения грунтов при разрушении подземного продуктопровода на примере характерных почв и грунтов в коридорах прокладки трубопроводов западной части Оренбургской области.

На территории Оренбургской области почвенный покров представлен многими типами, подтипами, родами, видами и разновидностями и отличается большим разнообразием и сложностью. Преобладают три типа: черноземный, каштановый и солонцовый. Черноземы четырех подтипов: выщелоченные, типичные, обыкновенные и южные. Морфологический профиль чернозема слагается из пяти генетических горизонтов: А-АВ-В-ВС-С. Мощность почв такого типа составляет в среднем 1,1. ..2,0 м. Известно, что реальные пористые грунты имеют сложное строение пустотного пространства, обусловленное гранулометрическим составом, который образует поровые каналы, резко изменяющие диаметр и направление.

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

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

их укладки. Поскольку глубина заложения продуктопроводов колеблется от 1,0 до 2,0 м в зависимости от диаметра трубы, получение характеристик гранулометрического состава грунта возможно исходя из знания морфологических профилей конкретных типов почв (таблица 1).

Таблица 1 - Гранулометрический состав типичных черноземов

Глубина, Содержание фракций, %, с размерами частиц, мм

м 1,00...0,25 0,25.-0,05 0,05... 0,01 0,010...0,005 0,005... 0,001 <0,001 <0,01

0...0.1 2,2 2,8 34,2 12,3 15,4 33,1 60,8

0,2...0,3 2,3 3,1 32,2 11,1 17,3 34,0 62,4

0,4...0,5 2,1 3,9 30,8 12,1 15,9 35,3 63,3

0,6...0,7 2,7 4,9 28,9 11,8 15,4 36,4 63,5

0,8...0,9 3,3 3,5 30,0 13,0 13,4 36,7 63,1

1.0...1.1 4,7 5,0 23,4 11,6 14,5 40,8 66,9

1,2...1,3 3,8 5,2 24,6 11,2 17,1 38,3 66,5

1,4...1,5 4,5 4,7 24,4 11,3 15,0 40,2 66,4

Ср. состав 3,2 4,1375 28,5625 11,8 15,5 36,85 64,15

Принято, что смешанный грунт имеет одинаковый гранулометрический состав по всей глубине, но открытая пористость и проницаемость грунта будут уменьшаться ко дну траншеи из-за более плотной укладки частиц и, соответственно, увеличиваться плотность. В «Рекомендациях по методике лабораторных испытаний грунтов на водопроницаемость и суффозионную устойчивость» рекомендуется принимать значения пористости п = 0,26 для предельно плотной и п = 0,48 для рыхлой упаковок.

Коэффициент фильтрации (кф, м/с) с учетом природной влажности (W, д.е.) и коэффициент проницаемости (кп, м2) зернистых фунтов определяем по формулам:

5, -Ф-n-g 2 , =кф'М (2)

кф=-— "о > Кп ~ „ > к >

36-v g-p

где Ф - коэффициент формы сечений поровых каналов, Ф-0,45-(1+0,3-Б);

Б - балл окатанности частиц, составляющих грунт, по шкале Хабакова A.B.

(для частиц угловатой формы Б = 0); м ~ динамическая вязкость

нефтепродукта, Пас; р - плотность нефтепродукта, кг/м3; d0 - гидравлически

эквивалентный диаметр поровых каналов в разнозернистом грунте, мм, с10 =0,455- с • й?17; V- кинематическая вязкость, см /с; Бг - коэффициент водонасыщения, д.е., Бг = (1- п)-\¥-ргр/(п -рж)\ с - коэффициент пористости, д.е., с = п /(1- п)\ п - пористость, д.е.; Т1Х = й?60 / ¿10 - коэффициент разнозернистости грунта; й?10, <Ии и й?60 - диаметры частиц, соответствующие 10, 17 и 60 % содержанию от суммарного относительного состава всех частиц, которые определяются по логарифмическому графику гранулометрического состава (кривой зернового состава грунта), мм.

Построение трехмерной модели подземного истечения жидких углеводородов основывается на методе конечных элементов (МКЭ), т.е. рассматриваемая область разбивается на элементарные объемы конечного размера. Дефектное отверстие на теле трубопровода является началом отсчета 3-мерной системы координат, рисунок 3. Таким образом, исходная информация представляет собой 3-мерный массив g . Размеры всех элементов одинаковы, каждый элемент массива g (г, j, к) несет в себе информацию о фильтрационных характеристиках элементарного объема грунта в зависимости от влажности и открытой пористости.

2

вО».Ь к) элементарный объем;

-««©ив направления распространения нефтепродуктов в грунте.

1....... г

Рисунок 3 - Модель подземного истечения жидких углеводородов

Предполагаем, что распространение нефтепродуктов в грунте происходит в нарушенной во время укладки трубопровода области, т.е. в траншее, размеры которой принимаем исходя из характеристик трубопровода по СНиП 2.05.06 -85* «Магистральные трубопроводы». На основе данных о трассе прокладки трубопровода определяемся с типами грунтов, характерных для данного района.

Из материалов инженерных изысканий получаем данные о профиле прохождения трассы трубопровода и морфологическом составе грунтов, при

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

Рисунок 4 - Схема траншеи для построения модели подземного истечения жидких углеводородов

Исходя из характеристик грунтов рассчитываем по формулам (2) их открытую пористость и коэффициенты проницаемости от поверхности до дна траншеи. Грунты вне траншеи считаем условно непроницаемыми для нефтепродуктов. Коэффициент проницаемости кп для них примем равным 1,0-КГ15 м2.

При нахождении аварийного отверстия на верхней образующей трубы в результате интенсивного суффозионного процесса и образования котлована возможен выход нефтепродукта на дневную поверхность. Длина врубовой щели (по Никонову Г.П.) находится из следующего соотношения:

АР

1щ=с1В\--0,2),

<7„

(3)

где 1Щ - глубина прорези, м; с1- эквивалентный диаметр дефектного отверстия, м; ДР - давление перед отверстием, Па; ас - прочность грунта на сжатие, Па (ас ~ 1,5 МПа для почвенных грунтов); В - экспериментальная константа, зависящая от грунта (В ~ 0,5 для несвязанных пылеватых грунтов).

Время на образование врубовой щели и размыв объема котлована определяем по формулам:

/ . = -

Щ

Q■pж■lщ

*к =

Ргр'ё'Ук

(4)

где Q — объемный расход нефтепродукта, определяется по зависимости Бернулли-Эйлера, м3/с; рж, ргр — плотности соответственно нефтепродукта и грунта, кг/м3; g — ускорение свободного падения, м/с2; Ущ — объем грунта в щели, м3, Ущ = luf-Scmp; Vk~ объем грунта в котловане, м3, VK = к-1щ/[4-(Н-hj}, Scmp - площадь струи, рассчитываемая исходя из формы и размеров дефектного отверстия, м2.

Движение жидкости в поровом пространстве вызывается перепадом (градиентом) давления. Определяем градиент давления на /«-итерации с помощью следующей формулы:

АРт=^-ДР, (5)

т

Sm = 2 • (2 ■ т+1)2 (bxby + bxbz + byb.) при (2 ■ m+ l)by <L или (2 ■ m+ l)b= < H,

Sm = 2 (L-H + bx-(2-m+l)-(H + L)) в остальных случаях. Здесь АР - давление перед дефектным отверстием, Па; АРт - градиент давления на тя-итерации, Па; Sm - площадь поверхности фильтрации на /я-итерации, м2; Some — площадь отверстия, м2; L,H— соответственно ширина и глубина траншеи, м; bx, Ьу, Ъг — линейные размеры элементарного объема, м.

Проекции расхода q на оси х,у, z в уравнении (1) являются скоростями в соответствующих направлениях. На их основе строим поле распределения скоростей и(х, у, z) в рассматриваемой области с использованием закона Дарси в постановке относительно скорости движения грунтового потока:

т _ I.y.z ,¿-1 ux,y,z --;-> w

"X,y,Z т

где bxy z — линейные размеры элементарного объема пространства g(i, j, k), м;

— коэффициент проницаемости грунта, м2; ц - динамическая вязкость нефти, Па-с; Д Р'„ — градиент давления, аналог напора p g-AH в законе Дарси, Па, с учетом силы тяжести и капиллярных сил, А = Р± р- g- A h-sin у + Д Ркап; у — угол уклона дна траншеи. Знак «±» означает движение нефтепродукта вниз или вверх.

Время движения через элементарный объем в направлениях х, у, z на /я-итерационном шаге рассчитываем по следующим формулам:

с=—. С=%- (7)

* итх у и; итх

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

т ш т

(8)

ООО

Выход из расчета происходит в случае достижения значения накопленного итерационного времени, равного значению времени, установленному в расчете (1рас). Причем выход из расчета происходит раздельно для каждого из направлений {х,у, г).

При выходе нефтепродуктов на поверхность происходит их растекание по рельефу. Интенсивность истечения (0л) будет обусловлена глубиной заложения трубопровода и плотностью нефтепродукта.

Массу нефти или нефтепродуктов, испарившихся с поверхности земли, определяем исходя из площади пропитанного грунта состава и свойств

жидкости, параметров окружающей среды (температуры воздуха, скорости ветра) по удельной величине выброса вещества с единицы площади за единицу времени (кг/(м2-с)). Численно площадь Брп равна площади проекции области нефтяного загрязнения на дневную поверхность на т-итерационном шаге:

(Ь-Ь \-L\-L

5=3,14

(9)

где Ly — полуширина области загрязнения, м; Lx — полудлина области загрязнения в сторону понижения трассы трубопровода, м; —Lx — полудлина области загрязнения в сторону повышения трассы трубопровода, м.

Площадь пропитанного грунта и количество выделившихся в атмосферу паров нефти или нефтепродуктов в результате разлива необходимы при расчетах поражающих факторов от возможных сценариев аварий, таких как взрыв топливно-воздушных смесей (TBC), «пожар пролива», «пожар -вспышка», и зон токсического заражения окружающего пространства. При этом используем утвержденные соответствующим образом модели и методики.

Уравнение (1) решаем прямым численным методом с привлечением математической системы MathCad и MS Excel, учитывая принятые допущения.

В вычислительном эксперименте приняты следующие исходные данные: вещество - метанол; плотность р = 792 кг/м3; ширина траншеи L = 1,0 м; глубина Н = 1,0 м; давление в трубопроводе — 2,5 МПа; площадь дефектного отверстия Some = 0,0002 м2; к = 4,7-10'11... 1,76-10"11 м2; ц = 0,58-10 "3 Па-с; п = 0,26...0,48. Количество итераций 3000, размер элемента - 0,001 м3 (линейный размер —0,1 м).

В результате расчета получено: объем загрязненного грунта — 78 м3, объем вышедших нефтепродуктов —12 м3.

Метод оценки последствий аварий, связанных с истечением нефтепродуктов через дефектные отверстия из подземных продуктопроводов, реализован в математической среде MathCad и MS Excel. Проект состоит из пяти блоков:

1. Ввод исходных данных с автоматизированным формированием матриц с фильтрационными характеристиками грунта;

2. Определение расчетным путем характеристик нарушенного во время укладки трубопровода грунта в среде MS Excel;

3. Построение трехмерной модели траншеи, формирование матриц фильтрационных характеристик грунта, полей градиентов давлений и скоростей в элементарных ячейках в среде MS Excel;

4. Расчет объема загрязненного грунта, площади контакта нефтяного загрязнения на границе «грунт — воздух» и количества испарившихся нефтепродуктов с поверхности пропитанного грунта для оценки опасных факторов при взрыве TBC и токсического заражения окружающего пространства;

5. Визуализация процесса загрязнения жидкими углеводородами толщи пористого грунта посредством возможностей математической среды MathCad (рисунок 5).

Остановка расчета происходит в случае достижения значения расчетного времени, равного значению заданного времени для моделирования.

При реализации метода применялись возможности символьной математики и языковые средства MathCad и MS Excel.

Рисунок 5 — Распространение метанола в грунте (расчет на 24 часа)

В четвертой главе диссертации представлены результаты апробации метода оценки последствий аварий, связанных с истечением нефтепродуктов через дефектные отверстия из подземных продуктопроводов, на примере инцидента, произошедшего на IV нитке конденсатопровода Оренбург -Салават - Уфа (на 57 км) в результате разрушения сварного шва.

Картина распространения конденсата в почвогрунтах на месте рассматриваемого инцидента обусловлена нарушением сплошности коренного грунта в период строительства III и IV ниток конденсатопровода. По этой причине конденсат в грунте распространялся не равномерно, а по наиболее проницаемым участкам, образованным при отрытии трассовой траншеи и последующей ее засыпке, расположение которых имеет очаговый характер. Данное обстоятельство привело к тому, что более тщательное и оперативное удаление из зоны инцидента загрязненного грунта потребовало также и удаления незагрязненного грунта, находящегося между этими скоплениями.

При моделировании последствий инцидента для расчета было принято: давление в трубе - 2,5 МПа; площадь отверстия - 0,0004 м2; динамическая вязкость конденсата - 0,2-10"3 Пас; проницаемость - (1,14...6,14)-10"п м2. В результате расчета были получены размеры зоны объемного загрязнения: ширина - 35,2 м; протяженность - 37,2 м; высота - около 5 м; объем загрязнения, с учетом грунта над пятном загрязнения, составил около 6 тыс. м3 (рисунок 6).

Рисунок 6 - Распространение конденсата в грунте (расчет на 24 часа)

Проведенное «шурфование» по контуру загрязнения с отбором проб и анализом на содержание конденсата показало, что основная часть конденсата находится на глубине 4,0...5,5 м, площадь распространения конденсата -1400 м2. Объем загрязненного грунта между III и IV нитками конденсатопровода, в полосе протяженностью 40 м и шириной 35 м, подлежащий утилизации, составил около 7,4 тыс. м3. Разница между расчетным объемом и объемом фактически подверженного утилизации грунта составила около 19 %.

На основе технологических параметров трубопроводов системы сбора продукции на Филипповской залежи Оренбургского НГКМ и характеристик нефти было смоделировано аварийное истечение нефти через дефектное отверстие эквивалентным диаметром 0,025 м из выкидного нефтепровода от СРВ - 6 до УКПГ - 15. Диаметр трубопровода 114 мм, давление до 4 МПа, транспортируемая продукция - пластовая нефть.

Длина полосы пропитанного на поверхности грунта составила около 25 м, ширина - около 5 м. Количество нефтепродуктов, испарившихся с площади 125 м2 за 3600 секунд, рассчитывалось по Приложению 1 к ПБ 09-540-03 «Общие правила взрывобезопасности для взрывопожароопасных химических, нефтехимических и нефтеперерабатывающих производств» и составило около 650 кг. Количество участвующего во взрыве TBC вещества, согласно нормативно-методической литературе, принимаем равным 10 % от объема

испарившейся нефти - 65 кг. Расчеты зон нижнего концентрационного предела воспламенения (НКПВ), зон действия избыточного давления выполнены по ГОСТ Р 12.3.047-98.

Диаметр возможной зоны распространения «пожара-вспышки» составляет 44 м, высота - 3,0 м. Радиус пороговой зоны поражения человека -23 м. Зоны негативного воздействия опасных поражающих факторов от смоделированных аварий на примере сценариев «пожар-вспышка» и взрыв TBC представлены на рисунке 7.

Рисунок 7 - Зоны действия поражающих факторов «пожара-вспышки» и взрыва ТВС

ОСНОВНЫЕ ВЫВОДЫ

1. Выявлено, что наибольшая доля аварий, связанных с проливами нефтепродуктов, происходит через дефектные отверстия максимальными линейными размерами до 0,75 от диаметра продуктопровода. Проведен анализ существующих методов оценки объемов загрязненных почвогрунтов при аварийных разливах нефтепродуктов.

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

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

4. Проведены апробация программного средства и сопоставление расчетных результатов по предложенному методу оценки геометрических параметров области загрязнения с таковыми реализованной аварии, произошедшей на IV нитке конденсатопровода Оренбург — Салават - Уфа (расчетный объем загрязненного грунта по предложенному методу составил около 6 тыс. м3, по фактическим данным - около 7,4 тыс. м3, разница составила около 19 %). Метод опробован для получения исходных данных при оценке зон действия опасных факторов аварийных сценариев (теплового воздействия от пожаров пролива, барического воздействия от взрывов топливно-воздушных смесей и токсического загрязнения окружающего пространства сероводородом) на выкидных трубопроводах системы сбора продукции в проектах обустройства Филипповской и Среднекаменноугольной залежей Оренбургского НГКМ.

Основные результаты работы опубликованы в следующих научных трудах:

Ведущие рецензируемые научные журналы

1. Киселев С.Ю., Клейменов A.B., Гендель Г.Л. Методика моделирования аварийного истечения жидких углеводородов из магистральных и промысловых продуктопроводов // Нефтепромысловое дело. - 2006. - № 12.-С. 27-33.

2. Киселев С.Ю., Клейменов A.B., Гендель Г.Л. Оценка площади разлива нефтепродуктов при их истечении из наземного трубопровода // Защита окружающей среды в нефтегазовом комплексе. — 2007. — № 6. — С. 30-33.

3. Киселев С.Ю., Клейменов A.B., Гендель Г.Л. Построение модели аварийного истечения жидких углеводородов из продуктопроводов // Нефтепромысловое дело. — 2007. - № 12. - С. 84-85.

4. Киселев С.Ю., Клейменов A.B., Гендель Г.Л. К вопросу моделирования аварийного истечения жидких углеводородов из подземных продуктопроводов // Нефтепромысловое дело. — 2009. - № 8. — С. 73-75.

5. Киселев С.Ю., Клейменов A.B., Гендель Г.Л. Основные задачи при моделировании аварийного истечения жидких углеводородов из подземных продуктопроводов // Защита окружающей среды в нефтегазовом комплексе. — 2010. -№ 12.-С. 81-83.

6. Короленок А.М., Киселев С.Ю., Клейменов A.B. Модель оценки объема замазученных почвогрунтов при аварийном истечении нефтепродуктов из подземных продуктопроводов // Нефть, газ и бизнес. — 2011. — № 7. - С. 44-46.

Прочие печатные издания

7. Гендель Г.Л., Киселев С.Ю., Клейменов A.B. Моделирование процесса аварийного истечения жидких углеводородов из магистральных и промысловых продуктопроводов // Современные наукоёмкие технологии. — 2005. —№ 11, —С. 31-32.

8. Киселев С.Ю., Гендель Г.Л., Клейменов A.B. Методика оценки объема загрязнения грунта при истечении жидких углеводородов из подземного трубопровода (диаметр дефектного отверстия один дюйм) // Современные наукоёмкие технологии. — 2006. — № 5. — С. 36-37.

9. Клейменов A.B., Киселев С.Ю. Моделирование процесса аварийного истечения жидких углеводородов из магистральных и промысловых продуктопроводов // Актуальные проблемы состояния и развития нефтегазового комплекса России. Тез. докл. 7-ой научн.-техн. конф. 29-30 января 2007 г. - М.: РГУ нефти и газа, 2007. - С. 487.

10. Киселев С.Ю., Клейменов A.B., Гендель Г.Л. К вопросу моделирования аварийного истечения жидких углеводородов из подземных продуктопроводов // Основные проблемы освоения и обустройства нефтегазовых месторождений и пути их решения. Тез. докл. III научн.-техн. конф. с международным участием 21-22 мая 2009 г. - Оренбург: ООО ВУНИПИГАЗ, 2009. - С. 78-79.

11. Киселев С.Ю., Клейменов A.B., Гендель Г.Л. К вопросу моделирования аварийного истечения жидких углеводородов из подземных продуктопроводов // Актуальные проблемы состояния и развития

нефтегазового комплекса России. Тез. докл. VIII Всеросс. научн.-техн. конф., посвященной 80-летию Российского государственного университета нефти и газа имени И.М. Губкина, 1-3 февраля 2010 г. - М.: ИЦ РГУ НиГ, 2010. -С. 94-95.

12. Киселев С.Ю., Клейменов A.B., Гендель Г.Л. Основные задачи при моделировании аварийного истечения жидких углеводородов из подземных продуктопроводов // Основные проблемы освоения и обустройства нефтегазовых месторождений и пути их решения. Тез. докл. IV научн.-техн. конф. с международным участием 19-20 августа 2010 г. - Оренбург: ООО ВУНИПИГАЗ, 2010. - С. 56-57.

13. Киселев С.Ю., Клейменов A.B., Гендель Г.Л. Методика оценки площади загрязнения почв при истечении жидких углеводородов из наземного трубопровода // Сб. научн. тр. учеников научной школы Генделя Г.Л. / ООО «ВолгоУралНИПИгаз». - Оренбург: Оренбургская губерния, 2010. - С. 91-97.

14. Киселев С.Ю., Клейменов А.П. Оценка опасных факторов аварий, связанных с истечением нефтепродуктов из подземных трубопроводов // Энергоэффективность. Проблемы и решения. Матер. XI Всеросс. научн.-практ. конф. 19 октября 2011 г. - Уфа, 2011. - С. 232-235.

Фонд содействия развитию научных исследований. Подписано к печати 16.08.2012 г. Бумага писчая. Заказ № 208. Тираж 100 экз. Ротапринт ГУП «ИПТЭР». 450055, г. Уфа, пр. Октября, 144/3.

ВВЕДЕНИЕ.

ГЛАВА 1 МЕТОДЫ ПРОГНОЗИРОВАНИЯ ОБЪЕМОВ ЗАГРЯЗНЕННЫХ ПОЧВОГРУНТОВ ПРИ АВАРИЙНОМ ИСТЕЧЕНИИ НЕФТЕПРОДУКТОВ ИЗ ПОДЗЕМНЫХ ПРОДУКТОПРОВОДОВ.

1.1 Аварийность на трубопроводном транспорте нефти и нефтепродуктов.

1.1.1 Описание аварий на нефтепродуктопроводах Российской Федерации.

1.2 Теоретические и практические результаты изучения процесса фильтрации жидкости в пористой среде.

1.2.1 Уравнение неразрывности.

1.2.2 Упругий режим фильтрации.

1.2.4 Капиллярные силы.

1.2.5 Суффозионная устойчивость несвязанных грунтов.

1.3 Методы определения объемов загрязненных почвогрунтов.

1.4 Выводы.

ГЛАВА 2. РАЗРАБОТКА МЕТОДОЛОГИЧЕСКОЙ ОСНОВЫ ПРОГНОЗИРОВАНИЯ ОБЪЕМОВ ЗАГРЯЗНЕНИЯ ПОЧВОГРУНТОВ ПРИ АВАРИЯХ НА ПОДЗЕМНЫХ ПРОДУКТОПРОВОДАХ.

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

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

2.4 Выводы.

ГЛАВА 3. МЕТОД ОЦЕНКИ ПОСЛЕДСТВИЙ ЗАГРЯЗНЕНИЯ ПОЧВОГРУНТОВ ПРИ РАЗРУШЕНИИ ПОДЗЕМНОГО ПРОДУКТОПРОВОДА.

3.1 Характеристика типичных почвогрунтов Оренбургской области.

3.1.1 Описание почвогрунтов, их распространение по площади и глубине.

3.1.2. Фильтрационные характеристики грунта со структурой, нарушенной во время строительства продуктопровода.

3.2. Модель распространения углеводородов в толще грунта из поврежденного подземного трубопровода.

3.2.1. Математическая модель проникновения углеводородов в грунтах различной плотности и гранулометрического состава.

3.2.2 Укрупненный алгоритм и программный комплекс для реализации методики оценки последствий загрязнения почвогрунтов.

3.3. Результаты вычислительных экспериментов с применением полученных характеристик.

3.4 Выводы.

ГЛАВА 4. АПРОБАЦИЯ И ВНЕДРЕНИЕ МЕТОДИКИ ПРОГНОЗИРОВАНИЯ ОБЪЕМОВ ЗАГРЯЗНЕННЫХ ПОЧВОГРУНТОВ ПРИ АВАРИЙНОМ ИСТЕЧЕНИИ НЕФТЕПРОДУКТОВ ИЗ ПОДЗЕМНЫХ

ПРОДУКТОПРОВОДОВ.

4.1. Оценка последствий разрушения сварного шва конденсатопровода Оренбург-Салават-Уфа на 57 км.

4.2. Расист зон действия поражающих факторов аварий на примере проектных решений при обустройстве Филипповской залежи Оренбургского НГКМ.

При эксплуатации промысловых и магистральных продуктопроводов одной из актуальных проблем является обеспечение безопасной эксплуатации линейной части, которая решается, главным образом, за счет своевременного проведения восстановительных и ремонтных работ на основе данных диагностики. Однако проводимые мероприятия не могут полностью исключить возможность возникновения аварийных ситуаций. Тяжесть последствий аварий зависит от оперативности принятия мер по локализации и ликвидации разлившихся нефтепродуктов, а также готовности эксплуатирующей организации к таким авариям. Для этих целей осуществляется прогноз масштабов возможных аварийных ситуаций, их последствий для социальной сферы и окружающей среды.

Существующие модели и методики позволяют оценить последствия от уже произошедших аварий или, применительно к проливам нефтепродуктов, дать примерные характеристики геометрии, площади разлива и объемов загрязненных грунтов. В настоящее время отсутствует доступный инженерный метод моделирования разлива жидких углеводородов при истечении нефтепродуктов из подземных трубопроводов через дефектные отверстия максимальными линейными размерами до 0,75 от диаметра трубопровода. При нарушении целостности заглубленного трубопровода существенно затруднена количественная оценка ареала распространения нефтепродукта в грунте, поскольку данный вид аварий характеризуется продолжительностью во времени от момента начала утечки до момента ее обнаружения и остановки транспортировки.

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

Цель работы.

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

Основные задачи работы:

-сбор и систематизация характеристик аварий, анализ существующих методов и подходов к оценке объемов загрязненных почвогрунтов и методов расчета массового расхода жидких углеводородов из типичных дефектных отверстий в теле трубопроводов при аварийном истечении;

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

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

- проверка адекватности разработанного метода на примере реально произошедших аварий.

Методы решения поставленных задач.

Применялись методы аналитической оценки движения жидких углеводородов в пористой среде, скорости их фильтрации и действительной скорости движения жидкости в фунте; линейный закон Дарси; методы определения фильтрационных параметров пористых сред — коэффициентов фильтрации и проницаемости; методы построения объективных зависимостей между емкостными и фильтрационными параметрами пористых сред (между коэффициентами пористости и проницаемости); математические методы моделирования процессов фильтрации и процессов формирования опасных факторов техногенных аварий, а также численное решение комплексов математических уравнений и их практическая проверка на основе материалов реализованных аварий.

Научная новизна работы:

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

- определены условия перехода подземного истечения на стадию поверхностного истечения;

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

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

На защиту выносятся следующие основные положения:

- усовершенствованная теоретическая модель подземного фильтрационного движения жидких нефтепродуктов в месте прокладки продуктопроводов;

- модель подземного истечения жидких углеводородов из продуктопроводов, разработанная с учетом закономерностей строения и особенностей грунтов в районе их прокладки;

- метод оценки последствий аварийных истечений нефтепродуктов через дефектные отверстия максимальными линейными размерами до 0,75 от диаметра трубы;

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

Практическая значимость результатов работы.

Возможность использования метода при оценке рисков объектов трубопроводного транспорта нефтепродуктов и адекватном прогнозировании потенциального ущерба, нанесенного социальной сфере и окружающей природной среде при проектировании объектов нефтегазового комплекса. Метод опробован при разработке деклараций промышленной безопасности к проектам обустройства Филипповской и Среднекаменноугольной залежей Оренбургского нефтегазоконденсатного месторождения ООО «Газпром добыча Оренбург».

Результаты оценки последствий распространения нефтепродуктов в толще грунта позволили уточнить масштабы возможных. загрязнений почвогрунтов, разработать эффективные Планы локализации и ликвидации аварийных разливов нефти (ПЛАРН), обосновать состав сил и средств для проведения работ по локализации и ликвидации аварийных разливов нефти и нефтепродуктов для нефтегазодобывающих предприятий.

Апробация работы.

Основные научные положения и практические результаты работы неоднократно доложены, обсуждены, одобрены и рекомендованы к применению на международных и российских научно- технических конференциях и симпозиумах, включая: VI Всероссийскую научно-техническую конференцию «Актуальные проблемы развития нефтегазового комплекса России», РГУ имени И.М. Губкина (Москва, 2008); VIII Всероссийскую научно-техническую конференцию, посвященную 80-летию Российского государственного университета нефти и газа имени И.М. Губкина «Актуальные проблемы развития нефтегазового комплекса России» (Москва, 2010); III и IV научно-технические конференции с международным участием «Основные проблемы поиска, освоения и обустройства нефтегазовых месторождений и пути их решения» (Оренбург, 2009, 2010).

Публикации.

По теме диссертации опубликовано 14 печатных работ, в т.ч. 6 статей в изданиях, рекомендованных ВАК Министерства образования и науки РФ для опубликования основных результатов диссертационных работ.

За помощь и консультации при подготовке диссертационной работы и помощь в работе автор выражает глубокую благодарность коллективу ООО «ВолгоУралНИПИгаз»: заведующему сектором охраны водных, почвенных ресурсов и геологической среды Донецковой Анастасии Акимовне; заведующему сектором программных средств Рахману Герасиму Самуиловичу; заведующему сектором анализа промышленных рисков Даньшиной Ольге Анатольевне.

ОСНОВНЫЕ ВЫВОДЫ И РЕЗУЛЬТАТЫ РАБОТЫ

1. Выявлено, что наибольшая доля аварий, связанных с проливами нефтепродуктов, происходит через дефектные отверстия максимальными линейными размерами до 0,75 от диаметра продуктопровода. Проведен анализ существующих методов оценки объемов загрязненных почвогрунтов при аварийных разливах нефтепродуктов.

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

3. Разработан метод оценки загрязнения почвогрунтов и возможных последствий с точки зрения пожарной и промышленной безопасности при разрушении подземного продуктопровода, построен алгоритм и разработано программное средство для реализации метода фильтрационного движения нефти и нефтепродуктов в толще грунта, выполнены вычислительные эксперименты в математической среде МаШСас!.

4. Проведены апробация программного средства и сопоставление расчетных результатов по предложенному методу оценки геометрических параметров области загрязнения с таковыми реализованной аварии, произошедшей на IV нитке конденсатопровода Оренбург - Салават - Уфа (расчетный объем загрязненного грунта по предложенному методу составил около 6 тыс. м3, по о фактическим данным - около 7,4 тыс. м , разница составила около 19 %). Метод опробован для получения исходных данных при оценке зон действия опасных факторов аварийных сценариев (теплового воздействия от пожаров пролива, барического воздействия от взрывов топливно-воздушных смесей и токсического загрязнения окружающего пространства сероводородом) на выкидных трубопроводах системы сбора продукции в проектах обустройства Филипповской и Среднекаменноугольной залежей Оренбургского НГКМ.

1. Отчет о деятельности Федеральной службы по экологическому, технологическому и атомному надзору за 2009 год Федеральная служба по экологическому, технологическому и атомному надзору. - М.: 2010. - 460 с.

2. Отчет о деятельности Федеральной службы по экологическому, технологическому и атомному надзору за 2003 год Федеральная служба по экологическому, технологическому и атомному надзору. - М.: 2004. - 340 с.

3. Отчет о деятельности Федеральной службы по экологическому, технологическому и атомному надзору за 2004 год Федеральная служба по экологическому, технологическому и атомному надзору. - М.: 2005. - 344 с.

4. Отчет о деятельности Федеральной службы по экологическому, технологическому и атомному надзору за 2005 год Федеральная служба по экологическому, технологическому и атомному надзору. - М.: 2006. - 354 с.

5. Отчет о деятельности Федеральной службы по экологическому, технологическому и атомному надзору за 2006 год Федеральная служба по экологическому, технологическому и атомному надзору. - М.: 2007. - 508 с.

6. Отчет о деятельности Федеральной службы по экологическому, технологическому и атомному надзору за 2007 год Федеральная служба по экологическому, технологическому и атомному надзору. - М.: 2008. - 548 с.

7. Отчет о деятельности Федеральной службы по экологическому, технологическому и атомному надзору за 2008 год Федеральная служба по экологическому, технологическому и атомному надзору. - М.: 2009. - 447 с.

8. CONservation of Clean Air and Water in Europe (CONCAWE) Европейская ассоциация нефтяных компаний (http://www.concawe.beV

9. Office of Pipeline Safety under the Pipeline and Hazardous Materials Safety Administration (OPS PHMSA) Бюро трубопроводной безопасности Управления по безопасности трубопроводов и опасным веществам Министерства транспорта США (http://www.phmsa.dot.gov).

10. Прусенко Б.Е., Мартынюк В.Ф. Анализ аварий и несчастных случаев на трубопроводном транспорте. Учебное пособие. М.: ООО «Анализ опасностей», 2003.-351 с.

11. Басниев К. С., Кочина Н. И., Максимов М. В. Подземная гидромеханика. -М.: Недра, 1993.-416 с.

12. Кременецкий H.H., Штеренлихт Д.В., Алышев В.М., Яковлева JI.B. Гидравлика. Учебник. М.: Энергия, 1980. - 384 с.

13. Полубаринова-Кочина П. Я. Теория движения грунтовых вод. М.: Наука, 1952.-677 с.

14. Ентов В.М. Теория фильтрации // Соросовский образовательный журнал. -1998.-№2.-С. 121-128.

15. Максимов В.М. Справочное руководство гидрогеолога. Справочник. Л.: Недра, 1967.-592 с.

16. Баренблатт Г.И., Ентов В.М., Рыжик В.М. Движение жидкостей и газов в природных пластах. М.: Недра, 1984.-208 с.

17. Дерягин Б.В., Чураев Н.В., Мюллер В.М. Поверхностные силы. М.: Наука,1985.-398 с.

18. Хейфец Л. И., Неймарк А. В. Многофазные процессы в пористых средах. -М.: Химия, 1982.-320 с.

19. Адамсон А. Физическая химия поверхностей. М.: Мир, 1979. 568 с.

20. Жуковский Н.Е. Гидродинамика. М.: Государственное издательство технико-технической литературы, 1949. - 764 с.

21. Русанов А.И., Гудрич Ф.Ч. Современная теория капиллярности. Л.: Химия, 1980.- 344 с.

22. Роулинсон Дж., Уидом Б. Молекулярная теория капиллярности. М.: Мир,1986.-376 с.

23. Кондратьев В. Н. Фильтрация и механическая суффозия в несвязных грунтах. Симферополь, Крымиздат, 1958. - 75 с.

24. Коста Ж., Санглера Г. Механика грунтов. Практический курс (перевод с французского). М.: Стройиздат, 1981.- 451 с.

25. Постановление Правительства РФ «О неотложных мерах по предупреждению и ликвидации аварийных разливов нефти и нефтепродуктов» № 613 от 21.08.2000 г. 7 с.

26. Методическое руководство по оценке степени риска аварий на магистральных нефтепроводах. Серия 27. Выпуск 1. М.: Государственное предприятие Научно-технический центр по безопасности в промышленности Госгортехнадзора России, 2000. 96 с.

27. Бракоренко H.H., Емельянова Т.Я. Типизация грунтовых толщ по составу и проницаемости в связи с прогнозом загрязнения их нефтепродуктами (на примере г. Томска) // Вестник Томского государственного университета. 2007. — № 302. -С. 215-218.

28. Середин В.В. Санация территорий, загрязненных нефтью и нефтепродуктами // Геоэкология. 2000. - № 6. - С. 525-540.

29. Сергеев Е.М. Грунтоведение. М.: МГУ, 1983. - 392 с.

30. Мироненко В.А., Румынии В.Г. Проблемы гидрогеоэкологии: В 3 т. М.: МГУ, 2002. - 732 с.

31. Котляровский В.А., Ларионов В.И., Сущев С.П. Энциклопедия безопасности: строительство, промышленность, экология: в 3 т. Т. 1 . М.: Наука, 2005. - 696 с.

32. НПБ 105-03 Определение категорий помещений, зданий и наружных установок по взрывопожарной и пожарной опасности, утверждены приказом № 314 МЧС РФ от 18.06.2003. 30 с.

33. РД 153-39.4-060-00 Методика расчета ущерба от криминальных врезок в нефтепродуктопроводы, утверждены Приказом Минэнерго РФ № 167 от 06.06.2001.-30 с.

34. Климентьев А.И., Чибилев A.A., Блохин Е.В., Грошев И.В. Красная книга почв Оренбургской области. Екатеринбург, УрО РАН, 2001. 295 с.

35. Балыхин Г.А., Нзамуйе И. Взаимодействие гидромоторной струи с горным массивом // Тр.Университета дружбы народов. Вып. 155. - 1992. - С. 208-114.

36. Топчиев А. В., Ведерников В. И., Горные машины. Справочник. М.: 1960. -347 с.

37. Шорохов С. М., Разработка россыпных месторождений и основы проектирования, М., 1963;

38. Нурок Г. А., Гидромеханизация открытых разработок, М., 1970.

39. Шкундин Б.М. Оборудование гидромеханизации земляных работ. М.: Энергия, 1970.- 240 с.

40. Никонов Г.П., Кузьмич И.А., Гольдин Ю.А. Разрушение горных пород струями воды высокого давления. М.: Недра, 1986. - 142 с.

41. Павчич М. П., Балыков Б. И. Методы определения коэффициента фильтрации грунтов. М.: Энергия, 1976. - 115 с.

42. Веригин H.H. Гидродинамические и физико-химические свойства горных пород. М.: Недра, 1977.-271 с.

43. Маскет М. Течение однородных жидкостей в пористой среде. М.: Институт компьютерных исследований, 2004. - 640 с.

44. Бэр Я., Заславский Д., Ирмей С. Физико-математические основы фильтрации воды. М.: Мир, 1971.-452 с.

45. Терцаги К., Пек Р. Механика грунтов в инженерной практике / Пер. с английского. М.: Госстройиздат, 1958.-608 с.

46. Чаповский Е. Г. Лабораторные работы по грунтоведению и механике грунтов. М.: Недра, 1966.- 136 с.

47. Ломтадзе В. Д. Методы лабораторных исследований механических свойств горных пород. Л.: Недра, 1972. - 171 с.

48. Энгельгардт В. Поровое пространство осадочных пород / Перевод с немецкого. М.: Недра, 1964.-232 с.

49. П 12-83 Рекомендации по методике лабораторных испытаний грунтов на водопроницаемость и суффозионную устойчивость, утверждены ВНИИГ имени Б. Е. Веденеева решением № 12 от 24.03 1983. 55 с.

50. ГОСТ 12536-79 Грунты. Методы лабораторного определения гранулометрического (зернового) и микроагрегатного состава.

51. ГОСТ 23278-78 Грунты. Методы полевых испытаний проницаемости.

52. ГОСТ 25100-95 Грунты. Классификация.

53. ГОСТ 25584-90 Грунты. Методы лабораторного определения коэффициента фильтрации.

54. ГОСТ 23001-90 Грунты. Методы лабораторных определений плотности и влажности.

55. СНиП 2.05.06-85* Магистральные трубопроводы, утвержденным постановлением Минстроя России № 18-78 от 10.11.1996. 81 с.

56. СНиП 11-02-96 Инженерные изыскания для строительства. Основные положения, введен в действие постановлением Минстроя России № 18-77 от 29.10.1996.-85 с.

57. СП 11-105-97 Инженерно-геологические изыскания для строительства, одобрено письмом Госстроя России № 9-4/116 от 14.10.1997. 78 с.

58. Кульчицкий Л.И., Усьяров О.Г. Физико-химические основы формирования свойств глинистых пород. М., Недра, 1981.

59. Ломтадзе A.M., Ломтадзе В.Д. Руководство по инженерно-геологическому изменению глинистых пород. М„ Энергия, 1970.

60. Бондарик Г.К., Царева A.M., Пономарев В.В. Текстура и деформация глинистых пород. М., Недра, 1975.

61. Гуменский Б.М. Основы физикохимии глинистых грунтов и их использование в строительстве. Л.-М., Стройиздат, 1965. 421 с.

62. Деклу Ж. Метод конечных элементов. М.: Мир, 1976. - 325 с.

63. Бенерджи П., Баттерфилд Р. Методы граничных элементов в прикладных науках. М.: Мир, 1984. - 653 с.

64. Зенкевич О. Метод конечных элементов в технике. М.: Мир, 1975. - 348 с.

65. Галлагер Р. Метод конечных элементов. Основы. М.: Мир, 1984 - 415 с.

66. Коннор Дж., Бреббиа К. Метод конечных элементов в механике жидкости. -JI.: Судостроение, 1979. 543 с.

67. Оден Дж. Конечные элементы в нелинейной механике сплошных сред. М.: Мир. 1976.-612 с.

68. ПБ 09-540-03 Общие правила взрывобезопасности для взрывопожароопасных химических, нефтехимических и нефтеперерабатывающих производств, утверждены Постановлением Госгортехнадзора России № 29 от 05.05.2003. 14 с.

69. ГОСТ 12.1. 004-91. ССБТ. Пожарная безопасность, (с изменением № 1, утв. в октябре 1993 г.)

70. ГОСТ Р 12.3.047-98. Пожарная безопасность технологических процессов. Общие требования. Методы контроля.

71. ГОСТ 12.1.010-76* ССБТ. Взрывобезопасность. Общие требования, (с изменениями № 1, утв. 02.1983 г.)

72. РД 03-409-01 Методика оценки последствий аварийных взрывов топливно-воздушных смесей, утверждены Постановлением Госгортехнадзора России № 25 от 26.06.2001.-21 с.

73. РД 03-496-02 Методические рекомендации по оценке ущерба от аварий на опасных производственных объектах, утверждены Постановлением Госгортехнадзора России № 63 от 29.10.2002. 18 с.

74. Бесчастнов М.В. Промышленные взрывы. Оценка и предупреждение. М.: Химия, 1991.-432 с.

75. В. Маршалл. Основные опасности химических производств. М.: Мир, 1989. -551 с.

76. Абрамов С.А., Гнездилова Г.Г., Капустина E.H., Селюн М.И. Задачи по программированию. М.: Наука, 1988.-596 с.

77. Бентли Б. Жемчужины творчества программистов. М.: Радио и связь, 1990. -269 с.

78. Кормен Т.Х., Лейзерсон Ч.И., Ривест Р.Л., Штайн К. Алгоритмы: Построение и анализ, 2-е издание. М.: Вильяме, 2005. - 1296 с.

79. Вирт Н. Алгоритмы + структуры данных = программы. . М.: Мир, 1985. -406 с.

80. Вирт Н. Алгоритмы и структуры данных. М.: Мир, 1989. - 360 с.

81. Голицына О.Л., Попов И.И. Основы алгоритмизации и программирования. -М.: Форум, 2008.-432 с.

82. Грин Д., Кнут Д. Математические методы анализа алгоритмов. М.: Мир, 1982.- 120 с.

83. Заварыкин В.М., Житомирский В.Г., Лапчик М.П. Основы информатики и вычислительной техники. М.: Просвещение, 1988. - 176 с.

84. Могилев A.B., Пак Н.И., Хеннер Е.К. Информатика: Учеб. пособие для студ. педвузов. М.: ACADEMIA, 1999. - 848 с.

85. Лингер Р., Миллс X., Уитт Б. Теория и практика структурного программирования. М.: Мир, 1982.-406 с.

86. Лэнгсам И., Огенстайн М., Тененбаум А. Структуры данных для персональных ЭВМ. М.: Мир, 1989. - 568 с.

87. Острейковский В.А. Информатика: Учебник для вузов. М.: Высш. шк., 2000.-511 с.

88. Сегерлнид Л. Применение метода конечных элементов. М.: Мир, 1979. -571 с.

89. Норрн Д., де Фриз Ж. Введение в метод конечных элементов. М.: Мир, 1981.-481 с.

90. MathCAD 6.0 Plus. Руководство пользователя / перевод с английского. М.: Филинъ. 1996. -712 с.

91. Дьяконов В.П. Справочник по MathCAD 6.0 Plus Pro. М.: СК Пресс, 1997. -336 с.

92. Дьяконов В.П., Абраменкова И.В. MathCAD 7.0 в математике, физике и в Internet. М.: Нолидж, 1998.- 352 с.

93. Очков В. Ф. MathCAD 8 PRO для студентов и инженеров. М.: Компьютер Пресс, 1999. - 523 с.97.0чков В. Ф. MathCAD 14 для студентов и инженеров. СПб.: БХВ -Петербург, 2009. -512 с.

94. Макаров Е. Инженерные расчеты в Mathcad 15. Учебный курс. СПб.: ПИТЕР, 2011.-400 с.

95. СТО 2-2.3-351-2009 Методические указания по проведению анализа риска для опасных производственных объектов газотранспортных предприятий ОАО «Газпром». М.: ООО «Газпром экспо», 2009. - 377 с.

96. СТО Газпром 2-2.3-400-2009 Методика анализа риска для опасных производственных объектов газодобывающих предприятий ОАО «Газпром». -М.: ООО «Газпром экспо», 2009. 300 с.

97. Руководство по оценке пожарного риска для промышленных предприятий. -М.: ФГУ ВНИИПО МЧС России, 2006. 45 с.

98. ГОСТ Р 12.3.047-98 Пожарная безопасность технологических процессов Приложение Г (Рекомендуемое).

99. РД 03-26-2007 Методические указания по оценке последствий аварийных выбросов опасных веществ, утверждены Постановлением Госгортехнадзора России № 261 от 29.10.2007. 34 с.