Managing Sustainable Tourism（优秀范文一：Butler 生命周期 + 承载力退化反馈动力系统）

2025 美赛 MCM B 题 · 优秀范文一。视角：把「网红旅游地为何兴起—巅峰—衰退」讲透——到访量 $V(t)$ 受逻辑斯蒂增长驱动、被口碑下滑拖拽，而环境承载力并非常数：累积退化 $D(t)$ 以指数方式侵蚀有效承载力 $C_{\mathrm{eff}}(D)=C_0e^{-\kappa D}$ ，形成「客流→退化→承载力缩水→客流回落」的闭环反馈。用三方程 ODE 刻画 Butler 生命周期的动力学本质，识别真超载窗口（ $V>C_{\mathrm{eff}}$ ）与社会不满越限窗口，并以足迹强度 $\alpha$ 与生态韧性 $\beta$ 两组扫描给出命运分岔图。正文 / 配图 / 附录 / 真源四路数字一致（确定性 ODE，无随机漂移）。

摘要

旅游地兴衰的本质是正负反馈的赛跑。本文把目的地抽象为两个状态变量——年到访量 $V(t)$ （万人次）与累积环境退化 $D(t)$ ，建立三方程动力系统： $\dot V=rV(1-V/C_{\mathrm{eff}})-\delta D$ （逻辑斯蒂增长减口碑下滑的绝对损失）， $\dot D=\alpha V-\beta D$ （客流足迹生成退化、生态按恢复率 $\beta$ 自愈）， $C_{\mathrm{eff}}=C_0e^{-\kappa D}$ （承载力随退化指数缩水）。取 $r=0.22/$ 年、 $C_0=250$ 万人次、 $\kappa=0.60$ 、 $\delta=1.5$ 、 $\alpha=0.0012$ 、 $\beta=0.10$ 的标称参数，欧拉法（ $\Delta t=0.01$ 年、仿真 120 年）解得经典 Butler 轨迹：第 24.0 年冲上峰值 131.9 万人次，随后回落 19.2%，收敛到长期稳态 106.6 万人次；退化峰值 $D_{\max}=1.30$ 使有效承载力最低压到 114.5 万人次，第 26–32 年出现长达 6.25 年的「真超载」窗口（ $V>C_{\mathrm{eff}}$ ，利用率峰值 1.019）；居民不满 $U(V)=0.4(V/100)^2$ 在 $V>122.5$ 万人次时越过社会阈值 0.6，累计越限 9.50 年（占 7.9%）——社会约束比生态超载更早、更久。参数扫描揭示两条管理铁律：足迹强度 $\alpha$ 翻 4 倍（0.0006→0.0024），长期稳态从 143.2 万腰斩至 73.0 万（−49%）；生态恢复率 $\beta$ 从 0.06 提到 0.14，稳态从 81.3 万升至 124.2 万（+53%）。可持续经营的目标不是追逐峰值，而是把 $V$ 压在 $0.9\times C_{\mathrm{eff}}$ 预警线之下的可持续走廊内。

一、问题重述

题目要求回答四件事：①刻画目的地旅游地生命周期——到访量如何随口碑、承载力与退化反馈演化；②量化旅游的环境与社会压力，识别承载力阈值与超载时段；③预测不同情境下客流与影响的演化路径；④设计管理政策使系统长期可持续。本文聚焦①—③构建统一动力系统并完成情景扫描，④的政策优化留给范文二、范文三分别从多目标规划与蒙特卡洛鲁棒性视角展开。核心是把三句定性常识翻译成方程：「游客越多退化越快」「退化越重能承载的游客越少」「口碑崩了客流就跌」。

二、模型假设

