A power grid engineering preliminary work cost automatic measurement method

Abstract

The application discloses a kind of power grid engineering early work cost automatic measurement method, comprising: receiving three-dimensional geographic space data is carried out grid processing, generates undirected scalar cost density;Extract topographic features and with cost density is carried out tensor operation, constructs anisotropic cost impedance tensor field;Extract blocking element and construct topological boundary layer and distribute scattering cost matrix;Planning line is dispersed into differential line segment and extracts unit direction tangent vector;Through cost impedance tensor matrix and the quadratic form integral operation of tangent vector, calculate continuous type cost, and combine space crossing angle to calculate the discrete step cost of penetrating topological boundary layer;Cost data fusion, rendering output cost heat map.The application realizes the accurate quantification of topographic fluctuation direction resistance and node step cost, provides intuitive data support for engineering personnel to optimize route selection path.

Description

Technical Field

[0001] This invention relates to the field of power grid engineering cost assessment and planning technology, specifically to an automatic calculation method for the preliminary work costs of power grid engineering. Background Technology

[0002] Power grid construction is a large-scale infrastructure project involving a wide range of aspects. In the preliminary planning and design phase, accurate calculation of preliminary work costs and rational planning of route selection directly affect the overall investment efficiency and feasibility of the project. Power grid lines typically traverse complex terrains and diverse geological environments; therefore, scientifically assessing the difficulty and resource consumption of exploration operations in different spatial sections is crucial for obtaining reliable project cost data.

[0003] In existing power grid engineering cost estimation and route planning practices, cost assessment generally relies on basic environmental data provided by geographic information systems (GIS) and nationally mandated cost quota tables. Cost estimators or conventional GIS systems typically divide the planned route selection area into multiple plots or grids. Based on surface and geological parameters such as land use type, soil bearing capacity, or rock hardness within each plot, they consult engineering quota standards to statically assign a fixed basic work cost per unit area or unit length to that plot. Subsequently, by extracting the crossing length of the planned power grid line within each plot and linearly multiplying and summing it with the corresponding basic work cost, the total preliminary cost of the project is obtained.

[0004] However, existing technologies for spatial cost estimation generally employ undirected scalar cost models, which can only roughly reflect the impact of terrain through discrete terrain category coefficients. They struggle to continuously quantify the directional cost differences resulting from the coupling effect of exploration operation direction and terrain gradient. In actual engineering sites, even in areas with identical geological features, there are significant cost differences in physical resistance and resource consumption faced by engineering machinery operating laterally along contour lines versus climbing uphill perpendicular to contour lines. Existing scalar calculation methods cannot accurately distinguish the cost changes caused by different exploration operation directions, leading to discrepancies between estimated costs and actual project expenditures in complex mountainous or undulating terrain, failing to meet the technical requirements for refined cost accounting. Summary of the Invention

[0005] To address the problems mentioned in the background art, this invention provides an automatic calculation method for the preliminary work costs of power grid projects, comprising:

[0006] Step S100: The basic data discretization mapping module receives the three-dimensional geospatial data of the planned route selection area, performs grid-based discretization processing on the continuous geospatial space, generates a three-dimensional grid set, extracts the geological and soil parameters and land use parameters of the corresponding spatial coordinate positions, matches the geological features with the land attributes to convert them into a benchmark operation unit price, and assigns the benchmark operation unit price to the corresponding three-dimensional grid unit to form an undirected scalar cost density.

[0007] Step S200: The anisotropic tensor field construction module generates an anisotropic cost impedance tensor field based on the three-dimensional mesh set, extracts the three-dimensional maximum slope direction vector of the ground surface at the location of each mesh cell, and performs tensor construction operation on the three-dimensional maximum slope direction vector and the undirected scalar cost density in the corresponding mesh cell to generate a second-order cost impedance tensor matrix.

[0008] Step S300: The topological boundary collision detection module extracts the blocking element layer within the selected line area to construct the topological boundary layer, and assigns corresponding scattering cost matrices to the topological boundary layers with different attributes.

[0009] Step S400: The spatial path calculus module imports the three-dimensional planned power grid line data to be measured, discretizes the continuous line space curve into a set of differential line segments containing directional attributes, and calculates the unit direction tangent vector of each differential line segment in three-dimensional space.

[0010] Step S500 involves performing continuous spatial path integration and topological boundary collision detection to classify and calculate continuous and discrete cost data. Specifically, this includes: the spatial path calculus module extracting the cost impedance tensor matrix corresponding to the coordinates of the differential line segment, performing a quadratic multiplication operation between the cost impedance tensor matrix and the unit direction tangent vector to obtain the cost value in a specific work advancement direction, and accumulating the values ​​along the entire line to obtain the continuous preliminary work cost; the topological boundary collision detection module simultaneously monitoring the spatial relationship between the differential line segment and the topological boundary layer, and when it detects that the differential line segment penetrates the topological boundary layer, analyzing the spatial crossing angle between the line tangent vector and the boundary surface normal vector, combining the preset scattering cost matrix, and calculating the discrete step cost generated when penetrating the topological boundary layer based on the spatial crossing angle and the attribute potential difference on both sides of the topological boundary layer.

[0011] In step S600, the cost data fusion and rendering module performs an arithmetic summation operation on the continuous preliminary work cost and the discrete step cost to generate a detailed table of the total preliminary work cost of the power grid project. It also iterates through the quadratic multiplication results in the path calculation process and selects the grid cells whose values ​​are greater than the preset cost threshold to generate a cost heat map.

[0012] Preferably, in step S100, the continuous geographic space is discretized into a grid to generate a three-dimensional grid set. This includes: the three-dimensional grid partitioning unit reads the boundary parameters of the three-dimensional bounding box of the planned route selection area, obtains the longitude span, latitude span, and elevation extreme values ​​of the area, and determines the coordinates of the minimum pole point in the three-dimensional space of the bounding box; using a preset three-dimensional space side length as the partitioning step size, the space within the three-dimensional bounding box is divided into mutually orthogonal and equal-volume cubic grids to generate the three-dimensional grid set; and the specific three-dimensional content of each cubic grid is calculated by step-accumulating the coordinate axis directions of the minimum pole point coordinates in the three-dimensional space and the three-dimensional space side length. The attribute parameter mapping unit vertically projects the three-dimensional center coordinates of the cube mesh onto the geological and soil distribution vector map and the land use status vector map, reads the layer polygonal region at the location of the projection point, and extracts the soil and soil mechanical feature value and land use code carried by the layer polygonal region. When the projection of the three-dimensional center coordinates of the cube mesh coincides with the boundary line of two adjacent polygonal regions, the attribute parameter mapping unit calculates the area ratio of the cube mesh cut by multiple polygonal regions on the horizontal projection plane, and selects the soil and soil mechanical feature value and land use code corresponding to the polygonal region with the largest coverage area ratio as the unique attribute parameters of the cube mesh.

[0013] Preferably, in step S100, matching geological features with land attributes to convert them into a benchmark operating unit price, and assigning the benchmark operating unit price to the corresponding three-dimensional grid cell to form an undirected scalar cost density, includes: the attribute parameter mapping unit inputting the land type code as a query keyword into the cost quota mapping database to perform precise matching retrieval, retrieving the benchmark unit price for land acquisition compensation at the spatial location corresponding to the grid cell; inputting the geotechnical mechanical feature value as a query keyword into the cost quota mapping database to perform precise matching retrieval, retrieving the benchmark unit price for basic exploration operations at the spatial location corresponding to the grid cell; and the attribute parameter mapping unit reading the side length of the grid cell and calculating the inner line of the grid cell based on the preset route corridor width. The horizontal projected area of ​​the road corridor is calculated, and the exploration space volume of the road corridor within the grid unit is calculated based on the preset geological sampling parameters of the road cross-section. The land compensation benchmark unit price is multiplied by the horizontal projected area of ​​the road corridor to obtain the land compensation cost of the road corridor within the grid unit. The basic exploration operation benchmark unit price is multiplied by the exploration space volume to obtain the exploration operation cost of the road corridor within the grid unit. The land compensation cost and the exploration operation cost are arithmetically summed to obtain the comprehensive benchmark cost of the road corridor within the grid unit. The comprehensive benchmark cost is divided by the extension length of the road corridor within the grid unit to calculate the undirected scalar foundation cost density of the grid unit.

[0014] Preferably, in step S200, extracting the three-dimensional maximum slope direction vector of the ground surface at the location of each grid cell includes: the terrain gradient analysis unit traverses the three-dimensional grid set, extracts the center elevation value of the target grid cell and the center elevation values ​​of the adjacent grid cells adjacent to the target grid cell, and calculates the spatial difference of elevation between the target grid cell and the adjacent grid cells in the X-axis and Y-axis directions of the three-dimensional rectangular coordinate system; divides the spatial difference of elevation by the fixed side length value of the target grid cell, performs finite difference calculation in the spatial domain, and obtains the first partial derivative of the elevation change of the target grid cell in the X-axis direction and the second partial derivative of the elevation change in the Y-axis direction; sets the first partial derivative as the X-axis vector component of the three-dimensional rectangular coordinate system, and sets the second partial derivative as the Y-axis vector component of the three-dimensional rectangular coordinate system. The Z-axis dynamic component is calculated based on the root mean square of the elevation difference between adjacent grid cells within the target grid cell. This component is set as the Z-axis vector component in a three-dimensional rectangular coordinate system, ensuring that the Z-axis vector component value is zero in flat areas without elevation undulations. The three-dimensional terrain gradient vector of the target grid cell is generated based on the X-axis vector component, the Y-axis vector component, and the Z-axis vector component. The impedance tensor generation unit calls the geological resistance penalty mapping table and uses the geotechnical mechanics feature value as a matching index to extract the geological resistance directional coupling coefficient corresponding to the target grid cell. When extracting the geological resistance directional coupling coefficient, the geological resistance directional coupling coefficient corresponding to the bearing capacity of conventional leveled soil is set as a reference value of 1. When the geotechnical mechanics feature value of the target grid cell is characterized as hard rock or soft, subsiding soil, a geological resistance directional coupling coefficient greater than the reference value of 1 is extracted.

[0015] Preferably, in step S200, the tensor construction operation is performed between the three-dimensional maximum slope direction vector and the undirected scalar cost density in the corresponding grid cell to generate a second-order cost impedance tensor matrix. This includes: the impedance tensor generation unit performing scalar multiplication of the undirected scalar foundation cost density with a third-order identity matrix to calculate and obtain an isotropic foundation cost matrix; the impedance tensor generation unit performing a transpose operation on the three-dimensional terrain gradient vector to generate a corresponding three-dimensional row vector, and performing matrix multiplication of the three-dimensional terrain gradient vector with the transposed three-dimensional row vector to calculate and obtain a gradient outer product matrix.

[0016] The geological resistance directional coupling coefficient is multiplied by the gradient outer product matrix using a scalar multiplication to generate an anisotropic penalty matrix. The impedance tensor generation unit performs matrix addition on the isotropic foundation cost matrix and the anisotropic penalty matrix to construct the cost impedance tensor matrix of the target mesh element. The specific mathematical expression of the cost impedance tensor matrix is ​​as follows:

[0017]

[0018] In the formula, The cost impedance tensor matrix represents the cost of the target mesh element; The undirected scalar base cost density of the target mesh cell; Represents a third-order identity matrix; The geological resistance directional coupling coefficient of the target grid cell is a dimensionless parameter. This represents the three-dimensional terrain gradient vector of the target mesh cell; This represents the transpose of the three-dimensional terrain gradient vector.

[0019] Preferably, the process of topological boundary collision detection in steps S300 and S500 includes: the boundary layer parsing unit performing conditional filtering on the vector feature dataset using feature classification codes to separate the centerline of navigable rivers, the roadbed boundary line of highways, and the boundary lines of administrative divisions at the municipal level or above; extracting the elevation maxima and minima of the three-dimensional mesh set, and uniformly stretching the separated linear and planar data along the vertical direction within the interval from the elevation minima to the elevation maxima, converting them into a spatial three-dimensional surface model, and constructing the topological boundary layer independent of the three-dimensional mesh set; when the cross potential difference is calculated... When the unit determines that the differential line segment intersects the specific boundary surface in space, it calculates the topological intersection coordinates of the differential line segment and the specific boundary surface in three-dimensional space using a three-dimensional linear interpolation algorithm, and extracts the surface unit normal vector of the specific boundary surface at the topological intersection coordinates. Based on the coordinates of the two endpoints of the differential line segment and the direction of advancement of the planned route, it calculates the unit direction tangent vector of the differential line segment at the topological intersection coordinates. The cross potential difference calculation unit calculates the dot product value of the surface unit normal vector and the unit direction tangent vector, and obtains the spatial crossing angle when the planned route penetrates the specific boundary surface through the inverse cosine function operation.

[0020] Preferably, in step S500, calculating the discrete step cost generated when penetrating the topological boundary layer based on the spatial crossing angle and the attribute potential difference on both sides of the topological boundary layer includes: the cross potential difference calculation unit reads the geographical attribute data or administrative attribute data on both sides of the topological intersection coordinates, calls the attribute weight mapping table, and matches and converts the ecological sensitivity identifier in the geographical attribute data or the administrative level identifier in the administrative attribute data into a digital weight value to extract the first attribute weight value at the front end of the boundary surface and the second attribute weight value at the rear end of the boundary surface; calculates the absolute difference between the first attribute weight value and the second attribute weight value to generate the attribute potential difference parameter; the cross potential difference calculation unit 32 calls the corresponding scattering cost matrix according to the classification attribute of the specific boundary surface, the scattering cost matrix... The matrix is ​​a pre-constructed second-order numerical square matrix. The scattering cost matrix is ​​constructed as a positive definite diagonal matrix with the direction of the surface unit normal vector of the boundary surface as the main axis. The main diagonal elements corresponding to the normal vector direction are limited to the smallest positive real number to represent the baseline step cost, and the main diagonal elements corresponding to the tangent plane direction are limited to real numbers greater than the baseline step cost to represent the crossing angle penalty coefficient, which is used to quantify the step cost coupling strength in the unit normal vector direction. The cross potential difference calculation unit substitutes the unit direction tangent vector, the scattering cost matrix, the surface unit normal vector, and the attribute potential difference parameter into the discrete step cost calculation model, performs a multiplication operation of spatial vector and matrix, and calculates the discrete step cost generated when the planned line penetrates the topological boundary layer. The mathematical expression of the discrete step cost calculation model is:

[0021]

[0022] In the formula, Indicates that the planned route crosses the first Discrete step cost triggered when a boundary surface is encountered; Indicates the first The property potential difference parameter on both sides of the boundary surface; This represents the unit direction tangent vector of the planned route at the topological intersection point. The transpose of the unit direction tangent vector; Indicates the first The scattering cost matrix corresponding to each boundary surface; Indicates the first The surface unit normal vector of a boundary surface at the topological intersection coordinates.

[0023] Preferably, in step S400, the continuous line space curve is discretized into a set of differential line segments containing directional attributes, and the unit direction tangent vector of each differential line segment in three-dimensional space is calculated. This includes: the curve differential slicing unit takes the starting point of the three-dimensional continuous curve equation as the parameter zero point, calculates the arc length parameter sequentially according to the preset reference deviation length, and substitutes the calculated arc length parameter into the three-dimensional continuous curve equation to obtain the corresponding three-dimensional discrete node coordinates; the adjacent three-dimensional discrete node coordinates are sequentially connected to generate the discrete differential line segment set, and the value of the reference deviation length is set to be smaller than the side length of the grid cell in the three-dimensional grid set to ensure that a single differential line segment in the differential line segment set is completely located inside a single grid cell; the curve differential slicing unit extracts the starting three-dimensional coordinates and the ending three-dimensional coordinates of a single differential line segment, subtracts the starting three-dimensional coordinates from the ending three-dimensional coordinates to obtain a coordinate difference vector, divides the coordinate difference vector by the reference deviation length for numerical normalization, and generates the unit direction tangent vector of the differential line segment.

[0024] Preferably, in step S500, the cost impedance tensor matrix and the unit direction tangent vector are multiplied by a quadratic multiplication operation to obtain the cost value in a specific work advancement direction, and the values ​​are accumulated along the entire route to obtain the continuous preliminary work cost. This includes: the quadratic integration unit obtains the three-dimensional coordinates of the center point of the differential line segment, uses the three-dimensional coordinates of the center point as a spatial retrieval index, queries the target grid cell enclosing the three-dimensional coordinates of the center point in the anisotropic cost impedance tensor field, and extracts the cost impedance tensor matrix corresponding to the target grid cell; the quadratic integration unit performs a multiplication operation of the transpose of the unit direction tangent vector, the cost impedance tensor matrix, and the unit direction tangent vector to obtain the quadratic directional cost density; the quadratic directional cost density is multiplied by the baseline deviation to obtain the local cost value of a single differential line segment, and the local cost values ​​of all differential line segments in the set of differential line segments are arithmetically summed to output the continuous preliminary work cost of the planned route. The mathematical expression of the continuous cost path integral model is:

[0025]

[0026] In the formula, This represents the total cost of continuous preliminary work for the planned route. A three-dimensional continuous curve representing the planned route; The arc length parameter represents a continuous curve in three-dimensional space; The parameter representing the arc length on a continuous curve in three-dimensional space is: Tangent vector in the unit direction at; The transpose of the unit direction tangent vector; The arc length parameter is The cost impedance tensor matrix corresponding to the spatial coordinates at that location; This represents the differential parameter of the scalar arc length. In computer numerical computation, the differential parameter of the scalar arc length is specifically replaced by the baseline discrete step length, which is solved by discrete accumulation.

[0027] Preferably, in step S600, the continuous preliminary work cost and the discrete step cost are arithmetically summed to generate a detailed list of total preliminary work costs for the power grid project. The quadratic multiplication results in the path calculation process are iterated through to select grid cells with values ​​greater than a preset cost threshold, generating a cost heatmap. This includes: a total cost calculation unit using the starting point coordinates of the planned route as the baseline zero point to calculate the spatial cumulative arc length from the three-dimensional continuous curve to each differential segment and topological intersection coordinate; mapping the spatial cumulative arc length to the corresponding engineering mileage station; and spatially associating the total preliminary work cost of the power grid project according to the engineering mileage station to generate a formatted engineering cost list data table; and a three-dimensional heatmap generation unit comparing the quadratic directional cost density with the continuous cost threshold to select a set of target differential segments with values ​​greater than the continuous cost threshold, and comparing the discrete step cost with the node cost threshold to select target segments with values ​​greater than the node cost threshold. The topological intersection coordinate set; the 3D heat map generation unit calculates the difference between the quadratic directional cost density and the continuous cost threshold, and performs a normalized division operation on the difference with the preset maximum alarm cost value to obtain a first proportional coefficient. The first proportional coefficient is input into a preset RGB color mapping function to generate corresponding RGB channel values. Based on the RGB channel values, a continuous tubular color rendering material is generated for the target differential line segment set; the discrete step cost is calculated to be different from the node cost threshold, and the difference is performed a normalized division operation on the difference with the maximum alarm cost value to obtain a second proportional coefficient. The second proportional coefficient is input into the RGB color mapping function to generate a dotted color rendering material for the target topological intersection coordinate set; the 3D heat map generation unit overlays and displays the 3D digital elevation model and basic geographic base map of the planned route selection area in the rendering interface, and combines the continuous tubular color rendering material and the dotted color rendering material to output a 3D heat map showing the spatial distribution of costs.

[0028] The present invention, by adopting the above technical solution, can bring the following beneficial effects:

[0029] 1. This invention constructs an anisotropic cost resistance tensor field and generates a second-order tensor matrix by operating the three-dimensional maximum slope direction vector and the undirected scalar cost density within the grid cell. This enables the off-diagonal elements of the matrix to quantify the directional operation resistance cost determined by the terrain undulation characteristics. This solves the technical defect of traditional scalar models that are difficult to accurately distinguish the cost difference between parallel contour lines and vertical slopes. It also enables accurate calculation and automatic output of the physical cost changes caused by different exploration operation directions under the same geological conditions.

[0030] 2. This invention constructs a topological boundary layer by extracting discretely distributed blocking elements and simultaneously monitors the spatial intersection relationship of differential line segments. When penetration is detected, it analyzes the spatial crossing angle between the line tangent vector and the boundary surface normal vector. Combined with the preset scattering cost matrix and the attribute potential difference on both sides of the topological boundary layer, it is substituted into the cost model for calculation. This overcomes the defect of existing technologies that make it difficult to quantify the cost of specific nodes in the boundary area, and realizes the refined accounting of the step coordination cost generated when power grid lines cross physical barriers or policy isolation zones.

[0031] 3. This invention discretizes the continuous line spatial curve into a set of differential line segments and extracts the corresponding cost impedance tensor matrix. It uses the quadratic multiplication operation of the tensor matrix and the unit direction tangent vector to calculate the cost value in the direction of operation and accumulates it along the line trend. At the same time, it maps the total cost details after fusion calculation with the three-dimensional coordinates of the high cost grid to the geographic interactive interface for rendering. This realizes the transformation of complex working condition characteristics into a physical space cost heat map, providing intuitive data support for planners to avoid high cost areas and optimize route selection. Attached Figure Description

[0032] Figure 1 This is a flowchart illustrating the overall workflow of an automatic calculation method for preliminary work costs in power grid engineering, according to an embodiment of the present invention.

[0033] Figure 2 This is an architecture diagram of an automatic calculation system for preliminary work costs of power grid projects according to an embodiment of the present invention;

[0034] Figure 3 This is a flowchart illustrating the basic geospatial data parsing and gridding construction of one embodiment of the present invention;

[0035] Figure 4 This is a flowchart illustrating the quantitative mapping of undirected scalar basic cost density according to an embodiment of the present invention.

[0036] Figure 5 This is a flowchart illustrating the differentiation of spatial terrain gradient and extraction of geological feature coefficients according to an embodiment of the present invention;

[0037] Figure 6This is a flowchart illustrating the tensor matrix construction and anisotropic impedance coupling logic of one embodiment of the present invention;

[0038] Figure 7 This is a flowchart of a topological boundary layer extraction and collision detection algorithm according to an embodiment of the present invention;

[0039] Figure 8 This is a flowchart illustrating the scattering cost matrix model and step cost calculation according to an embodiment of the present invention;

[0040] Figure 9 This is a flowchart of the differential slicing and continuous integral quadratic form operation of the planned route according to an embodiment of the present invention;

[0041] Figure 10 This is a flowchart illustrating the total cost data fusion and 3D heatmap rendering output of one embodiment of the present invention. Detailed Implementation

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

[0043] Example 1

[0044] Please refer to the appendix. Figure 1 This invention provides an automatic calculation method for the preliminary work costs of power grid projects. This method relies on an automatic calculation system for the preliminary work costs of power grid projects. Through the interaction and calculation of multi-dimensional spatial data, it achieves the analysis and aggregation of project costs. The automatic calculation system for the preliminary work costs of power grid projects performs data interaction and business logic operations according to the following steps:

[0045] In step S100, the basic data discretization mapping module 10 receives the three-dimensional geospatial data of the planning and route selection area, performs grid discretization processing on the continuous geospatial space, and generates a three-dimensional grid set. The three-dimensional geospatial data includes digital elevation model data, geological and soil distribution data, and land use status data. The discretization process is to divide the continuous geospatial space into three-dimensional grid units of equal volume according to the set spatial coordinate resolution. The basic data discretization mapping module 10 extracts the geological and soil parameters and land use parameters of the corresponding spatial coordinate positions, calls the cost quota mapping table stored in the database, matches the geological features with the land attributes and converts them into a benchmark operation unit price, and assigns the benchmark operation unit price to the corresponding three-dimensional grid unit to form an undirected scalar cost density.

[0046] In step S200, the anisotropic tensor field construction module 20 generates an anisotropic cost impedance tensor field based on the three-dimensional mesh set. The anisotropic tensor field construction module 20 reads the elevation values ​​of each three-dimensional mesh cell, analyzes the elevation difference between adjacent mesh cells using the finite difference algorithm, and extracts the three-dimensional maximum slope direction vector of the ground surface at the location of each mesh cell. The anisotropic tensor field construction module 20 performs tensor construction operation on the three-dimensional maximum slope direction vector and the undirected scalar cost density in the corresponding mesh cell. Specifically, the undirected scalar cost density is used as the diagonal basic element of the tensor matrix, and the three-dimensional maximum slope direction vector is used as the off-diagonal penalty element of the tensor matrix to generate a second-order cost impedance tensor matrix. The diagonal elements of the tensor matrix quantify the basic cost without direction attributes, and the off-diagonal elements quantify the directional operation resistance cost determined by the terrain undulation characteristics. The tensor matrices of all three-dimensional mesh cells in the planning and selection area together constitute the anisotropic cost impedance tensor field.

[0047] In step S300, the topology boundary collision detection module 30 extracts the blocking feature layer within the selected area to construct a topology boundary layer. The topology boundary collision detection module 30 separates discrete linear or area features from the input geographic information system, specifically including navigable river boundaries, highway control lines, and administrative jurisdiction boundaries. The topology boundary collision detection module 30 converts the aforementioned separated features into a topology boundary layer independent of the three-dimensional mesh set. The topology boundary collision detection module 30 assigns corresponding scattering cost matrices to topology boundary layers with different attributes. The scattering cost matrix is ​​used to record the discrete cost association calculation rules triggered when the power grid line crosses the corresponding physical barrier or policy barrier.

[0048] In step S400, the spatial path calculus module 40 imports the three-dimensional planned power grid line data to be measured, discretizes the continuous line spatial curve into a set of differential line segments containing directional attributes, and cuts the overall spatial curve into differential line segments of a preset reference length along the geometric trend of the line design. The spatial path calculus module 40 calculates the unit direction tangent vector of each differential line segment in three-dimensional space. The unit direction tangent vector represents the direction of exploration equipment entry and the direction of surveying operation in the early stage of the project.

[0049] Step S500 involves performing continuous spatial path integration and topological boundary collision detection to classify and calculate continuous and discrete cost data. Specifically, this includes: the spatial path calculus module 40 extracting the cost impedance tensor matrix corresponding to the coordinates of the differential line segment, performing a quadratic multiplication operation between the cost impedance tensor matrix and the unit direction tangent vector to obtain the cost value in a specific work advancement direction, and accumulating the values ​​along the entire line to obtain the continuous preliminary work cost; the topological boundary collision detection module 30 synchronously monitors the spatial relationship between the differential line segment and the topological boundary layer. When it detects that the differential line segment penetrates the topological boundary layer, it analyzes the spatial crossing angle between the line tangent vector and the boundary surface normal vector, combines the preset scattering cost matrix, and calculates the discrete step cost generated when penetrating the topological boundary layer based on the spatial crossing angle and the attribute potential difference on both sides of the topological boundary layer.

[0050] In step S600, the cost data fusion and rendering module 50 fuses multi-source cost data and renders and outputs spatial distribution data of engineering costs. The cost data fusion and rendering module 50 performs an arithmetic summation operation on the continuous preliminary work costs obtained by integral accumulation and the discrete step costs calculated by collision triggering, generating a detailed table of total preliminary work costs for power grid projects. The cost data fusion and rendering module 50 traverses the quadratic multiplication operation results in the path calculation process, filters out the three-dimensional coordinates of grid units with values ​​greater than the preset cost threshold, maps the filtered three-dimensional coordinates to the three-dimensional geographic information interactive interface, assigns specific color rendering, and generates a cost heat map showing the spatial distribution of high-cost sections.

[0051] Example 2

[0052] Based on the above-described Embodiment 1, please refer to the appendix. Figure 2 This embodiment provides a detailed description of the specific architecture of the automatic calculation system for the preliminary work costs of power grid projects that implements this method. The automatic calculation system for the preliminary work costs of power grid projects is deployed in a computing device 200. The computing device 200 includes a processor 210 and a memory 220. The processor 210 is a hardware chip with concurrent data operation and matrix processing capabilities, such as a central processing unit, a graphics processing unit, or a tensor processor. The memory 220 is a non-volatile storage medium such as a random access memory, a read-only memory, or a solid-state drive. The memory 220 is used to store multi-source geographic information data and computer program instructions for performing cost calculations. The processor 210 calls the program instructions in the memory 220 to run the calculation logic.

[0053] The automatic calculation system for preliminary work costs of power grid projects includes a basic data discretization and mapping module 10, an anisotropic tensor field construction module 20, a topology boundary collision detection module 30, a spatial path calculus module 40, and a cost data fusion and rendering module 50. The basic data discretization and mapping module 10, the anisotropic tensor field construction module 20, the topology boundary collision detection module 30, the spatial path calculus module 40, and the cost data fusion and rendering module 50 combine software instructions with underlying hardware resources to achieve structured processing of spatial data.

[0054] The basic data discretization mapping module 10 is used to perform structured slicing and parameter binding on continuous geographic spatial entities. The basic data discretization mapping module 10 includes a three-dimensional grid division unit 11 and an attribute parameter mapping unit 12. The three-dimensional grid division unit 11 divides the continuous three-dimensional geographic space containing the digital elevation model into a set of grids of equal volume according to a preset spatial resolution. The attribute parameter mapping unit 12 reads the vector graphic data in the land use status layer and the geological and soil distribution layer, and extracts the soil type parameters and geological bearing capacity parameters of the corresponding coordinate positions in the grid set. In order to realize the conversion of scalar data, the attribute parameter mapping unit 12 calls the cost quota mapping database pre-stored in the memory 220, matches the corresponding basic operation unit price in the cost quota mapping database according to the soil type parameters and geological bearing capacity parameters, converts the basic operation unit price into a scalar cost density, and assigns the scalar cost density to the corresponding grid cell in the grid set.

[0055] The anisotropic tensor field construction module 20 is used to calculate and generate matrix data representing the difficulty of traversing the space. The anisotropic tensor field construction module 20 includes a terrain gradient analysis unit 21 and an impedance tensor generation unit 22. The terrain gradient analysis unit 21 uses the finite difference algorithm to analyze the elevation difference between adjacent grid cells and extracts the three-dimensional maximum slope direction vector of the grid cell. The impedance tensor generation unit 22 generates a second-order tensor matrix based on the extracted three-dimensional maximum slope direction vector and the scalar cost density in the corresponding grid cell. The diagonal elements of the tensor matrix carry the basic operation cost without directional differences, and the off-diagonal elements carry the additional operation resistance cost caused by the interaction between the three-dimensional maximum slope direction vector and the approach direction.

[0056] The topology boundary collision detection module 30 is used to identify discrete blocking elements in the planning area and quantify the node cost brought by discrete blocking elements. The topology boundary collision detection module 30 includes a boundary layer parsing unit 31 and a cross potential difference calculation unit 32. The boundary layer parsing unit 31 extracts the center line of navigable rivers, the roadbed boundary line of highways and the administrative division boundary line from the geographic information system to form a topology boundary layer independent of the grid set. The cross potential difference calculation unit 32 is used to determine the intersection state between the line space curve and the topology boundary layer, and when an intersection point is detected, it calls the preset scattering cost matrix in the memory 220 to calculate the single discrete cost required to cross the topology boundary layer.

[0057] The spatial path calculus module 40 is used to perform continuous cost calculation along the geometric trend of the route. The spatial path calculus module 40 includes a curve differential slicing unit 41 and a quadratic integration unit 42. The curve differential slicing unit 41 reads the three-dimensional coordinate data of the planned route and cuts the continuous route equation into differential line segments with directional tangent vectors. The quadratic integration unit 42 extracts the tensor matrix in the grid cell where the differential line segment is located and performs tensor quadratic multiplication operation on the tensor matrix and the directional tangent vector of the differential line segment to obtain the continuous cumulative cost of the planned route in a specific direction of advancement.

[0058] The cost data fusion and rendering module 50 is used to stitch together multi-source calculation results and generate a visual presentation of engineering data. The cost data fusion and rendering module 50 includes a total cost calculation unit 51 and a three-dimensional heat map generation unit 52. The total cost calculation unit 51 performs arithmetic summation of the continuous cost data output by the quadratic integral unit 42 and the single discrete cost output by the cross potential difference calculation unit 32 to generate detailed data of the preliminary work cost of the project. The three-dimensional heat map generation unit 52 filters out the three-dimensional coordinates of grid cells whose quadratic multiplication operation results are greater than the preset high cost threshold, and performs color rendering of the filtered grid cell three-dimensional coordinates in the three-dimensional geographic information interface.

[0059] Example 3

[0060] Based on the above embodiments one and two, the specific implementation process of each step in the automatic calculation method for the preliminary work costs of power grid projects is further explained in detail.

[0061] Please refer to the appendix. Figure 3 The basic data discretization mapping module 10 calls upon the processor hardware resources within the computing device to perform structured slicing and parameter binding on the continuous spatial entities of the planned route selection area. The basic data discretization mapping module 10, relying on its internal 3D mesh generation unit 11 and attribute parameter mapping unit 12, performs the following processing steps:

[0062] Step S101: The basic data discretization mapping module 10 receives multi-source geographic information data, including three-dimensional digital elevation model data, geological and soil distribution vector map and land use status vector map. The three-dimensional grid division unit 11 performs coordinate system alignment processing on the multi-source geographic information data and uniformly converts it to the national geodetic coordinate system or the independent spatial coordinate system specified by the project.

[0063] In step S102, the three-dimensional mesh division unit 11 performs mesh slicing on the aligned continuous three-dimensional geographic space based on the set spatial resolution. The three-dimensional mesh division unit 11 reads the boundary parameters of the three-dimensional bounding box of the planned route selection area, obtains the longitude span, latitude span, and elevation extreme values ​​of the area, and determines the coordinates of the minimum pole point of the three-dimensional bounding box. The three-dimensional mesh division unit 11 uses the preset three-dimensional space side length as the division step size to divide the space within the three-dimensional bounding box into mutually orthogonal and equal-volume cubic meshes, generating a three-dimensional mesh set. Each cubic mesh in the three-dimensional mesh set has independent three-dimensional center coordinates, which are used to uniquely identify the spatial position of the cubic mesh in physical space. The three-dimensional mesh division unit 11 calculates and generates the specific three-dimensional center coordinates of each cubic mesh by stepping and accumulating the coordinate axis directions of the minimum pole point coordinates and the three-dimensional space side length.

[0064] In step S103, the attribute parameter mapping unit 12 performs spatial intersection operation on the multi-source geographic information data, extracts the geographic attribute parameters corresponding to the cube grid, and projects the three-dimensional center coordinates of the cube grid onto the geological and soil distribution vector map and the land use status vector map vertically. The attribute parameter mapping unit 12 reads the layer polygon area at the location of the projection point, and extracts the soil and rock mechanical characteristic values ​​and land use codes carried by the layer polygon area. The soil and rock mechanical characteristic values ​​represent the soil bearing capacity or rock hardness index at the location of the cube grid, and the land use codes represent the land ownership classification and current use status at the location of the cube grid.

[0065] In step S104, the attribute parameter mapping unit 12 performs attribute cleaning on the data overlapping mesh generated in the spatial intersection operation. When the three-dimensional center coordinate projection of the cube mesh coincides with the boundary line of two adjacent polygonal regions, the attribute parameter mapping unit 12 calculates the area ratio of the cube mesh cut by multiple polygonal regions on the horizontal projection plane. The attribute parameter mapping unit 12 selects the geotechnical characteristic value and land use code corresponding to the polygonal region with the largest coverage area ratio as the unique attribute parameter of the cube mesh, so as to ensure that each cube mesh in the three-dimensional mesh set has a single and definite geographical attribute data, and eliminates multi-value conflicts in the subsequent cost matrix operation.

[0066] Further details are provided in the appendix. Figure 4The attribute parameter mapping unit 12, based on the geographic attribute parameters obtained from the aforementioned spatial intersection operation, performs the binding and quantitative transformation of numerical cost features on the grid cells in the three-dimensional grid set, specifically including the following steps:

[0067] In step S105, the attribute parameter mapping unit 12 calls the cost quota mapping database that has been pre-built and stored in the memory 220. The cost quota mapping database stores a data table of the correspondence between geographical attributes and statutory engineering cost quota values. The data table includes the benchmark unit price of land acquisition compensation mapped by different land type codes, and the benchmark unit price of basic exploration operation mapped by different geotechnical mechanical characteristic values.

[0068] Step S106: The attribute parameter mapping unit 12 traverses each grid cell in the three-dimensional grid set, extracts the land use code and geotechnical mechanical feature value attached to the grid cell, inputs the land use code as the query keyword into the cost quota mapping database to perform precise matching retrieval, retrieves the land acquisition compensation benchmark unit price of the corresponding spatial location of the grid cell, and inputs the geotechnical mechanical feature value as the query keyword into the cost quota mapping database to perform precise matching retrieval, retrieves the basic exploration operation benchmark unit price of the corresponding spatial location of the grid cell, and the basic exploration operation benchmark unit price quantitatively expresses the unit physical cost of geological survey and foundation excavation under terrain without elevation undulation.

[0069] Step S107: The attribute parameter mapping unit 12 reads the side length of the grid cell, calculates the horizontal projected area of ​​the route corridor within the grid cell based on the preset route corridor width, and calculates the exploration space volume of the route corridor within the grid cell based on the preset geological sampling parameters of the route cross-section. The attribute parameter mapping unit 12 multiplies the land acquisition compensation benchmark unit price by the horizontal projected area of ​​the route corridor to obtain the land compensation cost of the route corridor within the grid cell. The attribute parameter mapping unit 12 multiplies the basic exploration operation benchmark unit price by the exploration space volume to obtain the exploration operation cost of the route corridor within the grid cell. The attribute parameter mapping unit 12 then compares the land compensation cost with the exploration cost. The operation cost is arithmetically summed to obtain the comprehensive baseline cost of the line corridor within the grid cell. The attribute parameter mapping unit 12 divides the comprehensive baseline cost by the extension length of the line corridor within the grid cell to calculate the undirected scalar foundation cost density of the grid cell. The undirected scalar foundation cost density represents the fixed foundation cost required per unit length of the engineering line to penetrate the grid cell under static conditions that do not consider changes in terrain slope and the direction of exploration. The attribute parameter mapping unit 12 assigns the undirected scalar foundation cost density and writes it into the data structure of the grid cell as the baseline input parameter for the subsequent construction of the diagonal elements of the cost impedance tensor matrix.

[0070] Further details are provided in the appendix. Figure 5The anisotropic tensor field construction module 20, by calling the hardware resources of the processor 210 within the computing device 200, analyzes the elevation changes and geological attribute differences within the three-dimensional mesh set. Relying on its internal terrain gradient analysis unit 21 and impedance tensor generation unit 22, the anisotropic tensor field construction module 20 performs the following parameter extraction and mathematical differentiation steps:

[0071] Step S201: The terrain gradient analysis unit 21 traverses the three-dimensional mesh set, extracts the center elevation values ​​of the target mesh cell and the center elevation values ​​of adjacent mesh cells. The terrain gradient analysis unit 21 calculates the spatial difference in elevation between the target mesh cell and its adjacent mesh cells in the X-axis and Y-axis directions of the three-dimensional Cartesian coordinate system. The terrain gradient analysis unit 21 divides the spatial difference in elevation by the fixed side length of the target mesh cell, performs finite difference calculation in the spatial domain, and obtains the first partial derivative of the elevation change of the target mesh cell in the X-axis direction and the second partial derivative of the elevation change in the Y-axis direction. Yuan 21 sets the first partial derivative as the X-axis vector component of the three-dimensional rectangular coordinate system, sets the second partial derivative as the Y-axis vector component of the three-dimensional rectangular coordinate system, and calculates the Z-axis dynamic component based on the root mean square of the elevation difference between adjacent grid cells. It sets this component as the Z-axis vector component of the three-dimensional rectangular coordinate system, so that the Z-axis vector component value is zero in flat areas without elevation undulations. Based on the X-axis vector component, Y-axis vector component, and Z-axis vector component, the three-dimensional terrain gradient vector of the target grid cell is generated. The three-dimensional terrain gradient vector accurately records the spatial direction and steepness of the surface tilt of the grid cell through independent three-dimensional coordinate components.

[0072] In step S202, the impedance tensor generation unit 22 reads the geotechnical eigenvalues ​​associated with the target mesh element. The impedance tensor generation unit 22 then calls the geological resistance penalty mapping table pre-stored in the memory 220. This geological resistance penalty mapping table contains a pre-defined correspondence between rock hardness indices, soil bearing capacity classifications, and corresponding amplification coefficients. The impedance tensor generation unit 22 uses the geotechnical eigenvalues ​​as matching indices to query and compare them in the geological resistance penalty mapping table, extracting the geological resistance directional coupling coefficient corresponding to the target mesh element. The impedance tensor generation unit 22 then extracts the geological resistance directional coupling coefficient. When setting the coupling coefficient, the geological resistance direction coupling coefficient corresponding to the bearing capacity of conventional leveled soil is set to the baseline value of 1. When the geotechnical characteristic value of the target grid unit is characterized as hard rock or soft settlement soil, the impedance tensor generation unit 22 extracts the geological resistance direction coupling coefficient which is greater than the baseline value of 1. The geological resistance direction coupling coefficient is a non-negative pure numerical scalar. The value of the geological resistance direction coupling coefficient increases positively with the increase of the hardness of the underlying rock or the settlement risk coefficient of the soft soil, which is used to quantify the increase in the working resistance caused by specific geological conditions to the forward movement of engineering machinery.

[0073] In step S203, the impedance tensor generation unit 22 associates and stores the extracted three-dimensional terrain gradient vector with the geological resistance direction coupling coefficient, and writes it into the extended storage field of the corresponding target grid cell. The three-dimensional terrain gradient vector and the geological resistance direction coupling coefficient serve as pre-calculation parameters, providing independent input variables with off-diagonal elements for the subsequent construction of a cost impedance tensor matrix with directional heterogeneity.

[0074] Further details are provided in the appendix. Figure 6 The anisotropic tensor field construction module 20 calls upon the matrix operation hardware resources of the computing device 200 to construct a mathematical tensor model representing the physical resistance to spatial passage for the mesh cells in the three-dimensional mesh set based on the preprocessed multidimensional basic parameters. The anisotropic tensor field construction module 20 performs the following tensor construction steps based on the internal impedance tensor generation unit 22:

[0075] In step S204, the impedance tensor generation unit 22 reads the undirected scalar basic cost density of the target grid cell and instantiates a third-order identity matrix in the memory of the computing device 200. The impedance tensor generation unit 22 performs scalar multiplication between the undirected scalar basic cost density and the third-order identity matrix to calculate and obtain the isotropic basic cost matrix. The isotropic basic cost matrix is ​​a diagonal matrix, and the values ​​of the elements on the diagonal of the isotropic basic cost matrix are all equal to the undirected scalar basic cost density. The isotropic basic cost matrix is ​​used to quantify the inherent basic operation cost within the target grid cell without considering the directional influence of terrain undulations.

[0076] In step S205, the impedance tensor generation unit 22 reads the three-dimensional terrain gradient vector and the geological resistance directional coupling coefficient of the target mesh cell. The three-dimensional terrain gradient vector is represented as a three-dimensional column vector containing vector components of the X-axis, Y-axis and Z-axis of the three-dimensional rectangular coordinate system. The impedance tensor generation unit 22 performs a transpose operation on the three-dimensional terrain gradient vector to generate the corresponding three-dimensional row vector. The impedance tensor generation unit 22 performs matrix multiplication operation on the three-dimensional terrain gradient vector and the transposed three-dimensional row vector to calculate and obtain the gradient outer product matrix. The impedance tensor generation unit 22 performs scalar multiplication operation on the geological resistance directional coupling coefficient and the gradient outer product matrix, and then performs scalar multiplication operation on the undirected scalar foundation cost density to generate an anisotropy penalty matrix.

[0077] In step S206, the impedance tensor generation unit 22 performs matrix addition on the isotropic foundation cost matrix and the anisotropic penalty matrix to construct the cost impedance tensor matrix of the target mesh element. The specific mathematical expression of the cost impedance tensor matrix is ​​as follows:

[0078]

[0079] In the formula, The cost impedance tensor matrix representing the target mesh element; The undirected scalar base cost density of the target mesh cell; Represents a third-order identity matrix; The geological resistance directional coupling coefficient of the target grid cell is a dimensionless parameter. This represents the three-dimensional terrain gradient vector of the target mesh cell; Represents the transpose of the three-dimensional terrain gradient vector;

[0080] In step S207, the impedance tensor generation unit 22 traverses the three-dimensional mesh set and calculates the cost impedance tensor matrix of all mesh elements in the three-dimensional mesh set. The impedance tensor generation unit 22 stores the calculated cost impedance tensor matrix into the database according to the spatial coordinate mapping relationship, forming an anisotropic cost impedance tensor field covering the planned route selection area. The cost impedance tensor matrix serves as the pre-spatial spatial environment data structure, providing underlying data support for the subsequent spatial path calculus module to perform quadratic multiplication operations.

[0081] In the construction logic of the cost resistance tensor matrix, the anisotropic penalty matrix endows the cost resistance tensor matrix with direction-aware characteristics. When the direction of the engineering operation is perpendicular to the three-dimensional terrain gradient vector, the projection value of the operation direction on the gradient outer product matrix is ​​zero, and the target grid cell only outputs the undirected scalar foundation cost density on the diagonal. When the direction of the engineering operation is parallel to the three-dimensional terrain gradient vector, the projection value of the operation direction on the gradient outer product matrix deviates from zero, and the target grid cell superimposes the output of the directional operation resistance cost determined by the geological resistance direction coupling coefficient. The cost resistance tensor matrix realizes the parameter coupling of terrain gradient and geological features through matrix addition, quantifying the change in physical cost caused by different exploration operation directions under the same geological conditions.

[0082] Further details are provided in the appendix. Figure 7 The topology boundary collision detection module 30 calls upon the hardware resources of the processor 210 within the computing device 200 to identify the spatial distribution of discrete blocking elements within the planning area, and calculates the spatial topological intersection state between the planned route and the discrete blocking elements. The topology boundary collision detection module 30, relying on its internal boundary layer parsing unit 31 and cross potential difference calculation unit 32, performs the following processing steps:

[0083] Step S301: Boundary layer parsing unit 31 reads the vector feature dataset from the geographic information system. The vector feature dataset includes linear water system features, areal highway features, and administrative division boundary features. Boundary layer parsing unit 31 filters the vector feature dataset using feature classification codes to separate the centerline of navigable rivers, the roadbed boundary line of expressways, and the boundary lines of administrative divisions at the municipal level and above. Boundary layer parsing unit 31 extracts the elevation maxima and minima of the three-dimensional grid set and stretches the separated linear and areal data along the vertical direction within the range of elevation minima to elevation maxima, converting them into a spatial three-dimensional surface model. A topological boundary layer independent of the three-dimensional grid set is constructed. The topological boundary layer contains multiple sets of discrete boundary surfaces, and each set of boundary surfaces corresponds to a unique policy isolation zone attribute or physical barrier attribute.

[0084] In step S302, the cross potential difference calculation unit 32 reads the three-dimensional coordinate data of the planned route and discretizes the planned route along the geometric trend into a set of differential line segments composed of multiple nodes. The cross potential difference calculation unit 32 traverses each differential line segment in the set of differential line segments and compares the coordinates of the two endpoints of the differential line segment with the spatial position of the boundary surface in the topological boundary layer. The cross potential difference calculation unit 32 determines whether the two endpoints of the differential line segment are located on both sides of the boundary surface.

[0085] Step S303: When the cross potential difference calculation unit 32 determines that the differential line segment intersects with the specific boundary surface in space, the cross potential difference calculation unit 32 calculates the topological intersection coordinates of the differential line segment and the specific boundary surface in three-dimensional space through a three-dimensional linear interpolation algorithm. The cross potential difference calculation unit 32 extracts the surface unit normal vector of the specific boundary surface at the topological intersection coordinates. Based on the coordinates of the two endpoints of the differential line segment and the advancing direction of the planned route, the cross potential difference calculation unit 32 calculates the unit direction tangent vector of the differential line segment at the topological intersection coordinates. The surface unit normal vector represents the spatial orientation of the boundary surface, and the unit direction tangent vector represents the crossing advancing direction in the actual early exploration process of the project.

[0086] In step S304, the cross potential difference calculation unit 32 calculates the dot product of the unit normal vector and the unit direction tangent vector of the surface, and obtains the spatial crossing angle when the planned route penetrates a specific boundary surface through the inverse cosine function. The spatial crossing angle is used to quantify the geometric difference when the route crosses an obstacle vertically or obliquely. The smaller the spatial crossing angle, the closer the route is to vertically crossing the topological boundary layer, and the lower the engineering crossing difficulty and the corresponding special coordination cost. The larger the spatial crossing angle, the longer the oblique crossing distance of the route, and the higher the additional crossing cost. The spatial crossing angle, as the underlying independent variable, provides mathematical operation parameters for subsequent calculation of discrete step cost by calling the scattering cost matrix.

[0087] Further details are provided in the appendix. Figure 8 After obtaining the coordinates of the topological intersection point and the corresponding spatial vector, the cross potential difference calculation unit 32 calls the preset numerical model in the memory 220 to perform the mathematical mapping and numerical calculation of discrete step cost. The cross potential difference calculation unit 32 performs the following processing steps:

[0088] In step S401, the cross potential difference calculation unit 32 calculates the attribute potential difference parameters on both sides of the boundary surface. The cross potential difference calculation unit 32 reads the geographic attribute data or administrative attribute data on both sides of the topological intersection coordinates. The cross potential difference calculation unit 32 calls the preset attribute weight mapping table in the memory 220 to match and convert the ecological sensitivity identifier in the geographic attribute data or the administrative level identifier in the administrative attribute data into a digital weight value, so as to extract the first attribute weight value at the front end of the boundary surface and the second attribute weight value at the back end of the boundary surface. The cross potential difference calculation unit 32 calculates the absolute difference between the first attribute weight value and the second attribute weight value to generate the attribute potential difference parameter. The attribute potential difference parameter is used to quantify the incremental management and coordination difficulty or the numerical change of ecological compensation level caused by the project crossing a specific boundary surface.

[0089] In step S402, the cross potential difference calculation unit 32 retrieves the corresponding scattering cost matrix from the memory 220 according to the classification attributes of the specific boundary surface. The classification attributes of the specific boundary surface include physical barrier attribute category and policy isolation zone attribute category. The scattering cost matrix is ​​a pre-constructed second-order numerical matrix. The scattering cost matrix is ​​constructed as a positive definite diagonal matrix with the direction of the surface unit normal vector of the boundary surface as the main axis. The main diagonal elements corresponding to the normal vector direction are limited to the smallest positive real number to represent the benchmark step cost. The main diagonal elements corresponding to the tangent plane direction are limited to real numbers greater than the benchmark step cost to represent the crossing angle penalty coefficient, which is used to quantify the step cost coupling strength in the unit normal vector direction.

[0090] In step S403, the cross-potential difference calculation unit 32 substitutes the unit direction tangent vector, scattering cost matrix, surface unit normal vector, and attribute potential difference parameters into the discrete step cost calculation model. The cross-potential difference calculation unit 32 performs a multiplication operation of spatial vectors and matrices to calculate the discrete step cost generated when the planned line penetrates the topological boundary layer. The mathematical expression of the discrete step cost calculation model is:

[0091]

[0092] In the formula, Indicates that the planned route crosses the first Discrete step cost triggered when a boundary surface is encountered; Indicates the first The property potential difference parameter on both sides of the boundary surface; This represents the unit direction tangent vector of the planned route at the topological intersection point coordinates, expressed as a column vector. The transpose of the unit direction tangent vector is represented as a row vector. Indicates the first The scattering cost matrix corresponding to each boundary surface; Indicates the first The surface unit normal vector of a boundary surface at the topological intersection point coordinates is represented as a column vector.

[0093] In step S404, in the underlying operation logic of the discrete step cost calculation model, the multiplication operation of the transpose of the unit direction tangent vector, the scattering cost matrix, and the surface unit normal vector constitutes a bilinear mapping. When the unit direction tangent vector and the surface unit normal vector are collinear, it means that the planned line penetrates the boundary surface vertically. The calculated value output by the multiplication operation is the smallest, reflecting that the engineering operation cost during vertical crossing is at the baseline state. When the unit direction tangent vector deviates from the surface unit normal vector, it means that the planned line penetrates the boundary surface obliquely. The calculated value output by the multiplication operation is nonlinearly amplified under the action of the off-diagonal elements of the scattering cost matrix. This is to quantify the additional penalty cost caused by the oblique crossing leading to an increase in the crossing span and the increase in the difficulty of special design. The cross potential difference calculation unit 32 associates and binds the calculated discrete step cost with the corresponding topological intersection coordinates and stores them in an independent data table to provide node cost data support for the fusion calculation of the total cost.

[0094] Further details are provided in the appendix. Figure 9 The spatial path calculus module 40 calls upon the hardware resources of the processor 210 within the computing device 200 to map the continuous engineering line onto an anisotropic cost impedance tensor field in order to extract cost data with directional attributes. The spatial path calculus module 40, relying on its internal curve differential slicing unit 41 and quadratic integration unit 42, performs the following processing steps:

[0095] Step S501: The curve differential slice unit (41) takes the starting point of the three-dimensional continuous curve equation as the parameter zero point, calculates the arc length parameter in sequence according to the preset reference deviation step length, and substitutes the calculated arc length parameter into the three-dimensional continuous curve equation to obtain the corresponding three-dimensional discrete node coordinates; connects the adjacent three-dimensional discrete node coordinates in sequence to generate the discrete differential line segment set, and sets the value of the reference deviation step length to be smaller than the side length of the grid cell in the three-dimensional grid set to ensure that a single differential line segment in the differential line segment set is completely located inside a single grid cell;

[0096] In step S502, the curve differential slicing unit 41 traverses the set of differential line segments, extracts the starting three-dimensional coordinates and ending three-dimensional coordinates of a single differential line segment, subtracts the starting three-dimensional coordinates from the ending three-dimensional coordinates to obtain a coordinate difference vector, divides the coordinate difference vector by the reference deviation step length for numerical normalization, and generates the unit direction tangent vector of the differential line segment. The unit direction tangent vector is stored in memory as a column vector data structure to record the spatial advancement direction of engineering exploration operations on local small road segments.

[0097] Step S503: The quadratic integration unit 42 obtains the three-dimensional coordinates of the center point of the differential line segment. The quadratic integration unit 42 uses the three-dimensional coordinates of the center point as a spatial retrieval index to query the target mesh cell with the three-dimensional coordinates of the envelope center point in the anisotropic cost impedance tensor field. The quadratic integration unit 42 extracts the cost impedance tensor matrix corresponding to the target mesh cell.

[0098] In step S504, the quadratic integration unit 42 substitutes the unit direction tangent vector and the cost impedance tensor matrix into the continuous cost path integration model. The quadratic integration unit 42 performs a multiplication operation in the processor 210 of the transpose of the unit direction tangent vector, the cost impedance tensor matrix, and the unit direction tangent vector to obtain the quadratic directional cost density. The quadratic integration unit 42 multiplies the quadratic directional cost density by the baseline dispersion step size to obtain the local cost value of a single differential segment. The quadratic integration unit 42 performs an arithmetic summation of the local cost values ​​of all differential segments in the differential segment set, outputting the continuous preliminary work cost of the planned route. The mathematical expression of the continuous cost path integration model is:

[0099]

[0100] In the formula, This represents the total cost of continuous preliminary work for the planned route. A three-dimensional continuous curve representing the planned route; The arc length parameter represents a continuous curve in three-dimensional space; The parameter representing the arc length on a continuous curve in three-dimensional space is: The unit direction tangent vector at a given location is represented as a column vector. The transpose of the unit direction tangent vector is represented as a row vector. The arc length parameter is The cost impedance tensor matrix corresponding to the spatial coordinates at that location; Represents the differential parameter of scalar arc length.

[0101] In the low-level numerical computation of the computer, the scalar arc length differential parameter is specifically replaced by the baseline discrete step length, and the continuous cost path integral model is solved in the processor 210 by a discrete accumulation method.

[0102] In step S505, in the operation logic of the continuous cost path integral model, the quadratic directional cost density is obtained by multiplying the cost impedance tensor matrix by the unit directional tangent vector. When the unit directional tangent vector is parallel to the local terrain gradient direction, the projection value of the unit directional tangent vector on the off-diagonal elements of the cost impedance tensor matrix reaches a maximum value, and the local cost value output by the operation is the highest. When the unit directional tangent vector is perpendicular to the local terrain gradient direction, the projection value of the unit directional tangent vector on the off-diagonal elements of the cost impedance tensor matrix is ​​zero, and the local cost value output by the operation is the lowest. The quadratic integration unit 42 stores the calculated continuous preliminary work costs in an independent data table to provide continuous linear cost data for the data assembly of the total cost.

[0103] Further details are provided in the appendix. Figure 10 The cost data fusion and rendering module 50 calls upon the hardware resources of the processor 210 within the computing device 200 to aggregate the cost details generated in the previous calculations across various dimensions, and outputs the project cost list and spatial distribution visualization results. The cost data fusion and rendering module 50, relying on its internal total cost calculation unit 51 and 3D heatmap generation unit 52, performs the following processing steps:

[0104] Step S601: The total cost calculation unit 51 extracts the continuous preliminary work cost of the planned line recorded in the independent data table in the memory 220, as well as the discrete step cost associated with all topological intersection coordinates. The total cost calculation unit 51 performs an arithmetic summation operation on the continuous preliminary work cost and all discrete step costs to generate the total preliminary work cost of the power grid project. The total cost calculation unit 51 uses the starting point coordinate of the planned line as the reference zero point to calculate the spatial cumulative arc length from the three-dimensional continuous curve to each differential line segment and the topological intersection coordinate. The total cost calculation unit 51 maps the spatial cumulative arc length to the corresponding engineering mileage station number, and associates the total preliminary work cost of the power grid project with spatial segments according to the engineering mileage station number to generate a formatted engineering cost list data table.

