Three-dimensional geological attribute model dynamic visualization method and system based on element universe, medium and equipment

Abstract

The invention relates to the field of three-dimensional geological modeling and element universe, and discloses a three-dimensional geological attribute model dynamic visualization method, system, medium and equipment based on element universe, and the method comprises the steps: reading top and bottom surface coordinates, depth coordinates, an ACTNUM table and an attribute table in a loaded three-dimensional geological attribute model, calculating eight vertex coordinates, an activation state and an attribute table of each grid; judging whether permeability layer-by-layer visualization is carried out or not according to a calculation result; and importing the processed model file into Unity3D, setting a dynamic visualization function of permeability or saturation, exporting a dynamic scene as an APK from the Unity3D, installing the APK to MR equipment, and clicking two dynamic visualization buttons on the MR equipment to dynamically characterize attributes. According to the method, the three-dimensional geological attribute model is dynamically displayed in the meta-universe environment, the visualization effect and understanding of the underground geological structure are improved, and a more scientific and efficient solution is provided for exploration and development of oil and gas fields.

Description

Technical Field

[0001] The present invention relates to the technical fields of three-dimensional geological modeling and the metaverse, and particularly to a method, system, medium and device for dynamically visualizing a three-dimensional geological attribute model based on the metaverse. Background Art

[0002] With the continuous development of the global economy, the demand for energy continues to grow. As the main energy resources, the exploration, development and production of oil and gas have become increasingly prominent. However, the geological structure of oil and gas fields is complex and changeable, and the distribution of reservoirs has a high degree of uncertainty. In order to more effectively explore and develop oil and gas resources, establishing an accurate geological attribute model has become the key. Through the geological attribute model, important parameters such as the underground geological structure, rock properties, porosity, and permeability can be quantitatively described and analyzed, providing a scientific basis for the evaluation of oil and gas reservoirs, the formulation of development plans, and the optimization of production processes.

[0003] In recent years, the metaverse has rapidly emerged as a new technology. The metaverse is a digital world that combines virtual and real worlds. It provides users with an immersive interactive experience through various technical means such as virtual reality, augmented reality, and artificial intelligence. In the energy field, the metaverse technology has great application potential. Through the metaverse platform, complex geological data can be transformed into intuitive three-dimensional visualization models, enabling geologists and engineers to more intuitively observe and analyze the underground geological structure. At the same time, the interactive characteristics of the metaverse can enable multi-person collaborative work, improving work efficiency and the accuracy of decision-making. In addition, the metaverse can also combine big data analysis and artificial intelligence algorithms to deeply mine and predict geological data, providing a more intelligent solution for the exploration and development of oil and gas fields.

[0004] Grdecl is a file format used to store input data for geological models and reservoir simulations, and is widely used in the oil and gas field. The structure of the grdecl file usually consists of multiple keywords, each keyword representing a data block. These keywords define different data contents. For example, DIMENS defines the dimensions of the grid (NX, NY, NZ), indicating the grid scale of the model; COORD and ZCORN are used to define the coordinates of the grid, describing the geometry of the model. COORD defines the coordinates of the four corner points of the top and bottom surfaces of the columnar grid, and ZCORN provides the corner point coordinates of each grid to support irregular grids; ACTNUM is the activation array, used to indicate which grids in the grid are active (value is 1) or inactive (value is 0). Inactive grids will not participate in the simulation operation.

[0005] Although the grdecl format is suitable for reservoir simulation and geological modeling, it has significant limitations on real-time 3D visualization platforms such as Unity3D. First, there is an incompatibility between the grdecl file format and the standard model file formats of real-time 3D engines such as Unity3D (such as OBJ, FBX, GLTF). Directly reading the grdecl file and using it in Unity3D requires additional parsing and conversion steps. Second, the grdecl file contains multiple attributes, but Unity3D does not directly support this block-based attribute management. To implement the loading, mapping, and updating of these attributes, additional parsing and conversion code needs to be written. Third, the grdecl file is very large, while the real-time rendering and computing capabilities of Unity3D are limited. Especially on mobile devices or head-mounted display devices (such as MR headsets), the large amount of data will cause slow loading, insufficient memory, and rendering lag. Fourth, the application of Unity3D in the metaverse often requires real-time interaction, such as dynamically changing object attributes and dynamically loading resources. However, the grdecl data file requires complex parsing and preprocessing and is difficult to achieve fast real-time updates. Fifth, the file structure of grdecl does not support directly operating on the attributes of objects, such as the object hierarchy (GameObject, Transform, etc.) and script control methods in Unity3D. This kind of interaction support is very important in the metaverse, but grdecl is not suitable for directly mapping to the object system of Unity3D.

[0006] Therefore, how to achieve the dynamic visualization of the 3D geological attribute model based on the metaverse is of great significance. Summary of the Invention

[0007] Aiming at the above problems, the purpose of the present invention is to provide a method, system, medium, and device for dynamic visualization of a 3D geological attribute model based on the metaverse, which solves the incompatibility problem between the grdecl file format and real-time 3D engines such as Unity3D, so as to observe the real-time update and dynamic display of the 3D geological attribute model in the metaverse scenario.

