大型儲油罐變形數據處理與分析

賈勇帥,獨知行,郭巍,陳哲(山東科技大學測繪科學與工程學院,山東青島 266590)

大型儲油罐變形數據處理與分析

賈勇帥?,獨知行,郭巍,陳哲
(山東科技大學測繪科學與工程學院,山東青島 266590)

摘 要:通過對某油庫大型儲油罐的觀測實例,介紹了TM30在油罐測量的原理,并分析了確定圓心坐標的方法的精度,選取MATLAB進行程序設計對油罐的變形數據進行自動化處理,得出一系列幾何參數,該方法明顯提高了工作效率。通過這些幾何參數來分析該儲油罐幾何形體的形變情況,為人員生命安全和安全生產提供了依據。

關鍵詞:油罐;TM30;程序設計;數據處理與分析

1 引 言

當今世界大多數國家都采用鋼制儲油罐儲存石油,其中在我國應用最廣泛的是立式鋼制浮頂油罐。隨著科技水平的不斷發展,儲油罐結構的規模一直在向著大型甚至是超大型的方向發展。儲油罐的大型化不僅使我國的戰略石油儲備邁向了一個新的臺階,也給石油石化行業的原油的存儲、輸送帶來了極大的便利,但與此同時,儲油罐的大型化對儲油罐的安全性能也提出了更高的要求,這些更高的要求其中之一就體現在儲油罐的幾何形體變形檢測上[1]。

本文選取的目標儲油罐隸屬于某石油化工集團管道儲運公司,位于我國東部沿海某大型商業儲備油庫中。該油罐于2002年建造完成并投產使用,公稱半徑30 000 mm,屬于典型的雙盤型浮頂油罐。該罐建造年代較早,服役時間較長,受罐內液位的反復升降、地面沉降、沿海風力較大以及日照較強等因素的影響,該儲罐的形體可能已發生垂直傾斜和水平偏移等形體上的變化[2],因此,需對儲罐的幾何形體變形情況重新進行檢測、評估。

2 測量原理與儀器操作

測量選用的儀器是徠卡公司生產的TM30型電子全站儀,TM30是最新的能替代人力進行自動搜索、跟蹤并能精確照準目標同時獲取所需要的距離、三維坐標以及其他有用信息的電子全站儀,通常稱作”測量機器人”[3]。

根據不同的測量環境和條件,我們將檢測儲油罐形體變形的方法分為了罐內測量和罐外測量。本文所進行討論的是以罐內測量為基礎,罐內測量的工作原理是利用測量機器人(以徠卡TM30為例)根據儀器自帶的罐內測量程序自動進行全角度的掃描測點并記錄數據而進行的。

圖1　罐內測量原理示意圖

如圖1所示,測量時只需將測量機器人架設于浮頂中心位置,打開儀器中的罐內測量程序,并測定或設置以下幾個參數即可開始測量:

(1)壁板的圈數以及各圈板的高度;

(2)基準點和起始點位置及高程;

(3)進行偏心改正的隨機選取的罐壁三點坐標;

(4)鏡頭轉動的角度步進值(一般為7．5°);

參數設置完畢后,儀器將進行自檢,之后開始數據采集。我們將得到各圈壁板1/4和3/4處的均勻分布的所需測點坐標(x,y,z),這也就意味著每一圈鋼板的形體參數我們將用這圈鋼板的1/4和3/4處來共同描述,這樣做可以有效地減少偶然誤差。

3 程序設計

3．1確定圓心坐標方法的精度分析

儲罐形體偏差計算首先需要設定一個基準進行對比,通常情況下我們都是以起始圈板為基準圈板。從儀器中導出的數據是各圈板1/4和3/4處一系列離散的點的坐標。我們將首先利用起始圈所測得的離散點,這些點在理論上可以擬合成一個基準圓。從起始圈各測點坐標開始,擬合成基準圓并求出該基準圓圓心坐標,然后根據擬合圓心坐標擬合出空間直線進行傾斜度的計算。

利用已知點的坐標來擬合圓的方法非常多,可以直接分別計算48個離散點x,y坐標平均值進行簡單的圓心估算。

以及通過最小二乘擬合求圓心的坐標,圓的基礎方程為:

經過變換可得:r

根據最小二乘準則:

達到最小就能構造出一個圓,為了方便計算,特將公式改寫為:

