Stratum sequence modeling method for stratigraphic range constraint

Abstract

The invention discloses a stratigraphic range constrained sequence modeling method, which relates to the technical field of three-dimensional geologic modeling, and comprises the following steps: preparing data including drilling data, regional geologic map data and a regional DEM, standardizing the drilling data including a standard stratigraphic table, a drilling basic information table and a drilling layered information table, the format of the form is unified as. Xlsx and a standard stratum table, the standard stratum table records the lithology of the whole stratum, and arrangement is carried out from top to bottom according to the new-to-old relationship; according to the sequence modeling method based on the stratum range constraint, the modeling boundary of each stratum is subjected to forced constraint, sequence modeling is carried out by using the stratum range constraint, the distribution condition of each stratum in a modeling area can be subjected to forced constraint, and meanwhile, modeling is carried out by combining drilling data and area elevation data, so that the model is finer, and the modeling efficiency is improved. And the error between a section cut through three-dimensional analysis of a sectioning model and the like and an original section is small.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of three-dimensional geological modeling, in particular to a stratum range constrained sequence modeling method. BACKGROUND

[0002] Three-dimensional geological modeling is a technology that converts two-dimensional geological data into a three-dimensional geological model using modern computer technology. It has been widely used in geology, oilfield exploitation and many other fields. Compared with traditional two-dimensional maps, converting two-dimensional geological data into a three-dimensional geological model is more conducive to data storage and management.

[0003] The layer-by-layer controlled sequence modeling method is suitable for areas where strata have iterative relationships and faults and structures are less developed, such as engineering geology Quaternary modeling and built-up area modeling. When the stratum relationship is complex and faults, folds and other structures are more developed, the layer-by-layer controlled sequence modeling method may face greater challenges, such as decreased model accuracy and increased computational complexity. Therefore, the applicability of this method under complex geological conditions still needs to be improved. At the same time, this method also has high requirements for the identification and division of stratum interfaces, and requires the participation and judgment of professional geologists, which to some extent limits its automation and intelligence. Therefore, we propose a stratum range constrained sequence modeling method. SUMMARY

[0004] The purpose of the present application is to provide a stratum range constrained sequence modeling method to solve the problems raised in the background art.

[0005] To achieve the above purpose, the present application provides the following technical solution: a stratum range constrained sequence modeling method, comprising the following steps:

[0006] S1: data preparation, including drilling data, regional geological map data and regional DEM;

[0007] S2: drilling data standardization processing, including standard stratum table, drilling basic information table and drilling layering information table, the table format is unified as.xlsx;

[0008] S201: standard stratum table, recording the lithology of the whole area stratum, and arranging from top to bottom according to the relationship from new to old;

[0009] S202: drilling basic information table, recording the coordinate position, hole depth and elevation of all drilling hole openings in the modeling area;

[0010] S203: drilling layering table, recording the stratum conditions corresponding to different buried depths in a single drilling hole, filling in the drilling number in the drilling layering table, extracting the stratum number corresponding to layering from the drilling data and filling it into the layering table, and calculating the layer bottom depth and layer bottom elevation according to the layer thickness;

[0011] S3: regional geological map data processing includes geological map standardization processing and stratum constraint range extraction;

[0012] S301: geological map standardization processing: extracting regional stratum surface data, and assigning attributes to each stratum surface;

[0013] S302: stratum constraint range extraction: after the geological map is standardized into SHP format, the stratum range of each stratum in the modeling area can be extracted layer by layer;

[0014] S4: regional DEM data, in modeling, the three-dimensional space expression of elevation data supports multi-field ground feature analysis and simulation, and adding DEM data in geological modeling restores the ground features of the modeling area;

[0015] S5: In the ZGIS Mine system, new projects are created, the above-mentioned standard stratum table, drilling basic information table, and drilling layer information table are imported into the engineering database, and field mapping matching is performed on each field. The SHP format data imported from the outside is used to control the boundary of each layer modeling;

[0016] S6: Roof and floor modeling is performed by selecting the modeling range, layer top DEM, downward depth / layer bottom DEM, and model accuracy to delineate the spatial range of modeling;

[0017] S7: Layer modeling is performed on the basis of roof and floor modeling, by selecting the interpolation method, selecting the controlled roof and floor layer model, checking the stratum range constraint, and adding regional drilling for control, single-layer layer model and body model can be constructed, and the full-area model can be constructed.

[0018] Preferably, in step S201, the standard stratum and color are extracted according to the stratum column chart of the geological map, and sorted and arranged.

[0019] Preferably, in step 202, the X, Y, and Z coordinates of the drilling hole and the hole depth data are filled according to the logging data, wherein the Z coordinate represents the hole elevation.

[0020] Preferably, in step 203, the layer top depth of the first layer of each drilling hole is 0m.

