CN102262220A - 一种基于非线性拟合的微震源或声发射源的定位方法 - Google Patents

一种基于非线性拟合的微震源或声发射源的定位方法 Download PDF

Info

Publication number
CN102262220A
CN102262220A CN2011101093724A CN201110109372A CN102262220A CN 102262220 A CN102262220 A CN 102262220A CN 2011101093724 A CN2011101093724 A CN 2011101093724A CN 201110109372 A CN201110109372 A CN 201110109372A CN 102262220 A CN102262220 A CN 102262220A
Authority
CN
China
Prior art keywords
acoustic emission
emission source
sensor
microquake sources
source
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
Application number
CN2011101093724A
Other languages
English (en)
Other versions
CN102262220B (zh
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.)
Central South University
Original Assignee
Central South University
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 Central South University filed Critical Central South University
Priority to CN 201110109372 priority Critical patent/CN102262220B/zh
Publication of CN102262220A publication Critical patent/CN102262220A/zh
Application granted granted Critical
Publication of CN102262220B publication Critical patent/CN102262220B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开了一种基于非线性拟合的微震源或声发射源的定位方法,将多个声发射传感器安放于待测对象的附近,根据已知传感器所在位置的坐标值和P、S波到时,根据距离运算公式及P、S波速度与时间关系式,可通过加速单纯形法或麦夸特法等非线性最小二乘回归方法求出微震源或声发射源的位置坐标。本发明的基于非线性拟合的微震源或声发射源的定位方法不用测量或预先给定速度,能避免测量声速给定位造成的影响,提高了定位精度,在实际工程应用中较传统方法方便实用。

Description

