A method and system for spatial classification and partitioning of soil materials

Abstract

In order to solve the technical problem that traditional soil material partitioning methods are difficult to fully express the spatial distribution regularity characteristics of soil materials during spatial partitioning, an embodiment of the present invention provides a spatial classification and partitioning method and system for soil materials, including: establishing a soil material geological reservoir grid model based on soil material exploration test data; judging whether the soil material exploration test data in the grids of the soil material geological reservoir grid model meet the soil material quality technical indicators, and assigning values ​​to each grid in the soil material geological reservoir grid model according to the judgment result to obtain a first assigned soil material model; dividing each grid that meets the soil material quality technical indicators in the first assigned soil material model into several stratigraphic partitions and selecting useful layers from the several stratigraphic partitions; reclassifying each useful layer to obtain a useful layer area; analyzing the spatial position and volume proportion of the spatial classification area in the useful layer area according to the stratigraphic geological characteristics, generalizing the useful layer area, and obtaining the spatial classification and partitioning of soil materials.

Description

Technical Field

[0001] The invention relates to a method and system for spatial classification and partitioning of soil materials. Background Art

[0002] As a natural building material, soil can be used as soil or raw material for anti-seepage, contact clay and slot wall reinforcement in hydropower project construction, provided that the quality indicators are qualified or the requirements can be met through engineering measures.

[0003] In order to find out the spatial distribution, reserves, quality, mining conditions and other conditions of soil materials, exploration methods such as drilling and vertical shafts are often used to find out its basic geological conditions and stratigraphic divisions, and sampling tests are carried out in boreholes and vertical shafts to analyze and study the quality of each layer of soil.

[0004] The geological conditions and stratigraphic interfaces revealed by boreholes and shafts serve as spatial landmarks. In addition to the spatial X, Y, and Z coordinates, each test result also contains multiple test parameters (such as water content, plasticity index, particle size content, and other data).

[0005] Traditional soil field zoning methods use field surveys, statistical charts, bar graphs, cross-section plots, and spatial diagrams to analyze the test parameters of each layer based on geological stratification and perform planar zoning. However, each layer of soil also exhibits spatial zoning characteristics based on thickness and elevation, and planar zoning cannot fully represent the distribution characteristics of the soil. Traditional soil field zoning methods are difficult to perform spatial zoning, and spatial zoning is difficult to display in two-dimensional diagrams, making it difficult to fully express the spatial distribution characteristics of the soil. Furthermore, soil quality is limited by multiple indicators, and a large number of test parameters that characterize soil characteristics are not fully utilized. Summary of the Invention

[0006] In order to solve the technical problem that traditional soil material partitioning methods are difficult to fully express the spatial distribution characteristics of soil materials during spatial partitioning, the embodiments of the present invention provide a spatial classification and partitioning method and system for soil materials.

[0007] The embodiments of the present invention are implemented through the following technical solutions:

[0008] In a first aspect, an embodiment of the present invention provides a method for spatial classification and partitioning of soil materials, comprising:

[0009] Establish a soil geological reservoir grid model based on soil exploration test data;

[0010] Determining whether soil exploration test data in a grid of a soil geological reservoir grid model meets soil quality technical indicators, and assigning values ​​to each grid in the soil geological reservoir grid model according to the determination result to obtain a first assigned soil model;

[0011] Dividing each grid in the first assigned soil material model that meets the soil material quality technical index into a plurality of stratum partitions and selecting useful layers from the plurality of stratum partitions;

[0012] Reclassify each useful layer to obtain the useful layer area;

[0013] According to the geological characteristics of the strata, the spatial classification areas are analyzed, and the spatial positions and volume proportions of the useful layer areas are generalized to obtain the spatial classification zoning of the soil materials.

[0014] Furthermore, a soil geological reservoir grid model is established based on the soil exploration test data; including:

[0015] Obtain soil material test parameters;

[0016] Set the modeling scope, load soil material test parameters, construct soil material geological reservoir grid objects, and obtain regularized spatial grids;

[0017] Using soil exploration test data to assign values ​​to corresponding soil exploration test data grids in the regularized spatial grid;

[0018] Initialize the attribute assignment for the null value grid in the regularized space grid after the soil exploration test data is assigned;

[0019] The discrete data of the null value grid are smoothed to obtain the soil geological reservoir grid model.

[0020] Furthermore, the first assigned soil material model is divided into a plurality of stratum partitions and useful layers are selected from the plurality of stratum partitions; including:

