摘要：水平井組開發(fā)致密油藏時(shí)，孔隙壓力大幅下降會顯著改變井間地應(yīng)力的大小與方向，迫使加密井壓裂裂縫朝向老井壓力衰竭區(qū)擴(kuò)展，導(dǎo)致加密井產(chǎn)量受限。因此，掌握致密油藏井間地應(yīng)力場時(shí)空演化特征是加密井布井與壓裂時(shí)機(jī)優(yōu)選的重要基礎(chǔ)。針對新疆油田金龍2 井區(qū)致密油藏，基于巖石多孔彈性理論，建立了雙水平井組應(yīng)力滲流耦合模型，數(shù)值模擬得到隨孔隙壓力下降，井間地應(yīng)力場的動態(tài)變化規(guī)律。以水平主應(yīng)力方向發(fā)生反轉(zhuǎn)時(shí)的孔隙壓力變化值為評價(jià)指標(biāo)，開展了儲層物性和壓裂完井參數(shù)對地應(yīng)力時(shí)空演化特征的影響規(guī)律。結(jié)果表明，隨孔隙壓力下降，井間最大水平主應(yīng)力方向發(fā)生反轉(zhuǎn)；儲層基質(zhì)滲透率、井間距減小，單井單位儲層厚度日產(chǎn)量及壓裂縫間距增加會加快主應(yīng)力反轉(zhuǎn)，而初始儲層孔隙壓力對水平主應(yīng)力變化趨勢無影響，其中主要影響因素為井間距。

關(guān)鍵詞：致密油藏；水力壓裂；滲流；孔隙壓力；地應(yīng)力場

引言

水平井組多級壓裂被公認(rèn)為是高效開發(fā)致密油藏的關(guān)鍵技術(shù)之一，其有效提高了油藏控制面積和最終采收率[1 4]。然而，生產(chǎn)實(shí)踐表明，開采數(shù)年后，在水平井組間布置加密水平井，其壓裂增產(chǎn)效果普遍低于老井[5]。其原因是孔隙壓力大幅下降會顯著改變井間地應(yīng)力的大小與方向，迫使加密井壓裂裂縫朝向老井壓力衰竭區(qū)擴(kuò)展，導(dǎo)致加密井產(chǎn)量受限[6 7]。因此，有必要建立耦合孔隙壓力變化的地應(yīng)力計(jì)算模型，掌握致密油藏井間地應(yīng)力場時(shí)空演化特征，優(yōu)選加密井布井與壓裂時(shí)機(jī)。

