可視化地顯示地質(zhì)異常體的方法及裝置制造方法
【專利摘要】本發(fā)明提供了一種可視化地顯示地質(zhì)異常體的方法和裝置，所述方法包括：(a)獲取原始三維地震數(shù)據(jù)；(b)確定原始三維地震數(shù)據(jù)中的有效地質(zhì)異常體數(shù)據(jù)；(c)利用三維圖像種子點生長方法確定有效地質(zhì)異常體數(shù)據(jù)中的地質(zhì)異常體散點數(shù)據(jù)；(d)對所述地質(zhì)異常體散點數(shù)據(jù)進行插值計算，以得到插值后的地質(zhì)異常體散點數(shù)據(jù)；(e)對插值后的地質(zhì)異常體散點數(shù)據(jù)進行成面操作，以得到地質(zhì)異常體成面數(shù)據(jù)模型；(f)基于地質(zhì)異常體成面數(shù)據(jù)模型來可視化地顯示地質(zhì)異常體。所述方法和裝置能夠高效地從原始三維地震數(shù)據(jù)中提取希望得到的地質(zhì)異常體散點數(shù)據(jù)，并實現(xiàn)了對地質(zhì)異常體的可視化顯示。
【專利說明】可視化地顯示地質(zhì)異常體的方法及裝置

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明總體說來涉及地球物理勘探領(lǐng)域，更具體地講，涉及一種可視化地顯示地 質(zhì)異常體的方法和裝置。

【背景技術(shù)】
[0002] 地質(zhì)體一般來說是一種不規(guī)則的物體，它不僅具有自然、任意的形狀和大小，而且 包含多種地質(zhì)構(gòu)造類型，比如礦體、斷層、巖丘、砂體、珊瑚礁、溶洞等。地質(zhì)異常體的概念屬 于地質(zhì)體的范疇，是一種特殊的地質(zhì)體。地質(zhì)異常體不僅具有上述地質(zhì)體的特點，在地球物 理勘探的過程中，存在地質(zhì)異常體的區(qū)域會在很多方面表現(xiàn)出和其它區(qū)域較大的不同。比 如，剖面波形的波峰波谷形狀改變、振幅突變、電法異?；虼朋w異常等。河道砂體、地質(zhì)斷層 以及一些礦體等都屬于地質(zhì)異常體。
[0003] 地質(zhì)異常體的建模一般需要經(jīng)過三個步驟：屬性分析、數(shù)據(jù)點追蹤以及表面模型 生成。具體地講，首先，通過屬性分析求出地質(zhì)異常體的數(shù)據(jù)點的屬性值范圍，然后，利用數(shù) 據(jù)點追蹤找到符合條件的數(shù)據(jù)點，最后，利用符合條件的數(shù)據(jù)點生成表面模型。目前地質(zhì)異 常體的建模方法大體上主要涉及數(shù)字圖像分割和體透視兩類，但這兩類方法都是以標量形 式顯示地質(zhì)異常體且這兩類方法的效率比較低。所以，可以看出，現(xiàn)有的地質(zhì)異常體的建模 方法無法高效地得到地質(zhì)異常體并對地質(zhì)異常體進行顯示。


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

[0004] 本發(fā)明的示例性實施例在于提供一種可視化地顯示地質(zhì)異常體的方法和裝置，以 解決現(xiàn)有技術(shù)中獲得地質(zhì)異常體時效率低下以及無法直觀顯示地質(zhì)異常體的技術(shù)問題。
[0005] 根據(jù)本發(fā)明的一方面，提供一種可視化地顯示地質(zhì)異常體的方法，包括以下步驟： 一種可視化地顯示地質(zhì)異常體的方法，其特征在于，包括以下步驟：(a)獲取原始三維地震 數(shù)據(jù)；(b)確定原始三維地震數(shù)據(jù)中的有效地質(zhì)異常體數(shù)據(jù)；（C)利用三維圖像種子點生長 方法確定有效地質(zhì)異常體數(shù)據(jù)中的地質(zhì)異常體散點數(shù)據(jù)；（d)對所述地質(zhì)異常體散點數(shù)據(jù) 進行插值計算，以得到插值后的地質(zhì)異常體散點數(shù)據(jù)；（e)對插值后的地質(zhì)異常體散點數(shù) 據(jù)進行成面操作，以得到地質(zhì)異常體成面數(shù)據(jù)模型；（f)基于地質(zhì)異常體成面數(shù)據(jù)模型來 可視化地顯示地質(zhì)異常體。
[0006] 在所述方法中，步驟（b)可包括：根據(jù)針對三維地震數(shù)據(jù)的屬性值范圍參數(shù)、坐標 值范圍參數(shù)和采樣間隔參數(shù)中的至少一個確定原始三維地震數(shù)據(jù)中的有效地質(zhì)異常體數(shù) 據(jù)。
[0007] 在所述方法中，步驟（C)可包括：將有效地質(zhì)異常體數(shù)據(jù)中的數(shù)據(jù)點作為種子點， 利用泛洪法基于所述種子點進行三維圖像種子點生長，將經(jīng)過三維圖像種子點生長得到的 所有數(shù)據(jù)點作為地質(zhì)異常體散點數(shù)據(jù)。
[0008] 在所述方法中，在步驟（C)中，可利用六鄰域或二十六鄰域的泛洪法基于所述種 子點進行三維圖像種子點生長。
[0009] 在所述方法中，在步驟（d)中，可基于二值形態(tài)學(xué)對所述地質(zhì)異常體散點數(shù)據(jù)執(zhí) 行先膨脹后腐蝕的插值計算，或者，可基于二值形態(tài)學(xué)對所述地質(zhì)異常體散點數(shù)據(jù)執(zhí)行先 腐蝕后膨脹的插值計算。
[0010] 在所述方法中，原始三維地震數(shù)據(jù)中的每個數(shù)據(jù)點具有各自的坐標信息，所述坐 標信息包括X軸坐標值、y軸坐標值和z軸坐標值，并且，在步驟（d)中，可將所述地質(zhì)異常 體散點數(shù)據(jù)中的X軸坐標值、y軸坐標值分別相同的數(shù)據(jù)點作為同一道數(shù)據(jù)，基于至少具有 兩個數(shù)據(jù)點的每一道數(shù)據(jù)中z軸坐標值最大的數(shù)據(jù)點進行插值以針對每道數(shù)據(jù)得到z軸坐 標值最大的數(shù)據(jù)點，并基于至少具有兩個數(shù)據(jù)點的每一道數(shù)據(jù)中z軸坐標值最小的數(shù)據(jù)點 進行插值以針對每道數(shù)據(jù)得到z軸坐標值最小的數(shù)據(jù)點。
[0011] 在所述方法中，步驟（e)可包括：利用Marching cubes算法對插值后的地質(zhì)異常 體散點數(shù)據(jù)進行成面操作，以得到地質(zhì)異常體成面數(shù)據(jù)模型。
[0012] 在所述方法中，步驟（f)可還包括：在顯示地質(zhì)異常體的頁面中提供地質(zhì)異常體 的體積，其中，所述地質(zhì)異常體的體積基于地質(zhì)異常體成面數(shù)據(jù)模型而計算得到。
[0013] 根據(jù)本發(fā)明的另一方面，提供一種可視化地顯示異常體的方法，包括：(a)獲取原 始三維圖像數(shù)據(jù)；(b)確定原始三維圖像數(shù)據(jù)中的有效異常體數(shù)據(jù)；（c)利用三維圖像種子 點生長方法確定有效異常體數(shù)據(jù)中的異常體散點數(shù)據(jù)；（d)對所述異常體散點數(shù)據(jù)進行插 值計算，以得到插值后的異常體散點數(shù)據(jù)；（e)對插值后的異常體散點數(shù)據(jù)進行成面操作， 以得到異常體成面數(shù)據(jù)模型；（f)基于異常體成面數(shù)據(jù)模型來可視化地顯示異常體。
[0014] 根據(jù)本發(fā)明的另一方面，提供一種可視化地顯示地質(zhì)異常體的裝置，包括：獲取單 元，用于獲取原始三維地震數(shù)據(jù)；有效數(shù)據(jù)確定單元，用于確定原始三維地震數(shù)據(jù)中的有效 地質(zhì)異常體數(shù)據(jù)；散點確定單元，用于利用三維圖像種子點生長方法確定有效地質(zhì)異常體 數(shù)據(jù)中的地質(zhì)異常體散點數(shù)據(jù)；插值單元，用于對所述地質(zhì)異常體散點數(shù)據(jù)進行插值計算， 以得到插值后的地質(zhì)異常體散點數(shù)據(jù)；成面單元，用于對插值后的地質(zhì)異常體散點數(shù)據(jù)進 行成面操作，以得到地質(zhì)異常體成面數(shù)據(jù)模型；顯示單元，用于基于地質(zhì)異常體成面數(shù)據(jù)模 型來可視化地顯示地質(zhì)異常體。
[0015] 上述可視化地顯示地質(zhì)異常體的方法和裝置，能夠高效地從原始三維地震數(shù)據(jù)中 提取希望得到的地質(zhì)異常體散點數(shù)據(jù)，并實現(xiàn)了對地質(zhì)異常體的可視化顯示。

