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

ГОРОХОВ Николай Леонидович

РАЗРАБОТКА МАТЕМАТИЧЕСКОЙ МОДЕЛИ

СЕЙСМИЧЕСКОГО ВОЗДЕЙСТВИЯ НА ПОДЗЕМНЫЕ

ГАЗОПРОВОДЫ (НА ПРИМЕРЕ МАГИСТРАЛЬНОГО

ГАЗОПРОВОДА «ЛЕНИНГРАД-ВЫБОРГ-ГОСГРАНИЦА»)

25.00.20 –

Специальность Геомеханика, разрушение горных пород, рудничная аэрогазодинамика и горная теплофизика

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

Санкт-Петербург – 2013

Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования "Национальный минеральносырьевой университет "Горный".

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

доктор технических наук, профессор Господариков Александр Петрович

Официальные оппоненты:

Эквист Борис Владимирович доктор технических наук, профессор, ФГБОУ ВПО "Московский государственный горный университет", кафедра взрывного дела, профессор Вильчинский Владислав Борисович кандидат технических наук, ООО "Институт Гипроникель", горная лаборатория, заведующий

Ведущая организация – ФГБОУ "Институт проблем комплексного освоения недр Российской академии наук"

Защита диссертации состоится 28 июня 2013 г. в 11 ч на заседании диссертационного совета Д 212.224.06 при Национальном минерально-сырьевом университете "Горный" по адресу: 199106, Санкт-Петербург, 21-я линия, дом 2, ауд. 1166.

С диссертацией можно ознакомиться в библиотеке Национального минерально-сырьевого университета "Горный".

Автореферат разослан 28 мая 2013 г.

УЧЕНЫЙ СЕКРЕТАРЬ СИДОРОВ

диссертационного совета Дмитрий Владимирович.т.н., профессор Э.И. Богуславский

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

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

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

Сейсмическое воздействие взрывных работ на охраняемые объекты исследовалось в работах М.А. Садовского, В.В. Адушкина, Б.Н. Кутузова, С.В. Медведева, Б.В. Эквиста, В.Ф. Богацкого, П.С.

Миронова, Я.И. Цейтлин и др.

Воздействие ударных волн в различных грунтовых средах на стальные трубопроводы исследовалось в работах К.Е. Ращепкина, А.Г. Гумерова, В.С. Силина, А.Г. Горшкова, Р.Г. Якупова, Е.Г. Янютина, В.Г. Баженова, А.В. Кочеткова, Т.Р. Рашидова и др.

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

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

Цель диссертационной работы. Разработка математической модели сейсмического воздействия на подземные газопроводы для оценки сейсмобезопасных параметров буровзрывных работ (БВР), проводимых в зоне залегания последних.

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

Основные задачи исследований:

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

• разработка численных методик и алгоритмов, необходимых для прямого интегрирования разрешающих уравнений математической модели;

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

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

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

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

Научная новизна диссертационной работы заключается в следующем:

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

• получена система разрешающих уравнений, для которой построена схема метода конечных элементов, описывающая pmlслои для задач сейсмодинамики, позволяющая исключить «фиктивные» отражения волн от условного контура расчетной модели и исследовать волновые процессы в грунтовых средах;

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

Основные защищаемые положения:

1. Исследование процессов взаимодействия сейсмовзрывных волн с подземными газопроводами необходимо осуществлять на основе математической модели совместного динамического деформирования системы грунт – трубопровод с учетом нелинейных эффектов их контактного взаимодействия, необратимых объемных и сдвиговых деформаций, объемной вязкости грунта.

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

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

Практическая значимость работы:

• разработан комплекс вычислительных программ на языке Matlab, для решения задач взаимодействия сейсмических волн в грунтах с подземными газопроводами методом конечных элементов в плоской постановке;

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

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

