煤田开发规划的合理编制问题Word下载.docx
- 文档编号:5930242
- 上传时间:2023-05-05
- 格式:DOCX
- 页数:21
- 大小:163.19KB
煤田开发规划的合理编制问题Word下载.docx
《煤田开发规划的合理编制问题Word下载.docx》由会员分享,可在线阅读,更多相关《煤田开发规划的合理编制问题Word下载.docx(21页珍藏版)》请在冰点文库上搜索。
费煤比:
综合效益比:
5、模型的分析、建立与求解
问题一:
因此,本文在规划“十一五”期间煤井的增产措施之前必须对老井在规划期内的各指标值进行预测。
下面首先对老井在规划期内的产煤量、产水量、费用建立灰色预测模型来进行预测:
5.1灰色理论简介
因为产煤量与产水量是多种因素确定的,要对其准确定量地描述预测是极其困难的。
根据灰色系统理论,我们可以不去研究这复杂系统的内部因素及相互关系,而从“产煤量,产水量”这两灰色量本身去挖掘有用信息,利用它的动态记忆特性,建立灰色模型,并以此建立模型对未来状态作出预测。
灰色系统在预测时,必须采用一定的方式对原始数据进行生成处理,使生成数据序列变成有规则序列。
数据生成有两个目的:
为模型提供中间信息;
弱化原随机数列的随机性。
常用的数据生成方式有累加生成和累减生成,设原始序列为
,则
次累加生成的结果为:
,其中
;
灰色理论在累加生成的基础上,采用灰色微分方程模型得到生成模型,记为GM(n,h),n是微分方程的阶数,h是变量个数。
在GM(n,h)模型中,当h
2时,所建GM模型不能作预测用,只能用于分析因子之间的相互关系。
作预测用的GM模型一般为GM(n,1)模型,其中最重要的也是在实际中应用最多的是GM(1,1)模型。
下面为GM(1,1)模型原理:
其灰色微分方程为:
待定系数
和
分别称为发展灰数和内生控制灰数。
它们可以利用最小二乘法求解,其计算公式为:
,
在
的边界条件下,特解为
,
上式称为生成模型,同样当
时,也是1阶累加量的预测公式。
在此模型下,原始数据的预测公式为
,
5.2运用灰色预测对数据进行预测
5.2.1数据分析与预测的原则
通过对数据分析,发现一些井只有近五年的数据,新井在开采初期数据基数小,预测的稳定性和可信性较差,而且煤井产煤存在先增后减的趋势,如果在整个增减区间内使用灰色预测,将与灰色预测的累加指数规律存在偏差,造成精度的降低。
故使用近五年的产煤量来作为预测的已知数据。
“十五”煤井在有数据范围内处于产量增加状态,且产量和“九五”煤井初期完全吻合。
利用同一煤田的新井产煤量与开发年份的动态记忆特性,可将“九五”煤井2001~2005年产煤量的预测值作为“十五”“十五”井产水量由于已知数据少,且开采初期数据基数小,受外界因素波动较大(产水受外界因素影响更大,如雨水补充地下水),直接用其数据进行灰色预测发现相对误差均值大于0.2,不满足精度要求。
而理论上,“十五”产水应该与“九五”井有相关性,数据上“九五”井“十五”井前期产水数据也大致吻合。
考虑个井存在差异因素,根据数据吻合趋势关系,故采用“九五”井2002~2006产水量预测值作为“十五”井2006~2010产水量的预测值。
“九五”煤井与“十五”煤井在开采初期的产煤量对比
注:
绿色曲线为“十五”井在开采5年内产煤曲线
红色曲线为“九五”井在开采第二年后产煤曲线
“九五”煤井与“十五”煤井在开采初期的产水量对比
蓝色曲线为“十五”井在开采5年内产水量
红色曲线为“九五”井在开采第二年后产水量曲线
5.2.2数据预测
利用灰色预测算法,运用MATLAB程序进行数据预测结果如下:
各类井网的产煤量数据及其预测数据
2001
2002
2003
2004
2005
2006
2007
2008
2009
2010
‘七五’井
实际
500.6
442.4
428.6
370.1
343.1
预测
500.60
450.30
411.87
376.72
344.56
315.15
288.25
263.65
241.15
220.57
“八
五”
井
354.7
318.0
280.7
246.6
229.0
354.70
315.55
281.74
251.48
224.51
200.42
178.92
159.74
142.60
127.30
“九
579.8
547.5
527.0
492.3
437.0
579.80
556.13
517.49
481.53
448.07
416.93
387.96
361.00
335.92
312.58
“十
72.3
218.2
297.1
416.1
508.7
579.80
各类井网的产水量数据及其预测数据
“七五”井
3198.5
2976.3
2929.5
2561.6
2330.3
3198.5
3047.8
2801.9
2575.9
2368.1
2177.0
2001.4
1839.9
1691.5
1555.0
587.7
636.6
653.1
671.2
670.1
587.70
640.15
651.73
663.52
675.53
687.75
700.20
712.87
725.77
738.90
941.5
1050.2
1147.2
1319.7
1364.3
941.5
1059.5
1159.7
1269.4
1389.5
1521.0
1664.9
1822.4
1994.8
2183.5
84.8
297.8
354.4
612.5
964.3
5.3、优化问题的求解:
由于上述灰色预测模型已经成功的预测出了老井网(“七五”、“八五”、“九五”和“十五”井)在“十一五”期间各年产煤量和产水量,并且这些老井网在“十一五”期间仍将使用。
因此我们可以在求解这一优化问题中除去这些老井网的影响,只通过对“十一五”期间各年所采用7种增产措施的井口数量求解来完成优化任务。
下面以“十一五”期间各年年产煤量为例来进行优化:
单井产煤量矩阵新增井口数量产煤矩阵
其中
为第j年采用第i种增产措施所开凿的井口数量(j=1~5分别表示2006年~2010年;
i=1~7),一共35个变量,用MATLAB求解时
分别对应于变量
~
。
产煤量
可以通过“单井产煤量矩阵”和“新增井口数量矩阵
”相乘表示成含有变量
的一次35元多项式形式。
例如产煤矩阵次对角线:
表示7种措施在2006本年产煤量;
表示7种措施在2006本年产煤量(其中
为2006年采取的措施在2007年所产的煤量,因此二者之和为2007的产煤量,
为2007采取的措施在本年产煤量,因此二者之和为2007的产煤量);
以此类推,
为2010本年产煤量(2006~2009所采取的措施在2010年本年所产生的煤量和2010年采取措施产生的煤量之和)。
、
分别表示新增措施产煤量、产水量与费用矩阵。
的求解与
方式一致,在次不做详述。
在处理优化问题中:
我们首先采取3种方案,每种方案只是目标函数不同,约束方程均一致。
目标函数1(五年内产煤量最大):
目标函数2(五年内产水量最小):
目标函数3(五年内费用最小):
其中约束方程:
约束方程主要由三部分组成:
新增措施井口数量约束、产煤量约束、产水量约束与费用约束。
由于MATLAB不能直接处理约束条件,我们可先将此约束条件改为统一的约束方程,得到如下矩阵A(35行22列);
再次MATLAB不能直接处理最大值问题,所以应把产煤下限转化为上限只需将不等式两端同时乘以“-1”,改变不等式的方向即可:
矩阵b=[426;
390;
430;
374;
445;
534;
476;
-32。
6;
-125.9;
-207.8;
-267.4;
-349.7;
548.2;
918.5;
1556;
1994.3;
2975.9;
69719;
98623;
111129;
116018;
114307];
求解这一次35元优化方程,分别得到三种方案的最优解:
(见附录1)
方案1(五年内产煤量最大):
1
426
2
100
21
269
3
323
107
4
188
136
5
283
6
98
134
302
7
73
403
产煤量:
1175.8万吨产水量:
2690.5万吨费用509690万元
方案2(五年内产水量最小):
234
156
430
00
374
440
192
342
476
1005.3万吨产水量:
1922.7万吨费用469630万元
方案3(五年内费用最小):
425
390
68
306
240
534
417
33
25
1004.6万吨产水量:
2029.5万吨费用451500万元
很显然:
由于3种方案3均采用单一目标函数来限定,很可能在满足(煤、水、费)一方面最优的同时,其它两种情况可能会不理想。
因此在上诉求解基础之上,我们建立了较为合理的多目标规划方案,即建立目标函数,使它同时权衡上述3种方案为:
方案4的解相当于将上述三种方案综合,最终得到一个折中方法,而非只苛求某一情况,使煤、水和费达到临近最优。
多目标规划最终最优方案4的解为:
方案4(多目标规划):
141
196
53
305
117
145
331
1011.6万吨产水量:
1947.5万吨费用467862万元
为了评定不同方案的好坏,我们设立产煤达标指数、产水达标指数、费用达标指数、盈亏指数、以及费煤比、水煤比来衡量方案的各个方面,从而有效的比较不同方案的优劣,并给出量化的结果。
降低产水量对煤开采的降污控本有着重要的影响。
水煤比反映了单位产煤量下,产水量的多少。
比值越低越有利于煤的持续科学开采。
因此,在煤井规划中要尽量保持水煤比的平衡
反映单位产煤量下,费用的多少。
比值越低,说明产煤的经济效益越高。
在煤井规划中也必须注重该指标的变化趋势
鉴于水煤比和费煤比都仅仅反映了煤开采各自的一个方面。
对于煤开采的整体反映不够全面。
本文将两者结合,提出了综合效益比。
它反映的是产水量、及费用的增加或减少对产煤量的影响,值越大说明综合效益越高
表18:
两类方案比较
评判指标
调整前的方案
调整后的方案
产煤不达到指数
1.44%
1.51%
产水不达标指数
12.73%
9.85%
费用不达标指数
3.16%
1.07%
盈亏指数
-1.79%
1.54%
水煤比
2.97147
2.44058
费煤比
534.1749
447.2481
综合效益比
2139405.127
1225058.629
题目的要求是尽量满足参考指标,而从模型的评判指数来看,我们规划的指标已经达到了题目的要求。
六:
模型的评价
1.利用灰色系统理论预测未知数据,可以不去研究这复杂系统的内部因素及相互关系,而从“产煤量,产水量”的时间序列这个综合灰色量本身去挖掘有用信息,利用它的动态记忆特性,建立灰色模型,并以此建立模型对未来状态作出预测,较为精确。
2.优化模型采用的是多目标规划方法,这种方法实用准确,可以反映对不同重要程度的目标进行有效的优化;
七、模型的推广与改进
我们的模型还可以应用到其他同类型的问题求解上,如对矿产资源,天然气资源等非可再生资源的利用开发上,模型可移植性好。
改进的方向:
对于老井的预测可以运用其它的方法,因为灰色预测主要是针对数据具有指数变化规律的情况,因而并非最佳的选择,可以进一步采用神经网络、时间序列等方法来进行预测,可能可以得到更好的结果;
而对于规划方案的提出,我们采用的是多目标规划的方法,该方法操作简便、易于理解也有很好的效果,但若想更好的反映煤田开采过程的动态变化,则可以考虑使用动态规划的方法。
九:
参考文献
1.姜启源,谢金星,叶俊,数学模型(第三版),高等教育出版社,2003年2月;
2.傅立,灰色系统理论及其应用,科学技术文献出版社,1991年11月;
3.丁艳清张彩霞,煤田开发产量的预测,工程论坛,中国科技信息2006年第3期;
4.孙祥,徐流美,吴清,MATLAB7.0基础教程,清华大学出版社,2005年5月第一版。
附录
1.MATLAB实现灰色预测程序
function[y,p,e]=huise_1_1(X,k)%灰色模型的malab程序
%Example[y,p]=gm_1_1([200250300350],2)
%接口描述:
X的预测的初始数列,|X|>
4,K是指向后进行预测的个数
%命令格式:
程序保存的文件名,eg:
huise.m则命令是:
huise([579.8547.5527.0492.3437.0],5)
ifnargout>
3;
error('
Toomanyoutputargument.'
);
end
ifnargin==1,k=1;
x_orig=X;
elseifnargin==0|nargin>
Wrongnumberofinputarguments.'
end
x_orig=X;
predict=k;
%AGO处理,即是对初始数列进行一阶累加
x=cumsum(x_orig);
%计算系数(a和u)------------------------
n=length(x_orig);
%生成矩阵B
fori=1:
(n-1);
B(i)=-(x(i)+x(i+1))/2;
B=[B'
ones(n-1,1)];
%生成矩阵Y
y(i)=x_orig(i+1);
Y=y'
;
%计算系数a=au
(1)u=au
(2)
au=(inv(B'
*B))*(B'
*Y);
%--------------------------------------------------------
%把huise模型公式转换成符号
coef1=au
(2)/au
(1);
coef2=x_orig
(1)-coef1;
coef3=0-au
(1);
costr1=num2str(coef1);
costr2=num2str(abs(coef2));
costr3=num2str(coef3);
eq=strcat(costr1,'
+'
costr2,'
e^'
costr3,'
*(t-1))'
%计算每一个值
fort=1:
(n+predict)
mcv(t)=coef1+coef2*exp(coef3*(t-1));
x_mcv0=diff(mcv);
x_mcve=[x_orig
(1)x_mcv0]%输出图形中的各点
x_mcv=diff(mcv(1:
end-predict));
x_orig_n=x_orig(2:
end);
x_c_error=x_orig_n-x_mcv;
x_error=mean(abs(x_c_error./x_orig_n));
ifx_error>
0.2%相对误差的均值
disp('
modeldisqualification!
'
elseifx_error>
0.1
disp('
modelcheckout'
else
modelisperfect!
end
plot(1:
n,x_orig,'
o'
1:
n+predict,x_mcve);
p=x_mcve(end-predict+1:
%画出预测模型和初始数列的点
xlabel('
年份(从第一个数据年份起)'
ylabel('
产水量(万吨)'
title('
灰度模型GM(1,1)'
gridon
y=eq;
e=x_error;
2.MATLAB实现灰色预测的预测数据与实际数据对比及精度评价
“七五”井产煤量预测与相关参数
“八五”井产煤量预测与相关参数
“九五”井产煤量预测与相关参数
“七五”井产水量预测与相关参数
“八五”井产水量预测与相关参数
“九五”井产水量预测与相关参数
多目标规划主要程序段(Matlab):
M文件:
functionf=myfun(x)
f
(1)=-0.9944*x
(1)-0.7693*x
(2)-0.5659*x(3)-0.298*x(4)-0.082*x(5)-0.5044*x(6)-0.4163*x(7)-0.3159*x(8)-0.168*x(
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 煤田 开发 规划 合理 编制 问题
![提示](https://static.bingdoc.com/images/bang_tan.gif)