Lindsay 等率先統(tǒng)計(jì)北美主要頁巖油氣盆地加密井產(chǎn)量數(shù)據(jù)，發(fā)現(xiàn)約70%～80% 加密井的單位井筒長度產(chǎn)量不及老井產(chǎn)量的50%，引發(fā)國內(nèi)外學(xué)者關(guān)注，揭開了加密井地應(yīng)力場與油藏?cái)?shù)值模擬研究的序幕[5]。Roussel 等建立了水平井組開發(fā)致密儲層流固耦合油藏模型和壓裂裂縫擴(kuò)展模型，發(fā)現(xiàn)井間應(yīng)力場隨生產(chǎn)時(shí)間延長發(fā)生轉(zhuǎn)向，且局部水平主應(yīng)力差逐漸降低，迫使加密井裂縫發(fā)生非平面擴(kuò)展[8]。Guo 等針對北美Eagle Ford 頁巖油藏，建立了可描述油水兩相滲流與巖石變形多場耦合的三維地應(yīng)力有限元模型，發(fā)現(xiàn)水平井組生產(chǎn)12 個(gè)月后井間地應(yīng)力場顯著改變，18 個(gè)月后水平主應(yīng)力方向反轉(zhuǎn)90? [9]。Zhu 等針對涪陵頁巖氣藏，建立了考慮天然裂縫的加密井壓裂裂縫擴(kuò)展模型，發(fā)現(xiàn)生產(chǎn)約5 a 后井間水平主應(yīng)力差增加，加密井的壓裂裂縫復(fù)雜度相較老井大幅降低，并獲得微地震監(jiān)測信號的證實(shí)[10 11]。Yang 等基于多孔彈性理論建立了致密儲層應(yīng)力滲流耦合模型，同時(shí)，利用離散裂縫擴(kuò)展模型模擬加密井壓裂裂縫擴(kuò)展行為，發(fā)現(xiàn)井間水平主應(yīng)力下降源于加密井裂縫朝向老井壓力衰竭區(qū)擴(kuò)展的力學(xué)機(jī)制[12]。Rezaei 等基于所建立的瞬態(tài)全耦合孔隙彈性模型，研究了泵注壓力、井間距、壓裂縫間距及水平主應(yīng)力差等因素對加密井壓裂裂縫擴(kuò)展行為的影響，結(jié)果表明，初始水平主應(yīng)力差小、老井壓裂縫間距小及加密井泵注壓力低均有利于延緩加密井裂縫擴(kuò)展進(jìn)入老井壓力衰竭區(qū)[13]。Zheng 等針對北美Haynesville 頁巖氣藏，建立了可綜合描述地質(zhì)構(gòu)造、油藏應(yīng)力滲流耦合、水力壓裂裂縫擴(kuò)展等多物理屬性與過程的加密井油藏?cái)?shù)值模擬方法，證實(shí)了加密井單位井筒長度的產(chǎn)量低于老井，發(fā)現(xiàn)增加井間距可降低加密井與老井產(chǎn)量差距[14]。Sangnimnuan 等針對立體水平井組開發(fā)北美Permain 盆地頁巖油，建立了可描述地質(zhì)分層的三維油藏地質(zhì)力學(xué)模型，發(fā)現(xiàn)在未開發(fā)的層系布置加密井，可降低老井對加密井的應(yīng)力干擾和產(chǎn)量[15]。Kumar 等利用地質(zhì)力學(xué)油藏模擬器，耦合三維水力壓裂模擬器，模擬發(fā)現(xiàn)加密井裂縫在非均勻地應(yīng)力場作用下，呈現(xiàn)出顯著的非對稱擴(kuò)展特征[16]。張金才等基于地應(yīng)力與巖性分析，提出準(zhǔn)確確定水平主應(yīng)力的大小與方向并將水平井布設(shè)在有利方位與有利層位，可以減小地應(yīng)力對壓裂的不利影響[17]。位云生等分析了威遠(yuǎn)區(qū)塊壓裂平臺的6 口水平井生產(chǎn)數(shù)據(jù)，表明水平井壓裂段長度、改造段數(shù)及簇?cái)?shù)與頁巖氣產(chǎn)量呈明顯正相關(guān)關(guān)系[18]。趙金洲等采用數(shù)值模擬方法研究裂縫對原始地應(yīng)力產(chǎn)生的影響，結(jié)果表明，張開的水力裂縫會在其周圍產(chǎn)生誘導(dǎo)應(yīng)力，壓裂液的濾失則會導(dǎo)致地層孔隙壓力變化，裂縫內(nèi)凈壓力越高、壓裂裂縫越長，水平應(yīng)力差異系數(shù)越小[19 20]。曾青冬等采用數(shù)值模擬的方法研究壓裂裂縫對裂縫間的應(yīng)力干擾，認(rèn)為小裂縫間距、長裂縫和小主應(yīng)力差更有利于縫間主應(yīng)力反轉(zhuǎn)和網(wǎng)狀裂縫的形成[21 24]。郭旭洋等綜合儲層滲流、巖石骨架線彈性變形和現(xiàn)場資料建立了頁巖油儲層三維地應(yīng)力時(shí)空演化數(shù)學(xué)模型，結(jié)果表明水平井投產(chǎn)后，縱向上會影響3 倍儲層改造厚度內(nèi)儲層的地應(yīng)力大小和方向[25]。以往研究主要集中在加密井壓裂裂縫擴(kuò)展與油藏?cái)?shù)值模擬方面，仍缺乏井間水平主應(yīng)力大小與方向的時(shí)空演化特征研究。

本文結(jié)合新疆油田金龍2 井區(qū)致密油藏典型開發(fā)參數(shù)，建立了水平井組壓裂裂縫的應(yīng)力滲流耦合模型，重點(diǎn)研究水平主應(yīng)力差及其應(yīng)力方向隨孔隙壓力降低的變化規(guī)律及其主控因素，旨在闡釋水平井組開發(fā)致密油藏井間地應(yīng)力場時(shí)空演化特征，為加密井布井與壓裂時(shí)機(jī)優(yōu)選提供理論指導(dǎo)。

1 水平井組井間地應(yīng)力場模型

建立如圖1 所示的物理模型，兩口水平井平行布置開發(fā)致密油藏，密切割壓裂形成均勻分布的橫切裂縫。模型幾何參數(shù)包括壓裂縫半長、壓裂縫間距及井間距，以井間中心點(diǎn)C 點(diǎn)為原點(diǎn)，建立XY直角坐標(biāo)系，將模型中心點(diǎn)記為C，并在X 軸間隔25.00 m 選取3 個(gè)特征點(diǎn)，在Y 軸間隔20.00 m 選取3 個(gè)特征點(diǎn)，分別記為H1～H3、V1～V3，旨在描述井間地應(yīng)力場隨孔隙壓力變化的時(shí)空演化特征。合理簡化模型，假設(shè)巖石處于二維平面應(yīng)變狀態(tài)，忽略上覆巖層壓實(shí)作用所帶來的巖石骨架應(yīng)力變化；兩水平井以相同定產(chǎn)量生產(chǎn)，且壓裂縫長均等，半縫長均為70 m；基質(zhì)與壓裂裂縫構(gòu)成雙重滲流通道，考慮到基質(zhì)中的流動通道由孔隙和天然微裂縫構(gòu)成，采用等效滲透率表征基質(zhì)孔隙和天然微裂縫的滲流能力，且符合達(dá)西流動定律。

