[發明專利]基于地面接收的反射地震資料的Q值場建模方法

有效

• 申請號: 201611001706.5
• 申請日: 2016-11-14
• 公開（公告）號: CN106443786B
• 公開（公告）日: 2018-04-20
• 發明（設計）人: 張劍鋒;劉偉;王進;劉禮農
• 申請（專利權）人: 中國科學院地質與地球物理研究所
• 主分類號: G01V1/30
• 分類號: G01V1/30;G01V1/34
• 代理公司: 工業和信息化部電子專利中心11010
• 代理人: 焉明濤

1.一種基于地面接收的反射地震資料的Q值場建模方法，其特征在于：包括以下步驟：步驟A、針對由拖纜或測線在地表記錄的人工震源激發的反射地震資料，利用疊前時間偏移方法得到時間深度域偏移速度場和一組等間距并平行于測線方向的時間深度域偏移成像剖面；步驟B、對記錄的反射地震資料進行高頻噪音道剔除，剔除含有強高頻噪音的地震道；步驟C、確定Q值的取值范圍，按1/Q等間距選取系列Q值q_i，用每個q_i，對剔除強高頻噪音道后的反射地震資料做反Q濾波計算，得到一系列Q值相關的Q補償地震資料；步驟D、針對每個Q值的Q補償地震資料，利用疊前時間偏移方法，在步驟A得到的一組偏移成像剖面處，求得一系列Q值相關的粘彈性偏移結果；其中，每個q_i的粘彈性偏移結果包括三個時間深度域的偏移成像剖面，一是所有偏移距的偏移結果疊加的成像剖面，二是偏移距小于中間偏移距的部分偏移結果疊加的成像剖面，三是偏移距大于中間偏移距的部分偏移結果疊加的成像剖面；步驟E、針對步驟A得到的每個時間深度域偏移成像剖面，選取一組確定Q值的Q值窗；步驟F、在每個Q值窗，針對每個q_i的粘彈性偏移結果，確定粘彈性偏移的高截止頻率，計算其頻譜包絡；步驟G、利用每個Q值窗的系列Q值的頻譜包絡，確定Q值窗中心處的等效Q值；步驟H、利用每個Q值窗的等效Q值和步驟A得到的時間深度域偏移速度場，求得時間深度域三維成像體上的等效Q值場；步驟I、沿時間深度方向對等效Q值進行反演，得到每個時間深度處的層Q值；步驟J、利用成像射線將時間深度域的層Q值轉換到深度域，得到深度域層Q值場。

2.根據權利 要求1所述的方法，其特征在于：步驟J之后還包括：

步驟K、對速度橫向變化屬于預設不劇烈的地下構造，利用等效Q值場和粘彈性疊前時間偏移方法，得到地下構造圖像；在速度橫向變化屬于預設劇烈的地下構造，應用深度域層Q值場和粘彈性疊前深度偏移方法，得到地下構造圖像；兩種地下構造圖像指示中、深層構造的斷裂體系和地層沉積樣式，用于確定地下生、儲油構造和識別油氣儲層。

3.根據權利 要求1所述的方法，其特征在于：所述的對記錄的反射地震資料進行高頻噪音道剔除，剔除含有強高頻噪音的地震道是這樣實現的：將記錄的反射地震資料按偏移距大小分組，對每組地震資料循環；先對該組地震資料進行付氏變換，設頻率域地震道為F_i(f)，i代表地震道的序號，f代表頻率，單位赫茲，設記錄的反射地震資料的主頻為f_p，計算比值

上式中f_max是粘彈性偏移計算采用的最高截止頻率，f_min是記錄的反射地震資料的最低頻率；對5b_i取整，將取整后數值相同的地震道歸為同一類型；選取所含地震道最多的一個類型，計算該類型的全部地震道的b_i的平均值，記為b₀；計算方差

式中n是該偏移距組中所含地震道的道數；對該組的每一地震道循環，若b_i-b₀>3σ，則判定該道含有強高頻噪音，予以剔除；完成對全部偏移距組的循環，即實現了剔除反射地震資料中含有強高頻噪音的地震道。

4.根據權利 要求1所述的方法，其特征在于：所述的針對步驟A得到的每個時間深度域偏移成像剖面，選取一組確定Q值的Q值窗是這樣實現的：在每個成像剖面上選取5到10個能覆蓋整個剖面的CDP點，并根據成像剖面上地質構造的形態確定CDP點的間距；在每個CDP位置，沿時間深度選取4到6個確定Q值的二維窗口；選定的窗口應滿足如下條件：在時間深度上至少要包含5個同相軸，在橫向上應包括40到60個CDP點，不要跨過斷層，窗口內的同相軸要有粗細變化，同相軸沿橫向要有幅值變化；在每個窗口，計算不同Q值對應的、所有偏移距疊加的時間深度域成像剖面在該窗口內的疊加頻譜，觀察頻譜的包絡是否存在高頻部分隨Q值倒數的增加而增加的趨勢，若不存在該趨勢，移動窗口的位置和調整窗口的大小，直至對應的疊加頻譜存在這一趨勢；對每個窗口觀察頻譜包絡并調整，最終在各個成像剖面上，確定一組確定Q值的Q值窗。