CN108181649B - 一种介质分界面振幅补偿方法及装置 - Google Patents
一种介质分界面振幅补偿方法及装置 Download PDFInfo
- Publication number
- CN108181649B CN108181649B CN201711181682.0A CN201711181682A CN108181649B CN 108181649 B CN108181649 B CN 108181649B CN 201711181682 A CN201711181682 A CN 201711181682A CN 108181649 B CN108181649 B CN 108181649B
- Authority
- CN
- China
- Prior art keywords
- amplitude
- interface
- dielectric interface
- value
- seismic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 95
- 230000005540 biological transmission Effects 0.000 claims abstract description 79
- 238000012545 processing Methods 0.000 claims abstract description 50
- 238000003860 storage Methods 0.000 claims description 14
- 238000005070 sampling Methods 0.000 claims description 8
- 238000004454 trace mineral analysis Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 abstract description 21
- 238000004364 calculation method Methods 0.000 abstract description 13
- 238000010586 diagram Methods 0.000 description 21
- 230000006870 function Effects 0.000 description 13
- 238000004458 analytical method Methods 0.000 description 12
- 239000003245 coal Substances 0.000 description 10
- 230000006872 improvement Effects 0.000 description 8
- 238000004590 computer program Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 239000000047 product Substances 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000003780 insertion Methods 0.000 description 2
- 230000037431 insertion Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000011017 operating method Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请提供一种介质分界面振幅补偿方法及装置。所述方法包括：根据待分析地震数据，获取地震道的均方根振幅值；拾取介质分界面，根据地震道的均方根振幅值，计算介质分界面对应的分界面反射振幅值；根据分界面反射振幅值和预设的地震能量总值，计算介质分界面对应的分界面透射振幅值；根据分界面透射振幅值，计算分界面振幅补偿因子；根据分界面振幅补偿因子和地震道的均方根振幅值，进行介质分界面的振幅补偿操作。利用本申请中各个实施例，实现了介质分界面处振幅的补偿，使得振幅补偿后的地震数据更加真实的反映地下地质的变化。为后续地震资料的处理和解释提供了准确的数据基础，提高了地震数据处理结果的准确性。
Description
技术领域
本申请属于石油勘探技术领域，尤其涉及一种介质分界面振幅补偿方法及装置。
背景技术
随着科技的发展，通过测井技术和地震资料的分析处理，进行油气的勘探在石油勘探领域内应用的越来越多。通常可以采用人工方法激发地震波，在各个检波点处接收地震数据。通过对地震数据的处理，形成能够真实反映地下反射体的图像，通过图像分析处理，识别地下构造的形态和地层的波阻抗特征，判断是否储存有油气。因此，需要获得能够准确反映地质结构的地震数据。
由于地质结构比较复杂，地震波会因地质结构的不同导致地震能量的变化。尤其在介质分界面处，当激发的地震波在地下传播到介质分界面时，由于波阻抗的变化引起质点振动位移振幅的变化，能量在分界处重新分配。部分能量透射，部分能量反射，还有部分能量可能进行了转化。介质分界面处的地震能量发生了变化，振幅发生了畸变。接收到的地震数据不能准确的表示地质的真实含义，给地震解释和储层反演都造成很大困难。现有技术中，对于介质分界面处振幅畸变的处理，通常针对地表因素、接收和激发因素进行，如球面扩散补偿、地表一致性振幅补偿等。但是，没有考虑地震波在介质分界面处的能量分配问题，可能导致介质分界面处振幅补偿处理结果的准确性不高，影响地震数据的处理精度。因此，业内亟需一种能够进行介质分界面处振幅的补偿处理，进一步提高地震数据的处理精度的实施方案。
发明内容
本申请目的在于提供一种介质分界面振幅补偿方法及装置，基于介质分界面处地震能量的重新分配，进行介质分界面处振幅的补偿，可以沿介质分界面进行振幅的补偿。为后续地震资料的处理和解释提供了准确的数据基础，提高了地震数据处理结果的准确性。
一方面本申请提供了一种介质分界面振幅补偿方法，包括：
根据待分析地震数据，获取地震道的均方根振幅值；
拾取介质分界面，根据所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值；
根据所述分界面反射振幅值和预设的地震能量总值，计算所述介质分界面对应的分界面透射振幅值；
根据所述分界面透射振幅值，计算分界面振幅补偿因子；
根据所述分界面振幅补偿因子和所述地震道的均方根振幅值，进行所述介质分界面的振幅补偿操作。
进一步地，所述方法的另一个实施例中，所述根据待分析地震数据，获取地震道的均方根振幅值，包括：
根据所述待分析地震数据，利用计算所述地震道的均方根振幅值；
上式中，Aj表示第j个地震道的均方根振幅，j表示地震道序号，L表示第j个地震道分析时窗内的总样点数，Xi表示第j个地震道分析时窗内第i个样点的振幅值。
进一步地，所述方法的另一个实施例中，所述拾取介质分界面，根据所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值，包括：
拾取所述介质分界面，选取空间时窗；
在所述空间时窗内，统计所述介质分界面的振幅属性；
根据所述介质分界面的振幅属性，将所述介质分界面的振幅数据进行归一化处理；
根据所述归一化处理后的介质分界面的振幅数据和所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值。
进一步地，所述方法的另一个实施例中，所述根据所述分界面反射振幅值和预设的地震能量总值，计算所述介质分界面对应的分界面透射振幅值，包括：
根据所述分界面反射振幅值和预设的地震能量总值，利用A＝Ar+At计算所述分界面透射振幅值；
上式中，A表示所述地震能量总值，Ar表示所述分界面反射振幅值，At表示所述分界面透射振幅值。
进一步地，所述方法的另一个实施例中，所述根据所述分界面透射振幅值，计算分界面振幅补偿因子，包括：
根据所述界面透射振幅值，利用计算所述分界面振幅补偿因子；
上式中，St表示所述分界面振幅补偿因子，As表示指定的归一化振幅级别值，At表示所述分界面透射振幅值。
进一步地，所述方法的另一个实施例中，所述根据所述分界面振幅补偿因子和所述地震道的均方根振幅值，进行所述介质分界面的振幅补偿操作，包括：
根据所述分界面振幅补偿因子，利用计算所述介质分界面下伏层补偿振幅值，进行所述介质分界面的振幅补偿操作；
上式中，St表示所述分界面振幅补偿因子，A21表示所述介质分界面下伏层振幅输入值，A22表示所述介质分界面下伏层补偿振幅值。
另一方面，本申请提供了一种介质分界面振幅补偿装置，包括：
基础数据处理模块，用于根据待分析地震数据，获取地震道的均方根振幅值；
反射振幅值计算模块，用于拾取介质分界面，根据所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值；
透射振幅值计算模块，用于根据所述分界面反射振幅值和预设的地震能量总值，计算所述介质分界面对应的分界面透射振幅值；
振幅补偿因子计算模块，用于根据所述分界面透射振幅值，计算分界面振幅补偿因子；
振幅补偿模块，用于根据所述分界面振幅补偿因子和所述地震道的均方根振幅值，进行所述介质分界面的振幅补偿操作。
进一步地，所述装置的另一实施例中，所述基础数据处理模块具体用于：
根据所述待分析地震数据，利用计算所述地震道的均方根振幅值；
上式中，Aj表示第j个地震道的均方根振幅，j表示地震道序号，L表示第j个地震道分析时窗内的总样点数，Xi表示第j个地震道分析时窗内第i个样点的振幅值。
进一步地，所述装置的另一个实施例中，所述反射振幅值计算模块包括：
分界面拾取模块，用于拾取所述介质分界面，选取空间时窗；
振幅属性统计模块，用于在所述空间时窗内，统计所述介质分界面的振幅属性；
归一化处理模块，用于根据所述介质分界面的振幅属性，将所述介质分界面的振幅数据进行归一化处理；
振幅值计算模块，用于根据所述归一化处理后的介质分界面的振幅数据和所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值。
进一步地，所述装置的另一个实施例中，所述透射振幅值计算模块具体用于：
根据所述分界面反射振幅值和预设的地震能量总值，利用A＝Ar+At计算所述分界面透射振幅值；
上式中，A表示所述地震能量总值，Ar表示所述分界面反射振幅值，At表示所述分界面透射振幅值。
进一步地，所述装置的另一个实施例中，所述振幅补偿因子计算模块具体用于：
根据所述界面透射振幅值，利用计算所述分界面振幅补偿因子；
上式中，St表示所述分界面振幅补偿因子，As表示指定的归一化振幅级别值，At表示所述分界面透射振幅值。
进一步地，所述装置的另一个实施例中，所述振幅补偿模块具体用于：
根据所述分界面振幅补偿因子，利用计算所述介质分界面下伏层补偿振幅值，进行所述介质分界面的振幅补偿操作；
上式中，St表示所述分界面振幅补偿因子，A21表示所述介质分界面下伏层振幅输入值，A22表示所述介质分界面下伏层补偿振幅值。
再一方面，本申请还提供了一种介质分界面振幅补偿装置，包括：处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现上述介质分界面振幅补偿方法。
本申请提供的介质分界面振幅补偿方法及装置，地震波在介质分界面主要是反射和透射，预设介质分界面处的地震能量总值不变。根据待分析地震数据，通过介质分界面的拾取，确定出介质分界面处的分界面反射振幅值。根据介质分界面处的分界面反射振幅值和预先确定的介质分界面处的地震能量总值，可以计算出介质分界面处的分界面透射振幅值的变化。介质分界面处的振幅畸变主要是由于地震波在介质分界面的透射带来的，基于介质分界面地震能量的重新分配，计算分界面振幅补偿因子。根据分界面振幅补偿因子进行介质分界面处的振幅补偿操作，可以实现介质分界面沿程的振幅补偿。基于介质分界面处地震能量的重新分配，进行介质分界面处的振幅补偿，提高了介质分界面处振幅补偿的准确性，为后续地震资料的处理和解释提供了准确的数据基础。减少了介质分界面处振幅强弱特征变化不均匀的现象，使得补偿后的地震数据振幅强弱变化更加自然，进一步使得振幅补偿后的地震数据更加真实的反映地下地质的变化。为后续地震资料的处理和解释提供了准确的数据基础，提高了地震数据处理结果的准确性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本申请提供的一种介质分界面振幅补偿方法一个实施例的方法流程示意图；
图2是本申请一个实施例中介质分界面处的振幅曲线示意图；
图3是本申请一个实施例中归一化处理后介质分界面处的振幅曲线示意图；
图4(a)-图4(b)是本申请一个实施例中某工区三维补偿前后的剖面对比示意图；
图5是本申请一个实施例中某工区补偿前后补偿后的振幅曲线对比示意图；
图6是本申请提供的介质分界面振幅补偿装置一个实施例的模块结构示意图；
图7是本申请一个实施例中反射振幅值计算模块的结构示意图；
图8是本申请提供的另一种介质分界面振幅补偿装置实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
在进行油气的勘探时，通常可以采用人工方法激发地震波，在各个检波点处接收地震数据。通过对地震数据的处理，形成能够真实反映地下反射体的图像，通过图像处理分析，识别地下构造的形态和地层的波阻抗特征，判断是否储存有油气。但是，由于地质结构比较复杂，地震波会因地质结构的不同导致地震能量的变化。例如：有些地质中包括煤层结构，由于煤层存在透射损失，且煤层的透射损失是与煤层厚度成正比，使得有无煤层的地方对下覆层的振幅影响比较大。正是由于煤层横向变化比较大的地质特征，使得地震资料中煤层界面处的能量重新分配。地震剖面特征是在有无煤层的地方振幅差异非常大，有煤层的下伏层的能量比较弱，没有煤层的下伏层能量比较强，给地震解释和储层反演都造成很大困难。
在实施本申请的过程中，发现对含有介质分界面的地震数据，地震能量由于地质分界面的影响，在介质分界面上进行重新分配。分界面波阻抗比大时，反射能量强，透射能量弱；反之，分界面波阻抗比小时，透射能量强，反射能量弱。当波阻抗比等于1，反射能量为0，所有能量都透射，随着波阻抗的比值趋于无穷大，透射能量就会接近于零，而反射能量会接近于1。可以认为介质分界面处的地震能量的总值为常数，可以通过分析地震数据，获得地震分界面处的透射能量和反射能量之间的关系，进一步确定介质分界面处振幅的变化，进行介质分界面处的振幅补偿。
图1是本申请提供的一种介质分界面振幅补偿方法一个实施例的方法流程示意图，本申请提供的介质分界面振幅补偿方法包括：
S1、根据待分析地震数据，获取地震道的均方根振幅值。
具体地，在进行油气勘探时，可以采用人工方法激发地震波，在各个检波点处接收地震数据。在每个观测点上记录地震波，都必须经过检波器、放大系统和记录系统三个基本环节，它们连在一起可以总称为地震道。每个采样点采集到的地震数据中可以包括该采样点处的地震波的振幅值。本申请一个实施例可以根据每个采样点获得的地震数据中的振幅值，计算地震道的均方根振幅值。其中每一个检波点获得的地震数据，可以认为是经过地层结构反射获得的数据。因此，地震道对应的均方根振幅值也可以认为是地震道对应的反射振幅值。
S2、拾取介质分界面，根据所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值。
可以通过获得的待分析地震数据，拾取介质分界面。本申请一个实施例中可以根据待分析地震数据的解释分析，确定介质分界面的在待分析地震数据中的时空位置，即可以确定介质分界面在地震道中对应的时空位置。步骤S1获得的地震道对应的均方根振幅值，也可以认为是地震道对应的反射振幅值。确定出介质分界面在地震道中对应的时空位置后，进一步可以根据获得的地震道的均方根振幅值，确定出介质分界面对应的分界面反射振幅值。
S3、根据所述分界面反射振幅值和预设的地震能量总值，计算所述介质分界面对应的分界面透射振幅值。
具体地，由能量守恒原理，可以认为介质分界面处的地震能量总值是固定的，即可以认为地震波在经过介质分界面的上覆层和下伏层前后的总能量不变。地震能量总值的具体数值可以根据实际情况确定，本申请不作具体限定。例如：可以通过实验或历史的地震资料，确定在待分析地质处，地震能量总值的数值。本申请一个实施例中，可以认为地震波在介质分界面处经过透射和反射。振幅可以表示地震波的能量，确定出介质分界面处地震能量总值和介质分界面处的分界面反射振幅值后，可以通过实验获得分界面透射振幅值与介质分界面处地震能量总值、介质分界面处的分界面反射振幅值之间的函数关系，计算分界面透射振幅值。
本申请一个实施例中可以参考如下公式(1)计算介质分界面对应的分界面透射振幅值：
A＝Ar+At(1)
上式中，A可以表示介质分界面处的地震能量总值，Ar可以表示介质分界面处的分界面反射振幅值，At可以表示介质分界面处的分界面透射振幅值。
S4、根据所述分界面透射振幅值，计算分界面振幅补偿因子。
具体地，介质分界面的振幅畸变，可以认为是地震波在介质分界面处透射带来的影响。可以通过实验或根据历史的地震资料，分析地震波在介质分界面处的透射振幅值对最后输出的振幅值的影响。根据分析获得的透射对介质分界面处输出振幅值的影响，获得介质分界面处对应的介质分界面透射振幅值后，可以根据分界面透射振幅值确定出分界面振幅补偿因子。例如：可以通过分析实验数据或历史的地震资料，发现将分界面处的分界面透射振幅值的K倍与分界面处输入的均方根振幅值相乘，可以得到期望输出的分界面处的振幅值。则可以将分界面透射振幅值与K的乘积作为分界面振幅补偿因子。或者，可以预设分界面透射振幅值与分界面振幅补偿因子之间的函数关系，通过实验、数值模拟等方式验证预设的函数关系是否合理，并进行调整，确定合适的函数关系，进一确定分界面振幅补偿因子。确定分界面振幅补偿因子的具体方法，本申请不作具体限定。
S5、根据所述振幅补偿因子和所述地震道的均方根振幅值，进行所述介质分界面的振幅补偿操作。
具体地，计算出分界面振幅补偿因子后，可以根据分界面振幅补偿因子对地震道对应的均方根振幅值进行振幅补偿操作。本申请一个实施例中可以根据分界面振幅补偿因子与期望输出的分界面处的振幅值之间的关系，进行振幅补偿操作。例如：若步骤S4中确定的分界面振幅补偿因子与期望输出的振幅值成比例关系，则可以将分界面处对应的均方根振幅值(或者可以称为分界面反射振幅值)与分界面振幅补偿因子相乘，进行分界面处的振幅补偿操作。
本申请提供的介质分界面振幅补偿方法，根据待分析地震数据，确定地震道对应的均方根振幅值，通过介质分界面的拾取，确定出介质分界面处的分界面反射振幅值。预设介质分界面处的地震能量总值不变，根据介质分界面处的分界面反射振幅值和预先确定的介质分界面处的地震能量总值，可以计算出介质分界面处的分界面透射振幅值。进一步，根据介质分界面透射振幅值，计算分界面振幅补偿因子，根据分界面振幅补偿因子进行介质分界面处的振幅补偿操作。提供了一种基于介质分界面处能量重新分配的振幅补偿方法，可以实现介质分界面沿程的振幅补偿。基于介质分界面处地震能量的重新分配，提高了介质分界面处振幅补偿的准确性，为后续地震资料的处理和解释提供了准确的数据基础，提高了地震数据的处理精度。
在上述实施例的基础上，所述根据待分析地震数据，获取地震道的均方根振幅值，包括：
根据所述待分析地震数据，利用计算所述地震道的均方根振幅值；
上式中，Aj表示第j个地震道的均方根振幅，j表示地震道序号，L表示第j个地震道分析时窗内的总样点数，Xi表示第j个地震道分析时窗内第i个样点的振幅值。
具体地，可以通过获取个采样点处采集到的地震数据，根据如下公式(2)计算获得灭一个地震道对应的均方根振幅值：
上式中，Aj可以表示第j个地震道的均方根振幅，j可以表示地震道序号，L可以表示第j个地震道(当前地震道)分析时窗内的总样点数，Xi可以表示第j个地震道(当前地震道)分析时窗内第i个样点的振幅值。
本申请提供的介质分界面振幅补偿方法，根据采样点采集到的地震数据，可以准确的计算出地震道对应的均方根振幅值，为后续进行振幅补偿操作提供了准确的数据基础，提高了地震数据的处理精度。
在上述实施例的基础上，所述拾取介质分界面，根据所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值，包括：
拾取所述介质分界面，选取空间时窗；
在所述空间时窗内，统计所述介质分界面的振幅属性；
根据所述介质分界面的振幅属性，将所述介质分界面的振幅数据进行归一化处理；
根据所述归一化处理后的介质分界面的振幅数据和所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值。
具体地，根据采样点处采集到的地震数据计算出地震道对应的均方根振幅值后，可以通过拾取介质分界面，确定介质分界面处的反射振幅值。具体可以拾取介质分界面对应的时空位置，指定空间时窗，在指定的空间时窗内，统计介质分界面的振幅属性。空间时窗可以根据地震资料通过试验确定，能统计出介质分界面处的均方根振幅变化即可。振幅属性可以包括均方根振幅，本申请一个实施例中可以先把振幅值限制在指定数据范围内，但是振幅值的相对关系保持不变。即可以根据介质分界面的振幅属性对介质分界面处的振幅数据进行归一化处理，避免原始的地震数据中振幅值比较大，为后续处理带来不便。
根据归一化处理后的介质分界面的振幅数据和地震道对应的均方根振幅值，计算出介质分界面对应的介质分界面反射振幅值。例如：可以确定出介质分界面的时空位置后，将归一化处理后的介质分界面的振幅数据代入上述公式(2)，计算出的均方根振幅Aj即可以表示介质分界面处的分界面反射振幅值。
图2是本申请一个实施例中介质分界面处的振幅曲线示意图，图3是本申请一个实施例中归一化处理后介质分界面处的振幅曲线示意图，如图2和图3所示，经过归一化处理后，介质分界面处的振幅值控制在较小的范围内，方便了后续数据的处理和计算。
本申请提供的介质分界面振幅补偿方法，通过拾取介质分界面，并对介质分界面处的振幅数据进行归一化处理，方便后续介质分界面处振幅补偿数据的处理。根据归一化处理后介质分界面处的振幅数据，计算介质分界面处的反射振幅值。我后续介质分界面处振幅补偿提供了准确的数据基础，同时提高了地震数据处理的速度。
在上述实施例的基础上所述根据所述分界面透射振幅值，计算分界面振幅补偿因子，包括：
根据所述界面透射振幅值，利用计算所述分界面振幅补偿因子；
上式中，St表示所述分界面振幅补偿因子，As表示指定的归一化振幅级别值，At表示所述分界面透射振幅值。
具体地，本申请一个实施例中可以根据如下公式(3)计算分界面振幅补偿因子：
上式中，St可以表示介质分界面处的分界面振幅补偿因子。As可以表示指定的归一化振幅级别值，可以根据需要预先设置，具体可以根据介质分界面处的波阻抗比进行设置，本申请不作具体限定。At可以表示上述实施例计算获得的介质分界面处的分界面透射振幅值。
可以看出，将计算获得的介质分界面处的分界面透射振幅值At代入上述公式(3)，并根据需要设置归一化振幅级别值As后，可以根据公式(3)计算获得分界面振幅补偿因子St。
在上述实施例的基础上，所述根据所述分界面振幅补偿因子和所述地震道的均方根振幅值，进行所述介质分界面的振幅补偿操作，包括：
根据所述分界面振幅补偿因子，利用计算所述介质分界面下伏层补偿振幅值，进行所述介质分界面的振幅补偿操作；
上式中，St表示所述分界面振幅补偿因子，A21表示所述介质分界面下伏层振幅输入值，A22表示所述介质分界面下伏层补偿振幅值。
具体地，计算获得分界面振幅补偿因子St后，可以通过拾取介质分界面下伏层，确定介质分界面下伏层的时空位置。利用上述实施例中计算获得的地震道对应的均方根振幅值，根据下伏层在地震道中对应的时空位置，计算获得介质分界面下伏层对应的下伏层振幅输入值。根据分界面振幅补偿因子和介质分界面下伏层振幅输入值，可以利用如下公式(4)计算获得介质分界面下伏层补偿振幅值，完成介质分界面的振幅补偿。
上式中，St可以表示分界面振幅补偿因子，A21可以表示介质分界面下伏层振幅输入值，可以根据地震道对应的均方根振幅值计算获得，A22可以表示介质分界面下伏层补偿振幅值。
本申请提供的介质分界面振幅补偿方法，基于介质分界面处地震能量的重新分配，提供了计算介质分界面处分界面振幅补偿因子的准确方法。根据计算获得的分界面振幅补偿因子，进行介质分界面处振幅的补偿，可以沿介质分界面进行振幅的补偿，获得沿介质分界面空间位置的振幅补偿值。为后续地震资料的处理和解释提供了准确的数据基础，提高了地震数据处理结果的准确性。
在本申请一个示例中，介质分界面振幅补偿方法可以参考如下的实施过程：
1)计算地震道的均方根振幅。具体可以参考上述公式(2)计算获得每一个地震道的均方根振幅值。
2)对介质分界面进行均方根振幅进行统计，并对均方根振幅曲线归一化处理。可以通过拾取介质分界面，在指定的空间时窗内，对介质分界面进行均方跟振幅值的统计，统计结果可以参考图2所示。根据介质分界面处的振幅属性，对介质分界面处的振幅曲线进行归一化处理，方便后续的数据处理。根据归一化处理后的介质分界面处的振幅数据，利用上述步骤1)计算获得的地震道的均方根振幅值，可以计算获得介质分界面处的均方根振幅值，即分界面反射振幅值。
3)根据地震能量在介质分界面上重新分配，计算介质分界面处的分界面透射振幅值。具体可以参考上述实施例中的公式(1)，计算获得分界面透射振幅值。公式(1)中的地震能量总值A可以根据需要预先设定，分界面反射振幅值Ar可以根据步骤2)利用上述公式(2)计算获得。因此，根据公式(1)可以直接计算获得介质分界面处的分界面透射振幅值，获得介质分界面的透射振幅相对变化。
4)分界面振幅补偿因子计算。根据得到分界面透射振幅值的相对变化，对分界面的下伏层进行振幅畸变因子计算即计算分界面振幅补偿因子。具体可以参考上述实施例中公式(3)计算获得分界面振幅补偿因子。
5)介质分界面振幅补偿处理。根据获得的介质分界面处的分界面振幅补偿因子，对介质分界面进行沿程振幅补偿，具体可以参考上述实施例中的公式(4)计算介质分界面的下伏层补偿振幅值，实现介质分界面处的振幅补偿。
图4(a)-图4(b)是本申请一个实施例中某工区三维补偿前后的剖面对比示意图，图4(a)是本申请一个实施例中一个工区三维补偿前的地震剖面示意图，图4(b)是本申请一个实施例中图4(a)中的工区三维补偿后的地震剖面示意图，图5是本申请一个实施例中某工区补偿前后补偿后的振幅曲线对比示意图。如图4(a)-图4(b)所示，采用本申请的方法，通过剖面对比，经过基于介质分界面的振幅畸变计算和补偿后，介质分界面下伏层的反射振幅特征更加自然。如图5所示，利用本申请介质分界面处振幅补偿方法对介质分界面处的振幅进行补偿后，振幅异常的情况得到比较好的消除。采用本申请的方法，可以将介质分界面处的振幅进行补偿，为后续地震资料的解释和处理，提供了准确的数据基础。
需要说明的是，本申请中的各个公式只是一种实施例，还可以根据需要对公式进行变形或变换，本申请不作具体限定。
本申请提供的介质分界面振幅补偿方法，通过分析介质分界面处的振幅特性，基于地震波能量在介质分界面处的能量重新分配，计算介质分界面处的分界面反射振幅值和分界面透射振幅值。根据介质分界面处的分界面透射振幅值，计算介质分界面对应的分界面振幅补偿因子。根据介质分界面振幅补偿因子，进行介质分界面处的振幅补偿。提高了介质分界面处振幅补偿的准确性，减少了介质分界面处振幅强弱特征变化不均匀的现象，使得补偿后的地震数据振幅强弱变化更加自然，进一步使得振幅补偿后的地震数据更加真实的反映地下地质的变化。为后续地震资料的处理和解释提供了准确的数据基础，提高了地震数据处理结果的准确性。
基于上述所述的介质分界面振幅补偿方法，本说明书一个或多个实施例还提供一种介质分界面振幅补偿装置。所述的装置可以包括使用了本说明书实施例所述方法的系统(包括分布式系统)、软件(应用)、模块、组件、服务器、客户端等并结合必要的实施硬件的装置。基于同一创新构思，本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似，因此本说明书实施例具体的装置的实施可以参见前述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
具体地，图6是本申请提供的介质分界面振幅补偿装置一个实施例的模块结构示意图，如图6所示，本申请中提供的介质分界面振幅补偿装置包括：基础数据处理模块61、反射振幅值计算模块62，透射振幅值计算模块63、振幅补偿因子计算模块64、振幅补偿模块65。
基础数据处理模块61，可以用于根据待分析地震数据，获取地震道的均方根振幅值；
反射振幅值计算模块62，可以用于拾取介质分界面，根据所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值；
透射振幅值计算模块63，可以用于根据所述分界面反射振幅值和预设的地震能量总值，计算所述介质分界面对应的分界面透射振幅值；
振幅补偿因子计算模块64，可以用于根据所述分界面透射振幅值，计算分界面振幅补偿因子；
振幅补偿模块65，可以用于根据所述分界面振幅补偿因子和所述地震道的均方根振幅值，进行所述介质分界面的振幅补偿操作。
本申请提供的介质分界面振幅补偿装置，根据待分析地震数据，确定地震道对应的均方根振幅值，通过介质分界面的拾取，确定出介质分界面处的分界面反射振幅值。预设介质分界面处的地震能量总值不变，根据介质分界面处的分界面反射振幅值和预先确定的介质分界面处的地震能量总值，可以计算出介质分界面处的分界面透射振幅值。进一步，根据介质分界面透射振幅值，计算分界面振幅补偿因子，根据分界面振幅补偿因子进行介质分界面处的振幅补偿操作。提供了一种基于介质分界面处能量重新分配的振幅补偿方法，可以实现介质分界面沿程的振幅补偿，为后续地震资料的处理和解释提供了准确的数据基础，提高了地震数据的处理精度。
在上述实施例的基础上，所述基础数据处理模块具体用于：
根据所述待分析地震数据，利用计算所述地震道的均方根振幅值；
上式中，Aj表示第j个地震道的均方根振幅，j表示地震道序号，L表示第j个地震道分析时窗内的总样点数，Xi表示第j个地震道分析时窗内第i个样点的振幅值。
本申请提供的介质分界面振幅补偿装置，根据采样点采集到的地震数据，可以准确的计算出地震道对应的均方根振幅值，为后续进行振幅补偿操作提供了准确的数据基础，提高了地震数据的处理精度。
图7是本申请一个实施例中反射振幅值计算模块的结构示意图，如图7所示，在上述实施例的基础上，所述反射振幅值计算模块62包括：
分界面拾取模块71，可以用于拾取所述介质分界面，选取空间时窗；
振幅属性统计模块72，可以用于在所述空间时窗内，统计所述介质分界面的振幅属性；
归一化处理模块73，可以用于根据所述介质分界面的振幅属性，将所述介质分界面的振幅数据进行归一化处理；
振幅值计算模块74，可以用于根据所述归一化处理后的介质分界面的振幅数据和所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值。
本申请提供的介质分界面振幅补偿装置，通过拾取介质分界面，并对介质分界面处的振幅数据进行归一化处理，方便后续介质分界面处振幅补偿数据的处理。根据归一化处理后介质分界面处的振幅数据，计算介质分界面处的反射振幅值。我后续介质分界面处振幅补偿提供了准确的数据基础，同时提高了地震数据处理的速度。
在上述实施例的基础上，所述透射振幅值计算模块具体用于：
根据所述分界面反射振幅值和预设的地震能量总值，利用A＝Ar+At计算所述分界面透射振幅值；
上式中，A表示所述地震能量总值，Ar表示所述分界面反射振幅值，At表示所述分界面透射振幅值。
本申请提供的介质分界面振幅补偿装置，基于介质分界面处地震能量的重新分配，计算介质分界面处的分界面反射振幅值和分界面透射振幅值。为后续地震资料的处理和解释提供了准确的数据基础，提高了地震数据处理结果的准确性。
在上述实施例的基础上，所述振幅补偿因子计算模块具体用于：
根据所述界面透射振幅值，利用计算所述分界面振幅补偿因子；
上式中，St表示所述分界面振幅补偿因子，As表示指定的归一化振幅级别值，At表示所述分界面透射振幅值。
本申请提供的介质分界面振幅补偿装置，基于介质分界面处地震能量的重新分配，提供了计算介质分界面处分界面振幅补偿因子的准确方法。根据计算获得的分界面振幅补偿因子，进行介质分界面处振幅的补偿，可以沿介质分界面进行振幅的补偿，获得沿介质分界面空间位置的振幅补偿值。为后续地震资料的处理和解释提供了准确的数据基础，提高了地震数据处理结果的准确性。
在上述实施例的基础上，所述振幅补偿模块具体用于：
根据所述分界面振幅补偿因子，利用计算所述介质分界面下伏层补偿振幅值，进行所述介质分界面的振幅补偿操作；
上式中，St表示所述分界面振幅补偿因子，A21表示所述介质分界面下伏层振幅输入值，A22表示所述介质分界面下伏层补偿振幅值。
本申请提供的介质分界面振幅补偿装置，根据计算获得的分界面振幅补偿因子，进行介质分界面处振幅的补偿，可以沿介质分界面进行振幅的补偿，获得沿介质分界面空间位置的振幅补偿值。为后续地震资料的处理和解释提供了准确的数据基础，减少了介质分界面处振幅强弱特征变化不均匀的现象，使得补偿后的地震数据振幅强弱变化更加自然，进一步使得振幅补偿后的地震数据更加真实的反映地下地质的变化。为后续地震资料的处理和解释提供了准确的数据基础，提高了地震数据处理结果的准确性。
需要说明的是，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
本说明书实施例提供的上述介质分界面振幅补偿方法或装置可以在计算机中由处理器执行相应的程序指令来实现，如使用windows操作系统的c++语言在PC端实现、linux系统实现，或其他例如使用android、iOS系统程序设计语言在智能终端实现，以及基于量子计算机的处理逻辑实现等。本说明书提供的一种介质分界面振幅补偿装置的另一个实施例中，图8是本申请提供的另一种介质分界面振幅补偿装置实施例的模块结构示意图，如图8所示，本申请另一实施例提供的介质分界面振幅补偿装置可以包括处理器81以及用于存储处理器可执行指令的存储器82，
处理器81和存储器82通过总线83完成相互间的通信；
所述处理器81用于调用所述存储器82中的程序指令，以执行上述各介质分界面振幅补偿方法实施例所提供的方法，例如包括：根据待分析地震数据，获取地震道的均方根振幅值；拾取介质分界面，根据所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值；根据所述分界面反射振幅值和预设的地震能量总值，计算所述介质分界面对应的分界面透射振幅值；根据所述分界面透射振幅值，计算分界面振幅补偿因子；根据所述分界面振幅补偿因子和所述地震道的均方根振幅值，进行所述介质分界面的振幅补偿操作。
需要说明的是说明书上述所述的装置根据相关方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照方法实施例的描述，在此不作一一赘述。本申请中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
本说明书实施例并不局限于必须是符合行业通信标准、标准计算机数据处理和数据存储规则或本说明书一个或多个实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、存储、判断、处理方式等获取的实施例，仍然可以属于本说明书实施例的可选实施方案范围之内。
在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray，FPGA))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(Hardware Description Language，HDL)，而HDL也并非仅有一种，而是有许多种，如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等，目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。
控制器可以按任何适当的方式实现，例如，控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320，存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
虽然本说明书一个或多个实施例提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。第一，第二等词语用来表示名称，而并不表示任何特定的顺序。
为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书一个或多个时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
本发明是参照根据本发明实施例的方法、装置(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储、石墨烯存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
本领域技术人员应明白，本说明书一个或多个实施例可提供为方法、系统或计算机程序产品。因此，本说明书一个或多个实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书一个或多个实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书一个或多个实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本本说明书一个或多个实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书一个或多个实施例的实施例而已，并不用于限制本本说明书一个或多个实施例。对于本领域技术人员来说，本说明书一个或多个实施例可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在权利要求范围之内。
Claims (9)
1.一种介质分界面振幅补偿方法，其特征在于，所述方法包括：
根据待分析地震数据，获取地震道的均方根振幅值；
拾取介质分界面，根据所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值；
根据所述分界面反射振幅值和预设的地震能量总值，利用A＝Ar+At计算所述介质分界面对应的分界面透射振幅值，上式中，A表示所述地震能量总值，Ar表示所述分界面反射振幅值，At表示所述分界面透射振幅值；
根据所述分界面透射振幅值，利用计算分界面振幅补偿因子，上式中，St表示所述分界面振幅补偿因子，As表示指定的归一化振幅级别值，At表示所述分界面透射振幅值；
根据所述分界面振幅补偿因子和所述地震道的均方根振幅值，进行所述介质分界面的振幅补偿操作。
2.如权利要求1所述的一种介质分界面振幅补偿方法，其特征在于，所述根据待分析地震数据，获取地震道的均方根振幅值，包括：
根据所述待分析地震数据，利用计算所述地震道的均方根振幅值；
上式中，Aj表示第j个地震道的均方根振幅，j表示地震道序号，L表示第j个地震道分析时窗内的总样点数，Xi表示第j个地震道分析时窗内第i个样点的振幅值。
3.如权利要求1所述的一种介质分界面振幅补偿方法，其特征在于，所述拾取介质分界面，根据所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值，包括：
拾取所述介质分界面，选取空间时窗；
在所述空间时窗内，统计所述介质分界面的振幅属性；
根据所述介质分界面的振幅属性，将所述介质分界面的振幅数据进行归一化处理；
根据所述归一化处理后的介质分界面的振幅数据和所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值。
4.如权利要求1-3任一项所述的一种介质分界面振幅补偿方法，其特征在于，所述根据所述分界面振幅补偿因子和所述地震道的均方根振幅值，进行所述介质分界面的振幅补偿操作，包括：
根据所述分界面振幅补偿因子，利用计算所述介质分界面下伏层补偿振幅值，进行所述介质分界面的振幅补偿操作；
上式中，St表示所述分界面振幅补偿因子，A21表示所述介质分界面下伏层振幅输入值，A22表示所述介质分界面下伏层补偿振幅值。
5.一种介质分界面振幅补偿装置，其特征在于，包括：
基础数据处理模块，用于根据待分析地震数据，获取地震道的均方根振幅值；
反射振幅值计算模块，用于拾取介质分界面，根据所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值；
透射振幅值计算模块，用于根据所述分界面反射振幅值和预设的地震能量总值，利用A＝Ar+At计算所述介质分界面对应的分界面透射振幅值，上式中，A表示所述地震能量总值，Ar表示所述分界面反射振幅值，At表示所述分界面透射振幅值；
振幅补偿因子计算模块，用于根据所述分界面透射振幅值，利用计算分界面振幅补偿因子，上式中，St表示所述分界面振幅补偿因子，As表示指定的归一化振幅级别值，At表示所述分界面透射振幅值；
振幅补偿模块，用于根据所述分界面振幅补偿因子和所述地震道的均方根振幅值，进行所述介质分界面的振幅补偿操作。
6.如权利要求5所述的一种介质分界面振幅补偿装置，其特征在于，所述基础数据处理模块具体用于：
根据所述待分析地震数据，利用计算所述地震道的均方根振幅值；
上式中，Aj表示第j个地震道的均方根振幅，j表示地震道序号，L表示第j个地震道分析时窗内的总样点数，Xi表示第j个地震道分析时窗内第i个样点的振幅值。
7.如权利要求5所述的一种介质分界面振幅补偿装置，其特征在于，所述反射振幅值计算模块包括：
分界面拾取模块，用于拾取所述介质分界面，选取空间时窗；
振幅属性统计模块，用于在所述空间时窗内，统计所述介质分界面的振幅属性；
归一化处理模块，用于根据所述介质分界面的振幅属性，将所述介质分界面的振幅数据进行归一化处理；
振幅值计算模块，用于根据所述归一化处理后的介质分界面的振幅数据和所述地震道的均方根振幅值，计算所述介质分界面对应的分界面反射振幅值。
8.如权利要求5-7任一项所述的一种介质分界面振幅补偿装置，其特征在于，所述振幅补偿模块具体用于：
根据所述分界面振幅补偿因子，利用计算所述介质分界面下伏层补偿振幅值，进行所述介质分界面的振幅补偿操作；
上式中，St表示所述分界面振幅补偿因子，A21表示所述介质分界面下伏层振幅输入值，A22表示所述介质分界面下伏层补偿振幅值。
9.一种介质分界面振幅补偿装置，其特征在于，包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现如权利要求1至4中任意一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711181682.0A CN108181649B (zh) | 2017-11-23 | 2017-11-23 | 一种介质分界面振幅补偿方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711181682.0A CN108181649B (zh) | 2017-11-23 | 2017-11-23 | 一种介质分界面振幅补偿方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108181649A CN108181649A (zh) | 2018-06-19 |
| CN108181649B true CN108181649B (zh) | 2019-09-10 |
Family
ID=62545184
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711181682.0A Active CN108181649B (zh) | 2017-11-23 | 2017-11-23 | 一种介质分界面振幅补偿方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108181649B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111435174B (zh) * | 2019-01-15 | 2023-04-25 | 中国石油天然气集团有限公司 | 强反射地区地震资料振幅补偿方法及装置 |
| CN110095812B (zh) * | 2019-05-09 | 2020-06-23 | 中国海洋石油集团有限公司 | 一种针对浅层气及底辟微裂隙影响的地震振幅补偿方法 |
| CN112649864A (zh) * | 2019-10-12 | 2021-04-13 | 中国石油化工股份有限公司 | 地震剖面变密度分段色标显示方法及系统 |
| CN113267808B (zh) * | 2020-02-14 | 2023-04-25 | 中国石油天然气集团有限公司 | 振幅补偿方法及装置 |
| CN113466943B (zh) * | 2020-03-31 | 2023-04-07 | 中国石油天然气集团有限公司 | 变偏移距vsp地表一致性振幅补偿方法和装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5774417A (en) * | 1996-10-25 | 1998-06-30 | Atlantic Richfield Company | Amplitude and phase compensation in dual-sensor ocean bottom cable seismic data processing |
| CN103645503A (zh) * | 2013-12-17 | 2014-03-19 | 中国海洋石油总公司 | 一种三维时间域照明分析及振幅补偿方法 |
| CN106291691A (zh) * | 2016-08-22 | 2017-01-04 | 中国石油天然气集团公司 | 一种地震偏移成像方法及装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101182838B1 (ko) * | 2010-08-24 | 2012-09-14 | 서울대학교산학협력단 | 송신원 추정을 통한 주파수 영역 역시간 구조보정 방법 및 장치 |
-
2017
- 2017-11-23 CN CN201711181682.0A patent/CN108181649B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5774417A (en) * | 1996-10-25 | 1998-06-30 | Atlantic Richfield Company | Amplitude and phase compensation in dual-sensor ocean bottom cable seismic data processing |
| CN103645503A (zh) * | 2013-12-17 | 2014-03-19 | 中国海洋石油总公司 | 一种三维时间域照明分析及振幅补偿方法 |
| CN106291691A (zh) * | 2016-08-22 | 2017-01-04 | 中国石油天然气集团公司 | 一种地震偏移成像方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 纵波和转换波联合AVO反演方法研究;陈天胜 等;《中国石油大学学报（自然科学版）》;20061231;第30卷(第1期);第33-37页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108181649A (zh) | 2018-06-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108181649B (zh) | 一种介质分界面振幅补偿方法及装置 | |
| Trorey | Diffractions for arbitrary source-receiver locations | |
| CN109100796A (zh) | 一种天然气水合物地震数据处理方法及装置 | |
| CN108983285B (zh) | 一种基于矩张量的多种震源波场模拟方法及装置 | |
| CN103210323B (zh) | 处理地震数据的系统和方法 | |
| CN106873030B (zh) | 断层封堵性的确定方法和装置 | |
| CN109061734A (zh) | 一种地震数据的处理方法及装置 | |
| CN107728210A (zh) | 多仪器采集数据中缺失道的确定方法和装置 | |
| CN107797154A (zh) | 一种电成像测井图像刻度方法及装置 | |
| CN108089238A (zh) | 一种利用综合地球物理资料确定深层裂谷的方法及装置 | |
| CN108828668A (zh) | 一种叠前时间偏移数据处理方法及装置 | |
| Ma et al. | Wavefield extrapolation in pseudodepth domain | |
| CN108107471B (zh) | 一种分方位初至数据体的获取方法及装置 | |
| Olsen et al. | Constraints of crustal heterogeneity and Q (f) from regional (< 4 Hz) wave propagation for the 2009 North Korea nuclear test | |
| CN104459774B (zh) | 一种地质岩性差异识别方法及系统 | |
| CN107656315B (zh) | 储层流体因子的确定方法和装置 | |
| CN108873085A (zh) | 一种沉积岩薄层电磁数据处理方法及装置 | |
| CN108693557A (zh) | 一种地震数据处理方法及装置 | |
| Hidalgo-Gato et al. | The monogenic signal of potential-field data: A Python implementation | |
| CN108710149A (zh) | 一种薄储层的预测方法及装置 | |
| CN109100803A (zh) | 微断裂的确定方法和装置 | |
| CN107390271B (zh) | 地震资料的振幅补偿方法和装置 | |
| CN107807409B (zh) | 地层密度和电阻率关系的确定方法和装置 | |
| Christensen et al. | Seismically driven reservoir characterization using an innovative integrated approach: Syd Arne Field | |
| CN107991705B (zh) | 基于二维统计特征的测井曲线校正方法和装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |