高含水期油藏液流方向優(yōu)化及流線模擬

李才學(xué)，沈曦，賈衛(wèi)平，毛立華，朱黎明

（1.中國(guó)石化中原油田分公司采油六廠，山東　東明　274515；2.中國(guó)石化中原油田分公司勘探開(kāi)發(fā)研究院，河南　濮陽(yáng)　457001；3.中國(guó)石化中原油田分公司采油一廠，河南　濮陽(yáng)　457172）

高含水期油藏液流方向優(yōu)化及流線模擬

李才學(xué)1，沈曦2，賈衛(wèi)平3，毛立華2，朱黎明2

（1.中國(guó)石化中原油田分公司采油六廠，山東東明274515；2.中國(guó)石化中原油田分公司勘探開(kāi)發(fā)研究院，河南濮陽(yáng)457001；3.中國(guó)石化中原油田分公司采油一廠，河南濮陽(yáng)457172）

目前國(guó)內(nèi)大部分油田已經(jīng)進(jìn)入高含水開(kāi)發(fā)、剩余油零星分布階段，造成該現(xiàn)象的主要原因是由于儲(chǔ)層平面非均質(zhì)性、水驅(qū)波及程度的不同。為了提高水驅(qū)開(kāi)發(fā)效果，需要對(duì)液流方向優(yōu)化進(jìn)行研究，進(jìn)一步提高水驅(qū)采收率。該次研究以中原油田某高含水油藏為研究對(duì)象，將數(shù)值模擬和最優(yōu)理論相結(jié)合，建立了基于改進(jìn)的CRM模型為驅(qū)動(dòng)、EnOpt集合優(yōu)化為最優(yōu)化手段，以經(jīng)濟(jì)效益為目標(biāo)函數(shù)的液流方向自動(dòng)化優(yōu)化方法，從而量化油藏配產(chǎn)配注。以相控精細(xì)地質(zhì)模型為基礎(chǔ)，利用流線模擬器進(jìn)行開(kāi)發(fā)指標(biāo)預(yù)測(cè)，確定油井泄流面積、注水井水驅(qū)波及面積、水驅(qū)速度、井間儲(chǔ)層波及面積和井間儲(chǔ)層含水等參數(shù)。

液流方向優(yōu)化；流線模擬；配產(chǎn)配注

據(jù)統(tǒng)計(jì)，目前中國(guó)石化大部分油田已經(jīng)進(jìn)入了高含水率開(kāi)發(fā)階段，含水率大于90%的區(qū)塊地質(zhì)儲(chǔ)量占中國(guó)石化整個(gè)地質(zhì)儲(chǔ)量的43%，面臨遞減速度快、穩(wěn)產(chǎn)難度大的開(kāi)發(fā)難題。高含水油藏剩余油分布零散，下步開(kāi)采難度大，造成該現(xiàn)象的主要原因是儲(chǔ)層平面非均質(zhì)性、水驅(qū)波及程度不同，部分波及程度較低的區(qū)域成為剩余油富集區(qū)。如何提高油田水驅(qū)波及程度，制定配產(chǎn)配注的最優(yōu)化方案，成為提高油田水驅(qū)采收率的主要問(wèn)題［1-5］。

1　液流方向優(yōu)化模型建立及求解

新的液流優(yōu)化模型構(gòu)架見(jiàn)圖1，改進(jìn)的CRM模型為驅(qū)動(dòng)，EnOpt集合優(yōu)化為手段，凈現(xiàn)值最大為優(yōu)化方向。優(yōu)化目標(biāo)為優(yōu)化出分階段的最優(yōu)注入量、采液量［6-11］。

1.1改進(jìn)的CRM模型

原CRM模型下，油藏飽和度是常數(shù)，無(wú)法反映水驅(qū)后，含油飽和度變化對(duì)傳導(dǎo)率的影響。本方法創(chuàng)新之處在于，通過(guò)計(jì)算各時(shí)間步累注與累產(chǎn)數(shù)據(jù)，不斷修正模型各井間傳導(dǎo)率。

本質(zhì)上講，CRM是一個(gè)帶有壓縮系數(shù)的總質(zhì)量平衡方程。通常，每口生產(chǎn)井同時(shí)受效于多個(gè)注入井，依據(jù)每口生產(chǎn)井周圍的控制量，主要的物質(zhì)平衡可以由不同的方程表示為

式中：Ct為控制體積周圍的總壓縮系數(shù)，1/MPa；Vp為每口生產(chǎn)井周圍的孔隙體積，m3；pj為生產(chǎn)井采出液量的平均壓力，MPa；ik為每口注入井的注入速度，m3/d；qj為注入井的總生產(chǎn)速度，m3/d；λkj為注入井與生產(chǎn)井之間的連通系數(shù)；t為時(shí)間步數(shù)。

式中：qkj為注入井k對(duì)應(yīng)的生產(chǎn)井j的產(chǎn)液量，m3。

圖1　模型流程

依據(jù)儲(chǔ)層傳導(dǎo)性，λkj可解釋為

式中：Tkj為注入井k與生產(chǎn)井j之間的有效傳導(dǎo)系數(shù)。

有效傳導(dǎo)系數(shù)Tkj是儲(chǔ)層表征的一個(gè)因素，取值由儲(chǔ)層飽和點(diǎn)滲透率決定，因此，在本研究中，規(guī)定λkj是一個(gè)隨時(shí)間的變量。式（1）中，生產(chǎn)井總采液量受生產(chǎn)井周圍的平均壓力和總注入量劈分的注入量決定，但該狀態(tài)方程在假定注采井間壓縮系數(shù)很小以及其他井沒(méi)有流體流動(dòng)情況下才能成立。為了使方程完全依據(jù)流速，引入生產(chǎn)指數(shù)模型：

式中：pwf，j為生產(chǎn)井的井底壓力，MPa；Jj為生產(chǎn)井的采油指數(shù)，t/（MPa·d）。

忽略式（1）和式（4）中平均壓力，得到CRM的基本微分方程：

式中：τj為每口生產(chǎn)井驅(qū)替體積的時(shí)間常數(shù)。

同時(shí)，為λkj給定了一個(gè)m×n的參數(shù)，m是τj的參數(shù)，n是Jj的參數(shù)，可以表述為

式中：ξ為積分變量，其變化范圍為to～t。

1.1.1井間連通性與時(shí)間滯留隨時(shí)間的變化

有關(guān)CRM的前期研究都以井間連通性與時(shí)間滯留是時(shí)間的常數(shù)為前提。利用數(shù)值模擬方法研究后發(fā)現(xiàn)，均質(zhì)儲(chǔ)層中λkj不隨注入速度改變而改變。注采井之間有效傳導(dǎo)系數(shù)Tkj可以根據(jù)巖石滲透率、有效剖面、黏度（流體組分黏度之和）、注采井之間的距離、相對(duì)滲透率（流體各組分相對(duì)黏度之和）等地質(zhì)參數(shù)確定。

解析式（6），假定λkj，Tj，Jj是常數(shù)，在本研究中，利用帶有分段常數(shù)的動(dòng)力系統(tǒng)模擬了這些參數(shù)隨時(shí)間的可能變化。在這種情況下，式（6）可改寫(xiě)為

在式（7）中，假定τj（tn）和λkj（tn）在每個(gè)時(shí)間間隔（tn-1，tn）為常數(shù)，該解法導(dǎo)致一個(gè)時(shí)間離散，即在每個(gè)時(shí)間段，這些參數(shù)假定為常數(shù)，這可以使CRM的適用性擴(kuò)展至非成熟油藏，并能提供關(guān)于油藏更定量的信息。

1.1.2采油量模型

CRM模型表述了注采井之間總的流體動(dòng)力情況，在此模型中，不考慮油水之間的區(qū)別。不同現(xiàn)象模型可用來(lái)描述油井生產(chǎn)情況和液體特征之間的關(guān)系。在本研究中，應(yīng)用了一個(gè)冪律關(guān)系來(lái)描述在第j個(gè)生產(chǎn)井（Wj）的累計(jì)注水量與油水采出比之間的關(guān)系：

式中：aj為系數(shù)。

1.2整體為基礎(chǔ)的集合優(yōu)化（EnOpt模型）

以凈現(xiàn)值為目標(biāo)函數(shù)，應(yīng)用EnKF（卡爾曼濾波）方程進(jìn)行多參數(shù)自動(dòng)集合最優(yōu)化逐步迭代，得到最優(yōu)化的油水井注采量參數(shù)。

