毛細(xì)管模型中流體啟動(dòng)流動(dòng)力學(xué)分析及應(yīng)用

玄建

（勝利油田勘探開發(fā)研究院，山東東營(yíng)257000）

毛細(xì)管模型中流體啟動(dòng)流動(dòng)力學(xué)分析及應(yīng)用

玄建

（勝利油田勘探開發(fā)研究院，山東東營(yíng)257000）

毛細(xì)管模型中，根據(jù)Poiseuille導(dǎo)管流動(dòng)公式和牛頓第二定律建立單相流體的啟動(dòng)流動(dòng)方程，求解得到位移、速度等參數(shù)隨時(shí)間的變化，進(jìn)而得出了毛細(xì)管中單相流體的啟動(dòng)流動(dòng)過程中慣性力和滲流阻力之間的變化規(guī)律。根據(jù)毛細(xì)管模型中流動(dòng)，闡述了啟動(dòng)壓力梯度問題。

毛細(xì)管；啟動(dòng)壓力梯度；流體

對(duì)砂巖儲(chǔ)集巖，常常把微小的孔隙看成毛細(xì)管，由于毛細(xì)管模型簡(jiǎn)單，便于分析，因此研究毛細(xì)管模型有利于對(duì)滲流問題的深入認(rèn)識(shí)。自Hagen-Poiseuille發(fā)表了導(dǎo)管流動(dòng)公式以來，毛細(xì)管模型在環(huán)境科學(xué)、水利工程、石油工業(yè)、醫(yī)學(xué)等方面得到了廣泛的應(yīng)用。很多學(xué)者做了大量的毛細(xì)管模型研究，大部分基于穩(wěn)定流態(tài)，對(duì)啟動(dòng)流動(dòng)的研究較少。在很多實(shí)驗(yàn)中，得到了啟動(dòng)壓力梯度，這正是測(cè)量初始階段啟動(dòng)流動(dòng)的反應(yīng)。很多學(xué)者對(duì)啟動(dòng)壓力梯度持懷疑態(tài)度，也給出了合理的解釋。本文通過毛細(xì)管模型中啟動(dòng)流動(dòng)的力學(xué)研究，進(jìn)一步闡釋了啟動(dòng)壓力梯度問題。

1 毛細(xì)管模型與達(dá)西公式

Hagen-Poiseuille[1]（1842年）發(fā)表了導(dǎo)管流動(dòng)公式，模型假設(shè)流體為不可壓縮牛頓流體并在管中作層流的條件下，得到了經(jīng)典的Poiseuille方程式：

Henry Darcy[2]（1856年）研究水在砂濾器中流動(dòng)時(shí)得出了經(jīng)驗(yàn)達(dá)西公式，用以描述流體流經(jīng)特定形狀的多孔介質(zhì)系統(tǒng)的性質(zhì)。

在毛細(xì)管中的流體流動(dòng)既遵循Poiseuille公式，也遵循達(dá)西公式，Kozeny公式給出了兩個(gè)公式之間的關(guān)系式。

2 毛細(xì)管中啟動(dòng)流動(dòng)方程的建立

假設(shè)在較長(zhǎng)毛細(xì)管中有長(zhǎng)為L(zhǎng)的流體，在其兩端施加Δp壓力差，由靜止開始運(yùn)動(dòng)。在運(yùn)動(dòng)過程中，Δp既要克服滲流阻力使流體運(yùn)動(dòng)，又要使流體做加速運(yùn)動(dòng)?？朔B流阻力使流體運(yùn)動(dòng)符合Poiseuille公式，使流體做加速運(yùn)動(dòng)符合牛頓第二運(yùn)動(dòng)定律。

根據(jù)位移、速度、加速度的關(guān)系可以得到：

初始條件為：當(dāng)t=0時(shí)，s=0，s'=0，求解得：

進(jìn)而求得速度v和加速度a：

根據(jù)牛頓第二運(yùn)動(dòng)定律作用在流體上并使其加速運(yùn)動(dòng)的壓力為：

根據(jù)以上計(jì)算可以看出當(dāng)毛細(xì)管半徑在微米級(jí)時(shí)，流體在極短時(shí)間內(nèi)加速度消失，損耗在加速上的壓力Δp'也隨之消失，流體達(dá)到了穩(wěn)定流動(dòng)狀態(tài)。

