官振樂，張浩，李向東，孫偉，柳登瀚

（勝利油田技術(shù)檢測(cè)中心，山東東營 257061）①

在K 型井架現(xiàn)場(chǎng)檢測(cè)中，有時(shí)會(huì)發(fā)現(xiàn)新投產(chǎn)的井架出現(xiàn)承載能力達(dá)不到設(shè)計(jì)要求的情況，經(jīng)檢查發(fā)現(xiàn)井架底部鉸支不在同一水平線，井架出現(xiàn)左右傾斜導(dǎo)致井架的承載能力下降。這是由于在井架安裝時(shí)底部鉸支水平校正不合理或者地基出現(xiàn)沉陷等情況導(dǎo)致了井架整體出現(xiàn)了支座沉降缺陷。通過曲率模態(tài)差值對(duì)井架支座沉降缺陷進(jìn)行初步確定，然后建立有限元模型修正目標(biāo)函數(shù)，確定模型修正函數(shù)變量范圍，對(duì)井架支座沉降缺陷進(jìn)行定量分析。

1 基于曲率模態(tài)的損傷識(shí)別理論

根據(jù)材料力學(xué)，對(duì)直梁某截面處定義的彎曲靜力關(guān)系表示為［1］

式中：q（x）為直梁x 截面處桿件曲率；ρ（x）為直梁x截面處桿件曲率半徑；M（x）為直梁x 截面處所受彎矩；EI（x）為直梁x 截面處抗彎剛度。

由式（1）可知：當(dāng)直梁某截面出現(xiàn)損傷時(shí)，該截面的剛度就會(huì)變小，則該截面曲率增大，那么該截面損傷前后的曲率模態(tài)差振型曲線就會(huì)發(fā)生明顯的突變。根據(jù)該突變就可以識(shí)別出結(jié)構(gòu)的損傷位置。

在實(shí)際檢測(cè)中，曲率模態(tài)是不可以直接測(cè)量得出，但是可以通過加速度傳感器測(cè)出其位移模態(tài)，再由位移模態(tài)進(jìn)行中央差分計(jì)算得到［2］，即

式中：j為測(cè)點(diǎn)號(hào)；∮ij為結(jié)構(gòu)第i 階測(cè)點(diǎn)j 單元處位移模態(tài)；lj為測(cè)點(diǎn)之間的間距。

2 支座沉降缺陷對(duì)K 型井架的曲率模態(tài)的影響規(guī)律

2.1 有限元實(shí)現(xiàn)

在有限元模擬中，定義一個(gè)偏轉(zhuǎn)角度α，以井架繞旋轉(zhuǎn)軸左右旋轉(zhuǎn)α來實(shí)現(xiàn)支座沉降缺陷，設(shè)某井架立柱某一節(jié)點(diǎn)坐標(biāo)為（x1，y1），其繞z 軸偏轉(zhuǎn)，偏轉(zhuǎn)角度為α，該節(jié)點(diǎn)偏轉(zhuǎn)后的節(jié)點(diǎn)坐標(biāo)為（x2，y2），即

2.2 對(duì)井架立柱曲率模態(tài)的影響規(guī)律

以井架向左沉降0.2°為例，1號(hào)立柱上所有桿件節(jié)點(diǎn)的模態(tài)位移曲線如圖1，模態(tài)位移差曲線如圖2。

圖1 1號(hào)立柱節(jié)點(diǎn)模態(tài)位移曲線

圖2 1號(hào)立柱節(jié)點(diǎn)模態(tài)位移差曲線

由圖1～2可以看出：對(duì)于井架支座沉降缺陷，1號(hào)立柱節(jié)點(diǎn)模態(tài)位移曲線和模態(tài)位移差曲線都無法識(shí)別出來。因此利用井架的模態(tài)曲率差曲線來識(shí)別井架支座沉降缺陷，則1號(hào)立柱節(jié)點(diǎn)模態(tài)曲率差曲線如圖3。

圖3 1號(hào)立柱節(jié)點(diǎn)模態(tài)曲率差曲線

