阻抗扫描（一）——传统频率耦合导纳测量方法

cho_ey

已于 2023-10-17 11:26:45 修改

阅读量655

点赞数

文章标签：经验分享

于 2023-10-17 11:25:10 首次发布

本文链接：https://blog.csdn.net/weixin_45879109/article/details/133876968

版权

前言

介绍了传统频率耦合导纳测量方法。

一、频率耦合导纳模型

新能源并网发电设备的频率耦合导纳模型定义如下：

$-\begin{pmatrix} \mathbf{I_p}\\ \mathbf{I_{p1}} \end{pmatrix}=\begin{pmatrix} Y_{11} & Y_{12}\\Y_{21} & Y_{22} \end{pmatrix}\begin{pmatrix} \mathbf{V_p}\\ \mathbf{V_{p1}} \end{pmatrix}$ （1)

其中， $Y_{11}$ 、 $Y_{12}$ 、 $Y_{21}$ 、 $Y_{22}$ 分别为频率耦合导纳矩阵的四个元素； $\mathbf{I_p}$ 、 $\mathbf{I_{p1}}$ 分别为新能源并网发电设备输出端口处频率为 $f_p$ 和 $f_{p1}$ 的电流分量； $\boldsymbol{\mathbf{V_p}}$ 、 $\boldsymbol{\mathbf{V_{p1}}}$ 分别为新能源并网发电设备输出端口处频率为 $f_p$ 和 $f_{p1}$ 的电压分量； $f_p$ 为正序频率， $f_{p1}$ 为负序频率， $f_{p1}=f_p-2f_1$ 。

二、测量流程

测量新能源并网发电设备的频率耦合导纳时，将频率为 $f_p$ 或 $f_{p1}$ 的小扰动电压串联注入至新能源并网发电设备的外部端口处。然后分别提取端口输出电压 $v_a$ 、 $v_b$ 、 $v_c$ 和输出电流 $i_a$ 、 $i_b$ 、 $i_c$ 中扰动频率 $f_p$ 和耦合频率 $f_{p1}$ 的响应或扰动频率 $f_{p1}$ 和耦合频率 $f_p$ 的响应。由于式(1)中只有两组方程，但导纳矩阵中却有四个未知量，因此需要分别进行两次线性独立的扰动注入，得到式(2)。

$-\begin{bmatrix} \mathbf{I^{'}_p} & \mathbf{I^{''}_p}\\ \mathbf{I^{'}_{p1}} & \mathbf{I^{''}_{p1}} \end{bmatrix}=\begin{bmatrix} Y_{11} &Y_{12} \\ Y_{21} & Y_{22} \end{bmatrix}\begin{bmatrix}\mathbf{V^{'}_p} &\mathbf{V^{''}_p}\\ \mathbf{V^{'}_{p1}} & \mathbf{V^{''}_{p1}} \end{bmatrix}$ (2)

式中，上标“'”表示第一次扰动注入后被测对象的响应，上标“''”表示第二次扰动注入后被测对象的响应。

对式(2)进行矩阵求逆运算后，可得被测新能源并网发电设备的导纳矩阵各元素表达式为：

