地震剖面变密度分段色标显示方法及系统与流程

1.本发明涉及地震勘探领域，更具体地，涉及一种地震剖面变密度分段色标显示方法及系统。


背景技术：

2.地震波传播时，不同地层里的波速是不一样的，在收集的地震信号中，发射波的振幅反映了反射界面的反射能力，通过不同色标与地震波的强弱相联系，直观显示地层层位变化。目前的显示结果质量低，图像不清晰。因此，有必要开发一种地震剖面变密度分段色标显示方法及系统。
3.公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。


技术实现要素：

4.本发明提出了一种地震剖面变密度分段色标显示方法及系统，其能够通过分段色标显示，形象高效地绘制地震剖面图，通过地震能量补偿，增强图像对比性，使地震剖面图显示更加清晰明了，有利于地震处理解释工作的顺利开展，提高工作效率。
5.根据本发明的一方面，提出了一种地震剖面变密度分段色标显示方法。所述方法可以包括：读取地震数据，针对所述地震数据进行归一化处理，获得归一化地震数据；针对所述归一化地震数据的振幅进行能量补偿，获得补偿后振幅；根据所述补偿后振幅确定色标；根据所述归一化地震数据与所述色标，绘制地震剖面图。
6.优选地，通过以下步骤进行所述能量补偿：确定时窗，设定时窗的起始时间和时窗半径，进而计算归一化地震数据样点的均方根振幅；根据所述均方根振幅，计算剩余补偿因子；根据所述剩余补偿因子在道集内进行叠加，获得所述补偿后振幅。
7.优选地，通过公式(1)计算归一化地震数据样点的均方根振幅：
[0008][0009]
其中，amp
rms
为均方根振幅，t1为时窗初始时间，r是时窗半径，amp
i
是样点的振幅值。
[0010]
优选地，通过公式(2)计算所述剩余补偿因子：
[0011][0012]
其中，f为剩余补偿因子，ns＝采样点总数-当前采样点下标+1，δdt为采样间隔。
[0013]
优选地，根据所述补偿后振幅确定色标包括：确定振幅值与所述色标的对应关系；将每一个采样点的补偿振幅值代入所述对应关系，标定所述色标。
[0014]
根据本发明的另一方面，提出了一种地震剖面变密度分段色标显示系统，其特征
在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：读取地震数据，针对所述地震数据进行归一化处理，获得归一化地震数据；针对所述归一化地震数据的振幅进行能量补偿，获得补偿后振幅；根据所述补偿后振幅确定色标；根据所述归一化地震数据与所述色标，绘制地震剖面图。
[0015]
优选地，通过以下步骤进行所述能量补偿：确定时窗，设定时窗的起始时间和时窗半径，进而计算归一化地震数据样点的均方根振幅；根据所述均方根振幅，计算剩余补偿因子；根据所述剩余补偿因子在道集内进行叠加，获得所述补偿后振幅。
[0016]
优选地，通过公式(1)计算归一化地震数据样点的均方根振幅：
[0017][0018]
其中，amp
rms
为均方根振幅，t1为时窗初始时间，r是时窗半径，amp
i
是样点的振幅值。
[0019]
优选地，通过公式(2)计算所述剩余补偿因子：
[0020][0021]
其中，f为剩余补偿因子，ns＝采样点总数-当前采样点下标+1，δdt为采样间隔。
[0022]
优选地，根据所述补偿后振幅确定色标包括：确定振幅值与所述色标的对应关系；将每一个采样点的补偿振幅值代入所述对应关系，标定所述色标。
[0023]
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
[0024]
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
[0025]
图1示出了根据本发明的地震剖面变密度分段色标显示方法的步骤的流程图。
[0026]
图2示出了根据本发明的一个实施例的色标的示意图。
[0027]
图3示出了根据本发明的一个实施例的地震剖面图。
[0028]
图4示出了根据本发明的一个实施例的能量补偿后获得的地震剖面图。
具体实施方式
[0029]
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
[0030]
图1示出了根据本发明的地震剖面变密度分段色标显示方法的步骤的流程图。
[0031]
在该实施例中，根据本发明的地震剖面变密度分段色标显示方法可以包括：步骤101，读取地震数据，针对地震数据进行归一化处理，获得归一化地震数据；步骤102，针对归一化地震数据的振幅进行能量补偿，获得补偿后振幅；步骤103，根据补偿后振幅确定色标；步骤104，根据归一化地震数据与色标，绘制地震剖面图。
[0032]
在一个示例中，通过以下步骤进行能量补偿：确定时窗，设定时窗的起始时间和时窗半径，进而计算归一化地震数据样点的均方根振幅；根据均方根振幅，计算剩余补偿因子；根据剩余补偿因子在道集内进行叠加，获得补偿后振幅。
[0033]
在一个示例中，通过公式(1)计算归一化地震数据样点的均方根振幅：
[0034][0035]
其中，amp
rms
为均方根振幅，t1为时窗初始时间，r是时窗半径，amp
i
是样点的振幅值。
[0036]
在一个示例中，通过公式(2)计算剩余补偿因子：
[0037][0038]
其中，f为剩余补偿因子，ns＝采样点总数-当前采样点下标+1，δdt为采样间隔。
[0039]
在一个示例中，根据补偿后振幅确定色标包括：确定振幅值与色标的对应关系；将每一个采样点的补偿振幅值代入对应关系，标定色标。
[0040]
具体地，根据本发明的地震剖面变密度分段色标显示方法可以包括：
[0041]
读取地震数据，当前窗口显示的道数是tracenum，地震采样点数是sample，将地震数据归一化，从-1到1，存储在数组seis_data[tracenum*sample]中，获得归一化地震数据。
[0042]
确定时窗，设定时窗的起始时间和时窗半径，进而通过公式(1)计算归一化地震数据样点的均方根振幅；由于地震信号在地层传播中会随着深度增加而衰减，所以剩余补偿因子f采用线性值，根据均方根振幅，通过公式(2)计算剩余补偿因子；根据剩余补偿因子在cmp道集内进行叠加，获得补偿后振幅。
[0043]
色标colortab有n种颜色，确定振幅值与色标的对应关系，min(n)对应补偿后振幅最小值，max(n)对应补偿后振幅最大值；将每一个采样点的补偿振幅值代入对应关系，标定色标；将归一化地震数据内的所有补偿振幅值根据色标确定颜色，绘制地震剖面图。
[0044]
本方法通过分段色标显示，形象高效地绘制地震剖面图，通过地震能量补偿，增强图像对比性，使地震剖面图显示更加清晰明了，有利于地震处理解释工作的顺利开展，提高工作效率。
[0045]
应用示例
[0046]
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
[0047]
根据本发明的地震剖面变密度分段色标显示方法包括：
[0048]
读取地震数据，当前窗口显示的道数是tracenum，地震采样点数是sample，将地震数据归一化，从-1到1，存储在数组seis_data[tracenum*sample]中，获得归一化地震数据。
[0049]
确定时窗，设定时窗的起始时间和时窗半径，进而通过公式(1)计算归一化地震数据样点的均方根振幅；由于地震信号在地层传播中会随着深度增加而衰减，所以剩余补偿因子f采用线性值，根据均方根振幅，通过公式(2)计算剩余补偿因子；根据剩余补偿因子在cmp道集内进行叠加，获得补偿后振幅。
[0050]
图2示出了根据本发明的一个实施例的色标的示意图。
[0051]
色标colortab有n种颜色，屏幕像素数组为subbit[nwidth*nheight*4]，用rgb色彩组合表示颜色，包括颜色值r、颜色值g、颜色值b，则色标如图2所示。确定振幅值与色标的对应关系，min(n)对应补偿后振幅最小值，max(n)对应补偿后振幅最大值；将每一个采样点的补偿振幅值代入对应关系，标定色标；将归一化地震数据内的所有补偿振幅值根据色标确定颜色，绘制地震剖面图。
[0052]
图3示出了根据本发明的一个实施例的未经能量补偿获得的地震剖面图。
[0053]
图4示出了根据本发明的一个实施例的能量补偿后获得的地震剖面图。
[0054]
二者对比可知，经过能量补偿后获得的地震剖面图更加清晰，显示的数据更多。
[0055]
综上所述，本发明通过分段色标显示，形象高效地绘制地震剖面图，通过地震能量补偿，增强图像对比性，使地震剖面图显示更加清晰明了，有利于地震处理解释工作的顺利开展，提高工作效率。
[0056]
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
[0057]
根据本发明的实施例，提供了一种地震剖面变密度分段色标显示系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：读取地震数据，针对地震数据进行归一化处理，获得归一化地震数据；针对归一化地震数据的振幅进行能量补偿，获得补偿后振幅；根据补偿后振幅确定色标；根据归一化地震数据与色标，绘制地震剖面图。
[0058]
在一个示例中，通过以下步骤进行能量补偿：确定时窗，设定时窗的起始时间和时窗半径，进而计算归一化地震数据样点的均方根振幅；根据均方根振幅，计算剩余补偿因子；根据剩余补偿因子在道集内进行叠加，获得补偿后振幅。
[0059]
在一个示例中，通过公式(1)计算归一化地震数据样点的均方根振幅：
[0060][0061]
其中，amp
rms
为均方根振幅，t1为时窗初始时间，r是时窗半径，amp
i
是样点的振幅值。
[0062]
在一个示例中，通过公式(2)计算剩余补偿因子：
[0063][0064]
其中，f为剩余补偿因子，ns＝采样点总数-当前采样点下标+1，δdt为采样间隔。
[0065]
在一个示例中，根据补偿后振幅确定色标包括：确定振幅值与色标的对应关系；将每一个采样点的补偿振幅值代入对应关系，标定色标。
[0066]
本系统通过分段色标显示，形象高效地绘制地震剖面图，通过地震能量补偿，增强
图像对比性，使地震剖面图显示更加清晰明了，有利于地震处理解释工作的顺利开展，提高工作效率。
[0067]
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
[0068]
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。