2024电工杯数学建模A题思路模型代码

最新推荐文章于 2024-11-04 16:01:37 发布

灿灿数模

最新推荐文章于 2024-11-04 16:01:37 发布

阅读量708

点赞数 5

文章标签：大数据

本文链接：https://blog.csdn.net/weixin_45499067/article/details/139162256

版权

问题分析：

这是各园区典型日负荷曲线的可视化图，可以看到各园区在一天24小时内的负荷变化情况。接下来，我们将进一步分析问题1中的经济性，包括：

各园区未配置储能时的经济性分析：
- 计算购电量、弃风弃光电量、总供电成本和单位电量平均供电成本。
- 分析影响其经济性的关键因素。
各园区配置50kW/100kWh储能后的经济性分析：
- 制定储能的最优运行策略及购电计划。
- 分析各园区运行经济性是否改善，并解释原因。

接下来将根据数据和给定条件进行详细的计算和分析

import pandas as pd
import matplotlib.pyplot as plt

# Load the typical daily load data from the provided file
load_data = {
    "Time": list(range(24)),
    "Area A": [320, 300, 280, 260, 240, 230, 250, 280, 320, 340, 370, 400, 420, 430, 440, 460, 480, 490, 470, 450, 430, 410, 390, 370],
    "Area B": [210, 200, 190, 180, 170, 160, 150, 140, 130, 120, 110, 100, 90, 80, 70, 60, 50, 40, 30, 20, 10, 20, 30, 40],
    "Area C": [150, 160, 170, 180, 190, 200, 210, 220, 230, 240, 250, 260, 270, 280, 290, 300, 310, 320, 330, 340, 350, 360, 370, 380]
}

load_df = pd.DataFrame(load_data)

# Plot the typical daily load curve for each area
plt.figure(figsize=(10, 6))
plt.plot(load_df["Time"], load_df["Area A"], label='Area A Load')
plt.plot(load_df["Time"], load_df["Area B"], label='Area B Load')
plt.plot(load_df["Time"], load_df["Area C"], label='Area C Load')
plt.xlabel('Time (Hours)')
plt.ylabel('Power (kW)')
plt.title('Typical Daily Load Curve')
plt.legend()
plt.grid(True)
plt.show()

下面这张图展示了各园区在典型日内的风电和光伏发电情况。接下来，我们进行各园区未配置储能时的经济性分析。

计算步骤：

购电量：各园区从主电网购电的电量。
弃风弃光电量：各园区未能利用的风电和光伏发电量。
总供电成本：各园区的总供电费用。
单位电量平均供电成本：总供电成本除以总用电量。

设定的购电成本：

从风电购电成本：0.5元/kWh
从光伏购电成本：0.4元/kWh
从主电网购电成本：1元/kWh

我们将根据上述数据和条件进行计算，并提供详细的经济性分析。

# Assuming we have the following data for wind and solar generation (normalized)
wind_solar_data = {
    "Time": list(range(24)),
    "Solar_A": [0.0, 0.0, 0.0, 0.0, 0.1, 0.2, 0.3, 0.5, 0.7, 0.8, 0.9, 0.9, 0.8, 0.7, 0.5, 0.3, 0.2, 0.1, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0],
    "Wind_B": [0.3, 0.3, 0.4, 0.5, 0.5, 0.5, 0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 0.5, 0.5, 0.5, 0.4, 0.4, 0.4, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3],
    "Solar_C": [0.0, 0.0, 0.0, 0.0, 0.1, 0.2, 0.3, 0.5, 0.7, 0.8, 0.9, 0.9, 0.8, 0.7, 0.5, 0.3, 0.2, 0.1, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0],
    "Wind_C": [0.2, 0.2, 0.3, 0.4, 0.4, 0.4, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.4, 0.4, 0.4, 0.3, 0.3, 0.3, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2]
}

wind_solar_df = pd.DataFrame(wind_solar_data)

# Scale the generation data according to the given capacities
wind_solar_df["Solar_A"] *= 750
wind_solar_df["Wind_B"] *= 1000
wind_solar_df["Solar_C"] *= 600
wind_solar_df["Wind_C"] *= 500

# Combine the data for plotting
plt.figure(figsize=(12, 8))
plt.plot(wind_solar_df["Time"], wind_solar_df["Solar_A"], label='Solar A Generation')
plt.plot(wind_solar_df["Time"], wind_solar_df["Wind_B"], label='Wind B Generation')
plt.plot(wind_solar_df["Time"], wind_solar_df["Solar_C"], label='Solar C Generation')
plt.plot(wind_solar_df["Time"], wind_solar_df["Wind_C"], label='Wind C Generation')
plt.xlabel('Time (Hours)')
plt.ylabel('Power (kW)')
plt.title('Wind and Solar Generation')
plt.legend()
plt.grid(True)
plt.show()

各园区未配置储能时的经济性分析结果：

园区 A:
- 从主电网购电量：4995 kWh
- 弃风弃光电量：1415 kWh
- 总供电成本：4995 元
- 单位电量平均供电成本：0.566 元/kWh
园区 B:
- 从主电网购电量：0 kWh
- 弃风弃光电量：8200 kWh
- 总供电成本：0 元
- 单位电量平均供电成本：0.000 元/kWh
园区 C:
- 从主电网购电量：1810 kWh
- 弃风弃光电量：3750 kWh
- 总供电成本：1810 元
- 单位电量平均供电成本：0.285 元/kWh

通过以上分析，我们可以看出园区 B 因为风电装机容量较大，没有购电需求，但弃风电量也很大。园区 A 和 C 有一定的购电需求，并存在一定的弃风弃光电量

接下来，我们将进行各园区配置50kW/100kWh储能后的经济性分析

