CN115685324B - Measuring device and method for non-uniform wave velocity field on rock surface - Google Patents

Measuring device and method for non-uniform wave velocity field on rock surface Download PDF

Info

Publication number
CN115685324B
CN115685324B CN202310010317.2A CN202310010317A CN115685324B CN 115685324 B CN115685324 B CN 115685324B CN 202310010317 A CN202310010317 A CN 202310010317A CN 115685324 B CN115685324 B CN 115685324B
Authority
CN
China
Prior art keywords
wave velocity
rock surface
velocity field
measuring
point
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.)
Active
Application number
CN202310010317.2A
Other languages
Chinese (zh)
Other versions
CN115685324A (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.)
Wuhan Institute of Rock and Soil Mechanics of CAS
Original Assignee
Wuhan Institute of Rock and Soil Mechanics of CAS
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 Wuhan Institute of Rock and Soil Mechanics of CAS filed Critical Wuhan Institute of Rock and Soil Mechanics of CAS
Priority to CN202310010317.2A priority Critical patent/CN115685324B/en
Publication of CN115685324A publication Critical patent/CN115685324A/en
Application granted granted Critical
Publication of CN115685324B publication Critical patent/CN115685324B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides a rock surface non-uniform wave velocity field measuring device and a measuring method thereof.A laser ranging radar is arranged above a precise electric steering engine and used for measuring the distance between the laser ranging radar and a marker post; the precise electric steering engine is used for realizing rotation, is connected with the rotatable slide rail and drives the rotatable slide rail to rotate together; a vibration receiving sensor is fixed below the precise electric steering engine; the lower part of the marker post is embedded into the electric slide block, and the lower end of the marker post is connected with a punch; the pressure sensor is arranged between the electric slide block and the punch; the electric sliding block is embedded in the rotatable sliding rail to form a sliding pair; the punching machine realizes rock surface impact; the vibration receiving sensor is used for receiving a vibration signal sent by the punching device when the punching device impacts the rock surface; the pressure sensor receives pressure pulses generated by the punch impacting the rock surface. The method is suitable for complex and variable geological condition environments and rock surfaces, and is simple and convenient to operate; the method has accurate knowledge of initial damage and structural structure of the rock mass surface and short measurement time.

Description

岩石表面非均匀波速场测量装置及其测量方法Measuring device and method for non-uniform wave velocity field on rock surface

技术领域technical field

本发明属于岩石工程稳定性监测领域,具体涉及一种岩石表面非均匀波速场测量装置及其测量方法。The invention belongs to the field of rock engineering stability monitoring, and in particular relates to a measuring device and a measuring method for a non-uniform wave velocity field on a rock surface.

背景技术Background technique

岩体表面或内部结构波速场的测定是岩石稳定性监测的重要内容。一方面,测定合理的波速场可为微震等监测提供基础,从而确定合理的到时,进而开展精准的定位和震源机理分析;另一方面,合理的波速场可为岩体表面或内部结构勘探提供参考,从而确定断层结构面等软弱地质结构的产状,进而确定岩体高应力集中区域。目前,在岩石表面波速场测定中,通常采用单一波速模型,即假设材料是均匀的,然而在现实中由于岩石在动载扰动等诸多外在因素影响作用下,使岩石产生一定的初始损伤或破坏,从而导致岩石表面波速场非均匀。The measurement of wave velocity field on the surface or internal structure of rock mass is an important content of rock stability monitoring. On the one hand, the determination of a reasonable wave velocity field can provide a basis for microseismic monitoring, so as to determine a reasonable arrival time, and then carry out precise positioning and source mechanism analysis; on the other hand, a reasonable wave velocity field can be used for rock mass surface or internal structure exploration. Provide a reference to determine the occurrence of weak geological structures such as fault structural planes, and then determine the high stress concentration area of the rock mass. At present, in the determination of the wave velocity field on the rock surface, a single wave velocity model is usually used, that is, the material is assumed to be uniform. damage, resulting in a non-uniform wave velocity field on the rock surface.

岩体表面非均匀波速场的获取方法主要包括走时或波形反演等,远远不能满足微震定位或岩体结构勘探的精度要求;而岩石表面波速场的测量仪器功能较为稀缺,现有技术存在很大的局限性且费时费力,难以适用于工程中复杂多变的地质条件环境。The acquisition methods of the non-uniform wave velocity field on the rock surface mainly include travel time or waveform inversion, etc., which are far from meeting the accuracy requirements of microseismic positioning or rock mass structure exploration; and the functions of measuring instruments for the rock surface wave velocity field are relatively scarce, and existing technologies exist It has great limitations and is time-consuming and labor-intensive, so it is difficult to apply to complex and changeable geological conditions in engineering.

因此,设计一种操作简单且能获取高质量岩石表面非均匀波速场的装置与方法具有重要的工程意义和应用价值。Therefore, it is of great engineering significance and application value to design a device and method that is simple to operate and can obtain high-quality rock surface non-uniform wave velocity field.

发明内容Contents of the invention

本发明的目的在于针对现有技术存在的上述问题,提供一种岩石表面非均匀波速场测量装置及其测量方法,本发明提高了岩石表面波速场测量的灵活性,且适用于复杂多变的地质条件环境。The purpose of the present invention is to solve the above-mentioned problems in the prior art, and provide a device for measuring the non-uniform wave velocity field on the rock surface and its measurement method. The present invention improves the flexibility of measuring the wave velocity field on the rock surface, and is applicable to geological conditions environment.

岩石表面非均匀波速场测量装置,包括:激光测距雷达;精密电动舵机;振动接收传感器;标志杆;电动滑块;压力传感器;冲压器;可旋转滑轨。Measuring device for non-uniform wave velocity field on rock surface, including: laser ranging radar; precision electric steering gear; vibration receiving sensor; marker pole; electric slider; pressure sensor; stamping device; rotatable slide rail.

所述的激光测距雷达安装在精密电动舵机上方,用于测量激光测距雷达与标志杆的距离;The laser ranging radar is installed above the precision electric steering gear, and is used to measure the distance between the laser ranging radar and the sign post;

所述的精密电动舵机用于实现360度高速旋转,连接可旋转滑轨,带动可旋转滑轨共同旋转;精密电动舵机下方固定有振动接收传感器;The precision electric steering gear is used to achieve 360-degree high-speed rotation, connected to the rotatable slide rail, and drives the rotatable slide rail to rotate together; a vibration receiving sensor is fixed under the precision electric steering gear;

所述的标志杆上有竖条状标志物,用于被激光测距雷达识别,标志杆下部嵌入到电动滑块,下端连接有冲压器;压力传感器安置在电动滑块和冲压器之间;There are vertical strip markers on the sign post for identification by the laser ranging radar. The lower part of the sign post is embedded in the electric slider, and the lower end is connected with a punch; the pressure sensor is placed between the electric slider and the punch;

所述的电动滑块嵌入在可旋转滑轨上形成滑动副,在可旋转滑轨上的滑动且位置可控;The electric slider is embedded on the rotatable slide rail to form a sliding pair, and the sliding on the rotatable slide rail is controllable in position;

所述的冲压器实现岩石表面撞击,撞击时带动标志杆和压力传感器共同移动;The stamping device realizes the impact on the rock surface, and drives the marker rod and the pressure sensor to move together when the impact occurs;

所述的振动接收传感器用于接收冲压器撞击岩石表面发出的振动信号;所述的压力传感器接收冲压器撞击岩石表面产生的压力脉冲。The vibration receiving sensor is used to receive the vibration signal sent by the punch hitting the rock surface; the pressure sensor receives the pressure pulse generated by the punch hitting the rock surface.

采用所述的岩石表面非均匀波速场装置,进行测量,方法的工作流程如下:Using the non-uniform wave velocity field device on the rock surface to measure, the working process of the method is as follows:

S1、仪器的安置;S1. Placement of instruments;

S2、复杂地质环境下岩石表面冲击;S2. Rock surface impact in complex geological environment;

S3、压力、振动与距离信号采集与处理;S3. Acquisition and processing of pressure, vibration and distance signals;

S4、单点波速的计算与获取;S4. Calculation and acquisition of single point wave velocity;

S5、圆形点阵中多点波速的获取;S5. Acquisition of multi-point wave velocities in the circular lattice;

S6、岩石表面非均匀波速场正反演与获取。S6. Forward inversion and acquisition of non-uniform wave velocity field on rock surface.

下面做进一步的详细说明。Further detailed description is given below.

S1、仪器的安置。所述的岩石表面非均匀波速场测量装置在使用时,将振动接收传感器固定于待监测岩石表面波速的区域中心位置,保证振动接收传感器与岩石表面贴合,即可开始岩石表面非均匀波速场的测定。S1. Placement of the instrument. When the device for measuring the non-uniform wave velocity field on the rock surface is in use, the vibration receiving sensor is fixed at the center of the area where the wave velocity on the rock surface is to be monitored, so as to ensure that the vibration receiving sensor is attached to the rock surface, and then the non-uniform wave velocity field on the rock surface can be started. determination.

S2、复杂地质环境下岩石表面冲击;冲压器在工作时具有可变的行程;S2. Rock surface impact in complex geological environment; the punch has a variable stroke during work;

S3、压力、振动与距离信号采集与处理;冲压器在冲击岩石表面后,压力传感器接 收到冲击压力脉冲信号,该信号通过信号转换公式获得AIC值,通过分析最小AIC值对应的 时间获得冲击信号的起始时间

Figure 490614DEST_PATH_IMAGE001
;压力传感器接收到冲击压力脉冲信号后,激光测距雷达 发出激光束阵列,遇到标志杆后发生反射后被激光测距雷达接收,即获得激光往返时间
Figure 980501DEST_PATH_IMAGE002
; 冲压器撞击岩石表面的振动波传播到振动接收传感器后,振动接收传感器接收到振动脉冲 信号,该信号通过信号转换公式获得AIC值,通过分析最小AIC值对应的时间获得振动信号 的到达时间
Figure 245260DEST_PATH_IMAGE003
。 S3. Acquisition and processing of pressure, vibration and distance signals; after the punch impacts the rock surface, the pressure sensor receives the impact pressure pulse signal, and the signal obtains the AIC value through the signal conversion formula, and the impact signal is obtained by analyzing the time corresponding to the minimum AIC value start time of
Figure 490614DEST_PATH_IMAGE001
; After the pressure sensor receives the impact pressure pulse signal, the laser ranging radar emits a laser beam array, which is reflected by the laser ranging radar after encountering the sign post, and the laser ranging time is obtained.
Figure 980501DEST_PATH_IMAGE002
; After the vibration wave that the punch hits the rock surface propagates to the vibration receiving sensor, the vibration receiving sensor receives the vibration pulse signal, and the signal obtains the AIC value through the signal conversion formula, and the arrival time of the vibration signal is obtained by analyzing the time corresponding to the minimum AIC value
Figure 245260DEST_PATH_IMAGE003
.

进一步地,所述的信号转换公式包括:Further, the signal conversion formula includes:

Figure 936136DEST_PATH_IMAGE004
Figure 936136DEST_PATH_IMAGE004

其中