[0021] For each stratum partition, determine whether the ratio of the first value-assigned area that meets the quality technical indicators to the stratum partition containing the first value-assigned area that meets the quality technical indicators is greater than a specified value. If so, determine that the stratum partition is a useful layer.

[0022] Furthermore, each useful layer is reclassified to obtain a useful layer area; including:

[0023] According to the engineering treatment conditions, index data that meet the engineering treatment conditions are selected from the soil exploration test data, burial depth index data and thickness index data of the useful layer as classification index data;

[0024] Divide the data range of each classification indicator data into several sub-range intervals;

[0025] Determine whether the classification index data within the secondary range interval of the grid of each useful layer is within the corresponding data standard range required by the engineering treatment conditions. If the classification index data within the secondary range interval of the grid of a useful layer is within the corresponding data standard range required by the engineering treatment conditions, then determine that the classification index data within the secondary range interval of the grid of the useful layer is a useful layer area unit; the set of all useful layer area units is a useful layer area.

[0026] Furthermore, assigning a value to each grid in the soil geological reservoir grid model according to the judgment result includes:

[0027] If the judgment result of a grid is that it meets the technical indicators of soil quality, the first value of the grid is 1, otherwise it is 0.

[0028] Furthermore, the initialization attribute assignment includes: using a DSI algorithm to initialize the attributes of the soil exploration test data and then assigning the initial attributes to a null value grid.

[0029] Furthermore, soil exploration test data is used to assign values ​​to corresponding soil exploration test data grids in the regularized spatial grid; including:

[0030] Assigning discrete soil exploration test data to a grid having the same spatial location as the soil exploration test data using attribute rendering; wherein a point set is created during the assignment;

[0031] The discrete data of the null value grid is smoothed to obtain the soil geological reservoir grid model; including:

[0032] The discrete data of the null value grid is initialized and attribute interpolation iterative calculation is performed, and the point set is used as a constraint before the attribute interpolation iterative calculation.

[0033] Furthermore, the modeling scope is set, soil material test parameters are loaded, and soil material geological reservoir grid objects are constructed to obtain a regularized spatial grid; including:

[0034] The modeling range of the reservoir grid model is determined by taking the top surface as the maximum elevation where the drilling test data points exist, and the bottom surface as the lowest elevation where the drilling test data points exist;

[0035] Import soil material test parameters into the reservoir grid model;

[0036] Select the top and bottom surfaces of the reservoir grid model boundary to set the number of grid divisions and establish a regularized spatial grid.

[0037] Furthermore, soil material test parameters include maximum particle size index, content of crushed stone and gravel larger than 5 mm after compaction, content of particles smaller than 0.075 mm, content of clay smaller than 0.005 mm, plasticity index, permeability coefficient after compaction, natural moisture content, organic matter content, water-soluble salt content, silicon-iron-aluminum ratio and / or soil dispersibility.

[0038] In a second aspect, an embodiment of the present invention provides a spatial classification and zoning system for soil materials, including:

[0039] A model building unit is used to build a soil geological reservoir grid model based on soil exploration test data;

[0040] a judgment unit, configured to judge whether the soil exploration test data in the grids of the soil geological reservoir grid model meets the soil quality technical indicators, and assign values ​​to each grid in the soil geological reservoir grid model according to the judgment result to obtain a first assigned soil model;

[0041] a useful layer screening unit, configured to divide each grid in the first assigned soil material model that meets the soil material quality technical index into a plurality of stratum partitions and select useful layers from the plurality of stratum partitions;

[0042] A useful layer area unit, used for reclassifying each useful layer to obtain a useful layer area; and

[0043] The spatial classification and partitioning unit is used to analyze the spatial position and volume proportion of the spatial classification area in the useful layer area according to the geological characteristics of the stratum, generalize the useful layer area, and obtain the spatial classification and partitioning of the soil material.

[0044] Compared with the prior art, the embodiments of the present invention have the following advantages and beneficial effects:

[0045] A method for spatial classification and partitioning of soil materials in an embodiment of the present invention establishes a soil material geological reservoir grid model based on soil material exploration test data; determines whether the soil material exploration test data in the grids of the soil material geological reservoir grid model meet the soil material quality technical indicators, and assigns values ​​to each grid in the soil material geological reservoir grid model according to the determination result to obtain a first assigned soil material model; divides each grid that meets the soil material quality technical indicators in the first assigned soil material model into several stratigraphic partitions and selects useful layers from the several stratigraphic partitions; reclassifies each useful layer to obtain a useful layer area; analyzes the spatial position and volume proportion of the spatial classification area in the useful layer area according to stratigraphic geological characteristics, generalizes the useful layer area, and obtains the spatial classification and partitioning of the soil material, thereby solving the technical problem that the traditional soil material partitioning method is difficult to fully express the spatial distribution regularity characteristics of the soil material when performing spatial partitioning. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0047] Figure 1 This is a flowchart of the spatial classification and zoning method for soil materials.

[0048] Figure 2 This is a structural diagram of the spatial classification and zoning system for soil materials. DETAILED DESCRIPTION

[0049] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with examples and drawings. The exemplary embodiments of the present invention and their descriptions are only used to explain the present invention and are not intended to limit the present invention.

[0050] In the following description, numerous specific details are set forth to provide a thorough understanding of the present invention. However, it will be apparent to one skilled in the art that these specific details are not necessarily required to practice the present invention. In other embodiments, well-known structures, circuits, materials, or methods are not described in detail to avoid obscuring the present invention.

[0051] Throughout this specification, references to "one embodiment," "an embodiment," "an example," or "an example" mean that a particular feature, structure, or characteristic described in connection with the embodiment or example is included in at least one embodiment of the present invention. Therefore, appearances of the phrases "one embodiment," "an embodiment," "an example," or "an example" in various places throughout this specification are not necessarily all referring to the same embodiment or example. Furthermore, the particular features, structures, or characteristics may be combined in one or more embodiments or examples in any suitable combinations and / or subcombinations. Furthermore, it will be understood by those of ordinary skill in the art that the figures provided herein are for illustrative purposes only and are not necessarily drawn to scale. As used herein, the term "and / or" includes any and all combinations of one or more of the associated listed items.

[0052] In the description of the present invention, the terms "front", "back", "left", "right", "up", "down", "vertical", "horizontal", "high", "low", "inside", "outside", etc., indicating directions or positional relationships, are based on the directions or positional relationships shown in the accompanying drawings. They are only for the convenience of describing the present invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific direction, be constructed and operated in a specific direction. Therefore, they should not be understood as limiting the scope of protection of the present invention.

[0053] Example

[0054] In order to solve the technical problem that the traditional soil material partitioning method is difficult to fully express the spatial distribution characteristics of soil materials during spatial partitioning, the embodiment of the present invention provides a spatial classification and partitioning method for soil materials, referring to Figure 1 As shown, including:

[0055] S1. Establishing a soil geological reservoir grid model based on soil exploration test data;

[0056] S2 determines whether the soil exploration test data in the grid of the soil geological reservoir grid model meets the soil quality technical indicators, and assign values ​​to each grid in the soil geological reservoir grid model based on the judgment result to obtain a first assigned soil model;

[0057] S3. Dividing each grid in the first assigned soil material model that meets the soil material quality technical indicators into several stratum partitions and selecting useful layers from the several stratum partitions;

[0058] S4. reclassifying each useful layer to obtain a useful layer region;

[0059] S5. Analyze the spatial position and volume proportion of the spatial classification area in the useful layer area according to the geological characteristics of the strata, generalize the useful layer area, and obtain the spatial classification and partitioning of the soil material.

[0060] Therefore, the embodiment of the present invention establishes a soil geological reservoir grid model based on soil exploration test data; judges whether the soil exploration test data in the grids of the soil geological reservoir grid model meets the soil quality technical indicators, and assigns values ​​to each grid in the soil geological reservoir grid model according to the judgment result to obtain a first assigned soil model; divides each grid that meets the soil quality technical indicators in the first assigned soil model into several stratigraphic partitions and selects useful layers from the several stratigraphic partitions; reclassifies each useful layer to obtain a useful layer area; analyzes the spatial position and volume proportion of the spatial classification area in the useful layer area according to the stratigraphic geological characteristics, generalizes the useful layer area, and obtains the spatial classification area of ​​the soil, which solves the technical problem that the traditional soil partitioning method is difficult to fully express the spatial distribution law characteristics of the soil when performing spatial partitioning.

[0061] Furthermore, a soil geological reservoir grid model is established based on the soil exploration test data; including:

[0062] Obtain soil material test parameters;

