基于回歸分析法的計量安全檢測設(shè)備運行質(zhì)量分析

牛曉明

(中國鐵路北京局集團有限公司 貨運部，北京 100860)

計量安全檢測設(shè)備包括計量衡器和安全檢測監(jiān)控設(shè)施2大類。軌道衡屬于計量衡器，是鐵路標準計量器具。超偏載檢測裝置屬于安全檢測監(jiān)控設(shè)施，是“鐵路貨運計量安全檢測監(jiān)控系統(tǒng)”的重要組成部分，是檢測貨車超載、偏載、偏重的主要裝置之一，是應用新技術(shù)保證運輸安全的重要手段，是行車事故分析的主要依據(jù)之一，在確保運輸安全上有著重要作用。軌道衡作為國家強制檢定的計量器具，其檢測數(shù)據(jù)具備法律效力，其運行穩(wěn)定性、數(shù)據(jù)準確性較超偏載檢測裝置要高，因而以軌道衡設(shè)備檢測數(shù)值為標準，分析超偏載檢測裝置的檢測狀態(tài)，做好貨運計量安全檢測監(jiān)控系統(tǒng)中超偏載檢測裝置和軌道衡檢測數(shù)據(jù)的對比，對設(shè)備運行質(zhì)量進行分析。

1 基于回歸分析法的計量安全檢測設(shè)備運行質(zhì)量分析

1.1 設(shè)備運行質(zhì)量指標分析

理想情況下對同一車輛進行計量安全檢測，采用超偏載檢測裝置檢測的數(shù)值和軌道衡檢測的數(shù)值應相同，但是由于受到設(shè)備安裝環(huán)境、列車運行狀態(tài)、車輛貨物狀態(tài)變化等影響，二者的檢測數(shù)值并不完全相同，軌道衡作為標準計量設(shè)備，其計量性能要高于超偏載檢測裝置。將軌道衡檢測數(shù)值“視為”無誤差，即將軌道衡的檢測數(shù)據(jù)按照車輛實際總重進行考慮[1]。

超偏載檢測裝置檢測車輛總重可以表示為

式中：y為超偏載檢測裝置檢測車輛總重，t；x為車輛實際總重，t；μ為超偏載檢測裝置檢測數(shù)據(jù)隨車輛總重變化而變化的比率，簡稱斜率，在理想超偏載檢測裝置下，μ= 1；z為無檢測車輛時超偏載檢測裝置顯示的重量，t，簡稱零點，在理想超偏載檢測裝置下，z= 0。一般情況下設(shè)備運行質(zhì)量完全體現(xiàn)在準確性和穩(wěn)定性2個方面，因而在將軌道衡檢測數(shù)據(jù)視為準確值的前提下，對超偏載檢測裝置運行質(zhì)量指標分析如下。①準確性。超偏載檢測裝置檢測數(shù)值的準確性體現(xiàn)在與實際重量相比檢測數(shù)值的誤差上，具體體現(xiàn)在零點z和斜率μ上。零點z：傳感器檢測值實際上是惠斯登電橋輸出端口的輸出電壓值，通過模數(shù)轉(zhuǎn)換器轉(zhuǎn)換成計算機可用的數(shù)字信息，即檢測的壓力值，對應的是貨車的輪重信息，而傳感器安裝時的狀態(tài)不同，零點電壓也不同，因而每次安裝傳感器后需重新設(shè)置檢測的零點數(shù)值，即在軌道上沒有壓力時惠斯登電橋產(chǎn)生的電壓值，相應z值是此時的顯示重量。斜率μ：傳感器檢測到軌道壓力后通過壓電轉(zhuǎn)換變?yōu)殡娦盘?，?jīng)模數(shù)轉(zhuǎn)換為相應數(shù)值，在計算機系統(tǒng)中轉(zhuǎn)換為相應稱重數(shù)值，需要設(shè)置一個固定參數(shù)，用以描述實際壓力變化和稱重數(shù)值變化之間的比例關(guān)系，這里描述為斜率。②穩(wěn)定性。按照超偏載檢測裝置的檢測數(shù)值和車輛實際總重的近似線性關(guān)系，二者數(shù)值表現(xiàn)出的線性關(guān)系越好，表明設(shè)備穩(wěn)定性越好，反之亦然。在數(shù)據(jù)分析中采用“判定系數(shù)”r2表示數(shù)據(jù)相關(guān)程度[2]，即樣本數(shù)據(jù)的線性擬合程度，該數(shù)值取值0-1之間，采樣數(shù)據(jù)的線性擬合程度越高該數(shù)值越近似為1。