[0008] To achieve the above purpose, in the first aspect, the technical solution adopted by the present invention is: A method for dynamic visualization of a 3D geological attribute model based on the metaverse, which includes:

[0009] Read the top and bottom surface coordinates, depth coordinates, ACTNUM table, and attribute table in the loaded 3D geological attribute model;

[0010] According to the loaded information, calculate the eight vertex coordinates, activation status, and attribute table of each grid;

[0011] Select the visualization method according to different visualization requirements. If you want to view the inside of the model, select layer-by-layer visualization; otherwise, perform the visualization of the hollowed-out model:

[0012] If layer-by-layer visualization is selected, set the number of layers to be skipped for cross-layer extraction, obtain the indices of all the layers to be extracted, and obtain the grid coordinates and attribute information layer by layer according to the layer indices to generate a new OBJ model file and MTL material file; repeat this step until all the specified layers are completely processed;

[0013] If layer-by-layer visualization is not selected, first perform internal hollowing calculation to obtain a new ACTNUM table, and obtain all the grid indices to be extracted. Generate a new OBJ model file according to the grid indices and coordinates; then determine whether all time steps are traversed. If not, generate a new MTL material file according to the attributes of the current time step until all time steps are traversed;

[0014] Import the processed model file into Unity3D, set the dynamic visualization function of permeability or saturation, export the dynamic scene from Unity3D as an APK and install it on the MR device, and click the two dynamic visualization buttons on the MR device to perform dynamic property characterization.

[0015] Furthermore, the 3D geological attribute model is in the grdecl file format; the grdecl file stores the data of the 3D geological model in text format, including the grid structure and physical property parameters;

[0016] In the grdecl file, the 3D geological model is divided into regular or irregular 3D grids, called grids or cells; each grid is defined by eight vertices, and its coordinates are used to describe the spatial position of the grid to form the top and bottom surfaces;

[0017] The grdecl file includes the following data:

[0018] COORD stores the coordinate information of each vertex, which is used to define the geometric structure of the grids in the model;

[0019] ZCORN includes the Z coordinates of the vertices of each grid, which determines the thickness and shape of each cell;

[0020] ACTNUM is used to identify the validity of the grid. A value of 1 indicates that the grid is valid in the model, and a value of 0 indicates that the grid is invalid or a blank area;

[0021] Attribute data, which is used for reservoir simulation and geological analysis.

[0022] Furthermore, the attribute table of each grid includes permeability and saturation;

[0023] Permeability refers to the ease with which a fluid flows through a rock; porosity refers to the ratio of the pore volume to the total volume in a rock.

[0024] Furthermore, set the number of layers to be skipped for cross-layer extraction, obtain the indices of all the layers to be extracted, and layer by layer according to the layer indices, obtain the grid coordinates and attribute information, and generate a new OBJ model file, including:

[0025] Map the coordinates of the corner grid in the grdecl file format to the OBJ format. The data of the corner grid is stored layer by layer. During processing, convert it into a list and assign an index to each grid.

[0026] According to the index and the size of the three-dimensional geological attribute model, extract specific layers to obtain the eight vertex coordinates of each grid.

[0027] Furthermore, set the number of layers to be skipped for cross-layer extraction, obtain the indices of all the layers to be extracted, and layer by layer according to the layer indices, obtain the grid coordinates and attribute information, and generate an MTL material file, including:

[0028] Map the attribute data in the corner grid to colors. The attribute data is converted into reflection color and reflectance material data through the mapping process; the reflection color represents the color displayed by the attribute during rendering, while the reflectance represents the ability of the area to reflect light.

[0029] During the process of generating the material, the mapping operation makes different areas present different colors, forming a visual display effect of the attributes.

[0030] Furthermore, perform internal hollowing calculation to obtain a new ACTNUM table, including:

[0031] Take the three-dimensional geological attribute model V represented by three-dimensional grids as input data, and an actnum array containing grid indices;

[0032] Assume that the model V is a three-dimensional body, and each point in the model space is represented by three-dimensional coordinates (x, y, z). For each point, there is a corresponding activation value actnum(x, y, z);

[0033] Define the boundary condition of the model as the set of points B that satisfy the set conditions, and on these points, the activation value actnum(x, y, z) is equal to 1;

[0034] For each non-boundary point (x, y, z), define its neighborhood N(x, y, z). If each non-boundary point (x, y, z) is completely surrounded by the points in its neighborhood and is not connected to the boundary B, then update the activation value actnum(x, y, z) to 0;

[0035] Perform hollowing determination, set the set I of all points inside that are not connected to the boundary B, and update the activation value actnum(x, y, z):

[0036]

[0037] After hollowing determination, use connected components to verify whether all the outer shells are connected as a whole to ensure that the outer shell grids form a connected structure;

[0038] Output the updated actnum array, where the internal grids are 0 and the outer shell grids are 1.

[0039] Furthermore, set the condition to be one of the following conditions:

[0040] B = {(x, y, z) ∈ V | x = 0 or x = X max or y = 0 or y = Y max or z = 0 or z = Z max}

