劉海龍

(中國(guó)石油勘探開(kāi)發(fā)研究院，北京　海淀　100083)

?

低滲透矩形油藏不穩(wěn)定滲流壓力動(dòng)態(tài)分析

劉海龍

(中國(guó)石油勘探開(kāi)發(fā)研究院，北京海淀100083)

摘要：目前，低滲透矩形油藏地層滲流壓力求解的研究比較成熟，但是對(duì)于滲流壓力的動(dòng)態(tài)分析及滲流壓力影響因素分析的研究尚存在一些不足。針對(duì)低滲透矩形油藏難動(dòng)用的難題，基于兩參數(shù)連續(xù)模型，運(yùn)用穩(wěn)態(tài)逐次替換法，分別建立了定壓和定產(chǎn)兩種生產(chǎn)形式下的壓力動(dòng)態(tài)分析模型，并從地層物性和導(dǎo)壓系數(shù)兩個(gè)層面進(jìn)行了相關(guān)因素分析。分析認(rèn)為：無(wú)論在定產(chǎn)模型還是在定壓模型中，壓力傳播距離與時(shí)間均為二次三項(xiàng)式關(guān)系，且壓力傳播快慢與產(chǎn)量、滲透率呈正相關(guān)，與流體黏度、井底流壓、啟動(dòng)壓力梯度呈負(fù)相關(guān)。實(shí)例計(jì)算表明：所建立的模型可靠，為矩形油藏注水開(kāi)發(fā)，預(yù)計(jì)注水見(jiàn)效時(shí)間，提供一定的理論依據(jù)。

關(guān)鍵詞：矩形油藏；低滲透油藏；不穩(wěn)定滲流

隨著我國(guó)國(guó)民經(jīng)濟(jì)穩(wěn)定增長(zhǎng)，石油需求不斷增加，我國(guó)已經(jīng)成為僅次于美國(guó)的世界第二大石油消費(fèi)國(guó)[1]。但是由于中、高滲透油藏后備石油儲(chǔ)量不足，每年新增探明儲(chǔ)量中，低滲透油藏比例越來(lái)越大，因此加速開(kāi)發(fā)低滲透油藏，意義重大[2]。

目前，我國(guó)的老油田開(kāi)采面臨兩大難問(wèn)題：高含水和低滲透。針對(duì)于低滲透的研究，前人已經(jīng)做出了很多研究，且取得了重大進(jìn)展：如馬爾哈辛提出低滲透油藏存在啟動(dòng)壓力梯度，并從微觀(guān)角度解釋了其產(chǎn)生的機(jī)理[3]；Pascal 等人首次應(yīng)用數(shù)值積分方法——有限差分法[4-5]，在考慮啟動(dòng)壓力梯度條件下，求解流體滲流壓力分布。

宋付權(quán)等人結(jié)合室內(nèi)實(shí)驗(yàn)，建立低滲透油藏巖心一開(kāi)一關(guān)滲流數(shù)學(xué)模型，對(duì)一維壓力傳播邊界進(jìn)行數(shù)值求解[6]；李凡華等人考慮啟動(dòng)壓力梯度，建立了無(wú)限大和有界低滲透油藏不穩(wěn)定滲流試井分析模型[7]；鄧英爾等人在實(shí)驗(yàn)的基礎(chǔ)上，首次提出非達(dá)西滲流連續(xù)函數(shù)模型，并建立了非線(xiàn)性穩(wěn)態(tài)滲流的壓力和產(chǎn)量公式[8]；劉鵬程等人結(jié)合室內(nèi)實(shí)驗(yàn)和油田實(shí)際生產(chǎn)，提出壓敏分段變化理論，并在此基礎(chǔ)上給出了油氣井單井產(chǎn)能公式和產(chǎn)量預(yù)測(cè)方法[9-10]；郝明強(qiáng)等人考慮啟動(dòng)壓力梯度和壓敏效應(yīng)，通過(guò)擬壓力變換，推導(dǎo)了平面徑向流壓力分布公式[11]。朱圣舉采用穩(wěn)態(tài)逐次替換法，給出了不同滲流方式下，低滲透油藏孔隙中彈性流體低速不穩(wěn)定滲流壓力的傳播規(guī)律[12]。楊清立提出的兩參數(shù)非線(xiàn)性連續(xù)模型既反映出流體在低滲介質(zhì)中滲流時(shí)存在最小啟動(dòng)壓力梯度的現(xiàn)象， 又可以很好地

