A method and system for extracting and editing sand body connected framework of a three-dimensional geological model
By adjusting the sparsity values through seed neighborhood search and backpropagation algorithms, the sand body connectivity of the 3D geological model is optimized, overcoming the limitations of inter-well connectivity editing and improving the model's accuracy and adaptability.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- YANGTZE UNIVERSITY
- Filing Date
- 2025-03-24
- Publication Date
- 2026-05-05
AI Technical Summary
Existing 3D geological models cannot directly edit inter-well connectivity, which limits their ability to determine the connectivity of underground sand bodies.
The propagation path of seed points is constructed by seed neighborhood search algorithm, and the sparsity value is adjusted by backpropagation algorithm to control the density of propagation path. The sand body connectivity of the three-dimensional geological model is optimized by splitting or merging according to the topological structure of the connected skeleton.
This method effectively identifies and optimizes the connectivity of sand bodies between wells, improves the accuracy and adaptability of three-dimensional geological models, and overcomes the limitations of traditional methods in editing inter-well connectivity.
Smart Images
Figure CN120495547B_ABST
Abstract
Description
Technical Field
[0001] This invention belongs to the field of three-dimensional geological modeling technology, and in particular relates to a method and system for extracting and editing the connected skeleton of sand bodies in a three-dimensional geological model. Background Technology
[0002] Reservoir geological modeling refers to the use of computer modeling software to establish high-precision reservoir geological models and to perform detailed dissection of the internal structure of oil and gas reservoirs. Currently, the mainstream modeling methods are mainly divided into two types: 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 reservoir prediction methods include reservoir sedimentology methods, kriging methods, and reservoir seismic modeling. Uncertainty modeling, on the other hand, utilizes known information as a foundation, uses stochastic functions as a theoretical basis, and applies stochastic simulation methods to obtain selectable, equally probable, and high-precision models that reflect the spatial distribution of variables. Common stochastic simulation methods include truncated Gaussian simulation, sequential Gaussian simulation, and sequential indicator simulation. Additionally, various other modeling methods have emerged, such as 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 reservoir modeling. The various methods mentioned above all build three-dimensional geological models based on input data. However, when there is a disconnect between two wells, the models constructed by these modeling methods usually cannot directly edit the connectivity between the wells. Summary of the Invention
[0005] In view of this, embodiments of the present invention provide a method and system for extracting and editing the sand body connectivity skeleton of 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 the present invention, a method for extracting and editing the connected skeleton of sand bodies in a three-dimensional geological model is provided, comprising:
[0007] A three-dimensional geological model is obtained. Based on the sand body mesh attribute values in the three-dimensional geological model, a seed point propagation path in three-dimensional space is constructed through a seed neighborhood search algorithm. The seed point is the control point for sand body modeling in the three-dimensional geological model.
[0008] The sparsity value of the propagation path is obtained based on the backpropagation algorithm, and the density of the propagation path is adjusted by constraining the sparsity value.
[0009] After adjusting the density of the propagation path, the propagation path is used as a connected skeleton. Based on the positional relationship between the skeleton branches, the skeleton branches are split or merged to obtain the topology of the connected skeleton.
[0010] Based on the topological structure of the connected skeleton, the sand bodies in the three-dimensional geological model are edited and optimized.
[0011] In a second aspect of the present invention, a system for extracting and editing the connected skeleton of sand bodies in a three-dimensional geological model is provided, comprising:
[0012] The path construction module is used to obtain a three-dimensional geological model. Based on the sand body mesh attribute values in the three-dimensional geological model, the propagation path of seed points in three-dimensional space is constructed through a seed neighborhood search algorithm. The seed points are the control points for sand body modeling in the three-dimensional geological model.
[0013] The sparse constraint module is used 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.
[0014] The splitting and merging module is used to adjust the density of the propagation path, and then use the propagation path as a connected skeleton. Based on the positional relationship between the skeleton branches, the skeleton branches are split or merged to obtain the topology of the connected skeleton.
[0015] The editing and optimization module is used to edit and optimize sand bodies in a 3D geological model based on the topology of the connected skeleton.
[0016] In a third aspect of the present invention, an electronic device is provided, including a memory, a processor, and a computer program stored in the memory and executable by the processor, wherein the processor executes the computer program to implement the steps of the method as described in the first aspect of the present invention.
[0017] In a fourth aspect of the present invention, a computer-readable storage medium is provided, the computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of the method provided in the first aspect of the present invention.
[0018] In this embodiment of the invention, the propagation path of the seed point is determined by the seed neighborhood search algorithm, the sparsity of the propagation path is controlled by the sparsity value, and the branch paths are merged and split to construct the skeleton topology. Based on the skeleton topology, 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 the direct editing of the inter-well sand body connectivity 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. Attached Figure Description
[0019] To more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0020] Figure 1 A flowchart illustrating a method for extracting and editing the connected skeleton of a three-dimensional geological model sand body, provided in one embodiment of the present invention;
[0021] Figure 2 This is a schematic diagram of seed propagation provided in one embodiment of the present invention;
[0022] Figure 3 A schematic diagram of backpropagation provided for one embodiment of the present invention;
[0023] Figure 4 A schematic diagram of sparse-density constraints provided for one embodiment of the present invention;
[0024] Figure 5 A schematic diagram of skeleton disassembly provided according to an embodiment of the present invention;
[0025] Figure 6 This is a schematic diagram of skeleton merging provided in one embodiment of the present invention;
[0026] Figure 7 This is a schematic diagram of similarity search provided in one embodiment of the present invention;
[0027] Figure 8 This is a schematic diagram of the structure of a three-dimensional geological model sand body connectivity skeleton extraction and editing system provided in one embodiment of the present invention;
[0028] Figure 9 This is a schematic diagram of the structure of an electronic device provided in one embodiment of the present invention. Detailed Implementation
[0029] To make the objectives, features, and advantages of this invention more apparent and understandable, the technical solutions of the embodiments of this invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the embodiments described below are only some embodiments of this invention, and not all embodiments. Based on the embodiments of this invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this invention.
[0030] It should be understood that the terms "comprising" and other similar expressions in the specification, claims, and accompanying drawings of this invention are intended to cover a non-exclusive inclusion, such as a process, method, system, or apparatus that includes a series of steps or units and is not limited to the listed steps or units. Furthermore, "first" and "second" are used to distinguish different objects and are not intended to describe a specific order.
[0031] Please see Figure 1 The present invention provides a flowchart illustrating a method for extracting and editing the connected skeleton of a three-dimensional geological model sand body, comprising:
[0032] S101. Obtain a three-dimensional geological model. Based on the sand body mesh attribute values in the three-dimensional geological model, construct the propagation path of seed points in three-dimensional space through a seed neighborhood search algorithm. The seed points are the control points for sand body modeling in the three-dimensional geological model.
[0033] A three-dimensional geological model is a quantitative three-dimensional model generated by integrating geological, well logging, and geophysical data with various interpretations. It is generally represented by a three-dimensional grid, which is built upon surfaces, faults, and stratigraphic levels. Each grid contains attribute values such as porosity, permeability, and water saturation. A sand body, also known as a sandstone body, refers to an oil-bearing sand body. Some high-permeability sandstone layers are divided into several independent, unconnected sand layers by low-permeability rock layers; these can be called sand bodies. The attribute values of a sand body grid can refer to permeability, porosity, or oil saturation.
[0034] The seed neighborhood search algorithm is a type of search algorithm that searches the neighborhood grid of a seed point. It searches the attribute values of the neighboring grids of a seed point. For a sand body attribute grid, it propagates the attribute values from largest to smallest, obtaining all propagation paths for the seed point. Seed points are key points that identify the characteristics of a sand body; they can also be called control points in sand body modeling. They are generally determined using seismic data or other geological data. Based on seed points, the spatial distribution of sand bodies and oil and gas reservoirs can be determined.
[0035] In this process, the grid with the largest grid attribute value in the sand body is taken as the seed point. The neighboring grids are obtained iteratively through the seed neighborhood search algorithm. The propagation paths of the seed point in three-dimensional space are constructed according to the principle of grid attribute values from large to small.
[0036] For example, such as Figure 2 As shown, each grid cell represents its attribute value using numbers. In the initial stage, each number point represents an attribute value, and 15 is defined as a random seed point. The propagation rule is from largest to smallest attribute value, and the steps are as follows:
[0037] A searches for the neighbors of seed point (15) and sorts them by priority as (14), (12), (10), with priority value being (14); B searches for the neighbor values of (14), with neighbors being (15), (11), (8), and deletes the parent class (15), with priority value being (11); C searches for the neighbor values of (11), with neighbors being (14), (12), (7), and deletes the parent class (14), and deletes duplicates of (12), with priority value being (7) and stops at the boundary.
[0038] Find the neighbors of B's second priority value (8), (14), (7), delete the parent class (14), delete duplicates of (7), and stop when priority value (8) reaches the boundary. The above steps are completed. Figure 2 During propagation, the red path is searched. After traversing the path, the neighbors with the second and third priority values in A can be searched to obtain the blue and green paths, respectively. Therefore, the red, blue, and green paths can be used as the propagation paths of the seed point.
[0039] S102. Obtain the sparsity value of the propagation path based on the backpropagation algorithm, and adjust the density of the propagation path by constraining the sparsity value;
[0040] The backpropagation algorithm described is a method of backpropagating to the endpoint of the propagation path. Backpropagation determines the number of branches along the path and thus the sparsity of the propagation path. The more branches a path has, the higher its sparsity. By constraining the sparsity, the number of branches can be reduced, thus controlling the density of the propagation path.
[0041] The preset sparsity value of the propagation path is used to truncate branch paths in the propagation path that are smaller than the preset sparsity value.
[0042] By setting a sparsity value, paths in the propagation path that are smaller than the preset sparsity value can be truncated, thereby reducing the density of the propagation path.
[0043] For example, such as Figure 3 As shown, for the propagation path, it can be seen through back propagation 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), with a density of 1 for both; grid point 14 has subclasses (8), (7), and (11), with a density of 3; seed point (15) is the initial parent class, and seed is the seed point, with a density value of the sum of the path density values of each subclass, which is 7.
[0044] Based on the calculation of path sparsity values, constraints can be imposed on propagation paths, such as path sparsity density constraints. Figure 4 As shown, when the sparsity value is set to zero, the propagation path is... Figure 4 When the sparsity value is set to 1, the propagation path is: Figure 4 b.
[0045] S103. After adjusting the density of the propagation path, the propagation path is used as a connected skeleton. Based on the positional relationship between the skeleton branches in the connected skeleton, the skeleton branches are split or merged to obtain the topology of the connected skeleton.
[0046] The connected skeleton is also known as the propagation path. Since the propagation path is composed of multiple branch paths, the skeleton branches 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.
[0047] When skeleton branches intersect, they can be merged or split according to the positional relationship at the intersection and the level of the skeleton branches.
[0048] Specifically, the levels of skeleton branches are determined based on the grid attribute values of the branches in the connected skeleton;
[0049] When branches of different levels of skeleton intersect, if the lower-level skeleton branch is not connected to the endpoint of the higher-level skeleton branch, then the lower-level connected skeleton is split.
[0050] If the endpoints of a lower-level skeleton branch are connected to the endpoints of a higher-level skeleton branch, then the lower-level skeleton branch is merged into the higher-level skeleton branch.
[0051] However, branches with higher grid attribute values in the seed point propagation path can be set as high-level branches, while branches with lower grid attribute values can be set as low-level branches. When 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.
[0052] For example, such as Figure 5 If a lower-level skeleton branch b2 is not connected to the endpoint of a higher-level skeleton branch b1, then the lower-level skeleton branch b2 needs to be split into b3 and b4.
[0053] like Figure 6 As shown, when the endpoints of a lower-level skeleton branch b2 and a higher-level skeleton branch b1 are connected or overlap, the lower-level skeleton branch b2 needs to be merged into the higher-level skeleton branch b1. That is, when the lower-level b2 and the higher-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 lower-level b2 is indirectly merged into b1, and the shorter part can be retained as the lower-level b2.
[0054] S104. Based on the topological structure of the connected skeleton, the sand bodies in the three-dimensional geological model are edited and optimized.
[0055] Based on the topological structure of the connected skeleton, attribute editing is performed on the grids within a certain range around the connected skeleton to achieve the editing and optimization of the sand body model.
[0056] Specifically, obtain the grid attribute values in the connected skeleton and the grid attribute values around the connected skeleton;
[0057] By calculating the Euclidean distance between the grid attribute values around the connected skeleton and the grid attribute values in the connected skeleton, the expansion boundary of the connected skeleton is judged, and attribute editing is performed on the grids within the expansion boundary.
[0058] Calculate the similarity of the grid attribute values within a certain range around the connected skeleton branches. According to the similarity calculation results, expand the connected skeleton outward. When the similarity of the grid in the first range skeleton is lower than the similarity of the grid in the second range, the grid corresponding to the first range is used as the boundary, and the outward expansion is terminated. Then, reassign the values of the grids within the first range, and the grids in the second range can also be reassigned, so as to complete the distinction of the outward expansion boundary of the selected skeleton branch and achieve the editing of the corresponding original geological model.
[0059] The similarity metric is the Euclidean distance. Calculate the magnitude of the attribute vector of the outward-expanded grid and the current skeleton. The smaller the value, the more similar. As Figure 7 shown, assume that A is the grid point where the skeleton is located, and the attribute value is between 700 - 80. B is the grid point with an attribute value between 700 - 60 outside the skeleton. Perform similarity judgment. d(A,B) < d(A,C). It can be considered that the attribute value of the grid point in the C range mutates, reaching the termination condition for outward search, and it can be considered that the attribute value in the B range needs to be edited. In practice, the sand body attribute value in the C range can be modified to disconnected shale to achieve the purpose of editing and correcting the model attribute value according to the skeleton topological relationship as expected.
[0060] In this embodiment, based on the sand body attributes, three-dimensional connected skeleton extraction and model editing are performed, which can solve the limitations of traditional reservoir modeling in well - to - well connectivity editing. Determine the propagation path of the seed points through the seed neighborhood search algorithm, control the density of the propagation path based on the sparse value, and directly edit and optimize 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.
[0061] 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 according to its function and internal logic, and should not constitute any limitation to the implementation process of the embodiments of the present invention.
[0062] Figure 8A schematic diagram of a three-dimensional geological model sand body connectivity skeleton extraction and editing system provided in this embodiment of the invention. The system includes:
[0063] The path construction module 810 is used to obtain a three-dimensional geological model. Based on the sand body mesh attribute values in the three-dimensional geological model, the propagation path of the seed point in the three-dimensional space is constructed through the seed neighborhood search algorithm. The seed point is the control point for sand body modeling in the three-dimensional geological model.
[0064] The step of constructing the propagation path of the seed point in three-dimensional space using the seed neighborhood search algorithm includes:
[0065] The grid with the largest grid attribute value in the sand body is taken as the seed point. The neighboring grids are obtained iteratively through the seed neighborhood search algorithm. The propagation paths of the seed point in three-dimensional space are constructed according to the principle of grid attribute value from large to small.
[0066] The sparse constraint module 820 is used 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.
[0067] Optionally, the sparse constraint module 820 includes:
[0068] The branch truncation unit is used to preset the sparsity value of the propagation path and truncate the branch paths in the propagation path that are smaller than the preset sparsity value.
[0069] The splitting and merging module 830 is used to adjust the density of the propagation path, take the propagation path as a connected skeleton, and split or merge the skeleton branches according to the positional relationship between each skeleton branch in the connected skeleton to obtain the topology of the connected skeleton.
[0070] The splitting and merging module 830 includes:
[0071] The level setting unit is used to classify the levels of skeleton branches based on the mesh attribute values of the branches in the connected skeleton;
[0072] The split unit is used to split the lower-level connected skeleton when different levels of skeleton branches intersect, if the lower-level skeleton branch is not connected to the endpoint of the higher-level skeleton branch.
[0073] The merge unit is used to merge the lower-level skeleton branch into the higher-level skeleton branch when different levels of skeleton branches intersect. If the endpoints of the lower-level skeleton branch are connected to the higher-level skeleton branch, the lower-level skeleton branch will be merged into the higher-level skeleton branch.
[0074] The editing and optimization module 840 is used to edit and optimize sand bodies in a three-dimensional geological model based on the topology of the connected skeleton.
[0075] Specifically, this involves obtaining the mesh attribute values within the connected skeleton and the mesh attribute values surrounding the connected skeleton;
[0076] By calculating the Euclidean distance between the attribute values of the meshes surrounding the connected skeleton and the attribute values of the meshes within the connected skeleton, the extended boundary of the connected skeleton is determined, and the attributes of the meshes within the extended boundary are edited.
[0077] Those skilled in the art will understand that, for the sake of convenience and brevity, the specific working process of the system and modules described above can be referred to the corresponding process in the foregoing method embodiments, and will not be repeated here.
[0078] Figure 9 This is a schematic diagram of an electronic device according to an embodiment of the present invention. The electronic device is used for three-dimensional attribute model connectivity skeleton extraction and model editing and optimization. Figure 9 As shown, the electronic device 9 of this embodiment includes: a memory 910, a processor 920, and a system bus 930. The memory 910 includes an executable program 9101 stored thereon. As those skilled in the art will understand, Figure 9 The electronic device structure shown does not constitute a limitation on the electronic device and may include more or fewer components than shown, or combine certain components, or have different component arrangements.
[0079] The following is combined Figure 9 A detailed introduction to each component of the electronic device:
[0080] The memory 910 can be used to store software programs and modules. The processor 920 executes various functional applications and data processing of the electronic device by running the software programs and modules stored in the memory 910. The memory 910 may mainly include a program storage area and a data storage area. The program storage area may store the operating system, application programs required for at least one function (such as sound playback function, image playback function, etc.), etc.; the data storage area may store data created according to the use of the electronic device (such as cached data), etc. In addition, the memory 910 may include high-speed random access memory, and may also include non-volatile memory, such as at least one disk storage device, flash memory device, or other volatile solid-state storage device.
[0081] An executable program 9101 containing an interface generation method is stored in memory 910. This executable program 9101 can be divided into one or more modules / units, which are stored in memory 910 and executed by processor 920 to perform tasks such as extracting the connected skeleton of a sand body model and editing its attributes. Each module / unit can be a series of computer program instruction segments capable of performing specific functions, describing 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 and optimization module.
[0082] The processor 920 is the control center of the electronic device. It connects various parts of the electronic device via various interfaces and lines. By running or executing software programs and / or modules stored in the memory 910, and by calling data stored in the memory 910, it performs various functions 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 mainly handles the operating system, application programs, etc., and the modem processor mainly handles wireless communication. It is understood that the aforementioned modem processor may not be integrated into the processor 920.
[0083] The system bus 930 is used to connect various functional components within the computer, transmitting data, address, and control information. Its type can be, for example, a PCI bus, an ISA bus, or a CAN bus. Instructions from the processor 920 are transmitted to the memory 910 via the bus, and the memory 910 sends data back to the processor 920. The system bus 930 handles the data and instruction exchange 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.
[0084] In this embodiment of the invention, the executable program executed by the processing 920 included in the electronic device includes:
[0085] A three-dimensional geological model is obtained. Based on the sand body mesh attribute values in the three-dimensional geological model, a seed point propagation path in three-dimensional space is constructed through a seed neighborhood search algorithm. The seed point is the control point for sand body modeling in the three-dimensional geological model.
[0086] The sparsity value of the propagation path is obtained based on the backpropagation algorithm, and the density of the propagation path is adjusted by constraining the sparsity value.
[0087] After adjusting the density of the propagation path, the propagation path is used as a connected skeleton. Based on the positional relationship between the skeleton branches, the skeleton branches are split or merged to obtain the topology of the connected skeleton.
[0088] Based on the topological structure of the connected skeleton, the sand bodies in the three-dimensional geological model are edited and optimized.
[0089] Those skilled in the art will understand that, for the sake of convenience and brevity, the specific working processes of the systems, devices, and modules described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.
[0090] In the above embodiments, the descriptions of each embodiment have different focuses. For parts that are not described in detail or recorded in a certain embodiment, please refer to the relevant descriptions of other embodiments.
[0091] The above-described embodiments are only used to illustrate the technical solutions of the present invention, and are not intended to limit it. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate 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: A three-dimensional geological model is obtained. Based on the sand body mesh attribute values in the three-dimensional geological model, a seed point propagation path in three-dimensional space is constructed through a seed neighborhood search algorithm. The seed point is the control point for sand body modeling in the three-dimensional geological model. Among them, the sand body grid attribute values are porosity, permeability, water saturation or oil saturation; The sparsity value of the propagation path is obtained based on the backpropagation algorithm, and the density of the propagation path is adjusted by constraining the sparsity value. After adjusting the density of the propagation path, the propagation path is used as a connected skeleton. Based on the positional relationship between the skeleton branches, the skeleton branches are split or merged to obtain the topology of the connected skeleton. The step of splitting or merging skeleton branches based on the positional relationships between them to obtain the topological structure of the connected skeleton includes: Based on the grid attribute values of the branches in the connected skeleton, the skeleton branches are classified into different levels; When branches of different levels of skeleton intersect, if the lower-level skeleton branch is not connected to the endpoint of the higher-level skeleton branch, then the lower-level connected skeleton is split. If a lower-level skeleton branch is connected to the endpoint of a higher-level skeleton branch, then the lower-level skeleton branch is merged into the higher-level skeleton branch. Based on the topological structure of the connected skeleton, the sand bodies in the three-dimensional geological model are edited and optimized; Specifically, this involves obtaining the mesh attribute values within the connected skeleton and the mesh attribute values surrounding the connected skeleton; By calculating the Euclidean distance between the attribute values of the meshes surrounding the connected skeleton and the attribute values of the meshes within the connected skeleton, the extended boundary of the connected skeleton is determined, and the attributes of the meshes within the extended boundary are edited.
2. The method according to claim 1, characterized in that, The propagation path of seed points in three-dimensional space, constructed using a seed neighborhood search algorithm based on the sand body mesh 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 neighboring grids are obtained iteratively through the seed neighborhood search algorithm. The 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 adjustment of the density of the propagation path by constraining sparse values includes: The sparsity value of the preset propagation path is used to truncate branch paths in the propagation path that are smaller than the preset sparsity value.
4. A system for extracting and editing the interconnected skeleton of sand bodies in a three-dimensional geological model, characterized in that, include: The path construction module is used to obtain a three-dimensional geological model. Based on the sand body mesh attribute values in the three-dimensional geological model, the propagation path of seed points in three-dimensional space is constructed through a seed neighborhood search algorithm. The seed points are the control points for sand body modeling in the three-dimensional geological model. Among them, the sand body grid attribute values are porosity, permeability, water saturation or oil saturation; The sparse constraint module is used 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. The splitting and merging module is used to adjust the density of the propagation path, and then use the propagation path as a connected skeleton. Based on the positional relationship between the skeleton branches, the skeleton branches are split or merged to obtain the topology of the connected skeleton. The splitting and merging module includes: The level setting unit is used to classify the levels of skeleton branches based on the mesh attribute values of the branches in the connected skeleton; The split unit is used to split the lower-level connected skeleton when different levels of skeleton branches intersect, if the lower-level skeleton branch is not connected to the endpoint of the higher-level skeleton branch. The merge unit is used to merge the lower-level skeleton branch into the higher-level skeleton branch when different levels of skeleton branches intersect. If the endpoints of the lower-level skeleton branch are connected to the higher-level skeleton branch, the lower-level skeleton branch is merged into the higher-level skeleton branch. The editing and optimization module is used to edit and optimize sand bodies in a 3D geological model based on the topology of the connected skeleton. Specifically, this involves obtaining the mesh attribute values within the connected skeleton and the mesh attribute values surrounding the connected skeleton; By calculating the Euclidean distance between the attribute values of the meshes surrounding the connected skeleton and the attribute values of the meshes within the connected skeleton, the extended boundary of the connected skeleton is determined, and the attributes of the meshes within the extended boundary are edited.
5. The system according to claim 4, 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 paths in the propagation path that are smaller than the preset sparsity value.
6. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the method for extracting and editing the connected skeleton of a three-dimensional geological model sand body as described in any one of claims 1 to 3.
7. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed, it implements the steps of the method for extracting and editing the connected skeleton of a three-dimensional geological model sand body as described in any one of claims 1 to 3.