Method and system for extracting and editing sand body connected skeleton of three-dimensional geologic model

Abstract

The invention provides a method and a system for extracting and editing a sand body connected skeleton of a three-dimensional geologic model, which are used in the field of three-dimensional geologic modelling, and the method comprises the following steps: obtaining a three-dimensional geologic model, constructing a propagation path of a seed point in a three-dimensional space through a seed neighborhood search algorithm based on a sand body grid attribute value in the three-dimensional geologic model, and obtaining a sand body connected skeleton of the three-dimensional geologic model; the seed points are control points for sand body modeling in the three-dimensional geologic model; obtaining a sparse value of the propagation path based on a back propagation algorithm, and adjusting the density of the propagation path by constraining the sparse value; adjusting the density of the propagation path, taking the propagation path as a communication skeleton, and splitting or combining skeleton branches according to the position relationship among the skeleton branches in the communication skeleton to obtain a topological structure of the communication skeleton; and editing and optimizing the sand body in the three-dimensional geologic model based on the topological structure of the connected skeleton. According to the scheme, the sand body connectivity can be effectively identified, and editing optimization of the sand body connectivity of the three-dimensional geologic model is realized.

Description

Technical Field

[0001] The present invention belongs to the technical field of three-dimensional geological modeling, and in particular relates to a method and system for extracting and editing a connected skeleton of a sand body in a three-dimensional geological model. Background Art

[0002] Reservoir geological modeling refers to the use of computer modeling software to establish high-precision reservoir geological models and conduct detailed analysis of the internal structure of oil and gas reservoirs. Currently, mainstream modeling methods are mainly divided into two types, namely deterministic modeling and stochastic modeling.

[0003] Deterministic modeling refers to providing deterministic predictions for unknown areas between wells based on known deterministic data. Common deterministic modeling reservoir prediction methods include reservoir sedimentology, kriging, and reservoir seismic modeling. Uncertainty modeling utilizes known information, based on random functions, and applies stochastic simulation methods to generate selectable, equally probabilistic, and highly accurate models that reflect the spatial distribution of variables. Common stochastic simulation methods include truncated Gaussian simulation, sequential Gaussian simulation, and sequential indicator simulation. Furthermore, various other modeling methods have emerged, including multi-point geostatistical stochastic simulation, phase-controlled modeling, and well-seismic integration.

[0004] Determining the connectivity of underground sand bodies has always been a core research topic in the field of reservoir modeling. The various methods mentioned above all establish three-dimensional geological models based on input data. In the case of disconnection between two wells, the models constructed by these modeling methods are usually unable to directly edit the connectivity between the wells. Summary of the Invention

[0005] In view of this, an embodiment of the present invention provides a method and system for extracting and editing the connectivity skeleton of sand bodies in a three-dimensional geological model, which is used to solve the problem that the connectivity between wells in the currently constructed three-dimensional geological model cannot be edited.

[0006] In a first aspect of an embodiment of the present invention, a method for extracting and editing a connected skeleton of a sand body in a three-dimensional geological model is provided, comprising: Obtaining a three-dimensional geological model, and constructing a propagation path of seed points in three-dimensional space based on the grid attribute values of the sand bodies in the three-dimensional geological model using a seed neighborhood search algorithm. The seed points are control points for sand body modeling in the three-dimensional geological model. Obtaining a sparse value of the propagation path based on a back-propagation algorithm, and adjusting the density of the propagation path by constraining the sparse value; After adjusting the density of the propagation path, the propagation path is used as a connected skeleton. According to the positional relationship between the skeleton branches in the connected skeleton, the skeleton branches are split or merged to obtain the topological structure of the connected skeleton. Based on the topological structure of the connected skeleton, the sand bodies in the 3D geological model are edited and optimized.

