哈薩克斯坦進(jìn)口液化石油氣短重原因分析

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

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

哈薩克斯坦進(jìn)口液化石油氣短重原因分析

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

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

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

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

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

1 實(shí)驗(yàn)方案設(shè)計(jì)

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

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

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

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

表2 車輛行進(jìn)順序及操作

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

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

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

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

2.1實(shí)驗(yàn)結(jié)果

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

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

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

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

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

表3 實(shí)驗(yàn)車輛皮重、毛重和凈重的平均值

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

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

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

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

表5 方差齊性檢驗(yàn)

表6 單因子方差分析表

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

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

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

從表8數(shù)據(jù)可以看出'獨(dú)山子和阿拉山口過衡數(shù)據(jù)的差值很小'兩地間皮重、毛重和凈重的平均差異分別為0.008 3、-0.049 0和-0.049 8t'換算為短溢率則分別為0.31‰、0.95‰和2.03‰'均小于±3‰的行業(yè)約定短溢率。獨(dú)山子和阿拉山口過衡數(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ī)律。而獨(dú)山子和多斯特克過衡結(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 實(shí)驗(yàn)結(jié)論及建議

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

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

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

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

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

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

[2］SN/T 0188.3-2010進(jìn)出口商品衡器鑒重規(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)計(jì)處理和解釋正態(tài)樣本離群值的判斷和處理[S］.

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

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

[9］尹海潔'劉耳.社會(huì)統(tǒng)計(jì)軟件SPSS for Windows簡明教程[M］.北京:社會(huì)科學(xué)文獻(xiàn)出版社'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)主要從事石油及天然氣計(jì)量與檢測(cè)的相關(guān)研究工作。