[0063] Set the modeling scope, load soil material test parameters, construct soil material geological reservoir grid objects, and obtain regularized spatial grids;

[0064] Using soil exploration test data to assign values ​​to corresponding soil exploration test data grids in the regularized spatial grid;

[0065] Initialize the attribute assignment for the null value grid in the regularized space grid after the soil exploration test data is assigned;

[0066] The discrete data of the null value grid are smoothed to obtain the soil geological reservoir grid model.

[0067] Furthermore, the first assigned soil material model is divided into a plurality of stratum partitions and useful layers are selected from the plurality of stratum partitions; including:

[0068] For each stratum partition, determine whether the ratio of the first value-assigned area that meets the quality technical indicators to the stratum partition containing the first value-assigned area that meets the quality technical indicators is greater than a specified value. If so, determine that the stratum partition is a useful layer.

[0069] Furthermore, each useful layer is reclassified to obtain a useful layer area; including:

[0070] According to the engineering treatment conditions, index data that meet the engineering treatment conditions are selected from the soil exploration test data, burial depth index data and thickness index data of the useful layer as classification index data;

[0071] Divide the data range of each classification indicator data into several sub-range intervals;

[0072] For example, if the data range of a certain indicator is 1-100, this data range can be divided into ten sub-range intervals at intervals of 10, such as 1-9, 10-19, 20-29...

[0073] Determine whether the classification index data within the secondary range interval of the grid of each useful layer is within the corresponding data standard range required by the engineering treatment conditions. If the classification index data within the secondary range interval of the grid of a useful layer is within the corresponding data standard range required by the engineering treatment conditions, then determine that the classification index data within the secondary range interval of the grid of the useful layer is a useful layer area unit; the set of all useful layer area units is a useful layer area.

[0074] Continuing with the example above, determine whether the data range in each grid of the useful layer is within the ten sub-range intervals mentioned above, and compare the data range in the grid of the useful layer with the data standard range. For example, the classification indicator data in the range of 1-9 in all grids of the useful layer and within the corresponding data standard range required by the engineering treatment conditions are regarded as the useful layer area with data within the range of 1-9. Similarly, the useful layer areas within the ranges of 10-19, 20-29, and so on can be obtained in this way.

[0075] For ease of processing, the useful layer area of ​​each secondary range interval can be assigned a value label, such as 1-9 is labeled 1, 10-19 is labeled 2, ... 90-99 is labeled 10, etc. Useful layer area units with the same label constitute the useful layer area under the label.

[0076] Furthermore, assigning a value to each grid in the soil geological reservoir grid model according to the judgment result includes:

[0077] If the judgment result of a grid is that it meets the technical indicators of soil quality, the first value of the grid is 1, otherwise it is 0.

[0078] Furthermore, the initialization attribute assignment includes: using a DSI algorithm to initialize the attributes of the soil exploration test data and then assigning the initial attributes to a null value grid.

[0079] Furthermore, soil exploration test data is used to assign values ​​to corresponding soil exploration test data grids in the regularized spatial grid; including:

[0080] Assigning discrete soil exploration test data to a grid having the same spatial location as the soil exploration test data using attribute rendering; wherein a point set is created during the assignment;

[0081] The discrete data of the null value grid is smoothed to obtain the soil geological reservoir grid model; including:

[0082] The discrete data of the null value grid is initialized and attribute interpolation iterative calculation is performed, and the point set is used as a constraint before the attribute interpolation iterative calculation.

[0083] Furthermore, the modeling scope is set, soil material test parameters are loaded, and soil material geological reservoir grid objects are constructed to obtain a regularized spatial grid; including:

[0084] The modeling range of the reservoir grid model is determined by taking the top surface as the maximum elevation where the drilling test data points exist, and the bottom surface as the lowest elevation where the drilling test data points exist;

[0085] Import soil material test parameters into the reservoir grid model;

[0086] Select the top and bottom surfaces of the reservoir grid model boundary to set the number of grid divisions and establish a regularized spatial grid.

[0087] Furthermore, soil material test parameters include maximum particle size index, content of crushed stone and gravel larger than 5 mm after compaction, content of particles smaller than 0.075 mm, content of clay smaller than 0.005 mm, plasticity index, permeability coefficient after compaction, natural moisture content, organic matter content, water-soluble salt content, silicon-iron-aluminum ratio and / or soil dispersibility.

[0088] Exemplarily, the method for spatial classification and zoning of a soil material field according to an embodiment of the present invention includes:

