一种小断层超前三维可视化建模方法技术

本发明专利技术公开了一种小断层超前三维可视化建模方法，先建立煤巷掘进超前预报观测系统；然后对各个三分量检波器接收到的三分量Love型反射槽波信号，进行时频域极化偏移成像，获取反射槽波在本盘界面上反射点位置及断层走向；在接收到的三分量信号中确定绕射体波信号，然后对绕射体波信号采用时频域希尔伯特极化成像方法进行处理，获取断层对盘绕射点的位置，进行收敛成像，并获得断层落差；通过上述得到的反射点及绕射点几何位置关系，获取小断层的倾角参数；最后获得断层性质及倾向，根据上述断层参数即可建立小断层三维空间模型。因此本发明专利技术能精确探测出巷道前方断层的三维空间下的各要素信息，从而能实现超前探测小断层的三维可视化。

A 3D visualization modeling method for small fault

The invention discloses a three-dimensional visual modeling method for small fault advance, which first establishes a coal tunnel advance prediction and observation system; then carries out time-frequency polarization migration imaging for the three-component love reflection groove wave signal received by each three-component geophone, obtains the reflection point position and fault direction of the reflection groove wave on the interface of the plate; determines the diffraction in the received three-component signal Body wave signal, and then the diffraction body wave signal is processed by time-frequency Hilbert polarization imaging method, to obtain the location of the diffraction point of the fault to the disk, to perform convergence imaging, and to obtain the fault drop; through the geometric position relationship of the reflection point and the diffraction point obtained above, to obtain the dip angle parameter of the small fault; finally, to obtain the fault property and tendency, according to the above fault parameters can be Three dimensional space model of small fault is established. Therefore, the invention can accurately detect the information of each element under the three-dimensional space of the fault in front of the roadway, so as to realize the three-dimensional visualization of detecting small fault in advance.

