Method and system for automatically generating stratigraphic partition map

Abstract

The invention discloses an automatic generation method and system for a stratigraphic partition map. The method comprises the steps that a unified stratigraphic sequence is established, and all drill hole layers are unified; acquiring distribution of all sampling points and boundary points of each stratum, and determining whether the attribute information of each point is used for marking whether the current stratum is missing or not; carrying out triangulation on the distribution of the sampling points and the boundary points of each stratum, judging whether a current triangle needs to be extinguished or not according to attribute information of three points forming the triangle, and if the current triangle needs to be extinguished, carrying out extinguishment position calculation to determine a extinguishment position point; and adding the pinching position points to the sampling point and boundary point distribution of the current stratum, assigning the attributes of the pinching position points according to the condition that the current stratum is not missing, deleting the points with the missing current stratum in the sampling point and boundary point distribution, and carrying out triangulation on the left points again. The invention provides a set of practical stratum pinching position algorithm, and the algorithm is used for tracking the stratum subareas, so that the construction precision of the three-dimensional geologic model can be improved.

Description

Technical Field

[0001] The invention relates to the field of three-dimensional geological modeling, and in particular to a method and system for automatically generating a stratigraphic zoning map. Background Art

[0002] The three-dimensional geological model is a simulation of the complex geological structure of the urban underground space. Through the geological model, the structural structure of the urban underground space can be clearly seen. The three-dimensional geological model plays an increasingly important and in-depth role in the current geological needs and applications. Three-dimensional geological modeling is one of the key technologies for geological data visualization, geological space analysis and glass earth construction. Its main idea is to process a series of geological data such as geological drilling, geological profiles, seismic profiles, faults, lenses, etc. through computer technology, and use three-dimensional visualization technology to express geological bodies, geological phenomena and geological processes in three-dimensional space, so as to assist geological personnel in visual analysis and decision-making. Three-dimensional geological models are widely used in the field of geology in advanced geological prediction, subway line selection, ore reserve calculation, etc., helping workers in related fields to improve work efficiency and conduct in-depth research.

[0003] Geological environment evaluation, deep geological prospecting, and urban underground space development and safe utilization all require accurate exploration of underground geological structures as support. Three-dimensional geological modeling is an important means to simulate geological structural characteristics and geological environment evolution. Three-dimensional geological models play a very important role in underground space development. High-precision models can guide urban underground space development and calculate capital investment. Using low-precision models to guide development may cause safety problems or increase unnecessary capital investment. In actual needs, people have higher and higher requirements for the accuracy of models. Therefore, we need to conduct further in-depth research on maintaining the consistency between models and data and the correctness of model trends to ensure that three-dimensional geological models can provide correct guidance and support for people's production and life. In short, three-dimensional geological models can provide many conveniences for people's production and life through simulation of urban underground space. As a powerful weapon to guide the development of underground space, the three-dimensional geological model has higher and higher requirements for its accuracy. Constructing high-precision three-dimensional geological models has become a key content in geological informatization.

[0004] The pinch-out position of the stratum determines the distribution range of the stratum, and thus determines the accuracy of the 3D geological model. The traditional stratum pinch-out algorithm does not consider the factors that affect the pinch-out position. The use of a fixed pinch-out position makes the distribution range of the stratum inaccurate, often resulting in a large difference between the stratum distribution of the model and the actual one. Summary of the invention

[0005] The present application provides a method and system for automatically generating a stratigraphic zoning map, in order to solve the problem that a traditional stratigraphic pinch-out algorithm does not consider factors that affect the pinch-out position, uses a fixed pinch-out position, makes the distribution range of the stratigraphic area inaccurate, and causes a large difference between the stratigraphic distribution of the three-dimensional model and the actual situation.

[0006] According to the first aspect, an embodiment provides a method for automatically generating a stratigraphic zoning map, the method comprising:

[0007] Step S1, obtaining multiple drilling data of the target area and preprocessing and extracting information from the drilling data;

[0008] Step S2, establishing a unified stratigraphic sequence, dividing the boreholes into stratigraphic layers based on the unified stratigraphic sequence, making the layers of each borehole consistent with the layers of the unified stratigraphic sequence, and unifying the layers of all boreholes;

