巫光勝 ，張濤 ，錢(qián)真 ，韓琦

（1.中國(guó)石化西北油田分公司石油工程技術(shù)研究院，新疆 烏魯木齊 830011；2.西南石油大學(xué)石油與天然氣工程學(xué)院，四川 成都 610500）

0 引言

碳酸鹽巖縫洞型油藏在注水開(kāi)發(fā)過(guò)程中，注入水極易沿優(yōu)勢(shì)通道快速突破，導(dǎo)致注水效果不理想。通過(guò)向地層中注入調(diào)流劑顆粒，顆粒堆積可以有效封堵目標(biāo)縫洞，實(shí)現(xiàn)調(diào)整注水流動(dòng)通道的目的［1］。目前，固相顆粒在裂縫中的運(yùn)移行為在以水力壓裂中支撐劑輸送為代表的領(lǐng)域得到了廣泛的研究［2］，國(guó)內(nèi)外主要從室內(nèi)實(shí)驗(yàn)和數(shù)值模擬2個(gè)方面對(duì)顆粒在裂縫中的輸送和沉降行為進(jìn)行研究。Kern等［3］研究了不同攜帶液中的支撐劑的輸送，并提出了沉降式布砂理論和砂堤平衡高度的概念。Liu等［4-5］應(yīng)用平板裂縫模型，研究了液體流變性、顆粒濃度、裂縫寬度與粗糙度對(duì)低密度顆粒沉降規(guī)律的影響。溫慶志等［6］研究了高砂比（8%）條件下，中高密度（2 550～3 440 kg/m3）顆粒在大型平板裂縫模型中的鋪置規(guī)律。上述研究均為中高密度顆粒、高砂比的實(shí)驗(yàn)條件，對(duì)于低濃度、低密度顆粒在裂縫中的輸送行為則鮮有研究。

除實(shí)驗(yàn)外，數(shù)值模擬也廣泛運(yùn)用于固液兩相流的研究，主要包括歐拉-歐拉方法（雙流體模型）和歐拉-拉格朗日方法（離散相顆粒模型）［7-9］。歐拉-歐拉方法適用于高濃度顆粒的情況，將顆粒相看作擬流體，需要對(duì)顆粒相構(gòu)建復(fù)雜的本構(gòu)方程。當(dāng)顆粒濃度較低，可忽略粒間相互作用時(shí)，基于拉格朗日坐標(biāo)系下的CFD-DPM模型在追蹤顆粒軌跡方面表現(xiàn)出了計(jì)算量小、高效實(shí)用的特點(diǎn)［10-12］。本文通過(guò)大型平板實(shí)驗(yàn)展示了調(diào)流劑顆粒在裂縫中的沉降鋪置特征，綜合考慮雙向耦合與湍流效應(yīng)，建立了歐拉-拉格朗日坐標(biāo)系的CFD-DPM模型，模擬不同注入速度、顆粒直徑、攜帶液密度及黏度等參數(shù)下裂縫中的顆粒輸送行為，分析了調(diào)流劑顆粒在裂縫中的運(yùn)移特征及沉降百分比隨時(shí)間的變化規(guī)律。

1 調(diào)流劑顆粒輸送實(shí)驗(yàn)

1.1 實(shí)驗(yàn)內(nèi)容

大型平板裂縫可視實(shí)驗(yàn)裝置由6個(gè)部分組成，如圖1所示。

圖1 大型平板裂縫可視實(shí)驗(yàn)系統(tǒng)

裂縫模型長(zhǎng)×高×寬為4 m×60 cm×6 mm，入口為6個(gè)直徑15 mm的射孔眼，位于模型右側(cè)，出口位于模型左側(cè)頂部，直徑20 mm。

在裂縫模型中泵入一定量清水，使其循環(huán)并達(dá)到穩(wěn)定，再將均勻混合的漿體（調(diào)流劑顆粒和清水的混合物）泵入縫中；當(dāng)顆粒流經(jīng)觀察區(qū)域時(shí)，通過(guò)高清視頻采集系統(tǒng)獲取動(dòng)態(tài)資料，后期經(jīng)圖像分析，觀測(cè)顆粒的輸送與沉降狀態(tài)。需要說(shuō)明的是，當(dāng)排量較小時(shí)，大量顆粒會(huì)在泵中沉積導(dǎo)致卡泵，故實(shí)驗(yàn)中先采用高排量注入，再逐漸將排量減小，直至完成所有實(shí)驗(yàn)。

針對(duì)調(diào)流施工現(xiàn)場(chǎng)的待優(yōu)化施工參數(shù)，開(kāi)展敏感性分析。參數(shù)取值如表1所示。

1.2 結(jié)果分析

隨顆粒的注入，顆粒不斷發(fā)生沉降，并在裂縫底部形成砂丘。從底部進(jìn)入的顆粒和直徑較大的顆粒在重力作用下迅速沉降至底部并向前滑動(dòng)，當(dāng)動(dòng)能完全耗盡后，則在底部停止運(yùn)動(dòng)，其后進(jìn)入的顆粒則在其上堆積形成砂堤。隨注入時(shí)間的增加，砂堤高度逐漸增加，直到縫內(nèi)顆粒流動(dòng)達(dá)到平衡，此時(shí)砂堤高度不再發(fā)生變化。從進(jìn)口中上部注入的小顆粒則懸浮于流場(chǎng)中，最后流出裂縫。

表1 輸送實(shí)驗(yàn)參數(shù)

1.2.1 攜帶液密度的影響

圖2描述了在不同攜帶液密度條件下，顆粒沉降體積與縫內(nèi)流速的關(guān)系。由圖2可以看出，攜帶液密度越大，即與顆粒密度差越小，則顆粒沉降體積越小，顆粒懸浮能力越好。流速較小時(shí)，攜帶液密度從0.998 g/cm3增大到1.140 g/cm3，最終顆粒沉降體積明顯減少；流速較大時(shí)，攜帶液密度的影響減小，縫內(nèi)流速為0.88 m/s時(shí)，3種攜帶液密度統(tǒng)計(jì)得到的顆粒沉降體積相差不大，而流速較小時(shí)統(tǒng)計(jì)的沉降體積相差明顯。

圖2 不同攜帶液密度下沉降體積與縫內(nèi)流速的關(guān)系

1.2.2 攜帶液黏度的影響

圖3描述了不同攜帶液黏度條件下顆粒沉降體積與縫內(nèi)流速的關(guān)系。由圖3可以看出，攜帶液黏度越大，其攜帶性能越好，顆粒沉降體積越小，其結(jié)果與攜帶液密度影響結(jié)果類似。不同于圖2的是，不同流速下的沉降體積變化差別不大。

圖3 不同攜帶液黏度下沉降體積與縫內(nèi)流速的關(guān)系

2 液-固兩相流數(shù)值模型

2.1 數(shù)值模型

液相的質(zhì)量和動(dòng)量守恒方程分別為

式中：ρ為液相密度，kg/m3；t為時(shí)間，s；u 為液相流動(dòng)速度，m/s；p 為流場(chǎng)壓力，Pa；g 為重力加速度，m/s2；SM為液相與固相的動(dòng)量交換源項(xiàng)，N/m3；τ為液相黏性應(yīng)力張量，Pa；μ 為液相動(dòng)力黏度，Pa·s。

顆粒的牛頓第二運(yùn)動(dòng)定律為

式中：up為顆粒速度，m/s；FD，F(xiàn)vm，F(xiàn)p分別為液相阻力、虛擬質(zhì)量力、 壓力梯度力，N；ρp為顆粒密度，kg/m3；dp為顆粒直徑，m；CD為曳力系數(shù)；Res為顆粒雷諾數(shù)；a1，a2，a3為球形顆粒常數(shù)［13］。

2.2 物理模型

