CN114063148B - 一种基于贝叶斯判别的折射波初至优选方法及系统 - Google Patents

一种基于贝叶斯判别的折射波初至优选方法及系统 Download PDF

Info

Publication number
CN114063148B
CN114063148B CN202010756674.XA CN202010756674A CN114063148B CN 114063148 B CN114063148 B CN 114063148B CN 202010756674 A CN202010756674 A CN 202010756674A CN 114063148 B CN114063148 B CN 114063148B
Authority
CN
China
Prior art keywords
arrival
wave
judging
sample
refraction wave
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
Application number
CN202010756674.XA
Other languages
English (en)
Other versions
CN114063148A (zh
Inventor
李艳东
宋建勇
李文科
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202010756674.XA priority Critical patent/CN114063148B/zh
Publication of CN114063148A publication Critical patent/CN114063148A/zh
Application granted granted Critical
Publication of CN114063148B publication Critical patent/CN114063148B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2415—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate
    • G06F18/24155—Bayesian classification
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00—Machine learning

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Remote Sensing (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Geology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Evolutionary Biology (AREA)
  • Software Systems (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geophysics (AREA)
  • Acoustics & Sound (AREA)
  • Probability & Statistics with Applications (AREA)
  • Mathematical Physics (AREA)
  • Computing Systems (AREA)
  • Medical Informatics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提出了一种基于贝叶斯判别的折射波初至优选方法及系统，该方法包括：获取初至数据；对初至数据进行线性动校正处理，得到初至样本，并将初至样本分成训练样本和判别样本；根据训练样本生成贝叶斯分类器；利用贝叶斯分类器，判断判别样本中的初至信息是否为有效折射波初至信息，根据判断结果得到有效折射波初至信息；根据有效折射波初至信息进行折射波静校正计算，得到炮点和检波点的折射波静校正量；本方案可以有效地实现对大数据量、低信噪比地震资料的有效折射波初至信息的快速筛选以及静校正计算，能够大大提高初至判识的效率和准确度，提高折射波静校正的精度，具有良好的应用前景。

Description

一种基于贝叶斯判别的折射波初至优选方法及系统
技术领域
本发明涉及石油勘探技术领域，尤指一种基于贝叶斯判别的折射波初至优选方法及系统。
背景技术
目前，折射波静校正是解决复杂地区中长波长静校正问题的一种有效方法。在折射波静校正方法中，一个关键的技术环节是从初至信号中优选出与有效折射波对应的初至信号，进行计算得到折射波静校正量。目前大多采用人机交互方式进行初至优级选，剔除不合理的初至信息，因而存在剔除不合理、剔除时间较长并且在操作过程中容易出错等实际问题。
针对上述问题，目前尚未提出有效的解决方案。
发明内容
为解决上述问题，本发明提出了一种基于贝叶斯判别的折射波初至优选方法及系统。该方法及系统将机器学习中的贝叶斯判别方法引入折射波静校正的初至优选中，能够解决海量地震炮集数据中折射波初至筛选中人工交互方式低、操作过程易出错等实际问题，有效提高折射波初至筛选的效率及准确性。
在本发明实施例的第一方面，提出了一种基于贝叶斯判别的折射波初至优选方法，该方法包括：
获取初至数据；
对所述初至数据进行线性动校正处理，得到初至样本，并将所述初至样本分成训练样本和判别样本；
根据所述训练样本生成贝叶斯分类器；
利用所述贝叶斯分类器，判断所述判别样本中的初至信息是否为有效折射波初至信息，根据判断结果得到有效折射波初至信息；
根据所述有效折射波初至信息进行折射波静校正计算，得到炮点和检波点的折射波静校正量。
在本发明实施例的第二方面，提出了一种基于贝叶斯判别的折射波初至优选系统，该系统包括：
数据获取模块，用于获取初至数据；
线性动校正处理模块，用于对所述初至数据进行线性动校正处理，得到初至样本，并将所述初至样本分成训练样本和判别样本；
训练模块，用于根据所述训练样本生成贝叶斯分类器；
判断模块，用于利用所述贝叶斯分类器，判断所述判别样本中的初至信息是否为有效折射波初至信息，根据判断结果得到有效折射波初至信息；
静校正计算模块，用于根据所述有效折射波初至信息进行折射波静校正计算，得到炮点和检波点的折射波静校正量。
在本发明实施例的第三方面，提出了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现基于贝叶斯判别的折射波初至优选方法。
在本发明实施例的第四方面，提出了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现基于贝叶斯判别的折射波初至优选方法。
本发明提出的基于贝叶斯判别的折射波初至优选方法及系统，适用于存在稳定折射层的复杂地表地区，可以有效地实现对大数据量、低信噪比地震资料的有效折射波初至信息的快速筛选以及静校正计算，与现有技术相比，能够大大提高初至判识的效率和准确度，提高折射波静校正的精度，具有良好的应用前景。
附图说明
为了更清楚地说明本申请实施例技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。
图1是本发明一实施例的基于贝叶斯判别的折射波初至优选方法的流程示意图。
图2是本发明一具体实施例的生成贝叶斯分类器的详细流程示意图。
图3是本发明一具体实施例的判断有效折射波初至的详细流程示意图。
图4是本发明一实施例的基于贝叶斯判别的折射波初至优选系统的架构示意图。
图5是本发明一具体实施例的训练模块的架构示意图。
图6是本发明一具体实施例的判断模块的架构示意图。
图7是本发明一实施例的计算机设备结构示意图。
具体实施方式
下面将参考若干示例性实施方式来描述本发明的原理和精神。应当理解，给出这些实施方式仅仅是为了使本领域技术人员能够更好地理解进而实现本发明，而并非以任何方式限制本发明的范围。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。
本领域技术人员知道，本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此，本公开可以具体实现为以下形式，即：完全的硬件、完全的软件(包括固件、驻留软件、微代码等)，或者硬件和软件结合的形式。
根据本发明的实施方式，提出了一种基于贝叶斯判别的折射波初至优选方法及系统，该方法及系统将机器学习中的贝叶斯判别方法引入折射波静校正的初至优选中，解决海量地震炮集数据中折射波初至筛选中人工交互方式低、操作过程易出错等实际问题，有效提高折射波初至筛选的效率及准确性。
下面参考本发明的若干代表性实施方式，详细阐释本发明的原理和精神。
图1是本发明一实施例的基于贝叶斯判别的折射波初至优选方法的流程示意图。如图1所示，该方法包括：
步骤S101，获取初至数据；
步骤S102，对所述初至数据进行线性动校正处理，得到初至样本，并将所述初至样本分成训练样本和判别样本；
步骤S103，根据所述训练样本生成贝叶斯分类器；
步骤S104，利用所述贝叶斯分类器，判断所述判别样本中的初至信息是否为有效折射波初至信息，根据判断结果得到有效折射波初至信息；
步骤S105，根据所述有效折射波初至信息进行折射波静校正计算，得到炮点和检波点的折射波静校正量。
为了对上述基于贝叶斯判别的折射波初至优选方法进行更为清楚的解释，下面结合每一步骤进行详细说明。
步骤S101：
读取工区内微测井提供的高速层速度数据、地震炮集数据对应的初至数据。
步骤S102：
对所述初至数据进行线性动校正处理，得到初至样本，并将所述初至样本分成训练样本和判别样本。
具体的，对初至数据进行线性动校正处理，线性动校正速度为高速层的地震波传播速度，例如，可用微测井资料提供的高速层速度资料。
在线性动校正后的初至时间服从高斯分布(式1)或混合高斯分布(式2)：
其中，t为线性动校后的初至走时，如式3：
t0为原始初至走时，x为对应的偏移距，v是选用的线性动校正速度，通常选用微测井资料提供的高速层速度作为线性动校正速度，使得线性动校正后的初至时间大致服从零均值高斯分布；
在式1中，μ和σ分别为高斯模型的均值和方差。
在式2中，M为混合高斯模型的个数，λc是每个高斯模型所占的比例，μc和σc分别为每个高斯模型的均值和方差。
进一步的，可以判断线性动校后的折射波初至走时服从高斯分布或混合高斯分布；
如服从高斯分布，计算相应的均值和方差：
ti为高斯分布模型的初至走时，N是样本的个数。
如果服从混合高斯模型，则先判别混合高斯模型的个数，通常选取M＝2或3即可满足实际需求。利用AIC准则判断混合高斯模型的个数，即定阶。AIC准则的计算公式如下：
AIC＝2*NlogL+2*M； (6)
其中，N是样本的个数，L是样本的负对数似然函数，M是混合高斯模型的个数。
分别计算混合高斯模型个数M＝2和3对应的AIC值，如果AIC(M＝2)<AIC(M＝3)，则认为线性动校后的初至服从M＝2的混合高斯模型；如果AIC(M＝2)>AIC(M＝3)，认为初至服从M＝3的混合高斯模型。
步骤S103：
根据所述训练样本生成贝叶斯分类器。
参考图2，为本发明一具体实施例的生成贝叶斯分类器的详细流程示意图。如图2所示，步骤S103的详细过程为：
步骤S1031，在训练样本中拾取有效折射波初至(第一类样本)以及非有效折射波初至(第二类样本)，其中，有效折射波初至为信号，非有效折射波初至为噪声。
步骤S1032，设置信号和噪声的先验概率P(ω1)与P(ω2)及信号和噪声的条件概率密度P(t|ω1)和P(t|ω2)，在设置的条件下利用概率函数估计方法预估信号样本和噪声样本的概率密度函数。概率函数估计方法可以采用如Parzen窗法等，即可估计此条件下的概率密度函数(PDF，probability density function)。
步骤S1033，根据所述概率密度函数生成贝叶斯分类器。
步骤S104：
利用所述贝叶斯分类器，判断所述判别样本中的初至信息是否为有效折射波初至信息，根据判断结果得到有效折射波初至信息。
参考图3，为本发明一具体实施例的判断有效折射波初至的详细流程示意图。如图3所示，步骤S104的详细过程为：
步骤S1041，根据所述判别样本中的初至信息，在先验概率P(ω1)与P(ω2)及条件概率密度P(t|ω1)和P(t|ω2)的条件下，利用贝叶斯公式把先验概率转换为后验概率；其中，样本总体概率密度P(t)分布为：
P(t)＝P(t|ω1)P(ω1)+P(t|ω2)P(ω2)； (7)
后验概率为：
其中，P(ω1|t)为信号的后验概率，P(ω2|t)为噪声的后验概率；
步骤S1042，根据式7、8、9得到后验概率，根据后验概率判断所述判别样本中的初至信息是否为有效折射波初至信息；其中，
若P(ω1|t)＞P(ω2|t)，判断当前初至信息为有效折射波初至信息；
若P(ω1|t)≤P(ω2|t)，判断当前初至信息为非有效折射波初至信息(噪声)。
步骤S1043，根据判断结果得到有效折射波初至信息。
步骤S105：
根据所述有效折射波初至信息进行折射波静校正计算，得到炮点和检波点的折射波静校正量。
本发明提出的基于贝叶斯判别的折射波初至优选方法将机器学习中的贝叶斯判别方法引入折射波静校正的初至优选中，能够解决海量地震炮集数据中折射波初至筛选中人工交互方式低、操作过程易出错等实际问题，有效提高折射波初至筛选的效率及准确性。
需要说明的是，尽管在上述实施例及附图中以特定顺序描述了本发明方法的操作，但是，这并非要求或者暗示必须按照该特定顺序来执行这些操作，或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地，可以省略某些步骤，将多个步骤合并为一个步骤执行，和/或将一个步骤分解为多个步骤执行。
在介绍了本发明示例性实施方式的方法之后，接下来，参考图4至图6对本发明示例性实施方式的基于贝叶斯判别的折射波初至优选系统进行介绍。
基于贝叶斯判别的折射波初至优选系统的实施可以参见上述方法的实施，重复之处不再赘述。以下所使用的术语“模块”或者“单元”，可以是实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
基于同一发明构思，本发明还提出了一种基于贝叶斯判别的折射波初至优选系统，如图4所示，该系统包括：
数据获取模块410，用于获取初至数据；
线性动校正处理模块420，用于对所述初至数据进行线性动校正处理，得到初至样本，并将所述初至样本分成训练样本和判别样本；
训练模块430，用于根据所述训练样本生成贝叶斯分类器；
判断模块440，用于利用所述贝叶斯分类器，判断所述判别样本中的初至信息是否为有效折射波初至信息，根据判断结果得到有效折射波初至信息；
静校正计算模块450，用于根据所述有效折射波初至信息进行折射波静校正计算，得到炮点和检波点的折射波静校正量。
在一实施例中，所述数据获取模块410具体用于：
读取工区内微测井提供的高速层速度数据、地震炮集数据对应的初至数据。
在一实施例中，参考图5，为本发明一具体实施例的训练模块的架构示意图。如图5所示，所述训练模块430包括：
拾取单元431，用于在训练样本中拾取有效折射波初至以及非有效折射波初至，其中，有效折射波初至为信号，非有效折射波初至为噪声；
概率密度函数预估单元432，用于设置信号和噪声的先验概率P(ω1)与P(ω2)及信号和噪声的条件概率密度P(t|ω1)和P(t|ω2)，在设置的条件下利用概率函数估计方法预估信号样本和噪声样本的概率密度函数；
贝叶斯分类器生成单元433，用于根据所述概率密度函数生成贝叶斯分类器。
在一实施例中，参考图6，为本发明一具体实施例的判断模块的架构示意图。如图6所示，所述判断模块440包括：
转换单元441，用于根据所述判别样本中的初至信息，在先验概率P(ω1)与P(ω2)及条件概率密度P(t|ω1)和P(t|ω2)的条件下，利用贝叶斯公式把先验概率转换为后验概率；其中，样本总体概率密度P(t)分布为：
P(t)＝P(t|ω1)P(ω1)+P(t|ω2)P(ω2)；
后验概率为：
其中，P(ω1|t)为信号的后验概率，P(ω2|t)为噪声的后验概率；
判断单元442，用于判断所述判别样本中的初至信息是否为有效折射波初至信息；其中，
若P(ω1|t)＞P(ω2|t)，判断当前初至信息为有效折射波初至信息；
若P(ω1|t)≤P(ω2|t)，判断当前初至信息为非有效折射波初至信息；
有效信息获取单元443，用于根据判断结果得到有效折射波初至信息。
应当注意，尽管在上文详细描述中提及了基于贝叶斯判别的折射波初至优选系统的若干模块，但是这种划分仅仅是示例性的并非强制性的。实际上，根据本发明的实施方式，上文描述的两个或更多模块的特征和功能可以在一个模块中具体化。反之，上文描述的一个模块的特征和功能可以进一步划分为由多个模块来具体化。
基于前述发明构思，如图7所示，本发明还提出了一种计算机设备700，包括存储器710、处理器720及存储在存储器710上并可在处理器720上运行的计算机程序730，所述处理器720执行所述计算机程序730时实现前述基于贝叶斯判别的折射波初至优选方法。
基于前述发明构思，本发明提出了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现前述基于贝叶斯判别的折射波初至优选方法。
本发明提出的基于贝叶斯判别的折射波初至优选方法及系统，适用于存在稳定折射层的复杂地表地区，可以有效地实现对大数据量、低信噪比地震资料的有效折射波初至信息的快速筛选以及静校正计算，与现有技术相比，能够大大提高初至判识的效率和准确度，提高折射波静校正的精度，具有良好的应用前景。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。

Claims (8)

1.一种基于贝叶斯判别的折射波初至优选方法，其特征在于，该方法包括：
S101，获取初至数据；
S102，对所述初至数据进行线性动校正处理，得到初至样本，并将所述初至样本分成训练样本和判别样本；
S103，根据所述训练样本生成贝叶斯分类器；
S104，利用所述贝叶斯分类器，判断所述判别样本中的初至信息是否为有效折射波初至信息，根据判断结果得到有效折射波初至信息；
S105，根据所述有效折射波初至信息进行折射波静校正计算，得到炮点和检波点的折射波静校正量；
其中，S101，获取初至数据，包括：
读取工区内微测井提供的高速层速度数据、地震炮集数据对应的初至数据；
其中，S102，对所述初至数据进行线性动校正处理，得到初至样本，并将所述初至样本分成训练样本和判别样本，包括：
对初至数据进行线性动校正处理，线性动校正速度为高速层的地震波传播速度；
在线性动校正后的初至时间服从式1的高斯分布或式2的混合高斯分布：
式1中，μ和σ分别为高斯模型的均值和方差；
式2中，M为混合高斯模型的个数，λc是每个高斯模型所占的比例，μc和σc分别为每个高斯模型的均值和方差；
t为线性动校后的初至走时，如式3：
式3中，t0为原始初至走时，x为对应的偏移距，v是选用的线性动校正速度，选用微测井提供的高速层速度数据作为线性动校正速度，使得线性动校正后的初至时间大致服从零均值高斯分布；
进一步的，判断线性动校后的折射波初至走时服从高斯分布或混合高斯分布；
如服从高斯分布，计算相应的均值和方差：
式4及式5中，ti为高斯分布模型的初至走时，N是样本的个数；
如果服从混合高斯模型，则先判别混合高斯模型的个数，选取M＝2或3；利用AIC准则判断混合高斯模型的个数；AIC准则的计算公式如下：
AIC＝2*NlogL+2*M； (6)
式6中，N是样本的个数，L是样本的负对数似然函数，M是混合高斯模型的个数；
分别计算混合高斯模型个数M＝2和3对应的AIC值，如果M＝2的AIC<M＝3的AIC，则认为线性动校后的初至服从M＝2的混合高斯模型；如果M＝2的AIC>M＝3的AIC，则认为线性动校后的初至服从M＝3的混合高斯模型。
2.根据权利要求1所述的基于贝叶斯判别的折射波初至优选方法，其特征在于，根据所述训练样本生成贝叶斯分类器，包括：
在训练样本中拾取有效折射波初至以及非有效折射波初至，其中，有效折射波初至为信号，非有效折射波初至为噪声；
设置信号和噪声的先验概率P(ω1)与P(ω2)及信号和噪声的条件概率密度P(t|ω1)和P(t|ω2)，在设置的条件下利用概率函数估计方法预估信号样本和噪声样本的概率密度函数；
根据所述概率密度函数生成贝叶斯分类器。
3.根据权利要求2所述的基于贝叶斯判别的折射波初至优选方法，其特征在于，利用所述贝叶斯分类器，判断所述判别样本中的初至信息是否为有效折射波初至信息，根据判断结果得到有效折射波初至信息，包括：
根据所述判别样本中的初至信息，在先验概率P(ω1)与P(ω2)及条件概率密度P(t|ω1)和P(t|ω2)的条件下，利用贝叶斯公式把先验概率转换为后验概率；其中，样本总体概率密度P(t)分布为：
P(t)＝P(t|ω1)P(ω1)+P(t|ω2)P(ω2)；
后验概率为：
其中，P(ω1|t)为信号的后验概率，P(ω2|t)为噪声的后验概率；
判断所述判别样本中的初至信息是否为有效折射波初至信息；其中，
若P(ω1|t)＞P(ω2|t)，判断当前初至信息为有效折射波初至信息；
若P(ω1|t)≤P(ω2|t)，判断当前初至信息为非有效折射波初至信息；
根据判断结果得到有效折射波初至信息。
4.一种基于贝叶斯判别的折射波初至优选系统，其特征在于，该系统包括：
数据获取模块，用于获取初至数据；
线性动校正处理模块，用于对所述初至数据进行线性动校正处理，得到初至样本，并将所述初至样本分成训练样本和判别样本；
训练模块，用于根据所述训练样本生成贝叶斯分类器；
判断模块，用于利用所述贝叶斯分类器，判断所述判别样本中的初至信息是否为有效折射波初至信息，根据判断结果得到有效折射波初至信息；
静校正计算模块，用于根据所述有效折射波初至信息进行折射波静校正计算，得到炮点和检波点的折射波静校正量；
其中，所述数据获取模块具体用于：
读取工区内微测井提供的高速层速度数据、地震炮集数据对应的初至数据；
其中，所述线性动校正处理模块具体用于：
对初至数据进行线性动校正处理，线性动校正速度为高速层的地震波传播速度；
在线性动校正后的初至时间服从式1的高斯分布或式2的混合高斯分布：
式1中，μ和σ分别为高斯模型的均值和方差；
式2中，M为混合高斯模型的个数，λc是每个高斯模型所占的比例，μc和σc分别为每个高斯模型的均值和方差；
t为线性动校后的初至走时，如式3：
式3中，t0为原始初至走时，x为对应的偏移距，v是选用的线性动校正速度，选用微测井提供的高速层速度数据作为线性动校正速度，使得线性动校正后的初至时间大致服从零均值高斯分布；
进一步的，判断线性动校后的折射波初至走时服从高斯分布或混合高斯分布；
如服从高斯分布，计算相应的均值和方差：
式4及式5中，ti为高斯分布模型的初至走时，N是样本的个数；
如果服从混合高斯模型，则先判别混合高斯模型的个数，选取M＝2或3；利用AIC准则判断混合高斯模型的个数；AIC准则的计算公式如下：
AIC＝2*NlogL+2*M； (6)
式6中，N是样本的个数，L是样本的负对数似然函数，M是混合高斯模型的个数；
分别计算混合高斯模型个数M＝2和3对应的AIC值，如果M＝2的AIC<M＝3的AIC，则认为线性动校后的初至服从M＝2的混合高斯模型；如果M＝2的AIC>M＝3的AIC，则认为线性动校后的初至服从M＝3的混合高斯模型。
5.根据权利要求4所述的基于贝叶斯判别的折射波初至优选系统，其特征在于，所述训练模块包括：
拾取单元，用于在训练样本中拾取有效折射波初至以及非有效折射波初至，其中，有效折射波初至为信号，非有效折射波初至为噪声；
概率密度函数预估单元，用于设置信号和噪声的先验概率P(ω1)与P(ω2)及信号和噪声的条件概率密度P(t|ω1)和P(t|ω2)，在设置的条件下利用概率函数估计方法预估信号样本和噪声样本的概率密度函数；
贝叶斯分类器生成单元，用于根据所述概率密度函数生成贝叶斯分类器。
6.根据权利要求5所述的基于贝叶斯判别的折射波初至优选系统，其特征在于，所述判断模块包括：
转换单元，用于根据所述判别样本中的初至信息，在先验概率P(ω1)与P(ω2)及条件概率密度P(t|ω1)和P(t|ω2)的条件下，利用贝叶斯公式把先验概率转换为后验概率；其中，样本总体概率密度P(t)分布为：
P(t)＝P(t|ω1)P(ω1)+P(t|ω2)P(ω2)；
后验概率为：
其中，P(ω1|t)为信号的后验概率，P(ω2|t)为噪声的后验概率；
判断单元，用于判断所述判别样本中的初至信息是否为有效折射波初至信息；其中，
若P(ω1|t)＞P(ω2|t)，判断当前初至信息为有效折射波初至信息；
若P(ω1|t)≤P(ω2|t)，判断当前初至信息为非有效折射波初至信息；
有效信息获取单元，用于根据判断结果得到有效折射波初至信息。
7.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至3任一所述方法。
8.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现权利要求1至3任一所述方法。
CN202010756674.XA 2020-07-31 2020-07-31 一种基于贝叶斯判别的折射波初至优选方法及系统 Active CN114063148B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010756674.XA CN114063148B (zh) 2020-07-31 2020-07-31 一种基于贝叶斯判别的折射波初至优选方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010756674.XA CN114063148B (zh) 2020-07-31 2020-07-31 一种基于贝叶斯判别的折射波初至优选方法及系统

Publications (2)

Publication Number Publication Date
CN114063148A CN114063148A (zh) 2022-02-18
CN114063148B true CN114063148B (zh) 2023-09-26

Family

ID=80227636

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010756674.XA Active CN114063148B (zh) 2020-07-31 2020-07-31 一种基于贝叶斯判别的折射波初至优选方法及系统

Country Status (1)

Country Link
CN (1) CN114063148B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109344881A (zh) * 2018-09-11 2019-02-15 中国科学技术大学 一种基于时空连续性的扩展分类器
CN111239805A (zh) * 2020-02-13 2020-06-05 中国石油大学(北京) 基于反射率法的块约束时移地震差异反演方法及系统
CN111381275A (zh) * 2019-11-20 2020-07-07 中国海洋石油集团有限公司 一种地震数据的初至拾取方法和装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9971051B2 (en) * 2013-01-14 2018-05-15 Cgg Services Sas High-fidelity adaptive curvelet domain primary-multiple separation processing of seismic data

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109344881A (zh) * 2018-09-11 2019-02-15 中国科学技术大学 一种基于时空连续性的扩展分类器
CN111381275A (zh) * 2019-11-20 2020-07-07 中国海洋石油集团有限公司 一种地震数据的初至拾取方法和装置
CN111239805A (zh) * 2020-02-13 2020-06-05 中国石油大学(北京) 基于反射率法的块约束时移地震差异反演方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
三维折射波静校正;王彦春 等;《石油地球物理勘探》;第35卷(第1期);第13-19页 *
地面微地震资料弱信号提取方法研究;宋维琪 等;《石油物探》;第52卷(第2期);第131-136页 *

Also Published As

Publication number Publication date
CN114063148A (zh) 2022-02-18

Similar Documents

Publication Publication Date Title
Bokulich Using models to correct data: Paleodiversity and the fossil record
CN108805269B (zh) 一种基于lstm循环神经网络拾取震相到时的方法
US20210026031A1 (en) Dolomite reservoir prediction method and system based on well and seismic combination, and storage medium
Mahasin et al. Comparison of CSPDarkNet53, CSPResNeXt-50, and EfficientNet-B0 backbones on YOLO v4 as object detector
CN112990330B (zh) 用户用能异常数据检测方法及设备
CN113687433B (zh) 一种基于Bi-LSTM的大地电磁信号去噪方法及系统
CN108334987B (zh) 一种基于小波分解-神经网络的海浪波高预测方法
CN115408925A (zh) 用于隧道施工的岩体参数预测方法和装置
US20220092359A1 (en) Image data classification method, device and system
CN116626753B (zh) 一种基于多模态神经网络的微地震事件识别方法及系统
Niksejel et al. OBSTransformer: a deep-learning seismic phase picker for OBS data using automated labelling and transfer learning
CN115032682A (zh) 一种基于图论的多站台地震震源参数估计方法
CN109165603B (zh) 一种舰船检测方法及装置
CN114063148B (zh) 一种基于贝叶斯判别的折射波初至优选方法及系统
CN112882899B (zh) 一种日志异常检测方法及装置
CN114548161A (zh) 一种重介分选精煤灰分预测方法、装置、电子设备及介质
CN112735466B (zh) 一种音频检测方法及装置
CN118330471A (zh) 一种基于数据增强模型的电池eis数据增强方法及相关设备
CN116541668B (zh) 一种游泳划水次数确定方法、装置、设备及存储介质
Ren et al. Automated dispersion curve picking using multi-attribute convolutional-neural-network based machine learning
CN114494999A (zh) 一种双分支联合型目标密集预测方法及系统
Dubus et al. From citizen science to AI models: Advancing cetacean vocalization automatic detection through multi-annotator campaigns
CN113361579A (zh) 一种水下目标探测辨识方法、系统、设备及可读存储介质
Shi A Survey of Changepoint Techniques for Time Series Data
Armstrong et al. A deep‐learning phase picker with calibrated Bayesian‐derived uncertainties for earthquakes in the Yellowstone volcanic region

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