趙洪波，劉寶林，王建強，3，王 瑜

(1.中國地質(zhì)大學(xué)〈北京〉國土資源部深部地質(zhì)鉆探技術(shù)重點實驗室，北京100083;2.北京探礦工程研究所，北京100083;3.中煤科工集團西安研究院，陜西 西安710077)

在高溫高壓下，由于橡膠的耐用環(huán)境的局限性，渦輪鉆具成為了地質(zhì)超深鉆探下的井底鉆具重要研究方向之一。尤其在萬米深孔下，渦輪鉆具配合新型金剛石鉆頭成為地質(zhì)鉆探技術(shù)方案的重要方向。

作為渦輪鉆具的核心部件，適應(yīng)于地質(zhì)鉆探的渦輪葉片的型線設(shè)計對渦輪鉆具的性能影響很大。目前，葉片設(shè)計主要包括幾何法和參數(shù)法2種，如Bezier曲線法、B樣條曲線等?；赟OLIDWORKS、UG/GRIP、Pro/E以及BAZIER曲線等與參數(shù)化設(shè)計相結(jié)合，都能很好地對渦輪葉片的型線設(shè)計提供幫助［1］。作為水力性能的研究隨著計算機和新型測試技術(shù)的發(fā)展不斷深入，三元流理論的應(yīng)用將進一步提高國內(nèi)渦輪鉆具研究水平［2］。本文采用五次多項式法設(shè)計渦輪葉片，使用SOLIDWORKS三維建模，導(dǎo)入ANSYS CFD模塊，通過模擬優(yōu)化可以避免所設(shè)計的葉片產(chǎn)生脫流嚴重、流動損失大、葉柵效率不高等狀況。同時，對不同粘度下的渦輪葉片的速度場、壓力場分布進行了模擬，為渦輪鉆具水力性能的預(yù)測提供了依據(jù)。

1 渦輪鉆具渦輪葉片的型線設(shè)計

地質(zhì)鉆探上所用的渦輪鉆具的渦輪一般為軸流式渦輪，渦輪定、轉(zhuǎn)子葉片的截面形狀，在徑向方向上不發(fā)生改變。參考地質(zhì)鉆探用127 mm的渦輪鉆具的性能要求，將渦輪葉柵葉片幾何參數(shù)求解后，各參數(shù)如為:葉片高度H=10 mm;進口結(jié)構(gòu)角β1k=110°;葉片弦長b=12.25 mm;出口結(jié)構(gòu)角 β2k=50°;前緣小圓半徑r1=0.8 mm;安裝角 βm=60°;后緣小圓半徑r2=0.4 mm;葉片前錐角γ1=15°;葉型最大厚度dmax=1.2 mm;葉片后錐角度γ2=10°;最大厚度位置a=4.5 mm;葉片數(shù)z=32 mm;葉柵間距t=8.83 mm;葉柵相對節(jié)距t=0.7208。

筆者從葉片型線與幾何參數(shù)的關(guān)系出發(fā)來研究葉片型線的造型，認為葉片型線選取的原則之一是葉片型線應(yīng)具有連續(xù)曲率。五次多項式作為葉片壓力面和吸力面型線方程式，結(jié)合計算機輔助設(shè)計可以很好地完成葉片造型［3］。其求解思路為假設(shè)葉片的壓力面yp和吸力面ys型線分別為:

為了更好求解方程，首先要對特殊點進行假設(shè)，這些特殊點包括壓力面上第一點的坐標(xp1，yp1)、一階導(dǎo)數(shù)yp1'和二階導(dǎo)數(shù)yp1″，最后一點的坐標(xpn，ypn)、一階導(dǎo)數(shù)ypn'、二階導(dǎo)數(shù)ypn″。吸力面上第一點的坐標(xs1，ys1)、一階導(dǎo)數(shù)ys1'和二階導(dǎo)數(shù)ys1″，最后一點的坐標(xsn，ysn)、一階導(dǎo)數(shù)ysn'、二階導(dǎo)數(shù)ysn″。將以上參數(shù)代入式(1)和式(2)。通過mathematics軟件求解，即可確定葉片壓力面和吸力面型線方程為:

此外，前緣圓心及后緣圓心坐標分別為(0.8，4.415)和(9.6，9.5756)。

對設(shè)計出的葉片進行檢驗一些附加檢驗，修正后畫出的渦輪鉆具定、轉(zhuǎn)子葉片如圖1所示。

圖1 渦輪葉片二維圖

