CN113466930A - 微动勘探方法、装置及电子设备 - Google Patents
微动勘探方法、装置及电子设备 Download PDFInfo
- Publication number
- CN113466930A CN113466930A CN202110819615.7A CN202110819615A CN113466930A CN 113466930 A CN113466930 A CN 113466930A CN 202110819615 A CN202110819615 A CN 202110819615A CN 113466930 A CN113466930 A CN 113466930A
- Authority
- CN
- China
- Prior art keywords
- micro
- motion
- motion data
- measuring points
- exploration
- 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
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000001514 detection method Methods 0.000 claims abstract description 16
- 239000006185 dispersion Substances 0.000 claims description 19
- 239000011295 pitch Substances 0.000 claims description 7
- 238000010586 diagram Methods 0.000 description 9
- 238000000605 extraction Methods 0.000 description 8
- 238000004891 communication Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 239000000523 sample Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 241000282414 Homo sapiens Species 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/20—Arrangements of receiving elements, e.g. geophone pattern
-
- 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/301—Analysis for determining seismic cross-sections or geostructures
-
- 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/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/10—Aspects of acoustic signal generation or detection
- G01V2210/16—Survey configurations
-
- 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
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
技术领域
本发明涉及地质勘探技术领域,尤其是涉及一种微动勘探方法、装置及电子设备。
背景技术
目前,城市地面塌陷、堤坝溃坝等安全问题频发,威胁着人类生命与财产安全,迫切需要一种可以在城市等干扰因素较多的地方进行地质探测的方法。现有勘探技术中,地质雷达的探测深度较浅,不能满足深度要求;电法勘探在城市路面难以布置电极,无法实施。而微动勘探受外界干扰小、无需人工源激发,且测深较大,可以满足城市地质调查的需求。
以往利用微动勘探技术进行数据采集时,一次布设台阵仅能采集得到一个勘探点的微动数据;如果需要完成一条剖面的勘探或者进行大面积区域勘探,则需要多次布设台阵多次采集,这种勘探方式的工作量大且耗时长,勘探效率较低;另外,由于同一剖面的不同勘探点的采集时间有先后,而不同时间采集的微动数据可能是来源于不同震源,且可能受到不同干扰源影响,这就使得同一勘探剖面的不同勘探点采集的微动数据之间的背景场差异增加,以致降低利用剖面数据反演得到的地下地层的横向分辨率,也降低了勘探结果的准确性。
发明内容
有鉴于此,本发明的目的在于提供一种微动勘探方法、装置及电子设备,可以提高勘探效率,提高利用剖面数据反演得到的地下地层的横向分辨率,提升勘探结果的准确性。
第一方面,本发明实施例提供了一种微动勘探方法,包括:基于预先布设的微动观测装置采集目标测线的原始微动数据;其中,该微动观测装置按预设道间距布设检波器,且布设的采集道数远大于单个测点的基本道数,该基本道数为满足预设探测深度要求的道数;基于该基本道数,从上述原始微动数据中确定多个连续测点的微动数据,得到该多个连续测点中每个测点对应的微动数据;根据该多个连续测点中每个测点对应的微动数据,确定该多个连续测点所在区间的地质勘探成果。
在本发明较佳的实施例中,基于该基本道数,从该原始微动数据中确定多个连续测点的微动数据,得到该多个连续测点中每个测点对应的微动数据的步骤,包括:以该基本道数为滑动窗口,以检波器布设区间的第一端点作为该滑动窗口的起点,按预设步长移动该滑动窗口,抽取每个该滑动窗口对应的道集数据,直至该滑动窗口移动到该检波器布设区间的第二端点为止,得到多个道集数据;其中,该滑动窗口的中间点对应为当前测点,该滑动窗口对应的道集数据为该当前测点的微动数据。
在本发明较佳的实施例中,上述步长的取值为一个或多个道间距。
在本发明较佳的实施例中,上述采集道数为24道的正整数倍。
在本发明较佳的实施例中,上述基本道数根据下述公式确定:D=K×(N-1)×R,其中,D为该预设探测深度,N为该基本道数,R为该预设道间距,K为常系数,表示该预设探测深度与观测台阵之间的倍数。
在本发明较佳的实施例中,根据该多个连续测点中每个测点对应的微动数据,确定该多个连续测点所在区间的地质勘探成果的步骤,包括:根据该多个连续测点中每个测点对应的微动数据,提取每个该测点的频散曲线;根据提取到的该频散曲线,反演得到该多个连续测点所在区间的地质剖面图。
在本发明较佳的实施例中,该频散曲线为速度-深度域曲线。
在本发明较佳的实施例中,在根据提取到的该频散曲线,反演得到该多个连续测点所在区间的地质剖面图的步骤之后,该方法还包括:在该微动观测装置的预设显示终端显示该地质剖面图。
第二方面,本发明实施例提供了一种微动勘探装置,包括:原始微动数据采集模块,用于基于预先布设的微动观测装置采集目标测线的原始微动数据;其中,该微动观测装置按预设道间距布设检波器,且布设的采集道数远大于单个测点的基本道数,该基本道数为满足预设探测深度要求的道数;测点微动数据确定模块,用于基于该基本道数,从该原始微动数据中确定多个连续测点的微动数据,得到该多个连续测点中每个测点对应的微动数据;地质勘探成果确定模块,用于根据该多个连续测点中每个测点对应的微动数据,确定该多个连续测点所在区间的地质勘探成果。
第三方面,本发明实施例提供了一种电子设备,该电子设备包括处理器和存储器,该存储器存储有能够被该处理器执行的计算机可执行指令,该处理器执行该计算机可执行指令以实现上述微动勘探方法。
本发明实施例带来了以下有益效果:
本发明实施例提供的一种微动勘探方法、装置及电子设备,基于预先布设的微动观测装置采集目标测线的原始微动数据;其中,该微动观测装置按预设道间距布设检波器,且布设的采集道数远大于单个测点的基本道数,该基本道数为满足预设探测深度要求的道数;基于该基本道数,从上述原始微动数据中确定多个连续测点的微动数据,得到该多个连续测点中每个测点对应的微动数据;根据该多个连续测点中每个测点对应的微动数据,确定该多个连续测点所在区间的地质勘探成果。该方式中,通过一次性布设远大于单个测点基本道数的采集道数,以超常规道数的方式布设检波器,从而一次性采集到多个测点的微动数据,并通过抽道的方式提取出多个测点中每个测点对应的微动数据,从而提高了勘探效率;另一方面,由于提取出的多个测点的微动数据的采集时间一致,因而可以保证各个测点的微动数据来源于相同震源,且受到的干扰也较为一致,从而减小了不同测点采集的微动数据之间的背景场差异,进而提升反演结果对地下地层的横向分辨率,提升勘探结果的准确性。
本公开的其他特征和优点将在随后的说明书中阐述,或者,部分特征和优点可以从说明书推知或毫无疑义地确定,或者通过实施本公开的上述技术即可得知。
为使本公开的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种微动勘探方法的流程示意图;
图2为本发明实施例提供的一种微动观测装置的布置示意图;
图3为本发明实施例提供的一种通过抽道集的方式提取测点微动数据的示意图;
图4为本发明实施例提供的一种反演地质剖面的示意图;
图5为本发明实施例提供的一种微动勘探装置的结构示意图;
图6为本发明实施例提供的一种电子设备的结构示意图。
图标:21-主机;22-检波器;23-大线;24-连接线;25-拖拽带;26-机动车;51-原始微动数据采集模块;52-测点微动数据确定模块;53-地质勘探成果确定模块;61-处理器;62-存储器;63-总线;64-通信接口。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
地球表面的任何地方和时刻即使在无地震之时都处在一种微弱的振动状态下,振幅很小(位移约为10-5~10-3cm),地球表面的这种连续的微弱振动就称为微动。微动信号属于天然源信号,其来源于两个方面:一是来源于人类的日常活动,包括各种机械振动,道路交通等,这些活动产生的信号频率大于1HZ,属于高频信号源,通常,类微动信号被称为常时微动或地脉动;二是来源于各种自然现象,包括海浪对海岸的撞击、河水的流动、风、雨、气压的变化等,这些现象产生的信号频率小于1HZ,属于低频信号源,这类微动被称为长波微动。
微动没有特定的震源,振动波来自观测点的四面八方,是一种由体波(P波和S波)和面波(瑞利波和拉夫波)组成的复杂振动,其中面波的能量占信号总能量的70%以上,携带有大量的地球物理信息。微动信号的振幅和形态随时空变化而发生变化,但在一定时空范围内具有统计稳定性,可用时间和空间上的平稳随机过程描述。利用上述理论进行勘探的方法称为微动(microtremor)勘探,或称为被动源面波(PassiveSeurfaceWave)勘探,也可称为环境振动(AmbientVibratinom)勘探。微动勘探理论上是以平稳随机过程理论为依据,从微动信号中提取面波(Rayleigh波)的频散曲线,通过对频散曲线进行反演,获得地下介质的横波速度结构。
考虑到当前的微动勘探方式还存在数据采集耗时长、勘探效率低,不同测点常受到不同时间震源差异性的影响,导致勘探结果准确性降低的问题,本发明实施例提供了一种微动勘探方法、装置及电子设备,可以提高勘探效率,并提升勘探结果的准确性。其中,该技术可以应用于各类利用微动勘探进行地质任务探测的应用场景中。
为便于对本实施例进行理解,首先对本发明实施例所公开的一种微动勘探方法进行详细介绍。
实施例一
参见图1,所示为一种微动勘探方法的流程示意图,由图1可见,该方法包括下述步骤:
步骤102:基于预先布设的微动观测装置采集目标测线的原始微动数据;其中,该微动观测装置按预设道间距布设检波器,且布设的采集道数远大于单个测点的基本道数,该基本道数为满足预设探测深度要求的道数。
预先在工区的目标测线上布设检波器,其中,根据工区的勘探设计要求确定道间距,通常可以取1米道间距。为了满足工区勘探的深度要求,需要预先确定出单个测点的基本道数,例如,可以通过施工前的试验来确定合适的基本道数,使得基于该基本道数采集到的微动数据可以反演得到预设勘探深度的地质情况。
在其中一种可能的实施方式中,可以通过下述公式确定该基本道数:
D=K×(N-1)×R (1)
其中,D为预设探测深度,单位为米;N为该基本道数;R为该预设道间距,单位为米;K为常系数,表示该预设探测深度与观测台阵之间的倍数,通常取值3~7。常系数K基于工区的地质情况进行确定,例如,如果是沿海地区,因潮汐产生较强微动信号,则K取值可以偏大;相比之下,对于内陆地区(青藏等地),K取值则相对偏小。
假设当前工区为山东沿海某地某城市道路,工程要求探查道路30米以内的地层情况,并且,需要区分出直径1米左右的异常体(如空洞、凹陷等等),测点间距为1米,且要求微动勘探的道间距为1米。此时,假设常系数K取值为6,则可以基于上述公式(1)推测出单个测点的基本道数N为6道。
按照通常的微动勘探方式,如果基本道数为6道,则一次布设6道采集一个测点的微动数据,然后移动采集道进行下一个测点的数据采集,如此往复,直至采集完整条测线上的测点。而在本实施例中,不同于常规的这种数据采集方式,而是通过一次性布设远大于单个测点的采集道数,例如,上述例子中,可以布设采集道数为24道或者其他远大于6的采集道数,并基于这种超常规道数同时采集目标测线的原始微动数据。
在其中一种可能的实施方式中,如图2所示为一种微动观测装置的布置示意图,图2示出的方式中,该微动观测装置包括主机21、检波器22、大线23,其中,每个检波器22通过连接线24连接到大线23上,大线23连接到主机21,从而可以将各个检波器22检测到的微动信号传输到主机21上。这里,该微动观测装置中的主机21可以采用主流的观测主机,例如WD-1智能微动勘探仪;并且,检波器22是以超常规道数的方式布设,各个检波器22设置在拖拽带25上,拖拽带25通过绳索固定在机动车26上,当需要挪动检波器22到新的测点位置时,可以通过该机动车26带动拖拽带25,以移动检波器22。在其他可能的实施方式中,该机动车26可以是电动车、小型汽车,或者还可以通过非机动车或其他动力带动拖拽带25。
步骤104:基于该基本道数,从上述原始微动数据中确定多个连续测点的微动数据,得到该多个连续测点中每个测点对应的微动数据。
通过步骤S102中超常规道数布设检波器的方式,一次性布设检波器,可以采集到多个连续测点的微动数据,这里,从采集到的原始微动数据中通过抽道集的方式将各个测点的微动数据提取出来,具体的,可以以该工区的基本道数为滑动窗口,以检波器布设区间的第一端点作为该滑动窗口的起点,按预设步长移动该滑动窗口,抽取每个该滑动窗口对应的道集数据,直至该滑动窗口移动到该检波器布设区间的第二端点为止,得到多个道集数据;其中,该滑动窗口的中间点对应为当前测点,该滑动窗口对应的道集数据为该当前测点的微动数据。这样,即提取到多个连续测点的每个测点对应的微动数据。
为了更清楚理解上述抽道集的过程,这里通过实例展示。如图3所示,为一种通过抽道集的方式提取测点微动数据的示意图,图3示出的方式中,基本道数为10道,一次性布设的采集道数为24道(远大于基本道数),测点间距为1米,道间距为1米。当一次性采集上述24道的数据后,可以通过抽道集的方式,提取出15个测点的微动数据。具体方式为:
以基本道数10道为滑动窗口,以检波器布设区间的第一端点(图中标号为1)作为该滑动窗口的起点,按预设步长移动该滑动窗口,抽取每个该滑动窗口对应的道集数据,直至该滑动窗口移动到该检波器布设区间的第二端点(图中标号为24)为止,得到多个道集数据。其中,第一个滑动窗口对应为1至10道(为一个道集数据)、第二个滑动窗口对应为2-11道,依次类推,得到15个道集数据,每个道集数据对应一个测点的数据。这里,1至10道这一道集数据对应的测点位置为图中标号5和标号6的中间位置,在实际操作中,通常是预先设定测线中的测点位置,进而倒推出检波器的布设位置,从而可以以上述超常规道数的方式布设检波器,以实现数据的快速采集。
在其他可能的实施方式中,上述步长的取值通常为一个道间距,也可以取值为多个道间距。并且,上述采集道数的布设通常取值为24的正整数倍,例如取24道、48道、96道等等。
步骤106:根据该多个连续测点中每个测点对应的微动数据,确定该多个连续测点所在区间的地质勘探成果。
在其中一种可能的实施方式中,可以通过下述步骤11-12确定该多个连续测点所在区间的地质勘探成果:
(11)根据该多个连续测点中每个测点对应的微动数据,提取每个该测点的频散曲线。
该频散曲线可以是速度-深度域曲线。
(12)根据提取到的该频散曲线,反演得到该多个连续测点所在区间的地质剖面图。
为了现场及时判断采集数据的质量,在其中一种方式中,还在该微动观测装置的预设显示终端显示该地质剖面图。例如,如果采用WD-1智能微动勘探仪,则可以在该勘探仪的显示界面上呈现上述反演处理得到的地质剖面图,以便于现场工作人员直观判断采集数据的质量情况,决定是否重新采集数据,或者,是否要调整采集参数重采,这种方式可以进一步提高微动勘探的效率,提高数据采集的质量。
在上述步骤S102至S106介绍的微动勘探方式中,通过一次性布设远大于单个测点基本道数的采集道数,以超常规道数的方式布设检波器,从而一次性采集到多个测点的微动数据,并通过抽道的方式提取出多个测点中每个测点对应的微动数据,从而提高了勘探效率;另一方面,由于提取出的多个测点的微动数据的采集时间一致,因而可以保证各个测点的微动数据来源于相同震源,且受到的干扰也较为一致,从而减小了不同测点采集的微动数据之间的背景场差异,进而提升反演结果对地下地层的横向分辨率,提升勘探结果的准确性。
在实际操作中,基于施工条件,如果单条测线的测点数目不是太多的情况下,可以一次性布设覆盖整条测线的检波器,实现剖面式勘探。通过这种超常规道数采集的方式,一次性同时采集整条测线上各个测点的微动数据,并通过抽道集的方式,将各个测点的微动数据分离出来,以进行后期提取频散曲线和反演等处理,这种方式可以极大提高微动勘探效率,并且,由于各测点的微动数据是同一时间采集的,可以降低不同时间震源差异对地层差异性的影响,提高地层的横向分辨率。
在其中一个实施例中,以城市地面道路塌陷的情形为例,在塌陷区的旁侧按照上图2中所示车载方式进行观测装置的布设,布设完毕后再进行微动数据采集;将采集好的数据进行抽道集处理,获得多个微动数据;再对各个微动数据提取频散曲线;最后将各频散连接成如图4所示的地质剖面图。
本发明实施例提供的微动勘探方法,基于预先布设的微动观测装置采集目标测线的原始微动数据;其中,该微动观测装置按预设道间距布设检波器,且布设的采集道数远大于单个测点的基本道数,该基本道数为满足预设探测深度要求的道数;基于该基本道数,从上述原始微动数据中确定多个连续测点的微动数据,得到该多个连续测点中每个测点对应的微动数据;根据该多个连续测点中每个测点对应的微动数据,确定该多个连续测点所在区间的地质勘探成果。该方式提高了微动勘探的效率,提升了微动勘探反演结果对地下地层的横向分辨率,提升勘探结果的准确性。
实施例二
对应于图1中所示的微动勘探方法,本发明实施例还提供了一种微动勘探装置,如图5所示为本发明实施例提供的一种微动勘探装置的结构示意图,由图5可见,该装置包括依次连接的原始微动数据采集模块51、测点微动数据确定模块52和地质勘探成果确定模块53;其中,各个模块的功能如下:
原始微动数据采集模块51,用于基于预先布设的微动观测装置采集目标测线的原始微动数据;其中,该微动观测装置按预设道间距布设检波器,且布设的采集道数远大于单个测点的基本道数,该基本道数为满足预设探测深度要求的道数;
测点微动数据确定模块52,用于基于该基本道数,从该原始微动数据中确定多个连续测点的微动数据,得到该多个连续测点中每个测点对应的微动数据;
地质勘探成果确定模块53,用于根据该多个连续测点中每个测点对应的微动数据,确定该多个连续测点所在区间的地质勘探成果。
本发明实施例提供的微动勘探装置,通过一次性布设远大于单个测点基本道数的采集道数,以超常规道数的方式布设检波器,从而一次性采集到多个测点的微动数据,并通过抽道的方式提取出多个测点中每个测点对应的微动数据,从而提高了勘探效率;另一方面,由于提取出的多个测点的微动数据的采集时间一致,因而可以保证各个测点的微动数据来源于相同震源,且受到的干扰也较为一致,从而减小了不同测点采集的微动数据之间的背景场差异,进而提升反演结果对地下地层的横向分辨率,提升勘探结果的准确性。
在其中一种可能的实施方式中,上述测点微动数据确定模块52还用于:以该基本道数为滑动窗口,以检波器布设区间的第一端点作为该滑动窗口的起点,按预设步长移动该滑动窗口,抽取每个该滑动窗口对应的道集数据,直至该滑动窗口移动到该检波器布设区间的第二端点为止,得到多个道集数据;其中,该滑动窗口的中间点对应为当前测点,该滑动窗口对应的道集数据为该当前测点的微动数据。
在另一种可能的实施方式中,上述步长的取值为一个或多个道间距。
在另一种可能的实施方式中,上述采集道数为24道的正整数倍。
在另一种可能的实施方式中,上述基本道数根据下述公式确定:D=K×(N-1)×R,其中,D为该预设探测深度,N为该基本道数,R为该预设道间距,K为常系数,表示该预设探测深度与观测台阵之间的倍数。
在另一种可能的实施方式中,上述地质勘探成果确定模块53还用于:根据该多个连续测点中每个测点对应的微动数据,提取每个该测点的频散曲线;根据提取到的该频散曲线,反演得到该多个连续测点所在区间的地质剖面图。
在另一种可能的实施方式中,该频散曲线为速度-深度域曲线。
在另一种可能的实施方式中,该装置还包括显示模块,用于在该微动观测装置的预设显示终端显示该地质剖面图。
本发明实施例提供的微动勘探装置,其实现原理及产生的技术效果和前述微动勘探方法实施例相同,为简要描述,微动勘探装置的实施例部分未提及之处,可参考前述微动勘探方法实施例中相应内容。
本发明实施例还提供了一种电子设备,如图6所示,为该电子设备的结构示意图,其中,该电子设备包括处理器61和存储器62,该存储器62存储有能够被该处理器61执行的机器可执行指令,该处理器61执行该机器可执行指令以实现上述微动勘探方法。
在图6示出的实施方式中,该电子设备还包括总线63和通信接口64,其中,处理器61、通信接口64和存储器62通过总线连接。
其中,存储器62可能包含高速随机存取存储器(RAM,Random Access Memory),也可能还包括非不稳定的存储器(non-volatile memory),例如至少一个磁盘存储器。通过至少一个通信接口64(可以是有线或者无线)实现该系统网元与至少一个其他网元之间的通信连接,可以使用互联网,广域网,本地网,城域网等。总线可以是ISA总线、PCI总线或EISA总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示,图6中仅用一个双向箭头表示,但并不表示仅有一根总线或一种类型的总线。
处理器61可能是一种集成电路芯片,具有信号的处理能力。在实现过程中,上述方法的各步骤可以通过处理器61中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器61可以是通用处理器,包括中央处理器(Central Processing Unit,简称CPU)、网络处理器(Network Processor,简称NP)等;还可以是数字信号处理器(Digital SignalProcessing,简称DSP)、专用集成电路(Application Specific Integrated Circuit,简称ASIC)、现成可编程门阵列(Field-Programmable Gate Array,简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成,或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器,闪存、只读存储器,可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器,处理器61读取存储器62中的信息,结合其硬件完成前述实施例的微动勘探方法的步骤。
本发明实施例还提供了一种机器可读存储介质,该机器可读存储介质存储有机器可执行指令,该机器可执行指令在被处理器调用和执行时,该机器可执行指令促使处理器实现上述微动勘探方法,具体实现可参见前述方法实施例,在此不再赘述。
本发明实施例所提供的微动勘探方法、微动勘探装置和电子设备的计算机程序产品,包括存储了程序代码的计算机可读存储介质,所述程序代码包括的指令可用于执行前面方法实施例中所述的微动勘探方法,具体实现可参见方法实施例,在此不再赘述。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个处理器可执行的非易失的计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
另外,在本发明实施例的描述中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
在本发明的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (10)
1.一种微动勘探方法,其特征在于,包括:
基于预先布设的微动观测装置采集目标测线的原始微动数据;其中,所述微动观测装置按预设道间距布设检波器,且布设的采集道数远大于单个测点的基本道数,所述基本道数为满足预设探测深度要求的道数;
基于所述基本道数,从所述原始微动数据中确定多个连续测点的微动数据,得到所述多个连续测点中每个测点对应的微动数据;
根据所述多个连续测点中每个测点对应的微动数据,确定所述多个连续测点所在区间的地质勘探成果。
2.根据权利要求1所述的微动勘探方法,其特征在于,基于所述基本道数,从所述原始微动数据中确定多个连续测点的微动数据,得到所述多个连续测点中每个测点对应的微动数据的步骤,包括:
以所述基本道数为滑动窗口,以检波器布设区间的第一端点作为所述滑动窗口的起点,按预设步长移动所述滑动窗口,抽取每个所述滑动窗口对应的道集数据,直至所述滑动窗口移动到所述检波器布设区间的第二端点为止,得到多个道集数据;
其中,所述滑动窗口的中间点对应为当前测点,所述滑动窗口对应的道集数据为所述当前测点的微动数据。
3.根据权利要求2所述的微动勘探方法,其特征在于,所述步长的取值为一个或多个道间距。
4.根据权利要求1所述的微动勘探方法,其特征在于,所述采集道数为24道的正整数倍。
5.根据权利要求1-4任一项所述的微动勘探方法,其特征在于,所述基本道数根据下述公式确定:
D=K×(N-1)×R
其中,D为所述预设探测深度,N为所述基本道数,R为所述预设道间距,K为常系数,表示所述预设探测深度与观测台阵之间的倍数。
6.根据权利要求1所述的微动勘探方法,其特征在于,根据所述多个连续测点中每个测点对应的微动数据,确定所述多个连续测点所在区间的地质勘探成果的步骤,包括:
根据所述多个连续测点中每个测点对应的微动数据,提取每个所述测点的频散曲线;
根据提取到的所述频散曲线,反演得到所述多个连续测点所在区间的地质剖面图。
7.根据权利要求6所述的微动勘探方法,其特征在于,所述频散曲线为速度-深度域曲线。
8.根据权利要求6所述的微动勘探方法,其特征在于,在根据提取到的所述频散曲线,反演得到所述多个连续测点所在区间的地质剖面图的步骤之后,所述方法还包括:
在所述微动观测装置的预设显示终端显示所述地质剖面图。
9.一种微动勘探装置,其特征在于,包括:
原始微动数据采集模块,用于基于预先布设的微动观测装置采集目标测线的原始微动数据;其中,所述微动观测装置按预设道间距布设检波器,且布设的采集道数远大于单个测点的基本道数,所述基本道数为满足预设探测深度要求的道数;
测点微动数据确定模块,用于基于所述基本道数,从所述原始微动数据中确定多个连续测点的微动数据,得到所述多个连续测点中每个测点对应的微动数据;
地质勘探成果确定模块,用于根据所述多个连续测点中每个测点对应的微动数据,确定所述多个连续测点所在区间的地质勘探成果。
10.一种电子设备,其特征在于,所述电子设备包括处理器和存储器,所述存储器存储有能够被所述处理器执行的计算机可执行指令,所述处理器执行所述计算机可执行指令以实现权利要求1至8任一项所述的微动勘探方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110819615.7A CN113466930B (zh) | 2021-07-20 | 2021-07-20 | 微动勘探方法、装置及电子设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110819615.7A CN113466930B (zh) | 2021-07-20 | 2021-07-20 | 微动勘探方法、装置及电子设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113466930A true CN113466930A (zh) | 2021-10-01 |
CN113466930B CN113466930B (zh) | 2024-06-07 |
Family
ID=77881391
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110819615.7A Active CN113466930B (zh) | 2021-07-20 | 2021-07-20 | 微动勘探方法、装置及电子设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113466930B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114240073A (zh) * | 2021-11-25 | 2022-03-25 | 湖北省地震局(中国地震局地震研究所) | 微动勘探施工方法、系统、计算机设备及存储介质 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10319128A (ja) * | 1997-05-21 | 1998-12-04 | Hiroshi Okada | 受動型地下構造探査装置 |
US20050052951A1 (en) * | 2003-05-30 | 2005-03-10 | Ray Clifford H. | Method and apparatus for seismic data acquisition |
JP2006349491A (ja) * | 2005-06-15 | 2006-12-28 | Central Res Inst Of Electric Power Ind | 微動アレイ探査法を用いた傾斜基盤構造の同定方法 |
JP2017122624A (ja) * | 2016-01-06 | 2017-07-13 | 中部電力株式会社 | 観測地震の距離パラメータを推定する方法ならびに距離パラメータ推定プログラムおよび距離パラメータ推定プログラムを記録したコンピュータ読み取り可能な記録媒体 |
JP6596626B2 (ja) * | 2018-04-19 | 2019-10-30 | 株式会社グリーンデザインオフィス | 常時微動計測による推定柱状図作成方法、プログラム |
CN112505749A (zh) * | 2020-10-19 | 2021-03-16 | 中国地质调查局南京地质调查中心(华东地质科技创新中心) | 一种基于线形台阵多次覆盖的微动数据采集方法 |
CN112526597A (zh) * | 2020-11-09 | 2021-03-19 | 中铁第四勘察设计院集团有限公司 | 一种微动技术的勘探方法、装置和存储介质 |
CN112731551A (zh) * | 2020-12-15 | 2021-04-30 | 核工业北京地质研究院 | 一种砂岩铀矿天窗构造的微动探测方法及系统 |
CN213482468U (zh) * | 2020-10-20 | 2021-06-18 | 北京市水电物探研究所 | 一种多工作模式无线智能微动采集装置 |
RU2750701C1 (ru) * | 2020-10-26 | 2021-07-01 | Федеральное государственное автономное образовательное учреждение высшего образования "Российский университет дружбы народов (РУДН) | Способ геофизической разведки |
-
2021
- 2021-07-20 CN CN202110819615.7A patent/CN113466930B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10319128A (ja) * | 1997-05-21 | 1998-12-04 | Hiroshi Okada | 受動型地下構造探査装置 |
US20050052951A1 (en) * | 2003-05-30 | 2005-03-10 | Ray Clifford H. | Method and apparatus for seismic data acquisition |
JP2006349491A (ja) * | 2005-06-15 | 2006-12-28 | Central Res Inst Of Electric Power Ind | 微動アレイ探査法を用いた傾斜基盤構造の同定方法 |
JP2017122624A (ja) * | 2016-01-06 | 2017-07-13 | 中部電力株式会社 | 観測地震の距離パラメータを推定する方法ならびに距離パラメータ推定プログラムおよび距離パラメータ推定プログラムを記録したコンピュータ読み取り可能な記録媒体 |
JP6596626B2 (ja) * | 2018-04-19 | 2019-10-30 | 株式会社グリーンデザインオフィス | 常時微動計測による推定柱状図作成方法、プログラム |
CN112505749A (zh) * | 2020-10-19 | 2021-03-16 | 中国地质调查局南京地质调查中心(华东地质科技创新中心) | 一种基于线形台阵多次覆盖的微动数据采集方法 |
CN213482468U (zh) * | 2020-10-20 | 2021-06-18 | 北京市水电物探研究所 | 一种多工作模式无线智能微动采集装置 |
RU2750701C1 (ru) * | 2020-10-26 | 2021-07-01 | Федеральное государственное автономное образовательное учреждение высшего образования "Российский университет дружбы народов (РУДН) | Способ геофизической разведки |
CN112526597A (zh) * | 2020-11-09 | 2021-03-19 | 中铁第四勘察设计院集团有限公司 | 一种微动技术的勘探方法、装置和存储介质 |
CN112731551A (zh) * | 2020-12-15 | 2021-04-30 | 核工业北京地质研究院 | 一种砂岩铀矿天窗构造的微动探测方法及系统 |
Non-Patent Citations (6)
Title |
---|
刘云祯;梅汝吾;叶佩;金荣杰: "WD智能天然源面波数据采集处理系统及其应用试验", 物探与化探, vol. 40, no. 05, 31 October 2016 (2016-10-31), pages 1007 - 1015 * |
刘庆华;鲁来玉;王凯明;: "主动源和被动源面波浅勘方法综述", 地球物理学进展, no. 06, 15 December 2015 (2015-12-15), pages 2906 - 2922 * |
张建清;熊永红;李鹏;: "微动勘探技术在城市轨道交通勘察中的应用", 人民长江, no. 01, 14 January 2016 (2016-01-14), pages 39 - 42 * |
李晓斌;杨振威;云美厚;杨双安;: "矿山微地震四线交错观测系统与面波速度成像", 煤炭学报, no. 2, 15 December 2019 (2019-12-15), pages 643 - 649 * |
董耀;高鹏举;金路;李恒;肖娟;: "微动探测在城市地质勘查中的应用研究", 能源与环保, no. 12, pages 88 - 92 * |
贾辉;陈义军;王铁领;高宏伟;陈昌彦;: "浅部地层微动勘探关键影响因素探讨", 工程勘察, no. 09, pages 68 - 73 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114240073A (zh) * | 2021-11-25 | 2022-03-25 | 湖北省地震局(中国地震局地震研究所) | 微动勘探施工方法、系统、计算机设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113466930B (zh) | 2024-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101968550B (zh) | 基于阵列光纤传感器的岩层识别装置及方法 | |
CN107478196B (zh) | 岩土分层沉降测量方法及测量系统 | |
CN108957563B (zh) | 一种隧道施工超前地质探测系统以及探测方法 | |
CN112505749B (zh) | 一种基于线形台阵多次覆盖的微动数据采集方法 | |
Gabàs et al. | Combination of geophysical methods to support urban geological mapping | |
CN115097445B (zh) | 道路路基病害三维探地雷达探测方法、系统、设备及终端 | |
CN206378448U (zh) | 基于综合物探方法的污染土检测测线布置结构 | |
Baptista et al. | A precise and efficient methodology to analyse the shoreline displacement rate | |
CN106706715A (zh) | 基于三维高密度电阻率法的污染土检测方法 | |
CN101762347A (zh) | 一种用半波法测量多跨钢拉索索力的方法 | |
CN106706029A (zh) | 一种面向地下结构施工的土体性能监测装置及其工作方法 | |
CN103205972A (zh) | 一种分析基坑变形与坑外地面沉降关系的方法 | |
CN113419294A (zh) | 一种多维度岩溶特殊地质综合探测方法 | |
CN102928713B (zh) | 一种磁场天线的本底噪声测量方法 | |
CN110579760A (zh) | 一种矿区土壤物理性质的三维呈现方法及系统 | |
CN113466930B (zh) | 微动勘探方法、装置及电子设备 | |
CN102767366B (zh) | 一种高分辨率方位电阻率侧向测井仪及测井方法 | |
CN113433595A (zh) | 基于自然电场隧道裂隙水的超前预报方法 | |
JP6716710B2 (ja) | 鉄塔基礎諸元予測方法 | |
CN114692391B (zh) | 用于模拟岩溶地区地层的相似试验装置及其使用方法 | |
CN115508907A (zh) | 一种车载堤防险情隐患探测系统及预警方法 | |
CN111695783B (zh) | 基于北斗的施工安全信息网络系统 | |
CN114660662A (zh) | 一种油气井下实时监测压裂方法、系统及装置 | |
CN113406707A (zh) | 一种大地电磁多尺度、多时段探测方法 | |
CN113900145A (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 |