[0021] Preferably, in step 203, the calculation formula of the layer bottom depth and the layer bottom elevation is: layer bottom depth = layer top depth + layer thickness, and layer bottom elevation = layer top elevation - layer thickness.

[0022] Preferably, in step S302, professional geological knowledge is needed to assist in judgment before extraction.

[0023] Preferably, in the step S5, the drilling data checking tool is provided by the ZGIS Mine system, and the tool can be used for drilling data error checking after the data is imported.

[0024] Preferably, in the step S6, the spatial interpolation of the subsequent single stratum can only be performed in the roof and floor model.

[0025] Preferably, the main fields of the standard stratum table include stratum code, lithology description and color.

[0026] Preferably, the elevation of the bottom of the last layer of the drilling hole needs to be consistent with the hole depth in the drilling hole basic information table.

[0027] Compared with the prior art, the method has the beneficial effects that:

[0028] According to the method, the modeling boundary of each stratum is forcedly constrained, the stratum range constraint is used for sequence modeling, the distribution of each stratum in the modeling area can be forcedly constrained, the modeling is performed in combination with the drilling data and the regional elevation data, the model is more fine and more in line with the actual situation, and the error between the profile cut by the model three-dimensional analysis and the original profile is small. BRIEF DESCRIPTION OF DRAWINGS

[0029] Figure 1 It is a whole structure schematic diagram of the present application;

[0030] Figure 2 It is a regional geological map standardization processing schematic diagram of the present application;

[0031] Figure 3 It is a stratum constraint extraction schematic diagram of the present application;

[0032] Figure 4 It is a DEM schematic diagram of the present application;

[0033] Figure 5 It is a data import schematic diagram of the present application;

[0034] Figure 6 It is a drilling hole checking schematic diagram of the present application;

[0035] Figure 7 It is an added boundary constraint schematic diagram of the present application;

[0036] Figure 8 It is a roof and floor modeling schematic diagram of the present application;

[0037] Figure 9 It is a roof and floor model schematic diagram of the present application;

[0038] Figure 10 It is a layer surface modeling schematic diagram of the present application;

[0039] Figure 11 Fig. 1 is a schematic diagram of a geological body model according to the present application. DETAILED DESCRIPTION

[0040] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by a person of ordinary skill in the art without creative work fall within the protection scope of the present application.

[0041] Embodiment: please refer to Figures 1 to 11 The present application provides a technical solution: a stratum range constrained sequence modeling method, comprising the following steps:

[0042] S1: data preparation, the data including drilling data, regional geological map data and regional DEM;

[0043] S2: drilling data standardization processing, the standardization processing including a standard stratum table, a drilling basic information table and a drilling layering information table, the table format being unified as.xlsx;

[0044] S201: standard stratum table, the standard stratum table recording the lithology of the whole region, and being arranged from top to bottom according to the relationship from new to old;

[0045] S202: drilling basic information table, recording the coordinate position, hole depth and elevation of all drilling hole openings in the modeling area;

[0046] S203: drilling layering table, recording the stratum condition corresponding to different buried depths in a single drilling hole, filling in the drilling number in the drilling layering table, extracting the stratum number corresponding to layering from the drilling data, filling in the layering table, and calculating the layer bottom depth and layer bottom elevation according to the layer thickness;

[0047] S3: regional geological map data processing including geological map standardization processing and stratum constraint range extraction;

[0048] S301: geological map standardization processing: extracting the stratum surface data of the region, and giving attributes to each stratum surface;

[0049] S302: stratum constraint range extraction: after the geological map is standardized into SHP format, the stratum range of each stratum in the modeling area can be extracted layer by layer;

[0050] S4: regional DEM data, the three-dimensional space expression of elevation data supports multi-field ground feature analysis and simulation in modeling, and adding DEM data in geological modeling restores the ground features of the modeling area;

[0051] S5: In the ZGIS Mine system, a new project is established, the above-mentioned standard stratum table, drilling basic information table and drilling stratified information table are imported into the project database, and field mapping matching is performed on each field, and the SHP format data imported from the outside is used to control the boundary of modeling of each layer;

[0052] S6: The roof and floor modeling is performed by selecting the modeling range, layer top DEM, lower extension depth / layer bottom DEM and model precision to define the spatial range of modeling.

[0053] S7: The layer surface modeling is performed on the basis of the roof and floor modeling, the interpolation method is selected, the controlled roof and floor layer surface model is selected, the stratum range constraint is checked, the regional drilling is added for control, and the single layer surface model and body model can be constructed and the full-area model can be constructed.

[0054] As a further limitation of the present application,

[0055] In step S201, the standard stratum and color are extracted according to the stratum column chart of the geological map, and are sorted and arranged.

[0056] In step 202, the X, Y and Z coordinates of the drilling and the hole depth data are filled according to the logging data, wherein the Z coordinate represents the hole opening elevation.

