CN107367480B - 基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法 - Google Patents

基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法 Download PDF

Info

Publication number
CN107367480B
CN107367480B CN201710852763.2A CN201710852763A CN107367480B CN 107367480 B CN107367480 B CN 107367480B CN 201710852763 A CN201710852763 A CN 201710852763A CN 107367480 B CN107367480 B CN 107367480B
Authority
CN
China
Prior art keywords
sample
thermal infrared
infrared spectrum
point
content
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710852763.2A
Other languages
English (en)
Other versions
CN107367480A (zh
Inventor
熊宏启
孙厚广
潘鹏飞
徐冬林
毛亚纯
刘善军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ansteel Mining Co Ltd
Original Assignee
Ansteel Mining Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Ansteel Mining Co Ltd filed Critical Ansteel Mining Co Ltd
Priority to CN201710852763.2A priority Critical patent/CN107367480B/zh
Publication of CN107367480A publication Critical patent/CN107367480A/zh
Application granted granted Critical
Publication of CN107367480B publication Critical patent/CN107367480B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/3563Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing solids; Preparation of samples therefor

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明针对于现有技术中化学分析法测定鞍山式铁矿中二氧化硅含量存在的问题,提供了一种基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法,属于矿山地质与遥感应用技术领域。该方法包括采集矿样、制备样品、样品光谱测定、构建光谱吸收深度(SAD)指标,建立SAD与SiO2含量的线性拟合模型。该方法基于样品的热红外光谱特征,提出了一种新的指标——光谱吸收深度SAD,用于SiO2含量反演,实现了矿山现场铁矿石SiO2含量实时原位测定,为矿山采场矿体边界实时快速圈定与精准区划奠定基础。

Description

基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法
技术领域
本发明属于矿山地质与遥感应用技术领域,特别是涉及一种基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法。
背景技术
鞍山式铁矿在我国分布最广,资源开采量最大,而矿山确定铁矿石样品的品位是地质勘探以及采矿工程中很重要的一个环节,因此品位确定方法对于矿山的生产效率是非常重要的。由于这种矿石中主要含有3种矿物成分,即石英(化学式SiO2)、磁铁矿(Fe3O4)、赤铁矿(Fe2O3),要确定矿石的铁含量,必须首先要确定SiO2的含量。
目前,鞍山式铁矿的成分确定方法主要为化学分析法,这种方法程序多、过程复杂、化验周期较长,使得矿石品位测试结果相对配矿流程存在滞后效应,这将导致选厂指标波动,使得配矿质量和效果也会受到极大影响。同时,矿体测试样本数量也会受到限制,必然导致采场取样密度偏低,难以满足采场精细区划的要求。并且,化学分析法不能完成矿体品位的原位即时测定,因此其难以满足露天采场矿体实时精准圈定与区划的技术要求。
在地表岩石成分定量反演方面,国内外学者利用热红外光谱的方法对岩矿SiO2含量反演进行了大量研究。Copper、闫柏琨以光谱库为数据源研究了硅酸盐矿物SiO2含量与光谱特征的关系,并建立了发射率光谱克里斯琴森特征(CF)、透射特征(TF)、CF+TF特征与岩石SiO2含量指数的统计关系。Greenhagen B T以克里斯琴森特征(CF)处的光谱作为组成矿物的权重,探测月球的主要岩石。杨长保等通过建立SiO2含量的多元回归分析模型对岩石SiO2含量进行定量反演,反演结果与地质图和实地考察结果较相符。杨杭等通过建立基于9.38、11.18、12.36和12.82um波段的归一化SiO2光谱指数,实现了对SiO2含量的高精度提取。郭帮杰、曾江维等利用TASI数据中余辉带特征(Reststrahlen Features简称RF)的B5、B6、B12、B17四个波段进行了岩石中石英含量定量分析。但是上述研究多是基于样品光谱的10~14um波段进行的研究,且上述研究建立的反演预测模型多是基于多波段的多元回归模型,模型较繁琐。
并且,上述研究均是利用光谱特征对自然界岩石和沙子中的SiO2含量的反演研究,但未见对于鞍山式变质型铁矿中SiO2含量的反演研究,因此,亟需对鞍山式变质型铁矿矿石现场实时原位测试方法。
发明内容
本发明针对于现有技术中化学分析法测定鞍山式铁矿中二氧化硅含量存在的问题,提供了一种基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法。该方法基于样品的热红外光谱特征,提出了一种新的指标——光谱吸收深度SAD,用于SiO2含量反演,实现了矿山现场铁矿石SiO2含量实时原位测定,为矿山采场矿体边界实时快速圈定与精准区划奠定基础。
本发明的技术方案为,一种基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法,包括如下步骤:
1)在鞍山式铁矿所在地,采集矿样,处理成圆形薄片状样品;
较好的,圆形薄片状样品直径为6-10cm、厚度为0.8-1.5cm;
进一步的,当所述矿样内部分布有石英条带时,制作的样品的表面与内部分布的石英条带相垂直;
当所述矿样内部没有分布石英条带,而分布有铁矿条带时,制作的样品的表面与内部分布的铁矿条带相垂直;
当所述矿样内部即没有分布石英条带、也没有分布铁矿条带时,则直接处理成圆形薄片状样品;
2)利用热红外光谱辐射计对样品的两个表面分别进行光谱测定,取两个表面的热红外光谱曲线均值作为该样品的热红外光谱曲线;
较好的,所述样品光谱测定方法为:
辐射计的光谱测试范围为2.55~16um,光谱分辨率优于4cm-1,下行辐射的测试选用金板进行;将扫描镜头对准金板,镜头与金板成90°,距离在0.5-1m左右,获得下行辐射;然后将扫描镜头垂直对准样品表面,距离在0.5-1米,视场角以覆盖样品测试表面为准,采集光谱曲线;
3)构建光谱吸收深度(SAD)指标:
SAD=EB-EB
式中:SAD为光谱吸收深度,EB为该样品的热红外光谱曲线中B点发射率值,EB,为B’点发射率值;
其中,所述热红外光谱曲线中B点为波长8.55~8.70um间曲线的峰顶,其发射率值为EB,所对应的波长为B点波长;
热红外光谱曲线中点A为在热红外光谱曲线中B点左侧的第一个谷底,热红外光谱曲线中点C为在热红外光谱曲线中B点右侧的第一个谷底;将点A与点C连成一条直线LAC,在直线LAC上于B点波长处对应的点即为B’点,其发射率值为EB’,EB与EB’的差值即为该峰的吸收深度(SAD);
4)建立SAD与SiO2含量的线性拟合模型,作为SiO2含量的反演预测模型,其公式为:
SiO2=332.68·SAD+34.779;
通过该模型计算得到SiO2含量。
本发明的原理为:根据热红外光谱理论,石英矿物(化学式SiO2)在8~10um波长处出现双吸收谷,该特征也被称为石英矿物在热红外波段的余晖带特征,它是石英矿物的诊断性特征。根据样品的光谱和化学成分测试结果,可以发现铁矿石中SiO2含量越高,该吸收谷就越深,SiO2含量与吸收谷深度具有很高的相关性。本发明通过样品的热红外光谱特征,根据公式计算即可反演样品中的SiO2含量。
与现有技术相比,本发明的优势在于:
1、本发明的方法,能够实时、快速、经济的确定鞍山式铁矿中二氧化硅的含量,且平均误差小于6%,精度较高。
2、本方法实现了矿山现场铁矿石SiO2含量实时原位测定,为矿山采场矿体边界实时快速圈定与精准区划奠定基础。
附图说明
图1、不同SiO2含量铁矿石样品的热红外光谱。
图2、SAD指数建立示意图。
具体实施方式
实施例1
1、测试样品采集
在鞍山式铁矿石采场采集鞍山式铁矿石样品,然后进行钻孔、取芯、切割、磨平处理,形成直径6-10cm、厚度为1cm的圆形薄片状实验样品。
为保证光谱测试样品表面的成分与整个样品的一致性,要求在加工的过程中,当所述矿样内部分布有石英条带时,制作的圆形薄片样品的表面与内部分布的石英条带相垂直;
当所述矿样内部没有分布石英条带,而分布有铁矿条带时,制作的圆形薄片样品的表面与内部分布的铁矿条带相垂直;
当所述矿样内部即没有分布石英条带、也没有分布铁矿条带时,则直接处理成圆形薄片状样品。
2、测试设备
1)使用便携式红外光谱辐射计,对样品进行光谱测定,辐射计光谱测试范围为2.55~16um,光谱分辨率优于4cm-1
2)视场角以覆盖样品测试表面为准。
3)下行辐射的测试应选用漫反射好的金板进行。
3、测试方法
1)在进行目标地物波谱测量前,首先要进行金板波谱的采集,获得下行辐射,将扫描镜头对准金板,镜头与金板成90°,距离在0.5-1m左右。
2)金板测定完毕后对样本红外辐射进行测试,采集波谱曲线。测试时将扫描镜头垂直对准样品表面,距离在0.5-1米左右,要求测试范围基本覆盖样品的表面。
4、把测试过光谱的实验样品粉碎、研磨成粉,然后进行化学分析,确定其中SiO2成分的含量。
5、根据热红外光谱理论,石英矿物(化学式SiO2)在8~10um波长处出现双吸收谷,该特征也被称为石英矿物在热红外波段的余晖带特征,它是石英矿物的诊断性特征,见图1不同SiO2含量铁矿石样品的热红外光谱。对照样品的光谱和化学成分测试结果,可以发现铁矿石中SiO2含量越高,该吸收谷就越深,SiO2含量与吸收谷深度具有很高的相关性。
6、根据样品的热红外光谱曲线,计算于8~10um波长处的光谱吸收深度SAD指标,SAD计算公式为:
SAD=EB-EB',
如图2的SAD指数建立示意图所示,在光谱8~10um波长之间取8.62um处的峰顶作为点B,其发射率为EB,于峰顶的左侧第一个谷底即双吸收谷的左谷部位8.48um处取点A作为左肩部,其发射率为EA;于峰顶的右侧第一个谷底即双吸收谷的右谷部位8.75um处取点C作为右肩部,其发射率为EC。左肩部点A与右肩部点C连成一条直线LAC,在直线LAC上于8.62um处对应的点为点B’,其发射率为EB'。EB与EB'的差值即为该峰的吸收深度。
7、在确定SAD后,根据下面的反演预测模型来确定SiO2含量:
SiO2=332.68·SAD+34.779。
通过10组样品的化学实测和本发明反演计算的结果可以看出,二者对SiO2测定误差小于6%,具体见表1。
表1、样品的SiO2测定数据表
样品编号 实测SiO<sub>2</sub>(%) 反演SiO<sub>2</sub>(%) 差值(%)
1 53.39 54.02 0.63
2 49.96 54.02 4.06
3 58.85 52.07 -6.78
4 66.03 67.82 1.79
5 61.12 57.14 -3.98
6 53.28 55.59 2.31
7 54.07 53.73 -0.34
8 55.29 53.00 -2.29
9 72.62 66.98 -5.64
10 55.07 55.13 0.06

