A method for automatically updating a three-dimensional geological model

Abstract

The application provides a kind of three-dimensional geological model automatic updating method, the method comprises: S11, according to the distribution position of new geological exploration data, the local area to be updated of three-dimensional geological model is circled by polygon;S12, extract the geological exploration data of local area to be updated, construct the three-dimensional structure model of local area to be updated, grid subdivision is carried out to three-dimensional structure model, and grid model is obtained;S13, construct the three-dimensional attribute model of local area to be updated, obtain the parameter used by three-dimensional attribute model, extract the sample point required by three-dimensional attribute model, carry out interpolation calculation to grid model by interpolation method, and the construction of three-dimensional attribute model is completed;S14, the three-dimensional structure model and three-dimensional attribute model obtained by construction are replaced to current three-dimensional geological model.The application carries out local update to existing three-dimensional geological model using new geological exploration data, and solves the problem of synchronous update structure model and attribute model.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of geological exploration processing, and particularly relates to a three-dimensional geological model automatic updating method. BACKGROUND

[0002] With the continuous development of geological exploration work, new geological exploration data (such as drilling data and related test data) will be generated. These new data can be used to update the built three-dimensional geological model. The general method is to merge the new data with the old data and then reconstruct the three-dimensional geological model, and then replace the old model with the new model. In the scenario where the model covers a large area (such as a city-level three-dimensional geological model), the workload of building a new model is large and time-consuming. The new data may be in a local area, and it is not worth performing a whole model update. In addition, the update target of the general method is single: either only the structure model reflecting the geological structure condition is updated, or only the attribute model reflecting the geological attribute distribution state is updated, and all related models cannot be updated synchronously. SUMMARY

[0003] The present application aims to provide a three-dimensional geological model automatic updating method, which uses new geological exploration data to perform local update on the current three-dimensional geological model and solves the problem of synchronous update of structure model and attribute model.

[0004] To achieve the above-mentioned application purpose, the present application provides a three-dimensional geological model automatic updating method, which comprises the following steps:

[0005] S11. According to the distribution position of the new geological exploration data, a local area to be updated of the three-dimensional geological model is circled by a polygon;

[0006] S12. The geological exploration data of the local area to be updated is extracted, a three-dimensional structure model of the local area to be updated is constructed, the three-dimensional structure model is meshed, and a mesh model is obtained;

[0007] S13. A three-dimensional attribute model of the local area to be updated is constructed, parameters used by the three-dimensional attribute model are obtained, sample points required by the three-dimensional attribute model are extracted, the mesh model is calculated by interpolation, and the construction of the three-dimensional attribute model is completed;

[0008] S14. The three-dimensional structure model and the three-dimensional attribute model obtained by construction are replaced into the existing three-dimensional geological model.

[0009] Further, in step S11, the following steps are specifically included:

[0010] The range affected by the new geological exploration data is automatically circled in the form of a polygon on the three-dimensional geological model, which is the local area to be updated of the three-dimensional geological model.

[0011] Further, in step S12, the geological exploration data in the local region to be updated is extracted, specifically including:

[0012] All geological exploration data in the local region to be updated, including new data and historical data, is extracted.

[0013] Virtual geological exploration data along the boundary of the local region to be updated is extracted.

[0014] Further, in step S12, a three-dimensional structure model of the local region is constructed, specifically including:

[0015] The boundary of the polygon in the local region to be updated is taken as the boundary constraint of the three-dimensional structure model, and the extracted geological exploration data is taken as the input data of the three-dimensional structure model, and the three-dimensional structure model is constructed, the three-dimensional structure model is composed of a plurality of geological bodies, and each geological body records the code of the stratum to which it belongs.

[0016] Further, in step S14, the three-dimensional structure model is replaced into the existing three-dimensional geological model, specifically including:

[0017] The geological body expressed by vector data in the three-dimensional structure model is decomposed into a vector grid, and the vector grid is replaced into the corresponding grid in the three-dimensional geological model one by one.

[0018] Further, in step S12, the three-dimensional structure model is meshed, specifically including:

