致密儲(chǔ)層體積壓裂井?dāng)?shù)值試井方法研究

王牧 冉旺 曹鵬程(中石油長(zhǎng)慶油田分公司第十二采油廠， 甘肅 合水 745400)

致密儲(chǔ)層體積壓裂井?dāng)?shù)值試井方法研究

王牧 冉旺 曹鵬程(中石油長(zhǎng)慶油田分公司第十二采油廠， 甘肅 合水 745400)

長(zhǎng)慶油田大規(guī)模開(kāi)展了以封隔器滑套分層壓裂技術(shù)、雙封單壓分段壓裂技術(shù)為主體的體積壓裂工藝措施，體積壓裂工藝技術(shù)的進(jìn)步帶動(dòng)了試井方法理論的發(fā)展。文章以“單孔雙區(qū)”和“雙孔雙區(qū)”物理模型為基礎(chǔ)，基于PEBI網(wǎng)格技術(shù)分別建立了體積壓裂井的兩種數(shù)值試井模型，求解后獲得了兩種模型的雙對(duì)數(shù)壓力響應(yīng)特征曲線，并對(duì)兩種模型的數(shù)值解和解析解進(jìn)行了比對(duì)分析，高度吻合的結(jié)果驗(yàn)證了數(shù)值模型的準(zhǔn)確性，同時(shí)對(duì)現(xiàn)場(chǎng)測(cè)試資料未能顯示第二徑向流段做出了理論解釋，文章最后通過(guò)對(duì)體積壓裂井周圍三維壓力場(chǎng)分布特征的研究，確定了壓力波的傳播規(guī)律，得出了影響油井產(chǎn)量的一個(gè)主要因素是改造區(qū)內(nèi)壓力的擾動(dòng)幅度。礦場(chǎng)的擬合效果肯定了該方法的實(shí)用性。

體積壓裂；滲流理論；數(shù)值試井；三維壓力場(chǎng)；特征曲線

Stimulated reservoir volume(SRV) treatment,mainly the layered fracturing technology by packer sleeve and double seal single pressure staged fracturing technology,were widely used in Changqing Oil Field.The improvement of SRV leads to the development of well testing theory.Based on the physical model of radial composite and double-porosity double-zone,two kinds of numerical well testing model of SRV wells was established by PEBI grid technology in this paper.The log-log pressure response characteristic curves of two models could be observed through model solving.CoMParison of numerical and analytical solution between two kinds of models shows considerably similar and reliable results.Further,a theoretical explain was given for field testing data’s lack of the second radial flow.In the last section of paper,we analyzed the distributing characteristics of 3D pressure field surrounding the SRV wells.Expanding pattern of pressure was established.It is revealed that pressure disturbance in stimulating area is a main factor to well production.Perfect match on field confirmed practicability of this new method.

長(zhǎng)慶油田借鑒美國(guó)儲(chǔ)層改造經(jīng)驗(yàn)，首次將儲(chǔ)層改造技術(shù)引入國(guó)內(nèi)低滲透油藏，并賦予了新的名稱“體積壓裂”，作為一個(gè)新概念的引入，體積壓裂井滲流機(jī)理研究目前尚少，2012年姜瑞忠、蔡振華等人分別提出了低滲透油藏壓裂井生態(tài)動(dòng)態(tài)分析方法和致密氣藏壓裂井的產(chǎn)能評(píng)價(jià)方法[1-2]；2013年王建忠、姚軍等人建立了裂縫性低滲透油藏的試井模型[3]；隨后姚軍等人在Ozkan基本點(diǎn)源解的基礎(chǔ)上，推導(dǎo)出不同邊界條件下裂縫性油氣藏分段壓裂水平井單條及多條裂縫生產(chǎn)時(shí)的井底無(wú)因次壓降求解曲線[4]；同年，蔡田田、劉振宇等人采用有限元數(shù)值模擬方法，對(duì)體積壓裂復(fù)雜縫網(wǎng)進(jìn)行數(shù)值模擬研究，建立了壓裂井生產(chǎn)動(dòng)態(tài)預(yù)測(cè)的數(shù)學(xué)模型[5]；而劉雄王曉冬等人考慮啟動(dòng)壓力梯度的影響，利用分形體描述改造區(qū)的特征，提出了預(yù)測(cè)體積壓裂直井和水平井的穩(wěn)態(tài)/非穩(wěn)態(tài)產(chǎn)能的解析方法[6]；2014年陳曉明等人針對(duì)長(zhǎng)慶油田低滲透油藏老井體積壓裂情況，提出了描述體積壓裂直井的兩種不穩(wěn)定試井模型：?jiǎn)慰纂p區(qū)模型[7]和雙孔雙區(qū)模型[8]。這兩種模型在長(zhǎng)慶油田的老井體積壓裂測(cè)試分析過(guò)程中得到了很好的應(yīng)用，為進(jìn)一步拓展單孔雙區(qū)模型和雙孔雙區(qū)模型的礦場(chǎng)適用性，本文提出一種體積壓裂井的數(shù)值試井分析方法。

