陳 雷 石家琦 石太昆 張 婷

(*北京信息科技大學(xué)計(jì)算機(jī)學(xué)院 北京 100101)

(**國(guó)勘數(shù)字地球(北京)科技有限公司 北京 100096)

(***北京工業(yè)大學(xué)信息學(xué)部 北京 100124)

0 引言

地震數(shù)據(jù)剖面可視化技術(shù)對(duì)于地震勘探中的解釋環(huán)節(jié)有著關(guān)鍵性的作用[1-2]。地震數(shù)據(jù)直觀、正確、特征層次分明的可視化可以極大地方便地質(zhì)科研人員對(duì)相關(guān)數(shù)據(jù)做出準(zhǔn)確的分析、得出正確的結(jié)論,對(duì)油氣勘探、地質(zhì)勘察等具有重要意義。

目前,地震數(shù)據(jù)可視化技術(shù)主要有2 類(lèi):一類(lèi)是二維可視化[3-6],另一類(lèi)是三維可視化[1-2,7-16]。地震數(shù)據(jù)二維可視化是將地震剖面數(shù)據(jù)轉(zhuǎn)換成一幅二維圖像進(jìn)行顯示的技術(shù),又可分為波形顯示、波形+變面積顯示、灰度顯示和彩色變密度顯示等4 個(gè)類(lèi)別。這種傳統(tǒng)可視化方法一定程度上展示了地震數(shù)據(jù)的空間變化特征,其不足是動(dòng)態(tài)范圍小、細(xì)節(jié)丟失嚴(yán)重、不能直觀顯示微觀地質(zhì)結(jié)構(gòu)、不便于地震解釋和地質(zhì)分析等。近年來(lái),隨著計(jì)算機(jī)圖形學(xué)的發(fā)展和顯卡中圖形處理單元(graphics processing unit,GPU)計(jì)算能力的提高,地震數(shù)據(jù)三維可視化獲得了廣泛研究,取得了較大進(jìn)展。它是將地震數(shù)據(jù)在三維空間中進(jìn)行可視化的技術(shù),是為解決二維可視化的不足而提出的。但現(xiàn)有地震數(shù)據(jù)三維可視化基本都是對(duì)三維地震數(shù)據(jù)體進(jìn)行可視化,再按需抽取某一地震剖面的灰度或彩色變密度顯示[1-2,7-16],這種可視化方法對(duì)地震微觀結(jié)構(gòu)的揭示能力并沒(méi)有實(shí)質(zhì)提高。文獻(xiàn)[15]提出了將二維地震剖面進(jìn)行三維可視化的方法,提高了對(duì)微觀結(jié)構(gòu)的顯示能力,但在頂點(diǎn)鑲嵌方面存在鑲嵌方式不能根據(jù)地震樣點(diǎn)傾斜方式自適應(yīng)選擇的問(wèn)題。

針對(duì)文獻(xiàn)[15]中可視化方法的網(wǎng)格頂點(diǎn)鑲嵌問(wèn)題,本文提出了一種新的低傾斜相關(guān)自適應(yīng)鑲嵌方法,以盡可能使網(wǎng)格頂點(diǎn)根據(jù)相應(yīng)地震樣點(diǎn)的傾斜方式來(lái)自適應(yīng)鑲嵌。

1 算法分析

本節(jié)先對(duì)文獻(xiàn)[15]中的地震剖面三維可視化方法進(jìn)行分析,然后再對(duì)所提鑲嵌方法進(jìn)行論證。

1.1 地震剖面數(shù)據(jù)三維可視化方法分析

文獻(xiàn)[15]提出的二維地震剖面三維可視化方法具有以下特征。

(1) 對(duì)二維地震剖面數(shù)據(jù)S(x,y),x表示道號(hào),y表示時(shí)間深度值,以(x,y) 處的地震樣點(diǎn)幅度S(x,y) 作為三維空間中的z坐標(biāo),即地震樣點(diǎn)對(duì)應(yīng)的三維空間中的點(diǎn)為P(x,y,S(x,y))。

(2) 除要顯示的最后一道地震數(shù)據(jù)外,將每一道上相鄰2 個(gè)地震樣點(diǎn)和其下一道上時(shí)間深度相同的相鄰2 個(gè)樣點(diǎn)連接為四邊形,由同一道上3 個(gè)地震樣點(diǎn)和下一道上1 個(gè)樣點(diǎn)連接成的2 個(gè)三角形單元組成,以此將地震剖面數(shù)據(jù)鑲嵌為網(wǎng)格曲面。