描述非線(xiàn)性段特征[13]。姜瑞忠基于兩參數(shù)模型，采用數(shù)值模擬方法，建立了兩維兩相非線(xiàn)性滲流數(shù)值模型[14]。但是目前針對(duì)于壓力傳播規(guī)律的研究還很欠缺，且大多只考慮擬啟動(dòng)壓力梯度，認(rèn)為低滲透油藏滲流的啟動(dòng)壓力梯度與油藏本身的無(wú)關(guān)，其無(wú)法真實(shí)反映地下流體的滲流特征。或是基于一維徑向滲流[15]，或是研究定產(chǎn)量生產(chǎn)[16]，然而實(shí)際生產(chǎn)中，存在定井底流壓生產(chǎn)情況，且對(duì)于生產(chǎn)中的大型線(xiàn)性排狀注水，或是滲流物理模擬實(shí)驗(yàn)，流體會(huì)發(fā)生一維單向流動(dòng)，因而有必要對(duì)低滲透油藏不穩(wěn)定滲流壓力進(jìn)行動(dòng)態(tài)分析。

1物理模型

圖1為一平面帶狀等厚低滲油藏，地層流體微可壓縮，粘度為μ，油藏平均厚度為h，寬度為w，長(zhǎng)度為L(zhǎng)，地層外邊界壓力為pe，井底流壓為pwf，油藏中一口井以定產(chǎn)量或是定井底流壓方法生產(chǎn)，儲(chǔ)層中的流體以不穩(wěn)定滲流形式流向井筒。

圖1　低滲透油藏一維滲流模型

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

目前非線(xiàn)性滲流的數(shù)學(xué)模型描述方法很多[17-19]，如表1所示。

表1　非線(xiàn)性滲流數(shù)學(xué)模型

由表1知，本文選取兩參數(shù)連續(xù)模型是合理的，不考慮壓力敏感效應(yīng)的低滲油藏非線(xiàn)性滲流的兩參數(shù)模型表達(dá)式為

(1)

式中：v為滲流速度，am/s；k為絕對(duì)滲透率，mD；μ為流體黏度，cp；p為流體壓力，atm；x為滲流介質(zhì)長(zhǎng)度；a，b——兩參數(shù)連續(xù)模型系數(shù)，可通過(guò)室內(nèi)試驗(yàn)測(cè)定獲得；C為單位換算系數(shù)，當(dāng)所取單位為達(dá)西制單位時(shí)，C為1，當(dāng)所取單位為SI制單位時(shí)，C為86.4，本文取達(dá)西制單位進(jìn)行研究。

任意位置處的滲流速度為

(2)

式中：w為油藏寬度，cm；h為油藏有效厚度，cm；q為地層流體單位時(shí)間內(nèi)的流量，cm3/s；A為滲流面積，cm2。

聯(lián)立(1)、(2)式，可得

(3)

求解(3)得

當(dāng)b≠0時(shí)，

(4)

利用文獻(xiàn)12中的穩(wěn)態(tài)逐次替換法，可得

p-pe=ε(q)[x(t)-L]

(5)

p-pwf=ε(q)[x(t)-rw]

(6)

式中：pe為外邊界壓力，atm；pwf為井底壓力，atm；p為油藏x處的流體壓力，atm；x(t)為壓力傳播距離，cm；L為油藏長(zhǎng)度，cm；rw為井筒半徑，cm；t為時(shí)間，s。

