CN114791623A - 一种微动采集方法 - Google Patents

一种微动采集方法 Download PDF

Info

Publication number
CN114791623A
CN114791623A CN202210194880.5A CN202210194880A CN114791623A CN 114791623 A CN114791623 A CN 114791623A CN 202210194880 A CN202210194880 A CN 202210194880A CN 114791623 A CN114791623 A CN 114791623A
Authority
CN
China
Prior art keywords
measurement
point
micro
cross
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.)
Pending
Application number
CN202210194880.5A
Other languages
English (en)
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.)
Zhejiang Huadong Construction Engineering Co Ltd
Original Assignee
Zhejiang Huadong Construction Engineering Co Ltd
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 Zhejiang Huadong Construction Engineering Co Ltd filed Critical Zhejiang Huadong Construction Engineering Co Ltd
Priority to CN202210194880.5A priority Critical patent/CN114791623A/zh
Publication of CN114791623A publication Critical patent/CN114791623A/zh
Pending legal-status Critical Current

Links

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/18—Receiving elements, e.g. seismometer, geophone or torque detectors, for localised single point measurements
    • G01V1/181—Geophones
    • 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/24—Recording seismic data
    • G01V1/247—Digital recording of seismic data, e.g. in acquisition units or nodes
    • 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/288—Event detection in seismic signals, e.g. microseismics
    • 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
    • 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/64—Geostructures, e.g. in 3D data cubes

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种微动采集方法，包括如下步骤：S1：通过电缆把检波器连接到地震仪，按照设计，布置检波器；S2：启动进行微动数据测量，得到1个测点或1个排列的微动数据；S3：移动检波器到下一个测点或排列；S4：重复步骤S2和S3的测量；S5：收集每个检波器的测量信号，根据基点和测点的位置坐标计算基点与测点之间的距离；S6：基点和测点的测量数据在频率域做互相关，形成互相关后的道集；S7：对有n个检波器的排列，最少可以提取n‑3个排列的互相关道集；S8：利用互相关道集提取频散曲线，并进行横波波速反演等成果图绘制。本发明通过基点测量，实现工区内线性排列或台阵的互相关，提高台阵间频散曲线计算精度，达到精细探测的目的。

Description

