非常規(guī)油氣藏體積壓裂數(shù)值模擬新進(jìn)展

唐子春,王 朝,張子珂,胡 蝶,雷征東

(1.中國石油勘探開發(fā)研究院，北京 100083；2.中國石油大學(xué)(北京))

非常規(guī)油氣藏體積壓裂數(shù)值模擬新進(jìn)展

唐子春1,2,王 朝1,2,張子珂1,2,胡 蝶1,2,雷征東1

(1.中國石油勘探開發(fā)研究院，北京 100083；2.中國石油大學(xué)(北京))

頁巖氣藏的高效開發(fā)與體積壓裂的合理設(shè)計(jì)密切相關(guān)，在調(diào)研大量相關(guān)文獻(xiàn)的基礎(chǔ)上，綜述了近年來國外非常規(guī)油氣藏水平井體積壓裂數(shù)值模擬技術(shù)。從天然裂縫與人工縫的相互作用、應(yīng)力影以及支撐劑模擬等方面敘述了壓裂力學(xué)模型的發(fā)展歷程，介紹了非常規(guī)壓裂模型的原理，為該技術(shù)的深入研究提供參考;闡述了連續(xù)介質(zhì)模型和離散裂縫模型在非常規(guī)油氣藏體積壓裂數(shù)值模擬中的應(yīng)用，分析了其原理、數(shù)模思路以及適用性;總結(jié)了國外研究的主要方向，提出了針對長慶油田致密儲(chǔ)層體積壓裂數(shù)值模擬的發(fā)展方向。

長慶油田；非常規(guī)油氣藏；數(shù)值模擬；體積壓裂；離散裂縫模型

北美頁巖氣的勘探與開發(fā)至今已有百年的歷史，而真正進(jìn)入商業(yè)化生產(chǎn)僅有20多年。水平井開發(fā)、體積壓裂等技術(shù)的突破促進(jìn)了美國的頁巖氣革命，使非常規(guī)油氣在世界能源供給中占有越來越大的比例。從非常規(guī)油氣藏獲得工業(yè)性的經(jīng)濟(jì)產(chǎn)量有賴于水平井近井地帶的改造體積，而在體積壓裂過程中產(chǎn)生的水力裂縫會(huì)與天然裂縫相互交錯(cuò)，在水平井近井地帶形成復(fù)雜的裂縫網(wǎng)絡(luò)。因此，如何合理、準(zhǔn)確地描述復(fù)雜縫網(wǎng)系統(tǒng)成為非常規(guī)油氣藏中數(shù)值模擬的關(guān)鍵。近年來，國外非常規(guī)油氣藏體積壓裂的數(shù)值模擬技術(shù)有了飛躍式的發(fā)展，而目前國內(nèi)還沒有對此最新研究成果進(jìn)行系統(tǒng)論述。本文針對該類問題，大量調(diào)研了近年來的相關(guān)研究，總結(jié)了雙重介質(zhì)模型、離散裂縫模型的數(shù)值模擬方法，提出了相應(yīng)的數(shù)模思路以及對未來技術(shù)的展望。

1 復(fù)雜縫網(wǎng)模型

不同于PKN模型、擬三維模型等常規(guī)壓裂模型，復(fù)雜縫網(wǎng)模型需要考慮很多因素，例如天然裂縫和水力裂縫的相互作用、應(yīng)力影的影響等。Olson等人提出了一個(gè)縫網(wǎng)模型來預(yù)測水力裂縫的擴(kuò)展，該模型雖然考慮了人工縫與天然縫的相互作用，但不能模擬壓裂液流動(dòng)和支撐劑輸送；Rogers等人[1]提出了一個(gè)離散裂縫模型，簡化了天然裂縫與人工縫的關(guān)系，但沒有考慮壓裂縫之間的作用；Nagel等人[2,3]提出了一個(gè)擬三維縫網(wǎng)模擬器，可模擬天然裂縫與人工縫之間的相互作用，但該模型計(jì)算量大，不具有工程適用性；而最為經(jīng)典的是Xu等人[4,5]提出的線網(wǎng)模型和Weng等人[6]提出的非常規(guī)壓裂模型。線網(wǎng)模型(Wire-Mesh)將復(fù)雜裂縫網(wǎng)絡(luò)簡化為一個(gè)橢球體，該橢球體由兩組正交、等距分布的垂直裂縫面組成。線網(wǎng)模型應(yīng)用了巖石力學(xué)理論，計(jì)算了縫網(wǎng)隨著壓裂施工的擴(kuò)展，考慮了天然裂縫與水力裂縫的力學(xué)作用，模擬了支撐劑在縫網(wǎng)中的輸送。非常規(guī)壓裂模型(UFM)是Kresse.O[7]于2011年提出的新模型，模型主要考慮以下幾個(gè)方面：

