ISSN 1009-6248CN 61-1149/P 双月刊

主管单位:中国地质调查局

主办单位:中国地质调查局西安地质调查中心
中国地质学会

    • 中文核心期刊
    • CSCD收录期刊
    • 中国科技核心期刊
    • Scopus收录期刊
高级检索

白云鄂博矿区矿石中磁铁矿及萤石的X射线粉晶衍射半定量分析方法及矿物三维空间分布模型

孙宁岳, 闫国英, 李国武, 王昭静, 申俊峰, 徐渴鑫, 孟文祥, 李凤坤

孙宁岳,闫国英,李国武,等. 白云鄂博矿区矿石中磁铁矿及萤石的X射线粉晶衍射半定量分析方法及矿物三维空间分布模型[J]. 西北地质,2024,57(4):113−120. doi: 10.12401/j.nwg.2024021
引用本文: 孙宁岳,闫国英,李国武,等. 白云鄂博矿区矿石中磁铁矿及萤石的X射线粉晶衍射半定量分析方法及矿物三维空间分布模型[J]. 西北地质,2024,57(4):113−120. doi: 10.12401/j.nwg.2024021
SUN Ningyue,YAN Guoying,LI Guowu,et al. Semi-quantitative Analysis of Magnetite and Fluorite by X-ray Powder Crystal Diffraction and Three-dimensional Distribution Model of Minerals in Bayan Obo Mining Area[J]. Northwestern Geology,2024,57(4):113−120. doi: 10.12401/j.nwg.2024021
Citation: SUN Ningyue,YAN Guoying,LI Guowu,et al. Semi-quantitative Analysis of Magnetite and Fluorite by X-ray Powder Crystal Diffraction and Three-dimensional Distribution Model of Minerals in Bayan Obo Mining Area[J]. Northwestern Geology,2024,57(4):113−120. doi: 10.12401/j.nwg.2024021

白云鄂博矿区矿石中磁铁矿及萤石的X射线粉晶衍射半定量分析方法及矿物三维空间分布模型

基金项目: 包头钢铁(集团)有限责任公司矿山研究院(有限责任公司)2021年A类重点项目“白云鄂博稀土−铌−铁资源矿物标型学研究(BGKY-ZH-2021-Z-017)”资助。
详细信息
    作者简介:

    孙宁岳(1998−),男,博士研究生,研究方向矿物晶体结构与晶体化学。E−mail:sunningyue1998@163.com

    通讯作者:

    李国武(1964−),男,教授,博士,从事矿物晶体结构与晶体化学研究。E−mail:liguowu@cugb.edu.cn

  • 中图分类号: P575.5

Semi-quantitative Analysis of Magnetite and Fluorite by X-ray Powder Crystal Diffraction and Three-dimensional Distribution Model of Minerals in Bayan Obo Mining Area

  • 摘要:

    白云鄂博超大型铁-稀土-铌矿床,伴生大量萤石等资源。现阶段,矿山资源评价采用的是多元素定量分析方法,该方法存在不能客观的表征出可用资源的矿物组成和赋存状态的缺陷。为了较为准确地确定元素赋存形式及磁铁矿、萤石矿物含量,本研究选用粉晶X射线衍射K值法定量分析法,通过以白云鄂博白云石为参比物质分别调配1∶1比例产自白云鄂博的磁铁矿、萤石与白云石单矿物样品,获得以白云石为参比的磁铁矿参比强度KFe3O41 = 0.61、萤石参比强度KCaF21 = 2.51,同时,精选了白云鄂博其他常见共生矿物的K值,以实现对白云鄂博不同矿石类型中矿物含量的半定量分析,通过测试已知标样验证了改进优化K值的正确性和适用性。利用大量实际岩心矿石样品矿物定量分析数据,通过克里金插值法获得白云鄂博主矿磁铁矿和萤石矿物的空间分布三维模型,初步呈现了资源矿物的空间分布特征,推测在矿区深部仍存在巨大的找矿潜力。X射线粉晶衍射定量分析技术直接对可回收资源的矿物半定量分析研究,为矿山精准分采和资源综合利用探索了新的解决思路,也为选矿流程的优化提供的重要的技术参考。

    Abstract:

    Bayan obo super large Fe-RE-Nb deposit, associated with a large amount of fluorite and other resources. At present, the multi-element quantitative analysis method is used in the evaluation of mine resources, which can not objectively characterize the mineral composition and occurrence state of available resources. In order to determine the occurrence form of elements and the content of magnetite and fluorite minerals more accurately, this study uses the powder crystal X-ray diffraction K-value method for quantitative analysis, and uses dolomite as the reference material to prepare 1∶1 proportion of single mineral samples of magnetite, fluorite and dolomite produced in Bayan Obo, the reference intensity of magnetite with dolomite as the reference is KFe3O4 = 0.61, and the reference intensity of fluorite is KCaF2 = 2.51. At the same time, the K value of other common co-minerals in Bayan Obo is selected. In order to realize the semi-quantitative analysis of the mineral content in different ore types of Bayan Obo, the correctness and applicability of the improved optimization K are verified by testing the known standard samples. Based on the quantitative analysis data of a large number of actual core ore samples, a three-dimensional model of spatial distribution of magnetite and fluorite minerals in Bayan Obo main mine was obtained by Kriegin interpolation method, and the spatial distribution characteristics of resource minerals were preliminarily presented. It is speculated that there is still a huge prospecting potential in the deep part of the mine area. X-ray powder diffraction quantitative analysis technology can directly study the semi-quantitative analysis of minerals in recoverable resources, which explores a new solution for precise mining and comprehensive utilization of resources, and also provides an important technical reference for the optimization of mineral processing process.

  • 滑坡是中国分布最多、危害最大的地质灾害。目前对其分布特征、破坏方式、运动机理等方面已有大量的研究(谢婉丽等,20172018a2018b刘颖莹等,2018郭倩怡等,2021Liu et al.,2022)。滑坡往往受多种地质环境作用影响,难以准确预测,使人类生命财产和生存环境受到重大威胁,因而很有必要开展区域性滑坡易发性评价,选择合适预防措施,以减轻地质灾害造成的影响。

    目前,滑坡易发性评价主要包括启发式推断法、数理统计分析法和机器学习法等(王高峰等,2021)。启发式推断法如层次分析、专家打分(冯卫等,2021)等主要凭借主观经验确定因子权重;数理统计分析法如逻辑回归法(屠水云等,2022)、证据权法(杨华阳等,2020)、模糊综合评判法(谢婉丽等,2018b)、信息量法(阮沈勇等,2001)等主要依靠对调查数据进行数理统计分析建立易发性评价模型,可以显著降低评价过程中的主观性,但是其评价结果依赖于大量精确数据支持,对于高维空间的复杂问题,数理统计分析法不可避免的存在欠拟合、预测准确度不高的缺点。

    目前学术界对滑坡易发性评价更偏重于机器学习方法的研究。机器学习可以深挖数据本质,获取最精准的滑坡易发性评价结果,具有自学习、自适应、非线性映射能力强等优点,能充分考虑事件发生概率及其影响因素之间非线性关系(王高峰等,2021),可以很好解决传统概率统计模型的缺陷(Dickson et al., 2016)。常用的评价方法有人工神经网络(冯杭建等,2016田乃满,2020李泽群,2022)、随机森林法(林荣福等,2020刘睿等,2020吴润泽等,2021)、支持向量机(王念秦等,2019王倩等,2021赵铮等,2022)和决策树(杨永刚等,2019)等。但是截至目前,学术界对于最佳的滑坡易发性评价模型仍未有定论。

    为了探究滑坡易发性的最佳模型,笔者将MaxEnt模型引入滑坡易发性评价。$ {\text{MaxEnt}} $模型作为一种机器学习方法,它以最大熵思想为核心,具有精确、高效、样本量要求低、可避免模型过度拟合等优点,常被应用于生态学领域的物种潜在分布范围预测(唐兴港等,2021),其思想与滑坡易发性区划极为类似,都是基于当前已有数据分析研究对象未来的发展变化趋势。近年来,不少国外学者已采用该方法进行滑坡易发性研究(Suchita et al.,2016Maryam et al.,2019Kornejady et al.,2017)。例如,Felicísimo Ángel M等(2013)在西班牙德巴河谷滑坡预测研究中将$ {\text{MaxEnt}} $模型和其他3种机器学习方法进行比较,结果显示$ {\text{MaxEnt}} $模型精度明显高于其他模型。然而,国内基于$ {\text{MaxEnt}} $模型在灾害易发评价方面的研究相对较少(赵冬梅等,2020麦鉴锋等,2021屈新星等,2021)。

    王益区和印台区位于铜川市中部,是铜川市政治、经济、文化和商贸中心,区内沟壑纵横,地下水资源匮乏,地质环境条件脆弱,滑坡分布较多,对区内展开滑坡易发性工作尤为重要。因此,笔者基于铜川市中部地区滑坡分布现状,采用$ {\text{MaxEnt}} $模型,结合$ {\text{ArcGIS}} $空间分析模块,选取高程、坡度、坡向、曲率、距道路的距离、距水系的距离、地形地貌和岩土体类型8个环境因子对滑坡易发性进行研究,为铜川市中部地区防灾减灾与国土空间规划提供科学参考。

    研究区位于铜川市中部,包括王益区、印台区,北接宜君县,南连富平县,东同白水县及蒲城县接壤,西与耀州区毗邻。经纬度涵盖范围为E 108°51′~109°26′,N 34°59′~35°12',总面积为791.74 km2。研究区属半干旱大陆性季风气候,历年年降水量达334.6~879.3 mm,平均为584.5 mm,且多集中在夏季,以暴雨、连续降雨的形式出现。

    研究区位于汾渭地堑与黄土高原过渡区域,丘陵起伏,沟壑密布,属于山、塬、川并存的地貌类型。区内地貌可分为土石山地区、黄土残塬区和河谷阶地区。区内构造简单,未见大型褶皱与断裂。复杂的地貌、脆弱的地质环境,使区域内滑坡灾害极其发育,对研究区社会发展和居民生命与财产安全造成了很大威胁。滑坡主要分布在研究区中部及东部人类工程活动频繁的人口密集区(图1)。

    图  1  研究区地理位置及滑坡分布图
    Figure  1.  Geographical location and landslide distribution map of the study area

    本文所涉及数据包括:$ {\text{DEM}} $(数字高程模型)(下载自地理空间数据云)、岩土体类型、地貌类型(下载自91卫图)、河流、道路信息(下载自全国地理信息资源目录系统)和野外调查收集的滑坡点数(44个)。

    结合野外实际调查数据,研究区共确定滑坡44处,其中约70%以上的滑坡属于浅层牵引式黄土滑坡,65%以上的滑坡属于小型滑坡。研究区滑坡主要分布在道路、水系两侧,且多发生于人类工程活动剧烈的区域。

    根据研究区滑坡分布规律以及诱发因素,结合实际野外调查,经过筛选对比,最终确定高程、坡度、坡向、曲率、距道路的距离、距水系的距离、地形地貌和岩土体类型8个环境因子作为评价指标。

    (1)高程:高程对滑坡易发性影响主要体现在不同高程区域气候、降水、植被类型、人类工程活动均有不同,因此对滑坡产生的影响有所不同。

    (2)坡度:坡度代表了边坡陡倾程度,同时影响了斜坡应力大小与方向。坡度数据在ArcGIS中通过DEM数据提取。

    (3)坡向:坡向对滑坡易发性的影响体现在朝向不同,受到的日照时长与太阳辐射强度不同,阴坡与阳坡的温差较大,同时对降水与植被类型也有一定的影响;坡向数据在ArcGIS中通过DEM数据提取。

    (4)曲率:曲率代表了边坡的凹凸程度,曲率值大于0为凸形坡,等于0为直线形坡,小于0为凹形坡。曲率数据在ArcGIS中通过DEM数据提取。

    (5)距道路的距离:修建道路不可避免会对沿途边坡坡脚进行开挖,破坏边坡的结构,诱发了滑坡的形成。研究区道路密集,通过ArcGIS欧氏距离分析模块获取距道路的距离栅格数据。

    (6)距水系的距离:河水往往对河岸两侧边坡产生冲蚀作用,导致斜坡稳定性降低,往往距离河流越近,滑坡发生的频率越高。通过ArcGIS欧氏距离分析模块获取距水系的距离栅格数据。

    (7)地形地貌:地形地貌是滑坡发生的重要因素,根据收集到的研究区地貌资料,研究区分为黄土残塬、土石山地和河谷阶地3种地貌。

    (8)岩土体类型:研究区岩土体类型共分为坚硬、半坚硬层状碎屑岩,碳酸盐岩,黄土,砂砾石土共5类。研究区整体上呈现为上部为黄土,下部为基岩。研究区由北到南呈现出基岩切割深度逐渐减小的趋势。

    在易发性评价开始前在ArcGIS中预先将所有数据坐标系与栅格像元大小(30 m × 30 m)进行统一,将结果数据以ASC文件导出。本次易发性评价所用的各环境变量图层见图2

    图  2  研究区环境影响因子图
    Figure  2.  Study area environmental impact factor layer

    最大熵($ {\text{MaxEnt}} $)原理是在满足全部已知约束的基础上开展概率事件的合理预测,但是这一过程中对于未知的部分不做任何假设。只有这样,预测结果中的概率分布最为均匀,误判概率最低,因而最终得出的概率分布结果熵值也就最高。换句话说,最大熵理论代表着在满足既定条件的情况下,熵值越大,越接近它的自然现实状态(范亦嵩,2020)。

    $ {\text{MaxEnt}} $模型运行同时需要滑坡地理分布与环境影响因子数据。在本文中,将研究区划分有限像元集,称为X。每个像元xX被分配一个非负概率值px),用来表示发生滑坡的概率分布,像元概率总和为1。这样就构建了研究区滑坡的概率分布,称为$ \pi $

    $ \pi $是未知的,最大熵原理认为熵最大时,概率分布p最接近真实状态。为了求解p,我们需要先构建p分布的约束条件,关于未知的分布$ \pi $,我们已知的信息就是历史滑坡的发生信息和历史滑坡点的环境因子信息,特征函数具有多种表达形式,本次选用MaxEnt模型自带的表达形式,即式(1):

    $$ {f_x}(...) = {\lambda _1}{f_1} + {\lambda _2}{f_2} + {\lambda _3}{f_3} + {\lambda _4}{f_4} + ...... + {\lambda _i}{f_i} $$ (1)

    式中:$ \lambda $是一组参数,$ {f_i} $为像素x位置的第i个环境因子值。已知信息可以表达为所有已发生滑坡的特征函数的平均值,即式(2):

    $$ e = \frac{1}{m}\sum\limits_1^m {{f_{x{\rm{i}}}}(...)} $$ (2)

    式中:m是历史滑坡发生数目,$ {f_{x{\rm{i}}}}(...) $是第i个滑坡所在像元的特征函数。概率分布p的环境因子期望值可按公式(3)计算:

    $$ E={\displaystyle {\sum }_{x\in X}p(\text{x}){f}_{x}}(\mathrm{...}) $$ (3)

    p的环境因子期望值需要无限接近于e,所以分布p的约束条件为式(4):

    $$ |e-E|<\beta $$ (4)

    式中:$ \beta $为任意小的正常数。最大熵分布通常具有式(5)形式:

    $$ p(x)=\frac{e{f}_{x}(\mathrm{...})}{Z} $$ (5)

    式中:$ p(x) $是第$ x $个像元概率值;$ Z $是一个确保所有像元概率值之和为1的常数。通过正则化变换后,约束条件可以表达为使式(6)的值最小:

    $$ \mathrm{ln}{Z}_{\lambda }-\frac{1}{m}{\displaystyle \sum _{{\rm{i}}=1}^{m}{f}_{x{\rm{i}}}(\mathrm{...})}+{\displaystyle \sum _{{\rm{j}}}{\beta }_{\text{j}}\lambda } $$ (6)

    式中:m是历史滑坡数目,$ \lambda $是特征函数的一组参数,$ {\beta _{\text{j}}} $是任意小的常数,只要确定$ \lambda $,使得上式值最小,就可以确定最大熵的概率分布。

    为计算$ \lambda $,MaxEnt模型采用机器学习中的经典算法连续迭代算法,即初始输入$ \lambda $=(0,0,0…0),代表每个环境因子作用相同,根据不断调整$ \lambda $值,使式(6)值不断变小,直到迭代次数达到用户指定次数,或不在显著变小时停止迭代。此时的$ \lambda $值为最终结果中的$ \lambda $值。

    将滑坡点(44个)地理位置和各环境影响因素数据输入MaxEnt模型中,分别取不同比例的滑坡点数量作为训练数据用于建模,余下滑坡点作为测试数据,进行模型验证。模型结果选用Logistic形式,表示某个滑坡在整个模拟区域每个栅格上的存在概率(p)。在预试验中运行多次,使得出的AUC值相对稳定(±0.001)即可。

    选用ROC–AUC值和Kappa一致性检验2种方法对模型精度进行交叉检验。首先,利用模型“Jacknife”中的ROC曲线,对模型模拟的滑坡点潜在分布点预测成效进行评价。其次,通过Kappa值验证模型整体准确度,即通过模型运行结束后的研究区滑坡易发预测图与现状滑坡分布图进行Kappa一致验证。评价指标对应准确度见表1

    表  1  AUC值与Kappa值评价标准表
    Table  1.  Assessment standard of AUC value and kappa value
    精确度极差较差一般较好优秀
    AUC0.5~0.60.6~0.70.7~0.80.8~0.90.9~1
    Kappa0~0.20.2~0.40.4~0.550.55~0.70.7~1
    下载: 导出CSV 
    | 显示表格

    本次研究采用标准差验证模型模拟结果是否稳定,对不同训练比例构建模型运行多次得出的AUC值进行标准差计算,计算公式为:

    $$ \begin{split} \\ {\text{SD}} = \sqrt {\dfrac{{\displaystyle\sum\limits_{i = 1}^N {{{\left( {{X_i} - \bar X} \right)}^2}} }}{{N - 1}}} \end{split}$$ (7)

    式中SD为AUC的标准差,N为重复次数,$ {X_i} $为第i次AUC的大小,$ \bar X $N次的AUC平均值,SD值越小,模型模拟结果越稳定。

    以往研究中,关于模型最适宜训练比例的选取并未达成一致,多数研究者采用训练比例取值范围为70%~90%(Shrestha et al.,2019)。

    为研究不同随机训练比例对MaxEnt构建模型准确度的影响,本次研究选取70%、75%、80%、85%和90%作为随机训练样本点数据比例,按照上述步骤模型运行10次,对评价结果进行标准差计算(表2)。从表中可知,滑坡样本点80%时AUC值达到最大值,但样本点为75%时该模型的标准差最小,表明研究区样本点为75%时MaxEnt模型数据最稳定。因此,笔者采用75%训练数据作为模型运行基础,通过模型10次迭代运行后获得的ROC曲线(图3)。从图中可知,训练集10次模拟得到的AUC值能达到0.905,大于0.9,评价精度优秀,证明该模型可以精确模拟灾害点分布与环境因子之间的关系

    表  2  AUC均值/SD值与训练比例的关系表
    Table  2.  Relationship between AUC mean value/SD value and training proportion
    训练样本比例70%75%80%85%90%
    AUC平均值0.9020.9050.9090.9040.887
    标准差0.07630.06610.08550.08390.0565
    下载: 导出CSV 
    | 显示表格
    图  3  模型运行10次ROC曲线图
    Figure  3.  ROC curve obtained by running the model 10 times

    将模型运行结果提取至ArcGIS,根据Jenks法,将研究区滑坡易发性等级分为5级(图3):低易发(0~0.11)、较低易发(0.11~0.24)、中易发(0.24~>0.43)、较高易发(0.43~0.69)、高易发(0.69~1.00)。易发性等级由高到低,面积分布分别为34.52 km2(4.36%)、45.69 km2(5.77%)、105.16 km2(13.28%)、239.37 km2(30.23%)和367.10 km2(46.36%)。滑坡较高易发区和高易发区主要分布于研究区中部与东部,多位于水系与道路两侧,人口密集分布区域。

    将滑坡现状分布与图4中得出的滑坡潜在分布进行$ {\text{Kappa}} $一致性检验,结果显示$ {\text{Kappa}} $系数为0.76,表示评价分区与现状滑坡点分布十分符合。

    图  4  研究区易发性区划图
    Figure  4.  Study area susceptibility zoning

    采用MaxEnt模型的“Jacknife”模块,可以对每个环境因子进行滑坡易发性贡献率分析(图5),图中深蓝色条块代表无其他环境影响因素干扰下其对易发性分布的贡献率,浅蓝条块即去掉该因子后,其他的所有变量的总贡献率。

    图  5  基于Jacknife环境因子贡献值分析图
    Figure  5.  Analysis of environmental factor contribution based on jacknife

    图5显示,环境变量贡献大小顺序为:距水系的距离>地貌类型>距道路的距离>高程>岩土体类型>坡向>曲率>坡度。可以看出距水系的距离、地貌类型是贡献度排名前二的因子,距道路的距离、高程和岩土体类型也是影响滑坡易发性的重要因子,而坡度、坡向和曲率对滑坡易发性分布的影响较小。

    通过$ {\text{MaxEnt}} $模型所生成的环境因子响应曲线(图6),不仅能反映每个环境影响因素与滑坡易发性之间的关联,也可通过每个环境因子变化时对应的易发性分布概率来体现每个环境因子对最终的易发性评价影响。其中,x轴是环境变量的取值范围,y轴是模型的易发性指数。由于无法用曲线来表示2个分类型变量对易发性评价结果的影响规律,因此用矩形面积来表示岩土体类型、地貌类型对易发性结果的影响,面积越大,对易发性评价影响的能力就越高。

    图  6  环境因子的响应曲线图
    Figure  6.  Response curve of environmental factors

    由高程的响应曲线可以得出(图6a),滑坡易发性于高程850 m左右达到峰值,随后便随高程增加而降低,这是由于研究区居民多居住于海拔较低的河谷平原区,伴随着频繁发生的工程活动;通过坡度的响应曲线来看(图6b),坡度在不断增长时,滑坡易发性随之持续增长;坡向的响应曲线表现为随着坡向取值增大,易发指数呈先增大后减小的趋势(图6c),可能原因为阳坡滑坡发生概率大于阴坡,因为阳坡接受到的阳光和降水更为充沛,岩土体更易风化;从曲率的响应曲线来看(图6d),随着曲率值增大,滑坡易发性随之增大,总体上表现为凸形坡对滑坡易发性的影响大于凹形坡;从距道路的距离的响应曲线来看(图6e),滑坡易发性在距道路800 m范围内时较大,随距道路的距离增加滑坡易发性随之减小,原因就是道路分布的地带由于各种工程活动的破坏,导致环境条件变得恶劣,边坡失稳,灾害频发;距水系的距离曲线与距道路的距离曲线类似,距离水系越近,易发指数越高(图6f);从地貌类型来看(图6g),河谷阶地区对滑坡易发性的影响最大,黄土残塬区次之,土石山地区最小,证明河谷阶地区是最适宜滑坡发育的地貌;从岩土体类型来看(图6h),坚硬层状碎屑岩对易发性影响最低,而砂砾石土对易发性影响最高,这与其易风化、透水性强和工程性质差的性质有关。

    (1)基于地质环境条件选取了高程、坡度、坡向、曲率、距道路的距离、距水系的距离、地形地貌、岩土体类型等8个环境影响因子,建立了以铜川市中部王益区、印台区为研究区的评价指标体系。

    (2)运用$ {\text{ArcGIS}} $$ {\text{MaxEnt}} $模型对研究区滑坡易发性进行了评价。结果表明,高和较高易发区分别占研究区总面积的4.36%和5.77%。滑坡较高易发区和高易发区主要分布于中部和东部道路、水系两侧,是人口集中分布区域。得到的评价结果与实际滑坡分布情况相符,$ {\text{AUC}} $值达到了0.905,表明该模型的评价精度十分优秀,可以用于研究区滑坡易发性评价。

    (3)使用“Jackknife”模块分析环境影响因素对模型结果影响程度,结果显示距水系的距离和地貌类型是贡献度最大的环境因子;坡向、坡度和曲率贡献度最低,表明研究区滑坡对其变化并不敏感。环境影响因素的响应曲线也揭示了易发性随各个环境影响因素的分布规律,表明该模型对分类环境变量和连续环境变量均有较好的适用性。

  • 图  1   配制样品中矿物百分含量与K值法计算结果的线性相关分析(K值来源晶体学数据库)

    Figure  1.   Linear correlation analysis between the percentage of minerals in the prepared samples and the results calculated by the K value method (K-values sourced from the crystallography database)

    图  2   配制样品中矿物百分含量与K值法计算结果线性相关分析(K值来源本次研究实测)

    Figure  2.   Linear correlation analysis between the percentage of minerals in the prepared samples and the results calculated by the K value method (K-values derived from measurements in this study)

    图  3   采样钻孔位置在白云鄂博主矿区分布图 (据徐志豪等,2023修改)

    Figure  3.   Distribution map of sampling locations in the main pit of Bayan Obo

    图  4   钻孔8-01化学分析全铁含量与K值法半定量分析矿物含量协同图

    Figure  4.   Synergistic diagram of full iron content analysis and semi-quantitative analysis of mineral content by K-value method in drill hole 8-01

    图  5   K值法分析与化学元素全铁分析线性相关分析

    Figure  5.   Linear correlation analysis between K-value method analysis and total iron analysis of chemical elements

    图  6   白云鄂博主矿萤石及磁铁矿分布模型(坐标轴字母代表方位S:南 E:东)

    Figure  6.   Distribution map of fluorite and magnetite ore deposits in the main pit of Bayan Obo

    表  1   代表矿石类型样品中典型矿物种类

    Table  1   Typical mineral species represented in samples of different ore types

    勘探
    线号
    钻孔号 样品原
    始编号
    标高(m) 主要矿物种类(%)* 矿石类型
    3副 3F-6 45 1339 磁铁矿(9.98 51.81)、霓石(61.13 26.93)、萤石(6.59 6.85)、白云石(12.63 8.16)、方铅矿(9.98 6.25 霓石型磁铁矿矿石
    4 4-4 50 1364 磁铁矿(17.66 57.80)、萤石(34.81 22.83)、白云石(9.83 4.00)、
    重晶石(37.66 15.35
    萤石型磁铁矿矿石
    8 8-6 106 1210 磁铁矿(3.80 22.86)、萤石(11.33 13.65)、白云石(81.06 60.66)、重晶石(3.79 2.84 白云石型磁铁矿矿石
    9 9-5 30 1397 磁铁矿(18.89 62.08)、萤石(19.02 12.52)、白云石(15.26 6.24)、重晶石(46.84 19.16)、霓石、氟碳铈矿、独居石 萤石型磁铁矿矿石
    8 8-7 23 1416 萤石(36.04 30.33)、磁铁矿(7.29 35.4)、赤铁矿(19.46 11.77)、
    白云石(16.79 10.15)、重晶石(13.04 7.89)、石英(7.38 4.46
    萤石型磁(赤)铁矿矿石
    8 8-01 237 912 磁铁矿(7.28 38.08)、白云石(77.47 50.46)、萤石(3.8 3.98)、
    重晶石(11.84 7.48
    白云石型磁铁矿矿石
    8 8-01 141 1157 磁铁矿(1.79 11.53)、云母(24.8 19.9)、萤石(11.42 18.8)、
    白云石(5.78 4.64)、长石(56.24 45.13
    云母型磁铁矿矿石
    8 8-01 212 977 磁铁矿(21.65 54.73)、钠闪石(28.45 13.21)、萤石(21.78 19.01)、云母(24.38 11.32)、白云石(3.71 1.72 萤石石型磁铁矿矿石
    8 8-6 55 1348 磁铁矿(44.16 70.63)、钠闪石(15.18 4.46)、白云石(4.4 1.29)、
    长石(5.51 1.62)、萤石(9.44 15.76)、云母(21.22 6.23
    萤石石型磁铁矿矿石
    8副 8F-12 110 1218 磁铁矿(56.19 72.98)、萤石(22.3 21.88)、白云石(8.43 2.01)、
    云母(4.70 1.12)、钠闪石(8.37 2.00
    萤石型磁铁矿矿石
    9副 9F-12 120 1193 磁铁矿(11.27 47.92)、钠闪石(42.82 24.84)、萤石(1.72 1.61)、
    云母(15.34 8.9)、白云石(4.04 2.34)、重晶石(24.81 14.39
    钠闪石型磁铁矿矿石
    13 WK13-01 10 1521 磁铁矿(16.47 43.27)、钠闪石(8.09 3.90)、萤石(33.09 32.39)、
    白云石(25.33 12.22)、黄铁矿(13.78 6.65)、云母(3.24 1.56
    萤石石型磁铁矿矿石
      注:*未加粗的数值代表使用相关文献中的萤石与磁铁矿K值进行矿物含量半定量计算,而加粗显示的数值则基于本研究所获得的K值测定。鉴于目前尚未获取适宜于稀土矿物的K值,这些矿物未被纳入本研究的定量分析范畴。因此,含有稀土矿物的样品所得到的矿物半定量分析结果,实为除稀土矿物外其他组分的归一化结果。
    下载: 导出CSV

    表  2   适用于白云鄂博矿区常见矿物的K值参数

    Table  2   Selection parameters of different minerals’ K-values in Bayan Obo district

    矿物种类 hkl 2θ (°) D(nm) K 来源
    白云石 104 30.95 0.2886 2.64 SY/T 5163-2010
    萤石 111 28.28 0.3150 2.51 本次实测
    磁铁矿 113 35.44 0.2530 0.61 本次实测
    石英 100 20.86 0.4255 0.91 SY/T 5163-2010
    方解石 104 29.42 0.3034 2.86 SY/T 5163-2010
    铁白云石 104 30.84 0.2894 2.60 SY/T 5163-2010
    重晶石 211 31.56 0.2833 0.87 SY/T 5163-2010
    黄铁矿 200 33.00 0.2712 2.06 SY/T 5163-2010
    方铅矿 200 30.09 0.2968 4.95 PDF 65-0241
    赤铁矿 104 33.09 0.7000 2.40 PDF 33-0664
    霓(辉)石 −221 31.02 0.2985 0.80 PDF 41-1370
    闪石(族) 110 0.82-0.85 1.36 PDF 89-7282
    云母(族) 001 0.98-1.00 3.4 PDF 42-1437
    长石(族) 0.32-0.33 1.22 俞旭等,1984
      注:SY/T5163-2010来源中华人民共和国行业标准;PDF卡片来源为粉晶衍射数据库(ICDD);2θ (°) 角度对应Cu靶数据。
    下载: 导出CSV
  • 迟广成, 肖刚, 汪寅夫, 等. 铁矿石矿物组分的X射线粉晶衍射半定量分析[J]. 冶金分析, 2015, 35(1): 38−44.

    CHI Guangcheng, XIAO Gang, WANG Yinfu, et al. Semi-quantitative analysis of the mineral components of iron ores by X-ray powder diffraction[J]. Metallurgical Analysis,2015,35(1):38−44.

    付伟, 彭召, 曾祥伟, 等. 基于XRD-Rietveld全谱拟合技术定量分析花岗岩风化壳中矿物组成[J]. 光谱学与光谱分析, 2018, 38(7): 2290−2295.

    FU Wei, PENG Zhao, ZENG Xiangwei, et al. Quantitative Analysis of Mineral Composition in Granite Regolith Based on XRD-Rietveld Full-Spectrum Fitting Method[J]. Spectroscopy and Spectral Analysis,2018,38(7):2290−2295.

    贾建业. 黄铁矿的X射线衍射谱及其找矿意义[J]. 西北地质, 1996(3): 38−45.

    JIA Jianye. X-ray diffraction spectrum of pyrite and its prospecting significance[J]. Northwestern Geology,1996(3):38−45.

    焦景慧, 陈天明. X射线法粉煤灰形成莫来石影响因素的研究[J]. 粉煤灰综合利用, 1996(3): 6−11.

    JIAO Jinghui, CHEN Tianming. Research on the Transformation of Fly Ash into Mullite with X-Ray Diffractometry[J]. Fly Ash Comprehensive Utilization,1996(3):6−11.

    柯昌辉, 李以科, 李立兴, 等. 白云鄂博矿区赋矿“白云岩”地质特征与成因再认识[J]. 中南大学学报(自然科学版), 2021, 52(9): 3047−3063. doi: 10.11817/j.issn.1672-7207.2021.09.007

    KE Changhui, LI Yike, LI Lixing, et al. Petrogenesis of ore-bearing “dolostone” in Bayan Obo deposit, Inner Mongolia, China: insights from geological features[J]. Journal of Central South University (Science and Technology),2021,52(9):3047−3063. doi: 10.11817/j.issn.1672-7207.2021.09.007

    李丹煜, 杨莉, 王金龙, 等. 白云鄂博萤石分布特征及放射性钍元素的影响[J]. 包钢科技, 2020, 46(3): 6−9. doi: 10.3969/j.issn.1009-5438.2020.03.003

    LI Danyu, YANG Li, WANG Jinlong, et al. Distribution Characteristics of Fluorite and Effects of Radioactive Element Thorium in Bayan Obo Deposit[J]. Science and Technology of Baotou Steel,2020,46(3):6−9. doi: 10.3969/j.issn.1009-5438.2020.03.003

    李强, 杨占峰. 白云鄂博主矿各矿石类型稀土配分特征研究[J]. 稀土, 2021, 42(5): 36−42.

    LI Qiang, YANG Zhanfeng. Study on REE Distribution Characteristics of Different Ore Types in Bayan Obo Main Orebody[J]. Chinese Rare Earths,2021,42(5):36−42.

    刘玉龙, 杨刚, 陈江峰, 等. 白云鄂博超大型稀土-铌-铁矿床黄铁矿Re-Os定年[J]. 科学通报, 2005, 50(2): 172−175. doi: 10.3321/j.issn:0023-074X.2005.02.012

    LIU Yulong, YANG Gang, CHEN Jiangfeng, et al. Re-Os Dating of pyrite from Baiyunebo Super Large REE-Nb-Fe Deposit[J]. Chinese Science Bulletin,2005,50(2):172−175. doi: 10.3321/j.issn:0023-074X.2005.02.012

    潘小菲, 张天福, 李岩, 等. 德兴斑岩铜(钼金)矿床蚀变岩石的X粉晶衍射分析及地质意义[J]. 矿床地质, 2012, 31(S1): 335−336.

    PAN Xiaofei, ZHANG Tianfu, LI Yan, el al. X-ray powder diffraction analysis of altered rocks and its geological significance of Tongchang Cu(-Mo-Au) porphyric deposit, Jiangxi Province[J]. Mineral Deposits,2012,31(S1):335−336.

    秦玉芳, 李娜, 王其伟, 等. 白云鄂博选铁尾矿稀土的工艺矿物学研究[J]. 中国稀土学报, 2021, 39(5): 796−804.

    QIN Yufang, LI Na, WANG Qiwei, et al. Technological Mineralogy of Rare Earth in Bayan Obo Iron Tailings[J]. Journal of the Chinese Society of Rare Earths,2021,39(5):796−804.

    王凯怡, 张继恩, 方爱民, 等. 白云鄂博矿床成因——矿体内霓长岩化成矿作用与赋矿白云岩的联系[J]. 岩石学报, 2018, 34(3): 785−798.

    WANG Kaiyi, ZHANG Ji'en, FANG Aimin, et al. Genesis of the Bayan Obo deposit, Inner Mongolia: The f enitized mineralization in the ore bodies and its relation to the ore-bearing dolomite[J]. Acta Petrologica Sinica,2018,34(3):785−798.

    王维维, 李二斗, 金海龙, 等. 白云鄂博萤石型稀土-铁矿石工艺矿物学研究[J]. 有色金属(选矿部分), 2020(6): 14−18. doi: 10.3969/j.issn.1671-9492.2020.06.003

    WANG Weiwei, LI Erdou, JIN Hailong, et al. Study on the Technological Mineralogy of Fluorite Type REE-Fe Ores from Bayan Obo Mine[J]. Nonferrous metals (beneficiation part),2020(6):14−18. doi: 10.3969/j.issn.1671-9492.2020.06.003

    徐志豪, 闫国英, 杨宗锋, 等. 白云鄂博矿床磁铁矿成分标型与深部富铁矿体预测[J]. 地学前缘, 2023, 30(2): 426−439.

    XU Zhihao, YAN Guoying, YANG Zongfeng, et al. Typomorphic characteristics of magnetite and prediction of deep iron-rich orebody in the Bayan Obo ore deposit[J]. Earth Science Frontiers,2023,30(2):426−439.

    薛雍, 江向峰, 钟玉锋. 标准曲线法X射线粉晶衍射直接分析滑石中微量石棉[J]. 岩矿测试, 2010, 29(3): 322−324. doi: 10.3969/j.issn.0254-5357.2010.03.028

    XUE Yong, JIANG Xiangfeng, ZHONG Yufeng. Quantitative Determination of Trace Asbestos inTalcum by X-ray Powder Diffraction with External Standardization[J]. Rock and Mineral Analysis,2010,29(3):322−324. doi: 10.3969/j.issn.0254-5357.2010.03.028

    杨波, 杨莉, 孟文祥, 等. 利用探针片进行X射线粉晶衍射分析在白云鄂博矿床中的应用[J]. 有色金属(选矿部分), 2021(6): 34−42.

    YANG Bo, YANG Li, MENG Wenxiang, et al. Application of X-ray Powder Diffraction Analysis in Bayan Obo Deposit with Microprobe Slice[J]. Nonferrous Metals (Beneficiation Part),2021(6):34−42.

    杨波, 于俊芳, 杨莉, 等. 白云鄂博主矿和东矿不同类型矿石中萤石矿物学特征[J]. 稀土, 2022, 43(1): 90−97.

    YANG Bo, YU Junfang, YANG Li, et al. Mineralogical Characteristics of Fluorite in Different Kinds of Ores of Main and East Orebodies of Bayan Obo Deposit[J]. Chinese Rare Earths,2022,43(1):90−97.

    俞旭, 江超华. 现代海洋沉积矿物及其X射线衍射研究[M].北京: 科学出版社, 1984.

    YU Xu, JANG Chaohua. Modern Marine sedimentary minerals and their X-ray diffraction studies[M]. Beijing:Science Press, 1984.

    于俊芳, 沈茂森, 郭爱芳. 白云鄂博矿萤石浮选分析[J]. 现代矿业, 2018, 593(9): 245−246.

    YU Junfang, SHEN Maosen, GUO Aifang. Flotation Analysis of Fluorite in Baiyunebo Mine[J]. Modern Mining,2018,593(9):245−246.

    于俊芳, 沈茂森, 杨波, 等. 白云鄂博东矿白云石型矿石特征分析[J]. 包钢科技, 2022, 48(3): 10−13.

    YU Junfang, SHEN Maosen, YANG Bo, et al. Analysis on Characteristics of Dolomite Type Orein East Mine of Bayan Obo[J]. Science and Technology of Baotou Steel,2022,48(3):10−13.

    袁珂, 廖立兵, 万红波, 等. 膨润土中方石英和α-石英的定量相分析——X射线衍射外标法和K值法的对比[J]. 硅酸盐学报, 2011, 39(2): 377−382.

    YUAN Ke, LIAO Libing, WAN Hongbo, et al. Quantitative Analysis of Cristobalite and α-Quartz in Bentonite by X-Ray Powder Diffraction-Comparison Between External Standard and K-Value Method[J]. Journal of the Chinese Ceramic Society,2011,39(2):377−382.

    Liu Yulong, Ling Mingxing, Williams I S, et al. The formation of the giant Bayan Obo REE-Nb-Fe deposit, North China, Mesoproterozoic carbonatite and overprinted Paleozoic dolomitization[J]. Ore Geology Reviews,2018,92:73−83. doi: 10.1016/j.oregeorev.2017.11.011

    Yang Kuifeng, Fan Hongrui, Pirajno Franco, et al. The Bayan Obo (China) giant REE accumulation conundrum elucidated by intense magmatic differentiation of carbonatite[J]. Geology,2019,47(12):1198−1202. doi: 10.1130/G46674.1

图(6)  /  表(2)
计量
  • 文章访问数:  309
  • HTML全文浏览量:  15
  • PDF下载量:  71
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-09-17
  • 修回日期:  2024-01-15
  • 录用日期:  2024-01-31
  • 网络出版日期:  2024-04-17
  • 刊出日期:  2024-08-19

目录

/

返回文章
返回