摘 要：附壁液滴的變形和運(yùn)移與殘余油的微觀驅(qū)油效率密切相關(guān)，鑒于聚合物驅(qū)替液及稠油均可表現(xiàn)出黏彈性流變性質(zhì)，闡明兩相流體黏度等物性條件對(duì)黏彈性流體系統(tǒng)附壁液滴行為的影響機(jī)制十分必要。采用黏彈性流體兩相流格子玻爾茲曼方法，結(jié)合接觸角幾何公式處理潤(rùn)濕邊界條件，開展了線性剪切流中Oldroyd-B黏彈性環(huán)境流體內(nèi)牛頓液滴（N/V系統(tǒng)）及牛頓環(huán)境流體內(nèi)Oldroyd-B黏彈性液滴（V/N系統(tǒng)）在兩相流體黏度比m、黏彈性流體溶劑黏度比β影響下的穩(wěn)態(tài)變形及運(yùn)移規(guī)律研究。結(jié)果表明：N/V系統(tǒng)聚合物的拉伸在界面尖端外側(cè)及前進(jìn)接觸線附近最為顯著，形成彈性拉力促進(jìn)鋪展；隨β降低（彈性增強(qiáng)），潤(rùn)濕面積增加，表面積減小；隨m增加，潤(rùn)濕面積增加，表面積增大；N/V液滴移動(dòng)速度在mlt;1.5時(shí)，隨β降低而升高；但m≥1.5時(shí)，液滴所受總驅(qū)動(dòng)力在β較小時(shí)因液滴高度顯著降低而減小，導(dǎo)致其移動(dòng)速度隨β降低而減慢。V/N液滴內(nèi)聚合物拉伸在后退接觸線附近更強(qiáng)，抑制液滴鋪展；m≤1時(shí)潤(rùn)濕面積受m及β影響不大；mgt;1時(shí)潤(rùn)濕面積隨m增加、β升高而增大。V/N液滴表面積受m和β的影響與N/V液滴趨勢(shì)相一致，且兩系統(tǒng)內(nèi)液滴移動(dòng)速度均隨m增加而減小，但V/N液滴在各m條件下始終在β較小時(shí)移動(dòng)速度略高。

關(guān)鍵詞：聚合物驅(qū)油；黏彈性流體；液滴；格子玻爾茲曼方法；變形和運(yùn)移

中圖分類號(hào)：TE319" 文獻(xiàn)標(biāo)志碼：A文章編號(hào)：1000-4939（2025）01-0212-11

Numerical study on deformation and motion of a wall-attached droplet in linear shear flow with viscoelastic rheology

WANG Ningning，NI Wanglai，WANG Dong，LIU Haihu

（School of Energy and Power Engineering，Xi’an Jiaotong University，710049 Xi’an，China）

Abstract：The deformation and movement of wall-attached droplets are closely related to the microscale displacement efficiency of residual oil.Since both the displacing liquid in polymer flooding and the heavy oil can exhibit viscoelasticity，it is necessary to clarify the influence of the viscosity-related fluid properties on the behavior of wall-attached droplets in viscoelastic fluid systems.The two-phase viscoelastic lattice Boltzmann method，combined with the geometric wetting boundary condition，is used to investigate the steady-state deformation and motion of a droplet under shear flow，where either the matrix phase （N/V system） or the droplet phase （V/N system） is the Oldroyd-B viscoelastic fluid while the other phase is the Newtonian fluid.The effects of the two-phase viscosity ratio m and the solvent viscosity ratio β of the viscoelastic fluid are studied.Results show that the stretching of polymer molecules in the N/V system is the most pronounced outside the droplet tip and near the advancing contact lines，which acts to promote droplet spreading.The wetting area of the Newtonian droplet in the Oldroyd-B matrix is increased by either decreasing β （increasing elasticity） or increasing m.The droplet surface area in the N/V system also expands with increased m， but it reduces with the decrease of β.In the N/V system，the sliding velocity of the droplet is enhanced with decreased β for mlt;1.5，whereas for m≥1.5 the trend is reversed owing to the reduced overall driving force caused by the much lower droplet height at smaller β.In the V/N system，the polymer molecules inside the droplet are stretched more significantly near the receding contact lines，which inhibits the droplet’s spreading.When m≤1，the wetting area is hardly affected by m or β，but the wetting area increases with increased m or β for mgt;1.Under the effects of m and β，the variations of the droplet surface area in the V/N system are consistent in trend with those in the N/V system.In both N/V and V/N systems，the droplet slides more slowly at a higher m.Unlike the droplet in the N/V system，the droplet in the V/N system with a smaller β always moves slightly faster under a given m.

Key words：polymer flooding;viscoelastic fluid;droplet;lattice Boltzmann method;deformation and movement

