2022加油站地埋罐容积标定加油机法

加油站地埋罐容积标定加油机法目 次前言 II引言 III范围 5规范引文件 5术语定义 5一般求 6标定备 7标定程 7标定果理判定 9标定积不定度 9附录A（料附）加机标记表 11附录B（料附）三非匀B样插法 13附录C（料附）加机标地罐的不定评定 17IPAGEPAGE5PAGEPAGE10加油站地埋罐容积标定加油机法范围本文件规定了在加油站的加油作业中，以加油机、液位计及温度计为主计量设备，并以它们同步测量的加油量、液位和温度为基础数据，在线标定地埋罐容积的原理、方法和步骤。本文件适用于地埋水平圆筒形罐的容积标定及容积表在使用中的期间核查。下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中，注日期的引用文件，仅该日期对应的版本适用于本文件；不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。GB/T8927石油和液体石油产品温度测量手工法GB/T13236石油和液体石油产品储罐液位手工测量设备GB/T13894石油和液体石油产品液位测量法手工法GB/T21117磁致伸缩液位计GB/T21451.11GB/T21451.44GB50156 JJG266 JJG443 JJF1059.1测量不确定度评定与表示下列术语和定义适用于本文件。3.1加油站fuellingstation包括加油加气合建站、加油加氢合建站和加油加气加氢合建站。3.2油罐容积volumeoftank地埋罐中液位高度对应的液体体积。3.3容积表volumetable地埋罐高度和容积对应关系的表格。3.4检尺点内竖直径insideverticaldiameteratdip-hatch从量油孔测量的地埋罐直圆筒部分的竖向内直径。3.5安全高度safetyheight地埋罐容积90%所对应的液位高度。3.6预设液位presetlevel在地埋罐安全高度至实际使用的最低液位之间，按照合理间隔预先设定的记录数据的液位。3.7记录液位recordlevel标定过程中实际记录数据的液位。3.8记录区间recordsegment两个相邻记录液位之间的区域。3.9记录范围recordrange最高记录液位至最低记录液位之间的区域。3.10罐枪关系tank-nozzlerelationship地埋罐和与之连接的油枪的对应关系。一个地埋罐可能对应单把油枪，也可能对应多把油枪，多于一把油枪的加油机应标注油枪序号，一把油枪只能连接一个地埋罐。3.11泵码pumpcode加油机记录的流量计读数累计值。3.12区间温度temperatureintheinterval同一记录区间内两个记录点所对应的地埋罐内油品温度的算术平均值。GB501561mm。50mm～100mm；10mm100mm8mm～10mm。100，且±0.5概要本方法的标定设备包括加油机、液位测量设备和地埋罐内油品温度测量设备。加油机JJG443JJG443液位测量设备可采用液位手工测量设备或地埋罐自动液位计。液位手工测量设备应符合GB/T13236，并应在检定或校准的有效期内使用。GB/T21451.1油温测量设备可采用便携式电子温度计（PET）或地埋罐自动温度计（ATT）。PETGB/T8927，应在检定或校准的有效期内使用。ATTGB/T21451.4进行选择、安装、校准和检验，其固有误差不超过±0.25℃，采用手工法的检验误差不超过±0.5℃。）GB/T21117𝒊𝑯1150mm10个。𝒊5JJG443𝜹𝒔,𝒌和重复性误差𝑬𝒏𝒔,𝒌，应符合第5章的要求，并进行数据记录。注：𝑘代表连接至待标定地埋罐对应加油机油枪的序号；符号s是开始标定前。ALGGB/T13894）ℎ𝑖，其中33𝐻𝑖；ATTPET（GB/T8927）1100mm𝑡𝑖；𝑣𝑘,𝑖。表1不同液位高度下地埋罐内油品温度的测量点数要求液位高度H/mm油品温度测量点数个H＜6001600≤H＜110021100≤H＜160031600≤H＜21004H≥21005（10m内6.1第n）JJG443𝜹𝒆,𝒌和𝑬𝒏𝒆,𝒌5注：𝑘代表连接至待标定地埋罐对应加油机油枪的序号；符号e是结束标定后。A中表A.1、表A.2。𝒋j𝑽（1）（3）𝒋=𝑣𝑘,𝑖−𝑣𝑘,(𝑖−1)…………（1）𝑘,𝑗𝑉′=𝑉𝑘,𝑗×(1−𝛿𝑘𝑘,𝑗

)………（2）𝑉′=∑𝑁

