天然气产量预测.docVIP

天然气产量预测模型的研究 摘 要 针对于本文天然气产量的预测问题，考虑到要求出前24个月的总产量，从而得到一个累加模型。对于求解该模型可以采用Matlab软件进行累加实现，由于定积分在一定程度上可以将问题从不连续转化为连续，大大减少了计算的步骤，使模型的求解变得更加的容易，所以本文也采用定积分的计算方法，对该模型进行了求解。在比较两种方法后，分析两者相对误差仅为0.013%。 将得到的两个结果取平均值，同时在参考假设的条件下，通过取平均值估计出了前24个月的总产量为：。 关键词 预测 累加 定积分 总产量 1 问题的提出 工程师们已经开始从墨西哥的一个新井开采天然气，根据初步的试验和以往的经验，他们预计天然气开采的第个月的月产量的函数为：（百万立方米），若能估计出前24个月的总产量，这对于发展经济及建设具有极其重要的作用。 2 问题的分析 针对于本问题，考虑到要求出前24个月的总产量，这就将问题转化为了每个月的产量的累加和，从而得到一个累加模型。对于求解该模型可以采用Matlab软件实现，同时也可以采用定积分或其它算法，对该模型进行求解。在比较几种方法后得到更好的结果，分析误差和进行模型的推广等等。 3 模型的假设 1 假设天然气的开采量在预计值的误差范围内； 2 假设天然气在24个月内不会开采完； 3 假设天然气在当地区的开采是连续性的。 4 符号说明 符 号 说 明 当地牌子的每听售价 外地牌子的每听售价 果汁的总利润 果汁的总收益 果汁的总成本 当地牌子果汁的总收益 外地牌子果汁的总收益 当地牌子果汁的总成本 外地牌子果汁的总成本 5 模型的建立和求解 5模型的建立与求解 5.1 模型的建立 针对于本问题，根据工程师的预测[1]，可以得到： ， （1） 其中：天然气开采的第个月的月产量的函数； 根据总产量为24个月产量和的原理，得： ， （2） 图 5-1 天然气产量每月的产量图 5.2 模型的求解 5.2.1 用Matlab[2]求解产量值： 根据， （3） 用Matlab软件（源代码见附录一）求解，从而求得。 5.2.2 用定积分求解： 令 （4） 则 （5） 且，从而得到为递增函数。 对于（4）式，在内进行定积分，并利用Matlab软件[3]（源代码见附录二），得到。 5.3 问题的结论： 通过比较用累加法求解和定积分求解得到结果，发现两结果相对误差为0.013%，这在误差范围内是允许的，于是本文采用平均值的方法求解。 本文在参考假设的条件下，通过取平均值估计出了前24个月的总产量为： 6模型的优缺点分析 本文仅从工程师的估计中得到天然气的开采量，这在实际情况下过于片面，同时天然气在地下会存在一定的不可探测性，而这在文章和现有资料中是很难考虑的，这也是本文的缺点所在。 在本文中存在的优点，主要体现在： 一、此Matlab程序求累加值在大部分公式中存在普遍性，对于同类问题可以在稍加修改后将模型[4]进行应用；二、定积分在一定程度上可以将问题从不连续转化为连续，大大减少了计算的步骤，使本文的求解变得更加的容易。 参考文献 [1]姜启源,谢金星,叶俊.数学模型.北京:高等教育出版社,2003. [2]王沫然,MATLAB 5.X 与科学计算,北京:清华大学出版社,2000. [3]陈桂明,戚红雨,潘伟编著,Matlab 数理统计（6.X）,北京:科学出版社,2002. [4]雷功炎,数学模型讲义,北京:北京大学出版社,1999. 附 录 Matlab累加求解源程序 p=0; for k=1:24 p=p+0.0849*k*exp(-0.02*k); end p Matlab函数积分源程序 clear clc syms k p=0.0849*k*(exp(-0.02*k)); int(p,k,0.5,24.5)