[0105] In step S602, the three-dimensional heat map generation unit 52 reads the quadratic direction cost density corresponding to each differential line segment stored by the spatial path calculus module 40 during the multiplication operation, and the single-cross discrete step cost output by the cross potential difference calculation unit 32. The three-dimensional heat map generation unit 52 retrieves the preset continuous cost threshold and node cost threshold in the memory 220. The three-dimensional heat map generation unit 52 compares the quadratic direction cost density with the continuous cost threshold and filters out the target differential line segment set whose value is greater than the continuous cost threshold. The three-dimensional heat map generation unit 52 compares the discrete step cost with the node cost threshold and filters out the target topological intersection coordinate set whose value is greater than the node cost threshold.

[0106] In step S603, the 3D heat map generation unit 52 extracts the spatial coordinate interval corresponding to the target differential line segment set and the spatial position corresponding to the target topological intersection point coordinate set. The 3D heat map generation unit 52 transmits the spatial coordinate interval and spatial position parameters to the rendering engine of the 3D geographic information interactive interface. The 3D heat map generation unit 52 calculates the difference between the quadratic direction cost density and the continuous cost threshold, and performs a normalized division operation on the difference with the preset maximum alarm cost value to obtain the first proportional coefficient. The 3D heat map generation unit 52 inputs the first proportional coefficient into the preset RGB color mapping function to generate the corresponding RGB channel values. Based on the RGB channel values, a continuous tubular color rendering material is generated for the target differential line segment set. The 3D heat map generation unit 52 calculates the difference between the discrete step cost and the node cost threshold, and performs a normalized division operation on the difference with the maximum alarm cost value to obtain the second proportional coefficient. The 3D heat map generation unit 52 inputs the second proportional coefficient into the RGB color mapping function to generate a point-like color rendering material for the target topological intersection point coordinate set.

[0107] In step S604, the rendering engine of the 3D geographic information interactive interface receives continuous tubular color rendering material and point color rendering material, overlays and displays the 3D digital elevation model and basic geographic base map of the planned route selection area in the rendering interface, and outputs a 3D heat map that presents the spatial distribution of costs. The 3D heat map objectively presents the high-cost section identifiers obtained by tensor integral calculation in a 3D spatial rectangular coordinate system, and provides 3D visualization graphic data for avoiding high-cost areas and iteratively optimizing the route selection path.

Claims

1. A method for automatically calculating the preliminary work costs of power grid projects, characterized in that, include: Step S100: The basic data discretization mapping module (10) receives the three-dimensional geospatial data of the planning route selection area, performs grid-based discretization processing on the continuous geospatial space, generates a three-dimensional grid set, extracts the geological and soil parameters and land use parameters of the corresponding spatial coordinate positions, matches the geological features with the land attributes to convert them into a benchmark operation unit price, and assigns the benchmark operation unit price to the corresponding three-dimensional grid unit to form an undirected scalar cost density. Step S200: The anisotropic tensor field construction module (20) generates an anisotropic cost impedance tensor field based on the three-dimensional mesh set, extracts the three-dimensional maximum slope direction vector of the ground surface at the location of each mesh unit, and performs tensor construction operation on the three-dimensional maximum slope direction vector and the undirected scalar cost density in the corresponding mesh unit to generate a second-order cost impedance tensor matrix. Step S300, the topological boundary collision detection module (30) extracts the blocking element layer in the selected line area to construct the topological boundary layer, and assigns the corresponding scattering cost matrix to the topological boundary layer with different attributes; Step S400: The spatial path calculus module (40) imports the three-dimensional planned power grid line data to be measured, discretizes the continuous line space curve into a set of differential line segments containing directional attributes, and calculates the unit direction tangent vector of each differential line segment in three-dimensional space. Step S500: Perform continuous spatial path integration and topological boundary collision detection, and classify and calculate continuous cost data and discrete cost data. Specifically, the following steps are taken: the spatial path calculus module (40) extracts the cost impedance tensor matrix corresponding to the coordinates of the differential line segment, performs a quadratic multiplication operation between the cost impedance tensor matrix and the unit direction tangent vector, calculates the cost value in the specific operation advancement direction, and accumulates the values ​​along the entire line trend to obtain the continuous preliminary work cost; the topological boundary collision detection module (30) synchronously monitors the spatial relationship between the differential line segment and the topological boundary layer. When the differential line segment of the line is detected to penetrate the topological boundary layer, the spatial crossing angle between the line tangent vector and the boundary surface normal vector is analyzed. Combined with the preset scattering cost matrix, and based on the spatial crossing angle and the attribute potential difference on both sides of the topological boundary layer, the discrete step cost generated when penetrating the topological boundary layer is calculated. Step S600: The cost data fusion rendering module (50) performs an arithmetic summation operation on the continuous preliminary work cost and the discrete step cost to generate a detailed table of total preliminary work costs for power grid projects. It also iterates through the quadratic multiplication results in the path calculation process and selects grid units with values ​​greater than the preset cost threshold to generate a cost heat map.

2. The method for automatically calculating the preliminary work costs of power grid projects according to claim 1, characterized in that, In step S100, the continuous geographic space is discretized into a grid to generate a three-dimensional grid set, including: The three-dimensional mesh division unit (11) reads the boundary parameters of the three-dimensional bounding box of the planned route selection area, obtains the longitude span, latitude span and elevation extreme value of the area, and determines the coordinates of the minimum pole point in the three-dimensional space of the three-dimensional bounding box; using the preset three-dimensional space side length as the division step size, the space inside the three-dimensional bounding box is divided into mutually orthogonal and equal volume cube meshes, generating the three-dimensional mesh set, and calculating the specific three-dimensional center coordinates of each cube mesh by stepping and accumulating the coordinate axis directions of the minimum pole point coordinates in the three-dimensional space and the three-dimensional space side length; The attribute parameter mapping unit (12) vertically projects the three-dimensional center coordinates of the cube mesh onto the geological and soil distribution vector map and the land use status vector map, reads the layer polygon area at the location of the projection point, and extracts the soil and soil mechanical feature values ​​and land use codes carried by the layer polygon area. When the three-dimensional center coordinate projection of the cube mesh coincides with the boundary line of two adjacent polygonal regions, the attribute parameter mapping unit (12) calculates the area ratio of the cube mesh cut by multiple polygonal regions on the horizontal projection plane, and selects the geotechnical characteristic value and the land use code corresponding to the polygonal region with the largest coverage area ratio as the unique attribute parameter of the cube mesh.

3. The method for automatically calculating the preliminary work costs of power grid projects according to claim 2, characterized in that, In step S100, geological features and land attributes are matched and converted into a benchmark operation unit price, and the benchmark operation unit price is assigned to the corresponding three-dimensional grid cell to form an undirected scalar cost density, including: The attribute parameter mapping unit (12) inputs the land type code as a query keyword into the cost quota mapping database to perform precise matching retrieval, retrieves the land acquisition compensation benchmark unit price of the corresponding spatial location of the grid unit, and inputs the rock and soil mechanical characteristic value as a query keyword into the cost quota mapping database to perform precise matching retrieval, retrieves the basic exploration operation benchmark unit price of the corresponding spatial location of the grid unit; The attribute parameter mapping unit (12) reads the side length of the grid unit, calculates the horizontal projected area of ​​the line corridor within the grid unit according to the preset line corridor width, and calculates the exploration space volume of the line corridor within the grid unit according to the preset line cross-section geological sampling parameters; multiplies the land acquisition compensation benchmark unit price by the horizontal projected area of ​​the line corridor to obtain the land compensation cost of the line corridor within the grid unit, multiplies the basic exploration operation benchmark unit price by the exploration space volume to obtain the exploration operation cost of the line corridor within the grid unit; performs an arithmetic summation operation on the land compensation cost and the exploration operation cost to obtain the comprehensive benchmark cost of the line corridor within the grid unit, and divides the comprehensive benchmark cost by the extension length of the line corridor within the grid unit to calculate the undirected scalar foundation cost density of the grid unit.

4. The method for automatically calculating the preliminary work costs of power grid projects according to claim 1, characterized in that, In step S200, the direction vector of the maximum three-dimensional slope of the ground surface at the location of each grid cell is extracted, including: The terrain gradient analysis unit (21) traverses the three-dimensional mesh set, extracts the center elevation value of the target mesh cell and the center elevation value of the adjacent mesh cells adjacent to the target mesh cell, and calculates the spatial difference of elevation between the target mesh cell and the adjacent mesh cells in the X-axis and Y-axis directions of the three-dimensional rectangular coordinate system; divides the spatial difference of elevation by the fixed side length value of the target mesh cell, performs finite difference calculation in the spatial domain, and obtains the first partial derivative of the elevation change of the target mesh cell in the X-axis direction and the elevation change in the Y-axis direction. The second partial derivative is converted; the first partial derivative is set as the X-axis vector component of the three-dimensional rectangular coordinate system, the second partial derivative is set as the Y-axis vector component of the three-dimensional rectangular coordinate system, and the Z-axis dynamic component is calculated based on the root mean square of the elevation difference between adjacent grid cells. This component is set as the Z-axis vector component of the three-dimensional rectangular coordinate system, so that the Z-axis vector component value is zero in flat areas without elevation undulations. The three-dimensional terrain gradient vector of the target grid cell is generated based on the X-axis vector component, the Y-axis vector component, and the Z-axis vector component. The impedance tensor generation unit (22) calls the geological resistance penalty mapping table and uses the geotechnical mechanics feature value as the matching index to extract the geological resistance directional coupling coefficient corresponding to the target grid unit. When extracting the geological resistance directional coupling coefficient, the geological resistance directional coupling coefficient corresponding to the bearing capacity of conventional leveled soil is set as the reference value 1. When the geotechnical mechanics feature value of the target grid unit is characterized as hard rock or soft settlement soil, a geological resistance directional coupling coefficient greater than the reference value 1 is extracted.

