Electrical method data mapping system, method and software based on cross-section and plane linkage

Abstract

The present application belongs to the field of electrical exploration, and relates to an electrical data mapping system, method and software based on cross section and plane linkage. The method comprises the following steps: receiving original measurement data, preliminary inversion data and target layer bottom elevation graphic data; extracting and generating the elevation file of the ground surface and the target layer bottom; converting the apparent resistivity-burial depth data into elevation domain apparent resistivity data based on an elevation reference surface through elevation domain conversion; then automatically generating the apparent resistivity cross section graph of each survey line in batches; extracting the in-seam apparent resistivity data corresponding to the target layer bottom elevation and generating the in-seam apparent resistivity plan in batches; linking the cross section graph and the plan; generating the single-point decay curve graph and the multi-channel graph; and finally automatically compiling all the graphs into a result atlas. The present application realizes the full-process automation from the data pre-processing to the generation of the result atlas, greatly improves the work efficiency and the graph accuracy, and solves the problems of complicated manual operation and easy errors.

Description

Technical Field

[0001] This invention belongs to the field of electrical exploration and relates to an electrical data mapping system, method and software based on cross-section and plane linkage. Background Technology

[0002] Mining voids and water accumulation in mining voids are major hidden geological factors causing disasters in coal mines. Once an underground water inrush accident occurs, it will cause extremely heavy economic losses to coal enterprises. Electrical exploration, as a geophysical method, is an effective means to find out the extent of mining voids and water accumulation in mining voids, and it is also an exploration method required by relevant hidden disaster regulations.

[0003] Currently, data obtained through electrical resistivity tomography is processed to generate maps, such as cross-sectional maps, plan maps, and multi-track maps, which provide key data support for subsequent determination of whether there is goaf or goaf water accumulation in the coal mining area.

[0004] Currently, there are three main processing methods for mapping electromagnetic data from electrical exploration:

[0005] First, electromagnetic data processing software includes a mapping module. Inputting data in a fixed format allows for the display of cross-sectional views, plan views, and multi-track maps. This type of electromagnetic data processing software has a built-in Surfer program. Inputting fixed-format data, or after data processing, it can generate a single transient electromagnetic cross-sectional view, plan view, or multi-track map. The output format is an image, which cannot be further modified. Modification refers to adding legends, survey line directions, target layer contour lines, and uniformly adjusting color codes to achieve a professional electrical resistivity tomography map that meets the required standards.

[0006] Secondly, professional mapping software such as Surfer and Grapher are used for map production. The data required for mapping needs to be manually edited or modified into the required format using software such as Python. Inputting data generates single transient electromagnetic cross-sectional views, plan views, or multi-track maps in SRF or GRP format. These SRF or GRP formats can be modified and have additional elements added, such as coal seam floor information, legends, survey line directions, and color scale adjustments, within the corresponding software. For batch generation of cross-sectional views, plan views, and multi-track maps, the built-in scripter plugin in Surfer is used for code editing. Technicians need to modify the code each time based on different content.

[0007] Third, software such as MATLAB and Python are used to call Surfer software to generate maps in batches using code. This requires writing code in the corresponding languages ​​of MATLAB and Python, and the functionality is limited. If batch generation of cross-sectional views is needed, new code needs to be written; if batch generation of planar views is required, new code needs to be written. Furthermore, the data needed for map generation needs to be modified to the required format using software such as Python or manually.

[0008] Overall, at present, traditional electrical resistivity tomography (EPT) data processing software (such as Surfer, Grapher, Voxler, etc.) requires manual switching, manual parameter setting, and repetitive operations, resulting in low efficiency; the preliminary data preparation work is inefficient and has a high error rate; cross-sectional views and layer-by-layer plan views cannot be linked, and after the cross-sectional view is modified, the plan view results need to be obtained again and the plan view needs to be created; the format and style of the results maps are inconsistent, and it takes a long time to organize and compile them into a map book. Summary of the Invention

[0009] To overcome the shortcomings of the aforementioned related technologies, this invention proposes an electrical resistivity data mapping system based on cross-sectional and planar linkage. This system can process electromagnetic data in batches, perform elevation correction, and simultaneously generate cross-sectional views, planar views, attenuation curves, and multi-track maps in batches. It integrates data processing and mapping processes in an automated, batch-processed, and unified manner.

[0010] The cross-sectional view can be linked with the plan view to quickly generate the plan view and improve the efficiency of graphic analysis; and an image conversion function can be added to it to batch produce standardized atlases.

[0011] To achieve the above technical objectives, the present invention provides an electrical resistivity data mapping system based on cross-sectional and planar linkage, which is applicable to data obtained by electrical resistivity exploration.

[0012] The described electrical resistivity data processing and mapping system is applicable to data obtained using electrical resistivity exploration. It includes a memory, a processor, and a computer program stored in the memory. The computer program has built-in geographic data visualization and analysis software, and is executed by the processor to perform the following steps:

[0013] S1: Data receiving steps: Receive measurement data, electromagnetic data, and apparent resistivity-depth data obtained from electrical exploration; Receive graphic data files with elevation values ​​of the target stratum bottom plate in the area to be surveyed.

[0014] S2: Preliminary data processing steps: including extracting the target layer floor and surface elevation data; generating multiple elevation files for the target layer floor based on the graphic data file of the target layer floor with elevation values ​​and the measurement data; and generating surface elevation files for multiple survey lines based on the measurement data.

[0015] S3: Elevation domain conversion step: Based on the apparent resistivity-burial depth data and measurement data, convert the burial depth data corresponding to the apparent resistivity value into an elevation domain apparent resistivity value based on the elevation datum.

[0016] S4: Batch generation steps of cross-section diagrams: Based on the apparent resistivity values ​​of the elevation domain, generate GRD grid files for multiple survey lines; repeatedly call the geographic data visualization software and read the GRD grid files in sequence to automatically generate apparent resistivity cross-section diagrams for each survey line.

[0017] S5: Stratified Information Extraction Step: Based on the target layer bottom plate graphic data file, extract the apparent resistivity values ​​of multiple measuring points corresponding to the elevation of the target layer bottom plate from the elevation domain apparent resistivity data to generate the target layer in-stratified apparent resistivity data.

[0018] S6: Batch generation steps of plan view: Generate a grd grid file based on the in-layer apparent resistivity data, and cyclically call the geographic data visualization software to generate an in-layer apparent resistivity plan view of the target layer bottom plate.

[0019] S7: Attenuation curve and multi-track generation steps: Based on the data collected by electrical resistivity tomography, generate single-point attenuation curves and multi-track maps, and output them as images;

[0020] S8: Atlas compilation steps: Convert the generated cross-sectional views, plan views, attenuation curves and multi-track maps into image formats, import and format them in batches, and automatically generate the result atlas.

[0021] Preferably, the computer program includes built-in geographic data visualization and analysis software, including Surfer software, and the processor achieves automated mapping by calling Surfer's API interface and Scripter plugin.

[0022] Preferably, generating multiple target layer bottom plate elevation files based on the target layer bottom plate graphic data file with target layer elevation values ​​and measurement data includes:

[0023] Convert the graphic data file of the target layer's base plate with elevation to a grd file and set the outer expansion range;

[0024] Identify the line number of the survey line, the point number of the measuring point on the survey line, the X-axis and Y-axis data corresponding to the measuring point, and the elevation value of each measuring point in the measurement data;

[0025] Extract the location points corresponding to the measuring points in the grd file, and obtain the base plate value of the target layer corresponding to the location points. Form a file including line number, point number, coordinate data of the location points corresponding to the measuring points in the grd file and base plate value. Then, collect the coordinate data and coordinate values ​​of the location points corresponding to multiple measuring points on each measuring line in the grd file into a target layer base plate elevation file.

[0026] When other target layers do not have base plate graphic data files, the base plate of the target layer can be pushed up or moved down based on the existing base plate of the target layer to obtain the base plate elevation files of the target layers in other parallel layers.

[0027] The process of generating surface elevation files for multiple survey lines based on measurement data includes: identifying the survey line number, point number, x-coordinate, y-coordinate, and elevation value in the measurement data; and generating the surface elevation file and whitening file required for subsequent cross-section mapping based on the blind zone range and the highest elevation value.

[0028] Preferably, the steps implemented by the computer program executed by the processor further include elevation correction, the elevation correction comprising:

[0029] Based on the apparent resistivity data and burial depth data generated after the preliminary inversion, and combined with the apparent resistivity data from multiple measuring points, elevation correction and apparent resistivity correction are performed. These include topographic correction method, fixed value correction method, and adaptive correction method.

[0030] Preferably, the method of cyclically running the geographic data visualization and analysis software and sequentially calling the cross-section mapping script to generate apparent resistivity cross-section maps of the survey lines includes:

[0031] Based on the obtained elevation file of the target layer's base plate, the corresponding surface elevation file, and the apparent resistivity value of the elevation domain, the Scripter plugin in the Surfer software is invoked to generate a map using built-in functions and scripting language. Specific steps include:

[0032] First, the processor retrieves the apparent resistivity file of the survey line after elevation correction and identifies the survey line number.

[0033] Second, the processor calls the Scripter plugin to input the elevation-corrected apparent resistivity file into the Surfer software, and sets the filtering parameters and GRD grid parameters to form the GRD grid file of the survey line.

[0034] Second, call Surfer's base map function to add multiple target layer bottom plate elevation files corresponding to the survey line to the grd grid file of a survey line, and set the line type characteristics;

[0035] Third, call Surfer's base map function, then add the corresponding surface elevation file to the grd grid file of the survey line, and set the corresponding line type features;

[0036] Fourth, call Surfer's base map function, add a whitening file corresponding to the survey line to the corresponding survey line's grd mesh file, fill it with white, and set the line type to invisible;

[0037] Fifth, arrange the base layers in order, with the whitened file being the second to last and the corresponding survey line's GRD mesh file being the last.

[0038] Sixth, set uniform color codes, scales, font sizes, survey line names, and uniform scales to form the apparent resistivity cross-sectional diagrams of the corresponding survey lines;

[0039] Seventh, the processor acquires the apparent resistivity cross-sectional diagram and GRD grid file of the survey line generated by Surfer software and stores them in the memory;

[0040] Eighth, after the processor acquires the apparent resistivity cross-sectional diagram and grd mesh file of each survey line, it will repeat steps one to seven. In step one, the processor sequentially calls the apparent resistivity files of multiple survey lines after elevation correction to form the apparent resistivity cross-sectional diagrams of multiple survey lines.

[0041] Preferably, the method for extracting bedding information and cyclically running the geographic data visualization and analysis software to generate a bedding apparent resistivity plane map includes:

[0042] First, the processor calls the grd mesh file of the multiple survey lines generated during the cross-section mapping process, as well as the elevation file and measurement file of the target layer bottom plate corresponding to the multiple survey lines;

[0043] Second, the processor identifies the survey line number corresponding to the grd grid file, and calls the elevation file and measurement data of the target layer bottom plate of the same survey line according to the identified survey line number;

[0044] Third, the processor extracts the apparent resistivity values ​​of multiple survey lines on the target layer bottom plate from the grd mesh file of multiple survey lines, and then matches the extracted apparent resistivity values ​​of multiple survey lines on the target layer bottom plate with the measurement data to form a file with line number, point number, target layer bottom plate elevation, and apparent resistivity value, thus forming the original data of the apparent resistivity value of the target layer in tandem.

[0045] Fourth, the processor calls the Scripter plugin of the Surfer software, and through its built-in functions and scripting language, it forms the GRD mesh file of the target layer from the raw data of the apparent resistivity values ​​of the target layer, and sets the outer expansion range, mesh spacing, and range of discarded values.

[0046] Fifth, the processor calls the substrate map function in the Surfer software to add the substrate map and color mark of the target layer base plate to the GRD mesh file of the target layer in sequence, and then sets the scale and font size to form a planar map of the apparent resistivity value of the target layer in sequence.

[0047] Sixth, the processor acquires the planar resistivity value of the target layer generated by Surfer software and stores it in the memory;

[0048] Seventh, after the processor obtains a plan view of the apparent resistivity value of the target layer in line, it will repeat steps one to five. In step one, the processor sequentially calls multiple target layer base plates corresponding to multiple survey lines to form multiple plan views of the apparent resistivity value of the target layer in line.

[0049] Preferably, the method for linking the cross-section and the plane includes:

[0050] Select the modified GRD mesh file of the multiple survey lines generated during the above cross-section mapping process. After selecting the other parameters according to the previous planar mapping parameters, the target layer in-line planar map after cross-section modification can be regenerated.

[0051] Preferably, the method for generating single-point attenuation curves and multi-track maps based on data collected by electrical resistivity tomography includes:

[0052] Input the original acquisition file and identify the three columns of data: measurement point, time, and voltage value.

[0053] Generate a single-point attenuation curve, with time on the horizontal axis and voltage value on the vertical axis. Each measurement point forms an attenuation curve. By setting the selected time channel, the number of measurement points, and the smoothing parameters, attenuation curves of all measurement points on each measurement line are generated.

[0054] Generate a multi-track map and establish a plane coordinate system for the multi-track map. The X-axis of the plane coordinate system represents the measurement points, and the Y-axis represents the voltage values. Connect the voltage values ​​of each measurement point at the same time to form a line, thus forming a multi-track map.

[0055] Preferably, the computer program is further executed by the processor in the following steps:

[0056] First, based on step S2, according to the apparent resistivity-burial depth data and measurement data, the burial depth data corresponding to the apparent resistivity value is converted into the apparent resistivity value of the elevation domain corresponding to the bottom plate of multiple target layers.

[0057] Second, based on the graphic data files of multiple target layer base plates and the apparent resistivity values ​​of the elevation domain corresponding to multiple target layer base plates, generate the grd grid file of the corresponding survey line on each target layer base plate.

[0058] Third, the geographic data visualization software is called in a loop, and the corresponding survey line's grd grid file is read in sequence to automatically generate the in-stratum apparent resistivity planar map of the bottom plate of each target layer.

[0059] Fourth, based on the corresponding grd grid files of the survey lines on the base plates of multiple target layers, extract the data of the survey lines from the corresponding grd grid files on each target layer base plate, and generate the grd grid files of the survey lines.

[0060] Fifth, when the spacing between multiple target layer base plates is greater than or equal to the required spacing, the apparent resistivity values ​​of the elevation domain at corresponding positions between the multiple target layer base plates are obtained by interpolation and collected into the grd grid file of the survey line.

[0061] Alternatively, when the spacing between multiple target layer base plates is greater than or equal to the distance between them, based on the apparent resistivity-burial depth data and the measurement data, the burial depth data corresponding to the apparent resistivity value is converted into the elevation domain apparent resistivity value of the corresponding layer between the multiple target layer base plates, and then collected into the grd grid file of the survey line.

[0062] Sixth, call the geographic data visualization software and read the grd grid file of the survey line and the apparent resistivity cross-section of the survey line in sequence.

[0063] On the other hand, the present invention also provides electrical resistivity data processing software based on cross-section and plane linkage, which is applicable to the above-mentioned electrical resistivity data mapping system based on cross-section and plane linkage. The electrical resistivity data processing software based on cross-section and plane linkage includes: a base plate extraction module, an elevation file creation module, an elevation correction module, a cross-section mapping module, a plane mapping module, an attenuation curve and multi-track module, and a processing module.

[0064] The base plate extraction module is configured to extract the base plate elevation of each survey line's target layer and generate the bln file required for subsequent mapping. The elevation file creation module is configured to generate bln format topographic lines and whitened bln files for each survey line. The correction module is configured to convert the original depth domain to the elevation domain and adjust the apparent resistivity values. The cross-section mapping module incorporates Surfer software and uses the Scripter plugin to call built-in functions to batch-process apparent resistivity cross-section mapping.

