低滲透矩形油藏不穩定滲流壓力動態分析

劉海龍

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

關鍵詞：矩形油藏；低滲透油藏；不穩定滲流

中圖分類號：TE1211 文獻標志碼：A文章編號：1672-1098（2016）01-0075-08

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

隨著我國國民經濟穩定增長，石油需求不斷增加，我國已經成為僅次于美國的世界第二大石油消費國[1]。但是由于中、高滲透油藏后備石油儲量不足，每年新增探明儲量中，低滲透油藏比例越來越大，因此加速開發低滲透油藏，意義重大[2]。

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

宋付權等人結合室內實驗，建立低滲透油藏巖心一開一關滲流數學模型，對一維壓力傳播邊界進行數值求解[6]；李凡華等人考慮啟動壓力梯度，建立了無限大和有界低滲透油藏不穩定滲流試井分析模型[7]；鄧英爾等人在實驗的基礎上，首次提出非達西滲流連續函數模型，并建立了非線性穩態滲流的壓力和產量公式[8]；劉鵬程等人結合室內實驗和油田實際生產，提出壓敏分段變化理論，并在此基礎上給出了油氣井單井產能公式和產量預測方法[9-10]；郝明強等人考慮啟動壓力梯度和壓敏效應，通過擬壓力變換，推導了平面徑向流壓力分布公式[11]。朱圣舉采用穩態逐次替換法，給出了不同滲流方式下，低滲透油藏孔隙中彈性流體低速不穩定滲流壓力的傳播規律[12]。楊清立提出的兩參數非線性連續模型既反映出流體在低滲介質中滲流時存在最小啟動壓力梯度的現象， 又可以很好地描述非線性段特征[13]。姜瑞忠基于兩參數模型，采用數值模擬方法，建立了兩維兩相非線性滲流數值模型[14]。但是目前針對于壓力傳播規律的研究還很欠缺，且大多只考慮擬啟動壓力梯度，認為低滲透油藏滲流的啟動壓力梯度與油藏本身的無關，其無法真實反映地下流體的滲流特征。或是基于一維徑向滲流[15]，或是研究定產量生產[16]，然而實際生產中，存在定井底流壓生產情況，且對于生產中的大型線性排狀注水，或是滲流物理模擬實驗，流體會發生一維單向流動，因而有必要對低滲透油藏不穩定滲流壓力進行動態分析。

1物理模型

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

目前非線性滲流的數學模型描述方法很多[17-19]，如表1所示。

dpdx=-（a-1-bN）+（a-1-bN）2+4abN2b=

-（a-1-bN）+（a-1+bN）2+4bN2b（4）

對（4）式進行討論，當a→∞，b→0時，（4）式可簡化為dpdx=N=μqwkh，此解與達西定律對應，驗證了解析解是正確的。

利用文獻12中的穩態逐次替換法，可得

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

對于實際油田開發而言，要提高壓力傳播速度，必須借助外來能量來補充地層能量的不足，如注水或注氣，從而提高儲層孔隙流體壓力，減小或消除啟動壓力梯度的影響，進而使得油田正常開發。采用壓裂改造低滲透儲層，在儲層中建立“流動網絡”，聯通儲層更多的滲流通道，增大儲層暴露的滲流面積，加上外部注水或注氣，提前補充地層能量，增加儲層壓力，減小儲層流體滲流阻力，使得低滲透儲層啟動壓力梯度減小或“消失”。

1） 模型應用。取文獻16的參數：低滲透油藏地層原始壓力為30 MPa，孔隙度為012，滲透率為12 mD，流體黏度為0256 mPa.s，啟動壓力梯度為002 MPa/m，流量為10 m3/d，矩形油藏長、寬、高分別為100 m、20 m、5 m，a、b取值為0908、105。將上述數據代入式（12），并與表1中的擬壓力梯度模型和達西模型作對比，結果如圖2所示。

t/d

1. 擬啟動壓力梯度模型；2. 擬線性模型；3. 達西模型

圖2顯示了不同滲流模型下壓力傳播的范圍和深度，三種模型的趨勢大致相同，在井筒附近壓降梯度最快，遠離井筒處的壓降梯度慢慢變緩。當壓力傳播距離相同時，非線性模型所用時間比擬啟動壓力梯度模型的小，比達西模型的大，這是因為在低滲透油藏中，當壓力達到最小啟動壓力梯度時，地層中的流體就開始流動，而啟動壓力梯度模型中，只有當壓力克服擬啟動壓力梯度時，流體才可以流動，故在傳播相同距離時，所用時間要長。相反，由于達西模型不考慮啟動壓力梯度的影響，因而在傳播相同距離時，用的時間最少。各模型計算結果相差很大，低滲透、特低滲透油藏由于具有明顯的非線性滲流特征，采用線性模型分析問題，必然會導致較大的誤差，因此建議使用非線性模型進行相關問題的分析和計算。

隨著油田實際生產的進行，壓力傳播的動邊界逐漸往外擴大，分別取x=20 m，40 m，60 m，80 m，100 m，將數據帶入式（12），觀察其壓力隨時間的變化規律，如圖（3）所示。再取t=60 d，70 d，80 d，90 d，觀察各時間地層中的壓力分布規律，如圖（4）所示。

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

2） 單因素分析。壓力傳播快慢受控于地層物性、溫壓力系統等因素，考慮實際生產需要，本文主要研究產量、滲透率、黏度等對壓力傳播快慢的影響。分別取圖5表明：隨著滲透率的增大，壓力傳播的越來越快（動邊界往外擴散速度增大），因為滲透率增大，地層流體流動能力增強，相應的地層阻力減小，壓力波更容易向外傳播，進而動邊界向往移動速度增大。

圖6顯示：隨著黏度的增大，壓力傳播的越來越慢（動邊界往外擴散速度減小），這是因為地層流體黏度越大，流動阻力就越大，需要更多的時間克服阻力流動，壓力傳播速度變慢，動邊界往外擴散速度降低。

由圖7知，對于同一時間，隨著采油量的增加，壓力傳播的越遠，即壓力傳播的越快。因為在相同時間內，增加采油量，就必須提高地層流體的流速，則就需要提高壓力梯度，但是對于定邊界壓力的低滲透油藏，相同距離的壓力梯度是一樣的，為達到壓力梯度在時間上的不一致，就必須使壓力傳播的更快，才能滿足油井產油量增加。反過來，采油量的增大，生產壓差也會在一定程度上增加（以增大生產壓差提高油井產量），需要壓力擴散到更大的空間范圍，增加壓力激動區的面積，即表現出，壓力傳播相同時間時，隨著油井產量的提高，壓力傳播的越快、越遠。

式（17）定壓條件下，時間與壓力傳播距離關系式。與定產模型呈現類似的規律。式（17）表明：定井底流壓條件下，矩形油藏壓力開始降落，呈“漏斗狀”向外擴展，只是壓力在近井壁處近似為一定值，當壓力傳播到邊界時，由于邊界壓力一定，這時為保持這種定井底流壓的情況，必須對地層補充能量，此時油井的產量主要是兩部分：一部分是邊界進入地層的流體，q1，并且流量逐漸增加；另一部分是邊界內部地層依靠彈性能量膨脹產生的流體q2，但是卻逐漸減少。當油井產量為q1時，由不穩態滲流轉變為穩態滲流。

1） 模型應用。由式（17）知：時間與壓力傳播距離為二次三項式，若井底流壓為75 MPa時，數據相關數據代入式（28），并依此取啟動壓力梯度為002 MPa/m，007 MPa/m，020 MPa/m和滲透率為025 mD，050 mD，10 mD，20 mD得到定壓生產10 d的壓力傳播規律，如圖8所示。

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

2） 單因素分析。同理，如212所述，本文主要研究井底流壓、滲透率等對壓力傳播快慢的影響。分別取pwf=12 MPa，14 MPa，16 MPa，18 MPa觀察壓力

由圖9知：圖形斜率逐漸變小，壓力傳播速度減小，直到壓力波及到邊界。主要是因為傳播距離越遠，耗散的能量越大，若得不到外界能量的及時補給，則地層流體無法克服地層阻力，從而使得油井產量下降，這也是諸多低滲透油藏采用注水開發，補給地層能量的原因。當傳播距離不大時，井底流壓的影響較小，隨著距離的增大，井底流壓的影響增大。

減小井底流壓，放大生產壓差，可以有效提高油井產量，但是并非井底流壓越小越好，而是油井以最合理的井底流壓生產，可實現油井高效、穩定開發。（14）式解出x（t），并代入（17）式，求解得合理井底壓力為

pwf=pe-（a2-a）2b[（bfa-1）+

（bfa-1）2+8b2a2μBfkφ0ct（1-a）]（22）

f=（1-a）2[a+ab-（b-1）a2] （23）

將文獻17數據代入式（22）得合理井底壓力為1031 MPa，對應的最大產油量326 m3/d。

3實例分析

我國部分地區由于長期受到地質活動的影響，導致這些地區出現了很多的復雜斷塊油藏，蘇北的臺興油田是一個典型的狹長型低滲透復雜斷塊油田[22-23]。該油田由于受到多級斷層相互作用，使得臺興油田發育Ⅲ-Ⅴ級斷層11條，斷層將整個臺興油田分為11個含油斷塊。每個含油斷塊都有獨立的溫度、壓力系統，每個小油藏形狀近似矩形，寬度較小，長度較長，符合本文模型的應用條件。下面以文獻23中的區塊為例，進行實例分析。

目前數值模型的理論是建立在經典達西滲流基礎上的，而Eclipse2010中的E300模塊是針對低滲透油藏，考慮了啟動壓力梯度，因此將本文的定產模型計算的某時刻的井底流壓（定產解）與數值模擬結果（數值解）對比，如表2所示。由表2可知，相對誤差基本控制在8%內，說明定流壓模型是可靠的。

表2定產模型實例計算對比表

時間

/d定產解

/（m3·d-1）數值解

/（m3·d-1）相對誤差

/% 1017492161717552 2017202158837668 3016951156317787 4016729154127873 5016532152157966 6016355150438022 7016194148848089 8016048147388163 9015913146098195

由于定流壓模型中，設定井底流壓為常數，因此將式（17）帶入式（14），求得不同時刻的產量（定壓解），并與數值模擬（數值解）對比，如表3所示。由表3可知，相對誤差控制在77%內，則定流壓模型是可靠的。

表3定流壓模型實例計算對比表

時間

/d定壓解

/（m3·d-1）數值解

/（m3·d-1）相對誤差

/% 10271225097485 20285527134974 30327930726313 40313933807678 50365938966477 60423845216678 70560159295856 80691772464756 4結論

1） 定產模型中，時間與壓力傳播距離為二次三項式關系，且在壓力傳播速度上，兩參數連續模型比達西模型慢，比擬啟動壓力梯度模型要快，與流體黏度呈負相關，與產量、滲透率呈正相關；

2） 定壓模型中，時間與壓力傳播距離亦為二次三項式關系，但壓力傳播速度與滲透率呈正相關，與井底流壓、啟動壓力梯度呈負相關。

3） 實例計算表面，建立的定產模型、定壓模型是可靠的。

參考文獻：

[1]胡文瑞.中國低滲透油氣的現狀與未來[J].中國工程科學，2009，11（8）：29-37.

[2]王光付，廖榮鳳，李江龍，等.中國石化低滲透油藏開發狀況及前景[J].油氣地質與采收率，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，FARUK C.Mddification of Darcys law for the threshold pressure gradient，[J].Petrol.Sci.Eng.1999，22（4）：237-240.

[6]宋付權，劉慈群，李凡華.低滲透介質含啟動壓力梯度一維瞬時壓力分析[J].應用數學和力學，1999，20（1）：25-33.

[7]李凡華，劉慈群.含啟動壓力梯度的不定常滲流的壓力動態分析[J].油氣井測試，1997，6（1）：1-4.

[8]鄧英爾，劉慈群.低滲油藏非線性滲流規律數學模型及其應用[J].石油學報，2001，22（4）：72-77.

[9]劉鵬程，王曉冬，李素珍，等.地層壓敏對低滲透氣井產能影響研究 [J].西南石油學院學報，2004，26（5）：37-41.

[10]劉鵬程，王曉冬，萬玉金.三區復合油藏有限導流垂直裂縫井井底壓力動態分析[J].油氣井測試，2004，13（1）：4-7.

[11]郝明強，胡永樂，李凡華.微裂縫性特低滲透油藏產能研究[J].石油天然氣學報，2009，31（2）：100-104.

[12]朱圣舉.低滲透油藏的壓力波傳播規律[J].新疆石油地質，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]史瑞娜，王曉冬.低滲透地層壓力擾動傳播規律[J].特種油氣藏，2011， 18（4）：80-82.

[16]黃爽英，陳祖華.引入啟動壓力梯度計算低滲透砂巖油藏注水見效時間[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]李松泉，程林松，李秀生，等.特低滲透油藏非線性滲流模型[J].石油勘探與開發，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.