CN110907897A - 一种适用于含孔洞岩石的声发射源定位方法 - Google Patents
一种适用于含孔洞岩石的声发射源定位方法 Download PDFInfo
- Publication number
- CN110907897A CN110907897A CN201911334569.0A CN201911334569A CN110907897A CN 110907897 A CN110907897 A CN 110907897A CN 201911334569 A CN201911334569 A CN 201911334569A CN 110907897 A CN110907897 A CN 110907897A
- Authority
- CN
- China
- Prior art keywords
- acoustic emission
- hole
- grid
- sensor
- arrival time
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000012360 testing method Methods 0.000 claims abstract description 33
- 239000011159 matrix material Substances 0.000 claims description 40
- 239000013598 vector Substances 0.000 claims description 39
- 238000004458 analytical method Methods 0.000 claims description 7
- BJAJDJDODCWPNS-UHFFFAOYSA-N dotp Chemical compound O=C1N2CCOC2=NC2=C1SC=C2 BJAJDJDODCWPNS-UHFFFAOYSA-N 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000006378 damage Effects 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001687 destabilization Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000035882 stress Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
- G01S5/26—Position of receiver fixed by co-ordinating a plurality of position lines defined by path-difference measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6161—Seismic or acoustic, e.g. land or sea measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Geophysics (AREA)
- Radar, Positioning & Navigation (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明提供的一种适用于含孔洞岩石的声发射源定位方法,解决了板状岩石试件中存在孔洞条件下,声发射波绕射孔洞导致传播路径增加引起声发射源定位误差加大的难点问题,能够在确定孔洞半径的基础上,结合声发射波传播路径的解析分析,求解含孔洞岩石中声发射波的传播路径长度,结合声发射波到时、网格搜索及误差权重系数法方法可较为精确的确定声发射源的位置。
Description
技术领域
本发明属于岩石力学中声发射源定位方法技术领域,涉及一种适用于适用于含孔洞岩石的声发射源定位方法。
背景技术
在荷载等因素作用下,岩石内部会形成破裂面并同时伴随应力波的传播,这种现象被称为声发射。声发射是岩石破裂过程的伴生现象,与岩石的物理力学行为具有密切的关系,因此,可通过对声发射监测信息的分析来推断岩石的受力状态及破裂程度,进而预警与控制岩石体的失稳性破坏。
近年来声发射监测技术已逐渐被世界各国所重视,并在国内外大量室内岩石实验中得到应用。声发射源的定位是该技术的核心功能之一,能否进行精确的定位是声发射监测系统是否有效发挥作用的关键评价指标。声发射定位主要通过时差测量方法实现。在岩石力学实验中,为了模拟隧道开挖,往往在岩石试件中加工出一个或多个孔洞来研究孔洞周围岩石介质的受力、损伤及破坏过程,并使用声发射系统来定位损伤及破裂的位置。目前广泛应用的定位算法没有考虑空区对波传播路径的影响。
由于孔洞的存在,大量岩石内部的破裂(声发射源)位置与声发射传感器之间存在孔洞阻隔,声发射源诱发的声发射波必然要绕过孔洞才能传播至声发射传感器位置被采集到,波传播路径及传播路程的改变会导致波传播时长的改变,忽略这一影响因素会导致定位结果误差的增加。
发明内容
本发明提供的一种适用于含孔洞岩石的声发射源定位方法,解决了板状岩石试件中存在孔洞条件下,声发射波绕射孔洞导致传播路径增加引起声发射源定位误差加大的难点问题,能够在确定孔洞半径的基础上,结合声发射波传播路径的解析分析,求解含孔洞岩石中声发射波的传播路径长度,结合声发射波到时、网格搜索及误差权重系数法方法可较为精确的确定声发射源的位置。
本发明的一种适用于含孔洞岩石的声发射源定位方法,其特征在于,包括如下步骤:
步骤1、建立含孔洞及传感器位置的网格模型:建立岩石试件模型并确定声发射传感器在岩石试件上的位置,划分网格,去除网格中心点位于孔洞内的网格;
步骤2、确定声发射波的传播路径:判断网格中心及声发射传感器的连线是否切割孔洞,若不切割则该线段即为声发射波传播路径,若切割则计算利用解析方法计算声发射波绕射孔洞的传播路径,进而建立声发射波自网格中心到声发射传感器的路径长度矩阵;
步骤3、声发射源的定位:基于岩石波速测试结果及步骤2建立的路径长度矩阵确定理论到时差矩阵,利用声发射传感器接受的真实声发射波到时建立真实到时差矩阵/向量,匹配搜索理论到时差矩阵与真实到时差矩阵,进而确定声发射源的位置。
作为本发明的进一步优化,所述步骤1中建立含孔洞及传感器位置的网格模型具体包括:
步骤1.1、根据岩石试件的形状及尺寸建立岩石试件模型,根据孔洞中心的位置及半径建立孔洞模型,并将两个模型做差集布尔运算,完成含孔洞岩石试件模型的建立,并确定声发射传感器在岩石试件上的坐标;
步骤1.2、对含空区的数值模型进行网格划分,网格尺度满足以下条件:
式中,LM为网格尺度,单位mm;VP为P波在岩体中的传播速度,单位mm/s;Sf为声发射系统的采样频率,单位Hz;emax为满足定位要求的最大误差,单位mm;
步骤1.3、计算所有网格中心点的坐标,并将中心点与孔洞中心点间距小于孔洞半径的网格去除:
式中,xi、yi为网格i的中心点坐标,单位mm;i=l,2,…,N, N为划分的网格总数;Xc、Yc为孔洞中心点的坐标,单位mm;R为孔洞的半径,单位mm。
作为本发明的进一步优化,所述步骤2中确定声发射波的传播路径:
步骤2.1、依次连接各网格中心点与声发射传感器坐标,该线段所在直线的表达式为:
(yi-Ysj)x+(Xsj-xi)y+xiYsj-yiXsj=0
式中,Xsj、Ysj为声发射传感器j的坐标,单位mm;j=l,2,…, n,n为声发射传感器的数量。
并计算该线段所在直线距孔洞中心点的距离dij:
式中,dij为该线段所在直线距孔洞中心点的距离,单位mm;
步骤2.2、若该线段所在直线距孔洞中心点的距离大于或等于孔洞半径(dij≥R),则连接网格中心点i及声发射传感器j的线段不切割孔洞,该线段即为声发射波自网格中心点i到声发射传感器j的传播路径:
式中,lij为声发射波自网格中心点i到声发射传感器j的传播路径,单位mm;
若该线段所在直线距孔洞中心点的距离小于孔洞半径(dij<R),则分别计算以网格中心点i为起点、以孔洞中心及声发射传感器j为终点的两向量的内积,以及以声发射传感器j为起点、以孔洞中心及网格中心点i为终点的两向量的内积:
dotp1ij=(Xc-xi)(Xsj-xi)+(Yc-yi)(Ysj-yi)
dotp2ji=(Xc-Xsj)(xi-Xsj)+(Yc-Ysj)(yi-Ysj)
步骤2.3、若该线段所在直线距孔洞中心点的距离小于孔洞半径 (dij<R)且步骤2.2中计算的两内积(dotp1ij,dotp2ij)不全为正,则连接网格中心点i及声发射传感器j的线段不切割孔洞,该线段即为声发射波自网格中心点i到声发射传感器j的传播路径;
若该线段所在直线距孔洞中心点的距离小于孔洞半径(dij<R)且步骤2.2中计算的两内积(dotp1ij,dotp2ij)全为正,则连接网格中心点i及声发射传感器j的线段切割孔洞,自网格中心点i到声发射传感器j的声发射波的传播路径会绕射孔洞;步骤2.4、利用解析方法计算自网格中心点i到声发射传感器j的声发射波绕射孔洞的传播路径,分别计算网格中心点i、声发射传感器j至孔洞的切线长度:
式中,l1ij、l2ji分别为网格中心点i、声发射传感器j至孔洞的切线长度,单位mm;
然后计算两个相切点间圆弧的长度:
式中,lcij为两个相切点间圆弧的长度,单位mm;
两切线长度与圆弧长度之和即为自网格中心点i到声发射传感器 j的声发射波绕射孔洞的传播路径:
lij=lcij+l1ij+l2ji
步骤2.5、依次计算各个网格中心至各个声发射传感器的传播路径长度,并建立声发射波自各网格中心点到各声发射传感器的路径长度矩阵Lij。
作为本发明的进一步优化,所述步骤3中声发射源的定位,具体包括以下步骤:
步骤3.1、根据路径长度矩阵Lij及P波在岩体中的传播速度确定理论走时矩阵:
Tij=Lij/VP
式中,Tij为声发射波自网格中心点i传播至声发射传感器j的理论走时,单位为s;
步骤3.2、根据理论计算得到的理论走时矩阵,计算依次从各个网格中心点传播至各个传感器位置的理论到时差矩阵:
ΔTij=Tij-Timin
式中,ΔTij为理论到时差矩阵,单位为s;Timin为理论走时矩阵 Tij中第i行中最小的元素,单位为s;
步骤3.3、真实到时差向量的建立:利用在岩石试件上安装的声发射传感器采集波形数据,对传感器接收到的波形进行人工或自动到时拾取工作,并根据下式计算真实到时差向量:
ΔTjreal=Tjreal-min(Tjreal)
式中,ΔTjreal为真实到时差向量,单位s;Tjreal为传感器j的真实到时,单位s;
步骤3.4、匹配搜索及声发射源定位:将真实到时差向量与理论到时差矩阵中的行向量进行匹配搜索,如果有偏差量绝对值之和为0 的行向量,则直接确定声发射源的坐标与该行向量对应的网格中心点坐标相同,完成定位;如果没有,则确定与真实到时差向量的总偏差量最小的4个行向量,将偏差量的大小作为权重系数来确定最终的定位点:
式中,Loc为声发射源的定位坐标,单位mm;ek为真实到时差与理论到时差间的总偏差量,单位s;P1、P2、P3、P4分别为与真实到时差向量最接近的4个网格中心点坐标。
本发明的一种适用于含孔洞岩石的声发射源定位方法,至少具有以下有益效果:
该方法解决了板状岩石试件中存在孔洞条件下,声发射波绕射孔洞导致传播路径增加引起声发射源定位误差加大的难点问题。将已存在的孔洞作为已知条件,结合声发射波传播路径的解析分析将孔洞对波传播路径的改变及由此造成的波传播路径长度的变化、波到达传感器的时间加以量化,显著降低了因孔洞存在诱发的声发射源定位误差。另外,在定位算法中联合使用了网格搜索及误差权重系数法,既避免了迭代算法中过分依赖初始迭代值和不收敛的缺点,又突破了网格尺度对定位结果的限制。
附图说明
图1为建立的含孔洞及传感器位置的网格模型,其中图1a为含孔洞、声发射传感器的岩石试件的模型,图1b为划分网格后的岩石试件模型,图1c为去除中心点在孔洞内的网格后的岩石试件模型。
图2为自网格中心点到声发射传感器的线段与孔洞间的相对位关系,其中图2a为线段及其所在的直线均与孔洞相离,图2b为线段与孔洞相离但线段所在的直线与孔洞相交,图2c为线段及其所在的直线均与孔洞相交。
图3为实测声发射波形示意图,其中图3a为传感器S1接收到的声发射波形,图3b为传感器S2接收到的声发射波形,图3c为传感器S3接收到的声发射波形,图3d为传感器S4接收到的声发射波形。
具体实施方式
以下结合附图和具体实施例对本发明作进一步的介绍。
本发明涉及的一种适用于含孔洞岩石的声发射源定位方法,包括如下步骤:
步骤1、建立含孔洞及传感器位置的网格模型:建立岩石试件模型并确定声发射传感器在岩石试件上的位置,划分网格,去除网格中心点位于孔洞内的网格;
所述步骤1中建立含孔洞及传感器位置的网格模型具体包括:
步骤1.1、根据岩石试件的形状及尺寸建立岩石试件模型,根据孔洞中心的位置及半径建立孔洞模型,并将两个模型做差集布尔运算,完成含孔洞岩石试件模型的建立,并确定声发射传感器在岩石试件上的坐标,如图1a所示;
步骤1.2、对含空区的数值模型进行网格划分,可将靠近孔洞的网格尺度划分的较小,而越远离孔洞的网格尺度划分的较大,如图 1b所示,孔洞周围的网格尺度为:
式中,LM为网格尺度,单位mm;VP为P波在岩体中的传播速度,单位mm/s;Sf为声发射系统的采样频率,单位Hz;emax为满足定位要求的最大误差,单位mm。
步骤1.3、计算所有网格中心点的坐标,并将中心点与孔洞中心点间距小于孔洞半径的网格去除,如图1c所示:
式中,xi、yi为网格i的中心点坐标,单位mm;i=l,2,…,N,N为划分的网格总数;Xc、Yc为孔洞中心点的坐标,单位mm;R为孔洞的半径,单位mm。
步骤2、确定声发射波的传播路径:判断网格中心及声发射传感器的连线是否切割孔洞,若不切割则该线段即为声发射波传播路径,若切割则计算利用解析方法计算声发射波绕射孔洞的传播路径,进而建立声发射波自网格中心到声发射传感器的路径长度矩阵;
所述步骤2中确定声发射波的传播路径具体包括:
步骤2.1、利用线段依次连接各网格中心点与声发射传感器坐标,该线段所在直线的表达式为:
(yi-Ysj)x+(Xsj-xi)y+xiYsj-yiXsj=0
式中,Xsj、Ysj为声发射传感器j的坐标,单位mm;j=l,2,…, n,n为声发射传感器的数量。
并计算该线段所在直线距孔洞中心点的距离dij:
式中,dij为该线段所在直线距孔洞中心点的距离,单位mm;
步骤2.2、若该线段所在直线距孔洞中心点的距离大于或等于孔洞半径(dij≥R),如图2a所示,则连接网格中心点i及声发射传感器j的线段不切割孔洞,该线段即为声发射波自网格中心点i到声发射传感器j的传播路径(图2a中粗线所示):
式中,lij为声发射波自网格中心点i到声发射传感器j的传播路径,单位mm;
若该线段所在直线距孔洞中心点的距离小于孔洞半径(dij<R),则分别计算以网格中心点i为起点、以孔洞中心及声发射传感器j为终点的两向量的内积,以及以声发射传感器j为起点、以孔洞中心及网格中心点i为终点的两向量的内积:
dotp1ij=(Xc-xi)(Xsj-xi)+(Yc-yi)(Ysj-yi)
dotp2ji=(Xc-Xsj)(xi-Xsj)+(Yc-Ysj)(yi-Ysj)
步骤2.3、若该线段所在直线距孔洞中心点的距离小于孔洞半径 (dij<R)且步骤2.2中计算的两内积(dotp1ij,dotp2ij)不全为正,如图2b所示,则连接网格中心点i及声发射传感器j的线段不切割孔洞,该线段即为声发射波自网格中心点i到声发射传感器j的传播路径,
如图2b中粗线所示。若该线段所在直线距孔洞中心点的距离小于孔洞半径(dij<R)且步骤2.2中计算的两内积(dotp1ij,dotp2ij)全为正,如图2c所示,则连接网格中心点i及声发射传感器j的线段切割孔洞,自网格中心点i到声发射传感器j的声发射波的传播路径会绕射孔洞。
步骤2.4、利用解析方法计算自网格中心点i到声发射传感器j 的声发射波绕射孔洞的传播路径,分别计算网格中心点i、声发射传感器j至孔洞的切线长度:
式中,l1ij、l2ji分别为网格中心点i、声发射传感器j至孔洞的切线长度,单位mm;
然后计算两个相切点间圆弧的长度:
式中,lcij为两个相切点间圆弧的长度,单位mm;
两切线长度与圆弧长度之和即为自网格中心点i到声发射传感器 j的声发射波绕射孔洞的传播路径(图2c中粗线所示):
lij=lcij+l1ij+l2ji
步骤2.5、依次计算各个网格中心至各个声发射传感器的传播路径长度,并建立声发射波自各网格中心点到各声发射传感器的路径长度矩阵Lij。
步骤3、声发射源的定位:基于岩石波速测试结果及步骤2建立的路径长度矩阵确定理论到时差矩阵,利用声发射传感器接受的真实声发射波到时建立真实到时差矩阵/向量,匹配搜索理论到时差矩阵与真实到时差矩阵,进而确定声发射源的位置。
所述步骤3中声发射源的定位具体包括:
步骤3.1、根据路径长度矩阵Lij及P波在岩体中的传播速度确定理论走时矩阵:
Tij=Lij/VP
式中,Tij为声发射波自网格中心点i传播至声发射传感器j的理论走时,单位为s。
步骤3.2、根据理论计算得到的理论走时矩阵,计算依次从各个网格中心点传播至各个传感器位置的理论到时差矩阵:
ΔTij=Tij-Timin
式中,ΔTij为理论到时差矩阵,单位为s;Timin为理论走时矩阵 Tij中第i行中最小的元素,单位为s。
步骤3.3、真实到时差向量的建立:利用在岩石试件上安装的声发射传感器采集波形数据(如图3所示),对传感器接收到的波形进行人工或自动到时拾取工作,并根据下式计算真实到时差向量:
ΔTjreal=Tjreal-min(Tjreal)
式中,ΔTjreal为真实到时差向量,单位s;Tjreal为传感器j的真实到时,单位s;
步骤3.4、匹配搜索及声发射源定位:将真实到时差向量与理论到时差矩阵中的行向量进行匹配搜索,如果有偏差量绝对值之和小于满足定位要求最大误差的行向量,则直接确定声发射源的坐标与该行向量对应的网格中心点坐标相同,完成定位;如果没有,则确定与真实到时差向量最接近(总偏差量最小)的4个行向量,将偏差量的大小作为权重系数来确定最终的定位点:
式中,Loc为声发射源的定位坐标,单位mm;ek为真实到时差与理论到时差间的总偏差量,单位s;P1、P2、P3、P4分别为与真实到时差向量最接近的4个网格中心点坐标。
实施例:
1、根据实际岩石试件的形状及尺寸建立岩石试件模型,岩石试件尺度为300×300mm,孔洞半径为50mm,如图1a所示,假定有4 个传感器S1、S2、S3、S4,传感器位置坐标如表1所示(岩石试件的左下角为坐标原点)。
表1
传感器 | x(mm) | y(mm) |
S1 | 70 | 70 |
S2 | 230 | 70 |
S3 | 70 | 230 |
S4 | 230 | 230 |
2、对含孔洞的岩石试件模型进行网格划分,如图1b所示。
3、去除网格中心点位于孔洞内的网格,如图1c所示。
4、依次计算网格中心点到各个声发射传感器的波传播路径长度,在有4个传感器的情况下,每个网格中心点会有4个直达传感器的传播路径,若有N个网格中心点就会形成N×4列的路径长度矩阵。以坐标为(78.75,116.25)的网格中心点为例,其到各个声发射传感器的波传播路径长度矩阵,如表2所示。
表2
传感器 | 波传播路径长度(mm) |
S1 | 37.72 |
S2 | 154.02 |
S3 | 104.77 |
S4 | 199.26 |
5、测量岩石试件的波速为3680000mm/s,结合波传播路径长度确定理论走时矩阵,通过减去走时最少的值可以得到理论到时差矩阵,如表3所示。
表3
6、在坐标(78.75,116.25)位置进行人工断铅,利用在岩石试件上安装的声发射传感器采集波形数据(如图3所示),对传感器接收到的波形进行人工到时拾取工作,并计算真实到时差,如表4所示。
传感器 | 实际走时(s) | 实际到时差(s) |
S1 | 0.0000102 | 0 |
S2 | 0.0000418 | 0.0000316 |
S3 | 0.0000284 | 0.0000182 |
S4 | 0.0000542 | 0.000044 |
7、将真实到时差与理论到时差进行匹配,发现真实的到时差与模拟到时差时的偏差量绝对值之和小于预设的最大误差5×10-7,因此确定实际人工断铅产生的声发射源的坐标为(78.75,116.25)。
以上所述仅为本发明的较佳实施例,并不用以限制本发明的思想,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种适用于含孔洞岩石的声发射源定位方法,其特征在于,包括如下步骤:
步骤1、建立含孔洞及传感器位置的网格模型:建立岩石试件模型并确定声发射传感器在岩石试件上的位置,划分网格,去除网格中心点位于孔洞内的网格;
步骤2、确定声发射波的传播路径:判断网格中心及声发射传感器的连线是否切割孔洞,若不切割则该线段即为声发射波传播路径,若切割则计算利用解析方法计算声发射波绕射孔洞的传播路径,进而建立声发射波自网格中心到声发射传感器的路径长度矩阵;
步骤3、声发射源的定位:基于岩石波速测试结果及步骤2建立的路径长度矩阵确定理论到时差矩阵,利用声发射传感器接受的真实声发射波到时建立真实到时差矩阵/向量,匹配搜索理论到时差矩阵与真实到时差矩阵,进而确定声发射源的位置。
2.如权利要求1所述的一种适用于含孔洞岩石的声发射源定位方法,其特征在于,所述步骤1中建立含孔洞及传感器位置的网格模型具体包括:
步骤1.1、根据岩石试件的形状及尺寸建立岩石试件模型,根据孔洞中心的位置及半径建立孔洞模型,并将两个模型做差集布尔运算,完成含孔洞岩石试件模型的建立,并确定声发射传感器在岩石试件上的坐标;
步骤1.2、对含空区的数值模型进行网格划分,网格尺度满足以下条件:
式中,LM为网格尺度,单位mm;VP为P波在岩体中的传播速度,单位mm/s;Sf为声发射系统的采样频率,单位Hz;emax为满足定位要求的最大误差,单位mm;
步骤1.3、计算所有网格中心点的坐标,并将中心点与孔洞中心点间距小于孔洞半径的网格去除:
式中,xi、yi为网格i的中心点坐标,单位mm;i=l,2,…,N,N为划分的网格总数;Xc、Yc为孔洞中心点的坐标,单位mm;R为孔洞的半径,单位mm。
3.如权利要求1所述的一种适用于含孔洞岩石的声发射源定位方法,其特征在于,所述步骤2中确定声发射波的传播路径:
步骤2.1、依次连接各网格中心点与声发射传感器坐标,该线段所在直线的表达式为:
(yi-Ysj)x+(Xsj-xi)y+xiYsj-yiXsj=0
式中,Xsj、Ysj为声发射传感器j的坐标,单位mm;j=l,2,…,n,n为声发射传感器的数量。
并计算该线段所在直线距孔洞中心点的距离dij:
式中,dij为该线段所在直线距孔洞中心点的距离,单位mm;
步骤2.2、若该线段所在直线距孔洞中心点的距离大于或等于孔洞半径(dij≥R),则连接网格中心点i及声发射传感器j的线段不切割孔洞,该线段即为声发射波自网格中心点i到声发射传感器j的传播路径:
式中,lij为声发射波自网格中心点i到声发射传感器j的传播路径,单位mm;
若该线段所在直线距孔洞中心点的距离小于孔洞半径(dij<R),则分别计算以网格中心点i为起点、以孔洞中心及声发射传感器j为终点的两向量的内积,以及以声发射传感器j为起点、以孔洞中心及网格中心点i为终点的两向量的内积:
dotp1ij=(Xc-xi)(Xsj-xi)+(Yc-yi)(Ysj-yi)
dotp2ji=(Xc-Xsj)(xi-Xsj)+(Yc-Ysj)(yi-Ysj)
步骤2.3、若该线段所在直线距孔洞中心点的距离小于孔洞半径(dij<R)且步骤2.2中计算的两内积(dotp1ij,dotp2ij)不全为正,则连接网格中心点i及声发射传感器j的线段不切割孔洞,该线段即为声发射波自网格中心点i到声发射传感器j的传播路径;
若该线段所在直线距孔洞中心点的距离小于孔洞半径(dij<R)且步骤2.2中计算的两内积(dotp1ij,dotp2ij)全为正,则连接网格中心点i及声发射传感器j的线段切割孔洞,自网格中心点i到声发射传感器j的声发射波的传播路径会绕射孔洞;
步骤2.4、利用解析方法计算自网格中心点i到声发射传感器j的声发射波绕射孔洞的传播路径,分别计算网格中心点i、声发射传感器j至孔洞的切线长度:
式中,l1ij、l2ji分别为网格中心点i、声发射传感器j至孔洞的切线长度,单位mm;
然后计算两个相切点间圆弧的长度:
式中,lcij为两个相切点间圆弧的长度,单位mm;
两切线长度与圆弧长度之和即为自网格中心点i到声发射传感器j的声发射波绕射孔洞的传播路径:
lij=lcij+l1ij+l2ji
步骤2.5、依次计算各个网格中心至各个声发射传感器的传播路径长度,并建立声发射波自各网格中心点到各声发射传感器的路径长度矩阵Lij。
4.如权利要求1所述的一种适用于含孔洞岩石的声发射源定位方法,其特征在于,所述步骤3中声发射源的定位,具体包括以下步骤:
步骤3.1、根据路径长度矩阵Lij及P波在岩体中的传播速度确定理论走时矩阵:
Tij=Lij/VP
式中,Tij为声发射波自网格中心点i传播至声发射传感器j的理论走时,单位为s;
步骤3.2、根据理论计算得到的理论走时矩阵,计算依次从各个网格中心点传播至各个传感器位置的理论到时差矩阵:
ΔTij=Tij-Timin
式中,ΔTij为理论到时差矩阵,单位为s;Timin为理论走时矩阵Tij中第i行中最小的元素,单位为s;
步骤3.3、真实到时差向量的建立:利用在岩石试件上安装的声发射传感器采集波形数据,对传感器接收到的波形进行人工或自动到时拾取工作,并根据下式计算真实到时差向量:
式中,ΔTjreal为真实到时差向量,单位s;Tjreal为传感器j的真实到时,单位s;
步骤3.4、匹配搜索及声发射源定位:将真实到时差向量与理论到时差矩阵中的行向量进行匹配搜索,如果有偏差量绝对值之和为0的行向量,则直接确定声发射源的坐标与该行向量对应的网格中心点坐标相同,完成定位;如果没有,则确定与真实到时差向量的总偏差量最小的4个行向量,将偏差量的大小作为权重系数来确定最终的定位点:
式中,Loc为声发射源的定位坐标,单位mm;ek为真实到时差与理论到时差间的总偏差量,单位s;P1、P2、P3、P4分别为与真实到时差向量最接近的4个网格中心点坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911334569.0A CN110907897B (zh) | 2019-12-23 | 2019-12-23 | 一种适用于含孔洞岩石的声发射源定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911334569.0A CN110907897B (zh) | 2019-12-23 | 2019-12-23 | 一种适用于含孔洞岩石的声发射源定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110907897A true CN110907897A (zh) | 2020-03-24 |
CN110907897B CN110907897B (zh) | 2023-09-15 |
Family
ID=69827082
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911334569.0A Active CN110907897B (zh) | 2019-12-23 | 2019-12-23 | 一种适用于含孔洞岩石的声发射源定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110907897B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111398433A (zh) * | 2020-04-17 | 2020-07-10 | 中南大学 | 一种基于线性加权最小二乘法的声发射源定位方法及系统 |
CN116593295A (zh) * | 2023-07-19 | 2023-08-15 | 北京科技大学 | 利用岩石各向异性波速提高声发射定位精度的方法及装置 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101035681A (zh) * | 2004-10-01 | 2007-09-12 | 拉伯赛特股份有限公司 | 与声发射有关的流体特性的声学评诂 |
DE102008017001B3 (de) * | 2008-04-03 | 2009-07-16 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Vorrichtung zur Bestimmung der akustischen Position von Mikrofonen im Raum |
US20090245019A1 (en) * | 2008-03-31 | 2009-10-01 | Jon Falkenberg | Method and system for determining geodetic positions of towed marine sensor array components |
CN101710105A (zh) * | 2009-09-24 | 2010-05-19 | 山东大学 | 多层包扎容器缺陷的声发射确定方法 |
CN102680579A (zh) * | 2012-06-01 | 2012-09-19 | 中国计量学院 | 基于声发射传感器阵列的混凝土内部缺陷检测方法 |
CN102928818A (zh) * | 2012-10-18 | 2013-02-13 | 东南大学 | 一种基于近场波束形成的碰摩声发射源的定位方法 |
CN103033565A (zh) * | 2012-12-07 | 2013-04-10 | 山东大学 | 岩样破裂过程的声发射与电阻率联合监测装置及监测方法 |
CN104914168A (zh) * | 2015-06-23 | 2015-09-16 | 中国科学院上海硅酸盐研究所 | 一种用于热障涂层高温声发射信号检测的工装 |
CN108646221A (zh) * | 2018-06-14 | 2018-10-12 | 河北工业大学 | 一种基于声源阵列的空间麦克风定位方法 |
CN109738940A (zh) * | 2019-01-18 | 2019-05-10 | 东北大学 | 一种存在空区条件下的声发射/微震事件定位方法 |
-
2019
- 2019-12-23 CN CN201911334569.0A patent/CN110907897B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101035681A (zh) * | 2004-10-01 | 2007-09-12 | 拉伯赛特股份有限公司 | 与声发射有关的流体特性的声学评诂 |
US20090245019A1 (en) * | 2008-03-31 | 2009-10-01 | Jon Falkenberg | Method and system for determining geodetic positions of towed marine sensor array components |
DE102008017001B3 (de) * | 2008-04-03 | 2009-07-16 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Vorrichtung zur Bestimmung der akustischen Position von Mikrofonen im Raum |
CN101710105A (zh) * | 2009-09-24 | 2010-05-19 | 山东大学 | 多层包扎容器缺陷的声发射确定方法 |
CN102680579A (zh) * | 2012-06-01 | 2012-09-19 | 中国计量学院 | 基于声发射传感器阵列的混凝土内部缺陷检测方法 |
CN102928818A (zh) * | 2012-10-18 | 2013-02-13 | 东南大学 | 一种基于近场波束形成的碰摩声发射源的定位方法 |
CN103033565A (zh) * | 2012-12-07 | 2013-04-10 | 山东大学 | 岩样破裂过程的声发射与电阻率联合监测装置及监测方法 |
CN104914168A (zh) * | 2015-06-23 | 2015-09-16 | 中国科学院上海硅酸盐研究所 | 一种用于热障涂层高温声发射信号检测的工装 |
CN108646221A (zh) * | 2018-06-14 | 2018-10-12 | 河北工业大学 | 一种基于声源阵列的空间麦克风定位方法 |
CN109738940A (zh) * | 2019-01-18 | 2019-05-10 | 东北大学 | 一种存在空区条件下的声发射/微震事件定位方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111398433A (zh) * | 2020-04-17 | 2020-07-10 | 中南大学 | 一种基于线性加权最小二乘法的声发射源定位方法及系统 |
CN116593295A (zh) * | 2023-07-19 | 2023-08-15 | 北京科技大学 | 利用岩石各向异性波速提高声发射定位精度的方法及装置 |
CN116593295B (zh) * | 2023-07-19 | 2023-10-03 | 北京科技大学 | 利用岩石各向异性波速提高声发射定位精度的方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN110907897B (zh) | 2023-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109738940B (zh) | 一种存在空区条件下的声发射/微震事件定位方法 | |
CN108710732B (zh) | 一种盾构刀具服役期磨损的评估和预测方法 | |
CN102129063B (zh) | 一种微震源或声发射源的定位方法 | |
Li et al. | Elastic mechanics-based fixturing scheme optimization of variable stiffness structure workpieces for surface quality improvement | |
CN110907897A (zh) | 一种适用于含孔洞岩石的声发射源定位方法 | |
CN109141266B (zh) | 一种钢结构测量方法及系统 | |
Rui et al. | A novel AE source localization method using clustering detection to eliminate abnormal arrivals | |
CN112924065A (zh) | 一种基于盲孔法测曲面残余应力的测量方法 | |
CN106093207A (zh) | 一种基于非线性无味卡尔曼滤波算法的兰姆波损伤定位方法 | |
CN110082559A (zh) | 一种喷丸过程中的速度测量装置及速度测量方法 | |
CN110083988B (zh) | 一种船舶水下辐射噪声评估方法 | |
CN110018062B (zh) | 一种直剪试验中岩石结构面剪切破坏位置定位方法 | |
CN105022091B (zh) | 一种无预测速的远场震源快速定位方法 | |
CN110580391B (zh) | 一种柔性结构的基频模态测量方法 | |
Ozevin | Geometry-based spatial acoustic source location for spaced structures | |
CN109085248B (zh) | 承压管道壁冲击源的定位方法、装置及系统 | |
CN105004795B (zh) | 伪缺陷信号识别及利用其提高管道无损检测精度的方法 | |
CN105548356A (zh) | 基于t模态导波检测带环焊缝小弯曲半径弯管缺陷的方法 | |
CN109933921B (zh) | 一种滚石灾害风险评估方法、装置、系统及存储介质 | |
CN112925018B (zh) | 一种异常压力成因判别方法及装置 | |
CN114578197A (zh) | 基于平面uhf传感器的变压器套管局放三维定位方法及系统 | |
CN114813962A (zh) | 一种基于gclm-网格搜索算法的声发射源定位方法 | |
CN108692806A (zh) | 固体中功率超声强度分布的测量方法 | |
CN104462790A (zh) | 疲劳耐久分析的自由表面方法 | |
CN113866827A (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 |