[0019] The three-dimensional structure model is meshed using the same division parameters as the current three-dimensional geological model, and a grid model composed of cuboid grids is obtained, and each grid unit records the stratum code to which it belongs.

[0020] Further, in step S13, the parameters used by the three-dimensional attribute model are obtained, specifically including:

[0021] According to the target attribute of the three-dimensional attribute model and the stratum code of the grid model, the Kriging parameters used for each target attribute to perform interpolation are obtained.

[0022] Further, in step S13, the sample points required by the three-dimensional attribute model in the local region to be updated are extracted, specifically including:

[0023] Take the maximum search ellipsoid radius of the Kriging parameters used by the three-dimensional attribute model, denoted as R;

[0024] The local region to be updated is denoted as P', and a buffer zone with a distance of R is added to the periphery of the local region to be updated to obtain a new region, denoted as P;

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

[0026] Further, in step S13, the grid model is calculated by interpolation method, specifically including:

[0027] By obtaining the Kriging parameters and extracting the sample points, the grid model belonging to the stratum is traversed, and the Kriging interpolation is performed on the grid model respectively, so as to complete the construction of the three-dimensional attribute model.

[0028] Further, in step S14, the three-dimensional attribute model of the local to-be-updated region is replaced to the existing three-dimensional geological model, specifically including:

[0029] The grid in the three-dimensional attribute model is directly replaced to the grid at the corresponding position in the existing three-dimensional geological model.

[0030] Compared with the prior art, the beneficial effects of the present application are:

[0031] The three-dimensional geological model automatic updating method provided by the present application divides the to-be-updated region first, and only the structure model and the attribute model in the local region are updated, instead of comprehensively updating the entire geological model, so that the calculation amount and the data processing workload are reduced; the updating process is divided into two steps of structure model updating and attribute model updating, and the local updating manner can more targetedly process the changed region; the updated three-dimensional structure model and attribute model are replaced to the original three-dimensional geological model as a whole, so that the continuity and consistency of the three-dimensional geological model between the updated region and the non-updated region are ensured. Compared with the traditional method, the present method has the capability of performing local model updating on demand, has small workload and high efficiency, and supports synchronous updating of the associated structure model and attribute model. BRIEF DESCRIPTION OF DRAWINGS

[0032] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the drawings needed to be used in the embodiments or the prior art description will be briefly introduced. Obviously, the drawings in the following description are only embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor on the basis of the provided drawings.

[0033] Figure 1 A three-dimensional geological model automatic updating method flowchart provided by the embodiment of the present application;

[0034] Figure 2 A two-dimensional profile decomposition into a vector grid diagram of the three-dimensional geological model automatic updating method provided by the embodiment of the present application;

[0035] Figure 3 A vector grid edge ratio comparison example diagram of the three-dimensional geological model automatic updating method provided by the embodiment of the present application;

[0036] Figure 4 A three-dimensional geological model automatic updating method provided by an embodiment of the present application extracts sample points required by a three-dimensional attribute model. DETAILED DESCRIPTION

[0037] The present application will be further described below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely intended for the purpose of interpretation of the present application and are not limiting of the present application. In addition, it should be noted that only parts related to the present application are shown in the accompanying drawings for the purpose of description.

[0038] Reference Figure 1 The present embodiment provides a three-dimensional geological model automatic updating method, which comprises the following steps:

[0039] S11. According to the distribution position of the newly added geological exploration data, a local region to be updated of the three-dimensional geological model is circled by a polygon.

[0040] S12. Geological exploration data of the local region to be updated is extracted, a three-dimensional structure model of the local region to be updated is constructed, the three-dimensional structure model is meshed, and a mesh model is obtained.

[0041] S13. A three-dimensional attribute model of the local region to be updated is constructed, parameters used by the three-dimensional attribute model are obtained, sample points required by the three-dimensional attribute model are extracted, the mesh model is calculated by an interpolation method, and the construction of the three-dimensional attribute model is completed.

[0042] S14. The three-dimensional structure model and the three-dimensional attribute model obtained by the construction are replaced into an existing three-dimensional geological model.

