发明名称 基于滤波与相位协调的静态励磁系统同步信号检测方法
摘要 本发明公开一种基于滤波与相位协调的静态励磁系统同步信号检测方法。该方法通过高压调理变压器及低压调理变压器,将阳极电压调理到同步电压检测板要求的范围内,通过数字采集单元将电压信号离散,进行带通滤波消除阳极电压的换相缺口,再以超前环节对引入滤波所带来的同步信号相位延迟进行补偿,通过多路相电压的相位协调,实现了多路冗余的同步信号检测,在该基础上进行过零检测,输出可控硅触发电路的同步信号,从而实现了同步信号的数字化检测,避免了采用模拟电路存在的检测精度低问题,避免了采用过零检测法存在的误判风险,保证了同步信号的实时性和准确性。
申请公布号 CN104009614A 申请公布日期 2014.08.27
申请号 CN201410280473.1 申请日期 2014.06.20
申请人 哈尔滨电机厂有限责任公司 发明人 赵梦莹;石雨涛;李显彤;解植莹
分类号 H02M1/08(2006.01)I;H02P9/30(2006.01)I 主分类号 H02M1/08(2006.01)I
代理机构 代理人
主权项 一种基于滤波与相位协调的静态励磁系统同步信号检测方法,其特征是:该方法包括如下步骤:1)高压调理变压器(1)接入励磁系统整流桥(4)的阳极电压,高压调理变压器(1)的输出接于低压调理变压器(2)的输入;低压调理变压器(2)接入同步信号检测板(3);2)对待检测的功率整流装置及其组成环节进行编号:将六个可控硅分别编号为T<sub>1</sub>、T<sub>2</sub>、T<sub>3</sub>、T<sub>4</sub>、T<sub>5</sub>和T<sub>6</sub>,整流装置的三相交流输入信号分别用UA、UB、UC表示;3)选择数字信号处理器:根据功率整流装置阳极信号的工作频率范围、检测同步信号算法的计算时间参数确定处理器的采样频率,应保证f<sub>s</sub>≥5f<sub>u</sub>和f<sub>s</sub>≥1/T<sub>proc</sub>,其中,f<sub>s</sub>为数字信号处理器的采样频率,f<sub>u</sub>为阳极信号的最大工作频率,T<sub>proc</sub>为算法的计算时间;4)计算确定测量装置变压器的型号和参数:将整流装置的交流输入信号UA、UB、UC引出,分别接于变压器的源端,利用变压器的连接方式和变比设定,使其输出幅值为±10V的电压信号U<sub>ab</sub>、U<sub>bc</sub>、U<sub>ca</sub>,变压器输出电压应满足数字采样板的输入范围要求,该要求根据数字信号处理器的技术参数确定;5)将变压器的输出电压信号U<sub>ab</sub>、U<sub>bc</sub>、U<sub>ca</sub>接在数字信号处理器的三个输入端口,通过数据采集端口的数字采样功能,处理得到离散化的电压信号,将其标记为<img file="FDA0000524261450000021.GIF" wi="282" he="79" />并且存在以下关系<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>U</mi><mo>&OverBar;</mo></mover><mi>ab</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>0</mn></mrow><mo>&infin;</mo></munderover><msub><mi>U</mi><mi>ab</mi></msub><mrow><mo>(</mo><mi>n</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><mi>&delta;</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>n</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000524261450000022.GIF" wi="1072" he="134" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>U</mi><mo>&OverBar;</mo></mover><mi>ca</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>0</mn></mrow><mo>&infin;</mo></munderover><msub><mi>U</mi><mi>ca</mi></msub><mrow><mo>(</mo><mi>n</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><mi>&delta;</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>n</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000524261450000023.GIF" wi="1067" he="134" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mover><mi>U</mi><mo>&OverBar;</mo></mover><mi>bc</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>0</mn></mrow><mo>&infin;</mo></munderover><msub><mi>U</mi><mi>bc</mi></msub><mrow><mo>(</mo><mi>n</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><mi>&delta;</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>n</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000524261450000024.GIF" wi="1067" he="134" /></maths>其中,T<sub>s</sub>为信号的采样周期,δ(t‑nT<sub>s</sub>)为采样函数,且满足<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>&delta;</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>n</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>t</mi><mo>=</mo><mi>n</mi><msub><mi>T</mi><mi>s</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>t</mi><mo>&NotEqual;</mo><mi>n</mi><msub><mi>T</mi><mi>s</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000524261450000025.GIF" wi="1044" he="157" /></maths>6)滤除<img file="FDA0000524261450000026.GIF" wi="257" he="84" />三路离散相电压信号中的高次谐波:根据阳极电压信号的工作频率范围选择带通滤波器的通带宽度,滤波器选用压控电压源二阶带通滤波器,其传递函数表达式为<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>Y</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mi>X</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mfrac><mi>sRC</mi><mrow><mn>1</mn><mo>+</mo><mo>[</mo><mn>3</mn><mo>-</mo><mi>X</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>]</mo><mi>sRC</mi><mo>+</mo><msup><mrow><mo>(</mo><mi>sRC</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000524261450000027.GIF" wi="1197" he="140" /></maths>其中,X(s)为三相电压信号的象函数,也即X(s)为滤波后信号的象函数,滤波器的中心频率为<img file="FDA0000524261450000028.GIF" wi="260" he="130" />根据滤波器的中心频率和通带范围,确定相应的滤波器参数R和C,上限频率和下限频率的计算公式为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>f</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub><mo>=</mo><mfrac><msub><mi>f</mi><mn>0</mn></msub><mn>2</mn></mfrac><mo>[</mo><msqrt><msup><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mi>K</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn></msqrt><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mi>K</mi><mo>)</mo></mrow><mo>]</mo></mtd></mtr><mtr><mtd><msub><mi>f</mi><mrow><mi>p</mi><mn>2</mn></mrow></msub><mo>=</mo><mfrac><msub><mi>f</mi><mn>0</mn></msub><mn>2</mn></mfrac><mo>[</mo><msqrt><msup><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mi>K</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn></msqrt><mo>+</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mi>K</mi><mo>)</mo></mrow><mo>]</mo></mtd></mtr><mtr><mtd><mi>K</mi><mo>=</mo><mn>1</mn><mo>+</mo><mfrac><msub><mi>R</mi><mi>f</mi></msub><mi>R</mi></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000524261450000029.GIF" wi="1130" he="419" /></maths>其中,f<sub>p1</sub>为下限频率,f<sub>p2</sub>为上限频率,K为放大倍数,R<sub>f</sub>为反馈电阻值,R为前馈电阻值;7)将步骤6中的S域传递函数,转化为Z域上的传递函数<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>M</mi></munderover><msub><mi>b</mi><mi>k</mi></msub><msup><mi>z</mi><mrow><mo>-</mo><mi>k</mi></mrow></msup></mrow><mrow><mn>1</mn><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>a</mi><mi>k</mi></msub><msup><mi>z</mi><mrow><mo>-</mo><mi>k</mi></mrow></msup></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>Y</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow><mrow><mi>X</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00005242614500000210.GIF" wi="1052" he="277" /></maths>将其表示为差分方程为<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>y</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>N</mi></munderover><msub><mi>a</mi><mi>k</mi></msub><mi>y</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>M</mi></munderover><msub><mi>b</mi><mi>x</mi></msub><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000524261450000031.GIF" wi="1091" he="150" /></maths>其中,x(n)是数字采样板采集得到的三路离散相电压信号序列,y(n)是经过滤波后的三路离散相电压信号序列,a<sub>k</sub>和b<sub>k</sub>为相应的系数,N和M为差分方程的阶次,阶次的确定以及系数的设计需要根据不同的滤波器设计方法得到;8)设计数字滤波器:a.将数字滤波器的技术指标,转换为相应的模拟滤波器技术指标;b.按照转换后的的技术指标,设计得到模拟带通滤波器H(s);c.将H(s)转换为H(z),将经过数字滤波后的信号记为U<sub>ab</sub>,U<sub>ca</sub>,U<sub>bc</sub>;9)对经过滤波后的三路相电压信号U<sub>ab</sub>,U<sub>ca</sub>,U<sub>bc</sub>进行相位补偿:补偿采用数字式超前环节,其相应的模拟传递函数为<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mi>G</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>1</mn><mo>+</mo><mi>a</mi><msub><mi>T</mi><mi>c</mi></msub><mi>s</mi></mrow><mrow><mn>1</mn><mo>+</mo><msub><mi>T</mi><mi>c</mi></msub><mi>s</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000524261450000032.GIF" wi="1011" he="147" /></maths>其中,a和T<sub>c</sub>为待设计的参数,根据相位补偿多少来确定,将传递函数转化为z域的传递函数,得到<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>M</mi></munderover><msub><mi>b</mi><mi>k</mi></msub><msup><mi>z</mi><mrow><mo>-</mo><mi>k</mi></mrow></msup></mrow><mrow><mn>1</mn><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>a</mi><mi>k</mi></msub><msup><mi>z</mi><mrow><mo>-</mo><mi>k</mi></mrow></msup></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>Y</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>X</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000524261450000033.GIF" wi="1156" he="277" /></maths>进一步将其转化为差分方程,可得<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msub><mi>y</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>N</mi></munderover><msub><mi>a</mi><mi>k</mi></msub><msub><mi>y</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>M</mi></munderover><msub><mi>b</mi><mi>x</mi></msub><msub><mi>x</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000524261450000034.GIF" wi="1173" he="141" /></maths>其中,y<sub>c</sub>(n)为超前环节的输出序列,x<sub>c</sub>(n)为超前环节的输入序列,并将补偿后的三路相电压信号记为<img file="FDA0000524261450000041.GIF" wi="491" he="91" />和<img file="FDA0000524261450000047.GIF" wi="260" he="104" />10)将补偿后的两路相电压信号进行相位协调,实现冗余的信号检测,具体公式如下:<img file="FDA0000524261450000042.GIF" wi="1769" he="562" />其中,y<sub>o</sub>(n)是相位协调后的输出信号序列;N[·]为故障判断函数,当其输入的信号存在故障时,其输出为0,否则输出为1,T<sub>d1</sub>和T<sub>d2</sub>分别是信号<img file="FDA0000524261450000043.GIF" wi="220" he="85" />和<img file="FDA0000524261450000044.GIF" wi="208" he="85" />相对于<img file="FDA0000524261450000045.GIF" wi="214" he="85" />的超前相位,T<sub>s</sub>为数字信号处理器的采样时间,T<sub>d1</sub>和T<sub>d2</sub>应为T<sub>s</sub>的整数倍;11)对协调后的相电压输出信号序列y<sub>o</sub>(n)进行过零检测,得到同步信号y<sub>out</sub>(n):过零检测环节的表达式为<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><msub><mi>y</mi><mi>out</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>y</mi><mi>max</mi></msub></mtd><mtd><mo>|</mo><msub><mi>y</mi><mi>o</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi></mi><mo>|</mo><mo>&le;</mo><msub><mi>y</mi><mrow><mi>o</mi><mi>min</mi></mrow></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mo>|</mo><msub><mi>y</mi><mi>o</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi></mi><mo>|</mo><mo>></mo><msub><mi>y</mi><mrow><mi>o</mi><mi>min</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000524261450000046.GIF" wi="1143" he="166" /></maths>其中,y<sub>max</sub>表示信号过零时的输出,y<sub>omin</sub>为设定的一个参数,当y<sub>o</sub>(n)的绝对值小于它时,便认为y<sub>o</sub>(n)已小于0,可根据情况设定,通常建议设置为0;12)控制可控硅整流装置的门极电压:在检测得到的同步信号y<sub>out</sub>(n)基础上,叠加由励磁系统自动电压调节器所计算得到的控制角,作为可控硅整流装置的触发脉冲信号,用以控制可控硅的门极电压。
地址 150040 黑龙江省哈尔滨市香坊区三大动力路99号哈尔滨电机厂有限责任公司技术管理部