(1)基本的控制方程。①流動(dòng)方程：一維泊肅葉圓管層流，假定流體為冪律流體；②質(zhì)量守恒；③考慮液體濾失效應(yīng)：限制條件是射孔點(diǎn)流速之和等于總排量。

(2)高度生長模型。主要基于擬三維壓裂模型計(jì)算，考慮流體壓力、地應(yīng)力、斷裂韌性、彈性模量等因素，在裂縫尖端的上部、下部分別計(jì)算斷裂韌性，判斷裂縫高度延伸的可能方位。

(3)裂縫相互作用模型[8]。該方法是最早實(shí)現(xiàn)水力裂縫與天然裂縫相互作用的模型之一，模型本身基于線彈性斷裂力學(xué)，沒有考慮液體影響。水力裂縫默認(rèn)垂直于最小水平主應(yīng)力方向起裂，單一裂縫時(shí)，非常規(guī)壓裂模型退化為擬三維壓裂模型。Warpinski[9]提出，當(dāng)水力裂縫遇到天然裂縫時(shí)，會(huì)出現(xiàn)穿過和滑移兩種情況，而Blaton[10]的觀點(diǎn)是會(huì)出現(xiàn)張開、穿過、捕捉三種情況，但二者本質(zhì)相同，分別由圖1b和c所示。

圖1 水力裂縫遇到天然裂縫的可能情況

(4)應(yīng)力影模型。模型使用邊界元(BEM)方法，該方法考慮了三維應(yīng)力修正項(xiàng)，模擬結(jié)果與CSIRO十分接近[11]。式(1)和式(2)分別代表第i個(gè)裂縫單元在受到應(yīng)力影作用后的法向與切向應(yīng)力值。

(1)

(2)

式中n，s分別代表張開、剪切應(yīng)力；Dn，Ds代表張開、剪切位移；N代表裂縫體內(nèi)總單元數(shù)；Cij是二維平面應(yīng)變影響因子，該因子是用于作用于計(jì)算來自相鄰裂縫單元的額外影響應(yīng)力值，具體表達(dá)式參見Crouch[12](1976)；Aij是三維修正系數(shù)，修正因子Cij。當(dāng)距離增加時(shí)，三維條件下有限的裂縫高度會(huì)削弱裂縫單元間的影響，需要引入該系數(shù)進(jìn)行修正。具體表達(dá)式參見Olson[13](2004)。

非常規(guī)壓裂模型在考慮應(yīng)力干擾時(shí)有兩大假定：①裂縫寬度等于法向位移不連續(xù)度；②裂縫表面剪切應(yīng)力等于0?；谶@兩個(gè)假設(shè)才能計(jì)算切向位移不連續(xù)度及法向的應(yīng)力干擾值。應(yīng)力影模型對裂縫傳遞過程的影響體現(xiàn)在兩個(gè)方面：一是影響井場所有裂縫單元的初始地應(yīng)力大小與方向，改變了裂縫的壓力與寬度；二是影響了裂縫尖端的應(yīng)力場，可能導(dǎo)致裂縫沿其他方向擴(kuò)展。多條裂縫的相互干擾模擬結(jié)果：當(dāng)裂縫間距低于10 m時(shí)，水力裂縫之間表示出強(qiáng)烈的排斥現(xiàn)象，而間距40 m以上排斥現(xiàn)象變得比較微弱。