一种微动采集方法
技术领域
本发明涉及一种微动采集方法，主要应用于对地下介质之间存在明显密度、波速差异的地质异常体的探测，适合对地下波速异常体的测量，适用于采空区、断层、空洞、孤石等等浅部空间的精细探测，适应对地下结构、构造探测速度、深度和精要求较高的项目。
背景技术
微动方法是用地震仪采集来自于天然地震、海浪、风、潮汐变化、雷电等自然现象或矿山开采、机械运转等人类活动产生的震动信号，并利用微动信号开展地下结构探测的地球物理方法。
微动探测是通过地震仪记录微动信号的垂直分量计算频率与相速度（频散曲线），并通过对频散曲线的反演获得地下介质的波速结构，达到地下介质的速度重建，随着探测方法和技术的不断发展，已经广泛应用于地质结构探测、地热调查、采空区、岩溶等多个领域的探测，并取得了大量成功的案例。
目前常用的微动观测装置主要有以下几种：嵌套三角形、十字形、L形、圆形、菱形等台阵和直线型。主要存在如下问题和缺点：（1）、台站装置复杂，数据采集效率低下；（2）、台站在地形复杂地区布设困难；（3）、抗固定震动干扰能力差，场地附近的高频振动或噪声对微动原始数据影响较大；（4）、由于微动探测深度与台阵边长大小成正比，地形复杂地区台阵的探测深度受限。
由于装置本身的特点，采用台阵和直线型装置经简单处理频散曲线提取困难，难以获得真实重构地下介质的速度结构。
发明内容
本发明的目的在于，针对现有技术中存在的不足，提供一种微动采集方法。
为此，本发明的上述目的通过如下技术方案实现：
一种微动采集方法，其特征在于：所述微动采集方法包括如下步骤：
S1：通过电缆把检波器连接到地震仪，按照设计，布置检波器；
S2：启动进行微动数据测量，得到1个测点或1个排列的微动数据；
S3：移动检波器到下一个测点或排列；
S4：重复步骤S2和S3的测量；
S5：收集每个检波器的测量信号，根据基点和测点的位置坐标计算基点与测点之间的距离Δx1、Δx2、Δx3、Δx4、…、Δxn，其中：n为检波点编号；
S6：基点和测点的测量数据在频率域做互相关，形成互相关后的道集；
S7：对有n个检波器的排列，最少可以提取n-3个排列的互相关道集；
S8：利用互相关道集提取频散曲线，并进行横波波速反演等成果图绘制。
在采用上述技术方案的同时，本发明还可以采用或者组合采用如下技术方案：
作为本发明的优选技术方案：采用直线型等间距排列，在二维空间中，简化计算，提高频散曲线质量。
作为本发明的优选技术方案：采用节点式地震仪，以提高施工速度，避免数据畸变。
作为本发明的优选技术方案：增加排列中节点的数量，以通过组合不同间距的节点地震仪数据获得不同探测深度的数据，提高成果的精度。
本发明提供一种微动采集方法，通过基点测量（虚拟震源），实现工区内线性排列或台阵的互相关，提高台阵间频散曲线计算精度，达到精细探测的目的。本发明所提供的微动采集方法能够适用于地下空洞、断层等和周围介质有明显波速异常地质体的精细探测。
附图说明
图1为本发明所提供的微动采集方法所涉及的物理架构图。
图2为山区陡坡微动测点布置示意图。
图3为检波器1-3的实际道集数据。
图4为以检波器1为基点，检波器2/3的互相关道集。
图5为以检波器1为基点，检波器2/3的互相关道集所提取的频散曲线。
具体实施方式
本发明提供一种微动采集方法，包括如下技术方案：
1）、采用设置基站的方式进行背景噪声数据的采集；
2）、通过排列与基点的互相关计算，压制固定震动干扰信号；
3）、在台站布设时对每个台站进行精确点位测量，获得相对坐标用于后继台站间频散曲线的提取；
4）、可以根据测深要求和台站数量进行台站滚动测量，在每个时间段测量中，基点与测量点之间进行互相关计算，因此，在一个工区内，基点可以随排列移动，也可以固定不动。
具体地，参照附图和具体实施例对本发明作进一步详细地描述。
一种微动采集方法，其特征在于：所述微动采集方法包括如下步骤：
S1：通过电缆把检波器连接到地震仪，按照设计，布置检波器（或节点地震仪）；
S2：启动进行微动数据测量，得到1个测点或1个排列的微动数据；
S3：移动检波器（或节点地震仪）到下一个测点或排列；
S4：重复步骤S2和S3的测量；
S5：收集每个检波器（或节点地震仪）的测量信号，根据基点和测点的位置坐标计算基点与测点之间的距离Δx1、Δx2、Δx3、Δx4、…、Δxn，其中：n为检波点编号；
S6：基点和测点的测量数据在频率域做互相关，形成互相关后的道集；
图2为山区陡坡上的排列布置示意图，图中可以用检波器1作为基点与检波器2和检波器3做互相关，得到相对于基点检波器1的2个道集，用于提取频散曲线；用检波器2作为基点与检波器3和检波器4做互相关，得到相对于基点检波器2的2个道集，用于提取频散曲线；依次类推，最后以检波器11作为基点与检波器12和检波器13做互相关，得到相对于基点检波器11的2个道集，用于提取频散曲线；
S7：对有n个检波器（或节点地震仪）的排列，最少可以提取n-3个排列的互相关道集；
S8：利用互相关道集提取频散曲线，并进行横波波速反演等成果图绘制。
作为优选，建议采用直线型等间距排列，在二维空间中，简化计算，提高频散曲线质量。
作为优选，建议采用节点式地震仪，可以提高施工速度，避免因地震电缆漏电等引起的数据畸变。
作为优选，建议增加排列中节点的数量，可以通过组合不同间距的节点地震仪数据获得不同探测深度的数据，提高成果的精度。
图3的横坐标为时间轴，纵坐标为检波器编号。图3中波形图为采用本发明装置选取的3个连续布置的检波器的道集数据。图3中波形图上连续波形是用于数据分析的基础，对n个检波器的排列，最少可以提取n-3个这样的道集数据。
图4的横坐标为时间轴，纵坐标为检波器编号。图4是图3数据采用互相关算法计算得到的道集数据。以第1道作为基准，第2道与第1道对比做相关，第3道与第1道对比做相关，得到的相关数据的波形图。
图5是频散曲线图，图中横坐标为相速度，纵坐标为提取频率。图5是对图4道集进行空间互相关计算的频散曲线。
上述具体实施方式用来解释说明本发明，仅为本发明的优选实施例，而不是对本发明进行限制，在本发明的精神和权利要求的保护范围内，对本发明做出的任何修改、等同替换、改进等，都落入本发明的保护范围。

