唐 琳,廖先莉,劉星月,趙永鑫,李躍鵬,余松科,3
1. 成都大學(xué)信息科學(xué)與工程學(xué)院,四川 成都 610106 2. 成都大學(xué)模式識別與智能信息處理四川省高校重點(diǎn)實(shí)驗(yàn)室,四川 成都 610106 3. 數(shù)學(xué)地質(zhì)四川省重點(diǎn)實(shí)驗(yàn)室,成都理工大學(xué),四川 成都 610059
能量色散X熒光光譜測量系統(tǒng)通常由半導(dǎo)體探測器、前端電子學(xué)電路以及后端數(shù)字信號處理電路組成。硅漂移探測器以其較高的能量分辨率和較好的性能被廣泛應(yīng)用[1-4]。探測器輸出信號幅度較小,為了便于傳輸絕大部分測量系統(tǒng)都需要有相應(yīng)的前端電子學(xué)電路對探測器輸出的電信號進(jìn)行放大,從而得到更高的能量分辨率[5-8],而開關(guān)復(fù)位型前放以其良好的噪聲新能被廣泛采用。該類型的前放輸出信號是一系列幅度不斷堆積上升的階躍脈沖,上升到一定程度當(dāng)開關(guān)復(fù)位時就跳變到零并開始下一輪堆積。由于開關(guān)復(fù)位的時間是不固定的,因此每次脈沖跳變的時刻也是隨機(jī)的,并且每一次跳變都有可能產(chǎn)生一個保持時間不夠的脈沖,這樣的脈沖在經(jīng)過CR微分電路轉(zhuǎn)換、后端放大電路放大、數(shù)字脈沖成形之后將會產(chǎn)生一個脈沖幅度和寬度受損的三角脈沖[5],本文將這樣的脈沖定義為突變脈沖。在前期的研究結(jié)果中,我們已經(jīng)證實(shí)如果不對突變脈沖進(jìn)行處理,將會在特征峰前面形成一個影子峰,也提出了通過剔除突變脈沖來消除影子峰的脈沖剔除技術(shù)[9]。經(jīng)過進(jìn)一步研究,我們發(fā)現(xiàn)剔除突變脈沖的方法雖然可以達(dá)到消除影子峰的目的,但也存在損失計(jì)數(shù)率的缺陷。
以往對計(jì)數(shù)率校正的研究大多是從死時間等方面著手,國內(nèi)外的研究學(xué)者們也對在計(jì)數(shù)率研究和校正領(lǐng)域取得了一些研究成果[10-12]??紤]到影子峰的本質(zhì)是計(jì)數(shù)率較高的元素特征峰因脈沖幅度受損造成的特征峰漂移,也可以稱為特征峰的影子峰。本文提出兩種特征峰漂移校正算法對突變負(fù)指數(shù)脈沖進(jìn)行修復(fù),該算法的可行性和準(zhǔn)確性通過模擬和實(shí)驗(yàn)進(jìn)行驗(yàn)證,結(jié)果表明特征峰漂移校正算法可以有效地將影子峰區(qū)域的計(jì)數(shù)率校正到元素特征峰中,既消除了影子峰,也保證了計(jì)數(shù)率不受損失,對獲取精細(xì)X射線譜具有重大意義。
開關(guān)復(fù)位型前置放大電路(下文簡稱前放)及其后端電路,如圖1所示,電路工作原理實(shí)際上就是探測器輸出的弱電流在反饋電容Cf上進(jìn)行一定時間的積分,并將積分結(jié)果存儲在采樣保持電容上由測量電路進(jìn)行測量。
圖1 前放電路以及其輸出信號Fig.1 Preamplifier circuit and its output signal
由于前放輸出信號是不斷上升且頻繁跳變的,如果直接對該信號進(jìn)行放大、數(shù)字成形將會造成脈沖幅度溢出,損失計(jì)數(shù)率。因此FPGA對前放輸出的階躍信號進(jìn)行脈沖成形處理時需要先將該信號通過CR微分電路轉(zhuǎn)變?yōu)樨?fù)指數(shù)信號。如果前放輸出的最后一個階躍脈沖剛好是一個脈沖寬度不夠的畸形脈沖,對這樣的脈沖進(jìn)行CR微分處理得到的負(fù)指數(shù)脈沖也將會是不完整的。對這樣的脈沖進(jìn)行三角成形得到的結(jié)果如圖2中pulse2所對應(yīng)的成形結(jié)果所示,可以看出成形后的脈沖寬度和脈沖幅度都受到了影響。
圖2 影子峰形成原理Fig.2 Principle of shadow peaks formation
如前文所述,經(jīng)CR微分電路輸出的負(fù)指數(shù)信號在放大和數(shù)字化后再進(jìn)行數(shù)字三角成形,當(dāng)負(fù)指數(shù)信號在下降沿出現(xiàn)了突變時,最終成形結(jié)果的脈沖幅度和寬度也會受到損失。
本文所采用的數(shù)字多道譜儀處理的核脈沖信號經(jīng)主放電路放大之后幅度范圍在0~2 V,對應(yīng)2 048個道址,在模擬的時候可以近似認(rèn)為脈沖幅度每增加1 mV則道址數(shù)加一[5]。當(dāng)出現(xiàn)圖2中pulse2這樣的突變脈沖時,該脈沖成形后的幅度大幅降低,進(jìn)行多道成譜時,原本元素特征峰所處的道址范圍應(yīng)該在1 000左右,但突變脈沖因幅度受損導(dǎo)致特征峰的計(jì)數(shù)左移,當(dāng)突變脈沖數(shù)量較大時就會在特征峰的前面形成一個影子峰。因此可以得出,影子峰的本質(zhì)其實(shí)就是特征峰中一部分計(jì)數(shù)率左移。由于總的脈沖數(shù)量是恒定的,所以當(dāng)特征峰區(qū)域的計(jì)數(shù)往左漂移時,影子峰計(jì)數(shù)升高必然也就會導(dǎo)致元素特征峰的計(jì)數(shù)降低,這對我們獲取精確的X射線譜是不利的,因此在下文中將提出兩種特征峰漂移校正方法將影子峰區(qū)域中的計(jì)數(shù)盡可能校正到特征峰中,既可以消除影子峰,也能夠達(dá)到提升特征峰計(jì)數(shù)的目的。
特征峰漂移產(chǎn)生的根源在于突變脈沖損失了部分采樣點(diǎn)導(dǎo)致三角成形結(jié)果的幅度受損,因此特征峰漂移校正的實(shí)質(zhì)就是對突變脈沖的修復(fù)。突變脈沖的修復(fù)在FPGA中實(shí)現(xiàn),其算法的核心在于修復(fù)條件和修復(fù)算法的確定??紤]到脈沖突變部分所有采樣點(diǎn)瞬間跳變成零的特點(diǎn),本文對需要修復(fù)的采樣點(diǎn)的定位采用判零法。在所有脈沖序列中,只要出現(xiàn)為零的采樣點(diǎn)都用校正的修復(fù)方法對其進(jìn)行修復(fù),下一個點(diǎn)若還是為零則繼續(xù)迭代,直到出現(xiàn)不為零的采樣點(diǎn)為止,完成修復(fù)后再對負(fù)指數(shù)脈沖序列進(jìn)行三角成形。而對修復(fù)算法的選擇,本文提出快校正和慢校正兩種方法,當(dāng)修復(fù)出的曲線衰減速度很快時將特征峰漂移校正方法稱之為快校正,當(dāng)衰減速度較慢并且接近于原始脈沖本身的衰減速度時,將這種校正方法稱之為慢校正,下文將對快校正和慢校正兩種方法得到的測量結(jié)果進(jìn)行詳細(xì)的介紹。
快校正的關(guān)鍵在于選擇最恰當(dāng)?shù)膬蓚€采樣點(diǎn)并根據(jù)其坐標(biāo)計(jì)算出修復(fù)直線的表達(dá)式。依然假定負(fù)指數(shù)信號v(n)損失了n+1之后的所有采樣點(diǎn),根據(jù)v(n)的峰值采樣點(diǎn)坐標(biāo)和突變前最后一個非零采樣點(diǎn)的坐標(biāo)來確定一條修復(fù)直線,根據(jù)求解出的直線表達(dá)式計(jì)算出損失的采樣點(diǎn),迭代修復(fù)直到查詢到下一個非零的采樣點(diǎn)為止。
文獻(xiàn)[9]已經(jīng)詳細(xì)討論過只有脈沖突變的時刻在三角成形的上升時間tup之前時該脈沖成形結(jié)果的幅度和脈沖寬度才會受到影響,并在最終測量得到的譜圖上以影子峰形式呈現(xiàn)。因此,此處對突變脈沖的校正就選擇突變時刻在tup之前的。
在本章所指定的脈沖序列中,tup時刻對應(yīng)著第200個采樣點(diǎn)。通過上文所述的直線修復(fù)法的執(zhí)行過程,設(shè)修復(fù)直線的表達(dá)式為y=kx+b。首先取突變脈沖的峰值坐標(biāo)(101,2 000),再取突變前的最后一點(diǎn)的坐標(biāo)(148,1 250),根據(jù)兩點(diǎn)坐標(biāo)求出修復(fù)直線的表達(dá)式如式(1)所示。
y=-16x+3 616
(1)
由式(1)可恢復(fù)出突變損失掉的部分采樣點(diǎn),直線修復(fù)法的優(yōu)點(diǎn)是計(jì)算簡單,缺點(diǎn)就在于它的衰減速度快,只能修復(fù)出一部分采樣點(diǎn),當(dāng)修復(fù)結(jié)果到0以后,其余的采樣點(diǎn)也只能以0來補(bǔ)充,修復(fù)結(jié)果如圖3中紅色曲線所示,修復(fù)后突變脈沖的三角成形結(jié)果如圖3中的藍(lán)色曲線所示。
圖3 快校正的修復(fù)結(jié)果及成形結(jié)果Fig.3 Repaired result by fast correction and shaping result
從圖中可以看出,修復(fù)后的脈沖成形結(jié)果幅度依然有所損失,在三角成形末端有一個較大的俯沖,導(dǎo)致在大于零的區(qū)域內(nèi)三角成形的脈沖寬度小于完整的負(fù)指數(shù)脈沖成形結(jié)果的寬度,這樣的修復(fù)結(jié)果對最終測量譜圖的影響將在下文的實(shí)驗(yàn)驗(yàn)證環(huán)節(jié)中進(jìn)行詳細(xì)分析。
慢校正修復(fù)算法的關(guān)鍵在于修復(fù)結(jié)果的衰減趨勢趨近于原始脈沖的衰減趨勢,該算法通過負(fù)指數(shù)脈沖的表達(dá)式迭代實(shí)現(xiàn),具體執(zhí)行過程如下所述。
假定負(fù)指數(shù)信號的表達(dá)方式如式(2)所示
(2)
式(2)中,A為負(fù)指數(shù)脈沖的幅值,τ為衰減時間常數(shù),Tclk為采樣周期。
發(fā)生脈沖突變時,可視為負(fù)指數(shù)信號v(n)損失了n+1之后的所有采樣點(diǎn),以表達(dá)式(2)為基礎(chǔ),可推導(dǎo)出慢校正的修復(fù)表達(dá)式如式(3)所示
(3)
以式(3)為理論基礎(chǔ)可以修復(fù)出突變脈沖丟失的采樣點(diǎn),修復(fù)結(jié)果如圖4所示。未經(jīng)處理的原始突變脈沖如圖4的黑色曲線所示,在第256個采樣點(diǎn)時脈沖幅度直接從428跳變到零,調(diào)用慢校正算法對突變脈沖進(jìn)行修復(fù)后得到如圖4藍(lán)色曲線所示的慢校正修復(fù)結(jié)果。
圖4 慢校正的修復(fù)結(jié)果及成形結(jié)果Fig.4 Repaired result by slow correction and shaping result
通過圖3和圖4的對比可以得出,快校正修復(fù)后的曲線衰減太快,脈沖下降沿的衰減趨勢不符合原始脈沖的衰減趨勢,而慢校正修復(fù)后的曲線衰減趨勢幾乎與原始脈沖的衰減趨勢保持一致,修復(fù)效果較好。下文將通過實(shí)際樣品的測量結(jié)果對快校正和慢校正更加詳細(xì)的對比。
測量系統(tǒng)的結(jié)構(gòu)圖如圖5所示,探測器、數(shù)字脈沖處理器(DPP)和pc軟件。
圖5 測量系統(tǒng)結(jié)構(gòu)圖Fig.5 Structure chart of the detection system
探測器(XR-100SDD)是AMPTEK生產(chǎn)的高性能探測器。DPP主要包括三部分:前端電路、ADC和FPGA的脈沖處理單元。經(jīng)前端電路放大后的脈沖由ADC進(jìn)行采樣,ADC以20 Msps的速度工作,分辨率為12位。FPGA選用Xilinx公司生產(chǎn)的xc3s400芯片,對數(shù)字脈沖序列進(jìn)行逐一篩選。當(dāng)脈沖被標(biāo)記為突變脈沖時,調(diào)用校正算法對其進(jìn)行校正,反之,若未標(biāo)記則直接進(jìn)入下一個單元的處理流程。
探測器FAST-SDD的有效探測面積為25 mm2,厚度為500 μm,鈹窗為0.5 mil,激發(fā)源采用科頤維KYW2000A型X光管,ADC采樣頻率20 MHz,采樣周期50 ns。為了便于驗(yàn)證各個能量段元素的特征峰漂移校正,選取大邑豐山鐵礦作為測試對象,該樣品在各個能量段中都含有某種代表性元素,測量得到的X射線光譜作為測試譜,譜圖中包含Ti,F(xiàn)e,Sr和Sn等元素,其中主元素為Fe和Sr,其他元素成分含量比較低。下文將分別通過譜圖對比和計(jì)數(shù)分析兩種方式對不同特征峰漂移校正方法進(jìn)行驗(yàn)證。
本文提出的特征峰漂移校正方法就是通過修復(fù)突變脈沖來達(dá)到降低影子峰計(jì)數(shù),提升元素特征峰計(jì)數(shù)的目的。實(shí)驗(yàn)選取Fe元素和Sr元素特征峰及其影子峰所在的道址區(qū)間為分析對象,將采用不同特征峰漂移校正方法得到的測量結(jié)果與不校正得到的結(jié)果進(jìn)行對比。結(jié)果如圖6所示,以Fe元素特征峰為分析對象時,從圖6(a)的放大區(qū)域可以看出,快校正和慢校正兩種方法都能夠讓Feshadow的計(jì)數(shù)大幅度降低;以Sr元素特征峰為分析對象時,從圖6(b)的放大區(qū)域可以觀察出快校正得到的Sr元素特征峰底部被展寬。
圖6 特征峰漂移校正前后對照圖(a):Fe元素特征峰;(b):Sr元素特征峰Fig.6 Characteristic peak contrast map before and after correction(a):Characteristic peak of Fe;(b):Characteristic peak of Sr
從計(jì)數(shù)的角度上分析可以得出,由于快校正對突變脈沖衰減部分的修復(fù)不完整,導(dǎo)致修復(fù)后脈沖成形結(jié)果的幅度依然有一定的損失,從而造成Srshadow區(qū)域的計(jì)數(shù)轉(zhuǎn)移到了Sr元素特征峰底部,最終可以得出慢校正相比于快校正更容易得到精確的X射線光譜。
取Fe元素和Sr元素特征峰及其影子峰所在的道址區(qū)間為分析對象,對特征峰漂移校正前后每個區(qū)間的計(jì)數(shù)率之和進(jìn)行分析,如表1所示。
表1 特征峰漂移校正數(shù)據(jù)對比Table 1 Comparison of characteristic peak drift correction data
表1中,Corigin表示原始計(jì)數(shù)率,Cfast和Cslow分別表示采用快校正法、慢校正法校正特征峰漂移后得到的計(jì)數(shù)率,Dfast和Dslow表示快校正和慢校正前后的計(jì)數(shù)率差值,E表示校正效率,其計(jì)算公式如式(4)所示
Efast=Dfast/Corigin
Eslow=Dslow/Corigin
(4)
從表1可以看出,采用不同方法校正后Feshadow和Srshadow所在道址區(qū)間的計(jì)數(shù)相比于不校正時有了明顯的降低,與此同時,兩個特征峰所在道址區(qū)間的計(jì)數(shù)相比于不校正則有了明顯的提高。因此Feshadow和Srshadow在校正前后的計(jì)數(shù)率差值始終為負(fù),而元素特征峰的計(jì)數(shù)率差值則始終為正。理想情況下,同一種元素在影子峰區(qū)域計(jì)數(shù)率的減小值與在特征峰區(qū)域計(jì)數(shù)率的增加值在數(shù)值上應(yīng)該趨于一致,表1中Feshadow和鐵元素特征峰所在區(qū)間的Dfast和Dslow的值基本符合這一趨勢,但Srshadow和鍶元素特征峰所在區(qū)間的Dfast值相差較大,不符合影子峰計(jì)數(shù)減小值即為特征峰計(jì)數(shù)增加值的規(guī)律。根據(jù)前文的譜圖對比分析可知造成這種結(jié)果的根本原因在于快校正對突變脈沖的修復(fù)不完整,而慢校正可以較好地實(shí)現(xiàn)所有采樣點(diǎn)的修復(fù),最后得出的修復(fù)效率也表明對于同樣的區(qū)間,慢校正法得到的修復(fù)效率更高,對特征峰漂移的校正效果更好。
在實(shí)際測量中,采用慢校正的方法可以解決開關(guān)復(fù)位型前放頻繁復(fù)位產(chǎn)生突變脈沖,并因此造成的特征峰漂移問題。針對上述問題,本文提出特征峰漂移校正方法,該方法由電路描述、原理、方法和實(shí)驗(yàn)驗(yàn)證四個部分組成,方法部分得出式(3)作為慢校正脈沖修復(fù)的理論依據(jù),并采用MATLAB對不同校正方法得到的修復(fù)結(jié)果進(jìn)行模擬對比。實(shí)驗(yàn)驗(yàn)證部分以鐵礦樣品為測量對象對兩種校正方法做出了驗(yàn)證,實(shí)驗(yàn)結(jié)果表明雖然快校正和慢校正都可以有效地降低影子峰計(jì)數(shù),但是快校正在消除影子峰的同時也會導(dǎo)致特征峰底部展寬,而慢校正得到的測量結(jié)果中,影子峰區(qū)域被降低的計(jì)數(shù)幾乎都校正到了元素特征峰的計(jì)數(shù)中,有效地實(shí)現(xiàn)了對特征峰漂移的校正。