張憲法，王繼偉，魏 馳，楊旭剛，李俊儒，陳付真

1 研究背景

目前我國(guó)部分油田進(jìn)入中高含水期，由水竄導(dǎo)致的含水率快速上升、采出程度低等問(wèn)題突出，因此調(diào)剖堵水作為一項(xiàng)經(jīng)濟(jì)有效的提高采收率手段得到了廣泛應(yīng)用。常用的堵水方法包括機(jī)械堵水和化學(xué)堵水。機(jī)械堵水工具主要包含堵水管柱以及相應(yīng)的封隔器，其中堵水管柱有可鉆式、整體式、懸掛式以及平衡式四種類型。油田實(shí)際開采進(jìn)入到中高含水期之后，機(jī)械堵水技術(shù)逐漸向可調(diào)層技術(shù)以及液壓可調(diào)層技術(shù)兩方面發(fā)展[1-2]?；瘜W(xué)堵水是通過(guò)將化學(xué)堵劑注入出水層位，利用堵劑的化學(xué)性質(zhì)或化學(xué)反應(yīng)物的物理化學(xué)性質(zhì)達(dá)到封堵地層出水孔道、降低油井綜合含水率的目的[3]。目前常用的化學(xué)堵劑主要有水泥類堵劑、樹脂類堵劑、無(wú)機(jī)鹽沉淀類堵劑、水溶性聚合物凍膠類堵劑、顆粒類堵劑以及泡沫類堵劑等[4]。

目前油田化學(xué)堵水大量使用由聚合物和交聯(lián)劑交聯(lián)形成的絮凝體作為堵劑，但絮凝體不耐沖刷、容易被突破等問(wèn)題嚴(yán)重影響其堵水效果[5-9]。通過(guò)在絮凝體中加入彈性纖維，形成纖維-絮凝體二元堵水體系，為解決上述問(wèn)題提供了新的思路。在該二元體系中，纖維在孔隙中與壁面接觸，搭建形成二元體系骨架，受力后產(chǎn)生彈性反力；絮凝體充填骨架體系，與孔隙接觸產(chǎn)生摩擦力、粘滯力，同時(shí)與纖維協(xié)同作用，提高二元體系強(qiáng)度。趙修太等人進(jìn)行了纖維-凍膠的聯(lián)合堵水實(shí)驗(yàn)，對(duì)比分析了凍膠體系與纖維-凍膠二元體系在不同溫度與不同pH值影響下的堵水效果，結(jié)果顯示兩者在同等 pH值和溫度條件下，纖維-凍膠二元體系在封堵率、突破壓力、孔隙體積減小程度上均優(yōu)于凍膠體系，同時(shí)實(shí)驗(yàn)結(jié)果也表明二元堵水體系封堵率比絮凝體封堵率提高31.68%[10]。

目前纖維-絮凝體二元堵水體系研究主要集中在室內(nèi)實(shí)驗(yàn)方面，得到的結(jié)論主要為宏觀統(tǒng)計(jì)規(guī)律，微觀機(jī)理尚不明確，纖維和絮凝體聯(lián)合作用的理論研究尚不完善。為此，本文開展了纖維-絮凝體二元堵水體系理論研究。首先借鑒材料力學(xué)的純彎曲梁變形特性[11]，并結(jié)合泰勒公式變換以及概率論中常用概率模型等數(shù)學(xué)工具[12]，構(gòu)建了纖維在堵水時(shí)的阻力計(jì)算模型；其次基于非牛頓流體的流體力學(xué)理論，給出了考慮屈服應(yīng)力的絮凝體堵水阻力模型[13]；與此同時(shí)，根據(jù)相界面法及摩擦阻力計(jì)算方法[14]，給出了通過(guò)相界面表示的絮凝體堵水阻力的計(jì)算模

2 纖維阻力計(jì)算模型

2.1 纖維單體阻力計(jì)算

型；最后將纖維和絮凝體的阻力計(jì)算模型進(jìn)行耦合，建立了一維纖維和絮凝體二元體系聯(lián)合堵水時(shí)的堵水阻力模型，并根據(jù)實(shí)驗(yàn)條件的不同進(jìn)行了優(yōu)化。在纖維阻力計(jì)算模型中，首先研究單根纖維在等直徑毛細(xì)管中所產(chǎn)生的阻力。假設(shè)長(zhǎng)度為L(zhǎng)的單根纖維垂直置于等直徑毛細(xì)管中，毛細(xì)管直徑為 D，外力為F，F(xiàn)為均勻力場(chǎng)，纖維變形如圖1所示。