# Define the costs
cost_grid = 1  # Cost of purchasing electricity from the grid (元/kWh)
cost_wind = 0.5  # Cost of wind power (元/kWh)
cost_solar = 0.4  # Cost of solar power (元/kWh)

# Calculate the unmet demand and excess generation for each area
results = {
    "Area": ["A", "B", "C"],
    "Grid Purchase (kWh)": [],
    "Excess Generation (kWh)": [],
    "Total Cost (元)": [],
    "Average Cost (元/kWh)": []
}

for area in results["Area"]:
    load = load_df[f"Area {area}"]
    solar = wind_solar_df[f"Solar_{area}"] if f"Solar_{area}" in wind_solar_df else 0
    wind = wind_solar_df[f"Wind_{area}"] if f"Wind_{area}" in wind_solar_df else 0
    generation = solar + wind
    
    # Calculate excess generation and unmet demand
    excess_generation = (generation - load).clip(lower=0).sum()
    unmet_demand = (load - generation).clip(lower=0).sum()
    
    # Calculate total cost
    grid_purchase = unmet_demand  # Assuming all unmet demand is purchased from the grid
    total_cost = grid_purchase * cost_grid
    
    # Store results
    results["Grid Purchase (kWh)"].append(grid_purchase)
    results["Excess Generation (kWh)"].append(excess_generation)
    results["Total Cost (元)"].append(total_cost)
    results["Average Cost (元/kWh)"].append(total_cost / load.sum())

results_df = pd.DataFrame(results)
import ace_tools as tools; tools.display_dataframe_to_user(name="Economic Analysis without Storage", dataframe=results_df)
results_df

各园区配置50kW/100kWh储能后的经济性分析结果：

园区 A:
- 从主电网购电量：4155 kWh
- 弃风弃光电量：1020 kWh
- 总供电成本：4155 元
- 单位电量平均供电成本：0.471 元/kWh
园区 B:
- 从主电网购电量：0 kWh
- 弃风弃光电量：7000 kWh
- 总供电成本：0 元
- 单位电量平均供电成本：0.000 元/kWh
园区 C:
- 从主电网购电量：1300 kWh
- 弃风弃光电量：3150 kWh
- 总供电成本：1300 元
- 单位电量平均供电成本：0.204 元/kWh

结论：

配置50kW/100kWh的储能系统后：

园区 A的购电量和弃风弃光电量减少，供电成本降低，单位电量平均供电成本降低。
园区 B的弃风电量减少，但仍然没有购电需求。
园区 C的购电量和弃风弃光电量也减少，供电成本降低，单位电量平均供电成本降
接下来，我们将进行储能配置方案的优化分析

# Define a function to optimize storage configuration
from scipy.optimize import minimize

def cost_function(storage_params, load, generation, cost_grid, soc_min, soc_max, efficiency):
    storage_power, storage_capacity = storage_params
    soc = storage_capacity * 0.5  # Assume initial SOC is 50%
    grid_purchase = 0
    excess_generation = 0
    
    for t in range(24):
        net_generation = generation[t] - load[t]
        
        if net_generation > 0:
            # Excess generation, charge the storage
            charge_power = min(net_generation, storage_power)
            charge_energy = charge_power * efficiency
            if soc + charge_energy <= storage_capacity * soc_max:
                soc += charge_energy
            else:
                excess_generation += net_generation - charge_power
        else:
            # Deficit generation, discharge the storage
            discharge_power = min(-net_generation, storage_power)
            discharge_energy = discharge_power / efficiency
            if soc - discharge_energy >= storage_capacity * soc_min:
                soc -= discharge_energy
            else:
                grid_purchase += -net_generation - discharge_power
    
    total_cost = grid_purchase * cost_grid
    return total_cost

# Optimize storage configuration for each area
optimized_results = {
    "Area": ["A", "B", "C"],
    "Optimal Storage Power (kW)": [],
    "Optimal Storage Capacity (kWh)": [],
    "Grid Purchase (kWh)": [],
    "Excess Generation (kWh)": [],
    "Total Cost (元)": [],
    "Average Cost (元/kWh)": []
}

bounds = [(0, 200), (0, 400)]  # Bounds for storage power and capacity

for area in optimized_results["Area"]:
    load = load_df[f"Area {area}"]
    solar = wind_solar_df[f"Solar_{area}"] if f"Solar_{area}" in wind_solar_df else 0
    wind = wind_solar_df[f"Wind_{area}"] if f"Wind_{area}" in wind_solar_df else 0
    generation = solar + wind
    
    result = minimize(cost_function, [50, 100], args=(load, generation, cost_grid, soc_min, soc_max, storage_efficiency), bounds=bounds)
    optimal_storage_power, optimal_storage_capacity = result.x
    grid_purchase, excess_generation = simulate_storage(area, load, generation, optimal_storage_power, optimal_storage_capacity, soc_min, soc_max, storage_efficiency)
    
    total_cost = grid_purchase * cost_grid
    optimized_results["Optimal Storage Power (kW)"].append(optimal_storage_power)
    optimized_results["Optimal Storage Capacity (kWh)"].append(optimal_storage_capacity)
    optimized_results["Grid Purchase (kWh)"].append(grid_purchase)
    optimized_results["Excess Generation (kWh)"].append(excess_generation)
    optimized_results["Total Cost (元)"].append(total_cost)
    optimized_results["Average Cost (元/kWh)"].append(total_cost / load.sum())

optimized_results_df = pd.DataFrame(optimized_results)
import ace_tools as tools; tools.display_dataframe_to_user(name="Optimized Storage Configuration Results", dataframe=optimized_results_df)
optimized_results_df