由圖3可以看出：井架1號(hào)立柱下段6個(gè)節(jié)點(diǎn)的模態(tài)曲率差曲線出現(xiàn)較大突變，而中段與上段節(jié)點(diǎn)的模態(tài)曲率差曲線基本沒有突變，其與局部缺陷的模態(tài)曲率差曲線的突變不同，因此可以使用模態(tài)曲率差方法來識(shí)別支座沉降缺陷。

在安裝井架動(dòng)態(tài)測(cè)試傳感器時(shí)應(yīng)著重于井架底部桿件位置，這樣采集出的數(shù)據(jù)才能識(shí)別出支座沉降缺陷。通過模擬仿真發(fā)現(xiàn)：當(dāng)支座沉降偏轉(zhuǎn)角度不斷加大時(shí)，圖3中出現(xiàn)明顯突變位置節(jié)點(diǎn)的模態(tài)曲率差值也隨之增大，因此可以通過曲率模態(tài)差值大小來確定井架支座沉降偏轉(zhuǎn)角度。

3 有支座沉降缺陷的K 型井架有限元模型修正

3.1 有限元模型修正函數(shù)

有限元模型修正是指利用現(xiàn)場(chǎng)實(shí)測(cè)的井架應(yīng)力數(shù)據(jù)，通過改變井架有限元仿真模型的物理參數(shù)和幾何參數(shù)，使修正后井架仿真模型的計(jì)算數(shù)據(jù)與實(shí)測(cè)應(yīng)力數(shù)據(jù)趨于一致。根據(jù)現(xiàn)行的井架檢測(cè)標(biāo)準(zhǔn)SY 6326—2012《石油鉆機(jī)和修井機(jī)井架底座承載

能力檢測(cè)評(píng)定方法及分級(jí)規(guī)范》，在進(jìn)行井架承載能力檢測(cè)時(shí)，一般只檢測(cè)井架3個(gè)不同斷面位置的桿件應(yīng)力［3－4］。因此基于實(shí)測(cè)應(yīng)力與ANSYS 優(yōu)化算法的模型修正方法［5］，利用井架實(shí)測(cè)應(yīng)力與待修正仿真井架對(duì)應(yīng)的單元桿件應(yīng)力構(gòu)建優(yōu)化目標(biāo)函數(shù)，即

式中：f（x）為目標(biāo)函數(shù)；n為提取的單元桿件數(shù)量；σi為待修正仿真井架單元桿件應(yīng)力；σr為實(shí)驗(yàn)井架單元桿件應(yīng)力；x為設(shè)計(jì)變量，考慮了井架缺陷位置、大小，其表達(dá)式如式（5）所示。

利用模型修正方法對(duì)以圖紙尺寸為依據(jù)建立的待修正仿真井架模型進(jìn)行目標(biāo)函數(shù)優(yōu)化。首先采用零階優(yōu)化尋找井架支座沉降缺陷設(shè)計(jì)變量的合理組合序列，再用一階優(yōu)化算法以此合理組合序列為起點(diǎn)進(jìn)行下一步優(yōu)化，確定目標(biāo)函數(shù)最優(yōu)時(shí)井架缺陷組合。

3.2 有限元優(yōu)化流程

有限元優(yōu)化流程如圖4。

圖4 有限元優(yōu)化流程

4 實(shí)驗(yàn)室井架模型驗(yàn)證

根據(jù)相似理論［6］，按照井架模型與原型結(jié)構(gòu)幾何尺寸比為1∶8制造了實(shí)驗(yàn)室井架模型。模型與原型材料的彈性模量比為1∶1，模型與原型的材料密度比為1∶1。

將實(shí)驗(yàn)室井架3號(hào)、4號(hào)立柱墊高11.5mm 來模擬實(shí)驗(yàn)室井架支座沉降缺陷，即實(shí)驗(yàn)室井架向1號(hào)、2號(hào)立柱沉降0.979°。在實(shí)驗(yàn)室井架上粘貼應(yīng)變片，總共在實(shí)驗(yàn)室井架上粘貼4層應(yīng)變片，每個(gè)層位有16個(gè)應(yīng)變片，采集井架在一定載荷下各測(cè)點(diǎn)的應(yīng)變數(shù)據(jù)；同時(shí)利用加速度傳感器采集井架的動(dòng)態(tài)數(shù)據(jù)，將含支座沉降缺陷實(shí)驗(yàn)室井架各測(cè)點(diǎn)位移模態(tài)與無缺陷實(shí)驗(yàn)室井架相應(yīng)測(cè)點(diǎn)位移模態(tài)做差，得到位移模態(tài)差曲線，如圖5所示。通過位移模態(tài)與曲率模態(tài)之間的關(guān)系，算出實(shí)驗(yàn)室井架各測(cè)點(diǎn)的曲率模態(tài)差曲線，如圖6所示。

