2023年第十五届电工杯A题-Python代码（个人观点）

  新能源主导电力体系是解决温控型负荷电器数量众多的重要举措。合理调控温控型

负荷用电方式，可提供电力系统调节能力并降低用电成本。本文研究电采暖负荷参与电

力系统功率调节的技术经济分析，建立了对室外、室内、墙体系统热模型，利用迭代法

求解集总参数的常微分方程，探究了不同条件下电采暖稳态运行的状态参数及影响。

  问题一：通过数值方法求解微分方程稳态解的性质，运用 Euler 方法迭代方法更新

温度，建立了典型住户温变模型，并评估了室外温度对其用电行为的影响。最后调整模

型参数实现对室内温度和室外温度的控制。针对第二小问，调整电采暖设备开关状态和

室外温度来研究影响其运行特性和耗电量的因素。

  问题二：要求分析电采暖负荷作为调度资源参与功率调节的技术潜力。通过计算不同

条件下可持续调节功率(上调/下调)的时间,评价电采暖负荷的调节能力及其与室外温

度的依赖关系。

  问题三：本文研究了多个电采暖负荷整体作为调度资源参与功率调节的潜力。通过计

算不同室外温度下 6 个典型住户的正常用电功率和可调节功率，本文分析温度对总调节

能力的影响，并为后续的技术经济分析提供了基础。本文选定多个代表性初始状态进行

计算，考虑多户情况下温度变化的随机性，评价了多个电采暖负荷作为调度资源的特征。

  问题四：分析了住宅区较大规模电采暖负荷整体参与系统调节的技术特征。需选择代

表性条件，确定各住户初始状态，在不同室外温度下计算小区日内室温变化、设备开关

状态及相应的用电功率，评价电采暖负荷的调节潜力。选取代表性初始状态，评估其作

为系统调度资源的潜力。本文从多个典型住户出发，考虑随机性选取初始状态，分析住

宅区大容量电采暖负荷整体参与系统调节的技术特征。

  问题五：要保证居民基本舒适度，计算最大可持续调节功率，用削峰或填谷的方式参

与系统功率调节。评估方案对室温和舒适度的影响，并根据影响大小进行调整。评估方

案的经济效益，包括住宅区日常电费和提供调节服务的收益。需考虑电价机制和政策环

境，并进行分析比较。

  问题六：分析了大面积省级电采暖负荷和典型城市空调负荷参与调节的潜力，提出了

制定技术标准、政策法规，选择适宜方式和机制等建议。从宏观上阐述了温控型负荷大

规模参与电网调节的展望，指出了潜力和问题所在，并给出了推进的建议和路径选择。

第一问代码：

import numpy as np

import matplotlib.pyplot as plt

from scipy.integrate import solve_ivp

from matplotlib.font_manager import FontProperties

# 参数

Cin = 1.1e6

Cwall = 1.86e8

R1 = 1.2e-3

R2 = 9.2e-3

# 微分方程

def dTdt(t, y, Pheat_func):

qin, qwall = y

Pheat = Pheat_func(t)

dqin_dt = -(qwall - qin) / (R1 * Cin) + Pheat / Cin

dqwall_dt = (qin - qwall) / R1 / Cwall - (qwall - qout) / R2 / Cwall

return [dqin_dt, dqwall_dt]

# 初始条件

qin0 = 18 # 初始室内温度，单位：摄氏度

qwall0 = 18 # 初始墙体温度，单位：摄氏度

qout = 0 # 室外温度，单位：摄氏度

# 时间区间

t_span = (0, 3600) # 以秒为单位的时间区间

t_eval = np.linspace(t_span[0], t_span[1], 1000) # 评估的时间点

# 制热功率函数

def Pheat_func(t):

if t < 3600: # 前 1 小时制热功率为 0

return 0

else:

return 8 # 之后的制热功率为 8000 瓦特

# 求解微分方程

sol = solve_ivp(dTdt, t_span, [qin0, qwall0], args=(Pheat_func,), t_eval=t_eval)

# 使用中文字体

font = FontProperties(fname=r"c:\windows\fonts\simsun.ttc", size=14)

# 绘制结果

plt.plot(sol.t / 3600, sol.y[0], label='室内温度 qin(t)')

plt.plot(sol.t / 3600, sol.y[1], label='墙体温度 qwall(t)')

plt.xlabel('时间 (小时)', fontproperties=font)

plt.ylabel('温度 (°C)', fontproperties=font)

plt.legend(prop=font)

plt.title('温度 vs. 时间', fontproperties=font)

plt.show()

其余代码文字量过大，需要的宝宝们关注并私信我😀