ε(q)=

(7)

2.1定產(chǎn)生產(chǎn)模型

以地下定流量q生產(chǎn)，則

Np=ρ0qt

(8)

式中：ρ0為流體密度，t/cm3；Np為t時(shí)間內(nèi)的累積地下流量，t。

考慮井筒存儲(chǔ)效益，由質(zhì)量守恒定律得

whrw[(ρφ)e-(ρφ)w]

(9)

由文獻(xiàn)20可得

ρφ=ρ0φ0[1+ct(p-p0)]

(10)

式中：ρ0為原始地層壓力的流體密度，t/cm3；φ0、φ為原始地層壓力、地層壓力下的流體孔隙度，小數(shù)；p0為原始地層壓力，atm；ct為綜合壓縮系數(shù)，atm-1。

由(6)-(10)式可得

(11)

式(11)即為定產(chǎn)條件下，時(shí)間與壓力傳播距離關(guān)系式。將式(11)帶入式(5)，即可得出定產(chǎn)條件下的壓力傳播規(guī)律。

(12)

由式(12)可知，壓力傳播與時(shí)間成冪指數(shù)關(guān)系，這與中高滲油藏的一維壓力傳播與時(shí)間為線(xiàn)性關(guān)系不同。原因在于低滲透油藏本身孔喉小，在壓力敏感效應(yīng)下，孔喉直徑縮短為原始的70%[21]。且存在啟動(dòng)壓力梯度，尤其是儲(chǔ)層能力不足時(shí)(壓力系數(shù)太小，小于1時(shí))，使得孔隙流體壓力不足以傳播到井筒，啟動(dòng)壓力梯度影響更為明顯。低滲透油藏在兩種機(jī)制作用下，壓力傳播速度變慢，傳播規(guī)律發(fā)生變化。

對(duì)于實(shí)際油田開(kāi)發(fā)而言，要提高壓力傳播速度，必須借助外來(lái)能量來(lái)補(bǔ)充地層能量的不足，如注水或注氣，從而提高儲(chǔ)層孔隙流體壓力，減小或消除啟動(dòng)壓力梯度的影響，進(jìn)而使得油田正常開(kāi)發(fā)。采用壓裂改造低滲透儲(chǔ)層，在儲(chǔ)層中建立“流動(dòng)網(wǎng)絡(luò)”，聯(lián)通儲(chǔ)層更多的滲流通道，增大儲(chǔ)層暴露的滲流面積，加上外部注水或注氣，提前補(bǔ)充地層能量，增加儲(chǔ)層壓力，減小儲(chǔ)層流體滲流阻力，使得低滲透儲(chǔ)層啟動(dòng)壓力梯度減小或“消失”。

1) 模型應(yīng)用。取文獻(xiàn)16的參數(shù)：低滲透油藏地層原始?jí)毫?0 MPa，孔隙度為0.12，滲透率為1.2 mD，流體黏度為0.256 mPa.s，啟動(dòng)壓力梯度為0.02 MPa/m，流量為1.0 m3/d，矩形油藏長(zhǎng)、寬、高分別為100 m、20 m、5 m，a、b取值為0.908、10.5。將上述數(shù)據(jù)代入式(12)，并與表1中的擬壓力梯度模型和達(dá)西模型作對(duì)比，結(jié)果如圖2所示。

t/d1. 擬啟動(dòng)壓力梯度模型；2. 擬線(xiàn)性模型；3. 達(dá)西模型圖2　計(jì)算結(jié)果對(duì)比圖

圖2顯示了不同滲流模型下壓力傳播的范圍和深度，三種模型的趨勢(shì)大致相同，在井筒附近壓降梯度最快，遠(yuǎn)離井筒處的壓降梯度慢慢變緩。當(dāng)壓力傳播距離相同時(shí)，非線(xiàn)性模型所用時(shí)間比擬啟動(dòng)壓力梯度模型的小，比達(dá)西模型的大，這是因?yàn)樵诘蜐B透油藏中，當(dāng)壓力達(dá)到最小啟動(dòng)壓力梯度時(shí)，地層中的流體就開(kāi)始流動(dòng)，而啟動(dòng)壓力梯度模型中，只有當(dāng)壓力克服擬啟動(dòng)壓力梯度時(shí)，流體才可以流動(dòng)，故在傳播相同距離時(shí)，所用時(shí)間要長(zhǎng)。相反，由于達(dá)西模型不考慮啟動(dòng)壓力梯度的影響，因而在傳播相同距離時(shí)，用的時(shí)間最少。各模型計(jì)算結(jié)果相差很大，低滲透、特低滲透油藏由于具有明顯的非線(xiàn)性滲流特征，采用線(xiàn)性模型分析問(wèn)題，必然會(huì)導(dǎo)致較大的誤差，因此建議使用非線(xiàn)性模型進(jìn)行相關(guān)問(wèn)題的分析和計(jì)算。

隨著油田實(shí)際生產(chǎn)的進(jìn)行，壓力傳播的動(dòng)邊界逐漸往外擴(kuò)大，分別取x=20 m，40 m，60 m，80 m，100 m，將數(shù)據(jù)帶入式(12)，觀(guān)察其壓力隨時(shí)間的變化規(guī)律，如圖(3)所示。再取t=60 d，70 d，80 d，90 d，觀(guān)察各時(shí)間地層中的壓力分布規(guī)律，如圖(4)所示。

t/d1. x=20 m；2. x=30 m；3. x=60 m；4. x=80 m；5. x=100 m圖3　距井底不同位置壓力隨時(shí)間變化規(guī)律

x/m1. t=60 d；2. t=70 d；3. t=80 d；4. t=90 d圖4　不同時(shí)刻地層各點(diǎn)壓力分布圖

由圖3知：隨著油井的采油，壓力逐漸向邊界擴(kuò)展，當(dāng)動(dòng)邊界擴(kuò)展到20 m，40 m，60 m，80 m，100 m時(shí)，壓力降落漏斗急劇增大，尤其在井筒附近表現(xiàn)最為突出，地層壓力不斷降低。圖4顯示：不同時(shí)間，壓力波及的范圍不一樣，地下壓力存在壓力波及區(qū)和未波及區(qū)，井筒附近的壓力變化最為明顯，越偏離井筒，壓力降落速度偏低。

2) 單因素分析。壓力傳播快慢受控于地層物性、溫壓力系統(tǒng)等因素，考慮實(shí)際生產(chǎn)需要，本文主要研究產(chǎn)量、滲透率、黏度等對(duì)壓力傳播快慢的影響。分別取k=0.1 mD，0.25 mD，0.5 mD，10 mD， 50 mD， 分別取μ=0.5 MPa·s， 3 MPa·s，10 MPa·s，20 MPa·s分別取Q=1 m3/d，3 m3/d，5 m3/d，7 m3/d，9 m3/d，觀(guān)察壓力與時(shí)間的關(guān)系，如圖5～圖7所示。

t/d1. k=50.0 mD；2. k=10.0 mD；3. k=0.50 mD；4. k=0.25 mD；5. k=0.10 mD圖5　滲透率對(duì)壓力傳播距離的影響

圖5表明：隨著滲透率的增大，壓力傳播的越來(lái)越快(動(dòng)邊界往外擴(kuò)散速度增大)，因?yàn)闈B透率增大，地層流體流動(dòng)能力增強(qiáng)，相應(yīng)的地層阻力減小，壓力波更容易向外傳播，進(jìn)而動(dòng)邊界向往移動(dòng)速度增大。

t/d1. μ=0.50 MPa·s；2. μ=3.00 MPa·s；3. μ=10.0 MPa·s；4. μ=20.0 MPa·s；圖6　黏度對(duì)壓力傳播距離的影響

圖6顯示：隨著黏度的增大，壓力傳播的越來(lái)越慢(動(dòng)邊界往外擴(kuò)散速度減小)，這是因?yàn)榈貙恿黧w黏度越大，流動(dòng)阻力就越大，需要更多的時(shí)間克服阻力流動(dòng)，壓力傳播速度變慢，動(dòng)邊界往外擴(kuò)散速度降低。

t/d1. Q=9 m3/d；2. Q=7 m3/d；3. Q=5 m3/d；4. Q=3 m3/d；5. Q=1 m3/d圖7　日產(chǎn)油量對(duì)壓力傳播距離的影響

由圖7知，對(duì)于同一時(shí)間，隨著采油量的增加，壓力傳播的越遠(yuǎn)，即壓力傳播的越快。因?yàn)樵谙嗤瑫r(shí)間內(nèi)，增加采油量，就必須提高地層流體的流速，則就需要提高壓力梯度，但是對(duì)于定邊界壓力的低滲透油藏，相同距離的壓力梯度是一樣的，為達(dá)到壓力梯度在時(shí)間上的不一致，就必須使壓力傳播的更快，才能滿(mǎn)足油井產(chǎn)油量增加。反過(guò)來(lái)，采油量的增大，生產(chǎn)壓差也會(huì)在一定程度上增加(以增大生產(chǎn)壓差提高油井產(chǎn)量)，需要壓力擴(kuò)散到更大的空間范圍，增加壓力激動(dòng)區(qū)的面積，即表現(xiàn)出，壓力傳播相同時(shí)間時(shí)，隨著油井產(chǎn)量的提高，壓力傳播的越快、越遠(yuǎn)。

2.2定壓生產(chǎn)模型

t時(shí)間內(nèi)累積產(chǎn)量為

Np=ρ0∫t0BQdt

(13)

由式(6)、(7)解得Q為

(14)

由式(5)、(6)相比得：

(15)

聯(lián)立(9)、(14)、(15)、(16)式并忽略無(wú)窮小項(xiàng)得

(16)

求解(16)式得

t=F[Cx(t)2+Dx(t)+F]

(17)

(18)

(19)

(20)

(21)

式(17)定壓條件下，時(shí)間與壓力傳播距離關(guān)系式。與定產(chǎn)模型呈現(xiàn)類(lèi)似的規(guī)律。式(17)表明：定井底流壓條件下，矩形油藏壓力開(kāi)始降落，呈“漏斗狀”向外擴(kuò)展，只是壓力在近井壁處近似為一定值，當(dāng)壓力傳播到邊界時(shí)，由于邊界壓力一定，這時(shí)為保持這種定井底流壓的情況，必須對(duì)地層補(bǔ)充能量，此時(shí)油井的產(chǎn)量主要是兩部分：一部分是邊界進(jìn)入地層的流體，q1，并且流量逐漸增加；另一部分是邊界內(nèi)部地層依靠彈性能量膨脹產(chǎn)生的流體q2，但是卻逐漸減少。當(dāng)油井產(chǎn)量為q1時(shí)，由不穩(wěn)態(tài)滲流轉(zhuǎn)變?yōu)榉€(wěn)態(tài)滲流。

1) 模型應(yīng)用。由式(17)知：時(shí)間與壓力傳播距離為二次三項(xiàng)式，若井底流壓為7.5 MPa時(shí)，數(shù)據(jù)相關(guān)數(shù)據(jù)代入式(28)，并依此取啟動(dòng)壓力梯度為0.02 MPa/m，0.07 MPa/m，0.20 MPa/m和滲透率為0.25 mD，0.50 mD，1.0 mD，2.0 mD得到定壓生產(chǎn)10 d的壓力傳播規(guī)律，如圖8所示。

t/d1. 啟動(dòng)壓力梯度為0.02 MPa/m；2. 啟動(dòng)壓力梯度為0.02 MPa/m；3. 啟動(dòng)壓力梯度為0.02 MPa/m

