魏臣興，練章華，丁亮亮，郭衍茹，韓 濤

（西南石油大學(xué)“油氣藏地質(zhì)及開發(fā)工程”國(guó)家重點(diǎn)實(shí)驗(yàn)室）

分支井滲流-應(yīng)力耦合場(chǎng)分析

魏臣興，練章華，丁亮亮，郭衍茹，韓 濤

（西南石油大學(xué)“油氣藏地質(zhì)及開發(fā)工程”國(guó)家重點(diǎn)實(shí)驗(yàn)室）

在多孔介質(zhì)的滲流過(guò)程中，存在著極強(qiáng)的流固耦合作用。針對(duì)廣泛使用的分支井型，在利用COMSOL Multiphysics軟件時(shí)充分考慮流固耦合作用，分析耦合滲流場(chǎng)和應(yīng)力場(chǎng)的分布特點(diǎn)，同時(shí)通過(guò)運(yùn)用三維庫(kù)倫準(zhǔn)則計(jì)算失效系數(shù)fail值來(lái)判斷井壁巖石的危險(xiǎn)截面，并進(jìn)一步計(jì)算出不同造斜率情況下失效系數(shù)fail的最小值和失效系數(shù)fail≤-2.5 MPa區(qū)域的體積積分值，從而研究采用不同造斜率側(cè)鉆分支井眼對(duì)井壁失效區(qū)域大小的影響。結(jié)果表明，在單滲流場(chǎng)分析中采用定壓縮系數(shù)方法得到的孔隙度其變化微乎其微，而耦合滲流場(chǎng)中的孔隙度、滲透率在井底附近區(qū)域變化明顯。因此，考慮流固耦合作用下的滲流場(chǎng)能更真實(shí)地反映地下滲流情況。在耦合應(yīng)力場(chǎng)分析中，Von Mises等效應(yīng)力的最大值和失效系數(shù)fail最小值均出現(xiàn)在側(cè)鉆分支井窗口處，即窗口處井壁巖石為危險(xiǎn)截面，進(jìn)而得出在取造斜率為15°/30 m進(jìn)行側(cè)鉆分支井眼的情況下井壁巖石失效范圍最小的結(jié)論，為分支井的設(shè)計(jì)和施工提供了理論依據(jù)。

分支井；滲流場(chǎng)；應(yīng)力場(chǎng)；流固耦合；數(shù)值模擬；井壁穩(wěn)定

0 引言

油藏?cái)?shù)值模擬分析方法經(jīng)歷了從單滲流場(chǎng)研究到耦合滲流場(chǎng)研究的歷程。單滲流場(chǎng)研究只考慮滲流場(chǎng)而忽略了應(yīng)力場(chǎng)的作用，隨著計(jì)算機(jī)技術(shù)的發(fā)展，耦合研究特別是流固耦合作用越來(lái)越受到重視。經(jīng)流固耦合研究發(fā)現(xiàn)，在油藏的開采過(guò)程中，地層中各種力處于一種動(dòng)態(tài)平衡狀態(tài)，地層孔隙壓力發(fā)生變化，引起巖石骨架有效應(yīng)力的改變，巖石骨架發(fā)生變形，進(jìn)而導(dǎo)致了孔隙度、滲透率等物性參數(shù)的改變；同時(shí)，滲流場(chǎng)中物性參數(shù)的變化又會(huì)影響地層孔隙壓力的重新分布，最終將影響地層流體的滲流和開采動(dòng)態(tài)。有研究表明，油藏開采過(guò)程中存在極強(qiáng)的流固耦合作用，是不可忽略的［1］。