(3) 創(chuàng)建一幅由一種顏色漸變到另一種顏色的調(diào)色板紋理,將經(jīng)過(guò)光柵化插值之后的地震幅度值S(x,y) 通過(guò)線性函數(shù)映射為紋理坐標(biāo),在GPU 像素著色器中根據(jù)此紋理坐標(biāo)采樣調(diào)色板紋理為網(wǎng)格頂點(diǎn)和三角形單元內(nèi)部的像素點(diǎn)著色。

(4) 將相鄰道上的4 個(gè)地震樣點(diǎn)鑲嵌為四邊形時(shí),鑲嵌有2 種方式,即上傾斜鑲嵌和下傾斜鑲嵌;不正確的鑲嵌方式將產(chǎn)生錯(cuò)誤的結(jié)果(脊或者槽)。對(duì)此問(wèn)題,采用了低傾斜相關(guān)自適應(yīng)鑲嵌,在GPU像素著色器中計(jì)算以當(dāng)前樣點(diǎn)為中心的某一長(zhǎng)度的樣點(diǎn)序列與下一道上一定偏移范圍內(nèi)的同樣長(zhǎng)度序列的相關(guān)系數(shù),將最大相關(guān)系數(shù)所對(duì)應(yīng)的偏移值作為當(dāng)前樣點(diǎn)的傾斜方式,通過(guò)渲染到紋理將數(shù)值化的傾斜方式存儲(chǔ)到紋理中,用來(lái)定義網(wǎng)格頂點(diǎn)的索引值。

(5) 將坐標(biāo)變換、光照計(jì)算、紋理坐標(biāo)映射等內(nèi)存消耗較大的計(jì)算都轉(zhuǎn)移到GPU 中進(jìn)行,將地震剖面數(shù)據(jù)填充到一幅浮點(diǎn)紋理中,在GPU 中通過(guò)采樣浮點(diǎn)紋理實(shí)現(xiàn)對(duì)地震剖面數(shù)據(jù)的訪問(wèn)。

1.1.1 網(wǎng)格頂點(diǎn)著色及顏色插值

在固定渲染管線和可編程渲染管線中,頂點(diǎn)數(shù)據(jù)經(jīng)過(guò)初始裝配之后的處理如圖1 所示。

圖1 渲染管線中經(jīng)過(guò)初始裝配之后的頂點(diǎn)處理過(guò)程

初始裝配之后的頂點(diǎn)數(shù)據(jù)需進(jìn)行光柵化與插值。固定渲染管線中,光柵化與插值輸出的像素經(jīng)過(guò)光柵器后直接送往幀緩沖等待顯示;可編程渲染管線中,光柵化與插值輸出的像素先經(jīng)過(guò)可編程像素著色器處理,再送往光柵器直至顯示。在像素著色器中,可根據(jù)頂點(diǎn)著色器輸出的信息對(duì)頂點(diǎn)顏色、紋理坐標(biāo)、光照等進(jìn)行處理,然后渲染輸出。

地震剖面三維可視化中,對(duì)不同幅度的樣點(diǎn)所對(duì)應(yīng)的網(wǎng)格頂點(diǎn)及其之間的像素用不同的顏色進(jìn)行著色,可有效提高對(duì)地震事件的分辨能力。著色有2 種方式:一種是將顏色作為頂點(diǎn)屬性分量為頂點(diǎn)及其之間的像素著色,可稱為頂點(diǎn)屬性著色方式;另一種是將光柵化與插值輸出的像素傳送到GPU 可編程像素著色器中,在像素著色器中根據(jù)頂點(diǎn)著色器輸出的地震幅度對(duì)調(diào)色板紋理進(jìn)行采樣,在其中,先將線性插值后的地震樣點(diǎn)幅度線性映射為紋理坐標(biāo),然后再采樣調(diào)色板紋理,以實(shí)現(xiàn)網(wǎng)格頂點(diǎn)及頂點(diǎn)之間三角形單元內(nèi)部像素的正確著色,這種方式可稱為線性漸進(jìn)調(diào)色板著色方式。文獻(xiàn)[15]采用第2種著色方式有效避免了地震數(shù)據(jù)可視化中網(wǎng)格頂點(diǎn)之間像素的顏色交叉問(wèn)題。

1.1.2 網(wǎng)格頂點(diǎn)低傾斜相關(guān)自適應(yīng)鑲嵌

網(wǎng)格頂點(diǎn)鑲嵌是地震剖面數(shù)據(jù)三維可視化中著色之外的另一個(gè)關(guān)鍵問(wèn)題,即如何在三維空間中將地震樣點(diǎn)對(duì)應(yīng)的頂點(diǎn)連接成三角形單元,以形成網(wǎng)格曲面。文獻(xiàn)[15]所用方法是除要顯示的最后一道地震數(shù)據(jù)外,將每一道上相鄰2 個(gè)地震樣點(diǎn)和其下一道上時(shí)間深度相同的相鄰2 個(gè)樣點(diǎn)連接為四邊形(此四邊形由同一道上2 個(gè)地震樣點(diǎn)和下一道上1個(gè)樣點(diǎn)連接成的2 個(gè)三角形單元組成),以將地震剖面鑲嵌為網(wǎng)格曲面。對(duì)以上4 個(gè)頂點(diǎn),不同鑲嵌方式會(huì)產(chǎn)生不同的效果,如圖2 所示。

圖2 相同4 個(gè)頂點(diǎn)的2 種鑲嵌方式

采用圖2(a)所示的上傾斜鑲嵌方式會(huì)產(chǎn)生“脊”結(jié)構(gòu)如圖2(c)所示;采用圖2(b)所示的下傾斜鑲嵌方式會(huì)產(chǎn)生“槽”結(jié)構(gòu)如圖2(d)所示;若統(tǒng)一采用上傾斜鑲嵌,則對(duì)下傾斜地震事件不合適,會(huì)在事件頂部產(chǎn)生“鋸齒”現(xiàn)象,反之亦然。

為解決傾斜事件頂部的“鋸齒”問(wèn)題,文獻(xiàn)[15]提出了一種低傾斜相關(guān)自適應(yīng)鑲嵌方式,在GPU 像素著色器中計(jì)算以當(dāng)前地震樣點(diǎn)為中心的某一長(zhǎng)度的樣點(diǎn)序列與下一道上同樣長(zhǎng)度序列在某一滑動(dòng)窗口(從上向下滑動(dòng))內(nèi)的相關(guān)系數(shù),并找出最大相關(guān)系數(shù)對(duì)應(yīng)的滑動(dòng)值。如果此滑動(dòng)值大于等于0,則采用下傾斜鑲嵌,用-1 表示,即以當(dāng)前頂點(diǎn)為控制點(diǎn)(四邊形左上角的頂點(diǎn))的四邊形采用下傾斜鑲嵌;反之,如果此滑動(dòng)值小于0,則以當(dāng)前頂點(diǎn)為控制點(diǎn)的四邊形采用上傾斜方式鑲嵌,用+1 表示。計(jì)算出鑲嵌方式之后,再把數(shù)值化的鑲嵌方式通過(guò)渲染到紋理存儲(chǔ)到紋理元素中,用以填充網(wǎng)格頂點(diǎn)索引。低傾斜是指在求2 個(gè)樣點(diǎn)序列的互相關(guān)函數(shù)值時(shí),只允許滑動(dòng)窗口有上下一個(gè)樣點(diǎn)間隔的滑動(dòng)范圍,這樣作為四邊形控制點(diǎn)的當(dāng)前道上的地震樣點(diǎn)只能與下一道上上傾斜一個(gè)采樣間隔或下傾斜一個(gè)采樣間隔的樣點(diǎn)相連。沒(méi)有采用高傾斜相關(guān)自適應(yīng)鑲嵌方式的原因是高傾斜情況下,鑲嵌方式過(guò)于復(fù)雜,很難處理自適應(yīng)鑲嵌帶來(lái)的網(wǎng)格漏洞和三角形單元重疊問(wèn)題。

文獻(xiàn)[15]所提低傾斜相關(guān)自適應(yīng)鑲嵌方式的不足是:只有當(dāng)四邊形控制點(diǎn)的傾斜方式為0(此時(shí)采用下傾斜鑲嵌)或-1(下傾斜)時(shí),組成此四邊形的2 個(gè)三角形單元的連接方式才與控制點(diǎn)的傾斜方式一致;如果控制點(diǎn)的傾斜方式為1,那么此控制點(diǎn)所控制的4 個(gè)頂點(diǎn)都要按照同一地震道上控制點(diǎn)之下一個(gè)樣點(diǎn)的上傾斜方式來(lái)鑲嵌,無(wú)論下一個(gè)樣點(diǎn)的傾斜方式是上傾斜、不傾斜還是下傾斜。這樣很大一部分網(wǎng)格頂點(diǎn)的鑲嵌方式不能與地震樣點(diǎn)的傾斜方式相一致。如圖3 所示。

圖3 低傾斜相關(guān)自適應(yīng)鑲嵌方式的不足說(shuō)明

文獻(xiàn)[15]中的方法以四邊形為基本單元鑲嵌地震數(shù)據(jù)網(wǎng)格,如圖3 所示,以四邊形ABCD為例,A點(diǎn)為四邊形的控制點(diǎn),四邊形的鑲嵌方式由A點(diǎn)的傾斜方式?jīng)Q定。當(dāng)A為下傾斜(或不傾斜)時(shí)采用圖3(a)所示的下傾斜鑲嵌方式鑲嵌四邊形ABCD,此時(shí)鑲嵌方式與A點(diǎn)的傾斜方式一致;當(dāng)A為上傾斜時(shí)采用圖3(b)所示的上傾斜方式鑲嵌四邊形ABCD,而此時(shí)按相關(guān)系數(shù)最大化原則,A點(diǎn)與N點(diǎn)相關(guān)度最大,A點(diǎn)應(yīng)與N點(diǎn)相連,生成四邊形MACN。而文獻(xiàn)[15]中方法生成了四邊形ABCD,這等于忽略了A點(diǎn)的上傾斜方式,同時(shí)將B點(diǎn)的傾斜方式強(qiáng)制設(shè)為上傾斜,當(dāng)B點(diǎn)的傾斜方式不是上傾斜時(shí),相當(dāng)于A、B2 個(gè)點(diǎn)都沒(méi)有根據(jù)其傾斜方式進(jìn)行鑲嵌,故文獻(xiàn)[15]所提低傾斜相關(guān)自適應(yīng)鑲嵌存在鑲嵌錯(cuò)誤問(wèn)題,需進(jìn)一步改進(jìn)。

1.2 本文算法

為解決1.1.2 小節(jié)中低傾斜相關(guān)自適應(yīng)鑲嵌方式存在的問(wèn)題,盡可能使得每一個(gè)頂點(diǎn)都按照相應(yīng)樣點(diǎn)的傾斜方式來(lái)鑲嵌(對(duì)于邊界頂點(diǎn)和產(chǎn)生漏洞的區(qū)域,為了保證網(wǎng)格的完整性與連續(xù)性,頂點(diǎn)鑲嵌方式可能與地震樣點(diǎn)的傾斜方式不一致),避免文獻(xiàn)[15]中低傾斜相關(guān)自適應(yīng)鑲嵌方法的不足。本文提出了新的低傾斜相關(guān)自適應(yīng)鑲嵌方法,具體如下。

首先在GPU 像素著色器中計(jì)算相鄰地震道上2個(gè)樣點(diǎn)序列的互相關(guān)函數(shù)值,以求取地震剖面中各地震樣點(diǎn)的傾斜方式及相應(yīng)的鑲嵌方式。GPU 像素著色器中2 個(gè)樣點(diǎn)序列的互相關(guān)函數(shù)值的計(jì)算原理,如圖4 所示。

圖4 計(jì)算互相關(guān)函數(shù)值原理圖

互相關(guān)函數(shù)和相關(guān)系數(shù)的計(jì)算公式分別為

其中,式(1)是計(jì)算序列X和Y的互相關(guān)函數(shù)值RXY(τ) 的公式,下標(biāo)XY表示2 個(gè)序列,X、Y的順序表示序列X未移動(dòng),序列Y移動(dòng)τ個(gè)單位;式(2)中表示序列X和Y的相關(guān)系數(shù)大小。τ是移位參數(shù),T為序列周期,N表示序列X和Y中元素編號(hào)的最大值,k表示元素序號(hào),k∈[0,N]。計(jì)算出相關(guān)系數(shù)后,找出最大值,根據(jù)其對(duì)應(yīng)的τ值來(lái)判斷當(dāng)前地震樣點(diǎn)是上傾斜還是下傾斜或者不傾斜,上傾斜用1 表示,下傾斜用-1 表示,傾斜幅度較小時(shí)視為不傾斜,用0 表示。也就是說(shuō)將最大相關(guān)系數(shù)對(duì)應(yīng)的τ值取反就可得地震樣點(diǎn)的傾斜方式。然后,使用渲染到紋理這一操作,把代表傾斜方式的數(shù)值存儲(chǔ)到紋理中,因紋理元素有R、G、B、A等多個(gè)分量,這里只用R分量來(lái)存儲(chǔ)傾斜方式數(shù)值,其他分量均為0。

獲得每個(gè)地震樣點(diǎn)的傾斜方式之后,根據(jù)傾斜方式的不同對(duì)地震樣點(diǎn)進(jìn)行自適應(yīng)鑲嵌,鑲嵌時(shí)分為3 種情況,具體如下。

(1) 頂部邊緣樣點(diǎn)均采用下傾斜鑲嵌,如圖5所示。圖5(a)為上邊緣點(diǎn)鑲嵌示意圖,圖5(b)為下邊緣點(diǎn)鑲嵌示意圖,T和B為當(dāng)前要鑲嵌的地震樣點(diǎn),它們分別處于地震剖面的上邊緣(時(shí)間深度值最小)和下邊緣(時(shí)間深度值最大)。對(duì)邊緣點(diǎn)采用統(tǒng)一的鑲嵌方式是為了避免當(dāng)同一地震道上頂部2 個(gè)或者底部2 個(gè)地震樣點(diǎn)的傾斜方式相反時(shí),由于其之間的矩形區(qū)域內(nèi)沒(méi)有鑲嵌三角形而形成漏洞,接下來(lái)就可以對(duì)其余非上下邊緣點(diǎn)進(jìn)行自適應(yīng)鑲嵌。

圖5 邊緣點(diǎn)鑲嵌示意圖

(2) 當(dāng)非上下邊緣點(diǎn)的傾斜方式為上傾斜時(shí),為避免同一矩形內(nèi)的4 個(gè)樣點(diǎn)被重復(fù)鑲嵌而產(chǎn)生三角形交叉重疊現(xiàn)象,先查看其上面相鄰地震樣點(diǎn)的傾斜方式是否為下傾斜,如果不是,則進(jìn)行上傾斜鑲嵌,如果是,則不進(jìn)行鑲嵌,如圖6 所示。

圖6 非上下邊緣點(diǎn)上傾斜鑲嵌示意圖

圖6 中,B點(diǎn)為當(dāng)前樣點(diǎn),圖6(a)為當(dāng)前一個(gè)樣點(diǎn)A的傾斜方式為下傾斜時(shí)產(chǎn)生的交叉現(xiàn)象。此時(shí)不對(duì)B點(diǎn)進(jìn)行鑲嵌,即不給B、D和C分配頂點(diǎn)索引以形成三角形BDC;圖6(b)為A點(diǎn)為上傾斜時(shí)的情況。

(3) 當(dāng)非上下邊緣點(diǎn)的傾斜方式為下傾斜時(shí),先進(jìn)行下傾斜鑲嵌。然后,為了避免產(chǎn)生漏洞現(xiàn)象,查看其上面相鄰樣點(diǎn)的傾斜方式是否為上傾斜,如果是,則對(duì)上面地震樣點(diǎn)和當(dāng)前樣點(diǎn)之間的矩形區(qū)域進(jìn)行鑲嵌,以對(duì)漏洞進(jìn)行縫補(bǔ);如果上面樣點(diǎn)的傾斜方式不是上傾斜,則不用進(jìn)行縫補(bǔ),如圖7 所示。

圖7 非上下邊緣點(diǎn)下傾斜鑲嵌示意圖

圖7 中,E為當(dāng)前樣點(diǎn),圖7(a)為當(dāng)B點(diǎn)為上傾斜時(shí)產(chǎn)生的漏洞現(xiàn)象,圖7(b)和(c)是可以采用的2 種漏洞縫補(bǔ)方法(此2 種方法任選一種即可)。

本文提出的鑲嵌方法,與文獻(xiàn)[15]每個(gè)四邊形單元中的4 個(gè)樣點(diǎn)的鑲嵌方式都由其左上角樣點(diǎn)的傾斜方式?jīng)Q定的方法不同,可最大程度地根據(jù)每一個(gè)地震樣點(diǎn)的傾斜方式對(duì)其進(jìn)行真正的自適應(yīng)鑲嵌。

2 實(shí)驗(yàn)結(jié)果及分析

