视电阻率深度快速计算方法技术

本发明专利技术涉及一种视电阻率深度快速计算方法，包括：a.选取已有的航空瞬变电磁成果数据库，对dB/dt数组数据列进行减道和限值；b.在数据库中置入时间数组列和深度系数数组列；c.基于视电阻率计算公式在数据库中计算出对应的视电阻率；d.对步骤c得到的视电阻率数据进行限值；e.基于深度计算公式在数据库中计算出对应的深度；f将上述步骤编写为gs文件，并在Oasismontaj的数据库中调用该gs文件，得到包含视电阻率深度数据的数据库。本发明专利技术的所有步骤均在Oasismontaj软件的数据库中完成，不需要将数据导出到其他软件中进行计算然后再将计算结果导入原数据库中的过程，克服了不同软件之间接口的问题，简化了处理的步骤，提高了现场数据处理效率。现场数据处理效率。现场数据处理效率。

【技术实现步骤摘要】
视电阻率深度快速计算方法


[0001]本专利技术涉及一种视电阻率计算方法，具体地说是一种视电阻率深度快速计算方法。

技术介绍

[0002]航空电磁成像是将观测数据转换为表征地下介质电性分布特征的中间参数，如视电阻率、视深度等。成像算法速度快，能从海量航空电磁数据中快速提取地下电性主要信息，适用于现场快速数据处理，也可为复杂的航空电磁反演提供初始模型。国内航空电磁系统多数还处于研发试验阶段，未有大量的实测数据。虽然电阻率深度成像方法多、成像理论较为成熟，但是基本上是基于模拟数据开展的。常规的做法是将模拟或实测航空电磁数据导入到商业软件或自编的软件中进行计算，再将计算好的视电阻率深度数据导出到专业的数据处理和成图软件中绘图。常规的视电阻率计算方法需要将数据在多个软件中相互转换，这会带来一定的误差，也使得处理流程相对繁琐。
[0003]Oasismontaj软件是一种用于地球科学成图及处理的软件，是航空物探领域常用的软件，航空物探数据一般都在Oasismontaj软件中进行处理。通过引进并使用先进、成熟的航空电磁系统，能够快速获得大量的实测数据，但如何在勘测现场快速开展视电阻率深度计算，将获得的实测数据快速地转换成所需的探测成果，是现在亟需解决的问题。

技术实现思路