Апробация работы. Содержание и основные положения диссертационной работы докладывались и обсуждались на симпозиуме «Неделя горняка-2011» (Москва, 2011 г.), на международной молодежной научной конференции во Фрайбергской горной академии (Германия, г. Фрайберг, 2011 г.), на ежегодной научной конференции молодых ученых и студентов «Полезные ископаемые России и их освоение» (Санкт-Петербург, 2010 г.), на заседаниях кафедры безопасности производств и разрушения горных пород и кафедры взрывного дела.

Личный вклад автора:

• построение математической модели взаимодействия сейсмовзрывных волн с подземным газопроводом;

• разработка численных методик и алгоритмов прямого интегрирования разрешающих уравнений;

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

• разработка метода расчета сейсмобезопасной массы заряда при ведении взрывных работ вблизи действующих газопроводов;

• построение математической модели и численной схемы метода конечных элементов для pml-слоев в задачах сейсмодинамики.

Реализация результатов работы. Разработанная методика определения сейсмобезопасной массы одновременно взрываемых зарядов применяется для расчета параметров БВР при проходке траншеи под «Североевропейский газопровод» вблизи действующего Магистрального газопровода «Ленинград-Выборг-Госграница».

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

Объем и структура работы. Диссертация состоит из введения, четырех глав, заключения и библиографического списка, изложенных на 142 страницах машинописного текста, содержит 53 рисунка, 19 таблиц и список литературы из 140 наименований.

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

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

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

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

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

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

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

В четвертой главе представлены результаты инструментальных измерений параметров сейсмовзрывных волн при проведении полигонных экспериментальных взрывов, которые позволили определить параметры вмещающих грунтов, находящихся в зоне залегания действующего газопровода; описан метод определения сейсмобезопасной массы заряда при ведении взрывных работ вблизи подземного газопровода, использующий результаты численного моделирования; определена сейсмобезопасная масса одновременно взрываемых зарядов при проведении буровзрывных работ вблизи действующего магистрального газопровода «Ленинград-ВыборгГосграница»; с учетом полученной величины сейсмобезопасной массы заряда был произведен расчет параметров буровзрывных работ для проведения траншеи.

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

Основные результаты исследований отражены в следующих защищаемых положениях:

1. Исследование процессов взаимодействия сейсмовзрывных волн с подземными газопроводами необходимо осуществлять на основе математической модели совместного динамического деформирования системы грунт – трубопровод с учетом нелинейных эффектов их контактного взаимодействия, необратимых объемных и сдвиговых деформаций, объемной вязкости грунта.

Расчетная схема, применяемая для решения задачи моделирования взаимодействия сейсмовзрывных волн с подземными газопроводами выглядит в виде, представленном на рисунке 1.

Рисунок 1 – Расчетная схема задачи: x1, x2 – координатные оси; стрелками обозначено направление, вдоль которого происходит поглощение волн в pml-слоях Разработанная математическая модель, описывающая расчетную схему, включает в себя систему разрешающих уравнений и подмоделей:

• вариационную постановку задачи динамического деформирования сплошной среды;

• вариационную постановку задачи динамического деформирования трубных оболочек;

• модель динамического деформирования грунта;

• модель поглощающих граничных условий (ПГУ), pmlслои;

• модель контакта системы «грунт-трубопровод».

Из различных формулировок задач исследования динамического деформирования механических систем с помощью МКЭ, как показала практика, наиболее эффективным является использование вариационного принципа возможных перемещений Лагранжа в сочетании с принципом Даламбера.

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

где ij – компоненты тензора напряжений, Па; ij – компоненты тензора деформаций; u i – упругие перемещения частицы деформируемого тела по направлению координатной оси xi, м; FV i – массовые силы, Н; FS i – поверхностные силы, заданные на внешней поверхности S, окружающей объем пространства V, Н; – плотность деформируемого тела, кг/м3; t – время, с.

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

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

Для мягких грунтов приняты известные уравнения состояния академика С.С. Григоряна (рисунок 2).

