离散资源依赖型开采模型中状态依赖承载力的分岔动力学与混沌控制

《Journal of Biological Dynamics》：Bifurcation dynamics and chaos control in a discrete stock-dependent resource-extraction model

　　

专业术语翻译摘要 本研究考察了一个具有状态依赖承载力的离散非线性资源-努力系统，其中最大开采能力自适应于资源可用性。与经典固定承载力模型不同，所提出的框架引入了资源与开采变量之间的自适应反馈，为分析非线性离散系统中的复杂动力学提供了一种广义方法。利用分段常数参数法（piecewise constant argument），将连续模型转化为非线性离散映射。分析了包括正性、有界性、平衡点及局部稳定性在内的基本性质。通过flip分岔（倍周期分岔）和Neimark–Sacker分岔探讨了内平衡点处稳定性的丧失，进而导致倍周期级联、准周期振荡与混沌。应用正规形理论（normal form theory）确定分岔方向与稳定性。利用分岔图、相图、时间序列及最大Lyapunov指数（maximum Lyapunov exponents）进行的数值模拟证实了复杂的转变过程。进一步提出了一种反馈控制策略，以稳定不稳定动力学并扩大稳定区域。

论文解读文章

一、研究背景与问题

自然资源开采的数学建模长期以来受到广泛关注，其核心在于资源可用性、经济决策、技术发展与环境约束之间的复杂交互。与可再生资源不同，矿产资源具有有限性和不可再生性，开采后面临品位下降、开采成本上升、开采深度增加等一系列问题，这使得采矿系统呈现出独特的动力学特征。经典资源经济学模型以Hotelling模型为代表，其假设开采努力为外部给定的控制变量，并基于理性预期和均衡假设进行优化分析。然而，实际开采过程中，开采努力水平并非外部可控变量，而是依赖于剩余资源存量的动态变量——高品位矿床的枯竭会导致矿石品位降低、能耗增加、开采成本上升，从而制约未来的开采能力。

在经典stock-effort模型（如Gordon-Schaefer模型）中，开采努力通常被视为与资源存量耦合的状态变量，但其承载力被假设为常数，忽略了资源存量对开采能力的反馈作用。现有采矿模型大多集中于优化框架，将开采水平视为可控因素而非动态变量，且多数模型假设承载力恒定。数学上，针对开采努力水平变化的离散采矿模型的分岔与混沌动力学研究仍十分有限。因此，本研究旨在构建一个将状态依赖的开采努力承载力引入离散资源开采框架的模型，系统分析其稳定性、分岔结构、复杂动力学及混沌控制特性。

二、研究内容与结论

研究人员提出了一种二维stock-effort模型，其中资源存量x(t)与开采努力y(t)作为相互作用的动态变量。模型的关键创新在于开采努力的承载力并非恒定，而是依赖于资源存量，即y_max = mx。这一机制建立了资源可用性与开采能力之间的双向反馈：资源枯竭降低可持续开采能力，而开采活动进一步减少资源存量。利用分段常数参数法，研究人员将连续模型转化为离散非线性映射，并系统分析了其动力学性质。

研究得出以下主要结论：系统存在两个不动点，即边界平衡点(K,0)和内平衡点(Kr/(r+αKm), Kmr/(r+αKm))；正象限是系统的不变区域；通过特征方程分析，确定了内平衡点失稳的条件，分别对应flip分岔（倍周期分岔）和Neimark–Sacker分岔的发生条件；正规形理论确定了分岔方向与稳定性；数值模拟证实了系统从稳定平衡到倍周期级联、准周期振荡直至混沌的复杂转变过程；提出的反馈控制策略能够有效稳定不稳定动力学并扩大稳定区域。

三、主要技术方法

本研究采用的主要技术方法包括：分段常数参数法（piecewise constant argument, PCA）用于将连续模型离散化，该方法不同于传统数值离散化，能够保持初始非线性系统的特征并允许在原始框架内分析不动点稳定性及各类分岔；局部分岔理论用于分析flip分岔和Neimark–Sacker分岔的条件与方向；正规形理论（normal form theory）用于确定分岔的方向与稳定性；数值模拟方法包括分岔图、相图、时间序列和最大Lyapunov指数（maximum Lyapunov exponents）分析；反馈控制策略用于稳定混沌和不稳定动力学。数值模拟中使用的参数为代表性无量纲值，而非针对特定矿区或经验数据校准的参数。

四、研究结果

4.1 模型建立与基本性质

研究人员建立了连续时间采矿资源开采模型，其中资源存量遵循logistic增长并受到开采活动的消耗，开采努力的演化受状态依赖承载力约束。通过分段常数参数法得到离散系统。证明了系统正性保持性质，即正象限为不变区域，确保资源存量与开采努力始终非负，这一性质对模型的经济学解释至关重要。

4.2 不动点与稳定性分析

系统存在两个不动点：边界平衡点(K,0)对应无开采活动的情形，内平衡点对应资源与开采共存的稳态。通过Jacobian矩阵特征值分析，研究人员给出了边界平衡点的稳定性分类：当h > 2/r时为源（source），当h < 2/r时为鞍点（saddle），当h = 2/r时为非双曲（nonhyperbolic）。对内平衡点，通过特征方程F(λ) = λ2 + pλ + q，利用Jury判据给出了稳定性条件，并确定了flip分岔与Neimark–Sacker分岔发生的参数条件集合。

4.3 Flip分岔（倍周期分岔）

当Jacobian矩阵的一个特征值穿过单位圆达到?1时，系统发生flip分岔。研究人员以离散化步长h为分岔参数，通过小扰动分析得到了扰动系统。利用正规形理论确定了分岔的方向与稳定性。数值模拟显示，随着h增大，系统从稳定不动点进入周期二振荡，随后经历倍周期级联，最终进入混沌状态。分岔图清晰地展示了这一转变过程。

4.4 Neimark–Sacker分岔

当Jacobian矩阵的两个复共轭特征值的模同时达到1时，系统发生Neimark–Sacker分岔，对应离散系统中的Hopf分岔。此时稳定不动点失稳，系统出现不变闭曲线上的准周期运动。数值模拟中的相图显示了不变闭曲线的形成，时间序列呈现准周期振荡特征，最大Lyapunov指数在混沌区域变为正值，证实了混沌动力学的存在。

4.5 混沌控制

针对系统中出现的混沌和不稳定动力学，研究人员提出了反馈控制策略。该策略通过引入控制项来稳定不稳定的周期轨道和混沌行为，从而扩大系统的稳定区域。数值模拟验证了控制策略的有效性，表明适当的反馈控制能够将系统从混沌状态引导回稳定或周期状态，为资源管理中的决策调整提供了理论依据。

五、讨论与结论

研究表明，离散化步长h作为采矿决策或调整的时间间隔，对系统稳定性具有显著影响——调整频率过高或过低均可能导致系统失稳。状态依赖承载力机制使该模型区别于经典固定承载力模型，能够产生更丰富的动力学行为。参数r、K、α、m、s的变化均可能触发分岔，导致系统在稳定、周期、准周期和混沌状态之间转变。这些复杂动力学行为对应于不稳定的生产动态、资源过度开采或意外的经济波动，因此对混沌行为的控制具有重要的实际意义。研究结论指出，该模型为理解资源可用性、开采效率、采矿调整与状态依赖开采能力如何交互塑造资源管理动力学提供了数学基础。未来研究需要对特定矿区进行经验参数校准，并进一步探索控制策略在实际资源管理中的应用。