專利名稱：：地震層位骨架化的制作方法
技術(shù)領(lǐng)域：
：本發(fā)明一般涉及地球物理和地質(zhì)勘探領(lǐng)域，更具體地涉及地震數(shù)據(jù)的分析。具體地，本發(fā)明是為了突出可能包含烴的區(qū)域產(chǎn)生目標(biāo)諸如表面和地質(zhì)體(geobody)并且自動(dòng)分析它們的方法。本發(fā)明一種具體的實(shí)施方式是從地震數(shù)據(jù)體(seismicdatavolume)中同時(shí)產(chǎn)生并分析許多地層學(xué)一致的表面。
背景技術(shù)：
：在地震數(shù)據(jù)處理和解釋中將地震數(shù)據(jù)體減少至其基于內(nèi)反射的表面或?qū)游皇怯欣摹＿@些表面共同形成地震體(seismicvolume)的骨架(skeleton)。已經(jīng)描述了許多通過地震數(shù)據(jù)體一次地抽取或跟蹤一個(gè)層位或表面的方法。這些方法中的大多數(shù)都產(chǎn)生最終相互重疊的表面。因此，同一表面可具有與同一空間位置相關(guān)的多個(gè)深度(或反射時(shí)間)。一些方法通過放棄每個(gè)位置除了一個(gè)以外的所有值防止多值表面。一般地，它們只存儲(chǔ)在執(zhí)行該方法期間遇到的第一個(gè)值，只是不記錄后面的值。而且，如果多個(gè)表面被跟蹤，一個(gè)表面可能在一個(gè)位置處與另一表面相重疊，而在另一位置出現(xiàn)相反的關(guān)系。這些情況可總稱為拓?fù)洳灰恢?。迄今為止公布的方法——其中一些在下面進(jìn)行了概述——很大程度上忽視了拓?fù)湟恢滦?topologyconsistency)。在“TheBinaryConsistencyCheckingSchemeandItsApplicationstoSeismicHorizonDetection,"IEEETransactionsonPatternAnalysisandMachineIntelligence,11,439-447(1989)中，Cheng和Lu描述了從二維數(shù)據(jù)中抽取地震骨架的方法。第三維引入的問題既沒討論也沒解決。該程序采用迭代方法，其中強(qiáng)的層位被首先跟蹤，而較弱的層位在后面的迭代中進(jìn)行跟蹤。在任意迭代時(shí)，跟蹤被局限于較早迭代中已經(jīng)跟蹤的層位所描繪的區(qū)域。通過將多個(gè)鄰近的跡線同時(shí)關(guān)聯(lián)進(jìn)行跟蹤。結(jié)合兩種方法使得將地質(zhì)組構(gòu)(geologicfabric)合并成結(jié)果。該方法還被描述在〃AnIterativeApproachtoSeismicSkeletonization,"LuandCheng,Geophysics55,1312-1320(1990)中。在"SeismicSkeletonization:ANewApproachtoInterpretationofSeismicReflectionData,"JournalofGeophysicalResearch-SolidEarth102,8427-8445(1997)中，LLVasudevan和Cook描述了利用地震骨架解釋地震數(shù)據(jù)的應(yīng)用。地震骨架是兩維的，并且當(dāng)層位裂開時(shí)，跟隨哪一個(gè)分支的決定沒有地質(zhì)學(xué)根據(jù)。相反，該方法嘗試以使傾角變化最小化的方式使三條鄰近跡線的事件相關(guān)聯(lián)。該方法僅包括層位的迭代生長。此夕卜，“AdaptationofSeismicSkeletonizationforOtherGeoscienceApplications,“Vasudevan,Eaton,andCook,GeophysicalJournalInternational162,975-993(2005)是早期工作的繼續(xù)，實(shí)現(xiàn)了骨架化(剖面線條化，skeletonization)具有除地震處理和解釋以外的地球科學(xué)應(yīng)用。在"BranchAndBoundSearchForAutomaticLinkingProcessOfSeismicHorizons,“Huang,PatternRecognition23,657-667(1990)中，Huang公開了層位生長的兩維方法，其允許層位相互交叉和穿透，這違反了地質(zhì)地層不交叉的地層學(xué)范例。該方法只揭示了地層的產(chǎn)生，這是通過挑選事件例如峰，在這些事件之間建立所有可能連接的樹，然后選擇產(chǎn)生最線性層位的事件。對(duì)連接樹的分支進(jìn)行選擇以使層位非線性的價(jià)格函數(shù)(costfunction)最小化。“HowToCreateAndUse3DWheelerTransformedSeismicVolumes,"deGroot,deBruin,andHemstra,SEG2006公開了一種解釋方法，所述方法通過依照局部的傾角和走向用重采樣分辨(二次抽樣決定，sub-sampleresolution)內(nèi)插層位，將這些層位以連續(xù)的順序組織，并且通過沿著層位將層位或?qū)傩泽w(attributevolume)整平使沉積域中這些層位或其上的屬性可視。具體地，該算法要求輸入需要用可選方法諸如手工挑選法挑選的主要層位。在主要層位包括(bracket)的層段內(nèi)，次要層位或者與頂部層位或底部層位平行地內(nèi)插，在其間線性地內(nèi)插，或者沿著從地震屬性估計(jì)的局部傾角和走向。通過構(gòu)建，內(nèi)插的次要層位沒有相互交叉。在為2008年6月9-12日意大利羅馬舉行的70thEAGE(EuropeanAssociationofGeoscientistsandEngineers)ConferenceandExhibition提交白勺且從2008年5^26HJf^nT/Awww.earthdoc.orgTicl^llW^^^"AnApproachofSeismicInterpretationBasedonCognitiveVision"的論文中，Verney等公開通過采用基于“認(rèn)知視覺(cognitivevision)”的人工智能工具以地質(zhì)學(xué)為基礎(chǔ)解釋地震數(shù)據(jù)的方法。利用地震數(shù)據(jù)中的體素連通性(voxelconnectivity)檢測(cè)第一順序反射層連續(xù)性。然后，進(jìn)行視覺表征步驟。例如，基于反射層位于彼此的上方或者下方建立地質(zhì)年代順序關(guān)系。最后，通過融合所有節(jié)點(diǎn)——(a)共享類似的視覺屬性(振幅、厚度、傾角)，并且(b)位于與至少一個(gè)其它反射層相似的距離處，從反射層中確定地質(zhì)層位。結(jié)果是一組按年代順序排列的層位。授予Dunn和Czernuszenko的美國專利號(hào)7，024，021，“MethodforPerformingStratigraphically-BasedSeedDetectionina3-DSeismicDataVolume，，公開了一種三維地質(zhì)體挑選器和分析器。在該專利中，挑選幾種精選的地質(zhì)體，其可包括屬性值在指定范圍內(nèi)的地質(zhì)體或者鄰近某屬性值的地質(zhì)體。在挑選期間，地質(zhì)體采用圖示標(biāo)準(zhǔn)(mapviewcriteria)進(jìn)行分析以檢測(cè)和消除自身重疊地質(zhì)體，并且代替地產(chǎn)生復(fù)合地質(zhì)體(compositegeobody)。復(fù)合地質(zhì)體至少滿足沒有自身重疊的拓?fù)錀l件，但地質(zhì)體之間的邊界通過體素(voxel)被檢測(cè)到的順序進(jìn)行確定。在“SystemandMethodforDisplayingSeismicHorizonswithAttributes"(PCT專利申請(qǐng)公布號(hào)WO2007046107)中，James公開了一種地震自動(dòng)挑選器，其產(chǎn)生單值層位并且當(dāng)層位裂開時(shí)通常選取正確的分支。解釋程序通過在地震數(shù)據(jù)體積內(nèi)手工選擇一個(gè)或多個(gè)種子點(diǎn)(seedpoint)將該方法初始化。該算法利用種子點(diǎn)從鄰近跡線中挑選一組次要點(diǎn)——其然后被作為新的種子點(diǎn)進(jìn)行處理，并且重復(fù)該程序。導(dǎo)致自身重疊的次要挑選被舍棄，但是沒有顯示與其它層位的拓?fù)湟恢滦?。該算法主要基于受控的前進(jìn)式(marching)。授予Cacas的美國專利號(hào)7,257,488(“MethodofSedimentologicInterpretationbyEstimationofVariousChronologicalScenariosofSedimentaryLayersDeposition")公開了一種采用上/下關(guān)系將地震和地質(zhì)層位組織成層次以便于其地層學(xué)解釋的方法。該方法通過利用對(duì)沉積層沉積的實(shí)際地質(zhì)年代情況的估計(jì)自動(dòng)從地震數(shù)據(jù)中抽取沉積學(xué)解釋有關(guān)的信息。該算法通過設(shè)定地震數(shù)據(jù)閾值并采用形態(tài)學(xué)細(xì)化以產(chǎn)生單個(gè)層位而開始。如果多個(gè)層位交叉，那么將最線性的一對(duì)組合而將其它的明確分離。該方法然后迭代地估計(jì)沉積層沉積的第一和第二地質(zhì)年代情況，分別假定各個(gè)反射層沉淀在沉積層沉積過程期間的最早和最晚可能時(shí)刻。以參考層位開始，該算法基本上列舉了之上和之下的層位以建立相對(duì)順序。最終進(jìn)行這兩個(gè)地質(zhì)年代情況的解釋以便重建沉積層的沉積情況。相對(duì)順序上的差異被用于估計(jì)情況(scenario)不確定性。授予Cacas的英國專利號(hào)2,444,167(“MethodforStratigraphicInterpretationofSeismicImages")公開了一種地震圖象的地質(zhì)學(xué)解釋以便確定地下的沉積歷史的方法。該方法涉及自動(dòng)跟蹤產(chǎn)生至少一個(gè)層位的事件，選擇具有在所述層位處或附近的窗口抽取的類似地震屬性的層位，以及沿所選擇的層位將地震體整平。授予Borgos的美國專利號(hào)7，248，539(〃ExtremaClassification")公開了一種通過波形簇和碎片特性中共同的成員(commonmembership)進(jìn)行層位碎片(horizonpatch)形成和合并的方法。該方法通過抽取例如所有的峰挑選層位，但是通過將波形簇集使它們相關(guān)聯(lián)。屬于同一簇的挑選被用于定義按照特性諸如簇指標(biāo)、位置或地震屬性合并成較大層位的層位碎片。具體地，方法通過3D地震輸入體(inputvolume)的極值表示法用二次抽樣精密度定義了地震層位的位置。對(duì)于各個(gè)極值，它得自表示極值位置附近地震波形的形狀的系數(shù)，并且通過利用這些系數(shù)無人監(jiān)督的或有監(jiān)督的分類將極值位置分成具有相似波形的組。它然后抽取原始表面作為表面層段，其沿著地震體的極值空間上連續(xù)并且在分類體的分類指標(biāo)中連續(xù)。通過對(duì)與每個(gè)碎片相連的特性諸如分類指標(biāo)、位置、屬性值等進(jìn)行過濾，一組碎片可被組合成最終的層位解釋。揭示了原始表面的三個(gè)基本應(yīng)用將原始表面組合成完整層位進(jìn)行解釋；將地震體內(nèi)的閉合容積(closedvolume)定義為垂直排列的原始表面的閉包(closure)；或者基于原始表面估計(jì)斷層位移。Monsen等(〃Geologic-process-controlledinterpretationbasedon3DWheelerdiagramgeneration,"SEG2007)通過抽取碎片的上/下關(guān)系將授予的Borgos美國專利號(hào)7，248,539進(jìn)行了延伸，并且通過應(yīng)用拓?fù)浞诸惱眠@些關(guān)系得出滿足這些限制的碎片的相對(duì)順序。整平的層位然后被以該相對(duì)順序放置以允許在沉積Wheeler域中解釋。SEG摘要是在2008年6月12日公布的美國專利申請(qǐng)公布號(hào)US2008/0140319的基礎(chǔ)。授予Pedersen的英國專利號(hào)2，375，448(〃ExtractingFeaturesfromanImagebyAutomaticSelectionofPixelsAssociatedwithaDesiredFeature,Pedersen")公開了從幾個(gè)精選種子點(diǎn)構(gòu)造表面的方法，諸如層位和斷層。該方法通過產(chǎn)生慢慢會(huì)聚成線(兩維)或表面(三維)的許多路徑在種子點(diǎn)之間內(nèi)插或者從種子點(diǎn)進(jìn)行外推。該方法基于螞蟻離開群體搜尋食物的方式。最初，它們的路徑幾乎是任意的，但是每個(gè)螞蟻都留下信息素痕跡。螞蟻跟隨各自的氣味，并且隨著時(shí)間的過去，出現(xiàn)了短的成功路徑。該策略適于層位跟蹤，其中成功由沿著路徑地震數(shù)據(jù)的一致性限定。對(duì)于斷層挑選，成功似乎由沿著該路徑的不連貫限定。隨著時(shí)間過去，單個(gè)區(qū)段生長，并且一些可以合并以形成更大的表面。在后續(xù)步驟中，片段根據(jù)其方向和預(yù)計(jì)軌道進(jìn)行連接。授予Viswanathan的美國專利號(hào)5，570，106(“MethodandApparatusforCreatingHorizonsfrom3-DSeismicData")公開了一種計(jì)算機(jī)輔助的層位挑選方法，其通過允許用戶刪除部分層位并且用剩余的層位作為種子點(diǎn)進(jìn)行自動(dòng)挑選進(jìn)行的。授予Sitoh的美國專利號(hào)5，537，365(〃ApparatusandMethodforEvaluationofPickingHorizonsin3-DSeismicData")公開了一種評(píng)價(jià)層位挑選質(zhì)量的方法，其通過應(yīng)用不同的挑選策略和參數(shù)以允許結(jié)果的相互校驗(yàn)進(jìn)行的。授予Stark的美國專利號(hào)6，853，922("SystemForInformationExtractionFromGeologicTimeVolumes")公開了將地震數(shù)據(jù)轉(zhuǎn)化成沉積相對(duì)地質(zhì)時(shí)間域的方法。該方法基于地震瞬間相數(shù)據(jù)的展開。授予Keskes等的美國專利號(hào)6，771，800("MethodofChrono-StratigraphicInterpretationofASeismicCrossSectionOrBlock")公開了將地震數(shù)據(jù)轉(zhuǎn)換成沉積或年代地層域的方法。它們構(gòu)建虛擬反射層，使地震剖面或地震體離散，計(jì)算每個(gè)像素或體素中虛擬反射層的個(gè)數(shù)，并且重正化(renormalizing)柱狀圖。通過對(duì)每條跡線執(zhí)行該程序，它們產(chǎn)生剖面或體，其中每個(gè)水平薄片近似為指示在一個(gè)時(shí)間沉積的地質(zhì)層的層位。該剖面或體然后被用于將數(shù)據(jù)轉(zhuǎn)換成沉積或年代地層域。然而，該參考文獻(xiàn)沒有公開表面的產(chǎn)生，也沒有公開表面的破裂或合并，還沒有公開拓?fù)鋵W(xué)或拓?fù)湟恢滦?。需要的是從地?或?qū)傩?數(shù)據(jù)或任何地球物理學(xué)數(shù)據(jù)產(chǎn)生拓?fù)湟恢碌姆瓷鋵游?topologicallyconsistentreflectionhorizon)的方法,優(yōu)選的是同時(shí)產(chǎn)生多個(gè)層位的方法。本發(fā)明實(shí)現(xiàn)了該需要。
發(fā)明內(nèi)容在一種實(shí)施方式中，本發(fā)明可以是將地震或地震屬性數(shù)據(jù)體(seismicattributedatavolume)中確定的表面合并以形成表示地下地質(zhì)結(jié)構(gòu)或物質(zhì)的地球物理狀態(tài)的較大表面的方法，其包括以拓?fù)湟恢碌姆绞胶喜⑧徑砻妗Ｔ谝恍?shí)施方式中，拓?fù)湟恢驴杀欢x為驗(yàn)證表面滿足下列中每一個(gè)(i)沒有自身重疊；(ii)局部一致性；和(iii)整體一致性。在更詳細(xì)的實(shí)施方式中，方法可以是將地震測(cè)量中獲得的地震數(shù)據(jù)體轉(zhuǎn)換成相應(yīng)的數(shù)據(jù)體的計(jì)算機(jī)執(zhí)行的方法，所述相應(yīng)的數(shù)據(jù)體當(dāng)視覺顯示時(shí)顯示出地下反射層表面的圖，反射層表面通過反射地震波產(chǎn)生數(shù)據(jù)，其中所述方法包括(a)從數(shù)據(jù)體中挑選地震反射，并從挑選中產(chǎn)生初始表面；(b)將表面破裂成拓?fù)湟恢抡贾鲗?dǎo)的較小部分(“碎片(patches)“)；(c)以拓?fù)湟恢碌姆绞胶喜⑧徑槠?，由此從地震?shù)據(jù)體中抽取基于拓?fù)湟恢路瓷涞谋砻妫缓?d)顯示抽取的表面(即骨架)以視覺檢查或解釋，或者將其數(shù)字表示保存到計(jì)算機(jī)存儲(chǔ)器或數(shù)據(jù)存儲(chǔ)器。任選地，將一次迭代的步驟(c)的表面用在下一次迭代的步驟(b)中，重復(fù)步驟(b)-(c)可以至少一次。在上面的步驟(a)中，通過使地震數(shù)據(jù)體中鄰近跡線之間的反射事件相關(guān)，可以挑選地震反射。相關(guān)可以用交叉事件相似系數(shù)或相關(guān)系數(shù)作為相關(guān)量度連接峰和谷，其中如果相關(guān)量度大于預(yù)先選擇的閾值，則連接是接受的，但是如果小于該閾值，則連接被舍棄。在本發(fā)明的一些實(shí)施方式中，只有獨(dú)特的相關(guān)性被接受?？蛇x地，可被確定而且也被公認(rèn)為多重相關(guān)連接，其特征在于從一個(gè)峰、谷或零交叉的兩個(gè)或更多個(gè)相關(guān)，都超過所述閾值。在步驟(C)中將鄰近碎片合并之前，碎片可針對(duì)拓?fù)湟恢滦赃M(jìn)行編輯并且拓?fù)洳灰恢碌乃槠杀粍h除，或者引起不一致的數(shù)據(jù)體素可被刪除。在上面的步驟(b)中，將表面破裂成碎片可通過下列步驟完成將初始表面收縮成線，在線中去除接點(diǎn)以形成更多單個(gè)線，將單個(gè)線收縮成單體素點(diǎn)(特征點(diǎn))，和通過添加鄰近的體素沿著初始表面?zhèn)鞑ヌ卣鼽c(diǎn)以形成體素碎片。野火般的傳播(wildfirepropagation)可被用于沿著初始表面?zhèn)鞑c(diǎn)，例如在每個(gè)特征點(diǎn)周圍圓周添加一個(gè)體素厚的依次更大的層，每個(gè)傳播被限于相應(yīng)特征點(diǎn)從其收縮的表面。連續(xù)的圓周添加體素可在不同碎片相遇的地方停止，因此防止任何體素屬于超過一個(gè)的碎片。傳播可以受到限制，從而任何碎片中的所有體素在收縮之前可追蹤回到相同的初始表面?？梢圆煌绞竭M(jìn)行收縮，例如通過形態(tài)學(xué)細(xì)化。線收縮成點(diǎn)可以通過同時(shí)從各端以相同的速度收縮所述線而實(shí)現(xiàn)。面收縮成線可以通過中間軸變換進(jìn)行。如果在點(diǎn)傳播期間，由于缺乏拓?fù)湟恢滦?，點(diǎn)被舍棄以添加到碎片，其可以被指定另外的特征點(diǎn)。在更一般的實(shí)施方式中，本發(fā)明可以是探測(cè)烴的方法，包括(a)從地震測(cè)量中獲得產(chǎn)生的地震或地震屬性數(shù)據(jù)的數(shù)據(jù)體；(b)將數(shù)據(jù)體細(xì)分成部分，稱為目標(biāo)(任選地，該步驟可通過前述段落的骨架化方法進(jìn)行)；(c)形成一個(gè)或多個(gè)目標(biāo)的區(qū)域；(d)開發(fā)或選擇量度，以便按照表示地質(zhì)體、界面或這些的交叉點(diǎn)，或其它物理地質(zhì)結(jié)構(gòu)或物質(zhì)的地質(zhì)物理狀態(tài)——其指示烴沉積——的可能性將區(qū)域歸類；和(e)利用該量度把區(qū)域區(qū)分優(yōu)先次序，然后利用所述區(qū)分優(yōu)先次序評(píng)定所述體的烴可能性。在另一實(shí)施方式中，本發(fā)明可以是從地下區(qū)域生產(chǎn)烴的方法。方法包括(a)獲得表示該地下區(qū)域的地震數(shù)據(jù)體；(b)至少部分基于通過上述骨架化方法從地震數(shù)據(jù)體中抽取的基于拓?fù)湟恢路瓷涞谋砻娅@得烴在地下區(qū)域中積聚的可能性的預(yù)測(cè)；和(c)響應(yīng)烴可能性的肯定預(yù)測(cè)，在該地下區(qū)域中鉆井并且生產(chǎn)烴。此外，方法的一種或多種實(shí)施方式可包括利用基于拓?fù)湟恢路瓷涞谋砻骖A(yù)測(cè)或分析烴積聚的可能性；其中拓?fù)湟恢率侵赶铝兄兄辽僖环N(i)沒有自身重疊；(ii)局部一致性，例如一個(gè)表面不可能在一個(gè)位置在第二表面上方但在另一位置在其下面；和(iii)整體一致性，是指例如對(duì)于三個(gè)表面A、B和C，如果A在B上面并且B在C上面，則在任意位置C不可能在A上面；其中拓?fù)湟恢率侵?i)、(ii)和(iii)所有三種情況；其中地震反射通過使地震數(shù)據(jù)體中鄰近跡線之間的反射事件相關(guān)而進(jìn)行挑選；其中相關(guān)可以用交叉事件相似系數(shù)或相關(guān)系數(shù)作為相關(guān)量度連接數(shù)據(jù)峰和谷，其中如果相關(guān)量度大于預(yù)先選擇的閾值則連接是接受的，但是如果小于該閾值則連接被舍棄；其中利用計(jì)算機(jī)自動(dòng)化挑選；并且其中在合并鄰近的碎片之前，碎片針對(duì)拓?fù)湟恢滦赃M(jìn)行編輯，并且拓?fù)洳灰恢碌乃槠粍h除，或者引起不一致的數(shù)據(jù)體素被刪除。而且，該方法的一種或多種實(shí)施方式可包括其中將表面破裂成碎片，包括將初始表面收縮成線，在線中去除接點(diǎn)以形成更多單個(gè)線，將單個(gè)線收縮成單體素點(diǎn)(特征點(diǎn))，通過添加鄰近的體素沿著初始表面?zhèn)鞑ヌ卣鼽c(diǎn)以形成體素碎片；其中各個(gè)特征點(diǎn)用不同的標(biāo)記進(jìn)行標(biāo)記，并且標(biāo)記應(yīng)用于在特征點(diǎn)周圍形成的碎片，因此提供了當(dāng)不同的碎片通過傳播延伸時(shí)追蹤不同的碎片的方法；其中野火般的傳播被用于沿著初始表面?zhèn)鞑c(diǎn)，其包括在每個(gè)特征點(diǎn)周圍圓周添加一個(gè)體素厚的依次更大的層，每個(gè)傳播被限于相應(yīng)特征點(diǎn)從其收縮的表面；其中連續(xù)的圓周添加體素在不同碎片相遇的地方停止，因此防止任何體素屬于超過一個(gè)的碎片；其中傳播可以受到限制，從而任何碎片中的所有體素在收縮之前可追蹤回到相同的初始表面；其中受控的前進(jìn)式被用于沿著初始表面?zhèn)鞑c(diǎn)；其中初始表面收縮成線包括從表面外圍接連去除一個(gè)體素厚的層直至產(chǎn)生單個(gè)體素的連續(xù)線；進(jìn)一步包括從線中刪除接點(diǎn)體素以形成更多的線，然后將線收縮成點(diǎn)；其中線收縮成點(diǎn)通過同時(shí)從各端以相同的速度收縮所述線而實(shí)現(xiàn)；其中收縮通過形態(tài)學(xué)細(xì)化進(jìn)行；其中表面收縮成線通過中間軸變換進(jìn)行；其中在點(diǎn)傳播期間要執(zhí)行拓?fù)湟恢滦裕黄渲杏捎谌狈ν負(fù)湟恢滦员簧釛壧砑拥剿槠狞c(diǎn)被指定另外的特征點(diǎn)。另外，方法的一種或多種實(shí)施方式可包括其中以拓?fù)湟恢碌姆绞胶喜⑧徑乃槠ㄟ^下列步驟進(jìn)行開發(fā)碎片的重疊和鄰居列表(overlapandneighbortable)，通過將重疊和鄰居列表分類產(chǎn)生合并對(duì)候選者的次序，利用重疊和鄰居列表校驗(yàn)候選者合并的拓?fù)湟恢滦?，以及接受拓?fù)湟恢碌暮喜ⅲ黄渲朽従恿斜淼姆诸惔涡蚧卩徑槠膸缀涡螤罨蛑g的幾何形狀差別，或者基于從分配有碎片的地震數(shù)據(jù)中抽取的一種或多種屬性的統(tǒng)計(jì)學(xué)特性或之間的差別；其中只有獨(dú)特的相關(guān)性被接受；確定并且也認(rèn)可多重相關(guān)連接，特征在于從單個(gè)峰、谷或零交叉的兩個(gè)或更多個(gè)相關(guān)，其都超過所述閾值；利用基于拓?fù)湟恢路瓷涞谋砻鎸⒒谕負(fù)湟恢路瓷涞谋砻婵臻g整平成表示沉積順序的次序，和利用整平的表面預(yù)測(cè)或分析烴積聚的可能性；整平其中存在基于拓?fù)湟恢路瓷浔砻娴南嚓P(guān)地震數(shù)據(jù)；其中地震數(shù)據(jù)整平通過地震數(shù)據(jù)的非線性伸長或者通過剪貼方法進(jìn)行；其中利用計(jì)算機(jī)自動(dòng)化進(jìn)行每一步；將一次迭代的步驟(c)的表面用在下一次迭代的步驟(b)中，重復(fù)步驟(b)-(c)至少一次；產(chǎn)生視覺表示法(即樹)，其顯示基于拓?fù)湟恢路瓷涞谋砻娴某练e次序或?qū)哟危焕盟鰳溥x擇一個(gè)或多個(gè)表面可視化；利用碎片將地震數(shù)據(jù)體分割成三維體或面間包(inter-surfacepackages)，其表示在共同層段內(nèi)沉積的地質(zhì)單元，利用它們分析烴可能性；分析基于拓?fù)湟恢路瓷涞谋砻娴倪吘壓徒K止點(diǎn)的位置和特征并且用其輔助預(yù)測(cè)或分析烴積聚的可能性；分析基于拓?fù)湟恢路瓷涞谋砻娴膶傩院蛶缀翁匦院?或所述表面位置處相關(guān)地震數(shù)據(jù)以輔助預(yù)測(cè)或分析烴積聚的可能性，其利用碎片或基于拓?fù)湟恢路瓷涞谋砻嬉来螠p少地震數(shù)據(jù)體中包含的信息量，從而減少地震數(shù)據(jù)后續(xù)數(shù)據(jù)處理的存儲(chǔ)或計(jì)算效率要求；以及其中合并鄰近碎片限于收縮前追蹤回同一初始表面的碎片。通過參考下面的具體實(shí)施方式和附圖將會(huì)更好地理解本發(fā)明及其優(yōu)點(diǎn)，其中圖1是準(zhǔn)備用于例如三維地震層位跟蹤的解釋的地震振幅數(shù)據(jù)體的計(jì)算機(jī)顯示；圖2顯示了四百五十(450)個(gè)對(duì)應(yīng)于峰和谷反射層位的表面，其通過本發(fā)明方法從圖1的地震數(shù)據(jù)體抽取，所有表面都是拓?fù)湟恢碌?；圖3A-3C圖解一個(gè)或多個(gè)層或表面之間的三種拓?fù)洳灰恢拢粓D4是一個(gè)流程圖，其顯示用于地震數(shù)據(jù)體拓?fù)涔羌芑谋景l(fā)明方法的一種實(shí)施方式；圖5A-D是圖解圖4中步驟的示意圖6顯示通過跨過鄰近跡線跟蹤峰獲得的地震反射表面；圖7顯示利用通過本發(fā)明方法的一種實(shí)施方式產(chǎn)生的一致表面組建立其整體次序并且將地震數(shù)據(jù)(在數(shù)據(jù)體內(nèi))重新組織成該次序用于解釋的方法步驟；圖8是一個(gè)流程圖，其顯示圖4方法的具體實(shí)施方式中基本步驟；圖9A-9F圖解將多值表面轉(zhuǎn)化成一致的表面的圖8流程圖的示例性應(yīng)用；圖10A-10C圖解事件跟蹤以及事件跟蹤后通過填充表面中的間隙進(jìn)行編輯；圖11圖解通過收縮(細(xì)化)從圖示中示意性的原始的、潛在多值的表面到線然后到特征點(diǎn)；接著點(diǎn)標(biāo)記并且將標(biāo)記傳播回到表面上的進(jìn)程；圖12A是一個(gè)流程圖，其顯示在本發(fā)明方法的一種實(shí)施方式中鄰近碎片的拓?fù)浜喜⒊蓪?duì)的基本步驟；圖12B是一個(gè)流程圖，其顯示在本發(fā)明方法的另一種實(shí)施方式中鄰近碎片的拓?fù)浜喜⒊蓪?duì)的基本步驟；圖13圖解由事件之間低的相關(guān)性或者多個(gè)好的但不明確的連接的存在所引起的孔，兩者都可在編輯步驟中進(jìn)行修整；圖14是將地震數(shù)據(jù)體的時(shí)間/深度垂直刻度轉(zhuǎn)換成推斷的地層學(xué)布置和沉積的水平次序的基本步驟的流程圖；圖15顯示基于圖13的表面的拓?fù)浯涡蚝退降膶?shí)例；圖16A顯示圖2數(shù)據(jù)體的四百五十個(gè)表面的表面水平樹，圖16B顯示所述樹的四個(gè)水平部分的放大圖，以及圖16C顯示與四個(gè)連續(xù)水平相關(guān)的圖2的所有表面；圖17是圖16A的表面水平樹中表面的拓?fù)洳淮_定性的圖示；圖18顯示通過本發(fā)明方法確定其表面然后從地球物理時(shí)間域轉(zhuǎn)化成(拓?fù)?水平域之后圖2的數(shù)據(jù)體；圖19圖解將地震體轉(zhuǎn)換成水平或次序域的方法；圖20顯示圖1的水平轉(zhuǎn)換的地震數(shù)據(jù)；圖21是一個(gè)流程圖，其顯示用于高級(jí)化(high-grading)地質(zhì)目標(biāo)的本發(fā)明方法的一種實(shí)施方式中的基本步驟；以及圖22A-B顯示沿著表面抽取的地震振幅上的兩個(gè)表面的深度等高線。本發(fā)明將結(jié)合實(shí)例實(shí)施方式進(jìn)行描述。對(duì)于下列具體實(shí)施方式針對(duì)本發(fā)明的特定實(shí)施方式或特定應(yīng)用來說明確的程度，這僅僅意欲是例證性的，不解釋為限定本發(fā)明的范圍。相反，它意欲覆蓋可能包括在如所附權(quán)利要求所限定的本發(fā)明范圍內(nèi)的所有改變、修飾或等價(jià)形式。具體實(shí)施例方式為了尋找地球中的烴積聚，地球?qū)W家正采用遙感方法以查看地球表面以下。常規(guī)使用的技術(shù)是地震反射方法，其中人造聲波在所述表面附近產(chǎn)生。聲音傳播到地球中，并且無論何時(shí)聲音從一個(gè)巖石層傳到另一巖石層，一小部分聲音被反射回到記錄它的表面。一般地，采用數(shù)百個(gè)至數(shù)千個(gè)記錄器械。聲波在許多不同位置被相繼激發(fā)。從所有這些記錄中，兩維QD)或三維(3D)的地下圖像可以在數(shù)據(jù)處理后獲得。地震解釋經(jīng)常涉及挑選表面以表征地下以便描繪與烴的探測(cè)、鑒定和生產(chǎn)相關(guān)的地下特征。本發(fā)明描述了同時(shí)挑選多個(gè)表面的方法。也就是說，本發(fā)明方法的實(shí)施方式可被用于同時(shí)挑選這些表面中的多個(gè)或所有表面。同時(shí)挑選多個(gè)表面的能力(即，骨架化地震數(shù)據(jù)的能力)能使圖像識(shí)別或機(jī)器學(xué)習(xí)方法(machinelearningmethod)搜尋直接表明烴或烴系統(tǒng)的要素諸如儲(chǔ)層、封閉層、源、熟化和遷移的地質(zhì)或地球物理數(shù)據(jù)以測(cè)定和描繪可能的烴積聚。在地球物理或地質(zhì)解釋的應(yīng)用中，在術(shù)語“層位”和“表面”之間經(jīng)常進(jìn)行區(qū)分。如本文所用，表面和層位可以互換使用。本發(fā)明是同時(shí)產(chǎn)生多個(gè)表面，同時(shí)使單個(gè)表面為單值的和使所有表面拓?fù)湟恢碌姆椒?。采用傳統(tǒng)方法的為多值的或者拓?fù)洳灰恢碌谋砻嬗靡唤M較小碎片代替，其中每一碎片是單值的并且與所有其它表面是拓?fù)湟恢碌?。該方法產(chǎn)生表示地震數(shù)據(jù)體中包含的多個(gè)或所有反射表面的表面。它產(chǎn)生地震數(shù)據(jù)的骨架化表示，其大大減少了數(shù)據(jù)量。有利地，它以地質(zhì)學(xué)直觀的方式組織和呈現(xiàn)地震數(shù)據(jù)，其促進(jìn)地震解釋和地下表征以及因此描繪與烴的探測(cè)和生產(chǎn)相關(guān)的地下特征。圖1呈現(xiàn)了地震振幅體(振幅體積，amplitudevolume)的實(shí)例。使從一條跡線至下一條的峰(亮的)或谷(暗的)相關(guān)聯(lián)允許限定表面。采用本發(fā)明方法的一種實(shí)施方式，對(duì)于圖1的這個(gè)實(shí)例體可抽取圖2中顯示的四百五十個(gè)表面。圖1和2和其它類似圖上顯示的數(shù)字網(wǎng)格表示地震測(cè)量的離散坐標(biāo)，其由源和接受器位置確定。許多地震表面對(duì)應(yīng)于地下巖層之間的界面。各個(gè)層是在大約同一時(shí)間沉積的一堆巖石。假定兩個(gè)并列的層，較深的層產(chǎn)生得較早，并且較淺的層較晚。地層學(xué)，即巖石層順序的科學(xué)，表明這種關(guān)系空間上持續(xù)。如果在一個(gè)位置一層覆蓋在另一層上，那么它在其存在的各個(gè)地方都覆蓋這一層。主要的例外是由結(jié)構(gòu)復(fù)雜性諸如掩沖斷層、逆斷層或倒轉(zhuǎn)褶皺引起的。在本發(fā)明的至少一種實(shí)施方式中，拓?fù)湟恢率侵戈P(guān)于巖層的幾何學(xué)排列滿足下列三個(gè)條件。1.巖石層不可以自身重疊。如果一個(gè)層自身重疊，則它同時(shí)比自身和夾在之間的巖石更晚和更早(更年輕和更老)。這種表述可被稱為圖3A中所圖解的沒有自身重疊的條件。2.兩個(gè)層不可以顛倒它們的沉積關(guān)系。一個(gè)層不可以在一個(gè)位置處位于另一層的上方并且在另一位置處位于它的下方。否則，一個(gè)層既比另一層更早又更晚。這種表述可被稱為圖3B中所圖解的局部一致性的條件。3.層系列必須保持傳遞性。上/下或更晚/更早是傳遞性關(guān)系的實(shí)例。如果層一在層二上方，并且層二在層三上方，那么層三必須在層一之下。否則，層一既比層三更早又比層三更晚。這種表述可被稱為圖3C中所圖解的整體一致性的條件?？梢宰⒁?，沒有自身重疊條件是局部一致性條件的特殊情形，并且局部一致性條件是整體一致性條件的特殊情形。然而，第一條件核對(duì)起來比另兩個(gè)條件更容易，并且第二條件核對(duì)起來比第三條件更容易。出于計(jì)算效率，分開處理所有三個(gè)條件是有用的，即使第三條件實(shí)際上合并了其它條件?？蛇x地，可以限定沒有自身重疊條件以便它適用于一個(gè)表面，可以限定局部一致性條件以便它僅適用于涉及兩個(gè)不同表面的時(shí)候，并且可以限定整體一致性條件==以便從而它僅適用于涉及至少三個(gè)不同表面的時(shí)候，其中三個(gè)條件是相斥的。如果地震反射事件是由聲波從一個(gè)層傳到另一層引起的，并且因此通常與巖層之間的界面相關(guān)，那么地震反射表面也需要滿足這三個(gè)條件。對(duì)地震數(shù)據(jù)的任何相旋轉(zhuǎn)形式(phaserotatedversion)來說也是如此，盡管在這種數(shù)據(jù)中的反射事件不必需與邏輯學(xué)界面相關(guān)。對(duì)于一組具有關(guān)聯(lián)的上/下關(guān)系的表面，三個(gè)上述條件可用于檢查這些表面的整體一致性。違反所述條件的表面或者不是由巖層引起的，或者被不正確地跟蹤。與層不相關(guān)的表面包括斷層、流體接觸、或者薄層反射混合的事件。跟蹤錯(cuò)誤可涉及噪聲、地震采集和處理人為因素或薄層調(diào)整。對(duì)于給定系列的層(或表面)，上/下(或更晚/更早)關(guān)系的集成限定了它們的拓?fù)浣Y(jié)構(gòu)。滿足三個(gè)條件的至少一個(gè)，優(yōu)選滿足所有三個(gè)條件的一系列層被稱為拓?fù)湟恢?。在下面給出的實(shí)例實(shí)施方式的討論中，其中在它有關(guān)的上下文中，拓?fù)湟恢率侵笣M足了所有的三個(gè)條件。對(duì)于拓?fù)湟恢碌膶酉盗?，不同事件的整體次序可通過針對(duì)這些關(guān)系進(jìn)行拓?fù)浞诸愡M(jìn)行限定(例如，Skiena,TheAlgorithmDesignManual,Springer,273-274(1998))。一般地，不應(yīng)用本發(fā)明方法的實(shí)施方式，由于層(或表面)之間沖突的關(guān)系，所以建立表面次序是有問題的和/或不可能的。這些拓?fù)洳灰恢乱话阋鹜負(fù)浞诸愂　＿@個(gè)論點(diǎn)可以被反過來測(cè)試拓?fù)湟恢滦匀绻⑶抑挥腥绻負(fù)浞诸惓晒?，表面才是一致的。本發(fā)明的目標(biāo)之一是建立表面之間的一致性。如果拓?fù)浞诸惓晒?，那么表面是拓?fù)湟恢碌摹Ｈ绻負(fù)浞诸愂?，那么表面是拓?fù)洳灰恢碌摹４送?，拓?fù)浞诸愃惴ㄨb定引起不一致的表面。一致性不暗示著所得到的表面次序是唯一的。例如，兩個(gè)小的、鄰近的但不重疊的表面是拓?fù)湟恢碌牟⑶覍?dǎo)致成功的拓?fù)浞诸?。然而，所得的線性分類次序是不唯一的，即任意表面可以首先列出而不違背任何的上/下限制或條件。許多小表面比幾個(gè)大的表面更有可能是拓?fù)湟恢碌?。在小尺寸限制下，每個(gè)表面在垂向和橫向上僅延伸一個(gè)點(diǎn)，并且因此，通過構(gòu)造，這些單點(diǎn)表面是拓?fù)湟恢碌?。本發(fā)明方法的實(shí)施方式部分基于該觀察。圖4是流程圖，其顯示在地震數(shù)據(jù)體骨架化的本發(fā)明方法的一種實(shí)施方式中的基本步驟。一般地，地震體是全疊加振幅體，但是可以使用任何地震或地震屬性體。在步驟41，跟蹤地震反射表面穿過地震體以找到原始的、潛在多值的表面。在該背景下，地震事件是峰(屬性最大值)、谷(屬性最小值)、從峰到谷交叉的零、或者從谷到峰交叉的零。一種或多種的所有事件被挑選并且與鄰近地震跡線上的事件相關(guān)。在本實(shí)例中，峰和谷都被挑選。因此，這個(gè)步驟涉及從數(shù)據(jù)體中挑選地震反射，并且從挑選中產(chǎn)生初始表面。在步驟42中，通過步驟41產(chǎn)生的表面被破裂成一組小的碎片。這些碎片優(yōu)選地足夠小，以致于它們相互拓?fù)湟恢抡贾鲗?dǎo)，并且拓?fù)洳灰恢碌哪切┧槠赏ㄟ^擦去幾個(gè)單個(gè)點(diǎn)(即數(shù)據(jù)體素或像素)或者甚至刪除整個(gè)小碎片——其產(chǎn)生拓?fù)洳灰恢隆菀椎刂苽?。因此，該步驟涉及將表面破裂成拓?fù)湟恢抡贾鲗?dǎo)的較小部分(“碎片”)。在步驟43，由多個(gè)小碎片通過合并鄰近的碎片產(chǎn)生較大的表面。如步驟42中所提供的，所有碎片是拓?fù)湟恢碌?。在本?shí)例實(shí)施方式中，對(duì)于每一碎片測(cè)定它重疊哪些碎片以及它是在這些碎片中每一個(gè)的上方還是下方。此外，對(duì)于每一碎片，其鄰居(即，相似水平上的碎片，其包含與正在分析的碎片鄰近的跡線)被鑒定。如果所得組合沒有引起拓?fù)洳灰恢拢敲脆徑槠赡軐儆谙嗤谋砻娌⑶冶缓喜?。該步驟可被稱為拓?fù)浜喜⒊绦?。因此，該步驟涉及以拓?fù)湟恢碌姆绞胶喜⑧徑槠虼藦牡卣饠?shù)據(jù)體中抽取基于拓?fù)湟恢路瓷涞谋砻?。在一個(gè)、多個(gè)或所有鄰近碎片被拓?fù)浜喜⒑螅Y(jié)果是構(gòu)建一組拓?fù)湟恢碌谋砻?。它們被?chǔ)存在計(jì)算機(jī)存儲(chǔ)器中以用于解釋和表征地下。在地震數(shù)據(jù)體的一些區(qū)域中，跟蹤程序(優(yōu)選地自動(dòng)的)可能使跡線之間的事件誤相關(guān)(miscorrelate)。在其它區(qū)域，差的數(shù)據(jù)質(zhì)量可防止地震事件跟蹤器關(guān)聯(lián)某些事件。最后，一些關(guān)聯(lián)可能如此不明確的以致于它們不可能歸屬于單個(gè)表面。在這些情況的每一種情況中，由周圍一致表面所提供的局部構(gòu)造可有助于解決這些問題。誤相關(guān)可被糾正，噪音區(qū)域中差的關(guān)聯(lián)可變得可接受，或者可使多重關(guān)聯(lián)消除歧義。來自步驟43的一致性表面組可允許改進(jìn)地震事件跟蹤并且，如果期望，進(jìn)一步通過工作流程(圖4中由虛線迭代箭頭所指示的)可以進(jìn)行以填充孔并產(chǎn)生較少的更廣闊的表面(即，改良的骨架)。圖5A-5D是圖解圖4的步驟的示意圖。在圖5A中，地震數(shù)據(jù)體中的峰事件被跟蹤，并且(圖5B)發(fā)現(xiàn)形成多值表面。圖5C顯示破裂成十六個(gè)相互拓?fù)湟恢碌男∷槠?1的表面。在圖5D中，鄰近碎片被合并成較大的碎片，除非這引起拓?fù)洳灰恢隆Ｗ詈蟮慕Y(jié)果52是一組四個(gè)拓?fù)湟恢碌谋砻?，每個(gè)由不同的交叉影線顯示。圖6顯示通過跨過鄰近跡線跟蹤峰獲得的地震反射表面的實(shí)例。在左側(cè)，表面是單值的，但是在右側(cè)，表面明顯是多值的并且自身重疊至少兩次。許多現(xiàn)有的地震自動(dòng)跟蹤器或者產(chǎn)生這種多值表面，或者僅僅返回每個(gè)位置的不同可能性之一，一般是首先找到的可能性，并且因此不一定是地質(zhì)相關(guān)的可能性。圖7呈現(xiàn)了本發(fā)明方法的一種實(shí)施方式的應(yīng)用，其中采用一組拓?fù)湟恢碌谋砻?，諸如本發(fā)明方法所產(chǎn)生的，重新組織地震屬性體。因?yàn)楸砻媸且恢碌?，所以至少有一種遵從(honor)單個(gè)上/下關(guān)系的次序。如果表面對(duì)應(yīng)于地質(zhì)地層之間的邊界，那么這種次序表示它們沉積的順序。一般地，次序不是唯一的，因?yàn)樾〉奶卣骺赡軝M向上不連接而沒有重疊，并且因此它們精確的次序不可能被建立。以這樣的方式垂直地扭曲地震數(shù)據(jù)(例如，整平地震表面)，以致于相應(yīng)的地震表面以這樣的次序排列以允許分析器以地質(zhì)地層可能已經(jīng)沉積的次序分析地震數(shù)據(jù)，這有利于烴的探測(cè)和生產(chǎn)。接下來，本發(fā)明方法被更詳細(xì)地解釋，如圖4方法的特定實(shí)施方式所圖解的，其被圖解在圖8的流程圖和圖9中該流程圖的示例性應(yīng)用中。在圖8中，“階段1”是指圖4的步驟41，“階段2”是指步驟42，以及“階段3”是指步驟43。圖9A-9F圖解圖8中流程圖的實(shí)例應(yīng)用并且可通過同時(shí)觀察圖8更好地理解。在圖9A中，多值表面通過跟蹤地震事件構(gòu)建(階段1)。在階段2(步驟42)中，表面被簡化成線(圖9B中顯示)，線接點(diǎn)被去除(圖9C中顯示)，并且線被簡化成特征點(diǎn)(圖9D中顯示)。剩余的點(diǎn)被標(biāo)記并且標(biāo)記被傳播回到表面上以構(gòu)建碎片(圖9E中顯示)。所得碎片在圖8的階段3(步驟43)中被拓?fù)浜喜?，產(chǎn)生一致的表面(圖9F中顯示)。階段1(步驟41)步驟81的第一部分是事件跟蹤。在本發(fā)明的這種實(shí)施方式中，所有事件的跟蹤涉及相關(guān)鄰近事件并且編輯空隙和誤相關(guān)。相關(guān)通過跨過所有跡線抽取所有期望的地震事件或反射表面開始。期望的事件可包括峰、谷、任一種零交叉(+/_或_/+)。經(jīng)驗(yàn)已經(jīng)表明抽取峰和谷(但不是零交叉)可能是使事件總數(shù)最小化(即，使計(jì)算效率最大化)和使所得表面質(zhì)量最大化之間良好的折衷。使用超過一種的事件減少了事件相關(guān)和拓?fù)浜喜⒅械牟淮_定性，因?yàn)榉濉⒐群土憬徊嫔⒉荚谡麄€(gè)地震數(shù)據(jù)體中。圖IOA圖解在本發(fā)明的這種實(shí)施方式中進(jìn)行的事件跟蹤。如圖的左邊所顯示，對(duì)于每條地震跡線，局部最小值被抽取以限定谷(虛線箭頭)，而局部最大值限定峰(實(shí)線箭頭)。由括號(hào)指示的地震跡線窗口101在各個(gè)事件中間，并且被用于不同跡線之間的事件相關(guān)。附圖的右手邊部分圖解不同跡線之間的事件窗口相關(guān)，并且因此構(gòu)建原始表面。一維的(導(dǎo)向的(Pilot))數(shù)據(jù)信息包102A(例如，峰的中心)與鄰近跡線中的其它信息包進(jìn)行比較。一些事件呈現(xiàn)大的類似性并且是獨(dú)特相關(guān)的(實(shí)線箭頭)。其它事件可與鄰近跡線上的超過一個(gè)的事件很好地相關(guān)(103箭頭)并且將被稱為多重相關(guān)。例如，在拓?fù)浜喜⒑蠡蛘咴诘诙谓?jīng)過工作流程時(shí)為了消除模凌兩可情況儲(chǔ)存兩種相關(guān)的位置和特性之后顯示兩種相關(guān)(在本發(fā)明的這種實(shí)施方式中)，而不是相比另一種選擇一種有效的相關(guān)。一些相關(guān)可能是差的(104箭頭)。通過考慮它們的背景和周圍的局部地震構(gòu)造在拓?fù)浜喜⒑?，僅僅具有差的相關(guān)的事件可被分配到表面。跡線間相關(guān)可采用許多方法諸如交叉相關(guān)或相似性分析數(shù)學(xué)上進(jìn)行測(cè)量。好的相關(guān)優(yōu)選地被定義為超過預(yù)先確定的閾值的相關(guān)，而差的相關(guān)是沒有超過預(yù)先確定的閾值的相關(guān)。另外的標(biāo)準(zhǔn)，諸如鄰近事件之間的垂直距離(滯后)，也可被用在相關(guān)過程期間。如果這種滯后超過預(yù)先確定的閾值，那么這兩個(gè)事件最有可能屬于不同的表面，并且沒有進(jìn)行連接(即，它們的相關(guān)被排除)。有利地，這可被用于防止周波跳躍。更一般地，跡線間相關(guān)可被計(jì)算為跡線間公制的結(jié)果。這可由限定計(jì)算多維向量空間中的距離的函數(shù)組成。例如，挑選在兩條跡線中每一條上居中的兩個(gè)窗口限定由那個(gè)窗口內(nèi)的值組成的多維向量。這些值可以是在各個(gè)體素上記錄的振幅或者從那些振幅計(jì)算的多個(gè)特征(例如，統(tǒng)計(jì)，諸如平均、方差和更高的矩；通過傅里葉變換的頻率分解；等)。比較兩個(gè)向量的函數(shù)可以是歐幾里得距離函數(shù)、1-范數(shù)、豪斯多夫距離等。實(shí)際上，兩個(gè)信息包通常不直接連接，這是因?yàn)樗鼈兊南嚓P(guān)性差或者它們的垂直差別超過一定的閾值。然而，它們可如圖10B-10C中所圖解以間接方式進(jìn)行明確的連接。圖IOB顯示具有幾個(gè)缺少的連接的表面。如果不考慮的話，這種缺少的連接引起許多不必要的碎片，其增加拓?fù)浜喜⒌挠?jì)算成本。然而，空隙可進(jìn)行修整，在那里已經(jīng)暗含連接。在本發(fā)明的一種實(shí)施方式中，當(dāng)事件可被唯一連接時(shí)，盡管以間接方式(即，沒有滿足先前討論的相關(guān)標(biāo)準(zhǔn))，可以明確進(jìn)行直接連接以關(guān)閉空隙并且因此防止不必要的碎片的產(chǎn)生。該編輯步驟(步驟82)依賴于這樣的事實(shí)以一個(gè)方向包圍空隙在鄰近跡線上導(dǎo)致與另一方向包圍空隙相同的點(diǎn)，這表明表面局部上是簡單的并且既沒裂開也沒螺旋。例如，考慮連接路徑105和106，其每個(gè)經(jīng)過繞圖IOA中空隙的相反路線。路徑105在同一位置結(jié)束，其暗示著那一點(diǎn)和路徑中的點(diǎn)之間唯一的連接。那些缺少的連接通過圖IOC中添加的兩個(gè)新單元邊界線(較粗的線)顯示。相反，路徑106顯示這些缺少的連接是不明確的，并且因此在圖IOC中的位置沒有進(jìn)行改變。階段2(步驟4第二階段是產(chǎn)生拓?fù)湟恢碌乃槠?。通過跟蹤由峰、谷和/或零交叉限定的反射事件在階段1中獲得的原始表面一般不是拓?fù)湟恢碌?。它們通?)自身重疊，幻在一個(gè)位置處存在于另一表面上方，但在不同位置處在同一表面下方(局部不一致)，或者；3)是在它們的上/下關(guān)系中包含回路的一組表面的一部分(整體不一致)。與幾個(gè)大的碎片相比，許多較小的碎片更有可能是拓?fù)湟恢碌?。事?shí)上，如果所有碎片僅僅是面積延伸中的一個(gè)樣品，那么通過構(gòu)建，它們是拓?fù)湟恢碌摹Ｒ虼?，該階段的目標(biāo)是將原始的、潛在多值的表面破裂成較小的、拓?fù)湟恢碌乃槠＿@通過應(yīng)用中間軸變換或形態(tài)學(xué)細(xì)化首先簡化(收縮)表面成拓?fù)漕愃频木€而進(jìn)行(步驟83)(例如參見HaralickandShapiro,ComputerandRobotVision,Vol.1,Chapter5，Addison-Wesley(1992))。因?yàn)橐?-連接指向0-connectedsense)應(yīng)用細(xì)化，線段之間的接點(diǎn)的特征在于具有至少三個(gè)直接的鄰居。去除接點(diǎn)，接著第二次應(yīng)用形態(tài)學(xué)細(xì)化，簡化(收縮)原始表面成幾個(gè)離散的特征點(diǎn)，其容易被給予獨(dú)特的標(biāo)識(shí)符或標(biāo)記。在步驟84，分配的標(biāo)記然后被傳播回到原始表面——該過程可更加描述性地稱為后向傳播，但是也可簡單地稱為傳播。圖11顯示通過形態(tài)學(xué)細(xì)化從圖示(111)中示意性的多值的表面到線(112)并且去除線接點(diǎn)(113)，通過形態(tài)學(xué)細(xì)化從線到點(diǎn)(114)，點(diǎn)標(biāo)記(115)，并且將標(biāo)記傳播回表面(116)的行程。這組具有相同標(biāo)記的事件定義碎片。圖11中的結(jié)果是步驟116中一組八個(gè)小碎片，其對(duì)應(yīng)于步驟115中八個(gè)特征點(diǎn)。通過構(gòu)建，碎片等于或小于它們的母體表面，因?yàn)槊總€(gè)特征點(diǎn)與附近的特征點(diǎn)競(jìng)爭連接成員。標(biāo)記的后向傳播可例如采用簡單的野火算法(wildfirealgorithm)或者采用受控前進(jìn)式算法(controlledmarchingalgorithm)-其受例如事件相關(guān)性、滯后(事件的垂直接近程度)或者表面的局部曲率指導(dǎo)——進(jìn)行。與簡單的野火算法相比，受控前進(jìn)式的優(yōu)點(diǎn)是它能更快地傳播標(biāo)記跨過較平的區(qū)域而較慢地移動(dòng)跨過復(fù)雜區(qū)域，由此產(chǎn)生更均質(zhì)的碎片。可以想象其他傳播方法，其在本發(fā)明方法的范圍內(nèi)。傳播標(biāo)記后，所得碎片盡管大體一致但不保證拓?fù)湟恢隆榱烁玫剡M(jìn)行階段3中的拓?fù)浜喜?，這些碎片優(yōu)選地被調(diào)整成拓?fù)湟恢碌摹ｈb定拓?fù)洳灰恢碌膬?yōu)選方式是通過構(gòu)建重疊表格(步驟85)開始的，以記錄哪些碎片重疊其它碎片以及它們重疊的方式。在步驟86時(shí)，鑒定不一致。在表格構(gòu)建或隨后的檢查期間，自身重疊的表面是很明顯的。從表中，鑒定出具有沖突的上-下關(guān)系(即局部不一致)的碎片對(duì)。最后，通過嘗試將重疊表格中剩余條目拓?fù)浞诸?，技術(shù)人員可發(fā)現(xiàn)上下關(guān)系是圓形(整體不一致)的三個(gè)或更多個(gè)碎片的組。如果沒有圓形關(guān)系存在，那么拓?fù)浞诸惓晒?。如果這種整體不一致存在，那么拓?fù)浞诸愂遣豢赡艿牟⑶掖娴胤祷鼐哂胁灰恢玛P(guān)系的碎片列表。步驟86的最后部分是編輯鑒定的拓?fù)洳灰恢碌乃槠?。最簡單的編輯方法是刪除不一致的碎片。另外的切除方法是通過僅僅去除沖突重疊來源來剪除這些碎片。該方法需要一些勤奮，因?yàn)橐恍┧槠ㄟ^該過程可能變得不連接，并且可能需要重新標(biāo)記所得片。另一編輯方法可以是反復(fù)將不一致的碎片裂開成更小的碎片直到消除了所有的不一致。實(shí)際上，簡單地刪除不一致的碎片似乎很奏效，這是因?yàn)橐恢碌乃槠h(yuǎn)多于不一致的碎片，并且不一致的碎片一般小很多并且經(jīng)常位于邊緣區(qū)域。編輯不一致的表面碎片后，優(yōu)選的是重新構(gòu)建重疊表格以說明這些編輯變化。階段3(步驟4第三階段涉及在合并的表面保持拓?fù)湟恢碌臈l件下合并鄰近碎片成較大的碎片。第一任務(wù)是測(cè)定哪些片段被連接(即，在數(shù)據(jù)體中以某些方式相互鄰接，但被不同標(biāo)記)。這些碎片被稱為鄰居，并且可被記錄(步驟87)在鄰居列表中作為拓?fù)浜喜⒊奢^大碎片的候選者，最終產(chǎn)生一個(gè)表面。例如，單獨(dú)的碎片通過細(xì)化(例如簡化、收縮)產(chǎn)生并且簡化成不同的特征點(diǎn)。如果表面是完美的長方形，其中在長方形內(nèi)所有方向上具有完美的連接，那么細(xì)化很可能產(chǎn)生五個(gè)特征點(diǎn)，并且傳播后由此產(chǎn)生五個(gè)碎片。僅僅因?yàn)樗鼈兪遣煌乃槠?，并不暗示著它們沒有很好相關(guān)性地相互連接。正要合并的大部分碎片是良好相關(guān)的表面的一部分。一般地，有許多碎片和許多對(duì)鄰居。所得的拓?fù)湟恢碌谋砻娴臄?shù)量、形狀和質(zhì)量取決于合并候選者被評(píng)價(jià)的次序。以重疊的兩個(gè)碎片以及與兩者都相鄰的第三碎片為例。第三碎片不可能與兩者都合并，這是因?yàn)樗煤喜⒌谋砻孀陨碇丿B。因此，它僅僅與一個(gè)合并。特定的選擇指示隨后合并的成功或失敗。繼續(xù)步驟87，鄰居列表中的鄰居對(duì)被優(yōu)選排列成進(jìn)行合并嘗試的次序。普通的次序是僅僅數(shù)字漸增的標(biāo)記的次序(即，遇到鄰居的次序)。更成熟的排序可以合并碎片特性，諸如鄰近碎片中事件之間的相關(guān)系數(shù)和碎片方位間的相似性。后者是建立合并次序的優(yōu)選方法。方位類似的鄰近碎片首先被合并，這是因?yàn)樗鼈兏锌赡鼙憩F(xiàn)共同的地質(zhì)地層，而方位很不同的鄰近碎片最后被合并，因?yàn)樗鼈兛赡芘c噪聲假象或非地層學(xué)事件諸如斷層或流體接觸相關(guān)。甚至更高級(jí)的排序可基于在碎片位置處或其附近抽取的次要地震屬性之間的統(tǒng)計(jì)學(xué)相似性。在合并次序通過一種方法或另一種方法建立的情況下，可以采取拓?fù)浜喜?步驟88)。本發(fā)明一種實(shí)施方式的方法被詳細(xì)描繪在圖12A的流程圖中，并且接下來被描述；第二種實(shí)施方式被描述在圖12B中并且在下面被進(jìn)一步描述。在步驟121，一對(duì)鄰近碎片被選為合并候選者，并且被假定構(gòu)成或者是一個(gè)表面的一部分，這是指一個(gè)碎片的重疊關(guān)系適用于另一個(gè)(反之亦然)。如果該行為產(chǎn)生拓?fù)洳灰恢?，那么該合并必須被舍棄。否則，合并被接受，并且重疊和鄰居列表通過用另一碎片的標(biāo)記替換一個(gè)碎片的標(biāo)記而進(jìn)行調(diào)整。假定的合并后評(píng)估三種一致性條件的計(jì)算成本很不同。自身重疊校驗(yàn)起來快且容易。這在圖12A中顯示為步驟122。局部一致性檢查需要檢查整個(gè)重疊列表(步驟123)。然而，整體一致性檢查需要拓?fù)浞诸?步驟124)，其是計(jì)算上昂貴的。在圖12A的實(shí)施方式中，三種一致性檢查以數(shù)字成本增加的次序級(jí)聯(lián)。更昂貴的檢查只對(duì)經(jīng)過較少成本檢查的合并候選者進(jìn)行。如果拓?fù)浞诸惓晒?步驟125)，那么合并的碎片整體一致，并且因此拓?fù)湟恢?。然后，接受該假定，并且因此修改列?步驟126)。然后，對(duì)下一對(duì)合并候選者重復(fù)這個(gè)程序(步驟121)。如果分類或任何其它測(cè)試失敗，那么(步驟128)假定被舍棄，并且所述程序應(yīng)用于下一對(duì)(步驟121)。甚至級(jí)聯(lián)三種一致性檢查計(jì)算成本高，這是因?yàn)橥負(fù)浞诸愋枰粓?zhí)行許多次。不對(duì)每一對(duì)鄰近碎片都執(zhí)行拓?fù)浞诸?，拓?fù)渌槠喜⑺惴杀伙@著加速。該算法的一個(gè)修正是引入隊(duì)列。經(jīng)過第一和第二測(cè)試的鄰居對(duì)(步驟122和12被置于一隊(duì)列中，而不立刻通過第三測(cè)試(步驟124)進(jìn)行評(píng)價(jià)。一旦隊(duì)列達(dá)到用戶指定的大小，例如四對(duì)，該重疊列表被復(fù)制，所有建議的合并被應(yīng)用于該副本，并且進(jìn)行拓?fù)浞诸悺Ｈ绻诸惓晒?，那么所有四種建議的合并是整體一致的并且可接受的。如果分類失敗，那么在建議的合并中一定有至少一個(gè)不一致。為了找到該不一致，原始重疊列表被再次復(fù)制，但只有前兩個(gè)合并被應(yīng)用到該副本中。剩余對(duì)只是被儲(chǔ)存在保持隊(duì)列(holdingqueue)中。如果分類成功，那么前兩對(duì)的合并是可接受的，并且知道在后兩對(duì)中存在不一致。如果分類失敗，就知道在前兩對(duì)中存在不一致。對(duì)包含不一致的那一組再次重復(fù)所述程序，但這次只評(píng)價(jià)一對(duì)。拓?fù)浞诸惡?，就立刻知道可能哪一?duì)導(dǎo)致不一致并且應(yīng)當(dāng)被舍棄。此時(shí)，再次進(jìn)行分類之前通過再充填隊(duì)列重復(fù)該循環(huán)。換句話說，在隊(duì)列中積聚的建議合并中發(fā)現(xiàn)不一致后，隊(duì)列被切成兩份，并對(duì)每一份進(jìn)行分類直到發(fā)現(xiàn)不一致，同時(shí)接受成功的合并。一般地，隊(duì)列應(yīng)當(dāng)不限于四對(duì)，而是幾百對(duì)或幾千對(duì)。而且，隊(duì)列大小可允許動(dòng)態(tài)變化。如果分類失敗，則隊(duì)列大小被減少，但是如果它成功了，那么隊(duì)列大小被增加用于接下來拓?fù)浞诸惖脑u(píng)價(jià)。在N對(duì)之間尋找一個(gè)不一致可以進(jìn)行Iog2N次分類而不是N次分類。對(duì)于具有一千零二十四個(gè)要素的隊(duì)列，一個(gè)不一致可以在至多十次分類中找到，其導(dǎo)致計(jì)算成本大大降低。拓?fù)浜喜⒌牡诙N實(shí)施方式被顯示在圖12B中，其中細(xì)節(jié)被呈現(xiàn)在表1中。本發(fā)明的這種可選實(shí)施方式在不一致檢查進(jìn)行的方式上不同于前面的實(shí)施方式。第一種方法檢查在合并兩個(gè)表面碎片后是否引入了有向循環(huán)(directedcycle)。相比之下，可選實(shí)施方式預(yù)測(cè)合并是否產(chǎn)生有向循環(huán)而不是檢查循環(huán)。這是計(jì)算量較不密集的多的任務(wù)，其不但執(zhí)行相同的功能，而且是也更有效。圖12B的方法輸入的是初始碎片、它們的排序(非循環(huán)有向圖)以及合并次序列表(鄰近碎片對(duì))。輸出的是較大的碎片，并且最終是表面。只要一個(gè)碎片在另一個(gè)上方，合并兩個(gè)表面碎片后就產(chǎn)生循環(huán)(即拓?fù)洳灰恢?。因此，為了避免引入不一致，它有能力檢查一個(gè)碎片是否在另一個(gè)碎片上方。提供這種關(guān)系方便表示的數(shù)據(jù)結(jié)構(gòu)是一種有向圖所述體內(nèi)表面碎片通過節(jié)點(diǎn)描繪，并且如果一個(gè)碎片在另一個(gè)上方，則節(jié)點(diǎn)之間的定向連接或邊緣就存在。因此，問題簡化成具體的圖遍歷問題(graphtraversalproblem)，其中所述問題是兩個(gè)節(jié)點(diǎn)(表面碎片)之間的路徑是否存在。圖遍歷問題可采用標(biāo)準(zhǔn)深度優(yōu)先搜尋(cbpth-firstsearch,DFS)算法(例如，參見IntroductiontoAlgorithms,Cormen,LeisersonandRiverst,MITPressandMcGrawHill,477-485(section23.3，"Depth-firstsearch")(1990))進(jìn)行解決。對(duì)這種一般算法執(zhí)行下列修正達(dá)到實(shí)際上更好的計(jì)算效率。首先，增加數(shù)據(jù)結(jié)構(gòu)、有向圖，其中在每個(gè)節(jié)點(diǎn)的另外屬性，u，表示為DEPTHATT(u)，其跟蹤碎片的絕對(duì)深度。第二，引入幾何深度特性(GeometricDepthProperty,⑶P)并且改進(jìn)遍歷算法以便確保⑶P在所有時(shí)間被保持和使用(圖12B中的步驟1201以及表1中的程序1)。該特性(GDP)要求當(dāng)沿圖中的定向邊緣時(shí)深度屬性單調(diào)增加。換句話說，如果在體中碎片a重疊碎片b，那么碎片a的深度屬性必須小于或等于碎片b的深度屬性。然后，在步驟1202中，從合并列表中選擇一對(duì)碎片，并且在步驟1203中，兩個(gè)所選的碎片的合并采用GDP檢查拓?fù)湟恢滦砸栽黾有?表1中程序2)。如果檢查是肯定的，那么在步驟1204中兩個(gè)碎片按照表1的程序3進(jìn)行合并。該方法是有效的，因?yàn)閮蓚€(gè)節(jié)點(diǎn)之間路徑的搜尋被限于圖的一小部分而不是整個(gè)結(jié)構(gòu)待合并的表面碎片具有在有限數(shù)值范圍內(nèi)的深度屬性值，并且搜尋只探測(cè)深度屬性在那范圍內(nèi)的節(jié)點(diǎn)。⑶P保證這是充分的。表.權(quán)利要求1.將地震測(cè)量中獲得的地震數(shù)據(jù)體轉(zhuǎn)換成相應(yīng)的數(shù)據(jù)體的計(jì)算機(jī)執(zhí)行的方法，所述相應(yīng)的數(shù)據(jù)體當(dāng)視覺顯示時(shí)顯示地下反射層表面的圖，所述地下反射層表面通過反射地震波產(chǎn)生數(shù)據(jù)，所述方法包括(a)從所述數(shù)據(jù)體中挑選地震反射，并從所述挑選中產(chǎn)生初始表面；(b)將表面破裂成拓?fù)湟恢抡贾鲗?dǎo)的較小部分(“碎片")；(c)以拓?fù)湟恢碌姆绞胶喜⑧徑槠?，由此從所述地震?shù)據(jù)體中抽取基于拓?fù)湟恢路瓷涞谋砻?；?d)顯示抽取的表面以視覺檢查或解釋，或者將它們的數(shù)字表示保存到計(jì)算機(jī)存儲(chǔ)器或數(shù)據(jù)存儲(chǔ)器。2.權(quán)利要求1的方法，進(jìn)一步包括利用所述基于拓?fù)湟恢路瓷涞谋砻骖A(yù)測(cè)或分析烴積聚的可能性。3.權(quán)利要求1的方法，其中拓?fù)湟恢掳?yàn)證表面滿足下列中至少一個(gè)(i)沒有自身重疊；(ii)局部一致性；和(iii)整體一致性。4.權(quán)利要求1的方法，其中所述地震反射通過使所述地震數(shù)據(jù)體中鄰近跡線之間的反射事件相關(guān)聯(lián)而挑選。5.權(quán)利要求4的方法，其中所述挑選是利用計(jì)算機(jī)自動(dòng)進(jìn)行的。6.權(quán)利要求1的方法，其中將表面破裂成碎片包括將初始表面收縮成線，去除所述線中的接點(diǎn)以形成更多單個(gè)線，將單個(gè)線收縮成單體素點(diǎn)(特征點(diǎn))，和通過添加鄰近的體素沿著所述初始表面?zhèn)鞑ニ鎏卣鼽c(diǎn)以形成體素碎片。7.權(quán)利要求6的方法，其中各個(gè)特征點(diǎn)用不同的標(biāo)記進(jìn)行標(biāo)記，并且所述標(biāo)記應(yīng)用于在所述特征點(diǎn)周圍形成的碎片，因此提供了當(dāng)不同的碎片通過傳播延伸時(shí)追蹤不同的碎片的方法。8.權(quán)利要求6的方法，其中受控的前進(jìn)式被用于沿初始表面?zhèn)鞑c(diǎn)。9.權(quán)利要求6的方法，其中初始表面收縮成線包括從表面外圍接連去除一個(gè)體素厚的層直至產(chǎn)生單個(gè)體素的連續(xù)線。10.權(quán)利要求6的方法，進(jìn)一步包括在將線收縮成點(diǎn)之前從線中刪除接點(diǎn)體素以形成更多的線。11.權(quán)利要求6的方法，其中在點(diǎn)傳播期間執(zhí)行拓?fù)湟恢滦浴?2.權(quán)利要求1的方法，其中以拓?fù)湟恢碌姆绞胶喜⑧徑乃槠ㄟ^下列步驟進(jìn)行開發(fā)所述碎片的重疊和鄰居列表，通過將所述重疊和鄰居列表分類產(chǎn)生合并對(duì)候選者的次序，利用所述重疊和鄰居列表校驗(yàn)候選者合并的拓?fù)湟恢滦裕约敖邮芡負(fù)湟恢碌暮喜ⅰ?3.權(quán)利要求12的方法，其中所述鄰居列表的分類次序基于所述鄰近碎片的幾何形狀或之間的幾何形狀差別，或者基于從分配有所述碎片的地震數(shù)據(jù)中抽取的一種或多種屬性的統(tǒng)計(jì)學(xué)特性或之間的差別。14.權(quán)利要求1的方法，進(jìn)一步包括利用所述基于拓?fù)湟恢路瓷涞谋砻鎸⑺龌谕負(fù)湟恢路瓷涞谋砻婵臻g整平成表示沉積順序的次序，和利用所述整平的表面預(yù)測(cè)或分析烴積聚的可能性。15.權(quán)利要求14的方法，其進(jìn)一步包括整平其中存在所述基于拓?fù)湟恢路瓷涞谋砻娴南嚓P(guān)地震數(shù)據(jù)。216.權(quán)利要求15的方法，其中所述地震數(shù)據(jù)整平通過所述地震數(shù)據(jù)的非線性伸長或者通過剪貼方法進(jìn)行。17.權(quán)利要求1的方法，進(jìn)一步包括產(chǎn)生視覺表示，其顯示所述基于拓?fù)湟恢路瓷涞谋砻娴某练e次序或?qū)哟巍?8.權(quán)利要求17的方法，其中所述視覺表示是樹，并且進(jìn)一步包括利用所述樹選擇一個(gè)或多個(gè)表面可視化。19.權(quán)利要求1的方法，進(jìn)一步包括利用所述碎片將所述地震數(shù)據(jù)體分割成三維體或面間包，其表示在共同層段內(nèi)沉積的地質(zhì)單元，并且利用它們分析烴可能性。20.權(quán)利要求2的方法，進(jìn)一步包括分析所述基于拓?fù)湟恢路瓷涞谋砻娴倪吘壓徒K止點(diǎn)的位置和特征并且用其輔助預(yù)測(cè)或分析烴積聚的可能性。21.權(quán)利要求2的方法，進(jìn)一步包括分析所述基于拓?fù)湟恢路瓷涞谋砻娴膶傩院蛶缀翁匦院?或所述表面位置處的相關(guān)地震數(shù)據(jù)以輔助預(yù)測(cè)或分析烴積聚的可能性。22.權(quán)利要求1的方法，其進(jìn)一步包括利用所述碎片或基于拓?fù)湟恢路瓷涞谋砻嬉来螠p少地震數(shù)據(jù)體中包含的信息量，從而減少所述地震數(shù)據(jù)后續(xù)數(shù)據(jù)處理的存儲(chǔ)或計(jì)算效率要求。23.權(quán)利要求1的方法，其中合并鄰近碎片限于收縮前追蹤回同一初始表面的碎片。24.權(quán)利要求12的方法，其中利用深度限制搜尋方法在合并鄰近碎片中執(zhí)行拓?fù)湟恢?，所述方法包?a)基于捕獲所述數(shù)據(jù)體中所述碎片的相對(duì)位置的重疊列表產(chǎn)生圖結(jié)構(gòu)；(b)給各個(gè)碎片賦予深度屬性，以便任意兩個(gè)碎片的深度屬性的比較表示出所述碎片之一是否在另一個(gè)上面；(c)利用所述圖結(jié)構(gòu)和所述深度屬性檢查基于所述鄰居列表提出的合并的拓?fù)湟恢滦裕缓?d)當(dāng)碎片合并被接受時(shí)更新所述深度屬性和圖結(jié)構(gòu)。25.權(quán)利要求1的方法，其中所述抽取的表面被顯示或保存為地球模型。26.計(jì)算機(jī)程序產(chǎn)品，其包括計(jì)算機(jī)可用介質(zhì)，其中包含有計(jì)算機(jī)可讀程序代碼，所述計(jì)算機(jī)可讀程序代碼適于被實(shí)行以執(zhí)行將地震數(shù)據(jù)體簡化成基于反射的表面的方法，所述方法包括(a)從所述數(shù)據(jù)體中挑選地震反射，并從所述挑選中產(chǎn)生初始表面；(b)將表面破裂成拓?fù)湟恢滦哉贾鲗?dǎo)的較小部分(“碎片")；和(c)以拓?fù)湟恢碌姆绞胶喜⑧徑槠?，由此從所述地震?shù)據(jù)體中抽取基于拓?fù)湟恢路瓷涞谋砻妗?7.從地下區(qū)域生產(chǎn)烴的方法，所述方法包括(a)獲得表示所述地下區(qū)域的地震數(shù)據(jù)體；(b)至少部分基于通過權(quán)利要求1中所述的方法——其通過引用在此并入——從所述地震數(shù)據(jù)體中抽取的基于拓?fù)湟恢路瓷涞谋砻娅@得在所述地下區(qū)域中烴積聚的可能性的預(yù)測(cè)；和(c)響應(yīng)烴可能性的肯定預(yù)測(cè)，在所述地下區(qū)域中鉆井并且生產(chǎn)烴。28.將在地震或地震屬性數(shù)據(jù)體中鑒定的表面合并以形成表示地下地質(zhì)結(jié)構(gòu)或物質(zhì)的地球物理狀態(tài)的較大表面的方法，其包括以拓?fù)湟恢碌姆绞胶喜⑧徑砻妗?9.探測(cè)烴的方法，包括(a)從地震測(cè)量中獲得產(chǎn)生的地震或地震屬性數(shù)據(jù)的數(shù)據(jù)體；(b)將所述數(shù)據(jù)體細(xì)分成部分，稱為目標(biāo)；(c)形成一個(gè)或多個(gè)目標(biāo)區(qū)域；(d)開發(fā)或選擇量度，以便按照表示地質(zhì)體、界面或這些的交叉點(diǎn)，或其它物理地質(zhì)結(jié)構(gòu)或指示烴沉積的物質(zhì)的地質(zhì)物理狀態(tài)的可能性將區(qū)域歸類；和(e)利用所述量度把區(qū)域區(qū)分優(yōu)先次序，然后利用所述區(qū)分優(yōu)先次序評(píng)定所述體的烴可能性。30.權(quán)利要求四的方法，其中各個(gè)目標(biāo)包含利用基于其數(shù)據(jù)或?qū)傩曰蛘咂渌锢砗侠順?biāo)準(zhǔn)的一個(gè)或多個(gè)標(biāo)準(zhǔn)分類在一起的單元。31.權(quán)利要求四的方法，進(jìn)一步包括利用所述區(qū)域區(qū)分優(yōu)先次序?qū)⑺鰯?shù)據(jù)體轉(zhuǎn)化成地球物理地球模型，并且利用所述地球模型評(píng)定所述體的烴可能性。32.權(quán)利要求四的方法，其中(b)利用權(quán)利要求1中所述的方法——其通過引用在此并入——進(jìn)行。33.權(quán)利要求四的方法，其中(d)中的量度包括直接烴指示(DHI)。全文摘要通過產(chǎn)生表面或地質(zhì)體并分析它們的烴指示來分析地下區(qū)域的烴可能性的方法?；诜瓷涞谋砻婵梢砸酝?fù)湟恢碌姆绞阶詣?dòng)產(chǎn)生，其中單個(gè)表面不自身重疊并且多個(gè)表面的組與地層學(xué)重疊原理一致。初始表面從地震數(shù)據(jù)中挑選(41)，然后破裂成拓?fù)湟恢抡贾鲗?dǎo)的小部分(“碎片”)(42)，在此基礎(chǔ)上鄰近碎片以拓?fù)湟恢碌姆绞胶喜?43)以形成廣闊并一致的表面組(“骨架”)。由此抽取的表面或地質(zhì)體可以基于選定的量度(213)諸如一個(gè)或多個(gè)直接烴指示(“DHI”)例如AVO分類進(jìn)行自動(dòng)分析和評(píng)估(214)。一個(gè)或多個(gè)表面的拓?fù)湟恢滦钥杀欢x為沒有自身重疊加上在多個(gè)表面中局部以及整體一致性(52)。文檔編號(hào)G01V1/40GK102066980SQ200980118674公開日2011年5月18日申請(qǐng)日期2009年4月24日優(yōu)先權(quán)日2008年5月22日發(fā)明者D·G·吉拉德,F·W·施羅德,K·古瑪南,M·J·特雷爾,M·英霍夫,P·迪米特洛夫,S·哈森諾德申請(qǐng)人:?？松梨谏嫌窝芯抗?