設(shè)y為該系統(tǒng)狀態(tài)變量。本次研究在每個(gè)時(shí)間點(diǎn)的y包含：動(dòng)態(tài)參數(shù)設(shè)置，每口井的總產(chǎn)量qj、每口井的累計(jì)注入量Wj、采油速度qo，j；靜態(tài)參數(shù)設(shè)置，｛λkj｝，｛Tj｝，｛Jj｝，｛aj｝，｛bj｝，以及可觀測(cè)的｛qj｝，｛qoj｝。在EnKF方程中，可觀測(cè)的參數(shù)包含在動(dòng)態(tài)參數(shù)之中，但作為可觀測(cè)的數(shù)據(jù)，也具分離的特點(diǎn)；因此，觀測(cè)參數(shù)可隨系統(tǒng)發(fā)展自動(dòng)更新，但也可通過(guò)單次投影所模擬的動(dòng)力變量得到。狀態(tài)向量變量的整體被存為Y=（y1，y2，y3，…，yNe）。ENKF程序提供了一個(gè)更新?tīng)顟B(tài)整體的方法，給出了觀測(cè)數(shù)dobs。

式中：dobs，j為每個(gè)整體的擾動(dòng)觀測(cè)量；yu為更新的整體；yf為通過(guò)數(shù)模預(yù)測(cè)的整體組分狀態(tài)；CD為測(cè)量噪音的協(xié)方差；CY為狀態(tài)向量變量的協(xié)方差矩陣；矩陣H為相態(tài)向量得到的觀測(cè)變量。

式中：〈Yf〉為狀態(tài)向量。

要優(yōu)化就必須存儲(chǔ)一個(gè)控制變量的集合，設(shè)X為描述控制變量的向量，X=［x1，x2，x3，…，xNx］，Nx為控制變量的總數(shù)，凈現(xiàn)值（NPV）是本次優(yōu)化的目標(biāo)函數(shù)，由控制向量X和油藏狀態(tài)向量Y表示。

式中：Nt為時(shí)間步總數(shù)；rΔt為時(shí)步Δt基礎(chǔ)上的折現(xiàn)率。

為確定年利率，可設(shè)定rΔt=exp（rΔt），其中r是年利率，時(shí)間步Δt是連續(xù)2個(gè)控制步之間的時(shí)間，控制步相當(dāng)于可以隨時(shí)間及時(shí)調(diào)整的注入速度這一點(diǎn)。例如，總開(kāi)發(fā)時(shí)間1000d，就可得到20個(gè)控制步，相當(dāng)于50d為一個(gè)時(shí)間間隔，這就可以根據(jù)優(yōu)化程序?qū)ψ⑷胨俣冗M(jìn)行調(diào)整。常數(shù)vo和vw分別為原油銷售和水處理的單價(jià)。Qo，t|x，y與Qw，t|x，y為所有井在控制變量x與狀態(tài)變量y所表述的t時(shí)間點(diǎn)，所有井的累計(jì)采油量和產(chǎn)水量。

利用求平均值的方式，在模型的不確定性中優(yōu)化控制變量。其公式表述如下：

參數(shù)y的每個(gè)選擇都可使目標(biāo)函數(shù)（x，yj）得到不同的值，下標(biāo)y表示NPV是依據(jù)于由EnKF所得的時(shí)變模型，模擬狀態(tài)向量Y在優(yōu)化部分設(shè)為不變量，NPV的優(yōu)化受控于控制變量x，在每一優(yōu)化步，更新的控制變量x1由式（14）得到：

式中：l為優(yōu)化迭代步；a為決定步長(zhǎng)的調(diào)優(yōu)參數(shù)；Cx，g（x）為控制變量與g（x）間的互協(xié)方差。

式中：xl，j為控制變量攝動(dòng)值的整體（控制變量是指具有平均值＜xl＞和前協(xié)方差的大小Ne）。

本研究用0平均值作為攝動(dòng)，對(duì)xl進(jìn)行高斯隨機(jī)域進(jìn)行校正。

1.3模擬分析

模擬油田包括5口注入井和4口生產(chǎn)井（見(jiàn)圖2）。油田大?。?5 m×15 m×2 m，網(wǎng)格數(shù)為31×31×5。油水流度比為1，油、水、巖石壓縮系數(shù)分別為5×10-6，1×10-6，1×10-6，注入井由注入速度控制，其井底壓力固定為0.1 MPa。儲(chǔ)層孔隙設(shè)定為24%，平均滲透率100×10-3μm2。儲(chǔ)層存在2個(gè)高滲透優(yōu)勢(shì)通道：注入井I1與生產(chǎn)井P1之間，優(yōu)勢(shì)通道滲透率1 000×10-3μm2；注入井I3與生產(chǎn)井P4之間，其滲透率500×10-3μm2。

圖2　概念模型

方案的擬合精度表明，EnKF-CRM的結(jié)合可以成功預(yù)測(cè)CRM在儲(chǔ)層研究中的未知參數(shù)。作為實(shí)際情況來(lái)看，EnOpt-CRM程序計(jì)算中，目標(biāo)函數(shù)的評(píng)價(jià)需要單獨(dú)的數(shù)值模擬來(lái)計(jì)算g（xl+1，yref），當(dāng)存在真實(shí)儲(chǔ)層的實(shí)際注入速度時(shí)，該計(jì)算由CRM來(lái)執(zhí)行。

針對(duì)不同優(yōu)化迭代步的注入速度與控制步之間關(guān)系的過(guò)程中，各個(gè)注入速度相關(guān)的變化趨勢(shì)是不同的，注入井I1和I3顯示的減低、增加趨勢(shì)與生產(chǎn)井P2，P4，P5剛好相反，I1和I3相對(duì)I2，I4，I5來(lái)說(shuō)，具有更低的注入速度，其優(yōu)化趨勢(shì)傾向于改善油藏驅(qū)替效率。I1和I3的注入速度低于其他井，這是因?yàn)镮1-1井組與I3-4井組內(nèi)更高的滲透率通道造成的。為了得到更高的凈現(xiàn)值，靠近高滲透帶井的注入速度需要降低，以使整個(gè)油藏可以得到更好的驅(qū)替。

2　流線模擬技術(shù)研究

以相控精細(xì)地質(zhì)模型為基礎(chǔ)，利用流線模擬器進(jìn)行開(kāi)發(fā)指標(biāo)預(yù)測(cè)，研究油井泄流面積、注水井水驅(qū)波及面積、水驅(qū)速度、井間儲(chǔ)層波及面積和井間儲(chǔ)層含水等參數(shù)，開(kāi)展高含水期油藏的注水利用效率研究［11-14］。

2.1配產(chǎn)配注量?jī)?yōu)化

應(yīng)用流線方法描述油水兩相驅(qū)替過(guò)程，要建立油水兩相流線數(shù)學(xué)模型。在黑油模型與流線方法相結(jié)合的基礎(chǔ)上，盡可能多地考慮地層流體滲流機(jī)理。

基本假設(shè)：1）忽略氣相的存在，流體為油、水兩相；2）油藏中流體的滲流是等溫的；3）油藏中流體和巖石均不可壓縮；4）油藏中流體的流動(dòng)符合達(dá)西滲流定律。

優(yōu)化流程可以歸納為以下步驟：1）利用流線法對(duì)油藏進(jìn)行數(shù)值模擬，并進(jìn)行歷史擬合，重點(diǎn)擬合現(xiàn)場(chǎng)井間受效情況，擬合油水井的連通情況；2）利用流線模擬器，計(jì)算每口油井連通的注水井，計(jì)算每口注水井注水量平面劈分量，并計(jì)算對(duì)周圍油井產(chǎn)油量的貢獻(xiàn)量；3）根據(jù)注入水利用率的計(jì)算方法計(jì)算每口注水井的注水利用率；4）根據(jù)注水井的注水利用率分配每口注水井的注水量，同時(shí)考慮注水井的注入能力，最大注水量和注水井的最大井底流壓；5）按照修改后的配產(chǎn)配注量進(jìn)行數(shù)值模擬，進(jìn)行下一步優(yōu)化，返回步驟2）。

2.2理論模型優(yōu)化

2.2.1流線優(yōu)化效果模擬

