锅炉正反平衡计算热效率

杨铮...

已于 2024-05-14 09:57:32 修改

阅读量2.1k

点赞数 1

分类专栏：算法文章标签： c++ 算法人工智能

于 2022-11-04 20:08:43 首次发布

本文链接：https://blog.csdn.net/weixin_37928884/article/details/127695046

版权

算法专栏收录该内容

37 篇文章 1 订阅

订阅专栏

在这里插入图片描述

热平衡方程：Qr=Q1+Q2+Q3+Q4+Q5+Q6

Qr—送入锅炉的热量

Q1—有效利用热

Q2—排烟热损失

Q3—气体不完全燃烧损失

Q4—固体不完全燃烧损失

Q5—散热损失

Q6—其他热损失

正平衡计算

热效率 ={ 过热蒸汽量*(过热蒸汽焓-给水焓)+排污量*(饱和水焓-给水焓)+再热蒸汽量*(再热蒸汽出口焓-再热蒸汽进口焓)}/(燃料消耗量*燃料热值)

一般锅炉排污量很少可忽略不计，并且没有再热器设备。
上述公式简化：

热效率 ={ 过热蒸汽量*(过热蒸汽焓-给水焓)}/(燃料消耗量*燃料热值)

这里会涉及到温压焓值表，大家可以从网上找找，应该有，没有的话去下载【其实代码中给出了220吨锅炉的温压焓值表，如果你研究的锅炉属于超高压，不懂如何下手分析，我很乐意解答】：https://download.csdn.net/download/weixin_37928884/85505889

过热蒸汽焓值二值插值计算程序如下：

