СПОСОБ ОПРЕДЕЛЕНИЯ ТЕКУЩИХ ПАРАМЕТРОВ ЭЛЕКТРИЧЕСКОГО РЕЖИМА ЛИНИИ ЭЛЕКТРОПЕРЕДАЧИ ДЛЯ ПОСТРОЕНИЯ ЕЕ Г-ОБРАЗНОЙ АДАПТИВНОЙ МОДЕЛИ Российский патент 2006 года по МПК G01R27/04 G01R27/08 

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

Известен способ определения текущих параметров электрического режима, например, для построения адаптивной модели ЛЭП [Мурзин П.В., Суворов А.А. Алгоритмы формирования параметров электрического режима в адаптивной модели ЛЭП // В 38 Вестник УГТУ-УПИ. Энергосистема: управление, качество, конкуренция: Сборник докладов II Всероссийской научно-технической конференции. Екатеринбург: ГОУ ВПО УГТУ-УПИ, 2004. №12 (42). с.424-428], выбранный в качестве прототипа, заключающийся в том, что в начале и в конце ЛЭП с помощью регистратора аварийных событий (РАС) получают синхронизированные во времени массивы мгновенных значений токов и напряжений, с помощью или оптоволоконного или высокочастотного канала связи передают второй массив с конца ЛЭП в ее начало. Для определения параметров модели ЛЭП необходимы следующие параметры электрического режима, полученные в начале и конце линии: токи и напряжения первой гармоники, и фазовые углы между ними.

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

где a_k - коэффициент при функции косинуса;

b_k - коэффициент при функции синуса;

k - номер гармоники;

y_i - мгновенное значение сигнала для i-го отсчета;

t_i - интервал от начала периода до (i+1)-го отсчета (в единицах отсчетов);

n - количество отсчетов в периоде;

Δt_i - шаг интегрирования.

Амплитуду и фазу k-й гармоники сигнала вычисляют следующим образом:

где M_k - амплитуда k-й гармоники сигнала;

ψ_k - фаза k-й гармоники, рад.

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

Находят амплитуды и фазы первых гармоник напряжения и тока в начале и конце линии, по которым производят построение адаптивной модели ЛЭП.

Известный способ определения текущих параметров электрического режима ЛЭП с помощью дискретного преобразования Фурье (ДПФ) для построения адаптивной модели линии электропередачи приведен на фиг.1, где КС - канал связи.

Недостатками известного способа являются многоэтапность и сложность его реализации из-за необходимости использования линейной или параболической аппроксимации мгновенных значений при контроле границ интервала интегрирования.

Задачей изобретения является создание простого, точного, информативного способа определения текущих параметров электрического режима линии электропередачи для построения ее Г-образной адаптивной модели.

Это достигается тем, что в способе определения текущих параметров электрического режима ЛЭП для построения ее Г-образной адаптивной модели, включающем так же, как в прототипе, измерение мгновенных значений сигналов напряжений и токов в начале и в конце линии и передачу массивов с конца линии в ее начало по каналу связи (КС).

Согласно изобретению, по массивам отсчетов мгновенных значений тока и напряжения в начале , и в конце линии электропередачи , , полученным в одни и те же моменты времени t_j=t₁, t₂, ..., t_N, с шагом ,

где Т - период сигнала тока (напряжения),

N - число отсчетов на периоде,

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

Полученные значения R₁ и X₁, R₂ и Х₂ являются исходными данными при создании Г-образной адаптивной модели ЛЭП.

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

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

На фиг.1 изображена структурная схема реализации способа определения текущих параметров электрического режима ЛЭП для построения ее адаптивной модели, основанного на дискретном преобразовании Фурье (ДПФ).

На фиг.2 изображена структурная схема реализации способа определения текущих параметров электрического режима ЛЭП для построения ее Г-образной адаптивной модели.

На фиг.3 изображена Г-образная схема замещения линии электропередачи.

На фиг.4 изображена аппаратная схема блока устройства, реализующего рассматриваемый способ определения текущих параметров ЛЭП для построения ее Г-образной адаптивной модели (блок расчета R₁ и Х₁).

На фиг.5 изображена аппаратная схема блока устройства, реализующего рассматриваемый способ определения текущих параметров ЛЭП для построения ее Г-образной адаптивной модели (блок расчета R₂ и X₂).

В табл.1 приведены цифровые отсчеты мгновенных значений напряжений и токов в начале и в конце линии , , , для фазы А.

В табл.2 приведены результаты расчета параметров Г-образной адаптивной модели ЛЭП.

Способ может быть осуществлен с помощью устройства, представленного на фиг.2. В начале и в конце ЛЭП установлены регистраторы аварийных ситуаций (РАС, на фиг.2 не показаны) для создания массивов мгновенных значений токов начала и конца линии и напряжений начала и конца линии электропередачи с шагом дискретизации Δt. Регистраторы аварийных ситуаций через каналы связи связаны системой сбора и обработки информации, которая обычно расположена в начале ЛЭП. Устройство для определения текущих параметров линии состоит из двух параллельно включенных блоков расчета, входы которых связаны с началом линии электропередачи и через канал связи (КС) с ее концом, а выходы блоков расчета подключены к ЭВМ: блока расчета R₁, X₁ и блока расчета R₂, Х₂. Массивы отсчетов мгновенных значений напряжений и токов поступают на соответствующие входные шины блоков.

Блок расчета R₁, X₁ (фиг.4) состоит из первого 1 (УВХ 1) и второго 2 (УВХ 2) устройств выборки и хранения, входы которых подключены к входу устройства. К первому устройству выборки-хранения 1 (УВХ 1) последовательно подключены первый инвертор 3, первый сумматор 4. Ко второму устройству выборки-хранения 2 (УВХ 2) последовательно подключены первый сумматор 4, третье устройство выборки-хранения 5 (УВХ 3), четвертое устройство выборки-хранения 6 (УВХ 4), второй инвертор 7, второй сумматор 8, выход которого связан с первым перемножителем 9. К первому перемножителю 9 последовательно подключены первый интегратор 10, первый перемножитель-делитель 11, выход которого подключен к ЭВМ. Вход пятого устройства выборки-хранения 12 (УВХ 5) подключен к входу устройства. Кроме того, к выходу третьего устройства выборки и хранения 5 (УВХ 3) подсоединен второй вход второго сумматора 8. К пятому устройству выборки-хранения 12 (УВХ 5) последовательно подключены шестое устройство выборки-хранения 13 (УВХ 6), третий сумматор 14, выход которого связан с первым перемножителем 9. К каждому устройству выборки-хранения подключен тактовый генератор 15 (ТГ). Входы третьего 5 (УВХ 3) и пятого 12 (УВХ 5) устройств выборки-хранения связаны со вторым перемножителем 16. Выход второго перемножителя 16 связан со вторым интегратором 17, выход которого соединен со вторым перемножителем-делителем 18, подключенным к ЭВМ. К выходу пятого устройства выборки-хранения 12 (УВХ 5) также подключены третий сумматор 14 и преобразователь действующих значений 19 (ПДЗ), выходы которого подключены к третьему перемножителю 20, связанному с входами первого 11 и второго 18 перемножителей-делителей.

Блок расчета R₂, Х₂ (фиг.5) по структуре аппаратной схемы является аналогичным, но выход первого сумматора 4 связан с пятым устройством выборки-хранения 12 (УВХ 5) и со вторым перемножителем 16, а вход третьего 5 (УВХ 3) устройства выборки и хранения подключен к входу устройства.

Первое 1 (УВХ 1), второе (УВХ 2), третье 5 (УВХ 3), четвертое 6 (УВХ 4), пятое 12 (УВХ 5), шестое 13 (УВХ 6) устройства выборки-хранения могут быть реализованы на микросхемах 1100СК2. Первый инвертор 3 и второй инвертор 7 могут быть реализованы на микросхеме 140УД17А. Первый 4, второй 8 и третий 14 сумматоры могут быть реализованы на операционных усилителях 140УД17А. В качестве первого перемножителя 9, второго перемножителя 16, третьего перемножителя 20, первого перемножителя-делителя 11 и второго перемножителя-делителя 18 может быть использована микросхема 525ПСЗ. Первый интегратор 10, второй интегратор 17 могут быть реализованы на операционном усилителе 140УД17А. Тактовый генератор 15 (ТТ) может быть реализован на микроконтроллере АТ80С2051. Программатор действующих значений 19 (ПДЗ) может быть выполнен на микроконтроллере серии 51 производителя atmel AT89S53.

Для исследования была выбрана воздушная линия электропередачи 500 кВ протяженностью 300 км, выполненная проводом АС-300/66.

На входы блока расчета R₁, X₁ и блока расчета R₂, X₂ устройства, реализующего рассматриваемый способ определения текущих параметров ЛЭП для построения ее Г-образной адаптивной модели, подавали одновременно следующие сигналы:

1) , , на входные шины блока расчета R₁, X₁

