地震反演方法概述

1、地震反演方法概述地震反演：由地震信息得到地質(zhì)信息的過程。地震反射波法勘探的基礎(chǔ)在于：地下不同地層存在波阻抗差異，當(dāng)?shù)卣鸩▊鞑ビ胁ㄗ杩?差異的地層分界面時，會發(fā)生反射從而形成地震反射波。地震反射波等于反射系數(shù)與地 震子波的褶積，而某界面的法向入射發(fā)射系數(shù)就等于該界面上下介質(zhì)的波阻抗差與波阻 抗和之比。也就是說，如果已知地下地層的波阻抗分布，我們可以得到地震反射波的分 布，即地震反射剖面。即由地層波阻抗剖面得到地震反射波剖面的過程稱為地震波阻抗 正演，反之，由地震反射剖面得到地層波阻抗剖面的過程稱為地震波阻抗反演。疊前反演 主要是指AVC反演，通過AVC反演，可以獲得全部的巖石參數(shù)，女口：巖石密度

2、、 縱橫波速度、縱橫波阻抗、泊松比等。疊前反演與疊后反演的根本區(qū)別在于疊前反演使 用了未經(jīng)疊加的地震資料。多道疊加雖然能夠改善資料的品質(zhì)，提高信噪比，但是另一 方面，疊加技術(shù)是以東校正后的地震反射振幅、波形等特征不隨炮檢距變化的假設(shè)為基 礎(chǔ)的。實際上，來自同一反射點的地震反射振幅在不同炮檢距上是不同的，并且反射波 形也隨炮檢距的變化而發(fā)生變化。這種地震反射振幅、波形特征隨炮檢距的變化關(guān)系很 復(fù)雜，主要原因就在于不同炮檢距的地震波經(jīng)過的地層結(jié)構(gòu)、彈性性質(zhì)、巖性組合等許 多方面都是不同的。疊加破壞了真實的振幅關(guān)系，同時損失了橫波信息。疊前反演通過 疊前地震信息隨炮檢距的變化特征，來揭示巖性和油氣的

3、關(guān)系。疊前反演的理論基礎(chǔ)是 地震波的反射和透射理論。理論上講，利用反射振幅隨入射角的變化規(guī)律可以實現(xiàn)全部 巖性參數(shù)的反演，提取縱波速度、橫波速度、縱橫波速度比、巖石密度、泊松比、體積 模量、剪切模量等參數(shù)。疊后地震剖面相當(dāng)于零炮檢距的自激自收記錄。與疊前反演不同，疊后反演只能得到縱波阻抗。雖然疊后反演與疊前反演想必有很多不足之處，但由于其技術(shù)方法成熟完 備，到目前為止，疊后反演仍然是主流的反演類型，是儲層預(yù)測的核心技術(shù)。介紹幾種疊后反演方法：1）道積分：利用疊后地震資料計算地層相對波阻抗（速度）的直接反演方法。因為它 是在地層波阻抗隨深度連續(xù)可微的條件下推導(dǎo)出來的，因而又稱為連續(xù)反演。原理簡述

4、：上述公式表示，反射系數(shù)的積分正比于波阻抗Z的自然對數(shù)，這是一種簡單的相對波阻抗概念。適用條件及優(yōu)缺點與絕對波阻抗反演相比，道積分的優(yōu)點：1.遞推時累積誤差較?。?.計算簡單，不需要反射系數(shù)標(biāo)定；3.無需鉆井控制，在勘探儲氣即可推廣使用。缺點：1.由于這種方法受到地震固有頻寬的限制，分辨率低，無法適用于薄層解釋的需 要；2.需要地震記錄經(jīng)過子波零相位化處理；3.無法求得地層的絕對波阻抗和絕對速度，不能用于定量計算儲層參數(shù)；4.這種方法在處理過程中不能用地質(zhì)或測井資料對其進(jìn)行 約束控制，因而結(jié)果比較粗略。2）遞推反演方法：根據(jù)反射系數(shù)進(jìn)行遞推計算地層波阻抗或?qū)铀俣龋潢P(guān)鍵在于由原 始地震記錄估算

5、反射系數(shù)和波阻抗，測井資料不直接參入反演，只起到標(biāo)定和質(zhì)量控制 的作用。因此又稱為直接反演。原理簡述：利用以上公式，可以從聲波時差曲線及密度曲線上（沒有密度曲線時可以利用Gard nar公式進(jìn)行換算）選擇標(biāo)準(zhǔn)層波阻抗作為基準(zhǔn)波阻抗，將反褶積得到的反射系數(shù)轉(zhuǎn)為波阻 抗。遞推反演是對地震資料的處理過程，其結(jié)果的分辨率、信噪比以及可靠程度主要依賴于 地震資料本身的品質(zhì)，因此用于反演的地震資料應(yīng)具有較寬的頻帶，較低的噪聲、相對 振幅保持和準(zhǔn)確成像。反演之前，應(yīng)對聲波測井和密度測井曲線進(jìn)行校正。遞推反演的核心技術(shù)在于由地震資料正確的估算地層的反射系數(shù)（或消除地震子波的影響），比較典型的實現(xiàn)方法有：基于地