[0089] 1. The soil layers are stratified according to the geological conditions revealed by the vertical drilling shaft, and the stratification interfaces of each layer of soil are constructed through the stratification marking points of the vertical drilling shaft.

[0090] 2. Let D be a non-empty set of n-element ordered arrays (n≥1), and f be a certain corresponding rule (which can be implemented in C language).

[0091] If for every ordered array (x1, x2, …, xn) ∈ D, there is a unique real number y corresponding to it according to the correspondence rule f, the function expression is y = f(x1, x2, …, xn), (x1, x2, …, xn) ∈ D. The variables x1, x2, …, xn are the test data of the soil material, and y is the technical indicator of the soil material quality.

[0092] In order to facilitate understanding, the technical scheme for partitioning test data is explained using the technical indicators of anti-seepage soil quality as an example.

[0093] Specifically, establishing a soil geological reservoir grid model based on soil exploration test data includes:

[0094] 1) Obtain soil material test parameters

[0095] In the "Regulations for the Investigation of Natural Building Materials for Hydropower Projects" (NB / T 10235-2019), Appendix F Technical Requirements for Soil Quality, maximum particle size index x1, crushed stone and gravel content larger than 5mm after compaction x2, particle content less than 0.075mm x3, clay content less than 0.005mm x4, plasticity index x5, permeability coefficient after compaction x6, natural moisture content x7, organic matter content x8, water-soluble salt content x9, silicon-iron-aluminum ratio x10, soil dispersibility x11.

[0096] 2) Set the modeling scope and load discrete test parameters

[0097] The plane range includes all the borehole shafts to be analyzed. The top plane is the highest elevation where there are borehole test data points, and the bottom plane is the lowest elevation where there are borehole test data points. Various test parameters are imported and attached to the exploration object (stored in the record of the borehole object).

[0098] The technical indicator model of anti-seepage soil quality includes the entire range of drilled shafts that need to be analyzed, and all exploration and test data obtained during exploration within the model range (x1, x2, ..., x11) are loaded.

[0099] 3) Constructing a geological reservoir grid (Voxet) object

[0100] Select the top and bottom surfaces of the model boundary, set the number of mesh divisions in three directions based on the computer hardware limitations and evaluation accuracy, and create a regularized spatial grid. No values ​​are assigned to this grid.

[0101] 4) Geological reservoir grid (Voxet) assignment

[0102] The discrete test data (attached to the exploration object) is assigned to the grid where it is located (same spatial position) using the attribute rendering method. If you choose to create a point set when assigning values, the assigned point set can be fixed as a hard constraint and used during interpolation. Multiple data may exist in the same grid, and the assignment will be calculated according to the selected method (net-to-gross ratio, homogenization, percentage points, other calculations). The technical indicators of the anti-seepage soil quality at this location are processed by arithmetic homogenization (in addition, power, geometric, harmonic, and inverse distance are also used). After the assignment is completed, the discrete data value is only matched to the grid where it is located. The grid that does not contain measured discrete data will be empty.

[0103] 5) Initialize properties

[0104] Using the DSI algorithm, the test data attributes are initialized and assigned to the empty value grids that do not contain measured discrete data. At this point, all grids in the model are matched with non-empty attribute values.

[0105] 6) Attribute object interpolation operation

[0106] Using the DSI algorithm, iterative calculations are performed on the initialized test data. The purpose of this iterative attribute interpolation calculation is to smooth the initialized discrete data. It is important to note that the discrete data point set created when assigning values ​​to the model grid before interpolation is used as a constraint to ensure that the representation of known geological conditions is not affected. Furthermore, as further exploration data is collected, new discrete data can be assigned to the corresponding grid and used as new constraint points for interpolation calculations to adjust the model.

[0107] 7) Assigning the y attribute value of the technical indicator of anti-seepage soil material quality (the quality of contact clay material and slot hole solid wall material is similar)

[0108] Create the y attribute object in the established Voxet object. The attribute is still empty. The judgment standard in the "Regulations for the Investigation of Natural Building Materials for Hydropower Projects" (NB / T 10235-2019) is based on C language programming. The test data (x1, x2, ..., x11) is judged by the judgment statement to form the attribute assignment of y (if the technical indicators of anti-seepage soil quality are met, y = 1; if the technical indicators of anti-seepage soil quality are not met, y = 0).

[0109] Quality of anti-seepage soil

