[发明专利]一种获取地震波传播方向的方法及装置有效
申请号: | 201510661856.8 | 申请日: | 2015-10-14 |
公开(公告)号: | CN105372704B | 公开(公告)日: | 2018-06-01 |
发明(设计)人: | 戴南浔;章威;吴晓丰;武威;刘磊;李文卉 | 申请(专利权)人: | 中国石油天然气集团公司;中国石油集团东方地球物理勘探有限责任公司 |
主分类号: | G01V1/30 | 分类号: | G01V1/30 |
代理公司: | 北京三友知识产权代理有限公司 11127 | 代理人: | 党晓林 |
地址: | 100007 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 地震波传播 地震波 处理效率 传播方向 地震数据 构造结构 波场 单频 申请 | ||
1.一种获取地震波传播方向的方法,其特征在于,包括:
利用所获取的地震波的单频波场,构造结构张量,所述结构张量中的每个元素均是三维数据体;
计算所述结构张量的特征值;
利用所述结构张量的特征值,确定所述地震波的传播方向,
其中,计算所述结构张量的特征值包括:
对所述结构张量进行平滑;
计算平滑后的结构张量的特征值;
其中,对所述结构张量进行平滑包括:
将所述结构张量的一个三维数据体中的一坐标点作为当前坐标点;
以所述当前坐标点为中心,计算所述当前坐标点在第二预设范围内的平均值;
将所得到的平均值作为所述当前坐标点的数值;
按照上述计算一坐标点的数值的方法,依次计算所述三维数据体中剩余N-1个坐标点的数值,N为一个三维数据体中坐标点的个数;
按照上述计算一个三维数据体中所有坐标点的数值的方法,依次计算所述结构张量中剩余M-1个三维数据体中所有坐标点的数值,M为三维数据体的个数。
2.根据权利要求1所述的方法,其特征在于,所述单频波场通过以下方式获取:
获取目标区域中第一预设范围内的地震波场以及地震波速度数据;
利用所获取的地震波场以及所述地震波速度数据,获取预设频率所对应的单频波场。
3.根据权利要求2所述的方法,其特征在于,所述获取预设频率所对应的单频波场包括:
利用单程波延拓方法或双程波延拓方法,计算所述第一预设范围内地下不同深度处的单频波场。
4.根据权利要求1所述的方法,其特征在于,所述利用所述结构张量的特征值,确定所述地震波的传播方向包括:
从所述结构张量的多个特征值中选取绝对值最大的特征值;
将所选取的特征值所对应的特征向量确定为所述地震波波前的法线方向。
5.一种获取地震波传播方向的装置,其特征在于,包括:
构造单元,用于利用所获取的所述地震波的单频波场,构造结构张量,所述结构张量中的每个元素均是三维数据体;
计算单元,用于计算所述结构张量的特征值;
确定单元,用于利用所述结构张量的特征值,确定所述地震波的传播方向,
其中,所述计算单元包括:
平滑子单元,用于对所述结构张量进行平滑;
计算子单元,用于计算平滑后的结构张量的特征值;
其中,所述平滑子单元具体用于:
将所述结构张量的一个三维数据体中的一坐标点作为当前坐标点;
以所述当前坐标点为中心,计算所述当前坐标点在第二预设范围内的平均值;
将所得到的平均值作为所述当前坐标点的数值;
按照上述计算一坐标点的数值的方法,依次计算所述三维数据体中剩余N-1个坐标点的数值;
按照上述计算一个三维数据体中所有坐标点的数值的方法,依次计算所述结构张量中剩余M-1个三维数据体中所有坐标点的数值。
6.根据权利要求5所述的装置,其特征在于,所述确定单元包括:
选取子单元,用于从所述结构张量的多个特征值中选取绝对值最大的特征值;
确定子单元,用于将所选取的特征值所对应的特征向量确定为所述地震波波前的法线方向。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油天然气集团公司;中国石油集团东方地球物理勘探有限责任公司,未经中国石油天然气集团公司;中国石油集团东方地球物理勘探有限责任公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201510661856.8/1.html,转载请声明来源钻瓜专利网。