于倩男*，劉義坤，劉 學(xué)，姚 迪，于 洋

1.東北石油大學(xué)石油工程學(xué)院，黑龍江 大慶 163318；2.“提高油氣采收率”教育部重點(diǎn)實(shí)驗(yàn)室·東北石油大學(xué)，黑龍江 大慶 163318；3.中國(guó)石油大慶油田有限責(zé)任公司，黑龍江 大慶 163001

引 言

隨著中國(guó)油氣開(kāi)發(fā)進(jìn)程的深入和相關(guān)理論與技術(shù)的發(fā)展，高含水油田水驅(qū)開(kāi)發(fā)的主要對(duì)象已由物性較好的均質(zhì)儲(chǔ)層轉(zhuǎn)變?yōu)榉蔷|(zhì)性較強(qiáng)的低滲透、特低滲透儲(chǔ)層。當(dāng)前開(kāi)發(fā)條件下非均質(zhì)低滲透油藏的整體動(dòng)用狀況仍有待提高，因此，流體在非均質(zhì)低滲透油藏中的滲流特征是需要解決的基礎(chǔ)問(wèn)題。滲流特征是承接滲流機(jī)理與生產(chǎn)狀況的重要中間過(guò)程，既是滲流機(jī)理具象直觀的表征，又將直接影響生產(chǎn)動(dòng)態(tài)。不同非均質(zhì)狀況低滲透儲(chǔ)層滲流特征的定性描述及定量表征可為深入研究滲流機(jī)理并指導(dǎo)生產(chǎn)實(shí)踐提供必要的研究基礎(chǔ)。

目前多孔介質(zhì)中流體滲流特征的相關(guān)研究已有較多，但非均質(zhì)低滲透儲(chǔ)層的滲流特征及機(jī)理的研究還不夠系統(tǒng)和完善，滲流特征測(cè)試中的實(shí)驗(yàn)參數(shù)多限于滲流壓差流量關(guān)系[1-5]或天然巖芯相對(duì)滲透率曲線的測(cè)定[6-8],而表征平面滲流特征的壓力梯度分布情況這一重要參數(shù)的相關(guān)文獻(xiàn)并不多見(jiàn)。在物理模擬實(shí)驗(yàn)手段方面，常規(guī)一維物理模型應(yīng)用廣泛，具有成本低廉和測(cè)試方便的優(yōu)勢(shì)，但其僅可模擬一維流動(dòng)，無(wú)法模擬不同注采關(guān)系的影響，因此在滲流機(jī)理角度存在缺陷，且與礦場(chǎng)實(shí)踐認(rèn)識(shí)差異較大[9-11]。三維的平板模型在壓力、壓力梯度及飽和度的分布狀況測(cè)試中更具優(yōu)勢(shì)，并在單相流體滲流、儲(chǔ)層應(yīng)力狀況測(cè)定和井網(wǎng)適應(yīng)性評(píng)價(jià)等實(shí)驗(yàn)中得到了較多的應(yīng)用[12-14]。非均質(zhì)低滲透平板模型可在實(shí)現(xiàn)儲(chǔ)層整體物性差模擬的同時(shí)表征注采關(guān)系及非均質(zhì)特性的影響，更好地模擬非均質(zhì)低滲透儲(chǔ)層的滲流特征。

本實(shí)驗(yàn)在明確平板模型相似理論基礎(chǔ)上，設(shè)計(jì)制作非均質(zhì)低滲透平板模型，并對(duì)其進(jìn)行評(píng)價(jià)，利用平板模型內(nèi)部對(duì)稱布置的壓力傳感器獲取的模型平面上的壓力數(shù)據(jù)，繪制滲流特征圖件，分析不同性質(zhì)非均質(zhì)低滲透儲(chǔ)層的滲流特征。

1 實(shí) 驗(yàn)

1.1 實(shí)驗(yàn)原理

相似理論是非均質(zhì)低滲透平板模型滲流特征實(shí)驗(yàn)的基礎(chǔ)原理，幾何相似、物性相似和制度相似是必要的相似條件。幾何相似指幾何尺寸等比例、形狀相似及注采關(guān)系的對(duì)應(yīng)；物性相似指滲透性等基本物性特征的相似；制度相似指邊界屬性和注采制度的相似。

平板模型物理模擬實(shí)驗(yàn)中，模型尺寸和注采參數(shù)等物理量應(yīng)依據(jù)相似條件縮小到合理尺度[15-17],將壓力、長(zhǎng)度、時(shí)間和其他流體及儲(chǔ)層性質(zhì)等參數(shù)的對(duì)應(yīng)比例關(guān)系帶入油水兩相滲流方程中，簡(jiǎn)化得到平板模型及模擬儲(chǔ)層中各物理量的相似比例尺。

