王小魯，賈永祿，王本成，王少飛，聶仁仕

（1.油氣藏地質(zhì)及開發(fā)工程國家重點實驗室（西南石油大學(xué)），成都610500；2.青海油田 天然氣開發(fā)公司，青海 格爾木816000；3.中國石化 西南油氣分公司 博士后科研工作站、勘探開發(fā)研究院，成都610041；4.中國石油 長慶油田分公司，西安710000）

隨著油氣藏的開發(fā)，動靜態(tài)結(jié)合的多井系統(tǒng)油田開發(fā)的試井過程中，后期試井資料受鄰井干擾的現(xiàn)象日趨增加［1］。然而，國內(nèi)外所應(yīng)用的試井分析方法和解釋軟件大多是基于單井系統(tǒng)的試井解釋模型，忽略了鄰井干擾的影響，難以得到正確的分析結(jié)果。

干擾試井是多井試井中歷史最長的一種方式，Theis率先給出了在均質(zhì)無限大油藏中由其他點流速變化引起的壓力變化的解［2，3］。Tongpenyai提出了考慮激動井表皮效應(yīng)和井筒儲存效應(yīng)的均質(zhì)油藏干擾試井模型［4］。Ogbe等進一步給出了注采井同時存在的均質(zhì)多井油藏中的一口生產(chǎn)井的壓力恢復(fù)分析理論，并建立了一口注水井的壓力降落分析方法［5］。張秀華等研究了考慮井筒儲集和表皮系數(shù)影響的復(fù)合油藏干擾試井模型，并繪制了相應(yīng)的樣板曲線［6］。張德志等研究了三重介質(zhì)油藏干擾試井壓力動態(tài)變化，并分析了井筒儲存、表皮系數(shù)、竄流系數(shù)以及彈性儲容比對觀測井井底壓力的影響［7］。謝林峰等利用數(shù)值模擬的方法開展了氣藏水平井干擾試井參數(shù)敏感性研究，定性地揭示了這些參數(shù)或參數(shù)組合與干擾壓力的關(guān)系，但不能用于干擾試井解釋與評價［8］。

目前國內(nèi)外許多文獻［9－14］對干擾試井問題作了一定的論述，建立了相應(yīng)的理論模型和數(shù)學(xué)模型，部分還作出了相應(yīng)的理論曲線，在油田中的應(yīng)用也比較普遍；但大都針對直井模型，水平井模型研究甚少。針對試井模型求解方法而言，大都采用Laplace變換求解［15］。本文采用一種新的水平井試井求解方法，獲得了三維實空間任意一點的壓力精確解，利用新研究的水平井干擾試井分析模型和試井曲線，有效地確定地層的連通情況以及準(zhǔn)確獲取地層動態(tài)參數(shù)信息，為油田高效開發(fā)提供保障。

1 物理模型

圖1 水平井干擾試井滲流模型示意圖Fig.1 Schematic diagram of the physical model of the horizontal well interference testing

無限大均質(zhì)油藏中，考慮兩水平井A和B，A為觀測井，B為激動井，其滲流物理模型如圖1所示。且基于如下物理模型假設(shè)：（1）油藏厚度為h，水平井距油藏底面距離為Zw，水平井長度為L；（2）油藏均質(zhì)且各向異性，水平等厚，頂面和底面為封閉邊界，外邊界為無限大；（3）激動井開井前整個油藏保持均一的原始壓力pi；（4）地層巖石和流體均微可壓縮，且壓縮系數(shù)為常數(shù)，原油黏度不變；（5）忽略重力、毛管力的影響，流體在地層中作達西流動。

2 數(shù)學(xué)模型

在三維笛卡爾坐標(biāo)系下，多段壓裂水平井滲流控制微分方程為

若將井作為點源列入方程，式（1）可改寫為

初始條件

外邊界條件

底面封閉

頂面封閉

3 無因次數(shù)學(xué)模型

無因次定義如下

無因次壓力

無因次時間

無因次流量

無因次距離

在三維笛卡爾坐標(biāo)系下，水平井無因次點源滲流控制微分方程為

初始條件

外邊界條件

底面封閉

頂面封閉

4 數(shù)學(xué)模型求解

根據(jù)水平井無因次點源滲流控制微分方程和邊界條件（式（8）－式（13））的特征，該數(shù)學(xué)問題為偏微分方程組的初邊值問題。本文利用正交變換，將此問題轉(zhuǎn)化為常微分方程的初值問題，再根據(jù)矩陣微分方程的基礎(chǔ)理論便可以求出本問題的精確解。本文考慮的三維特征值問題如下

求解式（14）的三維特征值λ和特征函數(shù)，首先需分別解出x，y和z方向的3個一維特征值方程

解出特征值和特征函數(shù)為

x方向：特征值為λβ＝β2，對應(yīng)的特征函數(shù)為Eβ（x）＝eiβx，β∈（－∞，＋∞）；

y方向：特征值為λγ＝γ2，對應(yīng)的特征函數(shù)為Eγ（y）＝eiγy，γ∈（－∞，＋∞）；

