劉禮軍， 姚 軍， 孫 海， 白玉湖， 徐兵祥， 陳 嶺

(1.中國(guó)石油大學(xué)(華東)油氣滲流研究中心，青島 266580；2.中海油研究總院，北京 100028)

?油氣開發(fā)?

考慮啟動(dòng)壓力梯度和應(yīng)力敏感的頁(yè)巖油井產(chǎn)能分析

劉禮軍1， 姚 軍1， 孫 海1， 白玉湖2， 徐兵祥2， 陳 嶺2

(1.中國(guó)石油大學(xué)(華東)油氣滲流研究中心，青島 266580；2.中海油研究總院，北京 100028)

為研究啟動(dòng)壓力梯度及應(yīng)力敏感效應(yīng)對(duì)頁(yè)巖油井產(chǎn)能的影響規(guī)律，基于油氣水三相滲流模型，建立了考慮啟動(dòng)壓力梯度和應(yīng)力敏感效應(yīng)的頁(yè)巖油滲流數(shù)學(xué)模型，給出了數(shù)值求解方法，并編制了頁(yè)巖油藏?cái)?shù)值模擬器。利用頁(yè)巖油藏?cái)?shù)值模擬器分析了啟動(dòng)壓力梯度和應(yīng)力敏感效應(yīng)對(duì)多級(jí)壓裂頁(yè)巖油水平井產(chǎn)能的影響規(guī)律。頁(yè)巖油藏?cái)?shù)值模擬器退化計(jì)算結(jié)果與Eclipse軟件的計(jì)算結(jié)果基本吻合，驗(yàn)證了模擬器的準(zhǔn)確性。數(shù)值模擬結(jié)果表明：當(dāng)啟動(dòng)壓力梯度達(dá)到0.5 MPa/m時(shí)，生產(chǎn)10 000 d的累計(jì)產(chǎn)油量為無(wú)啟動(dòng)壓力梯度時(shí)的24.7%；當(dāng)應(yīng)力敏感系數(shù)由0.1 MPa-1升至0.5 MPa-1時(shí)，生產(chǎn)10 000 d時(shí)儲(chǔ)層近井地帶的滲透率降低近一個(gè)數(shù)量級(jí)，累計(jì)產(chǎn)油量為無(wú)應(yīng)力敏感效應(yīng)時(shí)的36.2%；同時(shí)考慮啟動(dòng)壓力梯度和應(yīng)力敏感效應(yīng)時(shí)，累計(jì)產(chǎn)油量更低。研究結(jié)果表明，啟動(dòng)壓力梯度和應(yīng)力敏感效應(yīng)對(duì)頁(yè)巖油井的產(chǎn)能均有較大的抑制作用，如不考慮啟動(dòng)壓力梯度和應(yīng)力敏感效應(yīng)會(huì)過(guò)高估計(jì)頁(yè)巖油井的產(chǎn)能，頁(yè)巖油藏?cái)?shù)值模擬器可為頁(yè)巖油藏的開發(fā)提供指導(dǎo)。

頁(yè)巖油；啟動(dòng)壓力梯度；應(yīng)力敏感；產(chǎn)能

頁(yè)巖油儲(chǔ)層發(fā)育大量的納米級(jí)孔隙，微裂縫也非常發(fā)育，屬于超低滲透的致密儲(chǔ)層，需要通過(guò)大規(guī)模壓裂形成納米級(jí)孔隙-微裂縫-大裂縫共存的多尺度頁(yè)巖儲(chǔ)層[1]，才能有效開發(fā)。頁(yè)巖油儲(chǔ)層具有超低滲透特性，經(jīng)大規(guī)模壓裂后又存在多尺度孔隙空間，其滲流機(jī)理和流動(dòng)特征均與常規(guī)油氣藏不同[2]，因此，明確頁(yè)巖油藏的復(fù)雜滲流機(jī)理對(duì)頁(yè)巖油開發(fā)具有重要意義。國(guó)內(nèi)外學(xué)者對(duì)頁(yè)巖油藏的滲流機(jī)理進(jìn)行了多方面的研究：高英等人[3]認(rèn)為頁(yè)巖油滲流為低速非達(dá)西滲流，具有啟動(dòng)壓力梯度，滲透率越小啟動(dòng)壓力梯度越大，并基于分形理論建立了頁(yè)巖油藏壓裂井的產(chǎn)能解析模型；郭肖等人[4]在建立頁(yè)巖滲透率模型時(shí)考慮了應(yīng)力敏感效應(yīng)的影響，分析了儲(chǔ)層壓力對(duì)滲透率的影響，發(fā)現(xiàn)滲透率隨儲(chǔ)層壓力降低而降低；S.S.Chhatre等人[5]通過(guò)穩(wěn)態(tài)試驗(yàn)發(fā)現(xiàn)頁(yè)巖的滲透率隨凈圍壓升高而降低，且兩者呈指數(shù)關(guān)系。針對(duì)頁(yè)巖油儲(chǔ)層中存在的微裂縫和人工大裂縫，傳統(tǒng)的數(shù)值模擬方法通常采用雙重介質(zhì)模型[6]和網(wǎng)格加密等方式處理，但雙重介質(zhì)模型在處理頁(yè)巖油藏中的非穩(wěn)態(tài)竄流時(shí)通常誤差較大[7]，網(wǎng)格加密會(huì)增加計(jì)算量，且難以處理具有復(fù)雜分布特征的大裂縫。MINC(multiple interacting continua)模型[8]將基質(zhì)網(wǎng)格劃分為嵌套網(wǎng)格，解決了基質(zhì)與微裂縫間的非穩(wěn)態(tài)竄流問題。嵌入式離散裂縫模型[9-10]將裂縫單元直接嵌入結(jié)構(gòu)化網(wǎng)格中，有效減少了網(wǎng)格量，且易于模擬復(fù)雜的裂縫。

截至目前，許多學(xué)者分析研究了啟動(dòng)壓力梯度和應(yīng)力敏感對(duì)低滲透油氣藏的影響[11-17]，但都存在不足或不完善之處，如：楊仁鋒等人[16]考慮啟動(dòng)壓力梯度和壓敏效應(yīng)，對(duì)特低滲透油藏的非線性滲流進(jìn)行了數(shù)值模擬，取得了一定的認(rèn)識(shí)，但其以直井為研究對(duì)象，未考慮多級(jí)壓裂及微裂縫發(fā)育的情況；熊健等人[17]考慮啟動(dòng)壓力梯度、滑脫及應(yīng)力敏感等效應(yīng)，推導(dǎo)了非線性滲流下低滲氣藏壓裂井的產(chǎn)能方程，但其建立的物理模型較為簡(jiǎn)單，也未考慮大量微裂縫發(fā)育的情況。為此，筆者建立了考慮啟動(dòng)壓力梯度和應(yīng)力敏感效應(yīng)的頁(yè)巖油藏滲流數(shù)學(xué)模型，采用有限體積法對(duì)其進(jìn)行離散求解，并采用MINC模型和嵌入式離散裂縫模型處理頁(yè)巖儲(chǔ)層中的微裂縫和大裂縫，有效解決了多尺度孔隙特征下頁(yè)巖油藏中的油氣滲流問題，并在此基礎(chǔ)上分析了啟動(dòng)壓力梯度和應(yīng)力敏感效應(yīng)對(duì)頁(yè)巖油井產(chǎn)能的影響程度。

1 頁(yè)巖油藏滲流數(shù)學(xué)模型的建立

頁(yè)巖油密度低，一般含有溶解氣，且地層水可能參與流動(dòng)，因此建立頁(yè)巖油藏油、氣、水三相滲流數(shù)學(xué)模型時(shí)，根據(jù)文獻(xiàn)[18]作以下基本假設(shè)：1)頁(yè)巖油藏中的滲流為等溫滲流；2)油相和水相的滲流考慮啟動(dòng)壓力梯度和應(yīng)力敏感，氣相滲流考慮應(yīng)力敏感；3)頁(yè)巖油藏?zé)N類只含油和氣2種組分，油組分只存在于油相中，氣組分可存在于氣相和油相中；4)油、氣不溶于水。

1.1啟動(dòng)壓力梯度

頁(yè)巖油儲(chǔ)層的滲透率極低，一般為納達(dá)西量級(jí)，具有明顯的非線性滲流特征[3]?？紤]啟動(dòng)壓力梯度的滲流運(yùn)動(dòng)方程為：

式中：v為滲流速度，m/s；Kr為相對(duì)滲流率；K為滲透率，D；μ為黏度，Pa·s；G為啟動(dòng)壓力梯度，Pa/m；ψ為勢(shì)梯度，Pa/m。

1.2應(yīng)力敏感效應(yīng)

在頁(yè)巖油開采過(guò)程中，儲(chǔ)層壓力下降使巖石所受有效應(yīng)力增加，從而導(dǎo)致儲(chǔ)層滲透率降低。頁(yè)巖滲透率與有效應(yīng)力之間的關(guān)系式為[19-20]：

K=K0e-c(p-p0)

(2)

式中：K0為原始狀態(tài)下的滲透率，D；p0為原始油藏壓力，Pa；c為應(yīng)力敏感系數(shù)，Pa-1。

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

頁(yè)巖油藏中油、氣、水三相的滲流數(shù)學(xué)模型為：

-·(ρovo)+qo=

(3)

(4)

(5)

式中：vo，vg和vw分別為油、氣、水相的滲流速度，m/s；ρo和ρgo分別為油相中油組分和溶解氣組分的密度，kg/m3；ρg和ρw分別為氣相和水相的密度，kg/m3；qo，qg和qw分別為油、氣、水組分的源匯項(xiàng)，即單位時(shí)間內(nèi)單位地層體積的產(chǎn)出或注入量，kg/(m3·s)；φ為孔隙度；So，Sg和Sw為油、氣、水相的飽和度。

此外，為了使方程組封閉，還需要如下輔助方程：

So+Sg+Sw=1

(6)

pw=po-pcow

(7)

pg=po-pcog

(8)

式中：po，pg和pw分別為油、氣、水相的壓力，Pa；pcow和pcog為油水和油氣的毛管力，Pa。

由式(1)—式(8)組成的方程組即為頁(yè)巖油藏滲流數(shù)學(xué)模型。

2 頁(yè)巖油藏滲流模型數(shù)值求解

2.1數(shù)值離散格式

采用有限體積法對(duì)頁(yè)巖油藏滲流數(shù)學(xué)模型進(jìn)行數(shù)值離散，即對(duì)滲流基本方程式(3)—(5)進(jìn)行積分，然后進(jìn)行積分離散，并將運(yùn)動(dòng)方程式(1)代入可得到滲流數(shù)學(xué)模型的數(shù)值離散格式：

(9)

(10)

(11)

ψ=p-ρgD

(12)

根據(jù)離散格式的方程式(9)—(11)，得到方程的殘差格式為：

(13)

(14)

(15)

式中：Ro,n，Rg,n和Rw,n分別為油、氣、水滲流方程的殘差。

式(13)—式(15)構(gòu)成了頁(yè)巖油藏滲流數(shù)學(xué)模型的全隱式數(shù)值格式，可通過(guò)牛頓迭代法對(duì)方程組進(jìn)行求解，將求解過(guò)程編制為計(jì)算程序就形成了頁(yè)巖油藏?cái)?shù)值模擬器。

2.2微裂縫及大裂縫的處理

采用MINC模型和嵌入式離散裂縫模型處理儲(chǔ)層中存在的微裂縫和大裂縫。上述2種模型的核心為網(wǎng)格連接信息的處理。MINC模型是在雙重介質(zhì)的基礎(chǔ)上，將基質(zhì)網(wǎng)格進(jìn)一步細(xì)分為嵌套網(wǎng)格，并基于網(wǎng)格劃分情況計(jì)算網(wǎng)格間的連接信息，如網(wǎng)格間中心距離、網(wǎng)格間界面面積，具體計(jì)算可參考文獻(xiàn)[8]。嵌入式離散裂縫模型是將大裂縫直接嵌入到規(guī)則的基質(zhì)網(wǎng)格中，計(jì)算大裂縫網(wǎng)格與基質(zhì)網(wǎng)格間的連接信息[10]。MINC模型及嵌入式離散裂縫模型劃分所得的網(wǎng)格見圖1，對(duì)應(yīng)的網(wǎng)格連接關(guān)系見圖2。采用MINC模型和嵌入式離散裂縫模型對(duì)微裂縫和大裂縫的處理形成網(wǎng)格前處理模塊，其與頁(yè)巖油藏?cái)?shù)值模擬器結(jié)合，便可對(duì)采用多級(jí)壓裂水平井開發(fā)的頁(yè)巖油藏進(jìn)行模擬分析。

圖1 網(wǎng)格示意Fig.1 Schematic diagram of grid block

圖2 網(wǎng)格連接關(guān)系示意Fig.2 Schematic diagram of grid connections

3 模型驗(yàn)證

商業(yè)模擬器在功能上存在局限，且啟動(dòng)壓力梯度及應(yīng)力敏感的數(shù)學(xué)表達(dá)式有多種形式，因此對(duì)筆者所編制的模擬器進(jìn)行了退化驗(yàn)證(c=0Pa-1，G=0Pa/m)。采用筆者所編模擬器和Eclipse軟件對(duì)一個(gè)二維頁(yè)巖油藏模型(見圖3)進(jìn)行了模擬，并對(duì)比了模擬結(jié)果。二維頁(yè)巖油藏模型的外邊界封閉，水平井以井底流壓25MPa進(jìn)行定壓生產(chǎn)，原始油藏壓力為31.27MPa，網(wǎng)格數(shù)為100×51×1，3個(gè)方向上網(wǎng)格尺寸均為10m，滲透率為0.1μD，孔隙度為0.08，地層水的黏度為1.000mPa·s，原油黏度為0.132mPa·s，束縛水飽和度為0.22，水力裂縫的滲透率為4D，開度為0.001m。

圖3 二維頁(yè)巖油藏的物理模型Fig.3 Physical model of 2D shale oil reservoir

圖4所示為分別用Eclipse軟件和筆者所編模擬器模擬該油藏生產(chǎn)10000d時(shí)的壓力分布。從圖4可以看出，筆者所編模擬器模擬得到的壓力分布與Eclipse軟件模擬得到的壓力分布基本吻合。

圖4 壓力分布模擬結(jié)果Fig.4 Simulation results of pressure distribution

圖5所示為用Eclipse軟件和筆者所編模擬器計(jì)算得到的生產(chǎn)10000d的日產(chǎn)油量和累計(jì)產(chǎn)油量。從圖5可以看出，筆者所編模擬器計(jì)算得到的日產(chǎn)油量和累計(jì)產(chǎn)油量與Eclipse軟件計(jì)算的結(jié)果基本吻合。

通過(guò)對(duì)比筆者所編模擬器與Eclipse軟件的模擬結(jié)果，驗(yàn)證了筆者所建立數(shù)學(xué)模型的準(zhǔn)確性。

4 產(chǎn)能影響因素分析

4.1啟動(dòng)壓力梯度對(duì)產(chǎn)能的影響

為了研究在定壓衰竭開采條件下啟動(dòng)壓力梯度對(duì)頁(yè)巖油井產(chǎn)能的影響，利用模擬器對(duì)不同啟動(dòng)壓力梯度下頁(yè)巖油井的產(chǎn)能進(jìn)行了模擬，模擬條件與模型驗(yàn)證時(shí)相同。此外，在上述油藏物理模型的基礎(chǔ)上加入微裂縫，微裂縫孔隙度為0.002%，滲透率為0.1mD，微裂縫間距為10m，井底壓力為21MPa。啟動(dòng)壓力梯度分別取0.1和0.5MPa/m[21]。

圖5 產(chǎn)油量和累計(jì)產(chǎn)油量的模擬結(jié)果Fig.5 Simulation results of daily oil production and cumulative oil production

圖6所示為不同啟動(dòng)壓力梯度下生產(chǎn)10000d時(shí)的壓力分布。由圖6可知，啟動(dòng)壓力梯度的存在抑制了壓力波的傳播，隨著頁(yè)巖油藏啟動(dòng)壓力梯度的增大，壓力波的擴(kuò)散范圍縮小。

圖6 不同啟動(dòng)壓力梯度下的壓力分布Fig.6 Pressure distribution under different threshold pressure gradient

圖7為不同啟動(dòng)壓力梯度下生產(chǎn)10000d的日產(chǎn)油量和累計(jì)產(chǎn)油量。從圖7可以看出，隨啟動(dòng)壓力梯度增大，頁(yè)巖油井的產(chǎn)油量降低，當(dāng)啟動(dòng)壓力梯度升至0.5MPa/m時(shí)，生產(chǎn)10000d時(shí)的累計(jì)產(chǎn)油量降為無(wú)啟動(dòng)壓力梯度的24.7%。啟動(dòng)壓力梯度的影響與生產(chǎn)壓差有關(guān)，隨著生產(chǎn)壓差增大，啟動(dòng)壓力梯度的影響逐漸減弱，但總的來(lái)說(shuō)，啟動(dòng)壓力梯度對(duì)頁(yè)巖油井產(chǎn)能具有較大的抑制作用。

圖7 不同啟動(dòng)壓力梯度下的頁(yè)巖油井的產(chǎn)油量Fig.7 Shale oil production under different threshold pressure gradient

4.2應(yīng)力敏感效應(yīng)對(duì)產(chǎn)能的影響

為了研究在定壓衰竭開采條件下應(yīng)力敏感效應(yīng)對(duì)頁(yè)巖油井產(chǎn)能的影響，對(duì)不同應(yīng)力敏感系數(shù)下頁(yè)巖油井的產(chǎn)能進(jìn)行了模擬，模擬條件與模型驗(yàn)證時(shí)相同。假設(shè)基質(zhì)、微裂縫和大裂縫的應(yīng)力敏感系數(shù)相同，并分別取0.1，0.3和0.5MPa-1[22]。

