環空渦動流場特性分析范文

本站小編為你精心準備了環空渦動流場特性分析參考范文，愿這些范文能點燃您思維的火花，激發您的寫作靈感。歡迎深入閱讀并收藏。

《石油鉆采工藝雜志》2016年第二期

摘要：

為了深入研究環空渦動時的流場特性，根據流體力學理論，以連續性方程和N-S方程為控制方程，利用計算流體力學技術對鉆柱渦動時環空赫巴流體的流動進行了系統的數值模擬，研究了不同公轉方向下鉆柱偏心度的變化對環空切向速度剖面與合速度剖面的影響。通過對模擬數據的對比分析，發現不同公轉方向下環空流場分布截然不同，正向公轉時切向速度在環空寬間隙處隨著公轉速度、自轉速度和偏心度增大正向增大；反向公轉時會出現二次流，切向速度在環空寬間隙處隨偏心度的減小整體反向減小，同時二次流趨勢越明顯，摩阻壓耗越大。合理應用這些規律有助于完善現有鉆井水力學理論，更好地揭示井下環空流場特性，并為鉆井水力參數優化設計提供有效的理論指導。

關鍵詞：

渦動；環空；赫巴流體；摩阻壓耗；數值模擬；流場

鉆進過程中鉆柱在井筒內一般會出現渦動，即在地面轉盤帶動下順時針自轉的同時還繞井眼軸線公轉。渦動的出現以及公轉的半徑（即渦動時鉆柱的偏心度）和方向（順時針方向、逆時針方向）就目前鉆井水平而言無法控制且很難預測，通過實驗難以揭示其內部規律，最可行的方法是借助數值模擬技術進行研究［1-5］。因此研究渦動時的偏心度對由自轉旋流場、公轉旋流場和軸向流場疊加而成的環空真實流場［6-12］的影響有助于完善現有鉆井水力學理論，更好地揭示井下環空流場特性，并為鉆井水力參數優化設計提供理論指導［13-16］。

1環空流動的建模與模擬

1.1鉆井液模型為更好地描述非牛頓不可壓縮流體的流變特性，采用了更為合理的赫巴流變模式，該模式也可看作修正的冪律模式。

1.2幾何模型數值模擬所涉及到的幾何尺寸均采用現場實際尺寸，鉆柱外徑139.7mm，套管內徑244.0mm，環空長度15m，其他參數見表1。選取15m長的環空鉆井液作為研究對象，采用基于交錯網格的SIMPLEC算法對連續性方程和N-S方程進行離散化，動量方程的對流項進行離散化時為保證計算精度，采用二階格式；建立的網格幾何模型如圖1。設置3類邊界類型：速度入口邊界，壓力出口邊界以及旋轉壁面邊界，設定該邊界繞其自身軸線自轉的同時繞井壁軸線公轉，井壁為靜止壁面，所有壁面都簡化為無滑移水力光滑的壁面，且近壁面處的網格加密劃分，采用柱坐標系，對稱軸為z軸，按照速度分量方式設定速度入口值［20］。

1.3流動模型將流體視為不可壓縮流體，且具有一定的黏度值和密度值，不考慮傳熱和重力的影響［21］。用連續性方程和N-S方程共同描述環空流動，基于有限體積法的思想進行離散。由于渦動時的環空流場屬于復雜的旋流場，流體質點脈動程度較高，因此采用湍流兩方程模型中的Realizable模型，該模型能夠有效地將湍動黏度模型中的系數與應變率耦合起來，能夠更為精確地描述包括旋轉均勻剪切流、包含有射流和混合流的自由流、管道內流動、邊界層流動以及帶有分離的流動。因此文中采用Realizable模型［22］。

1.4環空流動的數值模擬渦動時自轉和公轉坐標系都遵循右手準則，規定轉盤自轉為順時針方向，以此作為旋轉正方向，鉆柱在井眼內的公轉根據已有研究的描述有順時針和逆時針2個方向，當公轉為順時針時，與自轉方向一致，根據右手準則，為正向渦動，反之為反向渦動。為更好地研究不同偏心度下流場的變化規律，設置其余鉆井參數在合理范圍，如表2。取15m長偏心環空鉆井液距入口9m處橫截面上沿對稱軸的剖面進行分析，橫坐標0表示井筒中軸線。

2鉆柱渦動時環空流場特性分析

設井筒圓心為O1，鉆柱圓心為O2，井筒半徑為R，鉆柱半徑為r，偏心距為h，公轉速度為ω1，自轉速度為ω2，P1為最寬間隙處對應的鉆柱外壁面位置，P2為最窄間隙處對應的鉆柱外壁面位置。設圖2中時刻的自轉線速度方向為正方向，根據分析可知：

