閆相禎, 董 衛(wèi), 樊 恒, 閆怡飛, 許建國,4

(1.中國石油大學(xué)油氣CAE技術(shù)研究中心,山東青島 266580; 2.中國石油大學(xué)機(jī)電工程學(xué)院,山東青島 266580;3.西安石油大學(xué)電子工程學(xué)院,陜西西安 710000; 4.中國石油吉林油田分公司,吉林松原 138000)

基于位移不連續(xù)模型的多裂縫應(yīng)力干擾極限距離計算

閆相禎1, 董 衛(wèi)2, 樊 恒3, 閆怡飛2, 許建國1,4

(1.中國石油大學(xué)油氣CAE技術(shù)研究中心,山東青島 266580; 2.中國石油大學(xué)機(jī)電工程學(xué)院,山東青島 266580;3.西安石油大學(xué)電子工程學(xué)院,陜西西安 710000; 4.中國石油吉林油田分公司,吉林松原 138000)

利用MTS試驗(yàn)裝置對“R13區(qū)塊”的致密砂巖巖心進(jìn)行小試樣試驗(yàn),獲得“R13區(qū)塊”致密砂巖的巖石力學(xué)參數(shù)(彈性模量、泊松比等)。忽略地層內(nèi)水力裂縫面的不規(guī)則性,將其簡化為互相平行的光滑裂縫面?？紤]到地層在裂縫面的不連續(xù)性,將裂縫視為連續(xù)地層的內(nèi)部不連續(xù)邊界,采用位移不連續(xù)理論建立邊界元模型,研究多裂縫系統(tǒng)周圍的地層應(yīng)力分布和不同簇數(shù)、不同布縫方式下的裂縫群縫間應(yīng)力干擾的極限距離。結(jié)果表明,人工多裂縫的出現(xiàn)明顯改變了裂縫周圍應(yīng)力場的分布,裂縫群簇數(shù)越多影響距離越大,但簇數(shù)的增加對影響距離的改變逐漸趨于平緩,呈對數(shù)增長,而段內(nèi)裂縫間距(20～30 m)的改變基本不影響應(yīng)力干擾的極限距離?；谖灰撇贿B續(xù)理論的邊界元模型能夠準(zhǔn)確地計算裂縫群周圍地層應(yīng)力場分布和裂縫群縫間應(yīng)力干擾極限距離。

多裂縫; 應(yīng)力干擾; 位移不連續(xù)理論; 極限距離

作為低滲透油氣井增產(chǎn)的主要措施,水平井分段多簇壓裂技術(shù)越來越受到重視。近幾年川慶鉆探、勘探開發(fā)研究院、吉林油田自主研發(fā)了水平井裸眼分隔器分段壓裂工具,現(xiàn)場成功試驗(yàn)百余口井,已發(fā)展成為水平井分段壓裂改造的主體技術(shù)[1];同時,對水平井分段壓裂技術(shù)機(jī)制的研究也成為了重要課題。大量的現(xiàn)場和室內(nèi)試驗(yàn)研究表明[2-10],油氣井壓裂后,人工裂縫的存在能改變近井地帶地應(yīng)力的大小和方向。隨著水平井壓裂段數(shù)和簇數(shù)的增加,縫間應(yīng)力干擾受到重視,合理地利用縫間應(yīng)力干擾可以改變裂縫起裂方向,提高儲層改造效果,但同時也會因?yàn)榭p間的應(yīng)力干擾形成局部高應(yīng)力,影響儲層改造[11-13]。為優(yōu)化水平井段的段間距和簇間距,須對地層人工裂縫群周圍應(yīng)力場的分布進(jìn)行細(xì)致研究。筆者通過對R13區(qū)塊致密砂巖巖心進(jìn)行力學(xué)試驗(yàn),獲得R13區(qū)塊巖石的力學(xué)參數(shù),采用位移不連續(xù)理論計算該區(qū)塊分段壓裂產(chǎn)生的裂縫群對地應(yīng)力分布影響的極限距離。

1 巖心力學(xué)試驗(yàn)

圖1 MTS巖石試驗(yàn)系統(tǒng)Fig.2 MTS rock-testing system

