Data processing method based on power transmission and transformation project geological map drawing

Abstract

The invention discloses a data processing method based on power transmission and transformation project geological map drawing, and relates to the technical field of data processing, and the method comprises the following steps: fusing the obtained multi-source geological data of a power transmission and transformation project region with planning line information, and dividing the project region into different types of geological structure units based on the fusion result; aiming at the geologic structure unit, predicting engineering geological parameters of the engineering area by adopting a space prediction algorithm corresponding to the geologic structure unit; performing engineering risk assessment based on the predicted engineering geological parameters and the power transmission and transformation project load information; and generating a power transmission and transformation project geological map according to a project risk assessment result. According to the method, the precision and reliability of geological parameter prediction can be improved, efficient fusion of multi-source data is realized, and a scientific decision basis is provided for line optimization, tower footing site selection and construction safety of a power transmission and transformation project.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of data processing, more particularly, the present application relates to a data processing method based on geological map drawing of power transmission and transformation engineering. BACKGROUND

[0002] Currently, in the field of power transmission and transformation engineering survey and design, the drawing of geological maps is crucial for line path optimization, tower foundation site selection and stability evaluation. The current mainstream drawing method has shifted from traditional manual drawing to a technical solution based on geographic information systems (GIS), that is, using general GIS software (such as ArcGIS, QGIS, etc.) for assisted drawing. Although this solution has achieved digital management and spatial visualization of geological data (such as remote sensing images, digital elevation models, survey point locations), its core processing flow still highly depends on manual intervention.

[0003] To further improve efficiency, secondary development is carried out on the general GIS platform, integrating some professional functions (such as tower foundation coordinate management, line buffer zone analysis, etc.), or building a preliminary professional geological database. This solution has optimized the workflow to some extent, but it still fails to fundamentally solve the core problem. Firstly, it lacks the ability to deeply integrate and intelligently understand multi-source heterogeneous data (especially unstructured data such as survey reports, geotechnical parameter tables), and still requires a large amount of manual data extraction and entry. Secondly, its geological analysis functions are mostly based on pre-set, fixed rule bases, lacking self-adaptation and learning ability, and unable to meet the intelligent identification and evaluation needs under complex geological conditions, making it difficult to meet the fine and dynamic evaluation requirements of power transmission and transformation engineering for geological safety. SUMMARY

[0004] In order to overcome the above-mentioned defects of the prior art, the embodiments of the present application provide a data processing method based on geological map drawing of power transmission and transformation engineering, which integrates multi-source geological data and planned line information, divides geological structure units, predicts engineering geological parameters, evaluates engineering risks, and generates a visual geological map, to solve the problems of insufficient multi-source data integration capability, limited geological parameter prediction accuracy, and difficult-to-visualize engineering risk evaluation results in the prior art.

[0005] To achieve the above-mentioned purpose, the present application provides the following technical solutions: The data processing method based on geological map drawing of power transmission and transformation engineering comprises the following steps: integrating the obtained multi-source geological data of the power transmission and transformation engineering area with the planned line information, dividing the engineering area into different types of geological structure units based on the integration result; for the geological structure units, using a spatial prediction algorithm corresponding to the geological structure units to predict the engineering geological parameters of the engineering area; based on the predicted engineering geological parameters and the power transmission and transformation engineering load information, performing engineering risk evaluation; and generating a power transmission and transformation engineering geological map according to the engineering risk evaluation result.

[0006] In a preferred embodiment, the division of the engineering area into different types of geological structural units specifically comprises: generating a structural stress field of the engineering area by processing multi-source geological data through numerical simulation; coupling the vector data of the planned route with the structural stress field for analysis to calculate the mechanical interaction index of each calculation point of the route; based on the mechanical interaction index of all calculation points, using a clustering algorithm to divide the route into sections, defining each section type obtained by clustering as a geological structural unit, and outputting a spatial distribution map of the geological structural units.

