存在部分連通斷層的條帶狀油藏試井解釋新模型

曾 楊， 康曉東， 謝曉慶， 石 爻， 張烈輝

(1.海洋石油高效開發國家重點實驗室， 北京 100028； 2.中海油研究總院， 北京 100028； 3.西南石油大學 油氣藏地質及開發工程國家重點實驗室， 四川 成都 610500)

存在部分連通斷層的條帶狀油藏試井解釋新模型

曾 楊1,2， 康曉東1,2， 謝曉慶1,2， 石 爻1,2， 張烈輝3

(1.海洋石油高效開發國家重點實驗室， 北京 100028； 2.中海油研究總院， 北京 100028； 3.西南石油大學 油氣藏地質及開發工程國家重點實驗室， 四川 成都 610500)

條帶狀油藏是斷塊油田中最常見的構造形態之一，國內外學者對非封閉斷層的研究大都集中于均質油藏和復合油藏.因此，用現有的試井解釋模型對具有存在非封閉斷層的條帶狀油藏試井資料進行解釋所得到的結果往往并不理想.通過引入界面表皮的概念，將部分連通斷層視為一個無限薄的表皮邊界，建立了條帶狀油藏中存在部分連通斷層的試井解釋新模型，并利用Fourier指數變換和Laplace變換等方法進行求解，繪制了井底壓力響應特征曲線.分析表明：界面表皮越大，曲線上翹越多，當它足夠大時，表現為封閉斷層的特征；邊界距離對典型曲線形態的影響顯著；流度比、厚度比、導壓系數比主要影響壓力導數曲線上升和下降.

部分連通斷層； 條帶狀油藏； 界面表皮； 試井模型； 典型曲線

0 引言

國內外學者對封閉斷層[1-4]的研究較為系統和成熟.一般的，常用于判別封閉斷層的方法就是Horner半對數曲線法[5].1975年，Prasad和Tiab等[6-8]將存在一條封閉斷層邊界的模型擴展到多條交叉或者平行的封閉斷層模型.然而，目前對非封閉斷層的研究大都集中于均質油藏和復合油藏.Bixel等[9]首次提出將斷層視為非封閉邊界，但他所建立模型僅僅考慮巖石和流體性質在平面上發生突然變化的情形，Kuchuk等[10]后來對該模型的求解方法進行了完善.1984年，Stewart和Streltsova等[11，12]在Bixel研究的基礎上通過數值模擬研究了部分連通斷層對干擾試井的影響.Yaxely[13]、李傳亮[14]推導了無限大均質油藏中含有部分連通斷層模型的解析解.在此基礎上，Ambastha和張望明等[15，16]進一步研究了復合油藏中存在部分連通斷層邊界的井底壓力響應特征.

上述對條帶狀油藏中存在非封閉斷層的情況并沒有研究，而實際油藏由于沉積環境的影響，斷塊油田中最常見的構造形態之一為條帶狀油藏，它直接影響到后續壓力波傳播到斷層邊界的反應.

為此，本文通過引入Everdingeen[17]和Hurst[18]提出的界面表皮的概念，將部分連通斷層視為一個無限薄的表皮邊界，建立了存在部分連通斷層的條帶狀油藏試井解釋新模型，繪制了模型的井底壓力響應特征曲線并對相關的影響因素進行了分析.

1 滲流物理模型

針對兩條平行不滲透邊界的條帶狀油藏中存在部分連通斷層的情形，如圖1所示.假設斷層兩側的巖石特性、儲層厚度及滲透率不同，但同一區域內為均質油藏且各向同性.激動井為定產量線源，且流體為微可壓縮等溫滲流，流動服從達西定律.引入界面表皮的概念考慮兩個儲層區域之間斷面的傳導率，忽略重力和毛管力的影響.

圖1 存在部分連通斷層的條帶狀油藏示意圖

2 數學模型的建立及求解

2.1 數學模型的建立

根據圖1所建立的坐標系和上述假設條件，以基本的滲流力學理論為依據，根據前面的假設條件，激動井為定產量線源，且以產量q生產，開井前地層各處壓力相等，均為pi，根據點源函數理論，可推導出條帶狀油藏中存在部分連通斷層的無因次試井解釋模型.

