控制所述組分?jǐn)?shù)值模擬子單 元基于所構(gòu)建的組分?jǐn)?shù)值模擬模型，利用所述先驗(yàn)控制參數(shù)向量集合進(jìn)行組分?jǐn)?shù)值模擬， 以得到當(dāng)前時(shí)間步的所述多個(gè)相對(duì)滲透率先驗(yàn)?zāi)Ｐ退鶎?duì)應(yīng)的先驗(yàn)動(dòng)態(tài)數(shù)據(jù)預(yù)測(cè)值向量，控 制所述第一計(jì)算子單元計(jì)算每個(gè)所述先驗(yàn)動(dòng)態(tài)數(shù)據(jù)預(yù)測(cè)值向量與所述動(dòng)態(tài)數(shù)據(jù)觀測(cè)值向 量之間的差值，以得到校正控制參數(shù)向量集合；以及控制所述第二計(jì)算子單元計(jì)算所述校 正控制參數(shù)向量集合所對(duì)應(yīng)的先驗(yàn)動(dòng)態(tài)數(shù)據(jù)預(yù)測(cè)值向量的平均值，以得到當(dāng)前動(dòng)態(tài)數(shù)據(jù)預(yù) 測(cè)值向量等操作；
[0046] 第二控制子單元，用于控制所述第三計(jì)算子單元以及所述判斷子單元的操作，直 到得到滿(mǎn)足迭代收斂條件的目標(biāo)函數(shù)。
[0047] 在至少一實(shí)施例中，所述迭代計(jì)算單元中所構(gòu)建的目標(biāo)函數(shù)的表達(dá)式如下：
[0048]
[0049] 其中，0(m)為目標(biāo)函數(shù);m為所述表征模型中的控制參數(shù)向量;T為表征向量或矩陣 轉(zhuǎn)置的符號(hào);cUs為動(dòng)態(tài)數(shù)據(jù)觀測(cè)值向量;g(m)為動(dòng)態(tài)數(shù)據(jù)預(yù)測(cè)值向量;CD為動(dòng)態(tài)數(shù)據(jù)測(cè)量 誤差的協(xié)方差矩陣;mpricir為所述控制參數(shù)向量m的先驗(yàn)?zāi)Ｐ托畔⒌钠骄?Cm為先驗(yàn)?zāi)Ｐ偷?協(xié)方差矩陣。
[0050] 在至少一實(shí)施例中，所述第三獲取單元包括：
[0051] 確定子單元，用于根據(jù)所述符合迭代收斂條件的目標(biāo)函數(shù)，確定出所述控制參數(shù) 向量中各控制參數(shù)的當(dāng)前數(shù)值；
[0052]第四計(jì)算子單元，用于根據(jù)所述控制參數(shù)向量中各控制參數(shù)的當(dāng)前數(shù)值，計(jì)算油-水兩相系統(tǒng)的油相相對(duì)滲透率、水相相對(duì)滲透率以及油-氣兩相系統(tǒng)的氣相相對(duì)滲透率以 及油相相對(duì)滲透率；
[0053] 第五計(jì)算子單元，用于利用所計(jì)算的油-水兩相系統(tǒng)的油相相對(duì)滲透率、水相相對(duì) 滲透率以及油-氣兩相系統(tǒng)的氣相相對(duì)滲透率以及油相相對(duì)滲透率，計(jì)算三相滲流條件下 的油相相對(duì)滲透率。
[0054] 由以上本申請(qǐng)實(shí)施例提供的技術(shù)方案可見(jiàn)，本申請(qǐng)實(shí)施例通過(guò)基于目標(biāo)區(qū)域的 氣-水交替非混相驅(qū)替實(shí)驗(yàn)，獲取不同時(shí)刻的驅(qū)替壓差、含水率以及氣油比，所述不同時(shí)刻 的驅(qū)替壓差、含水率以及氣油比構(gòu)成了動(dòng)態(tài)數(shù)據(jù)觀測(cè)值向量;基于所構(gòu)建的表征模型所表 征的先驗(yàn)三相相對(duì)滲透率曲線以及所構(gòu)建的組分?jǐn)?shù)值模擬模型，獲取當(dāng)前動(dòng)態(tài)數(shù)據(jù)預(yù)測(cè)值 向量;利用所述動(dòng)態(tài)數(shù)據(jù)觀測(cè)值向量和所述當(dāng)前動(dòng)態(tài)數(shù)據(jù)預(yù)測(cè)值向量，對(duì)所構(gòu)建的目標(biāo)函 數(shù)進(jìn)行迭代計(jì)算，直到所述目標(biāo)函數(shù)的迭代計(jì)算滿(mǎn)足預(yù)設(shè)的迭代收斂條件;根據(jù)滿(mǎn)足所述 迭代收斂條件的目標(biāo)函數(shù)，獲取目標(biāo)區(qū)域的油-氣-水三相相對(duì)滲透率曲線，從而實(shí)現(xiàn)了獲 取適用于氣-水交替非混相驅(qū)油技術(shù)的油-氣-水三相相對(duì)滲透率曲線的目的。
【附圖說(shuō)明】
[0055] 為了更清楚地說(shuō)明本申請(qǐng)實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn) 有技術(shù)描述中所需要使用的附圖作簡(jiǎn)單地介紹，顯而易見(jiàn)地，下面描述中的附圖僅僅是本 申請(qǐng)中記載的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來(lái)講，在不付出創(chuàng)造性勞動(dòng)性的前提 下，還可以根據(jù)這些附圖獲得其他的附圖。
[0056] 圖1是本申請(qǐng)實(shí)施例提供的獲取三相相對(duì)滲透率曲線的方法的流程圖；
[0057]圖2是步驟S2所包括的子步驟流程圖；
[0058]圖3是驅(qū)替壓差及含水率的擬合效果圖；
[0059] 圖4是累積生產(chǎn)氣油比的擬合效果圖；
[0060] 圖5是步驟S3所包括的子步驟流程圖；
[0061] 圖6為油-水兩相系統(tǒng)的油相相對(duì)滲透率曲線以及水相相對(duì)滲透率曲線；
[0062] 圖7為油-氣兩相系統(tǒng)的油相相對(duì)滲透率曲線以及氣相相對(duì)滲透率曲線；
[0063]圖8為三相滲流條件下的油相相對(duì)滲透率曲線。
[0064]圖9是一種獲取三相相對(duì)滲透率曲線的裝置的模塊示意圖。
【具體實(shí)施方式】
[0065]本申請(qǐng)實(shí)施例提供一種獲取三相相對(duì)滲透率曲線的方法及裝置。
[0066]為了使本技術(shù)領(lǐng)域的人員更好地理解本申請(qǐng)中的技術(shù)方案，下面將結(jié)合本申請(qǐng)實(shí) 施例中的附圖，對(duì)本申請(qǐng)實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施 例僅僅是本申請(qǐng)一部分實(shí)施例，而不是全部的實(shí)施例。基于本申請(qǐng)中的實(shí)施例，本領(lǐng)域普通 技術(shù)人員在沒(méi)有作出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都應(yīng)當(dāng)屬于本申請(qǐng)保護(hù) 的范圍。
[0067] 下面結(jié)合附圖對(duì)本申請(qǐng)所述的獲取三相相對(duì)滲透率曲線的方法進(jìn)行詳細(xì)的說(shuō)明。 雖然本申請(qǐng)?zhí)峁┝巳缦率鰧?shí)施例或流程圖所述的方法操作步驟，但基于常規(guī)或者無(wú)需創(chuàng)造 性的勞動(dòng)在所述方法中可以包括更多或者更少的操作步驟。在邏輯性上不存在必要因果關(guān) 系的步驟中，這些步驟的執(zhí)行順序不限于本申請(qǐng)實(shí)施例提供的執(zhí)行順序。所述的方法的在 實(shí)際中的裝置或終端產(chǎn)品執(zhí)行時(shí)，可以按照實(shí)施例或者附圖所示的方法順序執(zhí)行或者并行 執(zhí)行。
[0068] 如圖1所示，本申請(qǐng)實(shí)施例提供了一種獲取三相相對(duì)滲透率曲線的方法，該方法包 括：
[0069] S1:基于目標(biāo)區(qū)域的氣-水交替非混相驅(qū)替實(shí)驗(yàn)，獲取不同時(shí)刻的驅(qū)替壓差、含水 率以及氣油比，所述不同時(shí)刻的驅(qū)替壓差、含水率以及氣油比構(gòu)成了動(dòng)態(tài)數(shù)據(jù)觀測(cè)值向量。
[0070] 可以采用真實(shí)的油藏巖樣，通過(guò)室內(nèi)實(shí)驗(yàn)確定巖石、流體物性參數(shù)以及兩相滲流 時(shí)的飽和度端點(diǎn)數(shù)據(jù)，包括:巖石的孔隙度、滲透率、壓縮系數(shù);油-水兩相初始含油飽和度、 油-水兩相殘余油飽和度、氣-液兩相臨界含氣飽和度、氣-液兩相殘余油飽和度等。在確定 出上述數(shù)據(jù)后，可以模擬油藏條件，針對(duì)目標(biāo)區(qū)域開(kāi)展氣-水交替非混相巖心驅(qū)替實(shí)驗(yàn)，例 如注烴氣或二氧化碳驅(qū)替，并實(shí)時(shí)記錄不同注入孔隙體積倍數(shù)(PV)下的驅(qū)替壓差、驅(qū)油效 率、含水率和氣油比等實(shí)驗(yàn)數(shù)據(jù);然后從所記錄的實(shí)驗(yàn)數(shù)據(jù)中獲取不同時(shí)刻的驅(qū)替壓差、含 水率以及氣油比(采出的原油量與天然氣的體積量之間的比值)等數(shù)據(jù)，并將所獲取的不同 時(shí)刻的驅(qū)替壓差、含水率以及氣油比構(gòu)成動(dòng)態(tài)數(shù)據(jù)觀測(cè)值向量。
[0071] 所述目標(biāo)區(qū)域可以是指整個(gè)勘探區(qū)域，也可以是指勘探區(qū)域中的部分區(qū)域。
[0072] S2:基于所構(gòu)建的表征模型所表征的先驗(yàn)三相相對(duì)滲透率曲線以及所構(gòu)建的組分 數(shù)值模擬模型，獲取當(dāng)前動(dòng)態(tài)數(shù)據(jù)預(yù)測(cè)值向量。
[0073] 所述表征模型可以用于表征油-水兩相相對(duì)滲透率曲線、油-氣兩相相對(duì)滲透率曲 線和/或油-氣-水三相相對(duì)滲透率曲線的形狀，其可以是三次B樣條表征模型，但并不限于 該模型。其構(gòu)建方法可以參考現(xiàn)有技術(shù)中的相關(guān)方法，在此不再贅敘。
[0074] 所述組分?jǐn)?shù)值模擬模型可以是基于PVT相態(tài)特征分析以及最小混相壓力測(cè)試所確 定的流體物性參數(shù)分布特征來(lái)構(gòu)建的。具體的，
[0075] 可以利用高壓物性?xún)x和高壓落球式粘度計(jì)的聯(lián)合裝置，開(kāi)展單次閃蒸、恒質(zhì)膨脹、 微分釋放和注氣膨脹實(shí)驗(yàn)，測(cè)定注氣摩爾分?jǐn)?shù)取不同值時(shí)，油藏溫度條件下的飽和壓力、流 體壓力和體積之間的關(guān)系、原油密度的變化規(guī)律以及粘度的變化規(guī)律。
[0076] 可以基于長(zhǎng)細(xì)管驅(qū)替實(shí)驗(yàn)，在油藏溫度和高于泡點(diǎn)壓力條件下，測(cè)定注入氣體的 最小混相壓力。通常情況下，可以將氣體突破時(shí)采收程度達(dá)到80%或最終采收率達(dá)到90% ~95%這兩個(gè)采收率水平作為判定驅(qū)替是否混相的標(biāo)準(zhǔn)。
[0077] 進(jìn)行PVT相態(tài)特征分析以及最小混相壓力測(cè)試可以有利于提高表征氣-水交替非 混相驅(qū)油技術(shù)如何提高采收率機(jī)理的準(zhǔn)確性，也有利于確保所構(gòu)建組分?jǐn)?shù)值模擬模型與實(shí) 際實(shí)驗(yàn)裝置的一致性。
[0078] 在獲得上述流體物性參數(shù)分布特征后，可以基于飽和壓力、原油密度的變化規(guī)律、 粘度的變化規(guī)律以及所測(cè)定的最小混相壓力來(lái)構(gòu)建組分?jǐn)?shù)值模擬模型。所述組分?jǐn)?shù)值模擬 模型可以是用于對(duì)流體進(jìn)行組分?jǐn)?shù)值模擬。所述組分?jǐn)?shù)值模擬是三相流體數(shù)值模擬的一種 形式，其可以涉及氣體相態(tài)變化、油相密度、粘度以及體積系數(shù)變化等，其主要通過(guò)降低氣 液相間性質(zhì)的差異、改善水油流度比以及提高驅(qū)油效率等來(lái)達(dá)到提高原油采收率的效果。
[0079] 在構(gòu)建所述組分?jǐn)?shù)值模擬模型后，可以利用所述表征模型所表征的先驗(yàn)三相相對(duì) 滲透率曲線