油氣資源是國(guó)民經(jīng)濟(jì)領(lǐng)域重要的能源基礎(chǔ)和化工原料［1］。鑒于我國(guó)大部分陸相油田非均質(zhì)性十分嚴(yán)重，且已動(dòng)用油藏經(jīng)過多年的注水開發(fā)普遍處于高含水開采階段，提高采收率技術(shù)的發(fā)展面臨更艱巨的挑戰(zhàn)。聚合物驅(qū)作為一種提高采收率技術(shù)，在我國(guó)已形成了一定規(guī)模的工業(yè)化生產(chǎn)能力，成為油田增儲(chǔ)上產(chǎn)的重要措施［2］。然而，聚合物溶液驅(qū)殘余油的兩相流動(dòng)過程受到流體屬性、儲(chǔ)層條件、殘余油分布狀態(tài)、注采條件等多重因素的影響［3-5］，目前對(duì)聚合物驅(qū)油機(jī)理的認(rèn)識(shí)仍存在爭(zhēng)議，無(wú)法回答聚合物驅(qū)在部分礦場(chǎng)應(yīng)用中采收率偏低的問題［6-7］。微觀尺度上殘余油常常以油膜狀、盲端型、柱狀或者簇狀等分散態(tài)存在［8］，聚合物溶液驅(qū)油過程實(shí)質(zhì)上是多孔介質(zhì)內(nèi)復(fù)雜的多相界面動(dòng)力學(xué)流動(dòng)過程，闡明微觀尺度兩相流體的相互作用機(jī)制，是深入認(rèn)識(shí)聚合物驅(qū)油機(jī)理，進(jìn)而提高殘余油微觀驅(qū)油效率的重要手段之一。

在各種類型的殘余油中，附著于固壁的油膜、油滴狀殘余油是一種典型且廣泛的存在狀態(tài)［8］。外部驅(qū)動(dòng)流體在壁面附近形成的流場(chǎng)可近似看作線性剪切流，驅(qū)動(dòng)壁面油滴發(fā)生小幅振蕩、穩(wěn)定的移動(dòng)、部分分離或完全分離等［9-10］。目前，針對(duì)壁面附著液滴在線性剪切流條件下的動(dòng)力學(xué)行為研究已經(jīng)開展了較多工作，發(fā)現(xiàn)液滴行為與雷諾數(shù)、毛細(xì)數(shù)、兩相黏度比m、壁面潤(rùn)濕性等因素有關(guān)［11］。針對(duì)兩相流體屬性的影響，LIU等［12］研究發(fā)現(xiàn)小毛細(xì)數(shù)條件下液滴穩(wěn)態(tài)移動(dòng)速度隨黏度比的增加而單調(diào)降低。DIMITRAKOPOULOS 等［13］研究發(fā)現(xiàn)液滴在垂直于流動(dòng)方向的鋪展隨液滴黏度的增加而增大，并解釋為液滴界面張力與增大的水動(dòng)力學(xué)力相平衡的結(jié)果。

結(jié)合聚合物驅(qū)油過程，聚合物溶液實(shí)質(zhì)上具有非牛頓黏彈性性質(zhì)［14］；同時(shí)，稠油在不同的組分和溫度條件下，也可具有黏彈性流變性質(zhì)［15］。在流體黏彈性流變屬性的影響下，不僅液滴直接受到彈性力的作用，而且流場(chǎng)分布也受到彈性力的影響［16-17］。劉麗麗等［18］和朱光普等［19］數(shù)值研究均發(fā)現(xiàn)，聚合物溶液的黏彈性導(dǎo)致壁面液滴所受的法向偏應(yīng)力增大，促進(jìn)液滴變形；LIU等［20］數(shù)值研究發(fā)現(xiàn)高濃度或高分子質(zhì)量的高彈性聚合物溶液，促進(jìn)油滴變形、運(yùn)移和剝離。VARAGNOLO等［21］實(shí)驗(yàn)研究了傾斜壁面上黏彈性液滴的滑移，發(fā)現(xiàn)與穩(wěn)態(tài)移動(dòng)速度相關(guān)的毛細(xì)數(shù)與邦德數(shù)之間不再遵循牛頓流體系統(tǒng)中的線性關(guān)系，而是在邦德數(shù)較大時(shí)受彈性力的影響而趨近于常數(shù)。目前，不僅圍繞黏彈性環(huán)境流體內(nèi)牛頓液滴的研究工作不夠充分［19］，而且考慮油滴黏彈性時(shí)，即使對(duì)于更基礎(chǔ)水驅(qū)條件下黏彈性油滴的變形和運(yùn)動(dòng)規(guī)律也有待進(jìn)一步探索，尤其尚未闡明兩相流體黏度、黏彈性流體濃度（或溶劑黏度比）等重要物性參數(shù)的影響機(jī)制。

