[發(fā)明專利]一種基于匹配追蹤的超低信噪比高精度速度譜生成方法有效
| 申請?zhí)枺?/td> | 202010068543.2 | 申請日： | 2020-01-21 |
| 公開（公告）號： | CN111239816B | 公開（公告）日： | 2021-06-25 |
| 發(fā)明（設(shè)計）人： | 謝玉洪;赫建偉;黎孝璋;裴健翔;鄧勇;任婷;張文祥;魯統(tǒng)祥;彭海龍;王瑞敏 | 申請（專利權(quán)）人： | 中國海洋石油集團有限公司;中海油海南能源有限公司;中海石油（中國）有限公司湛江分公司;中海石油（中國）有限公司海南分公司 |
| 主分類號： | G01V1/30 | 分類號： | G01V1/30 |
| 代理公司： | 廣州市南鋒專利事務(wù)所有限公司 44228 | 代理人： | 李慧 |
| 地址： | 100010 北*** | 國省代碼： | 北京;11 |
| 權(quán)利要求書： | 查看更多 | 說明書： | 查看更多 |
| 摘要： | |||
| 搜索關(guān)鍵詞： | 一種 基于 匹配 追蹤 超低信噪 比高 精度 速度 生成 方法 | ||
本發(fā)明公開了一種基于匹配追蹤的超低信噪比高精度速度譜生成方法，該方法包括以下步驟：步驟(1)向計算機輸入計算參數(shù)，曲波變換的最大頻率fmax，有效信號曲波域分布范圍Ω，濾波閾值百分比ε，增益參數(shù)s，收斂系數(shù)α；步驟(2)輸入某一待處理CMP地震道集d；步驟(3)利用(1)中輸入的曲波變換的參數(shù)生成完備的原子庫A；步驟(4)令k＝1，設(shè)置初始迭代結(jié)果設(shè)置殘差rk＝d，迭代分解有效信號，保留Ω內(nèi)系數(shù)百分比ε范圍的曲波系數(shù)并重建得到對有效信號進(jìn)行增強及殘差以計算高精度速度譜，迭代終止條件為如果滿足，輸出步驟(5)設(shè)定速度計算范圍v∈[v1,v2]，速度間隔Δv，利用計算速度譜常規(guī)雙曲線振幅疊加及高精度加權(quán)系數(shù)計算高精度速度譜Sembhi＝Semb·S，sa為指數(shù)因子，輸出速度譜Sembhi。
技術(shù)領(lǐng)域
本發(fā)明涉及地震勘探技術(shù)領(lǐng)域，具體為一種基于匹配追蹤的超低信噪比高精度速度譜生成方法。
背景技術(shù)
地震速度分析是地震數(shù)據(jù)處理中非常重要的一個流程，地震速度的精度直接影響地震成像的好壞，準(zhǔn)確的地震速度可以使復(fù)雜構(gòu)造、低信噪比等難以成像的地區(qū)得以較好地成像，反之，如果成像速度與真實速度偏差較大，這些地質(zhì)情況下，地震數(shù)據(jù)則很難成像，從而嚴(yán)重制約地下地質(zhì)構(gòu)造認(rèn)識。目前，所有的地震速度分析均是從時間域速度譜拾取開始的，時間域速度通常作為后續(xù)深度域速度反演的初始速度，如果初始速度與真實速度偏離較大，后續(xù)速度反演則會出現(xiàn)不收斂或者往錯誤的解收斂等問題。常規(guī)的速度分析方法普遍采用相干性度量方法進(jìn)行速度分析，也有引入了相位信息的更高精度速度分析方法，但是這些方法針對信噪比極低的數(shù)據(jù)，例如多次波能量非常強、完全掩蓋有效信號的沉積破碎帶地區(qū)的海洋地震資料，均無法凸顯有效信號的速度特征。海洋地震資料最主要的干擾波就是多次波，多次波振幅相對于中深層有效反射能量相對較強，在沒有壓制多次波或者多次波無法壓制的情況下，道集速度譜往往凸顯的是多次波的能量，一次波的能量較弱甚至無法識別。針對以上所述，那么如何發(fā)明出一種針對超低信噪比、能識別有效弱反射的高精度速度譜生成方法，這成為我們需要解決的問題。
發(fā)明內(nèi)容
本發(fā)明的目的在于提供一種基于匹配追蹤的超低信噪比高精度速度譜生成方法，解決了背景技術(shù)中所提出的問題。
為解決上述問題，本發(fā)明提供如下技術(shù)方案：
一種基于匹配追蹤的超低信噪比高精度速度譜生成方法，其總的計算流程如圖1所示，該方法包括以下步驟：
步驟(1)向計算機輸入計算參數(shù)，曲波變換的最大頻率fmax，有效信號曲波域分布范圍Ω，濾波閾值百分比ε，增益參數(shù)s，收斂系數(shù)α；
步驟(2)輸入某一待處理CMP地震道集d；
步驟(3)利用(1)中輸入的曲波變換的參數(shù)生成完備的原子庫A，原子庫尺度層數(shù)為log2(Nt)-2，Nt為曲波變換時間樣點數(shù)且必須為2的指數(shù)，如果Nt不等于2的指數(shù)則取值2的最相近指數(shù)；最內(nèi)層與最外層分解角度數(shù)為1，其它層分解角度數(shù)依次為floor[[3:1:log2(Nt)-1]/2]+1，其中floor表示向下取整，a:b:c表示a、c之間步長為b的所有整數(shù)；
步驟(4)令k＝1，設(shè)置初始迭代結(jié)果設(shè)置殘差rk＝d，迭代計算：
步驟(4.1)利用A對rk進(jìn)行曲波域分解得到Ck；
步驟(4.2)保留Ω內(nèi)系數(shù)百分比ε范圍曲波系數(shù)參數(shù)ε為當(dāng)前迭代Ck在Ω內(nèi)最大能量值的百分比，控制稀疏性提高迭代的穩(wěn)定性和抗噪性，參數(shù)ε保證了每次迭代只保留Ω內(nèi)能量相對較大的有效信號；
步驟(4.3)利用對數(shù)據(jù)進(jìn)行重建得到
該專利技術(shù)資料僅供研究查看技術(shù)是否侵權(quán)等信息，商用須獲得專利權(quán)人授權(quán)。該專利全部權(quán)利屬于中國海洋石油集團有限公司;中海油海南能源有限公司;中海石油（中國）有限公司湛江分公司;中海石油（中國）有限公司海南分公司，未經(jīng)中國海洋石油集團有限公司;中海油海南能源有限公司;中海石油（中國）有限公司湛江分公司;中海石油（中國）有限公司海南分公司許可，擅自商用是侵權(quán)行為。如果您想購買此專利、獲得商業(yè)授權(quán)和技術(shù)合作，請聯(lián)系【客服】
本文鏈接：http://www.szxzyx.cn/pat/books/202010068543.2/2.html，轉(zhuǎn)載請聲明來源鉆瓜專利網(wǎng)。
