【技术实现步骤摘要】
本专利技术属于勘探地球物理，具体而言，涉及一种地震资料保幅程度定量评价方法及装置。
技术介绍
1、现有的地震勘探资料保幅程度评价主要是定性评价，即处理人员凭借经验，通过视觉观察保幅处理前后剖面，或者将资料推进到叠前时间偏移,再对比剖面与测井曲线吻合度等方法评价保幅程度，这个方法需要处理人员有较丰富的经验，带有一定的主观性，且较耗费工时。
技术实现思路
1、为了解决保幅模块处理缺乏统一的定量评价标准的问题，本申请提供了一种定量计算地震勘探资料保幅处理程度的方法及装置，通过分析对比保幅处理前后的叠加剖面相干性，能够为评价分析保幅处理程度提供一种判定依据。
2、第一方面，本申请提供了一种地震资料保幅程度定量评价方法，包括：
3、s1，获取地震勘探资料保幅处理前后的叠加道集，控制两个叠加道集具有相同程度的振幅增益；
4、s2，根据步骤s1中获得的两个叠加道集，从中选取一个共中心点，分别取出该点在保幅处理前后相同位置上对应的叠加地震道，并将两个地震道作为输入数据；
5、s3，在地震道上设定时间段[a,b]，在该时间段内对步骤s2中获得的两个地震道f(t)、g(t)分别进行连续小波变换；
6、s4，根据步骤s3中获得的两个连续小波变换wf、wg，计算他们之间的交叉小波谱以及自功率谱；
7、s5，根据步骤s4中获得的交叉小波谱和自功率谱，计算他们之间的小波相干谱；
8、s6，根据步骤s5得到的小波相干谱，计算超过
9、s7、遍历叠加道集上的每一个共中心点，每一个点重复执行s2到s6的步骤，计算得到所有共中心点的小波相干谱以及其相干阈值面积百分比数值；
10、s8，根据步骤s7获得的所有叠加道集共中心点的小波相干阈值面积百分比数值，统计所有共中心点的百分比数值的算数平均值，所述算术平均值反映了地震记录的保幅处理程度。
11、在一些实施例中，s1，获取地震勘探资料保幅处理前后的叠加道集，控制两个叠加道集具有相同程度的振幅增益，包括：
12、获取地震勘探资料保幅处理前的叠加道集，设定道集的振幅增益值为a；获取地震勘探资料保幅处理后的叠加道集，设定道集的振幅增益值为a。
13、在一些实施例中，s3，在地震道上设定时间段[a,b]，在该时间段内对步骤s2中获得的两个地震道f(t)、g(t)分别进行连续小波变换，包括：
14、在地震道上设定一个时间段[a,b]，在该时间段内对步骤s2中获得的两个地震道f(t)、g(t)分别进行连续小波变换：
15、
16、其中，f(t)表示单个地震道输入；表示基本小波簇的复共轭，s是基本小波拉伸或压缩的尺度量，τ代表基本小波的平移量，它表示基本小波在地震道上滑动的位置；基本小波采用适用于地震数据的morlet函数，对其拉伸和平移后的表达式为：
17、ψs，τ(t)＝s-1/2*exp[-(t-τ)2/2s2]*exp[-jw0(t-τ)/s]*π-1/4
18、其中，j是虚数单位，w0代表motlet函数的中心频率。
19、在一些实施例中，s4，根据步骤s3中获得的两个连续小波变换wf、wg，计算他们之间的交叉小波谱以及自功率谱，包括：
20、交叉小波谱的表达式为：
21、
22、自功率谱分别为：
23、
24、
25、其中w*f(s，τ)、w*g(s，τ)分别代表地震道f(t)、g(t)经过连续小波变换后的复共轭。
26、在一些实施例中，s5，根据步骤s4中获得的交叉小波谱和自功率谱，计算他们之间的小波相干谱，包括：相干谱的表达式为：
27、
28、其中smooth()函数表示时间和尺度上的平滑函数。
29、在一些实施例中，s6，根据步骤s5得到的小波相干谱，计算超过设定阈值的显著相干面积百分比，包括：
30、设定一个频带范围bw，在步骤s5得到的小波相干谱中，筛选出频率在bw范围内的相干点，再设定一个阈值th，计算显著相干面积百分比，计算公式为：
31、
32、其中wc(bw)代表在相干谱中，频带bw范围内的相干点，公式表示小波相干谱中某一频带的面积占整个小波相干谱面积的比例。
33、在一些实施例中，s8，根据步骤s7获得的所有叠加道集共中心点的小波相干阈值面积百分比数值，统计所有共中心点的百分比数值的算数平均值，所述算术平均值反映了地震记录的保幅处理程度，包括：
34、根据步骤s7获得的所有叠加道集共中心点的小波相干阈值面积百分比数值，统计所有共中心点的百分比数值的算数平均值，公式为：
35、
36、其中n代表共中心点的个数；再将算数平均值大小作为输出值，该输出值反映了地震记录的保幅处理程度，数值越大，保幅处理程度越高。
37、第二方面，本申请提供了一种地震资料保幅程度定量评价装置，该装置用于：
38、获取地震勘探资料保幅处理前后的叠加道集，控制两个叠加道集具有相同程度的振幅增益；
39、根据获得的两个叠加道集，从中选取一个共中心点，分别取出该点在保幅处理前后相同位置上对应的叠加地震道，并将两个地震道作为输入数据；
40、在地震道上设定时间段[a,b]，在该时间段内对获得的两个地震道f(t)、g(t)分别进行连续小波变换；
41、根据获得的两个连续小波变换wf、wg，计算他们之间的交叉小波谱以及自功率谱；
42、根据获得的交叉小波谱和自功率谱，计算他们之间的小波相干谱；
43、根据得到的小波相干谱，计算超过设定阈值的显著相干面积百分比；
44、遍历叠加道集上的每一个共中心点，每一个点重复执行上述步骤，计算得到所有共中心点的小波相干谱以及其相干阈值面积百分比数值；
45、根据获得的所有叠加道集共中心点的小波相干阈值面积百分比数值，统计所有共中心点的百分比数值的算数平均值，所述算术平均值反映了地震记录的保幅处理程度。
46、第三方面，本申请提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述任一项所述方法的步骤。
47、第四方面，本申请提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述任一项所述方法的步骤。
48、本申请地震资料保幅程度定量评价方法及装置具有如下有益效果：
49、本申请能够定量的评价经过保幅处理后的地震资料保幅程度，形成数值评价标准，能够为评价分析保幅处理程度提供一种判定依据。
本文档来自技高网...【技术保护点】
1.一种地震资料保幅程度定量评价方法，其特征在于，包括：
2.根据权利要求1所述地震资料保幅程度定量评价方法，其特征在于，S1，获取地震勘探资料保幅处理前后的叠加道集，控制两个叠加道集具有相同程度的振幅增益，包括：
3.根据权利要求2所述地震资料保幅程度定量评价方法，其特征在于，S3，在地震道上设定时间段[a,b]，在该时间段内对步骤S2中获得的两个地震道f(t)、g(t)分别进行连续小波变换，包括：
4.根据权利要求1-3任一项所述地震资料保幅程度定量评价方法，其特征在于，S4，根据步骤S3中获得的两个连续小波变换Wf、Wg，计算他们之间的交叉小波谱以及自功率谱，包括：
5.根据权利要求1-3任一项所述地震资料保幅程度定量评价方法，其特征在于，S5，根据步骤S4中获得的交叉小波谱和自功率谱，计算他们之间的小波相干谱，包括：相干谱的表达式为：
6.根据权利要求1-3任一项所述地震资料保幅程度定量评价方法，其特征在于，S6，根据步骤S5得到的小波相干谱，计算超过设定阈值的显著相干面积百分比，包括：
7.根据权利
8.一种地震资料保幅程度定量评价装置，其特征在于，，所述装置用于：
9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现所述权利要求1-7中任一项所述方法的步骤。
10.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现所述权利要求1-7中任一项所述方法的步骤。
...【技术特征摘要】
1.一种地震资料保幅程度定量评价方法，其特征在于，包括：
2.根据权利要求1所述地震资料保幅程度定量评价方法，其特征在于，s1，获取地震勘探资料保幅处理前后的叠加道集，控制两个叠加道集具有相同程度的振幅增益，包括：
3.根据权利要求2所述地震资料保幅程度定量评价方法，其特征在于，s3，在地震道上设定时间段[a,b]，在该时间段内对步骤s2中获得的两个地震道f(t)、g(t)分别进行连续小波变换，包括：
4.根据权利要求1-3任一项所述地震资料保幅程度定量评价方法，其特征在于，s4，根据步骤s3中获得的两个连续小波变换wf、wg，计算他们之间的交叉小波谱以及自功率谱，包括：
5.根据权利要求1-3任一项所述地震资料保幅程度定量评价方法，其特征在于，s5，根据步骤s4中获得的交叉小波谱和自功率谱，计算他们之间的小波相干谱，包括：相干谱的表达式为：
...
【专利技术属性】
技术研发人员：付嘉玮，廖文婷，郭廷超，鲍伟，李红彩，胡国斌，
申请(专利权)人：中国石油化工股份有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。