Hydraulic engineering geological data acquisition method and system

Abstract

The invention relates to the technical field of electric digital data processing, in particular to a hydraulic engineering geological data acquisition method and system. The method comprises the steps of determining a plurality of target points according to surveying and mapping reference points, obtaining various types of data of the target points, and then obtaining interference degrees of the target points according to fluctuation of geological data; constructing a feature space according to the geological attributes in the geological data of the target point, and obtaining an abnormal factor of the target point according to the distribution of data points in the feature space and the interference degree of the target point; obtaining the abnormal degree of the target point according to the terrain complexity and the abnormal factor of the target point; according to the abnormal degree of the target point and the actual distance between the target point and the to-be-interpolated point, obtaining the distance between the target point and the to-be-interpolated point after weighting adjustment; and finally, performing interpolation processing on the distance between the target point and the to-be-interpolated point after weighting adjustment to obtain data of all geological points. According to the invention, the accuracy of geological data acquisition can be improved.

Description

Technical Field

[0001] The invention relates to the technical field of electric digital data processing, and in particular to a method and system for collecting geological data of water conservancy projects. Background Art

[0002] The collection of geological data for water conservancy projects is a key activity in the construction process of water conservancy projects. It aims to provide basic data and scientific basis for engineering design, construction and management through detailed investigation and analysis of the geological environment. The collection of geological data for water conservancy projects plays a key role in risk assessment, design optimization, construction guidance and other aspects of water conservancy projects.

[0003] When surveying and mapping the geology in water conservancy projects, geological profiles are obtained by collecting various types of information from geological points. However, when obtaining various information data from geological points, there are great differences between the different soil compositions in different places and the corresponding terrains, which makes it easy to be interfered by noise when sampling some soils or terrains, making the obtained geological point information data not accurate. This leads to errors when obtaining all geological point information data through inverse distance weighted interpolation of some geological points, thereby reducing the accuracy of geological data collection. Summary of the invention

[0004] In order to solve the technical problem of low accuracy of geological data collection in the prior art, the purpose of the present invention is to provide a method and system for collecting geological data of water conservancy projects. The technical scheme adopted is as follows: In a first aspect, a method for collecting water conservancy engineering geological data is provided, the method comprising: Step S1: determining multiple target points according to the surveying and mapping reference points, and obtaining various types of data of the target points; Step S2: obtaining the interference degree of the target point according to the fluctuation of the various types of data; Step S3: constructing a feature space based on the geological attributes in various types of data of the target point, and obtaining an abnormal factor of the target point according to the distribution of data points in the feature space and the interference degree of the target point; Step S4: Obtain the abnormality degree of the target point according to the terrain complexity and abnormality factor of the target point; Step S5: according to the abnormality degree of the target point, a weighted distance between the target point and the point to be interpolated is obtained, and various types of data of the point to be interpolated are obtained by interpolation according to the weighted distance.

[0005] Furthermore, the step S1 specifically includes: Determine the surveying and mapping benchmarks, and determine multiple target points with a preset interval distance threshold by triangulation or grid method; Use a total station to obtain the coordinates, height, angle and distance of the target point, and use geological exploration to obtain the geological stratum, geological properties, meteorological data and groundwater level data of the target point; among them, the geological properties include the density, porosity and permeability of the soil or rock.

[0006] Furthermore, the step S2 specifically includes: According to a type of data of the target point, the mean of the type of data of all neighboring target points corresponding to the target point, and the standard deviation of the type of data of all neighboring target points corresponding to the target point, the disturbance degree of the type of data is obtained; According to the disturbance degree of all types of data of the target point, the initial disturbance degree of the target point is obtained; The initial interference degree of the target point is corrected by the distribution of the target point's neighborhood target points in the feature space to obtain the interference degree of the target point.

[0007] Furthermore, the initial interference degree of the target point is corrected by the distribution of the target point's neighborhood target points in the feature space, and the interference degree of the target point is obtained as follows: The product of the normalized data of the porosity, permeability and thickness of the geological layer of the rock or soil at each target point is recorded as the high abnormal dimension feature of each target point, and the product of the normalized data of the density of each target point is recorded as the low abnormal dimension feature of each target point; The feature space of the target point is constructed with the high abnormal dimension feature of the target point as the vertical axis and the low abnormal dimension feature as the horizontal axis, and all the target points are mapped in the feature space; The interference degree of the target point is obtained according to the distance between a target point and a corresponding neighborhood target point in the feature space, the number of neighborhood target points, and the initial interference degree of the target point.