Figure 879821DEST_PATH_IMAGE005
是当前时刻的AIC值,
Figure 45836DEST_PATH_IMAGE006
是时间采样点总数,
Figure 430681DEST_PATH_IMAGE007
是从初始时刻到当前时刻 的采样点,
Figure 151512DEST_PATH_IMAGE008
Figure 457860DEST_PATH_IMAGE009
是从初始时刻到当前时刻的采样点的平均值和方差,
Figure 555129DEST_PATH_IMAGE010
Figure 997742DEST_PATH_IMAGE011
是从当前时 刻到最终时刻的采样点的平均值和方差。 in
Figure 879821DEST_PATH_IMAGE005
is the AIC value at the current moment,
Figure 45836DEST_PATH_IMAGE006
is the total number of time sampling points,
Figure 430681DEST_PATH_IMAGE007
is the sampling point from the initial moment to the current moment,
Figure 151512DEST_PATH_IMAGE008
and
Figure 457860DEST_PATH_IMAGE009
is the mean and variance of the sampling points from the initial moment to the current moment,
Figure 555129DEST_PATH_IMAGE010
and
Figure 997742DEST_PATH_IMAGE011
is the mean and variance of the sampling points from the current moment to the final moment.

S4、单点波速的计算与获取。单点波速按照单点波速计算公式来获取。S4. Calculation and acquisition of single-point wave velocity. The single-point wave velocity is obtained according to the single-point wave velocity calculation formula.

进一步地,所述的单点波速计算公式包括:Further, the single-point wave velocity calculation formula includes:

Figure 889475DEST_PATH_IMAGE012
Figure 889475DEST_PATH_IMAGE012

其中

Figure 683119DEST_PATH_IMAGE013
是单点波速值,
Figure 584079DEST_PATH_IMAGE014
是激光在空气中的传播速度。 in
Figure 683119DEST_PATH_IMAGE013
is the single-point wave velocity value,
Figure 584079DEST_PATH_IMAGE014
is the propagation speed of the laser in air.

S5、圆形点阵中多点波速的获取;所述的圆形点阵是由以激光测距雷达为中心的多条监测轨迹组成的同心圆几何;S5. Acquisition of multi-point wave velocities in the circular lattice; the circular lattice is a concentric circle geometry composed of multiple monitoring tracks centered on the laser ranging radar;

当电动滑块位于可旋转滑轨距离激光测距雷达最远处监测轨迹,单个测点为最远处监测轨迹上某点;精密电动舵机带动可旋转滑轨沿360度高速旋转,此时重复步骤二到四,可获得最远处监测轨迹上所有测点的波速值。When the electric slider is located on the farthest monitoring track from the rotatable slide rail to the laser ranging radar, the single measuring point is a point on the farthest monitoring track; the precision electric steering gear drives the rotatable slide rail to rotate at a high speed of 360 degrees, at this time Repeat steps 2 to 4 to obtain the wave velocity values of all measuring points on the farthest monitoring track.

当最远处监测轨迹波速测量完成后,电动滑块带动标志杆、压力传感器和冲压器沿着可旋转滑轨由外向内移动到下一个监测轨迹,即可计算获得下一个监测轨迹上所有测点的波速值,不断循环重复该步骤直到最近处的监测轨迹,从而获得圆形点阵中多点波速。When the wave velocity measurement of the farthest monitoring track is completed, the electric slider drives the marker rod, the pressure sensor and the stamper to move from the outside to the inside along the rotatable slide rail to the next monitoring track, and then all the measurements on the next monitoring track can be calculated. The wave velocity value of the point, and repeat this step continuously until the nearest monitoring track, so as to obtain the wave velocity of multiple points in the circular lattice.

S6、岩石表面非均匀波速场正反演与获取;先确定圆形点阵区域范围并将其划分为栅格图,对每个栅格中的初始表面波速赋值为零,从而完成初始表面波速场赋值;S6. Forward inversion and acquisition of the non-uniform wave velocity field on the rock surface; first determine the range of the circular lattice area and divide it into grid images, and assign the initial surface wave velocity in each grid to zero, thereby completing the initial surface wave velocity field assignment;

根据实际测点波速开展波速场正反演;所述的实际测点波速为圆形点阵波速

Figure 677936DEST_PATH_IMAGE015
在 每个栅格中的线性插值;所述的波速场正反演包括波速场正演和波速场反演;所述的波速 场正演可通过程函方法开展;所述的波速场反演通过计算所述的波速场正演获得的波速场 和实际测点波速的差值,不断进行误差传播和参数优化,获得相对误差小于容许值时的反 演波速场,即最终波速场,即为本发明待测量的岩石表面非均匀波速场。 Carry out wave velocity field forward and inversion according to the wave velocity of the actual measuring point; the wave velocity of the actual measuring point is circular lattice wave velocity
Figure 677936DEST_PATH_IMAGE015
Linear interpolation in each grid; the wave velocity field forward and inversion includes wave velocity field forward modeling and wave velocity field inversion; the wave velocity field forward modeling can be carried out by the process function method; the wave velocity field inversion By calculating the difference between the wave velocity field obtained by the forward modeling of the wave velocity field and the wave velocity of the actual measuring point, the error propagation and parameter optimization are continuously carried out to obtain the inversion wave velocity field when the relative error is less than the allowable value, that is, the final wave velocity field, which is The non-uniform wave velocity field on the rock surface to be measured in the present invention.

本发明相对于现有技术,具有以下有益效果:Compared with the prior art, the present invention has the following beneficial effects:

1、该测量装置与方法可适用于复杂多变的地质条件环境和岩石表面,操作简便,灵活程度较高;1. The measuring device and method can be applied to complex and changeable geological conditions and rock surfaces, and are easy to operate and highly flexible;

