傾角傳播法層位自動(dòng)追蹤方法
【專利摘要】本發(fā)明涉及一種地震資料構(gòu)造解釋技術(shù)，是一種傾角傳播法層位自動(dòng)追蹤方法，根據(jù)地震資料的特點(diǎn)和地震層位追蹤要求對(duì)經(jīng)典的“洪水填充”圖像處理技術(shù)進(jìn)行改造，構(gòu)建“逆向洪水填充”算法，采用反向種子點(diǎn)搜索、設(shè)置地層阻尼系數(shù)和斷層阻擋等技術(shù)實(shí)現(xiàn)層位的自動(dòng)追蹤。本發(fā)明可增強(qiáng)地震層位追蹤的連續(xù)性，提高斷層附近層位追蹤的準(zhǔn)確性，適應(yīng)高精度層位自動(dòng)追蹤和地震資料高效構(gòu)造解釋的要求，既能適應(yīng)二維工區(qū)的層位追蹤，又能適應(yīng)三維工區(qū)的層位追蹤。
【專利說明】?jī)A角傳播法層位自動(dòng)追蹤方法

【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及石油天然氣勘探與開發(fā)技術(shù)，屬于地震資料構(gòu)造解釋中一種在地層傾角控制下的層位自動(dòng)追蹤方法。

【背景技術(shù)】
[0002]地震資料的構(gòu)造解釋與分析是地震勘探工程的重要階段，在這一階段，既可利用地震資料搞清地下構(gòu)造特征，又為尋找有利的油氣儲(chǔ)層和油氣藏提供基礎(chǔ)資料。基于地震資料的構(gòu)造特征分析不但是儲(chǔ)層特征分析的基礎(chǔ)而且直接影響儲(chǔ)層特征分析的質(zhì)量。
[0003]地震層位是指與地下介質(zhì)中某一地層所對(duì)應(yīng)的地震反射界面。一般情況下，在同一斷塊中，一個(gè)地震反射界面表現(xiàn)為一個(gè)可連續(xù)追蹤的地震反射同相軸。
[0004]地震資料的構(gòu)造解釋特別是層位解釋是一項(xiàng)非常費(fèi)時(shí)、費(fèi)力的工作，在過去的幾十年中，此項(xiàng)工作常常占到整個(gè)地震資料解釋與分析工作的50%-80%。因此，基于計(jì)算機(jī)智能技術(shù)的層位自動(dòng)識(shí)別與解釋(簡(jiǎn)稱為層位自動(dòng)追蹤或?qū)游蛔粉?非常重要。
[0005]地震層位自動(dòng)追蹤是現(xiàn)代地震資料解釋中最關(guān)鍵的技術(shù)之一，它直接體現(xiàn)了一個(gè)解釋系統(tǒng)在構(gòu)造解釋中的能力和效率。一個(gè)好的層位自動(dòng)追蹤方法可以大幅度提高構(gòu)造解釋的效率，提高解釋精度，將地質(zhì)綜合研究人員的精力從繁瑣的層位拾取中解放出來，并將研究重點(diǎn)集中到儲(chǔ)層預(yù)測(cè)中。
[0006]在地震資料解釋中，由于地震剖面上的地震同相軸是某一地層的地震響應(yīng)或某些地層的綜合地震響應(yīng)，地震同相軸與地層具有一定的對(duì)應(yīng)關(guān)系，因此在實(shí)際地震剖面上對(duì)同相軸進(jìn)行連續(xù)性追蹤是層位解釋的先決條件，地震同相軸上的傾角和方位角變化特征是地震層位追蹤的基礎(chǔ)。
[0007]多年來，層位追蹤方法一直在研究和發(fā)展過程中，并出現(xiàn)了眾多的層位追蹤技術(shù)。其中以相關(guān)法應(yīng)用最為廣泛，在很多的解釋系統(tǒng)中都采用了這種技術(shù)，但相關(guān)法存在精度低及易陷入局部極值等問題。為了解決互相關(guān)法所存在的問題，多種新的方法被提出并進(jìn)行了嘗試，如波形特征法(Waveform Controlled,)、最大能量法(Maximum Energy)、神經(jīng)網(wǎng)絡(luò)法(Neural Network)、線性插值法(Linear Interpolat1n)、高階量統(tǒng)計(jì)法(High-randStatistics)、邊緣探測(cè)法(Edge Detect1n)及圖像對(duì)比法等。在這些方法中，或多或少存在著問題，如精度問題、效率問題及地震資料的適應(yīng)性等問題。
[0008]地震同相軸是地層界面的地震響應(yīng)，傾角和方位角是地震同相軸的本質(zhì)特征。在傳統(tǒng)層位自動(dòng)追蹤技術(shù)中，盡管或多或少地利用了地震同相軸的傾角信息作為層位自動(dòng)追蹤的依據(jù)，但并未明確提出利用傾角信息控制種子點(diǎn)的傳播，并未強(qiáng)調(diào)傾角信息特別是不同方位的傾角信息在層位自動(dòng)追蹤過程中的重要性，導(dǎo)致每一種層位自動(dòng)追蹤方法的實(shí)現(xiàn)過程不一致，層位自動(dòng)追蹤精度差異很大，且嚴(yán)重受制于傾角估算的精度。因此，明確提出以地震同相軸傾角信息作為層位自動(dòng)追蹤的主要控制因素，可將不同類型的層位自動(dòng)追蹤技術(shù)統(tǒng)一起來，并根據(jù)傾角估算技術(shù)的特點(diǎn)進(jìn)一步提高層位自動(dòng)追蹤的精度。
[0009]從本質(zhì)上說，在地層的沉積過程中，地層及地層界面的形成是受流體運(yùn)動(dòng)所控制的，這種地質(zhì)運(yùn)動(dòng)過程可在一定程度上用流體力學(xué)來描述。在傳統(tǒng)層位自動(dòng)追蹤技術(shù)中，為節(jié)約計(jì)算機(jī)資源，常常沿固定方向?qū)崿F(xiàn)種子點(diǎn)的擴(kuò)散。盡管這種技術(shù)效率較高，但這并不符合流體力學(xué)的規(guī)律，會(huì)導(dǎo)致層位追蹤結(jié)果與實(shí)際層位存在較大的誤差。此外，一些層位自動(dòng)追蹤算法特別是三維層位自動(dòng)追蹤算法雖然采用了圖像處理中的洪水填充(Flood Fill)思想及實(shí)現(xiàn)思路，但一般使用遞歸技術(shù)實(shí)現(xiàn)洪水填充，但由于地震資料的特殊性，基于遞歸技術(shù)的的經(jīng)典洪水填充算法具有很多缺點(diǎn)，難以獲得好的追蹤效果或耗費(fèi)非常高的計(jì)算機(jī)資源。
[0010]與二維圖像或三維圖像相比，地震資料的特殊性主要表現(xiàn)在兩個(gè)方面，首先，地震資料是由一系列地震波形序列組成的，相鄰的地震波形序列在橫向上常常表現(xiàn)為地震同相軸。與圖像中的區(qū)域邊界相比，地震同相軸具有一定的寬度；其次，由于地層中常常存在斷層，導(dǎo)致斷層附近的地震同相軸發(fā)生錯(cuò)斷，如果對(duì)這種情況不加以處理，層位追蹤常常會(huì)發(fā)生“串軸”現(xiàn)象，導(dǎo)致層位追蹤不符合地層沉積特征。
[0011]洪水填充(Flood fill)，又稱種子點(diǎn)填充，是一種在多維空間中確定連通區(qū)域的算法，在圖像處理中應(yīng)用非常廣泛。它從一種子點(diǎn)出發(fā)，采用一定的規(guī)則判斷種子點(diǎn)周圍的像素點(diǎn)是否與種子點(diǎn)一致或類似，以確定該點(diǎn)是否屬于種子點(diǎn)所在的填充區(qū)域。
[0012]典型的洪水填充有二連通、四連通和八連通兩種算法。二連通是指從圖像中的某一點(diǎn)出發(fā)從左、右兩個(gè)方向?qū)δ繕?biāo)區(qū)域進(jìn)行擴(kuò)展；四連通算法是指對(duì)種子點(diǎn)的四個(gè)方向即上、下、左、右的像素點(diǎn)進(jìn)行判斷以確定是否對(duì)目標(biāo)區(qū)域進(jìn)行擴(kuò)展。八連通算法與四連通算法類似，可以通過上、下、左、右、左上、右上、左下、右下這8個(gè)方向?qū)崿F(xiàn)目標(biāo)區(qū)域的擴(kuò)展。
[0013]傳統(tǒng)的洪水填充采用遞歸算法。以四連通算法為例，采用遞歸算法可使洪水填充的實(shí)現(xiàn)非常簡(jiǎn)單:
[0014]如果目標(biāo)點(diǎn)顏色等于目標(biāo)顏色，則:
[0015](I)將目標(biāo)點(diǎn)顏色設(shè)置為種子點(diǎn)顏色；
[0016](2)遞歸調(diào)用Flood-fill，填充左側(cè)點(diǎn)；
[0017](3)遞歸調(diào)用Flood-fill，填充右側(cè)點(diǎn)；
[0018](4)遞歸調(diào)用Flood-fill，填充上側(cè)點(diǎn)；
[0019](5)遞歸調(diào)用Flood-fill，填充下側(cè)點(diǎn)；
[0020]遞歸法洪水填充在填充區(qū)域較大時(shí)，遞歸層次增多，由于函數(shù)的反復(fù)調(diào)用會(huì)使操作系統(tǒng)堆棧溢出，特別是由于地震數(shù)據(jù)的數(shù)據(jù)量常常很大，會(huì)使堆棧溢出問題非常嚴(yán)重。其次，遞歸算法效率較低，多數(shù)節(jié)點(diǎn)需要被訪問四次，遞歸調(diào)用函數(shù)時(shí)間和空間占用較大。
[0021]使用隊(duì)列代替遞歸算法可避免洪水填充算法的以上問題。隊(duì)列(Queue)是一種計(jì)算機(jī)數(shù)據(jù)結(jié)構(gòu)，是一種“先進(jìn)先出”的線性表。在具體應(yīng)用中通常用鏈表或者數(shù)組來實(shí)現(xiàn)。隊(duì)列只允許在后端進(jìn)行插入操作，在前端進(jìn)行刪除操作。隊(duì)列的操作方式和堆棧類似，唯一的區(qū)別在于隊(duì)列只允許新數(shù)據(jù)在后端進(jìn)行添加。


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

