An interactive Quaternary stratigraphic stratification method based on statistical laws

Abstract

The present invention relates to the field of geotechnical engineering geological survey, and provides an interactive Quaternary stratigraphic stratification method based on statistical laws, which consists of several links in two aspects: statistical stratification and three-dimensional correction. The statistical stratification is classified and counted based on the borehole catalog data (including at least the rock and soil name and depth) to obtain the statistical characteristics of the spatial distribution of each rock and soil layer; then the large layers are preliminarily set in order from high to low certainty, and the rock and soil layers that do not belong to the large layers are allowed to be specified at the same time, and the statistical results are corrected based on this constraint condition. The interactive operation is repeated in this way to complete the statistical stratification of all rock and soil layers and obtain standard strata that conforms to the distribution law of Quaternary stratigraphic layers. The three-dimensional correction is used to quickly check the problems in the connection relationship between boreholes, automatically calculate the area exceeding the allowable inclination angle, guide the user to check and correct, and correct the unreasonable connection relationship in an interactive manner.

Description

Technical Field

[0001] The present invention relates to the field of geotechnical engineering geological survey, and provides a fast and interactive stratification method that allows intervention based on professional knowledge, based on the basic data obtained from drilling exploration of Quaternary strata to generate a two-dimensional geological map or create a three-dimensional geological model based on the spatial stratigraphic sequence relationship of the strata, thereby improving the efficiency and quality of the stratification work, and the stratification results meet the requirements of three-dimensional geological modeling. Background Art

[0002] Quaternary stratification refers to assigning correct stratigraphic codes to each rock and soil layer exposed by drilling. This code can not only correctly reflect the stratigraphic relationship between strata within a single hole, but also reasonably describe the connection relationship between strata between adjacent holes, ensuring the professional rationality of the spatial relationship of the stratification results.

[0003] Quaternary stratification is one of the basic tasks of geotechnical engineering geological survey. In essence, it is necessary to infer the three-dimensional spatial form of the Quaternary strata based on one-dimensional drilling data. So far, Quaternary stratification mainly relies on the rationality check of the generated two-dimensional geological profiles, which can only ensure the rationality of the connection relationship between the strata in each profile, and it is difficult to ensure the correctness of the connection relationship between profiles from a three-dimensional perspective. The spatial rationality of the stratification results is the basis of three-dimensional geological modeling, which directly affects the practical feasibility of the promotion and application of three-dimensional geological modeling technology.

[0004] The difficulty of Quaternary stratification lies in the spatial variability of geological conditions, including property changes and spatial morphological changes. The former refers to the fact that the properties of the same rock and soil layer vary with location, which makes the results revealed in different parts different, leading to differences in the cataloging and test results of the same rock and soil layer; the latter refers to the spatial distribution, in which multiple rock and soil layers with the same properties are formed by alternating deposition, which are repeatedly exposed in drilling holes and have complex connection relationships between adjacent holes.

[0005] At present, two-dimensional stratification is mainly used in production practice, which often refers to ensuring the rationality of the stratification results and connection relationships on a given section, ignoring the relationship between the sections. In fact, the three-dimensional rationality of the stratification results cannot be guaranteed, and the requirements of three-dimensional geological modeling cannot be met. For this reason, when three-dimensional geological modeling technology is applied to the field of geotechnical engineering geological survey, a variety of three-dimensional stratification methods have emerged, such as a stratification method and device proposed in the publication number CN110119587A, and an automatic stratum comparison stratification method proposed in the publication number CN118296393A. Their common feature is that the names of the rocks and soils revealed by each borehole and their corresponding depths are used to complete the stratification at one time using mathematical methods. The common problem of these methods is that they are completely dependent on the data revealed by the boreholes, cannot effectively utilize professional knowledge, and the stratum division does not have the necessary interaction and iteration process, which is not enough to deal with the uncertainty problem of the spatial distribution of the actual Quaternary strata. In view of the limited drilling data and the uncertainty of geological bodies, professional knowledge and manual judgment cannot be completely replaced in actual work. Manual intervention, trial and error, and iteration are often indispensable links in address work. None of the above-mentioned layered methods that ignore professional intervention are accepted by engineering practice. The interaction and essence of the present invention is the necessary manual intervention in the computer-aided operation process to solve the problem of geological body uncertainty. Summary of the invention

