CN110488349A - 基于微动三分量谱比vhsr的无损探测方法及应用 - Google Patents
基于微动三分量谱比vhsr的无损探测方法及应用 Download PDFInfo
- Publication number
- CN110488349A CN110488349A CN201910769450.XA CN201910769450A CN110488349A CN 110488349 A CN110488349 A CN 110488349A CN 201910769450 A CN201910769450 A CN 201910769450A CN 110488349 A CN110488349 A CN 110488349A
- Authority
- CN
- China
- Prior art keywords
- vhsr
- component
- curve
- fine motion
- detection method
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 claims abstract description 42
- 238000010276 construction Methods 0.000 claims abstract description 24
- 230000000694 effects Effects 0.000 claims abstract description 17
- 230000002159 abnormal effect Effects 0.000 claims abstract description 10
- 239000011248 coating agent Substances 0.000 claims abstract description 6
- 238000000576 coating method Methods 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 4
- 206010044565 Tremor Diseases 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 1
- 239000010410 layer Substances 0.000 description 16
- 239000002689 soil Substances 0.000 description 13
- 230000004888 barrier function Effects 0.000 description 8
- 239000004927 clay Substances 0.000 description 6
- 230000002547 anomalous effect Effects 0.000 description 5
- 239000004576 sand Substances 0.000 description 4
- 241001270131 Agaricus moelleri Species 0.000 description 3
- 238000009412 basement excavation Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 239000000843 powder Substances 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 208000002177 Cataract Diseases 0.000 description 1
- 208000010412 Glaucoma Diseases 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000002247 constant time method Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000010291 electrical method Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 235000013312 flour Nutrition 0.000 description 1
- 239000010438 granite Substances 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000011499 joint compound Substances 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000003466 welding Methods 0.000 description 1
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
-
- 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
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
本发明为一种基于微动三分量谱比VHSR的无损探测方法及应用。方法包括以下步骤：1、采集二维阵列的地表微震动信号；2、分别计算各台三分量检波器所采集数据的VHSR曲线；3、用自研MATLAB程序将各测点的VHSR曲线形成二维VHSR曲线等值线图或者三维切片图；4、根据二维VHSR曲线等值线图或者三维切片图的异常特征推断解释异常目标体。本方法用于在城市工程建设中浅地表或覆盖层深度范围内桩基、承台、管道、人工填石、空洞、软弱夹层等对工程建设有影响的地下障碍物或不良地质的探测。本发明能准确查明地下障碍物和不良地质的分布范围，为盾构、顶管等城市工程建设施工提前发现潜在风险，避免发生重大施工事故，作用重大。
Description
技术领域
本发明涉及一种工程地球物理探测方法，尤其涉及城市工程建设中浅地表或覆盖层深度范围内桩基、承台、管道、人工填石、空洞、软弱夹层等对工程建设有影响的地下障碍物或不良地质探测的基于微动三分量谱比VHSR的无损探测方法及应用。
背景技术
城市工程建设中浅地表或覆盖层深度范围内桩基、承台、管道、人工填石、空洞、软弱夹层等经常对工程建设造成影响或危害。比如盾构机掘进中如果遇到未探明的桩基，经常会严重损坏盾构机，甚至会造成喷涌、塌方等意外情况；软弱夹层或者空洞可能造成盾构机发生“磕头”现象，使盾构机严重偏离预定轨道，无法正常掘进，甚至造成严重施工事故；未探明的承台、管道、人工填石等也会妨碍顶管施工，会造成顶管施工无法进行，重新设计施工线路无疑会产生巨大的成本等。被动的处理这些障碍物也会造成破坏环境、延误工期、堵塞路面交通，增加投资控制的不确定性等影响，造成的损失往往是巨大的。因此，在施工开始之前探明这些障碍物是非常有必要的。
地下障碍物或者不良地质的探测一般是直接开挖、钻探和地球物理探测等。障碍物位置不明的情况下采用直接开挖或者钻探，成本高、效果不好，特别在城市复杂环境下，基本不具备开挖和钻探施工条件且风险巨大。因此，寻找一种适应城市复杂环境的地球物理方法来探测这些障碍物或不良地质势在必行。
物探方法种类繁多，电磁波类方法受城市复杂电磁背景干扰严重，常规地震方法及直流电法类受城市狭小地面空间、复杂地表情况等影响较大。探测障碍物常用的如探地雷达，在南方地下水位较浅的城市，信号衰减较快，探测深度受限制，同时也受电磁干扰，效果往往不理想；各类跨孔CT方法有一定效果，但孔间不同方法的效果也不尽相同,且受到孔间距的影响，不仅成本高，工期长。总体上，诸多物探方法在城市复杂的环境下均达不到预期效果。
发明内容
本发明的目的在于提供一种基于微动三分量谱比VHSR的无损探测方法及应用。
本发明的目的通过如下技术方案实现：一种基于微动三分量谱比VHSR的无损探测方法，它包括以下步骤：
步骤1：采用二维小间距的方式布置三分量检波器，并采集地表微震动信号，采集记录时间不少于5min；
步骤2：分别计算各台三分量检波器所采集数据的VHSR曲线；它包括以下流程：
流程(1).对采集的原始数据进行基线漂移校正预处理；
流程(2).去除短时强干扰信号，获取相对平稳的微动信号；
流程(3).对信号进行分窗处理；
流程(4).利用快速傅里叶变换得到分窗的傅里叶振幅谱数据；
流程(5).对分窗处理的傅里叶振幅谱数据，使用统计的方式分别得到各三分量检波器的垂直分量和N-S向与E-W向两个水平分量的傅里叶振幅谱数据；
流程(6).对N-S向与E-W向两个水平分量的傅里叶振幅谱数据进行计算得到合并的水平分量傅里叶振幅谱数据；
流程(7).对垂直分量和合并的水平分量的傅里叶振幅谱数据进行直接比值计算，获得VHSR曲线。
步骤3：用自研的MATLAB程序将各测点的VHSR曲线形成二维VHSR曲线等值线图或者三维切片图；
步骤4：根据二维VHSR曲线等值线图或者三维切片图的异常特征推断解释异常目标体。
所述基于微动三分量谱比VHSR的无损探测方法的应用，用于在城市工程建设中浅地表或覆盖层深度范围内桩基、承台、管道、人工填石、空洞、软弱夹层等对工程建设有影响的地下障碍物或不良地质的探测。
较之现有技术而言，本发明的优点在于：
1、本发明采用了人类活动产生的震动作为震源，主要利用其中人类活动产生的震动包括车辆移动、工厂机械运行，甚至人的行走等频率大于1Hz的震源，城市环境中这类震源相当丰富。人类活动密集的城市环境不仅不会形成干扰，相反为浅层微动勘探提供了丰富的高频信号源；通过以日常人类活动产生的“噪音”作为震源，使得谱比VHSR在城市复杂环境中具备强抗干扰能力，而且采用谱比VHSR，能够放大土层中细微的信号异常，从而达到分辨出障碍物的目的；微动勘探不需要人工激发的震源，不会对周边环境产生影响，对地表条件要求高，不会妨碍交通，只需在测试时进行灵活的临时交通布控，完成测试后，可立即撤走布控设施；
2、采用VHSR谱比法，并且可选择将VHSR量值进行取对数计算，能够放大土层中细微的信号异常，更有效的反映出土层中障碍物或不良地质的位置；微动垂直-水平谱比(VHSR曲线)是一个随频率变化的函数，与台阵下方浅层岩土结构有密切的关系。对于均匀层状介质地层，纵向随深度变化的分层模型来说，这一函数只与基阶瑞利面波所到达位置的地下介质力学参数有关，与其他参数(如震源)无关。VHSR曲线函数对浅层结构尤其是横向不均匀异常体更加敏感，如覆盖层中若有桩基、箱涵，震动传播时水平、垂直的能量会发生变化，在埋深对应地频率段，VHSR曲线会产生明显的异常变化。
3、采用二维小间距的方式布置三分量检波器，可以绘制三维的YZ向、XZ向和YZ向的等频率切片图，能够直观呈现出障碍物的三维分布范围。VHSR曲线的本质为震动在不同岩土层结构中传播时垂直方向和水平方向的能量分配占比会发生变化，不同的土层结构观测到的VHSR曲线的形态量值各不相同，即能量的占比分配能够直接反映地下岩土层的结构变化，这也是VHSR曲线能够用于探测地下岩土层中异常体的直接依据。工程实例表明，一般震动传播遇到较硬的介质，垂直分量的能量会明显的减弱，导致VHSR的量值变小；遇到较软的介质，垂直分量的能量会明显的增强，导致VHSR的量值变大，这种现象可以用于判断覆盖层中的异常体，尤其是土层中的桩基、箱涵等横向变化对比明显的异常地质体。
附图说明
图1是本发明步骤2中获取VHSR曲线流程图。
图2是实施案例一三分量检波器布置示意图。
图3是实施案例一VHSR剖面图。
图4是实施案例一VHSR等频率切片图(2-18Hz)。
图5是实施案例一VHSR等频率切片图(20-45Hz)。
图6是实施案例二VHSR剖面图。
具体实施方式
下面结合说明书附图和实施案例对本发明内容进行详细说明：
一种基于微动三分量谱比VHSR的无损探测方法，它包括以下步骤：
步骤1：采用二维小间距的方式布置三分量检波器，并采集地表微震动信号，采集记录时间不少于5min；这里的三分量检波器采用三分量高精度一体化数字地震仪。
根据掌握的探测目标体的大致情况，在地表上以一定的间距布置二维的三分量检波器，确保布置的三分量检波器能覆盖到探测目标体(可分几次布置三分量检波器，能适应需要较大面积的探测，且三分量检波器数量不足的情况)，并采集地表微震动信号。
二维布置可以采用等间距的规则长方形或者不等间距的长方形，应考虑方便记录各三分量检波器的安放位置。
二维小间距的方式是指根据探测目标体的尺寸采用50公分到5米的三分量检波器间距，确保三分量检波器覆盖到地下目标体，并按方便记录各三分量检波器安放位置的二维平面布置三分量检波器。
步骤2：分别计算各台三分量检波器所采集数据的VHSR曲线；它包括以下流程：
流程(1).对采集的原始数据进行基线漂移校正预处理，具体顺序如下：
顺序1.将原始数据按每2000个采样点分为一段；顺序2.分别对每段数据进行polyfit多项式拟合，得到拟合结果；顺序3.将每一段的原始数据减去该段的拟合结果，得到基线漂移校正后的结果数据。
流程(2).去除短时强干扰信号，获取相对平稳的微动信号；
流程(3).对信号进行分窗处理，具体如下：
分窗的时间长度T，设置为周期的整数倍，其中周期为频率的倒数；
式中，n为整数(一般取10～30倍)，f为频率。
因此不同的频率分窗时间长度是不一样的，同一时间长度的信号，不同频率的分窗个数也不一样。
流程(4).利用快速傅里叶变换得到分窗的傅里叶振幅谱数据；
流程(5).对分窗处理的傅里叶振幅谱数据，使用统计的方式分别得到各三分量检波器的垂直分量和N-S向与E-W向两个水平分量的傅里叶振幅谱数据；
流程(6).对N-S向与E-W向两个水平分量的傅里叶振幅谱数据进行计算得到合并的水平分量傅里叶振幅谱数据，具体采用以下公式：
式中，H为合并的水平分量，NS为南北向的水平分量，EW为东西向的水平分量，f为频率。
流程(7).对垂直分量和合并的水平分量的傅里叶振幅谱数据进行直接比值计算，获得VHSR曲线。
步骤3：用自研的MATLAB程序将各测点的VHSR曲线形成二维VHSR曲线等值线图或者三维切片图，具体如下：用自研的MATLAB程序将各测点的VHSR曲线形成二维VHSR曲线等值线图或者三维切片图时，可以选择对VHSR量值进行取对数计算，VHSR量值取对数后成图可以突出异常范围，使异常的轮廓边界更清楚，有助于准确分析障碍物的位置。
VHlog(f)＝log10(VH(f))
式中，VHlog为取对数后的VHSR的量值，VH为取对数之前的VHSR的量值，f为频率。
步骤4：根据二维VHSR曲线等值线图或者三维切片图的异常特征推断解释异常目标体。
本发明所述的方法主要应用于在城市工程建设中浅地表或覆盖层深度范围内障碍物或不良地质的探测。探测内容包括桩基、承台、管道、人工填石、空洞、软弱夹层等。
下面结合实施案例对本发明进行说明
实施例1：对某地地铁盾构区间桩基位置及桩长VHSR法探测
工程概况：障碍物调查资料显示某地地铁盾构区间某天桥(地面以上部分早已拆除)桩基位于盾构隧道洞身范围内，但存有疑问，因此采用VHSR法进行探测查明。附近详勘钻孔MDZ3-29-34终孔47.8m未揭示基岩面，自上而下主要土层为杂填土、粉质粘土、淤泥、粉质粘土、淤泥、粉质粘土、砂土强风化花岗岩。
步骤1：本次VHSR法探测采用三分量高精度一体化数字地震仪,按0.5～1.5m的间距布置二维三分量检波器阵列(见图2)，采集地表微震动信号10～15min；
步骤2：计算各台三分量检波器所采集数据的VHSR曲线；
步骤3：用自研的MATLAB程序将各测点的VHSR曲线形成二维VHSR曲线等值线图和三维切片图(见图3～5)；
步骤4：探测成果分析解释：
如图3所示：图中横坐标为距离；纵坐标为频率(Hz)，高频反映浅部的信息，低频反映深部的信息；不同灰度代表VHSR量值的大小。剖面中间垂直方框圈定位置VHSR量值小的范围为桩基的位置。
如图4和5所示：横坐标为X方向的距离，纵坐标为Y方向的距离，不同灰度代表VHSR量值的大小，每张图为1个频率切片，高频反映浅部的信息，低频反映深部的信息。因此2～45Hz的频率切片反映了深部到浅部的信息，共20张等频率切片图(出2Hz的切片外，深度较深)中坐标(0,0)位置附近显示出相对较小的VHSR量值，反映的是桩基的位置；可以看出3Hz较低频时，桩基位置仍然有明显的反映，表明桩长较长，桩底埋藏较深，根据目标体的反映深度为半波长来计算，当3Hz频点相速度取值为300m/s时，估算桩底埋深为50m，这与设计的桩长基本吻合。
实施例2：对某建筑场地基坑支护范围人工填石范围VHSR法探测
工程概况：勘察资料揭露测区场地主要岩土层分布特征如下：(1)杂填土(局部粉质粘土)、淤泥(质土)夹砂、粉砂、粉质粘土、淤泥(质土)、粉质粘土、圆角砾、碎卵石、淤泥质土、砂土状强风化、碎块状强风化和中风化凝灰岩等，基岩埋深45m左右。测区边上为规划道路，施工时局部填石可能侵入基坑支护围护桩的施工范围内，会影响到后期的基坑支护围护桩施工作业，为查明基坑支护围护桩施工范围内可能存在的填石等不良地质体，采用无损的VHSR法进行探测。
步骤1：本次VHSR法探测采用三分量高精度一体化数字地震仪,沿着基坑支护围护桩施工范围以5m的台阵间距布置2.5m半径的圆形台阵，共布置34个圆形台阵，采集时间20min；
步骤2：计算各台三分量检波器所采集数据的VHSR曲线；
步骤3：用自研的MATLAB程序将各测点的VHSR曲线形成二维VHSR曲线等值线图(见图6)；
步骤4：探测成果分析解释：
如图3所示：图中横坐标为距离；纵坐标为频率(Hz)，高频反映浅部的信息，低频反映深部的信息；不同灰度代表VHSR量值的大小。剖面中W2～W7测点下方30～45Hz范围VHSR量值明显偏小，推断为填石的位置，估算埋深在0.5-2m左右；W21～W25测点下方15～35Hz范围VHSR量值明显偏小，推断为填石的位置，估算埋深在1.6-4.3m左右。施工方在基坑支护围护桩施工之前，直接进行开挖，证实查明填石位置准确。
Claims (8)
1.一种基于微动三分量谱比VHSR的无损探测方法，其特征在于：它包括以下步骤：
步骤1：采用二维小间距的方式布置三分量检波器，并采集地表微震动信号，采集记录时间不少于5min；
步骤2：分别计算各台三分量检波器所采集数据的VHSR曲线；它包括以下流程：
流程(1).对采集的原始数据进行基线漂移校正预处理；
流程(2).去除短时强干扰信号，获取相对平稳的微动信号；
流程(3).对信号进行分窗处理；
流程(4).利用快速傅里叶变换得到分窗的傅里叶振幅谱数据；
流程(5).对分窗处理的傅里叶振幅谱数据，使用统计的方式分别得到各三分量检波器的垂直分量和N-S向与E-W向两个水平分量的傅里叶振幅谱数据；
流程(6).对N-S向与E-W向两个水平分量的傅里叶振幅谱数据进行计算得到合并的水平分量傅里叶振幅谱数据；
流程(7).对垂直分量和合并的水平分量的傅里叶振幅谱数据进行直接比值计算，获得VHSR曲线。
步骤3：用MATLAB程序将各测点的VHSR曲线形成二维VHSR曲线等值线图或者三维切片图；
步骤4：根据二维VHSR曲线等值线图或者三维切片图的异常特征推断解释异常目标体。
2.根据权利要求1所述的基于微动三分量谱比VHSR的无损探测方法，其特征在于：步骤1中二维小间距的方式是指根据探测目标体的尺寸采用50公分到5米的三分量检波器间距，确保三分量检波器覆盖到地下目标体，并按方便记录各三分量检波器安放位置的二维平面布置三分量检波器，采集时间在5-15min。
3.根据权利要求1所述的基于微动三分量谱比VHSR的无损探测方法，其特征在于：流程(1)对原始数据进行基线漂移校正顺序如下：顺序1.将原始数据按每2000个采样点分为一段；顺序2.分别对每段数据进行polyfit多项式拟合，得到拟合结果；顺序3.将每一段的原始数据减去该段的拟合结果，得到基线漂移校正后的结果数据。
4.根据权利要求1所述的基于微动三分量谱比VHSR的无损探测方法，其特征在于：流程(3)中对信号进行分窗处理，分窗的时间长度T，设置为周期的整数倍，其中周期为频率的倒数；
式中，n为整数，f为频率，其中N取值为10～30。
5.根据权利要求1所述的基于微动三分量谱比VHSR的无损探测方法，其特征在于：流程(6)中对N-S向与E-W向两个水平分量的傅里叶振幅谱数据进行计算得到合并的水平分量傅里叶振幅谱数据，采用以下公式：
式中，H为合并的水平分量，NS为南北向的水平分量，EW为东西向的水平分量，f为频率。
6.根据权利要求1所述的基于微动三分量谱比VHSR的无损探测方法，其特征在于：步骤3中用自研的MATLAB程序将各测点的VHSR曲线形成二维VHSR曲线等值线图或者三维切片图时，选择对VHSR量值进行取对数计算，VHSR量值取对数后成图突出异常范围，使异常的轮廓边界更清楚，有助于准确分析障碍物的位置。
VHlog(f)＝log10(VH(f))
式中，VHlog为取对数后的VHSR的量值，VH为取对数之前的VHSR的量值，f为频率。
7.根据权利要求1所述的基于微动三分量谱比VHSR的无损探测方法，其特征在于：所述三分量检波器采用三分量高精度一体化数字地震仪。
8.根据权利要求1-7中任意一项所述基于微动三分量谱比VHSR的无损探测方法的应用，其特征在于：用于在城市工程建设中浅地表或覆盖层深度范围内桩基、承台、管道、人工填石、空洞、软弱夹层等对工程建设有影响的地下障碍物或不良地质的探测。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910769450.XA CN110488349B (zh) | 2019-08-20 | 2019-08-20 | 基于微动三分量谱比vhsr的无损探测方法及应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910769450.XA CN110488349B (zh) | 2019-08-20 | 2019-08-20 | 基于微动三分量谱比vhsr的无损探测方法及应用 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110488349A true CN110488349A (zh) | 2019-11-22 |
| CN110488349B CN110488349B (zh) | 2021-11-02 |
Family
ID=68552296
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910769450.XA Active CN110488349B (zh) | 2019-08-20 | 2019-08-20 | 基于微动三分量谱比vhsr的无损探测方法及应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110488349B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112558149A (zh) * | 2020-11-12 | 2021-03-26 | 中铁第四勘察设计院集团有限公司 | 基于非均匀介质的微动探测方法、装置、设备及存储介质 |
| CN114137609A (zh) * | 2021-11-09 | 2022-03-04 | 长江地球物理探测（武汉）有限公司 | 线性微动数据校正方法及装置 |
| CN118501939A (zh) * | 2024-05-20 | 2024-08-16 | 中国地质科学院 | 风化壳稀土资源评估方法、系统、设备、介质及程序产品 |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100027377A1 (en) * | 2006-02-24 | 2010-02-04 | Hannes Georges Zuercher | Locating oil or gas actively by exciting a porous oil and gas saturated system to give off its characteristic resonance response, with optional differentiation of oil, gas and water |
| CN101950034A (zh) * | 2010-08-16 | 2011-01-19 | 中国石油天然气股份有限公司 | 一种用于地震反演的自然电位曲线基线漂移校正方法 |
| CN102539467A (zh) * | 2010-12-10 | 2012-07-04 | 中国科学院上海微系统与信息技术研究所 | 一种分析相变材料结晶速率和结晶温度的方法 |
| CN102663768A (zh) * | 2012-05-10 | 2012-09-12 | 合肥工业大学 | 偏振视觉系统中自适应确定偏振镜的最佳透振方向的方法 |
| CN102737246A (zh) * | 2012-06-14 | 2012-10-17 | 公安部天津消防研究所 | 基于Canny算子的泡沫边界识别和粒度分析方法 |
| CN102508311B (zh) * | 2011-11-30 | 2013-10-16 | 河北煤炭科学研究院 | 巷道超前探测数据的多参数空间成图法 |
| CN103424777A (zh) * | 2013-07-01 | 2013-12-04 | 中国科学院地质与地球物理研究所 | 一种提高地震成像分辨率的方法 |
| CN104688234A (zh) * | 2015-03-17 | 2015-06-10 | 吉林大学 | 一种icg色素浓度谱的无创、抗扰动检测方法 |
| CN107065014A (zh) * | 2017-03-14 | 2017-08-18 | 中交公局厦门工程有限公司 | 一种基于微动探测结合地质钻探的地铁区间孤石群探测方法 |
| CN108152854A (zh) * | 2017-11-30 | 2018-06-12 | 福建省建筑设计研究院有限公司 | 一种基于微动功率谱密度的无损探测方法及其应用 |
| CN108318918A (zh) * | 2017-12-29 | 2018-07-24 | 福建省建筑设计研究院有限公司 | 基于微动频散曲线和h/v曲线的地下不良地质体无损探测方法及应用 |
| CN109257128A (zh) * | 2018-11-01 | 2019-01-22 | 南京邮电大学 | 一种基于傅里叶级数拟合去噪的频谱信号识别方法及系统 |
| CN109632016A (zh) * | 2019-02-20 | 2019-04-16 | 湖北理工学院 | 岩土体硐室开挖与围岩应力、应变监测试验装置及其方法 |
-
2019
- 2019-08-20 CN CN201910769450.XA patent/CN110488349B/zh active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100027377A1 (en) * | 2006-02-24 | 2010-02-04 | Hannes Georges Zuercher | Locating oil or gas actively by exciting a porous oil and gas saturated system to give off its characteristic resonance response, with optional differentiation of oil, gas and water |
| CN101950034A (zh) * | 2010-08-16 | 2011-01-19 | 中国石油天然气股份有限公司 | 一种用于地震反演的自然电位曲线基线漂移校正方法 |
| CN102539467A (zh) * | 2010-12-10 | 2012-07-04 | 中国科学院上海微系统与信息技术研究所 | 一种分析相变材料结晶速率和结晶温度的方法 |
| CN102508311B (zh) * | 2011-11-30 | 2013-10-16 | 河北煤炭科学研究院 | 巷道超前探测数据的多参数空间成图法 |
| CN102663768A (zh) * | 2012-05-10 | 2012-09-12 | 合肥工业大学 | 偏振视觉系统中自适应确定偏振镜的最佳透振方向的方法 |
| CN102737246A (zh) * | 2012-06-14 | 2012-10-17 | 公安部天津消防研究所 | 基于Canny算子的泡沫边界识别和粒度分析方法 |
| CN103424777A (zh) * | 2013-07-01 | 2013-12-04 | 中国科学院地质与地球物理研究所 | 一种提高地震成像分辨率的方法 |
| CN104688234A (zh) * | 2015-03-17 | 2015-06-10 | 吉林大学 | 一种icg色素浓度谱的无创、抗扰动检测方法 |
| CN107065014A (zh) * | 2017-03-14 | 2017-08-18 | 中交公局厦门工程有限公司 | 一种基于微动探测结合地质钻探的地铁区间孤石群探测方法 |
| CN108152854A (zh) * | 2017-11-30 | 2018-06-12 | 福建省建筑设计研究院有限公司 | 一种基于微动功率谱密度的无损探测方法及其应用 |
| CN108318918A (zh) * | 2017-12-29 | 2018-07-24 | 福建省建筑设计研究院有限公司 | 基于微动频散曲线和h/v曲线的地下不良地质体无损探测方法及应用 |
| CN109257128A (zh) * | 2018-11-01 | 2019-01-22 | 南京邮电大学 | 一种基于傅里叶级数拟合去噪的频谱信号识别方法及系统 |
| CN109632016A (zh) * | 2019-02-20 | 2019-04-16 | 湖北理工学院 | 岩土体硐室开挖与围岩应力、应变监测试验装置及其方法 |
Non-Patent Citations (3)
| Title |
|---|
| 刘伟 等: "高密度电阻率成像法与微动谱比法探测岩溶区塌陷的地质背景—以广东省高要市蛟塘镇塱下村塌陷区为例", 《地质与勘探》 * |
| 张振: "玉溪盆地三维速度结构建模及地脉动H/V谱模拟", 《中国优秀博硕士学位论文全文数据库(硕士) 基础科学辑》 * |
| 杨坤: "地铁盾构施工不良地质体微动探测技术研究及应用", 《轨道交通与地下工程》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112558149A (zh) * | 2020-11-12 | 2021-03-26 | 中铁第四勘察设计院集团有限公司 | 基于非均匀介质的微动探测方法、装置、设备及存储介质 |
| CN112558149B (zh) * | 2020-11-12 | 2023-04-25 | 中铁第四勘察设计院集团有限公司 | 基于非均匀介质的微动探测方法、装置、设备及存储介质 |
| CN114137609A (zh) * | 2021-11-09 | 2022-03-04 | 长江地球物理探测（武汉）有限公司 | 线性微动数据校正方法及装置 |
| CN114137609B (zh) * | 2021-11-09 | 2023-12-01 | 长江地球物理探测（武汉）有限公司 | 线性微动数据校正方法及装置 |
| CN118501939A (zh) * | 2024-05-20 | 2024-08-16 | 中国地质科学院 | 风化壳稀土资源评估方法、系统、设备、介质及程序产品 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110488349B (zh) | 2021-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Albarello et al. | The contribution of the ambient vibration prospecting in seismic microzoning: an example from the area damaged by the April 6, 2009 L’Aquila (Italy) earthquake | |
| Ivanov et al. | Delineating a shallow fault zone and dipping bedrock strata using multichannal analysis of surface waves with a land streamer | |
| CN108318918A (zh) | 基于微动频散曲线和h/v曲线的地下不良地质体无损探测方法及应用 | |
| Debeglia et al. | Karst investigations using microgravity and MASW; Application to Orléans, France | |
| CN110488349A (zh) | 基于微动三分量谱比vhsr的无损探测方法及应用 | |
| CN108152854A (zh) | 一种基于微动功率谱密度的无损探测方法及其应用 | |
| Pamuk et al. | 3D bedrock structure of Bornova plain and its surroundings (İzmir/western Turkey) | |
| Theodoulidis et al. | Basin effects on ground motion: the case of a high-resolution experiment in Cephalonia (Greece) | |
| Maghami et al. | Extracting the shear wave velocity structure of deep alluviums of “Qom” Basin (Iran) employing HVSR inversion of microtremor recordings | |
| Karastathis et al. | The contribution of geophysical techniques to site characterisation and liquefaction risk assessment: Case study of Nafplion City, Greece | |
| Hanson et al. | Advanced techniques in site characterization and mining hazard detection for the underground coal industry | |
| Macau et al. | Geophysical measurements for site effects characterisation in the urban area of Girona, Spain | |
| Zhang et al. | Multi-frequencies GPR measurements for delineating the shallow subsurface features of the Yushu strike slip fault | |
| Bertrand et al. | Site effects of the Roio basin, L’Aquila | |
| CN113050085A (zh) | 一种超前地质预报的方法 | |
| Di Giulio et al. | Using a vibratory source at Mt. Etna (Italy) to investigate the wavefield polarization at Pernicana Fault | |
| Durante et al. | A multidisciplinary approach to the seismic characterization of a mountain top (Monteluco, central Italy) | |
| Malik et al. | Ground-penetrating radar investigation along Pinjore Garden Fault: implication toward identification of shallow subsurface deformation along active fault, NW Himalaya | |
| Martino et al. | Integrated engineering-geological and numerical approach applied to the large Büyükçekmece (Turkey) landslide for evaluating earthquake-induced effects | |
| Montgomery et al. | Integrating surface-based geophysics into landslide investigations along highways | |
| Triantafyllidis et al. | Part I: theoretical site response estimation for microzoning purposes | |
| Schwardt et al. | A novel seismic full waveform inversion approach for assessing the internal structure of a medieval sea dike | |
| Cai et al. | Evaluate Geophysical Methods to Detect Underground Voids | |
| Civico et al. | Imaging the three-dimensional architecture of the Middle Aterno basin (2009 L’Aquila earthquake, Central Italy) using ground TDEM and seismic noise surveys: preliminary results | |
| Mondal et al. | Estimating Soil Thickness using the H/V Spectral Ratio Technique in CSIR-NGRI Campus, Hyderabad, India |
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 |