負壓振動篩鉆井液流動及氣液分離數值模擬

負壓振動篩是一種新型的固液分離設備，對鉆屑表現出了出色的隨鉆干燥效果。然而，對于該設備的鉆井液流動規律、空氣和鉆井液透過篩網后的氣液分離機理以及效率等基礎理論方面，迄今尚未進行深入的研究。為此，采用計算流體力學（CFD）方法，建立了負壓振動篩的三維瞬態兩相流篩分仿真模型；通過歐拉-歐拉多相流模型（VOF）和多孔介質模型，系統研究了負壓作用時篩網上鉆井液的流動規律和篩網下的壓力分布特性。同時，探討了負壓對振動篩的處理量、鉆井液透篩速度以及氣液分離效率的影響。研究結果表明：在相同振動條件下，負壓振動篩的透篩速度提高了87%，處理量增加了18%；在負壓為-10 kPa時，分離效率達到95.8%。隨著負壓值的增大，氣液分離效率下降；排氣管道高度升高，氣液分離效率增大。所得結論可為負壓鉆井液振動篩的研發與應用提供理論指導。

負壓振動篩；鉆井液流動；氣液分離；CFD；VOF；多孔介質

中圖分類號：TE927

文獻標識碼：A

DOI： 10.16082/j.cnki.issn.1001-4578.2024.12.005

基金項目：四川省科技計劃重點研發項目“沙漠鉆井液隨鉆高效回收關鍵裝備研究及應用”（2022YFQ0064）。

Numerical Simulation on Drilling Fluid Flow and Gas-Liquid

Separation of Negative Pressure Shale Shaker

Hou Yongjun1" Du Lei1" Sun Lin2

（1.School of Mechanical Engineering，Southwest Petroleum University；2.No.3 Drilling Engineering Company，CNPC Bohai Drilling Engineering Company Limited）

Negative pressure shale shaker is a new type of solid-liquid separation equipment，and it exhibits excellent drying while drilling effect on drilling cuttings.However，there has been no in-depth research on the basic theoretical aspects such as the drilling fluid flow law of the equipment and the gas-liquid separation mechanism and efficiency after the air and drilling fluid have passed through the screen.In this paper，the computational fluid dynamics （CFD） method was used to build a 3D transient two-phase flow screening simulation model of negative pressure shale shaker.Then，the Euler-Euler multiphase flow model （VOF） and porous media model were used to systematically study the flow pattern of drilling fluid on the screen and the pressure distribution characteristics below the screen under the action of negative pressure，and explore the influence of negative pressure on the treatment capacity of shale shaker，drilling fluid sieving speed and gas-liquid separation efficiency.The results show that under the same shaking conditions，the sieving speed of the negative pressure shale shaker is increased by 87%，and the treatment capacity is increased by 18%.When the negative pressure is -10 kPa，the separation efficiency reaches 95.8%.As the negative pressure increases，the gas-liquid separation efficiency decreases.As the height of the gas exhaust duct increases，the gas-liquid separation efficiency increases.The conclusions provide theoretical guidance for the development and application of negative pressure shale shakers.

negative pressure shale shaker；drilling fluid flow；gas-liquid separation；CFD；VOF；porous media

0" 引" 言

傳統的鉆井振動篩在常壓下工作，完全依靠振動作用實現鉆井液的透篩和巖屑在篩面的運移。在篩網孔徑大于0.075 mm后，振動篩的處理量很小，排出的鉆屑含液量高，高含液巖屑等廢棄物的直接排放將對地下水和周邊生態環境造成嚴重污染［1-3］。通過增大篩分面積和提高振動篩振動強度，可使傳統振動篩的處理能力增大，但受固控系統罐面空間、篩網壽命等限制，處理能力的增大幅度有限，尤其是鉆屑干燥能力很難增強。為此，國內外研究者開展了雙頻振動篩分［4］、真空過濾篩分［5］、負壓振動篩分［6-7］等鉆井液固液分離方法的研究。

侯勇俊，等：負壓振動篩鉆井液流動及氣液分離數值模擬

M.J.DERRICK等［8］提出在傳統振動篩的篩網下交替施加負壓抽吸和正壓反吹的鉆井液振動篩分方法，但正壓反吹時，高壓氣流只能對局部篩網吹掃，無法實現清掃篩網的作用，反而使系統復雜化。而采用壓縮空氣噴射形成真空時，能耗大，形成大量液霧，尤其是使用油基鉆井液時，具有較大的安全隱患［9-13］。

