Divided implicit three-dimensional stratum modeling method and device for drilling data of any scale

Abstract

The invention provides a divide-and-conquer implicit three-dimensional stratum modeling method and device for drilling data of any scale, and belongs to the technical field of three-dimensional stratum modeling.The method comprises the steps that region division is carried out based on the drilling data of any data scale, and a plurality of modeling sub-regions are obtained; extracting control points; according to the control points, independently resolving a stratum implicit interface function of each modeling sub-region, and constructing a symbolic distance function based on the boundary of the modeling sub-region; in each modeling sub-region, carrying out implicit Boolean operation on the stratum implicit interface function and the sign distance function to obtain a constrained implicit geological interface in the modeling sub-region; and according to the three-dimensional geological interface function, carrying out implicit fusion on the constraint implicit geological interfaces in the plurality of modeling sub-regions with the common points and the overlapping regions, and constructing a global implicit three-dimensional geological model of the region to be modeled. According to the method, rapid and accurate three-dimensional geologic model construction of drilling data of any number and scale can be realized.

Description

Technical Field

[0001] The present invention belongs to the technical field of three-dimensional stratum modeling, and in particular relates to a divide-and-conquer implicit three-dimensional stratum modeling method and device for drilling data of any scale. Background Art

[0002] Most of the urban 3D geological modeling based on borehole data is for small-volume boreholes in local engineering areas. These 3D geological models serve and solve engineering scale problems. They are difficult to reflect the overall stratigraphic development and soil distribution characteristics in a large urban area, and cannot assist in decision-making on the deployment and implementation of urban underground space construction projects. With the gradual increase in the scale of borehole data, the amount of data has gradually become huge, and the amount of newly added borehole data can reach tens of thousands or even hundreds of thousands. In the work of urban 3D geological modeling using borehole data, the commonly used 3D geological modeling methods are mainly divided into two categories: explicit modeling and implicit modeling. The former includes classic modeling methods such as parametric surface method and GTP (Generalized Tri-Prism, based on generalized triangular prism). The latter mainly includes radial basis functions and their derived implicit surface functions, including Hermite radial basis function (HRBF, High-Resolution Basis Function), AdaHRBF (gradient-adaptive Hermite RBF), Fast-RBF, locally compactly supported RBF (CSRBF, Compactly Supported RadialBasis Function), etc., which are the implicit surfaces commonly used in the current implicit three-dimensional geological modeling methods.

[0003] Explicit 3D geological modeling is a WYSIWYG 3D geological modeling method. The 3D geological modeling of massive borehole data in cities has a huge amount of interaction and is prone to topological contradictions. The late verification and interactive correction of the model make the modeling efficiency low, and the singleness of the grid scale also makes the reconstruction and hierarchical visualization of large-scale 3D geological models difficult. Implicit 3D geological modeling calculates the implicit equations of the isosurface representing the geometric shape of the geological body and uses a series of implicit function visualization methods to establish a 3D geological model. For the 3D geological modeling of limited borehole data in local areas, the efficiency, automation and spatial analysis capabilities of implicit 3D geological modeling can meet the needs of 3D geological modeling and spatial analysis at the engineering scale. However, for massive borehole data, the ability to perform global solution and modeling based on a single implicit function in the entire study area is ultimately limited. When the borehole data involved in the modeling continues to increase to more than 10,000 and the scope of urban modeling is further expanded, the modeling efficiency is often low, which greatly increases the time cost and cannot meet the real-time requirements of urban 3D geological models and spatial analysis. At the same time, in the context of massive drilling data, it is difficult to ensure the details of the geological body, resulting in a decrease in the accuracy and reliability of the established model, affecting subsequent analysis and application. It is highly dependent on hardware and requires more powerful computing power and storage resources. Ordinary hardware and software equipment are difficult to meet the needs, thus limiting the widespread application of 3D geological models in underground space analysis at the urban or even larger scale level. Summary of the invention

[0004] In response to the shortcomings of the prior art, the present application proposes a divide-and-conquer implicit three-dimensional stratigraphic modeling method and device for borehole data of any scale, which realizes fast, stable and accurate implicit three-dimensional geological modeling for boreholes of any scale (including tens of thousands and hundreds of thousands or above) in a wide range, and reduces the three-dimensional geological modeling method of large-scale borehole data's excessive reliance on hardware requirements.

[0005] In the first aspect, the present application proposes a divide-and-conquer implicit three-dimensional stratigraphic modeling method for borehole data of any scale, comprising:

[0006] Step S1: for drilling data of any data scale, the area to be modeled is divided based on Thiessen polygons, and clustered with Thiessen polygon units to obtain multiple modeling sub-areas;

[0007] Step S2: extracting control points of implicit stratum interfaces at different modeling sub-regions and neighborhood boundaries after clustering Thiessen polygon units;

[0008] Step S3: solving the implicit interface function of the stratum of each modeling sub-area separately according to the control points, and constructing a signed distance function based on the boundary of the modeling sub-area;

[0009] Step S4: in each modeling sub-region, performing implicit Boolean operation on the formation implicit interface function and the signed distance function to obtain the constrained implicit geological interface in the modeling sub-region;