[0008] Furthermore, the distance between the target point and a corresponding neighborhood target point in the feature space and the initial interference degree of the target point are both positively correlated with the interference degree of the target point.

[0009] Furthermore, the step S3 specifically includes: The abnormality factor of a target point is obtained according to the interference degree of the target point, the difference between the high abnormality dimension feature of the target point and the minimum high abnormality dimension feature, and the difference between the low abnormality dimension feature of the target point and the maximum low abnormality dimension feature.

[0010] Furthermore, the interference degree of a target point, the difference between the high abnormal dimension feature of the target point and the minimum high abnormal dimension feature, and the difference between the low abnormal dimension feature of the target point and the maximum low abnormal dimension feature are all positively correlated with the abnormal factor of the target point.

[0011] Furthermore, the step S4 specifically includes: According to the cosine value of the angle between a target point and the inclination vector of a corresponding neighboring target point and the number of neighboring target points, the terrain complexity of the target point is obtained; According to the terrain complexity and abnormal factor of each target point, the abnormal degree of each target point is obtained.

[0012] Furthermore, step S5 specifically includes: According to the abnormality degree of the target point and the actual distance between the target point and the point to be interpolated, a weighted adjusted distance between the target point and the point to be interpolated is obtained; According to the weighted distance between the target point and the point to be interpolated, the inverse distance weighted interpolation method is used to interpolate and obtain various types of data of the interpolation point.

[0013] In another aspect, the present invention provides a water conservancy engineering geological data acquisition system, the system comprising: A data acquisition module is used to determine multiple target points based on surveying and mapping benchmark points and acquire various types of data of the target points; A disturbance degree acquisition module, used for acquiring the disturbance degree of a target point according to the fluctuation of the geological data; The abnormal factor acquisition module is used to construct a feature space based on the geological attributes in various types of data of the target point, and obtain the abnormal factor of the target point according to the distribution of data points in the feature space and the interference degree of the target point; The abnormality degree acquisition module is used to obtain the abnormality degree of the target point according to the terrain complexity and abnormality factor of the target point; A distance acquisition module is used to acquire a weighted distance between a target point and a point to be interpolated according to the abnormality degree of the target point and the actual distance between the target point and the point to be interpolated; The interpolation point data acquisition module is used to perform interpolation processing on the weighted adjusted distance between the target point and the point to be interpolated, and obtain various types of data of the interpolation point; The interference degree acquisition module specifically includes: obtaining the disturbance degree of a type of data of the target point, the mean of the type of data of all neighboring target points corresponding to the target point, and the standard deviation of the type of data of all neighboring target points corresponding to the target point; obtaining the initial interference degree of the target point according to the disturbance degree of all types of data of the target point; and correcting the initial interference degree of the target point by the distribution of the neighboring target points of the target point in the feature space to obtain the interference degree of the target point.

[0014] The present invention has the following beneficial effects: obtaining the initial interference degree of each target point through data fluctuations; constructing a feature space of geological attributes of the target point, obtaining the abnormal factor of each target point through the distribution of data points in the feature space and the initial interference degree, reducing the influence of geological stability, and improving the accuracy of geological stability analysis; obtaining the abnormal degree of each target point according to the abnormal factor of each target point and the degree of consistency of the tendency of the geological target point, reducing the interference of complex terrain on data collection, and improving the accuracy of terrain complexity analysis; weighting the distance between the target point and the point to be interpolated by the abnormal degree and interpolating, completing the data collection of all geological points, and effectively improving the accuracy of geological data collection. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] In order to more clearly illustrate the technical solutions and advantages in the embodiments of the present invention or the prior art, the drawings required for use in the embodiments or the prior art descriptions are briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0016] Figure 1 A flow chart of a method for collecting geological data of water conservancy projects provided by one embodiment of the present invention; Figure 2 A block diagram of a water conservancy engineering geological data acquisition system provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0017] In order to further explain the technical means and effects adopted by the present invention to achieve the predetermined invention purpose, the following is a detailed description of a method and system for collecting geological data of water conservancy engineering proposed by the present invention, its specific implementation method, structure, characteristics and effects, in combination with the accompanying drawings and preferred embodiments. In the following description, different "one embodiment" or "another embodiment" does not necessarily refer to the same embodiment. Unless otherwise defined, all technical and scientific terms used in this article have the same meaning as those commonly understood by technicians in the technical field of the present invention.

