A method and system for drawing geological distribution maps of power transmission and transformation projects based on multi-source data

Abstract

The present invention provides a method and system for drawing geological distribution maps of power transmission and transformation projects based on multi-source data. Based on multi-source data such as digital geological maps, terrain elevation models, and satellite images, the digital regional geological map is initialized and partitioned, the partitions are corrected according to the satellite images and terrain elevation models, and further partitioned according to the micro-topography. After considering whether to filter and discard the data based on the degree of landform damage and the digital geological map information, the coverage layers of the exploration points in the partitions are counted, and assignment processing is performed based on the partition results. The assignment calculation is performed according to the partition categories in order of the geological age of the partitions from old to new. Compared with the existing technology, the drawing area is pre-partitioned and processed based on the boundaries of the multi-source data, and the cover layer thickness parameter interpolation calculation is performed in a targeted manner, which effectively solves the problem of imperfect drawing of geological distribution maps of power transmission and transformation projects under certain accuracy when the data points are limited and unevenly distributed.

Description

Technical Field

[0001] The present invention relates to the technical field of geological distribution map drawing, and in particular to a method and system for drawing geological distribution maps of power transmission and transformation projects based on multi-source data. Background Art

[0002] The geological distribution map of power transmission and transformation projects is a digital result that reflects geological conditions for power transmission and transformation projects, using survey data of power transmission and transformation projects and adopting unified drawing standards. The results can assist in management work such as mechanized construction and selection of environmentally friendly foundations. It mainly includes information such as the thickness of the cover layer, the basic quality grade of the bedrock, and topographic and geographical information. The "Standards for Drawing Geological Distribution Maps of Power Transmission and Transformation Projects (2021 Edition)" of the State Grid Infrastructure Department has detailed provisions for the geological distribution map of power transmission and transformation projects: the cover layer thickness is divided into 7 intervals, corresponding to the cover layer 0-2m (inclusive), 2-4m (inclusive), 4-6m (inclusive), 6-8m (inclusive), 8-10m (inclusive), 8-12m (inclusive), and >12m. The "Standards for Drawing Geological Distribution Maps of Power Transmission and Transformation Projects (2021 Edition)" provides a standardized drawing method as follows: ① Data screening. 1. Eliminate duplicate point data; when the thickness of the overburden at adjacent points varies significantly, analyze the data's credibility based on the origin of the strata, topography, and regional engineering geological characteristics, and eliminate erroneous data. 2. Data reclassification. Prioritize geological data obtained from transmission and transformation projects, followed by extracting relevant information from various nationally released vector graphics. If neither of these methods is available, interpolation methods such as kriging and inverse distance are used to calculate data at the target location. 3. Map-based structured data processing. Prioritize automatic calculation and direct extraction of information such as overburden thickness and basic bedrock quality grade. Transmission and transformation project geological distribution maps can be used for selecting environmental protection foundations for transmission projects and for survey and planning of substation projects.

[0003] Due to the complexity of geological conditions, the above standards do not specify the method for drawing the cover thickness. In actual drawing, some methods directly use interpolation methods such as Kriging to calculate and draw within the drawing area; some divide the area into bedrock area and soil area based on the characteristics of geological conditions, and set the cover thickness in the bedrock area to 0-2m. The remaining areas use power transmission and transformation project data for interpolation; due to the spatial unevenness and finiteness of the distribution of power transmission and transformation projects and the complexity of geological conditions, the calculation results of the above methods are insufficient. The two methods do not take into account that a single interpolation method is only intra-regional interpolation, while the geological conditions are extremely complex in spatial distribution, with regional laws and microscopic complexity. The first method is only applicable to the simple distribution of a single stratum. The second method distinguishes between bedrock areas and soil areas, but is also only applicable to specific simple geological conditions. It is not universal and has a limited scope of application. For general sites, the error is extremely large and the actual situation is not taken into account at all. Summary of the Invention

[0004] In view of this, the purpose of the embodiments of the present invention is to provide a relatively reliable and general method that is suitable for the construction and management needs of power transmission and transformation projects. It uses existing survey data and collected data to solve the contradiction between the limited number and distribution of existing exploration point data and the complexity of geological conditions, thereby effectively solving the problem of drawing geological distribution maps of power transmission and transformation projects with a certain accuracy when the data points are limited and unevenly distributed.

[0005] A first aspect of the present invention provides a method for drawing a geological distribution map of a power transmission and transformation project based on multi-source data, the method comprising a preprocessing step, a pre-assignment step, and an assignment calculation step; wherein:

[0006] The preprocessing step includes taking a single regional geological map as a unit, extracting the rock and soil age genesis and lithology of the digital regional geological map based on multi-source data such as the digital geological map, terrain elevation model, and satellite imagery, performing initial zoning, revising the zoning based on satellite imagery and terrain elevation model, and further zoning based on micro-relief zoning;

