1.一種基于廣義S變換的疊前Q值反演方法,其特征在于,包括:
獲得地震波層狀模型的疊前CMP道集;
采用廣義S變換對所述疊前CMP道集逐道進(jìn)行時頻分析,得到目標(biāo)層的振幅譜;
基于目標(biāo)層的振幅譜,求取目標(biāo)層每個地震道的頻譜比斜率;
根據(jù)頻譜比斜率與偏移距的線性關(guān)系,求取零偏移距的射線平均品質(zhì)因子;以及
根據(jù)零偏移距的射線平均品質(zhì)因子,估算地震波層狀模型的層間Q值。
2.根據(jù)權(quán)利要求1所述的疊前Q值反演方法,其特征在于,所述獲得地震波層狀模型的疊前CMP道集,包括:
采用頻率域粘聲波有限差分波動方程對地震波層狀模型進(jìn)行正演,得到地震波層狀模型的炮記錄,對炮記錄進(jìn)行預(yù)處理和動校正,獲得疊前CMP道集。
3.根據(jù)權(quán)利要求1所述的疊前Q值反演方法,其特征在于,所述應(yīng)用廣義S變換對所述疊前CMP道集逐道進(jìn)行時頻分析,得到目標(biāo)層的振幅譜,具體包括:
引入時窗調(diào)節(jié)因子控制S變換的窗函數(shù),再通過S變換對所述疊前CMP道集逐道進(jìn)行時頻分析,獲得時頻分辨率可調(diào)的時頻譜,基于該時頻譜計算出目標(biāo)層的振幅譜。
4.根據(jù)權(quán)利要求1所述的疊前Q值反演方法,其特征在于,還包括:對得到的目標(biāo)層的振幅譜進(jìn)行動校拉伸校正,根據(jù)動校拉伸校正后的振幅譜,求取目標(biāo)層每個地震道的頻譜比斜率。
5.根據(jù)權(quán)利要求1或4所述的疊前Q值反演方法,其特征在于,采用以下公式求取目標(biāo)層每個地震道的頻譜比斜率K:
其中,S1(f)和S2(f)分別為t1和t2時刻的振幅譜,f為頻率,P和G分別為能量分配因子和幾何擴(kuò)散因子。
6.根據(jù)權(quán)利要求1所述的疊前Q值反演方法,其特征在于,所述頻譜比斜率與偏移距的線性關(guān)系表達(dá)為:
式中,Kx是偏移距為x時的頻譜比斜率,Q0為零偏移距的Q值,Δt0為零偏移距地震波通過參考軸到目標(biāo)軸所用的雙程旅行時差,v為速度,t0是零偏移距的雙程旅行時。
7.根據(jù)權(quán)利要求6所述的疊前Q值反演方法,其特征在于,所述求取零偏移距的射線平均品質(zhì)因子,具體包括:對零偏移距的Q值Q0進(jìn)行最小平方線性回歸處理,得到零偏移距的射線平均品質(zhì)因子Q。
8.根據(jù)權(quán)利要求1所述的疊前Q值反演方法,其特征在于,采用以下公式估算地震波層狀模型的層間Q值Qi:
Qi=[tn-tn-1]/[tn/Qn-tn-1/Qn-1]
其中,Qn表示第n層的射線平均品質(zhì)因子,Qn-1表示第n-1層的射線平均品質(zhì)因子,tn和tn-1分別表示地震波在第n層和第n-1層的時刻。
9.一種基于廣義S變換的疊前Q值反演系統(tǒng),其特征在于,包括:
CMP道集獲取模塊,用于獲得地震波層狀模型的疊前CMP道集;
振幅譜計算模塊,用于采用廣義S變換對所述疊前CMP道集逐道進(jìn)行時頻分析,得到目標(biāo)層的振幅譜;
頻譜比斜率計算模塊,用于基于目標(biāo)層的振幅譜,求取目標(biāo)層每個地震道的頻譜比斜率;
平均Q值計算模塊,用于根據(jù)頻譜比斜率與偏移距的線性關(guān)系,求取零偏移距的射線平均品質(zhì)因子;以及
層間Q值計算模塊,用于根據(jù)零偏移距的射線平均品質(zhì)因子,估算地震波層狀模型的層間Q值。
10.根據(jù)權(quán)利要求9所述的疊前Q值反演系統(tǒng),其特征在于,還包括振幅譜動校模塊,用于對得到的目標(biāo)層的振幅譜進(jìn)行動校拉伸校正。