前人研究應(yīng)力場(chǎng)的作用時(shí)，雖通過(guò)采用隨定值壓縮系數(shù)變化的孔隙度來(lái)表征應(yīng)力場(chǎng)對(duì)滲流場(chǎng)的影響［2］，但本文計(jì)算結(jié)果表明采用此方法孔隙度變化甚微，近似為常數(shù)，而實(shí)際地層中的孔隙度和滲透率應(yīng)為地層有效應(yīng)力的函數(shù)。筆者采用廣泛使用的分支井井型［3］，并充分考慮滲流場(chǎng)與應(yīng)力場(chǎng)耦合作用下滲流參數(shù)將隨巖石應(yīng)力狀態(tài)不同而發(fā)生改變，來(lái)分析側(cè)鉆分支井開窗處的滲流場(chǎng)和應(yīng)力場(chǎng)特點(diǎn)，同時(shí)應(yīng)用三維庫(kù)倫準(zhǔn)則，判斷井壁巖石易于失穩(wěn)的危險(xiǎn)截面，以及從主井眼側(cè)鉆分支井眼的造斜率對(duì)井壁危險(xiǎn)截面的影響。

1 數(shù)學(xué)模型

1.1 滲流場(chǎng)數(shù)學(xué)模型

為了研究流固耦合滲流場(chǎng)的特點(diǎn)以及耦合場(chǎng)中井壁巖石的應(yīng)力分布情況，假設(shè)巖石孔隙為不可壓縮單相流體所飽和，整個(gè)滲流過(guò)程是穩(wěn)態(tài)的且等溫的，建立起流固耦合滲流的數(shù)學(xué)模型為［1，4-6］

為了研究側(cè)鉆分支井眼處的流場(chǎng)，根據(jù)模型的具體情況給定不同的邊界條件：

供給邊界

井底

對(duì)稱邊界和上、下邊界

1.2 應(yīng)力場(chǎng)數(shù)學(xué)模型

作用于地層巖石上的體積力主要是重力和滲透體積力（?p/?x，?p/?y，?p/?z）［2］。 這些力在壓降沒(méi)有傳遞到供給邊界期間是隨時(shí)間發(fā)生變化的，但壓降傳遞到供給邊界后則慢慢趨于穩(wěn)定。分析過(guò)程中描述地層巖石骨架變形的基本方程為

根據(jù)Shaybah油田所鉆裸眼分支井情況［3］，模型邊界條件如下：

供給邊界

對(duì)稱邊界

上邊界

下邊界

井底

1.3 庫(kù)倫判斷準(zhǔn)則

在滲流-應(yīng)力耦合場(chǎng)中，采用3D庫(kù)倫準(zhǔn)則，充分考慮3個(gè)主應(yīng)力和地層孔隙壓力的影響，來(lái)判斷危險(xiǎn)截面，見(jiàn)式（13）

其中

利用式（13）得到fail值，進(jìn)而判斷井壁巖石的狀態(tài)，其中fail為正值時(shí)表明巖石處于穩(wěn)定狀態(tài)；fail為負(fù)值時(shí)表明巖石易于失效，且負(fù)值越大，失效的可能性越大。

2 耦合場(chǎng)分析

COMSOL Multiphysics是一款基于有限元理論的直接以偏微分方程為研究對(duì)象的多物理場(chǎng)耦合軟件，具有求解耦合問(wèn)題的強(qiáng)大能力。將式（1）與式（7）2種微分方程在該軟件中進(jìn)行耦合計(jì)算，分析垂深1500 m處側(cè)鉆分支井眼區(qū)域在生產(chǎn)壓差4.7 MPa的條件下的滲流場(chǎng)分布（圖1）。圖1中流線顯示滲流方向，流線上顏色代表壓力的分布情況，切面上顏色代表壓力梯度的變化，在地層流體向井底滲流的過(guò)程中，壓力從供給邊界到井底的變化表現(xiàn)為明顯的“壓降漏斗”形式，這一趨勢(shì)在圖2的切線圖中表現(xiàn)得更為明顯。圖2是以側(cè)鉆分支井眼的窗口處為0 m，向井深方向，間隔20 m，取出井眼剖面上橫向的壓力分布情況，后續(xù)的切線圖也是在相同位置處取值。圖1中切面圖上壓力梯度的分布和圖2壓力切線圖中的“壓降漏斗”表明壓力主要消耗在井底附近，這與數(shù)值求解得到的結(jié)論一致［2］，說(shuō)明此模型的研究方法是正確合理的。

