一種對(duì)復(fù)雜地質(zhì)構(gòu)造層序建模的方法
【專利摘要】公開(kāi)了一種對(duì)復(fù)雜地質(zhì)構(gòu)造層序建模的方法，該方法包括：步驟A：通過(guò)對(duì)地震解釋的層位、斷層散點(diǎn)數(shù)據(jù)進(jìn)行第一組段，來(lái)找出并記錄與各個(gè)層位段相關(guān)聯(lián)的斷層；步驟B：根據(jù)相關(guān)聯(lián)的層位、斷層建立塊狀模型；步驟C：通過(guò)對(duì)在步驟B劃分的小塊進(jìn)行處理來(lái)獲得網(wǎng)格化后的小層結(jié)果；步驟D：根據(jù)在步驟C獲得的小層結(jié)果輸出整體層序建模結(jié)果。
【專利說(shuō)明】一種對(duì)復(fù)雜地質(zhì)構(gòu)造層序建模的方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及油氣勘探開(kāi)發(fā)領(lǐng)域，更具體地說(shuō)，涉及一種基于測(cè)井層位(分層)和斷層雙重約束對(duì)復(fù)雜地質(zhì)構(gòu)造層序建模的方法，所述雙重約束指的是同時(shí)考慮層位信息和斷層Ih息。
【背景技術(shù)】
[0002]油藏地質(zhì)建模技術(shù)是地震勘探開(kāi)發(fā)中一項(xiàng)重要的油藏描述技術(shù)，對(duì)于儲(chǔ)層空間分布區(qū)域的識(shí)別、油藏參數(shù)的表征具有非常重要的作用。儲(chǔ)層地質(zhì)模型是油藏描述的關(guān)鍵，是儲(chǔ)層特征及其非均質(zhì)性在三維空間上變化和分布的表征。儲(chǔ)層建模貫穿油田開(kāi)發(fā)的全過(guò)程，不同的開(kāi)發(fā)階段，獲得的資料基礎(chǔ)和地質(zhì)認(rèn)識(shí)不同，儲(chǔ)層建模的重點(diǎn)內(nèi)容和精度不同，其所采用的儲(chǔ)層建模技術(shù)和方法也有很大差別。
[0003]一般而言，建立儲(chǔ)層模型有兩種基本途徑，即確定性建模和隨機(jī)建模。確定性建模是對(duì)井間未知區(qū)給出確定性的預(yù)測(cè)結(jié)果，確定性建模的思路是試圖從具有確定性資料的控制點(diǎn)(井點(diǎn))出發(fā)，推測(cè)出井間確定的、唯一的、真實(shí)的儲(chǔ)層參數(shù)。隨機(jī)建模，是指以己知信息為基礎(chǔ)，以隨機(jī)函數(shù)為理論，應(yīng)用隨機(jī)模擬方法產(chǎn)生可選的、等概率的儲(chǔ)層模型的方法。這種方法承認(rèn)控制點(diǎn)(井點(diǎn))以外的儲(chǔ)層參數(shù)具有一定的不確定性，即具有一定的隨機(jī)性。隨機(jī)建模不僅能反映觀測(cè)數(shù)據(jù)中包含的總體趨勢(shì)，還能反映局部變化。所建立的模型更加符合地質(zhì)體的實(shí)際情況和物性分布規(guī)律，從而克服了確定性建模帶來(lái)的“平滑”效應(yīng)。因此，隨機(jī)建?？蓪?duì)儲(chǔ)層非均質(zhì)性進(jìn)行高分辨率的表征。
[0004]層序建模方法是后續(xù)的屬性建模方法的核心，它的主要功能是通過(guò)某種方法在兩個(gè)已知的層位之間模擬出其中未知空間的大體構(gòu)造。即通過(guò)已知層位信息估算在無(wú)層位的位置空間位置的可能層位走勢(shì)，這些層位既受到其上下兩個(gè)已知層位的影響，又滿足與之相關(guān)的所有的斷層信息的約束。它模擬出了兩個(gè)層位之間的實(shí)際的地質(zhì)構(gòu)造走勢(shì)，為后續(xù)的模型構(gòu)建提供依據(jù)。層序建模是儲(chǔ)層建模過(guò)程中的一個(gè)關(guān)鍵環(huán)節(jié)，它決定了最終的油藏地質(zhì)模型的準(zhǔn)確性和精度。而對(duì)于復(fù)雜地質(zhì)構(gòu)造的油氣藏，準(zhǔn)確建立層序模型是一個(gè)難道很大的難題。因此，需要一種能夠方便快捷地對(duì)復(fù)雜地質(zhì)構(gòu)造層序建模的方法。

【發(fā)明內(nèi)容】

[0005]為了解決含逆斷層的復(fù)雜地質(zhì)構(gòu)造在層序建模過(guò)程中的小層劃分不清、小斷層層序走向不確定和網(wǎng)格化不清楚、尖滅點(diǎn)、上超點(diǎn)無(wú)法確定的問(wèn)題，本發(fā)明提出了采用測(cè)井地質(zhì)層位(分層)和斷層雙重約束的方法，建立準(zhǔn)確、合理的層序模型，提高了層序建模的效率，節(jié)省了油藏建模工作者的工作時(shí)間，并且得到的結(jié)果更加合理。
[0006]根據(jù)本發(fā)明的一方面，提供了一種對(duì)復(fù)雜地質(zhì)構(gòu)造層序建模的方法，該方法包括:步驟A:通過(guò)對(duì)地震解釋的層位散點(diǎn)數(shù)據(jù)、斷層散點(diǎn)數(shù)據(jù)進(jìn)行第一組段，來(lái)找出并記錄各個(gè)層位的關(guān)聯(lián)斷層，其中，所述關(guān)聯(lián)斷層是將所述各個(gè)層位中的相應(yīng)層位分割為兩個(gè)單獨(dú)層位段的斷層，并且所述第一組段操作包括對(duì)斷層散點(diǎn)數(shù)據(jù)進(jìn)行第一組段的操作和對(duì)層位散點(diǎn)數(shù)據(jù)進(jìn)行第一組段的操作；步驟B:根據(jù)彼此關(guān)聯(lián)的層位和斷層建立塊狀模型，其中，通過(guò)找出兩個(gè)層位之間的所有大斷層并將其排序，來(lái)根據(jù)排好序的大斷層以及工區(qū)邊界信息將組段后的層位段劃分成小塊，所述大斷層是與上下兩個(gè)層位相關(guān)聯(lián)的斷層；步驟C:通過(guò)對(duì)在步驟B劃分的小塊進(jìn)行處理來(lái)獲得網(wǎng)格化后的小層結(jié)果；步驟D:根據(jù)在步驟C獲得的小層結(jié)果輸出整體層序建模結(jié)果。
[0007]優(yōu)選地，所述步驟C包括:步驟1:判斷小塊內(nèi)部是否包括斷層；步驟2:如果小塊內(nèi)部沒(méi)有斷層，則對(duì)小塊進(jìn)行網(wǎng)格化；步驟3:如果小塊內(nèi)部包括斷層，則通過(guò)抹去內(nèi)部斷層操作及層位拉平操作，將內(nèi)部有斷層的小塊近似轉(zhuǎn)化為內(nèi)部無(wú)斷層的小塊；對(duì)小塊進(jìn)行網(wǎng)格化；對(duì)網(wǎng)格化后的小塊進(jìn)行劃分；返回步驟1，直到所有小塊中都不包括斷層；步驟4:輸出網(wǎng)格化后的小層結(jié)果。
[0008]優(yōu)選地，對(duì)小塊進(jìn)行網(wǎng)格化的步驟包括:根據(jù)小塊的上下層位段信息，得到小塊的上下層位段的起始及結(jié)束位置；根據(jù)所述起始及結(jié)束位置，將上下層位段按等比例劃分為相同段數(shù)的小段；分別連接上下層位段上相同比例處的分割點(diǎn)，連接上下層位段上相同比例處的分割點(diǎn)的直線是參考線；根據(jù)小層個(gè)數(shù)，將每條參考線進(jìn)行相應(yīng)的等比例劃分；通過(guò)連接相應(yīng)的參考線上的分割點(diǎn)，來(lái)獲得所需的小層。
[0009]優(yōu)選地，對(duì)小塊進(jìn)行劃分的步驟包括:斷層范圍計(jì)算、關(guān)鍵斷層選取、約束小層構(gòu)建和小塊構(gòu)建。
[0010]優(yōu)選地，對(duì)層位數(shù)據(jù)進(jìn)行第一組段的步驟包括:確定與各個(gè)層位相關(guān)聯(lián)的斷層；遍歷每個(gè)層位點(diǎn)，找出沿上下方向離該層位點(diǎn)最近的斷層點(diǎn)，獲取關(guān)聯(lián)斷層信息并進(jìn)行標(biāo)記；對(duì)確定了關(guān)聯(lián)斷層信息的層位點(diǎn)進(jìn)行第二組段，所述第二組段操作包括組小段和組大段；確定層位與斷層的關(guān)聯(lián)關(guān)系信息，并根據(jù)所述關(guān)聯(lián)關(guān)系信息對(duì)層位段進(jìn)行排序；其中，組小段操作包括:在所有層位點(diǎn)中遍歷并尋找到具有相同屬性的點(diǎn)；并將相同屬性的點(diǎn)組合為一個(gè)小段；組大段操作包括:如果小段的尾點(diǎn)與另一小段的首點(diǎn)相鄰，則將所述兩個(gè)小段組合為一個(gè)大段。
[0011]優(yōu)選地，所述第二組段操作還包括:當(dāng)尋找到特定小段時(shí)，遍歷其他小段，根據(jù)該特定小段的共深度點(diǎn)CDP信息判斷其他小段是否可能與該特定小段進(jìn)行組段；如果未找到能夠與當(dāng)前小段進(jìn)行組段的小段，則當(dāng)前大段組段結(jié)束，并開(kāi)始下一大段的組段工作；如果找到一個(gè)小段能夠與當(dāng)前的特定小段進(jìn)行組段，則將這兩個(gè)小段合并為一個(gè)大段；如果找到能與當(dāng)前小段進(jìn)行組段的多個(gè)小段，則首先讀取出當(dāng)前小段的尾點(diǎn)以及所有符合條件的小段的首點(diǎn)；通過(guò)對(duì)所有小段的時(shí)間值或深度值進(jìn)行比較，找出首點(diǎn)距離與當(dāng)前小段尾點(diǎn)最近的小段，并將該小段與當(dāng)前的特定小段進(jìn)行組段。
[0012]優(yōu)選地，根據(jù)關(guān)聯(lián)關(guān)系信息對(duì)層位段進(jìn)行排序的步驟包括:根據(jù)關(guān)聯(lián)關(guān)系信息對(duì)層位段從左到右進(jìn)行排序，其中，首先找到左關(guān)聯(lián)斷層為空的層位段，之后循環(huán)插值左關(guān)聯(lián)斷層與上一層位段的右關(guān)聯(lián)斷層相同的層位段，直到找到一個(gè)右關(guān)聯(lián)斷層為空的層位段為止。
[0013]根據(jù)本發(fā)明的另一方面，提供一種對(duì)復(fù)雜地質(zhì)構(gòu)造三維層序建模的方法，包括:將整個(gè)工區(qū)按共深度點(diǎn)CDP坐標(biāo)劃分為二維剖面；在每個(gè)剖面上單獨(dú)進(jìn)行如權(quán)利要求1所述的層序建模的方法；將每個(gè)剖面的層序建模結(jié)果組合成三維空間中的建模結(jié)果?！緦＠綀D】