侯勇俊等［12-13］提出了采用負壓風機抽吸形成篩網下的負壓，并在篩網下的負壓室中直接進行氣液分離的負壓振動篩；現場試驗和應用表明其處理量大、鉆屑干燥效果好、能耗低、無安全隱患?，F有報道中，以各種負壓振動篩分方法的現場試驗較多［14-15］。對負壓振動篩的基礎理論研究方面，侯勇俊等［16］基于傳統的單顆粒模型研究了負壓作用下篩網上的固相運移規律。李文霞［17］提出了鉆井液在負壓鉆井液振動篩上的流動的數學模型。YIN P.等［18-19］研究了脈動負壓作用下篩網上固相的運動規律。周思柱等［20］研究了正負交替的脈沖壓力作用下篩網上固相的運動特性。雷廳［21］對循環篩網負壓振動篩的鉆井液流動開展了數值模擬研究。

已開展的負壓振動篩的鉆井液流動數值模擬，將篩網下方視為一個壓力均勻分布的出口，透過篩網的氣液均從該出口流出，這與實際工況差異較大。筆者針對負壓氣液分離一體化的振動篩，建立了三維仿真模型，考慮負壓不均勻分布的情況，對鉆井液的流動狀態、負壓腔中的氣液分離效率等開展了數值模擬研究，以期為負壓振動篩的研發和推廣應用提供理論指導。

1" 負壓振動篩工作原理

負壓振動篩工作原理如圖1所示。負壓風機在篩網下的真空盤中形成一定的負壓，篩網上方的鉆井液與空氣的混合物在壓力的作用下快速透過篩網進入真空盤。進入真空盤中的混合物進行氣液分離，分離后的空氣從篩箱側板上的通道通過柔性軟管，進入液霧分離器后排入大氣（路徑如圖1中綠色箭頭所示）；分離后的鉆井液在重力作用下沉降到下部，通過自動排液裝置排出（路徑如圖1中藍色箭頭所示）。

2" 數值模擬

圖2為負壓振動篩流體域模型。該模型由5個流體域構成，分別是篩網上部的篩箱、篩網、篩網下部的常壓腔室和負壓腔室，以及負壓腔兩側的排氣管。根據相仿原理將篩網縮小至實際尺寸的1/4，其尺寸分別為：長500 mm、寬250 mm、高4 mm。

由于鉆井振動篩采用的篩網孔徑普遍達到0.075 mm及以上，鉆井篩網的鋼絲編織結構和幾何形狀復雜且數量龐大，要建立篩網的實際幾何模型并劃分網格難度較大。筆者采用了多孔介質模型模擬篩網，篩面上的巖屑對鉆井液的透篩流動具有一定的干涉作用，因此假設其在篩面上均勻分布，并將其對鉆井液的透篩流動的影響合并到篩網的多孔介質模型中考慮。

多孔介質模型實際上是一種分布阻力模型，通過在動量方程中引入阻力源項，使得流體在流經多孔介質區域時遇到流動阻力，從而產生壓力損失。根據Darcy定律［22］，多孔介質區域動量的源項Si為：

Si=-μαui+12C2 ρuui（1）

式中：Si為動量源項，Pa/m；α為滲透率，μm2;C2為慣性阻力系數，mm-1；μ為黏度，Pa·s；ui為速度，m/s；u為速度的模，m/s；ρ為密度，kg/m3。

當流動情況為層流時才會考慮滲透率的影響，根據經驗公式［23］可得：

C2=10.98（Ap/Af）2-1t（2）

式中：t為多孔板的厚度，mm;Ap/Af為多孔板的面積與孔的面積之比，無量綱。

選擇將鉆井液中的固相顆粒視為平鋪在篩網上的均布顆粒層。這種簡化方法避免了復雜的實際篩網模型和網格劃分的挑戰，同時可更專注地研究氣液兩相在篩網上的流動行為。顆粒層的阻力系數可以用半經驗公式確定，黏性阻力系數和慣性阻力系數分別為：

1α=150（1-ε）2ε3Dp2

C2=3.5（1-ε）ε3Dp（3）

式中：Dp為顆粒粒徑，取0.1 mm；ε為篩網的孔隙率，無量綱，取0.34。

2.1" 邊界條件和求解設置