1.2 基于回歸分析法的計量安全檢測設(shè)備運行質(zhì)量分析

回歸分析法是研究一個隨機變量與一個或幾個可控變量之間相關(guān)關(guān)系的統(tǒng)計方法。只有一個自變量的回歸分析法稱為一元回歸分析法。自變量與因變量具有線性關(guān)系的稱為一元線性回歸分析法。由于超偏載檢測獲得的數(shù)據(jù)和軌道衡檢測獲得的數(shù)據(jù)之間具有顯著的線性關(guān)系，因而可以使用一元線性回歸分析法對超偏載檢測數(shù)據(jù)和軌道衡檢測數(shù)據(jù)進行分析。具體步驟如下。

（1）提取樣本數(shù)據(jù)。從“貨運計量安全檢測監(jiān)控系統(tǒng)”的“數(shù)據(jù)比對”模塊中提取同一車輛經(jīng)過超偏載檢測裝置和軌道衡檢測的總重數(shù)據(jù)作為分析設(shè)備運行質(zhì)量的樣本數(shù)據(jù)。在樣本數(shù)據(jù)選取時需掌握以下原則：①同一臺超偏載檢測裝置對應的軌道衡不多于3臺，以免由于軌道衡運行質(zhì)量不穩(wěn)定導致的數(shù)據(jù)雜亂、相關(guān)度太低影響設(shè)備質(zhì)量的分析；②避免貨物狀態(tài)變化前后的數(shù)據(jù)混用，如夏季降雨造成的貨物增載前后分別通過不同的設(shè)備檢測的數(shù)據(jù)，不應作為設(shè)備運行質(zhì)量分析的樣本；③剝離個別離散度太大的異常數(shù)據(jù)，此類數(shù)據(jù)有可能是設(shè)備突發(fā)的異常值，如匹配錯誤、干擾等原因造成的異常，此類異常數(shù)據(jù)可另作分析[3]。④剝離異常數(shù)據(jù)的過程可能在數(shù)據(jù)分析的過程之中，由于在雜亂的數(shù)據(jù)中很難發(fā)現(xiàn)其中的個別異常數(shù)據(jù)，而在分析過程中能夠很容易發(fā)現(xiàn)與預期不相符的數(shù)據(jù)。

（2）描繪散點圖。設(shè)總量為n的樣本觀測值X和Y為一維變量xi∈X，yi∈Y，yi(i= 1，2，…，n)為x取固定值x=xi時y的觀測值，即超偏載檢測裝置的第i個檢測數(shù)值和對應的軌道衡第i個檢測數(shù)值，則(xi，yi)為樣本數(shù)據(jù)的第i對數(shù)據(jù)，可在平面直角坐標系中表示一個點。將樣本數(shù)據(jù)(x1，y1)，(x2，y2)，…，(xn，yn)在平面直角坐標系中描繪出相應的散點圖，觀察散點分布是否大致在一條直線附近，如果在一條直線附近，則可進行下一步分析。

（3）設(shè)備準確性分析。估計樣本數(shù)據(jù)的一元線性回歸方程y=μx+z的各個參數(shù)，進行準確性分析。在回歸分析中，斜率μ又稱為回歸系數(shù)，可以表示為

零點z又稱為回歸常數(shù)，可以表示為

（4）設(shè)備穩(wěn)定性分析。使用判定系數(shù)r對線性關(guān)系的顯著性進行檢驗[2]，進行穩(wěn)定性分析。判定系數(shù)r可以表示為

