Earthwork analysis method and system based on comprehensive geophysical prospecting

Abstract

The invention relates to the technical field of geotechnical engineering geological investigation, in particular to an earthwork analysis method and system based on comprehensive geophysical prospecting, and the method comprises the steps: obtaining basic topographic data, and obtaining the stratum elevation data of a to-be-analyzed airport through a geophysical prospecting method; carrying out terrain design on the to-be-analyzed airport, and obtaining excavation and filling elevation data of the to-be-analyzed airport; basic topographic data, stratum elevation data obtained by geophysical prospecting and excavation and filling elevation data after terrain design are imported into an ArcGIS platform, and triangulated irregular network (TIN) models are generated respectively; constructing a digital elevation model (DEM) according to the triangulated irregular network (TIN); calculating the earth-rock volume in combination with the shrinkage-expansion square coefficient; and dynamically adjusting a digital elevation model (DEM) and an earth-rock volume calculation result according to actual conditions in the construction process. According to the method, the stratum information and the digging and filling coefficient are comprehensively considered, so that the earth-rock calculation result is more scientific and reasonable, three-dimensional display and dynamic adjustment can be performed on the conditions before and after earth-rock engineering construction, the visualization effect is stronger, and calculation is flexible.

Description

Technical Field

[0001] The present application relates to the technical field of geotechnical engineering geological survey, and in particular to an earthwork analysis method and system based on comprehensive geophysical prospecting. Background Art

[0002] Mountain airports have complex terrain. To ensure clear airspace for aircraft, obstacles near the proposed airport, especially on the take-off and landing routes, must be cleared. Earthwork is crucial for mountain airports, and its construction costs may exceed 60% of the airport construction costs. In addition, due to the different cut-and-fill ratios of earthwork and stonework, and the large difference in construction costs between earthwork and stonework, the calculation and cost control of earthwork and stonework are greatly affected. Therefore, the rational calculation of earthwork and stonework volume for each rock and soil type (cover layer, strongly weathered rock and soil layer, and moderately weathered rock and soil layer) is particularly important for the construction of mountain airport projects.

[0003] However, traditional methods such as the grid method and the cross-section method have low calculation accuracy under complex terrain conditions (such as high-fill airports in mountainous areas). The excavation and filling ratios of earthwork and stonework are different, and there are also large differences in construction costs. Traditional methods fail to fully consider this difference, resulting in inaccurate project cost estimates, which has a significant impact on cost control. In addition, the terrain of mountain airports is very undulating, and drilling work is extremely difficult. Without drilling to understand the lithology of the strata in the project area, the presence of weathered layers and overburden layers has a significant impact on the calculation of earthwork quantities. During the construction process, earthwork will change in volume due to factors such as compaction and expansion. Traditional calculation methods do not take into account the shrinkage coefficient and cannot accurately reflect the actual project quantity, which directly affects the site selection and construction decision-making of mountain airports. Summary of the Invention

[0004] To achieve the above objectives, this application provides the following technical solutions:

[0005] According to a first aspect of the present invention, the present invention claims protection for a method for analyzing earthwork based on comprehensive geophysical prospecting, comprising:

[0006] Obtain basic terrain data and use geophysical methods to obtain the ground elevation data of the airport to be analyzed;

[0007] Performing terrain design on the airport to be analyzed and obtaining cut-and-fill elevation data of the airport to be analyzed;

[0008] Import basic terrain data, stratum elevation data obtained from geophysical exploration, and cut-and-fill elevation data after terrain design into the ArcGIS platform, and generate irregular triangulated network (TIN) models for each.

[0009] Constructing a digital elevation model (DEM) based on the triangulated irregular network (TIN);

[0010] In combination with the shrinkage coefficient, the earthwork volume is calculated based on the area of each pixel in the digital elevation model (DEM) and the corresponding elevation change;

[0011] The construction is carried out based on the earthwork volume and according to the actual situation during the construction process, the digital elevation model DEM and the earthwork volume calculation result are dynamically adjusted.

