基于分頻數(shù)據(jù)的地震反演方法及其裝置的制造方法
【技術(shù)領域】
[0001] 本發(fā)明設及地球物理勘探技術(shù)領域，特別設及一種基于分頻數(shù)據(jù)的地震反演方法 及其裝置。
【背景技術(shù)】
[0002] 作為地下儲層空間展布形態(tài)描述的常用技術(shù)，地震反演長期在生產(chǎn)應用中發(fā)揮重 要作用。地震反演是利用地表觀測地震資料，W已知地質(zhì)規(guī)律和鉆井、測井資料為約束，對 地下巖層空間結(jié)構(gòu)和物理性質(zhì)進行成像的過程。而分頻反演技術(shù)是目前地震反演中使用的 一項新技術(shù)，它通過對地震數(shù)據(jù)進行分頻處理，提取不同頻率數(shù)據(jù)的地震屬性體，利用不同 頻率數(shù)據(jù)的地震屬性體擬合目標曲線，進而得到最終的儲層預測結(jié)果。本發(fā)明的發(fā)明人發(fā) 現(xiàn)常規(guī)的分頻反演技術(shù)利用地震屬性體擬合測井波阻抗曲線，由于地震屬性體只反映地震 數(shù)據(jù)單一屬性特征，其不能完全反映振幅、波形、相位等特征，進而無法反映地下具體地質(zhì) 情況，只能對地下儲層進行定性預測，常規(guī)的分頻反演技術(shù)的定性預測使其喪失了常規(guī)地 震反演定量化預測的功能，其具有一定的局限性。

【發(fā)明內(nèi)容】

[0003] 為了克服現(xiàn)有技術(shù)的上述缺陷，本發(fā)明實施例提供了一種基于分頻數(shù)據(jù)的地震反 演方法及其裝置，其能對地下儲層空間進行定量化預測。
[0004] 本發(fā)明實施例的具體技術(shù)方案是： 陽〇化]一種基于分頻數(shù)據(jù)的地震反演方法，它包括W下步驟：基于收集的區(qū)域內(nèi)的地震 和測井數(shù)據(jù)獲得地震反演初始模型和不同頻率的地震數(shù)據(jù)體；基于不同頻率的地震數(shù)據(jù)體 和地震反演初始模型獲得不同頻率的地震數(shù)據(jù)體的反演結(jié)果；基于不同頻率的地震數(shù)據(jù)體 的反演結(jié)果得到區(qū)域的定量化儲層預測結(jié)果。
[0006] 一種基于分頻數(shù)據(jù)的地震反演裝置，它包括：
[0007] 地震反演初始模型建立模塊，其用于根據(jù)區(qū)域內(nèi)的測井數(shù)據(jù)和地震數(shù)據(jù)目的層的 層位解釋建立地震反演初始模型；地震數(shù)據(jù)體提取模塊，其用于對收集區(qū)域內(nèi)的地震數(shù)據(jù) 按照一定間隔進行提取從而獲得不同頻率的地震數(shù)據(jù)體；地震數(shù)據(jù)反演模塊，其用于將不 同頻率的地震數(shù)據(jù)體的反演結(jié)果合并得到區(qū)域的反演剖面；地震反演數(shù)據(jù)目標擬合模塊， 其用于將不同頻率的地震數(shù)據(jù)體的反演結(jié)果擬合目標曲線得到區(qū)域的定量化儲層預測結(jié) 果。
[0008] 本發(fā)明實施例中基于分頻數(shù)據(jù)的地震反演方法及其裝置將不同分頻數(shù)據(jù)體通過 地震反演獲得的反演結(jié)果，用于擬合目標曲線，有效提升低頻、高頻信息在地震反演結(jié)果中 的作用，同時將原有分頻反演定性預測提高到定量化預測精度，進而得到定量化儲層預測 結(jié)果，為后期地質(zhì)分析和井位部署提供更加可靠的依據(jù)。
【附圖說明】
[0009] 在此描述的附圖僅用于解釋目的，而不意圖w任何方式來限制本發(fā)明公開的范 圍。另外，圖中的各部件的形狀和比例尺寸等僅為示意性的，用于幫助對本發(fā)明的理解，并 不是具體限定本發(fā)明各部件的形狀和比例尺寸。本領域的技術(shù)人員在本發(fā)明的教導下，可 W根據(jù)具體情況選擇各種可能的形狀和比例尺寸來實施本發(fā)明。
[0010] 圖1為本發(fā)明實施例中基于分頻數(shù)據(jù)的地震反演方法的流程圖。
[0011] 圖2為本發(fā)明實施例中基于分頻數(shù)據(jù)的地震反演裝置的結(jié)構(gòu)圖。
[0012] 圖3為本發(fā)明實施例中地震反演初始模型的示意圖。
[0013] 圖4a-4e為本發(fā)明實施例中不同頻率的地震數(shù)據(jù)振幅歸一化處理和去噪處理后 的地震剖面。
[0014] 圖5為本發(fā)明實施例中地震反演結(jié)果得到的儲層預測結(jié)果示意圖。
[0015] 圖6為基于常規(guī)地震反演方法得到的儲層預測結(jié)果示意圖。
【具體實施方式】
[0016] 結(jié)合附圖和本發(fā)明【具體實施方式】的描述，能夠更加清楚地了解本發(fā)明的細節(jié)。但 是，在此描述的本發(fā)明的【具體實施方式】，僅用于解釋本發(fā)明的目的，而不能W任何方式理解 成是對本發(fā)明的限制。在本發(fā)明的教導下，技術(shù)人員可W構(gòu)想基于本發(fā)明的任意可能的變 形，運些都應被視為屬于本發(fā)明的范圍。
[0017] 圖1為本發(fā)明實施例中基于分頻數(shù)據(jù)的地震反演方法的流程圖，如圖1所示，本 發(fā)明實施例基于分頻數(shù)據(jù)的地震反演方法包括W下步驟：
[0018]S101:收集區(qū)域內(nèi)的地震和測井數(shù)據(jù)得到收集的區(qū)域內(nèi)的地震和測井數(shù)據(jù)。
[0019] 收集區(qū)域內(nèi)的地震和測井數(shù)據(jù)得到收集的區(qū)域內(nèi)的地震和測井數(shù)據(jù)，收集的區(qū)域 內(nèi)的地震和測井數(shù)據(jù)包括區(qū)域內(nèi)的地震數(shù)據(jù)、測井數(shù)據(jù)、地震數(shù)據(jù)目的層的層位解釋。該收 集的區(qū)域內(nèi)的地震和測井數(shù)據(jù)可W通過保幅處理W增加地震數(shù)據(jù)的質(zhì)量。
[0020]S102 :基于收集的區(qū)域內(nèi)的地震和測井數(shù)據(jù)獲得地震反演初始模型和不同頻率的 地震數(shù)據(jù)體。
[0021] 基于收集的區(qū)域內(nèi)的地震和測井數(shù)據(jù)中區(qū)域內(nèi)的測井數(shù)據(jù)和地震數(shù)據(jù)目的層的 層位解釋建立地震反演初始模型。具體為首先利用目的層頂、底界的層位解釋搭建構(gòu)造模 型；然后利用標定后的測井數(shù)據(jù)對構(gòu)造模型內(nèi)部空間進行插值得到地震反演初始模型。
[0022] 分析地震數(shù)據(jù)中的頻譜信息獲得地震頻帶寬度，在有效地震頻帶寬度內(nèi)按照一定 間隔提取W獲得不同頻率的地震數(shù)據(jù)體。具體可W為分析原始地震數(shù)據(jù)的頻譜信息，用 W得到地震頻帶寬度，在有效地震頻帶范圍（fmi。，fmJ內(nèi)按照一定間隔提取不同頻率帶寬 (帶寬b可W為lHz、2Hz、…、a)的地震數(shù)據(jù)得到不同頻率的地震數(shù)據(jù)體Dfi，i= 1，2,…， k,義-二-K心二-，帶寬可W為1化、2化、…、曰，同時町1可W是二維數(shù)據(jù)體或者S維數(shù)據(jù) :煉 體。
[0023]S103:對不同頻率的地震數(shù)據(jù)體進行振幅歸一化處理獲得振幅歸一化處理后的不 同頻率的地震數(shù)據(jù)體。
[0024] 對各個分頻地震數(shù)據(jù)體Dfi可W進行振幅歸一化處理，振幅歸一化處理可W至少 包括最大振幅歸一化法和方差歸一化法其中之一。最大振幅歸一化法的具體計算過程如 下： 陽0巧]
。）
[00%] 其中，max值fi)表示一個頻率的地震數(shù)據(jù)體的最大振幅值，巧;表示對應的振幅歸 一化分頻地震數(shù)據(jù)體，Dfi表示地震數(shù)據(jù)體。
[0027] 方差歸一化法的具體計算過程如下：
[00測
巧)
[0029] 其中，std值fi)表示單個不同頻率的地震數(shù)據(jù)體中最大振幅方差值，表示對應 的振幅歸一化分頻地震數(shù)據(jù)體，Dfi表示地震數(shù)據(jù)體。
[0030] 本基于分頻數(shù)據(jù)的地震反演方法對獲取的地震數(shù)據(jù)體進行分頻處理獲得不同頻 率的地震數(shù)據(jù)體，進一步開展振幅歸一化處理，避免了不同頻率的地震數(shù)據(jù)由于振幅差異 太大從而在非線性擬合中引起的不確定性。
[0031] S104:對振幅歸一化處理后的不同頻率的地震數(shù)據(jù)體進行去噪處理獲得去噪處理 后的不同頻率的地震數(shù)據(jù)體。
[0032] 對振幅歸一化處理后的不同頻率的地震數(shù)據(jù)體分別做去噪處理W提高信噪比，該 步驟主要用于去除隨機噪音，去噪處理可W至少包括中值濾波法和均值濾波法其中之一。
[0033] 具體可W如下處理：在振幅歸一化處理后的不同頻率的地震數(shù)據(jù)體中定義一個時 窗，包含2j+l個樣本點，則在中值濾波法中第j+1個樣本點的值表示成：
[0034]
(3)
[0035] 其中：Sr""為采用中值濾波法獲得的時窗中間第j+1個樣本點的數(shù)值。
[0036] 在均值濾波法中第j+1個樣本點的值表示成：<