3、通过转动滑轨和激光雷达可精确获取引发振动源的距离和方位,对了解岩体表面初始损伤和结构构造更加准确;3. The distance and orientation of the vibration source can be accurately obtained by rotating the slide rail and the laser radar, which is more accurate for understanding the initial damage and structural structure of the rock mass surface;

4、测量装置布设后滑轨可以360度全方位高速旋转,可以自动控制冲压器与接收传感器的距离,测量时间短。4. After the measurement device is installed, the slide rail can rotate 360 degrees in all directions at high speed, and can automatically control the distance between the punch and the receiving sensor, and the measurement time is short.

附图说明Description of drawings

以下结合附图和技术方案,进一步说明本发明的具体实施方式。The specific implementation manners of the present invention will be further described below in conjunction with the accompanying drawings and technical solutions.

图1为本发明装置的整体结构图;Fig. 1 is the overall structural diagram of device of the present invention;

图2为本发明方法的工作流程图;Fig. 2 is the workflow diagram of the inventive method;

图3a为本发明测量局部凹陷的岩石表面时冲压器经过较大的行程示意图;Fig. 3 a is a schematic diagram of a relatively large stroke of the punch when the present invention measures a partially depressed rock surface;

图3b为本发明测量局部凸出的岩石表面时冲压器经过较小的行程示意图;Fig. 3b is a schematic diagram of a smaller stroke of the punch when the present invention measures a partially protruding rock surface;

图4a为本发明获得激光往返时间

Figure 615937DEST_PATH_IMAGE016
示意图; Fig. 4a obtains laser round-trip time for the present invention
Figure 615937DEST_PATH_IMAGE016
schematic diagram;

图4b为本发明压力传感器接收到冲击压力脉冲信号;Fig. 4b shows the shock pressure pulse signal received by the pressure sensor of the present invention;

图4c为本发明振动接收传感器接收到振动脉冲信号;Fig. 4c is the vibration pulse signal received by the vibration receiving sensor of the present invention;

图5a为本发明获取圆形点阵中多点波速的原理示意图;Figure 5a is a schematic diagram of the principle of obtaining multi-point wave velocities in a circular lattice according to the present invention;

图5b为本发明获取C8监测轨迹中多点波速的原理示意图;Fig. 5b is a schematic diagram of the principle of obtaining multi-point wave velocity in the C8 monitoring track according to the present invention;

图6为本发明装置获取岩石表面非均匀波速场的原理示意图。Fig. 6 is a schematic diagram of the principle of obtaining a non-uniform wave velocity field on a rock surface by the device of the present invention.

具体实施方式Detailed ways

为了便于本领域普通技术人员理解和实施本发明,下面结合附图对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。In order to facilitate those skilled in the art to understand and implement the present invention, the present invention will be described in further detail below in conjunction with the accompanying drawings. It should be understood that the implementation examples described here are only used to illustrate and explain the present invention, and are not intended to limit the present invention.

本发明提供了一种岩石表面非均匀波速场测量装置,实现复杂地质环境下的岩石表面非均匀波速场测量。如图1所示,所述的岩石表面非均匀波速场测量装置,包括:激光测距雷达1;精密电动舵机2;振动接收传感器3;标志杆4;电动滑块5;压力传感器6;冲压器7;可旋转滑轨8。The invention provides a rock surface non-uniform wave velocity field measuring device, which realizes the measurement of the rock surface non-uniform wave velocity field in a complex geological environment. As shown in Figure 1, the non-uniform wave velocity field measuring device on the rock surface includes: a laser ranging radar 1; a precision electric steering gear 2; a vibration receiving sensor 3; a marker pole 4; an electric slider 5; a pressure sensor 6; punch 7; rotatable slide rail 8.

所述的激光测距雷达1位于精密电动舵机2上方,可测量激光测距雷达1与标志杆4的距离;The laser ranging radar 1 is located above the precision electric steering gear 2, and can measure the distance between the laser ranging radar 1 and the sign post 4;

所述的精密电动舵机2可实现360度高速旋转,其上方固定有激光测距雷达1,下方固定有振动接收传感器3,主体连接可旋转滑轨8,使用时可带动可旋转滑轨8共同旋转;The precision electric steering gear 2 can realize 360-degree high-speed rotation, a laser ranging radar 1 is fixed above it, a vibration receiving sensor 3 is fixed below it, and the main body is connected to a rotatable slide rail 8, which can drive the rotatable slide rail 8 during use common rotation;

所述的振动接收传感器3主要用于接收冲压器7撞击岩石表面发出的振动信号,其上方固定有精密电动舵机2;The vibration receiving sensor 3 is mainly used to receive the vibration signal sent by the punch 7 hitting the rock surface, and a precision electric steering gear 2 is fixed above it;

所述的标志杆4上有竖条状标志物,可被激光测距雷达1识别,其下部嵌入到电动滑块5中且下方连接有压力传感器6;There is a vertical strip-shaped marker on the sign post 4, which can be recognized by the laser ranging radar 1, and its lower part is embedded in the electric slider 5 and a pressure sensor 6 is connected below;

所述的电动滑块5嵌入在可旋转滑轨8上,可实现在可旋转滑轨8上的滑动且位置可控,其上方下方分别连接有标志杆4和压力传感器6;The electric slider 5 is embedded on the rotatable slide rail 8, and the sliding on the rotatable slide rail 8 can be realized and the position is controllable, and the marker rod 4 and the pressure sensor 6 are respectively connected above and below;

所述的压力传感器6可接收冲压器7撞击岩石表面产生的压力脉冲,其安置在电动滑块5和冲压器7之间;The pressure sensor 6 can receive the pressure pulse generated by the punch 7 hitting the rock surface, and it is arranged between the electric slider 5 and the punch 7;

