机器学习

电力系统——基于10机39节点的电力系统仿真(Matlab)

本文主要是介绍电力系统——基于10机39节点的电力系统仿真(Matlab),对大家解决编程问题具有一定的参考价值,需要的程序猿们随着小编来一起学习吧!
电力系统——基于10机39节点的电力系统仿真(Matlab)  

1 引言

目前,随着科学技术的发展和电能需求量的日益增长,电力系统规模越来越庞大,超高压远距离输电、大容量发电机组、各种新型控制装置得到了广泛的应用,这对于合理利用能源,充分挖掘现有的输电潜力和保护环境都有重要意义。另一方面,随着国民经济的高速发展,以城市为中心的区域性用电增长越来越快,大电网负荷中心的用电容量越来越大,长距离重负荷输电的情况日益普遍,电力系统在人们的生活和工作中担任重要角色,电力系统的稳定运行直接影响着人们的日常生活。从技术和安全上考虑直接进行电力试验可能性很小,因此迫切要求运用电力仿真来解决这些问题。 电力系统仿真是将电力系统的模型化、数学化来模拟实际的电力系统的运行,可以帮助人们通过计算机手段分析实际电力系统的各种运行情况,从而有效的了解电力系统概况。本文根据电力系统的特点,利用MATLAB的动态仿真软件Simulink搭建了无穷大电源的系统仿真模型,得到了在该系统主供电线路电源端发生三相短路接地故障并由故障器自动跳闸隔离故障的仿真结果,并分析了这一暂态过程。通过仿真结果说明MATLAB电力系统工具箱是分析电力系统的有效工具。 此次仿真实验在MATLAB/Simulink平台上搭建10机39节点系统配电网的仿真模型,并对搭建仿真模型进行潮流计算。阐述其原理及功能,并结合仿真实验求解出现故障时的功角电压稳定的问题,进一步分析系统的暂态稳定性来验证其功能。

2 案例仿真

2.1 负荷参数 2.2 线路、变压器参数 2.3 发电机参数 2.4 励磁参数

3 10机39节点的仿真

3.1 建立Simulink模型 如图1所示,十机39节点系统共有10台发电机,39个节点,共46条支路,其中5条为联络开关支路、1个电源作为平衡节点。从拓扑图中可以看出,0至17节点作为干线,由干线分别从1、2、5节点引出分支。并根据网上找的资料参考进行各元件(发电机、变压器、线路和负荷)的参数设置,以便运行该模型。 图1 10机39节点模型图 图2 子模块图1 图3 子模块图2
3.2 MATLAB程序实现
% clear
% clc
% Definition for the columns
% system base =100 MVA
 
%% Line data Format (line)
% All values are given on the same system base MVA
% 1: From bus  
% 2: To bus   
% 3: Resistance (pu)
% 4: Reactance  (pu)
% 5: Charge     (pu)
% 6: Transformer Tap Amplitute
% 7: base MVA
% 8: Nomonal Voltage (KV) 
 
%   1   2    3        4       5     6    7   8    
line=[...
  0.0035  0.0411  0.6987  0   100 345
 0.001   0.025   0.75    0   100 345
  0.0013  0.0151  0.2572  0   100 345
 0.007   0.0086  0.146   0   100 345
 0       0.0181  0     1.025 100 22
  0.0013  0.0213  0.2214  0   100 345
 0.0011  0.0133  0.2138  0   100 345
  0.0008  0.0128  0.1342  0   100 345
 0.0008  0.0129  0.1382  0   100 345
  0.0008  0.0112  0.1476  0   100 345
  0.0002  0.0026  0.0434  0   100 345
  0.0006  0.0092  0.113   0   100 345
 0.0007  0.0082  0.1389  0   100 345
  0.0004  0.0046  0.078   0   100 345
  0.0023  0.0363  0.3804  0   100 345
 0.001   0.025   1.2     0   100 345
 0.0004  0.0043  0.0729  0   100 345
 0.0004  0.0043  0.0729  0   100 345
 0       0.02    0     1.07  100 22
 0.0016  0.0435  0     1.006 100 345
 0.0016  0.0435  0     1.006 100 345
 0.0009  0.0101  0.1723  0   100 345
 0.0018  0.0217  0.366   0   100 345
 0.0009  0.0094  0.171   0   100 345
 0.0007  0.0089  0.1342  0   100 345
 0.0016  0.0195  0.304   0   100 345
 0.0008  0.0135  0.2548  0   100 345
 0.0003  0.0059  0.068   0   100 345
 0.0007  0.0082  0.1319  0   100 345
 0.0013  0.0173  0.3216  0   100 345
 0.0007  0.0142  0     1.07  100 22
 0.0007  0.0138  0     1.06  100 345
 0.0009  0.018   0     1.009 100 22
 0.0008  0.014   0.2565  0   100 345
 0.0006  0.0096  0.1846  0   100 345
 0       0.0143  0     1.025 100 22
 0.0022  0.035   0.361   0   100 345
 0.0005  0.0272  0       1   100 22
 0.0032  0.0323  0.513   0   100 345
 0.0006  0.0232  0     1.025 100 22
 0.0014  0.0147  0.2396  0   100 345
 0.0043  0.0474  0.7802  0   100 345
 0.0057  0.0625  1.029   0   100 345
 0.0014  0.0151  0.249   0   100 345
 0.0008  0.0156  0     1.025 100 22
  0       0.025   0       1   100 22];
 
