王 蓓 ，劉向君 ，司馬立強(qiáng) ，徐 偉 ，李 騫 ，梁 瀚

（1.西南石油大學(xué)地球科學(xué)與技術(shù)學(xué)院，成都610500；2.中國石油西南油氣田分公司勘探開發(fā)研究院，成都610041）

0 引言

裂縫建模是一種常用的研究裂縫空間展布的手段，可分為確定性建模和隨機(jī)建模。確定性建模是根據(jù)已知信息建立確定的裂縫模型，如利用地震資料解釋出規(guī)模較大的裂縫，該方法一般不適用于規(guī)模較小的裂縫，且不能較好地綜合利用多種資料；隨機(jī)建模是利用裂縫的先驗(yàn)信息，通過隨機(jī)模擬方式生成可選的相同概率裂縫模型。上述建模方法不僅滿足已知點(diǎn)的裂縫統(tǒng)計(jì)學(xué)特征，而且承認(rèn)未知區(qū)域裂縫發(fā)育的隨機(jī)性，較好地尊重了裂縫模擬不確定性的客觀事實(shí)［1-3］。目前，國內(nèi)外主要的隨機(jī)裂縫建模方法大致可分為5類，即基于空間剖分的裂縫建模、離散裂縫網(wǎng)絡(luò)建模、基于變差函數(shù)的裂縫建模、基于多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)的裂縫建模以及基于分形特征迭代的裂縫建模［4-5］，然而，在實(shí)際應(yīng)用中，這些方法所需要的裂縫產(chǎn)狀、長(zhǎng)度、寬度等幾何特征的真實(shí)概率分布函數(shù)獲取困難，所建立的裂縫模型與實(shí)際地層中裂縫的發(fā)育情況差異較大。

四川盆地磨溪龍王廟組氣藏是迄今為止國內(nèi)外寒武系已探明縫洞型碳酸鹽巖儲(chǔ)層儲(chǔ)集規(guī)模最大的氣藏，具有油氣資源豐富、構(gòu)造低緩、孔洞縫配置關(guān)系復(fù)雜、儲(chǔ)層非均質(zhì)性較強(qiáng)、低孔-中高滲等特點(diǎn)。該氣藏 2016年初建成年產(chǎn) 90億m3的生產(chǎn)能力，開發(fā)井型均為定向井，隨著生產(chǎn)的平穩(wěn)推進(jìn)，氣藏南北兩翼邊水水侵情況凸顯［6］。為了及時(shí)掌握水體向氣藏內(nèi)部侵入的方向和方式，合理調(diào)整氣井產(chǎn)量，實(shí)現(xiàn)氣藏科學(xué)開發(fā)，亟須開展裂縫建模研究，精細(xì)描述儲(chǔ)層裂縫發(fā)育情況，明確氣藏高、低滲區(qū)域分布以及優(yōu)勢(shì)水侵方向與水侵方式，進(jìn)而優(yōu)化氣藏開發(fā)方式。

由于Petrel軟件將裂縫作為片元進(jìn)行處理，難以獲取裂縫在不同空間位置的產(chǎn)狀數(shù)據(jù)，且該軟件給出的片元過于理想，建立的角點(diǎn)網(wǎng)格模型可信度較低［7］，較難運(yùn)用于研究區(qū)的數(shù)值模擬。相較于角點(diǎn)網(wǎng)格模型，離散裂縫網(wǎng)格模型能更好地描述裂縫的復(fù)雜性和非均質(zhì)性，被廣泛運(yùn)用于氣藏模擬［4-5］。因此，擬以定向井的FMI成像測(cè)井裂縫識(shí)別技術(shù)為約束，進(jìn)行多尺度的DFM（discrete fracture model）非結(jié)構(gòu)化網(wǎng)格離散裂縫建模，解決隨機(jī)裂縫建模中裂縫參數(shù)獲取困難的問題，以構(gòu)建更加符合地質(zhì)實(shí)際的裂縫模型。

1 地質(zhì)概況