[0043] In the present embodiment, the newly added geological exploration data (drilling data and related test data) is used to perform local updating on the existing three-dimensional geological model, and the problem of synchronous updating of the structure model and the attribute model is solved. The precondition for using the method is that the global stratigraphic sequence of the region is determined, i.e., which strata are divided in total, the new and old coverage relationship of the strata is determined, and the newly added data also conforms to the stratigraphic division scheme.

[0044] As a preferred embodiment, step S11 specifically comprises the following steps:

[0045] The range affected by the newly added geological exploration data is automatically circled in the form of a polygon on the three-dimensional geological model, serving as the local region to be updated of the three-dimensional geological model.

[0046] In this embodiment, the region to be updated of the three-dimensional geological model is circled according to the distribution position of the newly added geological exploration data. The system displays the corresponding position of the newly added geological exploration data on the three-dimensional geological model, and the range affected by the new data is automatically circled by the software system or manually circled by the user, as a local region to be updated of the three-dimensional geological model, which is expressed by a polygon.

[0047] As a preferred embodiment, in step S12, the geological exploration data of the local region to be updated is extracted, specifically including:

[0048] All the geological exploration data in the local region to be updated, including the newly added data and the historical data, is extracted.

[0049] The virtual geological exploration data along the boundary of the local region to be updated is extracted.

[0050] In this embodiment, the geological exploration data used for the local update of the three-dimensional geological model is extracted. The geological exploration data is composed of two parts. One is all the geological exploration data (including the newly added data and the historical data) in the polygon output from the system library; the other is the virtual geological exploration data extracted from the existing three-dimensional geological model along the boundary of the polygon. These virtual geological exploration data can control the new three-dimensional geological model to maintain good data consistency with the old three-dimensional geological model at the boundary. When the virtual geological exploration data is extracted along the boundary, the interval needs to be specified, and the value can be equal to the side length of the grid used for modeling.

[0051] As a preferred embodiment, in step S12, the three-dimensional structure model of the local region is constructed, specifically including:

[0052] The boundary of the polygon in the local region to be updated is taken as the boundary constraint of the three-dimensional structure model, and the extracted geological exploration data is taken as the input data of the three-dimensional structure model, so as to construct the three-dimensional structure model. The three-dimensional structure model is composed of a plurality of geological bodies, and each geological body records the code of the stratum to which it belongs.

[0053] In this embodiment, the three-dimensional geological model of the local region, i.e., the structure model, is constructed by using the data obtained by the above steps. The polygon constructed above is the boundary constraint of modeling, and the geological exploration data extracted for the local update of the three-dimensional geological model is the input data of modeling. The modeling method is an automatic modeling method based on geological exploration data provided by the software system, and the output is a geological body expressed by vector data. A plurality of geological bodies constitute the three-dimensional structure model, and each geological body records the code of the stratum to which it belongs.

[0054] As a preferred embodiment, in step S14, the three-dimensional structure model is replaced to the existing three-dimensional geological model, specifically including:

[0055] The geological body expressed by vector data in the three-dimensional structure model is decomposed into vector grids, and the vector grids are replaced one by one into the grids at the corresponding positions in the three-dimensional geological model.

[0056] In this embodiment, the new local area three-dimensional geological model is integrated (or replaced) into the existing three-dimensional geological model, and the three-dimensional geological model is updated. The three-dimensional structure model and the three-dimensional attribute model are integrated using different methods. The method of integrating the structure model is as follows: The geological body expressed by vector data in the new three-dimensional structure model is decomposed into "vector grids". A vector grid is a "column" with a square cross section. Decomposing the three-dimensional structure model into vector grids is to divide it into a set of closely connected "columns" (see Figure 2 The schematic diagram of decomposing a two-dimensional profile into a vector grid). Each "column" on the plane has 2 edges, and each "column" in the three-dimensional environment has 4 edges. The edges record the positions of the stratigraphic division, which is equivalent to a virtual geological exploration data. Because the stratigraphic sequence is determined, the number and sequence of the stratigraphic layers recorded on each edge are the same, and the only difference is the position of the stratigraphic layering. When decomposing the vector grid, the same parameters as the existing three-dimensional geological model should be used (including the edge length of the grid cross section, the reference point of the grid subdivision, etc.).