[0057] In step 203, the layer top depth of the first layer of each drilling is 0 m.

[0058] In step 203, the calculation formula of the layer bottom depth and the layer bottom elevation is: layer bottom depth = layer top depth + layer thickness, and layer bottom elevation = layer top elevation - layer thickness.

[0059] In step S302, the professional geological knowledge of the geological personnel needs to be added for auxiliary judgment and extraction.

[0060] In step S5, the drilling data checking tool provided by the ZGIS Mine system is used to check the drilling data error after the data is imported.

[0061] In step S6, the spatial interpolation of the subsequent single stratum can only be performed in this roof and floor model.

[0062] The main fields referred to in the standard stratum table include stratum code, lithology description and color.

[0063] The layer bottom elevation of the last layer of the drilling needs to be consistent with the hole depth in the drilling basic information table.

[0064] The specific implementation of the embodiment is:

[0065] The original drilling data usually comes from the mine logging data, which cannot be directly used due to the difference in the data field content used by each modeling software for modeling. After standardization processing, it can be imported into the database of the modeling software to participate in subsequent modeling. The standardization processing mainly includes the following three tables. All table formats are unified as ".xlsx".

[0066] Standard stratigraphic table: The standard stratigraphic table records the lithology of the strata in the whole area. From top to bottom, it is arranged according to the relationship from new to old. The standard strata and colors can be extracted from the stratigraphic column of the geological map and sorted. The following is the format reference for the standard stratigraphic table. The main fields include stratigraphic code, lithology description, color, etc. as shown in the following table:

[0067] Standard stratigraphic table

[0068] Field name Type Explanation Stratum number Text Rock name Text Stratum color Text Hexadecimal code Geological description Text ……

[0069] Drilling basic information table: The drilling basic information table records the coordinate position, hole depth, elevation, etc. of all drilling hole openings in the modeling area, which is convenient for the system to identify drilling information. When arranging, the X, Y, Z coordinates (Z coordinate is the hole elevation) and hole depth data of the drilling should be filled according to the logging data. The arrangement format is as follows:

[0070] Drilling basic information table

[0071] Field name Type Explanation Hole number Text Orifice elevation (m) Number Hole depth (m) Number X (m) Number Y (m) Number ……

[0072] Drilling stratification information table: The drilling stratification table records the strata corresponding to different buried depths in a single drilling hole. The drilling stratification information table, drilling basic information table and standard stratigraphic table are in a correlated relationship and complement each other. The drilling number is filled in the drilling stratification table. The stratigraphic number corresponding to the stratification is extracted from the logging data and filled in the table. The top depth of the first layer of each drilling hole is 0m. The bottom depth and bottom elevation can be calculated according to the layer thickness (bottom depth = top depth + layer thickness; bottom elevation = top elevation - layer thickness). The following table shows the arrangement format:

[0073] Drilling stratification table

[0074] Field name Type Explanation Hole number Text Stratum number Text Depth of layer top Number From Depth of layer bottom Number To Layer thickness Number Geological description Text ……

[0075] Note: The bottom elevation of the last layer of the drilling hole should be consistent with the hole depth in the drilling basic information table.

[0076] Geological map standardization processing: The original geological map data is large and rich, containing many data that are not needed for modeling. The geological map data needs to be simplified to extract the data that can be used for this modeling. The main task is to extract the stratum surface data of the area and assign attributes to each stratum surface, such as shown in the following table. Figure 2 ​

[0077] Stratigraphic constraint range extraction: The highlight of this method is to add stratigraphic range to impose constraints on the modeling boundary. After standardizing the geological map into SHP format, the stratigraphic range of each stratigraphic layer in the modeling area can be extracted by stratification. In this process, geologists need to add professional geological knowledge to assist in judgment and extraction. It has certain requirements for the professional technical knowledge of modeling personnel, such as Figure 3 .

[0078] DEM (Digital Elevation Model) supports multi-field surface feature analysis and simulation through three-dimensional spatial expression of elevation data in modeling. Adding DEM data in geological modeling can accurately restore the surface features of the modeling area. DEM can be generated by interpolation method according to the terrain line of the region, or open source DEM data can be downloaded through 91 map, geological cloud and other software, such as Figure 4 .

[0079] After all the data preparation is completed, the construction of the sequence model can be started. In the ZGIS Mine system, a new project is created, and the above-mentioned standard stratigraphic table, drilling basic information table, and drilling stratification information table are imported into the project database. The fields are matched by field mapping, such as Figure 5 .

[0080] The correctness of the modeling data directly affects the modeling progress and model effect. Therefore, data quality inspection before modeling is necessary. The system provides a drilling data checking tool. After importing the data, the tool can be used to check the drilling data errors, such as Figure 6 .

[0081] Adding range constraints can achieve the purpose of accurately controlling the modeling boundary of each layer. In the standard stratigraphic management function, the method of adding boundary constraints can be selected through advanced settings. Generally, SHP format data is imported from the outside for constraint, such as Figure 7 .

[0082] The top and bottom plate modeling circumscribes the spatial range of modeling by selecting the modeling range, layer top DEM, downward depth / layer bottom DEM, and model precision. Subsequent spatial interpolation of individual stratigraphic layers can only be performed within this top and bottom plate model. It is the basis of sequence modeling, such as Figure 8 and Figure 9 .

[0083] Layer modeling is based on top and bottom plate modeling. By selecting the appropriate interpolation method, selecting the control top and bottom plate layer model, checking the stratigraphic range constraints, and adding regional drilling for control, single layer layer model and body model can be constructed. Full-area model can also be quickly constructed, such as Figure 10 and Figure 11 .

[0084] It has to be noted that, in the present document, the terms "first", "second", etc. merely serve to identify a subject or action, without necessarily requiring or implying any such actual relationship or order between such subjects or actions. Moreover, the terms "comprises", "comprising", or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but can include other elements not expressly listed or inherent to such process, method, article, or apparatus.

[0085] While embodiments of the present application have been shown and described, it is to be understood that various modifications, substitutions, combinations, and variations of the embodiments can be undertaken without departing from the spirit and scope of the present application, which is defined by the appended claims and their equivalents.

Claims

1. A method of stratigraphically-constrained sequence modeling, characterized by: The method comprises the following steps: S1: data preparation, the data including drilling data, regional geological map data and regional DEM; S2: drilling data standardization processing, the standardization processing including a standard stratum table, a drilling basic information table and a drilling layering information table, and the table format being unified as.xlsx; S201: the standard stratum table, the standard stratum table recording the lithology of the whole region, and being arranged from top to bottom according to the relationship from new to old; S202: the drilling basic information table, recording the coordinate position, hole depth and elevation of all drilling orifices in the modeling area; S203: the drilling layering table, recording the stratum condition corresponding to different buried depths in a single drilling, filling in the drilling number in the drilling layering table, extracting the stratum number corresponding to layering from the drilling data, filling in the layering table, and the layer bottom depth and layer bottom elevation being calculated according to the layer thickness; S3: regional geological map data processing including geological map standardization processing and stratum constraint range extraction; S301: the geological map standardization processing: extracting the stratum surface data of the region, and giving attributes to each stratum surface; S302: the stratum constraint range extraction: after the geological map is standardized into the SHP format, the stratum range of each stratum in the modeling area can be extracted layer by layer; S4: regional DEM data, the three-dimensional space expression of elevation data supports multi-field ground feature analysis and simulation in modeling, and adding DEM data in geological modeling restores the ground features of the modeling area; S5: in the ZGIS Mine system, a new project is newly built, the above-mentioned standard stratum table, drilling basic information table and drilling layering information table are imported into the project database, and field mapping matching is performed on each field, and the SHP format data imported from outside is used to control the boundary of each layer modeling; S6: the roof and floor modeling is performed by selecting the modeling range, layer top DEM, downward depth / layer bottom DEM and model precision, and the spatial range of modeling is circled; S7: the layer surface modeling is performed on the basis of the roof and floor modeling, the interpolation method is selected, the controlled roof and floor layer surface model is selected, the stratum range constraint is checked, the regional drilling is added for control, and the single layer surface model and body model can be constructed, and the whole area model can be constructed.

2. The method of claim 1, wherein: In the step S201, the standard stratum and color are extracted according to the stratum column chart of the geological map, and are sorted and arranged.

3. The method of claim 1, wherein: In the step 202, the X, Y and Z coordinates of the drilling and the hole depth data are filled according to the logging data, wherein the Z coordinate represents the orifice elevation.

4. The method of claim 1, wherein: In the step 203, the layer top depth of the first layer of each drilling is 0 m.

5. The method of claim 1, wherein: In the step 203, the calculation formula of the layer bottom depth and layer bottom elevation is: layer bottom depth = layer top depth + layer thickness, and layer bottom elevation = layer top elevation - layer thickness.

6. The method of claim 1, wherein: In the step S302, the professional geological knowledge of the geological personnel needs to be added for auxiliary judgment before extraction.

7. The method of claim 1, wherein: In the step S5, the drilling data checking tool is provided by the ZGIS Mine system, and the drilling data error checking can be performed by using the tool after the data is imported.

8. The method of claim 1, wherein: In the step S6, the spatial interpolation of a single stratum can only be performed in the roof and floor model.

9. The method of claim 1, wherein: The main fields referenced in the standard stratigraphic table include stratigraphic code, lithology description and color.

10. The method of claim 1, wherein: The elevation of the bottom of the last layer of the borehole should be consistent with the hole depth in the borehole basic information table.

Images