建立一個(gè)與裂縫模型大小相同的模型，在模型左側(cè)設(shè)置6個(gè)16 mm×6 mm的進(jìn)口，出口設(shè)置在模型右側(cè)頂部，長(zhǎng)×寬為20 mm×6 mm。為準(zhǔn)確模擬模型中的流場(chǎng)，采用正六面體結(jié)構(gòu)化網(wǎng)格對(duì)模型進(jìn)行網(wǎng)格劃分，網(wǎng)格總數(shù)目為129 600個(gè)：長(zhǎng)度方向上的網(wǎng)格數(shù)目為160個(gè)，寬度方向上為6個(gè)；在進(jìn)出口處對(duì)網(wǎng)格進(jìn)行加密處理，入口和出口的高度方向上的網(wǎng)格數(shù)目各7個(gè)，入口間壁面在高度方向上為10個(gè)。

進(jìn)口條件給定進(jìn)口速度及顆粒體積分?jǐn)?shù)，出口條件將標(biāo)準(zhǔn)大氣壓設(shè)為模擬壓力。液相在固體壁面上采用無(wú)滑移邊界條件；對(duì)于固相，底部邊界條件設(shè)置為“捕獲”，即顆粒運(yùn)動(dòng)到該邊界即被“吸收”。除底部以外的壁面邊界設(shè)置為“反彈”，其法向反射系數(shù)en和切向反射系數(shù)et計(jì)算公式分別為

式中：α為顆粒表面的入射角，（°）。

用有限體積法對(duì)流體相控制方程進(jìn)行離散，用相耦合的SIMPLE算法進(jìn)行迭代求解，收斂標(biāo)準(zhǔn)為各項(xiàng)殘差小于10-4，時(shí)間步長(zhǎng)0.002 s。采用CFD-DPM模型分析不同因素對(duì)顆粒的運(yùn)輸和沉降影響規(guī)律。各參數(shù)設(shè)置如表2所示。

表2 模擬實(shí)驗(yàn)參數(shù)

CFD-DPM模型中的流量Qm計(jì)算公式為

式中：v 為進(jìn)口速度，m/s；S 為進(jìn)口面積，m2；φ 為顆粒體積分?jǐn)?shù)。

實(shí)驗(yàn)中顆粒體積分?jǐn)?shù)為2%，當(dāng)進(jìn)口速度分別為1.5，2.5，4.0，5.5 m/s時(shí)，可計(jì)算得到流量分別為 0.021，0.034，0.055，0.075 kg/s。

3 結(jié)果與討論

3.1 顆粒分布規(guī)律

為了展示顆粒在裂縫中的沉降動(dòng)力學(xué)行為，圖4對(duì)比了相同進(jìn)口速度（4.0 m/s）、不同時(shí)間（2，5，10，30 s）縫內(nèi)的顆粒分布，圖5對(duì)比了50 s內(nèi)、不同進(jìn)口速度條件下的顆粒沉降特征。

在2 s時(shí)，顆粒流前端未表現(xiàn)出十分明顯的沉降特征。隨著顆粒向前推進(jìn)，在5 s時(shí)，顆粒到達(dá)裂縫出口端，由于顆粒沉降，出口端的上部空間形成了“純水區(qū)”，部分顆粒在液相攜帶下從下方的出口端流出。隨注入時(shí)間的增加，縫內(nèi)流動(dòng)逐漸達(dá)到動(dòng)態(tài)平衡，“純水區(qū)”的面積基本不再發(fā)生變化。需要說(shuō)明的是，由于底部邊界設(shè)置為“捕獲”，沉降至裂縫底部的顆粒不顯示在圖中。

當(dāng)v=1.5 m/s時(shí)，大量顆粒在裂縫入口附近沉積，顆粒能夠到達(dá)的最遠(yuǎn)距離為2.6 m。無(wú)顆粒從出口流出，出現(xiàn)大面積的“純水區(qū)”。當(dāng)v=2.5 m/s時(shí)，顆粒的沉降趨勢(shì)仍較為明顯，但有部分顆粒從出口流出。隨進(jìn)口速度的增加，從出口流出的顆粒數(shù)增多，“純水區(qū)”的面積逐漸減小。

圖4 不同時(shí)間內(nèi)縫內(nèi)顆粒的分布情況

圖5 不同進(jìn)口速度條件下縫內(nèi)顆粒的分布情況

