王明波 王瑞和

中國石油大學（華東），青島，266580

0 引言

低滲油氣藏高效開發(fā)、原油采收率提高和低煤階煤層氣徑向水平井技術，是當前具有重要戰(zhàn)略意義的熱點課題，也是石油工程技術發(fā)展的重要方向。利用磨料射流高效破巖來鉆進徑向水平井是近年來逐漸成熟、具有鮮明技術特點的新方法。對于老井或近井地帶污染嚴重的油井來說，這一方法可以有效改善油藏的壓力分布，提高殘余油的采收率和單井的油氣產量；對于占我國煤層氣資源量比重40%以上的低煤階煤層氣資源來說，這一方法可以依靠多條水平井眼溝通煤層內的裂縫和裂隙，減小流體流動阻力，從而有效提高煤層氣開發(fā)效率。

目前，對磨料水射流中磨料顆粒受力的研究比較少。文獻［1］采用量綱比較和典型函數法對磨料顆粒受力進行了分析比較。作為典型的液固兩相湍流流動，流體湍流對磨料顆粒的拖曳力一直是學者們研究的重點，但流體湍流對磨料顆粒所受拖曳力的影響規(guī)律則一直沒有確定的結論。文獻［2－3］研究了湍流對顆粒受力的影響規(guī)律，但該實驗數據分散性較大，并且不同實驗的顆粒雷諾數（以顆粒直徑為長度特征尺度計算得到的雷諾數）Rep也不盡相同，實驗結果不具有可比性。

1 數值計算

計算是在恒定密度、等溫、不可壓縮和湍流流動條件下進行的。采用直接求解雷諾時均輸運方程的方法來處理湍流流動中的雷諾應力項。流體流動控制方程的具體表達式和各物理量的具體含義見文獻［4－5］。

計算域如圖1a所示。隨著Rep的增大，顆粒附近的流場不再呈現對稱性質。考慮到磨料顆粒附近流場的復雜性，將整個磨料顆粒置于流場中進行計算，這同文獻［6］中的選取有所不同，本文的方法能夠最大限度地減小流場尺寸對計算結果的影響?，F場常用的磨料顆粒粒徑為300～400μm，本次計算的磨料顆粒直徑dp＝350μm。現場常用的磨料顆粒體積分數在8%左右，在此體積分數范圍附近，磨料顆粒的間距較大，粒間碰撞幾率較小，采用單個磨料顆粒所得結論可以直接應用于現場實際。當Rep小于2500時，來流速度小于150m/s。顆粒表面邊界條件采用標準壁面函數處理。顆粒附近網格劃分情況參見圖1b。最小網格尺寸為50μm，最大網格尺寸為200μm。網格總數為1.54×106。為提高計算精度，采用具有3階精度的守恒律的單調迎風中心格式（monotone upstream －centered schemes for conservation laws，MUSCL）進行控制方程中動量方程與湍流方程的離散。程序收斂判據為所有控制方程的殘差絕對值的和小于10－6。

圖1 計算域與網格劃分情況

表征湍流的參數有很多，本文選取了湍流強度和湍流黏度比這兩個參數來開展研究。湍流強度I表征湍流脈動的激烈程度，是湍流脈動速度與平均速度的比值。湍流黏度比R是湍流黏度同流體分子黏度的比值，它同湍流雷諾數成正比。在高雷諾數邊界層流動、剪切層流動以及充分發(fā)展管流流動中，R可高達100甚至達到1000；在大多數外部流動的自由流動邊界上，R則較小，一般在1～10之間。本次計算中，入口平面上的湍流強度I在5%～70%之間變化，湍流黏度比R在1～1200之間變化。給定進口的湍流強度和湍流黏度比后，可根據湍流理論換算出該處的湍動能和湍動能耗散率，然后進行流場中湍流的模擬計算。磨料顆粒的受力信息通過積分顆粒表面的黏性力和壓力來獲取。

2 計算結果與討論

為了檢驗各種湍流模型對計算結果的影響，本文采用標準k－ε模型、RNGk－ε模型、Realizablek－ε模型和雷諾應力模型對受限湍流沖擊射流進行計算，并將計算結果與實驗結果［7］進行比較。圖2a給出了計算域，流體介質為FC－77流體［7］。圖2b為采用雷諾應力模型計算出的沖擊射流流場內的流線分布圖，從中可以得到再附點的位置。表1給出了采用不同湍流模型的計算結果與實驗結果的，其中，b＝6.35mm。

圖2 計算域及流場中的流線分布

表1 不同湍流模型計算結果與實驗數據的對比