float ahanValue[31][12]={ //      
{   0.0,	   0.1,	   0.5,	   1.0,	   3.0,	   5.0,	   7.1,	 10.1,	  14.1,	  20.1,	  25.1,	  30.0},
{  42.0,	  42.1,	  42.5,	  43.0,	  44.9,	  46.9,	  48.8,	 51.7,	  55.6,	  61.3,	  66.1,	  70.8},
{  83.9,	  84.0,	  84.3,	  84.8,	  86.7,	  88.6,	  90.4,	 93.2,	  97.0,	 102.5,	 107.1,	 111.7},
{ 167.4,	 167.5,	 167.9,	 168.3,	 170.1,	 171.9,	 173.6,	 176.3,	 179.8,	 185.1,	 189.4,	 193.8},
{2611.3,	 251.2,	 251.2,	 251.9,	 253.6,	 255.3,	 256.9,	 259.4,	 262.8,	 267.8,	 272.0,	 276.1},
{2649.3,	 335.0,	 335.3,	 335.7,	 337.3,	 338.8,	 340.4,	 342.8,	 346.0,	 350.8,	 354.8,	 358.7},
{2687.3,	2676.5,	 419.4,	 419.7,	 421.2,	 422.7,	 424.2,	 426.5,	 429.5,	 434.0,	 437.8,	 441.6},
{2725.4,	2716.8,	 503.9,	 504.3,	 505.7,	 507.1,	 508.5,	 510.6,	 513.5,	 517.7,	 521.3,	 524.9},
{2763.6,	2756.6,	 589.2,	 589.5,	 590.8,	 592.1,	 593.4,	 595.4,	 598.0,	 602.0,	 605.4,	 603.1},
{2802.0,	2796.2,	2767.3,	 675.7,	 676.9,	 678.0,	 679.2,	 681.0,	 683.4,	 687.1,	 690.2,	 693.3},
{2840.6,	2835.7,	2812.1,	2777.3,	 764.1,	 765.2,	 766.2,	 767.8,	 769.9,	 773.1,	 775.9,	 778.7},
{2879.3,	2875.2,	2855.5,	2827.5,	 853.0,	 853.8,	 854.6,	 855.9,	 857.7,	 860.4,	 862.8,	 856.2},
{2918.3,	2914.7,	2898.0,	2874.9,	 943.9,	 944.4,	 945.0,	 946.0,	 947.2,	 949.3,	 951.2,	 953.1},
{2957.4,	2954.3,	2939.9,	2920.5,	2823.0,	1037.8,	1038.0,	1038.4,	1039.1,	1040.3,	1041.5,	1024.8},
{2996.8,	2994.1,	2981.5,	2964.8,	2885.5,	1135.0,	1134.7,	1134.3,	1134.1,	1134.0,	1134.3,	1134.8},
{3036.5,	3034.0,	3022.9,	3008.3,	2941.8,	2857.0,	1236.7,	1235.2,	1233.5,	1231.6,	1230.5,	1229.9},
{3076.3,	3074.1,	3064.2,	3051.3,	2994.2,	2925.4,	2839.2,	1343.7,	1339.5,	1334.6,	1331.5,	1329.0},
{3177.0,	3175.3,	3167.6,	3157.7,	3115.7,	3069.2,	3017.0,	2924.2,	2753.5,	1648.4,	1626.4,	1611.3},
{3279.4,	3278.0,	3217.8,	3264.0,	3231.6,	3196.9,	3159.7,	3098.5,	3004.0,	2820.1,	2583.2,	2159.1},
{3321.0,	3319.7,	3313.8,	3306.6,	3276.9,	3245.4,	3211.0,	3156.0,	3072.7,	2917.0,	2730.8,	2424.7},
{3362.5,	3361.4,	3355.9,	3349.3,	3321.9,	3293.2,	3262.3,	3213.5,	3141.4,	3013.9,	2878.3,	2690.3},
{3383.3,	3382.2,	3377.1,	3370.7,	3344.4,	3316.8,	3288.0,	3242.2,	3175.8,	3062.4,	2952.1,	2823.1},
{3404.4,	3403.3,	3398.3,	3392.1,	3366.8,	3340.4,	3312.4,	3268.6,	3205.2,	3098.0,	2994.7,	2875.3},
{3446.7,	3445.6,	3440.9,	3435.1,	3411.6,	3387.2,	3361.3,	3321.3,	3264.1,	3169.1,	3079.8,	2979.6},
{3488.9,	3487.9,	3483.7,	3478.3,	3456.4,	3433.8,	3410.2,	3374.1,	3323.0,	3240.2,	3165.0,	3083.9},
{3531.8,	3530.9,	3526.9,	3521.9,	3501.3,	3480.1,	3458.6,	3425.1,	3378.4,	3303.7,	3237.0,	3166.1},
{3574.7,	3573.9,	3570.1,	3565.4,	3546.2,	3526.4,	3506.4,	3475.4,	3432.5,	3364.6,	3304.7,	3241.7},
{3593.2,	3595.4,	3591.7,	3587.2,	3568.6,	3549.6,	3530.2,	3500.4,	3459.2,	3394.3,	3337.3,	3277.7},
{3618.0,	3617.2,	3613.6,	3609.2,	3591.2,	3572.8,	3554.1,	3525.4,	3485.8,	3423.6,	3369.2,	3312.6},
{3661.6,	3660.9,	3657.5,	3653.3,	3636.3,	3619.1,	3601.6,	3574.9,	3538.2,	3480.9,	3431.2,	3379.8},
{3705.2,	3704.5,	3701.4,	3697.4,	3681.5,	3665.4,	3649.0,	3624.0,	3589.8,	3536.9,	3491.2,	3444.2}}; 

float axlabel[31]={0,10,20,40,60,80,100,120,140,160,180,200,220,240,260,280,300,350,400,420,440,450,460,480,500,520,540,550,560,580,600};
float aylabel[12]={0.01,0.1,0.5,1,3,5,7,10,14,20,25,30};

float aXX = 0;
float aYY = 0;
float aZZ = 0;
int ainum = 0;
int ajnum = 0;
//X：温度
//Y：压力
float asecondctl(float X,float Y)                        
{ 
	float Z1;
	float Z2;
	int i,j;
	
	//找横坐标
	for (i = 0;i < 31;i++)
	{
		if((axlabel[i] < X) && (axlabel[i+1] >= X))
		{
			//aXX = axlabel[i];
			//ainum = i;
			aXX = axlabel[i+1];
			ainum = i+1;
			break;
		}
	}
	//找纵坐标
	for (j = 0;j < 12;j++)
	{
		if((aylabel[j] < Y) && (aylabel[j+1] >= Y))
		{
			//aYY = aylabel[j];
			//ajnum = j;
			//0-HOLLiASMACSV6.5.2功能块说明—参考76～78 下角标从1开始
			aYY = aylabel[j+1];
			ajnum = j+1;
			break;
		}
	}

	Z1 = ahanValue[ainum-1][ajnum-1] + (ahanValue[ainum-1][ajnum]-ahanValue[ainum-1][ajnum-1])*(aYY-aylabel[ajnum-1])/(aylabel[ajnum]-aylabel[ajnum-1]);
	Z2 = ahanValue[ainum][ajnum-1] + (ahanValue[ainum][ajnum]-ahanValue[ainum][ajnum-1])*(aYY-aylabel[ajnum-1])/(aylabel[ajnum]-aylabel[ajnum-1]);
	
	aZZ = Z1 + (Z2-Z1)*(X-axlabel[ainum-1])/(axlabel[ainum]-axlabel[ainum-1]);
	return aZZ;

}