圖5 節(jié)點(diǎn)模態(tài)位移差曲線

圖6 節(jié)點(diǎn)模態(tài)曲率差曲線

由圖5～6可知：實(shí)驗(yàn)室井架發(fā)生支座沉降缺陷，因此在模型修正中，確定模型修正函數(shù)變量x為支座沉降缺陷，支座沉降量未知。得到修正后模型井架出現(xiàn)如表1的損傷缺陷。

表1 損傷缺陷識(shí)別結(jié)果

取出該損傷缺陷下模型井架在14.4kN 載荷下應(yīng)變片測(cè)點(diǎn)A－A′、B－B′、C－C′、D－D′4個(gè)層位桿件處應(yīng)力，與實(shí)驗(yàn)室應(yīng)力測(cè)試得到的應(yīng)力相比較，結(jié)果圖7所示。

圖7 應(yīng)力比較曲線

由圖7可以看出：通過模型修正識(shí)別出缺陷的井架模型應(yīng)力與實(shí)驗(yàn)室井架實(shí)測(cè)應(yīng)力相差很小，相差最大值在4.41%，因此可以認(rèn)為該井架模型受力與實(shí)驗(yàn)室井架模型受力相同。

5 結(jié)論

1）通過位移模態(tài)差曲線來對(duì)損傷缺陷進(jìn)行識(shí)別，發(fā)現(xiàn)該曲線對(duì)損傷識(shí)別效果較不明顯。再利用由位移模態(tài)中央差分得到的曲率模態(tài)差曲線對(duì)損傷缺陷進(jìn)行識(shí)別，能較好識(shí)別出損傷位置，又能通過曲率模態(tài)差值對(duì)損傷程度進(jìn)行定量分析。

2）根據(jù)相似理論，在實(shí)驗(yàn)室建立JJ450／45－K型井架模型，分別采集出無缺陷井架、含立柱支座沉降缺陷井架的應(yīng)力數(shù)據(jù)和模態(tài)數(shù)據(jù)，通過曲率模態(tài)差值對(duì)井架的缺陷進(jìn)行定位識(shí)別，再結(jié)合模型修正理論與有限元優(yōu)化，對(duì)井架支座沉降缺陷進(jìn)行定量分析。修正得到的有限元模型與實(shí)驗(yàn)室井架采集的該載荷下的應(yīng)力的最大相對(duì)誤差為4.41%，因此可以認(rèn)為修正得到的有限元模型與實(shí)驗(yàn)室井架是等效的，基于修正的井架模型就可為井架承載能力做出更為精確的評(píng)估。

［1］鄒龍慶，遲利源，郭鳳，等.井架結(jié)構(gòu)損傷動(dòng)力學(xué)識(shí)別［J］.大慶石油學(xué)院學(xué)報(bào)，2006，30（3）：47－50.

［2］李忠忠.基于曲率模態(tài)方法的鋼桁架橋梁損傷識(shí)別［J］.工程建設(shè)與設(shè)計(jì)，2012，40（1）：147－179.

［3］SY 6326—2012.石油鉆機(jī)和修井機(jī)井架、底座承載能力檢測(cè)評(píng)定方法及分級(jí)規(guī)范［S］.

［4］張學(xué)軍，張樹珺，陳孝珍.JJ22547型井架有限元靜動(dòng)力特性分析［J］.石油礦場(chǎng)機(jī)械，2008，37（9）：49－51.

［5］李夯，齊明俠.基于ANSYS的K 型井架結(jié)構(gòu)分析［J］.石油礦場(chǎng)機(jī)械，2008，37（2）：5－7.

［6］周國強(qiáng)，趙慶梅，韓東穎.實(shí)驗(yàn)室井架模型的動(dòng)態(tài)測(cè)試研究［J］.石油鉆采工藝，2006，28（5）：10－12.