[0009] Step S3, taking each drilling position as a sampling point, extracting all sampling points of each stratum after division, taking into account the regional boundary constraints, obtaining all sampling points and boundary point distribution of each stratum, and determining the attribute information of each point, the attribute information of each point is used to mark whether the current stratum is missing;

[0010] Step S4, triangulate the sampling points and boundary points of each stratum to obtain an irregular triangulated network, traverse each triangle, and determine whether the current triangle needs to be pinched out according to the attribute information of the three points that make up the triangle. If pinching out is required, the pinching out position is calculated to determine the pinching out position point;

[0011] Step S5, adding the pinch-out location point to the sampling point and boundary point distribution of the current stratum, and assigning the attributes of the pinch-out location point as if there is no missing current stratum, deleting the points in the sampling point and boundary point distribution that are missing from the current stratum according to the attribute information of each point, re-triangulating the remaining points and obtaining the stratigraphic plane partition map of the current stratum.

[0012] Furthermore, the step S1 specifically includes:

[0013] Collect original drilling data and profile data, and convert profile data into drilling data through data extraction;

[0014] Check and modify the obtained drilling data, and screen out the drilling holes with missing formations, formation inversions and formation duplications;

[0015] Information is extracted based on the collected drilling data, and the extracted information is used for statistics and analysis to obtain the parameter information required for pinch-out position calculation.

[0016] Furthermore, the step S2 specifically includes:

[0017] Within the stratum, according to the preset stratigraphic division order, combined with the original stratification information of the boreholes within the area, the geology and the overall geological structure characteristics of the region, the strata are subdivided into similar levels to establish a unified stratigraphic sequence.

[0018] Furthermore, the step S2 specifically includes:

[0019] After the unified stratigraphic sequence is determined, for the original borehole strata that are inconsistent with the unified stratigraphic sequence, the corresponding zero-thickness layer is inserted to make the strata of each borehole consistent with the strata of the unified stratigraphic sequence.

[0020] Furthermore, the step S3 specifically includes:

[0021] For each stratum, each sampling point may contain only two situations: containing the corresponding stratum or not containing the corresponding stratum. The attribute value of the borehole sampling point that does not contain the corresponding stratum, i.e., the stratum is missing, is assigned a value of 1; the attribute value of the borehole sampling point that contains the corresponding stratum, i.e., the stratum is not missing, is assigned a value of 0;

[0022] The attribute value of the boundary point is determined by interpolating the sampling points of the existing attribute value, or by finding the nearest sampling point in the vicinity and taking the attribute value of the nearest sampling point in the vicinity as the attribute value of the current boundary point.

[0023] Furthermore, the step S4 specifically includes:

[0024] If the attribute values ​​of the three points that make up the triangle are all 0 or all 1, there is no need to determine the pinch-out position;

[0025] If the attributes of the three points that make up the triangle contain both 1 and 0, then it is necessary to perform a pinch-out process, and there must be two edges in the triangle where the pinch-out points appear.

[0026] Furthermore, the step S4 specifically includes:

[0027] For the points that need to be pinched out, the pinch-out coefficient on the boundary of the existing strata is calculated according to the pinch-out coefficient calculation method of the angular unconformity strata, and the pinch-out point position is calculated according to the pinch-out coefficient.

[0028] Furthermore, the step S5 specifically includes:

[0029] Assign the attribute value of the pinch-out position point to 1, filter out the points with an attribute value of 0, and leave the set of points with an attribute value of 1. Use the triangulation algorithm again to triangulate the points with an attribute value of 1 to obtain an irregular triangulated network, which is the stratigraphic partition of the corresponding stratigraphic layer.

[0030] The boundary extracted from the irregular triangulated network is the outer boundary of the stratigraphic partition; according to the stratigraphic partition boundary conditions, the stratigraphic partition boundary is smoothed using a smoothing algorithm.

[0031] Furthermore, the method further comprises:

[0032] The triangulation is a triangulation algorithm that combines the principle of Delaunay triangulation and triangulates adjacent rings of polygons from outside to inside based on the relationship between polygon vertices, while taking "island polygons" into consideration.

[0033] According to the second aspect, an embodiment provides a system for automatically generating a stratigraphic zoning map, the system comprising:

[0034] A data acquisition module is used to obtain multiple borehole data in the target area and perform preprocessing and information extraction on the borehole data;

