矩陣?yán)碚撟鳂I(yè)

矩陣?yán)碚撛阢@柱力學(xué)分析方面的應(yīng)用摘要：鉆柱力學(xué)是井眼軌道設(shè)計和控制、鉆柱設(shè)計及鉆井參數(shù)優(yōu)選的基礎(chǔ)。本文主要從鉆柱力學(xué)與井眼軌跡控制出發(fā)，以彈塑性力學(xué)為基礎(chǔ)，對多穩(wěn)定器旋轉(zhuǎn)鉆下部鉆具、帶彎外殼下部鉆具以及導(dǎo)向鉆具力學(xué)性能進(jìn)行了力學(xué)分析，并結(jié)合矩陣?yán)碚摰姆椒▽栴}簡化，最后對方程求解，并對計算結(jié)果進(jìn)行了分析。引言：鉆柱力學(xué)是指應(yīng)用數(shù)學(xué)、力學(xué)等基礎(chǔ)理論和方法，結(jié)合實驗以及井場資料等數(shù)據(jù)綜合研究受井眼約束的鉆柱的力學(xué)行為的工程科學(xué)。開展鉆柱力學(xué)研究，對鉆柱進(jìn)行系統(tǒng)、全面、準(zhǔn)確的力學(xué)分析，在井眼軌道設(shè)計與控制、鉆柱強度校核、鉆柱結(jié)構(gòu)和鉆井參數(shù)優(yōu)化等方面都具有重要意義。1下部鉆具力學(xué)分析方法鉆柱力學(xué)研究從最初的解決防斜打直問題，發(fā)展到解決定向井軌跡控制問題，從一維、二維發(fā)展到三維，從靜態(tài)發(fā)展到動態(tài)。最終形成了集中比較典型的研究方法，即：微分方程法、能量法、有限差分法、縱橫彎曲連續(xù)梁法和有限元法。1.1微分方程法經(jīng)典微分方程法是鉆柱力學(xué)中應(yīng)用最早的研究方法。該方法要求在滿足經(jīng)典材料力學(xué)的基本假設(shè)的前提下，建立鉆柱線彈性的經(jīng)典微分方程并求解。這種方法在考慮因素較多時，建立的微分方程很復(fù)雜，用經(jīng)典微分方程法求解比較困難。1.2能量法能量法是一種求解簡單的彈性力學(xué)問題的方法。它要求勢能函數(shù)不僅要滿足彈性力學(xué)的控制方程，而且要滿足邊界條件，通過解的形式的假設(shè)及有關(guān)參數(shù)的確定，可得到問題的解答。由于滿足以上2個條件是一件非常困難的事情。因此，這一方法的應(yīng)用受到了限制。1.3有限差分法有限差分法是一種近似方法。是通過對鉆柱進(jìn)行力學(xué)分析得到鉆柱微分方程式，再通過適當(dāng)?shù)牟罘洲D(zhuǎn)換將位移控制方程轉(zhuǎn)化為差分的形式求解。由于差分方程的系數(shù)是可變的，因此可以很容易考慮非線性的影響;同時，由于差分區(qū)間可以減小，可以比較容易考慮井眼的約束。但是要得到精確的解，答，差分區(qū)間必須取得很小，這樣就使矩陣的維數(shù)增加，降低了計算速度。對于鉆柱力學(xué)來說，有限差分法是一種有效的近似計算方法。1.4縱橫彎曲連續(xù)梁法縱橫彎曲連續(xù)梁法是一種精確解法，這種方法是將鉆柱視為相互聯(lián)系的縱橫彎曲的連續(xù)梁，應(yīng)用材料力學(xué)中的三彎矩方程建立一組非線性代數(shù)方程，該方程物理概念清楚，計算簡單，且速度較快。由于這種方法是將三維空間問題分解成2個獨立的二維問題求解，力學(xué)模型簡化得太多，忽略了扭矩及可能的力和變形的耦合問題。這種方法在國內(nèi)得到了推廣和應(yīng)用。1.5有限元法有限元法是一種近似數(shù)值計算方法，這種方法是通過將鉆柱分解為有限的離散梁單元，再通過適當(dāng)?shù)暮铣煞椒▽⑦@些單元組合成一個整體，用以代表原來的鉆柱狀態(tài)，并最終得到一組以節(jié)點位移為未知量的代數(shù)方程組。有限元法的物理概念清楚、簡單，實用性強。不限制鉆柱的材料和幾何形狀，且對單元尺寸也無嚴(yán)格的要求;又可以較容易地考慮非線性的影響。目前發(fā)展的接觸有限元法，考慮了鉆柱、穩(wěn)定器與井壁之間的初始接觸摩擦力，力學(xué)模型比較準(zhǔn)確，考慮因素較多，解題的速度雖然是這幾種方法中最慢的，但也可滿足需要。2鉆柱力學(xué)研究的有限元法2.1基本假設(shè)①鉆柱為三維彈性梁，離散后同一單元的幾何尺寸和材料性質(zhì)為常數(shù)。②井壁是剛性的，井眼是圓的。③鉆具上端與下井壁相切。④穩(wěn)定器與井壁之間為點接觸。⑤鉆頭處為鉸支，鉆頭與地層之間無力矩作用。2.2力學(xué)模型的建立2.2.1BHA的離散化分成三維梁柱單元，一般1m長為一個單元，理論上單元越短，計算結(jié)果越精確，劃分的原則是，單元內(nèi)外徑和內(nèi)徑是相同的。2.2.2建立單元的力學(xué)平衡方程主要包括單元的廣義節(jié)點力、廣義節(jié)點位移、單元的剛度矩陣（彈性剛度矩陣和幾何剛度矩陣）將下部鉆具組合離散化成兩節(jié)點三維直梁單元，單元位移列矩陣為：φ單元載荷列矩陣為：R由于作用在鉆柱上的力除橫向載荷外，還承受很大的軸向力，因此，需要考慮梁單元的幾何剛度矩陣KGK其中：KK式中：KEKG2.2.3等效節(jié)點力①單元所收到的外力主要是鉆具本身的自重；②自重分解為qx、qy、qz；③把qx、qy、qz通過虛功原理等效移植到兩個端點上。qx的移植為：Qqy的移植為：Q等效節(jié)點力計算公式為：R'按與廣義節(jié)點位移相對應(yīng)的等效節(jié)點力：R'單元平衡方程為：K2.2.4彎接頭和彎外殼的處理對于鉆柱中包含有彎接頭、彎外殼的情形，每一單元可作坐標(biāo)平移：xy在過渡坐標(biāo)系內(nèi)，原本不在同一條直線上的鉆柱變?yōu)橹钡你@柱，只不過此時在彎接頭處必須附加一個側(cè)向力。在過渡坐標(biāo)系下的方程為：K在局部坐標(biāo)系下：K2.2.5坐標(biāo)變換由于每一個單元的井斜和方位不同，有必要把單個局部坐標(biāo)系（x，y，z）下的單元平衡方程轉(zhuǎn)換到整體坐標(biāo)系下。整體坐標(biāo)系的定義為大地坐標(biāo)，X、Y、Z分別為北坐標(biāo)、東坐標(biāo)的垂深坐標(biāo)變換矩陣T應(yīng)滿足：x坐標(biāo)變換矩陣T應(yīng)滿足：T式中：T為轉(zhuǎn)換矩陣；θ為井斜角；φ為方位角；在整體坐標(biāo)系下的單元廣義位移和廣義力為：δR在整體坐標(biāo)系下的單元剛度矩陣為：KBHA平衡方程為：KK2.2.6方程的求解鉆柱動力學(xué)控制方程有多種求解方法，包括整體剛度矩陣有限元法和節(jié)點迭代有限元法。整體剛度矩陣有限元法：根據(jù)單元矩陣集成剛度，并根據(jù)邊界約束條件修正剛度、質(zhì)量、阻尼、作用力矩陣，然后求解方程。整體剛度矩陣法的特點是流程簡單，邊界處理方便。但在解決非線性問題時，整體剛度矩陣法需要多次求解大型矩陣，運算時間較長，因此，整體剛度矩陣法適合處理單元較少的線性問題。3鉆柱力學(xué)研究中存在的主要問題⑴現(xiàn)有的研究工作大多著重于理論上定性分析，與實際的鉆井工況還存在一定差距，定量分析尚有欠缺。⑵鉆柱力學(xué)基礎(chǔ)研究已日趨成熟，但是關(guān)于實際應(yīng)用研究的探索還不夠。4鉆柱力學(xué)的未來發(fā)展趨勢未來的鉆井技術(shù)必然走向自動化、智能化，鉆柱力學(xué)研究是前提和基礎(chǔ)。與現(xiàn)代鉆井技術(shù)的發(fā)展相適應(yīng)，基礎(chǔ)理論研究、系統(tǒng)控制論理論與方法、動態(tài)特性測試技術(shù)及系統(tǒng)虛擬仿真技術(shù)將成為鉆柱力學(xué)未來發(fā)展的主要方向。參考文獻(xiàn)[1]李子豐.鉆柱力學(xué)研究現(xiàn)狀及進(jìn)展[J].石油鉆采工藝,2008,30（2）[2]胡以寶.基于實際井眼軌跡的鉆柱動力學(xué)特性有限元分析[D].上海：上海大學(xué)，2011