CN101482623A - 矿山震动中震源的定位方法 - Google Patents
矿山震动中震源的定位方法 Download PDFInfo
- Publication number
- CN101482623A CN101482623A CNA2009100776492A CN200910077649A CN101482623A CN 101482623 A CN101482623 A CN 101482623A CN A2009100776492 A CNA2009100776492 A CN A2009100776492A CN 200910077649 A CN200910077649 A CN 200910077649A CN 101482623 A CN101482623 A CN 101482623A
- Authority
- CN
- China
- Prior art keywords
- focus
- coordinate
- equations
- sensor
- sensors
- 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.)
- Pending
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供了一种矿山震动中震源的定位方法,该方法属于采矿领域,该方法包括:在矿山上至少4个不同位置设置传感器,并获取所述至少4个传感器的坐标和所述至少4个传感器检测出的ti、t0;通过获取的所述至少4个传感器的坐标和至少4个传感器检测出的ti、t0计算得出震源的坐标。本发明具体实施方式提供的矿山震动中震源的定位方法具有计算过程简单,计算量小,适用矿山震源的快速定位的特点。
Description
技术领域
本发明涉及采矿领域,尤其涉及一种矿山震动中震源的定位方法。
背景技术
采矿技术为国家大力发展的技术,采矿技术的高低与国民经济的高低有直接的关系。在采矿过程中通常伴随着,如矿震引起的震动或人工放炮引起的震动等,而采矿过程中的震动容易造成各种矿山事故,对生命和财产安全造成重大威胁,所以对采矿过程中的震动研究有利于维护生命和财产安全,能够最大限度的减少损失。
在矿山震动研究中,震源的定位是极为重要的一个环节,现有技术通常采用的穷举法或双差法等定位方式进行震源的定位。穷举法计算每个可能震中位置的目标函数,取目标函数最小值为震中,并以此分层计算震源深度及发震时刻;双差法引入“事件对”i,j及双重残差,对所有台站和事件对反演得到震源的绝对位置。
在实现本发明的过程中,发明人发现现有技术存在如下问题:
现有技术的方案是针对地震等远距离震源采用的穷举法或双差法,进行震源定位的,所以其计算过程复杂,计算量大,并不适用于矿山震源的快速定位。
发明内容
本发明的具体实施方式提供一种矿山震动中震源的定位方法,该方法具有计算过程简单,计算量小,适用矿山震源的快速定位的特点。
本发明具体实施方式提供一种矿山震动中震源的定位方法,该方法包括:
A、在矿山上至少4个不同位置设置传感器,并获取所述至少4个传感器的坐标和所述至少4个传感器检测出的ti、t0;其中ti为震源到第i个传感器时间,t0为测量起始时间;
B、通过获取的所述至少4个传感器的坐标和至少4个传感器检测出的ti、t0计算得出震源的坐标;计算方法具体为:
将获取的所述至少4个传感器的坐标和至少4个传感器检测出的ti、t0分别代入vi 2(ti-t0)2=(x-ai)2+(y-bi)2+(z-ci)2获得包括至少4个方程的方程组,将所述方程组求解得到震源的坐标;其中vi为震源的传播速度;x,y,z为所述震源的坐标;ai,bi,Ci为第i个传感器的坐标。
由上述本发明的具体实施例提供的技术方案可以看出,本发明的具体实施例所述的技术方案通过4个传感器的坐标和至少4个传感器检测出的ti、t0分别代入vi 2(ti-t0)2=(x-ai)2+(y-bi)2+(z-ci)2获得包括至少4个方程的方程组,将所述方程组求解得到震源的坐标,其计算方法简单,计算量小。
附图说明
图1为本发明具体实施方式提供的矿山震动中震源的定位方法的流程图。
具体实施方式
本发明具体实施方式提供一种矿山震动中震源的定位方法,该方法的原理为,通过定出震源空间位置的3个坐标值来确定震源的位置,解决这一问题的最简捷方法就是用标准的距离方程:
式(1)中:ai,bi,ci为第i个传感器的坐标;x,y,z为所求震源的坐标;di为震源至第i个传感器的距离。
设Vi为震源至第i个传感器的传播速度,ti为震源到第i个传感器时间,t0为测量起始时间,则有:di=vi(ti-t0) (2)
(1)式可以转化为:
vi 2(ti-t0)2=(x-ai)2+(y-bi)2+(z-ci)2(3)
为了求解出式(3)中的x,y,z,vi,则需要4个独立的方程组成方程组,这里假设震源至每个传感器的传播速度Vi均相同;则只需要4个不同位置的传感器的坐标以及该传感器检测出的ti、t0组成方程组即可求出震源的坐标。
通过上述说明,本发明具体实施方式提供的一种矿山震动中震源的定位方法如图1所示,包括如下步骤:
步骤11、在矿山上至少4个不同位置设置传感器,并获取至少4个传感器的坐标和至少4个传感器检测出的ti、t0;其中ti为震源到第i个传感器时间,t0为测量起始时间;
步骤12,通过获取的至少4个传感器的坐标和至少4个传感器检测出的ti、t0计算得出震源的坐标。
实现该步骤12的具体方法可以为,将至少4个传感器的坐标和至少4个传感器检测出的ti、t0分别代入上式(3)中,得到包括至少4个独立方程组成的方程组,将该方程组求解即得到震源的坐标。
可选的,由于在实际情况中,将至少4个传感器的坐标和至少4个传感器检测出的ti、t0分别代入上式(3)中得到的方程组为非线性方程,而求解该非线性方程组成的方程组十分的复杂,为了简化求解过程,在实际情况中,可以将传感器增加到至少5个,这样就可以得到至少5个独立的非线性方程组成的方程组,由于该5个独立的的非线性方程均具有相同的未知数,所以将该5个独立的的非线性方程组成的方程组简化成由4个的独立的线性方程组成的方程组,其具体的简化方法可以为,这里以5个方程为例来进行说明,假设这里5个独立的非线性方程分别为,方程1、方程2、方程3、方程4、方程5;方程1—方程2得到第一个线性方程,方程2—方程3得到第二个线性方程,方程3—方程4得到第三个线性方程,方程4—方程5得到第四个线性方程,当然在实际情况中,还可以有其他的得到线性方程的方法,如由方程1分别减去方程2、3、4、5也可以得到4个线性方程;本发明具体实施方式并不限制该简化的具体方法。将得到的4个线性方程组合起来即得到由4个线性方程组合得到的方程组,对该方程组求解即可得到震源的坐标。
本发明具体实施方式提供的一种矿山震动中震源的定位方法只需通过获取至少4个点的传感器的坐标以及传感器的检测数据计算出震源的位置,所以该方法计算方法简单,计算量小,适合矿山震源的快速定位。
本发明具体实施方式所述的技术方案具有简单,计算量小,适合矿山震源的快速定位的优点。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (2)
1、一种矿山震动中震源的定位方法,其特征在于,所述方法包括:
A、在矿山上至少4个不同位置设置传感器,并获取所述至少4个传感器的坐标和所述至少4个传感器检测出的ti、t0;其中ti为震源到第i个传感器时间,t0为测量起始时间;
B、通过获取的所述至少4个传感器的坐标和至少4个传感器检测出的ti、t0计算得出震源的坐标;计算方法具体为:
将获取的所述至少4个传感器的坐标和至少4个传感器检测出的ti、t0分别代入vi 2(ti-t0)2=(x-ai)2+(y-bi)2+(z-ci)2获得包括至少4个方程的方程组,将所述方程组求解得到震源的坐标;其中Vi为震源的传播速度;x,y,z为所述震源的坐标;ai,bi,ci为第i个传感器的坐标。
2、根据权利要求1所述的方法,其特征在于,所述至少4个传感器为大于或等于5个,当所述传感器大于或等于5个时,所述步骤B具体包括:
将获取的至少5个传感器的坐标和至少5个传感器检测出的ti、t0分别代入vi 2(ti-t0)2=(x-ai)2+(y-bi)2+(z-ci)2获得包括至少5个方程的方程组,并将包括至少5个方程的方程组简化成包括至少4个线性方程的方程组后,求解包括至少4个线性方程的方程组得到震源的坐标;其中Vi为震源的传播速度;x,y,z为所述震源的坐标;ai,bi,ci为第i个传感器的坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2009100776492A CN101482623A (zh) | 2009-02-10 | 2009-02-10 | 矿山震动中震源的定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2009100776492A CN101482623A (zh) | 2009-02-10 | 2009-02-10 | 矿山震动中震源的定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101482623A true CN101482623A (zh) | 2009-07-15 |
Family
ID=40879823
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2009100776492A Pending CN101482623A (zh) | 2009-02-10 | 2009-02-10 | 矿山震动中震源的定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101482623A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102096093A (zh) * | 2010-11-29 | 2011-06-15 | 杨本才 | 一种利用微震点作为震源计算矿区地震波传播速度的方法 |
CN102262220A (zh) * | 2011-04-28 | 2011-11-30 | 中南大学 | 一种基于非线性拟合的微震源或声发射源的定位方法 |
CN102279410A (zh) * | 2011-06-21 | 2011-12-14 | 北京蓝尊科技有限公司 | 矿山地下开采活动实时监测系统及其方法 |
CN102778668A (zh) * | 2012-07-23 | 2012-11-14 | 中煤科工集团西安研究院 | 矿山被动震源快速精确定位方法 |
CN104360391A (zh) * | 2014-12-02 | 2015-02-18 | 武汉科技大学 | 一种基于聚束阵列波形的微震源定位方法 |
CN104536035A (zh) * | 2015-01-16 | 2015-04-22 | 淮南矿业(集团)有限责任公司 | 煤系地层中震源位置获取方法 |
CN105607040A (zh) * | 2015-09-07 | 2016-05-25 | 中国神华能源股份有限公司 | 一种矿区防盗采监测定位方法及系统 |
CN108931816A (zh) * | 2018-08-17 | 2018-12-04 | 山东省科学院激光研究所 | 一种震源定位方法及装置 |
CN111007564A (zh) * | 2019-12-23 | 2020-04-14 | 辽宁工程技术大学 | 一种基于手机传感器的露天矿山爆破震动监测方法 |
CN111221036A (zh) * | 2020-01-21 | 2020-06-02 | 中南大学 | 一种含未知空洞的目标区域震源定位方法及系统 |
-
2009
- 2009-02-10 CN CNA2009100776492A patent/CN101482623A/zh active Pending
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102096093A (zh) * | 2010-11-29 | 2011-06-15 | 杨本才 | 一种利用微震点作为震源计算矿区地震波传播速度的方法 |
CN102096093B (zh) * | 2010-11-29 | 2013-04-10 | 杨本才 | 一种利用微震点作为震源计算矿区地震波传播速度的方法 |
CN102262220A (zh) * | 2011-04-28 | 2011-11-30 | 中南大学 | 一种基于非线性拟合的微震源或声发射源的定位方法 |
CN102262220B (zh) * | 2011-04-28 | 2013-07-17 | 中南大学 | 一种基于非线性拟合的微震源或声发射源的定位方法 |
CN102279410A (zh) * | 2011-06-21 | 2011-12-14 | 北京蓝尊科技有限公司 | 矿山地下开采活动实时监测系统及其方法 |
CN102778668A (zh) * | 2012-07-23 | 2012-11-14 | 中煤科工集团西安研究院 | 矿山被动震源快速精确定位方法 |
CN104360391A (zh) * | 2014-12-02 | 2015-02-18 | 武汉科技大学 | 一种基于聚束阵列波形的微震源定位方法 |
CN104360391B (zh) * | 2014-12-02 | 2017-02-22 | 武汉科技大学 | 一种基于聚束阵列波形的微震源定位方法 |
CN104536035A (zh) * | 2015-01-16 | 2015-04-22 | 淮南矿业(集团)有限责任公司 | 煤系地层中震源位置获取方法 |
CN104536035B (zh) * | 2015-01-16 | 2017-03-22 | 淮南矿业(集团)有限责任公司 | 煤系地层中震源位置获取方法 |
CN105607040A (zh) * | 2015-09-07 | 2016-05-25 | 中国神华能源股份有限公司 | 一种矿区防盗采监测定位方法及系统 |
CN108931816A (zh) * | 2018-08-17 | 2018-12-04 | 山东省科学院激光研究所 | 一种震源定位方法及装置 |
CN108931816B (zh) * | 2018-08-17 | 2020-01-21 | 山东省科学院激光研究所 | 一种震源定位方法及装置 |
CN111007564A (zh) * | 2019-12-23 | 2020-04-14 | 辽宁工程技术大学 | 一种基于手机传感器的露天矿山爆破震动监测方法 |
CN111221036A (zh) * | 2020-01-21 | 2020-06-02 | 中南大学 | 一种含未知空洞的目标区域震源定位方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101482623A (zh) | 矿山震动中震源的定位方法 | |
CN104536043B (zh) | 一种深度域整体速度模型融合方法及装置 | |
CN104199090B (zh) | 一种地面监测微地震定位的速度模型构建与求解方法 | |
CN103454646B (zh) | 用于处理基于北斗的高精度位移形变监测数据的方法 | |
CN101520516A (zh) | 一种三维地震记录的叠前似二维变换的方法 | |
CN106054214A (zh) | 一种可削弱多路径效应的gnss数据处理方法 | |
CN113064170A (zh) | 一种基于时序InSAR技术的膨胀土区域地表形变监测方法 | |
CN103713322A (zh) | 一种方位高阶相干振幅变化率裂缝检测方法 | |
CN105093319A (zh) | 基于三维地震数据的地面微地震静校正方法 | |
De Guidi et al. | Brief communication: Co-seismic displacement on 26 and 30 October 2016 (Mw= 5.9 and 6.5)–earthquakes in central Italy from the analysis of a local GNSS network | |
CN106767922A (zh) | 一种海洋卫星星载高度计在轨绝对定标系统及定标方法 | |
CN104567802B (zh) | 集成船载重力和gnss的测线式陆海高程传递方法 | |
CN110865338B (zh) | 一种基于水下机动平台的未知周期声信标高精度定位方法 | |
CN111651711A (zh) | 一种地质勘查钻孔地理空间数据坐标转换方法 | |
CN109782355B (zh) | Obs检波点漂移的检测方法及装置 | |
CN104297790A (zh) | 一种起伏海底地震观测系统中炮检距分布的定量评估方法 | |
CN107102342A (zh) | 基于共模误差的gps坐标时间序列不连续性的补足方法 | |
CN111060064B (zh) | 一种惯性导航结合三角定位的智能管道形变检测系统 | |
Bingley et al. | Monitoring the vertical land movement component of changes in mean sea level using GPS: results from tide gauges in the UK | |
CN102353992A (zh) | 基于中心点规则化的起伏地表动校正方法 | |
Yao et al. | Emerged fossil corals on the coast of northwestern Hainan Island, China: Implications for mid-Holocene sea level change and tectonic uplift | |
CN104077476A (zh) | 扫描仪坐标系与工程测量坐标系的变换方法 | |
CN107015251A (zh) | 一种全球导航卫星系统伪距单点定位误差改正的方法 | |
CN104634260A (zh) | 地表位移监测系统 | |
De Guidi et al. | Brief Communication Co-seismic displacement on October 26 and 30, 2016 (M w 5.9 and 6.5)–earthquakes in central Italy from the analysis of discrete GNSS network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Open date: 20090715 |