[0007] In a preferred embodiment, the mechanical interaction index specifically comprises: spatially registering the vector data of the planned route with the structural stress field; for each calculation point on the route, reading the maximum horizontal principal stress direction and the maximum horizontal principal stress amplitude at the point; according to the stress tensor data at the location, calculating the included angle between the route direction of the point and the maximum horizontal principal stress direction to obtain a direction coupling degree parameter; based on the design load of the point and the maximum horizontal principal stress amplitude, obtaining a strength coupling degree parameter; combining the direction coupling degree parameter and the strength coupling degree parameter to form the mechanical interaction index.

[0008] In a preferred embodiment, the definition method of the geological structural unit specifically comprises: for the route sections obtained by clustering, calculating the statistical characteristics of the direction coupling degree parameters and the strength coupling degree parameters of all calculation points in the sections, the statistical characteristics including the average value of the direction coupling degree parameters and the average value of the strength coupling degree parameters; based on the comparison of the statistical characteristics with preset threshold values, classifying each route section into different types of geological structural units; the types of geological structural units at least include stress synergistic stable sections and stress conflict high-risk sections.

[0009] In a preferred embodiment, for the geological structural units, a spatial prediction algorithm corresponding to the geological structural unit is used to predict the engineering geological parameters of the engineering area, specifically comprising: matching a corresponding spatial interpolation algorithm according to the type of the geological structural unit; using the known exploration point data in the target geological structural unit and adjacent units as boundary conditions for interpolation calculation to obtain an engineering geological parameter prediction map and a prediction variance map; identifying high-uncertainty areas with prediction variance higher than a preset variance threshold; for the high-uncertainty areas, automatically generating virtual exploration point data and re-interpolating until the prediction variance of all areas is lower than the preset variance threshold, and outputting an engineering geological parameter map.

[0010] In a preferred embodiment, the method for automatically generating virtual exploration point data specifically comprises: calculating a spatial gradient field of the uncertainty region based on the predicted variance map, and determining the direction in which the predicted variance increases fastest; taking the existing exploration point as the starting point, searching for a path in the high-uncertainty region along the predicted variance gradient direction, and generating an optimal path that penetrates the high-variance region; automatically laying out virtual exploration points along the optimal path at a preset interval, and assigning engineering geological parameter values to the virtual exploration points.

[0011] In a preferred embodiment, the method for performing engineering risk assessment comprises: calculating a time-varying safety factor sequence for each tower foundation point in each geological structure unit based on the predicted engineering geological parameters and the power transmission and transformation engineering load information; determining the reference risk level of each tower foundation point according to the minimum value in the time-varying safety factor sequence; and generating a static risk influence envelope by spatial simulation using an anisotropic risk attenuation model related to the main direction of the local tectonic stress field, taking each tower foundation point as a diffusion source.

[0012] In a preferred embodiment, the method for performing engineering risk assessment further comprises: superimposing the risk envelope lines of different tower foundation points in space to identify the overlapping area of the risk envelope lines as a cascading risk critical area; calculating the sum of the reference risk levels superimposed in each cascading risk critical area to obtain the superimposed risk intensity; and generating a comprehensive engineering risk assessment map based on the risk level of a single tower foundation point, the distribution of the cascading risk critical area, and the superimposed risk intensity.

[0013] In a preferred embodiment, the method for generating a power transmission and transformation engineering geological map based on the engineering risk assessment result specifically comprises: spatially registering and superimposing a geological structure unit spatial distribution map, an engineering geological parameter map, and an engineering risk assessment map to generate a comprehensive digital base map; generating a recommended path using a path search algorithm based on the comprehensive digital base map, with the starting point and the end point of the line planning as constraints, and with the minimization of the objective function composed of the superimposed risk intensity and the economic cost as the optimization goal; and obtaining the power transmission and transformation engineering geological map based on the comprehensive digital base map and the recommended path.