[0041] In the formula, X max represents the maximum x-direction index of the 3D model; Y max represents the maximum y-direction index of the 3D model; Z max represents the maximum depth z-direction index of the 3D model.

[0042] In a second aspect, the technical solution adopted by the present invention is: a dynamic visualization system for a 3D geological attribute model based on the metaverse, which includes:

[0043] A model reading module that reads the top and bottom surface coordinates, depth coordinates, ACTNUM table, and attribute table in the loaded 3D geological attribute model;

[0044] A data preparation module that calculates the eight vertex coordinates, activation status, and attribute table of each grid according to the loaded information;

[0045] A visualization processing selection module that selects a visualization method according to different visualization requirements. If you want to view the inside of the model, select layer-by-layer visualization; otherwise, perform visualization of the hollowed model:

[0046] If layer-by-layer visualization is selected, set the number of layers to be skipped for cross-layer extraction, obtain the indices of all the layers to be extracted, and obtain the grid coordinates and attribute information layer by layer according to the layer indices to generate a new OBJ model file and MTL material file; repeat this step until all the specified layers are completely processed;

[0047] If layer-by-layer visualization is not selected, first perform internal hollowing calculation to obtain a new ACTNUM table, and obtain all the grid indices to be extracted. Generate a new OBJ model file based on the grid indices and coordinates. Then, determine whether to traverse all time steps. If not completed, generate a new MTL material file according to the attributes of the current time step until all time steps are traversed.

[0048] The visualization module imports the processed model file into Unity3D, sets the dynamic visualization function of permeability or saturation, exports the dynamic scene from Unity3D as an APK and installs it on the MR device, and clicks the two dynamic visualization buttons on the MR device for dynamic property characterization.

[0049] In a third aspect, the technical solution adopted by the present invention is: a computer-readable storage medium storing one or more programs, the one or more programs including instructions which, when executed by a computing device, cause the computing device to execute any of the above methods.

[0050] In a fourth aspect, the technical solution adopted by the present invention is: a computing device, which includes: one or more processors, a memory, and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, and the one or more programs include instructions for executing any of the above methods.

[0051] Due to the above technical solutions adopted by the present invention, it has the following advantages:

[0052] 1. The dynamic visualization method of the 3D geological property model based on the metaverse of the present invention successfully solves the incompatibility problem between the grdecl file format and real-time 3D engines such as Unity3D. By converting the corner point grid of the grdecl file format into an OBJ model file and an MTL material file, the 3D geological property model can run on the Unity3D platform and is finally imported into the MR device for viewing and analysis. This greatly improves the visualization effect of the geological property model. Users can observe the 3D geological property model in a more intuitive and realistic way in the metaverse environment, so as to better understand the underground geological structure and provide a more accurate basis for the exploration and development of oil and gas fields.

[0053] 2. The dynamic visualization method of the 3D geological attribute model based on the metaverse of the present invention, on the one hand, realizes the interactive dynamic display of the 3D geological attribute model in the metaverse environment. Users can interact through MR devices to achieve the dynamic display of attributes, enhancing the interactivity between users and the model and improving work efficiency and decision-making accuracy. On the other hand, for the case where layer-by-layer visualization is not selected, by performing internal hollowing calculations, the storage and rendering pressure is reduced, and the visualization performance is optimized. At the same time, the list structure organization of OBJ files and the method of mapping the attribute values of MTL files to colors also provide a clearer and more understandable model display for users, further enhancing the interactive experience.

[0054] 3. The present invention can be popularized and applied in various software such as 3D geological information systems, geographic information systems, geological modeling and simulation systems, etc. BRIEF DESCRIPTION OF THE DRAWINGS

[0055] Figure 1 is the flowchart of the dynamic visualization method of the 3D geological attribute model based on the metaverse in the embodiment of the present invention;

[0056] Figure 2 is the schematic diagram of the structure data mapping in the embodiment of the present invention;

[0057] Figure 3 is the schematic diagram of the attribute data mapping in the embodiment of the present invention;

[0058] Figure 4 is the experimental case of the layer-by-layer dynamic visualization of the permeability attribute model in the embodiment of the present invention;

[0059] Figure 5 is the experimental case of the dynamic visualization of the saturation attribute model in the embodiment of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0060] In view of the urgent demand for the dynamic visualization of 3D geological attribute models in the metaverse scenario in the current energy field and the deficiencies of existing grdecl format models, it is of great significance to propose a dynamic visualization method of 3D geological attribute models based on the metaverse. This method combines metaverse technology with geological modeling and utilizes the powerful visualization and interaction functions of the metaverse to achieve the dynamic display and real-time analysis of geological attribute models. Through this method, complex geological data can be transformed into a realistic 3D scene, and users can freely view and analyze geological models in the metaverse environment, improving the understanding and recognition of underground geological structures, thereby promoting the efficient development and management of oil and gas resources and providing a more scientific and efficient solution for the exploration and development of oil and gas fields.

