CN115831276A - Space classification and partition method and system for soil materials - Google Patents
Space classification and partition method and system for soil materials Download PDFInfo
- Publication number
- CN115831276A CN115831276A CN202211377313.XA CN202211377313A CN115831276A CN 115831276 A CN115831276 A CN 115831276A CN 202211377313 A CN202211377313 A CN 202211377313A CN 115831276 A CN115831276 A CN 115831276A
- Authority
- CN
- China
- Prior art keywords
- soil
- grid
- useful layer
- classification
- test data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
In order to solve the technical problem that the traditional soil material partitioning method has great difficulty in fully expressing the characteristics of the space distribution rule of the soil material when the space partitioning is carried out, the embodiment of the invention provides a soil material space classification partitioning method and a soil material space classification partitioning system, wherein the method comprises the following steps: establishing a soil geological reservoir grid model according to the soil exploration test data; judging whether the soil exploration test data in the grids of the soil geological reservoir grid model meet the soil quality technical index or not, and assigning values to each grid in the soil geological reservoir grid model according to the judgment result to obtain a first assigned soil model; dividing each grid meeting the technical index of the soil quality in the first assigned soil model into a plurality of stratum subareas and selecting a useful layer from the plurality of stratum subareas; reclassifying each useful layer to obtain a useful layer area; and analyzing the spatial position and the volume ratio of the spatial classification region in the useful layer region according to the geological features of the stratum to perform generalized treatment on the useful layer region to obtain a spatial classification partition of the soil material.
Description
Technical Field
The invention relates to a space classification and partition method and system for soil materials.
Background
The soil material is used as a natural building material in hydropower engineering construction, and can be used as a soil material or a raw material for preventing seepage, contacting clay and fixing walls of slotted holes under the condition that the quality index is qualified or the requirement can be met through engineering measures.
In order to find out the conditions of spatial distribution, reserves, quality, exploitation and the like of soil materials, basic geological conditions and stratigraphic division of the soil materials are mostly found out by utilizing exploration measures such as drilling holes and vertical shafts, and the quality of each layer of soil is analyzed and researched by sampling tests in the drilling holes and the vertical shafts.
The geological conditions revealed by the drilling holes and the vertical shaft and the stratum interface are spatial mark points, and each test result has a plurality of test parameters (such as data of water content, plasticity index, particle content and the like) besides spatial X, Y and Z coordinates.
According to the traditional soil yard zoning method, test parameters of each layer are researched according to geological stratification through means of site survey, statistical charts, histograms, tangent maps, space maps and the like, and plane zoning is carried out. However, the thickness and the elevation of each layer of soil material also have partition characteristics in space, the distribution characteristics of the soil material cannot be completely represented by plane partitions, the difficulty of the traditional soil material partition method in space partition is high, the space partitions are difficult to display in a two-dimensional graph, and the space distribution rule characteristics of the soil material are difficult to fully express; and the quality of the soil is limited by a plurality of indexes, and a large number of test parameters for representing the soil characteristic indexes are not used up.
Disclosure of Invention
In order to solve the technical problem that the traditional soil material partitioning method is difficult to fully express the characteristics of the space distribution rule of the soil material when the space partitioning is carried out, the embodiment of the invention provides a soil material space classification partitioning method and a soil material space classification partitioning system.
The embodiment of the invention is realized by the following technical scheme:
in a first aspect, an embodiment of the present invention provides a method for spatially classifying and partitioning soil materials, including:
establishing a soil geological reservoir grid model according to the soil exploration test data;
judging whether the soil exploration test data in the grids of the soil geological reservoir grid model meet the soil quality technical index or not, and assigning values to each grid in the soil geological reservoir grid model according to the judgment result to obtain a first assigned soil model;
dividing each grid meeting the technical index of the soil quality in the first assigned soil model into a plurality of stratum subareas and selecting a useful layer from the plurality of stratum subareas;
reclassifying each useful layer to obtain a useful layer area;
and analyzing the spatial position and the volume ratio of the spatial classification region in the useful layer region according to the geological features of the stratum to perform generalized treatment on the useful layer region to obtain a spatial classification partition of the soil material.
Further, establishing a soil geological reservoir grid model according to the soil exploration test data; the method comprises the following steps:
acquiring soil test parameters;
setting a modeling range, loading soil test parameters, and constructing a soil geological reservoir grid object to obtain a regularized spatial grid;
assigning a value to the corresponding soil exploration test data grid in the regularized spatial grid by using the soil exploration test data;
performing initialized attribute assignment on null value grids in the regularized spatial grids after the earth material exploration test data assignment is adopted;
and smoothing the discrete data of the null value grid to obtain the soil geological reservoir grid model.
Further, dividing the soil model subjected to the first assignment into a plurality of stratum partitions and selecting a useful layer from the plurality of stratum partitions; the method comprises the following steps:
and for each stratum partition, judging whether the ratio of the first assignment region conforming to the quality technical index to the stratum partition of the first assignment region containing the quality technical index is greater than a specified value, and if so, judging the stratum partition as a useful layer.
Further, classifying each useful layer to obtain a useful layer area; the method comprises the following steps:
selecting index data which accords with engineering processing conditions from the soil exploration test data, the burial depth index data and the thickness index data of the useful layer according to the engineering processing conditions 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 in the secondary range section of the grid of each useful layer is in the corresponding data standard range required by the engineering processing condition, and if the classification index data in the secondary range section of the grid of one useful layer is in the corresponding data standard range required by the engineering processing condition, judging that the classification index data in the secondary range section of the grid of one useful layer is a useful layer area unit; the collection of all the useful layer area units is the useful layer area.
Further, assigning a value to each grid in the soil geological reservoir grid model according to a judgment result; the method comprises the following steps:
and if the grid judgment result is that the grid meets the technical indexes of the soil quality, the first assignment of the grid is 1, otherwise, the grid is 0.
Further, the initializing the attribute assignments comprises: and (4) initializing the attribute of the soil exploration test data by using a DSI algorithm and then assigning to a null value grid.
Further, assigning a value to the corresponding soil exploration test data grid in the regularized spatial grid by using soil exploration test data; the method comprises the following steps:
assigning discrete soil exploration test data to a grid at the same spatial position as the soil exploration test data by utilizing attribute rendering; wherein, a point set is created during assignment;
smoothing discrete data of the null value grid to obtain a soil geological reservoir grid model; the method comprises the following steps:
initializing the discrete data of the null value grid and performing attribute interpolation iterative computation, wherein the point set is used as constraint before the attribute interpolation iterative computation.
Further, setting a modeling range, loading soil test parameters, and constructing a soil geological reservoir grid object to obtain a regularized spatial grid; the method comprises the following steps:
determining a reservoir grid model modeling range by taking the top surface as the maximum elevation with the drilling test data points and the bottom surface as the minimum elevation with the drilling test data points;
importing soil test parameters into a reservoir grid model;
and selecting the top surface and the bottom surface of the boundary of the reservoir grid model, setting the grid subdivision quantity, and establishing a regularized space grid.
Further, the soil test parameters comprise the maximum particle size index, crushed stone gravel content larger than 5mm after compaction, particle content smaller than 0.075mm, clay content smaller than 0.005mm, plasticity index, permeability coefficient after compaction, natural water content, organic matter content, water-soluble salt content, silicon-iron-aluminum ratio and/or soil dispersibility.
In a second aspect, an embodiment of the present invention provides a soil material spatial classification partitioning system, including:
the model establishing unit is used for establishing a soil geological reservoir grid model according to the soil exploration test data;
the judging unit is used for judging whether the soil exploration test data in the grids of the soil geological reservoir grid model meet the soil quality technical index or not, and assigning values to each grid in the soil geological reservoir grid model according to the judging result to obtain a first assigned soil model;
the useful layer screening unit is used for dividing each grid which accords with the technical index of the soil quality in the soil model of the first assignment into a plurality of stratum subareas and selecting a useful layer from the stratum subareas;
a useful layer area unit for reclassifying each useful layer to obtain a useful layer area; and
and the space classification partition unit is used for analyzing the space position and the volume ratio of the space classification area in the useful layer area according to the geological features of the stratum to perform generalized processing on the useful layer area so as to obtain the space classification partition of the soil material.
Compared with the prior art, the embodiment of the invention has the following advantages and beneficial effects:
according to the space classification and partitioning method for the soil material, disclosed by the embodiment of the invention, a soil material geological reservoir grid model is established according to soil material exploration test data; judging whether the soil exploration test data in the grids of the soil geological reservoir grid model meet the soil quality technical index or not, and assigning values to each grid in the soil geological reservoir grid model according to the judgment result to obtain a first assigned soil model; dividing each grid meeting the technical index of the soil quality in the first assigned soil model into a plurality of stratum subareas and selecting a useful layer from the plurality of stratum subareas; reclassifying each useful layer to obtain a useful layer area; the method comprises the steps of analyzing the spatial position and the volume ratio of a spatial classification region in a useful layer region according to geological features of a stratum, carrying out generalized processing on the useful layer region to obtain a spatial classification partition of the soil material, and solving the technical problem that the traditional soil material partition method is difficult to fully express the spatial distribution rule features of the soil material when spatial partition is carried out.
Drawings
In order to more clearly illustrate the technical solutions of the exemplary embodiments of the present invention, the drawings that are required to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present invention and therefore should not be considered as limiting the scope, and that for those skilled in the art, other related drawings can be obtained from these drawings without inventive effort.
Fig. 1 is a schematic flow chart of a method for spatially classifying and partitioning soil materials.
Fig. 2 is a schematic structural diagram of a space classification partition system of soil.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to examples and accompanying drawings, and the exemplary embodiments and descriptions thereof are only used for explaining the present invention and are not meant to limit the present invention.
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention. However, it will be apparent to one of ordinary skill in the art that: it is not necessary to employ these specific details to practice the present invention. In other instances, well-known structures, circuits, materials, or methods have not been described in detail so as not to obscure the present invention.
Throughout the specification, reference to "one embodiment," "an embodiment," "one example" or "an example" means: the particular features, structures, or characteristics described in connection with the embodiment or example are included in at least one embodiment of the invention. Thus, the appearances of the phrases "one embodiment," "an embodiment," "one 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 any suitable combination and/or sub-combination in one or more embodiments or examples. Further, those of ordinary skill in the art will appreciate that the illustrations provided herein are for illustrative purposes 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.
In the description of the present invention, the terms "front", "rear", "left", "right", "upper", "lower", "vertical", "horizontal", "upper", "lower", "inner", "outer", etc. indicate orientations or positional relationships based on those shown in the drawings, and are only for convenience of description and simplicity of description, but do not indicate or imply that the device or element being referred to must have a particular orientation, be constructed in a particular orientation, and be operated, and therefore, should not be construed as limiting the scope of the present invention.
Examples
In order to solve the technical problem that the conventional soil material partitioning method has great difficulty in fully expressing the characteristics of the spatial distribution rule of the soil material when performing spatial partitioning, in a first aspect, an embodiment of the present invention provides a method for spatially classifying and partitioning a soil material, which is shown in fig. 1 and includes:
s1, establishing a soil geological reservoir grid model according to soil exploration test data;
s2, judging whether the soil exploration test data in the grids of the soil geological reservoir grid model meet the soil quality technical index or not, and assigning values to each grid in the soil geological reservoir grid model according to the judgment result to obtain a first assigned soil model;
s3, dividing each grid meeting the technical index of the soil quality in the first assigned soil model into a plurality of stratum subareas and selecting a useful layer from the plurality of stratum subareas;
s4, classifying each useful layer to obtain a useful layer area;
and S5, analyzing the spatial position and the volume ratio of the spatial classification area in the useful layer area according to the geological features of the stratum, and performing generalized treatment on the useful layer area to obtain the spatial classification subarea of the soil material.
Therefore, the embodiment of the invention establishes the soil geological reservoir grid model according to the soil exploration test data; judging whether the soil exploration test data in the grids of the soil geological reservoir grid model meet the soil quality technical index or not, and assigning values to each grid in the soil geological reservoir grid model according to the judgment result to obtain a first assigned soil model; dividing each grid meeting the technical index of the soil quality in the first assigned soil model into a plurality of stratum subareas and selecting a useful layer from the plurality of stratum subareas; reclassifying each useful layer to obtain a useful layer area; the method comprises the steps of analyzing the spatial position and the volume ratio of a spatial classification region in a useful layer region according to geological features of a stratum, carrying out generalized processing on the useful layer region to obtain a spatial classification partition of the soil material, and solving the technical problem that the traditional soil material partition method is difficult to fully express the spatial distribution rule features of the soil material when spatial partition is carried out.
Further, establishing a soil geological reservoir grid model according to the soil exploration test data; the method comprises the following steps:
acquiring soil test parameters;
setting a modeling range, loading soil test parameters, and constructing a soil geological reservoir grid object to obtain a regularized spatial grid;
assigning a value to the corresponding soil exploration test data grid in the regularized spatial grid by using the soil exploration test data;
performing initialized attribute assignment on null value grids in the regularized space grids after the earth exploration test data assignment is adopted;
and smoothing the discrete data of the null value grid to obtain the soil geological reservoir grid model.
Further, dividing the soil model subjected to the first assignment into a plurality of stratum partitions and selecting a useful layer from the plurality of stratum partitions; the method comprises the following steps:
and for each stratum partition, judging whether the ratio of the first assignment region conforming to the quality technical index to the stratum partition of the first assignment region containing the quality technical index is greater than a specified value, and if so, judging the stratum partition as a useful layer.
Further, classifying each useful layer to obtain a useful layer area; the method comprises the following steps:
selecting index data which accords with engineering processing conditions from the soil exploration test data, the burial depth index data and the thickness index data of the useful layer according to the engineering processing conditions as classification index data;
dividing the data range of each classification index data into a plurality of secondary range intervals;
for example, the data range of a certain index is 1-100, and the data range can be divided into ten sub-range intervals of 1-9, 10-19,20-29 \8230 \ 8230:, etc. at intervals of 10.
Judging whether the classification index data in the secondary range section of the grid of each useful layer is in the corresponding data standard range required by the engineering processing condition, and if the classification index data in the secondary range section of the grid of one useful layer is in the corresponding data standard range required by the engineering processing condition, judging that the classification index data in the secondary range section of the grid of one useful layer is a useful layer area unit; the collection of all the useful layer area units is the useful layer area.
The examples are given above. Judging whether the data range in the grid of each useful layer is in the ten secondary range intervals, and comparing the data range in the grid of the useful layer with the data standard range, for example, taking classification index data which is in the range of 1-9 in all the grids of the useful layer and is in the corresponding data standard range required by the engineering processing conditions as the useful layer area of which the data is in the range of 1-9; likewise, useful-layer regions in the interval 10-19,20-29, etc. can be obtained in this way.
The active layer areas of each sub-range region may be assigned to a value, such as 1-9, 2, 8230, 10-99, etc., for ease of processing. The useful layer area unit of the same mark constitutes the useful layer area under the mark.
Further, assigning a value to each grid in the soil geological reservoir grid model according to a judgment result; the method comprises the following steps:
and if the grid judgment result is that the grid meets the technical indexes of the soil quality, the first assignment of the grid is 1, otherwise, the grid is 0.
Further, the initializing the attribute assignments comprises: and (4) initializing the attribute of the soil exploration test data by using a DSI algorithm and then assigning to a null value grid.
Further, evaluating the corresponding soil exploration test data grids in the regular space grids by adopting the soil exploration test data; the method comprises the following steps:
assigning discrete soil exploration test data to a grid at the same spatial position as the soil exploration test data by utilizing attribute rendering; wherein, a point set is created during assignment;
smoothing discrete data of the null value grid to obtain a soil geological reservoir grid model; the method comprises the following steps:
initializing the discrete data of the null value grid and performing attribute interpolation iterative computation, wherein the point set is used as constraint before the attribute interpolation iterative computation.
Further, setting a modeling range, loading soil test parameters, and constructing a soil geological reservoir grid object to obtain a regularized spatial grid; the method comprises the following steps:
determining a reservoir grid model modeling range by taking the top surface as the maximum elevation with the drilling test data points and the bottom surface as the minimum elevation with the drilling test data points;
importing soil test parameters into a reservoir grid model;
and selecting the top surface and the bottom surface of the boundary of the reservoir grid model, setting the grid subdivision quantity, and establishing a regularized space grid.
Further, the soil test parameters comprise the maximum particle size index, crushed stone gravel content larger than 5mm after compaction, particle content smaller than 0.075mm, clay content smaller than 0.005mm, plasticity index, permeability coefficient after compaction, natural water content, organic matter content, water-soluble salt content, silicon-iron-aluminum ratio and/or soil dispersibility.
Illustratively, the method for classifying and partitioning the space of the earth material field in the embodiment of the invention comprises the following steps:
1. and layering the soil layers according to the geological conditions revealed by the drilling shaft, and constructing a layering interface of each layer of soil through layering mark points of the drilling shaft.
2. Let D be a non-empty set of n-ary ordered arrays (n ≧ 1), and f be some defined rule (which can be implemented by C language).
If for each ordered array (x 1, x2, \ 8230;, xn) ∈ D, there is a uniquely determined real number y corresponding to it by the corresponding rule f, the functional expression is y = f (x 1, x2, \ 8230;, xn), (x 1, x2, \ 8230;, xn) ∈ D. The variable x1, x2, \ 8230, xn is the test data of the soil material, and y is the technical index of the soil material quality.
The technical scheme of test data partitioning is explained for the convenience of understanding by taking the technical indexes of the quality of the impermeable soil as examples.
Specifically, the establishing of the soil geological reservoir grid model according to the soil exploration test data comprises the following steps:
1) Obtaining soil test parameters
In the survey regulations of natural building materials of hydroelectric engineering (NB/T10235-2019), the appendix F technical requirements on the quality of soil materials, the maximum particle size index x1, the crushed stone and gravel content x2 which are larger than 5mm after compaction, the particle content x3 which is smaller than 0.075mm, the clay content x4 which is smaller than 0.005mm, the plasticity index x5, the permeability coefficient x6 after compaction, the natural water content x7, the organic matter content x8, the water-soluble salt content x9, the silicon-iron-aluminum ratio x10 and the soil dispersibility x11.
2) Setting the modeling range and loading discrete test parameters
The planar range includes the range of all drilled wells to be analyzed, the top surface being the maximum elevation at which drilled test data points exist, and the bottom surface being the minimum elevation at which drilled test data points exist. Various test parameters are introduced and attached to the exploratory object (stored in a record of the borehole object).
The impervious soil quality technical index model comprises all the drilling shaft ranges to be analyzed, and test data (x 1, x2, \8230; x 11) which are loaded in the model range, are all explored and acquired along with exploration
3) Constructing geological reservoir grid (Voxet) objects
And selecting the top surface and the bottom surface of the model boundary, setting the mesh subdivision quantity in three directions according to the hardware limit and the evaluation precision of the counting machine, and establishing a regularized space mesh. No value is assigned within the grid.
4) Geological reservoir grid (Voxet) valuation
The discrete test data (which are attached to the exploration object) are assigned to the grid where they are located (in the same spatial position) by using an attribute rendering method. And when the creation point set is selected during assignment, the assignment point set can be fixed as a hard constraint and used during interpolation. There may be multiple data in the same grid and the assignment will be calculated according to the chosen method (net-to-gross, average, percentile, other calculations). The technical indexes of the quality of the impermeable earth materials are subjected to arithmetric homogenization (power, geometric, harmonic and inverse distance are also added). After assignment is completed, the discrete data value is only matched with the grid where the discrete data value is located, and the grid which does not contain the actually measured discrete data is a null value.
5) Initializing attributes
And (4) initializing the attribute of the test data by using a DSI algorithm, and assigning the attribute to a null value grid which does not contain the actually measured discrete data. So far, all grids of the model are matched with non-null attribute values.
6) Attribute object interpolation
And (4) performing iterative calculation on the initialized test data by using a DSI algorithm. The purpose of the attribute interpolation iterative computation is to smooth the initialized discrete data. It should be noted that the discrete data point set created when the model mesh is assigned is used as a constraint before the interpolation operation, so that the expression of the known geological conditions is not affected. And after further exploration data is collected, new discrete data can be assigned to the corresponding grids to serve as new constraint points for interpolation calculation, so that the model can be adjusted.
7) Assignment of technical index y of impermeable soil quality (similar to judgment of contact clay quality and slotted hole solid wall soil quality)
And (2) creating a y attribute object in the established voxel object, wherein the attribute is still a null value, judging a standard in a natural building material survey procedure of hydropower engineering (NB/T10235-2019), judging test data (x 1, x2, \8230; x 11) through a judgment statement based on C language programming, and forming attribute assignment of y (meeting the technical index y =1 of the impermeable earth material quality and not meeting the technical index y =0 of the impermeable earth material quality).
Quality of anti-seepage soil
(1) The fine-grained soil material homogenizing dam is programmed as follows: if (x 4> =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; }
(2) The fine-grained soil material homogenizing dam is programmed as follows: if (x 4> =15& & x4< =40& & x5> =10& & x5< =20& & x6<1/100000& & x7> = "optimum moisture content" -2& & x7< = "optimum moisture content" +3& & x8<2& & x9<3& & x10> =2& & x10< =4& & x11= "non-dispersive soil") { y =1; else { y =0; }
(3) The weathered soil is programmed as: if (x 1<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-dispersive soil") { y =1; else { y =0; }
(4) The crushed (gravel) stone type soil material is programmed as follows: if (x 1< =150& & x2> =20& & x2< =50& & x3>15& & x4> =8| x4> =6& & x5>6& & x6<1/100000& & x7> = "optimal water content" -2& & x7 = "optimal water content" +3& & x8<2& & x9<3& & x10> =2& & x10 & =4& & x11= "non-dispersive soil") { y =1; else { y =0; }
Contact clay mass
The soil-soluble organic fertilizer is characterized by comprising the following components, by weight, more than 5mm of particle content x1, less than 0.075mm of particle content x2, less than 0.005mm of particle content x3, plasticity index x4, maximum particle size x5, silicon-iron-aluminum ratio x6, permeability coefficient x7, allowable ratio drop x8, organic matter content x9, water-soluble salt content x10, natural water content x11 and soil dispersibility x12.
Programming is as follows: if (x 1<10& & x2>60& & x3> =20& & x4>10& & x5> =20& & x 5& & & 40& & x6> =2& & x 6& & =4& & x7<1/1000000& & x8>5& & x9<2& & x10<3& & x11> "optimum moisture content" & & x11 & "=" optimum moisture content "+3& & x12=" non-dispersible soil ") { y =1; else { y =0; }
Quality of slotted hole wall-fixing soil
The grain content is more than 5mm x1, less than 0.075mm x2, less than 0.005mm x3, plasticity index x4, silicon-iron-aluminum ratio x5, pH value x6, activity index x7, and organic matter content x8
Programming is as follows: if (x 1<10& & x2>30& & x3> =15& & x4>17& & x5> =3& & x5< =4& & x6>7& & x7<1& & x8< 1) { y =1; else { y =0; }.
3. The useful layer of the soil material is determined.
The whole Voxet model is a, the terrain surface is used as a top boundary, the area below the Voxet terrain surface is extracted as b, and the area b is partitioned according to the stratum interface (the layered interface constructed in the step 1) in the area (b 1, b2, 8230; bn). Because most of the earth materials in the actual engineering can not completely meet the requirements of the technical indexes of the quality of the impermeable earth materials, namely y =1, and the partially unsatisfied indexes can change the technical indexes of the quality through engineering measures, the attribute of the partial technical index of the quality of the impermeable earth materials is corrected, the attribute of the corrected technical index of the quality of the impermeable earth materials y2 is established, test data (x 1, x2, \8230;, x 11) is judged through distinguishing sentences based on C language programming, and attribute assignment of y2 is formed (the technical index of the quality of the corrected impermeable earth materials y2=1 is met, the technical index of the quality of the corrected impermeable earth materials y2=0 is not met), and the correction at the position refers to the relaxation limitation of the condition of the part 7 in the step 2. The region of y2=1 in each stratum partition (b 1, b2 \8230;, bn) is extracted, and the useful layer (b 1, b2 \8230;, bn) in (b 1, b2 \8230;, bn) is determined by a volume ratio (e.g.: the ratio of the region volume satisfying y2=1 in b1 to the total volume of b 1) > 70% (this value is adjustable).
4. Useful layer space classification
Reclassifying the mesh of the useful layer (b 1, b2 8230;)
Wherein the buried depth and the thickness are relatively key classification indexes.
Obtaining a burial Depth index Depth: the elevation coordinate (z 1) of the terrain surface is vertically projected to grids of a useful layer (bx 1, bx2 \8230;), each grid has an elevation coordinate (z 2), and a buried Depth index Depth is obtained through z1-z 2.
Thickness index, thick: and vertically projecting the elevation coordinates of the top surface and the bottom surface of the stratum of the useful layer to grids, and calculating the difference of the elevation coordinates of each grid to obtain the Thickness index Thickness of the useful layer.
The (x 1, x2, \ 8230;, (x 11, xn1, xn 2) index of the experimental data was selected as the index with large variability or critical index. And classifying all the selected indexes, wherein the classification includes the condition of summarizing all the indexes. According to the classification rule, based on C language programming, the attribute assignment of y3 (y 3 is the quantitative discrimination value of useful layer classification) is formed through a conditional statement, and y3= f (x 1, x2, \ 8230;, x 11) (x 1, x2, \ 8230;, xn) ∈ D. The variable x1, x2, \ 8230, xn is test data of the soil material, and y3 is a soil material partition distinguishing value. And (3) carrying out region extraction by using different y3 values to obtain useful region regions (y 31, y32, y33, \8230; y3 n) meeting different y3 values, wherein the finer the index classification is, the more the regions are divided, and the sum of all the y3 regions is equal to the classified useful region (y 31+ y32+ \8230; + y3n = b 1).
5. And (4) partitioning the space of the soil material.
And analyzing the spatial position and the volume proportion of the spatial classification regions (y 31, y32, y33, \8230; y3 n) in the useful layer region b1 through the volume ratio and the geological features of the stratum, generalizing partial regions in the region, generalizing from y31 downwards, combining partial regions of y32, y33, \8230; y3n by y31, and then combining partial regions of y33, y34, \8230; y3n by y32, thereby completing the generalization of all the regions, wherein the sum of the generalized regions is still equal to the useful layer region b1. The area after the generalization forms the space subarea of the soil material.
And (5) repeating the processes 3 and 4 in other useful layer areas bn to complete the spatial classification and partition of all the soil materials.
In a second aspect, an embodiment of the present invention provides a soil material spatial classification partitioning system, which is shown in fig. 2, and includes:
the model establishing unit is used for establishing a soil geological reservoir grid model according to the soil exploration test data;
the judging unit is used for judging whether the soil exploration test data in the grids of the soil geological reservoir grid model meet the soil quality technical index or not, and assigning values to each grid in the soil geological reservoir grid model according to the judging result to obtain a first assigned soil model;
the useful layer screening unit is used for dividing each grid which accords with the technical index of the soil quality in the first assigned soil model into a plurality of stratum subareas and selecting a useful layer from the stratum subareas;
a useful layer area unit for reclassifying each useful layer to obtain a useful layer area; and
and the space classification partition unit is used for analyzing the space position and the volume ratio of the space classification area in the useful layer area according to the geological features of the stratum to perform generalized processing on the useful layer area so as to obtain the space classification partition of the soil material.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are merely exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (10)
1. A space classification and partition method of soil materials is characterized by comprising the following steps:
establishing a soil geological reservoir grid model according to the soil exploration test data;
judging whether the soil exploration test data in the grids of the soil geological reservoir grid model meet the soil quality technical index or not, and assigning values to each grid in the soil geological reservoir grid model according to the judgment result to obtain a first assigned soil model;
dividing each grid meeting the technical index of the soil quality in the first assigned soil model into a plurality of stratum subareas and selecting a useful layer from the plurality of stratum subareas;
reclassifying each useful layer to obtain a useful layer area;
and analyzing the spatial position and the volume ratio of the spatial classification region in the useful layer region according to the geological features of the stratum to perform generalized treatment on the useful layer region to obtain a spatial classification partition of the soil material.
2. The method of claim 1, wherein the earth material geological reservoir grid model is created based on earth material exploration test data; the method comprises the following steps:
acquiring soil test parameters;
setting a modeling range, loading soil test parameters, and constructing a soil geological reservoir grid object to obtain a regularized spatial grid;
assigning a value to the corresponding soil exploration test data grid in the regularized spatial grid by using the soil exploration test data;
performing initialized attribute assignment on null value grids in the regularized space grids after the earth exploration test data assignment is adopted;
and smoothing the discrete data of the null value grid to obtain a soil geological reservoir grid model.
3. The method of claim 2, wherein the first assigned earth model is divided into a plurality of stratigraphic sections and the useful layers are selected from the plurality of stratigraphic sections; the method comprises the following steps:
and for each stratum partition, judging whether the ratio of the first assignment region conforming to the quality technical index to the stratum partition of the first assignment region containing the quality technical index is greater than a specified value, and if so, judging the stratum partition as a useful layer.
4. The method of claim 2, wherein each useful layer is reclassified to obtain a useful layer area; the method comprises the following steps:
selecting index data which accords with engineering processing conditions from the soil exploration test data, the burial depth index data and the thickness index data of the useful layer according to the engineering processing conditions 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 in the secondary range section of the grid of each useful layer is in the corresponding data standard range required by the engineering processing condition, and if the classification index data in the secondary range section of the grid of one useful layer is in the corresponding data standard range required by the engineering processing condition, judging that the classification index data in the secondary range section of the grid of one useful layer is a useful layer area unit; the collection of all the useful layer area units is the useful layer area.
5. The method of spatially classifying and partitioning earthen material as claimed in claim 2, wherein each of said grids in said earthen material geological reservoir grid model is assigned a value in accordance with a result of said determination; the method comprises the following steps:
and if the grid judgment result is that the grid meets the technical indexes of the soil quality, the first assignment of the grid is 1, otherwise, the grid is 0.
6. The method for spatial classification and partition of earth material according to any one of claims 2-5, wherein the initializing attribute assignment comprises: and (4) initializing the attribute of the soil exploration test data by using a DSI algorithm and then assigning to a null value grid.
7. The method for spatially classifying and partitioning earth material according to any one of claims 2 to 5, wherein earth exploration test data is used to assign a value to a corresponding earth exploration test data grid in said regularized spatial grid; the method comprises the following steps:
assigning discrete soil exploration test data to a grid at the same spatial position as the soil exploration test data by utilizing attribute rendering; wherein, a point set is created during assignment;
smoothing discrete data of the null value grid to obtain a soil geological reservoir grid model; the method comprises the following steps:
initializing the discrete data of the null value grid and performing attribute interpolation iterative computation, wherein the point set is used as constraint before the attribute interpolation iterative computation.
8. The method for spatial classification and partition of the soil material according to any one of claims 2 to 5, wherein a modeling range is set, soil material test parameters are loaded, and a soil geological reservoir grid object is constructed to obtain a regularized spatial grid; the method comprises the following steps:
determining a reservoir grid model modeling range by taking the top surface as the maximum elevation with the drilling test data points and the bottom surface as the minimum elevation with the drilling test data points;
importing soil test parameters into a reservoir grid model;
and selecting the top surface and the bottom surface of the boundary of the reservoir grid model, setting the grid subdivision quantity, and establishing a regularized space grid.
9. The method of any one of claims 2 to 5, wherein the soil testing parameters include maximum particle size index, crushed gravel content greater than 5mm after compaction, particle content less than 0.075mm, clay content less than 0.005mm, plasticity index, permeability coefficient after compaction, natural water content, organic matter content, water soluble salt content, silicon-iron-aluminum ratio and/or soil dispersibility.
10. A soil material spatial classification zoning system, comprising:
the model establishing unit is used for establishing a soil geological reservoir grid model according to the soil exploration test data;
the judging unit is used for judging whether the soil exploration test data in the grids of the soil geological reservoir grid model meet the soil quality technical index or not, and assigning values to each grid in the soil geological reservoir grid model according to the judging result to obtain a first assigned soil model;
the useful layer screening unit is used for dividing each grid which accords with the technical index of the soil quality in the first assigned soil model into a plurality of stratum subareas and selecting a useful layer from the stratum subareas;
a useful layer area unit for reclassifying each useful layer to obtain a useful layer area; and
and the space classification partition unit is used for analyzing the space position and the volume ratio of the space classification area in the useful layer area according to the geological features of the stratum to perform generalized processing on the useful layer area so as to obtain the space classification partition of the soil material.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211377313.XA CN115831276A (en) | 2022-11-04 | 2022-11-04 | Space classification and partition method and system for soil materials |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211377313.XA CN115831276A (en) | 2022-11-04 | 2022-11-04 | Space classification and partition method and system for soil materials |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115831276A true CN115831276A (en) | 2023-03-21 |
Family
ID=85526649
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211377313.XA Pending CN115831276A (en) | 2022-11-04 | 2022-11-04 | Space classification and partition method and system for soil materials |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115831276A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117808376A (en) * | 2024-02-29 | 2024-04-02 | 杭州钱塘产城开发集团有限公司 | Construction engineering quality detection supervision method and system |
-
2022
- 2022-11-04 CN CN202211377313.XA patent/CN115831276A/en active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117808376A (en) * | 2024-02-29 | 2024-04-02 | 杭州钱塘产城开发集团有限公司 | Construction engineering quality detection supervision method and system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Dell'Arciprete et al. | Comparison of three geostatistical methods for hydrofacies simulation: a test on alluvial sediments | |
| Mariethoz et al. | Reconstruction of incomplete data sets or images using direct sampling | |
| Hou et al. | Assessing quality of urban underground spaces by coupling 3D geological models: The case study of Foshan city, South China | |
| CN111766629B (en) | Method for identifying and describing deep carbonate karst structure | |
| Tinker | Building the 3-D jigsaw puzzle: Applications of sequence stratigraphy to 3-D reservoir characterization, Permian Basin | |
| GB2523869A (en) | Representation of unstructured grids | |
| CN106023313A (en) | 3D geologic digital reconnaissance report generating method | |
| CN110244021B (en) | Stratum layering method based on anisotropic interpolation | |
| CN105321204B (en) | A kind of three-dimensional geological grid body Model construction method | |
| Verbovšek et al. | Morphometric properties of dolines in Matarsko podolje, SW Slovenia | |
| Elzain et al. | Lithofacies modeling of Late Jurassic in upper Ulayyah reservoir unit at central Saudi Arabia with inference of reservoir characterization | |
| CN115831276A (en) | Space classification and partition method and system for soil materials | |
| CN115494560A (en) | Novel modeling method based on combination of CT (computed tomography) and drilling data | |
| Grunwald et al. | Soil layer models created with profile cone penetrometer data | |
| Yeste et al. | Integrating outcrop and subsurface data to improve the predictability of geobodies distribution using a 3D training image: A case study of a Triassic Channel–Crevasse-splay complex | |
| CN107507179B (en) | Rock-soil mass quantitative analysis method based on GOCAD | |
| Cuddy et al. | The application of fuzzy logic and genetic algorithms to reservoir characterization and modeling | |
| Abdelkarim et al. | High-resolution lithofacies and porosity modeling of the mixed siliciclastic–carbonate deposits of the Burdigalian Dam Formation, Eastern Saudi Arabia | |
| Mehdipour et al. | The Best Scenario for Geostatistical Modeling of Porosity in the Sarvak Reservoir in an Iranian Oil Field, Using Electrofacies, Seismic Facies, and Seismic Attributes | |
| Albarrán-Ordás et al. | The D i models method: geological 3-D modeling of detrital systems consisting of varying grain fractions to predict the relative lithological variability for a multipurpose usability | |
| Stafleu et al. | Models Created as 3‐D Cellular Voxel Arrays | |
| Ji et al. | An automated method to build 3D multi-scale geological models for engineering sedimentary layers with stratum lenses | |
| CN112292714B (en) | Grid partition based on fault radiation | |
| Erzeybek Balan | Characterization and modeling of paleokarst reservoirs using multiple-point statistics on a non-gridded basis | |
| Berrino et al. | 3D gravity inversion by growing bodies and shaping layers at Mt. Vesuvius (Southern Italy) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |