CN110426739A - 一种地质勘探检测方法、装置和存储介质 - Google Patents
一种地质勘探检测方法、装置和存储介质 Download PDFInfo
- Publication number
- CN110426739A CN110426739A CN201910710584.4A CN201910710584A CN110426739A CN 110426739 A CN110426739 A CN 110426739A CN 201910710584 A CN201910710584 A CN 201910710584A CN 110426739 A CN110426739 A CN 110426739A
- Authority
- CN
- China
- Prior art keywords
- wave
- velocity
- surface layer
- rectifier
- component
- 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
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. analysis, for interpretation, for correction
- 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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
-
- 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/624—Reservoir parameters
Abstract
本发明实施例公开了一种地质勘探检测方法、装置和存储介质，其中，所述方法应用于四分量采集设备中；所述四分量采集设备包括用于采集高频分量数据的高频检波器和用于采集低频三分量数据的低频检波器；所述方法包括：通过所述高频检波器获得高频分量数据，通过所述低频检波器获得低频三分量数据；所述高频分量数据和所述低频三分量数据基于激发的人工震源获得；所述低频三分量数据包括两个水平分量数据和一个垂直分量数据；通过所述高频分量数据确定表层速度；基于所述表层速度和所述低频三分量数据获得速度模型，所述速度模型表征地层的深度和速度的对应关系；基于所述速度模型确定测试区域的地层速度剖面。
Description
技术领域
本发明涉及地震勘探技术领域，尤其涉及一种地质勘探检测方法、装置和存储介质。
背景技术
目前，地震学利用实际观测记录数据的反演方法仍多用于深部结构的研究，如单点地震噪音成像方法，该方法通常使用三分量低频检波器采集数据，再利用采集的数据得到测量的水平与垂直谱比曲线(Horizontal-to-Vertical Spectral Ratio，HVSR)，进而计算地层速度，但是在该计算过程中严重依赖地下表层横波速度。而使用HVSR单点噪音成像方法时，很难确定准确的表层横波速度，通常是给定经验值，这样会导致比较大的计算偏差，且计算结果不稳定。而针对该问题，目前尚无有效解决方案。
发明内容
有鉴于此，本发明实施例期望提供一种地质勘探检测方法、装置和存储介质。
本发明实施例的技术实施例是这样实现的：
本发明实施例提供一种地质勘探检测方法，所述方法应用于四分量采集设备中；所述四分量采集设备包括用于采集高频分量数据的高频检波器和用于采集低频三分量数据的低频检波器；所述方法包括：
通过所述高频检波器获得高频分量数据，通过所述低频检波器获得低频三分量数据；所述高频分量数据和所述低频三分量数据基于激发的人工震源获得；所述低频三分量数据包括两个水平分量数据和一个垂直分量数据；
通过所述高频分量数据确定表层速度；
基于所述表层速度和所述低频三分量数据获得速度模型，所述速度模型表征地层的深度和速度的对应关系；
基于所述速度模型确定测试区域的地层速度剖面。
在上述方案中，所述通过所述高频分量数据确定表层速度，包括：
通过所述高频分量数据确定所述人工震源发出的震动波到达所述高频检波器的时间；
基于所述时间确定所述震动波的波速；
基于所述波速确定表层速度。
在上述方案中，所述震动波至少包括以下之一：面波、纵波、横波；
通过所述高频分量数据确定所述人工震源发出的震动波到达所述高频检波器的时间，其中，所述时间至少包括以下之一：所述面波到达所述高频检波器的时间、所述纵波到达所述高频检波器的时间、所述横波到达所述高频检波器的时间；
所述基于所述时间确定所述震动波的波速，至少包括以下之一：基于所述面波到达所述高频检波器的时间确定所述面波的波速；
基于所述纵波到达所述高频检波器的时间确定所述纵波的波速；
基于所述横波到达所述高频检波器的时间确定所述横波的波速。
在上述方案中，所述基于所述波速确定表层速度，包括：
当所述波速为横波的波速时，将所述横波的波速确定为表层速度；
当所述波速为纵波的波速时，基于所述纵波的波速和第一比值确定表层速度；其中，所述第一比值为横波的波速与纵波的波速的比值；
当所述波速为面波的波速时，基于所述面波的波速和第二比值确定表层速度；其中，所述第二比值为横波的波速与面波的波速的比值。
在上述方案中，所述基于所述表层速度和所述低频三分量数据获得速度模型，包括：
基于所述表层速度建立初始模型；所述初始模型表征地层的深度和速度的对应关系；
对所述初始模型进行正演，获得模拟水平与垂直谱比HVSR曲线；
通过所述低频三分量数据获得测量HVSR曲线；
确定所述模拟HVSR曲线和所述测量HVSR曲线的误差因子，基于所述误差因子更新所述初始模型，确定速度模型。
在上述方案中，所述基于所述表层速度建立初始模型，包括：
根据所述表层速度和预设深度确定所述预设深度的地层速度；
基于所述表层速度和所述预设深度的地层速度建立初始模型。
在上述方案中，所述通过所述低频三分量数据获得测量HVSR曲线，包括：
通过所述两个水平分量数据确定水平分量的傅里叶谱；
通过所述垂直分量数据确定垂直分量的傅里叶谱；
将所述水平分量的傅里叶谱比上所述垂直分量的傅里叶谱获得测量HVSR曲线。
在上述方案中，所述基于所述速度模型确定测试区域的地层速度剖面，包括：
基于所述速度模型和测试区域的空间位置，生成所述测试区域的地层速度剖面。
本发明实施例提供一种地质勘探检测装置，所述装置包括：第一获得单元、第一确定单元、第二获得单元和第二确定单元，其中：
所述第一获得单元，用于通过高频检波器获得高频分量数据，通过低频检波器获得低频三分量数据；所述高频分量数据和所述低频三分量数据基于激发的人工震源获得；所述低频三分量数据包括两个水平分量数据和一个垂直分量数据；
所述第一确定单元，用于基于所述第一获得单元获得的所述高频分量数据确定表层速度；
所述第二获得单元，用于基于所述第一获得单元获得的所述低频三分量数据和所述第一确定单元确定的表层速度获得速度模型，所述速度模型表征地层的深度和速度的对应关系；
所述第二确定单元，用于基于所述第二获得单元获得的所述速度模型确定测试区域的地层速度剖面。
在上述方案中，所述第一确定单元，还用于通过所述高频分量数据确定所述人工震源发出的震动波到达所述高频检波器的时间；基于所述时间确定所述震动波的波速；基于所述波速确定表层速度。
在上述方案中，所述震动波至少包括以下之一：面波、纵波、横波；
通过所述高频分量数据确定所述人工震源发出的震动波到达所述高频检波器的时间，其中，所述时间至少包括以下之一：所述面波到达所述高频检波器的时间、所述纵波到达所述高频检波器的时间、所述横波到达所述高频检波器的时间；
所述基于所述时间确定所述震动波的波速，至少包括以下之一：基于所述面波到达所述高频检波器的时间确定所述面波的波速；
基于所述纵波到达所述高频检波器的时间确定所述纵波的波速；
基于所述横波到达所述高频检波器的时间确定所述横波的波速。
在上述方案中，所述第一确定单元，还用于当所述波速为横波的波速时，将所述横波的波速确定为表层速度；当所述波速为纵波的波速时，基于所述纵波的波速和第一比值确定表层速度；其中，所述第一比值为横波的波速与纵波的波速的比值；当所述波速为面波的波速时，基于所述面波的波速和第二比值确定表层速度；其中，所述第二比值为横波的波速与面波的波速的比值。
在上述方案中，所述第二获得单元包括：建立模块、获得模块和更新模块，其中
所述建立模块，用于基于所述表层速度建立初始模型；所述初始模型表征地层的深度和速度的对应关系；
所述获得模块，用于对所述建立模块中建立的初始模型进行正演，获得模拟水平与垂直谱比HVSR曲线；还用于通过所述低频三分量数据获得测量HVSR曲线；
所述更新模块，用于确定所述获得模块中的所述模拟HVSR曲线和所述测量HVSR曲线的误差因子，基于所述误差因子更新所述初始模型，确定速度模型。
在上述方案中，所述建立模块，还用于根据所述表层速度和预设深度确定所述预设深度的地层速度；基于所述表层速度和所述预设深度的地层速度建立初始模型。
在上述方案中，所述获得模块，还用于通过所述两个水平分量数据确定水平分量的傅里叶谱；通过所述垂直分量数据确定垂直分量的傅里叶谱；将所述水平分量的傅里叶谱比上所述垂直分量的傅里叶谱获得测量HVSR曲线。
在上述方案中，所述第二确定单元，还用于基于所述速度模型和测试区域的空间位置，生成所述测试区域的地层速度剖面。
本发明实施例提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述所述方法的任一步骤。
本发明实施例提供的一种地质勘探检测方法、装置和存储介质，其中，包括：通过所述高频检波器获得高频分量数据，通过所述低频检波器获得低频三分量数据；所述高频分量数据和所述低频三分量数据基于激发的人工震源获得；所述低频三分量数据包括两个水平分量数据和一个垂直分量数据；通过所述高频分量数据确定表层速度；基于所述表层速度和所述低频三分量数据获得速度模型，所述速度模型表征地层的深度和速度的对应关系；基于所述速度模型确定测试区域的地层速度剖面。采用本发明实施例的技术方案，通过高频分量数据确定表层速度，基于所述表层速度和所述低频三分量数据获得速度模型，所述速度模型表征地层的深度和速度的对应关系；相比于现有技术中直接给定表层速度(经验值)，大大提升了表层速度的准确率，进而减少速度模型的偏差，实现了地层速度的时效性和正确性。
附图说明
图1为本发明实施例一种地质勘探检测方法实现流程示意图；
图2为本发明实施例一种地质勘探检测方法中的高频分量数据示意图；
图3为本发明实施例一种地质勘探检测方法又一实现流程示意图；
图4为本发明实施例一种地质勘探检测方法中的一个测量HVSR曲线；
图5为本发明实施例一种地质勘探检测方法中的速度模型；
图6为本发明实施例一种地质勘探检测方法中得到的联合速度剖面；
图7为本发明实施例地质勘探检测装置的组成结构示意图；
图8为本发明实施例中地质勘探检测的一种硬件实体结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对发明的具体技术方案做进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
本实施例提出一种地质勘探检测方法，该方法应用于四分量采集设备中；所述四分量采集设备包括用于采集高频分量数据的高频检波器和用于采集低频三分量数据的低频检波器。
图1为本发明实施例一种地质勘探检测方法实现流程示意图，如图1所示，该方法包括：
步骤S101：通过所述高频检波器获得高频分量数据，通过所述低频检波器获得低频三分量数据；所述高频分量数据和所述低频三分量数据基于激发的人工震源获得；所述低频三分量数据包括两个水平分量数据和一个垂直分量数据。
需要说明的是，在实际应用中，先要根据勘察要求布置测点，在所述测点处设置四分量采集设备，所述四分量采集设备包括高频检波器和低频检波器。再在距离所述四分量采集设备预设距离的地方激发人工震源，该预设距离可以根据实际情况进行确定，在此不做限定。作为一种示例，该预设距离可以在10米左右，如果在表层速度较低的地区也可以适当增大距离；该人工震源的激发方式也可以根据实际情况进行确定，作为一种示例，该人工震源的激发方式可以为通过大锤敲击铁板激发震源或通过大锤敲击枕木激发震源。该人工震源激发的地方可以设置设备，该设备与四分量采集设备可以通过无线连接的方式或有线连接的方式连接，当人工震源激发的地方激发人工震源时，该设备便可以通过无线连接的方式或有线连接的方式连接触发四分量采集设备采集数据，并记录该数据。作为一种示例，该有线连接的方式可以为计时线的连接。由于该四分量采集设备包括用于采集高频分量数据的高频检波器和用于采集低频三分量数据的低频检波器，当人工震源激发的震动信号到达测点时，所述高频检波器便可以采集高频分量数据，并记录该高频分量数据；所述低频检波器便可以采集低频三分量数据，并记录该低频三分量数据。作为一种示例，该震动信号可以为震动波。
这里，通过所述高频检波器获得高频分量数据，其中，所述高频分量数据可以包括人工震源发出的震动波到达所述高频检波器的时间，该震动波可以为各种类型的波，在此不做限定，作为一种示例，该震动波可以为面波、纵波、横波等。
通过所述低频检波器获得低频三分量数据；所述低频三分量数据包括两个水平分量数据和一个垂直分量数据；其中，所述三分量数据可以为笛卡尔坐标的X、Y、Z三个分量，所述两个水平分量数据可以为X、Y两个分量数据，所述垂直分量数据可以为Z分量数据。
步骤S102：通过所述高频分量数据确定表层速度。
在本发明一种可选实施例中，所述通过所述高频分量数据确定表层速度可以包括：通过所述高频分量数据确定所述人工震源发出的震动波到达所述高频检波器的时间；基于所述时间确定所述震动波的波速；基于所述波速确定表层速度。
这里，所述震动波可以为各种类型的波，在此不做限定，作为一种示例，该震动波可以为面波、纵波、横波等。通过所述高频分量数据可以确定所述人工震源发出的各种类型的波到达所述高频检波器的时间；由于是在距离所述四分量采集设备预设距离的地方激发人工震源，虽然该预设距离是根据实际情况进行确定，但该预设距离是固定值，便可以基于该固定值和所述人工震源发出的各种类型的波到达所述高频检波器的时间确定该各种类型波的波速；再根据各种类型波的波速确定表层速度。
在本发明一种可选实施例中，所述震动波可以至少包括以下之一：面波、纵波、横波；
通过所述高频分量数据确定所述人工震源发出的震动波到达所述高频检波器的时间，其中，所述时间可以至少包括以下之一：所述面波到达所述高频检波器的时间、所述纵波到达所述高频检波器的时间、所述横波到达所述高频检波器的时间；
所述基于所述时间确定所述震动波的波速，可以至少包括以下之一：基于所述面波到达所述高频检波器的时间确定所述面波的波速；
基于所述纵波到达所述高频检波器的时间确定所述纵波的波速；
基于所述横波到达所述高频检波器的时间确定所述横波的波速。
这里，为了方便理解，进行示例说明，假设激发人工震源的地方距离所述四分量采集设备的距离为10米，图2为本发明实施例一种地质勘探检测方法中的高频分量数据示意图；如图2所示，在图2中记录了所述人工震源发出的震动波到达所述高频检波器的时间，该震动波可以为各种类型的波，在图2中仅标注出了面波和纵波到达的时间，其它的波并没有一一标注。该纵波到达时间为19ms，该面波到达时为58ms。再根据波速等于波传播的距离除以时间，在此，各种类型的波传播的距离为10米。通过波速公式的计算，便可以得到该纵波的波速为526m/s，该面波的波速为172m/s。
对应的，所述基于所述波速确定表层速度可以包括：当所述波速为横波的波速时，将所述横波的波速确定为表层速度；当所述波速为纵波的波速时，基于所述纵波的波速和第一比值确定表层速度；其中，所述第一比值为横波的波速与纵波的波速的比值；当所述波速为面波的波速时，基于所述面波的波速和第二比值确定表层速度；其中，所述第二比值为横波的波速与面波的波速的比值。
这里，所述表层速度可以为横波速度，所述第一比值为横波的波速与纵波的波速的比值，该第一比值在一般情况下，可以认为是个固定值。作为一种示例，该第一比值可以为1.732，当所述波速为纵波的波速时，且该纵波的波速为526m/s，基于所述纵波的波速和第一比值便可以确定表层速度为911m/s。所述第二比值为横波的波速与面波的波速的比值，该第二比值在一般情况下，也可以认为是个固定值。作为一种示例，该第二比值可以为1.1，当所述波速为面波的波速时，且该纵波的波速为172m/s，基于所述面波的波速和第二比值便可以确定表层速度为189m/s。
步骤S103：基于所述表层速度和所述低频三分量数据获得速度模型，所述速度模型表征地层的深度和速度的对应关系。
这里，基于所述表层速度和所述低频三分量数据获得速度模型可以为基于所述表层速度建立初始模型；所述初始模型表征地层的深度和速度的对应关系；对所述初始模型进行正演，获得模拟水平与垂直谱比HVSR曲线；再通过所述低频三分量数据获得测量HVSR曲线；将所述模拟HVSR曲线和所述测量HVSR曲线进行拟合，利用所述模拟HVSR曲线和所述测量HVSR曲线两者的差值对所述初始模型中的参数进行更新。不断迭代上述过程，直到所述模拟HVSR曲线和所述测量HVSR曲线的误差满足拟合误差或迭代参数满足预设值为止。该更新后的初始模型即为速度模型。
步骤S104：基于所述速度模型确定测试区域的地层速度剖面。
这里，所述测试区域可以为测点对应的区域，基于所述速度模型确定测试区域的地层速度剖面可以为基于所述速度模型和测试区域的空间位置，生成所述测试区域的地层速度剖面，具体，可以基于所述速度模型和对应的测点的空间位置，生成地层速度剖面。
本发明实施例提供的一种地质勘探检测方法，通过所述高频检波器获得高频分量数据，通过所述低频检波器获得低频三分量数据；所述高频分量数据和所述低频三分量数据基于激发的人工震源获得；所述低频三分量数据包括两个水平分量数据和一个垂直分量数据；通过所述高频分量数据确定表层速度；基于所述表层速度和所述低频三分量数据获得速度模型，所述速度模型表征地层的深度和速度的对应关系；基于所述速度模型确定测试区域的地层速度剖面。采用本发明实施例的技术方案，通过高频分量数据确定表层速度，基于所述表层速度和所述低频三分量数据获得速度模型，所述速度模型表征地层的深度和速度的对应关系；相比于现有技术中直接给定表层速度(经验值)，大大提升了表层速度的准确率，进而减少速度模型的偏差，实现了地层速度的时效性和正确性。
本实施例提出又一种地质勘探检测方法，图3为本发明实施例一种地质勘探检测方法又一实现流程示意图，如图3所示，该方法包括：
步骤S201，通过所述高频检波器获得高频分量数据，通过所述低频检波器获得低频三分量数据；所述高频分量数据和所述低频三分量数据基于激发的人工震源获得；所述低频三分量数据包括两个水平分量数据和一个垂直分量数据。
本实施例中的步骤S201可参照前述实施例中的步骤S101中的描述，这里不再赘述。
步骤S202，通过所述高频分量数据确定表层速度。
本实施例中的步骤S202可参照前述实施例中的步骤S102中的描述，这里不再赘述。
步骤S203，基于所述表层速度建立初始模型；所述初始模型表征地层的深度和速度的对应关系。
在本发明的一种可选实施例中，所述基于所述表层速度建立初始模型可以包括：根据所述表层速度和预设深度确定所述预设深度的地层速度；基于所述表层速度和所述预设深度的地层速度建立初始模型。
这里，预设深度可以根据实际情况进行确定，在此不做限定。作为一种示例，可以假设预设深度为2m，4m，6m等等。
根据所述表层速度和预设深度确定所述预设深度的地层速度可以为根据所述表层速度、预设深度和给定的速度增量确定所述预设深度的地层速度。所述给定的速度增量可以根据实际情况进行确定，在此不做限定。基于所述表层速度和所述预设深度的地层速度建立初始模型可以为基于所述表层速度和该表层速度的深度以及预设深度的地层速度和该地层速度对应的深度建立初始模型。为了方便理解，这里示例说明，假设所述表层速度为189m/s，深度为0m，由表层面向下每隔2m的速度增量为10m，即所述预设深度的地层速度可以为：深度为2m的地层速度为199m/s、深度为4m的地层速度为209m/s、深度为6m的地层速度为219m/s等等，进而根据不同深度和不同深度分别对应的地层速度建立初始模型。即所述初始模型表征地层的深度和速度的对应关系。
步骤S204，对所述初始模型进行正演，获得模拟水平与垂直谱比HVSR曲线。
这里，可以通过算法或软件对所述初始模型进行正演，便可获得模拟水平与垂直谱比HVSR曲线。
步骤S205，通过所述低频三分量数据获得测量HVSR曲线。
在本发明的一种可选实施例中，所述通过所述低频三分量数据获得测量HVSR曲线，包括：通过所述两个水平分量数据确定水平分量的傅里叶谱；通过所述垂直分量数据确定垂直分量的傅里叶谱；将所述水平分量的傅里叶谱比上所述垂直分量的傅里叶谱获得测量HVSR曲线。
为了方便理解，这里示例说明，图4为本发明实施例一种地质勘探检测方法中的一个测量HVSR曲线，如图4所示，纵坐标H/V表示HV谱比值，该HV谱比值为所述水平分量的傅里叶谱比上所述垂直分量的傅里叶谱，横坐标为频率。
步骤S206，确定所述模拟HVSR曲线和所述测量HVSR曲线的误差因子，基于所述误差因子更新所述初始模型，确定速度模型。
这里，确定所述模拟HVSR曲线和所述测量HVSR曲线的误差因子可以为将所述模拟HVSR曲线和所述测量HVSR曲线进行拟合，得到所述模拟HVSR曲线和所述测量HVSR曲线两者的差值。其中，所述误差因子即为所述模拟HVSR曲线和所述测量HVSR曲线两者的差值。基于所述误差因子更新所述初始模型可以为基于两者的差值对所述初始模型中的参数进行更新，该过程可以不断迭代，直到所述模拟HVSR曲线和所述测量HVSR曲线的误差满足拟合误差或迭代参数满足预设值为止。更新后的初始模型即为速度模型。图5为本发明实施例一种地质勘探检测方法中的速度模型，如图5所示，从图中可以看出，利用本发明实施例得到的速度模型能准确的反映试测区域地层的空间分布特征和速度变化规律。
步骤S207，基于所述速度模型和测试区域的空间位置，生成所述测试区域的地层速度剖面。
这里，所述测试区域可以为测点对应的区域，基于所述速度模型和测试区域的空间位置，生成所述测试区域的地层速度剖面可以为基于所述速度模型和对应的测点的空间位置，生成地层速度剖面。
在实际应用中，可以在测试区域中分别设置多个测点测量，针对每个测点都获得一个速度模型，将多个速度模型和分别对应的测点的空间位置进行联合，以获取整个测试区域的联合速度剖面。图6为本发明实施例一种地质勘探检测方法中得到的联合速度剖面。从图6中可以看出，利用本发明实施例得到的速度模型反映了测试区域地层的空间分布特征和速度变化规律。
本发明实施例的一种地质勘探检测方法，通过所述高频检波器获得高频分量数据，通过所述低频检波器获得低频三分量数据；所述高频分量数据和所述低频三分量数据基于激发的人工震源获得；所述低频三分量数据包括两个水平分量数据和一个垂直分量数据；通过所述高频分量数据确定表层速度；基于所述表层速度建立初始模型；所述初始模型表征地层的深度和速度的对应关系；对所述初始模型进行正演，获得模拟水平与垂直谱比HVSR曲线；通过所述低频三分量数据获得测量HVSR曲线；确定所述模拟HVSR曲线和所述测量HVSR曲线的误差因子，基于所述误差因子更新所述初始模型；确定速度模型基于所述速度模型和测试区域的空间位置，生成所述测试区域的地层速度剖面。采用本发明实施例的技术方案，通过高频分量数据确定表层速度，基于所述表层速度和所述低频三分量数据获得速度模型，所述速度模型表征地层的深度和速度的对应关系；相比于现有技术中直接给定表层速度(经验值)，大大提升了表层速度的准确率，进而减少速度模型的偏差，实现了地层速度的时效性和正确性。
本实施例提出一种地质勘探检测装置，图7为本发明实施例地质勘探检测装置的组成结构示意图，如图7所示，所述装置300包括：第一获得单元301、第一确定单元302、第二获得单元303和第二确定单元304，其中：
所述第一获得单元301，用于通过高频检波器获得高频分量数据，通过低频检波器获得低频三分量数据；所述高频分量数据和所述低频三分量数据基于激发的人工震源获得；所述低频三分量数据包括两个水平分量数据和一个垂直分量数据。
所述第一确定单元302，用于基于所述第一获得单元获得的所述高频分量数据确定表层速度。
所述第二获得单元303，用于基于所述第一获得单元获得的所述低频三分量数据和所述第一确定单元确定的表层速度获得速度模型，所述速度模型表征地层的深度和速度的对应关系。
所述第二确定单元304，用于基于所述第二获得单元获得的所述速度模型确定测试区域的地层速度剖面。
在其他的实施例中，所述第一确定单元302，还用于通过所述高频分量数据确定所述人工震源发出的震动波到达所述高频检波器的时间；基于所述时间确定所述震动波的波速；基于所述波速确定表层速度。
在其他的实施例中，所述震动波至少包括以下之一：面波、纵波、横波；
通过所述高频分量数据确定所述人工震源发出的震动波到达所述高频检波器的时间，其中，所述时间至少包括以下之一：所述面波到达所述高频检波器的时间、所述纵波到达所述高频检波器的时间、所述横波到达所述高频检波器的时间；
所述基于所述时间确定所述震动波的波速，至少包括以下之一：基于所述面波到达所述高频检波器的时间确定所述面波的波速；
基于所述纵波到达所述高频检波器的时间确定所述纵波的波速；
基于所述横波到达所述高频检波器的时间确定所述横波的波速。
在其他的实施例中，所述第一确定单元302，还用于当所述波速为横波的波速时，将所述横波的波速确定为表层速度；当所述波速为纵波的波速时，基于所述纵波的波速和第一比值确定表层速度；其中，所述第一比值为横波的波速与纵波的波速的比值；当所述波速为面波的波速时，基于所述面波的波速和第二比值确定表层速度；其中，所述第二比值为横波的波速与面波的波速的比值。
在其他的实施例中，所述第二获得单元303包括：建立模块、获得模块和更新模块，其中
所述建立模块，用于基于所述表层速度建立初始模型；所述初始模型表征地层的深度和速度的对应关系；
所述获得模块，用于对所述建立模块中建立的初始模型进行正演，获得模拟水平与垂直谱比HVSR曲线；还用于通过所述低频三分量数据获得测量HVSR曲线；
所述更新模块，用于确定所述获得模块中的所述模拟HVSR曲线和所述测量HVSR曲线的误差因子，基于所述误差因子更新所述初始模型，确定速度模型。
在其他的实施例中，所述建立模块，还用于根据所述表层速度和预设深度确定所述预设深度的地层速度；基于所述表层速度和所述预设深度的地层速度建立初始模型。
在其他的实施例中，所述获得模块，还用于通过所述两个水平分量数据确定水平分量的傅里叶谱；通过所述垂直分量数据确定垂直分量的傅里叶谱；将所述水平分量的傅里叶谱比上所述垂直分量的傅里叶谱获得测量HVSR曲线。
在其他的实施例中，所述第二确定单元304，还用于基于所述速度模型和测试区域的空间位置，生成所述测试区域的地层速度剖面。
以上装置实施例的描述，与上述方法实施例的描述是类似的，具有同方法实施例相似的有益效果。对于本发明装置实施例中未披露的技术细节，请参照本发明方法实施例的描述而理解。
需要说明的是，本发明实施例中，如果以软件功能模块的形式实现上述的一种地质勘探检测方法，并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实施例的技术实施例本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台地质勘探检测(可以是个人计算机、服务器、或者网络设备等)执行本发明各个实施例所述方法的全部或部分。而前述的存储介质包括：U盘、移动硬盘、只读存储器(Read Only Memory，ROM)、磁碟或者光盘等各种可以存储程序代码的介质。这样，本发明实施例不限制于任何特定的硬件和软件结合。
对应地，本发明实施例提供一种地质勘探检测，包括存储器和处理器，所述存储器存储有可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述实施例提供的一种地质勘探检测方法中的步骤。
对应地，本发明实施例提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例提供的一种地质勘探检测方法中的步骤。
这里需要指出的是：以上存储介质和设备实施例的描述，与上述方法实施例的描述是类似的，具有同方法实施例相似的有益效果。对于本发明存储介质和设备实施例中未披露的技术细节，请参照本发明方法实施例的描述而理解。
需要说明的是，图8为本发明实施例中地质勘探检测的一种硬件实体结构示意图，如图8所示，该地质勘探检测400的硬件实体包括：处理器401和存储器403，可选地，所述地质勘探检测400还可以包括通信接口402。
可以理解，存储器403可以是易失性存储器或非易失性存储器，也可包括易失性和非易失性存储器两者。其中，非易失性存储器可以是只读存储器(ROM，Read Only Memory)、可编程只读存储器(PROM，Programmable Read-Only Memory)、可擦除可编程只读存储器(EPROM，Erasable Programmable Read-Only Memory)、电可擦除可编程只读存储器(EEPROM，Electrically Erasable Programmable Read-Only Memory)、磁性随机存取存储器(FRAM，ferromagnetic random access memory)、快闪存储器(Flash Memory)、磁表面存储器、光盘、或只读光盘(CD-ROM，Compact Disc Read-Only Memory)；磁表面存储器可以是磁盘存储器或磁带存储器。易失性存储器可以是随机存取存储器(RAM，Random AccessMemory)，其用作外部高速缓存。通过示例性但不是限制性说明，许多形式的RAM可用，例如静态随机存取存储器(SRAM，Static Random Access Memory)、同步静态随机存取存储器(SSRAM，Synchronous Static Random Access Memory)、动态随机存取存储器(DRAM，Dynamic Random Access Memory)、同步动态随机存取存储器(SDRAM，SynchronousDynamic Random Access Memory)、双倍数据速率同步动态随机存取存储器(DDRSDRAM，Double Data Rate Synchronous Dynamic Random Access Memory)、增强型同步动态随机存取存储器(ESDRAM，Enhanced Synchronous Dynamic Random Access Memory)、同步连接动态随机存取存储器(SLDRAM，SyncLink Dynamic Random Access Memory)、直接内存总线随机存取存储器(DRRAM，Direct Rambus Random Access Memory)。本发明实施例描述的存储器403旨在包括但不限于这些和任意其它适合类型的存储器。
上述本发明实施例揭示的方法可以应用于处理器401中，或者由处理器401实现。处理器401可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤可以通过处理器401中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器401可以是通用处理器、数字信号处理器(DSP，Digital Signal Processor)，或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。处理器401可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者任何常规的处理器等。结合本发明实施例所公开的方法的步骤，可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于存储介质中，该存储介质位于存储器403，处理器401读取存储器403中的信息，结合其硬件完成前述方法的步骤。
在示例性实施例中，地质勘探检测可以被一个或多个应用专用集成电路(ASIC，Application Specific Integrated Circuit)、DSP、可编程逻辑器件(PLD，ProgrammableLogic Device)、复杂可编程逻辑器件(CPLD，Complex Programmable Logic Device)、现场可编程门阵列(FPGA，Field-Programmable Gate Array)、通用处理器、控制器、微控制器(MCU，Micro Controller Unit)、微处理器(Microprocessor)、或其他电子元件实现，用于执行前述方法。
在本发明所提供的几个实施例中，应该理解到，所揭露的方法和装置，可以通过其他的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，如：多个单元或组件可以结合，或可以集成到另一个观测量，或一些特征可以忽略，或不执行。另外，所显示或讨论的各组成部分相互之间的通信连接可以是通过一些接口，设备或单元的间接耦合或通信连接，可以是电性的、机械的或其他形式的。
上述作为分离部件说明的单元可以是、或也可以不是物理上分开的，作为单元显示的部件可以是、或也可以不是物理单元，即可以位于一个地方，也可以分布到多个网络单元上；可以根据实际的需要选择其中的部分或全部单元来实现本实施例的目的。
本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成，前述的程序可以存储于计算机可读取存储介质中，该程序在执行时，执行包括上述方法实施例的步骤；而前述的存储介质包括：移动存储设备、只读存储器(ROM，Read-Only Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
或者，本发明实施例上述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实施例的技术实施例本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台地质勘探检测(可以是个人计算机、服务器、或者网络设备等)执行本发明各个实施例所述方法的全部或部分。而前述的存储介质包括：移动存储设备、ROM、磁碟或者光盘等各种可以存储程序代码的介质。
本发明是实例中记载的基于确定卫星观测量质量方法、装置和计算机存储介质只以本发明所述实施例为例，但不仅限于此，只要涉及到该基于确定卫星观测量质量方法、装置和计算机存储介质均在本发明的保护范围。
应理解，说明书通篇中提到的“一个实施例”或“一实施例”意味着与实施例有关的特定特征、结构或特性包括在本发明的至少一个实施例中。因此，在整个说明书各处出现的“在一个实施例中”或“在一实施例中”未必一定指相同的实施例。此外，这些特定的特征、结构或特性可以任意适合的方式结合在一个或多个实施例中。应理解，在本发明的各种实施例中，上述各过程的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本发明实施例的实施过程构成任何限定。上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。
以上所述，仅为本发明的实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。
Claims (17)
1.一种地质勘探检测方法，其特征在于，所述方法应用于四分量采集设备中；所述四分量采集设备包括用于采集高频分量数据的高频检波器和用于采集低频三分量数据的低频检波器；所述方法包括：
通过所述高频检波器获得高频分量数据，通过所述低频检波器获得低频三分量数据；所述高频分量数据和所述低频三分量数据基于激发的人工震源获得；所述低频三分量数据包括两个水平分量数据和一个垂直分量数据；
通过所述高频分量数据确定表层速度；
基于所述表层速度和所述低频三分量数据获得速度模型，所述速度模型表征地层的深度和速度的对应关系；
基于所述速度模型确定测试区域的地层速度剖面。
2.根据权利要求1所述的方法，其特征在于，所述通过所述高频分量数据确定表层速度，包括：
通过所述高频分量数据确定所述人工震源发出的震动波到达所述高频检波器的时间；
基于所述时间确定所述震动波的波速；
基于所述波速确定表层速度。
3.根据权利要求2所述的方法，其特征在于，所述震动波至少包括以下之一：面波、纵波、横波；
通过所述高频分量数据确定所述人工震源发出的震动波到达所述高频检波器的时间，其中，所述时间至少包括以下之一：所述面波到达所述高频检波器的时间、所述纵波到达所述高频检波器的时间、所述横波到达所述高频检波器的时间；
所述基于所述时间确定所述震动波的波速，至少包括以下之一：基于所述面波到达所述高频检波器的时间确定所述面波的波速；
基于所述纵波到达所述高频检波器的时间确定所述纵波的波速；
基于所述横波到达所述高频检波器的时间确定所述横波的波速。
4.根据权利要求3所述的方法，其特征在于，所述基于所述波速确定表层速度，包括：
当所述波速为横波的波速时，将所述横波的波速确定为表层速度；
当所述波速为纵波的波速时，基于所述纵波的波速和第一比值确定表层速度；其中，所述第一比值为横波的波速与纵波的波速的比值；
当所述波速为面波的波速时，基于所述面波的波速和第二比值确定表层速度；其中，所述第二比值为横波的波速与面波的波速的比值。
5.根据权利要求1-4任意一项所述的方法，其特征在于，所述基于所述表层速度和所述低频三分量数据获得速度模型，包括：
基于所述表层速度建立初始模型；所述初始模型表征地层的深度和速度的对应关系；
对所述初始模型进行正演，获得模拟水平与垂直谱比HVSR曲线；
通过所述低频三分量数据获得测量HVSR曲线；
确定所述模拟HVSR曲线和所述测量HVSR曲线的误差因子，基于所述误差因子更新所述初始模型，确定速度模型。
6.根据权利要求5所述的方法，其特征在于，所述基于所述表层速度建立初始模型，包括：
根据所述表层速度和预设深度确定所述预设深度的地层速度；
基于所述表层速度和所述预设深度的地层速度建立初始模型。
7.根据权利要求5所述的方法，其特征在于，所述通过所述低频三分量数据获得测量HVSR曲线，包括：
通过所述两个水平分量数据确定水平分量的傅里叶谱；
通过所述垂直分量数据确定垂直分量的傅里叶谱；
将所述水平分量的傅里叶谱比上所述垂直分量的傅里叶谱获得测量HVSR曲线。
8.根据权利要求1所述的方法，其特征在于，所述基于所述速度模型确定测试区域的地层速度剖面，包括：
基于所述速度模型和测试区域的空间位置，生成所述测试区域的地层速度剖面。
9.一种地质勘探检测装置，其特征在于，所述装置包括：第一获得单元、第一确定单元、第二获得单元和第二确定单元，其中：
所述第一获得单元，用于通过高频检波器获得高频分量数据，通过低频检波器获得低频三分量数据；所述高频分量数据和所述低频三分量数据基于激发的人工震源获得；所述低频三分量数据包括两个水平分量数据和一个垂直分量数据；
所述第一确定单元，用于基于所述第一获得单元获得的所述高频分量数据确定表层速度；
所述第二获得单元，用于基于所述第一获得单元获得的所述低频三分量数据和所述第一确定单元确定的表层速度获得速度模型，所述速度模型表征地层的深度和速度的对应关系；
所述第二确定单元，用于基于所述第二获得单元获得的所述速度模型确定测试区域的地层速度剖面。
10.根据权利要求9所述的装置，其特征在于，所述第一确定单元，还用于通过所述高频分量数据确定所述人工震源发出的震动波到达所述高频检波器的时间；基于所述时间确定所述震动波的波速；基于所述波速确定表层速度。
11.根据权利要求10所述的装置，其特征在于，所述震动波至少包括以下之一：面波、纵波、横波；
通过所述高频分量数据确定所述人工震源发出的震动波到达所述高频检波器的时间，其中，所述时间至少包括以下之一：所述面波到达所述高频检波器的时间、所述纵波到达所述高频检波器的时间、所述横波到达所述高频检波器的时间；
所述基于所述时间确定所述震动波的波速，至少包括以下之一：基于所述面波到达所述高频检波器的时间确定所述面波的波速；
基于所述纵波到达所述高频检波器的时间确定所述纵波的波速；
基于所述横波到达所述高频检波器的时间确定所述横波的波速。
12.根据权利要求11所述的装置，其特征在于，所述第一确定单元，还用于当所述波速为横波的波速时，将所述横波的波速确定为表层速度；当所述波速为纵波的波速时，基于所述纵波的波速和第一比值确定表层速度；其中，所述第一比值为横波的波速与纵波的波速的比值；当所述波速为面波的波速时，基于所述面波的波速和第二比值确定表层速度；其中，所述第二比值为横波的波速与面波的波速的比值。
13.根据权利要求9至12任一项所述的装置，其特征在于，所述第二获得单元包括：建立模块、获得模块和更新模块，其中
所述建立模块，用于基于所述表层速度建立初始模型；所述初始模型表征地层的深度和速度的对应关系；
所述获得模块，用于对所述建立模块中建立的初始模型进行正演，获得模拟水平与垂直谱比HVSR曲线；还用于通过所述低频三分量数据获得测量HVSR曲线；
所述更新模块，用于确定所述获得模块中的所述模拟HVSR曲线和所述测量HVSR曲线的误差因子，基于所述误差因子更新所述初始模型，确定速度模型。
14.根据权利要求13所述的装置，其特征在于，所述建立模块，还用于根据所述表层速度和预设深度确定所述预设深度的地层速度；基于所述表层速度和所述预设深度的地层速度建立初始模型。
15.根据权利要求13所述的装置，其特征在于，所述获得模块，还用于通过所述两个水平分量数据确定水平分量的傅里叶谱；通过所述垂直分量数据确定垂直分量的傅里叶谱；将所述水平分量的傅里叶谱比上所述垂直分量的傅里叶谱获得测量HVSR曲线。
16.根据权利要求9所述的装置，其特征在于，所述第二确定单元，还用于基于所述速度模型和测试区域的空间位置，生成所述测试区域的地层速度剖面。
17.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至8任一项所述方法中的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910710584.4A CN110426739B (zh) | 2019-08-02 | 2019-08-02 | 一种地质勘探检测方法、装置和存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910710584.4A CN110426739B (zh) | 2019-08-02 | 2019-08-02 | 一种地质勘探检测方法、装置和存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110426739A true CN110426739A (zh) | 2019-11-08 |
| CN110426739B CN110426739B (zh) | 2021-07-16 |
Family
ID=68413751
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910710584.4A Active CN110426739B (zh) | 2019-08-02 | 2019-08-02 | 一种地质勘探检测方法、装置和存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110426739B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111425183A (zh) * | 2020-02-24 | 2020-07-17 | 中铁第四勘察设计院集团有限公司 | 一种地质勘探孔位定位方法及定位机器人 |
| CN112379418A (zh) * | 2020-10-16 | 2021-02-19 | 山东大学 | 一种地震直达波波速计算方法及系统 |
| CN112904348A (zh) * | 2021-01-22 | 2021-06-04 | 中铁第四勘察设计院集团有限公司 | 一种三维探测方法、装置、设备和存储介质 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103499835A (zh) * | 2013-10-13 | 2014-01-08 | 中国石油集团西北地质研究所 | 初至波波形反演近地表速度模型方法 |
| CN104536043A (zh) * | 2014-12-26 | 2015-04-22 | 中国石油天然气股份有限公司 | 一种深度域整体速度模型融合方法及装置 |
| CN104898162A (zh) * | 2014-03-06 | 2015-09-09 | 中国石油化工股份有限公司 | 地质勘探中的裂缝检测方法 |
| US20160291177A1 (en) * | 2015-03-30 | 2016-10-06 | Chevron U.S.A. Inc. | System and method for saltsurface updating via wavefield redatuming |
| CN108318918A (zh) * | 2017-12-29 | 2018-07-24 | 福建省建筑设计研究院有限公司 | 基于微动频散曲线和h/v曲线的地下不良地质体无损探测方法及应用 |
| WO2018158378A1 (en) * | 2017-03-01 | 2018-09-07 | Repsol Exploración, S.A. | Computer implemented method for improving a velocity model for seismic imaging |
| CN108646288A (zh) * | 2018-02-27 | 2018-10-12 | 中国石油天然气集团有限公司 | 近地表模型的建立方法和装置 |
| CN109061731A (zh) * | 2018-09-18 | 2018-12-21 | 中国地震局地壳应力研究所 | 面波频散与体波谱比联合反演浅层速度的全局优化方法 |
| CN109541690A (zh) * | 2018-11-30 | 2019-03-29 | 中铁第四勘察设计院集团有限公司 | 一种浅层介质结构面松散程度评价方法 |
-
2019
- 2019-08-02 CN CN201910710584.4A patent/CN110426739B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103499835A (zh) * | 2013-10-13 | 2014-01-08 | 中国石油集团西北地质研究所 | 初至波波形反演近地表速度模型方法 |
| CN104898162A (zh) * | 2014-03-06 | 2015-09-09 | 中国石油化工股份有限公司 | 地质勘探中的裂缝检测方法 |
| CN104536043A (zh) * | 2014-12-26 | 2015-04-22 | 中国石油天然气股份有限公司 | 一种深度域整体速度模型融合方法及装置 |
| US20160291177A1 (en) * | 2015-03-30 | 2016-10-06 | Chevron U.S.A. Inc. | System and method for saltsurface updating via wavefield redatuming |
| WO2018158378A1 (en) * | 2017-03-01 | 2018-09-07 | Repsol Exploración, S.A. | Computer implemented method for improving a velocity model for seismic imaging |
| CN108318918A (zh) * | 2017-12-29 | 2018-07-24 | 福建省建筑设计研究院有限公司 | 基于微动频散曲线和h/v曲线的地下不良地质体无损探测方法及应用 |
| CN108646288A (zh) * | 2018-02-27 | 2018-10-12 | 中国石油天然气集团有限公司 | 近地表模型的建立方法和装置 |
| CN109061731A (zh) * | 2018-09-18 | 2018-12-21 | 中国地震局地壳应力研究所 | 面波频散与体波谱比联合反演浅层速度的全局优化方法 |
| CN109541690A (zh) * | 2018-11-30 | 2019-03-29 | 中铁第四勘察设计院集团有限公司 | 一种浅层介质结构面松散程度评价方法 |
Non-Patent Citations (3)
| Title |
|---|
| 周翼 等: ""细分表层速度场约束层析静校正技术在西昆仑山黄土区的应用"", 《非常规油气》 * |
| 杨奎: ""微动面波H/V谱探测浅部地层横波速度结构的应用研究"", 《中国优秀硕士学位论文全文数据库•基础科学辑》 * |
| 王振宇 等: "近地表速度建模研究现状及发展趋势", 《地球科学进展》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111425183A (zh) * | 2020-02-24 | 2020-07-17 | 中铁第四勘察设计院集团有限公司 | 一种地质勘探孔位定位方法及定位机器人 |
| CN111425183B (zh) * | 2020-02-24 | 2023-12-08 | 中铁第四勘察设计院集团有限公司 | 一种地质勘探孔位定位方法及定位机器人 |
| CN112379418A (zh) * | 2020-10-16 | 2021-02-19 | 山东大学 | 一种地震直达波波速计算方法及系统 |
| CN112379418B (zh) * | 2020-10-16 | 2022-06-17 | 山东大学 | 一种地震直达波波速计算方法及系统 |
| CN112904348A (zh) * | 2021-01-22 | 2021-06-04 | 中铁第四勘察设计院集团有限公司 | 一种三维探测方法、装置、设备和存储介质 |
| CN112904348B (zh) * | 2021-01-22 | 2023-10-03 | 中铁第四勘察设计院集团有限公司 | 一种三维探测方法、装置、设备和存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110426739B (zh) | 2021-07-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110426739A (zh) | 一种地质勘探检测方法、装置和存储介质 | |
| CN108008459B (zh) | 一种获得剩余重力异常的方法及装置 | |
| CN105550448B (zh) | 基于钻井轨迹设计参数的预钻井三维井眼建模方法及装置 | |
| CN104793247B (zh) | 近地表速度建模方法及建模装置 | |
| CN105301641B (zh) | 一种方位各向异性速度反演的方法及装置 | |
| JP6001740B1 (ja) | 構造物の伝達関数の高精度評価・地震応答予測・劣化診断システム及びその方法 | |
| CN107329171A (zh) | 深度域储层地震反演方法及装置 | |
| CN107065013B (zh) | 一种地震尺度下的层速度确定方法及装置 | |
| CN107045143A (zh) | 一种预测裂缝发育的方法及装置 | |
| CN105223617A (zh) | 利用地震属性控制初始波阻抗的建模方法及系统 | |
| WO2023098441A1 (zh) | 基于地层记录沉降反演被动陆缘地壳结构的方法及装置 | |
| CN108873063A (zh) | 一种微地震矩张量反演的方法及装置 | |
| CN106021840A (zh) | 一种用于反演横观各向同性岩体地应力的方法 | |
| CN110426740A (zh) | 一种地震噪音成像勘探方法、装置和存储介质 | |
| CN106383367A (zh) | 绝对重力测量方法和装置 | |
| CN110426741A (zh) | 一种地震噪音成像勘探方法、装置和存储介质 | |
| CN108594300B (zh) | 地貌成像方法、装置及计算机存储介质 | |
| CN112180444B (zh) | 一种地层速度结构的探测方法、装置以及存储介质 | |
| Tran et al. | Inversion of first-arrival time using simulated annealing | |
| CN106125133B (zh) | 一种基于气云区约束下的精细速度建模方法 | |
| AU2013334868B2 (en) | Distributing petrofacies using analytical modeling | |
| CN103513279B (zh) | 一种基于地震波波动方程的照明分析计算方法及计算装置 | |
| CN106772599B (zh) | 一种计算地层横波速度的方法及装置 | |
| CN106950598B (zh) | 一种偏移速度场可靠性评价方法 | |
| CN114200538B (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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Liu Tiehua Inventor after: Bian Youyan Inventor after: Zhang Bang Inventor after: Li Kai Inventor after: Xiao Lifeng Inventor after: Liu Jian Inventor after: Zhao Xiaobo Inventor after: Liu Wei Inventor after: Duan Shenglong Inventor after: Liu Qing Inventor after: Guang Zhenxiong Inventor after: Sun Honglin Inventor after: Li Shiliang Inventor after: Zhang Zhanrong Inventor after: Liu Tie Inventor after: Hua Xirui Inventor after: Liao Jinxing Inventor after: Cui Dehai Inventor before: Liu Tiehua Inventor before: Liu Jian Inventor before: Zhao Xiaobo Inventor before: Liu Wei Inventor before: Duan Shenglong Inventor before: Liu Qing Inventor before: Liu Tie Inventor before: Hua Xirui Inventor before: Liao Jinxing Inventor before: Cui Dehai Inventor before: Bian Youyan Inventor before: Zhang Bang Inventor before: Li Kai Inventor before: Xiao Lifeng |
|
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20191108 Assignee: NATIONAL SCIENCE (CHONGQING) INSTRUMENT Co.,Ltd. Assignor: CHINA RAILWAY SIYUAN SURVEY AND DESIGN GROUP Co.,Ltd. Contract record no.: X2023980042888 Denomination of invention: A geological exploration detection method, device, and storage medium Granted publication date: 20210716 License type: Common License Record date: 20231007 Application publication date: 20191108 Assignee: CHONGQING SHENDI TECHNOLOGY CO.,LTD. Assignor: CHINA RAILWAY SIYUAN SURVEY AND DESIGN GROUP Co.,Ltd. Contract record no.: X2023980042879 Denomination of invention: A geological exploration detection method, device, and storage medium Granted publication date: 20210716 License type: Common License Record date: 20231007 |