圖1 纖維受力計(jì)算模型

根據(jù)力的平衡原理可知纖維所產(chǎn)生的阻力fF即為纖維所受外力：fF=F×D；根據(jù)材料的純彎曲變形撓曲線方程關(guān)系可知：

由于纖維在毛細(xì)管中的變形近似為弧形，經(jīng)過(guò)幾何變換得：

將sinα用泰勒公式展開，經(jīng)變換得：

慣性矩為：

將公式（4）和（5）代入式（2）得纖維彈性阻力計(jì)算模型如下：

式中：F為纖維所受均勻力場(chǎng)，N/m；fF為纖維彈性阻力，N；α為圓心半角，Rad；D為毛細(xì)管直徑，m；L為纖維長(zhǎng)度，m；E為纖維彈性模量，Pa；d為纖維截面圓的直徑，m；I為纖維的慣性矩，m4。

根據(jù)公式（7）可知纖維的阻力受纖維彈性模量、纖維直徑、纖維長(zhǎng)度以及孔隙直徑影響[15]。假設(shè)纖維的彈性模量E=106Pa、孔隙直徑D=1cm，對(duì)纖維彈性阻力與纖維長(zhǎng)度和纖維直徑的關(guān)系進(jìn)行研究，得出纖維屬性參數(shù)對(duì)阻力影響規(guī)律如圖2所示。

由圖2可知，在纖維長(zhǎng)度與孔隙直徑之比一定時(shí)，纖維彈性阻力隨纖維直徑的增大而增長(zhǎng)，所以為提高堵水效果，應(yīng)盡可能選用直徑較大的纖維。纖維彈性阻力的增長(zhǎng)幅度隨纖維長(zhǎng)度與孔隙直徑之比的增大而增大。纖維長(zhǎng)度與孔隙直徑之比在 1～2之間時(shí)，彈性阻力增長(zhǎng)幅度變化比較明顯；纖維長(zhǎng)度與孔隙直徑之比大于3后彈性阻力增長(zhǎng)幅度不明顯，故纖維長(zhǎng)度與孔隙直徑之比存在最優(yōu)值。

圖2 纖維屬性參數(shù)對(duì)阻力的影響

圖3可得出，纖維彈性阻力隨纖維長(zhǎng)度與孔隙直徑比值的增加而增大，當(dāng)比值達(dá)到3時(shí)上升趨勢(shì)變緩，曲線逐漸平穩(wěn)。所以可以推斷孔隙直徑一定時(shí)，對(duì)應(yīng)纖維阻力纖維長(zhǎng)度有最佳值。從圖中得出結(jié)論，纖維長(zhǎng)度與孔隙直徑之比最優(yōu)值為3。

圖3 纖維屬性參數(shù)對(duì)阻力的影響

2.2 纖維概率分布模型

當(dāng)多根纖維同時(shí)注入孔隙中時(shí)，纖維在孔隙中分布具有不確定性，纖維與孔隙壁面接觸角度不同時(shí)其對(duì)有效阻力的貢獻(xiàn)程度也不同。因此，在計(jì)算纖維有效阻力時(shí)需要考慮纖維在孔隙中的分布狀態(tài)和分布規(guī)律。

假設(shè)單位時(shí)間內(nèi)有n根纖維注入注水井，若有n～ m根滯留在儲(chǔ)層中，則有效纖維比例為：

對(duì)于滯留在孔隙中的纖維，其分布主要為圖 4所示的三種形態(tài)。纖維兩端連線與孔壁夾角為γ，已知γ的取值范圍為：

由圖5可知，γ在0.5π附近時(shí)纖維易在孔隙中滯留，即纖維與孔隙壁面角度呈γ=0.5π時(shí)概率最大，在 角 度 γ1= a rcsin（ D / L）或時(shí)概率最小，近似為0。這一分布規(guī)律與正態(tài)分布相吻合，因此假設(shè)纖維在孔隙中的概率分布符合正態(tài)分布。

圖5 纖維分布概率

正態(tài)分布的概率分布密度為：

式中：μ為常數(shù)，等于0.5π；σ為常數(shù)，纖維以γ=0.5π角度滯留在孔隙中的概率。

根據(jù)歸一化條件，即：

通過(guò)式（11）可求得σ的值。

