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 数据清洗的暗黑技巧

python 复制代码

# 处理促销期间的异常峰值
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 时间序列特征自动化

python 复制代码

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的供应链优化

python 复制代码

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的注意力机制解析

python 复制代码

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值揭示促销失效场景

python 复制代码

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 流式特征管道

python 复制代码

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格式，部署在门店边缘设备：

python 复制代码

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检验监控特征分布：

python 复制代码

from scipy.stats import ks_2samp

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