摘 要:為揭示水平井分段改造中的應力陰影效應及其影響因素,基于水平井分段體積改造理論,借鑒Bunger的應力場分析假設條件,從水平井單一橫切縫入手,建立起了考慮應力陰影的應力場疊加數學模型并推導出數學公式,在此基礎上得出簡化的水平井多段應力陰影效應應力場模型。對應力陰影效應影響因素的分析結果表明:儲層地質參數中楊氏模量、水平主應力差對應力陰影效應的影響大,而泊松比對其影響則很?。皇┕抵惺┕づ帕繉﹃幱靶挠绊懨黠@,在一定儲層條件下,分段及簇間距越小,應力陰影效果對水平井裂縫附近其他裂縫的擴展產生應力陰影影響越大。利用該研究成果對現場實際水平井壓裂進行了優化設計,微地震監測結果表明:利用應力陰影效應只壓裂6段就實現了以往l0段壓裂的改造體積,較之以往減少壓裂4段,節約成本在30%以上。這對合理利用應力陰影效應進行水平井體積改造裂縫尺寸設置及工藝參數優化具有重要的指導意義。
關鍵詞:應力陰影 水平井 分段改造體積 橫切縫 楊氏模量 凈壓力 水平主應力差 泊松比 微地震 優化設計
Influence of stress shadow on horizontal well fracturing
Abstract:In order to reveal stress shadow effect and its influencing factors in the processing of segregated fracturing in horizontal wells,starting with an analysis of single transverse fracture of horizontal wells,we established a stress field superposition model and developed a mathematical formula considering stress shadow according to the theories of segregated SRV(stimulated reservoir volume)fracturing in horizontal wells and stress field analysis assumed conditions proposed by Bunger.On this basis,a simplified stress field model was established for the segregated stress shadow of a horizontal well.The influencing factor analysis shows that among reservoir geological parameters,Young's modulus and horizontal principal stress difference greatly impact stress shadow,while the Poisson's ratio slightly affects stress shadow;among construction parameters,the displacement obviously affects stress shadow;in certain reservoir conditions,the smaller the distance between the segments clusters,the stronger stress shadow on the extension of other fractures around horizontal well fractures.The research results were used to optimize the design of actual horizontal wells.Micro seismic monitoring indicated that stress shadow reduced three to four segments to achieve the same stimulated reservoir volulne,with the cost saved by more than 30%.Therefore,it is of an important guiding significance for engineers to rationally use stress shadow to optimize fracturing design and process parameters during horizontal well fracturing.
Keywords:stress shadow,horizontal well,staged-stimulated reservoir volume,transverse fracture,Young's modulus,net pressure,horizontal principal stress difference,Poisson’s ratio,micro seismic,optimal design
近年來,國內外非常規領域水平井及體積壓裂技術的進步與規模應用,使非常規油氣資源得以高效經濟開發并發揮出革命性作用[1-2]。但在水平井壓裂設計中,裂縫間距及分簇間距這一影響產量、采收率和經濟效益的重要因素仍尚不清楚,其核心在于水力裂縫產生的應力陰影效應(stress shadow)的表征與計算仍不明朗。自20世紀80年代,Warpinski和Teufel已開始研究平行裂縫之間的相互作用,并給出開啟一條恒定裂縫高度、有限裂縫長度的裂縫后,裂縫周圍產生應力場的解析解[3];2004年Fisher和2006年Mayethofer在微地震的礦場研究中證實這一現象;2009年Olson得出這一效應的一個明顯作用效果是在多條平行裂縫中,中間的裂縫的裂縫寬度最小,這是因為兩邊的裂縫對中間的裂縫都產生壓縮應力[4]。2010年Bunger等在研究Barnett頁巖氣時進一步證實這種效應的存在,并將這種效應稱為應力陰影效應[5];2011年Roussel和Sharma指出水平井改造中,一定裂縫間距內前一條支撐裂縫會導致下一條裂縫附近的應力轉向,在下一相鄰段裂縫中出現與主縫向切的裂縫系統[6]。國內學者趙金洲等也在研究我國頁巖氣壓裂設計中,指出壓裂在單一裂縫脆弱面上將產生誘導應力,并指出產生應力陰影的誘導應力可以改變最大與最小主應力的分布,使裂縫發生轉向[7]。前人多數研究集中在靜態條件下的兩條平行裂縫間的相互作用上,對多段多條裂縫下如何表征應力陰影及其影響因素仍是有待解決的問題,研究也進一步證實應力陰影效應及其影響因素是合理設計并利用人工裂縫最大限度發揮裂縫作用和效果的重要科學問題。因此,筆者在上述前人研究的基礎上進行了深入研究。
1 應力陰影效應
隨著工具設備及水平井改造技術的不斷進步,儲層多段改造的規模、段數等越來越多[8-9]。水力壓裂沒計師可利用分簇射孔將分段改造距離設置得越來越近,但往往過多的分簇與分段未必就帶來理想的效果,Roussel認為應力陰影效應的產生機理主要在于支撐劑的填充導致人工裂縫附近引起了應力場的變化,蔣廷學等也在研究我國頁巖氣壓裂設計中指出,壓裂在單一裂縫脆弱面上產生誘導應力的存在,并指出產生應力陰影的誘導應力可以改變最大與最小主應力的分布,使裂縫發生轉向[10]。述研究表明,應力陰影效應可引起人工裂縫從水平井的橫切縫中開始偏離朝向或遠離以前的裂縫方向延伸,使得下一裂縫發生轉向。Barnctt頁巖裂縫間距30m時水平井壓裂裂縫產生的應力陰影效應如圖1-a所示,可見從第二條裂縫開始裂縫出現了與主縫向切的裂縫系統,而Eagle Ford頁巖壓裂的裂縫監測情況,同樣是30m的裂縫間距沒有出現裂縫應力陰影效應(圖1-b)。因此研究陰影效應的產生機理,建立分析數學模型并進行其影響因素分析對水平井分段體積改造優化設計及效果的提升具有重要研究與應用價值。
2 數學模型
2.1 假設條件
假定流體流動滿足應力 滲流耦合連續性方程、運動方程和狀態方程[11]。裂縫的擴展與延伸模型復合經典的Perkins、Kern和Norgren(PKN)與Khristianovic、Geertsma和Dekerk(KGD)理論模型,具體假設條件為:①油氣藏為均質無限大且各向同性;②流動流體為牛頓流體;③不考慮地層溫度對裂縫內流體的影響;④暫不考慮裂縫系統內啟動壓力的影響;⑤巖石破裂滿足摩爾庫倫屈服準則;⑥巖石的變形為彈性變形,不考慮巖石蠕變。
2.2 單縫應力場模型
為能合理有效研究多條裂縫間應力陰影效應引起的誘導應力變化情況,首先建立以單一裂縫為主的均質,各向同性,彈性和有界的三維裂縫模型(圖2),在此基礎上,利用應力疊加原理進行多條裂縫模型研究。假設巖石上下隔層的楊氏模量和泊松比為Eb、ub,生產層油層的楊氏模量和泊松比為Eb、up。,裂縫高度并沒有失控(即hf<hp),裂縫的閉合壓力可以利用凈壓力和最小水平應力計算得到。即
因此對于單一裂縫,可通過上述方程計算裂縫誘導應力大小計算。
2.3 多縫應力場模型
基于第1條裂縫應力分析基礎,對下一條形成的裂縫進行應力分析,由于應力陰影效應,第2條裂縫凈壓力的模型為:
則第n條裂縫內的凈壓力為:
因此,上述力學模型中對第九條縫的應力陰影效應為前面裂縫的綜合誘導應力疊加,只需要求得前面的綜合誘導應力就可以得到鉀條縫內的凈壓力。
2.4 應力陰影求解
考慮裂縫平面上的應力變化,假設二維人工垂直裂縫(圖3)。
則在水平面上人工裂縫所引起的不同位置的誘導應力為:
由廣義Hook定律計算得:
同時,各幾何參數間存在以下關系:
利用式(1)~(6)可計算裂縫誘導應力大小,上述研究表明在地層最小水平主應力不變的情況下,誘導應力大小與裂縫內的凈壓力具有直接關系,應力陰影誘導應力對裂縫網絡擴展形態的影響具有雙重作用。苗+先,在壓力和寬度疊加的過程中,由于應力陰影的影響,附加的正應力改變每個裂縫單元原始的就地應力值,這直接影響裂縫內壓力和寬度的分布,導致裂縫的擴展發生改變。其次,考慮應力陰影誘導應力(正應力和剪應力)后,裂縫擴展端部的局部應力場發生改變,可引起局部主應力方向發生改變。局部主應力方向發生改變又會引起裂縫發生轉向,進一步影響裂縫網絡擴展形態。因此,研究應力陰影影響因素指導壓裂設計顯得尤為重要。
3 影響因素分析
以某口致密儲層X水平井為例,水平段長度1043m,儲層有效孔隙度8%,有效滲透率0.12mD,地層孔隙壓力29.0MPa,楊氏模量25000MPa,泊松比0.2,最小水平主應力43.2 MPa,抗張強度3.5MPa。分段改造設計6~10段。利用上述數學模型計算水平井10段改造下不同段內的凈壓力及裂縫形態表明:100m的段間距內裂縫將發牛轉向(圖4、5)。
3.1 楊氏模量和泊松比
為進一步研究儲層巖石力學參數對應力陰影的影響,設置不同泊松比和楊氏模量進行研究。考察泊松0.2、0.35、0.45比時,結果表明泊松比對應力陰影的影響較小,即便在泊松比為0.5的極端情況下,導致的凈壓力變化幅度仍不大,因此產生的應力影響很小(圖6)??疾鞐钍夏A?span lang="EN-US">20000MPa、25000MPa、30000MPa等3種情況,研究表明:楊氏模量對裂縫內凈壓力值的變化影響較大,產生的應力陰影現象越強烈,楊氏模量越高對裂縫應力陰影效果影響也較大(圖7)。因此,儲層巖石力學參數中楊氏模量對應力陰影現象的影響最敏感。
3.2 水平主應力差
在給定儲層條件下變化水平主應力差值,考察水平主應力差值1MPa、3MPa、6MPa等3種條件對應力陰影的影響。結果表明:隨著應力差的增大凈壓力值變小,所形成的應力陰影效應也逐級變小,儲層應力差的越小,形成的誘導應力范圍越大,對鄰近段裂縫改造的陰影效應越大,導致裂縫間的相互作用越強,促使裂縫發生轉向,形成復雜裂縫網絡。
3.3 施工排量與砂比
考察同等施工規模條件下,施工排量對應力陰影效應的影響程度,施工排量在5~12m3/min變化時,結果表明隨著排量的增大凈壓力增大,所形成的應力陰影效應的也顯著增加,因此施工排量對應力陰影效應的影響敏感性強,同樣的研究表明施工砂比對應力陰影效應的影響敏感性較弱。
3.4 分段間距的影響
在水平井單段同等規模、同等施工參數條件下,考察水平井裂縫間距為100m、125m、200m等3種不同分段裂縫間距下對應力陰影效應的影響。結果表明:在同等條件下,隨著裂縫間距的增大,形成的應力陰影效應顯著變小,對鄰段的改造裂縫的陰影效應減小,利于裂縫長度的延伸但不利十復雜裂縫網絡的形成,但過近的裂縫間距可使得縫內誘導應力過大,對施工安全不利。
4 礦場試驗
在上述優化研究的基礎上,利用文中模型及研究結淪,對本井進行分段設計與工藝參數優化,按照設計參數對本井進行現場6段大型壓裂改造施工,注入總液量4480m3、加砂320m3。微地震監測結果表明:利用人工裂縫應力陰影效應進行分段改造后,人工裂縫出現多處轉向裂縫,說明利用應力陰影效應實現了6段壓裂下最優改造體積,較以往未考慮龐力陰影下節約3~4段改造井段,降低成本30%以上,為水平井多段的體積改造優化提供重要技術指導。
5 結論
筆者借鑒Bunger提出的應力場分析假設條件,從水平井單一橫切縫入手,建立考慮應力陰影的應力場疊加數學模型并推出數學公式,在此基礎上得出簡化的水平井多段應力陰影效應應力場模型。影響因素分析表明:儲層地質參數中楊氏模量,水平主應力差對應力陰影效應的影響大,而泊松比對應力陰影效應的影響幾乎很小;施工參數中施工排量對應力陰影效應的影響明顯,在一定儲層條件下,分段及簇間距越小,應力陰影效果對水平井裂縫附近其他裂縫的擴展產生應力陰影影響越大。
利用文中研究結果對致密層水平井進行優化設計與現場施工。微地震表明結果:利用應力陰影效應只壓裂6段就實現了以往10段壓裂的改造體積,較以往未考慮應力陰影下節約3~4段,降低30%以上成本,對合理利用應力陰影效應進行水平井體積改造裂縫尺寸設置及工藝參數優化具有重要指導意義。
符號說明
pnet為裂縫內凈壓力,MPa;pc為裂縫內壓力,MPa;shmin為最小水平主壓力,MPa;sf為巖石抗張強度,MPa;sy誘導為人工裂縫引起的誘導應力,MPa;wo為裂縫縫寬;Ep為油層的楊氏模量,MPa;up為儲層的泊松比;hf為裂縫高度,m;p2net為第2條裂縫內凈壓力,MPa;s1y誘導為第1條裂縫引起的誘導應力,MPa;pnnet為第n條裂縫內凈壓力,MPa;i為變量;C為裂縫半縫高度,m;x為縫長方向距離,m;y為縫寬方向距離,m;如q、q1和q2為負,則用q+180°、q1+180°和q2+180°代替。
參考文獻
[1]趙政璋,杜金虎.致密油氣[M].北京:石油工業出版社,2012.
ZHAO Zhengzhang,DU Jinhu.Tight oil and gas[M].Beijing:Petroleum Industry Press,2012.
[2]賈承造,張永峰,趙霞.中國天然氣工業發展前景與挑戰[J].天然氣工業,2014,34(2):1-11.
JIA Chengzao,ZHANG Yongfeng,ZHAO Xia.Prospects of and challenges to natural gas industry development in China[J].Natural Gas Industry,2014,34(2):1-11.
[3]潘林華,張士誠,程禮軍,等.水平井“多段分簇”壓裂簇間干擾的數值模擬[J].天然氣工業,2014,34(1):74-79.
PAN Linhua,ZHANG Shicheng,CHENG Lijun,et al.A numerical simulation of the intor-cluster interference in multi-cluster staged fracking for horizontal wells[J].Natural Gas Industry,2014,34(1):74-79.
[4]OISON J E,DAHI T A.Modeling simultaneous growth of multiple hydraulic fractures and their interaction with natural fractures[C]//paper 119739 presented at the SPE Hydraulic Fracturing Technology Conference,19-21 January 2009,Woodlands,Texas,USA.New York:SPE,2009.
[5]BUNGER A P,ZHANG X,JEFFREY R G.Parameters affecting the interaction among closely spaced hydraulic fractures[C]//paper 140426 presented at the SPE Hydraulic Fracturing Technology Conference,24-26 January 2011,the Woodlands,Texas,USA.New York:SPE,2011.
[6]ROUSSEL N P,SHARMA M M.Strategies to minimize frac spacing and stimulate natural fractures in horizontal completions[C]//paper 146104 presented at the SPE Annual Technical Conference and Exhibition.30 Octoher-2 November 2011,Denver,Colorado,USA.New York:SPE,2011.
[7]趙金洲,任嵐,胡永全.頁巖儲層壓裂縫成網延伸的受控因素分析[J].西南石油大學學報:自然科學版,2013,35(1):1-7.
ZHAO Jinzhou,REN Lan,HU Yongquan.Controlling factors of hydraulic fractures extending into network in shale formations[J].Journal of Southwest Petroleum University:Science&Technology Edition,2013,35(1):1-7.
[8]曾凡輝,郭建春,劉恒,等.北美頁巖氣高效壓裂經驗及對中國的啟示[J].西南石油大學學報:自然科學版,2013,35(6):90-98.
ZENG Fanhui,GUO Jianchun,LIU Heng,et al.Experience of efficient fracturing of shale gas in North America and enlightenment to China[J].Journal of Southwest Petroleum University:Science&Technology Edition,2013,35(6):90-98.
[9]張濤,郭建春,劉偉.清水壓裂中支撐劑輸送沉降行為的CFD模擬[J].西南石油大學學報:自然科學版,2014,36(1):74-82.
ZHANG Tao,GUO Jianchun,LIU Wei.CFD Simulation of proppant transportation and settling in water fracture treatments[J].Journal of Southwest Petroleum University:Science&Technology Edition,2014,36(1):74-82.
[10]蔣廷學,賈長貴,王海濤,等.頁巖氣網絡壓裂設計方法研究[J].石油鉆探技術,2011,39(3):36-40.
JIANG Tingxue,JIA Changgui,WANG Haitao,et al.Study on network fracturing design method in shale gas[J].Petroleum Drilling Techniques,2011,39(3):36-40.
[11]葉源新,劉光廷.巖石滲流應力耦合特性研究[J].巖石力學與工程學報,2005,24(14):2518-2523.
YE Yuanxin,LIU Guangting.Research on coupling characteristics of fluid flow and stress within rock[J].Chinese Journal of Rock Mechanics and Engineering,2005,24(14):2518-2523.
本文作者:才博 唐邦忠 丁云宏 楊振周 何春明 唐術華 陳春霞
作者單位:中國石油勘探開發研究院
中國石油勘探丌發研究院廊坊分院
中國石油華北油田公司
中國石油吐哈油出公司溫米采油廠
中國石油吐哈油田公司部善采油廠
您可以選擇一種方式贊助本站
支付寶轉賬贊助
微信轉賬贊助