[0035] A stratigraphic division module is used to establish a unified stratigraphic sequence, divide the boreholes into stratigraphic layers based on the unified stratigraphic sequence, make the layers of each borehole consistent with the layers of the unified stratigraphic sequence, and unify the layers of all boreholes;

[0036] The attribute determination module is used to take each drilling position as a sampling point, extract all sampling points of each stratum after division, and consider the regional boundary constraints at the same time, obtain all sampling points and boundary point distribution of each stratum, and determine the attribute information of each point. The attribute information of each point is used to mark whether there is a missing current stratum;

[0037] The pinch-out position determination module is used to triangulate the sampling points and boundary points of each stratum to obtain an irregular triangulated network, traverse each triangle, and determine whether the current triangle needs to be pinched out based on the attribute information of the three points that make up the triangle. If pinch-out is required, the pinch-out position calculation is performed to determine the pinch-out position point;

[0038] The stratigraphic zoning map generation module is used to add the pinch-out location points to the sampling points and boundary point distribution of the current stratigraphic layer, and assign the attributes of the pinch-out location points as if there is no missing current stratigraphic layer, delete the points in the sampling points and boundary point distribution that are missing from the current stratigraphic layer according to the attribute information of each point, re-triangulate the remaining points and obtain the stratigraphic plane zoning map of the current stratigraphic layer.

[0039] The present application provides a method and system for automatically generating a stratigraphic zoning map, which has the following beneficial effects:

[0040] (1) The traditional formation pinch-out algorithm does not consider the factors that affect the pinch-out position. The use of a fixed pinch-out position makes the distribution range of the formation inaccurate, which often leads to a large difference between the formation distribution of the model and the actual situation. The present invention fully analyzes the factors that affect the formation pinch-out position, based on a formation pinch-out position algorithm that conforms to reality, considers factors such as formation thickness, borehole spacing, and single-hole control range, realizes reasonable pinch-out position calculation, and tracks formation partitions based on the algorithm, which can help improve the construction accuracy of the three-dimensional geological model;

[0041] (2) The present invention unifies the stratigraphic sequence. For the original borehole strata that are inconsistent with the unified stratigraphic sequence, the corresponding zero-thickness layer is inserted, so that the strata of each borehole are consistent with the strata of the unified stratigraphic sequence. Finally, each borehole is divided into 38 strata, ensuring that the stratigraphic division adopts a unified standard and the data of each stratum is accurate.

[0042] (3) The present invention assigns a value to the attributes of each stratum sampling point according to whether there is a current stratum missing, and considers boundary constraints at the same time. Through triangulation, the attribute information of the three points that make up the mesh triangle is used to determine whether each mesh triangle needs to be pinched out. The pinch-out position judgment based on the attribute information can achieve fast and accurate pinch-out position processing and calculation. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] Figure 1 A flow chart of a method for automatically generating a stratigraphic zoning map provided by one embodiment of the present invention;

[0044] Figure 2 A specific implementation flow chart of a method for automatically generating a stratigraphic zoning map provided by an embodiment of the present invention;

[0045] Figure 3 A schematic diagram of unified stratigraphic sequence coding in a method for automatically generating a stratigraphic zoning map provided by one embodiment of the present invention;

[0046] Figure 4 A schematic diagram of the distribution of stratum sampling points in a method for automatically generating a stratum zoning map provided by an embodiment of the present invention;

[0047] Figure 5 A schematic diagram of boundary point distribution in a method for automatically generating a stratigraphic zoning map provided by an embodiment of the present invention;

[0048] Figure 6 A schematic diagram of triangulation results in a method for automatically generating a stratigraphic zoning map provided by an embodiment of the present invention;

[0049] Figure 7 A schematic diagram of determining the direction of formation pinch-out in a method for automatically generating a formation partition map provided by an embodiment of the present invention;

[0050] Figure 8 A stratigraphic zoning map in a method for automatically generating a stratigraphic zoning map provided by one embodiment of the present invention;