【技术实现步骤摘要】
一种小断层超前三维可视化建模方法
本专利技术涉及一种矿井地球物理探测方法，具体是一种小断层超前三维可视化建模方法。
技术介绍
随着煤矿开采深度的不断增加，三维地震探测断层的难度增加，导致目前对断层的探查程度不够，同时煤巷掘进过程中由于断层的影响造成的煤与瓦斯突出事故和突水频发，给煤矿生产带来巨大的安全隐患，因此积极发展井下探查技术是目前煤矿生产迫切需要。煤巷掘进超前探测断层主要通过矿井物探方法(反射槽波超前探、瞬变电磁法，直流电法等)和钻探技术，钻探技术虽然能够直观的反应前方的构造情况，但存在“一孔之径”的问题所导致的探测范围的有限，且钻探成本较高，目前地震反射槽波最适合超前探测断层，目前针对回采工作面内断层的探测，前人做了大量工作，而针对掘进迎头前方断层探测较少。专利申请号“201711007670.6”中杨思通提出厚煤层沿底板或顶板掘进煤巷反射槽波超前探测方法，利用Rayleigh槽波进行超前探测，但Rayleigh槽波形成条件较为苛刻，而Love型槽波更具普遍性。专利申请号“201711391055.X”中王保利提出了基于孔巷联合的井下独头巷道槽波超前探测方法及系统，利用偏振滤波对地震记录进行区噪处理，采用散射成像方法对去噪后的地震记录进行偏移成像；综上所述，针对煤巷掘进迎头前方断层研究上取得了一定的成果，但均没有解决如何可视化探测断层的几何要素(如断层面的具体走向、倾向、倾角等)的问题。因此如何专利技术一种小断层超前三维可视化建模方法是本行业的研究方向。
技术实现思路
针对上述现有技术存在的问题，本专利技术提供一种小断层超前三维可视化建模方法，能精确探测出巷道前方断层，在三维空间中的各要素参数信息，从而能实现超前探测断层的三维可视化。为了实现上述目的，本专利技术采用的技术方案是：一种小断层超前三维可视化建模方法，该方法的具体步骤为：步骤一：在煤巷左右帮的中心高度分别布置n个三分量检波器接收地震信号，道间距(即相邻三分量检波器之间的距离)为N米，从迎头反向依此布置；步骤二：采用CO2炮在煤巷左右帮煤层中各激发一次，分别产生震源S1和震源S2；(震源产生的原理为CO2炮激发后对煤岩体压裂过程中会产生震动)步骤三：沿煤巷建立三维坐标系，X方向为指向迎头，Y方向为垂直于煤壁方向，Z方向为垂直于煤层底板，以模型边界上煤巷中心为原点建立三维空间坐标系，将上述震源位置与三分量检波器位置归入到空间坐标系中，建立三维煤巷地震超前预报观测系统；步骤四：分析震源S1和S2的地震记录，选取信噪比高的地震记录进行处理分析；步骤五：在地震记录中确定来自本盘的反射Love型槽波信号，进行时频域极化偏移成像，获取本盘断层的位置及走向，并获得最大振幅反射点位置；步骤六：在地震记录中确定来自对盘的绕射体波信号，然后对绕射体波信号采用时频域希尔伯特极化成像方法进行处理，获取对盘绕射点的位置，进行收敛成像；步骤七：根据步骤六中对盘绕射点位置求取断层的落差；步骤八：结合步骤五、步骤六中反射点位置和绕射点位置，通过几何位置关系求取断层的倾角、性质、倾向；步骤九：根据步骤五至八获取的断层的位置、走向、倾角、性质及倾向的参数信息，建立断层三维空间模型。进一步，所述步骤五具体步骤为：(1)在由t组成的时窗t1范围内的反射Love型数据Ai(t)采用解析信号ui(t)来表示：式中，表示任意采样点优势频率，argci(t)表示求取Ai(t)的幅角，表示复数共轭；(2)根据上述时频特征，构建协方差矩阵：式中表示t时刻的l分量时窗叠加振幅，第l分量在时窗内的均值μl，Tlh(t)表示t时刻的自适应时窗长度，h为1、2、3；(3)通过上述协方差矩阵求取的最大特征值对应的特征向量，进行归一化后表示为[x1(t),y1(t),z1(t)]；(4)主极化方向的极化优势倾角表示为：其中即主极化方向与XOZ面之间的夹角；(5)根据以上参数进行时频域极化成像，其原理如下所示：其中Ω代表所有地震道，A(t,f)代表标量振幅的叠加函数，Ai代表振幅能量，Ls+Lr代表地震波传播的总路径，v代表地震波速；主极化方向的权重函数Pc表达如下：(6)通过上述(1)至(5)，得到断层本盘界面，计算得出断层界面走向为α1和断层距离坐标原点的水平长度为L；(7)在断层本盘界面上选取最大振幅反射点位置，设其为A点，三维空间坐标为(xi，yi，zi)。进一步，所述步骤六的具体过程为：(1)在地震记录中确定断层对盘断点绕射体波波列，选定时窗t2；(2)从时窗t2中选取负相位最大振幅采样点，对该采样点的三分量信号进行时频域希尔伯特极化分析，获得主极化方向的方位角及倾角参数；(3)通过n个检波器之间的组合计算关系，结合上步骤方位角及倾角参数，获得射线交汇点的空间分布；(4)将上步骤交汇点的位置归入到三维空间坐标系内，然后将三维空间坐标系以1m为单位划分三维空间网格；(5)统计每个三维空间网格内落入的交汇点，将落入交汇点最多的三维空间网格确定为对盘绕射点位置，设其点为B，坐标为(xj，yj，zj)；进一步，所述步骤七的具体过程为：(1)将步骤六中对盘断层绕射点B投影到XOZ面上，即可获得B点坐标为(xj，zj)；(2)根据几何位置关系即可获得落差H为：H＝zj进一步，所述步骤八的具体过程为：(1)由于断层本盘反射点A的坐标为(xi,yi,zi)以及断层对盘绕射点B坐标为(xj,yj,zj)；将A、B两点分别投影到XOZ面上，即坐标分别为(xi，zi)，(xj，zj)，根据坐标关系求取AB两点连线投影角度即：θ即为断层倾角。(2)在XOZ面上，根据几何位置坐标关系，判定断层的性质，具体步骤为：根据上述A、B两点的位置坐标关系，即可获得断层的性质，可分为以下四种情况：①当坐标关系满足zj＞zi，xj＞xi时，为正断层；②当坐标关系满足zj<zi，xj<xi时，为逆断层；③当坐标关系满足zj＜zi，xj＞xi时，为正断层；④当坐标关系满足zj＞zi，xj＜xi时，为逆断层；(3)根据上述步骤(2)得到的断层性质，结合确定的断层走向，即获得断层的倾向。与现有技术相比，本专利技术先建立煤巷掘进超前预报观测系统；然后对各个三分量检波器接收到的三分量Love型反射槽波信号，进行时频域极化偏移成像，获取反射槽波在断层本盘界面上反射点位置及走向、倾向的大小；在接收到的三分量信号中确定绕射体波信号，然后对绕射体波信号进行时频域希尔伯特极化成像方法，获取断层对盘绕射点的位置，进行收敛成像，并获得小断层落差；通过上述得到的反射点及绕射点位置关系，通过几何关系，获取小断层的倾角大小；最后根据获取的断层参数，建立小断层三维空间模型。因此本专利技术能精确本文档来自技高网...