(5)支撐劑傳輸模型。為提高計(jì)算效率，水平方向(沿縫長)設(shè)定為一維模型，垂向(沿縫高)設(shè)定為三層模型，從底部到上部分別為砂堤、混砂液、無砂液，同時(shí)考慮支撐劑的沉降與侵蝕作用。Cohen等人[14]對非常規(guī)壓裂模型的各個(gè)參量進(jìn)行了1800多次的擬合試驗(yàn)，結(jié)果表明，支撐劑尺寸為40～70目時(shí)，可流動(dòng)面積比30～50目大了3倍以上，分布也更為均勻。模型計(jì)算過程考慮了各個(gè)模型的相互作用，對每個(gè)時(shí)步，該模型基本計(jì)算流程如圖2所示。

圖2 UFM模型計(jì)算流程

非常規(guī)壓裂模型最大的優(yōu)勢在于，除了考慮常規(guī)壓裂軟件的基本功能外，還考慮了多縫間應(yīng)力干擾，包括單井多段間應(yīng)力干擾與井間應(yīng)力干擾，在現(xiàn)有的壓裂模型中有獨(dú)特的優(yōu)勢。該模型還可以模擬與天然裂縫的相互作用，從而進(jìn)行多井壓裂等復(fù)雜的壓裂設(shè)計(jì)，為壓裂領(lǐng)域較為前沿的壓裂設(shè)計(jì)模型。但是，模型本身的天然裂縫與水力裂縫的擴(kuò)展準(zhǔn)則是基于線彈性斷裂力學(xué)理論，沒有考慮流體的影響，可能會(huì)造成較大的偏差；地層天然裂縫識(shí)別困難，給建立天然裂縫網(wǎng)絡(luò)帶來一定的不確定性。目前，模型還處于進(jìn)一步完善之中。

2 體積壓裂數(shù)值模擬

2.1 連續(xù)介質(zhì)模型

連續(xù)介質(zhì)模型包括雙孔模型、雙孔雙滲模型以及多重孔隙介質(zhì)模型。傳統(tǒng)的雙孔模型假設(shè)基質(zhì)到裂縫的流動(dòng)為擬穩(wěn)定流，而在非常規(guī)油氣藏中，流動(dòng)并不是瞬間發(fā)生，需要多級基質(zhì)網(wǎng)格來模擬致密儲(chǔ)層中的微孔隙、裂縫和有機(jī)質(zhì)區(qū)域，獲取過渡流特征。Pruess[15]提出的MINC模型解決了這個(gè)問題，該模型將基質(zhì)單元剖分成一系列次級單元，這些次級單元形成獨(dú)立的一維流動(dòng)系統(tǒng)，并與對應(yīng)的裂縫單元相連，后續(xù)研究又將模型中次級單元與裂縫表面的距離改進(jìn)為對數(shù)變化。Zhang[16]使用ECLIPSE中的多重孔隙結(jié)構(gòu)來描述體積壓裂區(qū)域基質(zhì)和裂縫的關(guān)系，提出了一個(gè)綜合模擬方法，研究油藏和水力壓裂參數(shù)對頁巖氣藏生產(chǎn)動(dòng)態(tài)的影響，模擬結(jié)果表明，該方法很好地獲取了基質(zhì)到裂縫之間的過渡流動(dòng)特征。