[0010] Step S5: According to the three-dimensional geological interface function, the constrained implicit geological interfaces in multiple modeling sub-areas with common points and overlapping areas are implicitly merged to construct a global implicit three-dimensional geological model of the area to be modeled.

[0011] For drilling data of any data scale, the area to be modeled is divided based on Thiessen polygons, and clustered with Thiessen polygon units to obtain multiple modeling sub-areas, including:

[0012] Step S1.1: In the area to be modeled, grid points are evenly added according to the set grid accuracy, and the grid points and the drilling position points are used as the center points of the Thiessen polygons to construct the Thiessen polygons of the area to be modeled;

[0013] Step S1.2: assigning a flag value to each Thiessen polygon, wherein the flag value is used to distinguish the Thiessen polygon corresponding to the drilling location point from the Thiessen polygon corresponding to the grid point;

[0014] Step S1.3: Select the Thiessen polygons corresponding to the ungrown partitioned drilling location points as seed points of the new engineering drilling data set and put them into the stack;

[0015] Step S1.4: When the stack is not empty, the seed points in the stack are popped out one by one, and the Thiessen polygons corresponding to the currently popped seed points are added to the new engineering drilling data set and marked as grown;

[0016] Step S1.5: traverse the neighborhood Thiessen polygons corresponding to the current seed point, and search the neighborhood of the neighborhood Thiessen polygons until the search depth reaches a preset search depth, and obtain the neighborhood polygon set of the current seed point;

[0017] Step S1.6: Traverse the neighborhood polygon set. If the current polygon is the Thiessen polygon corresponding to the drilling location point, push the current polygon into the stack as a new seed point, and repeat steps S1.2 to S1.4 until the stack is empty, completing the clustering of the new engineering drilling data set.

[0018] The extraction of control points of implicit stratum interfaces at different modeling sub-areas and neighborhood boundaries after clustering Thiessen polygon units includes:

[0019] Step S2.1: The method of extracting control points of implicit stratum interfaces at different modeling sub-areas and neighborhood boundaries is divided into Thiessen polygon areas at drilling locations and Thiessen polygon areas at non-drilling locations;

[0020] Step S2.2: In the Thiessen polygon area of ​​the drilling position, extract the bottom points of the strata from the drilling data as the control points of the implicit surface to form a first control point set C1;

[0021] Step S2.3: In the Thiessen polygon area of ​​the drilling location, the midpoints and elevation values ​​on each side of the Thiessen polygon corresponding to the drilling location point are taken to form a second control point set C2. The first control point set C1 and the second control point set C2 together constitute the control points of the implicit geological interface function in the Thiessen polygon area of ​​the drilling location;

[0022] Step S2.4: In the Thiessen polygon area at the non-drilling location, the midpoint of the Thiessen polygon boundary corresponding to the grid point and the center point of the Thiessen polygon figure corresponding to the grid point are selected as the third control point set C3. The third control point set C3 constitutes the control point of the implicit geological interface function in the Thiessen polygon area at the non-drilling location.

[0023] The elevation value is calculated as follows:

[0024]

[0025] Among them, H is the elevation value, h max is the maximum surface elevation of the current modeling sub-area, h0 is the average elevation of the bottom of the surrounding drilling data, L max The maximum growth search depth d is the cluster center in the Thiessen polygon area of ​​the drilling location max The distance from the center of the Thiessen polygon, L i It is the distance from the cluster center in the Thiessen polygon area of ​​the drilling location to the center of the Thiessen polygon area of ​​the current non-drilling location.

[0026] According to the control points, the implicit interface function of the stratum in each modeling sub-area is solved separately, and the calculation formula is as follows:

[0027]

[0028] Among them, α i is the first parameter of the implicit HRBF function, β i is the second parameter of the implicit HRBF function, is the basis function of the implicit geological interface, H is the Hermite matrix, x j is the center control point of the current spatial implicit 3D surface, x i is the control point of the i-th spatial implicit 3D surface, n j is the normal vector at the surface control point, i is the i-th point, j is the j-th point, n is the number of points in the region, f(x j ) Implicit interface function of the formation.

[0029] The signed distance function is constructed based on the boundary of the modeling sub-region, and the calculation formula is as follows:

[0030]

[0031] Among them, Ω is a subset of the metric space U, d is the metric, Ω represents the boundary. For any u∈U, inf represents the infimum. u is the current point for calculating the signed distance function value, and y is the region boundary point.

[0032] In each modeling sub-region, an implicit Boolean operation is performed on the stratum implicit interface function and the signed distance function to obtain a constrained implicit geological interface in the modeling sub-region, including:

[0033] Step S4.1: Substitute the control point into the formation implicit interface function to calculate the formation implicit interface function value;

[0034] Step S4.2: Calculate the signed distance function value according to the signed distance function;

[0035] Step S4.3: Calculate the intersection of the implicit interface function value and the signed distance function value in space, constrain the implicit geology within the interface modeling sub-region, and take a negative value outside the modeling sub-region to obtain the constrained implicit geology interface within the modeling sub-region.

[0036] According to the three-dimensional geological interface function, the constrained implicit geological interfaces in multiple modeling sub-areas with common points and overlapping areas are implicitly fused to construct a global implicit three-dimensional geological model of the area to be modeled. The calculation formula is as follows:

[0037]

[0038] R p :

[0039] Among them, f1 and f2 are implicit geological interface functions of any two regions, p is the order of smooth fusion, R p is the global implicit 3D geological model of the area to be modeled, G i′ is the implicit constrained geological interface function obtained by Boolean intersection of the implicit geological interface function in the modeling sub-region of the i′th region to be fused and the signed distance function of the boundary, V i′ is the fusion result of the implicit constraint geological interface function between regions of the i′th region to be fused, V i′-1is the fusion result of the implicitly constrained geological interface function between regions of the i′-1th region to be fused, n′ is the number of regions to be fused, V1 is the fusion result of the implicitly constrained geological interface function between regions of the first region to be fused, and G1 is the implicitly constrained geological interface function obtained by the Boolean intersection of the implicit geological interface function in the modeling sub-region of the first region to be fused and the signed distance function of the boundary.

[0040] In the second aspect, the present application proposes a divide-and-conquer implicit three-dimensional stratigraphic modeling device for borehole data of any scale, comprising:

[0041] The regional clustering module is used to divide the modeling area based on Thiessen polygons for drilling data of any data scale, and cluster them with Thiessen polygon units to obtain multiple modeling sub-areas;

[0042] A control point extraction module is used to extract control points from implicit stratum interfaces at different modeling sub-areas and neighborhood boundaries after Thiessen polygon unit clustering;

[0043] A distance function construction module is used to solve the implicit interface function of the stratum of each modeling sub-area according to the control points, and to construct a signed distance function based on the boundary of the modeling sub-area;

[0044] A geological interface calculation module is used to perform implicit Boolean operations on the stratum implicit interface function and the signed distance function in each modeling sub-area to obtain a constrained implicit geological interface in the modeling sub-area;

[0045] The three-dimensional geological model output module is used to implicitly merge the constrained implicit geological interfaces in multiple modeling sub-areas with common points and overlapping areas according to the three-dimensional geological interface function, and construct a global implicit three-dimensional geological model of the area to be modeled.

[0046] In a third aspect, the present application proposes an electronic device comprising: one or more processors, and a memory, wherein the memory is used to store instructions, and when the instructions are executed by the one or more processors, the one or more processors execute the divide-and-conquer implicit three-dimensional stratigraphic modeling method for drilling data of any scale.

[0047] In a fourth aspect, the present application proposes a computer-readable storage medium storing executable instructions, which, when executed, enable a processor to perform the divide-and-conquer implicit three-dimensional stratigraphic modeling method for drilling data of any size.

[0048] Beneficial effects:

[0049] The present application proposes a method and device for dividing and conquering implicit three-dimensional stratigraphic modeling of borehole data of any scale, which realizes the rapid and refined three-dimensional geological modeling of tens of thousands of borehole data. The actual application results show that the method of the present application has a great advantage in modeling efficiency compared with the modeling results obtained by the implicit three-dimensional geological modeling method of global overall solution, while ensuring the high precision of the model and the ability to restore local details; it is suitable for implicit three-dimensional geological modeling tasks in urban areas with various geological conditions, and the stratigraphic coverage and exposure range, soil volume, etc. in the geological model are all in line with geological semantics; implicit three-dimensional geological interface functions can be constrained modeling in different modeling areas; it reduces the dependence on hardware, so that ordinary computing equipment can also be competent for the task of massive borehole three-dimensional geological modeling in large cities. When performing three-dimensional geological modeling on dozens or one or two hundred evenly distributed borehole data in a small area, the method of the present application has no difference in modeling results compared with the commonly used implicit three-dimensional geological modeling method, and the modeling efficiency is still improved. Therefore, the method of the present application is theoretically applicable to implicit three-dimensional geological modeling tasks of borehole data of any scale in a large area. BRIEF DESCRIPTION OF THE DRAWINGS

[0050] Figure 1 A flow chart of a divide-and-conquer implicit three-dimensional stratigraphic modeling method for borehole data of any scale according to an embodiment of the present invention;

[0051] Figure 2 A schematic flow chart of a divide-and-conquer implicit three-dimensional stratigraphic modeling method for borehole data of any scale according to an embodiment of the present invention;

[0052] Figure 3 Principle block diagram of a divide-and-conquer implicit three-dimensional stratigraphic modeling device for drilling data of any scale in an embodiment of the present application. DETAILED DESCRIPTION

[0053] The specific implementation of the present application is further described in detail below in conjunction with the drawings and examples.

[0054] Embodiment 1:

[0055] This embodiment proposes a divide-and-conquer implicit 3D stratigraphic modeling method for borehole data of any scale, such as Figure 1 , Figure 2 As shown, including:

[0056] In the first aspect, the present application proposes a divide-and-conquer implicit three-dimensional stratigraphic modeling method for borehole data of any scale, comprising:

[0057] Step S1: For drilling data of any data scale, the area to be modeled is divided based on Thiessen polygons, and clustered with Thiessen polygon units to obtain multiple modeling sub-areas, including:

[0058] Step S1.1: In the area to be modeled, grid points are evenly added according to the set grid accuracy, and the grid points and the drilling position points are used as the center points of the Thiessen polygons to construct the Thiessen polygons of the area to be modeled;

[0059] In this embodiment, based on the urban drilling data, the study area is divided into regions and Voronoi units are clustered based on the Vorronoi diagram (Thyssen polygon) to form several modeling sub-areas. First, grid points are evenly added in the study area according to the set grid accuracy. The grid points and the drilling location points together constitute the Voronoi diagram center point set P, and the Voronoi diagram of the study area is constructed based on this center point set P;

[0060] Step S1.2: assigning a flag value to each Thiessen polygon, wherein the flag value is used to distinguish the Thiessen polygon corresponding to the drilling location point from the Thiessen polygon corresponding to the grid point;

[0061] In this embodiment, each polygonal unit in the Voronoi diagram of the study area is assigned a flag value to distinguish the Voronoi corresponding to the borehole (Borehole-Voronoi) from the Voronoi corresponding to the ground elevation point (Grid-Voronoi).

[0062] Step S1.3: Select the Thiessen polygons corresponding to the ungrown partitioned drilling location points as seed points of the new engineering drilling data set and put them into the stack;

[0063] In this example, the ungrown Borehole-Voronoi is selected as the new engineering drilling data set Project i The seed points (i=1,2,3,...i∈R) are pushed into the stack;

[0064] Step S1.4: When the stack is not empty, the seed points in the stack are popped out one by one, and the Thiessen polygons corresponding to the currently popped seed points are added to the new engineering drilling data set and marked as grown;

[0065] In this embodiment, if the stack is not empty, the seed points in the stack are traversed, the current seed point is popped out of the stack, and its corresponding Voronoi is added to the Project i (i=1,2,3,...i∈R) and marked as grown.

[0066] Step S1.5: traverse the neighborhood Thiessen polygons corresponding to the current seed point, and search the neighborhood of the neighborhood Thiessen polygons until the search depth reaches a preset search depth, and obtain the neighborhood polygon set of the current seed point;

[0067] In this embodiment, the maximum search depth d is given max, traverse the neighborhood Voronoi polygons of the current seed point Voronoi, and continue searching the neighborhood of the neighborhood polygon until the search depth d reaches the maximum search depth d max , forming the neighborhood polygon set n of the current seed point;

[0068] Step S1.6: Traverse the neighborhood polygon set. If the current polygon is the Thiessen polygon corresponding to the drilling location point, push the current polygon into the stack as a new seed point, and repeat steps S1.2 to S1.4 until the stack is empty, completing the clustering of the new engineering drilling data set.

[0069] In this embodiment, the neighborhood polygon set n is traversed, and if the current Voronoi polygon is a Borehole-Voronoi, it is stacked as a new seed point. Steps S1.2 to S1.4 are repeated until the stack is empty, and the clustering of the new engineering borehole data set is completed.

[0070] During the growth and clustering process, the Borehole-Voronoi and Grid-Voronoi involved in each new borehole dataset are stored. These Voronoi polygons are merged and their boundaries extracted, and the neighborhood topology is reconstructed. Finally, the modeling area is re-divided into several modeling sub-areas based on Voronoi considering the distribution of borehole data.

[0071] Step S2: Extracting control points of implicit stratum interfaces at different modeling sub-regions and neighborhood boundaries after Thiessen polygon unit clustering, including:

[0072] Step S2.1: The method of extracting control points of implicit stratum interfaces at different modeling sub-areas and neighborhood boundaries is divided into Thiessen polygon areas at drilling locations and Thiessen polygon areas at non-drilling locations;

[0073] In this embodiment, the extraction methods of all modeling sub-regions and their internal implicit surface control points are summarized into two categories: Borehole Voronoi Region (Borehole Voronoi Region) and Non-Borehole Voronoi Region (Non-Borehole Voronoi Region).

[0074] Step S2.2: In the Thiessen polygon area of ​​the drilling position, extract the bottom points of the strata from the drilling data as the control points of the implicit surface to form a first control point set C1;

[0075] In this embodiment, in the Borehole Voronoi Region, the bottom points of the strata are extracted from the boreholes as the control points of the implicit surface to form a point set C1.

[0076] Step S2.3: In the Thiessen polygon area of ​​the drilling location, the midpoints and elevation values ​​on each side of the Thiessen polygon corresponding to the drilling location point are taken to form a second control point set C2. The first control point set C1 and the second control point set C2 together constitute the control points of the implicit geological interface function in the Thiessen polygon area of ​​the drilling location;

[0077] The elevation value is calculated as follows:

[0078]

[0079] Among them, H is the elevation value, h max is the maximum surface elevation of the current modeling sub-area, h0 is the average elevation of the bottom of the surrounding drilling data, L max The maximum growth search depth d is the cluster center in the Thiessen polygon area of ​​the drilling location max The distance from the center of the Thiessen polygon, L i The distance from the cluster center in the Thiessen polygon area of ​​the drilling location to the center of the Thiessen polygon area of ​​the current non-drilling location. After clustering the Thiessen polygon units, multiple different modeling areas will be formed. Control points will be extracted from these areas one by one and implicit geological models will be constructed.

[0080] In this embodiment, in order to ensure the continuity of the geological interface between regions and to ensure that the composite implicit geological interface function can transition as smoothly as possible when merging, the midpoints on each side of the Voronoi polygon are taken, and the elevation value is the weighted average of the elevation of the bottom position of the strata in the adjacent regions to form the control point set C2. The point sets C1 and C2 together constitute the control points of the implicit geological interface function of the region.

[0081] Step S2.4: In the Thiessen polygon area at the non-drilling location, the midpoint of the Thiessen polygon boundary corresponding to the grid point and the center point of the Thiessen polygon figure corresponding to the grid point are selected as the third control point set C3. The third control point set C3 constitutes the control point of the implicit geological interface function in the Thiessen polygon area at the non-drilling location.

[0082] In this embodiment, the midpoint of the polygon boundary and the center point of the Voronoi diagram are selected in the Non-Borehole Voronoi Region as the control points of the implicit geological interface to form a point set C3. In order to ensure that the stratigraphic extension within the Borehole-Voronoi and its surrounding neighborhood is reasonably estimated during the modeling process, the elevation value of the geological interface control point is calculated according to formula (1).

[0083] Step S3: According to the control points, the implicit interface function of the stratum of each modeling sub-area is solved separately, and the signed distance function is constructed based on the boundary of the modeling sub-area, including:

[0084] Step S3.1: According to the control points, the implicit interface function of the stratum in each modeling sub-area is solved separately, and the calculation formula is as follows:

[0085]

[0086] Among them, α i is the first parameter of the implicit HRBF function, β i is the second parameter of the implicit HRBF function, is the basis function of the implicit geological interface, H is the Hermite matrix, x j is the center control point of the current spatial implicit 3D surface, x i is the control point of the i-th spatial implicit 3D surface, n j is the normal vector at the surface control point, i is the i-th point, j is the j-th point, n is the number of points in the region, f(x j ) Implicit interface function of the formation. Here, each point needs to calculate the distance function with all other points, so x j The current point is the central control point and all other points are traversed.

[0087] In this embodiment, the implicit geological interface control points of the strata are extracted from the modeling sub-area, and the implicit HRBF three-dimensional geological interface function is constructed according to the calculation formula (2);

[0088] Step S3.2: Construct a signed distance function based on the boundary of the modeled sub-region, and the calculation formula is as follows:

[0089]

[0090] Among them, Ω is a subset of the metric space U, d is the metric, Ω represents the boundary. For any u∈U, inf represents the infimum. u is the current point for calculating the signed distance function value, and y is the region boundary point.

[0091] In this embodiment, the boundary of each modeling sub-region is constructed based on formula (3) to construct a signed distance function SDF, and each modeling sub-region is divided according to the signed distance function.

[0092] Step S4: in each modeling sub-region, performing implicit Boolean operation on the formation implicit interface function and the signed distance function to obtain the constrained implicit geological interface in the modeling sub-region;

[0093] In this embodiment, implicit Boolean operations are used in each modeling sub-region to perform implicit Boolean operations on the implicit function of the stratum interface and the boundary signed distance function, thereby completing the constrained implicit three-dimensional geological modeling in the sub-region.

[0094] The specific steps are as follows:

[0095] Step S4.1: Substitute the control point into the formation implicit interface function to calculate the formation implicit interface function value;

[0096] In this embodiment, the implicit function value of the constrained geological interface of the stratum in the modeling area consists of two parts. First, the point is brought into the HRBF implicit surface function of the stratum in the modeling sub-area to calculate the HRBF function value F.

[0097] Step S4.2: Calculate the signed distance function value according to the signed distance function;

[0098] In this embodiment, the SDF function value G is calculated based on the signed distance function SDF constructed based on the modeling area boundary;

[0099] Step S4.3: Calculate the intersection of the implicit interface function value and the signed distance function value in space, constrain the implicit geology within the interface modeling sub-region, and take a negative value outside the modeling sub-region to obtain the constrained implicit geological interface within the modeling sub-region.

[0100] In this embodiment, based on the calculation formula (4), the part where F and G intersect in space is calculated by implicit Boolean operation, that is, the space where the two implicit surfaces take positive values. Through the constraint of SDF, the original formation implicit HRBF interface is constrained within the modeling sub-region, and the values ​​outside the region are all negative, thereby avoiding the influence on the irrelevant grids outside the region;

[0101] f1∩1 f2=min(f1,f2) (4)

[0102] Step S5: According to the three-dimensional geological interface function, the constrained implicit geological interfaces in multiple modeling sub-areas with common points and overlapping areas are implicitly merged to construct a global implicit three-dimensional geological model of the area to be modeled.

[0103] In this embodiment, based on the general form of R-Function, that is, calculation formula (5), the constrained implicit HRBF geological interfaces in multiple modeling sub-areas with common points and overlapping areas are implicitly fused to construct a globally unique implicit three-dimensional geological interface function of the modeling area.

[0104] The main steps are as follows:

[0105] Step S5.1: Obtain the constrained implicit three-dimensional geological interface function R of several modeling sub-areas through step S4;

[0106] Step S5.2: Based on the calculation formula (5), the implicit three-dimensional geological interface fusion recursive formula of several regions is derived, as shown in formula (6), and the constrained implicit three-dimensional geological interface functions R of several modeling sub-regions are implicitly fused in combination with the neighborhood topological relationship of the modeling sub-regions to obtain a globally unique implicit three-dimensional geological interface function in the entire modeling region;

[0107]

[0108] Among them, f1 and f2 are implicit geological interface functions of any two regions, p is the order of smooth fusion, R p is the global implicit 3D geological model of the area to be modeled, G i′ is the implicit constrained geological interface function obtained by Boolean intersection of the implicit geological interface function in the modeling sub-region of the i′th region to be fused and the signed distance function of the boundary, V i′ is the fusion result of the implicit constraint geological interface function between regions of the i′th region to be fused (it can be understood that: V i′ is the fusion result of this iteration), V i′-1 is the fusion result of the implicit constraint geological interface function between regions of the i′-1th region to be fused (V i′-1 is the fusion result of the previous iteration), n′ is the number of regions to be fused, V1 is the fusion result of the implicit constrained geological interface function between regions of the first region to be fused, and G1 is the implicit constrained geological interface function obtained by Boolean intersection of the implicit geological interface function and the signed distance function of the boundary in the modeling sub-region of the first region to be fused.

[0109] Step S5.3: Based on step S5.1 and step S5.2, the moving tetrahedron method is used to make the implicit function of the grid field of the entire tetrahedron area explicit, so as to complete the rapid implicit three-dimensional geological modeling of large-scale drilling data in a large area.

[0110] The present embodiment proposes a divide-and-conquer implicit three-dimensional stratigraphic modeling method for borehole data of any scale. By partitioning the large-scale borehole data of a large-scale study area based on the Voronoi diagram, the large-scale borehole data is divided into spatial regions that conform to geological and geographical semantics, and then the regional clustering is completed by the growth algorithm, thereby generating several modeling sub-regions. For each modeling sub-region, the implicit Boolean operation is used to constrain the shape and range of the HRBF implicit three-dimensional geological interface function in the region through the signed distance function SDF to construct the constrained implicit three-dimensional geological interface function of each modeling sub-region. The implicit three-dimensional geological interface functions of several modeling sub-regions are fused through the implicit geological interface fusion recursion of the R-Function to generate the unique implicit three-dimensional geological interface function of each stratum in the study area. Finally, the implicit function explicit method is used to complete the rapid implicit three-dimensional geological modeling of large-scale borehole data in a large area. The present embodiment is simple to implement, has significant effects, and meets the requirements of practical engineering applications.

[0111] Embodiment 2:

[0112] This embodiment proposes a divide-and-conquer implicit three-dimensional stratigraphic modeling device for drilling data of any scale, such as Figure 3 As shown, it includes: a regional clustering module, a control point extraction module, a distance function construction module, a geological interface calculation module and a three-dimensional geological model output module;

[0113] The regional clustering module, the control point extraction module, the distance function construction module, the geological interface calculation module and the three-dimensional geological model output module are connected in sequence;

[0114] The regional clustering module is used to divide the modeling area based on Thiessen polygons for drilling data of any data scale, and cluster them with Thiessen polygon units to obtain multiple modeling sub-areas;

[0115] A control point extraction module is used to extract control points from implicit stratum interfaces at different modeling sub-areas and neighborhood boundaries after Thiessen polygon unit clustering;

[0116] A distance function construction module is used to solve the implicit interface function of the stratum of each modeling sub-area according to the control points, and to construct a signed distance function based on the boundary of the modeling sub-area;

[0117] A geological interface calculation module is used to perform implicit Boolean operations on the stratum implicit interface function and the signed distance function in each modeling sub-area to obtain a constrained implicit geological interface in the modeling sub-area;

[0118] The three-dimensional geological model output module is used to implicitly merge the constrained implicit geological interfaces in multiple modeling sub-areas with common points and overlapping areas according to the three-dimensional geological interface function, and construct a global implicit three-dimensional geological model of the area to be modeled.

[0119] Embodiment 3:

[0120] This embodiment proposes an electronic device, including: one or more processors, and a memory, wherein the memory is used to store instructions, and when the instructions are executed by the one or more processors, the one or more processors execute the divide-and-conquer implicit three-dimensional stratigraphic modeling method for drilling data of any scale.

[0121] The electronic device may be a mobile phone, a computer or a tablet computer, etc., including a memory and a processor, wherein a computer program is stored in the memory, and when the computer program is executed by the processor, the divide-and-conquer implicit three-dimensional stratigraphic modeling method for borehole data of any scale as described in the embodiment is implemented. It is understood that the electronic device may also include an input / output (I / O) interface and a communication component.

[0122] The processor is used to execute all or part of the steps in the divide-and-conquer implicit three-dimensional stratigraphic modeling method for borehole data of any scale as described in the above embodiment. The memory is used to store various types of data, which may include, for example, instructions of any application or method in the electronic device, and data related to the application.

[0123] The processor can be an application specific integrated circuit (ASIC), a digital signal processor (DSP), a programmable logic device (PLD), a field programmable gate array (FPGA), a controller, a microcontroller, a microprocessor or other electronic components, and is used to execute the divide-and-conquer implicit three-dimensional formation modeling method for drilling data of any scale described in the above embodiments.

[0124] Embodiment 4:

[0125] This embodiment provides a computer-readable storage medium storing executable instructions. When the instructions are executed and implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium.

[0126] The computer software product is stored in a storage medium and includes a number of instructions for enabling a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the divide-and-conquer implicit three-dimensional stratigraphic modeling method for drilling data of any scale as described in the various embodiments of the present application.

[0127] The aforementioned storage media include: flash memory, hard disk, multimedia card, card-type memory (for example, SD (Secure Digital Memory Card) or DX (Memory Data Register, MDR abbreviation, memory data register) memory, etc.), random access memory (RAM), static random access memory (SRAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), programmable read-only memory (PROM), magnetic memory, disk, CD, server, APP (Application, abbreviation of application software) application store and other media that can store program verification codes, on which computer programs are stored. When the computer program is executed by the processor, the various steps of the divide-and-conquer implicit three-dimensional stratigraphic modeling method for drilling data of any scale described above can be implemented.

[0128] The various embodiments in the present application are described in a progressive manner, and the same or similar parts between the various embodiments can be referenced to each other, and each embodiment focuses on the differences from other embodiments.

[0129] The protection scope of the present application is not limited to the above-mentioned embodiments. Obviously, those skilled in the art can make various changes and modifications to the present disclosure without departing from the scope and spirit of the present disclosure. If these changes and modifications fall within the scope of the claims of the present disclosure and their equivalents, the intention of the present disclosure also includes these changes and modifications.

Claims

1. A divide-and-conquer implicit 3D stratigraphic modeling method for borehole data of any scale, characterized in that: include: Step S1: for drilling data of any data scale, the area to be modeled is divided based on Thiessen polygons, and clustered with Thiessen polygon units to obtain multiple modeling sub-areas; Step S2: extracting control points of implicit stratum interfaces at different modeling sub-regions and neighborhood boundaries after clustering Thiessen polygon units; Step S3: solving the implicit interface function of the stratum of each modeling sub-area separately according to the control points, and constructing a signed distance function based on the boundary of the modeling sub-area; Step S4: in each modeling sub-region, performing implicit Boolean operation on the formation implicit interface function and the signed distance function to obtain the constrained implicit geological interface in the modeling sub-region; Step S5: According to the three-dimensional geological interface function, the constrained implicit geological interfaces in multiple modeling sub-areas with common points and overlapping areas are implicitly merged to construct a global implicit three-dimensional geological model of the area to be modeled.

2. The divide-and-conquer implicit 3D stratigraphic modeling method for borehole data of any scale according to claim 1, characterized in that: For drilling data of any data scale, the area to be modeled is divided based on Thiessen polygons, and clustered with Thiessen polygon units to obtain multiple modeling sub-areas, including: Step S1.1: In the area to be modeled, grid points are evenly added according to the set grid accuracy, and the grid points and the drilling position points are used as the center points of the Thiessen polygons to construct the Thiessen polygons of the area to be modeled; Step S1.2: assigning a flag value to each Thiessen polygon, wherein the flag value is used to distinguish the Thiessen polygon corresponding to the drilling location point from the Thiessen polygon corresponding to the grid point; Step S1.3: Select the Thiessen polygons corresponding to the ungrown partitioned drilling location points as seed points of the new engineering drilling data set and put them into the stack; Step S1.4: When the stack is not empty, the seed points in the stack are popped out one by one, and the Thiessen polygons corresponding to the currently popped seed points are added to the new engineering drilling data set and marked as grown; Step S1.5: traverse the neighborhood Thiessen polygons corresponding to the current seed point, and search the neighborhood of the neighborhood Thiessen polygons until the search depth reaches a preset search depth, and obtain the neighborhood polygon set of the current seed point; Step S1.6: Traverse the neighborhood polygon set. If the current polygon is the Thiessen polygon corresponding to the drilling location point, push the current polygon into the stack as a new seed point, and repeat steps S1.2 to S1.4 until the stack is empty, completing the clustering of the new engineering drilling data set.

3. The divide-and-conquer implicit 3D stratigraphic modeling method for borehole data of any scale according to claim 1, characterized in that: The extraction of control points of implicit stratum interfaces at different modeling sub-areas and neighborhood boundaries after clustering Thiessen polygon units includes: Step S2.1: The method of extracting control points of implicit stratum interfaces at different modeling sub-areas and neighborhood boundaries is divided into Thiessen polygon areas at drilling locations and Thiessen polygon areas at non-drilling locations; Step S2.2: In the Thiessen polygon area of ​​the drilling position, extract the bottom points of the strata from the drilling data as the control points of the implicit surface to form a first control point set C1; Step S2.3: In the Thiessen polygon area of ​​the drilling location, the midpoints and elevation values ​​on each side of the Thiessen polygon corresponding to the drilling location point are taken to form a second control point set C2. The first control point set C1 and the second control point set C2 together constitute the control points of the implicit geological interface function in the Thiessen polygon area of ​​the drilling location; Step S2.4: In the Thiessen polygon area at the non-drilling location, the midpoint of the Thiessen polygon boundary corresponding to the grid point and the center point of the Thiessen polygon figure corresponding to the grid point are selected as the third control point set C3. The third control point set C3 constitutes the control point of the implicit geological interface function in the Thiessen polygon area at the non-drilling location.

4. The divide-and-conquer implicit 3D stratigraphic modeling method for borehole data of any scale according to claim 3, characterized in that: The elevation value is calculated as follows: Among them, H is the elevation value, h max is the maximum surface elevation of the current modeling sub-area, h0 is the average elevation of the bottom of the surrounding drilling data, L max The maximum growth search depth d is the cluster center in the Thiessen polygon area of ​​the drilling location max The distance between the centers of Thiessen polygons, L i It is the distance from the cluster center in the Thiessen polygon area of ​​the drilling location to the center of the Thiessen polygon area of ​​the current non-drilling location.

5. The divide-and-conquer implicit 3D stratigraphic modeling method for borehole data of any scale according to claim 1, characterized in that: According to the control points, the implicit interface function of the stratum in each modeling sub-area is solved separately, and the calculation formula is as follows: Among them, α i is the first parameter of the implicit HRBF function, β i is the second parameter of the implicit HRBF function, is the basis function of the implicit geological interface, H is the Hermite matrix, x j is the center control point of the current spatial implicit 3D surface, x i is the control point of the ith spatial implicit 3D surface, n j is the normal vector at the surface control point, i is the i-th point, j is the j-th point, n is the number of points in the region, f(x j ) Implicit interface function of formation; The signed distance function is constructed based on the boundary of the modeling sub-region, and the calculation formula is as follows: Among them, Ω is a subset of the metric space U, d is the metric, Ω represents the boundary. For any u∈U, inf represents the infimum, u is the current point for calculating the signed distance function value, and y is the region boundary point.

6. The divide-and-conquer implicit 3D stratigraphic modeling method for borehole data of any scale according to claim 1, characterized in that: In each modeling sub-region, an implicit Boolean operation is performed on the stratum implicit interface function and the signed distance function to obtain a constrained implicit geological interface in the modeling sub-region, including: Step S4.1: Substitute the control point into the formation implicit interface function to calculate the formation implicit interface function value; Step S4.2: Calculate the signed distance function value according to the signed distance function; Step S4.3: Calculate the intersection of the implicit interface function value and the signed distance function value in space, constrain the implicit geology within the interface modeling sub-region, and take a negative value outside the modeling sub-region to obtain the constrained implicit geology interface within the modeling sub-region.

7. The divide-and-conquer implicit 3D stratigraphic modeling method for borehole data of any scale according to claim 1, characterized in that: According to the three-dimensional geological interface function, the constrained implicit geological interfaces in multiple modeling sub-areas with common points and overlapping areas are implicitly fused to construct a global implicit three-dimensional geological model of the area to be modeled. The calculation formula is as follows: According to the three-dimensional geological interface function, the constrained implicit geological interfaces in multiple modeling sub-areas with common points and overlapping areas are implicitly fused to construct a global implicit three-dimensional geological model of the area to be modeled. The calculation formula is as follows: Among them, f1 and f2 are implicit geological interface functions of any two regions, p is the order of smooth fusion, R p is the global implicit 3D geological model of the area to be modeled, G i , is the implicit constrained geological interface function obtained by Boolean intersection of the implicit geological interface function in the modeling sub-region of the i′th region to be fused and the signed distance function of the boundary, V i′ is the fusion result of the implicit constraint geological interface function between regions of the i′th region to be fused, V i′-1 is the fusion result of the implicitly constrained geological interface function between regions of the i′-1th region to be fused, n′ is the number of regions to be fused, V1 is the fusion result of the implicitly constrained geological interface function between regions of the first region to be fused, and G1 is the implicitly constrained geological interface function obtained by the Boolean intersection of the implicit geological interface function in the modeling sub-region of the first region to be fused and the signed distance function of the boundary.

8. A divide-and-conquer implicit three-dimensional stratigraphic modeling device for borehole data of any scale, characterized in that: include: The regional clustering module is used to divide the modeling area based on Thiessen polygons for drilling data of any data scale, and cluster them with Thiessen polygon units to obtain multiple modeling sub-areas; A control point extraction module is used to extract control points from implicit stratum interfaces at different modeling sub-areas and neighborhood boundaries after Thiessen polygon unit clustering; A distance function construction module is used to solve the implicit interface function of the stratum of each modeling sub-area according to the control points, and to construct a signed distance function based on the boundary of the modeling sub-area; A geological interface calculation module is used to perform implicit Boolean operations on the stratum implicit interface function and the signed distance function in each modeling sub-area to obtain a constrained implicit geological interface in the modeling sub-area; The three-dimensional geological model output module is used to implicitly merge the constrained implicit geological interfaces in multiple modeling sub-areas with common points and overlapping areas according to the three-dimensional geological interface function, and construct a global implicit three-dimensional geological model of the area to be modeled.

9. An electronic device, characterized in that: include: One or more processors, and a memory, wherein the memory is used to store instructions, and when the instructions are executed by the one or more processors, the one or more processors execute the divide-and-conquer implicit three-dimensional stratigraphic modeling method for borehole data of any scale as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: It stores executable instructions, which, when executed, enable the processor to execute the divide-and-conquer implicit three-dimensional stratigraphic modeling method for borehole data of any scale as described in any one of claims 1 to 7.

Images