[0065] The apparent resistivity values ​​at the bottom plate height of the target layer where each survey line is located are extracted and merged to form a layer-in-layer apparent resistivity value file with line number, point number, survey point coordinate data, and apparent resistivity value. Through the Scripter plugin, built-in functions are called to batch generate layer-in-layer apparent resistivity planar maps.

[0066] The multi-channel module is configured to draw attenuation curves and multi-channel diagrams based on the original acquired data.

[0067] The processing module is configured to batch convert apparent resistivity cross-sectional diagrams and in-layer apparent resistivity diagrams into image format and save them as a booklet.

[0068] Furthermore, the present invention also provides an electrical resistivity data processing method based on cross-sectional and planar linkage, wherein the electrical resistivity data processing method based on cross-sectional and planar linkage performs the steps implemented in the above-mentioned electrical resistivity data mapping system based on cross-sectional and planar linkage.

[0069] The beneficial effects of this invention are as follows:

[0070] This invention processes raw data from electrical resistivity tomography (EPT) to batch generate base elevation files, surface elevation files, whitening files, etc., required for cross-section mapping. It can also batch generate cross-section maps, plan maps, attenuation curve maps, and multi-track maps through built-in software. This automates, batches, and integrates the scattered, tedious, and manual data processing and mapping process, reducing operational difficulty and error rate.

[0071] The linkage between cross-sectional and plan views improves the efficiency of map generation and has the advantage of improving the efficiency of graphic analysis.

[0072] Furthermore, the generation of various maps in this invention can ultimately lead to the automatic generation of a standardized atlas, greatly improving the efficiency and standardization of deliverables. Attached Figure Description

[0073] To more clearly illustrate the technical solutions in the embodiments of the present invention or related technologies, the drawings used in the description of the embodiments or related technologies will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0074] Figure 1 This is a system flowchart of the present invention;

[0075] Figure 2 This is a cross-sectional view of the apparent resistivity of the present invention;

[0076] Figure 3 This is a plan view of the apparent resistivity along the layers of the present invention;

[0077] Figure 4 This is the multi-track diagram of the present invention;

[0078] Figure 5 This is a single-point attenuation curve diagram of the present invention. Detailed Implementation

[0079] To make the above-mentioned objects, features, and advantages of the present invention more apparent and understandable, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0080] In the description of this invention, it should be understood that the terms "center", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", etc., indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings. They are only for the convenience of describing this invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation. Therefore, they should not be construed as limitations on this invention.

[0081] The terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of this invention, unless otherwise stated, "a plurality of" means two or more.

[0082] Some embodiments of the present invention provide one aspect: the present invention proposes an electrical resistivity data mapping system based on cross-section and planar linkage, which can batch process the preliminary data of electromagnetic methods, perform elevation correction, and simultaneously batch draw cross-section maps, planar maps, attenuation curve maps and multi-track maps, thus automating, batching and integrating the data processing and mapping process.

[0083] The cross-sectional view can be linked with the plan view to quickly generate the plan view and improve the efficiency of graphic analysis; and an image conversion function can be added to it to batch produce standardized atlases.

[0084] To achieve the above technical objectives, the present invention provides an electrical resistivity data mapping system based on cross-sectional and planar linkage, which is applicable to data obtained by electrical resistivity exploration.

[0085] The described electrical resistivity data processing and mapping system is applicable to data obtained using electrical resistivity exploration. It includes a memory, a processor, and a computer program stored in the memory. The computer program has built-in geographic data visualization and analysis software, and is executed by the processor to perform the following steps:

[0086] S1: Data receiving steps: Receive measurement data, electromagnetic data, and apparent resistivity-depth data obtained from electrical exploration; Receive graphic data files with elevation values ​​of the target stratum bottom plate in the area to be surveyed.

[0087] S2: Preliminary data processing steps: including extracting the target layer floor and surface elevation data; generating multiple elevation files for the target layer floor based on the graphic data file of the target layer floor with elevation values ​​and the measurement data; and generating surface elevation files for multiple survey lines based on the measurement data.

[0088] S3: Elevation domain conversion step: Based on the apparent resistivity-burial depth data and measurement data, convert the burial depth data corresponding to the apparent resistivity value into an elevation domain apparent resistivity value based on the elevation datum.

[0089] S4: Batch generation steps of cross-section diagrams: Based on the apparent resistivity values ​​of the elevation domain, generate GRD grid files for multiple survey lines; repeatedly call the geographic data visualization software and read the GRD grid files in sequence to automatically generate apparent resistivity cross-section diagrams for each survey line.

[0090] S5: Stratified Information Extraction Step: Based on the target layer bottom plate graphic data file, extract the apparent resistivity values ​​of multiple measuring points corresponding to the elevation of the target layer bottom plate from the elevation domain apparent resistivity data to generate the target layer in-stratified apparent resistivity data.

[0091] S6: Batch generation steps of plan view: Generate a grd grid file based on the in-layer apparent resistivity data, and cyclically call the geographic data visualization software to generate an in-layer apparent resistivity plan view of the target layer bottom plate.

[0092] S7: Attenuation curve and multi-track generation steps: Based on the data collected by electrical resistivity tomography, generate single-point attenuation curves and multi-track maps, and output them as images;

[0093] S8: Atlas compilation steps: Convert the generated cross-sectional views, plan views, attenuation curves and multi-track maps into image formats, import and format them in batches, and automatically generate the result atlas.