黏彈性流體彈性應(yīng)力的演化常通過本構(gòu)方程進(jìn)行描述，其中Oldroyd-B本構(gòu)模型作為一種典型的黏彈性流體本構(gòu)模型，在石油領(lǐng)域流動(dòng)研究中應(yīng)用廣泛［4］，其所表征的黏彈性流體在線性剪切流動(dòng)中黏度始終為常數(shù)、具有非零的第一法向應(yīng)力差?；陴椥粤黧w系統(tǒng)壁面附著液滴動(dòng)力學(xué)行為的研究現(xiàn)狀，本研究采用Oldroyd-B本構(gòu)模型描述流體的黏彈性流變性質(zhì)，基于格子玻爾茲曼方法（lattice Boltzmann method，LBM）分別求解兩相流動(dòng)控制方程和黏彈性流體本構(gòu)方程，結(jié)合DING等［22］所發(fā)展的潤(rùn)濕邊界條件，分別開展液滴黏彈性和環(huán)境流體黏彈性兩相流系統(tǒng)中兩相流體黏度比、黏彈性流體溶劑黏度比對(duì)壁面附著液滴變形及運(yùn)動(dòng)影響機(jī)制的研究，以期加深對(duì)流體物性條件影響下驅(qū)油機(jī)理的認(rèn)識(shí)，并可為化學(xué)、生物等領(lǐng)域非牛頓流體系統(tǒng)液滴微流動(dòng)控制中的共性問題研究提供理論參考。

1 壁面液滴運(yùn)動(dòng)模型

1.1 物理模型

如圖1所示，半徑R=25的半球形液滴置于靜止的下壁面上，計(jì)算域長(zhǎng)、寬、高設(shè)置為L(zhǎng)×W×H=11.2R×8R×2R，上平板以速度uw沿著x正方向移動(dòng)，計(jì)算域四周為周期邊界，流場(chǎng)特征剪切率為γ·=uw/H。采用ρR、ρB和μR、μB分別表示兩相流體的密度和零剪切黏度，其中下標(biāo)“R”、“B”分別指代紅色液滴相流體和藍(lán)色環(huán)境負(fù)載流體，σ表示兩相流體間的界面張力系數(shù)，則剪切流場(chǎng)中液滴運(yùn)動(dòng)的主要無(wú)量綱控制參數(shù)［9，23-24］可表示為：毛細(xì)數(shù)Ca=μBγ·R/σ，定義為黏性力與界面張力之比；雷諾數(shù)Re=ρBγ·R2/μB，定義為慣性力與黏性力之比；魏森伯格數(shù)Wi=λγ·，定義為黏彈性流體的松弛時(shí)間λ與流動(dòng)特征時(shí)間1/γ·之比；液滴與環(huán)境負(fù)載流體間的兩相黏度比m=μR/μB。此外，以液滴相為黏彈性流體時(shí)為例，液滴零剪切黏度μR可視為由該流體的溶質(zhì)聚合物黏度μRp和溶劑牛頓黏度μRs兩部分組成，定義溶劑黏度比β=μRs/μR反映黏彈性流體的彈性濃度［25］。β在0～1.0之間變化；當(dāng)β=0時(shí)，表示純彈性流體；當(dāng)β=1.0時(shí)，表示牛頓流體。為簡(jiǎn)單起見，本研究忽略兩相流體的密度差異。

1.2 數(shù)值模型及邊界條件

液滴及其環(huán)境流體組成的黏彈性兩相流體系統(tǒng)，其流場(chǎng)演化遵循Navier-Stokes方程，即

給定計(jì)算域x×y×z為199×10×35，界面張力σ=0.001，潤(rùn)濕角θ=60°，藍(lán)色流體零剪切黏度μB固定為0.1；黏彈性流體β=0.5，λ=10000。進(jìn)、出口采用修正周期邊界條件［32］以保證自吸過程中紅色流體從進(jìn)口連續(xù)吸入，藍(lán)色流體從出口連續(xù)流出，且進(jìn)、出口位置處壓力相等。在兩相黏度比m分別為0.125、1.0、4.0的條件下，分別模擬了純牛頓流體系統(tǒng)（N/N）、被驅(qū)替相為Oldroyd-B流體（N/V），以及驅(qū)替相為Oldroyd-B流體（V/N）時(shí)的自發(fā)吸入過程。模擬所得的界面位置時(shí)間演化曲線如圖3所示，N/N、N/V及V/N流體系統(tǒng)在所測(cè)試的3種黏度比條件下，界面抵達(dá)出口前（l≤190）與相應(yīng)解析解間的最大相對(duì)誤差分別N/N≤1.7%，N/V≤4.0%，V/N≤2.0%，驗(yàn)證了所發(fā)展模型能夠較好地模擬移動(dòng)接觸線行為。

