馮江榮 ，趙圣賢 ，夏自強(qiáng) ，李志宏 ，劉永旸 ，何沅翰 ，高攀 ，王高翔

（1.中國(guó)石油西南油氣田分公司頁(yè)巖氣研究院，四川 成都 618300；2.北京源烴泰克科技有限公司，北京 100102）

0 引言

在頁(yè)巖氣的勘探開(kāi)發(fā)過(guò)程中，數(shù)值模擬法已廣泛應(yīng)用于水平井水力壓裂模擬中，為水平井的設(shè)計(jì)、壓裂參數(shù)選擇和壓裂改造效果的最終評(píng)估提供了有效技術(shù)支撐。數(shù)值模擬法主要包括有限元法（FEM）、基于非連續(xù)性介質(zhì)力學(xué)的位移不連續(xù)方法（DDM）、基于連續(xù)性介質(zhì)力學(xué)的離散元法（DEM）等。其中，F(xiàn)EM是目前的主流方法［1-3］，但作為有網(wǎng)格的數(shù)值模擬方法，它無(wú)法消除裂尖應(yīng)力奇異性的作用，從而增大計(jì)算結(jié)果的誤差，在實(shí)際應(yīng)用中模擬應(yīng)力場(chǎng)精度無(wú)法滿足生產(chǎn)需要。

為了規(guī)避壓裂模擬過(guò)程中網(wǎng)格畸變?cè)斐傻挠?jì)算誤差，作為無(wú)網(wǎng)格數(shù)值模擬方法，物質(zhì)點(diǎn)法（MPM）越來(lái)越受到學(xué)者們關(guān)注。Aimene等［4］針對(duì)北美Marcellus和Eagle Ford頁(yè)巖氣儲(chǔ)層，利用MPM對(duì)水平井進(jìn)行了壓裂模擬，應(yīng)力場(chǎng)模擬與微地震監(jiān)測(cè)結(jié)果基本一致。Chip等［5］利用特拉華盆地公開(kāi)的斷層數(shù)據(jù)，采用MPM模擬得到了該盆地最大水平主應(yīng)力方向的變化規(guī)律，與實(shí)鉆井情況基本一致。Han等［6］針對(duì)新疆瑪湖油田百口泉組砂礫巖儲(chǔ)層，采用MPM模擬了不同壓裂段長(zhǎng)度情況下改造效果差異，優(yōu)選50 m段長(zhǎng)進(jìn)行完井壓裂方案設(shè)計(jì)及產(chǎn)量預(yù)測(cè)。

本文針對(duì)渝西地區(qū)五峰組—龍馬溪組一段 （簡(jiǎn)稱龍一段，下同）開(kāi)展了基于物質(zhì)點(diǎn)法的頁(yè)巖儲(chǔ)層壓裂模擬（簡(jiǎn)稱MPM壓裂模擬，下同），并與微地震監(jiān)測(cè)資料對(duì)比，驗(yàn)證了其適用性。在此基礎(chǔ)上提出了研究區(qū)水平井軌跡部署方案，應(yīng)用效果顯著。

1 MPM基本原理

MPM是Sulsky等［7］提出的一種數(shù)值模擬方法，它采用攜帶材料信息的物質(zhì)點(diǎn)離散材料區(qū)域，表征其運(yùn)動(dòng)和變形狀態(tài)；采用規(guī)則的歐拉背景計(jì)算空間導(dǎo)數(shù)和動(dòng)量方程，避免網(wǎng)格畸變，實(shí)現(xiàn)質(zhì)點(diǎn)間的相互作用與聯(lián)系。該方法非常適合處理瞬間施壓固體材料而導(dǎo)致其大幅變形的模擬問(wèn)題。

1.1 控制方程

物質(zhì)點(diǎn)法的控制方程包括質(zhì)量守恒以及動(dòng)量守恒方程：

式中：ρ為密度，g/cm3；t為時(shí)間，s；v 為速度，m/s；a 為重力加速度，m/s2；σ 為 Cauchy 應(yīng)力張量，Pa；b 為質(zhì)量力，m/s2。

在物質(zhì)點(diǎn)法的模擬中，空間連續(xù)分布的基質(zhì)要被離散成具有集中質(zhì)量的質(zhì)點(diǎn)。由于每個(gè)質(zhì)點(diǎn)的質(zhì)量在整個(gè)計(jì)算過(guò)程中保持不變，自動(dòng)滿足質(zhì)量守恒方程，質(zhì)點(diǎn)的質(zhì)量密度ρ( x，t)可表示為

式中：Np為質(zhì)點(diǎn)的總數(shù)，個(gè)；MP為質(zhì)點(diǎn)的質(zhì)量，kg；δ為Dirac δ函數(shù)；為質(zhì)點(diǎn)p在t時(shí)的坐標(biāo)，m。

為了確保求解方程過(guò)程中的動(dòng)量守恒，需要布置固定于空間的背景網(wǎng)格，利用節(jié)點(diǎn)基函數(shù)反映質(zhì)點(diǎn)p和背景網(wǎng)格節(jié)點(diǎn)i之間信息的映射關(guān)系，即：

滿足質(zhì)量守恒和動(dòng)量守恒的節(jié)點(diǎn)離散模擬方程為

1.2 工作流程

在物質(zhì)點(diǎn)法計(jì)算過(guò)程中，質(zhì)點(diǎn)攜帶物體的所有物質(zhì)信息，應(yīng)用在非常規(guī)油氣藏的壓裂模擬中，這些物質(zhì)信息包括彈性模量、泊松比、密度等儲(chǔ)層巖石物理參數(shù)以及斷層（裂縫）等構(gòu)造參數(shù)。背景網(wǎng)格節(jié)點(diǎn)根據(jù)質(zhì)點(diǎn)的空間分布自動(dòng)形成，但不包含任何物質(zhì)信息。在求解式（5）時(shí)，要將質(zhì)點(diǎn)在當(dāng)前時(shí)刻的質(zhì)量、動(dòng)量、應(yīng)力等信息，通過(guò)函數(shù)映射至背景網(wǎng)格節(jié)點(diǎn)進(jìn)行求解；之后，將網(wǎng)格節(jié)點(diǎn)的速度、位置增量映射回質(zhì)點(diǎn)，使得質(zhì)點(diǎn)速度、位置得到更新。

1）數(shù)據(jù)網(wǎng)格化處理，并將網(wǎng)格化后的數(shù)據(jù)離散為質(zhì)點(diǎn)，賦予物質(zhì)信息，再根據(jù)質(zhì)點(diǎn)分布，自動(dòng)建立計(jì)算用的背景網(wǎng)格。

2）對(duì)網(wǎng)格施加應(yīng)力邊界條件，建立網(wǎng)格節(jié)點(diǎn)位置的動(dòng)量方程并映射回質(zhì)點(diǎn)，得到質(zhì)點(diǎn)的速度梯度、應(yīng)變?cè)隽?；更新質(zhì)點(diǎn)的體積、應(yīng)力偏量、內(nèi)能，并考慮質(zhì)點(diǎn)的巖石物理、斷層（裂縫）等參數(shù)計(jì)算局部應(yīng)力，計(jì)算背景網(wǎng)格節(jié)點(diǎn)的應(yīng)力和應(yīng)變，進(jìn)而根據(jù)邊界條件修正、更新背景網(wǎng)格節(jié)點(diǎn)動(dòng)量；將背景網(wǎng)格節(jié)點(diǎn)位置、速度變化量映射回質(zhì)點(diǎn)，計(jì)算應(yīng)變量；放棄已變形的網(wǎng)格，重新產(chǎn)生背景網(wǎng)格，進(jìn)行迭代直至結(jié)果收斂。

3）通過(guò)上述步驟，可得應(yīng)力邊界條件下質(zhì)點(diǎn)的位移量、方向和應(yīng)變量等。各質(zhì)點(diǎn)位置的總位移方向代表局部最大水平主應(yīng)力方向。水平主應(yīng)力差異系數(shù)為最大、最小水平主應(yīng)力方向的應(yīng)變量之比。在射孔位置進(jìn)行水力壓裂時(shí)的施工壓力作為應(yīng)力邊界條件。模擬結(jié)果中應(yīng)變量代表的是所有應(yīng)力作用下的地層變形量，可表示水力壓裂的受效范圍。

2 MPM壓裂模擬

2.1 地質(zhì)背景

研究區(qū)位于四川盆地川東南高陡構(gòu)造帶陽(yáng)高寺構(gòu)造群，西北部為新店子構(gòu)造，中部為來(lái)蘇向斜，東南部為黃瓜山構(gòu)造。區(qū)內(nèi)主要發(fā)育北東向逆斷層，呈條帶狀分布；同時(shí)還發(fā)育北北東、近南北向逆斷層。規(guī)模較大的逆斷層均發(fā)育于構(gòu)造軸部及兩翼，對(duì)構(gòu)造形態(tài)起控制作用［8-13］。來(lái)蘇向斜內(nèi)部斷層欠發(fā)育，斷距較小。與斷裂相伴生的裂縫主要發(fā)育在研究區(qū)東南部及東北部構(gòu)造主體部位，且呈北東向展布；而中部裂縫欠發(fā)育。

本次壓裂的主要目的層為五峰組—龍一段，黑色頁(yè)巖富含有機(jī)質(zhì)，優(yōu)質(zhì)頁(yè)巖孔隙度介于4.3%～4.4%，脆性指數(shù)均大于45。氣藏埋深在3 500～4 500 m，屬深層頁(yè)巖氣藏［14］，其鉆井、壓裂等作業(yè)成本高［15-27］。為實(shí)現(xiàn)效益開(kāi)發(fā)，亟需開(kāi)展科學(xué)的壓裂模擬，為井位部署、完井設(shè)計(jì)以及壓裂設(shè)計(jì)優(yōu)化提供依據(jù)。

2.2 MPM壓裂模擬

1）收集整理目的層構(gòu)造圖和鉆、測(cè)井獲取的彈性模量、泊松比、密度等巖石物理和裂縫預(yù)測(cè)資料，質(zhì)點(diǎn)化并賦予相應(yīng)數(shù)值（見(jiàn)圖1。圖1a中空白處為數(shù)據(jù)缺失區(qū)，圖1a，1c中黑線為斷層，圖1b，1c中藍(lán)點(diǎn)為無(wú)斷層質(zhì)點(diǎn)，1c中紅點(diǎn)為斷層質(zhì)點(diǎn)）。網(wǎng)格面元為20 m×20 m，每個(gè)面元離散為64個(gè)質(zhì)點(diǎn)，除在質(zhì)點(diǎn)上賦予巖石物理等信息外，還要標(biāo)記其是否為斷層（裂縫）質(zhì)點(diǎn)。

2）根據(jù)研究區(qū) H202，H203，H204 等井區(qū)導(dǎo)眼井測(cè)井所獲得的初始地應(yīng)力參數(shù)（見(jiàn)表1），結(jié)合鄰區(qū)數(shù)據(jù)和世界地應(yīng)力數(shù)據(jù)庫(kù)（WSM2016），令區(qū)域最大水平主應(yīng)力方向?yàn)?10°，并給定應(yīng)力邊界條件。

表1 初始地應(yīng)力參數(shù)

3）將準(zhǔn)備好的數(shù)據(jù)和邊界條件輸入物質(zhì)點(diǎn)法數(shù)值模擬器，進(jìn)行模擬運(yùn)算，MPM壓裂模擬結(jié)果見(jiàn)圖2—5（圖2中彩色線為斷層；圖3—5中黑線為水平井軌跡，井筒周圍的紅色區(qū)域代表高剪切應(yīng)變區(qū)，其包絡(luò)代表壓裂影響范圍，散點(diǎn)是各井實(shí)測(cè)的微地震信號(hào)）。

4）應(yīng)用實(shí)際資料對(duì)模擬結(jié)果進(jìn)行評(píng)價(jià)。由圖2a可知，研究區(qū)模擬地應(yīng)力方向與成像測(cè)井解釋的基本一致，最大水平主應(yīng)方向以北西—南東向?yàn)橹鳎軘鄬佑绊懢植坑猩僭S轉(zhuǎn)向，如YY1井周圍。由表2可知，水平主應(yīng)力差異系數(shù)模擬結(jié)果與實(shí)驗(yàn)室數(shù)據(jù)誤差在8%以內(nèi)。H202井區(qū)水平主應(yīng)力差異系數(shù)主要分布在0.10～0.14，這有利于水力壓裂時(shí)形成復(fù)雜縫網(wǎng)（見(jiàn)圖3—5）。

表2 研究區(qū)水平主應(yīng)力差異系數(shù)對(duì)比

由圖3可知：MPM壓裂模擬的高剪切應(yīng)變區(qū)與微地震監(jiān)測(cè)結(jié)果基本吻合；H202井區(qū)不同井軌跡方向的高應(yīng)變區(qū)，局部呈現(xiàn)左右不對(duì)稱性，其右翼范圍大于左翼，整體與微地震分布特征基本一致。MPM壓裂模擬過(guò)程中，微地震信號(hào)分布不參與約束，僅作為后驗(yàn)。研究區(qū)5口壓裂施工井微地震事件分布與模擬結(jié)果吻合，進(jìn)一步驗(yàn)證了該方法的適用性和可靠性。

3 MPM壓裂模擬結(jié)果應(yīng)用

基于可靠的MPM壓裂模擬，在鉆井或壓裂前，可預(yù)測(cè)應(yīng)力方向、應(yīng)力差異系數(shù)分布、應(yīng)變范圍等。因此，結(jié)合壓裂模擬結(jié)果，在鉆井前可優(yōu)化水平井的軌跡部署，在壓裂前優(yōu)化壓裂施工設(shè)計(jì)參數(shù)，從而提高勘探開(kāi)發(fā)效益。模擬的應(yīng)變范圍還可以替代微地震監(jiān)測(cè)，從而降低工程成本，提高開(kāi)發(fā)效益。

本研究針對(duì)H202井區(qū)，分別選擇裂縫發(fā)育帶（裂縫方向與最大水平主應(yīng)力方向近垂直，下同）和網(wǎng)狀裂縫發(fā)育區(qū)，通過(guò)設(shè)置一系列井軌跡方位，開(kāi)展MPM壓裂模擬（見(jiàn)圖4—7），對(duì)優(yōu)化水平井井軌跡提出了合理化建議。由圖 4—7（圖 6b，6c，7b中角度為負(fù)，表示最大水平主應(yīng)力方向在水平井井軌跡左側(cè)；角度為正，則表示最大水平主應(yīng)力方向在水平井井軌跡右側(cè)）可以看出：裂縫發(fā)育帶，當(dāng)井軌跡方位與最大水平主應(yīng)力夾角在70°左右時(shí)，模擬改造面積最大；網(wǎng)狀裂縫發(fā)育區(qū)，當(dāng)井軌跡方位與最大水平主應(yīng)力近垂直時(shí)，模擬改造面積最大。該研究成果為研究區(qū)水平井井軌跡部署提供了依據(jù)。

4 結(jié)論

1）相比常規(guī)有限元法，基于MPM的數(shù)值模擬在理論上更適合頁(yè)巖儲(chǔ)層的壓裂模擬，H202井區(qū)的壓裂模擬結(jié)果與成像測(cè)井得到的應(yīng)力方向等實(shí)驗(yàn)數(shù)據(jù)、微地震信號(hào)分布吻合較好，驗(yàn)證了該方法在五峰組—龍一段優(yōu)質(zhì)頁(yè)巖儲(chǔ)層開(kāi)展壓裂模擬的適用性和可靠性。

2）MPM壓裂模擬研究得出不同方位水平井井軌跡的壓裂模擬改造面積，為水平井軌跡優(yōu)化部署及壓裂設(shè)計(jì)提供了依據(jù)。

3）MPM壓裂模擬為鉆井、壓裂前預(yù)測(cè)以及壓裂后評(píng)估提供了有效技術(shù)手段，為非常規(guī)油氣藏效益開(kāi)發(fā)提供技術(shù)支持，可推廣應(yīng)用。