專利名稱：頻散avo反演方法
技術領域：
本發明涉及地震勘探技術領域，是一種新的地震屬性，用于描述地震信號的頻散程度，具體講是一種頻散AVO (Amplitude Versus Offset)反演方法。
背景技術：
目前的縱波AVO反演技術僅僅考慮了地震波的振幅信息隨偏移距的變化情況，認為儲層中流體性質和流體飽和度的變化導致地層中波阻抗的差異，從而引起地震波反射振幅的變化，而這種變化是與頻率無關的。然而無論砂巖儲層還是碳酸鹽巖儲層，都是由固體和液體兩種部分組成的雙相介質，流體的存在導致地震波發生不同程度的頻散和衰減。目前的AVO分析忽略了振幅隨頻率的變化特征，而這一特征受儲層中所含流體性質的影響， 是描述儲層中流體性質的依據之一。Zo印pritz方程是AVO技術的理論基礎。當彈性波傳播到彈性分界面時，要發生波型轉換和能量重新分配，在不同介質分界面上，根據位移和應力的邊界條件，通過解波動方程，可以得到彈性分界面能量分配關系。Knott和heppritz利用反射界面兩側位移和應力的連續性作為邊界條件，得到了反射系數和透射系數，它們是入射角和介質彈性參數(密度、體積模量和剪切模量)的函數，即Zc^ppritz方程
權利要求
1.一種頻散AVO反演方法，其特征在于，所述頻散AVO反演方法包括 獲取疊前道集，并對所述疊前道集進行球面擴散補償和動校正；對所述進行球面擴散補償和動校正的疊前道集進行時頻分析，生成不同頻率下的振幅譜，所述振幅譜為采樣時間、接收道號和頻率的函數； 對所述不同頻率下的振幅譜進行振幅譜均衡處理；獲取所述疊前道集中的采樣點的疊加速度，并根據所述疊加速度生成所述采樣點的層速度；根據所述譜均衡處理后的振幅譜和所述采樣點的層速度，采用最小二乘反演，生成設定頻率下的縱橫波速度變化率和所述采樣點的頻散程度，所述頻散程度為縱橫波速度變化率隨頻率變化的快慢；將所述生成的頻散程度按照道集序號進行排列，生成頻散程度屬性剖面。
2.如權利要求1所述的頻散AVO反演方法，其特征在于，所述對進行球面擴散補償和動校正的疊前道集進行時頻分析，生成不同頻率下的振幅譜，包括所述疊前道集有η個接收道，表示為矩陣形式為s (t，η)，采用短時傅里葉變換或者小波變換的時頻分析方法，進行頻譜分解得到不同頻率f下的振幅譜S (t，η, f)。
3.如權利要求1所述的頻散AVO反演方法，其特征在于，所述對不同頻率下的振幅譜進行振幅譜均衡處理，包括對不同頻率的振幅譜通過加權函數ω進行譜均衡Sb(t,n,f) = S(t，n，f) ω (f)，其中，S (t，n，f)為不同頻率下的振幅譜，&(t，n，f)為進行譜均衡處理后的振幅譜，ω (f)為進行譜均衡處理的加權函數。
4.如權利要求1所述的頻散AVO反演方法，其特征在于，所述根據所述疊加速度生成所述采樣點的層速度，包括獲取所述疊前道集中的采樣點的疊加速度，根據Dix公式
5.如權利要求3所述的頻散AVO反演方法，其特征在于，所述根據所述譜均衡處理后的振幅譜和所述采樣點的層速度，根據最小二乘反演，生成設定頻率下的縱橫波速度變化率和所述采樣點的頻散程度，包括根據公式下述兩個公式獲取所述縱橫波速度變化率和頻散程度
6.如權利要求1所述的頻散AVO反演方法，其特征在于，所述疊前道集包括共反射點道集CRP、共中心點道集CMP和共深度點道集⑶P。
全文摘要
本發明公開了一種頻散AVO反演方法，包括獲取疊前道集，并對疊前道集進行球面擴散補償和動校正；對疊前道集進行瞬時頻譜分析，生成不同頻率下的振幅譜；對振幅譜進行振幅譜均衡處理；獲取疊前道集中的采樣點的疊加速度，并根據疊加速度生成采樣點的層速度；根據譜均衡處理后的振幅譜和采樣點的層速度，根據最小二乘反演，生成設定頻率下的縱橫波速度變化率和所述采樣點的頻散程度；將生成的頻散程度按照道集序號進行排列，生成頻散程度屬性剖面。本發明實施例的頻散AVO反演方法，消除了由于彈性界面波阻抗差異導致的強震幅能量團的影響，也可避免疊后地震資料頻譜分析時疊加處理可能導致的虛假“頻率異常”，使得到的結果更加可靠。
文檔編號G01V1/28GK102478667SQ201010567538
公開日2012年5月30日 申請日期2010年11月30日 優先權日2010年11月30日
發明者吳小羊, 張峰, 李向陽, 陳雙全 申請人:中國石油大學(北京), 中國石油天然氣集團公司