瞬变电磁数据小波数据的圆滑处理方法
【技术领域】
[0001] 本发明属于地球物理勘探技术领域，涉及一种瞬变电磁数据小波数据的圆滑处理 方法。
【背景技术】
[0002] 瞬变电磁原始数据的圆滑处理，往往使用人为手工进行后期校正，这样由于一些 复杂的原因，使瞬变电磁原始数据在数据处理中可能出现更大的误差，产生多解性问题。在 这种情况下，数据处理工作者关心的瞬变电磁原始数据的圆滑质量得不到较好的效果，因 此，需要研究突出异常的解决办法。
【发明内容】
[0003] 有鉴于此，本发明的主要目的在于提供一种瞬变电磁数据小波数据的圆滑处理方 法。
[0004] 为达到上述目的，本发明的技术方案是这样实现的：
[0005] 本发明实施例提供一种瞬变电磁数据小波数据的圆滑处理方法，该方法为：对原 始瞬变电磁数据进行线性插值，对所述线性插值后获得的数据进行分解并且确定每个测点 在不同测道的斜率，当任意一个测点的斜率的方差值小于等于5时，根据数值定义校正系 数对该测点的原始瞬变电磁数据进行圆滑调节。
[0006] 上述方案中，所述对原始瞬变电磁数据进行线性插值，具体为：对原始瞬变电磁数 据沿X、Y方向按克里格网络插值法插值，插值后得到X、Y方向上需要进行计算的若干条测 线数据。
[0007] 上述方案中，所述对所述线性插值后获得的数据进行分解并且确定每个测点在不 同测道的斜率，具体为：对线性插值后的数据进行小波基分解，分解后得到该数据的双对数 坐标下的曲线值，计算分解后的每个数据不同测道的斜率；
[0008] 沿某一选择方向（如X方向或Y方向，二者选择一个即可）对第i条测线数据采
算衰减曲线在每个测点不同测道的斜率。
[0009] 上述方案中，所述当任意一个测点的斜率的方差值小于等于5时，根据数值定义 校正系数对该测点的原始瞬变电磁数据进行圆滑调节，具体为：计算出各测点斜率的平
测点的斜率的方差值小于等于5时，则该测点的方差值没有超限，根据数值定义校正系数 对该测点的原始瞬变电磁数据进行圆滑调节；反之，当任意一个测点的斜率的方差值大于 5时，则该测点的方差值超限，更换测点进行计算。
[0010] 上述方案中，所述数值定义校正系数分为第一数值定义校正系数和第二数值定义 校正系数；当前时间道时间与前一道时间差值平方的倒数的第一数值定义校正系数，记为 r 1;当前时间道时间与后一道时间差值平方的倒数的第二数值定义校正系数，记为r2，通过
[0011] 与现有技术相比，本发明的有益效果：
[0012] 本发明能够通过用小波变换方法将瞬变电磁原始数据进行计算，通过小波计算这 种数学计算方法对数据进行圆滑，达到消除人为因素的影响，突出数据原有基本特征的目 的。
【附图说明】
[0013] 图1为本发明实施例提供一种瞬变电磁数据小波数据的圆滑处理方法的总流程 图；
[0014] 图2为瞬变电磁圆滑前数据曲线；
[0015] 图3为瞬变电磁圆滑后数据曲线；
[0016] 图4为瞬变电磁数据圆滑前成果图；
[0017] 图5为瞬变电磁数据圆滑后成果图。
【具体实施方式】
[0018] 下面结合附图和【具体实施方式】对本发明进行详细说明。
[0019] 本发明实施例提供一种瞬变电磁数据小波数据的圆滑处理方法，该方法为：对原 始瞬变电磁数据进行线性插值，对所述线性插值后获得的数据进行分解并且确定每个测点 在不同测道的斜率，当任意一个测点的斜率的方差值小于等于5时，根据数值定义校正系 数对该测点的原始瞬变电磁数据进行圆滑调节。
[0020] 所述对原始瞬变电磁数据进行线性插值，具体为：对原始瞬变电磁数据沿X、Y方 向按克里格网络插值法插值，插值后得到X、Y方向上需要进行计算的若干条测线数据。
[0021] 所述对所述线性插值后获得的数据进行分解并且确定每个测点在不同测道的斜 率，具体为：对线性插值后的数据进行小波基分解，分解后得到该数据的双对数坐标下的曲 线值，计算分解后的每个数据不同测道的斜率；
[0022] 沿某一选择方向（如X方向或Y方向，二者选择一个即可）对第i条测线数据采
算衰减曲线在每个测点不同测道的斜率。
[0023] 所述当任意一个测点的斜率的方差值小于等于5时，根据数值定义校正系数对该
其中S2为方差值，X i为斜率值，n为需要计算的斜率的个数；当任意一个测点的斜率的方差 值小于等于5时，则该测点的方差值没有超限，根据数值定义校正系数对该测点的原始瞬 变电磁数据进行圆滑调节；反之，当任意一个测点的斜率的方差值大于5时，则该测点的方 差值超限，更换测点进行计算。
[0024] 所述数值定义校正系数分为第一数值定义校正系数和第二数值定义校正系数；当 前时间道时间与前一道时间差值平方的倒数的第一数值定义校正系数，记为r 1;当前时间 道时间与后一道时间差值平方的倒数的第二数值定义校正系数，记为r2，通过 ri、r2与它们
[0025] 实施例：
[0026] 本发明实施例提供一种瞬变电磁数据小波数据的圆滑处理方法，如图1所示，实 现步骤为：
[0027] (1)对数据沿X、Y方向按克里格网络插值法插值；
[0028] (2)沿选择方向（如X方向）对第i条测线数据（网格插值后数据）采用小波基
压值，t为时间）计算衰减曲线在每个测点不同测道的斜率；
的斜率的个数。
[0031] (5)判断各点方差是否超限，如果超限，更换测点进行计算，方差值要求为< 5的 数值，以此作为圆滑的标准；
[0032] (6)给出两个校正系数ri (当前时间道时间与前一道时间差值平方的倒数）和 r2(当前时间道时间与与后一道时间差值平方的倒数），通过这两个系数与它们和的比值，
[0034] (8)通过测线上测点数是否为最大判断是否对该方向所有测线均完成以上步骤， 如没有，更换测线，重复步骤（1) 一（7);
[0035] (9)通过（7)计算得到的衰减电压值绘图，得到最终结果。
[0036] 下面以实测数据为例说明本发明效果：
[0037] 从图2和图3可以看出，由于一些复杂的地质原因，使图2上的瞬变电磁原始采集 数据干扰很大，难以对数据进行分析处理，该干扰影响了对数据的正确认识和后期处理。经 过小波变换曲线圆滑后，图3反映出的圆滑后瞬变电磁曲线单调衰减，曲线圆滑，这与常规 瞬变电磁曲线，以及相邻点数据曲线匹配较好。
[0038] 从图4和图5对比可以看出，经过小波变换曲线圆滑后，图5反映出的圆滑后瞬变 电磁成果图曲线圆滑，这与常规瞬变电磁成果匹配较好。说明该方法使用是正确且有效的。 经过计算后与实际巷道异常位置吻合很好。
[0039] 以上所述，仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。
【主权项】
1. 一种瞬变电磁数据小波数据的圆滑处理方法，其特征在于，该方法为：对原始瞬变 电磁数据进行线性插值，对所述线性插值后获得的数据进行分解并且确定每个测点在不同 测道的斜率，当任意一个测点的斜率的方差值小于等于5时，根据数值定义校正系数对该 测点的原始瞬变电磁数据进行圆滑调节。2. 根据权利要求1所述的瞬变电磁数据小波数据的圆滑处理方法，其特征在于，所述 对原始瞬变电磁数据进行线性插值，具体为：对原始瞬变电磁数据沿X、Y方向按克里格网 络插值法插值，插值后得到X、Y方向上需要进行计算的若干条测线数据。3. 根据权利要求1或2所述的瞬变电磁数据小波数据的圆滑处理方法，其特征在于，所 述对所述线性插值后获得的数据进行分解并且确定每个测点在不同测道的斜率，具体为： 对线性插值后的数据进行小波基分解，分解后得到该数据的双对数坐标下的曲线值，计算 分解后的每个数据不同测道的斜率； 沿某一选择方向（如X方向或Y方向，二者选择一个即可）对第i条测线数据采用小 波基进行分解，分解方法为'其中Φ (t)为分解后结果，t为时间；在双对数 坐标中，首先通过公式，kS斜率，U为每个测点衰减电压值，t为时间；计算衰 减曲线在每个测点不同测道的斜率。4. 根据权利要求3所述的瞬变电磁数据小波数据的圆滑处理方法，其特征在于，所述 当任意一个测点的斜率的方差值小于等于5时，根据数值定义校正系数对该测点的原始瞬 变电磁数据进行圆滑调节，具体为：计算出各测点斜率的平均值其中X1为斜率 值，η为需要计算的斜率的个数，计算出各点斜率的方差值其中S2为方差 值，Χι为斜率值，η为需要计算的斜率的个数；当任意一个测点的斜率的方差值小于等于5 时，则该测点的方差值没有超限，根据数值定义校正系数对该测点的原始瞬变电磁数据进 行圆滑调节；反之，当任意一个测点的斜率的方差值大于5时，则该测点的方差值超限，更 换测点进行计算。5. 根据权利要求4所述的瞬变电磁数据小波数据的圆滑处理方法，其特征在于：所述 数值定义校正系数分为第一数值定义校正系数和第二数值定义校正系数；当前时间道时间 与前一道时间差值平方的倒数的第一数值定义校正系数，记为r 1;当前时间道时间与后一 道时间差值平方的倒数的第二数值定义校正系数，记为r2，通过^、^与它们和的比值，即?进行圆滑调节。
【专利摘要】本发明公开了一种瞬变电磁数据小波数据的圆滑处理方法，该方法为：对原始瞬变电磁数据进行线性插值，对所述线性插值后获得的数据进行分解并且确定每个测点在不同测道的斜率，当任意一个测点的斜率的方差值小于等于5时，根据数值定义校正系数对该测点的原始瞬变电磁数据进行圆滑调节。本发明能够通过用小波变换方法将瞬变电磁原始数据进行计算，通过小波计算这种数学计算方法对数据进行圆滑，达到消除人为因素的影响，突出数据原有基本特征的目的。
【IPC分类】G01V3/38
【公开号】CN105068129
【申请号】CN201510390409
【发明人】张军
【申请人】中煤科工集团西安研究院有限公司
【公开日】2015年11月18日
【申请日】2015年7月6日