蔡明金，張福祥，楊向同，彭建新，牛新年*

（1.西南石油大學博士后流動站，四川 成都 610500；2.中油塔里木油田分公司博士后工作站，新疆 庫爾勒 841000；3.中油塔里木油田分公司，新疆 庫爾勒 841000）

引 言

塔里木油田奧陶系碳酸鹽巖油氣藏發(fā)育隨機大型洞穴，與國內(nèi)外其他碳酸鹽巖油氣藏相比，具有其獨特性，其儲層具有極強的非均質性，是以儲集體和縫洞系統(tǒng)控藏為主的縫洞型油氣藏。目前，國內(nèi)外學者[1-11]在洞穴型碳酸鹽巖儲層試井解釋理論方面已做了大量研究工作，姚軍[1-2]、吳玉樹[3]等基于連續(xù)介質滲流理論建立了孔、洞、縫三重介質試井解釋模型，該模型是基于孔洞縫連續(xù)分布的，在隨機分布的大尺度洞穴儲層中無法得到很好的應用；修乃嶺[7]、林加恩[8]等研究了洞穴管流與滲流的機理，指出管流與滲流形式是統(tǒng)一的，證明應用滲流方法來研究大尺度洞穴是可行的；陳方方[9]基于純滲流理論將洞穴假定為一個區(qū)，滲透率無限大，建立了大尺度洞穴碳酸鹽巖均質復合油氣藏試井解釋模型。實際上洞穴是有充填的，滲透率無限大與地質情況不符，且由于洞穴形狀不確定，也無法解決生產(chǎn)中最關心的洞穴體積大小問題；彭小龍[10]、張利軍[11]等基于空腔流與管流理論建立了洞穴型碳酸鹽巖油氣藏試井解釋模型，該模型將地層中隨機分布大尺度溶洞抽象為幾個洞穴，雖然可以解釋洞穴的體積，但無法解釋洞穴距離，且除洞穴外的裂縫孔洞型儲層究竟有多大貢獻也不得而知，因此在實際試井解釋中也難以應用。塔里木油田大型洞穴碳酸鹽巖油氣藏多年的開發(fā)實踐表明，洞穴外的裂縫孔洞型儲層對生產(chǎn)也有貢獻，因此，在試井解釋中不但要解釋出生產(chǎn)中最關心的洞穴距離及大小，而且要研究洞穴外裂縫孔洞型儲層純滲流部分對產(chǎn)能的貢獻。因此，在前人研究成果[1-14]、[16-21]的基礎上，創(chuàng)新性地將裂縫孔洞型儲層滲流和洞穴空腔流動相結合，將油藏假定為內(nèi)區(qū)（大尺度洞穴發(fā)育區(qū)）雙重介質和外區(qū)（洞穴欠發(fā)育區(qū)）均質的兩區(qū)復合油藏，整個油藏滲流遵從達西定律，將地層純管流部分抽象為一個大型洞穴（由于洞穴在空間上形狀難以確定，忽略平面大小，只考慮洞穴位置和體積大?。?，洞穴內(nèi)為“空腔流”，這既可解釋出生產(chǎn)中關心的洞穴距離及大小問題，也可解釋出洞穴外儲層的滲流面積。

1 試井模型建立

1.1 物理模型的建立

模型基本假定如下：

（1）考慮圓形地層各向同性、水平等厚，地層厚度為h（m），外邊界距離為（m）。

（2）考慮地層呈雙重介質-均質兩區(qū)復合油藏特性分布（1區(qū)為縫洞發(fā)育高滲帶，2區(qū)為縫洞欠發(fā)育低滲帶，見圖1）。1區(qū)半徑為（m）；地層2區(qū)滲透率為K2（μ m2），孔隙度為?？紤]地層1區(qū)存在一個大型洞穴，處于位置A處，與井筒距離為rv（m），洞穴體積為V（m3），綜合壓縮系數(shù)為（MPa-1），洞穴到井筒的流動能量完全來自流體及巖石的彈性膨脹，流體通過長度為rv（m）、滲透率為Kf（μ m2）、寬度為（m）的裂縫直接流入井筒，洞穴壓力等于地層A點處壓力。

（3）除洞穴外，假定地層1區(qū)由天然裂縫和基質組成（圖2），天然裂縫是主要的流動通道，而基質作為流體的主要儲集空間，作為“源”，不斷地向天然裂縫供給，考慮基質向天然裂縫的竄流為Warren-Root擬穩(wěn)態(tài)竄流。基質、天然裂縫滲透率分別為Km、KNf；基質形狀因子為（m-2）。

（4）假定地層1區(qū)基質、天然裂縫的壓力系統(tǒng)相對獨立，儲層打開前，各處的壓力均為原始地層壓力，且考慮一種介質的孔隙度與另一種介質的壓力變化相對獨立，孔隙度分別為、。

（7）油井半徑為rw（m），以定產(chǎn)量q（m3/d）生產(chǎn)。

