地震剖面顯示的方法
【專利摘要】本發(fā)明提供了一種地震剖面顯示的方法，包括以下步驟：對(duì)于屏幕的每個(gè)像素，將其覆蓋范圍從像素空間映射到道集空間，獲得每個(gè)像素所覆蓋的地震道及采樣點(diǎn)；根據(jù)每個(gè)像素所覆蓋的地震道及采樣點(diǎn)對(duì)應(yīng)生成位于每個(gè)像素內(nèi)的波紋圖像；獲取每個(gè)像素內(nèi)的各個(gè)波紋圖像對(duì)本像素的亮度貢獻(xiàn)值，并將所述每個(gè)像素內(nèi)的各個(gè)波紋圖像對(duì)本像素的亮度貢獻(xiàn)值累加后作為本像素的亮度值；將各個(gè)像素的亮度值緩存于顯存；將所述顯存內(nèi)緩存的各個(gè)像素的亮度值繪制至屏幕顯示。本發(fā)明避免了由于抽道顯示導(dǎo)致的假頻現(xiàn)象。
【專利說明】地震剖面顯示的方法

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及地震剖面處理【技術(shù)領(lǐng)域】，尤其是涉及一種地震剖面顯示的方法。

【背景技術(shù)】
[0002] 地震數(shù)據(jù)尤其是疊加剖面數(shù)據(jù)的可視化是地震數(shù)據(jù)處理中非常重要的輔助手段， 從野外采集的地震數(shù)據(jù)經(jīng)過處理之后，需要借助于計(jì)算機(jī)技術(shù)以圖像的形式將其顯示出 來，從而為數(shù)據(jù)處理員、地質(zhì)分析人員提供一種可視化的判斷依據(jù)，提高石油勘探的精度和 效率。
[0003] 地震數(shù)據(jù)都是離散采樣，并按照道進(jìn)行排列，實(shí)際上就是空間中一系列的離散數(shù) 值，其可視化的主要形式是以波紋線、變密度的形式將數(shù)據(jù)繪制出來。當(dāng)以波紋線方式繪制 的時(shí)候，傳統(tǒng)的繪制方法是將離散的采樣點(diǎn)用直線連接起來，形成一條波形曲線，為了顯式 地表示波形振幅的正負(fù)，需要對(duì)地震波的某一個(gè)方向（一般是正向）進(jìn)行填充，這樣在地震 反射截面處就會(huì)有較清晰的反射同相軸。用戶縮小圖像的時(shí)候，由于現(xiàn)有計(jì)算機(jī)屏幕分辨 率的限制，當(dāng)縮小到一定程度之后，相鄰兩條波紋線之間的距離已經(jīng)小于一個(gè)像素，即使按 照一個(gè)像素進(jìn)行繪制，也會(huì)導(dǎo)致整個(gè)屏幕是全部涂黑的。一般的解決方法是對(duì)地震數(shù)據(jù)進(jìn) 行抽道處理，即每隔n條地震道抽取一條地震道作為代表進(jìn)行繪制。這種方法雖然實(shí)現(xiàn)簡 單，但是會(huì)產(chǎn)生假頻，這種現(xiàn)象會(huì)嚴(yán)重干擾解釋人員對(duì)地質(zhì)情況整體的把握。


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

[0004] 本發(fā)明的目的在于提供一種地震剖面顯示的方法，以解決地震剖面縮小顯示時(shí)出 現(xiàn)的假頻現(xiàn)象。
[0005] 為達(dá)到上述目的，本發(fā)明提供了一種地震剖面顯示的方法，包括以下步驟：
[0006] (SI)、對(duì)于屏幕的每個(gè)像素，將其覆蓋范圍從像素空間映射到道集空間，獲得每個(gè) 像素所覆蓋的地震道及采樣點(diǎn)；
[0007] (S2)、根據(jù)每個(gè)像素所覆蓋的地震道及采樣點(diǎn)對(duì)應(yīng)生成位于每個(gè)像素內(nèi)的波紋圖 像；
[0008] (S3)、獲取每個(gè)像素內(nèi)的各個(gè)波紋圖像對(duì)本像素的亮度貢獻(xiàn)值，并將所述每個(gè)像 素內(nèi)的各個(gè)波紋圖像對(duì)本像素的亮度貢獻(xiàn)值累加后作為本像素的亮度值；
[0009] (S4)、將各個(gè)像素的亮度值緩存于顯存；
[0010] (S5)、將所述顯存內(nèi)緩存的各個(gè)像素的亮度值繪制至屏幕顯示。
[0011] 本發(fā)明的地震剖面顯示的方法，步驟（SI)至（S4)由圖形處理器（Graphics Processing Unit，GPU)紋理生成線程完成；步驟（S5)由GPU紋理繪制線程完成。
[0012] 本發(fā)明的地震剖面顯示的方法，當(dāng)收到圖像縮放指令時(shí)，所述GPU紋理繪制線程 將所述屏幕當(dāng)前顯示的地震剖面進(jìn)行對(duì)應(yīng)的縮放處理生成臨時(shí)地震剖面并輸出顯示，并同 時(shí)觸發(fā)所述GPU紋理繪制線程，待所述GPU紋理繪制線程生成地震剖面時(shí)，將所述屏幕當(dāng)前 顯示的臨時(shí)地震剖面替換成所述GPU紋理繪制線程生成的地震剖面。
[0013] 本發(fā)明的地震剖面顯示的方法，所述屏幕由若干個(gè)固定大小的紋理單元所覆蓋， 所述每個(gè)紋理單元配置有一個(gè)獨(dú)立的紋理生成線程，且各個(gè)獨(dú)立的紋理生成線程并行生成 本紋理單元的地震剖面紋理。
[0014] 本發(fā)明的地震剖面顯示的方法，相鄰紋理單元的邊緣相互重疊。
[0015] 本發(fā)明的地震剖面顯示的方法，相鄰紋理單元的邊界相互重疊 2個(gè)像素。
[0016] 本發(fā)明的地震剖面顯示的方法，在所述將所述顯存內(nèi)緩存的各個(gè)像素的亮度值繪 制至屏幕顯示過程中，相鄰紋理單元在邊界重疊處減少一個(gè)像素。
[0017] 本發(fā)明的地震剖面顯示的方法，所述對(duì)于屏幕的每個(gè)像素，將其覆蓋范圍從像素 空間映射到道集空間，具體包括：
[0018] 將每個(gè)像素從屏幕坐標(biāo)轉(zhuǎn)換到世界坐標(biāo):W = P'M^S，其中F1為投影矩陣逆矩 陣，M4為模型矩陣逆矩陣，S為像素的屏幕坐標(biāo)，W為像素對(duì)應(yīng)的世界坐標(biāo)。
[0019] 本發(fā)明的地震剖面顯示的方法，所述獲取每個(gè)像素內(nèi)的各個(gè)波紋圖像對(duì)本像素的 亮度貢獻(xiàn)值，并將所述每個(gè)像素內(nèi)的各個(gè)波紋圖像對(duì)本像素的亮度貢獻(xiàn)值累加后作為本像 素的亮度值，具體包括：
[0020] 將每個(gè)像素的世界坐標(biāo)W均勻劃分成若干格子G，置每個(gè)格子Gi的初值為0 ;
[0021] 計(jì)算出與每個(gè)像素的世界坐標(biāo)W相交的波紋圖像，并將這些波紋圖像的覆蓋范圍 與本像素的每個(gè)格子Gi做相交測(cè)試，并將被覆蓋的格子Gi的值置1 ;
[0022] 計(jì)算每個(gè)像素中其值為1的格子的數(shù)目，并將該數(shù)目置為當(dāng)前像素的亮度值