[0094] Preferably, the computer program includes built-in geographic data visualization and analysis software, including Surfer software, and the processor achieves automated mapping by calling Surfer's API interface and Scripter plugin.

[0095] Preferably, generating multiple target layer bottom plate elevation files based on the target layer bottom plate graphic data file with target layer elevation values ​​and measurement data includes:

[0096] Convert the graphic data file of the target layer's base plate with elevation to a grd file and set the outer expansion range;

[0097] Identify the line number of the survey line, the point number of the measuring point on the survey line, the X-axis and Y-axis data corresponding to the measuring point, and the elevation value of each measuring point in the measurement data;

[0098] Extract the location points corresponding to the measuring points in the grd file, and obtain the base plate value of the target layer corresponding to the location points. Form a file including line number, point number, coordinate data of the location points corresponding to the measuring points in the grd file and base plate value. Then, collect the coordinate data and coordinate values ​​of the location points corresponding to multiple measuring points on each measuring line in the grd file into a target layer base plate elevation file.

[0099] When other target layers do not have base plate graphic data files, the base plate of the target layer can be pushed up or moved down based on the existing base plate of the target layer to obtain the base plate elevation files of the target layers in other parallel layers.

[0100] The process of generating surface elevation files for multiple survey lines based on measurement data includes: identifying the survey line number, point number, x-coordinate, y-coordinate, and elevation value in the measurement data; and generating the surface elevation file and whitening file required for subsequent cross-section mapping based on the blind zone range and the highest elevation value.

[0101] Preferably, the steps implemented by the computer program executed by the processor further include elevation correction, the elevation correction comprising:

[0102] Based on the apparent resistivity data and burial depth data generated after the preliminary inversion, and combined with the apparent resistivity data from multiple measuring points, elevation correction and apparent resistivity correction are performed. These include topographic correction method, fixed value correction method, and adaptive correction method.

[0103] Preferably, the method of cyclically running the geographic data visualization and analysis software and sequentially calling the cross-section mapping script to generate apparent resistivity cross-section maps of the survey lines includes:

[0104] Based on the obtained elevation file of the target layer's base plate, the corresponding surface elevation file, and the apparent resistivity value of the elevation domain, the Scripter plugin in the Surfer software is invoked to generate a map using built-in functions and scripting language. Specific steps include:

[0105] First, the processor retrieves the apparent resistivity file of the survey line after elevation correction and identifies the survey line number.

[0106] Second, the processor calls the Scripter plugin to input the elevation-corrected apparent resistivity file into the Surfer software, and sets the filtering parameters and GRD grid parameters to form the GRD grid file of the survey line.

[0107] Second, call Surfer's base map function to add multiple target layer bottom plate elevation files corresponding to the survey line to the grd grid file of a survey line, and set the line type characteristics;

[0108] Third, call Surfer's base map function, then add the corresponding surface elevation file to the grd grid file of the survey line, and set the corresponding line type features;

[0109] Fourth, call Surfer's base map function, add a whitening file corresponding to the survey line to the corresponding survey line's grd mesh file, fill it with white, and set the line type to invisible;

[0110] Fifth, arrange the base layers in order, with the whitened file being the second to last and the corresponding survey line's GRD mesh file being the last.

[0111] Sixth, set uniform color codes, scales, font sizes, survey line names, and uniform scales to form the apparent resistivity cross-sectional diagrams of the corresponding survey lines;

[0112] Seventh, the processor acquires the apparent resistivity cross-sectional diagram and GRD grid file of the survey line generated by Surfer software and stores them in the memory;

[0113] Eighth, after the processor acquires the apparent resistivity cross-sectional diagram and grd mesh file of each survey line, it will repeat steps one to seven. In step one, the processor sequentially calls the apparent resistivity files of multiple survey lines after elevation correction to form the apparent resistivity cross-sectional diagrams of multiple survey lines.

[0114] Preferably, the method for extracting bedding information and cyclically running the geographic data visualization and analysis software to generate a bedding apparent resistivity plane map includes:

[0115] First, the processor calls the grd mesh file of the multiple survey lines generated during the cross-section mapping process, as well as the elevation file and measurement file of the target layer bottom plate corresponding to the multiple survey lines;

[0116] Second, the processor identifies the survey line number corresponding to the grd grid file, and calls the elevation file and measurement data of the target layer bottom plate of the same survey line according to the identified survey line number;

[0117] Third, the processor extracts the apparent resistivity values ​​of multiple survey lines on the target layer bottom plate from the grd mesh file of multiple survey lines, and then matches the extracted apparent resistivity values ​​of multiple survey lines on the target layer bottom plate with the measurement data to form a file with line number, point number, target layer bottom plate elevation, and apparent resistivity value, thus forming the original data of the apparent resistivity value of the target layer in tandem.

[0118] Fourth, the processor calls the Scripter plugin of the Surfer software, and through its built-in functions and scripting language, it forms the GRD mesh file of the target layer from the raw data of the apparent resistivity values ​​of the target layer, and sets the outer expansion range, mesh spacing, and range of discarded values.

[0119] Fifth, the processor calls the substrate map function in the Surfer software to add the substrate map and color mark of the target layer base plate to the GRD mesh file of the target layer in sequence, and then sets the scale and font size to form a planar map of the apparent resistivity value of the target layer in sequence.

[0120] Sixth, the processor acquires the planar resistivity value of the target layer generated by Surfer software and stores it in the memory;