目的地同质单一景区，年到访量 $V$ 为唯一客流状态量；初始 $V_0=5$ 万人次（新兴期）、 $D_0=0.05$ 。
客流内在增长率 $r=0.22/$ 年恒定（营销投入、交通改善等外生因素不显式建模）；无退化时环境承载力 $C_0=250$ 万人次/年。
每单位到访量每年生成 $\alpha=0.0012$ 单位退化（足迹强度，涵盖垃圾、植被践踏、水质下降的综合折算）；生态系统以速率 $\beta=0.10/$ 年自然恢复。
口碑下滑对客流的拖拽正比于当期退化水平 $D$ ，系数 $\delta=1.5$ 万人次/年·单位退化——退化通过差评、负面舆情转化为绝对客流损失。
居民不满只依赖拥挤度： $U=u_0(V/V_{\mathrm{ref}})^{\gamma}$ ， $u_0=0.4$ 、 $V_{\mathrm{ref}}=100$ 、 $\gamma=2$ （凸性：越挤越敏感），社会阈值 $U>0.6$ 即视为社区抵触显性化。

三、符号说明

符号	含义
$V(t)$	年到访量（万人次）
$D(t)$	累积环境退化指数（无量纲）
$C_{\mathrm{eff}}(D)$	有效承载力 $=C_0e^{-\kappa D}$ （万人次/年）
$r,\ \delta$	客流内在增长率（/年）、口碑下滑系数
$\alpha,\ \beta$	足迹强度（单位退化/万人次）、生态恢复率（/年）
$\kappa$	承载力退化弹性
$u(t)=V/C_{\mathrm{eff}}$	承载力利用率； $u>1$ 为真超载
$U(V)$	居民不满指数

四、建模过程

4.1 从 Butler 叙事到状态方程

Butler 生命周期曲线把旅游地分为探索—参与—发展—巩固—停滞—衰退/复苏六阶段，但那是描述性的。给它装上动力学引擎只需三个方程：

$\dot V=rV\Big(1-\frac{V}{C_{\mathrm{eff}}(D)}\Big)-\delta D,\qquad \dot D=\alpha V-\beta D,\qquad C_{\mathrm{eff}}=C_0e^{-\kappa D}.$

第一式是「带时变天花板」的逻辑斯蒂增长再减一项口碑拖拽——注意拖拽项 $\delta D$ 是绝对量而非比例项，这正是峰值能够触及甚至略微越过承载力线的数学前提；第二式中退化由客流驱动、又按指数恢复，天然形成滞后响应（时间尺度 $1/\beta=10$ 年，远慢于客流的 $1/r\approx4.5$ 年），滞后正是过冲的来源；第三式让天花板本身随退化下沉。图 1 给出标称参数下的完整生命周期：前 15 年近乎指数兴起，第 24.0 年登顶 131.9 万人次，随后在退化反馈下拉回落 19.2%，约 60 年后稳定在 106.6 万人次——兴起、巅峰、回落、稳态四段俱全。

相空间视角更能看清机制。图 2 的 $(D,V)$ 轨迹从左下角出发先向右上冲刺（客流先行、退化跟涨），过峰后沿退化轴小幅回撤，最终收敛到稳态点 $(1.28,\ 106.6)$ （图中红点）。轨迹不绕圈、不振荡，说明该参数域内系统是单调收敛的——衰退以「一次深呼吸」而非反复痉挛的方式发生。

4.2 承载力退化反馈与真超载窗口

把 $V(t)$ 与 $C_{\mathrm{eff}}(t)$ 叠在一起（图 3），故事的关键角色登场：客流冲峰时退化尚未到位， $C_{\mathrm{eff}}$ 还维持在高位；峰值过后退化继续爬升到 $D_{\max}=1.30$ ，把有效承载力压到全程最低的 114.5 万人次——而客流此时恰好也回到这个量级，两条线在第 26–32 年交叉穿插，形成 6.25 年的真超载窗口（25 个采样点 $V>C_{\mathrm{eff}}$ ，占仿真时长 5.2%，利用率峰值 1.019）。这回答了「超载发生在何时」：不在客流最高的巅峰时刻，而在峰值之后、退化最深的回落段——最危险的时刻恰是看起来「形势已经好转」的时刻，因为承载力比客流掉得更快。

