馬應(yīng)奎，張錫偉

中國檢驗認(rèn)證集團新疆有限公司(新疆烏魯木齊830011)

哈薩克斯坦進口液化石油氣短重原因分析

馬應(yīng)奎，張錫偉

中國檢驗認(rèn)證集團新疆有限公司(新疆烏魯木齊830011)

為了尋找哈薩克斯坦進口液化石油氣(LPG)短重率較高的原因，在資料審核無果的情況下，對58車次LPG分別在裝車站多斯特克、進境站阿拉山口和到貨站獨山子進行皮重和毛重的衡重，對各站間衡重所得結(jié)果的差異值進行對比分析。采用數(shù)理統(tǒng)計假設(shè)檢驗中方差分析的方法對皮重、毛重和凈重共9組差異值進行差異顯著性判定，得到皮重、毛重和凈重3組衡重數(shù)據(jù)間均存在顯著性差異。對各組數(shù)據(jù)進一步進行多重比較，得出了裝車站多斯特克的過衡結(jié)果與進境站阿拉山口和到貨站獨山子的過衡結(jié)果存在顯著性差異的結(jié)論，認(rèn)為導(dǎo)致進口LPG持續(xù)短重的原因應(yīng)是裝車站多斯特克的衡器存在系統(tǒng)偏差。對多斯特克的電子汽車衡進行重新檢定后，短重率大幅下降至貿(mào)易各方可以接受的3‰以內(nèi)，解決了持續(xù)困擾各方的短重問題。

液化石油氣;電子汽車衡;短重原因;方差分析

某公司自2013年10月開始采用汽運罐車進口液化石油氣(LPG)'裝車站為哈薩克斯坦共和國多斯特克'經(jīng)中國新疆阿拉山口口岸進境'到達(dá)新疆克拉瑪依市獨山子區(qū)。運距全程約280km'多斯特克～阿拉山口10km'阿拉山口～獨山子270km。質(zhì)量鑒定方式為電子汽車衡'在裝車站多斯特克使用發(fā)貨方的汽車衡衡重'在到貨站獨山子使用經(jīng)過商檢部門確認(rèn)的汽車衡進行衡重，LPG進口商委托的獨立第三方檢驗機構(gòu)在發(fā)貨站和到貨站分別進行過衡過程的見證。自業(yè)務(wù)開展以來'嚴(yán)重的短重問題一直困擾各相關(guān)方'短重率持續(xù)高達(dá)10‰左右[1-2］。到貨站查驗鉛封完好。多斯特克發(fā)貨方和獨山子收貨方均堅稱所使用的衡器經(jīng)過法定部門檢定合格'且在有效期內(nèi)'雙方衡重結(jié)果沒有問題。雙方提供的衡器檢定證書也驗證了雙方所使用衡器確實均為檢定有效。問題陷入僵局'進口商和收貨方的利益持續(xù)遭受損失'并導(dǎo)致貿(mào)易雙方產(chǎn)生對對方的不信任'影響了貿(mào)易的進一步擴大。為了尋找短重原因'由進口商、獨立檢驗方、承運方、收貨方四方經(jīng)過共同討論'制定了進口LPG短重原因查找實驗方案'并將方案提交發(fā)貨方'得到了對方的認(rèn)可。

1 實驗方案設(shè)計

1)選取一定數(shù)量的車輛在裝貨地多斯特克、入境口岸阿拉山口、到貨地獨山子三地分別對車輛進行皮重和毛重過衡'然后對3組過衡數(shù)據(jù)進行分析對比'尋找各組數(shù)據(jù)之間的差別。為了避免由于車輛燃油變化對實驗的干擾'在每一地過衡前均將車輛油箱加滿[3］。

2)衡器采用:獨山子、阿拉山口、多斯特克三地均為數(shù)字型汽車衡'并經(jīng)所在地計量部門檢定合格'最大秤量值分別為100、120和80t'分度值均為20kg。所秤量范圍均在衡器的正常秤量范圍內(nèi)。三地所用衡器的詳細(xì)情況見表1。

3)過衡順序:以獨山子作為車輛出發(fā)和最終到達(dá)點'在各點的操作見表2。

表1 實驗所用電子汽車衡詳細(xì)數(shù)據(jù)

表2 車輛行進順序及操作

4)抽樣量:計劃抽取月計劃進口量12 000t(約500車)的10%'50車。

5)實驗過程見證:獨山子由收貨方、承運方、進口商、獨立檢驗方共同派員見證;阿拉山口由承運方、進口方委托貨代公司、獨立檢驗方共同派員見證;多斯特克由承運方負(fù)責(zé)見證。運輸車輛在每次稱量前'按照方案要求對油箱和鉛封狀態(tài)進行查驗。

6)最終對比數(shù)據(jù)由獨立檢驗方采用數(shù)理統(tǒng)計工具進行匯兌分析[4-5］。

2 實驗結(jié)果及數(shù)據(jù)分析

2.1實驗結(jié)果

經(jīng)過協(xié)調(diào)國內(nèi)外各方資源'于2014年6月分4批次安排了對比實驗'實際參加對比實驗車輛58車次'實驗車輛分別在三地過衡所得的皮重、毛重及計算凈重的平均值見表3。

2.2 數(shù)據(jù)處理及分析

1)分別計算出獨山子-多斯特克、阿拉山口-多斯特克和獨山子-阿拉山口3組數(shù)據(jù)間的皮重差、毛重差和凈重差。