[0007] The pre-assignment step includes projecting the collected exploration data into the partitions, considering whether to filter and discard the data based on the degree of landform damage and digital geological map information, counting the coverage layers of the exploration points in the partitions, and performing assignment processing based on the partition results;

[0008] The steps of the value assignment calculation include performing value assignment calculations according to the order of the geological age of the partition from old to new, and according to the partition category. The specific rules are: first, interpolation calculations are performed on the interpolation area with exploration points, and ordinary kriging method or triangular grid method is used for interpolation within the same partition range. According to the interpolation results, the partition is divided into single-valued partitions, and the spatial distribution, geographic information, micro-topography and other information of the partition are classified and statistically analyzed to determine whether the above interpolation results can be updated through the classification standard.

[0009] Furthermore, the further partitioning according to micro-topography includes:

[0010] Based on the corresponding regional geological survey data, consider whether to add orientation, slope, and spatial distribution factors for further zoning; the attributes of the zoning include geological age, geological causes, and micro-geomorphological factors.

[0011] Furthermore, in the pre-assignment step, the assignment process includes:

[0012] In the case of a single partition with an exploration point in the same type of partition, the value is directly assigned and all partitions of this type are marked as assignment areas;

[0013] When a single zone in the same type of zone has a large range or a large span of cover thickness, the user shall review the micro-topography and spatial distribution factors to determine whether to add classification criteria for further zone division. If further zone division is required, the preprocessing step shall be repeated. Otherwise, it shall be marked as an interpolation area. The zone without exploration points that is consistent with the attributes of the assignment zone according to the classification criteria shall be marked as the assignment zone, and the remaining areas without exploration points shall be marked as interpolation areas. The zone without exploration points that is consistent with the attributes of the assignment zone according to the classification criteria shall be marked as the assignment zone, and the remaining areas without exploration points shall be marked as interpolation areas.

[0014] Furthermore, the assignment calculation step further includes:

[0015] For partitions without exploration points, preprocessing should be performed again when the classification criteria are updated;

[0016] For the assigned areas of the same type without exploration points, the inverse distance-based interpolation method is used for the assignment calculation. The inverse square of the distance between the partition to be calculated and the surrounding similar assigned areas is used as the weight for interpolation calculation. To reduce the amount of calculation, the four partitions closest to each other within a fixed distance are generally used for calculation; wherein, the fixed distance is selected according to the accuracy of the mapping;

[0017] For the interpolation area without exploration points, professionals will perform the assignment process based on the above assignment results;

[0018] For similar zones where the distribution of cover thickness at exploration points is irregular or where geotechnical professional analysis shows that the cover thickness of such zones is difficult to estimate, the values ​​shall be comprehensively determined according to the standards that affect the application and selection of rock anchors.

[0019] Furthermore, the method also includes an inspection step: experts inspect and process the above-mentioned assignment calculation results and make adjustments to individual areas.

[0020] In addition, the second aspect of the present invention further provides a system for drawing geological distribution maps of power transmission and transformation projects based on multi-source data, the system comprising a pre-processing module, a pre-assignment module, and an assignment calculation module;

[0021] The preprocessing module is used to extract the geotechnical age and lithology of the digital regional geological map based on a single regional geological map, terrain elevation model, satellite imagery and other multi-source data, perform initial zoning, perform zoning correction based on satellite imagery and terrain elevation model, and further zoning based on micro-geomorphological zoning;

[0022] The pre-assignment module is used to project the collected exploration data into the partitions, consider whether to filter and discard the data based on the degree of landform damage and digital geological map information, calculate the coverage of the exploration points in the partitions, and perform assignment processing based on the partition results;

[0023] The assignment calculation module is used to perform assignment calculations according to the geological age of the partition from old to new and according to the partition category. The specific rules are: first perform interpolation calculations on the interpolation area with exploration points, and use ordinary kriging or triangular grid method to interpolate within the same partition range. According to the interpolation results, the partition is divided into single-assignment partitions, and the spatial distribution, geographic information, micro-topography and other information of the partition are classified and statistically analyzed to determine whether the above interpolation results can be updated through the classification standard.

[0024] Furthermore, the pre-processing module is also used to: consider whether to add orientation, slope, and spatial distribution factors for further partitioning based on the corresponding regional geological survey data; the attributes of the partition include geological age, geological origin, and micro-geomorphological factors.

[0025] Furthermore, the pre-assignment module is further configured to: directly assign a value to a single partition of the same type with an exploration point, and identify all partitions of this type as assignment areas;

