基于數(shù)據(jù)擬合的RM 215測厚儀測量曲線標定研究

王小波

（湖南華菱漣源鋼鐵有限公司信息化中心，湖南 婁底 417009）

金屬板帶加工過程中，板帶厚度和板帶尺寸控制是極為重要的控制目標，其中金屬板帶的厚度控制往往和板帶加工的質量直接掛鉤，具有高精度厚度控制的板帶產(chǎn)品不僅在市場競爭中處于有利地位，同時也可以增加產(chǎn)品的附加值。而板帶厚度控制中極為重要的一環(huán)就是加工過程中的板帶厚度檢測問題，精確的板帶厚度檢測是高精度厚度控制的基礎[1]。

X射線測厚儀由于具有非接觸、連續(xù)、快捷、無損傷檢測，且能夠適用于高油污、高溫惡劣環(huán)境的優(yōu)點，被廣泛應用于板帶軋制過程。但在實際使用過程中，X射線測厚儀的檢測精度受現(xiàn)場的X射線源穩(wěn)定性、合金成分、帶鋼位置、油污水分等的影響。因此需要對X射線的測量曲線進行校正或者補償[2]。

針對板帶加工測量過程中的X射線測厚儀補償和校正問題，為提高檢測的精度和穩(wěn)定性，一些學者和工程師做了相關的研究工作。文獻[3]系統(tǒng)地分析了X射線測厚儀射線強度、統(tǒng)計誤差以及被測物的材質、溫度、表面附著水層、傾斜角度等對測厚精度的影響并給出了理論上的精度修正方法。其后，在文獻[4]中給出了X射線測厚儀應用實例以及在國內(nèi)的應用情況。文獻[5]結合X射線測厚儀測量原理，著重分析了射線穩(wěn)定性、合金材質、帶鋼位置以及環(huán)境對測量精度的影響并給出了理論補償方法。文獻[6]和文獻[7]從具體的應用出發(fā)，介紹了工業(yè)應用現(xiàn)場中，對X射線測厚儀的精度補償方法。而針對具體使用過程中的精度維護過程，文獻[8]介紹了了實際生產(chǎn)過程中對X射線測厚儀的維護方法，以保證X射線測厚儀滿足精度要求。

從上述文獻中可以得知，在實際應用過程中，要完成對X射線測厚儀的理論測量曲線繪制和測量精度補償是一項較為復雜的過程。另外，由于使用場景的不同以及環(huán)境的影響，一次的理論計算并不能夠長期使用。因此需要尋求一種快捷方便而又具備一定程度精度保證的測量曲線繪制與精度補償方法。

為解決上述問題，本文結合漣鋼RM215X射線測厚儀在冷軋過程的具體應用實例，提出一種基于數(shù)據(jù)的X射線測厚儀測量曲線繪制與精度補償方法，并依據(jù)現(xiàn)場的實際檢測數(shù)據(jù)給出具體的應用過程。

1 基于數(shù)據(jù)的X射線測厚儀檢測模型

X射線測厚儀的基本原理為：當X射線穿過物體時光譜區(qū)輻射將被改物體的材料削弱。如圖1所示，X射線測厚儀發(fā)射輻射強度為I0，通過厚度為x的材料后的輻射強度為I，對于特定能量或波長，發(fā)射的輻射量和通過材料的的輻射量用下式給出：

式中：I0為X射線的入射輻射強度；I為通過厚度為x的樣片后的輻射強度；u為樣片材料的線性吸收系數(shù)；p為樣片材料密度；x為樣片厚度。其中，吸收系數(shù)u是材料對X射線進行衰減的能力的量度，是一個復雜的函數(shù)，并且隨著能量的衰減而非線性地減小。它取決于兩個因素：材料的化學組成以及X射線的能量。

圖1 X射線測厚儀測量社示意圖

根據(jù)測厚儀的測量原理，測厚儀傳感器的電壓信號強度V與X射線接收強度I有如下式所示的關系。

將式（2）代入式中可以得到樣片厚度與測厚儀傳感器電壓信號之間的關系式。

在V0已知下，可得到如式（4）所示的一般表達式

由于u是與發(fā)射管電壓和材質原子序數(shù)相關的復雜的非線性表達式，所以a也是一個復雜的非線性表達式。

根據(jù)文獻[3]的理論分析，在冷軋現(xiàn)場，安裝在固定位置，選取合適的入射強度，配置穩(wěn)定的電源之后，X射線測厚儀的檢測精度還受統(tǒng)計誤差、材質吸收系數(shù)u以及環(huán)境的溫度、濕度、帶材表面附著物的影響。其中，統(tǒng)計誤差與材質吸收系數(shù)、板材厚度具有復雜的非線性函數(shù)關系，材質吸收系數(shù)u反比于管電壓V的3次方，正比于被測物原子序數(shù)的3次方，而環(huán)境對檢測精度的影響，理論分析則更為困難。

依據(jù)上述分析，要完成理論上的X射線測厚儀理論測量曲線繪制和精度修正是極為困難的。實際上，即使完成理論上的計算與修正，在實際應用時還需要根據(jù)具體檢測情況進行精度的補償。在這種情況下，考慮到在設備完成安裝后，上述對X射線測厚儀精度影響的各個因素在相當長的一個時期內(nèi)都不會發(fā)生大的變化，因此可以使用基于標準檢測數(shù)據(jù)的數(shù)據(jù)建模方法來實現(xiàn)X射線測厚儀的測量曲線繪制工作。

上述分析表明在公式表達式a為一個三階以上非線性表達式，而常量b為與具體安裝環(huán)境相關的補償閾值。經(jīng)過數(shù)據(jù)與機理的綜合分析，這些復雜非線性關系可以采用5階多項式來描述：

其中，z=-1（V）.

通過機理分析，得到上述X射線檢測儀檢測曲線模型的模型結構，但是其中的相關參數(shù)還需要通過現(xiàn)場的標準檢測數(shù)據(jù)進行辨識，通過數(shù)據(jù)的訓練得到完整的檢測曲線模型。

2 基于最小二乘的多項式系數(shù)辨識方法

對上一節(jié)建立的X射線測厚儀檢測模型結構分析表明，模型為典型的多項式系數(shù)辨識問題。針對多項式系數(shù)辨識問題，最小二乘多項式擬合方法計算相對較為簡單，并且能夠得到較為精確的結果，因此本節(jié)采用最小二乘多項式擬合方法辨識檢測模型中的系數(shù)。

2.1 最小二乘法基本原理

從整體上考慮近似函數(shù)p（x）同所給數(shù)據(jù)點（xi，y）i（i=0，1，…，m）誤差ri=p（x）i-yi（i=0，1，…，m）的大小，常用的方法有以下三種：一是誤差ri=p（x）i-yi（i=0，1，… ，m）絕 對 值 的 最 大 值，即誤差向量 r=（r1，r2，…，rm）T的∞-范數(shù)；二是誤差絕對值的和，即誤差向量r

數(shù)據(jù)擬合的具體做法是：對給定數(shù)據(jù)（xi，yi）（i=0，1，…，m），在取定的函數(shù)類 Φ 中，求 p（x）∈Φ，使誤差ri=p（xi）-yi（i=0，1，…，m）的平方和最小，即

從幾何意義上講，就是尋求與給定點（xi，yi）（i=0，1，…，m）的距離平方和為最小的曲線 y=p（x），如圖2所示。函數(shù)p（x）稱為擬合函數(shù)或最小二乘解，求擬合函數(shù)p（x）的方法稱為曲線擬合的最小二乘法。

在曲線擬合中，函數(shù)類可有不同的選取方法。

圖2 曲線擬合示意圖

2.2 最小二乘多項式擬合

假設給定數(shù)據(jù)點（xi，yi）（i=0，1，…，m），Φ為所有次數(shù)不超過n（n≤m）的多項式構成的函數(shù)類，現(xiàn)求一使得

當擬合函數(shù)為多項式時，稱為多項式擬合，滿足式（7）的pn（x）稱為最小二乘擬合多項式。特別地，當n=1時，稱為線性擬合或直線擬合。

顯然式為a0，a1，…an的多元函數(shù)，因此上述問題即為求I=I（a0，a1，…an）的極值問題。由多元函數(shù)求極值的必要條件，得式（9）為關于 a0，a1，…an的線性方程組，用矩陣表示為：

該式（9）稱為正規(guī)方程組或法方程組。

可以證明，方程組的系數(shù)矩陣是一個對稱正定矩陣，故存在唯一解。從式（10）中解出ak（k=0，1，…，n），從而可得多項式可以證明，式（11）中的p（nx）滿足式（7），即p（nx）為所求的擬合多項式。我們把稱為最小二乘擬合多項式p（nx）的平方誤差，記作

由式（8）可得：

在實際應用中，n＜m或n≤m；當n＝m時所得的擬合多項式就是拉格朗日或牛頓插值多項式。通過上述多項式擬合方法，在給定標準測量數(shù)據(jù)基礎上，即能夠計算得到確定的X射線檢測曲線。

2.3 最小二乘擬合多項式的存在唯一性

為論證最小二多項式擬合的解具有唯一性，給出如下兩個定理。

定理1設節(jié)點x0，x1，…，xn互異，則法方程組（10）的解存在唯一。

證：由克萊姆法則，只需證明方程組的系數(shù)矩陣非奇異即可。

用反證法，設方程組（10）的系數(shù)矩陣奇異，則其所對應的齊次方程組為：