通過表1比較可知采取最小二乘擬合圓求得相應的擬合圓心精度較高[4],因此選取最小二乘擬合法確定圓心坐標。

圓心坐標精度對比結果 表1

3．2幾何參數數學模型

徑向偏差反映了儲油罐在水平方向上偏轉變形程度,若某個方向上的徑向偏差超過一定的限差,則該油罐存在安全隱患。各離散點到擬合圓心的半徑與標準半徑的差即徑向偏差。

圈板橢圓度的計算是以徑向偏差為基礎的,兩個對徑點到擬合圓心距離之和可以構成這圈壁板擬合圓的一條直徑。圈板橢圓度實際上是徑向偏差的進一步推算,規定中并沒有單獨列出橢圓度的限差值,只需作為徑向偏差的一個參考項即可。

3．3MATLAB數據處理程序設計

TM30采集得到的各離散點坐標數量為48?n個(n為掃描罐體的圈數),通過Excel也能進行數據的處理得到需要的幾何參數,但是由于數據量過大,需要進行反復的人工計算,在計算的過程中容易由于人為的原因產生錯誤,耗時耗力。本文通過MATLAB[7,8]進行程序設計預期達到的功能如圖2所示:

在理想狀態下,儲油罐是完全垂直于地面的。當儲罐發生傾斜時,其各圈板的擬合圓心也會發生改變,圓心之間的連線也就不再垂直于地面。若將各圈圈板所構成的擬合圓的圓心擬合成一條空間直線[5,6],求出直線與平面的夾角就是傾斜度。

zi-zo

圖2　程序設計流程圖

由于主要是對數據進行處理,所以可視化界面選取的控件為兩個Push Button分別為“數據處理”、“數據檢測”和兩個Edit Text控件分別用來“輸入n”和“顯示空間直線方程”。

[filename,pathname] = uigetfile('?．xlsx','選擇數據');%檢測測量數據是否有遺漏或者重復。

fid=[pathname,filename]

n=get(handles．edit8,'string');n=str2num(n) for i=48:48:48?n

C=['c',num2str(i),':','c',num2str(i)] str1=xlsread(fid,C);str1=num2str(str1) a=strcmp(str1,'48') if a= =0

msgbox(strcat('第',num2str(i-48),'到',num2str (i),'行出錯了'))

end end n=get(handles．edit1,'string');%然后輸入n來確定掃描的圈數

n=str2num(n)

[filename,pathname] =uigetfile('?．xlsx','選擇數據')

fid=[pathname,filename]

K=1

for i=1:48:48?(n-1)+1; % TM30采集的點的坐標數據是48?n個(n為對油罐進行了掃描圈數的次數),所以需要反復進行n次數據的讀取。

x=['d',num2str(i),':','d',num2str(i+47)] x=xlsread(fid,x);x=[x]'

y=['e',num2str(i),':','e',num2str(i+47)] y=xlsread(fid,y);y=[y]'

h=['f',num2str(i),':','f',num2str(i+47)]

z=xlsread(fid,h);z=[z]'

z=mean(z) [x0,y0] =circ(x,y);yx=[x0,y0] for j=1:24%計算所需幾何參數

r(j)= sqrt((x(j)-A)^2+(y(j)-B)^2)

r(j+24)= sqrt((x(j+24)-A)^2+(y(j+24)-B)^2) d(j)= r(j)+r(j+24) end

cha=max(d)-min(d);du=cha/ mean(d) k=k+1

F=['f',num2str(k)]

j=['a',num2str(i),':','b',num2str(i)]

[data,str] =xlsread(fid,j);data=num2str(data) bt=str{2};bt=strcat({data},{bt}) ma=max(r)-30;mi=min(r)-30; if abs(ma)＞abs(mi)%判斷并找出半徑的最大測量誤差。

xlswrite('C:UsersAdministratorDesktop1．xlsx',1000?ma,'sheet1',F); else

xlswrite('C:UsersAdministratorDesktop1．xlsx',1000?mi,'sheet1',F);

end end

4 實例數據分析

由國家能源局發布的SY/ T5921-2011《立式圓筒形鋼制焊接油罐操作維護修理規程》中第5．9．2．3條規定了大型立式常壓儲罐底圈板和其余圈板內表面上任意點半徑的允許偏差(徑向偏差),其中底圈板以上各圈板的半徑允許偏差為底圈板半徑允許偏差的3倍。其中45 m

