CN116736383A - 地震波速度模型的更新方法、装置、电子设备及存储介质 - Google Patents
地震波速度模型的更新方法、装置、电子设备及存储介质 Download PDFInfo
- Publication number
- CN116736383A CN116736383A CN202310730120.6A CN202310730120A CN116736383A CN 116736383 A CN116736383 A CN 116736383A CN 202310730120 A CN202310730120 A CN 202310730120A CN 116736383 A CN116736383 A CN 116736383A
- Authority
- CN
- China
- Prior art keywords
- seismic
- velocity model
- time
- updating
- arrival time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 82
- 238000003860 storage Methods 0.000 title claims abstract description 15
- 238000005520 cutting process Methods 0.000 claims abstract description 44
- 238000012937 correction Methods 0.000 claims abstract description 16
- 238000003384 imaging method Methods 0.000 claims description 15
- 238000013508 migration Methods 0.000 claims description 14
- 230000005012 migration Effects 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 13
- 238000001914 filtration Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 12
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 claims description 9
- 238000010586 diagram Methods 0.000 claims description 7
- 238000007781 pre-processing Methods 0.000 claims description 7
- 239000002360 explosive Substances 0.000 abstract description 11
- 238000009412 basement excavation Methods 0.000 abstract description 3
- 238000012544 monitoring process Methods 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 7
- 239000003245 coal Substances 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000005065 mining Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000004807 localization Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000005019 pattern of movement Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
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
本发明提供的一种地震波速度模型的更新方法、装置、电子设备及存储介质，通过获取检波器模块采集的微震事件，所述微震事件包括微震事件的发生坐标和微震事件的截割时间；根据所述微震事件，获取所述微震事件的P波到时和S波到时；根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，更新预先构建的地震波速度模型；通过更新后的地震波速度模型，预测地震波的旅行时。本发明无需在井下使用***激发产生地震波来校正速度模型，能够通过挖采工作中产生的固有微震事件，对预先构建的地震波速度模型进行校正，是一种可代替***震源且安全环保、操作简单的速度模型校正方法。
Description
技术领域
本发明涉及微地震技术领域，尤其涉及一种地震波速度模型的更新方法、装置、电子设备及存储介质。
背景技术
目前矿井地层井下定向孔压裂多在裸眼孔中进行后退式压裂，并且没有射孔，所以无法利用射孔信号来进行速度模型的校正；虽然可以通过***激发，基于标定炮来进行速度模型校正。但是，随着煤矿对***等火攻品管控的更加严格，在压裂区域多次放炮，不仅污染井下环境、成本高，而且还存在安全隐患，不符合绿色矿山建设的要求。
综上所述，现有技术中存在的问题亟需得到解决。
发明内容
本发明提供一种地震波速度模型的更新方法、装置、电子设备及存储介质，用以解决现有技术中基于标定炮来进行速度模型校正的缺陷，实现简单、绿色而且成本很低的速度模型校正。
本发明提供一种地震波速度模型的更新方法，包括：
获取检波器模块采集的微震事件，所述微震事件包括微震事件的发生坐标和微震事件的截割时间；
根据所述微震事件，获取所述微震事件的P波到时和S波到时；
根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，更新预先构建的地震波速度模型；
通过更新后的地震波速度模型，预测地震波的旅行时。
根据本发明提供的一种地震波速度模型的更新方法，根据所述发生坐标和所述截割时间，获取所述微震事件的P波到时和S波到时，具体包括：
对所述微震事件进行数据预处理，所述预处理包括去噪处理、滤波处理和增益校正处理；
对进行预处理后的微震事件进行初至拾取，以获取所述微震事件的P波到时和S波到时。
根据本发明提供的一种地震波速度模型的更新方法，根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，更新预先构建的地震波速度模型，具体包括：
将所述发生坐标输入至所述速度模型中，得到地震波的预测旅行时；
根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，得到地震波的真实旅行时；
根据所述预测旅行时和所述真实旅行时，对所述速度模型进行更新。
根据本发明提供的一种地震波速度模型的更新方法，根据所述预测旅行时和所述真实旅行时，对所述速度模型进行更新，具体包括：
根据所述预测旅行时和所述真实旅行时，通过Levenberg-Marquardt算法对所述速度模型进行更新。
根据本发明提供的一种地震波速度模型的更新方法，获取检波器模块采集的微震事件，具体包括：
按照预设的时间间隔获取检波器模块采集的微震事件。
根据本发明提供的一种地震波速度模型的更新方法，在通过更新后的地震波速度模型，预测震源到所述检波器模块的旅行时后，所述方法还包括：
根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，进行叠前偏移成像，以获得所述微震事件的偏移成像结果图。
本发明还提供一种地震波速度模型的更新装置，包括：
获取单元，用于获取检波器模块采集的微震事件，所述微震事件包括微震事件的发生坐标和微震事件的截割时间；
计算单元，用于根据所述微震事件，获取所述微震事件的P波到时和S波到时；
更新单元，用于根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，更新预先构建的地震波速度模型；
预测单元，用于通过更新后的地震波速度模型，预测地震波的旅行时。
本发明还提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上述任一种所述地震波速度模型的更新方法。
本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如上述任一种所述地震波速度模型的更新方法。
本发明还提供一种计算机程序产品，包括计算机程序，所述计算机程序被处理器执行时实现如上述任一种所述地震波速度模型的更新方法。
本发明提供的一种地震波速度模型的更新方法、装置、电子设备及存储介质，通过获取检波器模块采集的微震事件，所述微震事件包括微震事件的发生坐标和微震事件的截割时间；根据所述微震事件，获取所述微震事件的P波到时和S波到时；根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，更新预先构建的地震波速度模型；通过更新后的地震波速度模型，预测地震波的旅行时。本发明无需在井下使用***激发产生地震波来校正速度模型，能够通过挖采工作中产生的固有微震事件，对预先构建的地震波速度模型进行校正，是一种可代替***震源且安全环保、操作简单的速度模型校正方法。
附图说明
为了更清楚地说明本发明或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明提供的地震波速度模型的更新方法的流程示意图；
图2是本发明提供的检波器模块设置示意图；
图3是本发明提供的地震波速度模型的更新装置的结构示意图；
图4是本发明提供的电子设备的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
目前矿井地层井下定向孔压裂多在裸眼孔中进行后退式压裂，并且没有射孔，所以无法利用射孔信号来进行速度模型的校正；虽然可以通过***激发，基于标定炮来进行速度模型校正。但是，随着煤矿对***等火攻品管控的更加严格，在压裂区域多次放炮，不仅污染井下环境、成本高，而且还存在安全隐患，不符合绿色矿山建设的要求。
为了解决上述问题，本发明提出了一种地震波速度模型的更新方法，该方法包括但不限于步骤110-140：
步骤110、获取检波器模块采集的微震事件，所述微震事件包括微震事件的发生坐标和微震事件的截割时间。
在步骤110中，需要通过检波器模块采集待检测区域的微震事件，以获取微震事件的发生坐标和微震事件的截割时间。
检波器是一种用于检测和测量信号的设备或装置。在地震监测中，检波器也被称为地震仪，用于测量地震波的振动。检波器模块是设置在待检测区域的一组检波器，如图2所示。通过设置检波器模块，可以采集待检测区域的微震事件。
微震事件是指地球上非常小幅度的地震活动。它们通常由地壳中的岩石断裂或移动引起，但能量释放的规模非常小。微震事件通常在地震监测网络中被记录下来。这些事件在地震研究中具有重要意义，可以帮助了解地壳的运动和地震活动的模式。而在本发明中的微震事件指的是采煤机截割产生的微震事件，采煤机截割产生的微震事件对于水力压裂微震监测来说是一种干扰信号。但是，本发明可以将采煤机截割产生的固有微震事件利用起来，无需井下使用***激发产生地震波来校正速度模型。
步骤120、根据所述微震事件，获取所述微震事件的P波到时和S波到时。
在步骤120中，需要根据获取到的微震事件，得到P波到时和S波到时。微震事件的P波到时和S波到时是测量地震波传播时间的重要参数。P波(纵波)是一种沿着介质传播的压缩性波动，而S波(横波)是一种沿着介质传播的剪切性波动。P波到时指的是从地震事件发生到P波首次到达地震检测点的时间。它是地震波在地壳中传播的初始阶段，通常是最先被记录到的波。P波的到时信息可以用于确定地震的发生时间以及地震震源的大致位置。S波到时指的是从地震事件发生到S波首次到达地震检测点的时间。S波相对于P波而言传播速度较慢，因此在地震波到达的顺序中稍晚一些。S波的到时信息可以帮助确定地震的震级和震源的更准确位置。通过测量多个地震检测点上的P波和S波到时，可以进行三角定位或其他地震定位算法，以确定地震的震源位置和深度。这对于地震研究和地震监测非常重要，并可用于地震预警***的运作。
步骤130、根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，更新预先构建的地震波速度模型。
在步骤130中，需要将通过步骤110和步骤120测得的真实旅行时，与预先构建的地震波速度模型预测的预测旅行时进行对比，并采用Levenberg-Marquard(LM)算法不断迭代更新速度模型，直到观测的旅行时和正演的旅行时之间的差异最小。
具体地，本发明射线追踪正演方法得到震源到各检波器的旅行时，即采煤机截割产生的固有微震事件到各检波器的旅行时。射线追踪正演方法是一种常用的地震波传播模拟方法，可以用来计算地震波从震源到各个检波器的旅行时间。具体步骤包括：
构建速度模型：首先需要建立一个地震波速度模型，描述地球内部的速度分布。这个模型通常是基于地震观测数据或地球物理勘探数据获得的。
发射射线：从震源开始，发射一束射线(即地震波路径)向各个方向传播。射线的路径会根据地震波速度模型中不同介质的速度变化而弯曲、折射或反射。
计算旅行时间：根据射线路径和速度模型，计算每条射线的旅行时间。旅行时间取决于射线的长度和经过的不同介质的速度。
汇总旅行时间：对于每个检波器，将所有射线到达该检波器的旅行时间进行汇总，得到震源到各检波器的总旅行时间。
通过射线追踪正演方法得到的旅行时间信息可以用于地震定位、地震波形模拟、地震成像和地震速度模型验证等应用。这种方法在地震学研究和地震监测中具有广泛的应用价值。
步骤140、通过更新后的地震波速度模型，预测地震波的旅行时。
在步骤140中，可以通过更新后的地震波速度模型，预测地震波的旅行时。可以理解的是，更新后的地震波速度模型比未更新的地震波速度模型更加准确。
本发明提供的一种地震波速度模型的更新方法、装置、电子设备及存储介质，通过获取检波器模块采集的微震事件，所述微震事件包括微震事件的发生坐标和微震事件的截割时间；根据所述微震事件，获取所述微震事件的P波到时和S波到时；根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，更新预先构建的地震波速度模型；通过更新后的地震波速度模型，预测地震波的旅行时。本发明无需在井下使用***激发产生地震波来校正速度模型，能够通过挖采工作中产生的固有微震事件，对预先构建的地震波速度模型进行校正，是一种可代替***震源且安全环保、操作简单的速度模型校正方法。
根据本发明提供的一种地震波速度模型的更新方法，根据所述发生坐标和所述截割时间，获取所述微震事件的P波到时和S波到时，具体包括：
对所述微震事件进行数据预处理，所述预处理包括去噪处理、滤波处理和增益校正处理；
对进行预处理后的微震事件进行初至拾取，以获取所述微震事件的P波到时和S波到时。
本实施例中，为了获取更准确的微震事件的P波到时和S波到时，需要先对进行微震事件数据预处理，以减少噪音的影响，然后进行初至拾取，以确定地震波的P波和S波的到时。
对微震事件进行数据预处理通常包括去噪处理、滤波处理和增益校正处理。这些步骤有助于提取微震信号并减少噪音的影响。
去噪处理：去噪处理的目标是减少背景噪音对微震信号的干扰。常见的去噪方法包括中值滤波、小波去噪和频率域滤波等。这些方法可以根据具体情况选择最适合的方法来降低噪音水平。
滤波处理：滤波处理用于选择性地增强微震信号的频率成分并减少其他频率成分的影响。一般会采用带通滤波器来滤除低频和高频噪声，以突出微震事件的频率特征。常用的滤波方法包括低通滤波、高通滤波和带通滤波。
增益校正处理：增益校正处理旨在调整微震信号的振幅以实现更好的可视化效果和更准确的初至拾取。常用的增益校正方法包括线性增益校正和对数增益校正。这些方法可以根据信号的动态范围进行调整，以突出微震信号并保持合适的对比度。
对预处理后的微震事件进行初至拾取可以通过人工或自动的方式来进行。
人工初至拾取：在人工初至拾取中，地震学专家会仔细观察预处理后的微震信号，并标记出P波和S波的到时。这通常涉及对信号的到达时间进行目测，以确定地震波的P波和S波的到时。
自动初至拾取：自动初至拾取是使用计算机算法来自动识别和提取微震事件的到时。这可以通过应用信号处理和机器学习技术来实现。自动初至拾取方法可以提高效率，并在大量数据中进行快速的初至提取。
初至拾取的目标是确定微震事件的P波和S波到时，从而提供有关地震的时间和位置信息。这对于地震监测和地震学研究具有重要意义。
根据本发明提供的一种地震波速度模型的更新方法，根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，更新预先构建的地震波速度模型，具体包括：
将所述发生坐标输入至所述速度模型中，得到地震波的预测旅行时；
根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，得到地震波的真实旅行时；
根据所述预测旅行时和所述真实旅行时，对所述速度模型进行更新。
在本实施例中，更新地震波速度模型的方法可以按照以下步骤进行：
将发生坐标输入至速度模型中：根据发生坐标，将其作为输入参数输入到预先构建的地震波速度模型中。这可以用来计算地震波在该位置上的预测旅行时。
根据发生坐标、截割时间、P波到时和S波到时获取真实旅行时：利用发生坐标、截割时间、P波到时和S波到时的数据，可以计算出地震波的真实旅行时。这些数据可以基于地震观测数据、地震台网或其他相关的测量方法得到。
更新速度模型：通过对比预测旅行时和真实旅行时，可以评估速度模型的准确性和误差。根据比较结果，可以进行速度模型的更新。一种常见的方法是使用反演算法，通过最小化预测旅行时和真实旅行时之间的差异来调整速度模型中的参数。
这个更新方法允许根据实际观测数据来改进地震波速度模型，从而更准确地描述地球内部的速度结构。这对于地震学研究、地震监测和地震预测具有重要意义，并有助于提高地震定位和地震波形模拟的准确性。
根据本发明提供的一种地震波速度模型的更新方法，根据所述预测旅行时和所述真实旅行时，对所述速度模型进行更新，具体包括：
根据所述预测旅行时和所述真实旅行时，通过Levenberg-Marquardt算法对所述速度模型进行更新。
本实施例中，地震波速度模型的更新方法可以使用Levenberg-Marquardt算法对速度模型进行更新。Levenberg-Marquardt算法是一种非线性最小二乘优化算法，常用于参数估计和模型拟合。具体更新步骤如下：
首先，根据预测旅行时和真实旅行时之间的差异，构建一个目标函数，该目标函数用于衡量速度模型的拟合程度。然后，使用Levenberg-Marquardt算法来最小化目标函数。该算法将在速度模型的参数空间中搜索最优解，以使目标函数达到最小值。
在Levenberg-Marquardt算法的迭代过程中，会根据预测旅行时和真实旅行时的差异来更新速度模型的参数。通过迭代优化，可以逐步调整速度模型，使其与实际观测数据更好地吻合。
Levenberg-Marquardt算法的优点是可以处理非线性问题，并且具有快速收敛的特性。它在地震波速度模型的更新和反演中广泛应用，用于改进地球内部的速度结构描述。使用Levenberg-Marquardt算法对速度模型进行更新可以提高地震波传播模拟的准确性，有助于改善地震定位、地震波形模拟和地震预测等应用的结果。
根据本发明提供的一种地震波速度模型的更新方法，获取检波器模块采集的微震事件，具体包括：
按照预设的时间间隔获取检波器模块采集的微震事件。
在本实施例中，按照预设的时间间隔获取检波器模块采集的微震事件。
具体地，确定时间间隔：根据需要和实际情况，设定一个合适的时间间隔，表示想要获取微震事件的频率。时间间隔可以根据研究需求、监测目的或数据采集***的要求来确定。
设置数据采集***：将检波器模块配置为按照预设的时间间隔进行采集。这可以通过相应的软件或硬件设置来完成。确保检波器模块在每个时间间隔内进行数据采集。
获取微震事件数据：在预设的时间间隔结束后，从检波器模块中获取采集到的微震事件数据。这些数据通常以数字形式存储，可以包括微震信号的振幅、时间戳和其他相关参数。
根据本发明提供的一种地震波速度模型的更新方法，在通过更新后的地震波速度模型，预测震源到所述检波器模块的旅行时后，所述方法还包括：
根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，进行叠前偏移成像，以获得所述微震事件的偏移成像结果图。
本实施例中，在预测震源到检波器模块的旅行时后，可以进行叠前偏移成像以获得微震事件的偏移成像结果图。具体步骤如下：
准备叠前偏移数据：将预测的震源到检波器模块的旅行时数据作为输入，并结合发生坐标、截割时间、P波到时和S波到时等信息，准备进行叠前偏移的数据。
进行叠前偏移成像：根据叠前偏移算法，对准备好的数据进行处理，生成微震事件的偏移成像结果图。叠前偏移是一种常用的地震成像技术，通过将数据在地震波速度模型中反向传播和叠加，来还原地下的地震反射界面和体波散射。
分析偏移成像结果：通过分析偏移成像结果图，可以获得关于地下结构的信息，例如地下反射界面的位置、形态和反射强度等。这有助于地质解释、地下构造分析和勘探目标识别等应用。
通过进行叠前偏移成像，可以将微震事件的数据转化为地下结构的图像，提供更直观和详细的地质信息。这可以帮助地震学研究人员和勘探地球物理学家深入理解地下构造，进一步完善地震波速度模型，并促进地震监测和地下资源勘探的发展。
参考图3，下面对本发明提供的地震波速度模型的更新装置进行描述，下文描述的地震波速度模型的更新装置与上文描述的地震波速度模型的更新方法可相互对应参照。
本发明还提供一种地震波速度模型的更新装置，包括：
获取单元310，用于获取检波器模块采集的微震事件，所述微震事件包括微震事件的发生坐标和微震事件的截割时间；
计算单元320，用于根据所述微震事件，获取所述微震事件的P波到时和S波到时；
更新单元330，用于根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，更新预先构建的地震波速度模型；
预测单元340，用于通过更新后的地震波速度模型，预测地震波的旅行时。
图4示例了一种电子设备的实体结构示意图，如图4所示，该电子设备可以包括：处理器(processor)410、通信接口(Communications Interface)420、存储器(memory)430和通信总线440，其中，处理器410，通信接口420，存储器430通过通信总线440完成相互间的通信。处理器410可以调用存储器430中的逻辑指令，以执行地震波速度模型的更新方法，该方法包括：
获取检波器模块采集的微震事件，所述微震事件包括微震事件的发生坐标和微震事件的截割时间；
根据所述微震事件，获取所述微震事件的P波到时和S波到时；
根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，更新预先构建的地震波速度模型；
通过更新后的地震波速度模型，预测地震波的旅行时。
此外，上述的存储器430中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
另一方面，本发明还提供一种计算机程序产品，所述计算机程序产品包括计算机程序，计算机程序可存储在非暂态计算机可读存储介质上，所述计算机程序被处理器执行时，计算机能够执行上述各方法所提供的地震波速度模型的更新方法，该方法包括：
获取检波器模块采集的微震事件，所述微震事件包括微震事件的发生坐标和微震事件的截割时间；
根据所述微震事件，获取所述微震事件的P波到时和S波到时；
根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，更新预先构建的地震波速度模型；
通过更新后的地震波速度模型，预测地震波的旅行时。
又一方面，本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各方法提供的地震波速度模型的更新方法，该方法包括：
获取检波器模块采集的微震事件，所述微震事件包括微震事件的发生坐标和微震事件的截割时间；
根据所述微震事件，获取所述微震事件的P波到时和S波到时；
根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，更新预先构建的地震波速度模型；
通过更新后的地震波速度模型，预测地震波的旅行时。
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种地震波速度模型的更新方法，其特征在于，包括：
获取检波器模块采集的微震事件，所述微震事件包括微震事件的发生坐标和微震事件的截割时间；
根据所述微震事件，获取所述微震事件的P波到时和S波到时；
根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，更新预先构建的地震波速度模型；
通过更新后的地震波速度模型，预测地震波的旅行时。
2.根据权利要求1所述的地震波速度模型的更新方法，其特征在于，根据所述发生坐标和所述截割时间，获取所述微震事件的P波到时和S波到时，具体包括：
对所述微震事件进行数据预处理，所述预处理包括去噪处理、滤波处理和增益校正处理；
对进行预处理后的微震事件进行初至拾取，以获取所述微震事件的P波到时和S波到时。
3.根据权利要求1所述的地震波速度模型的更新方法，其特征在于，根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，更新预先构建的地震波速度模型，具体包括：
将所述发生坐标输入至所述速度模型中，得到地震波的预测旅行时；
根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，得到地震波的真实旅行时；
根据所述预测旅行时和所述真实旅行时，对所述速度模型进行更新。
4.根据权利要求3所述的地震波速度模型的更新方法，其特征在于，根据所述预测旅行时和所述真实旅行时，对所述速度模型进行更新，具体包括：
根据所述预测旅行时和所述真实旅行时，通过Levenberg-Marquardt算法对所述速度模型进行更新。
5.根据权利要求1所述的地震波速度模型的更新方法，其特征在于，获取检波器模块采集的微震事件，具体包括：
按照预设的时间间隔获取检波器模块采集的微震事件。
6.根据权利要求1所述的地震波速度模型的更新方法，其特征在于，在通过更新后的地震波速度模型，预测震源到所述检波器模块的旅行时后，所述方法还包括：
根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，进行叠前偏移成像，以获得所述微震事件的偏移成像结果图。
7.一种地震波速度模型的更新装置，其特征在于，包括：
获取单元，用于获取检波器模块采集的微震事件，所述微震事件包括微震事件的发生坐标和微震事件的截割时间；
计算单元，用于根据所述微震事件，获取所述微震事件的P波到时和S波到时；
更新单元，用于根据所述发生坐标、所述截割时间、所述P波到时和所述S波到时，更新预先构建的地震波速度模型；
预测单元，用于通过更新后的地震波速度模型，预测地震波的旅行时。
8.一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1至6任一项所述地震波速度模型的更新方法。
9.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至6任一项所述地震波速度模型的更新方法。
10.一种计算机程序产品，包括计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至6任一项所述地震波速度模型的更新方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310730120.6A CN116736383A (zh) | 2023-06-19 | 2023-06-19 | 地震波速度模型的更新方法、装置、电子设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310730120.6A CN116736383A (zh) | 2023-06-19 | 2023-06-19 | 地震波速度模型的更新方法、装置、电子设备及存储介质 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116736383A true CN116736383A (zh) | 2023-09-12 |
Family
ID=87907801
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310730120.6A Pending CN116736383A (zh) | 2023-06-19 | 2023-06-19 | 地震波速度模型的更新方法、装置、电子设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116736383A (zh) |
-
2023
- 2023-06-19 CN CN202310730120.6A patent/CN116736383A/zh active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR101548976B1 (ko) | 지진 표면파들의 파형들을 사용하는 토양 특성들의 추정 | |
| RU2457513C2 (ru) | Способы и системы для обработки микросейсмических данных | |
| US10073184B2 (en) | Sensor system of buried seismic array | |
| RU2518577C2 (ru) | Непрерывный адаптивный анализ поверхностных волн в случае трехмерных сейсмических данных | |
| KR20200014387A (ko) | 지하 구조물의 검출 | |
| US20130215716A1 (en) | Integrated Passive and Active Seismic Surveying Using Multiple Arrays | |
| EP2972502B1 (en) | System for seismic surveying of a subsurface volume | |
| US10345463B2 (en) | Methods and systems for using known source events in seismic data processing | |
| CN101086535A (zh) | 表征油藏随时间演变的方法和程序 | |
| WO2017035104A1 (en) | Velocity model seismic static correction | |
| Bellino et al. | An automatic method for data processing of seismic data in tunneling | |
| EA032186B1 (ru) | Сейсмическая адаптивная фокусировка | |
| CN111123359B (zh) | 随钻测井与地层格架约束的井周地震成像探测方法及装置 | |
| US20160161619A1 (en) | Systems and methods for aligning a monitor seismic survey with a baseline seismic survey | |
| US20150268365A1 (en) | Method to characterize geological formations using secondary source seismic data | |
| US20180164454A1 (en) | Determining residual statics of survey receivers | |
| Liu et al. | Detection of karst voids at pile foundation by full-waveform inversion of single borehole sonic data | |
| EP3359983B1 (en) | Seismic polynomial filter | |
| CN110579799B (zh) | 一种等旅行时间间隔的地震采集观测方法及*** | |
| Liu et al. | An inverted heterogeneous velocity model for microseismic source location in deep buried tunnels | |
| US20180335531A1 (en) | Salt proximity imaging using reverse time migration of transmitted vertical seismic profile data | |
| CN116736383A (zh) | 地震波速度模型的更新方法、装置、电子设备及存储介质 | |
| CN110579798B (zh) | 一种等反射角度间隔的地震采集观测方法及*** | |
| US20120143511A1 (en) | Removing ground roll from geophysical data | |
| CN113267810B (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 |