圖1 20°/30 m分支井滲流場(chǎng)分布圖Fig.1 The seepage field distribution of 20°/30 m

圖2 分支井壓力切線圖Fig.2 Pressure distribution of lateral wells

在圖2中，沿井深方向距側(cè)鉆點(diǎn)20 m，40 m和60 m處，由于主井眼和分支井眼橫向距離的加大及井間干擾作用的減弱，兩井眼之間的壓力也在逐漸恢復(fù)；在分支井中存在分支井眼的一側(cè)，由于供給邊界的限制，分支井眼的“壓降漏斗”更為陡峭。在不同深度條件下，分支井中不存在分支井眼的一側(cè)，由于距離分支井眼較遠(yuǎn)，其橫向分布的壓力受井間干擾引起的壓降不明顯，即側(cè)鉆分支井對(duì)驅(qū)動(dòng)該區(qū)域的地層流體意義不大，分支井眼僅對(duì)側(cè)鉆分支井眼的一側(cè)作用效果較好。因此，在地層壓力衰竭的油田或薄油層開展老井側(cè)鉆，最好是采用雙分支井或多分支井，以提高油井控制區(qū)域的壓力波及能力，從而實(shí)現(xiàn)增產(chǎn)增注。

在耦合滲流場(chǎng)分析中，巖石骨架顆粒在滲透體積力和重力的作用下其力學(xué)行為發(fā)生改變，在滲流場(chǎng)中反映為地層的孔隙度和滲透率隨著應(yīng)力場(chǎng)作用程度的不同而發(fā)生相應(yīng)的變化［7-11］，見(jiàn)圖3、圖4所示孔隙度、滲透率變化曲線（在分析中所采用初始孔隙度為25%，初始滲透率為10 mD）。結(jié)合圖2中的壓力分布情況對(duì)圖3、圖4中的物性參數(shù)分布趨勢(shì)進(jìn)行分析。在供給邊界附近，由于地層孔隙壓力變化不大，因而所引起的巖石變形相對(duì)也比較小，在達(dá)到穩(wěn)定滲流狀態(tài)時(shí)，巖石的孔隙度和滲透率變化不大，與初始孔隙度和滲透率相當(dāng)；而在靠近井底附近區(qū)域，孔隙壓力急劇變化，所引起的巖石變形也相當(dāng)顯著，對(duì)巖石的物性參數(shù)影響十分明顯?？偟膩?lái)說(shuō)，在分析區(qū)域內(nèi)，孔隙度、滲透率變化趨勢(shì)與“壓降漏斗”相對(duì)應(yīng)。而采用不考慮應(yīng)力場(chǎng)影響的單滲流場(chǎng)分析時(shí)，在分析程序中嵌入孔隙度隨定值的孔隙壓縮系數(shù)變化的表達(dá)式如下

經(jīng)計(jì)算可知，孔隙度值的改變微乎其微，基本上為初始常數(shù)，見(jiàn)圖3中虛線所示。將其與耦合滲流場(chǎng)分析得到的參數(shù)進(jìn)行對(duì)比，可見(jiàn)耦合滲流場(chǎng)分析得到的孔隙度、滲透率要小于單滲流場(chǎng)分析得到的參數(shù)，特別是在井底附近。由于孔隙度、滲透率的減小，整個(gè)地層的壓力分布將大于僅考慮滲流場(chǎng)時(shí)的壓力分布，且滲流速度也相應(yīng)地減小，最終將影響整個(gè)油藏的開發(fā)動(dòng)態(tài)［12］。

圖3 耦合滲流場(chǎng)中孔隙度分布圖Fig.3 Porosity distribution in coupling seepage field

圖4 耦合滲流場(chǎng)中滲透率分布圖Fig.4 Permeability distribution in coupling seepage field