數(shù)值模擬的模型使用的是表征河流相沉積的非均質(zhì)模型，地質(zhì)模型表征多期河道相互疊置的模型，在平面上沿河道沉積的方法滲透率較高，與河道垂直方向滲透率較低。使用非均值模型在ECLIPSE平臺(tái)上進(jìn)行數(shù)值模擬，并研究?jī)?yōu)化方法對(duì)提高采收率的效果［15-17］。以五點(diǎn)法進(jìn)行部署井位，共13口井，其中4口注水井，9口生產(chǎn)井。由于采用的是非均質(zhì)性儲(chǔ)層，井間的連通性差別較大（見(jiàn)圖3、圖4）。如I2和P5之間的連通性較高，而P3和I2的連通性則較低。

根據(jù)水井的注水利用率，I3的注水利用率高，應(yīng)優(yōu)先考慮增加I3的注水量，其次是I4的注水利用率；I2的注水利用率在初期就出現(xiàn)了急劇下降的現(xiàn)象，因?yàn)镮2和P2，P5之間的井間連通程度高，P2和P5很快就見(jiàn)水，含水急劇上升，造成I2和P2，P5之間的無(wú)效注水循環(huán)。

圖3　地質(zhì)模型孔隙度分布

圖4　地質(zhì)模型滲透率分布

對(duì)比優(yōu)化前后的流場(chǎng)，優(yōu)化前井間存在多個(gè)水驅(qū)不到或水驅(qū)程度低的區(qū)域，如井I3和P5，I3和P8，P2和I1中間的儲(chǔ)層水驅(qū)程度低；而優(yōu)化之后這些井間的流線較優(yōu)化前密集，這些區(qū)域的剩余油飽和度較優(yōu)化之前低。優(yōu)化之后注入水的利用率有了很大的提升，說(shuō)明注入同樣的水，驅(qū)動(dòng)的油量更多，波及體積越大。對(duì)比兩次模擬結(jié)果，優(yōu)化后的采出程度有了很大的提升，含水也有較大幅度的降低（見(jiàn)圖5—8）。

圖5　優(yōu)化前含油飽和度分布

圖6　優(yōu)化后含油飽和度分布

圖7　優(yōu)化前流線

圖8　優(yōu)化后流線

2.2.2優(yōu)化時(shí)機(jī)

優(yōu)化時(shí)機(jī)對(duì)該方法的影響程度也非常明顯，進(jìn)行不同含水階段的優(yōu)化效果對(duì)比，分別模擬了含水率30%，50%，70%，80%，90%和95%時(shí)的方案。對(duì)比不同時(shí)期的優(yōu)化效果，進(jìn)行優(yōu)化時(shí)機(jī)的分析（見(jiàn)圖9）。

不同含水階段優(yōu)化結(jié)果對(duì)比表明，早期優(yōu)化的效果最優(yōu)，優(yōu)化時(shí)間越晚，效果越差；當(dāng)含水率90%以后再優(yōu)化效果會(huì)劇烈下降。因此，最佳的優(yōu)化時(shí)機(jī)是在含水率80%階段，即使是含水率在95%時(shí)進(jìn)行優(yōu)化，優(yōu)化結(jié)果對(duì)比無(wú)優(yōu)化的方案仍然采收率是提高的。由此可以得出在高含水期進(jìn)行流線優(yōu)化仍然能夠提高水驅(qū)效果。

圖9　不同時(shí)期優(yōu)化采收率效果

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

選擇中原油田某區(qū)塊進(jìn)行試驗(yàn)區(qū)的液流方向優(yōu)化研究。對(duì)目標(biāo)油藏進(jìn)行流線模擬評(píng)價(jià)開(kāi)發(fā)指標(biāo)，識(shí)別高效驅(qū)替方向和低效驅(qū)替方向；通過(guò)液流方向優(yōu)化方法優(yōu)化目標(biāo)油藏單井配注量配產(chǎn)量，進(jìn)行井組的挑選和液流方向優(yōu)化礦場(chǎng)應(yīng)用以及效果分析。

3.1開(kāi)發(fā)歷史注采優(yōu)化

運(yùn)用該方法，從1990年開(kāi)始進(jìn)行優(yōu)化（見(jiàn)圖10）。優(yōu)化前后整個(gè)區(qū)塊的注入量和產(chǎn)液量完全相同，不同的是每口注水井的注入量，同樣每口油井的產(chǎn)液量會(huì)有差別，但是為了便于優(yōu)化前后進(jìn)行對(duì)比，假設(shè)整個(gè)區(qū)塊的產(chǎn)液量不變。

圖10　油田注水利用率對(duì)比

3.2預(yù)測(cè)方案優(yōu)化

儲(chǔ)層含水級(jí)別統(tǒng)計(jì)顯示，儲(chǔ)層水淹及動(dòng)用程度較高，儲(chǔ)層含水率95%～98%是油層的主體，低含水率儲(chǔ)層分布較少，如含水率小于60%的儲(chǔ)層只有4%。不同含水率儲(chǔ)層注水分配見(jiàn)圖11，高含水率儲(chǔ)層注水分配程度較高，而低含水率儲(chǔ)層的水驅(qū)分配程度非常低。如含水率小于90%的儲(chǔ)層注入水還不到5%，對(duì)于高含水方向，尤其是含水率大于95%的儲(chǔ)層，注入水分配程度超過(guò)了70%，這些水必然形成無(wú)效的水循環(huán)，造成油井含水率高，開(kāi)發(fā)效果差。優(yōu)化后含水率大于98%的儲(chǔ)層注入水分配率由33.3%下降到20.9%，而含水率在60%～90%的儲(chǔ)層由4.6%上升到15.5%。擴(kuò)大了對(duì)低含水率儲(chǔ)層的水洗強(qiáng)度（見(jiàn)圖12、圖13）。

圖11　優(yōu)化前注入水分布直方圖

圖12　優(yōu)化后注入水分布直方圖

圖13　優(yōu)化前后采出程度

通過(guò)流線優(yōu)化對(duì)試驗(yàn)區(qū)塊進(jìn)行流線優(yōu)化模擬。模擬油藏從2013年到2023年的開(kāi)發(fā)動(dòng)態(tài)指標(biāo)，這期間注入水量427×104m3，優(yōu)化前累產(chǎn)水396×104m3，累產(chǎn)油22×104t；優(yōu)化前后的整個(gè)區(qū)塊的注水量相同，優(yōu)化后累產(chǎn)水394×104m3，累產(chǎn)油24×104t，增油2×104t。

4　結(jié)論

1）建立了基于改進(jìn)的CRM模型為驅(qū)動(dòng)、EnOpt集合優(yōu)化為最優(yōu)化手段、以經(jīng)濟(jì)效益為目標(biāo)函數(shù)的液流方向自動(dòng)化優(yōu)化方法，來(lái)求解油藏最優(yōu)配產(chǎn)配注量，得到每個(gè)油水井不同時(shí)期的最優(yōu)注采參數(shù)，可以改善油藏驅(qū)替效率，提高油藏最終采收率10%以上。

2）根據(jù)剩余油的分布規(guī)律，建立流線模擬評(píng)價(jià)模型，對(duì)目標(biāo)油藏進(jìn)行模擬開(kāi)發(fā)指標(biāo)研究，可以識(shí)別高效驅(qū)替方向和低效驅(qū)替方向，優(yōu)化目標(biāo)油藏單井配注量配產(chǎn)量，實(shí)現(xiàn)有效注水，提高開(kāi)發(fā)調(diào)整效果。

［1］李東玻.國(guó)內(nèi)外高含水油田特高含水期主要技術(shù)措施及啟示［J］.當(dāng)代石油石化，2013，26（10）：13-15.

［2］王元搏.特高含水期水驅(qū)剩余油挖潛技術(shù)研究［J］.中國(guó)石油和化工標(biāo)準(zhǔn)與質(zhì)量，2011，31（6）：160.

［3］王華，邴紹獻(xiàn)，張海燕，等.特高含水期水驅(qū)油效率計(jì)算新方法［J］.斷塊油氣田，2013，20（2）：201-203.

［4］王德龍，郭平，喻高明.雙河油田注采井組剩余油分布規(guī)律研究［J］.石油天然氣學(xué)報(bào)，2012，34（2）：128-130.

［5］鄧新穎，殷旭東，尹承棣，等.特高含水期剩余油分析方法［J］.斷塊油氣田，2004，11（5）：54-57.

［6］俞啟泰，謝緒權(quán)，羅洪.改變液流方向可行性評(píng)價(jià)公式［J］.斷塊油氣田，2000，7（3）：29-32.

［7］計(jì)秉玉，楊際平，呂志國(guó).改變液流方向的數(shù)值模擬計(jì)算［J］.大慶石油地質(zhì)與開(kāi)發(fā)，1994，13（2）：73-74.