[0014] In a preferred embodiment, the comprehensive digital base map comprises: a simulation reinforcement operation of a user on any tower foundation point on the comprehensive digital base map; real-time updating of the reference risk level and the risk influence envelope of the tower foundation point based on the simulation reinforcement operation; and recalculating and visualizing the range and superimposed risk intensity of the cascading risk critical area related to the updated risk influence envelope.

[0015] The present application can realize accurate prediction of key engineering geological parameters under complex geological conditions by fusing multi-source geological data with planning route information and selecting corresponding spatial prediction algorithm based on geological structure unit type. Compared with the traditional method which relies on a single data source or a single interpolation method, the present application can make full use of the complementarity between different geological information, improve the spatial resolution and reliability of the prediction results, and effectively reduce the uncertainty caused by the error of geological parameter estimation to the engineering design and construction. By further combining the load information of the power transmission and transformation project on the basis of predicting the engineering geological parameters, the engineering risk is evaluated, and the evaluation results are intuitively mapped into a geological map, realizing the organic integration of data acquisition, parameter prediction, risk quantification and graphical expression. The method not only can accurately identify high-risk areas and potential engineering hazards, but also provides intuitive decision-making basis for site selection optimization, design adjustment and construction control of the power transmission and transformation project, and improves the scientificity and safety of the overall engineering construction. BRIEF DESCRIPTION OF DRAWINGS

[0016] Figure 1 FIG. 1 is a flowchart of the data processing method for drawing the geological map of the power transmission and transformation project according to the present application.

[0017] Figure 2 FIG. 2 is a flowchart of the method for predicting the engineering geological parameter map according to the present application. DETAILED DESCRIPTION

[0018] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor are within the scope of protection of the present application.

[0019] Embodiment 1, Figure 1 The data processing method for drawing the geological map of the power transmission and transformation project according to the present application is given, which comprises: S1, fuse the obtained multi-source geological data of the power transmission and transformation project area with the planning route information, and divide the engineering area into different types of geological structure units based on the fusion result; In the present embodiment, the engineering area is divided into different types of geological structure units, specifically: The multi-source geological data is processed by numerical simulation to generate the structure stress field of the engineering area, specifically: The multi-source geological data is acquired, including topography and geomorphology data, stratum lithology and geological structure data, in-situ stress measurement data and rock mass physical and mechanical parameter data, the topography and geomorphology data is acquired through satellite remote sensing or aerial photogrammetry, the distribution of regional geomorphic units (such as watershed, valley, slope) is clear, the stratum lithology and geological structure data is acquired through field geological surveying and mapping, drilling core logging and regional geological map digitization, the in-situ stress measurement data is acquired through field in-situ stress testing by hydraulic fracturing or stress relief method, and the rock mass physical and mechanical parameter data is acquired by engineering geological drilling to obtain core samples, and rock mechanics tests are carried out in the laboratory. Before numerical simulation, the multi-source geological data needs to be spatially registered and standardized to ensure that all data have a unified geographic coordinate system and data format. Common numerical simulation methods include finite element method, finite difference method and discrete element method, in the application, the finite element method or the finite difference method is preferentially adopted for simulation, firstly, based on the topography and geomorphology, stratum lithology and geological structure data in the multi-source geological data, a professional numerical simulation software (such as MIDAS GTS, FLAC3D, etc.) is used to establish a three-dimensional geomechanical model of the engineering area, the model should accurately depict the main topographic relief, the spatial distribution of different lithological strata and the geometric shape of the main faults and other structures. According to the lithology division result, the corresponding mechanical parameters are given to different rock mass units in the model, these parameters are determined by laboratory test results, a suitable constitutive model is selected for the model, and the Mohr-Coulomb elastoplasticity model is preferably selected to simulate the plastic deformation behavior of the rock mass after reaching the yield strength. The setting of the model boundary condition is the key to generating a real tectonic stress field. The bottom of the model is set as a fixed constraint, and the normal displacement constraint is applied to the side of the model. The initial stress field is superimposed by the self-weight stress field and the tectonic stress field. Firstly, the stress distribution of the model under the action of gravity is calculated by applying the gravity acceleration. Then, the equivalent tectonic force is applied to the boundary of the model (for example, the horizontal displacement or stress boundary is applied to the side of the model) to simulate the influence of regional geological structure movement. The stress state of the key points (especially the direction and size of the maximum horizontal principal stress) obtained by simulation calculation is repeatedly compared with the in-situ stress measurement data. The size and direction of the boundary tectonic force are adjusted, and the iterative calculation is carried out until the simulation result is best fitted with the measurement data within the preset error tolerance (for example, the direction error is less than 15°, and the size error is less than 20%). After the model is verified by inversion and reaches the equilibrium state, the required tectonic stress field can be extracted from the model, including the size and three-dimensional spatial vector direction of the maximum principal stress, the intermediate principal stress and the minimum principal stress.