(1)滲流微分方程

Ⅰ區，對xDgt;0：

(1)

Ⅱ區，對xDlt;0：

(2)

式(1)～(2)中，無因次變量分別為：

其中：p1D、p2D-Ⅰ區、Ⅱ區的無因次壓力；wD、xD、yD、aD、bD-無因次距離；tD-無因次時間；M-流度比；hD-厚度比；ηD-導壓系數比；p1、p2-Ⅰ區、Ⅱ區壓力；k1、k2-Ⅰ區、Ⅱ區地層滲透率；h1、h2-Ⅰ區、Ⅱ區地層厚度；μ1、μ2-Ⅰ區、Ⅱ區流體粘度；φ1、φ2-Ⅰ區、Ⅱ區孔隙度；Ci1、Ci2-Ⅰ區、Ⅱ區地層綜合壓縮系數；x、y-橫、縱坐標；a、b-井點橫、縱坐標；δ-δ函數表示定產量線源井.

(2)初始條件

p1D(xD,yD,0)=p2D(xD,yD,0)=0

(3)

(3)邊界條件

x方向的邊界條件：

p1D(∞D,yD,tD)=p2D(∞D,yD,tD)=0

(4)

y方向的邊界條件：

(5)

(4)連接條件

(6)

(7)

2.2 數學模型的求解

對上述無因次試井模型式(1)～(7)取基于yD的Fourier指數變換和基于tD的Laplace變換并化簡，可得：

(8)

式(8)中：

求解可得：

xDgt;0

(9)

xDlt;0

(10)

要求生產井的井底壓力，令xD=aD-1，yD=bD，得：

xDgt;0

(11)

3 典型曲線的繪制及特征分析

3.1 邊界距離的影響

圖2所示為井位于條帶狀地層中部(bD/wD=0.5)，aD、bD對典型曲線的影響.從圖2可以看出，無論aD、bD取值如何，典型曲線早期都表現出井儲階段的特征，壓力與壓力導數曲線重合，二者均為斜率為1的直線.隨著壓力波的不斷向外傳播，當tD較小時，壓力波尚未傳播到任何斷層邊界或界面，地層中表現出均質儲層滲流特征，壓力導數曲線表現為0.5水平線，該階段持續時間的長短取決于aD和bD中的較小值.

當aD/bDlt;1時，如圖2中的aD/wD=0.05情況，隨著tD不斷增大，壓力波首先傳播到兩區分界面處，由于此時考慮分界面的界面表皮S=106，這已接近于封閉斷層，因此壓力導數曲線上翹并出現第二個水平段，該水平段值趨于1，反映出分界面處為封閉斷層的特征，曲線上翹值的大小取決于界面表皮S值的大小，S越大，上翹越多，并最終趨于封閉斷層的特征.該水平段在壓力波傳播到平行斷層邊界時結束，其持續時間取決于井到平行斷層邊界的距離bD.由于井位于條帶狀地層中部(bD/wD=0.5)，故壓力波同時傳到兩平行斷層邊界處，之后地層中流動變為等效均質儲層中的線性流，壓力導數曲線表現為斜率為1/2的直線.

當aD/bD=1時，如圖2中的aD/wD=0.5情況，經過早期井儲流動和Ⅰ區均質儲層徑向流階段后，壓力波同時傳到兩平行斷層邊界和區域交界面處.由于交界面處S=106，已接近與封閉斷層，此時會出現一短暫的過渡階段，過渡段結束后，地層中流動進入等效均質儲層線性流階段，壓力導數曲線表現為斜率為1/2的直線.

圖2 aD和bD對典型曲線的影響(井位 于條帶狀地層中部)

當aD/bDgt;1時，如圖2中的aD/wD=2情況，經過早期井儲流動和Ⅰ區均質儲層徑向流階段后，隨著tD不斷增大，壓力波首先傳播到平行斷層邊界，地層中出現均質儲層線性流，壓力導數曲線表現為斜率1/2的直線，該階段一直持續到壓力波傳播到區域界面aD處為止.當壓力波繼續向外傳播至兩區界面時，由于界面處S=106，此時典型曲線上會出現短暫的過渡段，過渡段結束后，地層中流動進入等效均質儲層線性流階段，壓力導數曲線又出現斜率為1/2的直線.需要注意的是，第二次出現的1/2斜率直線段反映的是等效均質儲層中的線性流，故與之前出現的1/2斜率直線段并不重合，兩直線截距之差的大小與S值、流度比和導壓系數比等相關.

3.2 界面表皮S的影響

圖3顯示了當井靠近條帶狀地層的一條斷層邊界時，界面表皮S對井底壓力動態曲線的影響.由于aD/bDlt;1，壓力波首先傳到兩區分界面，界面表皮S主要影響壓力導數曲線上第二個水平段位置的高低，S越大，說明Ⅰ區與Ⅱ區的連通性越弱，流體在地層中流動的壓力損失就越大，壓力導數曲線由0.5水平線開始上翹，隨著S值的不斷增大，上翹越多.當S=106時，兩區分界面趨于封閉斷層，壓力導數曲線表現出值為1的水平線，隨后，當壓力波繼續傳到其中一條斷層時，此時形成角度為90 °的夾角斷層，壓力導數曲線表現出值為2的水平線，最后出現的反映等效均質儲層線性流的1/2斜率壓力導數曲線.當S=103時，兩區分界面部分連通，當壓力傳播到此分界面，由于界面表皮的存在，流體流動的壓力損失增大，壓力導數曲線逐漸從0.5水平線開始上翹，隨著Ⅱ區流體開始流動壓力有所下降，隨后，當壓力波繼續傳到其中一條斷層時，流動阻力再一次增加，壓力導數曲線出現上翹，直到壓力波傳播到另外一條斷層邊界，最后出現反映等效均質儲層線性流的1/2斜率壓力導數曲線.

圖3 界面表皮對典型曲線 的影響(aD/bDlt;1)

3.3 流度比、厚度比、導壓系數比的影響

圖4顯示了流度比M對井底壓力動態曲線的影響.當aD/bDlt;1時，流度比M主要影響壓力導數曲線上第二個水平段位置的高低，M越大，說明Ⅱ區物性越好，流體在地層中流動的壓力損失就越小，壓力導數曲線就下降越多，最后出現的反映等效均質儲層線性流的1/2斜率壓力導數曲線位置也越靠下.

厚度比hD與導壓系數比ηD對井底壓力動態的影響與流度比M對井底壓力動態的影響類似，這里就不在一一贅述.

圖4 流度比對典型曲線 的影響(aD/bDlt;1)

4 結論

(1)條帶狀油藏邊界直接影響到后續壓力波傳播到斷層后的反應，首次建立了條帶狀油藏中存在部分連通斷層邊界的試井解釋新模型，對具有類似邊界條件的油藏研究和應用具有一定的借鑒意義.

(2)邊界距離對存在部分連通斷層條帶狀油藏典型曲線形態的影響顯著；當界面表皮S值足夠大時，表現為封閉斷層的特征；流度比主要影響壓力波傳到Ⅱ區之后的壓力動態，流度比越大，壓力及壓力導數曲線的位置越靠下，厚度比與導壓系數比對井底壓力動態的影響與流度比類似.

(3)本文提出的斷層模型可用于雙重介質油藏、復合油藏等的試井解釋及壓力分布規律研究.

[1] 王 奎，楊 靜，侯光宗，等.斷層封閉性研究現狀及其發展趨勢[J].斷塊油氣田，2008，15(3)：43-45.

[2] 周 路，王麗君，羅曉容，等.斷層連通概率計算及其應用[J].西南石油大學學報(自然科學版)，2010，32(3)：11-18.

[3] 葉雙江，姜漢橋，孫傳宗，等.斷塊油藏水平井產能研究及斷層封堵影響分析[J].西南石油大學學報(自然科學版)，2011，33(5)：131-136.