【附圖說(shuō)明】
[0014]通過(guò)下面結(jié)合附圖進(jìn)行的詳細(xì)描述，本發(fā)明的上述和其它目的、特點(diǎn)和優(yōu)點(diǎn)將會(huì)變得更加清楚，其中:
[0015]圖1示出根據(jù)本發(fā)明實(shí)施例的層序建模的方法的流程圖；
[0016]圖2示出根據(jù)本發(fā)明實(shí)施例的組段的方法的流程圖；
[0017]圖3示出根據(jù)本發(fā)明實(shí)施例的關(guān)聯(lián)斷層的示意圖；
[0018]圖4示出根據(jù)本發(fā)明實(shí)施例的組段原理的示意圖；
[0019]圖5示出根據(jù)本發(fā)明實(shí)施例的層位組段的特殊情況的示意圖；
[0020]圖6示出根據(jù)本發(fā)明示例性實(shí)施例的確定關(guān)聯(lián)關(guān)系的示意圖；
[0021]圖7示出根據(jù)本發(fā)明實(shí)施例的斷層與層位之間的關(guān)系的特殊情況的示意圖；
[0022]圖8示出根據(jù)本發(fā)明實(shí)施例的成塊原理的示意圖；
[0023]圖9示出根據(jù)本發(fā)明實(shí)施例的小塊網(wǎng)格化的遞歸方法的流程圖；
[0024]圖10示出根據(jù)本發(fā)明實(shí)施例的小層劃分操作的示意圖；
[0025]圖11示出根據(jù)本發(fā)明實(shí)施例的層位拉平操作的示意圖；
[0026]圖12示圖根據(jù)本發(fā)明實(shí)施例的忽略斷層約束引起的兩種問(wèn)題；
[0027]圖13示出根據(jù)本發(fā)明實(shí)施例的小層與斷層閉合問(wèn)題的解決方法的流程圖；
[0028]圖14示出根據(jù)本發(fā)明實(shí)施例的小層穿越斷層的四種情況；
[0029]圖15示出根據(jù)本發(fā)明實(shí)施例的小層閉合的示意圖；
[0030]圖16示出根據(jù)本發(fā)明實(shí)施例的無(wú)斷層區(qū)域以及有斷層區(qū)域的示意圖；
[0031]圖17示出根據(jù)本發(fā)明實(shí)施例的小塊劃分的方法的流程圖；
[0032]圖18示出根據(jù)本發(fā)明實(shí)施例的斷層范圍統(tǒng)計(jì)的示意圖；
[0033]圖19示出根據(jù)本發(fā)明實(shí)施例的關(guān)鍵斷層與區(qū)域邊界的示意圖；
[0034]圖20A至圖20C分別示出根據(jù)本發(fā)明實(shí)施例的上有斷層、下有斷層、上下均有斷層的小塊構(gòu)造的三種情況；
[0035]圖21示出根據(jù)本發(fā)明實(shí)施例的上下均有斷層的特殊情況；
[0036]圖22示出根據(jù)本發(fā)明實(shí)施例的測(cè)井分層約束的層序模型構(gòu)建。
【具體實(shí)施方式】
[0037]現(xiàn)在，將參照附圖更充分地描述不同的示例實(shí)施例，其中，一些示例性實(shí)施例在附圖中示出。
[0038]圖1示出根據(jù)本發(fā)明實(shí)施例的層序建模的方法的流程圖。
[0039]參照?qǐng)D1，在步驟S110，通過(guò)對(duì)地震解釋的層位、斷層散點(diǎn)數(shù)據(jù)進(jìn)行第一組段，來(lái)找出并記錄每個(gè)層位段相關(guān)聯(lián)的斷層，其中，層位由斷層劃分為層位段，一個(gè)層位由多個(gè)層位段組成。由于數(shù)據(jù)庫(kù)中所保存的層位信息按散點(diǎn)方式存儲(chǔ)，因此將從數(shù)據(jù)庫(kù)中提取的層位數(shù)據(jù)稱為層位散點(diǎn)數(shù)據(jù)。對(duì)于散點(diǎn)數(shù)據(jù)，能夠方便的獲取每個(gè)點(diǎn)的CDP (共深度點(diǎn))坐標(biāo)和時(shí)間(深度)值，然而散點(diǎn)數(shù)據(jù)本身缺少了層位點(diǎn)之間的關(guān)聯(lián)關(guān)系以及層位點(diǎn)與斷層之間的關(guān)聯(lián)關(guān)系等關(guān)鍵信息，因此難以用于后續(xù)的方法處理。通過(guò)組段方法，可以尋找到層位點(diǎn)之間以及層位與斷層之間的關(guān)聯(lián)關(guān)系。
[0040]在步驟S120，通過(guò)根據(jù)相關(guān)聯(lián)的層位、斷層段建立塊狀模型(B卩，成塊)。建立塊狀模型是生成小塊的過(guò)程，是小塊網(wǎng)格化的前提。通過(guò)尋找同時(shí)劃分上下兩個(gè)層位的斷層，并記錄兩兩斷層之間的層位段，從而形成小塊，來(lái)建立塊狀模型。
[0041]在步驟S130，通過(guò)對(duì)在步驟120獲得的小塊進(jìn)行處理來(lái)獲得網(wǎng)格化后的小層結(jié)果。在步驟S140，根據(jù)小層結(jié)果輸出層序建模結(jié)果。下面將參照附圖描述操作SllO至操作S130的詳細(xì)操作。
[0042]圖2示出根據(jù)本發(fā)明實(shí)施例的第一組段的詳細(xì)方法的流程圖。
[0043]組段的實(shí)際工作分為斷層組段以及層位組段兩部分。然而，由于斷層結(jié)構(gòu)簡(jiǎn)單、情況單一，且每個(gè)⑶P坐標(biāo)上均不存在重值點(diǎn)，因此斷層組段方法簡(jiǎn)單，僅需讀取出相鄰的斷層點(diǎn)將其組合即可，因此，此處不再累述。以下將詳細(xì)描述層位數(shù)據(jù)的組段方法。
[0044]首先，在操作S210，確定關(guān)聯(lián)斷層。關(guān)聯(lián)斷層，即將一個(gè)層位段分割為兩個(gè)單獨(dú)層位段的斷層。參照?qǐng)D3，對(duì)于某一斷層，如果當(dāng)前層位存在這樣的點(diǎn)，其CDP坐標(biāo)與某些斷層點(diǎn)的CDP坐標(biāo)相等，而其時(shí)間(深度)值小于對(duì)應(yīng)斷層點(diǎn)的時(shí)間(深度)值，則將其稱為層位相對(duì)于此斷層的上盤點(diǎn)；類似地，如果所述點(diǎn)的時(shí)間(深度)值大于對(duì)應(yīng)斷層點(diǎn)的時(shí)間(深度)值，則將其稱為層位相對(duì)于此斷層的下盤點(diǎn)。根據(jù)圖3可以看出，如果當(dāng)前層位對(duì)于某一斷層既存在上盤點(diǎn)，也存在下盤點(diǎn)，則說(shuō)明當(dāng)前斷層將層位分割，即當(dāng)前斷層為層位的關(guān)聯(lián)斷層。
[0045]因此，針對(duì)每一個(gè)層位點(diǎn)，尋找其上(下)離其最近的斷層點(diǎn)并標(biāo)記。我們可以大體認(rèn)為當(dāng)前層位點(diǎn)是其上(下)斷層的下(上)盤層位點(diǎn)。最后，統(tǒng)計(jì)并找出既存在上盤層位點(diǎn)，也存在下盤層位點(diǎn)的斷層，并將其標(biāo)記為當(dāng)前層位的關(guān)聯(lián)斷層。
[0046]尋找關(guān)聯(lián)斷層的意義在于我們所處理的工區(qū)之中可能存在數(shù)量眾多的斷層數(shù)據(jù)，然而只有關(guān)聯(lián)斷層會(huì)將一個(gè)層位分割為多個(gè)段數(shù)據(jù)。因此，尋找關(guān)聯(lián)斷層能使我們?cè)诮M段時(shí)忽略其他無(wú)效斷層的影響。同時(shí)，由于每個(gè)關(guān)聯(lián)斷層都將層位分割為兩段，因此根據(jù)關(guān)聯(lián)斷層的數(shù)量還可以預(yù)先得到組段之后層位段的數(shù)量。
[0047]在操作S220，計(jì)算層位點(diǎn)(B卩，構(gòu)成層位的點(diǎn))的屬性，即，計(jì)算層位點(diǎn)的上下關(guān)聯(lián)斷層，也就是說(shuō)，確定沿上下方向離層位點(diǎn)最近的兩個(gè)關(guān)聯(lián)的斷層點(diǎn)，即，關(guān)聯(lián)斷層信息。所述點(diǎn)屬性指的是該層位點(diǎn)上下的關(guān)聯(lián)斷層點(diǎn)所屬的斷層，和關(guān)聯(lián)斷層信息相關(guān)。此外，一個(gè)層位點(diǎn)的屬性與其上方和下方的兩個(gè)斷層點(diǎn)信息相關(guān)。通過(guò)遍歷每一個(gè)層位點(diǎn)，找出沿上下方向離該層位點(diǎn)最近的斷層點(diǎn)，并獲取相關(guān)的斷層信息和標(biāo)記，來(lái)確定到特定層位點(diǎn)上下離其最近的關(guān)聯(lián)斷層信息。由于斷層將層位劃分為上盤與下盤兩部分，因此上盤點(diǎn)與下盤點(diǎn)具有不同的屬性值，可以根據(jù)屬性值將上盤點(diǎn)和下盤點(diǎn)分割到不同的段中。
[0048]在操作S230，對(duì)確定了屬性的層位點(diǎn)進(jìn)行第二組段。根據(jù)輸出結(jié)果的不同，可以大體上將組段劃分為組小段與組大段。
[0049]組小段是根據(jù)點(diǎn)屬性的不同，將相同屬性的點(diǎn)進(jìn)行初步聚合的操作。簡(jiǎn)而言之，在所有層位點(diǎn)中遍歷并尋找到具有相同屬性的點(diǎn)，并將其組合為一個(gè)小段。參照?qǐng)D4，根據(jù)屬性值的不同，可以將當(dāng)前層位劃分為7個(gè)小段。
[0050]在某些情況下，僅僅根據(jù)點(diǎn)屬性來(lái)判定小段并非十分準(zhǔn)確。如圖4中示出的第1、
4、7小段，在此三段中，每個(gè)層位的點(diǎn)的上下關(guān)聯(lián)斷層均為空，因此這三段中每個(gè)點(diǎn)的屬性值均相同。然而從圖4中可以看出，這些層位點(diǎn)分屬于三個(gè)不同的小段。因此，在組小段時(shí)，某些情況下還需要考慮層位點(diǎn)之間的位置關(guān)系，通過(guò)點(diǎn)屬性與位置關(guān)系才能準(zhǔn)確的進(jìn)行小段組合。
[0051]從圖4中可以看出，根據(jù)點(diǎn)屬性組合出的7個(gè)小段并非最終的層位段，因此，還需要將組合好的小段合并為最終的大段。采用CDP號(hào)作為判斷大段的標(biāo)準(zhǔn)，也就是說(shuō)，如果某一個(gè)小段的尾點(diǎn)與另一個(gè)小段的首點(diǎn)相鄰，則將其組合為同一個(gè)大段。根據(jù)此方法，可以將圖4中的第1、2段，第3、4、5段，第6、7段分別組合為三個(gè)大段。
[0052]然而，由于目前所使用的層位數(shù)據(jù)情況較為復(fù)雜，同時(shí)層位與斷層之間的閉合計(jì)算可能出現(xiàn)誤差，因此可能出現(xiàn)一種特殊情況導(dǎo)致此判斷方法失效，需要引進(jìn)其他判斷準(zhǔn)則來(lái)進(jìn)行更為準(zhǔn)確的判斷。
[0053]從圖5中可以看出，應(yīng)該將1、2組合為一個(gè)大段，而將3作為一個(gè)單獨(dú)的大段。然而，由于小段3與斷層并未閉合，僅僅依靠CDP信息來(lái)判斷小段的組合關(guān)系可能會(huì)出現(xiàn)錯(cuò)誤，因此，在這種情況下，需要引入其他的條件進(jìn)行輔助判斷。
[0054]具體地說(shuō)，當(dāng)尋找到特定小段時(shí)，遍歷其他小段，根據(jù)特定小段的CDP信息判斷其他小段是否可能與該特定小段進(jìn)行組段。如果未找到能夠與當(dāng)前小段進(jìn)行組段的小段，則當(dāng)前大段組段結(jié)束，并開(kāi)始下一大段的組段工作；如果找到一個(gè)小段能夠與當(dāng)前的特定小段進(jìn)行組段，則將這兩個(gè)小段合并為一個(gè)大段；然而，如果尋找到多個(gè)小段可能與當(dāng)前小段進(jìn)行組段，需要首先讀取出當(dāng)前小段的尾點(diǎn)以及所有符合條件的小段的首點(diǎn)，之后，通過(guò)對(duì)所有小段的時(shí)間值或深度值進(jìn)行比較，找出首點(diǎn)距離與當(dāng)前小段尾點(diǎn)最近的小段，并將其與當(dāng)前小段進(jìn)行組段。
[0055]在完成組段之后，還需要在操作S240確定層位與斷層的關(guān)聯(lián)關(guān)系信息，并根據(jù)此信息對(duì)層位段進(jìn)行排序，以供后續(xù)使用。
[0056]參照?qǐng)D6，在操作S240，針對(duì)每個(gè)層位段，尋找其左右斷點(diǎn)被哪一個(gè)關(guān)聯(lián)斷層所分害I]。首先，判斷當(dāng)前層位段的左右端點(diǎn)是否處于工區(qū)邊界處，如果特定層位段上存在一點(diǎn)位于工區(qū)邊界處，則將其對(duì)應(yīng)關(guān)聯(lián)斷層標(biāo)記為空，如圖6中示出的從左起第一段的左關(guān)聯(lián)斷層及第三段的右關(guān)聯(lián)斷層。
[0057]之后，如果某斷點(diǎn)不在工區(qū)邊界處，則取出其所作⑶P坐標(biāo)上所有的層位與斷層信息，并從中提取出與該層位相關(guān)聯(lián)的斷層的斷層點(diǎn)信息。然后，比較當(dāng)前層位端點(diǎn)與所有關(guān)聯(lián)斷層點(diǎn)之間的距離，并找出與其距離最小的一個(gè)斷層，將其標(biāo)記為該段的左(右)關(guān)聯(lián)斷層。
[0058]然而，由于工區(qū)與數(shù)據(jù)的相互關(guān)系，可能出現(xiàn)如圖7所示的導(dǎo)致關(guān)聯(lián)斷層判斷較為困難的特殊情況。
[0059]由于在實(shí)際的使用過(guò)程中，為了滿足生產(chǎn)需要，可能會(huì)將實(shí)際使用的各工區(qū)數(shù)據(jù)按一定需求進(jìn)行裁剪，導(dǎo)致層位與斷層關(guān)聯(lián)關(guān)系難以判斷。如圖7所示，由于在實(shí)際情況中層位與斷層數(shù)據(jù)被工區(qū)邊界裁剪掉，導(dǎo)致此處的部分關(guān)鍵信息丟失，現(xiàn)有信息量不足，從而導(dǎo)致無(wú)法判斷實(shí)際的層位與斷層之間的關(guān)聯(lián)關(guān)系。因此，將這種情況下的層位與斷層之間的關(guān)系視為以下情形:上下兩個(gè)層位段被中間的斷層分割，并分別與中間的斷層相互關(guān)聯(lián)。
[0060]因此，將在完成基本的組段方法之后查找此種情況并對(duì)組段結(jié)果進(jìn)行修改。具體地說(shuō)，遍歷組段好的所有層位段，尋找位于工區(qū)左(右)邊界的段，判斷此段左(右)關(guān)聯(lián)斷層是否為空。如果存在位于左(右)邊界上的多段層位段，則尋找其兩兩層位段之間的斷層。對(duì)于其中任意兩段相鄰層位段，如果其中存在斷層，則判斷兩層位段的右(左)關(guān)聯(lián)斷層是否為此斷層(即，所述存在的斷層)，如果所述兩層位段中的一段的右(左)關(guān)聯(lián)斷層為此斷層，則將所述兩層位段中的另一層位段的左(右)關(guān)聯(lián)斷層設(shè)置為此斷層；如果沒(méi)有層位右(左)關(guān)聯(lián)此斷層且兩層位之間無(wú)其他斷層，或者兩層位之間無(wú)斷層，則層位信息出錯(cuò)，結(jié)束處理過(guò)程。
[0061]至此，已完成層位由散點(diǎn)到段結(jié)構(gòu)的組段過(guò)程。由于后續(xù)需要，需要將層位段從左至右進(jìn)行排序。具體地說(shuō)，首先找到左關(guān)聯(lián)斷層為空的層位段，之后循環(huán)查找左關(guān)聯(lián)斷層與上一層位段的右關(guān)聯(lián)斷層相同的層位段，一直到找到一個(gè)右關(guān)聯(lián)斷層為空的層位段為止。
[0062]通過(guò)圖2示出的組段方法，可以將特定層位上的離散的點(diǎn)信息組合成為由斷層分割的段信息，并得到段信息與斷層之間的關(guān)聯(lián)關(guān)系。
[0063]圖8示出根據(jù)本發(fā)明實(shí)施例的成塊原理的示意圖。
[0064]在層序建模中，實(shí)際所需的最終結(jié)構(gòu)是通過(guò)層位與斷層約束所模擬出的兩個(gè)層位之間的空間位置中的構(gòu)造結(jié)構(gòu)。參照?qǐng)D8，當(dāng)兩個(gè)層位之間的特定區(qū)域左右均被斷層或工區(qū)邊界從整體區(qū)域中整體劃分出來(lái)，則可認(rèn)為其內(nèi)部構(gòu)造結(jié)構(gòu)不受其他區(qū)域構(gòu)造結(jié)構(gòu)的影響，因此可以作為層序建模中的一個(gè)單獨(dú)的構(gòu)造區(qū)域，所以，在進(jìn)行實(shí)際的層序建模之前，需要首先對(duì)兩個(gè)層位之間的區(qū)域進(jìn)行分塊。
[0065]參照?qǐng)D8，將同時(shí)與上下兩個(gè)層位相關(guān)聯(lián)的斷層稱為大斷層，而將只與其中之一關(guān)聯(lián)的斷層稱為小斷層。大斷層與工區(qū)邊界構(gòu)成了用于劃分塊區(qū)域的標(biāo)準(zhǔn)。因此，通過(guò)找出兩個(gè)層位之間的所有大斷層，并將其排序，來(lái)根據(jù)排好序的大斷層以及工區(qū)邊界信息將組好的層位段劃分成小塊，所述小塊指的是由層位與斷層為邊界構(gòu)成的封閉的全連通區(qū)域。
[0066]成塊的過(guò)程對(duì)后續(xù)所需要處理的區(qū)域進(jìn)行了進(jìn)一步劃分，構(gòu)成了后續(xù)處理所使用的基本單位，所述基本單位指的是使得兩兩小塊之間的建模過(guò)程互不影響的基本單元。因此，只需要考慮單獨(dú)考慮塊內(nèi)的信息，而不會(huì)受到其他部分的影響。
[0067]小塊網(wǎng)格化，又稱為小層劃分，是層序建模的基本步驟。小塊網(wǎng)格化根據(jù)當(dāng)前塊中的層位信息和斷層信息，構(gòu)造出經(jīng)過(guò)所述層位信息和斷層信息約束之后的扭曲的空間網(wǎng)格。圖10示出根據(jù)本發(fā)明實(shí)施例的小層劃分的示意圖。圖10中的橫向網(wǎng)格線稱為小層，即它是根據(jù)上下層位信息在當(dāng)前位置模擬的出的一些虛擬的層位信息。
[0068]小塊網(wǎng)格化的過(guò)程根據(jù)小塊的結(jié)構(gòu)可以分為兩種情況，即無(wú)內(nèi)部斷層與有內(nèi)部斷層的情況。如果小塊無(wú)內(nèi)部斷層(如圖8中的塊I與塊3)，則由于小塊為一個(gè)簡(jiǎn)單的多邊形模型，不需要進(jìn)一步劃分，因此可以直接進(jìn)行網(wǎng)格化操作，計(jì)算當(dāng)前小塊內(nèi)部的層序模型，并輸出層序建模結(jié)果，該層序建模結(jié)果如圖10中示出的橫向的網(wǎng)格線。然而，如果小塊中內(nèi)部存在斷層(如圖8中的塊2)，則無(wú)法在此小塊的基礎(chǔ)上直接計(jì)算出層序模型，而需要對(duì)其進(jìn)行一些預(yù)處理，并將其交給小塊劃分模塊進(jìn)行進(jìn)一步的小塊細(xì)分，最后使用遞歸的方法處理完所有小塊。
[0069]圖9示出根據(jù)本發(fā)明實(shí)施例的小塊網(wǎng)格化的遞歸方法的流程圖。
[0070]在步驟S910，從存儲(chǔ)小塊信息的棧頂取出一個(gè)小塊，對(duì)其進(jìn)行相應(yīng)處理。在步驟S920，判斷小塊內(nèi)部是否存在斷層。如果小塊內(nèi)部無(wú)斷層，則在步驟S930，對(duì)小塊進(jìn)行網(wǎng)格化。具體地說(shuō)，首先，根據(jù)當(dāng)前小塊的上下層位段信息，可以得到其上下層位段的起始及結(jié)束位置。根據(jù)此信息，可以將上下層位段按等比例劃分為相同段數(shù)的小段。之后，分別連接上下層位段上相同比例處的分割點(diǎn)，如圖10所示，可以得到圖中豎直方向的直線部分。[0071]這些直線就是后續(xù)劃分小層的依據(jù)，被稱為小層劃分的參考線。根據(jù)預(yù)先設(shè)定好的參數(shù)，可以獲得用戶所要求的小層個(gè)數(shù)。之后，根據(jù)小層個(gè)數(shù)，將每條參考線進(jìn)行相應(yīng)的等比例劃分，最后，連接相應(yīng)的參考線上的分割點(diǎn)，即得到所需的小層。
[0072]參照?qǐng)D10，在網(wǎng)格化分之后，可以得到如圖10中所示的所有交點(diǎn)的坐標(biāo)信息。之后，在步驟S940，可以根據(jù)網(wǎng)格信息(B卩，網(wǎng)格中每個(gè)頂點(diǎn)的實(shí)際空間坐標(biāo))輸出小層信息(即，關(guān)于兩個(gè)層位之間估計(jì)所得的其他層位的信息)，或者直接將網(wǎng)格信息輸出用于后續(xù)的確定建模處理。由于網(wǎng)格點(diǎn)坐標(biāo)為浮點(diǎn)數(shù)坐標(biāo)，因此在輸出小層信息時(shí)，需要將小層點(diǎn)投影至整數(shù)坐標(biāo)上，再進(jìn)行小層輸出。
[0073]通過(guò)以上描述可知，小塊由斷層和層位段構(gòu)成。然而，由于斷層原本為曲線，而在上述處理中將其視為直線進(jìn)行處理。因此，小塊網(wǎng)格化中作為基準(zhǔn)的斷層是一條擬合的直線，這導(dǎo)致網(wǎng)格化結(jié)果不準(zhǔn)確，從而產(chǎn)生閉合問(wèn)題。需要通過(guò)小層閉合操作對(duì)該問(wèn)題進(jìn)行修正。因此，在實(shí)際的小層輸出之前，仍然需要對(duì)小層進(jìn)行優(yōu)化處理，該優(yōu)化處理將在之后進(jìn)行詳細(xì)介紹。
[0074]如果在步驟S920確定小塊內(nèi)部存在斷層，則由于其內(nèi)部構(gòu)造不僅僅受上下層位的約束，還要受到內(nèi)部斷層的約束，因此不能直接使用上下層位信息進(jìn)行直接的網(wǎng)格劃分。因此，為了保證小層劃分的準(zhǔn)確性，需要對(duì)其進(jìn)行一定的處理，并交由小塊劃分模塊進(jìn)行進(jìn)一步的小塊細(xì)分，將其劃分為多個(gè)內(nèi)部無(wú)斷層小塊之后，再分別進(jìn)行處理。
[0075]為了給小塊劃分模塊提供一定的劃分依據(jù)，需要對(duì)當(dāng)前小塊進(jìn)行預(yù)處理，并進(jìn)行網(wǎng)格化。此時(shí)，網(wǎng)格化的結(jié)果并不作為實(shí)際的小層劃分依據(jù)，它與最終的小層信息沒(méi)有實(shí)際的關(guān)聯(lián)。
[0076]由于小塊內(nèi)部存在斷層，其上下層位段可能包含多段，因此在對(duì)其進(jìn)行網(wǎng)格化時(shí)，難以將其層位段進(jìn)行等比例劃分。因此，在進(jìn)行步驟S970的小塊網(wǎng)格化之前，需要在步驟S960，對(duì)包括多于一個(gè)層位段的層位進(jìn)行層位拉平的處理，將其近似為一個(gè)整體的層位段，將此步操作稱為層位拉平。
[0077]層位拉平即將兩段或者多段由斷層分割的層位段通過(guò)近似的方法還原為一個(gè)層位段。圖11是示出根據(jù)本發(fā)明示例性實(shí)施例的層位拉平的操作的示意圖。首先，由于單一層位段之中每個(gè)CDP坐標(biāo)點(diǎn)上均不存在重值點(diǎn)，因此需要先刪除其所有重值點(diǎn)信息。之后，統(tǒng)計(jì)出剩余層位部分在空白區(qū)域處的端點(diǎn)位置，并根據(jù)線性插值的方法將空白區(qū)域補(bǔ)齊。最后，為了保證最終小塊劃分的準(zhǔn)確性及美觀性，需要使用線平滑的方式將層位進(jìn)行平滑。
[0078]線平滑主要采用加權(quán)平均的方式對(duì)層位段進(jìn)行平滑。其權(quán)重分配使用1:2:1的分配方式，即當(dāng)前點(diǎn)標(biāo)量值=(前一點(diǎn)標(biāo)量值+后一點(diǎn)標(biāo)量值+當(dāng)前點(diǎn)標(biāo)量值*2)/4。通過(guò)此種方式可以逐漸的將層位進(jìn)行平滑處理，在實(shí)際應(yīng)用中，所需要的平滑次數(shù)約為500次。
[0079]通過(guò)抹去內(nèi)部斷層及層位拉平等操作，可以將一個(gè)內(nèi)部有斷層的小塊近似的轉(zhuǎn)化為一個(gè)內(nèi)部無(wú)斷層的小塊。之后，可以使用內(nèi)部無(wú)斷層小塊網(wǎng)格化方法對(duì)其進(jìn)行網(wǎng)格化操作。然而，由于忽略了內(nèi)部斷層的影響，此處的網(wǎng)格化操作結(jié)果極不準(zhǔn)確，無(wú)法直接用于小層結(jié)果的計(jì)算，只能用作在步驟S980的小塊劃分操作的數(shù)據(jù)基礎(chǔ)。將在下面參照?qǐng)D17詳細(xì)描述小塊劃分的方法的流程圖。在步驟S990，在完成小塊劃分處理之后，將處理后的小塊信息壓棧。之后處理返回步驟S910，直到在步驟S950確定棧為空時(shí)，結(jié)束處理。這里，由于層序建模最終的輸出結(jié)果為小塊網(wǎng)格，具有內(nèi)部斷層的小塊雖然進(jìn)行了層位拉平處理了消除了內(nèi)部斷層，但其建立的小塊網(wǎng)格并不準(zhǔn)確，因此只能作為小塊劃分的標(biāo)準(zhǔn)。為了獲得準(zhǔn)確的結(jié)果，需要將一個(gè)內(nèi)部含有斷層的小塊劃分為多個(gè)內(nèi)部不含斷層的小塊，之后再對(duì)這些小塊進(jìn)行準(zhǔn)確的小塊網(wǎng)格化。
[0080]在以上所介紹的網(wǎng)格化方法中，為了保證方法的高效性與簡(jiǎn)便性，采用上下層位的端點(diǎn)連線來(lái)模擬邊界斷層的方法。這種方法實(shí)際上忽略了邊界斷層對(duì)小塊網(wǎng)格化過(guò)程的影響，雖然簡(jiǎn)化了方法，卻會(huì)導(dǎo)致包括小層穿越及小層不足在內(nèi)的一些問(wèn)題，如圖12所示。為此，在進(jìn)行最終的小層結(jié)果輸出之前，需要先對(duì)小層進(jìn)行優(yōu)化，使其滿足小塊的邊界斷層約束，使網(wǎng)格化的結(jié)果更加的準(zhǔn)確，其中，各個(gè)單個(gè)小塊的小層輸出結(jié)果可以組成整體層序建豐吳結(jié)果。
[0081]小層與斷層閉合的問(wèn)題主要包括如圖12所示的兩種。為了解決這兩種問(wèn)題，本方法所使用的方法是根據(jù)邊界斷層對(duì)小層進(jìn)行優(yōu)化，具體方法為先將小層穿越斷層的部分進(jìn)行裁剪，再將小層不足的部分進(jìn)行補(bǔ)齊，使其最終與斷層閉合，將在以下參照?qǐng)D13對(duì)該操作進(jìn)行詳細(xì)描述。
[0082]在步驟S1310，進(jìn)行小層裁剪。小層裁剪的主要目的是將小層穿越斷層的部分進(jìn)行裁剪，其主要難點(diǎn)在于判斷小層是否穿越斷層，并找出所有穿越斷層的小層點(diǎn)。
[0083]判斷小層點(diǎn)是否穿越斷層點(diǎn)可以通過(guò)小層點(diǎn)與對(duì)應(yīng)⑶P坐標(biāo)上的斷層點(diǎn)的時(shí)間(深度)值的比較來(lái)確定。圖14示出根據(jù)本發(fā)明示例性實(shí)施例的小層穿越的四種情況。首先，根據(jù)斷層起始點(diǎn)和結(jié)束點(diǎn)的時(shí)間(深度)值以及其為小塊的左或右邊界等信息，可以計(jì)算出其參考線的傾斜方向。之后，根據(jù)小層點(diǎn)的時(shí)間(深度)值與斷層點(diǎn)的位置關(guān)系，以及斷層的傾斜方向，可以判斷出當(dāng)前點(diǎn)是否屬于小層穿越斷層部分的點(diǎn)。例如，如果特定斷層為某小塊的右邊界斷層，且其起始點(diǎn)的時(shí)間(深度)值小于結(jié)束點(diǎn)的時(shí)間(深度)值，則可以判斷出該斷層的參考線方向如圖14的左上圖所示。因此，凡是時(shí)間(深度)值小于同CDP坐標(biāo)下斷層點(diǎn)的小層點(diǎn)，均屬于穿越斷層的小層點(diǎn)。因此，可以從右至左遍歷整個(gè)小層，并將穿越斷層的小層點(diǎn)刪除。由于經(jīng)過(guò)小塊網(wǎng)格化所計(jì)算出的小層點(diǎn)坐標(biāo)為浮點(diǎn)數(shù)坐標(biāo)，因此在實(shí)際判斷時(shí)需要對(duì)斷層進(jìn)行線性插值以計(jì)算相應(yīng)坐標(biāo)點(diǎn)上的時(shí)間(深度)值。
[0084]在完成了小層裁剪之后，可以保證所有的小層均不穿越斷層，因此，可以在步驟S1320對(duì)進(jìn)行小層補(bǔ)齊操作，即，對(duì)所有的小層進(jìn)行小層與斷層閉合的操作而不用判斷小層是否穿越斷層。
[0085]圖15示出根據(jù)本發(fā)明示例性實(shí)施例的小層閉合的示意圖。參照?qǐng)D15，小層閉合采用了相同時(shí)間(深度)值延長(zhǎng)的方法。即針對(duì)每一個(gè)待處理的小層，取出其端點(diǎn)位置的時(shí)間(深度)值，并按此時(shí)間(深度)值對(duì)小層做延長(zhǎng)線，計(jì)算出此延長(zhǎng)線與小塊邊界斷層的交點(diǎn)。之后，根據(jù)斷層交點(diǎn)以及小層端點(diǎn)的距離，判斷是否需要進(jìn)行小層閉合操作?？赏ㄟ^(guò)使用當(dāng)前小層的端點(diǎn)和斷層點(diǎn)之間距離的絕對(duì)值來(lái)判斷小層是否穿越斷層。如果小層需要閉合操作，則從小層端點(diǎn)處出發(fā)，按一定的步長(zhǎng)在小層延長(zhǎng)線上進(jìn)行采樣，并將采樣點(diǎn)添加到小層之中。在實(shí)際使用中，為了保證后續(xù)確定建模等方法效果，優(yōu)選采樣步長(zhǎng)小于1/2個(gè)CDP坐標(biāo)標(biāo)準(zhǔn)單位長(zhǎng)度。采樣一直到當(dāng)前采樣點(diǎn)穿越斷層為止，至此，可保證小層與其關(guān)聯(lián)斷層完全符合。
[0086]如果層序建模的結(jié)果用作隨機(jī)建模或?qū)傩越５妮斎?，則由于方法的特殊要求，需要保證每個(gè)小層上的點(diǎn)數(shù)相等。因此，在完成小層閉合操作之后，需要在步驟S1330，統(tǒng)計(jì)出當(dāng)前小塊中小層的最大點(diǎn)數(shù)，并在其他小層首尾添加重值點(diǎn)，使所有小層長(zhǎng)度一致。之后，在步驟S1340，輸出小層結(jié)果。
[0087]圖17示出根據(jù)本發(fā)明示例性實(shí)施例的小塊劃分的方法的流程圖，包括斷層范圍計(jì)算(操作S1710)、關(guān)鍵斷層選取(操作S1720)、約束小層構(gòu)建(操作S1730)和小塊構(gòu)建(S1740)四個(gè)步驟。
[0088]小塊劃分，即將內(nèi)部存在斷層的小塊劃分為多個(gè)更小的小塊。由于內(nèi)部斷層會(huì)影響小塊內(nèi)部層序構(gòu)造，使小塊難以構(gòu)建正確的層序模型，因此需要盡可能的消除內(nèi)部斷層，從而將其轉(zhuǎn)化為更容易處理的邊界斷層。小塊劃分通過(guò)將小塊劃分為更小的塊來(lái)將其內(nèi)部斷層轉(zhuǎn)換為劃分后部分小塊的邊界斷層，從而減少了內(nèi)部斷層。再加上使用了遞歸調(diào)用的方法，因此小塊劃分最終能夠?qū)⑺械男K劃分成為一系列無(wú)內(nèi)部斷層的小塊，為后續(xù)的層序建模工作提供保證。
[0089]小塊劃分的主要依據(jù)是進(jìn)行層位拉平之后所建立的初步的網(wǎng)格模型。雖然經(jīng)過(guò)層位拉平之后忽略了內(nèi)部斷層對(duì)小塊內(nèi)部結(jié)構(gòu)的影響，無(wú)法準(zhǔn)確的描述小塊內(nèi)部構(gòu)造，但由于小塊內(nèi)部的主要構(gòu)造按層劃分，所以對(duì)于無(wú)斷層約束的部分，此種網(wǎng)格化方法仍具有一定的可靠性，因此，可以根據(jù)此網(wǎng)格化所得到的近似的小層結(jié)果將原小塊劃分為有斷層約束區(qū)域及無(wú)斷層約束區(qū)域兩部分，如圖16所示，以方便進(jìn)一步的小塊劃分工作。
[0090]根據(jù)這一原理，小塊劃分可以將具有內(nèi)部斷層的小塊進(jìn)行進(jìn)一步的細(xì)分，最終得到多個(gè)無(wú)內(nèi)部斷層的小塊。
[0091]參照?qǐng)D17，在操作S1710，由于需要使用根據(jù)層位拉平之后的結(jié)果所建立的網(wǎng)格模型(即，近似網(wǎng)格模型)進(jìn)行判斷，并在其中找出有斷層區(qū)域及無(wú)斷層區(qū)域的分割線，因此需要首先統(tǒng)計(jì)出每個(gè)內(nèi)部斷層所覆蓋的小層范圍。該步驟分為上盤斷層及下盤斷層兩種情況，其中，上盤斷層為只關(guān)聯(lián)上層位段的斷層，相應(yīng)地，下盤斷層為只關(guān)聯(lián)下層位段的斷層。由于兩種情況類似，因此此處將以下盤斷層為例，介紹相關(guān)方法。
[0092]圖18示出根據(jù)本發(fā)明實(shí)施例的斷層范圍統(tǒng)計(jì)的示意圖。參照?qǐng)D18，由于下盤斷層在小塊內(nèi)部的部分從下層位處開(kāi)始，因此，通過(guò)統(tǒng)計(jì)下盤斷層穿過(guò)的近似網(wǎng)格模型中的小層，來(lái)尋找其上端點(diǎn)所在位置的小層號(hào)。然而，由于小層線是按大量的離散點(diǎn)連接而成的折線，因此直接通過(guò)小層線來(lái)計(jì)算端點(diǎn)位置較為困難。因此，首先尋找離斷層端點(diǎn)最近的一根參考線。由于參考線實(shí)質(zhì)上是上層位段與下層位段上兩點(diǎn)的連線，因此可以很方便的計(jì)算出該參考線的直線表示式。之后，將斷層端點(diǎn)的時(shí)間(深度)值代入此表達(dá)式，可得出該參考線與此平面相交的CDP坐標(biāo)。之后，可通過(guò)尋找到在CDP坐標(biāo)上與該斷層端點(diǎn)最接近的一個(gè)交點(diǎn)，來(lái)獲得離該斷層端點(diǎn)最近的一條參考線。而小層點(diǎn)在參考線上按等比例排布，因此可通過(guò)計(jì)算出交點(diǎn)在參考線上的比例位置，來(lái)計(jì)算出小塊內(nèi)部的斷層端點(diǎn)所在的小層位置。
[0093]由于小塊內(nèi)部的實(shí)際構(gòu)造根據(jù)斷層和層位情況發(fā)生一定的扭曲，因此，此處所計(jì)算出的斷層影響范圍并非以斷層長(zhǎng)度為準(zhǔn)，而是考慮當(dāng)小塊內(nèi)部區(qū)域被層序化之后，該斷層影響的小層數(shù)。
[0094]在計(jì)算出每個(gè)斷層的覆蓋范圍之后，在步驟S1720，從所述覆蓋范圍中選取實(shí)際對(duì)區(qū)域劃分產(chǎn)生影響的斷層。由于影響范圍較小的斷層的影響區(qū)域會(huì)被影響范圍較大的斷層的影響區(qū)域所覆蓋，導(dǎo)致其無(wú)法約束區(qū)域邊界。因此，如圖19所示，在所有的內(nèi)部斷層中，通過(guò)從上盤斷層和下盤斷層中分別找出影響范圍最大的斷層，來(lái)確定實(shí)際的區(qū)域邊界。
[0095]在確定了離斷層端點(diǎn)最近(S卩，劃分區(qū)域邊界所使用)的小層后，通過(guò)前述小層輸出處理，將該小層輸出為一個(gè)新的層位段。由于對(duì)于劃分后得到的小塊，其上層為或者下層為可能不是原有的層位，而是通過(guò)前述小層輸出處理而得到的一個(gè)新的層位，因此所述小層可能作為新小塊的上層位或者下層位，從而需要對(duì)該小層的數(shù)據(jù)格式進(jìn)行一定的修正。
[0096]小層數(shù)據(jù)與層位數(shù)據(jù)最主要的區(qū)別在于其坐標(biāo)位置。小層數(shù)據(jù)的坐標(biāo)使用了浮點(diǎn)數(shù)坐標(biāo)，以便提升計(jì)算的精度。而層位坐標(biāo)采用的是以CDP坐標(biāo)為準(zhǔn)的整數(shù)坐標(biāo)，以方便數(shù)據(jù)的存儲(chǔ)與使用。為了保證方法的統(tǒng)一，需要將所有的輸入小塊的層位和斷層信息進(jìn)行統(tǒng)一。因此，在進(jìn)行后續(xù)的小塊構(gòu)建之前，需要先將其中作為約束線(即用于小塊劃分的約束小層)的小層數(shù)據(jù)重新采樣成層位格式。小層數(shù)據(jù)的重新采樣所使用的方法較為簡(jiǎn)單，即根據(jù)相鄰兩點(diǎn)在對(duì)應(yīng)的整數(shù)坐標(biāo)位置進(jìn)行線性插值。
[0097]在完成上述操作之后，在步驟S1740，根據(jù)層位、斷層、塊內(nèi)斷層端點(diǎn)以及約束小層信息，將當(dāng)前塊進(jìn)一步劃分為多個(gè)小塊。其中，關(guān)鍵斷層將成為部分小塊的邊界斷層，而約束線將成為部分小塊的邊界層位段。至此，當(dāng)前塊中的關(guān)鍵斷層將不再出現(xiàn)在任何小塊的內(nèi)部，從而不會(huì)再影響小塊網(wǎng)格化的過(guò)程。通過(guò)如圖9所述的遞歸調(diào)用的方法，可以逐漸的消除所有的內(nèi)部斷層的影響，最終將所有分塊劃分成無(wú)內(nèi)部斷層的簡(jiǎn)單小塊，提高層序建模的準(zhǔn)確性。
[0098]圖20A至圖20C分別示出上有斷層、下有斷層、上下均有斷層的小塊構(gòu)造的三種情況。參照?qǐng)D20A至圖20B，將小塊劃分為三塊，以上有斷層的情況為例，將小塊約束線以上的部分劃分為兩塊，分別以關(guān)鍵斷層為其左邊界斷層和右邊界斷層。而約束線以下自成一塊，以約束線所構(gòu)造的層位為其上層位。
[0099]參照?qǐng)D20C，對(duì)于上下均有斷層的情況，由于存在兩天約束線，則需將小層劃分為五塊。即上部?jī)蓧K，下部?jī)蓧K以及中部一塊。之后，在依據(jù)對(duì)應(yīng)的方式計(jì)算其上下層位及左右關(guān)聯(lián)斷層即可。在以上情況中，上下關(guān)鍵斷層可能存在多個(gè)，然而其處理方法與其對(duì)應(yīng)的基本處理方法類似，此處不再累述。
[0100]而對(duì)于上下均有斷層的情況，存在一種如圖21所示的特殊情況。在這種情況下，上下關(guān)鍵斷層覆蓋區(qū)域相互重疊，導(dǎo)致上下分塊區(qū)域相互交疊，且中間區(qū)域消失。在這種情況下，如果仍然按基本方法處理，則會(huì)導(dǎo)致交疊區(qū)域的小塊出現(xiàn)重值點(diǎn)，導(dǎo)致層序建模結(jié)果錯(cuò)誤。為此，目前所采用的方法是估計(jì)虛擬約束線的方法。即根據(jù)上下兩條約束線重新估算出一條中間約束線，并根據(jù)此中間約束線將原始分塊劃分為4塊，如圖21所示。
[0101]圖22示出根據(jù)本發(fā)明實(shí)施例的測(cè)井分層約束的層序模型構(gòu)建。在發(fā)明中，首先在不考慮斷層影響的情況下進(jìn)行層序模型及屬性模型的構(gòu)建。參照?qǐng)D22，由于每口井上每個(gè)TOP的厚度信息已知，因此可以使用克里金插值方法查找出每一道上每個(gè)TOP點(diǎn)的厚度，再根據(jù)等比例壓縮的方法將其壓縮到一個(gè)可接受的范圍內(nèi)。如果插值所得的厚度為O或者負(fù)值，則說(shuō)明該小層在此處尖滅，則將其厚度設(shè)為O。
[0102]根據(jù)本發(fā)明實(shí)施例的層序建模方法是一種基于二維剖面的層序建模方法，即將整個(gè)工區(qū)按CDP坐標(biāo)劃分為二維剖面，在每個(gè)剖面上單獨(dú)進(jìn)行層序建模方法。將每個(gè)剖面的層序建模結(jié)果組合成三維空間中的建模結(jié)果。
[0103]為了解決含逆斷層的復(fù)雜地質(zhì)構(gòu)造在層序建模過(guò)程中的小層劃分不清、小斷層層序走向不確定和網(wǎng)格化不清楚、尖滅點(diǎn)、上超點(diǎn)無(wú)法確定的問(wèn)題，本發(fā)明提出了采用測(cè)井地質(zhì)層位(分層)和斷層雙重約束的方法，建立準(zhǔn)確、合理的層序模型，提高了層序建模的效率，節(jié)省了油藏建模工作者的工作時(shí)間，并且得到的結(jié)果更加合理。
[0104]通過(guò)本發(fā)明的實(shí)施例，對(duì)于含逆斷層的層序模型的建立，能夠節(jié)省大量的工作時(shí)間，得到的層序模型邊界與實(shí)際地震解釋結(jié)果相一致，不需要任何簡(jiǎn)化，層序模型的尖滅、上超等地質(zhì)情況與實(shí)際情況更加吻合，解決了目前油藏建模大量運(yùn)用簡(jiǎn)化手段、得到的模型結(jié)果與實(shí)際地質(zhì)分布規(guī)律脫節(jié)的問(wèn)題。
[0105]盡管已經(jīng)參照其示例性實(shí)施例具體顯示和描述了本發(fā)明，但是本領(lǐng)域的技術(shù)人員應(yīng)該理解，在不脫離權(quán)利要求所限定的本發(fā)明的精神和范圍的情況下，可以對(duì)其進(jìn)行形式和細(xì)節(jié)上的各種改變。
【權(quán)利要求】
1.一種對(duì)復(fù)雜地質(zhì)構(gòu)造層序建模的方法，該方法包括: 步驟A:通過(guò)對(duì)地震解釋的層位散點(diǎn)數(shù)據(jù)、斷層散點(diǎn)數(shù)據(jù)進(jìn)行第一組段，來(lái)找出并記錄各個(gè)層位的關(guān)聯(lián)斷層，其中，所述關(guān)聯(lián)斷層是將所述各個(gè)層位分割為兩個(gè)單獨(dú)層位段的斷層，并且所述第一組段操作包括對(duì)斷層散點(diǎn)數(shù)據(jù)進(jìn)行第一組段的操作和對(duì)層位散點(diǎn)數(shù)據(jù)進(jìn)行第一組段的操作； 步驟B:根據(jù)彼此關(guān)聯(lián)的層位和斷層建立塊狀模型，其中，通過(guò)找出兩個(gè)層位之間的所有大斷層并將其排序，來(lái)根據(jù)排好序的大斷層以及工區(qū)邊界信息將組段后的層位段劃分成小塊，所述大斷層是與上下兩個(gè)層位相關(guān)聯(lián)的斷層； 步驟C:通過(guò)對(duì)在步驟B劃分的小塊進(jìn)行處理來(lái)獲得網(wǎng)格化后的小層結(jié)果； 步驟D:根據(jù)在步驟C獲得的小層結(jié)果輸出整體層序建模結(jié)果。
2.如權(quán)利要求1所述的方法，其中，步驟C包括: 步驟1:判斷小塊內(nèi)部是否包括斷層； 步驟2:如果小塊內(nèi)部不包括斷層，則對(duì)小塊進(jìn)行網(wǎng)格化； 步驟3:如果小塊內(nèi)部包括斷層，則通過(guò)抹去內(nèi)部斷層操作及層位拉平操作，將內(nèi)部有斷層的小塊近似轉(zhuǎn)化為內(nèi)部無(wú)斷層的小塊；對(duì)小塊進(jìn)行網(wǎng)格化；對(duì)網(wǎng)格化后的小塊進(jìn)行劃分；返回步驟I，直到所有小塊中都不包括斷層； 步驟4:輸出網(wǎng)格化后的 小層結(jié)果。
3.如權(quán)利要求2所述的方法，其中，對(duì)小塊進(jìn)行網(wǎng)格化的步驟包括: 根據(jù)小塊的上下層位段信息，得到小塊的上下層位段的起始及結(jié)束位置； 根據(jù)所述起始及結(jié)束位置，將上下層位段按等比例劃分為相同段數(shù)的小段； 分別連接上下層位段上相同比例處的分割點(diǎn)，連接上下層位段上相同比例處的分割點(diǎn)的直線是參考線； 根據(jù)小層個(gè)數(shù)，將每條參考線進(jìn)行相應(yīng)的等比例劃分； 通過(guò)連接相應(yīng)的參考線上的分割點(diǎn)，來(lái)獲得所需的小層。
4.如權(quán)利要求2所述的方法，其中，對(duì)小塊進(jìn)行劃分的步驟包括:斷層范圍計(jì)算、關(guān)鍵斷層選取、約束小層構(gòu)建和小塊構(gòu)建。
5.如權(quán)利要求1所述的方法，其中，對(duì)層位數(shù)據(jù)進(jìn)行第一組段的步驟包括: 確定與各個(gè)層位相關(guān)聯(lián)的斷層； 遍歷每個(gè)層位點(diǎn)，找出沿上下方向離該層位點(diǎn)最近的斷層點(diǎn)，獲取關(guān)聯(lián)斷層信息并進(jìn)行標(biāo)記； 對(duì)確定了關(guān)聯(lián)斷層信息的層位點(diǎn)進(jìn)行第二組段，所述第二組段操作包括組小段和組大段； 確定層位與斷層的關(guān)聯(lián)關(guān)系信息，并根據(jù)所述關(guān)聯(lián)關(guān)系信息對(duì)層位段進(jìn)行排序； 其中， 組小段操作包括:在所有層位點(diǎn)中遍歷并尋找到具有相同屬性的點(diǎn)；并將相同屬性的點(diǎn)組合為一個(gè)小段； 組大段操作包括:如果小段的尾點(diǎn)與另一小段的首點(diǎn)相鄰，則將所述兩個(gè)小段組合為一個(gè)大段。
6.如權(quán)利要求5所述的方法，其中，所述第二組段操作還包括:當(dāng)尋找到特定小段時(shí)，遍歷其他小段，根據(jù)該特定小段的共深度點(diǎn)CDP信息判斷其他小段是否可能與該特定小段進(jìn)行組段； 如果未找到能夠與當(dāng)前小段進(jìn)行組段的小段，則當(dāng)前大段組段結(jié)束，并開(kāi)始下一大段的組段工作； 如果找到一個(gè)小段能夠與當(dāng)前的特定小段進(jìn)行組段，則將這兩個(gè)小段合并為一個(gè)大段； 如果找到能與當(dāng)前小段進(jìn)行組段的多個(gè)小段，則首先讀取出當(dāng)前小段的尾點(diǎn)以及所有符合條件的小段的首點(diǎn)；通過(guò)對(duì)所有小段的時(shí)間值或深度值進(jìn)行比較，找出首點(diǎn)距離與當(dāng)前小段尾點(diǎn)最近的小段，并將該小段與當(dāng)前的特定小段進(jìn)行組段。
7.如權(quán)利要求5所述的方法，其中，根據(jù)關(guān)聯(lián)關(guān)系信息對(duì)層位段進(jìn)行排序的步驟包括:根據(jù)關(guān)聯(lián)關(guān)系信息對(duì)層位段從左到右進(jìn)行排序，其中，首先找到左關(guān)聯(lián)斷層為空的層位段，之后循環(huán)插值左關(guān)聯(lián)斷層與上一層位段的右關(guān)聯(lián)斷層相同的層位段，直到找到一個(gè)右關(guān)聯(lián)斷層為空的層位段為止。
8.一種對(duì)復(fù)雜地質(zhì)構(gòu)造三維層序建模的方法，包括: 將整個(gè)工區(qū)按共深度點(diǎn)CDP坐標(biāo)劃分為二維剖面； 在每個(gè)剖面上單獨(dú) 進(jìn)行如權(quán)利要求1所述的層序建模的方法； 將每個(gè)剖面的層序建模結(jié)果組合成三維空間中的建模結(jié)果。
【文檔編號(hào)】G06T17/05GK103632396SQ201310598209
【公開(kāi)日】2014年3月12日 申請(qǐng)日期:2013年11月22日 優(yōu)先權(quán)日:2013年11月22日
【發(fā)明者】王頎, 洪余剛, 魯才, 陶正喜, 鄒文, 張洞君, 黃東山, 周晶晶, 錢津 申請(qǐng)人:中國(guó)石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司