[0121] Seventh, after the processor obtains a plan view of the apparent resistivity value of the target layer in line, it will repeat steps one to five. In step one, the processor sequentially calls multiple target layer base plates corresponding to multiple survey lines to form multiple plan views of the apparent resistivity value of the target layer in line.

[0122] Preferably, the method for linking the cross-section and the plane includes:

[0123] Select the modified GRD mesh file of the multiple survey lines generated during the above cross-section mapping process. After selecting the other parameters according to the previous planar mapping parameters, the target layer in-line planar map after cross-section modification can be regenerated.

[0124] Preferably, the method for generating single-point attenuation curves and multi-track maps based on data collected by electrical resistivity tomography includes:

[0125] Input the original acquisition file and identify the three columns of data: measurement point, time, and voltage value.

[0126] Generate a single-point attenuation curve, with time on the horizontal axis and voltage value on the vertical axis. Each measurement point forms an attenuation curve. By setting the selected time channel, the number of measurement points, and the smoothing parameters, attenuation curves of all measurement points on each measurement line are generated.

[0127] Generate a multi-track map and establish a plane coordinate system for the multi-track map. The X-axis of the plane coordinate system represents the measurement points, and the Y-axis represents the voltage values. Connect the voltage values ​​of each measurement point at the same time to form a line, thus forming a multi-track map.

[0128] Preferably, the computer program is further executed by the processor in the following steps:

[0129] First, based on step S2, according to the apparent resistivity-burial depth data and measurement data, the burial depth data corresponding to the apparent resistivity value is converted into the apparent resistivity value of the elevation domain corresponding to the bottom plate of multiple target layers.

[0130] Second, based on the graphic data files of multiple target layer base plates and the apparent resistivity values ​​of the elevation domain corresponding to multiple target layer base plates, generate the grd grid file of the corresponding survey line on each target layer base plate.

[0131] Third, the geographic data visualization software is called in a loop, and the corresponding survey line's grd grid file is read in sequence to automatically generate the in-stratum apparent resistivity planar map of the bottom plate of each target layer.

[0132] Fourth, based on the corresponding grd grid files of the survey lines on the base plates of multiple target layers, extract the data of the survey lines from the corresponding grd grid files on each target layer base plate, and generate the grd grid files of the survey lines.

[0133] Fifth, when the spacing between multiple target layer base plates is greater than or equal to the required spacing, the apparent resistivity values ​​of the elevation domain at corresponding positions between the multiple target layer base plates are obtained by interpolation and collected into the grd grid file of the survey line.

[0134] Alternatively, when the spacing between multiple target layer base plates is greater than or equal to the distance between them, based on the apparent resistivity-burial depth data and the measurement data, the burial depth data corresponding to the apparent resistivity value is converted into the elevation domain apparent resistivity value of the corresponding layer between the multiple target layer base plates, and then collected into the grd grid file of the survey line.

[0135] Sixth, call the geographic data visualization software and read the grd grid file of the survey line and the apparent resistivity cross-section of the survey line in sequence.

[0136] On the other hand, the present invention also provides electrical resistivity data processing software based on cross-section and plane linkage, which is applicable to the above-mentioned electrical resistivity data mapping system based on cross-section and plane linkage. The electrical resistivity data processing software based on cross-section and plane linkage includes: a base plate extraction module, an elevation file creation module, an elevation correction module, a cross-section mapping module, a plane mapping module, an attenuation curve and multi-track module, and a processing module.

[0137] The base plate extraction module is configured to extract the base plate elevation of each survey line's target layer and generate the bln file required for subsequent mapping. The elevation file creation module is configured to generate bln format topographic lines and whitened bln files for each survey line. The correction module is configured to convert the original depth domain to the elevation domain and adjust the apparent resistivity values. The cross-section mapping module incorporates Surfer software and uses the Scripter plugin to call built-in functions to batch-process apparent resistivity cross-section mapping.

[0138] The apparent resistivity values ​​at the bottom plate height of the target layer where each survey line is located are extracted and merged to form a layer-in-layer apparent resistivity value file with line number, point number, survey point coordinate data, and apparent resistivity value. Through the Scripter plugin, built-in functions are called to batch generate layer-in-layer apparent resistivity planar maps.

[0139] The multi-channel module is configured to draw attenuation curves and multi-channel diagrams based on the original acquired data.

[0140] The processing module is configured to batch convert apparent resistivity cross-sectional diagrams and in-layer apparent resistivity diagrams into image format and save them as a booklet.

[0141] Furthermore, the present invention also provides an electrical resistivity data processing method based on cross-sectional and planar linkage, wherein the electrical resistivity data processing method based on cross-sectional and planar linkage performs the steps implemented in the above-mentioned electrical resistivity data mapping system based on cross-sectional and planar linkage.

[0142] In the description of this specification, specific features, structures, materials, or characteristics may be combined in any suitable manner in one or more embodiments or examples.

