水力采煤生产管理专家模拟系统

1、 1999 年 11 月 系统工程理论与实践 第 11 期 水力采煤生产管理专家模拟系统孟祥瑞( 淮南工业学院资源开发与管理工程系, 安徽 淮南 232001)摘要: 结合大屯煤电公司孔庄煤矿的具体条件, 在全面系统分析的基础上, 总结大量专家经验, 以 高产高效为总体目标, 将生产管理中的主要内容建成知识库, 利用面向对象的方法, 将专家系统与模 拟系统集成构建专家模拟系统. 根据不同生产目标, 得出可靠性较高的决策值, 部分结果应用于生产 实际, 效果明显.关键词: 水力采煤; 高产高效; 生产管理; 专家模拟系统E xp e r t S im u la t io n Sy s tem o

2、 f H yd rau lic M in ingP ro duc t io n M anagem en tM EN G X ian g 2ru i(H ua inan In st itu te o f T ech no lo gy, H ua inan 232001)A bstrac t: Co n side r ing rea lit ie s o f Ko ngzh uang Co a l M ine, D a tun Co a l E lec t r ic ity Com p any, th e p ap e r sum s up a la rge am o un t o f exp e

3、 r t s exp e r ience o n th e ba sis o f o ve r2 a ll sy stem a t ic ana ly sis, tak ing h igh deg ree o f p ro duc t io n and eff ic iency a s th e gene ra l go a l, deve lop ing k now ledge2ba se abo u t co n ten t s o f p ro duc t io n m anagem en t. M ean2 w h ile, by m ean s o f o b jec t2o r i

4、en ted p ro g ramm ing , exp e r t sim u la t io n sy stem is in teg ra ted by exp e r t sy stem and sim u la t io n sy stem. M uch m o re re liab le dec isio n s a re m ade acco rd2 ing to d iffe ren t p ro duc t io n o b jec t s. A nd po sit ive effec t s a re m ade af te r som e re su lt s be2 in

5、g app lied in p rac t ice.Keywords: h yd rau lic m in ing; h igh p ro duc t io n and h igh eff ic iency; p ro duc t io n m anage2m en t; exp e r t sim u la t io n sy stem1引言水力采煤生产管理问题是一个复杂的、非线性、半结构化或非结构化问题. 作为面对复杂多变的客观 地质条件, 工作地点分散且移动频繁, 时空关系转变很快, 机械化程度较高、生产连续性很强的水采工艺来 讲, 其系统的综合生产能力最终将受到薄弱环节的制约. 在矿井

6、这一复杂的大系统中, 采掘生产系统是其 核心部分, 它的正常运转关系到煤炭生产正常、稳定和持续地进行. 对于综采、高档普采等机械化采煤工 艺, 国内外许多学者已从不同侧面对采掘生产管理进行了大量研究, 并把计算机技术引进推广, 优化理论 与方法亦应用到设计研究和生产管理中. 但对水力采煤生产系统研究却很少, 对诸多系统参数的确定, 系 统分析和生产计划制订等, 长期以来凭借经验类比和传统的决策方法.通过紧密结合大屯煤电公司孔庄矿生产实际, 在进行全面系统分析的基础上, 总结大量专家经验, 以 高产高效为总体目标, 将生产管理中主要内容建成知识库, 利用面向对象的方法, 将专家系统和模拟系统 集

7、成为专家模拟系统, 解决水采领域半结构化或非结构化问题. 根据不同产量目标, 得出可靠性较高的决策值, 为现场应用提供有价值的参考依据. 收稿日期: 199820420442系统工程理论与实践1999 年 11 月2水力采煤生产管理专家模拟系统 ESS 的建造211 系统结构本系统采用面向对象程序设计思想开发而成, 以水力采煤高产高效为父类对象, 包括采垛参数、落煤 时间、落煤能力、产量进尺预测、顶板及瓦斯管理五个子类对象. 对象既是知识的基本元素, 又是问题求解 的独立单元. 各种类型的求解机制分布于各对象, 通过对象之间的信息传递完成整个问题求解过程. 在专 家系统运行中, 调用模拟系统,