2)對計算出的各組差值數(shù)據(jù)利用格拉布斯準(zhǔn)則(Grubbs)進行離群值的判斷'排除離群值;取置信度α=0.05'樣本量n=58'查表可知T0(58，0.05)=3.02。有若計算出某個值Tx＞3.02，則該值為離群值'不參加后續(xù)數(shù)據(jù)處理。經(jīng)計算共排除離群值7個'其中獨山子-多斯特克毛重差離群值1個'獨山子-多斯特克凈重差離群值5個'獨山子-阿拉山口凈重差離群值1個。

3)使用IBM SPSS Statistics專業(yè)統(tǒng)計軟件進行處理'計算出各組差值數(shù)據(jù)的均值、標(biāo)準(zhǔn)差、最大值、最小值、置信區(qū)間等'進行方差分析'尋找3組數(shù)據(jù)間是否存在顯著性差異和組間差異[6］。三地間衡重結(jié)果差的均值、標(biāo)準(zhǔn)差、95%置信區(qū)間和極值見表4'方差齊性檢驗結(jié)果見表5'單因子方差分析結(jié)果見表6'各組數(shù)據(jù)間的兩兩對比結(jié)果見表7。

表3 實驗車輛皮重、毛重和凈重的平均值

從表4的統(tǒng)計量可以看出'獨山子-阿拉山口間的衡重差皮重、毛重和凈重均值都遠(yuǎn)小于獨山子-多斯特克和阿拉山口-多斯特克間的差。

從表5的檢驗結(jié)果可知'3組數(shù)據(jù)的差異的顯著性水均小于設(shè)定顯著性水平因子0.05'該方差為非齊性。

從表6計算結(jié)果可知'皮重、毛重和凈重3組差異數(shù)據(jù)的組間均存在顯著性差異'想知道到底是哪一組數(shù)據(jù)與其他兩組數(shù)據(jù)存在顯著性差異還是3組數(shù)據(jù)間均存在顯著性差異需要進行各組間數(shù)據(jù)的多重比較。

表4 差值描述性統(tǒng)計量

表5 方差齊性檢驗

表6 單因子方差分析表

表7顯示了進一步多重比較的結(jié)果'因為方差為非齊性'困此以Tamhane模式計算結(jié)果為準(zhǔn)。計算結(jié)果顯示:①皮重衡重結(jié)果間的3組差異中'獨山子-多斯特克和阿拉山口-多斯特克與獨山子-阿拉山口差值間均存在統(tǒng)計學(xué)意義上的顯著性差異;阿拉山口-多斯特克與獨山子-多斯特克間差值不存在統(tǒng)計學(xué)意義上的顯著性差異;②毛重衡重結(jié)果間的3組差異中'獨山子-阿拉山口、獨山子-多斯特克、阿拉山口-多斯特克間均存在統(tǒng)計學(xué)意義上的顯著性差異;③凈重衡重結(jié)果間的3組差異中'獨山子-多斯特克和阿拉山口-多斯特克與獨山子-阿拉山口差值間均存在統(tǒng)計學(xué)意義上的顯著性差異;阿拉山口-多斯特克與獨山子-多斯特克間差值不存在統(tǒng)計學(xué)意義上的顯著性差異;④以數(shù)據(jù)較接近的獨山子-阿拉山口間的差異值為控制組'則表6中Dunnett t(雙側(cè))模式計算結(jié)果顯示'皮重、毛重和凈重3組數(shù)據(jù)中'獨山子-多斯特克和阿拉山口-多斯特克差值均與獨山子-阿拉山口對照組存在顯著性差異[7］。

表7 多重比較結(jié)果

4)為了將差異換算為貿(mào)易習(xí)慣性的千分率的表達(dá)方式'分別取表3計算的58車實驗車輛的皮重、毛重和凈重的平均數(shù)作為計算千分率的基數(shù)[8］。58個實驗車輛在三地過衡的皮重、毛重和凈重的平均值見表7。將表7計算出的均值、95%置信區(qū)間的值分別轉(zhuǎn)換為千分率后所得結(jié)果見表8。

從表8數(shù)據(jù)可以看出'獨山子和阿拉山口過衡數(shù)據(jù)的差值很小'兩地間皮重、毛重和凈重的平均差異分別為0.008 3、-0.049 0和-0.049 8t'換算為短溢率則分別為0.31‰、0.95‰和2.03‰'均小于±3‰的行業(yè)約定短溢率。獨山子和阿拉山口過衡數(shù)據(jù)的差值的95%置信區(qū)間'皮重為-0.98‰～1.59‰'毛重為-1.33‰～-0.57‰'凈重為-2.95‰～-0.70‰'均小于±3‰。阿拉山口和多斯特克間的過衡結(jié)果的差值也較大'兩地間皮重、毛重和凈重的平均差異分別為-0.203 1、-0.456 2和-0.253 1t'換算為短溢率則分別為-7.53‰、-8.87‰和-10.33‰'遠(yuǎn)高于±3‰的行業(yè)約定短溢率。差值的95%置信區(qū)間呈現(xiàn)與均值同樣的規(guī)律。而獨山子和多斯特克過衡結(jié)果的差值明顯要大'兩地間皮重、毛重和凈重的平均差異分別達(dá)到了-0.194 8、-0.511 6和-0.286 0t'換算為短溢率則分別為-7.22‰、-9.94‰和-11.68‰'遠(yuǎn)高于±3‰的行業(yè)約定短溢率。

表8 差值的均值、95%置信區(qū)間的轉(zhuǎn)換為千分率的結(jié)果

3 實驗結(jié)論及建議

1)目前世界各國計量體系的量值溯源系統(tǒng)'任何與重量(質(zhì)量)有關(guān)的標(biāo)準(zhǔn)器都要最終溯源到國際計量局(BIPM)的國際基準(zhǔn)千克原器'因而由于兩國計量體系的差異導(dǎo)致衡重結(jié)果產(chǎn)生較大差異的原因可以排除。

2)通過統(tǒng)計學(xué)的計算和現(xiàn)實差異率的計算'可以得出結(jié)論:獨山子-多斯特克差值、阿拉山口-多斯特克差值與獨山子-阿拉山口間差值在皮重、毛重和凈重上均存在顯著性差異'可以得出多斯特克過衡結(jié)果明顯異常于獨山子和阿拉山口的結(jié)果。

3)雖然毛重計算結(jié)果中獨山子-多斯特克差值與阿拉山口-多斯特克差值間也存在統(tǒng)計學(xué)意義上的顯著性差異'似乎獨山子與阿拉山口所過衡的毛重間也存在較大差異'但實際上獨山子與阿拉山口間毛重間差異平均值為50kg'差異率僅0.95‰'遠(yuǎn)小于行業(yè)慣例可以接受的±3‰的差異率'因而這種差異只存在統(tǒng)計學(xué)意義'現(xiàn)實中是完全可以接受的。

4)通過皮重和毛重的差異對比可以發(fā)現(xiàn)'貨物的重量越大'多斯特克的過衡結(jié)果與阿拉山口、獨山子的過衡結(jié)果差異也越大(凈重值為計算值'不是直接秤量值)。阿拉山口-多斯特克過衡的差異均值皮重差為-0.203 1t'-7.53‰'毛重差異則擴大至-0.456 2t'-8.87‰。獨山子-多斯特克過衡的差異均值皮重差-0.194 8t'-7.22‰'毛重差異則上升至-0.511 6t'-9.94‰。由此推斷多斯特克的衡器可能存在秤量值偏大的系統(tǒng)偏差'且該偏差隨著所秤量貨物重量的增加而增大。

5)建議協(xié)調(diào)有關(guān)方面'對裝車站多斯特克的電子汽車衡重新進行科學(xué)嚴(yán)謹(jǐn)?shù)臋z定'修正存在的偏差。

[1］SN/T 0188.2-2010進出口商品衡器鑒重規(guī)程第2部分:衡器鑒重通則[S］.

[2］SN/T 0188.3-2010進出口商品衡器鑒重規(guī)程第3部分:汽車衡器鑒重[S］.

[3］ГОСT 8.453-82靜態(tài)秤重用衡器檢定方法和器具[S］.

[4］ГОСT 29329-92靜態(tài)秤重用衡器總要求[S］.

[5］JJG 539-1997數(shù)字指示秤[S］.

[6］GB/T 4883-2008數(shù)據(jù)的統(tǒng)計處理和解釋正態(tài)樣本離群值的判斷和處理[S］.

[7］GB/T 4086.4-1983統(tǒng)計分布數(shù)值表F分布[S］.

[8］GB/T 10092-2009數(shù)據(jù)的統(tǒng)計處理和解釋測試結(jié)果的多重比較[S］.

[9］尹海潔'劉耳.社會統(tǒng)計軟件SPSS for Windows簡明教程[M］.北京:社會科學(xué)文獻出版社'2003.

The weight of the liquefied petroleum gas(LPG)imported from Kazakhstan is lacking for a long time.For finding out the causes，the tare and gross weight of 58 cars weighed at Duositeke station，Alashankou station and Dushanzi station separately.The difference in 9 group data of tare，gross and net weight of 58 cars is analyzed using the variance analysis of mathematical statistics，and the results show that the difference of 9 group data is all significant.By comparing the data，it is held that the less weight of LPG imported from Kazakhstan is caused by the systemic error of the electronic car weigher at Duositeke station.By the correction of the weighter，the less weight rate is decreased to within 3‰，which is acceptable to two parties.

liquefied petroleum gas;electronic car weigher;cause of less weight;variance analysis

左學(xué)敏

2014-11-17

馬應(yīng)奎(1976-)'男'工程師'現(xiàn)主要從事石油及天然氣計量與檢測的相關(guān)研究工作。