[0026] When a single zone in the same type of zone has a large range or a large span of cover thickness, the user shall review the micro-topography and spatial distribution factors to determine whether to add classification criteria for further zone division. If further zone division is required, the preprocessing step shall be repeated. Otherwise, it shall be marked as an interpolation area. The zone without exploration points that is consistent with the attributes of the assignment zone according to the classification criteria shall be marked as the assignment zone, and the remaining areas without exploration points shall be marked as interpolation areas. The zone without exploration points that is consistent with the attributes of the assignment zone according to the classification criteria shall be marked as the assignment zone, and the remaining areas without exploration points shall be marked as interpolation areas.

[0027] Furthermore, the assignment calculation module is further configured to: for partitions without exploration points, re-preprocess when the classification criteria are updated;

[0028] For the assigned areas of the same type without exploration points, the inverse distance-based interpolation method is used for the assignment calculation. The inverse square of the distance between the partition to be calculated and the surrounding similar assigned areas is used as the weight for interpolation calculation. To reduce the amount of calculation, the four partitions closest to each other within a fixed distance are generally used for calculation; wherein, the fixed distance is selected according to the accuracy of the mapping;

[0029] For the interpolation area without exploration points, professionals will perform the assignment process based on the above assignment results;

[0030] For similar zones where the distribution of cover thickness at exploration points is irregular or where geotechnical professional analysis shows that the cover thickness of such zones is difficult to estimate, the values ​​shall be comprehensively determined according to the standards that affect the application and selection of rock anchors.

[0031] Furthermore, the system also includes a checking module for checking and processing the above-mentioned assignment calculation results according to experts and adjusting individual areas.

[0032] In the solution of the present invention, a single regional geological map is used as a unit. Based on multi-source data such as a digital geological map, a terrain elevation model, and satellite images, the digital regional geological map is extracted to initialize the zoning by geotechnical age and lithology. The zoning is then corrected based on the satellite image and terrain elevation model, and further zoned based on micro-geomorphological zoning. After the collected exploration data is projected onto the zones, the data is screened and discarded based on the degree of geomorphological damage and the digital geological map information. The coverage of the exploration points within the zones is statistically analyzed, and a value assignment process is performed based on the zone results. The value assignment calculation is performed sequentially according to the zone categories in descending order of geological age. The specific rules are as follows: interpolation calculation is first performed on the interpolation area with exploration points. Ordinary kriging or triangulated meshing is used to interpolate within the same zone. Based on the interpolation results, the zone is divided into single-valued zones. Classification and statistical analysis of the spatial distribution, geographic information, micro-geomorphological information, and other information of the zones is performed to determine whether the interpolation results can be updated according to the classification criteria. Compared with the existing technology, by introducing multi-source data such as digital geological maps, terrain elevation models, satellite images, etc., the drawing area is pre-partitioned and processed based on the boundaries of the multi-source data, and on this basis, the cover layer thickness parameter interpolation calculation is carried out in a targeted manner, which effectively solves the problem of drawing geological distribution maps of power transmission and transformation projects with a certain accuracy when the data points are limited and unevenly distributed. BRIEF DESCRIPTION OF THE DRAWINGS

[0033] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the following briefly introduces the drawings required for use in the embodiments. It should be understood that the following drawings only illustrate certain embodiments of the present invention and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other relevant drawings can be obtained based on these drawings without paying any creative work.

[0034] Figure 1 This is a flow chart of a method for drawing a geological distribution map of a power transmission and transformation project based on multi-source data disclosed in Example 1 of the present invention;

[0035] Figure 2 It is a structural schematic diagram of the geological distribution map drawing system for power transmission and transformation projects based on multi-source data disclosed in Example 2 of the present invention. DETAILED DESCRIPTION

[0036] Example embodiments will now be described more fully with reference to the accompanying drawings. However, example embodiments can be implemented in many forms and should not be construed as limited to the examples set forth herein; rather, these embodiments are provided so that this application will be thorough and complete and will fully convey the concepts of the example embodiments to those skilled in the art.

[0037] In addition, described feature, structure or characteristic can be combined in one or more embodiments in any suitable manner.In the following description, many specific details are provided so as to provide a full understanding of the embodiments of the present application. However, it will be appreciated by those skilled in the art that the technical scheme of the present application can be put into practice without one or more of the specific details, or other methods, components, devices, steps etc. can be adopted. In other cases, known methods, devices, implementations or operations are not shown or described in detail to avoid blurring the various aspects of the application.

[0038] The block diagrams shown in the accompanying drawings are merely functional entities and do not necessarily correspond to physically separate entities. That is, these functional entities may be implemented in software, in one or more hardware modules or integrated circuits, or in different networks and / or processor devices and / or microcontroller devices.

[0039] The flowcharts shown in the accompanying drawings are for illustrative purposes only and do not necessarily include all contents and operations / steps, nor must they be executed in the order described. For example, some operations / steps may be decomposed, while others may be combined or partially combined. Therefore, the actual execution order may vary depending on the actual situation.