3 計(jì)算結(jié)果與討論

針對(duì)線性剪切流條件下的壁面附著液滴運(yùn)動(dòng)，采用Re=1.0表征慣性可忽略的低速運(yùn)動(dòng)情況；在彈性數(shù)的選擇方面，已有文獻(xiàn)較常選用Wi=0～1.0［19-20］，本研究采用Wi=1.0體現(xiàn)彈性的作用；壁面接觸角設(shè)置為90°，對(duì)應(yīng)為中性潤(rùn)濕條件。同時(shí)，由牛頓流體系統(tǒng)壁面附著液滴運(yùn)動(dòng)的研究可知［9，12］，Ca較小時(shí)液滴僅發(fā)生變形運(yùn)動(dòng)，而Ca超過某一臨界值后，液滴可發(fā)生破裂、甚至從壁面脫離。鑒于目前對(duì)黏彈性流體系統(tǒng)中黏度比、溶劑黏度比等參數(shù)影響下的液滴受力特點(diǎn)認(rèn)識(shí)不十分清晰，且變形運(yùn)動(dòng)機(jī)制是理解脫離、破裂等行為的重要前提，本研究中固定Ca=0.15以模擬液滴的穩(wěn)態(tài)變形和運(yùn)移。

3.1 黏彈性流體剪切流中的牛頓液滴（N/V系統(tǒng)）

在β=0.25、0.5、0.75的條件下，通過改變液滴相流體的黏度值，模擬了兩相黏度比m為0.125、0.5、1、1.5、2和4時(shí)，牛頓液滴在黏彈性流體剪切流（N/V系統(tǒng)）中的變形與運(yùn)移行為。在所模擬的各參數(shù)條件下，液滴均在外部剪切流的作用下發(fā)生變形和運(yùn)動(dòng)，液滴形態(tài)從初始時(shí)刻起發(fā)生持續(xù)的變形直至穩(wěn)定，并最終以恒定的速度ucl在壁面上移動(dòng)。為描述液滴的穩(wěn)態(tài)變形程度，定義液滴潤(rùn)濕面積、高度、及表面積的相對(duì)變化參數(shù)分別為Ar=（A－A0）/A0，hr=（h－h(huán)0）/h0，Sr=（S－S0）/S0，其中A、h、S和A0、h0、S0分別為液滴穩(wěn)態(tài)和初始潤(rùn)濕面積、高度及表面積。圖4給出了穩(wěn)態(tài)時(shí)Ar、hr、Sr在不同β條件下隨兩相黏度比m的變化曲線，以及β=0.25時(shí)x-y底面和x-z中截面上的界面形態(tài)。研究發(fā)現(xiàn)：隨著兩相黏度比m的增加，液滴潤(rùn)濕面積Ar和表面積Sr均隨m的增加而增大，而液滴高度hr隨m的增加而減小；各β條件下液滴隨m的變形趨勢(shì)相一致，與純牛頓流體系統(tǒng)中黏度比的影響相一致［12］。

從液滴受力的角度，黏彈性流體系統(tǒng)中液滴所受驅(qū)動(dòng)力主要包括黏性力、彈性力和壓力3部分，液滴穩(wěn)態(tài)時(shí)x方向的驅(qū)動(dòng)力與壁面摩擦力相平衡。界面上各驅(qū)動(dòng)力x方向分量在兩相流體界面SI上的合力的計(jì)算式分別為：黏性力∫SIn·τsxδΓdS、彈性力∫SIn·τpxδΓdS和壓力∫SIn·（－pI）xδΓdS。其中，δΓ為狄拉克函數(shù)，對(duì)于格子玻爾茲曼顏色模型［12］，狄拉克函數(shù)可取為δΓ=SymbolQC@ρN/2?；谝陨狭Φ挠?jì)算公式，采用穩(wěn)態(tài)時(shí)的模擬數(shù)據(jù)驗(yàn)證發(fā)現(xiàn)，各驅(qū)動(dòng)力之和與壁面摩擦阻力相平衡，表明以上各作用力計(jì)算式的合理性。

通過μBγ·πR2對(duì)界面所受黏性力、彈性力和壓力進(jìn)行無(wú)量綱化處理，并分別用符號(hào)n·τs′x、n·τp′x、－n·pI′x表示，各分力隨兩相黏度比的變化如圖5所示。聚合物構(gòu)型張量的跡trA，定義為構(gòu)型張量A主對(duì)角線上的元素之和，其大小與聚合物分子的拉伸程度正相關(guān)［25］，故從trA的分布能夠直觀反映聚合物的拉伸情況（對(duì)于牛頓流體，trA=3.0）。后文將從液滴界面所受x方向驅(qū)動(dòng)力隨控制參數(shù)的變化，結(jié)合聚合物拉伸情況討論液滴的變形機(jī)制。