[0018] The scenario targeted by the present invention is that in the process of collecting geological data for water conservancy projects, due to the complexity of the terrain and the different stability of the geology in different regions, data interference analysis is performed through terrain characteristics.

[0019] The specific scheme of a water conservancy engineering geological data acquisition method and system provided by the present invention is described in detail below with reference to the accompanying drawings.

[0020] First, see Figure 1, which shows a flow chart of a method for collecting geological data of water conservancy projects provided by an embodiment of the present invention, the method comprising the following steps: Step S1: Determine multiple target points based on the surveying and mapping reference points, and obtain various types of data of the target points.

[0021] Among them, step S1 specifically includes: determining the surveying and mapping benchmark point, determining multiple target points by using triangulation or grid method with a preset interval distance threshold, using a total station to obtain the coordinates, height, angle, distance and other types of data of the target points, and obtaining the geological layer, geological attributes, meteorological data and groundwater level data of the target points through geological exploration. Among them, the geological attributes include density, porosity and permeability of soil or rock.

[0022] More specifically, first determine the surveying and mapping benchmark points, obtain all coordinate data and geological data of each benchmark point, and then evenly determine multiple target points with a preset interval distance threshold. For example, the preset distance threshold can be 100 meters, that is, the interval between target points is 100 meters. Then, use the total station to use the coordinates of the benchmark points as a reference to use triangulation or grid method to obtain various types of data such as coordinates, heights, angles and distances of all target points; then explore and obtain the geological strata (thickness and arrangement order of each geological stratum), geological properties (density, porosity, permeability), meteorological data and groundwater level data of each benchmark point and target point. Among them, the geological strata and geological properties are obtained by drilling and sampling at the target point and sending them to the laboratory for measurement and analysis, and the meteorological data and groundwater level data are collected and obtained using corresponding sensors.

[0023] When collecting various types of data at the target point, information such as the water content and temperature in the soil affects the propagation and reflection characteristics of electromagnetic waves, and the differences in the geological properties of the soil and the complexity of the terrain make the stability of the soil different. Therefore, it is necessary to assist in analyzing and judging the degree of interference when collecting data at different target points based on geological data such as geological layers, geological properties, meteorological data, and groundwater level data. Based on this, this embodiment further sets the following steps.

[0024] Step S2: Obtain the interference degree of the target point according to the fluctuations of the various types of data.

[0025] Among them, step S2 specifically includes: according to a type of data of the target point, the mean of the type of data of all neighboring target points corresponding to the target point, and the standard deviation of the type of data of all neighboring target points corresponding to the target point, obtaining the disturbance degree of the data of this type; according to the disturbance degree of all types of data of the target point, obtaining the initial interference degree of the target point; correcting the initial interference degree of the target point by the distribution of the neighboring target points of the target point in the feature space, and obtaining the interference degree of the target point.

[0026] When obtaining the initial interference degree of the target point, the moisture content, temperature and other information in the general soil affect the propagation and reflection characteristics of electromagnetic waves, resulting in large fluctuations in the various data types of the collected geological points. Therefore, the degree of interference is analyzed by the difference between the data of each geological point and the neighboring geological points.

[0027] More specifically, first, obtain all neighboring target points (including the target point itself) within a radius of 1000 meters for each target point; calculate the mean of each type of data corresponding to all neighboring target points and standard deviation ; Preset a fluctuation threshold ; When the value of each type of data at each target point is within the fluctuation range, it is considered that there is no abnormality, and the interference degree is considered to be 0; when the value of each type of data at each target point is not within the fluctuation range, it is considered that there is an abnormality.

[0028] The mathematical calculation formula for the disturbance degree of each type of data at each target point constructed in this embodiment is as follows: ; In the formula, Represents the first Types of data; Indicates that each target point corresponds to the first The mean of the data types; Indicates that each target point corresponds to the first The standard deviation of the data of each type; since the data of all neighboring target points corresponding to each target point cannot be exactly the same, obviously, ; Represents the first The degree of disturbance of the data type.