%% Machine Data Format (mac_con)
% 1.  Machine Number
% 2.  Bus Number
% 3.  Base MVA
% 4.  Leakage Reactance x_l(pu)
% 5.  Resistance r_a(pu)
% 6.  d-axis sychronous reactance x_d(pu)
% 7.  d-axis transient reactance x_d(pu)
% 8.  d-axis subtransient reactance x"_d(pu)
% 9.  d-axis open-circuit time constant T_do(sec),
% 10. d-axis open-circuit subtransient time constant T"_do(sec)
% 11. q-axis sychronous reactance x_q(pu)
% 12. q-axis transient reactance x_q(pu)
% 13. q-axis subtransient reactance x"_q(pu)
% 14. q-axis open-circuit time constant T_qo(sec)
% 15. q-axis open circuit subtransient time constant % T"_qo(sec)
% 16. inertia constant H(sec)
% 17. damping coefficient d_o(pu)
% 18. dampling coefficient d_1(pu)
% 19. bus number
% Note: all the following machines use sub-transient model
% 1  2     3      4      5      6     7    8      9    10    11    12    13    14    15    16    17    18  19
mac_con=[
 1000.0 0.030 0.0010  0.200 0.060 0.01  7.000 0.003 0.190 0.080 0.03  1.500 0.005 50.00 0.000 0.00 39 ;
 1000.0 0.350 0.0270  2.950 0.697 0.01  6.560 0.003 2.820 1.7   0.03  1.500 0.005 3.030 0.000 0.00 31 ;
 1000.0 0.304 .00386  2.495 0.531 0.01  5.700 0.003 2.370 0.876 0.03  1.500 0.005 3.580 0.000 0.00 32 ;
 1000.0 0.295 .00222  2.620 0.436 0.01  5.690 0.003 2.580 1.66  0.03  1.500 0.005 2.860 0.000 0.00 33 ;
 1000.0 0.540 0.0014  6.700 1.320 0.01  5.400 0.003 6.200 1.66  0.03  0.440 0.005 2.600 0.000 0.00 34 ;
 1000.0 0.224 0.0615  2.540 0.500 0.01  7.300 0.003 2.410 0.814 0.03  0.400 0.005 3.480 0.000 0.00 35 ;
 1000.0 0.322 .00268  2.950 0.490 0.01  5.660 0.003 2.920 1.86  0.03  1.500 0.005 2.640 0.000 0.00 36 ;
 1000.0 0.280 .00686  2.900 0.570 0.01  6.700 0.003 2.800 0.911 0.03  0.410 0.005 2.430 0.000 0.00 37 ;
 1000.0 0.298 0.0030  2.106 0.570 0.01  4.790 0.003 2.050 0.587 0.03  1.960 0.005 3.450 0.000 0.00 38 ;
 1000.0 0.125 0.0014  1.000 0.310 0.01  10.20 0.003 0.690 0.08  0.03  1.500 0.005 4.200 0.000 0.00 30 ];
 
mac_con(:,5)=0;
 