8、 通过系统模拟检验专家系统的输出结果, 从而构成整个专家模拟系统 ( E x2 p e r t S im u la t io n Sy stem , 简称“E SS”).1) 对象的表示整个系统以对象的层次结构为中心, 该结构既是知识源, 也是工作存储区, 同时又是问题求解单元的 集合. 系统将问题求解中涉及的概念、实体等表示为对象, 如图 1 所示.水力采煤高产高效类名采垛参数落煤能力落煤时间超类子类 类1 实例类2槽规 则图1 对象的层次结构与对象的组成方法各对象以它们之间的超类、子类和实例的关系形成一个层次网络. 对象的各个“槽”记录着对象的有关 属 性, 涉及该对象的操作记录在“方法”

9、中, 对象的有关规则记录在“规则”中. “槽”和“方法”中的内容都具 有继承性. 系统采用多重继承的方式, 即一个对象可以有多个超类, 采用多重继承可以更大程度地提高信 息共享. 例如, 水采生产管理中很重要的一个内容是确定合理的采垛参数, 它既需要定量计算, 也需要有关专家经验, 能够根据煤层条件及设备参数选择合理的采垛尺寸. 对象的层次结构表示如图 2 所示.采垛参数冲采角最小掐枪步距开采眼间距60 6565 7070 754m 5m 6m 7m12 1414 1616 1818 20图2 采垛参数对象的层次结构每个对象的方法和规则都封装在对象之中, 对象中方法的启动只能通过对象之间的信息

10、传递. 对象中 涉及的方法可以是规则推理, 也可以是任何求解功能, 如对关系数据库 dBA SE 或 Fo xba se 的访问, 模拟神 经网络, 以及用户编程定义的过程等.由于对象的封装性, 各类型的求解操作不会互相干扰.2) 方法和规则的表示 方法是一个程序块封装在对象之中, 其编写格式为:第 11 期水力采煤生产管理专家模拟系统43func t io n: (类型) (方法) 参数语句end (方法名)方法的类型分“P ub lic”和“P r iva te”. 前者表示该方法可被子类对象调用, 后者只能由本对象调用. 语句 部分就是方法的定义体. 为定义方法, 系统提供了一种程序设计

11、语言 O ECL , 以供使用.系统开始工作时, 在第一个对象中查找方法“m a in ”并执行之. 如果有多个方法需要执行, 则方法名应 为“m a in 1”,“m a in 2”, . . . , 系统按顺序执行之. 方法类型“P r iva te”一般与“a_ dem o n ”或“m _ dem o n ”槽对应, 由方法求值后再返回槽为“a_ dem o n ”的属性.从目前用于专家系统的知识表达的方法来看, 基于规则的表示是一种最常用的方式. 一个对象中可以 列有多条规则, 规则的形式为:ru le: 规则名 if 前提子句 th en 结论子句其中, 规则名可以省略, 不同的规

12、则也可以有相同的规则名. 前提子句有两种型式:型式 1: and (对象) is (对象名)型式 2: and 属性 o f (对象) (关系) (值)型式 1 的前提子句中, 对象为规则变量, 前提子句的功能是将对象名赋给该变量. 规则变量以问号“? ” 开头的是局部变量, 其值只在该规则的执行过程中有效. 型式 2 的前提子句中, 对象可以是具体的对象名, 也可以是规则变量, 该变量的值由型式 1 的前提子句赋予. 关系可以是“is”或“isn t”表示是与不是.若 对象的属性值为数字, 关系还可为“= ”、“ ”、“ ”、“ ”. 值可以有多个, 彼此之间用“o r”连接.值既可以是数字

13、和词汇, 也可以是规则变量.3) 数据库的建立及访问专家系统中最重要的组成部分之一是数据库, 它既可作为独立的对象, 也可用来收集关于咨询过程获 取的数据信息, 记录中间结果. 利用 dBA SE 或 Fo xba se, 首先建立了原始生产记录表 ( 见表 1) , 故障时间表 (见表 2) 和统计表 (见表 3) 及可靠度计算表 (见表 4).表 1 原始生产记录表日期班次开始时间停枪时间停枪原因煤层号rqbcsta r t im ef in ish t im ecau sem listDC 2N 5. 2N 5. 2C 6N 2 3 D 2da te 日期型; C 2 ch a rac

14、t ic 字符型; N 2 num be r 数字型, 以下同.表 2 故障时间表高压泵溜槽水枪振动筛煤水仓信号电器顶板供水管路yb tlc tsq tzc tm c txh tdb tg ltN 6. 2N 6. 2N 6. 2N 6. 2N 6. 2N 6. 2N 6. 2N 6. 2表 3 统计表班内检修时间落煤生产时间故障时间改面时间交接时间总计jx tlm tgztgm tjjtto ta lN 7. 2N 7. 2N 7. 2N 7. 2N 7. 2N 8. 2表 4 可靠度计算表名称平均故障时间平均运行时间可靠度nam eM TBMM TB FA iC 6N 8. 4N 8. 4

15、N 6. 444系统工程理论与实践1999 年 11 月对数据库进行访问操作, 则由以下函数实现:1) op enf ile ( 文件,类型)打开文件, 文件的存取方式由类型指定.2) c lo sef ile ( 文件)关闭文件, 结束访问.3) ge tda ta ( 关系数据库文件,字段名,记录号)读取文件中由字段名和记录号指定的值, 并返回之.4) labe l( 关系数据库文件,记录号)若关系数据库文件中记录号指定的记录已被删除, 则返回真, 否则返回假.5) listdb ( 文件名,字段名,记录号)在窗口 tex t0 的当前光标位置上, 显示文件中字段名和记录号指定的值. 若字

16、段名为 a ll, 则列出所 有字段的值.6) db leng th (关系数据库文件)返回 文件中的记录数目.212 系统可靠度计算系统可靠性分析是生产管理中重要的一环, 是进行有效落煤时间参数分析和提高水枪落煤能力的基 础. 其可靠度计算以对“生产原始记录表”数据库进行访问为基础, 根据各种停枪原因及开、停枪时间点, 分 别自动统计出交接班时间、改面时间、检修时间及故障时间. 然后计算出部件可靠度、系统有效度, 并找出 关键部件.水力采煤正规循环应是每班冲采一个或两个采垛或者两班冲采一个采垛, 但由于采垛参数的不合理 或因顶板冒落中途停采, 造成临时换面冲采. 对于改面时间的统计, 只要停

17、枪原因是由于改面或换面, 就可 以将下次开枪时间与本次停枪时间的差值累计. 但要注意两点: 一是两次停、开枪时间是否为同一班次, 如 果为不同班次, 则为正常的两面交替开采, 不属于因改面而造成的故障时间; 二是在生产过程中, 由于本班 第一次开枪时间较晚, 偶尔发生直到下班人员已到工作面但本班尚未停采的现象, 这种现象亦算作改换 面, 只是不必停枪停泵. 这一般发生在三个交接班的时间点, 孔庄矿为 0 点、8 点和 16 点. 因此, 统计改面时 间的方法如下:o b jec t: “y s. dbf”func t io n: p ub lic m a in op enf ile ( se

18、lf, db )op enf ile “(t jb. dbf”, db )p u t db leng th ( se lf) in to num rep ea t w ith i = 1 to nump u t ge tda ta ( se lf, 班次, i ) in to bcp u t ge tda ta ( se lf, 停枪时间, i ) in to f tp u t ge tda ta ( se lf, 停枪原因, i ) in to cau seif (cau se = 改换面 and f t 0. 0 and f t 8. 0 and f t 16. 0) th enp u t

19、 i in to nex tp u t (nex t + 1) in to jp u t ge tda ta ( se lf, 班次, j ) in to bc1p u t ge tda ta ( se lf, 开始时间, j ) in to st if (bc1 = bc) th enp u t ( st - f t) in to gm t end if第 11 期水力采煤生产管理专家模拟系统45end ifp u t (gm t + sum ) in to sum p u t ( k + 1 ) in to kend rep ea tp u tda ta “(p u tda ta “(t

20、jb. dbf”, sum , 改面时间, 1)t jb. dbf”, k , 改面次数, 2)c lo sef ile ( se lf)c lo sef ile “(end m a int jb. dbf”)213 产量进尺预测生产管理中一项重要内容是制订作业计划, 而对产量与进尺的科学预测是制订生产作业计划的基础. 根据研究结果, 采用灰色马尔可夫预测模型对水采产量、进尺的预测较为适宜. 为了将该方法应用到专家 模拟系统中, 首先将预测程序编译为执行文件 GM M . EX E , 而后通过专家系统调用. 生产历史数据记录在 数据文件中, 执行文件运行过程中访问该数据文件, 而后输出相应预