2 CFD計算模型邊界條件的建立及網(wǎng)格的劃分

2.1 建立跨葉片的CFD計算模型

常用的CFD分析軟件通常具有全面的幾何建模能力，直接建立點、線、面、體等。本文采用ANSYS程序中的 FLOTRAN CFD進行分析。通過SOLIDWORKS三維造型軟件建立實體模型，再導(dǎo)入ANSYS進行分析比直接在ANSYS上建模更為簡單。用SOLIDWORKS三維專業(yè)造型軟件創(chuàng)建了定、轉(zhuǎn)子葉片的三維實體模型如圖2所示。

圖2 葉片三維實體模型圖

模型導(dǎo)入之后，往往所導(dǎo)入的不是實體形狀，應(yīng)在ANSYS里面首先進行處理，形成體。創(chuàng)建單周期渦輪定、轉(zhuǎn)子流道計算模型，為了使流場計算能夠得到穩(wěn)定解，需要將定子的入口處和轉(zhuǎn)子的出口處分別向上和向下延伸一定的距離，延伸距離取為3倍長度較為合適［4］，如圖3所示。

圖3 跨葉片流道的CFD模型圖

2.2 邊界條件

在建立邊界條件時，首先要選擇控制方程。渦輪葉片在正常工作時，可以認為通過葉片的流場是連續(xù)的，且所涉及的場(比如速度場、壓力場)是可微的。因此，筆者將 Navier－Stokes方程選擇為CFD分析計算的控制方程。雷諾時均方程(即NS方程)作為流場平均變量的控制方程，其相關(guān)的模擬理論都是基于湍流模式理論［5］。使用雷諾時均方程分析內(nèi)流場應(yīng)滿足以下條件［6］。

式中:ui——速度分量(i=x，y，z);ρ——密度;p——靜壓力;Rij——雷諾應(yīng)力張量，其計算公式為:

式中:μi——湍流粘性系數(shù);ui'，uj'——脈動速度;k——湍流動能;δij——克羅內(nèi)克爾函數(shù);Sij——變形率張量，其計算公式如下:

在CFD求解過程中，首先要對邊界條件進行定義，主要包括以下幾個方面［7］:

(1)取Y軸負方向為轉(zhuǎn)子旋轉(zhuǎn)方向，涉及到的邊界條件主要包括:轉(zhuǎn)子周期邊界條件、無滑移壁面邊界條件、出口邊界條件等(參見圖3);

(2)在模擬過程中將定子進口邊界條件設(shè)定為速度進口，速度為軸向平均流速;轉(zhuǎn)子出口邊界條件設(shè)定為壓力出口，將轉(zhuǎn)子出口條件設(shè)定為一個標準大氣壓(即參考壓力為14.7 PSI);

(3)定、轉(zhuǎn)子的周期邊界定義為周期循環(huán)(periodic);

(4)渦輪內(nèi)流場計算時，往往涉及到定子和轉(zhuǎn)子之間的相互作用，這是因為在定子出口與轉(zhuǎn)子進口之間，定子不旋轉(zhuǎn)，而轉(zhuǎn)子繞軸旋轉(zhuǎn)。因此將定、轉(zhuǎn)子結(jié)合面定義為frozen rotor;

(5)定、轉(zhuǎn)子葉片壁面及定、轉(zhuǎn)子的前緣和后緣均定義為無滑移壁面邊界;

(6)采用湍流模型中的k－ε模型，湍流強度(turbulence intensity)定義為medium。

2.3 網(wǎng)格劃分

通過分析，可以將渦輪一組定、轉(zhuǎn)子葉片及葉柵間空間當做一個流道分析，并且各組定、轉(zhuǎn)子葉片的流道可看做是相同的，可以認為葉柵間的流動狀態(tài)為周期相似。FLOTRAN CFD分析功能在完成跨葉片流道模型后，根據(jù)以上分析，其網(wǎng)格劃分使用ANSYS FLOTRAN CFD中用于三維分析的FLUID 142單元，進而對處于流道之間的定、轉(zhuǎn)子葉片的流場空間進行劃分。在網(wǎng)格劃分時，應(yīng)進行合理的細劃、加密，以提高求解準確性。網(wǎng)格劃分后的三維模型如圖4所示。

圖4 利用FLUID 142單元進行的網(wǎng)格劃分圖

3 CFD求解過程及計算結(jié)果分析

3.1 FLOTRAN CFD 求解過程

采用ANSYS程序中的FLOTRAN CFD模塊對定、轉(zhuǎn)子跨葉片的流道進行分析求解的過程可以用圖5表示。

圖5 渦輪葉片CFD分析流程圖

渦輪鉆具用于地質(zhì)超深鉆探時，為了模擬出流體(鉆井液、泥漿)的粘度變化對渦輪葉片的輸出性能會產(chǎn)生一定的影響，對127 mm的軸流式渦輪鉆具進行了多種粘度分類模擬，在流量的計算時，考慮容積的效率損失，近似容積率取0.8，這里取流量為14 L/s。在超過3000 m的深部鉆探時，懸屑和護壁的視粘度可調(diào)范圍在5～35 mPa·s之間，本文分別取粘度系數(shù)為5、20、35 mPa·s三種情況進行模擬。

3.2 CFD求解結(jié)果分析

對不同情況進行模擬后得到的定、轉(zhuǎn)子表面速度場分布圖和壓力分布圖如圖6～8所示。

從不同的粘度下的模擬云圖可以得出以下結(jié)論。

(1)通過渦輪的流體進入定子之后，分為2部分，一部分沖擊定子吸力面，另一部分對壓力面造成沖擊;進入渦輪葉片定子后的速度變化不大，只是在定子后緣兩邊速度開始變大，在后緣的末尾處速度近似為零。流體流過定子后，進入轉(zhuǎn)子，流體同樣分成兩部分，通過壓力面的流體速度逐漸變大，并在后緣處形成最大，進入后一組定、轉(zhuǎn)子葉片;

(2)通過吸力面的流體速度低于壓力面，從壓力場分布中得知，壓力面壓力也明顯大于吸力面，可以看出輸出扭矩較好，減少了水力損失。從葉片模擬的整體情況上看，沒有出現(xiàn)脫流現(xiàn)象，說明葉型選擇合理，可以認為所設(shè)計的葉片較好地滿足了設(shè)計要求;

圖6 5 mPa·s時，跨葉片流道速度場與壓力場分布圖

圖7 20 mPa·s時，跨葉片流道速度場與壓力場分布圖

(3)結(jié)合圖6～8，可以看出當流量一定時，流經(jīng)渦輪葉片的流體(鉆井液、泥漿)的粘度越大，其出口速度越小，渦輪定子葉片入口處的壓力值越大，定、轉(zhuǎn)子之間的壓差較明顯。

4 結(jié)論

圖8 35 mPa·s時，跨葉片流道速度場與壓力場分布圖

(1)基于五次多項式來構(gòu)造渦輪定、轉(zhuǎn)子葉柵葉片型線可以很好地滿足地質(zhì)鉆探用渦輪鉆具的設(shè)計要求。對模擬結(jié)果分析，可快速對葉片進行優(yōu)化設(shè)計。形成了設(shè)計——模擬——優(yōu)化的一整套渦輪鉆具渦輪葉片設(shè)計方案;

(2)通過三維造型軟件繪制葉片三維模型，并導(dǎo)入ANSYS CFD分析模塊，建立三維的跨葉片流道模型，改變不同情況下的鉆井液粘度來模擬出葉片的速度場和壓力場變化，可以為渦輪鉆具的實際應(yīng)用提供參考依據(jù);

(3)通過模擬，為下一步渦輪葉片的水力性能預(yù)測提供了依據(jù)。

［1］孟軍強，鐘易成.基于NURBS的渦輪葉片設(shè)計及其性能分析［R］.江蘇南京:中國航空學(xué)會第六屆動力年會，2006:215－221.

［2］姚堅毅，劉寶林，王瑜.渦輪鉆具水力設(shè)計與分析方法應(yīng)用現(xiàn)狀研究［J］.石油礦場機械，2012，(3):4 －7.

［3］馮進，符達良.渦輪鉆具渦輪葉片造型設(shè)計新方法［J］.石油機械，2000，28(11):9 －12.

［4］萬邦列，李繼志.石油工程流體機械［M］.北京:石油工業(yè)出版社，1999.

［5］金忠青.N－S方程的數(shù)值解和紊流模型［M］.江蘇南京:河海大學(xué)出版社，1987.23 －66.

［6］張慢來，馮進，丁凌云，等.一種軸流式葉輪的全三維優(yōu)化設(shè)計［J］.長江大學(xué)學(xué)報(自科版)，2006，(6):83 －87.

［7］劉孝光，潘培道，胡昌軍.渦輪鉆具葉柵水力性能仿真優(yōu)化技術(shù)研究［J］.冶金設(shè)備，2007，(1):21－24.