CN111897001A - 波场确定方法与装置 - Google Patents

波场确定方法与装置 Download PDF

Info

Publication number
CN111897001A
CN111897001A CN202010568542.4A CN202010568542A CN111897001A CN 111897001 A CN111897001 A CN 111897001A CN 202010568542 A CN202010568542 A CN 202010568542A CN 111897001 A CN111897001 A CN 111897001A
Authority
CN
China
Prior art keywords
target
slowness
seismic wave
seismic
wave data
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
Application number
CN202010568542.4A
Other languages
English (en)
Other versions
CN111897001B (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.)
Chengdu Wushi Technology Co.,Ltd.
Original Assignee
Chengdu Univeristy of Technology
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 Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN202010568542.4A priority Critical patent/CN111897001B/zh
Publication of CN111897001A publication Critical patent/CN111897001A/zh
Application granted granted Critical
Publication of CN111897001B publication Critical patent/CN111897001B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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/364—Seismic filtering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Geophysics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Acoustics & Sound (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明实施例公开了一种波场确定的方法与装置，所述方法包括：基于预先构建的慢度谱确定方法，确定目标地震波数据在τ‑p域上的目标慢度谱；基于目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数；基于目标地震波的慢度拟合直线参数，确定与目标地震波对应的目标时移量，并根据目标时移量，对目标地震波数据进行排齐处理，以使排齐处理后的目标地震波数据中的同相轴平行排齐，目标地震波为不同类型的地震波中的任意一种类型的地震波；对排齐处理后的目标地震波数据进行过滤处理，并基于目标时移量，对过滤处理后的目标地震波数据进行反排齐处理，得到与目标地震波对应的目标波场。通过本方法，可以确定地震波对应的波场。

Description

波场确定方法与装置
技术领域
本发明涉及地震勘探领域，尤其涉及一种波场确定方法与装置。
背景技术
地震勘探方法是人类获取地下空间信息的重要手段，其中，对垂直地震剖面(Vertical Seismic Profiling，VSP)的地震波数据中包含的不同类型的地震波波场进行分离，是对地震波数据进行处理的关键步骤。
目前，可以通过f-k滤波方法，对地震波数据中的不同类型的地震波波场进行分离，但是，由于不同f-k象限之间存在能量谱的重叠，所以会导致分离出的地震波波场的准确性差，因此，需要提供一种准确性更高的波场确定方案。
发明内容
本发明实施例的目的是提供一种波场确定方法与装置，以解决现有技术中在从地震波数据中确定地震波波场时，存在的波场确定准确性差的问题。
为解决上述技术问题，本发明实施例是这样实现的：
第一方面，本发明实施例提供的一种波场确定方法，所述方法包括：
基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱，所述慢度谱确定方法是根据检波器的走时、所述检波器的深度、所述检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，确定慢度谱的方法；
基于所述目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数；
基于目标地震波的所述慢度拟合直线参数，确定与所述目标地震波对应的目标时移量，并根据所述目标时移量，对所述目标地震波数据进行排齐处理，以使排齐处理后的所述目标地震波数据中的同相轴平行排齐，所述目标地震波为所述不同类型的地震波中的任意一种类型的地震波；
对所述排齐处理后的目标地震波数据进行过滤处理，并基于所述目标时移量，对过滤处理后的所述目标地震波数据进行反排齐处理，得到与所述目标地震波对应的目标波场。
可选地，在所述基于所述目标时移量，对过滤处理后的所述目标地震波数据进行反排齐处理，得到与所述目标地震波对应的目标波场之后，还包括：
基于所述目标地震波数据和所述目标波场，确定第一地震波数据，所述第一地震波为所述不同类型的地震波中，与所述目标地震波的类型不同的任意一种类型的地震波；
基于第一地震波的所述慢度拟合直线参数，确定与所述第一地震波对应的第一时移量，并根据所述第一时移量，对所述第一地震波数据进行排齐处理，以使排齐处理后的所述第一地震波数据中的同相轴平行排齐；
对所述排齐处理后的第一地震波数据进行过滤处理，并基于所述第一时移量，对过滤处理后的所述第一地震波数据进行反排齐处理，得到与所述第一地震波对应的第一波场。
可选地，所述基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱，包括：
将每个所述检波器的走时、所述检波器的深度、所述检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，代入公式
Figure BDA0002548677570000021
得到所述目标地震波数据在τ-p域上的目标慢度谱，其中，v(p,τ)为所述目标地震波数据在τ-p域上的目标慢度谱，p为慢度值，τ为所述截距时间，M为所述预设时窗长度，hi为第i个所述检波器的深度，t为第i个所述检波器的走时，
Figure BDA0002548677570000022
为第i个所述检波器对应的叠加路径内振幅序列的方差，Nx为所述检波器的数量，m＝-M/2，u(hi,t)为第i个所述检波器在地震域基于所述检波器的走时和所述检波器的深度确定的坐标。
可选地，在所述基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱之前，还包括：
基于所述目标地震波数据，获取每个所述检波器的截距时间、慢度值以及深度；
将每个所述检波器的截距时间、慢度值以及深度代入公式
Figure BDA0002548677570000023
得到每个所述检波器的走时，其中，t为所述检波器的走时，τ为所述检波器的截距时间，p为所述检波器的慢度值，hi为所述检波器的深度。
可选地，所述基于所述目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数，包括：
基于所述不同类型的地震波对应的预设慢度值范围，确定所述目标慢度谱中与每个类型的所述地震波对应的慢度点；
基于所述每个类型的地震波对应的慢度点的慢度值，确定所述每个类型的地震波对应的慢度拟合直线参数。
可选地，所述基于目标地震波的所述慢度拟合直线参数，确定与所述目标地震波对应的目标时移量，包括：
将所述目标地震波的所述慢度拟合直线参数，代入公式
Figure BDA0002548677570000031
得到与所述目标地震波对应的目标时移量，其中，i为第i个所述检波器，Δti为第i个所述检波器的目标时移量，ti为第i个所述检波器的走时，Δh＝hi-hi-1，hi为第i个所述检波器的深度，an为所述慢度拟合直线参数中的第一慢度拟合直线参数，bn为所述慢度拟合直线参数中的第二慢度拟合直线参数。
可选地，所述对所述排齐处理后的目标地震波数据进行过滤处理，包括：
基于均值滤波算法，对所述排齐处理后的目标地震波数据进行过滤处理。
第二方面，本发明实施例提供了一种波场确定装置，所述装置包括：
慢度谱确定模块，用于基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱，所述慢度谱确定方法是根据检波器的走时、所述检波器的深度、所述检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，确定慢度谱的方法；
参数确定模块，用于基于所述目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数；
第一确定模块，用于基于目标地震波的所述慢度拟合直线参数，确定与所述目标地震波对应的目标时移量，并根据所述目标时移量，对所述目标地震波数据进行排齐处理，以使排齐处理后的所述目标地震波数据中的同相轴平行排齐，所述目标地震波为所述不同类型的地震波中的任意一种类型的地震波；
第一波场确定模块，用于对所述排齐处理后的目标地震波数据进行过滤处理，并基于所述目标时移量，对过滤处理后的所述目标地震波数据进行反排齐处理，得到与所述目标地震波对应的目标波场。
可选地，所述装置，还包括：
数据确定模块，用于基于所述目标地震波数据和所述目标波场，确定第一地震波数据，所述第一地震波为所述不同类型的地震波中，与所述目标地震波的类型不同的任意一种类型的地震波；
第二确定模块，用于基于第一地震波的所述慢度拟合直线参数，确定与所述第一地震波对应的第一时移量，并根据所述第一时移量，对所述第一地震波数据进行排齐处理，以使排齐处理后的所述第一地震波数据中的同相轴平行排齐；
第二波场确定模块，用于对所述排齐处理后的第一地震波数据进行过滤处理，并基于所述第一时移量，对过滤处理后的所述第一地震波数据进行反排齐处理，得到与所述第一地震波对应的第一波场。
可选地，所述慢度谱确定模块，用于：
将每个所述检波器的走时、所述检波器的深度、所述检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，代入公式
Figure BDA0002548677570000041
得到所述目标地震波数据在τ-p域上的目标慢度谱，其中，v(p,τ)为所述目标地震波数据在τ-p域上的目标慢度谱，p为慢度值，τ为所述截距时间，M为所述预设时窗长度，hi为第i个所述检波器的深度，t为第i个所述检波器的走时，
Figure BDA0002548677570000042
为第i个所述检波器对应的叠加路径内振幅序列的方差，Nx为所述检波器的数量，m＝-M/2，u(hi,t)为第i个所述检波器在地震域基于所述检波器的走时和所述检波器的深度确定的坐标。
可选地，所述装置，还包括：
获取模块，用于基于所述目标地震波数据，获取每个所述检波器的截距时间、慢度值以及深度；
走时确定模块，用于将每个所述检波器的截距时间、慢度值以及深度代入公式
Figure BDA0002548677570000043
得到每个所述检波器的走时，其中，t为所述检波器的走时，τ为所述检波器的截距时间，p为所述检波器的慢度值，hi为所述检波器的深度。
可选地，所述参数确定模块，用于：
基于所述不同类型的地震波对应的预设慢度值范围，确定所述目标慢度谱中与每个类型的所述地震波对应的慢度点；
基于所述每个类型的地震波对应的慢度点的慢度值，确定所述每个类型的地震波对应的慢度拟合直线参数。
可选地，所述第一确定模块，用于：
将所述目标地震波的所述慢度拟合直线参数，代入公式
Figure BDA0002548677570000051
得到与所述目标地震波对应的目标时移量，其中，i为第i个所述检波器，Δti为第i个所述检波器的目标时移量，ti为第i个所述检波器的走时，Δh＝hi-hi-1，hi为第i个所述检波器的深度，an为所述慢度拟合直线参数中的第一慢度拟合直线参数，bn为所述慢度拟合直线参数中的第二慢度拟合直线参数。
可选地，所述第一波场确定模块，用于：
基于均值滤波算法，对所述排齐处理后的目标地震波数据进行过滤处理。
第三方面，本发明实施例提供一种电子设备，包括处理器、存储器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述计算机程序被所述处理器执行时实现上述第一方面提供的波场确定方法的步骤。
第四方面，本发明实施例提供一种计算机可读存储介质，所述计算机可读存储介质上存储计算机程序，所述计算机程序被处理器执行时实现上述第一方面提供的波场确定方法的步骤。
由以上本发明实施例提供的技术方案可见，本发明实施例通过基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱，慢度谱确定方法是根据检波器的走时、检波器的深度、检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，确定慢度谱的方法，基于目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数，基于目标地震波的慢度拟合直线参数，确定与目标地震波对应的目标时移量，并根据目标时移量，对目标地震波数据进行排齐处理，以使排齐处理后的目标地震波数据中的同相轴平行排齐，目标地震波为不同类型的地震波中的任意一种类型的地震波，对排齐处理后的目标地震波数据进行过滤处理，并基于目标时移量，对过滤处理后的目标地震波数据进行反排齐处理，得到与目标地震波对应的目标波场。这样，通过τ-p域上的目标慢度谱，可以准确的获取不同类型的地震波对应的慢度拟合直线参数，在根据目标地震波的慢度拟合直线参数，对目标地震波数据进行排齐处理、过滤处理以及反排齐处理，以提高目标波场从目标地震波数据中分离出来的准确性，同时也可以抑制目标地震波数据中的噪声的影响，即可以提高目标地震波的目标波场的确定准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明一种波场确定方法的流程示意图；
图2为本发明另一种波场确定方法的流程示意图；
图3为本发明一种波场确定装置的结构示意图；
图4为本发明一种电子设备的结构示意图。
具体实施方式
本发明实施例提供一种波场确定方法与装置。
为了使本技术领域的人员更好地理解本发明中的技术方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
实施例一
如图1所示，本说明书实施例提供一种波场确定方法，该方法的执行主体可以为服务器，该服务器可以是独立的服务器，也可以是由多个服务器组成的服务器集群。该方法具体可以包括以下步骤：
在S102中，基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱。
其中，慢度谱确定方法可以是根据检波器的走时、检波器的深度、检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，确定慢度谱的方法，目标地震波数据可以为通过相关的勘探器件进行实地勘探时记录的勘探资料中获取的已知数据。
在实施中，地震勘探方法是人类获取地下空间信息的重要手段，其中，对垂直地震剖面(VerticalSeismicProfiling，VSP)的地震波数据中包含的不同类型的地震波波场进行分离，是对地震波数据进行处理的关键步骤。
目前，可以通过f-k滤波方法，对地震波数据中的不同类型的地震波波场进行分离，但是，由于不同f-k象限之间存在能量谱的重叠，所以会导致分离出的地震波波场的准确性差，因此，需要提供一种准确性更高的波场确定方案。为此，本发明实施例提供一种能够解决上述问题的技术方案，具体可以参见下述内容：
服务器可以从预先存储的实际勘探资料中查找相关信息，通过查找的相关信息可以得到目标地震波数据。
服务器可以根据预先构建的慢度谱确定方法，将地震域中的目标地震波数据，转换为τ-p域上的目标慢度谱，通过慢度谱确定方法中的预设时长长度，可以将目标地震波数据转换到叠加范围更大的τ-p域的时窗内，转换得到的目标慢度谱可以显示不同截距时间和慢度值对应的叠加路径内振幅相似性以及地震波能量的强弱，同时，目标慢度谱上的能量极值点也是计算地震波波场的合理慢度的重要依据。
在S104中，基于目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数。
其中，不同类型的地震波可以包括上行P波、上行S波、下行P波以及下行S波等，不同类型的地震波在目标慢度谱中对应的慢度值的范围可以不同。
在实施中，不同类型的地震波，在目标慢度谱上的变化情况时不同的，假设目标慢度谱上的地震波的慢度值是随截距时间的变化满足双曲线特征的，服务器可以在τ-p域上拾取不同类型的地震波的对应的慢度值，并根据不同类型的地震波的慢度值，确定不同类型的地震波的慢度拟合直线参数。
例如，服务器可以在目标慢度谱，拾取针对不同类型的地震波对应的慢度点(pi,τi)，其中，pi为第i个慢度点对应的慢度值，τi为第i个慢度点对应的截距时间，且这些慢度点满足：
Figure BDA0002548677570000081
其中，n为第n个类型的地震波，pnmin为第n个类型的地震波对应的慢度值中的最小值，pnmax为第n个类型的地震波对应的慢度值中的最大值，pn(τ)为第n个类型的地震波的第τ个截距时间对应的慢度值。
在S106中，基于目标地震波的慢度拟合直线参数，确定与目标地震波对应的目标时移量，并根据目标时移量，对目标地震波数据进行排齐处理，以使排齐处理后的目标地震波数据中的同相轴平行排齐。
其中，目标地震波可以为不同类型的地震波中的任意一种类型的地震波。
在实施中，由于地震波的速度和地震波的到达方向会影响检波器的时移，而上行波和下行波中包含来自不同界面的反射波，即使是同一类型地震波，其传播方向也并不完全一致。这就导致静态时移并不能准确地将同相轴排齐，因此，需要通过确定的目标时移量，对目标地震波进行排齐处理。
在S108中，对排齐处理后的目标地震波数据进行过滤处理，并基于目标时移量，对过滤处理后的目标地震波数据进行反排齐处理，得到与目标地震波对应的目标波场。
在实施中，对排齐处理后的目标地震波数据进行过滤处理的方法可以有多个，可以根据实际应用场景的不同而有所不同，本发明实施例对此不作具体限定。
可以基于上述得到目标地震波对应的目标波场的方法，确定其他类型的地震波对应的波场，以实现对目标地震波数据的波场分离。
本发明实施例提供一种波场确定方法，本发明实施例通过基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱，慢度谱确定方法是根据检波器的走时、检波器的深度、检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，确定慢度谱的方法，基于目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数，基于目标地震波的慢度拟合直线参数，确定与目标地震波对应的目标时移量，并根据目标时移量，对目标地震波数据进行排齐处理，以使排齐处理后的目标地震波数据中的同相轴平行排齐，目标地震波为不同类型的地震波中的任意一种类型的地震波，对排齐处理后的目标地震波数据进行过滤处理，并基于目标时移量，对过滤处理后的目标地震波数据进行反排齐处理，得到与目标地震波对应的目标波场。这样，通过τ-p域上的目标慢度谱，可以准确的获取不同类型的地震波对应的慢度拟合直线参数，在根据目标地震波的慢度拟合直线参数，对目标地震波数据进行排齐处理、过滤处理以及反排齐处理，以提高目标波场从目标地震波数据中分离出来的准确性，同时也可以抑制目标地震波数据中的噪声的影响，即可以提高目标地震波的目标波场的确定准确性。
实施例二
如图2所示，本发明实施例提供一种波场确定方法，该方法的执行主体可以为服务器，该服务器可以是独立的服务器，也可以是由多个服务器组成的服务器集群。该方法具体可以包括以下步骤：
在S202中，基于目标地震波数据，获取每个检波器的截距时间、慢度值以及深度。
其中，目标地震波数据对应的检波器的数量可以大于预设检波器数量阈值(如可以为100个)。
在S204中，将每个检波器的截距时间、慢度值以及深度代入公式
Figure BDA0002548677570000091
得到每个检波器的走时。
其中，t为检波器的走时，τ为检波器的截距时间，p为检波器的慢度值，hi为第i个检波器的深度。
在实施中，当检波器的数量较多(即大于预设检波器数量阈值)时，目标地震波数据中包含的地震波可以为曲面波，并满足双曲线关系，所以，可以通过上述公式，确定每个检波器的走时。
在S206中，将每个检波器的走时、检波器的深度、检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，代入公式
Figure BDA0002548677570000092
得到目标地震波数据在τ-p域上的目标慢度谱。
其中，v(p,τ)为目标地震波数据在τ-p域上的目标慢度谱，p为慢度值，τ为截距时间，M为预设时窗长度，hi为第i个检波器的深度，t为第i个检波器的走时，
Figure BDA0002548677570000093
为第i个检波器对应的叠加路径内振幅序列的方差，Nx为检波器的数量，m＝-M/2，u(hi,t)为第i个检波器在地震域基于检波器的走时和检波器的深度确定的坐标。
在S208中，基于不同类型的地震波对应的预设慢度值范围，确定目标慢度谱中与每个类型的地震波对应的慢度点。
在实施中，不同类型的地震波对应的预设慢度值范围不同，例如，地震波可以包括上行P波、下行P波、上行S波和下行S波。其中，上行P波和上行S波对应的预设慢度值可以大于下行P波和下行S波对应的预设慢度值范围。
服务器可以根据预设慢度值范围，确定目标慢度谱中与每个类型的地震波对应的慢度点。
在S210中，基于每个类型的地震波对应的慢度点的慢度值，确定每个类型的地震波对应的慢度拟合直线参数。
在实施中，可以在获取到每个类型的地震波对应的慢度点的慢度值后，对这些慢度点进行线性拟合，并得到针对每个类型的地震波的慢度拟合直线参数。
在S212中，将目标地震波的慢度拟合直线参数，代入公式
Figure BDA0002548677570000101
得到与目标地震波对应的目标时移量。
其中，i为第i个检波器，Δti为第i个检波器的目标时移量，ti为第i个检波器的走时，Δh＝hi-hi-1，hi为第i个检波器的深度，an为慢度拟合直线参数中的第一慢度拟合直线参数，bn为慢度拟合直线参数中的第二慢度拟合直线参数。
在实施中，可以根据下述公式：
Δti＝ti-t1＝(i-1)*Δh*pn,
获得每个检波器的时移量，其中，pn为第n个类型的检波器的慢度值。
而在τ-p域内，不同检波器在不同截距时间的时移量满足下述公式：
Δtn(i,τ)＝(i-1)Δhpn＝(i-1)Δh(an+bnτ)＝ti-τ。
而将上述两个公式进行推导消除处理后，可以得到公式：
Figure BDA0002548677570000102
通过该公式可以看出，检波器的目标时移量是随着检波器的不同以及时间的不同而动态变化的。
在S214中，根据目标时移量，对目标地震波数据进行排齐处理，以使排齐处理后的目标地震波数据中的同相轴平行排齐。
在S216中，基于均值滤波算法，对排齐处理后的目标地震波数据进行过滤处理。
在实施中，可以根据均值滤波器：
Figure BDA0002548677570000111
对排齐处理后的目标地震波数据进行过滤处理，其中，i为第i个检波器，k为预设参数，且k小于检波器数量的一半，ti为第i个检波器的走时，d(ti)为第i个检波器在时间t的振幅。
在S218中，基于目标时移量，对过滤处理后的目标地震波数据进行反排齐处理，得到与目标地震波对应的目标波场。
在S220中，基于目标地震波数据和目标波场，确定第一地震波数据。
其中，第一地震波可以为不同类型的地震波中，与目标地震波的类型不同的任意一种类型的地震波。
在实施中，可以从目标地震波数据中减去目标波场，以得到第一地震波数据。
在S222中，基于第一地震波的慢度拟合直线参数，确定与第一地震波对应的第一时移量，并根据第一时移量，对第一地震波数据进行排齐处理，以使排齐处理后的第一地震波数据中的同相轴平行排齐。
在S224中，基于第一时移量，对过滤处理后的第一地震波数据进行反排齐处理，得到与第一地震波对应的第一波场。
在实施中，在得到第一波场后，可以继续从第一地震波数据中减去第一波场，得到与第二地震波对应的第二地震波数据，第二地震波为不同类型的地震波中，与目标地震波以及第一地震波的类型不同的任意一种类型的地震波，并根据上述获取第一地震波的第一波场相同的方法，确定第二地震波对应的第二波场。
可以基于上述获取第二地震波的方法，继续获取目标地震波数据中其他类型的地震波对应的波场，直到目标地震波数据中不包含有效波场，即已经获取了目标地震波数据中包含的所有类型的地震波对应的波场，实现了对目标地震波数据的波场分离，同时还可以输出目标地震波数据中包含的噪声数据。
此外，目标地震波的能量值可以大于第一地震波的能量值，第一地震波的能量值可以大于第二地震波的能量值，即可以根据地震波能量值的大小，依次从目标地震波数据中，分离出对应的波场。
本发明实施例提供一种波场确定方法，本发明实施例通过基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱，慢度谱确定方法是根据检波器的走时、检波器的深度、检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，确定慢度谱的方法，基于目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数，基于目标地震波的慢度拟合直线参数，确定与目标地震波对应的目标时移量，并根据目标时移量，对目标地震波数据进行排齐处理，以使排齐处理后的目标地震波数据中的同相轴平行排齐，目标地震波为不同类型的地震波中的任意一种类型的地震波，对排齐处理后的目标地震波数据进行过滤处理，并基于目标时移量，对过滤处理后的目标地震波数据进行反排齐处理，得到与目标地震波对应的目标波场。这样，通过τ-p域上的目标慢度谱，可以准确的获取不同类型的地震波对应的慢度拟合直线参数，在根据目标地震波的慢度拟合直线参数，对目标地震波数据进行排齐处理、过滤处理以及反排齐处理，以提高目标波场从目标地震波数据中分离出来的准确性，同时也可以抑制目标地震波数据中的噪声的影响，即可以提高目标地震波的目标波场的确定准确性。
实施例三
以上为本说明书实施例提供的波场确定方法，基于同样的思路，本说明书实施例还提供一种波场确定装置，如图3所示。
该波场确定装置包括：慢度谱确定模块301、参数确定模块302、第一确定模块303和第一波场确定模块304，其中：
慢度谱确定模块301，用于基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱，所述慢度谱确定方法是根据检波器的走时、所述检波器的深度、所述检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，确定慢度谱的方法；
参数确定模块302，用于基于所述目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数；
第一确定模块303，用于基于目标地震波的所述慢度拟合直线参数，确定与所述目标地震波对应的目标时移量，并根据所述目标时移量，对所述目标地震波数据进行排齐处理，以使排齐处理后的所述目标地震波数据中的同相轴平行排齐，所述目标地震波为所述不同类型的地震波中的任意一种类型的地震波；
第一波场确定模块304，用于对所述排齐处理后的目标地震波数据进行过滤处理，并基于所述目标时移量，对过滤处理后的所述目标地震波数据进行反排齐处理，得到与所述目标地震波对应的目标波场。
本发明实施例中，所述装置，还包括：
数据确定模块，用于基于所述目标地震波数据和所述目标波场，确定第一地震波数据，所述第一地震波为所述不同类型的地震波中，与所述目标地震波的类型不同的任意一种类型的地震波；
第二确定模块，用于基于第一地震波的所述慢度拟合直线参数，确定与所述第一地震波对应的第一时移量，并根据所述第一时移量，对所述第一地震波数据进行排齐处理，以使排齐处理后的所述第一地震波数据中的同相轴平行排齐；
第二波场确定模块，用于对所述排齐处理后的第一地震波数据进行过滤处理，并基于所述第一时移量，对过滤处理后的所述第一地震波数据进行反排齐处理，得到与所述第一地震波对应的第一波场。
本发明实施例中，所述慢度谱确定模块301，用于：
将每个所述检波器的走时、所述检波器的深度、所述检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，代入公式
Figure BDA0002548677570000131
得到所述目标地震波数据在τ-p域上的目标慢度谱，其中，v(p,τ)为所述目标地震波数据在τ-p域上的目标慢度谱，p为慢度值，τ为所述截距时间，M为所述预设时窗长度，hi为第i个所述检波器的深度，t为第i个所述检波器的走时，
Figure BDA0002548677570000132
为第i个所述检波器对应的叠加路径内振幅序列的方差，Nx为所述检波器的数量，m＝-M/2，u(hi,t)为第i个所述检波器在地震域基于所述检波器的走时和所述检波器的深度确定的坐标。
本发明实施例中，所述装置，还包括：
获取模块，用于基于所述目标地震波数据，获取每个所述检波器的截距时间、慢度值以及深度；
走时确定模块，用于将每个所述检波器的截距时间、慢度值以及深度代入公式
Figure BDA0002548677570000141
得到每个所述检波器的走时，其中，t为所述检波器的走时，τ为所述检波器的截距时间，p为所述检波器的慢度值，hi为所述检波器的深度。
本发明实施例中，所述参数确定模块302，用于：
基于所述不同类型的地震波对应的预设慢度值范围，确定所述目标慢度谱中与每个类型的所述地震波对应的慢度点；
基于所述每个类型的地震波对应的慢度点的慢度值，确定所述每个类型的地震波对应的慢度拟合直线参数。
本发明实施例中，所述第一确定模块303，用于：
将所述目标地震波的所述慢度拟合直线参数，代入公式
Figure BDA0002548677570000142
得到与所述目标地震波对应的目标时移量，其中，i为第i个所述检波器，Δti为第i个所述检波器的目标时移量，ti为第i个所述检波器的走时，Δh＝hi-hi-1，hi为第i个所述检波器的深度，an为所述慢度拟合直线参数中的第一慢度拟合直线参数，bn为所述慢度拟合直线参数中的第二慢度拟合直线参数。
本发明实施例中，所述第一波场确定模块304，用于：
基于均值滤波算法，对所述排齐处理后的目标地震波数据进行过滤处理。
本发明实施例提供一种波场确定方法，本发明实施例通过基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱，慢度谱确定方法是根据检波器的走时、检波器的深度、检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，确定慢度谱的方法，基于目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数，基于目标地震波的慢度拟合直线参数，确定与目标地震波对应的目标时移量，并根据目标时移量，对目标地震波数据进行排齐处理，以使排齐处理后的目标地震波数据中的同相轴平行排齐，目标地震波为不同类型的地震波中的任意一种类型的地震波，对排齐处理后的目标地震波数据进行过滤处理，并基于目标时移量，对过滤处理后的目标地震波数据进行反排齐处理，得到与目标地震波对应的目标波场。这样，通过τ-p域上的目标慢度谱，可以准确的获取不同类型的地震波对应的慢度拟合直线参数，在根据目标地震波的慢度拟合直线参数，对目标地震波数据进行排齐处理、过滤处理以及反排齐处理，以提高目标波场从目标地震波数据中分离出来的准确性，同时也可以抑制目标地震波数据中的噪声的影响，即可以提高目标地震波的目标波场的确定准确性。
实施例四
图4为实现本发明各个实施例的一种电子设备的硬件结构示意图，
该电子设备400包括但不限于：射频单元401、网络模块402、音频输出单元403、输入单元404、传感器405、显示单元406、用户输入单元407、接口单元404、存储器409、处理器410、以及电源411等部件。本领域技术人员可以理解，图4中示出的电子设备结构并不构成对电子设备的限定，电子设备可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。在本发明实施例中，电子设备包括但不限于手机、平板电脑、笔记本电脑、掌上电脑、车载终端、可穿戴设备、以及计步器等。
其中，处理器410，用于：基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱，所述慢度谱确定方法是根据检波器的走时、所述检波器的深度、所述检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，确定慢度谱的方法；基于所述目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数；基于目标地震波的所述慢度拟合直线参数，确定与所述目标地震波对应的目标时移量，并根据所述目标时移量，对所述目标地震波数据进行排齐处理，以使排齐处理后的所述目标地震波数据中的同相轴平行排齐，所述目标地震波为所述不同类型的地震波中的任意一种类型的地震波；对所述排齐处理后的目标地震波数据进行过滤处理，并基于所述目标时移量，对过滤处理后的所述目标地震波数据进行反排齐处理，得到与所述目标地震波对应的目标波场。
所述处理器410，还用于：基于所述目标地震波数据和所述目标波场，确定第一地震波数据，所述第一地震波为所述不同类型的地震波中，与所述目标地震波的类型不同的任意一种类型的地震波；基于第一地震波的所述慢度拟合直线参数，确定与所述第一地震波对应的第一时移量，并根据所述第一时移量，对所述第一地震波数据进行排齐处理，以使排齐处理后的所述第一地震波数据中的同相轴平行排齐；对所述排齐处理后的第一地震波数据进行过滤处理，并基于所述第一时移量，对过滤处理后的所述第一地震波数据进行反排齐处理，得到与所述第一地震波对应的第一波场。
所述处理器410，还用于：将每个所述检波器的走时、所述检波器的深度、所述检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，代入公式
Figure BDA0002548677570000161
得到所述目标地震波数据在τ-p域上的目标慢度谱，其中，v(p,τ)为所述目标地震波数据在τ-p域上的目标慢度谱，p为慢度值，τ为所述截距时间，M为所述预设时窗长度，hi为第i个所述检波器的深度，t为第i个所述检波器的走时，
Figure BDA0002548677570000162
为第i个所述检波器对应的叠加路径内振幅序列的方差，Nx为所述检波器的数量，m＝-M/2，u(hi,t)为第i个所述检波器在地震域基于所述检波器的走时和所述检波器的深度确定的坐标。
所述处理器410，还用于：基于所述目标地震波数据，获取每个所述检波器的截距时间、慢度值以及深度；将每个所述检波器的截距时间、慢度值以及深度代入公式
Figure BDA0002548677570000163
得到每个所述检波器的走时，其中，t为所述检波器的走时，τ为所述检波器的截距时间，p为所述检波器的慢度值，hi为所述检波器的深度。
所述处理器410，还用于：基于所述不同类型的地震波对应的预设慢度值范围，确定所述目标慢度谱中与每个类型的所述地震波对应的慢度点；基于所述每个类型的地震波对应的慢度点的慢度值，确定所述每个类型的地震波对应的慢度拟合直线参数。
所述处理器410，还用于：将所述目标地震波的所述慢度拟合直线参数，代入公式
Figure BDA0002548677570000164
得到与所述目标地震波对应的目标时移量，其中，i为第i个所述检波器，Δti为第i个所述检波器的目标时移量，ti为第i个所述检波器的走时，Δh＝hi-hi-1，hi为第i个所述检波器的深度，an为所述慢度拟合直线参数中的第一慢度拟合直线参数，bn为所述慢度拟合直线参数中的第二慢度拟合直线参数。
所述处理器410，还用于：基于均值滤波算法，对所述排齐处理后的目标地震波数据进行过滤处理。
本发明实施例提供一种电子设备，本发明实施例通过基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱，慢度谱确定方法是根据检波器的走时、检波器的深度、检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，确定慢度谱的方法，基于目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数，基于目标地震波的慢度拟合直线参数，确定与目标地震波对应的目标时移量，并根据目标时移量，对目标地震波数据进行排齐处理，以使排齐处理后的目标地震波数据中的同相轴平行排齐，目标地震波为不同类型的地震波中的任意一种类型的地震波，对排齐处理后的目标地震波数据进行过滤处理，并基于目标时移量，对过滤处理后的目标地震波数据进行反排齐处理，得到与目标地震波对应的目标波场。这样，通过τ-p域上的目标慢度谱，可以准确的获取不同类型的地震波对应的慢度拟合直线参数，在根据目标地震波的慢度拟合直线参数，对目标地震波数据进行排齐处理、过滤处理以及反排齐处理，以提高目标波场从目标地震波数据中分离出来的准确性，同时也可以抑制目标地震波数据中的噪声的影响，即可以提高目标地震波的目标波场的确定准确性。
应理解的是，本发明实施例中，射频单元401可用于收发信息或通话过程中，信号的接收和发送，具体的，将来自基站的下行数据接收后，给处理器410处理；另外，将上行的数据发送给基站。通常，射频单元401包括但不限于天线、至少一个放大器、收发信机、耦合器、低噪声放大器、双工器等。此外，射频单元401还可以通过无线通信系统与网络和其他设备通信。
电子设备通过网络模块402为用户提供了无线的宽带互联网访问，如帮助用户收发电子邮件、浏览网页和访问流式媒体等。
音频输出单元403可以将射频单元401或网络模块402接收的或者在存储器409中存储的音频数据转换成音频信号并且输出为声音。而且，音频输出单元403还可以提供与电子设备400执行的特定功能相关的音频输出(例如，呼叫信号接收声音、消息接收声音等等)。音频输出单元403包括扬声器、蜂鸣器以及受话器等。
输入单元404用于接收音频或视频信号。输入单元404可以包括图形处理器(Graphics Processing Unit，GPU)4041和麦克风4042，图形处理器4041对在视频捕获模式或图像捕获模式中由图像捕获装置(如摄像头)获得的静态图片或视频的图像数据进行处理。处理后的图像帧可以显示在显示单元406上。经图形处理器4041处理后的图像帧可以存储在存储器409(或其它存储介质)中或者经由射频单元401或网络模块402进行发送。麦克风4042可以接收声音，并且能够将这样的声音处理为音频数据。处理后的音频数据可以在电话通话模式的情况下转换为可经由射频单元401发送到移动通信基站的格式输出。
电子设备400还包括至少一种传感器405，比如光传感器、运动传感器以及其他传感器。具体地，光传感器包括环境光传感器及接近传感器，其中，环境光传感器可根据环境光线的明暗来调节显示面板4061的亮度，接近传感器可在电子设备400移动到耳边时，关闭显示面板4061和/或背光。作为运动传感器的一种，加速计传感器可检测各个方向上(一般为三轴)加速度的大小，静止时可检测出重力的大小及方向，可用于识别电子设备姿态(比如横竖屏切换、相关游戏、磁力计姿态校准)、振动识别相关功能(比如计步器、敲击)等；传感器405还可以包括指纹传感器、压力传感器、虹膜传感器、分子传感器、陀螺仪、气压计、湿度计、温度计、红外线传感器等，在此不再赘述。
显示单元406用于显示由用户输入的信息或提供给用户的信息。显示单元406可包括显示面板4061，可以采用液晶显示器(Liquid Crystal Display，LCD)、有机发光二极管(Organic Light-Emitting Diode,OLED)等形式来配置显示面板4061。
用户输入单元407可用于接收输入的数字或字符信息，以及产生与电子设备的用户设置以及功能控制有关的键信号输入。具体地，用户输入单元407包括触控面板4071以及其他输入设备4072。触控面板401，也称为触摸屏，可收集用户在其上或附近的触摸操作(比如用户使用手指、触笔等任何适合的物体或附件在触控面板4071上或在触控面板4071附近的操作)。触控面板4071可包括触摸检测装置和触摸控制器两个部分。其中，触摸检测装置检测用户的触摸方位，并检测触摸操作带来的信号，将信号传送给触摸控制器；触摸控制器从触摸检测装置上接收触摸信息，并将它转换成触点坐标，再送给处理器410，接收处理器410发来的命令并加以执行。此外，可以采用电阻式、电容式、红外线以及表面声波等多种类型实现触控面板4071。除了触控面板4071，用户输入单元407还可以包括其他输入设备4072。具体地，其他输入设备4072可以包括但不限于物理键盘、功能键(比如音量控制按键、开关按键等)、轨迹球、鼠标、操作杆，在此不再赘述。
进一步的，触控面板4071可覆盖在显示面板4061上，当触控面板4071检测到在其上或附近的触摸操作后，传送给处理器410以确定触摸事件的类型，随后处理器410根据触摸事件的类型在显示面板4061上提供相应的视觉输出。虽然在图4中，触控面板4071与显示面板4061是作为两个独立的部件来实现电子设备的输入和输出功能，但是在某些实施例中，可以将触控面板4071与显示面板4061集成而实现电子设备的输入和输出功能，具体此处不做限定。
接口单元408为外部装置与电子设备400连接的接口。例如，外部装置可以包括有线或无线头戴式耳机端口、外部电源(或电池充电器)端口、有线或无线数据端口、存储卡端口、用于连接具有识别模块的装置的端口、音频输入/输出(I/O)端口、视频I/O端口、耳机端口等等。接口单元408可以用于接收来自外部装置的输入(例如，数据信息、电力等等)并且将接收到的输入传输到电子设备400内的一个或多个元件或者可以用于在电子设备400和外部装置之间传输数据。
存储器409可用于存储软件程序以及各种数据。存储器409可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序(比如声音播放功能、图像播放功能等)等；存储数据区可存储根据手机的使用所创建的数据(比如音频数据、电话本等)等。此外，存储器409可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。
处理器410是电子设备的控制中心，利用各种接口和线路连接整个电子设备的各个部分，通过运行或执行存储在存储器409内的软件程序和/或模块，以及调用存储在存储器409内的数据，执行电子设备的各种功能和处理数据，从而对电子设备进行整体监控。处理器410可包括一个或多个处理单元；优选的，处理器410可集成应用处理器和调制解调处理器，其中，应用处理器主要处理操作系统、用户界面和应用程序等，调制解调处理器主要处理无线通信。可以理解的是，上述调制解调处理器也可以不集成到处理器410中。
电子设备400还可以包括给各个部件供电的电源411(比如电池)，优选的，电源411可以通过电源管理系统与处理器410逻辑相连，从而通过电源管理系统实现管理充电、放电、以及功耗管理等功能。
优选的，本发明实施例还提供一种电子设备，包括处理器410，存储器409，存储在存储器409上并可在所述处理器410上运行的计算机程序，该计算机程序被处理器410执行时实现上述波场确定方法实施例的各个过程，且能达到相同的技术效果，为避免重复，这里不再赘述。
实施例五
本发明实施例还提供一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述波场确定方法实施例的各个过程，且能达到相同的技术效果，为避免重复，这里不再赘述。其中，所述的计算机可读存储介质，如只读存储器(Read-Only Memory，简称ROM)、随机存取存储器(Random Access Memory，简称RAM)、磁碟或者光盘等。
本发明实施例提供一种计算机可读存储介质，通过基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱，慢度谱确定方法是根据检波器的走时、检波器的深度、检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，确定慢度谱的方法，基于目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数，基于目标地震波的慢度拟合直线参数，确定与目标地震波对应的目标时移量，并根据目标时移量，对目标地震波数据进行排齐处理，以使排齐处理后的目标地震波数据中的同相轴平行排齐，目标地震波为不同类型的地震波中的任意一种类型的地震波，对排齐处理后的目标地震波数据进行过滤处理，并基于目标时移量，对过滤处理后的目标地震波数据进行反排齐处理，得到与目标地震波对应的目标波场。这样，通过τ-p域上的目标慢度谱，可以准确的获取不同类型的地震波对应的慢度拟合直线参数，在根据目标地震波的慢度拟合直线参数，对目标地震波数据进行排齐处理、过滤处理以及反排齐处理，以提高目标波场从目标地震波数据中分离出来的准确性，同时也可以抑制目标地震波数据中的噪声的影响，即可以提高目标地震波的目标波场的确定准确性。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本发明的实施例可提供为方法、系统或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
以上所述仅为本发明的实施例而已，并不用于限制本发明。对于本领域技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。

Claims (10)

1.一种波场确定的方法，其特征在于，包括：
基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱，所述慢度谱确定方法是根据检波器的走时、所述检波器的深度、所述检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，确定慢度谱的方法；
基于所述目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数；
基于目标地震波的所述慢度拟合直线参数，确定与所述目标地震波对应的目标时移量，并根据所述目标时移量，对所述目标地震波数据进行排齐处理，以使排齐处理后的所述目标地震波数据中的同相轴平行排齐，所述目标地震波为所述不同类型的地震波中的任意一种类型的地震波；
对所述排齐处理后的目标地震波数据进行过滤处理，并基于所述目标时移量，对过滤处理后的所述目标地震波数据进行反排齐处理，得到与所述目标地震波对应的目标波场。
2.根据权利要求1所述的方法，其特征在于，在所述基于所述目标时移量，对过滤处理后的所述目标地震波数据进行反排齐处理，得到与所述目标地震波对应的目标波场之后，还包括：
基于所述目标地震波数据和所述目标波场，确定第一地震波数据，所述第一地震波为所述不同类型的地震波中，与所述目标地震波的类型不同的任意一种类型的地震波；
基于第一地震波的所述慢度拟合直线参数，确定与所述第一地震波对应的第一时移量，并根据所述第一时移量，对所述第一地震波数据进行排齐处理，以使排齐处理后的所述第一地震波数据中的同相轴平行排齐；
对所述排齐处理后的第一地震波数据进行过滤处理，并基于所述第一时移量，对过滤处理后的所述第一地震波数据进行反排齐处理，得到与所述第一地震波对应的第一波场。
3.根据权利要求2所述的方法，其特征在于，所述基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱，包括：
将每个所述检波器的走时、所述检波器的深度、所述检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，代入公式
Figure FDA0002548677560000021
得到所述目标地震波数据在τ-p域上的目标慢度谱，其中，v(p,τ)为所述目标地震波数据在τ-p域上的目标慢度谱，p为慢度值，τ为所述截距时间，M为所述预设时窗长度，hi为第i个所述检波器的深度，t为第i个所述检波器的走时，
Figure FDA0002548677560000022
为第i个所述检波器对应的叠加路径内振幅序列的方差，Nx为所述检波器的数量，m＝-M/2，u(hi,t)为第i个所述检波器在地震域基于所述检波器的走时和所述检波器的深度确定的坐标。
4.根据权利要求3所述的方法，其特征在于，在所述基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱之前，还包括：
基于所述目标地震波数据，获取每个所述检波器的截距时间、慢度值以及深度；
将每个所述检波器的截距时间、慢度值以及深度代入公式
Figure FDA0002548677560000023
得到每个所述检波器的走时，其中，t为所述检波器的走时，τ为所述检波器的截距时间，p为所述检波器的慢度值，hi为所述检波器的深度。
5.根据权利要求4所述的方法，其特征在于，所述基于所述目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数，包括：
基于所述不同类型的地震波对应的预设慢度值范围，确定所述目标慢度谱中与每个类型的所述地震波对应的慢度点；
基于所述每个类型的地震波对应的慢度点的慢度值，确定所述每个类型的地震波对应的慢度拟合直线参数。
6.根据权利要求5所述的方法，其特征在于，所述基于目标地震波的所述慢度拟合直线参数，确定与所述目标地震波对应的目标时移量，包括：
将所述目标地震波的所述慢度拟合直线参数，代入公式
Figure FDA0002548677560000031
得到与所述目标地震波对应的目标时移量，其中，i为第i个所述检波器，Δti为第i个所述检波器的目标时移量，ti为第i个所述检波器的走时，Δh＝hi-hi-1，hi为第i个所述检波器的深度，an为所述慢度拟合直线参数中的第一慢度拟合直线参数，bn为所述慢度拟合直线参数中的第二慢度拟合直线参数。
7.根据权利要求1所述的方法，其特征在于，所述对所述排齐处理后的目标地震波数据进行过滤处理，包括：
基于均值滤波算法，对所述排齐处理后的目标地震波数据进行过滤处理。
8.一种波场确定的装置，其特征在于，包括：
慢度谱确定模块，用于基于预先构建的慢度谱确定方法，确定目标地震波数据在τ-p域上的目标慢度谱，所述慢度谱确定方法是根据检波器的走时、所述检波器的深度、所述检波器对应的叠加路径内振幅序列的方差、截距时间以及预设时窗长度，确定慢度谱的方法；
参数确定模块，用于基于所述目标慢度谱，确定不同类型的地震波对应的慢度拟合直线参数；
第一确定模块，用于基于目标地震波的所述慢度拟合直线参数，确定与所述目标地震波对应的目标时移量，并根据所述目标时移量，对所述目标地震波数据进行排齐处理，以使排齐处理后的所述目标地震波数据中的同相轴平行排齐，所述目标地震波为所述不同类型的地震波中的任意一种类型的地震波；
第一波场确定模块，用于对所述排齐处理后的目标地震波数据进行过滤处理，并基于所述目标时移量，对过滤处理后的所述目标地震波数据进行反排齐处理，得到与所述目标地震波对应的目标波场。
9.根据权利要求8所述的装置，其特征在于，所述装置，还包括：
数据确定模块，用于基于所述目标地震波数据和所述目标波场，确定第一地震波数据，所述第一地震波为所述不同类型的地震波中，与所述目标地震波的类型不同的任意一种类型的地震波；
第二确定模块，用于基于第一地震波的所述慢度拟合直线参数，确定与所述第一地震波对应的第一时移量，并根据所述第一时移量，对所述第一地震波数据进行排齐处理，以使排齐处理后的所述第一地震波数据中的同相轴平行排齐；
第二波场确定模块，用于对所述排齐处理后的第一地震波数据进行过滤处理，并基于所述第一时移量，对过滤处理后的所述第一地震波数据进行反排齐处理，得到与所述第一地震波对应的第一波场。
10.一种电子设备，其特征在于，包括处理器、存储器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述计算机程序被所述处理器执行时实现如权利要求1至7中任一项所述的波场确定方法的步骤。
CN202010568542.4A 2020-06-19 2020-06-19 波场确定方法与装置 Active CN111897001B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010568542.4A CN111897001B (zh) 2020-06-19 2020-06-19 波场确定方法与装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010568542.4A CN111897001B (zh) 2020-06-19 2020-06-19 波场确定方法与装置

Publications (2)

Publication Number Publication Date
CN111897001A true CN111897001A (zh) 2020-11-06
CN111897001B CN111897001B (zh) 2021-04-13

Family

ID=73206897

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010568542.4A Active CN111897001B (zh) 2020-06-19 2020-06-19 波场确定方法与装置

Country Status (1)

Country Link
CN (1) CN111897001B (zh)

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107144880B (zh) * 2017-05-12 2019-07-05 招商局重庆交通科研设计院有限公司 一种地震波波场分离方法

Also Published As

Publication number Publication date
CN111897001B (zh) 2021-04-13

Similar Documents

Publication Publication Date Title
CN109558512B (zh) 一种基于音频的个性化推荐方法、装置和移动终端
CN107576986B (zh) 一种确定海底地震反射波旅行时的方法和装置
CN110674112A (zh) 一种数据查询的方法、装置及电子设备
CN110988993B (zh) 一种偏移成像方法、装置及电子设备
CN107505653A (zh) 一种确定叠前偏移时间结果的方法和装置
CN108828613B (zh) 一种去除噪点的方法和危化品仓储激光扫描装置
CN110389380B (zh) 地震剖面的同相轴自动拾取方法、装置及存储介质
CN110808872A (zh) 一种流量实验的实现方法、装置及电子设备
CN110348275A (zh) 手势识别方法、装置、智能设备及计算机可读存储介质
CN107784298B (zh) 一种识别方法及装置
CN112231144A (zh) 一种数据处理方法、装置及电子设备
CN112558144B (zh) 基于零偏vsp数据确定裂缝方位的方法、装置及电子设备
CN111897001B (zh) 波场确定方法与装置
CN111964773B (zh) 一种基于多个交叠线性调频组的振动信号提取方法及装置
CN115935816A (zh) 钻井参数确定方法、装置、设备及存储介质
CN111505717B (zh) 一种三维vsp共面元成像方法与装置
CN113341456B (zh) 一种基于非规则网格的地震偏移方法、装置及电子设备
CN111666421A (zh) 数据处理方法、装置及电子设备
CN111551987A (zh) 利用变偏移距vsp数据确定岩石弹性性质的方法与装置
CN111601302A (zh) 一种数据传输方法及电子设备
CN111238495A (zh) 一种定位车辆的方法和终端设备
CN110658388A (zh) 一种接地电阻的检测方法、装置及电子设备
CN112698391B (zh) 砂体刻画方法、装置、设备及存储介质
CN113535693B (zh) 一种移动平台的数据真值确定方法、装置及电子设备
CN110058306B (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230613

Address after: Room 522, Zhiyuan Building, Chengdu University of Technology, No. 1, Erxianqiao East Third Road, Chenghua District, Chengdu 610059, Sichuan

Patentee after: Chengdu Wushi Technology Co.,Ltd.

Address before: Three road 610051 Sichuan city of Chengdu Province, No. 1 East

Patentee before: Chengdu University of Technology