автореферат диссертации по энергетике, 05.14.02, диссертация на тему:Анализ динамических свойств электроэнергетической системы Вьетнама при малых возмущениях на основе энергетического подхода

МОСКОВСКИЙ ЭНЕРГЕТИЧЕСКИЙ ИНСТИТУТ (ТЕХНИЧЕСКИ УНИВЕРСИТЕТ)

на правах рукописи

ю таи зиеу хонг

АНАЛИЗ ДИНШИЕСКИХ СВОЙСТВ аЛЕКТРОЗНЕРГЕТИЧЕСКОИ СИСТЕМЫ ВЬЕТНАМА ПРИ Ш ВОЗМУЩЕНИЯХ НД. ОСНОВЕ ЭНЕРГЕТКЧЕС1ЮГ0 ПОДХОД А

■ ■ Специальность 05.14.02 -Электрические станции (электрическая часть), сэтн, электроэнергетические системы и управление ими.

АВТОРЕвЕР1Т

, диссертации на соискание ученой степени кандидата технических наук

Москва - 1994

- г -

Работа выполнена в иосковском Энергетическом Институте.

Научный руководитель : доктор технических наук. ;

профессор Зелэнохат Н.И.

Официальные оппонента : Акздемик дэн. доктор технических наук, профессор Шакарян В.Г.

кандидат технических наук,;, профессор Зотов В.И.

Ведущее предприятие : Всесоюзный Государственный проектный и научно-исследовательский

Институт "Энергосвтьцроект" "

. ^ ■■■ 1 •

Защита состоится " 10 " ишя 1994 г. в 15 час ' 00 мин.

на заседании специализированного совета К.053.16.17 в Московском Энергетическом Института в вуд. Г.201 .. ' .

Отзывы, заверенные печатью,просим высылать по адресу : 105835, ГСП, Москва, Е-250, Красноказарменная ул., д.14 Ученый Совет МЭИ.

С диссертацией можно ознакомиться в библиотеке МЭИ.

Автореферат разослан мая 1994 г.

' Ученый секретарь специализированного совета, кандидат технических наук, доцент

0Я4АЛ ХАРАКТЕРИСТИКА РАБСУШ

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

В нзетояцео время две части ЭЭО, расположенные соответственно на северной и гхаой половинах территории Вьетнама, работают раздельно. СеЕ&рная часть ЭЭС Вьетнама достаточно хорошо исследована. Интерес представляет южнзл часть, являвшаяся достаточно концентрированной системой с развитой электрической сетью. При исследовании стационарных рекимов для нее актуальной задачей является обеспечение затухания взаимных качаний роторов синхронных генераторов, шявленио условий, при которых сохраняется статическая устойчивость системы.

Решении подобного рода задач посвящены многочисленные научно-исследовательские работы, в первую очередь, заполненные на кафедре "Электроэнергетические система" МЭИ.

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

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

ботаншй на кафедре "Элоктроэноргетпчоскиэ системы" МЭИ энергетический подход к исслэдовангао устойчивости сложных ЭЗО.

Целью данной работы является реаение комплекса задач, связанных с применением энергетического подхода к анализу динамических свойств слсияоЯ ЗЭС и оценке эффективности припеве нкя различных законов регулирования возбуждения сишрошшх генераторов для улучшения качества электромеханически пере-' ходаых процессов в ЭЭС.

Основными задачами, решаемыми в диссертационной работа, являются :

1. Анализ особенностей развития ЭЭО Вьетнама я выявление, условий возможного нарушения ее устойчивости при малых возмущениях. " .

2. Разработка математической модели сложной ЭЭО в малых отклонениях и исследование динамических свойств ЭЭО на основе структурного анализа действующих в ней сил. .

3. Исследование динамических свойств сложной ЭЭО на основе внэргэплоского подхода и анализ влияния на дох мероприятий по ослабления действия веконсэрвативных и усиленно действия диссипативных сил.

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

Метода к средства проведения исследований.! При анализе динамических свойств сложной ЭЭО использоваятсь метод малых колебаний и критерии статической устойчивости, в большей мере энергетические критерии. При выполнении исследований использовались численные методы определения критериев устойчивости электроэнергетических систем. - -. Исследования проводились аналитически и численно с применением ЭВМ.

К аввргте представляются I

1.Научное обоснование целесообразности проведения исследований по у лучше шя> динамических свойств ЭЭО Вьэтнамз при малых возмущениях ее роима.