[0012] Furthermore, the method of obtaining basic terrain data and obtaining stratum elevation data of the airport to be analyzed by using geophysical exploration methods further includes:

[0013] Arrange a physical detection line in a local area of the airport to be analyzed, and use seismic exploration and geological radar methods to obtain layered data of the weathering layer and the overburden layer;

[0014] Obtaining the resistivity distribution of the weathered layer and the cover layer, and combining the on-site geological analysis of the airport to be analyzed to obtain the soil composition, resistivity, and thickness of the cover layer and the weathered layer of the airport to be analyzed;

[0015] A first local terrain elevation point is generated based on the elevation of the measuring point of the physical detection line and the terrain data measured on site.

[0016] Furthermore, the performing of terrain design on the airport to be analyzed and obtaining elevation data of cut and fill of the airport to be analyzed further includes:

[0017] Based on the buried depth of the cover layer and the severe weathering limit determined by the geophysical exploration results, the cover layer thickness and the severe weathering limit thickness of the area of the airport to be analyzed are determined by spatial interpolation, and the second elevation point of the cover layer and the third elevation point of the severe weathering limit layer are prepared;

[0018] Determine the excavation elevation point through terrain design.

[0019] Furthermore, the step of importing basic terrain data, stratum elevation data obtained through geophysical exploration, and cut-and-fill elevation data after terrain design into the ArcGIS platform to generate triangulated irregular network (TIN) models for each of them further includes:

[0020] The original terrain TIN model is constructed using the first local terrain elevation point, the cover layer boundary TIN model is constructed using the second elevation point, the weathering boundary TIN model is constructed using the third elevation point, and the design elevation TIN model is determined using the terrain design cut elevation.

[0021] Furthermore, the construction of a digital elevation model (DEM) based on the triangulated irregular network (TIN) further includes:

[0022] The original terrain TIN model is converted into an original terrain DEM model, the cover boundary TIN model is converted into a cover boundary DEM model, the weathering boundary TIN model is converted into a weathering boundary DEM model, and the design elevation TIN model is converted into a design elevation DEM model.

[0023] Furthermore, the calculation of earthwork volume based on the area of each pixel in the digital elevation model (DEM) and the corresponding elevation change in combination with the shrinkage coefficient also includes:

[0024] The earthwork volume is calculated based on the following formula:

[0025] V = [cell (area)] × ΔZ;

[0026] Where: V is the calculated earthwork volume m 3 ; cell(area) is the area of a pixel to be calculated in DEM data m 2 ; ΔZ is the elevation change m corresponding to the calculated pixel before and after the proposed design.

[0027] Furthermore, the construction based on the earthwork volume and dynamically adjusting the digital elevation model DEM and earthwork volume calculation results according to the actual situation during the construction process also includes:

[0028] Based on the original terrain DEM model and the above-mentioned cover layer boundary DEM model, ArcGIS cut and fill analysis is used to calculate the cover layer cut volume;

[0029] Based on the cover boundary DEM model and the weathering boundary DEM model, ArcGIS cut and fill analysis is used to calculate the excavation volume of the strongly weathered rock layer;

[0030] Based on the weathering limit DEM model and the design elevation DEM model, the cut and fill volume of the medium-weathered rock layer is calculated using ArcGIS cut and fill analysis;

[0031] The volume of abandoned earth and filled earth of the airport to be analyzed is obtained based on the volume of excavation of the overburden layer, the volume of excavation of the strongly weathered rock layer and the volume of excavation of the moderately weathered rock layer, combined with the expansion coefficients of different rock layers.

[0032] According to a second aspect of the present invention, the present invention claims protection for an earthwork analysis system based on comprehensive geophysical exploration, comprising:

[0033] one or more processors;

[0034] A memory having one or more programs stored thereon, wherein when the one or more programs are executed by the one or more processors, the one or more processors implement the earthwork analysis method based on comprehensive geophysical exploration.

[0035] The present application relates to the field of geotechnical engineering geological survey technology, and in particular to a method and system for analyzing earthwork based on comprehensive geophysical prospecting, which includes obtaining basic terrain data, obtaining stratum elevation data of an airport to be analyzed by using geophysical prospecting methods, performing terrain design on the airport to be analyzed, and obtaining cut and fill elevation data of the airport to be analyzed, importing basic terrain data, stratum elevation data obtained by geophysical prospecting, and cut and fill elevation data after terrain design into an ArcGIS platform, and generating an irregular triangulated network (TIN) model for each, constructing a digital elevation model (DEM) based on the irregular triangulated network (TIN), calculating earthwork volume by combining shrinkage and expansion coefficients, and dynamically adjusting the digital elevation model (DEM) and earthwork volume calculation results according to the actual situation during the construction process. The present invention comprehensively considers stratum information and cut and fill coefficients, making the earthwork calculation results more scientific and reasonable, and can perform three-dimensional display and dynamic adjustment of the conditions before and after the construction of the earthwork project, with stronger visualization and flexible calculation. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] Figure 1 A flowchart of a method for analyzing earthwork based on comprehensive geophysical exploration as claimed in an embodiment of the present application;

[0037] Figure 2 A geophysical exploration result diagram of an earthwork analysis method based on comprehensive geophysical exploration as claimed in an embodiment of the present application;

[0038] Figure 3 A schematic diagram of geophysical detection lines and original terrain elevation points for an earthwork analysis method based on comprehensive geophysical exploration as claimed in an embodiment of the present application;

[0039] Figure 4 A schematic diagram of an original terrain TIN model of an earthwork analysis method based on comprehensive geophysical exploration claimed in an embodiment of the present application;

[0040] Figure 5 A schematic diagram of a TIN model of a cover boundary of an earthwork analysis method based on comprehensive geophysical exploration as claimed in an embodiment of the present application;

[0041] Figure 6 A schematic diagram of a weathering limit TIN model for an earthwork analysis method based on comprehensive geophysical exploration as claimed in an embodiment of the present application;

[0042] Figure 7 A schematic diagram of a design elevation TIN model of an earthwork analysis method based on comprehensive geophysical exploration as claimed in an embodiment of the present application;

[0043] Figure 8 A schematic diagram of an original terrain DEM model of an earthwork analysis method based on comprehensive geophysical exploration claimed in an embodiment of the present application;

[0044] Figure 9 A schematic diagram of a DEM model of the overburden boundary of an earthwork analysis method based on comprehensive geophysical exploration as claimed in an embodiment of the present application;

[0045] Figure 10 A schematic diagram of a weathering limit DEM model of an earthwork analysis method based on comprehensive geophysical exploration as claimed in an embodiment of the present application;

[0046] Figure 11 This is a schematic diagram of a design elevation DEM model for an earthwork analysis method based on comprehensive geophysical exploration, as claimed in an embodiment of the present application. DETAILED DESCRIPTION

[0047] The following will be combined with the drawings in the embodiments of this application to clearly and completely describe the technical solutions in the embodiments of this application. Obviously, the embodiments described are only part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of this application.

[0048] The terms "first", "second" and "third" in this application are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Therefore, a feature defined as "first", "second" and "third" may explicitly or implicitly include at least one of such features. In the description of this application, "multiple" means at least two, for example, two, three, etc., unless otherwise clearly and specifically defined. All directional indications in the embodiments of this application (such as up, down, left, right, front, back...) are only used to explain the relative positional relationship, movement, etc. between the components under a certain specific posture (as shown in the accompanying drawings). If the specific posture changes, the directional indication will also change accordingly. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or units is not limited to the listed steps or units, but may optionally also include steps or units that are not listed, or may optionally also include other steps or units inherent to these processes, methods, products or devices.