Рисунок 2 – Диаграмма сжатия в модели академика С.С. Григоряна Диаграмма сжатия имеет начальный линейно-упругий участок (р< pe), при котором деформации грунта являются чисто упругими; pe и e – предельные значения давления р и плотности 7 на упругом участке деформирования среды;. При этом давление и плотность в диапазоне изменения давления от – pe до pe связаны линейной зависимостью вида:

где p 0 – начальное давление в грунте, Па; 0 – начальная плотность в грунте, кг/м3; С – числовая константа.

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

где 1, 2, 3 – объемные концентрации свободной пористости (защемленного воздуха), жидкости (воды) и твердого компонента, E1, E 2, E 3 – модули объемного сжатия компонентов, Па; 1, 2, 3.числовые коэффициенты.

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

где С1 и С2 - скорости звука, определяющие наклоны к оси первого и второго звеньев ломаной, соответственно; p00 = p c характеризует отношение длин участков ломаных; ( 00, p 00 ) – точка излома кривой разгрузки в координтаах (, p ) ; параметр c задает отношение С1 к C2 при = lim ; lim – предельное значение плотности 7, после которого кривая разгрузки совпадает с кривой нагрузки.

Сдвиговое упругопластическое деформирование грунтовой среды в допредельном состоянии описывается моделью линейно-упругой среды. В зоне упругопластического деформирования зависимость p() описывается схемой Прандтля-Рейса с условием пластичности Мизеса.

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

где S ij – компоненты девиатора тензора напряжений, Па; G – модуль сдвига, Па; ij – компоненты девиатора тензора деформаций;

ij – компоненты девиатора тензора скоростей деформаций;

S = ii – среднее напряжение, Па; = ii – средняя деформация; K – модуль всестороннего сжатия, Па; – коэффициент динамической вязкости, Па·с.

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

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

В результате численной реализации МКЭ расчетная область V разбивается на N конечных элементов, объемом Ve каждый. Внутри элементов искомые функции напряжений ij и перемещений u i аппроксимируются с помощью линейных интерполяционных полиномов, записанных в матрицах функций формы для напряжений [N G ] и перемещений [NU ] относительно значений напряжений {G e } и перемещений {U e } в узлах элементов. Далее, для каждого конечного элемента составляются дискретные аналоги исходных дифференциальных уравнений относительно узловых значений. После суммирования этих уравнений по всем элементам в сетке конечных элементов, уравнения динамики механической системы переводятся в матричные дифференциальные уравнения относительно глобальных векторов неизвестных {G}, {U }. В результате получаются конечно-элементные аналоги всех необходимых разрешающих уравнений математической модели совместного динамического деформирования грунта и трубопровода.

Для исключения отражения волн от условного контура используется вид поглощающих граничных условий – метод введения идеально согласованных слоев (pml) на фиктивных границах расчетной области. Pml-слои представляют собой специальные зоны поглощения, которые отлавливают и ослабляют приходящие в них волны, при этом не давая отражений в расчетную область. Определяющие уравнения, описывающие движение волн напряжений в области pml в плоской постановке, имеют вид:

где { } – вектор напряжений, Па; {u} – вектор перемещений, м; [D] – матрица упругих характеристик среды; µ – коэффициент демпфирования, Па·с; 0 – опорная частота, с-1; cs – скорость поперечной волны, м/с; точками сверху обозначены производные по времени.

Матрицы F1e, F2e, F1p и F2p имеют вид:

где f ie – профиль демпфирования для волн, распространяющихся перпендикулярно протяженной координате ~i ; f i p – профиль демпx фирования для волн, распространяющихся параллельно протяженной координате ~i.

Выражения f m, f c и f k определяются по формулам:

С помощью разработанной численной модели получены решения задачи распространения сейсмической волны в безграничной среде сначала без использования pml-слоев (рисунок 3), а затем с использованием последних (рисунок 4).