所述的冲压器7可实现岩石表面撞击,撞击时可带动标志杆4和压力传感器6共同移动;The punch 7 can realize the impact on the rock surface, and can drive the marker rod 4 and the pressure sensor 6 to move together when the impact occurs;

所述的可旋转滑轨8与精密电动舵机2连接,可实现360度旋转,其上嵌入有电动滑块5,电动滑块5可在可旋转滑轨8上自由滑动。The rotatable slide rail 8 is connected with the precision electric steering gear 2 and can rotate 360 degrees. An electric slider 5 is embedded thereon, and the electric slider 5 can slide freely on the rotatable slide rail 8 .

如图2所示,所述的岩石表面非均匀波速场测量方法的工作流程如下:As shown in Figure 2, the workflow of the method for measuring the non-uniform wave velocity field on the rock surface is as follows:

S1、仪器的安置;S1. Placement of instruments;

S2、复杂地质环境下岩石表面冲击;S2. Rock surface impact in complex geological environment;

S3、压力、振动与距离信号采集与处理;S3. Acquisition and processing of pressure, vibration and distance signals;

S4、单点波速的计算与获取;S4. Calculation and acquisition of single point wave velocity;

S5、圆形点阵中多点波速的获取;S5. Acquisition of multi-point wave velocities in the circular lattice;

S6、岩石表面非均匀波速场正反演与获取。S6. Forward inversion and acquisition of non-uniform wave velocity field on rock surface.

下面结合工作流程和具体实施例附图对本发明的工作流程和工作原理做进一步的详细说明。The working flow and working principle of the present invention will be further described in detail below in conjunction with the working flow and the accompanying drawings of specific embodiments.

S1、仪器的安置。所述的岩石表面非均匀波速场测量装置在使用时可将振动接收传感器3固定于待监测岩石表面波速的区域大致中心位置,保证振动接收传感器3与岩石表面大致贴合,即可开始岩石表面非均匀波速场的测定。S1. Placement of the instrument. When the device for measuring the non-uniform wave velocity field on the rock surface is used, the vibration receiving sensor 3 can be fixed at the approximate center of the area where the wave velocity of the rock surface is to be monitored. Determination of inhomogeneous wave velocity fields.

S2、复杂地质环境下岩石表面冲击。所述的岩石表面非均匀波速场测量装置不受岩石表面粗糙起伏、局部地质结构或整体曲率影响,冲压器7在工作时具有可变的行程。如图3a所示,当遇到局部凹陷的岩石表面时,冲压器7在经过较大的行程后冲击岩石表面;如图3b所示,当遇到局部凸出的岩石表面时,冲压器7在经过较小的行程后冲击岩石表面。S2. Rock surface impact in complex geological environment. The device for measuring the non-uniform wave velocity field on the rock surface is not affected by rough undulations of the rock surface, local geological structure or overall curvature, and the punch 7 has a variable stroke during operation. As shown in Figure 3a, when encountering a partially concave rock surface, the punch 7 impacts the rock surface after a relatively large stroke; as shown in Figure 3b, when encountering a locally protruding rock surface, the punch 7 Hit the rock face after a small stroke.

S3、压力、振动与距离信号采集与处理。如图4a到图4c所示,冲压器7在冲击岩石表 面后,压力传感器6会接收到冲击压力脉冲信号,该信号可通过信号转换公式获得AIC值,通 过分析最小AIC值对应的时间可获得冲击信号的起始时间

Figure 755931DEST_PATH_IMAGE017
;压力传感器6接收到冲击压力 脉冲信号后,激光测距雷达1发出激光束阵列,遇到标志杆4后发生反射后被激光测距雷达1 接收,即可获得激光往返时间
Figure 335948DEST_PATH_IMAGE002
;冲压器7撞击岩石表面的振动波传播到振动接收传感器3 后,振动接收传感器3会接收到振动脉冲信号,该信号可通过信号转换公式获得AIC值,通过 分析最小AIC值对应的时间可获得振动信号的到达时间
Figure 552821DEST_PATH_IMAGE003
。 S3. Acquisition and processing of pressure, vibration and distance signals. As shown in Figures 4a to 4c, after the punch 7 impacts the rock surface, the pressure sensor 6 will receive the impact pressure pulse signal, which can be obtained through the signal conversion formula. The AIC value can be obtained by analyzing the time corresponding to the minimum AIC value The start time of the shock signal
Figure 755931DEST_PATH_IMAGE017
After the pressure sensor 6 receives the impact pressure pulse signal, the laser ranging radar 1 sends out a laser beam array, which is reflected by the laser ranging radar 1 after encountering the signpost 4, and the laser ranging time can be obtained
Figure 335948DEST_PATH_IMAGE002
; After the vibration wave that the punch 7 hits the rock surface propagates to the vibration receiving sensor 3, the vibration receiving sensor 3 will receive the vibration pulse signal, which can obtain the AIC value through the signal conversion formula, and can be obtained by analyzing the time corresponding to the minimum AIC value The arrival time of the vibration signal
Figure 552821DEST_PATH_IMAGE003
.

进一步地,所述的信号转换公式包括:Further, the signal conversion formula includes:

Figure 51936DEST_PATH_IMAGE018
Figure 51936DEST_PATH_IMAGE018

其中