[0029] In the mathematical calculation formula for the disturbance degree of each type of data at each target point constructed above, It represents the process of dividing the difference between a number and the mean by the standard deviation, that is, the calculation of the standard score, and the result can reflect the first Types of data The relative standard distance of the distance average can reflect the The degree of disturbance of the data of each type; when the value of each type of data at each target point is If the value of each type of data at each target point is within the range of If it is within the range, it is considered to be abnormal.

[0030] Then, according to the disturbance degree of all types of data of the target point, the initial disturbance degree of the target point is obtained. The mathematical calculation formula of the initial disturbance degree of each target point constructed in this embodiment is as follows: ; In the formula, Indicates the number of all types of data. represents the initial interference degree of each target point, Represents the normalization function.

[0031] In the mathematical calculation formula for the initial interference degree of each target point constructed above, It represents the product of the normalized disturbance degree of all types of data. The greater the fluctuation difference between each type of data of each target point and the mean of each type of data of the corresponding neighborhood target point, the greater the interference. The greater the interference of all types of data, the greater the initial interference degree of the target point.

[0032] When obtaining the initial interference degree of the target point, because of the differences between different soil types, the stability of the physical state affecting the soil is different, and the interference during data collection is different. For example, when the high-density soil type means that the soil is more compact, the electromagnetic wave propagation is relatively stable, and the interference degree may be low, while the low-density soil may cause more pores and air, and the electromagnetic wave propagation is unstable, increasing the interference; the soil with high porosity has more water and air, and the electromagnetic wave may be more scattered and attenuated, thereby increasing the data interference; the soil with high permeability is easy for water to flow, and the water content changes frequently, which may cause the data fluctuation to increase during measurement, while the stability of the soil corresponding to low porosity and low permeability is higher, that is, the interference degree is smaller; the thick layer of soil causes multiple reflections and refraction phenomena, which increases the uncertainty and interference degree of the measurement, and the thin layer of soil has fewer layer changes, the electromagnetic wave propagation is relatively simple, and the interference degree is low. Therefore, according to the characteristics of the collected soil type, the initial interference degree is corrected to obtain the interference degree of each target point after correction.

[0033] Among them, the initial interference degree of the target point is corrected by the distribution of the target point's neighborhood target points in the feature space, and the interference degree of the target point is obtained as follows: the product of the standardized data of the porosity, permeability and thickness of the geological layer of each target point rock or soil is recorded as the high anomaly dimensional feature of each target point, and the product of the standardized data of the density of each target point is recorded as the low anomaly dimensional feature of each target point; the feature space of the target point is constructed with the high anomaly dimensional feature of the target point as the vertical axis and the low anomaly dimensional feature as the horizontal axis, and all target points are mapped in the feature space; the interference degree of the target point is obtained according to the distance between a target point and a corresponding neighborhood target point in the feature space, the number of neighborhood target points, and the initial interference degree of the target point.

[0034] More specifically, first, the density, porosity, permeability and thickness of the geological layer are standardized respectively; the product of the standardized data of the porosity, permeability and thickness of the geological layer of each target point is recorded as the high abnormal dimension feature of each target point, and the product of the standardized data of the density of each target point is recorded as the low abnormal dimension feature of each target point. Then, the feature space of the target point is constructed by taking the high abnormal dimension feature of the target point as the vertical axis and the low abnormal dimension feature as the horizontal axis; all target points are mapped in the feature space. Since the closer the geological characteristic attributes corresponding to the neighboring target points of each target point are, the smaller the distance corresponding to the neighboring target points of each target point in the feature space is, which means that the uniformity and stability of the geological environment corresponding to the target point are better, that is, the less interference is received. Therefore, the initial interference degree of each target point is corrected by the distribution of the neighboring target points of each target point in the feature space, and the corrected interference degree of each target point is obtained.

[0035] In this embodiment, the mathematical calculation formula for the interference degree of the target point after correction is constructed is as follows; ; In the formula, Indicates The target point and the corresponding The distance between neighboring target points in the feature space; Represents the number of neighborhood target points. Since neighborhood target points must exist, ; Indicates The initial interference degree of the target point is Indicates The interference degree of each target point after correction.

[0036] In the mathematical calculation formula for the interference degree of the target point after correction constructed above, Indicates traversing all neighborhood target points, The average distance between each target point and the corresponding neighboring target point in the feature space. The closer the distance between each target point and all neighboring target points in the feature space, the better the uniformity and stability of the geological environment around the target point, that is, the less interference, so its value is positively correlated with the first The interference degree of the target point after correction; obviously, The initial interference degree of the target point is also positively correlated with the The interference degree of each target point after correction.