如圖2所示，將計算模型的鉆井液入口定義為速度入口velocity-inlet，速度大小為1.9 m/s。模型的所有出口定義為壓力出口pressure-outlet，其中出口1、2和3的工作壓力（表壓）為0，出口4為負壓出口（0～3 s表壓為0，3 s之后為-10 kPa）。鉆井液在篩面上的流動為自由液面，采用VOF作為計算模型［24-25］。首相設置為空氣；次相設置為鉆井液，其密度為1 180 kg/m3，黏度為0.015 Pa·s。采用瞬態仿真，設置時間步長為1×10-4 s，收斂精度為1×10-5，壓力速度耦合算法為PISO，壓力離散格式為PRESTO！，動量、湍流動能以及湍流耗散率均采用二階迎風格式，體積分數采用幾何重構格式。

2.2" 篩網運動方程

振動篩的篩箱整體做往復運動，其運動方程為：

X=Asin（2πωt）

Y=Acos（2πωt）（4）

式中：A為篩網運動的振幅，取0.002 m；ω為篩網運動的角速度，取157 rad/s；t為運動時間，s。

2.3" 網格劃分與驗證

網格采用多面體網格，多面體網格通常具有更好的計算效率和精度［26］。面網格尺寸選擇4 mm。在多孔介質（篩網）附近的流體域和排氣管內的流體域，由于尺寸對流動的影響較大，并且流場也比較混亂，為了提高計算精度，對這2部分網格進行局部加密，面網格尺寸選擇1 mm。

網格的無關性驗證結果如表1所示。當網格數量超過535 284時，篩網上下面的壓降趨于一個穩定值，網格數量的增加不影響計算結果，因此后續的分析中采用網格數量為535 284。

李文霞［17］推導了鉆井液在篩面上的流動方程，并利用特征線法得到了鉆井液在篩面上各個位置的沿篩面流動的速度。這里通過仿真模擬得到的鉆井液沿篩面流動速度與文獻［17］具有一樣的趨勢。

3" 結果分析

3.1" 處理量

圖3為各出口的鉆井液流量情況。從圖3觀察到，在0～3 s階段，即振動篩的所有篩網均處于常壓工作狀態時，存在鉆井液“跑漿”的現象。這表明在這個時間段，進入振動篩的鉆井液流量超過了振動篩的處理能力。在3 s后，施加負壓，此時產生抽吸作用，導致出口2對應篩網的透篩量增加，“跑漿”現象消失。

在t=3.8 s時，篩網上的自由液面呈現明顯的“弓形效應”，篩網兩側壁面位置的鉆井液流程最長，但距離出口3還有一段距離，這確保了在振動速度周期性變化的過程中，振動篩不會出現“跑漿”的情況。在負壓狀態下，振動篩的處理量為入口進入的鉆井液流量，為4.5 L/s；而在常壓下，相同振動篩的極限處理量為3.8 L/s?？梢钥闯?，負壓作用后，振動篩的處理量增加了18%，這表明負壓對振動篩的性能有顯著的改善作用。

3.2" 鉆井液分布狀態

圖4為鉆井液在不同時刻沿篩面的自由液面分布圖。從圖4可以看到，隨著時間的推移，篩面上的鉆井液逐漸增多，并最終達到一個穩定的分布狀態。在3 s之前，即振動篩處于無負壓狀態，相當于傳統模式的振動篩。此時，鉆井液完全覆蓋了整張篩網。在3 s之后，由于負壓的作用，“跑漿”的鉆井液逐漸減少，最終在3.8 s時達到了穩定的分布狀態。

圖5為在1.4、1.8、2.8和3.8 s時刻的中間截面上，沿著篩面長度方向鉆井液層的厚度變化。從圖5可以發現：在篩網上方的流體域中，鉆井液從入口進入后，在重力和初始速度的作用下沖擊篩網，形成了一個鉆井液層較薄的區域；隨后，在其后方形成了一個鉆井液層較厚的區域。圖5的點A（140.1，21.4）、B（122.9，23.6）、C（126.0，26.8）、D（125.0，25.0）分別代表了相應時刻篩網上鉆井液厚度最厚的點。即使在同一位置，鉆井液層的厚度也隨著時間不斷變化。

圖6為負壓腔內鉆井液在垂直于x方向截面上的分布。在負壓腔內，鉆井液受到兩側負壓的抽吸，導致其向兩側聚集。然而，在排氣管道前的截面（x=300 mm和x=320 mm），負壓腔內鉆井液向兩側聚集較少;而在排氣管的中間截面（x=375 mm和x=385 mm），該聚集現象較為明顯。這一觀察表明負壓系統的設計和排氣管的位置對于負壓腔內鉆井液的分布產生了影響：在排氣管道前，由于可能存在較小的負壓，鉆井液在這個區域的聚集較為有限；在排氣管的中間位置，負壓更為明顯，造成鉆井液更集中地向兩側聚集。這些結果對于理解和優化負壓振動篩的性能具有重要意義。

圖7為鉆井液沿著篩面流動的速度曲線。從圖7可以看到：隨著沿著篩網的距離增加，在3 s內，即傳統振動篩上，鉆井液的流動速度逐漸減??；而在3.8 s時，即轉化為負壓模式后，鉆井液的流動速度先減小再增大。這是因為鉆井液抵達負壓區域附近后，篩網下方已經形成了低壓區，增大了這一區域的鉆井液流動速度，表明負壓系統對篩網上鉆井液的流動產生了顯著的影響。

3.3" 鉆井液透篩流動速度

為研究負壓對鉆井液透篩的影響，在中間面選擇了出口對應篩網上的監測點1（260，0，-125）、2（300，0，-125）、3（340，0，-125）。在有負壓作用和無負壓作用下，透篩速度（垂直方向）隨時間的變化關系如圖8所示。從圖8可見，鉆井液的透篩流動呈周期性變化。

在常壓時（見圖8a），3個監測點的最大透篩速度幾乎相同，為0.24 m/s（波峰值）。這是因為篩網的振動為平動，監測點的液面高度較小，鉆井液的透篩速度主要由篩網的振動決定，各點的自由液面高度差異對透篩速度影響很小。

在負壓作用下（見圖8b），鉆井液的最大透篩速度為0.45 m/s（波谷值），與無負壓作用相比提高了87%。此外，越靠近負壓腔室空氣出口的區域，鉆井液的透篩速度越大。值得注意的是，圖8顯示鉆井液的透篩速度出現了大于0的情況，這是因為鉆井液附著在篩網上，篩網的振動不足以克服這種黏附作用。當振動篩向上運動時，鉆井液具有正向速度，此時尚未透過篩網。特別是在常壓振動篩上，這種現象最為明顯，而負壓振動篩的負壓作用加大了克服黏附作用的能力，這種現象并不顯著。

3.4" 壓力分布

圖9為流場穩定后篩網下的壓力分布情況（t=5 s）。其中0～250 mm部分對應第一張篩網（下部為出口1），250～500 mm對應第二張篩網（下部為出口2）。當篩網向下運動時，篩網的速度大于鉆井液自由落體的速度，可能導致鉆井液和篩網之間的短暫分離。對于細目篩網或者致密的多孔介質，篩網下的液體可能暫時堵住篩網網孔，形成第一張篩網下的局部負壓，而第二張篩網下的負壓主要分布在篩網兩側。

圖10展示了第二張篩網下方的壓力隨時間的變化情況。在負壓系統剛開始作用的較短時間內，篩網下的靜壓力會發生一個突變，最大值可達到-635 Pa。隨后，負壓逐漸減少，并最終維持在-205 Pa左右。這是因為負壓作用后覆蓋在篩網上的鉆井液逐漸減少，空氣透過篩網的區域增大，所以該篩網下的負壓呈現這樣的變化規律。

3.5" 鉆井液分離效率

圖11為出口4在不同負壓下，液相體積分數所占的百分比。從圖11可以看到，鉆井液從出口4被吸出是隨機的，但隨著施加的負壓增大，出口4處鉆井液的含量會增加。在-10 kPa下，出口4處液相體積分數的平均值為4.2%（即分離效率為95.8%），而在-15 kPa時則增至9.0%。

圖12展示了不同壓力和排氣管道高度下的鉆井液分離效率。其中定義出口2的液相體積占出口2與4液相體積之和的百分比為鉆井液的分離效率。從圖12可以看到，隨著負壓的增大，鉆井液的分離效率降低。然而，隨著排氣管道高度的增加，鉆井液的分離效率卻提高。因此，在一定負壓下，通過適當增加排氣管道的高度可以提高鉆井液的分離效率。這些結果對于優化負壓振動篩的設計和操作參數具有指導意義。