式中：為樣本數(shù)據(jù)X的平均值；為樣本數(shù)據(jù)Y的平均值；r為判定系數(shù)。在已知樣本數(shù)據(jù)確定為線性關(guān)系的前提下，可用|r|反映系統(tǒng)的穩(wěn)定性，即系統(tǒng)受到其他因素干擾后偏離線性關(guān)系的程度，為便于分析和計算，也可用r2代替|r|進行分析，0≤|r|≤1，|r|越接近1，則樣本X與Y的線性相關(guān)程度越密切，|r|越接近0，則樣本X與Y的線性相關(guān)程度越小。

2 案例分析

以2017年8月19日至20日某一臺超偏載檢測裝置和對應的一臺軌道衡的檢測數(shù)據(jù)為例，應用回歸分析法對計量安全檢測設(shè)備運行質(zhì)量進行分析。

2.1 樣本數(shù)據(jù)的提取和處理

按照樣本數(shù)據(jù)選取的原則從“貨運計量安全檢測監(jiān)控系統(tǒng)”中提取相應的數(shù)據(jù)并處理后，獲得超偏載檢測裝置和軌道衡的檢測數(shù)據(jù)樣本如表1所示。

表1 超偏載檢測裝置與軌道衡的檢測數(shù)據(jù)樣本Tab.1 Sample detectiondata of overloading and unbalanced loading detection device and rail weigher

表1直接列出超偏載檢測和軌道衡檢測的車輛總重數(shù)據(jù)。目前“貨運計量安全檢測監(jiān)控系統(tǒng)”中“數(shù)據(jù)比對”功能通過對該表中數(shù)據(jù)的簡單計算求取超偏載檢測裝置和軌道衡檢測數(shù)據(jù)的平均誤差，公式為

式中：w為超偏載檢測裝置相對于軌道衡的檢測數(shù)值的平均誤差。

平均值計算能夠在一定程度上反映超偏載檢測裝置的誤差?！皵?shù)據(jù)比對”功能使用了該方法，以天為單位對每臺超偏載檢測裝置檢測的數(shù)據(jù)對應軌道衡檢測數(shù)據(jù)進行比對，計算設(shè)備的誤差[4]。這種計算方式僅能反映超偏載檢測裝置通過此樣本計算出的誤差平均值，但是在計算過程丟失了檢測數(shù)據(jù)中反映設(shè)備穩(wěn)定性的信息[5]。因此，這種方法對數(shù)據(jù)的分析過于簡單，未能充分發(fā)掘數(shù)據(jù)中存在的更多的信息。

2.2 繪制散點圖

把樣本值 (x1，y1)，(x2，y2)，…，(xn，yn) 作為平面直角坐標系的n個點描繪出來，這n個點構(gòu)成的圖像叫做散點圖。超偏載檢測裝置的檢測數(shù)值散點圖如圖1所示。

圖1 超偏載檢測裝置的檢測數(shù)值散點圖Fig.1 Scatterplot of the detection values of the overloading and unbalanced loading detection device

如果散點圖中的n個點近似在一條直線的附近，則可以認為該樣本適合用一元線性回歸模型來分析。圖1所示超偏載檢測裝置的檢測數(shù)據(jù)和軌道衡的檢測數(shù)據(jù)具有明顯的線性特征，適合此模型進行分析。

2.3 線性擬合分析

按照公式 ⑴ 提出的描述超偏載檢測裝置運行質(zhì)量的數(shù)學模型，運用一元線性回歸分析的方法對數(shù)據(jù)進行分析。

在Excel中運用一元線性回歸的參數(shù)計算函數(shù)LINEST()計算樣本數(shù)據(jù)X和Y對應的擬合直線的相關(guān)參數(shù)數(shù)組，應用INDEX()函數(shù)將數(shù)組中的相應參數(shù)提取出來，回歸系數(shù)μ= INDEX (LINEST (Y，X，1，1)，1， 1)，回歸常數(shù)z= INDEX (LINEST (Y，X，1，1)，1，2)，判定系數(shù)r2= INDEX (LINEST (Y，X， 1，1)，3，1)。將計算得到的參數(shù)μ和z的值代入公式 ⑴，得到反映超偏載檢測裝置的檢測重量y和軌道衡檢測重量x線性關(guān)系的一元線性回歸方程。將方程表示的直線在散點圖所在的直角坐標系中描繪出來，稱為“趨勢線”。超偏載檢測裝置的檢測數(shù)值散點圖和趨勢線如圖2所示。

