CN115437009A - 钻孔数据处理方法、装置、电子设备和计算机可读介质 - Google Patents

钻孔数据处理方法、装置、电子设备和计算机可读介质 Download PDF

Info

Publication number
CN115437009A
CN115437009A CN202211373224.8A CN202211373224A CN115437009A CN 115437009 A CN115437009 A CN 115437009A CN 202211373224 A CN202211373224 A CN 202211373224A CN 115437009 A CN115437009 A CN 115437009A
Authority
CN
China
Prior art keywords
drilling
shear wave
wave velocity
model
sample 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
CN202211373224.8A
Other languages
English (en)
Other versions
CN115437009B (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.)
INSTITUTE OF GEOPHYSICS CHINA EARTHQUAKE ADMINISTRATION
Original Assignee
INSTITUTE OF GEOPHYSICS CHINA EARTHQUAKE ADMINISTRATION
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 INSTITUTE OF GEOPHYSICS CHINA EARTHQUAKE ADMINISTRATION filed Critical INSTITUTE OF GEOPHYSICS CHINA EARTHQUAKE ADMINISTRATION
Priority to CN202211373224.8A priority Critical patent/CN115437009B/zh
Publication of CN115437009A publication Critical patent/CN115437009A/zh
Application granted granted Critical
Publication of CN115437009B publication Critical patent/CN115437009B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/303Analysis for determining velocity profiles or travel times
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282Application of seismic models, synthetic seismograms
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29Geographical information databases
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/622Velocity, density or impedance
    • G01V2210/6222Velocity; travel time

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Geophysics (AREA)
  • Databases & Information Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本公开的实施例公开了钻孔数据处理方法、装置、电子设备和计算机可读介质。该方法的一具体实施方式包括:获取钻孔数据集合;对于钻孔数据集合中的每个钻孔数据,执行以下剪切波速确定步骤:获取钻孔数据对应的钻孔参数信息,其中,钻孔参数信息包括以下中的至少一项:坡度钻孔参数、地质时代钻孔参数和土类钻孔参数;根据钻孔参数信息,从剪切波速模型中选取出对应钻孔参数信息的剪切波速子模型作为目标剪切波速子模型;将钻孔数据输入至目标剪切波速子模型,得到目标剪切波速;根据所得到的各个目标剪切波速,生成剪切波速地图;控制相关联的显示设备显示剪切波速地图。该实施方式提高了剪切波速地图表征的30米剪切波速参数的准确率。

Description

钻孔数据处理方法、装置、电子设备和计算机可读介质
技术领域
本公开的实施例涉及计算机技术领域,具体涉及钻孔数据处理方法、装置、电子设备和计算机可读介质。
背景技术
地表以下30米深度范围内的等效剪切波速(以下简称为30米剪切波速参数)是地震工程学中的重要参数,在地震区划、地震保险、地震应急响应等多个方面具有重要作用。目前,在实际应用30米剪切波速参数前,通常会先生成30米剪切波速参数,再根据30米剪切波速参数生成相应的剪切波速地图,以供在各个领域中进一步应用。例如,可通过双深度参数外推模型生成30米剪切波速参数。
然而,当采用上述方式生成30米剪切波速参数时,经常会存在如下技术问题:
生成的30米剪切波速参数的精度较低,导致剪切波速地图表征30米剪切波速参数的准确率较低。
发明内容
本公开的内容部分用于以简要的形式介绍构思,这些构思将在后面的具体实施方式部分被详细描述。本公开的内容部分并不旨在标识要求保护的技术方案的关键特征或必要特征,也不旨在用于限制所要求的保护的技术方案的范围。
本公开的一些实施例提出了钻孔数据处理方法、装置、电子设备和计算机可读介质,来解决以上背景技术部分提到的技术问题中的一项或多项。
第一方面,本公开的一些实施例提供了一种钻孔数据处理方法,该方法包括:获取钻孔数据集合;对于上述钻孔数据集合中的每个钻孔数据,执行以下剪切波速确定步骤:获取上述钻孔数据对应的钻孔参数信息,其中,上述钻孔参数信息包括以下中的至少一项:坡度钻孔参数、地质时代钻孔参数和土类钻孔参数;根据上述钻孔参数信息,从剪切波速模型中选取出对应上述钻孔参数信息的剪切波速子模型作为目标剪切波速子模型;将上述钻孔数据输入至上述目标剪切波速子模型,得到目标剪切波速;根据所得到的各个目标剪切波速,生成剪切波速地图;控制相关联的显示设备显示上述剪切波速地图。
第二方面,本公开的一些实施例提供了一种钻孔数据处理装置,装置包括:获取单元,被配置成获取钻孔数据集合;确定单元,被配置成对于上述钻孔数据集合中的每个钻孔数据,执行以下剪切波速确定步骤:获取上述钻孔数据对应的钻孔参数信息,其中,上述钻孔参数信息包括以下中的至少一项:坡度钻孔参数、地质时代钻孔参数和土类钻孔参数;根据上述钻孔参数信息,从剪切波速模型中选取出对应上述钻孔参数信息的剪切波速子模型作为目标剪切波速子模型;将上述钻孔数据输入至上述目标剪切波速子模型,得到目标剪切波速;生成单元,被配置成根据所得到的各个目标剪切波速,生成剪切波速地图;控制单元,被配置成控制相关联的显示设备显示上述剪切波速地图。
第三方面,本公开的一些实施例提供了一种电子设备,包括:一个或多个处理器;存储装置,其上存储有一个或多个程序,当一个或多个程序被一个或多个处理器执行,使得一个或多个处理器实现上述第一方面任一实现方式所描述的方法。
第四方面,本公开的一些实施例提供了一种计算机可读介质,其上存储有计算机程序,其中,程序被处理器执行时实现上述第一方面任一实现方式所描述的方法。
本公开的上述各个实施例具有如下有益效果:通过本公开的一些实施例的钻孔数据处理方法,可以提高剪切波速地图表征的30米剪切波速参数的准确率。具体来说,造成剪切波速地图表征30米剪切波速参数的准确率较低的原因在于:生成的30米剪切波速参数精度较低,导致剪切波速地图表征该30米剪切波速参数的准确率较低。基于此,本公开的一些实施例的地图显示方法,首先,获取钻孔数据集合。然后,对于钻孔数据集合中的每个钻孔数据,执行以下剪切波速确定步骤:获取钻孔数据对应的钻孔参数信息。由此,可以得到包括钻孔的相关参数的钻孔参数信息。根据钻孔参数信息,从剪切波速模型中选取出对应钻孔参数信息的剪切波速子模型作为目标剪切波速子模型。由此,可以得到对应上述钻孔参数信息的目标剪切波速子模型。将钻孔数据输入至目标剪切波速子模型,得到目标剪切波速。由此,可以得到对应上述钻孔数据的目标剪切波速。其次,根据所得到的各个目标剪切波速,生成剪切波速地图。由此,可以得到表征各个钻孔的剪切波速的剪切波速地图。最后,控制相关联的显示设备显示剪切波速地图。由此,可以显示上述剪切波速地图。因为通过从剪切波速模型中选取出对应钻孔参数信息的剪切波速子模型作为目标剪切波速子模型,得到目标剪切波速,以生成剪切波速地图,提高了生成的30米剪切波速参数的精度,从而提高了剪切波速地图表征的30米剪切波速参数的准确率。
附图说明
结合附图并参考以下具体实施方式,本公开各实施例的上述和其他特征、优点及方面将变得更加明显。贯穿附图中,相同或相似的附图标记表示相同或相似的元素。应当理解附图是示意性的,元件和元素不一定按照比例绘制。
图1是根据本公开的钻孔数据处理方法的一些实施例的流程图;
图2是根据本公开的钻孔数据处理装置的一些实施例的结构示意图;
图3是适于用来实现本公开的一些实施例的电子设备的结构示意图。
具体实施方式
下面将参照附图更详细地描述本公开的实施例。虽然附图中显示了本公开的某些实施例,然而应当理解的是,本公开可以通过各种形式来实现,而且不应该被解释为限于这里阐述的实施例。相反,提供这些实施例是为了更加透彻和完整地理解本公开。应当理解的是,本公开的附图及实施例仅用于示例性作用,并非用于限制本公开的保护范围。
另外还需要说明的是,为了便于描述,附图中仅示出了与有关发明相关的部分。在不冲突的情况下,本公开中的实施例及实施例中的特征可以相互组合。
需要注意,本公开中提及的“第一”、“第二”等概念仅用于对不同的装置、模块或单元进行区分,并非用于限定这些装置、模块或单元所执行的功能的顺序或者相互依存关系。
需要注意,本公开中提及的“一个”、“多个”的修饰是示意性而非限制性的,本领域技术人员应当理解,除非在上下文另有明确指出,否则应该理解为“一个或多个”。
本公开实施方式中的多个装置之间所交互的消息或者信息的名称仅用于说明性的目的,而并不是用于对这些消息或信息的范围进行限制。
下面将参考附图并结合实施例来详细说明本公开。
图1示出了根据本公开的钻孔数据处理方法的一些实施例的流程100。该钻孔数据处理方法,包括以下步骤:
步骤101,获取钻孔数据集合。
在一些实施例中,钻孔数据处理方法的执行主体(例如计算设备)可以通过有线连接方式或者无线连接方式从存储有钻孔数据集合的终端获取钻孔数据集合。其中,上述钻孔数据集合中的钻孔数据可以为钻孔的钻孔标识、钻孔坐标以及钻孔深度。上述钻孔标识可以唯一标识钻孔。上述钻孔可以为利用钻探设备向地下钻的柱状圆孔。例如,上述钻孔坐标可以为国家大地坐标系下上述钻孔的横纵坐标。上述钻孔深度可以为上述钻孔低于地面的深度,该深度的单位可以为米。需要指出的是,上述无线连接方式可以包括但不限于3G/4G连接、WiFi连接、蓝牙连接、WiMAX连接、Zigbee连接、UWB(ultra wideband)连接、以及其他现在已知或将来开发的无线连接方式。
需要说明的是,上述计算设备可以是硬件,也可以是软件。当计算设备为硬件时,可以实现成多个服务器或终端设备组成的分布式集群,也可以实现成单个服务器或单个终端设备。当计算设备体现为软件时,可以安装在上述所列举的硬件设备中。其可以实现成例如用来提供分布式服务的多个软件或软件模块,也可以实现成单个软件或软件模块。在此不做具体限定。应该理解,根据实现需要,可以具有任意数目的计算设备。
步骤102,对于钻孔数据集合中的每个钻孔数据,执行以下剪切波速确定步骤:
步骤1021,获取钻孔数据对应的钻孔参数信息。
在一些实施例中,上述执行主体可以获取上述钻孔数据对应的钻孔参数信息。其中,上述钻孔参数信息包括以下中的至少一项:坡度钻孔参数、地质时代钻孔参数和土类钻孔参数。其中,上述钻孔参数信息可以为包括钻孔的相关参数的信息。上述坡度钻孔参数可以为钻孔的坡度。上述地质时代钻孔参数可以为钻孔对应的地质年代。作为示例,上述地质年代钻孔参数可以为以下中的一项:古近纪之前、古近纪、新近纪、第四纪。其中,上述第四纪可以包括早中更新世、晚更新世、全新世。上述土类钻孔参数可以为钻孔的土类。上述土类钻孔参数可以为以下中的一项:岩土类和土类。其中,上述土类可以包括:“碎石土类”和“砂土/粉/黏土类”。上述“砂土/粉/黏土类”表征钻孔对应的土类钻孔参数为砂土、粉、黏土中的一项。实践中,上述执行主体可以通过有线方式或无线方式从存储有钻孔参数信息的终端获取钻孔参数信息。由此,可以得到包括钻孔的相关参数的钻孔参数信息。
步骤1022,根据钻孔参数信息,从剪切波速模型中选取出对应钻孔参数信息的剪切波速子模型作为目标剪切波速子模型。
在一些实施例中,上述执行主体可以根据上述钻孔参数信息,从剪切波速模型中选取出对应上述钻孔参数信息的剪切波速子模型作为目标剪切波速子模型。其中,上述剪切波速模型可以为预先训练的、包括至少一个剪切波速子模型的模型。上述剪切波速模型包括的剪切波速子模型可以为用于输入钻孔参数信息,以输出钻孔对应的剪切波速的模型。实践中,上述执行主体可以通过上述钻孔参数信息对应的坡度钻孔参数、地质时代钻孔参数以及土类钻孔参数,从剪切波速模型中选取出对应的坡度范围覆盖上述坡度钻孔参数、且对应的地质时代信息包括上述地质时代钻孔参数、对应的土类信息包括上述土类钻孔参数的剪切波速子模型作为目标剪切波速子模型。由此,可以得到对应上述钻孔参数信息的目标剪切波速子模型。
可选地,上述剪切波速模型可以是通过以下步骤训练得到的:
第一步,获取钻孔样本数据集,以及上述钻孔样本数据集中每个钻孔样本数据对应的坡度样本参数、地质时代样本参数和土类样本参数,得到坡度样本参数集合、地质时代样本参数集合和土类样本参数集合。实践中,可以通过有线连接或无限连接的方式从存储有钻孔样本数据集、坡度样本参数集合、地质时代样本参数集合和土类样本参数集合的终端获取上述钻孔样本数据集、上述坡度样本参数集合、上述地质时代样本参数集合和上述土类样本参数集合。
第二步,根据上述坡度样本参数集合、上述地质时代样本参数集合和上述土类样本参数集合,对上述钻孔样本数据集中的各个钻孔样本数据进行分组,得到钻孔样本数据组集。实践中,根据上述坡度样本参数集合、上述地质时代样本参数集合和上述土类样本参数集合,可以采用各种方式对上述钻孔样本数据集中的各个钻孔样本数据进行分组,得到钻孔样本数据组集。
可选地,在上述第二步之前,首先,上述执行主体可以将上述钻孔样本数据集中包括的钻孔位置信息满足预设钻孔位置条件的钻孔样本数据从上述钻孔样本数据集中删除,以对上述钻孔样本数据集进行更新。其中,上述钻孔位置信息可以为表征钻孔位置的横纵坐标。上述预设钻孔位置条件可以为位置精度大于1千米的钻孔样本数据。然后,可以将更新后的钻孔样本数据集中包括的钻孔地域信息满足预设钻孔地域条件的钻孔样本数据从更新后的钻孔样本数据集中删除,以对更新后的钻孔样本数据集再次更新。其中,上述钻孔地域信息可以为表征钻孔所属地域的信息。上述钻孔地域信息可以为以下中的一项:陆地区域、填海区域。上述预设钻孔地域条件可以为钻孔地域信息为填海区域的钻孔样本数据。由此,可以提高钻孔样本数据集的精度。
可选地,首先,上述执行主体可以根据坡度样本梯度范围信息和上述坡度样本参数集合,对上述钻孔样本数据集中的各个钻孔样本数据进行分组,得到至少一个坡度钻孔样本数据组。其中,上述坡度样本梯度范围信息可以表征样本的坡度的梯度范围的信息。作为示例,上述坡度样本梯度范围信息可以包括小于0.001,0.001-0.005,0.005-0.01,0.01-0.05,0.05-0.15和大于0.15这六个坡度样本梯度范围。实践中,可以根据上述坡度样本梯度范围信息和每个钻孔样本数据对应的坡度样本参数,按照该坡度样本参数在上述坡度样本梯度范围信息包括的坡度样本梯度范围,对相应的钻孔样本数据进行分组,以得到至少一个坡度钻孔样本数据组。作为示例,上述坡度样本梯度范围信息可以包括六个坡度样本梯度范围信息,则上述至少一个坡度钻孔样本数据组中坡度钻孔样本数据组的数量可以为六。其次,可以根据地质时代类别信息和上述地质时代样本参数集合,对上述钻孔样本数据集中的各个钻孔样本数据进行分组,得到至少一个地质钻孔样本数据组。其中,上述地质时代类别信息可以表征样本的地质时代。作为示例,上述地质时代类别信息可以包括古近纪之前类别、古近纪类别、新近纪类别和第四纪类别这四个地质时代类别。实践中,可以根据地质时代类别信息和每个钻孔样本数据对应的地质时代样本参数,按照该地质时代样本参数在上述地质时代类别信息包括的各个地质时代类别,对相应的钻孔样本数据进行分组,以得到至少一个地质钻孔样本数据组。作为示例,上述地质时代类别信息可以包括四个地质时代类别,则上述至少一个地质钻孔样本数据组中地质钻孔样本数据组的数量可以为四。然后,可以根据土类别信息和上述土类样本参数集合,对上述钻孔样本数据集中的各个钻孔样本数据进行分组,得到至少一个土类钻孔样本数据组。其中,上述土类别信息可以表征样本的土地质地的类别。作为示例,上述地质时代类别信息可以包括岩土类和土类。上述土类可以包括:“碎石土类”和“砂土/粉/黏土类”。实践中,可以根据土类别信息和每个钻孔样本数据对应的土类样本参数,按照该土类样本参数在上述地质时代类别信息包括的各个土地质地的类别,对相应的钻孔样本数据进行分组,以得到至少一个土类钻孔样本数据组。作为示例,上述土类别信息可以包括三个土地质地的类别,则上述至少一个土类钻孔样本数据组中土类钻孔样本数据组的数量可以为三。最后,可以将所得到的各个坡度钻孔样本数据组、各个地质钻孔样本数据组和各个土类钻孔样本数据组确定为钻孔样本数据组,得到钻孔样本数据组集。
第三步,根据上述钻孔样本数据组集中的每个钻孔样本数据组,对初始子模型进行模型训练,以将模型训练后的初始子模型作为剪切波速子模型,得到剪切波速子模型集合。实践中,根据上述钻孔样本数据组集中的每个钻孔样本数据组,可以采用各种方式对初始子模型进行模型训练,以将模型训练后的初始子模型作为剪切波速子模型,得到剪切波速子模型集合。由此,可以实现对于剪切波速子模型集合的训练。
可选地,上述初始子模型可以是通过以下步骤确定的:
第一子步骤,确定备选子模型集合。其中,上述备选子模型集合中的备选子模型对应有模型参数信息。其中,上述模型参数信息可以为包括用于生成剪切波速的模型系数的信息。用于生成剪切波速的模型系数可以包括对应钻孔深度的第一模型系数和第二模型系数。每个备选子模型对应的剪切波速生成公式如下式所示:
Figure 492637DEST_PATH_IMAGE001
Figure 803533DEST_PATH_IMAGE002
其中,
Figure DEST_PATH_IMAGE003
表示第一模型系数。
Figure 261059DEST_PATH_IMAGE004
表示第二模型系数。
Figure DEST_PATH_IMAGE005
表示钻孔深度(单位为米)。
Figure 622901DEST_PATH_IMAGE006
表示钻孔深度为
Figure 710943DEST_PATH_IMAGE005
的剪切波速(单位为m/s)。
Figure DEST_PATH_IMAGE007
表示备选子模型输出的剪切波速。
Figure 774714DEST_PATH_IMAGE008
表示以10为底的对数运算。
实践中,可以将各个现有的30米剪切波速(VS30)外推模型确定为备选子模型,得到备选子模型集合。现有的30米剪切波速(VS30)外推模型可以为已发表的论文中的剪切波速外推模型。作为示例,现有的30米剪切波速(VS30)外推模型可以为2015年Wills等人基于2000年发表的“A site-conditions map for California based on geology and shear-wave velocity”论文中的模型进行优化后所得到的模型。
第二子步骤,确定子模型样本数据集合。其中,上述子模型样本数据集合中的子模型样本数据可以为用于训练备选子模型集合中的备选子模型集合的样本数据。上述子模型样本数据可以包括但不限于:钻孔标识、钻孔坐标、钻孔深度、剪切波速、坡度样本参数、地质时代样本参数和土类样本参数。
第三子步骤,根据上述子模型样本数据集合,确定上述备选子模型集合中每个备选子模型的修正自由度信息,得到修正自由度信息集合。其中,上述修正自由度信息可以为包括校正决定系数的信息。
第四子步骤,根据上述子模型样本数据集合,确定上述备选子模型集合中每个备选子模型的均方根误差信息,得到均方根误差信息集合。
可选地,上述第四子步骤,包括:
第一,上述执行主体可以初始化迭代次数。例如,可以初始化迭代次数为1。
第二,对上述备选子模型集合可以执行以下信息生成步骤:
第一小步,根据预设迭代增长值,可以对迭代次数进行更新。其中,上述预设迭代增长值可以为预先设定的迭代的次数增长的值。例如,上述预设迭代增长值可以为1。实践中,可以将迭代次数与预设迭代增长值的和确定为更新后的迭代次数。
第二小步,根据上述子模型样本数据集合和预设数据集比例,确定子模型样本训练数据集合和子模型样本测试数据集合。其中,上述预设数据集比例可以为预先设定的子模型样本训练数据集合和子模型样本测试数据集合的比例。例如,上述预设数据集比例可以为9:1。上述执行主体可以将上述子模型样本数据集合中90%的子模型样本数据确定为子模型样本训练数据,得到子模型样本训练数据集合,以及将将上述子模型样本数据集合中10%的子模型样本数据确定为子模型样本测试数据,得到子模型样本测试数据集合。
第三小步,根据上述子模型样本训练数据集合,对上述备选子模型集合中的各个备选子模型进行模型训练,得到备选训练子模型集合。实践中,可以将子模型样本训练数据集合中的子模型样本训练数据代入每个备选子模型对应的剪切波速生成公式,以修正公式中的第一模型系数和第二模型系数,并将修正完成的备选子模型作为备选训练子模型,以得到备选训练子模型集合。其中,上述备选训练子模型对应有坡度范围、地质时代信息和土类信息。上述备选训练子模型对应的坡度范围可以为由训练上述备选训练子模型的各个子模型样本训练数据包括的坡度样本参数中的最大值和最小值所组成的取值范围。上述地质时代信息可以为表征训练上述备选训练子模型的各个子模型样本训练数据包括的地质时代样本参数的信息。上述土类信息可以为表征训练上述备选训练子模型的各个子模型样本训练数据包括的土类样本参数的信息。由此,可以通过修正第一模型系数和第二模型系数,提高各个备选子模型的模型精度,从而提高了模型输出的剪切波速的精度。此外,因为对上述备选子模型集合中的各个备选子模型进行模型训练,避免了直接对备选子模型进行筛选导致的可使用的备选子模型的丢弃,扩大了后续得到的初始子模型和剪切波速模型的可选择范围,使得剪切波速模型覆盖的情况更加全面,进一步提高了后续剪切波速模型输出的精度。
第四小步,根据上述子模型样本测试数据集合,对于上述备选训练子模型集合中的每个备选训练子模型,生成备选均方根误差。其中,上述备选均方根误差可以为备选训练子模型的均方根误差(RMSE,Root Mean Square Error)。
第五小步,响应于更新后的迭代次数小于预设迭代次数,再次执行上述信息生成步骤。其中,上述预设迭代次数可以为预先设置的次数。例如,上述预设迭代次数可以为10。
第三,响应于更新后的迭代次数大于等于上述预设迭代次数,对于上述备选训练子模型集合中的每个备选训练子模型,根据上述备选训练子模型对应的各个备选均方根误差,生成均方根误差信息。实践中,对于上述备选训练子模型集合中的每个备选训练子模型,可以将上述备选训练子模型对应的各个备选均方根误差的均值确定为均方根误差信息。
第五子步骤,根据上述修正自由度信息集合和上述均方根误差信息集合,将上述备选子模型集合中满足预设精度条件的备选子模型确定为初始子模型。其中,上述预设精度条件可以为备选子模型对应的修正自由度信息所表征的修正自由度大于等于预设修正自由度阈值,且对应的均方根误差信息表征的均方根误差小于等于预设均方根误差阈值。上述预设修正自由度阈值可以为预先设置的修正自由度的阈值。上述预设均方根误差阈值可以为预先设置的均方根误差的阈值。由此,可以根据精度进行筛选,得到精度更高的初始子模型。
第四步,对上述剪切波速子模型集合中的各个剪切波速子模型进行组合,得到剪切波速模型。实践中,可以将对应的坡度样本参数属于的坡度范围,对上述剪切波速子模型集合中的各个剪切波速子模型进行分组,并将每组中的各个剪切波速子模型确定为剪切波速模型所包括的各个剪切波速子模型,以得到剪切波速模型。其中,上述坡度范围可以为坡度样本参数所在的范围。例如,上述坡度范围可以包括:小于0.001、[0.0001-0.0005]、[0.0005-0.01]、[0.01-0.05]、大于0.05。
步骤1023,将钻孔数据输入至目标剪切波速子模型,得到目标剪切波速。
在一些实施例中,上述执行主体可以将上述钻孔数据输入至上述目标剪切波速子模型,得到目标剪切波速。由此,可以得到对应上述钻孔数据的目标剪切波速。
步骤103,根据所得到的各个目标剪切波速,生成剪切波速地图。
在一些实施例中,上述执行主体可以根据所得到的各个目标剪切波速,生成剪切波速地图。实践中,可以根据所生成的各个目标剪切波速属于的区间范围对应的颜色,对目标地图区域进行渲染。其中,上述目标地图区域可以为包括上述钻孔数据集合中各个钻孔数据在地图上的位置的区域。例如,可以将对应的值位于[100,200]区间的目标剪切波速在目标地图区域中渲染为红色,将对应的值位于[200,300]区间的目标剪切波速在目标地图区域中渲染为橙色。由此,可以得到表征各个钻孔的剪切波速的剪切波速地图。
步骤104,控制相关联的显示设备显示剪切波速地图。
在一些实施例中,上述执行主体可以控制相关联的显示设备显示上述剪切波速地图。上述相关联的显示设备可以为通信连接的智能设备。例如,上述显示设备可以为智能电脑。由此,可以显示上述剪切波速地图。
可选地,上述钻孔数据集合对应有目标区域。上述目标区域可以为上述目标地图区域。
可选地,首先,响应于检测到对应上述目标区域的地震信号,根据上述地震信号和预设地震动预测模型,上述执行主体可以生成至少一个基岩地震动反应谱。其中,上述至少一个基岩地震动反应谱中的基岩地震动反应谱对应上述目标区域。上述至少一个基岩地震动反应谱中的基岩地震动反应谱可以对应于目标点。上述至少一个基岩地震动反应谱中的基岩地震动反应谱对应有放大系数和特征周期等指标。上述预设地震动预测方程可以为预先设定的地震动预测方程(GMPE,Ground-Motion Prediction Equations)。上述地震信号可以表征发生地震的地震波。上述地震信号可以包括但不限于:地震波形参数、地震震级参数和震源深度参数。实践中,上述执行主体可以将上述地震信号包括的地震波形参数、地震震级参数和震源深度参数等各项参数代入上述预设地震动预测方程,并将预先获取的对应目标点的震中距离或断层距离依次输入已代入参数的预设地震动预测方程,得到至少一个基岩地震动反应谱。然后,根据所得到的各个目标剪切波速,可以对上述至少一个基岩地震动反应谱进行调整,以确定场地地震动反应谱。实践中,第一步,对于每个目标剪切波速,可以根据NEHRP Recommended Seismic Provisions for New Buildings and OtherStructures(以下简称为:NEHRP建议抗震规范)包括的场地类别和对应场地类别的剪切波速区间范围,以及上述目标剪切波速的大小,确定目标剪切波速对应的目标场地类别。第二步,可以将NEHRP建议抗震规范中上述目标场地类别对应的放大系数和特征周期等指标确定为目标指标。第三步,可以将上述目标指标确定为上述目标剪切波速对应的基岩地震动反应谱所对应的指标,以对基岩地震动反应谱对应的放大系数和特征周期等指标进行更新。第四步,可以将更新处理后的各个基岩地震动反应谱对应的指标的均值确定为场地地震动反应谱对应的指标,得到场地地震动反应谱。其次,根据上述场地地震动反应谱和预设预警等级,可以确定是否执行地震报警操作。其中,上述预设预警等级可以为预先设置的等级。实践中,第一步,可以根据上述场地地震动反应谱,生成地震强度。这里,对于生成地震强度的具体方法,不做限定。例如,可以将上述场地地震动反应谱中的峰值确定为地震强度。第二步,响应于上述地震强度大于等于上述预设预警等级,确定执行地震报警操作。第三步,响应于上述地震强度小于上述预设预警等级,确定不执行上述地震报警操作。最后,响应于确定执行上述地震报警操作,可以控制相关联的终端设备执行上述地震报警操作。实践中,可以控制通信连接的显示屏幕显示地震预警信息。上述地震预警信息可以为提示用户发生地震且提示用户地震等级的信息。由此,可以进行地震预警。
本公开的上述各个实施例具有如下有益效果:通过本公开的一些实施例的钻孔数据处理方法,可以提高剪切波速地图表征的30米剪切波速参数的准确率。具体来说,造成剪切波速地图表征30米剪切波速参数的准确率较低的原因在于:生成的30米剪切波速参数精度较低,导致剪切波速地图表征该30米剪切波速参数的准确率较低。基于此,本公开的一些实施例的地图显示方法,首先,获取钻孔数据集合。然后,对于钻孔数据集合中的每个钻孔数据,执行以下剪切波速确定步骤:获取钻孔数据对应的钻孔参数信息。由此,可以得到包括钻孔的相关参数的钻孔参数信息。根据钻孔参数信息,从剪切波速模型中选取出对应钻孔参数信息的剪切波速子模型作为目标剪切波速子模型。由此,可以得到对应上述钻孔参数信息的目标剪切波速子模型。将钻孔数据输入至目标剪切波速子模型,得到目标剪切波速。由此,可以得到对应上述钻孔数据的目标剪切波速。其次,根据所得到的各个目标剪切波速,生成剪切波速地图。由此,可以得到表征各个钻孔的剪切波速的剪切波速地图。最后,控制相关联的显示设备显示剪切波速地图。由此,可以显示上述剪切波速地图。因为通过从剪切波速模型中选取出对应钻孔参数信息的剪切波速子模型作为目标剪切波速子模型,得到目标剪切波速,以生成剪切波速地图,提高了生成的30米剪切波速参数的精度,从而提高了剪切波速地图表征的30米剪切波速参数的准确率。
继续参考图2,作为对上述各图所示方法的实现,本公开提供了一种钻孔数据处理装置的一些实施例,这些装置实施例与图1所示的那些方法实施例相对应,该装置具体可以应用于各种电子设备中。
如图2所示,一些实施例的钻孔数据处理装置200包括:获取单元201、确定单元202、生成单元203和控制单元204。其中,获取单元201被配置成获取钻孔数据集合;确定单元202被配置成对于上述钻孔数据集合中的每个钻孔数据,执行以下剪切波速确定步骤:获取上述钻孔数据对应的钻孔参数信息,其中,上述钻孔参数信息包括以下中的至少一项:坡度钻孔参数、地质时代钻孔参数和土类钻孔参数;根据上述钻孔参数信息,从剪切波速模型中选取出对应上述钻孔参数信息的剪切波速子模型作为目标剪切波速子模型;将上述钻孔数据输入至上述目标剪切波速子模型,得到目标剪切波速;生成单元203被配置成根据所得到的各个目标剪切波速,生成剪切波速地图;控制单元204被配置成控制相关联的显示设备显示上述剪切波速地图。
可以理解的是,该装置200中记载的诸单元与参考图1描述的方法中的各个步骤相对应。由此,上文针对方法描述的操作、特征以及产生的有益效果同样适用于装置200及其中包含的单元,在此不再赘述。
下面参考图3,其示出了适于用来实现本公开的一些实施例的电子设备(例如计算设备)300的结构示意图。图3示出的电子设备仅仅是一个示例,不应对本公开的实施例的功能和使用范围带来任何限制。
如图3所示,电子设备300可以包括处理装置(例如中央处理器、图形处理器等)301,其可以根据存储在只读存储器(ROM)302中的程序或者从存储装置308加载到随机访问存储器(RAM)303中的程序而执行各种适当的动作和处理。在RAM 303中,还存储有电子设备300操作所需的各种程序和数据。处理装置301、ROM 302以及RAM 303通过总线304彼此相连。输入/输出(I/O)接口305也连接至总线304。
通常,以下装置可以连接至I/O接口305:包括例如触摸屏、触摸板、键盘、鼠标、摄像头、麦克风、加速度计、陀螺仪等的输入装置306;包括例如液晶显示器(LCD)、扬声器、振动器等的输出装置307;包括例如磁带、硬盘等的存储装置308;以及通信装置309。通信装置309可以允许电子设备300与其他设备进行无线或有线通信以交换数据。虽然图3示出了具有各种装置的电子设备300,但是应理解的是,并不要求实施或具备所有示出的装置。可以替代地实施或具备更多或更少的装置。图3中示出的每个方框可以代表一个装置,也可以根据需要代表多个装置。
特别地,根据本公开的一些实施例,上文参考流程图描述的过程可以被实现为计算机软件程序。例如,本公开的一些实施例包括一种计算机程序产品,其包括承载在计算机可读介质上的计算机程序,该计算机程序包含用于执行流程图所示的方法的程序代码。在这样的一些实施例中,该计算机程序可以通过通信装置309从网络上被下载和安装,或者从存储装置308被安装,或者从ROM 302被安装。在该计算机程序被处理装置301执行时,执行本公开的一些实施例的方法中限定的上述功能。
需要说明的是,本公开的一些实施例中记载的计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质或者是上述两者的任意组合。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。计算机可读存储介质的更具体的例子可以包括但不限于:具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机访问存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本公开的一些实施例中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。而在本公开的一些实施例中,计算机可读信号介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式,包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读信号介质还可以是计算机可读存储介质以外的任何计算机可读介质,该计算机可读信号介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。计算机可读介质上包含的程序代码可以用任何适当的介质传输,包括但不限于:电线、光缆、RF(射频)等等,或者上述的任意合适的组合。
在一些实施方式中,客户端、服务器可以利用诸如HTTP(HyperText TransferProtocol,超文本传输协议)之类的任何当前已知或未来研发的网络协议进行通信,并且可以与任意形式或介质的数字数据通信(例如,通信网络)互连。通信网络的示例包括局域网(“LAN”),广域网(“WAN”),网际网(例如,互联网)以及端对端网络(例如,ad hoc端对端网络),以及任何当前已知或未来研发的网络。
上述计算机可读介质可以是上述电子设备中所包含的;也可以是单独存在,而未装配入该电子设备中。上述计算机可读介质承载有一个或者多个程序,当上述一个或者多个程序被该电子设备执行时,使得该电子设备:获取钻孔数据集合;对于上述钻孔数据集合中的每个钻孔数据,执行以下剪切波速确定步骤:获取上述钻孔数据对应的钻孔参数信息,其中,上述钻孔参数信息包括以下中的至少一项:坡度钻孔参数、地质时代钻孔参数和土类钻孔参数;根据上述钻孔参数信息,从剪切波速模型中选取出对应上述钻孔参数信息的剪切波速子模型作为目标剪切波速子模型;将上述钻孔数据输入至上述目标剪切波速子模型,得到目标剪切波速;根据所得到的各个目标剪切波速,生成剪切波速地图;控制相关联的显示设备显示上述剪切波速地图。
可以以一种或多种程序设计语言或其组合来编写用于执行本公开的一些实施例的操作的计算机程序代码,上述程序设计语言包括面向对象的程序设计语言—诸如Java、Smalltalk、C++,还包括常规的过程式程序设计语言—诸如“C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)——连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
附图中的流程图和框图,图示了按照本公开各种实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分,该模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个接连地表示的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或操作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
描述于本公开的一些实施例中的单元可以通过软件的方式实现,也可以通过硬件的方式来实现。所描述的单元也可以设置在处理器中,例如,可以描述为:一种处理器包括获取单元、确定单元、生成单元和控制单元。其中,这些单元的名称在某种情况下并不构成对该单元本身的限定,例如,获取单元还可以被描述为“获取钻孔数据集合的单元”。
本文中以上描述的功能可以至少部分地由一个或多个硬件逻辑部件来执行。例如,非限制性地,可以使用的示范类型的硬件逻辑部件包括:现场可编程门阵列(FPGA)、专用集成电路(ASIC)、专用标准产品(ASSP)、片上系统(SOC)、复杂可编程逻辑设备(CPLD)等等。
以上描述仅为本公开的一些较佳实施例以及对所运用技术原理的说明。本领域技术人员应当理解,本公开的实施例中所涉及的发明范围,并不限于上述技术特征的特定组合而成的技术方案,同时也应涵盖在不脱离上述发明构思的情况下,由上述技术特征或其等同特征进行任意组合而形成的其它技术方案。例如上述特征与本公开的实施例中公开的(但不限于)具有类似功能的技术特征进行互相替换而形成的技术方案。

Claims (10)

1.一种钻孔数据处理方法,包括:
获取钻孔数据集合;
对于所述钻孔数据集合中的每个钻孔数据,执行以下剪切波速确定步骤:
获取所述钻孔数据对应的钻孔参数信息,其中,所述钻孔参数信息包括以下中的至少一项:坡度钻孔参数、地质时代钻孔参数和土类钻孔参数;
根据所述钻孔参数信息,从剪切波速模型中选取出对应所述钻孔参数信息的剪切波速子模型作为目标剪切波速子模型;
将所述钻孔数据输入至所述目标剪切波速子模型,得到目标剪切波速;
根据所得到的各个目标剪切波速,生成剪切波速地图;
控制相关联的显示设备显示所述剪切波速地图。
2.根据权利要求1所述的方法,其中,所述钻孔数据集合对应有目标区域;以及
所述方法还包括:
响应于检测到对应所述目标区域的地震信号,根据所述地震信号和预设地震动预测模型,生成至少一个基岩地震动反应谱,其中,所述至少一个基岩地震动反应谱中的基岩地震动反应谱对应所述目标区域;
根据所得到的各个目标剪切波速,对所述至少一个基岩地震动反应谱进行调整,以确定场地地震动反应谱;
根据所述场地地震动反应谱和预设预警等级,确定是否执行地震报警操作;
响应于确定执行所述地震报警操作,控制相关联的终端设备执行所述地震报警操作。
3.根据权利要求1所述的方法,其中,所述剪切波速模型是通过以下步骤训练得到的:
获取钻孔样本数据集,以及所述钻孔样本数据集中每个钻孔样本数据对应的坡度样本参数、地质时代样本参数和土类样本参数,得到坡度样本参数集合、地质时代样本参数集合和土类样本参数集合;
根据所述坡度样本参数集合、所述地质时代样本参数集合和所述土类样本参数集合,对所述钻孔样本数据集中的各个钻孔样本数据进行分组,得到钻孔样本数据组集;
根据所述钻孔样本数据组集中的每个钻孔样本数据组,对初始子模型进行模型训练,以将模型训练后的初始子模型作为剪切波速子模型,得到剪切波速子模型集合;
对所述剪切波速子模型集合中的各个剪切波速子模型进行组合,得到剪切波速模型。
4.根据权利要求3所述的方法,其中,所述初始子模型是通过以下步骤确定的:
确定备选子模型集合,其中,所述备选子模型集合中的备选子模型对应有模型参数信息;
确定子模型样本数据集合;
根据所述子模型样本数据集合,确定所述备选子模型集合中每个备选子模型的修正自由度信息,得到修正自由度信息集合;
根据所述子模型样本数据集合,确定所述备选子模型集合中每个备选子模型的均方根误差信息,得到均方根误差信息集合;
根据所述修正自由度信息集合和所述均方根误差信息集合,将所述备选子模型集合中满足预设精度条件的备选子模型确定为初始子模型。
5.根据权利要求4所述的方法,其中,所述确定所述备选子模型集合中每个备选子模型的均方根误差信息,包括:
初始化迭代次数;
对所述备选子模型集合执行以下信息生成步骤:
根据预设迭代增长值,对迭代次数进行更新;
根据所述子模型样本数据集合和预设数据集比例,确定子模型样本训练数据集合和子模型样本测试数据集合;
根据所述子模型样本训练数据集合,对所述备选子模型集合中的各个备选子模型进行模型训练,得到备选训练子模型集合;
根据所述子模型样本测试数据集合,对于所述备选训练子模型集合中的每个备选训练子模型,生成备选均方根误差;
响应于更新后的迭代次数小于预设迭代次数,再次执行所述信息生成步骤;
响应于更新后的迭代次数大于等于所述预设迭代次数,对于所述备选训练子模型集合中的每个备选训练子模型,根据所述备选训练子模型对应的各个备选均方根误差,生成均方根误差信息。
6.根据权利要求3所述的方法,其中,所述钻孔样本数据集中的钻孔样本数据包括钻孔位置信息和钻孔地域信息;以及
在所述根据所述坡度样本参数集合、所述地质时代样本参数集合和所述土类样本参数集合,对所述钻孔样本数据集中的各个钻孔样本数据进行分组,得到钻孔样本数据组集之前,训练所述剪切波速模型的步骤还包括:
将所述钻孔样本数据集中包括的钻孔位置信息满足预设钻孔位置条件的钻孔样本数据从所述钻孔样本数据集中删除,以对所述钻孔样本数据集进行更新;
将更新后的钻孔样本数据集中包括的钻孔地域信息满足预设钻孔地域条件的钻孔样本数据从更新后的钻孔样本数据集中删除,以对更新后的钻孔样本数据集再次更新。
7.根据权利要求3所述的方法,其中,所述根据所述坡度样本参数集合、所述地质时代样本参数集合和所述土类样本参数集合,对所述钻孔样本数据集中的各个钻孔样本数据进行分组,得到钻孔样本数据组集,包括:
根据坡度样本梯度范围信息和所述坡度样本参数集合,对所述钻孔样本数据集中的各个钻孔样本数据进行分组,得到至少一个坡度钻孔样本数据组;
根据地质时代类别信息和所述地质时代样本参数集合,对所述钻孔样本数据集中的各个钻孔样本数据进行分组,得到至少一个地质钻孔样本数据组;
根据土类别信息和所述土类样本参数集合,对所述钻孔样本数据集中的各个钻孔样本数据进行分组,得到至少一个土类钻孔样本数据组;
将所得到的各个坡度钻孔样本数据组、各个地质钻孔样本数据组和各个土类钻孔样本数据组确定为钻孔样本数据组,得到钻孔样本数据组集。
8.一种钻孔数据处理装置,包括:
获取单元,被配置成获取钻孔数据集合;
确定单元,被配置成对于所述钻孔数据集合中的每个钻孔数据,执行以下剪切波速确定步骤:获取所述钻孔数据对应的钻孔参数信息,其中,所述钻孔参数信息包括以下中的至少一项:坡度钻孔参数、地质时代钻孔参数和土类钻孔参数;根据所述钻孔参数信息,从剪切波速模型中选取出对应所述钻孔参数信息的剪切波速子模型作为目标剪切波速子模型;将所述钻孔数据输入至所述目标剪切波速子模型,得到目标剪切波速;
生成单元,被配置成根据所得到的各个目标剪切波速,生成剪切波速地图;
控制单元,被配置成控制相关联的显示设备显示所述剪切波速地图。
9.一种电子设备,包括:
一个或多个处理器;
存储装置,其上存储有一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如权利要求1-7中任一所述的方法。
10.一种计算机可读介质,其上存储有计算机程序,其中,所述程序被处理器执行时实现如权利要求1-7中任一所述的方法。
CN202211373224.8A 2022-11-04 2022-11-04 钻孔数据处理方法、装置、电子设备和计算机可读介质 Active CN115437009B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211373224.8A CN115437009B (zh) 2022-11-04 2022-11-04 钻孔数据处理方法、装置、电子设备和计算机可读介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211373224.8A CN115437009B (zh) 2022-11-04 2022-11-04 钻孔数据处理方法、装置、电子设备和计算机可读介质

