Automatic updating method for three-dimensional geologic model

Abstract

The invention provides a three-dimensional geological model automatic updating method, and the method comprises the steps: S11, delineating a local to-be-updated region of a three-dimensional geological model through a polygon according to the distribution position of newly added geological exploration data; s12, extracting geological exploration data of the local to-be-updated area, constructing a three-dimensional structure model of the local to-be-updated area, and performing mesh generation on the three-dimensional structure model to obtain a mesh model; s13, constructing a three-dimensional attribute model of the local to-be-updated area, obtaining parameters used by the three-dimensional attribute model, extracting sample points required by the three-dimensional attribute model, and performing interpolation calculation on the grid model through an interpolation method to complete construction of the three-dimensional attribute model; and S14, replacing the three-dimensional structure model and the three-dimensional attribute model obtained by construction with the current three-dimensional geologic model. According to the method, the existing three-dimensional geological model is locally updated by using the newly added geological exploration data, and the problem of synchronously updating the structure model and the attribute model is solved.

Description

Technical Field

[0001] The present invention relates to the field of geological exploration processing technology, and in particular to a three-dimensional geological model automatic updating method. Background Art

[0002] As geological exploration continues, new geological exploration data (such as drilling data and related test data) is generated. This new data can be used to update existing 3D geological models. The general approach is to merge the new data with the old data to rebuild the 3D geological model, and then replace the old model with the new model. In scenarios where the model covers a large area (such as a city-level 3D geological model), the workload of building a new model as a whole is large and time-consuming. The new data may be localized, making it unworthy to perform a comprehensive model update. In addition, the general update method has a single update target: either only updating the structural model that reflects the geological tectonic conditions, or only updating the attribute model that reflects the distribution status of geological attributes, and it is not possible to update all related models simultaneously. Summary of the Invention

[0003] The purpose of the present invention is to provide a method for automatically updating a three-dimensional geological model, which uses newly added geological exploration data to perform local updates on the current three-dimensional geological model and solves the problem of synchronously updating the structural model and the attribute model.

[0004] To achieve the above-mentioned object, the present invention provides a method for automatically updating a three-dimensional geological model, the method comprising: S11, based on the distribution of the newly added geological exploration data, delineate the local area to be updated of the three-dimensional geological model using polygons; S12, extracting geological exploration data of the local area to be updated, constructing a three-dimensional structural model of the local area to be updated, and performing grid division on the three-dimensional structural model to obtain a grid model; S13, constructing a three-dimensional attribute model of the local area to be updated, obtaining parameters used by the three-dimensional attribute model, extracting sample points required for the three-dimensional attribute model, and performing interpolation calculation on the grid model by an interpolation method to complete the construction of the three-dimensional attribute model; S14. Replace the constructed three-dimensional structural model and three-dimensional attribute model with the existing three-dimensional geological model.

[0005] Furthermore, step S11 specifically includes: The range affected by the newly added geological exploration data is automatically delineated in the form of polygons on the 3D geological model as the local area to be updated in the 3D geological model.

[0006] Furthermore, in step S12, the geological exploration data of the local area to be updated is extracted, which specifically includes: Extract all geological exploration data in the local area to be updated, including new data and historical data; Extract virtual geological exploration data along the boundary of the local area to be updated.

[0007] Furthermore, in step S12, a three-dimensional structural model of the local area is constructed, specifically including: The boundaries of the polygons in the local area to be updated are used as boundary constraints of the three-dimensional structural model, and the extracted geological exploration data are used as input data of the three-dimensional structural model to construct a three-dimensional structural model. The three-dimensional structural model is composed of several geological bodies, which are expressed as vector data. Each geological body records the code of the stratum to which it belongs.

[0008] Furthermore, in step S14, the 3D structural model is replaced with the existing 3D geological model, specifically including: The geological body expressed by vector data in the three-dimensional structural model is decomposed into vector grids, and the vector grids are replaced one by one with the grids at the corresponding positions in the three-dimensional geological model.

[0009] Furthermore, in step S12, the three-dimensional structural model is meshed, specifically including: The three-dimensional structural model is meshed using the division parameters consistent with the current three-dimensional geological model to obtain a grid model composed of rectangular grids, and each grid unit records the stratigraphic code to which it belongs.