[0022]本發(fā)明目地在于提供一種具有高精度、高效率、可滿足地震資料高效構(gòu)造解釋的要求的傾角傳播法層位自動(dòng)追蹤方法。
[0023]本發(fā)明通過以下步驟實(shí)現(xiàn):
[0024](I)根據(jù)地質(zhì)研究目標(biāo)在地震剖面上確定地震層位并拾取種子點(diǎn)，創(chuàng)建種子點(diǎn)庫；
[0025]所述的創(chuàng)建種子點(diǎn)庫是，首先針對(duì)不同的地質(zhì)研究目標(biāo)，選擇要追蹤的目標(biāo)地質(zhì)層位，利用測(cè)井、鉆井等資料對(duì)地質(zhì)層位進(jìn)行標(biāo)定，確定層位的地震反射相位類型，然后在地震剖面上選擇與目標(biāo)地質(zhì)層位對(duì)應(yīng)的地震同相軸，在該同相軸上拾取種子點(diǎn)，并組成種子點(diǎn)庫。
[0026]所述的種子點(diǎn)庫中的種子點(diǎn)數(shù)量和分布特征可根據(jù)地震數(shù)據(jù)的信噪比和構(gòu)造復(fù)雜度確定。如果地震數(shù)據(jù)的信噪比較低或構(gòu)造復(fù)雜，則種子點(diǎn)數(shù)量要多一些，基本原則是在一個(gè)斷塊內(nèi)要設(shè)置至少一個(gè)種子點(diǎn)。當(dāng)工區(qū)中存在斷距較大的斷層時(shí)，斷層的上下盤均應(yīng)設(shè)置種子點(diǎn)。
[0027]所述的種子點(diǎn)庫中至少要有一個(gè)種子點(diǎn)。如果種子點(diǎn)庫中沒有種子點(diǎn)存在，則跳轉(zhuǎn)到步驟(10)，停止追蹤過程。
[0028](2)對(duì)種子點(diǎn)庫中的種子點(diǎn)進(jìn)行排序，將其組成一個(gè)種子點(diǎn)隊(duì)列；
[0029]所述的隊(duì)列是指一種“先進(jìn)先出的”數(shù)據(jù)結(jié)構(gòu)，其特點(diǎn)是只允許在后端進(jìn)行插入操作，在前端進(jìn)行刪除操作。隊(duì)列的操作方式與堆棧類似，其區(qū)別在于隊(duì)列只允許新數(shù)據(jù)在后端進(jìn)行添加。
[0030]所述的種子點(diǎn)隊(duì)列是指由所有種子點(diǎn)按一定順序所組成的一個(gè)對(duì)列。
[0031]所述的對(duì)種子點(diǎn)庫中的種子點(diǎn)進(jìn)行排序，是指將所有種子點(diǎn)按拾取順序或坐標(biāo)位置進(jìn)行排序，并將其存儲(chǔ)到事先預(yù)置的可包含所有種子點(diǎn)的隊(duì)列中。
[0032](3)從種子點(diǎn)隊(duì)列中提取當(dāng)前種子點(diǎn)；
[0033]所述的提取當(dāng)前種子點(diǎn)是從種子點(diǎn)隊(duì)列中選擇最前端的種子點(diǎn)作為當(dāng)前種子點(diǎn)。
[0034](4)在地震道相似性特征控制下，利用洪水填充(Flood Fill)技術(shù)確定層位追蹤目標(biāo)點(diǎn)；
[0035]所述的確定層位追蹤目標(biāo)點(diǎn)是采用經(jīng)典的二連通或四連通洪水填充(FloodFill)算法，根據(jù)種子點(diǎn)確定目標(biāo)點(diǎn)的機(jī)制，以當(dāng)前種子點(diǎn)為中心將種子點(diǎn)周圍左和右兩個(gè)方向或上、下、左和右四個(gè)方向且相鄰的測(cè)量點(diǎn)作為目標(biāo)點(diǎn)。
[0036]所述的洪水填充(Flood Fill)技術(shù)是一種在多維空間中確定連通區(qū)域的算法，是一種經(jīng)典的圖像處理方法，主要用于圖像中的區(qū)域填充或區(qū)域顏色替換；它從一個(gè)種子點(diǎn)出發(fā)，采用一定的規(guī)則判斷種子點(diǎn)周圍的像素點(diǎn)是否與種子點(diǎn)一致或類似，以確定該點(diǎn)是否屬于種子點(diǎn)所在的填充區(qū)域。
[0037]所述的洪水填充(Flood Fill)技術(shù)具有二連通、四連通和八連通三種算法。
[0038]所述的二連通是指從圖像中的某一點(diǎn)出發(fā)從左、右兩個(gè)方向?qū)δ繕?biāo)區(qū)域進(jìn)行擴(kuò)展；二連通算法用于二維地震數(shù)據(jù)。
[0039]所述的四連通算法是指從圖像中某一點(diǎn)出發(fā)，可通過對(duì)四個(gè)方向即上、下、左、右的像素點(diǎn)進(jìn)行判斷以確定是否對(duì)目標(biāo)區(qū)域進(jìn)行擴(kuò)展；四連通算法用于三維地震數(shù)據(jù)。
[0040]所述的八連通算法是通過上、下、左、右、左上、右上、左下、右下這8個(gè)方向?qū)崿F(xiàn)目標(biāo)區(qū)域的擴(kuò)展。八連通算法用于三維地震數(shù)據(jù)。
[0041]所述的地震道相似性特征，是指在利用洪水填充(Flood Fill)技術(shù)確定層位追蹤目標(biāo)點(diǎn)時(shí)，求取過目標(biāo)點(diǎn)的地震道與過種子點(diǎn)的地震道的相似性。如果兩個(gè)地震道的相似性小于設(shè)定的門檻值，認(rèn)為種子點(diǎn)和目標(biāo)點(diǎn)之間存在斷層，種子點(diǎn)不應(yīng)擴(kuò)散到斷層的另一邊，則將該目標(biāo)點(diǎn)剔除，反之則將其作為有效目標(biāo)點(diǎn)。
[0042]所述的利用洪水填充(Flood Fill)技術(shù)確定層位追蹤目標(biāo)點(diǎn)，是利用互相關(guān)分析、方差等算法計(jì)算過目標(biāo)點(diǎn)的地震道與過種子點(diǎn)的地震道的相似性。
[0043]所述的門檻值為可接受的最小相似系數(shù)。歸一化的相似系數(shù)的值域范圍為0.0?
1.0 ;一般情況下門檻值要大于0.5，缺省門檻值為0.85。
[0044](5)分別在目標(biāo)點(diǎn)的周圍反向?qū)ふ矣行ХN子點(diǎn)；
[0045]所述的反向?qū)ふ矣行ХN子點(diǎn)，是將地層傾角視作地層阻尼系數(shù)，在一設(shè)定范圍內(nèi)沿一固定方向?qū)ふ易枘嵯禂?shù)最小的種子點(diǎn)作為有效種子點(diǎn)。
[0046]所述的反向?qū)ふ矣行ХN子點(diǎn)，是將從種子點(diǎn)到目標(biāo)點(diǎn)方向視作正方向，將從目標(biāo)點(diǎn)到種子點(diǎn)方向視作反方向，從目標(biāo)點(diǎn)開始，反方向?qū)ふ易枘嵯禂?shù)最小的種子點(diǎn)作為有效種子點(diǎn)。
[0047]所述的地層阻尼系數(shù)，是與地層傾角成正比的一種表征流體在地層中運(yùn)動(dòng)的參數(shù)，用于控制層位追蹤方向和范圍。在本專利中，地層阻尼系數(shù)是歸一化的地震同相軸傾角。
[0048]所述的利用地層阻尼系數(shù)控制層位追蹤方向和范圍，是指當(dāng)?shù)貙幼枘嵯禂?shù)大于設(shè)定的門檻值時(shí)，兩個(gè)相鄰地震道的地層傾角過大，在兩個(gè)相鄰地震道之間存在斷層，不存在有效種子點(diǎn)，此時(shí)需放棄尋找有效種子點(diǎn)，使該方向的種子點(diǎn)擴(kuò)散停止。
[0049](6)估算目標(biāo)點(diǎn)相對(duì)于有效種子點(diǎn)的地震同相軸傾角；
[0050]所述的地震同相軸傾角，在層位追蹤中，表現(xiàn)為某一地震同相軸上的種子點(diǎn)和目標(biāo)點(diǎn)之間的時(shí)差。
[0051]所述的估算地震同相軸傾角是指在設(shè)定的時(shí)窗范圍內(nèi)利用地震道的波形特征對(duì)相鄰地震道的相似性進(jìn)行逐點(diǎn)掃描，并在設(shè)定的時(shí)窗范圍內(nèi)選擇相似系數(shù)最大的時(shí)差作為該點(diǎn)的地震同相軸傾角。
[0052]所述的估算地震同相軸傾角，對(duì)于二維地震數(shù)據(jù)，沿測(cè)線方向估算一個(gè)方位角的地震同相軸傾角；對(duì)于三維地震數(shù)據(jù),沿主測(cè)線和聯(lián)絡(luò)線兩個(gè)方位估算兩個(gè)地震同相軸傾角。
[0053]所述的估算地震同相軸傾角，是指利用互相關(guān)、方差分析等數(shù)學(xué)算法沿某一傾向求取兩個(gè)地震道或多個(gè)地震道的相關(guān)系數(shù)或方差，相關(guān)系數(shù)與相似性成正比，方差與相似性成反比。
[0054](7)將層位值從有效種子點(diǎn)傳播到目標(biāo)點(diǎn)，并將其標(biāo)記為有效目標(biāo)點(diǎn)。
[0055]所述的種子點(diǎn)傳播是指，將目標(biāo)點(diǎn)相對(duì)于有效種子點(diǎn)的地震同相軸傾角所對(duì)應(yīng)的時(shí)移量與有效種子點(diǎn)的層位值相加，并將修改后的層位值視作目標(biāo)點(diǎn)的層位值。
[0056]所述的有效目標(biāo)點(diǎn)是指，如果目標(biāo)點(diǎn)相對(duì)于有效種子點(diǎn)的地震同相軸傾角存在，則目標(biāo)點(diǎn)的層位值存在，該目標(biāo)點(diǎn)為有效目標(biāo)點(diǎn)，反之則為無效目標(biāo)點(diǎn)。
[0057](8)將有效目標(biāo)點(diǎn)追加到種子點(diǎn)隊(duì)列的后端；
[0058]所述的將目標(biāo)點(diǎn)追加到種子點(diǎn)隊(duì)列的后端，是將目標(biāo)點(diǎn)作為新的種子點(diǎn)，利用洪水填充算法的特點(diǎn)，將新的種子點(diǎn)納入種子點(diǎn)擴(kuò)散過程，實(shí)現(xiàn)全工區(qū)遍歷。
[0059](9)待全部目標(biāo)點(diǎn)處理完成后，隊(duì)列前移一位，如果前移后的種子點(diǎn)隊(duì)列為空，則停止追蹤，否則返回步驟(3)繼續(xù)進(jìn)行下一個(gè)種子點(diǎn)的追蹤；
[0060](10)當(dāng)種子點(diǎn)隊(duì)列清空后，將所有追蹤到的目標(biāo)點(diǎn)的層位值與原始種子點(diǎn)組合，輸出層位，完成層位追蹤。
[0061]本發(fā)明利用逆向洪水填充算法在地層傾角控制下實(shí)現(xiàn)層位自動(dòng)追蹤，通過反向種子點(diǎn)搜索、設(shè)置地層阻尼系數(shù)和斷層阻擋等技術(shù)增強(qiáng)地震層位追蹤的連續(xù)性，提高斷層附近層位追蹤的準(zhǔn)確性，適應(yīng)高精度層位自動(dòng)追蹤和地震資料高效構(gòu)造解釋的要求。由于本發(fā)明所述的傾角傳播技術(shù)包含了二維層位傾角傳播和三維層位傾角傳播技術(shù)，因此本發(fā)明既能適應(yīng)二維工區(qū)的層位追蹤，又能適應(yīng)三維工區(qū)的層位追蹤。

【專利附圖】

【附圖說明】
[0062]結(jié)合下面的附圖及其詳細(xì)說明，本領(lǐng)域的技術(shù)人員可更加容易地理解本發(fā)明的目的、意義和特點(diǎn)，其中:
[0063]圖1為傾角傳播法層位自動(dòng)追蹤基本流程。
[0064]圖2為隊(duì)列管理示意圖?；陉?duì)列的“先進(jìn)先出”特點(diǎn)，隊(duì)列元素只能從尾端進(jìn)頂端出。
[0065]圖3為基于遞歸算法的洪水填充示意圖。左圖中的右部區(qū)域?yàn)榇畛鋮^(qū)域，在第4行第三列的方格上設(shè)置一個(gè)種子點(diǎn)。右圖為像素填充效果圖，方格中的數(shù)字表示填充順序。
[0066]圖4為基于隊(duì)列算法的洪水填充示意圖。左圖中的右部區(qū)域?yàn)榇畛鋮^(qū)域，在第4行第三列的方格上設(shè)置一個(gè)種子點(diǎn)。右圖為像素填充效果圖，方格中的數(shù)字表示填充順序。與圖3相比，填充效果一致，但填充順序不同。
[0067]圖5表示逆向洪水填充算法。逆向洪水填充算法的核心是首先將種子點(diǎn)上、下、左、右四個(gè)方向上相鄰的測(cè)量點(diǎn)作為目標(biāo)點(diǎn)(圖5a)，然后從目標(biāo)點(diǎn)出發(fā)，沿種子點(diǎn)所在的方向在一定范圍內(nèi)尋找有效種子點(diǎn)(圖5b，以目標(biāo)點(diǎn)I為例)。
[0068]圖6表示根據(jù)地震波形特征估算地震同相軸傾角。在地震剖面上，地震同相軸的傾角表現(xiàn)為相鄰兩個(gè)地震道上同一個(gè)地震同相軸的傾向與水平方向的夾角。
[0069]圖7為二維傾角傳播方法示意圖。二維傾角傳播根據(jù)地震同相軸的傾角從一個(gè)種子點(diǎn)向兩邊擴(kuò)散，既能實(shí)現(xiàn)單種子點(diǎn)傳播也能實(shí)現(xiàn)多種子點(diǎn)傳播。
[0070]圖8為三維傾角傳播方法示意圖。三維傾角傳播根據(jù)地震同相軸的傾角從一個(gè)種子點(diǎn)向四周擴(kuò)散。與二維傾角傳播方法一樣，三維傾角傳播既能實(shí)現(xiàn)單種子點(diǎn)傳播也能實(shí)現(xiàn)多種子點(diǎn)傳播。
[0071]圖9表示一個(gè)基于傾角傳播算法的二維層位自動(dòng)追蹤的實(shí)施例。在本實(shí)施例中，利用波形特征法估算地震同相軸傾角，從3個(gè)種子點(diǎn)出發(fā)，利用傾角傳播實(shí)現(xiàn)整條二維地震剖面的層位追蹤。采用傾角傳播技術(shù)可保證斷層附近層位追蹤的正確性。
[0072]圖10表示一個(gè)基于傾角傳播的三維層位自動(dòng)追蹤的實(shí)施例。與圖9中的二維層位追蹤一樣，種子點(diǎn)也為3個(gè)。在本實(shí)施例中，利用波形特征法估算地震同相軸傾角，利用傾角傳播技術(shù)實(shí)現(xiàn)三維工區(qū)的地震層位追蹤(圖10b)，同時(shí)追蹤完成后可獲得沿層傾角信息(圖10a)。

【具體實(shí)施方式】
[0073]本發(fā)明根據(jù)地震資料的特點(diǎn)和地震層位追蹤要求對(duì)經(jīng)典的“洪水填充(FloodFill)”圖像處理技術(shù)進(jìn)行改造，構(gòu)建“逆向洪水填充(Inverse Flood Fill)"算法，采用多種技術(shù)增強(qiáng)地震層位追蹤的連續(xù)性和提高斷層附近層位追蹤的準(zhǔn)確性。本發(fā)明可應(yīng)用于二維地震層位自動(dòng)追蹤和三維地震層位自動(dòng)追蹤。
[0074]根據(jù)本發(fā)明的
【發(fā)明內(nèi)容】
，本發(fā)明具有以下實(shí)施步驟(圖1):
[0075](I)根據(jù)地質(zhì)研究目標(biāo)在地震剖面上確定地震層位并拾取種子點(diǎn)，創(chuàng)建種子點(diǎn)庫；
[0076]所述的創(chuàng)建種子點(diǎn)庫是，首先針對(duì)不同的地質(zhì)研究目標(biāo)，選擇要追蹤的目標(biāo)地質(zhì)層位，利用測(cè)井、鉆井等資料對(duì)地質(zhì)層位進(jìn)行標(biāo)定，確定層位的地震反射相位類型，然后在地震剖面上選擇與目標(biāo)地質(zhì)層位對(duì)應(yīng)的地震同相軸，在該同相軸上拾取種子點(diǎn)，并組成種子點(diǎn)庫。
[0077]所述的種子點(diǎn)庫中至少要有一個(gè)種子點(diǎn)。如果種子點(diǎn)庫中沒有種子點(diǎn)存在，則跳轉(zhuǎn)到步驟(10)，停止追蹤過程。
[0078]根據(jù)地震工區(qū)類型的不同，種子點(diǎn)的坐標(biāo)類型有所不同。對(duì)于圖9所表示的二維地震測(cè)線層位自動(dòng)追蹤實(shí)施例，僅需要地震道號(hào)和時(shí)間值即可確定種子點(diǎn)的空間位置；對(duì)于圖10所表示的三維地震工區(qū)，則需要測(cè)線號(hào)、地震道號(hào)和時(shí)間值來確定種子點(diǎn)的空間位置。
[0079](2)對(duì)種子點(diǎn)庫中的種子點(diǎn)進(jìn)行排序，將其組成一個(gè)種子點(diǎn)隊(duì)列；
[0080]所述的隊(duì)列是指，隊(duì)列是一種“先進(jìn)先出的”數(shù)據(jù)結(jié)構(gòu)，其特點(diǎn)是只允許在后端進(jìn)行插入操作，在前端進(jìn)行刪除操作，如圖2所示。隊(duì)列的操作方式與堆棧類似，其區(qū)別在于隊(duì)列只允許新數(shù)據(jù)在后端進(jìn)行添加。
[0081]在將種子點(diǎn)組成種子點(diǎn)隊(duì)列時(shí)，將所有種子點(diǎn)按拾取順序或坐標(biāo)位置進(jìn)行排序，并將其存儲(chǔ)到事先預(yù)置的可包含所有種子點(diǎn)的隊(duì)列中。
[0082](3)從種子點(diǎn)隊(duì)列中提取當(dāng)前種子點(diǎn)；
[0083]所述的提取當(dāng)前種子點(diǎn)是從種子點(diǎn)隊(duì)列中選擇最前端的種子點(diǎn)作為當(dāng)前種子點(diǎn)，如圖2中的元素ei。
[0084](4)在地震道相似性特征控制下，利用洪水填充(Flood Fill)技術(shù)確定層位追蹤目標(biāo)點(diǎn)；
[0085]所述的利用洪水填充技術(shù)確定層位追蹤目標(biāo)點(diǎn)是采用經(jīng)典的二連通、四連通或八連通洪水填充(Flood Fill)算法。根據(jù)地震工區(qū)類型的不同，所采用的洪水填充算法也不同。對(duì)于二維地震測(cè)線僅采用二連通算法，對(duì)于三維地震工區(qū)，可采用四連通和八連通算法，一般常用四連通算法。
[0086]當(dāng)采用四連通算法時(shí)，根據(jù)種子點(diǎn)確定目標(biāo)點(diǎn)的機(jī)制，以當(dāng)前種子點(diǎn)為中心將種子點(diǎn)周圍左和右兩個(gè)方向或上、下、左和右四個(gè)方向且相鄰的測(cè)量點(diǎn)作為目標(biāo)點(diǎn)。
[0087]在實(shí)現(xiàn)洪水填充時(shí)，可采用遞歸算法或隊(duì)列算法實(shí)現(xiàn)種子點(diǎn)的擴(kuò)散和區(qū)域填充。圖3為遞歸算法示意圖，圖4為隊(duì)列算法示意圖。在圖3和圖4中，左圖中的右部區(qū)域?yàn)榇畛鋮^(qū)域，在第4行第三列的方格上設(shè)置一個(gè)種子點(diǎn)。右圖為像素填充效果圖，方格中的數(shù)字表示填充順序。無論采用哪一種算法，洪水填充的基本原理是，從一個(gè)種子點(diǎn)出發(fā)，采用一定的規(guī)則判斷種子點(diǎn)周圍的像素點(diǎn)是否與種子點(diǎn)一致或類似，以確定該點(diǎn)是否屬于種子點(diǎn)所在的填充區(qū)域。
[0088]在將洪水填充技術(shù)用于層位自動(dòng)追蹤時(shí)，要用地震道相似性特征代替圖像區(qū)域填充中的像素差異性特征來控制種子點(diǎn)的擴(kuò)散方向和范圍。如果兩個(gè)地震道的相似性小于設(shè)定的門檻值，認(rèn)為種子點(diǎn)和目標(biāo)點(diǎn)之間存在斷層，種子點(diǎn)不應(yīng)擴(kuò)散到斷層的另一邊，將目標(biāo)點(diǎn)剔除，反之則將其作為有效目標(biāo)點(diǎn)。
[0089]地震道相似性是利用互相關(guān)分析、方差等算法對(duì)過目標(biāo)點(diǎn)的地震道與過種子點(diǎn)的地震道計(jì)算得到的，相似性門檻值為可接受的最小相似系數(shù)。歸一化的相似系數(shù)的值域范圍為0.0?1.0 ;—般情況下門檻值要大于0.5，缺省值為0.85。
[0090](5)分別在目標(biāo)點(diǎn)的周圍反向?qū)ふ矣行ХN子點(diǎn)；
[0091]所述的反向?qū)ふ矣行ХN子點(diǎn)，是將地層傾角視作地層阻尼系數(shù)，在設(shè)定范圍內(nèi)沿某一個(gè)固定方向?qū)ふ易枘嵯禂?shù)最小的種子點(diǎn)作為有效種子點(diǎn)，如圖5所示，首先將種子點(diǎn)上、下、左、右四個(gè)方向上相鄰的測(cè)量點(diǎn)作為目標(biāo)點(diǎn)(圖5a)，然后從目標(biāo)點(diǎn)出發(fā)，沿種子點(diǎn)所在的方向在一定范圍內(nèi)尋找有效種子點(diǎn)(圖5b，以目標(biāo)點(diǎn)I為例)。
[0092]在本專利中，所述的地層阻尼系數(shù)，是與地層傾角成正比的一種表征流體在地層中運(yùn)動(dòng)的參數(shù)，用于控制層位追蹤方向和范圍。當(dāng)?shù)貙幼枘嵯禂?shù)大于設(shè)定的門檻值時(shí)，兩個(gè)相鄰地震道的地層傾角過大，在兩個(gè)相鄰地震道之間存在斷層，不存在有效種子點(diǎn)，此時(shí)需放棄尋找有效種子點(diǎn)，使該方向的種子點(diǎn)擴(kuò)散停止。
[0093](6)估算目標(biāo)點(diǎn)相對(duì)于有效種子點(diǎn)的地震同相軸傾角；
[0094]所述的地震同相軸傾角，在層位追蹤中，表現(xiàn)為某一地震同相軸上的種子點(diǎn)和目標(biāo)點(diǎn)之間的時(shí)差，如圖6所示。對(duì)于二維地震數(shù)據(jù)，沿測(cè)線方向估算一個(gè)方位角的地震同相軸傾角；對(duì)于三維地震數(shù)據(jù)，沿主測(cè)線和聯(lián)絡(luò)線兩個(gè)方向估算兩個(gè)地震同相軸傾角。
[0095]所述的估算地震同相軸傾角是指在設(shè)定的時(shí)窗范圍內(nèi)利用地震道的波形特征對(duì)相鄰地震道的相似性進(jìn)行逐點(diǎn)掃描，并在設(shè)定的時(shí)窗范圍內(nèi)選擇相似系數(shù)最大的時(shí)差作為該點(diǎn)的地震同相軸傾角。
[0096]在本專利的具體實(shí)施步驟中，該步驟可簡(jiǎn)稱為“傾角估算”。傾角估算的精度會(huì)影響層位追蹤的準(zhǔn)確性。但由于地震同相軸傾角估算方法很多，且為使本專利所述的傾角傳播能夠適應(yīng)多種傾角估算方法，因此未對(duì)特定的傾角估算方法進(jìn)行描述，也未限制使用任何傾角估算方法。實(shí)際上，目前常用的一些傾角估算技術(shù)，如波形特征法、相關(guān)法、最大能量法、平面波分解法等，都可以與傾角傳播法配合形成可滿足不同類型、不同品質(zhì)地震資料的層位自動(dòng)追蹤軟件產(chǎn)品。
[0097](7)將層位值從有效種子點(diǎn)傳播到目標(biāo)點(diǎn)，并將其標(biāo)記為有效目標(biāo)點(diǎn)。
[0098]所述的種子點(diǎn)傳播是指，將目標(biāo)點(diǎn)相對(duì)于有效種子點(diǎn)的地震同相軸傾角所對(duì)應(yīng)的時(shí)移量與有效種子點(diǎn)的層位值相加，并將修改后的層位值視作目標(biāo)點(diǎn)的層位值，種子點(diǎn)傳播是依賴地震同相軸傾角而實(shí)現(xiàn)的，如圖7和圖8所示。其中，圖7為二維傾角傳播示意圖，圖8為三維傾角傳播示意圖。
[0099]在本專利的具體實(shí)施步驟中，該步驟可簡(jiǎn)稱為“傾角傳播”。該步驟是本專利的核心，根據(jù)地震同相軸傾角信息控制種子點(diǎn)的傳播，強(qiáng)調(diào)傾角信息特別是不同方位的傾角信息在層位自動(dòng)追蹤過程中的重要性，可將不同類型的層位自動(dòng)追蹤技術(shù)統(tǒng)一起來，并根據(jù)傾角估算技術(shù)的特點(diǎn)進(jìn)一步提聞層位自動(dòng)追蹤的精度。
[0100](8)將有效目標(biāo)點(diǎn)追加到種子點(diǎn)隊(duì)列的后端；
[0101]所述的將目標(biāo)點(diǎn)追加到種子點(diǎn)隊(duì)列的后端，是將目標(biāo)點(diǎn)作為新的種子點(diǎn)，利用洪水填充算法的特點(diǎn)，將新的種子點(diǎn)納入種子點(diǎn)擴(kuò)散過程，實(shí)現(xiàn)全工區(qū)遍歷，如圖2中的元素ei所示。
[0102](9)待全部目標(biāo)點(diǎn)處理完成后，隊(duì)列前移一位，如果前移后的種子點(diǎn)隊(duì)列為空，則停止追蹤，否則返回步驟3)繼續(xù)進(jìn)行下一個(gè)種子點(diǎn)的追蹤；
[0103](10)當(dāng)種子點(diǎn)隊(duì)列清空后，將所有追蹤到的目標(biāo)點(diǎn)的層位值與原始種子點(diǎn)組合，輸出層位，完成層位追蹤，如圖9和圖10所示。其中，圖9為一個(gè)基于傾角傳播算法的二維層位自動(dòng)追蹤的實(shí)施例，在本實(shí)施例中，利用波形特征法估算地震同相軸傾角，從3個(gè)種子點(diǎn)出發(fā)，利用傾角傳播實(shí)現(xiàn)整條二維地震剖面的層位追蹤。采用傾角傳播技術(shù)可保證斷層附近層位追蹤的正確性。圖10為一個(gè)基于傾角傳播的三維層位自動(dòng)追蹤的實(shí)施例。與圖9中的二維層位追蹤一樣，種子點(diǎn)也為3個(gè)。在本實(shí)施例中，利用波形特征法估算地震同相軸傾角，利用傾角傳播技術(shù)實(shí)現(xiàn)三維工區(qū)的地震層位追蹤(圖10b)，同時(shí)追蹤完成后可獲得沿層傾角信息(圖10a)。以上兩個(gè)實(shí)施例表明，傾角傳播法層位自動(dòng)追蹤方法具有高精度、聞效率等特點(diǎn)，可滿足地震資料聞效構(gòu)造解釋的要求。
【權(quán)利要求】
1.一種傾角傳播法層位自動(dòng)追蹤方法，特點(diǎn)是包括以下步驟: (1)根據(jù)地質(zhì)研究目標(biāo)在地震剖面上確定地震層位并拾取種子點(diǎn)，創(chuàng)建種子點(diǎn)庫；所述的種子點(diǎn)庫中至少要有一個(gè)種子點(diǎn)，如果種子點(diǎn)庫中沒有種子點(diǎn)存在，則跳轉(zhuǎn)到步驟(10)，停止追蹤過程； (2)對(duì)種子點(diǎn)庫中的種子點(diǎn)進(jìn)行排序，將其組成一個(gè)種子點(diǎn)隊(duì)列； (3)從種子點(diǎn)隊(duì)列中提取當(dāng)前種子點(diǎn)； (4)在地震道相似性特征控制下，利用洪水填充(FloodFill)技術(shù)確定層位追蹤目標(biāo)占.(5)分別在目標(biāo)點(diǎn)的周圍反向?qū)ふ矣行ХN子點(diǎn)； (6)估算目標(biāo)點(diǎn)相對(duì)于有效種子點(diǎn)的地震同相軸傾角； (7)將層位值從有效種子點(diǎn)傳播到目標(biāo)點(diǎn)，并將其標(biāo)記為有效目標(biāo)點(diǎn)； (8)將有效目標(biāo)點(diǎn)追加到種子點(diǎn)隊(duì)列的后端； (9)待全部目標(biāo)點(diǎn)處理完成后，隊(duì)列前移一位，如果前移后的種子點(diǎn)隊(duì)列為空，則停止追蹤，否則返回步驟(3)繼續(xù)進(jìn)行下一個(gè)種子點(diǎn)的追蹤； (10)當(dāng)種子點(diǎn)隊(duì)列清空后，將所有追蹤到的目標(biāo)點(diǎn)的層位值與原始種子點(diǎn)組合，輸出層位，完成層位追蹤。
2.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(I)所述的創(chuàng)建種子點(diǎn)庫是，首先針對(duì)不同的地質(zhì)研究目標(biāo)，選擇要追蹤的目標(biāo)地質(zhì)層位，利用測(cè)井、鉆井等資料對(duì)地質(zhì)層位進(jìn)行標(biāo)定，在地震剖面上確定對(duì)應(yīng)的地震反射同相軸，并根據(jù)地震同相軸的相位類型定義地震層位，然后在該地震層位對(duì)應(yīng)的地震同相軸上拾取種子點(diǎn)，并組成種子點(diǎn)庫。
3.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(I)所述的種子點(diǎn)庫中的種子點(diǎn)數(shù)量和分布特征可根據(jù)地震數(shù)據(jù)的信噪比和構(gòu)造復(fù)雜度確定。如果地震數(shù)據(jù)的信噪比較低或構(gòu)造復(fù)雜，則種子點(diǎn)數(shù)量要多一些，基本原則是在一個(gè)斷塊內(nèi)要設(shè)置至少一個(gè)種子點(diǎn)。當(dāng)工區(qū)中存在斷距較大的斷層時(shí)，斷層的上下盤均應(yīng)設(shè)置種子點(diǎn)。
4.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(2)所述的隊(duì)列是指一種先進(jìn)先出的數(shù)據(jù)結(jié)構(gòu)，只允許在后端進(jìn)行插入操作，在前端進(jìn)行刪除。
5.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(2)所述的種子點(diǎn)隊(duì)列是指由所有種子點(diǎn)按一定順序所組成的一個(gè)對(duì)列。
6.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(2)所述的對(duì)種子點(diǎn)庫中的種子點(diǎn)進(jìn)行排序，是指將所有種子點(diǎn)按拾取順序或坐標(biāo)位置進(jìn)行排序，并將其存儲(chǔ)到事先預(yù)置的可包含所有種子點(diǎn)的隊(duì)列中。
7.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(3)所述的提取當(dāng)前種子點(diǎn)是從種子點(diǎn)隊(duì)列中選擇最前端的種子點(diǎn)作為當(dāng)前種子點(diǎn)。
8.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(4)所述的確定層位追蹤目標(biāo)點(diǎn)是采用經(jīng)典的二連通或四連通洪水填充算法，根據(jù)種子點(diǎn)確定目標(biāo)點(diǎn)的機(jī)制，以當(dāng)前種子點(diǎn)為中心將種子點(diǎn)周圍左和右兩個(gè)方向或上、下、左和右四個(gè)方向且相鄰的測(cè)量點(diǎn)作為目標(biāo)點(diǎn)。
9.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(4)所述的洪水填充技術(shù)是一種在多維空間中確定連通區(qū)域的算法，是一種經(jīng)典的圖像處理方法，主要用于圖像中的區(qū)域填充或區(qū)域顏色替換；它從一個(gè)種子點(diǎn)出發(fā)，采用一定的規(guī)則判斷種子點(diǎn)周圍的像素點(diǎn)是否與種子點(diǎn)一致或類似，以確定該點(diǎn)是否屬于種子點(diǎn)所在的填充區(qū)域，實(shí)現(xiàn)圖像的區(qū)域填充或擴(kuò)展。
10.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(4)所述的洪水填充技術(shù)具有二連通、四連通和八連通三種算法。
11.根據(jù)權(quán)利要求10的方法，特點(diǎn)是所述的二連通是指從圖像中的某一點(diǎn)出發(fā)從左、右兩個(gè)方向?qū)δ繕?biāo)區(qū)域進(jìn)行擴(kuò)展；二連通算法用于二維地震數(shù)據(jù)。
12.根據(jù)權(quán)利要求10的方法，特點(diǎn)是所述的四連通算法是指從圖像中某一點(diǎn)出發(fā)，可通過對(duì)四個(gè)方向即上、下、左、右的像素點(diǎn)進(jìn)行判斷以確定是否對(duì)目標(biāo)區(qū)域進(jìn)行擴(kuò)展；四連通算法用于三維地震數(shù)據(jù)。
13.根據(jù)權(quán)利要求10的方法，特點(diǎn)是所述的八連通算法是通過上、下、左、右、左上、右上、左下、右下8個(gè)方向?qū)崿F(xiàn)目標(biāo)區(qū)域的擴(kuò)展；八連通算法用于三維地震數(shù)據(jù)。
14.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(4)所述的地震道相似性特征，是指在利用洪水填充技術(shù)確定層位追蹤目標(biāo)點(diǎn)時(shí)，求取過目標(biāo)點(diǎn)的地震道與過種子點(diǎn)的地震道的相似性；如果兩個(gè)地震道的相似性小于設(shè)定的門檻值，認(rèn)為種子點(diǎn)和目標(biāo)點(diǎn)之間存在斷層，種子點(diǎn)不應(yīng)擴(kuò)散到斷層的另一邊，則將該目標(biāo)點(diǎn)剔除，反之則將其作為有效目標(biāo)點(diǎn)。
15.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(4)所述的利用洪水填充技術(shù)確定層位追蹤目標(biāo)點(diǎn)，是利用互相關(guān)分析、方差等算法計(jì)算過目標(biāo)點(diǎn)的地震道與過種子點(diǎn)的地震道的相似性。
16.根據(jù)權(quán)利要求14的方法，特點(diǎn)是所述的門檻值為可接受的最小相似系數(shù)；歸一化的相似系數(shù)的值域范圍為0.0?1.0 ;—般情況下門檻值要大于0.5，缺省門檻值為0.85。
17.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(5)所述的反向?qū)ふ矣行ХN子點(diǎn)，是將地層傾角視作地層阻尼系數(shù)，在一設(shè)定范圍內(nèi)沿一固定方向?qū)ふ易枘嵯禂?shù)最小的種子點(diǎn)作為有效種子點(diǎn)。
18.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(5)所述的反向?qū)ふ矣行ХN子點(diǎn)，是將從種子點(diǎn)到目標(biāo)點(diǎn)方向視作正方向，將從目標(biāo)點(diǎn)到種子點(diǎn)方向視作反方向，從目標(biāo)點(diǎn)開始，反方向?qū)ふ易枘嵯禂?shù)最小的種子點(diǎn)作為有效種子點(diǎn)。
19.根據(jù)權(quán)利要求17的方法，特點(diǎn)是所述的地層阻尼系數(shù)，是與地層傾角成正比的一種表征流體在地層中運(yùn)動(dòng)的參數(shù)，用于控制層位追蹤方向和范圍。在本專利中，地層阻尼系數(shù)是歸一化的地震同相軸傾角；地震同相軸傾角越大，地層阻尼系數(shù)越大，反之則越小。
20.根據(jù)權(quán)利要求19的方法，特點(diǎn)是所述的利用地層阻尼系數(shù)控制層位追蹤方向和范圍，是指當(dāng)?shù)貙幼枘嵯禂?shù)大于設(shè)定的門檻值時(shí)，兩個(gè)相鄰地震道的地層傾角過大，在兩個(gè)相鄰地震道之間存在斷層，不存在有效種子點(diǎn)，此時(shí)需放棄尋找有效種子點(diǎn)，使該方向的種子點(diǎn)擴(kuò)散停止。
21.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(6)所述的地震同相軸傾角，在層位追蹤中，表現(xiàn)為某一地震同相軸上的種子點(diǎn)和目標(biāo)點(diǎn)之間的時(shí)差。
22.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(6)所述的估算地震同相軸傾角是指在設(shè)定的時(shí)窗范圍內(nèi)利用地震道的波形特征對(duì)相鄰地震道的相似性進(jìn)行逐點(diǎn)掃描，并在設(shè)定的時(shí)窗范圍內(nèi)選擇相似系數(shù)最大的時(shí)差作為該點(diǎn)的地震同相軸傾角。
23.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(6)所述的估算地震同相軸傾角，對(duì)于二維地震數(shù)據(jù)，沿測(cè)線方向估算一個(gè)方位角的地震同相軸傾角；對(duì)于三維地震數(shù)據(jù)，沿主測(cè)線和聯(lián)絡(luò)線兩個(gè)方向估算兩個(gè)地震同相軸傾角。
24.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(6)所述的估算地震同相軸傾角，是指利用互相關(guān)、方差分析等數(shù)學(xué)算法沿某一傾向求取兩個(gè)地震道或多個(gè)地震道的相關(guān)系數(shù)或方差，相關(guān)系數(shù)與相似性成正比，方差與相似性成反比。
25.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(7)所述的種子點(diǎn)傳播是指，將目標(biāo)點(diǎn)相對(duì)于有效種子點(diǎn)的地震同相軸傾角所對(duì)應(yīng)的時(shí)移量與有效種子點(diǎn)的層位值相加，并將修改后的層位值視作目標(biāo)點(diǎn)的層位值。
26.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(7)所述的有效目標(biāo)點(diǎn)是指，如果目標(biāo)點(diǎn)相對(duì)于有效種子點(diǎn)的地震同相軸傾角存在，則目標(biāo)點(diǎn)的層位值存在，該目標(biāo)點(diǎn)為有效目標(biāo)點(diǎn)，反之則為無效目標(biāo)點(diǎn)。
27.根據(jù)權(quán)利要求1的方法，特點(diǎn)是步驟(8)所述的將目標(biāo)點(diǎn)追加到種子點(diǎn)隊(duì)列的后端，是將目標(biāo)點(diǎn)作為新的種子點(diǎn)，利用洪水填充算法的特點(diǎn)，將新的種子點(diǎn)納入種子點(diǎn)擴(kuò)散過程，實(shí)現(xiàn)全工區(qū)遍歷。
【文檔編號(hào)】G01V1/28GK104375175SQ201310356318
【公開日】2015年2月25日 申請(qǐng)日期:2013年8月15日 優(yōu)先權(quán)日:2013年8月15日
【發(fā)明者】陳茂山, 詹仕凡, 于海生, 郝彥國(guó), 趙海珍, 左紅光 申請(qǐng)人:中國(guó)石油天然氣集團(tuán)公司, 中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司