[0061] To solve the problems that the grdecl file format is incompatible with real-time 3D engines such as Unity3D, making it difficult to apply to the metaverse scenario, and that the 3D geological attribute model is difficult to update and display dynamically in real time, the present invention provides a method, system, medium and device for dynamic visualization of 3D geological attribute models based on the metaverse. The corner point grid in the grdecl file format is converted into an OBJ model file and an MTL material file, and then imported into the MR headset through the Unity3D platform to view and analyze the dynamic 3D geological attribute model, aiming to improve the visualization effect of the geological attribute model, enhance the interactivity and improve the calculation efficiency, and provide a more scientific and efficient solution for the exploration and development of oil and gas fields. The present invention loads the 3D geological attribute model, extracts the top and bottom surface coordinates, depth coordinates, ACTNUM table and attribute table to obtain the basic data, and calculates the eight vertex coordinates, activation status and attribute information of each grid based on this as the preparatory work for visualization. According to the requirements, the user can choose whether to perform layer-by-layer visualization of permeability. If layer-by-layer visualization is selected, the cross-layer extraction interval is set, the index of the required layer is obtained, and the grid coordinates and attribute information of each layer are extracted according to the index to generate new OBJ and MTL files until all specified layers are processed. During the generation of the OBJ file, the corner point grid coordinates are processed layer by layer into a list structure, and an index is assigned to each grid. For example, the 0th grid is stored as list[0], and the 10,000th grid is stored as list[10,000], thus organizing into an OBJ model file. The generation of the MTL file maps the corner point grid attribute values in the grdecl file format to colors, which is achieved through an RGB mapping table. For example, an attribute value of 1 is mapped to yellow, and an attribute value of 100 is mapped to blue, and each color corresponds to a specific attribute value. If layer-by-layer visualization is not selected, a new ACTNUM table is obtained through hollowing calculation, all grid indexes are extracted, and a new OBJ file is generated. The hollowing process is to reduce the storage and rendering pressure and optimize the visualization performance. The internal grids are hollowed out through the hollowing algorithm to form a shell. Thereafter, an attribute MTL file is generated according to the current time step, and after traversing all time steps, it is imported into Unity3D to achieve dynamic visualization of permeability or saturation. Finally, the dynamic scene is exported as an APK and deployed to the MR device, and the dynamic display of attributes is realized through device interaction. The present invention realizes the dynamic display of 3D geological attribute models in the metaverse environment, improves the visualization effect and the understanding of underground geological structures, and provides a more scientific and efficient solution for the exploration and development of oil and gas fields.

[0062] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Apparently, the described embodiments are some, but not all, of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the described embodiments of the present invention fall within the scope of protection of the present invention.

[0063] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present invention. As used herein, unless the context clearly indicates otherwise, the singular form is also intended to include the plural form. In addition, it should be understood that when the terms "comprise" and / or "include" are used in this specification, they indicate the presence of features, steps, operations, devices, components, and / or combinations thereof.

[0064] In one embodiment of the present invention, a method for dynamically visualizing a three-dimensional geological attribute model based on the metaverse is provided. In this embodiment, as Figure 1 shown, the method includes the following steps:

[0065] 1) Read the top and bottom surface coordinates, depth coordinates, ACTNUM table, and attribute table in the loaded three-dimensional geological attribute model to provide basic data for subsequent processing.

[0066] 2) Calculate the eight vertex coordinates, activation status, and attribute table of each grid according to the loaded information to prepare data for visualization.

[0067] 3) Select a visualization method according to different visualization requirements. If you want to view the inside of the model, select layer-by-layer visualization; otherwise, perform visualization of the hollowed-out model.

[0068] 4) If layer-by-layer visualization is selected in step 3), set the number of layers to be skipped for cross-layer extraction, obtain the indices of all layers to be extracted, and obtain the grid coordinates and attribute information layer by layer according to the layer indices to generate a new OBJ model file and MTL material file; repeat this step until all specified layers are completely processed.

[0069] 5) If layer-by-layer visualization is not selected in step 3), first perform internal hollowing calculation to obtain a new ACTNUM table, and obtain all the grid indices to be extracted. Generate a new OBJ model file according to the grid indices and coordinates. Then, determine whether all time steps have been traversed. If not, generate a new MTL material file according to the attributes of the current time step until all time steps have been traversed.

[0070] 6) Import the processed model file into Unity3D, set the dynamic visualization function of permeability or saturation, export the dynamic scene from Unity3D as an APK and install it on the MR device, and click the two dynamic visualization buttons on the MR device to perform dynamic characterization of properties.

[0071] Thus, the interactive dynamic display of the 3D geological property model in the metaverse environment is realized.

[0072] In the above step 1), the loaded 3D geological property model is in the grdecl file format; the grdecl file stores the data of the 3D geological model in text format, including the grid structure and physical property parameters.

[0073] In the grdecl file, the 3D geological model is divided into regular or irregular 3D grids, called grids or cells; each grid is defined by eight vertices, and its coordinates are used to describe the spatial position of the grid to form the top and bottom surfaces.

[0074] The grdecl file includes the following types of data:

[0075] COORD stores the coordinate information of each vertex, which is used to define the geometric structure of the grids in the model;

[0076] ZCORN includes the Z coordinates of the vertices of each grid, which determines the thickness and shape of each cell;

[0077] ACTNUM is used to identify the validity of the grid, a value of 1 indicates that the grid is valid in the model, and a value of 0 indicates that the grid is invalid or a blank area;