由表1可知，雷諾應力模型的計算結果與實驗結果最接近。為此，本文選用雷諾應力模型對球形磨料顆粒附近的流動進行數值模擬，在Rep＜2500的條件下，研究了不同湍流強度和湍流黏度比對顆粒受力的影響。圖3給出了在Rep＝2100，R＝200條件下，湍流強度對顆粒受力的影響。隨著湍流強度的增大，顆粒所受的曳力也越來越大，相應的曳力系數CD也逐步增大。在時均速度一定的情況下，湍流強度增大意味著湍流脈動速度增大，即流場中湍流小渦團的比重增加。在顆粒表面粗糙度保持恒定的前提下，流經顆粒表面的湍流小渦團數目增加，湍流渦團耗散成機械能的數目也增大，對顆粒表面做功增加，造成的顆粒所受曳力增大，相應的其曳力系數也增大。文獻［2］采用實驗方法對顆粒所受曳力問題進行了研究，得到了湍流強度和積分長度尺度對顆粒所受曳力的影響規(guī)律。將本文計算結果同文獻［2］中的實驗結果對比，發(fā)現二者一致。

圖3 湍流強度對顆粒曳力的影響

湍流黏度比對顆粒所受曳力的影響見圖4。從圖4可以看出，隨著湍流黏度比的增大，顆粒所受曳力逐漸增大，相應的曳力系數也逐漸增大。若湍動能k和湍流黏度比R已知，則可計算得到湍動能耗散率ε

式中，Cμ為經驗常數，Cμ＝0.99。

圖4 湍流黏度比對顆粒曳力的影響

從式（5）可以看出，隨著湍流黏度比的增大，流場中的湍動能耗散率下降。湍動能耗散率是湍流中單位質量流體脈動動能的耗散率，即各向同性的小尺度渦的機械能轉化為熱能的速率。隨著小尺度渦所含機械能轉化為熱能速率的降低，含能大渦熱能轉化的速率也逐漸降低，含能大渦對固相顆粒做功逐漸增加，這表現為顆粒所受曳力逐漸增大，相應的其曳力系數也逐漸增大。

圖5 顆粒附近流體流速分布情況

為了分析湍流對顆粒受力的影響規(guī)律，圖5給出了圖1a中過顆粒球心向上直線上的速度分布。從圖5可知：對于湍流流動而言，在顆粒附近存在一個非常大的速度尖峰，然后流速快速降低并達到某一穩(wěn)定值。也就是說，在入口平面上的來流速度恒定的條件下，湍流黏性比的增大等同于湍流黏性的增大，而湍流黏性的增大勢必會影響到固相顆粒附近繞流情況的變化，這種變化所導致的固相顆粒近壁處較大的流速梯度變化勢必會對流體產生較大的阻力，根據作用力與反作用力的關系，固相顆粒所受曳力也勢必會增大。

3 結語

本文采用數值模擬的方法研究了磨料射流湍流流場中靜止磨料顆粒所受到的曳力受湍流參量的影響。數值計算中首先考察了不同湍流模型（標準k－ε模型、RNGk－ε模型、Realizablekε模型和雷諾應力模型）對淹沒沖擊射流流動計算結果的影響，從中選取表現比較優(yōu)越的雷諾應力模型對磨料顆粒（直徑dp＝350μm）附近的湍流流動進行了數值模擬。計算中，湍流強度在5%～70%之間，湍流黏度比為1～1200。由計算結果發(fā)現：在本文計算范圍內，隨著湍流強度的增強，固相顆粒所受曳力逐漸增大，顆粒的曳力系數增大；隨著湍流黏度比的增大，固相顆粒所受的曳力逐漸增大，相應的曳力系數也增大。

［1］王明波，王瑞和.磨料水射流中磨料顆粒的受力分析［J］.中國石油大學學報（自然科學版），2006，30（4）：47－49.

［2］Moradian N，Ting D S K，Cheng Shaohong.The Effects of Freestream Turbulence on the Drag Coefficient of a Sphere［J］.Experimental Thermal and Fluid Science，2009，33（3）：450－471.

［3］Warnica W D，Renksizbulut M，Strong A B.Drag Coefficients of Spherical Liquid Droplets，Part 2：Turbulent Gaseous Fields［J］.Experiments in Fluids，1995，18（4）：265－276.

［4］陶文銓.數值傳熱學［M］.2版.西安：西安交通大學出版社，2001.

［5］Hinze J O.Turbulence［M］.New York：McGraw－Hill Publishing Companies，1975.

［6］由長福，祁海鷹，徐旭常.湍流對氣固兩相流動中顆粒受力的影響［J］.清華大學學報（自然科學版），2002，42（10）：1357－1360.

［7］Fitzgerald J A，Garimella S V.A Study of the Flow Field of a Confined and Submerged Impinging Jet［J］.International Journal of Heat and Mass Transfer，1998，41（8/9）：1025－1034.