1 物理模型基礎(chǔ)

圖 1 雙孔雙區(qū)模型 2D Map

圖 2 雙孔雙區(qū)模型 2D Map

吳奇、陳作等人針對(duì)美國(guó)頁(yè)巖氣體積改造技術(shù)的內(nèi)涵以及在我國(guó)的應(yīng)用等方面做了深入的探討，從儲(chǔ)層巖性、裂縫破壞方式、天然裂縫溝通情況以及壓裂方式等四個(gè)方面提出了體積壓裂技術(shù)的關(guān)鍵[9-11]。以此為指導(dǎo)，長(zhǎng)慶油田開(kāi)展了以“封隔器滑套分層壓裂”為主體技術(shù)的直井體積壓裂措施和以“雙封單壓分段壓裂”為主體技術(shù)的水平井體積壓裂措施?，F(xiàn)場(chǎng)測(cè)試資料表明，體積壓裂后的測(cè)試曲線主要存在兩種特征響應(yīng)段：“雙線性流段”和“基質(zhì)竄流段”[12]，以這兩種特征響應(yīng)為依據(jù)，陳曉明等人提出了單孔雙區(qū)和雙孔雙區(qū)的概念。

針對(duì)陳曉明等人提出的與兩種儲(chǔ)層破壞結(jié)構(gòu)形式相對(duì)應(yīng)的滲流物理模型，借助于PEBI網(wǎng)格劃分方法分別建立了與其相對(duì)應(yīng)的單井?dāng)?shù)值模型。圖1為單孔雙區(qū)模型的平面圖，其中測(cè)試井為體積壓裂措施井，呈現(xiàn)“兩翼縫”特征，近井地帶是儲(chǔ)層改造區(qū)，儲(chǔ)層是單一介質(zhì)；遠(yuǎn)井地帶(錨區(qū))是儲(chǔ)層非改造區(qū)，儲(chǔ)層也是單一介質(zhì)。圖2為雙孔雙區(qū)模型的平面圖，其中測(cè)試井為體積壓裂措施井，呈現(xiàn)“縫網(wǎng)”特征，近井地帶是儲(chǔ)層改造區(qū)，儲(chǔ)層是雙重介質(zhì)；遠(yuǎn)井地帶(錨區(qū))是儲(chǔ)層非改造區(qū)，儲(chǔ)層是單一介質(zhì)。為確保數(shù)值模擬的結(jié)果更貼近實(shí)際，模型設(shè)計(jì)的參數(shù)在此取長(zhǎng)慶油田低滲透儲(chǔ)層的平均值，各個(gè)參數(shù)的具體取值如表1和表2。

表1 單孔雙區(qū)模型試井設(shè)計(jì)參數(shù)一覽表

表2 雙孔雙區(qū)模型試井設(shè)計(jì)參數(shù)一覽表

2 數(shù)值模型結(jié)果分析

