CN113496006B - 一种冲击地压矿井矿震空间集中程度的计算方法 - Google Patents
一种冲击地压矿井矿震空间集中程度的计算方法 Download PDFInfo
- Publication number
- CN113496006B CN113496006B CN202110878942.XA CN202110878942A CN113496006B CN 113496006 B CN113496006 B CN 113496006B CN 202110878942 A CN202110878942 A CN 202110878942A CN 113496006 B CN113496006 B CN 113496006B
- Authority
- CN
- China
- Prior art keywords
- events
- mine
- earthquake
- subset
- ore
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种冲击地压矿井矿震空间集中程度的计算方法,首先初步筛选工作面发生的矿震事件;设S为筛选后工作面矿震事件的集合,将矿震事件划分为j个子集,各子集中分别包含n个矿震事件;选择子集中任意两个矿震事件求解欧式距离;对子集中的矿震事件两两之间的欧式距离求和后并取平均值,得到矿震事件平均距离D;子集内矿震事件个数n与矿震事件平均距离D之比即可求得所有子集内矿震空间集中程度值P;通过比较P值,即可判断出工作面中矿震空间集中程度最大的区域。该方法原理简单,易于编程实现,便于现场推广应用。
Description
技术领域
本发明属于煤岩动力灾害防控领域,具体涉及一种冲击地压矿井矿震空间集中程度的计算方法。
背景技术
矿震即为矿山地震,是井下煤岩体在采掘扰动下释放弹性能产生震动的一种现象。随着我国煤矿采深逐年增加,冲击地压灾害愈演愈烈,同时矿震现象也十分普遍。目前微震监测系统是定位矿震最有效的手段,经标波定位可以确定出矿震发生的时间及空间坐标。已有研究表明,矿震的发生与煤岩体的应力集中相关,矿震发生密集,表明该区域煤岩体内应力集中程度高,发生大尺度破坏的可能性大,易发生具有诱冲能力的高能级矿震或者直接引起冲击地压。因此,计算矿震空间集中程度对预测高能矿震发生区域和冲击风险区域具有重要意义。
目前已有的矿震空间集中程度计算方法原理复杂,程序实现难度大,在煤矿现场技术人员水平较低的情况下不易推广应用,阻碍了煤矿现场人员对矿震数据的深度分析与挖掘,不利于开展高能级矿震的预测与防范工作。
发明内容
本发明的目的在于提供一种冲击地压矿井矿震空间集中程度的计算方法,该方法原理简单,易于编程实现,便于现场推广应用。
为了实现上述目的,本发明提供了一种冲击地压矿井矿震空间集中程度的计算方法,该方法基于欧式距离计算原理,包括如下步骤:
(1)初步筛选工作面发生的矿震事件;
(2)设S为筛选后工作面矿震事件的集合,将矿震事件划分为j个子集,各子集可依次命名为S1、S2、S3、……、Sj-2、Sj-1、Sj,即S={S1,S2,S3,S4,……,Sj-2,Sj-1,Sj};
各子集中分别包含n个矿震事件,每个矿震事件可表示为tn,即S1={t1,t2,t3,t4,t5,……,tn-2,tn-1,tn};每个矿震事件的空间位置在笛卡尔直角坐标系中表示为(xn,yn,zn);
(3)选择子集中任意两个矿震事件求解欧式距离,tn-1、tn的欧式距离表示为d(n-1)-n=,tn-2、tn的欧式距离表示为按以下顺序可求取子集中任意两个矿震事件的欧式距离:d1-2、d1-3、d1-4、……、d1-n,d2-3、d2-4、d2-5、……、d2-n,d3-4、d3-5、d3-6、……、d3-n,……,d(n-2)-(n-1)、d(n-2)-n,d(n-1)-n;
(4)对子集中的矿震事件两两之间的欧式距离求和:对d1-2、d1-3、d1-4、……、d1-n累计求和,表示为对d2-3、d2-4、d2-5、……、d2-n累计求和,表示为/>对d3-4、d3-5、d3-6、……、d3-n累计求和,表示为/>依次类推,直至d(n-1)-n;
(5)对子集中的矿震事件两两之间的欧式距离取平均值,得到矿震事件平均距离并设为D,计算方法如下:
;
(6)子集内矿震事件个数n与矿震事件平均距离D之比即可求得所有子集内矿震空间集中程度值,即P1、P2、P3、....、Pj,即:
进一步的,步骤(1)中,筛选准则为:矿震事件垂直标高位于开采煤层顶板上方20m至开采煤层底板下方20m范围内。
进一步的,步骤(2)中,将矿震事件按空间坐标或者时间顺序划分为j个子集。
本发明基于欧式距离提出了一种矿震空间集中程度的计算方法,原理简单、可操作性强,方法中涉及的计算过程易于用编程实现,便于现场推广使用;本方法可提高冲击地压矿井对矿震空间集中程度的量化评估水平,对提升冲击地压矿井的监测预警水平具有重要的作用。
附图说明
图1是实施例一工作面按垂直标高筛选后全部的矿震事件空间投影图;
图2是实施例一矿震事件集合按空间区域划分示意图;
图3是实施例一矿震事件中某个子集的欧式距离表示方式示意图;
图4是实施例一工作面各空间区域矿震空间集中程度值P的统计图;
图5是实施例二工作面按垂直标高筛选后5月份的矿震事件空间投影图;
图6是实施例二5月21日至5月26日的矿震事件分布图;(A)-(F)分别为5月21日、5月22日、5月23日、5月24日、5月25日、5月26日当日矿震事件分布图;
图7是实施例二工作面每日矿震空间集中程度值P的统计图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明。
实施例一
一种冲击地压矿井矿震空间集中程度的计算方法,该方法基于欧式距离计算原理,包括如下步骤:
(1)如图1所示为本实施例初步筛选工作面发生的矿震事件,对矿震事件进行震源空间集中程度分析,工作面煤层底板标高为+130m,煤层顶板标高为+140m,根据矿震事件的垂直坐标进行初步筛选,则保留垂直坐标处于+110m~+160m的矿震事件,剔除其余垂直标高距离开采煤层较远的矿震事件;
(2)设S为筛选后工作面矿震事件的集合,集合S共有7066个矿震事件;如图2所示,在xy平面将矿震事件沿工作面走向即x轴方向用矩形划分,矩形宽度30m,矩形长度可将顺槽两侧的矿震事件包含即可,工作面走向长度即总宽度为1050m,则该工作面可用35个矩形分割;将矿震事件按空间坐标划分为35个子集,各子集可依次命名为S1、S2、S3、……、S33、S34、S35,即S={S1,S2,S3,……,S33,S34,S35};
各子集中分别包含n个矿震事件,每个矿震事件可表示为tn,即S1={t1,t2,t3,t4,t5,……,tn-2,tn-1,tn};每个矿震事件的空间位置在笛卡尔直角坐标系中表示为(xn,yn,zn);
例如,如图2所示,在笛卡尔直角坐标系中,S1为空间坐标(41~71,23~323,110~160)所包围的柱体空间,则判断矿震事件tn是否处于S1空间区域的准则为41≤xn<71且23≤yn≤323且110≤zn≤160;S2为空间坐标(71~101,23~323,110~160)所包围的柱体空间,则判断矿震事件tn是否处于S2的准则为71≤xn<101且23≤yn≤323且110≤zn≤160;.........;S35为空间坐标(1020~1050,23~323,110~160)所包围的柱体空间,则判断矿震事件tn是否处于S35的准则为1020≤xn<1050且23≤yn≤323且110≤zn≤160,从将矿震事件集合S划分到35个矩形区域内,具体划分结果为:S1~S35矿震事件数目n分别为:5、18、10、16、33、30、34、37、63、36、51、58、68、66、79、92、128、216、258、326、385、518、505、428、331、255、234、212、264、278、360、400、490、601、181;
(3)如图3所示,选择S1子集中任意两个矿震事件求解欧式距离,t1、t2的欧式距离表示为t1、t3的欧式距离表示为按以下顺序可求取子集中任意两个矿震事件的欧式距离:d1-2、d1-3、d1-4、d1-5,d2-3、d2-4、d2-5,d3-4、d3-5,d4-5;
(4)对子集中的矿震事件两两之间的欧式距离求和:对d1-2、d1-3、d1-4、d1-5累计求和,表示为对d2-3、d2-4、d2-5累计求和,表示为/>对d3-4、d3-5累计求和,表示为/>依次类推,直至d4-5;
(5)对S1子集中的矿震事件两两之间的欧式距离取平均值,得到矿震事件平均距离并设为D1,计算方法如下:
同上述步骤(3)至步骤(5)的原理和过程,对S2至S35子集中任意两个矿震事件求解欧式距离、再对欧式距离进行求和、最后计算出矿震事件平均距离,即为D2、D3、…、D33、D34、D35;
S1~S35矿震事件平均距离Dj分别为:106.49、87.99、116.51、74.31、99.62、97.62、63.44、119.92、112.14、129.56、112.64、120.71、114.99、102.29、106.73、92.91、90.22、89.44、80.29、79.89、77.78、76.54、74.77、74.47、77.67、80.59、82.15、82.66、76.71、74.38、76.28、66.38、73.7、76.9、73.63;
(6)子集内矿震事件个数n与矿震事件平均距离D之比即可求得所有子集内矿震空间集中程度值,即P1、P2、P3、....、P35,即:
计算结果如图4所示,S1~S35区域的Pj分别为:0.05、0.2、0.09、0.22、0.33、0.31、0.54、0.31、0.56、0.28、0.45、0.48、0.59、0.65、0.74、0.99、1.42、2.41、3.21、4.08、4.95、6.77、6.75、5.75、4.26、3.16、2.85、2.56、3.44、3.74、4.72、6.03、6.65、7.81、2.46;
如图4所示,通过比较P1、P2、P3、....、P35,即可判断出工作面中矿震空间集中程度最大的区域为S22,该区域P值为6.77。
实施例二
一种冲击地压矿井矿震空间集中程度的计算方法,该方法基于欧式距离计算原理,包括如下步骤:
(1)如图5所示为本实施例初步筛选工作面发生的矿震事件,包含了该工作面从5月1日至5月31日的全部矿震事件,对矿震事件进行震源空间集中程度分析,工作面煤层底板标高为+130m,煤层顶板标高为+140m,根据矿震事件的垂直坐标进行初步筛选,则保留垂直坐标处于+110m~+160m的矿震事件,剔除其余垂直标高距离开采煤层较远的矿震事件;
(2)设S为筛选后工作面的5月份全部矿震事件的集合,集合S共有1623个矿震事件;如图6所示为5月21日、5月22日、5月23日、5月24日、5月25日、5月26日当日矿震事件分布图;将矿震事件按时间顺序划分为31个子集,各子集可依次命名为S1、S2、S3、……、S29、S30、S31,即S={S1,S2,S3,……,S29、S30、S31};
各子集中分别包含n个矿震事件,每个矿震事件可表示为tn,即S1={t1,t2,t3,t4,t5,……,tn-2,tn-1,tn};每个矿震事件的空间位置在笛卡尔直角坐标系中表示为(xn,yn,zn);
按时间顺序发生矿震事件统计结果为:5月1日至5月31日每个子集合Sj的矿震数目n分别为:16、15、18、28、19、21、22、31、34、5、25、53、34、33、35、39、72、42、79、103、61、59、92、73、67、51、78、72、156、56、134;
(3)选择S1子集中任意两个矿震事件求解欧式距离,t1、t2的欧式距离表示为t1、t3的欧式距离表示为按以下顺序可求取子集中任意两个矿震事件的欧式距离:d1-2、d1-3、d1-4、……、d1-16,d2-3、d2-4、d2-5、……、d2-16,d3-4、d3-5、d3-6、……、d3-16,……,d14-15、d14-16,d15-16;
(4)对子集中的矿震事件两两之间的欧式距离求和:对d1-2、d1-3、d1-4、……、d1-16累计求和,表示为对d2-3、d2-4、d2-5、……、d2-16累计求和,表示为/>对d3-4、d3-5、d3-6、……、d3-16累计求和,表示为/>依次类推,直至d15-16;
(5)对S1子集中的矿震事件两两之间的欧式距离取平均值,得到矿震事件平均距离并设为D,计算方法如下:
同上述步骤(3)至步骤(5)的原理和过程,对S2至S31子集中任意两个矿震事件求解欧式距离、再对欧式距离进行求和、最后计算出矿震事件平均距离,即为D2、D3、…、D29、D30、D31;
S1~S31矿震事件平均距离Dj分别为:151.23、173.73、142.73、155.5、144.06、188.07、150.94、96.19、99.2、180.6、174.76、189.38、186.51、142.16、108.74、106.01、109.03、112.53、127.4、151.78、128.38、120.06、115.77、123.93、177.1、135.21、157.59、114.56、116.68、101.9、102.95;
(6)子集内矿震事件个数n与矿震事件平均距离D之比即可求得所有子集内矿震空间集中程度值,即P1、P2、P3、....、P31,即:
计算结果如图7所示,S1~S31区域的Pj分别为:0.11、0.09、0.13、0.18、0.13、0.11、0.15、0.32、0.34、0.03、0.14、0.28、0.18、0.23、0.32、0.37、0.66、0.37、0.62、0.68、0.48、0.49、0.79、0.59、0.38、0.38、0.49、0.63、1.34、0.55、1.30;
通过比较P1、P2、P3、....、P31,即可判断出工作面中矿震空间集中程度最大的区域为S29,5月份该工作面矿震空间集中程度最大的时间为5月29日,当日P值为1.34。
Claims (3)
1.一种冲击地压矿井矿震空间集中程度的计算方法,其特征在于,该方法基于欧式距离计算原理,包括如下步骤:
(1)初步筛选工作面发生的矿震事件;
(2)设S为筛选后工作面矿震事件的集合,将矿震事件划分为j个子集,各子集依次命名为S1、S2、S3、……、Sj-2、Sj-1、Sj,即S={S1,S2,S3,S4,……,Sj-2,Sj-1,Sj};
各子集中分别包含n个矿震事件,每个矿震事件表示为tn,即S1={t1,t2,t3,t4,t5,……,tn-2,tn-1,tn};每个矿震事件的空间位置在笛卡尔直角坐标系中表示为(xn,yn,zn);
(3)选择子集中任意两个矿震事件求解欧式距离,tn-1、tn的欧式距离表示为 tn-2、tn的欧式距离表示为/> 按以下顺序求取子集中任意两个矿震事件的欧式距离:d1-2、d1-3、d1-4、……、d1-n,d2-3、d2-4、d2-5、……、d2-n,d3-4、d3-5、d3-6、……、d3-n,……,d(n-2)-(n-1)、d(n-2)-n,d(n-1)-n;
(4)对子集中的矿震事件两两之间的欧式距离求和:对d1-2、d1-3、d1-4、……、d1-n累计求和,表示为对d2-3、d2-4、d2-5、……、d2-n累计求和,表示为/>对d3-4、d3-5、d3-6、……、d3-n累计求和,表示为/>依次类推,直至d(n-1)-n;
(5)对子集中的矿震事件两两之间的欧式距离取平均值,得到矿震事件平均距离并设为D,计算方法如下:
;
(6)子集内矿震事件个数n与矿震事件平均距离D之比即求得所有子集内矿震空间集中程度值,即P1、P2、P3、…、Pj,即:
2.根据权利要求1所述的一种冲击地压矿井矿震空间集中程度的计算方法,其特征在于,步骤(1)中,筛选准则为:矿震事件垂直标高位于开采煤层顶板上方20m至开采煤层底板下方20m范围内。
3.根据权利要求1或2所述的一种冲击地压矿井矿震空间集中程度的计算方法,其特征在于,步骤(2)中,将矿震事件按空间坐标或者时间顺序划分为j个子集。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110878942.XA CN113496006B (zh) | 2021-08-02 | 2021-08-02 | 一种冲击地压矿井矿震空间集中程度的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110878942.XA CN113496006B (zh) | 2021-08-02 | 2021-08-02 | 一种冲击地压矿井矿震空间集中程度的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113496006A CN113496006A (zh) | 2021-10-12 |
CN113496006B true CN113496006B (zh) | 2023-08-01 |
Family
ID=77995904
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110878942.XA Active CN113496006B (zh) | 2021-08-02 | 2021-08-02 | 一种冲击地压矿井矿震空间集中程度的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113496006B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20100117344A (ko) * | 2009-04-24 | 2010-11-03 | 삼성전자주식회사 | 단일 반송파 주파수 분할 다중 접속 시스템을 위한 수신 장치 및 방법 |
CN202939619U (zh) * | 2012-10-31 | 2013-05-15 | 中国矿业大学(北京) | 一种基于图像离散多小波变换的煤岩识别系统 |
CN110761841A (zh) * | 2019-10-28 | 2020-02-07 | 煤炭科学技术研究院有限公司 | 一种基于微震监测的矿井群工作面开采互扰范围计算方法 |
CN112377258A (zh) * | 2020-11-04 | 2021-02-19 | 中国矿业大学 | 一种临近采区同采工作面间采动影响程度的评价方法 |
-
2021
- 2021-08-02 CN CN202110878942.XA patent/CN113496006B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20100117344A (ko) * | 2009-04-24 | 2010-11-03 | 삼성전자주식회사 | 단일 반송파 주파수 분할 다중 접속 시스템을 위한 수신 장치 및 방법 |
CN202939619U (zh) * | 2012-10-31 | 2013-05-15 | 中国矿业大学(北京) | 一种基于图像离散多小波变换的煤岩识别系统 |
CN110761841A (zh) * | 2019-10-28 | 2020-02-07 | 煤炭科学技术研究院有限公司 | 一种基于微震监测的矿井群工作面开采互扰范围计算方法 |
CN112377258A (zh) * | 2020-11-04 | 2021-02-19 | 中国矿业大学 | 一种临近采区同采工作面间采动影响程度的评价方法 |
Non-Patent Citations (2)
Title |
---|
东滩煤矿地垒断层型矿震发生机制及监测分析;王世娟;第05卷(第《中国优秀硕士学位论文全文数据库 工程科技I辑》期);B021-33 * |
煤矿冲击矿压前兆信息识别及预警模型研究;韩泽鹏;《中国优秀硕士学位论文全文数据库 工程科技I辑》(第01期);B021-278 * |
Also Published As
Publication number | Publication date |
---|---|
CN113496006A (zh) | 2021-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cheng et al. | A zoning model for coal mining-induced strata movement based on microseismic monitoring | |
CN109798149B (zh) | 一种坚硬顶板厚煤层采场来压强度分级预测方法 | |
Li et al. | Rockburst occurrences and microseismicity in a longwall panel experiencing frequent rockbursts | |
CN113514886B (zh) | 一种砂岩型铀矿成矿有利部位地质-地震三维预测方法 | |
CN110135030A (zh) | 一种采空区地表沉降的预测方法 | |
CN111581834B (zh) | 一种基于多源信息融合的煤层顶板溃水溃沙危险性评价方法 | |
CN111222254A (zh) | 基于应力叠加法的工作面冲击地压危险等级划分方法和系统 | |
CN116050171B (zh) | 一种掘进巷道人造解放层防冲方法 | |
Ma et al. | Fracture mechanism of roof key strata in Dongjiahe coal mine using microseismic moment tensor | |
Ma et al. | Microseismic monitoring, analysis and early warning of rockburst | |
CN104564069A (zh) | 一种基于方格网法的地面动态沉陷预测与复垦方法 | |
CN113496006B (zh) | 一种冲击地压矿井矿震空间集中程度的计算方法 | |
Cieślik et al. | Geomechanical analysis of location and conditions for mining-induced tremors in LGOM copper mines | |
CN115045660B (zh) | 一种大倾角回采工作面顶板深孔爆破致裂方法及系统 | |
CN112377258A (zh) | 一种临近采区同采工作面间采动影响程度的评价方法 | |
Gai et al. | A new method for evaluating floor spatial failure characteristics and water inrush risk based on microseismic monitoring | |
Oparin et al. | On possible causes of increase in seismic activity of mine fields in the Oktyabrsky and Taimyrsky mines of the Norilsk deposit in 2003. Part I: Seismic regime | |
CN113434819B (zh) | 工作面采动对采空区矿震活动影响时间和距离的确定方法 | |
CN112379460B (zh) | 基于应力场和震动波场的合理停采线位置确定的方法 | |
Westley-Hauta et al. | Characterisation of seismic activity at a kimberlite block caving operation in a complex geological setting in Quebec, Canada | |
CN115263427A (zh) | 覆岩破坏高度确定方法、装置、电子设备及存储介质 | |
Xue et al. | Mechanism and prevention practice of coal burst in deep gradual residual coal pillar area: a case study | |
CN114063152A (zh) | 一种基于矿震统计特征的冲击地压主控因素确定方法 | |
Ma et al. | Monitoring and evaluation of disaster risk caused by linkage failure and instability of residual coal pillar and rock strata in multi-coal seam mining | |
Heal et al. | Seismic hazard at Agnico-Eagle's Laronde mine using MS-RAP |
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 |