單井含水率，考慮到了壓力等因素，使得單井含水率的運(yùn)算結(jié)果更加準(zhǔn)確；
[0058] 4、在計(jì)算各個(gè)井點(diǎn)的壓力時(shí)，以連通單元為對(duì)象，基于物質(zhì)守恒定律，建立物質(zhì)守 恒方程，可以快速地計(jì)算出各個(gè)井點(diǎn)的壓力，提高了單井含水率的運(yùn)算速度，進(jìn)而可以提高 油藏井間連通關(guān)系的建立速度；
[0059] 5、在計(jì)算各個(gè)井點(diǎn)的含水率導(dǎo)數(shù)時(shí)，基于飽和度前緣推進(jìn)理論，整個(gè)飽和度追蹤 計(jì)算都是采用半解析方法，即僅采用某井點(diǎn)的上游井點(diǎn)進(jìn)行求解，可以快速、穩(wěn)定地計(jì)算出 各個(gè)井點(diǎn)的含水率導(dǎo)數(shù)，提高了單井含水率的運(yùn)算速度，進(jìn)而可以提高油藏井間連通關(guān)系 的建立速度。
[0060] 本發(fā)明的其它特征和優(yōu)點(diǎn)將在隨后的說明書中闡述，并且，部分地從說明書中變 得顯而易見，或者通過實(shí)施本發(fā)明而了解。本發(fā)明的目的和其他優(yōu)點(diǎn)可通過在說明書、權(quán)利 要求書以及附圖中所特別指出的結(jié)構(gòu)來實(shí)現(xiàn)和獲得。
【附圖說明】
[0061] 附圖用來提供對(duì)本發(fā)明的進(jìn)一步理解，并且構(gòu)成說明書的一部分，與本發(fā)明的實(shí) 施例共同用于解釋本發(fā)明，并不構(gòu)成對(duì)本發(fā)明的限制。在附圖中：
[0062] 圖1示出了本發(fā)明提供的油藏井間連通關(guān)系的建立方法的流程圖；
[0063] 圖2示出了本發(fā)明提供的油藏井間連通關(guān)系的簡(jiǎn)化示意圖；
[0064] 圖3示出了本發(fā)明的實(shí)施例2提供的單井含水率計(jì)算步驟的流程圖；
[0065] 圖4示出了含水率導(dǎo)數(shù)與含水飽和度以及含水飽和度與含水率的函數(shù)關(guān)系曲線。
【具體實(shí)施方式】
[0066] 以下將結(jié)合附圖及實(shí)施例來詳細(xì)說明本發(fā)明的實(shí)施方式，借此對(duì)本發(fā)明如何應(yīng)用 技術(shù)手段來解決技術(shù)問題，并達(dá)成技術(shù)效果的實(shí)現(xiàn)過程能充分理解并據(jù)以實(shí)施。需要說明 的是，只要不構(gòu)成沖突，本發(fā)明中的各個(gè)實(shí)施例以及各實(shí)施例中的各個(gè)特征可以相互結(jié)合， 所形成的技術(shù)方案均在本發(fā)明的保護(hù)范圍之內(nèi)。
[0067] 圖1示出了本發(fā)明的提供的油藏井間連通關(guān)系的建立方法的流程圖。該油藏井間 連通關(guān)系的建立方法包括以下步驟：
[0068] S100、連通單元?jiǎng)澐植襟E：將油藏井間連通區(qū)域劃分成多個(gè)連通單元，并使每?jī)蓚€(gè) 井點(diǎn)之間的連接區(qū)域?qū)?yīng)于一個(gè)連通單元，且各個(gè)連通單元的特征參數(shù)包括平均傳導(dǎo)率和 孔隙體積。
[0069] 圖2示出了本發(fā)明提供的油藏井間連通關(guān)系的簡(jiǎn)化示意圖。如圖2所示，將每一 油井視為一個(gè)井點(diǎn)，每?jī)蓚€(gè)井點(diǎn)之間的連接區(qū)域視為一個(gè)連通單元，且每一連通單元包括 平均傳導(dǎo)率和孔隙體積兩個(gè)特征參數(shù)。諸如對(duì)于井點(diǎn)i和井點(diǎn)j來說，兩者之間的連通單 元包括平均傳導(dǎo)率！和孔隙體積Vpu兩個(gè)特征參數(shù)。
[0070] 現(xiàn)有的油藏井間連通關(guān)系通常僅考慮了井點(diǎn)之間的連通關(guān)系，并未考慮井間連通 區(qū)域的實(shí)際地質(zhì)特征；此外，現(xiàn)有的油藏井間連通關(guān)系通常僅考慮了注采井之間的相互作 用，并未考慮生產(chǎn)井之間的相互作用，沒有考慮到停井和油井轉(zhuǎn)注等情況。
[0071] 本發(fā)明將井間連通區(qū)域劃分成多個(gè)連通單元，并使各個(gè)連通單元的特性參數(shù)包括 平均傳導(dǎo)率和孔隙體積。本發(fā)明采用平均傳導(dǎo)率來表征連通單元中油水的流動(dòng)特性，能夠 反映出各個(gè)井間(包括注采井和生產(chǎn)井）的相互作用，考慮到了停井和油井轉(zhuǎn)注等情況；采 用孔隙體積來表征連通單元的實(shí)際地質(zhì)特征，進(jìn)而考慮到了整個(gè)油藏井間的實(shí)際地質(zhì)特 征。
[0072] S200、單井含水率計(jì)算步驟，根據(jù)各個(gè)所述連通單元的所述平均傳導(dǎo)率和所述孔 隙體積，計(jì)算出各個(gè)所述井點(diǎn)的單井含水率。
[0073] S300、生產(chǎn)動(dòng)態(tài)指標(biāo)計(jì)算步驟，根據(jù)各個(gè)所述井點(diǎn)的所述單井含水率，計(jì)算出生產(chǎn) 動(dòng)態(tài)指標(biāo)。
[0074] S400、參數(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ù)修正。當(dāng)計(jì)算出的生產(chǎn)動(dòng)態(tài)指標(biāo)與其實(shí)際的生產(chǎn)動(dòng)態(tài)指標(biāo)不一致時(shí)，則需要對(duì)該連通 單元的平均傳導(dǎo)率和孔隙體積進(jìn)行參數(shù)修正。
[0075] 需要注意的是：在將計(jì)算出的各個(gè)井點(diǎn)的生產(chǎn)動(dòng)態(tài)指標(biāo)與其實(shí)際的生產(chǎn)動(dòng)態(tài)指標(biāo) 進(jìn)行擬合時(shí)，可以采用但不限于人工試算法、輔助優(yōu)化等方法。
[0076] 現(xiàn)有的油藏井間連通關(guān)系在模型參數(shù)修正過程中，通常采用諸如注采速度等特征 參數(shù)，計(jì)算出生產(chǎn)動(dòng)態(tài)數(shù)據(jù)，通過將計(jì)算出的生產(chǎn)動(dòng)態(tài)數(shù)據(jù)與實(shí)際的生產(chǎn)動(dòng)態(tài)數(shù)據(jù)進(jìn)行擬 合，對(duì)模型參數(shù)進(jìn)行修正。
[0077] 相對(duì)地，本發(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)行修正。
[0078] 由于單井含水率參數(shù)依賴于連通單元的平均傳導(dǎo)率和孔隙體積計(jì)算而來，平均傳 導(dǎo)率能夠反映出井間的相互作用，考慮到了停井和油井轉(zhuǎn)注等情況，而孔隙體積能夠反映 出油藏的實(shí)際地質(zhì)特征。由此可見，相比于注采速度等特征參數(shù)，單井含水率參數(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)確、可靠。
[0079] 圖3示出了本發(fā)明的提供的單井含水率計(jì)算步驟的流程圖。如圖3所示，上述單 井含水率計(jì)算步驟S200包括：
[0080] S210、以連通單元為對(duì)象，建立物質(zhì)守恒方程，計(jì)算出各個(gè)井點(diǎn)的壓力。
[0081] S220、根據(jù)各個(gè)連通單元的平均傳導(dǎo)率以及各個(gè)井點(diǎn)的壓力，計(jì)算出各個(gè)連通單 元的流量。
[0082] S230、根據(jù)各個(gè)連通單元的孔隙體積和流量，計(jì)算出各個(gè)連通單元的無因次累積 流量。
[0083] S240、根據(jù)各個(gè)連接單元的無因次累積流量，基于飽和度前緣推進(jìn)理論，計(jì)算出各 個(gè)井點(diǎn)的含水率導(dǎo)數(shù)。
[0084] S250、根據(jù)各個(gè)井點(diǎn)的含水率導(dǎo)數(shù)，計(jì)算出各個(gè)井點(diǎn)的含水率。
[0085] S260、根據(jù)各個(gè)井點(diǎn)的含水率以及各個(gè)連通單元的流量，計(jì)算出各個(gè)井點(diǎn)的單井 含水率。
[0086] 本發(fā)明在計(jì)算各個(gè)井點(diǎn)的單井含水率時(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)的單井含水率，考慮到了壓力等因素，使得單井含水率的運(yùn)算結(jié)果更 加準(zhǔn)確。
[0087] 在計(jì)算各個(gè)井點(diǎn)的壓力時(shí)，以連通單元為對(duì)象，基于物質(zhì)守恒定律，建立物質(zhì)守恒 方程，可以快速地計(jì)算出各個(gè)井點(diǎn)的壓力，提高了單井含水率的運(yùn)算速度，進(jìn)而可以提高油 藏井間連通關(guān)系的建立速度。
[0088] 在計(jì)算各個(gè)井點(diǎn)的含水率導(dǎo)數(shù)時(shí)，基于飽和度前緣推進(jìn)理論，可以快速、穩(wěn)定地計(jì) 算出各個(gè)井點(diǎn)的含水率導(dǎo)數(shù)，提高了單井含水率的運(yùn)算速度，進(jìn)而可以提高油藏井間連通 關(guān)系的建立速度。
[0089] 在本發(fā)明的一個(gè)實(shí)施例中，上述步驟S210包括以下步驟：
[0090] 步驟1、以連通單元為對(duì)象，建立物質(zhì)守恒方程：
[0091]
[0092] 其中，μ 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ù)，i和j為 正整數(shù)。
[0093] 步驟2、由于連通單元的平均傳導(dǎo)率戈
可將上述物質(zhì)守恒方程整理成：
[0094]
[0095] 其中，Tu表示井點(diǎn)i與井點(diǎn)j之間的連通單元的平均傳導(dǎo)率。
[0096] 步驟3、將整理后的所述物質(zhì)守恒方程進(jìn)行隱式差分運(yùn)算，得出壓力方程：
[0097]
[0098] 其中，pt+1表示在t+Ι時(shí)刻時(shí)井點(diǎn)i的壓力，Pt表示在t時(shí)刻時(shí)井點(diǎn)i的壓力，表 示在t+Ι時(shí)刻時(shí)井點(diǎn)j的壓力，
Λ t表示t到 t+Ι時(shí)刻的時(shí)間步長(zhǎng)。
[0099] 步驟4、將上述壓力方程整理成矩陣方程：
[0100]
，從而計(jì)算出各個(gè)井點(diǎn)的壓 力。
[0101] 從上述方程可以看出，壓力方程的個(gè)數(shù)與井點(diǎn)的個(gè)數(shù)相同，通過以連通單元為對(duì) 象，建立物質(zhì)守恒方程，可以快速地計(jì)算出各個(gè)井點(diǎn)的壓力，提高了單井含水率的運(yùn)算速 度，進(jìn)而可以提高油藏井間連通關(guān)系的建立速度。
[0102] 在本發(fā)明的另一實(shí)施例中，上述步驟S220包括：
[0103] 根據(jù)公式Q (i，j)=Tij (Pj-Pi)計(jì)算出各個(gè)所述連通單元的流量，
[0104] 其中，Q(i, j)表示井點(diǎn)i與井點(diǎn)j之間的連通單元的流量,且井點(diǎn)j為井點(diǎn)i的 上游井點(diǎn)，Tij表示井點(diǎn)i與井點(diǎn)j之間的連通單元的平均傳導(dǎo)率,Pi表示井點(diǎn)i的壓力，P j表示井點(diǎn)j的壓力。
[0105] 在本發(fā)明的另一實(shí)施例中，上述步驟S230包括：
[0106] 根據(jù)公式
t算出各個(gè)所述連通單元的所述無因次累積流量， 其中，Cv(i, j)表示井點(diǎn)j與井點(diǎn)i之間的