【技术实现步骤摘要】
本申请属于地球物理勘探，具体涉及一种基于阈值约束的构造信息提取方法及装置。
技术介绍
1、目前常规的深度域层析速度建模方法通常不考虑构造信息，反演速度模型缺乏地质意义，反演精度较低。针对该问题，如何在层析速度建模过程中加入地质构造信息进行约束是解决该问题的关键。而在这个技术流程过程中，如何高效、高质量的提取构造信息是首先需要解决的问题。目前提取构造的方法是以斜率或者局部能量叠加的方式获得构造信息，该种方法通常计算较为繁琐，且计算结果不够稳定，对于提取的构造信息无法通过某种参数定量控制，需要额外采用滤波等技术进行修饰和修改，在实际应用过程中，考虑到处理周期和复杂性等问题，应用较少。
技术实现思路
1、基于以上技术问题，本申请提出一种基于阈值约束的构造信息提取方法及装置。
2、第一方面，本申请提出一种基于阈值约束的构造信息提取方法，包括：
3、获取三维地震数据；
4、取所述三维地震数据的一个点，计算所述点在第一方向上的第一加权差分算子以及在第二方向上的第二加权差分算子；
5、根据所述第一加权差分算子以及第二加权差分算子，得到初始构造信息；
6、根据所述初始构造信息计算得到阈值；
7、基于所述阈值，对所述初始构造信息进行归一化处理，得到最终的构造信息。
8、所述三维地震数据的x轴方向为地震主测线方向的数据，y轴方向为地震联络测线方向的数据，z轴方向为地震深度方向的数据所述三维地震数据的一个点表示为w(i,
9、所述在x轴方向上与点w(i,j,k)邻近的n个点为：w(i-1,j,k)…
10、以及
11、所述在y轴方向上与点w(i,j,k)邻近的n个点为：w(i,j-1,k)…
12、以及
13、所述第一加权差分算子，计算式如下：
14、
15、其中，cx(i,j,k)为第一加权差分算子，q1为第一设置权重，q2为第设置权重，q3为第设置权重，w(i-1,j,k)为x轴负方向上与点w(i,j,k)邻近的第一个点，为x轴负方向上与点w(i,j,k)邻近的第个点，为x轴负方向上与点w(i,j,k)邻近的第个点，w(i+1,j,k)为x轴正方向上与点w(i,j,k)邻近的第一个点、为x轴正方向上与点w(i,j,k)邻近的第个点，为x轴正方向上与点w(i,j,k)邻近的第个点。
16、所述第二加权差分算子，计算式如下：
17、
18、其中，cy(i,j,k)为第二加权差分算子，q1为第一设置权重，q2为第设置权重，q3为第设置权重，w(i,j-1,k)为在y轴负方向上与点w(i,j,k)邻近的第一个点，为在y轴负方向上与点w(i,j,k)邻近的第个点，为在y轴负方向上与点w(i,j,k)邻近的第个点，w(i,j+1,k)为在y轴正方向上与点w(i,j,k)邻近的第一个点，为在y轴正方向上与点w(i,j,k)邻近的第个点，为在y轴正方向上与点w(i,j,k)邻近的第个点。
19、所述根据所述第一加权差分算子以及第二加权差分算子，得到初始构造信息，计算式如下：
20、
21、其中，s(i,j,k)为初始构造信息，cx(i,j,k)为第一加权差分算子，cy(i,j,k)为第二加权差分算子。
22、所述阈值，计算式如下：
23、
24、其中，thre为阈值，s(i,j,k)为初始构造信息，ε为设定的控制系数，n1为三维地震数据x轴的数据个数，n2三维地震数据y轴的数据个数，n3为三维地震数据z轴的数据个数。
25、所述将所述初始构造信息通过阈值进行归一化处理，得到最终的构造信息，计算式如下：
26、
27、其中，s(i,j,k)为初始构造信息，thre为阈值，为最终的构造信息。
28、第二方面，本申请提出一种基于阈值约束的构造信息提取装置，包括：
29、数据获取模块，用于获取三维地震数据；
30、差分算子计算模块，用于取所述三维地震数据的一个点，计算所述点在第一方向上的第一加权差分算子以及在第二方向上的第二加权差分算子；
31、初始构造信息计算模块，用于根据所述第一加权差分算子以及第二加权差分算子，得到初始构造信息；
32、阈值计算模块，用于根据所述初始构造信息计算得到阈值；
33、最终构造信息输出模块，用于将所述初始构造信息通过阈值进行归一化处理，得到最终的构造信息。
34、第三方面，本申请提出一种电子设备，包括：存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，执行所述的基于阈值约束的构造信息提取方法。
35、第四方面，本申请提出一种计算机可读存储介质，其存储有可执行指令，所述指令当被执行时使得处理器执行所述的基于阈值约束的构造信息提取方法。
36、有益效果：
37、本申请提出一种基于阈值约束的构造信息提取方法及装置，利用加权差分算子计算得到初始构造信息，然后基于初始构造信息，计算得到阈值，通过阈值约束得到最终的构造信息，本申请的方法剔除了噪声或对层析反演贡献较小的小尺度构造等的影响，仅留下对层析速度建模有贡献的构造信息，实现了面向深度域层析速度建模的构造信息提取。
本文档来自技高网...【技术保护点】
1.一种基于阈值约束的构造信息提取方法，其特征在于，包括：
2.根据权利要求1所述的基于阈值约束的构造信息提取方法，其特征在于，所述三维地震数据的x轴方向为地震主测线方向的数据，y轴方向为地震联络测线方向的数据，z轴方向为地震深度方向的数据，所述三维地震数据的一个点表示为W(i,j,k)，则所述在第一方向上的第一加权差分算子为在x轴方向上与点W(i,j,k)邻近的N个点加权差分的结果；所述在第二方向上的第二加权差分算子为在y轴方向上与点W(i,j,k)邻近的N个点加权差分的结果，其中，N为2的整数倍。
3.根据权利要求2所述的基于阈值约束的构造信息提取方法，其特征在于，所述在x轴方向上与点W(i,j,k)邻近的N个点为：以及
4.根据权利要求1所述的基于阈值约束的构造信息提取方法，其特征在于，所述第一加权差分算子，计算式如下：
5.根据权利要求1所述的基于阈值约束的构造信息提取方法，其特征在于，所述第二加权差分算子，计算式如下：
6.根据权利要求1所述的基于阈值约束的构造信息提取方法，其特征在于，所述根据所述第一加权差分
7.根据权利要求1所述的基于阈值约束的构造信息提取方法，其特征在于，所述阈值，计算式如下：
8.根据权利要求1所述的基于阈值约束的构造信息提取方法，其特征在于，所述将所述初始构造信息通过阈值进行归一化处理，得到最终的构造信息，计算式如下：
9.一种基于阈值约束的构造信息提取装置，其特征在于，包括：
10.一种电子设备，其特征在于，包括：存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，执行权利要求1-8任一项所述的基于阈值约束的构造信息提取方法。
11.一种计算机可读存储介质，其特征在于，其存储有可执行指令，所述指令当被执行时使得处理器执行权利要求1-8任一项所述的基于阈值约束的构造信息提取方法。
...【技术特征摘要】
1.一种基于阈值约束的构造信息提取方法，其特征在于，包括：
2.根据权利要求1所述的基于阈值约束的构造信息提取方法，其特征在于，所述三维地震数据的x轴方向为地震主测线方向的数据，y轴方向为地震联络测线方向的数据，z轴方向为地震深度方向的数据，所述三维地震数据的一个点表示为w(i,j,k)，则所述在第一方向上的第一加权差分算子为在x轴方向上与点w(i,j,k)邻近的n个点加权差分的结果；所述在第二方向上的第二加权差分算子为在y轴方向上与点w(i,j,k)邻近的n个点加权差分的结果，其中，n为2的整数倍。
3.根据权利要求2所述的基于阈值约束的构造信息提取方法，其特征在于，所述在x轴方向上与点w(i,j,k)邻近的n个点为：以及
4.根据权利要求1所述的基于阈值约束的构造信息提取方法，其特征在于，所述第一加权差分算子，计算式如下：
5.根据权利要求1所述的基于阈值约束的构造信息提取方法，其特征在于，所述第二加权差分算...
【专利技术属性】
技术研发人员：郑浩，陶永慧，刘立民，
申请(专利权)人：中国石油化工股份有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。