4.3 命运分岔：足迹强度 α 决定长期命运

同一目的地，游客行为模式不同，命运可以完全不同。图 4 把足迹强度从 $\alpha=0.0006$ （轻徒步型）扫到 $0.0024$ （大巴打卡型）：四条曲线起点相同、早期几乎重叠，却在 20 年后分道扬镳——低足迹曲线稳稳站上高平台，高足迹曲线冲峰后深度回落。图 5 汇总成命脉图： $\alpha$ 翻 4 倍，峰值仅从 163.6 万降到 99.3 万（−39%），长期稳态却从 143.2 万腰斩到 73.0 万（−49%），且回落幅度从 12% 扩大到 26%。直觉解释：峰值由增长期的「惯性」决定，对 $\alpha$ 不太敏感；稳态却由「退化—承载力」的长期平衡决定，对 $\alpha$ 高度敏感。短期看客流、长期看足迹——这是本模型最有政策含金量的结论。

4.4 生态韧性与社会阈值：两道隐形的闸门

图 6 回答「如果生态更脆弱/更有韧性会怎样」：恢复率 $\beta$ 从 0.06 提到 0.14，长期稳态从 81.3 万线性升至 124.2 万（+53%）。 $\beta$ 是管理者最容易忽视的杠杆——它不由门票决定，而取决于生态修复投入、轮休制与保护区的实质性维护。每提升 0.01 的恢复率约换回 5.4 万人次的可持续客流，「养生态就是养客流」在此有了定量表达。

社会维度同样是一道闸门。图 7 画出居民不满 $U(V)=0.4(V/100)^2$ 的凸性上升曲线与阈值线 0.6 的交点： $V=122.5$ 万人次。标称轨迹上有 38 个采样点越限，累计 9.50 年（占 7.9%）——比真超载窗口（6.25 年）更长且更早开始。含义直白：在本文参数下，社区居民的忍耐比生态环境先绷不住；若无视社会阈值，抵触情绪（本地反对新增酒店、抵制游客）将转化为额外口碑损失，相当于有效 $\delta$ 被动放大，进一步压低稳态。

五、可持续走廊：管理目标的图形化定义

把三张图合成一张管理地图（图 8）：红线是动态承载力 $C_{\mathrm{eff}}$ ，橙线是留有 10% 安全余量的预警线 $0.9\times C_{\mathrm{eff}}$ 。标称轨迹在冲峰段整体贴线运行并短暂越线——这正是「不加管理的自由生长」的典型形态。可持续经营的定义由此图形化：通过限流、预约、动态定价与淡旺季调节，让 $V$ 长期位于橙线之下。贴红线经营看似多做了一成一成的生意，代价却是超载段的体验滑坡与退化加速——模型显示这段透支要用之后数十年的低稳态来偿还。

六、灵敏度分析

增长率 $r\pm20\%$ ：主要改变冲峰速度与峰值时点（第 21–27 年间移动），对长期稳态影响不足 ±3%——什么时候到顶不重要，顶有多高之后的落点才重要。
承载力弹性 $\kappa\pm0.1$ ： $\kappa=0.7$ 时稳态降至约 97 万、超载窗口拉长； $\kappa=0.5$ 时稳态升至约 115 万。弹性刻画「环境多娇气」，宜用历史水质/植被数据实测校准。
口碑系数 $\delta$ ：拖拽翻倍（ $\delta=3$ ）使峰值降约 9%、稳态基本不变——口碑影响的是「跌得有多疼」，不是「最后剩多少」。
初值 $V_0, D_0$ ：稳态与初值无关（单调收敛系统），只有到达时间受影响；这说明稳态是结构参数（ $\alpha,\beta,\kappa,\delta,r,C_0$ ）的属性，管理必须作用于结构，不能指望「换个起点」。

七、模型评价与改进