[0049] References to "embodiments" herein mean that a particular feature, structure, or characteristic described in connection with the embodiments may be included in at least one embodiment of the present application. The appearance of the phrase in various places in the specification does not necessarily refer to the same embodiment, nor does it constitute an independent or alternative embodiment that is mutually exclusive of other embodiments. It is understood, both explicitly and implicitly, by those skilled in the art that the embodiments described herein may be combined with other embodiments.

[0050] According to the first embodiment of the present invention, the present invention claims a method for analyzing earthwork based on comprehensive geophysical exploration, referring to Figure 1 ,include:

[0051] Obtain basic terrain data and use geophysical methods to obtain the ground elevation data of the airport to be analyzed;

[0052] Performing terrain design on the airport to be analyzed and obtaining cut-and-fill elevation data of the airport to be analyzed;

[0053] Import basic terrain data, stratum elevation data obtained from geophysical exploration, and cut-and-fill elevation data after terrain design into the ArcGIS platform, and generate irregular triangulated network (TIN) models for each.

[0054] Constructing a digital elevation model (DEM) based on the triangulated irregular network (TIN);

[0055] In combination with the shrinkage coefficient, the earthwork volume is calculated based on the area of each pixel in the digital elevation model (DEM) and the corresponding elevation change;

[0056] The construction is carried out based on the earthwork volume and according to the actual situation during the construction process, the digital elevation model DEM and the earthwork volume calculation result are dynamically adjusted.

[0057] Furthermore, the method of obtaining basic terrain data and obtaining stratum elevation data of the airport to be analyzed by using geophysical exploration methods further includes:

[0058] Arrange a physical detection line in a local area of the airport to be analyzed, and use seismic exploration and geological radar methods to obtain layered data of the weathering layer and the overburden layer;

[0059] Obtaining the resistivity distribution of the weathered layer and the cover layer, and combining the on-site geological analysis of the airport to be analyzed to obtain the soil composition, resistivity, and thickness of the cover layer and the weathered layer of the airport to be analyzed;

[0060] A first local terrain elevation point is generated based on the elevation of the measuring point of the physical detection line and the terrain data measured on site.

[0061] In this embodiment, a survey line is laid out in a local area of a mountain airport, and the geophysical prospecting results (inverse flux transient electromagnetic method) are as follows: Figure 2As shown. From the results diagram, it can be seen that the resistivity of the profile varies greatly, ranging from 10Ω·m to 500Ω·m, and is unevenly distributed along the survey line. In the vertical direction, the resistivity tends to gradually increase from shallow to deep. From the distribution of resistivity combined with the analysis of on-site geological conditions, the cover layer of this profile is mainly composed of shallow clay, silt and gravel, etc., with a resistivity generally less than 150Ω·m and a thickness of less than 5m. Local bedrock is exposed, and the overall thickness of the cover layer does not change much, and the distribution is relatively uniform. From the distribution of resistivity combined with the previous data and on-site geological conditions, below the bottom boundary of the cover layer, the resistivity is between 150Ω·m and 300Ω·m. It is speculated to be a strongly weathered layer of mud dolomite and mud limestone, with a thickness of generally 20m to 30m, with little change in thickness and a relatively uniform distribution. Reference Figure 3 , based on the elevation of the geophysical detection line measurement points and the terrain data measured on site, a local terrain elevation point 1 is generated.

[0062] Furthermore, the performing of terrain design on the airport to be analyzed and obtaining elevation data of cut and fill of the airport to be analyzed further includes:

[0063] Based on the buried depth of the cover layer and the severe weathering limit determined by the geophysical exploration results, the cover layer thickness and the severe weathering limit thickness of the area of the airport to be analyzed are determined by spatial interpolation, and the second elevation point of the cover layer and the third elevation point of the severe weathering limit layer are prepared;

[0064] Determine the excavation elevation point through terrain design.