t/d1. k=2.00 mD；2. k=1.00 mD；3. k=0.50 mD；4. k=0.25 mD

由圖8知：?jiǎn)?dòng)壓力梯度和滲透率均能在一定程度上減緩壓力傳播速度。啟動(dòng)壓力梯度大，地層阻力大；滲透率大，地層阻力小。兩種因素相互制約，由于啟動(dòng)壓力梯度是低滲透油藏儲(chǔ)層本身的屬性，通過(guò)外部很難改變，但是地層滲透率卻是可以改造的。目前工藝上主要采用壓裂、酸化進(jìn)行，低滲透油藏更多采用壓裂進(jìn)行增產(chǎn)。

2) 單因素分析。同理，如2.1.2所述，本文主要研究井底流壓、滲透率等對(duì)壓力傳播快慢的影響。分別取pwf=12 MPa，14 MPa，16 MPa，18 MPa觀(guān)察壓力與時(shí)間的關(guān)系，如圖9所示。

t/d1. Pwf=18MPa mD；2. Pwf=16MPa mD；3. Pwf=14MPa mD；4. Pwf=12MPa mD圖9　井底流壓對(duì)壓力傳播規(guī)律的影響

由圖9知：圖形斜率逐漸變小，壓力傳播速度減小，直到壓力波及到邊界。主要是因?yàn)閭鞑ゾ嚯x越遠(yuǎn)，耗散的能量越大，若得不到外界能量的及時(shí)補(bǔ)給，則地層流體無(wú)法克服地層阻力，從而使得油井產(chǎn)量下降，這也是諸多低滲透油藏采用注水開(kāi)發(fā)，補(bǔ)給地層能量的原因。當(dāng)傳播距離不大時(shí)，井底流壓的影響較小，隨著距離的增大，井底流壓的影響增大。

減小井底流壓，放大生產(chǎn)壓差，可以有效提高油井產(chǎn)量，但是并非井底流壓越小越好，而是油井以最合理的井底流壓生產(chǎn)，可實(shí)現(xiàn)油井高效、穩(wěn)定開(kāi)發(fā)。(14)式解出x(t)，并代入(17)式，求解得合理井底壓力為

(22)

(23)

將文獻(xiàn)17數(shù)據(jù)代入式(22)得合理井底壓力為10.31 MPa，對(duì)應(yīng)的最大產(chǎn)油量3.26 m3/d。

3實(shí)例分析

我國(guó)部分地區(qū)由于長(zhǎng)期受到地質(zhì)活動(dòng)的影響，導(dǎo)致這些地區(qū)出現(xiàn)了很多的復(fù)雜斷塊油藏，蘇北的臺(tái)興油田是一個(gè)典型的狹長(zhǎng)型低滲透復(fù)雜斷塊油田[22-23]。該油田由于受到多級(jí)斷層相互作用，使得臺(tái)興油田發(fā)育Ⅲ-Ⅴ級(jí)斷層11條，斷層將整個(gè)臺(tái)興油田分為11個(gè)含油斷塊。每個(gè)含油斷塊都有獨(dú)立的溫度、壓力系統(tǒng)，每個(gè)小油藏形狀近似矩形，寬度較小，長(zhǎng)度較長(zhǎng)，符合本文模型的應(yīng)用條件。下面以文獻(xiàn)23中的區(qū)塊為例，進(jìn)行實(shí)例分析。

目前數(shù)值模型的理論是建立在經(jīng)典達(dá)西滲流基礎(chǔ)上的，而Eclipse2010中的E300模塊是針對(duì)低滲透油藏，考慮了啟動(dòng)壓力梯度，因此將本文的定產(chǎn)模型計(jì)算的某時(shí)刻的井底流壓(定產(chǎn)解)與數(shù)值模擬結(jié)果(數(shù)值解)對(duì)比，如表2所示。由表2可知，相對(duì)誤差基本控制在8%內(nèi)，說(shuō)明定流壓模型是可靠的。