[0006] The present invention aims at solving the technical problems existing in the prior art and provides an interactive Quaternary stratigraphic layering method based on statistical laws, comprising: S100. Based on the borehole cataloging data of the Quaternary rock and soil layers (at least only the rock and soil name and depth), the rock and soil layers with the same name in the selected boreholes are counted by depth to obtain preliminary statistical stratification data, and it is allowed to be modified when problems are found in the subsequent statistical stratification and three-dimensional correction process, forming part of the interactive operation. The rock and soil layers with the same name refer to the basic units with completely identical characteristic indicators such as rock and soil name, genesis and properties, which are determined by means of on-site cataloging, in-situ testing and indoor testing.

[0007] S200. Based on the statistical results of rock and soil layers, the rock and soil layers exposed by different boreholes and with similar characteristic indicators are merged, and the rock and soil layers with smaller numbers are merged into other rock and soil layers in accordance with professional requirements. In this way, the designated rock and soil layers in all the boreholes participating in the statistics are merged in batches at one time to improve the efficiency of basic data processing.

[0008] S300. The statistical stratification results are arranged in order of average bottom depth from small to large, reflecting the upper and lower layer relationship of rock and soil layers from a statistical perspective. Geological personnel divide the upper and lower layer relationship of each rock and soil layer to the greatest extent based on the information revealed by the statistical results and combined with professional knowledge, and assign stratigraphic codes. The specific implementation principle is achieved through the interactive operation of the two links of "include" and "exclude".

[0009] S400, statistical stratification does not consider the connection relationship between strata, and cannot fully guarantee the correctness of the stratification results. It needs to be further checked and corrected through three-dimensional stratification. According to the statistical stratification results, the three-dimensional connection relationship between the rock and soil layers with the same code between each borehole is established in the order of large layer first and sub-layer, and the triangulation technology is used to connect them into a spatial surface, and the average inclination of the surface is calculated. Then the allowable deviation value of the inclination is preset, and the triangulated network that exceeds the allowable value is clearly displayed through the color cloud map. The necessary graphic interaction operations, modification of the connection relationship and even stratification are performed to form a correct layer relationship, eliminate the abnormal fluctuation of the ground surface, and obtain the final stratification result that meets the requirements.

[0010] Preferably, the S300 includes the following steps: S301. Geologists divide large layers according to the order of certainty from high to low based on preliminary statistical stratification data and professional knowledge, and designate a certain rock and soil layer as the bottom of a certain large layer. All rock and soil layers from the bottom of this rock and soil layer to the bottom of the previous large layer (when the designated rock and soil layer is the bottom of the first large layer, the bottom of the previous large layer is the ground surface) are sublayers of this large layer and are included in this large layer, that is, the "include" link.

[0011] S302. After designating a certain rock and soil layer as the bottom of a large layer, you can also designate a certain rock and soil layer below it (with a greater depth) as a "removal" marker layer, that is, the rock and soil layers of the marker layer and below it in each borehole do not belong to the designated large layer, so as to remove the rock and soil layers that are "mixed" into the designated large layer due to the same name and do not belong to the large layer in terms of stratigraphic relationship.

[0012] S303, using the manually set "include" and "exclude" as conditions, updating the statistical results, thereby obtaining a list of different rock and soil layers, i.e., statistical stratification data. "Include" and "exclude" are set in turn according to the large layers to be divided, and the statistical stratification data are continuously updated until the operation of all rock and soil layers is completed.

[0013] S304. In the process of updating the statistical stratification list in sequence according to the large layer, it is agreed to assign the rock and soil layer stratum code in the form of Nj, where N is the large layer and j is the sublayer, including but not limited to Arabic numerals, English letters, and Roman characters. All rock and soil layers with the same N value belong to the same large layer, and the one with the largest j value is the bottom of the large layer.

[0014] S305. According to the assigned stratum code, all rock and soil layers of all boreholes are traversed. When errors that do not comply with professional rules can be clearly identified, corrections are made according to the requirements of professional rules. For strata suspected of having problems, prompt information is given to guide manual processing.

[0015] Preferably, the rock and soil layer refers to the basic unit with similar rock and soil properties. In addition to the required rock and soil name, it can also be a combination of any one or more indicators such as genesis, age, plasticity, density and weathering degree. If the rock and soil names are the same but other indicators are different, they can be regarded as different rock and soil layers.

[0016] Preferably, the preliminary statistical stratification data is calculated according to the rock and soil layer division indicators (rock and soil name, age, genesis, plasticity, density and weathering degree, etc.), and the average bottom depth, average thickness, number of rock and soil layers, the total number of drill holes revealing the rock and soil layers, etc. are used as the initial data for stratification.

