沈文潔，趙輝，劉偉，許凌飛，廖茂林

（長江大學石油工程學院，湖北 武漢 430100）

0 引言

油藏優(yōu)勢竄流通道識別對油層注水的開發(fā)效果起直接作用。根據(jù)以往地質(zhì)學家采樣實驗結(jié)果分析，縫洞型碳酸鹽巖油藏的儲存滲流空間主要為規(guī)模不同的溶洞、縫洞和裂縫等，構(gòu)造一般為多縫洞分布，縫洞大小規(guī)則不一，迂曲度較高［1］，因此，該類油藏井間連通方式較為復雜［2-3］。目前，塔河縫洞型油藏仍以注水開發(fā)為主要驅(qū)油方式，沿裂縫水淹水竄是注水開發(fā)中存在的主要問題。傳統(tǒng)的竄流通道識別方法［4-6］有含水指數(shù)特征曲線方法、霍爾曲線方法、灰色關(guān)聯(lián)分析方法、油藏動態(tài)關(guān)聯(lián)法以及模糊綜合評判法等。這些方法對于早期油藏開發(fā)有一定的借鑒指導作用，但后期油田為控水穩(wěn)油，實施了酸化壓裂、摻稀降黏、調(diào)剖堵水等措施，破壞了地層的滲流走向，改變了地層導流能力，使得計算結(jié)果產(chǎn)生了不同程度的誤差，降低了油藏竄流通道識別的準確性。本文針對該問題提出了基于可模擬油水動態(tài)的連通性模型［7-11］，這是一種不依賴于復雜地質(zhì)建模的模型。其主要思想是：利用油水井的日常注采動態(tài)數(shù)據(jù)，基于簡化表征的油藏注采系統(tǒng)，通過動態(tài)參數(shù)反演出傳導率和連通體積；結(jié)合示蹤劑資料更新不同階段地層流動特征參數(shù)，得到實時反映油藏模型的連通關(guān)系場圖；同時考慮井間干擾和注采受效時滯性等因素，對油藏區(qū)塊進行井史和注采關(guān)系分析，通過對比結(jié)果，最終實時反映油藏井間竄流通道及注水驅(qū)油效率。

1 連通性模型

井間連通性概念早在1997年就已提出，其計算模型主要包括多元回歸模型［8］、電容模型［9］和系統(tǒng)分析模型［10］。這些模型計算簡單，但考慮因素較少，特別是無法計算轉(zhuǎn)注后的產(chǎn)液參數(shù)。針對該問題，本文提出了一種基于趙輝等［10］研究的連通性模型。該模型主要以傳導率和連通體積來表征單元流動能力和物質(zhì)基礎，具體求解計算參見文獻［7，11］，本文僅詳述含水率的求解過程。

1.1 飽和度追蹤-含水率求解

基于建立的連通性模型，可以得到連通單元物質(zhì)平衡方程：

式中：Tij為第 i井和第 j井間的傳導率，m3/（s·MPa）；m為第i井的所有上游井點總數(shù)；pi，pj分別為第i井和第j井泄油區(qū)內(nèi)的平均壓力，MPa；α為單位換算系數(shù)（取值 9.8×10-6）；ρl為流體密度，kg/m3；g為重力加速度（取值9.8 m/s2）；Dij為第 i井和第 j井的中部深度差，m；dt為生產(chǎn)時間變化量，d；qi為第 i井流量（注入為正，產(chǎn)出為負），m3/s；Ct為綜合壓縮系數(shù)，MPa-1；Vi為第i井網(wǎng)格孔隙體積，m3。

對式（1）隱式差分離散，可得到壓力求解方程［7］。求出各井點壓力后，井點間連通單元內(nèi)流量為

式中：qij為從第j井流向第i井的流量，m3/s。

考慮油藏存在邊底水，結(jié)合前緣推進理論［11］進行推導。取x的某一上游位置xu（即xu

式中：f′w（Swu）為 xu處含水率導數(shù)；QDn為 n 時刻流入 xu到x形成單元的無因次累積注入量；f′w（Sw）為油層x位置的含水率導數(shù)。

進一步整理得到n時刻第j井追蹤至下游的第i井含水率導數(shù) f′w（Swnij）：

