《平均品位计算方法》由会员分享，可在线阅读，更多相关《平均品位计算方法（6页珍藏版）》请在金锄头文库上搜索。
1、平均品位计算新方法张亮（江西贵溪银矿，江西贵溪流3353）摘要:本文采用三次多项式曲线拟合方法,对单位工程品位指标值多边形图进行分段拟合处理，使得多边形图成为一条连续光滑的复合曲线，并对拟合的复合曲线进行积分，求出品位平均值，并给出计算编程计算程序及使用说明。关键词：单工程,三次多项式,拟合，平均品位计算，新方法。1、 引言如图1所示，沿某一直线方向品位的直方图和多边形图，设图1中的直方图间距分别为L、L2、Ln、且L=L1+L2+ Ln（L取样品跨度）,与其相对应的品位指标什值分别为、C2、n.现按常规的加权平均值法计算品位平均值，有又设SiLi，则上式变为 (） P2P1P0 P-1 X从
2、（)式计算公式,可以看出单工程平均品位加权平均值的计算方法实质上是在品位指标的直方图上进行了的，而品位指标值的直方图并不能反映品位指标值的变化情况。由于矿体指标值实际上大多具有连续和光滑的特点，也就是说,指标值曲线形变化,更符合、更接近矿体的实际变化,因此，采用平均值的计算方法计算品位平均值，其计算结果显然不能较好的反映出平均品位指标值的真实情况。为此,欲使计算结果更接近实际结果，本文针对形图为一条连续光滑的曲线，最后运用积分方法求出品位指标值的平均值。2、拟合方法图中，设X为取样点至原点（起点)的距离，Y为其相对应的品位指标什。其中P0、P1、n+1、Pn+2各点的品位均为表外矿品位，点及点
3、为最低可采品位点（边界品位点)。若P0、P1、2点必须可导，具有一阶导数。现将经过P0,1，P，n,+各点的曲线分成PP1，P1P,n1 Pn,PPn+段，每段利用三次多项式曲线来拟合,各分段曲线间光滑连接，使得整个曲线成为一条由多段三次多段三次多项式曲线组成的连续光滑的复合曲线。设每段三次多项式为y=0a1+a2x2+a3x （2)显然,每段三次多项式曲线本身是连续光滑的,欲使各分段曲线间的连接点光滑连接，则曲线间连接点处必须具有一阶导数.现设曲线连接点（如图1所示)处的一阶导数i，根据曲线函数一阶导数的几何意义,i处的导数ti可以用如下公式来拟合确定： ()式中x1、y1、i、yi、xi+
4、、i+1分别是取样点i-1、pi、i+1的坐标。用pi点相邻邻两点的直线斜率之和之半来拟合pi点的一阶导数,这样处理，其拟合曲线的变化趋势与矿体品位指标变化的实际情形基本一致。点i处(i=2，3,n-）导数ti可由（3）式计算可得，i点、n点的导数计算必须助于p0点的坐标和+处的导数计算可借助力n1点的坐标利用(3)式进行计算，同理，点p0和p+1处的导数计算可借助力于p-1点和pn2点的坐标进行计算.上述p1、p0、p1、pn+2各点坐标仅供曲线拟合计算之用，不参与平均品位的计算。3、拟合曲线方程系数求解如图2，以曲线边ii+1为例，由于三次多项式曲线通过pi点和pi+1点,将pi点、pi点
5、的坐标值代入（2)式有y=a+a1x+a22i+3xi （）yi+=a0a11a2x2i1+a3x3i (5)又由于曲线在pi点、pi+处有一阶导数t和i+1,则有 t=l+2a2+3a3x2i (6) t+1a1+22xi+3a3x2i+1 (7）由()、(5)、(6）、(7）解方程组可得 （8） 4、 平均值的计算品位指标值平均也即为图2中复印合曲线函数的平均值,而复y合曲线函数的平均平均等于函数曲线与坐标轴所围成的面积除以函数曲线跨度。如图2所示，复合函数曲线跨度（即取样总长)G=L，P1，。，Pn为边形图转折点；点为0P1段曲线与Y轴的交点，其坐标为H（0,yh）；点为PnPn+1段曲
6、线与取样边界线（边界品位线）GK的交点，交点坐标为K(L，yk)。设OEP1所围成的面积为a,EP1PiPnF所围成的面积为S，PKG所围成的面积为Sc，则复合曲线所围成的总面积为:S=+Sbc .（9)由区边PP所确定的面积为：dX= 。.。（10)利用（0）式可以计算每分段的面积,且(10）式中各分段曲线方程系数am(m0，1，2，3）由(8)式所得。由Pn所围成的总面积Sb为：S= 。.（11）同理可得:s= 。(12） = .(3）由(11）式，(1）式，（13）式可得总面积,故平均品位C为:CS/L 。.。(14)上述解算虽然简单，但计算点较我时，其计算工作量较大，如果采用计算机编程
7、计算，那么上述计算十分方便、快捷。5、计算程序设计计算程序由于BASC语言是最基本且通用性较强的语言,故此，本文采用BASIC语言编制计算程序，程序如下：10 RNT“取样个数，取样长度分别为：”：20 INPT N,L30 DIM X（+），Y(+），T(+1）,A（3）40 FOR I= T0 N+250 READ X(），（1）0 NX I65 FORI=IO N+170 T(I)1/2*（Y(I）-Y（I1)）/(X(I）-X（I）)+（Y(I1） Y（)）X(N+2),Y（N+2）存多边形图中P各点的坐标值T(N1)存多边形图中P0P+1各点的导数值A（3）存每间分段三次多项式曲线系
8、数1，Y1为1点的坐标值(2）运行过程将多边形图中各点（P0Pn+)横纵坐标依次置入DATA语句中,运行该程序，输入N，L入-1点的坐标1,!的值后，转入子程序50语句计算P0P分段三次多项式系数，接着,转入子程序60语句计算HP1分段的面积Sa。自20语入8语句,计算P1K段的面积（b+Sc），最后计算并打印出品位平均值.1、 结束语在进行三次多项式曲线拟合之前，应对特高品位进行处理.采用文中方法计算品位平均值，其结果更符合实际情况，更具有科学性，理论依据充分.而且很适用于计算机编程计算,方便快捷。由于此方法拟合的曲线严格通过了各点，因此在取样的偶然性的确定品位指标值的误差范围内，以及矿体连续较好的单工程平均品位的计算，采用本文计算方法优为方便、实用。参考文献朱晓岚，何新义，陈于恒.矿体几何学，徐州：中国矿业学院出版社,19872同济大学教研室主编,高等数学，北京:高等教育出版社，1985作者：张亮 男 工程师地址：江西省贵溪银矿采掘坑口 邮编：313电话：07073329文中如有不足，请您指教！ /
