如何在IBM CPLEX中定义集合并编写配电网韧性优化约束

一、集合与数据结构定义 #

首先用OPL的tuple类型定义支路的节点信息，再基于此定义各类支路集合和时间集合：

// 定义支路结构：包含起始节点i和终止节点j
tuple Branch {
    int i;
    int j;
}

// 1. 所有支路集合（共32条）
{Branch} L = ...; 
// 2. 配备远程控制开关（RCS）的支路集合（共3条）
{Branch} L_switch = ...; 
// 3. 时间时段集合（根据研究时段范围定义，比如1..T）
{int} Time = ...; 
// 4. 各时段受损未修复的支路集合（随时间t变化，初始有9条受损）
{Branch} L_t_damaged[Time] = ...; 

二、参数与决策变量定义 #

定义初始状态参数和支路连接状态变量：

// 支路初始连接状态参数（二进制：1连通，0断开）
float alpha0[L] = ...; 
// 支路在t时刻的连接状态决策变量（二进制：1连通，0断开）
dvar boolean alpha[L][Time]; 

三、核心约束编写 #

针对未受损且无RCS的支路，通过forall循环结合集合过滤条件实现约束逻辑，对应公式αij,t = α0_ij ∀(i,j) ∈ L \ {L_t_damaged, L_switch}, ∀t ∈ T：

// 约束：未受损且无RCS的支路在所有时段保持初始状态
forall(t in Time, b in L: b not in L_switch && b not in L_t_damaged[t]) {
    alpha[b][t] == alpha0[b];
}

四、数据赋值示例（参考） #

在.dat文件中填充实际配电网数据，示例如下：

// 所有支路集合（示例，需替换为你的32条支路节点对）
L = {<1,2>, <1,3>, <2,4>, <3,4>, ...}; 
// 带RCS的支路集合（示例3条）
L_switch = {<5,6>, <7,8>, <9,10>}; 
// 时间时段（假设研究3个时段）
Time = {1,2,3}; 
// 各时段受损未修复支路（示例第1时段9条受损，后续时段可根据修复情况调整）
L_t_damaged[1] = {<11,12>, <13,14>, <15,16>, <17,18>, <19,20>, <21,22>, <23,24>, <25,26>, <27,28>};
L_t_damaged[2] = {<11,12>, <13,14>, ...}; // 第2时段未修复的支路（部分已修复）
L_t_damaged[3] = {...}; 
// 初始状态参数（对应每条支路的初始连通状态）
alpha0 = [1, 1, 0, 1, ...]; 

逻辑说明 #

约束中的过滤条件b not in L_switch && b not in L_t_damaged[t]等价于集合运算L \ {L_t_damaged, L_switch}，确保仅对既不在RCS支路集合、也不在当前时段受损未修复集合的支路，强制其各时段连接状态等于初始状态。

内容的提问来源于stack exchange，提问作者VANDANA KUMARI