【專利附圖】

【附圖說明】
[0016] 通過下面結(jié)合附圖對本發(fā)明的示例性實施例進行的描述，本發(fā)明示例性實施例的 上述和其他目的和特點將會變得更加清楚，其中：
[0017] 圖1示出根據(jù)本發(fā)明示例性實施例的可視化地顯示地質(zhì)異常體的方法的流程圖；
[0018] 圖2A和圖2B分別為六鄰域種子點生長模型和二十六鄰域種子點生長模型；
[0019] 圖3示出根據(jù)本發(fā)明示例性實施例的對地質(zhì)異常體散點數(shù)據(jù)執(zhí)行先腐蝕后膨脹 的插值計算的步驟的流程圖；
[0020] 圖4示出根據(jù)本發(fā)明示例性實施例的對地質(zhì)異常體散點數(shù)據(jù)執(zhí)行先膨脹后腐蝕 的插值計算的步驟的流程圖；
[0021] 圖5示出根據(jù)本發(fā)明示例性實施例的基于上下表面曲面插值對地質(zhì)異常體散點 數(shù)據(jù)執(zhí)行插值計算的步驟的流程圖；
[0022] 圖6示出根據(jù)本發(fā)明示例性實施例的成面操作的步驟的流程圖；
[0023] 圖7示出根據(jù)本發(fā)明示例性實施例的可視化地顯示異常體的方法的流程圖；
[0024] 圖8示出根據(jù)本發(fā)明示例性實施例的可視化地顯示地質(zhì)異常體的裝置的框圖。