z方向：特征值為λn＝（nπ）2，對應(yīng)的特征函數(shù)為En（z）＝cosnπzD，（n＝0，1，2，…）。

因此，上述三維特征值λ和特征函數(shù)分別為

特征函數(shù)系：Λ＝ ｛Eβγn（xD，yD，zD），（β，γ∈（－∞，＋∞）；n＝0，1，2，…）｝，構(gòu)成了三維空間上的完備正交系，則

式中

利用特征函數(shù)系的完備正交性便可引入相應(yīng)的正交變換。

首先引入正交變換的定義

式（21）為pD（xD，yD，zD，tD）的正交變換，記為

如果右端的廣義積分是收斂的，那么根據(jù)式（20）就可以得到逆變換公式

對式（8）—式（13）進行正交變換，最后整理化簡得到下面的常微分方程組

解得

代入逆變換公式（22），得到水平井井底壓力的無因次精確解，即全空間Ω的點源解

因此，考慮整個機動水平井對觀察點（xD，yD，zD）的壓力影響，我們則需將所求得的點源解進行擴展到整個水平井，即對點源解（式（25））進行線積分，則觀察點的壓力表達式為

5 樣板曲線

圖2為激動井和觀察井均為水平井的干擾試井井底壓力雙對數(shù)樣板曲線。從圖中可以看出，激動井忽略井筒儲集與表皮效應(yīng)的影響，其雙對數(shù)曲線主要可以分為3個流動階段：早期徑向流、中期線性流和晚期徑向流。由于激動井早期徑向流階段發(fā)生在近井筒地帶，且時間較短，因此，對應(yīng)的觀察井主要反映為2個流動階段：線性流和徑向流，出現(xiàn)的時間均晚于激動井。

圖2 水平井干擾試井理論樣板曲線Fig.2 The type curves of the well test theory for the horizontal well interference testing

6 敏感性分析

圖3與圖4為觀察井位置（x方向和y方向）對井底壓力動態(tài)樣板曲線的影響。從圖中可以看出，觀察井位置與激動井位置越近，觀察井井底壓力受激動井干擾越早出現(xiàn)，且觀察井的位置變化只影響壓力和壓力導(dǎo)數(shù)曲線的前期，不影響壓力導(dǎo)數(shù)的后期，即后期徑向流壓力導(dǎo)數(shù)值不變。

圖3 觀察井x方向距離對試井樣板曲線的影響Fig.3 Effect of the distance for observation well in xdirection on the type curves of horizontal well interference testing

圖4 觀察井y方向距離對干擾試井樣板曲線的影響Fig.4 Effect of the distance of the observation well in ydirection on the type curves of the horizontal well interference testing

圖5、圖6和圖7為觀察井和激動井不同的堆疊方式對干擾試井曲線的影響。對比3種堆疊方式，激動井的曲線沒有變化，均沒反映出干擾的特征，而圖6觀察井雙對數(shù)曲線比圖5和圖7更早出現(xiàn)線性流階段，并且壓力導(dǎo)數(shù)曲線線性流段的斜率更小。圖5與圖7表現(xiàn)出“水平并排”和“水平直線并排”這兩種堆疊方式對觀察井的雙對數(shù)曲線影響不大，其樣板曲線變化趨勢一致。

圖5 觀察井與激動井水平并排干擾試井樣板曲線Fig.5 The type curves of the horizontal wells in the interference testing at the same horizontal level

圖6 觀察井與激動井上下垂直并排干擾試井樣板曲線Fig.6 The type curves of the stacked laterals horizontal wells in the interference testing

圖7 觀察井與激動井水平直線并排干擾試井樣板曲線Fig.7 The typical curves of the horizontal opposing dual laterals horizontal wells in interference testing

7 結(jié)論

a.建立了水平井干擾試井滲流模型，并利用現(xiàn)代數(shù)學(xué)分析方法進行了求解，獲得了三維實空間下任意一點壓力的精確解，以及新型的水平井干擾試井樣板曲線。

b.獲得了水平井干擾試井中觀察井的2個典型的流動特征：線性流和徑向流，且出現(xiàn)時間均晚與激動井。

c.分析并獲得了觀察井位置以及堆疊方式對試井樣板曲線的影響規(guī)律：觀察井離激動井越近，雙對數(shù)曲線線性流出現(xiàn)的時間越早；堆疊方式不同，雙對數(shù)曲線線性流出現(xiàn)的時間以及斜率均不相同。

符號說明

Ct為綜合壓縮系數(shù)（MPa－1）；E（）為特征函數(shù)；F（）為正交變換；h為油層厚度（m）；Kx，Ky，Kz分別為x，y，z方向的滲透率；L為水平井長度（m）；LD為無因次水平井長度；p為壓力（MPa）；pi為原始地層壓力（MPa）；pw為井底壓力（MPa）；pD為無因次壓力；pwD為無因次井底壓力；~Q為點源產(chǎn)量（m3/d）；Q為總產(chǎn)量（m3/d）；QD為無因次產(chǎn)量；rw為井半徑（m）；t為時間（h）；tD為無因次時間；x，y，z為三維空間坐標(biāo)（m）；xw，yw，zw為三維空間點源坐標(biāo)（m）；δ（）為δ函數(shù)；ξx，ξy和ξz分別為δ函數(shù)在x，y，z方向上的函數(shù)變量；η為黏度（MPa·s）；q為孔隙度；λβ，λγ，λn分別為x，y，z方向特征值；λβγn為三維特征值；An為中間代換變量；τ和σ均為積分變量。

［1］鄧英爾，劉樹根，麻翠杰.井間連通性的綜合分析方法［J］.斷塊油氣田，2004，10（5）：50－53.Deng Y E，Liu S G，Ma C J.Aggregate analysis method of continuity of formation between wells［J］.Fault-Block Oil ＆ Gas Field，2004，10（5）：50－53.（In Chinese）

［2］Morris E E，Tracy G W.Determination of pore volume in a naturally fractured reservoir［C］// SPE 1185.Paper presented at Fall Meeting of the Society of Petroleum Engineers of AIME held in Colorado，Denver，1965.

［3］Jargo J R.Effect of well bore storage and wellbore damage at the active well on interference test analysis［J］.Journal of Petroleum Technology，1976：28（8）：51－58.

［4］Tongpenyai Y，Raghavan R.The effect of wellbore storage and skin on interference test data［J］.Journal of Petroleum Technology，1981，33（1）：151－160.

［5］Ogbe D O，Brigham W E.A correlation for interference testing with wellbore-storage and skin effects［J］.SPE Formation Evaluation，1989，4（3）：391－396.

［6］張秀華，賈永祿，劉平，等.考慮井儲和表皮的復(fù)合油藏干擾試井分析［J］.西南石油學(xué)院學(xué)報，2005，27（1）：57－63.Zhang X H，Jia Y L，Liu P，etal.Composite reservoir interference well testing considering wellbore storage and skin factor［J］.Journal of Southwest Petroleum Institute，2005，27（1）：57－63.（In Chinese）

［7］張德志，王子勝，姚軍.三重介質(zhì)油藏干擾試井壓力動態(tài)特征［J］.油氣地質(zhì)與采收率，2006，13（1）：74－76.Zhang D Z，Wang Z S，Yao J.Pressure behavior character of interference well test in tri－media oil reservoir［J］.Petroleum Geology and Recvery Efficiency，2006，13（1）：74－76.（In Chinese）

［8］謝林峰，甘柏松.氣藏水平井干擾試井參數(shù)敏感性研究［J］.石油天然氣學(xué)報，2012，34（9）：117－120.Xie L F，Gan B S.Parameter sensitivities analysis for horizontal well interference testing in gas reservoirs［J］.Journal of Oil and Gas Technology，2012，34（9）：117－120.（In Chinese）

［9］Kamal M M，Hegeman P S.New developments in multiple-well testing［J］.SPE formation evaluation，1988，3（1）：159－168.

［10］高承泰.串漏垂向干擾試井的特征及解釋方法［J］.石油學(xué)報，1994，15（2）：65－75.Gao C T.The features and interpretation of vertical interference well testing with serial leakage［J］.Acta Pertolei Sinica，1994，15（2）：65－75.（In Chinese）

［11］董明，席裕庚.基于遺傳算法的干擾試井解釋方法［J］.石油大學(xué)學(xué)報：自然科學(xué)版，1997，21（4）：33－36.Dong M，Xi Y G.A new method of pressure analysis for interference test based on the genetic algorithm［J］.Journal of the University of Petroleum，1997，21（4）：33－36.（In Chinese）

［12］張繼芬，衛(wèi)秀芬，賈莉卿.確定夾層穩(wěn)定性垂向干擾試井方法［J］.石油學(xué)報，1999，20（2）：62－66.Zhang J F，Wei X F，Jia L Q.A vertical in terference test method recognizing an unstable impermeable interlayer in multi-layer reservoir［J］.Acta Pertolei Sinica，1999，20（2）：62－66.（In Chinese）

［13］張同義，楊永智，杜鵑.多井試井資料在油田開發(fā)中的應(yīng)用［J］.測井技術(shù)，2004，28（B02）：65－68.Zhang T Y，Yang Y Z，Du J.Applications of multiwell test data in oil field developments［J］.Well Logging Technology，2004，28（B02）：65－68.（In Chinese）

［14］劉啟國，陳彥麗，張烈輝，等.不關(guān)井干擾試井模型和解釋方法研究［J］.油氣井測試，2006，15（1）：10－12.Liu Q G，Chen Y L，Zhang L H，etal.Study for disturbing test model with the well opening and its interpreting method［J］.Well Testing，2006，15（1）：10－12.（In Chinese）

［15］王本成，石國新，路建國，等.考慮二次梯度的非線性雙重介質(zhì)滲流模型［J］.成都理工大學(xué)學(xué)報：自然科學(xué)版，2011，38（2）：140－146 Wang B C，Shi G X，Lu J G，etal.Study on quadratic pressure gradient terms of the nonlinear filtering flow model in a dual medium reservoir［J］.Journal of Chengdu University of Technology（Science ＆ Technology Edition），2011，38（2）：140－146.（In Chinese）