CN113051848B - 一种风化壳淋积型稀土矿溶浸过程数值模拟方法 - Google Patents
一种风化壳淋积型稀土矿溶浸过程数值模拟方法 Download PDFInfo
- Publication number
- CN113051848B CN113051848B CN202110319561.8A CN202110319561A CN113051848B CN 113051848 B CN113051848 B CN 113051848B CN 202110319561 A CN202110319561 A CN 202110319561A CN 113051848 B CN113051848 B CN 113051848B
- Authority
- CN
- China
- Prior art keywords
- rare earth
- earth ore
- leaching process
- weathering crust
- crust elution
- 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 137
- 238000000034 method Methods 0.000 title claims abstract description 133
- 150000002910 rare earth metals Chemical class 0.000 title claims abstract description 125
- 238000002386 leaching Methods 0.000 title claims abstract description 113
- 230000008569 process Effects 0.000 title claims abstract description 99
- 238000004088 simulation Methods 0.000 title claims abstract description 20
- 238000010828 elution Methods 0.000 claims abstract description 25
- 239000007787 solid Substances 0.000 claims abstract description 23
- 230000008878 coupling Effects 0.000 claims abstract description 21
- 238000010168 coupling process Methods 0.000 claims abstract description 21
- 238000005859 coupling reaction Methods 0.000 claims abstract description 21
- 238000012546 transfer Methods 0.000 claims abstract description 15
- 230000005514 two-phase flow Effects 0.000 claims abstract description 15
- 238000013178 mathematical model Methods 0.000 claims abstract description 10
- 239000011148 porous material Substances 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 claims description 9
- -1 rare earth ions Chemical class 0.000 claims description 9
- 239000003795 chemical substances by application Substances 0.000 claims description 7
- 238000009792 diffusion process Methods 0.000 claims description 5
- 239000007788 liquid Substances 0.000 claims description 5
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 239000000376 reactant Substances 0.000 claims description 3
- 238000005325 percolation Methods 0.000 claims description 2
- 230000007246 mechanism Effects 0.000 abstract description 8
- 238000004458 analytical method Methods 0.000 abstract description 6
- 238000011161 development Methods 0.000 description 5
- 230000001808 coupling effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004925 denaturation Methods 0.000 description 1
- 230000036425 denaturation Effects 0.000 description 1
- 238000003795 desorption Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000005342 ion exchange Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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)
- Physics & Mathematics (AREA)
- Materials Engineering (AREA)
- Mechanical Engineering (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Metallurgy (AREA)
- Manufacturing & Machinery (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Organic Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Geochemistry & Mineralogy (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Fluid Mechanics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种风化壳淋积型稀土矿溶浸过程数值模拟方法，其包括如下步骤：构建风化壳淋积型稀土矿溶浸过程的物理模型；根据风化壳淋积型稀土矿溶浸过程的物理模型构建风化壳淋积型稀土矿溶浸过程的数学模型；对稀土矿溶浸过程液固两相流多场耦合流动方程、稀土矿层骨架变形方程以及稀土矿溶浸过程传质方程进行求解以获取风化壳淋积型稀土矿溶浸过程的数据信息。本发明可以快速地描述风化壳淋积型稀土矿溶浸过程中渗流场、应力场、浓度场这三个物理场耦合作用机理，解决了传统的风化壳淋积型稀土矿溶浸过程分析方法存在时间周期长、难以描述风化壳淋积型稀土矿溶浸过程中多个物理场的耦合机理的问题。
Description
技术领域
本发明涉及计算机数值模拟技术领域，具体而言，涉及一种风化壳淋积型稀土矿溶浸过程数值模拟方法。
背景技术
风化壳淋积型稀土矿是我国宝贵的矿产资源，含有丰富的多种稀土元素。随着科技的发展，世界各国对稀土资源，特别是风化壳淋积型稀土矿中的多种稀土元素的需求量日益增长。而长期以来，由于风化壳淋积型稀土矿的粗放式开采，导致稀土资源的利用率不高，并且引发了一系列的环境问题，成为制约稀土产业可持续发展的重要因素。如何最大限度地开采出稀土矿中的稀土元素，是可持续发展的要求，也是缓解我国供需矛盾的战略要求。为从根本上解决开采稀土能力有限和安全环境隐患问题，急需开展风化壳淋积型稀土矿溶浸稀土过程的基础理论体系与关键技术研究。
风化壳淋积型稀土矿的溶浸过程是一种典型的固液两相流过程，涉及到质量扩散、异相化学反应等一系列复杂反应的过程。浸出过程中，溶液流动、液体压力、目的稀土离子的溶浸剂之间的化学反应以及目的稀土离子的吸附和解吸将导致多孔介质的弹性变性，同时也影响着溶浸过程进行的快慢程度。传统的风化壳淋积型稀土矿溶浸过程分析方法主要是基于物理实验来进行研究分析，其存在时间周期长、难以描述风化壳淋积型稀土矿溶浸过程中多个物理场的耦合机理的问题。
发明内容
基于此，为了解决传统的风化壳淋积型稀土矿溶浸过程分析方法存在时间周期长、难以描述风化壳淋积型稀土矿溶浸过程中多个物理场的耦合机理的问题，本发明提供了一种风化壳淋积型稀土矿溶浸过程数值模拟方法，其具体技术方案如下：
一种风化壳淋积型稀土矿溶浸过程数值模拟方法，其包括如下步骤：
构建风化壳淋积型稀土矿溶浸过程的物理模型；
根据所述风化壳淋积型稀土矿溶浸过程的物理模型构建风化壳淋积型稀土矿溶浸过程的数学模型；
其中，所述数学模型包括稀土矿溶浸过程液固两相流多场耦合流动方程
稀土矿层骨架变形方程
以及稀土矿溶浸过程传质方程
Xp、Xf分别为骨架孔隙变形系数和流动变形系数，p为液体压力梯度，e为水力坡度，k为渗透率，η为溶液粘度，g为重力加速度，ρf为溶液密度，εv为应力体积比，Qs为源项，θ为固体孔隙率，E为杨氏模量，S为位移变量，v为泊松比，σ为应力矩阵，εv为体积应变量，ε为各向同性弹塑性固体体积应变量，εij为柯西应变张量，εxx、εyy、εzz分别为ε沿坐标轴方向的二阶偏导数，y为轴坐标，t为时间，C1、C2分别为溶浸剂和浸出稀土离子的浓度，u为孔隙中的溶液流速，D为水动力弥散系数，b为孔隙张开度，R为阻滞系数，Cmax为溶液中稀土离子质量浓度的最大值，β为化学反应式中反应物的计量系数，x为扩散厚度，G为矿石品位；
对所述稀土矿溶浸过程液固两相流多场耦合流动方程、稀土矿层骨架变形方程以及稀土矿溶浸过程传质方程进行求解以获取风化壳淋积型稀土矿溶浸过程的数据信息。
所述风化壳淋积型稀土矿溶浸过程数值模拟方法通过根据所述风化壳淋积型稀土矿溶浸过程的物理模型构建风化壳淋积型稀土矿溶浸过程的数学模型，可以快速地描述风化壳淋积型稀土矿溶浸过程中渗流场、应力场、浓度场这三个物理场耦合作用机理，解决了传统的风化壳淋积型稀土矿溶浸过程分析方法存在时间周期长、难以描述风化壳淋积型稀土矿溶浸过程中多个物理场的耦合机理的问题，为稀土高效开发提供一种有力工具。
进一步地，通过流体数值模拟软件Fluent对所述稀土矿溶浸过程液固两相流多场耦合流动方程、稀土矿层骨架变形方程以及稀土矿溶浸过程传质方程进行求解以获取风化壳淋积型稀土矿溶浸过程的数据信息。
进一步地，所述稀土矿溶浸过程液固两相流多场耦合流动方程、稀土矿层骨架变形方程以及稀土矿溶浸过程传质方程进行求解方法为有限体积法。
进一步地，所述物理模型为二维浸出柱模型。
进一步地，所述二维浸出柱模型填充有稀土矿固体颗粒。
进一步地，所述稀土矿溶浸过程为等温渗流过程。
进一步地，所述β设定为1。
进一步地，所述有限体积法的边界条件为速度入口以及压力出口。
进一步地，所述二维浸出柱模型的长度为50CM，宽度为20CM。
相应地，本发明提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，当所述计算机程序被处理器执行时实现如上述所述的风化壳淋积型稀土矿溶浸过程数值模拟方法。
附图说明
从以下结合附图的描述可以进一步理解本发明。图中的部件不一定按比例绘制，而是将重点放在示出实施例的原理上。在不同的视图中，相同的附图标记指定对应的部分。
图1是本发明一实施例中一种风化壳淋积型稀土矿溶浸过程数值模拟方法；
图2是本发明一实施例中一种风化壳淋积型稀土矿溶浸过程数值模拟方法的物理模型的结构示意图。
具体实施方式
为了使得本发明的目的、技术方案及优点更加清楚明白，以下结合其实施例，对本发明进行进一步详细说明。应当理解的是，此处所描述的具体实施方式仅用以解释本发明，并不限定本发明的保护范围。
需要说明的是，当元件被称为“固定于”另一个元件，它可以直接在另一个元件上或者也可以存在居中的元件。当一个元件被认为是“连接”另一个元件，它可以是直接连接到另一个元件或者可能同时存在居中元件。本文所使用的术语“垂直的”、“水平的”、“左”、“右”以及类似的表述只是为了说明的目的，并不表示是唯一的实施方式。
除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施方式的目的，不是旨在于限制本发明。本文所使用的术语“及/或”包括一个或多个相关的所列项目的任意的和所有的组合。
本发明中所述“第一”、“第二”不代表具体的数量及顺序，仅仅是用于名称的区分。
如图1所示，本发明一实施例中的一种风化壳淋积型稀土矿溶浸过程数值模拟方法，其包括如下步骤：
构建风化壳淋积型稀土矿溶浸过程的物理模型；
根据所述风化壳淋积型稀土矿溶浸过程的物理模型构建风化壳淋积型稀土矿溶浸过程的数学模型；
其中，所述数学模型包括稀土矿溶浸过程液固两相流多场耦合流动方程
稀土矿层骨架变形方程
以及稀土矿溶浸过程传质方程
Xp、Xf分别为骨架孔隙变形系数和流动变形系数，p为液体压力梯度，e为水力坡度，k为渗透率，η为溶液粘度，g为重力加速度，ρf为溶液密度，εv为应力体积比，Qs为源项，θ为固体孔隙率，E为杨氏模量，S为位移变量，v为泊松比，σ为应力矩阵，εv为体积应变量，ε为各向同性弹塑性固体体积应变量，εij为柯西应变张量，εxx、εyy、εzz分别为ε沿坐标轴方向的二阶偏导数，y为沿浸出柱的轴坐标，t为时间，C1、C2分别为溶浸剂和浸出稀土离子的浓度，u为孔隙中的溶液流速，D为水动力弥散系数，b为孔隙张开度，R为阻滞系数，Cmax为溶液中稀土离子质量浓度的最大值，β为化学反应式中反应物的计量系数，x为扩散厚度，G为矿石品位，Si、Sj分别表示i轴方向以及j轴方向上的位移变量；
对所述稀土矿溶浸过程液固两相流多场耦合流动方程、稀土矿层骨架变形方程以及稀土矿溶浸过程传质方程进行求解以获取风化壳淋积型稀土矿溶浸过程的数据信息。
所述稀土矿溶浸过程液固两相流多场耦合流动方程、稀土矿层骨架变形方程以及稀土矿溶浸过程传质方程通过稀土矿多孔介质孔隙率耦合，组成了溶浸过程中渗流场、应力场以及浓度场的耦合模型。
所述风化壳淋积型稀土矿溶浸过程数值模拟方法通过根据所述风化壳淋积型稀土矿溶浸过程的物理模型构建风化壳淋积型稀土矿溶浸过程的数学模型，可以快速地描述风化壳淋积型稀土矿溶浸过程中渗流场、应力场、浓度场这三个物理场耦合作用机理，解决了传统的风化壳淋积型稀土矿溶浸过程分析方法存在时间周期长、难以描述风化壳淋积型稀土矿溶浸过程中多个物理场的耦合机理的问题，为稀土高效开发提供一种有力工具。
在其中一个实施例中，在稀土矿溶浸过程中，稀土矿可视为多孔介质，溶浸剂在渗流过程与稀土矿之间发生离子交换传递过程，稀土矿层的骨架结构会发生相应变化，溶浸剂和稀土离子利用化学反应速率的连续性耦合。
孔隙中的溶质迁移方程为
以及
其中，s为单位长度孔隙表面积吸附的溶质质量，Ri为化学反应速率。
所述化学反应速率Ri表示为
假设溶浸剂在孔隙表面为线性平衡等温吸附，溶解相与吸附相之间的关系可以表示为
故
其中，kf为分布系数。
考虑扩散质量通量(单位时间内单位面积上扩散的物质量)Jd，由Fick第一定律可得
其中，
为浓度梯度。
所述阻滞系数
在其中一个实施例中，通过流体数值模拟软件Fluent对所述稀土矿溶浸过程液固两相流多场耦合流动方程、稀土矿层骨架变形方程以及稀土矿溶浸过程传质方程进行求解以获取风化壳淋积型稀土矿溶浸过程的数据信息。
在其中一个实施例中，所述稀土矿溶浸过程液固两相流多场耦合流动方程、稀土矿层骨架变形方程以及稀土矿溶浸过程传质方程进行求解方法为有限体积法。
在其中一个实施例中，如图2所示，所述物理模型为二维浸出柱模型。
在其中一个实施例中，所述二维浸出柱模型填充有稀土矿固体颗粒。
在其中一个实施例中，所述稀土矿溶浸过程为等温渗流过程。
在其中一个实施例中，所述β设定为1。
在其中一个实施例中，所述有限体积法的边界条件为速度入口以及压力出口。通过采用速度入口以及压力出口的边界调节对所述数学模型进行封闭求解，即可得出整个稀土矿溶浸过程中的渗流场、应力场和浓度场的所有信息。
在其中一个实施例中，所述二维浸出柱模型的长度为50cm，宽度为20cm。
在其中一个实施例中，本发明提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，当所述计算机程序被处理器执行时实现如上述所述的风化壳淋积型稀土矿溶浸过程数值模拟方法。
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种风化壳淋积型稀土矿溶浸过程数值模拟方法，其特征在于，包括如下步骤：
构建风化壳淋积型稀土矿溶浸过程的物理模型；
根据所述风化壳淋积型稀土矿溶浸过程的物理模型构建风化壳淋积型稀土矿溶浸过程的数学模型；
其中，所述数学模型包括稀土矿溶浸过程液固两相流多场耦合流动方程
稀土矿层骨架变形方程
以及稀土矿溶浸过程传质方程
与
Xp、Xf分别为骨架孔隙变形系数和流动变形系数，p为液体压力梯度，e为水力坡度，k为渗透率，η为溶液粘度，g为重力加速度，ρf为溶液密度，εv为应力体积比，Qs为源项，θ为固体孔隙率，E为杨氏模量，S为位移变量，v为泊松比，σ为应力矩阵，εv为体积应变量，ε为各向同性弹塑性固体体积应变量，εij为柯西应变张量，εxx、εyy、εzz分别为ε沿坐标轴方向的二阶偏导数，y为轴坐标，t为时间，C1、C2分别为溶浸剂和浸出稀土离子的浓度，u为孔隙中的溶液流速，D为水动力弥散系数，b为孔隙张开度，R为阻滞系数，Cmax为溶液中稀土离子质量浓度的最大值，β为化学反应式中反应物的计量系数，x为扩散厚度，G为矿石品位；
对所述稀土矿溶浸过程液固两相流多场耦合流动方程、稀土矿层骨架变形方程以及稀土矿溶浸过程传质方程进行求解以获取风化壳淋积型稀土矿溶浸过程的数据信息。
2.如权利要求1所述的一种风化壳淋积型稀土矿溶浸过程数值模拟方法，其特征在于，通过流体数值模拟软件Fluent对所述稀土矿溶浸过程液固两相流多场耦合流动方程、稀土矿层骨架变形方程以及稀土矿溶浸过程传质方程进行求解以获取风化壳淋积型稀土矿溶浸过程的数据信息。
3.如权利要求2所述的一种风化壳淋积型稀土矿溶浸过程数值模拟方法，其特征在于，所述稀土矿溶浸过程液固两相流多场耦合流动方程、稀土矿层骨架变形方程以及稀土矿溶浸过程传质方程进行求解方法为有限体积法。
4.如权利要求1所述的一种风化壳淋积型稀土矿溶浸过程数值模拟方法，其特征在于，所述物理模型为二维浸出柱模型。
5.如权利要求4所述的一种风化壳淋积型稀土矿溶浸过程数值模拟方法，其特征在于，所述二维浸出柱模型填充有稀土矿固体颗粒。
6.如权利要求1所述的一种风化壳淋积型稀土矿溶浸过程数值模拟方法，其特征在于，所述稀土矿溶浸过程为等温渗流过程。
7.如权利要求1所述的一种风化壳淋积型稀土矿溶浸过程数值模拟方法，其特征在于，所述β设定为1。
8.如权利要求3所述的一种风化壳淋积型稀土矿溶浸过程数值模拟方法，其特征在于，所述有限体积法的边界条件为速度入口以及压力出口。
9.如权利要求5所述的一种风化壳淋积型稀土矿溶浸过程数值模拟方法，其特征在于，所述二维浸出柱模型的长度为50CM，宽度为20CM。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机程序，当所述计算机程序被处理器执行时实现如权利要求1-9中任意一项所述的风化壳淋积型稀土矿溶浸过程数值模拟方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110319561.8A CN113051848B (zh) | 2021-03-25 | 2021-03-25 | 一种风化壳淋积型稀土矿溶浸过程数值模拟方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110319561.8A CN113051848B (zh) | 2021-03-25 | 2021-03-25 | 一种风化壳淋积型稀土矿溶浸过程数值模拟方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113051848A CN113051848A (zh) | 2021-06-29 |
| CN113051848B true CN113051848B (zh) | 2022-09-02 |
Family
ID=76515726
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110319561.8A Active CN113051848B (zh) | 2021-03-25 | 2021-03-25 | 一种风化壳淋积型稀土矿溶浸过程数值模拟方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113051848B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113553787B (zh) * | 2021-07-22 | 2023-04-25 | 江西理工大学 | 一种离子型稀土矿搅拌浸出过程的数值模拟方法 |
| CN118112195A (zh) * | 2024-01-02 | 2024-05-31 | 中国恩菲工程技术有限公司 | 基于数字孪生的特征污染物监控方法及装置、介质、终端 |
| CN117807681B (zh) * | 2024-01-02 | 2024-05-17 | 中国恩菲工程技术有限公司 | 数字孪生矿山建模方法及装置、介质、设备 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104867294A (zh) * | 2015-05-27 | 2015-08-26 | 江西理工大学 | 一种离子型稀土原地浸矿滑坡在线监测预警方法 |
| CN106204312A (zh) * | 2016-07-21 | 2016-12-07 | 中国地质科学院矿产资源研究所 | 离子吸附型稀土矿快速经济评价方法 |
| CN108427836A (zh) * | 2018-03-02 | 2018-08-21 | 江西理工大学 | 一种定水头条件下离子型稀土单井注液影响范围的测定方法 |
| CN109207722A (zh) * | 2018-11-22 | 2019-01-15 | 江西理工大学 | 一种确定离子型稀土浸矿过程源汇项的方法 |
-
2021
- 2021-03-25 CN CN202110319561.8A patent/CN113051848B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104867294A (zh) * | 2015-05-27 | 2015-08-26 | 江西理工大学 | 一种离子型稀土原地浸矿滑坡在线监测预警方法 |
| CN106204312A (zh) * | 2016-07-21 | 2016-12-07 | 中国地质科学院矿产资源研究所 | 离子吸附型稀土矿快速经济评价方法 |
| CN108427836A (zh) * | 2018-03-02 | 2018-08-21 | 江西理工大学 | 一种定水头条件下离子型稀土单井注液影响范围的测定方法 |
| CN109207722A (zh) * | 2018-11-22 | 2019-01-15 | 江西理工大学 | 一种确定离子型稀土浸矿过程源汇项的方法 |
Non-Patent Citations (2)
| Title |
|---|
| Research progress on leaching technology and theory of weathered crust elution-deposited rare earth ore;Wenrui Nie 等;《Hydrometallurgy》;20200531;第193卷;1-11 * |
| 风化壳淋积型稀土矿浸出过程中离子交换模型与溶质运移方程;曹小晶;《中国优秀博硕士学位论文全文数据库(硕士) 基础科学辑》;20191215(第12期);A002-130 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113051848A (zh) | 2021-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113051848B (zh) | 一种风化壳淋积型稀土矿溶浸过程数值模拟方法 | |
| Luo et al. | Review on the development and utilization of ionic rare earth ore | |
| Van Orman et al. | Diffusive fractionation of trace elements during production and transport of melt in Earth’s upper mantle | |
| Huber et al. | A physical model for three‐phase compaction in silicic magma reservoirs | |
| Cugliandolo | Topics in coarsening phenomena | |
| Lin et al. | Modeling of 3D rock porous media by combining X-ray CT and Markov chain Monte Carlo | |
| Natarajan et al. | Evolution of fracture permeability due to co-colloidal bacterial transport in a coupled fracture-skin-matrix system | |
| CN117807681B (zh) | 数字孪生矿山建模方法及装置、介质、设备 | |
| CN117828732B (zh) | 基于数字孪生的边坡稳定性确定方法及系统、介质、终端 | |
| Ilankoon | Hydrodynamics of unsaturated particle beds pertaining to heap leaching | |
| Hilpert et al. | A binomial modeling approach for upscaling colloid transport under unfavorable conditions: Emergent prediction of extended tailing | |
| Altree-Williams et al. | Exact solution for coupled reactive flow and dissolution with porosity changes | |
| Saeedmonir et al. | A multiscale approach in modeling of chemically reactive porous media | |
| Mellado et al. | An analytical–numerical method for solving a heap leaching problem of one or more solid reactants from porous pellets | |
| Ambekar et al. | Particle-resolved simulations of local liquid spreading in packed beds: Effect of wettability at varying particle size | |
| Meirmanov et al. | Mesoscopic dynamics of solid-liquid interfaces. A general mathematical model | |
| Wu et al. | Simulation of coupled flowing-reaction-deformation with mass transfer in heap leaching processes | |
| Yang et al. | Simulation of pore scale fluid flow of granular ore media in heap leaching based on realistic model | |
| LU et al. | Kinetics of weathered-crust elution-deposited rare-earth ore in a leaching process | |
| Han et al. | An Experimental Study on the Kinetics of Leaching Ion-Adsorbed REE Deposits with Different Concentrations of Magnesium Sulfate | |
| Maaløe | Melt dynamics of a partially molten mantle with randomly oriented veins | |
| Ding et al. | A kinetic model for heap leaching of uranium ore considering variation of model parameters with depth of heap | |
| Warner et al. | Chemical reaction between recharge water and aquifer water | |
| Chen et al. | Pertechnetate removal by macroporous polymer impregnated with 2-nitrophenyl octyl ether (NPOE) | |
| Marion | The geochemistry of natural waters: Surface and groundwater environments |
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 |