一种基于范围约束的层析反演方法技术

本发明专利技术公开了一种基于范围约束的层析反演方法，先采用无线电波坑透获取工作面的原始数据；将工作面进行网格离散化，每个网格作为一个像素，利用直射线理论获得每条射线在工作面内所有网格的截距；根据原始无线电波坑透数据、网格截距以及直射线理论，构建无线电波层析反演矩阵方程组；通过已知的反演方法结合反演矩阵方程组，求得每个像素的反演吸收系数值；最后将已揭露的多源先验数据作为范围约束的数据引入到层析反演计算过程中，不断修正得出的每个像素的反演吸收系数值，最终将修正后的数据通过图像显示，从而得出工作面内的异常分布情况。该方法能提高反演得出的工作面内异常分布情况的准确性。

【技术实现步骤摘要】
一种基于范围约束的层析反演方法
本专利技术涉及一种基于范围约束的层析反演方法，属于煤矿工作面内构造反演

技术介绍
煤炭是我国的主体能源和基础产业，而煤矿工作面内部断层、陷落柱等地质异常的存在直接影响着工作面回采方案制定、回采效率以及井下工作人员的生命安全。以后随着煤矿开采深度的增加，所面临的开采条件将变得更加复杂。因此，对工作面内地质异常构造的精准反演与预测对煤矿的安全高效生产具有十分重要的现实意义。井下无线电波坑透技术作为一种非接触式探测方法，由于其施工便捷、效率高、分辨能力强、成本较低，而且还不存在破坏性等优点，使其成为目前国内外煤矿工作面内部地质异常探测和反演的主要方法。在通过无线电波坑透技术对工作面内地质异常反演结果进行分析解释的过程中，往往还需要结合巷探、钻探甚至回采揭露等多源数据进行综合判断。然而，这些已揭露的多源先验数据只是用于辅助结果的定性分析，并未直接定量参与到层析反演的计算过程中。不仅对预测结果的准确性产生一定的影响，同时还产生一定的数据资源浪费。
技术实现思路
针对上述现有技术存在的问题，本专利技术提供一种基于范围约束的层析反演方法，将已揭露的多源先验数据作为约束引入到层析反演计算过程中，不断修正反演计算过程，从而能提高反演得出的工作面内异常分布情况的准确性。为了实现上述目的，本专利技术采用的技术方案是：一种基于范围约束的层析反演方法，具体步骤为：Ⅰ、采用无线电波坑透获取工作面的原始数据；Ⅱ、将工作面进行网格离散化，利用直射线理论获得每条射线在工作面内所有网格的截距；其中，网格数目为B，射线条数为A；每个网格称为一个像素，并用xj代表第j个网格的真实吸收系数，j＝1,2,3…B；假设任一条射线yi从发射端穿过工作面达到接收端，该射线在第j个网格的截距为dij。其中若设定射线的条数为A，则在第i条射线路径上有：i＝1,2,3…A；由于已知：Hi＝H0-20yilge-20lgri(2)将公式(1)和(2)联立，即：上式中：Hi为针对第i条射线采集到的无线电波场强数据；H0为无线电波的发射初始场强值；ri为第i条射线的长度，初始通过计算得出的已知值；Ⅲ、根据原始无线电波坑透数据、网格截距以及直射线理论，构建无线电波层析反演矩阵方程组；由于煤矿无线电波坑透探测的整个过程是利用多个发射点产生无线电波，其中针对每一个发射点激发产生的无线电波均有多个接收点进行数据的采集，则公式(1)能进一步推广为：简化后，得出：[D][X]＝[Y]式中：[D]为A*B阶系数矩阵，即每条射线通过每个网格的截距，其中射线条数为A，网格数目为B；[X]为B*1阶未知系数矩阵，为每个像素的无线电波吸收系数；[Y]为A*1阶常数矩阵，各观测方式下与实测场强有关的常数，均为已知数据；其中，[X]为需要反演的目标值，即根据结合上述公式采用已知的反演方法求得每个像素的反演吸收系数值x'j；Ⅳ、在通过无线电波坑透技术对工作面内传播介质的吸收系数反演过程中，对得到的每个网格的反演吸收系数值x'j进行约束处理：①在煤矿工作面已揭露无地质异常区域范围内，若x'j小于等于当前反演得到的所有网格的平均吸收系数值xavg，则x'j并不改变；若x'j大于xavg，则通过在反演变量的下限lb与xavg之间随机产生一个数值来代替x'j；②在煤矿工作面已揭露存在地质异常区域范围内，若x'j大于等于xavg，则x'j并不发生改变；若x'j小于xavg，则通过在反演变量的上限ub与xavg之间随机产生一个数值来代替x'j；具体公式为：lb≤x'j≤ub其中：x'j为工作面第j个网格对应的反演吸收系数值；xavg为当前反演的所有网格吸收系数平均值；Q1为在煤矿工作面已揭露无地质异常区域范围内大于xavg的网格序号集合；Q2为在煤矿工作面已揭露存在地质异常区域范围内小于xavg的网格序号集合。ub为反演变量的上限，根据历史数据和多次测试得出的吸收系数值进行设定；lb为反演变量的下限，根据历史数据和多次测试得出的吸收系数值进行设定；δ为0至1之间的随机数；通过上述范围约束对反演得到的每个像素的吸收系数值进行处理，最终将处理后的数据通过图像显示，从而得出工作面内的异常分布情况。与现有技术相比，本专利技术先采用无线电波坑透获取工作面的原始数据；将工作面进行网格离散化，每个网格作为一个像素，利用直射线理论获得每条射线在工作面内所有网格的截距；根据原始无线电波坑透数据、网格截距以及直射线理论，构建无线电波层析反演矩阵方程组；通过已知的反演方法结合反演矩阵方程组，求得每个像素的反演吸收系数值；最后将已揭露的多源先验数据作为范围约束的数据引入到层析反演计算过程中，不断修正得出的每个像素的反演吸收系数值，最终将修正后的数据通过图像显示，从而得出工作面内的异常分布情况。因此本专利技术能提高反演得出的工作面内异常分布情况的准确性。附图说明图1是本专利技术中无线电波层析成像的网格离散化示意图。具体实施方式下面将对本专利技术作进一步说明。如图1所示，本专利技术的具体步骤为：Ⅰ、采用无线电波坑透获取工作面的原始数据；Ⅱ、将工作面进行网格离散化，利用直射线理论获得每条射线在工作面内所有网格的截距；其中，网格数目为B，射线条数为A；每个网格称为一个像素，并用xj代表第j个网格的真实吸收系数，j＝1,2,3…B；假设任一条射线yi从发射端穿过工作面达到接收端，该射线在第j个网格的截距为dij。其中若设定射线的条数为A，则在第i条射线路径上有：i＝1,2,3…A；由于已知：Hi＝H0-20yilge-20lgri(2)将公式(1)和(2)联立，即：上式中：Hi为针对第i条射线采集到的无线电波场强数据；H0为无线电波的发射初始场强值；ri为第i条射线的长度，初始通过计算得出的已知值；Ⅲ、根据原始无线电波坑透数据、网格截距以及直射线理论，构建无线电波层析反演矩阵方程组；由于煤矿无线电波坑透探测的整个过程是利用多个发射点产生无线电波，其中针对每一个发射点激发产生的无线电波均有多个接收点进行数据的采集，则公式(1)能进一步推广为：简化后，得出：[D][X]＝[Y]式中：[D]为A*B阶系数矩阵，即每条射线通过每个网格的截距，其中射线条数为A，网格数目为B；[X]为B*1阶未知系数矩阵，为每个像素的无线电波吸收系数；[Y]为A*1阶常数矩阵，各观测方式下与实测场强有关的常数，均为已知数据；其中，[X]为需要反演的目标值，即根据结合上述公式采用已知的反演方法求得每个像素的反演吸收系本文档来自技高网...

