Содержание к диссертации
Введение
Глава 1. Постановка задачи 20
1.1. Основные закономерности установившегося течения жидкости 20
1.2. Математическая модель потокораспределения 25
1.3. Внешнее и внутреннее представление модели гидравлической сети...27
Глава 2. Расчет стационарного потокораспределения в многокольцевых гидравлических сетях 33
2.1. Метод последовательных приближений 33
2.2. Решение системы уравнений (2.1) на каждом шаге итерации 38
2.3. Алгоритм основного цикла решения задачи потокораспределения... 45
2.4. Процедуры для учета регуляторов расхода, давления и обратных клапанов ...49
2.4.1. Вычисление сопротивлений обратных клапанов 50
2.4.2. Вычисление сопротивлений регуляторов расходов 51
2.4.3. Вычисление сопротивлений регуляторов давления 53
2.5. Формирование расчетной схемы на основе принципиальной схемы... 55
2.6. Методика расчета гидравлических сопротивлений участков расчетной схемы 58
2.6. 1. Расчет сопротивлений труб и арматуры 58
2.6. 2. Расчет сопротивлений источников напора 62
2.6. 3. Расчет гидравлического сопротивлений потребителей 62
Глава 3. Автоматическое эквивалентирование схем 66
3.1. Эквивалентирование при отключении отдельных участков 66
3.2. Эквивалентирование последовательно соединенных участков 67
3.3. Эквивалентирование сетей, связанных с магистральным трубопроводом одной врезкой 69
3.4. Эквивалентирование при решении задачи наладки тепловых сетей 72
3.5. Эквивалентирование при совместной работе источников на общую теплосеть с несколькими точками подпитки 76
Глава 4. Расчет тепловых потерь и температуры теплоносителя 80
4.1. Расчет температуры теплоносителя 80
4.2. Расчет тепловых потерь 82
4.2.1. Условные обозначения 82
4.2.2. Расчет тепловых потерь для участков трубопровода надземной прокладки 83
4.2.3. Расчет для участков трубопровода с прокладкой в помещениях и тоннелях 84
4.2.4. Расчет для участков трубопровода подземной канальной прокладки 84
4.2.5. Расчет для участков трубопровода подземной бесканальной прокладки 85
4.2.6. Реализация в системе выбора поправочных коэффициентов и констант расчета 86
4.3. Расчет тепловыделения в системах отопления 87
4.4. Расчет тепловыделения в водоподогревателях ГВС 89
Глава 5. Внутренняя модель принципиальной гидравлической схемы теплосети 90
5.1. Структура объектов БД ПГСТ 91
5.1.1. Абстрактный базовый класс 91
5.1.2. Виртуальные методы абстрактного базового класса 92
5.1.3. Классы и структуры данных, используемые методами абстрактного базового класса 94
5.1.4. Основные геометрические классы 97
5.1.5. Классы гидравлических элементов 100
5.2. Структура базы данных потребителей 108
Глава 6. Моделирование и расчет гидравлических режимов нефтепроводов . 112
6.1. Математическая модель потокораспределения нефтепровода 112
6.2. Укрупненная расчетная схема нефтепровода 112
6.3. Принципиальная технологическая схема нефтепровода 115
6.4. Внутренняя модель гидравлической схемы нефтепровода 116
6.5. Автоматическое эквиваленч ирование принципиальной технологической схемы нефтепровода 120
6.6. Поиск оптимальных режимов нефтепроводов 124
Заключение 129
Литература 130
- Основные закономерности установившегося течения жидкости
- Решение системы уравнений (2.1) на каждом шаге итерации
- Эквивалентирование сетей, связанных с магистральным трубопроводом одной врезкой
- Расчет температуры теплоносителя
Введение к работе
Возросшие на данный момент времени требования к экономии энергоресурсов делают очевидным необходимость применения эффективных алгоритмов расчета и оптимизации на ЭВМ для обеспечения оптимального функционирования сложных энергетических систем, таких как теплосети и нефтепроводы. Для сложных энергетических систем характерны большие объемы информации и сложность взаимосвязей между элементами систем, поэтому для применения на практике инженерами эффективных алгоритмов расчета и оптимизации требуется программное обеспечение, позволяющее инженерам моделировать энергосистемы в привычном для них представлении, в терминах их предметной области.
В связи с этим для современных программных продуктов важны эффективные алгоритмы реализации расчетных задач потокораспределения, адекватные модели гидросистем и эффективные среды интерфейса с пользователем, позволяющие инженерам работать с моделью гидросистемы, не выходя за пределы предметной области.
Наиболее сложными для моделирования являются тепловые сети.
По мере укрупнения ТЭЦ, увеличения радиуса транспортировки энергоносителя повышаются требования к надежности теплоснабжения, т. е. к бесперебойной подаче теплоты потребителю. Поскольку при современных конструкциях тепловых сетей и качестве их выполнения не исключается аварийное или плановое отключение отдельных участков сети, основным методом повышения надежности в настоящее время является сооружение резервных связей между тепломагистралями одной или нескольких ТЭЦ. В результате прокладки перемычек между магистралями тепловые сети развиваются как сложные многокольцевые гидравлические системы.
По условиям надежности многокольцевыми являются не только тепловые сети, но и электрические, водопроводные, нефтепроводные. газовые и вентиляционные сети. Однако гидравлический режим тепловых сетей значительно сложнее, чем режим водопроводных или газовых сетей. Во-первых, потому что тепловые сети являются двухтрубными. Обе магистрали (подающая и обратная), гидравлически связаны между собой и сотнями перемычек - потребителями тепла. Во-вторых, возможность вскипания воды и непосредственное присоединение к тепловой сети тысяч местных систем отопления, обладающих невысокой механической прочностью, жестко ограничивают максимальные и минимальные допустимые давления и напоры в подающей и обратной магистралях. Основные ограничения состоят в следующем:
• Напор в местных системах не должен превышать допустимого, равного для чугунных отопительных приборов 60 м, для конвекторных систем -100 м;
• Напоры во всех точках тепловой сети и местных систем должны быть избыточными, равными не менее 5 м во избежание подсоса воздуха;
• Напор в подающей линии должен обеспечивать невскипание сетей воды и не должен превышать расчетного для трубопроводов и вспомогательного оборудования и арматуры.
В схемах тепловых сетей крупных городов имеются гидравлические связи между ТЭЦ, обслуживающими соседние тепловые районы. Эти связи позволяют переключать отдельные участки сети на питание от той или иной ТЭЦ, а в отдельных случаях при достаточной пропускной способности гидравлической связи осуществлять и параллельную работу станции на общие тепловые сети. Следует отметить, что тепловая сеть является многокольцевой даже при отсутствии перемычек между соседними магистралями или соседними источниками теплоты. Кольца в этом случае образуются подающими и обратными линиями в сети и соединяющими их насосами и потребителями тепла.
Приведенные данные показывают, насколько сложны современные схемы тепловых сетей, а следовательно, и решение вопросов, связанных с их гидравлическими режимами. В практике проектирования, наладки и эксплуатации возникают следующие основные задачи, требующие расчета гидравлических режимов тепловых сетей:
• Трассировка тепловых сетей, оптимизация их диаметров, выбор мест установки подсосных подстанций;
• Проверка тидравлического режима проектируемой или эксплуатируемой сети при сезонных и суточных изменениях расхода воды;
• Проверка возможностей присоединения к определенным точкам действующей тепловой сети дополнительных тепловых нагрузок и необходимых для этой цели мероприятий;
• Разработка мероприятий по аварийному резервированию действующих и проектируемых сетей;
• Разработка схем автоматизации и защиты тепловых сетей;
• Анализ гидравлических режимов действующих сетей при необходимых в практике эксплуатации переключениях и выбора оптимальною варианта переключений;
• Разработка программы переключений для обнаружения места аварии и его локализации;
• Разработка режимов параллельной работы нескольких источников теплоснабжения на общие тепловые сети.
С точки зрения математического моделирования любая инженерная сеть является в первую очередь направленным математическим графом, в котором места соединения и деления потоков заменены узловыми точками (вершинами графа), а участки трубопроводной сети - ветвями (ребрами графа). Трубопроводные системы характеризуются общностью структуры своих схем, а движение жидкости в них подчиняется общим законам течения и двум сетевым законам Кирхгофа(закон сохранения массового баланса и закон сохранения энергии).
Основной задачей при моделировании трубопроводных систем является определение потокораспределения (т.е. определение расхода на каждом участке) и напоров (давлений) в узловых точках. Определение потокораспределеня основано на решении замкнутых систем нелинейных алгебраических уравнений.
Вопросам решения задачи потокораспределения посвящено большое количество работ.
Для многих гидросистем характерны простые схемы с последовательно-параллельным соединением участков и одним источником напора. К методам расчета таких «цепочечных схем» относятся: метод «эквивалентных отверстий», предложенный в 1873 г. Д. Мюргом [75], метод «перемещения единицы объема», предложенный П.Н. Каменевым [19], а также метод суммирования сопротивлений и проводимостей, впервые примененный Б.Л. Шифринсоном для расчета переменных режимов в тепловых сетях[65, 67]. В алгебраическом смысле эти методы эквивалентны и основаны на свойстве цепочечной схемы, позволяющем сводить ее матрицу к единственной строке (контуру).
Большой цикл работ посвящен графическим и графоаналитическим методам. Среди графических методов различают две группы построений. Первая из них основана на выполнении расчетов путем построения сопрягающих кривых в системе координат «напоры-расходы» [31, 43, 66]. Вторая группа сводится к построению пьезометрических графиков в системе координат «напоры(высоты) -длины участков сетей». Несомненна важная иллюстративная роль этих методов при анализе результатов расчета.
Наиболее важны работы посвященные аналитическим методам расчета, затраг ивающие вопросы количественного регулирования, совместной работы ряда источников на общей сети, режимов открытых систем теплоснабжения и т.п. [3, 4, 5, 8, 11, 13, 14, 22, 23, 27, 28, 56, 70, 68, 69] . Применение аналитических методов позволило получить ряд важных решений, способствовавших созданию рациональных систем теплоснабжения.
Наиболее эффективными методами расчета потокораспределения в кольцевых трубопроводных сетях стали итеративные, гак называемые «увязочные» методы, которые в математическом отношении являются покоординатными релаксационными методами последовательных приближений. «Увязка» сети с заданными сопротивлениями, нагрузками и действующими напорами имеет своей целью найти такие значения расходов на всех участках и давлений в узлах, которые с наперед заданной точностью удовлетворяли бы обоим законам Кирхгофа.
Увязка напоров в кольцах сети осуществляется путем циклического выполнения следующих операций:
1. По данным о нагрузках (расходах у потребителя) выбирается произвольное начальное приближение для расходов на всех участках так, чтобы первый закон Кирхгофа (материальный баланс) соблюдался во всех узлах системы.
2. По этим расходам и данным о гидравлических сопротивлениях вычисляются потери давления (напоры) на участках и суммарные невязки напоров во всех независимых контурах.
3. Для каждого контура ищется так называемый увязочный расход, отвечающий невязке напора.
4. Полученные увязочные расходы поконтурно алгебраически суммируются с расходами, принятыми в начальном приближении или на предыдущей итерации. Новые расходы используются в качестве очередного приближения для следующей итерации (см. п2).
5. Расчет прерывается, когда невязки напоров в любом из контуров перестают превышать заданное значение.
Для вычисления увязочного расхода АХс по данным о невязке напора
S:X? в контуре с М.М. Андрияшевым [2] и В.ГЛобачевым [25] была предложена известная формула
Аналогичная формула, широко цитируемая в зарубежной литературе, была дана и Харди Кроссом [70]. В этой же работе наряду с рассматриваемой здесь поконтурной увязкой напоров им был предложен и метод поузловой увязки расходов.
Увязочные методы, получившие широкое распространение еще в условиях ручного счета, стали и с пользоваться как основной метод первых расчетов потокораспределения на ЭВМ [1, 6, 71, 73, 76 и др.]. В числе первых были работы лаборатории теплофикации ВТИ им. Ф.Э. Дзержинского по применению ЭВМ для анализа гидравлических режимов в сложных тепловых сетях [15, 16, 18].
С появлением в 1958 - 1959 гг. первых универсальных ЭВМ в ВТИ была начата работа по их применению для гидравлических расчетов тепловых сетей. В 1960 году были разработаны алгоритм и программа расчета многокольцевых гидравлических сетей. Программа была составлена для наиболее распространенных в то время ЭВМ «Урал» [18]. Наряду с указанной программой, предусматривающей независимость расхода воды у абонентов от располагаемого напора в сети, в 1961 году была составлена также программа, учитывающая эту зависимость, имеющую место в реальных тепловых сетях. Начиная с 1960 года, применение ЭВМ для исследования гидравлических режимов тепловых сетей и решения, связанных с этим практических задач для тепловых сетей различных городов получило широкое развитие в ВТИ.
Разработка теорий гидравлических цепей и ее применение для расчета гидравлических режимов тепловых сетей с помощью современных ЭВМ развивается в Сибирском энергетическом институте СО АН СССР [61 ].
Большая работа по разработке и совершенствованию программ расчета гидравлических расчетов тепловых сетей проводится в последние годы в институте ВНИПИэнергопром и Карагандинском политехническом институте [21].
В эти же годы сложились основные требования к программным продуктам, моделирующим работу сложных гидросистем:
• Программа расчета на ЭВМ должна быть универсальной, то есть одна и та же программа должна допускать возможность расчета закрытой или открытой тепловой сети любой конфигурации, в одно - или двухлинейном изображении, регуляторами расхода у абонентов или без них, с регуляторами давления в сети и т. д.
• Программа должна допускать расчет сети с большим числом колец, участков, насосов, абонентов, соответствующим реальным сетям.
• Система подготовки исходных данных должна быть достаточно простой. В программе должен производиться контроль за правильной подготовкой и записью исходных данных.
• Расчет должен производиться с заданной степенью точности. Время счета должно быть невелико и должно измеряться минутами.
• Результаты расчета должны печататься в форме, удобной для их обработки и анализа, в частности, в виде пьезометрического графика.
Как отмечалось выше, система нелинейных уравнений, определяющая потокораспределения в сети, решается методом итераций путем определения поправочных расходов АХс для всех независимых контуров сети. Как показал опыт расчета тепловых сетей, особенно в двухлинейном изображении, когда сопротивление участков контура отличается на несколько порядков, серьезной проблемой является обеспечение сходимости процесса итеративного счета. Выяснилось, что определение поправочного расхода АХс по формуле (1), т. е. по методу Лобачева - Кросса, во многих случаях не обеспечивает сходимости процесса, особенно при неплоских схемах, когда один участок входит в 3 кольца и более, а так же когда процессе счета сопротивления некоторых участков, на которых установлены регуляторы расхода или давления, изменяются для поддержания заданных параметров, и появляется второй цикл итерации.
Начали разрабатываться математические методы расчета и приемы их использования, которые с максимальной вероятностью обеспечивают в этих условиях сходимость итерационного процесса.
Вначале предлагалось снижать коэффициент К в формуле (1) для уменьшения АХс [17], затем была предложена формула для определения поправочного контура расхода, использующая модификацию Л. Выханду метода Ньютона- Рафсона [61]: где AfyN) - невязка напора в очередном контуре С, N - номер итерации. Si -гидравлическое сопротивление і - го участка, Хі - расход на і - том участке. ДЛ(. - падение напора на контуре С, Не -действующий напор на контуре С.
Но и применение модификации Л. Выханду метода Ньютона - Рафсона не всегда обеспечивает сходимость вычислительного процесса. Безусловная сходимость обеспечивается для плоских схем, чем дальше схема отходит от плоской, тем хуже сходимость вычислительного процесса[36, Зингер]. В Сибирском энергетическом институте СО АН СССР была проведена большая работа по разработке строгих общих методов расчета потокораспределения [61, 32, 33, 34, 35, 57, 58, 59, 60 и др.], был разработан обобщенный метод контурных расходов [61. 36].
Обобщенный метод контурных расходов (МКР).
Этот метод - аналог известного метода контурных токов для расчета линейных электрических цепей в сочетании с методом Ньютона [9, 61 ] для решения системы нелинейных алгебраических уравнений.
Рассмотрим сеть имеющую п участков (ветвей), m узлов (вершин), к линейно независимых контуров (k— n-m+1). На схеме выделяется некоторое дерево, связывающее все ее m узлов, в результате все участки разобьются на (т-1) участков дерева и к участков, не вошедших в это дерево, которые называются хордами. Каждая хорда замыкает какую-то последовательность участков дерева и однозначно определяет контур, который фиксируется соответствующей строкой матрицы В.
Исходной для МКР является система (4) - (5) относительно вектора расходовл. Л - матрица размером (т—1)х/? соединений т-1 линейно независимых узлов сети и ее п участков, однозначно отображающая структуру сети и ориентацию ее участков; JC=( Xi - jc) " вектор расходов на всех участках; Q=( Q Q) " вектор узловых расходов в линейно независимых узлах; В - матрица совпадений выбранной системы из к линейно независимых контуров и участков сети с размерами fcxn і S - диагональная матрица из величин Si (i-1, ..., п) гидравлических сопротивлений участков; Х- диагональная матрица из абсолютных значений расходов на участках, т.е. из величин х(/=1,...,w) ; .//=( //( •"//„) " вектор действующих напоров; Процесс последовательных приближений имеет вид: (Л -Н) (.V) . ІЛЧ1) , ,\ X = Х + АХ ІР) где порядковый номер итерации. Процесс последовательных приближений осуществляется в МКР таким образом, что для любого X выполняются условия (4) материального баланса в узлах. Для этого достаточно, чтобы этим условиям удовлетворяло начальное приближение Axm = Q О) а все поправки (приращения расходов) Аг=(АХ,...,Ахл7) к расходам имели бы нулевые балансы: AAx=A,Ax + AJIAx,l = 0. (8)
Из (8) следует, что в качестве независимых (свободных) переменных можно взять приращения к расходам па хордах Дд;г, тогда соответствующие
им поправки к расходам на участках дерева Д будут определяться
следующим образом: Поскольку [34 ] то имеем
Ах=(Ах[,&х У=(Е,ВУАх,=В Ах, О О
Приращения к расходам на хордах Дд выбираются из условия уменьшения невязок, которые получаются при подстановке в (5) очередного приближения х }: яяГУ д/Г, (12) где Д/;т - вектор-столбец невязок напоров, вычисленных для контуров по данным N-ro приближения. Для пассивных цепей д/Г=1д/Г. (13) [ Для контуров, содержащих активные участки с источниками напора, д/Г=1д/Г-#о (14) с является алгебраической суммой напоров, действующих на участках, входящих в контур С. Приращения расходов должны обеспечивать выполнение первого закона Кирхгофа (условие материального баланса в узлах), т.е. выполнения равенства BSiX +AX X x AW v- 0. (16) Раскрывая скобки, пренебрегая членами, содержащими Дх и подставляя вместо Ахего выражение (11) через Ах , придем с учетом (12) к системе линейных выражений лГ ЛхГ— л/Г, (і?) которую надо решать на каждом шаге вычислительного процесса. В выражении (17) матрица Д/ Л) размера fcxk (по числу линейно-независимых контуров) есть симметрическая матрица M BSX B1. О») Решая систему уравнений (17) каким-либо из методов линейной алгебры [9, 24, 37, 49, 55] определим приращения расходов на хордах Дх +" а по ним по формулам (9) и (10) вычисляются поправки Дх +1 к расходам на участках дерева. Это позволяет получить новое приближение для вектора (6) расходов. Далее производится анализ сходимости вычислительного процесса: lx(V+, -Jt " ! . (19) Здесь S - вектор из допустимых погрешностей в значениях расходов. По новому приближению х Л+1 опять формируется и решается система уравнений (17) и т.д.
Матрица М является симметрической. Если приравнять нулю все недиагональные элементы, то вместо системы уравнений (17) получим набор формул для «увязочных расходов» в методе Андрияшева - Лобачева -Кросса. При этом безусловная сходимость лишь в случае преобладания диагональных элементов. Это условие выполняется в случаях, когда каждый участок сети входит не более чем в два кон тура, т.е. для плоских схем.
Таким образом, методы с определением «увязочных расходов» по формуле Андрияшева - Лобачева - Кросса (1) или формуле, использующая модификацию Л. Выханду метода Ньютона Рафсона (2), очень выгодны с точки зрения вычислительных затрат, но на реальных схемах со сложной топологией и сопротивлениями участков контура, огличающимися на несколько порядков, во многих случаях не обеспечивает сходимости итерационного процесса, особенно для сетей, на которых установлены регуляторы расхода или давления.
В общем случае система линейных уравнений (17) размера fcxk Д°лжна решаться в полном виде, что влечет за собой вычислительные затраты пропорциональные А на каждой итерации [7].
Усложнение схем тепловых сетей, требование практики о расширении круга решаемых задач, широкое распространение программ расчета гидравлических режимов тепловых сетей в вычислительных центрах энергосистем и других организаций приводит к непрерывному совершенствованию программ. В настоящее время такие работы ведут в основном в следующих направлениях: совершенствование методов счета; упрощение подготовки исходной информации и контроля за её качеством; повышение наглядности результатов расчета и облегчение их использования.
Методы гидравлического расчета были рассмотрены выше. Рассмотрим общие подходы к разработке эффективных сред интерфейса с пользователем, упрощающие подготовку исходной информации, контроль за её качеством: повышающими наглядность результатов расчета и облегчение их использования.
Эффективность сред (оболочек) интерфейса с пользователем во многом определяется выбранной моделью гидравлических сечей, ее структурой (набор составляющих сеть элементов и взаимосвязи между ними). Внутренняя структура модели отражается и на ее внешнем представлении, с которым работает конечный пользователь.
По мере развития аппаратных и программных средств вычислительной техники - менялись внутренняя структура данных моделей и внешнее отображение объектов моделирования.
Программные продукты написанные для ЭВМ «Урал», «БЭСМ», «ЕС», работавшие в пакетном режиме с исходными данными на перфокартах или лентах, не имели интерфейса пользователя в современном понимании, но в тоже время хорошо продуманные формы ввода исходной информации облегчали использование программ расчета инженерам-гидравликам, не являющимися программистами. К таким продуктам относятся например: «АСИГР», «СОСНА» [61, 38J, разработанные в Сибирском энергетическом институте СО АН СССР и др.
Программные продукты, разработанные для персональных компьютеров. имеют более дружественный интерфейс пользователя; современные программы, как правило, имеют графический интерфейс.
В различных системах применялись модели с различным внутренним и внешним представлением в однолинейном [21, 29] и двухлиненейном исполнении [61], с различной степенью детализации элементов, составляющих модель гидросеги.
Разрабатывались системы на основе геоинформационных систем (ГИС) [29] и с собственным графическим представлением схемы [21].
Уровень детализации модели также различен у различных авторов, в зависимости от задач, которые они перед собой ставили.
В представленной работе рассмотрены алгоритмы расчета и математические модели сложных гидравлических систем, ориентированные для реализации на ЭВМ, и разработанное на их основе программное обеспечение моделирования для двух энергетических систем: сетей теплотрасс большого города и многоветочного нефтепровода. Программное обеспечение ориентировано на инженеров разработчиков и инженеров, эксплуатирующих трубопроводные системы.
Материал диссертации расположен следующим образом:
• В главе 1 приведена постановка задачи проектирования системы моделирования сложных гидравлических сетей. Здесь рассмотрены : математическая модель потокораспределения; внутренняя и внешняя модели гидравлической сети для реализации на ЭВМ.
• В главе 2 рассмотрены алгоритмы решения задачи потокораспределения с учетом наличия в сети регуляторов расходов и напоров; рассмотрен метод последовательных приближений решения задачи потокораспределения сложных многокольцевых сетей.
• В главе 3 рассмотрено автоматическое эквивалентирование гидравлических схем, снижающее размерность решаемой системы нелинейных алгебраических уравнений.
• В главе 4 приведены алгоритмы расчета тепловых потерь, теплоотдачи на потребителях тепла и температур теплоносителя на участках тепловой сети.
• В главе 5 рассмотрена внутренняя модель принципиальной гидравлической схемы теплосети; основные структуры и объекты моделирующие элементы теплосети, их иерархия и базы данных модели.
• В главе 6 рассмотрены вопросы моделирования нефтепроводов и поиск оптимальных режимов для них.
• В заключении содержатся основные выводы работы.
Настоящая диссертация выполнена на кафедре «Информационной безопасности» Омского государственного университета под руководством Р.Т. Файзуллина.
Вошедшие в нее результаты опубликованы в работах [10. 12, 26. 64, 72 ] и докладывались на международной конференции «Методы оптимизации и экономические приложения» (Омск, 1997), IV международной научно-технической конференции «Динамика систем, механизмов и машин» (Омск, 2002), XI 1-й всероссийской конференции «Математическое программирование и приложения» (Екатеринбург. 2003).
Разработанные методы и алгоритмы реализованы в программных продуктах - системах моделирования теплосетей и нефтепроводов. Система моделирования тепловых сетей эксплуатируется в АК «ОМСКЭНЕРГО» и в НИУ КТИ ВТ СО РАН; система моделирования потокораспределения и система диспетчерского контроля и управления с выбором оптимальных режимов, эксплуатируется в ОАО «Транснефть».
В заключении введения автор выражает благодарность научному руководителю д.т.н. Р.Т. Файзуллину и всем сотрудникам кафедры «Информационной безопасности» ОмГУ за постоянное внимание к работе.
Основные закономерности установившегося течения жидкости
Матрица М является симметрической. Если приравнять нулю все недиагональные элементы, то вместо системы уравнений (17) получим набор формул для «увязочных расходов» в методе Андрияшева - Лобачева -Кросса. При этом безусловная сходимость лишь в случае преобладания диагональных элементов. Это условие выполняется в случаях, когда каждый участок сети входит не более чем в два кон тура, т.е. для плоских схем.
Таким образом, методы с определением «увязочных расходов» по формуле Андрияшева - Лобачева - Кросса (1) или формуле, использующая модификацию Л. Выханду метода Ньютона Рафсона (2), очень выгодны с точки зрения вычислительных затрат, но на реальных схемах со сложной топологией и сопротивлениями участков контура, огличающимися на несколько порядков, во многих случаях не обеспечивает сходимости итерационного процесса, особенно для сетей, на которых установлены регуляторы расхода или давления.
В общем случае система линейных уравнений (17) размера fcxk Длжна решаться в полном виде, что влечет за собой вычислительные затраты пропорциональные А на каждой итерации [7]. Усложнение схем тепловых сетей, требование практики о расширении круга решаемых задач, широкое распространение программ расчета гидравлических режимов тепловых сетей в вычислительных центрах энергосистем и других организаций приводит к непрерывному совершенствованию программ. В настоящее время такие работы ведут в основном в следующих направлениях: совершенствование методов счета; упрощение подготовки исходной информации и контроля за её качеством; повышение наглядности результатов расчета и облегчение их использования.
Методы гидравлического расчета были рассмотрены выше. Рассмотрим общие подходы к разработке эффективных сред интерфейса с пользователем, упрощающие подготовку исходной информации, контроль за её качеством: повышающими наглядность результатов расчета и облегчение их использования.
Эффективность сред (оболочек) интерфейса с пользователем во многом определяется выбранной моделью гидравлических сечей, ее структурой (набор составляющих сеть элементов и взаимосвязи между ними). Внутренняя структура модели отражается и на ее внешнем представлении, с которым работает конечный пользователь.
По мере развития аппаратных и программных средств вычислительной техники - менялись внутренняя структура данных моделей и внешнее отображение объектов моделирования.
Программные продукты написанные для ЭВМ «Урал», «БЭСМ», «ЕС», работавшие в пакетном режиме с исходными данными на перфокартах или лентах, не имели интерфейса пользователя в современном понимании, но в тоже время хорошо продуманные формы ввода исходной информации облегчали использование программ расчета инженерам-гидравликам, не являющимися программистами. К таким продуктам относятся например: «АСИГР», «СОСНА» [61, 38J, разработанные в Сибирском энергетическом институте СО АН СССР и др. Программные продукты, разработанные для персональных компьютеров. имеют более дружественный интерфейс пользователя; современные программы, как правило, имеют графический интерфейс.
В различных системах применялись модели с различным внутренним и внешним представлением в однолинейном [21, 29] и двухлиненейном исполнении [61], с различной степенью детализации элементов, составляющих модель гидросеги.
Разрабатывались системы на основе геоинформационных систем (ГИС) [29] и с собственным графическим представлением схемы [21].
Уровень детализации модели также различен у различных авторов, в зависимости от задач, которые они перед собой ставили.
В представленной работе рассмотрены алгоритмы расчета и математические модели сложных гидравлических систем, ориентированные для реализации на ЭВМ, и разработанное на их основе программное обеспечение моделирования для двух энергетических систем: сетей теплотрасс большого города и многоветочного нефтепровода. Программное обеспечение ориентировано на инженеров разработчиков и инженеров, эксплуатирующих трубопроводные системы.
Материал диссертации расположен следующим образом: В главе 1 приведена постановка задачи проектирования системы моделирования сложных гидравлических сетей. Здесь рассмотрены : математическая модель потокораспределения; внутренняя и внешняя модели гидравлической сети для реализации на ЭВМ. В главе 2 рассмотрены алгоритмы решения задачи потокораспределения с учетом наличия в сети регуляторов расходов и напоров; рассмотрен метод последовательных приближений решения задачи потокораспределения сложных многокольцевых сетей. В главе 3 рассмотрено автоматическое эквивалентирование гидравлических схем, снижающее размерность решаемой системы нелинейных алгебраических уравнений. В главе 4 приведены алгоритмы расчета тепловых потерь, теплоотдачи на потребителях тепла и температур теплоносителя на участках тепловой сети.
Решение системы уравнений (2.1) на каждом шаге итерации
Для расчета потокораспределения в сетях, оснащенных регуляторами расхода, давления и обратными клапанами в процедуре решения уравнений потокораспре деления (1.22) предусмотрены действия по изменению (вычислению) гидравлического сопротивления участков, соответствующих регуляторам. Эти действия выполняются при выполнении определенных условий, например при достижении невязки величины, которая меньше заранее установленной величины. После получения вектора решений (расходов) на итерации уточняется сопротивления регулирующих элементов.
Обратный клапан - это регулятор пропускающий расход жидкости только в одном направлении, если жидкость течет в обратном направлении, то обратный клапан закрывается. Типичная схема включения обратного клапана показана на Рис. 2.8.
Схема включения обратного клапана. Закрытие обратного клапана имитируется увеличением его сопротивления до величины тах=1Е+06. Однако следует иметь в виду, что скачкообразное изменение сопротивления регулятора с номинального сопротивления до максимального на больших схемах со сложной топологией приводит к нарушению сходимости вычислительного процесса. Поэтому в предлагаемом алгоритме сопротивление при закрытии обратного клапана увеличивается постепенно (в 100 раз на каждом шаге итерации). Регулирование заканчивается когда сопротивление достигает величины 5тях это соответствует полностью закрытому клапану. В случае если направление течения жидкости правильное, сопротивлению присваивается величина Smin.
Теперь рассмотрим алгоритм вычисления нового сопротивления. Пусть Sm-,a - сопротивление полностью открытого клапана. S\,ias -сопротивление полностью закрытого клапана, Scur - текущее сопротивление в процессе регулирования (первоначально равно Smia), X — расход через клапан, RegFlag - состояние процесса регулирования.
Следующие действия проводятся на каждой итерации после получения решения (расходов). Код дан на языке С. Вычисление сопротивлений регуляторов расходов. Вычисление нового сопротивления регуляторов расхода проводится после получения решения (расходов) на каждой итерации, где значения невязки напоров меньше 10 м.
Пусть Smin - сопротивление регулятора при полностью открытом клапане, Smax - сопротивление реіулятора при полностью закрытом клапане, Snew -новое вычисленное сопротивление, Scur - текущее сопротивление с предыдущей итерации (первоначально равно Smin), X - расход через регулятор, Q - расход, который должен поддерживать регулятор, RegFIag -состояние процесса регулирования, Nr - заданная допустимая невязка по расходам.
Эквивалентирование сетей, связанных с магистральным трубопроводом одной врезкой
Квартальные сети теплосетей, как правило, присоединяются к магистральному трубопроводу только одной врезкой (Рис. 3.3), т.е. образуют ответвление, которое можно назвать «листом». Общепринято такие «листы» заменять обобщенным потребителем, с эквивалентным сопротивлением. которое равно сопротивлению «листа». Такой подход укрупняет схему, а следовательно, уменьшает размерность решаемой системы нелинейных уравнений и улучшает сходимость вычислительного процесса. После решения задачи иотокораспределения и определения напоров в узлах определяется разность напоров на врезке, что позволяет затем решить задачу потокораспределения внутри квартальной сети («листа»).
Отметим, что такой подход правомерен, если внутри листа отсутствуют источники напоров и регуляторы расхода и напора, т.е. элементы с изменяющимся сопротивлением.
Кроме того, такая замена абсолютна эквивалентна, если внутри листа отсутствует открытый водоразбор. При наличии открытого водоразбора, не превышающего 10-20% от циркуляционного расхода, погрешность решения незначительна и данное эквивалентирование можно применять с приемлемой для практики эксплуатации точностью.
Рассмотри алгоритм, решения задачи потокораспределения, использующий данный способ эквивалентирования (термины: обобщенный потребитель; «лист»; квартальная сель - в данном контексте синонимы). Шаг 1. Определение эквивалентного сопротивления «листа» (Рис. 3.4). К двум узлам врезки присоединяется фиктивный источник напора с некоторым действующим напором Нд. В узле врезки, инцидентном с обратным трубопроводом, устанавливается некоторый произвольный напор ho. При наличии открытого водоразбора весь водоразбор «листа» присваивается данному узлу (подпитка). Составляется расчегная схема «листа» сети, которая включает в себя все участки «листа» и фиктивный источник напора. Решается задача потокораспределения «листа», в результате получаем расход X через фиктивный источник, который будет равен расходу через врезку, соединяющую «лист» с магистралью. Падение напора па врезке Ah соответственно равно действующему напору Нл фиктивного источника напора. Зная расход через врезку X и падение напора на врезке Д/ї определяем эквивалентное сопротивление обобщенного потребителя « А/? («листа») Ь=— , которое будет использовано при расчете X2 магистральной сети. Шаг 2. Решается задача потокораспределения и определения напоров в узлах для магистральной сети, с обобщенными потребителями, сопротивление которых найдено при выполнении шага 1. Для магистральной сети - это окончательное решение. Шаг 3. Решение задачи потокораспределения и определения напоров в узлах для сети «листа» (квартальной сети). После расчета магистральной сети мы имеем значения падения напоров Д# на всех «листах», присоединенных к магистральному трубопроводу, т.е. напоры в узлах всех интересующих нас врезках. 7! Фиктивным источникам присваиваются уже не произвольные действующие напоры, а реальные Н — /S.h, Составляется расчетная схема «листа» сети, которая включает в себя все участки «листа» и фиктивный источник напора. Окончательно решается задача потокораспределения «листа» и определения напоров в узлах. 72 Практическая реализация вышеизложенного алгоритма предполагает следующие решения: принципиальная схема «листа» физически находится в отдельном файле (или в том же файле, но в отдельном слое), что облегчает выделение его участков (ребер графа) из графа общей схемы; на врезке со стороны магистрального трубопровода и «листа» вводятся в схему специальные (фиктивные) элементы, в которых содержится информация о связи «листа» с магистралью (Рис. 3.4). Задача наладки тепловых сетей - это задача нахождения потокораспределения и напоров в узлах при заданных напорах на источниках и заданных расходах на потребителях, и подбора сопротивлений дросселирующих устройств (шайб и сопел элеваторов) у потребителей, которые обеспечили бы распределение расходов теплоносителя между потребителями в точном соответствии с заданными нагрузками. Иначе эту задачу еще называют - расчет с фиксированными нагрузками потребителей, в отличие от расчета с фиксированными сопротивлениями потребителей.
Одним из возможных решений является включение в схему вместо потребителей «фиктивных» регуляторов расхода, настроенных на требуемый расход потребителей [61]. Получаемые в результате расчета сопротивления этих регуляторов оказываются равными тем, которые требуются потребителям. По полученным сопротивлениям рассчитываются диаметры шайб или сопел элеваторов.
Расчет температуры теплоносителя
В системе уравнений должен соблюдаться баланс притоков и оттоков по узлам гидросети, поэтому потокораспределение однозначно определено при наличии одной точка подпитки, которая одновременно является точкой с заданным напором, от которого рассчитываются напоры во всех остальных узлах при известных расходах на участках. На эту точку подпитки собирается весь водоразбор, что обеспечивает баланс притоков и отгоков.
Задача расчета потокораспределения и в частности распределения подпиточных расходов для сети с несколькими точками подпитки и несколькими точками с фиксированным напором решается эквивалентированием гидросхемы. Рассмотрим схему, включающую в себя 2 источника напора и 2-е точки подпитки (Рис. 3.7). При этом единым режимом оказываются связанными и коммуникации подпиточной воды обоих источников напора.
В работе [61 ] предложен метод эквивалентирования схемы. Принципиальная расчетная схема совместной работы источников при заданных давлениях на обратных коллекторах приведена на Рис. 3,8.
Как видно из схемы вводятся дополнительные источники напоров (подпорные насосы), соединенные общей точкой подпитки и регуляторы давления «после себя». Для задания параметров участков 1-2 и 1-3 не обязательно иметь полную информацию о всех эквивалентируемых ими коммуникациях. Достаточно задать где h2 и A j - напоры, поддерживаемые регуляторами под низки; h/ -пьезометрическая отметка водоема. Полученные затем в результате расчета сопротивления S]_2 и Si-з автоматически обеспечат правильное эквивалентирование схемы. Однако данный метод, вводя дополнительные регуляторы, увеличивает число итераций вычислительного процесса и вносит погрешность в результат, так как напоры, заданные на регуляторах, при расчете обеспечиваются с некоторой заданной точностью. В данной работе предлагается способ эквивалентирования описанный ниже. Во внешней модели схемы (т.е. то с чем работает пользователь) имеется специальный элемент, размещением которого на схеме указывается точка подпитки с заданным давлением, на Рис. 3.7 это элементы Ш и П2. При составлении внутреннего представления расчетной схемы выполняются следующие шаги: каждый такой элемент заменяется парой «источник напора - обратный клапан», начальная точка одной из точек подпитки соединяется с начальными точками всех остальных точек подпитки, участком с действующим напором равным разности геодезических отметок между этими точками. Обратный клапан включается, чтобы исключить обратный ток воды, что фактически соответствует реальной практике включения насосов в сеть. Предлагаемая принципиальная расчетная схема совместной работы источников при заданных давлениях на обратных коллекторах приведена на Рис. 3.9.
