CC BY
64
УДК 528.482.3
ПОСТРОЕНИЕ ПРОГНОЗНОЙ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ПРОЦЕССА ПЕРЕМЕЩЕНИЙ ПЛОТИНЫ САЯНО-ШУШЕНСКОЙ ГЭС
Наталья Николаевна Кобелева
Сибирский государственный университет геосистем и технологий, 630108, Россия, г. Новосибирск, ул. Плахотного, 10, старший преподаватель кафедры физической геодезии и дистанционного зондирования, тел. (383)343-29-11, e-mail: n.n.kobeleva@mail.ru
Валерий Степанович Хорошилов
Сибирский государственный университет геосистем и технологий, 630108, Россия, г. Новосибирск, ул. Плахотного, 10, доктор технических наук, профессор кафедры физической геодезии и дистанционного зондирования, тел. (383)343-29-11, e-mail: Khoroshilovvs@mail.ru
В статье рассматривается подход к построению прогнозной математической модели для изучения процесса перемещений контролируемых точек гребня плотины Саяно-Шушенской ГЭС с целью изучения качественных изменений наблюдаемых деформационных процессов сооружения. Показана последовательность выполнения этапов оценивания при построении прогнозной модели в зависимости от характера прогностической задачи.
Ключевые слова: геодезические данные, прогнозная математическая модель, прогнозирование, перемещения контролируемых точек, деформации сооружения.
BUILD A PREDICTIVE MATHEMATICAL MODEL OF MOTION OF THE DAM SAYANO-SHUSHENSKAYA HPP
Natalia N. Kobeleva
Siberian State University of Geosystems and Technologies, 630108, Russia, Novosibirsk, 10 Plakhotnogo St., a senior lecturer in physical geodesy and remote sensing, tel. (383)343-29-11, e-mail: n.n.kobeleva@mail.ru
Valery S. Khoroshilov
Siberian State University of Geosystems and Technologies, 630108, Russia, Novosibirsk, 10 Plakhotnogo St., Ph. D., Prof. of Department of physical geodesy and remote sensing, tel. (383)343-29-11, e-mail: khoroshilovvs@mail.ru
In article approach to creation of expected mathematical model for studying of process of movements of controlled points of a crest of a dam of Sayano-Shushenskaya hydroelectric power station for the purpose of studying of high-quality changes of observed deformation processes of a construction is considered. The sequence of performance of stages of estimation at creation of expected model depending on character of a predictive task is shown.
Key words: geodetic data, expected mathematical model, forecasting, movements of controlled points, deformations of a construction.
Разработка прогнозных математических моделей является обязательным условием диагностического контроля плотины СШГЭС в соответствии с принятой «Методикой определения критериев безопасности гидротехнических со-
оружении» [1] и рекомендациями Федерального закона [2]. Существенный вклад в разработку прогнозных математических моделей отражен в работах [4,5,7], методологический подход в [3,6,8,9].
Осенью 1990 г. водохранилище Саяно-Шушенской ГЭС было впервые наполнено до проектной отметки (НПУ = 540 м) и с 1991 года ГЭС и её сооружения вышли на проектный эксплуатационный режим [7]. В этом режиме для плотины характерно циклическое изменение основной нагрузки - гидростатического давления верхнего бьефа (УВБ) на верховую грань плотины и дно водохранилища. Каждый цикл имеет две ветви - ветвь нагружения (наполнение водохранилища: май-октябрь); ветвь разгрузки (сработка водохранилища: ноябрь-апрель).
В соответствии с принятыми положениями [1] в качестве одного из диагностических показателей контроля Саяно-Шушенской ГЭС выбраны радиальные перемещения гребней плотины: левобережной секции 18, ключевой секции 33 и правобережной секции 45, измеренные отвесами и соответствующие им значения параметров УВБ и температура бетона Тниз ,ТверХ в базовых точках
[7], а также представлены прогнозные математические модели регрессионного типа для различных периодов работы плотины: 2004-2007; 2007-2009 и 20092014 годы.
В качестве исходных данных для построения динамической модели были использованы данные измерений за 2004-2007 годы, полученные на основании графика (рис.1) из работы [7] (необходимые данные натурных перемещений гребня плотины отсутствуют в открытой печати). В связи с этим, авторы не претендуют на точность прогнозной математической модели, скорее больший интерес представляет методический подход к её построению.
Рис. 1. Изменение диагностического показателя (радиального перемещения гребня плотины) и параметров-аргументов: УВБ и температура бетона в базо-
Т
вой точке низ во времени: — верхний бьеф; — температура бетона (секция
33); — радиальные перемещения.
Одним из вариантов для построения динамической модели может служить рекуррентное уравнение 1-го порядка процесса перемещения, происходящего под воздействием одного из основных факторов, а остаточная часть процесса описывается дополнительно моделью шумовой компоненты, т.е. путем расширения вектора состояний [10-12].
Стационарную динамическую модель представим в виде выражения:
хк = фхк-\ +(3ик + Ск ■ (1)
Заменим в выражении (1) значения х^ и и^ их центрированными по времени значениями х&= х^ - х и и^ - и (х и и - средние величины перемещения и входного параметра: уровень верхнего бьефа на годичном интервале, взятом за период основания прогноза). Это позволяет не только определить &, как будет показано ниже, но и упрощает вычисления, повышая тем самым степень обусловленности системы нормальных уравнений, решаемых для оценивания параметров.
При построении модели в качестве основного входного фактора использовался уровень верхнего бьефа (гидростатическое давление) за период с января 2004 г. по январь 2005 г. (табл. 1). Все исходные данные (уровень верхнего бьефа и результаты перемещений гребня плотины) получены на середину каждого месяца.
Таблица 1
Исходные данные для построения прогнозной модели
Исходные данные математической Остаточные
модели ошибки
Дата, УВБ, м хвых, 8
месяц мм
1 2 3 4
Янв. 2004 528,8 128,15 -
февраль 523,6 119,26 +5,917
март 516,4 100,74 +0,163
апрель 505,0 76,67 -2,720
май 503,6 66,67 -7,435
июнь 532,2 120,00 +7,242
июль 537,2 125,18 -6,798
август 538,4 125,91 -9,007
сентябрь 538,6 127,41 -7,968
октябрь 539,2 135,92 -0,688
ноябрь 537,0 138,52 +3,347
декабрь 533,6 134,44 +3,801
Янв. 2005 528,6 127,04 +4,806
На первом этапе построения динамических моделей методом МНК оценивались параметры по результатам наблюдений за входом |и^ | и выходом |х^ | на периоде основания прогноза £=1,2,...,N. Для этого минимизировался функ-
N ( Л2
ционал ^ (ф,Р)= X - ; через обозначено условное мате-
^ ' у
матическое ожидание уравнения (1) определяющееся выражением:
М|Хк 'Хк-1 >ик1 = *к/к-1 = Vхк-1 + Рик ■ (2)
Оценки параметров находили из решения полученных систем нормальных уравнений.
6041,2816(3 + 1695,98493 = 3873,7474;
Модель: (3)
1695,9849(3 + 1838,75043 = 3133,9478. Из решения системы (3) были получены параметры:
Л Л Л
(3 = 0,2196; у? = 1,5018. После вычисления оценок (,/? находили оценку Х^ из
/V
выражения Х^ = Х (1-()- ; оно получилось равным Х^ =-701,141мм.
Остаточные ошибки (табл.1) характеризуют, с одной стороны: - корректность построения модели с точки зрения структурной идентификации; с другой, - это характеристика свойств шумовой компоненты , т.е. используя остаточные ошибки £^= х^-^ = ус^ ^ можно произвести математическое описание процесса шума известными моделями авторегрессии 1-го или 2-го порядков [13].
Для определения порядка модели авторегрессии и выполнения второго этапа оценивания параметров вычислялись асимптотически несмещенные оценки корреляционной функции [14] остаточных ошибок по формуле:
1 N-m
КЛ т ] = IV £ £кек+т • (4)
где временной сдвиг т =0, 1, 2. М < N.
Построенный график корреляционной функции (рис. 2) свидетельствует о том, что описание процесса шума для модели следует производить моделью авторегрессии 2-го порядка. Шумовая компонента модели авторегрессии 2-го порядка [14] имеет вид:
юк = Ю-\ +Ю-2 + *к, (5)
где - оцениваемые параметры.
т
т
Рис. 2. График корреляционной функции
Оценка этих параметров производилась путём минимизации функционала:
М 2
^ [й,ц) = X (К£[т\- 1лК8\_т-1]-цК^[т-2]) ,т.е. ¡,ц находились из ре-
2 т=1 ]
шения соответствующей системы нормальных уравнений, а оценка коэффициента у осуществлялась по формуле [8]:
у,
1
Кс[ 0]
К([0]
(6)
где Кю[°] = -
1 -Ц
2 '
(1+ Ц) (1 -Ц)
Для оценивания параметров АР2 были получены следующие системы нормальных уравнений, по которым были вычислены оценки параметров ¡й;Ц;у.
1260,9663¡й+ 514,8224ц = 289,5179;
514,8224й +1307,1876ц = 49,6265. ^
й = ■+0,2551; Ц = -0,0625; у = 5,5246 Таким образом, были получены все оценки параметров модели, которые представляются в следующем виде:
х, = 0,2196х, ..+1,5018^ -701,141 + 5,5321®,; к к-1 к к
( = 0,2551® 0,0625® ~ к к-1 к-2
+ 1
(8)
к
Условные математические ожидания в соответствии с (7 и 8) имеют вид:
х, = 0,2196хЛГ . 1/лг+1,50Шлг, .-701,141+5,532(лг . .т; к N+1-И N ' N+1 ' N+1-И N
В^т 1/л7= 0,2551ВЛ7 . 1/Л7--0,0625ВЛГ . л/лг N+1/^ ' N+1-1/N ' N+1-2/N
При х^цу- = = 134,44мм;
, л ос с 1 +3,801мм 347мм
^ / N + 'Нр = 0'2551 5,5321 - 0'0625 5,5321 = 0,1375мм
Результаты прогнозирования методом инверсной верификации представлены в табл. 2.
Таблица 2
Результаты прогнозирования процесса перемещения гребня плотины
Х перемещ.2005 мм Дата прогноза, месяц Прогнозная модель
Хпрогн ., мм Ак 2, мм
1 2 5 6
127,04 Янв.2005 129,14 -2,10
120,26 февраль 113,16 +7,10
100,51 март 93,14 +7,37
75,64 апрель 63,21 +12,43
66,15 май 80,57 -14,42
91,80 июнь 101,77 -9,97
104,62 июль 113,40 -8,78
116,67 август 121,52 -4,85
130,00 сентябрь 134,03 -4,03
135,64 октябрь 140,98 -5,34
138,21 ноябрь 140,91 -2,70
134,36 Дек.2005 129,03 +5,33
Выводы.
1. Построенная прогнозная модель позволяют отслеживать основные закономерности развития процесса деформации гребня плотины.
2. Большинство ошибок прогноза перемещений Ак в сравнении с самими перемещениями за исключением отдельных месяцев лежат в пределах графической точности полученных исходных данных по обеим построенным моделям (неточность исходных данных с использованием рис. 1 вносит неточность в результаты прогнозирования).
3. Большие по величине ошибки прогнозирования приходятся на время конца срабатывания и начала нагружения плотины (апрель-июнь); очевидно, что инъецирование трещин [7] изменило характер работы плотины, что препятствует возврату плотины в исходное состояние. Представляется, что для данного периода работы плотины необходимо строить отдельные прогнозные модели.
4. Повысить точность прогнозирования возможно за счет увеличения периода основания прогноза при построении модели (включение нескольких лет натурных наблюдений в период основания прогноза), что планируется показать в последующем при построении прогнозной модели 2007-2009 годов.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Методика определения критериев безопасности гидротехнических сооружений. РД 153-34.2-21.342-00. - М. 2001. - 80 с.
2. Ф1.22. Федеральный закон «О безопасности гидротехнических сооружений / Собрание законодательства Российской Федерации. - 1997. - № 30. - Ст. 3589.
3. Дурчева В. Н., Пучкова С. М., Загрядский И. И. Учет сезонных изменений схемы работы бетонных плотин при анализе данных натурных измерений // Известия ВНИИГ им.Б. Е. Веденеева. - 2000. - № 237. - С. 45-53.
4. Гордон Л. А., Соколовский И. К., Цовикян Л. Х. Прогноз перемещений арочной плотины на основе идентифицированной прогнозной модели // Известия ВНИИГ им. Б. Е. Веденеева. - 1989. - Т. 214. - С. 124-132.
5. Александров Ю. Н. Использование расчетной модели плотины Саяно-Шушенской ГЭС для оценки и прогнозирования ее состояния // Гидротехническое строительство. -2008. - № 11. - С. 64-69.
6. Параметрическая идентификация расчетных моделей гидротехнических сооружений / Д. А. Ивашинцов, А. С. Соколов, С. Г. Шульман, А. М. Юделевич. - СПб: Изд-во ОАО «ВНИИГ им. Б. Е. Веденеева». 2001. - 432 с.
7. Вульфович Н. А., Гордон Л. А., Стефаненко Н. И. Арочно-гравитационная плотина Саяно-Шушенской ГЭС (Оценка технического состояния по данным натурных наблюдений). - СПб.: Изд-во ОАО «ВНИИГ им. Б.Е. Веденеева», 2012. - 107 с.
8. Гуляев Ю. П. Прогнозирование деформации сооружений на основе результатов геодезических наблюдений. - Новосибирск: СГГА, 2008. - 256 с.
9. Гуляев Ю. П., Хорошилов В. С., Лисицкий Д. В. О корректном подходе к математическому моделированию деформационных процессов инженерных сооружений по геодезическим данным // Изв. вузов. Геодезия и аэрофотосъемка. - 2014. - № 4С. - С. 22-30.
10. Хорошилов В.С., Павловская О.Г. Математическое моделирование оползневых процессов по геодезическим данным в условиях проведения взрывных работ и вывоза грунта // Интерэкспо ГЕО-Сибирь-2014. Х Междунар. науч. конгр. : Междунар. науч. конф. «Геодезия, геоинформатика, картография, маркшейдерия» : сб. материалов в 2 т. (Новосибирск, 8-18 апреля 2014 г.). - Новосибирск: СГГА, 2014. Т. 1. - С. 156-161.
11. Гуляев Ю. П., Каленицкий А. И. Структурная идентификация прогнозных моделей локальной техногенной геодинамики и деформаций // Совр. проблемы геодезии и оптики: Тез. докл. - Новосибирск, СГГА, 1998. - С. 121.
12. Павловская О. Г., Хорошилов В. С. Методика выделения однородных оползневых зон по результатам геодезических наблюдений // ГЕО-Сибирь-2011. VII Междунар. науч. конгр. : сб. материалов в 6 т. (Новосибирск, 19-29 апреля 2011 г.). - Новосибирск: СГГА, 2011. Т. 1, ч. 1. - С. 84-88.
13. Гуляев Ю. П., Хорошилов В. С. Математическое моделирование. Прогнозирование деформаций сооружений гидроузлов по геодезическим данным. - Новосибирск: СГГА, 2014. - 78 с.
14. Бокс Дж., Дженкинс Г. Анализ временных рядов. Прогноз и управление. - М.: Мир, 1974. Вып.1. - 405 с. Вып.2. - 197 с.
© Н. Н. Кобелева, В. С. Хорошилов, 2015