2. Предложенный подход к форштрованию математической модэли сложной ЭЭС в малых отклонениях и структурщй анализ ее динамических свойств с выделением "дейстЕУПзис-в система сил.' .

3. Методика анализа динамических свойств сложной ЭЭО на основе энергетического подхода с оценкой влияния действующих в системе сил.

4. Результаты исследования «¡фэктивности применения структурного метода и энергетического подхода к анализу динамических свойств ЭЭС Вьетнама при магах возмущениях и ьыбору мероприятий по улучшению качества переходных процессов.

Научная новизна.

На основе применения отдельных положений теоретической Кб ханша! и предложенного математического описания ЭЭС в малых отклонениях разработан к реализован при расчетах метод структурного анализа влияния дайствунцих б ЭЭС сил на ее динамические свойства. /.■.:...:

2. Разработан математический аппарат для анализа динами; ческих свойств слоеной ЭЗС с применением энергетического

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

3. Разработано программное обеспечение для выполнения расчетов и определения настроечных параметров АРВ синхронных генераторов о использованием разработанного подхода к анализу динамических свойств сложной ЭЭС.

4. Разработаны рекомендации го анализу динамических ' свойств сложных ЭЗС и оценке влияния па них выбранных^ законов регулирования АРВ синхронных, генерато^хш и использованы при исследовании режимов ЭЭС Вьетнама.'

Практическая ценность. Разработанный математический

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

Разработанные математические модели и программы, позволяющие оценить эффективность выбранной настройки АРВ по ! условию демпфероЕания взаимных качаний роторов, могут бить использованы при проведении исследований и улучшении управления режимами эксплуатируемых электроэнергетических систем.

- б -

Публикации. По материалам диссертационной работы

опубликованы тезисы доклада на Всесоюзном семинаре "Кибернетика елекроэнергетических систем" в Челябинске, 1990 г.

Апробация работы. Материалы диссертационной работы

докладывались и обсуждались на заседании кафедры "Элекроэнер-гетические системы" МЭИ, на Есесогсшом семинаре ."Кибернетика елекроэнергетачэских систем" в Челябинске, 1990 г.

Объем работ. Диссертация состоит из введения, четы- ;

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

содагаша РАБОТЫ Во введения обоснована актуальность темы,* диссертационной работы, в краткой форме изложено содержание проведенных исследований.

В первой.главе дан анализ современного состояния и

шрсшктив развития электроэнергетики Вьетаамз, рассмотрена проблема создания.в вей единой злакррэнаргетической системы.

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

Во второй главе дан структурный анализ динамических

свойств сложной ЭЭС при малых возмущениях ее режима.

Для выявления динамических свойств ЭЭО обычно выполняется известная процедура с использованием статической устойчивости: составляется математическое описание ЭЭО в малых отклонениях, составляется характеристический определитель

т -

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

В этом случае необходимо проводить полнее исследование статической устойчивости ЭЭС, чтобы охарактеризовать ее динамические свойства, что является весьма трудоемкой задачей. Для реиеняя такой задачи на первой стадии исследования предлагается использовать применяемые в теоретической механике понятия дкссилативных, консервативных, неконсерянтивных, гироскопических сил и математические приемы для их определения на основе системы дифференциальных уравнений ЭЭС :

- - р£1 - X + -

• 114.

П п п п

- «Т,с/Тй)[. - $ Р„ * 5 I )' <1)

[=» 1*1 1=11/1

отратвщкх относительное движение ротора генератора I относи- 1 тельно траектории движения ЭЭО как целого (центра инерции) с постоянной инерции :

п

В (1) обозначено: (1>4Э - скорость относительного движения ротора генератора I;

Р„ - г^а!^,,; « В.у^соаа,,; Г?, -'Е^у^аНа^.

. Остальные обозначения соответствуют общепринятым.

В малых отклонениях вта система уравнений в матричной форме имеет вид : ,

ЮС + ВХ + СХ = 0 (2)

где В, В, С - квадратные матрицы порядка гс»л;

X, X, X - столбцовые матрицы переменных : йЛш,

X «

"1 э

бг

йШ

«к

гэ

ййш

пЭ

<П

Лагэ

X = ; х » •

* • . •

причем А0И - А0| - АОд, где

1 П ■'■■:•■■"

«, - — Б Элемента матрк-л О, В, 0 определятся по формулам

А< - : А, - А« - 0 :

л

ъ» - •

№

иг »/с

1=1

- - 2(ТИЛЛ)( 2 ,] - Р?,еовВ„ -

Для анализа, структуры сил, действущих в ЭЭС используется следу падл форма записи системы (2) :

ОХ + Всх + СХ + Ссх + АХ = О

С

гдэ В , 0_ - гаммэтричнае матрицы порядка п-п, а частности

°с "

сс72 ссгг

°сгп

Спя

в, А - кососикмэтрнчные матрицы порядка пчг, в частности

А -

12

12

Чп

2п

-а

О

ск

»«1

иСС|

II

»п 2г»

Элементы матриц 1^,0

II

□с, А определяется по формулам:

Л'

е,

си

«V" Л»«* 1

Действующие в ЭЭС силы в механике принято назвать : ВСХ - дассипативные ; СХ - гироскопически»;

. 0ЙХ - консервативные ; АХ - неконсервативные.

Проведено исследование влияния действу ццих в ЭЭС сил на ее динамические свойства. При втом выявлено, что неконгар-ваташша силы оказывают дестабилизирующее влияние на колеба-тольнул устойчивость системы, способствуя раскачивании роторов ее генераторов. Применением метода Гурвица показ&яо, что их дестабялизирупцеё действие о увеличением диссгагативных сил убывает и вследствие этого улучшаются условия затухания качаний роторов генераторов.

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

АРВ составляющих, как функций параметров режима - в зависимости от того, на усиление каких именно сил в ЭЭС направлено действие АРВ. ■ , ■ л

Показано, что можно ставить задачей осуществление такого управления, при котором компенсировалось Сы действие неконсерватишшх сил, как дестабилизирующих режим ЭЭС, а можно за счет управления создавать или усиливать действие дассипа-' тивннх сил, увеличивая тем самим интенсивность затухания взаимных качаний роторов генераторов.

В третьей главе на основе энергетического подхода,

развиваемого на кафедре "Электроэнергетические система" МЭИ,-разработан матоматическнй аппарат, позеоляпций проводить анализ динамических свойств сложной ЭЭС. Црименовмо энергетического подхода отдано предпочтение после проведанного обзора методов енвлиза динамических свойств ЭЭО, используемых при исследовании статической устойчивости ЭЭО (критерии. Гурвица, Рауса, второй метод Ляпунова и др.).

разработаны математические модели для анализа влияния действупцих в ЭЭС сод на качество переходных процессов и для оценки эффективности применения различного рода мероприятий, улучшающих динамические свойства системы, в частности, учитывающие действие АРВ с законами регулирования, ооеспечивающими компенсацию действия неконсэрзаткишх'. сил-, усиление действия дассипативких сил с учетом выполнения основной функции - поддержание изданного уровня напряжения на : выводах синхронных генераторов в стационарных режимах.

В работе доказано, что последнее уравнение (1) при удовлетворении (п-1) предыдущих автоматически удовлетворяется, т.е. система (1) линейно зависима, поэтому необходимо исключить 'одну переменную и последнее уравнение и получить систему (п-1) независимых уравнений :

П'Х + В'Х + О'Х -О (4)

С учетом этого в соответствии с энергетическим подходом составляются определители Д(О') и Д(АОВ), причэм

Л(40В)

А

где А'г - транспонированная квадратная матрица А'; О - нулевая матрица.

Условие устойчивости сводится к требованию положи--тельностя всех главных диагональных миноров матрицы о'порядка

А® - с',, > 0; Лс2

сг» сгг

> 0; ...; Ас , - Д(С) > 0, (5)

П-1

я главных диагональных миноров матрицы АОВ' порядка (2п-г)*(гп-2) : .

.ЛОВ . п. „АО

А, = Ьс„ > О; Д2

ьси V« ьс2т ьсгг

> 0;

2п-2

ьсп ьс,г ьс/.„-» 0 а,2

ъс1г ъсгг " \z.n-i 'а\г 0

ьс 1.п-1 ьсг.п-1"' ьсп-1.п-1 ,п~1 ~аг.л-1

гг

12 О

Г.п-1

... -а

г.п-1

О О

*1.п-1 *2.п-1 — 0

О"

о

1.П-1 а2,п-(

... О

. □ . а

о

>о

: .' . • (б) В (4), (5) и (6) матрицы В^, А', С* квадратичные порядка (п-1)«). Их элементы определяется гараж?ш:я?я1 :

С11 - Си " ~ ;

' »Гц :

си - си -

Ч?

•т

Чп '

Jтг

-

I т^ т^ -............- *■

аи - - - 7" Кс« - — с1п ] - [сл ~ %)] 5

С 1¿п . ■ ^п

ьси -Чм - т - —)+ Ь«" т" И ■

04 Л»

Критерии устойчивости (5), (6), кок и обычно получа емые на основе второго метода Ляпунова, дают лишь достаточные, но не необходимые условия статической устойчивости.

В работе показано, что затухание качаний роторов генераторов в ЭЭС возможно лшь под действием диссипативных, а не каких либо других силг и сделан вывод, что для обеспечения демпферования качаний роторов действием регулирующих устойств должны создаваться силы, эквивалентные диссипативным. Если компенсировать действие неконсервативаых сил, то в ЭЭО при наличии диссипативных сил оСэспачиваатся затухвнг переходного процесса, а при их отсутствии - устанавливается автоколеба-' тельный процесс.

С учетам этого разработаны математические модели, учитывающие действие АРВ с различными законами регулирования, которые можно объединить и записать в обобщенной форме :

Е*с0 + ^го^+ <7>

где 01Э - относительный угол сдвига ротора генератора I.

В работе принято, что введением жесткой отрицательной обратной связи по току ротора компенсирована электромагнитная инерционность роторной цепи (Е^ - ), как это делается у

асинхроннизированных синхронных машин.

С учетом этого а (?) уравнение (2) принимает вид : ,

Й+В^+С^Д-О (8)

гдэ обозначено :

В^ - В + В;; С^ - 0 + С, + с2

в, - кн"'кш; о, - шг'а ; ог = иг'к3

Элементы матрица II порядка (п«п) определяются по формулам :

п

ш„ - 2Е,Уиз1т1{ + I Е^аШО,, - -

- <2ТЛ/ТЛ> ( В.*1€в1лас- I Т 1

2 Е^в^созо,,] . .

1-1

Здесь матрица. Н~' являбтся обратной матрицей Н порядка (п«л).о элементами (при упрощен в виде иг^€«-и^ () :

- 1 + 1 - ^„соэа,, )

" Фй17исоа<й1Г «V

К? и К" - диагональные матрицы коэффициентов усиления Я®, . соответственно.

Для матрицы 0 порядка (п»п) :

п

Если не вводить-упрощение Огч<а.иг{, то выражения для

: элементов -матриц .Н и принимают более слокяую форму, показанную в диссертации. ■ ' •

Переход к линейно независимой системе уреЕнешй осу-; . пдествляется.пересчётом коэф^щиентсв (8). ;

- H -

ЬШ ' bEU - (Tj/TJn>bEta : V

°zu - cEt 1 - <WTJn>czm ' 1

На основе общей математической модели ЭЭО в катя, отклонениях с учетом регулирования получены модели, отвечающие упрощенным законам регулирования. :

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

Разработан математический аппарат для определения условий создания дкссипатййпнт сил в ЭХ и определения величин коэффициентов усиления АРВ, на основе которого аналитически доказано, что в этом случав в законе регулирования возбуждения необходимо включение составляющей, изменяющейся пропорционально отклонению относительной скорости вращения ро ; >ра генератора.

- Разработана методика определения настроечных . параметров АРВ на основе энергетического подхода для общего случая, когда действием АРВ обеспечивается одновременное . создание диссипатшшых сил и компенсация неконсервтивных сил, действующих в-ЭЭО. . ....

В четвертой главе проведено исследование эффективности применения разработанного математического и программного-обеспечения для анализа динамических свойств ЭЭО Вьетнама. Исследования прозедены применительно к че тире хмашинной схеме замещения (рис.7) с учетом действия АРВ с использованием структурного анализа и энергетического подхода при различных законах регулирования дг.я режимов максимальных, средних и минимальных нагрузок ЭЭО Вьатнамз. ■

ПроЕедв!Шыо с пртубненийм энергетического подхода исследования ствтичослсой устс!^шодстя для этих режимов пока- > 'зали, что в- .режимах яакскязлышх• ТГ^бдаих нагрузок можно

-■»et-

обеспечить сохранение статической устойсиеости ЭЭО посредством—-компенсации неконсервативных сил, но при этом ухудшается апериодическая устойчивость, для ев улучшения необходимо увеличивать коэффициенты усиления по отклонению напряжения генератора (pua.г, 3); в режиме минимальных нагрузок компенса- • ция неконсервативных сил также ведет к предотвращению саморас- _ качивзния. но не обеспечивает сохранение апериодической устойчивости. Для сохранения устойчивости ЭЭО Вьетнама в режиме минимальных нагрузок необходимо действием AFB способствовать созданию диссипативных сил в системе, что улучшает также статическую устойчивость системы в режимах максимальных и средних нагрузок.

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

Обращено внимание на то, что хотя алгоритм выбора настроечных коэффициентов Я® при компенсации неконсервативных сил'прощэ, чем при создании диссипативных сил, что .позволяет сократить затрачиваемое время при расчетах на ЭВМ, однако включение составлящей К^(0{Э-04ЭО) в закон (7) к- может быть рекомендовано, так как практическая реализация регулирования • го углу связана с большими трудностями технического характера.

В итоге обоснована целесообразность регулирования . возбуждения генераторов в сложной ЭЭО по закону : .

Eqe„t- Eci0 + Wo. - иг«> + *>»э г «.э0> ~ <*>

Проведено сопоставление результатов исследования статической устойчивости но энергетическому критерию и критерию Гурвица, что позволило подтвердить достоверность исследований с применением энергетического подхода, и тот факт, что определяемая область устойчивости ш энергетическому критерию является подобластью, расположенной внутри области устойчивости по критерию Гурвица, т.е. если ЭЭС устойчива от энергетическому критерию, то сгна обязательно устойчива да критерию . Гурвица. ■

В результате сопоставления выявлены преимущества а недостатки энергетического подхода^чс^анализу статической : устойчивости и оценка динамических свойствЭЭЦ. а также к вы-

гооо ¿¡со 2000 1500

юса

500

к? / к*

1

а; \\\ к?

кЧ

° 20 АО СО во <00 «о

Рис.г. Область устойчивости по энергетическому критерии при регулировании по закону (7), если К^ = О и К^ измаюпвтся в соответствии с графиком рис.3 (рапа* максимальных нагрузок).

«о

100

50

-го

-•(00

: !: - ■ 1 :. Г ••

1 1

-' !

к5 ! г -——

к5" 4 ———1Г~'—* ' 1 К1,

(ООО «00 (200 НОО 4400 1ГС0

Рис.3. Зависимости коэфЗзс-у.еитов Я^ ((-2+4) от коэффициента - .- -(закон (7) при Я" = 0, ре угол максимальных нагрузок).

бору волгшн настроечных параметров АРВ синхронных генераторов^ ОСНОВНЫЕ РЕЗУЛЬТАТЫ ДИССЕРТАОДОННОИ РАБОТЫ

1. Проведенными исследованиями обоснована и подтверздена расчетами необходимость улучпения динамических свойств электроэнергетической системы Вьетнама при малых возмущениях ее режима в условиях изменения суммарной нагрузки системы. - . •

2. Разработано теоретическое и методическое обеспечение для осуществления структурного анализа действующих в сложной электроэнергетической системе сил и для характеристики'es динамических свойств, эфЕектшость которого подтверждена выполненными расчетами примэнительно к многомашинной схеме замещения ЗЭС Вьетнама.

3. Разработано методическое и программное оба спечете для проведения исследования динамических свойств регулируемой ЭЭС с применением энергетического подхода и на его осноге дана оценка эффективности функционирования различных структур автоматических регуляторов еозОуждэеия с определением их настройки по условию обеспечения интенсивного двмпфврования Езаимяых качаний роторов синхронных генераторов в сложной 330.

4. Сопоставлена энергетический подход и метод Гурвица теоретически и расчетами на ЭВМ по специально разработанным программам; выявлены преимущества энергетического подхода при решении задачи выбора мероприятий по улучшению 'динамических свойств сложной ЭЭС.

Б. Проведенными исследованиями установлено, что для ' обеспечения интенсивного демпферования , взаимных качаний роторов синхронных генераторов в ЭЭС Вьетнама нет необходимости . усложнять законы . регулирования на АРВ генераторов всех электростанций, достаточный это сделать лишь " на генераторах одной - двух электростанций, если настройку АРВ осуществить в соответствии с разработанными рекомендациями.

Основные результаты диссертационной работа'отражены ' в публикации : " . ... .

1. Зеленохат Н.И., Поляков М.А., Хо Тхи Зиву ХЬнг. Анализ устойчивости слоюой электреаоергетической"системы // Все союз- -ный семинар "Кибернетика электроендэтетических систем" : Тез.докл. -Челябинск, 19Э0. -С.18-19. .