中国有色金属学报 2003,(04),1005-1010 DOI:10.19476/j.ysxb.1004.0609.2003.04.037
铝酸钠溶液碳酸化分解的热力学
李小斌 刘祥民 苟中入 彭志宏 刘桂华 周秋生 丁安平 李明 刘业翔
中南大学冶金科学与工程学院,中南大学冶金科学与工程学院,中南大学冶金科学与工程学院,中南大学冶金科学与工程学院,中南大学冶金科学与工程学院,中南大学冶金科学与工程学院,中国铝业股份有限公司山西分公司,中国铝业股份有限公司山西分公司,中南大学冶金科学与工程学院 长沙410083 ,长沙410083 ,长沙410083 ,长沙410083 ,长沙410083 ,长沙410083 ,河津043300 ,河津043300 ,长沙410083
摘 要:
为对高浓度铝酸钠溶液碳酸化分解过程中产品杂质含量的控制提供理论依据,分析了铝酸钠溶液碳酸化分解过程的热力学,认为该过程中Al(OH)3析出的真正机理是过饱和铝酸钠溶液的自发结晶;而基于分解过程中平衡浓度的热力学计算表明,在碳酸化分解过程中SiO2的变化规律与公认的三段变化规律一致,说明该过程产品中SiO2的含量取决于它在溶液中的平衡浓度,而非吸附所致;同时还分析了丝钠铝石的形成热力学。
关键词:
铝酸钠溶液 ;碳酸化分解 ;热力学 ;机理 ;
中图分类号: TF821
作者简介: 李小斌(1962),男,教授,博士生导师,博士,电话:07318830454;Email:X.B.Li@mail.csu.edu.cn;
收稿日期: 2002-07-29
基金: 国家重点基础研究发展规划资助项目(G1999064910);
Thermodynamics of carbonization of aluminate solution
Abstract:
The mechanism of carbonization process of aluminate solution was analyzed thermodynamically. The results show that the internal reason of precipitation of Al(OH)3 in the process is its spontaneous nucleation caused by supersaturation of aluminate solution. The behaviors of SiO2 in the process are in accordance with that from published experiment, which implies that the content of product SiO2 from carbonization depends on its solubility but is not the result of absorption. The thermodynamic calculation for forming Na2O·Al2O3·2CO2·nH2O was achieved.
Keyword:
aluminate solution; carbonization; thermodynamics; mechanism;
Received: 2002-07-29
烧结法氧化铝生产流程中, 碳酸化分解过程是包含液-固-气多相反应的复杂过程。 对该过程的机理解释存在多种不同的观点。
1)氢氧化铝的析出机制 认为可能存在两种机制, 即溶液稳定性破坏机制和酸碱中和直接作用机制。 所谓溶液稳定性的破坏机制, 即认为铝酸钠溶液中的游离苛性碱被通入的CO2 按反应(1)所中和, 导致体系苛性比降低, 溶液过饱和度增大, 其稳定性降低, 然后按反应式(2)析出Al(OH)3 沉淀。 该过程可用反应式(3)来描述。
2OH- +CO2 =CO
2 ? 3
3
2
-
+H2 O (1)
Al(OH)- 4 =Al(OH)3 +OH- (2)
Al(OH)- 4 +OH- +CO2 =Al(OH)3 +CO
2 ? 3
3
2
-
+
H2 O (3)
而酸碱中和直接作用机制则认为Al(OH)3 的析出是碱性的NaAl(OH)4 与酸性的CO2 直接作用的结果, 铝酸钠溶液碳酸化分解发生的反应为
2Al(OH)- 4 +CO2 =2Al(OH)3 +CO
2 ? 3
3
2
-
+H2 O (4)
2)SiO2 在碳分过程中的行为 普遍认为
[1 ]
, 在碳分过程中, SiO2 的析出分3个阶段: 第一阶段发生在分解初期, 溶液中SiO2 浓度下降较快; 第二阶段则主要析出Al(OH)3 , 溶液中SiO2 浓度变化较小; 第三阶段则由于溶液中苛性碱浓度的降低, SiO2 在溶液中平衡浓度大幅减小而基本上全部析出, 溶液中SiO2 浓度下降较快, 此阶段发生的反应如反应式(5)所示。
2Na+ +2Al(OH)- 4 +2SiO2 (OH)
2 ? 2
2
2
-
=
Na2 O·Al2 O3 ·2SiO2 ·2H2 O+4OH- +
2H2 O (5)
3)丝钠铝石的形成 丝钠铝石的形成主要是HCO- 3 与Al(OH)- 4 作用的结果。 在碳分末期或体系局部过度碳酸化时, 可能发生如式(6)~(8)的反应。 但在碳分初期, 由于游离苛性碱的存在, 丝钠铝石可以按反应(9)被其溶解。 但对丝钠铝石的形成条件,文献
[
2 ,
3 ]
报道的结果并不一致。
2Al(OH)- 4 +2Na+ +4HCO- 3 =
Na2 O·Al2 O3 ·2CO2 ·n H2 O+
2CO
2 ? 3
3
2
-
+(6-n )H2 O (6)
2Al(OH)- 4 +2Na+ +2CO- 3 =
Na2 O·Al2 O3 ·2CO2 ·n H2 O+
4OH- +(2-n )H2 O (7)
2Al(OH)- 4 +2Na+ +2HCO- 3 =
Na2 O·Al2 O3 ·2CO2 ·n H2 O+
2OH- +(4-n )H2 O (8)
Na2 O·Al2 O3 ·2CO2 ·n H2 O+4OH- =
2Al(OH)- 4 +2CO
2 ? 3
3
2
-
+2Na+ +(n -2)H2 O (9)
上述观点主要是针对较低浓度的铝酸钠溶液碳酸化分解过程提出的, 很少涉及高浓度铝酸钠溶液。 而烧结法一个主要的发展方向是大幅度提高原液中Al2 O3 浓度以提高溶液的产出率。 本研究在实验测定分解过程苛性碱浓度变化的基础上, 对碳酸化分解过程进行了热力学分析, 并从理论上加以解释, 以期实现高浓度铝酸钠溶液的碳酸化分解工艺优化。
1 实验
1.1 实验方法
以工业纯Al(OH)3 、 NaOH、 Na2 CO3 配制成铝酸钠溶液, 工业熟料溶出后得到含一定浓度的SiO2 粗液, 将两种溶液配制成所要求的分解原液。 准确移取2.0 L原液至反应槽中, 密封后采用JBV-Ⅲ变频调速搅拌器搅拌, 分别用CYES-ⅡO2 /CO2 气体测定仪和LZB-6玻璃转子流量计控制CO2 气体浓度和流量, 控温精度为±1 ℃。 在规定时间取样, 溶液经离心分离后分析相关组分的浓度。
1.2 实验结果
表1列出了高浓度铝酸钠溶液碳酸化分解母液中Al2 O3 和Na2 Ok 浓度变化和分解产物中Na2 O含量变化情况, 图1所示为根据表1计算的分解过程母液苛性比的变化。
由图1可见, 高浓度铝酸钠溶液碳酸化分解中溶液苛性比随分解时间的变化均近似“W”型曲线, 隐含了氧化铝和苛性碱等物质浓度以及所涉及化学反应的反应速度等大量相关信息。
2 热力学分析
2.1 Al(OH)3析出的热力学分析
为了确定碳酸化分解过程Al(OH)3 的析出机理, 结合实验条件和结果(表1)计算了反应(3)、
(4)的反应自由能。 由于碳酸化分解的反应均在非标准状态下进行, 溶液中主要物质的浓度高且变化较快, 按标准状态计算的结果与实际状态存在很大差距。 因此, 在进行Al(OH)3 析出的热力学计算时, 离子的活度系数采用Robinson等修正的Debye-hückel方程式计算
[4 ,5 ,6 ]
:
-lgf i =0.511I /(1+1.50I )-0.20I (10)
式中 I 为离子强度,
I = ∑ i x i z 2 i / 2 ( x i
Ι
=
∑
i
x
i
z
i
2
/
2
(
x
i
为离子i 的摩尔分数); f i 为离子i 的活度系数。
而化学反应Gibbs自由能可以通过式(11)计算
[7 ]
:
ΔG ? T =ΔG ? 298 -(T -298)ΔS ? 298 +
T ∫
T 2 9 8
2
9
8
Τ
Δ cp d T/T+∫
T 2 9 8
2
9
8
Τ
Δ cp d T (11)
cp =a+b×10-3 T+c×105 T-2 (12)
Δ GT =Δ G? T +RTln KT (13)
式中 a、 b、 c分别为热容系数, KT 为反应平衡常数。 反应(3)和(4)的KT 分别为:
图1 碳酸化分解过程中溶液αk随时间的变化
Fig.1 Changes of caustic ratio in carbonization process of pregnant solutions
表1 高浓度铝酸钠溶液碳酸化分解溶液的组成变化
Table 1 Content change in carbonization of high concentration aluminate solution(g/L)
Time/min
Exp.1
Exp.2
Exp.3
Exp.4
Al2 O3
Na2 Ok
Al2 O3
Na2 Ok
Al2 O3
Na2 Ok
Al2 O3
Na2 Ok
0
170.91
163.88
183.66
168.58
187.50
176.24
192.50
183.28
30
163.89
136.89
177.90
140.39
176.50
131.59
174.50
136.29
60
99.12
89.29
132.50
112.79
144.37
110.07
128.45
112.40
80
60.68
68.14
100.81
98.69
115.44
96.34
99.64
96.34
100
33.88
47.10
64.26
68.14
75.33
70.49
61.79
68.14
115
16.45
21.15
45.92
51.96
42.12
44.65
29.33
39.95
130
6.67
9.40
37.46
39.95
21.69
28.20
19.43
21.15
145
-
-
20.04
14.90
12.47
7.64
8.40
9.94
160
-
-
10.70
9.40
2.70
3.40
2.10
4.10
Notes: temperature 85 ℃; w (CO2 )=34%~37%; seed ratio(K r )=0(Exp.1 and Exp.2), 0.2 (Exp.3), 0.12 (Exp.4).
K T (3) =x (CO
2 ? 3
3
2
-
)·F (CO
2 ? 3
3
2
-
)·p? /
[x(Al (OH )- 4 )·f(Al (OH )- 4 )·
x(OH - )·f(OH - )·p(CO 2 )] (14)
KT(4) =x(CO
2 ? 3
3
2
-
)·F(CO
2 ? 3
3
2
-
)·p? /
[x2 (Al (OH )- 4 )·f2 (Al (OH )- 4 )·
p(CO 2 )] (15)
实验数据的计算在“氧化铝生产热力学数据库软件”
[8 ]
上完成, 计算结果如图2, 3所示。 从图2和3可见, 在碳酸化分解工艺条件下, 析出氢氧化铝的两个反应的ΔG T 均小于-18 kJ/mol。 但反应(3)的ΔG T 在-200 kJ/mol以下, 说明从热力学上考虑, 碳酸化分解过程中Al(OH)3 的析出主要原因是由过饱和度的升高而产生的自发分解为主, 酸碱中和直接作用析出Al(OH)3 是次要的。 根据这一结果, 为在高浓度铝酸钠溶液碳酸化分解过程中控制产品碱含量而制订的相应措施能取得良好的实用效果。
图2 反应式(3)在不同温度下的反应自由能ΔGT
Fig.2 Gibbs free energy changes of reaction(3) at different temperatures 1—ρ (Al2 O3 )=10.70 g/L; 2—ρ (Al2 O3 )=20.00 g/L; 3—ρ (Al2 O3 )=100.80 g/L; 4—ρ (Al2 O3 )=183.66 g/L
图3 反应式(4)在不同温度下的反应自由能ΔGT
Fig.3 Gibbs free energy changes of reaction(4) at different temperatures 1—ρ (Al2 O3 )=10.70 g/L; 2—ρ (Al2 O3 )=20.00 g/L; 3—ρ (Al2 O3 )=100.80 g/L; 4—ρ (Al2 O3 )=183.66 g/L
2.2SiO2在碳酸化分解过程中的热力学分析
文献
[
9 ,
10 ,
11 ,
12 ,
13 ,
14 ]
研究了SiO2 在铝酸钠溶液中的平衡浓度、 水合铝硅酸钠的析出动力学及其影响因素。 但对碳酸化分解过程中水合铝硅酸钠形成的热力学研究未见报道, 可能是因为缺乏相应的热力学数据。
文献
[
15 ,
16 ,
17 ]
计算了Na2 O-Al2 O3 -SiO2 -H2 O系中各物质的ΔG ? 298 和恒压热容系数a 、 b 、 c 的值。 根据Na2 O·Al2 O3 ·2SiO2 ·2H2 O的吉布斯自由能ΔG ? 298 的计算值, 结合反应式(5), 可计算铝酸钠溶液中SiO2 的平衡浓度。 由于铝酸钠溶液结构极其复杂, 离子种类及结构至今不明确, 为简化计算, 若以离子浓度代替活度, 则反应(5)的反应平衡常数为
K T = [ O H ? ] 4 [ N a + ] 2 [ A l ( O H ) ? 4 ] 2 [ S i O 2 ( O H ) 2 ? 2 ] 2 ? ? ? ( 1 6 )
Κ
Τ
=
[
Ο
Η
-
]
4
[
Ν
a
+
]
2
[
A
l
(
Ο
Η
)
4
-
]
2
[
S
i
Ο
2
(
Ο
Η
)
2
2
-
]
2
?
?
?
(
1
6
)
而α k -1=[OH- ]/[Al(OH)- 4 ], 因此铝酸钠溶液中SiO2 (OH)
2 ? 2
2
2
-
的浓度为
[ S i O 2 ( O H ) 2 ? 2 ] = ( α k ? 1 ) [ O H ? ] [ N a + ] ? exp ( Δ G ? T / ( 2 R T ) ) ? ? ? ( 1 7 )
[
S
i
Ο
2
(
Ο
Η
)
2
2
-
]
=
(
α
k
-
1
)
[
Ο
Η
-
]
[
Ν
a
+
]
?
exp
(
Δ
G
Τ
?
/
(
2
R
Τ
)
)
?
?
?
(
1
7
)
由式(17)表示的SiO2 (OH)
2 ? 2
2
2
-
浓度可以转化成SiO2 的浓度。 在不同温度下铝酸钠溶液中SiO2 的平衡质量浓度随碳酸化分解的分解率变化的计算结果如图4所示。
图4 碳酸化分解溶液中SiO2平衡质量浓度的变化
Fig.4 Changes of equilibrium mass concentration of silica in spent solution of carbonization with precipitation rate
由图4和式(17)可以看出, 高浓度铝酸钠溶液碳酸化分解过程中SiO2 的平衡浓度与溶液温度、 苛性比、 OH- 浓度和Na+ 浓度等有关, 并且随分解率变化存在较明显的阶段性特征, 与文献报道的实验研究结果完全相符, 说明所进行的热力学分析是正确的。
图1和图4表明: 在分解初期, 苛性碱与CO2 之间的中和反应速度较快, OH- 浓度降低, 虽然Al2 O3 的过饱和度增大, 但Al(OH)3 析出并不显著, 因此溶液的α k 和OH- 浓度迅速降低, 导致SiO2 的平衡浓度快速下降; 分解中期, 铝酸钠溶液因稳定性急剧下降而快速分解, 按反应式(2)析出Al(OH)3 , 并产生OH- , 游离苛性碱得到一定程度的补充, α k 逐步回升, 但Na+ 浓度因浓缩也有所增加, 最后导致溶液中SiO2 的平衡浓度变化不大; 而后, 随着分解过程的进行, Al(OH)3 析出量减少, 反应(2)中OH- 生成量比反应(1)中OH- 的消耗量要少得多, α k 再次降低, 且Na+ 浓度因溶液浓缩而进一步升高, SiO2 的平衡浓度急剧下降。 图1中分解末期α k 的第二次增大现象是因为该阶段减少甚至停止通气, 由反应(2)造成的。
这种规律与文献
[
3 ]
报道的实验结果完全一致。 而文献
[
3 ]
对此的解释是: 分解初期, 母液中SiO2 浓度的下降是由于新析出的Al(OH)3 表面活性大, 吸附了溶液中的硅酸根离子所致; 而分解中期则是因为随着分解过程的进行, Al(OH)3 表面钝化, 吸附硅酸根离子的能力减小, 只有在分解末期, 才因为硅酸根在溶液中平衡浓度的减小而以水合铝硅酸钠的形式析出。 但从热力学的角度, 可以认为SiO2 在分解过程中所表现出的阶段性特征的真正原因应该是其平衡浓度的变化。
2.3 丝钠铝石形成的热力学分析
在铝酸钠溶液碳酸化分解中, 丝钠铝石(Na2 O·Al2 O3 ·2CO2 ·n H2 O)的析出将严重影响碳酸化分解产品的质量。 因此, 对它的形成进行热力学分析是很有意义的。
按文献
[
8 ]
的拆分和计算方法, 计算出在343~373 K内丝钠铝石的生成自由能ΔG ? T , 计算结果可以经过线性回归成ΔG ? T =A +BT 形式:
Δ G ? T = ? 3 3 2 0 . 5 ? 0 . 6 9 5 1 T ? ( n = 1 ) Δ G ? T = ? 3 6 0 9 . 5 ? 0 . 7 5 7 6 T ? ( n = 2 ) } ? ? ? ( 1 8 )
Δ
G
Τ
?
=
-
3
3
2
0
.
5
-
0
.
6
9
5
1
Τ
?
(
n
=
1
)
Δ
G
Τ
?
=
-
3
6
0
9
.
5
-
0
.
7
5
7
6
Τ
?
(
n
=
2
)
}
?
?
?
(
1
8
)
然后, 根据表1所列Exp.2的数据, 按式(10)计算出的实际分解条件下参与反应离子的活度系数, 再由式(11)和(12), 结合式(13)和(19)计算出反应(6)的ΔG ? T 和ΔG T , 其中反应平衡常数K T 为
K T =x 2 (CO
2 ? 3
3
2
-
)·f 2 (CO
2 ? 3
3
2
-
)/
[x 2 (Al(OH)- 4 )·f 2 (Al(OH)- 4 )·
x 4 (HCO- 3 )·f 4 (HCO- 3 )·
x 2 (Na+ )·f 2 (Na+ )] (19)
以同样方法可计算出反应(7)和(8)的ΔG ? T 和ΔG T , 结果分别如图5~7所示。
图5 不同温度下按反应(6)析出 Na2O·Al2O3·2CO2·nH2O的 ΔGT与Al2O3质量浓度的关系
Fig.5 Relationship between ΔG T and mass concentration of Al2 O3 for reaction (6)
图6 不同温度下按反应(7)析出 Na2O·Al2O3·2CO2·nH2O的 ΔGT与Al2O3质量浓度的关系
Fig.6 Relationship between ΔG T and mass concentration of Al2 O3 for reaction (7)
图7 不同温度下按反应(8)析出 Na2O·Al2O3·2CO2·nH2O的 ΔGT与Al2O3质量浓度的关系
Fig.7 Relationship between ΔG T and mass concentration of Al2 O3 for reaction (8)
从图5~7可以看出, 温度、 分解深度(Al2 O3 浓度)、 丝钠铝石中化合水数量以及反应途径都是影响丝钠铝石析出的因素。 温度升高时, 式(6)~(8)3个反应的反应自由能ΔG T 都增大, 即在热力学上提高分解温度有利于减少该物质析出, 这一结论与文献
[
2 ]
的实验结论一致, 但是同文献
[
3 ]
相悖, 这可能与丝钠铝石析出动力学有关; 分解深度越高(母液中Al2 O3 浓度越小), 丝钠铝石越容易形成; 在碳酸化分解过程中, 由HCO- 3 比由CO
2 ? 3
3
2
-
更易形成丝钠铝石; 析出丝钠铝石的自由能还与化合水数量密切相关, n =1的丝钠铝石比n =2的更难析出。 但n 与分解温度有关, 一般来说, 温度越高, n 越小, 结合热力学计算和生产实践经验, 工业分解条件下n 应该为1。
就形成Na2 O·Al2 O3 ·2CO2 ·H2 O而言, 从图5~7可以看出, 在实际分解温度下(80~90 ℃), 由于反应(6)、 (7)的ΔG T >-20 kJ/mol, 而反应(8)的ΔG T <-96 kJ/mol, 从而热力学上可以认为, 该反应可能是丝钠铝石在分解过程中析出的主要途径。
以上分析表明, 在进行碳酸化分解作业时, 为了避免丝钠铝石析出而污染产品, 应该控制较高的分解温度, 同时还要控制分解深度和加强搅拌, 以保证体系Al2 O3 浓度均匀和避免过度碳酸化而产生大量HCO- 3 。
3 结论
1) 在碳酸化分解工艺条件下, 溶液的苛性比值随分解时间呈一条“W”字型曲线, 意味着该过程是OH- 被CO2 中和以及NaAl(OH)4 水解析出两个反应的共同作用。 热力学计算也证明了铝酸钠溶液的碳酸化分解的实质就是这种共同作用, 而不是CO2 与NaAl(OH)4 直接作用的结果。
2) 热力学分析表明, 铝酸钠溶液中SiO2 的平衡浓度与溶液温度、 苛性比、 苛性碱浓度和Na+ 浓度等有关, SiO2 在分解过程中的行为主要取决于它在溶液中的平衡浓度, 而非吸附作用。
3) 分解温度、 体系中Al2 O3 和HCO- 3 浓度是影响丝钠铝石析出的重要因素, 且丝钠铝石主要通过反应式(8)形成。
参考文献
[1] 杨重愚.轻金属冶金学[M].北京:冶金工业出版社,1996.5.YANGZhong yu.LightMetalsMetallurgy[M].Beijing:MetallurgicalIndustryPress,1996.5.
[2] HanamirovaAA.Aluminaandrouteofdecreasingitsimpurities[M].Armenia:ArmeniaAcademicPress(inSovietUnion),1983.116158.(inRussian)
[3] 杨重愚.氧化铝生产工艺学(修订版)[M].北京:冶金工业出版社,1993.YANGZhong yu.TechnologyofAluminaProduction(revised)[M].Beijing:MetallurgicalIndustryPress,1993.
[4] DeanJA.兰氏化学手册[M].第3版.尚方久译.北京:科学出版社,19991.329.DeanJA.LANGE SHandbookofChemistry[M].3rded.SHANGFang jiutransl.Beijing:SciencePress,1991.3
[5] 黄子卿.电解质溶液理论导论(修订版)[M].北京:科学出版社,1983.818.HUANGZi qing.IntroductionofElectrolyteTheory(revised)[M].Beijing:SciencePress,1983.818.
[6] 李以圭.金属溶剂萃取热力学[M].北京:清华大学出版社,1988.5.LIYi gui.ExtractionThermodynamicofMetal ioninSolution[M].Beijing:TsinghuaUniversityPress,1988.5.
[7] 傅崇说.有色冶金原理[M].北京:冶金工业出版社,1993.154157.FUChong shuo.PrincipleofNonferrousMetallurgy[M].Beijing:MetallurgicalIndustryPress,1993.154157.
[8] 李永芳.氧化铝生产热力学数据库的研究与开发[D].长沙:中南大学,2001.LIYong fang.ResearchonThermodynamicDataforAluminaProduction[D].Changsha:CentralSouthUniversity,2001.
[9] OkuT,YamadaK.ThedissolutionratioofquartzandtherateofdesilicationinBayerliquor[J].LightMetals,1971:
[10] AdamsonAN,BlooreEJ,CarrAR.BasicprinciplesofBayerprocessdesign[J].JExtrMetofAluminum,1963,1:2325.
[11] MurakamiM.Aspectofdesilicationkineticswithrecycledseed[J].LightMetals,1992:113.
[12] CresswellJ.FactorsaffectingdesilicationofBayerprocessliquor[J].Chemeca,1984(48):285.
[13] JamialahmadiM,SteinhagenHM.DeterminingsilicasolubilityinBayerprocessliquor[J].JOM,1998,11:4449.
[14] GersonAR,ZhengK,etal.TheprecipitationmechanismofsodiumaluminosilicatescaleinBayerplants[J].LightMetals,1999:2828.
[15] 刘桂华,李小斌,李永芳,等.复杂无机化合物组成与热力学数据间的线性关系及其初步应用[J].科学通报,2000,45(13):13861392.LIUGui hua,LIXiao bin,LIYong fang,etal.Thelinearrelationshipofcomponentsofcomplicatecompoundandther modynamicdataanditsapplication[J].ChineseScienceBulletin,2000,45(13):13861392.
[16] 李小斌,李永芳,刘桂华,等.复杂化合物硅酸盐吉布斯自由能和焓的一种近似算法[J].硅酸盐学报,2001,3:4549.LIXiao bin,LIYong fang,LIUGui hua,etal.Acalcula tionmethodofGibbsfreeenergyofcomplicatecompoundcontainingsilica[J].JournalofCeramicSociety,2001,3:4549.
[17] 李小斌,刘桂华,彭志宏,等.铝酸钠溶液中二氧化硅平衡浓度的研究[J].东北大学学报,2002,23(3):251254.LIXiao bin,LIUGui hua,PENGZhi hong,etal.Equilib riumconcentrationofsilicainaluminatesolution[J].JournalofNortheasternUniversity,2002,23(3):251254.