基于孔隙彈性理論，建立孔隙-裂基于孔隙彈性理論，建立孔隙裂縫滲流與巖石變形耦合模型，其控制方程包括孔隙裂縫流體流動連續(xù)性方程、巖石應(yīng)力平衡方程，控制方程如下。

1）孔隙-裂縫流體流動連續(xù)性方程

假設(shè)孔隙-裂縫內(nèi)流體為單相飽和流體，孔隙介質(zhì)孔隙度為?。如式（1）所示，考慮控制體體積V表面積為S ，連續(xù)性方程為控制體體積V 內(nèi)的質(zhì)量增量與表面質(zhì)量增量相等。流體流動符合達(dá)西定律，如式（2）所示

式中：

ρ—孔隙流體密度，kg/m3；

?—孔隙度，%；

t—時(shí)間，s；

V—控制體體積，m3；

S—控制體表面積，m2；

n—控制體表面法向張量；

v—流體滲流速度張量，m/s；

K—基質(zhì)滲透率，mD；

μ—流體黏度，Pa·s；

p—流動壓力，Pa。

2）巖石應(yīng)力平衡方程

儲層由巖石骨架與流體組成，流體分布在骨架孔隙和裂縫中?；谔郴行?yīng)力原理，確定任意時(shí)刻巖石應(yīng)力平衡方程為

式中：

σ—有效應(yīng)力，Pa；

pp—孔隙壓力，Pa；

I—單位矩陣；

δ\"—虛應(yīng)變率，s?1；

g—單位面積的牽引力，N；

δV—虛應(yīng)變，Pa；

f—單位體積的體積力，N。

模型采用有限元法求解，儲層巖石基質(zhì)選用四邊形網(wǎng)格劃分，壓裂裂縫選用寬度為6 mm 的矩形網(wǎng)格表征。對網(wǎng)格節(jié)點(diǎn)施加原始地應(yīng)力場，模型邊界固定法向位移且無滲透性，布置壓裂裂縫前使得儲層巖石處于應(yīng)力平衡狀態(tài)，從而考慮了壓裂裂縫對原始地應(yīng)力場的擾動作用。模型參數(shù)取值依據(jù)新疆油田金龍2 區(qū)塊致密油藏水平井組開發(fā)參數(shù)設(shè)置，如表1 所示。

在基本算例中，初始孔隙壓力45 MPa，基質(zhì)滲透率1.00 mD，原始水平最大、最小水平主應(yīng)力分別為70 MPa 和60 MPa。模擬井距400 m，縫間距25 m，單井單位厚度儲層的日產(chǎn)油量30 t；假設(shè)各條壓裂裂縫產(chǎn)量均衡，將單井日產(chǎn)量平均分配至每條裂縫。改變儲層初始孔隙壓力、基質(zhì)滲透率、井間距、縫間距和產(chǎn)量等參數(shù)，開展地應(yīng)力時(shí)空演化影響因素敏感性分析。

儲層計(jì)算域設(shè)置長960.00 m× 寬570.00 m，為消除邊界效應(yīng)，始終保持計(jì)算域外邊界距離壓裂裂縫尖端210 m，即當(dāng)增加井間距或壓裂縫間距時(shí)，同步增大計(jì)算域尺寸。經(jīng)網(wǎng)格無關(guān)性檢驗(yàn)，網(wǎng)格最大步長取2.00 m，網(wǎng)格數(shù)量123 579 個(gè)。數(shù)值模擬得到孔隙壓力由初始孔隙壓力降低至0 過程中的地應(yīng)力場變化規(guī)律。

2 水平井組井間地應(yīng)力場時(shí)空演化特征

2.1 井間主應(yīng)力方向演化反轉(zhuǎn)特征

圖2 為4 個(gè)不同井底壓力條件下儲層孔隙壓力分布和水平最大主應(yīng)力場，選取4 個(gè)典型的井底壓力，分別為初始孔隙壓力的100%（初始生產(chǎn)時(shí)刻）、70%、35% 和0%。由圖2 可知，初始生產(chǎn)時(shí)，井間地應(yīng)力場已受到壓裂裂縫引起的應(yīng)力干擾的影響，不再是原始地應(yīng)力場。

圖2a 與圖2b 分別展示了井底壓力由初始孔隙壓力100%pp0=45.000 MPa 降低到0%pp0=0 MPa 時(shí)，孔隙壓力、主應(yīng)力大小及方向分布情況。由圖2b可知，由于存在張開裂縫，地應(yīng)力場在水平井投產(chǎn)前已發(fā)生擾動，隨著孔隙壓力下降，水平最大主應(yīng)力大小及方向持續(xù)改變，壓裂半縫長1 倍以內(nèi)區(qū)域偏轉(zhuǎn)顯著高于其他區(qū)域。

圖3 為井間中心線水平最大主應(yīng)力方位角與孔隙壓力間的變化關(guān)系曲線。隨著孔隙壓力的降低C、V1、H1 處水平最大主應(yīng)力方位角呈加速上升趨勢，當(dāng)水平最大主應(yīng)力方位角上升至最大值時(shí)，水平最大主應(yīng)力方位角由正值突變?yōu)樨?fù)值，水平最大主應(yīng)力方向順時(shí)針轉(zhuǎn)動90? ，即水平最大主應(yīng)力方向由y 方向變?yōu)閤 方向（或稱為水平主應(yīng)力方向反轉(zhuǎn)），例如，C 點(diǎn)處最大主應(yīng)力方向由豎直變?yōu)樗?。主?yīng)力方向反轉(zhuǎn)后，水平最大主應(yīng)力方位角仍隨孔隙壓力下降而下降，最終穩(wěn)定在0 值附近。

從圖3 中可知，H2、H3 在本文模擬條件下未發(fā)生反轉(zhuǎn)，且在水平井組投產(chǎn)時(shí)刻V2、V3 點(diǎn)處水平最大主應(yīng)力已發(fā)生反轉(zhuǎn)，故圖3b 中該兩點(diǎn)處方位角未發(fā)生明顯變化。

2.2 井間主應(yīng)力差演化及其方向反轉(zhuǎn)規(guī)律

圖4 為井間中心線X 方向4 個(gè)特征點(diǎn)處水平主應(yīng)力和水平主應(yīng)力差與孔隙壓力間的變化關(guān)系曲線。沿井間中心線X 方向，隨油藏孔隙壓力下降，水平最大主應(yīng)力逐漸降低，而水平最小主應(yīng)力逐漸升高。水平主應(yīng)力差先降低，當(dāng)達(dá)到臨界孔隙壓力時(shí)，水平主應(yīng)力差降為零，水平主應(yīng)力方向反轉(zhuǎn)，而后水平主應(yīng)力差又逐漸增大。各點(diǎn)水平主應(yīng)力差為0 時(shí)對應(yīng)的孔隙壓力與圖3 中主應(yīng)力夾角波動時(shí)對應(yīng)的孔隙壓力相同。如：井間中心點(diǎn)（C 點(diǎn)）達(dá)到臨界孔隙壓力32.62 MPa 時(shí)，水平主應(yīng)力差變?yōu)榱?，水平主?yīng)力方向發(fā)生反轉(zhuǎn)。且不同位置處的水平主應(yīng)力差值和地應(yīng)力方向變化趨勢均不相同。特征點(diǎn)越靠近井間中心點(diǎn)，水平主應(yīng)力方向發(fā)生反轉(zhuǎn)時(shí)的臨界孔隙壓力值越高，相較遠(yuǎn)離井間中心點(diǎn)位置水平主應(yīng)力反轉(zhuǎn)時(shí)間越早。

圖5 為井間中心線Y 方向4 個(gè)特征點(diǎn)處水平主應(yīng)力和水平主應(yīng)力差與孔隙壓力間的變化關(guān)系曲線。從圖5 可以看出，隨孔隙壓力下降，水平最大主應(yīng)力先減小后增加，水平最小主應(yīng)力先增加后減小。其水平主應(yīng)力差存在臨界孔隙壓力，對應(yīng)的水平主應(yīng)力差降至0，而后水平主應(yīng)力方向發(fā)生反轉(zhuǎn)。越遠(yuǎn)離井間中心點(diǎn)處，臨界孔隙壓力值越高，即發(fā)生水平主應(yīng)力方向反轉(zhuǎn)時(shí)間越早。特別的是，靠近裂縫處的V2、V3 點(diǎn)處，水平主應(yīng)力在壓裂投產(chǎn)前已發(fā)生了反轉(zhuǎn)，且水平主應(yīng)力差隨孔隙壓力下降一直增加。

綜上所述，針對致密油藏水平井組壓裂開發(fā)模式，油藏孔隙壓力下降及其非均勻分布是誘發(fā)井間地應(yīng)力時(shí)空演化的根本原因，其作用機(jī)理符合多孔介質(zhì)孔隙彈性理論。水平主應(yīng)力差先減后增、水平主應(yīng)力方向反轉(zhuǎn)是井間地應(yīng)力時(shí)空演化的兩個(gè)重要特征，是優(yōu)選加密井布井與壓裂時(shí)機(jī)的關(guān)鍵參數(shù)，這與Zhu 等[10] 的研究結(jié)論一致。

3 地應(yīng)力時(shí)空演化影響因素敏感性分析

3.1 儲層物性對地應(yīng)力時(shí)空演化的影響

儲層初始孔隙壓力和基質(zhì)滲透率是儲層物性的兩個(gè)重要參數(shù)。圖6 對比了不同初始孔隙壓力和基質(zhì)滲透率條件下，井間中心點(diǎn)處水平主應(yīng)力差與孔隙壓力間的變化關(guān)系曲線圖。由圖6 可知，對不同的初始孔隙壓力，水平主應(yīng)力差下降速率均為0.45，當(dāng)達(dá)到臨界孔隙壓力時(shí)，水平主應(yīng)力差為0，其水平主應(yīng)力方向發(fā)生反轉(zhuǎn)，并且隨著初始孔隙壓力增加，其臨界孔隙壓力也增加。不同的是，基質(zhì)滲透率對井間地應(yīng)力演化有顯著影響，根據(jù)達(dá)西定律可知，定產(chǎn)量生產(chǎn)條件下，滲透率越低壓降越大，故基質(zhì)滲透率越低，水平主應(yīng)力差下降越快，水平主應(yīng)力方向反轉(zhuǎn)時(shí)間越早。如：當(dāng)基質(zhì)滲透率為5.00 mD 時(shí)，孔隙壓力降至23.37 MPa 時(shí)，井間中心點(diǎn)處水平主應(yīng)力方向發(fā)生反轉(zhuǎn)；但當(dāng)基質(zhì)滲透率變?yōu)?.25 mD 時(shí)，孔隙壓力降至33.11 MPa時(shí)，井間中心點(diǎn)處水平主應(yīng)力方向發(fā)生反轉(zhuǎn)。說明儲層特低滲透性是促進(jìn)水平主應(yīng)力方向反轉(zhuǎn)的重要原因之一，高滲透油藏發(fā)生水平主應(yīng)力方向反轉(zhuǎn)相對較晚。

3.2 布井方案對地應(yīng)力時(shí)空演化的影響

單井配產(chǎn)和井間距是水平井組壓裂完井開發(fā)致密油藏的兩個(gè)關(guān)鍵參數(shù)。圖7 為不同單井日產(chǎn)量和井間距條件下，井間中心點(diǎn)處水平主應(yīng)力差與孔隙壓力間的變化關(guān)系曲線圖。由圖7 可知，降低單井日產(chǎn)量，可延緩水平主應(yīng)力方向反轉(zhuǎn)時(shí)間，如：當(dāng)單井日產(chǎn)量由30.0 t 減少至7.5 t 時(shí)，臨界孔隙壓力由32.62 MPa 降至26.53 MPa 時(shí)水平主應(yīng)力發(fā)生反轉(zhuǎn)。增大井間距，會顯著減小水平主應(yīng)力下降速率，可有效延緩水平主應(yīng)力方向反轉(zhuǎn)時(shí)間，如：當(dāng)井間距為400.00 m 時(shí)，當(dāng)孔隙壓力下降至32.62 MPa 水平主應(yīng)力發(fā)生反轉(zhuǎn)；但當(dāng)井間距增大至450.00 m 時(shí)，在本文模擬條件下水平主應(yīng)力方向未發(fā)生反轉(zhuǎn)。值得注意的是，井間距過小會導(dǎo)致壓裂投產(chǎn)前井間中心點(diǎn)水平主應(yīng)力方向已經(jīng)反轉(zhuǎn)。

3.3 壓裂縫間距對地應(yīng)力時(shí)空演化的影響

壓裂縫間距是重要壓裂完井參數(shù)之一。由圖8a可知，縮小壓裂縫間距有利于推遲井間中心處水平主應(yīng)力方向反轉(zhuǎn)，如：當(dāng)壓裂縫間距由25.00 m增至30.00 m 時(shí)，臨界孔隙壓力由32.62 MPa 上升為41.00 MPa。然而，當(dāng)壓裂縫間距增大至40 m 以上時(shí)，水平主應(yīng)力差隨著孔隙壓力減小而增加，井間中心點(diǎn)水平主應(yīng)力方向在壓裂投產(chǎn)前就發(fā)生了反轉(zhuǎn)，其原因是壓裂縫整體對水平最大主應(yīng)力影響范圍隨著壓裂縫間距增大而增大。如圖8b 所示，與裂縫間距為10.00 m 相比，裂縫間距為40.00 m 的裂縫整體對壓裂投產(chǎn)前的水平最大主應(yīng)力影響已波及到井間中心點(diǎn)，使得井間中心點(diǎn)處的水平最大主應(yīng)力由垂向變?yōu)樗较颉?/p>

3.4 敏感程度分析

以標(biāo)準(zhǔn)組與相鄰組井間中心點(diǎn)處，水平主應(yīng)力方向發(fā)生反轉(zhuǎn)時(shí)的孔隙壓力相對變化值表征各因素對水平主應(yīng)力反轉(zhuǎn)影響程度[式（4）]，考察單井單位儲層厚度日產(chǎn)量、儲層基質(zhì)滲透率、壓裂縫間距、井間距和儲層初始孔隙壓力等參數(shù)的影響敏感性。如圖9 所示，數(shù)值0 表示該因素對地應(yīng)力方向反轉(zhuǎn)無影響，正數(shù)、負(fù)數(shù)分別代表地應(yīng)力方向反轉(zhuǎn)與該因素呈正相關(guān)、負(fù)相關(guān)，正負(fù)數(shù)值越大表明地應(yīng)力方向反轉(zhuǎn)對該因素越敏感。各因素對水平主應(yīng)力方向反轉(zhuǎn)的敏感性排序：井間距gt; 壓裂縫間距gt; 單井單位儲層厚度日產(chǎn)量gt; 儲層基質(zhì)滲透率gt; 儲層初始孔隙壓力。

式中：Δp1—標(biāo)準(zhǔn)組井間中心點(diǎn)應(yīng)力反轉(zhuǎn)所需降低的孔隙壓力值，MPa；

Δp2—相鄰組井間中心點(diǎn)應(yīng)力反轉(zhuǎn)所需降低的孔隙壓力值，MPa；

X1—標(biāo)準(zhǔn)組自變量值；

X2—相鄰組自變量值。

自變量包含單井單位儲層厚度日產(chǎn)量、儲層基質(zhì)滲透率、壓裂縫間距、井間距和儲層初始孔隙壓力。

4 結(jié)論

1）井間水平主應(yīng)力方向和數(shù)值變化是水平井組開發(fā)致密油藏典型特征之一。開發(fā)期內(nèi)井間水平主應(yīng)力方向發(fā)生反轉(zhuǎn)，水平主應(yīng)力差呈現(xiàn)先減小后增大的趨勢；投產(chǎn)前水平主應(yīng)力方向已反轉(zhuǎn)，水平主應(yīng)力差呈現(xiàn)一直增大的趨勢。

2）水平主應(yīng)力方向反轉(zhuǎn)時(shí)間隨時(shí)間和空間動態(tài)變化，靠近井間中心點(diǎn)和壓裂裂縫區(qū)域，水平主應(yīng)力方向反轉(zhuǎn)時(shí)間越早。

3）儲層基質(zhì)滲透率越低、縮小井間距、提高單井日產(chǎn)量、增加壓裂縫間距均會加速水平主應(yīng)力方向反轉(zhuǎn)。其敏感性排序?yàn)榫g距gt; 壓裂縫間距gt; 單井單位儲層厚度日產(chǎn)量gt; 儲層基質(zhì)滲透率gt; 儲層初始孔隙壓力。

參考文獻(xiàn)

[1] 馬新華，李熙喆，梁峰，等. 威遠(yuǎn)頁巖氣田單井產(chǎn)能主控因素與開發(fā)優(yōu)化技術(shù)對策[J]. 石油勘探與開發(fā)，2020，47（3）：555-563. doi：10.11698/PED.2020.03.11

MA Xinhua， LI Xizhe， LIANG Feng， et al. Dominatingfactors on well productivity and development strategiesoptimization in Weiyuan shale gas play， Sichuan Basin，SW China[J]. Petroleum Exploration and Development，2020， 47（3）： 555–563. doi： 10.11698/PED.2020.03.11

[2] 王永輝，盧擁軍，李永平，等. 非常規(guī)儲層壓裂改造技術(shù)進(jìn)展及應(yīng)用[J]. 石油學(xué)報(bào)，2012，33（S1）：149-158.doi：10.7623/syxb2012S1018

WANG Yonghui， LU Yongjun， LI Yongping， et al. Progressand application of hydraulic fracturing technologyin unconventional reservoir[J]. Acta Petrolei Sinica， 2012，33（S1）： 149–158. doi： 10.7623/syxb2012S1018

[3] 吳寶成，李建民，鄔元月，等. 準(zhǔn)噶爾盆地吉木薩爾凹陷蘆草溝組頁巖油上甜點(diǎn)地質(zhì)工程一體化開發(fā)實(shí)踐[J]. 中國石油勘探，2019，24（5）：679-690. doi：10.-3969/j.issn.1672-7703.2019.05.014

WU Baocheng， LI Jianmin， WU Yuanyue， et al. Developmentpractices of geology-engineering integration onupper sweet spots of Lucaogou Formation shale oil inJimsar Sag， Junggar Basin[J]. China Petroleum Exploration，2019， 24（5）： 679–690. doi： 10.3969/j.issn.1672-7703.2019.05.014

[4] ZHANG Hao， SHENG J. Optimization of horizontal wellfracturing in shale gas reservoir based on stimulated reservoirvolume[J]. Journal of Petroleum Science and Engineering，2020， 190： 107059. doi： 10.1016/j.petrol.2020.-107059

[5] LINDSAY G， MILLER G， XU T， et al. Production performanceof infill horizontal wells vs. pre-existing wells inthe major US unconventional basins[C]. SPE 189875-MS，2018. doi： 10.2118/189875-MS

[6] 王俊超，賈俊山，李佳琦，等. 老井開發(fā)對加密井水力壓裂裂縫擴(kuò)展的影響規(guī)律[J]. 科學(xué)技術(shù)與工程，2022，22（26）：11356-11363. doi：10.3969/j.issn.1671-1815.2022.26.012

WANG Junchao， JIA Junshan， LI Jiaqi， et al. Influencesof development of old wells on fracture propagation fromhydraulic fracturing of infill wells[J]. Science Technologyand Engineering， 2022， 22（26）： 11356–11363. doi： 10.-3969/j.issn.1671-1815.2022.26.012

[7] 秦勇，李保柱，胡水清，等. 瑪湖凹陷致密礫巖油藏四維地應(yīng)力場模擬研究與應(yīng)用[J]. 石油科技論壇，2022，41（2）：23-31. doi：10.3969/j.issn.1002-302x.-2022.02.003

QIN Yong， LI Baozhu， HU Shuiqing， et al. Numerical simulationof four-dimensional stress field for tight glutenitereservoir in Mahu Sag， Junggar Basin[J]. Petroleum Scienceand Technology Forum， 2022， 41（2）： 23–31. doi： 10.-3969/j.issn.1002-302x.2022.02.003

[8] ROUSSEL N P， FLOREZ H A， RODRIGUEZ A A.Hydraulic fracture propagation from infill horizontalwells[C]. SPE 166503-MS， 2013. doi： 10.2118/166503-MS

[9] GUO Xuyang， WU Kan， AN Cheng， et al. Numericalinvestigation of effects of subsequent parent-well injectionon interwell fracturing interference using reservoirgeomechanics-fracturing modeling[J]. SPE Journal， 2019，24（4）： 1884–1902. doi： 10.2118/195580-PA

[10] ZHU Haiyan， TANG Xuanhe， LIU Qingyou， et al. Complexfractures propagations of infill well based on reservoirstress evolution after long-time shale gas production[C].New York： 53rd US Rock Mechanics/Geomechanics Symposium，2019.

[11] 朱海燕，宋宇家，唐煊赫，等. 頁巖氣藏加密井壓裂時(shí)機(jī)優(yōu)化——以四川盆地涪陵頁巖氣田X1 井組為例[J]. 天然氣工業(yè)，2021，41（1）：154-168. doi：10.-3787/j.issn.1000-0976.2021.01.014

ZHU Haiyan， SONG Yujia， TANG Xuanhe， et al. Optimizationof fracturing timing of infill wells in shale gasreservoirs： A case study on Well Group X1 of Fuling ShaleGas Field in the Sichuan Basin[J]. Natural Gas Industry，2021， 41（1）： 154–168. doi： 10.3787/j.issn.1000-0976.-2021.01.014

[12] YANG L， ZHANG F， DUAN Y， et al. Effect of depletioninducedstress reorientation on infill well fracture propagation[C]. New York： 53rd US Rock Mechanics/GeomechanicsSymposium， 2019.

[13] REZAEI A， DINDORUK B， SOLIMAN M Y. On parametersaffecting the propagation of hydraulic fracturesfrom infill wells[J]. Journal of Petroleum Science and Engineering，2019， 182： 106255. doi： 10.1016/j.petrol.2019.-106255

[14] ZHENG Wei， XU Tao， BAIHLY J， et al. Advanced modelingof production induced pressure depletion impact oninfill well using cloud computation in the Haynesville[C].IPTC 19460-MS， 2019. doi： 10.2523/IPTC-19460-MS

[15] SANGNIMNUAN A， LI Jiawei， WU Kan， et al. Impact ofparent well depletion on stress changes and infill well completionin multiple layers in Permian Basin[C]. Denver：Proceedings of the 7th Unconventional Resources TechnologyConference， 2019. doi： 10.15530/urtec-2019-972

[16] KUMAR A， SHRIVASTAVA K， ELLIOTT B， et al. Effectof parent well production on child well stimulation andproductivity[C]. SPE 199700-MS， 2020. doi： 10.2118/199700-MS

[17] 張金才，亓原昌. 地應(yīng)力對頁巖儲層開發(fā)的影響與對策[J]. 石油與天然氣地質(zhì)，2020，41（4）：776-783，799.doi：10.11743/ogg20200411

ZHANG Jincai， QI Yuanchang. Impact of in-situ stresseson shale reservoir development and its countermeasures[J]. Oil amp; Gas Geology， 2020， 41（4）： 776–783， 799.doi： 10.11743/ogg20200411

[18] 位云生，齊亞東，賈成業(yè)，等. 四川盆地威遠(yuǎn)區(qū)塊典型平臺頁巖氣水平井動態(tài)特征及開發(fā)建議[J]. 天然氣工業(yè)，2019，39（1）：81-86. doi：10.3787/j.issn.1000-0976.-2019.01.009

WEI Yunsheng， QI Yadong， JIA Chengye， et al. Productionperformance of and development measures for typicalplatform horizontal wells in Weiyuan Shale Gas Field，Sichuan Basin[J]. Natural Gas Industry， 2019， 39（1）： 81–86. doi： 10.3787/j.issn.1000-0976.2019.01.009

[19] 趙金洲，許文俊，李勇明，等. 低滲透油氣藏水平井分段多簇壓裂簇間距優(yōu)化新方法[J]. 天然氣工業(yè)，2016，36（10）：63 69. doi：10.3787/j.issn.1000-0976.2016.10.-008

ZHAO Jinzhou， XU Wenjun， LI Yongming， et al. A newmethod for cluster spacing optimization of multi-clusterstaged fracturing in horizontal wells of low-permeabilityoil and gas reservoirs[J]. Natural Gas Industry， 2016，36（10）： 63–69. doi： 10.3787/j.issn.1000-0976.2016.10.-008

[20] 郭建春，李根，周鑫浩. 頁巖氣藏縫網(wǎng)壓裂裂縫間距優(yōu)化研究[J]. 巖土力學(xué)，2016，37（11）：3123-3129. doi：10.16285/j.rsm.2016.11.011

GUO Jianchun， LI Gen， ZHOU Xinhao. Optimization offracture spacing in fracture network of shale gas reservoir[J]. Rock and Soil Mechanics， 2016， 37（11）： 3123–3129. doi： 10.16285/j.rsm.2016.11.011

[21] 曾青冬，佟穎，姚軍. 考慮應(yīng)力干擾的多簇壓裂水平井產(chǎn)能分布規(guī)律[J]. 中國石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2019，43（1）：99-107. doi：10.3969/j.issn.1673-5005.2019.01.012

ZENG Qingdong， TONG Ying， YAO Jun. Production distributionin multi-cluster fractured horizontal wells accountingfor stress interference[J]. Journal of China Universityof Petroleum， 2019， 43（1）： 99–107. doi： 10.3969/j.-issn.1673-5005.2019.01.012

[22] 劉乃震，張兆鵬，鄒雨時(shí)，等. 致密砂巖水平井多段壓裂裂縫擴(kuò)展規(guī)律[J]. 石油勘探與開發(fā)，2018，45（6）：1059-1068. doi：10.11698/PED.2018.06.14

LIU Naizhen， ZHANG Zhaopeng， ZOU Yushi， et al. Propagationlaw of hydraulic fractures during multi-stagedhorizontal well fracturing in a tight reservoir[J]. PetroleumExploration and Development， 2018， 45（6）： 1059–1068.doi： 10.11698/PED.2018.06.14

[23] 李旺，李連崇，唐春安. 水平井平行裂縫間誘導(dǎo)應(yīng)力干擾機(jī)制的數(shù)值模擬研究[J]. 天然氣地球科學(xué)，2016，27（11）：2043-2053. doi：10.11764/j.issn.-1672-1926.2016.11.2043

LI Wang， LI Lianchong， TANG Chun′an. Numerical simulationresearch on mechanism of induced stress perturbationbetween parallel fractures in horizontal wells[J].Natural Gas Geoscience， 2016， 27（11）： 2043–2053. doi：10.11764/j.issn.1672-1926.2016.11.2043

[24] 李士斌，官兵，張立剛，等. 水平井壓裂裂縫局部應(yīng)力場擾動規(guī)律[J]. 油氣地質(zhì)與采收率，2016，23（6）：112-119. doi：10.13673/j.cnki.cn37-1359/te.2016.06.019

LI Shibin， GUAN Bing， ZHANG Ligang， et al. Localstress field disturbance law of horizontal well fracturing[J]. Petroleum Geology and Recovery Efficiency，2016， 23（6）： 112–119. doi： 10.13673/j.cnki.cn37-1359/te.-2016.06.019

[25] 郭旭洋，金衍，林伯韜. 頁巖油立體開發(fā)誘發(fā)的儲層地應(yīng)力動態(tài)響應(yīng)特征[J]. 石油鉆采工藝，2020，42（6）：738-744. doi：10.13639/j.odpt.2020.06.012

GUO Xuyang， JIN Yan， LIN Botao. Dynamic responsecharacteristics of reservoir in-site stress induced by thetridimensional development of shale oil[J]. Oil Drilling amp;Production Technology， 2020， 42（6）： 738–744. doi： 10.-13639/j.odpt.2020.06.012

作者簡介

張羽鵬，1996 年生，男，漢族，黑龍江大慶人，工程師，主要從事非常規(guī)油氣壓裂改造理論與應(yīng)用方面的工作。E-mail：Zhangyp2021@petrochina.com.cn

呂振虎，1990 年生，男，漢族，新疆奇臺人，工程師，主要從事油氣田儲層改造研究。E-mail：lvzhenhu2016@petrochina.com.cn

李嘉成，1994 年生，男，漢族，新疆克拉瑪依人，工程師，主要從事非常規(guī)油藏改造及地質(zhì)工程一體化研究工作。E-mail：lijc831@petrochina.com.cn

陳小璐，1995 年生，女，滿族，黑龍江虎林人，工程師，主要從事儲層壓裂改造技術(shù)研究工作。E-mail：chenxiaolu@petrochina.com.cn

盛茂，1985 年生，男，漢族，安徽池州人，教授，博士研究生導(dǎo)師，主要從事非常規(guī)油氣壓裂改造、壓裂人工智能等基礎(chǔ)理論與應(yīng)用研究。E-mail：shengmao@cup.edu.cn

編輯：牛靜靜

基金項(xiàng)目：中國石油科技創(chuàng)新基金（2018D 5007 0308）