3 啟動(dòng)壓力梯度探討

有的學(xué)者通過實(shí)驗(yàn)方法提出啟動(dòng)壓力梯度是存在的（見圖1）[3-10]，有的學(xué)者認(rèn)為不存在，認(rèn)為啟動(dòng)壓力梯度只是一個(gè)實(shí)驗(yàn)假象[11]。從毛細(xì)管理論模型中流體啟動(dòng)流動(dòng)研究來看，只要壓力大于毛管力，毛細(xì)管中的流體就會(huì)流動(dòng)，因此啟動(dòng)壓力梯度是不存在的。而在實(shí)驗(yàn)過程中經(jīng)常得到一個(gè)啟動(dòng)壓力梯度，即在壓力梯度較小時(shí)，滲流規(guī)律并不符合Poiseuille公式和達(dá)西公式，也與毛細(xì)管模型中啟動(dòng)流動(dòng)公式不一致。

圖1 實(shí)驗(yàn)室得到的滲流指示曲線示意圖

實(shí)驗(yàn)中用的巖心可以看成有很多個(gè)不等徑毛細(xì)管組成的毛細(xì)管模型，由于在油藏條件下，油相基本上是非潤(rùn)濕相[12，13]。假設(shè)該模型中流體為非潤(rùn)濕相，那么毛管力就是流動(dòng)阻力。由于毛管力的作用，在較低的壓力梯度驅(qū)替時(shí)，大部分毛細(xì)管中的流體不流動(dòng)，只有較大的毛細(xì)管中的流體流動(dòng)。而這部分較大毛細(xì)管的橫截面積僅占整個(gè)巖心橫截面積的一部分，隨著驅(qū)替壓力的逐漸增加，參與流動(dòng)的毛細(xì)管越來越多，直至驅(qū)替壓力能夠驅(qū)動(dòng)絕大部分毛細(xì)管中的流體。

假設(shè)有一巖心由r1、r2（r1＞r2）兩種毛細(xì)管組成，毛細(xì)管中的流體為非潤(rùn)濕相，當(dāng)壓力達(dá)到pc1＜pr1＜pc2時(shí)，r1毛細(xì)管中流體流動(dòng)，流動(dòng)符合：

當(dāng)壓力達(dá)到pr2＞pc2時(shí)，r1、r2毛細(xì)管中流體流動(dòng)，流動(dòng)符合：

合并得：

因?yàn)閞1＞r2，所以，也就是說，pc1＜p1＜pc2時(shí)的斜率大于pr2＞pc2時(shí)的斜率，而在實(shí)驗(yàn)測(cè)量時(shí)，得到的斜率是逐漸增大的，這是因?yàn)閷?shí)驗(yàn)時(shí)采用的速度不是真實(shí)速度，而是視速度，即流量除以巖心的橫截面積得到，那么斜率就變成了，這樣，與實(shí)驗(yàn)結(jié)果一致。

同理可得，當(dāng)存在三種半徑r1、r2、r3（r1＞r2＞r3）或者更多種毛細(xì)管時(shí)，視速度V與關(guān)系圖中，斜率逐漸增大，與實(shí)驗(yàn)得出數(shù)據(jù)一致。

4 結(jié)論

（1）從流體在毛細(xì)管模型中由靜止開始運(yùn)動(dòng)過程的理論推導(dǎo)來看，毛細(xì)管中的流體在極短的時(shí)間內(nèi)就能達(dá)到穩(wěn)定流態(tài)，慣性力在極短的時(shí)間內(nèi)消失，壓力只用于克服滲流阻力。

（2）根據(jù)理想毛細(xì)管模型來看，實(shí)驗(yàn)測(cè)量的滲流指示曲線中出現(xiàn)的低速滲流曲線段和高速滲流直線段是合理，這是由于巖心中存在不同半徑的孔隙具有不同的毛管力導(dǎo)致的。即在較低的壓力梯度時(shí)，僅克服較大孔道的毛管力使其中的流體流動(dòng)，隨著壓力梯度的不斷增加，越來越多的更小半徑的孔隙中流體參與流動(dòng)。當(dāng)絕大部分孔道中的流體流動(dòng)時(shí)，滲流指示曲線呈現(xiàn)直線段。

符號(hào)注釋：