5. The method for automatically calculating the preliminary work costs of power grid projects according to claim 4, characterized in that, In step S200, the three-dimensional maximum slope direction vector and the undirected scalar cost density within the corresponding grid cell are subjected to tensor construction operations to generate a second-order cost resistance tensor matrix, including: The impedance tensor generation unit (22) performs scalar multiplication on the undirected scalar foundation cost density and the third-order identity matrix to calculate and obtain the isotropic foundation cost matrix; The impedance tensor generation unit (22) performs a transpose operation on the three-dimensional terrain gradient vector to generate a corresponding three-dimensional row vector. It performs matrix multiplication operation on the three-dimensional terrain gradient vector and the transposed three-dimensional row vector to calculate and obtain the gradient outer product matrix. It performs scalar multiplication operation on the geological resistance directional coupling coefficient and the gradient outer product matrix to generate an anisotropic penalty matrix. The impedance tensor generation unit (22) performs matrix addition on the isotropic basic cost matrix and the anisotropic penalty matrix to construct the cost impedance tensor matrix of the target mesh element. The specific mathematical expression of the cost impedance tensor matrix is ​​as follows: In the formula, The cost impedance tensor matrix representing the target mesh element; The undirected scalar base cost density of the target mesh cell; Represents a third-order identity matrix; The geological resistance directional coupling coefficient of the target grid cell is a dimensionless parameter. This represents the three-dimensional terrain gradient vector of the target mesh cell; This represents the transpose of the three-dimensional terrain gradient vector.

6. The method for automatically calculating the preliminary work costs of power grid projects according to claim 1, characterized in that, The process involving topological boundary collision detection in steps S300 and S500 includes: The boundary layer parsing unit (31) filters the vector feature dataset by using feature classification codes to separate the center line of navigable rivers, the roadbed boundary line of expressways, and the boundary line of administrative divisions above the municipal level; it extracts the elevation maxima and minima of the three-dimensional grid set, and stretches the separated linear and planar data along the vertical direction within the range from the elevation minima to the elevation maxima, converting them into a spatial three-dimensional surface model, and constructing the topological boundary layer independent of the three-dimensional grid set; When the cross potential difference calculation unit (32) determines that the differential line segment intersects with the specific boundary surface in space, it calculates the topological intersection point coordinates of the differential line segment and the specific boundary surface in three-dimensional space through a three-dimensional linear interpolation algorithm, and extracts the surface unit normal vector of the specific boundary surface at the topological intersection point coordinates; based on the coordinates of the two endpoints of the differential line segment and the advancing direction of the planned route, it calculates the unit direction tangent vector of the differential line segment at the topological intersection point coordinates. The cross potential difference calculation unit (32) calculates the dot product value of the unit normal vector of the surface and the unit direction tangent vector, and obtains the spatial crossing angle when the planned route penetrates a specific boundary surface through the inverse cosine function operation.

7. The method for automatically calculating the preliminary work costs of power grid projects according to claim 6, characterized in that, In step S500, based on the spatial crossing angle and the potential difference of properties on both sides of the topological boundary layer, the discrete step cost generated when penetrating the topological boundary layer is calculated, including: The cross potential difference calculation unit (32) reads the geographic attribute data or administrative attribute data on both sides of the topological intersection coordinates, calls the attribute weight mapping table, matches and converts the ecological sensitivity identifier in the geographic attribute data or the administrative level identifier in the administrative attribute data into a digital weight value, so as to extract the first attribute weight value at the front end of the boundary surface and the second attribute weight value at the back end of the boundary surface; calculates the absolute difference between the first attribute weight value and the second attribute weight value, and generates the attribute potential difference parameter; The cross potential difference calculation unit (32) calls the corresponding scattering cost matrix according to the classification attribute of a specific boundary surface. The scattering cost matrix is ​​a pre-constructed second-order numerical matrix. The scattering cost matrix is ​​constructed as a positive definite diagonal matrix with the direction of the surface unit normal vector of the boundary surface as the main axis. The main diagonal elements corresponding to the normal vector direction are limited to the smallest positive real number to represent the benchmark step cost. The main diagonal elements corresponding to the tangent plane direction are limited to real numbers greater than the benchmark step cost to represent the crossing angle penalty coefficient, which is used to quantify the step cost coupling strength in the unit normal vector direction. The cross-potential difference calculation unit (32) substitutes the unit direction tangent vector, the scattering cost matrix, the surface unit normal vector, and the attribute potential difference parameter into the discrete step cost calculation model, performs a multiplication operation of spatial vectors and matrices, and calculates the discrete step cost generated when the planned line penetrates the topological boundary layer. The mathematical expression of the discrete step cost calculation model is: In the formula, Indicates that the planned route crosses the first Discrete step cost triggered when a boundary surface is encountered; Indicates the first The property potential difference parameter on both sides of the boundary surface; This represents the unit direction tangent vector of the planned route at the topological intersection point. The transpose of the unit direction tangent vector; Indicates the first The scattering cost matrix corresponding to each boundary surface; Indicates the first The surface unit normal vector of a boundary surface at the topological intersection coordinates.

8. The method for automatically calculating the preliminary work costs of power grid projects according to claim 1, characterized in that, In step S400, the continuous line space curve is discretized into a set of differential line segments containing directional attributes, and the unit direction tangent vector of each differential line segment in three-dimensional space is calculated, including: The curve differential slice unit (41) takes the starting point of the three-dimensional continuous curve equation as the parameter zero point, calculates the arc length parameter in sequence according to the preset reference deviation step length, and substitutes the calculated arc length parameter into the three-dimensional continuous curve equation to obtain the corresponding three-dimensional discrete node coordinates; connects the adjacent three-dimensional discrete node coordinates in sequence to generate the discrete differential line segment set, and sets the value of the reference deviation step length to be smaller than the side length of the grid cell in the three-dimensional grid set to ensure that a single differential line segment in the differential line segment set is completely located inside a single grid cell; The curve differential slicing unit (41) extracts the starting three-dimensional coordinates and the ending three-dimensional coordinates of a single differential line segment, subtracts the starting three-dimensional coordinates from the ending three-dimensional coordinates to obtain the coordinate difference vector, divides the coordinate difference vector by the reference deviation step length for numerical normalization, and generates the unit direction tangent vector of the differential line segment.

9. The method for automatically calculating the preliminary work costs of power grid projects according to claim 8, characterized in that, In step S500, the cost impedance tensor matrix and the unit direction tangent vector are multiplied by a quadratic multiplication to obtain the cost value in a specific work advancement direction. The cost is then accumulated along the entire route to obtain the continuous preliminary work cost, including: The quadratic integration unit (42) obtains the three-dimensional coordinates of the center point of the differential line segment, uses the three-dimensional coordinates of the center point as a spatial retrieval index, queries the target grid cell that encloses the three-dimensional coordinates of the center point in the anisotropic cost impedance tensor field, and extracts the cost impedance tensor matrix corresponding to the target grid cell. The quadratic integration unit (42) performs a multiplication operation of the transpose of the unit direction tangent vector, the cost impedance tensor matrix, and the unit direction tangent vector to obtain the quadratic directional cost density; multiplies the quadratic directional cost density by the reference deviation step size to obtain the local cost value of a single differential line segment; performs an arithmetic summation of the local cost values ​​of all differential line segments in the set of differential line segments to output the continuous preliminary work cost of the planned route. The mathematical expression of the continuous cost path integral model is: In the formula, This represents the total cost of continuous preliminary work for the planned route. A three-dimensional continuous curve representing the planned route; The arc length parameter represents a continuous curve in three-dimensional space; The parameter representing the arc length on a continuous curve in three-dimensional space is: Tangent vector in the unit direction at; The transpose of the unit direction tangent vector; The arc length parameter is The cost impedance tensor matrix corresponding to the spatial coordinates at that location; This represents the differential parameter of the scalar arc length. In computer numerical computation, the differential parameter of the scalar arc length is specifically replaced by the baseline discrete step length, which is solved by discrete accumulation.

10. The method for automatically calculating the preliminary work costs of power grid projects according to claim 1, characterized in that, In step S600, the continuous preliminary work cost and the discrete step cost are arithmetically summed to generate a detailed table of total preliminary work costs for the power grid project. The results of quadratic multiplication operations during the path calculation process are then iterated through to select grid cells with values ​​greater than a preset cost threshold, generating a cost heatmap based on their three-dimensional coordinates. This includes: The total cost calculation unit (51) uses the starting point coordinates of the planned line as the reference zero point, calculates the spatial cumulative arc length of the three-dimensional continuous curve to each differential line segment and topological intersection coordinate, maps the spatial cumulative arc length to the corresponding engineering mileage station number, and associates the total cost of the preliminary work of the power grid project with the spatial segment according to the engineering mileage station number to generate a formatted engineering cost list data table. The three-dimensional heat map generation unit (52) compares the quadratic directional cost density with the continuous cost threshold, filters out the target differential line segment set whose value is greater than the continuous cost threshold, and compares the discrete step cost with the node cost threshold, filters out the target topological intersection coordinate set whose value is greater than the node cost threshold. The three-dimensional heat map generation unit (52) calculates the difference between the quadratic directional cost density and the continuous cost threshold, and performs a normalized division operation on the difference with the preset maximum alarm cost value to obtain a first proportional coefficient. The first proportional coefficient is input into the preset RGB color mapping function to generate the corresponding RGB channel value. Based on the RGB channel value, a continuous tubular color rendering material is generated for the target differential line segment set. The unit calculates the difference between the discrete step cost and the node cost threshold, and performs a normalized division operation on the difference with the maximum alarm cost value to obtain a second proportional coefficient. The second proportional coefficient is input into the RGB color mapping function to generate a point-like color rendering material for the target topological intersection coordinate set. The three-dimensional heat map generation unit (52) overlays the three-dimensional digital elevation model and basic geographic base map of the planned route selection area in the rendering interface, and outputs a three-dimensional heat map that presents the spatial distribution of the cost by combining the continuous tubular color rendering material and the point color rendering material.

Images