2) , , на входные шины блока расчета R₂, X₂,

На блок расчета R₁, X₁ на вход первого устройства выборки-хранения 1 (УВХ 1) поступает сигнал u₂(t_j) на вход второго устройства выборки-хранения 2 (УВХ 2) сигнал u₁(t_j), а на вход пятого устройства выборки-хранения 12 (УВХ 5) сигнал i₁(t_j),

где t_j=t₁, t₂, ..., t_N - моменты времени,

- число разбиений на периоде Т,

Δt=1·10^-3 с - шаг дискретизации массивов мгновенных значений токов и напряжений в начале и в конце ЛЭП.

Массивы значений сигналов с ЛЭП представлены в табл.1. Значения сигналов записывали в блоки выборки-хранения 1 (УВХ 1), 2 (УВХ 2) и 12 (УВХ 5) и хранили там, как текущие, затем с выхода устройства выборки-хранения 1 (УВХ 1) сигнал u₂(t_j) поступал на инвертор 3. С помощью инвертора 3 отрицательное значение предыдущего сигнала u₂(t_j) преобразовывалось в положительное. С выхода инвертора 3 значение сигнала u₂(t_j) поступало на вход сумматора 4. В то же время с выхода устройства выборки-хранения 2 (УВХ 2), значение сигнала u₁(t_j) поступало на второй вход сумматора 4. С помощью сумматора 4 определяли разность значений сигналов u₁(t_j)-u₂(t_j). С выхода сумматора 4 разность значений сигналов u₁(t_j)-u₂(t_j) поступала в устройство выборки-хранения 5 (УВХ 3) и на вход перемножителя 16. Одновременно значение сигнала i₁(t_j) поступало в блок выборки-хранения 12 (УВХ 5) и на второй вход перемножителя 16. Значения сигналов, записанные в блоки выборки-хранения 5 (УВХ 3) и 12 (УВХ 5), хранили там как текущие. С выхода устройства выборки-хранения 5 (УВХ 3) сигнал u₁(t_j)-u₂(t_j) поступал на вход сумматора 8 и в устройство выборки-хранения 6 (УВХ 4), в котором становился предыдущим значением, а с выхода устройства выборки-хранения 12 (УВХ 5), значение сигнала i₁(t_j) поступало на первый и второй входы программатора действующих значений 19 (ПДЗ), затем поступало в устройство выборки-хранения 13 (УВХ 6) и становилось предыдущим значением. На выходах программатора действующих значений 19 (ПДЗ) получали дважды действующее значение сигнала