[0037] Therefore, the distance between a target point and a corresponding neighborhood target point in the feature space and the initial interference degree of the target point are both positively correlated with the interference degree of the target point.

[0038] The data anomaly of the target point is not limited to the interference degree of the target point. In the geological feature space, the higher the value of the high anomaly dimension feature, the more abnormal it is, while the lower the value of the low anomaly dimension feature, the more abnormal it is. Therefore, it is necessary to adjust the interference degree of the target point by the height of the two feature values ​​corresponding to each target point. Therefore, this embodiment further sets the following steps.

[0039] Step S3: construct a feature space based on the geological attributes in various types of data of the target point, and obtain the abnormal factor of the target point according to the distribution of data points in the feature space and the interference degree of the target point.

[0040] Among them, step S3 specifically includes: obtaining the abnormal factor of the target point according to the interference degree of a target point, the difference between the high abnormal dimensional feature of the target point and the minimum high abnormal dimensional feature, and the difference between the low abnormal dimensional feature of the target point and the maximum low abnormal dimensional feature.

[0041] In this embodiment, the mathematical calculation formula of the abnormal factor of the target point is constructed as follows: ; In the formula, Indicates The abnormal factor of each target point; Indicates The interference degree of each target point after correction; , represents the weight coefficient; Indicates Highly abnormal dimensional features of target points; Indicates Low abnormal dimensional features of target points; Represents the minimum value of the high abnormal dimension features of all target points; Indicates the maximum value of the high abnormal dimension features of all target points; Represents the maximum value among the low-anomaly dimensional features of all target points; Represents the minimum value among the low-anomaly dimensional features of all target points.

[0042] In the mathematical calculation formula of the abnormal factor of the target point constructed above, The interference degree of each target point after correction Positively correlated with The abnormal factor of each target point; It represents the difference between the maximum and minimum values ​​in the high anomaly dimensional features of all target points. Its value is obviously not 0. The difference between each target point and the smallest high abnormal dimension feature is normalized again. The larger the high abnormal feature of the target point, the larger the result, which means the more abnormal it is. It represents the difference between the maximum and minimum values ​​in the low-anomaly dimensional features of all target points. Its value is obviously not 0. The difference between each target point and the largest low-abnormal dimensional feature is normalized. The smaller the low-abnormal feature of the target point, the larger the result, which means the more abnormal it is. The high-abnormal dimensional feature of the target point is the product of the standardized data of the porosity, permeability and thickness of the geological layer of each target point, which has three features, and the low-abnormal dimensional feature of the target point is the product of the standardized data of the density of each target point, which has one feature. Therefore, the weight between the two is 3 to 1, so the definition is , .

[0043] Therefore, the degree of interference of a target point, the difference between the high abnormal dimension feature of the target point and the smallest high abnormal dimension feature, and the difference between the low abnormal dimension feature of the target point and the largest low abnormal dimension feature are all positively correlated with the abnormal factor of the target point.

[0044] When faults, folds, rock layer contacts and other terrains appear in the terrain, the terrain is relatively complex; and because the height difference between such terrains is large, the consistency of the direction and tendency of the target point is poorer. Therefore, the degree of abnormality of each target point can be obtained through the degree of consistency of the direction and the abnormality of the target point. Accordingly, the present embodiment further sets the following steps.

[0045] Step S4: Obtain the abnormality degree of the target point according to the terrain complexity and abnormality factor of the target point.

[0046] Among them, step S4 specifically includes: obtaining the terrain complexity of the target point according to the cosine value of the angle between a target point and the inclination vector of a corresponding neighboring target point and the number of neighboring target points; obtaining the abnormality degree of each target point according to the terrain complexity and abnormality factor of each target point.

[0047] More specifically, first, the height vector between each target point and each neighboring target point is obtained. The magnitude of the vector is the height difference between the two points, and the direction is from the point with a higher altitude to the point with a lower altitude. The vector sum operation is performed on the vectors between each target point and all neighboring target points to obtain the inclination vector of each target point. Similarly, the inclination vectors of the neighboring target points of the target point are obtained. Then, the terrain complexity of each target point is obtained by the difference between the inclination vectors of each target point and the corresponding neighboring target point.