6、層反褶積方法、稀稀疏脈沖反演法、測井控制地 震反演法、頻域反褶積法等。 地層反褶積法：很據(jù)已有測井資料（聲波和密度）與井旁地震記錄，利用最小平方法估算數(shù)學(xué)意義上的“最佳”子波和反射系數(shù)。優(yōu)點：把子波求解的“欠定”問題變成確 定問題，再井點已有測井段范圍內(nèi)可獲得與測井最吻合的反演結(jié)果；局限性：1.完全忽略了測井誤差和地震噪音，尤其是側(cè)井誤差的客觀存在使“子波”確定更加困難；2.地層反褶積因子的估算是在計算時窗內(nèi)數(shù)學(xué)上的最佳逼近，實際處理范圍與該時窗的不同已超出該方法的適用范圍，即便是在井點位置，得到的反演結(jié)果已不可能是“誤差最小” 稀疏脈沖反演法：基于稀疏脈沖反褶積基礎(chǔ)上的遞推反演法，主要包括最

7、大似然反褶積、L1模反褶積和最小熵反褶積。這類方法針對地震記錄的欠定問題，提出了地層反射系數(shù)為一系列疊加于高斯背景下的強(qiáng)軸的基本假設(shè)，在此條件下以不同的方法估算地下“強(qiáng)”反射系數(shù)和地震子波。優(yōu)點：無需考慮鉆井資料，直接由地震記錄計算反射系數(shù)， 實現(xiàn)遞推反演；其缺陷在于很難得到與測井曲線相吻合的最終結(jié)果。 基于頻域反褶積與相位校正的遞推反演法：從方法實現(xiàn)上回避了計算子波或反射系數(shù)的欠定問題，以井旁反演結(jié)果與實際測井曲線的吻合程度作為參數(shù)優(yōu)選的基本依據(jù)，從 而保證了反演資料的可信度和可解釋性，是遞推反演的主導(dǎo)方法。其主要技術(shù)關(guān)鍵有： 恢復(fù)地層反射系數(shù)振幅譜的頻域反褶積、使井旁反演道與測井最佳吻合的

8、相位校正以及反映地層波阻抗變化趨勢的低頻模型技術(shù)。適用條件及優(yōu)缺點應(yīng)用領(lǐng)域較寬：勘探初期鉆井較少的情況下，通過反演資料進(jìn)行巖相分析確定地層的沉 積體系，根據(jù)鉆井揭示儲集層特征進(jìn)行橫向預(yù)測，確定評價井位；開發(fā)前期，在儲層較 厚的條件下，遞推反演資料可為地質(zhì)建模提供較可靠的構(gòu)造、厚度和物性信息，優(yōu)化方 案設(shè)計；在油藏檢測階段，通過時延地震反演速度差異分析，可幫助確定儲集層壓力、 物性的空間變化，進(jìn)而推斷油氣前緣。缺點：1.由于受地震頻帶寬度的限制，遞推反演資料的分辨率相對較低，不能滿足薄儲 集層研究的需要；2.算法相對復(fù)雜，在具體實現(xiàn)過程中存在一些難點：反射剖面的極 性問題：地震反射波的極性是正還

9、是負(fù)直接影響反演波阻抗后速度變高還是變低；標(biāo) 定問題：地震反演中對反射系數(shù)的標(biāo)定，通常是根據(jù)井中反射系數(shù)來標(biāo)定反褶積后的振 幅值的。但是，求波阻抗是一個積分的過程，反褶積后的地震道振幅實際上還不是反射 系數(shù)，而是相對于反射系數(shù)再褶積的一個剩余子波。這個剩余子波一般在淺層主頻高些，在深層主頻低些。頻率低的波積分后數(shù)值偏大，會使深層產(chǎn)生較大的波阻抗值。因此在 標(biāo)定時，除了要考慮時變的振幅因素外，還要考慮時變的主頻變化；低頻分量的補(bǔ)償 問題：在有井的情況下，以井為控制，能夠得到該點的低頻分量，但是在井與井之間低 頻分量的內(nèi)插又是一個難題，簡單的內(nèi)插只有在地層等厚且產(chǎn)狀平穩(wěn)時才行。即使利用 地層產(chǎn)狀起

10、伏控制內(nèi)插，還有高低、頻帶的一個銜接問題，因為低頻成分一定要與子波 的譜“互補(bǔ)”。在無井區(qū)，波阻抗反演往往要從疊加速度譜中提取低頻分量，又存在速 度譜的質(zhì)量與分辨率問題。3）基于模型反演法基本思路：先建立一個初始地層波阻抗模型，然后由此模型進(jìn)行正演，求得地震合成記錄，將合成地震記錄與實際地震記錄相比較，根據(jù)比較結(jié)果修改地下波阻抗模型的速度、密度、深度值及子波，再正演求取合成地震記錄，與實際地震記錄進(jìn)行比較，繼續(xù)修改 波阻抗模型，如此多次反復(fù)，從而不斷地通過迭代修改，直至合成地震記錄與實際地震 記錄最接近，最終得到地下的波阻抗模型?；谀Ｐ头囱莘椒ㄖ饕幸韵聨追N： 測井約束反演：是目前生產(chǎn)中廣泛

11、采用的基于模型的地震反演方法。將地震與測井有 機(jī)結(jié)合，突破傳統(tǒng)意義上地震分辨率的限制，理論上可以得到與測井資料相同的分辨率（問題實質(zhì)在于怎么減少多解性）。 地震巖性模擬： 西方地球物理公司的 SLIM （地震巖性模擬） 程序也可以把地震剖面反 演成很詳細(xì)的波阻抗剖面。它是將模型正演的結(jié)果與實際地震記錄作比較，然后根據(jù)比較的結(jié)果，反復(fù)修改地下波阻抗模型的速度、密度和深度數(shù)值（同時也修改子波），從 而不斷的通過迭代修改，找到一個詳細(xì)的地下波阻抗模型。該方法避免了一般反褶積算 法最子波的最小相位假設(shè)，也不需要假設(shè)反射系數(shù)是白噪，但是反演結(jié)果只有一種可能 的解答，不一定真實。 廣義線性反演：通過模型正

12、演與實際地震剖面作比較，根據(jù)誤差的情況，在最小二乘 意義上，或者在誤差絕對值之和最小的意義上，最佳逼近實際數(shù)據(jù)，從而迭代反復(fù)修改 模型，直到滿意為止。 多道反演法：分為無井多道反演和有井多道反演。有井多道反演是在無井多道反演的 基礎(chǔ)上，結(jié)合已知井的資料建立初始模型和提取子波，將多道反演的思路應(yīng)用于寬帶約 束反演，該方法保持了寬帶高分辨率的特點，與測井的吻合性較好，可以較好的壓制隨 即噪聲，但是最規(guī)則噪聲不適用。 地質(zhì)統(tǒng)計學(xué)反演：在地質(zhì)和地層模型中對一個三維地震數(shù)據(jù)體進(jìn)行轉(zhuǎn)換，得到一些儲層尺度的波阻抗數(shù)據(jù)體， 并且通過這些三維數(shù)據(jù)體進(jìn)行統(tǒng)計學(xué)計算，來量化其不確定性。地質(zhì)統(tǒng)計學(xué)反演首先在地震時間域

13、內(nèi)建立儲層的地質(zhì)模型，層面由拾取的地震層位決 定，地層網(wǎng)格的結(jié)構(gòu)（上超，剝蝕）取決于地質(zhì)情況，并將井位出的原始地質(zhì)波阻抗曲 線放置于地層網(wǎng)格內(nèi)。利用井和地震數(shù)據(jù)來決定地質(zhì)統(tǒng)計學(xué)參數(shù)，然后開始地質(zhì)學(xué)統(tǒng)計 反演過程。模擬過程沿著一個隨機(jī)路徑進(jìn)行，并且在每一個隨機(jī)拉伸道位置，通過序貫 高斯模擬產(chǎn)生波阻抗值，并計算出相應(yīng)的反射系數(shù)。反射系數(shù)與子波褶積后，與實際地 震資料擬合最好的波阻抗道被保留，并且與井?dāng)?shù)據(jù)及以前的模擬波阻抗道合并。適用于 各類復(fù)雜儲層的地震預(yù)測和描述，尤其鉆井資料較多、需要進(jìn)行精細(xì)儲層描述地區(qū)；該 算法運算量大，速度慢。 波阻抗多尺度反演：采用小波變化，把目標(biāo)函數(shù)分成不同尺度的分量，

14、根據(jù)不同尺度上目標(biāo)函數(shù)的特征逐步搜索全局最小點。一般情況下，大尺度上，目標(biāo)函數(shù)的極值點較 少，且分得開，用通常的線性化方法很容易所搜到該尺度。在相對較小尺度上，目標(biāo)函 數(shù)極值點較多，直接尋找全局極值點比較困難。但是，如果以大尺度上搜索到的總體背 景上“全局極小點”為起始點，則能很容易地在其附近搜索到對應(yīng)尺度上的“全局極小點”。最后，當(dāng)尺度降至目標(biāo)函數(shù)的原始尺度時，對應(yīng)搜索出的“全局極小點”就是目 標(biāo)函數(shù)的全局最小點。優(yōu)點：反演穩(wěn)定，反演結(jié)果不受選定的初始點影響，從而避免其 后的反演落入錯誤的領(lǐng)域，并且收斂速度加快。缺點：對目標(biāo)函數(shù)由大到小的多尺度分解過程中，總是假定上一尺度（較大尺度）在迭代終止點就是下一尺度的“全局極