(東北石油大學(xué) 機(jī)械科學(xué)與工程學(xué)院， 大慶 163318)

大型鋼制儲(chǔ)罐因具有薄壁、耐高壓和良好的密閉性等特性，而廣泛應(yīng)用于石油石化行業(yè)中[1]。由于大型儲(chǔ)罐大多建造在沿海等比較松軟的土地上，地基的不穩(wěn)定易使儲(chǔ)罐發(fā)生不同程度的沉降。這種不均勻沉降會(huì)影響儲(chǔ)罐的正常運(yùn)行，嚴(yán)重時(shí)會(huì)帶來(lái)重大的安全隱患[2-4]。儲(chǔ)罐在不均勻沉降過(guò)程中，罐底板與罐壁板會(huì)發(fā)生相應(yīng)的應(yīng)力集中和徑向位移變化，從而使罐底板產(chǎn)生扭曲以及罐壁板產(chǎn)生內(nèi)凹或外凸的現(xiàn)象，也會(huì)使罐底板與罐壁板的焊縫處發(fā)生撕裂[5-6]。MALIK[7]提出使用傅里葉諧波法對(duì)儲(chǔ)罐的不均勻沉降進(jìn)行模擬。JONAIDI等[8-9]對(duì)不同類型的大型鋼制浮頂型儲(chǔ)罐進(jìn)行了線性的靜力學(xué)研究。

浙江大學(xué)空間結(jié)構(gòu)研究中心分析了不同類型鋼儲(chǔ)罐在罐周不均勻沉降下的結(jié)構(gòu)響應(yīng)，并且指出目前的沉降控制標(biāo)準(zhǔn)不合理，應(yīng)該充分考慮地基變形對(duì)儲(chǔ)罐結(jié)構(gòu)的影響[10]。在儲(chǔ)罐檢測(cè)方面，東北石油大學(xué)李偉、戴光等[11]利用聲發(fā)射技術(shù)對(duì)儲(chǔ)罐底板狀態(tài)進(jìn)行在線檢測(cè)，取得了一定的研究結(jié)果。大量研究表明，儲(chǔ)罐發(fā)生不均勻沉降后會(huì)加速罐底板的腐蝕。因此加強(qiáng)對(duì)儲(chǔ)罐不均勻沉降和儲(chǔ)罐底板腐蝕情況的檢測(cè)，對(duì)維護(hù)儲(chǔ)罐的安全運(yùn)行有非常重要的作用。

筆者針對(duì)儲(chǔ)罐不均勻沉降條件下底板的受力狀態(tài)和腐蝕速率變化進(jìn)行數(shù)值模擬和試驗(yàn)研究，結(jié)合聲發(fā)射技術(shù)研究不同應(yīng)力條件下底板腐蝕速率與聲發(fā)射特征之間的對(duì)應(yīng)關(guān)系，并利用均方差得到對(duì)于不均勻應(yīng)力腐蝕最為敏感的聲發(fā)射參量，以為相應(yīng)條件下聲發(fā)射在線監(jiān)測(cè)提供技術(shù)依據(jù)。

圖1 不同沉降面積的儲(chǔ)罐應(yīng)力分布圖

1 儲(chǔ)罐不均勻沉降數(shù)值模擬

應(yīng)用ABAQUS有限元軟件進(jìn)行儲(chǔ)罐不均勻沉降應(yīng)力狀態(tài)數(shù)值模擬，分析在不同沉降面積的條件下，儲(chǔ)罐罐底板與罐壁板的應(yīng)力變化以及應(yīng)力集中部位。筆者選擇30°，45°，60°扇形等3個(gè)不同沉降面積進(jìn)行儲(chǔ)罐有限元數(shù)值模擬分析(見(jiàn)圖1)，液位為滿液位狀態(tài)。

從圖1可看出當(dāng)滿液位時(shí)，沉降面積為30°扇形面積的罐壁最大MISES應(yīng)力(米塞斯應(yīng)力)主要分布在底部靠近焊縫處，且罐底邊緣板部分的應(yīng)力明顯大于中部未發(fā)生沉降部位的。沉降面積為45°扇形面積的罐壁板最大應(yīng)力位置與30°的類似，只是分布范圍有所增大，邊緣板應(yīng)力也呈現(xiàn)增大趨勢(shì)；沉降面積為60°扇形面積的罐壁板的最大應(yīng)力分布區(qū)域較大，但最大應(yīng)力值為1.478 MPa，較之前的應(yīng)力值減小。由此可以看出，受不均勻沉降影響，罐底邊緣板和罐壁處均會(huì)出現(xiàn)應(yīng)力增大趨勢(shì)，且主要分布位置集中于罐壁與罐底交界處以及邊緣板附近。

2 不均勻沉降時(shí)的應(yīng)力與腐蝕速率試驗(yàn)