[0057] The vector grids are replaced one by one into the grids at the corresponding positions in the existing three-dimensional geological model, and the three-dimensional geological model is updated. When replacing, there may be a situation that the positions of the stratigraphic division recorded on the edges of the adjacent grids are inconsistent, such as stratigraphic pinchout (stratigraphic thickness is 0) and stratigraphic layering position deviation. In the two-dimensional plane, each edge is shared by two adjacent columns; in the three-dimensional environment, each edge is shared by four adjacent columns. Under the condition that the stratigraphic sequence is determined, by comparing the elevations of the stratigraphic layering points in the new edge and the old edge, when there is a deviation, the average elevation of the corresponding two points is taken as the final stratigraphic layering position, to ensure the data consistency of the new and old grids (such as Figure 3 Only the grids at the boundary of the new three-dimensional geological model need to perform the operations of comparison and calculation with the old grids in the replacement process, and the grids inside the new three-dimensional geological model can directly replace the old grids.

[0058] As a preferred embodiment, in step S12, the three-dimensional structure model is subjected to grid subdivision, specifically including:

[0059] The three-dimensional structure model is subjected to grid subdivision using the same division parameters as the current three-dimensional geological model, to obtain a grid model composed of cuboid grids, and each grid unit records the stratigraphic code to which it belongs.

[0060] In this embodiment, the grid subdivision is performed on the three-dimensional structure model to obtain a grid model composed of cuboid grids (also referred to as regular grid model). Each grid cell also records the stratum code to which it belongs. The parameters used in the grid subdivision (including the reference point of the subdivision, the size of the grid cell, etc.) are consistent with the existing three-dimensional geological model. The grid model is the basis for subsequent three-dimensional geological model updating.

[0061] As a preferred embodiment, in step S13, the parameters used in the three-dimensional attribute model are obtained, specifically including:

[0062] According to the target attribute of the three-dimensional attribute model and the stratum code of the grid model, the Kriging parameters used in the interpolation of each target attribute are obtained.

[0063] In this embodiment, the parameters used in the three-dimensional attribute modeling are queried. The Kriging interpolation algorithm is used in the three-dimensional attribute modeling. The system records the Kriging parameters used in the interpolation of each attribute in advance. If necessary, different parameters can also be used in different strata. Therefore, the grid model needs to record the stratum code of each grid. According to the target attribute and the stratum code, the corresponding Kriging parameters can be queried from the system.

[0064] As a preferred embodiment, in step S13, the sample points required by the three-dimensional attribute model in the local region to be updated are extracted, specifically including:

[0065] The maximum search ellipsoid radius of the Kriging parameters used in the three-dimensional attribute model is taken as R.

[0066] The local region 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 region, denoted as P.

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

[0068] In this embodiment, the sample points required by the three-dimensional attribute model in the local region 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 region 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 region (denoted as P). Then, the new sample points and the existing sample points located in the region P are extracted as the sample points of the subsequent attribute model. The purpose of adding the buffer zone with a distance of the maximum search ellipsoid radius to extract the sample points is to make the Kriging interpolation result at the boundary of the updated region continuous with the existing model. Figure 4 Taking the rectangle P' as an example, a new polygon P with a distance of R is added to the buffer zone to extract the sample points for the Kriging interpolation of the three-dimensional attribute model.

[0069] As a preferred embodiment, in step S13, the grid model is calculated by interpolation, specifically including:

[0070] By obtaining the Kriging parameters and the extracted sample points, the grid model is traversed to the formation to which the grid model belongs, and Kriging interpolation is performed on the grid model, thereby completing the construction of the three-dimensional attribute model.

[0071] In this embodiment, the three-dimensional geological model of the local region to be updated is calculated by interpolation: the three-dimensional attribute model. The obtained sample points and the corresponding Kriging parameters are used to perform Kriging interpolation on the grid model, thereby obtaining the updated three-dimensional attribute model. Different formations have different attribute interpolation parameters, so this step needs to traverse the formation to which the grid belongs and perform interpolation respectively.

[0072] As a preferred embodiment, in step S14, the three-dimensional attribute model of the local region to be updated is replaced into the existing three-dimensional geological model, specifically including:

[0073] The grid in the three-dimensional attribute model is directly replaced into the grid at the corresponding position in the existing three-dimensional geological model.

[0074] In this embodiment, the new local region three-dimensional geological model is integrated (or replaced) into the existing three-dimensional geological model, thereby completing the three-dimensional geological model update. The three-dimensional structure model and the three-dimensional attribute model use different methods to perform integration. The three-dimensional attribute model integration only needs to replace the grid in the existing three-dimensional geological model with the grid in the new three-dimensional geological model.

[0075] The above only describes the preferred embodiments of the present application and is not intended to limit the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application shall be included in the protection scope of the present application.

Claims

1. A method for automatically updating a three-dimensional geological model, characterized in that, The method comprises: S11, according to the distribution position of the newly added geological exploration data, a local area to be updated of the three-dimensional geological model is circled by a polygon; S12, geological exploration data of the local area to be updated is extracted, and a three-dimensional structure model of the local area to be updated is constructed, and the three-dimensional structure model is meshed to obtain a grid model; The three-dimensional structure model is meshed, specifically comprising: Using the same division parameters as the current three-dimensional geological model, the three-dimensional structure model is meshed to obtain a grid model composed of cuboid grids, and each grid unit records the stratigraphic code to which it belongs; S13, a three-dimensional attribute model of the local area to be updated is constructed, parameters used by the three-dimensional attribute model are obtained, sample points required by the three-dimensional attribute model are extracted, and the grid model is calculated by interpolation to complete the construction of the three-dimensional attribute model; The parameters used by the three-dimensional attribute model are obtained, specifically comprising: According to the target attribute of the three-dimensional attribute model and the stratigraphic code of the grid model, the Kriging parameters used by each target attribute for interpolation are obtained; The sample points required by the three-dimensional attribute model of the local area to be updated are extracted, specifically comprising: Taking the maximum search ellipsoid radius of the Kriging parameters used by the three-dimensional attribute model as R; The local area to be updated is denoted as P', a buffer zone with a distance of R is added to the periphery of the local area to be updated to obtain a new area, denoted as P; The new sample points and the existing sample points located in the area P are extracted as the sample points of the three-dimensional attribute model; S14, the three-dimensional structure model and the three-dimensional attribute model obtained by construction are replaced into the existing three-dimensional geological model.

2. The method of claim 1, wherein, In step S11, specifically comprising: The range affected by the newly added geological exploration data is automatically circled in the form of a polygon on the three-dimensional geological model as the local area to be updated of the three-dimensional geological model.

3. The method of claim 1, wherein, In step S12, the geological exploration data of the local area to be updated is extracted, specifically comprising: All geological exploration data in the local area to be updated, including new data and historical data, is extracted; Virtual geological exploration data is extracted along the boundary of the local area to be updated.

4. The method of claim 1, wherein, In step S12, the three-dimensional structure model of the local area is constructed, specifically comprising: The boundary of the polygon in the local area to be updated is taken as the boundary constraint of the three-dimensional structure model, and the extracted geological exploration data is taken as the input data of the three-dimensional structure model to construct the three-dimensional structure model, wherein the three-dimensional structure model is composed of a plurality of geological bodies, and each geological body records the code of the stratum to which it belongs.

5. The method of claim 4, wherein, In step S14, the three-dimensional structure model is replaced into the existing three-dimensional geological model, specifically comprising: The geological bodies expressed in vector data in the three-dimensional structure model are decomposed into vector grids, and the vector grids are replaced into the corresponding grids in the three-dimensional geological model one by one.

6. The method of claim 1, wherein, In step S13, the grid model is calculated by interpolation, specifically comprising: Through the obtained Kriging parameters and the extracted sample points, the strata to which the grid model belongs are traversed, and Kriging interpolation is performed on the grid model respectively to complete the construction of the three-dimensional attribute model.

7. The method of claim 1, wherein, In step S14, the three-dimensional attribute model of the local area to be updated is replaced into the existing three-dimensional geological model, specifically comprising: Directly replacing the grid in the three-dimensional attribute model to the grid in the corresponding position of the existing three-dimensional geological model.

Images