長(zhǎng)為L(zhǎng)的纖維傾斜在孔隙中產(chǎn)生的阻力，相當(dāng)于纖維投影在豎直方向上纖維長(zhǎng)度在孔隙中產(chǎn)生的阻力，即相當(dāng)于長(zhǎng)為 sinL γ的纖維垂直與孔隙壁面時(shí)在孔隙中產(chǎn)生的阻力?；诶w維彈性阻力計(jì)算模型，以任意形狀滯留在毛管中的纖維產(chǎn)生的阻力為：

對(duì)于注入的n根纖維，產(chǎn)生的總阻力tfF為：

上式即為假設(shè)纖維在孔隙中的分布形態(tài)符合正態(tài)分布時(shí)，纖維總阻力計(jì)算模型。

3 絮凝體堵水阻力模型

3.1 屈服應(yīng)力法

絮凝體是由一種線型高分子經(jīng)過(guò)交聯(lián)后產(chǎn)生的具有網(wǎng)狀結(jié)構(gòu)的化學(xué)物質(zhì)，具有一定的連續(xù)性和穩(wěn)定性。在調(diào)剖堵水過(guò)程中，絮凝體能夠滯留在孔隙中有效封堵水竄通道。從流變模式上說(shuō)，絮凝體是一種非牛頓流體（塑性流體），它具有與牛頓流體不同的流變特性，絮凝體流動(dòng)存在啟動(dòng)屈服應(yīng)力[16]，該啟動(dòng)屈服應(yīng)力也就是絮凝體能夠承受的最大應(yīng)力。當(dāng)流體壓力超過(guò)啟動(dòng)屈服應(yīng)力時(shí)，絮凝體產(chǎn)生形變開始流動(dòng)，因而失去封堵效果。絮凝體產(chǎn)生的堵水阻力為：

式中：jF為絮凝體產(chǎn)生的阻力，N；sτ為絮凝體塑性流體的啟動(dòng)屈服應(yīng)力，可由實(shí)驗(yàn)確定，Pa；A為絮凝體與孔隙的接觸面積，可根據(jù)絮凝體的體積和孔隙的半徑求得，m2。

這種計(jì)算絮凝體阻力的方法，是根據(jù)實(shí)驗(yàn)現(xiàn)象推導(dǎo)出來(lái)的，是從宏觀角度考慮問(wèn)題得到的絮凝體阻力計(jì)算模型。

3.2 相界面法

注入水與絮凝體之間存在相界面，相界面上的分子受兩相的引力大小不同，產(chǎn)生了界面張力。對(duì)于相界面的曲界面來(lái)說(shuō)，其兩側(cè)的壓力不同，此曲界面壓差是堵水阻力的一部分。

根據(jù)界面化學(xué)理論，水與絮凝體產(chǎn)生的曲界面壓差為：

式中：cP為曲界面壓差，Pa；σ為水與絮凝體的界面張力，N/m；θ為水對(duì)絮凝體的潤(rùn)濕角，Rad；r為毛細(xì)管孔隙的半徑，m。

已知當(dāng)0θ=即cos1θ=時(shí)，cP取得最大值，為：

此時(shí)，由相界面產(chǎn)生的絮凝體阻力最大，為：

式中：pA是毛細(xì)管孔隙的橫截面積，m2。

除了相界面壓差外，水與絮凝體相界面還存在靜摩擦力，在絮凝體未被突破前，絮凝體與孔隙壁之間的靜摩擦力大小為：

式中：N為正壓力，N；η為摩擦系數(shù)。

式（18）中正壓力N可根據(jù)水與絮凝體相界面張力求得：

將式（19）代入式（18）可得絮凝體產(chǎn)生的阻力為：

當(dāng)sin1=θ時(shí)，絮凝體產(chǎn)生的阻力最大，為：

將上述絮凝體曲界面壓差和靜摩擦力進(jìn)行耦合，可得絮凝體產(chǎn)生的阻力為：

相界面法計(jì)算絮凝體堵水阻力，是從微觀角度考慮絮凝體阻力問(wèn)題，與屈服應(yīng)力法并不矛盾，這兩者可以相互補(bǔ)充，同時(shí)也可以相互驗(yàn)證。

4 二元體系堵水阻力模型

纖維與絮凝體在儲(chǔ)層中產(chǎn)生復(fù)合堵水網(wǎng)絡(luò)體系，二元體系在孔隙中的分布形態(tài)如圖6所示。

將纖維阻力模型與絮凝體阻力模型進(jìn)行耦合得到纖維-絮凝體復(fù)合網(wǎng)絡(luò)的阻力計(jì)算模型為：

式中：a為纖維疊加阻力系數(shù)；b為絮凝體疊加阻力系數(shù)。a、b可根據(jù)實(shí)際條件的不同，通過(guò)室內(nèi)實(shí)驗(yàn)進(jìn)行求得。

圖6 纖維-絮凝體復(fù)合堵水網(wǎng)絡(luò)體系示意圖

最后，通過(guò)VB軟件進(jìn)行了計(jì)算模型的編譯，實(shí)現(xiàn)了單根纖維在堵水時(shí)的阻力、絮凝體單獨(dú)堵水時(shí)的阻力以及纖維-絮凝體二元網(wǎng)絡(luò)堵水體系堵水阻力的計(jì)算，從而提高了纖維-絮凝體二元體系優(yōu)化的效率。

5 結(jié)論

（1） 建立了纖維在一維二元網(wǎng)絡(luò)堵水體系中的阻力計(jì)算模型，分析了纖維屬性參數(shù)對(duì)纖維阻力影響規(guī)律，得到了纖維長(zhǎng)度與孔隙直徑之比最優(yōu)值為3。

（2） 基于屈服應(yīng)力法和相界面法分別建立了絮凝體在一維二元網(wǎng)絡(luò)堵水體系中的阻力計(jì)算模型，兩模型可以相互補(bǔ)充，同時(shí)也可以相互驗(yàn)證。

（3） 研究了纖維在多孔介質(zhì)中的分布規(guī)律，建立了纖維分布概率模型。在此基礎(chǔ)上實(shí)現(xiàn)了纖維與絮凝體阻力模型的耦合，得到纖維-絮凝體一維二元網(wǎng)絡(luò)聯(lián)合堵水體系阻力計(jì)算模型。

參考文獻(xiàn)：

[1] 章秀平.堵水調(diào)剖技術(shù)在油田發(fā)展中的應(yīng)用[J].化工管理，2016，29（6）：128.

[2] 尤建華，王衛(wèi)真.油田機(jī)械堵水技術(shù)發(fā)展研究[J].科技風(fēng)，2016，29（15）：141.

[3] 吳玉柱，嚴(yán)慧呈，楊小東.油田化學(xué)堵水工藝技術(shù)[J].化學(xué)工程與裝備，2016，45（7）：88–89.

[4] 柯耀斌.論油田化學(xué)中堵水調(diào)剖的開發(fā)及應(yīng)用[J].中國(guó)石油和化工標(biāo)準(zhǔn)與質(zhì)量，2017，37（3）：62–63.

[5] 劉文龍，張杰，趙英杰，等.油田化學(xué)耐鹽堵劑研究進(jìn)展[J].內(nèi)蒙古石油化工，2011，21（18）：30–34.

[6] 劉玉章.聚合物驅(qū)油提高油田采收率技術(shù)[M].北京：石油工業(yè)出版社，2006.

[7] D D WHITNEY，徐文芹.聚合物凝膠體系的堵水試驗(yàn)[J].國(guó)外油田工程，1999，14（7）：25–29.

[8] 周彬，王慧玲，樊理山，等.含膠類纖維成膜裝置的設(shè)計(jì)及實(shí)驗(yàn)[J].上海紡織科技，2013，40（10）：21–22.

[9] 任建曉，林濤，殷學(xué)風(fēng)，等.細(xì)小纖維與CPAM共絮凝體的性質(zhì)[J].湖南造紙，2013，32（4）：3–11.

[10] 趙修太，董林燕.纖維–聚合物凍膠復(fù)合堵水體系研究[J].精細(xì)石油化工，2013，29（4）：18–21.

[11] 范欽珊，唐靜靜.工程力學(xué)（靜力學(xué)和材料力學(xué)）[M].北京：高等教育出版社，2007：7–10.

[12] 王清河，常照光，曹曉敏，等.隨機(jī)數(shù)據(jù)處理方法[M].東營(yíng)：中國(guó)石油大學(xué)出版社，2011.

[13] 祝連慶，董明利，丁裕棟.非牛頓流體流變學(xué)特性測(cè)試與建模[N].北京機(jī)械工業(yè)學(xué)院學(xué)報(bào)，1997，12（1）：21–24.

[14] 曾作祥，孫莉.界面現(xiàn)象[M].上海：華東理工大學(xué)出版社，2016.

[15] 任志華，于永玲.影響亞麻纖維彈性的因素及其提高方法[J].大連輕工業(yè)學(xué)院學(xué)報(bào)，2002，22（3）：228–231.

[16] 張國(guó)忠，劉剛.大慶膠凝原油啟動(dòng)屈服應(yīng)力研究[J].石油大學(xué)學(xué)報(bào)，2005，47（6）：91–93.