Publications (2)

Publication Number Publication Date
CN115437009A true CN115437009A (zh) 2022-12-06
CN115437009B CN115437009B (zh) 2023-03-24

Family

ID=84252898

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211373224.8A Active CN115437009B (zh) 2022-11-04 2022-11-04 钻孔数据处理方法、装置、电子设备和计算机可读介质

Country Status (1)

Country Link
CN (1) CN115437009B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4575830A (en) * 1982-10-15 1986-03-11 Schlumberger Technology Corporation Indirect shearwave determination
CN106094022A (zh) * 2016-06-06 2016-11-09 中国地震局地球物理勘探中心 一种土层等效剪切波速和土层反射界面埋深的测量方法
CN111208559A (zh) * 2020-01-09 2020-05-29 中国地震灾害防御中心 一种地震安全性评价计算装置和计算方法
CN113267814A (zh) * 2021-07-05 2021-08-17 中国地震局地球物理研究所 场地剪切波速测算方法及装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4575830A (en) * 1982-10-15 1986-03-11 Schlumberger Technology Corporation Indirect shearwave determination
CN106094022A (zh) * 2016-06-06 2016-11-09 中国地震局地球物理勘探中心 一种土层等效剪切波速和土层反射界面埋深的测量方法
CN111208559A (zh) * 2020-01-09 2020-05-29 中国地震灾害防御中心 一种地震安全性评价计算装置和计算方法
CN113267814A (zh) * 2021-07-05 2021-08-17 中国地震局地球物理研究所 场地剪切波速测算方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
何仲太 等: "衡水市岩土体剪切波速与土层深度的关系", 《科学技术与工程》 *
周健: "《中国博士学位论文全文数据库(电子期刊)》", 15 February 2022 *