$\left\{\begin{matrix} Y_{11}=\frac{\mathbf{I^{'}_p}\times\mathbf{V^{''}_{p1}}-\mathbf{V^{'}_{p1}} \times \mathbf{I^{''}_{p}}}{\mathbf{V^{''}_{p}}\times\mathbf{V^{'}_{p1}}-\mathbf{V^{'}_{p}}\times\mathbf{V^{''}_{p1}}} & Y_{12}=\frac{\mathbf{V^{'}_p}\times\mathbf{I^{''}_{p}}-\mathbf{I^{'}_{p}} \times \mathbf{V^{''}_{p}}}{\mathbf{V^{''}_{p}}\times\mathbf{V^{'}_{p1}}-\mathbf{V^{'}_{p}}\times\mathbf{V^{''}_{p1}}}\\ Y_{21}=\frac{\mathbf{I^{'}_{p1}}\times\mathbf{V^{''}_{p1}}-\mathbf{I^{''}_{p1}} \times \mathbf{V^{'}_{p1}}}{\mathbf{V^{''}_{p}}\times\mathbf{V^{'}_{p1}}-\mathbf{V^{'}_{p}}\times\mathbf{V^{''}_{p1}}} & Y_{22}=\frac{\mathbf{V^{'}_p}\times\mathbf{I^{''}_{p1}}-\mathbf{I^{'}_{p1}} \times \mathbf{V^{''}_{p}}}{\mathbf{V^{''}_{p}}\times\mathbf{V^{'}_{p1}}-\mathbf{V^{'}_{p}}\times\mathbf{V^{''}_{p1}}} \end{matrix}\right.$ (3)

传统频率耦合导纳测量方法的主要步骤如图1所示。根据式(3)，可以得到频率fp处被测新能源并网发电设备的频率耦合导纳测量结果，然后，通过扫频法可以逐点测量出被测新能源并网发电设备的宽频带频率耦合导纳特性。

图1 传统频率耦合导纳测量方法

在采用该方法测量新能源并网发电设备的频率耦合导纳模型时，可采用A相采集的电压电流分量。设第一次注入频率为 $f_p$ 的正序电压扰动和第二次注入频率为 $f_{p1}$ 的负序电压扰动后，A相并网电压和并网电流的时域表达式为：

$\left\{\begin{matrix} v^{'}_{a}=V_1\cos (2\pi f_1t^{'}_0)+V^{'}_p\cos (2\pi f_pt^{'}_0+\varphi^{'}_{vp} )+V^{'}_{p1}\cos (2\pi f_{p1}t^{'}_0+\varphi^{'}_{vp1} )\\ i^{'}_{a}=I_1\cos (2\pi f_1t^{'}_0+\varphi_{i1})+I^{'}_p\cos (2\pi f_pt^{'}_0+\varphi^{'}_{ip} )+I^{'}_{p1}\cos (2\pi f_{p1}t^{'}_0+\varphi^{'}_{ip1} ) \\ v^{''}_{a}=V_1\cos (2\pi f_1t^{''}_0)+V^{''}_p\cos (2\pi f_pt^{''}_0+\varphi^{''}_{vp} )+V^{''}_{p1}\cos (2\pi f_{p1}t^{''}_0+\varphi^{''}_{vp1} ) \\ i^{''}_{a}=I_1\cos (2\pi f_1t^{''}_0+\varphi_{i1})+I^{''}_p\cos (2\pi f_pt^{''}_0+\varphi^{''}_{ip} )+I^{''}_{p1}\cos (2\pi f_{p1}t^{''}_0+\varphi^{''}_{ip1} ) \end{matrix}\right.$ (4)

式中，上标“'”和“''”分别表示第一次和第二次注入扰动后的变量；V1、Vp、Vp1与I1、Ip、Ip1分别为电网基波频率f1、频率fp、频率fp1处的电压与电流分量幅值；φvp、φvp1为fp、fp1处的电压分量初相位；φi1、 φip、φip1为f1、fp、fp1处的电流分量初相位。

因此，在一次宽频带频率耦合导纳测量过程中，采样和FFT分析后得到信号的相位是一个与时间t0相关的变量。在某个频率下进行频率耦合导纳测量时，第一次注入频率为fp的正序电压扰动后，信号采样从t'0时刻开始，该情况下采样后的A相并网电压和电流中提取的频率为fp及fp1的分量可表示为相量的形式：