[0007] In a second aspect of an embodiment of the present invention, a system for extracting and editing a connected skeleton of a sand body in a three-dimensional geological model is provided, comprising: A path construction module is used to obtain a three-dimensional geological model and construct a propagation path of seed points in three-dimensional space based on the grid attribute values of the sand body in the three-dimensional geological model through a seed neighborhood search algorithm. The seed points are control points for sand body modeling in the three-dimensional geological model. A sparse constraint module, configured to obtain a sparse value of the propagation path based on a back propagation algorithm, and adjust the density of the propagation path by constraining the sparse value; The splitting and merging module is used to adjust the density of the propagation path, use the propagation path as a connected skeleton, and split or merge the skeleton branches according to the positional relationship between the skeleton branches in the connected skeleton to obtain the topological structure of the connected skeleton; The editing and optimization module is used to edit and optimize the sand bodies in the three-dimensional geological model based on the topological structure of the connected skeleton.

[0008] In a third aspect of an embodiment of the present invention, an electronic device is provided, comprising a memory, a processor, and a computer program stored in the memory and executable by the processor, wherein the processor implements the steps of the method described in the first aspect of the embodiment of the present invention when executing the computer program.

[0009] In a fourth aspect of an embodiment of the present invention, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the method provided in the first aspect of the embodiment of the present invention are implemented.

[0010] In an embodiment of the present invention, a seed neighborhood search algorithm is used to determine the propagation path of a seed point, a sparse value is used to control the sparseness of the propagation path, and each branch path is merged and split to construct a skeleton topology structure. Based on the skeleton topology structure, the connectivity of the sand body model in the three-dimensional geological model is edited and optimized, thereby effectively identifying the connectivity of the sand body model and realizing direct editing of the connectivity of the sand bodies between wells in the three-dimensional geological model, optimizing the spatial distribution of the sand body model, and improving the accuracy and adaptability of the three-dimensional geological model. BRIEF DESCRIPTION OF THE DRAWINGS

[0011] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the embodiments or descriptions of the prior art. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0012] Figure 1A schematic flow chart of a method for extracting and editing a connected skeleton of a sand body in a three-dimensional geological model provided by one embodiment of the present invention; Figure 2 A schematic diagram of seed propagation provided by one embodiment of the present invention; Figure 3 A schematic diagram of back propagation provided for one embodiment of the present invention; Figure 4 A schematic diagram of density constraints provided by one embodiment of the present invention; Figure 5 A schematic diagram of skeleton disassembly provided in one embodiment of the present invention; Figure 6 A schematic diagram of skeleton merging provided in one embodiment of the present invention; Figure 7 A schematic diagram of similarity search provided by one embodiment of the present invention; Figure 8 A schematic structural diagram of a system for extracting and editing a connected skeleton of a three-dimensional geological model sand body provided by one embodiment of the present invention; Figure 9 The present invention provides a schematic structural diagram of an electronic device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0013] In order to make the purpose, features, and advantages of the present invention more obvious and easy to understand, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the embodiments described below are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0014] It should be understood that the terms "including" and similar expressions in the specification, claims, and drawings of the present invention are intended to cover non-exclusive inclusions. For example, a process, method, system, or apparatus comprising a series of steps or units is not limited to the listed steps or units. Furthermore, the terms "first" and "second" are used to distinguish between different objects and are not intended to describe a specific order.

[0015] See also Figure 1 , a schematic flow chart of a method for extracting and editing a connected skeleton of a sand body in a three-dimensional geological model provided by an embodiment of the present invention, comprising: S101, obtaining a three-dimensional geological model, and constructing a propagation path of seed points in three-dimensional space using a seed neighborhood search algorithm based on sand body grid attribute values in the three-dimensional geological model, wherein the seed points are control points for sand body modeling in the three-dimensional geological model; A 3D geological model is a quantitative 3D model generated by integrating geological, well logging, and geophysical data with various interpretations. It is generally represented by a 3D grid. These grids are based on surfaces, faults, and horizons, and each grid contains attribute values such as porosity, permeability, and water saturation. A sand body, also known as a sandstone body, refers to oil-bearing sandstone. Some highly permeable sandstone layers are separated by less permeable rock layers into several independent, disconnected sand layers, which are referred to as sand bodies. The attributes of the sand body grid can refer to permeability, porosity, or oil saturation.

[0016] The seed neighborhood search algorithm is a grid-based search algorithm for seed points. It searches the attribute values of the grids surrounding the seed points. For a sand body attribute grid, the algorithm propagates the attribute values of the neighboring grids from largest to smallest, obtaining all propagation paths for the seed points. Seed points are key points that identify sand body characteristics and can also be called control points for sand body modeling. They are typically determined from seismic data or other geological data. Based on these seed points, the spatial distribution of sand bodies and oil and gas reservoirs can be determined.

[0017] Among them, the grid with the largest grid attribute value in the sand body is used as the seed point, and the grid of the neighborhood is iteratively obtained through the seed neighborhood search algorithm. According to the principle of grid attribute value from large to small, all the propagation paths of the seed point in three-dimensional space are constructed.

[0018] For example, Figure 2 As shown in the figure, each grid represents its attribute value through a number. In the initial stage, each digital point is the attribute value, 15 is defined as a random seed point (seed), and the propagation rule is from large to small attribute value. The steps are: A searches for the neighbors of seed point (15) and sorts them by priority as (14), (12), (10), with a priority of (14); B searches for the neighbor value of (14), with the neighbors being (15), (11), (8), deletes the parent class (15), and the priority value is (11); C searches for the neighbor value of (11), with the neighbors being (14), (12), (7), deletes the parent class (14), deletes the duplicate (12), and the priority value is (7) and stops when it reaches the boundary.

[0019] Find the neighbors of B's second priority value (8), (14), (7), delete the parent class (14), delete the duplicate (7), and stop when the priority value (8) reaches the boundary. The above steps are completed. Figure 2 During the propagation, the red path is searched. After the traversal is completed, the neighbors with the second and third priority values in A can be searched to obtain the blue path and green path respectively. Therefore, the red path, blue path and green path can be used as the propagation path of the seed point.

[0020] S102, obtaining a sparse value of the propagation path based on a back propagation algorithm, and adjusting the density of the propagation path by constraining the sparse value; The backpropagation algorithm performs backpropagation on the endpoints of a propagation path. Backpropagation can determine the number of branches along the propagation path and the sparsity of the propagation path. The more branches there are, the higher the sparsity of the propagation path. By constraining the sparsity value, the number of branches can be reduced, controlling the density of the propagation path.

[0021] The sparsity value of the propagation path is preset, and the branch paths in the propagation path that are smaller than the preset sparsity value are truncated.

[0022] By setting the sparseness value, the paths in the propagation path that are smaller than the preset sparseness value can be truncated to reduce the density of the propagation path.

[0023] For example, Figure 3 As shown in the figure, for the propagation path, through back propagation, we can know that grid points (7), (8), (9), and (10) have no subclasses, and the path density is 0; grid point (11) has subclass (7), and grid point (12) has subclass (9), and the density is 1; grid point 14 has subclasses (8), (7), and (11), and the density is 3; seed point (15) is the initial parent class, seed is the seed point, and the density value is the sum of the path density values of each subclass, so it is 7.

[0024] Based on the calculation of the path sparsity value, the propagation path can be constrained. The path sparsity density constraint is as follows: Figure 4 As shown, when the sparsity value is set to zero, the propagation path is Figure 4 a, when the sparse value is set to 1, the propagation path is Figure 4 b.

[0025] S103, after adjusting the density of the propagation path, using the propagation path as a connected skeleton, and splitting or merging the skeleton branches according to the positional relationship between the skeleton branches in the connected skeleton to obtain a topological structure of the connected skeleton; The connected skeleton is also the propagation path. Since the propagation path is composed of multiple branch paths, the skeleton branches also correspond to the branch paths in the propagation path. When the branch paths or branch skeletons intersect, the intersecting skeleton branches need to be split or merged.