[0051] Fig. 9 A schematic diagram of the logical structure of a system for automatically generating a stratigraphic zoning map provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0052] The present invention is further described in detail below by specific embodiments in conjunction with the accompanying drawings. Wherein similar elements in different embodiments adopt associated similar element numbers. In the following embodiments, many detailed descriptions are for making the present application better understood. However, those skilled in the art can easily recognize that some features can be omitted in different situations, or can be replaced by other elements, materials, methods. In some cases, some operations related to the present application are not shown or described in the specification, this is to avoid the core part of the present application being overwhelmed by too much description, and for those skilled in the art, it is not necessary to describe these related operations in detail, and they can fully understand the related operations according to the description in the specification and the general technical knowledge in the art.

[0053] In addition, the features, operations or characteristics described in the specification can be combined in any appropriate manner to form various implementations. At the same time, the steps or actions in the method description can also be interchanged or adjusted in a manner that is obvious to those skilled in the art. Therefore, the various sequences in the specification and the drawings are only for the purpose of clearly describing a certain embodiment and are not meant to be a required sequence, unless otherwise specified that a certain sequence must be followed.

[0054] The first embodiment of the present invention provides a method for automatically generating a stratigraphic zoning map. Figure 1 and Figure 2 Provide detailed explanation.

[0055] like Figure 1 As shown, in step S1, multiple drilling data of the target area are acquired and the drilling data are preprocessed and information extracted.

[0056] The above steps specifically include:

[0057] S11. Preprocessing of drilling data includes checking and modifying drilling data. The drilling data is checked by software to ensure the data format is standardized, and the stratification information of the drilling holes is checked and screened to extract the drilling holes with missing strata, strata inversion and strata duplication. For the information in the inspection results that cannot be modified by the software, manual modification is performed.

[0058] S12. Mark the missing strata of each borehole, and detect the missing strata information according to the strata missing discrimination rules; combine the three-dimensional geological modeling idea of ​​"first large layer, then sub-layer, and then sub-layer again", and mark the missing strata of each borehole step by step and graded.

[0059] S13. Utilize the borehole database to extract information that can be used for calculations, such as stratum thickness, borehole spacing, and stratum burial depth, and calculate and analyze the average thickness of each stratum, the median stratum thickness, the average controlled area of ​​the borehole, and other parameter information that can be used for pinch-out position calculations.

[0060] The main data used in this embodiment are drilling data and profile data. The drilling data mainly comes from engineering survey reports such as foundation, roadbed, pipeline, etc., and a total of 186 boreholes in Beihai City have been collected. The profile data includes 3 cross sections and 3 longitudinal sections, a total of 6 profile data, which are used to constrain stratum modeling.

[0061] The drilling data was checked using Smart3dMap modeling software and found that there were 174 boreholes with missing strata, 32 boreholes with inverted strata, and 57 boreholes with repeated strata. The profile data was automatically processed using Smart3dMap modeling software, and a borehole layer data was extracted every 1 meter from the profile, and the profile data was converted into drilling data for participation.

[0062] like Figure 1 As shown, in step S2, a unified stratigraphic sequence is established, and the boreholes are stratigraphically divided based on the unified stratigraphic sequence, so that the stratification of each borehole is consistent with the stratification of the unified stratigraphic sequence, and the stratification of all boreholes is unified.

[0063] In this embodiment, within the stratum, according to the order of the stratum system, group, and lithology, combined with the original layering information of the boreholes within the modeling range, the geology, and the overall geological structure characteristics of the region, the stratum near level is subdivided to establish a standard stratigraphic sequence table, as shown in Table 1. When connecting strata, when there is a missing stratum, a stratum duplication, or a stratum inversion, the insertion position of the zero thickness layer is determined by referring to the standard stratigraphic sequence table.

[0064] Table 1 Part of the Quaternary standard stratigraphic sequence in the region

[0065]

[0066]

[0067] Specifically, Figure 3, the unified stratigraphic sequence method was used to analyze and count the borehole strata within the modeling range, and a total of 38 unified stratigraphic sequences were established using the unified stratigraphic sequence coding method (the original standard stratigraphic sequences were 21). After the unified stratigraphic sequence was determined, for the original borehole strata that were inconsistent with the unified stratigraphic sequence, the corresponding zero-thickness layer was inserted to make the stratification of each borehole consistent with the stratification of the unified stratigraphic sequence, that is, each borehole was divided into 38 strata (missing strata were represented by zero-thickness layers).

[0068] like Figure 1 As shown, in step S3, each drilling position is taken as a sampling point, and all sampling points of each divided stratum are extracted. At the same time, the regional boundary constraints are considered, and the distribution of all sampling points and boundary points of each stratum is obtained, and the attribute information of each point is determined. The attribute information of each point is used to mark whether there is a missing current stratum.

[0069] Specifically, according to the above content, each borehole is divided into 38 strata, and all the borehole points belonging to the same stratum are extracted, and each borehole position is used as a sampling point, such as Figure 4 As shown. The sampling points are marked according to whether they contain the current stratum. Each sampling point can only contain the stratum or not contain the stratum. The missing boreholes are assigned an attribute value of 0, and the non-missing boreholes are assigned an attribute value of 1.

[0070] Since stratigraphic partitioning is used to constrain modeling, the modeling boundary points also need to be used as sampling points for stratigraphic partitioning calculations, so the attributes of the modeling boundary points need to be determined. The attribute values ​​of the boundary points are determined by interpolating the sampling points with existing attribute values ​​or by finding the nearest sampling points around them and taking the attribute values ​​of the nearest sampling points around them as the attribute values ​​of the current boundary points, such as Figure 5 shown.

[0071] like Figure 1 As shown, in step S4, the sampling points and boundary points of each stratum are triangulated to obtain an irregular triangulated network, and each triangle is traversed to determine whether the current triangle needs to be pinched out based on the attribute information of the three points that make up the triangle. If pinching out is required, the pinching out position is calculated to determine the pinching out position point.

[0072] Specifically, all sampling points are triangulated to obtain an irregular triangulated network, such as Figure 6 For strata with missing parts, loop through each triangle and choose whether to pinch out based on the attribute values ​​of the three points of the current triangle. If the attribute values ​​of the three points are all 0 or all 1, there is no need to determine the pinch-out position. If the attributes of the three points are both 1 and 0, pinch-out processing is required, and there must be two edges in the triangle with pinch-out points.

[0073] In this embodiment, for the points that need to be pinched out, the pinch-out coefficient on the boundary of the existing strata is calculated according to the pinch-out coefficient calculation method of the angular unconformity strata, and the pinch-out point position is accurately determined according to the pinch-out coefficient.

[0074] If pinch-out is required, the pinch-out algorithm for unconformity strata is used, taking into account the constraints of the plane boundary. Calculate the average drilling distance D in the modeling area h By comparing the actual distance between the drill holes to be pinched out (drill holes without missing strata at a certain stratum level) and the pinched-out drill holes (with missing strata at a certain stratum level) and the average drill hole distance in the modeling area, the pinch-out calculation method of the stratum to be pinched out is determined. The specific process of calculating the pinch-out position is as follows:

[0075] a. The strata pinch out step by step in the order of "large layer first, then sublayer, and then sublayer again", taking into account the constraints of the top and bottom DEM interfaces of the model, the constraints of the modeling plane boundaries, and the constraints of the bottom DEM surfaces of each large layer.

[0076] b. Calculate the average drilling distance D in the modeling area h Assuming the area of ​​the modeling range is A and the number of holes is N, the control area of ​​a single hole is:

[0077]

[0078] The diameter of a single drill hole (applicable to the case where the drill holes are evenly distributed), that is, the average distance between two drill holes:

[0079]

[0080] c. Select the drill holes to be pinched out (drill holes with no missing strata at a certain level) and the pinched-out drill holes (drill holes with missing strata at a certain level), and determine the pinch-out calculation method for the formations to be pinched out by comparing the actual distance and average distance between the two holes.

[0081] d. If the actual distance between the two holes is D ac Greater than the average distance D h , then the pinch-out point is selected at a distance D from the drill hole to be pinched out h Position. If the actual distance between the two holes is D ac Less than or equal to the average distance D h , then go to step e.

[0082] e. Set the pinch-out coefficient PinchCoef (pc), and the minimum thickness of the pinch-out stratum in the modeling area is t min , the maximum thickness is t max , the average layer thickness is thick ave (t ave ), the median layer thickness is t m , compare t ave and tm , according to the size of the two, we get two values ​​t a and t b (t a <t b ). Assuming that the current formation pinch-out thickness is t, the pinch-out coefficient calculation is divided into two cases:

[0083] like The value range of pc is For example, when n=3,

[0084] If t a ≤t≤t max ,but The value range of pc is

[0085] The corresponding positions to be pinched out in the above two cases are respectively Place.