$\left\{\begin{matrix} \mathbf{V^{'}_{p}}=V^{'}_{p}\angle (2\pi f_{p}t^{'}_0+\varphi ^{'}_{vp})\\ \mathbf{V^{'}_{p1}}=V^{'}_{p1}\angle (2\pi f_{p1}t^{'}_0+\varphi ^{'}_{vp1}) \end{matrix}\right.$ (5)

$\left\{\begin{matrix} \mathbf{I^{'}_{p}}=I^{'}_{p}\angle (2\pi f_{p}t^{'}_0+\varphi ^{'}_{ip})\\ \mathbf{I^{'}_{p1}}=I^{'}_{p1}\angle (2\pi f_{p1}t^{'}_0+\varphi ^{'}_{ip1}) \end{matrix}\right.$ (6)

同样，第二次注入频率为fp1的负序电压扰动后，信号采样从t''0时刻开始，采样后的A相并网电压和电流中提取的频率为fp及fp1的分量可表示为：

$\left\{\begin{matrix} \mathbf{V^{''}_{p}}=V^{''}_{p}\angle (2\pi f_{p}t^{''}_0+\varphi ^{''}_{vp})\\ \mathbf{V^{''}_{p1}}=V^{''}_{p1}\angle (2\pi f_{p1}t^{''}_0+\varphi ^{''}_{vp1}) \end{matrix}\right.$ (7)

$\left\{\begin{matrix} \mathbf{I^{''}_{p}}=I^{''}_{p}\angle (2\pi f_{p}t^{''}_0+\varphi ^{''}_{ip})\\ \mathbf{I^{''}_{p1}}=I^{''}_{p1}\angle (2\pi f_{p1}t^{''}_0+\varphi ^{''}_{ip1}) \end{matrix}\right.$ (8)

将式(5)-(8)代入式(3)，可得新能源并网发电设备在fp处的导纳如式(7)所示。因此，在一次宽频带频率耦合导纳测量过程中，采样FFT分析后信号的相位是与第一次扰动注入后采样起始时刻t'0 和第二次扰动注入后采样起始时刻t''0相关的变量。

$\left\{\begin{matrix} Y_{11}=\frac{I^{'}_{p}\angle(2\pi f_{p}t^{'}_0+\varphi^{'}_{ip})\cdot V^{''}_{p1}\angle(2\pi f_{p1}t^{''}_0+\varphi^{''}_{vp1})-V^{'}_{p1}\angle(2\pi f_{p1}t^{'}_0+\varphi^{'}_{vp1})\cdot I^{''}_{p}\angle(2\pi f_{p}t^{''}_0+\varphi^{''}_{ip})}{ V^{''}_{p}\angle(2\pi f_{p}t^{''}_0+\varphi^{''}_{vp})\cdot V^{'}_{p1}\angle(2\pi f_{p1}t^{'}_0+\varphi^{'}_{vp1})-V^{'}_{p}\angle(2\pi f_{p}t^{'}_0+\varphi^{'}_{vp})\cdot V^{''}_{p1}\angle(2\pi f_{p1}t^{''}_0+\varphi^{''}_{vp1})}\\ Y_{12}=\frac{V^{'}_{p}\angle(2\pi f_{p}t^{'}_0+\varphi^{'}_{vp})\cdot I^{''}_{p1}\angle(2\pi f_{p1}t^{''}_0+\varphi^{''}_{ip1})-I^{'}_{p}\angle(2\pi f_{p}t^{'}_0+\varphi^{'}_{ip})\cdot V^{''}_{p}\angle(2\pi f_{p}t^{''}_0+\varphi^{''}_{vp})}{ V^{''}_{p}\angle(2\pi f_{p}t^{''}_0+\varphi^{''}_{vp})\cdot V^{'}_{p1}\angle(2\pi f_{p1}t^{'}_0+\varphi^{'}_{vp1})-V^{'}_{p}\angle(2\pi f_{p}t^{'}_0+\varphi^{'}_{vp})\cdot V^{''}_{p1}\angle(2\pi f_{p1}t^{''}_0+\varphi^{''}_{vp1})} \\ Y_{21}=\frac{I^{'}_{p1}\angle(2\pi f_{p1}t^{'}_0+\varphi^{'}_{ip1})\cdot V^{''}_{p1}\angle(2\pi f_{p1}t^{''}_0+\varphi^{''}_{vp1})-I^{''}_{p1}\angle(2\pi f_{p1}t^{''}_0+\varphi^{''}_{ip1})\cdot V^{'}_{p1}\angle(2\pi f_{p1}t^{'}_0+\varphi^{'}_{vp1})}{ V^{''}_{p}\angle(2\pi f_{p}t^{''}_0+\varphi^{''}_{vp})\cdot V^{'}_{p1}\angle(2\pi f_{p1}t^{'}_0+\varphi^{'}_{vp1})-V^{'}_{p}\angle(2\pi f_{p}t^{'}_0+\varphi^{'}_{vp})\cdot V^{''}_{p1}\angle(2\pi f_{p1}t^{''}_0+\varphi^{''}_{vp1})} \\ Y_{22}=\frac{V^{'}_{p}\angle(2\pi f_{p}t^{'}_0+\varphi^{'}_{vp})\cdot I^{''}_{p1}\angle(2\pi f_{p1}t^{''}_0+\varphi^{''}_{ip1})-I^{'}_{p1}\angle(2\pi f_{p1}t^{'}_0+\varphi^{'}_{ip1})\cdot V^{''}_{p}\angle(2\pi f_{p}t^{''}_0+\varphi^{''}_{vp})}{ V^{''}_{p}\angle(2\pi f_{p}t^{''}_0+\varphi^{''}_{vp})\cdot V^{'}_{p1}\angle(2\pi f_{p1}t^{'}_0+\varphi^{'}_{vp1})-V^{'}_{p}\angle(2\pi f_{p}t^{'}_0+\varphi^{'}_{vp})\cdot V^{''}_{p1}\angle(2\pi f_{p1}t^{''}_0+\varphi^{''}_{vp1})} \end{matrix}\right.$ (9)

分析

采用上述流程，可测量计算获得传统的频率耦合导纳阻抗。对（9）进行化简，可得其表达式为：

$\left\{\begin{matrix} Y_{11}=\frac{I^{'}_{p}\cdot V^{''}_{p1} \angle(-4\pi f_1t^{''}_0+\varphi^{'}_{ip}+\varphi^{''}_{vp1})-V^{'}_{p1}\cdot I^{''}_{p} \angle(-4\pi f_1t^{'}_0+\varphi^{'}_{vp1}+\varphi^{''}_{ip})} {V^{''}_{p}\cdot V^{''}_{p1} \angle(-4\pi f_1t^{'}_0+\varphi^{''}_{vp}+\varphi^{'}_{vp1})-V^{'}_{p}\cdot V^{''}_{p1} \angle(-4\pi f_1t^{''}_0+\varphi^{'}_{vp}+\varphi^{''}_{vp1})}\\ Y_{12}=\frac{V^{'}_{p}\cdot I^{''}_{p} \angle( \varphi^{'}_{vp}+\varphi^{''}_{ip})-I^{'}_{p}\cdot V^{''}_{p} \angle( \varphi^{'}_{ip}+\varphi^{''}_{vp})} {V^{''}_{p}\cdot V^{''}_{p1} \angle(-4\pi f_1t^{'}_0+\varphi^{''}_{vp}+\varphi^{'}_{vp1})-V^{'}_{p}\cdot V^{''}_{p1} \angle(-4\pi f_1t^{''}_0+\varphi^{'}_{vp}+\varphi^{''}_{vp1})}\\ Y_{21}=\frac{[I^{'}_{p1}\cdot V^{''}_{p1} \angle( \varphi^{'}_{ip1}+\varphi^{''}_{vp1})-I^{''}_{p1}\cdot V^{'}_{p1} \angle( \varphi^{''}_{ip1}+\varphi^{'}_{vp1})]\angle[-4\pi f_1(t^{'}_0+t^{''}_0))]} {V^{''}_{p}\cdot V^{''}_{p1} \angle(-4\pi f_1t^{'}_0+\varphi^{''}_{vp}+\varphi^{'}_{vp1})-V^{'}_{p}\cdot V^{''}_{p1} \angle(-4\pi f_1t^{''}_0+\varphi^{'}_{vp}+\varphi^{''}_{vp1})}\\ Y_{22}=\frac{V^{'}_{p}\cdot I^{''}_{p1} \angle(-4\pi f_1t^{''}_0+\varphi^{'}_{vp}+\varphi^{''}_{ip1})-I^{'}_{p1}\cdot V^{''}_{p} \angle(-4\pi f_1t^{'}_0+\varphi^{'}_{ip1}+\varphi^{''}_{vp})} {V^{''}_{p}\cdot V^{''}_{p1} \angle(-4\pi f_1t^{'}_0+\varphi^{''}_{vp}+\varphi^{'}_{vp1})-V^{'}_{p}\cdot V^{''}_{p1} \angle(-4\pi f_1t^{''}_0+\varphi^{'}_{vp}+\varphi^{''}_{vp1})} \end{matrix}\right.$ (10)