【權(quán)利要求】
1. 一種地震剖面顯示的方法，其特征在于，包括以下步驟： (51) 、對(duì)于屏幕的每個(gè)像素，將其覆蓋范圍從像素空間映射到道集空間，獲得每個(gè)像素 所覆蓋的地震道及采樣點(diǎn)； (52) 、根據(jù)每個(gè)像素所覆蓋的地震道及采樣點(diǎn)對(duì)應(yīng)生成位于每個(gè)像素內(nèi)的波紋圖像； (53) 、獲取每個(gè)像素內(nèi)的各個(gè)波紋圖像對(duì)本像素的亮度貢獻(xiàn)值，并將所述每個(gè)像素內(nèi) 的各個(gè)波紋圖像對(duì)本像素的亮度貢獻(xiàn)值累加后作為本像素的亮度值； (54) 、將各個(gè)像素的亮度值緩存于顯存； (55) 、將所述顯存內(nèi)緩存的各個(gè)像素的亮度值繪制至屏幕顯示。
2. 根據(jù)權(quán)利要求1所述的地震剖面顯示的方法，其特征在于，步驟（S1)至（S4)由圖形 處理器GPU紋理生成線程完成；步驟（S5)由GPU紋理繪制線程完成。
3. 根據(jù)權(quán)利要求2所述的地震剖面顯示的方法，其特征在于，當(dāng)收到圖像縮放指令時(shí)， 所述GPU紋理繪制線程將所述屏幕當(dāng)前顯示的地震剖面進(jìn)行對(duì)應(yīng)的縮放處理生成臨時(shí)地 震剖面并輸出顯示，并同時(shí)觸發(fā)所述紋理繪制線程，待所述GPU紋理繪制線程生成地震 剖面時(shí)，將所述屏幕當(dāng)前顯示的臨時(shí)地震剖面替換成所述GPU紋理繪制線程生成的地震剖 面。
4. 根據(jù)權(quán)利要求2或3所述的地震剖面顯示的方法，其特征在于，所述屏幕由若干個(gè)固 定大小的紋理單元所覆蓋，所述每個(gè)紋理單元配置有一個(gè)獨(dú)立的紋理生成線程，且各個(gè)獨(dú) 立的紋理生成線程并行生成本紋理單元的地震剖面紋理。
5. 根據(jù)權(quán)利要求4所述的地震剖面顯示的方法，其特征在于，相鄰紋理單元的邊緣相 互重疊。
6. 根據(jù)權(quán)利要求5所述的地震剖面顯示的方法，其特征在于，相鄰紋理單元的邊界相 互重疊2個(gè)像素。
7. 根據(jù)權(quán)利要求6所述的地震剖面顯示的方法，其特征在于，在所述將所述顯存內(nèi)緩 存的各個(gè)像素的亮度值繪制至屏幕顯示過程中，相鄰紋理單元在邊界重疊處減少一個(gè)像 素。
8. 根據(jù)權(quán)利要求1所述的地震剖面顯示的方法，其特征在于，所述對(duì)于屏幕的每個(gè)像 素，將其覆蓋范圍從像素空間映射到道集空間，具體包括： 將每個(gè)像素從屏幕坐標(biāo)轉(zhuǎn)換到世界坐標(biāo)：w= 其中r1為投影矩陣逆矩陣，ivr1 為模型矩陣逆矩陣，s為像素的屏幕坐標(biāo)，w為像素對(duì)應(yīng)的世界坐標(biāo)。
9. 根據(jù)權(quán)利要求8所述的地震剖面顯示的方法，其特征在于，所述獲取每個(gè)像素內(nèi)的 各個(gè)波紋圖像對(duì)本像素的亮度貢獻(xiàn)值，并將所述每個(gè)像素內(nèi)的各個(gè)波紋圖像對(duì)本像素的亮 度貢獻(xiàn)值累加后作為本像素的亮度值，具體包括： 將每個(gè)像素的世界坐標(biāo)W均勻劃分成若干格子G，置每個(gè)格子Gi的初值為0 ; 計(jì)算出與每個(gè)像素的世界坐標(biāo)W相交的波紋圖像，并將這些波紋圖像的覆蓋范圍與本 像素的每個(gè)格子h做相交測(cè)試，并將被覆蓋的格子匕的值置1 ; 計(jì)算每個(gè)像素中其值為1的格子的數(shù)目，并將該數(shù)目置為當(dāng)前像素的亮度值
10. 根據(jù)權(quán)利要求9所述的地震剖面顯示的方法，其特征在于，將每個(gè)像素的世界坐標(biāo) W均勻劃分成16X16個(gè)格子G。
11. 根據(jù)權(quán)利要求4所述的地震剖面顯示的方法，其特征在于，每個(gè)紋理單元覆蓋 512X512 或者 1024X1024 個(gè)像素。
【文檔編號(hào)】G01V1/34GK104407384SQ201410743516
【公開日】2015年3月11日 申請(qǐng)日期:2014年12月8日 優(yōu)先權(quán)日:2014年12月8日
【發(fā)明者】梁兼棟, 呂偉偉 申請(qǐng)人:中國石油天然氣集團(tuán)公司, 中國石油集團(tuán)東方地球物理勘探有限責(zé)任公司