[0065] In this embodiment, based on the buried depth of the overburden and the severe weathering limit determined by geophysical exploration results, spatial interpolation is used to determine the overburden thickness and the severe weathering limit thickness in the region. Elevation point 2 for the overburden and elevation point 3 for the severe weathering limit are generated. Furthermore, excavation elevation point 4 is determined through terrain design.

[0066] Furthermore, the step of importing basic terrain data, stratum elevation data obtained through geophysical exploration, and cut-and-fill elevation data after terrain design into the ArcGIS platform to generate triangulated irregular network (TIN) models for each of them further includes:

[0067] The original terrain TIN model is constructed using the first local terrain elevation point, the cover layer boundary TIN model is constructed using the second elevation point, the weathering boundary TIN model is constructed using the third elevation point, and the design elevation TIN model is determined using the terrain design cut elevation.

[0068] In this embodiment, the TIN model is constructed based on point and line data with elevation attributes. The TIN model is then generated using the Create TIN tool in ArcToolbox, using 3D Analyst Tools → Data Management → TIN. In the dialog box that pops up, set the path and name for the output TIN file, select the input feature class, and specify the elevation field as the height source. For example, if the input data is contour lines, select the contour elevation field; if it is point elevations, select the point elevation attribute field.

[0069] Reference Figure 4-7 , use the original terrain elevation point data 1 to construct the original terrain TIN model 1, use the elevation point data 2 to construct the cover layer boundary TIN model 2, use the elevation point data 3 to construct the weathering boundary TIN model 3, and use the terrain design excavation elevation to determine the design elevation TIN model 4.

[0070] Furthermore, the construction of a digital elevation model (DEM) based on the triangulated irregular network (TIN) further includes:

[0071] The original terrain TIN model is converted into an original terrain DEM model, the cover boundary TIN model is converted into a cover boundary DEM model, the weathering boundary TIN model is converted into a weathering boundary DEM model, and the design elevation TIN model is converted into a design elevation DEM model.

[0072] In this embodiment, when converting a TIN model to a DEM model, the TIN model is converted to DEM data in raster format through the [3D Analyst Tools] → [Conversion] → [From TIN] → [TIN to Raster] tool in ArcToolbox. In the dialog box, select the input TIN file, set the pixel size (i.e., resolution) of the output raster, and the path and name of the output file. The pixel size should be determined based on actual needs and data accuracy. A smaller pixel size can provide higher resolution, but it will increase the amount of data and calculation time.

[0073] Reference Figure 8-11 , the original terrain TIN model 1 is converted to the original terrain DEM model 1, the cover boundary TIN model 2 is converted to the cover boundary DEM model 2, the weathering boundary TIN model 3 is converted to the weathering boundary DEM model 3, and the design elevation TIN model 4 is converted to the design elevation DEM model 4.

[0074] Furthermore, the calculation of earthwork volume based on the area of each pixel in the digital elevation model (DEM) and the corresponding elevation change in combination with the shrinkage coefficient also includes:

[0075] The earthwork volume is calculated based on the following formula:

[0076] V = [cell (area)] × ΔZ;

[0077] Where: V is the calculated earthwork volume m 3 ; cell(area) is the area of a pixel to be calculated in DEM data m 2 ; ΔZ is the elevation change m corresponding to the calculated pixel before and after the proposed design.

[0078] Furthermore, the construction based on the earthwork volume and dynamically adjusting the digital elevation model DEM and earthwork volume calculation results according to the actual situation during the construction process also includes:

[0079] Based on the original terrain DEM model and the above-mentioned cover layer boundary DEM model, ArcGIS cut and fill analysis is used to calculate the cover layer cut volume;

[0080] Based on the cover boundary DEM model and the weathering boundary DEM model, ArcGIS cut and fill analysis is used to calculate the excavation volume of the strongly weathered rock layer;

[0081] Based on the weathering limit DEM model and the design elevation DEM model, the cut and fill volume of the medium-weathered rock layer is calculated using ArcGIS cut and fill analysis;