[0078] Property data, such as permeability, porosity, and saturation, etc., are used for reservoir simulation and geological analysis.

[0079] In the above step 2), the property table of each grid includes permeability and saturation; permeability refers to the ease with which a fluid (such as water, oil, or natural gas) flows through a rock; it is a measure of the ability of a rock to allow a fluid to pass through under a pressure difference, usually measured in millidarcies (mD).

[0080] Porosity is the ratio of the pore volume to the total volume in a rock, usually expressed as a percentage. It describes the amount of voids inside the rock and is an important indicator for measuring the fluid storage capacity of a rock.

[0081] In the above step 4), set the number of layers to be spaced for cross-layer extraction, obtain the indices of all the layers to be extracted, and obtain the grid coordinates and property information layer by layer according to the layer indices to generate a new OBJ model file, as Figure 2 shown, including the following steps:

[0082] 4.1) Map the coordinates of the corner grid in the grdecl file format to the OBJ format. The data of the corner grid is stored layer by layer. During processing, convert it into a list and assign an index to each grid;

[0083] For example, the 0th grid corresponds to list[0], and the 10000th grid corresponds to list

[10000] .

[0084] 4.2) Extract specific layers according to the index and the size of the three-dimensional geological attribute model. Through such index processing, the eight vertex coordinates of each grid can be obtained.

[0085] In step 4) above, set the number of layers to be skipped for cross-layer extraction, obtain the indexes of all layers to be extracted, obtain the grid coordinates and attribute information layer by layer according to the layer indexes, and generate an MTL material file. As Figure 3 shown, it includes the following steps:

[0086] 4.3) Map the attribute data in the corner grid to colors. The attribute data is converted into material data such as reflection color and reflectivity through the mapping process. For example, the attribute value can correspond to an RGB color table, and the attribute data is assigned a specific color during this mapping process to visually display the characteristics of different attribute regions; The reflection color represents the color displayed by the attribute during rendering, while the reflectivity represents the light reflection ability of this region.

[0087] 4.4) During the process of generating materials, the mapping operation makes different regions present different colors, forming a visual attribute display effect.

[0088] In step 5) above, the hollowing operation is required to optimize the storage and visualization performance. Traditional MR devices are difficult to support offline loading and visualization of large-scale data, which easily leads to problems such as overheating or freezing of the device. Performing the hollowing operation can reduce the pressure of storage and rendering.

[0089] Specifically, perform internal hollowing calculation to obtain a new ACTNUM table, including the following steps:

[0090] 5.1) Take the three-dimensional geological attribute model V represented by three-dimensional grids as input data, including the actnum array containing grid indexes;

[0091] 5.2) Initialization: Set the model V as a three-dimensional body. Each point in the model space is represented by three-dimensional coordinates (x, y, z), and for each point, there is a corresponding activation value actnum(x, y, z);

[0092] At the beginning:

[0093]

[0094] 5.3) Marking the boundary: Define the boundary condition of the model as the point set B that satisfies the set conditions, and activate the numerical value actnum(x, y, z) equal to 1 at these points, i.e.: actnum(x, y, z) = 1,

[0095] where the set condition is one of the following conditions:

[0096] B = {(x, y, z) ∈ V | x = 0 or x = X max or y = 0 or y = Y max or z = 0 or z = Z max}

[0097] In the formula, X max represents the maximum x-direction index value of the three-dimensional model; Y max represents the maximum y-direction index value of the three-dimensional model; Z max represents the maximum depth z-direction index value of the three-dimensional model.

[0098] 5.4) Judgment of internal grids: For each non-boundary point (x, y, z), define its neighborhood N(x, y, z). If each non-boundary point (x, y, z) is completely surrounded by the points in its neighborhood and is not connected to the boundary B, then update the activation numerical value actnum(x, y, z) to 0;

[0099] In this embodiment, the neighborhood N(x, y, z) is:

[0100] N(x, y, z) = {(x ′ , y ′ , z ′ ) || |x ′ - x| ≤ 1, |y ′ - y| ≤ 1, |z ′ - z| ≤ 1, (x ′ , y ′ , z ′ ) ≠ (x, y, z)}

[0101] 5.5) Perform hollowing determination, set the set I of all points inside that are not connected to the boundary B, and update the activation numerical value actnum(x, y, z):

[0102]

[0103] Among them, the set I is defined as:

[0104] I = {(x, y, z) ∈ V | if (x, y, z) is not connected to B and it is surrounded by N(x, y, z)}

[0105] 5.6) Connectivity verification: After the hollowing determination, the connected components are used to verify whether all the outer shells are connected as a whole to ensure that the lattice S of the outer shell forms a connected structure;

[0106] S = {(x, y, z) ∈ V | actnum(x, y, z) = 1 and connected to B}

[0107] 5.7) Output the updated actnum array, where the internal lattice is 0 and the outer shell lattice is 1.