[0086] When all triangles are traversed, the pinch-out positions are calculated for all places with missing strata within the area to be partitioned.

[0087] Subsequent modeling will involve the judgment of the pinch-out direction. i , the method for determining the direction of pinch-out is as follows:

[0088] There are strata D1…D from top to bottom in the adjacent drill holes to be pinched out. i-1 , D i , D i+1 …D n , if i = 1, the pinch-out direction is the direction of the formation D in the hole to be pinched out. i The base plate pinches out upward to the calculated pinch-out position, e.g. Figure 7 a; if i = n, the pinch-out direction is the direction of the formation D in the hole to be pinched out. i The top plate pinches out downward to the calculated pinch-out position, such as Figure 7 b; if i≠1 and i≠n, the pinch-out direction is the direction of the formation D in the hole to be pinched out. i The top plate and bottom plate are pushed flat and intersect at the calculated pinch-out position, such as Figure 7 c.

[0089] like Figure 1 As shown, in step S5, the pinch-out location point is added to the sampling point and boundary point distribution of the current stratum, and the attributes of the pinch-out location point are assigned according to the absence of the current stratum missing, and the points in the sampling point and boundary point distribution that are missing the current stratum are deleted according to the attribute information of each point, and the remaining points are re-triangulated to obtain the stratigraphic plane zoning map of the current stratum.

[0090] Specifically, add the pinch-out position to the sampling point set, and assign the attribute value of the pinch-out position point to 1. Filter the sampling points with an attribute value of 0, leaving the sampling point set with an attribute value of 1, and use the triangulation algorithm again to triangulate the sampling points with an attribute value of 1. The irregular triangulated network obtained by triangulation is the stratigraphic partition of the stratum, such as Figure 8 As shown in the figure, the boundary extracted from the irregular triangulated network is the outer boundary of the stratigraphic partition. According to the boundary of the stratigraphic partition, the stratigraphic partition is smoothed using a smoothing algorithm.

[0091] The triangulation algorithm described in this embodiment is a triangulation algorithm that combines the principle of Delaunay triangulation and triangulates adjacent rings of polygons from the outside to the inside based on the relationship between polygon vertices, while taking into account the triangulation algorithm of "island polygons".

[0092] An automatic generation method of a stratigraphic zoning map provided in an embodiment of the present invention fully analyzes the factors affecting the stratigraphic pinch-out position, and from the data level, considers factors such as stratigraphic thickness, borehole spacing, and single-hole control range, proposes a set of stratigraphic pinch-out position algorithms that conform to reality, and uses the algorithm to track stratigraphic zoning, thereby improving the accuracy of model construction.

[0093] Corresponding to the above disclosed method for automatically generating a stratigraphic zoning map, the embodiment of the present invention further discloses a system for automatically generating a stratigraphic zoning map, such as Fig. 9 As shown, it specifically includes:

[0094] A data acquisition module is used to obtain multiple borehole data in the target area and perform preprocessing and information extraction on the borehole data;

[0095] A stratigraphic division module is used to establish a unified stratigraphic sequence, divide the boreholes into stratigraphic layers based on the unified stratigraphic sequence, make the layers of each borehole consistent with the layers of the unified stratigraphic sequence, and unify the layers of all boreholes;

[0096] The attribute determination module is used to take each drilling position as a sampling point, extract all sampling points of each stratum after division, and consider the regional boundary constraints at the same time, obtain all sampling points and boundary point distribution of each stratum, and determine the attribute information of each point. The attribute information of each point is used to mark whether there is a missing current stratum;

[0097] The pinch-out position determination module is used to triangulate the sampling points and boundary points of each stratum to obtain an irregular triangulated network, traverse each triangle, and determine whether the current triangle needs to be pinched out based on the attribute information of the three points that make up the triangle. If pinch-out is required, the pinch-out position calculation is performed to determine the pinch-out position point;

[0098] The stratigraphic zoning map generation module is used to add the pinch-out location points to the sampling points and boundary point distribution of the current stratigraphic layer, and assign the attributes of the pinch-out location points as if there is no missing current stratigraphic layer, delete the points in the sampling points and boundary point distribution that are missing from the current stratigraphic layer according to the attribute information of each point, re-triangulate the remaining points and obtain the stratigraphic plane zoning map of the current stratigraphic layer.