徑向偏差誤差對比結果 表2

在表2上可以明顯看出罐體半徑的最大測量偏差在誤差允許范圍內,但是隨著高度增加罐體的最大測量偏差從外傾向內傾變化。在圖3上可以明顯看出紅色標記的直徑中上部罐體存在較大的變形,所以有可能是因為儲罐在長期高罐位貯油的情況下,油壓力對中部圈板的壓迫導致了中上部圈板的輕微變形。

為了對油罐的形體變形情況進行全面的了解,還需要對油罐的傾斜度進行分析,利用圈板擬合圓心數據進行空間直線的擬合,進一步得到傾斜度。

從表3上我們可以看出,x有增大的趨勢,y有增大的趨勢,若將基準圓圓心作為參考,其他各圈板1/4處和3/4處的圓心都位于基準圓心的東北面,說明該罐有向東北方向傾斜的趨勢,我們可以更直觀的在圖4上看出。

圖3　罐體俯視圖

圈板擬合圓心數據 表3

圖4　擬合圓心空間圖

根據MATLAB擬合得到空間直線方程,主要代碼如下:

隨機選取兩點坐標( -0．0096, 0．0058,0),( -0．0089,0．0065, 1 ),利用式( 9 )可得傾斜度為0．056 7°,遠低于JJG168-2005《立式金屬罐容量檢定規程》第5．3條規定:“罐體傾斜度不得超過1°”。

T=[z;ones(1,length(z))]; %擬合空間直線方程

N=T?A';O=T?B';M=T?T';a=(MN)';b=(MO)' a=['x=',num2str(a(1)),'?z+',num2str(a(2))]

b=['y=',num2str(b(1)),'?z+',num2str(b(2))] set(handles．edit2,'string',{a,b})

最后得到擬合圓心在空間的直線方程:

5 結 語

本文通過MATLAB進行了數據處理,相比Excel處理數據,通過本程序不但工作效率有了明顯提升,而且通過對數據分析可以看出罐體存在微小的幾何形體變形,徑向偏差較小,符合規定,有較輕微的傾斜趨勢,傾斜方向為東偏北,傾斜度為0．056 7°,也在誤差允許范圍內,證明此罐的狀況良好。

參考文獻

[1] 武銅柱．大型立式油罐發展綜述[J]．石油化工設備技術,2004(3):56～59．

[2] 張玉平,陶彬,郎需慶等．大型儲罐變形評估技術研究及應用[J]．安全、健康和環境,2014(5):43～47,53．

[3] 劉朋姣,侯東亞．徠卡TM30在罐體掃描檢測中的開發及應用[J]．測繪通報,2010(10):73～74．

[4] 張應芳,張愛國,蘇積聰等．觀測圓上多個任意點及其圓心坐標擬合[J]．黑龍江工程學院學報,2014(2):27～31．

[5] 襲楊．空間直線擬合的一種方法[J]．齊齊哈爾大學學報,2009(2):64-68．

[6] 姚宜斌,黃書華,孔建等．空間直線擬合的整體最小二乘算法[J]．武漢大學學報·信息科學版,2014(5):571～574．

[7] 羅華飛．MATLAB GUI設計學習手記(第二版)[M]．北京:北京航空航天大學出版社,2011．

[8] 阮沈勇,王永利,桑群芳．MATLAB程序設計[M]．北京:電子工業出版社,2003．

Processing and Analysis of Deformation Data of Large Oil Storage Tank

Jia Yongshuai,Du Zhixing,Guo Wei,Chen Zhe
(Shandong University of Science and Technology,Qingdao 266590,China)

Abstract:the principle of TM30 in the oil tank measurement is introduced through a observational instances of large storage tanks in this paper, analyze the Precision of the Method which determines the center of the circle and a series of geometric parameters are obtained by using the program design of MATLAB．The method has obviously improved the working efficiency．In order to ensure the life safety and safety in production the deformation of the geometrical shape of the oil tank is analyzed．

Key words:oil tank; TM30; program design; data processing and analysis

文章編號:1672-8262(2015)06-117-04中圖分類號:P258

文獻標識碼:B

收稿日期:?2015—09—25

作者簡介:賈勇帥(1991—),男,碩士研究生,研究方向變形監測。