CN109407144A - 一种基于多波的单孔孤石立体探测方法 - Google Patents
一种基于多波的单孔孤石立体探测方法 Download PDFInfo
- Publication number
- CN109407144A CN109407144A CN201811476965.2A CN201811476965A CN109407144A CN 109407144 A CN109407144 A CN 109407144A CN 201811476965 A CN201811476965 A CN 201811476965A CN 109407144 A CN109407144 A CN 109407144A
- Authority
- CN
- China
- Prior art keywords
- wave
- boulder
- single hole
- signal
- waves
- 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 57
- 239000000523 sample Substances 0.000 title claims abstract description 31
- 238000001514 detection method Methods 0.000 claims abstract description 24
- 238000003384 imaging method Methods 0.000 claims abstract description 16
- 230000005012 migration Effects 0.000 claims abstract description 12
- 238000013508 migration Methods 0.000 claims abstract description 12
- 238000000926 separation method Methods 0.000 claims abstract description 9
- 230000001133 acceleration Effects 0.000 claims abstract description 6
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 230000005284 excitation Effects 0.000 claims description 24
- 238000005553 drilling Methods 0.000 claims description 18
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 239000011148 porous material Substances 0.000 claims description 6
- 238000003786 synthesis reaction Methods 0.000 claims description 4
- 239000007787 solid Substances 0.000 claims description 3
- 125000004122 cyclic group Chemical group 0.000 claims description 2
- 238000010892 electric spark Methods 0.000 claims description 2
- 239000002131 composite material Substances 0.000 claims 1
- 238000001308 synthesis method Methods 0.000 claims 1
- 238000010276 construction Methods 0.000 abstract description 10
- 238000011835 investigation Methods 0.000 abstract description 5
- 230000007812 deficiency Effects 0.000 abstract description 2
- 239000000203 mixture Substances 0.000 abstract description 2
- 230000005641 tunneling Effects 0.000 abstract description 2
- 239000002689 soil Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000002547 anomalous effect Effects 0.000 description 5
- 229910052704 radon Inorganic materials 0.000 description 4
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 4
- 230000001360 synchronised effect Effects 0.000 description 4
- 230000011514 reflex Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000010438 granite Substances 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000003325 tomography 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/02—Generating seismic energy
- G01V1/157—Generating seismic energy using spark discharges; using exploding wires
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/003—Seismic data acquisition in general, e.g. survey design
-
- 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/189—Combinations of different types of receiving elements
-
- 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
-
- 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
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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
- G01V1/364—Seismic filtering
-
- 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/12—Signal generation
- G01V2210/129—Source location
- G01V2210/1299—Subsurface, e.g. in borehole or below weathering layer or mud line
-
- 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/14—Signal detection
- G01V2210/142—Receiver location
- G01V2210/1429—Subsurface, e.g. in borehole or below weathering layer or mud line
-
- 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
- G01V2210/163—Cross-well
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
- G01V2210/24—Multi-trace filtering
- G01V2210/244—Radon transform
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
- G01V2210/512—Pre-stack
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
技术领域
本发明涉及工程地质勘查技术领域，具体涉及一种基于多波的单孔孤石立体探测方法。
背景技术
孤石是指花岗岩和凝灰岩的不均匀风化所致而形成的孤立岩石体，其形状各异，粒径从亚米到米级不等。孤石在我国南方地区普片存在，广东、福建更为集中。孤石的存在不仅对隧道盾构掘进造成困难，而且会对盾构机刀片产生严重磨损，同时对施工作业区域地层造成强扰动，影响施工安全。此外，孤石的存在还容易造成边坡失稳、建筑基础差异沉降及桩基沉桩断桩等工程地质问题。因此查明孤石的空间分布对工程建设十分重要。
根据孤石与周围介质的物性参数差异，孤石探测方法主要有电阻率CT(Computertomograpy，层析成像)，弹性波(或声波)CT，电磁波CT和微动探测法。由于孤石和土层的速度、密度差异很大，因此，相比其他方法，弹性波(声波)探测法具有较好的物性基础。运用这些方法，既可以从隧道里进行超前探测，也可以从地面上进行探测。在地面上探测具有不中断盾构施工的优点，但由于地面干扰强，浅地表结构复杂，在地面上基于钻孔的探测方法的分辨率更高。相比两孔或多孔的孔间弹性波探测，单孔弹性波探测不需要大量钻孔，具有探测成本低，施工准备快的特点。近几年在声波测井技术的基础上单孔弹性波探测孔旁地质异常体的方法得到了快速发展，主要有管波探测、共振波探测和声波远探测等方法。
管波探测是利用一种沿着钻孔传播的瑞雷面波对孔旁溶洞、软弱夹层等不良地质体有效探测的方法。该方法只能对井旁极小小范围地质情况进行探测，不能对远孔区域进行探测，且该方法无法确定异常方位。共振波探测主要依赖地质异常体受激诱发的共振信号进行目标体的定位，因共振信号产生的条件较为严苛，且其能量极弱，很难在钻孔中接收并提取出共振信号。因此，这种方法尚未见较好的实际应用效果的文献报道。声波远探测方法主要用于油气田测井中，受不能破坏油气井井壁的限制，其采用压电陶瓷换能器发射的声波功率不能太高，且频率较高。由于孤石主要分布于沉积土层中，属于慢速地层，高频信号衰减快，探测距离不大。此外，该方法运用的信号源是反射波，根据其探测原理可知，其无法对近孔区域异常体成像不佳。
专利号为CN1245637C(03/2006)、CN102565848B(11/2015)、CN103558637B(01/2016)、CN105589103B(02/2018)、CN208000384U(10/2018)、CN105604557B(03/2018)、CN106324683A(01/2017)、CN105804763B(03/2017)、CN103790594B(02/2016)、CN108693249A(10/2018)、CN107065014A(08/2017)的中国专利。
李学文，等.管波及其工程应用.《物探与化探》，2005，第29卷(第5期)，第463-466页.
张玉池，等.管波探测方法在岩土工程中的应用.《工程勘察》，2010，增刊(第s1期)，第640-645页.
唐晓明，等.偶极横波远探测测井技术进展及其应用.《测井技术》，2013，第37卷(第4期)，第333-339页.
李国英，等.方位远探测反射声波成像测井仪器.《测井技术》，2018，第42卷(第2期)，第221-226页.
上述专利或文献集中介绍了孤石探测的方法种类，以及单孔弹性波(管波、共振波、反射波)地质异常体探测方法特点。其现有单孔弹性波探测方法存在显著不足之处在于：单一方法只重点采集和利用单一类型波信号，使得成像范围受限。本发明介绍了一种基于单个钻孔，选用中低频激发源，多方位接收，同步采集全波信号，并利用管波和反射波联合解释，实现基于多波单孔孤石立体全域准确探测的方法。
发明内容
本发明的目的是针对上述不足，提出了一种基于单个钻孔，选用中低频脉冲振动波激发源，多方位接收，同步采集全波信号，并利用管波和反射波联合解释，实现基于多波单孔孤石立体全域准确探测的方法。
本发明具体采用如下技术方案：
一种基于多波的单孔孤石立体探测方法，具体包括：
a.选用中低频脉冲波激发仪作为震源，八分量加速度传感器接和姿态传感器组成接收探头，以及信号记录仪为主组成单孔多波信号采集***；
b.固定激发和接收探头距离，在同一钻孔中，由深至浅按固定步距提升，完成不同深度的全波信号数据采集和探头姿态信息记录；
c.根据姿态传感器提供的不同采集深度对应的接收探头方位，对获取的八分量信息进行旋转并合成得到固定坐标系下的十六分量信息；
d.选排组合得到固定坐标系下不同方位的深度-时间波形剖面，对管波和反射波波场分离；
e.对分离得到的反射波剖面，采用绕射速度分析方法获取不同方位、不同深度的纵波速度；
f.利用分析所得的速度场信息，对反射波剖面深度-时间波形记录实施偏移成像；
g.分别基于管波信号和反射波成像结果联合解释得到以钻孔为中心的圆柱体空间中，孤石的尺度及空间展布信息。
进一步地，所述单孔多波信号采集***中的脉冲波激发仪为激发频率范围在几百赫兹到几千赫兹、能量大于1万焦耳的电火花类型震源，八分量加速度传感器为环状45度等方位间隔分布，姿态传感器能测出探头方位角、俯仰角和横滚角，信号记录仪采样率最高采样率可达0.1ms。
进一步地，所述固定激发和接收探头距离取值范围0.5m-2m，固定步距取值范围0.25m-1m。
进一步地，所述全波信号至少包含有管波和反射波信息。
进一步地，所述的旋转合成为校正探头的俯仰角和横滚角，并采用相控阵信号合成方法合成固定坐标系下指定方位的信号。
进一步地，所述固定坐标系为地球正北向为0度，按逆时针递增，十六分量信息为在固定坐标系下方位间隔为22.5度情况下的分量信息。
进一步地，所述的联合解释是利用管波信息解释近孔区域孤石分布情况，同时利用反射波成像结果解释远孔区域孤石分布情况，综合二者解释成果实现近孔和远孔全域立体孤石探测。
优选地，所述管波和反射波波场分离采用拉动变换。
优选地，所述反射波偏移成像采用绕射叠前偏移成像。
本发明所述的一种基于多波的单孔孤石立体探测方法具有的有效果是：
(1)选用中低频脉冲波激发仪作为震源，八分量加速度传感器接和姿态传感器组成接收探头，以及信号记录仪为主组成单孔多波信号采集***，能够一次激发、多方位多类型波信号同步接收，适合低速土层中的孤石探测的情况，具有远距离多方位探测的能力。
(2)能够基于管波信号和反射波成像结果联合解释，解决了只利用单一信号解释存在探测盲区的缺点，可实现以钻孔为中心的圆柱体空间的全域孤石探测，获得其探测区孤石尺度及空间展布信息。
(3)该方法数据采集装置便携、施工简单、工期短、成本低、精度高，能够为隧道盾构施工及其它建筑工程建设中排除孤石影响，提供准确的地质信息。
附图说明
图1本方明方法所提供的单孔立体多波探测孤石方法流程图；
图2本发明方法所提供的单孔立体多波探测孤石的数据采集示意图；
图3实施例一中单一孤石空间分布方位模型图；
图4实施例一模型模拟得到的孤石深度不同方位深度-时间波形剖面图；
图5实施例二中孤石群空间分布、管波和反射波全波剖面及联合解释成果图。
图中：1-地面信号记录仪；2-地面震源控制器；3-激发电缆；4-信号电缆；5-井口滑轮；6-；7-接收探头；8-姿态传感器；9-信号传感器；10-震源发射头；11-钻孔；12-孤石；13-反射波；14-管波。
具体实施方式
下面结合附图和实施例对本发明做进一步的详细说明，但不作为对本发明进行任何限制的依据。
实施例一：
在实施例一中，探测对象为如图3所示的一个单一孤石，且远离钻孔的地质模型。孤石位于固定坐标系下方位角约为90度，埋深10m，距孔15m的位置，大小约为直径2m。采用高阶有限差分弹性波动方程正演模拟算法进行模拟，按图1所示本发明方法所提供的单孔立体多波探测孤石方法流程，具体包括：
a.按图2所示本发明方法所提供的单孔立体多波探测孤石的数据采集示意图，将震源头[10]和接收探头[7]，通过井口滑轮[5]输送至同一钻孔[11]中，部分激发电缆[3]和信号电缆[4]一并入井。控制地面震源控制器[2]，发出指令使震源发射头[10]放电激发出脉冲声波。利用接收信号传感器[9]同步接收来自不同方位反射纵波[13]信息、以及沿着井孔传播的管波[14]信号，通过信号电缆[4]，由地面信号记录仪[1]进行数据采集、模数转换并存储。同时，信号记录仪也将姿态传感器[8]记录的井中探头方位角、俯仰角和横滚角等信息记录存储。
b.按图2所示本发明方法所提供的单孔立体多波探测孤石的数据采集示意图，固定激发和接收探头距离1m，在同一钻孔中，由深(17m)至浅(2m)按固定步距(1m)提升，完成不同深度的全波信号数据采集和探头姿态信息记录。
c.根据姿态传感器提供的不同采集深度对应的接收探头方位，对获取的八分量信息进行旋转并合成得到固定坐标系下的十六分量信息。根据姿态传感器提供的不同采集深度对应的接收探头方位，对获取的八分量信息进行旋转，合成得到固定坐标系下的十六方位信息。对于某一深度h记录的数据，根据姿态传感记录的探头俯仰角θ，横滚角φ，对每一分量数据(i是分量序号，为传感器定义的0度方位与地球正北向之间的夹角，这里定义为传感器的方位角)都按下式做倾角校正：
根据姿态传感记录的接收传感器方位角采用相控阵接收波形合成法(该方法为本领域成熟方法)合成十六分量波形信号，因此本发明方位分辨率可达22.5度，图4所示为其中八个分量的波形剖面(方位角分别为0，45，90，135，180，225，270，315)。通过改用更高分量接收传感器，方位分辨率可进一步提高。按照上述方法，对所有深度测量所得的数据，都进行倾角校正和相控阵波形合成。值得注意的是：在实施例一中，所有数据的方位角、俯仰角和横滚角均为0，且八方位信号是由模拟所得水平XY分量信号合成的。
图4所示为不同方位深度-时间波形剖面中，孤石的反射波呈绕射双曲线形态，其中能量最强的绕射波位于90度方位角，绕射波顶点与孤石埋深一致。孤石偏离钻孔的距离为r＝vt0/2，v为土层速度，t0为绕射点自激自收旅行时。因此可以通过绕射波的能量和形态特征直接解释孤石的空间位置。也可进一步做偏移成像，进行孤石空间分布精确定位。图4所示剖面中无管波存在，其原因是由于钻孔旁无地质异常体存在，且模拟过程中未考虑地面和钻孔底部的管波反射问题。
d.对步骤c所合成的数据进行选排，组合得到不同方位的深度-时间波形剖面，使用拉东变换进行管波和反射纵波波场分离。拉东变换将时间-深度域内的波场沿着不同斜率、不同截距进行数据叠加，得到斜率-截距(τ-p)域内的波场，对于具有不同同相轴的波场经过变换以不同形式分布在τ-p域内不同位置对τ-p域内分离的波场进行拾取，再反变换即可得到分离的波场，实现波场分离。本发明中输入信号用表示，其中为方位角，z为深度，t为时间，拉东变换结果为：
变换到τ-p域中，其中τ为截距，p为斜率。根据初至波、管波和反射波在τ-p域中的分布特征，分别确定管波分离滤波器Ft(τ，p)和反射波分离滤波器Fr(τ，p)，并进行如下反变换：
从而实现管波和反射波提取。
e.对步骤d分离得到的反射波剖面，采用绕射速度分析方法获取不同方位不同深度介质的纵波速度。由于孤石尺度及形态特征，其反射信息主要以绕射被形式存在。这里假设主要以土层为主，其纵波速度为常数，绕射双曲线方程为：
其中z为激发接收中点与绕射点的垂直高度差，d为激发接收间距，v为土层速度，t0为绕射点自激自收旅行时，tz为激发接收中点偏离绕射点垂直高度差为z时的旅行时。按照上述公式(5)，即可进行基于绕射双曲线方程的速度扫描分析。
f.利用步骤e分析所得的速度场信息，对各方位二维深度-时间反射波记录实施偏移处理。首先把满足方程(5)的曲线上的地震能量(即采样点振幅)全部送到对应的绕射双曲线的顶点上去叠加，其顶点对应时间为
在此基础上，需要把叠加能量输送到零激发接收间距绕射双曲线的顶点上，从而实现绕射叠前偏移成像。因为非零激发接收间距绕射波顶点和零激发接收间距绕射双曲线顶点在深度上相同，上述叠加能量的只需在时间上做时移，时移计算公式如下：
按照上述步骤，对深度-时间剖面所有样点全部实施偏移，将偏移后的数据叠加即完成一个二维剖面的绕射叠前偏移处理。对所有方位反射波数据按此方法处理，孤石反射能量即可归位，恢复出其真实的空间形态。也可按照绕射双曲面做三维偏移，可以获得更高空间分辨率。
实施例二：
在实施例二中，探测对象为如图5(a)所示左侧孤石群模型，其中含2个孤石，孤石1靠近钻孔，孤石2远离钻孔，分布在同一方位不同深度的钻孔同一侧。孤石1位于固定坐标系下方位角为90度，埋深17m，距孔0.2m的位置，直径约2m。孤石2位于固定坐标系下方位角为90度，埋深23m，距孔12m的位置，直径约2m。采用高阶有限差分弹性波动方程正演模拟算法进行模拟，按图1所示本发明方法所提供的单孔立体多波探测孤石方法流程，具体包括：
a.按图2所示本发明方法所提供的单孔立体多波探测孤石的数据采集示意图，布设数据观测方式(与实施例一相同，不再赘述)。
b.按图2所示本发明方法所提供的单孔立体多波探测孤石的数据采集示意图，固定激发和接收探头距离1m，在同一钻孔中，由深(40m)至浅(1m)按固定步距(1m)提升，完成不同深度的全波信号数据采集和探头姿态信息记录。
c.根据姿态传感器提供的不同采集深度对应的接收探头方位，对获取的八分量信息进行旋转并合成得到固定坐标系下的十六分量信息(具体实施措施与实施例一相同)。图5(b)所示为合成方位角为90度所对应深度-时间全波波形剖面。图中含清晰可见的初至波(直达波)、反射波和管波。管波的特征表现为在孤石1界面上的反射，因此管波从0时刻就开始发育，且孤石上下均有反射现象，在波形剖面上呈V字形分布。V字形交叉处对应孤石1的空间位置。除此之外，管波还表现出频散特性。远离钻孔的孤石2表现出明显的反射现象，呈反射双曲线，无频散特性。反射双曲线的顶点对应深度为孤石2的埋深。孤石偏离钻孔的距离为r＝vt0/2，v为土层速度，t0为绕射点自激自收旅行时。
d.对步骤c所合成的数据进行选排，组合得到不同方位的深度-时间波形剖面，使用拉东变换进行管波和反射纵波波场分离(具体实施方案与实施例一相同，不再赘述)。
e.对步骤d分离得到的反射波剖面，采用绕射速度分析方法获取不同方位不同深度介质的纵波速度(具体实施方案与实施例一相同，不再赘述)。
f.利用步骤e分析所得的速度场信息，对反射波剖面深度-时间波形记录实施绕射叠前偏移成像。
g.分别基于管波信号和反射波成像结果联合解释得到以钻孔为中心的圆柱体空间中，不同孤石的尺度及空间展布信息。图5(c)所示为方位角为90度的综合成像结果，孤石群的分布特征与模型完全吻合。
由实施例一和实施例二可以看出，利用本发明方法可以综合利用管波和反射波，进行联合解释成像，既可以利用管波对近孔区域探测，也可同时利用反射波对远孔区域探测，避免了只利用单一信号解释存在探测盲区的缺点。
当然，上述说明并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换，也应属于本发明的保护范围。
Claims (7)
1.一种基于多波的单孔孤石立体探测方法，其特征在于，具体包括：
a.选用中低频脉冲波激发仪作为震源，八分量加速度传感器接和姿态传感器组成接收探头，以及信号记录仪为主组成单孔多波信号采集***；
b.固定激发和接收探头距离，在同一钻孔中，由深至浅按固定步距提升，完成不同深度的全波信号数据采集和探头姿态信息记录；
c.根据姿态传感器提供的不同采集深度对应的接收探头方位，对获取的八分量信息进行旋转并合成得到固定坐标系下的十六分量信息；
d.选排组合得到固定坐标系下不同方位的深度-时间波形剖面，对管波和反射波波场分离；
e.对分离得到的反射波剖面，采用绕射速度分析方法获取不同方位、不同深度的纵波速度；
f.利用分析所得的速度场信息，对反射波剖面深度-时间波形记录实施偏移成像；
g.分别基于管波信号和反射波成像结果联合解释得到以钻孔为中心的圆柱体空间中，孤石的尺度及空间展布信息。
2.如权利要求1所述的一种基于多波的单孔孤石立体探测方法，其特征在于，所述单孔多波信号采集***中的脉冲波激发仪为激发频率范围在几百赫兹到几千赫兹、能量大于1万焦耳的电火花类型震源，八分量加速度传感器为环状45度等方位间隔分布，姿态传感器能测出探头方位角、俯仰角和横滚角，信号记录仪采样率最高采样率可达0.1ms。
3.如权利要求1所述的一种基于多波的单孔孤石立体探测方法，其特征在于，所述固定激发和接收探头距离取值范围0.5m-2m，固定步距取值范围0.25m-1m。
4.如权利要求1所述的一种基于多波的单孔孤石立体探测方法，其特征在于，所述全波信号至少包含有管波和反射波信息。
5.如权利要求1所述的一种基于多波的单孔孤石立体探测方法，其特征在于，所述的旋转合成为校正探头的俯仰角和横滚角，并采用相控阵信号合成方法合成固定坐标系下指定方位的信号。
6.如权利要求1或5所述的一种基于多波的单孔孤石立体探测方法，其特征在于，所述固定坐标系为地球正北向为0度，按逆时针递增，十六分量信息为在固定坐标系下方位间隔为22.5度情况下的合成信息。
7.如权利要求1所述的一种基于多波的单孔孤石立体探测方法，其特征在于，所述的联合解释是利用管波信息解释近孔区域孤石分布情况，同时利用反射波成像结果解释远孔区域孤石分布情况，综合二者解释成果实现近孔和远孔全域立体孤石探测。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811476965.2A CN109407144B (zh) | 2018-12-05 | 2018-12-05 | 一种基于多波的单孔孤石立体探测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811476965.2A CN109407144B (zh) | 2018-12-05 | 2018-12-05 | 一种基于多波的单孔孤石立体探测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109407144A true CN109407144A (zh) | 2019-03-01 |
| CN109407144B CN109407144B (zh) | 2020-02-18 |
Family
ID=65457224
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811476965.2A Expired - Fee Related CN109407144B (zh) | 2018-12-05 | 2018-12-05 | 一种基于多波的单孔孤石立体探测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109407144B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109917450A (zh) * | 2019-04-10 | 2019-06-21 | 长江水利委员会长江科学院 | 一种全断面网络化布设的土石坝位移监测***及方法 |
| CN110670569A (zh) * | 2019-10-17 | 2020-01-10 | 华侨大学 | 一种在冲孔桩成孔过程中对地层情况进行勘察的装置及方法 |
| CN110988980A (zh) * | 2019-12-23 | 2020-04-10 | 山东大学 | 一种盾构机搭载的相控阵声波超前地质探测***及方法 |
| CN113279435A (zh) * | 2021-01-13 | 2021-08-20 | 浙江省建设工程质量检验站有限公司 | 基于旁孔绕射波分析的基桩质量无损检测方法 |
| CN113917561A (zh) * | 2021-09-17 | 2022-01-11 | 东华理工大学 | 一种矿区三维地质建模中目标地质体的确定方法及*** |
| CN115822579A (zh) * | 2022-12-13 | 2023-03-21 | 北京华晖探测科技股份有限公司 | 一种分方位孔中不良地质体探测技术 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH10121452A (ja) * | 1996-10-17 | 1998-05-12 | Raito Kogyo Co Ltd | 地盤改良範囲の判定方法 |
| CN1544934A (zh) * | 2003-11-25 | 2004-11-10 | 饶其荣 | 管波探测法 |
| CN102692650A (zh) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | 一种具有假频压制功能的井筒波分离方法 |
| CN104133247A (zh) * | 2014-07-31 | 2014-11-05 | 中国石油天然气集团公司 | 垂直地震剖面数据中套管波的压制方法及装置 |
| CN105604557A (zh) * | 2015-12-30 | 2016-05-25 | 福建工程学院 | 基于地震ct的盾构施工孤石探测方法 |
| CN208000384U (zh) * | 2018-02-24 | 2018-10-23 | 广州市城市建设工程监理公司 | 一种孤石探测装置 |
-
2018
- 2018-12-05 CN CN201811476965.2A patent/CN109407144B/zh not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH10121452A (ja) * | 1996-10-17 | 1998-05-12 | Raito Kogyo Co Ltd | 地盤改良範囲の判定方法 |
| CN1544934A (zh) * | 2003-11-25 | 2004-11-10 | 饶其荣 | 管波探测法 |
| CN102692650A (zh) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | 一种具有假频压制功能的井筒波分离方法 |
| CN104133247A (zh) * | 2014-07-31 | 2014-11-05 | 中国石油天然气集团公司 | 垂直地震剖面数据中套管波的压制方法及装置 |
| CN105604557A (zh) * | 2015-12-30 | 2016-05-25 | 福建工程学院 | 基于地震ct的盾构施工孤石探测方法 |
| CN208000384U (zh) * | 2018-02-24 | 2018-10-23 | 广州市城市建设工程监理公司 | 一种孤石探测装置 |
Non-Patent Citations (1)
| Title |
|---|
| 刘宏岳: "地震反射波CDP叠加技术在海域花岗岩孤石探测中的应用", 《工程地球物理学报》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109917450A (zh) * | 2019-04-10 | 2019-06-21 | 长江水利委员会长江科学院 | 一种全断面网络化布设的土石坝位移监测***及方法 |
| CN109917450B (zh) * | 2019-04-10 | 2024-02-06 | 长江水利委员会长江科学院 | 一种全断面网络化布设的土石坝位移监测***及方法 |
| CN110670569A (zh) * | 2019-10-17 | 2020-01-10 | 华侨大学 | 一种在冲孔桩成孔过程中对地层情况进行勘察的装置及方法 |
| CN110670569B (zh) * | 2019-10-17 | 2024-06-07 | 华侨大学 | 一种在冲孔桩成孔过程中对地层情况进行勘察的装置及方法 |
| CN110988980A (zh) * | 2019-12-23 | 2020-04-10 | 山东大学 | 一种盾构机搭载的相控阵声波超前地质探测***及方法 |
| WO2021129381A1 (zh) * | 2019-12-23 | 2021-07-01 | 山东大学 | 一种盾构机搭载的相控阵声波超前地质探测***及方法 |
| US11640007B2 (en) | 2019-12-23 | 2023-05-02 | Shandong University | System and method for phased array sound wave advanced geological exploration for shield tunneling machine |
| CN113279435A (zh) * | 2021-01-13 | 2021-08-20 | 浙江省建设工程质量检验站有限公司 | 基于旁孔绕射波分析的基桩质量无损检测方法 |
| CN113917561A (zh) * | 2021-09-17 | 2022-01-11 | 东华理工大学 | 一种矿区三维地质建模中目标地质体的确定方法及*** |
| CN113917561B (zh) * | 2021-09-17 | 2023-06-20 | 东华理工大学 | 一种矿区三维地质建模中目标地质体的确定方法及*** |
| CN115822579A (zh) * | 2022-12-13 | 2023-03-21 | 北京华晖探测科技股份有限公司 | 一种分方位孔中不良地质体探测技术 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109407144B (zh) | 2020-02-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109407144A (zh) | 一种基于多波的单孔孤石立体探测方法 | |
| US7245560B2 (en) | Acoustic source for infrasonic electromagnetic wave exploration using induced electrokinetic effect | |
| US10073184B2 (en) | Sensor system of buried seismic array | |
| US20040044479A1 (en) | Method and apparatus for interferometry, spectral analysis, and three-dimensional holographic imaging of hydrocarbon accumulations and buried objects | |
| CA2905892C (en) | Sensor system of buried seismic array | |
| CN102866417A (zh) | 一种地下溶洞地震跨孔ct探测及层析成像装置及方法 | |
| CN112068206B (zh) | 一种基于井孔与弹性波相互作用理论的邻井探测方法 | |
| CN202837558U (zh) | 一种地下溶洞地震跨孔ct探测及层析成像装置 | |
| CN109001839A (zh) | 一种城市车载综合地球物理探测***及工作方法 | |
| CN109782360A (zh) | 一种低阻覆盖区深部砂体探测方法 | |
| US6476608B1 (en) | Combining seismic waves with seismoelectrics to perform prospecting and measurements | |
| CN110187394A (zh) | 双场源电磁测深法获取地层电阻率各向异性的方法及装置 | |
| CN101100940B (zh) | 一种阵列化声信号检测***及其工程应用 | |
| CN105093314B (zh) | 一种测定微地震震源的方法 | |
| AU2022200417B2 (en) | A method for seismic frequency resonance exploration technology | |
| CN207829030U (zh) | 桩基浇注前桩底隐患探查装置 | |
| US9551803B2 (en) | Geological medium exploration | |
| CN113777653B (zh) | 一种合成孔径式浅层地震剖面勘探方法和*** | |
| CN114779346A (zh) | 一种基于正交偶极四分量测量精准探测井外地质构造的声波测井方法及装置 | |
| CN109375251B (zh) | 利用城市既有地下空间与地表的探测方法及*** | |
| Juhlin et al. | High‐resolution reflection seismics applied to detection of groundwater resources in glacial deposits, Sweden | |
| Wang et al. | Retrieving drill bit seismic signals using surface seismometers | |
| CN110780345A (zh) | 一种隧道超前地震探测地震数据的三维速度分析方法 | |
| Cheng et al. | 3D trajectory inversion of an adjacent well using scattered P-wave | |
| CN108732625A (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200218 Termination date: 20201205 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |