0 引 言
1 集成框架与边界设定
1.1 集成框架
1.2 变量与映射关系
表1 植物工厂经济脆弱性诊断的变量与映射关系表Table 1 Variable definitions and mapping relationships for economic vulnerability assessment of plant factories with artificial lighting |
| 模块 | 核心变量 | 操作性定义 | 主要映射关系 | 植物工厂经济脆弱性诊断中的角色 |
|---|---|---|---|---|
| 生物学层 | , | 作物原型与栽培管理 | (1) | 决定毛产量上限 |
| 生物学层 | , | 商品率与品质一致性 | (2) | 将产品转化为商品 |
| 数字-物理控制层 | 控制成熟度 | 、 、 、 均为 的函数 | 连接控制升级与经济结果 | |
| 数字-物理控制层 | , , , | 补光、温湿调控、自动化和数字化 | (3) (4) | 决定能耗与劳动替代 |
| 技术经济层 | (元/m2), (%), (年), (%) | 单位资本性支出(capital expenditure, CAPEX)、贴现率、寿命期、维护率 | (5) (6) | 构成资本年化压力 |
| 技术经济层 | 元/kWh, (元/(人·年)), (元/kg), (元/kg) | 电价、工资、耗材与其他运营成本 | , , | 构成OPEX主体 |
| 技术经济层 | (元/kg), (元/kg), (元/(m2·年)) | 售价、单位成本、单位面积利润 | (7) (8) | 识别盈亏边界 |
| 脆弱性层 | , , | 盈利缺口、能源暴露度、碳约束暴露度 | (9) | 对情景稳健性排序 |
| 输出层 | TP | 转型路径类型 | 结果反馈 | 形成效率优先、能源耦合、价值升级和综合转型四类路径 |
|
1.3 研究对象与系统边界
1.4 数据来源与参数标定
表2 植物工厂生菜经济脆弱性模型参数设定与说明Table 2 Parameter settings and description for the economic vulnerability model of plant factories with artificial lighting lettuce |
| 参数 | 符号 | 取值 | 来源/说明 |
|---|---|---|---|
| 基础控制下毛产量/(kg/(m2·年)) | 80 | 介于元分析年化代理值(62.8)与先进商业值(115)之间 | |
| 控制成熟度对毛产量提升系数 | 0.437 5 | 使C=1时达到115 kg/(m2·年) | |
| 基础控制下商品率/% | 0.88 | 闭环智能情景提升至 0.95 | |
| 控制成熟度对商品率提升系数 | 0.07 | 模型设定,使C=1时达到0.95 | |
| 基础控制下单位商品电耗/(kWh/kg) | 12.0 | 位于公开文献当前区间内 | |
| 控制成熟度对单位商品电耗降低系数 | 0.30 | 使C=1时 kWh/kg | |
| 基准单位CAPEX/(元/m2) | 8 000 | 位于公开PFAL投资区间内部 | |
| 控制升级导致CAPEX增加系数 | 0.15 | 模型设定,反映自动化与控制升级成本 | |
| 贴现率/% | 5 | 沿用PFAL规模经济研究设定[23] | |
| 寿命期/年 | 15 | 沿用PFAL规模经济研究设定[23] | |
| 维护费率/% | 1.5 | 沿用PFAL规模经济研究设定[23] | |
| 生产制造及有关人员平均工资/(元/(人·年)) | 78 561 | 国家统计局2024年数据 | |
| 基础人工需求系数/(人·年/m2) | 1/250 | 模型设定,人工成本核算 | |
| 基础自动化替代率 | 0.10 | 模型设定,人工替代基准 | |
| 控制成熟度对自动化提升系数 | 0.35 | 使C=1时β=0.45 | |
| 营养液、种子和卫生材料成本/(元/kg) | 1.6 | 模型设定 | |
| 其他可变运营成本 | 0.4元/kg | 模型设定 |
2 情景设计与核心模型
2.1 控制成熟度定义
2.2 能源及市场模型情景
2.3 情景组合设计
表3 植物工厂生菜控制、能源、市场情景组合与取值依据Table 3 Control-Energy-Market scenario combinations and value basis for PFAL lettuce |
| 维度 | 符号 | 取值 | 含义 | 设定依据 |
|---|---|---|---|---|
| 控制成熟度 | C0 | 0.0 | 基础控制 | 规则控制 有限自动化 |
| 控制成熟度 | C1 | 0.5 | 增强控制 | 传感、执行与调度增强 |
| 控制成熟度 | C2 | 1.0 | 闭环智能 | 闭环调控 较高自动化 |
| 能源情景 | E1 | 元/kWh kg CO2/kWh | 能源富裕型 | 低碳低价电力锚点 |
| 能源情景 | E2 | 元/kWh; kg CO2/kWh | 常规电网型 | 全国平均电网碳因子锚点 |
| 能源情景 | E3 | 元/kWh; kg CO2/kWh | 高价高碳型 | 不利能源条件 |
| 市场情景 | M1 | 元/kg | 常规鲜食市场 | 官方周均蔬菜批发价上沿 |
| 市场情景 | M2 | 元/kg | 品牌溢价市场 | 中等价格阈值情景 |
| 市场情景 | M3 | 元/kg | 高值直供市场 | 高价格阈值情景 |
2.4 核心模型
2.4.1 作物生产模块
2.4.2 能源与自动化模块
2.4.3 年化资本成本
2.4.4 运营成本与利润函数
2.4.5 碳强度与经济脆弱性
2.5 模型运行与结果输出
2.6 局部归一化弹性的定义与计算
2.7 政策支持与渠道成本扩展情景
3 结果与分析
3.1 控制成熟度对产量、成本与利润的影响
表4 E2+M3情景下不同控制成熟度的产量、成本、利润与碳强度结果Table 4 Yield, cost, profit and carbon intensity under different control maturity levels in the E2+M3 Scenario |
| 控制情景 | 毛产量 /(kg/(m2·年)) | 商品率 | 商品产量 /(kg/(m2·年)) | 单位商品电耗 /(kWh/kg) | (元/kg) 单位成本 / | 利润 /(元/(m2·年)) | 碳强度 /(kg CO2/kg) |
|---|---|---|---|---|---|---|---|
| 基础控制 | 80.0 | 0.880 | 70.40 | 12.0 | 27.67 | -258.36 | 6.37 |
| 增强控制 | 97.5 | 0.915 | 89.21 | 10.2 | 22.94 | 94.83 | 5.41 |
| 闭环智能 | 115.0 | 0.950 | 109.25 | 8.4 | 19.26 | 518.04 | 4.46 |
3.2 盈亏平衡售价与电价边界
3.3 27个情景组合下的利润空间
3.4 27个情景组合下的经济脆弱性
3.5 敏感度与转型路径识别
表5 基于利润、EVI与局部弹性的植物工厂生菜结果驱动的转型路径识别Table 5 Transformation path identification for plant factories with artificial lighting lettuce based on profit, EVI and local elasticity |
| 典型情景 | 结果特征 | 主导约束 | 对应路径 | 优先措施 |
|---|---|---|---|---|
| 增强控制+常规电网型+高值直供市场 | 利润已转正但边际有限,EVI中等 | 商品产量、资本强度 | 效率优先型 | 提高商品率、降低单位CAPEX、优化作业节拍 |
| 闭环智能+高价高碳型+高值直供市场 | 利润为正但能源暴露度和碳约束暴露度偏高 | 电价、电网碳因子 | 能源耦合型 | 分时电价响应、绿电采购、储能或源网荷储协同 |
| 闭环智能+任意能源+品牌/常规市场 | 技术成熟但售价不足以覆盖总成本 | 市场实现 | 价值升级型 | 品牌化、净菜化、活体供应、合同直供 |
| 基础控制+常规鲜食市场+常规/高价高碳电网 | 利润和EVI均处于不利区间 | 控制、能源、市场三重约束 | 综合转型型 | 同步升级控制、能源获取与销售模式 |
3.6 政策支持与渠道成本敏感性分析
表6 政策支持与渠道成本对E2+M3代表性情景利润的影响Table 6 Effects of policy support and channel cost on profits under the E2+M3 representative scenario |
| 扩展情景 | 基础控制C0/(元/(m2·年)) | 增强控制C1/(元/(m2·年)) | 闭环智能C2/(元/(m2·年)) |
|---|---|---|---|
| 基准情景 | -258.36 | 94.83 | 518.04 |
| 投资补贴20% | -104.21 | 260.54 | 695.31 |
| 电价优惠20% | -131.64 | 231.32 | 655.70 |
| 投资补贴20%+电价优惠20% | 22.51 | 397.03 | 832.97 |
| 渠道成本增加1元/kg | -328.76 | 5.62 | 408.79 |
| 渠道成本增加2元/kg | -399.16 | -83.60 | 299.54 |
| 渠道成本增加3元/kg | -469.56 | -172.81 | 190.29 |
| 投资补贴20%+电价优惠20%+渠道成本增加1元/kg | -47.89 | 307.82 | 723.72 |
3.7 中国情境下的典型案例对照分析
表 7 中国情境下典型植物工厂案例与模型路径对应关系Table 7 Typical PFAL cases in China and their correspondence with model pathways |
| 案例 | 公开资料显示的主要特征 | 对应模型路径 | 对本研究结论的启示 |
|---|---|---|---|
| 京东植物工厂 | 自有品牌“京觅”;线上京东生鲜与线下七鲜销售;冷链覆盖区域市场;年产量约300 t;强调安全、标准化和渠道能力 | 价值升级型;综合转型型 | 高值直供市场需要品牌、冷链、销售渠道和终端信任支撑;高价格情景不能只看售价,也要计入渠道成本 |
| 中科三安植物工厂 | LED光谱、环境控制和智能装备;叶菜、药用植物、茄果类等多品类布局;规模化产业基地;生产周期缩短 | 效率优先型;能源耦合型 | 控制升级和规模化可提高单位面积产出,但同时带来资本强度和能源依赖,需与能源条件和高附加值作物组合协同 |
4 讨 论
4.1 PFAL经济可行性是典型的“窄窗口”问题
4.2 控制成熟度同时改变收入项、成本项与单位成本的分母
4.3 脆弱性评价能够补充利润视角,但应避免过度解读
4.4 技术发展趋势对敏感性结构的影响
4.5 参数真实性、适用边界与局限性
5 结 论
附录
附表1 关键映射系数±10%扰动结果Appendix Table 1 Robustness under ±10% perturbations of key mapping coefficients |
| 扰动对象 | 扰动幅度 | 正利润情景数 | 最佳EVI情景 | 最差EVI情景 | 最高利润情景 |
|---|---|---|---|---|---|
| -10% | 5 | C2-E1-M3 | C0-E3-M1 | C2-E1-M3 | |
| 10% | 5 | C2-E1-M3 | C0-E3-M1 | C2-E1-M3 | |
| -10% | 5 | C2-E1-M3 | C0-E3-M1 | C2-E1-M3 | |
| 10% | 5 | C2-E1-M3 | C0-E3-M1 | C2-E1-M3 | |
| -10% | 5 | C2-E1-M3 | C0-E3-M1 | C2-E1-M3 | |
| 10% | 5 | C2-E1-M3 | C0-E3-M1 | C2-E1-M3 | |
| -10% | 5 | C2-E1-M3 | C0-E3-M1 | C2-E1-M3 | |
| 10% | 5 | C2-E1-M3 | C0-E3-M1 | C2-E1-M3 | |
| -10% | 5 | C2-E1-M3 | C0-E3-M1 | C2-E1-M3 | |
| 10% | 5 | C2-E1-M3 | C0-E3-M1 | C2-E1-M3 |
附表2 27个情景组合下的主要计算结果Appendix Table 2 Main computational results under 27 scenario combinations |
| 情景 | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|
| C0-E1-M1 | 70.4 | 12.0 | 25.27 | -1 356.60 | -1.06 | 1.60 | 0.763 | 0.261 | 0.068 | 0.364 |
| C0-E1-M2 | 70.4 | 12.0 | 25.27 | -723.00 | -0.31 | 1.60 | 0.406 | 0.261 | 0.068 | 0.245 |
| C0-E1-M3 | 70.4 | 12.0 | 25.27 | -89.40 | 0.44 | 1.60 | 0.050 | 0.261 | 0.068 | 0.127 |
| C0-E2-M1 | 70.4 | 12.0 | 27.67 | -1 525.56 | -1.06 | 6.37 | 0.783 | 0.325 | 0.746 | 0.618 |
| C0-E2-M2 | 70.4 | 12.0 | 27.67 | -891.96 | -0.31 | 6.37 | 0.458 | 0.325 | 0.746 | 0.510 |
| C0-E2-M3 | 70.4 | 12.0 | 27.67 | -258.36 | 0.44 | 6.37 | 0.133 | 0.325 | 0.746 | 0.401 |
| C0-E3-M1 | 70.4 | 12.0 | 30.07 | -1 694.52 | -1.06 | 8.16 | 0.800 | 0.379 | 1.000 | 0.727 |
| C0-E3-M2 | 70.4 | 12.0 | 30.07 | -1 060.92 | -0.31 | 8.16 | 0.501 | 0.379 | 1.000 | 0.627 |
| C0-E3-M3 | 70.4 | 12.0 | 30.07 | -427.32 | 0.44 | 8.16 | 0.202 | 0.379 | 1.000 | 0.527 |
| C1-E1-M1 | 89.21 | 10.2 | 20.9 | -1 329.00 | -0.91 | 1.36 | 0.713 | 0.268 | 0.034 | 0.338 |
| C1-E1-M2 | 89.21 | 10.2 | 20.9 | -526.09 | -0.03 | 1.36 | 0.282 | 0.268 | 0.034 | 0.195 |
| C1-E1-M3 | 89.21 | 10.2 | 20.9 | 276.82 | 0.85 | 1.36 | 0.00 | 0.268 | 0.034 | 0.101 |
| C1-E2-M1 | 89.21 | 10.2 | 22.94 | -1 511.00 | -0.91 | 5.41 | 0.738 | 0.334 | 0.610 | 0.561 |
| C1-E2-M2 | 89.21 | 10.2 | 22.94 | -708.08 | -0.03 | 5.41 | 0.346 | 0.334 | 0.610 | 0.43 |
| C1-E2-M3 | 89.21 | 10.2 | 22.94 | 94.83 | 0.85 | 5.41 | 0.000 | 0.334 | 0.610 | 0.315 |
| C1-E3-M1 | 89.21 | 10.2 | 24.98 | -1 692.99 | -0.91 | 6.93 | 0.760 | 0.388 | 0.826 | 0.658 |
| C1-E3-M2 | 89.21 | 10.2 | 24.98 | -890.08 | -0.03 | 6.93 | 0.399 | 0.388 | 0.826 | 0.538 |
| C1-E3-M3 | 89.21 | 10.2 | 24.98 | -87.16 | 0.85 | 6.93 | 0.039 | 0.388 | 0.826 | 0.418 |
| C2-E1-M1 | 109.25 | 8.4 | 17.58 | -1 264.92 | -0.83 | 1.12 | 0.659 | 0.263 | 0.000 | 0.307 |
| C2-E1-M2 | 109.25 | 8.4 | 17.58 | -281.67 | 0.24 | 1.12 | 0.147 | 0.263 | 0.000 | 0.136 |
| C2-E1-M3 | 109.25 | 8.4 | 17.58 | 701.58 | 1.31 | 1.12 | 0.000 | 0.263 | 0.000 | 0.088 |
| C2-E2-M1 | 109.25 | 8.4 | 19.26 | -1 448.46 | -0.83 | 4.46 | 0.688 | 0.327 | 0.474 | 0.497 |
| C2-E2-M2 | 109.25 | 8.4 | 19.26 | -465.21 | 0.24 | 4.46 | 0.221 | 0.327 | 0.474 | 0.341 |
| C2-E2-M3 | 109.25 | 8.4 | 19.26 | 518.04 | 1.31 | 4.46 | 0.000 | 0.327 | 0.474 | 0.267 |
| C2-E3-M1 | 109.25 | 8.4 | 20.94 | -1 632.00 | -0.83 | 5.71 | 0.713 | 0.381 | 0.652 | 0.582 |
| C2-E3-M2 | 109.25 | 8.4 | 20.94 | -648.75 | 0.24 | 5.71 | 0.284 | 0.381 | 0.652 | 0.439 |
| C2-E3-M3 | 109.25 | 8.4 | 20.94 | 334.50 | 1.31 | 5.71 | 0.000 | 0.381 | 0.652 | 0.344 |
附表3 EVI权重稳健性检验结果Appendix Table 3 Robustness test of EVI weighting schemes |
| 权重方案 | 权重设置( , , ) | 最佳情景 | 最差情景 | 前5重合数 | 后5重合数 | Spearman | Kendall |
|---|---|---|---|---|---|---|---|
| 等权重 | 0.333, 0.333, 0.333 | C2-E1-M3 | C0-E3-M1 | 5 | 5 | 1.000 | 1.000 |
| 盈利缺口加权 | 0.500, 0.250, 0.250 | C2-E1-M3 | C0-E3-M1 | 4 | 4 | 0.944 | 0.829 |
| 能源暴露加权 | 0.250, 0.500, 0.250 | C2-E1-M3 | C0-E3-M1 | 5 | 5 | 0.999 | 0.989 |
| 碳约束加权 | 0.250, 0.250, 0.500 | C2-E1-M3 | C0-E3-M1 | 5 | 4 | 0.976 | 0.886 |
| 剔除 | 0.000, 0.500, 0.500 | C2-E1-M1/M2/M3并列 | C0-E3-M1 | 3 | 3 | 0.842 | 0.682 |
| 剔除 | 0.500, 0.000, 0.500 | C2-E1-M3 | C0-E3-M1 | 5 | 5 | 0.991 | 0.943 |
| 剔除 | 0.500, 0.500, 0.000 | C2-E1-M3 | C0-E3-M1 | 3 | 4 | 0.782 | 0.647 |





