Файл: Резниковский, А. Ш. Управление режимами водохранилищ гидроэлектростанций.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 22.10.2024
Просмотров: 61
Скачиваний: 0
П Р И Л О Ж Е Н И Е 1
ОСНОВЫ МЕТОДА СТОХАСТИЧЕСКОГО УПРАВЛЕНИЯ
В |
основ е |
о п и сы ва е м о го |
м етод а у п р а в л е н и я |
л е ж и т |
векто р но е |
|||||
сто х а с ти че ск о е д и ф ф е р е н ц и а л ьн о е |
у р а в н е н и е в о д н о го б а л а н с а , |
з а п и |
||||||||
са н н о е |
в |
ф орм е, п р и в е д е н н о й |
в р аб о те [ Л . 3 7 ] : |
|
|
|
|
|||
|
|
|
|
V = A ,Q e + B ,Q H.6 + С г7). |
|
|
|
( 1 |
||
З д е с ь |
V — ( Г , ; |
. . . ;Vn) — |
векто р |
объем о в н а п о л н е н и я |
вод о хр ани - |
|||||
л и щ ; п — к о л и че с тв о с т у п е н е й в к а с к а д е ; <Зн .б = |
(Ф н .б ь Qн .б п ) — |
|||||||||
ве кто р |
о тд а ч |
Г Э С |
к а с к а д а ; |
Qe — |
(Q s il ••• 1 *3бк)— векто р |
боковой |
п р и - |
|||
то чн о сти |
к во д о хр а н и ли щ ам ; k — чи сл о п р и то ко в ; |
v) = ( yj, ; |
. . . |
; т|4) |
— |
о б о б щ е н н а я с л у ч а й н а я в е к т о р -ф у н к ц и я тр е б о в а н и й р а з л и ч н ы х во д о п о л ь зо в а те л е й и во д о п о тр е б и те л е й ; s — ч и с л о в о д о п о л ь зо в а те л е й и
в о д о л о тр е б и те л е й ; |
A i,Bi, С 2 — п о с т о я н н ы е м а тр и ц ы . |
|
|
|||
С т о х а с т и ч е с к а я |
с у т ь |
в о д о б а ла н со в о го |
у р а в н е н и я |
о б у с л о в л е н а |
||
в е р о я тн о стн о й |
п р и р о д о й |
е сте ств е н н о го |
п р и т о к а |
к |
в о д о х р а н и л и - |
|
|
|
|
|
■+ |
|
|
щ а м Qo и с л у ч а й н ы м х а р а к те р о м п о тр е б л е н и я во д ы т). О ч е в и д н о , чт о
т о м у и д р у го м у н е о бход им о |
д а т ь п о д х о д я щ е е |
а н а л и ти ч е с к о е |
о п и с а |
||||||||||||
ние. С тр о го е |
м а те м а ти че ск о е |
з а д а н и е |
сто к о во го |
пр о ц е сс а |
тр е б уе т |
||||||||||
п р и в л е ч е н и я |
со вр е м е н н о го |
а п п а р а т а |
те о р и и |
с л у ч а й н ы х |
пр оц ессо в. |
||||||||||
Р е ч н о й |
сто к (б о к о в у ю |
п р и т о ч н о с т ь ) |
буд е м |
з а д а в а т ь |
т -м е р н ы м |
||||||||||
м а р к о в ск и м |
э р го д и че ск и м |
га р м о н и зу е м ы м |
|
п р оц ессо м |
|
с н е п р е р ы вн ы м |
|||||||||
врем ене м и |
п е р ио д о м |
в |
о д и н |
год . |
В ы б о р |
н е п р е р ы вн о го |
врем ени |
||||||||
в у р а в н е н и и |
в о д н о го |
б а л а н с а |
о б усл о в л е н |
н а л и ч и е м |
р а з р а б о та н н о го |
||||||||||
а н а л и ти ч е с к о го |
а п п а р а т а , |
п о зв о л я ю щ е го н а х о д и т ь |
в |
это м |
с л у ч а е пр и |
о п р е д е л е н н ы х д о п у щ е н и я х и н те р е с ую щ и е н а с у п р а в л я ю щ и е ф у н к
ц и и . К р о м е |
то го , те о р и я р е гу л и р о в а н и я |
сто к а , о п и р а ю щ а я с я |
н а п р е д |
|||||||
с т а в л е н и е |
о |
речно м |
сто к е |
к а к |
о н е п р е р ы вн о м |
с л у ч а й н о м |
процессе, |
|||
б л ц ж е |
все го |
со о тв е т с т в у е т |
и ф и зи че с к о й |
стор он е |
пр о б ле м ы |
с т о х а с т и |
||||
ческо го |
у п р а в л е н и я |
(п р а в о м е р н о ст ь |
п о л о ж е н и я |
о ст а т и с ти ч е с к о й |
||||||
о д н о р о д н о ст и |
в н у т р и го д о в ы х |
с т о к о в ы х |
в е л и ч и н р а с с м а тр и в а е т с я |
в гл . 3 ) .
В к а ч е с т в е а п п р о к с и м и р у ю щ е й м о д ел и д л я к о м п о н е н т к о р р е л и -
р о в а н н о го ве кто р а с т о к а Qo возьм е м с л е д у ю щ у ю а вто р е гр е сс и о н н у ю с х е м у т-го п о р я д к а :
9 (m ,( Q e ) - = a i( 0 < P ( Q e ) + ••• + а т ( О ф ( т - ,) (< Э о )+ < т .( О г . |
(2) |
159
З д е с ь ф — н е ко то р а я |
га у с с о в с к а я |
тр а н с ф о р м а ц и я |
пр о ц е сса ; |
||
v — га у с с о в с к и й « б е л ы й |
ш у м » с M |
v = 0 , |
C o v [ v ( / j ) ; |
v ( ^6(^2)] —= A ) ; |
|
6( r ) — д е л ь т а -ф у н к ц и я Д и р а к а . |
|
|
|
|
|
П о л о ж и в Ф(<2б)=<2б, |
п о л у ч и м , |
чт о |
у р а в н е н и е (2) |
о п и сы вЬ е т |
о б ы ч н ы й га у с с о в с к и й п р о ц е сс; п о л о ж и в - ф ( ф б) = l n Q e и л и , п о п р е д
л о ж е н и ю Н . |
А . |
К а р т в е л и ш в и л и |
[ Л . |
|
37], |
ф ( ( 2 б ) = а 1 п |
In |
( 1 4 - Q s ) -+- |
||||||||||||||
+ Ып Q b +c V Qo, |
п о л у ч и м , |
чт о |
|
у р а в н е н и е 2)( |
б у д е т |
о п и с ы в а т ь |
||||||||||||||||
в |
перво м |
с л у ч а е |
|
л о гн о р м а л ь н ы й |
пр о ц есс, |
а |
во |
второ м |
— пр о ц есс |
|||||||||||||
с п л о т н о с т ь ю |
р аспр е д е ле н и я |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
Р (1, |
Q«) = ^ = X |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
I |
|
|
[« In In (1 +Q6) + |
MnQ6+ |
c VQbY l. |
|
(3) |
|
||||||||||
|
|
X exp ^ ------------------------- 2---------------------f |
|
|
||||||||||||||||||
гд е a, b, c — |
ф у н к ц и и |
врем ени . |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
т |
И зв е с тн о |
[ Л . |
|
6 9 ], ч т о |
т а к н а зы в а е м ы е |
ко эф ф и ц и е н т ы |
|
сноса |
||||||||||||||
&i (t) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
и |
|
д иф ф узи0и ,( 0 в |
ур а в н е н и и |
2)( |
о п р е д е ля ю тс я |
у с л о в - |
|||||||||||||
i = 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
н ы м м а те м а ти че ск и м |
о ж и д а н и е м |
пр о ц е сс а vVf,p (^) = |
А1 [ ф6)( Q/ф (Q c ) |
= |
|
|||||||||||||||||
= ф |
н] , |
гд е |
ф н — |
н а ч а л ь н о е |
со сто я н и е |
пр о ц е сса , |
и его |
а в т о к о р р е л я |
||||||||||||||
ц и о н н о й |
ф у н к ц и е й |
A'v(tu t2). М о ж н о |
п о к а з а т ь , |
ч т о |
|
и м е ю т |
м есто |
с л е |
||||||||||||||
д у ю щ и е ф о р м у л ы ; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
м м |
( О |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
а , |
(0 |
= |
{ |
"Мф |
(t) |
д л я |
г' = |
1 |
|
<сч н та е м М9 > ° > : |
|
|
|
|
||||||
|
|
|
|
|
|
|
О |
|
д л я1 > |
1 . |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4) |
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В \ т ( О |
|
|
дГт—1 |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
Н апр им ер , |
д л я |
т = |
1 п о л учае м |
|
В |
|
2 |
|
|
|
2 а , Д , , . |
||||||||||
|
а , = |
|
51и = |
|
|
|||||||||||||||||
З д е с ь н е о бход им о |
|
о тм е ти ть , |
ч т о |
ф ( £ ) |
|
в у р а в н е н и и |
(2) ц е н тр и р о в а н о |
|||||||||||||||
свои м б е з усл о в н ы м |
м а те м а ти че ск и м |
|
о ж и д а н и е м |
|
М [ ф ( | ) ] . |
П е р е х о д |
||||||||||||||||
н а |
о б щ и й |
|
с л у ч а й |
|
ф о р м а ль н о о с у щ е с тв л я е тс я за м е н о й ф ( £ ) |
н а |
||||||||||||||||
ф ( | ) — М [ ф ( | ) ] . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
Т р а н с ф о р м а ц и я ж е Н . |
А . |
К а р т в е л и ш в и л и н о р м и р о ва н а и ц е н тр и |
|||||||||||||||||||
р о в а н а , |
к о р р е л я ц и о н н ы е за в и с и м о с ти , |
а т а к ж е п а р а м е тар,ы b и с |
|
|||||||||||||||||||
р е ко м е н д уе тс я н а х о д и т ь т а к , к а к э т о п р е д л о ж е н о в [ Л . 3 7 ] . |
|
|
||||||||||||||||||||
|
Д л я н е к о то р ы х ф у н к ц и й ф с у щ е с т в у е т в о з м о ж н о с т ь а н а л и т и ч е |
|||||||||||||||||||||
ск и о п р е д е л я т ь |
и х б е з усл о в н ы е |
ч и с л о в ы е |
х а р а к т е р и с т и к и : |
М Г ф ], |
||||||||||||||||||
D [ф ] и / С , (см . |
[Л . |
2 1 ] ) . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
Р а с ч е т ы п о к а з ы в а ю т , |
чт о д л я ц е ле й п р и к л а д н о г о х а р а к те р а , |
в и |
|||||||||||||||||||
д и м о , д о с та то чн о |
|
б р а т ьт=3. |
В |
это м |
с л у ч а е |
у р а в н е н и е (2 ) |
п р и |
|||||||||||||||
м ет |
вид : |
|
|
|
|
|
ф (3)( < 2 б ) = П 1 ф ( < 2 б ) + с г г у , |
|
|
|
|
|
(5) |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
гд е а , и 0! о п р е д е л я ю тс я п р ит=3 по (4 ).
160
Это уравнение стандартной заменой переменных можно Свести к системе линейных относительно Zi = (p(i_1>(Q6) стохастических диф ференциальных уравнений:
|
z l |
— z2; |
|
|
|
|
|
|
|
z2 = |
z3; |
|
|
|
|
|
|
|
z3 = |
a 3 l z , - + |
a33v |
|
|
|
||
или в векторной форме |
|
|
|
|
|
|
|
|
|
z -= A2 (t) |
z + |
a (t ) |
v. |
|
|
||
Здесь |
0 |
1 |
0 |
|
|
0 |
0 |
0 |
|
|
|
||||||
а 2 (0 - |
0 |
0 |
1 ; а ( 0 = |
0 0 |
0 |
|||
|
Д31 |
0 |
0 |
|
|
0 |
0 |
°зэ |
|
|
0 s i = f li ; |
еУзз— :CTi. |
|
|
|
( 6)
(6')
Уравнение (6') легко обобщается на векторный случай, т. е. на
случай, |
когда у (Q6) = {у! (Q6, y |
ft(Q6. h)}, |
причем каждая ком |
понента |
<ft может удовлетворять |
уравнению |
(2). Характеристиками |
стохастической зависимости между компонентами этого вектора будут
являться взаимокорреляцконные функции процесса <р(Q6) : |
^ |
(П. t 2)' |
|||||
Под компонентами вектора у (Q6) |
подразумеваются расходы |
воды в |
|||||
разных створах одной или нескольких |
рек. |
Матрицы |
Аг(/) |
и ar(t) |
|||
будут иметь при этом следующую структуру: |
|
|
|
|
|||
0 |
1 |
0 . . |
0 |
|
00 |
|
|
0 |
0 |
1 . . |
0 |
|
00 |
|
|
^31 |
0 |
0 . . |
Q |
|
00 |
|
|
(<>“ |
0 |
0 . - |
0 |
|
9 |
|
|
0 |
|
10 |
|
|
|||
0 |
0 |
0 . . |
0 |
|
01 |
|
|
0 |
0 |
0 . . |
a 3k, |
3ft 2 |
00 |
|
|
|
|
|
|
|
|||
0 . . |
|
|
|
|
|
|
|
0 . . |
|
|
|
|
|
|
|
0 . о 83 |
|
• ° 3 1 |
|
• |
®«, 3ft |
|
|
о |
а |
6 • ®зк,з •
• ■ |
• |
at, 3ft |
• °з&, i |
• |
»3fc »3ft |
Элементы матрицы А2(/) находятся, как и ранее, элементы же симметричной матрицы a ( t ) определяются в соответствии со следую-
11—712 |
161 |
fiuiMii общего видя соотношениями:
3ft
|
|
bqf — B qf — 2 |
(a Z i B f i |
+ a H ^ i i ) • |
(7) |
|
где |
|
i= 1 |
|
|
|
|
|
|
|
|
|
|
|
3ft |
|
|
дч+f - z K |
(ti. ^2) |
||
|
|
Bqf (t) |
||||
bqj |
|
|
' q - b |
|
|
|
|
O t 4 - l d t i - |
|
|
|||
1=1s |
|
|
|
|||
|
|
|
|
|
||
Отсюда и из (4) вытекает, |
что матрица в |
(/) |
является функцией |
|||
авто (взаимно) |
корреляционных |
матриц |
Кф (Кф |
|
) рассматриваемого |
процесса.
Заметим, что последние как непрерывные функции двух незави симых аргументов могут быть получены двумерной аппроксимацией исходных 12-мерных теплицевых матриц коэффициентов корреляции между среднемесячными расходами с последующим переходом к корреляции между функциями ср; и ср3-. Можно также их опреде лять методом наибольшего правдоподобия [Л. 37], а затем произво дить двумерную аппроксимацию соответствующих им матриц.
Элементы вектора математических ожиданий стока и его авто(взаимно) корреляционных матриц аппроксимируются тригоно метрическими полиномами с периодом в один год, что и обеспечи-
вает гармонизуемость процесса Q с, а следовательно, и ф(Qe) - Уравнение (6 ') совместно с уравнением (1) и заданными исход
ными значениями дает исчерпывающее описание той или иной управляемой стохастической водохозяйственной системы.
Первое из этих уравнений можно разрешить относительно Qe, а затем, воспользовавшись идеей статистической линеаризации коэф фициентов сноса и диффузии, привести его к виду:
|
|
|
3k |
|
|
|
|
|
|
|
Qe = А3 (/) Q6 + 2 |
<*** (0 |
*v. |
|
|
(8) |
|
|
|
|
i=i |
|
|
|
|
|
Очевидно, |
что |
уравнения |
(8) и (1) путем соответствующей |
за |
||||
мены обозначений |
могут быть объединены следующим образом: |
|
||||||
|
|
|
|
I |
|
|
|
|
|
Х - = А ( t ) x + |
BQH>6 + |
Y i <*t (0 |
|
|
(9) |
||
где A(f), В и O i ( t ) |
|
|
i= l |
|
|
|
|
|
— некоторые матрицы. |
|
|
|
|
||||
Размерность системы (9) |
равна l = n + 3 k + s , где п — число водо^ |
|||||||
хранилищ, k — число притоков к ним, s — размерность |
функции |
т), |
||||||
которая, как и Q e , |
может быть включена в общее уравнение |
(9). |
|
|||||
Управляющей |
функцией |
(управлением) в |
уравнении (9) |
удобнее |
||||
всего принять |
отдачу QH.e, |
т. е. положить |
-у —^ |
^ |
■У |
|
||
U ( х , t ) — |
QH.6 (•*. О- |
Это управление по определению марковское, так как оно зависит от состояний системы на данный момент времени.
Под состояниями системы будут пониматься величины наполне ния всех водохранилищ каскада и величины боковой прпточности к ним в данный момент времени.
162