21、测结果. 同时, 考虑在煤矿生产中, 地 质条件的不确定性及对产量的决定性影响, 在产量预测中会根据地质条件变化有一定的调整系数.灰色马尔可夫预测程序框图如图 3 所示.数据文件根据不同采掘队分别建立, 分别进行预测. 一 般的灰色预测是按现实时刻 t = n 及过去的全体数据建模, 其模型是连续的时间函数. 从理论上讲, 该模型可从初值一 直延伸到未来任何一个时刻. 不过对于现实的本征性灰系 统, 随着时间的推移, 未来的一些扰动因素, 如采矿中的煤 层地质条件、设备参数变化等, 将不断地相继进入系统造成 影响, 因此需动态地将每一个新得到的数据记入数据文件, 重新建模重新预测, 这便是所谓

22、信息模型. 但是, 随着时间 的变化, 信息越来越多 (如月产量月进尺每月一个信息) , 这 样便要求计算机内存容量随着时间推移不断扩大, 运算量 不断增加, 这从经济和技术上看显然是不合理的. 此外, 随着时间的推移, 老的数据将越来越不适应新的情况, 或者说老数据的信息意义将随时间推移而降低, 因此每补充一个 信息便去掉一个最老的数据, 以维持数据的个数, 显然是合理的. 因此, 所采用的数据模型是这种新陈代谢模型, 具体开 始读取数据文件累加生成及矩阵运算X(1) (0)(k + 1) 生成模型及还原模型 X(k + 1)残 差 检 验误差 e (k ) TF建立残差灰色模型状 态 划

23、分计算状态转移概率矩阵输 出 预 测 值结 束图 3 灰色马尔可夫预测程序框图体现在每个数据文件中, 就是补充一个新数据, 去掉一个老数据, 用户操作较为简单方便, 无须重复输入.同时, 通过模糊评判矩阵计算, 地质条件对产量的影响程度分为四个等级: 好, 一般, 差, 很差. 其调整 系数分别为 1. 2, 1. 0, 0. 9 和 0. 75.214 采垛参数及巷道布置优化优化采垛参数及巷道布置是以煤层条件、供水落煤设备参数和顶板稳定性为基础的, 采垛参数的内容 包括冲采角, 最小移枪步距和工作面长度, 其推理路径如图 4 所示.215落煤能力计算水枪生产落煤能力主要取决于现场设备参数对煤

24、层地质条件的适应程度, 其推理路径如图 5 所示. 由 此, 可根据一定的产量目标及工艺水平, 确定经济合理的落煤能力; 亦可依当前设备状况和煤层条件, 优化 水枪落煤参数, 以达到最佳落煤能力.216 有效落煤时间参数分析影响有效落煤时间的因素主要有两大方面: 一是生产系统设备可靠性; 二是工艺之间的平衡. 应从以46系统工程理论与实践1999 年 11 月采 垛 参 数冲采角最小移枪步距回采巷间距落煤方法煤层条件管子长度顶板允许暴露面积有效射程工作射程小阶段漏斗式倾角煤厚采深顶板喷嘴直径水枪压力煤层裂隙图4 采垛参数对象推理路径图落煤能力水枪出口压力冲采流量采场运输条件泵压管阻煤的水力破碎

25、强度喷嘴直径煤层条件煤水比矿压作用煤硬度裂隙发育脆韧性夹矸煤厚采深采场周围情况顶板岩性图5 落煤能力对象推理路径图下三个方面考虑提高水枪有效落煤时间: 提高每个部件可靠度从而提高整个系统有效度; 研究煤水平衡, 确保煤水流畅通; 保证水采采掘平衡, 从而保证退采正常进行. 该对象层次结构如图 6 所示.落煤时间采掘平衡系统可靠性煤水平衡采后掘采前掘水平垂高煤水泵流量煤水仓容积退采速度掘进速度煤层条件冲采浆量掘进浆量图6 落煤时间对象的层次结构图217顶板管理和瓦斯管理根据分析, 顶板是目前孔庄矿水采生产管理最薄弱环节, 瓦斯是影响水采安全生产的最不利因素. 为 此特在专家系统中建立一个对象,