С выходов программатора действующих значений 19 (ПДЗ) действующие значения сигналов I₁ и I₁ поступали на входы перемножителя 20. С помощью третьего перемножителя 20 значения сигналов I₁ и I₁ перемножали и подавали на входы первого 11 и второго 18 перемножителей-делителей. С выхода четвертого устройства выборки-хранения 6 (УВХ 4) предыдущее значение сигнала u₁(t_j)-u₂(t_j) поступало во второй инвертор 7, с помощью которого отрицательное значение предыдущего сигнала u₁(t_j)-u₂(t_j) преобразовывалось в положительное. С выхода второго инвертора 7 значение сигнала u₁(t_j)-u₂(t_j) поступало на вход второго сумматора 8. В то же время с выхода третьего устройства выборки-хранения 5 (УВХ 3) текущее значение сигнала u₁(t_j)-u₂(t_j) поступало на вход второго сумматора 8, с помощью которого определяли разность текущего и предыдущего значений сигнала u₁(t_j)-u₂(t_j). Одновременно с описанным выше процессом с выхода шестого устройства выборки-хранения 13 (УВХ 6) предыдущее значение сигнала i₁(t_j) поступало на вход третьего сумматора 14 и с выхода пятого устройства выборки-хранения 12 (УВХ 5) текущее значение сигнала i₁(t_j) поступало на вход третьего сумматора 14. С помощью третьего сумматора 14 определяли сумму текущего и предыдущего значений сигнала i₁(t_j). С выхода второго сумматора 8 разность текущего и предыдущего значений сигнала u₁(t_j)-u₂(t_j) поступала на вход первого перемножителя 9, а с выхода третьего сумматора 14 сумма текущего и предыдущего значений сигнала i₁(t_j) поступала на вход первого перемножителя 9. С помощью первого перемножителя 9 значения разности и суммы сигналов перемножали и подавали на вход первого интегратора 10. С помощью первого интегратора 10 суммировали произведения разности и суммы сигналов и определяли значение потери реактивной мощности