[0020] The vector data of the planned route is coupled with the tectonic stress field for analysis, and the mechanical interaction index of each calculation point of the route is calculated, specifically: The planned route vector data undergoes a topology check to ensure it is a continuous linear feature. Subsequently, it is spatially registered with the tectonic stress field with high precision. To accurately reflect the interaction between the route and the stress field, an adaptive strategy is adopted for the layout of route calculation points, rather than simple equal-spaced placement. Basic calculation points are placed along the route centerline at a certain initial interval (e.g., 50 to 100 meters). For route curves or arcs, the spacing needs to be increased according to the radius of curvature. Preferably, when the radius of curvature R < 500m, the spacing should be shortened to 1 / 2 or even 1 / 5 of the initial interval to ensure sufficient sampling points in areas with drastic directional changes. Near known fault zones, weak interlayers, and other geologically complex areas, the spacing also needs to be increased, with intervals shortened to 10-25 meters, ultimately forming a set of route calculation points. .

[0021] The constructed stress field is essentially the stress tensor data for each grid node, and the stress tensor is symmetric. A matrix is ​​usually represented as: For each line calculation point The coordinates of the point are used to calculate the stress tensor at that point from the stress tensor data of its corresponding mesh element using bilinear interpolation. Eigenvalue decomposition is then performed on this stress tensor to obtain the direction of the maximum horizontal principal stress. and maximum horizontal principal stress amplitude Calculation points based on the line Centered on a point, take three points before and three points before it, for a total of seven points. Use the least squares method to fit a straight line, and the direction of this line is the calculation point for the route. Route at the location .

[0022] Calculate the route to this point With respect to the direction of maximum horizontal principal stress The angle between them yields the directional coupling parameter. Specifically:

[0023]

[0024] in, This represents the absolute angle difference. Directional coupling parameter. The smaller the value, the more coordinated the route direction is with the direction of ground stress.

[0025] Calculate the ratio of the design load amplitude to the maximum horizontal principal stress amplitude at this point to obtain the strength coupling parameter, specifically: Obtain line calculation points Design load of tower foundation Based on the basic design dimensions, such as the expanded foundation area converted into equivalent design stress :

[0026] intensity coupling degree parameter The calculation formula is:

[0027] Combining the directional coupling degree parameter and the intensity coupling degree parameter, a parameter group of each line calculation point is formed as a mechanical interaction index.

