ICS.73.100 D.96 MT 备案号：1069-1998 中华人民共和国煤行业标准 MT/T 145—1997 评定煤用重选设备工艺性能的 计算机算法 Gravity separating equipment for coal-Performance evaluationComputer algorithm 1997-12-12发布 1998-07-01实施 中华人民共和国煤炭工业部发布 MT/T145-1997 前 言 本标准是配合国家标准GB/T15715~1996《煤用重选设备工艺性能评定方法》使用的一项行业标 准，主要供算法编程人员使用。使用本标准时务必仔细阅读GB/T15715一1995。 本标准于1986年首次发布，本次为第一次修订。 本次修订除在格式方面按GB/T1.1-1993做了修改以外，主要在技术内容方面依据 GB/T15715--1995做了相应的修改，同时在算法规定方面注意到与计算机技术的进步相适应。 本标准从实施之时起代替MT145一86。 本标准的附录A、附录B和附录C都是提示的附录。 本标准由煤炭工业部科技教育司提出。 本标准由全国煤炭标准化技术委员会归口。 本标准由煤炭科学研究总院唐山分院起草并负责解释。 本标准主要起草人：李学琨。 1 中华人民共和国煤炭行业标准 评定煤用重选设备工艺性能的 MT/T145—1997 计算机算法 代警MT145-86 Gravity separating equipmentforcoal-Performance evaluation--Computeralgorithm 1范围 本标准规定了计算煤用重选设备各项工艺性能指标的计算机算法。这些指标包括可能偏差、不完善 度、数量效率和错配物总量等。 本标准适用于煤炭的重介质和水介质分选作业中各种重选设备工艺性能的评定。 本标准不适用于选煤厂月综合资料的处理。 2引用标准 下列标准所包含的条文，通过在本标准中引用面构成为本标准的条文。本标准出版时，所示版本均 为有效。所有标准都会被修订，使用本标准的各方应探讨使用下列标准最新版本的可能性。 GB/T6949-86煤炭视比重测定方法 GB/T157151995煤用重选设备工艺性能评定方法 3分选产品产率的计算及数据检验 3.1在无实测数据的情况下，分选产品的产率可按GB/T15715附录A提供的方法计算。 3.2在完成产品产率、“计算原煤”的可选性数据和重产品分配率的计算以后，应按GB/T15715中4.2 至4.4的规定对计算所得的数据进行检验。对于判定为不合格的数据，非经必要的处理不能继续进行计 算。 4各密度级平均密度的确定 4.1中间密度级 4.1.1原煤和分选产品浮沉组成中第二个密度级至倒数第二个密度级的平均密度，根据质量守恒关 系，按式(1)确定： Y, -- Y.-I ·( 1 ) 式中：p,——第：个密度级的平均密度，kg/L； Y,第，个密度级的浮物累计产率，%； Y,-1一第(t一1)个密度级的浮物累计产率，%； 一密度的倒数，可视为浮物累计产率的函数。 p 中华人民共和国煤炭工业部1997-12-12批准 1998-07-01实施 1 MT/T145-1997 4.1.2函数关系 (2) 可用分段插值或分段拟合的方法确定。 4.2端部密度级 4.2.1原煤和分选产品浮沉组成中第一个密度级和倒数第一个密度级的平均密度可按GB6949测 定。 4.2.2在不具备测定值的情况下，上述两个密度值可利用“计算原煤”的可选性数据，通过中间各密度 级的平均密度对平均灰分的线性回归方程外推获得。 4.2.3如果第一个密度级平均密度的外推值高于浮沉试验的第个密度时，第一个密度级平均密度可 按下述原则取值： 当P≥1.3时，取P=P-0.02 (3) 当P<1.3时，取=p:-0.01 ·(4 ) 上两式中：P:—一浮沉试验的第一个密度，kg/L； Pi—第个密度级的平均密度，kg/L。 4.2.4对于按上述方法推算获得的端部密度级平均密度，如果发生歧义，可按GB/T6949补充测定； 或由合同双方依据经验协商取值。 5曲线绘制的一般原则 5.1用以评定工艺性能的分配曲线，可选性曲线和错配物曲线的绘制应符合GB/T15715第5章之规 定。 5.2上述曲线的绘制，可采用插值方法处理，也可采用拟合方法处理。附录A和附录B分别介绍了一 种使用效果较好的插值方法和非线性拟合方法，可供选用。 5.3无论是采用插值方法还是采用拟合方法处理曲线，都要求所选用的数学模型能保证获得的曲线连 续、光滑、没有不适当的拐点，具有合理的性态。 5.4在采用拟合方法处理曲线时，允许对于同一组型值点更换若干不同的数学模型进行拟合运算，并 以拟合误差比较拟合效果的优劣。在曲线性态符合要求的前提下，应按拟合误差最小的原则选定模型。 拟合误差按下式计算： F (5) VN 式中：F。--拟合误差； N-密度级数； E,一第1个型值点的纵坐标； 3第：个型值点的拟合值。 5.5在采用拟合方法处理曲线时，选定模型的拟合误差应小于曲线纵坐标全尺度的5%，否则可认为 原始数据不合格或模型选用不当。 6曲线绘制的具体规定 6.1分配曲线 6.1.1分配曲线的型值点分为基本型值点和虚拟型值点。 6.1.2基本型值点的横坐标为按本标准第4章所述方法确定的各密度级的平均密度；纵坐标为相应各 密度级的入料在重产品中的分配率。基本型值点应通过GB/T15715一1995中4.3的规定检验。 6.1.3当基本型值点中分配率为0的型值点或分配率为100%的型值点多于两个时，为了保持拟合曲 2 MT/T 145—1997 线的性态，可以在拟合时舍去外侧的型值点。 6.1.4如基本型值点中缺少分配率小于2%或大于90%的型值点，为了保持拟合曲线的性态，可考虑 增设相应的虚拟型值点，虚拟型值点的横坐标分别设置于最低密度物的虚拟密度（pm）处或最高密度物 的虚拟密度（pmax）处，其分配率分别取为0或1(100%）。 6.1.5最低密度物的虚拟密度（pmm）和最高密度物的虚拟密度（pmax）的取值，可依据经验，也可按第 个密度级中和最后一个密度级中物料对密度呈线性分布的假定，分别用以下两个公式推算： Pmn= 2prPl (6) Pmax = 2p, - pa-1 (7 ) 上两式中：p。一最后一个密度级的平均密度，kg/L； 一浮沉试验的最后一个密度，kg/L。 Pr-1- 6.1.6附录C中给出了两个经验模型，可供拟合分配曲线时参考。 6.2可选性曲线 当采用β曲线时，对应于浮物累计产率为零的型值点的横坐标（即最低密度物的灰分）应谨慎推 算。在不影响评定指标计算的情况下，曲线端部可略去不画。 注：如采用M曲线，则端部型值点可以准确地获得。 6.3损失曲线和污染曲线 根据损失曲线和污染曲线的交点，确定等误密度和等误密度下的错配物总量。如果曲线的型值点分 布过于离散，可以仅做数学处理，而不绘制曲线图形。 MT/T145-1997 附录A （提示的附录） 分段连续的三次多项式插值方法 本方法的基本原则是： a）在相邻的两个型值点之间建立一个三次多项式。 b）每个型值点左右两侧的三次多项式在该型值点处具有相同的函数值和一阶导数值。 c）三次多项式在型值点处的导数，指定取为过该型值点及左右相邻型值点所做二次抛物线在该型 值点处的导数值。 实际计算可分为两步进行： a）确定插值曲线在各个型值点处的斜率 设有（N+1)个型值点（r，y)=0,1，，N)，则对于1≤j≤（N1)诸型值点可以做出通过该型 值点及其左右相邻型值点的二次抛物线。抛物线方程的系数p，q.r可由下述方程组中解出： [j-1 +q,1 +r = y,-1 .(A1 ) ++y [++q,++=+ 根据所得到的抛物线方程的系数，就可以定出插值曲线在第」个型值点处的斜率M，即： M,=2pr,+q (A2 ) b）确定每两个相邻型值点之间的插值多项式 根据上文所给出的基本原则，可以对1≤j≤(N一2)诸型值点为起点的插值区间建立三次插值多项 式: y=a+br+cr+d ·(A3) 这个多项式的系数可以根据第，个型值点和第(十1)个型值点的函数值及导数值来确定。多项式 的系数a、b、c、d可由下述方程组解出： [ar,+br,+ca,+d=y, axj+ +brj++ +ct,+1 +d = y,+1 ·(A4) 3az+2br,+c=M, [3azj+1+2bt,+1+c=M+) 附录B （提示的附录） 估计模型参数的阻尼最小二乘法 设有N个型值点（，y）（=1，2，，N)，并准备用数学模型 y=f(r,bo,bi,"bm) (B1 ) 去进行拟合。 式(B1)中，ba，b1，b.可以用6.(1=0，1，，m）表示，它们是非线性模型f的待定参数。这些参数需 要用逐步逼近的方法去估计，现介绍如下： 先给b.一个初始值，记为b，并记初始值与真值（这是未知的)之差为A，于是可以写出： b,=b+A(=0,1,m) ..(B2) 这样，就可以把确定b,的问题化为确定修正量△的问题。为确定△，可将函数(B1)在b附近做台 4 MT/T145-1997 劳级数展开，并略去4,的二次及二次以上的项，得： f(tu,bo,b,...bm) af af 4。+ (B3) abo ab. 式中：f=f(x，bo,,b%) afi. af(r,ba,br,..bm) 'qe ab, (B4 ) bo=b8 6m=b% af 当给定后f和 f(1k,b,,bi,",bm)) Q： af (B5 ) abo ab,