圖5 圖1所標(biāo)示區(qū)域1處應(yīng)力和失效系數(shù)局部分布Fig.5 Local distribution of stress and fail coefficient of region 1 showed in Fig.1

圖5為圖1中所標(biāo)示區(qū)域1處應(yīng)力場(chǎng)和失效系數(shù)fail的局部分布圖。根據(jù)材料的第四強(qiáng)度理論計(jì)算的Von Mises等效應(yīng)力最大值22.5 MPa出現(xiàn)在窗口附近的井壁處；同樣失效系數(shù)fail值在井壁附近負(fù)值較為集中，其最小值處也就是最易發(fā)生井壁不穩(wěn)定的截面，同時(shí)也正是在側(cè)鉆分支井眼的窗口處。在不考慮井壁巖石內(nèi)化學(xué)作用的前提下，單從力學(xué)方面分析，側(cè)鉆分支井眼窗口附近的井壁巖石極易垮塌。因此除非巖石結(jié)構(gòu)極為穩(wěn)定，否則在側(cè)鉆分支井中，應(yīng)對(duì)側(cè)鉆分支井眼的窗口附近加強(qiáng)完井質(zhì)量，將其周圍封固好，最好將分支井的完井級(jí)別提高到3級(jí)以上，以便可用井下連接工具實(shí)現(xiàn)分支井眼和主井眼的機(jī)械連接，從而提高不穩(wěn)定地層的連接強(qiáng)度，防止窗口附近區(qū)域于生產(chǎn)過(guò)程中在滲流場(chǎng)和應(yīng)力場(chǎng)耦合作用下發(fā)生事故，提高安全生產(chǎn)時(shí)間，延長(zhǎng)油井壽命［13］。

3 造斜率對(duì)井壁穩(wěn)定性的影響

在研究主井眼中側(cè)鉆分支井眼的造斜率對(duì)井壁穩(wěn)定性的影響時(shí)，考慮到分支井眼造斜率一般要求在（10°～20°）/30 m，并結(jié)合國(guó)內(nèi)實(shí)施的分支井情況，分別取造斜率為 5°/30 m，10°/30 m，15°/30 m，20°/30 m 及 25°/30 m 來(lái)加以研究［14-16］。

造斜率對(duì)失效系數(shù)fail的影響如圖6所示。隨著造斜率的增加，fail系數(shù)最小值是曲折變化的，在模型考慮的造斜率范圍內(nèi)，造斜率為15°/30 m時(shí)取得較大值。也就是在地層因素相同的條件下，其它造斜率相對(duì)于15°/30 m時(shí)的造斜率更易出現(xiàn)井壁不穩(wěn)定問(wèn)題。取fail系數(shù)≤-2.5 MPa的嚴(yán)重失效區(qū)域進(jìn)行體積積分，可以得到各個(gè)造斜率條件下井壁不穩(wěn)定的范圍大小。如圖6所示，在造斜率為15°/30 m時(shí)的體積積分取得較小值。因此，考慮到現(xiàn)有側(cè)鉆分支井的工具設(shè)備以及井眼軌跡控制能力，在滿足采油工藝的前提下，對(duì)于主井眼中側(cè)鉆分支井最好采取15°/30 m的造斜率為宜。

圖6 造斜率對(duì)fail系數(shù)和fail≤-2.5的區(qū)域體積的影響Fig.6 The impact of build-up rate on fail coefficient and volume integration in the condition of fail≤-2.5

4 結(jié)論

（1）耦合滲流場(chǎng)的壓力分布與常規(guī)數(shù)值求解結(jié)果分布趨勢(shì)大致相似，但耦合滲流場(chǎng)中的孔隙度和滲透率是變化的，更能反映真實(shí)地層滲流。

（2）在應(yīng)力場(chǎng)作用下，地層巖石的孔隙度、滲透率急劇減小，進(jìn)而影響滲流場(chǎng)和油藏的開采動(dòng)態(tài)，其中耦合滲流場(chǎng)各點(diǎn)處的壓力大于單滲流場(chǎng)分析得到的壓力，耦合滲流場(chǎng)的速度小于單滲流場(chǎng)的滲流速度。

（3）分支井窗口處的井壁最不穩(wěn)定，最好能將分支井的完井級(jí)別提高到3級(jí)以上，以便提高不穩(wěn)定地層的連接強(qiáng)度，避免井下事故的發(fā)生。

（4）主井眼中側(cè)鉆分支井眼的造斜率不同，對(duì)井壁的危險(xiǎn)程度影響較大?？紤]到現(xiàn)階段的鉆井能力，側(cè)鉆分支井眼時(shí)取造斜率為15°/30 m最為合適。

符號(hào)說(shuō)明：

p——地層孔隙壓力，MPa；

pr——供給邊界壓力，MPa；

Δp——地層孔隙壓力變化值，MPa；

pw——井底壓力，MPa；

p0——上覆地層壓力，MPa；

p1——模型頂界處骨架應(yīng)力，MPa；

p2——井底巖石所承受的外壓，MPa；

η——?jiǎng)恿︷ざ龋琺Pa·s；

D——相對(duì)于基準(zhǔn)面，計(jì)算區(qū)域的位置水頭，m；

φ——孔隙度，%；

k——滲透率，mD；

φ0——初始孔隙度，%；

k0——初始滲透率，mD；

ε——巖石的體積應(yīng)變；

Cp——巖石壓縮系數(shù)，取1×10-6，MPa-1；

ρf——地層流體密度，kg/m3；

n——邊界的單位法向向量；

u——模型中x方向位移，m；

v——模型中y方向位移，m；

w——模型中z方向位移，m；

u——位移矢量u，v和w；

F——力矢量；

c——系數(shù)；

α——有效應(yīng)力系數(shù)；

θ——巖石內(nèi)摩擦角，（°）；

S0——巖石內(nèi)聚力，MPa；

σ1——第一主應(yīng)力，MPa；

σ2——第二主應(yīng)力，MPa；

σ3——第三主應(yīng)力，MPa；

fail——判斷井壁是否失效系數(shù)，MPa。

［1］胡永樂(lè)，冉啟全，孫建平.流固耦合油藏?cái)?shù)值模擬［M］.北京：石油工業(yè)出版社，2007：3-98.

［2］李曉平.地下油氣滲流力學(xué)［M］.北京：石油工業(yè)出版社，2008：35-59.

［3］Dossary A S，Saudi Aramco，Mahgoub A A.Challenges and Achievements of Drilling Maximum Reservoir Contact（MRC） Wells in Shaybah Field［R］.SPE 85307，2003.

［4］李允.油藏模擬［M］.東營(yíng)：石油大學(xué)出版社，1999：37-49.

［5］趙延林，曹平，趙陽(yáng)升，等.雙重介質(zhì)溫度場(chǎng)-滲流場(chǎng)-應(yīng)力場(chǎng)耦合模型及三維數(shù)值研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2007，26（增2）：4024-4031.

［6］盛金昌.多孔介質(zhì)流-固-熱三場(chǎng)全耦合數(shù)學(xué)模型及數(shù)值模擬［J］.巖石力學(xué)與工程學(xué)報(bào)，2006，25（增 1）：3028-3033.

［7］樊懷才，李曉平，竇天財(cái)，等.應(yīng)力敏感效應(yīng)的氣井流量動(dòng)態(tài)特征研究［J］.巖性油氣藏，2010，22（4）：130-134.

［8］張宇焜，汪偉英，周江江.注水壓力對(duì)低滲透儲(chǔ)層滲流特征的影響［J］.巖性油氣藏，2010，22（2）：120-122.

［9］傅春梅，唐海，鄒一鋒，等.應(yīng)力敏感對(duì)蘇里格致密低滲氣井廢棄壓力及采收率的影響研究［J］.巖性油氣藏，2009，21（4）：96-98.