[0028] All mechanical interaction indexes of the line calculation points are taken as classification features, and a K-means clustering algorithm is used to divide the line into sections, and each section type obtained by clustering is defined as a geological structure unit, the geological structure unit at least includes a stress coordination stable section and a stress conflict high-risk section, which is obtained by comparing statistical features of the directional coupling degree parameter and the intensity coupling degree parameter of all calculation points in the geological structure unit with preset thresholds, the statistical features include an average value of the directional coupling degree parameter and an average value of the intensity coupling degree parameter, the preset thresholds include a first threshold, a second threshold, a third threshold, a fourth threshold and a fifth threshold, the line section with the average value of the directional coupling degree parameter less than the first threshold and the average value of the intensity coupling degree parameter between the second threshold and the third threshold is the stress coordination stable section, and the line section with the average value of the directional coupling degree parameter greater than the fourth threshold and the average value of the intensity coupling degree parameter greater than the fifth threshold is the stress conflict high-risk section. The first threshold , the second threshold , the third threshold , the fourth threshold , the fifth threshold The thresholds are set based on rock mechanics theory, underground engineering practical experience and rock mass failure criterion.

[0029] Finally, a geological structure unit spatial distribution map is output.

[0030] S2, for the geological structure unit, a spatial prediction algorithm corresponding to the geological structure unit is used to predict the engineering geological parameters of the engineering area; In the embodiment, for the geological structure unit, a spatial prediction algorithm corresponding to the geological structure unit is used to predict the engineering geological parameters of the engineering area, specifically: ​According to the type of geological structural unit, a corresponding spatial interpolation algorithm is matched, anisotropic Kriging interpolation algorithm is matched for the stress conflict high-risk section, the variogram model of the algorithm is selected as a spherical model, the main direction of the variogram of the algorithm is aligned with the principal stress direction of the local tectonic stress field, and the initial value of the anisotropy ratio (i.e. the ratio of the main direction range to the range perpendicular to the main direction) is preferably set to 2.0. The ordinary Kriging interpolation algorithm (isotropy) is configured for the stress cooperative stable section or other section, and the variogram function is also selected as a spherical model.

[0031] In the interpolation calculation of any geological structural unit, not only the known exploration point data in it is used, but also the known exploration point data in the adjacent geological structural unit is used as a boundary constraint condition to construct and participate in the construction and solution of the Kriging equation set. After the interpolation is completed, the initial engineering geological parameter prediction map of the engineering geological parameter and its corresponding prediction variance (Kriging variance) map are output, an operable variance threshold is set to identify the high uncertainty area, and a preferred embodiment of the threshold is to take the arithmetic mean of the prediction variance of all grid points in the entire engineering area plus one standard deviation. All continuous grid areas in the prediction variance map that satisfy the prediction variance greater than the variance threshold are identified as high uncertainty areas. For the high uncertainty area, virtual exploration point data is automatically generated and supplemented to the original data set, and the interpolation calculation is re-executed to generate a new prediction surface and a prediction variance map. The iteration process is repeated until the prediction variance of all areas is lower than the preset variance threshold, and the engineering geological parameter map is output.

[0032] The method of automatically generating virtual exploration point data is specifically: For the identified high uncertainty area, the spatial gradient field of the prediction variance map is calculated to determine the direction in which the prediction variance grows fastest. The existing exploration point closest to the edge of the high uncertainty area and having a prediction variance lower than the variance threshold is taken as the starting point of the path. Starting from the starting point, the Dijkstra shortest path algorithm in graph theory is used to search within the high uncertainty area. In the search process, the passing cost (Cost) of each grid cell is defined as the reciprocal of the prediction variance of the cell. The algorithm will automatically find a path with the lowest total cost. This path is essentially inclined to pass through the area with the highest prediction variance, thereby achieving the goal of "maximum range penetration of high prediction variance area". This path is determined as the optimal path for laying virtual exploration points. Along the optimal path, a series of virtual exploration points are automatically generated at a preset interval, and each generated virtual exploration point is assigned an engineering geological parameter value. The engineering geological parameter value of the virtual exploration point is determined by the following method: The engineering geological parameter value of the virtual exploration point specifically refers to the characteristic value of the foundation bearing capacity.

[0033] The trained Kriging interpolation algorithm is used to obtain the measurement value of the virtual exploration point by weighted average according to the measurement values of the known exploration points within a certain range around the virtual exploration point.