【具體實施方式】
[0025] 現(xiàn)將詳細參照本發(fā)明的實施例，所述實施例的示例在附圖中示出，其中，相同的標 號始終指的是相同的部件。
[0026] 下面將結(jié)合本發(fā)明的具體示例性實施例來對本發(fā)明的可視化地顯示地質(zhì)異常體 的方法進行詳細描述。
[0027] 圖1是根據(jù)本發(fā)明示例性實施例的可視化地顯示地質(zhì)異常體的方法的流程圖。
[0028] 在步驟S100,獲取原始三維地震數(shù)據(jù)。作為示例，原始三維地震數(shù)據(jù)可包括多個數(shù) 據(jù)點，每個數(shù)據(jù)點可具有各自的坐標信息和屬性值?？蛇x地，所述屬性值可為能夠表明數(shù)據(jù) 點的某種特定屬性的值，例如地震屬性值等。
[0029] 在步驟S200,確定原始三維地震數(shù)據(jù)中的有效地質(zhì)異常體數(shù)據(jù)。作為示例，可根據(jù) 針對三維地震數(shù)據(jù)的屬性值范圍參數(shù)、坐標值范圍參數(shù)和采樣間隔參數(shù)中的至少一個來確 定原始三維地震數(shù)據(jù)中的有效地質(zhì)異常體數(shù)據(jù)。具體說來，可預(yù)先獲取三維地震數(shù)據(jù)的參 數(shù)列表，其中，所述參數(shù)列表包括針對三維地震數(shù)據(jù)的屬性值范圍參數(shù)、坐標值范圍參數(shù)和 采樣間隔參數(shù)。這里，三維地震數(shù)據(jù)的參數(shù)列表中的參數(shù)可由用戶設(shè)定。
[0030] 例如，可將原始地震數(shù)據(jù)中滿足屬性值范圍參數(shù)的數(shù)據(jù)點確定為有效地質(zhì)異常體 數(shù)據(jù)?；蛘撸蓪⒃嫉卣饠?shù)據(jù)中滿足坐標值范圍參數(shù)的數(shù)據(jù)點確定為有效地質(zhì)異常體數(shù) 據(jù)。或者，可將對原始地震數(shù)據(jù)根據(jù)采樣間隔參數(shù)進行采樣后得到的數(shù)據(jù)點確定為有效地 質(zhì)異常體數(shù)據(jù)?；蛘撸蓪⑸鲜鰲l件進行組合以確定有效地質(zhì)異常體數(shù)據(jù)。
[0031] 在步驟S300,利用三維圖像種子點生長方法確定有效地質(zhì)異常體數(shù)據(jù)中的地質(zhì)異 常體散點數(shù)據(jù)。
[0032] 作為示例，可將有效地質(zhì)異常體數(shù)據(jù)中的數(shù)據(jù)點作為種子點，利用泛洪法基于所 述種子點進行三維圖像種子點生長，將經(jīng)過三維圖像種子點生長得到的所有數(shù)據(jù)點作為地 質(zhì)異常體散點數(shù)據(jù)。
[0033] 進一步地，作為示例，可利用六鄰域（如圖2A所示的立方體的六個面的中心點） 或二十六鄰域（如圖2B所示的立方體的六個面的中心點、八個頂點和十二條棱的中點）的 泛洪法基于所述種子點進行三維圖像種子點生長。
[0034] 在步驟S400,對地質(zhì)異常體散點數(shù)據(jù)進行插值計算，以得到插值后的地質(zhì)異常體 散點數(shù)據(jù)。
[0035] 作為示例，可基于二值形態(tài)學(xué)對地質(zhì)異常體散點數(shù)據(jù)執(zhí)行先膨脹后腐蝕的插值計 算?；蛘?，可基于二值形態(tài)學(xué)對地質(zhì)異常體散點數(shù)據(jù)執(zhí)行先腐蝕后膨脹的插值計算?；蛘?， 可基于上下表面曲面插值對地質(zhì)異常體散點數(shù)據(jù)執(zhí)行插值計算。
[0036] 在步驟S500,對插值后的地質(zhì)異常體散點數(shù)據(jù)進行成面操作，以得到地質(zhì)異常體 成面數(shù)據(jù)模型。
[0037] 作為示例,可利用Marching cubes算法對插值后的地質(zhì)異常體散點數(shù)據(jù)進行成面 操作，以得到地質(zhì)異常體成面數(shù)據(jù)模型。這里，地質(zhì)異常體成面數(shù)據(jù)模型可由多個等值面組 成，每個等值面對應(yīng)一個預(yù)定閾值。等值面可表示為{(x，y，z)|f(x，y，z) =c}，即，等值面 是空間中所有具有相同預(yù)定閾值的數(shù)據(jù)點的集合，其中，c為所述預(yù)定閾值。
[0038] 在步驟S600,基于地質(zhì)異常體成面數(shù)據(jù)模型來可視化地顯示地質(zhì)異常體。
[0039] 進一步地，作為示例，在顯示地質(zhì)異常體的頁面中可提供地質(zhì)異常體的體積，地質(zhì) 異常體的體積可基于地質(zhì)異常體成面數(shù)據(jù)模型而計算得到。例如，可通過將地質(zhì)異常體劃 分為多個小立方體并將小立方體的體積求和得到地質(zhì)異常體的體積，或者將地質(zhì)異常體分 割為不同空間區(qū)域并將每個空間區(qū)域進行積分運算的結(jié)果求和以得到地質(zhì)異常體的體積。 這里，可以看出，通過地質(zhì)異常體成面數(shù)據(jù)模型可有效地估算地質(zhì)異常體的體積，為后續(xù)對 地質(zhì)異常體進行其他的處理提供便利。
[0040] 由此可見，通過上述步驟能夠高效地從原始三維地震數(shù)據(jù)中提取希望得到的地質(zhì) 異常體散點數(shù)據(jù)，并實現(xiàn)了對地質(zhì)異常體的可視化顯示。
[0041] 以下，將詳細描述步驟S300的實現(xiàn)示例。具體說來，步驟S300可包括以下步驟：
[0042] 步驟一，從有效地質(zhì)異常體數(shù)據(jù)中選取任意一個數(shù)據(jù)點。
[0043] 步驟二，以選取的數(shù)據(jù)點為中心確定一個立方體，將選取的數(shù)據(jù)點壓入第一堆棧。
[0044] 這里，所述立方體的六個面的中心點位于以選取的數(shù)據(jù)點為原點的三維坐標系的 三個坐標軸上，且立方體的邊長為預(yù)先設(shè)定的值（例如，可由用戶根據(jù)需要預(yù)先設(shè)定立方 體的邊長）。
[0045] 步驟三，判斷位于立方體上的預(yù)定點的數(shù)據(jù)點中是否存在屬于有效地質(zhì)異常體數(shù) 據(jù)的數(shù)據(jù)點。作為示例，位于立方體上的預(yù)定點可包括：立方體的六個面的中心點（如圖 2A所示），或者，立方體的六個面的中心點、立方體八個頂點以及立方體的十二條棱的中點 (如圖2B所示）。
[0046] 步驟四，如果不存在屬于有效地質(zhì)異常體數(shù)據(jù)的數(shù)據(jù)點，則首先計算第一堆棧中 的數(shù)據(jù)點的個數(shù)，并將計算出的個數(shù)與預(yù)定個數(shù)閾值（可由用戶設(shè)定）進行比較，如果計算 出的個數(shù)大于所述預(yù)定個數(shù)閾值，則將第一堆棧中的數(shù)據(jù)點壓入第二堆棧，并且，將第一堆 棧中的數(shù)據(jù)點從有效地質(zhì)異常體數(shù)據(jù)中刪除；如果計算出的個數(shù)大于所述預(yù)定個數(shù)閾值， 則直接將第一堆棧中的數(shù)據(jù)點從有效地質(zhì)異常體數(shù)據(jù)中刪除。然后，清空第一堆棧，從有效 地質(zhì)異常體數(shù)據(jù)的剩余數(shù)據(jù)點中選取任意一個數(shù)據(jù)點，繼續(xù)執(zhí)行步驟二，直到有效地質(zhì)異 常體數(shù)據(jù)中沒有剩余的數(shù)據(jù)點為止。
[0047] 步驟五，如果存在屬于有效地質(zhì)異常體數(shù)據(jù)的數(shù)據(jù)點，則針對每個屬于有效地質(zhì) 異常體數(shù)據(jù)的數(shù)據(jù)點再次開始執(zhí)行步驟二，直到步驟四被執(zhí)行為止。
[0048] 這里，當所有步驟執(zhí)行完畢時（即，當有效地質(zhì)異常體數(shù)據(jù)中沒有剩余的數(shù)據(jù)點 時），第二堆棧中存儲的數(shù)據(jù)點即為希望得到的地質(zhì)異常體散點數(shù)據(jù)。
[0049] 關(guān)于圖1所示的步驟S400,根據(jù)本發(fā)明的示例性實施例，對地質(zhì)異常體散點數(shù)據(jù) 進行插值計算的方法包括但不限于：基于二值形態(tài)學(xué)對地質(zhì)異常體散點數(shù)據(jù)執(zhí)行先腐蝕后 膨脹的插值計算、基于二值形態(tài)學(xué)對地質(zhì)異常體散點數(shù)據(jù)執(zhí)行先膨脹后腐蝕的插值計算以 及基于上下表面曲面插值對地質(zhì)異常體散點數(shù)據(jù)進行插值計算。
[0050] 這里，二值形態(tài)學(xué)主要指腐蝕和膨脹運算，其中，腐蝕和膨脹為對偶運算。腐蝕是 一種消除邊界點，以使得邊界向內(nèi)部收縮的過程，根據(jù)本發(fā)明的示例性實施例，可利用腐蝕 的操作來消除無意義的地質(zhì)異常體散點數(shù)據(jù)。膨脹是將與地質(zhì)異常體接觸的所有有效地質(zhì) 異常體數(shù)據(jù)中除了地質(zhì)異常體散點數(shù)據(jù)之外的數(shù)據(jù)點合并到地質(zhì)異常體中，以使得邊界向 外部擴張的過程，根據(jù)本發(fā)明的示例性實施例，可利用膨脹的操作來填補地質(zhì)異常體的空 洞。
[0051] 具體說來，對地質(zhì)異常體散點數(shù)據(jù)進行腐蝕的運算公式如公式（1)所示。
[0052] A QB = {.X： B + X a A} (1)
[0053] 其中，A代表地質(zhì)異常體散點數(shù)據(jù)，B代表三維結(jié)構(gòu)元素，用于掃描A所代表的地質(zhì) 異常體散點數(shù)據(jù)。A0B表示A與B做腐蝕運算的結(jié)果，X表示腐蝕運算結(jié)果中的任一元素， 這里，X應(yīng)滿足下面的條件：將B與X進行"與"運算，運算的結(jié)果仍包含在A內(nèi)。關(guān)于三維 結(jié)構(gòu)元素B，其具有立方體形狀，包含中心元素和周圍元素，例如，中心元素為立方體的中心 點，周邊元素可包括立方體的六個面的中心點，或者，周邊元素可包括立方體的六個面的中 心點、立方體八個頂點以及立方體的十二條棱的中點，中心元素和周邊元素的值在不同的 運算過程中可被設(shè)定為不同的值。此外，立方體的邊長為預(yù)先設(shè)定的值（例如，可由用戶根 據(jù)需要來預(yù)先設(shè)定所述邊長）。
[0054] 對地質(zhì)異常體散點數(shù)據(jù)進行膨脹的運算公式如公式（2)所示。
[0055] A?B= {x:B+x f A} (2)
[0056] 其中，A代表地質(zhì)異常體散點數(shù)據(jù)，B代表三維結(jié)構(gòu)元素，用于掃描A所代表的地質(zhì) 異常體散點數(shù)據(jù)。A ? B表示A與B做膨脹運算的結(jié)果。X表示膨脹運算結(jié)果中的任一元 素，這里，X應(yīng)滿足下面的條件：將B與X進行"或"運算，運算的結(jié)果擊中A。
[0057] 以下，將結(jié)合圖3、圖4和圖5對步驟S400的具體實現(xiàn)流程進行詳細描述。
[0058] 圖3示出根據(jù)本發(fā)明示例性實施例的對地質(zhì)異常體散點數(shù)據(jù)執(zhí)行先腐蝕后膨脹 的插值計算的步驟的流程圖。
[0059] 在步驟S411，將組成原始三維地震數(shù)據(jù)的所有數(shù)據(jù)點存入第一數(shù)組。
[0060] 在步驟S412,將第一數(shù)組中的組成地質(zhì)異常體散點數(shù)據(jù)的數(shù)據(jù)點賦予值1，將第 一數(shù)組中的其他數(shù)據(jù)點賦予值0。
[0061] 在步驟S413,按照公式（1)所示來執(zhí)行腐蝕插值。具體說來，用三維結(jié)構(gòu)元素來 掃描第一數(shù)組中的所有數(shù)據(jù)點，即，將三維結(jié)構(gòu)元素遍歷地與第一數(shù)組中的數(shù)據(jù)點被賦予 的值進行"與"運算，其中，三維結(jié)構(gòu)元素的中心元素的元素值為〇,周圍元素的元素值為1。 運算的結(jié)果作為第一數(shù)組中相應(yīng)數(shù)據(jù)點的更新值。
[0062] 在步驟S414,按照公式（2)所示來對步驟S413中得到的更新值執(zhí)行膨脹插值。具 體說來，用三維結(jié)構(gòu)元素來掃描步驟S413中更新后的第一數(shù)組中的所有數(shù)據(jù)點，S卩，將三 維結(jié)構(gòu)元素遍歷地與更新后的第一數(shù)組中的數(shù)據(jù)點被賦予的值進行"或"運算，其中，三維 結(jié)構(gòu)元素的中心元素的元素值為1，周圍元素的元素值為0。運算的結(jié)果作為第一數(shù)組中相 應(yīng)數(shù)據(jù)點的更新值。
[0063] 在步驟S415,將在步驟S414中得到的第一數(shù)組中具有值1的數(shù)據(jù)點確定為插值后 的地質(zhì)異常體散點數(shù)據(jù)。
[0064] 如圖3的流程圖所示，對地質(zhì)異常體散點數(shù)據(jù)執(zhí)行了先腐蝕后膨脹的插值計算。 可以看出，上述處理在不明顯改變地質(zhì)異常體的體積的情況下先消除了地質(zhì)異常體中的無 意義的數(shù)據(jù)點，并隨后進行了相應(yīng)的平滑處理，從而確保了地質(zhì)異常體散點數(shù)據(jù)的插值效 果。
[0065] 圖4示出根據(jù)本發(fā)明示例性實施例的對地質(zhì)異常體散點數(shù)據(jù)執(zhí)行先膨脹后腐蝕 的插值計算的步驟的流程圖。
[0066] 在步驟S421，將組成原始三維地震數(shù)據(jù)的所有數(shù)據(jù)點存入第一數(shù)組。
[0067] 在步驟S422,將第一數(shù)組中的組成地質(zhì)異常體散點數(shù)據(jù)的數(shù)據(jù)點賦予值1，將第 一數(shù)組中的其他數(shù)據(jù)點賦予值0。
[0068] 在步驟S423,按照公式（2)所示來執(zhí)行膨脹插值。具體說來，用三維結(jié)構(gòu)元素來 掃描第一數(shù)組中的所有數(shù)據(jù)點，即，將三維結(jié)構(gòu)元素遍歷地與第一數(shù)組中的數(shù)據(jù)點被賦予 的值進行"或"運算，其中，三維結(jié)構(gòu)元素的中心元素的元素值為1，周圍元素的元素值為0。 運算的結(jié)果作為第一數(shù)組中相應(yīng)數(shù)據(jù)點的更新值。
[0069] 在步驟S424,按照公式（1)所示來對步驟S423中得到的更新值執(zhí)行腐蝕插值。具 體說來，用三維結(jié)構(gòu)元素來掃描步驟S423中更新后的第一數(shù)組中的所有數(shù)據(jù)點，S卩，將三 維結(jié)構(gòu)元素遍歷地與更新后的第一數(shù)組中的數(shù)據(jù)點被賦予的值進行"與"運算，其中，三維 結(jié)構(gòu)元素的中心元素的元素值為0,周圍元素的元素值為1。運算的結(jié)果作為第一數(shù)組中相 應(yīng)數(shù)據(jù)點的更新值。
[0070] 在步驟S425,將在步驟S424中得到第一數(shù)組中具有值為1的數(shù)據(jù)點確定為插值后 的地質(zhì)異常體散點數(shù)據(jù)。
[0071] 如圖4的流程圖所示，對地質(zhì)異常體散點數(shù)據(jù)執(zhí)行了先膨脹后腐蝕的插值計算。 可以看出，上述處理在不明顯改變地質(zhì)異常體的體積的情況下先平滑地質(zhì)異常體的邊界， 而后消除地質(zhì)異常體中的無意義的數(shù)據(jù)點。
[0072] 可以看出，通過利用二值形態(tài)學(xué)中的膨脹和腐蝕對地質(zhì)異常體散點數(shù)據(jù)進行插值 處理，可達到有效地填補地質(zhì)異常體散點數(shù)據(jù)中的空洞、平滑地質(zhì)異常體散點數(shù)據(jù)的邊緣 的目的，經(jīng)過插值處理后的插值后的地質(zhì)異常體散點數(shù)據(jù)更加飽滿、緊湊且完整。這樣的插 值處理為后續(xù)的成面操作提供有效的成面數(shù)據(jù)。
[0073] 圖5示出根據(jù)本發(fā)明示例性實施例的基于上下表面曲面插值對地質(zhì)異常體散點 數(shù)據(jù)執(zhí)行插值計算的步驟的流程圖。
[0074] 具體地講，原始三維地震數(shù)據(jù)中的每個數(shù)據(jù)點的坐標信息可包括X軸坐標值、y軸 坐標值和z軸坐標值。在步驟S431，將地質(zhì)異常體散點數(shù)據(jù)中的X軸坐標值、y軸坐標值分 別相同的數(shù)據(jù)點作為同一道數(shù)據(jù)。
[0075] 在步驟S432,對至少具有兩個數(shù)據(jù)點的每一道數(shù)據(jù)中的z軸坐標值最大的數(shù)據(jù)點 形成的上表面進行曲面插值，使得每一道數(shù)據(jù)都存在z軸坐標值最大的數(shù)據(jù)點。這里，所插 入的數(shù)據(jù)點可以為原始三維地震數(shù)據(jù)中的數(shù)據(jù)點。具體地講，基于至少具有兩個數(shù)據(jù)點的 每一道數(shù)據(jù)中z軸坐標值最大的數(shù)據(jù)點進行插值以得到所有道數(shù)據(jù)中每一道數(shù)據(jù)z軸坐標 值最大的數(shù)據(jù)點。
[0076] 在步驟S433,對至少具有兩個數(shù)據(jù)點的每一道數(shù)據(jù)中的z軸坐標值最小的數(shù)據(jù)點 形成的下表面進行曲面插值，使得每一道數(shù)據(jù)都存在z軸坐標值最小的數(shù)據(jù)點。這里，所插 入的數(shù)據(jù)點可以為原始三維地震數(shù)據(jù)中的數(shù)據(jù)點。具體地講，基于至少具有兩個數(shù)據(jù)點的 每一道數(shù)據(jù)中z軸坐標值最小的數(shù)據(jù)點進行插值以得到所有道數(shù)據(jù)中每一道數(shù)據(jù)z軸坐標 值最小的數(shù)據(jù)點。
[0077] 在步驟S434,將經(jīng)過步驟S432和S433中的插值計算之后得到的所有數(shù)據(jù)點作為 插值后的地質(zhì)異常體散點數(shù)據(jù)。
[0078] 應(yīng)該理解，上述步驟S432和S433的執(zhí)行順序可互換，即先執(zhí)行步驟S433,再執(zhí)行 步驟S432,以得到插值后的地質(zhì)異常體散點數(shù)據(jù)。
[0079] 可以看出，通過上下表面曲面插值可得到異常體的外表面（即，上表面和下表面） 上的數(shù)據(jù)（即，插值后的地質(zhì)異常體散點數(shù)據(jù)）。這樣的插值處理為后續(xù)的成面操作提供有 效的成面數(shù)據(jù)。
[0080] 以下，將結(jié)合圖6對步驟S500的具體實現(xiàn)流程進行詳細描述。圖6示出根據(jù)本發(fā) 明示例性實施例的成面操作的步驟的流程圖。
[0081] 作為示例，可利用Marching cubes算法對插值后的地質(zhì)異常體散點數(shù)據(jù)進行成面 操作，其基本思路是：在插值后的地質(zhì)異常體散點數(shù)據(jù)所在的空間網(wǎng)格（體素立方體，該體 素立方體的棱長被設(shè)定為原始三維地震數(shù)據(jù)采樣率的正整數(shù)倍）中確定出與地質(zhì)異常體 的外表面（即，插值后的地質(zhì)異常體散點數(shù)據(jù)中位于最外層的數(shù)據(jù)點的集合所形成的面） 相交的體素立方體，對于每一個確定的相交體素立方體，插值計算出體素立方體的棱邊與 地質(zhì)異常體的外表面的交點。然后，根據(jù)體素立方體每一頂點與地質(zhì)異常體的外表面的相 對位置，將地質(zhì)異常體的外表面與體素立方體邊的交點按一定方式連接生成地質(zhì)異常體成 面數(shù)據(jù)模型的一部分，作為地質(zhì)異常體的外表面在體素立方體內(nèi)的逼近表不。
[0082] 具體說來，在步驟S501，設(shè)定預(yù)定閾值，所述預(yù)定閾值介于插值后的地質(zhì)異常體散 點數(shù)據(jù)的各個數(shù)據(jù)點的屬性值的最大值和最小值之間，且與每一個預(yù)定閾值對應(yīng)的所有數(shù) 據(jù)點將形成一個等值面。
[0083] 在步驟S502,采用以下公式計算得到與所述預(yù)定閾值對應(yīng)的所有數(shù)據(jù)點的坐標， 所述所有數(shù)據(jù)點中的任意一個數(shù)據(jù)點P的坐標為：
[0084] P = Pl+(isovalue-Vl) (P2-P1)/(V2-V1) (5)
[0085] 其中，isovalue為預(yù)定閾值，P2為插值后的地質(zhì)異常體散點數(shù)據(jù)中位于最外層的 數(shù)據(jù)點的坐標，且P2位于體素立方體的頂點上，Pl為原始三維地震數(shù)據(jù)中與P2的空間距離 最近的數(shù)據(jù)點的坐標，且位于與P2相同和不同的體素立方體的頂點上，Vl為坐標為Pl的 數(shù)據(jù)點的屬性值，V2為坐標為P2的數(shù)據(jù)點的屬性值。
[0086] 在步驟S503,根據(jù)用戶的需求判斷是否需要通過另一個預(yù)定閾值求出另一個等值 面。
[0087] 如果需要通過另一個預(yù)定閾值求出另一個等值面，則在步驟S504,設(shè)定另一預(yù)定 閾值，并返回步驟S502,直到獲得與期望的所有等值面對應(yīng)的所有對應(yīng)的數(shù)據(jù)點（即，在步 驟S503判斷不需要通過另一個預(yù)定閾值求出另一個等值面）時，結(jié)束成面操作。相應(yīng)地， 可將異常體成面數(shù)據(jù)模型構(gòu)建為所述所有對應(yīng)的數(shù)據(jù)點和插值后的地質(zhì)異常體散點數(shù)據(jù) 中位于最外層的所有數(shù)據(jù)點所形成的面。具體地講，作為示例，將上述數(shù)據(jù)點中每三個相鄰 的數(shù)據(jù)點形成一個三角面，所有三角面將形成最終的異常體成面數(shù)據(jù)模型。
[0088] 可以看出，通過上述成面操作可對地質(zhì)異常體的外表面進行近似描述。作為示例， 在步驟S600中，可將不同預(yù)定閾值（也即，屬性值）對應(yīng)的曲面采用不同的顏色進行顯示， 以實現(xiàn)地質(zhì)異常體的可視化。
[0089] 以上的描述主要是針對地質(zhì)領(lǐng)域的地質(zhì)異常體的可視化顯示，但本領(lǐng)域技術(shù)人員 應(yīng)可以理解，上述方法不限于僅在地質(zhì)領(lǐng)域可視化地顯示地質(zhì)異常體，也可以應(yīng)用于其他 領(lǐng)域（例如醫(yī)療衛(wèi)生領(lǐng)域等）的異常體的可視化顯示。以下，將結(jié)合圖7對可視化地顯示 異常體的方法進行詳細描述。圖7是根據(jù)本發(fā)明示例性實施例的可視化地顯示異常體的方 法的流程圖。
[0090] 在步驟S10,獲取原始三維圖像數(shù)據(jù)。作為示例，原始三維圖像數(shù)據(jù)可包括多個數(shù) 據(jù)點，每個數(shù)據(jù)點可具有各自的坐標信息和屬性值。可選地，所述屬性值可為能夠表明數(shù)據(jù) 點的某種特定屬性（例如，圖像三維屬性等）的值。
[0091] 在步驟S20,確定原始三維圖像數(shù)據(jù)中的有效異常體數(shù)據(jù)。作為示例，可根據(jù)針對 三維圖像數(shù)據(jù)的屬性值范圍參數(shù)、坐標值范圍參數(shù)和采樣間隔參數(shù)中的至少一個來確定原 始三維圖像數(shù)據(jù)中的有效異常體數(shù)據(jù)。具體說來，可預(yù)先獲取三維圖像數(shù)據(jù)的參數(shù)列表，其 中，所述參數(shù)列表包括針對三維圖像數(shù)據(jù)的屬性值范圍參數(shù)、坐標值范圍參數(shù)和采樣間隔 參數(shù)。這里，三維圖像數(shù)據(jù)的參數(shù)列表中的參數(shù)可由用戶設(shè)定。
[0092] 例如，可將原始圖像數(shù)據(jù)中滿足屬性值范圍參數(shù)的數(shù)據(jù)點確定為有效異常體數(shù) 據(jù)。或者，可將原始圖像數(shù)據(jù)中滿足坐標值范圍參數(shù)的數(shù)據(jù)點確定為有效異常體數(shù)據(jù)?；?者，可將對原始圖像數(shù)據(jù)根據(jù)采樣間隔參數(shù)進行采樣后得到對數(shù)據(jù)點確定為有效異常體數(shù) 據(jù)?；蛘?，可將上述條件進行組合以確定有效異常體數(shù)據(jù)。
[0093] 在步驟S30,利用三維圖像種子點生長方法確定有效異常體數(shù)據(jù)中的異常體散點 數(shù)據(jù)。
[0094] 作為示例，可將有效異常體數(shù)據(jù)中的數(shù)據(jù)點作為種子點，利用泛洪法基于所述種 子點進行三維圖像種子點生長，將經(jīng)過三維圖像種子點生長得到的所有數(shù)據(jù)點作為異常體 散點數(shù)據(jù)。
[0095] 進一步地，作為示例，可利用六鄰域（如圖2A所示的立方體的六個面的中心點） 或二十六鄰域（如圖2B所示的立方體的六個面的中心點、八個頂點和十二條棱的中點）的 泛洪法基于所述種子點進行三維圖像種子點生長。
[0096] 在步驟S40,對異常體散點數(shù)據(jù)進行插值計算，以得到插值后的異常體散點數(shù)據(jù)。
[0097] 作為示例，可基于二值形態(tài)學(xué)對異常體散點數(shù)據(jù)執(zhí)行先膨脹后腐蝕的插值計算。 或者，可基于二值形態(tài)學(xué)對異常體散點數(shù)據(jù)執(zhí)行先腐蝕后膨脹的插值計算?；蛘撸苫谏?下表面曲面插值對異常體散點數(shù)據(jù)進行插值計算。
[0098] 在步驟S50,對插值后的異常體散點數(shù)據(jù)進行成面操作，以得到異常體成面數(shù)據(jù)模 型。
[0099] 作為示例,可利用Marching cubes算法對插值后的異常體散點數(shù)據(jù)進行成面操 作，以得到異常體成面數(shù)據(jù)模型。
[0100] 在步驟S60,基于異常體成面數(shù)據(jù)模型來可視化地顯示異常體。
[0101] 進一步地，作為示例，在顯示異常體的頁面中可提供異常體的體積，異常體的體積 可基于異常體成面數(shù)據(jù)模型而計算得到。例如，可通過將異常體劃分為多個小立方體并將 小立方體的體積求和得到異常體的體積，或者將異常體分割為不同空間區(qū)域并將每個空間 區(qū)域進行積分運算的結(jié)果求和以得到異常體的體積。這里，可以看出，通過異常體成面模型 可有效地估算異常體的體積，為后續(xù)對異常體進行其他的處理提供便利。
[0102] 可以看出，可視化地顯示異常體的方法不僅適用于地質(zhì)領(lǐng)域，也可將其應(yīng)用到其 他領(lǐng)域（例如醫(yī)療衛(wèi)生領(lǐng)域等），對其他領(lǐng)域中的任意目標體作為上述步驟中的異常體實 施上述步驟，即，可實現(xiàn)三維目標體的數(shù)據(jù)提取和可視化。
[0103] 以下，將結(jié)合圖8來描述根據(jù)本發(fā)明示例性實施例的可視化地顯示地質(zhì)異常體的 裝置。所述裝置可實現(xiàn)圖1到圖6所示的方法。
[0104] 圖8是根據(jù)本發(fā)明示例性實施例的可視化地顯示地質(zhì)異常體的裝置的框圖。
[0105] 如圖8所示，根據(jù)本發(fā)明示例性實施例的可視化地顯示地質(zhì)異常體的裝置可包 括：獲取單元100、有效數(shù)據(jù)確定單元200、散點確定單元300、插值單元400、成面單元500 和顯示單元600。這些單元可由數(shù)字信號處理器、現(xiàn)場可編程門陣列等通用硬件處理器來實 現(xiàn)，也可通過專用芯片等專用硬件處理器來實現(xiàn)，還可完全通過計算機程序來以軟件方式 實現(xiàn)，例如，被實現(xiàn)為安裝在電子設(shè)備（例如，計算機）中用于可視化地顯示地質(zhì)異常體的 各個模塊。
[0106] 具體地講，獲取單元100用于獲取原始三維地震數(shù)據(jù)。作為示例，原始三維地震數(shù) 據(jù)可包括多個數(shù)據(jù)點，每個數(shù)據(jù)點可具有各自的坐標信息和屬性值。可選地，所述屬性值可 為能夠表明數(shù)據(jù)點的某種特定屬性的值，例如地震屬性值等。
[0107] 有效數(shù)據(jù)確定單元200用于確定原始三維地震數(shù)據(jù)中的有效地質(zhì)異常體數(shù)據(jù)。作 為示例，可根據(jù)針對三維地震數(shù)據(jù)的屬性值范圍參數(shù)、坐標值范圍參數(shù)和采樣間隔參數(shù)中 的至少一個來確定原始三維地震數(shù)據(jù)中的有效地質(zhì)異常體數(shù)據(jù)。具體說來，可預(yù)先獲取三 維地震數(shù)據(jù)的參數(shù)列表，其中，所述參數(shù)列表包括針對三維地震數(shù)據(jù)的屬性值范圍參數(shù)、坐 標值范圍參數(shù)和采樣間隔參數(shù)。這里，三維地震數(shù)據(jù)的參數(shù)列表中的參數(shù)可由用戶設(shè)定。
[0108] 例如，有效數(shù)據(jù)確定單元200可將原始地震數(shù)據(jù)中滿足屬性值范圍參數(shù)的數(shù)據(jù)點 確定為有效地質(zhì)異常體數(shù)據(jù)。或者，有效數(shù)據(jù)確定單元200可將原始地震數(shù)據(jù)中滿足坐標 值范圍參數(shù)的數(shù)據(jù)點確定為有效地質(zhì)異常體數(shù)據(jù)?；蛘撸行?shù)據(jù)確定單元200可將對原 始地震數(shù)據(jù)根據(jù)采樣間隔參數(shù)進行采樣后得到的數(shù)據(jù)點確定為有效地質(zhì)異常體數(shù)據(jù)。或 者，有效數(shù)據(jù)確定單元200可將上述條件進行組合以確定有效地質(zhì)異常體數(shù)據(jù)。
[0109] 散點確定單元300用于利用三維圖像種子點生長方法確定有效地質(zhì)異常體數(shù)據(jù) 中的地質(zhì)異常體散點數(shù)據(jù)。
[0110] 作為示例，可將有效地質(zhì)異常體數(shù)據(jù)中的數(shù)據(jù)點作為種子點，利用泛洪法基于所 述種子點進行三維圖像種子點生長，將經(jīng)過三維圖像種子點生長得到的所有數(shù)據(jù)點作為地 質(zhì)異常體散點數(shù)據(jù)。
[0111] 進一步地，作為示例，可利用六鄰域（如圖2A所示的立方體的六個面的中心點） 或二十六鄰域（如圖2B所示的立方體的六個面的中心點、八個頂點和十二條棱的中點）的 泛洪法基于所述種子點進行三維圖像種子點生長。
[0112] 插值單元400用于對地質(zhì)異常體散點數(shù)據(jù)進行插值計算，以得到插值后的地質(zhì)異 常體散點數(shù)據(jù)。
[0113] 作為示例，可基于二值形態(tài)學(xué)對地質(zhì)異常體散點數(shù)據(jù)執(zhí)行先膨脹后腐蝕的插值計 算?；蛘?，可基于二值形態(tài)學(xué)對地質(zhì)異常體散點數(shù)據(jù)執(zhí)行先腐蝕后膨脹的插值計算?；蛘?， 可基于上下表面曲面插值對地質(zhì)異常體散點數(shù)據(jù)進行插值計算。這里，關(guān)于二值形態(tài)學(xué)的 相關(guān)內(nèi)容已經(jīng)在對圖1的步驟S400的具體描述時進行了詳細的闡述，此處為避免重復(fù)，不 再贅述。
[0114] 成面單元500用于對插值后的地質(zhì)異常體散點數(shù)據(jù)進行成面操作，以得到地質(zhì)異 常體成面數(shù)據(jù)模型。
[0115] 作為示例，可利用Marching cubes算法對插值后的地質(zhì)異常體散點數(shù)據(jù)進行成面 操作，以得到地質(zhì)異常體成面數(shù)據(jù)模型。這里，地質(zhì)異常體成面數(shù)據(jù)模型可由多個等值面組 成，每個等值面對應(yīng)一個預(yù)定閾值。
[0116] 顯示單元600用于基于地質(zhì)異常體成面數(shù)據(jù)模型來可視化地顯示地質(zhì)異常體。
[0117] 進一步地，作為示例，在顯示地質(zhì)異常體的頁面中可提供地質(zhì)異常體的體積，地質(zhì) 異常體的體積可基于地質(zhì)異常體成面數(shù)據(jù)模型而計算得到。例如，可通過將地質(zhì)異常體劃 分為多個小立方體并將小立方體的體積求和得到地質(zhì)異常體的體積，或者將地質(zhì)異常體分 割為不同空間區(qū)域并將每個空間區(qū)域進行積分運算的結(jié)果求和以得到地質(zhì)異常體的體積。 這里，可以看出，通過地質(zhì)異常體成面數(shù)據(jù)模型可有效地估算地質(zhì)異常體的體積，為后續(xù)對 地質(zhì)異常體進行其他的工作提供便利。
[0118] 由此可見，通過上述裝置能夠高效地從原始三維地震數(shù)據(jù)中提取希望得到的地質(zhì) 異常體散點數(shù)據(jù)，并實現(xiàn)了對地質(zhì)異常體的可視化顯示。
[0119] 上述可視化地顯示地質(zhì)異常體的方法和裝置，能夠高效地從原始三維地震數(shù)據(jù)中 提取希望得到的地質(zhì)異常體散點數(shù)據(jù)，并實現(xiàn)了對地質(zhì)異常體的可視化顯示。
[0120] 上面已經(jīng)結(jié)合具體示例性實施例描述了本發(fā)明，但是本發(fā)明的實施不限于此。在 本發(fā)明的精神和范圍內(nèi)，本領(lǐng)域技術(shù)人員可以進行各種修改和變型，這些修改和變型將落 入權(quán)利要求限定的保護范圍之內(nèi)。
【權(quán)利要求】
1. 一種可視化地顯示地質(zhì)異常體的方法，其特征在于，包括以下步驟： (a) 獲取原始三維地震數(shù)據(jù)； (b) 確定原始三維地震數(shù)據(jù)中的有效地質(zhì)異常體數(shù)據(jù)； (c) 利用三維圖像種子點生長方法確定有效地質(zhì)異常體數(shù)據(jù)中的地質(zhì)異常體散點數(shù) 據(jù)； (d) 對所述地質(zhì)異常體散點數(shù)據(jù)進行插值計算，以得到插值后的地質(zhì)異常體散點數(shù) 據(jù)； (e) 對插值后的地質(zhì)異常體散點數(shù)據(jù)進行成面操作，以得到地質(zhì)異常體成面數(shù)據(jù)模 型； (f) 基于地質(zhì)異常體成面數(shù)據(jù)模型來可視化地顯示地質(zhì)異常體。
2. 如權(quán)利要求1所述的方法，其特征在于，步驟（b)包括：根據(jù)針對三維地震數(shù)據(jù)的屬 性值范圍參數(shù)、坐標值范圍參數(shù)和采樣間隔參數(shù)中的至少一個確定原始三維地震數(shù)據(jù)中的 有效地質(zhì)異常體數(shù)據(jù)。
3. 如權(quán)利要求1所述的方法，其特征在于，步驟（c)包括：將有效地質(zhì)異常體數(shù)據(jù)中的 數(shù)據(jù)點作為種子點，利用泛洪法基于所述種子點進行三維圖像種子點生長，將經(jīng)過三維圖 像種子點生長得到的所有數(shù)據(jù)點作為地質(zhì)異常體散點數(shù)據(jù)。
4. 如權(quán)利要求3所述的方法，其特征在于，在步驟（c)中，利用六鄰域或二十六鄰域的 泛洪法基于所述種子點進行三維圖像種子點生長。
5. 如權(quán)利要求1所述的方法，其特征在于，在步驟（d)中，基于二值形態(tài)學(xué)對所述地質(zhì) 異常體散點數(shù)據(jù)執(zhí)行先膨脹后腐蝕的插值計算，或者，基于二值形態(tài)學(xué)對所述地質(zhì)異常體 散點數(shù)據(jù)執(zhí)行先腐蝕后膨脹的插值計算。
6. 如權(quán)利要求1所述的方法，其特征在于，原始三維地震數(shù)據(jù)中的每個數(shù)據(jù)點具有各 自的坐標信息，所述坐標信息包括x軸坐標值、y軸坐標值和z軸坐標值，并且，在步驟（d) 中，將所述地質(zhì)異常體散點數(shù)據(jù)中的x軸坐標值、y軸坐標值分別相同的數(shù)據(jù)點作為同一道 數(shù)據(jù)，基于至少具有兩個數(shù)據(jù)點的每一道數(shù)據(jù)中z軸坐標值最大的數(shù)據(jù)點進行插值以針對 每道數(shù)據(jù)得到z軸坐標值最大的數(shù)據(jù)點，并基于至少具有兩個數(shù)據(jù)點的每一道數(shù)據(jù)中z軸 坐標值最小的數(shù)據(jù)點進行插值以針對每一道數(shù)據(jù)得到z軸坐標值最小的數(shù)據(jù)點。
7. 如權(quán)利要求1所述的方法，其特征在于，步驟（e)包括：利用Marchingcubes算法對 插值后的地質(zhì)異常體散點數(shù)據(jù)進行成面操作，以得到地質(zhì)異常體成面數(shù)據(jù)模型。
8. 如權(quán)利要求1所述的方法，其特征在于，步驟（f)還包括：在顯示地質(zhì)異常體的頁面 中提供地質(zhì)異常體的體積，其中，所述地質(zhì)異常體的體積基于地質(zhì)異常體成面數(shù)據(jù)模型而 計算得到。
9. 一種可視化地顯示異常體的方法，其特征在于，包括以下步驟： (a) 獲取原始三維圖像數(shù)據(jù)； (b) 確定原始三維圖像數(shù)據(jù)中的有效異常體數(shù)據(jù)； (c) 利用三維圖像種子點生長方法確定有效異常體數(shù)據(jù)中的異常體散點數(shù)據(jù)； (d) 對所述異常體散點數(shù)據(jù)進行插值計算，以得到插值后的異常體散點數(shù)據(jù)； (e) 對插值后的異常體散點數(shù)據(jù)進行成面操作，以得到異常體成面數(shù)據(jù)模型； (f) 基于異常體成面數(shù)據(jù)模型來可視化地顯示異常體。
10. -種可視化地顯示地質(zhì)異常體的裝置，包括： 獲取單元，用于獲取原始三維地震數(shù)據(jù)； 有效數(shù)據(jù)確定單元，用于確定原始三維地震數(shù)據(jù)中的有效地質(zhì)異常體數(shù)據(jù)； 散點確定單元，用于利用三維圖像種子點生長方法確定有效地質(zhì)異常體數(shù)據(jù)中的地質(zhì) 異常體散點數(shù)據(jù)； 插值單元，用于對所述地質(zhì)異常體散點數(shù)據(jù)進行插值計算，以得到插值后的地質(zhì)異常 體散點數(shù)據(jù)； 成面單元，用于對插值后的地質(zhì)異常體散點數(shù)據(jù)進行成面操作，以得到地質(zhì)異常體成 面數(shù)據(jù)模型； 顯示單元，用于基于地質(zhì)異常體成面數(shù)據(jù)模型來可視化地顯示地質(zhì)異常體。
【文檔編號】G01V1/34GK104360392SQ201410719825
【公開日】2015年2月18日 申請日期:2014年12月1日 優(yōu)先權(quán)日:2014年12月1日
【發(fā)明者】王頎, 鄒文, 何光明, 陳愛萍, 陶正喜, 夏銘, 陳小二 申請人:中國石油集團川慶鉆探工程有限公司地球物理勘探公司