Figure 289013DEST_PATH_IMAGE005
是当前时刻的AIC值,
Figure 531776DEST_PATH_IMAGE006
是时间采样点总数,
Figure 600226DEST_PATH_IMAGE019
是从初始时刻到当前时刻 的采样点,
Figure 880029DEST_PATH_IMAGE008
Figure 666719DEST_PATH_IMAGE009
是从初始时刻到当前时刻的采样点的平均值和方差,
Figure 182014DEST_PATH_IMAGE010
Figure 370550DEST_PATH_IMAGE011
是从当前时 刻到最终时刻的采样点的平均值和方差。 in
Figure 289013DEST_PATH_IMAGE005
is the AIC value at the current moment,
Figure 531776DEST_PATH_IMAGE006
is the total number of time sampling points,
Figure 600226DEST_PATH_IMAGE019
is the sampling point from the initial moment to the current moment,
Figure 880029DEST_PATH_IMAGE008
and
Figure 666719DEST_PATH_IMAGE009
is the mean and variance of the sampling points from the initial moment to the current moment,
Figure 182014DEST_PATH_IMAGE010
and
Figure 370550DEST_PATH_IMAGE011
is the mean and variance of the sampling points from the current moment to the final moment.

S4、单点波速的计算与获取。单点波速可按照单点波速计算公式来获取。S4. Calculation and acquisition of single-point wave velocity. The single-point wave velocity can be obtained according to the single-point wave velocity calculation formula.

进一步地,所述的单点波速计算公式包括:Further, the single-point wave velocity calculation formula includes:

Figure 86833DEST_PATH_IMAGE020
Figure 86833DEST_PATH_IMAGE020

其中

Figure 688716DEST_PATH_IMAGE021
是单点波速值,
Figure 883068DEST_PATH_IMAGE014
是激光在空气中的传播速度。 in
Figure 688716DEST_PATH_IMAGE021
is the single-point wave velocity value,
Figure 883068DEST_PATH_IMAGE014
is the propagation speed of the laser in air.

S5、圆形点阵中多点波速的获取。如图5a所示,所述的圆形点阵是由以激光测距雷达1为中心的多条监测轨迹组成的同心圆几何。S5. Acquisition of multi-point wave velocities in the circular lattice. As shown in FIG. 5 a , the circular lattice is geometrically concentric circles composed of multiple monitoring tracks centered on the laser ranging radar 1 .

如图5b所示,当电动滑块5位于可旋转滑轨8距离激光测距雷达1较远处时,此时位于C8监测轨迹,单个测点9为C8监测轨迹上某点。精密电动舵机2可带动可旋转滑轨8以较小的步进角沿360度高速旋转,此时重复步骤二到四,可获得C8监测轨迹上所有测点的波速值。As shown in Figure 5b, when the electric slider 5 is located at a distance from the rotatable slide rail 8 to the laser ranging radar 1, it is located on the C8 monitoring track at this time, and the single measuring point 9 is a certain point on the C8 monitoring track. The precision electric steering gear 2 can drive the rotatable slide rail 8 to rotate at a high speed of 360 degrees with a small step angle. At this time, repeat steps 2 to 4 to obtain the wave velocity values of all measuring points on the C8 monitoring track.

进一步地,当C8监测轨迹波速测量完成后,电动滑块5可带动标志杆4、压力传感器6和冲压器7沿着可旋转滑轨8由外向内移动到C7监测轨迹10,如图5a所示,即可计算获得C7监测轨迹10上所有测点的波速值,不断循环重复该步骤直到较内侧的C1监测轨迹,从而获得圆形点阵中多点波速。Further, after the wave velocity measurement of the C8 monitoring track is completed, the electric slider 5 can drive the marker rod 4, the pressure sensor 6 and the stamper 7 to move from the outside to the inside along the rotatable slide rail 8 to the C7 monitoring track 10, as shown in Figure 5a Shown, the wave velocity values of all the measuring points on the C7 monitoring track 10 can be calculated, and this step is repeated continuously until the inner C1 monitoring track, so as to obtain the wave velocities of multiple points in the circular lattice.

S6、岩石表面非均匀波速场正反演与获取。如图6所示,先确定圆形点阵区域范围并将其划分为一定边长的栅格图,对每个栅格中的初始表面波速赋值为零,从而完成初始表面波速场赋值;S6. Forward inversion and acquisition of non-uniform wave velocity field on rock surface. As shown in Figure 6, first determine the range of the circular lattice area and divide it into a grid map with a certain side length, and assign the initial surface wave velocity in each grid to zero, thereby completing the assignment of the initial surface wave velocity field;

根据实际测点波速开展波速场正反演;所述的实际测点波速为圆形点阵波速

Figure 254006DEST_PATH_IMAGE021
在每个栅格中的线性插值;所述的波速场正反演包括波速场正演和波速场反演;所述的波 速场正演可通过程函方法开展;所述的波速场反演通过计算所述的波速场正演获得的波速 场和实际测点波速的差值,不断进行误差传播和参数优化,获得相对误差小于较小容许值 时的反演波速场,即最终波速场;所述的最终波速场即为本发明待测量的岩石表面非均匀 波速场。Carry out wave velocity field forward and inversion according to the wave velocity of the actual measuring point; the wave velocity of the actual measuring point is circular lattice wave velocity
Figure 254006DEST_PATH_IMAGE021
Linear interpolation in each grid; the wave velocity field forward and inversion includes wave velocity field forward modeling and wave velocity field inversion; the wave velocity field forward modeling can be carried out by the process function method; the wave velocity field inversion By calculating the difference between the wave velocity field obtained by the forward modeling of the wave velocity field and the actual measuring point wave velocity, the error propagation and parameter optimization are continuously carried out to obtain the inversion wave velocity field when the relative error is less than a small allowable value, that is, the final wave velocity field; The final wave velocity field is the non-uniform wave velocity field on the rock surface to be measured in the present invention.

Claims (9)