為準(zhǔn)確獲取R13區(qū)塊致密砂巖的巖石力學(xué)參數(shù),利用MTS試驗(yàn)系統(tǒng)(圖1)對R13區(qū)塊的致密砂巖巖心進(jìn)行了常規(guī)力學(xué)試驗(yàn)。根據(jù)國際巖石力學(xué)學(xué)會(ISRM)推薦的方法,將鉆取巖心制成直徑為25 mm,高為50 mm的圓柱型試件,如圖2所示。其端面不平行度不超過0.05 mm,直徑誤差不超過0.3 mm,尺寸測量精度為0.1 mm。根據(jù)油藏埋深地應(yīng)力選取圍壓為35 MPa。繪制R13區(qū)塊巖心應(yīng)力-應(yīng)變曲線如圖3所示。

圖2 砂巖巖心試件Fig.2 Specimens of sandstone cores

圖3 砂巖巖心應(yīng)力-應(yīng)變曲線Fig.3 Stress-strain curves of sandstone

圖4 徑向-軸向應(yīng)變曲線Fig.4 Radial-axial strain curves

通過巖心力學(xué)試驗(yàn),對試驗(yàn)結(jié)果進(jìn)行線性擬合,試驗(yàn)結(jié)果如圖3、 4所示。圖中藍(lán)色線表示試驗(yàn)數(shù)據(jù),紅色直線為線性擬合結(jié)果。R13區(qū)塊的巖心彈性模量為34.6 GPa,泊松比為0.14。通過力學(xué)試驗(yàn)得到R13區(qū)塊的巖石力學(xué)參數(shù),為后續(xù)裂縫應(yīng)力干擾極限的研究提供了基礎(chǔ)數(shù)據(jù)。

2 位移不連續(xù)理論基本解

位移不連續(xù)法是由美國學(xué)者Crouch于1973年在國際巖石力學(xué)會議上第一次提出,并給出了位移不連續(xù)理論的二維解[14-16]。迄今為止,已有大量的學(xué)者運(yùn)用位移不連續(xù)理論進(jìn)行巖石力學(xué)問題的研究[17-18]。Kresse和Wu等[19-20]采用改進(jìn)的二維位移不連續(xù)法研究裂縫擴(kuò)展問題。

當(dāng)進(jìn)行水力壓裂時,地層中會出現(xiàn)水力裂縫,裂縫的上下表面間會發(fā)生相對位移和錯動,即為所說的位移不連續(xù)量。故而將裂縫的兩個面視為地層內(nèi)部的不連續(xù)邊界。建立長為2l的不連續(xù)單元如圖5所示。裂縫面x與y方向的不連續(xù)量分別用Dx、Dy表示,相對的單個面的位置則用u(x,0-)、u(x,0+)、u(y,0-) 、u(y,0+)表示。裂縫面的位移不連續(xù)量可以表示為

Dx=ux(x,0-)-ux(x,0+),

Dy=uy(y,0-)-uy(y,0+),|x|≤l.

(1)

式中,正號表示上表面,負(fù)號表示下表面;裂縫張開Dy為負(fù)。

圖5 二維位移不連續(xù)單元Fig.5 Two-dimensional displacement discontinuity element

Crouch和Starfied給出了計算位移和地應(yīng)力的計算公式[15]如下:

(2)

式中,ux為水平位移;uy為垂直位移;Dx為切向不連續(xù)位移;Dy為法向不連續(xù)位移。

f是與點(diǎn)的位置坐標(biāo)有關(guān)的函數(shù),表述為

(3)

(4)

(5)

(6)

(7)

(8)

(9)

則x與y方向的應(yīng)力及剪應(yīng)力表達(dá)式為

(10)

3 R13區(qū)塊裂縫群間應(yīng)力影響范圍分析

依據(jù)位移不連續(xù)原理,采用R13區(qū)塊的巖心試驗(yàn)結(jié)果,對R13區(qū)塊的裂縫間應(yīng)力干擾進(jìn)行了研究。

3.1 模型建立

水平井每段簇數(shù)分別為2、3、4、5、6簇。圖6為多裂縫壓裂過程的產(chǎn)層受力示意圖。模擬計算中,產(chǎn)層模型尺寸取800 m×400 m,裂縫長度為100 m,井眼位置位于裂縫中間;AB、CD加載最大水平主應(yīng)力σH;AB、BC邊加載最小主應(yīng)力σh;初始地層壓力為p0。

圖6 地層受力示意圖Fig.6 Diagram of strata stress