對于水力裂縫的模擬，Kalantari Dahaghi[17]研究認(rèn)為，對數(shù)間距局部網(wǎng)格加密對水平井體積壓裂區(qū)域有較好的模擬效果。在此基礎(chǔ)上，Cipolla和CMG公司的Rubin[18-19]基于雙重介質(zhì)模型提出了一種叫做“雙滲模型-對數(shù)間距變化-局部網(wǎng)格加密(DK-LS-LGR)”的方法，結(jié)合微地震監(jiān)測結(jié)果，對致密儲(chǔ)層中的復(fù)雜縫網(wǎng)進(jìn)行整體建模。該方法在體積改造區(qū)域內(nèi)部采用對數(shù)間距變化、局部網(wǎng)格加密的雙滲模型，在外部采用普通的雙滲網(wǎng)格。復(fù)雜縫網(wǎng)的建模首先通過微地震監(jiān)測數(shù)據(jù)篩選體積改造區(qū)域的裂縫網(wǎng)絡(luò)，之后，通過對數(shù)間距局部網(wǎng)格加密模擬水平井的體積改造區(qū)域。結(jié)果表明，該方法歷史擬合效果好，可以較好反映頁巖的生產(chǎn)特征。Du等人[20-21]提出了一個(gè)綜合的工作流對Barnett頁巖氣藏進(jìn)行模擬。在地質(zhì)模型的建立過程中，天然裂縫建模基于井筒成像、測井和微地震監(jiān)測的解釋結(jié)果，水力裂縫建模依賴于壓裂施工參數(shù)、支撐劑的分布以及微地震監(jiān)測的結(jié)果，采用Oda方法將油藏地質(zhì)模型和離散裂縫模型粗化為雙孔模型?，F(xiàn)場應(yīng)用結(jié)果表明，這個(gè)雙孔模型工作流可以較好地模擬非常規(guī)油氣藏的生產(chǎn)，使人們對致密儲(chǔ)層的生產(chǎn)有了更深的認(rèn)識(shí)，同時(shí)，也為非常規(guī)油氣藏一體化模擬工作流的快速發(fā)展奠定了基礎(chǔ)。

2.2 離散裂縫模型

離散裂縫模型(DFM)是指用一系列指定網(wǎng)格代表裂縫，使裂縫離散分布在模型中，利用網(wǎng)格之間的傳導(dǎo)率控制基質(zhì)與裂縫之間的流動(dòng)。根據(jù)求解方法的不同，可將該方法分為非結(jié)構(gòu)化離散裂縫模型(UDFM)和嵌入式離散列分模型(EDFM)[22-23]。

2.2.1 非結(jié)構(gòu)化離散裂縫模型

常規(guī)的離散裂縫模型有賴于非結(jié)構(gòu)化網(wǎng)格以及局部網(wǎng)格加密實(shí)現(xiàn)復(fù)雜縫網(wǎng)的精細(xì)描述，Karimi-Fard[24]、Fu[25]、Marcondes[26]等人應(yīng)用有限元方法建立了離散裂縫模型，采用非結(jié)構(gòu)化網(wǎng)格對裂縫進(jìn)行了顯性剖分。

近年來，通過將地震、測井、微地震監(jiān)測以及巖心分析結(jié)合到復(fù)雜縫網(wǎng)模型，進(jìn)行不同方式的網(wǎng)格化，實(shí)現(xiàn)了離散裂縫模型到油藏?cái)?shù)值模型的耦合。Cipolla等人[27-29]提出了地震-數(shù)模一體化工作流,該工作流包括以下步驟：①通過地震、測井解釋以及巖心分析手段，提供必要的巖石力學(xué)性質(zhì)、地應(yīng)力分布、斷層以及天然裂縫的分布，建立離散裂縫模型和地質(zhì)力學(xué)模型，并運(yùn)用微地震監(jiān)測進(jìn)行約束；②基于已建立的地質(zhì)模型，使用非常規(guī)壓裂模型(UFM)輸入壓裂施工參數(shù)和巖石力學(xué)性質(zhì)，通過計(jì)算耦合人工裂縫和天然裂縫，生成復(fù)雜縫網(wǎng)系統(tǒng)；③通過自動(dòng)非結(jié)構(gòu)網(wǎng)格化算法進(jìn)行局部網(wǎng)格加密，生成油藏?cái)?shù)值模型，由下一代數(shù)值模擬器INTERSECT進(jìn)行計(jì)算。圖3為頁巖氣藏水平井衰竭式開發(fā)10年以后壓力分布模擬結(jié)果，從圖中可以看出，該工作流可以精細(xì)地描述水平井體積壓裂區(qū)域的開發(fā)特征。

圖3 體積壓裂水平井衰竭開采10年后的壓力分布

