1999年11月系统工程理论与实践第11期水力采煤生产管理专家模拟系统α孟祥瑞(淮南工业学院资源开发与管理工程系,安徽淮南232001摘要:结合大屯煤电公司孔庄煤矿的具体条件,在全面系统分析的基础上,总结大量专家经验,以高产高效为总体目标,将生产管理中的主要内容建成知识库,利用面向对象的方法,将专家系统与模拟系统集成构建专家模拟系统.根据不同生产目标,得出可靠性较高的决策值,部分结果应用于生产实际,效果明显.关键词:水力采煤;高产高效;生产管理;专家模拟系统Expert Si m u lati on System of H ydrau lic M in ingP roducti on M anagem en tM EN G X iang2ru i(H uainan In stitu te of T echno logy,H uainan232001Abstract:Con sidering realities of Kongzhuang Coal M ine,D atun Coal E lectricityCompany,the paper sum s up a large amoun t of expert’s experience on the basis of over2all system atic analysis,tak ing h igh degree of p roducti on and efficiency as the generalgoal,develop ing know ledge2base abou t con ten ts of p roducti on m anagem ent.M ean2w h ile,by m ean s of ob ject2o rien ted p rogramm ing,expert si m u lati on system is in tegratedby expert system and si m u lati on system.M uch mo re reliab le decisi on s are m ade acco rd2ing to differen t p roducti on ob jects.A nd po sitive effects are m ade after som e resu lts be2ing app lied in p ractice.Keywords:hydrau lic m in ing;h igh p roducti on and h igh efficiency;p roducti on m anage2m en t;expert si m u lati on system1引言水力采煤生产管理问题是一个复杂的、非线性、半结构化或非结构化问题.作为面对复杂多变的客观地质条件,工作地点分散且移动频繁,时空关系转变很快,机械化程度较高、生产连续性很强的水采工艺来讲,其系统的综合生产能力最终将受到薄弱环节的制约.在矿井这一复杂的大系统中,采掘生产系统是其核心部分,它的正常运转关系到煤炭生产正常、稳定和持续地进行.对于综采、高档普采等机械化采煤工艺,国内外许多学者已从不同侧面对采掘生产管理进行了大量研究,并把计算机技术引进推广,优化理论与方法亦应用到设计研究和生产管理中.但对水力采煤生产系统研究却很少,对诸多系统参数的确定,系统分析和生产计划制订等,长期以来凭借经验类比和传统的决策方法.通过紧密结合大屯煤电公司孔庄矿生产实际,在进行全面系统分析的基础上,总结大量专家经验,以高产高效为总体目标,将生产管理中主要内容建成知识库,利用面向对象的方法,将专家系统和模拟系统集成为专家模拟系统,解决水采领域半结构化或非结构化问题.根据不同产量目标,得出可靠性较高的决策值,为现场应用提供有价值的参考依据.α2水力采煤生产管理专家模拟系统ESS 的建造211系统结构本系统采用面向对象程序设计思想开发而成,以水力采煤高产高效为父类对象,包括采垛参数、落煤时间、落煤能力、产量进尺预测、顶板及瓦斯管理五个子类对象.对象既是知识的基本元素,又是问题求解的独立单元.各种类型的求解机制分布于各对象,通过对象之间的信息传递完成整个问题求解过程.在专家系统运行中,调用模拟系统,通过系统模拟检验专家系统的输出结果,从而构成整个专家模拟系统(Ex 2pert Si m u lati on System ,简称“ESS ”.1对象的表示整个系统以对象的层次结构为中心,该结构既是知识源,也是工作存储区,同时又是问题求解单元的集合.系统将问题求解中涉及的概念、实体等表示为对象,如图1所示.图1对象的层次结构与对象的组成类1类2实例方法规则槽子类超类类名落煤时间采垛参数落煤能力水力采煤高产高效各对象以它们之间的超类、子类和实例的关系形成一个层次网络.对象的各个“槽”记录着对象的有关属性,涉及该对象的操作记录在“方法”中,对象的有关规则记录在“规则”中.“槽”和“方法”中的内容都具有继承性.系统采用多重继承的方式,即一个对象可以有多个超类,采用多重继承可以更大程度地提高信息共享.例如,水采生产管理中很重要的一个内容是确定合理的采垛参数,它既需要定量计算,也需要有关专家经验,能够根据煤层条件及设备参数选择合理的采垛尺寸.对象的层次结构表示如图2所示.图2采垛参数对象的层次结构18~2016~1814~1612~147m 6m 5m 4m 70~7565~7060~65开采眼间距冲采角最小掐枪步距采垛参数每个对象的方法和规则都封装在对象之中,对象中方法的启动只能通过对象之间的信息传递.对象中涉及的方法可以是规则推理,也可以是任何求解功能,如对关系数据库dBA SE 或Foxbase 的访问,模拟神经网络,以及用户编程定义的过程等.由于对象的封装性,各类型的求解操作不会互相干扰.2方法和规则的表示方法是一个程序块封装在对象之中,其编写格式为:24系统工程理论与实践1999年11月functi on :(类型(方法{参数}{语句}end (方法名方法的类型分“Pub lic ”和“P rivate ”.前者表示该方法可被子类对象调用,后者只能由本对象调用.语句部分就是方法的定义体.为定义方法,系统提供了一种程序设计语言O ECL ,以供使用.系统开始工作时,在第一个对象中查找方法“m ain ”并执行之.如果有多个方法需要执行,则方法名应为“m ain 1”,“m ain 2”,...,系统按顺序执行之.方法类型“P rivate ”一般与“a _demon ”或“m _demon ”槽对应,由方法求值后再返回槽为“a _demon ”的属性.从目前用于专家系统的知识表达的方法来看,基于规则的表示是一种最常用的方式.一个对象中可以列有多条规则,规则的形式为:ru le :[规则名]if 前提子句then 结论子句其中,规则名可以省略,不同的规则也可以有相同的规则名.前提子句有两种型式:型式1:[and ](对象is (对象名型式2:[and ]属性of (对象(关系(值型式1的前提子句中,对象为规则变量,前提子句的功能是将对象名赋给该变量.规则变量以问号“?”开头的是局部变量,其值只在该规则的执行过程中有效.型式2的前提子句中,对象可以是具体的对象名,也可以是规则变量,该变量的值由型式1的前提子句赋予.关系可以是“is ”或“isn’t ”表示是与不是.若对象的属性值为数字,关系还可为“=”、“>”、“<”、“<>”.值可以有多个,彼此之间用“o r ”连接.值既可以是数字和词汇,也可以是规则变量.3数据库的建立及访问专家系统中最重要的组成部分之一是数据库,它既可作为独立的对象,也可用来收集关于咨询过程获取的数据信息,记录中间结果.利用dBA SE 或Foxbase ,首先建立了原始生产记录表(见表1,故障时间表(见表2和统计表(见表3及可靠度计算表(见表4.表1原始生产记录表日期班次开始时间停枪时间停枪原因煤层号rq bc starti m e finish ti m e cause m list DC 2N 5.2N 5.2C 6N 23D 2date 日期型;C 2charactic 字符型;N 2num ber 数字型,以下同. 表2故障时间表高压泵溜槽水枪振动筛煤水仓信号电器顶板供水管路ybt lct sqt zct m ct xh t dbt glt N 6.2 N 6.2N 6.2N 6.2N 6.2N 6.2N 6.2N 6.2表3统计表班内检修时间落煤生产时间故障时间改面时间交接时间总计jxt l m t gzt gm t jjt to tal N 7.2N 7.2N 7.2N 7.2N 7.2N 8.2表4可靠度计算表名称平均故障时间平均运行时间可靠度nam e M TBM M TBF A i C 6N 8.4N 8.4N 6.434第11期水力采煤生产管理专家模拟系统44系统工程理论与实践1999年11月对数据库进行访问操作,则由以下函数实现:1openfile(〈文件〉,〈类型〉打开文件,文件的存取方式由类型指定.2clo sefile(〈文件〉关闭文件,结束访问.3getdata(〈关系数据库文件〉,〈字段名〉,〈记录号〉读取文件中由字段名和记录号指定的值,并返回之.4label(〈关系数据库文件〉,〈记录号〉若关系数据库文件中记录号指定的记录已被删除,则返回真,否则返回假.5listdb(〈文件名〉,〈字段名〉,〈记录号〉在窗口tex t0的当前光标位置上,显示文件中字段名和记录号指定的值.若〈字段名〉为all,则列出所有字段的值.6db length(〈关系数据库文件〉返回〈文件〉中的记录数目.212系统可靠度计算系统可靠性分析是生产管理中重要的一环,是进行有效落煤时间参数分析和提高水枪落煤能力的基础.其可靠度计算以对“生产原始记录表”数据库进行访问为基础,根据各种停枪原因及开、停枪时间点,分别自动统计出交接班时间、改面时间、检修时间及故障时间.然后计算出部件可靠度、系统有效度,并找出关键部件.水力采煤正规循环应是每班冲采一个或两个采垛或者两班冲采一个采垛,但由于采垛参数的不合理或因顶板冒落中途停采,造成临时换面冲采.对于改面时间的统计,只要停枪原因是由于改面或换面,就可以将下次开枪时间与本次停枪时间的差值累计.但要注意两点:一是两次停、开枪时间是否为同一班次,如果为不同班次,则为正常的两面交替开采,不属于因改面而造成的故障时间;二是在生产过程中,由于本班第一次开枪时间较晚,偶尔发生直到下班人员已到工作面但本班尚未停采的现象,这种现象亦算作改换面,只是不必停枪停泵.这一般发生在三个交接班的时间点,孔庄矿为0点、8点和16点.因此,统计改面时间的方法如下:ob ject:“ys.dbf”functi on:pub lic m ainopenfile(self,dbopenfile(“tjb.dbf”,dbpu t db length(selfin to numrepeat w ith i=1to numpu t getdata(self,班次,iin to bcpu t getdata(self,停枪时间,iin to ftpu t getdata(self,停枪原因,iin to cau seif(cau se=改换面and ft<>0.0and ft<>8.0and ft<>16.0thenpu t‘i’in to nex tpu t(nex t+1in to jpu t getdata(self,班次,jin to bc1pu t getdata(self,开始时间,jin to stif(bc1=bcthenpu t(st-ftin to gm tend ifend ifpu t (gm t +sum in to sum pu t (k +1in to kend repeat pu tdata (“tjb .dbf ”,sum ,改面时间,1pu tdata (“tjb .dbf ”,k ,改面次数,2clo sefile (self clo sefile (“tjb .dbf ”end m ain213产量进尺预测生产管理中一项重要内容是制订作业计划,而对产量与进尺的科学预测是制订生产作业计划的基础.根据研究结果,采用灰色马尔可夫预测模型对水采产量、进尺的预测较为适宜.为了将该方法应用到专家模拟系统中,首先将预测程序编译为执行文件G MM .EXE ,而后通过专家系统调用.生产历史数据记录在数据文件中,执行文件运行过程中访问该数据文件,而后输出相应预测结果.同时,考虑在煤矿生产中,地质条件的不确定性及对产量的决定性影响,在产量预测中会根据地质条件变化有一定的调整系数.灰色马尔可夫预测程序框图如图3所示.开始读取数据文件累加生成及矩阵运算Xδ(1(k +1生成模型及还原模型X δ(0(k +1残差检验误差e (k <ΕTF建立残差灰色模型状态划分计算状态转移概率矩阵输出预测值结束图3灰色马尔可夫预测程序框图数据文件根据不同采掘队分别建立,分别进行预测.一般的灰色预测是按现实时刻t =n 及过去的全体数据建模,其模型是连续的时间函数.从理论上讲,该模型可从初值一直延伸到未来任何一个时刻.不过对于现实的本征性灰系统,随着时间的推移,未来的一些扰动因素,如采矿中的煤层地质条件、设备参数变化等,将不断地相继进入系统造成影响,因此需动态地将每一个新得到的数据记入数据文件,重新建模重新预测,这便是所谓信息模型.但是,随着时间的变化,信息越来越多(如月产量月进尺每月一个信息,这样便要求计算机内存容量随着时间推移不断扩大,运算量不断增加,这从经济和技术上看显然是不合理的.此外,随着时间的推移,老的数据将越来越不适应新的情况,或者说老数据的信息意义将随时间推移而降低,因此每补充一个信息便去掉一个最老的数据,以维持数据的个数,显然是合理的.因此,所采用的数据模型是这种新陈代谢模型,具体体现在每个数据文件中,就是补充一个新数据,去掉一个老数据,用户操作较为简单方便,无须重复输入.同时,通过模糊评判矩阵计算,地质条件对产量的影响程度分为四个等级:好,一般,差,很差.其调整系数分别为1.2,1.0,0.9和0.75.214采垛参数及巷道布置优化优化采垛参数及巷道布置是以煤层条件、供水落煤设备参数和顶板稳定性为基础的,采垛参数的内容包括冲采角,最小移枪步距和工作面长度,其推理路径如图4所示.215落煤能力计算水枪生产落煤能力主要取决于现场设备参数对煤层地质条件的适应程度,其推理路径如图5所示.由此,可根据一定的产量目标及工艺水平,确定经济合理的落煤能力;亦可依当前设备状况和煤层条件,优化水枪落煤参数,以达到最佳落煤能力.216有效落煤时间参数分析影响有效落煤时间的因素主要有两大方面:一是生产系统设备可靠性;二是工艺之间的平衡.应从以54第11期水力采煤生产管理专家模拟系统图4采垛参数对象推理路径图最小移枪步距回采巷间距工作射程有效射程顶板允许暴露面积管子长度煤层裂隙水枪压力喷嘴直径顶板采深煤厚倾角漏斗式小阶段煤层条件落煤方法冲采角采垛参数图5落煤能力对象推理路径图顶板岩性采场周围情况采深煤厚夹矸脆韧性裂隙发育煤硬度矿压作用煤水比煤层条件喷嘴直径煤的水力破碎强度管阻泵压采场运输条件冲采流量水枪出口压力落煤能力下三个方面考虑提高水枪有效落煤时间:①提高每个部件可靠度从而提高整个系统有效度;②研究煤水平衡,确保煤水流畅通;③保证水采采掘平衡,从而保证退采正常进行.该对象层次结构如图6所示.图6落煤时间对象的层次结构图煤层条件掘进速度退采速度水平垂高采前掘采后掘掘进浆量冲采浆量煤水泵流量煤水仓容积煤水平衡采掘平衡系统可靠性落煤时间217顶板管理和瓦斯管理根据分析,顶板是目前孔庄矿水采生产管理最薄弱环节,瓦斯是影响水采安全生产的最不利因素.为此特在专家系统中建立一个对象,专门阐述顶板管理和瓦斯管理的定性分析、定量计算和专家经验.3“ESS ”组成及其功能311ESS 系统组成该专家模拟系统是水力采煤生产管理专家系统和模拟系统的集成,采用面向对象的程序设计的思想开发而成.程序O EC .EXE 、知识库编译文件及模拟系统构成了ESS 的基本成分.ESS 包括的主要文件如表5所示.表5O KD .EXE 知识库开发系统O EC.EXE 推理系统W S .COM 文本编辑程序ANNB .EXE 神经网络模拟程序G MM .EXE 灰色马尔可夫预测程序JJ .DA T 掘进进尺数据文件CM .EXE 采煤产量数据文件H IGH PRO .CM P 水采高产高效知识库编译文件SYST E M .CM P 系统可靠度计算编译文件FORCA ST .CM P 产量进尺预测知识库文件YS .DBF 生产原始记录表数据库文件H EAD .DA T 标题文件GRA PH .DA T显示库文件3.2ESS 系统运行环境1I BM -PC 类微机及其兼容机;2U CDO S ,CCDO S ,WM 等中文操作系统.313ESS 系统功能“ESS ”可实现以下功能:1通过对“生产原始记录”数据库文件的访问,可自动计算系统部件可靠度和系统有效度,找出系统关键环节,进行系统有效度按重要度重新分配;2根据煤层条件及开采方法,优化采垛参数和采面巷道布置;3以水力采煤高产高效为总体目标,进行有效落煤时间参数分析,进行采掘平衡和煤水平衡的定量计算和定性分析,提出有效的平衡措施,确定合理的煤水仓容积和煤水泵流量,以及在一定掘进速度下的采掘面头比、合理退采速度等;4通过系统向用户提问,根据水采生产系统的设备参数及环节能力,核定水枪落煤能力,推荐合理落煤压力和合理流量值;5调用外部文件执行过程,运用灰色马尔柯夫预测方法,对产量、进尺进行精度较高的预测,为生产作业计划的制订提供科学依据;6对目前水采生产系统的最薄弱环节顶板的管理提出专家级的管理措施;7在对瓦斯涌出定量计算的基础上,确定其对产量、水枪有效落煤时间的约束程度,提出综合防治方法.4水力采煤生产管理专家模拟系统ESS 应用实例根据大屯煤电公司孔庄矿生产实际,应用ESS 进行水采高产高效的生产管理决策咨询.首先系统向用户提问,获取基本煤层地质条件和生产技术条件后,即可利用知识库中专家经验进行定性分析和定量计算.系统依次提问如表6所示.1请问目前开采深度为多少米?选择:<200,200-400,400-6002请选择所咨询的煤层号:选择:7#,8#,其它3您对水采年产量的期望是:选择:50万t,60万t,70万t,80万t4该矿水采面落煤方法是:选择:走向小阶段式,倾斜漏斗式5请问水采面回采顺槽长度为多少米?选择:80-100,100-120,60-80,120-140 6请问工作面煤水溜槽波度为多少?输入一个数.7该煤层裂隙发育情况是:选择:中等发育;发育;不发育8该煤层的脆韧性是:选择:较脆;脆;韧9请问煤层硬度普氏系数是多少?输入一个值.10煤层有无夹矸:选择:有;没有11请选择煤层顶板类别:选择:中等稳定;坚硬顶板;松软顶板12请问煤水提升方式是:选择:一次提升;分级提升在不同年产量(50万t a、60万t a、70万t a、80万t a的目标下,根据孔庄矿生产实际情况,利用“人2机”交互方式逐一回答系统提问,分别对7#、8#煤层的合理采垛参数、落煤能力及落煤时间提出建议.经过规则匹配和推理机推理,最后系统得出结论如表7所示.表7不同产量目标下的落煤时间及落煤能力分析计算表产量目标煤层落煤能力落煤压力合理流量落煤时间退采速度面头比煤水泵煤水仓507#1501453305.511.432.68#1301602703.07.411.2全矿1408.53.8510280607#160145350613.723.38#14018029048.881.5全矿150104.8530300707#1701453707.515.693.78#1451803005.010.161.8全矿16012.55.5560310807#1851603708.017.784.18#1501903005.511.121.9全矿17019013.56.0600330对于系统可靠度计算,需首先按照“生产原始记录表”逐日逐班输入计算机,以数据库文件形式保存.“ESS”系统通过对该数据库的访问,即可自动计算部件、系统可靠度,找出系统薄弱环节,并将部件故障时间、故障次数记入数据库文件中,以供查询.例如,我们将95年7月水采生产记录输入计算机,进入专家模拟系统后,操作“O EC PS”,经运算几秒后,给出结果为:1系统目前的最薄弱环节是顶板,其有效度是:0.8821;2系统有效度为:0.7561,其它各部件有效度分别如表8.部件名称有效度故障时间故障次数顶板0.8821131438高压泵0.9912795煤水仓0.96433315水枪0.9908814瓦斯0.9982131振动筛0.95654025煤水溜槽0.9913703信号电器0.9955312供水管路0.958841393非正常改面的时间及次数是:67m in,3次;4操作因素故障时间是:756m in;5本月工作天数是29d,有效落煤时间为:178.48h.现生产条件下合理采垛参数决策数据如表9,系统改造前后两时期可靠度计算对比如表10所示.表11为95年1~12月水采进尺预测值及拟合值,表12为95年1~12月水采产量预测值及拟合值.表9目前设备参数条件下合理采垛参数煤层冲采角(°最小移枪步距(m采垛斜长(m7#70~756.017.08#70~755.016.0表10系统改进前后两期可靠度计算表项目1994年6月~1995年6月1995年7月~1995年12月故障时间m in 次数可靠度故障时间m in 次数可靠度顶板7743 2490.8876155 2960.8999瓦斯157 40.997186 70.9971振动筛6012 1360.90572783 620.9630溜槽3107 1020.9481625 400.9778水枪523 290.9825872 270.9878煤水仓2643 360.97412701 440.9638信号287 80.9962392 120.9942高压泵3688 1060.97033230 980.9575供水管路1083 280.98281759 440.9761系统合计0.69310.7436表111995年1~12月水采进尺预测值及拟合值项目123456789101112预测值(m639788666813694839723866879768799937实际值(m630806634760634874768913854770774871相对误差%1.4-2.25.06.99.4-4.0-5.9-5.22.9-0.33.27.6表121995年1~12月水采产量预测值及拟合值项目123456789101112预测值(万t2.832.893.443.113.193.663.713.853.934.314.284.05实际值(万t2.992.633.563.273.373.573.953.673.824.394.143.76相对误差%-5.49.9-3.4-4.9-5.32.5-6.14.92.9-6.13.47.7参考文献:[1]洪允和1水力采煤1北京:煤炭工业出版社,19881[2]朱克强,黎志成1C I M S环境下综合生产计划专家模拟系统的研究1系统工程理论与实践,1994,14(3:30~351[3]邓聚龙1灰色预测与决策1武汉:华中理工大学出版社,19931(上接第40页从文献[1]表1、表4、表5及计算过程看,因素B i的层次总排序U(B i值的计算过程已考虑B i的上层V21、V2的权重因素.因而B i的层次总排序U(B i是B i 对最高层(目标层V的权重,而U(B i×F(B i则反映的是因素B1,B2,…,B6之中的某元素B i对最高层(目标层V的贡献.仍以B i为例(文献[1]表5.B1:0.060×0.70=0.042所得0.042是元素B1对最高层(目标层V的贡献,而不是对其上层元素V21的贡献.B i:66(U(B i×F(B i=0.216i=1所得到的0.216并不是V21的评价值F(V21,不可用于从V2i到V2的几何综合.对A i,C i,D i亦类似.文献[1]给出的表1科技农业示范区开发效益综合评价指标体系,既有线性加权成分又有几何综合成分,是一种混合综合评价体系.其层次总排序U(i是下层元素对最高层(目标层的权重,在文献[1]设定的条件下,并不可直接用于计算中间层的评价值(文献[5].上层某元素V对其下层次相关子因素V i的层次单排序W(V i,是其子因素V1,V2,…,V n对上层某元素V的权重.下层次元素的层次单排序也不直接用于最高层(目标层的评价值计算.在文献[1]公式(3、(4中,用层次总排序U(i计算中间层次评价值,这是不正确的.513在几何综合法和混合综合法中,层次总排序U(i计算较复杂,若用层次单排序W(i从最底层起逐层向上计算,则步骤清晰.参考文献:[1]石金贵1科技农业示范区开发效益的综合评判方法1中国农业大学学报,1996,1(4:8~141[2]朱永达1农业系统工程1北京:农业出版社,19931[3]张崇甫等1统计分析方法及其应用1重庆:重庆大学出版社,19951[4]陈义华1数学模型1重庆:重庆大学出版社,19951[5]李德等1运筹学1北京:清华大学出版社,19821。
