王俊超, 賈俊山, 李佳琦, 佟亮, 李嘉成
(1.新疆油田公司工程技術(shù)研究院, 克拉瑪依 834000; 2.北京遠(yuǎn)望景盛科技發(fā)展有限公司, 北京 100083)
中國(guó)頁巖油氣資源豐富,近年來實(shí)現(xiàn)了頁巖氣的商業(yè)化開采,形成了“水平井+體積壓裂”的施工工藝[1-2]?;谒毫褜?shí)驗(yàn)[3-4]和數(shù)值模擬[5-6],目前人們對(duì)水力裂縫的擴(kuò)展和延伸已經(jīng)形成了相對(duì)成熟的認(rèn)識(shí)。由于中國(guó)初期對(duì)頁巖氣開發(fā)認(rèn)識(shí)不足,現(xiàn)有水力壓裂水平井間距較大,進(jìn)而導(dǎo)致井間資源動(dòng)用難度大的問題[7]。有資料表明,目前中國(guó)蜀南地區(qū)頁巖氣藏水平井井距介于 400~500 m,涪陵地區(qū)初始井距約 600 m[8]。為了降低水平井的間距,提高頁巖儲(chǔ)層資源的動(dòng)用率,降低產(chǎn)能衰竭,需要對(duì)已有井網(wǎng)部署加密井[9-10]。
水平井加密時(shí)應(yīng)考慮已有水平井(也稱為老井)開發(fā)的影響。老井長(zhǎng)期的開發(fā)會(huì)導(dǎo)致井周孔隙壓力和地層應(yīng)力發(fā)生變化。而加密井壓裂時(shí)水力裂縫的擴(kuò)展路徑與應(yīng)力狀態(tài)密切相關(guān),因此在考慮加密井布置時(shí),不可忽略老井持續(xù)開采對(duì)地層應(yīng)力的干擾,進(jìn)而影響到對(duì)加密井裂縫的形態(tài)。針對(duì)油氣開發(fā)對(duì)地層應(yīng)力的影響問題,前人開展了一系列的研究。趙凱等[11]基于均質(zhì)各向同性彈性儲(chǔ)集層,認(rèn)為儲(chǔ)層壓力衰竭將導(dǎo)致地層水平應(yīng)力降低。孟憲波等[12]基于多場(chǎng)耦合方法,計(jì)算了致密儲(chǔ)層開采過程中的地應(yīng)力變化規(guī)律。Yang等[13]和Guo等[14]采用有限元方法,研究了頁巖氣井生產(chǎn)過程中不同水力裂縫參數(shù)及巖石力學(xué)參數(shù)對(duì)儲(chǔ)層應(yīng)力變化的影響。綜上可知隨著水平井內(nèi)頁巖氣的衰竭開發(fā),水平井周圍的地應(yīng)力會(huì)發(fā)生變化。有限元法[15]或邊界元方法[16]均證明均質(zhì)氣藏多裂縫水平井生產(chǎn)過程中會(huì)對(duì)地層應(yīng)力場(chǎng)產(chǎn)生影響,且地層應(yīng)力的變化會(huì)導(dǎo)致水力裂縫的偏轉(zhuǎn)。因此在部署加密井進(jìn)行水力壓裂作業(yè)時(shí)不能忽略老井開采的影響。
目前人們對(duì)水平井開采過程中的應(yīng)力變化規(guī)律取得了初步的認(rèn)識(shí),認(rèn)為頁巖儲(chǔ)層衰竭開發(fā)會(huì)導(dǎo)致地層的水平地應(yīng)力降低。而地層應(yīng)力是控制水力裂縫擴(kuò)展的關(guān)鍵因素。學(xué)者們對(duì)老井開發(fā)對(duì)加密井水力裂縫擴(kuò)展進(jìn)行了探究性工作。Ashish等[17]采用位移不連續(xù)法分析了老井生產(chǎn)對(duì)加密井水力壓裂裂縫的影響,結(jié)果表明,加密時(shí)機(jī)越晚,加密井水力裂縫非對(duì)稱性越嚴(yán)重??梢悦鞔_的是老井生產(chǎn)會(huì)造成地應(yīng)力的擾動(dòng),且地應(yīng)力的變化影響水力裂縫的擴(kuò)展。老井生產(chǎn)引起的地應(yīng)力演化問題屬于典型的流固耦合問題,前人對(duì)此取得了較為成熟的認(rèn)識(shí)。但是針對(duì)水力裂縫在老井生產(chǎn)引起的地應(yīng)力變化條件下的擴(kuò)展規(guī)律的研究尚不成熟。
因此,現(xiàn)采用塊體離散元方法,基于流固耦合理論,對(duì)不同地應(yīng)力梯度條件下的水力裂縫擴(kuò)展行為進(jìn)行模擬計(jì)算,分析地應(yīng)力變化對(duì)水力裂縫的影響規(guī)律,明確不同因素(施工參數(shù)和地應(yīng)力條件)對(duì)水力裂縫擴(kuò)展的影響。基于此嘗試建立老井生產(chǎn)引起的地應(yīng)力變化對(duì)加密井水力壓裂的影響特點(diǎn)。最后將該認(rèn)識(shí)與現(xiàn)場(chǎng)實(shí)例進(jìn)行對(duì)比。以期為頁巖氣開發(fā)過程中的水平井井網(wǎng)加密提供理論支撐和技術(shù)參考。
在巖土計(jì)算力學(xué)中,由于離散單元法在表達(dá)真實(shí)節(jié)理巖體中具有優(yōu)勢(shì),可以更好地表征節(jié)理發(fā)育巖體的非連續(xù)性,因此廣泛應(yīng)用于邊坡處理和節(jié)理巖體的滲流問題。在塊體離散元方法中可以通多離散隨機(jī)裂縫(discrete fracture network,DFN)模塊生成隨機(jī)裂縫網(wǎng)絡(luò),可以更真實(shí)地描述地下的復(fù)雜結(jié)構(gòu)。Zhang等[18]和鄭永香[19]采用3DEC對(duì)頁巖的水力壓裂開展了一系列研究,驗(yàn)證了該方法在水力壓裂中的適用性,因此本文不再討論該方法模擬水力壓裂的可行性,而是直接采用該方法進(jìn)行水力壓裂的模擬。
1.1.1 裂縫破壞的判定
塊體離散元方法中需要通過預(yù)置節(jié)理面來實(shí)現(xiàn)裂縫的模擬。裂縫的破壞則表示裂縫的開啟。在塊體離散元方法中,節(jié)理通過接觸來進(jìn)行描述。在塊體離散元方法中節(jié)理的基本模型為庫(kù)倫滑移節(jié)理模型,該模型在剛性塊體和柔性塊體的接觸中的原理相似。在模型中均考慮了剪切和張拉失效,同時(shí)也考慮了節(jié)理的膨脹。在彈性階段,接觸之間通過法向剛度與切向剛度來描述。節(jié)理的法向力學(xué)行為[20]可以表示為
ΔFn=-KnΔUnAc
(1)
其剪切行為可以表示為
(2)
一個(gè)原始節(jié)理(未發(fā)生滑移或開裂)的法向拉力的最大值為
Tmax=-TAc
(3)
式(3)中:T為抗拉強(qiáng)度。節(jié)理允許的最大剪切力為
(4)
式(4)中:c為節(jié)理的內(nèi)聚力;φ為摩擦角。一旦節(jié)理上的力超過其抗拉或抗剪強(qiáng)度后,接觸失效,此時(shí)節(jié)理的抗拉強(qiáng)度以及內(nèi)聚力變?yōu)?。此時(shí)節(jié)理面上的最大拉力和最大剪切力為
Tmax=0
(5)
(6)
剪切膨脹只發(fā)生在滑移模式。此時(shí),剪切位移增量可表示為
(7)
剪切位移會(huì)引起法向的位移變化,其關(guān)系可通過膨脹角ψ來表征,其關(guān)系為
ΔUn(dil)=ΔUstanψ
(8)
式(8)中:ΔUn(dil)表示由于膨脹引起的法向位移。
此時(shí)法向力應(yīng)該綜合考慮膨脹的影響,法向力變?yōu)?/p>
Fn:=Fn+KnAcΔUstanψ
(9)
1.1.2 節(jié)理內(nèi)流體流動(dòng)
節(jié)理內(nèi)的流體流動(dòng)符合N-S方程。當(dāng)兩個(gè)面近似平行的非滲透面,且流體為不可壓縮流體時(shí),N-S方程可以簡(jiǎn)化為雷諾方程:
(10)
式(10)中:u=u(xi)為平面上一點(diǎn)xi處的非滲透邊界的距離;φ為水頭高度;ρ為流體密度;μ為流體黏度;g為重力加速度;i為變分法中逗號(hào)約定的下標(biāo)。由式(10)可以得到流體的流速為
(11)
為探究地應(yīng)力對(duì)水力裂縫,建立如圖1所示的地層模型。模型中預(yù)置一條節(jié)理面作為水力裂縫的擴(kuò)展路徑。在預(yù)制裂縫兩側(cè)10 m范圍內(nèi)采用細(xì)化網(wǎng)格,在遠(yuǎn)處則采用粗化的網(wǎng)格,這樣既能保證計(jì)算精度又可提高計(jì)算效率。在預(yù)制裂縫的中心點(diǎn)設(shè)置注入點(diǎn),注入點(diǎn)位置如圖2(a)所示。為了模擬老井油氣開采引起的地應(yīng)力衰減對(duì)水力裂縫擴(kuò)展的影響,在注入點(diǎn)左側(cè)設(shè)置0.01 MPa/m的應(yīng)力梯度,右側(cè)保持原有地層應(yīng)力值,如圖2(b)所示。其他相關(guān)參數(shù)如表1所示。
圖1 地層模型圖Fig.1 Stratigraphic model
圖2 預(yù)置裂縫面示意圖Fig.2 Diagram of preset fracture surface
表1 模型參數(shù)取值表Table 1 Value of model parameter
根據(jù)確定的計(jì)算參數(shù),模擬地層內(nèi)水力裂縫的擴(kuò)展行為。根據(jù)計(jì)算過程,將所有計(jì)算步數(shù)平均分為8個(gè)階段,選擇8個(gè)階段的裂縫狀態(tài)以分析不同時(shí)間下裂縫的擴(kuò)展形態(tài),進(jìn)而分析裂縫擴(kuò)展的演化過程。計(jì)算結(jié)果如圖3所示。
圖3為隨注液時(shí)長(zhǎng)增加水力裂縫的形態(tài)演化過程。圖3中自左向右為階段1~階段8的裂縫形態(tài)。由圖3可知,水力裂縫的左側(cè)長(zhǎng)度大于裂縫右側(cè)長(zhǎng)度,這說明低地應(yīng)力區(qū)域?qū)λα芽p產(chǎn)生了吸引作用,水力裂縫向地應(yīng)力低的一側(cè)偏移。此外,由于在注入點(diǎn)左側(cè),隨著裂縫向左擴(kuò)展,地層的地應(yīng)力逐漸降低,因此裂縫向左擴(kuò)展的趨勢(shì)高于右側(cè)。圖4中分析了裂縫左右兩側(cè)長(zhǎng)度的差值,由圖3可知,在階段1時(shí),裂縫的左右長(zhǎng)度差值為0.9 m,隨著裂縫的擴(kuò)展,左右差值分別為1.8、3.9、4.7、6.0、7.5、8.4、9.7 m。隨著水力裂縫的擴(kuò)展,裂縫左右差值逐漸變大,裂縫向左偏移越大。為了表征裂縫向地應(yīng)力的偏移狀態(tài),定義水力裂縫的偏移率e,可用式(12)表示。偏心率表示裂縫中心的偏移量與裂縫總長(zhǎng)度的比值。
圖3 不同時(shí)間步裂縫的擴(kuò)展情況Fig.3 Fracture propagation in different time steps
圖4 裂縫兩側(cè)的偏移情況Fig.4 Offset on both sides of the fracture
(12)
式(12)中:e為水力裂縫的偏移率;PL為水力裂縫左側(cè)的尖端的橫坐標(biāo)點(diǎn);PR為水力裂縫右側(cè)的尖端的橫坐標(biāo)點(diǎn);Pin為注入點(diǎn)的橫坐標(biāo)點(diǎn)。
不同階段的水力裂縫偏移率如圖4所示,由圖4可知,隨著裂縫的擴(kuò)展,偏移率逐漸增加,說明裂縫的中心逐漸向地應(yīng)力區(qū)偏移。由圖4可知,在初始階段,水力裂縫兩側(cè)的應(yīng)力差值低,此時(shí)水力裂縫的左右兩側(cè)差值和裂縫偏移率較低,分別為0.9 m和0.98%。在最終階段,水力裂縫兩側(cè)尖端的應(yīng)力差增加,此時(shí)水力裂縫的左右兩側(cè)差值和裂縫偏移率較高,分別為9.7 m和5.28%。
綜上,由于地層應(yīng)力的差異,水力壓裂時(shí),水力裂縫將更趨于向低地應(yīng)力區(qū)域擴(kuò)展,水力裂縫將不再以注入點(diǎn)為中心向兩側(cè)對(duì)稱擴(kuò)展,水力裂縫的中心也會(huì)向低地應(yīng)力方向偏移。隨著已有水平井頁巖氣的開發(fā),地層地應(yīng)力發(fā)生變化。此時(shí)如果在老井周圍進(jìn)行加密,則需要關(guān)注老井開采引起的地應(yīng)力場(chǎng)擾動(dòng)對(duì)新井壓裂時(shí)水力裂縫形態(tài)的影響。為了明確不同條件下水力裂縫的擴(kuò)展規(guī)律,接下來探討壓裂工藝(注液速率與壓裂液黏度)和地層應(yīng)力狀態(tài)(應(yīng)力梯度與初始應(yīng)力值)對(duì)水力裂縫擴(kuò)展的影響規(guī)律。
壓裂施工參數(shù)是人為可控因素,主要包括壓裂液黏度和注液速率,兩者是控制裂縫擴(kuò)展速度和裂縫路徑的關(guān)鍵因素。為了分析這兩個(gè)因素對(duì)裂縫偏移率的影響,分別選取壓裂液黏度為20、40、50、100、200 cp,注液速率為0.01、0.02、0.03、0.04、0.05 m3/s。對(duì)兩個(gè)因素進(jìn)行正交設(shè)計(jì),共開展25組模擬實(shí)驗(yàn)。由于不同條件下裂縫左右尖端的位置不同,為了統(tǒng)一考察,選擇裂縫左側(cè)長(zhǎng)度達(dá)到50 m為終止條件。下述研究均基于此條件的結(jié)果展開討論。
圖5為注液速率為0.01 m3/s,應(yīng)力梯度為0.2 MPa/m時(shí),不同黏度條件下水力裂縫的擴(kuò)展情況。由圖5可知,水力裂縫均不同程度的向左發(fā)生了偏移。其中,壓裂液的黏度越低,水力裂縫向左偏移的情況越明顯。圖5中云圖的不同顏色表示裂縫的開裂寬度(其中由于計(jì)算方法的問題,設(shè)置裂縫的初始開度為0.01 mm,可視為裂縫未開裂。)由圖5可知,水力裂縫的最大寬度亦向左移動(dòng)。此外,隨著黏度的增大,水力裂縫的最大寬度值增加。綜上可得,增大壓裂液黏度可以緩解地應(yīng)力梯度引起的水力裂縫偏移。
圖5 不同黏度下水力裂縫擴(kuò)展云圖Fig.5 Hydraulic fracture propagation with different viscosities
一般地需要綜合考慮壓裂液黏度和注液速率兩個(gè)因素共同作用下裂縫的擴(kuò)展行為。計(jì)算中均以水力裂縫左側(cè)擴(kuò)展長(zhǎng)度達(dá)到50 m為終止條件,因此只需考慮水力裂縫右側(cè)的擴(kuò)展長(zhǎng)度即可。圖6(a)中展示了不同黏度下、不同注液速率情況下右側(cè)水力裂縫的擴(kuò)展情況。由圖6(a)可知,隨著注液速率的增加,右側(cè)的水力裂縫長(zhǎng)度增加。這意味著增大注液速率可以降低由于地應(yīng)力梯度導(dǎo)致的水力裂縫偏移。同樣地,增大壓裂液黏度也會(huì)導(dǎo)致右側(cè)水力裂縫長(zhǎng)度的增加。因此增加壓裂液黏度和注液速率可以降低由于地應(yīng)力梯度所導(dǎo)致的水力裂縫偏移。此外,由圖6(a)可知,當(dāng)黏度較低,注液速率較小時(shí),兩者對(duì)裂縫右側(cè)長(zhǎng)度的影響較為明顯;當(dāng)黏度和注液速率較大時(shí),其對(duì)水力裂縫右側(cè)擴(kuò)展長(zhǎng)度的影響較小。圖6(b)為不同壓裂液黏度和注液速率條件下水力裂縫在應(yīng)力梯度作用下的偏移率。由圖6可知,隨著注液速率的增大,水力裂縫的偏移率降低。從圖6中縱向分析,可知隨著壓裂液黏度的增加,水力裂縫的偏移率降低,同時(shí),在注液速率較低時(shí)(如0.01 m3/s),黏度對(duì)偏移率的影響程度高,當(dāng)注液速率較大時(shí)(如0.05 m3/s),黏度對(duì)偏移率的影響則十分有限。究其原因,主要是因?yàn)楫?dāng)壓裂液黏度增大,注液速率增加時(shí),裂縫內(nèi)的流體壓力增加,這在一定程度上抵消了由于地應(yīng)力差值造成的裂縫擴(kuò)展優(yōu)勢(shì)。綜上,通過黏度和注液速率可以有效調(diào)控地應(yīng)力對(duì)水力裂縫擴(kuò)展方向的影響,提高黏度和注液速率可在一定程度上降低地應(yīng)力不均勻?qū)λα芽p擴(kuò)展的影響。
圖6 黏度和注液速率對(duì)水力裂縫的影響規(guī)律Fig.6 Influence of viscosity and injection rate on hydraulic fracture
地應(yīng)力是影響水力裂縫擴(kuò)展的直接控制因素。地層的地應(yīng)力狀態(tài)包括初始條件下的初始地應(yīng)力和后續(xù)油氣開采造成的地應(yīng)力變化。為了探究規(guī)律,將水平井生產(chǎn)造成的地應(yīng)力變化用均勻的地應(yīng)力梯度表示。為了綜合考慮地應(yīng)力梯度和初始地應(yīng)力對(duì)水力裂縫擴(kuò)展的影響,分別選擇地應(yīng)力梯度為0、0.01、0.05、0.1、0.2 MPa/m,初始地應(yīng)力值分別選擇為10、15、20、25 MPa。對(duì)兩個(gè)因素進(jìn)行正交分析,分別開展了20組模擬實(shí)驗(yàn)。同樣地,模擬中選擇裂縫左側(cè)長(zhǎng)度達(dá)到50 m為終止條件。
圖7為不同地應(yīng)力梯度下水力裂縫的擴(kuò)展情況。由圖7可知,當(dāng)?shù)貞?yīng)力梯度為0時(shí),裂縫均勻擴(kuò)展,裂縫以注入點(diǎn)所在位置豎線為對(duì)稱軸對(duì)稱分布,右側(cè)裂縫長(zhǎng)度達(dá)到50 m。隨著地應(yīng)力梯度的增加,水力裂縫右側(cè)的擴(kuò)展長(zhǎng)度逐漸降低,且水力裂縫的最大寬度所在位置向左偏移,水力裂縫最大寬度增加。因此地應(yīng)力梯度越大,水力裂縫的偏移情況越明顯。
圖7 不同地應(yīng)力梯度下水力裂縫特征Fig.7 Hydraulic fracture propagation under different crustal stress gradients
除應(yīng)力梯度外,地應(yīng)力場(chǎng)的另外一個(gè)重要因素是初始地應(yīng)力。為了綜合考慮兩者的共同作用。將不同初始地應(yīng)力和不同應(yīng)力梯度條件下的20組模擬結(jié)果如圖8所示。圖8(a)為不同情況下裂縫右側(cè)長(zhǎng)度,由圖8(a)可知,隨應(yīng)力梯度增加,右側(cè)裂縫長(zhǎng)度減小。從縱向上來看,隨初始地應(yīng)力值增加,水力裂縫右側(cè)擴(kuò)展長(zhǎng)度降低。圖8(b)為不同情況下的裂縫偏移率,由圖8(b)可知,隨應(yīng)力梯度和初始地應(yīng)力的增加,水力裂縫的偏移率增加。因此同樣的應(yīng)力梯度條件下,高地應(yīng)力地層中的水力裂縫偏移率會(huì)高于地應(yīng)力地層。
圖8 地應(yīng)力梯度和初始地應(yīng)力對(duì)水力裂縫的影響規(guī)律Fig.8 Influence of crustal stress gradient and initial in-situ stress on hydraulic fracture
基于上述模擬,在某平臺(tái)的加密井壓裂過程的進(jìn)行了驗(yàn)證。該平臺(tái)存在兩口壓后生產(chǎn)三年的生產(chǎn)井,需綜合考慮老井對(duì)新井壓裂生產(chǎn)的影響。采用FracMan生產(chǎn)動(dòng)態(tài)模擬功能,模擬老井生產(chǎn)獲得孔隙壓力變化,結(jié)果如圖9所示。由圖9可知,在兩口老井周圍由于繼續(xù)開采,井周的地層孔隙壓力下降,經(jīng)計(jì)算目前老井周圍壓力下降約10 MPa。
圖9 孔隙壓力分布Fig.9 Dsitribution of pore pressure
結(jié)合巖石力學(xué)、初始地應(yīng)力場(chǎng),多場(chǎng)耦合獲得壓前地應(yīng)力場(chǎng),用于新井壓裂模擬,擬合得到的新井壓前最小主應(yīng)力場(chǎng)分布如圖10所示。由圖10可知,老井周圍的應(yīng)力場(chǎng)發(fā)生了變化,老井周圍的最小主應(yīng)力降低。
圖10 新井壓前水平最小主應(yīng)力Fig.10 Minimum horizontal principal stress before fracturing in new wells
基于新井壓前地應(yīng)力場(chǎng),模擬新井壓裂改造過程,模擬得到的結(jié)果如圖11所示。由圖11可知,老井壓后生產(chǎn)造成的應(yīng)力低值區(qū)對(duì)新井裂縫擴(kuò)展起到吸引作用。在圖11中選取部分典型特征裂縫進(jìn)行分析。不同裂縫的特征如下:①壓裂段距離應(yīng)力低值區(qū)較遠(yuǎn),兩側(cè)均衡擴(kuò)展;②壓裂段緊鄰應(yīng)力低值區(qū),壓裂縫網(wǎng)向低值區(qū)擴(kuò)展;③單簇射孔緊鄰應(yīng)力低值區(qū),單簇壓裂縫網(wǎng)向低值區(qū)擴(kuò)展;④受前段裂縫擴(kuò)展影響,壓裂縫遠(yuǎn)離鄰井壓裂段擴(kuò)展。
圖11 新井壓裂裂縫網(wǎng)絡(luò)分布情況Fig.11 Distribution of fracture network after fracturing in new wells
結(jié)合現(xiàn)場(chǎng)實(shí)例,可以取得以下認(rèn)識(shí):老井壓后生產(chǎn)后,會(huì)造成近壓裂縫網(wǎng)的應(yīng)力低值區(qū);新井壓裂改造過程,應(yīng)力低值區(qū)對(duì)新井裂縫擴(kuò)展起到吸引作用;老井生產(chǎn)造成應(yīng)力低值區(qū),多造成裂縫非對(duì)稱擴(kuò)展,造成改造不充分。因此在實(shí)際生產(chǎn)中應(yīng)該結(jié)合現(xiàn)場(chǎng)的地應(yīng)力條件,調(diào)控不同的壓裂施工參數(shù)(壓裂液黏度和注液速率),來降低非均勻地應(yīng)力對(duì)新井水力裂縫擴(kuò)展帶來的影響。
通過塊體離散元數(shù)值模擬方法,探究了非均勻地應(yīng)力條件下水力裂縫的擴(kuò)展規(guī)律。得出如下結(jié)論。
(1)應(yīng)力低值區(qū)會(huì)對(duì)新井壓裂的水力裂縫起到吸引作用,水力裂縫的擴(kuò)展向低應(yīng)力區(qū)域偏移。
(2)壓裂液黏度和注液速率增加,水力裂縫在非均勻地應(yīng)力條件下的偏移率降低。因此通過黏度和注液速率可以有效調(diào)控地應(yīng)力對(duì)水力裂縫擴(kuò)展方向的影響,提高黏度和注液速率可降低地應(yīng)力不均勻?qū)λα芽p擴(kuò)展的影響。
(3)地應(yīng)力對(duì)水力裂縫擴(kuò)展具有重要影響,地應(yīng)力梯度和初始地應(yīng)力的增大,水力裂縫的偏移率增加。因此在加密井壓裂時(shí)應(yīng)充分考慮初始地應(yīng)力場(chǎng)和由于老井開采引起的地應(yīng)力變化對(duì)新井壓裂時(shí)裂縫網(wǎng)絡(luò)形成的影響。