[0040] It should be noted that the “plurality” mentioned in this article refers to two or more.

[0041] The following is a detailed description of the implementation details of the technical solution of the embodiment of the present application:

[0042] Example 1

[0043] This embodiment is based on the construction characteristics of power transmission projects and the accuracy requirements of the results applied to the construction of power transmission and transformation projects. Its core principles are as follows: (1) The strata at the exploration points of the same geomorphic unit are representative of this type of geomorphic unit. The geomorphic unit is further refined by combining multi-source data with factors such as geotechnical age and lithology in the area based on the geomorphic unit to further enhance the representativeness; (2) The locations of the exploration points of the power transmission project are limited due to the requirements of the tower erection conditions. Most of the towers in mountainous and hilly areas are located on the ridges and mountain tops. Therefore, the mapping area can be simply merged; (3) In areas with specific geological conditions where the construction of the project leads to drastic changes in the geomorphic unit or the cover layer, the principle of covering layer value selection is simplified to the engineering service orientation, that is, whether it affects the application and selection of rock anchors as the accuracy requirement, avoiding the problem of value assignment accuracy.

[0044] like Figure 1 As shown, this embodiment provides a method for drawing a geological distribution map of a power transmission and transformation project based on multi-source data, the method comprising a pre-processing step S101, a pre-value assignment step S102, and a value assignment calculation step S103; wherein:

[0045] The S101 preprocessing step includes taking a single regional geological map as a unit, extracting the rock and soil age genesis and lithology of the digital regional geological map based on multi-source data such as the digital geological map, terrain elevation model, and satellite imagery, performing initial zoning, correcting the zoning based on the satellite imagery and terrain elevation model, and further zoning based on the micro-geomorphological zoning.

[0046] Furthermore, in the S101 preprocessing step, the further partitioning according to the micro-geomorphology includes: considering whether to add orientation, slope, and spatial distribution factors for further partitioning based on the corresponding regional geological survey data; the attributes of the partition include geological age, geological causes, and micro-geomorphology factors.

[0047] Specifically, in this embodiment, a single regional geological map is used as a unit, and based on multi-source data such as digital geological maps, terrain elevation models, and satellite images, the digital regional geological map is extracted to initialize the zoning according to the rock and soil age and lithology, and the zoning is corrected according to the satellite images and terrain elevation models, and further zoned according to the micro-topography. Professionals consider whether to add factors such as orientation, slope, and spatial distribution for further zoning based on the corresponding regional geological survey data. The zoning attributes include geological age, geological origin, micro-topography and other factors.

[0048] S102 is a pre-assignment step, which includes projecting the collected exploration data into partitions, considering whether to filter and discard the data based on the degree of landform damage and digital geological map information, counting the coverage layers of the exploration points in the partitions, and performing assignment processing based on the partition results.

[0049] Furthermore, in the S102 pre-assignment step, the assignment processing includes: when a single partition with exploration points in the same type of partition is assigned a single value, the assignment is performed directly, and all partitions of this type are identified as assignment areas; when the range of a single partition in the same type of partition is large or the span of the cover layer thickness is large, the user reviews the micro-topography and spatial distribution factors to determine whether to add classification standards for further partitioning. If further partitioning is required, the pre-processing step is repeated; otherwise, it is identified as an interpolation area. The partition without exploration points whose classification standards are consistent with the attributes of the assignment area is identified as the assignment area, and the remaining areas without exploration points are identified as interpolation areas; the partition without exploration points whose classification standards are consistent with the attributes of the assignment area is identified as the assignment area, and the remaining areas without exploration points are identified as interpolation areas.

[0050] Specifically, in this embodiment, after the collected exploration data is projected into subareas, professionals will consider whether to filter and discard the data based on the degree of landform damage and digital geological map information. Software is then used to perform coverage statistics of the exploration points within the subareas. Based on the calculation results, professionals perform the following processing based on the subarea results:

[0051] If a single zone with exploration points in the same type of zone can be assigned a single value, the value is assigned directly, and all zones of this type are marked as assigned zones. If a single zone in the same type of zone has a large range or a large span of cover thickness, professionals should review factors such as micro-topography and spatial distribution to study whether to add classification standards for further partitioning. If further partitioning is required, the above work should be repeated. Otherwise, it is marked as an interpolation zone. The zones without exploration points that are consistent with the classification standards and the attributes of the assigned zone are marked as assigned zones, and the remaining areas without exploration points are marked as interpolation zones. Taking a sandstone area as an example, most of the cover layer is less than 2m, and there may be a cover layer thickness greater than 2m in some areas, or most of the bedrock grade is high, and some areas have a low bedrock grade. Analysis of the reasons found that the bedrock distribution has changed due to the existence of structural faults. Therefore, separate zones are divided according to the range of meters around the fault structure (linear distribution). Taking a certain alluvial area as an example, there is an obvious mutation in the thickness of the cover layer and there is regularity in space. Analysis of the reasons found that the alluvial area is bounded by a river, and the cover layer thickness is thin in the hilly area south of the river. The reason was found to be that the river has little erosion and cutting effect on the hills on the south side, resulting in a thin soft soil layer, while the north side is flat land, the original bedrock has been severely eroded and damaged, and the thickness of the soft soil is large; therefore, it is necessary to re-zone it with the river as the boundary.

