[發明專利]基于激發主能量優化算法的全波形反演梯度算子提取方法

有效

• 申請號: 201410206451.0
• 申請日: 2014-05-16
• 公開（公告）號: CN105093278B
• 公開（公告）日: 2018-06-29
• 發明（設計）人: 慎國強;王玉梅;孟憲軍;王希萍;梁鍇;王延光;王興謀;苗永康;揭景榮;王振濤;劉立彬;張達;張洪;鈕學民;李燕;李海濤;李玉鳳;張睿璇
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勝利油田分公司物探研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 東營雙橋專利代理有限責任公司 37107
• 代理人: 侯華頌

本發明公開了一種基于激發主能量優化算法的全波形反演梯度算子的提取方法，通過建立全波形反演參數模型；采用常規的疊加前單炮數據作為觀測數據，建立正反演觀測系統；將震源置于地表炮點處，利用優化有限差分算法計算該震源產生的正傳波場；確定震源波場的正向傳播過程中地下各點震源波場最大值及其對應的時刻，并保存地下各點在激發主能量時窗內的波場；再利用優化有限差分模擬算法計算波場殘差逆時傳播波場，判斷地下各點不同時刻的波場殘差逆時傳播波場是否位于該點激發主能量時窗之內，通過基于激發主能量優化算法得到全波形反演梯度算子。本發明在保證精度的前提下，大大減少了I/O操作，提升了計算效率。

技術領域

本發明涉及石油地震勘探速度建模技術領域，尤其涉及一種基于激發主能量優化算法的全波形反演梯度算子的提取方法。

背景技術

地震波的傳播速度是地震勘探中的一個基本參數，在地震資料處理和地震解釋的諸多環節中扮演著重要的角色。影響地震波傳播速度的因素有很多，諸如地下介質的巖性、密度、孔隙度、埋藏深度等等，這些因素最終體現在地震波的走時、波形、振幅和相位等屬性上變化。只有充分了解速度及影響因素的相互關系，才有可能更為準確地估計出地下的結構和構造。我國大多數勘探區域越來越復雜，呈現出地表構造和地下構造雙復雜特征。傳統的速度分析方法已經很難達到所需的精度要求，也會直接影響成像質量。常用的速度分析方法有：疊加速度分析、偏移速度分析、層析速度反演和全波形反演。這四種速度分析方法的準確性依次升高，復雜度也依次升高。

全波形反演技術是目前精度很高的速度分析和反演工具，其利用疊前地震記錄與正演模擬炮記錄之間的尋求最佳匹配來獲取最優反演速度場。在反演過程中不僅考慮了地震波傳播的旅行時等運動學信息，而且加入了振幅、相位、波形等動力學信息，能夠適應強橫向變速介質和各向異性介質的速度反演。

全波形反演的發展大致可以分為三個階段，第一個階段是從20世紀80年代到90年代，這一時期是全波形反演理論形成和初步驗證階段，Lailly、Tarantola等人建立了基于波動方程的波形反演的理論構架，之后一批學者分別利用簡單模型對反射和透射波形反演效果進行了理論驗證。Gauthier等人率先對理論模型合成的多偏移距數據進行了聲波和彈性波全波形反演測試，反演結果既展示了全波形反演方法的應用潛力，也暴露出其根本缺陷——對初始模型的嚴重依賴性。

第二個階段是從上世紀90年代初到21世紀初，為全波形穩定發展階段，全波形反演方法形成了多個分支和變種，其中以Pratt等提出的頻率域波形反演方法最具代表性，Pratt提出頻率域聲波和彈性波波形反演理論，并利用實際井間地震數據對方法進行驗證。該方法利用稀疏矩陣直接解法求解頻率域的離散化聲波或彈性波波動方程，只要頻率與物性參數不改變，通過一次稀疏矩陣LU分解和多次高斯消元法中的前向、后向回代運算，即可完成所有炮點的單頻波場響應計算；應用研究方面以井間地震波形反演和海洋天然氣水合物的疊前波形反演最具影響力。此外為解決全波形反演因初始模型不準而陷入局部極小問題，Zhou等人先后給出了波動方程走時和波形聯合反演方法，以充分利用走時和波形兩類數據對不同尺度模型參數變化的敏感度的優勢，使井間數據全波形反演更為穩定。

自21世紀初至今，全波形反演達到了蓬勃發展的階段，全波形反演逐步從二維反演走向三維反演，從聲波走向彈性波反演，從單參數反演到多參數反演，從各向同性到各向異性反演。此外還有Laplace/Laplace-Fourier域全波形反演、各向異性全波形反演、先驗信息約束的全波形反演和基于超級炮和相位編碼的全波形反演。此外隨著計算機水平的進步，全波形反演的計算效率進一步得到了提升，其實現平臺逐步從單機單核走向了多核多節點并行平臺、GPU并行平臺。在一定程度上，全波形反演的計算效率完全取決于大數據的I/O交換。

發明內容