在此基礎(chǔ)上，Cipolla等人[30]提出了一個(gè)完井設(shè)計(jì)工具，將其耦合到了工作流中。該工具利用應(yīng)力場、巖石力學(xué)特性和成像技術(shù)，對頁巖孔隙度、黏土含量、有機(jī)碳含量(TOC)等參數(shù)進(jìn)行篩選，確定合理的壓裂段和射孔位置，進(jìn)行嚴(yán)格高效的水平井分段壓裂完井設(shè)計(jì)，判別標(biāo)準(zhǔn)如表1所示。Cohen等人[31]也將工作流做了一些改進(jìn)，包括壓裂段之間的應(yīng)力干擾、壓裂液溫度和流變性的耦合等因素。目前，該流程已在斯倫貝謝公司的Mangrove軟件中實(shí)現(xiàn)，工作流如圖4所示。

地震-數(shù)模一體化工作流實(shí)現(xiàn)了壓裂模型與油藏?cái)?shù)值模擬器之間的耦合，但是沒有有效的工具來模擬體積壓裂對非常規(guī)油氣藏單井產(chǎn)能的影響。Cohen等人[32]耦合了一個(gè)半解析的非常規(guī)產(chǎn)能模型(UPM)，進(jìn)行多參數(shù)敏感性分析，研究壓裂施工參數(shù)設(shè)計(jì)和油藏性質(zhì)與產(chǎn)能之間的關(guān)系。在此基礎(chǔ)上，Cohen等人[33]耦合了一個(gè)溫度模型，研究壓裂液體系和油藏溫度對產(chǎn)能的影響。

圖4 地震-數(shù)模一體化工作流

油藏性質(zhì)完井性質(zhì)總有機(jī)碳>=3%低應(yīng)力值天然氣含量>2.83m3/t電阻率>15Ω·m干酪根含量高黏土含量<40%泥巖孔隙度>4%楊氏模量>13.8GPa泊松比<0.2氣測滲透率>10-4μm2中子測井孔隙度<35%密度孔隙度>8%

2.2.2 嵌入式離散裂縫模型

雙重介質(zhì)模型建立的縫網(wǎng)形態(tài)相對理想，但不能正確描述改造區(qū)域的復(fù)雜情況；而常規(guī)離散裂縫模型依賴于局部網(wǎng)格加密或非結(jié)構(gòu)化網(wǎng)格，可以滿足高精度的要求，但網(wǎng)格數(shù)目增加，造成計(jì)算量呈數(shù)量級增長。為了解決這些問題，Li和Lee(2008)[34]首次提出了嵌入式離散裂縫模型(EDFM)，該模型使用結(jié)構(gòu)化網(wǎng)格描述基質(zhì)系統(tǒng)，通過計(jì)算裂縫與基質(zhì)網(wǎng)格的交會(huì)部分，顯性地引入裂縫系統(tǒng)控制體作為單獨(dú)的一部分，將裂縫平面嵌入到基質(zhì)網(wǎng)格中并進(jìn)行邊界元離散化。

Moinfar[35-36]系統(tǒng)地發(fā)展了嵌入式離散裂縫模型，增加了傾斜裂縫的特性，并將其應(yīng)用于全隱式組分模擬器中。由于兩個(gè)系統(tǒng)的物質(zhì)平衡方程之間沒有直接關(guān)聯(lián)，需要通過非相鄰連接(NNC)實(shí)現(xiàn)相互之間的流體流動(dòng)，即模型中需要計(jì)算的四種傳導(dǎo)率——裂縫系統(tǒng)中裂縫控制體之間的傳導(dǎo)率、基質(zhì)網(wǎng)格與裂縫控制體之間的傳導(dǎo)率、相交裂縫之間的傳導(dǎo)率以及裂縫和井筒之間的傳導(dǎo)率。通過定義多種類型的傳導(dǎo)率，嵌入式離散裂縫模型可以動(dòng)態(tài)地表征裂縫特征，在每個(gè)時(shí)間步末端更新裂縫的開度和滲透率，重新計(jì)算傳導(dǎo)率并應(yīng)用到下一個(gè)時(shí)間步的流動(dòng)方程中。Shakiba等人[37-38]利用微地震監(jiān)測對嵌入式離散裂縫模型進(jìn)行約束，同時(shí)應(yīng)用組分模擬器UTCOMP更加精準(zhǔn)而高效地對體積壓裂區(qū)域的縫網(wǎng)系統(tǒng)進(jìn)行了數(shù)值模擬。

