0
点赞
收藏
分享

微信扫一扫

AI驱动的供应链管理:需求预测实战指南

AI驱动的供应链管理:需求预测实战指南

从统计模型到深度学习的完整落地路径

1. 需求预测为何决定供应链生死

1.1 牛鞭效应与成本放大

在传统供应链中,终端需求5%的波动会在上游被放大到40%的库存偏差(MIT Beer Game实验数据)。需求预测误差每降低1%,库存周转率可提升7-10%,缺货率下降3-5%。

1.2 传统方法的三大瓶颈
  • 移动平均/指数平滑:无法捕捉非线性与多变量关系
  • ARIMA:对季节性处理粗糙,高阶差分导致信息损失
  • 人工规则:依赖专家经验,无法适应SKU级别的动态变化

2. AI预测的技术跃迁

2.1 机器学习vs深度学习

维度

LightGBM

Temporal Fusion Transformer (TFT)

数据量

1万-100万样本

100万+样本

特征工程

需手工构造

自动提取时序特征

解释性

SHAP值

内置注意力可视化

冷启动问题

可迁移学习

需预训练

2.2 特征工程的质变

传统方法仅用历史销量,AI引入:

  • 外部特征:气象(温度/降水)、节假日、宏观经济指标
  • 交互特征:促销折扣与天气的交叉效应(晴天促销销量提升1.8倍)
  • 图特征:供应商-仓库的拓扑关系影响补货提前期

3. 数据工程:从脏数据到特征仓库

3.1 数据清洗的暗黑技巧

# 处理促销期间的异常峰值
def winsorize_promo(df, promo_col='is_promo', sales_col='sales'):
    promo_mask = df[promo_col] == 1
    q99 = df[promo_mask][sales_col].quantile(0.99)
    df.loc[promo_mask & (df[sales_col] > q99), sales_col] = q99
    return df

# 动态缺失值填充:相似SKU的同期均值
def dynamic_impute(df, sku_col='sku', date_col='date'):
    df['imputed'] = df.groupby([sku_col, df[date_col].dt.month])[sales_col].transform(
        lambda x: x.fillna(x.median())
    )
    return df

3.2 时间序列特征自动化

from tsfresh import extract_features
from tsfresh.feature_selection.relevance import calculate_relevance_table

# 自动提取300+时序特征
features = extract_features(
    df, 
    column_id='sku', 
    column_sort='date',
    default_fc_parameters=EfficientFCParameters()
)

# 基于假设检验筛选显著特征
relevance = calculate_relevance_table(features, df['demand'])
selected_features = relevance[relevance['p_value'] < 0.01]['feature'].tolist()

4. 模型实战:LightGBM与TFT双案例

4.1 LightGBM的供应链优化

import lightgbm as lgb
from sklearn.metrics import mean_absolute_percentage_error as mape

# 构造滞后特征(过去7/14/30天销量)
def create_lags(df, lags=[7,14,30]):
    for lag in lags:
        df[f'sales_lag_{lag}'] = df.groupby('sku')['sales'].shift(lag)
    return df

# 自定义损失函数:加权MAPE(近期误差权重更高)
def weighted_mape(preds, train_data):
    labels = train_data.get_label()
    weights = np.power(0.9, np.arange(len(labels))[::-1])  # 指数衰减权重
    return 'weighted_mape', mape(labels, preds, sample_weight=weights), False

# 训练
train_data = lgb.Dataset(X_train, label=y_train)
model = lgb.train(
    params={
        'objective': 'regression_l1',
        'learning_rate': 0.03,
        'max_depth': 8,
        'num_leaves': 64,
        'feature_fraction': 0.8,
        'bagging_fraction': 0.8,
        'bagging_freq': 1
    },
    train_set=train_data,
    feval=weighted_mape,
    num_boost_round=1000,
    early_stopping_rounds=100
)

4.2 TFT的注意力机制解析

from pytorch_forecasting import TemporalFusionTransformer, Baseline

# 定义时间序列数据集
dataset = TimeSeriesDataSet(
    data=df,
    time_idx="time_idx",
    target="demand",
    group_ids=["sku", "warehouse"],
    min_encoder_length=30,
    max_encoder_length=60,
    min_prediction_length=7,
    max_prediction_length=14,
    static_categoricals=["category", "region"],
    time_varying_known_reals=["discount", "temperature"],
    time_varying_unknown_reals=["demand"],
    add_relative_time_idx=True,
    add_target_scales=True,
    add_encoder_length=True,
)

# 训练
tft = TemporalFusionTransformer.from_dataset(
    dataset,
    learning_rate=0.01,
    hidden_size=32,
    attention_head_size=4,
    dropout=0.1,
    loss=QuantileLoss(),
)
trainer.fit(
    tft,
    train_dataloaders=train_dataloader,
    val_dataloaders=val_dataloader,
)

# 注意力可视化:识别关键驱动因子
interpretation = tft.interpret_output(batch, reduction="sum")
tft.plot_interpretation(interpretation)

5. 模型解释性:让业务方信服

5.1 SHAP值揭示促销失效场景

import shap

# 计算全局特征重要性
explainer = shap.TreeExplainer(model)
shap_values = explainer.shap_values(X_test)

# 发现:当折扣>30%且温度<5℃时,促销反而抑制需求(shap值<-0.3)
shap.dependence_plot(
    "discount", 
    shap_values, 
    X_test,
    interaction_index="temperature"
)

5.2 TFT的变量选择网络

TFT的静态变量选择权重显示:

  • 品类权重:生鲜(0.82) > 日化(0.45) > 家电(0.21)
  • 区域权重:一线城市(0.75) > 三线以下(0.33)

6. 部署与实时推断

6.1 流式特征管道

from kafka import KafkaConsumer
import faust

# 实时温度数据接入
app = faust.App('feature-pipeline', broker='kafka://localhost:9092')
temperature_topic = app.topic('weather', value_type=float)

@app.agent(temperature_topic)
async def process_temperature(temperatures):
    async for temp in temperatures:
        redis_client.set(f"temp_{datetime.now().date()}", temp)

6.2 边缘计算优化

将LightGBM模型转换为ONNX格式,部署在门店边缘设备:

import onnxmltools
onnx_model = onnxmltools.convert_lightgbm(model)
onnxmltools.utils.save_model(onnx_model, 'demand_predict.onnx')

7. 效果评估与持续迭代

7.1 供应链KPI闭环
  • 预测准确率:MAPE从23%降至11%
  • 库存周转:核心SKU周转天数从45天降至28天
  • 缺货率:大促期间从8.5%降至2.3%
7.2 模型漂移检测

使用Kolmogorov-Smirnov检验监控特征分布:

from scipy.stats import ks_2samp

# 检测温度特征漂移
_, p_value = ks_2samp(
    reference_data['temperature'], 
    current_data['temperature']
)
if p_value < 0.05:
    trigger_model_retrain()

举报

相关推荐

0 条评论