С выхода первого интегратора 10 значение потери реактивной мощности поступало на вход первого перемножителя-делителя 11. В то же время с помощью второго перемножителя 16 определяли произведения текущих значений сигналов u₁(t_j)-u₂(t_j) и i₁(t_j), которые поступали на вход второго интегратора 17. С помощью второго интегратора 17 определяли потери активной мощности

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

Работа другого блока расчета R₂, Х₂ аналогична работе первого блока расчета R₁, X₁, но на вход первого устройства выборки-хранения 1 (УВХ 1) поступал сигнал i₂(t_j), на вход второго устройства выборки-хранения 2 (УВХ 2) сигнал i₁(t_j), а на вход третьего устройства выборки-хранения 5 (УВХ 3) сигнал u₂(t_j).

Значения сигналов записывали в блоки выборки-хранения 1 (УВХ 1), 2 (УВХ 2) и 5 (УВХ 3) и хранили там как текущие, затем с выхода устройства выборки-хранения 1 (УВХ 1) сигнал i₂(t_j) поступал на первый инвертор 3. С помощью первого инвертора 3 отрицательное значение предыдущего сигнала i₂(t_j) преобразовывалось в положительное. С выхода первого инвертора 3 значение сигнала i₂(t_j) поступало на вход первого сумматора 4. В то же время с выхода второго устройства выборки-хранения 2 (УВХ 2), значение сигнала i₁(t_j) поступало на второй вход первого сумматора 4. С помощью первого сумматора 4 определяли разность значений сигналов i₁(t_j)-i₂(t_j). C выхода первого сумматора 4 разность значений сигналов i₁(t_j)-i₂(t_j) поступала в устройство выборки-хранения 12 (УВХ 5) и на вход второго перемножителя 16. Одновременно значение сигнала u₂(t_j) поступало в третий блок выборки-хранения 5 (УВХ 3) и на второй вход второго перемножителя 16.

Значения сигналов, записанные в блоки выборки-хранения 5 (УВХ 3) и 12 (УВХ 5) хранили там как текущие. С выхода устройства выборки-хранения 5 (УВХ 3) сигнал u₂(t_j) поступал на вход сумматора 8 и в устройство выборки-хранения 6 (УВХ 4), в котором становился предыдущим значением, а с выхода устройства выборки-хранения 12 (УВХ 5), значение сигнала i₁(t_j)-i₂(t_j) поступало на первый и второй входы программатора действующих значений 19 (ПДЗ), затем поступало в устройство выборки-хранения 13 (УВХ 6) и становилось предыдущим значением. На выходах программатора действующих значений 19 (ПДЗ) получали дважды действующее значение сигнала (табл.2).

С выходов программатора действующих значений 19 (ПДЗ) действующие значения сигналов ΔI₁₂ и ΔI₁₂ поступали на входы перемножителя 20. С помощью третьего перемножителя 20 значения сигналов ΔI₁₂ и ΔI₁₂ перемножали и подавали на входы первого 11 и второго 18 перемножителей-делителей. С выхода четвертого устройства выборки-хранения 6 (УВХ 4) предыдущее значение сигнала u₂(t_j) поступало во второй инвертор 7, с помощью которого отрицательное значение предыдущего сигнала u₂(t_j) преобразовывалось в положительное. С выхода второго инвертора 7 значение сигнала u₂(t_j) поступало на вход второго сумматора 8. В то же время с выхода третьего устройства выборки-хранения 5 (УВХ 3) текущее значение сигнала u₂(t_j) поступало на вход второго сумматора 8, с помощью которого определяли разность текущего и предыдущего значений сигнала u₂(t_j). Одновременно с описанным выше процессом, с выхода шестого устройства выборки-хранения 13 (УВХ 6) предыдущее значение сигнала i₁(t_j)-i₂(t_j) поступало на вход третьего сумматора 14 и с выхода пятого устройства выборки-хранения 12 (УВХ 5) текущее значение сигнала i₁(t_j)-i₂(t_j) поступало на вход третьего сумматора 14. С помощью третьего сумматора 14 определяли сумму текущего и предыдущего значений сигнала i₁(t_j)-i₂(t_j). C выхода второго сумматора 8 разность текущего и предыдущего значений сигнала u₂(t_j) поступала на вход первого перемножителя 9, а с выхода третьего сумматора 14 сумма текущего и предыдущего значений сигнала i₁(t_j)-i₂(t_j) поступала на вход первого перемножителя 9. С помощью первого перемножителя 9 значения разности и суммы сигналов перемножали и подавали на вход первого интегратора 10. С помощью первого интегратора 10 суммировали произведения разности и суммы сигналов и определяли значение потери реактивной мощности (табл.2). С выхода первого интегратора 10 значение потери реактивной мощности поступало на вход первого перемножителя-делителя 11. В то же время с помощью второго перемножителя 16 определяли произведения текущих значений сигналов u₂(t_j) и i₁(t_j)-i₂(t_j), которые поступали на вход второго интегратора 17. С помощью второго интегратора 17 определяли потери активной мощности (табл.2). С выхода второго интегратора 17 значение потери активной мощности поступало на вход второго перемножителя-делителя 18. С помощью первого перемножителя-делителя 11 определяли значение поперечного реактивного сопротивления линии электропередачи (табл.2). С помощью второго перемножителя-делителя 18 определяли значение поперечного активного сопротивления линии (табл.2).

Таким образом, получен простой, точный и информативный способ определения текущих параметров электрического режима линии электропередачи для построения ее Г-образной адаптивной модели.

Табл.1Время t, сu_a1, кВi_a1, кАu_a2, кВi_a2, кА0402,5650,272318,224-0,4310,001453,4090,568370,863-0,6580,002458,1950,817385,817-0,840,003418,1180,976366,675-0,9310,004339,7581,044314,635-0,9310,005231,490,999230,892-0,8630,00698,6970,863125,017-0,7270,007-92,1170,6587,776-0,50,008-223,7140,386-111,259-0,2270,009-335,5710,091-219,5270,0230,01-417,519-0,204-304,4660,3180,011-456,4-0,522-361,2920,5680,012-454,606-0,863-388,8080,7490,013-410,94-0,999-376,2460,8850,014-326,598-1,044-326,5980,9310,015-202,778-0,999-246,4440,8630,016-66,995-0,84-138,1760,7270,01773,574-0,613-22,1320,5450,018209,358-0,31896,3050,1360,019325,402-0,023238,668-0,1360,02408,5470,272322,411-0,431

Табл.2ΔР₁, МВтΔQ₁, МВАрI₁, кАR₁, ОмX₁, ОмΔР₂, МВтΔQ₂, МВАрΔI₁₂, кАR₂, ОмХ₂, Ом2,556855,6850,73374,7491103,43331,86208,661,3962170,24107,04

Предложенное изобретение относится к области систем обработки информации и может быть использовано при управлении линией электропередачи. Задачей изобретения является создание простого, точного, информативного способа определения текущих параметров электрического режима линии электропередачи. Способ определения текущих параметров электрического режима линии электропередачи для построения ее Г-образной адаптивной модели заключается в измерении мгновенных значений сигналов напряжений и токов в начале и в конце линии и передаче массивов с конца линии в ее начало по каналу связи. При этом по массивам отсчетов мгновенных значений тока и напряжения определяют разности массивов мгновенных значений токов начала и конца линии электропередачи, а также разности массивов мгновенных значений напряжений начала и конца линии электропередачи; сохраняют массив значений каждой из указанных разностей, а также массив мгновенных значений токов начала линии и напряжений конца линии, как текущий и предыдущий; определяют значения потери реактивной мощности и потери активной мощности в продольной и поперечной ветвях моделируемой линии электропередачи; на основании значений потерь реактивной и активной мощностей в продольной и поперечной ветвях моделируемой линии электропередачи и полученных произведений действующих значений токов в начале линии, а также действующих значений разности токов начала и конца линии электропередачи определяют активное и реактивное сопротивления в продольной и поперечной ветвях моделируемой линии электропередачи. 2 табл., 5 ил.

Способ определения текущих параметров электрического режима линии электропередачи для построения ее Г-образной адаптивной модели, включающий измерение мгновенных значений сигналов напряжений и токов в начале и в конце линии и передачу массивов с конца линии в ее начало по каналу связи, отличающийся тем, что по массивам отсчетов мгновенных значений тока и напряжения в начале , и в конце линии электропередачи , , полученным в одни и те же моменты времени t_j=t₁, t₂, ..., t_N, с шагом ,

где Т - период сигнала тока (напряжения);

N - число отсчетов в периоде,

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