CN112379428A - 地震数据一致性处理方法及装置 - Google Patents
地震数据一致性处理方法及装置 Download PDFInfo
- Publication number
- CN112379428A CN112379428A CN202011202155.5A CN202011202155A CN112379428A CN 112379428 A CN112379428 A CN 112379428A CN 202011202155 A CN202011202155 A CN 202011202155A CN 112379428 A CN112379428 A CN 112379428A
- Authority
- CN
- China
- Prior art keywords
- quality
- low
- cannon
- section
- quality single
- 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
- 238000003672 processing method Methods 0.000 title claims description 25
- 238000012545 processing Methods 0.000 claims abstract description 94
- 238000012216 screening Methods 0.000 claims abstract description 43
- 238000000034 method Methods 0.000 claims abstract description 34
- 238000004590 computer program Methods 0.000 claims description 19
- 238000001914 filtration Methods 0.000 claims description 14
- 230000006872 improvement Effects 0.000 abstract description 16
- 230000000694 effects Effects 0.000 abstract description 12
- 238000010586 diagram Methods 0.000 description 23
- 230000005284 excitation Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 4
- 238000010276 construction Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 235000019738 Limestone Nutrition 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000006028 limestone Substances 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
- 230000010363 phase shift Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 238000006467 substitution reaction Methods 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
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
技术领域
本发明涉及石油地震资料处理技术领域，尤其涉及地震数据一致性处理方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
随着我国油气勘探进入大规模三维采集阶段，地震数据因工区的面积过大，出现了地表激发岩性差异大的问题。地震数据在进行一致性处理时，由低频干扰引起的横向特征不一致的现象比较严重。其主因在于地震数据的岩性激发品质差异大，且地震数据中相当数量的单炮记录，存在着频率低和信噪比低的问题。所以传统的地震数据一致性处理方法，对地震数据的改善程度低，制约了后续处理环节，严重影响构造解释、储层反演和油气预测的工作。
长期以来，传统的地震数据一致性处理方法，常采用地表一致性反褶积技术来改善激发岩性、激发效应和接收效应等造成的一致性问题，但此种方法对低频干扰重，且对子波差异大的地震数据改善程度小，对地震数据的横向特征一致性的改善效果较差，难以满足目前油气勘探工作的需求。
发明内容
本发明实施例提供一种地震数据一致性处理方法，用以提高地震数据一致性处理的改善效果，该方法包括：
求取地震数据中每一单炮的子波半周期长度；
根据单炮的子波半周期长度，从地震数据中筛选出低品质单炮和高品质单炮；所述低品质单炮为子波半周期长度大于等于第一预设长度的单炮；所述高品质单炮为子波半周期长度小于等于第二预设长度的单炮；第一预设长度大于第二预设长度；
对低品质单炮和高品质单炮进行叠加处理，得到品质改善后的低品质单炮；
以品质改善后的低品质单炮替换地震数据中的原始低品质单炮，得到品质改善后的地震数据；
对品质改善后的地震数据进行地表一致性处理。
本发明实施例还提供一种地震数据一致性处理装置，用以提高地震数据一致性处理的改善效果，该装置包括：
子波半周期长度求取模块，用于求取地震数据中每一单炮的子波半周期长度；
单炮筛选模块，用于根据单炮的子波半周期长度，从地震数据中筛选出低品质单炮和高品质单炮；所述低品质单炮为子波半周期长度大于等于第一预设长度的单炮；所述高品质单炮为子波半周期长度小于等于第二预设长度的单炮；第一预设长度大于第二预设长度；
叠加处理模块，用于对低品质单炮和高品质单炮进行叠加处理，得到品质改善后的低品质单炮；
数据替换模块，用于以品质改善后的低品质单炮替换地震数据中的原始低品质单炮，得到品质改善后的地震数据；
地表一致性处理模块，用于对品质改善后的地震数据进行地表一致性处理。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述地震数据一致性处理方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述地震数据一致性处理方法的计算机程序。
本发明实施例中，求取地震数据中每一单炮的子波半周期长度；根据单炮的子波半周期长度，从地震数据中筛选出低品质单炮和高品质单炮；所述低品质单炮为子波半周期长度大于等于第一预设长度的单炮；所述高品质单炮为子波半周期长度小于等于第二预设长度的单炮；第一预设长度大于第二预设长度；对低品质单炮和高品质单炮进行叠加处理，得到品质改善后的低品质单炮；以品质改善后的低品质单炮替换地震数据中的原始低品质单炮，得到品质改善后的地震数据；对品质改善后的地震数据进行地表一致性处理，从而可在地表一致性处理之前对地震数据进行品质改善，与现有技术对比，避免了直接对地震数据进行地表一致性处理时易出现的改善程度小的问题，提高了地震数据一致性处理的改善效果，可满足目前油气勘探工作的需求。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中一种地震数据一致性处理方法的流程示意图；
图2为本发明实施例中一种地震数据一致性处理方法的具体示例图；
图3为本发明实施例中一种地震数据一致性处理方法的具体示例图；
图4为本发明实施例中一种地震数据一致性处理方法的具体示例图；
图5为本发明实施例中一种地震数据一致性处理方法的具体示例图；
图6为本发明实施例中一种地震数据一致性处理方法的具体示例图；
图7为本发明实施例中一种地震数据一致性处理方法的具体示例图；
图8为本发明实施例中一种地震数据一致性处理方法的具体示例图；
图9为本发明实施例中一种地震数据一致性处理装置的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
图1为本发明实施例中一种地震数据一致性处理方法的流程示意图，如图1所示，本发明实施例提供的一种地震数据一致性处理方法，可以包括：
步骤101：求取地震数据中每一单炮的子波半周期长度；
步骤102：根据单炮的子波半周期长度，从地震数据中筛选出低品质单炮和高品质单炮；所述低品质单炮为子波半周期长度大于等于第一预设长度的单炮；所述高品质单炮为子波半周期长度小于等于第二预设长度的单炮；第一预设长度大于第二预设长度；
步骤103：对低品质单炮和高品质单炮进行叠加处理，得到品质改善后的低品质单炮；
步骤104：以品质改善后的低品质单炮替换地震数据中的原始低品质单炮，得到品质改善后的地震数据；
步骤105：对品质改善后的地震数据进行地表一致性处理。
本发明实施例中，求取地震数据中每一单炮的子波半周期长度；根据单炮的子波半周期长度，从地震数据中筛选出低品质单炮和高品质单炮；所述低品质单炮为子波半周期长度大于等于第一预设长度的单炮；所述高品质单炮为子波半周期长度小于等于第二预设长度的单炮；第一预设长度大于第二预设长度；对低品质单炮和高品质单炮进行叠加处理，得到品质改善后的低品质单炮；以品质改善后的低品质单炮替换地震数据中的原始低品质单炮，得到品质改善后的地震数据；对品质改善后的地震数据进行地表一致性处理，从而可在地表一致性处理之前对地震数据进行品质改善，与现有技术对比，避免了直接对地震数据进行地表一致性处理时易出现的改善程度小的问题，提高了地震数据一致性处理的改善效果，可满足目前油气勘探工作的需求。
具体实施时，首先求取地震数据中每一单炮的子波半周期长度。
实施例中，求取地震数据中每一单炮的子波半周期长度，可以包括：提取工区范围内所有单炮的子波半周期长度。
在上述实施例中，通过求取地震数据中每一单炮的子波半周期长度，有利于在后续步骤中对地震数据的品质进行改善。
具体实施时，在求取地震数据中每一单炮的子波半周期长度后，根据单炮的子波半周期长度，从地震数据中筛选出低品质单炮和高品质单炮；所述低品质单炮为子波半周期长度大于等于第一预设长度的单炮；所述高品质单炮为子波半周期长度小于等于第二预设长度的单炮；第一预设长度大于第二预设长度。
实施例中，可通过设计的程序自动搜索地震数据处理范围内的低品质单炮和高品质单炮。
在上述实施例中，可通过设计程序自动实现从地震数据中筛选出低品质单炮和高品质单炮，提高了工作效率。
具体实施时，根据单炮的子波半周期长度，从地震数据中筛选出低品质单炮的方法有多种，例如，可以包括：以子波半周期长度从小到大的顺序，提取预设数量的单炮；所述预设数量根据预设的占单炮总数的比例和单炮总数确定；计算提取的预设数量的单炮的子波半周期长度的平均值；将所述平均值的预设倍数作为第一预设长度，筛选出子波半周期长度大于第一预设长度，且总数小于所述预设数量的低品质单炮。
实施例中，首先对单炮以子波半周期长度从小到大的顺序进行排列；之后，提取预设数量的单炮，其中，预设数量可为占单炮总数的前百分之十的数量；其次，计算提取的该部分单炮的子波半周期长度的平均值；设定低品质单炮为子波半周期大于等于所述平均值的预设倍数，如2倍，且低品质单炮数量小于等于工区范围内所有单炮总炮数的百分之十，筛选出低品质单炮。
在上述实施例中，通过将所述平均值的预设倍数作为第一预设长度，筛选出子波半周期长度大于第一预设长度，且总数小于所述预设数量的低品质单炮，可有效地筛选出地震数据中的低品质单炮；通过控制低品质单炮的数量小于等于预设数量，可保证不破坏后续的地震数据地表一致性处理。
具体实施时，根据单炮的子波半周期长度，从地震数据中筛选出高品质单炮的方法有多种，例如，可以包括：对筛选出的低品质单炮进行叠加处理，得到工区的叠加剖面体；以预设的信噪比，从工区的叠加剖面体中抽取叠加剖面；计算从工区的叠加剖面体中抽取的叠加剖面的覆盖次数；以子波半周期长度从小到大的顺序，计算每一单炮的叠加剖面体；所述每一单炮的叠加剖面体根据该单炮以及子波半周期长度小于等于该单炮子波半周期长度的单炮，按顺序进行叠加处理得到；从每一单炮的叠加剖面体中，分别抽取每一单炮的叠加剖面，每一单炮的叠加剖面与从工区的叠加剖面体中抽取的叠加剖面地理位置相同；计算每一单炮的叠加剖面的覆盖次数；根据从工区的叠加剖面体中抽取的叠加剖面的覆盖次数，和每一单炮的叠加剖面的覆盖次数，从地震数据中筛选出高品质单炮。
实施例中，根据从工区的叠加剖面体中抽取的叠加剖面的覆盖次数，和每一单炮的叠加剖面的覆盖次数，从地震数据中筛选出高品质单炮的方法有多种，例如，可以包括：将每一单炮的叠加剖面的覆盖次数，与从工区的叠加剖面体中抽取的叠加剖面的覆盖次数进行匹配，得到每一单炮的匹配程度；所述匹配程度用于表征每一单炮的叠加剖面的覆盖次数，与从工区的叠加剖面体中抽取的叠加剖面的覆盖次数的比值；根据匹配程度最高的单炮的叠加剖面的覆盖次数，从地震数据中筛选出高品质单炮。
举一例，按以下步骤对地震数据的道头进行计算，可以从地震数据中筛选出高品质单炮：
1、对全工区地震数据中筛选出的低品质单炮进行叠加，得到全工区的叠加剖面体a；
2、从a中抽取一条叠加剖面A，并计算该叠加剖面A的覆盖次数，其中，覆盖次数等于叠加剖面A中参与进行叠加的地震道数；
3、以子波半周期长度从小到大的顺序，计算每一单炮的叠加剖面体，可以包括：对全工区的单炮，按子波半周期长度从小到大的顺序，根据该单炮以及子波半周期长度小于等于该单炮子波半周期长度的单炮，按顺序进行叠加处理，得到每一单炮的叠加剖面体，如b1、b2、b3…bn，bn表示子波半周期长度为nms的单炮的叠加剖面体；
例如，①、对于子波半周期长度为1毫秒的单炮，计算该单炮的叠加剖面体的方式为对自身进行叠加处理；②、对于子波半周期长度为2毫秒的单炮，计算该单炮的叠加剖面体的方式为将子波半周期长度为1毫秒的单炮和子波半周期长度为2毫秒的单炮进行叠加处理；③、对于子波半周期长度为n毫秒的单炮，计算该单炮的叠加剖面体的方式为将子波半周期长度分别1ms、2ms、…、nms的单炮进行叠加处理；
4、在b1、b2、b3…bn上分别抽取与该叠加剖面A相同地理位置的叠加剖面，记为B1、B2、B3…Bn，并计算B1、B2、B3…Bn的覆盖次数；
5、从B1、B2、B3…Bn中找出与叠加剖面A覆盖次数相同或者最接近的叠加剖面Bn，从而找到对应的高品质单炮，从地震数据中筛选出高品质单炮。如，在确定子波半周期长度为nms的单炮的匹配程度最高时，则将子波半周期长度为1ms、2ms、…、nms的单炮作为高品质单炮。
在上述实施例中，通过叠加处理，可自动从地震数据中筛选出高品质单炮，提高了运算效率。
具体实施时，在根据单炮的子波半周期长度，从地震数据中筛选出低品质单炮和高品质单炮后，对低品质单炮和高品质单炮进行叠加处理，得到品质改善后的低品质单炮。
实施例中，对低品质单炮和高品质单炮进行叠加处理，得到品质改善后的低品质单炮数据的方法有多种，例如，可以包括：对低品质单炮和高品质单炮进行叠加处理，分别得到低品质单炮的地震剖面数据和高品质单炮的地震剖面数据；从低品质单炮的地震剖面数据和高品质单炮的地震剖面数据中，抽取与从工区的叠加剖面体中抽取的叠加剖面地理位置相同的叠加剖面，分别得到低品质单炮的地震剖面和高品质单炮的地震剖面；将低品质单炮的地震剖面和高品质单炮的地震剖面进行互相关处理，得到滤波算子；所述滤波算子用于改善低品质单炮的品质；将滤波算子与低品质单炮进行褶积处理，得到品质改善后的低品质单炮。
举一例，将得到的低品质单炮和高品质单炮分别进行叠加，可分别得到低品质单炮的地震剖面数据d和高品质单炮的地震剖面数据地震剖面数据g，并分别抽取与从工区的叠加剖面体中抽取的叠加剖面地理位置相同的叠加剖面d1和g1，之后将叠加剖面g1、d1进行互相关处理，得到滤波算子。
在上述实施例中，在滤波算子的求取过程中，应尽量选择水平方向的叠加道作为目标，提取接近零相位的滤波算子，避免匹配处理后的单炮产生较大的相位移动，从而影响后续一致性处理。
具体实施时，在对低品质单炮和高品质单炮进行叠加处理，得到品质改善后的低品质单炮后，以品质改善后的低品质单炮替换地震数据中的原始低品质单炮，得到品质改善后的地震数据；对品质改善后的地震数据进行地表一致性处理。
实施例中，在将品质改善后的低品质单炮，替代地震数据中的原始低品质单炮，可得到品质改善后的地震数据，再进行常规的地表一致性处理。
在上述实施例中，对比现有技术中，在对地震数据进行处理时，未对高品质单炮和低品质单炮进行识别筛选，因此无论用匹配滤波还是子波整形的方法去消除低频干扰带来的横向特征差异，效果都不理想，此外，这样做会给高品质的单炮记录带来不良影响，破坏地震数据原有的地球物理特征。本发明实施例在对低品质单炮和高品质单炮进行叠加处理，得到品质改善后的低品质单炮后，以品质改善后的低品质单炮替换地震数据中的原始低品质单炮，得到品质改善后的地震数据；对品质改善后的地震数据进行地表一致性处理，仅对低品质单炮进行了处理，可以较好地解决上述问题。
下面给出一个具体实施例，来说明本发明的方法的具体应用。
图2为本发明实施例中一种地震数据一致性处理方法的具体示例图，如图2所示，本发明实施例提供的地震数据一致性处理方法的具体示例，可以包括：
(1)、首先加载地震数据：从地震处理软件中提取工区范围内的单炮子波，并求取各单炮子波的半周期长度；可按照炮点桩号、子波半周期长度的格式规范整理，将其放入预设置的执行上述地震数据一致性处理方法的计算机程序的计算路径下，可直接得出工区地震数据的单炮位置和子波半周期长度的统计结果；举一例，如可直接得出工区地震数据的单炮位置图和子波半周期长度的统计结果图，如图3展示工区地震数据的单炮位置图，图3采用BJ54坐标，图纵轴为X坐标，横轴为Y坐标，图中黑点表示工区内所有单炮位置；图4展示子波半周期长度的统计图，图4中左轴表示工区内该单炮的子波半周期长度，单位为毫秒；右轴表示该类型单炮在工区内的百分比；
(2)、运行预设置的执行上述地震数据一致性处理方法的计算机程序中的低品质单炮自动筛选程序(运行程序2)，求取低品质单炮的子波半周期长度平均值，实现对低品质单炮的筛选。如图5所示，图5中展示了提取的低品质单炮位置图与地表地质图，图5中筛选出的低品质单炮的位置轮廓范围与地质图上灰岩区(低品质单炮激发区)一致，由此说明上述低品质单炮自动筛选程序的有效性；
(3)、运行预设置的执行上述地震数据一致性处理方法的计算机程序中的高品质单炮自动筛选程序(运行程序1)，根据工区范围内高品质与低品质单炮覆盖次数对等的原则，由子波半周期从小到大对高品质单炮进行筛选，求取高品质单炮子波半周期值范围，实现对高品质单炮的筛选，如图6所示，图6展示提取的高品质单炮的位置图，图6采用BJ54坐标，图纵轴为X坐标，横轴为Y坐标，图中黑点表示从提取的高品质单炮位置；
(4)、将高品质单炮和低品质单炮选排成共中心点道集(CMP)，分别进行CMP域叠加，选择重合叠加段；在给定的时窗内提取相关子波，求取滤波算子，如图7所示，图7是高品质单炮和低品质单炮在相同位置的叠加图与提取的滤波算子图，其中a部分表示频率异常的单炮叠加道，展示的是提取的高品质单炮进行叠加的剖面重合段；b部分表示高品质单炮的叠加道，展示的是提取的低品质单炮进行叠加的剖面重合段；c部分表示提取的滤波算子，为表示该高品质单炮剖面与低品质剖面重合段进行互相关提取的滤波算子；
(5)、将求取的滤波算子在共炮集上应用于低品质单炮，获得品质改善后的低品质单炮；
(6)、将改善后的低品质单炮代替对应的低品质单炮，再进行常规的地表一致性处理；
(7)、如图8所示，在用户对种地震数据一致性处理的结果满意时，进行结果及图件输出，图8匹配(匹配是指上述的“以品质改善后的低品质单炮替换地震数据中的原始低品质单炮，得到品质改善后的地震数据”步骤)前后叠加剖面及低频干扰段频谱，图8展示了匹配前剖面上存在明显的低频干扰，造成了整个剖面特征不一致，匹配后剖面特征趋于一致，频谱拓宽，分辨率明显提高。
本发明实施例还提供了一种应用上述地震数据一致性处理方法的另一实例，在该实例中，在对四川盆地长宁、威远地区的页岩气连片开发处理中应用了上述方案，可使低信噪比段地震数据的横向一致性得到较好地改善，纵向分辨率提高了约30％，同时为后续构造解释、储层反演和油气预测奠定了良好的基础。
当然，可以理解的是，上述详细流程还可以有其他变化例，相关变化例均应落入本发明的保护范围。
本发明实施例中，求取地震数据中每一单炮的子波半周期长度；根据单炮的子波半周期长度，从地震数据中筛选出低品质单炮和高品质单炮；所述低品质单炮为子波半周期长度大于等于第一预设长度的单炮；所述高品质单炮为子波半周期长度小于等于第二预设长度的单炮；第一预设长度大于第二预设长度；对低品质单炮和高品质单炮进行叠加处理，得到品质改善后的低品质单炮；以品质改善后的低品质单炮替换地震数据中的原始低品质单炮，得到品质改善后的地震数据；对品质改善后的地震数据进行地表一致性处理，从而可在地表一致性处理之前对地震数据进行品质改善，与现有技术对比，避免了直接对地震数据进行地表一致性处理时易出现的改善程度小的问题，提高了地震数据一致性处理的改善效果，可满足目前油气勘探工作的需求。
如上所述，为解决复杂地表的激发条件导致的剖面特征问题，本发明实施例提供了一种可应用于复杂激发岩性区域的地震数据一致性处理方法，可用于改善叠加剖面横向一致性和提高纵向分辨率。通过在叠前处理前实施上述地震数据一致性处理方法，可以有效的解决上述难题。上述地震数据一致性处理方法运行效率高、操作简便，具有进行项目生产的良好条件，同时丰富了目前较为单一的地震数据一致性处理手段，可以为石油地震资料处理领域加强技术储备。本发明实施例提出了分离高品质单炮和低品质单炮的思路来求取滤波算子的方法，并可通过设计的程序自动实现，可以较大的提高工作效率。
本发明实施例中还提供了一种地震数据一致性处理装置，如下面的实施例所述。由于该装置解决问题的原理与地震数据一致性处理方法相似，因此该装置的实施可以参见地震数据一致性处理方法的实施，重复之处不再赘述。
图9为本发明实施例中一种地震数据一致性处理装置的结构示意图，如图9所示，本发明实施例提供的地震数据一致性处理装置，可以包括：
子波半周期长度求取模块01，用于求取地震数据中每一单炮的子波半周期长度；
单炮筛选模块02，用于根据单炮的子波半周期长度，从地震数据中筛选出低品质单炮和高品质单炮；所述低品质单炮为子波半周期长度大于等于第一预设长度的单炮；所述高品质单炮为子波半周期长度小于等于第二预设长度的单炮；第一预设长度大于第二预设长度；
叠加处理模块03，用于对低品质单炮和高品质单炮进行叠加处理，得到品质改善后的低品质单炮；
数据替换模块04，用于以品质改善后的低品质单炮替换地震数据中的原始低品质单炮，得到品质改善后的地震数据；
地表一致性处理模块05，用于对品质改善后的地震数据进行地表一致性处理。
在一个实施例中，单炮筛选模块，具体用于：
以子波半周期长度从小到大的顺序，提取预设数量的单炮；所述预设数量根据预设的占单炮总数的比例和单炮总数确定；
计算提取的预设数量的单炮的子波半周期长度的平均值；
将所述平均值的预设倍数作为第一预设长度，筛选出子波半周期长度大于等于第一预设长度，且总数小于所述预设数量的低品质单炮。
在一个实施例中，单炮筛选模块，具体用于：
对筛选出的低品质单炮进行叠加处理，得到工区的叠加剖面体；
以预设的信噪比，从工区的叠加剖面体中抽取叠加剖面；
计算从工区的叠加剖面体中抽取的叠加剖面的覆盖次数；
以子波半周期长度从小到大的顺序，计算每一单炮的叠加剖面体；所述每一单炮的叠加剖面体根据该单炮以及子波半周期长度小于等于该单炮子波半周期长度的单炮，按顺序进行叠加处理得到；
从每一单炮的叠加剖面体中，分别抽取每一单炮的叠加剖面，每一单炮的叠加剖面与从工区的叠加剖面体中抽取的叠加剖面地理位置相同；
计算每一单炮的叠加剖面的覆盖次数；
根据从工区的叠加剖面体中抽取的叠加剖面的覆盖次数，和每一单炮的叠加剖面的覆盖次数，从地震数据中筛选出高品质单炮。
在一个实施例中，单炮筛选模块，具体用于：
将每一单炮的叠加剖面的覆盖次数，与从工区的叠加剖面体中抽取的叠加剖面的覆盖次数进行匹配，得到每一单炮的匹配程度；所述匹配程度用于表征每一单炮的叠加剖面的覆盖次数，与从工区的叠加剖面体中抽取的叠加剖面的覆盖次数的比值；
根据匹配程度最高的单炮的叠加剖面的覆盖次数，从地震数据中筛选出高品质单炮。
在一个实施例中，叠加处理模块，具体用于：
对低品质单炮和高品质单炮进行叠加处理，分别得到低品质单炮的地震剖面数据和高品质单炮的地震剖面数据；
从低品质单炮的地震剖面数据和高品质单炮的地震剖面数据中，抽取与从工区的叠加剖面体中抽取的叠加剖面地理位置相同的叠加剖面，分别得到低品质单炮的地震剖面和高品质单炮的地震剖面；
将低品质单炮的地震剖面和高品质单炮的地震剖面进行互相关处理，得到滤波算子；所述滤波算子用于改善低品质单炮的品质；
将滤波算子与低品质单炮进行褶积处理，得到品质改善后的低品质单炮。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述地震数据一致性处理方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述地震数据一致性处理方法的计算机程序。
本发明实施例中，求取地震数据中每一单炮的子波半周期长度；根据单炮的子波半周期长度，从地震数据中筛选出低品质单炮和高品质单炮；所述低品质单炮为子波半周期长度大于等于第一预设长度的单炮；所述高品质单炮为子波半周期长度小于等于第二预设长度的单炮；第一预设长度大于第二预设长度；对低品质单炮和高品质单炮进行叠加处理，得到品质改善后的低品质单炮；以品质改善后的低品质单炮替换地震数据中的原始低品质单炮，得到品质改善后的地震数据；对品质改善后的地震数据进行地表一致性处理，从而可在地表一致性处理之前对地震数据进行品质改善，与现有技术对比，避免了直接对地震数据进行地表一致性处理时易出现的改善程度小的问题，提高了地震数据一致性处理的改善效果，可满足目前油气勘探工作的需求。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (12)
1.一种地震数据一致性处理方法，其特征在于，包括：
求取地震数据中每一单炮的子波半周期长度；
根据单炮的子波半周期长度，从地震数据中筛选出低品质单炮和高品质单炮；所述低品质单炮为子波半周期长度大于等于第一预设长度的单炮；所述高品质单炮为子波半周期长度小于等于第二预设长度的单炮；第一预设长度大于第二预设长度；
对低品质单炮和高品质单炮进行叠加处理，得到品质改善后的低品质单炮；
以品质改善后的低品质单炮替换地震数据中的原始低品质单炮，得到品质改善后的地震数据；
对品质改善后的地震数据进行地表一致性处理。
2.如权利要求1所述的方法，其特征在于，根据单炮的子波半周期长度，从地震数据中筛选出低品质单炮，包括：
以子波半周期长度从小到大的顺序，提取预设数量的单炮；所述预设数量根据预设的占单炮总数的比例和单炮总数确定；
计算提取的预设数量的单炮的子波半周期长度的平均值；
将所述平均值的预设倍数作为第一预设长度，筛选出子波半周期长度大于等于第一预设长度，且总数小于所述预设数量的低品质单炮。
3.如权利要求1所述的方法，其特征在于，根据单炮的子波半周期长度，从地震数据中筛选出高品质单炮，包括：
对筛选出的低品质单炮进行叠加处理，得到工区的叠加剖面体；
以预设的信噪比，从工区的叠加剖面体中抽取叠加剖面；
计算从工区的叠加剖面体中抽取的叠加剖面的覆盖次数；
以子波半周期长度从小到大的顺序，计算每一单炮的叠加剖面体；所述每一单炮的叠加剖面体根据该单炮以及子波半周期长度小于等于该单炮子波半周期长度的单炮，按顺序进行叠加处理得到；
从每一单炮的叠加剖面体中，分别抽取每一单炮的叠加剖面，每一单炮的叠加剖面与从工区的叠加剖面体中抽取的叠加剖面地理位置相同；
计算每一单炮的叠加剖面的覆盖次数；
根据从工区的叠加剖面体中抽取的叠加剖面的覆盖次数，和每一单炮的叠加剖面的覆盖次数，从地震数据中筛选出高品质单炮。
4.如权利要求3所述的方法，其特征在于，根据从工区的叠加剖面体中抽取的叠加剖面的覆盖次数，和每一单炮的叠加剖面的覆盖次数，从地震数据中筛选出高品质单炮，包括：
将每一单炮的叠加剖面的覆盖次数，与从工区的叠加剖面体中抽取的叠加剖面的覆盖次数进行匹配，得到每一单炮的匹配程度；所述匹配程度用于表征每一单炮的叠加剖面的覆盖次数，与从工区的叠加剖面体中抽取的叠加剖面的覆盖次数的比值；
根据匹配程度最高的单炮的叠加剖面的覆盖次数，从地震数据中筛选出高品质单炮。
5.如权利要求3所述的方法，其特征在于，对低品质单炮和高品质单炮进行叠加处理，得到品质改善后的低品质单炮数据，包括：
对低品质单炮和高品质单炮进行叠加处理，分别得到低品质单炮的地震剖面数据和高品质单炮的地震剖面数据；
从低品质单炮的地震剖面数据和高品质单炮的地震剖面数据中，抽取与从工区的叠加剖面体中抽取的叠加剖面地理位置相同的叠加剖面，分别得到低品质单炮的地震剖面和高品质单炮的地震剖面；
将低品质单炮的地震剖面和高品质单炮的地震剖面进行互相关处理，得到滤波算子；所述滤波算子用于改善低品质单炮的品质；
将滤波算子与低品质单炮进行褶积处理，得到品质改善后的低品质单炮。
6.一种地震数据一致性处理装置，其特征在于，包括：
子波半周期长度求取模块，用于求取地震数据中每一单炮的子波半周期长度；
单炮筛选模块，用于根据单炮的子波半周期长度，从地震数据中筛选出低品质单炮和高品质单炮；所述低品质单炮为子波半周期长度大于等于第一预设长度的单炮；所述高品质单炮为子波半周期长度小于等于第二预设长度的单炮；第一预设长度大于第二预设长度；
叠加处理模块，用于对低品质单炮和高品质单炮进行叠加处理，得到品质改善后的低品质单炮；
数据替换模块，用于以品质改善后的低品质单炮替换地震数据中的原始低品质单炮，得到品质改善后的地震数据；
地表一致性处理模块，用于对品质改善后的地震数据进行地表一致性处理。
7.如权利要求6所述的装置，其特征在于，单炮筛选模块，具体用于：
以子波半周期长度从小到大的顺序，提取预设数量的单炮；所述预设数量根据预设的占单炮总数的比例和单炮总数确定；
计算提取的预设数量的单炮的子波半周期长度的平均值；
将所述平均值的预设倍数作为第一预设长度，筛选出子波半周期长度大于等于第一预设长度，且总数小于所述预设数量的低品质单炮。
8.如权利要求6所述的装置，其特征在于，单炮筛选模块，具体用于：
对筛选出的低品质单炮进行叠加处理，得到工区的叠加剖面体；
以预设的信噪比，从工区的叠加剖面体中抽取叠加剖面；
计算从工区的叠加剖面体中抽取的叠加剖面的覆盖次数；
以子波半周期长度从小到大的顺序，计算每一单炮的叠加剖面体；所述每一单炮的叠加剖面体根据该单炮以及子波半周期长度小于等于该单炮子波半周期长度的单炮，按顺序进行叠加处理得到；
从每一单炮的叠加剖面体中，分别抽取每一单炮的叠加剖面，每一单炮的叠加剖面与从工区的叠加剖面体中抽取的叠加剖面地理位置相同；
计算每一单炮的叠加剖面的覆盖次数；
根据从工区的叠加剖面体中抽取的叠加剖面的覆盖次数，和每一单炮的叠加剖面的覆盖次数，从地震数据中筛选出高品质单炮。
9.如权利要求8所述的装置，其特征在于，单炮筛选模块，具体用于：
将每一单炮的叠加剖面的覆盖次数，与从工区的叠加剖面体中抽取的叠加剖面的覆盖次数进行匹配，得到每一单炮的匹配程度；所述匹配程度用于表征每一单炮的叠加剖面的覆盖次数，与从工区的叠加剖面体中抽取的叠加剖面的覆盖次数的比值；
根据匹配程度最高的单炮的叠加剖面的覆盖次数，从地震数据中筛选出高品质单炮。
10.如权利要求8所述的装置，其特征在于，叠加处理模块，具体用于：
对低品质单炮和高品质单炮进行叠加处理，分别得到低品质单炮的地震剖面数据和高品质单炮的地震剖面数据；
从低品质单炮的地震剖面数据和高品质单炮的地震剖面数据中，抽取与从工区的叠加剖面体中抽取的叠加剖面地理位置相同的叠加剖面，分别得到低品质单炮的地震剖面和高品质单炮的地震剖面；
将低品质单炮的地震剖面和高品质单炮的地震剖面进行互相关处理，得到滤波算子；所述滤波算子用于改善低品质单炮的品质；
将滤波算子与低品质单炮进行褶积处理，得到品质改善后的低品质单炮。
11.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5任一所述方法。
12.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至5任一所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011202155.5A CN112379428A (zh) | 2020-11-02 | 2020-11-02 | 地震数据一致性处理方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011202155.5A CN112379428A (zh) | 2020-11-02 | 2020-11-02 | 地震数据一致性处理方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112379428A true CN112379428A (zh) | 2021-02-19 |
Family
ID=74576101
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011202155.5A Pending CN112379428A (zh) | 2020-11-02 | 2020-11-02 | 地震数据一致性处理方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112379428A (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010014118A1 (en) * | 2008-08-01 | 2010-02-04 | Ikelle Luc T | Statistical decoding and imaging of multishot and single-shot seismic data |
| CN101900833A (zh) * | 2010-06-02 | 2010-12-01 | 西安石油大学 | 一种地震散射p-p波成像速度分析方法 |
| CN102590858A (zh) * | 2011-12-31 | 2012-07-18 | 中国石油集团西北地质研究所 | 基于宽频子波重构的双程波成像方法 |
| CN103777239A (zh) * | 2012-10-17 | 2014-05-07 | 中国石油化工股份有限公司 | 一种提高不同震源资料相位一致性的相位旋转方法 |
| US20180246241A1 (en) * | 2016-02-02 | 2018-08-30 | China University Of Mining & Technology, Beijing | Method and apparatus for identifying collapsed coal column |
| CN109541677A (zh) * | 2018-11-15 | 2019-03-29 | 中石化石油工程技术服务有限公司 | 一种保护激发频率的延时震源组合参数计算方法 |
| CN109884705A (zh) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | 双重约束时频域子波谱提高地震分辨率处理方法 |
-
2020
- 2020-11-02 CN CN202011202155.5A patent/CN112379428A/zh active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010014118A1 (en) * | 2008-08-01 | 2010-02-04 | Ikelle Luc T | Statistical decoding and imaging of multishot and single-shot seismic data |
| CN101900833A (zh) * | 2010-06-02 | 2010-12-01 | 西安石油大学 | 一种地震散射p-p波成像速度分析方法 |
| CN102590858A (zh) * | 2011-12-31 | 2012-07-18 | 中国石油集团西北地质研究所 | 基于宽频子波重构的双程波成像方法 |
| CN103777239A (zh) * | 2012-10-17 | 2014-05-07 | 中国石油化工股份有限公司 | 一种提高不同震源资料相位一致性的相位旋转方法 |
| US20180246241A1 (en) * | 2016-02-02 | 2018-08-30 | China University Of Mining & Technology, Beijing | Method and apparatus for identifying collapsed coal column |
| CN109541677A (zh) * | 2018-11-15 | 2019-03-29 | 中石化石油工程技术服务有限公司 | 一种保护激发频率的延时震源组合参数计算方法 |
| CN109884705A (zh) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | 双重约束时频域子波谱提高地震分辨率处理方法 |
Non-Patent Citations (6)
| Title |
|---|
| BERNARD RAYNAUD: "《Diffraction modelling of 3-D lower-crustal reflectors 》", GEOPHYSICAL JOURNAL, 31 December 1988 (1988-12-31), pages 149 - 161 * |
| 王功军等: "《波阻抗反演中地震子波的处理》", 长江大学学报(自然科学版), vol. 8, no. 2, 28 February 2011 (2011-02-28), pages 35 - 37 * |
| 翟桐立等: "《全方位高密度单点接收地震采集技术》", 石油学报, vol. 37, no. 2, 31 December 2016 (2016-12-31), pages 56 - 63 * |
| 赵桂玲等: "《叠前偏移能量一致性处理方法及应用》", 地球物理学进展, vol. 31, no. 3, 30 June 2016 (2016-06-30), pages 1295 - 1299 * |
| 龙浩等: "《一种新的子波匹配技术在复杂岩性激发地震数据处理中的应用——以四川盆地CN区块为例》", 2019油气田勘探与开发国际会议论文集, 31 October 2019 (2019-10-31), pages 2 - 4 * |
| 龙浩等: "一种新的子波匹配技术在复杂岩性激发地震数据处理中的应用——以四川盆地CN区块为例", 2019油气田勘探与开发国际会议论文集, 31 October 2019 (2019-10-31), pages 2 - 4 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105353407B (zh) | 一种叠后地震波阻抗反演方法 | |
| US7917295B2 (en) | Modeling and filtering coherent noise in seismic surveying | |
| CN110907995B (zh) | 井中vsp地震数据的逆时偏移方法及装置 | |
| CN109633742B (zh) | 一种全波形反演方法及装置 | |
| CN100561256C (zh) | 一种避免动校拉伸的动校正叠加方法 | |
| CN107728206A (zh) | 一种速度场建模方法 | |
| CN111983683A (zh) | 一种少井条件下盐下湖相灰岩储层预测方法和系统 | |
| CN113534259B (zh) | 一种可控震源高效采集实时叠前时间偏移成像方法 | |
| CN105911590A (zh) | 一种谐波噪声的压制处理方法及装置 | |
| CN112379428A (zh) | 地震数据一致性处理方法及装置 | |
| CN112198547A (zh) | 深层或超深层地震资料处理方法及装置 | |
| CN113391353A (zh) | 地震数据的处理方法及装置 | |
| CN112946742B (zh) | 一种拾取精确叠加速度谱的方法 | |
| CN105277981A (zh) | 基于波场延拓补偿的非一致性时移地震面元匹配方法 | |
| CN112882101B (zh) | 叠前地震数据的随机噪音衰减方法及装置 | |
| CN112462425A (zh) | 海底节点混采数据中交叉干扰源的识别方法及装置 | |
| CN111929726B (zh) | 地震相干数据体处理方法及装置 | |
| CN111965703B (zh) | 井中地震数据中单炮记录的下行波反褶积处理方法及系统 | |
| CN111965706B (zh) | 地震反演方法及装置 | |
| CN112799132B (zh) | 微局部线性噪声压制方法及装置 | |
| CN113138419B (zh) | 提取下行子波和衰减参数的方法、装置 | |
| CN113589384B (zh) | 基于信号随偏移距变化特征的叠前道集保幅去噪方法 | |
| CN112255689B (zh) | 多次波发育区地震数据保真速度分析方法 | |
| CN113835123B (zh) | 基于地质目标叠前偏移成像的地震采集参数分析方法 | |
| CN113917539A (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 |