［8］周書(shū)院，陳雷，解立春，等.改變液流方向技術(shù)在安塞油田王窯區(qū)中高含水期措施挖潛中的應(yīng)用［J］.石油天然氣學(xué)報(bào)，2007，29（3）：423-425.

［9］王健，謝華鋒.輪南油田中高含水期液流改向數(shù)值模擬研究［J］.石油地質(zhì)與工程，2011（3）：90-92.

［10］李偉才，崔連訓(xùn)，趙蕊.水動(dòng)力改變液流方向技術(shù)在低滲透油藏中的應(yīng)用：以新疆寶浪油田北區(qū)為例［J］.石油與天然氣地質(zhì)，2012，33（5）：796-801.

［11］NeliaJafroodi，DongxiaoZhang.Newmethodforreservoir characterization and optimization using CRM-EnOpt approach［J］. Journal of Petroleum Science and Engineering，2011，77：155-171.

［12］張婷.基于流線模擬的井組注采關(guān)系定量研究［J］.長(zhǎng)江大學(xué)學(xué)報(bào)：自然版，2014，11（20）：83-86.

［13］羅二輝，胡永樂(lè).流線數(shù)值模擬中的流線追蹤技術(shù)［J］.油氣井測(cè)試，2011，22（3）：10-13.

［14］王冰，王煥.流線模擬技術(shù)在油藏?cái)?shù)值模擬中的應(yīng)用［J］.內(nèi)蒙古石油化工，2012，（7）：80-82.

［15］王敏.泥頁(yè)巖地層孔隙流體實(shí)驗(yàn)數(shù)據(jù)分析［J］.石油實(shí)驗(yàn)地質(zhì)，2013，35（2）：199-201.

［16］崔傳智，張繼慶，楊勇，等.斷塊油藏典型井組特高含水期配產(chǎn)配注優(yōu)化研究［J］.科學(xué)技術(shù)與工程，2013，13（3）：697-700.

［17］寧士華，肖斐，束寧凱.特高含水開(kāi)發(fā)期曲流河儲(chǔ)層構(gòu)型深化研究及應(yīng)用［J］.斷塊油氣田，2013，20（3）：354-358.

（編輯楊會(huì)朋）

Fluid flow direction optimization and streamline simulation for high water reservoir

Li Caixue1，Shen Xi2，Jia Weiping3，Mao Lihua2，Zhu Liming2
（1.No.6 Oil Production Plant，Zhongyuan Oilfield Company，SINOPEC，Dongming 274515，China；2.Research Institute of Exploration and Development，Zhongyuan Oilfield Company，SINOPEC，Puyang 457001，China；3.No.1 Oil Production Plant，Zhongyuan Oilfield Company，SINOPEC，Puyang 457172，China）

Because of the planar heterogeneity and different water flooding degree，most oil fields in China have already entered into high water cut development stage，and the residual oil distributes sparsely.In order to optimize the water flooding development，the direction of fluid flow needs to be optimized to improve water drive recovery efficiency.Combining the numerical simulation and optimaltheory，based on improved CRM modeland EnOpt set optimization，taking economic benefit as the objective function，reservoir production and injection proration are quantified in a high water cut reservoir in Zhongyuan Oilfield.Development indexes are predicted by phased fine geological model and streamline simulator.The parameters about product well drainage area，water flooding sweep area of injection well，water flooding speed，inter-wellreservoir sweep area and inter-wellreservoir water are studied.

fluid flow direction optimization；streamline simulation；production and injection proration

TE312.

A

10.6056/dkyqt201505021

2015-03-01；改回日期：2015-07-22。

李才學(xué)，男，1968年生，高級(jí)工程師，碩士，1990年畢業(yè)于成都地質(zhì)學(xué)院石油地質(zhì)專業(yè)，現(xiàn)在從事開(kāi)發(fā)管理與研究工作。E-mail：lclcx001@163.com。

引用格式：李才學(xué)，沈曦，賈衛(wèi)平，等.高含水期油藏液流方向優(yōu)化及流線模擬［J］.斷塊油氣田，2015，22（5）：641-646.

Li Caixue，Shen Xi，Jia Weiping，et al.Fluid flow direction optimization and streamline simulation for high water reservoir［J］.Fault-Block Oil&Gas Field，2015，22（5）：641-646.