碳平衡法排放計(jì)算流程.doc

1. 發(fā)動(dòng)機(jī)試驗(yàn)狀況：實(shí)驗(yàn)室大氣系數(shù)fa試驗(yàn)時(shí)應(yīng)在參數(shù)0.96fa1.06之間進(jìn)行。（1）對(duì)于壓燃式帶進(jìn)氣中冷的渦輪增壓式發(fā)動(dòng)機(jī)，fa應(yīng)按下式進(jìn)行計(jì)算：fa=99ps0.7Ta2981.5Ta發(fā)動(dòng)機(jī)進(jìn)氣絕對(duì)溫度，K；ps發(fā)動(dòng)機(jī)進(jìn)氣干氣壓，kPa；psph-ehpa發(fā)動(dòng)機(jī)進(jìn)氣濕氣壓，kPa；ea發(fā)動(dòng)機(jī)進(jìn)氣蒸氣壓，kPa，進(jìn)氣蒸氣壓的計(jì)算見下面絕對(duì)濕度Ha的計(jì)算過程中。0.93fa1.07例，當(dāng)測(cè)得進(jìn)氣濕壓力pin=100.6kPa，tin31.4，相對(duì)濕度RHin54%時(shí)，可計(jì)算fa數(shù)值得，fa=99100.6-2.481610.7273.15+31.42981.51.022. 濕基排氣質(zhì)量流量qmew的確定：所有與體積有關(guān)的量綱均按標(biāo)準(zhǔn)狀況（0，101.325kPa）確定，CO2，O2，H2O，N2組分所用濃度單位為體積百分?jǐn)?shù)，其他用ppm。（1） 直接測(cè)量法本法需要對(duì)排氣流量直接進(jìn)行測(cè)量。排氣流量可用諸如下列設(shè)備直接測(cè)量：壓差計(jì)，如流量噴嘴（見ISO 5167）；超聲波流量計(jì)；渦街流量計(jì)。排氣流量的測(cè)量允許偏差為讀數(shù)的2.5%或發(fā)動(dòng)機(jī)最大排氣流量的1.5%，以較大值為準(zhǔn)。（2） 空氣和燃料測(cè)量法本法需要對(duì)空氣流量和燃料流量進(jìn)行測(cè)量。排氣流量計(jì)算公式如下：qmew=qmaw+qmfqmew濕基排氣質(zhì)量流量，kg/h；qmaw濕基進(jìn)氣空氣質(zhì)量流量，kg/h；qmf 燃料質(zhì)量流量，kg/h；空氣消耗量的允許偏差為讀數(shù)的2%或發(fā)動(dòng)機(jī)最大排氣流量的1%，以較大值為準(zhǔn)。燃料流量的允許偏差為發(fā)動(dòng)機(jī)最大燃料消耗量的2%。（3） 碳平衡：迭代計(jì)算法干排氣質(zhì)量流量qmed的計(jì)算需要有干排氣密度ed和干濕基修正系數(shù)kwr值，其本身與qmad值有關(guān)，從而就與qmed的計(jì)算結(jié)果有關(guān)。因此，在下述方法中將應(yīng)用迭代（或多步）計(jì)算法。利用ed和kwr的初始值（例如1.34kg/m3和1）計(jì)算qmed值，由此計(jì)算qmad值，再由這些值計(jì)算ed和kwr值。根據(jù)這些幾乎已經(jīng)很精確的ed和kwr值，再利用相同公式進(jìn)行下面迭代計(jì)算時(shí)，所有數(shù)據(jù)值就已足夠精確，因此在大多數(shù)情況下就不必再進(jìn)行第三次迭代計(jì)算。 步驟一：利用ed和kwr的初始值（例如1.34kg/m3和1）計(jì)算qmed值先根據(jù)碳平衡法用下列公式計(jì)算干排氣質(zhì)量流量qmed：qmedqmfBETed104ArCCCO2d-CCO2 ,a104VmCO2+CCOdVmCO11-prpd+CHCwVmHC+CCwArC1kwrCCO2d 原排氣中干CO2濃度，%； 測(cè)得CCO2a環(huán)境中干CO2濃度，%； 已知，0.04%CCOd 原排氣中干CO濃度，%； 測(cè)得CHC 原排氣中濕HC濃度，ppm； 測(cè)得CCw碳煙濃度，mg/m3，計(jì)算時(shí)，可以用波許煙度值SN計(jì)算得到碳煙濃度，在此計(jì)算時(shí)忽略了碳煙的濃度，即令CCw0。 步驟二：利用步驟一得到的干排氣質(zhì)量流量qmed計(jì)算干進(jìn)氣空氣質(zhì)量流量qmadqmadqmed-qmf1-ALF0.08936 步驟三：利用步驟二得到的干進(jìn)氣空氣質(zhì)量流量qmad和測(cè)得的燃料質(zhì)量流量qmf計(jì)算ed和kwredqmad+qmf1-ALF0.08936qmad1.293+ffdqmfkwr1-1.2442Ha+111.187ALFqmfqmad-773.4prpb773.4+1.2442Ha+qmfqmadffd1000其中，ffd-0.055593ALF+0.008002DEL+0.0070046EPS 步驟四：利用步驟三計(jì)算得到的ed和kwr值再迭代入步驟一，步驟二進(jìn)行計(jì)算得到qmad，然后再利用迭代后計(jì)算得到的qmad代入下式中進(jìn)行濕基排氣質(zhì)量流量qmew的計(jì)算qmew=qmad1+Ha1000+qmf將步驟一，步驟二，步驟四合并成一個(gè)公式，并進(jìn)行一些簡(jiǎn)化忽略未燃燒的碳煙和假設(shè)冷卻器溫度固定為4，即1/(1-pr/pb)=1.008，可以利用步驟三計(jì)算得到的ed和kwr值直接代入下面更容易使用的濕排氣質(zhì)量流量計(jì)算公式中計(jì)算qmew，qmewqmfBETedCCO2d-CCO2 ,a0.540+CCOd+CHCwkwr10000.00471+ALF0.08936-11+Ha1000+1（4） 碳平衡：一步計(jì)算法本法需要對(duì)燃料消耗、燃料組分和排氣組分濃度進(jìn)行測(cè)量。注：在本方法的推導(dǎo)中，假設(shè)的前提是：當(dāng)pb=1013mBar，pr=7.5mBar（冷卻器溫度為4），或1/1-prpb=1.008、kwr=0.93和CCw0qmew=qmfBETBET1.41.4BETfc+ALF0.08936-111.293+ffdfcfc+ALF0.08936 -11+Ha1000+1qmew濕基排氣質(zhì)量流量，kg/h； 計(jì)算量qmf 燃料質(zhì)量流量，kg/h； 測(cè)得與燃料組分有關(guān)，已知量BET燃料含C量，%質(zhì)量；ALF燃料含H量，%質(zhì)量；ffd 干基排氣流量計(jì)算用燃料特定系數(shù)；fc 碳系數(shù)；Ha 進(jìn)氣絕對(duì)濕度，又稱含濕量，即每千克干空氣含水量g，單位為g（水）/kg（干空氣）。 干基排氣流量計(jì)算用燃料特定系數(shù)ffd的計(jì)算：ffd=qved-qvadqmfqmf 燃料質(zhì)量流量，kg/h；qved 干排氣體積流量，m3/h；qvad 干進(jìn)氣空氣體積流量，m3/h；ffd 干基排氣流量計(jì)算用燃料特定系數(shù)；ffw 濕基排氣流量計(jì)算用燃料特定系數(shù)，m3，每kg燃料從燃燒空氣到濕排氣的體積變化。對(duì)于燃料組分為CHSNO的燃燒，一般燃燒后，由于會(huì)產(chǎn)生更多不同的氣體，所以體積會(huì)增大；對(duì)于像H2類的燃料，由于燃燒后，會(huì)產(chǎn)生水，體積有可能會(huì)減小。qved的計(jì)算應(yīng)從燃燒（qmfffw）所產(chǎn)生的總體積變化中除去燃燒產(chǎn)生的水蒸氣，產(chǎn)生的水蒸氣的體積計(jì)算如下：ALFqmfVmH2O1002ArHALF燃料含H量，%質(zhì)量；VmH2O水摩爾體積，22.414 L/mol；ArH 氫原子的原子質(zhì)量，g，ArH1.00794 g。則，qvedqvad+qmfffw-ALFqmfVmH2O1002ArH則，ffd=ffwALFVmH2O200ArHffwALF22.4142001.00794ffwALF0.11118又知，ffw0.055594ALF+0.0080021DEL+0.0070046EPS代入上式得，ffd=-0.055593ALF+0.0080021DEL+0.0070046EPS與燃料組分有關(guān)，已知量ALF燃料含H量，%質(zhì)量；DEL燃料含N量，%質(zhì)量；EPS燃料含O量，%質(zhì)量； fc的計(jì)算：fc=CCO2d-CCO2ad0.5441+CCOd18522+CHC17355CCO2d 原排氣中干CO2濃度，%； 測(cè)得CCO2ad環(huán)境中干CO2濃度，%； 已知，后面計(jì)算組分排放質(zhì)量流量時(shí)需要CCO2awCCOd 原排氣中干CO濃度，%； 測(cè)得CHC 原排氣中濕HC濃度，ppm； 測(cè)得對(duì)于燃燒的理想化學(xué)計(jì)算，假設(shè)干進(jìn)氣空氣的組分如下：A） 惰性氣體濃度inert=76.8% mass 79.0% vol惰性氣體中的CO2其質(zhì)量濃度為0.061%，體積濃度為0.04%B） 氧氣濃度OX=23.2% mass 21.0% vol Ha進(jìn)氣絕對(duì)濕度（又稱含濕量）的計(jì)算：一般在測(cè)量過程中，不會(huì)對(duì)進(jìn)氣絕對(duì)溫度進(jìn)行測(cè)量，而往往會(huì)測(cè)量進(jìn)氣口的溫度tin（），壓力pin（kPa），相對(duì)濕度RHin（%），由此三個(gè)已知量來推導(dǎo)進(jìn)氣絕對(duì)濕度。此外，Ha 在GBT8190.1 Pag32關(guān)于NOx的干濕基修正內(nèi)容中標(biāo)為空氣濾清器入口處的大氣濕度，即表明了Ha基本可以等價(jià)為大氣的濕度，并不是經(jīng)增壓冷卻后進(jìn)入穩(wěn)壓箱中的空氣的濕度。推導(dǎo)過程如下：A. 該溫度下的飽和蒸氣壓（又稱飽和水汽壓）es（kPa）飽和蒸氣壓的計(jì)算有多種方法：（）飽和蒸氣壓僅與溫度t有關(guān)，es=0.1es0expat273.15+t-bes00時(shí)的水氣壓，值為6.1078；t 水蒸氣的溫度，。乘以0.1的原因是單位轉(zhuǎn)換的需要，因?yàn)?.1后面計(jì)算得到的結(jié)果單位為hpa（百帕）。水面和冰面飽和水汽壓不同，a、b系數(shù)取值不同：t15水面：a=17.269，b=35.86-40t-15冰水混合：es=0.00280+2tes水-(30+2t)es冰t40冰面：a=21.8746，b=7.66（）根據(jù)GB/T 8190.1 A.15的簡(jiǎn)化公式：es=4.856884+0.2660089t+0.01688919t2-7.47712310-5t3+8.1052510-6t4-3.11522110-8t51013.2760文檔中所列公式為GB/T 8190.1 A.15，標(biāo)準(zhǔn)中約定的計(jì)算單位為kPa，但按該單位得到的飽和蒸氣濕度比實(shí)際大10倍，個(gè)人猜想標(biāo)準(zhǔn)中單位應(yīng)該是hPa（百帕），而不是kPa。（）采用Antoine公式（適用于0150）lgP=A-B/(t+C)P物質(zhì)的蒸氣壓，毫米汞柱，1mm汞柱=133.3Pa；t溫度，；A、B、C系數(shù)，查表可得。溫度t范圍（）ABC601507.966811668.210228.0000-608.107651750.286235.000經(jīng)整理可得，es=133.3100010A-B/(t+C)B. 該溫度下的蒸氣壓（又稱水汽壓）e（kPa）eRHes100%注：如果相對(duì)濕度測(cè)量值為45%則，上式中的RH應(yīng)為45，即e45es100%0.45esC. 該溫度下的進(jìn)氣絕對(duì)濕度Ha（g/kg）Ha=622ep-e注：一般情況下，由于柴油機(jī)進(jìn)氣口的溫度ta（），壓力pa（kPa），相對(duì)濕度RHa（%），與環(huán)境溫度、壓力、相對(duì)濕度非常接近，所以上述三個(gè)量也可以用環(huán)境溫度、壓力、相對(duì)濕度來代替，但在有條件的情況下，測(cè)量點(diǎn)應(yīng)盡量在進(jìn)氣口處進(jìn)行測(cè)量。對(duì)于有渦輪增壓并有中冷的柴油機(jī)的進(jìn)氣絕對(duì)濕度，由于經(jīng)過中冷器冷卻后，空氣一般都會(huì)成為飽和蒸氣，所以此時(shí)的進(jìn)氣絕對(duì)濕度應(yīng)為飽和絕對(duì)濕度，只通過測(cè)量穩(wěn)壓箱溫度、與穩(wěn)壓箱壓力， 不必利用相對(duì)濕度即可計(jì)算。（5） 氧平衡法與碳平衡法相同，氧平衡法需要對(duì)燃料消耗、燃料組分和排氣組分濃度進(jìn)行測(cè)量。與理論排氣質(zhì)量流量相比，氧平衡法的偏差要稍大一些（可達(dá)1%，而碳平衡法小于0.2%）。因此應(yīng)優(yōu)先采用碳平衡法。（6） 空氣流量和空燃比測(cè)量法利用空氣流量和空燃比來計(jì)算排氣質(zhì)量，瞬時(shí)排氣質(zhì)量流量的計(jì)算如下：qmew=qmaw1+1AFstqmew濕基排氣質(zhì)量流量，kg/h；qmaw濕基進(jìn)氣空氣質(zhì)量流量，kg/h； 理想空燃比AFst的計(jì)算：AFst138.0+4-2+12.001+1.00794+15.9994+14.0067+32.065AFst理想空燃比，kg/kg；（H/C）摩爾比；與燃料組分有關(guān)，已知量（C/C）摩爾比；（S/C）摩爾比；（N/C）摩爾比；（O/C）摩爾比；注：燃料組分為CHSNO時(shí)，1。對(duì)于不含C的燃料（如：氫），不能用上式計(jì)算。 過量空氣系數(shù)的計(jì)算：=100-CCO10-42-CHC10-4+41-2CCO10-43.5CCO21+CCO10-43.5CCO2-2-2CCO2+CCO10-44.764+4-2+CCO2+CCO10-4+CHC10-4 過量空氣系數(shù)； 測(cè)得CCO2 干CO2濃度，%； 測(cè)得CCO 干CO濃度，ppm； 測(cè)得CHC HC濃度，ppm； 測(cè)得注：傳感器應(yīng)直接安裝在排氣溫度足夠高的位置，以防產(chǎn)生凝結(jié)水。內(nèi)裝電子元器件的傳感器，其準(zhǔn)確度應(yīng)符合下列要求：2 讀數(shù)的3%25 讀數(shù)的5%5 讀數(shù)的10%3. 確定濕基排氣質(zhì)量流量中各組分的濕基質(zhì)量流量qmgas（1） 干/濕基修正如果測(cè)得的各排氣組分未按濕基測(cè)量，當(dāng)需要將其干基濃度轉(zhuǎn)換成濕基濃度時(shí)；如果測(cè)得（或計(jì)算得）的是干進(jìn)氣空氣的進(jìn)氣質(zhì)量流量，當(dāng)需要進(jìn)行干/基修正時(shí)；Cgasw=kwCgasdCgasw濕基排氣中各組分的濕基濃度，下標(biāo)“gas”是指各種氣體組分； 計(jì)算得Cgasd各組分的干基濃度； 測(cè)得kw干/濕基修正系數(shù)； 計(jì)算得干/濕基修正系數(shù)kw的計(jì)算，有以下幾種情況： 對(duì)于原排氣a) 當(dāng)燃料完全燃燒時(shí)本方法需要測(cè)得燃料質(zhì)量流量qmf，干進(jìn)氣空氣的質(zhì)量流量qmad。kwr1=1-1.2442Ha+111.19ALFqmfqmad773.4+1.2442Ha+qmfqmadffw10001-prpbqmf 燃料質(zhì)量流量，kg/h； 測(cè)得qmad干進(jìn)氣空氣質(zhì)量流量，kg/h； 測(cè)得ffw 濕基排氣流量計(jì)算用燃料特定系數(shù)； 計(jì)算得，僅與燃料組分有關(guān)ALF燃料含H量，%質(zhì)量； 僅與燃料組分有關(guān)Ha 進(jìn)氣絕對(duì)濕度，g（水）/kg（空氣）。 計(jì)算得pr冷卻槽后的水蒸氣分壓； 與分析儀的冷凝溫度設(shè)置有關(guān)prpb水蒸氣的摩爾分量（冷卻器后水蒸氣的體積分量）； 與冷凝溫度設(shè)置有關(guān)因?yàn)閷?duì)于CO、CO2、O2濃度的測(cè)量，CO和CO2的測(cè)量是利用不分光紅外線分析儀測(cè)得，如果有水蒸氣的存在，水蒸氣也會(huì)吸收部分紅外線，從而會(huì)對(duì)濃度的測(cè)量產(chǎn)生正干擾誤差，即測(cè)得的CO和CO2濃度比實(shí)際濃度大。而對(duì)于O2的測(cè)量由于是利用順磁法測(cè)得的，如果有水蒸氣的存在，同樣水蒸氣會(huì)產(chǎn)生類似于O2的性質(zhì)，而產(chǎn)生正干擾誤差。因而，當(dāng)對(duì)CO、CO2、O2濃度進(jìn)行測(cè)量時(shí)，需要將樣氣中的水蒸氣去除，具體方法即將樣氣進(jìn)行冷凝（一般冷凝單元的冷凝溫度在4，如CAI的SCU單元），當(dāng)冷凝溫度固定在4時(shí)，可求得11-prpb1.008。故，kwr1=1-1.2442Ha+111.19ALFqmfqmad773.4+1.2442Ha+qmfqmadffw10001.008ffw0.055594ALF+0.0080021DEL+0.0070046EPSb) 當(dāng)燃料不完全燃燒時(shí)（當(dāng)采用碳平衡法時(shí)常用）這種計(jì)算方法尤其適用于濃燃料空氣混合氣（當(dāng)存在相當(dāng)數(shù)量未燃燒的組分（CO，H2）時(shí)）的情況，而且也適用于不直接測(cè)量空氣流量的排放測(cè)量。應(yīng)使用下列公式求出（推導(dǎo)見附錄A.24-A.41）：kwr2=11+0.005CCO2d+CCOd-0.01CH2d+kw2-prpbkw2進(jìn)氣空氣中水分，由GBT 8190.1 A.39 計(jì)算得到，kw21.608Ha1000+1.608HaHa空氣濾清器入口處的大氣濕度，g水/kg干空氣。（Ha 在GBT8190.1 Pag32關(guān)于NOx的干濕基修正內(nèi)容中標(biāo)為空氣濾清器入口處的大氣濕度，即表明了Ha基本可以等價(jià)為大氣的濕度，并不是經(jīng)增壓冷卻后進(jìn)入穩(wěn)壓箱中的空氣的濕度。）根據(jù)SAE J 1088用下列公式由水氣平衡原理推算出氫的濃度：CH2d=0.5CCOd104CCO2d+CCOd104CCOd104+3CCO2dCCO2d 干CO2濃度，%； 測(cè)得CCOd 干CO濃度，ppm； 測(cè)得CH2d 干H2濃度，%； 計(jì)算得pr冷卻槽后的水蒸氣分壓； 與分析儀的冷凝溫度設(shè)置有關(guān)prpb水蒸氣的摩爾分量（冷卻器后水蒸氣的體積分量）； 與冷凝溫度設(shè)置有關(guān)當(dāng)冷凝溫度固定在4時(shí)，可求得11-prpb1.008。 對(duì)于稀釋空氣kw2=1.608Hd1-1D+Ha1D1000+1.608Hd1-1D+Ha1DD稀釋系數(shù)；Ha、Hd進(jìn)氣空氣和稀釋空氣的絕對(duì)濕度，單位為克每千克（g/kg）。 對(duì)于進(jìn)氣空氣（當(dāng)需要通過干進(jìn)氣空氣轉(zhuǎn)換為濕基進(jìn)氣空氣時(shí)常用）kwa=1-kw2其中，kw2=1.608Ha1000+1.608HaHa 進(jìn)氣絕對(duì)濕度，g（水）/kg（空氣）。注：Ha可用通用公式由相對(duì)濕度、露點(diǎn)、蒸氣壓或干/濕球的測(cè)量值求得，具體計(jì)算可看上述部分中關(guān)于Ha的詳細(xì)計(jì)算過程。 NOx的濕度和溫度修正（不同于CO、CO2、O2，對(duì)于NOx需單獨(dú)進(jìn)行干/濕基修正）因NOx排放與大氣狀況有關(guān)，應(yīng)采用下列公式中給出的系數(shù)對(duì)NOx濃度進(jìn)行大氣溫度和濕度修正，當(dāng)濕度在0g/kg和25g/kg之間時(shí)這些系數(shù)才有效。a) 對(duì)于增壓中冷壓燃式發(fā)動(dòng)機(jī)，可用下列公式：khd=11-0.012Ha-10.71-0.00275Ta-298+0.00285(TSC-TSCRef)Ha空氣濾清器入口處的大氣濕度，g水/kg干空氣；Ta空氣濾清器入口處的大氣溫度，K；TSC中冷空氣溫度，K；TSCRef制造商規(guī)定的中冷空氣基準(zhǔn)溫度，K。（2） 濕基排氣中各組分濕基質(zhì)量流量qmgas的計(jì)算應(yīng)根據(jù)污染物的原濃度、表7中的u值和按照上面求得的排氣質(zhì)量流量計(jì)算每工況的排放質(zhì)量流量。如果按干基測(cè)量濃度，則應(yīng)先對(duì)濃度值進(jìn)行干-濕基修正然后再作進(jìn)一點(diǎn)計(jì)算。如果試驗(yàn)所用燃料不