[0082] The volume of abandoned earth and filled earth of the airport to be analyzed is obtained based on the volume of excavation of the overburden layer, the volume of excavation of the strongly weathered rock layer and the volume of excavation of the moderately weathered rock layer, combined with the expansion coefficients of different rock layers.

[0083] In this embodiment, during the excavation process, the overburden layer generally needs to be cleared. In this calculation, the overburden layer clearing volume V1 = 488643.31m 3 .

[0084] According to the results of ArcGIS calculation, the excavation volume of strongly weathered dolomite and limestone V2 = 3049475.48m 3 According to the basic principle of earthwork balance, filling after excavation usually involves expansion. In this example, considering that the expansion coefficient of the strongly weathered layer is 1.1, if the excavated volume of strongly weathered dolomite and limestone is used for filling, the filling volume can reach 3049475.48×1.1=3354423.03m 3 .

[0085] According to the results of ArcGIS calculation, the excavation volume of medium-weathered dolomite and limestone V3 = 8523312.10m 3According to the basic principle of earthwork balance, filling after excavation usually involves expansion. In this example, considering that the expansion coefficient of the medium-weathered layer is 1.2, if the excavated volume of the strongly weathered dolomite and limestone is used for filling, the filling volume can reach 8523312.10×1.2=10227974.52m 3 .

[0086] Based on the above calculation results, it can be seen that the local area of the excavation area near an airport in the southwestern mountainous area involves abandoned earthwork V1 = 488643.31m 3 , the excavation volume is V2+V3=3049475.48+8523312.10=11572787.58m 3 Considering the expansion coefficient of different rock layers, the volume available for filling is (V2×1.1+V3×1.2)=3354423.03+10227974.52=13582397.55m 3 .

[0087] According to a second embodiment of the present invention, the present invention claims protection for an earthwork analysis system based on comprehensive geophysical exploration, comprising:

[0088] one or more processors;

[0089] A memory having one or more programs stored thereon, wherein when the one or more programs are executed by the one or more processors, the one or more processors implement the earthwork analysis method based on comprehensive geophysical exploration.

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

[0091] In addition, the functional units in the various embodiments of the present application can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The above-mentioned integrated units can be implemented in the form of hardware or in the form of software functional units. The above is only an implementation method of the present application and does not limit the patent scope of the present application. Any equivalent structure or equivalent process transformation made using the contents of the description and drawings of this application, or directly or indirectly used in other related technical fields, are also included in the patent protection scope of the present application.

[0092] The above detailed description of the specific embodiments of the invention is intended only as an example, and the present application is not limited to the specific embodiments described above. For those skilled in the art, any equivalent modifications or substitutions of the invention are also within the scope of the present application. Therefore, equivalent changes, modifications, and improvements made without departing from the spirit and scope of the present application should be included within the scope of the present application.

Claims

1. A soil and rock analysis method based on comprehensive geophysical exploration, characterized in that: include: Obtain basic terrain data and use geophysical methods to obtain the ground elevation data of the airport to be analyzed; Performing terrain design on the airport to be analyzed and obtaining cut-and-fill elevation data of the airport to be analyzed; Import basic terrain data, stratum elevation data obtained from geophysical exploration, and cut-and-fill elevation data after terrain design into the ArcGIS platform, and generate irregular triangulated network (TIN) models for each. Constructing a digital elevation model (DEM) based on the triangulated irregular network (TIN); In combination with the shrinkage coefficient, the earthwork volume is calculated based on the area of each pixel in the digital elevation model (DEM) and the corresponding elevation change; The construction is carried out based on the earthwork volume and according to the actual situation during the construction process, the digital elevation model DEM and the earthwork volume calculation result are dynamically adjusted.