[0048] In this embodiment, the mathematical calculation formula for constructing the terrain complexity of each target point is as follows: ; In the formula, Indicates The target point and the corresponding The cosine value of the angle between the inclination vectors of the neighboring target points; Represents the number of neighborhood target points. Since neighborhood target points must exist, ; Indicates The complexity of the terrain at each target point; Represents an exponential function with the natural constant e as the base.

[0049] In the mathematical calculation formula for the terrain complexity of each target point constructed above, The target point and the corresponding The cosine value of the angle between the inclination vectors of the neighboring target points , when the cosine value is larger, the angle is smaller, which means the terrain near the target point is less complex, the terrain complexity factor of the target point is smaller, and the data of the target point is less affected by the terrain complexity; in addition, the exponential function with natural constant as the base is convenient to limit the result to between.

[0050] Furthermore, in this embodiment, the mathematical calculation formula for the abnormality degree of each target point is constructed as follows: ; In the formula, Indicates The degree of abnormality of each target point; Indicates The complexity of the terrain at each target point; Indicates The abnormal factor of each target point; represents the linear normalization function.

[0051] In the mathematical calculation formula for the abnormality degree of each target point constructed above, The complexity of the terrain at each target point , No. The abnormal factor of the target point Both are positively correlated with The abnormality of the target point , when the terrain complexity of the target point is greater and the anomaly factor of the target point is greater, the anomaly degree of the target point is greater.

[0052] Step S5: according to the abnormality degree of the target point, a weighted distance between the target point and the point to be interpolated is obtained, and various types of data of the point to be interpolated are obtained by interpolation according to the weighted distance.

[0053] Among them, step S5 specifically includes: obtaining the weighted adjusted distance between the target point and the point to be interpolated according to the abnormality degree of the target point and the actual distance between the target point and the point to be interpolated; interpolating using the inverse distance weighted interpolation method according to the weighted adjusted distance between the target point and the point to be interpolated to obtain various types of data of the interpolation point.

[0054] The greater the abnormality of each target point, the lower the credibility of the data passing through the point, which means that the target point is farther away from the point to be interpolated. Therefore, the actual distance between the target point and the point to be interpolated can be weighted to ensure the accuracy of the difference data of the point to be interpolated.

[0055] In this embodiment, the mathematical calculation formula for constructing the weighted adjusted distance between the target point and the point to be interpolated is as follows: ; In the formula, Indicates The degree of abnormality of each target point; Represents the actual distance between the target point and the point to be interpolated. Represents the weighted distance between the target point and the point to be interpolated.

[0056] In the weighted distance mathematical calculation formula between the target point and the point to be interpolated constructed above, the actual distance between the target point and the point to be interpolated is As a benchmark value, Represents the weight. The greater the abnormality of the target point, the farther the target point is from the interpolation point. Therefore, the actual distance between the target point and the interpolation point is The basis needs to be increased the farther the distance.

[0057] Furthermore, according to the weighted distance between the target point and the point to be interpolated, the inverse distance weighted interpolation method is used to interpolate to obtain various types of data of the interpolation point; based on the various types of data of the interpolation point and the target point, the data of all geological points are obtained.

[0058] In the second aspect, this embodiment provides a water conservancy engineering geological data acquisition system, see Figure 2 , which shows a block diagram of a water conservancy engineering geological data acquisition system provided by an embodiment of the present invention, the system includes: The data acquisition module 101 is used to determine multiple target points according to the surveying and mapping reference points and acquire various types of data of the target points; The interference degree acquisition module 102 is used to acquire the interference degree of the target point according to the fluctuation of the geological data; The abnormal factor acquisition module 103 is used to construct a feature space based on the geological attributes in the geological data of the target point, and acquire the abnormal factor of the target point according to the distribution of the data points in the feature space and the interference degree of the target point; The abnormality degree acquisition module 104 is used to obtain the abnormality degree of the target point according to the terrain complexity and abnormality factor of the target point; The distance acquisition module 105 is used to acquire the weighted distance between the target point and the point to be interpolated according to the abnormality degree of the target point and the actual distance between the target point and the point to be interpolated; The interpolation point data acquisition module 106 is used to perform interpolation processing on the weighted adjusted distance between the target point and the to-be-interpolated point to acquire various types of data of the interpolation point.