[0052] S103 is a step of assigning values, which includes performing assignment calculations according to the geological age of the partitions, from old to new, and by partition category. The specific rules are as follows: first, interpolation calculations are performed on the interpolation area with exploration points, and ordinary kriging or triangular grid method is used for interpolation within the same partition. Based on the interpolation results, the partition is divided into single-valued partitions, and classified statistical analysis is performed on the spatial distribution, geographic information, micro-geomorphology and other information of the partitions to determine whether the above interpolation results can be updated through the classification standard.

[0053] Furthermore, the S103 value assignment calculation step also includes: for partitions without exploration points, re-preprocessing should be performed when the classification standard is updated; for similar assignment areas without exploration points, an inverse distance-based interpolation method is used for assignment calculation, and the inverse square of the distance between the partition to be calculated and the surrounding similar assigned areas is used as the weight for interpolation calculation. In order to reduce the amount of calculation, generally the four partitions closest to each other within a fixed distance are used for calculation; wherein the fixed distance is selected according to the accuracy of the mapping; for interpolation areas without exploration points, professional personnel perform assignment processing based on the aforementioned assignment results; for similar partitions where the thickness of the cover layer of the exploration points is irregularly distributed or the geotechnical professional analysis shows that the thickness of the cover layer of this type of partition is difficult to estimate, a comprehensive value is taken according to the standards that affect the application and selection of rock anchors.

[0054] Specifically, this embodiment performs value assignment calculations based on the geological age of the partitions, from oldest to newest, and by partition category. The specific rule is to first perform interpolation calculations on the interpolation area with exploration points. Within the same partition, ordinary kriging or triangular meshing is then used for interpolation. Kriging interpolation is a statistical interpolation method that essentially calculates using a two-dimensional continuous function to ensure that the function value of the existing data point is consistent with the actual value. Kriging interpolation is a mature calculation method. The inverse distance weighted interpolation method is based on the principle of proximity and similarity: the closer two objects are, the more similar their properties are; conversely, the farther they are, the less similar they are. It uses the distance between the interpolation point and the sample point as a weight for weighted averaging, with sample points closer to the interpolation point being assigned a greater weight.

[0055] Furthermore, based on the interpolation results, the partition is divided into single-valued partitions, and the spatial distribution, geographic information, micro-topography and other information of the partition are classified and statistically analyzed to determine whether the above interpolation results can be updated through the classification standard; for partitions without exploration points, preprocessing should be re-performed when the classification standard is updated; for similar assigned areas without exploration points, the inverse distance-based interpolation method is used for assignment calculation, and the inverse square of the distance between the partition to be calculated and the surrounding similar assigned areas is used as the weight for interpolation calculation. To reduce the amount of calculation, the four partitions closest to each other within a fixed distance are generally used for calculation; the fixed distance is selected based on the mapping accuracy.

[0056] For interpolation areas without exploration points, professionals will perform assignment processing based on the above-mentioned assignment results. In the above-mentioned assignment calculations, for areas with irregular distribution of cover thickness at exploration points in the same type of partitions, or areas where geotechnical professional analysis shows that the cover thickness of such partitions is difficult to estimate, comprehensive values ​​will be taken according to the standards that affect the application and selection of rock anchors.

[0057] For the interpolation area without exploration points, professionals will perform assignment processing based on the above-mentioned assignment results. Specifically, the core idea of ​​the 2+3+4 mapping adopted in this embodiment is that the attributes of the initial partition are mainly geological age, geological genesis, and rock and soil properties. On this basis, the attributes of the geomorphic unit are added according to satellite images and topographic maps. For the partitions with exploration points, when a single assignment is possible, it is directly assigned. If a single assignment is not possible, it is necessary to use the interpolation method for segmentation. The segmentation results are analyzed by geotechnical professionals to see whether the segmentation rules can be found from the distribution range, such as the influence of elevation and river distribution. There is a type of situation here where the distribution of cover thickness is irregular, which requires special marking, and a statistical method is used to calculate the possible maximum and minimum cover thickness of such areas. Cover thickness, average cover thickness. Since rock anchors are generally used in areas with cover thickness within 2m, 2-4m generally require testing, and areas below 4m are not considered. This type of partition can be assigned values ​​based on possible rock anchor application scenarios. For partitions without exploration points, first look for similar partitions (geological age, geological origin, lithology, geomorphic units). On this basis, the partitions are further refined into assignment areas based on the results of the previous similar partition assignments, and the inverse distance method is used for assignment. When similar partitions cannot be found in areas without exploration points, this type can only be assigned manually by geotechnical professionals using other collected geological data, and is mainly based on meeting the rock anchor selection criteria, without considering the actual value of the cover thickness. For areas with damaged landforms, there is no regularity based on the cover thickness.

