一種地層傾角檢測方法
【技術(shù)領域】
[0001] 本發(fā)明涉及油氣勘探技術(shù)領域，具體地說，涉及一種地層傾角檢測方法。
【背景技術(shù)】
[0002] 地震數(shù)據(jù)處理中深度域偏移剖面的傾角、方位角是三維地下地質(zhì)體的重要屬性。 作為中間結(jié)果，它可以用來進行深度域網(wǎng)格層析反演，提高深度域速度模型精度。
[0003] 現(xiàn)有的地層傾角檢測方法包括傾斜疊加法、平面波分解法、復地震道分析法、地震 波相干法、梯度結(jié)構(gòu)張量法和給予加權(quán)結(jié)構(gòu)張量法等方法。其中，利用傾斜疊加法能夠計算 二維地震剖面同相軸的傾角；利用平面波分解法、復地震道分析法和地震波相干法能夠獲 取得到局部地層傾角；利用梯度結(jié)構(gòu)張量法能夠?qū)θS地震數(shù)據(jù)的梯度向量進行分析，從 而得到局部地層傾角和方位角；利用給予加權(quán)結(jié)構(gòu)張量法能夠提高傾角估計的一致性。
[0004] 但是，現(xiàn)有的地層傾角檢測方法在對低信噪比數(shù)據(jù)時的地層傾角進行檢測時，會 出現(xiàn)檢測結(jié)果不準確的問題。同時，當?shù)貙颖容^稀疏時，現(xiàn)有的地層傾角檢測方法無法檢測 缺乏反射波同相軸的地層的地層傾角。
[0005] 基于上述情況，亟需一種準確、可靠地進行地層傾角檢測的方法。

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

[0006] 為解決上述問題，本發(fā)明提供了一種地層傾角檢測方法，所述方法包括：
[0007] 地震波形數(shù)據(jù)獲取步驟，根據(jù)當前道分析點和各個預設掃描傾角，確定對應于所 述各個預設掃描傾角的臨近道分析點，根據(jù)所述當前道分析點和臨近道分析點，分別獲取 當前道的地震波形數(shù)據(jù)和對應于所述各個預設掃描傾角的臨近道的地震波形數(shù)據(jù)；
[0008] 相似能量譜計算步驟，根據(jù)所述當前道的地震波形數(shù)據(jù)和對應于所述各個預設掃 描傾角的臨近道的地震波形數(shù)據(jù)，確定所述各個預設掃描傾角的相似能量譜；
[0009] 地層傾角確定步驟，從各個相似能量譜中提取最大相似能量譜，并將所述最大相 似能量譜所對應的預設掃描傾角作為地層傾角。
[0010] 根據(jù)本發(fā)明的一個實施例，在所述地震波形數(shù)據(jù)獲取步驟中，
[0011] 在當前道中，獲取以所述當前道分析點為中心的預設深度窗口的地震波形數(shù)據(jù)， 作為所述當前道的地震波形數(shù)據(jù)；
[0012] 在各個臨近道中，獲取以相應臨近道分析點為中心的預設深度窗口的地震波形數(shù) 據(jù)，作為所述臨近道的地震波形數(shù)據(jù)。
[0013] 根據(jù)本發(fā)明的一個實施例，所述當前道兩側(cè)的臨近道的數(shù)量相等。
[0014] 根據(jù)本發(fā)明的一個實施例，所述相似能量譜計算步驟包括：
[0015] 根據(jù)所述當前道的地震波形數(shù)據(jù)和對應于所述各個預設掃描傾角的臨近道的地 震波形數(shù)據(jù)，計算對應于所述各個預設掃描傾角的各個臨近道的地震波形數(shù)據(jù)的相關(guān)值；
[0016] 根據(jù)所述各個臨近道的地震波形數(shù)據(jù)的相關(guān)值，計算所述各個預設掃描傾角的相 似能量譜。
[0017] 根據(jù)本發(fā)明的一個實施例，根據(jù)如下公式計算所述各個臨近道的地震波形數(shù)據(jù)的 相關(guān)值：
[0019] 其中，Cini表示對應于掃描傾角m的第i個臨近道的地震波形數(shù)據(jù)的相關(guān)值，Wd表 示預設深度窗口所包含的分析點的總數(shù)，4表示所述當前道的預設深度窗口中第k個分析 點，b ink表示第i個臨近道中對應于傾角m的預設深度窗口中的第k個分析點，Dak和Dbink 分別表示分析點ak和binik的地震波形數(shù)據(jù)。
[0020] 根據(jù)本發(fā)明的一個實施例，根據(jù)如下公式計算所述各個預設掃描傾角的相似能量 譜：
[0022] 其中，Sm表示掃描傾角m的相似能量譜，η表示臨近道的總數(shù)。
[0023] 根據(jù)本發(fā)明的一個實施例，在所述地層傾角確定步驟中，還從各個相似能量譜中 提取與所述最大相似能量譜臨近的若干相似能量譜，
[0024] 其中，根據(jù)所述最大相似能量譜和所述若干臨近的能量譜確定地層傾角插值曲 線，并根據(jù)所述地層傾角插值曲線確定地層傾角。
[0025] 根據(jù)本發(fā)明的一個實施例，從各個相似能量譜中提取與所述最大相似能量譜相鄰 的兩個相似能量譜。
[0026] 根據(jù)本發(fā)明的一個實施例，在所述地層傾角確定步驟中，
[0027] 計算所述地層傾角插值曲線的一階導數(shù)；
[0028] 將所述一階導數(shù)等于零時所對應的傾角作為地層傾角。
[0029] 根據(jù)本發(fā)明的一個實施例，根據(jù)如下公式計算所述地層傾角：
[0031] 其中，X表示地層傾角，Y1表示最大相似能量譜，X1表示最大相似能量譜對應的預 設掃描傾角，^和y 2表示與所述最大相似能量譜相鄰的兩個相似能量譜。
[0032] 本發(fā)明所提供的地層傾角檢測方法通過利用預設窗口長度的深度窗口來其他區(qū) 相應的地震波形數(shù)據(jù)，這樣相較于現(xiàn)有的地層傾角檢測方法，能夠明顯減少同相軸稀疏性 的影響，能夠得到更為穩(wěn)定、精度更高的傾角剖面。
[0033] 同時，通過相關(guān)屬性的同相疊加，使得相較于現(xiàn)有的地層傾角檢測方法，本方法對 數(shù)據(jù)信噪比的依賴性更低，并具有更好的通用性。本方法流程簡單、易于實現(xiàn)。同時由于在 進行地層傾角檢測時，是針對地震剖面上行的一個具體的點來計算的，與其他點之間沒有 聯(lián)系，因此本方法的計算密度高、并行性好，適合大規(guī)模并行計算。
[0034] 此外，本方法還可以根據(jù)最大相似能量譜及其取值鄰近的相似能量譜來確定地層 傾角插值曲線，并根據(jù)該插值曲線來確定地層傾角。這樣所得到的檢測結(jié)果能夠更加接近 真實值。而通過地層傾角插值曲線能夠得到連續(xù)的地層傾角信息，還可以作為輸入信息進 行深度域?qū)游鏊俣确囱?，具有更加廣泛的應用價值。
[0035] 本發(fā)明的其它特征和優(yōu)點將在隨后的說明書中闡述，并且，部分地從說明書中變 得顯而易見，或者通過實施本發(fā)明而了解。本發(fā)明的目的和其他優(yōu)點可通過在說明書、權(quán)利 要求書以及附圖中所特別指出的結(jié)構(gòu)來實現(xiàn)和獲得。
【附圖說明】
[0036] 為了更清楚地說明本發(fā)明實施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對實施例或現(xiàn) 有技術(shù)描述中所需要的附圖做簡單的介紹：
[0037] 圖1是根據(jù)本發(fā)明一個實施例的地層傾角檢測方法的流程圖；
[0038] 圖2是根據(jù)本發(fā)明另一個實施例的地層傾角檢測方法的流程圖；
[0039] 圖3是根據(jù)本發(fā)明一個實施例的輸入的地震剖面圖；
[0040] 圖4是利用傾斜疊加法檢測得到的圖3的地層傾角剖面圖；
[0041] 圖5是根據(jù)本發(fā)明一個實施例的地層傾角檢測方法得到的圖3的地層傾角剖面 圖。
【具體實施方式】
[0042] 以下將結(jié)合附圖及實施例來詳細說明本發(fā)明的實施方式，借此對本發(fā)明如何應用 技術(shù)手段來解決技術(shù)問題，并達成技術(shù)效果的實現(xiàn)過程能充分理解并據(jù)以實施。需要說明 的是，只要不構(gòu)成沖突，本發(fā)明中的各個實施例以及各實施例中的各個特征可以相互結(jié)合， 所形成的技術(shù)方案均在本發(fā)明的保護范圍之內(nèi)。
[0043] 同時，在以下說明中，出于解釋的目的而闡述了許多具體細節(jié)，以提供對本發(fā)明實 施例的徹底理解。然而，對本領域的技術(shù)人員來說顯而易見的是，本發(fā)明可以不用這里的具 體細節(jié)或者所描述的特定方式來實施。
[0044] 另外，在附圖的流程圖示出的步驟可以在諸如一組計算機可執(zhí)行指令的計算機系 統(tǒng)中執(zhí)行，并且，雖然在流程圖中示出了邏輯順序，但是在某些情況下，可以以不同于此處 的順序執(zhí)行所示出或描述的步驟。
[0045] 實施例一:
[0046] 針對現(xiàn)有的地層傾角檢測方法所存在的缺陷，本實施例基于窗口內(nèi)地震波形的相 關(guān)值進行傾斜疊加分析，并根據(jù)疊加得到的能量最大值來確定最佳的地層傾角。
[0047] 圖1示出了本實施例所提供的地層傾角檢測方法的流程圖。
[0048] 本實施例首先在地震波形數(shù)據(jù)獲取步驟中獲取當前道和對應于各個預設掃描傾 角的臨近道的地震波形數(shù)據(jù)。具體地，如圖1所示，首先在步驟SlOl中根據(jù)當前道分析點 和各個預設掃描傾角，確定對應于各個預設掃描傾角的臨近道分析點。臨近道分布在當前 道的兩側(cè)且靠近當前道，也可以說是以當前道為中心向兩側(cè)延伸預設道數(shù)的地震道，優(yōu)選 地分布在當前道兩側(cè)的臨近道的數(shù)量相等。當然，在本發(fā)明的其他實施例中，分布在當前到 兩側(cè)的臨近道的數(shù)量也可以不等，本發(fā)明不限于此。
[0049] 當前道分析點是已知的，那么過當前道分析點以預設掃描傾角構(gòu)建一條直線，而 這條直線與各個臨近道的交點便是各個臨近道的臨近道分析點。對于不同的預設掃描傾 角，各個臨近道的臨近道分析點是不同的。所以在步驟S102中，根據(jù)當前道分析點和各個 臨近道分析點，便可以分別獲取當前道的地震波形數(shù)據(jù)和對應于各個預設掃描傾角的臨近 道的地震波形數(shù)據(jù)。
[0050] 本實施例中，對于當前道，獲取以當前道分析點為中心的預設深度窗口的地震波 形數(shù)據(jù)，并將獲取的地震波形數(shù)據(jù)作為當前道的地震波形數(shù)據(jù)。同樣，對于各個臨近道，獲 取以相應臨近道分析點為中心的預設深度窗口的地震波形數(shù)據(jù)，并將獲取的地震波形數(shù)據(jù) 作為相應臨近道的地震波形數(shù)據(jù)。
[0051] 需要說明的是，在本發(fā)明的其他實施例中，當前道和各個臨近道的地震波形數(shù)據(jù) 還可以通過其他合理方式獲得，例如對于當前道，將以當前道分析點為起點的預設深度窗 口的地震波形數(shù)據(jù)作為當前道的地震波形數(shù)據(jù)，而各個臨近道的地震波形數(shù)據(jù)也可以通過 同樣方式獲得，但本發(fā)明不限于此。
[0052] 本實施例所提供的地層檢測方法通過利用預設窗口長度的深度窗口來其他區(qū)相 應的地震波形數(shù)據(jù)，這樣相較于現(xiàn)有的地層傾角檢測方法，能夠明顯減少同相軸稀疏性的 影響，從而得到更為穩(wěn)定、精度更高的傾角剖面。
[0053] 獲得當前道和各個臨近道的地震波形數(shù)據(jù)后，本實施例在相似能量譜計算步驟 中，根據(jù)當前道的地震波形數(shù)據(jù)和對應于各個預設掃描傾角的臨近道的地震波形數(shù)據(jù)，來 確定各個預設掃描傾角的相似能量譜。
[0054] 如圖1所示，具體地，本實施例在步驟S103中，根據(jù)當前道的地震波形數(shù)據(jù)和對應 于各個預設掃描傾角的臨近道的地震波形數(shù)據(jù)，計算對應于各個預設掃描傾角的各個臨近 道的地震波形數(shù)據(jù)的相關(guān)值。本實施例中，根據(jù)如下公式計算各個臨近道的地震波形數(shù)據(jù) 的相關(guān)值：
[0056] 其中，Cini表示對應于傾角m的第i個臨近道的地震波形數(shù)據(jù)的相關(guān)值，Wd表示預 設深度窗口