CN109207722B - 一种确定离子型稀土浸矿过程源汇项的方法 - Google Patents
一种确定离子型稀土浸矿过程源汇项的方法 Download PDFInfo
- Publication number
- CN109207722B CN109207722B CN201811397711.1A CN201811397711A CN109207722B CN 109207722 B CN109207722 B CN 109207722B CN 201811397711 A CN201811397711 A CN 201811397711A CN 109207722 B CN109207722 B CN 109207722B
- Authority
- CN
- China
- Prior art keywords
- rare earth
- ions
- leaching
- concentration
- relation
- 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
Links
- 229910052761 rare earth metal Inorganic materials 0.000 title claims abstract description 404
- 150000002910 rare earth metals Chemical class 0.000 title claims abstract description 213
- 238000002386 leaching Methods 0.000 title claims abstract description 171
- 238000000034 method Methods 0.000 title claims abstract description 86
- BFNBIHQBYMNNAN-UHFFFAOYSA-N ammonium sulfate Chemical compound N.N.OS(O)(=O)=O BFNBIHQBYMNNAN-UHFFFAOYSA-N 0.000 claims abstract description 45
- 229910052921 ammonium sulfate Inorganic materials 0.000 claims abstract description 45
- 235000011130 ammonium sulphate Nutrition 0.000 claims abstract description 45
- 238000012360 testing method Methods 0.000 claims abstract description 42
- 238000005342 ion exchange Methods 0.000 claims abstract description 41
- -1 rare earth ions Chemical class 0.000 claims description 298
- 239000002245 particle Substances 0.000 claims description 90
- 230000014509 gene expression Effects 0.000 claims description 63
- 238000009792 diffusion process Methods 0.000 claims description 40
- 239000003795 chemical substances by application Substances 0.000 claims description 23
- 239000007791 liquid phase Substances 0.000 claims description 16
- 239000007788 liquid Substances 0.000 claims description 13
- 238000009795 derivation Methods 0.000 claims description 12
- 150000002500 ions Chemical class 0.000 claims description 12
- 238000011160 research Methods 0.000 claims description 10
- KCXVZYZYPLLWCC-UHFFFAOYSA-N EDTA Chemical compound OC(=O)CN(CC(O)=O)CCN(CC(O)=O)CC(O)=O KCXVZYZYPLLWCC-UHFFFAOYSA-N 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000004448 titration Methods 0.000 claims description 7
- 238000005119 centrifugation Methods 0.000 claims description 6
- 230000035484 reaction time Effects 0.000 claims description 6
- 238000000926 separation method Methods 0.000 claims description 6
- 239000007787 solid Substances 0.000 claims description 6
- 229910001404 rare earth metal oxide Inorganic materials 0.000 claims description 4
- 238000001179 sorption measurement Methods 0.000 claims description 4
- QGZKDVFQNNGYKY-UHFFFAOYSA-O Ammonium Chemical compound [NH4+] QGZKDVFQNNGYKY-UHFFFAOYSA-O 0.000 claims description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 2
- 239000011707 mineral Substances 0.000 claims description 2
- 150000003254 radicals Chemical class 0.000 claims 2
- 150000005837 radical ions Chemical class 0.000 claims 1
- 238000011065 in-situ storage Methods 0.000 abstract description 6
- 239000006185 dispersion Substances 0.000 abstract description 4
- 238000004088 simulation Methods 0.000 abstract 1
- 239000000243 solution Substances 0.000 description 69
- 238000013508 migration Methods 0.000 description 7
- 230000005012 migration Effects 0.000 description 7
- 238000002347 injection Methods 0.000 description 5
- 239000007924 injection Substances 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 239000002689 soil Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000005065 mining Methods 0.000 description 3
- 230000035515 penetration Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000013312 flour Nutrition 0.000 description 1
- 239000010413 mother solution Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000006228 supernatant Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- C—CHEMISTRY; METALLURGY
- C22—METALLURGY; FERROUS OR NON-FERROUS ALLOYS; TREATMENT OF ALLOYS OR NON-FERROUS METALS
- C22B—PRODUCTION AND REFINING OF METALS; PRETREATMENT OF RAW MATERIALS
- C22B3/00—Extraction of metal compounds from ores or concentrates by wet processes
- C22B3/04—Extraction of metal compounds from ores or concentrates by wet processes by leaching
- C22B3/12—Extraction of metal compounds from ores or concentrates by wet processes by leaching in inorganic alkaline solutions
- C22B3/14—Extraction of metal compounds from ores or concentrates by wet processes by leaching in inorganic alkaline solutions containing ammonia or ammonium salts
-
- C—CHEMISTRY; METALLURGY
- C22—METALLURGY; FERROUS OR NON-FERROUS ALLOYS; TREATMENT OF ALLOYS OR NON-FERROUS METALS
- C22B—PRODUCTION AND REFINING OF METALS; PRETREATMENT OF RAW MATERIALS
- C22B3/00—Extraction of metal compounds from ores or concentrates by wet processes
- C22B3/04—Extraction of metal compounds from ores or concentrates by wet processes by leaching
-
- C—CHEMISTRY; METALLURGY
- C22—METALLURGY; FERROUS OR NON-FERROUS ALLOYS; TREATMENT OF ALLOYS OR NON-FERROUS METALS
- C22B—PRODUCTION AND REFINING OF METALS; PRETREATMENT OF RAW MATERIALS
- C22B59/00—Obtaining rare earth metals
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P10/00—Technologies related to metal processing
- Y02P10/20—Recycling
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Metallurgy (AREA)
- Organic Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Manufacturing & Machinery (AREA)
- Materials Engineering (AREA)
- Mechanical Engineering (AREA)
- Geochemistry & Mineralogy (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Inorganic Chemistry (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Manufacture And Refinement Of Metals (AREA)
Abstract
本发明涉及一种确定离子型稀土浸矿过程源汇项的方法,适用于离子型稀土矿浸矿过程的模拟和最优浸矿剂的浓度的确定。本发明包含四个步骤:(1)确定稀土矿样的离子交换选择系数;(2)确定稀土矿样的稀土品位;(3)建立浸矿过程的源汇项模型;(4)确定源汇项模型中的参数。本发明结合对流‑弥散方程,可以模拟离子型稀土浸矿过程,可以确定最优浸矿剂的浓度,为原地浸矿确定浸矿剂的浓度提供依据。当采用最优浓度为12.0g/L的硫酸铵溶液进行柱浸试验,得到的稀土浸取率高达96.3%。
Description
技术领域
本发明属于原地浸矿领域,涉及一种确定离子型稀土浸矿过程源汇项的方法。
背景技术
南方离子型稀土富含中重稀土资源,中重稀土对高端科技产品的发展起着至关重要的作用。现阶段主要以硫酸铵为浸矿剂,采用原地浸矿的工艺开采离子型稀土矿。原地浸矿的工艺为在矿区按孔间距为1.0m~2.0m布置注液孔网,根据探矿资料确定硫酸铵单耗,往注液孔内注入一定浓度的硫酸铵溶液,溶液中的铵根离子与吸附在稀土矿样颗粒表面的稀土离子发生离子交换,解吸下来的稀土离子随溶液流入收液工程,待硫酸铵溶液注完后,向注液孔内注入上清液,待母液中的稀土离子小于0.01g/L后,完成开采(参见期刊《岩土工程学报》中文章“离子型稀土原地浸矿注液孔周含水率分布的计算模型”,2018第5期)。影响离子型稀土矿开采的主要因素为注液孔网的布置、硫酸铵溶液浓度、硫酸铵用量和收液工程的布置,其中硫酸铵溶液的浓度和硫酸铵用量是影响稀土浸取率的主要因素。
现在确定硫酸铵溶液的浓度和硫酸铵用量主要采用两种方式(见中国专利“离子型稀土原地浸矿的浸矿剂浓度计算方法”,申请号201810130892.5),一种是根据前期开采的经验,结合探矿资料,确定硫酸铵溶液的浓度和硫酸铵用量,由于地质条件的复杂性,各个矿山的地质条件有明显差异,通过经验确定的硫酸铵浓度往往不在合理的浓度范围内,若浓度过高,会使周边水体的氨氮超标,造成环境污染,若浓度过低,稀土浸出不充分,造成资源浪费,而且增加浸矿周期;另一种是现场取样进行室内柱浸试验,根据室内柱浸试验的结果为实际开采的硫酸铵溶液的浓度和硫酸铵用量提供参考(参见期刊《科学技术与工程》中文章“浸矿剂浓度对稀土浸取效果的影响”,2015第28期),由于实际矿山的地质条件复杂、稀土品位分布不均匀,室内柱浸试验难以较为准确的模拟实际矿山的地质条件和品位分布,因而室内柱浸试验的结果对准确确定实际浸矿过程的硫酸铵溶液的浓度和硫酸铵用量的指导意义很有限。
建立合适的数学模型,通过探矿资料获取实际矿山的水运动参数、溶质运移参数和稀土品位等的分布信息,采用数值计算的方法模拟浸矿过程是充分考虑实际矿山的地质条件和稀土品位分布,确定硫酸铵溶液的浓度和硫酸铵用量最直接有效的方法。离子型稀土浸矿过程可以分为外部溶液的铵根离子交换吸附在稀土矿样颗粒表面稀土离子的浸出过程和稀土离子与铵根离子随溶液向下流动的运移过程,由此可见,离子型稀土浸矿过程实质是铵根离子和稀土离子在矿体内的溶质运移过程(浸出过程为源汇项,运移过程主要受对流和弥散的影响),科学研究铵根离子和稀土离子在矿体内的溶质运移过程,建立相应的溶质运移模型、模拟浸矿过程,是从科学的角度确定硫酸铵溶液的浓度和硫酸铵用量。关于离子在矿体内的运移过程,已经做了大量研究,运移的机理和模型已相对成熟,而浸出过程的研究很少,提出确定离子型稀土浸矿过程源汇项的方法对科学确定硫酸铵溶液的浓度和硫酸铵用量、绿色开采离子型稀土矿具有重要意义。
发明内容
本发明的目的是提供一种确定离子型稀土浸矿过程源汇项的方法。
本发明的技术方案:
一种确定离子型稀土浸矿过程源汇项的方法,包括以下步骤:
第一步,确定稀土矿样的离子交换选择系数:
进行浸矿剂浓度为1.0g/L~20.0g/L浸取离子型稀土矿样的平衡杯浸试验,采用现有方法测试浸出液中稀土离子浓度,得到浸矿剂溶液与浸出液中稀土离子浓度的关系曲线,以离子交换选择系数为基本未知量,采用离子交换模型拟合浸矿剂浓度与浸出液中稀土离子浓度的关系曲线,确定稀土矿样的离子交换选择系数;
第二步,确定稀土矿样的稀土品位:
进行液固比为5:1mL/g、硫酸铵溶液的浓度为20g/L的平衡杯浸试验,对浸出液和稀土矿样进行固液分离,量取与杯浸试验相同体积的硫酸铵溶液对稀土矿样进行淋洗,测试浸出稀土离子的总质量,采用关系式(1)计算稀土矿样的品位:
关系式(1)中:η为稀土品位,mRE为浸出的稀土离子的总质量,ms为稀土矿样的质量,MREO和MRE分别为稀土氧化物和稀土离子的摩尔质量;
第三步,建立浸矿过程的源汇项模型:
离子型稀土矿样的稀土浸出包括三个过程:(1)铵根离子通过扩散层达到稀土矿样颗粒表面,(2)铵根离子与吸附在稀土矿样颗粒表面的稀土离子发生离子交换反应,(3)解吸下来的稀土离子通过扩散层进入外部溶液;
铵根离子和稀土离子通过扩散层的过程采用Fick定律描述,可以分别通过关系式(2)和关系式(3)进行计算;
关系式(2)和关系式(3)中:NRE为吸附在稀土矿样颗粒表面的稀土离子的摩尔数,t为时间,d为求导符号,ms为稀土矿样的质量,和分别为外部溶液中铵根离子和稀土离子的浓度,和分别为靠近稀土矿样颗粒表面的液相铵根离子和液相稀土离子的浓度,和αRE分别为与铵根离子和稀土离子在扩散层中的扩散系数相关的参数,和MRE分别为铵根离子和稀土离子的摩尔质量;
铵根离子与吸附在稀土矿样颗粒表面的稀土离子发生离子交换反应采用关系式(4)描述;
浸出过程中,以稀土矿样颗粒、扩散层和外部溶液所组成的体系为研究对象,在整个体系内,铵根离子和稀土离子应满足质量守恒,分别采用关系式(5)和关系式(6)描述;
关系式(5)和关系式(6)中:ms为稀土矿样的质量,VL为外部溶液体积,和分别为吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,和分别为外部溶液中铵根离子和稀土离子的浓度,和分别为浸矿前吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,和分别为浸矿前外部溶液中铵根离子和稀土离子的浓度;
以稀土矿样颗粒为研究对象,浸矿前后,稀土矿样颗粒表面的吸附位数量不发生变化,则浸矿前后吸附在稀土矿样颗粒表面的铵根离子和稀土离子应满足关系式(7);
把关系式(2)、关系式(3)和关系式(5)~关系式(7)代入关系式(4)中,可以得到关系式(8);
关系式(8)中:K为离子交换选择系数,和MRE分别为铵根离子和稀土离子的摩尔质量,和分别为浸矿前吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,和分别为浸矿前外部溶液中铵根离子和稀土离子的浓度,为外部溶液中稀土离子的浓度,为吸附在稀土矿样颗粒表面的稀土离子的浓度,和αRE分别为与铵根离子和稀土离子在扩散层中的扩散系数相关的参数,t为时间,d为求导符号,ms为稀土矿样的质量,VL为外部溶液体积;
关系式(8)中为吸附在稀土矿样颗粒表面的稀土离子浓度随时间的变化速率,即为浸出过程的源汇项,以为基本未知量,关系式(8)为关于的一元三次方程,采用现有的方法解关系式(8),即可得到离子型稀土浸矿过程的源汇项,结果为关系式(9),关系式(9)中的参数可以通过关系式(10)~关系式(15)确定;
关系式(9)~关系式(15)中:p、q、a、b、c和e为计算过程的临时参数,为吸附在稀土矿样颗粒表面的稀土离子的浓度,和MRE分别为铵根离子和稀土离子的摩尔质量,和αRE分别为与铵根离子和稀土离子在扩散层中的扩散系数相关的参数,和分别为浸矿前外部溶液中铵根离子和稀土离子的浓度,为外部溶液中稀土离子的浓度,和分别为浸矿前吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,K为离子交换选择系数,t为时间,d为求导符号,ms为稀土矿样的质量,VL为外部溶液体积;
第四步,确定源汇项模型中的参数:
进行稀土矿样的浸出动力学试验,测试不同反应时间浸出液中稀土离子的浓度,确定外部溶液中稀土离子的生成速率,通过关系式(16)得到稀土矿样颗粒表面吸附的稀土离子浓度随时间的变化速率的试验数据,
本发明基于Fick定律和离子交换模型,全面考虑浸出离子型稀土矿的物理过程,并给出了源汇项中参数的测试方法,通过探矿资料确定实际矿山的地形,以及水运动参数、溶质运移参数和稀土品位的分布信息,采用数值方法对实际矿山进行网格划分,通过随机有限元的方法将水运动参数、溶质运移参数和稀土品位的分布信息赋值给各个网格节点,本发明提出的源汇项,结合溶质运移模型,可以在充分考虑实际矿山复杂的地质条件和稀土品位分布不均匀的基础上,模拟离子型稀土浸矿过程,并结合优化算法,确定最优浸矿剂(硫酸铵溶液)的浓度。从实施效果可以看出,当采用最优浓度为12.0g/L的硫酸铵溶液进行柱浸试验,得到的稀土浸取率高达96.3%。
具体实施方式
应用本发明在实验室内进行柱浸试验,以硫酸铵作为浸矿剂,先通过平衡杯浸试验确定稀土矿样的离子交换选择系数和品位,根据离子型稀土矿浸出过程的物理过程建立浸出过程的数学模型,通过浸出动力学试验确定模型中的参数,具体的实施方案如下。
第一步,确定稀土矿样的离子交换选择系数:
以硫酸铵作为浸矿剂,取10份稀土矿样,每份20.0g,分别采用1.0g/L、1.5g/L、2.0g/L、3.0g/L、5.0g/L、7.0g/L、10.0g/L、13.0g/L、16.0g/L和20.0g/L的硫酸铵溶液按液固比为5:1mL/g进行浸矿,先振荡2.0h,后静置0.5h,采用中速滤纸进行固液分离,采用EDTA滴定法测试浸出液中的稀土离子浓度,得到浸矿剂溶液与浸出液中稀土离子浓度的关系曲线,以离子交换选择系数为基本未知量,采用离子交换模型拟合浸矿剂浓度与浸出液中稀土离子浓度的关系曲线,确定稀土矿样的离子交换选择系数为12.59×10-10L2/g2;
第二步,确定稀土矿样的稀土品位:
以硫酸铵作为浸矿剂,取20.0g稀土矿样,采用20.0g/L的硫酸铵溶液按液固比为5:1mL/g进行浸矿,先振荡2.0h,后静置0.5h,采用中速滤纸进行固液分离,再采用体积为100mL、浓度为20g/L的硫酸铵溶液淋洗尾矿,采用EDTA滴定法测试浸出稀土离子的总质量,采用关系式(1)计算稀土矿样的品位为1.01‰;
关系式(1):
关系式(1)中:η为稀土品位,mRE为浸出的稀土离子的总质量,ms为稀土矿样的质量,MREO和MRE分别为稀土氧化物和稀土离子的摩尔质量;
第三步,建立浸矿过程的源汇项模型:
离子型稀土矿样的稀土浸出包括三个过程:(1)铵根离子通过扩散层达到稀土矿样颗粒表面,(2)铵根离子与吸附在稀土矿样颗粒表面的稀土离子发生离子交换反应,(3)解吸下来的稀土离子通过扩散层进入外部溶液;
铵根离子和稀土离子通过扩散层的过程采用Fick定律描述,可以分别通过关系式(2)和关系式(3)进行计算;
关系式(2)和关系式(3):
关系式(2)和关系式(3)中:NRE为吸附在稀土矿样颗粒表面的稀土离子的摩尔数,t为时间,d为求导符号,ms为稀土矿样的质量,和分别为外部溶液中铵根离子和稀土离子的浓度,和分别为靠近稀土矿样颗粒表面的液相铵根离子和液相稀土离子的浓度,和αRE分别为与铵根离子和稀土离子在扩散层中的扩散系数相关的参数,和MRE分别为铵根离子和稀土离子的摩尔质量;
铵根离子与吸附在稀土矿样颗粒表面的稀土离子发生离子交换反应采用关系式(4)描述;
关系式(4):
浸出过程中,以稀土矿样颗粒、扩散层和外部溶液所组成的体系为研究对象,在整个体系内,铵根离子和稀土离子应满足质量守恒,分别采用关系式(5)和关系式(6)描述;
关系式(5)和关系式(6):
关系式(5)和关系式(6)中:ms为稀土矿样的质量,VL为外部溶液体积,和分别为吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,和分别为外部溶液中铵根离子和稀土离子的浓度,和分别为浸矿前吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,和分别为浸矿前外部溶液中铵根离子和稀土离子的浓度;
以稀土矿样颗粒为研究对象,浸矿前后,稀土矿样颗粒表面的吸附位数量不发生变化,则浸矿前后吸附在稀土矿样颗粒表面的铵根离子和稀土离子应满足关系式(7);
关系式(7):
把关系式(2)、关系式(3)和关系式(5)~关系式(7)代入关系式(4)中,可以得到关系式(8);
关系式(8):
关系式(8)中:K为离子交换选择系数,和MRE分别为铵根离子和稀土离子的摩尔质量,和分别为浸矿前吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,和分别为浸矿前外部溶液中铵根离子和稀土离子的浓度,为外部溶液中稀土离子的浓度,为吸附在稀土矿样颗粒表面的稀土离子的浓度,和αRE分别为与铵根离子和稀土离子在扩散层中的扩散系数相关的参数,t为时间,d为求导符号,ms为稀土矿样的质量,VL为外部溶液体积;
关系式(8)中为吸附在稀土矿样颗粒表面的稀土离子浓度随时间的变化速率,即为浸出过程的源汇项,以为基本未知量,关系式(8)为关于的一元三次方程,采用卡丹公式解关系式(8),即可得到离子型稀土浸矿过程的源汇项,结果为关系式(9),关系式(9)中的参数可以通过关系式(10)~关系式(15)确定;
关系式(9)~关系式(15):
关系式(9)~关系式(15)中:p、q、a、b、c和e为计算过程的临时参数,为吸附在稀土矿样颗粒表面的稀土离子的浓度,和MRE分别为铵根离子和稀土离子的摩尔质量,和αRE分别为与铵根离子和稀土离子在扩散层中的扩散系数相关的参数,和分别为浸矿前外部溶液中铵根离子和稀土离子的浓度,为外部溶液中稀土离子的浓度,和分别为浸矿前吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,K为离子交换选择系数,t为时间,d为求导符号,ms为稀土矿样的质量,VL为外部溶液体积;
第四步,确定源汇项模型中的参数:
取20.0g聚合稀土矿样装入250mL的离心瓶中,将离心瓶放置于磁力搅拌器上,倒入体积为100mL、浓度为20.0g/L的硫酸铵溶液,同时启动磁力搅拌器,转速为1000r/min;并计时,到达预定反应时间:0s、5s、10s、15s、20s、30s、60s、90s、150s后,迅速将离心盒置于离心机中进行离心,转速设置为3000r/min,离心时间为20s,待离心结束后,采用EDTA滴定法测试不同反应时间浸出液中稀土离子的浓度,确定外部溶液中稀土离子的生成速率,通过关系式(16)得到稀土矿样颗粒表面吸附的稀土离子浓度随时间的变化速率的试验数据,
关系式(16):
以和αRE为基本未知量,采用第三步建立的源汇项模型中关系式(9)拟合稀土矿样颗粒表面吸附的稀土离子浓度随时间的变化速率的试验数据,最终确定稀土浸矿过程的源汇项参数和αRE分别为3.7×103s/L和1.3×103s/L。
实施效果:
以硫酸铵作为浸矿剂,进行定水头的柱浸试验,柱浸的孔隙比设置为1.0,水头高度设置为5.0cm,柱长为1.0m,采用自动收液器每隔2.0h收集一次浸出液,采用EDTA滴定法测试浸出液中稀土离子的浓度,待浸出液中稀土离子小于0.01g/L后结束试验,得到稀土离子的穿透曲线。把源汇项代入对流-弥散方程中,模拟离子型稀土矿样的浸矿过程,计算穿透曲线与试验穿透曲线的相关系数达到0.91,改变硫酸铵溶液的浓度,采用本发明提出的源汇项的确定方法和结合对流-弥散方程,计算不同硫酸铵溶液的浓度对应的稀土浸取率,设定目标浸取率为95.0%,当计算浸取率大于目标浸取率时,确定最优硫酸铵溶液的浓度为12.0g/L,当采用浓度为12.0g/L的硫酸铵溶液进行柱浸试验,得到稀土浸取率为96.3%。
Claims (2)
1.一种确定离子型稀土浸矿过程源汇项的方法,其特征是,包括以下步骤:
第一步,确定稀土矿样的离子交换选择系数:
进行浸矿剂浓度为1.0g/L~20.0g/L浸取离子型稀土矿样的平衡杯浸试验,采用现有方法测试浸出液中稀土离子浓度,得到浸矿剂溶液与浸出液中稀土离子浓度的关系曲线,以离子交换选择系数为基本未知量,采用离子交换模型拟合浸矿剂浓度与浸出液中稀土离子浓度的关系曲线,确定稀土矿样的离子交换选择系数;
第二步,确定稀土矿样的稀土品位:
进行液固比为5:1mL/g、硫酸铵溶液的浓度为20g/L的平衡杯浸试验,对浸出液和稀土矿样进行固液分离,量取与杯浸试验相同体积的硫酸铵溶液对稀土矿样进行淋洗,测试浸出稀土离子的总质量,采用关系式(1)计算稀土矿样的品位:
关系式(1)中:η为稀土品位,mRE为浸出的稀土离子的总质量,ms为稀土矿样的质量,MREO和MRE分别为稀土氧化物和稀土离子的摩尔质量;
第三步,建立浸矿过程的源汇项模型:
离子型稀土矿样的稀土浸出包括三个过程:(1)铵根离子通过扩散层达到稀土矿样颗粒表面,(2)铵根离子与吸附在稀土矿样颗粒表面的稀土离子发生离子交换反应,(3)解吸下来的稀土离子通过扩散层进入外部溶液;
铵根离子和稀土离子通过扩散层的过程采用Fick定律描述,可以分别通过关系式(2)和关系式(3)进行计算;
关系式(2)和关系式(3)中:NRE为吸附在稀土矿样颗粒表面的稀土离子的摩尔数,t为时间,d为求导符号,ms为稀土矿样的质量,和分别为外部溶液中铵根离子和稀土离子的浓度,和分别为靠近稀土矿样颗粒表面的液相铵根离子和液相稀土离子的浓度,和αRE分别为与铵根离子和稀土离子在扩散层中的扩散系数相关的参数,和MRE分别为铵根离子和稀土离子的摩尔质量;
铵根离子与吸附在稀土矿样颗粒表面的稀土离子发生离子交换反应采用关系式(4)描述;
浸出过程中,以稀土矿样颗粒、扩散层和外部溶液所组成的体系为研究对象,在整个体系内,铵根离子和稀土离子应满足质量守恒,分别采用关系式(5)和关系式(6)描述;
关系式(5)和关系式(6)中:ms为稀土矿样的质量,VL为外部溶液体积,和分别为吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,和分别为外部溶液中铵根离子和稀土离子的浓度,和分别为浸矿前吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,和分别为浸矿前外部溶液中铵根离子和稀土离子的浓度;
以稀土矿样颗粒为研究对象,浸矿前后,稀土矿样颗粒表面的吸附位数量不发生变化,则浸矿前后吸附在稀土矿样颗粒表面的铵根离子和稀土离子应满足关系式(7);
把关系式(2)、关系式(3)和关系式(5)~关系式(7)代入关系式(4)中,可以得到关系式(8);
关系式(8)中:K为离子交换选择系数,和MRE分别为铵根离子和稀土离子的摩尔质量,和分别为浸矿前吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,和分别为浸矿前外部溶液中铵根离子和稀土离子的浓度,为外部溶液中稀土离子的浓度,为吸附在稀土矿样颗粒表面的稀土离子的浓度,和αRE分别为与铵根离子和稀土离子在扩散层中的扩散系数相关的参数,t为时间,d为求导符号,ms为稀土矿样的质量,VL为外部溶液体积;
关系式(8)中为吸附在稀土矿样颗粒表面的稀土离子浓度随时间的变化速率,即为浸出过程的源汇项,以为基本未知量,关系式(8)为关于的一元三次方程,采用现有的方法解关系式(8),即可得到离子型稀土浸矿过程的源汇项,结果为关系式(9),关系式(9)中的参数可以通过关系式(10)~关系式(15)确定;
关系式(9)~关系式(15)中:p、q、a、b、c和e为计算过程的临时参数,为吸附在稀土矿样颗粒表面的稀土离子的浓度,和MRE分别为铵根离子和稀土离子的摩尔质量,和αRE分别为与铵根离子和稀土离子在扩散层中的扩散系数相关的参数,和分别为浸矿前外部溶液中铵根离子和稀土离子的浓度,为外部溶液中稀土离子的浓度,和分别为浸矿前吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,K为离子交换选择系数,t为时间,d为求导符号,ms为稀土矿样的质量,VL为外部溶液体积;
第四步,确定源汇项模型中的参数:
进行稀土矿样的浸出动力学试验,测试不同反应时间浸出液中稀土离子的浓度,确定外部溶液中稀土离子的生成速率,通过关系式(16)得到稀土矿样颗粒表面吸附的稀土离子浓度随时间的变化速率的试验数据,
2.根据权利要求1所述的一种确定离子型稀土浸矿过程源汇项的方法,其特征是,具体实施为:
第一步,确定稀土矿样的离子交换选择系数:
以硫酸铵作为浸矿剂,取10份稀土矿样,每份20.0g,分别采用1.0g/L、1.5g/L、2.0g/L、3.0g/L、5.0g/L、7.0g/L、10.0g/L、13.0g/L、16.0g/L和20.0g/L的硫酸铵溶液按液固比为5:1mL/g进行浸矿,先振荡2.0h,后静置0.5h,采用中速滤纸进行固液分离,采用EDTA滴定法测试浸出液中的稀土离子浓度,得到浸矿剂溶液与浸出液中稀土离子浓度的关系曲线,以离子交换选择系数为基本未知量,采用离子交换模型拟合浸矿剂浓度与浸出液中稀土离子浓度的关系曲线,确定稀土矿样的离子交换选择系数为12.59×10-10L2/g2;
第二步,确定稀土矿样的稀土品位:
以硫酸铵作为浸矿剂,取20.0g稀土矿样,采用20.0g/L的硫酸铵溶液按液固比为5:1mL/g进行浸矿,先振荡2.0h,后静置0.5h,采用中速滤纸进行固液分离,再采用体积为100mL、浓度为20g/L的硫酸铵溶液淋洗尾矿,采用EDTA滴定法测试浸出稀土离子的总质量,采用关系式(1)计算稀土矿样的品位为1.01‰;
关系式(1):
关系式(1)中:η为稀土品位,mRE为浸出的稀土离子的总质量,ms为稀土矿样的质量,MREO和MRE分别为稀土氧化物和稀土离子的摩尔质量;
第三步,建立浸矿过程的源汇项模型:
离子型稀土矿样的稀土浸出包括三个过程:(1)铵根离子通过扩散层达到稀土矿样颗粒表面,(2)铵根离子与吸附在稀土矿样颗粒表面的稀土离子发生离子交换反应,(3)解吸下来的稀土离子通过扩散层进入外部溶液;
铵根离子和稀土离子通过扩散层的过程采用Fick定律描述,可以分别通过关系式(2)和关系式(3)进行计算;
关系式(2)和关系式(3):
关系式(2)和关系式(3)中:NRE为吸附在稀土矿样颗粒表面的稀土离子的摩尔数,t为时间,d为求导符号,ms为稀土矿样的质量,和分别为外部溶液中铵根离子和稀土离子的浓度,和分别为靠近稀土矿样颗粒表面的液相铵根离子和液相稀土离子的浓度,和αRE分别为与铵根离子和稀土离子在扩散层中的扩散系数相关的参数,和MRE分别为铵根离子和稀土离子的摩尔质量;
铵根离子与吸附在稀土矿样颗粒表面的稀土离子发生离子交换反应采用关系式(4)描述;
关系式(4):
浸出过程中,以稀土矿样颗粒、扩散层和外部溶液所组成的体系为研究对象,在整个体系内,铵根离子和稀土离子应满足质量守恒,分别采用关系式(5)和关系式(6)描述;
关系式(5)和关系式(6):
关系式(5)和关系式(6)中:ms为稀土矿样的质量,VL为外部溶液体积,和分别为吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,和分别为外部溶液中铵根离子和稀土离子的浓度,和分别为浸矿前吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,和分别为浸矿前外部溶液中铵根离子和稀土离子的浓度;
以稀土矿样颗粒为研究对象,浸矿前后,稀土矿样颗粒表面的吸附位数量不发生变化,则浸矿前后吸附在稀土矿样颗粒表面的铵根离子和稀土离子应满足关系式(7);
关系式(7):
把关系式(2)、关系式(3)和关系式(5)~关系式(7)代入关系式(4)中,可以得到关系式(8);
关系式(8):
关系式(8)中:K为离子交换选择系数,和MRE分别为铵根离子和稀土离子的摩尔质量,和分别为浸矿前吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,和分别为浸矿前外部溶液中铵根离子和稀土离子的浓度,为外部溶液中稀土离子的浓度,为吸附在稀土矿样颗粒表面的稀土离子的浓度,和αRE分别为与铵根离子和稀土离子在扩散层中的扩散系数相关的参数,t为时间,d为求导符号,ms为稀土矿样的质量,VL为外部溶液体积;
关系式(8)中为吸附在稀土矿样颗粒表面的稀土离子浓度随时间的变化速率,即为浸出过程的源汇项,以为基本未知量,关系式(8)为关于的一元三次方程,采用卡丹公式解关系式(8),即可得到离子型稀土浸矿过程的源汇项,结果为关系式(9),关系式(9)中的参数可以通过关系式(10)~关系式(15)确定;
关系式(9)~关系式(15):
关系式(9)~关系式(15)中:p、q、a、b、c和e为计算过程的临时参数,为吸附在稀土矿样颗粒表面的稀土离子的浓度,和MRE分别为铵根离子和稀土离子的摩尔质量,和αRE分别为与铵根离子和稀土离子在扩散层中的扩散系数相关的参数,和分别为浸矿前外部溶液中铵根离子和稀土离子的浓度,为外部溶液中稀土离子的浓度,和分别为浸矿前吸附在稀土矿样颗粒表面的铵根离子和稀土离子的浓度,K为离子交换选择系数,t为时间,d为求导符号,ms为稀土矿样的质量,VL为外部溶液体积;
第四步,确定源汇项模型中的参数:
取20.0g聚合稀土矿样装入250mL的离心瓶中,将离心瓶放置于磁力搅拌器上,倒入体积为100mL、浓度为20.0g/L的硫酸铵溶液,同时启动磁力搅拌器,转速为1000r/min;并计时,到达预定反应时间:0s、5s、10s、15s、20s、30s、60s、90s、150s后,迅速将离心盒置于离心机中进行离心,转速设置为3000r/min,离心时间为20s,待离心结束后,采用EDTA滴定法测试不同反应时间浸出液中稀土离子的浓度,确定外部溶液中稀土离子的生成速率,通过关系式(16)得到稀土矿样颗粒表面吸附的稀土离子浓度随时间的变化速率的试验数据,
关系式(16):
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811397711.1A CN109207722B (zh) | 2018-11-22 | 2018-11-22 | 一种确定离子型稀土浸矿过程源汇项的方法 |
US16/691,632 US11873541B2 (en) | 2018-11-22 | 2019-11-22 | Method for determining source sink term in ionic type rare earth ore leaching process |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811397711.1A CN109207722B (zh) | 2018-11-22 | 2018-11-22 | 一种确定离子型稀土浸矿过程源汇项的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109207722A CN109207722A (zh) | 2019-01-15 |
CN109207722B true CN109207722B (zh) | 2020-04-17 |
Family
ID=64994086
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811397711.1A Active CN109207722B (zh) | 2018-11-22 | 2018-11-22 | 一种确定离子型稀土浸矿过程源汇项的方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US11873541B2 (zh) |
CN (1) | CN109207722B (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110195158B (zh) * | 2019-06-26 | 2021-03-05 | 江西理工大学 | 以矿石体积为依据的离子型稀土分区注液方法 |
CN111690831B (zh) * | 2020-05-29 | 2022-05-24 | 江西理工大学 | 一种离子型稀土矿的注液工艺优化方法 |
CN112329972B (zh) * | 2020-07-22 | 2023-11-10 | 矿冶科技集团有限公司 | 一种基于水环境容量的离子型稀土矿山开采预测方法及系统 |
CN112198129B (zh) * | 2020-09-15 | 2023-08-01 | 首都师范大学 | 一种不同分子量溶解性有机碳与碱土金属结合常数测定方法 |
CN112699618B (zh) * | 2020-12-18 | 2023-01-17 | 赣江新区澳博颗粒科技研究院有限公司 | 一种离子型稀土矿原地浸矿过程数值模拟方法 |
CN112964611B (zh) * | 2021-02-01 | 2024-04-05 | 攀枝花学院 | 离子型稀土堆浸渗流效率实时监测方法 |
CN113051848B (zh) * | 2021-03-25 | 2022-09-02 | 赣江新区澳博颗粒科技研究院有限公司 | 一种风化壳淋积型稀土矿溶浸过程数值模拟方法 |
CN113553787B (zh) * | 2021-07-22 | 2023-04-25 | 江西理工大学 | 一种离子型稀土矿搅拌浸出过程的数值模拟方法 |
CN113679997B (zh) * | 2021-08-23 | 2022-07-26 | 湖南省环境保护科学研究院 | 用于处理砷渣的稳定化药剂的制备方法及处理砷渣的方法 |
CN113774238B (zh) * | 2021-09-15 | 2022-10-14 | 江西离子型稀土工程技术研究有限公司 | 一种串级淋洗离子型稀土尾矿中氨氮的工艺 |
CN114297951B (zh) * | 2021-12-15 | 2024-04-05 | 中冶长天国际工程有限责任公司 | 一种采用加压注液的离子型稀土开采方法 |
CN114444348B (zh) * | 2021-12-31 | 2024-04-16 | 海油发展珠海管道工程有限公司 | 一种螺旋列板涡激振动抑制装置的动力学设计方法 |
CN114544616B (zh) * | 2022-02-09 | 2024-08-06 | 江西省地质局第七地质大队(江西省地质局稀土应用研究所) | 一种快速选定离子型稀土找矿靶区的方法 |
CN114522969B (zh) * | 2022-02-11 | 2023-09-19 | 江西理工大学 | 一种离子型稀土尾矿淋洗方法 |
CN114959317B (zh) * | 2022-06-17 | 2023-10-10 | 江西理工大学 | 一种强化离子型稀土原地浸出渗流过程的装置及方法 |
CN115216653A (zh) * | 2022-08-04 | 2022-10-21 | 中国科学院赣江创新研究院 | 一种利用电场强化风化壳淋积型稀土矿浸出的方法 |
CN115389605B (zh) * | 2022-08-23 | 2024-07-26 | 江西省地质局第七地质大队(江西省地质局稀土应用研究所) | 一种快速评价离子吸附型稀土矿山复灌区资源的方法 |
CN115992318A (zh) * | 2023-02-17 | 2023-04-21 | 中国科学院赣江创新研究院 | 一种浸出离子型稀土矿中离子相稀土的方法及其应用 |
CN117808266B (zh) * | 2024-01-02 | 2024-06-21 | 中国恩菲工程技术有限公司 | 离子型稀土智能矿山管控系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108319791B (zh) | 2018-02-09 | 2021-05-11 | 江西理工大学 | 离子型稀土原地浸矿的浸矿剂浓度计算方法 |
-
2018
- 2018-11-22 CN CN201811397711.1A patent/CN109207722B/zh active Active
-
2019
- 2019-11-22 US US16/691,632 patent/US11873541B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
US11873541B2 (en) | 2024-01-16 |
US20200199708A1 (en) | 2020-06-25 |
CN109207722A (zh) | 2019-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109207722B (zh) | 一种确定离子型稀土浸矿过程源汇项的方法 | |
Bronswijk et al. | Field‐scale solute transport in a heavy clay soil | |
Chen et al. | Pore scale study of multiphase multicomponent reactive transport during CO2 dissolution trapping | |
CN104612675B (zh) | 一种碳酸盐岩地层随钻岩性快速识别方法 | |
CN108319791A (zh) | 离子型稀土原地浸矿的浸矿剂浓度计算方法 | |
CN103510933B (zh) | 一种渗吸剂渗吸采油效果的定量评价方法 | |
US11891678B2 (en) | Method for optimizing liquid injection process of ionic rare earth ore | |
CN106522935A (zh) | 实验确定碳酸盐岩油气藏酸压裂缝导流能力分布的方法 | |
CN110410054A (zh) | 一种煤层气井压裂裂缝变时空导流能力的预测方法 | |
CN102560099B (zh) | 串联柱浸试验装置 | |
Zhao et al. | Effect of reactive surface areas associated with different particle shapes on chemical-dissolution front instability in fluid-saturated porous rocks | |
CN107130951A (zh) | 一种监测蒸汽驱井间连通状况的方法及其系统 | |
WO2022227822A1 (zh) | 一种孔隙型可渗透岩石的仿真模拟方法及系统 | |
CN112464582B (zh) | 一种基于离子型稀土孔隙优先流的原地浸矿分析方法 | |
CN102564900B (zh) | 聚合物溶液在地层不同位置处渗流过程模拟测试方法 | |
CN106959261A (zh) | 一种预测沉积物颗粒分布与配比的方法 | |
Russo et al. | Water flow and solute transport in soils: Developments and applications in memoriam Eshel Bresler (1930–1991) | |
Stollenwerk | Molybdate transport in a chemically complex aquifer: Field measurements compared with solute‐transport model predictions | |
CN110195158B (zh) | 以矿石体积为依据的离子型稀土分区注液方法 | |
CN104809275A (zh) | 碳酸盐岩微观流动模拟方法及装置 | |
CN107858537B (zh) | 离子型稀土裸脚式矿山原地浸矿孔网参数设计方法 | |
CN107798211B (zh) | 地下烃源岩层中有机质的有机酸生成量的确定方法和装置 | |
CN105938084A (zh) | 一种化学渗吸剂渗透性能评价方法 | |
CN109030775A (zh) | 一种封闭系统下水流对可溶岩溶蚀的分析设备及分析方法 | |
CN103437745A (zh) | 模拟聚合物对储层出砂及防砂效果影响的实验装置及其填砂模型 |
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 |