表2　定產(chǎn)模型實(shí)例計(jì)算對(duì)比表

由于定流壓模型中，設(shè)定井底流壓為常數(shù)，因此將式(17)帶入式(14)，求得不同時(shí)刻的產(chǎn)量(定壓解)，并與數(shù)值模擬(數(shù)值解)對(duì)比，如表3所示。由表3可知，相對(duì)誤差控制在7.7%內(nèi)，則定流壓模型是可靠的。

表3　定流壓模型實(shí)例計(jì)算對(duì)比表

4結(jié)論

1) 定產(chǎn)模型中，時(shí)間與壓力傳播距離為二次三項(xiàng)式關(guān)系，且在壓力傳播速度上，兩參數(shù)連續(xù)模型比達(dá)西模型慢，比擬啟動(dòng)壓力梯度模型要快，與流體黏度呈負(fù)相關(guān)，與產(chǎn)量、滲透率呈正相關(guān)；

2) 定壓模型中，時(shí)間與壓力傳播距離亦為二次三項(xiàng)式關(guān)系，但壓力傳播速度與滲透率呈正相關(guān)，與井底流壓、啟動(dòng)壓力梯度呈負(fù)相關(guān)。

3) 實(shí)例計(jì)算表面，建立的定產(chǎn)模型、定壓模型是可靠的。

參考文獻(xiàn)：

[1]胡文瑞.中國(guó)低滲透油氣的現(xiàn)狀與未來(lái)[J].中國(guó)工程科學(xué)，2009，11(8)：29-37.

[2]王光付，廖榮鳳，李江龍，等.中國(guó)石化低滲透油藏開(kāi)發(fā)狀況及前景[J].油氣地質(zhì)與采收率，2007，14(3)：84-89.

[3]MILLER R J.Threshold gradient for water flow in clay system， Proc[J].Soil.Sci.Soc.Am.，1963，27：606-609.

[4]PASCAL.F.Consolidation with Threshold Gradient[J].Inter.J. for Numerical and Analytical Methods in Geomechanics.1980，5：247-261.

[5]ALVARO，P，F(xiàn)ARUK C.Mddification of Darcy's law for the threshold pressure gradient，[J].Petrol.Sci.Eng.1999，22(4)：237-240.

[6]宋付權(quán)，劉慈群，李凡華.低滲透介質(zhì)含啟動(dòng)壓力梯度一維瞬時(shí)壓力分析[J].應(yīng)用數(shù)學(xué)和力學(xué)，1999，20(1)：25-33.

[7]李凡華，劉慈群.含啟動(dòng)壓力梯度的不定常滲流的壓力動(dòng)態(tài)分析[J].油氣井測(cè)試，1997，6(1)：1-4.

[8]鄧英爾，劉慈群.低滲油藏非線(xiàn)性滲流規(guī)律數(shù)學(xué)模型及其應(yīng)用[J].石油學(xué)報(bào)，2001，22(4)：72-77.

[9]劉鵬程，王曉冬，李素珍，等.地層壓敏對(duì)低滲透氣井產(chǎn)能影響研究 [J].西南石油學(xué)院學(xué)報(bào)，2004，26(5)：37-41.

[10]劉鵬程，王曉冬，萬(wàn)玉金.三區(qū)復(fù)合油藏有限導(dǎo)流垂直裂縫井井底壓力動(dòng)態(tài)分析[J].油氣井測(cè)試，2004，13(1)：4-7.

[11]郝明強(qiáng)，胡永樂(lè)，李凡華.微裂縫性特低滲透油藏產(chǎn)能研究[J].石油天然氣學(xué)報(bào)，2009，31(2)：100-104.

[12]朱圣舉.低滲透油藏的壓力波傳播規(guī)律[J].新疆石油地質(zhì)，2007， 8(1)：85-87.

