Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

THE IMPACT OF ENVIRONMENT ON THERMAL STATE OF OFFSHORE PLATFORMS

Starokon I.V. 1
1 Russian State oil and gas university named after I.M. Gubkin
Offshore Platforms are exposed to solar radiation, which causes them to heat. The change in the thermal condition (TC) structural elements (FE) results in a different direction and magnitude of thermal stresses, whose analysis by the author in [4]. Previously, the author has developed a method [3], which allows to calculate the temperature of the elements without the influence of the environment. This paper presents the augmented technique that allows us to estimate the effect of this in terms of free and forced convection. In addition to the formula to estimate the surface temperature of the elements introduced a correction factor which takes into account their degree of blackness. On the basis of the developed technique is an example of the calculation of the thermal state of offshore platforms in the Black Sea. Revealed a difference between the TC columns, diagonals and horizontal elements.
marine
offshore platforms
solar radiation
temperature
exposure
thermal
flow
heat transfer

Морские стационарные платформы подвергаются различным воздействиям, влияющим на их напряженное состояние. Одним из таких воздействий является температурное, вызывающее переменные температурные напряжения в сечениях их конструктивных элементов. Ранее, в работе [3] автором было доказано это влияние и обоснована методика его численной оценки. Однако эта методика не учитывала теплообмена с окружающей средой. Целью настоящей статьи является исследование влияния окружающей среды на формирование температуры внешней стенки конструктивных элементов МСП.В качестве примера будут рассмотрены тепловые процессы опорных блоков морских стационарных платформ, расположенных на месторождениях нефти и газа Черного моря. Анализ проектов МСП показал, что конструктивно опорные блоки МСП выполнены из стальных труб различного диаметра от 325мм до 1020мм. Как это было показано в работах [1,3-5], опорный блок МСП подвергается воздействию солнечной радиации, что вызывает изменение его температурного состояния. Автором в работе [3] была предложена методика, позволяющая численно оценить величину этого воздействия, которая заключается в вычислении величины плотности теплового потока, воздействующего на горизонтальные, вертикальные и наклонные конструктивные элементы, и предложена новая формула расчета температуры внешней стенки конструктивных элементов(КЭ) МСП в зависимости от динамики изменения плотности теплового потока. Однако за рамками статьи [3] осталось оценка влияния окружающей среды на формирование температуры внешней стенки исследуемого объекта. В настоящей статье проводится изучение этого влияния. Из классической теории термодинамики нам известно, что при воздействии теплового потока на твердое тело, часть радиации поглощается и нагревает тело, а часть отражается. Соотношение между поглощенной и отраженной солнечной радиацией определяется так называемой «степенью черноты». Что же представляют собой КЭ МСП? Это стальные трубы, с нанесенным на него лакокрасочным покрытием, предохраняющим его от коррозионного воздействия и имеющим свою, в некоторых случаях весьма значительную, отражающую способность. Однако, как показала практика эксплуатации опорных блоков, фактически разрушение лакокрасочного покрытия (ЛКП) происходит уже через один, два года после его нанесения в зависимости от марки и завода изготовителя. В связи с тем, что процесс восстановления ЛКП является чрезвычайно дорогостоящим, его проводят один раз в 10 лет. Кроме того, на практике не в достаточной мере соблюдаются требования к предварительным операциям по подготовке поверхности КЭ МСП к нанесению ЛКП (например, зачистка поверхности), что также снижает стойкость лакокрасочного покрытия к разрушению. В свою очередь, некачественное ЛКП приводит к усилению коррозионных процессов. Фактическое состояние элементов опорного блока можно описать как «стальные трубы с сильно окисленной поверхностью», для которых приведенная степень черноты согласно справочным данным составляет 0,74, т.е. другими словами, при воздействии солнечной радиации на элементы опорного блока лишь 26% этой энергии отражаются, остальные же 74% – поглощаются и нагревают элементы опорного блока. В свою очередь нагретый конструктивный элемент излучает энергию в окружающую среду. Процесс теплоотдачи происходит за счет излучения и конвекции, которая подразделяется на –свободную и вынужденную. Все эти процессы характеризуются коэффициентами теплоотдачи. При теплообмене излучением, коэффициент теплоотдачи αл вычисляется по формуле:

(1)

где: ε - степень черноты КЭ МСП, принимаемая равной0,74; С - коэффициент излучения абсолютно черного тела, равный 5,67 Вт/(м2К4); Тст и Твозд температура стенки КЭ МСП и температура воздуха в градусах Цельсия. При свободной или вынужденной конвекции коэффициент теплоотдачи α можно вычислить по формуле:

,(2)

Где: α - коэффициент теплоотдачи, Вт/м2К; Nu-число Нуссельдта; D - диаметр конструктивного элемента МСП, м; λ - коэфициент теплопроводности воздуха, 10-2 (Вт/мК). В свою очередь, число Нуссельдта определяется по формуле:

(3),

где: С и n значения, выбираемые из таблицы 1; Gr-число Грасгофа; Pr-число Прандтля при определенной температуре воздуха; Prст-число Прандтля воздуха, при температуре стенки КЭ МСП.

Таблица 1

Значение коэффициента С и показателя степени

Режим

Gr·Pr

C

n

Ламинарный

1·103 ÷ 5·102

1,18

0,125

Переходный

5·102 ÷ 2·107

0,54

0,25

Турбулентный

2·107 ÷ 1·1012

0,185

0,33

В свою очередь число Грасгофа определяется по формуле:

(4),

где g – ускорение свободного падения, g=9,81 м/с2; D – диаметр КЭ МСП в метрах; β – температурный коэффициент объемного расширения воздуха, равный ; tc и tv – температуры стенки КЭ МСП и воздуха. Число Прандтля является справочной величиной. Процесс вынужденного конвективного теплообмена имеет свои особенности. Анализ данных из справочников по климату Черного моря позволят утверждать, что наиболее сильные штормовые ветра, способные серьезно повлиять на формирование температуры КЭ МСП, приходятся с ноября по март в период минимальной солнечной активности. Так, например, в прибрежной зоне СВЧЧМ диапазон значений средней скорости ветра изменяется от 2,7–2,8 м/с на ЮБК до 6–7 м/с, в районах Керченского пролива и Кавказа. С другой стороны, с апреля по ноябрь происходит ослабление ветровой активности, скорость ветра уменьшается до 1,9–2,4 м/с и 3,3–5,3 м/с соответственно, при этом значительно усиливается интенсивность солнечного излучения. Исходя из этого, автором для рассмотрения был выбран период с апреля по ноябрь, т.к. именно в этот период происходят наиболее значительные изменения температурного состояния КЭ МСП. Совершенно очевидно, что при обтекании элементов МСП ветровым потоком происходит сложный теплообмен, складывающийся из двух составляющих: вынужденный конвективный теплообмен и теплообмен излучением. Большое значение на формирование температуры КЭ МСП имеет ориентация ветрового потока по отношению к ним. Рассмотрим случай поперечного обтекания ветровым потоком элементов КЭ МСП. Основным критерием, определяющим интенсивность теплообмена, является коэффициент теплоотдачи α. Для вычисления этого коэффициента необходимо знать следующие исходные данные: диаметр конструктивного элемента МСП, скорость ветрового потока, его направление по отношению к рассматриваемому элементу. При расчете числа Нусельдта большое значение имеет ориентация потока, определяемая углом µ. Среднее по периметру число Нусельдта определяется с учетом угла атаки Φ при помощи умножения на поправочный коэффициент . Для рассматриваемых нами условий обтекания ветром элементов МСП число Нусельдта определяется из условий

1) При Re=10…103 2) При Re=103…2*105 3) При Re более 2*105 . В этих условиях приняты следующие обозначения: Re-число Рейнольдса, Nu-число Нуссельдта, -коэффициент, зависящий от направления ветра и выбираемый по таблице 1, С=0,023 и m=0,8 коэффициенты, принимаемые в соответствии с работой [2], Pr-число Прандтля, принимаемое в соответствии с таблицей 3.

Таблица 2

Значение коэффициента в зависимости от направления ветра

Угол ?

10

20

30

40

50

60

70

80

90

0,42

0,52

0,67

0,78

0,88

0,94

0,98

1

1

В свою очередь, число Рейнольдса определяется по формуле:

, (5)

Где: w– скорость ветра, м/с;D – диаметр конструктивного элемента МСП, м; v – кинематическая вязкость воздуха м2/c, которая определяется из таблицы3.

Таблица 3

Некоторые параметры свойств воздуха при различных температурах

Температура воздуха

0

10

20

30

40

50

60

Кинематическая вязкость, 10-6м2/с

13,28

14,16

15,06

16

16,96

17,95

18,97

Число Прандтля

0,715

0,714

0,713

0,712

0,711

0,710

0,709

Коэфициент теплопроводности воздуха λ, 10-2 (Вт/мК)

2,44

2,51

2,59

2,67

2,76

2,83

2,9

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

,(6)

Затем следует вычислить плотность теплового потока от стенки трубы к более холодной среде набегающего ветрового потока по формуле:

(7)

q – плотность теплового потока от стенки трубы к набегающего ветровому потоку, Вт/м2; α – коэффициент теплоотдачи, Вт/м2К; температура внешней стенки КЭ МСП; – температура набегающего ветрового потока; коэффициент, учитывающий что в рассматриваемом случае воздействию прямого солнечного излучения подвергается только половина трубы, коэффициент принят равным 0,5.В статье [3] автором предлагалась методика, позволяющая рассчитать температуру поверхности КЭ МСП в зависимости от плотности теплового потока, возникающего от солнечного излучения:

(8)

где: Qкэ – суммарный тепловой поток солнечного излучения, R и L – внешний радиус и длина конструктивного элемента МСП в метрах, Sнакл – плотность теплового потока прямого солнечного излучения, па­дающего на наклонную поверхность конструктивного элемента МСП, Вт/м2; Dнакл – плотность теплового потока рассеянной сол­нечной энергии, падающей на наклонную поверхность конструктивного элемента МСП, Вт/м2; К – поправочный коэффициент, учитывающий отраженный от поверхности морской воды тепловой поток. Методика расчета этих значений приведена в статье [3].Помимо этого, автором разработана формула, позволяющая определить температуру внешней стенки КЭ МСП, подвергающейся воздействию солнечного излучения (9):

, (9)

где: Тi и Тi+1 – начальная и рассматриваемая в некоторый момент времени температура поверхности конструктивного элемента МСП, находящегося под действием прямого солнечного излучения; R – внешний радиус конструктивного элемента МСП, м; l – его длина, м; m-масса рассматриваемого участка; 0,74 – поправочный коэффициент, учитывающий степень черноты трубы, а остальные значения приведены в статье [3].В связи с тем, что начальная и рассматриваемая в некоторый момент времени температура поверхности конструктивного элемента МСП являются неизвестными величинами, автором предлагается следующая методика: 1) рассчитать по формуле (8) суммарные величины плотностей тепловых потоков от солнечного излучения, и определить по формуле (9) температуру внешней стенки КЭ МСП; 2) зная данные о параметрах ветрового потока (скорость, температура и др.), рассчитать плотность теплового потока q, создаваемого ветровым воздействием по формуле (7); 3) скорректировать величину плотностей тепловых потоков по формуле (9); 4) используя полученную величину , определить температуру поверхности конструктивного элемента МСП. Рассмотрим пример. Исследуем процесс нагрева и охлаждения МСП на примере 14 августа 2010г. Данные о климатических условиях получены из архива погоды на сайте meteo.ua. Рассмотрим горизонтальные конструктивные элементы МСП, выполненные согласно проекту из стальных труб диаметром 325х8мм (масса исследуемого элемента 312,5 кг) и 530х12мм (масса исследуемого элемента 767,7 кг). Также проанализируем вертикальные элементы МСП (колонны), выполненные согласно проекту из стальных труб диаметром 720х20мм (масса исследуемого элемента 1739 кг)и 1020х20мм (масса исследуемого элемента 2474 кг). Для сравнительного анализа примем длину всех элементов одинаковой и равной 10м.Процесс нагрева начинается после восхода солнца в 5:00. Примем коэффициент прозрачности атмосферы с=0,38, при этом облачность отсутствует. Создаваемая плотность теплового потока на горизонтальные элементы (ГЭ) к 5:00 имеет значения 13Вт\м2, что практически не влияет на формирование температуры поверхности КЭ МСП. Уже к 7:00 плотность теплового потока на ГЭ достигает значений 384 Вт\м2, и уже к 9:00 плотность теплового потока составляет 845 Вт\м2. Затем, постепенно возрастая в логарифмической зависимости плотность теплового потока достигает своего максимального значения 1137 Вт\м2 к 12:00 часам. После этого начинается уменьшение плотности теплового потока, которое к 15:00 снижается до 845 Вт\м2, к 17:00 составляет 384 Вт\м2,а к 19:00 составляет всего лишь 19 Вт\м2.Несколько иначе формируются плотности тепловых потоков от солнечного воздействия на вертикальные элементы (колонны). Следует отметить, что величины плотностей тепловых потоков при схожих условиях для колон практически в три раза меньше, чем для горизонтальных элементов. Так в момент максимальных значений тепловых потоков 12:00 в рассматриваемом примере для горизонтальных элементов это значение составляет 1137 Вт\м2, а для колон всего лишь 348 Вт\м2. Подставляя полученные значения в формулу (9) получим следующие значения (таблица 4).

Таблица 4

Расчетные значения температуры горизонтальных и вертикальных элементов опорного блока МСП при различных условиях окружающей среды

Часы

7:00

8:00

9:00

10:00

11:00

12:00

7:00

8:00

9:00

10:00

11:00

12:00

Плотность теплового потока солнечного излучения , Вт/м2

384

626

845

1014

1106

1137

190

307

367

389

360

348

Плотность теплового потока поглощенного конструктивными элементами МСП, Вт/м2

284

463

625

750

818

841

141

227

272

288

266

258

Температура воздуха, С0

28

29

30

30

31

32

28

29

30

30

31

32

Скорость ветра, м\с

1

2

1

2

1

1

1

2

1

2

1

1

Наименование параметра

Горизонтальный элемент МСП 325х8мм

Вертикальный элемент (колонна) МСП 720х20мм

Температура внешней поверхности КЭ МСП без учета влияния окружающей среды, С0

56

79

100

115

124

127

29

32

34

35

33

33

Коэффициент теплоотдачи излучением, αл Вт/(м2С0)

5,3

7,1

6,8

7,3

7,6

7,8

4,6

5,7

4,8

4,8

4,8

4,8

Коэффициент теплоотдачи при свободной конвекции, α Вт/(м2С0)

5,9

7,1

7,9

7,9

8,1

8,1

2,6

3,8

4,1

4,5

3,3

2,6

Коэффициент теплоотдачи при вынужденной конвекции, α Вт/(м2С0)

7

10

7

10

7

7

5

7,5

5

7,5

5

5

Плотность отводящего теплового потока при теплообмене излучением и свободной конвекцией

157

355

566

707

793

811

4

14

18

23

8

4

Температура внешней поверхности КЭ МСП с учетом теплоотдачи излучением и свободной конвекцией, С0

49

55

83

86

87

89

27

31

33

33

33

33

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

172

428

531

804

737

755

5

20

20

31

10

5

Температура внешней поверхности КЭ МСП с учетом теплоотдачи излучением и вынужденной конвекцией, С0

47

45

73

60

80

82

27

30

33

33

33

33

Наименование параметра

Горизонтальный элемент МСП 530х12мм

Вертикальный элемент (колонна) МСП 1020х20мм

Температура внешней поверхности КЭ МСП без учета влияния окружающей среды, С0

44

59

73

83

89

91

29

32

34

35

33

33

Коэффициент теплоотдачи излучением, αл Вт/(м2С0)

5,0

6,5

5,8

6,0

6,2

6,4

4,6

5,7

4,8

4,8

4,8

4,8

Коэффициент теплоотдачи при свободной конвекции, α Вт/(м2С0)

6,3

7,6

8,4

8,8

8,9

9,0

2,6

3,8

4,1

4,5

3,3

2,6

Коэффициент теплоотдачи при вынужденной конвекции, α Вт/(м2С0)

5,6

8,5

5,6

8,5

5,6

5,6

5

7,5

5

7,5

5

5

Плотность отводящего теплового потока при теплообмене излучением и свободной конвекцией

90

212

305

392

438

454

4

14

18

23

8

4

Температура внешней поверхности КЭ МСП с учетом теплоотдачи излучением и свободной конвекцией, С0

36

41

47

50

52

53

27

31

33

33

33

33

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

85

225

245

384

342

354

5

20

20

31

10

5

Температура внешней поверхности КЭ МСП с учетом теплоотдачи излучением и вынужденной конвекцией, С0

37

40

52

51

60

61

27

30

33

33

33

33

Как видно из таблицы, формирование температуры поверхности КЭ МСП зависит от многих факторов. На первый взгляд, казалось бы, что чем больше диаметр, а, следовательно, и площадь облучаемой поверхности, тем более высокие температуры должны формироваться на поверхности КЭ МСП. Однако с увеличением диаметра, как правило, увеличивается и толщина стенки КЭ МСП, а это приводит к увеличению массы элемента. А чем больше масса элемента, тем менее высоких температур он будет достигать при нагреве. Следует также отметить одну немаловажную особенность наклонных элементов, а именно – что при их отклонении по отношению к горизонту до угла 250происходит увеличение плотности теплового потока солнечного излучения. Так при одинаковых условиях плотность теплового потока горизонтально расположенного КЭ составляла бы 1137 Вт\м2, в то же время для элемента наклоненного под углом 250 эта плотность составляет уже 1254 Вт\м2. В проанализированных автором проектах наклонные элементы типа раскосов расположены под углом 570 по отношению к горизонту, в свою очередь это вызывает снижение максимальной плотности теплового потока до 1072 Вт\м2. Это явление влияния угла наклона плоскости на плотность теплового потока от солнечного излучения хорошо изучено в области проектирования солнечных энергоустановок и разработаны методы выбора оптимальных углов наклона плоскости, позволяющие аккумулировать максимально возможную плотность теплового потока от солнечного излучения. Изменяются и коэффициенты теплоотдачи (КТ). При скоростях ветрового потока приблизительно 1 м/с КТ сопоставимы с КТ при свободной конвекции, однако при увеличении скорости ветра до 2 м/с коэффициенты теплоотдачи существенно возрастают. В результате проведенных расчетов установлено, что в отличие от ГЭ для колонн температурное воздействие незначительное.

Рецензенты:

Бородавкин П.П., д.т.н., профессор, генеральный директор АО «Интергаз», г. Москва.

Литвин И.Е., д.т.н., генеральный директор ООО «СТД», г. Москва.