Also Published As

Publication number Publication date
CN115437009B (zh) 2023-03-24

Similar Documents

Publication Publication Date Title
Ming et al. Real‐time flood forecasting based on a high‐performance 2‐D hydrodynamic model and numerical weather predictions
Wellmann et al. Uncertainties have a meaning: Information entropy as a quality measure for 3-D geological models
Bourne et al. A Monte Carlo method for probabilistic hazard assessment of induced seismicity due to conventional natural gas production
Snow et al. A high‐resolution national‐scale hydrologic forecast system from a global ensemble land surface model
JP5502084B2 (ja) 地震損害予測及び回避システム並びにその方法
Bevilacqua et al. The effects of vent location, event scale, and time forecasts on pyroclastic density current hazard maps at Campi Flegrei caldera (Italy)
CN103339627B (zh) 评估表面数据
US20130218538A1 (en) Simulation model optimization
CN115186037B (zh) 地图显示方法、装置、电子设备和计算机可读介质
CN104854482A (zh) 根据取芯井对非取芯井中的岩相预测
US11846175B2 (en) Estimating reservoir production rates using machine learning models for wellbore operation control
Panhalkar et al. Assessment of spatial interpolation techniques for river bathymetry generation of Panchganga River basin using geoinformatic techniques
Mignan et al. Considering large earthquake clustering in seismic risk analysis
CN110827405A (zh) 一种数码遥感地质测绘方法及系统
CN115688046A (zh) 一种岩爆预测方法、装置及计算机设备
Nikakhtar et al. Application of ANN-PSO algorithm based on FDM numerical modelling for back analysis of EPB TBM tunneling parameters
CN115437009B (zh) 钻孔数据处理方法、装置、电子设备和计算机可读介质
CN114371504B (zh) 一种地震震中位置确定方法、装置、设备及可读存储介质
CN115203893A (zh) 基于贝叶斯定理的砂岩比例不确定性区间缩小方法及系统
CN116484576A (zh) 一种基于贝叶斯更新的隧道沉降预测方法、装置、设备及介质
RU2591239C1 (ru) Системы и способы для оценивания возможности коллекторной системы
US20210149077A1 (en) Optimized methodology for automatic history matching of a petroleum reservoir model with ensemble kalman filter (enkf)
Ren et al. Data assimilation for surface wave method by ensemble Kalman filter with random field modeling
CN117648873B (zh) 地面沉降预测方法、训练方法、装置、设备
CN116797755B (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