1.一种基于数学模型指导的粗铜吹炼工艺，其特征在于：包括以下步骤：
①先将冰铜、烟尘、石灰和石英沙混合均匀后置入闪速吹炼炉的冰铜喷嘴，再将炉料与工艺风由喷嘴一起置入吹炼炉的反应塔内迅速完成吹炼反应，为改变反应塔中间区域的氧势，利用反应塔内安装的中央氧枪吹入氧气，并向反应塔内通入压缩空气，通过压缩空气的气流改善炉料在反应塔内的分布状况；炉内的氧化程度和产出物的温度分别由氧/料比和工艺风的氧浓度控制，当炉内热量不足时，通过冰铜喷嘴周围的天然气烧咀补充热量，当热量过剩的情况下，则加入烟尘调节炉内的热平衡状况；
②上步反应生成的熔体在沉淀池内分离为粗铜和渣，产出的粗铜通过溜槽定期进入阳极炉进行精炼，吹炼渣经水淬磨细脱水后返回闪速熔炼炉继续熔炼，吹炼炉的烟气经余热锅炉和电收尘后送制酸系统制酸，烟气经余热锅炉沉降下来的烟尘和电收尘的烟尘经破碎筛分后返回闪速熔炼，电收尘的其他烟尘则返回吹炼炉作冷料；
③根据冶金学关于在密闭环境中火法冶金过程化学平衡和热力学平衡理论，建立以金属平衡和热平衡为基础的闪速炉粗铜吹炼数学模型；通过金属平衡建立一个11元线性方程组，以求取当粗铜含硫达到目标值时物料在炉内化学反应所需要的氧量，渣中四氧化三铁、硅品位达到目标值时的钙量和硅量；热平衡是求取渣温达到目标值时的氧气浓度或天然气量，模型既能把氧气浓度作为未知数，也能将天然气量作为未知数；
利用该数学模型，根据给定的目标参数和生产条件能够精确计算出需要的最佳操作参数，使粗铜吹炼工艺过程的关键控制指标：粗铜硫品位、渣中四氧化三铁、渣温保持稳定性，按所得的参数进行生产，具体计算的相关表格与计算公式如下：
一、金属平衡
(一)、金属平衡表
根据控制要求和投入产出条件，构造如下金属平衡表
备注：表中 G：给定值 C：解方程前计算 Xn：未知数 空白：解方程后计算--：不计算
(二)、金属平衡方程式
根据金属平衡原理，将金属平衡表中的未知量，构造成一个11元线性方程组，解此方程组，即可求得为达到粗铜吹炼目标参数所需的氧系数、石灰量、硅酸矿量的操作值，具体方程如下：
1、物料平衡
冰铜量+投入C烟尘量+石灰量+硅酸矿量+O量
＝粗铜量+渣量+产出B烟尘量+产出C烟尘量+烟气中的S量+烟气中的其它量冰铜量+投入C烟尘量+X1+X2+X11
＝X3+X4+X7+X8+X9+X10
2、Cu平衡
冰铜Cu量+投入C烟尘Cu量+石灰Cu量+硅酸矿Cu量
＝粗铜Cu量+渣中的Cu量+产出B烟尘Cu量+产出C烟尘Cu量冰铜Cu量+投入
C烟尘Cu量+X1*石灰Cu品位+X2*硅酸矿Cu品位
＝X3*粗铜Cu品位+X4*渣中的Cu品位+X7*产出B烟尘Cu品位+X8*产出C烟尘Cu品位
3、Fe平衡
冰铜Fe量+投入C烟尘Fe量+石灰Fe量+硅酸矿Fe量
＝粗铜Fe量+渣中的Fe量+产出B烟尘Fe量+产出C烟尘Fe量冰铜Fe量+投入
C烟尘Fe量+X1*石灰Fe品位+X2*硅酸矿Fe品位
＝X3*粗铜Fe品位+X5+X7*产出B烟尘Fe品位+X8*产出C烟尘Fe品位
4、S平衡
冰铜S量+投入C烟尘S量+石灰S量+硅酸矿S量
＝粗铜S量+渣中的S量+产出B烟尘S量+产出C烟尘S量+烟气中的S量冰铜S
量+投入C烟尘S量+X1*石灰S品位+X2*硅酸矿S品位
＝X3*粗铜S品位+X4*渣中的S品位+X7*产出B烟尘S品位+X8*产出C烟尘S品位+X9
5、SiO2平衡
冰铜SiO2量+投入C烟尘SiO2量+石灰SiO2量+硅酸矿SiO2量
＝粗铜SiO2量+渣中的SiO2量+产出B烟尘SiO2量+产出C烟尘SiO2量冰铜SiO2量+投入C烟尘SiO2量+X1*石灰SiO2品位+X2*硅酸矿SiO2品位
＝X3*粗铜SiO2品位+X4*渣中的SiO2品位+X7*产出B烟尘SiO2品位+X8*产出C烟尘SiO2品位
6、CaO平衡
冰铜CaO量+投入C烟尘CaO量+石灰CaO量+硅酸矿CaO量
＝粗铜CaO量+渣中的CaO量+产出B烟尘CaO量+产出C烟尘CaO量冰铜CaO量+投入C烟尘CaO量+X1*石灰CaO品位+X2*硅酸矿CaO品位
＝X3*粗铜CaO品位+X6+X7*产出B烟尘CaO品位+X8*产出C烟尘CaO品位
7、渣中的CaO/Fe平衡
X6＝X5*目标CaO/Fe
8、产出B烟尘量
X7＝(冰铜量+投入C烟尘量+X1+X2)*烟尘发生率*(1-C烟尘率)
9、产出C烟尘量
X8＝(冰铜量+投入C烟尘量+X1+X2)*烟尘发生率*C烟尘率
10、烟气中的其它量平衡
X10＝[(冰铜其它量+投入C烟尘其它量+X1*石灰其它品位+X2*硅酸矿其它品位)-(冰铜O2量+投入C烟尘O2量+硅酸矿O2量)]
*(1-其它系数)
其中：冰铜O2量＝Fe3O4的O2量+2FeO.SiO2的O2量+Cu2O的O2量
＝冰铜Fe3O4的Fe量*K7+冰铜2FeO.SiO2的Fe量*K14+冰铜Cu2O的Cu量*K16投入C烟尘O2量＝Cu2O的O2量+CuSO4的O2量+Fe3O4的O2量
＝投入C烟尘Cu2O的Cu量*K16+[FCF C烟尘的S量-Cu2S的Cu量*K1]*K17+FCF C烟尘Fe量*K7
硅酸矿O2量＝硅酸矿的Fe量*K11
11、O2平衡
X11＝(粗铜O2量+渣中的O2量+产出C烟尘的O2量+产出B烟尘的O2量)-(冰铜O2量+投入C烟尘O2量+硅酸矿O2量)
其中：粗铜O2量＝Cu2O的O量+Fe3O4的O量
＝粗铜Cu2O的Cu量*K16+粗铜中的Fe量*K7
渣中的O2量＝Cu2O的O2量+CuO的O2量+Fe2O3的O2量+Fe3O4的O2量+FeOX.CaO的O2量+2FeO.SiO2的O2量+CaSO4的O2量
＝渣中的Cu2O的Cu量*K16+渣中的CuO的Cu量*K8+渣中的Fe2O3的Fe量*K11+渣中的Fe3O4的Fe量*K7+渣中的FeOX.CaO的Fe量*K14+渣中的2FeO.SiO2的Fe量*K14+渣量*P10*K12*K15
产出B烟尘的O2量＝Cu2O的O量+Fe3O4的O量
＝[FCF B烟尘的Cu量-Cu2S的Cu量]*K16+FCF B烟尘Fe量*K7
其它各项的O2量见烟气中的其它量平衡，
计算常数
(三)氧气量计算
3 3
1、理论氧气量＝(烟气中的S量*22.4/32.06+投入产出的O2量*22.4/32)*10Nm，
100％O2
2、必要氧气量＝理论氧气量/反应塔氧效率；
二、热平衡
(一)热平衡表
1、反应塔热平衡
2、沉淀池热平衡
3、上升烟道热平衡
4、总热平衡
(二)热平衡计算
A、入热
1、S燃烧热＝烟气中的S量*反应塔S的氧化率*SO2生成热
2、FeO的生成热＝(渣中的Fe量-渣中FeS的Fe量-渣中Fe3O4的Fe量)*FeO反应热
其中：FeS的Fe量＝{渣中的S量*P30-[渣中的Cu量-Cu2O的Cu量-CuO的Cu量-金属Cu量]*k1-渣量*P10*K12}*P2
Fe3O4的Fe量＝渣中的Fe量*P9
3、Fe3O4含CuSO4、CaSO4生成热＝Fe3O4生成热+CuSO4生成热+CaSO4生成热其中：
·Fe3O4生成热＝(渣中Fe3O4的Fe量+产出C烟尘Fe3O4的Fe量+产出B烟尘Fe3O4的Fe量+粗铜Fe3O4的Fe量)*Fe3O4反应热
其中：渣中Fe3O4的Fe量＝渣中的Fe量*P9
产出C烟尘Fe3O4的Fe量＝产出C烟尘Fe量
B烟尘Fe3O4的Fe量＝FCF B烟尘Fe量
粗铜Fe3O4的Fe量＝粗铜中的Fe量
·CuSO4生成热＝产出C烟尘CuSO4的Cu量*CuSO4反应热
其中：CuSO4的Cu量＝(C烟尘的Cu量-C烟尘Cu2S的Cu量-C烟尘Cu2O的Cu量)·CaSO4生成热＝渣中CaSO4的Ca量*CaSO4反应热
其中：渣中CaSO4的CaO量＝渣量*P10*(40.08+16)/(40.08+32.06+16*4)
4、分解热＝-(Fe3O4分解热+投入C烟尘CuSO4分解热+FeS的分解热+Cu2S变成CuSO4分解热+硅酸矿分解热+石灰分解热)
其中：
·Fe3O4分解热＝(冰铜Fe3O4的Fe量+投入C烟尘的Fe3O4的Fe量)*Fe3O4反应热其中：冰铜Fe3O4的Fe量
＝冰铜Fe量-冰铜FeS的Fe量-冰铜2FeO.SiO2的Fe量投入C烟尘Fe3O4的Fe量＝投入C烟尘Fe量
·投入C烟尘CuSO4分解热＝投入C烟尘CuSO4的Cu量*CuSO4反应热
·FeS的分解热＝[(粗铜的Fe量+渣中Fe量+产出C烟尘的Fe量+B烟尘的Fe
量)-(冰铜Fe3O4的Fe量+冰铜2FeO.SiO2的Fe量+投入C烟尘Fe3O4的Fe量+硅酸矿的Fe量)]*FeS反应热
·Cu2S变成CuSO4分解热＝(产出C烟尘的Cu量-投入C烟尘的Cu量)*Cu2S反应热·硅酸矿分解热＝硅酸矿的Fe量*Fe2O3反应热
·石灰分解热＝CaCO3量*CaCO3反应热+MgCO3量*MgCO3反应热
5、炉料显热
[总装入量*炉料比热+(冰铜量*冰铜水分率+石灰量*石灰水分率+硅酸矿量*硅酸矿水分率)]*(炉料温度-基准温度)
6、沉淀池造渣热＝渣中2FeO.SiO2的Fe量*2FeO.SiO2生成热+渣中FeOX.CaO的Fe量*FeOX.CaO生成热
7、反应塔燃料燃烧热＝反应塔天然气量*天然气发热量
8、反应塔送风显热＝(O2量*O2单位热量+N2量*N2单位热量+H2O量*H2O单位热-3
量)*10 Mj/h
其中：单位热量计算通式：
2 -3 -1 3 3
单位热量＝(a*TG+b*TG*10 +c*TG *10+d)*1/22.4*4.187Kj/Nm
系数表：
成分 a b c d
SO2 10.38 1.27 1.42 -3683
N2 6.66 0.51 0 -2030
H2O 7.17 1.28 -0.08 -2223
TG：反应塔送风温度+217.16K
O2量＝反应塔工艺风O2量+反应塔燃烧风O2量
反应塔工艺风O2量＝必要矿石氧气量
反应塔燃烧风O2量＝反应塔燃烧风量*反应塔燃烧风氧浓度
N2量＝反应塔工艺风N2量+反应塔燃烧风N2量+反应塔天然气量*天然气含N2率
反应塔工艺风N2量＝必要矿石氧气量/工艺风富氧浓度*(1-工艺风富氧浓度)
反应塔燃烧风N2量＝反应塔燃烧风量湿基*(1-反应塔燃烧风氧浓度)H2O量＝反应塔工艺风H2O量+反应塔燃烧风H2O量
反应塔工艺风H2O量＝反应塔工艺空气量干基*(空气水分率)
反应塔燃烧风H2O量＝反应塔燃烧空气量*反应塔空气水分率
9、沉淀池S燃烧热＝烟气中的S量*(1-反应塔S的氧化率)*SO2生成热
10、沉淀池天然气燃烧热＝沉淀池天然气量*天然气发热量
11、沉淀池燃烧风显热＝(O2量*O2单位热量+N2量*N2单位热量+H2O量*H2O单位热-3
量)*1O O2量、N2量、H2O量及其单位热量参见第8项，
12、上升烟道天然气燃烧热＝上升烟道天然气量*天然气发热量
13、上升烟道燃烧风显热＝(O2量*O2单位热量+N2量*N2单位热量+H2O量*H2O单位热-3
量)*1O
O2量、N2量、H2O量及其单位热量参见第8项，
14、入热合计＝上述1～13项之和
常数、系数
B、出热
1、粗铜显热＝粗铜量*粗铜单位热量
粗铜单位热量＝0.1766*TB+5.29-(0.0885+0.0005*TB)*粗铜品位*100
TB：在反应塔中为反应塔气体温度K
在沉淀池中为粗铜温度K
2、渣显热＝渣量*渣单位热量
渣单位热量＝0.3*TS-131.9
TS：在反应塔中为反应塔气体温度K
在沉淀池中为粗铜温度K
3、烟气显热＝(SO2量*SO2单位热量+O2量*O2单位热量+NO2量*NO2单位热量+CO2量-3
*CO2单位热量+H2O量*H2O单位热量)*10
其中：单位热量计算通式：
2 -3 -1 3 3
单位热量＝(a*TG+b*TG*10 +C*TG *10+d)*1/22.4*4.187Kj/Nm
系数表：
成分 a b C d
SO2 10.38 1.27 1.42 -3683
O2 7.16 0.5 0.40 -2312
N2 6.66 0.51 0 -2030
CO2 10.55 1.08 2.04 -3926
H2O 7.17 1.28 -0.08 -2223
TG：反应塔热平衡时为反应塔烟气温度+217.16K
沉淀池热平衡时为沉淀池烟气温度+217.16K
上升烟道热平衡时为上升烟道烟气温度+217.16K
4、辐射热损＝上升烟道烟气显热*辐射损失率
5、反应塔烟尘显热＝(产出B烟尘量+产出C烟尘量+烟气中其他量)*烟尘单位热量烟尘单位热量＝0.3*TG-131.9
TG：反应塔热平衡时为反应塔烟气温度+217.16K
沉淀池热平衡时为沉淀池烟气温度+217.16K
上升烟道热平衡时为上升烟道烟气温度+217.16K
6、水分蒸发热＝(冰铜量*冰铜水分率+石灰*石灰水分率+硅酸矿*硅酸矿水分率)*(水分蒸发热+饱和水焓-装入水分焓)
7、反应塔散热＝给定值
8、沉淀池散热＝给定值
9、上升烟道散热＝给定值
10、出热合计＝上述1～9项之和
C、烟气成分计算
1、反应塔烟气成分量
3 3
·SO2量＝烟气中S量*(22.4/32.06)*10 Nm/h
·O2量＝反应塔工艺风O2量*(1-反应塔工艺风氧效率)+反应塔燃烧风O2量*(1-反应塔燃烧风氧效率)
反应塔工艺风O2量＝必要矿石氧气量
反应塔燃烧风O2量＝反应塔燃烧风量*反应塔燃烧风氧浓度
·N2量＝反应塔工艺风N2量+反应塔燃烧风N2量+反应塔天然气量*天然气含N2率反应塔工艺风N2量＝必要矿石氧气量/工艺风富氧浓度*(1-工艺风富氧浓度)反应塔燃烧风N2量＝反应塔燃烧风量湿基*(1-反应塔燃烧风氧浓度)
·CO2量＝天然气C品位*反应塔天然气量*22.4/12+[石灰CaO量*CaCO3比率/56.08+
3
石灰MgO量*MgCO3比率/40.32]*22.4*10+反应塔天然气量*天然气CO2比率·H2O量＝反应塔工艺风H2O量+反应塔燃烧风H2O量+炉料H2O量+天然气H2O量反应塔工艺风H2O量＝反应塔工艺空气量干基*空气水分率
反应塔燃烧风H2O量＝反应塔燃烧空气量*反应塔空气水分率
炉料H2O量＝[冰铜量*冰铜水分率+石灰量*石灰水分率+硅酸矿量*硅酸矿水分
3
率]*22.4/18*10
天然气H2O量＝天然气量*天然气H品位*22.4/2
2、沉淀池烟气成分量
·沉淀池SO2量＝反应塔SO2量
·沉淀池O2量＝反应塔O2量+沉淀池燃烧风O2量*(1-沉淀池氧效率)沉淀池燃烧风O2量＝沉淀池燃烧风量*沉淀池燃烧风氧浓度
·沉淀池N2量＝反应塔N2量+沉淀池燃烧风N2量+沉淀池天然气量*天然气含N2率沉淀池燃烧风N2量＝沉淀池燃烧风量湿基*(1-沉淀池燃烧风氧浓度)
·沉淀池CO2量＝反应塔CO2量+天然气C品位*沉淀池天然气量*22.4/12+沉淀池天然气量*天然气CO2比率
·沉淀池H2O量＝反应塔H2O量+沉淀池燃烧风H2O量+沉淀池天然气H2O量沉淀池燃烧风H2O量＝沉淀池燃烧空气量*沉淀池空气水分率沉淀池天然气H2O量＝沉淀池天然气量*天然气H品位*22.4/2
3、上升烟道烟气成分量
·上升烟道SO2量＝沉淀池SO2量
·上升烟道O2量＝沉淀池O2量+上升烟道燃烧风O2量*(1-上升烟道氧效率)上升烟道燃烧风O2量＝上升烟道燃烧风量*上升烟道燃烧风氧浓度
·上升烟道N2量＝沉淀池N2量+上升烟道燃烧风N2量+上升烟道天然气量*天然气含N2率
上升烟道燃烧风N2量＝上升烟道燃烧风量湿基*(1-上升烟道燃烧风氧浓度)·上升烟道CO2量＝沉淀池CO2量+天然气C品位*上升烟道天然气量*22.4/12+上升烟道天然气量*天然气CO2比率
·上升烟道H2O量＝沉淀池H2O量+上升烟道燃烧风H2O量+上升烟道天然气H2O量上升烟道燃烧风H2O量＝上升烟道燃烧空气量*上升烟道空气水分率沉淀池天然气H2O量＝沉淀池天然气量*天然气H品位*22.4/2。
2.根据权利要求1所述的一种基于数学模型指导的粗铜吹炼工艺，其特征在于：反应塔的工艺风是工艺空气和工业制备氧混合的常温富氧气体，富化率达80％。
3.根据权利要求1所述的一种基于数学模型指导的粗铜吹炼工艺，其特征在于：其配料过程是在吹炼炉炉顶部完成，吹炼炉炉顶分别安装了冰铜、烟尘、石灰和石英沙配料仓，各物料均配有独立的计量装置，经计量给料后的炉料进入两台埋刮板运输机，再由埋刮板运输机送入闪速吹炼炉的喷嘴。