4" 現場應用

圖13為2種條件下鉆井液振動篩應用效果。從圖13可以看到：沒有負壓作用時，篩面上排出的鉆屑中含有大量的鉆井液，并伴隨有部分“跑漿”現象；在有負壓作用時，篩面上排出的鉆屑較為干燥，鉆井液回收效果好?，F場應用的效果驗證了數值分析結果。

5" 結論及建議

（1）相較于常壓模式，負壓模式能夠有效防止鉆井液振動篩出現“跑漿”問題，并顯著提升振動篩對鉆井液的處理量，模型中負壓模式的處理量提高了18%。

（2）負壓振動篩的透篩速度呈現周期性變化，最大透篩速度比無負壓時顯著增加，增幅達87%，使振動篩具備更強的克服鉆井液黏附作用的能力。這對提高振動篩對排出鉆屑的干燥能力具有重要作用。

（3）篩網下方的壓力分布并不均勻。在開始施加負壓時，篩網下方觀測點的壓力會突然降低至-635 Pa而后隨著篩面上鉆井液逐漸透過篩面，負壓穩定在-205 Pa。

（4）系統的氣液分離效率會隨著負壓值的增大而減小，隨著排氣管高度的增加而增加。因此在負壓值一定時，可以考慮提高負壓系統的高度來增加氣液分離的效率。

（5）排氣管高度一定時，出口負壓在-10～-15 kPa時，排出的空氣中液相體積分數為4.2%～9.0%。建議在排出管出口設置氣霧分離裝置，尤其是使用油基鉆井液時，能減少出口外空氣中的油霧含量，確保井場安全。

［1］" 謝水祥，鄧皓，王蓉沙，等．鉆井環境污染過程控制技術綜述［J］．油氣田環境保護，2008，18（2）：38-40．

XIE S X，DENG H，WANG R S，et al.A review on environmental pollution process control technology in drilling［J］.Environmental Protection of Oil amp; Gas Fields，2008，18（2）：38-40.

［2］" 袁姍姍．鉆井生產過程中的環境污染及治理［J］．油氣田環境保護，2009，19（1）：47-50．

YUAN S S.Environmental pollution and control in drilling process［J］.Environmental Protection of Oil amp; Gas Fields，2009，19（1）：47-50.

［3］" 雷先革，牟長清，涂志威，等．鉆井固控環保一體化裝備的研制與應用［J］．石油機械，2017，45（12）：28-31．

LEI X G，MU C Q，TU Z W，et al.Development and application of integrated drilling solid control and environmental protection equipment［J］.China Petroleum Machinery，2017，45（12）：28-31.

［4］" 侯勇俊，方潘．一種雙頻振動篩：CN201120479461.3［P］．2011-11-28．

HOU Y J，FANG P.A double-frequency vibrating screen：CN201120479461.3［P］.2011-11-28.

［5］" COURTNEY D.Novel screening unit provides alternative to conventional shale shaker［J］.Oil amp; Gas Journal，1999，97（15）：40-48.

［6］" 史萬明，李志偉，萬國華，等．一種真空負壓固液分離振動篩：CN201520693912.1［P］．2015-09-09．

SHI W M，LI Z W，WAN G H，et al.A vacuum negative pressure solid-liquid separation shaker：CN201520693912.1［P］.2015-09-09.

［7］" 陳金瑞，馮群娣，周靜，等．真空吸附設備在油基鉆井液處理中的應用研究［J］．石油機械，2022，50（12）：46-50．

CHEN J R，FENG Q D，ZHOU J，et al.Application of vacuum adsorption equipment in treatment of oil-based drilling fluid［J］.China Petroleum Machinery，2022，50（12）：46-50.

［8］" DERRICK M J，DERRICK R G，MOONEY J A.Vibratory screening machine with suction and pressure and method for screening a slurry：US2003/0230541 Al［P］.2022-06-12.

［9］" POMERLEAU D G.Systems and methods for drying drill cuttings：US 9，359，840 B2［P］.2012-11-12.

［10］" POMERLEAU D G.Optimization of vacuum system and methods for drying drill cuttings：US 20150204151A1［P］.2015-03-27.

［11］" 侯勇俊，方潘，張麗萍，等．一種負壓鉆井液振動篩：CN201510618406.0［P］．2015-09-25．

HOU Y J，FANG P，ZHANG L P，et al.A negative pressure drilling fluid vibrating screen：CN201510618406.0［P］.2015-09-25.

［12］" 侯勇俊，王鈺文，方潘，等．負壓發生和氣液分離一體化的負壓鉆井振動篩：CN202110311176.9［P］．2021-03-24．

HOU Y J，WANG Y W，FANG P，et al.Negative pressure drilling vibrating screen integrated with negative pressure Generation and gas-liquid separation：CN202110311176.9［P］.2021-03-24.

［13］" 侯勇俊，吳先進，方潘，等．一種脈沖負壓鉆井振動篩：CN202111109944.9［P］．2021-09-23．

HOU Y J，WU X J，FANG P，et al.A pulse negative pressure drilling vibrating screen：CN202111109944.9［P］.2021-09-23.

［14］" 岳明，劉陽，李治衡，等．水基鉆井液鉆屑減量處理技術探討［J］．非常規油氣，2021，8（6）：1-6，39．

YUE M，LIU Y，LI Z H，et al.Discussion on cuttings reduction treatment technology of water-based drilling fluid［J］.Unconventional Oil amp; Gas，2021，8（6）：1-6，39.

［15］" 張羽臣，林家昱，謝濤，等．履帶式負壓振動篩在海洋平臺鉆井廢棄物處置中的應用［J］．能源化工，2021，42（1）：61-65．

ZHANG Y C，LIN J Y，XIE T，et al.Application of crawler negative pressure vibrating screen in offshore drilling waste disposal［J］.Energy Chemical Industry，2021，42（1）：61-65.

［16］" 侯勇俊，李文霞，吳先進，等．負壓鉆井液振動篩上固相顆粒運移規律研究［J］．工程設計學報，2018，25（4）：465-471．

HOU Y J，LI W X，WU X J，et al.Study on conveyance law of solid particle on negative pressure drilling fluid vibrating screen［J］.Journal of Engineering Design，2018，25（4）：465-471.

［17］" 李文霞．負壓振動篩篩分機理研究［D］．成都：西南石油大學，2018．

LI W X.Study on screening mechanism of negative pressure vibrating screen［D］.Chengdu：Southwest Petroleum University，2018.

［18］" YIN P，HOU Y J，WU X J.Simulation of particles screening in pulsating negative pressure shale shaker by coupling CFD and DEM［J］.Engineering Computations，2022，39（5）：1701-1722.

［19］" YIN P，HOU Y J.Dynamic study on particle throwing motion of pulsating negative pressure shale shaker［J］.Journal of Petroleum Exploration and Production Technology，2023，13（2）：753-762.

［20］" 周思柱，李旭東，呂志鵬，等．負壓鉆井振動篩的動力特性分析［J］．機械設計與制造，2018（10）：43-46．

ZHOU S Z，LI X D，L Z P，et al.Dynamic characteristics analysis of negative-pressure shale shaker［J］.Machinery Design amp; Manufacture，2018（10）：43-46.

［21］" 雷廳．循環篩網負壓振動篩鉆井液流動規律研究［D］．成都：西南石油大學，2018．

LEI T.Study on drilling fluid flow law of circulating screen negative pressure vibrating screen［D］.Chengdu：Southwest Petroleum University，2018.

［22］" MENTER F R.Two-equation eddy-viscosity turbulence models for engineering applications［J］.AIAA Journal，1994，32（8）：1598-1605.

［23］" SMITHJR P L，VAN WINKLE M.Discharge coefficients through perforated plates at reynolds numbers of 400 to 3，000［J］.AIChE Journal，1958，4（3）：266-268.

［24］" NICHOLS B D，HIRT C W.Improved free surface boundary conditions for numerical incompressible-flow calculations［J］.Journal of Computational Physics，1971，8（3）：434-448.

［25］" HIRT C W，NICHOLS B D.Volume of fluid （VOF）method for the dynamics of free boundaries［J］.Journal of Computational Physics，1981，39（1）：201-225.

［26］" WANG W，CAO Y，OKAZE T.Comparison of hexahedral，tetrahedral and polyhedral cells for reproducing the wind field around an isolated building by LES［J］.Building and Environment，2021，195：107717.

第一侯勇俊，教授，博士生導師，生于1967年，2002年畢業于西南石油大學機械設計及理論專業，現從事石油機械方面的研究工作。地址：（610500）四川省成都市。email：hyj2643446@126.com。

2024-06-23任" 武