由于軌道衡也存在一定誤差，同時為了簡化分析過程，當零點值很小時，可以將零點值近似地認為0，在Excel中將趨勢線的“截距”設(shè)置為0，LINEST ()函數(shù)的第3個參數(shù)調(diào)整為0，則回歸系數(shù)μ= INDEX(LINEST (Y，X，0，1)，1，1)， 判 定 系數(shù)r2= INDEX(LINEST (Y，X，0，1)，3，1)，得到只含有2個參數(shù)的趨勢線。超偏載檢測裝置的檢測數(shù)值散點圖(截距設(shè)置為0)如圖3所示。

圖2 超偏載檢測裝置的檢測數(shù)值散點圖和趨勢線Fig.2 Scatterplot and trendline of the detection values of the overloading and unbalanced loading detection device

圖3 得出趨勢線對應的線性公式為y= 0.986x，μ= 0.986直觀地表示出超偏載檢測裝置的誤差情況，即誤差百分比為(1 - 0.986)×100% = 1.4%；判定系數(shù)r2= 0.989，即該線性關(guān)系的可信度為0.989，在軌道衡設(shè)備穩(wěn)定的前提下，可以作為超偏載檢測裝置穩(wěn)定性的判定指標。

圖3 超偏載檢測裝置的檢測數(shù)值散點圖(截距設(shè)置為0)Fig.3 Scatterplot of the detection values of the overloading and unbalanced loading detection device (Set the intercept to 0)

2.4 設(shè)備檢測準確性驗證

使用回歸分析法驗證超偏載檢測裝置準確性時，可設(shè)定幾個x值代入得到的線性公式進行驗證，如代入案例中得到的超偏載檢測裝置特征線性公式y(tǒng)= 0.924x+ 4.903。按空車25 t，即x= 25 t時，得到y(tǒng)= 28 t，超偏載檢測裝置檢測誤差為28 - 25 = 3 t；當重車x= 65 t時，得到y(tǒng)= 64.95 t，超偏載檢測裝置檢測數(shù)值基本準確；當重車x= 90 t時，得到y(tǒng)=88.04 t，超偏載檢測裝置檢測數(shù)值的誤差為88.04 -90 = -1.16 t；當重車x= 100 t時，得到y(tǒng)= 97.28 t，超偏載檢測裝置檢測數(shù)值的誤差將達到97.28 - 100 =-2.72 t。

從以上計算可知，當被檢測車輛總重數(shù)值偏離65 t越大，檢測的誤差越大，這也是造成超偏載檢測裝置檢測的“空車不空”“檢測數(shù)值偏低”問題的重要原因。

3 結(jié)束語

在貨運計量安全檢測設(shè)備的日常管理工作中使用回歸分析法對檢測數(shù)據(jù)進行分析，能夠準確判定超偏載檢測裝置或軌道衡的運行質(zhì)量問題，為設(shè)備維修單位提供準確的維修建議[6-7]，對提高設(shè)備運行質(zhì)量起到有效的數(shù)據(jù)支撐作用。然而，在“貨運計量安全檢測監(jiān)控系統(tǒng)”中尚未提供回歸分析法，如果在系統(tǒng)中編寫相應的分析模塊并修改設(shè)備運行質(zhì)量考核標準，將為設(shè)備管理人員和設(shè)備維修單位提供更準確的判定依據(jù)，有利于提高設(shè)備運行質(zhì)量；還可以充分利用系統(tǒng)中大量數(shù)據(jù)從多個分析維度實現(xiàn)更深入的分析[8]。基于回歸分析法的貨運計量安全檢測設(shè)備運行質(zhì)量分析對構(gòu)建貨運計量安全檢測設(shè)備運行質(zhì)量評價體系提供了方向。