數(shù)據(jù)預(yù)測值向量。
[0113] S34:在得到所述校正控制參數(shù)向量集合所對應(yīng)的當(dāng)前動態(tài)數(shù)據(jù)預(yù)測值向量后，重 復(fù)步驟S31-S32,直到所述迭代計(jì)算滿足預(yù)設(shè)的迭代收斂條件，得到滿足所述迭代收斂條件 的目標(biāo)函數(shù)。
[0114] S4:根據(jù)滿足所述迭代收斂條件的目標(biāo)函數(shù)，獲取目標(biāo)區(qū)域的油-氣-水三相相對 滲透率曲線。
[0115] 在判斷出當(dāng)前迭代滿足預(yù)設(shè)的迭代收斂條件時，根據(jù)當(dāng)前迭代計(jì)算出的目標(biāo)函數(shù) 的當(dāng)前數(shù)值，獲取對應(yīng)于目標(biāo)區(qū)域的油-氣-水三相相對滲透率曲線。具體的，
[0116] 可以根據(jù)所述目標(biāo)函數(shù)的當(dāng)前數(shù)值，確定出所述控制參數(shù)向量中各控制參數(shù)的當(dāng) 前數(shù)值。然后，可以根據(jù)所述控制參數(shù)向量中各控制參數(shù)的當(dāng)前數(shù)值，計(jì)算油-水兩相系統(tǒng) 的油相相對滲透率、水相相對滲透率以及油-氣兩相系統(tǒng)的氣相相對滲透率以及油相相對 滲透率，計(jì)算公式可以表示如下：
[0117] (5)
[0?18]其中，krp為p相流體的相對滲透率;p = W，OW，g，Og ; Bj, 3 ( SpD )為四階(3次）的B樣條 基函數(shù)
[0119] 最后，可以利用油-水兩相系統(tǒng)的油相相對滲透率、水相相對滲透率以及油-氣兩 相系統(tǒng)的氣相相對滲透率以及油相相對滲透率，計(jì)算三相滲流條件下的油相相對滲透率， 其計(jì)算公式可以表示如下：
[0120]
(6)
[0121] 上式中，Kr。為三相滲流條件下的油相相對滲透率;Kr?為油-水兩相系統(tǒng)的油相相 對滲透率;K rw為油-水兩相系統(tǒng)的水相相對滲透率;Krcig為油-氣兩相系統(tǒng)的油相相對滲透 率;K rg為油-氣兩相系統(tǒng)的氣相相對滲透率。
[0122] 根據(jù)所計(jì)算出的三相滲流條件下的油相相對滲透率、水相相對滲透率以及氣相相 對滲透率，即可以得到油-氣-水三相相對滲透率曲線。
[0123] 圖6-圖8分別示出了所獲取的冀東油田中某一區(qū)塊的油-氣-水三相相對滲透率曲 線。其中，圖6為油-水兩相系統(tǒng)的油相相對滲透率曲線以及水相相對滲透率曲線；圖7為油-氣兩相系統(tǒng)的油相相對滲透率曲線以及氣相相對滲透率曲線；圖8為三相滲流條件下的油 相相對滲透率曲線。
[0124] 根據(jù)上述描述可以看出，本申請實(shí)施例通過獲取在目標(biāo)區(qū)域的氣-水交替非混相 驅(qū)替實(shí)驗(yàn)中所采集的不同時刻的驅(qū)替壓差、含水率以及氣油比等動態(tài)數(shù)據(jù)的觀測值以及基 于所構(gòu)建的組分?jǐn)?shù)值模擬模型所獲取的動態(tài)數(shù)據(jù)預(yù)測值，然后利用這些數(shù)據(jù)構(gòu)建目標(biāo)函 數(shù)，并對所構(gòu)建的目標(biāo)函數(shù)進(jìn)行迭代計(jì)算，得到滿足預(yù)設(shè)的迭代收斂條件的目標(biāo)函數(shù)，最后 根據(jù)所得到的目標(biāo)函數(shù)，獲取目標(biāo)區(qū)域的油-氣-水三相相對滲透率曲線，從而實(shí)現(xiàn)了提供 適用于氣-水交替非混相驅(qū)油技術(shù)的油-氣-水三相相對滲透率曲線的目的，這可以為描述 孔隙介質(zhì)中氣-水交替非混相驅(qū)替的多相滲流特征提供理論依據(jù)。
[0125] 本申請實(shí)施例還提供了一種獲取三相相對滲透率曲線的裝置，如圖9所示。該裝置 包括:第一獲取單元410、第二獲取單元420、迭代計(jì)算單元430以及第三獲取單元440。其中， 第一獲取單元410可以用于基于目標(biāo)區(qū)域的氣-水交替非混相驅(qū)替實(shí)驗(yàn)，獲取不同時刻的驅(qū) 替壓差、含水率以及氣油比，所述不同時刻的驅(qū)替壓差、含水率以及氣油比構(gòu)成了動態(tài)數(shù)據(jù) 觀測值向量;第二獲取單元420可以用于基于所構(gòu)建的表征模型所表征的先驗(yàn)三相相對滲 透率曲線以及所構(gòu)建的組分?jǐn)?shù)值模擬模型，獲取當(dāng)前動態(tài)數(shù)據(jù)預(yù)測值向量;迭代計(jì)算單元 430可以用于利用所述動態(tài)數(shù)據(jù)觀測值向量和所述當(dāng)前動態(tài)數(shù)據(jù)預(yù)測值向量，對所構(gòu)建的 目標(biāo)函數(shù)進(jìn)行迭代計(jì)算，直到所述目標(biāo)函數(shù)的迭代計(jì)算滿足預(yù)設(shè)的迭代收斂條件;第三獲 取單元440可以用于根據(jù)滿足所述迭代收斂條件的目標(biāo)函數(shù)，獲取目標(biāo)區(qū)域的油-氣-水三 相相對滲透率曲線。
[0126] 在至少一實(shí)施例中，第二獲取單元420可以包括(圖中未示出）：
[0127] 獲取子單元，可以用于獲取先驗(yàn)控制參數(shù)向量集合，所述先驗(yàn)控制參數(shù)向量集合 包括多個相對滲透率先驗(yàn)?zāi)Ｐ退鶎?yīng)的先驗(yàn)控制參數(shù)向量；
[0128] 組分?jǐn)?shù)值模擬子單元，可以用于基于所構(gòu)建的組分?jǐn)?shù)值模擬模型，利用所述先驗(yàn) 控制參數(shù)向量集合進(jìn)行組分?jǐn)?shù)值模擬，以得到當(dāng)前時間步的所述多個相對滲透率先驗(yàn)?zāi)Ｐ?所對應(yīng)的先驗(yàn)動態(tài)數(shù)據(jù)預(yù)測值向量；
[0129] 第一計(jì)算子單元，可以用于計(jì)算每個所述先驗(yàn)動態(tài)數(shù)據(jù)預(yù)測值向量與所述動態(tài)數(shù) 據(jù)觀測值向量之間的差值，以得到校正控制參數(shù)向量集合；
[0130] 第二計(jì)算子單元，可以用于計(jì)算所述校正控制參數(shù)向量集合所對應(yīng)的先驗(yàn)動態(tài)數(shù) 據(jù)預(yù)測值向量之間的平均值，以得到當(dāng)前動態(tài)數(shù)據(jù)預(yù)測值向量，
[0131] 相應(yīng)的，迭代計(jì)算單元430可以包括：
[0132] 第三計(jì)算子單元，可以用于利用所述動態(tài)數(shù)據(jù)觀測值向量和所述當(dāng)前動態(tài)數(shù)據(jù)預(yù) 測值向量，計(jì)算所述目標(biāo)函數(shù)的當(dāng)前數(shù)值；
[0133] 判斷子單元，可以用于判斷當(dāng)前迭代計(jì)算是否滿足迭代收斂條件；
[0134] 第一控制子單元，可以用于在判斷出當(dāng)前迭代計(jì)算不滿足迭代收斂條件時，將所 述校正控制參數(shù)向量集合作為先驗(yàn)控制參數(shù)向量集合，并控制所述組分?jǐn)?shù)值模擬子單元、 所述第一計(jì)算子單元以及所述第二計(jì)算子單元進(jìn)行相應(yīng)的操作，即控制所述組分?jǐn)?shù)值模擬 子單元基于所構(gòu)建的組分?jǐn)?shù)值模擬模型，利用所述先驗(yàn)控制參數(shù)向量集合進(jìn)行組分?jǐn)?shù)值模 擬，以得到當(dāng)前時間步的所述多個相對滲透率先驗(yàn)?zāi)Ｐ退鶎?yīng)的先驗(yàn)動態(tài)數(shù)據(jù)預(yù)測值向 量，控制所述第一計(jì)算子單元計(jì)算每個所述先驗(yàn)動態(tài)數(shù)據(jù)預(yù)測值向量與所述動態(tài)數(shù)據(jù)觀測 值向量之間的差值，以得到校正控制參數(shù)向量集合；以及控制所述第二計(jì)算子單元計(jì)算所 述校正控制參數(shù)向量集合所對應(yīng)的先驗(yàn)動態(tài)數(shù)據(jù)預(yù)測值向量的平均值，以得到當(dāng)前動態(tài)數(shù) 據(jù)預(yù)測值向量等操作；
[0135] 第二控制子單元，用于控制所述第三計(jì)算子單元以及所述判斷子單元的操作，具 體的，用于控制第三計(jì)算子單元利用所述動態(tài)數(shù)據(jù)觀測值向量和所述當(dāng)前動態(tài)數(shù)據(jù)預(yù)測值 向量來計(jì)算所述目標(biāo)函數(shù)的當(dāng)前數(shù)值，以及控制判斷子單元判斷當(dāng)前迭代計(jì)算是否滿足迭 代收斂條件，直到得到滿足迭代收斂條件的目標(biāo)函數(shù)。
[0136] 在至少一實(shí)施例中，第三獲取單元440可以包括(圖中未示出）：
[0137] 確定子單元，用于根據(jù)所述符合迭代收斂條件的目標(biāo)函數(shù)，確定出所述控制參數(shù) 向量中各控制參數(shù)的當(dāng)前數(shù)值；
[0138] 第四計(jì)算子單元，用于根據(jù)所述控制參數(shù)向量中各控制參數(shù)的當(dāng)前數(shù)值，計(jì)算油-水兩相系統(tǒng)的油相相對滲透率、水相相對滲透率以及油-氣兩相系統(tǒng)的氣相相對滲透率以 及油相相對滲透率；
[0139] 第五計(jì)算子單元，用于利用所計(jì)算的油-水兩相系統(tǒng)的油相相對滲透率、水相相對 滲透率以及油-氣兩相系統(tǒng)的氣相相對滲透率以及油相相對滲透率，計(jì)算三相滲流條件下 的油相相對滲透率。
[0140] 上述實(shí)施例闡明的系統(tǒng)、裝置或單元，具體可以由計(jì)算機(jī)芯片或?qū)嶓w實(shí)現(xiàn)，或者由 具有某種功能的產(chǎn)品來實(shí)現(xiàn)。
[0141]為了描述的方便，描述以上裝置時以功能分為各種單元分別描述。當(dāng)然，在實(shí)施本 申請時可以把各單元的功能在同一個或多個軟件和/或硬件中實(shí)現(xiàn)。
[0142] 本發(fā)明實(shí)施例中所描述的方法或算法的步驟可以直接嵌入硬件、處理器執(zhí)行的軟 件模塊、或者這兩者的結(jié)合。在一個或多個示例性的設(shè)計(jì)中，本發(fā)明實(shí)施例所描述的上述功 能可以在硬件、軟件、固件或這三者的任意組合來實(shí)現(xiàn)。
[0143] 本說明書中的各個實(shí)施例均采用遞進(jìn)的方式描述，各個實(shí)施例之間相同相似的部 分互相參見即可，每個實(shí)施例重點(diǎn)說明的都是與其他實(shí)施例的不同之處。尤其，對于系統(tǒng)實(shí) 施例而言，由于其基本相似于方法實(shí)施例，所以描述的比較簡單，相關(guān)之處參見方法實(shí)施例 的部分說明即可。
[0144] 雖然通過實(shí)施例描繪了本申請，本領(lǐng)域普通技術(shù)人員知道，本申請有許多變形和 變化而不脫離本申請的精神，希望所附的權(quán)利要求包括這些變形和變化而不脫離本申請的 精神。