[0059] The interference degree acquisition module 102 specifically includes: obtaining the disturbance degree of a type of data of the target point, the mean of the type of data of all neighboring target points corresponding to the target point, and the standard deviation of the type of data of all neighboring target points corresponding to the target point; obtaining the initial interference degree of the target point according to the disturbance degree of all types of data of the target point; and correcting the initial interference degree of the target point by the distribution of the neighboring target points of the target point in the feature space to obtain the interference degree of the target point.

[0060] Furthermore, the system also includes a terrain complexity acquisition module for acquiring the terrain complexity of a target point based on the cosine value of the angle between a target point and a corresponding inclination vector of a neighboring target point and the number of neighboring target points.

[0061] It should be further explained that the interference degree acquisition module 102 is used to obtain the interference degree of the target point according to the fluctuation of the geological data. Specifically, the following steps are taken: according to a type of data of the target point, the mean of the type of data of all neighboring target points corresponding to the target point, and the standard deviation of the type of data of all neighboring target points corresponding to the target point, the disturbance degree of the type of data is obtained; according to the disturbance degree of all types of data of the target point, the initial interference degree of the target point is obtained; the product of the porosity, permeability and thickness of the geological layer of the rock or soil of each target point is recorded as the high abnormal dimensional feature of each target point, and the product of the density standardized data of each target point is recorded as the low abnormal dimensional feature of each target point; the feature space of the target point is constructed with the high abnormal dimensional feature of the target point as the vertical axis and the low abnormal dimensional feature as the horizontal axis, and all target points are mapped in the feature space; the interference degree of the target point is obtained according to the distance between a target point and a corresponding neighboring target point in the feature space, the number of neighboring target points, and the initial interference degree of the target point.

[0062] The present embodiment provides a method and system for collecting geological data of water conservancy projects. Multiple target points are determined according to surveying and mapping benchmark points to obtain various types of data of the target points; the interference degree of the target points is obtained according to the fluctuation of the geological data; the feature space is constructed based on the geological attributes in the geological data of the target points, and the abnormal factors of the target points are obtained according to the distribution of data points in the feature space and the interference degree of the target points; the abnormal degree of the target points is obtained according to the terrain complexity and abnormal factors of the target points; the weighted adjusted distance between the target points and the points to be interpolated is obtained according to the abnormal degree of the target points and the actual distance between the target points and the points to be interpolated; the weighted adjusted distance between the target points and the points to be interpolated is interpolated to obtain the data of all geological points. This effectively improves the accuracy of geological data collection.

[0063] It should be noted that the sequence of the above embodiments of the present invention is only for description and does not represent the advantages and disadvantages of the embodiments. The processes depicted in the accompanying drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0064] The embodiments described above are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the aforementioned embodiments, a person skilled in the art should understand that the technical solutions described in the aforementioned embodiments may still be modified, or some of the technical features may be replaced by equivalents. Such modifications or replacements do not deviate the essence of the corresponding technical solutions from the scope of the technical solutions of the embodiments of the present application, and should all be included in the protection scope of the present application.

Claims

1. A method for collecting geological data of water conservancy projects, characterized in that: The method comprises: Step S1: determining multiple target points according to the surveying and mapping reference points, and obtaining various types of data of the target points; Step S2: obtaining the interference degree of the target point according to the fluctuation of the various types of data; Step S3: constructing a feature space based on the geological attributes in various types of data of the target point, and obtaining the abnormal factor of the target point according to the distribution of the data points in the feature space and the interference degree of the target point; Step S4: Obtain the abnormality degree of the target point according to the terrain complexity and abnormality factor of the target point; Step S5: according to the abnormality degree of the target point, obtaining the weighted adjusted distance between the target point and the point to be interpolated, and interpolating and obtaining various types of data of the point to be interpolated according to the weighted adjusted distance; The step S2 specifically includes: obtaining the disturbance degree of a type of data of the target point, the mean of the type of data of all neighboring target points corresponding to the target point, and the standard deviation of the type of data of all neighboring target points corresponding to the target point; obtaining the initial disturbance degree of the target point according to the disturbance degree of all types of data of the target point; and correcting the initial disturbance degree of the target point by the distribution of the neighboring target points of the target point in the feature space to obtain the disturbance degree of the target point.