上式有非零解。式（14）可寫為

將式（15）中第 j個方程乘以 aj（j=0，1，…，n），然后將新得到的n+1個方程左右兩端分別相加，得：

式中，pn（x）是次數(shù)不超過n的多項式，它有m+1>n個相異零點，由代數(shù)基本定理，必須有a0，a1，…an＝0，與齊次方程組有非零解的假設矛盾。因此正規(guī)方程組（14）必有唯一解。

定理 2：設 a0，a1，…an是正規(guī)方程組（14）的解，則是滿足式（7）的最小二乘擬合多項式。

證：只需證明，對任意一組數(shù)b0，b1，…bn組成的多項式恒有

因為 ak（k=0，1，…，n）是正規(guī)方程組（14）的解，所以滿足式（8），因此有

故pn（x）為最小二乘擬合多項式。

上述證明表明，在給定足夠的不具備線性相關辨識數(shù)據(jù)情況下，總能得到唯一的檢測模型唯一解。

3 基于數(shù)據(jù)的X射線測厚儀檢測模型應用

在漣鋼冷軋生產(chǎn)現(xiàn)場按照要求安裝完RM215型X射線測厚儀后，采用本文的方法繪制該測厚儀的測量曲線，具體的過程如下。

首先，根據(jù)九個材料樣片標準件插入到X射線源發(fā)射的光束中進行測量，得到表1所示的標準測量數(shù)據(jù)，后續(xù)將利用該標準測量數(shù)據(jù)對測量模型進行參數(shù)辨識。

表1 測厚儀標準樣片測量數(shù)據(jù)

根據(jù)測厚儀的測量原理，X射線的實際入射輻射量和穿過樣板后的實際接收輻射量加上X射線在空氣中的衰減輻射量等于X射線的入射輻射量與接收輻射量的測量值，所以把表1轉化成表2所示的X射線實際入射與接收輻射量。

表2 X射線實際入射與輻射強度

根據(jù)式把表2數(shù)據(jù)轉化成表3所示的數(shù)據(jù)。

表3 樣片標準厚度與射線強度比值對照

利用表3中的9組值辨識測試模型中的系數(shù)，可以得到如下式方程組：

對式中各項相加得到下式：

式（24）~（29）是由公式（23）線性變換得到的，所以公式（24）~（29）與公式（23）是線性相關的，由線性相關性定理知公式（24）~（29）的解也是公式（23）的解。

將表4中的數(shù)據(jù)代入公式~列表如下：

表4 方程求解值表

解得多項式：

式（31）即為擬合得到的測量曲線，調整完擬合曲線，生成擬合修正值，寫入對應程序參數(shù)當中。這樣測厚儀的測量曲線基本完成?，F(xiàn)場的應用表明，采用該方法得到的X射線測厚儀測量曲線能夠滿足現(xiàn)場的應用需求，生產(chǎn)的冷軋板材厚度控制精度高且穩(wěn)定。

4 結束語

本文在分析X射線測厚儀檢測原理基礎上，指出實際生產(chǎn)現(xiàn)場影響測厚儀檢測精度的影響因素。根據(jù)分析的結論，給出一種通過標準檢測數(shù)據(jù)建立測厚儀檢測模型的數(shù)據(jù)驅動建模方法。通過建立的模型能夠完成測量曲線的繪制與精度補償。最后，通過漣鋼冷軋廠的實際應用實例表明該方法具有可操作性，且得到的測量曲線能夠滿足現(xiàn)場的精度要求，具備實用價值。

[1]趙春紅，趙檢羅.X射線測厚儀的應用與改進[J].漣鋼科技與管理，2008 ，（4）：27-29.

[2]馬竹梧.穿透式X射線測厚儀及其在冶金工業(yè)中的應用[J].計量裝置及應用，2010，20（6）：25-25.

[3]馬竹梧.X射線測厚影響因素分析,技術進展及其在冶金工業(yè)中的應用（上）[J].冶金自動化，2011，35（1）:1-5.

[4]馬竹梧.X射線測厚影響因素分析,技術進展及其在冶金工業(yè)中的應用（下）[J].冶金自動化，2011，35（2）：1-3.

[5]魏 鵬，方偉嶄.X射線測厚儀測量精度影響因素及補償措施[J].自動化儀表，2011，32（1）：17-18.

[6]陳德斌，胥福順，欒爭位.影響X射線測厚儀測量精度的因素及維護方法[J].輕合金加工技術，2009，37（4）：38-40.

[7]張曉春.X射線測厚儀在熱軋板帶中的應用[J].自動化儀表,2013，35（2）：88-91.

[8]陳德斌，胥福順，欒爭位.影響X射線測厚儀測量精度的因素及維護方法[J].輕合金加工技術，2009，37（4）：38-40.