26、专门阐述顶板管理和瓦斯管理的定性分析、定量计算和专家经验.第 11 期水力采煤生产管理专家模拟系统473“ESS”组成及其功能311 ESS 系统组成该专家模拟系统是水力采煤生产管理专家系统和模拟系统的集成, 采用面向对象的程序设计的思想 开发而成. 程序 O EC. EX E、知识库编译文件及模拟系统构成了 E SS 的基本成分. E SS 包括的主要文件如表 5 所示.3. 2ESS 系统运行环境1) IBM - PC 类微机及其兼容机;2 ) U CDO S, CCDO S, WM 等 中 文 操 作 系 统.313 ESS 系统功能“E SS”可实现以下功能:1) 通过对“生产原始记录

27、”数据库文件的访 问, 可自动计算系统部件可靠度和系统有效度, 找出系统关键环节, 进行系统有效度按重要度 重新分配;2) 根据煤层条件及开采方法, 优化采垛参 数和采面巷道布置;3) 以水力采煤高产高效为总体目标, 进行 有效落煤时间参数分析, 进行采掘平衡和煤水平衡的定量计算和定性分析, 提出有效的平衡表 5O KD. EX E 知识库开发系统O EC. EX E 推理系统W S. COM 文本编辑程序A N N B. EX E 神经网络模拟程序GM M . EX E 灰色马尔可夫预测程序JJ. DA T 掘进进尺数据文件CM . EX E 采煤产量数据文件H IGH PRO. CM P

28、水采高产高效知识库编译文件SYST EM . CM P 系统可靠度计算编译文件 FO RCA ST. CM P 产量进尺预测知识库文件 YS. DB F 生产原始记录表数据库文件 H EA D. DA T 标题文件 GRA PH. DA T 显示库文件 措施, 确定合理的煤水仓容积和煤水泵流量, 以及在一定掘进速度下的采掘面头比、合理退采速度等;4) 通过系统向用户提问, 根据水采生产系统的设备参数及环节能力, 核定水枪落煤能力, 推荐合理落 煤压力和合理流量值;5) 调用外部文件执行过程, 运用灰色马尔柯夫预测方法, 对产量、进尺进行精度较高的预测, 为生产作 业计划的制订提供科学依据;6)

29、 对目前水采生产系统的最薄弱环节顶板的管理提出专家级的管理措施;7) 在对瓦斯涌出定量计算的基础上, 确定其对产量、水枪有效落煤时间的约束程度, 提出综合防治方法.4 水力采煤生产管理专家模拟系统 ESS 应用实例 根据大屯煤电公司孔庄矿生产实际, 应用 E SS 进行水采高产高效的生产管理决策咨询. 首先系统向用户提问, 获取基本煤层地质条件和生产技术条件后, 即可利用知识库中专家经验进行定性分析和定量计算. 系统依次提问如表 6 所示.48系统工程理论与实践1999 年 11 月表 61) 请问目前开采深度为多少米? 选择: 200, 200- 400, 400- 6002) 请选择所咨询

30、的煤层号: 选择: 7# , 8# , 其它3) 您对水采年产量的期望是: 选择: 50 万 t, 60 万 t, 70 万 t, 80 万 t4) 该矿水采面落煤方法是: 选择: 走向小阶段式, 倾斜漏斗式5) 请问水采面回采顺槽长度为多少米? 选择: 80- 100, 100- 120, 60- 80, 120- 1406) 请问工作面煤水溜槽波度为多少? 输入一个数.7) 该煤层裂隙发育情况是: 选择: 中等发育; 发育; 不发育8) 该煤层的脆韧性是: 选择: 较脆; 脆; 韧9) 请问煤层硬度普氏系数是多少? 输入一个值.10) 煤层有无夹矸: 选择: 有; 没有11) 请选择煤层顶

