以地質(zhì)信息、測井資料為約束條件,構(gòu)建寬頻帶初始波阻抗模型,研究參數(shù)置換法建立擬線性矩陣方程,對反演解加入硬約束條件,采用非線性預(yù)條件共軛梯度數(shù)值算法,改進(jìn)搜索方向,降低系數(shù)矩陣條件數(shù),提高收斂速度,實現(xiàn)從帶限的地震數(shù)據(jù)資料中恢復(fù)出寬頻帶波阻抗的空間分布特征。通過砂泥巖互層復(fù)雜速度模型數(shù)值反演,在5%的隨機(jī)噪聲干擾下,反演得到速度模型也能夠描述出儲層和逆斷層的位置,由于噪音強(qiáng)度分布不同,造成各道反演得到的絕對速度值存在偏差,但最大偏差值僅為真實值的±5%。進(jìn)一步驗證算法可行有效,反演精度高,具有一定的抗噪能力。

【文章頁數(shù)】：7 頁

【部分圖文】：

圖1寬帶約束波阻抗反演流程圖

(3)建立反演目標(biāo)函數(shù)，采用參數(shù)置換法對原始波阻抗進(jìn)行參數(shù)變換，得到式(6)近似線性系統(tǒng);采用更改的預(yù)處理共軛梯度法求解，對初始波阻抗模型進(jìn)行迭代修改，直到求解的模型波阻抗合成的地震記錄與原始地震數(shù)據(jù)最佳匹配為止，進(jìn)而獲得最優(yōu)的高分辨率波阻抗模型，寬帶約束波阻抗反演流程見圖1。4....

圖2十二層單道理論速度模型及其初始模型

設(shè)計一個12層單道理論速度模型及一個6層初始速度模型，相關(guān)系數(shù)僅為0.785(圖2)，采樣時間為2ms，密度常數(shù)1.0，采樣點數(shù)235個。選取地震子波為零相位阻尼余弦子波，采樣率2ms，子波的長度為40ms，主頻為60Hz，地震記錄是由地震子波與反射系數(shù)相褶積形成，加入....

圖3地震記錄中加入5%隨機(jī)噪聲反演結(jié)果對比

利用fortran語言編制計算機(jī)程序進(jìn)行波阻抗反演，其反演效果如圖3所示:由圖3可知，在原始地震記錄中加入5%的隨機(jī)噪聲，反演結(jié)果與理論模型相關(guān)系數(shù)可達(dá)0.997，剖面與絕對波阻抗值均與理論值大體相同;對比采用稀疏約束反褶積法[13-14]，由迭代求出的反射系數(shù)存在誤差，小的反射....

圖4砂泥巖互層巖性油氣藏模型

以一個61道砂泥巖互層理論速度模型模擬實際油氣儲層。該模型大致可分為5層:地表為速度較低的第四系覆蓋層;第二層為泥巖層;第三層為楔狀砂巖體儲層，局部含有油氣資源;而第四層泥巖層局部存在逆斷層，使得下部第五層含水砂巖向上侵入，各層的巖性及速度值見表1。選取的地震子波為阻尼余弦子波，....

本文編號：3922988