|
(21), (22) Заявка: 2006114899/28, 02.05.2006
(24) Дата начала отсчета срока действия патента:
02.05.2006
(43) Дата публикации заявки: 10.11.2007
(46) Опубликовано: 20.08.2008
(56) Список документов, цитированных в отчете о поиске:
SU 1478139 A1, 07.05.1989. SU 834614 A1, 30.05.1981. RU 2282201 C1, 20.08.2006. SU 813283 A1, 15.03.1981. FR 2731077 A1, 30.08.1996.
Адрес для переписки:
634050, г.Томск, пр. Ленина, 30, ТПУ, отдел интеллектуальной и промышленной собственности
|
(72) Автор(ы):
Джумик Дмитрий Валерьевич (RU), Гольдштейн Ефрем Иосифович (RU)
(73) Патентообладатель(и):
Государственное образовательное учреждение высшего профессионального образования Томский политехнический университет (RU)
|
(54) СПОСОБ ОПРЕДЕЛЕНИЯ НАПРЯЖЕНИЯ НЕЙТРАЛИ И ПОЛОЖЕНИЯ НУЛЕВОЙ ТОЧКИ
(57) Реферат:
Способ заключается в измерении массивов отсчетов мгновенных значений фазных напряжений в режиме несимметричной нагрузки фаз и симметричной составляющей сигнала фазы A полученной путем суммирования в дифференциальном усилителе сигнала фазы A uAO1(tj) и суммарного сигнала трех фаз с нейтрали звезды резисторов, полученных в одни и те же моменты времени tj=t1, t2,…, tN, с шагом где Т – период сигнала тока (напряжения); N – число отсчетов на периоде, одновременно определяют несимметричную UAO1 и симметричную UA составляющие сигнала фазы А. Определяют реактивную квазимощность напряжений фазы A перемножают текущие отсчеты сигналов и определяют их активную квазимощность напряжений фазы А Определяют косинус и синус угла F расхождения векторов между напряжением UA и между напряжением UAO1. Определяют напряжение нейтрали UN. Определяют синус угла расхождения векторов D между напряжением UAO1 и напряжением UN; и синус угла расхождения векторов Е между напряжением UA и напряжением UN. Технический результат заключается в повышении точности и простоты вычислений. 3 ил., 2 табл.
Изобретение относится к области систем обработки информации и может быть использовано при определении напряжения нейтрали и положения нулевой точки по массивам отсчетов мгновенных значений фазных напряжений.
Известен способ определения напряжения нейтрали и положения нулевой точки при отсутствии нейтрального провода [Основы теории цепей: Учебник для вузов / Г.В.Зевеке, П.А.Ионкин, А.В.Нетушил, С.В.Страхов. – 5-е изд., перераб. – М.: Энергоатомиздат, 1989. – 528 с.: ил.], заключающийся в том, что измеряют мгновенные значения линейных напряжений, а также проводимости отдельных фаз, а напряжение нейтрали и положение нулевой точки определяют по формуле
Недостатком известного способа является многоэтапность и сложность его реализации, необходимость в дополнительных операциях по измерению линейных напряжений.
Известен способ определения напряжения нейтрали и положения нулевой точки [Основы теории цепей: Учебник для вузов / Г.В.Зевеке, П.А.Ионкин, А.В.Нетушил, С.В.Страхов. – 5-е изд., перераб. – М.: Энергоатомиздат, 1989. – 528 с.: ил.], выбранный в качестве прототипа, заключающийся в том, что измеряют мгновенные значения фазных напряжений, а также проводимости отдельных фаз и нейтрального (нулевого) провода (при наличии такового ), а напряжение нейтрали и положение нулевой точки определяют по формуле
где – мгновенные значения фазных напряжений;
– проводимости отдельных фаз и нейтрального (нулевого) провода.
Недостатком известного способа являются многоэтапность и сложность его реализации.
Задачей изобретения является создание простого, точного, информативного способа определения напряжения нейтрали и положения нулевой точки по массивам отсчетов мгновенных значений фазных напряжений.
Это достигается тем, что способ определения напряжения нейтрали и положения нулевой точки по массивам отсчетов мгновенных значений фазных напряжений включает также как в прототипе измерение мгновенных значений сигналов фазных напряжений. Согласно изобретению по массивам отсчетов мгновенных значений напряжения в режиме несимметричной нагрузки фаз и симметричной составляющей сигнала фазы А полученной путем суммирования в дифференциальном усилителе сигнала фазы A uAO1(tj) и суммарного сигнала трех фаз с нейтрали звезды резисторов, полученных в одни и те же моменты времени tj=t1, t2,…, tN, с шагом
где Т – период сигнала напряжения,
N – число отсчетов на периоде,
одновременно определяют действующие значения напряжений в режиме несимметричной нагрузки фаз несимметричную UAO1 и симметричную UA составляющие сигнала фазы А. Сохраняют каждый цифровой отсчет как текущий и предыдущий, определяют разность и сумму каждой пары текущего и предыдущего значений, перемножают их разность и сумму. Суммируют эти произведения, затем одновременно определяют реактивную квазимощность напряжений фазы А . Далее перемножают текущие отсчеты сигналов и определяют их активную квазимощность напряжений фазы А Затем определяют косинус и синус угла F расхождения векторов между напряжением UA и между напряжением UAO1. Далее определяют действующее значение напряжения нейтрали UN. Затем одновременно определяют синусы углов: расхождения векторов D между напряжением UAO1 и напряжением UN, E между напряжением UA и напряжением UN.
Полученные значения D, Е, F, UN являются исходными данными при определении напряжения нейтрали и положения нулевой точки по массивам отсчетов мгновенных значений фазных напряжений.
Простота предложенного способа заключается в том, что нет необходимости измерения всех фазных напряжений и проводимостей отдельных фаз и нейтрального (нулевого) провода.
Точность предложенного способа заключается в том, что нет необходимости одновременного измерения проводимостей отдельных фаз и нейтрального (нулевого) провода. В настоящее время невозможно обеспечить одновременность измерения всех проводимостей фаз и проводимости нейтрального (нулевого) провода, ввиду чего возникает погрешность.
Предложенный способ является информативным за счет того, что позволяет определять углы расхождения векторов D между напряжением UN и напряжением UAO1, E между напряжением UA и напряжением UN, а также угол расхождения векторов F между напряжением UA и между напряжением UAO1, действующее значение напряжения нейтрали UN.
На фиг.1 приведена структурная схема реализации предложенного способа определения напряжения нейтрали и положения нулевой точки по массивам отсчетов мгновенных значений фазных напряжений.
На фиг.2 изображена аппаратная схема устройства, реализующего рассматриваемый способ определения напряжения нейтрали и положения нулевой точки по массивам отсчетов мгновенных значений фазных напряжений.
На фиг.3 приведена векторная диаграмма линейных и фазных напряжений.
В табл.1 приведены цифровые отсчеты мгновенных значений напряжения в режиме несимметричной нагрузки фаз и симметричной составляющей сигнала фазы А
В табл.2 приведены результаты расчета углов расхождения векторов D между напряжением UN и напряжением UAO1, Е между напряжением UA и напряжением UN, а также угол расхождения векторов F между напряжением UA и между напряжением UAO1, результаты расчета действующего значения напряжения нейтрали UN.
Способ может быть осуществлен с помощью устройства, представленного на фиг.1. В конце ЛЭП установлены регистраторы аварийных ситуаций (РАС, на фиг.1 не показаны) для создания массивов мгновенных значений напряжений конца линии электропередачи с шагом дискретизации t. Устройство для определения напряжения нейтрали и положения нулевой точки состоит из одного блока расчета, входы которого связаны с концом линии электропередачи, а выходы блока расчета подключены к ЭВМ.
Блок расчета D, E, F, UN (фиг.2) содержит два устройства выборки-хранения. К первому устройству выборки-хранения 1 (УВХ 1) последовательно подключены третье устройство выборки-хранения 3 (УВХ 3), инвертор 4, сумматор 5, перемножитель 6, интегратор 7, первый перемножитель-делитель 8. Ко второму устройству выборки-хранения 2 (УВХ 2) последовательно подключены четвертое устройство выборки-хранения 9 (УВХ 4), второй сумматор 10, выход которого подключен ко второму входу первого перемножителя 6. Кроме того, к выходу первого устройства выборки и хранения 1 (УВХ 1) подсоединен второй вход первого сумматора 5, вход преобразователя действующих значений 11 (ПДЗ), а к выходу второго устройства выборки и хранения 2 (УВХ 2) подсоединен второй вход второго сумматора 10, второй вход преобразователя действующих значений 11 (ПДЗ). К первому выходу преобразователя действующих значений 11 (ПДЗ) подключен второй перемножитель 12, к которому последовательно подключены третий сумматор 13, второй инвертор 14, третий перемножитель 15, четвертый сумматор 16, второй перемножитель-делитель 17. Выход второго перемножителя-делителя 17 подключен к ЭВМ, а также к входу третьего 18 и четвертого 19 перемножителей-делителей. Кроме того, к первому выходу преобразователя действующих значений 11 (ПДЗ) подсоединены первый и второй входы четвертого перемножителя 20, выход которого связан с пятым сумматором 21. Также к первому выходу преобразователя действующих значений 11 (ПДЗ) подключен второй вход третьего перемножителя-делителя 18. Ко второму выходу преобразователя действующих значений 11 (ПДЗ) подсоединены первый и второй входы пятого перемножителя 22, выход которого связан со вторым входом пятого сумматора 21, подключенным ко второму входу четвертого сумматора 16. Второй выход преобразователя действующих значений 11 (ПДЗ) связан со вторым входом второго перемножителя 12 и со вторым входом четвертого перемножителя-делителя 19. Входы первого 1 (УВХ 1) и второго (УВХ 2) устройств выборки-хранения связаны с шестым перемножителем 23. Выход шестого перемножителя 23 связан со вторым интегратором 24, выход которого соединен с пятым перемножителем-делителем 25, подключенным ко второму входу третьего перемножителя 15. К каждому устройству выборки-хранения подключен тактовый генератор 26 (ТГ). Кроме того, выход второго перемножителя 12 связан со вторым входом третьего сумматора 13, вторыми входами первого 8 и пятого 25 перемножителей-делителей. Выход первого перемножителя-делителя 8 подключен к ЭВМ, к третьим входам третьего 18 и четвертого 19 перемножителей-делителей, подключенных в свою очередь к ЭВМ. Блок расчета D, E, F, UN (фиг.2) также включает в себя резисторы, соединенные в звезду, входы которых подключены к входу устройства. К нейтрали звезды резисторов подключен второй вход дифференциального усилителя 27 (ДУ), выход которого соединен с входом второго устройства выборки-хранения 2 (УВХ 2). Второй вход дифференциального усилителя 27 (ДУ) подключен к входу устройства.
Все устройства выборки-хранения реализованы на микросхемах 1100СК2. В качестве резисторов могут быть использованы высокоомные активные сопротивления (например, 3 МОм и более). Программатор действующих значений 11 (ПДЗ) выполнен на микроконтроллере серии 51 производителя atmel AT89S53. Дифференциальный усилитель 27 (ДУ), инверторы, сумматоры и интеграторы реализованы на операционных усилителях 140УД17А. В качестве перемножителей и перемножителей-делителей может быть использована микросхема 525ПС3. Тактовый генератор 26 (ТГ) может быть реализован на микроконтроллере АТ80С2051.
Исследования были проведены на подстанции «Октябрьская» 110/35/10 кВ ОАО «Томскэнерго», 2-ая секция шин 10 кВ.
На входы блока расчета D, E, F, UN устройства, реализующего рассматриваемый способ определения напряжения нейтрали и положения нулевой точки по массивам отсчетов мгновенных значений фазных напряжений, подавали одновременно следующие сигналы:
1) – режим несимметричной нагрузки фазы А,
2) – режим несимметричной нагрузки фазы В,
3) – режим несимметричной нагрузки фазы С.
На блок расчета D, E, F, UA на вход первого устройства выборки-хранения 1 (УВХ 1) поступал сигнал uAO1(tj), на соответствующие резисторы поступали сигналы:
1) – на резистор фазы А,
2) – на резистор фазы В,
3) – на резистор фазы С,
на первый вход дифференциального усилителя (ДУ) сигнал uAO1(tj), где tj=t1, t2,…, tn – моменты времени,
– число разбиений на периоде Т,
t=1·10-3 с – шаг дискретизации массивов мгновенных значений токов и напряжений в начале и в конце ЛЭП.
Суммарный сигнал трех фаз с нейтрали звезды резисторов поступал на второй вход дифференциального усилителя 27 (ДУ). С помощью дифференциального усилителя 27 (ДУ) суммировали сигнал фазы А и суммарный сигнал трех фаз с нейтрали звезды резисторов. На выходе дифференциального усилителя 27 (ДУ) получали симметричную составляющую сигнала фазы A uA(tj). Значения сигналов подавали на соответствующие входы шестого перемножителя 23 и записывали в блоки выборки-хранения 1 (УВХ 1), 2 (УВХ 2) и хранили там, как текущие. Массивы значений сигналов с ЛЭП представлены в табл.1. С выхода устройства выборки-хранения 1 (УВХ 1) сигнал UAO1(tj) поступал на второй вход сумматора 5, на первый вход программатора действующих значений 11 (ПДЗ) и в устройство выборки-хранения 3 (УВХ 3), в котором становился предыдущим значением. В то же время с выхода устройства выборки-хранения 2 (УВХ 2), значение сигнала uA(tj) поступало на второй вход второго сумматора 10, на второй вход программатора действующих значений 11 (ПДЗ) и в устройство выборки-хранения 9 (УВХ 4), в котором становилось предыдущим значением. На выходах программатора действующих значений 11 (ПДЗ) получали действующие значения сигналов uAO1(tj) и UA(tj)
(табл.2),
(табл.2).
С выхода устройства выборки-хранения 3 (УВХ 3) предыдущее значение сигнала UAO1(tj) поступало в инвертор 4, в котором отрицательное значение предыдущего сигнала UAO1(tj) преобразовывалось в положительное. С выхода инвертора 4 значение сигнала UAO1(tj) поступало на вход сумматора 5. С помощью сумматора 5 определяли разность текущего и предыдущего значений сигнала UAO1(tj). С выхода сумматора 5 разность текущего и предыдущего значений сигнала UAO1(tj) поступала на вход перемножителя 6. В то же время с выхода устройства выборки-хранения 9 (УВХ 4) предыдущее значение сигнала и uA(tj) поступало на первый вход второго сумматора 10. С помощью второго сумматора 10 определяли разность текущего и предыдущего значений сигнала uA(tj). С выхода второго сумматора 10 разность текущего и предыдущего значений сигнала uA(tj) поступала на второй вход перемножителя 6. С помощью перемножителя 6 значения разности и суммы сигналов перемножали и подавали на вход первого интегратора 7. С помощью первого интегратора 7 суммировали произведения разности и суммы сигналов и определяли значение реактивной квазимощности
(табл.2).
С выхода первого интегратора 7 значение реактивной квазимощности поступало на вход первого перемножителя-делителя 8. Одновременно с описанным выше процессом в шестом перемножителе 23 сигналы uAO1(tj) и uA(tj) перемножали и подавали на вход второго интегратора 24. С помощью второго интегратора 24 суммировали произведения сигналов и определяли значение активной квазимощности
(табл.2).
С выхода второго интегратора 24 значение активной квазимощности поступало на вход пятого перемножителя-делителя 25. В то же время с первого выхода программатора действующих значений 11 (ПДЗ) действующее значение сигнала UAO1 поступало на первый вход второго перемножителя 12, на первый и второй входы четвертого перемножителя 20, на второй вход третьего перемножителя-делителя 18. Со второго выхода программатора действующих значений 11 (ПДЗ) действующее значение сигнала UA поступало на второй вход второго перемножителя 12, на первый и второй входы пятого перемножителя 22, на второй вход четвертого перемножителя-делителя 19. Значения сигналов UAO1 и UA перемножали с помощью четвертого 20 и пятого 22 перемножителей и подавали на первый и второй входы пятого сумматора 21, в котором суммировали произведения сигналов UAO1 и UA. С выхода пятого сумматора 21 сумма произведений сигналов UAO1 и UA поступала на второй вход четвертого сумматора 16. Одновременно с описанным выше процессом во втором перемножителе 12 сигналы UAO1 и UA перемножали и подавали на второй вход пятого перемножителя-делителя 25, на второй вход перемножителя-делителя 8, на первый и второй входы третьего сумматора 13. С помощью пятого перемножителя-делителя 25 определяли значение косинуса угла расхождения векторов F между напряжением UA и между напряжением UAO1, которое затем подавали на первый вход третьего перемножителя 15.
(табл.2).
С помощью первого перемножителя-делителя 8 определяли значение синуса угла расхождения векторов F между напряжением UA и между напряжением UAO1 (см. фиг.3), которое затем подавали на первые входы третьего 18 и четвертого 19 перемножителей-делителей.
(табл.2).
В то же время с помощью третьего сумматора 13 определяли сумму произведений сигналов UAO1 и UA. С выхода третьего сумматора 13 сумма произведений сигналов UAO1 и UA поступала во второй инвертор 14, с помощью которого отрицательное значение суммы произведений сигналов UАО1 и UA преобразовывалось в положительное. С выхода второго инвертора 14 значение суммы произведений сигналов UAO1 и UA поступала на второй вход третьего перемножителя 15. С помощью третьего перемножителя 15 значения сигналов перемножали и подавали на первый вход четвертого сумматора 16, в котором суммировали произведения сигналов. С выхода четвертого сумматора сигнал поступал на первый вход второго перемножителя-делителя 17. При этом третий вход второго перемножителя-делителя 17 имеет отрицательную обратную связь со своим выходом. С помощью второго перемножителя-делителя 17 определяли действующее значение напряжения нейтрали
(табл.2).
С выхода второго перемножителя-делителя 17 действующее значение напряжения нейтрали поступало на третьи входы третьего 18 и четвертого 19 перемножителей-делителей. С помощью третьего перемножителя-делителя 18 определяли значение синуса угла расхождения векторов D между напряжением UN и напряжением UAO1 (см. фиг.3)
(табл.2).
С помощью четвертого перемножителя-делителя 19 определяли значение синуса угла расхождения векторов E между напряжением UA и напряжением UN (см. фиг.3)
(табл.2).
Для наглядного представления полученных результатов на фиг.3 изображена векторная диаграмма линейных и фазных напряжений.
Таким образом, получен простой, точный и информативный способ определения напряжения нейтрали и положения нулевой точки по массивам отсчетов мгновенных значений фазных напряжений.
Табл.1 |
Время t, с |
uA(tj), кВ |
uAO1(tj), кВ |
0 |
0 |
2406,265 |
0,001 |
2523,126 |
4251,135 |
0,002 |
4799,25 |
5679,867 |
0,003 |
6605,596 |
6552,611 |
0,004 |
7765,335 |
6783,944 |
0,005 |
8164,966 |
6351,22 |
0,006 |
7765,335 |
5296,795 |
0,007 |
6605,596 |
3723,883 |
0,008 |
4799,25 |
1786,453 |
0,009 |
2523,126 |
-325,862 |
0,01 |
0 |
-2406,26 |
0,011 |
-2523,13 |
-4251,13 |
0,012 |
-4799,25 |
-5679,87 |
0,013 |
-6605,6 |
-6552,61 |
0,014 |
-7765,33 |
-6783,94 |
0,015 |
-8164,97 |
-6351,22 |
0,016 |
-7765,33 |
-5296,8 |
0,017 |
-6605,6 |
-3723,88 |
0,018 |
-4799,25 |
-1786,45 |
0,019 |
-2523,13 |
325,8617 |
0,02 |
0 |
2406,265 |
Табл.2 |
|
|
UA, кВ |
UAO1, кВ |
cos F |
sin F |
UN, кВ |
sin D |
sin E |
37648,55 |
14030,3 |
5,7735 |
4,8025 |
0,935135 |
0,348492 |
2,1057 |
0,95552 |
0,79482 |
Формула изобретения
Способ определения напряжения нейтрали и положения нулевой точки, включающий измерение массивов отсчетов мгновенных значений фазных напряжений, отличающийся тем, что по массивам отсчетов мгновенных значений напряжения в режиме несимметричной нагрузки фаз и симметричной составляющей сигнала фазы А , полученной путем суммирования в дифференциальном усилителе сигнала фазы A uAO1(tj) и суммарного сигнала трех фаз с нейтрали звезды резисторов, полученных в одни и те же моменты времени tj=t1, t2,…, tn, с шагом ,
где Т – период сигнала тока (напряжения);
N – число отсчетов на периоде,
одновременно определяют действующие значения напряжения в режиме несимметричной нагрузки фаз несимметричную UAO1 и симметричную UA составляющие сигнала фазы А, сохраняют каждый цифровой отсчет как текущий и предыдущий, определяют разность и сумму каждой пары текущего и предыдущего значений, перемножают их разность и сумму, суммируют эти произведения, затем одновременно определяют реактивную квазимощность напряжений фазы A , перемножают текущие отсчеты сигналов и определяют их активную квазимощность напряжений фазы А , затем определяют косинус и синус угла F расхождения векторов между напряжением UA и между напряжением UAO1, определяют действующее значение напряжения нейтрали UN, затем для определения положения нулевой точки одновременно определяют синус угла расхождения векторов D между напряжением UAO1 и напряжением UN; и синус угла расхождения векторов Е между напряжением UA и напряжением UN.
РИСУНКИ
MM4A – Досрочное прекращение действия патента СССР или патента Российской Федерации на изобретение из-за неуплаты в установленный срок пошлины за поддержание патента в силе
Дата прекращения действия патента: 03.05.2008
Извещение опубликовано: 20.03.2010 БИ: 08/2010
|
|