磨溪龍王廟組氣藏位于四川省遂寧市、資陽市及重慶市潼南縣境內(nèi)，構(gòu)造位置處于四川盆地川中古隆起平緩構(gòu)造區(qū)樂山—龍女寺古隆起東端。龍王廟組頂界構(gòu)造低緩、多高點(diǎn)，斷層發(fā)育程度較低（圖1），地層厚度為80～110 m。儲(chǔ)集層受有利亞相顆粒灘發(fā)育程度的控制［8-9］，連續(xù)性普遍較好，但局部區(qū)域物性相對(duì)較差；受地層埋深壓實(shí)和多期溶蝕成巖作用的改造，次生孔、洞、縫極為發(fā)育［10-11］。分析巖心、鑄體薄片和數(shù)字巖心等資料表明，研究區(qū)宏觀裂縫發(fā)育構(gòu)造縫、壓溶縫以及構(gòu)造溶蝕縫；構(gòu)造縫以高角度縫為主，斜交縫、水平縫以及網(wǎng)狀縫發(fā)育程度均較低；微裂縫較發(fā)育，發(fā)育頻率達(dá)到了40%，在部分薄片中可見瀝青、黃鐵礦充填或者半充填［12］。

圖1 磨溪龍王廟組氣藏相干地層切片F(xiàn)ig.1 Coherent stratigraphic slice of gas reservoir of Longwangmiao Formation in Moxi area

2 大斜度井和水平井裂縫識(shí)別技術(shù)

通過巖心觀察和薄片鑒定的裂縫僅是小尺度的、局部的、片面的存在，為了獲得實(shí)鉆井未取心井段的裂縫參數(shù)，約束裂縫模型的建立，提高裂縫預(yù)測(cè)的精度，須利用FMI成像測(cè)井技術(shù)對(duì)有效的裂縫進(jìn)行精細(xì)識(shí)別［13-15］。在以往的研究中，基于直井的裂縫已形成一套系統(tǒng)的識(shí)別和評(píng)價(jià)技術(shù)［14-17］，鑒于磨溪龍王廟組氣藏開發(fā)井均為大斜度井或水平井，著重針對(duì)該類定向井型中FMI成像測(cè)井裂縫識(shí)別技術(shù)和評(píng)價(jià)開展研究。

2.1 大斜度井和水平井裂縫的識(shí)別

定向井受井斜、地層傾角的影響，其井眼軌跡分為下穿地層和上穿地層2類。定向井FMI成像測(cè)井圖天然裂縫的識(shí)別模式與直井存在明顯差異。為了精細(xì)描述裂縫發(fā)育特征，須在排除巖性等其他地質(zhì)特征成像響應(yīng)的基礎(chǔ)上，結(jié)合常規(guī)測(cè)井資料及區(qū)域構(gòu)造資料，加以精細(xì)判別。

假設(shè)模型中地層傾角為2°～10°，按照裂縫與巖心中線垂直面的夾角分類，研究區(qū)發(fā)育高角度縫、斜交縫和低角度縫（表1）：①高角度縫傾角為75°～90°，與定向井下穿地層和上穿地層時(shí)垂直切割，切割井眼很短，在成像測(cè)井圖中表現(xiàn)為1組或2組呈對(duì)稱狀的幅度很低、傾向相差約90°的正弦曲線，該正弦曲線不能指示裂縫的最大主應(yīng)力方向；深淺雙側(cè)向測(cè)井曲線呈小幅度正差異，與鉆井誘導(dǎo)縫的成像響應(yīng)特征相似。直井中高角度縫在成像測(cè)井圖上表現(xiàn)出的特征與定向井中的低角度縫特征具有相似性。②斜交縫傾角為15°～75°，在定向井下穿地層和上穿地層時(shí)，隨著井斜的增加，與井眼斜交的面越大，成像測(cè)井圖中顯示的波谷狀高導(dǎo)異常也越高。隨著斜交縫角度及方位的變化，其特征也隨之變化。直井中斜交縫在成像測(cè)井曲線上表現(xiàn)為暗色正弦曲線。③低角度縫傾角為0°～15°，在定向井下穿地層和上穿地層時(shí)由于層界面和井眼斜交，切割井眼較長(zhǎng)，成像測(cè)井圖表現(xiàn)為幅度很高的波峰狀和波谷狀高導(dǎo)異常。在井軌跡平行于層界面時(shí)表現(xiàn)出變形的縱向彎曲曲線。常規(guī)深淺雙側(cè)向測(cè)井曲線表現(xiàn)為正差異，與層理成像響應(yīng)特征相似。直井中低角度縫在成像測(cè)井圖上表現(xiàn)為近直線狀互相平行的高導(dǎo)異常，深淺雙側(cè)向測(cè)井曲線呈小幅度正差異，與定向井中的高角度縫特征具有相似性。

表1 大斜度井和水平井FMI成像測(cè)井裂縫識(shí)別圖版Table 1 FMI imaging logging fracture recognition chart for highly deviated wells and horizontal wells

2.2 大斜度井和水平井裂縫的定量評(píng)價(jià)

根據(jù)大斜度井和水平井FMI成像測(cè)井裂縫識(shí)別圖版（表1），在定性識(shí)別裂縫單井發(fā)育位置的基礎(chǔ)上，可定量評(píng)價(jià)裂縫發(fā)育的產(chǎn)狀、開度、密度及孔隙度等參數(shù)，進(jìn)而綜合評(píng)價(jià)裂縫的發(fā)育程度，即巖石的破裂程度，以達(dá)到精細(xì)表征井點(diǎn)裂縫發(fā)育情況的目的，為裂縫建模奠定基礎(chǔ)。

對(duì)電成像圖像進(jìn)行處理，拾取裂縫，再計(jì)算裂縫特征參數(shù)，對(duì)于裂縫建模具有重要的約束作用。裂縫的發(fā)育程度是單位體積內(nèi)裂縫發(fā)育情況的綜合反映，裂縫的密度、張開度、孔隙度可通過電成像處理獲取，裂縫的長(zhǎng)度可通過雙側(cè)向測(cè)井等探測(cè)深度較深的測(cè)井信息獲取，裂縫的滲透性可通過陣列聲波斯通利波能量獲取。利用有效裂縫密度、有效裂縫孔隙度、有效裂縫張開度和有效裂縫長(zhǎng)度分別與其權(quán)重的乘積之和的方法，建立了一個(gè)適用于研究區(qū)的綜合裂縫指數(shù)計(jì)算模型。該模型中裂縫密度、孔隙度、張開度和長(zhǎng)度的下限均為定值；有效裂縫密度對(duì)裂縫滲透性的影響較大，經(jīng)過多次權(quán)重系數(shù)迭代擬合，確定該參數(shù)所占權(quán)重最大，為0.5；裂縫的徑向延伸度對(duì)裂縫發(fā)育程度的影響較大，但是該參數(shù)無法依據(jù)測(cè)井資料直接獲取，所以在該模型中加入了FMI成像測(cè)井解釋的裂縫長(zhǎng)度，有效裂縫長(zhǎng)度對(duì)裂縫發(fā)育程度的影響相對(duì)較小，所以經(jīng)驗(yàn)值擬合結(jié)果顯示該參數(shù)所占權(quán)重最小，為0.1。FI=（FD-FDM）/FD×W1+（FP-FPM）/FP×W2+

（FA-FAM）/FA×W3+（FL-FLM）/FL×W4 （1）式中：FI為裂縫發(fā)育指數(shù)；FD為計(jì)算裂縫密度（條/m）；FDM為裂縫密度下限（條/m）；FP為計(jì)算裂縫孔隙度，%；FPM為裂縫孔隙度下限，%；FA為計(jì)算裂縫張開度，μm；FAM為裂縫張開度下限，μm；FL為計(jì)算裂縫長(zhǎng)度，m/m；FLM為裂縫長(zhǎng)度下限，m/m；W1為裂縫密度指數(shù)權(quán)重；W2為裂縫孔隙度指數(shù)；W3為裂縫寬度指數(shù)權(quán)重；W4為裂縫長(zhǎng)度指數(shù)權(quán)重；W1+W2+W3+W4=1。

3 多尺度非結(jié)構(gòu)化網(wǎng)格裂縫建模技術(shù)

對(duì)于碳酸鹽巖縫洞型儲(chǔ)層發(fā)育的氣藏，裂縫的發(fā)育程度對(duì)氣藏連通性、儲(chǔ)集性的影響均較大。利用基于定向井的FMI成像測(cè)井裂縫識(shí)別技術(shù)，在單井有效儲(chǔ)集空間精細(xì)描述的基礎(chǔ)上，綜合測(cè)井資料及疊前地震各向異性裂縫預(yù)測(cè)和不連續(xù)性檢測(cè)等數(shù)據(jù)，結(jié)合三維地質(zhì)建模和離散裂縫建模技術(shù)，可實(shí)現(xiàn)多尺度縫洞型儲(chǔ)層地質(zhì)建模，形成多尺度DFM非結(jié)構(gòu)化網(wǎng)格離散裂縫模型，這對(duì)于刻畫氣藏高、低滲區(qū)域分布及優(yōu)勢(shì)水侵通道，優(yōu)化水侵方式具有重要的指導(dǎo)意義［18-21］。