由于頁(yè)巖油藏基質(zhì)的滲透率極低，對(duì)頁(yè)巖油流動(dòng)貢獻(xiàn)很小，而微裂縫為主要的頁(yè)巖油滲流通道，因此分析微裂縫滲透率在應(yīng)力敏感效應(yīng)下的響應(yīng)，對(duì)于認(rèn)識(shí)頁(yè)巖油井的產(chǎn)能規(guī)律具有重要的指導(dǎo)作用。圖8所示為不同應(yīng)力敏感系數(shù)下生產(chǎn)10000d時(shí)微裂縫的滲透率分布。從圖8可以看出，在開發(fā)頁(yè)巖油的過(guò)程中，應(yīng)力敏感效應(yīng)會(huì)使微裂縫的滲透率降低，尤其是近井地帶，而且隨著應(yīng)力敏感系數(shù)增大，微裂縫滲透率下降幅度增大，當(dāng)應(yīng)力敏感系數(shù)由0.1MPa-1升至0.5MPa-1時(shí)，生產(chǎn)10000d時(shí)近井處微裂縫的滲透率降低近一個(gè)數(shù)量級(jí)。

圖8 不同應(yīng)力敏感系數(shù)下的微裂縫滲透率分布Fig.8 Micro-fracture permeability distribution with different stress sensitivity coefficientst

圖9所示為不同應(yīng)力敏感系數(shù)下生產(chǎn)10000d時(shí)的日產(chǎn)油量和累計(jì)產(chǎn)油量。從圖9可以看出：應(yīng)力敏感效應(yīng)的影響主要體現(xiàn)在開發(fā)初期，這主要是由于開發(fā)初期儲(chǔ)層壓力下降快，大裂縫、微裂縫及基質(zhì)的滲透率下降快，從而導(dǎo)致開發(fā)初期油井的產(chǎn)油量大幅降低，當(dāng)進(jìn)入開發(fā)后期，儲(chǔ)層壓力變化平緩，產(chǎn)油量變平穩(wěn)；應(yīng)力敏感效應(yīng)對(duì)油井產(chǎn)能影響顯著，當(dāng)應(yīng)力敏感系數(shù)達(dá)到0.5MPa-1時(shí)，生產(chǎn)10000d的累計(jì)產(chǎn)油量為無(wú)應(yīng)力敏感效應(yīng)的36.2%。同樣，應(yīng)力敏感效應(yīng)的影響與生產(chǎn)壓差有關(guān)，隨著生產(chǎn)壓差增大，應(yīng)力敏感效應(yīng)的影響逐漸顯著?？傮w來(lái)說(shuō)，應(yīng)力敏感效應(yīng)對(duì)頁(yè)巖油井產(chǎn)能具有很大的抑制作用。

圖9 不同應(yīng)力敏感系數(shù)下的頁(yè)巖油井的產(chǎn)油量Fig.9 Shale oil production with different stress sensitivity coefficients

4.3應(yīng)力敏感系數(shù)及啟動(dòng)壓力梯度對(duì)產(chǎn)能的影響

對(duì)不同應(yīng)力敏感系數(shù)及啟動(dòng)壓力梯度參數(shù)組合下頁(yè)巖油井的產(chǎn)能進(jìn)行了模擬，模擬條件與模型驗(yàn)證時(shí)相同。應(yīng)力敏感系數(shù)和啟動(dòng)壓力分別取0.3MPa-1、0.3MPa/m和0.5MPa-1、0.5MPa/m。

圖10所示為不同應(yīng)力敏感系數(shù)和啟動(dòng)壓力梯度組合下生產(chǎn)10000d時(shí)的日產(chǎn)油量和累計(jì)產(chǎn)油量。從圖10可以看出，同時(shí)考慮應(yīng)力敏感和啟動(dòng)壓力梯度與不考慮或只考慮應(yīng)力敏感(啟動(dòng)壓力梯度)相比，頁(yè)巖油井產(chǎn)油量的遞減速度更快，當(dāng)應(yīng)力敏感系數(shù)和啟動(dòng)壓力梯度為0.5MPa-1和0.5MPa/m

圖10 考慮應(yīng)力敏感系數(shù)和啟動(dòng)壓力梯度時(shí)頁(yè)巖油井的產(chǎn)油量Fig.10 Shale oil production with consideration of stress sensitivity coefficient and threshold pressure gradients

時(shí)，生產(chǎn)10000d的累計(jì)產(chǎn)油量為不考慮應(yīng)力敏感和啟動(dòng)壓力梯度時(shí)的5.72%。由此可見，應(yīng)力敏感和啟動(dòng)壓力梯度對(duì)頁(yè)巖油井的產(chǎn)能具有很大的影響，不考慮應(yīng)力敏感和啟動(dòng)壓力梯度會(huì)極大地高估頁(yè)巖油井的產(chǎn)能，因此在模擬頁(yè)巖油藏開發(fā)產(chǎn)過(guò)程中，必須考慮應(yīng)力敏感系數(shù)和啟動(dòng)壓力梯度的影響。

5 結(jié) 論

1) 在考慮啟動(dòng)壓力梯度和應(yīng)力敏感效應(yīng)的基礎(chǔ)上，建立了頁(yè)巖油藏油、氣、水三相滲流數(shù)學(xué)模型，給出了模型的數(shù)值求解方法，并編制了頁(yè)巖油藏?cái)?shù)值模擬器。

2) 數(shù)值模擬結(jié)果表明：?jiǎn)?dòng)壓力梯度與應(yīng)力敏感效應(yīng)對(duì)頁(yè)巖油井的產(chǎn)能均有較大的抑制作用。啟動(dòng)壓力梯度減緩了壓力波的傳播，當(dāng)啟動(dòng)壓力梯度達(dá)到0.5MPa/m時(shí)，生產(chǎn)10000d時(shí)的累計(jì)產(chǎn)油量為無(wú)啟動(dòng)壓力梯度時(shí)的24.7%；應(yīng)力敏感效應(yīng)使開發(fā)中頁(yè)巖油儲(chǔ)層的滲透率降低，尤其是近井區(qū)域，當(dāng)應(yīng)力敏感系數(shù)由0.1MPa-1升至0.5MPa-1時(shí)，生產(chǎn)10000d時(shí)近井區(qū)域的滲透率降低近一個(gè)數(shù)量級(jí)，油井累計(jì)產(chǎn)油量為無(wú)應(yīng)力敏感時(shí)的36.2%；同時(shí)考慮啟動(dòng)壓力梯度和應(yīng)力敏感效應(yīng)時(shí)，頁(yè)巖油井的產(chǎn)能更低。

3) 頁(yè)巖油藏的滲流機(jī)理復(fù)雜，目前對(duì)頁(yè)巖油的滲流機(jī)理認(rèn)識(shí)不足。因此，下一步應(yīng)著重開展頁(yè)巖油藏的微觀滲流機(jī)理研究，完善頁(yè)巖油的滲流數(shù)學(xué)模型，從而指導(dǎo)頁(yè)巖油藏的開發(fā)。

References

[1] 盛湘，陳祥，章新文，等.中國(guó)陸相頁(yè)巖油開發(fā)前景與挑戰(zhàn)[J].石油實(shí)驗(yàn)地質(zhì)，2015，37(3)：267-271.

SHENG Xiang,CHEN Xiang,ZHANG Xinwen,et al. Prospects and challenges of continental shale oil development in China[J]. Petroleum Geology and Experiment,2015,37(3):267-271.

[2] 姚軍，孫致學(xué)，張凱，等.非常規(guī)油氣藏開采中的工程科學(xué)問題及其發(fā)展趨勢(shì)[J].石油科學(xué)通報(bào)，2016，1(1)：128-142.

YAO Jun,SUN Zhixue,ZHANG Kai,et al. Scientific engineering problems and development trends in unconventional oil and gas reservoirs[J]. Petroleum Science Bulletin,2016,1(1):128-142.

[3] 高英，朱維耀，岳明，等.體積壓裂頁(yè)巖油儲(chǔ)層滲流規(guī)律及產(chǎn)能模型[J].東北石油大學(xué)學(xué)報(bào)，2015，39(1)：80-86.

GAO Ying,ZHU Weiyao,YUE Ming,et al. Seepage law and productivity model of stimulated reservoir volume in shale oil reservoirs[J]. Journal of Northeast Petroleum University,2015,39(1):80-86.

[4] 郭肖，任影，吳紅琴.考慮應(yīng)力敏感和吸附的頁(yè)巖表觀滲透率模型[J].巖性油氣藏，2015，27(4)：109-112，118.

GUO Xiao,REN Ying,WU Hongqin. Apparent permeability model of shale gas considering stress sensitivity and adsorption[J]. Lithologic Reservoirs,2015,27(4):109-112,118.

[5] CHHARTE S S,SINHA S,BRAUN E M,et al. Effect ofstress,creep,and fluid type on steady state permeability measurements in tight liquid unconventional reservoirs[C]//Unconventional Resources Technology Conference,August25-27,Denver,Colorado,USA,2014.

[6] 孫海，姚軍，孫致學(xué)，等.頁(yè)巖氣數(shù)值模擬技術(shù)進(jìn)展及展望[J].油氣地質(zhì)與采收率，2012，19(1)：46-49.

SUN Hai,YAO Jun,SUN Zhixue,et al. Recent development and prospect on numerical simulation of shale gas reservoirs[J]. Petroleum Geology and Recovery Efficiency,2012,19(1):46-49.

