簡介本資源是一套面向地球物理專業(yè)研究生、油氣勘探工程師及地震數(shù)據(jù)處理從業(yè)者的疊前三參數(shù)反演MATLAB實現(xiàn)工具包聚焦于AVO反演核心問題——從疊前道集同步反演縱波速度、橫波速度與密度三個關(guān)鍵地質(zhì)參數(shù)顯著提升儲層流體識別與巖性解釋精度。壓縮包共12個文件10個.m主程序腳本2個.txt配置文件總大小僅9KB輕量緊湊其中包含正演建模syntrace.m、rflecoe.m、雅可比矩陣計算JacobiMatrix.m、確定性反演主流程deterministic_inversion.m及初始模型與子波定義initial.txt、seiswavelet.m等完整模塊結(jié)構(gòu)清晰、功能閉環(huán)。已有1193人學(xué)習(xí)下載適用于地震資料處理課程實踐、科研算法驗證或工業(yè)級反演流程快速原型開發(fā)。用戶可直接調(diào)用各函數(shù)模塊理解三參數(shù)耦合反演機(jī)制結(jié)合初始模型與實測道集開展迭代優(yōu)化獲得高分辨率地下彈性參數(shù)剖面。 在油氣勘探圈里摸爬滾打這么多年疊前反演始終是讓我覺得最有嚼頭的一項技術(shù)。很多做地震解釋的朋友一聽到“反演”兩個字就頭大覺得那是反演工程師的活兒自己只要拿到最終屬性體就行。但說句實在話如果你不理解疊前反演里的那些參數(shù)到底在反什么、為什么有時候結(jié)果靠譜有時候完全不靠譜那你拿到手的東西很可能就是一個披著華麗外衣的“陷阱”。我剛開始接觸three-parameter-inversion三參數(shù)反演的時候也走過不少彎路。最早自己以為疊前反演就是比疊后反演多了一個角度道集的事無非就是“把三個角度體的波阻抗都反出來”。真正做下去才發(fā)現(xiàn)這里面涉及的概念遠(yuǎn)沒有那么簡單從Zoeppritz方程出發(fā)到Aki-Richards近似再到截距、梯度、曲率屬性的提取和縱橫波阻抗、密度的同步反演每一步都對輸入數(shù)據(jù)的品質(zhì)、約束條件的合理性、初始模型的精度極其敏感。一個環(huán)節(jié)沒做好后面全白搭。這篇文章我打算把三參數(shù)疊前反演的整個邏輯鏈條和技術(shù)要點掰開揉碎講清楚包括它和AVO之間的關(guān)系、三種參數(shù)各自的地質(zhì)含義、實際反演流程中那些最容易出問題的環(huán)節(jié)以及我在Petrel、Hampson-Russell等軟件里反復(fù)實操后總結(jié)出來的一些避坑經(jīng)驗。內(nèi)容偏干貨盡量說人話希望能讓正在研究疊前反演或者正準(zhǔn)備從疊后反演往疊前轉(zhuǎn)的朋友少走點彎路。1. 三參數(shù)反演到底在反什么從Zoeppritz方程到A、B、C屬性1.1 為什么是“三參數(shù)”A、B、C背后的物理意義要聊三參數(shù)反演必須先搞清楚它反演的這三個參數(shù)到底是什么。很多初學(xué)者以為三參數(shù)就是縱波阻抗、橫波阻抗和密度這個理解不能算錯但不準(zhǔn)確因為這是“最終輸出”而反演過程的“中間產(chǎn)物”或者說“核心參數(shù)”是AVO屬性中的截距Intercept、梯度Gradient和曲率Curvature也就是常說的A、B、C三項。這三項來源于Zoeppritz方程。這個方程從彈性波理論出發(fā)描述了平面縱波入射到兩個彈性介質(zhì)分界面時反射縱波、透射縱波、反射橫波和透射橫波的振幅分配關(guān)系。完整版的Zoeppritz方程形式極其復(fù)雜涉及縱橫波速度、密度、入射角、透射角等多個變量直接用來做反演幾乎是不可能的數(shù)值穩(wěn)定性也極差。所以實際工作中用的都是各種近似公式最常用的就是Aki-Richards近似和Shuey近似。Shuey在1985年給出的三項近似式可以寫成這樣R(θ) ≈ A B·sin2θ C·tan2θ·sin2θ其中A就是AVO截距它代表零偏移距也就是垂直入射時的反射系數(shù)本質(zhì)上反映的是縱波阻抗差異B是AVO梯度它描述的是反射振幅隨入射角增大而減小的速率對縱橫波速度比的相對變化極其敏感是判別含氣異常的關(guān)鍵C是AVO曲率它控制大入射角時的振幅變化行為和密度項的貢獻(xiàn)高度相關(guān)。三參數(shù)反演就是利用不同偏移距或者不同角度地震道集上振幅的變化特征把這三個參數(shù)估計出來。有了A、B、C之后你可以再通過一組線性變換公式推導(dǎo)出縱波阻抗反射系數(shù)、橫波阻抗反射系數(shù)和密度反射系數(shù)再經(jīng)過遞推或約束反演最終得到縱波阻抗、橫波阻抗和密度三個絕對屬性體。這就是三參數(shù)反演最核心的技術(shù)邏輯。1.2 從A、B、C到巖性參數(shù)這些屬性到底有什么用很多解釋人員第一次拿到三參數(shù)反演結(jié)果時會問你給我縱波阻抗、橫波阻抗、密度我怎么用直接拿縱波阻抗看儲層疊后反演已經(jīng)做得很好了干嘛還要費(fèi)勁做疊前答案隱藏在“組合”里。單看縱波阻抗確實有時候和疊后反演差別不大但一旦你有了橫波阻抗和密度就可以計算出大量組合屬性比如縱橫波速度比Vp/Vs這是識別巖性和流體的利器。常規(guī)泥巖的Vp/Vs通常在1.8到2.2之間含氣砂巖往往顯著降低到1.5以下碳酸鹽巖則通常在1.7到1.9之間。這個參數(shù)比單純看縱波阻抗的區(qū)分度要高得多。泊松比Poissons Ratio由縱橫波速度直接換算而來是巖石力學(xué)和流體識別的經(jīng)典參數(shù)。含氣儲層的泊松比異常降低這個特征在很多盆地都得到了驗證。拉梅參數(shù)λρ和μρλρ反映的是流體的不可壓縮性μρ反映的是巖石骨架的剛性。這兩個參數(shù)組合做交會圖對巖性和流體的識別能力比單一彈性參數(shù)要強(qiáng)得多。流體因子ΔF基于Fatti近似的流體因子表達(dá)式可以直接指示流體異常。這些衍生屬性才是三參數(shù)反演真正值錢的地方。疊后反演你只能拿到一個縱波阻抗相當(dāng)于只看了一張黑白照片三參數(shù)反演則像是做了一次CT掃描能把巖石骨架和孔隙流體的影響相對分離出來。這也是為什么現(xiàn)在做儲層預(yù)測和流體檢測行業(yè)里都傾向于上疊前反演而不僅僅是疊后反演。1.3 三參數(shù)反演和AVO反演是不是一回事這個問題我經(jīng)常被問到包括很多干了五六年解釋工作的人也未必分得很清。我的理解是AVO反演和三參數(shù)反演在數(shù)學(xué)本質(zhì)上是一回事——都是利用AVO響應(yīng)特征來反演彈性參數(shù)但側(cè)重點和實現(xiàn)路徑略有不同。AVO反演更側(cè)重于利用A、B屬性以及它們的交會圖來進(jìn)行異常檢測核心工作是做AVO屬性分析和屬性交會通常不強(qiáng)調(diào)得到絕對彈性參數(shù)體。經(jīng)典的AVO分類一類到四類含氣砂巖就是基于截距和梯度的組合特征來劃分的。而三參數(shù)反演通常指的是“AVO加權(quán)疊加反演”或者“疊前同步反演”它不只是停留在屬性提取層面而是進(jìn)一步通過反演框架把地震數(shù)據(jù)、測井?dāng)?shù)據(jù)、地質(zhì)約束綜合在一起直接產(chǎn)出縱波阻抗、橫波阻抗和密度等彈性參數(shù)體。這個流程中A、B、C屬性的提取只是一個中間步驟后面還有大量的模型約束和優(yōu)化迭代過程。所以我一般跟別人解釋時說AVO反演是“看圖說話”從反射振幅特征看異常三參數(shù)反演是“解方程求解”從反射振幅特征反推彈性參數(shù)場。兩者同根同源但目標(biāo)不同三參數(shù)反演對數(shù)據(jù)和流程的要求更高結(jié)果價值也更大。2. 反演之前最重要的事道集品質(zhì)和測井約束決定反演上限2.1 疊前道集上的每一次馬虎最后都會報復(fù)在你的結(jié)果上這里我要非常嚴(yán)肅地說一句三參數(shù)反演這個技術(shù)對輸入道集數(shù)據(jù)的要求是你無法想象的嚴(yán)格。疊后反演對道集品質(zhì)的容忍度相對較高因為疊加本身就在一定程度上壓制了噪聲但疊前反演完全不一樣。它是利用不同偏移距上的振幅差異來做擬合的如果道集上有剩余時差、AVO振幅畸變、多次波殘留、各向異性效應(yīng)沒有校正干凈那同一個反射點的不同角度振幅就不滿足Zoeppritz方程描述的物理規(guī)律反演時硬擬合出來的A、B、C屬性必然被污染。我見過太多失敗的疊前反演項目最后排查來排查去問題都出在道集預(yù)處理上。比如有一個海上工區(qū)原始道集的遠(yuǎn)偏移距振幅明顯偏弱且存在頻率衰減團(tuán)隊沒有做遠(yuǎn)偏移距補(bǔ)償就進(jìn)了反演流程結(jié)果反出來的密度體在遠(yuǎn)道貢獻(xiàn)大的區(qū)域出現(xiàn)了一堆不合理的“低速異常”讓人誤以為是含氣響應(yīng)。后來補(bǔ)做了遠(yuǎn)偏移距Q補(bǔ)償和振幅恢復(fù)密度體的假象就消失了。所以我的建議是疊前反演之前至少要做以下幾項道集處理去噪包括異常振幅衰減、隨機(jī)噪聲壓制、多次波衰減尤其是多次波對疊前反演的破壞力極大。剩余時差校正道集上同相軸必須拉平到同一個零相位對齊基準(zhǔn)。剩余時差會引起振幅相位混疊破壞AVO趨勢。振幅隨偏移距角度變化的補(bǔ)償包括幾何擴(kuò)散校正、吸收衰減補(bǔ)償、透射損失補(bǔ)償?shù)?。這些處理在生產(chǎn)地震處理流程中通常做了但不同偏移距之間的殘余振幅差異需要專門檢查。各向異性校正如果工區(qū)存在明顯的VTI或HTI各向異性不做校正的話遠(yuǎn)偏移距道集的振幅和走時都會畸變直接導(dǎo)致大角度數(shù)據(jù)的AVO響應(yīng)失真。切除和角度范圍控制遠(yuǎn)偏移距的拉伸畸變數(shù)據(jù)、低信噪比區(qū)域要果斷切除寧可少用角度也不要用壞數(shù)據(jù)。這些處理做完之后最好做一個AVO正演模擬驗證。具體做法是在關(guān)鍵井位處利用測井曲線做Zoeppritz正演得到理論合成角道集再將實際地震道集與合成道集進(jìn)行對比。如果兩者在目的層段整體AVO趨勢一致說明道集品質(zhì)可以進(jìn)入反演如果差異很大說明還有處理問題沒解決。2.2 測井曲線整理橫波預(yù)測和密度曲線的鬼門關(guān)三參數(shù)反演需要用測井曲線來建立初始模型和約束條件其中最關(guān)鍵也最容易出問題的是橫波速度曲線和密度曲線。橫波測井并不是每口井都做的國內(nèi)很多老井只有縱波和密度。這時候就需要利用經(jīng)驗公式或者巖石物理模型來預(yù)測橫波速度。最常用的經(jīng)驗關(guān)系包括Castagna泥巖線、Gardner公式等但這些公式都有特定的適用區(qū)域。更可靠的做法是利用工區(qū)內(nèi)已有的實測橫波井通過多元回歸建立區(qū)域性的Vp-Vs關(guān)系再做橫波預(yù)測。我自己的經(jīng)驗是橫波預(yù)測的誤差控制必須嚴(yán)格特別是在目的層段。如果預(yù)測橫波和實際橫波偏差超過5%反演出來的Vp/Vs和泊松比屬性就會出現(xiàn)系統(tǒng)性偏差導(dǎo)致含油氣檢測結(jié)論不可靠。做橫波預(yù)測時一定要用至少一到兩口實測橫波井做盲測驗證看到吻合度達(dá)到90%以上才能批量外推。密度曲線的問題相對隱蔽。密度測井受井壁垮塌、泥餅、侵入帶影響很大經(jīng)常出現(xiàn)非真實的地層密度響應(yīng)。某些井段的密度曲線明顯跳變或者數(shù)值異常比如小于1.0 g/cc這種曲線直接用來做反演會直接污染低角度道集的權(quán)重。密度曲線要做環(huán)境校正同時要有經(jīng)驗地進(jìn)行深部區(qū)域編輯和擬合。還有一個很多人容易忽略的點測井曲線的極性、單位、深度基準(zhǔn)必須和地震數(shù)據(jù)完全一致。井震標(biāo)定做不好子波提取就做不準(zhǔn)反演結(jié)果是空中樓閣。我后來在Petrel里做地震反演項目時第一步永遠(yuǎn)是重新做井震標(biāo)定強(qiáng)制統(tǒng)一極性把時深關(guān)系做得嚴(yán)絲合縫才敢往下走。2.3 子波提取疊前反演的成敗手子波提取在疊前反演中的重要性怎么強(qiáng)調(diào)都不為過。疊后反演的子波提取相對簡單從井旁道提取統(tǒng)計子波就行疊前反演要考慮不同角度地震數(shù)據(jù)的子波一致性而且要在每個部分角度范圍或者每個入射角子集上分別提取子波。理論上Zoeppritz方程正演是確定性的如果子波是零相位的那么合成道集應(yīng)該和實際道集一致。但實際上處理流程中不同偏移距吸收衰減不同不同角度的子波往往存在頻率和相位差異導(dǎo)致實際的角道集上的子波是變化的。如果反演時用同一個子波去套不同角度的數(shù)據(jù)就會出現(xiàn)系統(tǒng)誤差。我的做法是將角道集分成幾個疊加范圍比如近、中、遠(yuǎn)、超遠(yuǎn)分別提取子波并檢查不同角度子波的振幅譜和相位譜是否穩(wěn)定。如果某個角度的子波與其他角度差異太大通常說明那個角度的地震數(shù)據(jù)品質(zhì)有問題需要返回去重新處理。子波相位一致性對A、B、C屬性提取的影響巨大。我記得有一次反演結(jié)果中梯度屬性與測井曲線完全不匹配處處都是“亮點”最后發(fā)現(xiàn)是遠(yuǎn)角度道集的子波相位比近角度道集偏了約30度導(dǎo)致梯度項被虛擬放大。后來做了相位校正G屬性馬上正常了。所以任何一個認(rèn)真做疊前反演的人都應(yīng)該把子波提取與一致性分析當(dāng)成關(guān)鍵質(zhì)量控制環(huán)節(jié)這一步?jīng)]做好后面所有反演結(jié)果的可信度都要打上問號。3. 實操流程拆解在Petrel中跑通三參數(shù)疊前同步反演3.1 Petrel地震反演模塊的定位與核心路線既然提到了Petrel我多說幾句。Petrel是斯倫貝謝旗下的一體化油藏描述平臺它內(nèi)置的地震反演模塊包括確定性反演和基于地質(zhì)統(tǒng)計學(xué)的隨機(jī)反演其中就包含疊前同步反演的功能。在Petrel中做三參數(shù)疊前反演通常有兩種可選路徑第一條路徑是“AVO屬性提取→屬性反演”先在Petrel中利用角道集提取A、B、C屬性體再對屬性體進(jìn)行模擬退火或者地質(zhì)統(tǒng)計學(xué)反演得到橫波阻抗、縱波阻抗和密度模型。第二條路徑是“同步反演”直接利用不同角度部分疊加數(shù)據(jù)體和測井?dāng)?shù)據(jù)在一個整體反演框架下同時反演縱波阻抗、橫波阻抗和密度。這條路徑在Hampson-Russell軟件里更成熟Petrel也支持但對參數(shù)設(shè)置和初始模型要求更高。我個人的使用感受是Petrel做三參數(shù)反演的優(yōu)勢在于它的一體化程度高。解釋層位、斷層模型、相建模、屬性建模都在同一個平臺構(gòu)建三維地質(zhì)框架非常順暢。如果你已經(jīng)習(xí)慣在Petrel里做解釋和地質(zhì)建模那直接在這個平臺里做疊前反演的效率會很高。但這里要提醒一點Petrel的地震反演模塊在國內(nèi)團(tuán)隊的普及度和熟練度整體上不如Hampson-Russell國內(nèi)很多反演工程師習(xí)慣HR或者Jason。如果你之前完全沒有用過Petrel的地震反演模塊建議先找個小工區(qū)練手不要拿重點項目直接上。3.2 數(shù)據(jù)加載與工區(qū)框架構(gòu)建基礎(chǔ)工作決定成敗在Petrel中做三參數(shù)疊前反演的第一步是把所有數(shù)據(jù)裝進(jìn)工區(qū)?？此坪唵蔚@一步的坑非常多。首先是角道集數(shù)據(jù)體的加載。不同處理軟件輸出的角道集格式多種多樣有分角度部分疊加的比如近、中、遠(yuǎn)三個數(shù)據(jù)體也有共角度道集的連續(xù)角度域數(shù)據(jù)比如每5度一個數(shù)據(jù)體的15個數(shù)據(jù)體。Petrel中要明確地震數(shù)據(jù)體的坐標(biāo)系統(tǒng)、測網(wǎng)密度和極性加載之后先做可視化檢查確認(rèn)空間位置和剖面響應(yīng)正確。然后是測井?dāng)?shù)據(jù)的加載。Petrel支持LAS、ASCII、DLIS等格式。加載時需要注意測井曲線的深度單位米還是英尺、曲線名稱不同公司命名的Vp、VP、P-wave、DT等千奇百怪、橫波是否存在、密度曲線是否完整。加載完成后要逐井做曲線完整性檢查尤其是目的層段范圍內(nèi)不能有大的盲區(qū)。接下來是構(gòu)造框架的建立。地震解釋層位是構(gòu)建低頻模型的骨架。在Petrel里解釋層位需要解釋到目的層段的頂?shù)缀蛢?nèi)部主要分界面而且層位的解釋方案要和測井分層的認(rèn)識一致。如果層位解釋有偏差低頻模型就會“跑偏”反演結(jié)果的背景趨勢就會錯誤。這一步也最容易出現(xiàn)“垃圾進(jìn)、垃圾出”的問題。我見過一個項目地震解釋層位在局部斷層附近解釋明顯不合理但直接用了這個層位來建低頻模型結(jié)果反演出的橫波阻抗在斷層附近出現(xiàn)了一長條不合理的低值帶讓人誤以為是斷層封堵形成的氣藏響應(yīng)。后來修正了解釋層位這個問題才消失。所以我的經(jīng)驗是在建立低頻模型之前至少花一個下午的時間把解釋層位和測井分層一張一張檢查確保每個井點處的地震層位和測井分層誤差小于半個子波長度大約5-10米左右這是低頻模型可信的前提。3.3 低頻模型構(gòu)建和約束條件設(shè)置的關(guān)鍵參數(shù)三參數(shù)反演的初始模型本質(zhì)上是把測井資料里的低頻信息向地震工區(qū)外推。為什么需要低頻模型因為地震數(shù)據(jù)本身是帶限的缺乏低頻成分通常低于5-10Hz的信息地震反射數(shù)據(jù)里沒有而測井?dāng)?shù)據(jù)有非常豐富的低頻信息。如果不補(bǔ)充低頻背景反演結(jié)果就會缺乏絕對數(shù)值的約束縱波阻抗、橫波阻抗和密度的絕對值都會“漂移”。在Petrel中構(gòu)建低頻模型時需要設(shè)置以下幾個關(guān)鍵參數(shù)層位約束方式通常選擇“層間插值”即每個地層單元內(nèi)部單獨(dú)插值模擬地層沉積的相對等時性。要避免跨層位插值那會讓不同地層單元的彈性參數(shù)互相混淆。橫向變差函數(shù)用于控制井間插值的橫向連續(xù)性。變差函數(shù)的主方向要和沉積體系方向一致變程通常取一個工區(qū)主河道寬度或者沉積體寬度的量級。變程太小井間模型破碎反演結(jié)果橫向變化劇烈變程太大模型過于平滑會丟失真實的橫向非均質(zhì)性?？v向變差函數(shù)縱向變程和垂向采樣率、地層厚度有關(guān)一般取2-6米。各參數(shù)之間的相關(guān)性設(shè)置縱波阻抗、橫波阻抗、密度之間存在一定的相關(guān)性在建立初始模型時可以設(shè)置這些參數(shù)之間的相關(guān)系數(shù)。但要注意如果相關(guān)系數(shù)設(shè)得過高會強(qiáng)制橫波阻抗和密度跟著縱波阻抗走失去獨(dú)立反演的意義設(shè)得過低約束太弱反演不穩(wěn)定。一般建議縱波阻抗和橫波阻抗相關(guān)系數(shù)設(shè)置在0.6-0.8之間密度獨(dú)立一些。對于Petrel中的約束條件設(shè)置最核心的是每個參數(shù)的取值范圍或者標(biāo)準(zhǔn)差。你可以從測井?dāng)?shù)據(jù)統(tǒng)計出目的層段各參數(shù)的分布范圍然后在反演時設(shè)置合理的上下限。上下限太寬反演會在噪聲中迷失太窄反演結(jié)果會被人為限制在狹窄區(qū)間里出現(xiàn)“假精細(xì)”。關(guān)于這一點我有一個非常深的教訓(xùn)。早年間我做一個碳酸鹽巖工區(qū)的反演把密度上下限設(shè)置得太窄反演結(jié)果密度體幾乎就是初始模型的復(fù)制品完全看不到地震數(shù)據(jù)的貢獻(xiàn)橫波阻抗倒是正常。后來放寬了密度的允許變化范圍密度體才出現(xiàn)了正確的異常響應(yīng)。所以約束條件的設(shè)置要適中不能太保守。3.4 反演參數(shù)與實際運(yùn)行角度子體重構(gòu)與質(zhì)量控制在Petrel中運(yùn)行疊前同步反演時需要定義參與反演的角度子體。通常是將角道集數(shù)據(jù)劃分為三個部分疊加體近角度比如0-15°、中角度15-30°、遠(yuǎn)角度30-45°。子體數(shù)與地質(zhì)目標(biāo)有關(guān)至少要用三個才有能力約束三個參數(shù)。這里有一個數(shù)學(xué)上的直覺值得說清楚。三參數(shù)反演是在擬合R(θ) A B·sin2θ C·tan2θ·sin2θ其中A主要受近角度約束B受中等角度約束C強(qiáng)烈依賴大角度數(shù)據(jù)。如果工區(qū)遠(yuǎn)角度信噪比很差大角度數(shù)據(jù)被切除了那么C項的約束就非常弱密度反演的不確定性很大。實際運(yùn)行前我會先在Petrel里生成AVO截距和梯度屬性和質(zhì)量控制面板對比關(guān)鍵井的AVO響應(yīng)。這個過程很快但能幫你發(fā)現(xiàn)道集中是否存在明顯的振幅異常模式。正式反演時通常會先做一次“快速測試運(yùn)行”用較粗的采樣步長比如4ms跑通全流程檢查結(jié)果是否穩(wěn)定如果穩(wěn)定再加密采樣比如2ms做正式反演。這一招能有效節(jié)省調(diào)試時間。Petrel的疊前同步反演計算量比較大尤其是多個角度子體參與時一個中等規(guī)模的工區(qū)幾百平方公里常常需要跑數(shù)小時到十幾小時。我一般建議夜間運(yùn)行第二天早上來檢查結(jié)果。運(yùn)行結(jié)束后先看幾條過井剖面。檢查合成地震記錄與原始地震數(shù)據(jù)的殘差殘差越小說明擬合越好。同時對比反演出的縱波阻抗、橫波阻抗、密度曲線與測井原始曲線的匹配程度。如果曲線匹配良好下一步才做平面成圖和交會圖分析。3.5 結(jié)果驗證和屬性導(dǎo)出別急著下結(jié)論反演完成之后千萬不要急著出圖寫報告。三參數(shù)反演的結(jié)果驗證是必須的一步而且要驗證的不是一條曲線而是多個維度的合理性。第一個驗證是井點驗證選擇參與反演的井檢查反演結(jié)果曲線與原始測井曲線的匹配度。再選擇一口或多口不參與反演的檢驗井做“盲井驗證”。盲井驗證通過說明反演結(jié)果在井間具有預(yù)測能力如果真的盲井處反演結(jié)果和測井曲線嚴(yán)重不匹配就要回頭審查初始模型、約束條件或者道集品質(zhì)。第二個驗證是地質(zhì)合理性驗證把反演得到的縱橫波速度比或泊松比體與已知油氣層分布進(jìn)行對比。比如工區(qū)內(nèi)有已鉆井證實的氣層在含氣區(qū)應(yīng)該看到Vp/Vs明顯降低的響應(yīng)。如果連已知井的已知油氣層都看不到響應(yīng)反演結(jié)果的可靠性就很值得懷疑。第三個驗證是空間分布合理性對反演體做切片和剖面展示檢查是否存在條帶狀、串珠狀的不合理異常體。特別是斷層附近、地震數(shù)據(jù)采集邊界處、測井資料稀疏區(qū)經(jīng)常會出現(xiàn)人為假象這些位置的結(jié)果要謹(jǐn)慎解釋。一旦驗證通過就可以把反演結(jié)果導(dǎo)出為體文件如SEGY格式或者ZGY格式供后續(xù)地質(zhì)建模或者儲層預(yù)測使用。在Petrel中可以直接屬性計算得到Vp/Vs、拉梅參數(shù)、流體因子等衍生屬性體這些屬性體還可以進(jìn)一步作為油氣檢測的直接證據(jù)。4. 實戰(zhàn)中那些讓人腦殼疼的問題常見故障與排查方案4.1 密度反演總是不受控大角度道集是命門三參數(shù)反演里密度參數(shù)的反演是最讓人頭疼的。前面說過密度項主要由C參數(shù)控制而C參數(shù)主要依賴大角度道集。但在實際地震數(shù)據(jù)中大角度道集通常意味著低信噪比、動校正拉伸畸變嚴(yán)重、各向異性影響顯著能“干凈”地到達(dá)35度以上偏移距的數(shù)據(jù)少之又少。如果反演出的密度體要么高到離譜、要么低到不可思議或者密度剖面出現(xiàn)“千層餅”樣的高頻變化你首先要檢查的是參與反演的角度子體是否包含了足夠多的大角度信息。我有一個典型的排查實驗先在近-中-遠(yuǎn)三個角度子體參與反演時密度反演結(jié)果和測井密度嚴(yán)重不匹配然后我把遠(yuǎn)角度子體進(jìn)一步細(xì)化替換為超遠(yuǎn)角度35-45度子體參與反演雖然超遠(yuǎn)角度信噪比略低但密度反演的分辨率和準(zhǔn)確性顯著提升。此外密度反演弱約束時的非唯一性問題也是常見問題。在反演時必須給密度設(shè)置合理的先驗約束否則結(jié)果極度發(fā)散。在Petrel中你可以定義密度和縱波阻抗之間的背景關(guān)系利用Gardner公式或工區(qū)統(tǒng)計關(guān)系作為先驗約束既能增加反演穩(wěn)定性又能保留異常信息。4.2 橫波速度估算不準(zhǔn)Vp/Vs屬性全線失真這本問題的根源大概率不是反演本身而是橫波速度預(yù)測的精度不夠。很多工區(qū)沒有實測橫波只能靠經(jīng)驗公式預(yù)測如果預(yù)測的橫波速度在目的層段與實際地層相差超過5%那么反演得到的Vp/Vs就會偏移油氣檢測結(jié)論就會失真。排查思路很簡單看反演結(jié)果在已鉆井處的Vp/Vs值與井上實測Vp/Vs或用實測橫波計算的差多少。如果偏差大直接把這個差值標(biāo)在圖上看是否有系統(tǒng)性的偏移。如果各井偏差方向一致說明橫波預(yù)測訓(xùn)練數(shù)據(jù)本身有問題如果偏差方向不一致、時正時負(fù)說明變差函數(shù)或者約束條件的非均質(zhì)刻畫有問題。我對橫波預(yù)測的建議永遠(yuǎn)是不要偷懶。如果工區(qū)內(nèi)有超過5口實測橫波井一定要建立工區(qū)專屬的Vp-Vs回歸公式而且要用未參與回歸的井做盲測。如果沒有足夠的實測橫波井可以考慮基于巖石物理模型的干巖框架法比如Xu-White模型雖然建模過程復(fù)雜但從機(jī)理上更可靠尤其是復(fù)雜巖性儲層。4.3 反演結(jié)果橫向不連續(xù)、出現(xiàn)大量假亮點怎么辦如果反演結(jié)果在中淺層出現(xiàn)大量零星分布的“假亮點”或者橫向連續(xù)性極差最可能的元兇是子波空間不穩(wěn)定性。前面說了不同角度的子波如果相位或頻率不一致反演時就會產(chǎn)生虛假的A、B、C異常而A、B異常又會被映射到縱波阻抗和橫波阻抗中造成反射剖面上到處是亮點。解決思路有兩個方向第一嚴(yán)格檢查各角度子波的統(tǒng)一性要求所有角度的地震數(shù)據(jù)經(jīng)過處理之后具有一致的子波特征如果做不到就用匹配濾波的方式把不同角度數(shù)據(jù)統(tǒng)一到相同子波第二在反演時引入橫向正則化約束讓相鄰道之間的反演結(jié)果具有一定連續(xù)性。在Petrel中你可以在反演參數(shù)里設(shè)置橫向平滑約束的權(quán)重。但注意這個約束不是越大越好。過大的橫向平滑約束會讓結(jié)果喪失橫向分辨率小斷層和河道邊界會被抹掉。我通常的做法是把橫向平滑權(quán)重從小往大調(diào)每調(diào)一次跑一個快速測試觀察目的層段的橫向連續(xù)性和細(xì)節(jié)保留之間的平衡點。4.4 快速排查清單按順序排除問題如果反演結(jié)果有嚴(yán)重問題我自己習(xí)慣按下面的順序排查效率和準(zhǔn)確性都還不錯排查次序檢查項主要現(xiàn)象常見對策1原始道集品質(zhì)AVO趨勢異常、遠(yuǎn)道信噪比差返回預(yù)處理做剩余校正、切除、去噪2井震標(biāo)定質(zhì)量井旁道殘差大、合成記錄匹配差重新標(biāo)定、校正時深、檢查極性3子波一致性不同角度子波相位譜差異大統(tǒng)一子波、相位校正、匹配濾波4橫波速度精度井點Vp/Vs偏差大重新構(gòu)建Vp-Vs關(guān)系、巖石物理預(yù)測5低頻模型約束背景趨勢偏離地質(zhì)認(rèn)識修正層位框架、調(diào)整變差函數(shù)6反演參數(shù)尺度密度發(fā)散或過于平滑調(diào)整約束范圍、松緊合適按照這個順序排查絕大部分三參數(shù)反演的翻車事故都能找到病根。反過來如果一上來就懷疑反演算法本身有問題往往會走很多彎路。絕大多數(shù)情況下算法本身沒有錯是喂給算法的數(shù)據(jù)或者約束不靠譜。5. 一些更進(jìn)階的思考三參數(shù)反演的前沿方向和適用邊界這些年疊前反演技術(shù)本身也在快速演進(jìn)。傳統(tǒng)三參數(shù)反演依賴疊前道集的分方位/分角度疊加對數(shù)據(jù)品質(zhì)要求極高。而現(xiàn)在比較前沿的技術(shù)比如基于深度學(xué)習(xí)的三參數(shù)直接反演、基于全波形反演FWI的彈性參數(shù)建模、多波多分量聯(lián)合反演等都在試圖突破傳統(tǒng)反演的瓶頸。但作為一線工程師我的態(tài)度一直比較務(wù)實在絕大多數(shù)工區(qū)條件下傳統(tǒng)三參數(shù)反演依然是性價比最高的選擇只要把數(shù)據(jù)品質(zhì)關(guān)、子波一致性和約束條件這幾個核心環(huán)節(jié)做到位結(jié)果可靠性已經(jīng)能夠滿足儲層預(yù)測和流體檢測的需求。深度學(xué)習(xí)和FWI雖然前景廣闊但它們的投入成本、算力要求、數(shù)據(jù)需求都遠(yuǎn)遠(yuǎn)高于傳統(tǒng)方法在常規(guī)勘探生產(chǎn)任務(wù)中很難全面鋪開。還有一個要提醒的點是工區(qū)適用性。三參數(shù)反演也不是萬能的。在復(fù)雜斷塊區(qū)、薄互層儲層、低信噪比地區(qū)反演結(jié)果的不確定性會大幅增加。在這些地區(qū)必須結(jié)合地質(zhì)認(rèn)識做人工干預(yù)有些情況下甚至不得不放棄定量反演退回定性AVO檢測的老路上來。所以我對想學(xué)三參數(shù)反演的朋友的建議是先啃透Zoeppritz方程及其近似的數(shù)學(xué)基礎(chǔ)再踏踏實實做一兩個工區(qū)的全流程實操最后一定要學(xué)會利用正演來驗證反演。正演是你最忠實的老師它能幫你理解地震數(shù)據(jù)中哪些響應(yīng)對應(yīng)的是真實地質(zhì)信息哪些是被處理流程改造過的偽影?；氐介_頭說的那句話三參數(shù)反演不是黑盒子魔術(shù)它是一套邏輯嚴(yán)密、但對每一步輸入都極其敏感的技術(shù)體系。你給它好數(shù)據(jù)、好約束、好模型它回報你可靠的彈性參數(shù)體你給它垃圾數(shù)據(jù)、敷衍的初始模型、馬虎的質(zhì)量控制它也會毫不客氣地還你一堆看似無誤其實全是坑的“精美”成果。做反演本質(zhì)上是在跟物理規(guī)律和數(shù)據(jù)品質(zhì)賽跑而這份嚴(yán)謹(jǐn)與耐心才是這個領(lǐng)域最值錢的經(jīng)驗。本文還有配套的精品資源點擊獲取