下面參考根據(jù)本發(fā)明的各實施例的方法、設(shè)備(系統(tǒng))和計算機程序產(chǎn)品的流程圖 圖示和/或框圖來描述本發(fā)明的各方面。應當理解的是，可W通過計算機程序指令實現(xiàn)流程 圖圖示和/或框圖的每個塊，W及流程圖圖示和/或框圖中的塊的組合?？蒞將運些計算機 程序指令提供至通用計算機、專用計算機的處理器或者其它可編程數(shù)據(jù)處理設(shè)備，W產(chǎn)生 機器，使得經(jīng)由計算機的處理器或其它可編程數(shù)據(jù)處理設(shè)備執(zhí)行的指令產(chǎn)生用于實現(xiàn)流程 圖和/或框圖的一個塊或多個塊中指定的功能/動作的裝置。
[0027] 運些計算機程序指令也可W存儲在能夠指揮計算機、其它可編程數(shù)據(jù)處理設(shè)備或 其它裝置W特定方式運行的計算機可讀介質(zhì)中，使得存儲在計算機可讀介質(zhì)中的指令產(chǎn)生 包括實現(xiàn)流程圖和/或框圖的一個塊或多個塊中指定的功能/動作的指令的制品。
[0028] 計算機程序指令也可W加載到計算機、其它可編程數(shù)據(jù)處理設(shè)備或其它裝置上， W致使在計算機、其它可編程設(shè)備或其它裝置上執(zhí)行一系列操作步驟來產(chǎn)生計算機實現(xiàn)的 過程，使得在計算機或其它可編程設(shè)備上執(zhí)行的指令提供用于實現(xiàn)流程圖和/或框圖的一 個塊或多個塊中指定的功能/動作的過程。
[0029] 現(xiàn)在轉(zhuǎn)到附圖，并且更具體地，圖1示出了根據(jù)本發(fā)明的優(yōu)選實施例的用于評估的 預期場的示例，該示例W匯編關(guān)于目標儲層的信息為開始。該信息可被可W根據(jù)該所采集 的信息提供數(shù)目η個相的技術(shù)人員理解。相的數(shù)目優(yōu)選地在1與8之間。圖1示出了在其表面 的特定部分下方具有域(D)的場，用字母A表示該區(qū)域。該圖示出了具有分布在域內(nèi)的Ξ個 相Fi，F(xiàn)2,的的域(D)。
[0030] 在缺少來源于地震或地圖的任何信息的情況下，儲層域(D)被視為六面體形狀?？?慮關(guān)于形狀和體積兩者的結(jié)構(gòu)不確定性。在目前的情形下，由于概念模型是規(guī)則形狀，所W 用儲層面積、厚度W及寬長比來表示不確定性。
[0031] 已經(jīng)研究出插值函數(shù)來根據(jù)類比數(shù)據(jù)庫估計儲層面積和厚度的聯(lián)合分布，然而， 假設(shè)寬長比參數(shù)在其各邊界之間均勻變化。數(shù)值分析已經(jīng)示出儲層面積和厚度可W近似于 具有高置信區(qū)間的二元對數(shù)正態(tài)分布。
[0032] 重要的是強調(diào)面積和厚度的聯(lián)合概率是關(guān)鍵，W保持根據(jù)類比物關(guān)于體積和儲量 給出的置信區(qū)間。如果存在關(guān)于儲層的石油的總量的信息或者可W估計該量，則必選選擇 大到足W容納石油的量的域。
[0033] 根據(jù)本發(fā)明的一個實施例，可W通過將目標信息與類比數(shù)據(jù)庫進行比較來得到相 的數(shù)目和其它特性。
[0034] 通過每個相的偏率flW及具有均值吳和用于預定闊值的截斷值咨'in和終-的孔隙 率的正態(tài)分布來模擬每個相Fi，i = l...n，其中，預定闊值使得正態(tài)分布的超過運樣的闊值 的尾部區(qū)域被認為是可忽略的。一旦預先確定了闊值，則使正態(tài)分布被限制在祭in和終?之 間，并且均值滿足在=(《胃-終li。)/2。偏率村馬足玄'=1乂 =1，其中，η是已固定的相的數(shù)目。
[0035] 如圖3所示，通過可W由技術(shù)人員提供的《尹、和義的值，或者通過可W從 類比數(shù)據(jù)庫得到的估計，或者更簡要地通過"類比物"來模擬整個域。
[0036] 下一步驟計算疼、終…、黎U% ?' = 1...?為最小化問題的解，其中，由W下表 達式定義成本函數(shù)：
[0037] >
[0038] 其中，由^下條件定義後^'和《識：
[0041 ] 并且，取《使得否 < 苗I /二i-U
[0042] 在一個實施例中，范數(shù)II . II已經(jīng)被實現(xiàn)為歐幾里得范數(shù)。要被最小化的目標函 數(shù)表示作為相類型的函數(shù)的所估計特性與例如通過類比物提供的值之間的誤差。
[0043] 圖3示出了四個相(n = 4)的情形，其中，根據(jù)成本函數(shù)的最小值來分布孔隙率的四 個正態(tài)分布;并且，同時確定了偏率fi，i = l. . .4。
[0044] 可W使用混合整數(shù)線性規(guī)劃(MILP)來求解該最小值問題。該求解方法聯(lián)合計算每 個相的相對比例W及每個相的特性的條件分布。該估計將與用作參考(ref)的:^、病 和充;f的最小值(min)、最大值(max) W及平均值一致。
[0045] 在本發(fā)明的一個實施例中，從類比數(shù)據(jù)庫(2)得到參考值每15"、#^和哀^。當匯編 關(guān)于目標（1)油氣儲層的信息時，所述信息至少具有類比數(shù)據(jù)庫的記錄中包括的材料特性 中的一個材料特性。W如下方式確定參考值：
[0046] 1.當提供包括已知儲層的記錄的類比數(shù)據(jù)庫（2)時，至少一個記錄還包括材料的 孔隙率，即孔隙率的最小值然:"，孔隙率的均值運。W及孔隙率的最大值cr;
[0047] 運Ξ個值允許假設(shè)諸如正態(tài)分布的特性分布的形狀?？紫堵蕦⒂糜谟嬎銋⒖贾担?作為具有運樣的特性的類比數(shù)據(jù)的記錄的子集的均值。
[0048] 2.當在步驟a)匯編關(guān)于目標油氣儲層的信息時，運樣的信息具有類比數(shù)據(jù)庫的記 錄所包括的材料特性中的至少一個材料特性；
[0049] 目標油氣儲層中W及類比數(shù)據(jù)庫中的共同特性允許將目標(1)信息與密切相關(guān)的 類比記錄(2.1)進行比較。共同特性可W與孔隙率不同。
[0050] 3.提供用于至少共享共同材料特性的儲層的比較的相似性準則；
[0051] 相似性準則是定義類比數(shù)據(jù)庫的哪些記錄接近于目標（1)的準則。該準則甚至可 W組合特性或者兩個或更多個特性之間的相關(guān)性。
[0052] 4.在類比數(shù)據(jù)庫(2)中選擇與目標儲層滿足所述相似性準則的m個記錄(2.1)的集 合，其中，所選記錄中的至少一個記錄包括材料的孔隙率；
[0053] 5.確定為；·'、和衣,/的值時，其中，針對在步驟4)中選擇的類比數(shù)據(jù)庫的記錄 的集合得到下述值：
[0056]
，針對具有孔隙率信息的所有記錄。
[0057] 在本發(fā)明的一個實施例中，得到具有Ξ維域的多個Ξ維場景，隨機生成所述域的 維度。如果類比物具有諸如儲層面積和厚度的體積數(shù)據(jù)、NTG(凈毛比）和飽和度信息，則計 算出在域內(nèi)存儲石油的最大容量。如果隨機生成的域(D)不大到足W存儲所計算的石油的 體積，則丟棄所述域(D)。
[005引例如，00IP(原始石油地質(zhì)儲量)可W被估計為：
[0化9]
[0060]其中，V是域的體積，Φ是孔隙率的均值，NTG是凈毛比的均值，S是飽和度，W及B是 當從儲層提取石油并W環(huán)境壓力存儲時石油的膨脹系數(shù)。
[0061 ] 一旦已經(jīng)生成了各相，并且確定了fi、充、病胃、：?Τ% ? = ι...Μ變量，則通過提 供預定縱向相關(guān)性方程和預定取向的地質(zhì)統(tǒng)計算法，基于孔隙率的統(tǒng)計分布 是、《nin、來空間分布孔隙率。在運種情形下，在多個場景的每個場景中定義孔隙率 的空間分布。
[0062] 在根據(jù)孔隙率的統(tǒng)計分布、縱向相關(guān)性方程W及取向填充孔隙率之后，針對每個 場景確定滲透率特性化）。
[0063] 為了執(zhí)行該進一步步驟，類比數(shù)據(jù)庫（2)的記錄必須包括材料的滲透率特性和孔 隙率。根據(jù)所選記錄，計算出滲透率與孔隙率之間的相關(guān)性，其中，在大多數(shù)情形下，根據(jù)卡 曼-科澤尼定律，兩個變量之間的定律被表示成滲透率與孔隙率之間的對數(shù)值。所選記錄是 在從類比數(shù)據(jù)庫(2)得到參考值探;*'、和哀/時選擇的那些記錄。
[0064] 根據(jù)卡曼-科澤尼定律，方程
[00 化]
[0066] 用于擬合來自于滲透率化)與孔隙率（Φ)之間的關(guān)系的散布。可W使用最小二乘 回歸來估計系數(shù)C。為了使分布的低散布成為可能，可W對類比結(jié)果應用另外的濾波器，諸 如利用具有如可能為深度的其它變量的可能的相關(guān)性。還使用將線性回歸模型用作孔隙率 的函數(shù)來用參數(shù)表示滲透性方差，W正確傳遞Ξ維實現(xiàn)的不確定性。為了得到可靠的模擬， 重要的是，分離用孔隙率確定的均值(趨勢)和通過地質(zhì)統(tǒng)計技術(shù)模擬的方差的作用