CN114545500B - 波速分布信息的确定方法以及装置 - Google Patents
波速分布信息的确定方法以及装置 Download PDFInfo
- Publication number
- CN114545500B CN114545500B CN202210110494.3A CN202210110494A CN114545500B CN 114545500 B CN114545500 B CN 114545500B CN 202210110494 A CN202210110494 A CN 202210110494A CN 114545500 B CN114545500 B CN 114545500B
- Authority
- CN
- China
- Prior art keywords
- waveform
- wave velocity
- distribution information
- velocity distribution
- field
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000003325 tomography Methods 0.000 claims description 7
- 239000000463 material Substances 0.000 description 11
- 238000004590 computer program Methods 0.000 description 10
- 238000004891 communication Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer 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/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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
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
技术领域
本申请涉及地质勘探技术领域，尤其涉及一种波速分布信息的确定方法以及装置。
背景技术
在地质勘探中，相关技术中通常是在待探测场地中的多个位置上设置震源设备，各个位置上的震源设备激发出波形相似且不可人为控制的应力波，并通过迭代反演的方式来确定该待探测场地的波速分布信息。然而，这种方式所需要的实地探测工作量比较大，并且需要在不同空间位置开辟规整空间，激发完成后需要转移震源设备以及重新落位，或者需要若干套同种震源设备，从而造成确定待探测场地的波速分布信息所需的时间较长、工程成本较高、需要足够的场地空间。因此，如何在受限空间快速高效确定待探测场地内部的波速分布信息是亟需解决的问题。
发明内容
本申请一方面实施例提出了一种波速分布信息的确定方法，所述方法包括：控制待探测场地中的可控波形震源设备发出多个不同的震源波形；获取所述待探测场地中的检波器对多个不同的震源波形进行检测所得到的多个波形数据；根据所述多个波形数据和所述多个不同的震源波形，确定所述待探测场地中内部的波速分布信息。
在本申请的一个实施例中，所述多个不同的震源波形为具有相同特征频率f，而波形不同的应力波形。
在本申请的一个实施例中，在所述控制待探测场地中的可控波形震源设备发出多个不同的震源波形之前，所述方法还包括：获取应力波形可分辨的最小障碍物尺度；获取所述待探测场地的平均波速；根据所述最小障碍物尺度和所述平均波速，确定所述特征频率f。
在本申请的一个实施例中，所述根据所述多个波形数据和所述多个不同的震源波形，确定所述待探测场地中内部的波速分布信息，包括：根据初始波速分布信息，对每个震源波形进行波场正演模拟，以得到每个震源波形的正演模拟波形；确定所述正演模拟波形和所述多个波形数据的数据残差；对所述数据残差进行全波形反演，以得到波速分布信息的更新梯度与更新步长；在所述数据残差大于或者等于预设阈值时，根据所述更新梯度与更新步长对所述初始波速分布信息进行更新，并跳转至根据初始波速分布信息，对每个震源波形进行波场正演模拟，以得到每个震源波形的正演模拟波形的步骤，直至所述数据残差小于所述预设阈值；将所述数据残差小于预设阈值时所对应的更新后的初始波速分布信息作为所述待探测场地中内部的波速分布信息。
在本申请的一个实施例中，在所述根据初始波速分布信息，对每个震源波形进行波场正演模拟，以得到每个震源波形的正演模拟波形之前，所述方法还包括：对每个波形数据进行初至走时层析反演，获取每个波形数据的初始波速场；对各个波形数据的初始波速场进行线性叠加，以得到所述初始波速分布信息。
本申请提出一种波速分布信息的确定方法，通过控制待探测场地中的可控波形震源设备发出多个不同的震源波形，并获取待探测场地中的检波器对多个不同的震源波形进行检测所得到的多个波形数据，以及根据多个波形数据和多个不同的震源波形，确定待探测场地中内部的波速分布信息。由此，实现了通过一个可控波形震源设备所发出的多个不同的震源波形，准确确定出了待测场地中内部的波速分布信息。
本申请另一方面实施例提出了一种波速分布信息的确定装置，所述装置包括：控制模块，用于控制待探测场地中的可控波形震源设备发出多个不同的震源波形；第一获取模块，用于获取所述待探测场地中的检波器对多个不同的震源波形进行检测所得到的多个波形数据；第一确定模块，用于根据所述多个波形数据和所述多个不同的震源波形，确定所述待探测场地中内部的波速分布信息。
在本申请的一个实施例中，所述多个不同的震源波形为具有相同特征频率f，而波形不同的应力波形。
在本申请的一个实施例中，所述装置还包括：第二获取模块，用于获取应力波形可分辨的最小障碍物尺度；第三获取模块，用于获取所述待探测场地的平均波速；第二确定模块，用于根据所述最小障碍物尺度和所述平均波速，确定所述特征频率f。
在本申请的一个实施例中，所述第二确定模块，具体用于：根据初始波速分布信息，对每个震源波形进行波场正演模拟，以得到每个震源波形的正演模拟波形；确定所述正演模拟波形和所述多个波形数据的数据残差；对所述数据残差进行全波形反演，以得到波速分布信息的更新梯度与更新步长；在所述数据残差大于或者等于预设阈值时，根据所述更新梯度与更新步长对所述初始波速分布信息进行更新，并跳转至根据初始波速分布信息，对每个震源波形进行波场正演模拟，以得到每个震源波形的正演模拟波形的步骤，直至所述数据残差小于所述预设阈值；将所述数据残差小于预设阈值时所对应的更新后的初始波速分布信息作为所述待探测场地中内部的波速分布信息。
在本申请的一个实施例中，所述装置还包括：走时层析模块，用于对每个波形数据进行初至走时层析反演，以得到每个波形数据的初始波速场；叠加模块，用于对各个波形数据的初始波速场进行线性叠加，以得到所述初始波速分布信息。
本申请提出一种波速分布信息的确定装置，通过控制待探测场地中的可控波形震源设备发出多个不同的震源波形，并获取待探测场地中的检波器对多个不同的震源波形进行检测所得到的多个波形数据，以及根据多个波形数据和多个不同的震源波形，确定待探测场地中内部的波速分布信息。由此，实现了通过一个可控波形震源设备所发出的多个不同的震源波形，准确确定出了待测场地中内部的波速分布信息。
本申请另一方面实施例提出了一种电子设备，包括：存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现波速分布信息的确定方法。
本申请另一方面实施例提出了一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现波速分布信息的确定方法。
本申请另一方面实施例提出了一种计算机程序产品，包括计算机程序，所述计算机程序被处理器执行时实现波速分布信息的确定方法。
上述可选方式所具有的其他效果将在下文中结合具体实施例加以说明。
附图说明
图1是本申请实施例所提供的一种波速分布信息的确定方法的流程示意图。
图2是本申请实施例所提供的另一种波速分布信息的确定方法的流程示意图。
图3是本申请实施例所提供的另一种波速分布信息的确定方法的流程示意图。
图4是待反演波速场的示例图。
图5是通过本实施例方法进行反演所得到的反演结果的示例图。
图6是本申请实施例的一个波速分布信息的确定装置的结构示意图。
图7是本申请实施例的另一个波速分布信息的确定装置的结构示意图。
图8是根据本申请一个实施例的电子设备的框图。
具体实施方式
下面详细描述本申请的实施例，实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本申请，而不能理解为对本申请的限制。
下面参考附图描述本申请实施例的波速分布信息的确定方法、装置和电子设备。
图1是本申请实施例所提供的一种波速分布信息的确定方法的流程示意图。
如图1所示，该波速分布信息的确定方法可以包括：
步骤101，控制待探测场地中的可控波形震源设备发出多个不同的震源波形。
其中，需要说明的是，本实施例提供的波速分布信息的确定方法的执行主体为波速分布信息的确定装置，该波速分布信息的确定装置可以由软件和/或硬件的方式实现，该实施例中的波速分布信息的确定装置可以配置电子设备中，本实施例中的电子设备可以包括终端设备、服务器等，该实施例对电子设备不作具体限定。
在一些实施例中，上述待探测场地可以空间受限的场地。由于待探测场地的空闲受限，可以在待探测场地中的指定位置上设置一个可控波形震源设备。
其中，可控波形震源设备为可对固体介质激发不同震源波形的设备。在一些实施例中，上述可控波形震源设备可以包括多个冲击弹头，其中，为了可以使得可控波形震源设备可以产生不同的震源波形，本实施例中的冲击弹头所对应的截面是不同的。
在一些实施例中，上述多个不同的震源波形可以为具有不同的特征频率，且波形不同的应力波形，为了使得可控波形震源设备可准确产生具有多个特征频率不同，且波形不同的应力波形，在控制待探测场地中的可控波形震源设备发出多个不同的震源波形之前，还可以获取所需要的多个特征频率f，并获取与多个特征频率f所对应的可控波形震源设备。
其中，该可控波形震源设备包括多个冲击弹头，其中，冲击弹头与上述特征频率f一一对应。
其中，可以理解的是，对于上述可控波形震源设备中的各个冲击弹头，该冲击弹头是基于对应特征频率f所对应的弹头长度L、材料密度ρ、材料单轴压缩模量E这三个参数所制成的。
在另一些实施例中，在上述多个不同的震源波形为具有相同特征频率f，而波形不同的应力波形的情况下，为了使得可控波形震源设备可准确产生具有相同特征频率f，而波形不同的应力波形，在控制待探测场地中的可控波形震源设备发出多个不同的震源波形之前，还可以获取应力波形可分辨的最小障碍物尺度，并获取待探测场地的平均波速，以及根据最小障碍物尺度和平均波速，确定特征频率f。
在一些实施例中，在实际应用中，用户可根据实际应用需求，来预先设置应力波形可分辨的最小障碍物尺度u，也可以将应力波形可分辨的最小障碍物尺度称为反演精度。
在一些实施例中，在获取根据最小障碍物尺度u和平均波速v之后，可通过如下的公式，计算出特征频率f，其中：
f＝10*v/u。
可以理解的是，在确定出特征频率f后，可根据特征频率f，获得弹头长度L、材料密度ρ、材料单轴压缩模量E，并获取基于弹头长度L、材料密度ρ、材料单轴压缩模量E这三个参数为基准所加工而成的可控波形震源设备。其中，该可控波形震源设备包括多个冲击弹头，各个冲击弹头的截面虽然不同，但是各个冲击弹头均是基于弹头长度L、材料密度ρ、材料单轴压缩模量E这三个参数为基准进行加工的。
基于上述描述，可以看出，本实施例中可控波形震源设备不仅可以对波形进行控制，还可以对波形的频率进行控制。
其中，可控波形震源设备的可控性，是通过对可控波形震源设备的冲击弹头材料分布函数、截面分布函数进行控制，以实现对震源波形的完全可控激发。
步骤102，获取待探测场地中的检波器对多个不同的震源波形进行检测所得到的多个波形数据。
在一些实施例中，为了可以准确对震源波形进行检测，可以在待探测场地中的不同位置上设置多个检波器。
在一些示例性的实施方式中，可在待探测场地中设置单排或者多排的检波器，其中，每排检波器包括多个检波器。其中，多个检波器之间间隔预设距离。
步骤103，根据多个波形数据和多个不同的震源波形，确定待探测场地中内部的波速分布信息。
本申请实施例的波速分布信息的确定方法，通过控制待探测场地中的可控波形震源设备发出多个不同的震源波形，并获取待探测场地中的检波器对多个不同的震源波形进行检测所得到的多个波形数据，以及根据多个波形数据和多个不同的震源波形，确定待探测场地中内部的波速分布信息。由此，实现了通过一个可控波形震源设备所发出的多个不同的震源波形，准确确定出了待测场地中内部的波速分布信息。
在本申请的一个实施例中，为了可以准确确定出待探测场地中内部的波速分布信息，上述根据多个波形数据和多个不同的震源波形，确定待探测场地中内部的波速分布信息的一种可能的实现方式，如图2所示，可以包括：
步骤201，根据初始波速分布信息，对每个震源波形进行波场正演模拟，以得到每个震源波形的正演模拟波形。
在一些实施例中，可采用有限差分法(Finite Difference Method-FDM)，利用初始波速分布信息，对每个震源波形进行波场正演模拟，以得到每个震源波形的正演模拟波形。
在一些示例性的实施例中，上述有限差分法可以为九点有限差分法。
在本申请的一个实施例中，上述初始波速分布信息可以是预先设定好的。
在本申请的另一些实施例中，上述初始波速分布信息可以是通过下述方式得到的：对每个波形数据进行初至走时层析反演，以得到每个波形数据的初始波速场；对各个波形数据的初始波速场进行线性叠加，以得到初始波速分布信息。
具体地，针对每个波形数据，可确定该波形数据的起始点，并利用初至走时层析技术对该波形数据的起始点进行处理，以得到该波形数据对应的初始波速场。对应地，在确定出每个波形数据的初始波速场后，可将每个波形数据的初始波速场进行线性叠加，以得到初始波速分布信息，以得到初始波速分布信息。
步骤202，确定正演模拟波形和多个波形数据的数据残差。
步骤203，对数据残差进行全波形反演，以得到波速分布信息的更新梯度与更新步长。
在一些实施例中，可采用伴随状态法，对数据残差进行全波形反演，以得到波速分布信息的更新梯度与更新步长。
步骤204，在数据残差大于或者等于预设阈值时，根据更新梯度与更新步长对初始波速分布信息进行更新，并跳转至步骤201，直至数据残差小于预设阈值。
步骤205，将数据残差小于预设阈值时所对应的更新后的初始波速分布信息作为待探测场地中内部的波速分布信息。
在一些实施例中，在得到每个震源波形的反演波形后，可计算每个震源波形的反演波形和多个波形数据之间的二范数残差，并确定该二范数残差是否小于预设阈值，当确定该二范数残差小于预设阈值的情况下，则停止反演迭代过程，并将数据残差小于预设阈值时所对应的更新后的初始波速分布信息作为待探测场地中内部的波速分布信息。
其中，上述预设阈值是在波速分布信息的确定装置中预先设置的值，在实际应用中可根据实际需求在波速分布信息的确定装置中预先设置该预设阈值的取值，该实施例对此不作具体限定。
为了可以清楚了解本申请，下面结合图3对该实施例的方法进行进一步描述。
图3是本申请实施例所提供的另一种波速分布信息的确定方法的流程示意图。
如图3所示，该方法可以包括：
步骤301，获取应力波形可分辨的最小障碍物尺度。
步骤302，获取待探测场地的平均波速。
步骤303，根据最小障碍物尺度和平均波速，确定特征频率f。
步骤304，根据特征频率f，确定设置待探测场地中的可控波形震源设备。
在一些实施例中，在确定出设置待探测场地中的可控波形震源设备后，可在真实的待探测场地中设置该可控波形震源设备。
步骤305，控制待探测场地中的可控波形震源设备发出多个不同的震源波形。
作为一种示例性的实施方式，在真实的待探测场地中设置好可控波形震源设备和检波器后，可控制待探测场地中的可控波形震源设备发出多个不同的震源波形。
步骤306，获取待探测场地中的检波器对多个不同的震源波形进行检测所得到的多个波形数据。
其中，可以理解的是，上述波速分布信息的确定装置可与可控波形震源设备和检波器进行通信。
步骤307，对每个波形数据进行初至走时层析反演，以得到每个波形数据的初始波速场。
步骤308，对各个波形数据的初始波速场进行线性叠加，以得到初始波速分布信息。
步骤309，根据初始波速分布信息，对每个震源波形进行波场正演模拟，以得到每个震源波形的正演模拟波形。
步骤310，确定正演模拟波形和多个波形数据的数据残差。
步骤311，对数据残差进行全波形反演，以得到波速分布信息的更新梯度与更新步长。
步骤312，在数据残差大于或者等于预设阈值时，根据更新梯度与更新步长对初始波速分布信息进行更新，并跳转至步骤309，直至数据残差小于预设阈值。
步骤313，将数据残差小于预设阈值时所对应的更新后的初始波速分布信息作为待探测场地中内部的波速分布信息。
在本申请的实施例中，通过一个可控波形震源设备发出多个不同的震源波形，并基于检波器所检测到的波形进行全波形反演，并基于反演出的反演波形和震源波形，确定出待探测场地内部的波速分布信息。
例如，待反演波速场如图4所示，也就是说，假设检波器所检测到的波形数据如图4所示，在通过可控波形震源设备发出多个不同的震源波形，并通过本实施例的方法对该待反演波速场进行反演时，基于可控波形震源设备所发出的多个不同的震源波形，并基于检波器所检测到的波形进行全波形反演所得到的反演结果的示例，如图5所示，通过将图4和图5进行比较，可以看出，本实施例不但能够给出够给出非常精细的波速场反演结果。
图6是本申请实施例的一个波速分布信息的确定装置的结构示意图。
如图6所示，该波速分布信息的确定装置600包括：
控制模块601，用于控制待探测场地中的可控波形震源设备发出多个不同的震源波形。
第一获取模块602，用于获取待探测场地中的检波器对多个不同的震源波形进行检测所得到的多个波形数据。
第一确定模块603，用于根据多个波形数据和多个不同的震源波形，确定待探测场地中内部的波速分布信息。
在本申请的一个实施例中，多个不同的震源波形为具有相同特征频率f，而波形不同的应力波形。
在本申请的一个实施例中，在图6所示的装置实施例的基础上，如图7所示，该装置还可以包括：
第二获取模块604，用于获取应力波形可分辨的最小障碍物尺度。
第三获取模块605，用于获取待探测场地的平均波速。
第二确定模块606，用于根据最小障碍物尺度和平均波速，确定特征频率f。
在本申请的一个实施例中，上述第二确定模块606，具体用于：根据初始波速分布信息，对每个震源波形进行波场正演模拟，以得到每个震源波形的正演模拟波形；确定正演模拟波形和多个波形数据的数据残差；对数据残差进行全波形反演，以得到波速分布信息的更新梯度与更新步长；在数据残差大于或者等于预设阈值时，根据更新梯度与更新步长对初始波速分布信息进行更新，并跳转至根据初始波速分布信息，对每个震源波形进行波场正演模拟，以得到每个震源波形的正演模拟波形的步骤，直至数据残差小于预设阈值；将数据残差小于预设阈值时所对应的更新后的初始波速分布信息作为待探测场地中内部的波速分布信息。
在本申请的一个实施例中，如图7所示，该装置还可以包括：
走时层析模块607，用于对每个波形数据进行初至走时层析反演，以得到每个波形数据的初始波速场；
叠加模块608，用于对各个波形数据的初始波速场进行线性叠加，以得到初始波速分布信息。
其中，需要说明的是，前述对波速分布信息的确定方法实施例的解释说明也适用于本实施例的波速分布信息的确定装置，该实施例对此不再赘述。
本申请提出一种波速分布信息的确定装置，通过控制待探测场地中的可控波形震源设备发出多个不同的震源波形，并获取待探测场地中的检波器对多个不同的震源波形进行检测所得到的多个波形数据，以及根据多个波形数据和多个不同的震源波形，确定待探测场地中内部的波速分布信息。由此，实现了通过一个可控波形震源设备所发出的多个不同的震源波形，准确确定出了待测场地中内部的波速分布信息。
根据本申请的实施例，本申请还提供了一种电子设备。
如图8所示，是根据本申请一个实施例的电子设备的框图。
如图8所示，该电子设备包括：
存储器801、处理器802及存储在存储器801上并可在处理器802上运行的计算机指令。
处理器802执行指令时实现上述实施例中提供的波速分布信息的确定方法。
进一步地，电子设备还包括：
通信接口803，用于存储器801和处理器802之间的通信。
存储器801，用于存放可在处理器802上运行的计算机指令。
存储器801可能包含高速RAM存储器，也可能还包括非易失性存储器(non-volatile memory)，例如至少一个磁盘存储器。
处理器802，用于执行程序时实现上述实施例的波速分布信息的确定方法。
如果存储器801、处理器802和通信接口803独立实现，则通信接口803、存储器801和处理器802可以通过总线相互连接并完成相互间的通信。总线可以是工业标准体系结构(Industry Standard Architecture，简称为ISA)总线、外部设备互连(PeripheralComponent，简称为PCI)总线或扩展工业标准体系结构(Extended Industry StandardArchitecture，简称为EISA)总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示，图8中仅用一条粗线表示，但并不表示仅有一根总线或一种类型的总线。
可选的，在具体实现上，如果存储器801、处理器802及通信接口803，集成在一块芯片上实现，则存储器801、处理器802及通信接口803可以通过内部接口完成相互间的通信。
处理器802可能是一个中央处理器(Central Processing Unit，简称为CPU)，或者是特定集成电路(Application Specific Integrated Circuit，简称为ASIC)，或者是被配置成实施本申请实施例的一个或多个集成电路。
本申请还提出一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述实施例的波速分布信息的确定方法。
本申请还提出一种计算机程序产品，包括计算机程序，该计算机程序被处理器执行时实现上述实施例的波速分布信息的确定方法。
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本申请的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本申请的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本申请的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本申请的限制，本领域的普通技术人员在本申请的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (8)
1.一种波速分布信息的确定方法，其特征在于，所述方法包括：
控制待探测场地中的可控波形震源设备发出多个不同的震源波形；
获取所述待探测场地中的检波器对多个不同的震源波形进行检测所得到的多个波形数据；
根据所述多个波形数据和所述多个不同的震源波形，确定所述待探测场地中内部的波速分布信息；
所述根据所述多个波形数据和所述多个不同的震源波形，确定所述待探测场地中内部的波速分布信息，包括：
根据初始波速分布信息，对每个震源波形进行波场正演模拟，以得到每个震源波形的正演模拟波形；
确定所述正演模拟波形和所述多个波形数据的数据残差；
对所述数据残差进行全波形反演，以得到波速分布信息的更新梯度与更新步长；
在所述数据残差大于或者等于预设阈值时，根据所述更新梯度与更新步长对所述初始波速分布信息进行更新，并跳转至根据初始波速分布信息，对每个震源波形进行波场正演模拟，以得到每个震源波形的正演模拟波形的步骤，直至所述数据残差小于所述预设阈值；
将所述数据残差小于预设阈值时所对应的更新后的初始波速分布信息作为所述待探测场地中内部的波速分布信息。
2.如权利要求1所述的方法，其特征在于，所述多个不同的震源波形为具有相同特征频率f，而波形不同的应力波形。
3.如权利要求2所述的方法，其特征在于，在所述控制待探测场地中的可控波形震源设备发出多个不同的震源波形之前，所述方法还包括：
获取应力波形可分辨的最小障碍物尺度；
获取所述待探测场地的平均波速；
根据所述最小障碍物尺度和所述平均波速，确定所述特征频率f。
4.如权利要求1所述的方法，其特征在于，在所述根据初始波速分布信息，对每个震源波形进行波场正演模拟，以得到每个震源波形的正演模拟波形之前，所述方法还包括：
对每个波形数据进行初至走时层析反演，以得到每个波形数据的初始波速场；
对各个波形数据的初始波速场进行线性叠加，以得到所述初始波速分布信息。
5.一种波速分布信息的确定装置，其特征在于，所述装置包括：
控制模块，用于控制待探测场地中的可控波形震源设备发出多个不同的震源波形；
第一获取模块，用于获取所述待探测场地中的检波器对多个不同的震源波形进行检测所得到的多个波形数据；
第一确定模块，用于根据所述多个波形数据和所述多个不同的震源波形，确定所述待探测场地中内部的波速分布信息；
所述第一确定模块，具体用于：
根据初始波速分布信息，对每个震源波形进行波场正演模拟，以得到每个震源波形的正演模拟波形；
确定所述正演模拟波形和所述多个波形数据的数据残差；
对所述数据残差进行全波形反演，以得到波速分布信息的更新梯度与更新步长；
在所述数据残差大于或者等于预设阈值时，根据所述更新梯度与更新步长对所述初始波速分布信息进行更新，并跳转至根据初始波速分布信息，对每个震源波形进行波场正演模拟，以得到每个震源波形的正演模拟波形的步骤，直至所述数据残差小于所述预设阈值；
将所述数据残差小于预设阈值时所对应的更新后的初始波速分布信息作为所述待探测场地中内部的波速分布信息。
6.如权利要求5所述的装置，其特征在于，所述多个不同的震源波形为具有相同特征频率f，而波形不同的应力波形。
7.如权利要求6所述的装置，其特征在于，所述装置还包括：
第二获取模块，用于获取应力波形可分辨的最小障碍物尺度；
第三获取模块，用于获取所述待探测场地的平均波速；
第二确定模块，用于根据所述最小障碍物尺度和所述平均波速，确定所述特征频率f。
8.如权利要求5所述的装置，其特征在于，所述装置还包括：
走时层析模块，用于对每个波形数据进行初至走时层析反演，以得到每个波形数据的初始波速场；
叠加模块，用于对各个波形数据的初始波速场进行线性叠加，以得到所述初始波速分布信息。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210110494.3A CN114545500B (zh) | 2022-01-29 | 2022-01-29 | 波速分布信息的确定方法以及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210110494.3A CN114545500B (zh) | 2022-01-29 | 2022-01-29 | 波速分布信息的确定方法以及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114545500A CN114545500A (zh) | 2022-05-27 |
| CN114545500B true CN114545500B (zh) | 2022-09-09 |
Family
ID=81673822
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210110494.3A Active CN114545500B (zh) | 2022-01-29 | 2022-01-29 | 波速分布信息的确定方法以及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114545500B (zh) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113156493A (zh) * | 2021-05-06 | 2021-07-23 | 中国矿业大学 | 一种使用归一化震源的时频域全波形反演方法及装置 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7725266B2 (en) * | 2006-05-31 | 2010-05-25 | Bp Corporation North America Inc. | System and method for 3D frequency domain waveform inversion based on 3D time-domain forward modeling |
| CN101770038B (zh) * | 2010-01-22 | 2012-08-22 | 中国科学院武汉岩土力学研究所 | 矿山微震源智能定位方法 |
| CN103207409B (zh) * | 2013-04-17 | 2016-01-06 | 中国海洋石油总公司 | 一种频率域全波形反演地震速度建模方法 |
| CN104656123A (zh) * | 2015-01-23 | 2015-05-27 | 四川大学 | 一种测定区域岩体等效波速的方法 |
| US10670751B2 (en) * | 2015-03-27 | 2020-06-02 | Cgg Services Sas | Full waveform inversion method for seismic data processing using preserved amplitude reverse time migration |
| WO2019157242A1 (en) * | 2018-02-08 | 2019-08-15 | Schlumberger Technology Corporation | Ultrasonic acoustic sensors for measuring formation velocities |
| WO2020009752A1 (en) * | 2018-07-02 | 2020-01-09 | Exxonmobil Upstream Research Company | Full wavefield inversion with an image-gather-flatness constraint |
| CN109270574B (zh) * | 2018-10-18 | 2019-09-03 | 国家海洋局第二海洋研究所 | 一种基于多种震源采集数据的联合反褶积方法 |
-
2022
- 2022-01-29 CN CN202210110494.3A patent/CN114545500B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113156493A (zh) * | 2021-05-06 | 2021-07-23 | 中国矿业大学 | 一种使用归一化震源的时频域全波形反演方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114545500A (zh) | 2022-05-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9311566B2 (en) | Method and system for direct strain imaging | |
| US4644487A (en) | Method and apparatus for verifying the design of digital electronic components | |
| JP2003521015A (ja) | 非線形モデルの発生及びそれを用いたシミュレーション試験用駆動信号の発生 | |
| CN110954950B (zh) | 地下横波速度反演方法、装置、计算设备及存储介质 | |
| Remedia et al. | An enhanced methodology for spacecraft correlation activity using virtual testing tools | |
| JP2020125911A (ja) | 家屋の耐震性評価方法 | |
| CN113487019A (zh) | 电路故障诊断方法、装置、计算机设备和存储介质 | |
| US11892577B2 (en) | Multi-scale photoacoustic detection method of geological structure around borehole and related devices | |
| CN114545500B (zh) | 波速分布信息的确定方法以及装置 | |
| CN112773396B (zh) | 基于全波形反演的医学成像方法、计算机设备及存储介质 | |
| RU2544761C1 (ru) | Устройство для моделирования каталога разведки разнотипных подвижных объектов | |
| US20170016861A1 (en) | System and method for providing simulated ultrasound porosity waveforms | |
| CN113608262A (zh) | 利用平动分量计算旋转分量的地震数据处理方法及装置 | |
| CN106324688A (zh) | 一种储层束缚水饱和度确定方法及装置 | |
| CN115421195A (zh) | 地震探测中速度场的生成方法、装置、设备及存储介质 | |
| CN110824559B (zh) | 共坐标点道集的生成方法及装置、计算机可读存储介质 | |
| CN115356767A (zh) | 基于区域分割提高微震定位与波速反演精度的方法及装置 | |
| JP2008151560A (ja) | 振動特性推定方法及び振動特性推定装置 | |
| US20210334435A1 (en) | Method and device for simulating a technical system | |
| CN112308206B (zh) | 基于深度学习的大地电磁-地震波初至时联合反演方法 | |
| CN108376123A (zh) | 一种卫星重力观测数据异常值提取方法、装置及电子设备 | |
| CN113704085A (zh) | 用于检查技术系统的方法和设备 | |
| Chakraborty et al. | Estimation of in-plane elastic parameters and stiffener geometry of stiffened plates | |
| US20220164496A1 (en) | Quantum simulation device and operating method of quantum simulation device | |
| CN116299374B (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 |