[发明专利]一种时频域真振幅恢复的方法和纵向能量补偿的方法有效
| 申请号： | 201711062296.X | 申请日： | 2017-11-02 |
| 公开（公告）号： | CN108303739A | 公开（公告）日： | 2018-07-20 |
| 发明（设计）人： | 曹中林;何光明;张新华;李忠;陈爱萍;陈丹;杨海涛;王珑;熊晶璇 | 申请（专利权）人： | 中国石油集团东方地球物理勘探有限责任公司 |
| 主分类号： | G01V1/32 | 分类号： | G01V1/32;G01V1/28 |
| 代理公司： | 成都中玺知识产权代理有限公司 51233 | 代理人： | 谭昌驰;张敏 |
| 地址： | 071000 河北*** | 国省代码： | 河北;13 |
| 权利要求书： | 查看更多 | 说明书： | 查看更多 |
| 摘要： | |||
| 搜索关键词： | 振幅恢复 地震数据 时频域 衰减 能量补偿 振幅曲线 加窗 分数阶傅里叶变换 相乘 傅里叶反变换 高阶多项式 比例因子 地层吸收 地震资料 二维数据 衰减函数 分数阶 时间点 拟合 频点 频域 吸收 恢复 | ||
1.一种时频域真振幅恢复的方法，其特征在于，所述方法包括以下步骤：
通过加窗分数阶傅里叶变换将地震数据S(x,t)中每一道变换到时频域，在所述时频域，所述每一道的数据都变成包括频率u和时间t的二维数据，其中，x为道数；
在时频域中，求取每个所述二维数据中的每个频点在不同时间点对应的振幅值，获得每个所述二维数据中的每个频点的第一时间振幅曲线Amp(t,u)；
采用高阶多项式拟合所述第一时间振幅曲线Amp(t,u)，得到吸收衰减函数du，所述吸收衰减函数du的各项系数通过非线性迭代方法求解；
对所述吸收衰减函数du进行变换，得到每个所述二维数据中的每个频点在不同时间点的衰减比例因子St，St＝exp(-du)；
将所述衰减比例因子St与所述第一时间振幅曲线Amp(t,u)相乘，得到每个所述二维数据中的每个频点的第二时间振幅曲线Amp'(t,u)；
对每个所述二维数据中的每个频点的第二时间振幅曲线Amp'(t,u)进行加窗分数阶傅里叶反变换，得到消除了吸收衰减之后的地震数据S'(x,t)。
2.根据权利要求1所述的时频域真振幅恢复的方法，其特征在于，所述二维数据为zp(t,u)，其中，
其中，Kp(t,u)为分数阶傅里叶变换的变换核，s(t)为地震数据S(x,t)中的一道，g(u-t)为高斯函数g(u)的积分表达式，p为分数阶傅里叶变换的阶数。
3.根据权利要求2所述的时频域真振幅恢复的方法，其特征在于，所述分数阶傅里叶变换的变换核Kp(t,u)为：
其中，p为分数阶傅里叶变换的阶数，α＝pπ/2，α为旋转角度，j为虚部。
4.根据权利要求2所述的时频域真振幅恢复的方法，其特征在于，所述高斯函数g(u)包括均值为零的高斯函数，g(u)为：
其中，σt为时域中标准偏差。
5.根据权利要求2所述的时频域真振幅恢复的方法，其特征在于，所述第一时间振幅曲线Amp(t,u)为：
Amp(t,u)＝abs(zp(t,u))，
其中，abs为绝对值的标识。
6.根据权利要求2所述的时频域真振幅恢复的方法，其特征在于，所述非线性迭代方法包括模拟退火算法。
7.一种地震数据中纵向能量的补偿方法，其特征在于，所述补偿方法包括采用权利要求1至6中任意一项所述的时频域真振幅恢复的方法进行能量补偿。
该专利技术资料仅供研究查看技术是否侵权等信息，商用须获得专利权人授权。该专利全部权利属于中国石油集团东方地球物理勘探有限责任公司，未经中国石油集团东方地球物理勘探有限责任公司许可，擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作，请联系【客服】
本文链接：http://www.vipzhuanli.com/pat/books/201711062296.X/1.html，转载请声明来源钻瓜专利网。