[0017] Preferably, when using "include" to set the large layer bottom, the differences in causes are considered first, and the rock and soil layers representing different causes are used as the large layer bottom; the second is to consider the exposure conditions, and the rock and soil layers with the same causes and relatively stable distribution are used as the large layer bottom.

[0018] Preferably, the use of "removal" is to specify the rock and soil layers that do not belong to the large layer (located below the bottom of the specified large layer), giving priority to the differences in causes and rock and soil types, and specifying the rock and soil layers to be removed, whose causes and types are different from the specified large layer.

[0019] Preferably, N and j in the formation code Nj are both Arabic numerals; As a preference, when checking and automatically correcting the statistical stratification results of each borehole one by one according to the standard stratum, the principle that sub-strata belonging to different large strata do not appear alternately is followed; The above process is implemented using computer software technology, wherein the drilling data is stored in a relational database, and a corresponding software system is compiled to complete statistical stratification and three-dimensional correction.

[0020] The present invention has the following beneficial effects: The present invention has three beneficial effects: upgrading the traditional operation mode, improving the operation efficiency, and revolutionary improving the quality of the result data. The upgrading of the traditional operation mode refers to upgrading the two-dimensional stratification used in the three-dimensional strata to three-dimensional stratification to ensure the rationality of the three-dimensional spatial relationship; improving the operation efficiency is manifested in statistics and batch operations of all boreholes and all strata, instead of the current one-by-one and repeated modification of a single borehole, which greatly reduces the amount of repeated work. The batch operation includes "merge", "include", "eliminate", "code assignment" and other links, and completes the corresponding operations for all boreholes in batches at one time, without finding each borehole and operating one by one. The verification example shows that for conventional sites with dozens of boreholes, the stratification time can be shortened from the current 2 to 3 days to less than 1 hour; and for large sites with hundreds or thousands of boreholes, the stratification time can be reduced from the current several months to less than 2 days, and the efficiency is more prominent. When the traditional two-dimensional stratification results are used to create a three-dimensional geological model, the modeling work often fails because the correctness of the three-dimensional spatial relationship cannot be guaranteed. The three-dimensional stratification of the present invention completely solves this pain point problem and is a key link in the promotion of three-dimensional modeling technology. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] Figure 1 It is a schematic flow chart of the interactive Quaternary stratigraphic layering method based on statistical laws of the present invention; Figure 2 A schematic diagram of batch merging of strata is provided for the present invention; Figure 3 It is a preliminary statistical layered schematic diagram of the present invention; Figure 4 This is a schematic diagram of refreshing the statistical stratification results of the present invention; Figure 5 This is a schematic diagram of the result of "removing" the rock and soil layer in the present invention; Figure 6 It is a schematic diagram of the statistical stratified assignment results of the present invention; Figure 7 This is a schematic diagram of the results of artificially correcting the drilling stratification of the present invention; Figure 8 It is a schematic diagram of the three-dimensional correction inspection of the adjacent hole layer relationship of the present invention; Fig. 9 It is a schematic diagram of the final layered result of the three-dimensional correction of the present invention. DETAILED DESCRIPTION

[0022] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention. In addition, the technical features in the various embodiments or single embodiments provided by the present invention can be arbitrarily combined with each other to form a feasible technical solution. This combination is not subject to the constraints of the sequence of steps and / or the structural composition mode, but must be based on the ability of ordinary technicians in this field to achieve. When the combination of technical solutions is contradictory or cannot be achieved, it should be considered that such a combination of technical solutions does not exist and is not within the scope of protection required by the present invention.

[0023] In view of the background technology, the present invention proposes an interactive Quaternary stratification method based on statistical laws in accordance with the long-standing work idea of ​​mutual verification of exploration data and professional knowledge in geological work. The core of the method is "interaction", which allows different stratifications to be tried according to professional knowledge, and the optimal results are obtained after necessary corrections are made according to the corresponding statistical rationality and spatial rationality, meeting the requirements of three-dimensional geological modeling. The method is implemented through software technology to ensure the speed of the stratification process. Even under relatively complex site conditions, the stratification work can still be completed within 1 hour, greatly improving work efficiency.