[0026] When skeleton branches intersect, they can be merged or split according to the positional relationship of the intersection and the level of the skeleton branches.

[0027] Specifically, the levels of the skeleton branches are divided according to the mesh attribute values of the branches in the connected skeleton; When skeleton branches of different levels intersect, if the lower-level skeleton branch is not connected to the endpoint of the higher-level skeleton branch, the lower-level connected skeleton is split; If a lower-level skeleton branch is connected to an endpoint of a higher-level skeleton branch, the lower-level skeleton branch is merged into the higher-level skeleton branch.

[0028] However, for the branch paths with higher mesh attribute values in the seed point propagation path, they can be set as high-level branches, and for the branch paths with lower mesh attribute values, they can be set as low-level branches. When the intersecting skeleton branches are low-level branches and the intersection is an endpoint, the low-level skeleton branches are merged. When the intersection is not an endpoint, the low-level skeleton branches need to be split.

[0029] For example, Figure 5 , when the low-level skeleton branch b2 is not connected to the endpoint of the high-level skeleton branch b1, it is necessary to split the low-level skeleton branch b2 and split it into b3 and b4; like Figure 6 As shown, when the endpoints of the low-level skeleton branch b2 and the high-level skeleton branch b1 are connected or overlapped, the low-level skeleton branch b2 needs to be merged into the high-level skeleton branch b1, that is, when the low-level b2 and the high-level b1 are connected end to end (overlap), they are directly merged; when b1 and b2 are not connected end to end (connected), the longer part of the low-level b2 is indirectly merged into b1, and the shorter part can be retained as the low-level b2.

[0030] S104. Edit and optimize the sand body in the three-dimensional geological model based on the topological structure of the connected skeleton.

[0031] Based on the topological structure of the connected skeleton, the properties of the grids within a certain range around the connected skeleton are edited to achieve editing optimization of the sand body model.

[0032] Specifically, obtain the mesh attribute values in the connected skeleton and the mesh attribute values around the connected skeleton; By calculating the Euclidean distance between the attribute values of the meshes around the connected skeleton and the attribute values of the meshes in the connected skeleton, the extended boundary of the connected skeleton is determined, and the attributes of the meshes within the extended boundary are edited.

[0033] The similarity of the mesh attribute values within a certain range around the connected skeleton branch is calculated, and the connected skeleton is expanded outward based on the similarity calculation results. When the similarity of the skeleton meshes in the first range is lower than that in the second range, the mesh corresponding to the first range is used as the boundary, the outward expansion is terminated, and the meshes in the first range are reassigned. The meshes in the second range can also be reassigned, thus completing the identification of the expansion boundary of the selected skeleton branch and enabling the editing of the corresponding original geological model.

[0034] The similarity metric is the Euclidean distance, which calculates the magnitude of the outer-expanded grid attributes and the current skeleton attribute vector. The smaller the value, the more similar. As Figure 7 shown, assume that A is the grid point where the skeleton is located, with attribute values between 700 and 80, and B is the grid point with attribute values between 700 and 60 outside the skeleton. For similarity judgment, d(A,B) < d(A,C). It can be considered that the attribute values of the grid points in the range of C have mutated, reaching the termination condition for outward search, and it can be considered that the attribute values in the range of B need to be edited. In practice, the sand body attribute values in the range of C can be modified to unconnected mudstone to achieve the purpose of editing and correcting the model attribute values based on the skeleton topological relationship as expected.

