Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Сооружение подводных трубопроводов

..pdf
Скачиваний:
30
Добавлен:
15.11.2022
Размер:
15.38 Mб
Скачать

подходах к реке можно выявить точки (узлы), через которые может быть проложен трубопровод (рис. 1.6). Эти точки могут быть соединены дугами, каждой из которых может быть постав­ лена в соответствие, после анализа условий строительства, стоимость сооружения соответствующего участка трубопровода Су, ще ij - номера точек (узлов) начала и конца дуги. За­

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

(//) lJ

по этому пути.

Задача является многовариантной и для ее решения можно использовать динамическое программирование, рекуррентные со­ отношения которого записываются в виде:

/„</> “ min {fnA(i) + С..},

(n - 2, 3, ...);

u ©

(1.15)

(1.16)

ще А - номер шага вычислений, соответствующий номеру сече­ ния графа.

Численные результаты могут быть получены с помощью на­ писанной на БЕЙСИКе программы для персональной ЭВМ. Пример расчета приводится в табл. 1.1.

Исходные данные:

Граф содержит 6 сечений и 20 узлов, число узлов, находя­ щихся:

в 1-м сечении равно 1; во 2-м 3; в 3-м 6; в 4-м - 6; в 5-м - 3; в 6-м - 1. В 1-м сечении находится узел графа с

номером 1, во 2-м

- с

номерами 5;

3; 2, в 3-м

- 4;

6;

/,

8;

9; 10, в 4-м - 11;

12;

13; 14; 15;

16, в 5-м -

П ;

18;

19;

в

6-м - с номером 20. Результаты расчета показывают, что уча­ сток и створ подводного перехода должны располагаться на

трассе,

проходящей через

узлы графа 1 - 5 - 6 - 1 2 - 1 7

-

20;

при

этом

значение

критерия

оптимальности равно 2,6

+

+ 2,8 +1,7 +

2,3 +

3,4

- 12,8 млн руб. Отметим,

что ви­

зуально

краткий путь,

проходящий

через узлы графа

1 - 3

-

7 -

13

- 18 -

20, характеризуется

большим значением крите­

рия:

3,9

+ 2,9

+ 1,6 + 2,4+ 3,5 - 14,3 млн руб.

 

 

При выполнении расчета персональная ЭВМ работает в инте­

рактивном режиме (режиме диалога

между человеком

и ЭВМ).

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

Стоимостные харакгеристакн строительства

 

 

Таблица 1.1

 

 

 

Но­

Номер

Номер

Значение

Но­

Номер

Номер

Значение

мер

началь­

конеч­

критерия

мер

началь­

конеч­

критерия

дуги

ного

ного

для дуги,

дуги

ного

ного

для дуги,

 

узла

узла

млн руб.

 

узла

узла

млн руб.

 

дуги

дуги

 

 

дуги

дуги

 

1

1

2

3,3

18

7

13

1,6

2

1

3

3,9

19

8

14

1,5

3

3

4

3,3

20

9

15

1.8

4

5

4

2,4

21

10

16

1,2

5

1

5

2,6

22

11

17

2,5

6

3

6

2,7

23

12

17

2,3

7

5

6

2,8

24

13

17

2,7

8

3

7

2,9

25

12

18

2,9

9

5

7

3

26

13

18

2,4

10

2

8

3,1

27

14

18

3

11

3

8

3,2

28

14

19

3,1

12

2

9

3,3

29

15

19

3,7

13

3

9

3,4

30

16

19

3,3

14

2

10

3,5

31

17

20

3,4

15

3

10

3,6

32

18

20

3,5

16

4

11

1,9

33

19

20

3,9

17

6

12

1,7

 

 

 

 

Отладка программы выполнена на мини-ЭВМ типа ’’Искра 226”. При необходимости расчеты могут быть выполнены на персональ­ ных и мини-ЭВМ других типов, операционные системы которых содержат БЕЙСИК.

1.4. ОПТИМИЗАЦИЯ ПРОФИЛЯ

Себестоимость строительно-монтажных работ существенно за­ висит не только от выбора участка и створа перехода, но и от заглубления трубопровода, определяющего объем выполняемых подводных земляных работ. Стоимость выполнения подводных земляных работ достигает 50% от общей стоимости строитель­ ства перехода. Причины этого кроются не только в значитель­ ной стоимости выполнения единицы объема подводных земляных работ, обусловленной применением дорогостоящей техники, но также и в необходимости выполнения этих работ в больших объ­ емах, что связано со значительной шириной траншеи по дну, достигающей 6-10 м, и незначительным заложением откосов, равным иноща 1/3. Объем земляных работ зависит не только от указанных параметров, но и от глубины подводной траншеи, ко­

торая может достигать на реках со

значительными сезонными

и многолетними переформированиями

русла 6 м и более.

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

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

3-271

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

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

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

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

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

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

Z(x)

-

M[Z(x)]

(1.17)

и корреляционная функция

 

К(х,

хд

- М {[Z(x)-Z(x)] [Z(Xl)-Z(Xl)]},

(1.18)

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

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

Практическое нахождение оценок математического ожидания и корреляционной функции может быть выполнено с использованием методов математической статистики. Если при русловых съемках получено п реализаций случайной функции, то ординаты Zi(x), Z2(х), ..., Zn(x) реализаций профиля дна реки в створе пере­

хода для значения абсциссы х можно рассматривать как найден­ ные из опыта значения случайной величины Z(x). Тоща оценки математического ожидания и корреляционной функции могут быть найдены по формулам

Z (х)

 

Z

Z.(x),

 

 

(1.19)

 

 

П 1-1

1

 

 

 

 

К(х,

х\)

 

1

Е

ZXx)Z.(xi) - Z(x)Z(xi)

( 1. 20)

л - 1

 

 

j - 1

1

1

 

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

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

3*

Рис. 1.7. План (а) и профили (б) створов подводного перехода через р.Кума нефтепровода Ходмогоры-Европейская часть России:

1-9

- створы; 10 - верхнее положение трубопровода по створу № 5; 11

бе­

рег

реки

 

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

У * >

ехр

f v n ( x ) - Z ( x ) 1 2

( 1.21)

2n<rz (x)

2crz (х)

 

 

 

где °'z (x)

среднеквадратичное отклонение

случайного

значения высотной отметки дна реки в точке с абсциссой х от математического ожидания

<rz (je)

-

V К(х, Xj);

 

 

(1.22)

сг^(дс)

 

 

среднеквадратичное

отклонение

угла

наклона

случайного

профиля дна реки

по створу перехода в

сечении

с абсциссой х

 

 

 

Ос)

-

/ ^

— 1d— К(х, Xj).

 

 

(1.23)

^d х 2

Расчет числа потенциальных размывов может быть выполнен на миниили персональной ЭВМ с помощью программы, разрабо­ танной на языке БЕЙСИК. На рис. 1.7 представлены исход­ ные данные для примера расчета. Ввиду небольшой длины про­ филь подводного трубопровода принят прямолинейным. Расчеты на ЭВМ показали, что при заглублении трубопровода на 0,5 м ниже дна реки следует ожидать 0,019 размывов, а при заглуб­

лении на

2 м

число размывов будет

равным 0,0011, т.ё.

в 17 раз меньше.

трубопровода вместе с

Однако

при

увеличении заглубления

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

в

некоторых

случаях

расходы не только на укрепление рус­

ла и

берегов, но и ремонт собственно трубопровода.

 

В

данной связи правомерной является постановка задачи

определения

оптимального профиля

подводного трубопровода.

В

качестве

критерия

оптимальности

рассматриваемой задачи

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

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

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

работкой статистических данных по эксплуатации подводных пе­ реходов.

Затраты на выполнение земляных работ, включающие рытье траншеи и ее обратную засыпку (замыв), могут быть определены с учетом глубины разрабатываемой траншеи, имеющей обычно в поперечном сечении трапецеидальную форму (см. рис. Г.7, 7). Обозначив через Z(x), уОс) высотные отметки соответ­ ственно дна реки и дна подводной траншеи, затраты на выпол­ нение земляных работ при сооружении участка подводного тру­ бопровода единичной длины можно определить по формуле:

Сз(у, х) - [(Z(x) - у(х))(А + m(Z(x) - у(*))]х

x(Ct + Сг) - Ci

,

(1.24)

где Ci, С2 стоимости соответственно разработки и обратного замыва единицы объема траншеи; Ь> т ~ ширина траншеи по дну,

коэффициент заложения откоса; De диаметр забалластирован­ ного трубопровода.

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

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

 

L

Сз(у,

x)dx

L

_

+ D6, x)dx ►,

(1.25)

min

« S

+ Ср S

v(y

 

о

 

 

О

 

 

 

ще

v(y

+ De,

х)

отнесенное

к единице длины

трубопровода

число размывов.

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

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

кривизна обратна радиусу изгиба,

это условие можно записать

в виде

 

 

«

— - ,

(1.26)

dx2I

Лм

 

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

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

EJ

+ Р(х) - д,

(1.27)

 

dx4

 

условие прилегания трубопровода к основанию и наличие отпора можно записать в виде

где

q

вес трубощювода под

водой;

Р(х)

отпор основания;

EJ -

жесткость трубопровода.

 

 

 

 

Производя вычитание неравенства (1.28) из (1.27), условие

прилегания представим окончательно в виде

 

^

-

< - 2 -

 

 

(1.29)

d x 4

EJ

 

 

 

Таким образом, постановка задачи оптимизации профиля под­

водного трубопровода состоит

из

критерия оптимальности

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

ний. Обозначив через

С* (у,

х) =

С3(у, х)

+ v(y + As, х),

введем функцию-свертку

 

ч

 

 

 

 

Л d l l < JZ_

X) ;

d h \

 

 

d x 2 '

* « J

d x 4

EJ

S(y, х) =• -

IA Л . d * 2 '

V

il±

^ a -

dx<

EJ)

 

(1.30)

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

min

L

S(y, x)dx.

(1.31)

J*

у (JC) 0

 

 

 

После введения функционала

 

W(y,

х)

- min

S S(y, x)dx

(1.32)

 

 

у (х)

0

 

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

W(y,

х) = min [W(y(x

- A), * - A)

+

x

y( JC-A)

 

 

S(y, x) , dx], x

(0, L ) ;

(1.33)

+ S

x- L