[发明专利]地震反演剥离煤层强反射的方法和装置

权利要求书

1.一种地震反演剥离煤层强反射的方法，该方法包括：

提取原始地震偏移成果数据体S(t)中每一道的统计子波W_n(f_m)以及部分频谱其中，原始地震偏移成果数据体S(t)包括N个道，下标n表示道号，n＝1,2,...,N，共M个频点，f_m＝(m-1)×Δf+f₁，m＝1,2,...,M，[f₁,f_M]是预设的频谱区间；

针对每一道，基于使目标函数J_a极小反演得到顶界面等效反射系数r_n,l,u、底界面等效反射系数r_n,l,d以及r_n,l,u和r_n,l,d分别对应的时间位置和

J_a＝‖Gα-D‖²，

其中，下标l表示煤层序号，l＝1,...,L，共L个煤层，G为M×2L的矩阵，G的第m行第2l-1列的元素第m行第2l列的元素α＝[r_n,1,u,r_n,1,d,...,r_n,L,u,r_n,L,d]^T，

基于r_n,l,u、r_n,l,d、得到煤层强反射数据体S^*(t)；

基于下式得到储层弱反射数据体S_w(t)：

S_w(t)＝S(t)-S^*(t)。

2.根据权利要求1所述的方法，其中，提取每一道的W_n(f_m)和包括：

结合地质的分层、测井曲线以及原始地震偏移成果数据体S(t)，进行井震标定，以提取每一道的W_n(f_m)和

3.根据权利要求1所述的方法，其中，采用非常快速模拟退火算法求解目标函数J_a以得到r_n,l,u、r_n,l,d、其中，在求解过程中，设置初始的时间位置分别是校正前相应的顶、底界面时间，并且设置第1道的初始的顶、底界面等效反射系数均为0，以及除第1道外每一道的煤层l的初始的顶、底界面等效反射系数为上一道求解得到的r_n-1,l,u和r_n-1,l,d。

4.根据权利要求3所述的方法，其中，针对每一道，时间位置与校正前相应的顶界面时间的差值小于5ms，时间位置与校正前的相应的底界面时间的差值小于5ms。

5.根据权利要求1所述的方法，其中，得到煤层强反射数据体S^*(t)包括：

针对每一道，基于下式得到该道的强地震响应

其中，符号“*”表示褶积运算，W_n(t)是W_n(f_m)，m＝1,2,...,M对应的时域信号，R_n(t)是基于r_n,l,u、r_n,l,d、l＝1,...,L得到的该道的反射系数；

通过合成所述N个道的强地震响应得到煤层强反射数据体S^*(t)。