TWI802434B - 粉塵濃度訊號處理裝置及其訊號處理方法 - Google Patents
粉塵濃度訊號處理裝置及其訊號處理方法 Download PDFInfo
- Publication number
- TWI802434B TWI802434B TW111121251A TW111121251A TWI802434B TW I802434 B TWI802434 B TW I802434B TW 111121251 A TW111121251 A TW 111121251A TW 111121251 A TW111121251 A TW 111121251A TW I802434 B TWI802434 B TW I802434B
- Authority
- TW
- Taiwan
- Prior art keywords
- dust concentration
- value
- signal processing
- dust
- unit
- Prior art date
Links
- 239000000428 dust Substances 0.000 title claims abstract description 217
- 238000012545 processing Methods 0.000 title claims description 54
- 238000000034 method Methods 0.000 title claims description 13
- 238000005070 sampling Methods 0.000 claims abstract description 34
- 238000001514 detection method Methods 0.000 claims abstract description 29
- 238000003672 processing method Methods 0.000 claims abstract description 26
- 238000004364 calculation method Methods 0.000 claims description 57
- 238000006243 chemical reaction Methods 0.000 claims description 38
- 238000012937 correction Methods 0.000 claims description 27
- 238000001914 filtration Methods 0.000 claims description 24
- 230000007613 environmental effect Effects 0.000 claims description 17
- 238000012360 testing method Methods 0.000 claims description 15
- 239000000523 sample Substances 0.000 description 17
- 238000010586 diagram Methods 0.000 description 15
- 238000005259 measurement Methods 0.000 description 11
- 239000002245 particle Substances 0.000 description 10
- 239000012530 fluid Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 238000000149 argon plasma sintering Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000007664 blowing Methods 0.000 description 1
- 239000013590 bulk material Substances 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000010881 fly ash Substances 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 230000031700 light absorption Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 231100000817 safety factor Toxicity 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000000779 smoke Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 238000004056 waste incineration Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/06—Investigating concentration of particle suspensions
- G01N15/075—Investigating concentration of particle suspensions by optical means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N2015/0042—Investigating dispersion of solids
- G01N2015/0046—Investigating dispersion of solids in gas, e.g. smoke
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Databases & Information Systems (AREA)
- Dispersion Chemistry (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Evolutionary Biology (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Probability & Statistics with Applications (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Sampling And Sample Adjustment (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measurement Of Radiation (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
一種粉塵濃度訊號處理方法,係接收粉塵濃度計的偵測訊號,且粉塵濃度訊號處理方法包括下列步驟:(a)對偵測訊號進行取樣而計算複數個特徵值,且基於該些特徵值判斷偵測訊號的訊號類型。(b)基於該些特徵值選擇特定運算法,且該些特徵值通過特定運算法計算出粉塵濃度估算值。(c)對粉塵濃度估算值進行校正補償而獲得粉塵濃度精確值。(d)對粉塵濃度精確值進行數值轉換,以提供特定資料格式的粉塵濃度。
Description
本發明係有關一種粉塵濃度訊號處理裝置及其訊號處理方法,尤指一種可降低雜訊影響的粉塵濃度訊號處理裝置及其訊號處理方法。
按,所謂總懸浮粒子濃度係指單位體積空氣中的粉塵含量。在對於廠房或其它較要求總懸浮粒子濃度排放規範,特別是排放到外部環境的情況下,有環保法規等要求。為了量測其數據是否符合標準,必須通過一定的標定程序及相關的量測設備為之,才能建立廠房設備在排煙道上安裝的量測不確定度與計量準確度。
傳統的光學式量測安裝於排煙道,是採用光學吸收(Light Absorption)或散射(Light Scattering)原理,偵測排煙到一端的光發射源與另一端的接收源的訊號,作為校準。然而由於在排煙道中往往會有雜物(例如但不限於樹葉、紙屑、蒸氣、飛灰等)入侵,造成光學路徑上雜訊訊號高於真實粉塵濃度訊號,必須維護設備的頻率必須以”小時”,重複校準,才能得到可靠的基準測量。
此外,排煙道的幾何形狀因子差異(圓形截面、方形截面),以及上游下游各種制動元件(閥件、節流閥、孔口版等),其流場分布與壓力差異,往往會造成排煙道內的器流不平均,容易產生流體渦流(Vortex Shedding)的狀
況。因此,判讀裝置所收到的訊號,參雜太多的即時雜訊,導致難以準確的判讀精確的粉塵濃度。所以,如何設計出一種可降低雜訊影響的粉塵濃度訊號處理裝置,以及即時訊號處理方法,改善排煙道渦流產生的雜訊影響,提供較為精確校正的粉塵濃度,乃為本案創作人所欲行研究的一大課題。
為了解決上述問題,本發明係提供一種可降低雜訊影響的粉塵濃度訊號處理方法,以克服習知技術的問題。因此,本發明粉塵濃度訊號處理方法,係對粉塵濃度計的偵測訊號進行訊號分析而計算出粉塵濃度,且粉塵濃度訊號處理方法包括下列步驟:(a)對偵測訊號進行取樣而計算複數個特徵值,且基於該些特徵值判斷偵測訊號的訊號類型。(b)基於該些特徵值選擇特定運算法,且該些特徵值通過特定運算法計算出粉塵濃度估算值。(c)對粉塵濃度估算值進行校正補償而獲得粉塵濃度精確值。(d)對粉塵濃度精確值進行數值轉換,以提供特定資料格式的粉塵濃度。
為了解決上述問題,本發明係提供一種可降低雜訊影響的粉塵濃度訊號處理裝置,以克服習知技術的問題。因此,本發明塵濃度訊號處理裝置係對粉塵濃度計的偵測訊號進行取樣而計算粉塵濃度,且粉塵濃度訊號處理裝置包括特徵處理模組、運算模組、校正補償模組及數值轉換模組。特徵處理模組對偵測訊號進行取樣而計算複數個特徵值,且基於該些特徵值判斷偵測訊號的訊號類型。運算模組基於該些特徵值選擇特定運算法,且該些特徵值通過特定運算法計算出粉塵濃度估算值。校正補償模組對粉塵濃度估算值進行校正補償而獲
得粉塵濃度精確值。數值轉換模組對粉塵濃度精確值進行數值轉換,以提供特定資料格式的粉塵濃度。
本發明之主要目的及功效在於,提供排煙道總懸浮粒子(TSP,Total Suspended Particulates)濃度標定訊號處理方法,其包含非同步訊號的線上取樣特徵與處理方法,對於流體渦流的變化響應,軟體訊號的取樣頻率、訊號做小波濾波或移動平均等處理,以達成改善排煙道渦流產生的雜訊影響,且提供較為精確的粉塵濃度之功效。
為了能更進一步瞭解本發明為達成預定目的所採取之技術、手段及功效,請參閱以下有關本發明之詳細說明與附圖,相信本發明之目的、特徵與特點,當可由此得一深入且具體之瞭解,然而所附圖式僅提供參考與說明用,並非用來對本發明加以限制者。
100:粉塵監控設備
1:粉塵濃度計
12:探棒
14:感測電路
142:放大電路
142-1:電荷放大器
142-2:可程式放大器
144:濾波電路
146:類比/數位轉換器
2:粉塵濃度訊號處理裝置
22:特徵處理模組
222:取樣單元
224:特徵計算單元
226:類型判斷單元
228:濾波單元
24:運算模組
242:第一運算單元
244:第二運算單元
26:校正補償模組
262:校正單元
264:補償單元
28:數值轉換模組
282:輸出處理單元
284:格式轉換單元
284-1:顯示轉換單元
284-2:類比轉換單元
284-3:閾值比較單元
284-4:RS485轉換單元
200:排煙道
A、B:點位
D:粉塵
Cd:粉塵濃度
Sa:放大訊號
Sf:濾波訊號
Sd:偵測訊號
Vs:取樣值
Es:特徵值
Cc:粉塵濃度估算值
Ccv:粉塵濃度校正值
Ce:粉塵濃度精確值
C1~C6:曲線
Fd:顯示輸出格式
Fa:類比輸出格式
Fr:繼電器輸出格式
Frs:RS485輸出格式
(S100)~(S280):步驟
圖1為本發明粉塵監控設備實施方式的示意圖;圖2為本發明感測電路的電路方塊示意圖;圖3為本發明粉塵濃度訊號處理裝置的電路方塊示意圖;圖4A為本發明粉塵連續排放的特徵示意圖;圖4B為本發明粉塵間歇性排放的特徵示意圖;圖5為本發明濾波單元濾波波形示意圖;圖6為本發明校正測試曲線示意圖;圖7為本發明數值轉換模組的電路方塊圖;圖8為本發明量測範圍曲線示意圖;及
圖9為本發明粉塵濃度訊號處理方法的方法流程圖。
茲有關本發明之技術內容及詳細說明,配合圖式說明如下:請參閱圖1為本發明粉塵監控設備實施方式的示意圖。粉塵監控設備100包括粉塵濃度計1與粉塵濃度訊號處理裝置2,且粉塵濃度計1耦接粉塵濃度訊號處理裝置2。排煙道200具有特定風速(例如但不限於3米、40米等風速)由A點往B點的方向吹動,以將空氣中的粉塵D由A點吹至B點。粉塵濃度計1包括探棒12與感測電路14,且探棒12插設於排煙道200。其中,探棒12插入排煙道200的深度以1/3~2/3的排煙道截面直徑為較佳。當粉塵D(帶電粒子)與探棒12接觸、撞擊或摩擦時,在探棒12上會產生電荷。因此,可通過感測電路14與粉塵濃度訊號處理裝置2將感應的電荷進行濾波、放大、分析及處理等程序,即可獲得粉塵濃度Cd。
進一步而言,由於風往探棒12吹時,實體的探棒12會對流體產生干擾,容易產生懸浮粒子分布不均或渦流現象,容易使所偵測到的訊號失真或產生不必要的雜訊(例如粉塵D反覆接觸、撞擊或摩擦探棒12或懸浮粒子恰巧繞過探棒12的狀況)。因此本發明之主要目的及功效在於,提供排煙道總懸浮粒子濃度標定訊號處理方法,其包含非同步訊號的線上取樣特徵與處理方法,對於流體渦流(Vortex Shedding)的變化響應,軟體訊號的取樣頻率、訊號做小波濾波或移動平均等處理,可調整1秒至24小時做平均後再輸出,以改善排煙道渦流產生的雜訊影響,且提供較為精確的粉塵濃度Cd。
請參閱圖2為本發明感測電路的電路方塊示意圖,復配合參閱圖1。感測電路14主要係由硬體電路所組成,且包括放大電路142、濾波電路144及類比/數位轉換器146。放大電路142耦接探棒12,且將探棒12上所產生電荷(交流反應,主要係因應懸浮粒子的多寡產生振幅的變化)放大為放大訊號Sa。放大電路142可以包括電荷放大器142-1與可程式放大器142-2,且電荷放大器142-1耦接可程式放大器142-2。電荷放大器142-1主要是對探棒12上所產生電荷進行訊號放大以及初階的濾波,且可程式放大器142-2具有良好的訊號雜訊(S/N)比,以特別針對所需要的訊號進行放大,降低雜訊的影響。濾波電路144耦接可程式放大器142-2,且將放大訊號Sa濾波為濾波訊號Sf。濾波電路144可以使用三階以上的濾波器衰減率,每倍頻從原本衰減6dB提升到18dB以上,以提高訊號雜訊比。其中,濾波電路144較佳的還增加了增加50-60Hz的濾波能力,改善一般供電造成的幅射干擾影響。類比/數位轉換器146耦接濾波電路144,且將濾波訊號Sf轉換為偵測訊號Sd。值得一提,於本發明之一實施例中,上述感測電路14的內部結構,僅為眾多感測電路14的其中一種,並不以此為限。因此,舉凡可用以感測懸浮粒子而相應的產生偵測訊號Sd的感測電路14,皆應包含在本實施例之範疇當中。
請參閱圖3為本發明粉塵濃度訊號處理裝置的電路方塊示意圖,復配合參閱圖1~2。粉塵濃度訊號處理裝置2可以為微控制器(MCU)等具有資料處理功能的控制器,但不排除部分可以使用硬體電路所組成具有相同功能的實體控制電路(例如但不限於由比較器、放大器及邏輯閘所組成的電路)。粉塵濃度訊號處理裝置2包括特徵處理模組22、運算模組24、校正補償模組26及數值轉換模組28,且上述模組可以由軟體程式所組成,或者由部分軟體搭配部分硬體
所組成。特徵處理模組22耦接類比/數位轉換器146而接收偵測訊號Sd,以對偵測訊號Sd進行取樣而計算複數個特徵值Es,且基於這些特徵值Es判斷偵測訊號Sd的訊號類型。
具體地,特徵處理模組22會基於預設的取樣率來對偵測訊號Sd進行取樣,合適的取樣率範圍在10Hz~10kHz為較佳,尤其又以100Hz為最佳。具體地,若取樣太低,則擷取的信號因時間間隔太久,除了會漏補捉實際排煙道200中瞬間發生粉塵提高的情形,還會造成即時性降低、誤差變大等狀況。若取樣太高,則會將排煙道200中亂流產生的特徵都取樣到,粉塵濃度訊號處理裝置2會無法去判斷出哪些訊號是真實的或是因亂流造成的。因此,使用合適的取樣率,可以排除取樣到排煙道200內部亂流的特徵,也有足夠的取樣訊號進行訊號的分析。
特徵處理模組22包括取樣單元222、特徵計算單元224及類型判斷單元226,且還可以選擇性的包括濾波單元228。取樣單元222耦接類比/數位轉換器146,且對偵測訊號Sd進行取樣而提供複數個取樣值Vs。特徵計算單元224耦接取樣單元222,且基於取樣值Vs計算特徵值Es。具體地,粉塵濃度訊號處理裝置2可以預設特定的單位時間(例如但不限於1秒),特徵計算單元224分別計算每個單位時間中的取樣值Vs來獲得一個特徵值Es,以於連續的單位時間(例如但不限於1分鐘、1小時等)形成複數個特徵值Es。其中,特徵計算單元224包括平均值計算法、均方根值計算法、標準差計算法及峰對峰值計算法,且其公式如下:平均值計算法如式1所示:
峰對峰值計算法如式5所示:Vpp=[陣列A(最大值)]-[陣列A(最小值)]...(5)
特徵計算單元224係可基於取樣值Vs的態樣(例如但不限於,最大值與最小值較為明顯可選擇峰對峰值計算法及/或標準差計算法),選擇平均值計算法、該均方根值計算法、該標準差計算法及該峰對峰值計算法的至少其中之一將取樣值Vs計算為特徵值Es。
類型判斷單元226耦接特徵計算單元224,且基於特徵值Es的特徵判斷取樣值Vs屬於粉塵連續排放或粉塵間歇性排放。具體地,如圖4A與4B所示分別為粉塵連續排放或粉塵間歇性排放的特徵示意圖。粉塵連續排放通常會出現在例如但不限於,火力發電廠、垃圾焚化廠、造紙等連續燃燒排放粉塵場所,粉塵的排放通常是連續性的,僅有在設備不工作時會暫時的降低粉塵濃度。反之,粉塵間歇性排放通常會出現在工廠等場所,粉塵的排放通常是非連續性的,僅有在例如但不限於,大宗物料(塑化、糧食、混凝土、鐵沙)車輛廢氣排放、或施工時等狀況,粉塵濃度才會飆高。
另外一方面,濾波單元228耦接取樣單元222與特徵計算單元224之間,且基於取樣值Vs的雜訊大小,選擇性的對該些取樣值Vs進行濾波。具
體地,如圖5所示為濾波單元濾波波形示意圖。當取樣值Vs所構成的波形如曲線C1所示時,代表即便經過感測電路14的實體電路濾波後,仍然有高頻諧波。因此,濾波單元228主要係以小波濾波法,對曲線C1分別進行一至五階的濾波而分別產生C2~C6的曲線。小波濾波法主要是用標準差的方式,將曲線C1上的高頻諧波濾除,預估將取三階小波濾波法做為上限值(即曲線C4)。若進行太多階的濾波,則曲線將會失真(如曲線C5~C6所示)。然而,若是經過感測電路14的實體電路濾波後,高頻諧波較少,甚至是未有高頻諧波時(即感測電路14的功能較好),則可以不需要另行使用濾波單元228對取樣值Vs所構成的波形如曲線進行濾波。
運算模組24耦接類型判斷單元226,且基於特徵值Es屬於粉塵連續排放或粉塵間歇性排放(即訊號類型)選擇特定運算法,使特徵值Es通過特定運算法計算出粉塵濃度估算值Cc。具體地,運算模組24包括第一運算單元242與第二運算單元244,且第一運算單元242與第二運算單元244耦接類型判斷單元226。當特徵值Es屬於粉塵連續排放時,第一運算單元242提供加權移動平均運算,以通過加權移動平均運算將特徵值Es計算為粉塵濃度估算值Cc。其中,類型判斷單元226例如但不限於,可以基於特徵值Es屬於粉塵連續排放或粉塵間歇性排放(即訊號類型)提供相應的代碼給運算模組24,以使運算模組24可以選擇第一運算單元242或第二運算單元244進行後續的運算。加權移動平均運算如式6~7所示:F t =w 1 A t-1 +w 2 A t-2 +w 3 A t-3 +…+w n A t-n ...(6)
其中,Ft代表當下的平均值,Wi代表每一筆的權重(不管有幾筆,總合為1),越接近當下時間的值權重越大為較佳。At-n設定做移動平均時間段的粉塵濃度資料筆數(例如但不限於1秒100筆,即100Hz),At-1代表當下數據、At-2代表前1秒數據、At-3代表前2秒數據,依此類推。舉例而言,例如取5秒進行加權移動平均運算,數據如下表所示:
經計算結果為(8*0.05)+(11*0.05)+(15*0.15)+(12*0.25)+(10*0.5)=11.2,11.2即經計算後所獲得的粉塵濃度估算值Cc。
當特徵值Es屬於粉塵間歇性排放時,第二運算單元244提供積分運算,以通過積分運算將特徵值Es計算為粉塵濃度估算值Cc。積分運算如式8所示:i s =t 1 A 1+t 2 A 2+t 3 A 3+…+t n A n ...(8)
其中,is代表設定時間段進行積分計算面積,且An代表各個單位時間當下的粉塵濃度大小。tn代表單位時間,較佳的可以取濃度值(即特徵值Es)較為突出時,其前後的時間做為一區段。
校正補償模組26耦接第一運算單元242與第二運算單元244,且對粉塵濃度估算值Cc進行校正補償而獲得粉塵濃度精確值Ce。具體地,校正補償模組26包括校正單元262與補償單元264,校正單元262耦接第一運算單元242與第二運算單元244,且補償單元264耦接校正單元262。校正單元262通過校正常數值對粉塵濃度估算值Cc進行計算而獲得粉塵濃度校正值Ccv。進一步而言,如圖6所示為校正測試曲線圖。校正單元262係通過多點校正做內差法來計算粉塵濃度校正值Ccv,主要是校正單元262係於實際運作前,預先於測試環境(例如但不限於,標準實驗室)進行校正測試而獲得校正常數值,其公式如下式9所示:
校正後可取得(x1,y1)與(x2,y2)的標定點,區間的範圍為類線性的特性,將校正後的標定點代入公式9,取得校正常數值x代表粉塵濃度估算值Cc,其為校正單元262實際運作時所計算出的實際值(單位是V),代入上述公式就可以轉成校正後的粉塵濃度值(即為粉塵濃度校正值Ccv,單位是mg/m3)。
補償單元264通過環境常數值對粉塵濃度校正值Ccv進行計算而獲得粉塵濃度精確值Ce,具體地為補償單元264係基於測試環境與實際運作環境的環境差異來設定環境常數值。例如但不限於,實際運作的環境只要與測試環境(即標準實驗室)相同,代表Kt值(即環境常數值)為1。環境常數值的組成包含粉塵粉徑差異、排煙道200形狀差異、管內風速、溫濕度、探棒12深入長度、實際現場的工安因素...等。環境常數值組成因環境差異影響,有可能為以下任一關係的組合但不限定(K1、K2、K3到Kn代表上述因素)所組成的公式10~14:K t =K 1+K 2+K 3…+K n ...(10)
K t =|K 1-K 2-K 3…-K n |...(11)
K t =K 1×K 2×K 3…×K n ...(12)
通過上述公式所取得的環境常數值與粉塵濃度校正值Ccv進行計算,即可獲得粉塵濃度精確值Ce。
配合參閱圖7為本發明數值轉換模組的電路方塊圖。數值轉換模組28耦接補償單元264,且對粉塵濃度精確值Ce進行數值轉換,以提供特定資料格式且較為精確的粉塵濃度Cd。具體而言,數值轉換模組28包括輸出處理單元282與格式轉換單元284。輸出處理單元282耦接補償單元264,且格式轉換單元284耦接輸出處理單元282。輸出處理單元282包括粉塵濃度百分比(%)顯示公式計算與輸出電流訊號公式計算,且格式轉換單元284通過格式轉換,可提供顯示輸出格式Fd的塵濃度精確值Ce、類比輸出格式Fa的塵濃度精確值Ce、繼電器輸出格式Fr的塵濃度精確值Ce或RS485輸出格式Frs的塵濃度精確值Ce。數值轉換模組28通過上述二者計算,再通過特定資料格式的資料轉換,將塵濃度精確值Ce轉換為特定格式的粉塵濃度Cd。值得一提,由於本發明的運算模組24可以利用加權移動平均運算或積分運算對一段特定時段(例如但不限於當前時間至先前12小時之間的時段)的粉塵濃度進行計算,因此粉塵濃度訊號處理裝置2可以提供非同步輸出的粉塵濃度Cd。
具體地,粉塵濃度百分比(%)顯示公式計算係依設定的量測範圍做計算(例如0.1~A mg/m3),且粉塵濃度百分比(%)=(即時值/A)x 100%。其中,如圖8所示為量測範圍曲線圖。粉塵濃度百分比(%)顯示公式計算係依設定的量
測範圍做計算(例如0.1~A mg/m3),且即時值即為粉塵濃度精確值Ce。A代表粉塵濃度訊號處理裝置2內設定的量測範圍,例如圖8的點位C1~C2分別為量測範圍的上下限值。舉例而言,C1~C2為0.1mg/m3~10mg/m3,則C1為0.1且C2為10,或0.1mg/m3~500mg/m3,則C1為0.1且C2為500(此範圍可設定)。輸出電流訊號(mA)公式計算符合業界40mA~20mA的輸出格式,其主要是將粉塵濃度精確值Ce轉換為類比的輸出電流訊號。具體地,輸出電流訊號公式計算係依設定的量測範圍做計算(例如0.1~A mg/m3),且輸出電流訊號(mA)=[(即時值/A)x16]+4。其中,即時值即為粉塵濃度精確值Ce,且A代表粉塵濃度訊號處理裝置2內設定的量測範圍(即4~20)。
格式轉換單元284可包括顯示轉換單元284-1、類比轉換單元284-2、閾值比較單元284-3及RS485轉換單元284-4,且上述四種不同的轉換單元耦接輸出處理單元282與後端所耦接的裝置(例如但不限於LCM螢幕、電腦、繼電器等)。顯示轉換單元係將輸出處理單元282所轉換出的數值轉換為顯示輸出格式Fd的粉塵濃度Cd,以提供較為精確的粉塵濃度Cd給液晶模組(LCM)顯示。類比轉換單元係將輸出處理單元282所轉換出的數值轉換為類比輸出格式Fa的粉塵濃度Cd,以提供較為精確的粉塵濃度Cd給後端的裝置。其中,類比轉換單元主要係將數值轉換為4mA~20mA的類比輸出電流訊號。
閾值比較單元係將輸出處理單元282所轉換出的數值轉換為繼電器輸出格式Fr的粉塵濃度Cd,以提供較為精確的粉塵濃度Cd給繼電器去做導通/關斷。即當粉塵濃度Cd大於閾值,則繼電器作動而表示粉塵濃度Cd過高,反之則繼電器不作動而表示粉塵濃度Cd未超標。RS485轉換單元係將輸出處理單元282所轉換出的數值轉換為RS485輸出格式Frs的粉塵濃度Cd,以提供較
為精確的粉塵濃度Cd給後端的裝置。值得一提,於本發明之一實施例中,數值轉換模組28可輸出的格式並不限定必須要上述四種,依實際需求可取上述四種的其中之一種輸出,亦或是額外加入其他種的輸出格式。此外,上述的四種輸出方式除了可提供非同步的粉塵濃度Cd外,還可以提供當下即時的粉塵濃度Cd,使使用者可以直觀的知悉當下的數值與非同步的數值差異。
請參閱圖9為本發明粉塵濃度訊號處理方法的方法流程圖,復配合參閱圖1~8。粉塵濃度訊號處理方法主要是對粉塵濃度計1的偵測訊號Sd進行取樣而計算較為精確的粉塵濃度Cd,以改善排煙道渦流產生的雜訊影響。具體地,粉塵濃度訊號處理方法包括,對偵測訊號進行取樣(S100)。利用較佳的實施方式為,特徵處理模組22會基於預設的取樣率來對偵測訊號Sd進行取樣,且合適的取樣率範圍在10Hz~10kHz為較佳,尤其又以100Hz為最佳。然後,是否進行濾波(S120)。當結果為”是”時,對取樣值進行濾波(S140)。較佳的實施方式為,利用濾波單元228主要係以小波濾波法,對取樣值Vs所構成的曲線C1進行(一至五階)的濾波,以曲線C1上的高頻諧波濾除。當結果為”否”時,則基於取樣值來計算特徵值(S160)。較佳的實施方式為,利用粉塵濃度訊號處理裝置2預設特定的單位時間(例如但不限於1秒),再利用特徵計算單元224分別計算每個單位時間中的取樣值Vs來獲得一個特徵值Es,以於連續的單位時間(例如但不限於1分鐘、1小時等)形成複數個特徵值Es。其中,若高頻諧波較少,甚至是未有高頻諧波時,則步驟(S120)~(S140)可以省略。
然後,判斷排放類型(S180)。較佳的實施方式為,利用類型判斷單元226基於特徵值Es的特徵判斷取樣值Vs屬於粉塵連續排放或粉塵間歇性排放。然後,提供加權移動平均運算(S200)。較佳的實施方式為,當取樣值Vs屬
於粉塵連續排放時,通過第一運算單元242提供加權移動平均運算,以通過加權移動平均運算將特徵值Es計算為粉塵濃度估算值Cc。然後,提供積分運算(S220)。較佳的實施方式為,當取樣值Vs屬於粉塵間歇性排放時,通過第二運算單元2442提供加積分運算,以通過積分運算將特徵值Es計算為粉塵濃度估算值Cc。
然後,對粉塵濃度估算值進行校正(S240)。較佳的實施方式為,利用校正單元262通過多點校正做內差法來計算粉塵濃度校正值Ccv。主要是在實際運作前,粉塵濃度訊號處理裝置2預先於測試環境(例如但不限於,標準實驗室)進行校正測試而獲得校正常數值,然後再利用校正常數值與粉塵濃度估算值Cc計算粉塵濃度校正值Ccv。然後,對粉塵濃度校正值進行補償(S260)。較佳的實施方式為,利用補償單元264通過環境常數值對粉塵濃度校正值Ccv進行計算,以獲得粉塵濃度精確值Ce。具體地,補償單元264係基於測試環境與實際運作環境的環境差異來設定環境常數值,且環境常數值與粉塵濃度校正值Ccv進行計算,即可獲得粉塵濃度精確值Ce。
最後,進行數值轉換而獲得粉塵濃度(S280)。較佳的實施方式為,利用數值轉換模組28對粉塵濃度精確值Ce進行計算與格式的轉換,以將塵濃度精確值Ce轉換為特定格式的粉塵濃度Cd。其中,數值轉換模組28至少可將粉塵濃度精確值Ce轉換為顯示輸出格式Fd、類比輸出格式Fa、繼電器輸出格式Fr或RS485輸出格式Frs的至少其中之一。值得一提,於本發明之一實施例中,上述的流程並非僅能使用圖3所出示的模組來達成,舉凡可達成相同功能之模組/電路/程式等,皆應包含在本實施例之範疇當中。
以上所述,僅為本發明較佳具體實施例之詳細說明與圖式,惟本發明之特徵並不侷限於此,並非用以限制本發明,本發明之
所有範圍應以下述之申請專利範圍為準,凡合於本發明申請專利範圍之精神與其類似變化之實施例,皆應包含於本發明之範疇中,任何熟悉該項技藝者在本發明之領域內,可輕易思及之變化或修飾皆可涵蓋在以下本案之專利範圍。
(S100)~(S280):步驟
Claims (12)
- 一種粉塵濃度訊號處理方法,係對一粉塵濃度計的一偵測訊號進行訊號分析而計算出一粉塵濃度,且該粉塵濃度訊號處理方法包括下列步驟:對該偵測訊號進行取樣而計算複數個特徵值,且基於該些特徵值判斷該偵測訊號的一訊號類型;基於該些特徵值選擇一特定運算法,且該些特徵值通過該特定運算法計算出一粉塵濃度估算值;對該粉塵濃度估算值進行一校正補償而獲得一粉塵濃度精確值;及對該粉塵濃度精確值進行一數值轉換,以提供一特定資料格式的該粉塵濃度;對該偵測訊號以一取樣率進行取樣而提供複數個取樣值;分別計算每個單位時間中的該些取樣值來獲得一特徵值,以於連續的單位時間形成該些特徵值;及基於該些特徵值的一特徵判斷,該些取樣值屬於一粉塵連續排放或一粉塵間歇性排放。
- 如請求項1所述之粉塵濃度訊號處理方法,該粉塵濃度訊號處理方法更包括下列步驟:基於該些取樣值的一雜訊大小,選擇性的對該些取樣值進行一濾波處理程序。
- 如請求項1所述之粉塵濃度訊號處理方法,該粉塵濃度訊號處理方法更包括下列步驟: 選擇平均值計算法、該均方根值計算法、該標準差計算法及該峰對峰值計算法的至少其中之一將該些取樣值計算為該些特徵值。
- 如請求項1所述之粉塵濃度訊號處理方法,該粉塵濃度訊號處理方法更包括下列步驟:基於該些特徵值屬於該粉塵連續排放而提供一加權移動平均運算,以將該些特徵值計算為該粉塵濃度估算值;或基於該些特徵值屬於該粉塵間歇性排放而提供一積分運算,以將該些特徵值計算為該粉塵濃度估算值。
- 如請求項1所述之粉塵濃度訊號處理方法,該粉塵濃度訊號處理方法更包括下列步驟:通過一校正常數值對該粉塵濃度估算值進行計算而獲得一粉塵濃度校正值;通過一環境常數值對該粉塵濃度校正值進行計算而獲得該粉塵濃度精確值;於實際運作前,預先於一測試環境進行一校正測試而獲得該校正常數值;及基於該測試環境與實際運作的環境的一環境差異來設定該環境常數值。
- 一種粉塵濃度訊號處理裝置,係對一粉塵濃度計的一偵測訊號進行取樣而計算一粉塵濃度,且該粉塵濃度訊號處理裝置包括:一特徵處理模組,對該偵測訊號進行取樣而計算複數個特徵值,且基於該些特徵值判斷該偵測訊號的一訊號類型,該特徵處理模組包括:一取樣單元,對該偵測訊號以一取樣率進行取樣而提供複數個取樣值; 一特徵計算單元,基於該些取樣值計算該些特徵值;及一類型判斷單元,基於該些特徵值的一特徵判斷該偵測訊號屬於一粉塵連續排放或一粉塵間歇性排放;一運算模組,基於該些特徵值選擇一特定運算法,且該些特徵值通過該特定運算法計算出一粉塵濃度估算值;一校正補償模組,對該粉塵濃度估算值進行一校正補償而獲得一粉塵濃度精確值;及一數值轉換模組,對該粉塵濃度精確值進行一數值轉換,以提供一特定資料格式的該粉塵濃度;其中,該特徵計算單元分別計算每個單位時間中的該些取樣值來獲得一特徵值,以於連續的單位時間形成該些特徵值。
- 如請求項6所述之粉塵濃度訊號處理裝置,其中該特徵處理模組更包括:一濾波單元,基於該些取樣值的一雜訊大小,選擇性的對該些取樣值進行一濾波。
- 如請求項6所述之粉塵濃度訊號處理裝置,其中該特徵計算單元包括一平均值計算法、一均方根值計算法、一標準差計算法及一峰對峰值計算法,該特徵計算單元選擇平均值計算法、該均方根值計算法、該標準差計算法及該峰對峰值計算法的至少其中之一將該些取樣值計算為該些特徵值。
- 如請求項6所述之粉塵濃度訊號處理裝置,其中該運算模組包括: 一第一運算單元,基於該些特徵值屬於該粉塵連續排放而提供一加權移動平均運算,以將該些特徵值計算為該粉塵濃度估算值;及一第二運算單元,基於該些特徵值屬於該粉塵間歇性排放而提供一積分運算,以將該些特徵值計算為該粉塵濃度估算值。
- 如請求項6所述之粉塵濃度訊號處理裝置,其中該校正補償模組包括:一校正單元,通過一校正常數值對該粉塵濃度估算值進行計算而獲得一粉塵濃度校正值;及一補償單元,通過一環境常數值對該粉塵濃度校正值進行計算而獲得該粉塵濃度精確值。
- 如請求項10所述之粉塵濃度訊號處理裝置,其中該校正單元係於實際運作前,預先於一測試環境進行一校正測試而獲得該校正常數值。
- 如請求項11所述之粉塵濃度訊號處理裝置,其中該補償單元係基於該測試環境與實際運作的環境的一環境差異來設定該環境常數值。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW111121251A TWI802434B (zh) | 2022-06-08 | 2022-06-08 | 粉塵濃度訊號處理裝置及其訊號處理方法 |
DE102022121355.2A DE102022121355A1 (de) | 2022-06-08 | 2022-08-24 | Staubkonzentrationssignal-Verarbeitungsvorrichtung und Signalverarbeitungsverfahren dafür |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW111121251A TWI802434B (zh) | 2022-06-08 | 2022-06-08 | 粉塵濃度訊號處理裝置及其訊號處理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
TWI802434B true TWI802434B (zh) | 2023-05-11 |
TW202348971A TW202348971A (zh) | 2023-12-16 |
Family
ID=87424379
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW111121251A TWI802434B (zh) | 2022-06-08 | 2022-06-08 | 粉塵濃度訊號處理裝置及其訊號處理方法 |
Country Status (2)
Country | Link |
---|---|
DE (1) | DE102022121355A1 (zh) |
TW (1) | TWI802434B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019074285A1 (ko) * | 2017-10-12 | 2019-04-18 | 주식회사 엑스엘 | 미세먼지 농도 연산장치 및 방법 |
TW202045909A (zh) * | 2019-02-01 | 2020-12-16 | 日商夏普股份有限公司 | 微粒子檢測感測器、灰塵感測器、空調設備、及微粒子檢測感測器的控制方法 |
TW202206803A (zh) * | 2020-08-03 | 2022-02-16 | 思維環境科技有限公司 | 提高環境感測裝置檢測精度的校正方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102009000904A1 (de) | 2008-03-03 | 2009-09-10 | Agilent Technologies Inc., Santa Clara | Verfahren und Systeme zum Berechnen einer Grössenverteilung von kleinen Partikeln |
GB201217228D0 (en) | 2012-09-26 | 2012-11-07 | Pepric Nv | Methods and systems for determining a particle distribution |
-
2022
- 2022-06-08 TW TW111121251A patent/TWI802434B/zh active
- 2022-08-24 DE DE102022121355.2A patent/DE102022121355A1/de active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019074285A1 (ko) * | 2017-10-12 | 2019-04-18 | 주식회사 엑스엘 | 미세먼지 농도 연산장치 및 방법 |
TW202045909A (zh) * | 2019-02-01 | 2020-12-16 | 日商夏普股份有限公司 | 微粒子檢測感測器、灰塵感測器、空調設備、及微粒子檢測感測器的控制方法 |
TW202206803A (zh) * | 2020-08-03 | 2022-02-16 | 思維環境科技有限公司 | 提高環境感測裝置檢測精度的校正方法 |
Also Published As
Publication number | Publication date |
---|---|
DE102022121355A1 (de) | 2023-12-14 |
TW202348971A (zh) | 2023-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105334147A (zh) | 基于β射线法和光散射法的颗粒物在线监测系统及方法 | |
CN101113947A (zh) | 测量车辆颗粒物排放的装置及其方法 | |
CN108827843B (zh) | 固定污染源颗粒物的质量浓度和粒度谱的测量装置及方法 | |
CN100460831C (zh) | 微电荷型在线粉尘及探头污染检测的方法 | |
CN103033399A (zh) | 一种pm2.5颗粒物监测仪用大气多参数变送器及其使用方法 | |
Hu et al. | On-line sizing of pneumatically conveyed particles through acoustic emission detection and signal analysis | |
CN107132166A (zh) | 一种直插式粉尘浓度测量仪及测量方法 | |
CN105488772B (zh) | 一种传感器信号峰值检测方法 | |
CN109443445A (zh) | 一种颗粒物在线监测装置与方法 | |
CN106323830A (zh) | 一种高湿烟气粉尘监测系统及其方法 | |
CN202661356U (zh) | 空气微小颗粒物监测仪器 | |
CN201788149U (zh) | 具有湿度校正的光散射式颗粒物在线监测仪 | |
TWI802434B (zh) | 粉塵濃度訊號處理裝置及其訊號處理方法 | |
CN105527203A (zh) | 烟尘浓度直读监测系统及方法 | |
CN205027638U (zh) | 一种pm2.5工业烟气的在线监测装置 | |
CN104459197B (zh) | 降低现场管道气体流速测量下限的方法及系统 | |
CN110220558A (zh) | 一种气固两相流参数检测方法 | |
CN102095773B (zh) | 烟气含氧量测量方法及系统 | |
CN105466826A (zh) | 颗粒物在线监测系统及方法 | |
CN210833647U (zh) | 一种输粉管内粉体速度、浓度和细度在线测量系统 | |
CN107764889B (zh) | 一种基于信号能量法的气固两相流粉尘含碳量测量方法 | |
CN106644867B (zh) | 气体中颗粒物的检测装置及方法 | |
CN109253953B (zh) | 一种可测量多种颗粒物质量浓度的方法和系统 | |
CN201795972U (zh) | 在线粉尘检测仪 | |
CN117309700A (zh) | 粉尘浓度信号处理装置及其信号处理方法 |