31、板类别: 选择: 中等稳定; 坚硬顶板; 松软顶板12) 请问煤水提升方式是: 选择: 一次提升; 分级提升 在不同年产量 ( 50 万 ta、60 万 ta、70 万 ta 、80 万 ta ) 的目标下, 根据孔庄矿生产实际情况, 利用 “人2机”交互方式逐一回答系统提问, 分别对 7# 、8# 煤层的合理采垛参数、落煤能力及落煤时间提出建 议. 经过规则匹配和推理机推理, 最后系统得出结论如表 7 所示.表 7 不同产量目标下的落煤时间及落煤能力分析计算表产量目标煤层落煤能力落煤压力合理流量落煤时间退采速度面头比煤水泵煤水仓507#1501453305. 511. 432. 651028

32、08#1301602703. 07. 411. 2全矿1408. 53. 8607#160145350613. 723. 35303008 881. 5全矿150104. 8707#1701453707. 515. 693. 75603108#1451803005. 010. 161. 8全矿16012. 55. 5807#1851603708. 017. 784. 16003308#1501903005. 511. 121. 9全矿17019013. 56. 0 对于系统可靠度计算, 需首先按照“生产原始记录表”逐日逐班输入计算机, 以数据库文件形式保存. “E S

33、S”系统通过对该数据库的访问, 即可自动计算部件、系统可靠度, 找出系统薄弱环节, 并将部件故障时 间、故障次数记入数据库文件中, 以供查询. 例如, 我们将 95 年 7 月水采生产记录输入计算机, 进入专家模 拟系统后, 操作“O EC P S”, 经运算几秒后, 给出结果为:1) 系统目前的最薄弱环节是顶板, 其有效度是: 0. 8821;2) 系统有效度为: 0. 7561, 其它各部件有效度分别如表 8.第 11 期水力采煤生产管理专家模拟系统49表 8部件名称 有效度故障时间故障次数 顶板0. 8821 1314 38 高压泵 0. 9912 79 5 煤水仓 0. 9643 33

34、1 5 水枪0. 9908 81 4瓦斯0. 9982 13 1 振动筛0. 9565 402 5 煤水溜槽0. 9913 70 3 信号电器0. 9955 31 2 供水管路0. 9588 413 93) 非正常改面的时间及次数是: 67m in , 3 次;4) 操作因素故障时间是: 756m in;5) 本月工作天数是 29d, 有效落煤时间为: 178. 48h.现生产条件下合理采垛参数决策数据如表 9, 系统改造前后两时期可靠度计算对比如表 10 所示. 表 11为 95 年 1 12 月水采进尺预测值及拟合值, 表 12 为 95 年 1 12 月水采产量预测值及拟合值.表 9 目

35、前设备参数条件下合理采垛参数煤 层冲采角 ( )最小移枪步距 (m )采垛斜长 (m )7#70 756. 017. 08#70 755. 016. 0表 10 系统改进前后两期可靠度计算表项 目1994 年 6 月 1995 年 6 月1995 年 7 月 1995 年 12 月故障时间 m in次数可靠度故障时间 m in次数可靠度顶 板77432490. 88761552960. 8999瓦 斯15740. 99718670. 9971振 动 筛60121360. 90572783620. 9630溜 槽31071020. 9481625400. 9778水 枪523290. 98258

36、72270. 9878煤 水 仓2643360. 97412701440. 9638信 号28780. 9962392120. 9942高 压 泵36881060. 97033230980. 9575供水管路1083280. 98281759440. 9761系统合计0. 69310. 7436表 11 1995 年 1 12 月水采进尺预测值及拟合值项 目123456789101112预测值 (m )639788666813694839723866879768799937实际值 (m )630806634760634874768913854770774871相对误差%1. 4- 2. 25. 06. 99. 4- 4. 0- 5. 9- 5. 22. 9- 0. 33. 27. 650系统工程理论与实践1999 年 11 月表 12 1995 年 1 12 月水采产量预测值及拟合值项 目123456789101112预测值 ( 万 t)2. 832. 893. 443. 113. 193. 663. 713. 853. 934. 314. 284. 05实际值 ( 万 t)2. 992. 633. 563. 273. 373. 573. 953. 673. 824. 394. 143. 76相对误差%- 5. 49. 9