3.1 大尺度裂縫建模

大尺度裂縫是指裂縫長(zhǎng)度在公里級(jí)的大尺度斷裂［19］。為了了解更加接近真實(shí)情況的強(qiáng)非均質(zhì)性裂縫系統(tǒng)，本次研究基于不連續(xù)性檢測(cè)體，通過橢圓或者線性擬合得到大尺度裂縫元的方位角，然后根據(jù)裂縫密度和方位角重構(gòu)大尺度離散裂縫元，再連接離散裂縫元，形成真實(shí)的裂縫網(wǎng)絡(luò)，進(jìn)而通過對(duì)大尺度裂縫的約束剖分非結(jié)構(gòu)化網(wǎng)格?；诓贿B續(xù)性檢測(cè)體的方位角擬合是根據(jù)局部區(qū)域內(nèi)裂縫密度分布，利用最小二乘法擬合目標(biāo)網(wǎng)格最有可能的方向。非結(jié)構(gòu)化網(wǎng)格技術(shù)中的四面體網(wǎng)格較雙重孔隙介質(zhì)建模過程中常用的角點(diǎn)網(wǎng)格更靈活，可以適應(yīng)復(fù)雜斷層或者裂縫系統(tǒng)的幾何約束，實(shí)現(xiàn)利用相對(duì)少的網(wǎng)格數(shù)精細(xì)表征裂縫，不僅提高了裂縫建模的準(zhǔn)確性，還大大縮短了數(shù)值模擬的運(yùn)算時(shí)間。

對(duì)于構(gòu)造裂縫或壓裂裂縫等大尺度裂縫，利用離散裂縫模型進(jìn)行模擬，關(guān)鍵在于確定其滲透率或?qū)Я飨禂?shù)。由于裂縫的開度遠(yuǎn)小于其長(zhǎng)度和高度，可采用無限平板中的流體流動(dòng)模型確定裂縫的滲透率。網(wǎng)格的傳導(dǎo)率與滲透率類似，可認(rèn)為是一個(gè)與流體性質(zhì)無關(guān)的屬性，因此，可將其假設(shè)為單相流動(dòng)問題，且忽略重力的影響。對(duì)于非結(jié)構(gòu)化網(wǎng)格，相鄰的2個(gè)網(wǎng)格可能正交，也可能不正交，網(wǎng)格中心與接觸面中點(diǎn)的連線并不一定垂直，每個(gè)網(wǎng)格的相鄰網(wǎng)格數(shù)量也不確定，因此，需要根據(jù)鄰接網(wǎng)格的數(shù)量，分別描述基質(zhì)與基質(zhì)、基質(zhì)與裂縫、裂縫與裂縫之間的關(guān)系（圖2）。

圖2 二維模型中的3種傳導(dǎo)率Fig.2 Three kinds of conductivity in two-dimensional model

3.2 中小尺度裂縫建模

中小尺度裂縫是指裂縫長(zhǎng)度在米級(jí)—百米級(jí)的裂縫［19］。中小尺度非結(jié)構(gòu)化離散裂縫建模的思路與大尺度裂縫建模的思路基本一致，但仍存在2個(gè)方面的差異：①模型建立基于不同的地震數(shù)據(jù)體；②獲取裂縫走向和密度的方法不同。首先，在獲取疊前地震各向異性裂縫預(yù)測(cè)體強(qiáng)度和方位角資料的基礎(chǔ)上，利用成像測(cè)井裂縫參數(shù)對(duì)地震預(yù)測(cè)的裂縫走向和密度進(jìn)行校正，得到符合井眼數(shù)據(jù)的裂縫方位角和密度數(shù)據(jù)體，再利用裂縫重構(gòu)算法提取中小尺度裂縫元；其次，根據(jù)局部裂縫單元的展布和密度的連續(xù)性，將最有可能處于一條裂縫上的裂縫單元在橫向或縱向上連接起來，裂縫橫向連接受相鄰裂縫元夾角和距離控制，縱向連接在算法上有2個(gè)條件，即裂縫之間的距離足夠接近，以及用于連接上下層之間的裂縫單元的傾角應(yīng)盡可能與被連接的裂縫單元一致；最終，對(duì)中小尺度裂縫作等效滲透率處理，以裂縫帶形式體現(xiàn)中小尺度裂縫的展布。對(duì)于中小尺度的裂縫一般采用等效處理方式，將裂縫介質(zhì)的滲流屬性考慮到基巖介質(zhì)中，其方法有實(shí)驗(yàn)測(cè)量法、數(shù)字巖心分析法以及數(shù)值等效法［22-24］。磨溪龍王廟組碳酸鹽巖儲(chǔ)層微小裂縫較發(fā)育，因此，主要采用數(shù)值等效法，通過將中小尺度裂縫滲透率等效為基質(zhì)滲透率的方法體現(xiàn)中小尺度裂縫對(duì)氣藏的貢獻(xiàn)。該方法將裂縫介質(zhì)的裂縫孔隙度和裂縫滲透率等效至基質(zhì)孔隙度和基質(zhì)滲透率中，其形式上與雙孔雙滲模型一致，區(qū)別在于雙孔雙滲模型計(jì)算得到的滲透率為裂縫網(wǎng)格滲透率，而微小裂縫等效計(jì)算得到的滲透率為基質(zhì)增強(qiáng)滲透率（圖3）。

圖3 微小裂縫滲透率等效計(jì)算示意圖Fig.3 Schematic diagram for equivalent calculation of permeability of micro-fracture

4 應(yīng)用實(shí)例

磨溪龍王廟組氣藏儲(chǔ)層非均質(zhì)性較強(qiáng)，縫洞配置關(guān)系復(fù)雜。在利用FMI成像測(cè)井進(jìn)行裂縫有效評(píng)價(jià)的基礎(chǔ)上，采用多尺度非結(jié)構(gòu)化網(wǎng)格裂縫建模技術(shù)，對(duì)研究區(qū)大尺度裂縫和中小尺度裂縫分別進(jìn)行了預(yù)測(cè)，所預(yù)測(cè)裂縫發(fā)育的高滲帶分布得到了后續(xù)開發(fā)井的證實(shí)，所預(yù)測(cè)的中小尺度裂縫與單井成像測(cè)井解釋結(jié)果吻合程度較高，為研究區(qū)高、低滲區(qū)域分布和水侵優(yōu)勢(shì)通道研究奠定了基礎(chǔ)。

4.1 裂縫模型可靠性評(píng)價(jià)

在磨溪龍王廟組氣藏基于疊前地震各向異性預(yù)測(cè)和基于廣義S變換譜分解的不連續(xù)性檢測(cè)的基礎(chǔ)上，以基于定向井的FMI成像測(cè)井裂縫識(shí)別技術(shù)為約束，三維重構(gòu)了大尺度裂縫和中小尺度裂縫的發(fā)育情況，結(jié)果表明研究區(qū)大尺度裂縫和中小尺度裂縫均較發(fā)育（圖4）。大尺度裂縫約束非結(jié)構(gòu)化剖分，計(jì)算大尺度裂縫與基質(zhì)傳導(dǎo)率，從而建立大尺度離散裂縫模型；對(duì)于中小尺度裂縫，通過將裂縫介質(zhì)滲透率等效到基質(zhì)滲透率的方法，建立離散裂縫模型。針對(duì)大尺度裂縫進(jìn)行離散裂縫建模，剖分的網(wǎng)格為非結(jié)構(gòu)化網(wǎng)格，總的網(wǎng)格數(shù)約為217.8萬個(gè)，其中基質(zhì)網(wǎng)格數(shù)約為190.5萬個(gè)，裂縫網(wǎng)格數(shù)約為27.3萬個(gè)；針對(duì)中小尺度裂縫進(jìn)行等效處理，將小尺度裂縫屬性等效到剖分的非結(jié)構(gòu)化基質(zhì)網(wǎng)格中。

圖4 磨溪龍王廟組氣藏多尺度裂縫模型Fig.4 Multi-scale fracture model for gas reservoir of Longwangmiao Formation in Moxi area

以位于MX8井區(qū)的MX009-8-X1井為例，在裂縫模型中可見該井周紅色線條表示的大尺度裂縫基本不發(fā)育，藍(lán)色線條表示的中小尺度裂縫相對(duì)欠發(fā)育［圖5（a）］；玫瑰花圖顯示裂縫方位為南北向和北西—南東向［圖5（b）］；成像測(cè)井響應(yīng)顯示該井縱向僅發(fā)育14條裂縫且多為高角度縫，溶蝕孔、洞較發(fā)育［圖5（c）］；試井壓力恢復(fù)雙對(duì)數(shù)曲線顯示該井具有視邊界特征［圖5（d）］，儲(chǔ)層滲透率為14.4 mD，具有中滲特征；酸化測(cè)試獲超過100萬m3/d的工業(yè)氣流，井產(chǎn)量穩(wěn)定中—高產(chǎn)。通過將MX009-8-X1井動(dòng)、靜態(tài)資料與實(shí)鉆井周緣裂縫模型對(duì)比分析表明，該井裂縫展布情況與生產(chǎn)動(dòng)態(tài)特征相吻合。

圖5 驗(yàn)證井MX009-8-X1及其周緣裂縫發(fā)育情況Fig.5 Fracture development in validation well MX009-8-X1 and adjacent areas

磨溪龍王廟組氣藏目前完鉆的定向井共計(jì)35口，其中34口井錄取了FMI成像測(cè)井資料?；诤罄m(xù)實(shí)鉆井的FMI成像測(cè)井裂縫響應(yīng)、裂縫參數(shù)、試井解釋等資料，證實(shí)了所建立磨溪龍王廟組多尺度裂縫模型的可靠性。

4.2 高、低滲區(qū)域分布及水侵通道研究

將常規(guī)測(cè)井儲(chǔ)層參數(shù)解釋、均方根振幅屬性圖和地震剖面等資料相結(jié)合，初步編繪了研究區(qū)高、低滲區(qū)域分布圖，再利用氣井試井解釋滲透率參數(shù)場(chǎng)等生產(chǎn)動(dòng)態(tài)資料，修正縫洞的規(guī)模、邊界和連通性，最終明確了研究區(qū)高、低滲區(qū)域分布（圖6）。磨溪龍王廟組氣藏高滲區(qū)呈連片狀廣泛分布，儲(chǔ)層連續(xù)性較好，地震響應(yīng)同相軸表現(xiàn)為連續(xù)的強(qiáng)波峰“亮點(diǎn)”反射，試井解釋滲透率高達(dá)500 mD；相對(duì)低滲區(qū)多呈條帶狀或片狀局限分布，儲(chǔ)層連續(xù)性較好，地震響應(yīng)同相軸表現(xiàn)為較連續(xù)的弱波峰反射或雜亂反射，試井解釋滲透率以大于0.5 mD為界限，與高滲區(qū)大面積接觸；低滲區(qū)受巖性影響呈局限的塊狀分布，儲(chǔ)層連續(xù)性相對(duì)較差，地震同相軸表現(xiàn)為斷續(xù)的雜亂反射或連續(xù)的波谷“暗點(diǎn)”反射。研究區(qū)生產(chǎn)井主要分布于氣藏統(tǒng)一氣水界面-4 385 m構(gòu)造等值線之上的高滲區(qū)。

在明確了研究區(qū)高、低滲區(qū)分布的基礎(chǔ)上，綜合巖心、鑄體薄片、測(cè)井、地震及生產(chǎn)動(dòng)態(tài)等資料，將研究區(qū)水侵通道具體劃分為4個(gè)方向的9條通道（圖6）。依據(jù)生產(chǎn)井水氣比等產(chǎn)水特征，將水侵模式歸納為沿裂縫水竄型和沿溶蝕孔洞均勻推進(jìn)型（圖 7）。

沿裂縫水竄型氣井下部強(qiáng)烈發(fā)育裂縫發(fā)育帶，以MX009-3-X2井為例，該井中小尺度裂縫重構(gòu)和等效滲透率均反映出井軌跡下部地層中存在發(fā)育的北西、北東向裂縫帶，大尺度裂縫重構(gòu)表明井軌跡下部地層中存在發(fā)育的北東向大裂縫帶與水體溝通，大尺度裂縫通過中小尺度裂縫帶與模型上層氣體溝通；地震疊前偏移時(shí)間剖面顯示該井儲(chǔ)層與氣藏北翼高滲水體儲(chǔ)層連續(xù)性好，裂縫FMI成像測(cè)井圖顯示裂縫非常發(fā)育，裂縫方位為北東—南西向，試井結(jié)果表現(xiàn)為高滲特征（滲透率大于60 mD），生產(chǎn)動(dòng)態(tài)表現(xiàn)為隨著產(chǎn)量增高，水氣比持續(xù)上升（表 2、表 3）。

圖6 磨溪龍王廟組氣藏優(yōu)勢(shì)水侵通道分布Fig.6 Distribution of dominant water invasion channel of gas reservoir of Longwangmiao Formation in Moxi area

圖7 磨溪龍王廟組氣藏優(yōu)勢(shì)水侵模式Fig.7 Water invasion patterns of gas reservoir of Longwangmiao Formation in Moxi area

沿溶蝕孔洞均勻推進(jìn)型氣井井軌跡下部地層中強(qiáng)烈發(fā)育溶蝕孔洞帶，裂縫相對(duì)欠發(fā)育。以MX8井為例，該井巖心、成像測(cè)井等資料以及大尺度和中小尺度裂縫模型均表現(xiàn)出裂縫欠發(fā)育的特征，并且沿裂縫發(fā)育方向南部的儲(chǔ)層分布較為局限，據(jù)靜態(tài)資料確認(rèn)該區(qū)域發(fā)生裂縫水竄的可能性較小，試井解釋表現(xiàn)出高滲的特征，表明邊水沿層推進(jìn)相對(duì)均勻（表 2、表 3）。

根據(jù)模型預(yù)測(cè)結(jié)果分析并判斷出氣藏邊水可能沿裂縫發(fā)育的高滲通道侵入氣藏，后續(xù)開發(fā)生產(chǎn)動(dòng)態(tài)證實(shí)了模型預(yù)測(cè)的水侵方向和通道與實(shí)際情況相符合。

表2 不同水侵模式典型井地球物理響應(yīng)特征Table 2 Geophysical characteristics of typical wells with different water invasion patterns

表3 不同水侵模式典型井試井特征Table 3 Well test characteristics of typical wells with different water invasion patterns

5 結(jié)論

（1）利用大斜度井和水平井裂縫識(shí)別技術(shù)，針對(duì)磨溪龍王廟組氣藏進(jìn)行了高角度縫、斜交縫、低角度縫在FMI成像測(cè)井響應(yīng)中井軌跡下穿地層和上穿地層2種情況下裂縫的識(shí)別及裂縫參數(shù)的定量評(píng)價(jià)，評(píng)價(jià)結(jié)果表明較直井FMI成像測(cè)井顯示可獲得更多的裂縫發(fā)育信息。

（2）利用地震疊后不連續(xù)性檢測(cè)體，根據(jù)裂縫密度和方位角可重構(gòu)大尺度裂縫；利用疊前方位角各向異性裂縫預(yù)測(cè)體，結(jié)合定向井FMI成像測(cè)井響應(yīng)，可重構(gòu)中小尺度裂縫?；诖笮倍染退骄瓼MI成像測(cè)井約束下的多尺度非結(jié)構(gòu)化網(wǎng)格裂縫模型，能更好地反映實(shí)際地層裂縫發(fā)育情況，實(shí)鉆結(jié)果證實(shí)了所建立的離散裂縫模型的可靠性。

（3）精細(xì)描述磨溪龍王廟組氣藏高滲、相對(duì)低滲和低滲區(qū)域分布范圍，刻畫出4個(gè)方向的9條水侵優(yōu)勢(shì)通道，建立了沿裂縫水竄型和沿溶蝕孔洞均勻推進(jìn)型2種水侵模式。研究結(jié)果不僅對(duì)該區(qū)氣藏產(chǎn)能補(bǔ)充井位部署、開發(fā)技術(shù)對(duì)策的優(yōu)化調(diào)整具有支撐作用，還對(duì)同類特大型有水深層碳酸鹽巖縫洞型儲(chǔ)層氣藏開展裂縫精細(xì)描述、水侵優(yōu)勢(shì)通道刻畫和水侵模式建立等研究具有借鑒意義。