對(duì)不同時(shí)間裂縫內(nèi)沉降與流出的顆粒數(shù)進(jìn)行統(tǒng)計(jì)，計(jì)算沉降顆粒占總顆粒數(shù)的百分比（見(jiàn)圖6）。由圖可知，顆粒流出裂縫前，顆粒沉降百分比增長(zhǎng)迅速。在前5 s內(nèi)，入口附近的湍流現(xiàn)象復(fù)雜，沉降顆粒數(shù)隨v的增大而增加；隨著縫內(nèi)流動(dòng)狀態(tài)逐漸穩(wěn)定，顆粒沉降百分比的增長(zhǎng)速度變緩；隨著v從1.5 m/s增大到5.5 m/s，漿體從入口到達(dá)出口的時(shí)間變短，縫內(nèi)流動(dòng)更快達(dá)到平衡狀態(tài)，顆粒沉降百分比達(dá)到穩(wěn)定的時(shí)間由30 s縮短至10 s，而顆粒沉降百分比由91.4%減小至50.2%。

圖6 不同進(jìn)口速度下的顆粒沉降百分比

3.2 顆粒沉降百分比的影響因素

3.2.1 顆粒直徑

從圖7可以看出：顆粒直徑越小，曲線的平緩段出現(xiàn)越早，說(shuō)明隨顆粒直徑的增大，顆粒沉降達(dá)到平衡需要的時(shí)間越長(zhǎng)；顆粒直徑從1.0 mm增大到2.0 mm，沉降的顆粒數(shù)增多，顆粒沉降百分比增大；在較短時(shí)間內(nèi)，3種直徑顆粒的沉降百分比差異不大，隨注入時(shí)間的增加，三者間的差異逐漸加大，直至趨于穩(wěn)定；注入時(shí)間為30 s時(shí)，直徑1.0 mm顆粒的沉降百分比為64.8%，而直徑2.0 mm顆粒的為36.9%。

圖7 顆粒直徑對(duì)顆粒沉降百分比的影響

3.2.2 攜帶液密度

從圖8可以看出，攜帶液密度越大，曲線平緩段出現(xiàn)越早。這是因?yàn)殡S攜帶液密度的增加，顆粒與攜帶液之間的重力分異減小，顆粒在液相中的沉降速度降低，則沉降達(dá)到動(dòng)態(tài)平衡所需的時(shí)間減少，同時(shí)顆粒沉降百分比也減小。

圖8 攜帶液密度對(duì)顆粒沉降百分比的影響

3.2.3 攜帶液黏度

注入時(shí)間為5 s時(shí)，隨著攜帶液黏度由1 mPa·s增加到3 mPa·s，攜帶液的攜帶性能增加，顆粒沉降百分比由29.9%減小至7.3%。隨注入時(shí)間的增加，顆粒沉降百分比逐漸增大，并最終趨于穩(wěn)定（見(jiàn)圖9）。

圖9 攜帶液黏度對(duì)顆粒沉降百分比的影響

4 結(jié)論

1）低密度調(diào)流劑顆粒在裂縫中的輸送實(shí)驗(yàn)顯示，顆粒會(huì)沉降在裂縫底部，且堆積高度最終會(huì)增大到一個(gè)平衡值。從顆粒沉降結(jié)果可以看出：攜帶液密度和黏度越大，最終堆積高度就越小，堆積位置越靠后；反之亦然。

2）隨注入時(shí)間的增加，顆粒沉降百分比迅速增長(zhǎng)并逐漸達(dá)到穩(wěn)定，其穩(wěn)定時(shí)間與沉降百分比隨進(jìn)口速度、攜帶液密度和黏度的增大而減小，隨顆粒直徑的增大而增大。

3）針對(duì)低密度、低濃度顆粒流建立的CFD-DPM模型能夠很好地揭示調(diào)流劑顆粒進(jìn)入裂縫后的不同時(shí)間內(nèi)顆粒的流動(dòng)及分布狀態(tài)，可定量描述顆粒沉降規(guī)律，為施工參數(shù)的設(shè)計(jì)與優(yōu)化提供指導(dǎo)。