[0143] The above are merely specific embodiments of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A mapping system for electrical resistivity data based on cross-section and planar linkage, applicable to data obtained through electrical resistivity exploration, characterized in that, It includes a memory, a processor, and a computer program stored on the memory, the computer program having built-in geographic data visualization and analysis software, the computer program being executed by the processor to perform the following steps: S1: Data receiving steps: Receive measurement data, electromagnetic data and apparent resistivity-depth data obtained from electrical exploration. Receive graphic data files with elevation values ​​of the target layer bottom plate of the area to be surveyed; S2: Preliminary data processing steps: including extracting the target layer floor and surface elevation data; generating multiple elevation files for the target layer floor based on the graphic data file of the target layer floor with elevation values ​​and the measurement data; and generating surface elevation files for multiple survey lines based on the measurement data. S3: Elevation domain conversion step: Based on the apparent resistivity-burial depth data and measurement data, convert the burial depth data corresponding to the apparent resistivity value into an elevation domain apparent resistivity value based on the elevation datum. S4: Batch generation steps of cross-section diagrams: Based on the apparent resistivity values ​​of the elevation domain, generate GRD grid files for multiple survey lines; repeatedly call the geographic data visualization software and read the GRD grid files in sequence to automatically generate apparent resistivity cross-section diagrams for each survey line. S5: Stratified Information Extraction Step: Extract the apparent resistivity values ​​of multiple measuring points corresponding to the elevation of the target layer's bottom plate from the apparent resistivity data in the elevation domain, and generate the apparent resistivity data of the target layer in stratified resistivity. S6: Batch generation steps of plan view: Generate a grd grid file based on the in-layer apparent resistivity data, and cyclically call the geographic data visualization software to generate an in-layer apparent resistivity plan view of the target layer bottom plate. S7: Attenuation curve and multi-track generation steps: Based on the data collected by electrical resistivity tomography, generate single-point attenuation curves and multi-track maps, and output them as images; S8: Atlas compilation steps: Convert the generated cross-sectional views, plan views, attenuation curves and multi-track maps into image formats, import and format them in batches, and automatically generate the result atlas.

2. The electrical resistivity tomography system based on cross-sectional and planar linkage according to claim 1, characterized in that, The computer program has built-in geographic data visualization and analysis software, including Surfer software. The processor achieves automated mapping by calling Surfer's API interface and Scripter plugin.

3. The electrical resistivity tomography system based on cross-sectional and planar linkage according to claim 1 or 2, characterized in that, The process of generating multiple target layer bottom slab elevation files based on the target layer bottom slab graphic data file with elevation values ​​and measurement data includes: Convert the graphic data file of the target layer's base plate with elevation to a grd file and set the outer expansion range; Identify the line number of the survey line, the point number of the measuring point on the survey line, the X-axis and Y-axis data corresponding to the measuring point, and the elevation value of each measuring point in the measurement data; Extract the location points corresponding to the measuring points in the grd file, and obtain the base plate value of the target layer corresponding to the location points. Form a file including line number, point number, coordinate data of the location points corresponding to the measuring points in the grd file and base plate value. Then, collect the coordinate data and coordinate values ​​of the location points corresponding to multiple measuring points on each measuring line in the grd file into a target layer base plate elevation file.

4. The electrical resistivity tomography system based on cross-sectional and planar linkage according to claim 3, characterized in that, The steps implemented by the computer program executed by the processor further include elevation correction, the elevation correction comprising: Based on the apparent resistivity data and burial depth data generated after the preliminary inversion, and combined with the apparent resistivity data from multiple measuring points, elevation correction and apparent resistivity correction are performed.

5. The electrical resistivity tomography system based on cross-sectional and planar linkage according to claim 2, characterized in that, The method of cyclically running the geographic data visualization and analysis software and sequentially calling the cross-section mapping script to generate apparent resistivity cross-section maps of the survey lines includes: Based on the obtained elevation file of the target layer's base plate, the corresponding surface elevation file, and the apparent resistivity value of the elevation domain, the Scripter plugin in the Surfer software is invoked to generate a map using built-in functions and scripting language. Specific steps include: First, the processor retrieves the apparent resistivity file of the survey line after elevation correction and identifies the survey line number. Second, the processor calls the Scripter plugin to input the elevation-corrected apparent resistivity file into the Surfer software, and sets the filtering parameters and GRD grid parameters to form the GRD grid file of the survey line. Second, call Surfer's base map function to add multiple target layer bottom plate elevation files corresponding to the survey line to the grd grid file of a survey line, and set the line type characteristics; Third, call Surfer's base map function, then add the corresponding surface elevation file to the grd grid file of the survey line, and set the corresponding line type features; Fourth, call Surfer's base map function, add a whitening file corresponding to the survey line to the corresponding survey line's grd mesh file, fill it with white, and set the line type to invisible; Fifth, arrange the base layers in order, with the whitened file being the second to last and the corresponding survey line's GRD mesh file being the last. Sixth, set uniform color codes, scales, font sizes, survey line names, and uniform scales to form the apparent resistivity cross-sectional diagrams of the corresponding survey lines; Seventh, the processor acquires the apparent resistivity cross-sectional diagram and GRD grid file of the survey line generated by Surfer software and stores them in the memory; Eighth, after the processor acquires the apparent resistivity cross-sectional diagram and grd mesh file of each survey line, it will repeat steps one to seven. In step one, the processor sequentially calls the apparent resistivity files of multiple survey lines after elevation correction to form the apparent resistivity cross-sectional diagrams of multiple survey lines.