[0035] In this embodiment, three-dimensional connectivity skeleton extraction and model editing based on sand body attributes can solve the limitations of traditional reservoir modeling in well - to - well connectivity editing. By using the seed neighborhood search algorithm to determine the propagation path of the seed points, controlling the density of the propagation path based on sparse values, and directly editing and optimizing the connectivity of the sand body model according to the skeleton topological structure. Thus, it can effectively identify the connectivity of the sand body model and optimize the spatial distribution of the unconnected sand body model, improve the accuracy and adaptability of three-dimensional reservoir modeling, and provide reference and support for reservoir connectivity analysis under complex geological conditions.

[0036] It should be understood that the sequence numbers of the steps in the above embodiments do not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation to the implementation process of the embodiments of the present invention.

[0037] Figure 8 FIG. 12 is a schematic structural diagram of a three - dimensional geological model sand body connectivity skeleton extraction and editing system provided by an embodiment of the present invention. The system includes: A path construction module 810, configured to obtain a three - dimensional geological model, and construct a propagation path of seed points in three - dimensional space based on the sand body grid attribute values in the three - dimensional geological model. The seed points are control points for sand body modeling in the three - dimensional geological model; Among them, constructing the propagation path of seed points in three - dimensional space by the seed neighborhood search algorithm includes: Taking the grid with the largest grid attribute value in the sand body as the seed point, and iteratively obtaining the grids in the neighborhood through the seed neighborhood search algorithm, and constructing all the propagation paths of the seed points in three - dimensional space according to the principle of arranging the grid attribute values from large to small.

[0038] A sparse constraint module 820, configured to obtain the sparse value of the propagation path based on the backpropagation algorithm, and adjust the density of the propagation path by constraining the sparse value; Optionally, the sparse constraint module 820 includes: The branch truncation unit is used to preset the sparsity value of the propagation path and truncate the branch path in the propagation path that is less than the preset sparsity value.

[0039] The splitting and merging module 830 is configured to adjust the density of the propagation path, use the propagation path as a connected skeleton, and split or merge the skeleton branches according to the positional relationship between the skeleton branches in the connected skeleton to obtain the topological structure of the connected skeleton; The splitting and merging module 830 includes: A level setting unit is used to classify the levels of the skeleton branches according to the mesh attribute values of the branches in the connected skeleton; A splitting unit is used to split the lower-level connected skeleton when different-level skeleton branches intersect, if the lower-level skeleton branch is not connected to the endpoint of the higher-level skeleton branch; The merging unit is used to merge the lower-level skeleton branch into the higher-level skeleton branch when skeleton branches of different levels intersect and if the lower-level skeleton branch is connected to the endpoint of the higher-level skeleton branch.

[0040] The editing and optimization module 840 is used to edit and optimize the sand bodies in the three-dimensional geological model based on the topological structure of the connected skeleton.

[0041] Among them, the mesh attribute values in the connected skeleton and the mesh attribute values around the connected skeleton are obtained; By calculating the Euclidean distance between the attribute values of the meshes around the connected skeleton and the attribute values of the meshes in the connected skeleton, the extended boundary of the connected skeleton is determined, and the attributes of the meshes within the extended boundary are edited.

[0042] Those skilled in the art will clearly understand that, for the convenience and brevity of description, the specific working processes of the above-described systems and modules can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0043] Figure 9 This is a schematic diagram of the structure of an electronic device provided by one embodiment of the present invention. The electronic device is used for connecting skeleton extraction and model editing optimization of three-dimensional attribute models. Figure 9 As shown, the electronic device 9 of this embodiment includes: a memory 910, a processor 920 and a system bus 930, wherein the memory 910 includes an executable program 9101 stored thereon. It can be understood by those skilled in the art that Figure 9 The electronic device structure shown in the figure does not constitute a limitation to the electronic device, and may include more or fewer components than shown in the figure, or combine certain components, or arrange the components differently.

