CN114460636A - 基于炮检对信噪比地震照明方法、装置、电子设备及介质 - Google Patents
基于炮检对信噪比地震照明方法、装置、电子设备及介质 Download PDFInfo
- Publication number
- CN114460636A CN114460636A CN202011132717.3A CN202011132717A CN114460636A CN 114460636 A CN114460636 A CN 114460636A CN 202011132717 A CN202011132717 A CN 202011132717A CN 114460636 A CN114460636 A CN 114460636A
- Authority
- CN
- China
- Prior art keywords
- signal
- noise ratio
- shot
- seismic
- gather
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 238000011156 evaluation Methods 0.000 claims abstract description 16
- 238000005286 illumination Methods 0.000 claims description 56
- 238000012545 processing Methods 0.000 claims description 24
- 230000008859 change Effects 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 5
- 238000007689 inspection Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 21
- 238000013461 design Methods 0.000 abstract description 20
- 238000005457 optimization Methods 0.000 abstract description 10
- 238000012937 correction Methods 0.000 abstract description 7
- 230000003068 static effect Effects 0.000 abstract description 7
- 238000011158 quantitative evaluation Methods 0.000 abstract description 3
- 230000005284 excitation Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 8
- 239000002360 explosive Substances 0.000 description 8
- 239000010410 layer Substances 0.000 description 7
- 238000004519 manufacturing process Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 235000019738 Limestone Nutrition 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000010459 dolomite Substances 0.000 description 2
- 229910000514 dolomite Inorganic materials 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 239000006028 limestone Substances 0.000 description 2
- 239000010453 quartz Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 206010013883 Dwarfism Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
Images
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
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Manufacturing & Machinery (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请公开了一种基于炮检对信噪比地震照明方法、装置、电子设备及介质。该方法可以包括：建立炮检对三维信噪比模型；确定标准道或标准道集；根据标准道或标准道集计算模型参数并进行评价。本发明适用于地震采集观测系统设计、采集效果评估和地震资料静校正、去噪处理效果评估等，为观测系统设计优化和资料处理效果的定量评价等提供了方法和依据，具有良好的应用前景。
Description
技术领域
本发明涉及地震勘探技术领域，更具体地，涉及一种基于炮检对信噪比地震照明方法、装置、电子设备及介质。
背景技术
应用反射地震方法获取地下地质结构及储层含油气性特征是油气勘探过程中最重要的手段之一，包括地震采集、地震资料处理和地震解释三个子过程。在这一过程中，对深层地质结构的客观性观测和储层的预测性描述会受到诸多因素影响，主要的影响来自四个方面：一是近地表结构对波传播的影响，复杂的近地表结构导致了波场的复杂，降低了地震采集资料信噪比；其次是低信噪比资料条件下，地下构造的复杂程度对速度估计与偏移成像精度的影响；三是储层含油气性的变化对油气预测准确度的影响；四是有限的经济投入等。这些影响，在很大程度上聚焦于地震采集工作，尤其是地震采集观测系统设计研究，一方面，相对于地震资料处理和地震解释，地震采集具有成本高和不可重复性的特点；另一方面，观测系统设计与优化的目标就是针对地震地质条件复杂性，获得消除或改善这些影响的条件与潜力，并有效地降低地震采集成本。
基于模型的地震照明分析是地震采集观测系统设计中一项关键的技术。技术的发展总体上经历了从简单二维水平层状介质模型到非均匀全三维模型；由简单射线、高斯射线束到单程波方程、全波动方程模拟；由采用CMP覆盖次数、CRP覆盖次数描述到照明能量描述等过程。90年代以前，地震采集观测系统设计一般基于均匀速度场、水平反射层和对称地震射线等假设，地震照明采用CMP覆盖次数描述。90年代后期，射线(束)方法和波动方程方法被逐步引入到地震照明分析中，基于非均匀介质全三维模型地震照明技术，及基于目的层面元的CRP覆盖次数与照明能量等参数描述，成为地震采集观测系统设计、评价和优化的主要方法手段，并取得一定的效果。对于一个已知的地震地质模型，地震照明为确定的观测系统提供了对深层地质结构观测能力的定量评价方法和优化的依据，广泛应用于地震采集观测系统设计与优化评价。
生产的实际需求推进了地震照明技术和软件的发展。在方法研究方面，除射线覆盖次数以外，单程波方程地震照明方法，全波方程地震照明方法，高斯射线束地震照明方法等，满足了实际生产的不同需求。
在实际生产应用中，采用射线照明与描述使得针对深层地质结构复杂、或浅表层激发岩性多变的观测系统设计，其结果缺乏严密性和可信度，对地震资料的处理效果缺乏预见性。
因此，有必要开发一种基于炮检对信噪比地震照明方法、装置、电子设备及介质。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种基于炮检对信噪比地震照明方法、装置、电子设备及介质，其能够适用于地震采集观测系统设计、采集效果评估和地震资料静校正、去噪处理效果评估等，为观测系统设计优化和资料处理效果的定量评价等提供了方法和依据，具有良好的应用前景。
第一方面，本公开实施例提供了一种基于炮检对信噪比地震照明方法，包括：
建立炮检对三维信噪比模型；
确定标准道或标准道集；
根据标准道或标准道集计算模型参数并进行评价。
优选地，模型参数为平均信噪比或去噪后与去噪前道集信号比的变化值。
优选地，计算所述平均信噪比包括：
计算叠后信噪比，进而计算叠加道集信号的平均振幅的叠后信噪比；
建立有效反射能量方程，计算所述平均信噪比。
优选地，通过公式(1)计算叠加道集信号的平均振幅的叠后信噪比：
其中，m为覆盖次数，Sa叠加道集信号的平均振幅，N为处理后噪声的振幅，噪音振幅为
优选地，所述有效反射能量方程为：
S＝m×Sna (4)
优选地，计算所述去噪后与去噪前道集信号比的变化值包括：
计算叠后信噪比，进而计算叠加道集信号的平均振幅的叠后信噪比；
计算去噪后与去噪前道集信号比的变化值。
优选地，通过公式(5)计算去噪后与去噪前道集信号比的变化值：
其中，K为去噪后与去噪前道集信号比的变化值，表示地震资料去噪处理能力。
作为本公开实施例的一种具体实现方式，
第二方面，本公开实施例还提供了一种基于炮检对信噪比地震照明装置，包括：
建模模块，建立炮检对三维信噪比模型；
标准模块，确定标准道或标准道集；
评价模块，根据标准道或标准道集计算模型参数并进行评价。
优选地，模型参数为平均信噪比或去噪后与去噪前道集信号比的变化值。
优选地，计算所述平均信噪比包括：
计算叠后信噪比，进而计算叠加道集信号的平均振幅的叠后信噪比；
建立有效反射能量方程，计算所述平均信噪比。
优选地，通过公式(1)计算叠加道集信号的平均振幅的叠后信噪比：
其中，m为覆盖次数，Sa叠加道集信号的平均振幅，N为处理后噪声的振幅，噪音振幅为
优选地，所述有效反射能量方程为：
S＝m×Sna (4)
其中，Sna为叠加道集的平均信噪比。
优选地，计算所述去噪后与去噪前道集信号比的变化值包括：
计算叠后信噪比，进而计算叠加道集信号的平均振幅的叠后信噪比；
计算去噪后与去噪前道集信号比的变化值。
优选地，通过公式(5)计算去噪后与去噪前道集信号比的变化值：
其中，K为去噪后与去噪前道集信号比的变化值，表示地震资料去噪处理能力。
第三方面，本公开实施例还提供了一种电子设备，该电子设备包括：
存储器，存储有可执行指令；
处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的基于炮检对信噪比地震照明方法。
第四方面，本公开实施例还提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的基于炮检对信噪比地震照明方法。
其有益效果在于：
1、建立了剖面品质、有效覆盖次数与道集信噪比的量化关系，理论基础扎实，物理意义明确，更符合地震勘探实际。
2、涵盖了射线覆盖次数和能量照明的内容，基于射线理论的观测系统设计，其覆盖次数是假设道集平均信噪比为1的地震信噪比照明；而高斯束和波动方程的能量照明，是假设噪声为1的地震信噪比照明，方法应用更具有广泛性。
3、实现了过程的量化设计与评估。通过炮点道集信噪比差异评价激发接收条件，优化观测系统设计；通过对比单炮去噪后与去噪前信噪比与剖面品质、有效覆盖次数量化关系，评价静校正、去噪等处理效果等，初步实现了地震采集和资料处理的量化设计与评估。
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的一个实施例的基于炮检对信噪比地震照明方法的步骤的流程图。
图2a和图2b分别示出了根据本发明的一个实施例的NJ地质图与三维地震单炮信噪比的对比示意图。
图3a和图3b分别示出了根据本发明的一个实施例的炸药震源与可控震源的单炮对比的示意图。
图4示出了根据本发明的一个实施例的一种基于炮检对信噪比地震照明装置的框图。
附图标记说明：
201、建模模块；202、标准模块；203、评价模块。
具体实施方式
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
本发明提供一种基于炮检对信噪比地震照明方法，包括：
建立炮检对三维信噪比模型；
确定标准道或标准道集；
根据标准道或标准道集计算模型参数并进行评价。
在一个示例中，模型参数为平均信噪比或去噪后与去噪前道集信号比的变化值。
在一个示例中，计算平均信噪比包括：
计算叠后信噪比，进而计算叠加道集信号的平均振幅的叠后信噪比；
建立有效反射能量方程，计算平均信噪比。
在一个示例中，通过公式(1)计算叠加道集信号的平均振幅的叠后信噪比：
其中，m为覆盖次数，Sa叠加道集信号的平均振幅，N为处理后噪声的振幅，噪音振幅为
在一个示例中，有效反射能量方程为：
S＝m×Sna (4)
其中，Sna为叠加道集的平均信噪比。
在一个示例中，计算去噪后与去噪前道集信号比的变化值包括：
计算叠后信噪比，进而计算叠加道集信号的平均振幅的叠后信噪比；
计算去噪后与去噪前道集信号比的变化值。
在一个示例中，通过公式(5)计算去噪后与去噪前道集信号比的变化值：
其中，K为去噪后与去噪前道集信号比的变化值，表示地震资料去噪处理能力。
具体地，以射线的地震照明为例，最理想的情况，假设信号S稳定不变，同相叠加，叠加后的振幅为mS，m为覆盖次数，对于平稳的随机噪声，N为原始噪声的振幅，计算叠加后信噪比以dB表示有：
等式右端的第一项是原始资料的信噪比分贝数，第二项是叠加后信噪比提高的分贝增量。覆盖次数m与剖面信噪比
原始道集信噪比
存在：
地震资料经过去噪、地表一致性补偿和自动增益处理后，道集中炮检对的信号S稳定不变，信噪比的变化是信号S振幅的变化，道集叠加后的振幅为mSa，m为覆盖次数，Sa叠加道集信号的平均振幅，噪音为平稳的随机噪声，N为处理后噪声的振幅，叠加道集信号的平均振幅的叠后信噪比以dB表示为公式(1)。
如果Sa相对标准信号S的振幅比为K，则有：
其中，去噪后与去噪前道集信号比的变化为公式(5)，公式(8)右端的第三项是去噪处理后剖面信噪比提高的分贝增量。地震剖面信噪比取决于单炮信噪比、有效覆盖次数和地震资料去噪处理的能力。对道集的有效反射或噪音进行归一化处理，不改变道集的信噪比。不失一般性假设，归一化处理后的N能量相等且等于1，叠加后道集的噪音能量：
则有效反射能量方程为公式(2)，叠加道集或剖面的信噪比Snp有：
从上述公式和讨论可以看出，基于射线理论的观测系统设计，其覆盖次数是假设道集平均信噪比为1的地震信噪比照明；而高斯束和波动方程能量照明，是假设噪音能量为1的地震信噪比照明。射线覆盖次数和能量照明是信噪比照明的特殊形式。覆盖次数、能量照明仍是地震照明的一种定量描述方式，与信噪比照明都是地震采集观测系统设计评价中不同的属性参数，且覆盖次数是能量照明和信噪比照明实现的基本条件。地震信噪比照明的优势在于用实际资料的信噪比衡量地震采集效果及资料处理能力，弥补了覆盖次数和能量照明法存在的不足。
建立炮检对三维信噪比模型；与地质模型对应的地震单炮记录不同岩性的炮检对-偏移距的信噪比值。包括去噪前和去噪后两类，去噪后的单炮记录计算信噪比，体现当前技术条件下地震资料处理的基本能力；去噪前炮检对-偏移距的信噪比，则概括了激发因素、接收因素及偏移距大小(衰减)的影响。在实际生产中，激发比接收对地震资料信噪比的影响更大，可采用单炮地震记录去噪后的信噪比代替叠加道集的信噪比，动态监测地震采集质量和评估地震采集效果，有利于采集参数的实时优化。建立炮检对信噪比三维模型，主要应用于观测系统设计。去噪后与去噪前的信噪比模型对比分析，用于评价静校正或去噪处理效果。
确定标准道或标准道集，获得评价道集准确的信噪比值，这是地震信噪比照明的重量依据；理论上，要求标准道无噪音。有声波测井资料，采用合成记录或单炮作为标准道或标准道集。无声波测井资料时，采用去噪处理后，区域信噪比高的地震记录或地震叠加道作为标准道集或标准道。
进而根据标准道或标准道集计算模型参数并进行评价。针对不同的目的有不同的应用。应用于观测系统评价与优化：按三维信噪比模型，计算面元的平均信噪比；计算拟定观测系统去噪前、后成像道集的信噪比，统计道间距与K值及覆盖次数与道集信噪比的变化关系；根据理论计算或实际需要，选择道集信噪比的阈值范围，例如8～10；对低信噪比(<8)照明区域，采用缩小接收道距或加密接收线或炮点，以增加去噪后和增加有效覆盖次数后剖面信噪比的分贝数，通常，当K值变化显著时，采用缩小接收道距的方法，反之，采用增加覆盖次数的方法。而高信噪比(>10)地震照明区域，可采取适当加大激发线距降低覆盖次数或减小激发井深药量、降低单炮记录信噪比等措施，降低地震勘探成本。激发/接收条件评估，计算炮集平均信噪比或K值。静校正或去噪处理效果评估，计算处理前后道集的平均信噪比或K值等。
本发明还提供一种基于炮检对信噪比地震照明装置，包括：
建模模块，建立炮检对三维信噪比模型；
标准模块，确定标准道或标准道集；
评价模块，根据标准道或标准道集计算模型参数并进行评价。
在一个示例中，模型参数为平均信噪比或去噪后与去噪前道集信号比的变化值。
在一个示例中，计算平均信噪比包括：
计算叠后信噪比，进而计算叠加道集信号的平均振幅的叠后信噪比；
建立有效反射能量方程，计算平均信噪比。
在一个示例中，通过公式(1)计算叠加道集信号的平均振幅的叠后信噪比：
其中，m为覆盖次数，Sa叠加道集信号的平均振幅，N为处理后噪声的振幅，噪音振幅为
在一个示例中，有效反射能量方程为：
S＝m×Sna (4)
其中，Sna为叠加道集的平均信噪比。
在一个示例中，计算去噪后与去噪前道集信号比的变化值包括：
计算叠后信噪比，进而计算叠加道集信号的平均振幅的叠后信噪比；
计算去噪后与去噪前道集信号比的变化值。
在一个示例中，通过公式(5)计算去噪后与去噪前道集信号比的变化值：
其中，K为去噪后与去噪前道集信号比的变化值，表示地震资料去噪处理能力。
具体地，以射线的地震照明为例，最理想的情况，假设信号S稳定不变，同相叠加，叠加后的振幅为mS，m为覆盖次数，对于平稳的随机噪声，N为原始噪声的振幅，计算叠加后信噪比以dB表示有公式(6)，等式右端的第一项是原始资料的信噪比分贝数，第二项是叠加后信噪比提高的分贝增量。覆盖次数m与剖面信噪比
原始道集信噪比
存在公式(7)。
地震资料经过去噪、地表一致性补偿和自动增益处理后，道集中炮检对的信号S稳定不变，信噪比的变化是信号S振幅的变化，道集叠加后的振幅为mSa，m为覆盖次数，Sa叠加道集信号的平均振幅，噪音为平稳的随机噪声，N为处理后噪声的振幅，叠加道集信号的平均振幅的叠后信噪比以dB表示为公式(1)。
如果Sa相对标准信号S的振幅比为K，则有公式(8)，其中，去噪后与去噪前道集信号比的变化为公式(5)，公式(8)右端的第三项是去噪处理后剖面信噪比提高的分贝增量。地震剖面信噪比取决于单炮信噪比、有效覆盖次数和地震资料去噪处理的能力。对道集的有效反射或噪音进行归一化处理，不改变道集的信噪比。不失一般性假设，归一化处理后的N能量相等且等于1，叠加后道集的噪音能量为公式(9)。则有效反射能量方程为公式(2)，叠加道集或剖面的信噪比Snp有公式(10)。
从上述公式和讨论可以看出，基于射线理论的观测系统设计，其覆盖次数是假设道集平均信噪比为1的地震信噪比照明；而高斯束和波动方程能量照明，是假设噪音能量为1的地震信噪比照明。射线覆盖次数和能量照明是信噪比照明的特殊形式。覆盖次数、能量照明仍是地震照明的一种定量描述方式，与信噪比照明都是地震采集观测系统设计评价中不同的属性参数，且覆盖次数是能量照明和信噪比照明实现的基本条件。地震信噪比照明的优势在于用实际资料的信噪比衡量地震采集效果及资料处理能力，弥补了覆盖次数和能量照明法存在的不足。
建立炮检对三维信噪比模型；与地质模型对应的地震单炮记录不同岩性的炮检对-偏移距的信噪比值。包括去噪前和去噪后两类，去噪后的单炮记录计算信噪比，体现当前技术条件下地震资料处理的基本能力；去噪前炮检对-偏移距的信噪比，则概括了激发因素、接收因素及偏移距大小(衰减)的影响。在实际生产中，激发比接收对地震资料信噪比的影响更大，可采用单炮地震记录去噪后的信噪比代替叠加道集的信噪比，动态监测地震采集质量和评估地震采集效果，有利于采集参数的实时优化。建立炮检对信噪比三维模型，主要应用于观测系统设计。去噪后与去噪前的信噪比模型对比分析，用于评价静校正或去噪处理效果。
确定标准道或标准道集，获得评价道集准确的信噪比值，这是地震信噪比照明的重量依据；理论上，要求标准道无噪音。有声波测井资料，采用合成记录或单炮作为标准道或标准道集。无声波测井资料时，采用去噪处理后，区域信噪比高的地震记录或地震叠加道作为标准道集或标准道。
进而根据标准道或标准道集计算模型参数并进行评价。针对不同的目的有不同的应用。应用于观测系统评价与优化：按三维信噪比模型，计算面元的平均信噪比；计算拟定观测系统去噪前、后成像道集的信噪比，统计道间距与K值及覆盖次数与道集信噪比的变化关系；根据理论计算或实际需要，选择道集信噪比的阈值范围，例如8～10；对低信噪比(<8)照明区域，采用缩小接收道距或加密接收线或炮点，以增加去噪后和增加有效覆盖次数后剖面信噪比的分贝数，通常，当K值变化显著时，采用缩小接收道距的方法，反之，采用增加覆盖次数的方法。而高信噪比(>10)地震照明区域，可采取适当加大激发线距降低覆盖次数或减小激发井深药量、降低单炮记录信噪比等措施，降低地震勘探成本。激发/接收条件评估，计算炮集平均信噪比或K值。静校正或去噪处理效果评估，计算处理前后道集的平均信噪比或K值等。
本发明还提供一种电子设备，电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述的基于炮检对信噪比地震照明方法。
本发明还提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述的基于炮检对信噪比地震照明方法。
为便于理解本发明实施例的方案及其效果，以下给出四个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
实施例1
图1示出了根据本发明的一个实施例的基于炮检对信噪比地震照明方法的步骤的流程图。
如图1所示，该基于炮检对信噪比地震照明方法包括：步骤101，建立炮检对三维信噪比模型；步骤102，确定标准道或标准道集；步骤103，根据标准道或标准道集计算模型参数并进行评价。
图2a和图2b分别示出了根据本发明的一个实施例的NJ地质图与三维地震单炮信噪比的对比示意图。
地震采集的基本参数：接收参数，检波器Dz10Hz×12，6串2并，沿测线线性组合，组内间距2米；激发参数，三叠系(T)灰岩、白云岩、泥灰岩等岩性，20m×16Kg×1；下侏罗统(J3)泥岩、中侏罗统(J1)石英砂岩及白垩系(K)砂泥岩地层，18m×12Kg×1。炮集信噪比的计算方法采用目的层时间2000-4500ms带通滤波(10，15，60，65)滤波后的平均振幅值与低通(10，15)滤波后的平均振幅值之比。炮集信噪比的总体分布与区域岩性的展布基本一致，灰岩、白云岩、泥灰岩与石英砂岩激发炮集信噪比相对偏低，泥岩激发炮集信噪比相对较高，激发岩性是影响炮集资料信噪比的主要因素。工区中部NJ县城及沿省道一线，采用增加井深和减少药量等工艺，炮集信噪比有明显提升，激发参数有进一步优化潜力。
从图2a、图2b可知，炮集信噪比分布规律与区域岩性的展布基本一致，激发岩性是影响炮集资料信噪比的主要因素。
图3a和图3b分别示出了根据本发明的一个实施例的炸药震源与可控震源的单炮对比的示意图。
某三维工区炸药震源与可控震源的单炮对比，勘探目的层相关信噪比：炸药1.15，可控震源0.78，信噪比之比1.47。高密度三维中试采用12.5m×12.5m面元，炸药震源激发，PSDM处理。中深层目的层强反射，当覆盖次数大于264次后，剖面信噪比的改善不明显；而对于目的层弱反射，覆盖次数达到396次，信噪比有显著提升。相同区域可控震源三维采集试验结论，同一目的层的覆盖次数分别为540次和大于900次，两者理论覆盖次数之比为2.05-2.27，接近炸药震源与可控震源信噪比比值1.47的平方。实际生产采集可控震源的覆盖次数约为炸药震源的3倍，取得了较好的勘探效果，如图3a、图3b所示。
实施例2
图4示出了根据本发明的一个实施例的一种基于炮检对信噪比地震照明装置的框图。
如图4所示，该基于炮检对信噪比地震照明装置，包括：
建模模块201，建立炮检对三维信噪比模型；
标准模块202，确定标准道或标准道集；
评价模块203，根据标准道或标准道集计算模型参数并进行评价。
作为可选方案，模型参数为平均信噪比或去噪后与去噪前道集信号比的变化值。
作为可选方案，计算平均信噪比包括：
计算叠后信噪比，进而计算叠加道集信号的平均振幅的叠后信噪比；
建立有效反射能量方程，计算平均信噪比。
作为可选方案，通过公式(1)计算叠加道集信号的平均振幅的叠后信噪比：
其中，m为覆盖次数，Sa叠加道集信号的平均振幅，N为处理后噪声的振幅，噪音振幅为
作为可选方案，有效反射能量方程为：
S＝m×Sna (4)
其中，Sna为叠加道集的平均信噪比。
作为可选方案，计算去噪后与去噪前道集信号比的变化值包括：
计算叠后信噪比，进而计算叠加道集信号的平均振幅的叠后信噪比；
计算去噪后与去噪前道集信号比的变化值。
作为可选方案，通过公式(5)计算去噪后与去噪前道集信号比的变化值：
其中，K为去噪后与去噪前道集信号比的变化值，表示地震资料去噪处理能力。
实施例3
本公开提供一种电子设备包括，该电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述基于炮检对信噪比地震照明方法。
根据本公开实施例的电子设备包括存储器和处理器。
该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(RAM)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(ROM)、硬盘、闪存等。
该处理器可以是中央处理单元(CPU)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元，并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中，该处理器用于运行该存储器中存储的该计算机可读指令。
本领域技术人员应能理解，为了解决如何获得良好用户体验效果的技术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些公知的结构也应包含在本公开的保护范围之内。
有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此不再赘述。
实施例4
本公开实施例提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的基于炮检对信噪比地震照明方法。
根据本公开实施例的计算机可读存储介质，其上存储有非暂时性计算机可读指令。当该非暂时性计算机可读指令由处理器运行时，执行前述的本公开各实施例方法的全部或部分步骤。
上述计算机可读存储介质包括但不限于：光存储介质(例如：CD－ROM和DVD)、磁光存储介质(例如：MO)、磁存储介质(例如：磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如：存储卡)和具有内置ROM的媒体(例如：ROM盒)。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种基于炮检对信噪比地震照明方法，其特征在于，包括：
建立炮检对三维信噪比模型；
确定标准道或标准道集；
根据标准道或标准道集计算模型参数并进行评价。
2.根据权利要求1所述的基于炮检对信噪比地震照明方法，其中，模型参数为平均信噪比或去噪后与去噪前道集信号比的变化值。
3.根据权利要求1所述的基于炮检对信噪比地震照明方法，其中，计算所述平均信噪比包括：
计算叠后信噪比，进而计算叠加道集信号的平均振幅的叠后信噪比；
建立有效反射能量方程，计算所述平均信噪比。
4.根据权利要求3所述的基于炮检对信噪比地震照明方法，其中，通过公式(1)计算叠加道集信号的平均振幅的叠后信噪比：
其中，m为覆盖次数，Sa叠加道集信号的平均振幅，N为处理后噪声的振幅，噪音振幅为
5.根据权利要求3所述的基于炮检对信噪比地震照明方法，其中，所述有效反射能量方程为：
S＝m×Sna (4)
其中，Sna为叠加道集的平均信噪比。
6.根据权利要求2所述的基于炮检对信噪比地震照明方法，其中，计算所述去噪后与去噪前道集信号比的变化值包括：
计算叠后信噪比，进而计算叠加道集信号的平均振幅的叠后信噪比；
计算去噪后与去噪前道集信号比的变化值。
7.根据权利要求6所述的基于炮检对信噪比地震照明方法，其中，通过公式(3)计算去噪后与去噪前道集信号比的变化值：
其中，K为去噪后与去噪前道集信号比的变化值，表示地震资料去噪处理能力。
8.一种基于炮检对信噪比地震照明装置，其特征在于，包括：
建模模块，建立炮检对三维信噪比模型；
标准模块，确定标准道或标准道集；
评价模块，根据标准道或标准道集计算模型参数并进行评价。
9.一种电子设备，其特征在于，所述电子设备包括：
存储器，存储有可执行指令；
处理器，所述处理器运行所述存储器中的所述可执行指令，以实现权利要求1-7中任一项所述的基于炮检对信噪比地震照明方法。
10.一种计算机可读存储介质，其特征在于，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现权利要求1-7中任一项所述的基于炮检对信噪比地震照明方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011132717.3A CN114460636A (zh) | 2020-10-21 | 2020-10-21 | 基于炮检对信噪比地震照明方法、装置、电子设备及介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011132717.3A CN114460636A (zh) | 2020-10-21 | 2020-10-21 | 基于炮检对信噪比地震照明方法、装置、电子设备及介质 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114460636A true CN114460636A (zh) | 2022-05-10 |
Family
ID=81404897
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011132717.3A Pending CN114460636A (zh) | 2020-10-21 | 2020-10-21 | 基于炮检对信噪比地震照明方法、装置、电子设备及介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114460636A (zh) |
-
2020
- 2020-10-21 CN CN202011132717.3A patent/CN114460636A/zh active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109425896B (zh) | 白云岩油气储层分布预测方法及装置 | |
| Rentsch et al. | Fast location of seismicity: A migration-type approach with application to hydraulic-fracturing data | |
| CN102176054B (zh) | 近地表综合信息处理解释方法 | |
| CN107065013B (zh) | 一种地震尺度下的层速度确定方法及装置 | |
| US11719840B2 (en) | Subsurface wave slowness prediction system | |
| US11209563B2 (en) | Data optimization method and integral prestack depth migration method | |
| CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
| CN111983683B (zh) | 一种少井条件下盐下湖相灰岩储层预测方法和系统 | |
| CN105425299A (zh) | 确定地层裂缝分布的方法和装置 | |
| Ktenidou et al. | Directional dependence of site effects observed near a basin edge at Aegion, Greece | |
| CN112114380A (zh) | 一种用于砂砾岩油藏储层的预测方法 | |
| CN111812716B (zh) | 页岩气储层总有机碳含量叠前定量预测方法、装置及设备 | |
| CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
| CN103513279A (zh) | 一种基于地震波波动方程的照明分析计算方法及计算装置 | |
| CN114460636A (zh) | 基于炮检对信噪比地震照明方法、装置、电子设备及介质 | |
| CN112213774B (zh) | 一种浅层q模型估计方法及装置 | |
| CN114755740A (zh) | 岩石分布确定方法、装置、设备及介质 | |
| CN113156494B (zh) | 面元等效覆盖次数的计算方法及装置 | |
| CN113806674A (zh) | 古河道纵向尺度的量化方法、装置、电子设备及存储介质 | |
| US10067254B2 (en) | Removal of an estimated acquisition effect from a marine survey measurement | |
| US10401515B2 (en) | Estimation of water properties from seismic data | |
| CN113267808B (zh) | 振幅补偿方法及装置 | |
| CN110850504B (zh) | 一种基于铀曲线拟阻抗约束的页岩密度参数叠前反演方法 | |
| RU2488145C1 (ru) | Способ построения сейсмических изображений геологической среды | |
| CN115639602A (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 |