[0099] It should be noted that, for the detailed description of the automatic generation system of a stratigraphic zoning map provided in an embodiment of the present invention, reference can be made to the relevant description of the automatic generation method of a stratigraphic zoning map provided in an embodiment of the present application, which will not be repeated here.

[0100] In addition, an embodiment of the present invention also provides an electronic device, which includes: a processor and a memory; the memory is used to store one or more program instructions; the processor is used to run one or more program instructions to execute the steps of a method for automatically generating a stratigraphic zoning map as described in any of the above items.

[0101] It should be noted that, for the detailed description of an electronic device provided in an embodiment of the present invention, reference can be made to the relevant description of a method for automatically generating a stratigraphic zoning map provided in an embodiment of the present application, which will not be repeated here.

[0102] In addition, an embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of the method for automatically generating a stratigraphic zoning map as described in any of the above items are implemented.

[0103] It should be noted that, for a detailed description of a computer-readable storage medium provided in an embodiment of the present invention, reference can be made to the relevant description of a method for automatically generating a stratigraphic zoning map provided in an embodiment of the present application, which will not be repeated here.

[0104] Those skilled in the art will appreciate that all or part of the functions of the various methods in the above-mentioned embodiments can be implemented by hardware or by computer programs. When all or part of the functions in the above-mentioned embodiments are implemented by computer programs, the program can be stored in a computer-readable storage medium, and the storage medium can include: read-only memory, random access memory, disk, optical disk, hard disk, etc., and the program is executed by a computer to implement the above-mentioned functions. For example, the program is stored in the memory of the device, and when the program in the memory is executed by the processor, all or part of the above-mentioned functions can be implemented. In addition, when all or part of the functions in the above-mentioned embodiments are implemented by computer programs, the program can also be stored in a storage medium such as a server, another computer, disk, optical disk, flash disk or mobile hard disk, and can be downloaded or copied and saved in the memory of the local device, or the system of the local device is updated, and when the program in the memory is executed by the processor, all or part of the functions in the above-mentioned embodiments can be implemented.

[0105] The above specific examples are used to illustrate the present invention, which is only used to help understand the present invention and is not intended to limit the present invention. For those skilled in the art, according to the concept of the present invention, some simple deductions, modifications or substitutions can be made.

Claims

1. A method for automatically generating a stratigraphic zoning map, characterized in that: The method comprises: Step S1, obtaining multiple drilling data of the target area and preprocessing and extracting information from the drilling data; Step S2, establishing a unified stratigraphic sequence, dividing the boreholes into stratigraphic layers based on the unified stratigraphic sequence, making the layers of each borehole consistent with the layers of the unified stratigraphic sequence, and unifying the layers of all boreholes; Step S3, taking each drilling position as a sampling point, extracting all sampling points of each stratum after division, taking into account the regional boundary constraints, obtaining all sampling points and boundary point distribution of each stratum, and determining the attribute information of each point, the attribute information of each point is used to mark whether the current stratum is missing; Step S4, triangulate the sampling points and boundary points of each stratum to obtain an irregular triangulated network, traverse each triangle, and determine whether the current triangle needs to be pinched out according to the attribute information of the three points that make up the triangle. If pinching out is required, the pinching out position is calculated to determine the pinching out position point; Step S5, adding the pinch-out location point to the sampling point and boundary point distribution of the current stratum, and assigning the attributes of the pinch-out location point as if there is no missing current stratum, deleting the points in the sampling point and boundary point distribution that are missing from the current stratum according to the attribute information of each point, re-triangulating the remaining points and obtaining the stratigraphic plane partition map of the current stratum.

2. The method for automatically generating a stratigraphic zoning map according to claim 1, characterized in that: The step S1 specifically includes: Collect original drilling data and profile data, and convert profile data into drilling data through data extraction; Check and modify the obtained drilling data, and screen out the drilling holes with missing formations, formation inversions and formation duplications; Information is extracted based on the collected drilling data, and the extracted information is used for statistics and analysis to obtain the parameter information required for pinch-out position calculation.

3. The method for automatically generating a stratigraphic zoning map according to claim 1, characterized in that: The step S2 specifically includes: Within the stratum, according to the preset stratigraphic division order, combined with the original stratification information of the boreholes within the area, the geology and the overall geological structure characteristics of the region, the strata are subdivided into similar levels to establish a unified stratigraphic sequence.