[4] 張 旭，劉啟國，張小龍，等.斷層板狀油藏斜井試井模型研究[J].天然氣與石油，2012，30(5)：60-63.

[5] Horner D R.Pressure build-up in wells[C]//Presented at the 3rd World Petroleum Congress.Hague：SPE,1951：503-521.

[6] Prasad R K.Pressure transient analysis in the presence of two intersecting boundaries[J].JPT,1975,27(1)：89-96.

[7] Tiab D,Kumar A.Detection and location of two parallel sealing faults around a well[J].JPT,1980,32(10)：1 701-1 708.

[8] Earlougher R C.Pressure distributions in a rectangular reservoirs[J].JPT,1968,20(2)：199-208.

[9] Bixel H C,Larkin B K,Van Poolen H K.Effect of linear discontinuity on pressure buildup and drawdown behavior[J].JPT,1963,15(8)：885-895.

[10] Kuchuk F J,Habashy T M.Pressure behavior of laterally composite reservoir[J].SPEFE,1997,12(1)：47-56.

[11] Stewart G,Gupta A,Westaway P.The interpretation of interference tests in a reservoir with a sealing and a partially communicating faults[C]//Presented at the 1984 European PetroIeum Conference.London：SPE,1984：411-429.

[12] Streltsova T D,McKinley R M.Effect of flow time duration on buildup pattern for reservoirs with heterogeneous properties[J].SPEJ,1984,24(3)：294-306.

[13] Yaxely L M.Effect of partially communicating fault on transient pressure behavior[J].SPEFE,1987,2(4)：590-598.

[14] 李傳亮.非封閉斷層的試井解釋方法[J].新疆石油地質，1997，18(4)：370-376.

[15] Ambastha A K,McLeroy P G,Grader AS.Effect of a partially communicating fault in a composite reservoir on transient pressure testing[J].SPEFE,1989,4(2)：210-218.

[16] 張望明，曾 萍.存在局部連通斷層的復合油藏不穩定試井分析及典型曲線擬合方法[J].石油勘探與開發，1992，19(2)：69-76.

[17] Van Everdingen A F.The skin effect and its influence on the productivity capacity of a well[J].JPT,1953,5(6)：171-176.

[18] Hurst W.Establishment of the skin effect and its impediment to fluid flow in a wellbore[J].Pet.Eng,1953,25(8)：6-16.

【責任編輯：蔣亞儒】

Transientpressurebehavioranalysisinabandedreservoirwithpartiallycommunicatingfaults

ZENG Yang1,2， KANG Xiao-dong1,2， XIE Xiao-qing1,2， SHI Yao1,2， ZHANG Lie-hui3

(1.State Key Laboratory of Offshore Oil Exploitation, Beijing 100028, China； 2.CNOOC Research Institute, Beijing 100028, China； 3.State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation， Southwest Petroleum University， Chengdu 610500, China)

Banded reservoir is one of the most common structures in fault block reservoir,most of the researches on non-closed faults are concentrated in homogeneous reservoirs and composite reservoirs.When we used the existing well testing interpretation models to interpret the actual banded reservoir well test data with partially communicating fault boundary,the obtained results are not always ideal.A well testing model for wells located near a partially communicating fault in a two-zone,banded reservoir was established by introduced the “boundary skin” concept,the solution was obtained with the Fourier and Laplace transformation.Type curves were plotted by the Stehfest numerical inversion method and this model was verified by the examples.The result showed that the bigger the boundary skin value,the more the curve upturned;the boundary distance has a significant effect on the typical curve shape;the mobility ratio,thickness ratio,pressure conductivity coefficient mainly affected the pressure derivative curve rise and fall..

partially communicating fault； banded reservoir； boundary skin； well testing model； type curve

2017-08-03

國家科技重大專項項目(2016ZX05025-003)； 中海石油有限公司綜合科研項目(YXKY-2014-ZY-03)

曾 楊(1987-)，女，四川廣安人，工程師，碩士，研究方向：油氣田開發及提高采收率

2096-398X(2017)06-0104-05

TE353

A