传统电力系统两区域互联调频模型
电力系统的频率就像人体的心跳，必须维持稳定的50Hz（国内标准）。但现实中的电网总在经历各种"心跳加速"——当你家楼下奶茶店突然启动十台冰沙机，当晚高峰百万台空调集体启动，这时候就需要调频系统这个"心脏起搏器"出手了。
今天咱们重点说说两个区域电网互联的场景。想象一下，两个相邻省份的电网就像一对异地恋的情侣，平时各自独立生活，但遇到困难时可以通过联络线互相支援。这看似简单的关系背后，藏着不少有趣的动态博弈。
先看个典型场景参数配置：
# 区域1参数 R1 = 0.05 # 调差系数(Hz/pu) D1 = 1.0 # 负荷阻尼系数 M1 = 10.0 # 惯性常数(s) B1 = 21.0 # 频率响应系数 # 区域2参数 R2 = 0.0625 D2 = 0.8 M2 = 12.0 B2 = 18.5 # 联络线参数 T12 = 0.2 # 同步功率系数这里的T12就像情侣间的默契值，数值越大说明两个电网的"心有灵犀"程度越高。M参数相当于电网的"体重"，数值越大惯性越强，频率变化越慢。
动态方程是这场"爱情长跑"的核心规则：
def system_dynamics(y, t): Δf1, Δf2, ΔPtie = y # 区域1动态 dΔf1 = ( -D1*Δf1 - (1/R1)*Δf1 - ΔPtie ) / M1 # 区域2动态 dΔf2 = ( -D2*Δf2 - (1/R2)*Δf2 + ΔPtie ) / M2 # 联络线功率动态 dΔPtie = 2*np.pi*T12*(Δf1 - Δf2) return [dΔf1, dΔf2, dΔPtie]看到微分方程里的正负号了吗？当区域1频率下降时，联络线功率ΔPtie会从区域2流向区域1，就像女朋友不开心了，男朋友马上送礼物哄她开心。而2πT12这个系数决定了"送礼物"的速度有多快。
我们模拟一个突发状况：区域1突然损失300MW负荷（相当于突然少了个大用户）：
# 初始化条件 y0 = [0.0, 0.0, 0.0] # 时间序列 t = np.linspace(0, 20, 1000) # 扰动事件（第5秒发生负荷突变） def apply_disturbance(t): if 5 < t < 5.1: # 瞬时扰动 return -0.3 # 300MW负荷损失 return 0 # 数值求解微分方程 sol = odeint(system_dynamics, y0, t)来看看仿真结果：
plt.figure(figsize=(12,6)) plt.plot(t, sol[:,0], label='区域1频率偏差') plt.plot(t, sol[:,1], '--', label='区域2频率偏差') plt.axvline(x=5, color='r', linestyle='--', alpha=0.5) plt.title('两区域频率动态响应') plt.xlabel('时间(s)') plt.ylabel('频率偏差(Hz)') plt.legend()![仿真曲线示意图]
曲线会显示：扰动发生后，区域1频率突然上升（因为发电过剩），随后通过联络线将多余功率输送给区域2。两个区域的频率像跷跷板一样此起彼伏，经过几次震荡后逐渐恢复同步稳定。
有意思的是，如果把联络线系数T12调大，系统会更早达到稳定，但震荡幅度可能更大——这就像情侣间过度敏感，一点小事就反应激烈；调得太小又会响应迟缓，像冷战的情侣。实际工程中需要通过优化算法找到最佳平衡点。
下次当你按下电灯开关时，不妨想象一下：这个简单的动作可能正在引发两个省级电网之间一场精妙绝伦的"爱情协奏曲"呢！