Рисунок 3 - Распространение сферической волны в безграничной среде без использования pml-слоев: ti – момент времени, c Рисунок 4 - Распространение сферической волны в безграничной среде с применением pml-слоев: ti – момент времени, c Как видно из результатов численного эксперимента, технология применения pml-слоев позволяет эффективно поглощать падающие волны. Для совместного решения уравнений динамики грунта и трубопровода необходимо обеспечить условия сопряжения контактирующих поверхностей между трубопроводом и окружающим его грунтом. В настоящей работе для этого используются контактные конечные элементы сопряжения конструкций (КЭСК), (рисунок 5).

Рисунок 5 Принципиальная схема КЭСК: А, В – узлы контектного При этом величина производной энергии контактного взаимодействия, определяющая собой нагрузку от контактных сил, включая силы предварительного сопряжения, используется при формировании глобального вектора сил. Одновременно, формирование этого вектора сопровождается необходимой модификацией в глобальной матрице жесткости, вызванной применением метода штрафных функций. Итерационный процесс решения нелинейной задачи организован относительно минимизации величины невязки (неравенства) вектора контактных перемещений в глобальной системе рассматриваемой деформируемой конструкции.

Для численного интегрирования уравнений динамики модели из конечных элементов в настоящей работе применяется метод прямого интегрирования по времени (схема Хаболта).

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

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

Для определения первого приближения к сейсмобезопасной массе заряда используются эмпирические формулы, приведенные в совместной работе профессора Г.П. Парамонова и доцентов кафедры взрывного дела Горного университета В.А. Артемова, В.Н. Ковалевского и Ю.И. Виноградова, на основе которых были определены нижняя Qmin и верхняя Qmax границы ориентировочной массы взрываемого заряда:

где R = 18 м – расстояние от линии зарядов до газопровода, м;

1=0,7вр – расчетное допустимое напряжение в стальной стенке трубопровода с учетом длительного срока эксплуатации газопровода, МПа; вр = 520 МПа – временное сопротивление стали разрыву;

0 = 251 МПа – кольцевые напряжения в стенке трубопровода от внутреннего давления P = 5 МПа; = 1, kв = 1, kр = 1, kз = 1, kн = 1, kч = 1, kt = 1 – числовые коэффициенты, зависящие от условий веления взрывных работ, группы разрыхляемых пород и условий эксплуатации трубопровода.

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

где Q – масса одновременно взрываемых зарядов ВВ, кг; Кc – коэффициент сейсмичности, отражающий свойства грунта в месте расположения охраняемого объекта, определяемый экспериментально; R – расстояние от места регистрации до очага взрыва, м; n - показатель степени затухания колебаний, зависящий от свойств пород по профилю распространения сейсмовзрывных волн.

Для определения эмпирических коэффициентов Кc и n использованы данные экспериментов доцента В.А. Артемова по измерениям скорости смещения в сейсмовзрывной волне при проведении опытных взрывов. В результате обработки данных экспериментов получены значения Кc=670, n = 2,2. Графики (экспериментальный и расчетный) изменения значения массовых скоростей смещений в сейсмовзрывной волне от относительного расстояния в диаметрах шпуров и скважин представлены на рисунке 6.

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

где KР – экспериментальный коэффициент; Qл – линейная масса скважинного заряда, кг/м; l – длина скважинного заряда, м.

Рисунок 6 – Зависимость изменения скорости смещения в сейсмовзрывной волне от относительного расстояния в диаметрах заряда Расчетная область разбивается сеткой конечных элементов и включает в себя подобласти оболочки, скального грунта, песчаной засыпки, а также подобласти pml-слоев и объемного генератора волн, в котором генерируется сейсмовзрывная нагрузка с заданными параметрами. Между областями засыпки и скального грунта введены переходные области, в которых деформационные параметры плавно меняются от области скального грунта к мягкому. Последнее необходимо для достижения устойчивости и сходимости решения в местах соединения областей с разными деформационными характеристиками (рисунок 7).

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

Параметры падающей сейсмовзрывной волны, в этом случае составляют: амплитуда скорости смещения v = 17,4 cм/c; период колебаний T=15 мс.