有限元模型尺寸為600 m×1 000 m,裂縫長度為100 m,最大、最小水平地應(yīng)力分別為24 和20 MPa,砂巖彈性模量為34.6 GPa,泊松比為 0.14,地層孔隙壓力為19 MPa,地層孔隙度為9%。結(jié)合水平井壓裂射孔工具,簇間裂縫間距取25 m。

3.2 計算結(jié)果分析

為更直觀地看出地層應(yīng)力改變區(qū)域,將計算結(jié)果制成云圖如圖7、8所示。結(jié)合工程實(shí)際,當(dāng)應(yīng)力改變值小于5%時,認(rèn)為地層應(yīng)力值未發(fā)生變化,即為圖中藍(lán)色區(qū)域。

圖7 不同簇數(shù)下地層最大主應(yīng)力云圖Fig.7 Maximum main stress nephogram under different clusters

圖8 不同簇數(shù)下地層最小主應(yīng)力云圖Fig.8 Minimum main stress nephogram under different clusters

從圖7和圖8可以得出,最大主應(yīng)力和最小主應(yīng)力的影響范圍均隨著壓裂簇數(shù)的增加而增加。將圖中數(shù)據(jù)提取作成曲線如圖9所示。由圖9可以看出,在裂縫簇數(shù)小于等于4時應(yīng)力影響范圍隨著簇數(shù)增加,增加較快,當(dāng)裂縫簇數(shù)大于4時增加較慢;垂直于裂縫方向的影響距離要大于平行于裂縫延伸方向的。改變簇間距(30 m)進(jìn)行了相同的處理方法,得到的曲線與圖9和圖10基本一致。

圖9 不同壓裂簇數(shù)下最大和最小主應(yīng)力影響范圍Fig.9 Reach of maximum and minimum principal stresses with different fracturing cluster number

4 不同布縫方式下的應(yīng)力干擾極限距離

針對油田現(xiàn)場的情況,采取了兩種布縫方式(第一種是垂直于裂縫延伸方向,第二種是平行于裂縫延伸方向)見圖10。

圖10 第一種和第二種布縫方式示意圖Fig.10 Diagrams under the first and the second decorate fractures

根據(jù)圖10的布縫方式,分別建立計算模型,將計算結(jié)果列于表1、2,表中數(shù)據(jù)為縫長的倍數(shù)。從表中可以看出,每個裂縫群簇數(shù)越多應(yīng)力影響極限距離越長,但增加逐漸放緩;最大主應(yīng)力要比最小主應(yīng)力的影響范圍大;水平井間距應(yīng)小于表中給出的極限距離,以利于形成縫網(wǎng)。

表1 第一種布縫方式下不同簇數(shù)下縫間應(yīng)力干擾極限距離Table 1 Seam stress interference limit distance of different clusters under the first decorate fractures

表2 第二種布縫方式下不同簇數(shù)下縫間應(yīng)力干擾極限距離Table 2 Seam stress interference limit distance of different clusters under the second decorate fractures

對比這兩種布縫方式的計算結(jié)果可以看出,第一種布縫方式應(yīng)力的影響極限距離明顯大于第二種布縫方式下的極限距離,因此在實(shí)際生產(chǎn)中建議采用第一種布縫方式,以利于形成縫網(wǎng)。

5 結(jié) 論

(1)由于人工裂縫的存在,明顯改變了裂縫周圍地層應(yīng)力場分布。采用位移不連續(xù)法得到的計算模型可以較好地計算裂縫地層應(yīng)力分布。

(2)每個裂縫群內(nèi)簇數(shù)越多應(yīng)力影響距離越長,但增加逐漸放緩,當(dāng)簇數(shù)大于4時,應(yīng)力影響的極限距離增加速率接近于0。

(3)最大主應(yīng)力比最小主應(yīng)力的影響范圍大;

簇間距對應(yīng)力干擾距離基本無影響。

(4)第一種布縫方式的影響極限距離大于第二種布縫方式的,實(shí)際生產(chǎn)中建議采用第一種布縫方式,以利于縫網(wǎng)形成。

[1] 王峰,劉長宇,張應(yīng)安,等.水平井叢書——吉林油田第四分冊[M]. 北京:石油工業(yè)出版社,2012.

[2] 岳迎春,郭建春.重復(fù)壓裂轉(zhuǎn)向機(jī)制流-固耦合分析[J].巖土力學(xué),2012,33(10):3189-3193.