由此可知，在對(duì)應(yīng)的時(shí)間節(jié)點(diǎn)上和對(duì)應(yīng)位置處模型和油藏中的和等4個(gè)組合物理量應(yīng)相等，才可以保證模型對(duì)油藏的滲流特征模擬的精確性。

根據(jù)平板模型的相似理論，由典型非均質(zhì)低滲透注采單元的參數(shù)，得到平板模型與注采單元參數(shù)對(duì)比表（表1，表2），據(jù)此設(shè)計(jì)平板模型并確定實(shí)驗(yàn)參數(shù)。由于模具尺寸標(biāo)準(zhǔn)及模型強(qiáng)度要求，平板模型厚度無(wú)法嚴(yán)格根據(jù)相似比例縮小。

表1 典型非均質(zhì)低滲透注采單元參數(shù)表Tab.1 Parameters of typical heterogeneous low-permeability reservoir

表2 平板模型參數(shù)表Tab.2 Parameters of plate model

1.2 平板模型

平板模型的設(shè)計(jì)、制作與評(píng)價(jià)是非均質(zhì)低滲透儲(chǔ)層滲流特征實(shí)驗(yàn)的基礎(chǔ)環(huán)節(jié)，以低滲透儲(chǔ)層滲透性特征為基準(zhǔn)，平板模型設(shè)計(jì)及制備需同時(shí)考慮平面非均質(zhì)特征和注采關(guān)系。

1.2.1 平板模型的設(shè)計(jì)

通常情況下，非均質(zhì)儲(chǔ)層中注采井間沉積相變化較大，根據(jù)注采端控制面積相等的原則，將平板模型分為填砂區(qū)域來(lái)模擬儲(chǔ)層的非均質(zhì)狀況。

壓力數(shù)據(jù)采集所需壓力傳感器布置如圖1所示。

平板模型的設(shè)計(jì)與制備環(huán)節(jié)中測(cè)壓點(diǎn)的設(shè)置不可過(guò)多，過(guò)多則影響平板模型整體承壓性，并擾動(dòng)流體的滲流狀態(tài)，過(guò)少則無(wú)法獲取足夠的特征值來(lái)研究平板模型內(nèi)部壓力分布規(guī)律?？紤]到獲取盡可能多的特征壓力值的同時(shí)減少其對(duì)滲流特征的影響，12個(gè)測(cè)壓點(diǎn)對(duì)稱分布在平板模型上，縱向上壓力傳感布置在模型中部位置。

圖1 非均質(zhì)低滲透平板模型設(shè)計(jì)圖Fig.1 Design of heterogeneous low-permeability artificial core plate models

1.2.2 平板模型的制備

平板模型的制備包括物料混合、石英砂填充、膠結(jié)物配制、模型成型和固化膠結(jié)等[18-19]。模型滲透率可通過(guò)調(diào)節(jié)石英砂粒度和膠結(jié)物含量控制，孔隙度則利用膠結(jié)物含量控制，孔隙結(jié)構(gòu)通過(guò)添加黏土礦物和天然巖芯碎屑模擬，進(jìn)而更精確地模擬非均質(zhì)低滲透儲(chǔ)層的物性特征。

基礎(chǔ)測(cè)試表明，依照設(shè)計(jì)方案及工藝流程制備得到的非均質(zhì)低滲透平板模型（圖2）的尺寸規(guī)格及密封性符合設(shè)計(jì)標(biāo)準(zhǔn)并滿足實(shí)驗(yàn)要求。

圖2 非均質(zhì)低滲透平板模型俯視圖及正視圖Fig.2 Top view and front view of some heterogeneous low-permeability artificial core plate model

1.2.3 平板模型的評(píng)價(jià)

不同壓力狀況下低滲透地層滲透性相應(yīng)發(fā)生變化，滲透率相似前提下，滲透性隨壓力的變化趨勢(shì)相似才能保證平板模型模擬的準(zhǔn)確性。

借鑒五點(diǎn)法井網(wǎng)達(dá)西滲流的產(chǎn)量公式，利用平板模型驅(qū)替實(shí)驗(yàn)的流量、流體性質(zhì)及地層參數(shù)，可推導(dǎo)得到表征模型整體允許流體通過(guò)的能力的滲透性參數(shù)，將此參數(shù)作為平板模型的等效滲透率

典型非均質(zhì)低滲透平板模型驅(qū)替實(shí)驗(yàn)結(jié)果表明，平板模型等效滲透率與注采壓差呈正相關(guān)關(guān)系。由等效滲透率與注采壓差關(guān)系曲線（圖3）可知，隨注采壓差的增加，等效滲透率呈急劇升高到緩慢升高最終逐漸平緩的變化趨勢(shì)。