[0010] Furthermore, in step S13, the parameters used by the three-dimensional attribute model are obtained, specifically including: According to the target attributes of the three-dimensional attribute model and the stratigraphic code of the grid model, the kriging parameters used to perform interpolation for each target attribute are obtained.

[0011] Furthermore, in step S13, the sample points required for the three-dimensional attribute model of the local area to be updated are extracted, which specifically includes: Take the maximum search ellipsoid radius of the Kriging parameters used in the three-dimensional attribute model, denoted as R; The local area to be updated is denoted as P'. A buffer zone with a distance R is added to the periphery of the local area to be updated to obtain a new area, denoted as P. New sample points and existing sample points located in the region P are extracted as sample points of the three-dimensional attribute model.

[0012] Furthermore, in step S13, interpolation calculation is performed on the grid model by an interpolation method, specifically including: By using the obtained Kriging parameters and extracted sample points, the stratum to which the grid model belongs is traversed, and Kriging interpolation is performed on the grid models respectively to complete the construction of the three-dimensional attribute model.

[0013] Furthermore, in step S14, the 3D attribute model of the local area to be updated is replaced with the existing 3D geological model, which specifically includes: The grid in the 3D attribute model is directly replaced with the grid at the corresponding position in the existing 3D geological model.

[0014] Compared with the current technology, the beneficial effects of the present invention are: The present invention provides a method for automatically updating a 3D geological model. By first dividing the area to be updated, the structural and attribute models are updated only in that local area, rather than updating the entire geological model comprehensively. This reduces computational and data processing workload. By breaking the update process into two steps—updating the structural model and updating the attribute model—this localized update approach allows for more targeted processing of the changed areas. The updated 3D structural and attribute models are then fully replaced with the original 3D geological model, ensuring continuity and consistency between the updated and unupdated areas. Compared to traditional methods, this method offers the ability to perform local model updates on demand, reduces workload, improves efficiency, and supports the simultaneous updating of associated structural and attribute models. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the current technology, the following is a brief introduction to the drawings required for use in the embodiments or the current technology description. Obviously, the drawings described below are merely embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without any creative work. Figure 1 A schematic flow chart of a method for automatically updating a three-dimensional geological model provided by an embodiment of the present invention; Figure 2 A schematic diagram of a two-dimensional section decomposition into a vector grid in a three-dimensional geological model automatic updating method provided by an embodiment of the present invention; Figure 3 A schematic diagram of an example of vector grid edge comparison in a three-dimensional geological model automatic updating method provided by an embodiment of the present invention; Figure 4 A schematic diagram of sample points required for extracting a three-dimensional attribute model is provided in an automatic updating method of a three-dimensional geological model provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0016] The present invention will be further described in detail below with reference to the accompanying drawings and examples. It will be understood that the specific embodiments described herein are intended only to illustrate the present invention and are not intended to limit the present invention. It should also be noted that, for ease of description, the accompanying drawings only illustrate portions relevant to the present invention, not all structures.

[0017] Reference Figure 1 This embodiment provides a method for automatically updating a three-dimensional geological model, the method comprising: S11. Based on the distribution of the newly added geological exploration data, a local area to be updated of the three-dimensional geological model is delineated by polygons.

[0018] S12. Extracting geological exploration data of the local area to be updated, constructing a three-dimensional structural model of the local area to be updated, and performing grid division on the three-dimensional structural model to obtain a grid model.

[0019] S13, constructing a three-dimensional attribute model of the local area to be updated, obtaining parameters used by the three-dimensional attribute model, extracting sample points required for the three-dimensional attribute model, performing interpolation calculation on the grid model through interpolation method, and completing the construction of the three-dimensional attribute model.

[0020] S14. Replace the constructed three-dimensional structural model and three-dimensional attribute model with the existing three-dimensional geological model.

[0021] In this example, newly added geological exploration data (drilling data and related test data) is used to locally update the existing 3D geological model and simultaneously update the structural and attribute models. This method assumes that the global stratigraphic sequence of the region is well defined—that is, the total number of strata and the relationship between the old and new strata are well defined—and that the newly added data conforms to this stratigraphic division scheme.