……………（3）式中：

𝑗

𝑘,𝑗𝑣𝑘,𝑖 ki，L；𝑣𝑘,(𝑖−1) ki-1L；kjL；k(𝛿𝑠𝑘+𝛿𝑒𝑘)/2，%；𝑘,𝑗𝑉′ kjL；𝑘,𝑗𝑗𝑉′ jL；𝑗∑𝑁∑𝑘=1

′𝑉𝑘,𝑗𝑉

第j个记录区间经过误差修正后的付油量之和，L；𝑘 (1N)；𝑁 𝑖 2、n)；𝑗 1、、m)，且m=N-1。j𝑻𝒋（4）=(𝑡𝑖+𝑡𝑖−1)/2…………（4）式中：𝑡𝑖 第i𝑡𝑖−1 i-1ji𝟐𝟎𝒊56（）（）𝑉′20,𝑖=𝑉20,𝑠+1⋅𝑉′

×[1+×−20)]+∑𝑚

×[1+𝛽2×(𝑇𝑗−20)]………（5）2 𝑗

𝑗=𝑖+1 𝑗△=𝑉′20,𝑖×[(𝑡𝑖+1−𝑡𝑖)×+(𝑡𝑖−𝑡𝑖+1)×𝛽2]…（6）=+△𝑉′𝑗)×[1+×−20)]………（7）𝑗=𝑖式中：𝑗=𝑖

=+

𝑉20,𝑗×[1+𝛽2×(𝑡𝑖−20)]…（8）′,𝑖 第i（，L；△jL；∑𝑉𝑚 ′∑𝑉𝑗=𝑖+1

第i+1个记录液位下的所有记录区间的容积之和，L；L；jL；∑𝑚∑𝑗=𝑖

第i个记录液位下的所有记录区间标准体积之和，L；-1；-1；𝑚 JJG266CBBB容积表的编制应符合下列规定：1mm1cmL；1JJG266H。概述JJF1059.1C确定度评定方法。依照测量液位高度方式和加油机不确定度分量不同，不确定度评定方式不同，此标准给定典型评定方式。加油机不确定度分量采用MPE：U=0.33%，k=2；：U=0.38%，k=2。加油机不确定分量采用示值误差和重复性误差加油机不确定分量具备明确示值误差及重复性误差，参见附录C评定方法进行评定，并给出不确定度。表A.1加油机法标定记录表表A.1加油机法标定记录表加油机法标定记录表1、基础数据公司标称容积L加油枪序号k平均示值误差𝛿差ത内竖直径mm加油站mm安全高度mm罐号开始高度mm油品记录区间mm2、数据记录液位高度序号i预设液时间位𝐻′𝑖罐内油品温度𝑡𝑖 油枪名油枪名加油机泵码值𝑣𝑘,𝑖L油枪名 油枪名 油枪油枪名mm人工计量mm液位计mm油枪名123456789…n1附录A（）11℃称:称:称:称:称:称:称:𝑣1,𝑖𝑣2,𝑖𝑣3,𝑖𝑣4,𝑖𝑣5,𝑖𝑣…,𝑖𝑣𝑁,𝑖1212校准或核查结果校准或核查校准或核查结果校准或核查加油机油品量器加油机油枪出口量器量器内实际单次示值周期阶段油枪序 名称 容积示值处油温刻度油温体积误差示值误差重复性误差（标定前、后）号LL℃mm℃L𝛿（%）𝛿（%）𝐸𝑛（%）kPAGEPAGE13PAGEPAGE14三次非均匀B