針對(duì)數(shù)值模擬結(jié)果，儲(chǔ)罐在出現(xiàn)不均勻沉降情況下，罐底邊緣板處會(huì)呈現(xiàn)不均勻的分布狀態(tài)，由于儲(chǔ)罐中的原油等液體介質(zhì)多具有腐蝕性，且研究表明，不均勻應(yīng)力往往會(huì)加速腐蝕的發(fā)展。因此為了得到在不均勻應(yīng)力狀態(tài)下儲(chǔ)罐底板的腐蝕變化情況，試驗(yàn)以常用的Q235材料立式儲(chǔ)罐為試件，設(shè)計(jì)受力模擬結(jié)構(gòu)，并在其中放入腐蝕液，以研究試件在不同濃度的酸性介質(zhì)腐蝕條件下應(yīng)力變化與腐蝕速率之間的對(duì)應(yīng)關(guān)系，得到腐蝕速率隨試件應(yīng)力變化的演化趨勢(shì)。分別用質(zhì)量分?jǐn)?shù)為2%，4%，6%，8%的酸性介質(zhì)進(jìn)行應(yīng)力變化與腐蝕速率變化試驗(yàn)。同時(shí)為了得到應(yīng)力及應(yīng)變變化，在每個(gè)試件上布置兩個(gè)光柵光纖應(yīng)變傳感器。

試驗(yàn)裝置示意如圖2所示：試件放置在螺旋下壓裝置和試件板支撐裝置之間。采用積分法建立試件彎曲變形與不均勻應(yīng)力之間的關(guān)系。試驗(yàn)系統(tǒng)實(shí)物圖片見(jiàn)圖3。

圖2 試驗(yàn)裝置示意

圖3 應(yīng)力變化與腐蝕速率變化試驗(yàn)系統(tǒng)實(shí)物

圖4 介質(zhì)濃度為2%時(shí)，試件應(yīng)變-腐蝕速率的關(guān)系曲線

圖5 介質(zhì)濃度為4%時(shí)，試件應(yīng)變-腐蝕速率的關(guān)系曲線

圖6 介質(zhì)濃度為6%時(shí)，試件應(yīng)變-腐蝕速率的關(guān)系曲線

圖7 介質(zhì)濃度為8%時(shí)，試件應(yīng)變-腐蝕速率的關(guān)系曲線

2.1 試驗(yàn)結(jié)果分析

不同濃度介質(zhì)時(shí)，試件的應(yīng)變-腐蝕速率的關(guān)系曲線如圖47所示(圖中濃度均為質(zhì)量分?jǐn)?shù)，下同)，從圖中可以看出，隨著螺旋下壓裝置向下位移(應(yīng)力變大)的不斷變化，試件應(yīng)變?cè)谥饾u增大，此時(shí)腐蝕速率也呈現(xiàn)增大趨勢(shì)，但增長(zhǎng)率有所減小。當(dāng)向下位移增加到16 mm以后，試件所受應(yīng)變有所減小并趨于平緩，判斷此時(shí)試件發(fā)生了塑性變形，而相應(yīng)腐蝕速率也處于下降趨勢(shì)。而不同介質(zhì)濃度下，腐蝕速率隨著介質(zhì)濃度的增大呈現(xiàn)增大趨勢(shì)。由此可以看出，應(yīng)力的變化與腐蝕速率的變化存在正相關(guān)的關(guān)系。

2.2 不均勻應(yīng)力下腐蝕速率與聲發(fā)射特性分析

在腐蝕試驗(yàn)的過(guò)程中，利用美國(guó)物理聲學(xué)公司的PCI-8聲發(fā)射監(jiān)測(cè)系統(tǒng)對(duì)整個(gè)腐蝕速率變化過(guò)程進(jìn)行動(dòng)態(tài)監(jiān)測(cè)，傳感器選用DP3I低頻傳感器，監(jiān)測(cè)頻率為20 kHz100 kHz，在裝置外側(cè)兩面均布4個(gè)傳感器，試驗(yàn)裝置實(shí)物如圖8所示。

圖8 聲發(fā)射監(jiān)測(cè)試驗(yàn)裝置實(shí)物

不同濃度介質(zhì)時(shí)，試件的聲發(fā)射參數(shù)變化曲線如圖912所示(圖中ASL表示有效值電壓，單位V；幅值單位為dB;計(jì)數(shù)和能量無(wú)量綱)。

