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 PDFInfo
- 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
Links
- 239000011435 rock Substances 0.000 title claims abstract description 87
- 238000000034 method Methods 0.000 title claims abstract description 27
- 239000003550 marker Substances 0.000 claims abstract description 14
- 238000005259 measurement Methods 0.000 claims abstract description 9
- 238000004080 punching Methods 0.000 claims abstract 6
- 238000012544 monitoring process Methods 0.000 claims description 28
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 5
- 238000005457 optimization Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims 2
- 230000003116 impacting effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 8
- 238000005070 sampling Methods 0.000 description 8
- 230000035939 shock Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000994 depressogenic effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
技术领域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值对应的 时间获得冲击信号的起始时间;压力传感器接收到冲击压力脉冲信号后,激光测距雷达 发出激光束阵列,遇到标志杆后发生反射后被激光测距雷达接收,即获得激光往返时间; 冲压器撞击岩石表面的振动波传播到振动接收传感器后,振动接收传感器接收到振动脉冲 信号,该信号通过信号转换公式获得AIC值,通过分析最小AIC值对应的时间获得振动信号 的到达时间。 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 ; 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. ; 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 .
进一步地,所述的信号转换公式包括:Further, the signal conversion formula includes:
其中是当前时刻的AIC值,是时间采样点总数,是从初始时刻到当前时刻 的采样点,和是从初始时刻到当前时刻的采样点的平均值和方差,和是从当前时 刻到最终时刻的采样点的平均值和方差。 in is the AIC value at the current moment, is the total number of time sampling points, is the sampling point from the initial moment to the current moment, and is the mean and variance of the sampling points from the initial moment to the current moment, and 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:
其中是单点波速值,是激光在空气中的传播速度。 in is the single-point wave velocity value, 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
当最远处监测轨迹波速测量完成后,电动滑块带动标志杆、压力传感器和冲压器沿着可旋转滑轨由外向内移动到下一个监测轨迹,即可计算获得下一个监测轨迹上所有测点的波速值,不断循环重复该步骤直到最近处的监测轨迹,从而获得圆形点阵中多点波速。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;
根据实际测点波速开展波速场正反演;所述的实际测点波速为圆形点阵波速在 每个栅格中的线性插值;所述的波速场正反演包括波速场正演和波速场反演;所述的波速 场正演可通过程函方法开展;所述的波速场反演通过计算所述的波速场正演获得的波速场 和实际测点波速的差值,不断进行误差传播和参数优化,获得相对误差小于容许值时的反 演波速场,即最终波速场,即为本发明待测量的岩石表面非均匀波速场。 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 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为本发明获得激光往返时间示意图; Fig. 4a obtains laser round-trip time for the present invention 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
所述的激光测距雷达1位于精密电动舵机2上方,可测量激光测距雷达1与标志杆4的距离;The
所述的精密电动舵机2可实现360度高速旋转,其上方固定有激光测距雷达1,下方固定有振动接收传感器3,主体连接可旋转滑轨8,使用时可带动可旋转滑轨8共同旋转;The precision
所述的振动接收传感器3主要用于接收冲压器7撞击岩石表面发出的振动信号,其上方固定有精密电动舵机2;The
所述的标志杆4上有竖条状标志物,可被激光测距雷达1识别,其下部嵌入到电动滑块5中且下方连接有压力传感器6;There is a vertical strip-shaped marker on the
所述的电动滑块5嵌入在可旋转滑轨8上,可实现在可旋转滑轨8上的滑动且位置可控,其上方下方分别连接有标志杆4和压力传感器6;The
所述的压力传感器6可接收冲压器7撞击岩石表面产生的压力脉冲,其安置在电动滑块5和冲压器7之间;The
所述的冲压器7可实现岩石表面撞击,撞击时可带动标志杆4和压力传感器6共同移动;The
所述的可旋转滑轨8与精密电动舵机2连接,可实现360度旋转,其上嵌入有电动滑块5,电动滑块5可在可旋转滑轨8上自由滑动。The
如图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
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
S3、压力、振动与距离信号采集与处理。如图4a到图4c所示,冲压器7在冲击岩石表
面后,压力传感器6会接收到冲击压力脉冲信号,该信号可通过信号转换公式获得AIC值,通
过分析最小AIC值对应的时间可获得冲击信号的起始时间;压力传感器6接收到冲击压力
脉冲信号后,激光测距雷达1发出激光束阵列,遇到标志杆4后发生反射后被激光测距雷达1
接收,即可获得激光往返时间;冲压器7撞击岩石表面的振动波传播到振动接收传感器3
后,振动接收传感器3会接收到振动脉冲信号,该信号可通过信号转换公式获得AIC值,通过
分析最小AIC值对应的时间可获得振动信号的到达时间。
S3. Acquisition and processing of pressure, vibration and distance signals. As shown in Figures 4a to 4c, after the
进一步地,所述的信号转换公式包括:Further, the signal conversion formula includes:
其中是当前时刻的AIC值,是时间采样点总数,是从初始时刻到当前时刻 的采样点,和是从初始时刻到当前时刻的采样点的平均值和方差,和是从当前时 刻到最终时刻的采样点的平均值和方差。 in is the AIC value at the current moment, is the total number of time sampling points, is the sampling point from the initial moment to the current moment, and is the mean and variance of the sampling points from the initial moment to the current moment, and 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:
其中是单点波速值,是激光在空气中的传播速度。 in is the single-point wave velocity value, 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
如图5b所示,当电动滑块5位于可旋转滑轨8距离激光测距雷达1较远处时,此时位于C8监测轨迹,单个测点9为C8监测轨迹上某点。精密电动舵机2可带动可旋转滑轨8以较小的步进角沿360度高速旋转,此时重复步骤二到四,可获得C8监测轨迹上所有测点的波速值。As shown in Figure 5b, when the
进一步地,当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
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;
根据实际测点波速开展波速场正反演;所述的实际测点波速为圆形点阵波速 在每个栅格中的线性插值;所述的波速场正反演包括波速场正演和波速场反演;所述的波 速场正演可通过程函方法开展;所述的波速场反演通过计算所述的波速场正演获得的波速 场和实际测点波速的差值,不断进行误差传播和参数优化,获得相对误差小于较小容许值 时的反演波速场,即最终波速场;所述的最终波速场即为本发明待测量的岩石表面非均匀 波速场。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 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)
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)
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)
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 |
-
2023
- 2023-01-05 CN CN202310010317.2A patent/CN115685324B/en active Active
Patent Citations (3)
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)
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 |