圖3 平板模型及對(duì)照巖芯的滲透率與注采壓差關(guān)系曲線Fig.3 Relationship curves between equivalent permeability and pressure difference of artificial core plate models and contrast cores from low-permeability reservior

選取滲透率與平板模型等效滲透率相近的低滲透儲(chǔ)層對(duì)照巖芯開(kāi)展對(duì)照實(shí)驗(yàn)測(cè)試，結(jié)果（圖3）表明，平板模型及對(duì)照巖芯中滲透率隨注采壓差的變化趨勢(shì)相似。

1.3 實(shí)驗(yàn)條件

1.3.1 實(shí)驗(yàn)設(shè)備

實(shí)驗(yàn)采用自主設(shè)計(jì)的平板模型實(shí)驗(yàn)裝置（圖4），裝置由注入控制系統(tǒng)、平板模型、壓力測(cè)量系統(tǒng)和采出計(jì)量系統(tǒng)4部分組成。

圖4 自主設(shè)計(jì)平板模型實(shí)驗(yàn)裝置Fig.4 The self-designed equipments for experiments on artificial core plate models

注入控制系統(tǒng)可提供低至0.001 MPa的穩(wěn)定壓力，系統(tǒng)包括高壓氣泵、穩(wěn)壓裝置和注入流體罐；壓力測(cè)量記錄系統(tǒng)由壓力巡檢儀和記錄系統(tǒng)組成，用于平板模型實(shí)驗(yàn)中壓力數(shù)據(jù)的獲??；采出計(jì)量系統(tǒng)由微流量計(jì)和電子天平組成，高精度微流量計(jì)和電子天平可實(shí)現(xiàn)采出流體的精確計(jì)量。

圖5 平板模型真空飽和裝置示意圖Fig.5 The vacuum and saturation equipments for artificial core plate models

非均質(zhì)低滲透平板模型體積較大，室內(nèi)物理模擬用于天然巖芯和一維人造巖芯的流體真空及飽和的常規(guī)方法已不再適用。利用平板模型設(shè)計(jì)中預(yù)留多孔的特點(diǎn)，改進(jìn)裝置并形成多點(diǎn)抽真空及飽和方法處理平板模型。平板模型封裝后注采端及預(yù)留測(cè)壓孔皆與外界連通（圖5），利用與外界連通多點(diǎn)抽真空，壓力表實(shí)時(shí)監(jiān)測(cè)模型不同位置的真空度變化。平板模型達(dá)到完全真空后，將非抽真空點(diǎn)連接流體罐，繼續(xù)飽和實(shí)驗(yàn)流體，當(dāng)全部壓力表數(shù)值等于大氣壓時(shí)平板模型完全飽和。

1.3.2 實(shí)驗(yàn)條件

為了探討不同非均質(zhì)條件低滲儲(chǔ)層的滲流特征，設(shè)計(jì)制備符合不同研究需求的注采單元平板模型，測(cè)試平板模型的的壓力分布特征。注入端滲透率設(shè)置為10 mD，控制采出端滲透率在10～100 mD，重點(diǎn)分析注采端滲透率極差為1、2、5和10的4種非均質(zhì)情況。安裝調(diào)試實(shí)驗(yàn)裝置，利用煤油處理地層油，過(guò)濾地層水，模擬定壓生產(chǎn)情況，恒定壓差0.35 MPa下注入地層水。為獲得注采過(guò)程中注采單元控制面積內(nèi)壓力傳播與分布狀態(tài)，實(shí)時(shí)測(cè)量平板模型各點(diǎn)的壓力值，并記錄采出端流量數(shù)據(jù)。當(dāng)采出端流量穩(wěn)定后，獲取平板模型各點(diǎn)穩(wěn)定狀態(tài)的壓力值，之后依照實(shí)驗(yàn)方案調(diào)整注采壓差或更換平板模型，重復(fù)實(shí)驗(yàn)流程。利用實(shí)驗(yàn)獲取的壓力值和流量等實(shí)驗(yàn)數(shù)據(jù)，繪制相應(yīng)圖件，分析平板模型的壓力梯度分布狀況，進(jìn)而研究非均質(zhì)低滲透儲(chǔ)層的滲流特征。

2 實(shí)驗(yàn)結(jié)果與討論

平板模型的各壓力測(cè)試點(diǎn)未在平板模型內(nèi)部均勻分布，造成部分未監(jiān)測(cè)部位數(shù)據(jù)的缺失，將平板模型內(nèi)部平面網(wǎng)格劃分后，依據(jù)測(cè)試點(diǎn)壓力值利用外推內(nèi)插算法填充數(shù)據(jù)，繪制不同平板模型的壓力梯度分布圖和滲流區(qū)域劃分圖件，進(jìn)而分析不同性質(zhì)非均質(zhì)低滲透儲(chǔ)層的滲流特征。

2.1 模型平面壓力梯度分布特征

儲(chǔ)層流體滲流是消耗能量使流體流動(dòng)的過(guò)程，即消耗壓力獲得流量的過(guò)程。壓力梯度是沿壓力傳播及流體流動(dòng)方向單位長(zhǎng)度上的壓力變化，非均質(zhì)低滲透平板模型的壓力梯度分布狀況可反映非均質(zhì)低滲透儲(chǔ)層的滲流特征。

利用代號(hào)為10-10 mD、10-20 mD、10-50 mD和10-100 mD的4組平板模型進(jìn)行滲流特征物理模擬實(shí)驗(yàn)，表征注入端滲透率控制為10 mD，采出端滲透率為10、20、50和100 mD的4種平面非均質(zhì)變化狀況。利用實(shí)驗(yàn)測(cè)試的壓力數(shù)據(jù)，求取壓力梯度分布數(shù)據(jù)，并繪制非均質(zhì)低滲透平板模型的壓力梯度分布圖，進(jìn)而定性分析非均質(zhì)低滲透儲(chǔ)層的滲流特征。為方便對(duì)比分析，不同圖件采用相同的色度標(biāo)尺。

由平板模型壓力梯度場(chǎng)分布圖（圖6）可知，不同非均質(zhì)狀況的平板模型注入端與采出端壓力梯度分布狀況相似，壓力梯度等值線形態(tài)為近似橢圓形，橢圓形圓心為注入點(diǎn)或采出點(diǎn)，橢圓形長(zhǎng)軸為注采端的連線。平板模型中注采點(diǎn)周圍的壓力梯度非常高，則注采單元中近井地帶的壓力消耗很大；平板模型中注采對(duì)角處壓力梯度最小，則注采單元中同號(hào)井連線中點(diǎn)處壓力變化幅度最小。

圖6 平板模型壓力梯度場(chǎng)分布圖Fig.6 Pressure gradient distribution maps of artificial core plate models

平板模型的采出端滲透性增強(qiáng)，注采端平面的非均質(zhì)程度加劇。結(jié)果表明，注入點(diǎn)周圍相同傳播距離的壓降增大，壓力梯度增大，且壓力梯度變化幅度隨非均質(zhì)性增強(qiáng)而增大。由于采出端滲透率的增加，采出點(diǎn)處等距離壓降減小，壓力梯度有所降低，壓力梯度等值線分布更加稀疏。此時(shí)，平板模型注采對(duì)角處的低壓力梯度區(qū)域形態(tài)有所改變，區(qū)域面積變小。即在非均質(zhì)低滲透注采單元中，注入井點(diǎn)處滲透率不變，采出井點(diǎn)滲透率增加，導(dǎo)致注入井區(qū)域壓力梯度增加，采出井區(qū)域壓力梯度減小，同號(hào)井連線中點(diǎn)處的低壓力梯度區(qū)域面積減小。

2.2 注采連線壓力梯度分布特征

非均質(zhì)平板模型注采井間連線的壓力梯度示意圖（圖7）中，縱軸為壓力梯度（G），橫軸為無(wú)因次距離（L），指注采連線上的任意點(diǎn)距注入井的距離與注采井距的比值。

從圖7可知，整體上注采點(diǎn)附近壓力梯度較高，注采連線中點(diǎn)處壓力梯度最低，注采井連線中點(diǎn)附近一定距離內(nèi)壓力梯度穩(wěn)定在較低值；從注入端出發(fā)，注采連線上壓力梯度呈先降后升的變化趨勢(shì)。

統(tǒng)計(jì)注入端、注采連線中點(diǎn)處和采出端的壓力梯度（表3），不同非均質(zhì)低滲透平板模型的注入端及采出端壓力梯度變化幅度較大，注采連線中點(diǎn)處壓力梯度幾乎沒(méi)有變化；對(duì)于注入端，10-10 mD平板模型的壓力梯度值最小，為0.866 9 MPa/m，10-100 mD平板模型的壓力梯度值達(dá)到了0.935 2 MPa/m，增加幅度為9.95%；采出端10-10 mD平板模型的壓力梯度值最大為0.866 4 MPa/m，隨滲透率的變化，壓力梯度逐漸大幅降低為0.825 8、0.756 2和0.693 6 MPa/m，下降幅度分別為5.01%、13.02%和20.22%；注采連線中點(diǎn)處壓力梯度在0.095 0 MPa/m左右，幾乎沒(méi)有變化。

圖7 注采連線壓力梯度示意圖Fig.7 Pressure gradient distribution map of connection line of injection and production points

表3 注采連線壓力梯度數(shù)據(jù)表Tab.3 Datasheet of pressure gradient of connection line between injection and production points