[0034] S3, based on the predicted engineering geological parameters and the power transmission engineering load information, engineering risk assessment is carried out; In the embodiment, the engineering risk assessment based on the predicted engineering geological parameters and the power transmission engineering load information is specifically: The engineering risk assessment includes: For each line tower foundation point in each geological structure unit , based on the predicted engineering geological parameters (foundation bearing capacity characteristic value) and the power transmission engineering load information of the point, the time-varying safety factor is calculated, the time-varying nature is realized by considering different working conditions, and a group of working conditions are defined , for example, Normal working condition, Rainstorm working condition, Seismic working condition, and then a time-varying safety factor sequence is obtained The calculation formula of the time-varying safety factor is:

[0035] Among them, is the foundation bearing capacity characteristic value of the line tower foundation point , and is the vertical design load.

[0036] For each tower foundation point, the minimum value in the time-varying safety factor sequence is taken as the basis for risk assessment, and the risk level threshold , is set to determine the reference risk level of the point , which is specifically:

[0037] Taking each tower foundation point as a diffusion source, the reference risk level of the tower foundation point is defined as the initial risk intensity of the source point , and an anisotropic risk attenuation model is established, the core of the model is to describe how the risk influence intensity decays with the distance d and direction θ away from the diffusion source point, the anisotropy refers to that the decay rate is related to the direction, and the elliptical decay model is used in the model, the long axis direction of the ellipse is the direction of the maximum horizontal principal stress at the tower foundation point Consistency indicates that the risk's influence extends further (attenuates more slowly) along the principal direction of tectonic stress. The minor axis being perpendicular to the major axis indicates that the risk attenuates faster in that direction. The intensity of the risk's influence in any direction θ... The decay formula with distance d can be expressed as:

[0038] Where m is the base index, specifying which base point the current calculation is for, and d is the distance between the calculation point and the base. The Euclidean distance between them, where θ is the calculation point relative to the base point of the tower. azimuth angle, This is a function of attenuation coefficients related to direction and tower base location.

[0039] The calculation is related to elliptic geometry, specifically:

[0040] in, Base point of the tower The direction of the maximum horizontal principal stress at that location, This represents the characteristic attenuation distance along the direction of maximum horizontal principal stress (major axis). The characteristic attenuation distance is perpendicular to the direction of the maximum horizontal principal stress (minor axis) and satisfies... .

[0041] Set a lower threshold for the intensity of risk impact. The static risk impact envelope is all conditions that satisfy the following: A closed figure formed by points on a plane.

[0042] Spatial superposition analysis of risk envelopes at different base points is performed to identify areas where risk envelopes overlap, which are defined as critical cascading risk zones. The sum of the baseline risk levels corresponding to the superimposed risk envelopes within each critical cascading risk zone is calculated to obtain the superimposed risk intensity of the zone. The larger the superimposed risk intensity value, the higher the potential risk of cascading failure in the zone.

[0043] Based on the risk level of a single tower base point, the distribution of critical risk zones in a chain reaction, and the intensity of superimposed risks, a comprehensive engineering risk assessment map is generated.

[0044] S4 generates a geological map of the power transmission and transformation project based on the project risk assessment results.

[0045] In this embodiment, generating a geological map of the power transmission and transformation project based on the engineering risk assessment results specifically involves: The process of generating a geological map of the power transmission and transformation project based on the engineering risk assessment results is as follows: The geological structure unit spatial distribution map, the engineering geological parameter map, and the engineering risk assessment atlas are spatially registered and superimposed, and are managed as different layers. For example, the geological structure unit map is taken as a base layer, the engineering geological parameter contour is taken as a semi-transparent cover layer, and the risk assessment atlas (especially the key area of chain risk) is superimposed on the uppermost layer in a high-light color. Thus, the integrated digital map containing multiple geological and risk information is the comprehensive digital base map, Based on the comprehensive digital base map, the starting point and the ending point of the line planning are taken as constraints, the minimization of the target function composed of the superimposed risk intensity and the economic cost of the path through area is taken as the optimization goal, and a path search algorithm (such as genetic algorithm) is used to generate one or more recommended paths. The comprehensive digital base map comprises: The user performs a simulation reinforcement operation on any tower foundation point on the comprehensive digital base map According to the simulation reinforcement operation, the reference risk level and the risk influence envelope of the tower foundation point are updated in real time, the range and the superimposed risk intensity of the key area of chain risk related to the updated risk influence envelope are recalculated and visualized, and the range and the superimposed risk intensity of the key area of chain risk related to the updated risk influence envelope are updated.