式中：f′w（Swnj），f′w（Swzi）分別為 n，z時刻第 j，i井處的含水率導數(shù)；QDnij，QDnji分別為n時刻第j井流向第i井的無因次累積流量和第i井流向第j井的無因次累積流量；QDzji為z時刻第i井流向第j井的無因次累積流量。

式（4）的修正保證了計算穩(wěn)定性［10］。通過式（4）插值求出第i井來自第j井方向的含水率，再依次求出各上游方向上的含水率，最終求出第i井的綜合含水率：

式中：fw（Swni）為第 i井的綜合含水率。

將連通性模型計算出的示蹤劑時刻劈分結(jié)果與現(xiàn)場給的劈分結(jié)果作對比。前者通過歷史擬合模型參數(shù)反演后計算得到，其值隨工作制度的改變而改變，反映了井間平面動態(tài)響應。假設第i井為水井，則其與周圍油井j間的流量分配系數(shù)αnij為

利用上述信息可進一步精確計算水井注水驅(qū)油效率，即注水井向周邊油井供水驅(qū)替出的原油總量與該井在該方向劈分過去的注水量比值：

式中：ηi為第i井注水效率；fwni為n時刻第i井含水率。

1.2 特征參數(shù)反演

基于隨機擾動和投影梯度類算法［10-11］，通過不斷優(yōu)化和調(diào)整油藏模型特征參數(shù)——傳導率和連通體積，使計算出的單井含水率、日產(chǎn)油量等結(jié)果與單井日常動態(tài)參數(shù)相吻合。換句話說，就是如何求解特征參數(shù)矩陣b，以使目標函數(shù)O（b） 取得最小值。

式中：b為油藏參數(shù)矩陣（這里主要指傳導率和連通體積）；br為先驗油藏模型估計；為模型參數(shù)的協(xié)方差逆矩陣；kobs為實際觀測數(shù)據(jù)；為先驗模型參數(shù)的協(xié)方差矩陣；h（b）為由數(shù)值模擬器計算得到的油藏觀測數(shù)據(jù)初值。

求解此類優(yōu)化問題，就是在滿足約束條件的同時，盡可能使目標函數(shù)值最小，最終得到對應的b。

為此，對約束條件采用了梯度投影方法進行迭代求解。求解公式為

式中：bl+1為第l+1步的迭代控制變量；ζ為搜索步長；T為Nu維投影矩陣；I為單位矩陣；B為由油藏參數(shù)b形成的矩陣；▽l（bl）為O（b）的隨機擾動梯度。

O（b）的隨機擾動梯度計算公式［10］為

式中：εl為擾動步長；Δl為 Nu維隨機擾動向量（其中所包含元素 Δlξ（ξ=1，2，…，Nu）為服從多元高斯分布的擾動向量）。

2 實例應用

塔河縫洞型油藏Z區(qū)塊屬于碳酸鹽巖油藏，該區(qū)塊2000年投產(chǎn)，油水井共計30口，其中關(guān)停井9口。截至2017年2月，累計產(chǎn)油量301.210×104t，累計注水量80.518×104t。目前采用傳統(tǒng)的注水驅(qū)油模式。該油藏中后期注水愈加頻繁，導致注入水無效循環(huán)嚴重，沿主流通道竄流，油井含水率快速上升，產(chǎn)量遞減快，整體采收率低。目前，該區(qū)塊5口井處于水淹狀態(tài)，7口井供液不足，9口井關(guān)停。在不改變主要開發(fā)方式和不影響高產(chǎn)穩(wěn)產(chǎn)井當前狀況的生產(chǎn)條件下，對區(qū)塊優(yōu)勢竄流通道進行判別，是該區(qū)塊亟待解決的問題。

應用井間連通性模型［11］，結(jié)合實際地質(zhì)參數(shù)對該油藏進行生產(chǎn)動態(tài)自動歷史擬合，反演結(jié)果如圖1所示。

為驗證該模型的可靠性，以D63井組和D12井組為例（見圖2）。

圖1 反演結(jié)果

圖2 井組井位示意

分析示蹤劑時刻注水劈分和本模型劈分結(jié)果，并以注水受效結(jié)果補充驗證。D63井示蹤劑測試時間為2007年7月。選取現(xiàn)場數(shù)據(jù)，基于連通性模型，計算示蹤劑時刻劈分系數(shù)（見圖3）?？梢钥闯?，B26CX，B35，B36H井示蹤劑時刻注水劈分和連通性模型注水劈分級別相當，且未劈分比例較多，存在外溢情況。對D12井組，采取了同樣的分析措施（見圖4）。

圖3 D63井組注水劈分結(jié)果對比

圖4 D12井組注水劈分結(jié)果對比

注水井 D12 井，周圍油井有 B15，B44，B47，B07井。D12井示蹤劑測試時間為2011年8月。由圖4可以看出，B47，B07，B15，B44 井數(shù)值相當，B11 井數(shù)值存在差異。原因是，在模型建立期間，只對Z區(qū)塊各井建模，而B11井不屬于該單元內(nèi)井，故該方向劈分系數(shù)無法給出。同時，結(jié)合注采井之間的動態(tài)響應，認為D12井的注水對B47，B07井的影響較大，劈分較多。

從2個井組中分別選取2口油井進行注采響應分析，并對歷史措施進行總結(jié)歸納（見表1）。

表1 各井組注水受效情況分析

上述井組示蹤劑分析結(jié)果及模型計算基本一致，說明該模型準確程度較高；同時，結(jié)合井史，對井組間注水受效進行分析，再一次驗證該模型的合理性?；谇捌趧討B(tài)響應分析，結(jié)合示蹤劑結(jié)果，利用連通性計算模型，最終得到表2中各單井數(shù)據(jù)，并據(jù)此給出潛力分析。

由表2分析可知，D63井注水，主要流向為B35，B26CX，B80，B11井，其次是 B36H，B14 井，最后是B25井。生產(chǎn)潛力分析認為：1）B14，B35井剩余油較多。B14井水淹與D63井注水沒有直接關(guān)系，目前處于低產(chǎn)、低含水率生產(chǎn)中，驅(qū)油效率高，建議提液；B35井間開生產(chǎn)，含水率下降，建議適當提液。2）B11，B36H井剩余油較多。D63井強注，使得2口井水淹，且B11井關(guān)井壓錐無效，屬于優(yōu)勢竄流通道，建議維持當前狀態(tài)生產(chǎn)或靶向調(diào)堵。3）B26CX，B25，B80井剩余油較少。D63井注水，使得B80井含水率異常，驅(qū)油效率較低，說明該方向存在優(yōu)勢竄流通道，建議維持當前生產(chǎn)；B26CX井未水淹，建議提液；B25井自身提液水竄，提液存在風險，建議維持當前生產(chǎn)。

D12井組主要流向B47，B07，B48井等，其次是B44，B15 井。生產(chǎn)潛力分析認為：1）B44，B47，B48 井剩余油較多。D12井注水，使得B44井水淹，該方向?qū)儆诟Z流通道，建議降液或維持當前生產(chǎn)或調(diào)剖堵水；B47井交替提液降液，含水率上升或下降，建議維持當前生產(chǎn)；B12井注水時，B48井交替提液降液，含水率波動較大，受水體和注入水共同作用，建議維持當前生產(chǎn)。2）B07井剩余油較高。D12井注水，未水淹，但驅(qū)油效率較低，建議維持當前狀態(tài)生產(chǎn)。3）B15井剩余油較少。B15井水淹與D12井注水沒有直接關(guān)系，可能是底水作用，目前轉(zhuǎn)層生產(chǎn)，驅(qū)油效率低，建議維持當前狀態(tài)生產(chǎn)。

表2 井組注水評價

3 結(jié)論

1）可模擬油水動態(tài)下的井間連通性模型，不及傳統(tǒng)地質(zhì)建模復雜，且在反演傳導率和連通體積等問題時，需要參數(shù)少，計算速度快，一定程度上可以反映某個方向上的滲流能力，為識別優(yōu)勢通道提供參考，有利于指導現(xiàn)場應用。

2）基于井間連通性模型，結(jié)合Z油藏開發(fā)動態(tài)特征，反演后的傳導率實時變化，與注采響應、示蹤劑結(jié)果相吻合，驗證了模型的準確性；在此基礎上計算出的注水井向周圍油井的注水劈分，反映了井間的注水流向，即優(yōu)勢通道，為下一步注水策略調(diào)整提供了依據(jù)。

3）采用該方法計算出的注水效率，能反映出該滲流通道的類型。注水效率高，說明驅(qū)油效果好，可以保持注水速度；反之，該方向易發(fā)生水竄，建議停注或減少注水量。