結(jié)合模型平面整體與注采連線上壓力分布分析可知，注采端物性的變化決定了平板模型的非均質(zhì)性，二者對(duì)平板模型的壓力梯度分布有一定影響，表現(xiàn)為采出端壓力梯度降低，而注入端壓力梯度則有所增加。注入端滲透率穩(wěn)定在較低水平時(shí)，采出端滲透率越高，平板模型注采端的非均質(zhì)性越強(qiáng)，隨滲透率變大物性變好，區(qū)域內(nèi)壓力梯度降低，等壓力的傳播距離增大，同時(shí)非均質(zhì)性增強(qiáng)，對(duì)壓力傳播有負(fù)面影響，流體流動(dòng)所需能量衰減變大，井筒周圍壓力梯度升高，單位距離壓降有所增加。

2.3 滲流特征區(qū)域的劃分

2.3.1 滲流特征區(qū)域劃分的基本原理

利用典型低滲透儲(chǔ)層取芯天然巖樣的滲流實(shí)驗(yàn)結(jié)果，繪制的滲流速度與壓力梯度關(guān)系曲線即為典型低滲透儲(chǔ)層滲流特征曲線（圖8）。

當(dāng)注采壓力梯度小于低滲透儲(chǔ)層啟動(dòng)壓力梯度時(shí)，低滲透儲(chǔ)層不啟動(dòng)，此時(shí)滲流速度為零；當(dāng)注采壓力梯度大于啟動(dòng)壓力梯度且小于臨界啟動(dòng)壓力梯度時(shí)，低滲透儲(chǔ)層啟動(dòng)，且流體滲流速度非線性增加；當(dāng)注采壓力梯度大于臨界啟動(dòng)壓力梯度時(shí)，滲流速度線性增加[20-21]。低滲儲(chǔ)層中的滲流速度關(guān)系為

圖8 典型低滲透儲(chǔ)層的滲流特征曲線Fig.8 Typical relation curves of seepage velocity and pressure gradient in low-permeability reservoir

利用根據(jù)低滲透儲(chǔ)層取芯天然巖樣的滲流實(shí)驗(yàn)中非線性流動(dòng)段端點(diǎn)的驅(qū)替壓力梯度、滲透性參數(shù)和注入流體黏度數(shù)據(jù)，繪制低滲透儲(chǔ)層的滲流特征劃分圖版（圖9）。

低滲透儲(chǔ)層中壓力梯度分布情況的差異導(dǎo)致某黏度的流體在儲(chǔ)層中有不流動(dòng)、非線性滲流和擬線性滲流等3種流態(tài)[22-23]，即可將平板模型劃分為不流動(dòng)區(qū)、非線性滲流區(qū)和擬線性滲流區(qū)等3個(gè)滲流特征區(qū)域。

2.3.2 平板模型滲流特征區(qū)域的劃分

根據(jù)平板模型滲透率和實(shí)驗(yàn)注入流體黏度特征，參照低滲透儲(chǔ)層的滲流特征劃分圖版，得到不流動(dòng)區(qū)、非線性滲流區(qū)和擬線性滲流區(qū)等3個(gè)滲流特征區(qū)域劃分的壓力梯度值為0.047 0和0.193 0 MPa/m，進(jìn)而可將10-10 mD、10-20 mD、10-50 mD和10-100 mD不同非均質(zhì)性低滲透平板模型劃分為不同的滲流區(qū)域。

圖9 典型低滲透儲(chǔ)層滲流特征曲線Fig.9 The typical seepage flow pattern chart of low-permeability reservoir

圖10 平板模型滲流區(qū)域劃分圖Fig.10 Seepage flow section maps of artificial core plate models

由低滲透非均質(zhì)平板模型滲流區(qū)域劃分示意圖（圖10）可知，注采端相對(duì)的平板模型對(duì)角處存在不流動(dòng)區(qū)；整體上低滲透平板模型中非線性流動(dòng)區(qū)域面積占比最大，不流動(dòng)域面積占比最??；隨著采出端滲透率提高，平板模型注采端非均質(zhì)性增強(qiáng)，各滲流特征區(qū)域形態(tài)及面積占比相應(yīng)變化。

非均質(zhì)低滲透平板模型滲流區(qū)域劃分?jǐn)?shù)據(jù)表（表4）給出了不同滲流區(qū)域面積占比，及相比均質(zhì)狀況面積占比的變化幅度。10-10 mD平板模型的不流動(dòng)區(qū)域面積最大，面積占比達(dá)到了15.749%，隨著采出端滲透率增加，不流動(dòng)區(qū)的面積逐漸減小，10-100 mD平板模型的不流動(dòng)區(qū)域面積占比最低為14.361%，下降幅度為8.813%。相應(yīng)地，平板模型中可流動(dòng)區(qū)域面積（非線性滲流區(qū)和擬線性滲流區(qū)）整體上有所增加，其中，非線性滲流區(qū)域面積占比由50.707%逐漸減小至50.292%、49.668%、49.504%，同時(shí)，平板模型中更有利于流體流動(dòng)的擬線性滲流區(qū)域面積占比則由33.544%逐漸增加至36.135%。