[0108] In step 6) above, the processed model file is imported into Unity3D, the dynamic visualization function of permeability or saturation is set, the dynamic scene is exported from Unity3D as an APK and installed on the MR device, and the two dynamic visualization buttons are clicked on the MR device for dynamic property characterization. As Figure 4 shown, it is a case of the layer-by-layer dynamic visualization experiment of the permeability property model. As Figure 5 shown, it is a case of the dynamic visualization experiment of the saturation property model. These two experimental cases show that the dynamic visualization method of the three-dimensional geological property model based on the metaverse provided by the present invention can more intuitively display the key properties in the geological model and provide strong support for geological analysis and research.

[0109] In an embodiment of the present invention, a dynamic visualization system for a three-dimensional geological property model based on the metaverse is provided, which includes:

[0110] A model reading module that reads the top and bottom surface coordinates, depth coordinates, ACTNUM table, and property table in the loaded three-dimensional geological property model;

[0111] A data preparation module that calculates the eight vertex coordinates, activation status, and property table of each lattice according to the loaded information;

[0112] A visualization processing selection module that selects a visualization method according to different visualization requirements. If you want to view the inside of the model, then select layer-by-layer visualization, otherwise perform the visualization of the hollowed-out model:

[0113] If layer-by-layer visualization is selected, set the number of layers to be skipped for cross-layer extraction, obtain the indices of all the layers to be extracted, obtain the lattice coordinates and property information layer by layer according to the layer indices, and generate new OBJ model files and MTL material files; repeat this step until all the specified layers are completely processed;

[0114] If layer-by-layer visualization is not selected, first perform internal hollowing calculation to obtain a new ACTNUM table, and obtain all the grid indices to be extracted. Generate a new OBJ model file based on the grid indices and coordinates. Then, determine whether all time steps are traversed. If not, generate a new MTL material file according to the attributes of the current time step until all time steps are traversed.

[0115] The visualization module imports the processed model file into Unity3D, sets the dynamic visualization function of permeability or saturation, exports the dynamic scene from Unity3D as an APK and installs it on the MR device, and clicks the two dynamic visualization buttons on the MR device for dynamic property characterization.

[0116] In the above embodiment, the 3D geological property model is in the grdecl file format; the grdecl file stores the data of the 3D geological model in text format, including the grid structure and physical property parameters.

[0117] In the grdecl file, the 3D geological model is divided into regular or irregular 3D grids, called grids or cells; each grid is defined by eight vertices, and its coordinates are used to describe the spatial position of the grid to form the top and bottom surfaces.

[0118] The grdecl file includes the following data:

[0119] COORD stores the coordinate information of each vertex, which is used to define the geometric structure of the grids in the model.

[0120] ZCORN includes the Z coordinates of the vertices of each grid, which determines the thickness and shape of each cell.

[0121] ACTNUM is used to identify the validity of the grid. A value of 1 indicates that the grid is valid in the model, and a value of 0 indicates that the grid is invalid or a blank area.

[0122] Property data, which is used for reservoir simulation and geological analysis.

[0123] In the above embodiment, the property table of each grid includes permeability and saturation; permeability refers to the ease of fluid flow in the rock; porosity refers to the ratio of the pore volume to the total volume in the rock.

[0124] In the above embodiment, set the number of layers to be skipped for cross-layer extraction, obtain the indices of all the layers to be extracted, and obtain the grid coordinates and property information layer by layer according to the layer indices to generate a new OBJ model file, including:

[0125] Map the coordinates of the corner point grids in the grdecl file format to the OBJ format. The data of the corner point grids is stored layer by layer. During processing, convert it into a list and assign an index to each grid.

[0126] Extract specific layers according to the index and the size of the 3D geological attribute model to obtain the eight vertex coordinates of each grid.

[0127] In the above embodiments, set the number of layers to be skipped for cross-layer extraction, obtain the indices of all layers to be extracted, and layer by layer according to the layer indices, obtain the grid coordinates and attribute information, and generate an MTL material file, including:

[0128] Map the attribute data in the corner point grid to colors. The attribute data is converted into reflection color and reflectivity material data through the mapping process; the reflection color represents the color displayed by the attribute during rendering, while the reflectivity represents the light reflection ability of the area;

[0129] During the process of generating the material, the mapping operation makes different areas present different colors, forming a visual attribute display effect.

[0130] In the above embodiments, perform internal hollowing calculation to obtain a new ACTNUM table, including:

[0131] Use the 3D geological attribute model V represented by 3D grids as input data, an actnum array containing grid indices;

[0132] Assume that the model V is a 3D body, and each point in the model space is represented by 3D coordinates (x, y, z). For each point, there is a corresponding activation value actnum(x, y, z);

[0133] Define the boundary condition of the model as the set of points B that satisfy the set conditions, and on these points, the activation value actnum(x, y, z) is equal to 1;

[0134] For each non-boundary point (x, y, z), define its neighborhood N(x, y, z). If each non-boundary point (x, y, z) is completely surrounded by the points in its neighborhood and is not connected to the boundary B, then update the activation value actnum(x, y, z) to 0;

[0135] Perform hollowing determination, set the set of points I inside that are not connected to the boundary B, and update the activation value actnum(x, y, z):

[0136]

[0137] After the hollowing determination, use the connected components to verify whether all the outer shells are connected as a whole to ensure that the outer shell grids form a connected structure;