[0046] The recommended path generated by the intelligent optimization is synthesized with the comprehensive digital base map, drawing elements such as legends and scales are added, and the final power transmission and transformation engineering geological map is output.

[0047] The above formulas are all dimensionless numerical calculations, the formulas are obtained by software simulation of a large amount of collected data to obtain a formula closest to the real situation, and the preset parameters in the formula are set by a person skilled in the art according to the actual situation.

[0048] The above embodiments can be realized wholly or partially by software, hardware, firmware, or any combination thereof. When realized by software, the above embodiments can be realized in the form of a computer program product wholly or partially.

[0049] Those skilled in the art can realize that the modules and algorithm steps of the examples described in combination with the embodiments disclosed herein can be realized by electronic hardware or a combination of computer software and electronic hardware. Whether the functions are realized in hardware or software depends on the specific application and design constraints of the technical solutions. A person skilled in the art can use different methods to realize the described functions for each specific application, but such implementation should not be considered beyond the scope of the present application.

[0050] In addition, the functional modules in each embodiment of the present application can be integrated in one processing module, or each module can be physically present alone, or two or more modules can be integrated in one module.

[0051] The above merely provides the specific implementation of the present application, but the protection scope of the present application is not limited to this. Any person skilled in the art can easily think of the changes or replacements within the technical range disclosed by the present application, which should be included in the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

[0052] Finally, the above merely provides the preferred embodiments of the present application, but is not used to limit the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the protection scope of the present application.

Claims

1. A data processing method based on a power transmission and transformation engineering geological map drawing, characterized in that, The method comprises the following steps: fusing the obtained multi-source geological data of the power transmission and transformation project area with the planned line information, dividing the project area into different types of geological structure units based on the fusion result; for the geological structure units, using a spatial prediction algorithm corresponding to the geological structure units to predict the engineering geological parameters of the project area; based on the predicted engineering geological parameters and the power transmission and transformation project load information, performing engineering risk assessment; generating a power transmission and transformation project geological map according to the engineering risk assessment result.

2. The data processing method based on the geological mapping of power transmission projects according to claim 1, characterized in that, The division of the project area into different types of geological structure units is specifically: processing the multi-source geological data through numerical simulation to generate a structural stress field of the project area; coupling the vector data of the planned line with the structural stress field to calculate the mechanical interaction index of each calculation point of the line; based on the mechanical interaction index of all calculation points, using a clustering algorithm to divide the line into sections, defining each section type obtained by clustering as a geological structure unit, and outputting a spatial distribution map of the geological structure units.

3. The data processing method based on the geological mapping of power transmission projects according to claim 2, characterized in that, The mechanical interaction index is specifically: spatially registering the vector data of the planned line with the structural stress field; for each calculation point on the line, reading the maximum horizontal principal stress direction and the maximum horizontal principal stress amplitude at the point; according to the stress tensor data of the position, calculating the included angle between the line direction of the point and the maximum horizontal principal stress direction to obtain a direction coupling degree parameter; based on the design load of the point and the maximum horizontal principal stress amplitude, obtaining a strength coupling degree parameter; combining the direction coupling degree parameter and the strength coupling degree parameter to form the mechanical interaction index.

