CC BY
64
Худякова А.А. и др. Модель деятельного слоя многолетнемерзлых пород...
УДК 622.692.4
Модель деятельного слоя многолетнемерзлых пород с морозобойными трещинами и усовершенствование методики расчёта напряженно-деформированного состояния грунтовых оснований
А. А. Худякова, М.Г. Губайдуллин, А.В. Конюхов
Северный (Арктический) федеральный университет, Институт нефти и газа, кафедра транспорта и хранения нефти и газа, г. Архангельск
Аннотация. На основе предложенной новой модели деятельного слоя многолетнемёрзлых пород усовершенствована методика расчёта напряженно-деформированного состояния промерзающих-оттаивающих оснований. Она учитывает влияние морозобойных трещин, вертикальных и горизонтальных составляющих тепломассообмена, что позволяет повысить достоверность результатов расчетов, надёжность, продолжительность эксплуатации надземных трубопроводов в северных регионах и позволяет более широко применять поверхностные фундаменты.
Abstract. On the basis of the offered new model of active layer permafrost the technique of calculation of intense-deformed condition of the freezing through-thawing bases has been advanced. It considers influence of cracks with ice, vertical and horizontal components of heat and mass exchange that allows to increase reliability of results of calculations, reliability, duration of elevated pipelines operation in the northern regions and allows to apply the superficial bases more widely.
Ключевые слова: многолетнемерзлый, морозобойные трещины, модель деятельного слоя, расчёт напряженно-деформированного состояния грунтовых оснований, надземный трубопровод, фундамент поверхностного типа
Key words: permafrost, cracks with ice, model of active layer, calculation of intense-deformed condition of soil bases, elevated oil pipeline, base of superficial type
1. Постановка задачи
Темп нефтегазового освоения северных регионов России существенно зависит от сроков возведения, эксплуатации трубопроводных систем. Трубопроводное строительство в Заполярье, и в частности в Ненецком автономном округе (НАО), в условиях повсеместного распространения многолетнемёрзлых пород (ММП), избыточной обводненности, засолённости грунтов на побережье Баренцева моря, ранимости тундровых земель от техногенных воздействий характеризуется значительными затратами.
В НАО с учётом мерзлотно-грунтовых условий территории основная часть нефтегазопромысловых сооружений возводится на заглублённых фундаментах - бурозабивных и буроопускных сваях. Их применение часто приводит к сверхнормативным вертикальным и горизонтальным деформациям, осложнениям при эксплуатации сооружений (Винокуров и др., 2008), что обусловливает дополнительные затраты на ремонт и ликвидацию аварийных ситуаций. Фундаменты поверхностного типа (плитные) в НАО применяются в ограниченном количестве из-за недостаточной разработанности и адаптированности модели деятельного слоя, теории и методики расчетов оснований к условиям развития ММП. В связи с этим представляется актуальным усовершенствование методики расчёта напряженно-деформированного состояния (НДС) промерзающих-оттаивающих грунтовых оснований, устойчивости фундаментов и трубопроводов на ММП с учетом состояния и изменений геокриологической среды с морозобойными трещинами на контакте с фундаментами и нефтепроводами применительно к природно-климатическим условиям НАО.
2. Модель деятельного слоя ММП с морозобойными трещинами, учитывающая наличие горизонтальных и вертикальных составляющих тепломассообмена
Геокриологическими исследованиями С.Е. Гречищева, Л.В. Чистотинова, Р.П. Петровой, С.Е. Мельникова и др. установлено, что под воздействием фазовых превращений влаги в покровных отложениях при высоком положении уровня грунтовых вод (УГВ) формируются морозобойные трещины. Однако с позиций влияния на водный и тепловой режимы и физико-механические свойства деятельного слоя ММП морозобойные трещины до настоящего времени не изучались.
Ниже приведены результаты геокриологических изысканий, проведённых в НАО на трассе нефтепровода Харьяга - Южное Хыльчую, на НПС "Инзырей" и "Восточное Сарутаю". На указанных
810
Вестник МГТУ, том 13, №4/1, 2010 г.
стр.810-815
объектах в 2007-2008 гг. изучались: геометрические размеры морозобойных трещин, процессы формирования в деятельном слое блочных структур, изменения водно-тепловых свойств покровных отложений и подстилаемых ММП. В плане трещины имели четырёхугольную и пятиугольную формы со сторонами 1,6...4,4 м при ширине раскрытия на поверхности 0,05...0,35 м. Глубина смыкания трещин соответствовала границам максимального оттаивания грунтов (0,5.2,1 м) или положению верхней плоскости подстилаемых слитых ММП. Наличие грунтовых блоков, разобщённых между собой трещинами (рис. 1), является определяющим фактором для оценки исходного и прогнозного состояний термовлажностного режима, физико-механических свойств покровных отложений.
Предлагаемая модель деятельного слоя ММП позволяет прогнозировать изменения капиллярной (вертикальной) и гидрологической (горизонтальной) составляющих водного режима грунтов на осваиваемых участках. Так, например, по данным наблюдений на площадке, имеющей морозобойные трещины шириной 0,25 м, глубиной 1,2 м, протяжённостью сторон в плане 2,4 м, было установлено следующее. Объём капиллярной влаги в одном грунтовом блоке, имеющем в плане размеры 2,4x2,4 м, УГВ = 1,2 м и осреднённую по высоте объёмную влажность w = 0,22, к концу весеннего периода составлял 1,2 м3. Содержание воды в трещинах в этот период было наибольшим - 1,44 м3. В меженный период года объём воды в трещинах составлял 30 % от максимального, в осенний период - 70 % от значений весеннего паводка. Объем капиллярной влаги внутри блоков был практически неизменным в течение всего периода. Содержание влаги (льда) в трещинах зависит от размеров трещин и положения УГВ в годовом режиме, а в разобщённых блоках - от атмосферных осадков, положения УГВ и объема капиллярной влаги. Водный режим в блоках можно прогнозировать относительно УГВ экспоненциальными зависимостями.
Наличие в деятельном слое блоков, разобщённых между собой трещинами, существенно отражается на формировании вертикальных и горизонтальных составляющих тепломассопереноса в грунтовых разрезах, фазовых превращений влаги в трещинах и грунтовых блоках. Наблюдениями за промерзанием-оттаиванием деятельного слоя было установлено превышение значений осадок (на 25.35 %) при оттаивании грунтов в разобщённых трещинами блоках в сравнении с деформациями их пучения. Это вызвано одновременными процессами миграции влаги из блоков к верхней границе промерзания и горизонтального оттока влаги из блоков в промерзающие трещины.
Приведенные выше данные подтверждают, что при определении НДС грунтовых оснований с морозобойными трещинами необходимо учитывать наличие вертикальных и горизонтальных составляющих тепломассообмена как непосредственно в грунтовых блоках, так и в морозобойных трещинах. Результаты проведённых изысканий позволяют повысить достоверность оценок исходных данных, прогноза изменений геокриологического состояния среды на контакте с сооружениями, получить данные для расчётов НДС грунтовых оснований, устойчивости нефтепроводов и других объектов на Крайнем Севере.
Рис. 1. Модель структуры деятельного слоя ММП с морозобойными трещинами
1 - грунтовые блоки, разобщённые морозобойными трещинами;
2 - водонепроницаемый подстилающий слой, верхняя граница которого соответствует глубине промерзания-оттаивания грунтов;
3 - морозобойные трещины, заполненные водой или льдом
3. Совершенствование методики расчёта напряженно-деформированного состояния промерзающих-оттаивающих грунтовых оснований, устойчивости фундаментов и нефтепроводов на ММП
В расчётах НДС талых и мёрзлых грунтовых оснований под сооружениями наиболее часто используют теорию линейно-деформируемых сред, позволяющую определять перемещения Ut точек в упругой среде грунтового полупространства, ограниченного плоскостью (x, y), на которую действует сила P(x, y). Тогда для силы
P = F-8(x) S(y),
где S - функция Дирака; F - сила, действующая на фундамент; с учётом работы Е.Ф. Винокурова можно получить численные значения перемещений Ux, Uy, Uz..
Ниже приведены результаты сравнения расчётных значений деформаций талых и мёрзлых грунтов в основании плитного фундамента, выполненные по теории линейно-деформируемых сред, с фактическими осадками плитного фундамента. Использование плитных фундаментов на деформируемых
811
Худякова А.А. и др. Модель деятельного слоя многолетнемерзлых пород...
засоленных грунтах позволяет снизить стоимость строительства до 32 %, повысить надежность эксплуатации нефтепроводов за счет возможности обеспечения горизонтальных перемещений и компенсации вертикальных подвижек (Худякова, 2010).
Расчёты деформаций неоднородных по глубине грунтов проводились из условий решения плоской задачи взаимодействия грунтового полупространства с пластинкой (шириной 0,9 м, длиной 10 м), на которую действует сила 123 кН/м. Изменения физико-механических свойств геологической среды в течение года учитывались соответствующими корректировками упругих констант грунта. Осадка плитного фундамента в натурных условиях была определена в ходе наблюдений за его перемещениями в процессе бурения скважины № 21 на Восточно-Сарутаюском месторождении НАО. Фактическая осадка оттаивающих грунтов в основании поверхностного фундамента составила 0,28 м и в 7 раз превысила расчётное значение. Такое несоответствие результатов объясняется тем, что НДС деформируемых грунтов под подошвой мелко заглублённых или плитных фундаментов нельзя рассчитывать по теории, не учитывающей нелинейную зависимость между напряжениями и деформациями, явления консолидации и вязкости деформируемых при промерзании-оттаивании грунтов деятельного слоя.
Полученные результаты подтвердили необходимость применения для решения подобных задач вязкопластической теории. Для определения НДС сильно деформируемых грунтов деятельного слоя под фундаментами были использованы следующие уравнения Пыжина:
где J\ - первый инвариант тензора напряжений; J22 - второй инвариант девиатора напряжений; д - упругий модуль сдвига; у - постоянная вязкости; к - предел текучести при простом сдвиге; а - постоянная величина, характеризующая скорость объемного расширения грунта; др - символ Кронекера, др = 1 при i = j и др = 0 при i Ф j. Функция Ф(М) определяется по результатам испытаний грунтов на динамические воздействия.
Уравнение (2) удовлетворительно описывает скорость объёмных деформаций грунтового полупространства. Постулируя определённую зависимость между вторым инвариантом тензора скорости неупругой деформации и скоростью продвижения границы промерзания-оттаивания, используя решения теории линейно-деформируемой среды в качестве начальных условий, можно из уравнений (1) и (2) найти зависимость перемещений щ точек основания во времени, определить площадь фундамента, в целом оценить устойчивость проектируемого сооружения. Результаты расчетов по вязкопластической теории не соответствовали деформациям грунтов под буровой установкой. Фактические осадки грунтов превышали расчётные значения более чем в 2,5 раза.
Совершенствование методики расчёта включало использование полученных данных по формированию морозобойных трещин, водного и теплового режимов в деятельном слое ММП. По методике с использованием известных теоретических положений определён радиус-вектор силы P в точке (x, y, z) грунта, составлен и реализован методами конечных элементов алгоритм определения НДС в грунтовых блоках с учетом наличия в них подвижной границы фазовых превращений влаги.
Расчет теплового режима в грунтовом блоке включал решение с учетом изменения физикомеханических свойств основания двумерного уравнения теплопроводности. При двух переменных нестационарное уравнение теплопроводности можно представить в виде:
где T = T(r,t) - температура в точке с радиус-вектором r(x,y) в момент времени t; c(T) - коэффициент теплоемкости; к(Т) - коэффициент теплопроводности.
Поставленная задача была реализована численным методом с учетом переменных направлений тепловых Q и водных W потоков в грунтовом блоке, загруженном усилием q (рис. 2).
Для нахождения численного решения уравнения (3) применительно к грунтовому блоку использовался метод контрольных объемов. Температура, теплофизические параметры породы рассчитывались в центрах объемов и принимались за характеристики всего объема. Предлагаемый подход позволил получить наиболее достоверные значения распределения тепловых потоков и НДС внутри грунтового блока под сооружениями. Реализация поставленной задачи была осуществлена с использованием данных инженерно-геологических изысканий. Результаты расчёта осадки фундамента при оттаивании грунтов в основании имели близкие к натурным значениям и равнялись соответственно 0,22 и 0,28 м.
(1)
(2)
c iL — jL (к ■dT) + -i- (к —)
C dt ~ dx \K dx dy \K by)
dyj ,
(3)
812
Вестник МГТУ, том 13, №4/1, 2010 г.
стр.810-815
Рис. 2. Схема расположения нагрузки, тепловых и водных потоков в грунтовом блоке
Для выполнения на основе теории вязкопластичности многочисленных расчётов НДС на неслитых ММП, расположенных под плитным фундаментом площадью 2,25 м2, были использованы следующие исходные данные: глубина сезонного оттаивания пылеватых суглинков Нот = 2,0 м; диаметр нефтепровода изменялся в пределах от 213 до 500 мм; температура оттаивающего грунта t^ под подошвой фундамента изменялась от -0,1 °С до +1,1 °С, а в подстилающих ММП на глубине 4,0 м от дневной поверхности температура tp изменялась от -2,2 °С до -2,52 °С; засолённость грунтов на глубине 4,0 м от дневной поверхности Zp от 0 до 1,05 %, а в деятельном слое Z0 от 0 до 0,18 %; расчётные значения нагрузки от трубопровода Рртр включали массу фундамента, трубы и прокачиваемого флюида. В расчётах учитывалось, что трещины в плане имели четырёхугольную форму со сторонами блока 3,6 м, ширину раскрытия трещин 0,25 м, глубину 1,8 м. Результаты расчётов, приведенные в таблице, подтверждают существенную зависимость напряжений, деформаций оттаивающих грунтовых оснований с морозобойными трещинами от нагрузок, температур и засоленности грунтов. Следует отметить, что максимальные значения деформаций в оттаивающих грунтовых основаниях (151...112 мм) практически не отражались на устойчивости и надёжности эксплуатации пространственно гибких надземных трубопроводов. Аналогичные расчёты НДС грунтовых оснований под плитными фундаментами были проведены и на участках со слитыми ММП. Расчётные деформации в оттаивающих основаниях слитых ММП были на 5.1 % меньше значений, приведенных в таблице.
Предлагаемый усовершенствованный метод расчёта НДС грунтовых оснований под мелко заглублёнными и плитными фундаментами, учитывающий наличие в деятельном слое морозобойных трещин, сезонные изменения водно-теплового режима, физико-механические свойства грунтов, позволяет получать более достоверные показатели несущей способности оснований, устойчивости сооружений на деформируемых, засоленных ММП.
Таблица. Значения напряжений асж, кПа (числитель), деформаций s, мм (знаменатель) оттаивающих, неслитых ММП под фундаментом
№ п/п Характеристика разрезов, сложенных пылеватыми суглинками с показателем текучести Ji , равным 0,15 Расчётные значения напряжений в грунте*
Деформации грунта под фундаментом**
Расчётные нагрузки на фундамент, Рртр, (тс)
1,25 2,42 5,92
1 Zо = 0,0 %, zro = 0,0 %, U = -0,1 °С, tro = -2,3 °С 62,5* / 13,5** 115,7 / 22,9 272,3 / 46,6
2 Zо = 0,1 %, zro = 0,5 %, U = -0,3 °С, tro = -2,4 °С 60,2 / 14,3 122,6 / 26,3 288,1 / 49,8
3 Zо = 0,2 %, zro = 1,05 %, U = -0,5 °С, tro = -2,7 °С 58,3 / 15,9 138,3 / 29,8 296,2 / 52,2
4 Zо = 0,0 %, zro = 0,0 %, U = +0,5 °С, tro = -2,0 °С 64,6 / 29,6 148,3 / 61,2 303,6 / 151,3
5 Zо = 0,1 %, zro = 0,5 %, U = +0,1 °С, tro = -2,2 °С 61,8 / 37,1 156,7 / 65,7 312,1 / 160,2
6 Zо = 0,2 %, zro = 1,05 %, U = 0,0 °С, tro = -2,4 °С 59,1 / 52,6 169,3 / 68,4 312,1 / 160,2
7 Zо = 0,0 %, ZгD = 0,0 %, tпф = +1,1 °С, tro = -2,4 °С 66,8 / 59,8 156,7 / 102,3 312,3 / 162,8
8 Zо = 0,1 %, zro = 0,5 %, tпф = -0,1 °С, tгD = -2,5 °С 63,5 / 62,5 159,5 / 108,2 315,1 / 168,9
9 Zо = 0,2 %, zro = 1,0 %, tпф = +0,9 °С, tro = -2,5 °С 61,2 / 69,5 162,7 / 113,6 319,8 / 172,4
Ниже приведен расчёт устойчивости надземного нефтепровода 120*12 мм, полностью заполненного сырой нефтью, оснащённого теплогидроизолированным покрытием толщиной 140 мм (рис. 3). Нефтепровод возведён на свайных (в местах расположения неподвижных опор) и плитных, имеющих возможность пространственных перемещений, фундаментах (1). Фундамент (1) выполнен из искусственно-упрочнённого грунта и оснащён опорной стойкой (3) и армирующими диафрагмами (4). Для проведения расчёта были использованы следующие сведения. Трубопровод расположен в зоне распространения слитых ММП. В составе расчётного сценария не учитывались силы смерзания его подошвы с грунтовым основанием в осенне-зимние периоды, удерживающие силы от свайных фундаментов и углов поворота по трассе трубопровода. Сценарий рассматривает возможность опрокидывания штормовыми порывами ветра надземного нефтепровода на тундровую поверхность.
813
Худякова А.А. и др. Модель деятельного слоя многолетнемерзлых пород...
Удельные нормативные ветровые нагрузки g,,e в примере приняты для шестого климатического района и соответствуют: ge=73,0 кгс/м2 [84]. Расчётная площадь трубопровода Прн=\,2 м2, воспринимающая ветровую нагрузку Feem, определялась произведением диаметра надземного нефтепровода (5), имеющего теплогидроизоляционное (6) и гололедное в осенне-зимние периоды покрытие, на шаг расстановки опор фундаментов L по длине осваиваемой трассы. Шаг расстановки опор фундаментов по длине осваиваемой трассы L принят равным 12 м. Тогда .Рет=1,2 м х 12 м = 14,4 м2. Следует отметить, что в расчёте на устойчивость трубопровода учитывалась и вертикальная площадь каждого из плитных фундаментов. Поверхностные фундаменты выполнены в виде трапеций, имеющих в плане размеры 2,2х1,1 м, 0,8х0,4 м, высоту 1,1 м. Площадь плитных фундаментов для расчета устойчивости трубопровода равна 0,825 м2 (F% = [(0,4+1,1)/2] х 1,1 м2).
Тогда суммарная площадь контакта надземного нефтепровода и плитного фундамента составит Х^рвт+-^реф=14,4+0,825=15,225 м2. Проверка на устойчивость нефтепровода под ветровым напором включала определение опрокидывающего момента Мнопр = 73,0 • 15,225 • 1,8 = 2000,565 кгсм.
Можно рассчитать значение опрокидывающего момента
Мроп = М опр тст тдин = 2000,565 • 1,1 • 1,4 = 3080,87 кгс-м.
Расчётное значение момента (Мруд), удерживающего надземный нефтепровод от опрокидывания (рис. 3), можно рассчитать по формуле:
Муд=Мрпф+ Мтр+ Мн+Мис, (4)
где Мрпф - удерживающий момент поверхностного фундамента; М^тр - удерживающий момент стальной трубы; Мн - удерживающий момент прокачиваемой нефти; Мр1с - удерживающий момент изоляционного слоя.
Расчётное значение удерживающих моментов зависит от веса соответствующих элементов
тнэл = Гэл • унэл, (5)
где Унэл - объем элемента; унэл - объемный вес элемента.
Подставляя в выражение (5) соответствующие значения получим:
тнпф = [(2,2+0,8)/2 • (1,1+0,4)/2] • 1,1 • 1,6 = 1,98 тс; тнтр = 0,2095 -12 = 2,514 тс; тнн = 0,348 12 = 4,176 тс;
тнис = (3,14 • 0,62 - 3,14 •0,362) • 0,05 • 12 = 0,4341 тс.
Расчетное значение удерживающего момента от опрокидывания надземного трубопровода на тундровый покров равно:
Мруд = Мнуд • кз, (6)
где к . - коэффициент запаса.
Рис. 3. Вертикальный разрез поверхностного фундамента для надземного нефтепровода
1 - поверхностного типа фундамент;
2 - песчано-гравийная подготовка;
3 - опорная стойка;
4 - армирующие диафрагмы;
5 - нефтепровод;
6 - теплогидроизоляционный слой;
7 - опорная плита стойки;
8 - стальной лист;
9 - опорная площадка трубы.
814
Вестник МГТУ, том 13, №4/1, 2010 г.
стр.810-815
Суммарное максимальное значение удерживающего момента от опрокидывания надземного трубопровода на тундровый покров равно:
Мрудт“ = [1,98+2,514+4,176+0,434] • 1,1 • 0,9 = 9,013 тсм.
Суммарное минимальное значение удерживающего момента от опрокидывания надземного трубопровода на тундровый покров равно:
Мрудт,т = [1,98+2,514+0,434] • 1,1 • 0,9 = 4,879 тсм.
Следовательно, Мроп = 3,081 тс-м < Мрудтт = 4,879 тс-м < Мрудтас = 9,013 тс-м.
Выполненный расчет подтверждает устойчивость поверхностных фундаментов в случаях воздействий на надземные нефтепроводы опрокидывающих усилий без использования на их трассах промежуточных свайных опор при сокращении числа компенсаторов. На деформируемых грунтах в качестве промежуточных опор для надземных трубопроводов следует применять поверхностные фундаменты с учётом возможностей их пространственных перемещений при сезонных вертикальных перемещениях поверхностных отложений, а также при температурных горизонтальных деформациях непосредственно надземных нефтепроводов. Целесообразность применения поверхностных фундаментов подтверждена и результатами экономического сравнения их стоимости с затратами на возведение свайных промежуточных подвижных опор. Расчёты показывают, что частичное (до 80 %) применение плитных фундаментов снижает стоимость строительства нефтепроводов более чем на 32 %.
4. Выводы
1. Наличие в покровных отложениях блоков, разобщённых между собой трещинами, существенно отражается на вертикальных и горизонтальных составляющих тепломассопереноса. Предложена модель деятельного слоя ММП, в котором процессы миграции влаги одновременно происходят внутри блоков не только в вертикальном, но и в горизонтальном направлении путем ее разгрузки в систему трещин.
2. Усовершенствована методика расчёта НДС промерзающих-оттаивающих оснований. Она учитывает влияние морозобойных трещин, вертикальных и горизонтальных составляющих тепломассообмена, что позволяет повысить достоверность результатов расчетов, надёжность, продолжительность эксплуатации трубопроводов в северных регионах и позволяет более широко применять поверхностные фундаменты.
3. В качестве подвижных опор надземных трубопроводов на деформируемых засоленных грунтах предлагается использовать плитные фундаменты, которые снижают стоимость строительства до 32 %, повышают надежность эксплуатации нефтепроводов за счет возможности обеспечения горизонтальных перемещений и компенсации вертикальных подвижек. Выполненный расчет на устойчивость подтверждает возможность и целесообразность их применения в НАО.
Литература
Винокуров Р.С., Конюхов А.В., Садриев Р.А. Анализ результатов геокриологического мониторинга по трассе нефтепровода Южно-Шапкино-Харьяга. Вестник АГТУ, серия Прикладная геоэкология, Архангельск, вып. 75, с.117-124, 2008.
Худякова А.А. Совершенствование строительства нефтепроводов в северных регионах. Автореферат дисс. ... канд. техн. наук, Уфа, ИПТЭР, 24 с., 2010.
815