4. The method for automatically generating a stratigraphic zoning map according to claim 3, characterized in that: The step S2 specifically includes: After the unified stratigraphic sequence is determined, for the original borehole strata that are inconsistent with the unified stratigraphic sequence, the corresponding zero-thickness layer is inserted to make the strata of each borehole consistent with the strata of the unified stratigraphic sequence.

5. The method for automatically generating a stratigraphic zoning map according to claim 1, characterized in that: The step S3 specifically includes: For each stratum, each sampling point may contain only two situations: containing the corresponding stratum or not containing the corresponding stratum. The attribute value of the borehole sampling point that does not contain the corresponding stratum, i.e., the stratum is missing, is assigned a value of 1; the attribute value of the borehole sampling point that contains the corresponding stratum, i.e., the stratum is not missing, is assigned a value of 0; The attribute value of the boundary point is determined by interpolating the sampling points of the existing attribute value, or by finding the nearest sampling point in the vicinity and taking the attribute value of the nearest sampling point in the vicinity as the attribute value of the current boundary point.

6. The method for automatically generating a stratigraphic zoning map according to claim 5, characterized in that: The step S4 specifically includes: If the attribute values ​​of the three points that make up the triangle are all 0 or all 1, there is no need to determine the pinch-out position; If the attributes of the three points that make up the triangle contain both 1 and 0, then it is necessary to perform a pinch-out process, and there must be two edges in the triangle where the pinch-out points appear.

7. The method for automatically generating a stratigraphic zoning map according to claim 6, characterized in that: The step S4 specifically includes: For the points that need to be pinched out, the pinch-out coefficient on the boundary of the existing strata is calculated according to the pinch-out coefficient calculation method of the angular unconformity strata, and the pinch-out point position is calculated according to the pinch-out coefficient.

8. The method for automatically generating a stratigraphic zoning map according to claim 5, characterized in that: The step S5 specifically includes: Assign the attribute value of the pinch-out position point to 1, filter out the points with an attribute value of 0, and leave the set of points with an attribute value of 1. Use the triangulation algorithm again to triangulate the points with an attribute value of 1 to obtain an irregular triangulated network, which is the stratigraphic partition of the corresponding stratigraphic layer. The boundary extracted from the irregular triangulated network is the outer boundary of the stratigraphic partition; according to the stratigraphic partition boundary conditions, the stratigraphic partition boundary is smoothed using a smoothing algorithm.

9. The method for automatically generating a stratigraphic zoning map according to claim 1, characterized in that: The method further comprises: The triangulation is a triangulation algorithm that combines the principle of Delaunay triangulation and triangulates adjacent rings of polygons from the outside to the inside on the basis of the relationship between polygon vertices, while taking into account "island polygons".

10. A system for automatically generating a stratigraphic zoning map, characterized in that: The system comprises: A data acquisition module is used to obtain multiple borehole data in the target area and perform preprocessing and information extraction on the borehole data; A stratigraphic division module is used to establish a unified stratigraphic sequence, divide the boreholes into stratigraphic layers based on the unified stratigraphic sequence, make the layers of each borehole consistent with the layers of the unified stratigraphic sequence, and unify the layers of all boreholes; The attribute determination module is used to take each drilling position as a sampling point, extract all sampling points of each stratum after division, and consider the regional boundary constraints at the same time, obtain all sampling points and boundary point distribution of each stratum, and determine the attribute information of each point. The attribute information of each point is used to mark whether there is a missing current stratum; The pinch-out position determination module is used to triangulate the sampling points and boundary points of each stratum to obtain an irregular triangulated network, traverse each triangle, and determine whether the current triangle needs to be pinched out based on the attribute information of the three points that make up the triangle. If pinch-out is required, the pinch-out position calculation is performed to determine the pinch-out position point; The stratigraphic zoning map generation module is used to add the pinch-out location points to the sampling points and boundary point distribution of the current stratigraphic layer, and assign the attributes of the pinch-out location points as if there is no missing current stratigraphic layer, delete the points in the sampling points and boundary point distribution that are missing from the current stratigraphic layer according to the attribute information of each point, re-triangulate the remaining points and obtain the stratigraphic plane zoning map of the current stratigraphic layer.

Images