给水焓值二值插值计算程序如下：

float bhanValue[31][12]={ //      
{   0.0,	   0.1,	   0.5,	   1.0,	   3.0,	   5.0,	   7.1,	 10.1,	  14.1,	  20.1,	  25.1,	  30.0},
{  42.0,	  42.1,	  42.5,	  43.0,	  44.9,	  46.9,	  48.8,	 51.7,	  55.6,	  61.3,	  66.1,	  70.8},
{  83.9,	  84.0,	  84.3,	  84.8,	  86.7,	  88.6,	  90.4,	 93.2,	  97.0,	 102.5,	 107.1,	 111.7},
{ 167.4,	 167.5,	 167.9,	 168.3,	 170.1,	 171.9,	 173.6,	 176.3,	 179.8,	 185.1,	 189.4,	 193.8},
{2611.3,	 251.2,	 251.2,	 251.9,	 253.6,	 255.3,	 256.9,	 259.4,	 262.8,	 267.8,	 272.0,	 276.1},
{2649.3,	 335.0,	 335.3,	 335.7,	 337.3,	 338.8,	 340.4,	 342.8,	 346.0,	 350.8,	 354.8,	 358.7},
{2687.3,	2676.5,	 419.4,	 419.7,	 421.2,	 422.7,	 424.2,	 426.5,	 429.5,	 434.0,	 437.8,	 441.6},
{2725.4,	2716.8,	 503.9,	 504.3,	 505.7,	 507.1,	 508.5,	 510.6,	 513.5,	 517.7,	 521.3,	 524.9},
{2763.6,	2756.6,	 589.2,	 589.5,	 590.8,	 592.1,	 593.4,	 595.4,	 598.0,	 602.0,	 605.4,	 603.1},
{2802.0,	2796.2,	2767.3,	 675.7,	 676.9,	 678.0,	 679.2,	 681.0,	 683.4,	 687.1,	 690.2,	 693.3},
{2840.6,	2835.7,	2812.1,	2777.3,	 764.1,	 765.2,	 766.2,	 767.8,	 769.9,	 773.1,	 775.9,	 778.7},
{2879.3,	2875.2,	2855.5,	2827.5,	 853.0,	 853.8,	 854.6,	 855.9,	 857.7,	 860.4,	 862.8,	 856.2},
{2918.3,	2914.7,	2898.0,	2874.9,	 943.9,	 944.4,	 945.0,	 946.0,	 947.2,	 949.3,	 951.2,	 953.1},
{2957.4,	2954.3,	2939.9,	2920.5,	2823.0,	1037.8,	1038.0,	1038.4,	1039.1,	1040.3,	1041.5,	1024.8},
{2996.8,	2994.1,	2981.5,	2964.8,	2885.5,	1135.0,	1134.7,	1134.3,	1134.1,	1134.0,	1134.3,	1134.8},
{3036.5,	3034.0,	3022.9,	3008.3,	2941.8,	2857.0,	1236.7,	1235.2,	1233.5,	1231.6,	1230.5,	1229.9},
{3076.3,	3074.1,	3064.2,	3051.3,	2994.2,	2925.4,	2839.2,	1343.7,	1339.5,	1334.6,	1331.5,	1329.0},
{3177.0,	3175.3,	3167.6,	3157.7,	3115.7,	3069.2,	3017.0,	2924.2,	2753.5,	1648.4,	1626.4,	1611.3},
{3279.4,	3278.0,	3217.8,	3264.0,	3231.6,	3196.9,	3159.7,	3098.5,	3004.0,	2820.1,	2583.2,	2159.1},
{3321.0,	3319.7,	3313.8,	3306.6,	3276.9,	3245.4,	3211.0,	3156.0,	3072.7,	2917.0,	2730.8,	2424.7},
{3362.5,	3361.4,	3355.9,	3349.3,	3321.9,	3293.2,	3262.3,	3213.5,	3141.4,	3013.9,	2878.3,	2690.3},
{3383.3,	3382.2,	3377.1,	3370.7,	3344.4,	3316.8,	3288.0,	3242.2,	3175.8,	3062.4,	2952.1,	2823.1},
{3404.4,	3403.3,	3398.3,	3392.1,	3366.8,	3340.4,	3312.4,	3268.6,	3205.2,	3098.0,	2994.7,	2875.3},
{3446.7,	3445.6,	3440.9,	3435.1,	3411.6,	3387.2,	3361.3,	3321.3,	3264.1,	3169.1,	3079.8,	2979.6},
{3488.9,	3487.9,	3483.7,	3478.3,	3456.4,	3433.8,	3410.2,	3374.1,	3323.0,	3240.2,	3165.0,	3083.9},
{3531.8,	3530.9,	3526.9,	3521.9,	3501.3,	3480.1,	3458.6,	3425.1,	3378.4,	3303.7,	3237.0,	3166.1},
{3574.7,	3573.9,	3570.1,	3565.4,	3546.2,	3526.4,	3506.4,	3475.4,	3432.5,	3364.6,	3304.7,	3241.7},
{3593.2,	3595.4,	3591.7,	3587.2,	3568.6,	3549.6,	3530.2,	3500.4,	3459.2,	3394.3,	3337.3,	3277.7},
{3618.0,	3617.2,	3613.6,	3609.2,	3591.2,	3572.8,	3554.1,	3525.4,	3485.8,	3423.6,	3369.2,	3312.6},
{3661.6,	3660.9,	3657.5,	3653.3,	3636.3,	3619.1,	3601.6,	3574.9,	3538.2,	3480.9,	3431.2,	3379.8},
{3705.2,	3704.5,	3701.4,	3697.4,	3681.5,	3665.4,	3649.0,	3624.0,	3589.8,	3536.9,	3491.2,	3444.2}}; 

float bxlabel[31]={0,10,20,40,60,80,100,120,140,160,180,200,220,240,260,280,300,350,400,420,440,450,460,480,500,520,540,550,560,580,600};
float bylabel[12]={0.01,0.1,0.5,1,3,5,7,10,14,20,25,30};

float bXX = 0;
float bYY = 0;
float bZZ = 0;
int binum = 0;
int bjnum = 0;
float bsecondctl(float X,float Y)                        
{ 
	float Z1;
	float Z2;
	int i,j;
	
	//找横坐标
	for (i = 0;i < 31;i++)
	{
		if((bxlabel[i] < X) && (bxlabel[i+1] >= X))
		{
			//bXX = bxlabel[i];
			//inum = i;
			//HOLLiASMACSV6.5.2功能块说明—参考76～78 下角标从1开始
			bXX = bxlabel[i+1];
			binum = i+1;
			break;
		}
	}
	//找纵坐标
	for (j = 0;j < 12;j++)
	{
		if((bylabel[j] < Y) && (bylabel[j+1] >= Y))
		{
			//bYY = bylabel[j];
			//bjnum = j;
			//0-HOLLiASMACSV6.5.2功能块说明—参考76～78 下角标从1开始
			bYY = bylabel[j+1];
			bjnum = j+1;
			break;
		}
	}

	Z1 = bhanValue[binum-1][bjnum-1] + (bhanValue[binum-1][bjnum]-bhanValue[binum-1][bjnum-1])*(bYY-bylabel[bjnum-1])/(bylabel[bjnum]-bylabel[bjnum-1]);
	Z2 = bhanValue[binum][bjnum-1] + (bhanValue[binum][bjnum]-bhanValue[binum][bjnum-1])*(bYY-bylabel[bjnum-1])/(bylabel[bjnum]-bylabel[bjnum-1]);
	
	bZZ = Z1 + (Z2-Z1)*(bXX-bxlabel[binum-1])/(bxlabel[binum]-bxlabel[binum-1]);
	return bZZ;

}

剩下的计算就简单了，这里留给你们。正平衡就是上面这点东西。

反平衡计算
反平衡计算稍微麻烦点，计算气体不完全燃烧损失中涉及到烟气量问题，这里的烟气量不是仪表测出的量，应该是通过煤气消耗量*单位煤气燃烧产生的烟气量计算出来的。这里的细节太多想了解更多，请移步
请参考：https://download.csdn.net/download/weixin_37928884/85505889

这里给和利时DCS点赞！《HOLLiASMACSV6.5.2功能块说明》帮助文档真的很人性化。