CN105765408A - 用于分析位于地下矿场巷道上方的层的地质结构以及相对应力变化的方法及系统 - Google Patents
用于分析位于地下矿场巷道上方的层的地质结构以及相对应力变化的方法及系统 Download PDFInfo
- Publication number
- CN105765408A CN105765408A CN201480002443.2A CN201480002443A CN105765408A CN 105765408 A CN105765408 A CN 105765408A CN 201480002443 A CN201480002443 A CN 201480002443A CN 105765408 A CN105765408 A CN 105765408A
- Authority
- CN
- China
- Prior art keywords
- seismic
- mining
- station
- record
- mining site
- 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
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000004458 analytical method Methods 0.000 title abstract description 3
- 230000003245 working effect Effects 0.000 title abstract 3
- 238000005065 mining Methods 0.000 claims abstract description 68
- 238000005259 measurement Methods 0.000 claims abstract description 49
- 239000011435 rock Substances 0.000 claims abstract description 22
- 238000003325 tomography Methods 0.000 claims abstract description 18
- 230000005540 biological transmission Effects 0.000 claims abstract description 11
- 230000011664 signaling Effects 0.000 claims abstract description 7
- 230000008859 change Effects 0.000 claims description 22
- 238000003384 imaging method Methods 0.000 claims description 11
- 230000005534 acoustic noise Effects 0.000 claims description 6
- 238000004587 chromatography analysis Methods 0.000 claims description 5
- 238000012800 visualization Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 238000010835 comparative analysis Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims 1
- 230000002596 correlated effect Effects 0.000 abstract 2
- 230000000875 corresponding effect Effects 0.000 abstract 1
- 238000005305 interferometry Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 239000006185 dispersion Substances 0.000 description 3
- 230000012447 hatching Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000002360 explosive Substances 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000002045 lasting effect Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 241000287219 Serinus canaria Species 0.000 description 1
- 206010044565 Tremor Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000003028 elevating effect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 230000003116 impacting effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
Classifications
-
- 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/303—Analysis for determining velocity profiles or travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/12—Signal generation
- G01V2210/123—Passive source, e.g. microseismics
- G01V2210/1236—Acoustic daylight, e.g. cultural noise
-
- 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/612—Previously recorded data, e.g. time-lapse or 4D
- G01V2210/6122—Tracking reservoir changes over time, e.g. due to production
-
- 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/612—Previously recorded data, e.g. time-lapse or 4D
- G01V2210/6122—Tracking reservoir changes over time, e.g. due to production
- G01V2210/6124—Subsidence, i.e. upwards or downwards
-
- 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/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
在所提出的分析位于地下矿场采矿巷道上方的层的地质结构以及相对应力变化的方法中，将通过在时间上密切相关地记录来自地面系统的低频地震噪声(Dn.cz.)以及由采矿产生的地震突发(Dw.cz.)而获得的来自移动测量数据记录器(3)和来自矿场地震系统的中心站(10)的数据传输到固定处理中心(1)。随后，使用用于所述噪声记录的地震干涉测量方法以及用于所述矿场突裂记录的无源速度和/或衰减层析成像，来处理优选地30秒的时间窗中的呈由采矿产生的低频地震噪声(Dn.cz.)和地震突发(Dw.cz.)的三轴记录形式的所述所记录测量数据。在此基础上，针对岩体(7)的所研究区域以所述无源速度和/或衰减层析成像方法确定横波速度的等值线以及纵波速度的等值线和/或衰减。它们最后将体现采矿巷道(B)上方的所述层中相对应力变化(ΔNP)的平均状态。在矿场突裂(W)的时刻，位置坐标(X、Y和Z)的参数以及计算出的突裂发生时间(To)与以下相关：矿场地面上记录的波的三轴低频测量站(5)的记录中由此岩层突裂产生的纵波首次进入的时间(Tp)，以及从所述纵波的所述首次输入直到所述岩层突裂(W)的所记录信号达到每一测量站(5)中的首次最大值为止所述信号的相对应上升时间。在所提出的测量系统中，所述固定处理中心(1)优选地经由GSM调制解调器(2)连接到所述移动测量数据记录器(3)以及到所述矿场地震系统的中心站(10)上，所述中心站连接到时钟(GPS)以及可视化信令模块(11)上，并且经由本质安全的数字传输系统(12)以及矿场远程传输网络(13)连接到至少四个地下测震站(14)和/或至少四个地下地震听地器站(15)上。
Description
技术领域
本发明涉及一种方法及系统，其用于分析位于地下矿场巷道上方的层的地质结构以及相对应力变化，以既定确定由于在此区域采矿造成岩层突裂的可能性而引起的地面上基础设施的危险。
背景技术
现有技术水平
目前，使用侵入性方法来进行采矿巷道上方岩层的地质结构的确定，所述侵入性方法包括激发由定位在矿场地面上的大功率致动器或爆炸物产生的在岩体中传播的地震波。通常使用地震层析成像工具。
低频无源地震探测(LFS)方法的应用以及与区域性质的噪声相关的现象(通常由例如地震、火山活动，大海或海洋的波浪、海流对大陆影响或气团的快速流动影响等自然过程引起)从以下出版物中已知：当热尔S.(DangelS.)的“油气层上观察到的类震颤信号的现象论(Phenomenologyoftremor-likesignalsobservedoverhydrocarbonreservoirs)”(火山学和地热学研究杂志(JournalofVolcanologyandGeothermalResearch)，2003年，128(1-3)，第135到158页)；乔巴帝科夫A.V.(GorbatikovA.V.)、加里宁娜A.V.(KalininaA.V.)、沃尔科夫V.A.(VolkovV.A.)等人的“在兰萨罗特岛、加那利、西班牙的微震勘测数据分析结果(ResultsofAnalysisofDataofMicroseismicSurveyatLanzaroteIsland,Canary,Spain)”(纯粹与应用地球物理(PureAppl.Geophys)，2004年，161，第1561到1578页)；或博乐格B.(BoullengerB.)的“用于CO2监测的时移地震噪声干涉测量的有限差可行性建模(FiniteDifferenceFeasibilityModellingofTime-lapseSeismicNoiseInterferometryforCO2monitoring)”(代尔夫特理工大学，硕士论文，2012年)。本地噪声是由通过的车辆的振动、装备、设备或人类的工作引起，或由采矿或其它活动产生的碰撞引起。所分析地震信号的频率通常在区域性噪声的情况下是从0.1Hz到3Hz，并且在本地噪声的情况下最多30Hz。与所分析地面地震波的频率相关联的最大深度范围可以从几十米到几千米。
低频无源地震探测(LFS)方法的发展现在可能通过使用现代技术产生宽带传感器并且通过使用平行计算提高计算能力。LFS方法使用多小时记录的数据，其迫使创建新的应用用于记录、处理以及解释这些数据。在LFS方法中人们可以区分两种基本的方法：微震探测(MS)方法和地震干涉测量(IS)方法。
当使用地震探测方法时，在多个移动测量站和固定基准站上进行记录，并且解释主要是基于地震噪声的垂直分量。在每一点处每至少1小时执行一次记录以获得其平稳性。
反过来，使用地震干涉测量方法时，连续进行记录，采样选择用于特定任务而无需移动传感器。所调查层的深度越浅，采样步骤就越少。可以使用瑞利类型的表面波和/或水平洛夫波的垂直分量来进行记录。记录长度取决于：所测试区域的波场的特性，包含主频和传播方向，以及所执行的任务的类型，例如：监视变化、定位弱化区域、确定介质结构。
国际专利申请WO2012044480(A2)揭示了地震干涉测量的数据处理方法以借助于由调查船沿海底拖曳的低频地震传感器来研究海床的地质结构，其允许海床以下地质结构的三维成像。所述方法使用格林函数来计算。
美国专利申请US2011069580(A1)描述了用于在确定地球的地下地质结构时修改地震干涉测量的方向性的装置。
另一美国专利申请US2010315902(A1)揭示了一种使用记录由在岩体中发生的地震现象产生的地震信号的地震干涉测量，用无源层析成像来成像地下结构的方法。此方法还包含适配所检测信号的过滤速度。
科扎阿尼R.(CzarnyR)在“地震干涉测量方法概述(Anoverviewofthemethodofseismicinterferometry)”(采矿综述(MiningReview)，2014年，第7期)中描述了地震干涉测量方法，所述方法在于用通过一对传感器记录的地震信号的相互相关或去卷积操作来映射这对传感器之间的介质的冲激响应(格林函数)。此方法提供了非常广泛范围的应用，从成像地球地壳的地下结构以及监视变化到地质工程应用。
马卡克H.(MarcakH.)、皮雷茨基Z.(PileckiZ.)、伊萨科夫Z.(IsakowZ.)和科扎阿尼R.(CzarnyR)在“在采矿中使用地震干涉测量方法的可能性(Possibilitiesofusingseismicinterferometrymethodinmining)”(采矿综述(MiningReview)，2014年，第7期)中描述了在采矿的地球物理学领域中使用此方法的趋势，并呈现了由采矿产生的噪声的分析结果，其可以在使用地震干涉测量方法的研究中使用。呈现了此方法的数学和物理基础。如上文已经提到的，在地震干涉测量方法中连续地执行记录，同时针对特定地质和采矿条件个别地调整采样。所调查层的地质结构越浅，采样步骤就越少。使用瑞利表面波和/或水平洛夫波的垂直分量来进行记录。记录长度取决于：所测试区域的波场的特性，即主频和传播方向，以及任务的类型，即监视变化、或定位弱化区域、或确定介质结构。
伊萨科夫Z.(IsakowZ.)、皮雷茨基Z.(PileckiZ.)和斯俄子奇P.(Sierodzki.P.)在“现代的低频无源地震LOFRES系统(ThemodernLOFRESsystemoflow-frequencypassiveseismics)”(采矿综述(MiningReview)，2014年，第7期)中示出了一种设计为使用地震噪声通过低频无源地震探测LFS方法测试地下地质层的系统LOFRES。此系统由于其测量功能而被设计成用来执行微震探测方法和地震干涉测量方法。所述系统由固定的中央处理单元构成，所述中央处理单元无线地连接到移动数据采集模块，且随后经由WI-FI自主接入点连接到测量站和基准站。
三轴低频测量站配备有经由模数转换器连接到微处理器上的三轴低频地震传感器，并且配备有大容量的内部存储器、GPS接收器、无线传输系统以及电池。测量站使用三轴低频地震传感器记录三个分量中的振动的速度。数据存储在传感器的具有几千兆字节容量的内部非易失性存储器中。通过附接到每一测量站上的GPS接收器提供时间同步。以WDS模式工作的自主接入点WI-FI以广泛区域增加WI-FI传输范围。
现有技术论述
以前使用的通过有源地震层析成像来分析位于采矿巷道上方的层的地质结构以及相对应力变化的方法及系统的主要缺点是它们的繁琐性。这与侵入性方法的使用相关联，特别是在所观察区域中矿场的地面用于住宅或工业建筑物和相关技术基础设施时的区域中。
由伊萨科夫Z.(IsakowZ.)、皮雷茨基Z.(PileckiZ.)和斯俄子奇P.(Sierodzki.P.)在上述文章“现代的低频无源地震LOFRES系统(ThemodernLOFRESsystemoflow-frequencypassiveseismics)”中所揭示，以地震探测MS的形式使用一种低频无源地震探测方法来识别位于所选择采矿巷道上方的层的地质结构以及相对应力变化。尽管是非侵入性的并且具有良好的水平分辨率(波长的大约3.5％)，但是所述方法在实际使用中低效且困难，因为它需要布置在不受采矿影响的区域中的基准站，并且对用于成像的低频噪声的非平稳性敏感。MS方法还具有垂直成像的低分辨率(波长的约30％)。虽然对噪声非平稳性不敏感且适于长期现场使用，但是地震干涉测量IS方法也以垂直成像的低分辨率为特征，并且在此已知的形式中不可用于分析位于地下矿场巷道上方的层的地质结构以及相对应力变化。
发明内容
发明目标
本发明的目标是提供用于通过低频无源地震探测非侵入性确定所选择矿场巷道上方的层的的地质结构以及相对应力变化的新的且更可靠的方法及系统。能够在周期性的或持续的基础上在本地构造冲击之前就地质异常状态和应力相对变化集中增大进行报警是必需的。
发明主旨
根据本发明的方法的特征在于，将通过在时间上严格相关地记录来自地面系统的低频地震噪声以及由采矿产生的地震碰撞而获得的来自移动测量数据记录器和来自矿场地震系统的中心站的数据传输到固定处理中心。随后，使用用于噪声记录的地震干涉测量方法以及用于冲击记录的无源速度和/或衰减层析成像，来处理优选地30秒的时间窗中的呈由采矿产生的低频地震噪声和地震突发的三分量记录形式的所记录测量数据。
在此基础上，针对岩体的所调查部分以无源速度和/或衰减层析成像方法确定横波速度的等值线以及纵波速度的等值线和/或衰减。最后，它们反映位于采矿巷道上方的层中相对应力变化的平均状态。在采矿冲击的时刻，定位其震源坐标的参数以及在岩层突裂中心中计算出的此事件的发生时间与以下相关：三轴低频纵波测量站的记录中由冲击产生的波的首次到达时刻，以及从纵波的首次到达至达到冲击记录的首次最大值所计算出的信号的相对应上升时间。上述纵波记录在矿场地面上。
在固定处理中心中，使用所采用的用于岩体所调查区域的相对应力变化的可接受阈值，对采矿巷道上方岩体的所研究区域的层中相对应力变化的状态持续进行比较分析。在当前测得的相对应力变化增大超过阈值的情况下，用信号发送已发生这种状况的位置。随后，当在可视化信令模块中观测到相对应力变化时，将地震干涉测量层析成像以及平均层析成像的空间结果从固定处理中心传输到矿场地震系统的中心站。
在根据本发明的系统中，固定处理中心优选地经由GSM调制解调器连接到移动测量数据记录器以及到矿场地震系统的中心站上。此站连接到时钟以及可视化信令模块上，并且经由本质安全的数字传输系统连接到具有至少四个地震地下测震站和/或具有至少四个地震听地器站的矿场远程传输网络上。
发明效果
允许通过无源低频地震干涉测量方法以及无源地震层析成像方法(所述方法使用由采矿产生的地震冲击)的相关使用来现行有效地分析采矿巷道上方的层的相对应力变化，这提高了地下矿场的指定采矿区域上进行的测量的功能、分辨率以及精确性。本发明使得能够先于区域性构造突裂较早地确定其中应力集中增大的位置。在矿场地面上存在住宅和/或工业建筑物以及技术基础设施的区域中，这对适合的预防性手段的应用极为重要。
分析相对应力变化使得能够对预期危险情况进行报警。本发明允许非侵入性穿透地质层，延伸到几百米的深度，而不需要应用例如使用爆炸物等的人工地震波激发方法。测量仅仅需要安装传感器并且创建用于测量与采矿地震系统相关联的低频无源地震探测的系统。重要的是通过使用地震干涉测量方法缩短连续分析之间的间隔，其中由于使噪声记录相关联因此不要求所记录噪声的平稳性。
附图说明
本发明以附图中的示例性实施例来说明，其中图1-示意性地示出测量系统连同在矿场的地面以及地下的其分量的相互位置，以及图2-示出测量系统的方块图。
具体实施方式
使用地震干涉测量方法执行根据本发明的方法，地震干涉测量方法使用优选地从测量开始处位于五个剖面线“k”中的四十个三轴低频测量站5到5i。地震干涉测量方法包括通过使用由一对三轴低频测量站5记录的地震信号的交叉相关或去卷积来映射岩体的冲激响应(所谓的格林函数)。
为了使用地震干涉测量方法连续监视，有必要确保从直流电源8连续供应具有安全电压的电源，并且使用具有本地自主WI-FI接入点4和移动测量数据记录器3的无线数据传输网络WI-FI来无线接收数据，移动测量数据记录器3记录数据并且有利地以块传输GPRS通过GSM调制解调器2将数据发送到固定处理中心1。使用地震干涉测量方法处理来处理优选地30秒的时间窗中来自三轴低频测量站5的呈三轴低频地震噪声Dn.cz.形式的所记录测量数据。
通过无源地震层析成像方法的使用来处理三轴低频测量站5中记录的由采矿产生的地震突发Dw.cz.，为此目的使用自动确定的地震波沿着地震射线从地震突发W的出现位置过渡到三轴低频测量站5中记录的它们在地面上的位置的时间。在此基础上确定岩体7的所研究区域的横波速度和纵波速度的等值线。它们体现位于岩体7所调查区域上方的矿场巷道B上的层中的相对应力变化ΔNP。
在突裂W的时刻，坐标X、Y和Z的位置参数以及在矿场突裂的震源中的计算出的时间To与以下相关：矿场A地面上三轴低频测量站5的记录中所记录的纵波的首次到达时间Tp。这样做是为了自动地确定地震波沿着地震射线路径的过渡时间。这些时间对确定波在这些方向上的传播速度是必需的。此外，为了允许另外使用空间的无源衰减层析成像方法(所述方法在三轴低频测量站5中识别突裂信号记录之后同样识别相对应力变化)，自动地确定从由突裂引起的纵波的首次到达时间Tp到所记录信号的首次最大值的信号上升时间。
在地震干涉测量方法中使用地面波，其承载有关如以色散曲线形式描述的地质介质的结构和特性的信息。在固定处理中心1中的处理包含：通过消除趋势和记录误差进行数据质量控制，地面波在0.2Hz到15Hz范围内的过滤，一位标准化，用于每对测量站的数据交叉相关性，选择具有最佳信噪比的相关性绘制，识别瑞利波的相速色散曲线，横波速度的1D模型的色散曲线的反转、横波速度场的2D模型的开发，以及使用来自三轴低频测量站5的多个测量剖面线“k”的数据开发横波速度场的空间(3D)模型。
具有专用算法的软件允许实时地进行测量，并使用计算机化的地震干涉测量方法对其进行处理。通过处理，以横波速度的等值线形式获得所控制区域的空间观测。在不太密集的地方观察到幅度的降低，横波速度应力较小，而在较密集的地方观察到幅度的较大应力增大。
在根据本发明的方法中，对于空间无源速度或衰减层析成像，使用已知坐标X、Y和Z的经定位矿场突裂和突裂震源的所计算时间To，以及在地面上测量网格P的节点处布置的测量站5的记录中所产生波首次到达的所识别时间Tp。此层析成像的结果与低频层析成像互补并且在使所获得的层析图相关联之后，它们增加分析的可靠性。
基于地面上记录的低频噪声Dn.cz.，从通过地震干涉测量方法获得的沿着特定剖面线“k”的低频层析成像的平均映射形成横向波传播速度的空间分布。随后，纵向波传播速度的空间分布从基于矿场突裂W的无源速度和/或衰减层析成像来形成，其体现矿场巷道B上方岩体7的所研究区域的层中的相对应力变化状态ΔN。
随后，在固定处理中心1中，使用所采用的用于岩体7所研究区域的相对应力变化的可接受阈值ΔNPgr，对相对应力变化ΔNP持续进行比较分析。在所测得的持续相对应力变化ΔNP增大超过阈值(ΔNP≥ΔNPgr)的情况下，启动报警并且将层析成像以及平均层析成像的空间结果从固定处理中心1发送到矿场地震系统的中心站10。在那里，在连接到所述站上的可视化信令模块11中检测到的相对应力变化ΔNP。
在图1和图2中示出的根据本发明布置中，固定处理中心1在WI-FI网络中经由GSM调制解调器2无线连接到移动测量数据记录器3，并且连续通过网状配置的两个自主WI-FI接入点4连接到四十个固定的三轴低频测量站5的数字输出端，所述三轴低频测量站在沿着剖面线“k”定位的五个测量面板6中分组。
在所调查采矿区域7上方布置具有三轴低频测量站5(SP)以及测量面板6-6k的测量网格的测量节点P。它们的数目和位置取决于此区域的大小在每种情况下个别地确定。三轴低频测量站5放置于地面中并且经由用避雷器固定的地下电信电缆线路供应恒定电压直流电源8。三轴低频测量站5中的每一者配备有三轴低频地震传感器5a，其经由模/数转换器5b连接到微处理器5c，具有大容量的内部非易失性存储器5d。微处理器5c连接到电池5e上，并且通过自动充电电路5f连接到直流电源8和GPS时钟。
固定处理中心1还在以太网网络9中连接到配备有GPS时钟和可视化信令模块11的矿场地震系统的中心站10。随后，矿场地震系统的中心站10经由本质安全的数字传输系统12，通过矿场数字电信网络13连接到地下测震站14和地震听地器站15的数字输出端。三轴低频测量站5连续地并且与GPS时钟同步记录来自三轴低频地震传感器5a的呈低频地震噪声Dn.cz.形式的测量数据的时间。在通过具有模数转换器5b的微处理器5c将这些数据处理至数字形式之后，所述数据存储在大容量的内部非易失性存储器5d中。
随后，三轴低频测量站5在数百米的距离上通过增大传输范围的自主接入点WI-FI4将上述测量数据传送到移动测量数据记录器3。所述数据也存储在那里，并且优选地通过GSM通信调制解调器2或替代地经由有线因特网网络周期性地或连续地传输到固定处理中心1。
Claims (3)
1.一种分析位于地下矿场采矿巷道上方的层的地质结构以及相对应力变化的方法，其包括使用岩体的所调查区域上方地面上安装的测量所述岩体振动的自主低频测量站，通过用于地震噪声三轴记录的地震干涉测量来测量这些层的结构特性，并且将测量数据传输到移动测量数据记录器，其特征在于，将通过在时间上密切相关地记录来自地面系统的低频地震噪声(Dn.cz.)以及由采矿产生的地震突发(Dw.cz.)而获得的来自移动测量数据记录器(3)和来自矿场地震系统的中心站(10)的数据传输到固定处理中心(1)，并且随后，使用用于所述噪声记录的地震干涉测量方法以及用于所述矿场突裂记录的无源速度和/或衰减层析成像，来处理优选地30秒的时间窗中的呈由采矿产生的低频地震噪声(Dn.cz.)和地震突发(Dw.cz.)的三轴记录形式的所述所记录测量数据，并在此基础上，针对岩体(7)的所研究区域以所述无源速度和/或衰减层析成像方法确定横波速度的等值线以及纵波速度的等值线和/或衰减，它们最后将体现采矿巷道(B)上方的所述层中相对应力变化(ΔNP)的平均状态，其中在矿场突裂(W)的时刻，位置坐标(X、Y和Z)的参数以及计算出的突裂发生时间(To)与以下相关：矿场地面上记录的波的三轴低频测量站(5)的记录中由此岩层突裂产生的纵波首次到达的时间(Tp)，以及从所述纵波的所述首次到达直到所述岩层突裂的所记录信号(W)达到每一测量站(5)中的首次最大值为止所述信号的相对应上升时间。
2.根据权利要求1所述的方法，其特征在于，在所述固定处理中心(1)中，使用用于所述岩体(7)所研究区域的相对应力变化的可接受阈值(ΔNPgr)，对所述矿场巷道(B)上方所述岩体(7)的所述所研究区域的所述层中相对应力变化(ΔNP)的状态持续进行比较分析，并且在当前测得的相对应力变化(ΔNP)增大超过所述阈值(ΔNP≥ΔNPgr)的情况下，用信号发送已发生这种状况的位置，并且随后，当在所述可视化信令模块(11)中观测到所测得的相对应力变化(ΔNP)时，将通过地震干涉测量执行的层析成像以及平均层析成像的空间结果从所述固定处理中心(1)传输到所述矿场地震系统的中心站(10)。
3.一种用于分析位于地下矿场采矿巷道上方的层的地质结构以及相对应力变化的系统，其包括无线地连接到移动测量数据记录单元上且随后经由WI-FI自主接入点连接到三轴低频测量站上的固定处理中心，所述三轴低频测量站配备有经由模数转换器连接到微处理器上的三轴低频地震传感器，具有大容量的内部非易失性存储器、GPS接收器、无线传输系统以及电池，对于如在权利要求1到2中界定的所述方法的使用，其特征在于，所述固定处理中心(1)优选地经由GSM调制解调器(2)连接到所述移动测量数据记录器(3)以及到所述矿场地震系统的中心站(10)上，所述中心站连接到时钟(GPS)以及所述可视化信令模块(11)上，并且经由本质安全的数字传输系统(12)以及矿场远程传输网络(13)连接到至少四个地下测震站(14)和/或至少四个地下地震听地器站(15)上。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PL409989A PL230219B1 (pl) | 2014-10-30 | 2014-10-30 | Sposób i układ do analizy struktury geologicznej i względnych zmian naprężeń w warstwach usytuowanych nad wyrobiskami górniczymi kopalni głębinowej |
| PLP.409989 | 2014-10-30 | ||
| PCT/PL2014/000128 WO2015002558A2 (en) | 2014-10-30 | 2014-11-07 | Method and system for analysis of geological structure and relative changes in stress in the layers located above the workings of an underground mine |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105765408A true CN105765408A (zh) | 2016-07-13 |
| CN105765408B CN105765408B (zh) | 2019-05-07 |
Family
ID=52004033
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201480002443.2A Active CN105765408B (zh) | 2014-10-30 | 2014-11-07 | 用于分析位于地下矿场巷道上方的层的方法及系统 |
Country Status (5)
| Country | Link |
|---|---|
| CN (1) | CN105765408B (zh) |
| PL (1) | PL230219B1 (zh) |
| RU (1) | RU2587521C1 (zh) |
| UA (1) | UA118543C2 (zh) |
| WO (1) | WO2015002558A2 (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110794460A (zh) * | 2019-11-15 | 2020-02-14 | 中国矿业大学 | 应力值变化方向约束下的二维矿震全波形反演方法 |
| CN113015926A (zh) * | 2018-11-14 | 2021-06-22 | Bp探索操作有限公司 | 无源地震成像 |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20180203144A1 (en) * | 2015-05-20 | 2018-07-19 | Optasense, Inc. | Interferometric Microseismic Imaging Methods and Apparatus |
| PL422137A1 (pl) * | 2017-07-10 | 2019-01-14 | Pytel Witold | Sposób prognozowania samoistnych zjawisk sejsmicznych indukowanych eksploatacją górniczą |
| CN112051548B (zh) * | 2020-08-11 | 2024-03-22 | 武汉工程大学 | 一种岩爆监测和定位方法、装置和系统 |
| US11323285B1 (en) | 2020-08-28 | 2022-05-03 | Earthsystems Technologies, Inc. | Architecture for a multichannel geophysical data acquisition system and method of use |
| CN112346123B (zh) * | 2020-11-06 | 2023-02-10 | 中国地震灾害防御中心 | 一种地震数据处理via双参数分析方法 |
| US11808797B1 (en) | 2021-03-19 | 2023-11-07 | Earthsystems Technologies, Inc. | Hemispherical dome electrode configuration and method of use |
| CN113404523B (zh) * | 2021-07-05 | 2023-11-10 | 淮北市平远软岩支护工程技术有限公司 | 一种基于卸压爆破的冲击地压监测系统 |
| CN113703046A (zh) * | 2021-08-31 | 2021-11-26 | 中煤科工集团重庆研究院有限公司 | 井下全空间地震波隐伏构造识别方法及系统 |
| CN113985482B (zh) * | 2021-10-28 | 2023-11-03 | 西安科技大学 | 基于煤矿井下通讯光缆的矿震震源定位方法 |
| CN114895353B (zh) * | 2022-05-27 | 2023-03-10 | 中国矿业大学 | 井地一体微震监测系统监测单元采集数据授时对齐方法 |
| CN115542381B (zh) * | 2022-09-26 | 2024-02-02 | 徐州弘毅科技发展有限公司 | 基于三分向监测仪的矿震井地一体化融合监测系统及方法 |
| CN115933803B (zh) * | 2023-01-09 | 2023-06-13 | 江苏东成工具科技有限公司 | 一种设备控制方法、设备及计算机可读介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101581789A (zh) * | 2009-06-23 | 2009-11-18 | 刘盛东 | 矿井工作面巷间震波ct探测方法 |
| CN102597808A (zh) * | 2009-11-02 | 2012-07-18 | 界标制图有限公司 | 采用倾斜的横向各向同性的3d逆向时间偏移的地震成像系统与方法 |
| CN102788995A (zh) * | 2012-08-02 | 2012-11-21 | 中煤科工集团西安研究院 | 以切割震动为地震信号的煤矿工作面探测方法 |
| US20130191044A1 (en) * | 2011-03-21 | 2013-07-25 | Schlumberger Technology Corporation | Method and system for presenting stress domain |
| CN103728655A (zh) * | 2013-12-24 | 2014-04-16 | 天地科技股份有限公司 | 一种回采工作面冲击危险性采前预评价方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7859942B2 (en) * | 2007-03-01 | 2010-12-28 | Christof Stork | Measuring and modifying directionality of seismic interferometry data |
| CA2747426C (en) * | 2009-01-09 | 2017-05-23 | Exxonmobil Upstream Research Company | Hydrocarbon detection with passive seismic data |
| US8218394B2 (en) * | 2009-06-16 | 2012-07-10 | Microseismic, Inc. | Method for imaging the earths subsurface using passive seismic interferometry and adaptive velocity filtering |
| US8838392B2 (en) * | 2009-10-05 | 2014-09-16 | Westerngeco L.L.C. | Noise attenuation in passive seismic data |
| US8737165B2 (en) * | 2010-10-01 | 2014-05-27 | Westerngeco L.L.C. | Interferometric seismic data processing for a towed marine survey |
| US20130003499A1 (en) * | 2011-06-28 | 2013-01-03 | King Abdulaziz City For Science And Technology | Interferometric method of enhancing passive seismic events |
-
2014
- 2014-10-30 PL PL409989A patent/PL230219B1/pl unknown
- 2014-11-07 UA UAA201501087A patent/UA118543C2/uk unknown
- 2014-11-07 CN CN201480002443.2A patent/CN105765408B/zh active Active
- 2014-11-07 WO PCT/PL2014/000128 patent/WO2015002558A2/en active Application Filing
- 2014-11-07 RU RU2015101763/28A patent/RU2587521C1/ru active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101581789A (zh) * | 2009-06-23 | 2009-11-18 | 刘盛东 | 矿井工作面巷间震波ct探测方法 |
| CN102597808A (zh) * | 2009-11-02 | 2012-07-18 | 界标制图有限公司 | 采用倾斜的横向各向同性的3d逆向时间偏移的地震成像系统与方法 |
| US20130191044A1 (en) * | 2011-03-21 | 2013-07-25 | Schlumberger Technology Corporation | Method and system for presenting stress domain |
| CN102788995A (zh) * | 2012-08-02 | 2012-11-21 | 中煤科工集团西安研究院 | 以切割震动为地震信号的煤矿工作面探测方法 |
| CN103728655A (zh) * | 2013-12-24 | 2014-04-16 | 天地科技股份有限公司 | 一种回采工作面冲击危险性采前预评价方法 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113015926A (zh) * | 2018-11-14 | 2021-06-22 | Bp探索操作有限公司 | 无源地震成像 |
| CN110794460A (zh) * | 2019-11-15 | 2020-02-14 | 中国矿业大学 | 应力值变化方向约束下的二维矿震全波形反演方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2015002558A3 (en) | 2015-09-11 |
| WO2015002558A2 (en) | 2015-01-08 |
| CN105765408B (zh) | 2019-05-07 |
| UA118543C2 (uk) | 2019-02-11 |
| PL409989A1 (pl) | 2016-05-09 |
| RU2587521C1 (ru) | 2016-06-20 |
| PL230219B1 (pl) | 2018-10-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105765408A (zh) | 用于分析位于地下矿场巷道上方的层的地质结构以及相对应力变化的方法及系统 | |
| US9523790B1 (en) | Hybrid sensing apparatus and method | |
| EP1946129B1 (en) | Continuous reservoir monitoring for fluid pathways using 3d microseismic data | |
| EP2255222B1 (en) | Method for imaging the earth's subsurface using passive seismic sensing | |
| CA2620819C (en) | Method and apparatus for imaging permeability pathways of geologic fluid reservoirs using seismic emission tomography | |
| EP2326795B1 (en) | Determining characteristics of a subterranean body using pressure data and seismic data | |
| CN106501848B (zh) | 一种隧道掘进过程中隐性断层超前物探方法 | |
| CN103336297B (zh) | 微破裂向量扫描方法 | |
| CN102879805A (zh) | 一种基于钻孔与地面相结合的地震波空间探测方法 | |
| NO20170663A1 (en) | Three-dimensional rock properties using cross well seismic | |
| CN113153430B (zh) | 巷道围岩破坏声发射定位与波速成像监测及灾变预警方法 | |
| Roux et al. | Microseismic activity within a serac zone in an alpine glacier (Glacier d’Argentiere, Mont Blanc, France) | |
| CA2961168A1 (en) | Integrating vertical seismic profile data for microseismic anisotropy velocity analysis | |
| Luo et al. | Tomographic imaging of rock conditions ahead of mining using the shearer as a seismic source—a feasibility study | |
| Xu et al. | Optimal design of microseismic monitoring networking and error analysis of seismic source location for rock slope | |
| Feng et al. | Rockfall localization from seismic polarization considering multiple triaxial geophones and frequency bands | |
| du Toit et al. | Can DAS be used to monitor mining induced seismicity? | |
| JP2008014830A (ja) | ハイドレートの存在領域探査方法及び探査システム | |
| Fang et al. | A single station approach for subsurface anomaly detection | |
| Wang et al. | Retrieving drill bit seismic signals using surface seismometers | |
| RU2645790C1 (ru) | Способ определения границ субвертикальных протяженных объектов в геологической среде | |
| Luo et al. | Sensing roof conditions ahead of a longwall mining using the shearer as a seismic source | |
| Zhou et al. | Feasibility of local anisotropy estimation using DAS VSP | |
| Henninges et al. | Exploration and monitoring with distributed acoustic sensing at the EGS Site Groß Schönebeck | |
| CN105929445A (zh) | 微破裂向量扫描方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 |