2.1偏心度對流速剖面的影響分析模擬數據發現，其他條件不變時，隨軸向流速增大，環空內寬窄間隙處合速度剖面變化規律與一般環空流動相同，由于徑向速度較小且變化不明顯，本研究未作討論。下面著重分析渦動時的偏心度及公轉正反方向的變化對環空氣場的影響。設定環空鉆井液密度為1.0g/cm3，軸向流速為2.0m/s，公轉速度為18.84rad/s，自轉速度為18.84rad/s，考查偏心度分別為0.32、0.64、0.95時各流動參數剖面的變化規律。切向速度剖面：偏心度增大，寬間隙處切向速度剖面變大變寬，窄間隙處的切向速度剖面整體增大，鉆柱外壁面在寬間隙處根據式（3）的推導切向速度越小，在窄間隙處根據式（4）的推導切向速度越大，如圖3a。合速度剖面：偏心度增大，寬間隙處合速度剖面整體增大，窄間隙處合速度剖面整體減小，鉆柱外壁面在寬窄間隙處的變化規律同切向速度的變化規律，如圖3b。設定環空鉆井液密度為1000kg/m3，軸向流速為2.0m/s，公轉速度為–18.84rad/s，自轉速度為18.84rad/s，考查偏心度分別為0.32、0.64、0.95時各流動參數剖面的變化規律。切向速度剖面：隨偏心度增大，寬間隙處的切向速度剖面整體反向增大，窄間隙處的切向速度剖面先增大后減小，鉆柱外壁面在寬間隙處根據式（5）的推導切向速度正向增大，在窄間隙處根據式（6）的推導切向速度減小，井壁處切向速度為0，如圖3c。合速度剖面：隨偏心度增大，寬間隙處合速度剖面整體增大，窄間隙處的合速度剖面整體減小，鉆柱外壁面寬窄間隙處的變化規律同切向速度的變化規律，如圖3d。

2.2二次流的分析反向公轉時，寬間隙鉆柱外壁面附近的流體質點跟隨鉆柱自轉呈正向旋轉，此處的切向速度為正值，但在遠離鉆柱外壁面的區域，由于鉆柱正向自轉的趨勢弱于反向公轉的趨勢，故此區域的流體質點跟隨鉆柱反向公轉，所以此處切向速度為負值，如圖3c所示。在寬間隙區域，同時出現正反2個方向的切向速度，兩股不同方向的旋流在環空窄間隙的兩側分離與匯合，使得環空寬間隙處出現一個明顯的大渦，環空窄間隙處由于鉆柱旋轉壁面帶動的旋流將本應對應出現的另一個不對稱的小渦淹沒，最終形成了不對稱的二次流，如圖4所示。反向公轉時，其他參數不變，偏心度越大，越不容易出現二次流。寬間隙處二次流產生的大渦越靠近旋轉鉆柱外壁面，同時二次流的分離與匯合點越靠近窄間隙處，當偏心度大到一定程度后，鉆柱的正向自轉趨勢被二次流在寬窄間隙處的反向流動趨勢完全包裹且基本抵消掉，此時狀態與偏心環空自轉的螺旋流場（自轉方向為反向，自轉速度大小同公轉最為接近。二次流的明顯出現使得環空流體質點的脈動趨勢進一步加強，隨著偏心度減小，鉆柱正向自轉趨勢越不容易被抵消，二次流趨勢越明顯，偏心環空流動的摩阻壓耗將進一步增大。

3結論

（1）偏心度對渦動流場的影響較為復雜，不能輕易忽略。宏觀上不論公轉方向如何，隨偏心度增大，寬間隙處的合速度剖面整體增大，窄間隙處的合速度剖面整體減小，但當具體分析不同公轉方向的環空流場時，又有各自截然不同的流動特性。（2）正向公轉時，鉆柱的自轉和公轉共同增強了寬窄間隙處的正向旋流趨勢，此時切向速度整體為正值，無二次流出現，偏心度增大，寬間隙處切向速度增大。（3）反向公轉時，切向速度同時出現正反向旋流趨勢，鉆柱的公轉流場與自轉流場出現相互抵消狀態，環空截面呈現不對稱的二次流。偏心度減小，寬間隙處切向速度反向越小，公轉流場與自轉流場相互抵消越小，二次流趨勢越明顯，摩阻壓耗越大。

作者：張晉凱 李根生 黃中偉 田守嶒 宋先知 王海柱 單位：油氣資源與探測國家重點實驗室 中國石化石油工程技術研究院