马应奎,张锡伟
中国检验认证集团新疆有限公司(新疆乌鲁木齐830011)
哈萨克斯坦进口液化石油气短重原因分析
马应奎,张锡伟
中国检验认证集团新疆有限公司(新疆乌鲁木齐830011)
为了寻找哈萨克斯坦进口液化石油气(LPG)短重率较高的原因,在资料审核无果的情况下,对58车次LPG分别在装车站多斯特克、进境站阿拉山口和到货站独山子进行皮重和毛重的衡重,对各站间衡重所得结果的差异值进行对比分析。采用数理统计假设检验中方差分析的方法对皮重、毛重和净重共9组差异值进行差异显著性判定,得到皮重、毛重和净重3组衡重数据间均存在显著性差异。对各组数据进一步进行多重比较,得出了装车站多斯特克的过衡结果与进境站阿拉山口和到货站独山子的过衡结果存在显著性差异的结论,认为导致进口LPG持续短重的原因应是装车站多斯特克的衡器存在系统偏差。对多斯特克的电子汽车衡进行重新检定后,短重率大幅下降至贸易各方可以接受的3‰以内,解决了持续困扰各方的短重问题。
液化石油气;电子汽车衡;短重原因;方差分析
某公司自2013年10月开始采用汽运罐车进口液化石油气(LPG)'装车站为哈萨克斯坦共和国多斯特克'经中国新疆阿拉山口口岸进境'到达新疆克拉玛依市独山子区。运距全程约280km'多斯特克~阿拉山口10km'阿拉山口~独山子270km。质量鉴定方式为电子汽车衡'在装车站多斯特克使用发货方的汽车衡衡重'在到货站独山子使用经过商检部门确认的汽车衡进行衡重,LPG进口商委托的独立第三方检验机构在发货站和到货站分别进行过衡过程的见证。自业务开展以来'严重的短重问题一直困扰各相关方'短重率持续高达10‰左右[1-2]。到货站查验铅封完好。多斯特克发货方和独山子收货方均坚称所使用的衡器经过法定部门检定合格'且在有效期内'双方衡重结果没有问题。双方提供的衡器检定证书也验证了双方所使用衡器确实均为检定有效。问题陷入僵局'进口商和收货方的利益持续遭受损失'并导致贸易双方产生对对方的不信任'影响了贸易的进一步扩大。为了寻找短重原因'由进口商、独立检验方、承运方、收货方四方经过共同讨论'制定了进口LPG短重原因查找实验方案'并将方案提交发货方'得到了对方的认可。
1)选取一定数量的车辆在装货地多斯特克、入境口岸阿拉山口、到货地独山子三地分别对车辆进行皮重和毛重过衡'然后对3组过衡数据进行分析对比'寻找各组数据之间的差别。为了避免由于车辆燃油变化对实验的干扰'在每一地过衡前均将车辆油箱加满[3]。
2)衡器采用:独山子、阿拉山口、多斯特克三地均为数字型汽车衡'并经所在地计量部门检定合格'最大秤量值分别为100、120和80t'分度值均为20kg。所秤量范围均在衡器的正常秤量范围内。三地所用衡器的详细情况见表1。
3)过衡顺序:以独山子作为车辆出发和最终到达点'在各点的操作见表2。
表1 实验所用电子汽车衡详细数据
表2 车辆行进顺序及操作
4)抽样量:计划抽取月计划进口量12 000t(约500车)的10%'50车。
5)实验过程见证:独山子由收货方、承运方、进口商、独立检验方共同派员见证;阿拉山口由承运方、进口方委托货代公司、独立检验方共同派员见证;多斯特克由承运方负责见证。运输车辆在每次称量前'按照方案要求对油箱和铅封状态进行查验。
6)最终对比数据由独立检验方采用数理统计工具进行汇兑分析[4-5]。
2.1实验结果
经过协调国内外各方资源'于2014年6月分4批次安排了对比实验'实际参加对比实验车辆58车次'实验车辆分别在三地过衡所得的皮重、毛重及计算净重的平均值见表3。
2.2 数据处理及分析
1)分别计算出独山子-多斯特克、阿拉山口-多斯特克和独山子-阿拉山口3组数据间的皮重差、毛重差和净重差。
2)对计算出的各组差值数据利用格拉布斯准则(Grubbs)进行离群值的判断'排除离群值;取置信度α=0.05'样本量n=58'查表可知T0(58,0.05)=3.02。有若计算出某个值Tx>3.02,则该值为离群值'不参加后续数据处理。经计算共排除离群值7个'其中独山子-多斯特克毛重差离群值1个'独山子-多斯特克净重差离群值5个'独山子-阿拉山口净重差离群值1个。
3)使用IBM SPSS Statistics专业统计软件进行处理'计算出各组差值数据的均值、标准差、最大值、最小值、置信区间等'进行方差分析'寻找3组数据间是否存在显著性差异和组间差异[6]。三地间衡重结果差的均值、标准差、95%置信区间和极值见表4'方差齐性检验结果见表5'单因子方差分析结果见表6'各组数据间的两两对比结果见表7。
表3 实验车辆皮重、毛重和净重的平均值
从表4的统计量可以看出'独山子-阿拉山口间的衡重差皮重、毛重和净重均值都远小于独山子-多斯特克和阿拉山口-多斯特克间的差。
从表5的检验结果可知'3组数据的差异的显著性水均小于设定显著性水平因子0.05'该方差为非齐性。
从表6计算结果可知'皮重、毛重和净重3组差异数据的组间均存在显著性差异'想知道到底是哪一组数据与其他两组数据存在显著性差异还是3组数据间均存在显著性差异需要进行各组间数据的多重比较。
表4 差值描述性统计量
表5 方差齐性检验
表6 单因子方差分析表
表7显示了进一步多重比较的结果'因为方差为非齐性'困此以Tamhane模式计算结果为准。计算结果显示:①皮重衡重结果间的3组差异中'独山子-多斯特克和阿拉山口-多斯特克与独山子-阿拉山口差值间均存在统计学意义上的显著性差异;阿拉山口-多斯特克与独山子-多斯特克间差值不存在统计学意义上的显著性差异;②毛重衡重结果间的3组差异中'独山子-阿拉山口、独山子-多斯特克、阿拉山口-多斯特克间均存在统计学意义上的显著性差异;③净重衡重结果间的3组差异中'独山子-多斯特克和阿拉山口-多斯特克与独山子-阿拉山口差值间均存在统计学意义上的显著性差异;阿拉山口-多斯特克与独山子-多斯特克间差值不存在统计学意义上的显著性差异;④以数据较接近的独山子-阿拉山口间的差异值为控制组'则表6中Dunnett t(双侧)模式计算结果显示'皮重、毛重和净重3组数据中'独山子-多斯特克和阿拉山口-多斯特克差值均与独山子-阿拉山口对照组存在显著性差异[7]。
表7 多重比较结果
4)为了将差异换算为贸易习惯性的千分率的表达方式'分别取表3计算的58车实验车辆的皮重、毛重和净重的平均数作为计算千分率的基数[8]。58个实验车辆在三地过衡的皮重、毛重和净重的平均值见表7。将表7计算出的均值、95%置信区间的值分别转换为千分率后所得结果见表8。
从表8数据可以看出'独山子和阿拉山口过衡数据的差值很小'两地间皮重、毛重和净重的平均差异分别为0.008 3、-0.049 0和-0.049 8t'换算为短溢率则分别为0.31‰、0.95‰和2.03‰'均小于±3‰的行业约定短溢率。独山子和阿拉山口过衡数据的差值的95%置信区间'皮重为-0.98‰~1.59‰'毛重为-1.33‰~-0.57‰'净重为-2.95‰~-0.70‰'均小于±3‰。阿拉山口和多斯特克间的过衡结果的差值也较大'两地间皮重、毛重和净重的平均差异分别为-0.203 1、-0.456 2和-0.253 1t'换算为短溢率则分别为-7.53‰、-8.87‰和-10.33‰'远高于±3‰的行业约定短溢率。差值的95%置信区间呈现与均值同样的规律。而独山子和多斯特克过衡结果的差值明显要大'两地间皮重、毛重和净重的平均差异分别达到了-0.194 8、-0.511 6和-0.286 0t'换算为短溢率则分别为-7.22‰、-9.94‰和-11.68‰'远高于±3‰的行业约定短溢率。
表8 差值的均值、95%置信区间的转换为千分率的结果
1)目前世界各国计量体系的量值溯源系统'任何与重量(质量)有关的标准器都要最终溯源到国际计量局(BIPM)的国际基准千克原器'因而由于两国计量体系的差异导致衡重结果产生较大差异的原因可以排除。
2)通过统计学的计算和现实差异率的计算'可以得出结论:独山子-多斯特克差值、阿拉山口-多斯特克差值与独山子-阿拉山口间差值在皮重、毛重和净重上均存在显著性差异'可以得出多斯特克过衡结果明显异常于独山子和阿拉山口的结果。
3)虽然毛重计算结果中独山子-多斯特克差值与阿拉山口-多斯特克差值间也存在统计学意义上的显著性差异'似乎独山子与阿拉山口所过衡的毛重间也存在较大差异'但实际上独山子与阿拉山口间毛重间差异平均值为50kg'差异率仅0.95‰'远小于行业惯例可以接受的±3‰的差异率'因而这种差异只存在统计学意义'现实中是完全可以接受的。
4)通过皮重和毛重的差异对比可以发现'货物的重量越大'多斯特克的过衡结果与阿拉山口、独山子的过衡结果差异也越大(净重值为计算值'不是直接秤量值)。阿拉山口-多斯特克过衡的差异均值皮重差为-0.203 1t'-7.53‰'毛重差异则扩大至-0.456 2t'-8.87‰。独山子-多斯特克过衡的差异均值皮重差-0.194 8t'-7.22‰'毛重差异则上升至-0.511 6t'-9.94‰。由此推断多斯特克的衡器可能存在秤量值偏大的系统偏差'且该偏差随着所秤量货物重量的增加而增大。
5)建议协调有关方面'对装车站多斯特克的电子汽车衡重新进行科学严谨的检定'修正存在的偏差。
[1]SN/T 0188.2-2010进出口商品衡器鉴重规程第2部分:衡器鉴重通则[S].
[2]SN/T 0188.3-2010进出口商品衡器鉴重规程第3部分:汽车衡器鉴重[S].
[3]ГОСT 8.453-82静态秤重用衡器检定方法和器具[S].
[4]ГОСT 29329-92静态秤重用衡器总要求[S].
[5]JJG 539-1997数字指示秤[S].
[6]GB/T 4883-2008数据的统计处理和解释正态样本离群值的判断和处理[S].
[7]GB/T 4086.4-1983统计分布数值表F分布[S].
[8]GB/T 10092-2009数据的统计处理和解释测试结果的多重比较[S].
[9]尹海洁'刘耳.社会统计软件SPSS for Windows简明教程[M].北京:社会科学文献出版社'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
左学敏
2014-11-17
马应奎(1976-)'男'工程师'现主要从事石油及天然气计量与检测的相关研究工作。