一種獲取三相相對滲透率曲線的方法及裝置的制造方法
【技術領域】
[0001 ]本申請涉及油氣藏開發(fā)技術領域，特別涉及一種獲取油-氣-水三相相對滲透率曲 線的方法及裝置。
【背景技術】
[0002] 相對滲透率曲線是油氣藏開發(fā)中的一項重要數(shù)據(jù)，它能夠精細描述孔隙介質中的 多相滲流特征。目前，油-水或氣-液兩相相對滲透率曲線主要通過巖心驅替實驗獲取，計算 方法主要以Johnson-Bossler-Naumann(JBN)等解析方法以及數(shù)值反演方法為主。而對獲取 油-氣-水三相相對滲透率曲線的方法的研究相對較少。
[0003] 隨著油藏縱向及平面非均質矛盾的加劇，中國注水開發(fā)油田主體已普遍進入高含 水、高采出程度的"雙高"開發(fā)階段，剩余油呈現(xiàn)"整體高度分散、局部相對富集"的分布格 局。室內實驗及礦場實踐表明，氣-水交替非混相驅是水驅油藏后期進一步提高原油采收率 的有效技術。因此，研究適用于氣-水交替非混相驅油技術的獲取油-氣-水三相相對滲透率 曲線的方法顯得意義重大，這可以為準確描述孔隙介質中氣-水交替非混相驅替的多相滲 流特征提供理論依據(jù)。

【發(fā)明內容】

