油藏井間連通關(guān)系的建立方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及油氣田開發(fā)領(lǐng)域，具體說，涉及一種油藏井間連通關(guān)系的建立方法。
【背景技術(shù)】
[0002] 縫洞型油藏具有儲(chǔ)集空間分布隨機(jī)性強(qiáng)、儲(chǔ)層非均質(zhì)性嚴(yán)重等特點(diǎn)，其流體流動(dòng) 特征及油水流動(dòng)規(guī)律相比常規(guī)砂巖油藏更為復(fù)雜，因此對(duì)該類油藏進(jìn)行連通性評(píng)價(jià)一直是 油田開發(fā)動(dòng)態(tài)分析和油藏描述工作的難點(diǎn)和重點(diǎn)。
[0003] 如何獲取該類型油藏井間連通關(guān)系，對(duì)于認(rèn)識(shí)儲(chǔ)層地質(zhì)特征和剩余油分布規(guī)律、 指導(dǎo)后期注采政策優(yōu)化和堵水調(diào)剖等措施調(diào)整、提高油藏經(jīng)營管理水平具有重要的作用。
[0004] 常用的井間連通關(guān)系研究方法歸納起來可以分為兩類：現(xiàn)場方法和基于生產(chǎn)動(dòng)態(tài) 的分析方法。現(xiàn)場方法包括示蹤劑測試、井間微地震、試井等，操作較復(fù)雜、解釋周期長且費(fèi) 用較高，因而適用范圍有限，無法大范圍推廣。
[0005] 基于生產(chǎn)動(dòng)態(tài)的分析方法主要是通過建立相關(guān)模型，利用生產(chǎn)動(dòng)態(tài)數(shù)據(jù)，進(jìn)行擬 合對(duì)比來修正模型參數(shù)。由于生產(chǎn)動(dòng)態(tài)數(shù)據(jù)獲取方便，該方法操作簡單、容易推廣。
[0006] 目前基于生產(chǎn)動(dòng)態(tài)的分析方法的模型主要包括多?；貧w模型、電容模型和系統(tǒng)分 析模型等。然而，現(xiàn)有的模型但普遍存在以下問題：（1)模型參數(shù)沒有明確的地質(zhì)意義，只 能反映井間相對(duì)的連通關(guān)系，無法反映井間真實(shí)地質(zhì)特征；（2)只考慮注采井之間的作用， 沒有考慮生產(chǎn)井之間的相互干擾，且模型本身無法考慮關(guān)停井及油井轉(zhuǎn)注等情況，因此反 演結(jié)果受油藏措施調(diào)整影響較大、可靠性差；（3)模型參數(shù)反演過程中主要使用注采速度 等，沒有結(jié)合含水率和壓力等動(dòng)態(tài)指標(biāo)、難以準(zhǔn)確反映實(shí)際油井生產(chǎn)狀況和儲(chǔ)層地質(zhì)特征。
【發(fā)明內(nèi)容】
[0007] 為了解決上述技術(shù)問題，本發(fā)明提供了一種油藏井間連通關(guān)系的建立方法，包 括：
[0008] 連通單元?jiǎng)澐植襟E，將所述油藏井間連通區(qū)域劃分為多個(gè)連通單元，并使每兩個(gè) 井點(diǎn)之間的連接區(qū)域?qū)?yīng)于一個(gè)連通單元，且各個(gè)所述連通單元的特性參數(shù)包括平均傳導(dǎo) 率和孔隙體積；
[0009] 單井含水率計(jì)算步驟，根據(jù)各個(gè)所述連通單元的所述平均傳導(dǎo)率和所述孔隙體 積，計(jì)算出各個(gè)所述井點(diǎn)的單井含水率；
[0010] 生產(chǎn)動(dòng)態(tài)指標(biāo)計(jì)算步驟，根據(jù)各個(gè)所述井點(diǎn)的所述單井含水率，計(jì)算出生產(chǎn)動(dòng)態(tài) 指標(biāo)；
[0011] 參數(shù)修正步驟，通過分別將計(jì)算出的各個(gè)所述井點(diǎn)的生產(chǎn)動(dòng)態(tài)指標(biāo)與其實(shí)際的生 產(chǎn)動(dòng)態(tài)指標(biāo)進(jìn)行擬合，對(duì)各個(gè)所述連通單元的所述平均傳導(dǎo)率和所述孔隙體積進(jìn)行參數(shù)修 正。
[0012] 根據(jù)本發(fā)明的一個(gè)實(shí)施例，所述單井含水率計(jì)算步驟包括：
[0013] 以所述連通單元為對(duì)象，建立物質(zhì)守恒方程，計(jì)算出各個(gè)井點(diǎn)的壓力；
[0014] 根據(jù)各個(gè)所述連通單元的所述平均傳導(dǎo)率以及各個(gè)所述井點(diǎn)的所述壓力，計(jì)算出 各個(gè)所述連通單元的流量；
[0015] 根據(jù)各個(gè)所述連通單元的所述孔隙體積和所述流量，計(jì)算出各個(gè)所述連通單元的 無因次累積流量；
[0016] 根據(jù)各個(gè)所述連接單元的所述無因次累積流量，基于飽和度前緣推進(jìn)理論，計(jì)算 出各個(gè)所述井點(diǎn)的含水率導(dǎo)數(shù)；
[0017] 根據(jù)各個(gè)所述井點(diǎn)的所述含水率導(dǎo)數(shù)，計(jì)算出各個(gè)所述井點(diǎn)的含水率；
[0018] 根據(jù)各個(gè)所述井點(diǎn)的所述含水率以及各個(gè)所述連通單元的所述流量，計(jì)算出各個(gè) 所述井點(diǎn)的單井含水率。
[0019] 根據(jù)本發(fā)明的一個(gè)實(shí)施例，所述以所述連通單元為對(duì)象，建立物質(zhì)守恒方程，計(jì)算 出各個(gè)井點(diǎn)的壓力，包括：
[0020] 以所述連通單元為對(duì)象，建立物質(zhì)守恒方程：
[0021]
[0022] 其中，μ i表示液相粘度，Iiij表示井點(diǎn)i與井點(diǎn)j之間的連通單元的平均滲透率， Aij表示井點(diǎn)i與井點(diǎn)j之間的連通單元的平均滲流截面積，Lij表示井點(diǎn)i與井點(diǎn)j之間 的距離，Pi表示井點(diǎn)i的壓力，Pj表示井點(diǎn)j的壓力，qi表示井點(diǎn)i的流量速度，Vpi表示井 點(diǎn)i的流量控制體積，Ct表示巖石綜合壓縮系數(shù)，t表示時(shí)刻，η表示井點(diǎn)的總數(shù)；
[0023] 由于所述連通單元的平均傳導(dǎo)率為
，將所述物質(zhì)守恒方程整理成：
[0024]
[0025] 其中，表示井點(diǎn)i與井點(diǎn)j之間的連通單元的平均傳導(dǎo)率；
[0026] 將整理后的所述物質(zhì)守恒方程進(jìn)行隱式差分運(yùn)算，得出壓力方程：
[0027]
[0028] 其中，p|+1表示在t+Ι時(shí)刻時(shí)井點(diǎn)i的壓力，P1表示在t時(shí)刻時(shí)井點(diǎn)i的壓力，P 1/1表 示在t+Ι時(shí)刻時(shí)井點(diǎn)j的壓力，
[0032] Λ t表示t到t+Ι時(shí)刻的時(shí)間步長；[0033] 將所述壓力方程整理成矩陣方程：
[0029]
[0030]
[0031]
[0034]
[0035] 從而計(jì)算出各個(gè)井點(diǎn)的壓力。
[0036] 根據(jù)本發(fā)明的一個(gè)實(shí)施例，所述根據(jù)各個(gè)所述連通單元的平均傳導(dǎo)率以及各個(gè)所 述井點(diǎn)的壓力，計(jì)算出各個(gè)所述連通單元的流量，包括：
[0037] 根據(jù)公式Q(i，W=Tij(Pj-Pi)計(jì)算出各個(gè)所述連通單元的流量，其中，Q(i，j)表示 井點(diǎn)i與井點(diǎn)j之間的連通單元的流量，Tij表示井點(diǎn)i與井點(diǎn)j之間的連通單元的平均傳 導(dǎo)率，Pi表示井點(diǎn)i的壓力，Pj表示井點(diǎn)j的壓力。
[0038] 根據(jù)本發(fā)明的一個(gè)實(shí)施例，所述根據(jù)各個(gè)所述連通單元的所述孔隙體積和所述流 量，計(jì)算出各個(gè)所述連通單元的無因次累積流量，包括：
[0039] 根據(jù)公式
計(jì)算出各個(gè)所述連通單元的所述無因次累積流量， 其中，Cv(i，j)表示井點(diǎn)j與井點(diǎn)i之間的連通單元的無因次累積流量，
[0040] Q(i, j)表示井點(diǎn)i與井點(diǎn)j之間的連通單元的流量，
[0041] Vpij表示井點(diǎn)i與井點(diǎn)j之間的連通單元的孔隙體積。
[0042] 根據(jù)本發(fā)明的一個(gè)實(shí)施例，所述根據(jù)各個(gè)所述連接單元的所述無因次累積流量， 基于飽和度前緣推進(jìn)理論，計(jì)算出各個(gè)所述井點(diǎn)的含水率導(dǎo)數(shù)，包括：
[0043] 根據(jù)公式
計(jì)算出各個(gè)所述井點(diǎn)的含水率導(dǎo)數(shù)，其 中，表示從井點(diǎn)j追蹤到井點(diǎn)i時(shí)所述井點(diǎn)i的含水率導(dǎo)數(shù)且井點(diǎn)j為井點(diǎn)i的 上游井點(diǎn)，供CCC/))表示井點(diǎn)j的含水率導(dǎo)數(shù)，Cv(i, j)表示井點(diǎn)j與井點(diǎn)i之間的連通單 元的無因次累積流量。
[0044] 根據(jù)本發(fā)明的一個(gè)實(shí)施例，所述根據(jù)各個(gè)所述井點(diǎn)的含水率導(dǎo)數(shù)，計(jì)算出各個(gè)所 述井點(diǎn)的含水率，包括：
[0045] 根據(jù)含水率導(dǎo)數(shù)與含水飽和度的函數(shù)關(guān)系，計(jì)算出各個(gè)所述井點(diǎn)的含水飽和度； 以及根據(jù)含水飽和度與含水率的函數(shù)關(guān)系，計(jì)算出各個(gè)所述井點(diǎn)的含水率。
[0046] 根據(jù)本發(fā)明的一個(gè)實(shí)施例，所述根據(jù)各個(gè)所述井點(diǎn)的含水率以及各個(gè)所述連通單 元的流量，計(jì)算出各個(gè)所述井點(diǎn)的單井含水率，包括：
[0047]根據(jù)公式.
?計(jì)算出各個(gè)所述井點(diǎn)的單井含水率，其中， fw(i)表示井點(diǎn)i的單井含水率，Q(i，j)表示井點(diǎn)i與井點(diǎn)j之間的連通單元的流量， fw(i，j)表示從井點(diǎn)j追蹤到井點(diǎn)i時(shí)所述井點(diǎn)i的追蹤含水率且井點(diǎn)j為井點(diǎn)i的上游 井點(diǎn)，n'表示井點(diǎn)i的上游井點(diǎn)的總數(shù)。
[0048] 根據(jù)本發(fā)明的一個(gè)實(shí)施例，所述生產(chǎn)動(dòng)態(tài)指標(biāo)計(jì)算步驟包括至少計(jì)算出以下一種 生產(chǎn)動(dòng)態(tài)指標(biāo)：
[0049] 根據(jù)公式q^qiX (l_fw(i))計(jì)算出單井日產(chǎn)油；
[0050] 根據(jù)公式qwi=qiXfw(i)計(jì)算出單井日產(chǎn)水；
[0051] 根據(jù)公式qf。= Σ qiX (l-fw(i))計(jì)算出區(qū)塊日產(chǎn)油；
[0052] 根據(jù)公式qfw= Σ qiXfw(i)計(jì)算出區(qū)塊日產(chǎn)水；
[0053] 其中，Qi表示井口 i的流量速度，fw(i)表示井點(diǎn)i的單井含水率。
[0054] 本發(fā)明的有益效果是：
[0055] 1、本發(fā)明將井間連通區(qū)域劃分成多個(gè)連通單元，并使各個(gè)連通單元的特性參數(shù)包 括平均傳導(dǎo)率和孔隙體積，采用平均傳導(dǎo)率來表征連通單元中油水的流動(dòng)特性，能夠反映 出各個(gè)井間(包括注采井和生產(chǎn)井）的相互作用，考慮到了停井和油井轉(zhuǎn)注等情況；采用孔 隙體積來表征連通單元的實(shí)際地質(zhì)特征，進(jìn)而考慮到了整個(gè)油藏井間的實(shí)際地質(zhì)特征；
[0056] 2、在修正模型參數(shù)時(shí)，本發(fā)明采用各個(gè)井點(diǎn)的單井含水率，計(jì)算出各個(gè)井點(diǎn)的生 產(chǎn)動(dòng)態(tài)指標(biāo)，通過分別將計(jì)算出的各個(gè)井點(diǎn)的生產(chǎn)動(dòng)態(tài)指標(biāo)與其實(shí)際生產(chǎn)動(dòng)態(tài)指標(biāo)進(jìn)行擬 合，對(duì)模型參數(shù)（即平均傳導(dǎo)率和孔隙體積）進(jìn)行修正，由于單井含水率參數(shù)依賴于連通單 元的平均傳導(dǎo)率和孔隙體積計(jì)算而來，單井含水率參數(shù)能夠更加準(zhǔn)確地反映出油藏的實(shí)際 生產(chǎn)狀態(tài)和地質(zhì)特征，采用含水率計(jì)算出的生產(chǎn)動(dòng)態(tài)數(shù)據(jù)也更加準(zhǔn)確，從而使得根據(jù)擬合 結(jié)果，修正獲得的模型參數(shù)更加準(zhǔn)確、可靠；
[0057] 3、在計(jì)算單井含水率時(shí)，首先，計(jì)算出各個(gè)井點(diǎn)的壓力；其次，根據(jù)各個(gè)連通單元 的平均傳導(dǎo)率以及各個(gè)井點(diǎn)的壓力，計(jì)算出各個(gè)連通單元的流量；再者，根據(jù)各個(gè)連通單元 的孔隙體積和流量，計(jì)算出各個(gè)連通單元的無因次累積流量；最后，根據(jù)各個(gè)連通單元的無 因次累積流量，計(jì)算出各個(gè)井點(diǎn)的含水率導(dǎo)數(shù)，并根據(jù)各個(gè)井點(diǎn)的含水率導(dǎo)數(shù)，計(jì)算出各個(gè) 井點(diǎn)的