Claims (4)

1.一种微动采集方法，其特征在于：所述微动采集方法包括如下步骤：
S1：通过电缆把检波器连接到地震仪，按照设计，布置检波器；
S2：启动进行微动数据测量，得到1个测点或1个排列的微动数据；
S3：移动检波器到下一个测点或排列；
S4：重复步骤S2和S3的测量；
S5：收集每个检波器的测量信号，根据基点和测点的位置坐标计算基点与测点之间的距离Δx1、Δx2、Δx3、Δx4、…、Δxn，其中：n为检波点编号；
S6：基点和测点的测量数据在频率域做互相关，形成互相关后的道集；
S7：对有n个检波器的排列，最少可以提取n-3个排列的互相关道集；
S8：利用互相关道集提取频散曲线，并进行横波波速反演等成果图绘制。
2.根据权利要求1所述的微动采集方法，其特征在于：采用直线型等间距排列，在二维空间中，简化计算，提高频散曲线质量。
3.根据权利要求1所述的微动采集方法，其特征在于：采用节点式地震仪，以提高施工速度，避免数据畸变。
4.根据权利要求1所述的微动采集方法，其特征在于：增加排列中节点的数量，以通过组合不同间距的节点地震仪数据获得不同探测深度的数据，提高成果的精度。
CN202210194880.5A 2022-03-01 2022-03-01 一种微动采集方法 Pending CN114791623A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210194880.5A CN114791623A (zh) 2022-03-01 2022-03-01 一种微动采集方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210194880.5A CN114791623A (zh) 2022-03-01 2022-03-01 一种微动采集方法

Publications (1)

Publication Number Publication Date
CN114791623A true CN114791623A (zh) 2022-07-26

Family

ID=82460584

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210194880.5A Pending CN114791623A (zh) 2022-03-01 2022-03-01 一种微动采集方法

Country Status (1)

Country Link
CN (1) CN114791623A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115755167A (zh) * 2022-12-09 2023-03-07 中铁第四勘察设计院集团有限公司 一种自适应瞬态面波勘探方法、设备及存储介质
CN116400406A (zh) * 2023-04-21 2023-07-07 中国地震局地球物理研究所 一种基于阵列的被动源多模式面波频散曲线提取方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110255371A1 (en) * 2009-01-09 2011-10-20 Charlie Jing Hydrocarbon Detection With Passive Seismic Data
CN111352153A (zh) * 2020-03-23 2020-06-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 核工业北京地质研究院 一种砂岩铀矿天窗构造的微动探测方法及系统
CN113900145A (zh) * 2021-09-30 2022-01-07 福建省交通规划设计院有限公司 一种基于微动信号的地下横波速度分布的推测方法和装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110255371A1 (en) * 2009-01-09 2011-10-20 Charlie Jing Hydrocarbon Detection With Passive Seismic Data
CN111352153A (zh) * 2020-03-23 2020-06-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 核工业北京地质研究院 一种砂岩铀矿天窗构造的微动探测方法及系统
CN113900145A (zh) * 2021-09-30 2022-01-07 福建省交通规划设计院有限公司 一种基于微动信号的地下横波速度分布的推测方法和装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘玉海等: "基于互相关函数法的地面微地震信号检测研究", 《石油物探》, vol. 51, no. 06, 30 November 2012 (2012-11-30), pages 633 - 637 *
张建清等: "微动勘探技术在城市轨道交通勘察中的应用", 《人民长江》, vol. 47, no. 01, 31 January 2016 (2016-01-31), pages 39 - 42 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115755167A (zh) * 2022-12-09 2023-03-07 中铁第四勘察设计院集团有限公司 一种自适应瞬态面波勘探方法、设备及存储介质
WO2024119966A1 (zh) * 2022-12-09 2024-06-13 中铁第四勘察设计院集团有限公司 一种自适应瞬态面波勘探方法、设备及存储介质
CN116400406A (zh) * 2023-04-21 2023-07-07 中国地震局地球物理研究所 一种基于阵列的被动源多模式面波频散曲线提取方法
CN116400406B (zh) * 2023-04-21 2023-12-19 中国地震局地球物理研究所 一种基于阵列的被动源多模式面波频散曲线提取方法

Similar Documents

Publication Publication Date Title
Zhu et al. Characterizing thunder‐induced ground motions using fiber‐optic distributed acoustic sensing array
Uebayashi Extrapolation of irregular subsurface structures using the horizontal-to-vertical spectral ratio of long-period microtremors
CN112505749B (zh) 一种基于线形台阵多次覆盖的微动数据采集方法
CN113900145B (zh) 一种基于微动信号的地下横波速度分布的推测方法和装置
CN102879805A (zh) 一种基于钻孔与地面相结合的地震波空间探测方法
CN114791623A (zh) 一种微动采集方法
CN111025380A (zh) 一种矿井微地震观测系统设计及面波提取方法
Coccia et al. Application of Refraction Microtremor (ReMi) technique for determination of 1-D shear wave velocity in a landslide area
CN114924315A (zh) 一种多态式地震勘探方法和系统
CN112180444A (zh) 一种地层速度结构的探测方法、装置以及存储介质
Asano et al. Surface wave group velocity in the Osaka sedimentary basin, Japan, estimated using ambient noise cross-correlation functions
CN110850469A (zh) 一种基于克希霍夫积分解的地震槽波深度偏移的成像方法
Jiang et al. Distributed acoustic sensing for shallow structure imaging using mechanical noise: A case study in Guangzhou, China
Liu et al. Incorporating geotechnical and geophysical investigations for underground obstruction detection: A case study
JP6531934B2 (ja) ハイブリッド表面波探査方法及びハイブリッド表面波探査システム
Guan et al. Improving the retrieval of high-frequency surface waves using convolution-based three-station interferometry for dense linear arrays
CN105093314A (zh) 一种测定微地震震源的方法
Wu et al. Ambient noise surface wave tomography of Quaternary structures derived from a high-density array in the central Hebei Depression, North China
Lin et al. Modelling of pulse-like velocity ground motion during the 2018 M w 6.3 Hualien earthquake, Taiwan
Ikeda et al. Development of surface-wave monitoring system for leaked CO2 using a continuous and controlled seismic source
Guan et al. 3-D tunnel seismic advance prediction method with wide illumination and high-precision
CN116338779B (zh) 用于密集线性台阵的交火式主被动源联合探测方法
US11933928B1 (en) Forward simulation-based irregular seismic data acquisition method
CN117687118A (zh) 一种基于光纤微动的地下空区综合定位勘察方法
Taipodia et al. A review of active and passive MASW techniques

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