[0004] 本申請實施例的目的是提供一種獲取三相相對滲透率曲線的方法及裝置，以實現(xiàn) 獲取適用于氣-水交替非混相驅油技術的油-氣-水三相相對滲透率曲線的目的。
[0005] 為解決上述技術問題，本申請實施例提供一種獲取三相相對滲透率曲線的方法及 裝置是這樣實現(xiàn)的：
[0006] 本申請實施例提供了一種獲取三相相對滲透率曲線的方法，包括：
[0007] S1，基于目標區(qū)域的氣-水交替非混相驅替實驗，獲取不同時刻的驅替壓差、含水 率以及氣油比，所述不同時刻的驅替壓差、含水率以及氣油比構成了動態(tài)數(shù)據(jù)觀測值向量；
[0008] S2，基于所構建的表征模型所表征的先驗三相相對滲透率曲線以及所構建的組分 數(shù)值模擬模型，獲取當前動態(tài)數(shù)據(jù)預測值向量；
[0009] S3,利用所述動態(tài)數(shù)據(jù)觀測值向量和所述當前動態(tài)數(shù)據(jù)預測值向量，對所構建的 目標函數(shù)進行迭代計算，直到所述目標函數(shù)的迭代計算滿足預設的迭代收斂條件；
[0010] S4,根據(jù)滿足所述迭代收斂條件的目標函數(shù)，獲取所述目標區(qū)域的油-氣-水三相 相對滲透率曲線。
[0011] 在至少一實施例中，所述組分數(shù)值模擬模型基于所確定的流體物性參數(shù)分布特征 來構建，所述流體物性參數(shù)特征包括飽和壓力、原油密度、粘度和體積系數(shù)的變化規(guī)律。
[0012] 在至少一實施例中，所述步驟S2包括：
[0013] S21，獲取先驗控制參數(shù)向量集合，所述先驗控制參數(shù)向量集合包括多個相對滲透 率先驗模型所對應的先驗控制參數(shù)向量；
[0014] S22,基于所構建的組分數(shù)值模擬模型，利用所述先驗控制參數(shù)向量集合進行組分 數(shù)值模擬，以得到當前時間步的所述多個相對滲透率先驗模型所對應的先驗動態(tài)數(shù)據(jù)預測 值向量；
[0015] S23,計算每個所述先驗動態(tài)數(shù)據(jù)預測值向量與所述動態(tài)數(shù)據(jù)觀測值向量之間的 差值或誤差平方和，以得到符合預設要求的校正控制參數(shù)向量集合；
[0016] S24,計算所述校正控制參數(shù)向量集合所對應的先驗動態(tài)數(shù)據(jù)預測值向量的平均 值，以得到當前動態(tài)數(shù)據(jù)預測值向量，
[0017]相應的，所述步驟S3包括：
[0018] S31，利用所述動態(tài)數(shù)據(jù)觀測值向量和所述當前動態(tài)數(shù)據(jù)預測值向量，計算所述目 標函數(shù)的當前數(shù)值；
[0019] S32,判斷當前迭代計算是否滿足迭代收斂條件；
[0020] S33,在判斷出當前迭代計算不滿足迭代收斂條件時，將所述校正控制參數(shù)向量集 合作為先驗控制參數(shù)向量集合，重復步驟S22-S24,得到當前動態(tài)數(shù)據(jù)預測值向量；
[0021 ] S34，在得到所述校正控制參數(shù)向量集合所對應的當前動態(tài)數(shù)據(jù)預測值向量后，重 復步驟S31-32,直到得到滿足所述迭代收斂條件的目標函數(shù)。
[0022] 在至少一實施例中，所構建的目標函數(shù)的表達式如下：
[0023]
[0024] 其中，0(m)為目標函數(shù);m為所述表征模型中的控制參數(shù)向量;T為表征向量或矩陣 轉置的符號;cUs為動態(tài)數(shù)據(jù)觀測值向量;g(m)為動態(tài)數(shù)據(jù)預測值向量;C D為動態(tài)數(shù)據(jù)測量 誤差的協(xié)方差矩陣;mpricir為所述控制參數(shù)向量m的先驗模型信息的平均值;C m為先驗模型的 協(xié)方差矩陣。
[0025] 在至少一實施例中，所述根據(jù)滿足所述迭代收斂條件的目標函數(shù)，獲取所述目標 區(qū)域的油-氣-水三相相對滲透率曲線包括：
[0026] 根據(jù)滿足所述迭代收斂條件的目標函數(shù)，確定出所述控制參數(shù)向量中各控制參數(shù) 的當前數(shù)值；
[0027] 根據(jù)所述控制參數(shù)向量中各控制參數(shù)的當前數(shù)值，計算油-水兩相系統(tǒng)的油相相 對滲透率、水相相對滲透率以及油-氣兩相系統(tǒng)的氣相相對滲透率以及油相相對滲透率；
[0028] 利用所計算的油-水兩相系統(tǒng)的油相相對滲透率、水相相對滲透率以及油-氣兩相 系統(tǒng)的氣相相對滲透率以及油相相對滲透率，計算三相滲流條件下的油相相對滲透率。 [0029]在至少一實施例中，通過以下公式來計算三相滲流條件下的油相相對滲透率：
[0030] Kro= (Krow+Krw) (Krog+Krg)_(Krw+Krg)
[0031] 式中，Kr。為三相滲流條件下的油相相對滲透率;Kr?為油-水兩相系統(tǒng)的油相相對 滲透率;Krw為油-水兩相系統(tǒng)的水相相對滲透率;K rcig為油-氣兩相系統(tǒng)的油相相對滲透率； Krg為油-氣兩相系統(tǒng)的氣相相對滲透率。
[0032]本申請實施例還提供了一種獲取三相相對滲透率曲線的裝置，包括：
[0033]第一獲取單元，用于基于目標區(qū)域的氣-水交替非混相驅替實驗，獲取不同時刻的 驅替壓差、含水率以及氣油比，所述不同時刻的驅替壓差、含水率以及氣油比構成了動態(tài)數(shù) 據(jù)觀測值向量；
[0034] 第二獲取單元，用于基于所構建的表征模型所表征的先驗三相相對滲透率曲線以 及所構建的組分數(shù)值模擬模型，獲取當前動態(tài)數(shù)據(jù)預測值向量；
[0035] 迭代計算單元，用于利用所述動態(tài)數(shù)據(jù)觀測值向量和所述當前動態(tài)數(shù)據(jù)預測值向 量，對所構建的目標函數(shù)進行迭代計算，直到所述目標函數(shù)的迭代計算滿足預設的迭代收 斂條件；
[0036] 第三獲取單元，用于根據(jù)滿足所述迭代收斂條件的目標函數(shù)，獲取所述目標區(qū)域 的油-氣-水三相相對滲透率曲線。
[0037] 在至少一實施例中，所述第二獲取單元包括：
[0038] 獲取子單元，用于獲取先驗控制參數(shù)向量集合，所述先驗控制參數(shù)向量集合包括 多個相對滲透率先驗模型所對應的先驗控制參數(shù)向量；
[0039] 組分數(shù)值模擬子單元，用于基于所構建的組分數(shù)值模擬模型，利用所述先驗控制 參數(shù)向量集合進行組分數(shù)值模擬，以得到當前時間步的所述多個相對滲透率先驗模型所對 應的先驗動態(tài)數(shù)據(jù)預測值向量；
[0040] 第一計算子單元，用于計算每個所述先驗動態(tài)數(shù)據(jù)預測值向量與所述動態(tài)數(shù)據(jù)觀 測值向量之間的差值或誤差平方和，以得到校正控制參數(shù)向量集合；
[0041] 第二計算子單元，用于計算所述校正控制參數(shù)向量集合所對應的先驗動態(tài)數(shù)據(jù)預 測值向量的平均值，以得到當前動態(tài)數(shù)據(jù)預測值向量，
[0042]相應的，迭代計算單元包括：
[0043]第三計算子單元，用于利用所述動態(tài)數(shù)據(jù)觀測值向量和所述當前動態(tài)數(shù)據(jù)預測值 向量，計算所述目標函數(shù)的當前數(shù)值；
[0044] 判斷子單元，用于判斷當前迭代計算是否滿足迭代收斂條件；
[0045] 第一控制子單元，用于在判斷出當前迭代計算不滿足迭代收斂條件時，將所述校 正控制參數(shù)向量集合作為先驗控制參數(shù)向量集合，并控制所述組分數(shù)值模擬子單元、所述 第一計算子單元以及所述第二計算子單元進行相應的操作，即