基于有限差分法的电磁波波场中瓦斯隧道掌子面空腔及边界条件的方程改进方法技术

本发明专利技术涉及基于有限差分法的电磁波波场中瓦斯隧道掌子面空腔及边界条件的方程改进方法，属隧道工程领域。本发明专利技术先根据探地雷达下正演的相关知识以及真实场景，运用非均匀网格技术，求解案例当前情况下的有限差分方程；其次引入加权平均值的概念，对网格中的各点周围的电导率进行面积加权平均，得到离散化后有限差分算法各个节点电导率值的抽象表达式；再次，监测此时矿井磁场中的各点电导率的稳定性，根据极值来判定监测的时间长短；然后，运用中心差分近似法对上方结果计算差分网度以及最大网度值；最后，根据巷道内满足的拉普拉斯方程对整个矿井巷道空间求体积分，运用时间线性插值进行最后的整理，即得到矿井巷道内可控制边界的相应所需方程。制边界的相应所需方程。制边界的相应所需方程。

【技术实现步骤摘要】
基于有限差分法的电磁波波场中瓦斯隧道掌子面空腔及边界条件的方程改进方法


[0001]本专利技术涉及基于有限差分法的电磁波波场中瓦斯隧道掌子面空腔及边界条件的方程改进方法，尤其涉及到时间限制以及有限采掘边界长度的电磁波强度限制的解决方法，属于隧道工程


技术介绍

[0002]传统的电磁波波场中，瓦斯隧道掌子面空腔以及异构体的断层结构只要传统的运用矿井直流电法勘探、矿井瞬变电磁法勘探、探底雷达和地震勘探类方法。
[0003]在本专利技术中主要依托于有限差分法的电磁波波场中瓦斯隧道掌子面空腔及边界条件进行研究，电磁波在地下岩层中传播时，由于各种岩石、矿电性的不同，它们所呈现出来的电磁波波场有一定的差异，且不同的电磁波波场对岩层的识别存在一定的限制，这就需要考虑到在探测过程中电磁波波场对巷道空腔的边界问题，通过改进有限差分方程，可以改进出有边界条件的巷道空腔电磁波探测方程。但是当电磁波在地下岩层传播时，如果存在着含水地段、陷落柱、断层、空洞或其它不均匀的地质构造，电磁波的波场能量就会被其吸收或者完全屏蔽，此时就会产生信号显著性异常，透视异常。
[0004]对于之前的有限差分相关算法，大多数学者都是基于其稳定性和精确性进行研究，少有考虑到精细网格细分改进差分方程的问题，所以本专利技术针对算法入手，对有限差分方程进行时间线性插值，以此实现减小误差的效果。此外，在重建的过程中，会出现电磁波之间不是平行束的情况，此时应考虑角度问题，本专利技术将运用对算法进行改进，以最大程度达到真值，实现对岩层的勘探实测目的。

技术实现思路

[0005]本专利技术提供了基于有限差分法的电磁波波场中瓦斯隧道掌子面空腔及边界条件的方程改进方法，以用于解决矿井探测中使用的有效差分法中的网格化不够精细、时间长，不能稳定的解出有限采掘边界长度的电磁波波场强度问题。
[0006]本专利技术的技术方案是：基于有限差分法的电磁波波场中瓦斯隧道掌子面空腔及边界条件的方程改进方法，所述方法的具体步骤如下：
[0007]Step1、根据探地雷达下正演的相关知识以及真实场景，运用非均匀网格技术，求解案例当前情况下的有限差分方程；
[0008]首先，根据背景知识，地质雷达主要通过发射天线M向探测的目标发射低频或高频电磁波，其通过高频电磁波在介质中的反射和折射来实现对地下介质的探测。在探测过程中，若探测到电性差异较大的界面时，会发生折射或投射等现象，之后被接天线N接收返回雷达主体，即地质雷达是通过雷达得到的电磁波信号分析波的形状、强度等特征，从而可以判断出地质体的结构。
[0009]本专利技术基于如图4所示的地质雷达探测图，对于隧道建设中出现的空洞进行质量
检测。电阻率、介电常数、磁导率等是影响雷达波在地下介质中传播的电性参数。如果在隧道中出现空洞，空洞内一般为空气，而空气(相对介电常数为1))的介电常数等电性参数与衬砌混凝土(相对介电常数6～10)有很大差异，电磁波在衬砌内传播时，会在空洞上下两个界面产生两次强反射，雷达剖面上会出现双曲线形态的强反射波，其同相轴与相邻道会发生错位。因此能够利用波形特征预计波场的分布情况来识别隧道中的空洞。
[0010]需要运用数值计算的方法对整个空间进行非均匀网格技术分割，将无限大场的分布空间截断为有限的计算空间，在分割的过程中主要采用的是引入合适的边界条件，使得电磁场扩散至边界时能没有明显的反射，运用数学公式得到每一个分割点的有限差分方程如下：
[0011][0012]其中σ表示介质电导率；σ
a
为异常体电导率；j
m
为磁流密度矢量；h
p
，h
s
为磁场强度；μ为介质磁导率。
[0013]Step2、引入加权平均值的概念，对网格中的各点周围的电导率进行面积加权平均，得到离散化后有限差分算法各个节点的电磁波场强的数学表达式；
[0014]所述Step2中包括的各个节点的电磁波场强的数学表达式如下：
[0015][0016]其中，n表示当前时刻，n+1表示下一时刻，i+1就是节点(i,j,k)x轴正方向上多一个单位的数学表示，i
‑
1就是节点(i,j,k)x轴负方向上多一个单位的数学表示，j+1就是节点(i,j,k)y轴正方向上多一个单位的数学表示，j
‑
1就是节点(i,j,k)y轴负方向上多一个单位的数学表示，k+1就是节点(i,j,k)z轴正方向上多一个单位的数学表示，k
‑
1就是节点(i,j,k)z轴负方向上多一个单位的数学表示，为可探测区域总场t＝(n+1)Δt时的磁场内中心节点的场强扩散方程；是x，y，z方向的平均网度；为x轴方向的网度；为y轴方向的网度；为z轴方向的网度；为节点(i,j,k)周围电导率的面积加权平均值；I(t)为关断电流函数；S为所在节点处双侧封闭曲面面积，μ为介质磁导率；
[0017]具体的操作步骤为：
[0018]给有限差分方程的两边在ΔV内积分，得到
[0019][0020]接下来，根据高斯公式对上述式子进行计算并整理，有
[0021][0022]其中ΔV表示前后两个节点移动形成的体积差；S
x
、S
y
和S
z
分别为ΔV与各个方向轴正割的曲面正面；S
‑
x
、S
‑
y
和S
‑
z
分别为ΔV与各个方向轴正割的曲面反面，σ表示介质电导率；σ
a
为异常体电导率；j
m
为磁流密度矢量；h
p
，h
s
为磁场强度；接下来对各积分得到积分结果后，进行面积加权平均计算，得到：
[0023][0024]为节点(i,j,k)周围电导率的面积加权平均，即为
[0025][0026]其中i,j,k分别表示x,y,z方向上的坐标表示；σ
i,j,k
为节点(i,j,k)周围电导率的面积加权值；σ
i+1,j,k
为节点(i+1,j,k)周围电导率的面积加权值，同理，变换下标就求得不同节点的面积加权值；Δx
i
为x轴上的网格距，Δy
j
为y轴上的网格距，Δz
k
为z轴上的网格距；
[0027]最后运用时间线性插值算法，得到各个节点电磁波场强数学表达式。
[0028]Step3、监测此时矿井磁场中的各点电导率的稳定性，根据极值来判定监测的时间长短；
[0029]作为本专利技术的进一步方案，所述Step3中包括具体如下：运用局部细分网格的思路找到时间步长与网格之间的关系来得到趋于稳定监测的时间区域，用时间插值公式表示为：
[0030][0031]其中，h为与全区域中网格不重合的局部网格内部强度，H为与局部网格重合的网格中以及全区域网格中的强度；v为时间步长；m为分的网格数量。
[0032]用本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于有限差分法的电磁波波场中瓦斯隧道掌子面空腔及边界条件的方程改进方法，其特征在于：所述方法的具体步骤如下：Step1、根据探地雷达下正演的相关知识以及真实场景，运用非均匀网格技术，求解案例当前情况下的有限差分方程；Step2、引入加权平均值的概念，对网格中的各点周围的电导率进行面积加权平均，得到离散化后有限差分算法各个节点的电磁波场强的数学表达式；Step3、监测此时矿井磁场中的各点电导率的稳定性，根据极值来判定监测的时间长短；Step4、运用中心差分近似法对Step2中的结果计算差分网度以及最大网度值；Step5、根据巷道内满足的拉普拉斯方程，对整个矿井巷道空间求体积分，运用时间线性插值进行最后的整理，即可以得到矿井巷道内可控制边界的相应所需方程。2.根据权利要求1所述的基于有限差分法的电磁波波场中瓦斯隧道掌子面空腔及边界条件的方程改进方法，其特征在于：所述Step1中包括：需要运用数值计算的方法对整个空间进行非均匀网格技术分割，将无限大场的分布空间截断为有限的计算空间，在分割的过程中采用的是引入合适的边界条件，使得电磁场扩散至边界时能没有明显的反射，运用数学公式得到每一个分割点的有限差分方程如下：其中σ表示介质电导率；σ
a
为异常体电导率；j
m
为磁流密度矢量；h
p
，h
s
为磁场强度；μ为介质磁导率。3.根据权利要求1所述的基于有限差分法的电磁波波场中瓦斯隧道掌子面空腔及边界条件的方程改进方法，其特征在于：所述Step2中包括的各个节点的电磁波场强的数学表达式如下：其中，n表示当前时刻，n+1表示下一时刻，i+1就是节点(i,j,k)x轴正方向上多一个单位的数学表示，i
‑
1就是节点(i,j,k)x轴负方向上多一个单位的数学表示，j+1就是节点(i,j,k)y轴正方向上多一个单位的数学表示，j
‑
1就是节点(i,j,k)y轴负方向上多一个单位的数学表示，k+1就是节点(i,j,k)z轴正方向上多一个单位的数学表示，k
‑
1就是节点(i,j,k)z轴负方向上多一个单位的数学表示，为可探测区域总场t＝(n+1)Δt时的磁场内中心节点的场强扩散方程；是x，y，z方向的平均网度；为x轴方向的网度；为y轴方向的网度；为z轴方向的网度；为节点(i,j,k)周围电导率的面
积加权平均值；I(t)为关断电流函数；S为所在节点处双侧封闭曲面面积，μ为介质磁导率；具体的操作步骤为：给有限差分方程的两边在ΔV内积分，得到接下来，根据高斯公式对上述式子进行计算并整理，有其中ΔV表示前后两个节点移动形成的体积差；S
x
、S
y
和S
z
分...

【专利技术属性】
技术研发人员：郑万波，董锦晓，谭春琳，冉丝蔓，朱榕，李磊，王耀，金正浩，钟佩玲，
申请(专利权)人：云南卫士盾科技有限公司，
类型：发明
国别省市：