Q、Q1、、Q2-流量，cm3/s；r、r1、r2、r3-毛細(xì)管半徑，μm；Δp-毛細(xì)管內(nèi)流體兩端壓力差，MPa；μ-流體黏度，mPa·s；L-毛細(xì)管內(nèi)流體長(zhǎng)度，m；K-巖心滲透率，μm2；A-巖心橫截面積，cm2；v、v1、v2-毛細(xì)管內(nèi)流體的速度，m/s；V-流體視速度，m/s；p1、p2、pr1、pr2-壓力，MPa；ρ-毛細(xì)管內(nèi)流體的密度，kg/m3；a-毛細(xì)管內(nèi)流體的加速度，m/s2；s-毛細(xì)管內(nèi)流體的位移，m；t-時(shí)間，s；b、c、d-常數(shù)；Δp'-損耗在流體加速上的壓力，MPa；pc1、pc2-毛管力，MPa；m1、m2-半徑為r1、r2的毛細(xì)管個(gè)數(shù)。

［1］秦積舜，李愛芬.油層物理學(xué)［M］.東營(yíng)：石油大學(xué)出版社，2001：147-149.

［2］楊勝來，魏俊之.油層物理學(xué)［M］.北京：石油工業(yè)出版社，2004：36-137.

［3］郭平，徐永高，陳召佑，等.對(duì)低滲氣藏滲流機(jī)理實(shí)驗(yàn)研究的新認(rèn)識(shí)［J］.天然氣工業(yè)，2007，27（7）：86-88.

［4］郝斐，程林松，李春蘭，等.特低滲透油藏啟動(dòng)壓力梯度研究［J］.西南石油學(xué)院學(xué)報(bào)，2006，28（6）：29-32.

［5］呂成遠(yuǎn)，王建，孫志剛.低滲透砂巖油藏滲流啟動(dòng)壓力梯度實(shí)驗(yàn)研究［J］.石油勘探與開發(fā)，2002，29（2）：86-88.

［6］王道成，李閩，喬國(guó)安，等.天然氣啟動(dòng)壓力梯度實(shí)驗(yàn)研究［J］.鉆采工藝，2007，30（5）：53-55.

［7］吳凡.氣體滲流特征及啟動(dòng)壓力規(guī)律的研究［J］.天然氣工業(yè)，2001，21（1）：82-84.

［8］李忠興.特低滲透油藏啟動(dòng)壓力梯度的求解方法及應(yīng)用［J］.石油勘探與開發(fā)，2004，31（3）：107-109.

［9］許建紅，程林松，周穎，等.一種求解低滲透油藏啟動(dòng)壓力梯度的新方法［J］.石油勘探與開發(fā)，2007，34（5）：594-597.

［10］李士倫，孫雷，杜建芬，等.低滲致密氣藏、凝析氣藏開發(fā)難點(diǎn)與對(duì)策［J］.新疆石油地質(zhì)，2004，25（2）：156.

［11］李傳亮，楊永全.啟動(dòng)壓力其實(shí)并不存在［J］.西南石油大學(xué)學(xué)報(bào)，2008，30（3）：167-170.

［12］周惠忠，王利群.兩維油藏物理模擬裝置［J］.清華大學(xué)學(xué)報(bào)（自然科學(xué)版），1994，34（3）：74-82.

［13］楊珂，徐守余.微觀剩余油實(shí)驗(yàn)方法研究［J］.斷塊油氣田，2009，16（4）：75-77.

The mechanical analysis and application of the flow launch in the capillary tube

XUAN Jian
（Exploration and Development Research Institute of Shengli Oilfield，Dongying Shandong 257000，China）

In the capillary tube,the equation of the single-phase instability flow is established on the basis of the Poiseuille equation and the Newton second law.The change of displacement and velocity with time is obtained.And the relationship of inertia force and seepage flow resistance is calculated.The capillary tube with flow is applied in the analysis of the threshold pressure gradient.

capillary；threshold pressure gradient；flow

TE312

A

1673-5285（2017）02-0077-03

10.3969/j.issn.1673-5285.2017.02.018

2016－12-08

玄建，男（1982-），勝利油田勘探開發(fā)研究院工程師，碩士，研究方向?yàn)橛蜌馓镩_發(fā)，郵箱：xuanjian0202@163.com。