通过总结全站仪坐标测量原理及误差分析,提出将其用于大型立式储油罐体变形检验的一种检测方法,总结了常见罐体圆度、垂直度、凹凸变形等尺寸误差检验的工作流程和数据分析方法,同时还讨论了浮顶变形相关的检测方法,并通过实际工程中的实例罐体检验验证,证实了全站仪坐标测量方法用于立式储罐变形检验的适用性和有效性,为今后圆筒状立式储油罐体的变形检验工作提供了一种便捷有效的检验方法。 By summarizing the principle of total station coordinate measurement and error analysis, a method of measuring deformation of large vertical oil storage tank is put forward. The workflow and data analysis methods of measuring dimensional errors, such as roundness, verticality, concave-convex deformation of tank, are summarized. Meanwhile, the detection methods related to floating roof deformation are also discussed. The applicability and validity of the total station coordinate measuring method for deformation inspection of vertical storage tanks are verified by a practical engineering example. It provides a convenient and effective method for deformation inspection of cylindrical vertical storage tanks in the future.
黄正均1,2*,秦琨1*,张磊1,2,刘钰1,2,张栋1,2,侍吉清3
1北京科技大学,土木与资源工程学院,北京
2城市地下空间工程北京市重点实验室,北京
3甘肃蓝科石化高新装备股份有限公司,甘肃 兰州
收稿日期:2018年10月22日;录用日期:2018年11月9日;发布日期:2018年11月16日
通过总结全站仪坐标测量原理及误差分析,提出将其用于大型立式储油罐体变形检验的一种检测方法,总结了常见罐体圆度、垂直度、凹凸变形等尺寸误差检验的工作流程和数据分析方法,同时还讨论了浮顶变形相关的检测方法,并通过实际工程中的实例罐体检验验证,证实了全站仪坐标测量方法用于立式储罐变形检验的适用性和有效性,为今后圆筒状立式储油罐体的变形检验工作提供了一种便捷有效的检验方法。
关键词 :坐标测量,储油罐体,变形检验,检测
Copyright © 2018 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
钢制立式储油罐体作为原油储备运输的重要设备,其具有省材料、占地少、投资省、储量大、易操作等优点,也是大多数国家广为采用的石油储存设备之一 [
现有的储罐变形检测主要针对外观质量检测,对其完整的变形情况往往无法进行准确、有效、经济高效的获取,又或需采用三维激光扫描等复杂技术才能得以实现,但设备成本过高、后期数据处理过于复杂,导致针对罐体的各项变形检验尚无快捷有效的方法 [
全站仪为当前工程测量所应用的主要设备之一,其在土木、交通、矿山、水利、能源等各个领域均应用广泛。其空间点位坐标测量原理为空间球坐标测量的原理,即已知一现有点坐标,假设全站仪对目标待测点的水平角、垂直角(天顶距)、斜距分别为(α, β, D),根据右手坐标系可得出待测点P的空间三维坐标为:
{ x = D sin β cos α y = D sin β sin α z = D cos β (1)
测量示意图如图1所示。此外,根据范百兴等 [
序号 | 测量距离/m | mpc/mm | mp/mm | mz/mm | 备注 |
---|---|---|---|---|---|
1 | 10 | 0.17 | 0.24 | 0.17 | |
2 | 30 | 0.18 | 0.25 | 0.17 | |
3 | 50 | 0.20 | 0.27 | 0.20 | |
4 | 120 | 0.32 | 0.41 | 0.25 | |
5 | 200 | 1.09 | 1.52 | 1.03 | |
6 | 300 | 1.60 | 2.23 | 1.52 |
表1. TS30全站仪坐标点位测量误差估算值表
注:其中mpc为平面坐标误差;mp为空间坐标误差;mz为Z方向垂直坐标误差。
大型储油罐体外观变形检验主要包括罐体几何圆度、罐壁垂直度、罐壁凹凸度、导向油管垂直度、以及浮顶相关的变形检测等内容。主要工作方式为根据全站仪坐标测量原理,对罐体各部分组件设立若干测点,并进行平面坐标测量和三角高程测量,再根据测点坐标,建立局部坐标系,并采用平差法进行坐标计算并绘图,以最终确定罐体各部分的几何变形情况。
测量工作可根据《立式圆筒形钢制焊接储罐施工规范》(GB50128-2014)的有关要求 [
图1. 观测坐标示意图
测量数据采用平差法进行计算,平差法为测量数据处理中最常采用的方法之一,包括条件平差、间接平差等,所有经典的平差计算均是基于最小二乘法原则对测量数据误差进行调整计算 [
如已知点坐标为: A ( X A , Y A ) , B ( X B , Y B ) ,待定点:P (需计算其坐标),观测数据:a, b (g = 180 − a − b),计算各边方位角:
α A B = tan − 1 y B − y A x B − x A α A P = α A B − α α B P = α B A + β (2)
其各边长度得:
c = ( x B − x A ) 2 + ( y B − y A ) 2 a = c ⋅ sin α sin γ , b = c ⋅ sin β sin γ (3)
则待测点P的坐标如下:
x P = x A + b ⋅ cos α A P y P = y A + b ⋅ sin α A P x P = x B + a ⋅ cos α B P y P = y B + a ⋅ sin α B P (4)
1) 圆心点坐标计算
根据圆的一般公式: x 2 + y 2 + D x + E y + F = 0 ,为达到测量数据与拟合函数之间误差的平方和为最小,可将其改写为:
M = ∑ i = 1 n ( x i 2 + y i 2 + D x i + E y i + F ) 2 (5)
令:
∂ M ∂ D = ∂ M ∂ E = ∂ M ∂ F = 0
则可将式(5)转换为线性方程组,整理合并为:
{ D ∑ i = 1 n x i 2 + E ∑ i = 1 n x i y i + F ∑ i = 1 n x i = − ∑ i = 1 n x i 3 − ∑ i = 1 n x i y i 2 D ∑ i = 1 n x i y i + E ∑ i = 1 n y i 2 + F ∑ i = 1 n y i = − ∑ i = 1 n x i 2 y i − ∑ i = 1 n y i 3 D ∑ i = 1 n x i + E ∑ i = 1 n y i + n F = − ∑ i = 1 n x i 2 − ∑ i = 1 n y i 2 (6)
求解方程组(6)可获得D、E、F的相应数值,并由此可以得出圆心坐标为 ( − D 2 , − E 2 ) 。
2) 圆度误差计算
根据前述圆心坐标计算中求解计算得出的圆心坐标为 ( − D 2 , − E 2 ) ,圆半径为 D 2 + E 2 − 4 F 4 ,并根
据拟合圆心和半径计算出圆心坐标至各测点的实测半径值,据此计算储油罐体测量断面各测点的半径偏差以及圆度误差,圆度误差计算方法为最大半径与最小半径之差和拟合半径的比值,以百分比表示,计算示意图如图2。
图2. 圆度误差计算示意图
图3. 局部凹凸变形测量及计算示意图
3) 罐壁局部凹凸变形计算
通过现场测量点进行拟合得出圆心坐标与半径后,计算罐壁各点对应半径,即可得出罐壁各点位置的局部凹凸变形值,计算示意图如图3所示。但是,由于罐壁现场测量点有限,现有的检测方法仅根据所测测点值给出有限点的变形情况,但是无法准确体现罐壁圆弧表面的全面变形情况,故可采用三次样条函数插值计算方法进行拟合计算圆弧表面的变形,计算可采用MATLAB编程进行,据此可得出整个罐壁圆弧表面的局部凹凸变形情况。
4) 导向油管和罐壁垂直度计算
导向油管和罐壁垂直度计算方法均利用投影在同一点上的罐顶观测点和罐底观测点坐标进行计算 [
Δ = ( x 1 − x 2 ) 2 + ( y 1 − y 2 ) 2 H × 100 % (7)
图4. 罐壁垂直度测量及计算示意图
5) 扶梯中心线与导轨中心线偏差计算
储油罐体内部设有通行扶梯,扶梯下方与位于浮顶上方的导轨相接,以供人员检修等使用,由于罐体体积大、高度深,扶梯在使用过程中会出现晃动、偏移等,导致其与下部导轨会出现偏差甚至脱离,故在进行变形检验时也需对其中心偏差进行检测校对。扶梯中心线与导轨中心偏差检验可采用多点测量的方式,即测量各点的坐标,再进行平面投影,以做图形比对,进而从成果图上量测其偏差值,如图5所示。
图5. 扶梯与导轨中心偏差测量示意图
6) 浮顶部分变形检验计算
立式储油罐体内设钢制浮顶,其结构分上下两层,中部设有浮舱,浮顶下盘设有钢制立柱,以便原油排空后浮顶下落支撑至罐底;浮顶上盘上部设有排水孔及浮舱人孔人若干。油罐在使用过程中往往因为各种因素会导致浮顶出现不同程度的变形,因此也需进行检验计算,现有的凹凸变形检测方法主要采用1 m长直式样板进行间歇测量,所测结果仅为相对变形值,对某一区域整体大面积变形或浮盘各区域准确变形值尚无法得出,本文提出的主要采用方法为下盘凹凸变形利用立柱支撑点坐标实测与设计值对比,上盘采用直接测量的方式进行计算分析,可得出浮盘整体的凹凸变形情况。
某储油基地43号成品油储罐为立式双盘内浮顶钢制圆柱形储罐,公称容积100,000 m3,直径为80 m,罐壁高度为21.8 m,由9层不同厚度钢板焊接而成,壁板公称厚度从下至上分别为:6~9#层12 mm、5#层15 mm、4#层18.5 mm、3#层21.5 mm、2#层27 mm和1#层32 mm,单层高度2.44 m (1#层除外)。内浮顶为上下两层双盘式结构,上浮盘中心与边缘高,中间圈层底,设计坡度1.5%;下浮盘为同一设计高度,浮盘与罐底设有8层支柱支撑,罐底中心高,外侧低,坡度为0.8%,其中下浮盘边缘距罐底上表面原设计高度1.8 m,如图6所示。自2009年开始使用起,经过8年多使用,期间未对其进行过任何变形及质量检测,2017年中由于罐内积水过多且浮盘排水孔堵塞导致浮盘发生整体倾覆,肉眼可见罐壁内侧有大量擦痕,且下浮盘有明显凹凸变形,如图7所示。故此,罐体各部件变形及是否会影响安全生产情况未知,需对该储罐进行全面变形检测,以掌握其整体变形情况,为储罐修复和安全生产提供技术支持和参考。拟对其进行检测的主要内容包括:1) 储罐不同壁厚处罐壁的几何圆度;2) 储罐罐壁及整体垂直度;3) 储罐罐壁局部凹凸变形;4) 导向油管的垂直度;5) 扶梯中心线与导轨中心线偏差;6) 浮顶下浮盘凹凸变形;7) 浮顶上部排水孔与倾斜变形部分浮舱人孔高差。
测量选用徕卡TS30全站仪,采用相对坐标系,原点设在罐底中心,以扶梯中心为北向。针对罐壁从底至顶布设了6个断面,每层24个测量点;扶梯从下往上,布设测点5个,导轨从扶梯连接处往罐内,布设测点3个;浮顶下盘凹凸变形测量采用沿支撑立柱布设测点的方式,共布置8圈188个测点,各圈测点起点均以导向油管及扶梯下方的出入口为起点(测量北向),同时针对南侧变形较大区域进行了加密测量,共布设加密点20个。根据设计原图,浮顶下盘原设计标高均为同一个标高水平,为一水平圆盘,故在对其凹凸变形检测时以中心高程为基准水平,进行比对计算;上浮盘检测排水孔2个,紧急排水孔4个以及浮舱人孔23个,主要针对浮顶倾斜变形区域进行检测。
图6. 10万立方储罐外观图片
图7. 储罐内部状态照片
根据前述检测和计算方法,获得了该立式储油罐体的各项变形情况,具体结果见表2、表3所示,其浮顶下盘凹凸变形检验结果及示意见图8。
测量 层位 | 最大 半径 偏差 /mm | 最大 半径 偏差 /mm | 圆度/% | 是否 符合 标准 | 最大 凹凸 偏差 /mm | 最大 凹凸 偏差 /mm | 是否 符合 标准 | 圆心 垂直 度(%) | 圆心 垂直 度误 差/mm | 偏移 方向 | 是否 符合 标准 | 北侧 导向 油管 垂直 度/% | 北侧导 向油管 垂直度 误差/mm | 偏移 方向 | 南侧 导向 油管 垂直 度/% | 南侧导 向油管 垂直度 误差/mm | 偏移 方向 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
第一层 | 50.369 | −35.502 | 0.215 | 否 | 12.298 | −17.449 | 否 | 0.112 | 20.396 | 北 | 是 | 0.174 | 35.355 | 东 | 0.075 | 15.133 | 东 |
第二层 | 66.807 | −24.850 | 0.229 | 否 | 4.074 | −10.952 | 否 | ||||||||||
第三层 | 58.507 | −60.245 | 0.297 | 否 | 10.954 | −10.638 | 是 | ||||||||||
第四层 | 39.186 | −66.129 | 0.263 | 否 | 8.029 | −15.375 | 否 | ||||||||||
第五层 | 35.460 | −40.726 | 0.191 | 否 | 8.485 | −7.388 | 是 | ||||||||||
第六层 | 61.213 | −2.403 | 0.159 | 否 | 0.502 | −12.788 | 是 |
表2. 43号立式储罐外观变形检测结果
测量 层位 | 最大 凹凸 变形/m | 最大 凹凸 变形/m | 平均/m | 人舱浮 孔与北 侧排水 孔高差 最小值/mm | 人舱浮 孔与北 侧排水 孔高差 最大值/mm | 人舱浮 孔与南 侧排水 孔高差 最小值/mm | 人舱浮 孔与南 侧排水 孔高差 最大值/mm | 扶梯与 导轨中 心线偏 差最小 值/mm | 扶梯与 导轨中 心线偏 差最大 值/mm |
---|---|---|---|---|---|---|---|---|---|
A | 0.041 | 0.008 | 0.024 | 245 | 302 | 476 | 493 | 20 | 60 |
B | 0.018 | −0.042 | 0.005 | ||||||
C | 0.012 | −0.051 | −0.015 | ||||||
D | 0.005 | −0.049 | −0.022 | ||||||
E | −0.002 | −0.049 | −0.028 | ||||||
F | −0.017 | −0.078 | −0.051 | ||||||
G | −0.072 | −0.124 | −0.093 | ||||||
H | −0.079 | −0.155 | −0.109 | ||||||
加密区 | 0.039 | −0.092 | −0.014 |
表3. 43号立式储罐浮顶变形检测结果
从表2、表3以及图8可以看出,储罐罐壁最大凹凸变形值分别为12.298 mm和−17.449 mm,罐体圆度误差最大为0.297%,罐体圆心垂直度为0.112%,导向油管垂直度分别为0.174%和0.075%,浮顶下盘倾斜变形区最大凹凸变形为41 mm和−155 mm,排水孔与人舱浮孔最小高差为245 mm。参考《立式圆筒形钢制焊接储罐施工规范》(GB50128-2014)的相关规定 [
1) 各层罐壁圆度误差均超出施工验收规范要求。
2) 储罐壁板局部凹凸偏差中,第一层(壁厚32 mm)、第二层(壁厚27 mm)及第四层(壁厚18.5 mm)超出验收规范要求,其余各层符合要求。
3) 罐壁垂直度在北、西方向略超出验收规范要求,储罐整体垂直度符合规范要求,偏移方向为北侧(即扶梯侧方向)。
4) 对比原设计资料,浮顶下盘整体下沉,且下沉量从内往外逐步增大;各测点高差负值表明浮盘相对原标准高度有下沉或凸出鼓起,高差为正值表明浮盘相对原高度上隆或向上凹陷隆起。其中A/B/C/D三层变形上凹变形量有超过原设计高度,整体各层最大凹凸变形差值约60 mm,大致圈定三个凹凸变形较大区域位置如图8示。
5) 各浮舱人孔与北侧(扶梯侧)排水孔的高差最小为0.245 m,最大为0.302 m;与南侧排水孔的高差最小为0.476 m,最大为0.493 m。各人孔与北侧和南侧排水孔的高差明显不一致,北侧低于南侧,说明此区域浮盘较南侧有明显变形,与浮顶下盘凹凸变形测量结果基本一致。
图8. 浮顶下盘凹凸变形检测示意图
本文通过分析全站仪坐标测量的原理及误差估算,提出基于坐标测量的储油罐体变形检验方法,总结出罐体圆度、垂直度、凹凸变形、浮顶变形等检测内容的工作流程和数据分析处理方法,并通过某一实际工程案例的立式储罐变形检验验证,获得了较为准确有效的储罐变形数据,并为其使用修复和安全生产提供了数据参考,同时证实了全站仪坐标测量方法用于立式储罐变形检验工作的适用性和有效性,也为今后类似各型储油罐体的变形检验工作提供了一种经济、便捷、有效的检验方法。
中央高校基本科研业务费(FRF-BD-17-006A);国家重点研发计划项目资助(2016YFC0600703)。
黄正均,秦 琨,张 磊,刘 钰,张 栋,侍吉清. 基于全站仪坐标测量的立式储罐变形检验方法研究Research on Deformation Detection Method of Vertical Storage Tank Based on Total Station Coordinate Measurement[J]. 土木工程, 2018, 07(06): 853-862. https://doi.org/10.12677/HJCE.2018.76103