一种基于非线性拟合的微震源或声发射源的定位方法
技术领域
[0001] 本发明涉及一种基于非线性拟合的微震源或声发射源的定位方法。 背景技术
[0002] 如何准确定位材料破损或岩爆前期伴随产生的微震/声发射源位置十分重要。微震/声发射的源定位技术是声发射技术研究的核心问题之一,其定位准确程度反映声源检测位置与真实位置的符合程度。
[0003] 许多科学工作者对声发射的源定位进行了大量的研究,提出了许多对声发射的源定位方法,如基于模态分析和小波变换的定位方法、基于神经元网络系统的定位方法等,但目前的声发射仪器进行定位仍普遍采用传统的第1次门槛跨越技术时差定位方法,该方法技术比较成熟,但需要预先给定声速或实测声速。波的传播速度受到传播介质的材质、尺寸及表面状态等因素影响,当输入的速度不同于被测物体中的真实速度时,特别对于岩石介质复杂性,将给系统定位带来误差,另外,实测声速受探头间距的大小影响也较大。
[0004] 总之,在目前声发射仪及微震定位系统中还不能消除速度测量偏差给定位造成的影响。为解决这一问题,有必要提出了一种不测速度的非线性定位的新方法,可以不用测量速度,能避免测量速度给定位造成的影响。
发明内容
[0005] 本发明所要解决的技术问题是提出一种基于非线性拟合的微震源或声发射源的定位方法,该基于非线性拟合的微震源或声发射源的定位方法不用测量或预先给定速度, 能避免测量声速给定位造成的影响,提高了定位精度,在实际工程应用中较传统方法方便实用。
[0006] 本发明的技术解决方案如下:
[0007] —种基于非线性拟合的微震源或声发射源的定位方法,在待测对象即微震源或声发射源的周边放置N个不处于同一平面的声发射传感器,N个声发射传感器的三维位置坐标已知,分另Ij为(X1, Y1, Z1)、(x2, J2, ζ2)、· · ·,(xN, yN, zN) ;N 大于等于 5 ;
[0008] 微震源或声发射源发声后,第i个声发射传感器感应到微震源或声发射源的发声讯号并分别记录下感应到该发声讯号的时刻,P波和S波到达该声发射传感器的时刻分别为tai,tM,共得到N组时间数据;
[0009] 对于公式:
[0010]
Figure CN102262220AD00031
其中,i和j取值为1,2,. . .,N,且i不等于j ; (x0, y0, Z0)为微震源或声发射源的坐标,参 αβ
数^^^^,其中α,β分别为P波、S波速度; α- β
[0011] 上式中,(x0, y0, Z0)和發为未知量,其余量为已知量;[0012] i和j不同组合形成C》组参数得到基于上式的C》个方程;G表示N个值中任2个值的组合数;
[0013] 根据该C》个方程,以
Figure CN102262220AD00041
为因变量,以(X(l,yQ,Z(l)及發为未知量, 采用非线性拟合方法得到震源微震源或声发射源的坐标(χο,y。,z0)。
[0014] 所述的非线性拟合方法为加速单纯形法或麦夸特法(Marqimrdt)。
[0015] 本发明的技术构思为:在待测对象的附近放置N个不处于同一平面的声发射传感器,N个声发射传感器的三维位置坐标已知,分别为(Xl,Z1)、(x2, y2, z2).... , (xN, yN, zN) ;N大于等于5 ;
[0016] 微震源或声发射源发声后,台站的第i个传感器感应到微震源或声发射源的发声讯号并分别记录下感应到该发声讯号的时刻,P波S波到时分别为tai,tM,震源到台站i 的距离为Ri。
[0017] 采用均勻介质模型,设P波、S波速度分别为a,β,假设参数
[0018]
Figure CN102262220AD00042
[0019]对传感器 i 有:(χ0 - X1 )2 + (y0 -兄.)2 + (Z0 - Z1 f = Rf
[0020]对传感器 j 有:(x。-x;)2+(y0-力)2 + (ζ。-z])2 = R2j
[0021] 将方程O),⑶求和,可得:
[0022]
Figure CN102262220AD00043
[0023] 将方程(4)用于5个以上的传感器(包括5个),即得到一因变量为 (tM-tai)2+(tej-taj)2,关于(Xtl,%,Ztl)及伶的非线性拟合公式,从而拟合得到震源位置。
[0024] 有益效果:
[0025] 本发明提出了一种基于非线性拟合的微震源或声发射源的定位方法,本发明根据声发射源或微震源信号传播至传感器的时刻及传感器所监测点的坐标值,巧妙地根据P、S 波触发传感器时刻差与声发射或微震源坐标值的关系推导得到声发射或微震源的位置坐标公式。将传感器安放于待测对象的附近,至少需要四个传感器,安放的传感器不要在同一平面上。通过传感器测量声发射源或微震源信号到达安放点的时刻(见附图1)。因此建立起声发射源或微震源坐标,传感器所在位置坐标,各传感器之触发的P、S波时间差之间的关系式,通过数据拟合,即可以得到声发射源或微震源的位置坐标。
[0026] 本发明突破了过去传统的给定或测量速度的方法。波的传播速度受到传播介质的材质、尺寸、表面状态及探头间距等因素影响,当输人的速度不同于被测物体中的真实速度时,将给系统定位带来误差。本发明提出了一种非线性定位方法,可以不用测量速度,既能避免测量速度给定位带来的误差,又能补充单纯P波定位信息不全面的缺陷,经多处矿山的实际验证,采用本发明方法的定位的平均误差为5米,而现有定位方法的平均误差为30米。因而,本发明为声发射源或微震源的定位提供了一种简便宜行、精确可靠、实用的新方法。
附图说明
[0027] 图1是三维条件下微震源或声发射源的微震源或声发射源的定位示意图。
[0028] 标号说明:1-第一传感器,2-第二传感器,3-第三传感器,4-第四传感器,5-微震源或声发射源。
具体实施方式
[0029] 以下将结合附图和具体实施例对本发明做进一步详细说明:
[0030] 实施例1 :
[0031] 如图1所示,预设一微震源/声发射源的位置为QOOm,300m,600m),5个传感器的坐标(X1, y” Z1)、(x2, Y2, z2)、(x3, y3, z3)、(x4, Yi, z4)与(x5, y5, z5)分另ll为(0,0,0), (800, 0,0), (800,800,0),(0,800,0),(0,0,800),单位均为 m。设 P、S 波的波速分别为 5600m/s, 3100m/s。P波触发以上传感器的时刻分别为0. 125,0. 175872,0. 143969,0. 073627秒;S波触发以上传感器的时刻分别为 0. 225806,0. 290323,0. 317705,0. 260073,0. 133003 秒。
[0032] 以本例对三维定位问题进行详细说明,实际定位时,已知量为五个传感器的坐标及传感器触发记录的时刻,微震源/声发射源的位置未知,在这里之所以给出,目的是通过本专利提出的方法进行验证。具体实施步骤如下:
[0033] (1)将5个传感器布置在待测对象附近,且5个传感器的不在同一平面上,5个传感器的坐标(X1, Υι» Z1)、(x2, Y2, z2, )、(x3, y3, z3)、(x4, y4,Z4)、(x5, y5, z5)分另ll为(0,0, 0),(800,0,0), (800,800,0),(0,800,0), (0,0,800) 0 P 波触发以上传感器的时刻分别为 0. 125s,0. 175872s,0. 143969s,0. 073627s ;S 波触发以上传感器的时刻分别为 0. 225806s, 0.290323s,0.317705s,0. 260073s,0. 133003s。
[0034] (2)求出5个传感器接收到S与P波到时时间差,进而求得因变量 (tM-tai)2+(t0J-taJ)2 的值,共有 10 组。
[0035] (3)将以上(1)与⑵求出的10组值与相应传感器坐标,采用加速单纯形法或麦夸特法(Marqimrdt)中至少一种方法拟合下式中的4个未知量xQ,y0, z0, Vo:
Figure CN102262220AD00051
[0036]通过拟合可知 x0, y0, z0,扒分别为 199. 9935m, 299. 9998m, 599. 9988m, 6943. 9935, x0,y。,z0的值(199. 9935m, 299. 9998m, 599. 9988m)即为定位结果,可见定位结果与预设的坐标QOOm,300m,600m)吻合较好,定位精度高。

Claims (2)

1. 一种基于非线性拟合的微震源或声发射源的定位方法,其特征在于,在待测对象即微震源或声发射源的周边放置N个不处于同一平面的声发射传感器,N个声发射传感器的三维位置坐标已知,分别为(X1, Y1, Z1)、(x2, y2, Z2)、. . .,(xN, yN,zN) ;N大于等于5 ;微震源或声发射源发声后,第i个声发射传感器感应到微震源或声发射源的发声讯号并分别记录下感应到该发声讯号的时刻,P波和S波到达该声发射传感器的时刻分别为 tai,tM,共得到N组时间数据; 对于公式:(x0 - X. )2 + (J0 - J,. )2 + (z0 - Zi )2 + (x0 - X . )2 + (J0 - J . )2 + (z0 - Z. )2 、2 SL-^-Mn.)+(、_g);“中,i和j取值为1,2,. . .,N,且i不等于j ; (x0, y0, Z0)为微震源或声发射源的坐标,参数 αβ其中α,β分别为P波、S波速度;上式中,(¾,〜Ztl)和發为未知量,其余量为已OC — μ知量;i和j不同组合形成C》组参数得到基于上式的C》个方程;表示N个值中任2个值的组合数;根据该个方程,以(tei-tai)2+(tej-taj)2为因变量,以(χ。,y。,z0)及發为未知量,采用非线性拟合方法得到震源微震源或声发射源的坐标(xo,y。,z0)。
2.根据权利要求1所述的基于非线性拟合的微震源或声发射源的定位方法,其特征在于,所述的非线性拟合方法为加速单纯形法或麦夸特法。
CN 201110109372 2011-04-28 2011-04-28 一种基于非线性拟合的微震源或声发射源的定位方法 Active CN102262220B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110109372 CN102262220B (zh) 2011-04-28 2011-04-28 一种基于非线性拟合的微震源或声发射源的定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110109372 CN102262220B (zh) 2011-04-28 2011-04-28 一种基于非线性拟合的微震源或声发射源的定位方法

Publications (2)

Publication Number Publication Date
CN102262220A true CN102262220A (zh) 2011-11-30
CN102262220B CN102262220B (zh) 2013-07-17