2.1 單孔雙區(qū)模型壓力響應(yīng)特征曲線分析

圖3是經(jīng)過(guò)數(shù)值計(jì)算得出的井底壓力及壓力導(dǎo)數(shù)雙對(duì)數(shù)曲線，可以將曲線形態(tài)劃分為以下四個(gè)流動(dòng)階段：

(1)第Ⅰ段是早期部分，壓力與壓力導(dǎo)數(shù)曲線合并成為45度直線，表明續(xù)流的影響。

(2)第Ⅱ段斜率近似為1/4，反映了裂縫-地層雙線性流，其持續(xù)時(shí)間長(zhǎng)短與裂縫半長(zhǎng)和導(dǎo)流能力有關(guān)。

(3) 第Ⅲ段呈現(xiàn)平臺(tái)狀，達(dá)到了第一徑向流階段，反映的是改造區(qū)的流動(dòng)系數(shù)。

(4)第Ⅳ段曲線開(kāi)始上翹，預(yù)示著壓力波開(kāi)始向外區(qū)過(guò)渡。

該模型的數(shù)值解中沒(méi)有表現(xiàn)出第二徑向流階段的特征，說(shuō)明了實(shí)際測(cè)試時(shí)間段內(nèi)，儲(chǔ)層中壓力降未能傳播到外區(qū)。(圖3)

將該模型的設(shè)計(jì)參數(shù)及數(shù)值解無(wú)因次化處理，然后將無(wú)因次化后的模型參數(shù)輸入單孔雙區(qū)解析模型的求解程序中，可得此模型的半解析解。對(duì)比數(shù)值解與半解析解如圖4所示，兩種算法下所得到的壓力及壓力導(dǎo)數(shù)曲線基本重合(誤差小于1%)，說(shuō)明了單孔雙區(qū)數(shù)值模型可靠。模型參數(shù)無(wú)因次化方法如下：

式中：Ct為壓縮系數(shù)，MPa-1； C為井筒存儲(chǔ)系數(shù)，m3/MPa；h為地層厚度，m； wf為裂縫寬度，m， xf為裂縫半長(zhǎng)，m； w1為內(nèi)區(qū)滲透率，mD；PwD為垂直裂縫井的井底壓力；q為井的產(chǎn)量，m3/D；t為時(shí)間，h；μ為流體的粘度，mPa*s。

2.2 雙孔雙區(qū)模型壓力響應(yīng)特征曲線分析

圖5是由數(shù)值模型計(jì)算得出的井底壓力及壓力導(dǎo)數(shù)雙對(duì)數(shù)曲線，曲線形態(tài)可以劃分為以下四個(gè)階段：

(1)第Ⅰ段是早期部分，壓力與壓力導(dǎo)數(shù)曲線合并成為45度直線，表明續(xù)流的影響，井筒儲(chǔ)集與表皮效應(yīng)結(jié)束后導(dǎo)數(shù)值向下傾斜。

(2)第Ⅱ段出現(xiàn)下凹段，反映了基質(zhì)-裂縫的竄流效應(yīng)，其出現(xiàn)的早晚以及持續(xù)的長(zhǎng)短與竄流系數(shù)和儲(chǔ)容比有關(guān)。

(3) 第Ⅲ段呈現(xiàn)平臺(tái)狀，達(dá)到了第一徑向流階段，反映的是改造區(qū)的流動(dòng)系數(shù)。

(4)第Ⅳ段曲線開(kāi)始上翹，預(yù)示著壓力波開(kāi)始向外區(qū)過(guò)渡。

該模型的數(shù)值解中同樣沒(méi)有表現(xiàn)出第二徑向流階段的特征，說(shuō)明了實(shí)際測(cè)試時(shí)間段內(nèi)，儲(chǔ)層中壓力降未能傳播到外區(qū)。

同理，將雙孔雙區(qū)數(shù)值模型的設(shè)計(jì)參數(shù)及數(shù)值解經(jīng)無(wú)因次化處理，然后將無(wú)因次化后的模型參數(shù)輸入雙孔雙區(qū)解析模型的求解程序中可得此模型的解析解。對(duì)比數(shù)值解與半解析解如圖6所示，兩種算法下所得到的壓力及壓力導(dǎo)數(shù)曲線基本重合(誤差小于1%)，說(shuō)明雙孔雙區(qū)數(shù)值模型也是可靠的。模型參數(shù)無(wú)因次化方法如下：

圖3 單孔雙區(qū)模型壓力及壓力導(dǎo)數(shù)曲線

圖4 單孔雙區(qū)模型壓力及壓力導(dǎo)數(shù)曲線驗(yàn)證

圖5 雙孔雙區(qū)模型壓力及壓力導(dǎo)數(shù)曲線

圖6 雙孔雙區(qū)模型壓力及壓力導(dǎo)數(shù)曲線驗(yàn)證

2.3 體積壓裂井壓力場(chǎng)特征分析

由圖4至圖6可知，數(shù)值解求得的壓力及壓力導(dǎo)數(shù)雙對(duì)數(shù)曲線并不能完整顯示出第二徑向流段，來(lái)自現(xiàn)場(chǎng)的測(cè)試曲線同樣很難出現(xiàn)第二徑向流特征段。本文通過(guò)分析壓力場(chǎng)的特征對(duì)這種現(xiàn)象做出解釋。

圖7是單孔雙區(qū)模型下單井生產(chǎn)了100小時(shí)、300小時(shí)、600小時(shí)后的三維空間壓力分布示意圖，從該圖中可以看出，改造區(qū)域內(nèi)的壓力衰竭幅度較大，而非改造區(qū)域內(nèi)的壓力衰竭幅度很小，這一方面說(shuō)明了儲(chǔ)層物性的改善使改造區(qū)域內(nèi)產(chǎn)生了較大的壓力擾動(dòng)，另一方面說(shuō)明了在低滲油藏原始儲(chǔ)層物性狀態(tài)下壓力傳播的速度是非常有限的，這也是導(dǎo)致油井低產(chǎn)的一個(gè)重要原因。(圖7)

圖8是雙孔雙區(qū)模型下單井生產(chǎn)160小時(shí)、320小時(shí)、480小時(shí)后的三維空間壓力分布示意圖，同樣在改造區(qū)域內(nèi)的壓力衰竭程度遠(yuǎn)遠(yuǎn)大于非改造區(qū)，說(shuō)明了體積壓裂后產(chǎn)量供給主要是來(lái)自改造區(qū)。(圖8)

3 礦場(chǎng)應(yīng)用

取長(zhǎng)慶油田XX井的測(cè)試數(shù)據(jù)進(jìn)行該方法的應(yīng)用。儲(chǔ)層的基礎(chǔ)參數(shù)為：儲(chǔ)層平均厚度15 m，平均孔隙度10%，井筒半徑0．6 m，原油粘度3．14 MPa·s，原油體積系數(shù)1．2，綜合壓縮系數(shù)15．38E-4 MPa-1，測(cè)試前井的產(chǎn)量10．8 m3/D。

圖7 單孔雙區(qū)模型壓力空間分布圖

圖8 雙孔雙區(qū)模型壓力空間分布圖

測(cè)試曲線顯示出“竄流”特征段，因此建立雙孔雙區(qū)數(shù)值模型，當(dāng)參數(shù)設(shè)置為以下值時(shí)計(jì)算測(cè)試井的井底響應(yīng)壓力值與實(shí)測(cè)值吻合較好：內(nèi)區(qū)滲透率K=2．5md，表皮系數(shù)S=-2，彈性儲(chǔ)能比Omega=0．3，竄流系數(shù)Lambda=5．0E-4，流度比M=8，內(nèi)區(qū)半徑R=60 m。擬合結(jié)果如圖9所示。

圖 9 測(cè)試井?dāng)?shù)值模擬結(jié)果與實(shí)測(cè)結(jié)果擬合圖

4 結(jié)論

(1)以單孔雙區(qū)滲流物理模型和雙孔雙區(qū)滲流物理模型為基礎(chǔ)分別建立了與其相對(duì)應(yīng)的數(shù)值試井模型，并求得了數(shù)值解。

(2)兩種模型數(shù)值解與解析解的高度吻合分別驗(yàn)證了單孔雙區(qū)不穩(wěn)定試井模型和雙孔雙區(qū)不穩(wěn)定試井模型的準(zhǔn)確性，為體積壓裂試井技術(shù)的發(fā)展奠定了理論基礎(chǔ)。

(3)數(shù)值模型的結(jié)果表明體積壓裂后，壓力波的擾動(dòng)主要發(fā)生在近井地帶，這對(duì)礦場(chǎng)測(cè)試不到第二徑向流的現(xiàn)象給出了合理解釋。

(4)數(shù)值模型可以彌補(bǔ)解析模型中無(wú)法考慮非均質(zhì)性、鄰井干擾等因素，拓展了“單孔雙區(qū)”與“雙孔雙區(qū)”模型在礦場(chǎng)的應(yīng)用性。

[1]姜瑞忠,楊明,王公昌,等.低滲透油藏壓裂井生產(chǎn)動(dòng)態(tài)分析[J].特種油氣藏,2013,20(1):52-55.

[2]蔡振華,廖新維.非常規(guī)氣藏壓裂井產(chǎn)能評(píng)價(jià)方法[J].特種油氣藏,2013,20(4):96-98.

[3]王建忠,姚軍.裂縫性低滲透油藏試井模型與非穩(wěn)態(tài)壓力特征[J].特種油氣藏,2013,20(2):69-71.

[4]姚軍,劉丕養(yǎng),吳明錄.裂縫性油氣藏壓裂水平井試井分析[J].中國(guó)石油大學(xué)學(xué)報(bào):自然科學(xué)版,2013,37(5):107-113.

[5]蔡田田.低滲透油藏體積壓裂數(shù)值模擬研究[D]:東北石油大學(xué),2013.

[6]劉雄.致密油儲(chǔ)層體積改造產(chǎn)能評(píng)價(jià)方法研究[D]:中國(guó)地質(zhì)大學(xué)(北京),2013.

[7]陳曉明,廖新維,趙曉亮等.直井體積壓裂不穩(wěn)定試井研究—單孔雙區(qū)模型[J].科學(xué)技術(shù)與工程,2014,14(26):45-49.

[8]陳曉明,廖新維,李東暉等.直井體積壓裂不穩(wěn)定試井研究—雙孔雙區(qū)模型[J].油氣井測(cè)試,2014,23(4):4-8.

[9]吳奇,胥云,劉玉章.美國(guó)頁(yè)巖氣體積改造技術(shù)現(xiàn)狀及對(duì)我國(guó)的啟示[J].石油鉆采工藝,2011,33(2):1-7.

[10]吳奇,胥云,王曉泉.非常規(guī)油氣藏體積改造技術(shù)——內(nèi)涵,優(yōu)化設(shè)計(jì)與實(shí)現(xiàn)[J].石油勘探與開(kāi)發(fā),2012,39(3):352-358.

[11]陳作,薛承瑾,蔣廷學(xué).頁(yè)巖氣井體積壓裂技術(shù)在我國(guó)的應(yīng)用建議[J].天然氣工業(yè),2010,30(10):30-32.

[12]趙曉亮,廖新維,王歡等.體積壓裂復(fù)雜裂縫形態(tài)壓力響應(yīng)特征分析[J].陜西科技大學(xué)學(xué)報(bào),2014,32(3):89-92.

The Numerical Well Testing Method of Volume Fracturing Wells in Tight Rerservoir

Wang Mu(Changqing Oilfield CoMPany,CNPC,Xi’an 710000,Shaanxi Province.)

王牧(1987- )，男， 中國(guó)石油長(zhǎng)慶油田分公司第十二采油廠，油氣田開(kāi)發(fā)專業(yè)。