以β=0.25為例，圖6（a）～圖6（c）給出了m=0.125、1.0、4.0時(shí)x-y底面和x-z中截面上的trA分布。黏彈性聚合物在界面尖端外側(cè)（圖6子圖中的虛線部分）及前進(jìn)接觸線附近拉伸最為顯著，形成的彈性拉力促進(jìn)液滴在壁面的鋪展并降低了液滴的高度，且其拉伸程度隨m的增加而增強(qiáng)，導(dǎo)致液滴界面所受的彈性力隨之增大，與圖5中的統(tǒng)計(jì)結(jié)果相一致。另外，液滴黏度的增加減慢了液滴的運(yùn)動(dòng)速度，導(dǎo)致液滴經(jīng)受更強(qiáng)的剪切作用，所以圖5中黏性驅(qū)動(dòng)力也隨兩相黏度比的增加而增大。

由圖5還可以看出，液滴所受的壓力隨兩相黏度比的增加由阻力、或影響不大轉(zhuǎn)變?yōu)轵?qū)動(dòng)力的作用?？傮w上，液滴界面所受黏性力、彈性力、壓力3者合力隨兩相黏度比增加而增大，促進(jìn)了液滴表面積的增加。

在各兩相黏度比m條件下，比較圖4環(huán)境流體溶劑黏度比β

對(duì)液滴穩(wěn)態(tài)變形的影響，發(fā)現(xiàn)液滴均在β最小時(shí)（β=0.25）更易在壁面上鋪展，對(duì)應(yīng)更大的潤(rùn)濕面積Ar和更低的液滴高度hr。圖6（c）～圖（e）給出了m=4.0時(shí)trA隨β的變化，雖然黏彈性聚合物的拉伸隨溶劑黏度比β的增加而增強(qiáng)，但由于溶質(zhì)黏度系數(shù)μp隨β的增加而降低，導(dǎo)致由式（3）計(jì)算所得的彈性力在β=0.25時(shí)更大，因此隨著β的減小，潤(rùn)濕面積Ar增大而液滴高度hr減小。由圖5（b）可知，在黏彈性環(huán)境流體的溶劑黏度μs和溶質(zhì)黏度μp相等時(shí)（β=0.5），x方向黏性力分量明顯高于彈性力分量。相同黏度比條件下，隨β減小（彈性增強(qiáng)），黏性力減小、液滴所受總的合力減小，但彈性力的增加促進(jìn)了液滴潤(rùn)濕面積的增大，表明液滴潤(rùn)濕面積的改變對(duì)于彈性力的變化更為敏感。另一方面，各m條件下液滴表面積均在β=0.75時(shí)最大，與液滴所受黏性力、以及總的合力略增大的趨勢(shì)相一致，由此可見液滴表面積的變化由黏性力和彈性力的變化共同決定。

圖7給出了N/V系統(tǒng)液滴穩(wěn)態(tài)移動(dòng)速度ucl/uw隨兩相黏度比m的變化，由于壁面摩擦阻力作用隨液滴黏度的增加而增強(qiáng)，各溶劑黏度比β條件下液滴穩(wěn)態(tài)移動(dòng)速度均隨m的增加而降低。由圖7可知：液滴移動(dòng)速度在兩相黏度比m＜1.5時(shí)，β越小移動(dòng)速度越快，其原因在于m較低時(shí)，液滴潤(rùn)濕面積及表面積均變形較小，環(huán)境流體彈性拉力在β較小時(shí)占比相對(duì)較大，有助于促進(jìn)液滴運(yùn)動(dòng)；而m≥1.5時(shí)，β越小潤(rùn)濕面積越大，液滴高度明顯降低，以圖5中m=3為例，β=0.75時(shí)（圖5c）的黏性力分量明顯大于β=0.25時(shí)（圖5a）的黏性力分量，且液滴所受合力隨β減小而略減弱，所以β越大液滴移動(dòng)速度越快。

3.2 牛頓流體剪切流中的黏彈性液滴（V/N系統(tǒng)）

圖8給出了V/N系統(tǒng)液滴穩(wěn)態(tài)時(shí)Ar、hr、Sr在不同β條件下隨兩相黏度比m的變化曲線，以及β=0.25時(shí)x-y底面及x-z中截面上的穩(wěn)定界面形態(tài)。

