nx-xT/nxtr)2+ ω 2 (yu/ny-yT/nytr)2+ ω 3 (zu/nz-zT/nztr)2；
[0039]6)選擇其中最小距離值所對應(yīng)的數(shù)據(jù)事件作為待估點處模擬結(jié)果，評估公式為:
[0040]Simulated (xu，yu, zu) = min{fi (xT, yT, zT), i = I, m}；
[0041]步驟四、整體替換數(shù)據(jù)事件:
[0042]將選擇的數(shù)據(jù)事件整體替換掉待估點處的數(shù)據(jù)事件，模擬指向下一個節(jié)點，重復步驟三和步驟四；
[0043]步驟五、建立模型:
[0044]直到所有待估點都完成預測，獲得對應(yīng)的模擬值，生成最終三維地質(zhì)模型。
[0045]實施例:
[0046]下面就某工區(qū)三角洲前緣基于位置的多點地質(zhì)統(tǒng)計學建模對本發(fā)明作進一步說明。
[0047]本實施例的數(shù)據(jù)為井數(shù)據(jù)以及網(wǎng)格化的工區(qū)。
[0048]工區(qū)規(guī)模:規(guī)格為columns (列數(shù))*rows(行數(shù))*layers(垂向網(wǎng)格)，具體為(columns = 199，rows = 199，layers = 10)，原始網(wǎng)格點數(shù)目為 396010 個。其 X 坐標值的間距為 cel I size IX (cel I size IX = 10)，Y 坐標值的間距為 cellsizelY (cellsizelY =10)，Z坐標值的間距為cellsizelZ(cellsizelZ = I)，而且X，Y，Z坐標值均只有一個起點值(5，5，0.5)0
[0049]網(wǎng)格按照標準的矩形網(wǎng)格體設(shè)計。通過起點和間隔就可以算出所有網(wǎng)格中心的三維坐標值。
[0050]1、工區(qū)地質(zhì)數(shù)據(jù)庫的建立:
[0051]參見圖3，根據(jù)研宄區(qū)井資料以及其范圍，將工區(qū)劃分成網(wǎng)格(199*199*10)的網(wǎng)格，并獲得每個網(wǎng)格點的三維坐標，以及每口井的沉積相類型。
[0052]2、訓練圖像的建立:
[0053]參見圖4，根據(jù)井資料沉積相數(shù)據(jù)計算了變差函數(shù)。利用截斷高斯模擬進行了沉積相空間分布的預測。然后結(jié)合手工編輯，形成研宄區(qū)三維訓練圖像。研宄區(qū)的三維訓練圖像與實際地形所編輯的二維圖像具有很好的一致性，反映了所建立的訓練圖像能夠有效反映儲層結(jié)構(gòu)特征。
[0054]3、相對距離計算:
[0055]參見圖5，首先是利用待估點周圍條件數(shù)據(jù)對訓練圖像進行掃描，確定滿足條件的數(shù)據(jù)事件，并確定每一個數(shù)據(jù)事件中心節(jié)點的位置。其次，選擇相對位置，一般我們選擇模擬區(qū)域和訓練圖像最大坐標點作為相對位置。在模擬區(qū)域中，最大坐標點坐標為(nx，ny,nz)，在訓練圖像中最大坐標點記為(nxtr，nytr，nztr)。這樣對于待估點坐標(xu，yu，zu)，則待估點在模擬區(qū)域的相對坐標為(xu/nx，yu/ny，zu/nz)，對于訓練圖像中的數(shù)據(jù)事件(χτ，yT，zT)，數(shù)據(jù)事件在訓練圖像中的相對坐標為(xT/nxtr，yT/nytr，zT/nztr)。隨后，計算待估點與數(shù)據(jù)事件的距離。其計算公式采用歐式距離:
[0056]f (xT, yT, ζτ) = ω J (xu/nx_xT/nxtr) 2+ω 2 (yu/ny_yT/nytr) 2+ω 3 (zu/nz_zT/nztr) 2O
[0057]最后，對所有滿足條件的數(shù)據(jù)事件的歐式距離做一個排序，選擇其中距離最小的數(shù)據(jù)事件作為模擬結(jié)果，具體公式如下:
[0058]Simulated (xu，yu, zu) = min {fj (xT, yT, zT), i = I, m}。
[0059]4、數(shù)據(jù)事件整體替換:
[0060]參見圖6，將選擇的數(shù)據(jù)事件投擲到待估點處，完成對待估點處儲層的預測。模擬將轉(zhuǎn)向下一個待估點。重復第3步和第4步。
[0061]5、模型建立:
[0062]將所有待估網(wǎng)格點模擬完成后，就可以將模擬結(jié)果進行顯示和輸出。如圖7所示，建立的地質(zhì)模型與觀測到現(xiàn)象較為接近，表明建立的模型能夠反映真實的地下地層情況。與傳統(tǒng)的多點統(tǒng)計建模方法相比，具有明顯優(yōu)勢。
[0063]顯然，本領(lǐng)域的技術(shù)人員可以對本發(fā)明進行各種改動和變型而不脫離本發(fā)明的精神和范圍。這樣，倘若本發(fā)明的這些修改和變型屬于本發(fā)明權(quán)利要求及其等同技術(shù)的范圍之內(nèi)，則本發(fā)明也意圖包含這些改動和變型在內(nèi)。
[0064]本說明書中未作詳細描述的內(nèi)容屬于本領(lǐng)域?qū)I(yè)技術(shù)人員公知的現(xiàn)有技術(shù)。
【主權(quán)項】
1.一種基于位置的多點地質(zhì)統(tǒng)計學建模方法，其特征在于:包括如下步驟: 步驟一、建立工區(qū)地質(zhì)數(shù)據(jù)庫； 步驟二、建立訓練圖像； 步驟三、計算相對距離； 步驟四、整體替換數(shù)據(jù)事件； 步驟五、建立模型。
2.根據(jù)權(quán)利要求1所述的基于位置的多點地質(zhì)統(tǒng)計學建模方法，其特征在于:所述步驟一中，根據(jù)地質(zhì)建模要求，將工區(qū)劃分成網(wǎng)格，獲得每一 口井在網(wǎng)格中的位置，以及其沉積相類型，建立地質(zhì)數(shù)據(jù)庫。
3.根據(jù)權(quán)利要求1所述的基于位置的多點地質(zhì)統(tǒng)計學建模方法，其特征在于:所述步驟二中，根據(jù)地質(zhì)知識庫的沉積相數(shù)據(jù)，計算各沉積相的變差函數(shù)，確定沉積相空間結(jié)構(gòu)特征和參數(shù)，然后，采用截斷高斯建模方法隨機插值，結(jié)合人工編輯獲得三維訓練圖像。
4.根據(jù)權(quán)利要求3所述的基于位置的多點地質(zhì)統(tǒng)計學建模方法，其特征在于:所述地質(zhì)知識庫為井資料。
5.根據(jù)權(quán)利要求1所述的基于位置的多點地質(zhì)統(tǒng)計學建模方法，其特征在于:所述步驟三中，過程如下: 1)從建模區(qū)域提取待估點的網(wǎng)格點坐標以及周圍條件數(shù)據(jù)； 2)利用條件數(shù)據(jù)掃描訓練圖像，尋找滿足條件數(shù)據(jù)的數(shù)據(jù)事件； 3)提取滿足條件數(shù)據(jù)的數(shù)據(jù)事件中心網(wǎng)格節(jié)點位置坐標； 4)分別在模擬區(qū)域和訓練圖像中選擇同一個參照點，獲得待估點在模擬區(qū)域相對位置，以及數(shù)據(jù)事件在訓練圖像中的相對位置； 5)采用歐式距離計算數(shù)據(jù)事件相對位置與待估點相對位置的距離； 6)選擇其中最小距離值所對應(yīng)的數(shù)據(jù)事件作為待估點處模擬結(jié)果。
6.根據(jù)權(quán)利要求5所述的基于位置的多點地質(zhì)統(tǒng)計學建模方法，其特征在于:所述步驟三第4)項中，模擬區(qū)域參照點的坐標為(nx，ny，nz)，所述參照點在訓練圖像中的坐標記為(nxtr，nytr, nztr)，則待估點的坐標為(xu，yu, zu)，所述待估點在模擬區(qū)域的相對坐標為(xu/nx，yu/ny，zu/nz)，對于訓練圖像中的數(shù)據(jù)事件(xT，yT，zT)，所述數(shù)據(jù)事件在訓練圖像中的相對坐標為(xT/nxtr，yT/nytr, zT/nztr)。
7.根據(jù)權(quán)利要求6所述的基于位置的多點地質(zhì)統(tǒng)計學建模方法，其特征在于:所述步驟三第5)項中，數(shù)據(jù)事件相對位置與待估點相對位置的歐氏距離為: f(xT，yT, ζτ) = oJxu/nx-XT/nxtrV+G^yu/ny-yT/nytrV+o^Zu/nz-ZT/nztr)、
8.根據(jù)權(quán)利要求7所述的基于位置的多點地質(zhì)統(tǒng)計學建模方法，其特征在于:所述步驟三第6)項中，評估公式為:Simulated (xu，yu, zu) = min {fj (xT，yT, zT), i = I, m}。
9.根據(jù)權(quán)利要求1所述的基于位置的多點地質(zhì)統(tǒng)計學建模方法，其特征在于:所述步驟四中，將選擇的數(shù)據(jù)事件整體替換掉待估點處的數(shù)據(jù)事件，模擬指向下一個節(jié)點，重復步驟三和步驟四。
10.根據(jù)權(quán)利要求1所述的基于位置的多點地質(zhì)統(tǒng)計學建模方法，其特征在于:所述步驟五中，直到所有待估點都完成預測，獲得對應(yīng)的模擬值，生成最終三維地質(zhì)模型。
【專利摘要】本發(fā)明公開了一種基于位置的多點地質(zhì)統(tǒng)計學建模方法，包括如下步驟：步驟一、建立工區(qū)地質(zhì)數(shù)據(jù)庫；步驟二、建立訓練圖像；步驟三、計算相對距離；步驟四、整體替換數(shù)據(jù)事件；步驟五、建立模型。本發(fā)明解決了訓練圖像與模擬區(qū)域尺度差異問題，方法運用靈活，提高了三角洲前緣儲層精細預測精度，獲得了更準確的三維儲層地質(zhì)模型，可以廣泛應(yīng)用于油氣勘探開發(fā)技術(shù)領(lǐng)域。
【IPC分類】G06F17-50, G06T17-05
【公開號】CN104850682
【申請?zhí)枴緾N201510184109
【發(fā)明人】尹艷樹, 張昌民, 王軍, 李少華
【申請人】長江大學
【公開日】2015年8月19日
【申請日】2015年4月17日