[0058] Furthermore, the method further includes a checking step 104: checking and processing the above-mentioned assignment calculation results by experts, and adjusting individual areas. Specifically, the above-mentioned assignment results are checked and processed by professionals, and individual areas are adjusted;

[0059] In this embodiment, the general geological map is drawn based on the topographic map to distinguish the landform units and combine the geological survey results such as the geological age of the rock and soil body, the lithology of the fault structure, etc. to draw the geological boundaries, and finally form a structure of points, lines and surfaces, wherein the surface is mainly the geological unit body, and the line is the geological boundary or the linear distribution of the structure or the vein rock. The principle of drawing the geological distribution map in this embodiment is the same as that in this embodiment. Its basic work is to partition the region according to the topography, geological age, genesis, etc. Since the digital regional geological map has preliminarily completed this part, it is only necessary to check and verify the geological unit body (i.e., partition) based on high-precision satellite images and DEM data to see whether it needs to be corrected and merged. The final drawing is to divide the region into multiple partitions, and assign values ​​to each partition (cover layer thickness and basic rock quality grade). The contribution point that distinguishes this embodiment from the prior art is to first partition and then assign values ​​to the partitions or to interpolate and then partition, and the drawing result is a surface.

[0060] This embodiment leverages existing survey and collected data to address the conflict between the limited quantity and distribution of existing exploration point data and the complexity of geological conditions, providing a relatively reliable and versatile method suitable for the construction and management of power transmission and transformation projects. Compared to existing technologies, this method incorporates multi-source data, such as digital geological maps, terrain elevation models, and satellite imagery. The mapping area is pre-partitioned and processed based on the boundaries of these multi-source data. On this basis, methods such as the inverse distance method and kriging interpolation are then applied to interpolate the overburden thickness parameter. This effectively addresses the challenge of mapping the geological distribution of power transmission and transformation projects with a certain level of accuracy when data points are limited and unevenly distributed.

[0061] Example 2

[0062] like Figure 2 As shown, this embodiment provides a system for drawing geological distribution maps of power transmission and transformation projects based on multi-source data. The system includes a pre-processing module 10 , a pre-assignment module 20 , and an assignment calculation module 30 .

[0063] The preprocessing module 10 is used to extract the geotechnical age and lithology of the digital regional geological map based on a single regional geological map, a terrain elevation model, a satellite image and other multi-source data, perform initial zoning, perform zoning correction based on the satellite image and the terrain elevation model, and further zoning based on the micro-geomorphological zoning;

[0064] The pre-assignment module 10 is used to project the collected exploration data into the partitions, consider whether to filter and discard the data based on the degree of landform damage and digital geological map information, calculate the coverage of the exploration points in the partitions, and perform assignment processing based on the partition results;

[0065] The assignment calculation module 30 is used to perform assignment calculations according to the geological age of the partitions from old to new, and according to the partition categories. The specific rules are as follows: first, interpolation calculations are performed on the interpolation area with exploration points, and ordinary kriging or triangulated grid method is used to interpolate within the same partition range. Based on the interpolation results, the partition is divided into single-valued partitions, and the spatial distribution, geographic information, micro-topography and other information of the partitions are classified and statistically analyzed to determine whether the above interpolation results can be updated according to the classification standards.

[0066] Furthermore, the pre-processing module 10 is also used to: consider whether to add orientation, slope, and spatial distribution factors for further partitioning based on the corresponding regional geological survey data; the attributes of the partition include geological age, geological origin, and micro-geomorphological factors.

[0067] Furthermore, the pre-assignment module 20 is further configured to: directly assign a value to a single partition of the same type with an exploration point, and identify all partitions of this type as assignment areas;

[0068] When a single zone in the same type of zone has a large range or a large span of cover thickness, the user shall review the micro-topography and spatial distribution factors to determine whether to add classification criteria for further zone division. If further zone division is required, the preprocessing step shall be repeated. Otherwise, it shall be marked as an interpolation area. The zone without exploration points that is consistent with the attributes of the assignment zone according to the classification criteria shall be marked as the assignment zone, and the remaining areas without exploration points shall be marked as interpolation areas. The zone without exploration points that is consistent with the attributes of the assignment zone according to the classification criteria shall be marked as the assignment zone, and the remaining areas without exploration points shall be marked as interpolation areas.