[0024] The present invention provides an interactive quaternary stratification method based on statistical laws. Figure 1 , including the following steps: S100: Based on the borehole catalog data of the Quaternary rock and soil layers (at least only the rock and soil name and depth), the rock and soil layers with the same name in the selected boreholes are counted by depth, such as Figure 3 As shown, the average bottom depth, average thickness, number of rock and soil layers, total number of drill holes revealing the rock and soil layers, etc., of miscellaneous fill, plain fill, silty clay, gravel and sand, are calculated as the initial data for stratification.

[0025] Specifically, rock and soil layers refer to the basic units with similar rock and soil properties. In addition to the required rock and soil names, they can also be a combination of any one or more indicators such as genesis, age, plasticity, density and degree of weathering. If the rock and soil names are the same but other indicators are different, they can be regarded as different rock and soil layers.

[0026] S200: Based on the statistical results of rock and soil layers, the rock and soil layers with similar properties revealed by different boreholes are merged, such as Figure 2As shown, a smaller amount of miscellaneous fill is incorporated into the plain fill, so that the miscellaneous fill of three boreholes can be incorporated into the plain fill in batches at one time, without having to find out each borehole and operate one by one, thereby improving the efficiency of basic data processing.

[0027] S300: The statistical stratification results are arranged in order of average bottom depth from small to large, reflecting the upper and lower layer relationship of rock and soil layers from a statistical perspective. Geological personnel divide the upper and lower layer relationship of each rock and soil layer to the greatest extent based on the information revealed by the statistical results and combined with professional knowledge, and assign stratigraphic codes. The specific implementation principle is achieved through the interactive operation of the two links of "include" and "exclude".

[0028] S301: Geologists divide the large layers according to their professional knowledge in descending order of certainty, e.g. Figure 3 As shown, plain fill, gravel sand, sandy clay soil and slightly weathered mixed granite are designated as the bottom of a large layer, and all rock and soil layers from the bottom of this rock and soil layer to the bottom of the previous large layer (when the designated rock and soil layer is the bottom of the first large layer, the bottom of the previous large layer is the ground surface) are sublayers of this large layer and are included in this large layer, that is, the "include" link.

[0029] S302: Figure 3 As shown in the figure, according to the average bottom depth obtained by statistical stratification, the rock and soil names designated as the bottom of the large layer are assigned large layer numbers in accordance with the principle of sorting from small to large, and are assigned integer values ​​of 1 to 4. They are manually selected and can be left blank, but the selected data are kept in ascending order, increasing from top to bottom (interval is greater than or equal to 1), and it is agreed that the rock and soil layer stratum code is assigned in the form of Nj, such as Figure 5 As shown, for example, miscellaneous fill 1-1 and plain fill 1-2 are respectively the first and second sublayers of the first large layer, and the sublayer numbers are arranged from top to bottom.

[0030] S303: Re-count the assigned stratum code in a large-layer-first manner, refresh the statistical stratification results, and then divide the upper and lower layers of each rock and soil layer to the greatest extent possible based on professional knowledge, and allow "elimination" of rock and soil layers that do not belong to the large layer, and use this as a constraint condition to correct the statistical results. Figure 4 As shown in the figure, after refreshing the statistical stratification results, the fully weathered mixed granite is mixed into the previous large layer 3. By "removing", it is specified that the fully weathered mixed granite is a rock and soil layer that does not belong to this large layer, such as Figure 5 As shown, after the designation, the fully weathered mixed granite in all the boreholes is removed from the large layer 3 and is located below the bottom of the designated large layer 3. Similar to "merging", removal is also a batch operation of the corresponding boreholes to improve the work efficiency.

[0031] Furthermore, the interactive operation is repeated in this way, and the bottom of the large layer is manually specified and the rock and soil layer is "eliminated" interactively until the rock and soil name record values ​​of all the boreholes are assigned stratigraphic codes with stratigraphic sequence relationships, completing the interactive statistical stratification, such as Figure 6 As shown, a standard stratum was obtained that conforms to the distribution law of Quaternary strata.

[0032] S304: According to the assigned stratum code, all rock and soil layers of all boreholes are traversed, and the statistical stratification results of each borehole are checked one by one according to the standard stratum. When it can be clearly determined that there are errors that do not comply with professional rules, the erroneous stratum that does not conform to the upper and lower layer relationship is automatically corrected; Furthermore, for the strata suspected of having problems, the computer will list the questions that cannot be judged, give prompt information, and guide manual processing. Figure 7 As shown, professionals can refer to the stratification results of adjacent holes with high reliability. When the catalog data of borehole ZK1 cannot clearly define the stratigraphic boundary within the depth range of 4.85-14.3m, after referring to the adjacent hole ZK2 with high reliability of stratification results on the left, when it is considered that the rock and soil layers can be further divided within the above-mentioned depth range, the in-situ test results and geotechnical test results of ZK1 are used to determine the need for further identification of the rock and soil names and to divide the lithological boundary between medium sand and fine sand (fine sand 14.2~16.5m, medium sand 16.5~19.1m).

[0033] S400, statistical stratification does not consider the three-dimensional spatial relationship of the strata, and does not guarantee the correctness of the results. It needs to be further checked and corrected through three-dimensional correction. According to the statistical stratification results, the three-dimensional connection relationship between the rock and soil layers with the same code between the boreholes is established in the order of large layers first and sub-layers, and the triangulation technology is used to connect them into a spatial surface, and the average inclination of the surface is calculated. Then the inclination allowable deviation value is preset, and the triangulated network that exceeds the allowable value is clearly displayed through the color cloud map. The necessary graphic interaction operations, modification of connection relationships and even stratification are carried out to form a correct layer relationship and eliminate abnormal fluctuations in the stratum surface until the stratification results meet the quality requirements of the three-dimensional geological modeling for the input data.

[0034] like Figure 8 As shown in the figure, the three-dimensional connection relationship between the rock and soil layers with the same code in each borehole is established in the order of large layer first and then sub-layer. The triangulation technology is used to connect them into spatial surfaces, and the average inclination angle of each surface is calculated. Then, the allowable deviation range of the inclination angle is preset, and the spatial fluctuation of the rock and soil layer 1-2 is clearly displayed through the color cloud map. The triangulated network that exceeds the allowable value range is regarded as a potential abnormal fluctuation (such as Figure 8The system quickly finds the problem and presents it prominently through cloud rendering. Then the borehole rock and soil layer 1-2 that causes abnormal layer connection relationship is removed, and the rock and soil layer 1-2 of the same name in the borehole is moved to the specified layer surface to form a correct layer connection relationship, eliminate the abnormal fluctuation of the layer surface, and obtain the final stratification result that meets the requirements, such as Fig. 9 shown.

[0035] Furthermore, the above process is implemented using computer software technology, wherein the drilling data is stored in a relational database, and a corresponding software system is compiled to complete statistical stratification and three-dimensional correction.

[0036] The embodiment of the present invention provides an interactive Quaternary stratification method based on statistical laws, the core of which is "interaction", which allows different stratifications to be tried according to professional knowledge, and the optimal result is obtained after necessary corrections are made according to the corresponding statistical rationality and spatial rationality, meeting the requirements of three-dimensional geological modeling. The method is implemented through software technology to ensure the speed of the stratification process. Even under relatively complex site conditions, the stratification work can still be completed in a relatively short time, greatly improving work efficiency.

[0037] It should be noted that in the above embodiments, the description of each embodiment has its own emphasis, and for parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0038] It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Furthermore, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0039] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded computer, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0040] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0041] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0042] Although the preferred embodiments of the present invention have been described, those skilled in the art may make other changes and modifications to these embodiments once they have learned the basic creative concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications that fall within the scope of the present invention.

[0043] Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include these modifications and variations.

Claims

1. An interactive Quaternary stratigraphic stratification method based on statistical laws, characterized in that: It consists of two aspects: statistical stratification and three-dimensional correction. The statistical stratification includes: Based on the borehole cataloging data of the Quaternary rock and soil layers (at least only the name and depth of the borehole rock and soil, and can also include the cause, properties, weathering degree, etc.), the rock and soil layers with the same name of the selected boreholes are classified and counted by depth to obtain the statistical characteristics of the spatial distribution of each rock and soil layer; then, based on professional knowledge and in order of certainty from high to low, the large layers are preliminarily set in sequence, and the rock and soil layers that do not belong to the large layer are allowed to be specified at the same time, and the statistical results are corrected based on this constraint condition. In this way, the statistical stratification of all rock and soil layers is completed by repeated interactive operations according to the large layers, and the standard strata that conforms to the distribution law of the Quaternary strata are obtained; then, the statistical stratification results of each borehole are checked one by one according to the standard strata, and the erroneous strata that do not conform to the upper and lower layer relationships are corrected, and the questions that the computer cannot judge are listed, and the professionals are guided to check and correct them to ensure that the stratification results of each borehole conform to the geological laws. The statistical stratification process can also merge the rock and soil layers with similar properties exposed by different boreholes according to the statistical results to improve the efficiency of basic data processing; Three-dimensional correction is used to check and correct the connection relationship between drill holes, automatically calculate the area that exceeds the allowable inclination angle, guide the user to check and correct, and correct unreasonable connection relationships in an interactive manner.

2. The interactive Quaternary stratigraphic layering method based on statistical laws according to claim 1, characterized in that: The rock and soil layers with the same name of the selected boreholes are classified and counted by depth to obtain preliminary statistical stratification data. The average bottom depth, average thickness, number of rock and soil layers, total number of boreholes revealing the rock and soil layers, etc. are counted according to the rock and soil layer division indicators (rock and soil name, age, genesis, plasticity, density and weathering degree, etc.). The statistical results are arranged in order from small to large according to the average bottom depth. As the initial data for stratification, the rock and soil layer refers to the basic unit with similar rock and soil properties. In addition to the rock and soil name that must be included, it can also be a combination of any one or more indicators such as genesis, age, plasticity, density and weathering degree. If the rock and soil names are the same but other indicators are different, they can be regarded as different rock and soil layers.

3. The interactive quaternary stratification method based on statistical laws according to claim 2 is characterized in that: The statistical stratification results are arranged in order of average bottom depth from small to large, reflecting the upper and lower layer relationship of the rock and soil layers from a statistical perspective. The technical personnel divide the upper and lower layer relationship of each rock and soil layer to the greatest extent based on the information revealed by the statistical results and combined with professional knowledge, and assign stratum codes. The specific implementation principle is realized by the interactive operation of the two links of "include" and "exclude", among which "include" is to designate a certain rock and soil layer as the bottom of a certain large layer according to professional knowledge, and all rock and soil layers from the bottom of this rock and soil layer to the bottom of the previous large layer (when the designated rock and soil layer is the bottom of the first large layer, the bottom of the previous large layer is the surface) are The sublayers belonging to the large layer are included in the large layer; "exclusion" means specifying the rock and soil layers that do not belong to the large layer (located below the bottom of the specified large layer). After the specification, the rock and soil layer and all the rock and soil layers below it in all boreholes do not belong to the specified large layer, so as to eliminate other rock and soil layers that are "mixed" into the large layer with the same name and different stratigraphic relationships in the statistical results, and update the statistical results. This is repeated, and the interactive operation method of manually specifying the bottom of the large layer and the rock and soil layers to be eliminated is used to complete the stratification of all rock and soil layers and assign stratigraphic codes, and obtain standard stratigraphic layers that conform to the distribution law of Quaternary stratigraphy.

4. The interactive Quaternary stratigraphic layering method based on statistical laws according to claim 3, characterized in that: According to the statistical results of the rock and soil layers, the rock and soil layers with small numbers can be merged into other rock and soil layers in accordance with professional requirements, so as to complete the merging processing of different borehole rock and soil layers in batches at one time and improve the working efficiency.

5. The interactive quaternary stratification method and system based on statistical laws according to claim 3 is characterized in that: It is agreed that the rock and soil layer stratigraphic codes are assigned in the form of Nj, where N is the main layer and j is the sublayer, including but not limited to Arabic numerals, English letters, and Roman characters. All rock and soil layers with the same N value belong to the same main layer, and the one with the largest j value is the bottom of the main layer.

6. The interactive Quaternary stratigraphic layering method based on statistical laws according to claim 5, characterized in that: According to the assigned stratigraphic code, all rock and soil layers of all boreholes are traversed. When it can be clearly determined that there are errors that do not comply with professional rules (such as stratigraphic intersections), they are corrected automatically in accordance with professional rules. For stratigraphic layers suspected of having problems, prompt information is given to guide manual processing.

7. The interactive quaternary stratification method based on statistical laws according to claim 1, characterized in that: According to the statistical stratification results, the three-dimensional connection relationship between the rock and soil layers with the same code in each borehole is established in the order of large layer first and sub-layer. The triangulation technology is used to connect them into a spatial surface, and the average inclination of the surface is calculated. Then the allowable deviation value of the inclination is preset, and the triangulated network that exceeds the allowable value is clearly displayed through a color cloud map. The necessary graphic interaction operations, modification of the connection relationship and even stratification are carried out to obtain the final result that meets the requirements.

8. The interactive Quaternary stratigraphic layering method based on statistical laws according to claim 1, characterized in that: The drilling data is stored in a relational database, and a corresponding software system is compiled to complete statistical stratification and three-dimensional correction.

Images