Claims (4)

1.一种基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法,其特征在于,包括如下步骤:
1)在鞍山式铁矿所在地,采集矿样,处理成圆形薄片状样品;
2)利用热红外光谱辐射计对样品的两个表面分别进行光谱测定,取两个表面的热红外光谱曲线均值作为该样品的热红外光谱曲线;
3) 计算光谱吸收深度指标:
SAD= EB- EB’
式中:SAD为光谱吸收深度,EB为该样品的热红外光谱曲线中B点发射率值, EB’为B’点发射率值;
其中,所述B点为热红外光谱曲线中波长8.55~8.70um间曲线的峰顶,其发射率值为EB,所对应的波长为B点波长;
所述热红外光谱曲线中B’点为:热红外光谱曲线中点A为在热红外光谱曲线中B点左侧的第一个谷底,热红外光谱曲线中点C为在热红外光谱曲线中B点右侧的第一个谷底;将所述点A与点C连成一条直线LAC,在直线LAC上于B点波长处对应的点即为B’点,其发射率值为E B’
4) 再根据下述的反演预测模型计算SiO2含量:
SiO 2=332.68•SAD+34.779。
2.根据权利要求1所述的一种基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法,其特征在于,所述圆形薄片状样品直径为6-10cm、厚度为0.8-1.5cm。
3.根据权利要求1或2所述的一种基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法,其特征在于,当所述矿样内部分布有石英条带时,制作的样品的表面与内部分布的石英条带相垂直;
当所述矿样内部没有分布石英条带,而分布有铁矿条带时,制作的样品的表面与内部分布的铁矿条带相垂直。
4.根据权利要求1所述的一种基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法,其特征在于,所述样品光谱测定方法为:
辐射计的光谱测定范围为2.55~16um,光谱分辨率优于4cm-1,下行辐射的测试选用金板;将扫描镜头对准金板,镜头与金板成90°,距离0.5-1m,获得下行辐射;然后将扫描镜头垂直对准样品表面,距离在0.5-1m,视场角以覆盖样品测试表面为准,采集光谱曲线。
CN201710852763.2A 2017-09-20 2017-09-20 基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法 Active CN107367480B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710852763.2A CN107367480B (zh) 2017-09-20 2017-09-20 基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710852763.2A CN107367480B (zh) 2017-09-20 2017-09-20 基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法

Publications (2)

Publication Number Publication Date
CN107367480A CN107367480A (zh) 2017-11-21
CN107367480B true CN107367480B (zh) 2020-07-28

Family

ID=60302956

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710852763.2A Active CN107367480B (zh) 2017-09-20 2017-09-20 基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法

Country Status (1)

Country Link
CN (1) CN107367480B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110161075B (zh) * 2019-05-05 2020-09-22 中国地质科学院矿产资源研究所 一种利用热红外波谱的石榴子石分析方法及其应用
CN110927087A (zh) * 2019-12-23 2020-03-27 鞍钢集团矿业有限公司 基于可见光-近红外光谱的鞍山式铁矿磁性率反演方法
CN110907379A (zh) * 2019-12-23 2020-03-24 鞍钢集团矿业有限公司 一种基于随机森林算法的鞍山式铁矿成分及含量测定方法
CN111691888B (zh) * 2020-06-12 2021-10-22 鞍钢集团矿业有限公司 一种露天矿山精确铲装系统及方法
CN113049521B (zh) * 2021-03-29 2023-09-26 自然资源实物地质资料中心 识别碳酸盐岩的方法、装置及设备、存储介质
CN113125353A (zh) * 2021-04-16 2021-07-16 自然资源实物地质资料中心 检测方解石含量的方法和预设模型的建设方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0869341A1 (en) * 1997-03-31 1998-10-07 Nec Corporation Thermal infrared detecting device
CN101520343A (zh) * 2009-03-25 2009-09-02 中国科学院上海技术物理研究所 热红外光谱成像系统装校装置及装校方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102879107B (zh) * 2012-09-25 2015-02-25 中国科学院地理科学与资源研究所 热红外下行辐射测定装置及其测定方法
CN106198437B (zh) * 2016-06-21 2018-10-26 中国土地勘测规划院 一种植被光合有效辐射吸收率的遥感估算方法及装置
CN106770054A (zh) * 2016-12-30 2017-05-31 黑龙江禾禾遥感科技有限公司 基于光谱吸收特征的黑土含水量速测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0869341A1 (en) * 1997-03-31 1998-10-07 Nec Corporation Thermal infrared detecting device
CN101520343A (zh) * 2009-03-25 2009-09-02 中国科学院上海技术物理研究所 热红外光谱成像系统装校装置及装校方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
沙质土壤热红外高光谱特征及其含沙量预测研究;黄启厅;《光谱学与光谱分析》;20110831;第31卷(第8期);第2195-2199页 *

Also Published As

Publication number Publication date
CN107367480A (zh) 2017-11-21

Similar Documents

Publication Publication Date Title
CN107367480B (zh) 基于热红外光谱的鞍山式铁矿中二氧化硅含量测定方法
CN107589094B (zh) 基于光谱特征的鞍山式铁矿石类型确定方法
CN104612675B (zh) 一种碳酸盐岩地层随钻岩性快速识别方法
McBride et al. Estimating forest soil quality from terrain measurements of apparent electrical conductivity
Dang et al. Shale gas potential of Lower Permian marine-continental transitional black shales in the Southern North China Basin, central China: Characterization of organic geochemistry
CN110333200B (zh) 一种基于短波红外光谱圈定矿化中心方法
Pommier Interpretation of magnetotelluric results using laboratory measurements
Gauthier et al. Assessing XRF for the geochemical characterization of radiolarian chert artifacts from northeastern North America
CN104749341B (zh) 碳酸盐岩类烃源岩有效性的判断方法
CN110967311B (zh) 基于红外光谱和磁化率测量的斑岩矿床蚀变分带识别方法
CN104142519B (zh) 一种泥岩裂缝油藏预测方法
CN106932837A (zh) 一种测试晶质铀矿形成年龄的方法
CN105445441A (zh) 储层含气饱和度确定方法及装置
CN104849409A (zh) 一种泥页岩成岩演化模拟实验方法
CN114646682B (zh) 一种基于绿帘石微量元素的找矿方法
Camargo et al. Pedotransfer functions to assess adsorbed phosphate using iron oxide content and magnetic susceptibility in an Oxisol
CN111999260B (zh) 一种热红外光谱识别含锂辉石伟晶岩的方法及热红外光谱的应用
CN103954679B (zh) 含铀脉体多点取样等时线定年方法
CN101581222A (zh) 一种石油钻井液核磁共振随钻分析方法
CN110132902A (zh) 一种下古生界海相页岩有机质成熟度的评价方法
CN110320571A (zh) 一种致密砂岩储层岩石脆性测井评价方法
Pitblado et al. Pilot study experiments sourcing quartzite, Gunnison Basin, Colorado
CN109709130B (zh) 一种测试全油基钻井液地层含油率的方法
Sharma et al. Multi-technique characterization of shale reservoir quality parameters
Fan et al. Geochemical characteristics of chlorite in the Luohe iron deposit in the middle-lower Yangtze metallogenic Belt, Eastern China

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant