砂巖氣藏水平井體積壓裂中Y341封隔器裂縫控制模擬

摘 要 使用Eclipse軟件建立了體積壓裂中Y341封隔器裂縫控制模型，在該模型中以不同段數(shù)、段長、裂縫夾角、裂縫半長、裂縫密度和裂縫導(dǎo)流能力作為控制因素，約束Y341封隔器對裂縫的控制能力。模擬結(jié)果顯示：在約束參數(shù)控制下，Y341封隔器可完成裂縫控制，提高致密砂巖氣藏水平井中的天然氣采集量。

關(guān)鍵詞 Y341封隔器 致密砂巖 氣藏水平井 體積壓裂 裂縫調(diào)控模擬

中圖分類號(hào) TE925.3" "文獻(xiàn)標(biāo)志碼 A" "文章編號(hào) 0254?6094（2024）05?0728?06

作者簡介：樊永濤（1989-），工程師，從事鄂爾多斯盆地延安氣田致密區(qū)域的鉆井、試氣、壓裂等研究工作，919977400@qq.com。

引用本文：樊永濤，梁東陽，張照鴻，等.砂巖氣藏水平井體積壓裂中Y341封隔器裂縫控制模擬［J］.化工機(jī)械，2024，51（5）：728-733.

致密巖層氣藏作為一種重要的天然氣資源儲(chǔ)存方式，已經(jīng)引起了天然氣行業(yè)研究者的廣泛關(guān)注。實(shí)際上，致密巖層氣藏是指以游離狀態(tài)存在于致密巖層中，但沒有長距離大規(guī)模運(yùn)移的聚集性天然氣［1］。目前，水平井體積壓裂是致密巖層氣藏開采中最常用的方法［2～4］。雖然20世紀(jì)末期就已經(jīng)對水平井體積壓裂進(jìn)行了相關(guān)研究，但當(dāng)時(shí)使用的試井模型只能假設(shè)單一主裂縫，無法模擬復(fù)雜的裂縫網(wǎng)絡(luò)，這導(dǎo)致最終研究結(jié)果與實(shí)際體積壓裂結(jié)果存在較大差異，獲取的有關(guān)裂縫網(wǎng)絡(luò)的具體參數(shù)失去了準(zhǔn)確性和研究價(jià)值［5］。水平井體積壓裂造成的裂縫網(wǎng)絡(luò)中的各個(gè)裂縫，無論是長度還是間距等參數(shù)均不相同。如果裂縫間距差異較大，會(huì)直接導(dǎo)致最終計(jì)算結(jié)果的偏差。由此可見，調(diào)控優(yōu)化裂縫相關(guān)參數(shù)是水平井體積壓裂技術(shù)中的一項(xiàng)重要工作［6～8］。

近年來，諸多研究者都對致密巖層氣藏水平井整體開發(fā)以及開裂進(jìn)行了研究。例如，胡納川等以提升開采率為目的，使用數(shù)值模擬的方式進(jìn)行模擬分析，該方法充分分析了氣藏情況以及多種開采方式對最終開采結(jié)果的影響，但其側(cè)重點(diǎn)并不是開壓裂縫方面，研究內(nèi)容也不夠充分［9］。張衍君等針對體積裂縫的漏失類型采取了對應(yīng)的控制措施，但該方法在壓裂參數(shù)優(yōu)選方面的效果較差，導(dǎo)致體積壓裂裂縫控制不佳，開采量不高［10］。余佩蓉等以瑪湖凹陷風(fēng)城組為對象分析了壓力裂縫擴(kuò)展情況，該方法結(jié)合Petrel軟件與Abaqus軟件模擬出了裂縫的擴(kuò)展情況，雖然獲得了較為準(zhǔn)確的裂縫擴(kuò)展分析結(jié)果，但并未給出針對這些裂縫的詳細(xì)控制方法［11］。綜上，裂縫的控制仍有較廣的擴(kuò)展研究空間。

筆者以Y341封隔器為研究對象，在體積壓裂過程中，完成裂縫控制過程建模；利用該模型模擬分析各裂縫參數(shù)變化下Y341封隔器對裂縫的控制過程，從而確定Y341封隔器應(yīng)用中的調(diào)控裂縫參數(shù)。

1 水平井體積壓裂中Y341封隔器裂縫調(diào)控模擬方法

1.1 研究區(qū)域概況

筆者選取我國中部某天然氣儲(chǔ)存區(qū)域作為研究區(qū)域，該區(qū)域地質(zhì)構(gòu)造復(fù)雜，面積約143.7平方千米，天然氣儲(chǔ)量約2178.53×104 t。該區(qū)域地下水資源豐富，分布有淺湖相水下分流河道、天然堤等結(jié)構(gòu)，地理分布呈現(xiàn)向西傾斜的單斜構(gòu)造特征。由于多期河道疊加，區(qū)域內(nèi)沙體呈現(xiàn)垂直向上的變化特點(diǎn)。儲(chǔ)層的主要地質(zhì)結(jié)構(gòu)由巖屑長石砂巖和長石砂巖組成，平均孔隙度為9.24%，具有較強(qiáng)的非均質(zhì)性儲(chǔ)層特點(diǎn)，該地層的平均滲透率為0.45 mD，礦化度約16 291 mg/L，主要水質(zhì)類型為CaCl2。水平井地質(zhì)參數(shù)如下：

水平井儲(chǔ)層壓力 70.6 MPa

水平井段長度 1 614 m

水平井儲(chǔ)層壓力系數(shù) 1.7

水平井泊松比 0.4

水平井最大水平主應(yīng)力 99.2 MPa

水平井最小水平主應(yīng)力 91.8 MPa

水平井地層溫度 128 ℃

水平井垂向水平主應(yīng)力 108.5 MPa

水平井彈性模量 46.7 GPa

1.2 Y341封隔器裂縫擴(kuò)展路徑控制模型設(shè)計(jì)

1.2.1 復(fù)雜裂縫網(wǎng)絡(luò)結(jié)構(gòu)數(shù)據(jù)

體積壓裂導(dǎo)致的氣藏區(qū)域裂縫網(wǎng)絡(luò)是一種隨機(jī)性強(qiáng)且形態(tài)各異的復(fù)雜網(wǎng)絡(luò)結(jié)構(gòu)，使用微地震監(jiān)測設(shè)備獲得的結(jié)果顯示，裂縫網(wǎng)絡(luò)可以分成中間窄兩端寬、等寬或兩端窄中間寬3種形態(tài)［12］。因此，在構(gòu)建裂縫網(wǎng)絡(luò)數(shù)據(jù)模型時(shí)，可以分別使用雙曲形、橢圓形或矩形等形狀來描述裂縫數(shù)據(jù)。按照主次級(jí)別也可以把裂縫網(wǎng)絡(luò)劃分成主裂縫、支裂縫和次裂縫3種類型，此時(shí)支裂縫與水平井段所在方向平行，另兩種裂縫則與水平井段所在方向垂直。

1.2.2 Y341封隔器裂縫擴(kuò)展路徑控制模型構(gòu)建

將待研究的砂巖區(qū)域看作數(shù)個(gè)基質(zhì)塊體單元，彈簧在各塊體單元之間發(fā)揮連接作用，模型中的彈簧出現(xiàn)斷裂則代表致密砂巖出現(xiàn)較為明顯的斷裂，即砂巖中出現(xiàn)裂縫［13，14］。壓裂液在每兩個(gè)塊體單元（Y341封隔器）之間的流動(dòng)情況使用一個(gè)裂縫單元表征，裂縫承載的外部荷載是壓裂液導(dǎo)致的Y341封隔器壓力，利用有限元法求解各塊體單元所發(fā)生的變形，使用離散元法計(jì)算彈簧的斷裂情況，最終構(gòu)建的Y341封隔器裂縫擴(kuò)展路徑控制模型如圖1所示。

在使用體積壓裂開采天然氣時(shí)，模型構(gòu)建不需要考慮超低滲透率的影響，Y341封隔器裂縫控制中的關(guān)鍵參數(shù)只考慮應(yīng)力、儲(chǔ)層壓力及天然裂縫等［15］。然后，對壓裂裂縫表征，評估裂縫的方向以及Y341封隔器分布法向應(yīng)力和切向應(yīng)力等，從而確定最佳的裂縫控制條件，提高油氣產(chǎn)能和油氣采收率。

1.2.3 模型裂縫控制過程模擬

裂縫單元上隨機(jī)一點(diǎn)記為S，兩個(gè)Y341封隔器塊體單元與裂縫單元相對的兩個(gè)面分別記為裂縫面1和裂縫面2，這兩個(gè)裂縫面上的隨機(jī)兩點(diǎn)記為S和S。裂縫面上任一點(diǎn)沿x方向和y方向的變形量分別記為Δu和Δv，用u、u與v、v分別表示這兩個(gè)變形在裂縫面1和裂縫面2上的響應(yīng)變量值，其關(guān)系如下：

點(diǎn)S在法向和切向兩個(gè)方向上的變形量用Δv和Δu表示，把坐標(biāo)整體向沿著邊界結(jié)構(gòu)建立的坐標(biāo)系轉(zhuǎn)換，則可以得到：

其中，θ表示坐標(biāo)系旋轉(zhuǎn)角度。

裂縫面隨機(jī)一點(diǎn)的法向應(yīng)力σ與切向應(yīng)力τ計(jì)算式為：

其中，k、k表示兩個(gè)方向上的抗裂系數(shù)；u、u表示裂縫面上的兩個(gè)切向應(yīng)力分量；v、v表示裂縫面上的兩個(gè)法向應(yīng)力分量。

以上完成裂縫控制路徑模擬。

1.3 模型模擬中的裂縫控制條件

通過上述模型設(shè)計(jì)，并結(jié)合地質(zhì)和巖石物理數(shù)據(jù)，進(jìn)行水平井體積壓裂的Y341封隔器裂縫控制，控制條件約束原則如下：

a. 裂縫路徑控制中的數(shù)量限制。致密砂巖氣藏水平井受到水平壓裂作用影響，地下內(nèi)部結(jié)構(gòu)中會(huì)產(chǎn)生較多的復(fù)雜裂縫，這些復(fù)雜裂縫會(huì)與自然環(huán)境中的現(xiàn)有裂縫相互連通，導(dǎo)致裂縫數(shù)量增多，使得致密砂巖氣藏發(fā)生較為嚴(yán)重的滲透［16］。使用Y341封隔器體積壓裂裂縫擴(kuò)展模型開展計(jì)算時(shí)，為了使計(jì)算結(jié)果更加精簡，按照裂縫情況根據(jù)Y341封隔器特點(diǎn)，將復(fù)雜裂縫網(wǎng)絡(luò)劃分成高滲透區(qū)域與主裂縫區(qū)域。按照水平井水平段特點(diǎn)，將水平井劃分成多個(gè)不同的段，每段按照裂縫的情況劃分成多個(gè)不同的簇，2簇、3簇、4簇的帶寬分別為60、70、80 m，經(jīng)過劃分后的高滲透區(qū)域滲透率為3×10-3 μm2。

b. 劃分裂縫網(wǎng)絡(luò)的區(qū)域。開展體積壓裂工作時(shí)，壓裂工具的具體布置方式是按照固定壓裂段長度，在水平井上均勻布置Y341封隔器，通常各Y341封隔器之間的距離設(shè)置為10 m。水平井的總長度為900 mm，在此基礎(chǔ)上，分析每2個(gè)簇至4個(gè)簇范圍內(nèi)，裂縫長度或壓裂段數(shù)對Y341封隔器壓裂效果的影響情況。

c. 裂縫控制中需要考慮水平井的水平段長度。Y341封隔器裂縫控制時(shí)，水平段長度與水平井產(chǎn)能呈現(xiàn)出正比例關(guān)系，但是相比直井，天然氣在水平井中產(chǎn)生的摩擦阻力更高，這種阻力會(huì)導(dǎo)致水平段長度增大，同時(shí)間接造成天然氣產(chǎn)量降低。因此，為實(shí)現(xiàn)水平井體積壓裂裂縫調(diào)控，在綜合考慮節(jié)約成本與優(yōu)化施工技術(shù)的前提下優(yōu)化水平井的水平段長度，從而實(shí)現(xiàn)Y341封隔器對體積壓裂裂縫的有效控制。設(shè)定水平段長度分別為450、550、650、750 m，在模型中利用Eclipse軟件模擬出不同水平段長度、體積壓力下的天然氣開采量變化情況。

d. 優(yōu)化體積壓力裂縫夾角。利用Y341封隔器壓裂形成的裂縫之間會(huì)形成夾角，在模型中使用Eclipse軟件自帶的LGR方法對裂縫之間形成的夾角開展計(jì)算，然后在模型中優(yōu)化水平井Y341封隔器體積壓裂中的裂縫夾角。

e. 優(yōu)化體積壓力裂縫長度。在模型中模擬調(diào)控裂縫的平均長度分別在110、160、210、260 m，最終確定出合適的調(diào)控優(yōu)化體積壓裂裂縫的裂縫長度。

f. 優(yōu)化體積壓力裂縫密度。統(tǒng)計(jì)1 000 m水平井中的裂縫數(shù)量，將該數(shù)量表征體積壓力裂縫密度，以4.5 a內(nèi)的天然氣開采量為評價(jià)標(biāo)準(zhǔn)，根據(jù)模型模擬的分析結(jié)果優(yōu)化裂縫密度。

g. 優(yōu)化體積壓力裂縫導(dǎo)流能力。導(dǎo)流能力是影響體積壓力裂縫、改變天然氣產(chǎn)量的關(guān)鍵因素。按照研究經(jīng)驗(yàn)，設(shè)定裂縫導(dǎo)流能力分別為15、25、35 D×m（1 D=9.869233×10-13 m2），在模型中模擬不同裂縫導(dǎo)流能力下的天然氣開采情況。

2 模擬實(shí)驗(yàn)與結(jié)果分析

2.1 模擬參數(shù)設(shè)置

Y341封隔器的模擬參數(shù)如下：

最大外徑 114 mm

最小內(nèi)徑 75 mm

坐封壓力 10～15 MPa

最大工作壓力 35 MPa

解封負(fù)荷 40～60 kN

工作溫度 不大于120 ℃

兩端連接螺紋尺寸 3 1/2UP（加厚螺紋，直徑88.9 mm）

Y341封隔器模型如圖2所示。

其他相關(guān)的水平井體積壓裂裂縫控制模擬參數(shù)如下：

滲透率 0.51×10-3 μm2

注液速度 20 m3/min

最大注液壓力 40 MPa

水平井半縫長 240 m

水平井帶寬 165 m

頁巖彈性模量 50 GPa

泊松比 0.25

在上述模型參數(shù)下，采用Eclipse軟件進(jìn)行Y341封隔器裂縫控制模擬。

2.2 模擬結(jié)果分析

2.2.1 水平井段數(shù)優(yōu)化效果

水平井段設(shè)置為2簇，水平井體積壓裂中Y341封隔器裂縫控制的段數(shù)分別設(shè)置為5段、10段、15段，模擬研究區(qū)域的天然氣采集量隨時(shí)間的變化關(guān)系如圖3所示。由圖3可知，隨著時(shí)間的增加，天然氣采集量呈上升趨勢，尤其在天然氣采集初期，各體積壓裂段數(shù)下均呈現(xiàn)出天然氣采集量迅速上升的趨勢，此后上升趨勢逐漸減緩，4.0 a后，天然氣采集量增幅平緩。在時(shí)間一定的情況下，隨著體積壓裂段數(shù)的增加，天然氣采集量增大，其中體積壓裂段數(shù)為15段時(shí)，天然氣采集量最大，5段時(shí)天然氣采集量最小。

2.2.2 水平井壓裂段長度優(yōu)化效果

不同水平井壓裂段長度下天然氣采集量隨時(shí)間的變化關(guān)系如圖4所示。由圖4可知，隨著開采時(shí)間的增加，天然氣采集量均呈現(xiàn)出持續(xù)上升的趨勢，說明無論使用何種壓裂段長，均能提升致密砂巖氣藏的開采量。4種段長相比而言，當(dāng)段長為750 m時(shí)天然氣開采量最大，因此筆者選擇750 m段長作為調(diào)控體積壓裂Y341封隔器裂縫控制的參數(shù)。

2.2.3 裂縫夾角優(yōu)化效果

不同裂縫夾角下天然氣采集量隨時(shí)間的變化關(guān)系如圖5所示。由圖5可知，當(dāng)裂縫夾角呈垂直狀態(tài)（90°）時(shí)天然氣采集量最大，而在裂縫夾角較小的情況下，研究區(qū)域的天然氣采集量明顯較低，說明較小的裂縫夾角不利于改善致密砂巖氣藏的開采，因此為實(shí)現(xiàn)水平井體積壓裂裂縫控制，宜選取裂縫夾角參數(shù)為90°。

2.2.4 裂縫半長優(yōu)化效果

不同裂縫半長（125、195、295、335 m）下天然氣采集量隨時(shí)間的變化關(guān)系如圖6所示?？梢钥闯?，當(dāng)裂縫半長較小（125 m）時(shí)，天然氣采集量的增長趨勢相對平緩，表明采用較小的裂縫半長，即使配合體積壓裂技術(shù)，也無法有效提高天然氣的采集量。當(dāng)裂縫半長為335 m時(shí)，天然氣采集量達(dá)到最大，說明可以通過Y341封隔器擴(kuò)大裂縫半長，從而顯著提高研究區(qū)域的天然氣采集量。因此，只要實(shí)際環(huán)境滿足要求，可以通過Y341封隔器擴(kuò)大裂縫半長來優(yōu)化體積壓裂技術(shù)中的裂縫參數(shù)，從而提高天然氣的采集量。

2.2.5 裂縫密度優(yōu)化效果

設(shè)定每1 000 m水平井分別存在10、15、20條裂縫（裂縫密度），得到不同裂縫密度下天然氣采集量隨時(shí)間的變化關(guān)系如圖7所示?？梢钥闯?，隨著裂縫密度的增加，天然氣采集量呈上升趨勢，說明在裂縫密度較大的情況下，能夠提升致密砂巖氣藏水平井的天然氣采集量。因此，在Y341封隔器裂縫控制中應(yīng)提高裂縫密度。

2.2.6 裂縫導(dǎo)流能力優(yōu)化效果

設(shè)定裂縫導(dǎo)流能力分別為15、25、35 D×m，不同裂縫導(dǎo)流能力下天然氣采集量隨時(shí)間的變化關(guān)系如圖8所示。由圖8可知，裂縫導(dǎo)流能力越大，研究區(qū)域的天然氣采集量越大，當(dāng)裂縫導(dǎo)流能力為35 D×m時(shí)，4.5 a后天然氣采集量幾乎突破1.1×104 t，與開采初期相比開采量的增長率得到了極大提升。因此，在調(diào)控致密砂巖氣藏水平井體積壓裂裂縫時(shí)，選擇35 D×m的裂縫導(dǎo)流能力能夠顯著提升研究區(qū)域的裂縫性能，最終達(dá)到提升天然氣采集量的目的。

3 結(jié)束語

筆者采用模擬軟件構(gòu)建研究區(qū)域體積壓裂范圍的區(qū)域模型，利用該模型分析Y341封隔器裂縫路徑控制中，不同裂縫參數(shù)變化下的天然氣采集量的變化情況，并根據(jù)分析結(jié)果確定體積壓裂Y341封隔器裂縫路徑控制中最終的參數(shù)優(yōu)化值。將模擬計(jì)算確定的裂縫參數(shù)應(yīng)用到實(shí)際的致密砂巖氣藏水平井體積壓裂施工中，可以綜合性提升研究區(qū)域的天然氣產(chǎn)量。

參 考 文 獻(xiàn)

［1］" "黃卓.“工廠化”壓裂多裂縫應(yīng)力干擾與延伸規(guī)律研究［J］.斷塊油氣田，2022，29（4）：572-576.

［2］" "付海峰，才博，庚勐，等.基于儲(chǔ)層縱向非均質(zhì)性的水力壓裂裂縫三維擴(kuò)展模擬［J］.天然氣工業(yè)，2022，42（5）：56-68.

［3］" "宋毅，林然，黃浩勇，等.深層頁巖氣水平井壓裂異步起裂裂縫延伸模擬與調(diào)控［J］.大慶石油地質(zhì)與開發(fā)，2022，41（5）：145-152.

［4］" "胡東風(fēng)，任嵐，李真祥，等.深層超深層頁巖氣水平井縫口暫堵壓裂的裂縫調(diào)控模擬［J］.天然氣工業(yè)，2022，42（2）：50-58.

［5］" "向洪，隋陽，王靜，等.勝北深層致密砂巖氣藏水平井細(xì)分切割體積壓裂技術(shù)［J］.石油鉆采工藝，2021，43（3）：368-373.

［6］" "肖陽，劉守昱，何永志，等.致密砂巖裂縫性氣藏縫網(wǎng)壓裂裂縫復(fù)雜程度評價(jià)方法［J］.特種油氣藏，2022，

29（2）：157-163.

［7］" "費(fèi)世祥，余浩杰，陳存良，等.致密砂巖氣藏水平井開發(fā)關(guān)鍵技術(shù)——以蘇里格氣田上古生界為例［J］.西安石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2022，37（4）：26-35.

［8］" "張慶輝，高發(fā)潤，譚武林，等.柴達(dá)木盆地英西地區(qū)水平井體積改造及縫網(wǎng)預(yù)測技術(shù)研究［J］.石油物探，2022，61（6）：1043-1052；1089.

［9］" " 胡納川，何東博，郭建林，等.致密砂巖氣藏水平井整體開發(fā)提高采收率——以長嶺氣田登婁庫組氣藏為例［J］.西安石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2021，36（5）：55-63.

［10］" "張衍君，曾會(huì)，陶秀娟，等.致密儲(chǔ)層體積壓裂裂縫漏失及控制綜述［J］.科學(xué)技術(shù)與工程，2022，22（26）：11277-11286.

［11］" "余佩蓉，鄭國慶，孫福泰，等.瑪湖凹陷風(fēng)城組頁巖油藏水平井壓裂裂縫擴(kuò)展模擬［J］.新疆石油地質(zhì)，2022，43（6）：750-756.

［12］" "劉傳喜，方文超，秦學(xué)杰.非常規(guī)油氣藏壓裂水平井動(dòng)態(tài)縫網(wǎng)模擬方法及應(yīng)用［J］.石油與天然氣地質(zhì)，2022，43（3）：696-702.

［13］" "王華，崔越華，劉雪玲，等.致密砂巖氣藏多層系水平井立體開發(fā)技術(shù)——以鄂爾多斯盆地致密氣示范區(qū)為例［J］.天然氣地球科學(xué)，2021，32（4）：472-480.

［14］" "康毅力，邵佳新，游利軍，等.深層致密砂巖氣藏干濕交替誘發(fā)氣井出砂實(shí)驗(yàn)?zāi)M［J］.油氣地質(zhì)與采收率，2021，28（2）：127-134.

［15］" "慕立俊，吳順林，徐創(chuàng)朝，等.基于縫網(wǎng)擴(kuò)展模擬的致密儲(chǔ)層體積壓裂水平井產(chǎn)能貢獻(xiàn)分析［J］.特種油氣藏，2021，28（2）：126-132.

［16］" "楊林，劉彧軒，向斌，等.金華—秋林致密氣藏直井層內(nèi)分簇壓裂裂縫擴(kuò)展規(guī)律及應(yīng)用［J］.鉆采工藝，2021，44（2）：49-51；57.

（收稿日期：2023-10-12，修回日期：2024-09-18）

Simulation of Y341 Packer Fracture Control in Volume Fracturing

of Horizontal Wells in Sandstone Gas Reservoirs

FAN Yong?tao， LIANG Dong?yang， ZHANG Zhao?hong， HUANG Xin?yu， WANG Xin

（Oil amp; Gas Exploration Company， Shaanxi Yanchang Petroleum （Group） Co.， Ltd.）

Abstract" "Making use of Eclipse software establish the crack control model of Y341 packer in volume fracturing was implemented. In which， the different stage number， stage length， fracture angle， fracture half? length， fracture density and fracture diversion ability were taken as the control factors to constrain control ability of Y341 packer to fractures. Simulation results show that， the Y341 packer can complete fracture control and improve gas collection in horizontal wells of tight sandstone gas reservoirs under the control of constraint parameters.

Key words" " Y341 packer， dense sandstone， gas reservoir horizontal wells， volume fracturing， crack control simulation