[0022] As a preferred embodiment, step S11 specifically includes: The range affected by the newly added geological exploration data is automatically delineated in the form of polygons on the 3D geological model as the local area to be updated in the 3D geological model.

[0023] In this embodiment, the area to be updated in the 3D geological model is delineated based on the distribution of the newly added geological exploration data. The system displays the corresponding location of the newly added geological exploration data on the 3D geological model. The software automatically or the user manually delineates the area affected by the new data as the local area to be updated in the 3D geological model, representing this area with a polygon.

[0024] As a preferred embodiment, in step S12, extracting geological exploration data of the local area to be updated specifically includes: Extract all geological exploration data in the local area to be updated, including new data and historical data; Extract virtual geological exploration data along the boundary of the local area to be updated.

[0025] In this embodiment, geological exploration data is extracted for local updates of the 3D geological model. This geological exploration data consists of two parts: First, all geological exploration data (including new and historical data) within the output polygon is extracted from the system library; second, virtual geological exploration data is extracted from the existing 3D geological model along the polygon's boundaries. This virtual geological exploration data ensures that the new 3D geological model maintains a good fit with the old 3D geological model at the boundaries. When extracting virtual geological exploration data along the boundaries, a spacing must be specified, which can be equal to the side length of the mesh used during modeling.

[0026] As a preferred embodiment, in step S12, constructing a three-dimensional structural model of the local area specifically includes: The boundaries of the polygons in the local area to be updated are used as boundary constraints of the three-dimensional structural model, and the extracted geological exploration data are used as input data of the three-dimensional structural model to construct a three-dimensional structural model. The three-dimensional structural model is composed of several geological bodies, which are expressed as vector data. Each geological body records the code of the stratum to which it belongs.

[0027] In this embodiment, the data obtained in the above steps is used to construct a three-dimensional geological model of a local area: a structural model. The polygons constructed above serve as the modeling boundary constraints, and the geological exploration data extracted for local updates of the three-dimensional geological model serves as the modeling input data. The modeling method is an automatic modeling method based on geological exploration data provided by the software system. The output is a geological body expressed as vector data. Several geological bodies constitute the three-dimensional structural model, and each geological body records the code of the stratum to which it belongs.

[0028] As a preferred embodiment, in step S14, replacing the 3D structural model with the existing 3D geological model specifically includes: The geological body expressed by vector data in the three-dimensional structural model is decomposed into vector grids, and the vector grids are replaced one by one with the grids at the corresponding positions in the three-dimensional geological model.