优点：三个方程各自对应一条可实证的证据链（客流序列拟合 $r$ 、环境监测反演 $\alpha,\beta$ 、舆情调查标定 $U$ ），参数少而机理完整；「真超载窗口」与「社会越限窗口」两个指标可直接翻译为管理动作； $\alpha$ / $\beta$ 双杠杆扫描给出了清晰的优先级排序。局限：①单景区同质假设忽略了时空分流的空间维度——热点瞬时超载比年均超载更致命；②口碑拖拽取线性 $\delta D$ ，真实口碑存在临界与雪崩效应；③未引入价格与政策变量，限流/定价只能事后解读而不能内生优化。改进方向：范文二将把限流比例、定价与分流作为控制变量做多目标优化；范文三将以蒙特卡洛注入需求冲击，检验政策鲁棒性。

八、结论

旅游地的兴衰不是玄学，而是三方程动力系统的确定演化：客流以 $r=0.22$ 的内在速度冲向一个正在下沉的天花板，第 24.0 年触顶 131.9 万人次后回落 19.2%，稳态锁定 106.6 万；真正的危险不在巅峰而在其后——第 26–32 年承载力被退化压到最低、客流反向越线的 6.25 年真超载窗口，叠加更早更久的 9.50 年社会不满越限。命运的分岔由两个「慢变量」执掌：足迹强度翻 4 倍、稳态腰斩；恢复率提高四成、稳态加五成。因此可持续经营的第一性原理是：少谈冲峰、多管足迹、养护韧性，把客流永久关进 $0.9\times C_{\mathrm{eff}}$ 的走廊里。

参考文献

COMAP. MCM Problem B: Managing Sustainable Tourism. 2025.
Butler R. W. The Concept of a Tourist Area Cycle of Evolution. Canadian Geographer, 1980.
Lodhi S., Mikic A., Ivanov I. Sustainable Tourism Development Model. 2016.
Hardin G. The Tragedy of the Commons. Science, 1968.
Dredge D., Jenkins J. Tourism Planning and Policy. Wiley, 2007.

附录：核心 Python 实现

import sys, os
_HERE = os.path.dirname(os.path.abspath(__file__))
sys.path.insert(0, os.path.abspath(os.path.join(_HERE, "..", "..", "..", "tools")))
import gen_mcm2025b as M   # 真源：Butler 生命周期 + 承载力退化反馈

# 1) 标称轨迹：峰值 / 稳态 / 回落幅度
ts, Vs, Ds, Ces = M.simulate()
pt, pv = M.peak_info(ts, Vs)
ss = M.steady_state(ts, Vs)
print("峰值: 第%.1f年 V=%.1f万人次  长期稳态=%.2f万人次  回落%.1f%%"
      % (pt, pv, ss, 100.0 * (pv - ss) / pv))

# 2) 真超载窗口与利用率峰值
ovl_n, ovl_pct, t0, t1 = M.overload_stats(ts, Vs, Ces)
print("真超载: %d个采样点=%.2f年(占%.1f%%) 第%.2f-%.2f年  利用率峰值=%.3f"
      % (ovl_n, ovl_n * 0.25, ovl_pct, t0, t1, M.utilization_peak(Vs, Ces)))

# 3) 社会阈值越限
v_thr, sb_n, sb_pct = M.social_stats(ts, Vs)
print("社会阈值=%.1f万人次  不满越限=%.2f年(占%.1f%%)"
      % (v_thr, sb_n * 0.25, sb_pct))

# 4) 两大杠杆：足迹强度 alpha 与生态恢复率 beta
fr = M.fate_vs_alpha()
print("足迹分岔: a=0.0006稳%.1f -> a=0.0024稳%.1f (降幅%.0f%%)"
      % (fr[0][2], fr[-1][2], 100.0 * (fr[0][2] - fr[-1][2]) / fr[0][2]))
bs = M.beta_scan()
print("韧性扫描: b=0.06稳%.1f -> b=0.14稳%.1f (+%.0f%%)"
      % (bs[0][1], bs[-1][1], 100.0 * (bs[-1][1] - bs[0][1]) / bs[0][1]))