DOI:10.19476/j.ysxb.1004.0609.2001.06.037
串级萃取最低成本优化萃取比方程
钟盛华
江西省稀土研究所 南昌330013
摘 要:
以寻求串级萃取最低生产成本的工艺条件为最优化目标 , 通过将萃取工艺参数与随之会发生变化的生产成本关联 , 建立了分馏萃取最低成本优化萃取比方程。为简化计算 , 推导出萃取段、洗涤段最低成本萃取比方程。如果不考虑酸碱消耗的成本变化 , 推导出的简略方程其优化目标与现有最优萃取比方程相一致 , 计算值相近。建立的最低成本优化萃取比方程及其优化目标发展了现有最优萃取比理论。
关键词:
串级萃取 ;最优化理论 ;成本 ;萃取比方程 ;稀土 ;
中图分类号: TF804.2
收稿日期: 2000-10-26
Equations of optimized extraction factor for minimum production cost in counter current extraction
Abstract:
Based on the optimization goal searching for extraction parameters of minimun production cost, a novel equation for calculating optimized extraction factor of minimum production cost in fractional extraction was proposed. For short cut calculation, new equations for determining the optimized extraction factors of minimum production cost in both the extraction and scrubbing sections were derived from above novel equation. These equations and theory have further developed the theory on equations of optimized extraction factors.
Keyword:
counter current extraction; theory of optimization; production cost; extraction factor; rare earth;
Received: 2000-10-26
1959年, Alders提出了分馏萃取理论
[1 ]
。 20世纪70年代国内开始有串级萃取理论的文献发表
[2 ,3 ,4 ]
, 80年代是我国串级萃取理论发展的旺盛期
[5 ,6 ,7 ,8 ]
。 北京大学徐光宪院士等创立的串级萃取理论在有色金属领域, 特别是在稀土分离工业中得到了广泛应用, 对稀土分离起到了重要的指导作用。 1978年徐光宪发表了串级萃取理论Ⅰ, 最优萃取比方程。 这个方程是以萃取段或洗涤段的萃取槽总体积不变, 分离出产品最多为优化目标推导出来的。 这个方程的建立, 为萃取分离工业选择好的工艺条件, 提高经济效益起到了很好作用。 最佳的萃取工艺应该是产生经济效益最好的工艺, 这是工业生产最需要的工艺。 然而, 它是一个十分复杂的问题。 从工艺条件角度讲, 生产合格产品的成本最低工艺就是经济效益最好的工艺。 作者从事过工艺、 设备的研究和设计
[9 ,10 ,11 ,12 ,13 ]
。 本文在串级萃取理论基础上, 将生产成本与萃取工艺参数关联, 建立了最低成本优化萃取比方程, 简化推导出萃取段和洗涤段最低成本萃取比方程。 它们在分别不考虑皂化碱消耗和洗涤酸消耗成本时, 其优化目标与徐光宪的最优萃取比方程相一致, 计算值相近。
1 分馏萃取最低成本优化萃取比方程
对于分馏萃取, 首先优选确定萃取体系、 有机相组成、 皂化率、 料液浓度和洗酸酸度等。 在确定原料组成和分离指标后, 即在两头都出合格产品情况下, 寻求该分馏萃取最低成本的萃取工艺参数, 也就是建立最低成本优化萃取比方程。
令Q 为这个分馏萃取工程的处理原料能力, 单位为t/月。
V
?
S , VF , VW 为每分钟进1 mol 原料时有机相、 料液、 洗酸的流量, 单位为mL /min 。 本文符号除文中标明外, 均与北京大学串级萃取理论中含意相同, 可参见文献
[
14 ]
。
进萃取槽有机相流量为
Q
?
S
/
(
L
?
min
-
1
)
=
1
0
6
?
Q
3
0
×
2
4
×
6
0
Μ
?
×
V
?
S
/
1
0
0
0
=
2
.
3
1
5
×
1
0
-
2
Q
Μ
?
×
V
?
S
?
?
?
(
1
)
进萃取槽洗酸流量为
Q
W
/
(
L
?
min
-
1
)
=
2
.
3
1
5
×
1
0
-
2
Q
Μ
?
×
V
W
?
?
?
(
2
)
进萃取槽料液流量为
Q
F
/
(
L
?
min
-
1
)
=
2
.
3
1
5
×
1
0
-
2
Q
Μ
?
×
V
F
?
?
?
(
3
)
单级萃取槽体积为
V单
=
2
.
3
1
5
×
1
0
-
2
Q
(
1
+
r
)
Μ
?
?
t
(
V
?
S
+
V
F
+
V
W
)
?
?
?
(
4
)
式中 t为混合室中混合时间, min ; r为澄清室与混合室的边长比。
由于原料组成、 分离指标确定, 有机相出口分数f′A 也就确定。 分离规模和萃取槽确定, 反萃段也就确定。 它不会随工艺参数的变化而改变, 推导中忽略反萃段, 只考虑萃取段、 洗涤段。 本文中只讨论水相进料。
随萃取工艺参数 (n, m,
V
?
S ∶VF ∶VW ) 变化, 会发生变动的生产成本, 以及它们之间的相互关系如下。
1) 萃取设备投资。 它与 (n+m) 级数和每级萃取槽体积乘积成正比。 其投资的折旧成本为
c1 =k1 (n+m) ·V单 (5)
式中 c1 表示生产成本, 元/月; k1 为系数, 以下类同。
2) 萃取设备所占厂房投资。 当萃取设备选用混合澄清萃取器, 并确定其边高比后, 其投资就与 (n+m) 级数和萃取槽体积的乘积成正比。 投资的折旧成本为
c2 =k2 (n+m) ·V单 (6)
3) 萃取槽中储存的萃取剂有机相及金属量所占用资金与萃取槽总体积成正比, 其利息为
c3 =k3 (n+m) ·V单 (7)
4) 萃取槽的动力消耗, 也即搅拌的电耗。 它与萃取槽混合室容积大小成正比, 与 (n+m) 也成正比。 当混合室与澄清室容积比确定, 其电耗成本为
c4 =k4 (n+m) ·V单 (8)
5) 皂化萃取剂消耗碱与有机相流量成正比。
c
5
=
k
5
?
Q
?
S
?
?
?
(
9
)
6) 洗涤消耗的酸与洗酸流量成正比。
c6 =k6 ·QW (10)
萃取的高低位储桶、 给料装置及其所占厂房等投资与萃取工艺参数的变动无关或关系不大。 有机相损耗、 维修费、 化验费、 工资等与萃取工艺参数变动关系也不大, 所以在推导方程中不加考虑。
将式 (5) ~ (10) 相加, 式 (1) ~ (3) 代入, 得到随萃取工艺参数变化发生变动的成本总和C:
C
=
2
.
3
1
5
×
1
0
-
2
Q
Μ
?
[
1
+
r
t
?
(
k
1
+
k
2
+
k
3
+
k
4
)
×
(
V
?
S
+
V
W
+
V
F
)
?
(
n
+
m
)
+
k
5
V
?
S
+
k
6
V
W
]
?
?
?
(
1
1
)
因
V
?
S
=
S
?
/
c
—
S
V
F
=
Μ
F
/
c
F
=
1
/
c
F
?
?
?
(
1
2
)
对于稀土萃取分离, 因稀土离子为+3价
VW =W·3/cW (13)
W
=
S
?
-
f
′
A
?
?
?
(
1
4
)
将 (12) , (13) , (14) 式代入 (11) 式, 得
C
=
2
.
3
1
5
×
1
0
-
2
Q
Μ
?
[
Κ
1
(
n
+
m
)
S
?
+
Κ
2
(
n
+
m
)
+
Κ
3
S
?
+
Κ
4
W
]
?
?
?
(
1
5
)
式中
Κ
1
=
(
k
1
+
k
2
+
k
3
+
k
4
)
?
1
+
r
t
(
1
/
c
—
S
+
3
/
c
W
)
Κ
2
=
(
k
1
+
k
2
+
k
3
+
k
4
)
?
1
+
r
t
(
1
/
c
—
F
-
3
f
′
A
/
c
W
)
Κ
3
=
k
5
/
c
—
S
Κ
4
=
3
k
6
/
c
W
对 (15) 式求导, 并令C′=0, 得
C
′
=
Κ
1
(
n
′
+
m
′
)
S
?
+
Κ
1
(
n
+
m
)
S
?
′
+
Κ
2
(
n
′
+
m
′
)
+
Κ
3
S
?
′
+
Κ
4
W
′
=
0
?
?
?
(
1
6
)
式 (16) 是最低成本工艺参数优化方程。
由北京大学串级萃取理论,
S
?
=
E
Μ
f
′
B
1
-
E
Μ
?
?
?
(
1
7
)
对式 (14) 求导, 式 (17) 对EM 求导, 得
W
′
=
S
?
′
=
f
′
B
(
1
-
E
Μ
)
2
?
?
?
(
1
8
)
将式 (17) 和 (18) 代入式 (16) , 得
C′=K1 (n′+m′) (EM -E
Μ
2
) +
K2 (n′+m′) (1-EM ) 2 /f′B ?+
K1 (n+m) +K3 +K4 =0 (19)
由北京大学串级萃取理论, 洗涤段萃取比E′M 与萃取段萃取比EM 关系为
E
′
Μ
=
E
Μ
f
′
B
E
Μ
-
f
′
A
?
?
?
(
2
0
)
n
=
lg
?
b
lg
?
β
E
Μ
?
?
?
(
2
1
)
m
=
lg
?
a
lg
(
β
′
/
E
′
Μ
?
)
-
1
=
lg
?
a
lg
β
′
(
E
Μ
-
f
′
A
?
)
E
Μ
f
′
B
-
1
?
?
?
(
2
2
)
式 (21) , 式 (22) 对EM 求导, 得
n
′
=
-
0
.
4
3
4
?
lg
?
b
(
lg
?
β
E
Μ
)
2
E
Μ
?
?
?
(
2
3
)
m
′
=
-
0
.
4
3
4
f
′
A
?
lg
?
a
(
E
Μ
-
f
′
A
?
)
E
Μ
?
[
lg
β
′
(
E
Μ
-
f
′
A
?
)
E
Μ
f
′
B
]
-
2
?
?
?
(
2
4
)
式 (19) 就是本文建立的最低成本优化萃取比方程。 将式 (21) ~ (24) 代入式 (19) , 得到只含未知数EM 的最低成本优化萃取比方程。 但是这个方程展开很复杂, 且是EM 的高次方程。 可采用数学计算方法来求解这个方程, 如可采用Muller 迭代法来求解。 若用计算机来求解这方程会方便很多。 限于篇幅, 本文不进行讨论。 求出的EM 就是最低成本优化萃取比。 再由这萃取比EM , 应用北京大学串级理论可计算出最低成本的萃取工艺参数 (n, m ,
V
?
S ∶VF ∶VW ) 。
2 分段计算的最低成本萃取比方程
将式 (15) 变为
C
=
2
.
3
1
5
×
1
0
-
2
Q
Μ
?
[
(
Κ
1
n
S
?
+
Κ
2
n
+
Κ
3
S
?
)
+
(
Κ
1
m
S
?
+
Κ
2
m
+
Κ
4
W
)
]
?
?
?
(
2
5
)
令
C
=
C
E
+
C
S
C
E
=
2
.
3
1
5
×
1
0
-
2
Q
Μ
?
(
Κ
1
n
S
?
+
Κ
2
n
+
Κ
3
S
?
)
?
?
?
(
2
6
)
C
S
=
2
.
3
1
5
×
1
0
-
2
Q
Μ
?
(
Κ
1
m
S
?
+
Κ
2
m
+
Κ
4
W
)
?
?
?
(
2
7
)
可以看出, CE 就是萃取段随萃取工艺参数变化发生变动的成本总和, CS 就是洗涤段随萃取工艺参数变化发生变动的成本总和。
2.1 萃取段最低成本萃取比方程
式 (26) 是萃取段萃取工艺参数与变动总成本关联式, 对其求导, 并令C′E =0, 得
Κ
1
n
′
S
?
+
Κ
1
n
S
?
′
+
Κ
2
n
′
+
Κ
3
S
?
′
=
0
?
?
?
(
2
8
)
由式 (21)
n
=
lg
?
b
lg
?
β
E
Μ
=
ln
?
b
ln
?
β
E
Μ
?
?
?
(
2
9
)
已知
ln
Ν
+
1
Ν
=
2
(
1
2
Ν
+
1
+
1
3
?
1
(
2
Ν
+
1
)
3
+
1
5
?
1
(
2
Ν
+
1
)
5
+
?
?
)
?
?
?
(
3
0
)
式中
Ν
+
1
Ν
>
0
,
式 (30) 的右边只取第一项
ln
Ν
+
1
Ν
≈
2
2
Ν
+
1
?
?
?
(
3
1
)
令
β
E
Μ
=
Ν
+
1
Ν
, 则
Ν
=
1
β
E
Μ
-
1
, 代入 (31) 式, 得
ln
β
E
Μ
≈
2
(
β
E
Μ
-
1
)
β
E
Μ
+
1
?
?
?
(
3
2
)
将式 (32) 代入式 (29) , 得
n
=
(
β
E
Μ
+
1
)
ln
?
b
2
(
β
E
Μ
-
1
)
?
?
?
(
3
3
)
将式 (33) 对E M 求导, 得
n
′
=
-
β
ln
?
b
(
β
E
Μ
-
1
)
2
?
?
?
(
3
4
)
将式 (17) , (18) , (33) 和 (34) 代入式 (28) 中, 得
pE
Μ
2
+qE M +r =0 (35)
式中
{
p
=
2
Κ
1
β
?
ln
?
b
+
Κ
1
β
2
ln
?
b
-
2
β
(
Κ
2
/
f
′
B
?
)
ln
?
b
+
?
?
2
Κ
3
β
2
q
=
4
β
(
Κ
2
/
f
′
B
?
)
ln
?
b
-
2
Κ
1
β
?
ln
?
b
-
4
Κ
3
β
r
=
2
Κ
3
-
Κ
1
?
ln
?
b
-
2
β
(
Κ
2
/
f
′
B
?
)
ln
?
b
?
?
?
(
3
6
)
E
Μ
=
-
q
±
q
2
-
4
p
r
2
p
?
?
?
(
3
7
)
式 (35) 就是萃取段最低成本萃取比方程。 由式 (37) 算出的E M 就是萃取段最低成本萃取比。
2.2 洗涤段最低成本萃取比方程
式 (27) 是洗涤段萃取工艺参数与变动总成本的关联式, 对其求导, 并令C ′S ?=0, 由式 (18) , 得
Κ
1
m
S
?
′
+
Κ
1
m
′
S
?
+
Κ
2
m
′
+
Κ
4
S
?
′
=
0
?
?
?
(
3
8
)
由式 (31) , 令β ′/E ′M ?
=
Ν
+
1
Ν
, 则
Ν
=
E
′
Μ
?
β
′
-
E
′
Μ
?
代入式 (31) , 得
ln (β ′/E ′M ?)
=
2
β
′
-
2
E
′
Μ
?
β
′
+
E
′
Μ
?
?
?
?
(
3
9
)
将 (39) 式代入 (22) 式中, 得
m
=
(
β
′
+
E
′
Μ
?
)
ln
?
a
2
β
′
-
2
E
′
Μ
?
-
1
?
?
?
(
4
0
)
将式 (40) 对E ′M ?求导, 得
m
′
=
β
′
ln
?
a
(
β
′
-
E
′
Μ
?
)
2
?
?
?
(
4
1
)
由串级理论,
S
?
=
E
′
Μ
?
f
′
A
?
E
′
Μ
?
-
1
?
?
?
(
4
2
)
将式 (42) 对E ′M ?求导, 得
S
?
′
=
-
f
′
A
?
(
E
′
Μ
?
-
1
)
2
?
?
?
(
4
3
)
将式 (40) ~ (43) 代入式 (38) , 得
p (E ′M ?) 2 +qE ′M ?+r =0 (44)
式中
{
p
=
2
Κ
1
β
′
ln
?
a
+
Κ
1
ln
?
a
+
2
(
Κ
2
/
f
′
A
?
)
β
′
ln
?
a
+
?
2
(
Κ
1
-
Κ
4
)
q
=
-
[
2
Κ
1
β
′
ln
?
a
+
4
(
Κ
2
/
f
′
A
?
)
β
′
ln
?
a
+
?
4
(
Κ
1
-
Κ
4
)
β
′
]
r
=
2
(
Κ
2
/
f
′
A
?
)
β
′
ln
?
a
-
Κ
1
(
β
′
)
2
ln
?
a
+
?
2
(
Κ
1
-
Κ
4
)
(
β
′
)
2
E ′M ?
=
-
q
±
q
2
-
4
p
r
2
p
?
?
?
(
4
5
)
式 (44) 就是洗涤段最低成本萃取比方程。
3 简略方程
由式 (12) , (13) , (14) , 得
V
?
S
+
V
W
+
V
F
=
(
V
?
S
+
S
?
?
3
/
c
W
)
+
(
V
F
-
3
f
′
A
?
/
c
W
)
?
?
?
(
4
6
)
有些情况下, 料液流量V F 比有机相流量
V
?
S 小得多, 比洗酸流量VW 也小许多。 在式 (46) 中 (VF -3f′A ?/cW ) 比
(
V
?
S
+
S
?
?
3
/
c
W
)
小很多, 可以将 (VF -3f′A ?/cW ) 略去不计。 从式 (15) 知, 此时K2 =0。 这样式 (26) 和式 (27) 简化为
C
E
=
2
.
3
1
5
×
1
0
-
2
Q
Μ
?
(
Κ
1
n
S
?
+
Κ
3
S
?
)
?
?
?
(
4
7
)
C
S
=
2
.
3
1
5
×
1
0
-
2
Q
Μ
?
(
Κ
1
m
S
?
+
Κ
4
W
)
?
?
?
(
4
8
)
3.1 萃取段简略方程
如果我们不考虑有机相皂化消耗碱的成本, 式 (47) 中K3 =0, (47) 式就为
C
E
=
2
.
3
1
5
×
1
0
-
2
Q
Μ
?
Κ
1
n
S
?
?
?
?
(
4
9
)
对 (49) 式求导, 并令C′E ?=0, 将式 (33) 和 (18) 代入, 得
β (β+2) E
Μ
2
-2βEM -1=0 (50)
因为EM 不为负值, 所以
E
Μ
=
β
+
2
β
2
+
2
β
β
(
β
+
2
)
?
?
?
(
5
1
)
从 (36) 式也可以得出 (50) 式, 因K2 =0, K3 =0, 代入 (36) 式, 得 (50) 式。
式 (51) 是不考虑有机相皂化消耗碱的成本变化, 且 (VF -3f′A ?/cW ) 忽略不计的情况下萃取段最低成本萃取比简略方程。
从式 (11) 和式 (15) 可以看出式 (49) 中CE 是与萃取槽总体积变化相关的成本近似值。 对式 (49) 求导, 并令C′E ?=0, 得到方程 (50) 。 其优化目标变为生产合格产品相同产量所用萃取段萃取槽总体积最小, 这与徐光宪的最优萃取比方程优化目标——萃取段萃取槽总体积不变产出产品B1 最多, 应该是一致的, 只是求解过程形式不一样。 表1列出了萃取段不同分离系数β用本文式 (51) 和采用徐光宪最优萃取比方程计算出的EM 值。 从表中看出, 计算的EM 值相近。
3.2 洗涤段简略方程
如果我们不考虑洗酸消耗的成本变化, 也即
表1 用本文 (51) 方程和最优萃取比方程计算萃取段不同β的EM值
Table 1 Calculated values of optimized extraction factors E M at various β in extraction section by different equations
Separation factor β
E
Μ
=
1
/
β
[
2
]
E
Μ
=
2
+
0
.
1
2
5
β
1
+
1
.
1
2
5
β
[
2
]
Eq. (51) in this paper
1.2
0.913
0.915
0.911
1.5
0.816
0.814
0.807
1.8
0.745
0, 736
0.727
2.0
0.707
0.692
0.683
2.3
0.659
0.638
0.627
2.6
0.620
0.592
0.579
3.0
0.577
0.543
0.527
K 4 =0, 同理可得
(2β ′+1) (E ′M ?) 2 -2β ′E ′M ?- (β ′) 2 =0 (52)
因为E ′M ?不为负值, 所以
E ′M ?
=
β
′
+
β
′
2
(
1
+
β
′
)
2
β
′
+
1
?
?
?
(
5
3
)
式 (53) 是在不考虑洗酸消耗成本变化, 且 (V F -3f ′A ?/c W ) 忽略不计的情况下洗涤段最低成本萃取比简略方程。
同理, 式 (52) 也与徐光宪最优萃取比方程的优化目标相一致, 只是推导过程及表达形式不一样。 表2列出了洗涤段不同分离系数β ′用本文式 (53) 和采用徐光宪最优萃取比方程计算出的E ′M ?值。 从表中看出, 计算结果相近。
表2 用式 (53) 和最优萃取比方程计算洗涤段不同β′的优化E′M?值
Table 2 Calculated values of optimized extraction factors E ′M ?? at various β ′ in scrubbing section by different equations
Separation factor β
E ′M ?
≈
β
′
[
2
]
E ′M ?
=
1
+
0
.
1
2
5
β
′
2
+
0
.
1
2
5
β
′
[
2
]
Eq. (53) in this paper
1.2
1.095
1.093
1.095
1.5
1.22
1.23
1.21
1.8
1.34
1.36
1.32
2.0
1.41
1.44
1.38
2.3
1.52
1.57
1.47
2.6
1.61
1.68
1.54
3.0
1.73
1.84
1.64
4 结论
1) 从串级萃取最佳工艺应该是生产成本最低的工艺出发, 将萃取工艺参数与随之发生变动的生产成本关联, 建立了最低成本优化萃取比方程。
2) 为了简化计算, 推导出萃取段和洗涤段最低成本萃取比方程。 由于忽略了萃取段与洗涤段的相互联系, 求得的萃取比不如最低成本优化萃取比方程的优化。
3) 推导了萃取段最低成本萃取比简略方程和洗涤段最低成本萃取比简略方程。
所建立的最低成本优化萃取比方程及优化目标发展了徐光宪的最优萃取比串级理论。
参考文献
[1] ] AldersL .LiquidLiquidExtraction[M ].Amsterdam:ElsevierPublishingCompany, 1959.
[2] XUGuang xian (徐光宪) .串级萃取理论 (Ⅰ) , (Ⅱ) [J].ActaScientiarumNoturaliumUniversitatisPekinen sis (北京大学学报 (自然科学版) ) , 1978 (1) :51-75.
[3] LUOZhen hua (罗振华) .溶剂萃取分离基本计算[J].CricularofChemistry (化学通报) , 1977 (2) :33-37.
[4] JINTian zhu (金天柱) , XUGuang xian (徐光宪) .稀土萃取串级理论—多组份串级萃取体系质量分布的计算方法[J].ChinJRareMetals (稀有金属) , 1979 (3) :1-14.
[5] XUGuang xian (徐光宪) , LIBiao guo (李标国) , NIYa ming (倪亚明) .串级萃取理论 (Ⅳ) —串级萃取参数的极值公式和工艺设计[J].ChinJRareMetals (稀有金属) , 1982 (1) :10-19.
[6] XUGuang xian (徐光宪) , LIBiao guo (李标国) , YANChun hua (严纯华) .串级萃取理论的进展及其在稀土工艺中的应用[J].ChinRareEarths (稀土) , 1985 (1) :56-67.
[7] LiBiao guo (李标国) , XUGuang xian (徐光宪) .串级萃取理论—动态平衡的计算和程序设计[J].ChinJRareMetals (稀有金属) , 1980 (1) :1-12.
[8] XUGuang xian (徐光宪) , LIBiao guo (李标国) .串级萃取理论发展的十五年[A].Procof2ndAnnualAca demicconfofCSRE (中国稀土学会第二届年会会议论文集) [C].Beijing, 1990 (1) :T1-T4.
[9] ZHONGCheng hua (钟盛华) .离子吸附型稀土萃取分离新工艺[J].TheChineseJournalofNonferrousMetals (中国有色金属学报) , 2000, 10 (2) :262-265.
[10] ZHONGCheng hua (钟盛华) , LIJi (李 机) .带支体工艺萃取法[P].CN 1044409A , 1990.
[11] ZHONGCheng hua (钟盛华) .Anewliquid liquidmix ersettlerextractor[A].LogsdailDH , SlaterMJeds.SolventExtractionintheProcessIndustries[C].Lon donandNewYork:ElsevierAppliedScience, 1993.111-118.
[12] ZHONGCheng hua.Uppersettlerextractor:Anewsolventextractionmixer settler[A].ShallcrossDC , PaiminR , PrvcicLM .ValueAddingthroughSolventExtraction[C].Australia, 1996, (2) :1221-1226.
[13] ZHONGCheng hua (钟盛华) .中混室多层澄清界面萃取器的研制[J].ChinJRareMetals (稀有金属) , 1996 (2) :100-104.
[14] XUGuang xian (徐光宪) .RareEarth (PartⅠ) (稀土, 上册) [M ].Beijing:MetallurgicalIndustryPress, 1995.709-720.