矿产资源的管理模型 《数学模型》(三版)电子姜启源

矿产资源的管理模型

基于优化商业利润与生产可持续的角度来研究主要讨论以下几个模型：1、开发会枯竭资源的泛函最大值模型2、矿业的非线性模型3、有采掘成本的最优控制模型重点确定最优开发速度以及税金等所建立的数学模型主要用到微分方程和泛函变分的最大值原理其建模思想和原理对其它可枯竭或可再生资源亦适用。1、开发会枯竭资源的泛函最大值模型令为t时刻的贮量，初始贮量为是已知的，表示采掘率，于是矿物贮量的数学模型是线性微分方程的哥西问题:其中(1)(2)设矿石价格为是已知的，经营者的目标是使收益现值最大，设贴现率为，则时刻t的收益满足其中即现值，解得，故值，对于上述开矿问题，其总现值PV为(3)我们需要确定

，以使总现值最大。这样把问题转化成了求泛函的最大值问题。考虑哥西问题（4）（5）其中是已知的三元连续函数，称为状态变量，（4）、（5）为状态方程.对于每个给定的函数，（4）、（5）的解称为的一个响应，这里是区间上的分段连续函数集合.其中是已知连续可微函数，最优控制理论中的核心问题是求使最大的所谓可行控制，如果这样的可行控制存在，则称其为最优控制。设，且设终点是固定的，而起点是可变的，设是最优控制，即函数称为最优获利函数。(6)（7）引入

的泛函在整个区间的最优控制策略对划分的各个子区间也分别是最优的对（8）式进行计算得即由（7）可得其中是与最优控制对应的最优轨线，公式（8）指出下面的最优性原理：（8）上式两端除以且令取极限得即考虑任意函数，定义控制

为其中是对应于起点的最优控制，是对

的响应，于是用代替变点得令固定,由(11)得(10)则有(11)(12)由(11),当时,达到它的最大值,在足够光滑的条件下令(13)(14)故由(12)得(15)其中右端的自变量皆为,由的定义称为哈密顿函数,为伴随变量,上述推导得到下面的重要结论于是有(16)

最大值原理：对于满足边条件的状态方程使目标泛函取最大值的最优控制和相应的最优响应满足伴随方程：(16)和(15)要确定的最优控制和最优响应可由及伴随方程和状态方程解出。例如对于矿产资源管理工作线性模型(1)(2)(3)，其哈密顿函数是由伴随方程得(17)所以，假设到时矿产采完，则,因此为使哈密顿函数达到最大，则得（18）式给出了为使采矿收益的现值最大采用的最优采掘率。即当矿产价格超过时，用采掘率去开采，否则停产，出现若干生产时期，则满足，其中，见图2.1(18)0tq图2.1元$0t2、矿业的非线性模型若价格受采掘率的影响，,则与（3）对应的目标泛函为以及条件哈密顿函数是设,对（22）又，故得（23）设价格随采掘率之提高而线性下降由最大值原理由伴随方程式，所以相应于最优控制采掘率，价格则由（20）得由最大值原理，为使取最大值，即T由（24）解出，由于，故当时单调上升，所以满足（24）的T，即采完日期是唯一存在的。采掘日起T满足（24）3、有采掘成本的最优控制模型最大设采掘成本，且是的增函数，即采掘越快，成本越高，这时要求最优控制，使得（25）由伴随方程得相应的哈密顿函数为对于任固定的，在为变量的坐标系上（见图），是（26）的唯一正数解，如果曲线与相交，这种正数解存在，注意这里的已固定，视为常数了。由最大值原理，，即(26)唯一的正数解是的斜率为点的横坐标，如果与只一个公共点，则如果多家竞争开发一种矿藏，设有家厂商组成的采矿业，第家的成本为是第家厂商的采掘率，则第家厂商的目标泛函为（27）于是有（28）是第厂商的开采上限；又由市场平衡知（30）(28)(29)(30)中共有2n+1个未知量有2n+1个方程，可从中解得各自的最优控制以及（29）设成本为是常数，时表示价格补贴率，时表示税率，则目标函数为（31）由最大值原理可以导出（32）最后讨论有价格补贴或税金的情形如果成本则由（32）得又由于到时刻T开采完毕，由得（33）（34）采完日期可由公式确定，如果价格则