[7] 嚴(yán)俠，黃朝琴，姚軍，等.裂縫性油藏改進(jìn)多重子區(qū)域模型[J].中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版)，2016，40(3)：121-129.

YAN Xia,HUANG Zhaoqin,YAO Jun,et al. An improved multiple sub-region method for flow simulation of fractured reservoirs[J]. Journal of China University of Petroleum (Edition of Natural Science),2016,40(3):121-129.

[8] PRUESS K,NARASIMHAN T N. Practical method for modeling fluid and heat flow in fractured porous media[J]. Society of Petroleum Engineers Journal,1985,25(1):14-26.

[9] LEE S H,LOUGH M F,JENSEN C L. Hierarchical modeling of flow in naturally fractured formations with multiple length scales[J]. Water Resources Research,2001,37(3):443-455.

[10] LI Liyong,LEE S H. Efficient field-scale simulation of black oil in a naturally fractured reservoir through discrete fracture networks and homogenized media[J]. SPE Reservoir Evaluation & Engineering,2008,11(4):750-758.

[11] 田冷，何順利，李秀生.低滲透氣田砂巖儲(chǔ)層應(yīng)力敏感試井模型研究[J].石油鉆探技術(shù)，2007，35(6)：89-92.

TIAN Leng,HE Shunli,LI Xiusheng. Study of well test of stress-sensitive sandstone in low permeability gas reservoir[J]. Petroleum Drilling Techniques,2007,35(6):89-92.

[12] 王敬，劉慧卿，劉仁靜，等.考慮啟動(dòng)壓力和應(yīng)力敏感效應(yīng)的低滲、特低滲油藏?cái)?shù)值模擬研究[J].巖石力學(xué)與工程學(xué)報(bào)，2013,32(增刊2)：3317-3327.

WANG Jing,LIU Huiqing,LIU Renjing,et al. Numerical simulation for low-permeability and extra-low permeability reservoirs with considering starting pressure and stress sensitivity effects[J]. Chinese Journal of Rock Mechanics and Engineering,2013,32(supplement2):3317-3327.

[13] 曲占慶，翟恒立，田相雷，等.考慮壓敏效應(yīng)的變啟動(dòng)壓力梯度試驗(yàn)研究[J].石油鉆探技術(shù)，2012，40(3)：78-82.

QU Zhanqing,ZHAI Hengli,TIAN Xianglei,et al. Experimental research on variable threshold pressure gradient considering pressure sensitive effect[J]. Petroleum Drilling Techniques,2012,40(3):78-82.

[14] 張慶輝，李相方，張磊，等.考慮啟動(dòng)壓力梯度的低滲底水氣藏見水時(shí)間預(yù)測(cè)[J].石油鉆探技術(shù)，2012，40(5)：96-99.

ZHANG Qinghui,LI Xiangfang,ZHANG Lei,et al. Prediction of water breakthrough time in low permeability gas reservoirs with bottom water in consideration of threshold pressure[J]. Petroleum Drilling Techniques,2012,40(5):96-99.

[15] 張楠，王曉琴，徐鋒，等.啟動(dòng)壓力梯度和應(yīng)力敏感效應(yīng)對(duì)低滲透油藏直井產(chǎn)能的影響[J].特種油氣藏，2012，19(1)：74-77.

ZHANG Nan,WANG Xiaoqin,XU Feng,et al. Effects of kickoff pressure gradient and stress sensitivity on the productivity of vertical wells in low permeability reservoirs[J]. Special Oil & Gas Reservoirs,2012,19(1):74-77.

[16] 楊仁鋒，姜瑞忠，劉世華，等.特低滲透油藏非線性滲流數(shù)值模擬[J].石油學(xué)報(bào)，2011,32(2)：299-306.

YANG Renfeng,JIANG Ruizhong,LIU Shihua,et al. Numerical simulation of nonlinear seepage in ultra-low permeability reservoirs[J]. Acta Petrolei Sinica,2011,32(2):299-306.

[17] 熊健，邱桃，郭平，等.非線性滲流下低滲氣藏壓裂井產(chǎn)能評(píng)價(jià)[J].石油鉆探技術(shù)，2012,40(3)：92-96.

XIONG Jian,QIU Tao,GUO Ping,et al. Production evaluation of fractured wells in low permeability reservoirs with nonlinear flow[J]. Petroleum Drilling Techniques,2012,40(3):92-96.

[18] 韓大匡，陳欽雷，閆存章.油藏?cái)?shù)值模擬基礎(chǔ)[M].北京：石油工業(yè)出版社，1993.

HAN Dakuang,CHEN Qinlei,YAN Cunzhang. Basis of reservoir numerical simulation[M]. Beijing:Petroleum Industry Press,1993.

[19] 張睿，寧正福，楊峰，等.微觀孔隙結(jié)構(gòu)對(duì)頁(yè)巖應(yīng)力敏感影響的實(shí)驗(yàn)研究[J].天然氣地球科學(xué)，2014，25(8)：1284-1289.

ZHANG Rui,NING Zhengfu,YANG Feng,et al. Experimental study on microscopic pore structure controls on shale permeability under compaction process[J]. Natural Gas Geoscience,2014,25(8):1284-1289.

[20] JU Binshan,WU Yushu,FAN Tailiang. Study on fluid flow in nonlinear elastic porous media:experimental and modeling approaches[J]. Journal of Petroleum Science and Engineering,2011,76(3/4):205-211.

[21] 楊瓊，聶孟喜，宋付權(quán).低滲透砂巖滲流啟動(dòng)壓力梯度[J].清華大學(xué)學(xué)報(bào)(自然科學(xué)版)，2004，44(12)：1650-1652.

YANG Qiong,NIE Mengxi,SONG Fuquan. Threshold pressure gradient of low permeability sandstone[J]. Journal of Tsinghua University(Science and Technology),2004,44(12):1650-1652.

[22] REYES L,OSISANYA S O. Empirical correlation of effective stress dependent shale rock properties[J]. Journal of Canadian Petroleum Technology,2002,41(12):47-53.

[編輯 劉文臣]

TheEffectofThresholdPressureGradientandStressSensitivityonShaleOilReservoirProductivity

LIULijun1,YAOJun1,SUNHai1,BAIYuhu2,XUBingxiang2,CHENLing2

(1.ResearchCentreofMultiphaseFlowinPorousMedia,ChinaUniversityofPetroleum(EastChina),Qingdao,Shandong,266580，China;2.CNOOCResearchInstitute,Beijing,100028,China)

In order to research on the effect of threshold pressure gradient and stress sensitivity on shale oil well productivity,a mathematical model that takes into consideration of the threshold pressure gradient and stress sensitivity for shale oil seepage was developed on the basis of oil-gas-water three-phase seepage model.A numerical solution method was offered and numerical simulator of shale oil reservoir was compiled.Patterns and mathematical rules for threshold pressure gradient and stress sensitivity on horizontal well productivity of multistage fracturing shale oil reservoir were analyzed by using a shale oil reservoir numerical simulator.The calculation results on the degradation of the simulator were basically identical with Eclipse software calculation results,which validated the correctness of the simulator.The numerical simulation results showed that the cumulative oil production rate for 10 000 days was 24.7% of that where there was no threshold pressure gradient when threshold pressure gradient reached 0.5 MPa/m. When the stress sensitivity coefficient rose from 0.1 MPa-1to 0.5 MPa-1,the reservoir permeability near borehole zones dropped approximately one order of magnitude and the cumulative oil production rate was 36.2% that of no stress sensitivity effect after 10 000 days’production.Moreover,the cumulative oil production was lower while simultaneously considering threshold pressure gradient and stress sensitivity effect.The results showed that the threshold pressure gradient and stress sensitivity effect both greatly inhibited shale oil productivity and the shale oil productivity might be overestimated without consideration of the threshold pressure gradient and stress sensitivity effect.Therefore a shale oil reservoir numerical simulator can provide guidance for the development of shale oil reservoirs.

shale oil;threshold pressure gradient;stress sensitivity;productivity

TE32+8

A

1001-0890(2017)05-0084-08

10.11911/syztjs.201705015

2017-05-19；改回日期2017-09-07。

劉禮軍(1993—)，男，山東海陽(yáng)人，2016年畢業(yè)于中國(guó)石油大學(xué)(華東)石油工程專業(yè)，在讀博士研究生，主要從事油氣田開發(fā)工程技術(shù)研究工作。E-mail：liulijunupc@163.com。

國(guó)家自然科學(xué)基金項(xiàng)目“頁(yè)巖油氣多尺度滲流特征與開采理論”(編號(hào)：51490654)、中海石油(中國(guó))有限公司綜合科研項(xiàng)目“海外頁(yè)巖油氣產(chǎn)能評(píng)價(jià)技術(shù)與方法研究”(編號(hào)：YXKY-2016-ZY-03)、中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金“頁(yè)巖油儲(chǔ)層多尺度數(shù)字巖心構(gòu)建及流動(dòng)機(jī)制研究”(編號(hào)：17CX02008A)及國(guó)家科技重大專項(xiàng)“基于數(shù)字巖心的頁(yè)巖氣藏微觀流動(dòng)模擬及產(chǎn)能評(píng)價(jià)方法研究”(編號(hào)：2016ZX05061-14)聯(lián)合資助。