為驗(yàn)證本文算法的有效性,在個(gè)人計(jì)算機(jī)上對(duì)新疆某工區(qū)的一條二維線中的標(biāo)準(zhǔn)segy 數(shù)據(jù)進(jìn)行了實(shí)驗(yàn)。運(yùn)行環(huán)境為Windows 7 旗艦版,處理器是i3-3240@3.40 GHz,顯卡為ATI Radeon HD 6450A/7450A(1 GB/聯(lián)想),內(nèi)存為4 GB(海力士DDR3L 1600 MHz),實(shí)驗(yàn)平臺(tái)是Microsoft Visual Studio 2010,圖形庫(kù)是Dirtect3D 10,著色器shader 4.0。

本文基于以上環(huán)境開(kāi)發(fā)的地震剖面三維可視化的系統(tǒng)流程圖,如圖8 所示。

圖8 本文地震剖面三維顯示系統(tǒng)流程圖

圖9～13 所示為不同鑲嵌方式對(duì)從標(biāo)準(zhǔn)segy剖面數(shù)據(jù)中選取的一塊含有下傾斜和上傾斜結(jié)構(gòu)的區(qū)域進(jìn)行可視化的結(jié)果,大小為200×80 樣點(diǎn)。

圖10 固定模式下傾斜鑲嵌

對(duì)比圖9～13 可知,固定模式鑲嵌沒(méi)有考慮相鄰道樣點(diǎn)之間的相關(guān)性,采用固定上傾斜或下傾斜鑲嵌方式鑲嵌頂點(diǎn)時(shí),分別會(huì)在下傾斜事件和上傾斜事件的頂部產(chǎn)生鋸齒(槽或脊),如圖9 和10 所示。文獻(xiàn)[15]采用低傾斜相關(guān)自適應(yīng)鑲嵌方法在一定程度上解決了此問(wèn)題,如圖11 所示。雖然相比固定鑲嵌模式,在消除鋸齒效果上有了一定的改進(jìn),但由于其會(huì)出現(xiàn)鑲嵌方式與樣點(diǎn)傾斜方式不一致的情況,存在鑲嵌錯(cuò)誤問(wèn)題,所以鋸齒現(xiàn)象仍然很明顯。如果根據(jù)其方法對(duì)每個(gè)樣點(diǎn)都根據(jù)其傾斜方式進(jìn)行自適應(yīng)鑲嵌,則會(huì)產(chǎn)生交疊和漏洞現(xiàn)象,如圖12所示。

圖11 文獻(xiàn)[15]低傾斜相關(guān)自適應(yīng)鑲嵌

圖12 對(duì)每個(gè)樣點(diǎn)都用文獻(xiàn)[15]方法進(jìn)行自適應(yīng)鑲嵌

圖13 是用本文提出的新的低傾斜相關(guān)自適應(yīng)鑲嵌方法得到的結(jié)果,該方法最大程度地保證了每個(gè)頂點(diǎn)的鑲嵌方式都與相應(yīng)樣點(diǎn)的傾斜方式一致,只有在產(chǎn)生三角形單元交疊或網(wǎng)格漏洞的情況下,少數(shù)頂點(diǎn)的鑲嵌方式才可能會(huì)與相應(yīng)樣點(diǎn)的傾斜方式不一致,而這里的不一致并不會(huì)造成傾斜事件頂部的鋸齒現(xiàn)象發(fā)生,因?yàn)槠洳⑽刺幱趦A斜事件上。

圖13 本文新的低傾斜相關(guān)自適應(yīng)鑲嵌

綜合以上對(duì)比分析可知,本文所提出的新的低傾斜相關(guān)自適應(yīng)鑲嵌方法效果較好,其有效性和可行性得到了驗(yàn)證。

3 結(jié)論

本文針對(duì)二維地震剖面數(shù)據(jù)的三維可視化問(wèn)題進(jìn)行研究,對(duì)于地震數(shù)據(jù)網(wǎng)格的頂點(diǎn)鑲嵌問(wèn)題,提出了一種新的低傾斜相關(guān)自適應(yīng)鑲嵌方法,并通過(guò)實(shí)驗(yàn)驗(yàn)證了算法的有效性和實(shí)用價(jià)值。進(jìn)一步研究的方向是對(duì)網(wǎng)格的細(xì)分方法進(jìn)行研究,以求將地震數(shù)據(jù)可視化的同時(shí),盡可能減少由鑲嵌帶來(lái)的鋸齒等現(xiàn)象。