夏　沖， 廖　旋， 劉荊成

(長江大學 石油工程學院，湖北 武漢 430100)

?

多分支縫水力壓裂暫堵起裂點位置研究

夏沖， 廖旋， 劉荊成

(長江大學 石油工程學院，湖北 武漢 430100)

針對多分支縫壓裂成縫起裂點不明確的問題，基于多孔介質(zhì)流固耦合理論模型，通過Comsol Multiphysics軟件進行有限元模擬分支縫壓裂過程，分析不同應力差條件下縫內(nèi)起裂點和井筒起裂點的集中應力變化規(guī)律，隨著縫內(nèi)凈壓力的增加，存在一個差別點，區(qū)分起裂點起裂的優(yōu)先順序。通過控制縫內(nèi)凈壓力變化，調(diào)整起裂優(yōu)先順序，控制分支縫延伸條數(shù)，保證了多分支縫水力壓裂縫網(wǎng)的復雜性。對海拉爾油田現(xiàn)場實例計算分析，調(diào)整縫內(nèi)凈壓力，施工結果與模擬結果相同，可以指導現(xiàn)場施工，提高多分支縫水力壓裂的成功率。

多分支縫；水力壓裂；起裂點；有限元；優(yōu)先起裂

低產(chǎn)井增產(chǎn)改造工藝主要是沿用了常規(guī)壓裂工藝[1-2](如重復壓裂、轉(zhuǎn)向壓裂、多裂縫壓裂等)，這類常規(guī)壓裂工藝技術針對性不強，人工裂縫波及儲層體積小，瀉油面積有限，措施效果差。多分支縫暫堵壓裂提高縫內(nèi)靜壓力[3]，暫堵條件下，多條裂縫分別從井筒及縫內(nèi)延伸，很好地解決了人工裂縫波及體積小的問題，能夠極大地提高產(chǎn)量。

多分支縫水力壓裂技術在部分油田成功實施，增產(chǎn)效果明顯。海拉爾油田興安嶺儲層多分支縫壓裂試驗平均單井增油3.9 t/d，與常規(guī)壓裂相比增產(chǎn)效果提高34.4%。國內(nèi)趙金洲、程遠方等[4-5]研究了縫內(nèi)存在天然裂縫情況下分支縫的起裂及延伸機理，然而多分支縫水力壓裂暫堵后裂縫形成機理并不明確，暫堵條件下難以確定分支縫的起裂位置。本文以多孔介質(zhì)流固耦合理論作為依據(jù)，建立多分支壓裂有限元模型，分析井筒起裂及縫內(nèi)起裂的壓力變化，確定分支縫起裂點位置，最后進行現(xiàn)場數(shù)據(jù)算例分析。

1　多孔介質(zhì)流固耦合理論模型

多孔介質(zhì)流固耦合理論模型包括基于單向流滲流方程的滲流場、線彈性本構方程的變形場、滲流與變形耦合場三部分組成[6]。

滲流場方程為：

(1)

變形場平衡方程為：

(2)

劉建軍等[7]通過實驗的方法得出低滲透儲層滲透率可以表示為：

(3)

式中：K為儲層滲透率，μm2；μ為油相黏度，Pa·s；Ct為綜合壓縮系數(shù)，1/Pa；p為油相壓力，Pa；εV為體積應變；t為時間，s；σij為有效應力分量，Pa；fi為體積分量；α為Biot系數(shù)；a為實驗擬合敏感性系數(shù)，10-3μm2；b為實驗擬合敏感性系數(shù)；σ為有效應力，MPa。

上述方程聯(lián)立構成流固耦合的基本控制方程，在相應邊界條件以及初始條件下，構成多孔介質(zhì)流固耦合理論模型?；谏鲜龆嗫捉橘|(zhì)流固耦合模型，利用Comsol Multiphysics軟件，建立多分支縫水力壓裂有限元模型。

2　近井區(qū)域單縫壓裂有限元模擬

選擇多孔彈性物理場(即流固耦合場),以致密砂巖儲層為例。在模型中心建立井筒，取地層彈性模量E為4 GPa，泊松比為0.26，最大、最小水平主應力分別為30、26 MPa，液柱壓力為30 MPa時，對有限元模型進行計算分析，結果如圖1、2所示。