[0044] The following combination Figure 9 A detailed introduction to the various components of electronic equipment: Memory 910 can be used to store software programs and modules. Processor 920 executes the software programs and modules stored in memory 910 to perform various functional applications and data processing of the electronic device. Memory 910 may primarily include a program storage area and a data storage area. The program storage area may store an operating system and at least one application required for a function (such as sound playback or image playback). The data storage area may store data generated based on the use of the electronic device (such as cached data). Memory 910 may also include high-speed random access memory and non-volatile memory, such as at least one disk storage device, flash memory device, or other volatile solid-state memory device.

[0045] The memory 910 includes an executable program 9101 for the interface generation method. The executable program 9101 can be divided into one or more modules / units, which are stored in the memory 910 and executed by the processor 920 to implement sand body model connectivity skeleton extraction and attribute editing, etc. The one or more modules / units can be a series of computer program instruction segments that can perform specific functions, and the instruction segments are used to describe the execution process of the executable program 9101 in the electronic device 9. For example, the executable program 9101 can be divided into functional modules such as a path construction module, a sparse constraint module, a splitting and merging module, and an editing optimization module.

[0046] The processor 920 is the control center of the electronic device. It connects the various parts of the entire electronic device using various interfaces and lines. By running or executing software programs and / or modules stored in the memory 910 and accessing data stored in the memory 910, it performs various functions of the electronic device and processes data, thereby monitoring the overall status of the electronic device. Optionally, the processor 920 may include one or more processing units; preferably, the processor 920 may integrate an application processor and a modem processor, wherein the application processor primarily processes the operating system, application programs, etc., and the modem processor primarily handles wireless communications. It is understood that the modem processor described above may not be integrated into the processor 920.

[0047] The system bus 930 connects the various functional components within the computer and can transmit data, address information, and control information. It can be a PCI bus, an ISA bus, a CAN bus, or other bus types. Instructions from the processor 920 are transmitted to the memory 910 via the bus, and the memory 910 feeds data back to the processor 920. The system bus 930 facilitates the exchange of data and instructions between the processor 920 and the memory 910. Of course, the system bus 930 can also connect to other devices, such as network interfaces and display devices.

[0048] In an embodiment of the present invention, the executable program executed by the processing 920 included in the electronic device includes: Obtaining a three-dimensional geological model, and constructing a propagation path of seed points in three-dimensional space based on the grid attribute values of the sand bodies in the three-dimensional geological model using a seed neighborhood search algorithm. The seed points are control points for sand body modeling in the three-dimensional geological model. Obtaining a sparse value of the propagation path based on a back-propagation algorithm, and adjusting the density of the propagation path by constraining the sparse value; After adjusting the density of the propagation path, the propagation path is used as a connected skeleton. According to the positional relationship between the skeleton branches in the connected skeleton, the skeleton branches are split or merged to obtain the topological structure of the connected skeleton. Based on the topological structure of the connected skeleton, the sand bodies in the 3D geological model are edited and optimized.

[0049] Those skilled in the art will clearly understand that, for the convenience and brevity of description, the specific working processes of the above-described systems, devices, and modules can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0050] In the above embodiments, the description of each embodiment has its own focus. For parts that are not described or recorded in detail in a certain embodiment, reference can be made to the relevant description of other embodiments.

[0051] As described above, the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit the same. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art should understand that the technical solutions described in the above embodiments can still be modified, or some of the technical features thereof can be replaced by equivalents. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for extracting and editing the connected skeleton of sand bodies in a three-dimensional geological model, characterized in that: include: Obtaining a three-dimensional geological model, and constructing a propagation path of seed points in three-dimensional space based on the grid attribute values of the sand bodies in the three-dimensional geological model using a seed neighborhood search algorithm. The seed points are control points for sand body modeling in the three-dimensional geological model. Obtaining a sparse value of the propagation path based on a back-propagation algorithm, and adjusting the density of the propagation path by constraining the sparse value; After adjusting the density of the propagation path, the propagation path is used as a connected skeleton. According to the positional relationship between the skeleton branches in the connected skeleton, the skeleton branches are split or merged to obtain the topological structure of the connected skeleton. Based on the topological structure of the connected skeleton, the sand bodies in the 3D geological model are edited and optimized.

2. The method according to claim 1, characterized in that The method of constructing a propagation path of a seed point in three-dimensional space by using a seed neighborhood search algorithm based on the sand body grid attribute values in the three-dimensional geological model includes: The grid with the largest grid attribute value in the sand body is taken as the seed point. The neighborhood grid is iteratively obtained through the seed neighborhood search algorithm. All propagation paths of the seed point in three-dimensional space are constructed according to the principle of grid attribute value from large to small.

3. The method according to claim 1, characterized in that The adjusting the density of the propagation path by constraining the sparsity value includes: The sparseness value of the propagation path is preset, and the branch paths in the propagation path that are less than the preset sparseness value are truncated.

4. The method according to claim 1, wherein The process of splitting or merging the skeleton branches according to the positional relationship between the skeleton branches in the connected skeleton to obtain the topological structure of the connected skeleton includes: According to the mesh attribute values of the branches in the connected skeleton, the levels of the skeleton branches are divided; When skeleton branches of different levels intersect, if the lower-level skeleton branch is not connected to the endpoint of the higher-level skeleton branch, the lower-level connected skeleton is split; If a lower-level skeleton branch is connected to an endpoint of a higher-level skeleton branch, the lower-level skeleton branch is merged into the higher-level skeleton branch.

5. The method according to claim 1, wherein The editing and optimization of the sand body in the three-dimensional geological model based on the topological structure of the connected skeleton includes: Get the mesh attribute values in the connected skeleton and the mesh attribute values around the connected skeleton; By calculating the Euclidean distance between the attribute values of the meshes around the connected skeleton and the attribute values of the meshes in the connected skeleton, the extended boundary of the connected skeleton is determined, and the attributes of the meshes within the extended boundary are edited.

6. A three-dimensional geological model sand body connected skeleton extraction and editing system, characterized by: include: A path construction module is used to obtain a three-dimensional geological model and construct a propagation path of seed points in three-dimensional space based on the grid attribute values of the sand body in the three-dimensional geological model through a seed neighborhood search algorithm. The seed points are control points for sand body modeling in the three-dimensional geological model. A sparse constraint module, configured to obtain a sparse value of the propagation path based on a back propagation algorithm, and adjust the density of the propagation path by constraining the sparse value; The splitting and merging module is used to adjust the density of the propagation path, use the propagation path as a connected skeleton, and split or merge the skeleton branches according to the positional relationship between the skeleton branches in the connected skeleton to obtain the topological structure of the connected skeleton; The editing and optimization module is used to edit and optimize the sand bodies in the three-dimensional geological model based on the topological structure of the connected skeleton.

7. The system according to claim 6, characterized in that The sparse constraint module includes: The branch truncation unit is used to preset the sparsity value of the propagation path and truncate the branch path in the propagation path that is less than the preset sparsity value.

8. The system according to claim 6, wherein: The splitting and merging module includes: A level setting unit is used to classify the levels of the skeleton branches according to the mesh attribute values of the branches in the connected skeleton; A splitting unit is used to split the lower-level connected skeleton when different-level skeleton branches intersect, if the lower-level skeleton branch is not connected to the endpoint of the higher-level skeleton branch; The merging unit is used to merge the lower-level skeleton branch into the higher-level skeleton branch when skeleton branches of different levels intersect and if the lower-level skeleton branch is connected to the endpoint of the higher-level skeleton branch.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the steps of the method for extracting and editing the connected skeleton of a three-dimensional geological model sand body are implemented as described in any one of claims 1 to 5.

10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed, the steps of the method for extracting and editing the connected skeleton of sand bodies in a three-dimensional geological model are implemented as described in any one of claims 1 to 5.

Images