一种基于纯矿物模拟的煤矸X射线识别系统

一种基于纯矿物模拟的煤矸x射线识别系统
技术领域
1.本发明属于煤和矸石智能识别技术领域，具体涉及一种基于纯矿物模拟的煤矸x射线识别系统。


背景技术：

2.随着煤炭开采量的逐年增加，优质品位的煤炭资源几乎被开采殆尽，煤炭资源已经向“贫、细、杂”方向发展。为了提高煤炭质量，目前选煤厂实现煤炭的全粒级洗选。当前主流的煤炭洗选方法分为湿法选煤和干法选煤两种。对于湿法选煤而言，重介质选煤成为主导工艺，其精度较高，但是其缺点也尤为明显：脱介脱水工艺复杂且用水量大，需要后续的煤泥水处理工艺。对于干法选煤而言，主要是空气重介质流化床干法选煤、风选、复合式干选以及射线选煤等，由于无需用水、工艺简单，因此一度成为研究热点。但是，干法选煤分选精度较低且空气污染严重，因此使用范围较小，仅限于西北缺水地区。后来，γ射线、x射线选煤技术开始进入选煤工作者的视线，特别是部分射线选煤设备已经开始在现场进行应用，虽然取得了一定的效果，但是其弊端也尤为明显：首先是设备大都选用γ射线作为射线，虽然穿透力强，但是射源管理是个难题，残余的γ射线会给检修工人的身体带来辐射伤害。二是部分设备为了射源安全和节约成本采用了伪双能x射线，但却会带来“厚度效应”，使得设备对煤和矸石的识别精度大大降低，使得精煤中夹带矸石的比率与矸石中夹带精煤的比率均达不到选煤厂对主选设备的精度需求，很难得到广泛应用和推广。此外，近些年的研究主要集中在煤和矸石的特征选择和算法选择及优化上，其中特征选择包括灰度值、灰度共生矩阵、表面纹理特征、图像熵等，均为基于可见光的特征，而对于不可见光的研究还比较少，亟待解决。


技术实现要素：

3.本发明的目的是克服上述现有技术的不足，提供一种基于纯矿物模拟的煤矸x射线识别系统，其具备高识别精度的优点，能实现对煤和矸石的高效率的稳定识别需求，尤其适用于不可见光基础下的干法选煤行业所使用。
4.为实现上述目的，本发明采用了以下技术方案：
5.一种基于纯矿物模拟的煤矸x射线识别系统，其特征在于包括以下步骤：
6.s1，建立标准煤和矸石的性质参数：
7.对矿区的原煤进行矿物组分的分析，得到包括碳在内的各矿物的含量区间；
8.s2，建立纯矿物灰度模型：
9.利用组成步骤s1所述原煤的各个矿物所对应的纯矿物进行组合，形成作为试验样品的复合矿物；改变复合矿物中不同纯矿物的含量，计算其像素体积的分子数，建立各矿物的分子数也即含量与灰度值之间的数学关系，也即灰度模型；
10.s3，煤和矸石的标准数据采集：
11.选择分选区间内每个粒度级的标准煤矸石样品，提取出每个标准厚度下的煤和矸
石的灰度，找出灰度重叠的范围；
12.s4，使用灰度模型进行分选：
13.灰度重叠范围内的煤和矸石使用s2步骤所述的灰度模型进行分选，也即：将待分选的煤和矸石的灰度输入灰度模型，反推得到待分选的煤和矸石内各矿物的分子数；
14.s5，结果的判断：
15.若步骤s4的数据中的碳的分子数占全部分子数的比值数落入矿区煤的碳含量范围或大于该范围，则判断为煤，否则为矸石。
16.优选的，所述矿物还包括高岭土、蒙脱石和石英。
17.优选的，所述步骤s2包括以下子步骤：
18.(1)将碳粉、高岭土粉、石英粉和蒙脱石粉均压制成厚度不同，但密度和质量均相等的长方体块；
19.(2)改变厚度，获得不同厚度下的复合矿物；该复合矿物通过将x厚的高岭土、y厚度的石墨、z厚度的石英与h厚度的蒙脱石叠放在一起形成；
20.(3)提取所有复合矿物的灰度值；
21.(4)计算每种复合矿物的一个像素体积的分子数；
22.(5)建立分子数与灰度之间的数学模型，得到灰度模型。
23.优选的，所述x，y，z、h均属于7到11之间的正整数，厚度的变化为整数变化。
24.优选的，子步骤(3)中，提取所有复合矿物的灰度值之前，需对图像进行中值滤波降噪处理，再对图像进行自适应图像分割，从而得到完整的矿物区域；再提取矿物区域的平均灰度值，减去背景的灰度值后，得到每种矿物的绝对灰度值，并将此作为复合矿物的灰度值。
25.优选的，子步骤(4)的操作方法如下：
26.首先计算出复合矿物中总的分子数量，再利用成像后图像中矿物区域所占总区域像素数的比例，将样品划分为一个个的像素体积，再以下式计算出每一个像素体积中的分子数代替物质的总分子数：
[0027][0028]
其中：
[0029]
n
e
‑
p
为样品的像素体积的分子数；
[0030]
m为样品的质量；
[0031]
n
a
为阿伏伽德罗常数；
[0032]
m为物质的量；
[0033]
s0为图像的面积；
[0034]
s为矿物区域的面积；
[0035]
α为图片的像素。
[0036]
优选的，在得到每种复合矿物的一个像素体积的分子数后，再以分子数为自变量，灰度值为因变量，进行数据拟合，并以拟合度最高为条件得到拟合方式，获得复合矿物中不同蒙脱石厚度下的拟合方程y＝f(m,g,k,q)，其中y为复合矿物的灰度值，g为复合矿物中石墨的分子数，m为复合矿物中蒙脱石的分子数，k为复合矿物中高岭土的分子数，q为复合矿
物中石英的分子数。
[0037]
优选的，所述步骤s3包括：
[0038]
在矿区选择标准厚度的煤和矸石样品，测得每一个粒度级标准厚度样品的灰度值用来划分灰度重叠的区域，对于区域之外的灰度值可以直接判断为煤和矸石，在区域之内的灰度值需要带入灰度模型之中，具体为：
[0039]
在对标准厚度的煤和矸石进行灰度提取后，得到煤和矸石的灰度，对于煤和矸石来说，总是有灰度范围相交的区域，无法区分是煤还是矸石，而此区域称作灰度重叠区域；
[0040]
提取待分选的煤和矸石的灰度并以矩阵形式处理，矩阵中最大的元素为图像中最大的灰度值，该值在灰度重叠区域外则可清楚的分辨出煤和矸石，在灰度重叠区域中则无法判断，只能得到是hx厚度的煤或者是hy厚度的矸石。
[0041]
优选的，将在灰度重叠区间的矿物，统计其灰度矩阵中的众数，将众数带入所述模型y＝f(m,g,k,q)中得到如下数据，
[0042][0043]
数据输出为矩阵，然后对矩阵作如下运算：
[0044][0045]
对矩阵中每一行做如下处理：
[0046][0047][0048][0049][0050]
最终结果为：
[0051]
[0052]
该矩阵的行表示的是在该灰度值下，蒙脱石，石英，高岭土和碳的含量的组合为a1,b1,c1,d1，接下来根据该矿区各矿物含量中蒙脱石，高岭土和石英的含量区间和该矩阵进行对比，从全部组合中选出满足要求的矿物含量矩阵，也即要求该矩阵中a，b，c的值与测得的矿物含量相同，然后再计算所选出的矩阵的d的平均值d
n
，该平均值代表了在假设的蒙脱石，高岭土，石英含量一定的情况下碳的含量，如果d
n
在煤的标准碳含量区间内或者大于该区间则为煤，否则为矸石。
[0053]
本发明的有益效果在于：
[0054]
在描述效果之前，此处先说明煤夹矸难以分选的原因，在于：(1)从表面的纹理特征来看，表面可能是煤而内部是矸石；(2)煤夹矸在x射线下成像后的灰度值是与纯煤或者纯矸石的灰度值极为相近难以识别。而本发明所确保煤夹矸能选的原因在于：煤夹矸内部的矿物组分不同，可能是碳含量居多，可能是其他矿物含量居多，如果能确定夹矸煤中的碳含量占比范围即能分选。传统射线分选方法是以煤和矸石的x射线吸收系数不同或者结合可见光来提取煤和矸石的表面纹理差异来识别煤和矸石，两者对纯煤和纯矸的分选效果显著，但是对煤夹矸分选效果极差。
[0055]
有鉴于此，本发明从煤和矸石的内部差异出发，即矿物组分的差异，从而抛弃了传统的基于可见光的特征算法，转而提供了一种基于不可见光的算法优化流程。也即本发明首先利用组成煤和矸石矿物组分所对应的纯矿物，形成作为试验样品的复合矿物，再改变复合矿物中不同矿物的含量，计算其像素体积的分子数，随之建立各矿物分子数与灰度值之间的数学关系，再将煤和矸石的灰度值带入模型，获得各矿物组分的分子数，作为判断是否为煤或者矸石的依据。具体而言：本发明通过对煤和矸石的矿物组分的模拟，得到煤和矸石中主要的四种矿物的分子数为自变量，灰度值为因变量的方程，通过该方程可以对无法用灰度值识别的夹矸煤进行处理，将灰度值带入方程，可以得到每个矿物成分的含量进而识别。
[0056]
显然，通过上述方式，本发明具备了以下好处：
[0057]
(1)本发明提供的识别方法从分子层面进行出发，使用纯矿物模拟煤和矸石中复杂的矿物组成，用试验的方法混合矿物中碳，高岭土，蒙脱石和石英四种矿物在不同含量下的灰度值模型，通过此模型可以得到煤夹矸中各矿物组分的相对分子数量，在进一步得到各矿物的相对含量，对比所选矿区标准煤样的矿物成分即可确定煤和矸石，即可确定煤和矸石。上述模拟矿物组成的方法避免了分选过程中外部因素对煤和矸石的影响，比如矸石表面的煤灰，煤和矸石在振动时的分裂而导致的厚度减小从而使煤和矸石的x射线吸收系数改变等，从而能得到精确度更高的灰度模型，能显著提高夹矸煤的识别率，解决了煤夹矸这一问题。
[0058]
(2)本发明的识别方法简单，识别精确度高，在实现过程中不需要复杂的算法和庞大的煤和矸石图片的数据库，只需要普通的工业计算机即可实现，从而能够增大代码运行速度，提高分选效率。
[0059]
(3)本发明原理简单，不需要复杂的辅助设备，也能有效的降低了设备和维护成本。
附图说明
[0060]
图1是本发明的工作流程图；
[0061]
图2是纯矿物的厚度与灰度之间的关系曲线图；
[0062]
图3是纯矿物像素体积分子数与灰度之间的关系曲线图；
[0063]
图4和图5为实施例1的两组灰度值图。
具体实施方式
[0064]
本发明的具体结构如图1
‑
3所示，其主要包括以下步骤：
[0065]
步骤1，建立标准煤和矸石的性质参数；
[0066]
对矿区进行煤质分析，得到该矿区煤和矸石的碳、高岭土、蒙脱石、石英的含量范围，记为标准煤和矸石的性质参数。
[0067]
步骤2，建立纯矿物灰度模型；
[0068]
其具体步骤为：
[0069]
步骤201：将碳粉也即石墨粉，高岭土粉，石英粉和蒙脱石粉压制成7mm
‑
11mm，以厚度为变化量，且密度和质量均相等的长方体块。
[0070]
步骤202：将x厚的高岭土，y厚度的石墨，z厚度的石英与7mm
‑
11mm的蒙脱石叠放在一起，其中x，y，z属于7到11之间的正整数。
[0071]
步骤203：提取步所有样品的灰度值；
[0072]
提取矿物灰度值之前对图像进行中值滤波降噪处理，在对图像进行自适应图像分割，得到完整的矿物区域，再提取矿物区域的平均灰度值，减去背景的灰度值后，得到每种矿物的绝对灰度值，并将此作为矿物的灰度值，样品厚度与灰度之间的关系如图2。
[0073]
步骤204：计算每种样品的一个像素体积的分子数；
[0074]
计算样品的像素体积的分子数，其具体方法如下：
[0075]
首先计算出样品中总的分子数量，再利用成像后图像中矿物区域所占总区域像素数的比例，将样品划分为一个个像素体积，再计算出每一个像素体积中的分子数，公式如下：
[0076][0077]
其中：
[0078]
n
e
‑
p
为样品的像素体积的分子数；
[0079]
m为样品的质量；
[0080]
n
a
为阿伏伽德罗常数；
[0081]
m为物质的量；
[0082]
s0为图像的面积；
[0083]
s为矿物区域的面积；
[0084]
α为图片的像素；
[0085]
计算出的单矿物像素体积的分子数与灰度之间的关系如图3。
[0086]
步骤205：建立灰度与分子数之间的数学模型；
[0087]
得到每种样品的像素体积的分子数之后，以分子数为自变量，灰度值为因变量，进
行数据拟合，并选择拟合度最高得到拟合方式，得拟合方程y＝f(m,g，k，q)。其中，y为复合矿物的灰度值，m为复合矿物中蒙脱石的分子数，g为复合矿物中石墨的分子数，k为复合矿物中高岭土的分子数，q为复合矿物中石英的分子数。显然，理论上，实际操作时，也可以通过复合矿物中的任何非碳的矿物进行厚度变化来获得上述模型，也能进行相应计算并获得相应的煤和矸石的在线分辨效果。
[0088]
步骤3，煤和矸石的标准数据采集；
[0089]
在矿区选择标准粒度的煤和矸石样品，测得每一个粒度级标准样品的灰度值用来划分灰度重叠的区域，对于区域之外的灰度值可以直接判断为煤和矸石，在区域之内的灰度值需要带入模型之中。
[0090]
步骤4，使用灰度模型进行分选。
[0091]
对于模型的使用具体分以下步骤：
[0092]
步骤401：图像灰度提取；
[0093]
步骤402：对灰度值进行初步分类；
[0094]
在对标准厚度的煤和矸石进行灰度提取后，得到煤和矸石的灰度方程y＝f(m，g，k，q)，自变量为各矿物的分子数，对于煤和矸石来说，此方程有不同的值域，而值域重叠的部分为灰度重叠区域；
[0095]
再提取待分选的煤和矸石的灰度并以矩阵形式处理，矩阵中最大的元素为图像中最大的灰度值，该值在灰度重叠区域外则可清楚的分辨出煤和矸石，在区域中则带入灰度方程，得到该灰度值下的矿物是a厚度的煤或者是b厚度的矸石。
[0096]
步骤403：将分类后的灰度值带入模型，具体如下：
[0097]
步骤404：对每一个灰度得到的结果进行统计，对比，分类。
[0098]
实际操作时，需将在灰度重叠区间的矿物，统计其灰度矩阵中的众数，将众数带入所述模型y＝f(m,g,k,q)中得到如下数据，
[0099][0100]
数据输出为矩阵，然后对矩阵作如下运算：
[0101][0102]
对矩阵中每一行做如下处理：
[0103]
[0104][0105][0106][0107]
最终结果为：
[0108][0109]
该矩阵的行表示的是在该灰度值下，蒙脱石，石英，高岭土和碳的含量的组合为a1,b1,c1,d1，接下来根据该矿区各矿物含量中蒙脱石，高岭土和石英的含量区间和该矩阵进行对比，从全部组合中选出满足要求的矿物含量矩阵，也即要求该矩阵中a，b，c的值与测得的矿物含量相同，然后再计算所选出的矩阵的d的平均值d
n
，该平均值代表了在假设的蒙脱石，高岭土，石英含量一定的情况下碳的含量，如果d
n
在煤的标准碳含量区间内护着大于该区间则为煤，否则为矸石。
[0110]
为进一步的描述本发明，此处结合图4
‑
5，给出以下实施例1：
[0111]
实施例1
[0112]
蒙脱石，石英，高岭土和碳的厚度为7mm
‑
11mm的5个厚度级，将四种矿物按照不同厚度的组合可得到5
×5×5×
5＝625个灰度数据，利用多元线性回归拟合方法得到四元线性回归方程：
[0113][0114]
经过拟合优度检验，f检验和t检验后表明模型符合要求。
[0115]
以淮南矿区的煤样做验证，经过对淮南矿区的矿物成分分析，得到煤中高岭土含量为12％
‑
16％，石英含量为9％
‑
14％，碳含量为65％
‑
75％，蒙脱石含量为5％
‑
10％。其10mm
‑
50mm标准煤样的灰度范围是＞50，标准矸石范围在＜70，所以其灰度重合范围是[50,70]。选择10mm
‑
50mm的200个煤和矸石，其中煤150个，矸石50个，其灰度值如图4所示。
[0116]
将灰度重叠范围内的灰度值输入到灰度模型中，以灰度值62为例，提取得到结果为一个216
×
4的矩阵；该矩阵的第一列为蒙脱石含量5％
‑
10％，6个量级；石英含量9％
‑
14％，6个量级；高岭土含量5％
‑
10％，6个量级；也即均满足前述a,b,c的选择要求。矩阵第四列d为碳含量，需判定。此时，计算216
×
4的矩阵中的全部第四列碳元素的平均值d
n
＝82.36，大于标准煤的碳含量范围，所以判断为煤，与实际复合。对其他在灰度重叠区间的夹矸煤做同样处理得到结果如图5所示。
[0117]
综上，该方法的灰度值取的是矿物灰度图像中矿物区域灰度矩阵中所有元素的众数，而灰度重叠区域是按照最大元素划分，操作简便快捷。
[0118]
当然，对于本领域技术人员而言，本发明不限于上述示范性实施例的细节，而还包括在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现的相同或类似
方法。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。
[0119]
此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。
[0120]
本发明未详细描述的技术部分均为公知技术。