2. The earthwork analysis method based on comprehensive geophysical exploration according to claim 1, characterized in that: The process of obtaining basic terrain data and obtaining stratum elevation data of the airport to be analyzed by using a geophysical prospecting method further includes: Arrange a physical detection line in a local area of the airport to be analyzed, and use seismic exploration and geological radar methods to obtain layered data of the weathering layer and the overburden layer; Obtaining the resistivity distribution of the weathered layer and the cover layer, and combining the on-site geological analysis of the airport to be analyzed to obtain the soil composition, resistivity, and thickness of the cover layer and the weathered layer of the airport to be analyzed; A first local terrain elevation point is generated based on the elevation of the measuring point of the physical detection line and the terrain data measured on site.

3. The earthwork analysis method based on comprehensive geophysical exploration according to claim 1, characterized in that: The performing of terrain design on the airport to be analyzed and obtaining elevation data of cut and fill of the airport to be analyzed further includes: Based on the buried depth of the cover layer and the severe weathering limit determined by the geophysical exploration results, the cover layer thickness and the severe weathering limit thickness of the area of the airport to be analyzed are determined by spatial interpolation, and the second elevation point of the cover layer and the third elevation point of the severe weathering limit layer are prepared; Determine the excavation elevation point through terrain design.

4. The earthwork analysis method based on comprehensive geophysical exploration according to claim 1 is characterized in that: The method of importing basic terrain data, stratum elevation data obtained by geophysical exploration, and cut-and-fill elevation data after terrain design into the ArcGIS platform to generate triangulated irregular network (TIN) models respectively includes: The original terrain TIN model is constructed using the first local terrain elevation point, the cover layer boundary TIN model is constructed using the second elevation point, the weathering boundary TIN model is constructed using the third elevation point, and the design elevation TIN model is determined using the terrain design cut elevation.

5. The earthwork analysis method based on comprehensive geophysical exploration according to claim 4 is characterized in that: The method of constructing a digital elevation model (DEM) based on the triangulated irregular network (TIN) further comprises: The original terrain TIN model is converted into an original terrain DEM model, the cover boundary TIN model is converted into a cover boundary DEM model, the weathering boundary TIN model is converted into a weathering boundary DEM model, and the design elevation TIN model is converted into a design elevation DEM model.

6. The earthwork analysis method based on comprehensive geophysical exploration according to claim 4 is characterized in that: The method of calculating the earthwork volume based on the area of each pixel in the digital elevation model (DEM) and the corresponding elevation change in combination with the shrinkage coefficient also includes: The earthwork volume is calculated based on the following formula: V = [cell (area)] × ΔZ; Where: V is the calculated earthwork volume m 3 ; cell(area) is the area of a pixel to be calculated in DEM data m 2 ; ΔZ is the elevation change m corresponding to the calculated pixel before and after the proposed design.

7. The earthwork analysis method based on comprehensive geophysical exploration according to claim 5 is characterized in that: The method of performing construction based on the earthwork volume and dynamically adjusting the digital elevation model (DEM) and the earthwork volume calculation result according to the actual situation during the construction process further includes: Based on the original terrain DEM model and the above-mentioned cover layer boundary DEM model, ArcGIS cut and fill analysis is used to calculate the cover layer cut volume; Based on the cover boundary DEM model and the weathering boundary DEM model, ArcGIS cut and fill analysis is used to calculate the excavation volume of the strongly weathered rock layer; Based on the weathering limit DEM model and the design elevation DEM model, the cut and fill volume of the medium-weathered rock layer is calculated using ArcGIS cut and fill analysis; The volume of abandoned earth and filled earth of the airport to be analyzed is obtained based on the volume of excavation of the overburden layer, the volume of excavation of the strongly weathered rock layer and the volume of excavation of the moderately weathered rock layer, combined with the expansion coefficients of different rock layers.

8. An earthwork analysis system based on comprehensive geophysical exploration, characterized in that: include: one or more processors; A memory having one or more programs stored thereon, wherein when the one or more programs are executed by the one or more processors, the one or more processors implement a soil and rock analysis method based on integrated geophysical exploration according to any one of claims 1 to 7.

Images