當(dāng)m≤1時(shí)，黏彈性液滴Ar及hr受m及β的影響不大；當(dāng)mgt;1時(shí)，Ar隨m的增大而增加，hr隨m的增大而減小，Ar及hr變化幅度均小于黏彈性環(huán)境流體中牛頓液滴相應(yīng)的變化幅度。表面積Sr總體隨m增大而增大；相同m條件下，Sr隨β增大而增大。

圖9～圖10分別給出了各β條件下，黏彈性液滴界面所受黏性力、彈性力、壓力x方向分量以及液滴穩(wěn)態(tài)移動(dòng)速度隨兩相黏度比m的變化。

隨著黏彈性液滴黏度的增加，液滴所受的壁面摩擦力增大，導(dǎo)致液滴移動(dòng)速度降低，液滴與環(huán)境流體間的速度梯度增大，因此總體上液滴黏性力和彈性力均隨黏度比的增加而增大；液滴所受的壓力總體占比較小。與N/V系統(tǒng)類似，液滴所受驅(qū)動(dòng)力合力的增大導(dǎo)致液滴表面積隨黏度比增加而增大。

另一方面，圖11給出了β=0.25時(shí)，V/N系統(tǒng)液滴穩(wěn)態(tài)時(shí)x-y底面及x-z中截面上的trA分布。

圖11（a）內(nèi)虛線框標(biāo)識(shí)了強(qiáng)彈性拉伸區(qū)域在m=0.125時(shí)，黏彈性聚合物的拉伸主要集中在壁面附近的流體上，以液滴后退接觸線區(qū)域及前進(jìn)接觸線包圍的內(nèi)部區(qū)域最為顯著（圖11a～b）。因此，盡管m≤1時(shí)液滴界面所受合力也隨m增加而增大，但液滴底部較強(qiáng)的彈性拉伸作用表現(xiàn)為抵抗液滴變形，導(dǎo)致潤(rùn)濕面積和液滴高度受m及β的影響并不明顯。mgt;1時(shí)，隨著m的增加，以圖11（c） 中m=4.0時(shí)為例，聚合物的拉伸明顯減弱且僅體現(xiàn)在液滴后退接觸線附近區(qū)域，界面總驅(qū)動(dòng)力增大，因此液滴潤(rùn)濕面積增大、高度降低；同時(shí)，彈性抑制鋪展的作用在β

較低時(shí)相對(duì)更強(qiáng)，所以mgt;1范圍內(nèi)，Ar隨β降低而減小，hr隨β降低而增大，與N/V系統(tǒng)中β的影響相反。VARAGNOLO等［21］的二維數(shù)值研究也發(fā)現(xiàn)了液滴內(nèi)聚合物拉伸類似的分布特點(diǎn)。

在對(duì)黏彈性液滴移動(dòng)速度的影響方面，由圖10知，V/N液滴在各m條件下的穩(wěn)態(tài)移動(dòng)速度ucl/uw始終在β較小時(shí)更高；在各β條件下均隨兩相黏度比m的增加而降低。其原因在于：β較小或兩相黏度比m較低時(shí)，液滴黏彈性抑制變形的能力較強(qiáng)，所以液滴底面積擴(kuò)張較小、壁面摩擦力小，液滴高度更高，承受的剪切作用較強(qiáng)，因此液滴移動(dòng)速度更高。

4 結(jié) 論

本研究在黏彈性兩相流顏色模型的基礎(chǔ)上，通過結(jié)合潤(rùn)濕邊界條件，發(fā)展了一種能夠同時(shí)處理黏彈性流變和移動(dòng)接觸線問題的兩相流格子玻爾茲曼方法，推導(dǎo)了Oldroyd-B 流體毛細(xì)自吸現(xiàn)象界面位置時(shí)間演化的解析解，開展了N/V和V/N兩類系統(tǒng)線性剪切流中壁面附著液滴在兩相流體黏度比m、黏彈性流體溶劑黏度比β影響下的變形和運(yùn)移規(guī)律研究，結(jié)論如下。

1）N/V系統(tǒng)聚合物拉伸以液滴尖端外側(cè)和前進(jìn)接觸線區(qū)域最為顯著，彈性拉力促進(jìn)液滴在壁面的鋪展；而V/N系統(tǒng)液滴的彈性拉伸主要集中在后退接觸線附近區(qū)域，表現(xiàn)為抑制液滴鋪展。

2）N/V和V/N液滴所受彈性力和黏性力均隨m增加而增大，N/V液滴隨m增大表現(xiàn)為潤(rùn)濕面積增加、高度降低及表面積增大。V/N液滴在m≤1時(shí)潤(rùn)濕面積及高度受m影響不大，在m更高時(shí)與N/V液滴變形趨勢(shì)相一致，但幅度更小。