[0138] Output the updated actnum array, where the internal grids are 0 and the outer shell grids are 1.

[0139] In this embodiment, the set condition is one of the following conditions:

[0140] B = {(x, y, z) ∈ V | x = 0 or x = X max or y = 0 or y = Y max or z = 0 or z = Z max}

[0141] In the formula, X max represents the maximum value of the x-direction index of the three-dimensional model; Y max represents the maximum value of the y-direction index of the three-dimensional model; Z max represents the maximum value of the depth z-direction index of the three-dimensional model.

[0142] The system provided in this embodiment is used to execute the above method embodiments. For the specific process and detailed content, please refer to the above embodiments and will not be elaborated here.

[0143] In an embodiment of the present invention, a computing device is provided. The computing device may be a terminal, which may include: a processor, a communications interface, a memory, a display screen, and an input device. Among them, the processor, the communications interface, and the memory complete communication with each other through a communication bus. The processor is used to provide computing and control capabilities. The memory includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. When the computer program is executed by the processor, the methods in the above embodiments are implemented; the internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The communications interface is used to communicate with an external terminal in a wired or wireless manner. The wireless manner can be implemented through WIFI, a carrier network, NFC (Near Field Communication), or other technologies. The display screen may be a liquid crystal display screen or an electronic ink display screen. The input device may be a touch layer covering the display screen, or a button, a trackball, or a touchpad provided on the housing of the computing device, or an external keyboard, touchpad, or mouse, etc. The processor can call the logical instructions in the memory.

[0144] In addition, when the logical instructions in the above-mentioned memory can be implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memories (ROMs), random access memories (RAMs), magnetic disks, or optical discs that can store program codes.

[0145] In an embodiment of the present invention, a computer program product is provided. The computer program product includes a computer program stored on a non-transitory computer-readable storage medium. The computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the methods provided in the above-mentioned method embodiments.

[0146] In an embodiment of the present invention, a non-transitory computer-readable storage medium is provided. The non-transitory computer-readable storage medium stores server instructions, and the computer instructions cause the computer to execute the methods provided in the above-mentioned embodiments.

[0147] For the computer-readable storage medium provided in the above-mentioned embodiment, its implementation principle and technical effects are similar to those of the above-mentioned method embodiment, and will not be elaborated here.