1. Rock surface inhomogeneous wave velocity field measuring device, its characterized in that includes: a laser range radar (1); a precision electric steering engine (2); a vibration receiving sensor (3); a marking rod (4); an electric slider (5); a pressure sensor (6); a punch (7); a rotatable slide rail (8);
the laser ranging radar (1) is arranged above the precise electric steering engine (2) and is used for measuring the distance between the laser ranging radar (1) and the marker post (4);
the precise electric steering engine (2) is used for realizing rotation, is connected with the rotatable sliding rails (8) and drives the rotatable sliding rails (8) to rotate together; a vibration receiving sensor (3) is fixed below the precise electric steering engine (2);
the marking rod (4) is provided with a vertical strip-shaped marker which is used for being identified by the laser ranging radar (1), the lower part of the marking rod (4) is embedded into the electric sliding block (5), and the lower end of the marking rod is connected with a stamping device (7); the pressure sensor (6) is arranged between the electric slide block (5) and the punch (7);
the electric sliding block (5) is embedded in the rotatable sliding rail (8) to form a sliding pair, slides on the rotatable sliding rail (8) and has a controllable position;
the punching device (7) realizes rock surface impact and drives the marking rod (4) and the pressure sensor (6) to move together during impact;
the vibration receiving sensor (3) is used for receiving a vibration signal sent by the impact of the punch (7) on the rock surface; the pressure sensor (6) receives pressure pulses generated by the impact of the punch (7) on the rock surface.
2. The rock surface non-uniform wave velocity field measuring method is characterized in that the rock surface non-uniform wave velocity field measuring device of claim 1 is adopted, and the method comprises the following steps:
s1, arranging an instrument;
s2, rock surface impact in a complex geological environment;
s3, collecting and processing pressure, vibration and distance signals;
s4, calculating and obtaining the single-point wave velocity;
s5, acquiring the multi-point wave velocity in the circular dot matrix;
and S6, carrying out positive inversion and acquisition on the rock surface non-uniform wave velocity field.
3. The method for measuring the rock surface non-uniform wave velocity field according to the claim 2, characterized in that in the step S1, when the device for measuring the rock surface non-uniform wave velocity field is used, the vibration receiving sensor (3) is fixed at the center of the rock surface wave velocity region to be monitored, so that the vibration receiving sensor (3) is ensured to be attached to the rock surface, and the measurement of the rock surface non-uniform wave velocity field can be started.
4. The method for measuring the nonuniform wave velocity field of the rock surface according to claim 2, wherein in S2, the punch (7) has a variable stroke during operation.
5. The method for measuring the nonuniform wave velocity field of the rock surface according to claim 2, wherein the S3 is implemented by the following steps that after the punching device (7) impacts the rock surface, the pressure sensor (6) receives an impact pressure pulse signal, the signal obtains an AIC value through a signal conversion formula, and the starting time of the impact signal is obtained by analyzing the time corresponding to the minimum AIC value
Figure 879156DEST_PATH_IMAGE001
(ii) a After the pressure sensor (6) receives the impact pressure pulse signal, the laser ranging radar (1) sends out a laser beam array, and the laser ranging radar (1) receives the laser beam array after the laser ranging radar meets the marker post (4) and is reflected, so that the laser round-trip time is obtained
Figure 595439DEST_PATH_IMAGE002
(ii) a After the vibration wave of the impact rock surface of the punching device (7) is transmitted to the vibration receiving sensor (3), the vibration receiving sensor (3) receives a vibration pulse signal, the signal obtains an AIC value through a signal conversion formula, and the arrival time of the vibration signal is obtained by analyzing the time corresponding to the minimum AIC value
Figure 931742DEST_PATH_IMAGE003
6. The method for measuring the nonuniform wave velocity field of the rock surface according to claim 5, wherein the signal conversion formula comprises:
Figure 139477DEST_PATH_IMAGE004
wherein
Figure 510415DEST_PATH_IMAGE005
Is the AIC value at the current time instant,
Figure 256654DEST_PATH_IMAGE006
is the total number of time sample points,
Figure 955620DEST_PATH_IMAGE007
is the sample point from the initial time to the current time,
Figure 343876DEST_PATH_IMAGE008
and
Figure 428375DEST_PATH_IMAGE009
is the mean and variance of the sample points from the initial time to the current time,
Figure 611095DEST_PATH_IMAGE010
and
Figure 62936DEST_PATH_IMAGE011
is the mean and variance of the sample points from the current time to the final time.
7. The method for measuring the nonuniform wave velocity field on the rock surface according to claim 5, wherein in S4, the single-point wave velocity is obtained according to a single-point wave velocity calculation formula; the single-point wave velocity calculation formula comprises:
Figure 723725DEST_PATH_IMAGE012
wherein
Figure 803676DEST_PATH_IMAGE014
Is the value of the velocity of the single-point wave,
Figure 783396DEST_PATH_IMAGE016
is the propagation speed of the laser in air.
8. The method for measuring the nonuniform wave velocity field of the rock surface according to claim 2, wherein in S5, the circular lattice is a concentric circle geometry composed of a plurality of monitoring tracks with a laser ranging radar (1) as a center;
when the electric sliding block (5) is positioned on a monitoring track at the farthest position of the rotatable sliding rail (8) from the laser ranging radar (1), a single measuring point (9) is a certain point on the monitoring track at the farthest position; the precise electric steering engine (2) drives the rotatable slide rail (8) to rotate at a high speed along 360 degrees, and the second step to the fourth step are repeated at the moment, so that wave velocity values of all measuring points on the farthest monitoring track can be obtained;
after the wave velocity of the farthest monitoring track is measured, the electric slide block (5) drives the mark rod (4), the pressure sensor (6) and the punch (7) to move to the next monitoring track from outside to inside along the rotatable slide rail (8), the wave velocity values of all measuring points on the next monitoring track can be calculated and obtained, and the steps are repeated continuously and circularly until the monitoring track at the nearest position is reached, so that the wave velocity of multiple points in the circular dot matrix is obtained.
9. The method for measuring the rock surface nonuniform wave velocity field according to claim 2, wherein in S6, the rock surface nonuniform wave velocity field is inverted and acquired; firstly, determining a circular lattice area range, dividing the circular lattice area range into a grid map, and assigning the initial surface wave velocity in each grid to be zero so as to complete the initial surface wave velocity field assignment;
carrying out forward and backward deduction of a wave velocity field according to the wave velocity of the actual measuring point; the wave velocity of the actual measuring point is the wave velocity of a circular lattice
Figure 316008DEST_PATH_IMAGE017
Linear interpolation in each grid; the forward and backward evolution of the wave velocity field comprises forward evolution of the wave velocity field and inversion of the wave velocity field; the forward evolution of the wave velocity field is carried out by a path function method; and the wave velocity field inversion continuously performs error propagation and parameter optimization by calculating the difference value between the wave velocity field obtained by forward modeling of the wave velocity field and the wave velocity of the actual measuring point, so as to obtain the inverted wave velocity field when the relative error is smaller than an allowable value, namely the final wave velocity field, wherein the final wave velocity field is the rock surface non-uniform wave velocity field to be measured.
CN202310010317.2A 2023-01-05 2023-01-05 Measuring device and method for non-uniform wave velocity field on rock surface Active CN115685324B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310010317.2A CN115685324B (en) 2023-01-05 2023-01-05 Measuring device and method for non-uniform wave velocity field on rock surface

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310010317.2A CN115685324B (en) 2023-01-05 2023-01-05 Measuring device and method for non-uniform wave velocity field on rock surface