圖1　初始裂縫起裂點

從圖1中可以看出，單縫壓裂時，裂縫沿最大主應力方向延伸。如圖2所示，縫內(nèi)壓力在井筒附近較大，隨著裂縫的延伸，縫內(nèi)應力逐漸變小，裂縫尖端存在應力集中，符合水力壓裂經(jīng)典理論及實驗[8-11]。因此本模型可以進行水力壓裂過程的分析。

圖2　單縫條件縫內(nèi)壓力分布

3　多分支縫起裂點應力分析

多分支縫水力壓裂是在封堵先導裂縫條件下，對裂縫再次壓裂，結果如圖3所示。從圖3中可以看出，分支縫二次壓裂條件下，裂縫區(qū)域存在井筒、縫內(nèi)、封端3個應力集中區(qū)域，然而由于暫堵的存在封端壓力不能傳播，因而存在縫內(nèi)、封端兩個起裂點的可能。

圖3　分支縫二次壓裂效果

在最大主應力為30 MPa，不同應力差條件下，對模型進行分支縫二次壓裂計算分析，結果如圖4、5所示。

圖4　縫內(nèi)集中應力隨凈壓力變化值

Fig.4Concentrated stress in the seam with the net pressure change values

圖4為不同應力差條件下，縫內(nèi)最大集中應力隨凈壓力值變化曲線。從圖4中可以看出，凈壓力越大，縫內(nèi)集中應力越明顯，越接近起裂。應力差越大，縫內(nèi)集中應力值越大?？梢猿醪脚袛鄳Σ钤酱蟮牡貙舆M行多分支壓裂時，初次初次裂縫縫內(nèi)更易起裂。

圖5　井筒集中應力隨凈壓力變化值

Fig.5Concentrated stress in wellbore with the net pressure change values

圖5為不同應力差條件下，井筒集中應力隨凈壓力變化曲線。從圖5中可以看出，隨著縫內(nèi)凈壓力的增大，井筒集中應力先減小后增大，當縫內(nèi)凈壓力大于0時，井筒集中應力和凈壓力變化趨勢呈正相關。與縫內(nèi)集中應力不同的是，應力差越大和井筒集中應力變化趨勢呈現(xiàn)負相關。容易得出應力差小的地層井筒應力集中大，分支縫容易從井筒起裂的結論。

4　起裂點起裂順序分析

由于多分支的起裂點存在兩個，因而不同起裂點的優(yōu)先起裂決定了不同的裂縫延伸情況。當縫內(nèi)優(yōu)先起裂時，雖然存在起裂時的卸壓，由于井筒的集中應力過大，井筒仍然存在起裂可能。如果井筒優(yōu)先起裂，導致初次壓裂裂縫壓力下降，縫內(nèi)起裂點便不能起裂。因此，根據(jù)地層特征，通過調(diào)整縫內(nèi)靜壓力，控制二次壓裂起裂點的優(yōu)先順序，便可以達到增加縫網(wǎng)復雜結構，擴大波及體積的目的。

縫內(nèi)起裂及井筒起裂的優(yōu)先順序可以簡單的通過起裂點的集中應力差來體現(xiàn)。相同條件下，集中應力越大的點優(yōu)先起裂。圖6為等效應力差隨凈壓力變化情況，等效應力差大于0時，縫內(nèi)優(yōu)先起裂，反之井筒優(yōu)先起裂。從圖6中可以看出，存在一個差別值，當縫內(nèi)壓裂小于此值時，縫內(nèi)優(yōu)先起裂，大于此值時井筒優(yōu)先起裂。等效應力差越大，差別值越大，因而等效應力差越大的地層，縫內(nèi)起裂的可能性越大。為了保證分支縫的復雜性，在保證達到縫內(nèi)集中應力達到破裂值的條件下，保持縫內(nèi)靜壓力不要高于差別值，令縫內(nèi)優(yōu)先起裂，進而增大流量，使縫內(nèi)井筒同時延伸，達到如圖7所示的分支縫效果。

圖6　等效應力差隨凈壓力變化

Fig.6The equivalent stress difference with the net pressure change

圖7　多分支縫裂縫形態(tài)

Fig.7Fracture morphology of multi-branch seam

5　現(xiàn)場實例

根據(jù)本文的有限元分析方法，按照興安嶺儲層特點及產(chǎn)生縫網(wǎng)壓裂的界限條件[12]，對貝28-X62-58井采用大規(guī)模多分支縫壓裂工藝。該層段的巖石力學參數(shù)和地應力數(shù)據(jù)見表1，壓裂施工曲線見圖8。

表1　貝28壓裂層段的巖石物理力學解釋數(shù)據(jù)

圖8　XⅡ16～XⅡ11-12層壓裂施工曲線圖

Fig.8The graph of fracturing in XⅡ16～ XⅡ11-12layers

根據(jù)壓裂井數(shù)據(jù)，應用有限元軟件對本井初始縫進行有限元數(shù)值模擬。有限元模擬結果如圖9所示，應力差隨凈壓力變化結果如圖10所示。

圖9　分支縫壓裂集中應力點

Fig.9Concentrated stress points of branch seam

從圖9和圖10中可以看出，縫內(nèi)凈壓力達到28 MPa時，才會從井筒優(yōu)先起裂，因此該次多分支縫壓裂縫內(nèi)凈壓力保持為2～3 MPa，預測從縫內(nèi)優(yōu)先起裂，然后縫內(nèi)和井筒同時延伸與現(xiàn)場檢測的實際情況相符。

圖10　應力差隨凈壓力變化

Fig.10The stress difference with the net pressure change

6　結論

(1) 多分支壓裂二次壓裂的集中應力點存在縫內(nèi)和井筒，起裂點的應力集中隨著縫內(nèi)凈壓力的增大而增大，兩個起裂點存在同時起裂的可能性?？p內(nèi)起裂的集中應力值和應力差呈現(xiàn)正相關；井筒起裂的集中應力值和應力差呈現(xiàn)負相關。

(2) 若縫內(nèi)優(yōu)先起裂，不影響井筒隨后起裂，兩條裂縫同時延伸；若井筒優(yōu)先起裂，則不會出現(xiàn)縫內(nèi)分支起裂情況，通過調(diào)整縫內(nèi)凈壓力可以控制分支縫起裂的優(yōu)先順序，保證分支縫網(wǎng)的復雜性。

(3) 根據(jù)該模型判斷海拉爾貝28井分支縫壓裂起裂點的優(yōu)先起裂順序，進一步預測其縫網(wǎng)形態(tài)，現(xiàn)場結果相同，驗證了模型的準確性，并且指導現(xiàn)場工作。

[1]董光，鄧金根，朱海燕，等.重復壓裂前的地應力場分析[J].斷塊油氣田，2012，19(4)：485-488；492.

Dong Guang，Deng Jingen，Zhu Haiyan，et al. Analysis of stress field before refracturetreatment[J]. Fault-Block Oil & Gas Field，2012，19(4)：485-488；492.

[2]羅天雨, 趙金洲, 王嘉淮,等. 復雜裂縫產(chǎn)生機理研究[J]. 斷塊油氣田, 2008, 15(3):46-48.

Luo Tianyu, Zhao Jinzhou,Wang Jiahuai, et al. Forming mechanism of complicated fractures in hydraulic fracturing[J]. Fault-Block Oil & Gas Field, 2008,15(3)：46-48.

[3]甘云雁，張士誠，劉書杰，等．整體壓裂井網(wǎng)與裂縫優(yōu)化設計新方法[J].石油學報，2011，32(2): 290-294.

Gan Yunyan，Zhang Shicheng，Liu Shujie，et al. A new method for well pattern optimization and integral fracturing design in low permeability reservoirs[J].Acta Petrolei Sinica，2011，32(2): 290-294.

[4]趙金洲，楊海，李勇明，等．水力裂縫逼近時天然裂縫穩(wěn)定性分析[J]．天然氣地球科學，2014,25(3)：402-408.

Zhao Jinzhou, Yang Hai, Li Yongming, et al. Stability of the natural fracture when the hydraulic fracture is approaching[J].Natural Gas Geoscience, 2014,25(3)：402-408.

[5]程遠方，常鑫，孫元偉，等．基于斷裂力學的頁巖儲層縫網(wǎng)延伸形態(tài)研究[J]．天然氣地球科學，2014,25(4):603-611.

Cheng Yuanfang, Chang Xin, Sun Yuanwei, et al. Research on fracture network propagation pattern of shale reservoir based on fracture mechanics[J].Natural Gas Geoscience, 2014,25(4):603-611.

[6]李小波，史英，張修明，等.水力壓裂裂縫數(shù)值模擬研究[J].西安石油大學學報(自然科學版)，2009,24(3)：52-55.

Li Xiaobo, Shi Ying, Zhang Xiuming, et al. Study on the simulation of hydraulic fracturing cracks[J].Journal of Xi’an Shiyou University(Natural Science Engineering),2009,24(3)：52-55.

[7]劉建軍，馮夏庭. 我國油藏滲流-溫度-應力耦合的研究進展[J].巖土力學，2003, 24(增刊):645-650.

Liu Jianjun, Feng Xiating. Advance of studies on thermo-hydro-mechanical interaction in oil reservoir in China[J]. Rock and Soil Mechanics, 2003, 24(Supplement):645-650.

[8]金衍, 陳勉, 周健, 等. 巖性突變體對水力裂縫延伸影響的實驗研究[J]. 石油學報, 2008, 29(2): 300-303.

Jin Yan, Chen Mian, Zhou Jian, et al. Experimental study on the effects of salutatory barrier on hydraulic fracture propagation of cement blocks[J]. Acta Petrolei Sinica, 2008,29(2): 300-303.

[9]程遠方，徐太雙，吳百烈，等．煤巖水力壓裂裂縫形態(tài)實驗研究[J]．天然氣地球科學，2013,24(1):134-137.

Cheng Yuanfang, Xu Taishuang, Wu Bailie,et al.Experimental study on the hydraulic fractures morphology of coal bed [J].Natural Gas Geoscience,2013,24(1):134-137.

[10]王鳳江，丁云宏，路勇．低滲透油田重復壓裂技術研究[J]．石油勘探與開發(fā)，1999，26(1): 71-73.

Wang Fengjiang, Ding Yunhong, Lu Yong. Study on refracturing technology in low permeability oil reservoir [J]．Petroleum Exploration and Development，1999，26(1): 71-73.

[11]Chen Yiheng, Hasebe N. A consistency check for strongly interacting multiple crack problems in isotropic, bimaterial and orthotropic bodies[J]. International Journal of Fracture, 1998, 89(4): 333-353.

(編輯宋官龍)

The Initial Cracking Point of Multi-Branch Hydraulic Fracture Temporary Plugging

Xia Chong， Liao Xuan， Liu Jingcheng

(CollegeofPetroleumEngineering,YangtzeUniversity,WuhanHubei430100,China)

The seam crack point of multiple branch seam fracturing is not clear. Based on the theory of porous media of fluid-solid coupling mode, Comsol software was used to simulate branch seam fracturing process using finite element simulation. The changing law of stress concentration of seam crack point and shaft crack point was analysized under the condition of different stress difference it was found that with the increase of net pressure on the seam, there was a difference which distinguished between crack priorities. By controlling the net pressure change within joint, the crack priorities was adjusted, the number of multi-branch fracture was controlled, and the multi-branch seam hydraulic pressure of the complexity of the fracture network was ensured. Calculation analysis of hailaer oil field on the spot was used to adjust the net pressure inside seam. The result of construction was same with simulation, and can be used to guide the construction to improve the success rate of multi-branch seam hydraulic fracturing.

Multi-branch; Hydraulic fracture; Crack point; Finite element; Crack priorities

1006-396X(2016)02-0044-05

投稿網(wǎng)址：http://journal.lnpu.edu.cn

2015-11-17

2016-01-14

國家自然科學基金項目(90210022)。

夏沖(1993-)，男，碩士研究生，從事石油與天然氣工程研究；E-mail：623064276@qq.com。

TE357

Adoi:10.3969/j.issn.1006-396X.2016.02.009