一種垂直地震數(shù)據(jù)成像方法及裝置的制造方法
【專利摘要】本申請實施例公開了一種垂直地震數(shù)據(jù)成像方法及裝置，所述方法包括：獲取采集到的原始垂直地震數(shù)據(jù)，對所述采集到的原始垂直地震數(shù)據(jù)進行預處理；采用預設的參數(shù)提取方法提取所述預處理后的垂直地震數(shù)據(jù)的參數(shù)對應的數(shù)據(jù)；所述參數(shù)包括：振幅、頻率和相位；采用第一成像方法將所述提取的每一種參數(shù)對應的數(shù)據(jù)分別轉換為參數(shù)數(shù)據(jù)體；根據(jù)所述參數(shù)數(shù)據(jù)體，利用第二成像方法得到垂直地震數(shù)據(jù)成像。本申請實施例提供的一種垂直地震數(shù)據(jù)成像方法及裝置，可以提高垂直地震數(shù)據(jù)的頻率和相位的保真性。
【專利說明】
一種垂直地震數(shù)據(jù)成像方法及裝置
技術領域
[0001] 本申請涉及地球物理勘探技術領域，特別涉及一種垂直地震數(shù)據(jù)成像方法及裝 置。
【背景技術】
[0002] 垂直地震剖面(VSP)方法是一種介于地震和測井方法之間的地球物理技術，其是 一種地震觀測方法。它與通常地面觀測的地震剖面相對應。垂直地震剖面方法是在地表附 近的一些點上激發(fā)地震波，在沿井孔不同深度布置的一些多級多分量的檢波點上進行觀 測。在垂直地震剖面中，因為檢波器通過井置于地層內部，所以不僅能接收到自下而上傳播 的上行縱波和上行轉換波，也能接收到自上而下傳播的下行縱波及下行轉換波，甚至能接 收到橫波。因此該地震方法比常規(guī)地震方法具有更高的信噪比和分辨率。然而VSP成像技術 卻一直制約著該方法的發(fā)展，其中兩個主要的原因分別是:該觀測方式中不能利用均勻的 多次覆蓋改善成像效果，以及轉換成像或偏移成像過程中產(chǎn)生了波場畸變。不能利用均勻 的多次覆蓋改善成像效果是由炮檢位置關系決定，難以通過室內處理改善提高。轉換成像 或偏移成像過程中產(chǎn)生了波場畸變是由于VSP成像處理過程中的動校正拉伸或偏移距拉伸 產(chǎn)生，可以通過一些保真處理手段進行改善。
[0003] 為了改善波場畸變，現(xiàn)有的垂直地震數(shù)據(jù)成像方法通常包括:對垂直地震數(shù)據(jù)進 行預處理操作，所述預處理操作可以包括：自動增益處理、道均衡處理或帶通波處理，再將 所述預處理后的垂直地震數(shù)據(jù)直接向地震共深度點道集(⑶P)進行行數(shù)據(jù)轉換，即VSP-CDP (垂直地震剖面-地震共深度點道集)成像。
[0004] 發(fā)明人發(fā)現(xiàn)現(xiàn)有技術中至少存在如下問題：現(xiàn)有技術中的垂直地震數(shù)據(jù)成像方 法，在利用VSP-CDP成像的過程中，先對垂直地震數(shù)據(jù)進行預處理操作，再將包括振幅、相 位、頻率的地震數(shù)據(jù)采用VSP-CDP技術進行成像操作，該成像過程可以保真振幅，但頻率和 相位仍然有較大的畸變，從而導致頻率和相位的保真性較低。

【發(fā)明內容】

[0005] 本申請實施例的目的是提供一種垂直地震數(shù)據(jù)成像方法及裝置，以提高垂直地震 數(shù)據(jù)的頻率和相位的保真性。
[0006] 為解決上述技術問題，本申請實施例提供一種垂直地震數(shù)據(jù)成像方法及裝置是這 樣實現(xiàn)的：
[0007] -種垂直地震數(shù)據(jù)成像方法，包括:獲取采集到的原始垂直地震數(shù)據(jù)，對所述采集 到的原始垂直地震數(shù)據(jù)進行預處理;采用預設的參數(shù)提取方法提取所述預處理后的垂直地 震數(shù)據(jù)的參數(shù)對應的數(shù)據(jù);所述參數(shù)包括:振幅、頻率和相位;采用第一成像方法將所述提 取的每一種參數(shù)對應的數(shù)據(jù)分別轉換為參數(shù)數(shù)據(jù)體;根據(jù)所述參數(shù)數(shù)據(jù)體，利用第二成像 方法得到垂直地震數(shù)據(jù)成像。
[0008] 優(yōu)選方案中，所述預設的參數(shù)提取方法，具體包括:確定滑動時間窗長度和所述滑 動時間窗的位置;在所述滑動時間窗長度內對所述垂直地震數(shù)據(jù)進行正弦擬合運算得到擬 合函數(shù);提取所述擬合函數(shù)中的所述參數(shù)對應的數(shù)據(jù)，將所述提取的數(shù)據(jù)作為所述滑動時 間窗中心位置處對應的參數(shù)數(shù)據(jù)。
[0009] 優(yōu)選方案中，所述滑動時間窗長度大于或等于一個波長對應的時間長度。
[0010] 優(yōu)選方案中，所述滑動時間窗的長度為35毫秒~50毫秒。
[0011] 優(yōu)選方案中，所述預設的參數(shù)提取方法還包括:從所述預處理后的垂直地震數(shù)據(jù) 的起始位置開始，每提取一次參數(shù)，滑動所述時間窗單位時長，直至所述預處理后的垂直地 震數(shù)據(jù)的終點位置。
[0012] 優(yōu)選方案中，所述擬合函數(shù)的表達式為：
[0013] x = Am\(2n: j) + q>)
[0014] 其中，x表示當前滑動時窗中擬合的函數(shù);A表示振幅，f表示頻率，t表示時間，9表 示相位;其中，A、f、<P為所述提取的參數(shù)對應的數(shù)據(jù)。
[0015]優(yōu)選方案中，所述采用第一成像方法具體包括:采用VSP-⑶P的成像方法，即垂直 地震剖面-地震共深度點道集的成像方法。
[0016] 優(yōu)選方案中，根據(jù)所述參數(shù)數(shù)據(jù)體，利用第二成像方法得到垂直地震數(shù)據(jù)成像，包 括:將相同位置處的三個參數(shù)數(shù)據(jù)體作為一組，帶入正弦函數(shù)中計算得到新的地震數(shù)據(jù)，所 有的地震數(shù)據(jù)的集合為垂直地震數(shù)據(jù)成像。
[0017] 優(yōu)選方案中，所述對采集到的原始垂直地震數(shù)據(jù)進行預處理，包括下述中的至少 一種:帶通濾波處理、相位轉換處理、旋轉定向處理、反褶積計算、波場分離處理。
[0018] 一種垂直地震數(shù)據(jù)成像裝置，包括:預處理模塊、參數(shù)提取模塊、參數(shù)分別轉換模 塊和成像模塊;其中，所述預處理模塊，用于獲取采集到的原始垂直地震數(shù)據(jù)，對所述采集 到的原始垂直地震數(shù)據(jù)進行預處理;所述參數(shù)提取模塊，用于采用預設的參數(shù)提取方法提 取所述預處理模塊預處理后的垂直地震數(shù)據(jù)的參數(shù)對應的數(shù)據(jù);所述參數(shù)包括:振幅、頻率 和相位;所述參數(shù)分別轉換模塊，用于采用第一成像方法將所述參數(shù)提取模塊提取的每一 種參數(shù)對應的數(shù)據(jù)分別轉換為參數(shù)數(shù)據(jù)體;所述成像模塊，用于根據(jù)所述參數(shù)分別轉換模 塊轉換的參數(shù)數(shù)據(jù)體，利用第二成像方法得到垂直地震數(shù)據(jù)成像。
[0019] 優(yōu)選方案中，所述參數(shù)提取模塊包括:確定模塊、擬合模塊和提取模塊;其中，所述 確定模塊，用于確定滑動時間窗長度和所述滑動時間窗的位置;所述擬合模塊，用于在所述 確定模塊，確定的滑動時間窗長度內對所述垂直地震數(shù)據(jù)進行正弦擬合運算得到擬合函 數(shù);所述提取模塊，用于提取所述擬合模塊，擬合的函數(shù)中的參數(shù)對應的數(shù)據(jù)，將所述提取 的數(shù)據(jù)作為所述滑動時間窗中心位置處對應的參數(shù)數(shù)據(jù);所述參數(shù)包括:振幅、頻率和相 位。
[0020] 優(yōu)選方案中，所述參數(shù)提取模塊還包括:滑動模塊，其用于從所述預處理后的垂直 地震數(shù)據(jù)的起始位置開始，每提取一次參數(shù)，滑動所述時間窗單位時長，直至所述預處理后 的垂直地震數(shù)據(jù)的終點位置。
[0021] 由以上本申請實施例提供的技術方案可見，本申請實施例提供的垂直地震數(shù)據(jù)成 像方法和裝置，先提取垂直地震數(shù)據(jù)中的振幅、頻率和相位三種參數(shù)，再基于VSP-CDP成像 方法將三種參數(shù)的數(shù)據(jù)分別轉換為數(shù)據(jù)體，再將同一位置的三個數(shù)據(jù)體進行合成，避免了 VSP-CDP成像過程中數(shù)據(jù)拉伸和畸變影響，可以使成像數(shù)據(jù)的振幅、頻率和相位得到最大限 度的保護，從而提高垂直地震數(shù)據(jù)的頻率和相位的保真性。同時，本申請對原始垂直地震數(shù) 據(jù)進行預處理的過程中，沒有采用破壞振幅、頻率、相位相對關系的處理操作，可以進一步 提高垂直地震數(shù)據(jù)的頻率和相位的保真性。
【附圖說明】
[0022] 為了更清楚地說明本申請實施例或現(xiàn)有技術中的技術方案，下面將對實施例或現(xiàn) 有技術描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本 申請中記載的一些實施例，對于本領域普通技術人員來講，在不付出創(chuàng)造性勞動性的前提 下，還可以根據(jù)這些附圖獲得其他的附圖。
[0023] 圖1是本申請一個垂直地震數(shù)據(jù)成像方法實施例的流程圖；
[0024] 圖2是本申請一個垂直地震數(shù)據(jù)成像裝置實施例的模塊圖；
[0025] 圖3是本申請裝置實施例中參數(shù)提取模塊的一個組成圖。
【具體實施方式】
[0026] 本申請實施例提供一種垂直地震數(shù)據(jù)成像方法及裝置。
[0027] 為了使本技術領域的人員更好地理解本申請中的技術方案，下面將結合本申請實 施例中的附圖，對本申請實施例中的技術方案進行清楚、完整地描述，顯然，所描述的實施 例僅僅是本申請一部分實施例，而不是全部的實施例?；诒旧暾堉械膶嵤├?，本領域普通 技術人員在沒有作出創(chuàng)造性勞動前提下所獲得的所有其他實施例，都應當屬于本申請保護 的范圍。
[0028] 圖1是本申請一個垂直地震數(shù)據(jù)成像方法實施例的流程圖。如圖1所示，所述垂直 地震數(shù)據(jù)成像方法可以包括：
[0029] S101:獲取采集到的原始垂直地震數(shù)據(jù)，對所述采集到的原始垂直地震數(shù)據(jù)進行 預處理。
[0030] 可以獲取采集到的原始垂直地震數(shù)據(jù)。所述原始垂直地震數(shù)據(jù)可以通過垂直地震 剖面(VSP)方法獲取到。
[0031 ]所述原始垂直地震數(shù)據(jù)中可以包括振幅、相位和頻率。
[0032]可以對所述采集到的原始垂直地震數(shù)據(jù)進行預處理。
[0033]所述預處理可以包括下述操作中的至少一種：帶通濾波處理、相位轉換處理、旋轉 定向處理、反褶積計算、波場分離處理。
[0034] 通過上述預處理操作，可以實現(xiàn)對原始垂直地震數(shù)據(jù)中的去噪處理，得到轉換成 像前的垂直地震數(shù)據(jù)。上述預處理的步驟一般不會破壞頻率、振幅和相位的相對關系，也不 會改變數(shù)據(jù)的連續(xù)性，因此可以較好地實現(xiàn)對原始垂直地震數(shù)據(jù)的保真。
[0035] S102:采用預設的參數(shù)提取方法提取所述預處理后的垂直地震數(shù)據(jù)的參數(shù)對應的 數(shù)據(jù)。
[0036] 可以采用預設的參數(shù)提取方法提取所述預處理后的垂直地震數(shù)據(jù)的參數(shù)對應的 數(shù)據(jù)，其中，所述參數(shù)可以包括:振幅、頻率和相位。那么，提取的數(shù)據(jù)可以分別是振幅數(shù)據(jù)、 頻率數(shù)據(jù)和相位數(shù)據(jù)。
[0037] 所述預設的參數(shù)提取方法，具體可以包括:確定滑動時間窗長度和所述滑動時間 窗的位置，在所述滑動時間窗長度內對所述垂直地震數(shù)據(jù)進行正弦擬合運算得到擬合函 數(shù)，提取所述擬合函數(shù)中的所述參數(shù)對應的數(shù)據(jù);將所述提取的數(shù)據(jù)作為所述滑動時間窗 中心位置處對應的參數(shù)數(shù)據(jù)。
[0038]其中，所述滑動時間窗的長度可以根據(jù)所述垂直地震數(shù)據(jù)的頻率來確定，具體地， 所述滑動時間窗的長度可以大于或等于一個波長對應的時間長度。
[0039]例如，垂直地震數(shù)據(jù)的頻率約為20-30HZ(赫茲），其對應的波長為33-50毫秒。那么 所述滑動時間窗的長度可以為33-50毫秒，一般可以選取為50毫秒。
[0040]進一步地，所述擬合函數(shù)的表達式為：
[0041 ] .二.4 si n (2;r./i十 p) (1)
[0042] 其中，x表示當前滑動時窗中擬合的函數(shù);A表示振幅，f表示頻率，t表示時間，$表 示相位;其中，A、f、f為所述提取的參數(shù)對應的數(shù)據(jù)。
[0043]需要說明的是，所述對滑動時間窗長度內所述垂直地震數(shù)據(jù)進行正弦擬合運算可 以通過現(xiàn)有的數(shù)學計算方法實現(xiàn)，本申請對此并不作出限定
[0044] 進一步地，所述預設的參數(shù)提取方法還可以包括:從所述預處理后的垂直地震數(shù) 據(jù)的起始位置開始，每提取一次參數(shù)，滑動所述時間窗單位時長，直至所述預處理后的垂直 地震數(shù)據(jù)的終點位置。所述單位時長可以是1毫秒。
[0045] 通過逐步滑動所述滑動時間窗，可以確定每個時間窗中心位置處的參數(shù)數(shù)據(jù)。
[0046] 所述確定的每個參數(shù)對應的數(shù)據(jù)的個數(shù)一般小于所述預處理后的垂直地震數(shù)據(jù) 的參數(shù)對應的數(shù)據(jù)個數(shù)。由于滑動時間窗具備一定的時間窗長度，時間窗的中心位置不能 到達所述預處理后的垂直地震數(shù)據(jù)的起始位置和終點位置，所以不能到達的位置處參數(shù)數(shù) 據(jù)被舍棄。舍棄的參數(shù)數(shù)據(jù)相比于所有地震數(shù)據(jù)數(shù)量較小，因此，不會影響地震數(shù)據(jù)的真實 性。
[0047] S103:采用第一成像方法將所述提取的每一種參數(shù)對應的數(shù)據(jù)分別轉換為參數(shù)數(shù) 據(jù)體。
[0048] 可以采用第一成像方法，將所述提取的每一種參數(shù)對應的數(shù)據(jù)分別轉換為參數(shù)數(shù) 據(jù)體。
[0049] 具體地，所述第一成像方法可以采用VSP-CDP的成像方法，可以實現(xiàn)振幅的保真 性。
[0050] 分別轉換成像的過程中，所述第一成像方法還可以采用現(xiàn)有技術中的其他方法， 本申請對此并不作出限定。
[0051] S104:根據(jù)所述參數(shù)數(shù)據(jù)體，利用第二成像方法得到垂直地震數(shù)據(jù)成像。
[0052]可以根據(jù)步驟S103中所述轉換的參數(shù)數(shù)據(jù)體，利用第二成像方法得到垂直地震數(shù) 據(jù)成像。所述第二成像方法具體可以包括:將相同位置處的三個參數(shù)數(shù)據(jù)體作為一組，帶入 正弦函數(shù)，即公式（1)，中計算得到新的地震數(shù)據(jù)，所有的地震數(shù)據(jù)的集合為垂直地震數(shù)據(jù) 成像。
[0053]上述實施例提供的垂直地震數(shù)據(jù)成像方法，先提取垂直地震數(shù)據(jù)中的振幅、頻率 和相位三種參數(shù)，再基于VSP-CDP成像方法將三種參數(shù)的數(shù)據(jù)分別轉換為數(shù)據(jù)體，再將同一 位置的三個數(shù)據(jù)體進行合成，避免了 VSP-CDP成像過程中數(shù)據(jù)拉伸和畸變影響，可以使成像 數(shù)據(jù)的振幅、頻率和相位得到最大限度的保護，從而提高垂直地震數(shù)據(jù)的頻率和相位的保 真性。同時，本申請對原始垂直地震數(shù)據(jù)進行預處理的過程中，沒有采用破壞振幅、頻率、相 位相對關系的處理操作，可以進一步提高垂直地震數(shù)據(jù)的頻率和相位的保真性。
[0054] 下面介紹本申請的垂直地震數(shù)據(jù)成像裝置實施例。
[0055] 圖2是本申請一個垂直地震數(shù)據(jù)成像裝置實施例的模塊圖。如圖2所示，所述垂直 地震數(shù)據(jù)成像裝置可以包括:預處理模塊201、參數(shù)提取模塊202、參數(shù)分別轉換模塊203和 成像模塊204。其中，
[0056]所述預處理模塊201，可以用于獲取采集到的原始垂直地震數(shù)據(jù)，對所述采集到的 原始垂直地震數(shù)據(jù)進行預處理。
[0057]所述參數(shù)提取模塊202,可以用于采用預設的參數(shù)提取方法提取所述預處理模塊 201預處理后的垂直地震數(shù)據(jù)的參數(shù)對應的數(shù)據(jù);所述參數(shù)包括:振幅、頻率和相位。
[0058]所述參數(shù)分別轉換模塊203,可以用于采用第一成像方法將所述參數(shù)提取模塊202 提取的每一種參數(shù)對應的數(shù)據(jù)分別轉換為參數(shù)數(shù)據(jù)體。
[0059]所述成像模塊204,可以用于根據(jù)所述參數(shù)分別轉換模塊203轉換的參數(shù)數(shù)據(jù)體， 利用第二成像方法得到垂直地震數(shù)據(jù)成像。
[0060] 圖3是本申請裝置實施例中參數(shù)提取模塊的一個組成圖。如圖3所示，所述參數(shù)提 取模塊202可以包括:確定模塊2021、擬合模塊2022和提取模塊2023。
[0061] 所述確定模塊2021，可以用于確定滑動時間窗長度和所述滑動時間窗的位置。所 述滑動時間窗的長度可以根據(jù)所述垂直地震數(shù)據(jù)的頻率來確定，具體地，所述滑動時間窗 的長度可以大于或等于一個波長對應的時間長度。例如，垂直地震數(shù)據(jù)的頻率約為20-30HZ (赫茲），其對應的波長為33-50暈秒。那么所述滑動時間窗的長度可以為33-50暈秒，一般可 以選取為50毫秒。
[0062]所述擬合模塊2022,可以用于在所述確定模塊2021確定的滑動時間窗長度內對所 述垂直地震數(shù)據(jù)進行正弦擬合運算得到擬合函數(shù)。
[0063]所述提取模塊2023,可以用于提取所述擬合模塊2022擬合的函數(shù)中的參數(shù)對應的 數(shù)據(jù)，將所述提取的數(shù)據(jù)作為所述滑動時間窗中心位置處對應的參數(shù)數(shù)據(jù)。所述參數(shù)可以 包括:振幅、頻率和相位。
[0064]進一步地，所述參數(shù)提取模塊202還可以包括:滑動模塊2024，其可以用于從所述 預處理后的垂直地震數(shù)據(jù)的起始位置開始，每提取一次參數(shù)，滑動所述時間窗單位時長，直 至所述預處理后的垂直地震數(shù)據(jù)的終點位置。
[0065]上述實施例公開的垂直地震數(shù)據(jù)成像裝置與本申請的垂直地震數(shù)據(jù)成像方法實 施例相對應，可以實現(xiàn)本申請的方法實施例并取得方法實施例的技術效果。
[0066]在20世紀90年代，對于一個技術的改進可以很明顯地區(qū)分是硬件上的改進(例如， 對二極管、晶體管、開關等電路結構的改進)還是軟件上的改進(對于方法流程的改進）。然 而，隨著技術的發(fā)展，當今的很多方法流程的改進已經(jīng)可以視為硬件電路結構的直接改進。 設計人員幾乎都通過將改進的方法流程編程到硬件電路中來得到相應的硬件電路結構。因 此，不能說一個方法流程的改進就不能用硬件實體模塊來實現(xiàn)。例如，可編程邏輯器件 (Programmable Logic Device,PLD)(例如現(xiàn)場可編程門陣列（Field Programmable Gate Array，F(xiàn)PGA))就是這樣一種集成電路，其邏輯功能由用戶對器件編程來確定。由設計人員 自行編程來把一個數(shù)字系統(tǒng)"集成"在一片PLD上，而不需要請芯片制造廠商來設計和制作 專用的集成電路芯片2。而且，如今，取代手工地制作集成電路芯片，這種編程也多半改用 "邏輯編譯器(logic compiler)"軟件來實現(xiàn)，它與程序開發(fā)撰寫時所用的軟件編譯器相類 似，而要編譯之前的原始代碼也得用特定的編程語言來撰寫，此稱之為硬件描述語言 (Hardware Description Language，HDL)，而HDL也并非僅有一種，而是有許多種，如ABEL (Advanced Boolean Expression Language)、AHDL(Altera Hardware Description Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL (Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby Hardware Description Language)等，目前最普遍使用的是VHDL(Very-High_Speed Integrated Circuit Hardware Description Language)與Verilog2〇本令頁域技術人員也 應該清楚，只需要將方法流程用上述幾種硬件描述語言稍作邏輯編程并編程到集成電路 中，就可以很容易得到實現(xiàn)該邏輯方法流程的硬件電路。
[0067] 控制器可以按任何適當?shù)姆绞綄崿F(xiàn)，例如，控制器可以采取例如微處理器或處理 器以及存儲可由該(微)處理器執(zhí)行的計算機可讀程序代碼(例如軟件或固件）的計算機可 讀介質、邏輯門、開關、專用集成電路(Application Specific Integrated Circuit， ASIC)、可編程邏輯控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制 器：ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320,存 儲器控制器還可以被實現(xiàn)為存儲器的控制邏輯的一部分。
[0068] 本領域技術人員也知道，除了以純計算機可讀程序代碼方式實現(xiàn)控制器以外，完 全可以通過將方法步驟進行邏輯編程來使得控制器以邏輯門、開關、專用集成電路、可編程 邏輯控制器和嵌入微控制器等的形式來實現(xiàn)相同功能。因此這種控制器可以被認為是一種 硬件部件，而對其內包括的用于實現(xiàn)各種功能的裝置也可以視為硬件部件內的結構。或者 甚至，可以將用于實現(xiàn)各種功能的裝置視為既可以是實現(xiàn)方法的軟件模塊又可以是硬件部 件內的結構。
[0069] 上述實施例闡明的系統(tǒng)、裝置、模塊或單元，具體可以由計算機芯片或實體實現(xiàn)， 或者由具有某種功能的產(chǎn)品來實現(xiàn)。
[0070] 為了描述的方便，描述以上裝置時以功能分為各種單元分別描述。當然，在實施本 申請時可以把各單元的功能在同一個或多個軟件和/或硬件中實現(xiàn)。
[0071]通過以上的實施方式的描述可知，本領域的技術人員可以清楚地了解到本申請可 借助軟件加必需的通用硬件平臺的方式來實現(xiàn)。基于這樣的理解，本申請的技術方案本質 上或者說對現(xiàn)有技術做出貢獻的部分可以以軟件產(chǎn)品的形式體現(xiàn)出來，在一個典型的配置 中，計算設備包括一個或多個處理器(CPU)、輸入/輸出接口、網(wǎng)絡接口和內存。該計算機軟 件產(chǎn)品可以包括若干指令用以使得一臺計算機設備(可以是個人計算機，服務器，或者網(wǎng)絡 設備等)執(zhí)行本申請各個實施例或者實施例的某些部分所述的方法。該計算機軟件產(chǎn)品可 以存儲在內存中，內存可能包括計算機可讀介質中的非永久性存儲器，隨機存取存儲器 (RAM)和/或非易失性內存等形式，如只讀存儲器(ROM)或閃存(flash RAM)。內存是計算機 可讀介質的示例。計算機可讀介質包括永久性和非永久性、可移動和非可移動媒體可以由 任何方法或技術來實現(xiàn)信息存儲。信息可以是計算機可讀指令、數(shù)據(jù)結構、程序的模塊或其 他數(shù)據(jù)。計算機的存儲介質的例子包括，但不限于相變內存(PRAM)、靜態(tài)隨機存取存儲器 (SRAM)、動態(tài)隨機存取存儲器（DRAM)、其他類型的隨機存取存儲器（RAM)、只讀存儲器 (ROM)、電可擦除可編程只讀存儲器(EEPROM)、快閃記憶體或其他內存技術、只讀光盤只讀 存儲器(CD-ROM)、數(shù)字多功能光盤(DVD)或其他光學存儲、磁盒式磁帶，磁帶磁磁盤存儲或 其他磁性存儲設備或任何其他非傳輸介質，可用于存儲可以被計算設備訪問的信息。按照 本文中的界定，計算機可讀介質不包括短暫電腦可讀媒體(transitory media)，如調制的 數(shù)據(jù)信號和載波。
[0072] 本說明書中的各個實施例均采用遞進的方式描述，各個實施例之間相同相似的部 分互相參見即可，每個實施例重點說明的都是與其他實施例的不同之處。尤其，對于系統(tǒng)實 施例而言，由于其基本相似于方法實施例，所以描述的比較簡單，相關之處參見方法實施例 的部分說明即可。
[0073] 本申請可用于眾多通用或專用的計算機系統(tǒng)環(huán)境或配置中。例如:個人計算機、月艮 務器計算機、手持設備或便攜式設備、平板型設備、多處理器系統(tǒng)、基于微處理器的系統(tǒng)、置 頂盒、可編程的消費電子設備、網(wǎng)絡PC、小型計算機、大型計算機、包括以上任何系統(tǒng)或設備 的分布式計算環(huán)境等等。
[0074] 本申請可以在由計算機執(zhí)行的計算機可執(zhí)行指令的一般上下文中描述，例如程序 模塊。一般地，程序模塊包括執(zhí)行特定任務或實現(xiàn)特定抽象數(shù)據(jù)類型的例程、程序、對象、組 件、數(shù)據(jù)結構等等。也可以在分布式計算環(huán)境中實踐本申請，在這些分布式計算環(huán)境中，由 通過通信網(wǎng)絡而被連接的遠程處理設備來執(zhí)行任務。在分布式計算環(huán)境中，程序模塊可以 位于包括存儲設備在內的本地和遠程計算機存儲介質中。
[0075] 雖然通過實施例描繪了本申請，本領域普通技術人員知道，本申請有許多變形和 變化而不脫離本申請的精神，希望所附的權利要求包括這些變形和變化而不脫離本申請的 精神。
【主權項】
1. 一種垂直地震數(shù)據(jù)成像方法，其特征在于，包括： 獲取采集到的原始垂直地震數(shù)據(jù)，對所述采集到的原始垂直地震數(shù)據(jù)進行預處理； 采用預設的參數(shù)提取方法提取所述預處理后的垂直地震數(shù)據(jù)的參數(shù)對應的數(shù)據(jù);所述 參數(shù)包括:振幅、頻率和相位； 采用第一成像方法將所述提取的每一種參數(shù)對應的數(shù)據(jù)分別轉換為參數(shù)數(shù)據(jù)體； 根據(jù)所述參數(shù)數(shù)據(jù)體，利用第二成像方法得到垂直地震數(shù)據(jù)成像。2. 如權利要求1所述的一種垂直地震數(shù)據(jù)成像方法，其特征在于，所述預設的參數(shù)提取 方法，具體包括： 確定滑動時間窗長度和所述滑動時間窗的位置； 在所述滑動時間窗長度內對所述垂直地震數(shù)據(jù)進行正弦擬合運算得到擬合函數(shù)； 提取所述擬合函數(shù)中的所述參數(shù)對應的數(shù)據(jù)，將所述提取的數(shù)據(jù)作為所述滑動時間窗 中屯、位置處對應的參數(shù)數(shù)據(jù)。3. 如權利要求2所述的一種垂直地震數(shù)據(jù)成像方法，其特征在于，所述滑動時間窗長度 大于或等于一個波長對應的時間長度。4. 如權利要求3所述的一種垂直地震數(shù)據(jù)成像方法，其特征在于，所述滑動時間窗的長 度為35毫秒~50毫秒。5. 如權利要求2所述的一種垂直地震數(shù)據(jù)成像方法，其特征在于，所述預設的參數(shù)提取 方法還包括:從所述預處理后的垂直地震數(shù)據(jù)的起始位置開始，每提取一次參數(shù)，滑動所述 時間窗單位時長，直至所述預處理后的垂直地震數(shù)據(jù)的終點位置。6. 如權利要求2所述的一種垂直地震數(shù)據(jù)成像方法，其特征在于，所述擬合函數(shù)的表達 式為：其中，X表示當前滑動時窗中擬合的函數(shù);A表示振幅，f表示頻率，t表示時間，表示相 位;其中，A、f、q>為所述提取的參數(shù)對應的數(shù)據(jù)。7. 如權利要求1所述的一種垂直地震數(shù)據(jù)成像方法，其特征在于，所述采用第一成像方 法具體包括:采用VSP-CDP的成像方法，即垂直地震剖面-地震共深度點道集的成像方法。8. 如權利要求1所述的一種垂直地震數(shù)據(jù)成像方法，其特征在于，根據(jù)所述參數(shù)數(shù)據(jù) 體，利用第二成像方法得到垂直地震數(shù)據(jù)成像，包括:將相同位置處的=個參數(shù)數(shù)據(jù)體作為 一組，帶入正弦函數(shù)中計算得到新的地震數(shù)據(jù)，所有的地震數(shù)據(jù)的集合為垂直地震數(shù)據(jù)成 像。9. 如權利要求1所述的一種垂直地震數(shù)據(jù)成像方法，其特征在于，所述對采集到的原始 垂直地震數(shù)據(jù)進行預處理，包括下述中的至少一種：帶通濾波處理、相位轉換處理、旋轉定 向處理、反權積計算、波場分離處理。10. -種垂直地震數(shù)據(jù)成像裝置，其特征在于，包括:預處理模塊、參數(shù)提取模塊、參數(shù) 分別轉換模塊和成像模塊;其中， 所述預處理模塊，用于獲取采集到的原始垂直地震數(shù)據(jù)，對所述采集到的原始垂直地 震數(shù)據(jù)進行預處理； 所述參數(shù)提取模塊，用于采用預設的參數(shù)提取方法提取所述預處理模塊預處理后的垂 直地震數(shù)據(jù)的參數(shù)對應的數(shù)據(jù);所述參數(shù)包括:振幅、頻率和相位； 所述參數(shù)分別轉換模塊，用于采用第一成像方法將所述參數(shù)提取模塊提取的每一種參 數(shù)對應的數(shù)據(jù)分別轉換為參數(shù)數(shù)據(jù)體； 所述成像模塊，用于根據(jù)所述參數(shù)分別轉換模塊轉換的參數(shù)數(shù)據(jù)體，利用第二成像方 法得到垂直地震數(shù)據(jù)成像。11. 如權利要求10所述的一種垂直地震數(shù)據(jù)成像裝置，其特征在于，所述參數(shù)提取模塊 包括:確定模塊、擬合模塊和提取模塊;其中， 所述確定模塊，用于確定滑動時間窗長度和所述滑動時間窗的位置； 所述擬合模塊，用于在所述確定模塊，確定的滑動時間窗長度內對所述垂直地震數(shù)據(jù) 進行正弦擬合運算得到擬合函數(shù)； 所述提取模塊，用于提取所述擬合模塊，擬合的函數(shù)中的參數(shù)對應的數(shù)據(jù)，將所述提取 的數(shù)據(jù)作為所述滑動時間窗中屯、位置處對應的參數(shù)數(shù)據(jù);所述參數(shù)包括:振幅、頻率和相 位。12. 如權利要求11所述的一種垂直地震數(shù)據(jù)成像裝置，其特征在于，所述參數(shù)提取模塊 還包括:滑動模塊，其用于從所述預處理后的垂直地震數(shù)據(jù)的起始位置開始，每提取一次參 數(shù)，滑動所述時間窗單位時長，直至所述預處理后的垂直地震數(shù)據(jù)的終點位置。
【文檔編號】G01V1/30GK105911586SQ201610228403
【公開日】2016年8月31日
【申請日】2016年4月13日
【發(fā)明人】蔡志東, 張慶紅, 趙建章, 王艷華, 王勇, 王沖
【申請人】中國石油天然氣集團公司, 中國石油集團東方地球物理勘探有限責任公司