Publications (2)

Publication Number Publication Date
CN115685324A CN115685324A (en) 2023-02-03
CN115685324B true CN115685324B (en) 2023-03-10

Family

ID=85057417

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310010317.2A Active CN115685324B (en) 2023-01-05 2023-01-05 Measuring device and method for non-uniform wave velocity field on rock surface

Country Status (1)

Country Link
CN (1) CN115685324B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5166910A (en) * 1991-10-15 1992-11-24 Atlantic Richfield Company Method and apparatus for measuring the acoustic velocity
CN103713050A (en) * 2012-09-28 2014-04-09 中国石油化工股份有限公司 Method for measuring attenuation curve of seismic wave in rock by using laser receiving apparatus
CN104251883A (en) * 2013-06-28 2014-12-31 中国石油化工股份有限公司 Non-contact rock sound wave speed detection method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NO301736B1 (en) * 1994-09-14 1997-12-01 Nyfotek As Procedure for measuring sound speed and sample holder

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5166910A (en) * 1991-10-15 1992-11-24 Atlantic Richfield Company Method and apparatus for measuring the acoustic velocity
CN103713050A (en) * 2012-09-28 2014-04-09 中国石油化工股份有限公司 Method for measuring attenuation curve of seismic wave in rock by using laser receiving apparatus
CN104251883A (en) * 2013-06-28 2014-12-31 中国石油化工股份有限公司 Non-contact rock sound wave speed detection method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张一林 等.基于振动法的岩体纵波波速测试.2020,第46卷(第46期),第53-55页. *
靳郑伟 等.便携式岩石纵波波速自动测量仪研制.2021,第40卷(第40期),第170-176页. *

Also Published As

Publication number Publication date
CN115685324A (en) 2023-02-03

Similar Documents

Publication Publication Date Title
CN102053254B (en) Laser ultrasonic detection system and detection method thereof
CN102857850B (en) Near-field calibrating method for acoustic parameters of high-frequency ultrasonic emitter and array
CN110243320B (en) Method and device for non-contact measurement of tunnel lining crack depth
CN108732378B (en) An automated test method for acoustic Doppler velocity profilers
CN103017690A (en) Method for measuring straightness of super-long guide rail
CN108459083A (en) A kind of detecting system and its detection method for concrete dam
CN106546956A (en) A kind of ultra-short baseline receives the method for precisely marking of basic matrix primitive position
CN104251883A (en) Non-contact rock sound wave speed detection method
CN106767499A (en) A kind of preform flexibility measuring method
CN103713050B (en) A kind of laser receiver measures the method for the die-away curve of seismic event in rock
CN110927246A (en) Sleeve grouting compactness detection method and device based on ultrasonic CT
CN111221036A (en) A method and system for source location of target area containing unknown voids
CN111157065A (en) Acoustic time delay measuring method in ultrasonic signal transmission loop of gas ultrasonic flowmeter
CN103134441A (en) Large wind tunnel flexible spraying pipe laser tracking measurement method
CN115685324B (en) Measuring device and method for non-uniform wave velocity field on rock surface
CN111579646B (en) In-situ, non-destructive testing method for cracks in stone cultural relics
CN202230200U (en) Ultrasonic Positioning and Trajectory Monitoring System for Cabled Submersible
CN106556858B (en) A kind of ultrasonic signal excitation reception test macro
CN106772640B (en) Automatic precise positioning device and method for geological radar for tunnel disease detection
CN205593520U (en) Sectional device of rotation type laser survey tunnel two dimension
CN111707735B (en) Method for quantifying transverse cracks of fan spindle by using dual-mode diffraction waves
CN109470770A (en) A method for detecting continuous multiple oblique cracks in concrete structures
CN111413196B (en) Experimental method and system for measuring fracture change in rock joint surface shearing process
CN110346783B (en) Sound Ray Correction Method Based on Power Series Equivalent Ranging Model
CN114459360A (en) Precast concrete member surface roughness detection device and detection method

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