[0069] Furthermore, the assignment calculation module 30 is further configured to: for partitions without exploration points, re-preprocess when the classification criteria are updated;

[0070] For the assigned areas of the same type without exploration points, the inverse distance-based interpolation method is used for the assignment calculation. The inverse square of the distance between the partition to be calculated and the surrounding similar assigned areas is used as the weight for interpolation calculation. To reduce the amount of calculation, the four partitions closest to each other within a fixed distance are generally used for calculation; wherein, the fixed distance is selected according to the accuracy of the mapping;

[0071] For the interpolation area without exploration points, professionals will perform the assignment process based on the above assignment results;

[0072] For similar zones where the distribution of cover thickness at exploration points is irregular or where geotechnical professional analysis shows that the cover thickness of such zones is difficult to estimate, the values ​​shall be comprehensively determined according to the standards that affect the application and selection of rock anchors.

[0073] Furthermore, the system further includes a checking module 40 for checking and processing the above-mentioned assignment calculation results according to experts, and adjusting individual areas.

[0074] This embodiment proposes an electronic device, which includes: one or more processors, and a memory, wherein the memory is used to store one or more computer programs; the computer program is configured to be executed by the one or more processors, and the program includes steps for executing the method for drawing a geological distribution map of a power transmission and transformation project based on multi-source data as described in Example 1 above.

[0075] In addition, on the other hand, this embodiment also proposes a computer-readable storage medium, characterized in that at least one instruction is stored in the storage medium, and the at least one instruction is loaded and executed by the processor to implement the steps of the method for drawing the geological distribution map of the power transmission and transformation engineering based on multi-source data as described in Example 1 above.

[0076] Those skilled in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the above description has generally described the composition and steps of each example according to function. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0077] In the several embodiments provided in this application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are merely schematic. For example, the division of the units is merely a logical function division. In actual implementation, there may be other division methods, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces, devices or units, or can be electrical, mechanical or other forms of connection.

[0078] The units described as separate components may or may not be physically separated. As a unit, a person of ordinary skill in the art can realize that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in terms of function in the above description. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0079] In addition, the functional units in the various embodiments of the present invention may be integrated into a single processing unit, each unit may exist physically separately, or two or more units may be integrated into a single unit. The aforementioned integrated units may be implemented in the form of hardware or software functional units.

[0080] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions for enabling a computer device (which can be a personal computer, server, or grid device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: various media that can store program code, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk.

[0081] The specific implementation methods described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific implementation method of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for drawing geological distribution maps of power transmission and transformation projects based on multi-source data, characterized in that: The method includes a pre-processing step, a pre-assignment step, and an assignment calculation step; wherein: The preprocessing step includes taking a single regional geological map as a unit, extracting the rock and soil age genesis and lithology of the digital regional geological map based on multi-source data including digital geological map, terrain elevation model, and satellite image, performing initial zoning, revising the zoning according to satellite image and terrain elevation model, and further zoning according to micro-relief zoning; The pre-assignment step includes projecting the collected exploration data into the partitions, considering whether to filter and discard the data based on the degree of landform damage and digital geological map information, counting the coverage layers of the exploration points in the partitions, and performing assignment processing based on the partition results; The steps of the value assignment calculation include performing value assignment calculations according to the geological age of the partitions, from old to new, and by partition category. The specific rules are as follows: first, interpolation calculations are performed on the interpolation area with exploration points, and ordinary kriging or triangular grid method is used to interpolate within the same partition. Based on the interpolation results, the partition is divided into single-valued partitions, and the spatial distribution, geographic information, micro-geomorphology and other information of the partitions are classified and statistically analyzed to determine whether the above interpolation results can be updated through the classification standard. In the pre-assignment step, the assignment process includes: In the case of a single partition with an exploration point in the same type of partition, the value is directly assigned and all partitions of this type are marked as assignment areas; When a single zone in the same type of zone has a large range or a large span of cover thickness, the user shall review the micro-topography and spatial distribution factors to determine whether to add classification criteria for further zone division. If further zone division is required, the preprocessing step shall be repeated. Otherwise, it shall be marked as an interpolation area. The zone without exploration points that is consistent with the attributes of the assignment zone according to the classification criteria shall be marked as the assignment zone, and the remaining areas without exploration points shall be marked as interpolation areas. The zone without exploration points that is consistent with the attributes of the assignment zone according to the classification criteria shall be marked as the assignment zone, and the remaining areas without exploration points shall be marked as interpolation areas.

2. The method for drawing geological distribution maps of power transmission and transformation projects based on multi-source data according to claim 1, characterized in that: The further division according to micro-topography includes: Based on the corresponding regional geological survey data, consider whether to add orientation, slope, and spatial distribution factors for further zoning; the attributes of the zoning include geological age, geological causes, and micro-geomorphological factors.

3. The method for drawing geological distribution maps of power transmission and transformation projects based on multi-source data according to claim 2, characterized in that: The assignment calculation step further includes: For partitions without exploration points, preprocessing should be performed again when the classification criteria are updated; For the assigned areas of the same type without exploration points, the inverse distance-based interpolation method is used for the assignment calculation. The inverse square of the distance between the partition to be calculated and the surrounding similar assigned areas is used as the weight for interpolation calculation. To reduce the amount of calculation, the four partitions closest to each other within a fixed distance are generally used for calculation; wherein, the fixed distance is selected according to the accuracy of the mapping; For the interpolation area without exploration points, professionals will perform the assignment process based on the above assignment results; For similar zones where the distribution of cover thickness at exploration points is irregular or where geotechnical professional analysis shows that the cover thickness of such zones is difficult to estimate, the values ​​shall be comprehensively determined according to the standards that affect the application and selection of rock anchors.

4. The method for drawing geological distribution maps of power transmission and transformation projects based on multi-source data according to claim 3 is characterized in that: The method further comprises a checking step: checking and processing the above-mentioned assignment calculation results by experts, and adjusting individual areas.

5. A system for drawing geological distribution maps of power transmission and transformation projects based on multi-source data, characterized in that: The system includes a pre-processing module, a pre-assignment module, and an assignment calculation module, and is characterized by: The preprocessing module is used to extract the geochronology and lithology of the digital regional geological map based on a single regional geological map, and initialize the zoning based on multi-source data including digital geological map, terrain elevation model, and satellite image. The zoning is corrected based on satellite image and terrain elevation model, and further zoned based on micro-geomorphological zoning. The pre-assignment module is used to project the collected exploration data into the partitions, consider whether to filter and discard the data based on the degree of landform damage and digital geological map information, calculate the coverage of the exploration points in the partitions, and perform assignment processing based on the partition results; The assignment calculation module is used to perform assignment calculations according to the geological age of the partitions from old to new, and by partition category. The specific rules are as follows: first, interpolation calculations are performed on the interpolation area with exploration points. Ordinary Kriging or triangular grid method is used for interpolation within the same partition. Based on the interpolation results, the partition is divided into single-assignment partitions, and classified statistical analysis is performed on the spatial distribution, geographic information, micro-geomorphology and other information of the partitions to determine whether the above interpolation results can be updated according to the classification standards. The pre-assignment module is further configured to: directly assign a value to a single partition of the same type where there is an exploration point and mark all partitions of this type as assignment areas; When a single zone in the same type of zone has a large range or a large span of cover thickness, the user shall review the micro-topography and spatial distribution factors to determine whether to add classification criteria for further zone division. If further zone division is required, preprocessing based on the preprocessing module shall be performed again. Otherwise, it shall be identified as an interpolation area. The zone without exploration points whose attributes are consistent with the assigned zone according to the classification criteria shall be identified as the assigned zone, and the remaining areas without exploration points shall be identified as the interpolation area. The zone without exploration points whose attributes are consistent with the assigned zone according to the classification criteria shall be identified as the assigned zone, and the remaining areas without exploration points shall be identified as the interpolation area.

6. The system for drawing geological distribution maps of power transmission and transformation projects based on multi-source data according to claim 5 is characterized in that: The pre-processing module is also used to consider whether to add orientation, slope, and spatial distribution factors for further partitioning based on the corresponding regional geological survey data; the attributes of the partition include geological age, geological origin, and micro-geomorphological factors.

7. The system for drawing geological distribution maps of power transmission and transformation projects based on multi-source data according to claim 6 is characterized in that: The assignment calculation module is further configured to: for partitions without exploration points, re-preprocess when the classification standard is updated; For the assigned areas of the same type without exploration points, the inverse distance-based interpolation method is used for the assignment calculation. The inverse square of the distance between the partition to be calculated and the surrounding similar assigned areas is used as the weight for interpolation calculation. To reduce the amount of calculation, the four partitions closest to each other within a fixed distance are generally used for calculation; wherein, the fixed distance is selected according to the accuracy of the mapping; For the interpolation area without exploration points, professionals will perform the assignment process based on the above assignment results; For similar zones where the distribution of cover thickness at exploration points is irregular or where geotechnical professional analysis shows that the cover thickness of such zones is difficult to estimate, the values ​​shall be comprehensively determined according to the standards that affect the application and selection of rock anchors.

8. The system for drawing geological distribution maps of power transmission and transformation projects based on multi-source data according to claim 7 is characterized in that: The system further comprises a checking module for checking and processing the above-mentioned assignment calculation results according to experts and making adjustments to individual areas.

Images