【技术保护点】
1.一种小断层超前三维可视化建模方法，其特征在于，该方法的具体步骤为：/n步骤一：在煤巷左右帮的中心高度分别布置n个三分量检波器接收地震信号，道间距为N米，从迎头反向依此布置；/n步骤二：采用CO

【技术特征摘要】
1.一种小断层超前三维可视化建模方法，其特征在于，该方法的具体步骤为：
步骤一：在煤巷左右帮的中心高度分别布置n个三分量检波器接收地震信号，道间距为N米，从迎头反向依此布置；
步骤二：采用CO2炮在煤巷左右帮煤层中各激发一次，分别产生震源S1和震源S2；
步骤三：沿煤巷建立三维坐标系，X方向为指向迎头，Y方向为垂直于煤壁方向，Z方向为垂直于煤层底板，以模型边界上煤巷中心为原点建立三维空间坐标系，将上述震源位置与三分量检波器位置归入到空间坐标系中，建立三维煤巷地震超前预报观测系统；
步骤四：分析震源S1和S2的地震记录，选取信噪比高的地震记录进行处理分析；
步骤五：在地震记录中确定来自本盘的反射Love型槽波信号，进行时频域极化偏移成像，获取本盘断层的位置及走向，并获得最大振幅反射点位置；
步骤六：在地震记录中确定来自对盘的绕射体波信号，然后对绕射体波信号采用时频域希尔伯特极化成像方法进行处理，获取对盘绕射点的位置，进行收敛成像；
步骤七：根据步骤六中对盘绕射点位置求取断层的落差；
步骤八：结合步骤五、步骤六中反射点位置和绕射点位置，通过几何位置关系求取断层的倾角、性质、倾向；
步骤九：根据步骤五至八获取的断层的位置、走向、倾角、性质及倾向的参数信息，建立断层三维空间模型。


2.根据权利要求1所述的一种小断层超前三维可视化建模方法，其特征在于，所述步骤五具体步骤为：
(1)在由t组成的时窗t1范围内的反射Love型数据Ai(t)采用解析信号ui(t)来表示：



式中，表示任意采样点优势频率，argci(t)表示求取Ai(t)的幅角，表示复数共轭；
(2)根据上述时频特征，构建协方差矩阵：






式中表示t时刻的l分量时窗叠加振幅，第l分量在时窗内的均值μl，Tlh(t)表示t时刻的自适应时窗长度，h为1、2、3；
(3)通过上述协方差矩阵求取的最大特征值对应的特征向量，进行归一化后表示为[x1(t),y1(t),z1(t)]；
(4)主极化方向的极化优势倾角表示为：



其中指主极化方向与XOZ面之间夹角；
(5)根据以上参数进行时频域极化成像，其原理如下所示：



其中Ω代表所有地震道，A(t,f)代表标量振幅的叠加函数，Ai代表振幅能量，Ls+Lr代表地震波传播的总路径，...

【专利技术属性】
技术研发人员：王勃，孙华超，刘盛东，张恒，李兴兴，邢世雨，龚震，王志丹，
申请(专利权)人：中国矿业大学，
类型：发明
国别省市：江苏;32