2. A method for collecting water conservancy engineering geological data according to claim 1, characterized in that: The step S1 specifically includes: Determine the surveying and mapping benchmarks, and determine multiple target points with a preset interval distance threshold by triangulation or grid method; Use a total station to obtain the coordinates, height, angle and distance of the target point, and use geological exploration to obtain the geological stratum, geological properties, meteorological data and groundwater level data of the target point; among them, the geological properties include the density, porosity and permeability of the soil or rock.

3. A method for collecting water conservancy engineering geological data according to claim 1, characterized in that: The initial interference degree of the target point is corrected by the distribution of the target point's neighborhood target points in the feature space, and the interference degree of the target point is obtained as follows: The product of the normalized data of the porosity, permeability and thickness of the geological layer of the rock or soil at each target point is recorded as the high abnormal dimension feature of each target point, and the product of the normalized data of the density of each target point is recorded as the low abnormal dimension feature of each target point; The feature space of the target point is constructed with the high abnormal dimension feature of the target point as the vertical axis and the low abnormal dimension feature as the horizontal axis, and all the target points are mapped in the feature space; The interference degree of the target point is obtained according to the distance between a target point and a corresponding neighborhood target point in the feature space, the number of neighborhood target points, and the initial interference degree of the target point.

4. A method for collecting water conservancy engineering geological data according to claim 3, characterized in that: The distance between the target point and a corresponding neighborhood target point in the feature space and the initial interference degree of the target point are both positively correlated to the interference degree of the target point.

5. A method for collecting water conservancy engineering geological data according to claim 1, characterized in that: The step S3 specifically includes: The abnormality factor of a target point is obtained according to the interference degree of the target point, the difference between the high abnormality dimension feature of the target point and the minimum high abnormality dimension feature, and the difference between the low abnormality dimension feature of the target point and the maximum low abnormality dimension feature.

6. A method for collecting water conservancy engineering geological data according to claim 5, characterized in that: The degree of interference of a target point, the difference between the high abnormal dimension feature of the target point and the smallest high abnormal dimension feature, and the difference between the low abnormal dimension feature of the target point and the largest low abnormal dimension feature are all positively correlated with the abnormal factor of the target point.

7. A method for collecting water conservancy engineering geological data according to claim 1, characterized in that: The step S4 specifically includes: According to the cosine value of the angle between a target point and the inclination vector of a corresponding neighboring target point and the number of neighboring target points, the terrain complexity of the target point is obtained; According to the terrain complexity and abnormal factor of each target point, the abnormal degree of each target point is obtained.

8. A method for collecting water conservancy engineering geological data according to claim 1, characterized in that: The step S5 specifically includes: According to the abnormality degree of the target point and the actual distance between the target point and the point to be interpolated, a weighted adjusted distance between the target point and the point to be interpolated is obtained; According to the weighted distance between the target point and the point to be interpolated, the inverse distance weighted interpolation method is used to interpolate and obtain various types of data of the interpolation point.

9. A water conservancy engineering geological data acquisition system, characterized in that: The system comprises: A data acquisition module is used to determine multiple target points based on surveying and mapping benchmark points and acquire various types of data of the target points; A disturbance degree acquisition module, used for acquiring the disturbance degree of a target point according to the fluctuation of the geological data; The abnormal factor acquisition module is used to construct a feature space based on the geological attributes in various types of data of the target point, and obtain the abnormal factor of the target point according to the distribution of data points in the feature space and the interference degree of the target point; The abnormality degree acquisition module is used to obtain the abnormality degree of the target point according to the terrain complexity and abnormality factor of the target point; A distance acquisition module is used to acquire a weighted distance between a target point and a point to be interpolated according to the abnormality degree of the target point and the actual distance between the target point and the point to be interpolated; The interpolation point data acquisition module is used to perform interpolation processing on the weighted adjusted distance between the target point and the point to be interpolated, and obtain various types of data of the interpolation point; The interference degree acquisition module specifically includes: obtaining the disturbance degree of a type of data of the target point, the mean of the type of data of all neighboring target points corresponding to the target point, and the standard deviation of the type of data of all neighboring target points corresponding to the target point; obtaining the initial interference degree of the target point according to the disturbance degree of all types of data of the target point; and correcting the initial interference degree of the target point by the distribution of the neighboring target points of the target point in the feature space to obtain the interference degree of the target point.

Images