（8）流體在地層中滲流滿足達西定律。

（9）不考慮重力、毛管力的影響，考慮井儲和表皮影響，井筒儲集系數(shù)為C（m/MPa），表皮系數(shù)為S。

圖1 大型洞穴油藏物理模型

圖2 Warren-Root模型示意圖

1.2 數(shù)學模型的建立及求解

形聲字是現(xiàn)代漢字的主體，其中形符的表義功能，有助于留學生理解漢字的含義；聲符的示音功能，利于留學生了解讀音。對外漢語教師運用形聲字的這些優(yōu)勢進行教學，不僅可以減輕學生們的負擔，而且能夠增強他們的信心，激發(fā)他們學習漢字的熱情。

采用上述無因次變量（其他無因次變量定義參見文獻[9]），可建立如下無因次數(shù)學模型：

儲層控制方程如下所示。

1區(qū)：

2區(qū)：

式中：η12為1、2區(qū)導壓系數(shù)比。

根據(jù)達西定律及物質平衡原理，洞穴向井筒補給流量滿足以下公式：

洞穴處，洞穴與天然裂縫壓力相等：

內(nèi)邊界條件為：

1、2區(qū)界面連接處流量連續(xù)條件為：

1、2區(qū)界面連接處壓力連續(xù)條件為：

外邊界條件為：

對式（1）～（11）組成的方程組進行拉普拉斯變換求解，再進行stefest數(shù)值反演[15]，最終可求得井底壓力pwD，具體求解過程略。

2 試井曲線特征分析

圖3是洞穴型油藏考慮井筒儲集和表皮效應影響的無因次雙對數(shù)樣版曲線圖，可劃分為8個流動段：第Ⅰ階段為井儲和表皮效應影響段；第Ⅱ階段為大型洞穴向井筒供液段，無因次壓力導數(shù)曲線出現(xiàn)第1個凹陷；第Ⅲ階段為裂縫流段，無因次壓力及其導數(shù)曲線呈現(xiàn)1/2或1/4斜率平行直線；第Ⅳ階段為大型洞穴“空腔流”向1區(qū)天然裂縫系統(tǒng)徑向流過渡段，無因次壓力及其導數(shù)曲線下掉；第Ⅴ階段為1區(qū)基質向裂縫竄流段，無因次壓力導數(shù)曲線出現(xiàn)第2個凹陷；第Ⅵ階段為內(nèi)外區(qū)過渡段，無因次壓力及其導數(shù)曲線急劇上翹；第Ⅶ階段為1、2區(qū)總系統(tǒng)徑向流段，出現(xiàn)0.5水平線；第Ⅷ階段為外邊界反應段，表現(xiàn)為：無限大邊界，無因次壓力導數(shù)呈0.5水平線；封閉邊界，無因次壓力導數(shù)上翹；定壓邊界，無因次壓力導數(shù)出現(xiàn)下掉。

圖3 大尺度洞穴油藏試井典型樣版曲線

圖4 實例井無因次壓力及壓力導數(shù)雙對數(shù)擬合圖

從典型曲線特征可以看出，其典型曲線特征與常規(guī)雙重-均質兩區(qū)復合油藏明顯不同，不同之處在于：大型洞穴油藏多出第Ⅱ、Ⅲ、Ⅳ流動段，在Ⅱ段，無因次壓力導數(shù)曲線會出現(xiàn)第1個凹陷，凹子的深淺受洞穴的大小及與井筒距離控制，洞穴越大、與井筒距離越近，凹子就越深；在第Ⅲ段，無因次壓力及其導數(shù)曲線呈現(xiàn)平行直線，表現(xiàn)為裂縫流特征。

3 應用實例

塔里木盆地某油田奧陶系碳酸鹽巖油藏某油井，其基本參數(shù)如表1所示。

表1 模擬油井基本參數(shù)

該井地震反射為強串珠，在鉆井過程中發(fā)生放空漏失，最終裸眼完井。從地質靜態(tài)認識來看，該井鉆遇大型溶洞。該井于2008年5月進行壓力恢復測試，關井295.37h，關井前生產(chǎn)630h，日產(chǎn)油為54.1t/d。采用建立的洞穴型碳酸鹽巖雙重介質-均質兩區(qū)復合油氣藏無限大邊界模型進行解釋（圖4），擬合效果較好。解釋結果為：井筒儲存系數(shù)為0.55m3/MPa，表皮系數(shù)為0.7，地層1區(qū)半徑為80m，1區(qū)滲透率為156×10-3μ m2，儲容比為0.15，竄流系數(shù)為1.5×10-6m-2，2區(qū)滲透率為1.25×10-3μ m2，1、2區(qū)導壓系數(shù)比為100，洞穴距井筒5.3m，洞穴體積為23×104m3。解釋結果與地質靜態(tài)認識及生產(chǎn)情況吻合，解釋結果合理。

4 結 論

（1）建立并求解了碳酸鹽巖洞穴型雙重介質-均質兩區(qū)復合油藏試井模型，繪制了樣版曲線，分析了樣版曲線特征，同時用實例井進行了驗證。

（2）碳酸鹽巖洞穴型油藏試井樣版曲線不同于常規(guī)雙重介質-均質兩區(qū)復合油藏，不同之處為：洞穴型油藏無因次壓力導數(shù)曲線早期會出現(xiàn)一個凹陷，凹子的深淺受洞穴的大小及與井筒距離控制，洞穴越大、與井筒距離越近，凹子就越深；之后，壓力及其導數(shù)曲線呈現(xiàn)平行直線，表現(xiàn)為裂縫流特征。

（3）將裂縫孔洞滲流與洞穴“空腔流”結合起來建立試井模型，彌補前人建立的洞穴型儲層試井模型的缺陷，可為碳酸鹽巖洞穴型儲層油氣藏勘探開發(fā)提供理論和技術支撐。

[1]楊堅，姚軍，王子勝.三重介質復合油藏壓力特征研究[J].水動力學研究與進展，2005，20（4）：418-425.

[2]趙冬梅，姚軍，王子勝.壓敏三重介質油藏壓力響應特征[J].計算物理，2005，22（5）：444-448.

[3]戴衛(wèi)國，朱玉國，王子勝.遺傳算法在三重介質試井分析中的應用[J].油氣井測試，2005，14（2）：14-16.

[4]張訓華，高惠梅，劉先貴，等.三重介質油藏壓力動態(tài)分析的邊界元法[J].西安石油大學學報，2006，21（6）：35-39.

[5]姚軍，侯力群，李愛芬.天然裂縫性碳酸鹽巖封閉油藏產(chǎn)量遞減規(guī)律研究與應用[J].油氣地質與采收率，2005，12（1）：56-58.

[6]Wu Yushu.A triple-continuum pressure-transient model for a naturally fractured vuggy reservoir[C].SPE110044，2007：1-10.

[7]修乃嶺，熊偉，班凡生，等.縫洞型油藏裂縫-溶洞中流體運動特征研究[J].重慶科技學院學報：自然科學版，2007，9（2）：14-16.

[8]林加恩，李亮，楊慧珠.管流與滲流耦合流動理論研究初探[J].西南石油大學學報：自然科學版，2007，22（2）：11-15.

[9]張福祥，陳方方，彭建新，等.井打在大尺度溶洞內(nèi)的縫洞型油藏試井模型[J].石油學報，2008，30（6）：912-915.

[10]彭小龍，杜志敏，劉學利，等.大尺度溶洞裂縫型油藏試井新模型[J].西南石油大學學報：自然科學版，2008，30（2）：74-77.

[11]張利軍，程時清，梁繼德，等.縫洞型油藏試井新模型探討[J].大慶石油地質與開發(fā)，2010，29（6）：145-148.

[12]李雪梅，李藝，楊翔.圓形封閉雙層油藏新型產(chǎn)量曲線研究[J].特種油氣藏，2005，12（6）：40-43.

[13]王東權，賈永祿，林濤，等.雙重介質重整產(chǎn)量遞減分析模型與應用研究[J].特種油氣藏，2008，15（5）：70-73.

[14]李慶輝，陳勉，金衍，等.新型壓裂技術在頁巖氣開發(fā)中的應用[J].特種油氣藏，2012，19（6）：3-7.

[15]賈永祿，趙必榮.拉普拉斯變換及數(shù)值反演在試井分析中的應用[J].天然氣工業(yè)，1992，12（1）：60-64.

[16]王怒濤，孫翠容，郭金城，等. 火山巖氣藏儲集空間類型的動態(tài)描述及開發(fā)對策［J］.大慶石油地質與開發(fā)，2011，30（2）：75-79.

[17]侯瑋,郭杰,賀成才.徐深氣田火山巖氣藏動態(tài)描述追蹤方法［J］.大慶石油地質與開發(fā),2011,30（4）:76-79.

[18]康志江,趙艷艷,張冬麗,等.大尺度縫洞型碳酸鹽巖油藏彈性開采特征分析［J］.大慶石油地質與開發(fā),2011,30（5）:72-75.

[19]李愛芬，張東，高成海.封閉定容型縫洞單元注水替油開采規(guī)律[J].油氣地質與采收率，2012，19（3）：94-97.

[20]柳明，張士誠，牟建業(yè).碳酸鹽巖酸化徑向蚓孔擴展形態(tài)研究[J].油氣地質與采收率，2012，19（2）：106-110.

[21]劉顯太，王玉芹，杜玉山，等.背斜構造裂縫性碳酸鹽巖油藏裂縫網(wǎng)絡發(fā)育模式——以勝利油區(qū)富臺油田車古201潛山為例[J].油氣地質與采收率，2011,18（6）：70-73.