Family

ID=45008923

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110109372 Active CN102262220B (zh) 2011-04-28 2011-04-28 一种基于非线性拟合的微震源或声发射源的定位方法

Country Status (1)

Country Link
CN (1) CN102262220B (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103206244A (zh) * 2012-01-12 2013-07-17 李成武 以太网总线式煤与瓦斯突出微震声发射集成监测系统
CN103389511A (zh) * 2012-05-10 2013-11-13 韩国地质资源研究院 利用声发射传感器的地层插入型力学行为探测装置及其制作方法
CN104166159A (zh) * 2014-07-15 2014-11-26 刘改成 四维微地震监测的裂缝形态处理方法和系统
CN104502964A (zh) * 2014-12-19 2015-04-08 桂林电子科技大学 一种基于空间几何关系的获得微震波速的方法
CN105022091A (zh) * 2015-08-07 2015-11-04 中国矿业大学 一种无预测速的远场震源快速定位方法
CN105022031A (zh) * 2015-07-03 2015-11-04 四川大学 一种区域岩体微震震源的分层速度定位方法
CN105093314A (zh) * 2015-07-10 2015-11-25 中联煤层气有限责任公司 一种测定微地震震源的方法
CN105866252A (zh) * 2016-06-07 2016-08-17 北京航空航天大学 一种用于中小型长方状箱体声发射源定位的方法
CN105929444A (zh) * 2016-04-08 2016-09-07 中国科学院地质与地球物理研究所 一种基于互相关偏移与最小二乘思想的微地震定位方法
CN106646378A (zh) * 2017-01-03 2017-05-10 中北大学 一种地下挖掘位置的定位方法
CN107727744A (zh) * 2017-10-19 2018-02-23 中南大学 用于岩石力学三轴试验的声发射源定位方法及系统
CN109239775A (zh) * 2018-09-18 2019-01-18 长沙迪迈数码科技股份有限公司 矿产资源被盗采追踪定位方法
CN110376290A (zh) * 2019-07-19 2019-10-25 中南大学 基于多维核密度估计的声发射源定位方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1999042857A1 (en) * 1998-02-23 1999-08-26 Exxonmobil Upstream Research Company Single well system for mapping sources of acoustic energy
CN101482623A (zh) * 2009-02-10 2009-07-15 中国安全生产科学研究院 矿山震动中震源的定位方法
CN101504455A (zh) * 2009-03-06 2009-08-12 沈阳化工学院 一种发射源的时差定位方法
CN101770038A (zh) * 2010-01-22 2010-07-07 中国科学院武汉岩土力学研究所 矿山微震源智能定位方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1999042857A1 (en) * 1998-02-23 1999-08-26 Exxonmobil Upstream Research Company Single well system for mapping sources of acoustic energy
CN101482623A (zh) * 2009-02-10 2009-07-15 中国安全生产科学研究院 矿山震动中震源的定位方法
CN101504455A (zh) * 2009-03-06 2009-08-12 沈阳化工学院 一种发射源的时差定位方法
CN101770038A (zh) * 2010-01-22 2010-07-07 中国科学院武汉岩土力学研究所 矿山微震源智能定位方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
唐绍辉等: "深井矿山地压灾害微震监测技术应用研究", 《岩石力学与工程学报》 *
程慧高: "基于三分量P S 波形的岩体声发射源定位算法", 《工业安全与环保》 *

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103206244A (zh) * 2012-01-12 2013-07-17 李成武 以太网总线式煤与瓦斯突出微震声发射集成监测系统
CN103389511A (zh) * 2012-05-10 2013-11-13 韩国地质资源研究院 利用声发射传感器的地层插入型力学行为探测装置及其制作方法
CN104166159A (zh) * 2014-07-15 2014-11-26 刘改成 四维微地震监测的裂缝形态处理方法和系统
CN104166159B (zh) * 2014-07-15 2015-10-21 刘改成 四维微地震监测的裂缝形态处理方法和系统
CN104502964A (zh) * 2014-12-19 2015-04-08 桂林电子科技大学 一种基于空间几何关系的获得微震波速的方法
CN105022031A (zh) * 2015-07-03 2015-11-04 四川大学 一种区域岩体微震震源的分层速度定位方法
CN105093314B (zh) * 2015-07-10 2017-09-22 中联煤层气有限责任公司 一种测定微地震震源的方法
CN105093314A (zh) * 2015-07-10 2015-11-25 中联煤层气有限责任公司 一种测定微地震震源的方法
CN105022091A (zh) * 2015-08-07 2015-11-04 中国矿业大学 一种无预测速的远场震源快速定位方法
CN105929444B (zh) * 2016-04-08 2018-06-22 中国科学院地质与地球物理研究所 一种基于互相关偏移与最小二乘思想的微地震定位方法
CN105929444A (zh) * 2016-04-08 2016-09-07 中国科学院地质与地球物理研究所 一种基于互相关偏移与最小二乘思想的微地震定位方法
CN105866252A (zh) * 2016-06-07 2016-08-17 北京航空航天大学 一种用于中小型长方状箱体声发射源定位的方法
CN105866252B (zh) * 2016-06-07 2018-09-21 北京航空航天大学 一种用于中小型长方状箱体声发射源定位的方法
CN106646378A (zh) * 2017-01-03 2017-05-10 中北大学 一种地下挖掘位置的定位方法
CN107727744A (zh) * 2017-10-19 2018-02-23 中南大学 用于岩石力学三轴试验的声发射源定位方法及系统
CN109239775A (zh) * 2018-09-18 2019-01-18 长沙迪迈数码科技股份有限公司 矿产资源被盗采追踪定位方法
CN109239775B (zh) * 2018-09-18 2020-08-14 长沙迪迈数码科技股份有限公司 矿产资源被盗采追踪定位方法
CN110376290A (zh) * 2019-07-19 2019-10-25 中南大学 基于多维核密度估计的声发射源定位方法

Also Published As

Publication number Publication date
CN102262220B (zh) 2013-07-17

Similar Documents

Publication Publication Date Title
CN102262220B (zh) 一种基于非线性拟合的微震源或声发射源的定位方法
CN102129063B (zh) 一种微震源或声发射源的定位方法
CN102435980B (zh) 一种基于解析求解的声发射源或微震源定位方法
CN103134473B (zh) 无线气压测高装置及测高方法
CN105022031B (zh) 一种区域岩体微震震源的分层速度定位方法
CN104656123A (zh) 一种测定区域岩体等效波速的方法
CN109061723B (zh) 一种隧洞岩爆孕育过程的微震源高精度定位方法及系统
CN102213769A (zh) 一种利用三维垂直地震剖面资料确定各向异性参数的方法
US20110095763A1 (en) Geologic mapping using lightning
CN105445756A (zh) 构造物的安全性诊断系统
CN104749630B (zh) 构建微地震监测速度模型的方法
CN106949844A (zh) 一种井筒井壁变形自动测量仪及其工作方法
Moriya et al. Delineation of large localized damage structures forming ahead of an active mining front by using advanced acoustic emission mapping techniques
CN105022091B (zh) 一种无预测速的远场震源快速定位方法
CN104237892A (zh) 一种风廓线雷达大气温度探测系统及方法
CN104101648A (zh) 基于李雅普诺夫指数的超声导波定位缺陷的方法
CN103733089B (zh) 用于包括不确定性估计的地下表征的系统和方法
CN103376443B (zh) 探地雷达探测地面干扰快速消除方法
RU2300122C1 (ru) Способ дистанционного определения параметров инфразвукового сигнала вблизи неопознанного источника сигнала
CN109597125B (zh) 一种基于p波到时与最大振幅波形的微震源定位方法
CN110954033A (zh) 混凝土裂缝深度检测方法及其系统
CN104635246A (zh) 卫星导航信号的动态范围检测系统及检测方法
CN104406681B (zh) 一种实时确定微震波速的测试方法
CN102939548A (zh) 确定放置于海底的探测器的位置的方法
FR2941782A1 (fr) Procede de localisation d'apparition d'un defaut dans un milieu par une onde.

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model