[0110] ① The programming for a homogeneous dam of fine-grained soil is: if(x4>=10&&x4<=30&&x5>=7&&x5<=17&&x6<1 / 10000&&x7>="optimal moisture content"-2&&x7<="optimal moisture content"+3&&x8<5&&x9<3&&x10>=2&&x10<=4&&x11="non-dispersed soil"){y=1;}else{y=0;}

[0111] ② The programming for the homogeneous dam of fine-grained soil is: if(x4>=15&&x4<=40&&x5>=10&&x5<=20&&x6<1 / 100000&&x7>="optimal moisture content"-2&&x7<="optimal moisture content"+3&&x8<2&&x9<3&&x10>=2&&x10<=4&&x11="non-dispersed soil"){y=1;}else{y=0;}

[0112] ③ The programming for weathered soil is: if(x1<150&&x2>=20&&x2<=50&&x3>15&&x4>8&&x5>8&&x6<1 / 100000&&x7>="optimal moisture content"-2&&x7<="optimal moisture content"+3&&x8<2&&x9<3&&x10>=2&&x10<=4&&x11="non-dispersed soil"){y=1;}else{y=0;}

[0113] ④ The programming for crushed (gravel) stone soil is: if(x1<=150&&x2>=20&&x2<=50&&x3>15&&x4>=8||x4>=6&&x5>6&&x6<1 / 100000&&x7>="optimal moisture content"-2&&x7<="optimal moisture content"+3&&x8<2&&x9<3&&x10>=2&&x10<=4&&x11="non-dispersed soil"){y=1;}else{y=0;}

[0114] Contact clay material quality

[0115] Content of particles larger than 5mm x1, content of particles smaller than 0.075mm x2, content of particles smaller than 0.005mm x3, plasticity index x4, maximum particle size x5, silicon-iron-aluminum ratio x6, permeability coefficient x7, allowable gradient x8, organic matter content x9, water-soluble salt content x10, natural moisture content x11, soil dispersibility x12.

[0116] Programming: if(x1<10&&x2>60&&x3>=20&&x4>10&&x5>=20&&x5<=40&&x6>=2&&x6<=4&&x7<1 / 1000000&&x8>5&&x9<2&&x10<3&&x11>"Optimal moisture content"&&x11<="Optimal moisture content"+3&&x12="Non-dispersible soil"){y=1;}else{y=0;}

[0117] Slot wall reinforcement soil quality

[0118] Content of particles larger than 5mm x1, content of particles smaller than 0.075mm x2, content of particles smaller than 0.005mm x3, plasticity index x4, silicon-iron-aluminum ratio x5, pH value x6, activity index x7, organic matter content x8

[0119] The programming is: if (x1<10&&x2>30&&x3>=15&&x4>17&&x5>=3&&x5<=4&&x6>7&&x7<1&&x8<1){y=1;}else{y=0;}.

[0120] 3. Determination of the useful layer of soil material.

[0121] The entire Voxet model is a, with the terrain surface as the top boundary. The area below the Voxet terrain surface is extracted as b. Within this area, region b is partitioned (b1, b2…, bn) based on the stratum interface (the stratum interface constructed in step 1). Because most soil materials in actual projects do not fully meet the anti-seepage soil quality technical index requirements, that is, y = 1, some of the unsatisfied indicators can be improved through engineering measures. Therefore, the attributes of some anti-seepage soil quality technical indicators y are modified to establish the modified anti-seepage soil quality technical indicators y2. Based on C language programming, the test data (x1, x2, …, x11) are judged through a discriminant statement to form the attribute assignment of y2 (if the modified anti-seepage soil quality technical indicators are met, y2 = 1; if the modified anti-seepage soil quality technical indicators are not met, y2 = 0). This modification refers to relaxing the restrictions on some of the conditions in step 2, section 7 above. Extract the area where y2=1 in each stratum partition (b1, b2…, bn), and determine the useful layers (b1, b2…, bn) in (b1, b2…, bn) through the volume ratio (such as: the ratio of the volume of the area in b1 that satisfies y2=1 to the total volume of b1) > 70% (this value can be adjusted).

[0122] 4. Useful layer space classification

[0123] Reclassify the grids of the useful layers (b1, b2...)

[0124] Among them, burial depth and thickness are relatively key classification indicators.

[0125] Depth index is obtained: the elevation coordinate (z1) of the terrain surface is vertically projected to the grid of the useful layer (bx1, bx2...). Each grid has an elevation coordinate (z2), and the depth index Depth is obtained by z1-z2.

[0126] Thickness index: The elevation coordinates of the top and bottom surfaces of the useful layer are vertically projected onto the grid, and the elevation coordinate difference of each grid is calculated to obtain the useful layer thickness index Thickness.

[0127] Select the most significant or critical indicators (x1, x2, …, x11, xn1, xn2) from the test data. Classify all selected indicators, ensuring that all indicators are fully covered. Based on the classification rules, conditional statements are used in C language programming to assign the attribute value y3 (y3 is the quantitative discriminant value for the useful layer classification): y3 = f(x1, x2, …, x11)(x1, x2, …, xn) ∈ D. The variables x1, x2, …, xn represent the test data for the soil material, and y3 is the discriminant value for the soil material partitioning. Regions are extracted using different y3 values, resulting in useful layer regions (y31, y32, y33, …, y3n) that meet these different y3 values. The finer the indicator classification, the more partitions there are. The sum of all y3 regions equals the classified useful layer region (y31 + y32 + … + y3n = b1).

[0128] 5. Spatial partitioning of soil materials.

[0129] The spatial location and volume proportions of the spatially classified regions (y31, y32, y33, ..., y3n) within the useful layer region b1 were analyzed using volume ratios and stratum geological characteristics. Some of these regions were generalized, starting with y31 and generalizing downwards. Y31 merged parts of y32, y33, ..., y3n, and then y32 merged parts of y33, y34, ..., y3n. This generalized region remained equal to the useful layer region b1. These generalized regions formed the spatial partitions of the soil material.

[0130] Repeat the process of 3 and 4 in other useful layer areas bn to complete the spatial classification and partitioning of all soil materials.

[0131] In the second aspect, the embodiment of the present invention provides a spatial classification and partitioning system for soil materials, referring to Figure 2 As shown, including:

[0132] A model building unit is used to build a soil geological reservoir grid model based on soil exploration test data;

[0133] a judgment unit, configured to judge whether the soil exploration test data in the grids of the soil geological reservoir grid model meets the soil quality technical indicators, and assign values ​​to each grid in the soil geological reservoir grid model according to the judgment result to obtain a first assigned soil model;

[0134] a useful layer screening unit, configured to divide each grid in the first assigned soil material model that meets the soil material quality technical index into a plurality of stratum partitions and select useful layers from the plurality of stratum partitions;

[0135] A useful layer area unit, used for reclassifying each useful layer to obtain a useful layer area; and

[0136] The spatial classification and partitioning unit is used to analyze the spatial position and volume proportion of the spatial classification area in the useful layer area according to the geological characteristics of the stratum, generalize the useful layer area, and obtain the spatial classification and partitioning of the soil material.

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

Claims

1. A method for spatial classification and partitioning of soil materials, characterized in that: include: Establish a soil geological reservoir grid model based on soil exploration test data; Determining whether soil exploration test data in a grid of a soil geological reservoir grid model meets soil quality technical indicators, and assigning values ​​to each grid in the soil geological reservoir grid model according to the determination result to obtain a first assigned soil model; Dividing each grid in the first assigned soil material model that meets the soil material quality technical index into a plurality of stratum partitions and selecting useful layers from the plurality of stratum partitions; Reclassifying each useful layer to obtain a useful layer region, including: selecting, according to the engineering treatment condition, index data that meets the engineering treatment condition from soil exploration test data, burial depth index data, and thickness index data of the useful layer as classification index data; dividing the data range of each classification index data into a plurality of secondary range intervals; judging whether the classification index data within the secondary range interval of the grid of each useful layer is within the corresponding data standard range required by the engineering treatment condition; if the classification index data within the secondary range interval of the grid of a useful layer is within the corresponding data standard range required by the engineering treatment condition, then judging the classification index data within the secondary range interval of the grid of the useful layer as a useful layer region unit; and the set of all useful layer region units is the useful layer region; According to the geological characteristics of the strata, the spatial classification areas are analyzed, and the spatial positions and volume proportions of the useful layer areas are generalized to obtain the spatial classification zoning of the soil materials.

2. The method for spatial classification and zoning of soil materials according to claim 1, characterized in that: Establish a soil geological reservoir grid model based on soil exploration test data; including: Obtain soil material test parameters; Set the modeling scope, load soil material test parameters, construct soil material geological reservoir grid objects, and obtain regularized spatial grids; Using soil exploration test data to assign values ​​to corresponding soil exploration test data grids in the regularized spatial grid; Initialization attribute assignment is performed on the null-value grid in the regularized spatial grid after the soil exploration test data is assigned; the discrete data of the null-value grid is smoothed to obtain the soil geological reservoir grid model.

3. The method for spatial classification and zoning of soil materials according to claim 2, characterized in that: Dividing the first assigned soil material model into a plurality of stratum partitions and selecting useful layers from the plurality of stratum partitions; comprising: For each stratum partition, determine whether the ratio of the first value-assigned area that meets the quality technical indicators to the stratum partition containing the first value-assigned area that meets the quality technical indicators is greater than a specified value. If so, determine that the stratum partition is a useful layer.

4. The method for spatial classification and zoning of soil materials according to claim 2, characterized in that: Assigning a value to each grid in the soil geological reservoir grid model according to the judgment result; including: If the judgment result of a grid is that it meets the technical indicators of soil quality, the first value of the grid is 1, otherwise it is 0.

5. The method for spatial classification and partitioning of soil materials according to any one of claims 2 to 4, characterized in that: The initialization attribute assignment includes: using the DSI algorithm to initialize the attributes of the soil exploration test data and then assigning them to a null value grid.

6. The method for spatial classification and partitioning of soil materials according to any one of claims 2 to 4, characterized in that: Using soil exploration test data to assign values ​​to corresponding soil exploration test data grids in the regularized spatial grid; including: Assigning discrete soil exploration test data to a grid having the same spatial location as the soil exploration test data using attribute rendering; wherein a point set is created during the assignment; The discrete data of the null value grid is smoothed to obtain the soil geological reservoir grid model; including: The discrete data of the null value grid is initialized and attribute interpolation iterative calculation is performed, and the point set is used as a constraint before the attribute interpolation iterative calculation.

7. The method for spatial classification and partitioning of soil materials according to any one of claims 2 to 4, characterized in that: Set the modeling scope, load soil material test parameters, construct soil material geological reservoir grid objects, and obtain regularized spatial grids; including: The modeling range of the reservoir grid model is determined by taking the top surface as the maximum elevation where the drilling test data points exist, and the bottom surface as the lowest elevation where the drilling test data points exist; Import soil material test parameters into the reservoir grid model; Select the top and bottom surfaces of the reservoir grid model boundary to set the number of grid divisions and establish a regularized spatial grid.

8. The method for spatial classification and partitioning of soil materials according to any one of claims 2 to 4, characterized in that: Soil material test parameters include maximum particle size index, content of crushed stone and gravel larger than 5 mm after compaction, content of particles smaller than 0.075 mm, content of clay smaller than 0.005 mm, plasticity index, permeability coefficient after compaction, natural moisture content, organic matter content, water-soluble salt content, silicon-iron-aluminum ratio and / or soil dispersibility.

9. A spatial classification and zoning system for soil materials, characterized in that: include: A model building unit is used to build a soil geological reservoir grid model based on soil exploration test data; a judgment unit, configured to judge whether the soil exploration test data in the grids of the soil geological reservoir grid model meets the soil quality technical indicators, and assign values ​​to each grid in the soil geological reservoir grid model according to the judgment result to obtain a first assigned soil model; a useful layer screening unit, configured to divide each grid in the first assigned soil material model that meets the soil material quality technical index into a plurality of stratum partitions and select useful layers from the plurality of stratum partitions; The useful layer region unit is used to reclassify each useful layer to obtain a useful layer region, specifically comprising: selecting index data that meets the engineering treatment conditions from soil exploration test data, burial depth index data, and thickness index data of the useful layer as classification index data according to the engineering treatment conditions; dividing the data range of each classification index data into a plurality of secondary range intervals; judging whether the classification index data within the secondary range interval of the grid of each useful layer is within the corresponding data standard range required by the engineering treatment conditions; if the classification index data within the secondary range interval of the grid of a useful layer is within the corresponding data standard range required by the engineering treatment conditions, then judging the classification index data within the secondary range interval of the grid of the useful layer as a useful layer region unit; the set of all useful layer region units is a useful layer region; The spatial classification and partitioning unit is used to analyze the spatial position and volume proportion of the spatial classification area in the useful layer area according to the geological characteristics of the stratum, generalize the useful layer area, and obtain the spatial classification and partitioning of the soil material.

Images