CN113109865A - 混源激发地震波的处理方法及装置 - Google Patents
混源激发地震波的处理方法及装置 Download PDFInfo
- Publication number
- CN113109865A CN113109865A CN202010021168.6A CN202010021168A CN113109865A CN 113109865 A CN113109865 A CN 113109865A CN 202010021168 A CN202010021168 A CN 202010021168A CN 113109865 A CN113109865 A CN 113109865A
- Authority
- CN
- China
- Prior art keywords
- seismic
- source
- explosive
- consistency
- controllable
- 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.)
- Granted
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 14
- 238000012545 processing Methods 0.000 claims abstract description 83
- 230000000694 effects Effects 0.000 claims abstract description 16
- 239000002360 explosive Substances 0.000 claims description 87
- 238000000034 method Methods 0.000 claims description 44
- 238000012937 correction Methods 0.000 claims description 38
- 238000001914 filtration Methods 0.000 claims description 35
- 238000007689 inspection Methods 0.000 claims description 22
- 238000006243 chemical reaction Methods 0.000 claims description 21
- 238000004590 computer program Methods 0.000 claims description 16
- 230000001427 coherent effect Effects 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000009792 diffusion process Methods 0.000 claims description 8
- 230000005284 excitation Effects 0.000 claims description 8
- 238000011835 investigation Methods 0.000 claims description 7
- 238000013459 approach Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims 4
- 238000010586 diagram Methods 0.000 description 15
- 230000006870 function Effects 0.000 description 12
- 238000004891 communication Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 230000004044 response Effects 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 230000003068 static effect Effects 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/02—Generating seismic energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/02—Generating seismic energy
- G01V1/104—Generating seismic energy using explosive charges
-
- 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
技术领域
本发明涉及石油勘探技术领域，具体涉及一种混源激发地震波的处理方法及装置。
背景技术
随着城市化进程的加快以及工农业的迅速发展，华北平原地区地震资料采集难度越来越大。其中最直接的影响就是城市周边及城市内的炮点布设问题。为了缓解施工生产与地方政府及百姓的工农矛盾，同时在国家大力倡导绿色勘探的背景要求下，近年来在华北平原部分区域(尤其是城区)采用可控震源代替炸药震源进行激发，大大缓解了工农矛盾同时保障了施工生产的顺利进行。
混源激发地震波是指既包含可控震源又包含炸药震源激发的地震资料。但由于同一工区内既存在炸药震源又存在可控震源，地震资料之间存在相位、振幅、频率等差异，不做任何处理将严重影响资料品质。因此，不同震源间的一致性处理就成为了混源激发地震波必须要解决的问题。
以往针对混源激发，处理上仅仅是做了时差校正，保证两种震源资料的相位基本一致。这样的处理成果在当今复杂岩性勘探阶段将严重影响地质认识，从而导致失利井出现。
发明内容
针对现有技术中的问题，本发明提供一种混源激发地震波的处理方法，包括：
对可控震源进行最小相位化处理，以使所述可控震源与炸药震源激发出的地震波相位相同；
根据激发出的炮点数量，从所述可控震源和所述炸药震源中确定出其中一个震源为一致性主体，另一个为非一致性主体；
将非一致性主体的震源的振幅级别调节至与一致性主体的震源相同的振幅级别；
对所述可控震源和炸药震源进行时差校正；
以所述一致性主体为基准，对已知时差的所述可控震源和炸药震源进行频率一致性匹配，得到处理后的所述可控震源和所述炸药震源，进而两个震源激发出的地震波的相位、振幅以及频率的差异低于设定值。
一些实施例中，所述对可控震源进行最小相位化处理，包括：
记录可控震源激发的地震波，生成地震波记录数据；
对所述地震波记录数据进行自相关，生成可控震源零相位子波；
利用所述可控震源零相位子波，计算将所述可控震源零相位子波转化为最小相位子波的转换算子；
利用所述转换算子将所述可控震源零相位子波转化为最小相位子波，生成最小相位子波的地震波记录数据。
一些实施例中，所述将非一致性主体的震源的振幅级别调节至与一致性主体的震源相同的振幅级别，包括：
按照所述可控震源与所述炸药震源的能量差异级别以及两者能量值的倍数，将非一致性主体的能量级别补偿至与所述一致性主体相同的能量级别；
通过几何扩散补偿的方式补偿两个震源纵向上的能量衰减；
通过地表一致性振幅补偿的方式补偿两个震源横向上的能量差异，进而使所述可控震源与所述炸药震源的振幅级别一致。
一些实施例中，所述对所述可控震源和炸药震源进行时差校正，包括：
选取重叠段进行时差调查，采用网格化互相关的方式求取时差；
将求取的时差插值成校正量库，进行时差校正，并通过叠加检查时差校正效果。
一些实施例中，所述对已知时差的所述可控震源和炸药震源进行频率一致性匹配，包括：
以一致性主体为目标地震道，以非一致性主体为原始地震道，对所述原始地震道进行匹配滤波使其逼近所述目标地震道；
应用托布里兹矩阵方程求解匹配滤波算子；
用所述匹配滤波算子对原始地震道进行滤波，进而使所述可控震源和炸药震源实现频率的一致性匹配。
一些实施例中，还包括：
对进行一致性匹配后的所述可控震源和所述炸药震源进行叠加检查。
一些实施例中，所述对进行一致性匹配后的所述可控震源和所述炸药震源进行叠加检查，包括：
将所述炸药震源和所述可控震源激发的地震道按照线、道以及偏移距进行分选，得到多个分选道集；
利用动校速度将共深度点的道集叠加成一道，形成剖地震面，进而根据形成的剖地震面；
根据所述剖地震面进行叠加检查操作。
本发明第二方面实施例提供一种混源激发地震波的处理装置，包括：
最小相位化处理模块，对可控震源进行最小相位化处理，以使所述可控震源与炸药震源激发出的地震波相位相同；
一致性主体确定模块，根据激发出的炮点数量，从所述可控震源和所述炸药震源中确定出其中一个震源为一致性主体，另一个为非一致性主体；
振幅级别统一模块，将非一致性主体的震源的振幅级别调节至与一致性主体的震源相同的振幅级别；
时差确定模块，对所述可控震源和炸药震源进行时差校正；
频率一致性匹配模块，以所述一致性主体为基准，对已知时差的所述可控震源和炸药震源进行频率一致性匹配，得到处理后的所述可控震源和所述炸药震源，进而两个震源激发出的地震波的相位、振幅以及频率的差异低于设定值。
一些实施例中，所述最小相位化处理模块，包括：
地震波记录数据生成单元，记录可控震源激发的地震波，生成地震波记录数据；
零相位子波生成单元，对所述地震波记录数据进行自相关，生成可控震源零相位子波；
转换算子计算单元，利用所述可控震源零相位子波，计算将所述可控震源零相位子波转化为最小相位子波的转换算子；
最小相位子波转换单元，利用所述转换算子将所述可控震源零相位子波转化为最小相位子波，生成最小相位子波的地震波记录数据。
一些实施例中，所述振幅级别统一模块，包括：
第一补偿单元，按照所述可控震源与所述炸药震源的能量差异级别以及两者能量值的倍数，将非一致性主体的能量级别补偿至与所述一致性主体相同的能量级别；
第二补偿单元，通过几何扩散补偿的方式补偿两个震源纵向上的能量衰减；
第三补偿单元，通过地表一致性振幅补偿的方式补偿两个震源横向上的能量差异，进而使所述可控震源与所述炸药震源的振幅级别一致。
一些实施例中，所述时差确定模块，包括：
时差求取单元，选取重叠段进行时差调查，采用网格化互相关的方式求取时差；
校正单元，将求取的时差插值成校正量库，进行时差校正，并通过叠加检查时差校正效果。
一些实施例中，所述频率一致性匹配模块，包括：
匹配滤波单元，以一致性主体为目标地震道，以非一致性主体为原始地震道，对所述原始地震道进行匹配滤波使其逼近所述目标地震道；
匹配滤波算子求解单元，应用托布里兹矩阵方程求解匹配滤波算子；
用所述匹配滤波算子对原始地震道进行滤波，进而使所述可控震源和炸药震源实现频率的一致性匹配。
一些实施例中，还包括：
叠加检查模块，对进行一致性匹配后的所述可控震源和所述炸药震源进行叠加检查。
一些实施例中，所述叠加检查模块，包括：
分选单元，将所述炸药震源和所述可控震源激发的地震道按照线、道以及偏移距进行分选，得到多个分选道集；
剖地震面形成单元，利用动校速度将共深度点的道集叠加成一道，形成剖地震面，进而根据形成的剖地震面；
叠加检查单元，根据所述剖地震面进行叠加检查操作。
第三方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述混源激发地震波的处理方法的步骤。
第四方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现所述混源激发地震波的处理方法的步骤。
由上述技术方案可知，本发明提供一种混源激发地震波的处理方法及装置，首先进行最小相位化处理，最小相位化处理可以保证处理后的可控震源子波就是最小相位，不会影响反褶积效果，然后进行振幅一致性处理，以及进行频率一致性匹配处理，考虑振幅和频率等一致性的影响，实施本发明的技术方案后，由于混源激发造成的子波一致性问题得到很好地解决，为后续处理奠定良好的资料基础。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的实施例中的混源激发地震波的处理方法的流程示意图。
图2示出了可控震源和炸药震源的子波对比示意图。
图3a示出了可控震源最小相位化前拼接剖面及可控震源子波示意图。
图3b示出了可控震源最小相位化后拼接剖面及可控震源子波示意图。
图4a示出了原始叠加剖面。
图4b示出了经过系列去躁，可控震源与炸药震源振幅级别统一后的剖面。
图4c示出了时间和空间上再次对震源进行上述能力级别调整后的剖面。
图5a示出了可控震源最小相位化处理之前的剖图。
图5b示出了可控震源最小相位化处理之后的剖图。
图5c示出了振幅级别统一后效果示意图。
图5d示出了实施本发明的方法之前某区资料叠加效果图。
图5e示出了实施本发明的方法之后某区资料叠加效果图。
图6为本发明的实施例中的混源激发地震波的处理装置的结构示意图。
图7为本发明的实施例中的电子设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的实施例提供一种混源激发地震波的处理方法的具体实施方式，参见图1，所述混源激发地震波的处理方法具体包括如下内容：
S1：对可控震源进行最小相位化处理，以使所述可控震源与炸药震源激发出的地震波相位相同；
S2：根据激发出的炮点数量，从所述可控震源和所述炸药震源中确定出其中一个震源为一致性主体，另一个为非一致性主体；
S3：将非一致性主体的震源的振幅级别调节至与一致性主体的震源相同的振幅级别；
S4：对所述可控震源和炸药震源进行时差校正；
S5：以所述一致性主体为基准，对已知时差的所述可控震源和炸药震源进行频率一致性匹配，得到处理后的所述可控震源和所述炸药震源，进而两个震源激发出的地震波的相位、振幅以及频率的差异低于设定值。
从上述描述可知，本发明的应用实例提供的混源激发地震波的处理方法，首先进行最小相位化处理，最小相位化处理可以保证处理后的可控震源子波就是最小相位，不会影响反褶积效果，然后进行振幅一致性处理，以及进行频率一致性匹配处理，考虑振幅和频率等一致性的影响，实施本发明的技术方案后，由于混源激发造成的子波一致性问题得到很好地解决，为后续处理奠定良好的资料基础。
具体的，步骤S1中，具体包括：
S11：记录可控震源激发的地震波，生成地震波记录数据；
S12：对所述地震波记录数据进行自相关，生成可控震源零相位子波；
S13：利用所述可控震源零相位子波，计算将所述可控震源零相位子波转化为最小相位子波的转换算子；
S14：利用所述转换算子将所述可控震源零相位子波转化为最小相位子波，生成最小相位子波的地震波记录数据。
图2示出了可控震源和炸药震源的子波对比示意图，在一些实施例中，可控震源激发的为零相位子波，而炸药震源为最小相位的，其相位不一致，而反褶积理论模型是针对子波为最小相位来计算的。通过其他手段，如子波整形处理，并不能保证处理后的可控震源子波就是最小相位，会影响反褶积效果，所以对可控震源资料首先进行最小相位化处理可以保证处理后的可控震源子波就是最小相位，不会影响反褶积效果。
具体的，在一些实施例中，从可控震源地震记录采集时使用的扫描信号参数出发，计算可控震源扫描信号，再将可控震源扫描信号进行自相关求出可控震源零相位子波。然后利用可控震源零相位子波，求出将零相位子波转化为最小相位子波的转换算子。最后利用该转换算子，将可控震源零相位子波转化为最小相位子波，实现了把可控震源零相位震源子波的地震记录转化为最小相位震源子波的地震记录的功能，图3a和图3b示出了可控震源最小相位化前后拼接剖面及可控震源子波对比，通过图3a可以看到可控震源区与周边区域存在相位差，图3b可以看出经可控震源最小相位化后，可控震源区与炸药震源区相位保持一致，均为最小相位，为城市与周边三维资料进行拼接奠定了基础。
在步骤S2中，确定一致性主体，哪种震源激发炮点多，将谁当成一致性主体，另一个为非一致性主体，在后续的步骤中，均以一致性主体为基准。
在确定一致性主体后，执行步骤S3，具体的，步骤3具体包括：
S31：按照所述可控震源与所述炸药震源的能量差异级别以及两者能量值的倍数，将非一致性主体的能量级别补偿至与所述一致性主体相同的能量级别；
S32：通过几何扩散补偿的方式补偿两个震源纵向上的能量衰减；
S33：通过地表一致性振幅补偿的方式补偿两个震源横向上的能量差异，进而使所述可控震源与所述炸药震源的振幅级别一致。
具体的，可控震源与井炮激发能量差异大，首先按照能量差异级别乘以差值的倍数，使其能量级别基本统一，然后通过几何扩散补偿补偿纵向上的能量衰减，通过地表一致性振幅补偿补偿横向上的能量差异，最后保证两种震源的振幅级别统一。
在一些实施例中，采用整炮方式统计振幅，根据振幅期望输出逐炮进行能量调整，并不改变炮内振幅关系，为保幅处理，通过此方法，可以将能量级别差异较大的可控震源单炮与炸药震源单炮调整至同一级别，以便于下一步振幅处理；时间与空间上的振幅处理：首先，充分考虑地震波在传播时球面发散及地层的吸收影响对振幅的衰减，采用球面扩散补偿来对时间上的衰减进行振幅恢复，使浅中深层能量保持一致；其次，考虑到由于地表条件以及激发接收条件的变化对炮点响应、接收点响应、炮检距响应以及共中心点响应的综合反映，进行四分量的地表一致性振幅补偿处理，充分解决地震资料能量在空间上的差异。
图4a、图4b以及图4c依次示出了原始叠加剖面、经过系列去躁，可控震源与炸药震源振幅级别统一后，以及时间和空间上再次对震源进行上述能量级别调整后的剖面。图4a可以看出可控震源与炸药震源差异化去噪后，消除了后期振幅处理中噪声的影响；图4b可以看出通过采用整炮方式统计振幅的思路，将能量级别差异较大的可控震源单炮与炸药震源单炮调整至同一级别；图4c可以看出通过时间及空间上振幅补偿后，充分解决地震资料能量在空间上的差异。
在步骤S4中，时差校正的具体步骤为：先选取重叠区位置叠加剖面，通过互相关方式进行时差调查，时差调查后插值成校正量库，按照校正量进行时差校正。
进一步的，进行时差校正后进行步骤S5，具体的，步骤S5包括：
S51：以一致性主体为目标地震道，以非一致性主体为原始地震道，对所述原始地震道进行匹配滤波使其逼近所述目标地震道；
S52：应用托布里兹矩阵方程求解匹配滤波算子；
S53：用所述匹配滤波算子对原始地震道进行滤波，进而使所述可控震源和炸药震源实现频率的一致性匹配。
具体而言，本发明通过具体研究发现，在基本解决了振幅问题和可控震源资料的相位问题后，可控震源与炸药震源的子波仍然存在一定的差异，这种差异更多地表现在频率及剖面特征上。
本发明以匹配滤波法解决子波的一致性问题。匹配滤波利用重复地震道(可分别称为原始道和目标道)设计匹配滤波算子，然后对原始道进行匹配滤波使其最大限度地接近目标道。设原始地震数据道的地震子波为x(t)，目标地震数据道的地震子波为y(t)，匹配滤波算子为a(t)，则应有以下关系：
x(t)×a(t)＝y(t)
转换到频率域，即为X(f)·A(f)＝Y(f)
应用最小平方法求解匹配滤波算子的托布里兹矩阵方程，即Rxx·A＝Rxy，式中，Rxx为原始地震数据道的自相关函数矩阵；Rxy为原始地震数据道x(t)与目标地震数据道y(t)的互相关函数矩阵；A为匹配滤波算子向量。求解即可得到匹配滤波算子A即a(t)，再用确定的匹配滤波算子对原始道数据滤波即可完成匹配滤波，从而实现子波一致性处理。通过上述匹配滤波计算方法，求取可控震源向炸药震源匹配的滤波算子，应用匹配滤波后可控震源与炸药震源的子波基本达到一致，满足了地震资料后续处理的要求，成功地解决了可控震源与炸药震源之间的能量差异、相位差异以及频率差异等，较好地保证了混源激发区地震资料的子波一致性。
进一步的，本发明还包括：
S6：对进行一致性匹配后的所述可控震源和所述炸药震源进行叠加检查。
在一些实施例中，步骤S6具体包括：
S61：将所述炸药震源和所述可控震源激发的地震道按照线、道以及偏移距进行分选，得到多个分选道集；
S62：利用动校速度将共深度点的道集叠加成一道，形成剖地震面，进而根据形成的剖地震面；
S63：根据所述剖地震面进行叠加检查操作。
该实施例中，叠加检查就是将cmp道集按照线、道、偏移距分选，然后利用动较速度将相同的cdp点叠加成一道，形成剖地震面检查效果。如图5a至图5e所示，通过本发明处理过程实施后，可控震源区与井炮区一致性得到明显改善，剖面成像质量显著提高。
从上述描述可以知晓，本发明提供一种混源激发地震波的处理方法，首先进行最小相位化处理，最小相位化处理可以保证处理后的可控震源子波就是最小相位，不会影响反褶积效果，然后进行振幅一致性处理，以及进行频率一致性匹配处理，考虑振幅和频率等一致性的影响，实施本发明的技术方案后，由于混源激发造成的子波一致性问题得到很好地解决，为后续处理奠定良好的资料基础。
本发明的实施例还提供一种用于实现混源激发地震波的处理方法中全部内容的混源激发地震波的处理装置的具体实施方式，参见图6，所述混源激发地震波的处理装置具体包括如下内容：
最小相位化处理模块1，对可控震源进行最小相位化处理，以使所述可控震源与炸药震源激发出的地震波相位相同；
一致性主体确定模块2，根据激发出的炮点数量，从所述可控震源和所述炸药震源中确定出其中一个震源为一致性主体，另一个为非一致性主体；
振幅级别统一模块3，将非一致性主体的震源的振幅级别调节至与一致性主体的震源相同的振幅级别；
时差确定模块4，对所述可控震源和炸药震源进行时差校正；
频率一致性匹配模块5，以所述一致性主体为基准，对已知时差的所述可控震源和炸药震源进行频率一致性匹配，得到处理后的所述可控震源和所述炸药震源，进而两个震源激发出的地震波的相位、振幅以及频率的差异低于设定值。
本申请提供的用于实现混源激发地震波的处理装置的实施例具体可以用于执行上述实施例中的用于实现混源激发地震波的处理方法的实施例的处理流程，其功能在此不再赘述，可以参照上述方法实施例的详细描述。
其中，所述最小相位化处理模块，包括：
地震波记录数据生成单元，记录可控震源激发的地震波，生成地震波记录数据；
零相位子波生成单元，对所述地震波记录数据进行自相关，生成可控震源零相位子波；
转换算子计算单元，利用所述可控震源零相位子波，计算将所述可控震源零相位子波转化为最小相位子波的转换算子；
最小相位子波转换单元，利用所述转换算子将所述可控震源零相位子波转化为最小相位子波，生成最小相位子波的地震波记录数据。
其中，所述振幅级别统一模块，包括：
第一补偿单元，按照所述可控震源与所述炸药震源的能量差异级别以及两者能量值的倍数，将非一致性主体的能量级别补偿至与所述一致性主体相同的能量级别；
第二补偿单元，通过几何扩散补偿的方式补偿两个震源纵向上的能量衰减；
第三补偿单元，通过地表一致性振幅补偿的方式补偿两个震源横向上的能量差异，进而使所述可控震源与所述炸药震源的振幅级别一致。
其中，所述时差确定模块，包括：
时差求取单元，选取重叠段进行时差调查，采用网格化互相关的方式求取时差；
校正单元，将求取的时差插值成校正量库，进行时差校正，并通过叠加检查时差校正效果。
其中，所述频率一致性匹配模块，包括：
匹配滤波单元，以一致性主体为目标地震道，以非一致性主体为原始地震道，对所述原始地震道进行匹配滤波使其逼近所述目标地震道；
匹配滤波算子求解单元，应用托布里兹矩阵方程求解匹配滤波算子；
用所述匹配滤波算子对原始地震道进行滤波，进而使所述可控震源和炸药震源实现频率的一致性匹配。
其中，所述装置还包括：
叠加检查模块，对进行一致性匹配后的所述可控震源和所述炸药震源进行叠加检查。
其中，所述叠加检查模块，包括：
分选单元，将所述炸药震源和所述可控震源激发的地震道按照线、道以及偏移距进行分选，得到多个分选道集；
剖地震面形成单元，利用动校速度将共深度点的道集叠加成一道，形成剖地震面，进而根据形成的剖地震面；
叠加检查单元，根据所述剖地震面进行叠加检查操作。
从上述描述可知，本发明的实施例提供的混源激发地震波的处理装置，首先进行最小相位化处理，最小相位化处理可以保证处理后的可控震源子波就是最小相位，不会影响反褶积效果，然后进行振幅一致性处理，以及进行频率一致性匹配处理，考虑振幅和频率等一致性的影响，实施本发明的技术方案后，由于混源激发造成的子波一致性问题得到很好地解决，为后续处理奠定良好的资料基础。
本申请的实施例还提供能够实现上述实施例中的混源激发地震波的处理方法中全部步骤的一种电子设备的具体实施方式，参见图7，所述电子设备具体包括如下内容：
处理器(processor)601、存储器(memory)602、通信接口(CommunicationsInterface)603和总线604；
其中，所述处理器601、存储器602、通信接口603通过所述总线604完成相互间的通信；所述通信接口603用于实现混源激发地震波的处理装置以及用户终端等相关设备之间的信息传输；
所述处理器601用于调用所述存储器602中的计算机程序，所述处理器执行所述计算机程序时实现上述实施例中的混源激发地震波的处理方法中的全部步骤。
本申请的实施例还提供能够实现上述实施例中的混源激发地震波的处理方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的混源激发地震波的处理方法的全部步骤。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
虽然本说明书实施例提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书实施例的实施例而已，并不用于限制本说明书实施例。对于本领域技术人员来说，本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书实施例的权利要求范围之内。
Claims (16)
1.一种混源激发地震波的处理方法，其特征在于，包括：
对可控震源进行最小相位化处理，以使所述可控震源与炸药震源激发出的地震波相位相同；
根据激发出的炮点数量，从所述可控震源和所述炸药震源中确定出其中一个震源为一致性主体，另一个为非一致性主体；
将非一致性主体的震源的振幅级别调节至与一致性主体的震源相同的振幅级别；
对所述可控震源和炸药震源进行时差校正；
以所述一致性主体为基准，对已知时差的所述可控震源和炸药震源进行频率一致性匹配，得到处理后的所述可控震源和所述炸药震源，进而两个震源激发出的地震波的相位、振幅以及频率的差异低于设定值。
2.根据权利要求1所述的混源激发地震波的处理方法，其特征在于，所述对可控震源进行最小相位化处理，包括：
记录可控震源激发的地震波，生成地震波记录数据；
对所述地震波记录数据进行自相关，生成可控震源零相位子波；
利用所述可控震源零相位子波，计算将所述可控震源零相位子波转化为最小相位子波的转换算子；
利用所述转换算子将所述可控震源零相位子波转化为最小相位子波，生成最小相位子波的地震波记录数据。
3.根据权利要求1所述的混源激发地震波的处理方法，其特征在于，所述将非一致性主体的震源的振幅级别调节至与一致性主体的震源相同的振幅级别，包括：
按照所述可控震源与所述炸药震源的能量差异级别以及两者能量值的倍数，将非一致性主体的能量级别补偿至与所述一致性主体相同的能量级别；
通过几何扩散补偿的方式补偿两个震源纵向上的能量衰减；
通过地表一致性振幅补偿的方式补偿两个震源横向上的能量差异，进而使所述可控震源与所述炸药震源的振幅级别一致。
4.根据权利要求2所述的混源激发地震波的处理方法，其特征在于，所述对所述可控震源和炸药震源进行时差校正，包括：
选取重叠段进行时差调查，采用网格化互相关的方式求取时差；
将求取的时差插值成校正量库，进行时差校正，并通过叠加检查时差校正效果。
5.根据权利要求2所述的混源激发地震波的处理方法，其特征在于，所述对已知时差的所述可控震源和炸药震源进行频率一致性匹配，包括：
以一致性主体为目标地震道，以非一致性主体为原始地震道，对所述原始地震道进行匹配滤波使其逼近所述目标地震道；
应用托布里兹矩阵方程求解匹配滤波算子；
用所述匹配滤波算子对原始地震道进行滤波，进而使所述可控震源和炸药震源实现频率的一致性匹配。
6.根据权利要求1所述的混源激发地震波的处理方法，其特征在于，还包括：
对进行一致性匹配后的所述可控震源和所述炸药震源进行叠加检查。
7.根据权利要求6所述的混源激发地震波的处理方法，其特征在于，所述对进行一致性匹配后的所述可控震源和所述炸药震源进行叠加检查，包括：
将所述炸药震源和所述可控震源激发的地震道按照线、道以及偏移距进行分选，得到多个分选道集；
利用动校速度将共深度点的道集叠加成一道，形成剖地震面，进而根据形成的剖地震面；
根据所述剖地震面进行叠加检查操作。
8.一种混源激发地震波的处理装置，其特征在于，包括：
最小相位化处理模块，对可控震源进行最小相位化处理，以使所述可控震源与炸药震源激发出的地震波相位相同；
一致性主体确定模块，根据激发出的炮点数量，从所述可控震源和所述炸药震源中确定出其中一个震源为一致性主体，另一个为非一致性主体；
振幅级别统一模块，将非一致性主体的震源的振幅级别调节至与一致性主体的震源相同的振幅级别；
时差确定模块，对所述可控震源和炸药震源进行时差校正；
频率一致性匹配模块，以所述一致性主体为基准，对已知时差的所述可控震源和炸药震源进行频率一致性匹配，得到处理后的所述可控震源和所述炸药震源，进而两个震源激发出的地震波的相位、振幅以及频率的差异低于设定值。
9.根据权利要求8所述的混源激发地震波的处理装置，其特征在于，所述最小相位化处理模块，包括：
地震波记录数据生成单元，记录可控震源激发的地震波，生成地震波记录数据；
零相位子波生成单元，对所述地震波记录数据进行自相关，生成可控震源零相位子波；
转换算子计算单元，利用所述可控震源零相位子波，计算将所述可控震源零相位子波转化为最小相位子波的转换算子；
最小相位子波转换单元，利用所述转换算子将所述可控震源零相位子波转化为最小相位子波，生成最小相位子波的地震波记录数据。
10.根据权利要求8所述的混源激发地震波的处理装置，其特征在于，所述振幅级别统一模块，包括：
第一补偿单元，按照所述可控震源与所述炸药震源的能量差异级别以及两者能量值的倍数，将非一致性主体的能量级别补偿至与所述一致性主体相同的能量级别；
第二补偿单元，通过几何扩散补偿的方式补偿两个震源纵向上的能量衰减；
第三补偿单元，通过地表一致性振幅补偿的方式补偿两个震源横向上的能量差异，进而使所述可控震源与所述炸药震源的振幅级别一致。
11.根据权利要求8所述的混源激发地震波的处理装置，其特征在于，所述时差确定模块，包括：
时差求取单元，选取重叠段进行时差调查，采用网格化互相关的方式求取时差；
校正单元，将求取的时差插值成校正量库，进行时差校正，并通过叠加检查时差校正效果。
12.根据权利要求8所述的混源激发地震波的处理装置，其特征在于，所述频率一致性匹配模块，包括：
匹配滤波单元，以一致性主体为目标地震道，以非一致性主体为原始地震道，对所述原始地震道进行匹配滤波使其逼近所述目标地震道；
匹配滤波算子求解单元，应用托布里兹矩阵方程求解匹配滤波算子；
用所述匹配滤波算子对原始地震道进行滤波，进而使所述可控震源和炸药震源实现频率的一致性匹配。
13.根据权利要求8所述的混源激发地震波的处理装置，其特征在于，还包括：
叠加检查模块，对进行一致性匹配后的所述可控震源和所述炸药震源进行叠加检查。
14.根据权利要求13所述的混源激发地震波的处理装置，其特征在于，所述叠加检查模块，包括：
分选单元，将所述炸药震源和所述可控震源激发的地震道按照线、道以及偏移距进行分选，得到多个分选道集；
剖地震面形成单元，利用动校速度将共深度点的道集叠加成一道，形成剖地震面，进而根据形成的剖地震面；
叠加检查单元，根据所述剖地震面进行叠加检查操作。
15.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至7任一项所述混源激发地震波的处理方法的步骤。
16.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至7任一项所述混源激发地震波的处理方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010021168.6A CN113109865B (zh) | 2020-01-09 | 2020-01-09 | 混源激发地震波的处理方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010021168.6A CN113109865B (zh) | 2020-01-09 | 2020-01-09 | 混源激发地震波的处理方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113109865A true CN113109865A (zh) | 2021-07-13 |
| CN113109865B CN113109865B (zh) | 2024-07-26 |
Family
ID=76709918
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010021168.6A Active CN113109865B (zh) | 2020-01-09 | 2020-01-09 | 混源激发地震波的处理方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113109865B (zh) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101545981A (zh) * | 2008-03-28 | 2009-09-30 | 中国石油天然气集团公司 | 可控震源地震数据零相位子波最小相位化方法 |
| CN101556338A (zh) * | 2008-04-10 | 2009-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | 一种可控震源自适应地表一致性反褶积方法 |
| US20100157733A1 (en) * | 2008-12-23 | 2010-06-24 | Mark Daly | Compensating seismic data for source variations |
| CN102798887A (zh) * | 2011-05-27 | 2012-11-28 | 中国石油天然气集团公司 | 一种提高可控震源记录初至质量的滤波方法 |
| CN103592683A (zh) * | 2013-10-17 | 2014-02-19 | 中国石油天然气集团公司 | 一种用于转换波地震资料的地表一致性振幅补偿方法及装置 |
| CN104536046A (zh) * | 2014-12-26 | 2015-04-22 | 吉林大学 | 基于地震记录的震源激发信号一致性评价方法 |
| CN105319587A (zh) * | 2014-07-04 | 2016-02-10 | 中国石油化工股份有限公司 | 一种基于初至子波提取的可控震源整形滤波方法 |
| WO2016141598A1 (zh) * | 2015-03-12 | 2016-09-15 | 中国矿业大学(北京) | 共姿态道集地表一致性振幅补偿方法及装置 |
| CN110261911A (zh) * | 2019-06-27 | 2019-09-20 | 中国石油化工股份有限公司 | 地震数据时差、相位差自动识别及校正方法 |
-
2020
- 2020-01-09 CN CN202010021168.6A patent/CN113109865B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101545981A (zh) * | 2008-03-28 | 2009-09-30 | 中国石油天然气集团公司 | 可控震源地震数据零相位子波最小相位化方法 |
| CN101556338A (zh) * | 2008-04-10 | 2009-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | 一种可控震源自适应地表一致性反褶积方法 |
| US20100157733A1 (en) * | 2008-12-23 | 2010-06-24 | Mark Daly | Compensating seismic data for source variations |
| CN102798887A (zh) * | 2011-05-27 | 2012-11-28 | 中国石油天然气集团公司 | 一种提高可控震源记录初至质量的滤波方法 |
| CN103592683A (zh) * | 2013-10-17 | 2014-02-19 | 中国石油天然气集团公司 | 一种用于转换波地震资料的地表一致性振幅补偿方法及装置 |
| CN105319587A (zh) * | 2014-07-04 | 2016-02-10 | 中国石油化工股份有限公司 | 一种基于初至子波提取的可控震源整形滤波方法 |
| CN104536046A (zh) * | 2014-12-26 | 2015-04-22 | 吉林大学 | 基于地震记录的震源激发信号一致性评价方法 |
| WO2016141598A1 (zh) * | 2015-03-12 | 2016-09-15 | 中国矿业大学(北京) | 共姿态道集地表一致性振幅补偿方法及装置 |
| CN110261911A (zh) * | 2019-06-27 | 2019-09-20 | 中国石油化工股份有限公司 | 地震数据时差、相位差自动识别及校正方法 |
Non-Patent Citations (10)
| Title |
|---|
| 向晓丽 等: "《混合震源交替激发地震资料处理方法》", 《物探与化探》, vol. 34, no. 6, pages 727 - 731 * |
| 曾维望 等: "《《山地复杂障碍区下煤层准确成像混合震源采集技术研究与实践》", 《中国煤炭地质》, no. 9, 30 September 2018 (2018-09-30), pages 76 - 81 * |
| 朱伟强 等: "《炸药震源和可控震源资料相位差异校正方法》", 《中国石油大学学报(自然科学版)》, vol. 30, no. 6, pages 31 - 35 * |
| 李建民: "《复杂地表条件下实施3D地震采集的基本对策――以涪陵义和3D采集项目为例》", 《江汉石油职工大学学报》, vol. 29, no. 6, 30 November 2016 (2016-11-30), pages 7 - 10 * |
| 杜娟: "《塔里木盆地二维地震老资料拼接处理技术及应用》", 《内蒙古石油化工》, no. 2, 28 February 2018 (2018-02-28), pages 92 - 98 * |
| 梁鸿贤: "《基于初至信息的可控震源和炸药震源地震资料匹配滤波技术》", 《石油物探》, vol. 57, no. 2, pages 231 - 236 * |
| 王敏: "《匹配滤波技术在江汉平原海相大剖面拼接处理中的应用》", 《工程地球物理学报》, vol. 14, no. 2, 31 March 2017 (2017-03-31), pages 165 - 173 * |
| 王翠华 等: "《相对振幅保持在连片处理中的应用》", 《天然气工业》, vol. 27, pages 228 - 231 * |
| 苏贵仕 等: "《地震数据处理中关于地震子波相位特性的探讨》", 《石油地球物理勘探》, vol. 43, no. 2, pages 121 - 122 * |
| 陈志卿 等: "《子波一致性校正方法在地震资料处理中的应用――以伊朗Kashan区块为例》", 《石油物探》, vol. 45, no. 2, pages 169 - 173 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113109865B (zh) | 2024-07-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Sun et al. | A convolutional neural network approach to deblending seismic data | |
| US9632192B2 (en) | Method of processing seismic data by providing surface offset common image gathers | |
| US4259733A (en) | Multi-dimensional seismic imaging | |
| US10557956B2 (en) | Method and system of processing seismic data by providing surface aperture common image gathers | |
| US9625593B2 (en) | Seismic data processing | |
| Bakulin et al. | Nonlinear beamforming for enhancement of 3D prestack land seismic data | |
| CN109633742B (zh) | 一种全波形反演方法及装置 | |
| CN108828668B (zh) | 一种叠前时间偏移数据处理方法及装置 | |
| CN110907995B (zh) | 井中vsp地震数据的逆时偏移方法及装置 | |
| CN111781635B (zh) | 海底四分量弹性波高斯束深度偏移方法和装置 | |
| CN113391351B (zh) | 一种基于被动源地震波场分析提取矿集区结构的方法 | |
| Li et al. | Finite difference seismic forward modeling method for fluid–solid coupled media with irregular seabed interface | |
| CN113109865B (zh) | 混源激发地震波的处理方法及装置 | |
| CN113743193B (zh) | 一种叠前地震资料线性干扰压制方法及系统 | |
| Corradini et al. | Investigating the influence of earthquake source complexity on back-projection images using convolutional neural networks | |
| CN112180446B (zh) | 三维地震数据叠前道集叠加方法及装置 | |
| CN114427452B (zh) | 微构造地质体的成像方法、装置、存储介质和计算机设备 | |
| CN103852786B (zh) | 一种应用于陆上地震数据的逆时偏移成像的方法及系统 | |
| CN113031071B (zh) | 地震波长波长静校正方法及装置 | |
| US20240319393A1 (en) | Systems and methods for seismic image generation using deep learning networks | |
| CN113126153B (zh) | 基于数据组合的叠前深度偏移方法及装置 | |
| CN114002741B (zh) | 叠前深度偏移方法及装置、计算机可读存储介质 | |
| CN112180448A (zh) | 剩余静校正方法、装置和电子设备 | |
| CN114063165B (zh) | 一种三维地震数据的拼接方法及装置 | |
| Mammadov et al. | Enhancement of Seismic Image Based on Novel Approach with Custom U-Net Architecture |
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 |