嵌入式離散裂縫模型繼承了雙重介質(zhì)模型的優(yōu)勢，不需要裂縫控制體與基質(zhì)網(wǎng)格相鄰，也不要求基質(zhì)和裂縫系統(tǒng)具有相同的網(wǎng)格尺寸。同時(shí)，該模型離散化表征了裂縫系統(tǒng)，不需要在模型中使用局部網(wǎng)格加密或非結(jié)構(gòu)網(wǎng)格化。上述優(yōu)勢使該模型可以在常規(guī)數(shù)值模擬器中應(yīng)用，大大提高了計(jì)算效率和精確度。

3 發(fā)展趨勢

國外頁巖氣藏開發(fā)具有獨(dú)特的商業(yè)化模式，需要在衰竭式開采條件下將采收率最大化，其研究內(nèi)容主要包括：①體積壓裂工藝參數(shù)對頁巖氣產(chǎn)能的影響，例如排量、支撐劑類型、壓裂液體系和溫度等；②完井工藝對頁巖氣藏生產(chǎn)動(dòng)態(tài)的影響，例如分段多簇壓裂中的壓裂段數(shù)、段間距以及每段中的射孔簇?cái)?shù)等因素，相應(yīng)的數(shù)值模擬技術(shù)的突破與創(chuàng)新也著眼于這些方面。Mangrove的一體化工作流主要為單井模擬，而近年來對頁巖儲(chǔ)層中井間干擾現(xiàn)象的關(guān)注度越來越高；同時(shí)，隨著拉鏈?zhǔn)綁毫押途S作業(yè)的推廣，基于多井模型的非常規(guī)油藏?cái)?shù)值模擬技術(shù)開始在國外發(fā)展。

不同于非常規(guī)頁巖儲(chǔ)層，長慶油田主要是致密砂巖儲(chǔ)層，在壓裂施工和射孔完井工藝上不需要進(jìn)行有機(jī)碳含量、黏土含量等因素的篩選，即"甜點(diǎn)區(qū)"的優(yōu)選；同時(shí)，長慶油田主要采用注水開發(fā)致密油，在開發(fā)方式上與國外存在根本性的差異。因此，體積壓裂數(shù)值模擬技術(shù)的發(fā)展，需要考慮超低滲、致密儲(chǔ)層中的應(yīng)力敏感、啟動(dòng)壓力梯度等特殊現(xiàn)象，并將研究重點(diǎn)放在注采井網(wǎng)的優(yōu)化上來。

4 結(jié)論

(1)復(fù)雜縫網(wǎng)模型逐漸向復(fù)雜化、系統(tǒng)化的方向發(fā)展。非常規(guī)壓裂模型可以較好地模擬壓裂過程中縫網(wǎng)的形成及縫間應(yīng)力干擾等問題，能夠并入油藏?cái)?shù)值模型進(jìn)行一體化分析，對現(xiàn)場問題提供了很好的解決方案。過于簡化的天然裂縫與水力裂縫的擴(kuò)展準(zhǔn)則，使裂縫擴(kuò)展與實(shí)際有一定偏差，且對天然裂縫的描述仍然停留在統(tǒng)計(jì)意義上，很難準(zhǔn)確描述地層裂縫系統(tǒng)，模擬的裂縫形態(tài)有很大不確定性。在以后的研究中，仍然需要大量的現(xiàn)場數(shù)據(jù)對非常規(guī)壓裂模型進(jìn)行完善與修改，使之更加精確與高效。

(2)對于體積壓裂水平井的數(shù)值模擬，首先需要地質(zhì)建模，通過力學(xué)模型與微地震手段的結(jié)合對復(fù)雜縫網(wǎng)進(jìn)行描述；對于流動(dòng)模型的選擇，如果選用雙重介質(zhì)模型，則需要一種合理的粗化技術(shù)；如果選擇離散裂縫模型，則需要一種高效的非結(jié)構(gòu)化網(wǎng)格技術(shù)以及局部網(wǎng)格加密方法對裂縫網(wǎng)絡(luò)準(zhǔn)確刻畫。

