华电MATLAB上机实习报告

[CO2 ( aq ) ] [1
3
K C1 K C1K C 2 1 ] [ H ] [ H ]2
（20）
[ H ] K C 2 1 [ HCO ] [1 ] K C1 [ H ] [ H ] [ H ]2 1 [CO ] [1 ] K C 2 K C1 K C 2
原题：利用VB或C语言或Matlab 编制电厂锅炉给水化学工况的选 择程序。要求根据电厂关键水化 学基础数据以及关键水化学表现 （腐蚀、沉积、材质特点等）， 利用该程序自动做出给水化学工 况的判定。参考热力设备水汽理 化过程教材第46页。
2015-2-22
Matlab第五题
• • • • • • • • • • • • • • • • • • • • disp(a) A=input(x); if A==1 disp(b) else disp(c) B=input(x); if B==1 disp(d) disp(j) disp(f) C=input(x); if C==0 while C==0 disp(d) disp(e)
disp(g) D=input(x); while D==0 disp(d) disp(e) disp(i) disp(f) disp(g) D=input(x); end if D==1 disp(h) disp(i) disp(j) disp(k) end end end end
2015-2-22
在 CO2-H2O 体系中的存在 H2CO3 的一级和二级电离反应以及水的离解反应：
H 2O CO2 ( aq ) H HCO3
HCO3 H CO32
H 2O H OH
（7） （8） （9）
•
得到
2015-2-22
Matlab第四题
则 CO2, HCO3-, CO32-在碳酸盐组分总量 CT 中所占的比例分别为：
专业软件编制实习报告
四大板块
Word and Excel
• 文字部分的加工 • Excel数据计算 • 曲线拟合
Visio
• 图1 • 图2
Origin
• 图3
• 图4
Matlab
• 第三题 • 第四题 • 第五题
2
Visio图一
CO2加压 封存
再生氨水
Hale Waihona Puke Baidu洁净烟气
氨法
直流 电源
填充床 电解吸
过滤器 脱碳后溶液
解： 第一小题：p=input('输入CO2浓度p=');
Kh=3.36*10^(-7); Ka=10^(-6.3); Ca=Kh*p*(10^(-6))*1.013*10^5;
原题：利用VB或C语言或Matlab编制程序计算酸雨 的定义限值，并做出CO2浓度330~660ppm范围内 （步距为30）降水的CO2浓度~pH图。酸雨的定义 限值系通过如下给定参数求得：温度25℃，CO2浓 度330ppm, 其他常数查相关手册或环境化学教材。
H2CO3=100./(1+kc1./h+kc1*kc2./(h.*h));
HCO3=100./(1+h./kc1+kc2./h); CO3=100./(1+h./kc2+h.*h./(kc1.*kc2)); plot(pH,H2CO3,'r',pH,HCO3,'b',pH,CO3,'y'); title('碳酸化合物物种分布图') %设置题目 xlabel('pH');%x ylabel('分布律'); %y % legend('CO2+H2CO3','HCO3','CO3'); %设置图注
碳捕集
模拟烟气
泵
解吸环节
2015-2-22
吸收环节
氨吸收剂
Visio图二
再生氨水
（屏蔽开关） 产氨水室 NH3· H2O
HCO3-
CO2
（屏蔽开关） 产氨水室 NH3· H2O 产CO2室
HCO3
-
气液 分离槽
产CO2室
HCO3-
OHH2O OH-
H+
OHH2O OH-
H+
HCO3-
H2O
H2O
HCO3HCO3-
•
• • • • • • • •
C=sqrt(Ka*Ca);
ph=-log10(C) 第二小题： p=5*[66:132]; Kh=3.36*10^(-7); Ka=10^(-6.3); C=sqrt(Ka*Kh*p*10^(-6)*1.013*10^5); ph=-log10(C); plot(p,ph)
2015-2-22
Matlab第五题
• • • • else disp(g) D=input(x); while D==0
if D==1
disp(h) disp(i) disp(j) disp(k) end end else disp(e) disp(i) disp(f)
•
• •
disp(d)
disp(e) disp(i)
2015-2-22
第二小题运行结果
2015-2-22
Matlab第四题
• 利用VB或C语言或Matlab编制程序计算碳酸水溶液体系诸 化学形态的浓度分布，做出pH~形态分率α的图示，并输 出pH值在4~11范围内，pH步距为0.1，CO2, HCO3-, CO32各自的形态分数数值列表。 • 解：通过
2015-2-22
第四题运行结果
碳酸化合物物种分布图 100 90 80 70 60
分布律
50 40 30 20 10 0
4
5
6
7 pH
8
9
10
11
2015-2-22
Matlab第五题
• x='若是，输入1，若否，输入0 ';
•
• • • • • • • • • •
a='目前给水处理方式有无问题、铁铜含量是否较低';
b='按此方式继续运行'; c='机组是否为无铜系统'; d='优化AVT（O）（pH）'; e='优化AVT（R）（pH溶解氧）'; f='给水氢电导率小于0.15uS/cm,精处理系统运行正常。'; g='通过加氧试验后水汽系统的铜含量合格，铁含量较小。'; h='OT'; i='监测基础数据'; j='技术、经济性评价'; k='制定该锅炉给水AVT（O）、AVT（R）或OT规范';
C=input(x); if C==0 while C==0 disp(d) disp(e)
•
• •
disp(f)
disp(g) D=input(x);
disp(i) disp(f) C=input(x);
•
2015-2-22
end
end
Matlab第五题
• • • • • • • • • • • • • • • • • disp(g) D=input(x); while D==0 disp(d) disp(e) disp(i) disp(f) disp(g) D=input(x); end if D==1 disp(h) disp(i) disp(j) disp(k) end else
2 3
根据上式可绘制碳酸水溶液中各碳酸盐组分浓度分率随 pH 值的变化关系图
• 由此写出程序
2015-2-22
Matlab第四题
• • • • pH=[4:0.1:11]; h=10.^(-pH); kc1=10.^(-6.35); kc2=10.^(-10.33);
•
• • • • • • •
稀 溶 液 负极
HCO3
-
正极
HCO3-
HCO3 OH
NH4+ AEM
-
-
HCO3-
HCO3OH
-
HCO3-
H
+
HCO3-
NH4+ AEM BPM
HCO3-
BPM
AEM
脱碳后溶液
2015-2-22
Origin图一
2015-2-22
Origin图二
2015-2-22
Matlab第三题
• •
• • •
disp(g) D=input(x); while D==0 disp(d) disp(e) disp(i) disp(f) disp(g) D=input(x); end if D==1 disp(h) disp(i)
disp(i)
disp(j)
disp(f) C=input(x); end
disp(k) end