6. The electrical resistivity tomography system based on cross-sectional and planar linkage according to claim 5, characterized in that, The method for extracting bedding information and cyclically running the geographic data visualization and analysis software to generate a bedding apparent resistivity plane map includes: First, the processor calls the grd mesh files of multiple survey lines, the elevation files and measurement files of the target layer bottom plate corresponding to the multiple survey lines; Second, the processor identifies the survey line number corresponding to the grd grid file, and calls the elevation file and measurement data of the target layer bottom plate of the same survey line according to the identified survey line number; Third, the processor extracts the apparent resistivity values ​​of multiple survey lines on the target layer bottom plate from the grd mesh file of multiple survey lines, and then matches the extracted apparent resistivity values ​​of multiple survey lines on the target layer bottom plate with the measurement data to form a file with line number, point number, target layer bottom plate elevation, and apparent resistivity value, thus forming the original data of the apparent resistivity value of the target layer in tandem. Fourth, the processor calls the Scripter plugin of the Surfer software, and through its built-in functions and scripting language, it forms the GRD mesh file of the target layer from the raw data of the apparent resistivity values ​​of the target layer, and sets the outer expansion range, mesh spacing, and range of discarded values. Fifth, the processor calls the substrate map function in the Surfer software to add the substrate map and color mark of the target layer base plate to the GRD mesh file of the target layer in sequence, and then sets the scale and font size to form a planar map of the apparent resistivity value of the target layer in sequence. Sixth, the processor acquires the planar resistivity value of the target layer generated by Surfer software and stores it in the memory; Seventh, after the processor obtains a plan view of the apparent resistivity value of the target layer in line, it will repeat steps one to five. In step one, the processor sequentially calls multiple target layer base plates corresponding to multiple survey lines to form multiple plan views of the apparent resistivity value of the target layer in line.

7. The electrical resistivity tomography system based on cross-sectional and planar linkage according to claim 2, characterized in that, The method for generating single-point attenuation curves and multi-track plots based on raw data includes: Input the original data acquisition file and identify the three columns of data: measurement point, time, and voltage value. Generate a single-point attenuation curve, with time on the horizontal axis and voltage value on the vertical axis. Each measurement point forms an attenuation curve. By setting the selected time channel, the number of measurement points, and the smoothing parameters, attenuation curves of all measurement points on each measurement line are generated. Generate a multi-track map and establish a plane coordinate system for the multi-track map. The X-axis of the plane coordinate system represents the measurement points, and the Y-axis represents the voltage values. Connect the voltage values ​​of each measurement point at the same time to form a line, thus forming a multi-track map.

8. The electrical resistivity tomography system based on cross-sectional and planar linkage according to claim 2, characterized in that, The computer program is further executed by the processor in the following steps: Based on step S2, according to the apparent resistivity-burial depth data and measurement data, the burial depth data corresponding to the apparent resistivity value is converted into the apparent resistivity value of the elevation domain corresponding to the bottom plate of multiple target layers. Based on the graphic data files of multiple target layer base plates and the apparent resistivity values ​​of the elevation domain corresponding to multiple target layer base plates, generate the grd grid file of the corresponding survey line on each target layer base plate for multiple survey lines; The system continuously calls the geographic data visualization software and reads the corresponding survey line's grd grid file in sequence to automatically generate the in-stratum apparent resistivity planar map of the bottom plate of each target layer. Based on the GRD mesh files of the corresponding survey lines on the base plates of multiple target layers, extract the data of the survey lines from the corresponding GRD mesh files on each target layer base plate, and generate the GRD mesh files of the survey lines. When the spacing between multiple target layer base plates is greater than or equal to the distance between them, the apparent resistivity values ​​of the elevation domain at the corresponding positions between the multiple target layer base plates are obtained by interpolation and collected into the grd grid file of the survey line. Alternatively, when the spacing between multiple target layer base plates is greater than or equal to the apparent resistivity-burial depth data and measurement data, the burial depth data corresponding to the apparent resistivity value is converted into the elevation domain apparent resistivity value of the corresponding layer between multiple target layer base plates, and collected into the grd grid file of the survey line. The geographic data visualization software is invoked, and the grd grid file of the survey line and the apparent resistivity cross-section of the survey line are read sequentially.

9. A software for processing electrical resistivity data based on cross-sectional and planar linkage, applicable to the electrical resistivity data mapping system based on cross-sectional and planar linkage as described in any one of claims 1 to 8, characterized in that, The electrical resistivity data processing software based on cross-sectional and planar linkage includes: The base plate extraction module is configured to extract the base plate elevation of each survey line target layer and generate a bln file required for subsequent mapping. The elevation file creation module is configured to generate bln format for each survey line and topographic line, as well as whitened bln format files. An elevation correction module is configured to convert the original depth domain into an elevation domain and adjust the apparent resistivity value. The cross-section mapping module has built-in Surfer software and uses the Scripter plugin to call built-in functions to batch complete the apparent resistivity cross-section mapping. The planar mapping module extracts the apparent resistivity values ​​of the target layer bottom plate height of each survey line and merges them to form a layer-by-layer apparent resistivity value file with line number, point number, survey point coordinate data, and apparent resistivity value. Through the Scripter plugin, built-in functions are called to batch complete the layer-by-layer apparent resistivity planar maps. Attenuation curve and multi-channel module, wherein the multi-channel module is configured to: draw attenuation curve and multi-channel graph based on the original acquired data; The processing module is configured to batch convert apparent resistivity cross-sectional diagrams and in-layer apparent resistivity diagrams into image format and save them as a booklet.

10. A method for processing electrical resistivity data based on cross-sectional and planar linkage, characterized in that, The electrical resistivity data processing method based on cross-section and planar linkage performs the steps implemented in the electrical resistivity data mapping system based on cross-section and planar linkage as described in any one of claims 1 to 8.

Images