【技术保护点】
1.一种基于范围约束的层析反演方法，其特征在于，具体步骤为：/nⅠ、采用无线电波坑透获取工作面的原始数据；/nⅡ、将工作面进行网格离散化，利用直射线理论获得每条射线在工作面内所有网格的截距；其中，网格数目为B，射线条数为A；每个网格称为一个像素，并用x

【技术特征摘要】
1.一种基于范围约束的层析反演方法，其特征在于，具体步骤为：
Ⅰ、采用无线电波坑透获取工作面的原始数据；
Ⅱ、将工作面进行网格离散化，利用直射线理论获得每条射线在工作面内所有网格的截距；其中，网格数目为B，射线条数为A；每个网格称为一个像素，并用xj代表第j个网格的真实吸收系数，j＝1,2,3…B；
假设任一条射线yi从发射端穿过工作面达到接收端，该射线在第j个网格的截距为dij；其中若设定射线的条数为A，则在第i条射线路径上有：i＝1,2,3…A；



由于已知：
Hi＝H0-20yilge-20lgri(2)
将公式(1)和(2)联立，即：



上式中：
Hi为针对第i条射线采集到的无线电波场强数据；
H0为无线电波的发射初始场强值；
ri为第i条射线的长度；
Ⅲ、根据原始无线电波坑透数据、网格截距以及直射线理论，构建无线电波层析反演矩阵方程组；
由于煤矿无线电波坑透探测的整个过程是利用多个发射点产生无线电波，其中针对每一个发射点激发产生的无线电波均有多个接收点进行数据的采集，则公式(1)能进一步推广为：



简化后，得出：
[D][X]＝[Y]
式中：
[D]为A*B阶系数矩阵，即每条射线通过每个网格的截距，其中射线条数为A，网格数目为B；
[X]为B*1阶未知系数矩阵，为每个像素的无线电波吸收系数；
[Y]为A*1阶常数矩阵，各观测方式下与实测场强有关的常数，均为...

【专利技术属性】
技术研发人员：郭昌放，武祥，
申请(专利权)人：徐州三叉戟信息科技有限公司，
类型：发明
国别省市：江苏;32