YUE Yingchun, GUO Jianchun. Fluid-solid coupling analysis of reorientation mechanism of refracturing[J]. Rock and Soil Mechanics,2012,33(10):3189-3193.

[3] 張廣清,陳勉,趙艷波.新井定向射孔轉(zhuǎn)向壓裂裂縫起裂與延伸機(jī)制研究[J].石油學(xué)報,2008,29(1):116-119.

ZHANG Guangqing, CHEN Mian, ZHAO Yanbo. Study on initiation and propagation mechanism of fractures in oriented perforation of new wells[J]. Acta Petrolei Sinica, 2008,29(1):116-119.

[4] 陳勉.頁巖氣儲層水力裂縫轉(zhuǎn)向擴(kuò)展機(jī)制[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2013,37(5):88-94.

CHEN Mian. Re-orientation and propagation of hydraulic fractures in shale gas reservoir[J]. Journal of China University of Petroleum(Edition of Natural Science), 2013,37(5):88-94.

[5] 周健,陳勉,金衍,等.壓裂中天然裂縫剪切破壞機(jī)制研究[J].巖石力學(xué)與工程學(xué)報,2008,27(S1):2637-2641.

ZHOU Jian, CHEN Mian, JIN Yan, et al. Mechanism study of shearing slippage damage of natural fracture in hydraulic fracturing[J]. Chinese Journal of Rock Mechanics and Engineering, 2008,27(S1):2637-2641.

[6] 陳勉,周健,金衍,等.隨機(jī)裂縫性儲層壓裂特征實(shí)驗(yàn)研究[J].石油學(xué)報,2008,29(3):431-434.

CHEN Mian, ZHOU Jian, JIN Yan, et al. Experimental study on fracturing features in naturally fractured reservoir[J]. Acta Petrolei Sinica, 2008,29(3):431-434.

[7] ROUSSEL N P, SHARMA M M. Quantifying transient effects in altered-stress refracturing of vertical wells[J]. SPE Journal, 2010,15(3):770-782.

[8] SIEBRITS E, ELBEL J L, HOOVER R S, et al. Refracture reorientation enhances gas production in Barnett shale tight gas wells[R]. SPE 63030, 2000.

[9] 曲占慶,田雨,李建雄,等.水平井多段分簇壓裂裂縫擴(kuò)展形態(tài)數(shù)值模擬[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2017,41(1):102-105.

QU Zhanqing, TIAN Yu, LI Jianxiong, et al. Numerical simulation study on fracture extension and morphology of multi-cluster staged fracturing for horizontal wells[J] Journal of China University of Petroleum(Edition of Natural Science), 2017,41(1):102-105.

[10] 閆相禎,李向陽.基于裂縫干涉模型的非常規(guī)油氣井壓裂優(yōu)化設(shè)計軟件的開發(fā)與應(yīng)用[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2013,37(5):120-128.

YAN Xiangzhen, LI Xiangyang. Development and application of unconventional oil & gas well fracturing optimization design software based on fracture network interference model[J]. Journal of China University of Petroleum(Edition of Natural Science), 2013,37(5):120-128.

[11] 趙金洲,許文俊,李勇明,等.低滲透油氣藏水平井分段多簇壓裂簇間距優(yōu)化新方法[J].天然氣工業(yè),2016,36(10):63-69.

ZHAO Jinzhou, XU Wenjun, LI Yongming, et al. A new method for cluster spacing optimization of multi-cluster staged fracturing in horizontal wells of low-permeability oil and gas reservoirs[J]. Natural Gas Industry, 2016,36(10):63-69.

[12] CIPOLLA C L, WARPINSKI N R, MAYERHOFER M J. Hydraulic fracture complexity: diagnosis, remediation, and explotation[R]. SPE 115771,2008.

[13] 趙金洲,陳曦宇,劉長宇,等.水平井分段多簇壓裂縫間干擾影響分析[J].天然氣地球科學(xué),2015,26(3):533-538.

ZHAO Jinzhou, CHEN Xiyu, LIU Changyu, et al. The analysis of crack interaction in multi-stage horizontal fracturing[J]. Natural Gas Geoscience, 2015,26(3):533-538.

[14] CROUCH S L. Solution of plane elasticity problems by the displacement discontinuity method: I infinite body solution[J]. International Journal for Numerical Methods in Engineering, 1976,10(2):301-343.

[15] CROUCH S L, STARFIELD A M, RIZZO F J. Boundary element methods in solid mechanics [M]. George: Allen & Unwin, 1983.

[16] 王飛,董越鵬,黃醒春.彈性層狀半平面的位移不連續(xù)基本解[J].上海交通大學(xué)學(xué)報,2009,43(10):1675-1680.

WANG Fei, DONG Yuepeng, HUANG Xingchun. Displacement discontinuity fundamental solutions for a layered elastics half-plane[J].Journal of Shanghai Jiaotong University, 2009,43(10):1675-1680.

[17] OLSON J E, TALEGHANI A D. Modeling simultaneous growth of multiple hydraulic fractures and their interaction with natural fractures[R]. SPE 163982, 2009.

[18] CHENG Y. Mechanical interaction of multiple fractures--exploring impacts of the selection of the spacing/number of perforation clusters on horizontal shale-gas wells[R]. SPE 125769, 2012.

[19] WU R, KRESS O, WENG X. Modeling of interaction of hydraulic fractures in complex fracture networks[R]. SPE 152052, 2012.

[20] KRESSE O, WENG X, GU H, et al. Numerical modeling of hydraulic fractures interaction in complex naturally fractured formations[J]. Rock Mechanics and Rock Engineering, 2013,46(3):555-568.

Criticaldistanceofstressshadowofmultiple-cracksbasedondisplacementdiscontinuitymodel

YAN Xiangzhen1, DONG Wei2, FAN Heng3, YAN Yifei2, XU Jianguo1,4

(1.CAETechnologyResearchCenter,ChinaUniversityofPetroleum,Qingdao266580,China;2.CollegeofMechanicalandElectronicEngineeringinChinaUniversityofPetroleum,Qingdao266580,China;3.CollegeofElectronicEngineering,XianShiyouUniversity,Xian710000,China;4.PetroChinaJilinOilfieldCompany,Songyuan138000,China)

A small specimen test is conducted in MTS to get the elastic modulus and Poissons ratio of tight sandstone core in R13 field. The irregularity of hydraulic crack surface is ignored, and the cracks are simplified as smooth and parallel to each other. Considering the discontinuity of the multi-cracks, they are treated as discontinuous boundaries of the continuous stratum. Then the displacement discontinuity method is used to construct a boundary element model, which is applicable to calculate the distribution of stress around the multi-crack system, and determine the critical distance of the crack swarmsstress shadow with different clusters and different ways of cloth sewn. The results show that the appearance of artificial multiple-cracks has significantly changed the distribution of stress field around the cracks. The more the cluster number, the longer the distance takes effect. But the influence of the cluster number increasing on the distance becomes smaller, taking on a logarithmic tendency. However, the change of crack spacing between 20 and 30 meters, in a single stage has little influence on the critical distance of stress shadow. This study indicates that the boundary element model based on the displacement discontinuity theory can accurately predict the distribution of the stress field around the multiple-cracks and the critical distance of the crack swarmsstress shadow.

multiple-cracks; stress shadow; displacement discontinuity method; critical distance

2017-03-28

國家自然科學(xué)基金項(xiàng)目(51374228); 中國博士后科學(xué)基金項(xiàng)目(2017M612375);國家科技重大專項(xiàng)(2016ZX05017-003-01);中國石油天然氣集團(tuán)公司重點(diǎn)實(shí)驗(yàn)室項(xiàng)目(2016A-3905)

閆相禎(1956-),男,教授,博士生導(dǎo)師,研究方向?yàn)橛蜌夤こ塘W(xué)。E-mail:yanxzh@163.com。

樊恒(1987-),男,講師,博士,研究方向?yàn)橛蜌夤こ塘W(xué)與安全。E-mail:fanheng1@126.com。

1673-5005(2017)06-0140-07

10.3969/j.issn.1673-5005.2017.06.017

TE 357

A

閆相禎,董衛(wèi),樊恒,等. 基于位移不連續(xù)模型的多裂縫應(yīng)力干擾極限距離計算[J].中國石油大學(xué)學(xué)報(自然科學(xué)版), 2017,41(6):140-146.

YAN Xiangzhen, DONG Wei, FAN Heng, et al. Critical distance of stress shadow of multiple-cracks based on displacement discontinuity model[J]. Journal of China University of Petroleum (Edition of Natural Science), 2017,41(6):140-146.

(編輯 沈玉英)