[0148] The present invention is described with reference to the flowcharts and / or block diagrams of methods, apparatuses (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, and the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices generate means for implementing the functions specified in one process Figure 1 one process or multiple processes and / or blocks Figure 1 one block or multiple blocks.

[0149] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to work in a particular manner, such that the instructions stored in the computer-readable memory produce a manufacture including an instruction means that implements the functions specified in one or more processes and / or blocks Figure 1 in the process or processes and / or blocks Figure 1 specified in the block or blocks.

[0150] These computer program instructions can also be loaded onto a computer or other programmable data processing apparatus, such that a series of operational steps are performed on the computer or other programmable apparatus to produce a computer-implemented process, whereby the instructions executed on the computer or other programmable apparatus provide steps for implementing the functions specified in one or more processes and / or blocks Figure 1 in the process or processes and / or blocks Figure 1 specified in the block or blocks.

[0151] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements for some of the technical features; and these modifications or replacements 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 dynamic visualization method of a three-dimensional geological attribute model based on a metaverse, characterized in that: include: Read the top and bottom surface coordinates, depth coordinates, ACTNUM table and attribute table in the loaded 3D geological attribute model; Based on the loaded information, the eight vertex coordinates, activation status and attribute table of each grid are calculated; Select the visualization method according to different visualization requirements. If you want to view the internal situation of the model, select layer-by-layer visualization. Otherwise, perform visualization of the hollow model: If you choose layer-by-layer visualization, set the number of layers required for cross-layer extraction, obtain the indexes of all layers to be extracted, obtain the grid coordinates and attribute information layer by layer according to the layer index, and generate new OBJ model files and MTL material files; repeat this step until all specified layers are processed; If layer-by-layer visualization is not selected, first perform internal hollowing calculation to obtain a new ACTNUM table, obtain all grid indexes that need to be extracted, and generate a new OBJ model file based on the grid index and coordinates; then determine whether to traverse all time steps. If not completed, generate a new MTL material file based on the attributes of the current time step until all time steps are traversed; Import the processed model file into Unity3D, set the dynamic visualization function of permeability or saturation, export the dynamic scene from Unity3D as APK and install it to the MR device, and click the two dynamic visualization buttons on the MR device to dynamically represent the attributes.

2. The method for dynamic visualization of a three-dimensional geological attribute model based on a metaverse according to claim 1, characterized in that: The three-dimensional geological attribute model is in the grdecl file format; the grdecl file stores the data of the three-dimensional geological model in a text format, including the grid structure and physical property parameters; In the grdecl file, the 3D geological model is divided into regular or irregular 3D grids, called cells or cells; each cell is defined by eight vertices, whose coordinates are used to describe the spatial position of the cell to form the top and bottom surfaces; The grdecl file includes the following data: COORD stores the coordinate information of each vertex and is used to define the geometric structure of the grid in the model; ZCORN includes the Z coordinates of the vertices of each grid, which determines the thickness and shape of each unit; ACTNUM is used to identify the validity of the grid. A value of 1 indicates that the grid is valid in the model, and a value of 0 indicates that the grid is invalid or a blank area. Attribute data for reservoir simulation and geological analysis.

3. The method for dynamic visualization of a three-dimensional geological attribute model based on a metaverse according to claim 1, characterized in that: The attribute table for each grid includes permeability and saturation; Permeability refers to the ease with which fluid flows in rock; porosity refers to the ratio of the pore volume to the total volume of rock.

4. The method for dynamic visualization of a three-dimensional geological attribute model based on a metaverse according to claim 1, characterized in that: Set the number of layers required for cross-layer extraction, obtain the indexes of all layers to be extracted, obtain the grid coordinates and attribute information layer by layer according to the layer index, and generate a new OBJ model file, including: Map the coordinates of the corner grid in the grdecl file format to the OBJ format. The data of the corner grid is stored layer by layer. When processing, it is converted into a list and an index is assigned to each grid. According to the index and the size of the 3D geological attribute model, a specific layer is extracted to obtain the eight vertex coordinates of each grid.

5. The method for dynamic visualization of a three-dimensional geological attribute model based on a metaverse according to claim 1, characterized in that: Set the number of layers required for cross-layer extraction, obtain the indexes of all layers that need to be extracted, obtain the grid coordinates and attribute information layer by layer according to the layer index, and generate the MTL material file, including: Map the attribute data in the corner grid to color. The attribute data is converted into reflection color and reflectivity material data through the mapping process. The reflection color indicates the color of the attribute when rendering, while the reflectivity indicates the ability of the area to reflect light. In the process of generating materials, the mapping operation makes different areas appear in different colors, forming a visual attribute display effect.

6. The method for dynamic visualization of a three-dimensional geological attribute model based on a metaverse according to claim 1, characterized in that: Perform internal hollowing calculations to obtain a new ACTNUM table, including: The three-dimensional geological attribute model V represented by a three-dimensional grid is used as input data, including the actnum array of grid index; Assume that the model V is a three-dimensional body, each point in the model space is represented by a three-dimensional coordinate (x, y, z), and each point has a corresponding activation value actnum(x, y, z); Define the boundary conditions of the model as a set of points B that meet the set conditions, and activate the value actnum(x,y,z) at these points equal to 1; For each non-boundary point (x, y, z), define its neighborhood N(x, y, z). If each non-boundary point (x, y, z) is completely surrounded by points in its neighborhood and is not connected to the boundary B, then update the activation value actnum(x, y, z) to 0; Perform hollowing judgment, set all internal points that are not connected to the boundary B as I, and update the activation value actnum(x,y,z): After the hollowing determination, the connected components are used to verify whether all the external shells are connected as a whole, ensuring that the external shell lattices form a connected structure; Output the updated actnum array, where the inner grid is 0 and the outer shell grid is 1.

7. The method for dynamic visualization of a three-dimensional geological attribute model based on a metaverse according to claim 6, characterized in that: Set the condition to one of the following: B={(x,y,z)∈V|x=0orx=X max Or y=0 or y=Y max Or z=0 or z=Z max } Where, X max Indicates the maximum x-direction index of the 3D model; Y max Indicates the maximum value of the y-direction index of the three-dimensional model; Z max Indicates the maximum value of the depth z-direction index of the 3D model.

8. A three-dimensional geological attribute model dynamic visualization system based on the metaverse, characterized in that: include: Model reading module, reading the top and bottom surface coordinates, depth coordinates, ACTNUM table and attribute table in the loaded three-dimensional geological attribute model; The data preparation module calculates the eight vertex coordinates, activation status and attribute table of each grid based on the loaded information; In the visualization processing selection module, select the visualization method according to different visualization requirements. If you want to view the internal situation of the model, select layer-by-layer visualization, otherwise perform visualization of the hollow model: If you choose layer-by-layer visualization, set the number of layers required for cross-layer extraction, obtain the indexes of all layers to be extracted, obtain the grid coordinates and attribute information layer by layer according to the layer index, and generate new OBJ model files and MTL material files; repeat this step until all specified layers are processed; If layer-by-layer visualization is not selected, first perform internal hollowing calculation to obtain a new ACTNUM table, obtain all grid indexes that need to be extracted, and generate a new OBJ model file based on the grid index and coordinates; then determine whether to traverse all time steps. If not completed, generate a new MTL material file based on the attributes of the current time step until all time steps are traversed; The visualization module imports the processed model file into Unity3D, sets the dynamic visualization function of permeability or saturation, exports the dynamic scene from Unity3D as APK and installs it to the MR device, and clicks two dynamic visualization buttons on the MR device to dynamically represent the attributes.

9. A computer-readable storage medium storing one or more programs, characterized in that: The one or more programs include instructions, which, when executed by a computing device, cause the computing device to perform any one of the methods of claims 1 to 7.

10. A computing device, characterized in that include: One or more processors, a memory, and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, and the one or more programs include instructions for executing any one of the methods described in claims 1 to 7.

Images