p0=[1000 520.81 650 632 508 650 560 540 830 250]./1000;%Active Power Generation of PV units
Pn(1:10,1)=mac_con(:,3)/1000; % Nominal Power
 
%% Power System Stabilizer Format (MB)
% Applied power system stabilizer is MBPSS with simplified settings
% Note: All machines use MBPSS with same configuration 
% 1: Global gain (G)
% 2: Frequency of low frequency band (FL) Hz
% 3: Gain of low frequency band (KL)
% 4: Frequency of intermediate frequency band (FI) Hz
% 5: Gain of intermediate frequency band (KI)
% 6: Frequency of high frequency band (FH) Hz
% 7: Gain of high frequency band (KH)
%   1    2  3     4   5      6  7
MB=[1   0.2 30   1.25 40    12 160];
 
%% Excitation System format (AVR_Data)
% All machines use IEEE type 1 synchronous machine voltage regulator combined to an exciter
% 1. Low pass filter time constant (Tr) sec
% 2. Regulator gain (Ka)
% 3. regulator time constant (Ta) sec
% 4. Lead-lag compensator time constant (Tb) sec
% 5. Lead-lag compensator time constant (Tc) sec
% 6. Terminal voltage (pu)
% 7. Lower limit for regulator output (Emin)   
% 8. Upper limit for regulator output (Emax)
%     1    2      3      4     5     6       7    8
  AVR_Data=[...
    0.01  200    0.015   10    1    1.03    -5    5
    0.01  200    0.015   10    1    1.03    -5    5
    0.01  200    0.015   10    1    1.03    -5    5
    0.01  200    0.015   10    1    1.03    -5    5
    0.01  200    0.015   10    1    1.03    -5    5
    0.01  200    0.015   10    1    1.03    -5    5
    0.01  200    0.015   10    1    1.03    -5    5
    0.01  200    0.015   10    1    1.03    -5    5
    0.01  200    0.015   10    1    1.03    -5    5
    0.01  200    0.015   10    1    1.03    -5    5];
 
C0=829.7e-9;
L0=3.220e-3;
R0=1;
Ns=120*60/(2); %Nominal speed of synchronous machines
s=10;
% PSSModel=1;%1:No pSS //1:MB 
%% Bus data (Bus)
% 1. Bus number
% 2. Nominal phase-phase voltage 
%    1    2
Bus = [ ...

  345

  345

  345

  345

  345

  230

  345

  345

  345

  345

  345

  345

  345

  345

   22

   22

   22

   22

   22
  345];
 
   zbase=(line(:,8).^2)./line(:,7);
   line(:,3)=line(:,3).*zbase;
   line(:,4)=line(:,4).*zbase/(120*pi);
   line(:,5)=line(:,5)./zbase/(120*pi);
     
PSSMODEL=1;
3.3 运行结果 根据所建立的Simulink仿真模型和程序,进行运行,所得结果如图4所示。 图4 Machines示波图 图5 System示波图 图6 1-10节点电压 图7 11-20节点电压 图8 21-30节点电压 图9 31-39节点电压
3.4 结果分析 由以上图中所得结果可以看出, 图3和图4所示波形图由一开始的波动逐渐趋于稳定,表明该系统是一个稳定的系统;图5至图8所示的39个节点的电值,都在0.93-1.1范围内波动,而且只有少数几个节点电压小于0.95,说明节点电压值处于正常电压水平。无故障情况下,该电力系统能够保持稳定行,当出现故障时,距离距离故障点越远母线越不容易受到干扰,距离故障点近,母线越不稳定,故障时间越短,各机组之间的功角差值越小,系统就越稳定。
4 总结 用Simulink构建10机39节点系统模型,并了解了相关的操作步骤。将Simulink教学和以上咱们学习的专业内容相结合,Simulink提供了各种仿真工具,尤其是它不断扩展的、内容丰富的模块率,为系统的仿真提供了极大便利。在Simulink平台上拖曳和连接典型模块就可以绘制仿真对象的模型框图,并对模型进行仿真分析。

5 资源下载

这篇关于电力系统——基于10机39节点的电力系统仿真(Matlab)的文章就介绍到这儿,希望我们推荐的文章对大家有所帮助,也希望大家多多支持为之网!