采出端滲透率增加，平板模型的滲流能力增強(qiáng)，同時(shí)，非均質(zhì)程度變大對(duì)壓力傳播及流體流動(dòng)則有負(fù)面影響。分析非均質(zhì)低滲透平板模型滲流區(qū)域劃分?jǐn)?shù)據(jù)可知，由于實(shí)驗(yàn)采用非均質(zhì)低滲透平板模型的整體滲透率較低，采出端物性變好的影響程度大于非均質(zhì)的負(fù)面影響。

表4 平板模型滲流區(qū)域數(shù)據(jù)表Tab.4 Datasheet of seepage flow sections in artificial core plate models

3 結(jié) 論

（1）平板模型的注入端和采出端壓力梯度分布狀況相似，壓力梯度等值線形態(tài)相似；整體上壓力梯度等值線以井筒為圓心，以注采端連線為長(zhǎng)軸，呈近似橢圓形分布；注采點(diǎn)周圍壓力梯度非常高，模型對(duì)角處壓力梯度最小。非均質(zhì)低滲透平板模型的壓力分布狀況受滲透率及平面非均質(zhì)性的影響，壓力梯度等值線的形態(tài)相應(yīng)發(fā)生變化。

（2）注入端滲透率保持不變，采出端滲透率增大時(shí)，注入端壓力梯度升高，注采兩端非均質(zhì)性的增強(qiáng)對(duì)壓力傳播有負(fù)面影響，流體流動(dòng)所需能量衰減變大，單位距離壓降增加；采出端區(qū)域內(nèi)壓力梯度有所降低，滲透率的增大使得壓力傳播范圍增大；當(dāng)儲(chǔ)層整體滲透率較低時(shí)，滲透率增大影響程度大于非均質(zhì)的負(fù)面影響。

（3）非均質(zhì)低滲透儲(chǔ)層可劃分為不流動(dòng)區(qū)、非線性滲流區(qū)和擬線性滲流區(qū)等3個(gè)滲流特征區(qū)域。注入端滲透率不變采出端滲透率增大，注采單元中不流動(dòng)區(qū)域面積占比有所下降，相應(yīng)的可流動(dòng)區(qū)域面積有所增加，其中更利于流體流動(dòng)的擬線性滲流區(qū)域面積占比增大。

符號(hào)說(shuō)明

pm—平板模型壓力，MPa；

pr—油藏壓力，MPa；

tm—平板模型模擬時(shí)間，d；

tr—油藏的生產(chǎn)時(shí)間，d；

Lm—平板模型尺寸（長(zhǎng)、寬、高），m；

Lr—油藏尺寸（長(zhǎng)、寬、高），m；

μm—平板模型中流體黏度，mPa·s；

μr—油藏中流體黏度，mPa·s；

ρm—平板模型中流體密度，g/cm3；

ρr—油藏中流體密度，g/cm3；

σm—平板模型界面張力，N/m；

σr—油藏界面張力，N/m；

θm—平板模型潤(rùn)濕角，（?）；

θr—油藏潤(rùn)濕角，（?）；

—平板模型等效滲透率，mD；

Q—采出量，m3；

h—地層厚度，m；

?p—注采壓差，MPa；

μ—流體黏度，mPa·s；

d—注采井距，m；

rw—井筒半徑，m；

υ—滲流速度，m/s；

a，b—地層及流體性質(zhì)相關(guān)的常數(shù)；

G—壓力梯度，MPa/m；

Gs—啟動(dòng)壓力梯度，MPa/m；

Gq—擬啟動(dòng)壓力梯度，MPa/m；

Gc—臨界啟動(dòng)壓力梯度，MPa/m。

[1]KORTEKAAS T.Water/oil displacement characteristics in crossbedded reservoir zones[J].Society of Petroleum Engineers Journal,1985,25(6):917–926.doi:10.2118/-12112-PA

[2]MALA G M,LI Dongqing.Flow characteristics of water in microtubes[J].International Journal of Heat and Fluid Flow,1999,20(2):142–148.doi:10.1016/S0142-727-X(98)10043-7

[3]曲志浩，孔令榮.低滲透油層微觀水驅(qū)油特征[J].西北大學(xué)學(xué)報(bào)（自然科學(xué)版），2002，32（4）：329-334.doi：10.3321/j.issn:1000-274X.2002.04.002 QU Zhihao,KONG Lingrong.A study of characteristics of low permeability formation with water flooding using micromodel experiment[J].Journal of Northwest University(Natural Science Edition),2002,32(4):329–334.doi:10.3321/j.issn:1000-274X.2002.04.002

[4]崔悅，石京平.大慶油田表外儲(chǔ)層滲流特征實(shí)驗(yàn)[J].大慶石油地質(zhì)與開(kāi)發(fā)，2017，36（2）：69-72.doi：10.3969/-j.issn.1000-3754.2017.02.011 CUI Yue,SHI Jingping.Experiment on the seepage features for the untabulated reservoirs in Daqing Oilfield[J].Petroleum Geology&Oilfield Development in Daqing,2017,36(2):69–72.doi:10.3969/j.issn.1000-3754.2017.-02.011

[5]李亞軍，李康寧，蘇玉亮，等.低滲多孔介質(zhì)單相液體穩(wěn)定流動(dòng)特征實(shí)驗(yàn)研究[J].科學(xué)技術(shù)與工程，2016，16（11）：49-53.doi：10.3969/j.issn.1671-1815.2016.11.-009 LI Yajun,LI Kangning,SU Yuliang,et al.Laboratory measurements of single-phase steady flow in lowpermeability rocks[J].Science Technology and Engineering,2016,16(11):49–53.doi:10.3969/j.issn.1671-1815.-2016.11.009

[6]張光明，廖銳全，徐永高，等.特低滲透油藏流體滲流特征試驗(yàn)研究[J].江漢石油學(xué)院學(xué)報(bào)，2004，26（1）：88-89，145.doi：10.3969/j.issn.1000-9752.2004.01.035 ZHANG Guangming,LIAO Ruiquan,XU Yonggao,et al.Experimental study of fluid percolating characteristics in extremely low permeable reservoirs[J].Journal of Jianghan Petroleum Institute,2004,26(1):88–89,145.doi:10.-3969/j.issn.1000-9752.2004.01.035

[7]李愛(ài)芬，劉敏，張少輝，等.特低滲透油藏滲流特征實(shí)驗(yàn)研究[J].西安石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2008，23（2）：35-39，117.doi：10.3969/j.issn.1673-064X.2008.02.009 LI Aifen,LIU Min,ZHANG Shaohui,et al.Experimental study on the percolation characteristic of extra lowpermeability reservoir[J].Journal of Xi′an Shiyou University(Natural Science Edition),2008,23(2):35–39,117.doi:10.3969/j.issn.1673-064X.2008.02.009

[8]XU Jianchun,GUO Chaohua,JIANG Ruizhong,et al.Study on relative permeability characteristics affected by displacement pressure gradient:Experimental study and numerical simulation[J].Fuel,2016,163(1):314–323.doi:10.1016/j.fuel.2015.09.049

[9]閆慶來(lái)，何秋軒，尉立崗，等.低滲透油層中單相液體滲流特征的實(shí)驗(yàn)研究[J].西安石油大學(xué)學(xué)報(bào)（自然科學(xué)版），1990，5（2）：1-6.YAN Qinglai,HE Qiuxuan,WEI Ligang,et al.A laboratory study on percolation characteristics of single phase flow in low-permeability reserviors[J].Journal of Xi′an Shiyou University(Natural Science Edition),1990,5(2):1–6.

[10]吳岱峰，魯曉兵，劉慶杰，等.低滲巖芯驅(qū)替過(guò)程中的壓力分布特征研究[J].力學(xué)與實(shí)踐，2012，34（5）：27-31，51.doi：10.6052/1000-0879-12-063 WU Daifeng,LU Xiaobing,LIU Qingjie,et al.Pressure distribution in flooding of low permeability reservoirs[J].Mechanics in Engineering,2012,34(5):27–31,51.doi:10.6052/1000-0879-12-063

[11]劉凱，文章，梁杏，等.一維低滲透介質(zhì)非達(dá)西滲流實(shí)驗(yàn)[J].水動(dòng)力學(xué)研究與進(jìn)展，2013，28（1）：81-87.doi：10.3969/j.issn1000-4874.2013.01.012 LIU Kai,WEN Zhang,LIANG Xing,et al.Onedimensional column test for non-darcy flow in low permeability media[J].Chinese Journal of Hydrodynamics,2013,28(1):81–87.doi:10.3969/j.issn1000-4874.2013.-01.012

[12]薛成國(guó).特低滲透儲(chǔ)層大型物理模擬實(shí)驗(yàn)研究[D].北京：中國(guó)科學(xué)院滲流流體力學(xué)研究所，2011.XUE Chengguo.Large physical simulation experimental study on the ultra-low permeability reservoir[D].Beijing:Institute of Porous Flow and Fluid Mechanics,Chinese Academy of Sciences,2011.

[13]徐軒，劉學(xué)偉，楊正明，等.特低滲透砂巖大型露頭模型單相滲流特征實(shí)驗(yàn)[J].石油學(xué)報(bào)，2012，33（3）：453-458.doi：10.7623/syxb201203016 XU Xuan,LIU Xuewei YANG Zhengming,et al.An experimental study on single-phase seepage characteristics with a large-scale model made of ultra-low permeability sandstone outcrops[J].Acta Petrolei Sinica,2012,33(3):453–458.doi:10.7623/syxb201203016

[14]于倩男.杏南開(kāi)發(fā)區(qū)精控壓裂薄差儲(chǔ)層滲流特征研究[D].大慶：東北石油大學(xué)，2017.YU Qiannan. Research on seepage flow patterns of fine controlled fractured thin and poor reservoir in Xingnan development aera[D]. Daqing: Northeast Petroleum University,2017.

[15]GEERTSMA J,CROES G A,SCHWARZ N.Theory of dimensionally scaled models of petroleum reservoirs[C].SPE-539-G,1955.

[16]孔祥言，陳峰磊，裴柏林.水驅(qū)油物理模擬理論和相似準(zhǔn)則[J].石油勘探與開(kāi)發(fā)，1997，24（6）：56-60.doi：10.3321/j.issn:1000-0747.1997.06.015 KONG Xiangyan.Similan theory,similar rule and data conversion software for water drive modeling[J].Petroleum Exploration and Development,1997,24(6):56–60.doi:10.3321/j.issn:1000-0747.1997.06.015

[17]滕起，楊正明，劉學(xué)偉，等.特低滲透油藏水驅(qū)油物理模擬相似準(zhǔn)則的推導(dǎo)和應(yīng)用[J].科技導(dǎo)報(bào)，2013，31（9）：40-45.TENG Qi,YANG Zhengming,LIU Xuewei,et al.Similar criteria derivation for the physical simulation of water flooding in the plate model of ultra-low permeability reservoir and its applications[J].Science&Technology Review,2013,31(9):40–45.

[18]皮彥夫.石英砂環(huán)氧樹(shù)脂膠結(jié)人造巖心的技術(shù)與應(yīng)用[J].科學(xué)技術(shù)與工程，2010，10（28）：6998-7000，7010.doi：10.3969/j.issn.1671-1815.2010.28.034 PI Yanfu.Technology and application of making artificial cores by the cementation of quartzite and colophony[J].Science Technology and Engineering,2010,10(28):6998–7000,7010.doi:10.3969/j.issn.1671-1815.2010.-28.034

[19]徐宏光，熊鈺，王永清，等.特高滲疏松砂巖人造巖心的制作及評(píng)價(jià)[J].石油鉆采工藝，2017，39（4）：477-483.doi：10.13639/j.odpt.2017.04.015 XU Hongguang,XIONG Yu,WANG Yongqing,et al.Preparation and evaluation on artificial core of extra high permeability unconsolidated sandstone[J].Oil Drilling&Production Technology,2017,39(4):477–483.doi:10.-13639/j.odpt.2017.04.015

[20]WANG Xiukun,Sheng J J.Discussion of liquid threshold pressure gradient[J].Petroleum,2017,3(2):232–236.doi:10.1016/j.petlm.2017.01.001

[21]呂成遠(yuǎn)，王建，孫志剛.低滲透砂巖油藏滲流啟動(dòng)壓力梯度實(shí)驗(yàn)研究[J].石油勘探與開(kāi)發(fā)，2002，29（2）：86-89.doi：10.3321/j.issn:1000-0747.2002.02.023 Lü Chengyuan,WANG Jian,SUN Zhigang.An experimental study on starting pressure gradient of fluids flow in low permeability sandstone porous media[J].Petroleum Exploration and Development,2002,29(2):86–89.doi:10.3321/j.issn:1000-0747.2002.02.023

[22]時(shí)佃海.低滲透砂巖油藏平面徑向滲流流態(tài)分布[J].石油勘探與開(kāi)發(fā)，2006，33（4）：491-494.doi：10.3321/j.-issn:1000-0747.2006.04.019 SHI Dianhai. Flow state distribution of areal radial flow in low permeability sandstone reservoir[J]. Petroleum Exploration and Development, 2006, 33(4): 491–494. doi: 10.-3321/j.issn:1000-0747.2006.04.019

[23]YU Qiannan,LIU Yikun,LIU Xue,et al.Experimental study on seepage flow patterns in heterogeneous lowpermeability reservoirs[J].Journal of Petroleum Exploration&Production Technology,17(6):1–8.doi:10.-1007/s13202-017-0354-y