4. The data processing method based on the power transmission engineering geological mapping according to claim 3, characterized in that, The definition method of the geological structure unit is specifically: for the line sections obtained by clustering, calculating the statistical characteristics of the direction coupling degree parameters and the strength coupling degree parameters of all calculation points in the sections, the statistical characteristics including the average value of the direction coupling degree parameters and the average value of the strength coupling degree parameters; based on the comparison of the statistical characteristics with preset threshold values, classifying the line sections into different types of geological structure units; the types of geological structure units at least include stress coordination stable sections and stress conflict high-risk sections.

5. The data processing method based on the power transmission engineering geological mapping according to claim 4, characterized in that, For the geological structure units, using a spatial prediction algorithm corresponding to the geological structure units to predict the engineering geological parameters of the project area is specifically: matching a corresponding spatial interpolation algorithm according to the type of the geological structure unit; using the known exploration point data in the target geological structure unit and the adjacent units as boundary conditions to perform interpolation calculation to obtain an engineering geological parameter prediction map and a prediction variance map; identifying a high-uncertainty area with a prediction variance higher than a preset variance threshold; for the high-uncertainty area, automatically generating virtual exploration point data and re-interpolating until the prediction variance of all areas is lower than the preset variance threshold, and outputting an engineering geological parameter map.

6. The data processing method based on the power transmission engineering geological mapping according to claim 5, characterized in that, The method for automatically generating virtual exploration point data is specifically: based on the prediction variance map, calculating the spatial gradient field of the uncertainty area to determine the direction in which the prediction variance grows fastest; taking an existing exploration point as a starting point, searching for an optimal path that penetrates the high-variance area in the high-uncertainty area along the direction of the prediction variance gradient; and The virtual exploration points are automatically arranged along the optimal path at preset intervals, and the engineering geological parameter values are assigned to the virtual exploration points.

7. The data processing method based on the power transmission engineering geological mapping according to claim 6, characterized in that, The engineering risk assessment includes: For each tower foundation point in each geological structure unit, a time-varying safety factor sequence is calculated based on the predicted engineering geological parameters and the power transmission and transformation engineering load information; According to the minimum value in the time-varying safety factor sequence, the reference risk level of each tower foundation point is determined; Taking each tower foundation point as a diffusion source, an anisotropic risk attenuation model related to the main direction of the local tectonic stress field is used for spatial simulation to generate a static risk influence envelope.

8. The data processing method based on the power transmission engineering geological mapping according to claim 7, characterized in that, The engineering risk assessment also includes: The risk envelope lines of different tower foundation points are spatially superimposed, and the overlapping area of the risk envelope lines is identified as a cascading risk critical area; The sum of the superimposed reference risk levels in each cascading risk critical area is calculated to obtain the superimposed risk intensity; Based on the risk level of a single tower foundation point, the distribution of the cascading risk critical area, and the superimposed risk intensity, a comprehensive engineering risk assessment map is generated.

9. The data processing method based on the power transmission engineering geological mapping according to claim 8, characterized in that, According to the engineering risk assessment results, a power transmission and transformation engineering geological map is generated, specifically: The geological structure unit spatial distribution map, the engineering geological parameter map, and the engineering risk assessment map are spatially registered and superimposed to generate a comprehensive digital base map; Based on the comprehensive digital base map, the starting point and the ending point of the line planning are constrained, the superimposed risk intensity and the economic cost constitute the objective function minimization as the optimization goal, and the path search algorithm is used to generate a recommended path; Based on the comprehensive digital base map and the recommended path, a power transmission and transformation engineering geological map is obtained.

10. The data processing method based on the power transmission engineering geological mapping according to claim 9, characterized in that, The comprehensive digital base map includes: In response to the user's simulation reinforcement operation on any tower foundation point on the comprehensive digital base map; Based on the simulation reinforcement operation, the reference risk level of the tower foundation point and its risk influence envelope are updated in real time; The range and superimposed risk intensity of the cascading risk critical area related to the updated risk influence envelope are recalculated and visualized.

Images