CC BY
23
shows that in flows with different resistance laws various «mechanisms» of turbulence take place. The content of the article makes it possible to expand the knowledge and understanding of some hydraulics fundamentals by specialists.
Dissipation, entropy, minimum, maximum, transition, resistance.
References
1. Klimontovich Yu.L. Turbulentnoe dvizhenie i struktura haosa. - М.: Nauka, 1990. - 320 s.
2. Zegzhda A.P. Gidravlicheskie poteri na trenie v kanalah i truboprovodah. - Ь.М.: Gosstrojizdat. 1957, 278 s.
3. Lyutikov A.V. Propusknaya sposob-nost i ustojchivost rusel krupnyh kanalov pri malyh uklonah. Avtoreferat cand. Diss. VNIIGiM. 1987.
4. Verbitskiy V.S., Hodzinskaya A.G. Uchet form dvizheniya nanosov pri opredele-
nii gidravlicheskih soprotivlenij ruslovyh potokov. // Gidrotehnicheskoe stroiteljstvo. -2005. - № 8. - S. 46-50.
The material was received at the editorial office
15.08.2017
Information about the author Vladimir Saulovich Verbitskiy, senior researcher FSRBI VNIIGiM named after A.N. Kostyakov; Moscow, ul. Bolj-shaya Academicheskaya, d. 44, cor. 2, tel.: 8(499)9770816, e-mail: asher5@mail.ru
УДК 502/504:532.5:626.4
Ю.Г. БУРКОВА, В.А. ФАРТУКОВ
Федеральное государственное бюджетное образовательное учреждение высшего образования
«Российский государственный аграрный университет — МСХА имени К. А. Тимирязева», г. Москва, Российская Федерация
Л.Б. БЕКИШЕВА
Акционерное общество «Концерн «Моринсис — Агат», г. Москва, Российская Федерация
ИМИТАЦИЯ ОТКАЗОВ ФУНКЦИОНИРОВАНИЯ НАПОРНОЙ ГИДРАВЛИЧЕСКОЙ СИСТЕМЫ
Представлены результаты имитационного моделирования работы напорной гидравлической системы, представляющей собой разветвленную (кольцевую) схему водоподачи. В основе имитационной модели лежит анализ многомерной случайной величины, базирующийся на методах теории временных рядов. Полученные результаты имитационного моделирования показали хорошее согласование с натурными экспериментами напорной гидравлической системы (НГС). При проектировании и реконструкции НГС необходимо учитывать фактор риска. Учет этого фактора может существенно снизить стоимость всей системы в целом. Имитационное моделирование НГС позволяет получить интервальные оценки надежности с учетом влияния отключения участков, снижения напоров в системе в расчетный час и часы максимального водопотребления, числа отказов за период эксплуатации. Полученные результаты могут служить основой оптимизации как вновь создаваемой, так и реконструируемой НГС с целью обеспечения заданной надежности подачи воды.
Напорная гидравлическая система, имитационное моделирование, факторный анализ, надежность водоподачи, оптимизация работы напорной гидравлической системы.
Введение. Процесс функционирования напорной гидравлической системы (НГС) можно представить в виде смены различных состояний, характеризуемых величиной водопотребления и техническим состоянием входящих в нее элементов. При этом имеющая место стохастичность ха-
рактеризуется не регулированным случайным процессом отбора воды потребителями, процессами износа и старения, аварийностью, отказами работы оборудования.
При моделировании гидравлических режимов водопроводных сетей в условиях аварийных отключений отдельных участков
необходимо учитывать сокращение водопо-требления, вызванное снижением напоров у ряда абонентов.
Целью работы является определение метода решения поставленной проблемы и частных задач работы и проектирования напорной гидравлической системы на основе проведения имитационного моделирования, основанного на теории надежности и учета влияния отключения каждого участка системы.
В соответствии с оценками, приведенными в [5] можно заключить, что при напорах у потребителей Н1, равных или больших требуемого И. треб, водопотребление 1 — го узла Q1 равно нормативному значению Q1 треб, а при снижении напора Н водопотребление уменьшается по закону:
=
<2^ треб ,
QIшреб ■ ^н~Тн~п О
при Н 1 > Н 1 треб ПриН 1 крыт < Н 1 < Н 1 треб
При Н 1 < Н 1 чрит
где И. крит — критический напор, соответствующий полному прекращению водоразбора.
Метод решения. При такой постановке задачи, обычные методы контурной увязки работают неустойчиво, что приводит к тому, что итерационный процесс расходится. Для устранения этого недостатка в аварийных ситуациях (при отключении отдельных участков), НГС моделируется системой линейных и нелинейных уравнений [3]:
^рп+(2т(Нг)=0 '=¡-т;
+ Н1)-+ Н])-■ Я 1} ■ Р = 0 = 1-т^ ¡-1ап -1 зд;
Hj + К- -нщ + —2■ qвj ■ =0 }=1-,к;
Н,2 - Нг1 - Н ст + -у^ qV ' Щ \ = 0 У ^1ап;
Нц -Нг-
g-ж-й
- qiy ■ =0 >Уе¡зд■
Решение этой системы уравнений осуществляется методом Ньютона. При моделировании случая отключения участка НГС от водопотребления, диаметр отключаемого участка при расчете снижается до величины, близкой к нулю, или увеличивается его гидравлическое сопротивление, что приводит к уменьшению пропускной способности моделируемого участка практически до нуля.
В качестве имитируемых были выбраны случайные факторы, оказывающие наибольшее воздействие на систему: — вектор узловых расходов водопотребления; - отказы и восстановления участков водопроводной сети. Узловые расходы имитируются с шагом квантования в 60 минут. По полученным результатам имитационного моделирования можно оценить возможность системы к саморегулированию, её резервы и "узкие места", скорректировать отдельные параметры и повторить имитацию.
Напорная гидравлическая система является восстанавливаемым объектом, работоспособность которого, в случае возникновения отказа, подлежит восстановлению. Рассмотрим процесс отказа и восстановления отдельных участков НГС. Процесс функционирования восстанавливаемого элемента представляет собой чередование периодов: исправной работы, отказа и восстановления, после чего снова наступает период исправной работы [4]. Продолжительность ремонта зависит от ряда операций, определяемых различными причинами, и считается случайной величиной.
На рисунке 1 приведен фактический временной ряд водопотребления жилого дома. Из него видно, что отбор воды имеет вид случайных периодических колебаний. В отдельные моменты времени он достигает нулевых значений, что соответствует полным отказным состояниям.
Процесс восстановления продолжается несколько часов, в течение которых водопотре-бление постоянно изменяется в соответствии с суточными колебаниями и стохастичностью процесса. Таким образом, при отключении участков в одних случаях будет происходить отказ узлов из-за нехватки напоров в них, в других случаях при уменьшении водопотре-бления отказа узлов не будет. Для уменьшения негативных явлений в послеаварийный период (до включения участка) можно управлять параметрами водопотребителей, например, путем включения дополнительных насосов или частотным регулированием оборотов электродвигателя насосной установки.
Исследования показали [1], [2] и др., что как процесс отказов, так и процесс восстановлений участков, подчиняются экспоненциальному закону распределения. В работе использовались величины удельных интенсивностей отказов и восстановлений, как функции диаметров и материалов труб, приведенные в таблице:
Таблица
Основные свойства случайных величин для имитации НГС [5]
Распределение М^ Функция распределения Плотность вероятности Генерирование случайных чисел
Экспоненциальное 1 A 1 Л2 1-е- М lxe-At -1 11 = ~Т 1п п л
Нормальное m с2 1 ге 1 2 (t-m)2 с42л e 2с2 Л к Ъ Г - 2 11 = <г—-- + т (к ту к > 12
Вейбулла ЬГ[(1+с)/с] Ь2Г[(с+2)/с]-{Г[(с+1)/с]}2 -(г 1 ь 1 - е -(b )с с^с-1/ Ьсе 11=Ь(-1п п)1/с
В результате приходим к следующей последовательности расчета:
- в начальный момент времени все участки считаются исправными. Формируются массивы времени отказа и восстановления каждого участка. Определяется время ближайшего отказа, которое переводится в календарную дату и час суток. По ним находятся номер графика водопотребления и расчетный час. По другой программной реализации алгоритма, выбирается наиболее неблагоприятный случай — час с максимальным в течение суток расходом.
Осуществляется имитация случайного вектора водопотребления и производится гидравлический расчет системы при отключенном аварийном участке по схеме с нефиксированными отборами. Определяется возможность включения дополнительных насосов и частотного регулирования.
Определяются узлы с напорами, меньше требуемых величин, а результаты сохраняются в файле вывода. Производится переход к следующему часу суток и сравнение времени после аварии с заданным временем восстановления. В случае продолжения ремонта для нового часа разыгрывается вектор водопотребле-ния и производится увязка сети. Расчет повторяется до истечения времени ликвидации ава-
рии, после чего участок считается работающим и для него определяются новое время отказа и восстановления. Моделируется время ближайшей аварии и расчет повторяется до окончания имитации. Производится обработка полученных данных, в результате чего для каждого узла определяется суммарное время отказа за период имитации и вклад аварий отдельных участков в отказ выбранного узла.
При разработке алгоритма были сделаны следующие допущения:
- отказы участков независимы между собой, при имитации узловых расходов не учитывалось изменение характера водо-потребления при аварии, связанное с перераспределением времени разбора воды. Узловые расходы подчиняются нормальному закону распределения.
В качестве примера рассмотрим кольцевую систему подачи и распределения воды с контррезервуаром. По существующим методикам строился ступенчатый график водопотребления. Часом максимального водопотребления является час (21-22). В утренние часы максимум разбора приходится на час (7-8). На них мы и будем ориентироваться при оценке надежности. При аварии узел с минимальным напором может сместиться и ориентироваться необходимо на узлы,
в которых снижение напора наиболее критично (например, крупные предприятия с непрерывным технологическим процессом). Кольцевая структура сети повышает ее надежность, но не гарантирует отсутствия отказов. На первом этапе проводился детерминированный расчет сети, и определялись параметры критического узла в час макси-
мального водопотребления при отключении каждого участка. Величина свободного напора в критическом узле в зависимости от отключаемого участка приведена на рисунке 2. Подобные расчеты могут служить основой для оптимизации системы с целью обеспечения снижения напоров в отдельных узлах, не ниже заданной величины.
Аварийный напор в узле 9 от отказавших участков в Час Qmax
Аварийный напор в узле 9 от отказавших участков в Час Отах
10 12 14 Участок
Рис. 2. Влияние отключения участков на напор в час Qmax (21-22), Qтреб = 14м.
(Гидравлическая имитация)
Для этой системы был проведен имитационный расчет по алгоритму имитационного моделирования [3] в течение 20 лет. В результате определялось число часов отказов каждого узла, а также число отказов критического узла в зависимости от отключаемого участка за весь период расчета, который принят равному нормативному времени эксплуатации. Результаты фикси-
ровались либо в час аварии с последующим изменением водопотребления за период ремонта, либо аварии приводились к часу максимального водопотребления по принципу наименее благоприятного случая. Результаты приведены на рисунке 3. Анализ подобных графиков может быть основанием для корректировки параметров системы и режимов ее эксплуатации.
Число часов отказа каждого узла за период расчёта
Рис. 3. Число часов отказов узлов за 20 лет эксплуатации (расчет в часы аварии).
(Гидравлическая имитация)
Выводы
Результаты имитационного моделирования хорошо согласуются с натурными экспериментами НГС, при анализе, как с помощью многомерной случайной величины, так и методов теории временных рядов.
При проектировании и реконструкции НГС необходимо учитывать фактор риска, что может существенно снизить стоимость системы.
На основе теории надежности, гидравлических расчетов и имитационного моделирования для НГС получены оценки:
- влияния отключения каждого участка на уменьшение напора выбранного узла в заданный час;
- числа часов отказов каждого узла за период эксплуатации;
- числа отказов выбранного узла за период эксплуатации;
- проведение расчетов в час аварии и час максимального водопотребления, что позволяет получить интервальные оценки надежности.
Полученные результаты могут служить основой оптимизации как вновь создаваемой, так и реконструируемой НГС с целью обеспечения заданной надежности подачи воды.
Имитационное моделирование является эффективным методом расчета динамики функционирования гидравлических систем, включающих насосные станции и напор-но-регулирующие емкости.
При проектировании гидравлической системы необходимо учитывать выбранную схему эксплуатации и процесс управления.
Библиографический список
1. Абрамов Н.Н. Надежность систем водоснабжения: (Надежность и качество) — М.: Стройиздат, 1979. - 231 с.
2. Эксплуатация систем водоснабжения, канализации и газоснабжения: Справочник / В.Д. Дмитриев, Д.А. Коровин, А.И. Кораблев и др. - Л.: Стройиздат, 1988. - 383с.; ил.
3. Карамбиров С. Н. Математическое моделирование систем подачи и распределения воды в условиях многорежимности и неопределенности: Монография. — М.: МГУП, 2004. - 197с.
4. Карамбиров С.Н., Манукьян Д. А., Бекишева Л.Б. Оценка надежности подачи воды системами водоснабжения. / Доклады Российской Академии сельскохозяйственных наук. - 2013. - № 6. - С. 63-65.
5. Эгильский И. С. Автоматизированные системы управления технологическими процессами подачи и распределения воды. -Л.: Стройиздат, 1988. - 216 с.: ил.
Материал поступил в редакцию 22.07.2017 г.
Сведения об авторах
Буркова Юлия Геннадьевна, кандидат технических наук, доцент. ФГБОУ ВО РГАУ-МСХА имени К.А. Тимирязева; 127550, г. Москва, Б. Академическая ул., д. 44; тел.: 8(499)1539766; e-mail: burkova.msuee@mail.ru
Фартуков Василий Александрович, кандидат технических наук, доцент. ФГБОУ ВО РГАУ-МСХА имени К.А. Тимирязева; 127550, г. Москва, Прянишникова ул., д. 19; тел.: 8(916)6531759, e-mail: vasfar@mail.ru
Бекишева Лаура Борисовна, аспирантка, менеджер АО «Концерн «Моринсис - Агат», 105275, г. Москва, шоссе Энтузиастов, 29; тел.: 8(917)5425746; e-mail: laura.lora@list.ru
YU.G. BURKOVA, V.A FARTUKOV
Federal State Budget Educational Institution of Higher Education "Russian State Agrarian University-MAA named after C.A. Timiryazev", Moscow, Russian Federation
L.B. BEKISHEVA
JSC "Concern "Sea information systems — Agat", Moscow, Russian Federation
SIMULATION OF FAILURES OF PRESSURE HYDRAULIC SYSTEM FUNCTIONING
The paper presents the results of simulation modeling of the operation of the pressure hydraulic system which is a branched (ring) water supply scheme. The simulation model is based on the analysis of a multidimensional random variable based on the methods of the theory of time series. The obtained results of the simulation modeling showed a good agreement with the field experiments of the pressure hydraulic system (PHS). When designing and reconstructing PHSs, it is necessary to take into account the risk factor. Accounting for this factor can significantly reduce the cost of the entire system as a whole. Simulation modeling of PHS allows obtaining interval
reliability estimates taking into consideration the effect of disconnection of sections, reduction of pressure in the system at the calculated hour and hours of maximum water consumption, number of failures during the operation period. The obtained results can serve as a basis for optimization of both newly created and reconstructed PHS in order to ensure a given reliability of water supply.
Рressure hydraulic system, simulation modeling, factor analysis, reliability of water supply, optimization of pressure hydraulic system.
References
1. Abramov N.N. Nadezhnostj system vo-dosnabzheniya: (Nadezhnostj i kachestvo) -M.: Strojizdat, 1979. - 231 s.
2. Expluatatsiya system vodosnabzheni-ya, kanalizatsii I gazosnabzheniya: Spravoch-nik / V.D. Dmitriev, D.A. Korovin, A.I. Kora-blev I dr. - L: Stroyizdat. 1988-383; il.
3. Karambirov S.N. Matematicheskoe mod-elirovanie system podachi I raspredeleniya vody v usloviyah mnogorezhimnosti i neopredelennos-ti: Monografiya. - M.: MGUP, 2004-197s.
4. Karambirov S.N., Manukyan D.A., Be-kisheva L.B. Otsenka nadezhnosti podachi vody sistemami vodosnabzheniya. // Doklady Rossi-jskoj Akademii seljskohozyajstvennyh nauk. -2013. - № 6. - S. 63-65.
5. Egilsky I. S Avtomatizirovannye siste-my upravleniya tehnologicheskimi protsessa-mi podachi i raspredeleniya vody. - L .: Stroyizdat, 1988. - 216 s .: ill.
The material was received at the editorial office
22.07.2017
Information about the authors Fartukov Vasilij Alexandrovich, associate professor, candidate of technical sciences, Federal state budgetary educational institution of higher education "Russian state agrarian University-MAA named after C.A. Timiryazev", Moscow, 125550, ul. Pryanishnikova, d. 19 tel.: 8(916)6531759; e-mail: vasfar@mail.ru
Burkova Yulia Gennadjevna, associate professor, candidate of technical sciences, Federal state budgetary educational institution of higher education "Russian state agrarian University-MAA named after C.A. Timiryazev", Moscow, 125550, ul. Pryanishnikova, d. 19 tel.: 8(499)1539766; e-mail: burkova.msuee@mail.ru
Bekisheva Laura Borisovna, post graduate student, manager of JSC "Concern "Sea information system - Agat", Moscow, 105275, Enthusiasts highway, 29 tel.: 8(917)5425746; e-mail: laura.lora@list.ru
УДК 502/504:627.8 Б01 10.26897/1997-6011/2018-1-18-24
В.А. ВОЛОСУХИН, Е.Н. БЕЛОКОНЕВ
Институт безопасности гидротехнических сооружений (ИБГТС), Ростовская область, г. Новочеркасск, Российская Федерация
М.А. ВОЛЫНОВ
Государственное бюджетное научное учреждение Всероссийский научно-исследовательский институт гидротехники и мелиорации им. А.Н. Костякова (ВНИИГиМ), г. Москва, Российская Федерация
ВОДОЗАБОРНЫЕ СООРУЖЕНИЯ ГАЭС
Целью исследований являлось экспериментальное подтверждение того факта, что интенсивность воронкообразования в усовершенствованном варианте водоприемника и входного оголовка существенно ниже, чем в проектном и разработка рекомендаций по внесению изменений в проект. В статье изложены результаты физического моделирования водоприемника и входного - выходного оголовка верхнего бассейна (ВБ) второй очереди Загорской гидроаккумулирующей станции ЗаГАЭС-2. Исследовались проектный и реконструированный варианты водоприемника и оголовка. Реконструкция водоприемника и оголовка была предложена с целью сведения до минимума интенсивности воронкообразования в акватории сооружений по результатам выполненного ранее 2Б и 3Б моделирования движения водного потока. В новой конструкции верхний бассейн переходит в водоприемник с уклоном дна 1:7. Перед входным оголовком дно водоприемника запроектировано с нулевым уклоном. Входные отверстия имеют раструбное очертание с вертикальными стенками, достигающими верха входного оголовка. Бычки короче по длине,
