[发明专利]近地表各向异性参数的测量方法、装置、地震计及介质有效
申请号: | 201810937297.2 | 申请日: | 2018-08-16 |
公开(公告)号: | CN108957544B | 公开(公告)日: | 2019-07-26 |
发明(设计)人: | 陶开;唐有彩;钮凤林 | 申请(专利权)人: | 中国石油大学(北京) |
主分类号: | G01V1/30 | 分类号: | G01V1/30 |
代理公司: | 北京同立钧成知识产权代理有限公司 11205 | 代理人: | 吴会英;黄健 |
地址: | 102249*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 各向异性参数 快波 慢波 互相关函数 归一化相关系数 近地表 连续背景 测量 噪音 预设时间段 天然地震 地震源 观测点 三分量 求解 频域 采集 申请 | ||
1.一种近地表各向异性参数的测量方法,其特征在于,包括:
利用单个三分量地震计在待测观测点采集预设时间段的连续背景噪音;
获取所述连续背景噪音中的快波数据和慢波数据;
分别获取所述快波数据和慢波数据在频域的互相关函数,所述快波数据和所述慢波数据的互相关函数中包括所述快波数据的各向异性参数;
获取所述快波数据的互相关函数和所述慢波数据的互相关函数间的波形归一化相关系数模型;
求解使所述波形归一化相关系数模型中所述波形归一化相关系数最大时所述快波数据的各向异性参数,所述快波数据的各向异性参数为近地表各向异性参数。
2.根据权利要求1所述的方法,其特征在于,所述波形归一化相关系数模型为:其中,
θFP表示快波数据的振幅方位角;a表示各向异性强度,且a=Vfast/Vslow-1,Vfast为快波的传播速度,Vslow为慢波的传播速度;Cθθ(ω)表示互相关函数;
相应地,所述求解使所述波形归一化相关系数模型中所述波形归一化相关系数最大时所述快波数据的各向异性参数,所述快波数据的各向异性参数为近地表各向异性参数,具体包括:
采用网格搜索方法计算中的快波数据的各向异性参数;
其中,所述快波数据的各向异性参数包括:快波数据的振幅方位角和各向异性强度。
3.根据权利要求2所述的方法,其特征在于,所述获取所述快波数据和慢波数据在频域的互相关函数,包括:
计算快波数据和慢波数据在频域的振幅谱归一化的平均互相关函数;
根据所述快波数据的振幅方位角和所述快波数据和慢波数据的振幅谱归一化的平均互相关函数获取所述快波数据和慢波数据在频域的互相关函数。
4.根据权利要求3所述的方法,其特征在于,所述计算快波数据和慢波数据在频域的振幅谱归一化的平均互相关函数,包括:
计算所述快波数据和慢波数据对应的水平分量波形X(t)和Y(t),所述X(t)和Y(t)正交;
对所述水平分量波形X(t)和Y(t)进行归一化处理;
将所述水平分量波形X(t)和Y(t)分割为连续不重叠的时间窗;
对每个时间窗内的水平分量波形X(t)和Y(t)进行频域变换,以获得X(ω)和Y(ω);
采用第一公式计算所述X(ω)和Y(ω)对应的振幅谱归一化的平均互相关函数;
所述第一公式为:
其中,CXX(ω)、CXY(ω)和CYY(ω)分别为所述X(ω)和Y(ω)对应的振幅谱归一化的平均互相关函数;X*(ω)和Y*(ω)分别为X(ω)和Y(ω)的复共轭函数;所述<·>表示对所有时间窗内的频域值取平均,
5.根据权利要求4所述的方法,其特征在于,所述对每个时间窗内的水平分量波形X(t)和Y(t)进行频域变换,以获得X(ω)和Y(ω)之后,还包括:
对所述X(ω)和Y(ω)进行预处理。
6.根据权利要求5所述的方法,其特征在于,所述对所述X(ω)和Y(ω)进行预处理,包括:
对所述X(ω)和Y(ω)进行带通滤波处理;
对带通滤波处理后的X(ω)和Y(ω)进行移动窗均衡化处理;
对所述移动窗均衡化处理后的X(ω)和Y(ω)进行谱白化处理。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油大学(北京),未经中国石油大学(北京)许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201810937297.2/1.html,转载请声明来源钻瓜专利网。
- 上一篇:采用地震法预测砂泥岩储层地层压力的新方法
- 下一篇:多次波成像预测方法及系统