供应链计划篇:实际案例计划精解
摘自《供应链管理:战略、规划与运作》(Supply Chain Management: Strategy, Planning, and Operation)
作者:Sunil Chopra 和 Peter Meindl
在供应链管理实践中,如何在需求波动显著的环境下制定成本最优的生产与库存计划,是企业面临的关键挑战之一。以“红西红柿”公司为例,这家园艺工具制造商在春季面临需求高峰,为确保供应响应能力和控制运营成本,必须在雇佣、解雇、加班、转包与库存之间做出权衡。
在供应链管理中,总体计划(Aggregate Planning)是连接长期战略规划与短期运营执行的桥梁,尤其适用于应对季节性需求、产能波动及资源配置问题。本文以“红西红柿”(Red Tomato's)公司为例,系统分析总体计划模型的构建过程、目标函数、限制条件及应对不确定性的优化策略,并进一步探讨在成本和需求变化下的敏感性分析。
一、问题背景与目标概述
“红西红柿”(Red Tomato's)公司是一家园艺工具生产商,生产流程中主要受到劳动力资源的限制。其产品具有强烈的季节性需求高峰,尤其在春季(3月-4月)需求激增。面对波动的市场需求,公司需要通过总体计划(Aggregate Planning)协调劳动力、库存、转包等资源,以最小化总成本,同时保证一定的服务水平和末期库存。计划周期设定为6个月,公司需保证:
- 1月初库存1000单位;
- 初始劳动力人数为80人;
- 最终(6月末)库存不少于500单位;
- 所有需求都应满足,允许短期积压。
| 1月 | 2月 | 3月 | 4月 | 5月 | 6月 | |
|---|---|---|---|---|---|---|
| 需求 | 1600 | 3000 | 3200 | 3800 | 2200 | 2200 |
二、模型设定
2.1 决策变量
令 $ t=1,2,...,6 $ 表示1月到6月,定义如下变量:
- $D_t$: \(t\)月的需求量
- $W_t$: \(t\)月的工人数
- $H_t$: \(t\)月新雇佣工人数
- $L_t$: \(t\)月解雇工人数
- $P_t$: \(t\)月生产量
- $I_t$: \(t\)月期末库存
- $S_t$: \(t\)月期末积压(缺货)量
- $C_t$: \(t\)月转包产量
- $O_t$: \(t\)月加班总小时数
2.2 参数设定
- 单位售价:40美元(不影响成本最小化)
- 原材料成本:10美元/单位
- 库存成本:2美元/单位/月(例2改为6美元)
- 缺货成本:5美元/单位/月
- 雇佣成本:300美元/人
- 解雇成本:500美元/人
- 正常工资:4美元/小时 × 8小时 × 20天 = 640美元/人/月
- 加班工资:6美元/小时
- 每个工人每月最多加班10小时
- 每单位产品生产需要4小时
2.3 目标函数(成本最小化)
2.4 约束条件
- 工人数量变动约束:
- 生产能力约束(40单位/工人/月 + 加班):
- 库存平衡约束:
- 加班时间限制:
- 非负约束:所有变量均大于等于0;$W_t$可取实数估计人力规模。
整合后的模型为
2.5 📦 平均库存与平均周转时间
在供应链总体计划中,平均库存(Average Inventory) 和 平均周转时间(Average Inventory Turnover Time) 是衡量库存效率和资金占用水平的重要指标。它们的计算依赖于各期的期末库存水平。
平均库存计算公式
若计划周期为 $ T $ 期,假定每期的平均库存量取为该期期初与期末库存的平均,即:
则整个计划周期的平均库存为:
其中:
- $ I_0 $:初始库存
- $ I_T $:最后一期(第 $T $ 期)末库存
- $ I_t $:第 $ t $ 期末库存(\(1 ≤ t ≤ T-1\))
平均库存周转时间计算公式
平均库存周转时间衡量库存单位从投入到售出的平均时间长度,计算公式如下:
其中:
- 分子为计划期内的平均库存;
- 分母是根据最小法则定义的平均需求率 = \(\frac{1}{T} \sum_{t=1}^{T} D_t\);
- 周转时间单位为“期”(如“月”)。
指标应用说明
- 平均库存越小,说明企业库存管理效率越高;
- 平均周转时间越短,表示库存周转速度更快,资金占用更少;
- 在实际计划优化中,这两个指标通常与服务水平、安全库存和成本目标协调考虑。
--
三、Python建模求解实现
import pulp
# 初始化模型
model = pulp.LpProblem("Aggregate_Planning", pulp.LpMinimize)
# 时间周期
T = 6
months = range(1, T + 1)
# 需求数据
D = {1: 1600, 2: 3000, 3: 3200, 4: 3800, 5: 2200, 6: 2200}
# 参数
wage = 640
ot_cost = 6
hire_cost = 300
layoff_cost = 500
inv_cost = 2
backlog_cost = 5
prod_cost = 10
subcontract_cost = 30
# 初始状态
I0 = 1000
S0 = 0
W0 = 80
# 定义变量
W = pulp.LpVariable.dicts("W", months, lowBound=0)
H = pulp.LpVariable.dicts("H", months, lowBound=0)
L = pulp.LpVariable.dicts("L", months, lowBound=0)
P = pulp.LpVariable.dicts("P", months, lowBound=0)
C = pulp.LpVariable.dicts("C", months, lowBound=0)
O = pulp.LpVariable.dicts("O", months, lowBound=0)
I = pulp.LpVariable.dicts("I", range(0, T + 1), lowBound=0)
S = pulp.LpVariable.dicts("S", range(0, T + 1), lowBound=0)
# 初始库存和积压设置
I[0] = I0
S[0] = S0
W_prev = W0
# 目标函数
model += pulp.lpSum([
wage * W[t] + ot_cost * O[t] + hire_cost * H[t] + layoff_cost * L[t] +
inv_cost * I[t] + backlog_cost * S[t] + prod_cost * P[t] + subcontract_cost * C[t]
for t in months
])
# 约束条件
for t in months:
if t == 1:
model += W[t] == W_prev + H[t] - L[t]
else:
model += W[t] == W[t - 1] + H[t] - L[t]
model += O[t] <= 10 * W[t]
model += P[t] <= 40 * W[t] + 0.25 * O[t]
model += I[t - 1] + P[t] + C[t] == D[t] + S[t - 1] + I[t] - S[t]
# 末期库存与积压要求
model += I[T] >= 500
model += S[T] == 0
# 求解模型
model.solve()
# 输出总成本
total_cost = pulp.value(model.objective)
print(f"总成本:{total_cost:.2f}")
# 输出每月详细信息
for t in months:
print(f"第{t}月: 工人={int(W[t].varValue)}, "
f"雇用={int(H[t].varValue)}, 解雇={int(L[t].varValue)}, "
f"产量={int(P[t].varValue)}, 转包={int(C[t].varValue)}, "
f"加班={int(O[t].varValue)}, 库存={int(I[t].varValue)}, "
f"积压={int(S[t].varValue)}")
# 平均库存计算
I_vals = [I[0]] + [I[t].varValue for t in months]
avg_inventory = ((I_vals[0] + I_vals[T]) / 2 + sum(I_vals[1:T])) / T
print(f"平均库存:{avg_inventory:.2f} 单位")
# 平均需求
avg_demand = sum(D[t] for t in months) / T
# 平均库存周转时间
turnover_time = avg_inventory / avg_demand
print(f"平均库存周转时间:{turnover_time:.4f} 月")
总成本:422275.00
第1月: 工人=64, 雇用=0, 解雇=15, 产量=2583, 转包=0, 加班=0, 库存=1983, 积压=0
第2月: 工人=64, 雇用=0, 解雇=0, 产量=2583, 转包=0, 加班=0, 库存=1566, 积压=0
第3月: 工人=64, 雇用=0, 解雇=0, 产量=2583, 转包=0, 加班=0, 库存=950, 积压=0
第4月: 工人=64, 雇用=0, 解雇=0, 产量=2583, 转包=0, 加班=0, 库存=0, 积压=266
第5月: 工人=64, 雇用=0, 解雇=0, 产量=2583, 转包=0, 加班=0, 库存=116, 积压=0
第6月: 工人=64, 雇用=0, 解雇=0, 产量=2583, 转包=0, 加班=0, 库存=500, 积压=0
平均库存:894.44 单位
平均库存周转时间:0.3354 月
四、例1:高季节性需求影响分析
在该例中,保持总需求16000单位不变,将其分配如下:
| 1月 | 2月 | 3月 | 4月 | 5月 | 6月 | |
|---|---|---|---|---|---|---|
| 需求 | 1000 | 3000 | 3800 | 4800 | 2000 | 1400 |
4.1 Python建模求解(季节性大)
只需替换需求数组 D:
import pulp
# 初始化模型
model = pulp.LpProblem("Aggregate_Planning", pulp.LpMinimize)
# 时间周期
T = 6
months = range(1, T + 1)
# ✅ 新需求数据
D = {1: 1000, 2: 3000, 3: 3800, 4: 4800, 5: 2000, 6: 1400}
# 参数
wage = 640
ot_cost = 6
hire_cost = 300
layoff_cost = 500
inv_cost = 2
backlog_cost = 5
prod_cost = 10
subcontract_cost = 30
# 初始状态
I0 = 1000
S0 = 0
W0 = 80
# 定义变量
W = pulp.LpVariable.dicts("W", months, lowBound=0)
H = pulp.LpVariable.dicts("H", months, lowBound=0)
L = pulp.LpVariable.dicts("L", months, lowBound=0)
P = pulp.LpVariable.dicts("P", months, lowBound=0)
C = pulp.LpVariable.dicts("C", months, lowBound=0)
O = pulp.LpVariable.dicts("O", months, lowBound=0)
I = pulp.LpVariable.dicts("I", range(0, T + 1), lowBound=0)
S = pulp.LpVariable.dicts("S", range(0, T + 1), lowBound=0)
# 初始库存和积压设置
I[0] = I0
S[0] = S0
W_prev = W0
# 目标函数
model += pulp.lpSum([
wage * W[t] + ot_cost * O[t] + hire_cost * H[t] + layoff_cost * L[t] +
inv_cost * I[t] + backlog_cost * S[t] + prod_cost * P[t] + subcontract_cost * C[t]
for t in months
])
# 约束条件
for t in months:
if t == 1:
model += W[t] == W_prev + H[t] - L[t]
else:
model += W[t] == W[t - 1] + H[t] - L[t]
model += O[t] <= 10 * W[t]
model += P[t] <= 40 * W[t] + 0.25 * O[t]
model += I[t - 1] + P[t] + C[t] == D[t] + S[t - 1] + I[t] - S[t]
# 末期库存与积压要求
model += I[T] >= 500
model += S[T] == 0
# 求解模型
model.solve()
# 输出总成本
total_cost = pulp.value(model.objective)
print(f"总成本:{total_cost:.2f}")
# 输出每月详细信息
for t in months:
print(f"第{t}月: 工人={int(W[t].varValue)}, "
f"雇用={int(H[t].varValue)}, 解雇={int(L[t].varValue)}, "
f"产量={int(P[t].varValue)}, 转包={int(C[t].varValue)}, "
f"加班={int(O[t].varValue)}, 库存={int(I[t].varValue)}, "
f"积压={int(S[t].varValue)}")
# 平均库存计算
I_vals = [I[0]] + [I[t].varValue for t in months]
avg_inventory = ((I_vals[0] + I_vals[T]) / 2 + sum(I_vals[1:T])) / T
print(f"平均库存:{avg_inventory:.2f} 单位")
# 平均需求
avg_demand = sum(D[t] for t in months) / T
# 平均库存周转时间
turnover_time = avg_inventory / avg_demand
print(f"平均库存周转时间:{turnover_time:.4f} 月")
4.2 结果分析:
总成本:432858.33
第1月: 工人=64, 雇用=0, 解雇=15, 产量=2583, 转包=0, 加班=0, 库存=2583, 积压=0
第2月: 工人=64, 雇用=0, 解雇=0, 产量=2583, 转包=0, 加班=0, 库存=2166, 积压=0
第3月: 工人=64, 雇用=0, 解雇=0, 产量=2583, 转包=0, 加班=0, 库存=950, 积压=0
第4月: 工人=64, 雇用=0, 解雇=0, 产量=2583, 转包=0, 加班=0, 库存=0, 积压=1266
第5月: 工人=64, 雇用=0, 解雇=0, 产量=2583, 转包=0, 加班=0, 库存=0, 积压=683
第6月: 工人=64, 雇用=0, 解雇=0, 产量=2583, 转包=0, 加班=0, 库存=500, 积压=0
平均库存:1075.00 单位
平均库存周转时间:0.4031 月
- 工人规模维持不变,依靠库存和短期积压完成供需协调。
- 总成本从原来的422275美元升至432858美元。
- 平均周转时间增加为0.40个月,库存管理压力增大。
五、例2:库存成本升高的情景模拟
将库存持有成本由2美元提高至6美元。
5.1 Python建模求解(库存成本上升)
在目标函数中调整库存成本系数:
import pulp
# 初始化模型
model = pulp.LpProblem("Aggregate_Planning", pulp.LpMinimize)
# 时间周期
T = 6
months = range(1, T + 1)
# ✅ 新需求数据
D = {1: 1600, 2: 3000, 3: 3200, 4: 3800, 5: 2200, 6: 2200}
# 参数(库存持有成本提高为6)
wage = 640
ot_cost = 6
hire_cost = 300
layoff_cost = 500
inv_cost = 6 # ✅ 修改在此(原为2)
backlog_cost = 5
prod_cost = 10
subcontract_cost = 30
# 初始状态
I0 = 1000
S0 = 0
W0 = 80
# 定义变量
W = pulp.LpVariable.dicts("W", months, lowBound=0)
H = pulp.LpVariable.dicts("H", months, lowBound=0)
L = pulp.LpVariable.dicts("L", months, lowBound=0)
P = pulp.LpVariable.dicts("P", months, lowBound=0)
C = pulp.LpVariable.dicts("C", months, lowBound=0)
O = pulp.LpVariable.dicts("O", months, lowBound=0)
I = pulp.LpVariable.dicts("I", range(0, T + 1), lowBound=0)
S = pulp.LpVariable.dicts("S", range(0, T + 1), lowBound=0)
# 初始库存和积压设置
I[0] = I0
S[0] = S0
W_prev = W0
# 目标函数
model += pulp.lpSum([
wage * W[t] + ot_cost * O[t] + hire_cost * H[t] + layoff_cost * L[t] +
inv_cost * I[t] + backlog_cost * S[t] + prod_cost * P[t] + subcontract_cost * C[t]
for t in months
])
# 约束条件
for t in months:
if t == 1:
model += W[t] == W_prev + H[t] - L[t]
else:
model += W[t] == W[t - 1] + H[t] - L[t]
model += O[t] <= 10 * W[t]
model += P[t] <= 40 * W[t] + 0.25 * O[t]
model += I[t - 1] + P[t] + C[t] == D[t] + S[t - 1] + I[t] - S[t]
# 末期库存与积压要求
model += I[T] >= 500
model += S[T] == 0
# 求解模型
model.solve()
# 输出总成本
total_cost = pulp.value(model.objective)
print(f"总成本:{total_cost:.2f}")
# 输出每月详细信息
for t in months:
print(f"第{t}月: 工人={int(W[t].varValue)}, "
f"雇用={int(H[t].varValue)}, 解雇={int(L[t].varValue)}, "
f"产量={int(P[t].varValue)}, 转包={int(C[t].varValue)}, "
f"加班={int(O[t].varValue)}, 库存={int(I[t].varValue)}, "
f"积压={int(S[t].varValue)}")
# 平均库存计算
I_vals = [I[0]] + [I[t].varValue for t in months]
avg_inventory = ((I_vals[0] + I_vals[T]) / 2 + sum(I_vals[1:T])) / T
print(f"平均库存:{avg_inventory:.2f} 单位")
# 平均需求
avg_demand = sum(D[t] for t in months) / T
# 平均库存周转时间
turnover_time = avg_inventory / avg_demand
print(f"平均库存周转时间:{turnover_time:.4f} 月")
总成本:441200.00
第1月: 工人=56, 雇用=0, 解雇=23, 产量=2266, 转包=0, 加班=0, 库存=1666, 积压=0
第2月: 工人=56, 雇用=0, 解雇=0, 产量=2266, 转包=0, 加班=0, 库存=933, 积压=0
第3月: 工人=56, 雇用=0, 解雇=0, 产量=2266, 转包=0, 加班=0, 库存=0, 积压=0
第4月: 工人=56, 雇用=0, 解雇=0, 产量=2266, 转包=1466, 加班=0, 库存=0, 积压=66
第5月: 工人=56, 雇用=0, 解雇=0, 产量=2266, 转包=0, 加班=0, 库存=0, 积压=0
第6月: 工人=56, 雇用=0, 解雇=0, 产量=2266, 转包=433, 加班=0, 库存=500, 积压=0
平均库存:558.33 单位
平均库存周转时间:0.2094 月
5.2 结果分析:
- 为降低库存持有成本,公司更倾向于通过解雇、转包等方式减少库存。
- 总体计划总成本 = 441200美元,低于沿用旧方案时的442742美元。
- 平均季节性库存下降至558单位,平均周转时间降至0.21个月。
六、三案例横向比较(表格与指标分析)
前面介绍了“红西红柿”公司在不同需求和成本条件下的总体计划模型求解结果,分别是基准需求(例0)、高波动需求(例1)和库存成本提升(例2)三种典型情景。为了更直观地理解不同情景对企业生产运营的影响,本节将基于这些案例,从关键运营指标层面对比分析,帮助管理者把握成本结构、库存水平和服务水平之间的平衡。
6.1 关键指标选取
选取以下关键指标进行对比:
| 指标名称 | 说明 |
|---|---|
| 总成本(美元) | 6个月计划期内的全部成本总和,包括劳动力、库存、缺货等成本 |
| 平均库存(单位) | 计划期内平均持有库存量,反映库存资金占用及潜在积压风险 |
| 库存缺货(累计单位) | 计划期内的累计缺货或积压量,体现服务水平及客户满意度 |
| 加班总小时数(小时) | 计划期内的累计加班小时,反映生产弹性及劳动力调度强度 |
| 转包总产量(单位) | 计划期内通过外包完成的产量,体现对内部产能不足的依赖程度 |
| 平均周转时间(月) | 反映库存周转速度,数值越小表示资金流动越快,库存管理越高效 |
6.2 三案例指标对比表
| 指标 | 初始案例(基准需求) | 例1(高波动需求) | 例2(库存成本提升) |
|---|---|---|---|
| 总成本(美元) | 422,275 | 432,858 | 441,200 |
| 平均库存(单位) | 895 | 1,075 | 558 |
| 库存缺货(累计单位) | 0 | 1,950 | 1,467 |
| 加班总小时数(小时) | 0 | 0 | 0 |
| 转包总产量(单位) | 0 | 0 | 433 |
| 平均周转时间(月) | 0.34 | 0.40 | 0.21 |
6.3 分析解读
-
成本层面
基准需求案例成本最低,表明生产与需求较为匹配,库存和缺货都控制良好。高波动需求导致库存缺货增加,虽然加班和转包未启用,但为保持产能稳定和避免服务水平下降,库存水平明显上升,推高了整体成本。库存成本提升案例中,企业采取减少库存并增加转包产量的策略,转包成本高于内部生产,但总成本仍略高于基准案例,显示出通过调整生产方式以适应成本结构变化的重要性。 -
库存与缺货
高波动需求下,库存和缺货双双上升,反映需求不确定性给供应链带来的压力。库存成本提升情景则有效降低库存量,但同时引入部分缺货,说明企业在成本压力下更倾向于以转包和适度缺货换取库存资金减少。 -
加班与转包
三个案例均未采用加班策略,说明在该模型和成本结构下,加班不具备成本优势。转包则在库存成本较高的情况下被激活,成为弥补内部产能不足的重要手段。 -
平均周转时间
周转时间反映库存资金占用效率。高波动需求导致周转时间最长,增加了资金压力。库存成本提升案例周转时间最短,表明减少库存有助于提升资金周转效率。
三案例横向比较揭示了需求波动和成本变化对生产计划的深刻影响。管理者应根据企业具体运营环境,合理调整库存水平、产能配置和外包策略,以平衡成本、服务和风险。动态的总体计划工具使企业能快速响应市场变化,避免盲目扩产或库存积压,提高整体供应链效率和竞争力。
七、不确定性应对策略
在实际供应链管理中,预测误差和需求不确定性是不可避免的风险因素。尽管总体计划模型能够基于预测进行最优化,但单纯依赖预测值很可能导致服务水平波动、库存成本激增或生产资源浪费。因此,企业需要引入多种策略,增强计划的鲁棒性和适应能力。
7.1 预测误差与风险来源
- 需求预测误差:销售数据、市场环境变化、新产品推广等均可能导致需求预测偏差。
- 供应风险:原材料延迟、供应商能力不足、运输中断等影响生产连续性。
- 产能不确定性:设备故障、工人流动、劳动法规变化等影响有效产能。
- 外部环境波动:政策调整、自然灾害、突发事件等均会扰动计划执行。
7.2 安全库存(Safety Stock)
安全库存是在预测需求基础上额外保留的一定量库存,用以缓冲需求和供应的波动,保证较高的客户服务水平。
7.2.1 计算方法
安全库存常基于需求标准差和服务水平(置信度)计算:
- \(S S\):安全库存量
- \(Z\):标准正态分布对应的服务水平系数(例如95%服务水平对应1.65)
- \(\sigma_{D}\):需求标准差
- \(L\):提前期(lead time)长度(期数)
7.2.2 应用要点
- 通过增加安全库存,企业能在需求突然增加时快速响应,避免缺货和客户流失。
- 安全库存增加会提升库存持有成本,因此需权衡服务水平和资金占用。
- 不同产品和供应链环节的安全库存应区别对待,重点保障关键产品或核心供应商。
7.3 安全产能(Safety Capacity)
安全产能指企业在正常产能基础上保留的额外生产能力,用以应对短期产能突增需求或生产中断。
7.3.1 实现方式
- 加班工时:安排弹性加班以应对短期峰值需求。
- 临时雇佣:通过临时工或短期劳务合同增加劳动力。
- 转包生产:将部分生产外包给第三方,缓解内部产能瓶颈。
- 多线生产:增加生产线冗余,保证设备维护时不影响产能。
- 灵活生产工艺:采用通用设备和工艺快速切换不同产品生产。
7.3.2 优劣权衡
- 保持安全产能提高供应链灵活性,缩短响应时间。
- 但增加产能会导致固定成本和管理复杂度上升。
- 过度依赖加班和临时劳动力可能影响员工士气和生产质量。
7.4 综合策略与总体计划结合
在总体计划模型中,可通过以下方式引入安全库存和安全产能:
- 调整库存约束:将安全库存作为最低库存要求,确保库存水平不低于安全阈值。
- 扩展产能约束:允许加班时数和转包量上限根据安全产能需求灵活调整。
- 场景模拟:基于不同预测误差和需求波动情景,反复求解总体计划,确定最优安全缓冲水平。
- 动态调整:结合滚动计划理念,周期性更新安全库存和产能参数,响应市场和运营变化。
7.5 现实应用案例
- 制造业:汽车厂商通常保持多日安全库存及多层供应商转包能力,保障产线连续性。
- 零售业:季节性商品设置较高安全库存,应对促销和节日需求波动。
- 电子行业:高科技产品加大安全产能,通过灵活用工和外包应对市场不确定。
面对复杂多变的供应链环境,仅依赖精准预测难以保证运营稳健。通过合理设计安全库存和安全产能,企业能有效缓冲不确定性,保障客户需求满足和生产平稳。结合动态总体计划模型,安全策略成为供应链管理中的关键杠杆,有助于实现成本控制与高服务水平的平衡,提升整体竞争力和抗风险能力。
八、管理启示总结
总体计划模型不仅为“红西红柿”公司提供了应对季节性需求波动的科学方法,还通过数学建模将生产、库存、劳动力和外包等多个因素系统化整合,实现成本最优化与服务水平平衡。通过不同需求和成本情景下的案例分析,模型展示了供应链管理中资源配置与风险权衡的重要性,为管理者提供了决策支持的量化依据。
此外,模型的灵活性允许根据企业实际情况调整约束和参数,支持多样化的生产策略,如加班、转包和库存管理等,为企业提升供应链韧性和响应速度提供可能。未来,可结合预测误差与不确定性因素,进一步引入安全库存、安全生产能力等缓冲机制,提升总体计划的鲁棒性和适应性。同时,结合机器学习预测与动态优化技术,将模型应用于更复杂、多变的供应链环境,实现智能化、自动化的运营决策支持。此类扩展将极大推动供应链管理理论与实践的深度融合,助力企业在激烈竞争中保持领先。
作者注:本文所有Python代码均可在
PuLP优化包中运行,建议搭配Jupyter Notebook使用,便于逐步调试与结果分析。

浙公网安备 33010602011771号