[0004]本专利技术的目的就是提供一种视电阻率深度快速计算方法，以解决现有对获得的航空物探数据进行数据处理、视电阻率深度计算流程复杂的问题。
[0005]本专利技术是这样实现的：一种视电阻率深度快速计算方法，包括以下步骤。r/>[0006]a.选取已有的航空瞬变电磁成果数据库，数据库中包含处理后的dB/dt数组数据列，对dB/dt数组数据列进行减道和限值，去除不参与成像和反演的数据以及噪声数据。
[0007]b.在数据库中置入与采样时间道一一对应的衰减采样时间数组列和深度系数数组列。
[0008]c.对步骤a得到的dB/dt数组数据列和步骤b中的衰减采样时间数组列，基于视电阻率计算公式在数据库中计算出对应的视电阻率。
[0009]d.对步骤c得到的视电阻率数据进行限值。
[0010]e.对步骤b中的衰减采样时间数组列和深度系数数组列以及步骤d中经过限值的视电阻率数据，基于深度计算公式在数据库中计算出对应的深度。
[0011]f将上述步骤编写为Oasismontaj软件的gs文件，并在Oasismontaj的数据库中调用该gs文件，得到包含视电阻率深度数据的数据库。
[0012]对于本专利技术的视电阻率深度快速计算方法，在步骤a中，将dB/dt数组数据列中用于计算噪声水平的最后一道去除，完成对数组数据列的减道，然后取dB/dt≥0.001pV/(A
·
m4)的数据，完成对数组数据列的限值。
[0013]对于本专利技术的视电阻率深度快速计算方法，在步骤b中，衰减采样时间数组列为航空电磁系统设置的固定值，深度系数数组列与采样时间道数存在如下三次多项式关系：k_depth=
‑
0.0043n3+0.567n2+7.8469n+213.32式中，k_depth为深度系数数组列；n为采样时间道数。
[0014]对于本专利技术的视电阻率深度快速计算方法，在步骤c中，视电阻率计算公式为：res=6.32
×
0.00001/(SF2
2/3
×
time
5/3
)式中，res为视电阻率；SF2为经过减道、限值后的dB/dt数组数据列；time为衰减采样时间数组列。
[0015]对于本专利技术的视电阻率深度快速计算方法，在步骤d中，对步骤c得到的视电阻率取≤8000Ω
·
m的数据，以进行视电阻率数据的限值。
[0016]对于本专利技术的视电阻率深度快速计算方法，在步骤e中，深度计算公式为：depth=k_depth
×
res1
1/2
×
time
1/2
式中，k_depth为深度系数数组列；res1为经过限值后的视电阻率；time为衰减采样时间数组列。
[0017]对于本专利技术的视电阻率深度快速计算方法，在步骤f中，gs文件通过调用编译好的gx插件来实现不同功能。
[0018]本专利技术通过编译视电阻率深度快速计算的gs文件，并在Oasismontaj数据库中调用该gs文件，实现在同一个软件中完成数据处理、视电阻率深度快速计算及成像。本专利技术的所有步骤均在Oasismontaj软件的数据库中完成，不需要将数据导出到其他软件中进行计算然后再将计算结果导入原数据库中的过程，克服了不同软件之间接口的问题，简化了处理的步骤，提高了现场数据处理效率。通过对航空瞬变电磁数据进行快速视电阻率深度计算并成像，能够在最短的时间内获得测区地质体电阻率信息，为地质解释提供了基础资料，具有广泛的应用前景，适于地质环境调查推广应用。
附图说明
[0019]图1是本专利技术的流程图。
[0020]图2为本专利技术对dB/dt(SF)数据减道限值gs文件内容。
[0021]图3为本专利技术置入时间time列的gs文件内容。
[0022]图4为本专利技术置入深度系数k_depth列的gs文件内容。
[0023]图5为本专利技术计算视电阻率res列的gs文件内容。
[0024]图6为本专利技术对计算视电阻率res列限值的gs文件内容。
[0025]图7为本专利技术计算深度depth列的gs文件内容。
[0026]图8为调用gs文件的示意图。
[0027]图9为通过gs文件计算出视电阻率
‑
深度等的数据库。
[0028]图10为根据计算数据绘制的电阻率深度成像断面等值图。
具体实施方式
[0029]如图1所示，本专利技术包括以下步骤。
[0030]a.选取已有的航空瞬变电磁成果数据库，数据库中包含处理后的dB/dt数组数据
列，对dB/dt数组数据列进行减道和限值，去除不参与成像和反演的数据以及噪声数据。
[0031]b.在数据库中置入与采样时间道一一对应的衰减采样时间数组列和深度系数数组列。
[0032]c.对步骤a得到的dB/dt数组数据列和步骤b中的衰减采样时间数组列，基于视电阻率计算公式在数据库中计算出对应的视电阻率。
[0033]d.对步骤c得到的视电阻率数据进行限值。
[0034]e.对步骤b中的衰减采样时间数组列和深度系数数组列以及步骤d中经过限值的视电阻率数据，基于深度计算公式在数据库中计算出对应的深度。
[0035]f将上述步骤编写为Oasismontaj软件的gs文件，并在Oasismontaj的数据库中调用该gs文件，得到包含视电阻率深度数据的数据库。
[0036]下面基于Oasismontaj软件的gs文件对各步骤进行详细说明。
[0037]步骤1：选取已有的航空瞬变电磁成果数据库，数据库中包含处理本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种视电阻率深度快速计算方法，其特征在于，包括以下步骤：a.选取已有的航空瞬变电磁成果数据库，数据库中包含处理后的dB/dt数组数据列，对dB/dt数组数据列进行减道和限值，去除不参与成像和反演的数据以及噪声数据；b.在数据库中置入与采样时间道一一对应的衰减采样时间数组列和深度系数数组列；c.对步骤a得到的dB/dt数组数据列和步骤b中的衰减采样时间数组列，基于视电阻率计算公式在数据库中计算出对应的视电阻率；d.对步骤c得到的视电阻率数据进行限值；e.对步骤b中的衰减采样时间数组列和深度系数数组列以及步骤d中经过限值的视电阻率数据，基于深度计算公式在数据库中计算出对应的深度；f将上述步骤编写为Oasismontaj软件的gs文件，并在Oasismontaj的数据库中调用该gs文件，得到包含视电阻率深度数据的数据库。2.根据权利要求1所述的视电阻率深度快速计算方法，其特征在于，在步骤a中，将dB/dt数组数据列中用于计算噪声水平的最后一道去除，完成对数组数据列的减道，然后取dB/dt≥0.001pV/(A
·
m4)的数据，完成对数组数据列的限值。3.根据权利要求1所述的视电阻率深度快速计算方法，其特征在于，在步骤b中，衰减采样时间数组列为航空电磁系统设置的固定值，深度系数数组列与采样时间道数存在如下三次多项式关系：k_dep...

【专利技术属性】
技术研发人员：骆燕，陈江源，彭莉红，孔志召，刘波，赵丛，程莎莎，王勇刚，杨波，张伟盟，陈伟，
申请(专利权)人：河北航遥科技有限公司，
类型：发明
国别省市：