［10］史英，顏菲，李小波，等.考慮應(yīng)力敏感疏松砂巖氣藏試井分析［J］.巖性油氣藏，2009，21（3）：114-117.

［11］李傳亮.巖石應(yīng)力敏感指數(shù)與壓縮系數(shù)之間的關(guān)系式［J］.巖性油氣藏，2007，19（4）：95-98.

［12］李傳亮，王雙才，周涌沂.巖石壓縮系數(shù)對(duì)油藏動(dòng)態(tài)儲(chǔ)量計(jì)算結(jié)果的影響［J］.大慶石油地質(zhì)與開發(fā)，2004，23（6）：31-32.

［13］張紹槐.多分支井鉆井完井技術(shù)新進(jìn)展［J］.石油鉆采工藝，2001，23（2）：1-3.

［14］邁克爾 R.錢伯斯.多分支井［M］.孫仁遠(yuǎn)，譯.北京：石油工業(yè)出版社，2006：77-172.

［15］石曉兵，喻著成，陳平，等.側(cè)鉆水平井、分支井井眼軌跡設(shè)計(jì)與控制理論［M］.北京：石油工業(yè)出版社，2009：1-62.

［16］趙軍，王淼，祁興中，等.輪西地區(qū)奧陶系地應(yīng)力方向及裂縫展布規(guī)律分析［J］.巖性油氣藏，2010，22（3）：95-99.

Analysis of seepage-stress coupling field for lateral wells

WEI Chen-xing，LIAN Zhang-hua，DING Liang-liang，GUO Yan-ruo，HAN Tao
（State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation，Southwest Petroleum University， Chengdu 610500， China）

There is very strong fluid-solid coupling during seeping in porous media.Taking full account of fluid-solid coupling in the COMSOL software for multilateral wells,the distribution characteristics of coupling seepage and stress field are analyzed,and the dangerous section ofsidewall is recognized accordingtocalculatingthe value offail coefficient by the Coulomb criterion.Through calculating the minimum value of fail coefficient and the volume integration in which fail coefficient is less than-2.5 MPa,the impact of the build-up rate on sidewall failure is studied.The results showthat the changes ofporosityis verysmall obtained in the simple seepage field byusingunique compressibilityfactor,while the porosityand permeability of the coupling seepage field change remarkably in the vicinity of the bottom.Therefore,fluidsolid coupling seeping is a true reflection of underground seepage.In the stress field,both the maximum valve of Von Mises equivalent stress and the minimum valve of coefficients fail emerge in the side tracking point of the lower lateral,which means that the sidetrack windowis the dangerous section of sidewall.Finally,the conclusion of smallest range of rock failure drilling lower lateral with build-up rate of 15°/30m is obtained,which provides a theoretical basis for the designingand buildingofthe multilateral wells.

lateral wells；seepage field；stress field；fluid-solid coupling；numerical simulation；borehole stability

TE243

A

2011-04-02；

2011-05-26

國(guó)家自然科學(xué)基金項(xiàng)目“基于數(shù)值模擬的復(fù)雜地層地應(yīng)力場(chǎng)反演研究”（編號(hào)：50774063）資助。

魏臣興，1979年生，男，西南石油大學(xué)在讀博士研究生，主要從事CAD/CAE/CFD、管柱力學(xué)、地應(yīng)力及巖石力學(xué)方面的研究。

地址：（610500）四川省成都市新都區(qū)西南石油大學(xué)博士 2009級(jí)。 電話：（028）83032210。E-mail：weichenxing@163.com

練章華，1964年生，男，博士，教授，博士生導(dǎo)師，現(xiàn)從事CAD/CAE/CFD、管柱力學(xué)、地應(yīng)力及巖石力學(xué)等方面的教學(xué)與科研工作。E-mail：cwctlzh@swpu.edu.cn

1673-8926（2011）04-0124-05

于惠宇）