(3)雙重介質(zhì)模型假設(shè)的縫網(wǎng)形態(tài)較為均勻，不能準(zhǔn)確描述實(shí)際情況；常規(guī)離散裂縫模型依賴于非結(jié)構(gòu)網(wǎng)格化技術(shù)，具有較高的精確度，但是大量的網(wǎng)格數(shù)目造成了巨大的計(jì)算量；嵌入式離散裂縫模型綜合了兩種模型的優(yōu)勢，較好地解決了一些復(fù)雜的問題。近年來發(fā)展較快的一體化工作流，基于連續(xù)介質(zhì)模型和離散裂縫模型，綜合應(yīng)用地質(zhì)研究、微地震監(jiān)測和壓裂參數(shù)等手段對體積壓裂水平井中的復(fù)雜縫網(wǎng)進(jìn)行了表征，表現(xiàn)出良好的效果。

[1] STEPHEN R, DAVIDE E, RORY D, et al. Understanding hydraulic fracture geometry and interactions in the horn river basin through DFN and numerical modeling[C].SPE 137488, 2010.

[2] NEAL B N, IVAN G, MARISELA S,et al.Simulating hydraulic fracturing in real fractured rocks-overcoming the limits of pseudo 3D models[C].SPE 140480, 2011.

[3] NEAL B N, IVAN G, MARISELA S, et al. Stress shadowing and microseismic events: A Numerical Evaluation[C].SPE 147363, 2011.

[4]WENYUE X,MARC J, UTPAL G,et al.Wiremesh:A novel shale fracturing simulator[C].SPE 132218, 2010.

[5] WENYUE X, JAMES L, MIKE D. Quick estimate of initial production from stimulated reservoirs with complex hydraulic fracture network[C].SPE 146753, 2011.

[6]XIAOWEI W, OLGA K, CHARLES E C, et al. Modeling of hydraulic fracture network propagation in a naturally fractured formation[C].SPE 140253, 2011.

[7] KRESSE O, COHEN C, WENG X, et al. Numerical modeling of hydraulic fracturing in naturally fractured formations[J].Rock Mechanics,2011,15(5)516-535.

[8] GU H, WENG X, LUND J B, et al. Hydraulic fracture crossing natural fracture at non-orthogonal angles, a criterion, its validation and applications[C].SPE 139984, 2011.

[9] WARPINSKI N R, TEUFEL LW. Influence of geologic discontinuities on hydraulic fracture propagation[C].SPE 13224,1987.

[10] BLANTON T L. An experimental study of interaction between hydraulically induced and pre-existing fractures[C].SPE 10847, 1982.

[11]KRESSE O, WENG X, WU R, et al. Numerical modeling of hydraulic fractures interaction in complex naturally fractured formations[J]. Rock Mechanics and Rock Engieering,2013,46(3):555-568.

[12] CROUCH S L. Solution of plane elasticity problems by the displacement discontinuity method[J]. International Journal for Numerical Methods in Engineering, 1976, 10(2):301-343.

[13] OISON J E. Predicting fracture swarms——the influence of subcritical crack growth and the crack-tip process zone on joint spacing in rock[J]. Geological Society London Special Publications, 2004, 231(1):73-88.

[14] COHEN C E, ABAD C, WENG X, et al. Analysis on the impact of fracturing treatment design and re-servoir properties on production from shale gas reservoirs[C].IPTC 16400, 2013.

[15] KARSTEN P. A practical method for modeling fluid and heat flow in fractured porous media[C].SPE 10509, 1985.

[16] XU Z, CHANGAN D, FRANZ D, et al. Sensitivity studies of horizontal wells with hydraulic fractures in shale gas reservoirs[C]. IPTC 13338, 2009.

[17]AMIRMASOUD K D. Numerical simulation and modeling of enhanced gas recovery and CO2sequestration in shale gas reservoirs: A Feasibility Study[C]. SPE 139701, 2010.

[18] CRAIG L C, ELYEZER L, JAMES C E, et al. Re-servoir modeling in shale-gas reservoirs[C].SPE 125530, 2009.

[19] BARRY R. Accurate simulation of non-darcy flow in stimulated fractured shale reservoirs[C].SPE 132093, 2010.

[20] CHANGAN D, XU Z, BRAD M, et al. A workflow for integrated Barnett shale gas reservoir modeling and simulation[C].SPE 122934, 2009.

[21] CHANGAN M D, XU Z, LANG Z, et al. Modeling hydraulic fracturing induced fracture networks in shale gas reservoirs as a dual porosity system[C].SPE 132180, 2010.

[22] HUI M H， MALLISON B T. System and method for predicting fluid flow characteristics within fractured subsurface reservoirs:US[P],2010.

[23] HUI M H， MALLISON B T. System and method of grid generation for discrete fracture modeling: US[P],2011.

[24] MOHAMMAD K F, ABBAS F. Numerical simulation of water injection in fractured media using the discrete-fracture model and the Galerkin method[C].SPE 83633, 2003.

[25] FU Y, YANG Y K, DEO M. Three-dimensional, three-phase discrete-fracture reservoir simulator based on control volume finite element formulation[C].SPE 93292, 2005.

[26] MARCONDES F. An element-based finite-volume method approach for naturally fractured compositional reservoir simulation.13th Brazilian Thermal Sciences Meeting, Uberlandia, Brazil, December.

[27] CRAIG L C, XIAOWEI W, MARK G M, et al. Integrating microseismic mapping and complex fracture modeling to characterize hydraulic fracture complexity[C]. SPE 140185, 2011.

[28] CRAIG L C, TONY F, MICHAEL J W, et al. Seismic-to-simulation for unconventional reservoir development[C].SPE 146876, 2011.

[29] MAHSANAM M, CRAIG L C. A workflow for modeling and simulation of hydraulic fractures in unconventional gas reservoirs[C].SPE 153022, 2012.

[30] CRAIG L C, XIAOWEI W, HITOSHI O, et al. New algorithms and integrated workflow for tight gas and shale completions[C].SPE 146872, 2011.

[31]HASHEMI R, CAOTHIEN S, ROVIG J. Field measurement of viscosified fluids[C]. SPE 17462, 1988.

[32] COHEN C E, WENYUE X, XIAOWEI W, et al. Production forecast after hydraulic fracturing in naturally fractured reservoirs: coupling a complex fracturing simulator and a semi-analytical production model[C].SPE 152541, 2012.

[33] COHEN C E, ABAD C, WENG X, et al. Analysis on the impact of fracturing treatment design and re-servoir properties on production from shale gas reservoirs[C].IPTC 16400, 2013.

[34] LIYONG L, SEONG H L. Efficient field-scale si-mulation for black oil in a naturally fractured reservoir via discrete fracture networks and homogenized media[C].SPE 103901, 2006.

[35] ALI M, ABDOLJALIL V, KAMY S, et al. Deve-lopment of a novel and computationally-efficient discrete-fracture model to study IOR processes in na-turally fractured reservoirs[C].SPE 154246, 2012.

[36] ALI M, KAMY S, RUSSELL T J, et al. Coupled geomechanics and flow simulation for an embedded discrete fracture model[C].SPE 163666, 2013.

[37] JOSESA C F, MAHMOOD S, ALI M, et al. Implementation of a preprocessor for embedded discrete fracture modeling in an IMPEC compositional reservoir simulator[C].SPE 173289, 2015.

[38] MAHMOOD S, KAMY S. Using embedded discrete fracture model and microseismic monitoring data to characterize the complex hydraulic fracture networks[C]. SPE 175142, 2015.

編輯：王金旗

1673-8217(2017)03-0108-06

2016-10-22

唐子春，1991年出生，2014年畢業(yè)于中國石油大學(xué)(北京)石油工程專業(yè)，現(xiàn)為中國石油大學(xué)(北京)石油與天然氣工程專業(yè)在讀碩士研究生，主要從事超低滲透油藏工程與數(shù)值模擬研究。

國家科技重大專項(xiàng)“低滲、特低滲油藏水驅(qū)擴(kuò)大波及體積方法與關(guān)鍵技術(shù)”(2017ZX05013-002)。

TE357.7

A