附录B（资料性附录）三次非均匀B样条插值法n+1𝑞𝑖(ℎ𝑖𝑣𝑖)𝑖(0,1B𝑞𝑣)0p(u)=∑n+2di(ℎ𝑖,𝑣𝑖)Ni,3(u) （B.1）0𝑑𝑖(ℎ𝑖𝑣𝑖)𝑖(0,1,𝑛+∈[0,1]为三次Bﻟ (𝑢−𝑢𝑖)3

𝑢≤𝑢＜𝑢I (𝑢𝑖+2−𝑢𝑖)(𝑢𝑖+1−𝑢𝑖)(𝑢𝑖+3−𝑢𝑖)

𝑖 𝑖+1I (𝑢−𝑢𝑖)2(𝑢𝑖+2−𝑢) + (𝑢−𝑢𝑖)(𝑢𝑖+3−𝑢)(𝑢−𝑢𝑖+1) I

(𝑢𝑖+3−𝑢𝑖+1)(𝑢𝑖+2−𝑢𝑖+1)(𝑢𝑖+3−𝑢𝑖)+I (𝑢−𝑢𝑖+1)2(𝑢𝑖+4−𝑢) +I

𝑢𝑖+1

≤𝑢＜𝑢𝑖+2Ni,3(u)=

(𝑢−𝑢𝑖)(𝑢𝑖+3−𝑢)2 + (𝑢𝑖+4−𝑢)(𝑢−𝑢𝑖+1)(𝑢𝑖+3−𝑢)

……………（B.2）❪ (𝑢𝑖+3−𝑢𝑖)(𝑢𝑖+3−𝑢𝑖+1)(𝑢𝑖+3−𝑢𝑖+2) (𝑢𝑖+4−𝑢𝑖+1)(𝑢𝑖+3−𝑢𝑖+1)(𝑢𝑖+3−𝑢𝑖+2)+I (𝑢𝑖+4−𝑢)2(𝑢−𝑢𝑖+2) +I I (𝑢𝑖+4−𝑢)3

𝑢𝑖+2

≤𝑢＜𝑢𝑖+3II (𝑢𝑖+4−𝑢𝑖+1)(𝑢𝑖+4−𝑢𝑖+2)(𝑢𝑖+4−𝑢𝑖+3) 𝑢𝑖+3≤𝑢＜𝑢𝑖+4I{0 其他3记录数据点𝑞𝑖(ℎ𝑖,𝑣𝑖)𝑖(0,1,⋯,𝑛−1)，n是数据记录点总数，采用积累弦长法参数化。𝑢′={𝑢′,𝑢′,𝑢′,𝑢′…𝑢′ } （B.3）0 1 2 3 𝑛−1𝑢′=𝑢′ +（B.4）𝑖 𝑖−1∆𝑝

………（B.5）=|𝑝𝑖−𝑝𝑖−1|=|√(ℎ𝑖−ℎ𝑖−1)2+(𝑣𝑖−𝑣𝑖−1)2|………（B.5）𝑖=1𝑖=1式中：

𝑖 𝑠𝑖 𝑖(1,⋯,𝑛−1)；

∑𝑛−1|√(ℎ𝑖−ℎ𝑖−1)2+(𝑣𝑖−𝑣𝑖−1)2|𝑆 ∆𝑝𝑖为向前差分矢量。0注1：𝑢′=00首尾端重复度为3，即分别为𝑢0=𝑢1=𝑢2=0,𝑢𝑖+3=𝑢𝑖+4=𝑢𝑖+5=1。𝑢𝑖={𝑢0,𝑢1,𝑢2,𝑢′,𝑢𝑖+3,𝑢𝑖+4,𝑢𝑖+5} （B.6）记录数据点𝑞𝑖(ℎ𝑖,𝑣𝑖)𝑖(0,1,⋯,𝑛−1)和𝑢𝑖(0,1,⋯,𝑛+6)代入式（B.7）求得系数矩阵对应元素𝑎𝑖、𝑏𝑖𝑐𝑖、𝑒𝑖，𝑎𝑖、𝑏𝑖𝑐𝑖、𝑒𝑖（B.8）𝑑𝑖(0,1𝑛1)。𝑎��I

=(𝑢𝑖+3−𝑢𝑖+2)2𝑢𝑖+3−𝑢𝑖(𝑢𝑖+3−𝑢𝑖+2)(𝑢𝑖+2−𝑢𝑖)

(𝑢𝑖+2−𝑢𝑖+1)(𝑢𝑖+4−𝑢𝑖+2)𝑖I𝑖=

𝑢𝑖+3−𝑢𝑖 (𝑢𝑖+2−𝑢𝑖+1)2𝑢𝑖+4−𝑢𝑖+1

𝑢𝑖+4−𝑢𝑖𝑒𝑖=(𝑢𝑖+3−𝑢𝑖+1)𝑞𝑖−1I❪𝑑0=𝑞0,𝑑𝑛+2=𝑞𝑛I3I𝑎1=1,𝑏1=𝑐1=0,𝑒1=𝑞0+I3

𝑢4−𝑢3

𝑢𝑛+3−𝑢𝑛+2

……………（B.7）𝑎𝑛+1=𝑛+1=0,+1=1,𝑛+1=𝑛−{𝑖∈{2,3,⋯𝑛}

3 F𝑎1

1F𝑑11

𝑒1F 1l𝑎2

Il𝑑2I lIl⋱ ⋱

Il⋮I=l⋮

I （B.8）l 𝑎𝑛 Il𝑑𝑛

l𝑒𝑛I

[ 𝑎𝑛+1 𝑏𝑛+1 𝑐𝑛+1][𝑑𝑛+1]

[𝑒𝑛+1]4先确定h在数据点列中的区间，若ℎ𝑖≤ℎ＜ℎ𝑖+1，则该数据点对应的节点矢量𝑢𝑖。代入式（B.1）中组成式（B.9）计算出该数值点对应的节点矢量𝑢𝑖的值。0即∑𝑛+2ℎ𝑖𝑁𝑖,3(𝑢𝑖)=ℎ （B.9）0式中ℎ𝑖 (0,1𝑛+2)。将𝑢𝑖（B.1）0𝑃(𝑢)=∑𝑛+2𝑑𝑖(ℎ𝑖,𝑣𝑖)𝑁𝑖,3(𝑢)=(ℎ,𝑣) （B.10）0v就是h例：设数值点为𝑞0(400𝑚𝑚,2300L),𝑞1(450𝑚𝑚,2500L),𝑞2(500𝑚𝑚,2700L),𝑞3(600𝑚𝑚,3200L),𝑞4(800𝑚𝑚,3900L),𝑞5(900𝑚𝑚4600L)𝑞6(950𝑚𝑚5200L)𝑞7(1050𝑚𝑚5200L)𝑞8(1200𝑚𝑚v。𝑞0~𝑞8代入式（B.5）分别求出个相邻两点欧氏距离：|𝑞1−𝑞0|=206|𝑞2−𝑞1|=206|𝑞3−𝑞2|=510|𝑞4−𝑞3|=728|𝑞5−𝑞4|=707|𝑞6−𝑞5|=602|𝑞7−𝑞6|=806|𝑞8−𝑞7|=10111计算𝑠=∑𝑛−1|𝑞𝑖−𝑞𝑖−1|=206+206+510+728+707+602+806+1011=47761计算𝑢′={0,206,412,922,1650,2357,2959,3765,1}，代入式（B.6）组成𝑢=4776477647764776477647764776{0,0,0,0,206,412,922,1650,2357,2959,3765,1,1,1,1}。4776477647764776477647764776计算系数矩阵：将所有u和数值点的横坐标代入式（B.7）中可以得到：𝑎1=1,𝑏1=𝑐1=0𝑒=

+𝑢4−𝑢3

=(400,2300)+(206⁄4776−0)⋅(1200,7000)=(417,2401)1 0 3 8 3𝑎9=𝑏9=0,𝑐9=1,𝑒=

−−

=(1200,7000)−(1−3765⁄47769 8 3 8

3 )⋅(1200,7000)=(1115,6506)𝑎2=

(𝑢5−𝑢4)2𝑢5−𝑢2

(412⁄4776−206⁄4776)2412⁄4776−0 =0.022𝑏=(𝑢5−𝑢4)(𝑢4−𝑢2)+(𝑢4−𝑢3)(𝑢6−𝑢4)=(412⁄4776−206⁄4776)(206⁄4776−0)2 𝑢5−𝑢2

𝑢6−𝑢2

412⁄4776−0+(206⁄4776−0)(922⁄4776−206⁄4776)=0.055922⁄4776−0𝑐2=

(𝑢4−𝑢3)2𝑢6−𝑢3

(206⁄4776−0)2922⁄4776−0=0.010𝑒2=(𝑢5−𝑢3)𝑞1=(412⁄4776−0)⋅(450,2500)=(39,216)按照𝑎2,𝑏2,𝑐2,𝑒2计算方式，计算𝑎𝑖,𝑏𝑖,𝑐𝑖,𝑒𝑖,𝑖∈{3,4,⋯8},即可得到系数矩阵d和矩阵e：F 1 1l0.022 0.0550.010 Ill𝑑=ll

0.0590.085 0.006 I0.077 0.1540.028 I0.054 0.1920.054 Il0.0380.1880.049Il0.0640.1990.032Il0.0890.2170.074I[1 ]417F

2401139l75ll

216405II829II𝑒=llll

240248280399

1170I1265I1534I2280I[11156506]由系数矩阵及矩阵e代入式（B.8）𝑑×ll[𝑥9

𝑦1𝑦21𝑦3I𝑦4I𝑦5I=𝑒I𝑦6II𝑦7I𝑦8I𝑦9]求解该方程可以得到控制点𝑑1到𝑑9，由𝑑0=𝑞0,𝑑10=𝑞8,求得控制点𝑑𝑖(ℎ𝑖,𝑣𝑖)𝑖(0,1,…,10)。控制点列：𝑑0(400,2300),𝑑1(417,2401),𝑑2(449,2485),𝑑3(529,2808),𝑑4(594,3270),𝑑5(829,3887),𝑑6(902,4568),𝑑7(947,5258),𝑑8(1070,6128),𝑑9(1115,6506),𝑑10(1200,7000)。h求高度h=410mm处容积v即𝑞′(410,𝑣)根据横坐标h=410mm，𝑞′∈{𝑞0,𝑞1}，𝑞′对应的𝑢′∈3,4B.1（B）将0,1…,6𝑢′=008728712𝑢′B.1到数据点𝑞′。P(0.008472260837713248)=(410,2352.93672238)𝑣=2352.9（。

附录C（资料性附录）加油机法标定地埋罐容积的不确定度评定ℎ𝑖0,𝑖6.3（5（（7（𝑢2(𝑉 )=𝑐2(𝑉′)𝑢2(𝑉′)+𝑐2(ℎ)𝑢2(ℎ)+𝑐2(𝑇)𝑢2(𝑇)+𝑐2(𝑡)𝑢2(𝑡)+𝑢2(𝑓)………（C.1）𝑟式中：

20,𝑖

𝑗 𝑟𝑗

𝑖 𝑟 𝑖

𝑗 𝑟𝑗

𝑖 𝑟𝑖 𝑟𝑗𝑉′ 𝑗L；𝑗ℎ𝑖 ，mm；𝑗𝑡𝑖 𝑖𝑢𝑟(𝑓)为非检定点插值引入的不确定度附加参量。灵敏度系数：𝑐(𝑉′)=1；𝑐(ℎ𝑖)=1；在此认为𝑉′≈𝑉20,𝑗适当的近似计算𝑐(𝑇𝑗)=−1,𝑐(𝑡𝑖)=𝑗 𝑗6.2.2𝑢𝑟(𝑓)C1𝑢2(𝑉 )=𝑢2(𝑉′)+𝑢2(ℎ)+𝑢2(𝑡)+𝑢2(𝑓) （C.2）𝑟 20,𝑖

𝑟

𝑟

𝑟𝑖 𝑟

表C.1标准不确定度分量表序号影响量灵敏系数标准不确定度备注1𝑢𝑟(ℎ𝑖)10.034%液位计方式测量0.085%手工测量2𝑢𝑟(𝑡𝑖)10.007%介质为汽油,材质为碳钢3𝑢𝑟(𝑓)10.050％三次非均匀B样条插值算法注1：采用液位仪测量高度：𝑢𝑟(ℎ𝑖)=0.034%。注2：采用手工测量液位高度：𝑢𝑟(𝐻𝑖)=0.085%。注3：介质温度测量引入的不确定度：𝑢𝑟(𝑡𝑖)=0.007%。注4：插值计算引入的不确定度：𝑢𝑟(𝑓)=0.050％。C.3C）𝑢𝑟(𝑉′)=√𝑢2(𝛿)+𝑢2(𝐸𝑛) （C.3）𝑗 𝑟 𝑟𝑟𝑢(𝐸𝑛)=ത𝑟式中：

………（C.4）√3𝑢𝑟(𝛿) 𝑢𝑟(𝐸𝑛)重复性误差分量；ത 𝑢𝑟(𝛿)当连接地埋罐有多个油枪时，采用所有参与测量的加油机不确定度𝑈𝑟(𝛿1),𝑈𝑟(𝛿2),…,𝑈𝑟