CN113985482B - 基于煤矿井下通讯光缆的矿震震源定位方法 - Google Patents
基于煤矿井下通讯光缆的矿震震源定位方法 Download PDFInfo
- Publication number
- CN113985482B CN113985482B CN202111262256.6A CN202111262256A CN113985482B CN 113985482 B CN113985482 B CN 113985482B CN 202111262256 A CN202111262256 A CN 202111262256A CN 113985482 B CN113985482 B CN 113985482B
- Authority
- CN
- China
- Prior art keywords
- optical cable
- communication optical
- communication
- deformation
- mine
- 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
- 230000003287 optical effect Effects 0.000 title claims abstract description 180
- 238000004891 communication Methods 0.000 title claims abstract description 100
- 239000003245 coal Substances 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 28
- 230000010287 polarization Effects 0.000 claims abstract description 37
- 238000005065 mining Methods 0.000 claims abstract description 25
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 16
- 239000011707 mineral Substances 0.000 claims abstract description 16
- 238000005259 measurement Methods 0.000 claims description 18
- 238000004458 analytical method Methods 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 8
- 239000011435 rock Substances 0.000 claims description 8
- 230000035945 sensitivity Effects 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 238000009826 distribution Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 241000084490 Esenbeckia delta Species 0.000 claims description 3
- 238000005452 bending Methods 0.000 claims description 3
- 230000035939 shock Effects 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 41
- 239000013307 optical fiber Substances 0.000 abstract description 4
- 238000001514 detection method Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011946 reduction process 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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
- G01B11/168—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge by means of polarisation
-
- 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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/22—Transmitting seismic signals to recording or processing apparatus
- G01V1/226—Optoseismic systems
-
- 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
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及矿震震源定位技术领域,公开了基于煤矿井下通讯光缆的矿震震源定位方法,包括以下步骤:采用矿井下通讯光缆为载体,在地面上将通讯光缆的一端与激光器连接,将通讯光缆的另一端与波分复用器连接,波分复用器的一端连接有偏振测量仪;通过控制激光器对通讯光缆发出具有偏振态的光信号,通讯光缆另一端的偏振测量仪对接收到的光信号进行解调,以获得通讯光缆上所有测点的第一斯托克斯参量s1的峰值Δs1max。本发明利用光纤偏振信息进行传感,创新性的利用煤矿井下既有的矿用通讯光缆替代价格昂贵、铺设困难的专用的应变传感光缆,可以有效解决煤矿现有微震监测系统布置位置单一、传感器覆盖范围小、监测范围有限等问题。
Description
技术领域
本发明涉及矿震震源定位技术领域,特别涉及基于煤矿井下通讯光缆的矿震震源定位方法。
背景技术
随着国家大力推动“互联网+”的行业解决方案,以及智慧化矿山的建设热潮,使得越来越多的先进智能设备及装置应用到煤矿井下,来满足不断增长的矿山安全监测需求,解决煤矿现有监测系统抗灾能力差、安装不便、监测范围小等问题,提高煤矿的安全生产和安全监测水平。特别是对于传统的矿震监测系统来说,微震传感器探头和数据采集主机采用有线连接,一般只能采集8到12通道的同一类传感数据,这种的方式不仅仅限制了采集传感器的种类,也限制了采集的数量,导致监测范围有限,监测精度难以提升的困境。
矿震是由采矿诱发的矿井地震,矿震发生时井下围岩迅速释放能量,往往造成井下巷道或采掘工作面突然破坏、地面震动、房屋损坏等,严重时造成人员伤亡,随着我国煤矿开采深度增加,矿震事故频发,严重威胁着我国煤矿的安全生产。目前矿震监测主要是由井下专业的微震监测系统来完成,矿震发生时,震动波检测台站将监测到的震动波数据记录发送至地面监测系统,通过提取矿震波信息计算得到发送的时间、地点和能量大小。现有矿震监测设备,设备价格昂贵,机构复杂,每个矿井布点很少,造成网络覆盖面积小,存在大量的监测盲区,难以及时、准确的进行监测定位并达到有效预警目的,为此,我们提出基于煤矿井下通讯光缆的矿震震源定位方法。
发明内容
本发明的主要目的在于提供基于煤矿井下通讯光缆的矿震震源定位方法,可以有效解决背景技术中传统的矿震监测系统微震传感器探头和数据采集主机采用有线连接,一般只能采集8到12通道的同一类传感数据,限制了采集传感器的种类,也限制了采集的数量,导致监测范围有限,监测精度难以提升的困境,同时现有矿震监测设备,设备价格昂贵,机构复杂,每个矿井布点很少,造成网络覆盖面积小,存在大量的监测盲区,难以及时、准确的进行监测定位并达到有效预警目的的问题。
为实现上述目的,本发明采取的技术方案为:
基于煤矿井下通讯光缆的矿震震源定位方法,包括以下步骤:
步骤一:采用矿井下通讯光缆为载体,在地面上将通讯光缆的一端与激光器连接,将通讯光缆的另一端与波分复用器连接,波分复用器的一端连接有偏振测量仪;
步骤二:通过控制激光器对通讯光缆发出具有偏振态的光信号,通讯光缆另一端的偏振测量仪对接收到的光信号进行解调,以获得通讯光缆上所有测点的第一斯托克斯参量s1的峰值Δs1max;
步骤三:对第一斯托克斯参量s1的峰值进行反演分析,利用通讯光缆的第一斯托克斯参量s1获得通讯光缆变形量,感知矿震所引起的光缆形变量Δx,其中,光缆形变量公式为:
式中,Δx为通讯光缆某一测点位置处光缆形变量;D光缆直径;a为光缆受力形变导致的折射率变化系数,矿用通讯光缆一般取0.11×1012;Δs1max1为第一斯托克斯参量s1的峰值变化量,由偏振仪测得;K为Δs1max与矿用光缆形变量的灵敏度系数,可以由室内试验测得;
步骤四:为了区分井下矿震引起的光缆变形和其他因素引起的光缆变形,将Δx≥1mm设置为判别条件,当Δx≥1mm时,认为是矿震引起光缆变形;
步骤五:按照上述方法对通讯光缆进行实时监测,当发生矿震时,测量获得通讯光缆上所有测量点的第一斯托克斯参量s1的峰值,对于测量点峰值频率进行分析计算,获得所有测量点的位移,当Δx≥1mm时,将这些测量点的位置及测量点监测到时直接上传至云端服务器进行分析处理,通过空间定位算法实现矿震震源定位;
步骤六:首先选取4个符合条件的测点为一个参考平面,在参考平面上将4个测点命名为:A、B、C、D,其中三角形的AB、AC边长分别为LAB=a,LAC=b,假设AB与AC垂直,D为AB中点,并选取A为空间直角坐标原点,震源位置为点P(x,y,z),P点到A、B、C、D四点的距离分别为Ra、Rb、Rc和Rd;
步骤七:当矿震在P点产生应力波引发的岩层变形进而作用到光缆上压力F致使光缆发生形变,通过偏振检测仪和矿震震源定位算法判别出符合条件的测量点,选取其中4个测点发生变形的时刻分别为ta,tb,tc和td,则4个点到震源P的距离有如下关系,
式中:t1=ta-tb;t2=ta-tc;t3=ta-td;
结合式(2)、式(3)可得式(4):
令Ra-Rb=R1,Ra-Rd=R2,得到Ra的一元一次方程:
可以求出Ra的值,代入到式(1)得到Rb和Rc的值,进而根据P点坐标得方程组:
由此可以得到震源P点的坐标;
步骤八:重复步骤四、步骤五、步骤六与步骤七,通过对多组符合要求的测量点进行反演,获得到更精确的矿震震源位置。
进一步的,步骤三中所述Δs1max与通讯光缆受力大小呈线性关系,一般为9×103,而通讯光缆受力大小与煤矿井下悬挂的矿用光缆存在一定的关系,其中矿用光缆受到一定的预紧力为σ,其次当矿震引发的岩层变形进而作用到光缆上压力F致使光缆发生形变,引起光缆的弯矩M发生变化。
进一步的,步骤三中所述光缆形变量Δx与光缆作用力之间关系可以由室内试验测试得出,对光缆施加恒定的力,然后记录其形变量变化,光缆受力与光缆型变量的分布曲线的斜率的倒数再乘以9×103即为Δs1max与矿用光缆形变量的灵敏度系数K。
进一步的,步骤四中由于考虑测量时间和测量精度的平衡,一般偏振测量仪的空间分辨率为1m,采样间隔为0.05m,也就是说对于50km长的通讯光缆,它具有1000000个测量点。
与现有技术相比,本发明具有如下有益效果:
利用光纤偏振信息进行传感,创新性的利用煤矿井下既有的矿用通讯光缆替代价格昂贵、铺设困难的专用的应变传感光缆,可以有效解决煤矿现有微震监测系统布置位置单一、传感器覆盖范围小、监测范围有限等问题,同时研发煤矿环境下高效精确的定位算法,对通讯光缆偏振态信号进行反演处理,以实现煤矿井下实时监测并精确定位震源,并且系统中使用到的都是高性价比、通用性较强的技术平台和设备,与现有的百万级矿震监测系统相比,价格上有显著优势。同时利用井下既有的矿用通讯光缆替代价格昂贵、铺设困难的专用的应变传感光缆,使得监测覆盖范围更大,更适用于全矿井范围的震源监测及定位,具有实时监测、全范围立体监测、空间定位、全数字化数据采集、存储和处理、远程监测和信息远传输送以及多用户计算机可视化监控与分析的特点。在煤矿井下中,通讯光缆不像地面上那样易受影响,使得基于偏振态的检测方式成为高效低成本的震源定位及监测预警方法。及时准确的获取,传递并传播预警信息,在减灾过程中,具有决定性作用。成功在不额外增加通信负担的前提下,通过测量光信号的偏振态,首次将传统通信光纤变为“震源监测定位网络”,实现了震源信号和定位信息的光速传递,为自救避难争取了宝贵时间。可以在不影响井下光缆通信功能的情况下,实现对煤矿井下矿震的探测定位,且不会对通信系统增加额外负载。
附图说明
图1为本发明基于煤矿井下通讯光缆的矿震震源定位方法的激光器、通讯光缆、波分复用器与偏振测量仪的连接结构示意图。
图2为本发明基于煤矿井下通讯光缆的矿震震源定位方法的矿井内部通讯光缆分布布置图。
图3为本发明基于煤矿井下通讯光缆的矿震震源定位方法中选取的4个符合条件的测点的参考平面坐标轴。
图4为本发明基于煤矿井下通讯光缆的矿震震源定位方法的通讯光缆受力大小线性图。
图中:1、激光器;2、通讯光缆;3、波分复用器;4、偏振测量仪。
具体实施方式
为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体实施方式,进一步阐述本发明。
如图1-4所示,基于煤矿井下通讯光缆的矿震震源定位方法,包括以下步骤:
步骤一:采用矿井下通讯光缆2为载体,在地面上将通讯光缆2的一端与激光器1连接,将通讯光缆2的另一端与波分复用器3连接,波分复用器3的一端连接有偏振测量仪4;
步骤二:通过控制激光器1对通讯光缆2发出具有偏振态的光信号,通讯光缆2另一端的偏振测量仪4对接收到的光信号进行解调,以获得通讯光缆2上所有测点的第一斯托克斯参量s1的峰值Δs1max;
步骤三:对第一斯托克斯参量s1的峰值进行反演分析,利用通讯光缆2的第一斯托克斯参量s1获得通讯光缆2变形量,感知矿震所引起的光缆形变量Δx,其中,光缆形变量公式为:
式中,Δx为通讯光缆2某一测点位置处光缆形变量;D光缆直径;a为光缆受力形变导致的折射率变化系数,矿用通讯光缆2一般取0.11×1012;Δs1max为第一斯托克斯参量s1的峰值变化量,由偏振仪测得;K为Δs1max与矿用光缆形变量的灵敏度系数,可以由室内试验测得;
步骤四:为了区分井下矿震引起的光缆变形和其他因素引起的光缆变形,将Δx≥1mm设置为判别条件,当Δx≥1mm时,认为是矿震引起光缆变形;由于考虑测量时间和测量精度的平衡,一般偏振测量仪4的空间分辨率为1m,采样间隔为0.05m,也就是说对于50km长的通讯光缆2,它具有1000000个测量点;
步骤五:按照上述方法对通讯光缆2进行实时监测,当发生矿震时,测量获得通讯光缆2上所有测量点的第一斯托克斯参量s1的峰值,对于测量点峰值频率进行分析计算,获得所有测量点的位移,当Δx≥1mm时,将这些测量点的位置及测量点监测到时直接上传至云端服务器进行分析处理,通过空间定位算法实现矿震震源定位;
步骤六:首先选取4个符合条件的测点为一个参考平面,在参考平面上将4个测点命名为:A、B、C、D,其中三角形的AB、AC边长分别为LAB=a,LAC=b,假设AB与AC垂直,D为AB中点,并选取A为空间直角坐标原点,震源位置为点P(x,y,z),P点到A、B、C、D四点的距离分别为Ra、Rb、Rc和Rd;
步骤七:当矿震在P点产生应力波引发的岩层变形进而作用到光缆上压力F致使光缆发生形变,通过偏振检测仪和矿震震源定位算法判别出符合条件的测量点,选取其中4个测点发生变形的时刻分别为ta,tb,tc和td,则4个点到震源P的距离有如下关系,
式中:t1=ta-tb;t2=ta-tc;t3=ta-td;
结合式(2)、式(3)可得式(4):
令Ra-Rb=R1,Ra-Rd=R2,得到Ra的一元一次方程:
可以求出Ra的值,代入到式(1)得到Rb和Rc的值,进而根据P点坐标得方程组:
由此可以得到震源P点的坐标;
步骤八:重复步骤四、步骤五、步骤六与步骤七,通过对多组符合要求的测量点进行反演,获得到更精确的矿震震源位置。
如图4所示,Δs1max与通讯光缆2受力大小呈线性关系,一般为9×103,而通讯光缆2受力大小与煤矿井下悬挂的矿用光缆存在一定的关系,其中矿用光缆受到一定的预紧力为σ,其次当矿震引发的岩层变形进而作用到光缆上压力F致使光缆发生形变,引起光缆的弯矩M发生变化;
光缆形变量Δx与光缆作用力之间关系可以由室内试验测试得出,对光缆施加恒定的力,然后记录其形变量变化,光缆受力与光缆型变量的分布曲线的斜率的倒数再乘以9×103即为Δs1max与矿用光缆形变量的灵敏度系数K。
需要说明的是,本发明为基于煤矿井下通讯光缆2的矿震震源定位方法,在使用时,现在矿井运输巷内构建如图1所示的检测系统,包括激光器1和偏振测量仪4,激光器1安装在通讯光缆2一端,通讯光缆2的另一端安装有波分复用器3,偏振测量仪4安装在波分复用器3的一端,其中激光器1实时发出具有偏振态的光信号到通讯光缆2内,通讯光缆2另一端输出光信号到偏振测量仪4;本发明利用光纤偏振信息进行传感,创新性的利用煤矿井下既有的矿用通讯光缆2替代价格昂贵、铺设困难的专用的应变传感光缆,可以有效解决煤矿现有微震监测系统布置位置单一、传感器覆盖范围小、监测范围有限等问题。
本发明采用矿井下通讯光缆2为载体,在地面上将通讯光缆2的一端与激光器1连接,将通讯光缆2的另一端与波分复用器3连接,波分复用器3的一端连接有偏振测量仪4;通过控制激光器1对通讯光缆2发出具有偏振态的光信号,通讯光缆2另一端的偏振测量仪4对接收到的光信号进行解调,以获得通讯光缆2上所有测点的第一斯托克斯参量s1的峰值Δs1max;然后对第一斯托克斯参量s1的峰值进行反演分析,利用通讯光缆2的第一斯托克斯参量s1获得通讯光缆2变形量,感知矿震所引起的光缆形变量Δx,其中,光缆形变量公式为:式中,Δx为通讯光缆2某一测点位置处光缆形变量;D光缆直径;a为光缆受力形变导致的折射率变化系数,矿用通讯光缆2一般取0.11×1012;Δs1max为第一斯托克斯参量s1的峰值变化量,由偏振仪测得;K为Δs1max与矿用光缆形变量的灵敏度系数,可以由室内试验测得;为了区分井下矿震引起的光缆变形和其他因素引起的光缆变形,将Δx≥1mm设置为判别条件,当Δx≥1mm时,认为是矿震引起光缆变形;由于考虑测量时间和测量精度的平衡,一般偏振测量仪4的空间分辨率为1m,采样间隔为0.05m,也就是说对于50km长的通讯光缆2,它具有1000000个测量点;按照上述方法对通讯光缆2进行实时监测,当发生矿震时,测量获得通讯光缆2上所有测量点的第一斯托克斯参量s1的峰值,对于测量点峰值频率进行分析计算,获得所有测量点的位移,当Δx≥1mm时,将这些测量点的位置及测量点监测到时直接上传至云端服务器进行分析处理,通过空间定位算法实现矿震震源定位;如图3所示,首先选取4个符合条件的测点为一个参考平面,在参考平面上将4个测点命名为:A、B、C、D,其中三角形的AB、AC边长分别为LAB=a,LAC=b,假设AB与AC垂直,D为AB中点,并选取A为空间直角坐标原点,震源位置为点P(x,y,z),P点到A、B、C、D四点的距离分别为Ra、Rb、Rc和Rd;当矿震在P点产生应力波引发的岩层变形进而作用到光缆上压力F致使光缆发生形变,通过偏振检测仪和矿震震源定位算法判别出符合条件的测量点,选取其中4个测点发生变形的时刻分别为ta,tb,tc和td,则4个点到震源P的距离有如下关系,
式中:t1=ta-tb;t2=ta-tc;t3=ta-td;
结合式(2)、式(3)可得式(4):
令Ra-Rb=R1,Ra-Rd=R2,得到Ra的一元一次方程:
[(2Ra-R1)2-a2]×[a2-R1 2]=4[(2Ra-R2)2-a2/4]×[a2/4-R2 2] (5)
可以求出Ra的值,代入到式(1)得到Rb和Rc的值,进而根据P点坐标得方程组:
由此可以得到震源P点的坐标;重复以上步骤,一次性对多个矿震震源位置进行定位,通过采用煤矿环境下高效精确的定位算法,对通讯光缆2偏振态信号进行反演处理,以实现煤矿井下实时监测并精确定位震源,并且系统中使用到的都是高性价比、通用性较强的技术平台和设备,与现有的百万级矿震监测系统相比,价格上有显著优势。
本发明利用井下既有的矿用通讯光缆2替代价格昂贵、铺设困难的专用的应变传感光缆,使得监测覆盖范围更大,更适用于全矿井范围的震源监测及定位,具有实时监测、全范围立体监测、空间定位、全数字化数据采集、存储和处理、远程监测和信息远传输送以及多用户计算机可视化监控与分析的特点。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (4)
1.基于煤矿井下通讯光缆的矿震震源定位方法,其特征在于,包括以下步骤:
步骤一:采用矿井下通讯光缆为载体,在地面上将通讯光缆的一端与激光器连接,将通讯光缆的另一端与波分复用器连接,波分复用器的一端连接有偏振测量仪;
步骤二:通过控制激光器对通讯光缆发出具有偏振态的光信号,通讯光缆另一端的偏振测量仪对接收到的光信号进行解调,以获得通讯光缆上所有测点的第一斯托克斯参量s1的峰值Δs1max;
步骤三:对第一斯托克斯参量s1的峰值进行反演分析,利用通讯光缆的第一斯托克斯参量s1获得通讯光缆变形量,感知矿震所引起的光缆形变量Δx,其中,光缆形变量公式为:
式中,Δx为通讯光缆某一测点位置处光缆形变量;D光缆直径;a为光缆受力形变导致的折射率变化系数,矿用通讯光缆取0.11×1012;Δs1max为第一斯托克斯参量s1的峰值变化量,由偏振仪测得;K为Δs1max与矿用光缆形变量的灵敏度系数,可以由室内试验测得;
步骤四:为了区分井下矿震引起的光缆变形和其他因素引起的光缆变形,将Δx≥1mm设置为判别条件,当Δx≥1mm时,认为是矿震引起光缆变形;
步骤五:当发生矿震时,测量获得通讯光缆上所有测量点的第一斯托克斯参量s1的峰值,对于测量点峰值频率进行分析计算,获得所有测量点的位移,当Δx≥1mm时,将这些测量点的位置及测量点监测到时直接上传至云端服务器进行分析处理,通过空间定位算法实现矿震震源定位;
步骤六:首先选取4个符合条件的测点为一个参考平面,在参考平面上将4个测点命名为:A、B、C、D,其中三角形的AB、AC边长分别为LAB=a,LAC=b,假设AB与AC垂直,D为AB中点,并选取A为空间直角坐标原点,震源位置为点P(x,y,z),P点到A、B、C、D四点的距离分别为Ra、Rb、Rc和Rd;
步骤七:当矿震在P点产生应力波引发的岩层变形进而作用到光缆上压力F致使光缆发生形变,通过偏振检测仪和矿震震源定位算法判别出符合条件的测量点,选取其中4个测点发生变形的时刻分别为ta,tb,tc和td,则4个点到震源P的距离有如下关系,
式中:t1=ta-tb;t2=ta-tc;t3=ta-td;
结合式(2)、式(3)可得式(4):
令Ra-Rb=R1,Ra-Rd=R2,得到Ra的一元一次方程:
可以求出Ra的值,代入到式(1)得到Rb和Rc的值,进而根据P点坐标得方程组:
由此可以得到震源P点的坐标;
步骤八:重复步骤四、步骤五、步骤六与步骤七,通过对多组符合要求的测量点进行反演,获得到更精确的矿震震源位置。
2.根据权利要求1所述的基于煤矿井下通讯光缆的矿震震源定位方法,其特征在于:步骤三中所述Δs1max与通讯光缆受力大小呈线性关系,为9×103,而通讯光缆受力大小与煤矿井下悬挂的矿用光缆存在一定的关系,其中矿用光缆受到一定的预紧力为σ,其次当矿震引发的岩层变形进而作用到光缆上压力F致使光缆发生形变,引起光缆的弯矩M发生变化。
3.根据权利要求1所述的基于煤矿井下通讯光缆的矿震震源定位方法,其特征在于:步骤三中所述光缆形变量Δx与光缆作用力之间关系可以由室内试验测试得出,对光缆施加恒定的力,然后记录其形变量变化,光缆受力与光缆型变量的分布曲线的斜率的倒数再乘以9×103即为Δs1max与矿用光缆形变量的灵敏度系数K。
4.根据权利要求1所述的基于煤矿井下通讯光缆的矿震震源定位方法,其特征在于:步骤四中由于考虑测量时间和测量精度的平衡,偏振测量仪的空间分辨率为1m,采样间隔为0.05m,也就是说对于50km长的通讯光缆,它具有1000000个测量点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111262256.6A CN113985482B (zh) | 2021-10-28 | 2021-10-28 | 基于煤矿井下通讯光缆的矿震震源定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111262256.6A CN113985482B (zh) | 2021-10-28 | 2021-10-28 | 基于煤矿井下通讯光缆的矿震震源定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113985482A CN113985482A (zh) | 2022-01-28 |
CN113985482B true CN113985482B (zh) | 2023-11-03 |
Family
ID=79743305
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111262256.6A Active CN113985482B (zh) | 2021-10-28 | 2021-10-28 | 基于煤矿井下通讯光缆的矿震震源定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113985482B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1096273A2 (en) * | 1999-10-29 | 2001-05-02 | Litton Systems, Inc. | Accoustic sensing systems for downhole seismic applications utilizing an array of fiber optic sensors |
CN101393269A (zh) * | 2008-11-06 | 2009-03-25 | 复旦大学 | 一种利用通信光缆监测地质状况的方法 |
CN201730613U (zh) * | 2010-05-18 | 2011-02-02 | 西部矿业股份有限公司 | 一种非煤矿山微震监测系统 |
CN103292889A (zh) * | 2013-05-23 | 2013-09-11 | 无锡聚为传感科技有限公司 | 一种分布式光纤振动传感器振源定位方法 |
WO2014209141A2 (en) * | 2014-10-30 | 2014-12-31 | Instytut Technik Innowacyjnych Emag | Method and system for assessing a risk of high-energy earth bursts generated by underground mining |
WO2015002558A2 (en) * | 2014-10-30 | 2015-01-08 | Instytut Technik Innowacyjnych Emag | Method and system for analysis of geological structure and relative changes in stress in the layers located above the workings of an underground mine |
CN111665568A (zh) * | 2020-07-09 | 2020-09-15 | 中油奥博(成都)科技有限公司 | 基于分布式光纤声波传感技术的微测井装置及测量方法 |
CN113417699A (zh) * | 2021-08-23 | 2021-09-21 | 辽宁工程技术大学 | 一种煤矿复合灾害分布式光纤动态监测系统及使用方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2411001B (en) * | 2004-02-10 | 2007-03-28 | Statoil Asa | Seismic exploration |
US8031555B2 (en) * | 2008-01-07 | 2011-10-04 | University Of Utah Research Foundation | Seismic location and communication system |
US9377551B2 (en) * | 2013-05-22 | 2016-06-28 | Schlumberger Technology Corporation | Method of borehole seismic surveying using an optical fiber |
-
2021
- 2021-10-28 CN CN202111262256.6A patent/CN113985482B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1096273A2 (en) * | 1999-10-29 | 2001-05-02 | Litton Systems, Inc. | Accoustic sensing systems for downhole seismic applications utilizing an array of fiber optic sensors |
CN101393269A (zh) * | 2008-11-06 | 2009-03-25 | 复旦大学 | 一种利用通信光缆监测地质状况的方法 |
CN201730613U (zh) * | 2010-05-18 | 2011-02-02 | 西部矿业股份有限公司 | 一种非煤矿山微震监测系统 |
CN103292889A (zh) * | 2013-05-23 | 2013-09-11 | 无锡聚为传感科技有限公司 | 一种分布式光纤振动传感器振源定位方法 |
WO2014209141A2 (en) * | 2014-10-30 | 2014-12-31 | Instytut Technik Innowacyjnych Emag | Method and system for assessing a risk of high-energy earth bursts generated by underground mining |
WO2015002558A2 (en) * | 2014-10-30 | 2015-01-08 | Instytut Technik Innowacyjnych Emag | Method and system for analysis of geological structure and relative changes in stress in the layers located above the workings of an underground mine |
CN111665568A (zh) * | 2020-07-09 | 2020-09-15 | 中油奥博(成都)科技有限公司 | 基于分布式光纤声波传感技术的微测井装置及测量方法 |
CN113417699A (zh) * | 2021-08-23 | 2021-09-21 | 辽宁工程技术大学 | 一种煤矿复合灾害分布式光纤动态监测系统及使用方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113985482A (zh) | 2022-01-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103398860B (zh) | 一种基于位移传感器的海洋平台安全监测方法 | |
US20180003550A1 (en) | Gauge length optimization in distributed vibration sensing | |
CA2878584C (en) | Cable comprising twisted sinusoid for use in distributed sensing | |
CN107144380A (zh) | 隧道施工过程中围岩扰动应力场的动态测试方法 | |
CN111442788A (zh) | 一种架空输电线路健康监测方法及系统 | |
CN103727980A (zh) | 一种用于边坡滑坡状况实时监测的光纤传感系统 | |
CN104729777A (zh) | 一种用于相似模拟试验的应力测试装置及方法 | |
CN209166366U (zh) | 一种基于光纤光栅传感技术的新型顶板离层仪 | |
CN105741278A (zh) | 一种基于计算机视觉的拉索分布应力在线监测方法 | |
CN110836690A (zh) | 基于变形量及含水率的高位远程滑坡预警装置及方法 | |
CN208780164U (zh) | 全深度阵列式测斜仪 | |
CN104501732B (zh) | 一种水工结构变形分布式光纤监测装置及方法 | |
CN103364823B (zh) | 震动源实时定位与分析系统 | |
CN111551985A (zh) | 一种基于桩锤激震和k近邻算法的地下溶洞探测方法 | |
CN204854709U (zh) | 一种滑坡的分布监测装置 | |
CN204286374U (zh) | 一种水工结构变形分布式光纤监测装置 | |
CN113985482B (zh) | 基于煤矿井下通讯光缆的矿震震源定位方法 | |
CN211291851U (zh) | 一种盾构隧道土体地震液化实时监测系统 | |
CN210719251U (zh) | 基于变形量及含水率的高位远程滑坡预警装置 | |
CN211147591U (zh) | 一种分布式光纤应变传感系统里程标定装置 | |
CN102346012A (zh) | 一种采空塌陷区油气管道管土相对位移监测方法和系统及系统的构建方法 | |
Karthik et al. | Review on low-cost wireless communication systems for slope stability monitoring in opencast mines | |
CN109323684A (zh) | 一种测斜系统及其测斜方法 | |
CN113639668A (zh) | 一种基于光纤感测技术的角度位移计及其安装方法 | |
CN108982001B (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 |