[0029] In this embodiment, the new local area 3D geological model is integrated (or replaced) into the existing 3D geological model to complete the 3D geological model update. The 3D structural model and the 3D attribute model are integrated using different methods. The method of structural model integration is as follows: the geological body expressed as vector data in the new 3D structural model is decomposed into "vector grids". A vector grid is a "pillar" with a square cross section. Decomposing the 3D structural model into vector grids is to divide it into a set of closely adjacent "pillars" (see Figure 2Schematic diagram of decomposing a 2D section into a vector grid. Each "column" in a two-dimensional environment has two edges, while in a 3D environment, each "column" has four edges. These edges record the locations of stratigraphic divisions, effectively creating virtual geological exploration data. Because the stratigraphic sequence is fixed, the number and order of stratigraphic layers recorded on each edge are identical; the only difference is the location of the layers. When decomposing the vector grid, parameters consistent with the existing 3D geological model (including mesh cross-section side lengths and mesh reference points) must be used.

[0030] Use vector grids to replace the grids at corresponding positions in the existing three-dimensional geological model one by one to complete the update of the three-dimensional geological model. During the replacement, there may be inconsistencies in the stratigraphic division positions recorded on the edges of adjacent grids, such as stratigraphic pinch-out (stratum thickness is 0) and deviations in stratigraphic stratification positions. On a two-dimensional plane, each edge is shared by two adjacent columns; in a three-dimensional environment, each edge is shared by four adjacent columns. Under the condition that the stratigraphic sequence is determined, compare the elevations of the stratigraphic stratification points in the new and old edges. If there is a deviation, take the average elevation of the corresponding two points as the final stratification position to ensure data consistency between the new and old grids (such as Figure 3 ). Only the grids at the boundaries of the new 3D geological model need to be compared and calculated with the old grid during the replacement process. The grids inside the new 3D geological model can directly replace the old grid.

[0031] As a preferred embodiment, in step S12, meshing the three-dimensional structure model specifically includes: The three-dimensional structural model is meshed using the division parameters consistent with the current three-dimensional geological model to obtain a grid model composed of rectangular grids, and each grid unit records the stratigraphic code to which it belongs.

[0032] In this embodiment, a grid model consisting of rectangular grid cells is generated on the 3D structural model (also known as a regular grid model). Each grid cell also records the stratigraphic code to which it belongs. The parameters used for gridding (including grid reference points and grid cell size) are consistent with those of the existing 3D geological model. This grid model serves as the basis for subsequent updates to the 3D geological model.

[0033] As a preferred embodiment, in step S13, obtaining parameters used by the three-dimensional attribute model specifically includes: According to the target attributes of the three-dimensional attribute model and the stratum code of the grid model, the kriging parameters used to perform interpolation for each target attribute are obtained.

[0034] In this embodiment, the parameters used for 3D attribute modeling are queried. The 3D attribute model is built using the Kriging interpolation algorithm. The system records the Kriging parameters used for interpolation of each attribute. If necessary, different parameters can be used for different strata. Therefore, the grid model needs to record the stratum code to which each grid belongs. Based on the target attribute and stratum code, the corresponding Kriging parameters can be queried from the system.

[0035] As a preferred embodiment, in step S13, extracting sample points required for the three-dimensional attribute model of the local area to be updated specifically includes: The maximum search ellipsoid radius of the Kriging parameters used in the three-dimensional attribute model is taken, denoted as R.

[0036] The local area to be updated is denoted as P', and a buffer zone with a distance R is added to the periphery of the local area to be updated to obtain a new area, which is denoted as P.

[0037] New sample points and existing sample points located in the region P are extracted as sample points of the three-dimensional attribute model.

[0038] In this embodiment, the sample points required for updating the three-dimensional attribute model of the local area to be updated are extracted. The maximum search ellipsoid radius (denoted as R) of the Kriging parameters used in the three-dimensional attribute model is taken, the local area to be updated is denoted as P', and a buffer zone with a distance of R is added to the periphery to obtain a new area (denoted as P). Afterwards, new sample points and existing sample points located in area P are extracted as sample points for subsequent updating of the attribute model. The purpose of adding a buffer zone with a distance of the maximum search ellipsoid radius to extract sample points is to ensure that the Kriging interpolation results at the boundary of the update area can maintain continuity with the existing model. Figure 4 , taking P' as a rectangle as an example, a buffer with a distance R is added to obtain a new polygon P, which is used to extract sample points for interpolation of the three-dimensional attribute model.

[0039] As a preferred embodiment, in step S13, interpolation calculation is performed on the grid model by an interpolation method, specifically including: By using the obtained Kriging parameters and extracted sample points, the stratum to which the grid model belongs is traversed, and Kriging interpolation is performed on the grid models respectively to complete the construction of the three-dimensional attribute model.

[0040] In this embodiment, interpolation is performed to calculate a 3D geological model of the local area to be updated: a 3D attribute model. Using the obtained sample points and the corresponding kriging parameters, kriging interpolation is performed on the grid model to obtain an updated 3D attribute model. Different strata have different attribute interpolation parameters, so this step involves traversing each stratum to which the grid belongs and performing interpolation on each stratum.

[0041] As a preferred embodiment, in step S14, replacing the 3D attribute model of the local area to be updated with the existing 3D geological model specifically includes: The grid in the 3D attribute model is directly replaced with the grid at the corresponding position in the existing 3D geological model.

[0042] In this embodiment, the new local 3D geological model is integrated (or replaced) into the existing 3D geological model to complete the 3D geological model update. The 3D structural model and 3D attribute model are integrated using different methods. 3D attribute model integration simply replaces the meshes at corresponding locations in the existing 3D geological model with those in the new 3D geological model.

[0043] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for automatically updating a three-dimensional geological model, characterized in that: The method comprises: S11, based on the distribution of the newly added geological exploration data, delineate the local area to be updated of the three-dimensional geological model using polygons; S12, extracting geological exploration data of the local area to be updated, constructing a three-dimensional structural model of the local area to be updated, and performing grid division on the three-dimensional structural model to obtain a grid model; S13, constructing a three-dimensional attribute model of the local area to be updated, obtaining parameters used by the three-dimensional attribute model, extracting sample points required for the three-dimensional attribute model, and performing interpolation calculation on the grid model by an interpolation method to complete the construction of the three-dimensional attribute model; S14. Replace the constructed three-dimensional structural model and three-dimensional attribute model with the existing three-dimensional geological model.

2. The automatic updating method of a three-dimensional geological model according to claim 1, characterized in that: Step S11 specifically includes: The range affected by the newly added geological exploration data is automatically delineated in the form of polygons on the 3D geological model as the local area to be updated in the 3D geological model.

3. The automatic updating method of a three-dimensional geological model according to claim 1, characterized in that: In step S12, the geological exploration data of the local area to be updated is extracted, which specifically includes: Extract all geological exploration data in the local area to be updated, including new data and historical data; Extract virtual geological exploration data along the boundary of the local area to be updated.

4. The automatic updating method of a three-dimensional geological model according to claim 1, characterized in that: In step S12, a three-dimensional structural model of the local area is constructed, which specifically includes: The boundaries of the polygons in the local area to be updated are used as boundary constraints of the three-dimensional structural model, and the extracted geological exploration data are used as input data of the three-dimensional structural model to construct a three-dimensional structural model. The three-dimensional structural model is composed of several geological bodies, which are expressed as vector data. Each geological body records the code of the stratum to which it belongs.

5. The automatic updating method of a three-dimensional geological model according to claim 4, characterized in that: In step S14, the 3D structural model is replaced with the existing 3D geological model, specifically including: The geological body expressed by vector data in the three-dimensional structural model is decomposed into vector grids, and the vector grids are replaced one by one with the grids at the corresponding positions in the three-dimensional geological model.

6. The method for automatically updating a three-dimensional geological model according to claim 1, characterized in that: In step S12, the three-dimensional structural model is meshed, specifically including: The three-dimensional structural model is meshed using the division parameters consistent with the current three-dimensional geological model to obtain a grid model composed of rectangular grids, and each grid unit records the stratigraphic code to which it belongs.

7. The automatic updating method of a three-dimensional geological model according to claim 6, characterized in that: In step S13, the parameters used by the three-dimensional attribute model are obtained, specifically including: According to the target attributes of the three-dimensional attribute model and the stratigraphic code of the grid model, the kriging parameters used to perform interpolation for each target attribute are obtained.

8. The method for automatically updating a three-dimensional geological model according to claim 7, characterized in that: In step S13, the sample points required for the three-dimensional attribute model of the local area to be updated are extracted, which specifically includes: Take the maximum search ellipsoid radius of the Kriging parameters used in the three-dimensional attribute model, denoted as R; The local area to be updated is denoted as P'. A buffer zone with a distance R is added to the periphery of the local area to be updated to obtain a new area, denoted as P. New sample points and existing sample points located in the region P are extracted as sample points of the three-dimensional attribute model.

9. The method for automatically updating a three-dimensional geological model according to claim 8, characterized in that: In step S13, interpolation calculation is performed on the grid model by using an interpolation method, specifically including: By using the obtained Kriging parameters and extracted sample points, the stratum to which the grid model belongs is traversed, and Kriging interpolation is performed on the grid models respectively to complete the construction of the three-dimensional attribute model.

10. The method for automatically updating a three-dimensional geological model according to claim 1, characterized in that: In step S14, the 3D attribute model of the local area to be updated is replaced with the existing 3D geological model, which specifically includes: The grid in the 3D attribute model is directly replaced with the grid at the corresponding position in the existing 3D geological model.

Images