На рисунках 8 - 14 представлены результаты моделирования воздействия сейсмовзрывной волны на газопровод при взрыве заряда массой 40 кг.

Максимальные суммарные напряжения в трубопроводе были достигнуты в момент времени t=17,6 мc и составили max = 356 МПа.

1 – трубопровод; 2 – песчаная засыпка; 3 – переходные области;

4 – скальный грунт;5 – объемный генератор волн; 6 – pml-слои Рисунок 9 Изополя напряжений в грунте и трубопроводе, t=10 мс: 11 – нормальные наРисунок 10 Изополя напряжений в грунте и трубопроводе, t=15 мс: 11 – нормальные напряпряжения вдоль оси x1, экв – эквивалентные напряжения по Мизесу Рисунок 11 Изополя напряжений в грунте и трубопроводе, t=20 мс: 11 – нормальные Рисунок 12 Изополя напряжений в грунте и трубопроводе, t=25 мс: 11 – нормальные напряжения вдоль оси x1, экв – эквивалентные напряжения по Мизесу напряжения вдоль оси x1, экв – эквивалентные напряжения по Мизесу Рисунок 13 Изополя напряжений в грунте и трубопроводе, t=30 мс: 11 – нормальные напряжения Увеличение напряжений в трубопроводе над уровнем начальных напряжений от действия внутреннего давления составило:

= max0 = 104 МПа.

Поскольку максимальные напряжения не превышают принятый уровень допустимых напряжений 1 = 364 МПа, то они являются безопасными для трубопровода и не могут разрушить его. Следовательно, заряды массой Q = 40 кг, приводящие к формированию сейсмовзрывной волны с амплитудой v = 17,4 cм/c на подходе к трубопроводу являются сейсмобезопасными.

Из результатов численного моделирования следует, что сейсмобезопасные массы зарядов, определяемые с помощью рекомендуемых приближенных формул (Q=14 кг и Q=18 кг), сильно занижены - трубопровод выдержит и больший по величине заряд.

ЗАКЛЮЧЕНИЕ

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

Основные результаты диссертационной работы позволяют сделать следующие выводы и рекомендации:

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

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

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

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

5. Предложен метод определения сейсмобезопасной массы заряда при ведении взрывных работ вблизи действующих газопроводов на основе разработанных алгоритмов и программ.

6. Рассчитана сейсмобезопасная масса одновременно взрываемых зарядов при проведении буровзрывных работ вблизи действующего Магистрального газопровода «Ленинград-ВыборгГосграница» (для расстояния 18 м от места взрыва до газопровода – Q=40кг при допустимой скорости смещения v = 17,4 см/с).

7. Из результатов численного моделирования следует, что сейсмобезопасные массы зарядов, определяемые с помощью нормативных приближенных формул (14 и 18 кг), сильно занижены - трубопровод выдержит и больший по величине заряд (40 кг).

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ:

1. Горохов Н.Л. Численное моделирование процессов воздействия сейсмовзрывных волн на подземные трубопроводы // Горный информационно-аналитический бюллетень. №3/2013.

М., 2013. – С. 211 – 217.

2. Горохов Н.Л., Господариков А.П. Оценка прочности газопровода, находящегося в грунте при ведении взрывных работ // Записки горного института. т. 195. С-Пб, 2012. – С. 89 – 94.

3. Господариков А.П., Горохов Н.Л. Динамический расчет трубопроводов на сейсмические воздействия // Записки горного института. т. 193. С-Пб, 2011. – С. 318 – 321.

4. Gorokhov N.L. The mathematical formulation and numerical implementation of dynamic problems of geomechanics using finite element method // Scientific reports on resource issues, vol 1. Frierberg, International University of resources, 2011. – p. 213 – 220.

Положительное решение № 2011124558, 2011 г. на Патент, «Способ рассредоточения заряда в скважине» / Т.М. Магомедов, Г.П. Парамонов, В.Н. Ковалевский, А.Е. Румянцев, Н.Л. Горохов, М.С. Шалаев.