專利名稱：對地震斷面或地震塊進行年代地層學分析的方法
技術領域：
本發(fā)明涉及一種對地震斷面或者地震塊，也就是說對地震斷面或者地震塊所包含的地質(zhì)記錄數(shù)據(jù)，進行年代地層學分析的方法。
本發(fā)明從石油勘探的角度進行說明，但可以從地球物理學領域轉(zhuǎn)用到地質(zhì)學領域。
現(xiàn)有技術本發(fā)明的方法應用于地震斷面或者地震塊。地震斷面是將稱為地震記錄跡的一維采樣信號并列在一個平面內(nèi)形成的。類似地，地震塊是將地震記錄跡并列在一個空間中形成的。術語“地震剖面”既指上述地震斷面，也指地震塊的某個切面。一個地震剖面提供了在該剖面平面內(nèi)所包含的并列地震記錄跡的一個視圖。這些視圖就是地震圖像，為了說明如何實施本發(fā)明的方法，下面將所述地震圖像稱為地震圖像剖面。在一個地震圖像中，像素的亮度與所述一維信號所代表的地震波幅成正比。
對地震斷面或者地震塊的年代地層學分析涉及地震斷面或者地震塊中地震地層的合成。在地層的合成方面，已有好幾種方法。這些方法的結(jié)果的好壞取決于用地震剖面為其提供地震圖像的地質(zhì)環(huán)境。因此，在地層以單斜為主的地區(qū)，通過測量相鄰地震記錄跡之間的相似性來進行地層的合成，可以得到較好的結(jié)果。另一方面，在地質(zhì)情況受到較多擾動的地區(qū)，最好首先計算出相鄰像素之間亮度的梯度矢量，然后通過對所計算出的梯度矢量的方向域求積分而完成地層的合成。
由Marc Donias于1999年1月28日提交給波爾多第一大學的題為《通過主分量分析和曲率估算表征方向域—在地震圖像中的應用》的論文詳細描述了上述地層合成的原理。
對地震圖像剖面的每個像素進行地層合成，會生成與圖像中的像素一樣多的地層。在局部地質(zhì)圖中體現(xiàn)為標志層的地震地層不能相交。另一方面，地震地層可以收斂、局部合并、匯合，或者甚至分叉。地層的合并導致了合成累積概念的產(chǎn)生。
為了進行地層合成的累積，要定義一個在大小上與地震圖像相同的矩陣。該矩陣的每個元素與地震圖像的一個像素關聯(lián)，并被賦以初始的零值。對于地震圖像的每一個像素，計算出一個連續(xù)曲線，對應于通過該像素的合成地層。所有的連續(xù)曲線都在圖像垂直尺度的橫向上。在計算連續(xù)曲線時，圖像中的像素每被一個連續(xù)曲線通過一次，與該像素關聯(lián)的矩陣元素就增加一個單位。
在給定的像素處，對圖像剖面垂直尺度的橫向上的連續(xù)曲線的計算是這樣進行的計算所述選定像素的一個鄰域中所有像素的亮度的梯度，然后從在所述鄰域上獲得的梯度值計算出一個局部梯度，再將該局部梯度賦給所述選定像素。然后，從所述選定像素開始，在所述局部梯度及其加性逆元所指示的兩個方向上，在橫向上一個像素一個像素地反復重復前述兩個步驟，直到圖像的側(cè)垂直邊界，從而形成連續(xù)曲線。
當?shù)卣饒D像剖面的所有像素都被掃描之后，進行上述合成累積的矩陣就表現(xiàn)為一個新的圖像，其中，每個像素具有與相應矩陣元素中的累計數(shù)成比例的亮度，該累計數(shù)至少為一。在該圖像上所看到的邊界可以很好地體現(xiàn)亞土層的地質(zhì)層位構造。
本發(fā)明的效果本發(fā)明方法利用通過合成累積獲得的圖像來測定地質(zhì)沉積物沉積時的狀態(tài)，而不是現(xiàn)在以地質(zhì)地層形式所觀察到的狀態(tài)。為此，本方法定義了一種從按照地震時序測量的地震剖面垂直尺度到按照地質(zhì)時期測量的地質(zhì)垂直尺度的轉(zhuǎn)換。該方法因此使得可以確定影響地質(zhì)地層沉積的沉積速率。特別地，該方法可以突出顯示地質(zhì)上的間斷，即侵蝕和間隙。
本發(fā)明概述本發(fā)明的主題是一種對地震圖像剖面S進行年代地層學分析的方法，所述地震圖像剖面S包括一個水平尺度或者說寬度，以及一個在亞土層方向的垂直尺度或者說高度，并由像素列構成，該方法包括—定義一個在大小上與圖像剖面S相同的矩陣M，該矩陣的每個元素與圖像剖面S的一個像素關聯(lián)，并被賦以初始的零值，—對于圖像剖面S的每一個像素i，計算出一個通過該像素、在圖像剖面S的垂直尺度的橫向上的連續(xù)曲線Ci，—圖像剖面S中的像素每被一個連續(xù)曲線Ci通過一次，與該像素關聯(lián)的矩陣M的元素就增加一個，所述方法的特征在于它還包括—對于矩陣M的每一列c，繪制一個由許多組構成的直方圖Hc，所述組的數(shù)目等于所述列c的元素的數(shù)目，每個組對應于所述列c的一個元素，并包括許多樣本，樣本數(shù)等于矩陣M的有關元素中的累計值，該值等于通過與所述元素關聯(lián)的像素的曲線數(shù)，為每一列繪制的所述直方圖中分布的樣本總數(shù)等于圖像剖面S中的像素總數(shù)，—對每一個直方圖Hc進行均衡處理，得到均衡直方圖Hc’，—定義一個空的圖像剖面S’，其按像素計算的寬度等于圖像剖面S按像素計算的寬度，其按像素計算的高度等于所述直方圖Hc’的組數(shù)，—通過將所述均衡直方圖Hc’的關聯(lián)組的內(nèi)容的基數(shù)分配給所述列c’的每一個像素，將Hc’所確定的分布賦給圖像剖面S’的每一列c’，—在圖像剖面S’中，確定含有樣本的鄰接像素的組，并標記所述每一個組，—將在圖像剖面S’中分配給所述像素組的標記分配給圖像剖面S的每一個像素，并顯示所述被標記的圖像剖面S。
根據(jù)本發(fā)明另一項特征，在給定的像素i處，對圖像剖面S垂直尺度的橫向上的連續(xù)曲線Ci的計算是這樣進行的—計算像素i的一個鄰域Vi中所有像素的亮度的梯度，—從在所述鄰域Vi上獲得的梯度值計算出一個局部梯度Gi，并將該梯度值Gi賦給像素i，—從像素i開始，在所述梯度Gi及其加性逆元-Gi所指示的兩個方向上，在橫向上一個像素一個像素地反復重復前述兩個步驟，直到圖像剖面S的側(cè)垂直邊界。
本發(fā)明也能應用于三維地震圖像塊。因此本發(fā)明的主題也是一種對包含兩個水平尺度即寬度和深度以及一個在亞土層方向的垂直尺度或者說高度、由像素列構成的地震圖像塊B進行年代地層學分析的方法，該方法包括—定義一個在大小上與圖像塊B相同的塊N，該塊的每個元素與圖像塊B的一個像素關聯(lián)，并被賦以初始的零值，—對于圖像塊B的每一個像素i，計算出一個通過該像素、在圖像塊B的垂直尺度的橫向上的連續(xù)表面Si，—圖像塊B中的像素每被一個連續(xù)表面Si穿過一次，與該像素關聯(lián)的塊N的元素就增加一個，所述方法的特征在于它還包括—對于塊N的每一列c，繪制一個由許多組構成的直方圖Hc，所述組的數(shù)目等于所述列c的元素的數(shù)目，每個組對應于所述列c的一個元素，并包括許多樣本，樣本數(shù)等于塊N的有關元素中的累計值，該值等于通過與所述元素關聯(lián)的像素的表面數(shù)，為每一列繪制的所述直方圖中分布的樣本總數(shù)等于圖像塊B中的像素總數(shù)，—對每一個直方圖Hc進行均衡處理，得到均衡直方圖Hc’，—定義一個空的圖像塊B’，其按像素計算的寬度和深度等于圖像塊B按像素計算的寬度和深度，其按像素計算的高度等于所述直方圖Hc’的組數(shù)，—通過將所述均衡直方圖Hc’的關聯(lián)組的內(nèi)容分配給所述列c’的每一個像素，將Hc’所確定的分布賦給圖像塊B’的每一列c’，—在圖像塊B’中，確定含有樣本的鄰接像素的組，并標記所述每一個組，—將在圖像塊B’中分配給所述像素組的標記分配給圖像塊B的每一個像素，并顯示所述被標記的圖像塊B。
根據(jù)本發(fā)明另一項特征，在給定的像素i處，對圖像塊B垂直尺度的橫向上的連續(xù)表面Si的計算是這樣進行的—計算像素列Ki中每一個像素i的一個鄰域中所有像素的亮度的梯度，—對像素列Ki的每一個像素i，對在像素i的鄰域Vi中計算出的梯度進行主分量分析，以確定在像素列Ki的像素i處與所述連續(xù)表面Si相切的平面所指示的兩個方向矢量，—從像素列Ki直到圖像塊B的側(cè)垂直邊界，同心地從一列到另一列，反復重復前述兩個步驟。


圖1a和圖1b按照不同的比例示出了應用本發(fā)明方法之前的同一個地震圖像剖面。
圖2a示出了一條與所述地震圖像剖面的一個像素相關的連續(xù)曲線，圖2b則是所有連續(xù)曲線的視圖。
圖3a示出一個與圖2b中一條垂線相關的直方圖，圖3b示出經(jīng)過均衡處理后的同一個直方圖。
圖4是一個均衡圖像剖面，稱為年代地層學圖像剖面。
圖5a是一個二元化的均衡圖像剖面，圖5b是經(jīng)過標記后的同一個圖像剖面。
圖6是應用了本發(fā)明方法后經(jīng)過標記的地震圖像剖面。
詳細說明本發(fā)明的方法是一種對地震圖像剖面進行自動年代地層學分析的方法。下面參照附圖給出一種該方法的實施例。
圖1a是一個地震圖像剖面S。該圖像剖面S表示的是從地震勘測得到的亞土層圖像。該圖像剖面S包括兩個尺度。它由沿一水平軸的水平尺度和沿一在亞土層方向的垂直軸的垂直尺度限定而成。圖像剖面S由像素組成，像素按照水平軸上的某個水平間距及垂直軸上的某個垂直間距均勻分布。圖像剖面S尤其包含有一系列像素列，列數(shù)等于圖像剖面的水平尺度除以所述水平間距的商，每像素列的像素數(shù)等于所述垂直尺度除以所述垂直間距的商。特別地，圖1a中的垂直粗實線10表示一個像素列，該像素列是下文對本發(fā)明的方法進行說明的各附圖的基礎。
為了實施本發(fā)明的方法，定義一個在大小上與圖像剖面S相同的矩陣M。該矩陣M的行數(shù)等于圖像剖面的一列所含的像素數(shù)，該矩陣M的列數(shù)等于圖像剖面S的列數(shù)，并等于所述圖像剖面中一行所包含的像素數(shù)。這樣，該矩陣M所包含的元素數(shù)就與所述圖像剖面S中的像素數(shù)一樣多，并且每個元素都與圖像剖面S的一個像素相關聯(lián)。矩陣M的所有元素均為整型，并被賦以初始的零值。
對于圖像剖面S的每一個像素i，計算出一個通過該像素、在圖像剖面S的垂直尺度的橫向上的連續(xù)曲線Ci。該曲線Ci的計算涉及像素i處的亮度局部梯度Gi的計算。
與某個像素相關聯(lián)的亮度，定義為某個地震參數(shù)、例如地震波信號的振幅按某種調(diào)色板的表達，例如按灰度值調(diào)色板的表達。在本例所采用的條件下，高振幅用亮像素表示，低振幅用暗像素表示。然后就可以計算出某個像素i及其相鄰像素之間的亮度梯度Gi了。該梯度Gi事實上是有關地震參數(shù)的梯度。將該梯度Gi賦給像素i。該梯度Gi包括一個水平分量和一個垂直分量，尤其是通過對在像素i鄰域Vi的所有像素上計算而得的所有梯度進行主分量分析而確定。所述鄰域Vi由一個以像素i為中心的窗口確定，這個窗口的大小最好為7像素乘以7像素。所述鄰域Vi中的梯度測量值的由第一慣性軸的方向給定的梯度測量束主方向使得確定所述梯度Gi的局部方向成為可能，該方向被分配給像素i。
圖2a示出了一條與某個像素i相關聯(lián)的連續(xù)曲線Ci。該曲線Ci是一條平滑化的連續(xù)曲線，因為該曲線是在所述圖像中從像素i向特性與像素i類似的像素延伸。該曲線Ci在示于與圖2a比例尺相同的圖1b中的所述圖像剖面S的垂直尺度的橫向上延伸，在該圖像剖面上，可以隱約辨認出由曲線Ci合成的地層。該曲線Ci是在相繼的連續(xù)計算的梯度Gi及其加性逆元-Gi所確定的相繼的方向上一個像素一個像素地處理而獲得的。結(jié)果，所述連續(xù)曲線Ci在構成該曲線的每一個像素點處的導數(shù)值都是在同一像素處計算而得的局部梯度值。
每有連續(xù)曲線Ci通過像素i，與像素i相應的矩陣M的元素就增加一個單位。矩陣M的元素就好比是與圖像剖面S的像素相關聯(lián)的計數(shù)器。
圖2b是所有橫向連續(xù)曲線的視圖。橫向曲線的線條越顯得粗實，疊加的橫向曲線數(shù)就越多。更具體地說，如果圖像剖面S表示的例如是地震波信號的振幅值，從像素i起始的曲線Ci將基于在像素i處觀察到的振幅相似性在圖像剖面S上尋找一個振幅連續(xù)值。如果像素i處的振幅值非常高，從該圖像中的像素i出發(fā)的曲線就可能是沿著一條地質(zhì)標志層。
在計算出所有連續(xù)曲線之后，矩陣M的元素就含有連續(xù)曲線通過圖像剖面S的每一個像素的嚴格按照元素-像素對應關系的計數(shù)結(jié)果。這個結(jié)果就是示于圖2b的矩陣M的精確內(nèi)容。
對于對應于從圖2b所示圖像中選取的某像素列的矩陣M的每一列，繪制一個直方圖Hc，就如圖3a所示。圖3a所示的直方圖與圖2b中的垂直黑實線所指示的列c相關聯(lián)。該直方圖所包含的組數(shù)等于所述列c中的像素數(shù)。直方圖的每個組與矩陣M的所述列中的一個元素相關聯(lián)，該元素對應于圖像剖面S中相應的像素列中的一個像素；通過該像素的連續(xù)曲線的數(shù)目越多，Hc的相應的組就越充實。因此，圖3的直方圖Hc展示了連續(xù)曲線在標于圖像剖面S上的垂直線10上的分布，這個分布精確存儲在與圖像剖面S關聯(lián)的矩陣M的相應列中。
與不同的列相關聯(lián)的不同的直方圖的組所包含的樣本數(shù)在均為正整數(shù)的一個最小值和一個最大值之間變化。最小值是一，因為總是至少有一條連續(xù)曲線通過每個像素。最大值是可變的，但受到連續(xù)曲線總數(shù)的限制，該連續(xù)曲線總數(shù)不外乎是圖像剖面S的總像素數(shù)。
然后用一種均衡算法對每一個直方圖Hc進行處理，生成相應的均衡直方圖Hc’。為此，例如可以利用由Addison-Wesley于1992年出版的R.Gonzalez和R.Woods的題為《數(shù)字圖像處理》的著作中所描述的算法，或者可以利用由Masson于1995年出版的J-P Coquerez和S.Philipp的題為《圖像分析濾波和分段》的著作中所描述的算法。對圖3a的直方圖實施的直方圖均衡算法使總體重新分布，使之按照圖3b的直方圖Hc’所示的新的分布進行排列。主要目的是提高圖像對比度的直方圖均衡算法在這里的目的是將所述排列總體以分布更加均衡的方式分為更多的組，從而重新分布所述排列總體。均衡直方圖Hc’的組數(shù)等于直方圖Hc的組數(shù)乘以依所需分辨率而定的擴展系數(shù)。
在實施本發(fā)明的方法時，定義一個新的圖像剖面S’，其按像素計的寬度等于圖像剖面按像素計的寬度，其按像素計的高度等于所述直方圖Hc’的組數(shù)。開始時圖像剖面S’是空的，也就是說其像素看起來是黑的。
將相應的均衡直方圖所確定的分布賦給圖像剖面S’的每一列，從而基于所有的均衡直方圖在所述圖像剖面S’中構造一個新的圖像。為此，用新分布到所述均衡直方圖的組中的樣本來計算與每個組相關聯(lián)的每個像素的亮度。圖像剖面S’中某個像素的亮度，例如在灰度值調(diào)色板中的亮度，取決于與該像素相關聯(lián)的組中所包含的樣本數(shù)。這樣，如果在與該像素相關聯(lián)的組中沒有樣本，該像素的亮度就為零，該像素在圖像剖面S’中就顯得是黑的。類似的，如果所述關聯(lián)組包含較多或較少的樣本數(shù)，亮度就是非零值，像素就是或亮或暗的灰色。
重要的是要注意到，對于所有的直方圖，所述排列總體都是一樣的，無論其是否被均衡處理過。每個直方圖表現(xiàn)出不同的分布，但對于所有的直方圖來說，所有組的內(nèi)容的總和是一樣的。就是這種不變性使得可以將圖像剖面S’視為年代地層學圖像剖面，也就是說，將所述圖像剖面所凸現(xiàn)出來的邊界視為地質(zhì)學意義上的等時線。
圖4示出了從均衡直方圖構造的年代地層學圖像剖面S’。在這里所作的變換，使得可以從按米或毫秒標度的深度方向的地球物理學垂直比例尺，轉(zhuǎn)換到按百萬年標度的地質(zhì)時代比例尺。特別地，可以注意到在圖4中劃定最初在圖1a中選定的像素列的垂直線11。
然后，將含有非零值的每一個像素配置為“一”值，而在無信息的像素處保持零值。在這樣二元化了的如圖5a所示的圖像中，在黑色的背景上呈現(xiàn)白色的區(qū)域。白色區(qū)域表示該區(qū)域的像素包含至少某些地震信息，而黑色背景則表示沒有信息，也就是說，在均衡直方圖的關聯(lián)組中沒有樣本。從地質(zhì)學的角度來看，黑色背景揭示了間斷期或者侵蝕期，因為，從年代地層學圖像剖面S’的垂直方向看，黑色段所指示的時期沒有沉積物的沉積。
標記所述二元化圖像，為相連的白色區(qū)域指定完全不同的圖案或者顏色，使獨立的白色區(qū)域之間決不具有相同的圖案或者顏色。圖5b示出了標記所述二元化圖像剖面5a的結(jié)果。
最后，對圖像剖面5b的每一列進行與直方圖均衡化相逆的變換，也就是說，被均衡化直方圖排列的所述樣本被恢復到原始直方圖的原始組。然后，為每一個像素分配上述標記過程所分配的顏色，從而重新形成原始剖面。所述逆變換的結(jié)果示于圖6。圖5b中從12到15的標記區(qū)對應于圖6中類似標記的區(qū)域。具體來說，在圖5b中可以注意到，區(qū)域12幾乎就被分割為兩個區(qū)域12和12’，但由于有一個小狹部w將這兩個區(qū)域連接起來，這兩個區(qū)域就使用同一個標記。在圖6中，也可以從12通到12’，因為，由于狹部w的存在，這兩個區(qū)域之間的分界是不連續(xù)的。
圖6非常接近于真正的地質(zhì)剖面。具有相同圖案的區(qū)域是沉積物的沉積比較規(guī)則和均一的區(qū)域。具有不同圖案又相互接觸的兩個區(qū)域往往表示有地質(zhì)構造的變化，也就是說存在新的沉積體系，例如從海進現(xiàn)象到海退現(xiàn)象的轉(zhuǎn)換所導致的新的沉積體系。而且，通過識別包含在示于圖6的經(jīng)過標記的圖像剖面中的地質(zhì)標志層，可以在空間中協(xié)調(diào)一致地確定所觀察到的所有地質(zhì)事件的年代。
本發(fā)明的整個方法也可以在三個尺度上應用于地震圖像塊。一個地震圖像塊由兩個水平軸和一個在亞土層方向的垂直軸確定。在這種實施方案中，在圖像塊垂直尺度的橫向上的連續(xù)表面相當于在所述圖像剖面的垂直尺度的橫向上的連續(xù)曲線，在將所述方法應用于圖像塊時，只有計算這些表面的步驟有所不同。
為了確定所述圖像塊的通過該圖像塊所有像素的連續(xù)表面，該圖像塊中所有的像素列都要被依次掃描。對于圖像塊中選定的某像素列Ki，在組成該像素列Ki的每個像素處與所述連續(xù)表面相切的所有平面同時被確定。為此，首先在像素列Ki的每一個像素i處，在該有關像素的一個鄰域Vi中，計算出一個亮度局部梯度。這里，所述鄰域Vi是一個以有關像素為中心的小立方體，例如是一個7×7×7像素的立方體。所述梯度是具有在所述圖像塊三個尺度方向的三個分量的矢量。在像素列Ki的每一個像素i處，所述梯度是這樣確定的對在像素i鄰域Vi中所有像素處計算出的所有梯度進行主分量分析。對在像素i鄰域Vi中計算出的梯度進行的三維主分量分析可以確定三個正交分量，即一個與第一慣性軸關聯(lián)的主矢量Fi1，一個與第二慣性軸相關聯(lián)的主矢量Fi2和一個與第三慣性軸相關聯(lián)的主矢量Fi3，所述第一慣性軸確定像素i的切平面的最大斜度線，所述第二慣性軸與第一慣性軸正交，并沿像素i的切平面方向，所述第三慣性軸正交于像素i的切平面。這樣，在像素列Ki的每個像素i處，就獲得了沿著與所要的連續(xù)表面相切并通過該像素的平面方向的一對方向矢量(Fi1，F(xiàn)i2)。反復重復上述過程，從像素列Ki直到圖像塊的側(cè)垂直邊界，同心地從一列到另一列進行所述主分量分析。
權利要求
1.一種對地震圖像剖面S進行年代地層學分析的方法，所述地震圖像剖面S包括一個水平尺度或者寬度，以及一個在亞土層方向的垂直尺度或者高度，并由像素列構成，該方法在于—定義一個在大小上與圖像剖面S相同的矩陣M，該矩陣的每個元素與圖像剖面S的一個像素關聯(lián)，并被賦以初始的零值，—對于圖像剖面S的每一個像素i，計算出一個通過該像素、在圖像剖面S的垂直尺度的橫向上的連續(xù)曲線Ci，—圖像剖面S中的像素每被一個連續(xù)曲線Ci通過一次，與該像素關聯(lián)的矩陣M的元素就增加一個單位，該方法的特征在于它還包括—對于矩陣M的每一列c，繪制一個由許多組構成的直方圖Hc，所述組的數(shù)目等于所述列c的元素的數(shù)目，每個組對應于所述列c的一個元素，并包括許多樣本，樣本數(shù)等于矩陣M的有關元素中的累計值，該值等于通過與所述元素關聯(lián)的像素的曲線數(shù)，為每一列繪制的所述直方圖中分布的樣本總數(shù)等于圖像剖面S中的像素總數(shù)，—對每一個直方圖Hc進行均衡處理，得到均衡直方圖Hc’，—定義一個空的圖像剖面S’，其按像素計算的寬度等于圖像剖面S按像素計算的寬度，其按像素計算的高度等于所述直方圖Hc’的組數(shù)，—通過將所述均衡直方圖Hc’的關聯(lián)組的內(nèi)容的基數(shù)分配給所述列c’的每一個像素，將Hc’所確定的分布賦給圖像剖面S’的每一列c’，—在圖像剖面S’中，確定含有樣本的鄰接像素的組，并標記所述每一個組，—將在圖像剖面S’中分配給所述像素組的標記分配給圖像剖面S的每一個像素，并顯示所述被標記的圖像剖面S。
2.根據(jù)權利要求1所述的方法，其特征在于，在給定的像素i處，對圖像剖面S垂直尺度的橫向上的連續(xù)曲線Ci的計算是這樣進行的—計算像素i的一個鄰域Vi中所有像素的亮度的梯度，—從在所述鄰域Vi上獲得的梯度值計算出一個局部梯度Gi，并將該梯度值Gi賦給像素i，—從像素i開始，在所述梯度Gi及其加性逆元-Gi所指示的兩個方向上，在橫向上一個像素一個像素地反復重復前述兩個步驟，直到圖像剖面S的側(cè)垂直邊界。
3.一種對地震圖像塊B進行年代地層學分析的方法，所述地震圖像塊包含兩個水平尺度、即寬度和深度，以及一個在亞土層方向的垂直尺度或者高度，并由像素列構成，該方法在于—定義一個在大小上與圖像塊B相同的塊N，該塊的每個元素與圖像塊B的一個像素關聯(lián)，并被賦以初始的零值，—對于圖像塊B的每一個像素i，計算出一個通過該像素、在圖像塊B的垂直尺度的橫向上的連續(xù)表面Si，—圖像塊B中的像素每被一個連續(xù)表面Si通過一次，與該像素關聯(lián)的塊N的元素就增加一個單位，該方法的特征在于它還包括—對于塊N的每一列c，繪制一個由許多組構成的直方圖Hc，所述組的數(shù)目等于所述列c的元素的數(shù)目，每個組對應于所述列c的一個元素，并包括許多樣本，樣本數(shù)等于塊N的有關元素中的累計值，該值等于通過與所述元素關聯(lián)的像素的表面數(shù)，為每一列繪制的所述直方圖中分布的樣本總數(shù)等于圖像塊B中的像素總數(shù)，—對每一個直方圖Hc進行均衡處理，得到均衡直方圖Hc’，—定義一個空的圖像塊B’，其按像素計算的寬度和深度等于圖像塊B按像素計算的寬度和深度，其按像素計算的高度等于所述直方圖Hc’的組數(shù)，—通過將所述均衡直方圖Hc’的關聯(lián)組的內(nèi)容分配給所述列c’的每一個像素，將Hc’所確定的分布賦給圖像塊B’的每一列c’，—在圖像塊B’中，確定含有樣本的鄰接像素的組，并標記所述每一個組，—將在圖像塊B’中分配給所述像素組的標記分配給圖像塊B的每一個像素，并顯示所述被標記的圖像塊B。
4.根據(jù)權利要求3所述的方法，其特征在于，在給定的像素i處，對圖像塊B垂直尺度的橫向上的連續(xù)表面Si的計算是這樣進行的—計算像素列Ki中每一個像素i的一個鄰域中所有像素的亮度的梯度，—對像素列Ki的每一個像素i，對在像素i的鄰域Vi中計算出的梯度進行主分量分析，以確定在像素列Ki的像素i處與所述連續(xù)表面Si相切的平面所指示的兩個方向矢量，—從像素列Ki直到圖像塊B的側(cè)垂直邊界，同心地從一列到另一列，反復重復前述兩個步驟。
全文摘要
本發(fā)明方法利用累積地震地層的連續(xù)代表曲線而獲得的圖像,來測定構成地震剖面的地層在沉積物沉積時的狀態(tài),而不是現(xiàn)在所觀察到的狀態(tài)。為此,本方法定義了一種從按照地震時序測量的地震剖面垂直尺度到按照地質(zhì)時期測量的地質(zhì)垂直尺度的轉(zhuǎn)換。利用均衡直方圖,本方法可以得到均衡地震剖面,再用來確定影響地質(zhì)地層沉積的沉積速率。特別地,該方法可以突出顯示地質(zhì)上的間斷,即侵蝕和間隙。
文檔編號G01V1/30GK1332378SQ0111712
公開日2002年1月23日 申請日期2001年4月26日 優(yōu)先權日2000年4月26日
發(fā)明者納曼·凱司克斯, 塞巴斯蒂安·吉龍, 馬爾克·多尼亞斯, 皮埃爾·貝龍, 法比安·波蓋特 申請人:埃勒夫勘探產(chǎn)品公司