由圖912可以看出，在不同腐蝕介質(zhì)濃度下，隨著螺旋下壓裝置向下位移的增加，相應(yīng)濃度腐蝕介質(zhì)條件下聲發(fā)射的相關(guān)參量也隨之發(fā)生變化。聲發(fā)射特征參量撞擊數(shù)增大趨勢(shì)明顯，幅值、能量變化幅度較小，幅值在5070之間浮動(dòng)，聲發(fā)射計(jì)數(shù)和能量值變化浮動(dòng)較小，用于表征腐蝕速率隨不均勻沉降的變化規(guī)律的效果較差；觀察不同濃度條件下的聲發(fā)射參量數(shù)值可以發(fā)現(xiàn)，當(dāng)腐蝕介質(zhì)質(zhì)量分?jǐn)?shù)變?yōu)?%時(shí)，撞擊數(shù)相較于2%質(zhì)量分?jǐn)?shù)時(shí)有所增大(由24增加到42)，且隨著質(zhì)量分?jǐn)?shù)的繼續(xù)增大而增大，這說(shuō)明腐蝕速率的增大會(huì)使得聲源信號(hào)明顯增多，而聲源信號(hào)的增多在聲發(fā)射檢測(cè)中更加明顯地體現(xiàn)在了撞擊數(shù)的經(jīng)歷變化上。

圖13 聲發(fā)射參量與腐蝕速率的方差值圖

圖9 介質(zhì)濃度為2%時(shí)，試件聲發(fā)射參數(shù)變化曲線

圖10 介質(zhì)濃度為4%時(shí)，試件聲發(fā)射參數(shù)變化曲線

圖11 介質(zhì)濃度為6%時(shí)，試件聲發(fā)射參數(shù)變化曲線

圖12 介質(zhì)濃度為8%時(shí)，試件聲發(fā)射參數(shù)變化曲線

為了進(jìn)一步研究聲發(fā)射不同參量與腐蝕速率的關(guān)聯(lián)關(guān)系，將4種不同腐蝕介質(zhì)濃度下的聲發(fā)射參量的撞擊數(shù)、幅值、計(jì)數(shù)、能量、ASL值以及腐蝕速率求平均值并取方差，計(jì)算出不同濃度與下壓裝置不同向下位移下聲發(fā)射參量與腐蝕速率的均方根差值，如式(1)所示。

(1)

式中：x為聲發(fā)射參量或腐蝕速率；為聲發(fā)射參量或腐蝕速率的平均值；n表示數(shù)量。

聲發(fā)射參量的方差值如圖13所示。由圖13(a)可以看出，在撞擊數(shù)的均方根差值圖中，隨著向下位移的增大，方差不斷增大，說(shuō)明不均勻應(yīng)力的大小對(duì)聲發(fā)射撞擊數(shù)的變化會(huì)產(chǎn)生重要的影響，同時(shí)隨著濃度的升高，撞擊數(shù)的差值增大，說(shuō)明腐蝕速率的增大也會(huì)對(duì)撞擊數(shù)產(chǎn)生重要影響。由圖13(c)(f)可以看出，隨著下壓裝置向下位移的變化以及腐蝕介質(zhì)濃度的增加，幅值、計(jì)數(shù)、能量和ASL的方差變化幅度較小，這說(shuō)明不均勻應(yīng)力狀態(tài)下的聲源信號(hào)種類較為單一，基本聲源都屬于同類聲源，對(duì)于同類聲源利用幅值、能量等參數(shù)來(lái)分析演化規(guī)律的效果較差，由此在不均勻應(yīng)力條件下進(jìn)行聲發(fā)射檢測(cè)可以選擇撞擊數(shù)作為分析的主要參量，通過(guò)其時(shí)間經(jīng)歷變化情況來(lái)分析儲(chǔ)罐相應(yīng)時(shí)刻的腐蝕狀態(tài)。

3 結(jié)論

(1) 隨著不均勻沉降面積的增加，罐底板與罐壁板的應(yīng)力集中區(qū)域會(huì)出現(xiàn)在儲(chǔ)罐罐壁與底板交接部位以及儲(chǔ)罐沉降位置附近的邊緣板部位。

(2) 當(dāng)應(yīng)力恒定時(shí)，隨著腐蝕介質(zhì)濃度的增加，試件的腐蝕速率增大，且應(yīng)力越大，試件的腐蝕速率越大。

(3) 當(dāng)應(yīng)力恒定時(shí)，隨著腐蝕介質(zhì)濃度的增大，聲發(fā)射參數(shù)特征參量也隨之增加；當(dāng)腐蝕介質(zhì)濃度恒定時(shí)，試件的腐蝕速率隨著應(yīng)力的增大而增大，聲發(fā)射參數(shù)的撞擊數(shù)也隨之增大，當(dāng)試件的應(yīng)力較大時(shí)，試件腐蝕速率的變化率要比應(yīng)力較小時(shí)的變化率大，撞擊數(shù)的變化率也較大，進(jìn)一步驗(yàn)證了用撞擊數(shù)來(lái)判定不均勻沉降狀態(tài)下儲(chǔ)罐底板腐蝕情況的有效性。