[13]楊清立，楊正明，王一飛，等.特低滲透油藏滲流理論研究[J].鉆采工藝，2007，30(6)：52-54.

[14]XU J.Non-Darcy flow numerical simulation for low permeability reservoirs，SPE 154890.

[15]史瑞娜，王曉冬.低滲透地層壓力擾動(dòng)傳播規(guī)律[J].特種油氣藏，2011， 18(4)：80-82.

[16]黃爽英，陳祖華.引入啟動(dòng)壓力梯度計(jì)算低滲透砂巖油藏注水見(jiàn)效時(shí)間[J]，河南石油，2011，15(5)：22-24.

[17]BAOQUAN ZENG，LINSONG CHENG，CHUNLAN LI.Low velocity non-linear flow in low permeability reservoir[J].Petrol.Sci.Eng.2011，80(1)：1-6.

[18]李松泉，程林松，李秀生，等.特低滲透油藏非線(xiàn)性滲流模型[J].石油勘探與開(kāi)發(fā)，2008，35(5)：606-612.

[19]XU J.Transient pressure behavior for dual porosity low peameability reservoir based on modified darcy.s equation[C]//SPE 153480.

[20]王曉冬.滲流力學(xué)基礎(chǔ)[M].北京：中國(guó)地質(zhì)大學(xué)出版社，2006：11-16.

[21]SHAOUL J ，ZELM L Damage mechanisms in unconventional gas well stimulation-A new look at an old problem[C]// SPE 142479.

[22]劉睿，姜漢橋，陳民峰，等.江蘇油田復(fù)雜小斷塊油藏分類(lèi)[J].新疆石油地質(zhì)，2009，30(6)：680-682.

[23]LIU RUI， JIANG HANQIAO， CHEN MINFENG， et，al. The classification of small complex fault block reservoirs， taking Jiangsu Oilfield for example [J]. Xinjiang Petroleum Geology， 2009，30 (6)： 680-682.

[24]黃偉.復(fù)雜小斷塊水平井整體開(kāi)發(fā)合理井網(wǎng)研究[J].石油天然氣學(xué)報(bào)，2009，31(3)：105-111.

[25]HUANG WEI. Study on reasonable well pattern in horizontal well integrated development of complex small fault block[J]. Journal of Oil and Gas Technology，2009，31(3)：105-111.

(責(zé)任編輯：李麗)

Dynamic Analysis of Unsteady Seepage Pressure in Low Permeability Rectangular Reservoirs

LIU Hai-long

(Research Institute of Petroleum Exploration and Development, PetroChina, Beijing 100083, China)

Abstract:At present, the research on solution of seepage pressure in low permeability rectangular reservoir formation is relatively mature. However, there are still some deficiencies in the study of dynamic analysis of the seepage pressure and the influencing factors on seepage pressure. Aiming at the difficult problem of low permeability rectangular reservoir, based on the two-parameter continuous model, the dynamic analysis models for constant pressure or fixed production were established separately by using the steady state successive substitution method, and the related factors analysis was carried out from the aspects of the formation characteristics and the pressure coefficient.The analysis showed that whether in the fixed production model or in the constant pressure model, the relationship between the pressure propagation distance and time is quadratic trinomial, and the speed of pressure propagation is positively correlated with the yield and the permeability, but negatively correlated with the viscosity of the fluid, the bottom hole flowing pressure and the starting pressure gradient.Example calculation showed that the models are reliable, which can provide some theoretical basis for the rectangular reservoir development by water injection and for the calculation of expected effective time of water injection.

Key words:rectangular reservoirs; low permeability reservoirs; unsteady seepage

收稿日期：2015-03-09

作者簡(jiǎn)介：劉海龍(1989-)，男，四川廣安人，在讀碩士，研究方向：油氣田開(kāi)發(fā)。

中圖分類(lèi)號(hào)：TE121.1

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1672-1098(2016)01-0075-08