一种地质勘探数据处理方法及系统

1.本发明涉及地质勘探数据处理技术领域，特别是涉及一种地质勘探数据处理方法及系统。


背景技术：

2.激发极化法是以岩矿石和地下水的激发极化效应差异为基础，用人工源进行激发，以某种极距的装置形式观测地下介质的极化率和视电阻率的变化，从而解决资源和矿产等地质问题的方法。作为地球物理勘探中最悠久和最成熟的方法之一，常规激电法在浅部水资源和矿产资源勘探中发挥着重要的作用，特别是针对含硫化物类的金属矿，有显著的效果。
3.wdjs
‑
3激电仪是目前国内实施激发极化法的主要实用仪器，仪器自带的软件只能进行原始数据的读取，对原始数据的进一步的处理需要人工手动完成，操作过程繁琐，并且一次只能针对一个文件，在数据量较大的时候会耗费大量时间，处理过程有很多重复性的工作，也会因人为因素产生大量错误，因此，迫切需要一种可以进行自动处理原始数据的方法。


技术实现要素：

4.本发明的目的是提供一种地质勘探数据自动处理方法及系统，可以实现数据的自动整理、自动处理和文件自动汇总等功能，为地球物理勘探和分析提供了数据基础，极大地提高了数据处理的效率和准确性。
5.为实现上述目的，本发明提供了如下方案：
6.一种地质勘探数据处理方法，所述处理方法包括：
7.获取激电仪参数、原始数据和原始数据文件；所述激电仪参数包括：距离am、距离an、距离bm和距离bn，所述距离am为第一发射电极a和第一接收电极m之间的距离，所述距离an为所述第一发射电极a和第二接收电极n之间的距离，所述距离bm为第二发射电极b和所述第一接收电极m之间的距离，所述距离bn为所述第二发射电极b和所述第二接收电极n之间的距离；所述原始数据由所述激电仪测得，所述原始数据包括：测点点号、测线号、发射电流和电位差；所述原始数据文件用于存储所述原始数据；所述激电仪参数、所述原始数据和所述原始数据文件数量相同且一一对应；
8.根据所述距离am、所述距离an、所述距离bm、所述距离bn、所述发射电流和所述电位差，利用公式和计算每个测点相应的k值和视电阻率；其中p
s
为视电阻率，v为电位差，i为发射电流；
9.将所述k值和所述视电阻率增加至所述原始数据文件中以得到更新后的原始数据文件；
10.汇总所有所述原始数据以得到原始数据汇总文件；
11.根据质检数据的测点信息与所述原始数据的测点进行匹配，以得到相应测点的k值，并根据所述质检数据相应的发射电流和电位差计算所述质检数据的视电阻率；判断所述质检数据的视电阻率与相应测点的视电阻率的绝对值是否超出第一阈值；若判断结果为否，则按照质检数据的测线号和质检数据的测点点号以从小到大的顺序对所述质检数据进行排列，生成质检数据汇总文件；所述质检数据的测点信息包括：质检数据的测线号和质检数据的测点点号；
12.将所述质检数据汇总文件中的质检数据与所述原始数据汇总文件中的原始数据匹配，并按照需要的格式存储在同一个数据文件里不同的表格中，以并生成质检汇总总文件；
13.分别为所述更新后的原始数据文件、所述汇总质检数据文件和所述质检汇总总文件添加平面坐标以得到相应的带坐标的更新后的原始数据文件、带坐标的汇总质检数据文件和带坐标的质检汇总总文件。
14.可选的，所述处理方法还包括：判断同一测线的测点点号是否单调递或单调递减，以确定是否有重复的测点点号。
15.可选的，所述处理方法还包括：对多个测线号进行单一化处理以提取唯一的线号，得到测线数量；所述测线数量与所述原始数据文件的数量相等。
16.可选的，所述原始数据还包括：距离ab、测程d
测程
、第一端头和第二端头；所述距离ab为所述第一发射电极a和所述第二发射电极b之间的距离，所述测程为所述第一端头和所述第二端头之间的距离；
17.所述处理方法还包括：判断起始点号和终止点号是否超出测点范围；若判断结果为是，则产生“错误”提示；所述测点范围为所述第一端头和所述第二端头之间的范围，并包括所述第一端头和所述第二端头；
18.所述起始点号为p
起始点号
＝(ab
‑
d
测程
)
×
0.5，所述终止点号为p
终止点号
＝p
起始点号
+d
测程
。
19.可选的，所述处理方法还包括：判断相邻测点点号的差值是否相等；若判断结果为否，则产生“数据遗漏”提示。
20.一种地质勘探数据处理系统，所述处理系统包括：
21.数据获取单元，用于获取激电仪参数、原始数据和原始数据文件；所述激电仪参数包括：距离am、距离an、距离bm和距离bn，所述距离am为第一发射电极a和第一接收电极m之间的距离，所述距离an为所述第一发射电极a和第二接收电极n之间的距离，所述距离bm为第二发射电极b和所述第一接收电极m之间的距离，所述距离bn为所述第二发射电极b和所述第二接收电极n之间的距离；所述原始数据由所述激电仪测得，所述原始数据包括：测点点号、测线号、发射电流和电位差；所述原始数据文件用于存储所述原始数据；所述激电仪参数、所述原始数据和所述原始数据文件数量相同且一一对应；
22.计算单元，用于根据所述距离am、所述距离an、所述距离bm、所述距离bn、所述发射电流和所述电位差，利用公式和计算每个测点相应的k值和视电阻率；其中p
s
为视电阻率，v为电位差，i为发射电流；
23.原始数据文件更新单元，用于将所述k值和所述视电阻率增加至所述原始数据文件中以得到更新后的原始数据文件；
24.原始数据汇总单元，用于汇总所有所述原始数据以得到原始数据汇总文件；
25.质检数据汇总文件生成单元，用于根据质检数据的测点信息与所述原始数据的测点进行匹配，以得到相应测点的k值，并根据所述质检数据相应的发射电流和电位差计算所述质检数据的视电阻率；判断所述质检数据的视电阻率与相应测点的视电阻率的绝对值是否超出第一阈值；若判断结果为否，则按照质检数据的测线号和质检数据的测点点号以从小到大的顺序对所述质检数据进行排列，生成质检数据汇总文件；所述质检数据的测点信息包括：质检数据的测线号和质检数据的测点点号；
26.质检汇总总文件生成单元，用于将所述质检数据汇总文件中的质检数据与所述原始数据汇总文件中的原始数据匹配，并按照需要的格式存储在同一个数据文件里不同的表格中，以并生成质检汇总总文件；
27.含坐标系形式文件生成单元，用于分别为所述更新后的原始数据文件、所述汇总质检数据文件和所述质检汇总总文件添加平面坐标以得到相应的带坐标的更新后的原始数据文件、带坐标的汇总质检数据文件和带坐标的质检汇总总文件。
28.可选的，所述处理系统还包括：测点点号判断单元，用于判断同一测线的测点点号是否单调递或单调递减，以确定是否有重复的测点点号。
29.可选的，所述处理系统还包括：测线数量获得单元，用于对多个测线号进行单一化处理以提取唯一的线号，得到测线数量；所述测线数量与所述原始数据文件的数量相等。
30.可选的，所述原始数据还包括：距离ab、测程d
测程
、第一端头和第二端头；所述距离ab为所述第一发射电极a和所述第二发射电极b之间的距离，所述测程为所述第一端头和所述第二端头之间的距离；
31.所述处理系统还包括：测点判断单元，用于判断起始点号和终止点号是否超出测点范围；若判断结果为是，则产生“错误”提示；所述测点范围为所述第一端头和所述第二端头之间的范围，并包括所述第一端头和所述第二端头；
32.所述起始点号为p
起始点号
＝(ab
‑
d
测程
)
×
0.5，所述终止点号为p
终止点号
＝p
起始点号
+d
测程
。
33.可选的，所述处理系统还包括：数据遗漏判断单元，用于判断相邻测点点号的差值是否相等；若判断结果为否，则产生“数据遗漏”提示。
34.根据本发明提供的具体实施例，本发明公开了以下技术效果：
35.(1)本发明实现了数据的自动整理、自动处理和文件自动汇总等功能，避免人为因素干扰带来的误差，为地球物理勘探和分析提供了数据基础，极大地提高了数据处理的效率和准确性。
36.(2)本发明可以对多个数据文件进行批量处理，大大提高了效率。
附图说明
37.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
38.图1为本发明提供的地质勘探数据处理方法流程图；
39.图2为本发明提供的地质勘探数据处理系统示意图。
具体实施方式
40.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
41.本发明的目的是提供一种地质勘探数据处理方法及系统，提高数据处理的效率和准确性。
42.为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
43.图1示出了本发明提供的地质勘探数据处理方法流程图，所述处理方法包括：
44.步骤101：获取激电仪参数、原始数据和原始数据文件；所述激电仪参数包括：距离am、距离an、距离bm和距离bn，所述距离am为第一发射电极a和第一接收电极m之间的距离，所述距离an为所述第一发射电极a和第二接收电极n之间的距离，所述距离bm为第二发射电极b和所述第一接收电极m之间的距离，所述距离bn为所述第二发射电极b和所述第二接收电极n之间的距离；所述原始数据由所述激电仪测得，所述原始数据包括：测点点号、测线号、发射电流和电位差；所述原始数据文件用于存储所述原始数据；所述激电仪参数、所述原始数据和所述原始数据文件数量相同且一一对应。
45.为了保证数据采集准确性，一般要进行1
‑
3次测量再保存数据。在重复的测点处，需标识“重复”。本发明可以自动判断m线的n点，和n+1点是否一样，如果一样，代表n点和n+1点是重复测点，则在n+1点对应列显示“重复”；如果不一样，代表n点和n+1点是不同的测点，则继续搜索n+2点和n+1点是否一样，以此类推。
46.步骤102：根据所述距离am、所述距离an、所述距离bm、所述距离bn、所述发射电流和所述电位差，利用公式和计算每个测点相应的k值和视电阻率；其中p
s
为视电阻率，v为电位差，i为发射电流。
47.步骤103：将所述k值和所述视电阻率增加至所述原始数据文件中以得到更新后的原始数据文件。
48.视电阻率代表了采集目标的导电特性，是评判目标体地球物理特征的重要指标，因此，本发明为地球物理勘探和分析提供了数据基础。
49.步骤104：汇总所有所述原始数据以得到原始数据汇总文件。
50.将所有的原始数据进行汇总，便于统一处理，提高数据处理效率。
51.步骤105：根据质检数据的测点信息与所述原始数据的测点进行匹配，以得到相应测点的k值，并根据所述质检数据相应的发射电流和电位差计算所述质检数据的视电阻率；判断所述质检数据的视电阻率与相应测点的视电阻率的绝对值是否超出第一阈值；若判断结果为否，则按照质检数据的测线号和质检数据的测点点号以从小到大的顺序对所述质检
数据进行排列，生成质检数据汇总文件；所述质检数据的测点信息包括：质检数据的测线号和质检数据的测点点号。
52.质检数据又称质量检查测量数据，用于针对测量数据进行的质量检查，即不同操作人员在不同时间用不同仪器设备在同一个点进行测量，得到质检数据，通过比较质检数据和原始测量数据是否一致，以判断原始数据是否正确，有利于提高数据处理的准确性。
53.步骤106：将所述质检数据汇总文件中的质检数据与所述原始数据汇总文件中的原始数据匹配，并按照需要的格式存储在同一个数据文件里不同的表格中，以并生成质检汇总总文件；
54.步骤107：分别为所述更新后的原始数据文件、所述汇总质检数据文件和所述质检汇总总文件添加平面坐标以得到相应的带坐标的更新后的原始数据文件、带坐标的汇总质检数据文件和带坐标的质检汇总总文件。
55.带坐标的文件便于测点位置查找，以减少时间消耗。
56.利用激电仪进行数据采集的时，测量点是放置发射电极的位置，例如在n点放置第一发射电极a，在n+1点放置第二发射电极b，而实际采集的是相邻测点的中间位置的数据，所以需要将在同一条测线上的n点和n+1点的平面坐标取平均值。当n点位最后一个测点，则不存在n+1点，数据处理过程会报错。因此，在平面坐标数据最后一行，加上了数据元素均为π的“托底数据体”，可在n点终止数据采集，即本发明自动搜索到表示为π的数据时，以n
‑
1点和n点的平均值作为终止坐标。
57.将平面坐标和数据文件的对应测点坐标匹配，最终得到带坐标的数据文件本发明设置了终止标识“π”，有利于在保证准确率的情况下完成数据的自动识别和处理。
58.本发明的一个实施例为，该处理方法还包括：判断同一测线的测点点号是否单调递或单调递减，以确定是否有重复的测点点号。在存在点号输入错误的情况下，输出错误提示，以及时修正测点点号。
59.本发明的一个实施例为，该处理方法还包括：对多个测线号进行单一化处理以提取唯一的线号，得到测线数量；所述测线数量与所述原始数据文件的数量相等。通过对测线号的单一化处理，可以排除重复数据，减小数据量，缩短数据处理时间。
60.本发明的一个实施例为，原始数据还包括：距离ab、测程d
测程
、第一端头和第二端头；所述距离ab为所述第一发射电极a和所述第二发射电极b之间的距离，所述测程为所述第一端头和所述第二端头之间的距离。
61.该处理方法还包括：判断起始点号和终止点号是否超出测点范围；若判断结果为是，则产生“错误”提示；所述测点范围为所述第一端头和所述第二端头之间的范围，并包括所述第一端头和所述第二端头。在起始点号和终止点号超出测点范围时，要舍弃相应的测点数据。
62.所述起始点号为p
起始点号
＝(ab
‑
d
测程
)
×
0.5，所述终止点号为p
终止点号
＝p
起始点号
+d
测程
。
63.本发明的一个实施例为，该处理方法还包括：判断相邻测点点号的差值是否相等；若判断结果为否，则产生“数据遗漏”提示。根据提示，共作人员可以及时补全遗漏的数据，有利于数据处理的准确性。
64.本发明还提供了一种地质勘探数据处理方法相应处理系统，如图2所示，该处理系统包括：数据获取单元201、计算单元202、原始数据文件更新单元203、原始数据汇总单元
204、质检数据汇总文件生成单元205、质检汇总总文件生成单元206和含坐标系形式文件生成单元207。
65.数据获取单元201用于获取激电仪参数、原始数据和原始数据文件；所述激电仪参数包括：距离am、距离an、距离bm和距离bn，所述距离am为第一发射电极a和第一接收电极m之间的距离，所述距离an为所述第一发射电极a和第二接收电极n之间的距离，所述距离bm为第二发射电极b和所述第一接收电极m之间的距离，所述距离bn为所述第二发射电极b和所述第二接收电极n之间的距离；所述原始数据由所述激电仪测得，所述原始数据包括：测点点号、测线号、发射电流和电位差；所述原始数据文件用于存储所述原始数据；所述激电仪参数、所述原始数据和所述原始数据文件数量相同且一一对应；
66.计算单元202用于根据所述距离am、所述距离an、所述距离bm、所述距离bn、所述发射电流和所述电位差，利用公式和计算每个测点相应的k值和视电阻率；其中p
s
为视电阻率，v为电位差，i为发射电流；
67.原始数据文件更新单元203用于将所述k值和所述视电阻率增加至所述原始数据文件中以得到更新后的原始数据文件；
68.原始数据汇总单元204用于汇总所有所述原始数据以得到原始数据汇总文件；
69.质检数据汇总文件生成单元205用于根据质检数据的测点信息与所述原始数据的测点进行匹配，以得到相应测点的k值，并根据所述质检数据相应的发射电流和电位差计算所述质检数据的视电阻率；判断所述质检数据的视电阻率与相应测点的视电阻率的绝对值是否超出第一阈值；若判断结果为否，则按照质检数据的测线号和质检数据的测点点号以从小到大的顺序对所述质检数据进行排列，生成质检数据汇总文件；所述质检数据的测点信息包括：质检数据的测线号和质检数据的测点点号；
70.质检汇总总文件生成单元206用于将所述质检数据汇总文件中的质检数据与所述原始数据汇总文件中的原始数据匹配，并按照需要的格式存储在同一个数据文件里不同的表格中，以并生成质检汇总总文件；
71.含坐标系形式文件生成单元207用于分别为所述更新后的原始数据文件、所述汇总质检数据文件和所述质检汇总总文件添加平面坐标以得到相应的带坐标的更新后的原始数据文件、带坐标的汇总质检数据文件和带坐标的质检汇总总文件。
72.本发明的一个实施例为，该处理系统还包括：用于判断同一测线的测点点号是否单调递或单调递减，以确定是否有重复的测点点号的测点点号判断单元。
73.本发明的一个实施例为，该处理系统还包括：用于对多个测线号进行单一化处理以提取唯一的线号，得到测线数量的测线数量获得单元；所述测线数量与所述原始数据文件的数量相等。
74.本发明的一个实施例为，所述原始数据还包括：距离ab、测程d
测程
、第一端头和第二端头；所述距离ab为所述第一发射电极a和所述第二发射电极b之间的距离，所述测程为所述第一端头和所述第二端头之间的距离。
75.该处理系统还包括：用于判断起始点号和终止点号是否超出测点范围的测点判断单元；若判断结果为是，则产生“错误”提示；所述测点范围为所述第一端头和所述第二端头
之间的范围，并包括所述第一端头和所述第二端头。
76.所述起始点号为p
起始点号
＝(ab
‑
d
测程
)
×
0.5，所述终止点号为p
终止点号
＝p
起始点号
+d
测程
。
77.本发明的一个实施例为，该处理系统还包括：用于判断相邻测点点号的差值是否相等；若判断结果为否，则产生“数据遗漏”提示的数据遗漏判断单元。
78.本发明的处理系统可以实现的技术效果与本发明处理方法实现的技术效果相同。
79.除非本文另有说明或与上下文明显矛盾，否则本文所述的所有方法都能够以任何合适的顺序进行。
80.本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
81.本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。