3）黏彈性流體溶劑黏度比β增加時(shí)，N/V液滴由于受到的彈性促進(jìn)作用減小，其潤(rùn)濕面積減小。V/N液滴潤(rùn)濕面積受β的影響在m≤1時(shí)不明顯，在mgt;1時(shí)表現(xiàn)為隨β增加（彈性抑制作用減弱）而增大。

4）N/V及V/N液滴穩(wěn)態(tài)移動(dòng)速度均隨兩相黏度比m的增加而降低。同時(shí)，V/N液滴始終在β較低時(shí)由于潤(rùn)濕面積偏小而速度略高。N/V液滴在m較小時(shí)，由于彈性力在β較低時(shí)的驅(qū)動(dòng)作用更強(qiáng)而運(yùn)動(dòng)速度更快；但m≥1.5時(shí)，N/V液滴因在β較低時(shí)的鋪展增強(qiáng)而移動(dòng)更慢。

參考文獻(xiàn)：

［1］ 馮丹，陸朝榮，董仕寶.我國(guó)國(guó)內(nèi)石油供給安全形勢(shì)分析［J］.廣州化工，2020，48（8）：21-23.

FENG Dan，LU Zhaorong，DONG Shibao.Analysis of Chinese domestic oil supply security situation［J］.Guangzhou chemical industry，2020，48（8）：21-23（in Chinese）.

［2］ 金亞杰.國(guó)外聚合物驅(qū)油技術(shù)研究及應(yīng)用現(xiàn)狀［J］.非常規(guī)油氣，2017，4（1）：116-122.

JIN Yajie.Progress in research and application of polymer flooding technology abroad［J］.Unconventional oil amp; gas，2017，4（1）：116-122（in Chinese）.

［3］ 王德民.黏彈性流體的特殊性對(duì)油藏工程、地面工程及采油工程的影響［J］.大慶石油學(xué)院學(xué)報(bào)，2001，25（3）：46-52.

WANG Demin.Effect of particularity of visco-elastic fluid on reservoir engineering and production engineering［J］.Journal of Daqing Petroleum Institute，2001，25（3）：46-52（in Chinese）.

［4］ 鐘會(huì)影，張偉東，劉義坤，等.黏彈性聚合物驅(qū)普通稠油微觀滲流數(shù)學(xué)模型［J］.地球科學(xué)，2017，42（8）：1364-1372.

ZHONG Huiying，ZHANG Weidong，LIU Yikun，et al.The micro-flow mathematical model study on viscoelastic polymer displacement viscous oil［J］.Earth science，2017，42（8）：1364-1372（in Chinese）.

［5］ 劉西洋，付濤濤，朱春英，等.微通道內(nèi)非牛頓流體中液滴生成機(jī)理研究進(jìn)展［J］.化工學(xué)報(bào)，2021，72（2）：772-782.

LIU Xiyang，F(xiàn)U Taotao，ZHU Chunying，et al.Progress on droplet formation mechanism in non-Newtonian fluids in microchannels［J］.CIESC journal，2021，72（2）：772-782（in Chinese）.

［6］ VERMOLEN E C，HAASTERECHT M J，MASALMEH S K.A systematic study of the polymer visco-elastic effect on residual oil saturation by core flooding［C］//SPE EOR Conference at Oil and Gas West Asia.Richardson，Texas，USA：OnePetro，2014：SPE-169681-MS.

［7］ SANDENGEN K，MELHUUS K，KRISTOFFERSEN A.Polymer “viscoelastic effect”; does it reduce residual oil saturation［J］.Journal of petroleum science and engineering，2017，153：355-363.

［8］ ZHAO H，YANG H B，KANG X，et al.Study on the types and formation mechanisms of residual oil after two surfactant imbibition［J］.Journal of petroleum science and engineering，2020，195：107904.

［9］ DING H，GILANI M N H，SPELT P D M.Sliding，pinch-off and detachment of a droplet on a wall in shear flow［J］.Journal of fluid mechanics，2010，644：217-244.

［10］SEEVARATNAM G K，DING H，MICHEL O，et al.Laminar flow deformation of a droplet adhering to a wall in a channel［J］.Chemical engineering science，2010，65（16）：4523-4534.

［11］WANG F L，LIU T，LEI W H，et al.Dynamic analysis of deformation and start-up process of residual-oil droplet on wall under shear flow［J］.Journal of petroleum science and engineering，2021，199：108335.

［12］LIU H H，ZHANG J G，BA Y，et al.Modelling a surfactant-covered droplet on a solid surface in three-dimensional shear flow［J］.Journal of fluid mechanics，2020，897：A33.

［13］DIMITRAKOPOULOS P，HIGDON J J L.On the displacement of three-dimensional fluid droplets from solid surfaces in low-Reynolds-number shear flows［J］.Journal of fluid mechanics，1998，377：189-222.

［14］ROCK A，HINCAPIE R E，TAHIR M，et al.On the role of polymer viscoelasticity in enhanced oil recovery：extensive laboratory data and review［J］.Polymers，2020，12（10）：2276.

［15］徐丕東，王麗霞，蘇玉亮.黏彈性稠油的驅(qū)替特征及其影響因素分析［J］.西南石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2008，30（6）：113-116.

XU Pidong，WANG Lixia，SU Yuliang.Displacement characteristics of viscoelastic heavy oil and the analysis on influencing factors［J］.Journal of Southwest Petroleum University（science amp; technology edition），2008，30（6）：113-116（in Chinese）.

［16］MUKHERJEE S，SARKAR K.Effects of viscosity ratio on deformation of a viscoelastic drop in a Newtonian matrix under steady shear［J］.Journal of non-newtonian fluid mechanics，2009，160（2/3）：104-112.

［17］WANG D，WANG N N，LIU H H.Droplet deformation and breakup in shear-thinning viscoelastic fluid under simple shear flow［J］.Journal of rheology，2022，66（3）：585-603.

［18］劉麗麗，楊樹人，王立輝，等.微孔道中殘余油膜受力與變形的數(shù)值計(jì)算［J］.油氣地質(zhì)與采收率，2012，19（2）：60-63.

LIU Lili，YANG Shuren，WANG Lihui，et al.The numerical calculation about the stress and deformation of residual oil in the micropore［J］.Petroleum geology and recovery efficiency，2012，19（2）：60-63（in Chinese）.

［19］朱光普，姚軍，孫海，等.基于界面追蹤的黏彈性流體微觀流動(dòng)模擬［J］.科學(xué)通報(bào)，2016，61（36）：3973-3981.

ZHU Guangpu，YAO Jun，SUN Hai，et al.Microscopic flow simulation of viscoelastic fluid based on interface tracking method［J］.Chinese science bulletin，2016，61（36）：3973-3981（in Chinese）.

［20］LIU Y，F(xiàn)AN J W，LIU L L，et al.Numerical simulation of residual oil flooded by polymer solution in microchannels［J］.Geofluids，2018，2018（1）：8947839.

［21］VARAGNOLO S，F(xiàn)ILIPPI D，MISTURA G，et al.Stretching of viscoelastic drops in steady sliding［J］.Soft matter，2017，13（17）：3116-3124.

［22］DING H，SPELT P D M.Wetting condition in diffuse interface simulations of contact line motion［J］.Physical review e，2007，75（4）：046708.

［23］AGGARWAL N，SARKAR K.Effects of matrix viscoelasticity on viscous and viscoelastic drop deformation in a shear flow［J］.Journal of fluid mechanics，2008，601：63-84.

［24］AGGARWAL N，SARKAR K.Deformation and breakup of a viscoelastic drop in a Newtonian matrix under steady shear［J］.Journal of fluid mechanics，2007，584：1-21.

［25］WANG D，TAN D S，KHOO B C，et al.A lattice Boltzmann modeling of viscoelastic drops’ deformation and breakup in simple shear flows［J］.Physics of fluids，2020，32（12）：123101.

［26］MALASPINAS O，F(xiàn)ITIER N，DEVILLE M.Lattice Boltzmann method for the simulation of viscoelastic fluid flows［J］.Journal of non-newtonian fluid mechanics，2010，165（23/24）：1637-1653.

［27］OLDROYD J G.On the formulation of rheological equations of state［J］.Proceedings of the Royal Society of London.series a.mathematical and physical sciences，1950，200（1063）：523-541.

［28］HALLIDAY I，LAW R，CARE C M，et al.Improved simulation of drop dynamics in a shear flow at low Reynolds and capillary number［J］.Physical review e，2006，73（5）：056708.

［29］LIU H H，JU Y P，WANG N N，et al.Lattice Boltzmann modeling of contact angle and its hysteresis in two-phase flow with large viscosity difference［J］.Physical review e，2015，92（3）：033306.

［30］LADD A J C.Numerical simulations of particulate suspensions via a discretized Boltzmann equation.Part 1.Theoretical foundation［J］.Journal of fluid mechanics，1994，271：285-309.

［31］WASHBURN E W.The dynamics of capillary flow［J］.Physical review，1921，17（3）：273-283.

［32］LIU H H，SUN S L，WU R，et al.Pore-scale modeling of spontaneous imbibition in porous media using the lattice Boltzmann method［J］.Water resources research，2021，57（6）：e2020WR029219.

（編輯 張璐）