Grid generation method and apparatus for oil reservoir model, and device, storage medium and program product

Abstract

A grid generation method and apparatus for an oil reservoir model, and a device, a storage medium and a program product, which relate to the technical field of oil exploration. The method comprises: acquiring an initial oil reservoir grid model, wherein the initial oil reservoir grid model is obtained by performing tetrahedral meshing on an initial oil reservoir model (S201); using a preset layer extraction strategy to extract sub-layers from the initial oil reservoir grid model, so as to generate an initial sub-layer model, wherein the initial sub-layer model is composed of a plurality of sub-layers, and each sub-layer comprises a plurality of triangular mesh elements (S202); using a preset mesh processing strategy to convert the triangular mesh elements in the plurality of sub-layers into rectangular mesh elements, so as to generate a target sub-layer model (S203); and on the basis of a mesh coordinate system and the target sub-layer model, generating a target oil reservoir grid model, wherein grids in the target oil reservoir grid model are all hexahedral grids (S204). The method of the present application improves the grid quality of the oil reservoir model, ensures the accuracy and reliability of subsequent analysis and research based on the oil reservoir grid model, and is applicable to oil reservoir models with non-transmissive faults and overthrust nappe structures, etc.

Description

Methods, apparatus, equipment, storage media, and program products for mesh generation of reservoir models

[0001] This application claims priority to Chinese Patent Application No. 202411853828.1, filed on December 16, 2024, entitled “Mesh Generation Method, Apparatus, Device, Storage Medium and Program Product for Reservoir Models”, the entire contents of which are incorporated herein by reference. Technical Field

[0002] This application relates to the field of petroleum exploration technology, and in particular to a method, apparatus, equipment, storage medium and program product for generating grids for reservoir models. Background Technology

[0003] Three-dimensional reservoir models have a wide range of applications and play an important role in all aspects of oil exploration and development, serving as a bridge and link for integrated exploration and development. Specifically, reservoir mesh modeling technology is one of the core technologies for reservoir modeling and simulation, and different mesh generation methods and mesh types have a significant impact on the accuracy of modeling and simulation.

[0004] Currently, the traditional Pillar-based corner mesh modeling approach is mainly suitable for thin reservoir mesh modeling with relatively simple fault relationships. However, when dealing with complex multi-level X-type or Y-type intersecting faults, as well as reservoir models with non-permeable faults, the generated mesh effect is poor, which affects the accuracy and reliability of subsequent analysis and processing based on the reservoir mesh model. Summary of the Invention

[0005] This application provides a method, apparatus, device, storage medium, and program product for generating meshes for reservoir models, in order to solve the problem that existing technologies produce poor mesh effects when dealing with complex multi-level X-type or Y-type cross-cutting layers and non-permeable faults in reservoir models.

[0006] In a first aspect, this application provides a method for generating a mesh for a reservoir model, including:

[0007] An initial reservoir mesh model is obtained, which is obtained by dividing the initial reservoir model into tetrahedral meshes.

[0008] A preset layer extraction strategy is used to extract sub-layers from the initial reservoir grid model to generate an initial sub-layer model; the initial sub-layer model is composed of multiple sub-layers, and each sub-layer includes multiple triangular mesh elements.

[0009] A preset mesh processing strategy is used to convert triangular mesh cells in multiple small layers into rectangular mesh cells to generate the target small layer model;

[0010] A target reservoir grid model is generated based on the grid coordinate system and the target small-level model; all grid types in the target reservoir grid model are hexahedral grids.

[0011] In one possible design, the step of extracting sub-layers from the initial reservoir grid model using a preset layer extraction strategy to generate an initial sub-layer model includes:

[0012] The sedimentary patterns corresponding to each stratum are determined using the strata represented in the initial reservoir grid model as the dimension;

[0013] Based on the sedimentary patterns of each stratum and the corresponding preset layer extraction interval values, small layers of each stratum are extracted along the height direction of the initial reservoir grid model.

[0014] The initial small-level model is generated based on the extracted small-levels of each stratum.

[0015] In one possible design, the step of using a preset mesh processing strategy to convert triangular mesh cells in multiple small layers into rectangular mesh cells to generate a target small layer model includes:

[0016] A preset mesh transformation strategy is used to convert triangular mesh units in multiple small layers into hybrid mesh units; the hybrid mesh units include rectangular mesh units far from the fault region and non-rectangular mesh units close to the fault region;

[0017] The non-rectangular mesh cells are processed using a preset mesh regularization strategy to generate a target small-level model; the multiple small levels of the target small-level model are all composed of rectangular mesh cells.

[0018] In one possible design, the adoption of a preset mesh transformation strategy to convert triangular mesh cells in multiple small layers into hybrid mesh cells includes:

[0019] The triangular mesh cells of the small layer are divided by a unit rectangular mesh based on the grid coordinate system to obtain a hybrid mesh cell; the unit rectangular mesh of the grid coordinate system is formed by the intersection of multiple horizontal grid lines and multiple horizontal grid lines.

[0020] In one possible design, the step of using a preset mesh processing strategy to perform mesh regularization processing on the non-rectangular mesh cells to generate the target small-scale model includes:

[0021] Determine whether the non-rectangular mesh cell covers the center point of the unit rectangular mesh in the mesh coordinate system;

[0022] In response to the non-rectangular mesh cell covering the center point of the unit rectangular mesh in the mesh coordinate system, the non-rectangular mesh cell is replaced with the unit rectangular mesh to generate the target small-scale model.

[0023] One possible design also includes:

[0024] In response to the non-rectangular mesh cell not covering the center point of the unit rectangular mesh of the mesh coordinate system, the non-rectangular mesh cell is deleted.

[0025] In one possible design, after replacing the non-rectangular mesh cells with the unit rectangular mesh, the method further includes:

[0026] Determine whether the vertex coordinates of adjacent sides of adjacent rectangular grid cells are consistent;

[0027] In response to the inconsistency of vertex coordinates of adjacent sides of adjacent rectangular grid cells, the vertices of adjacent sides of adjacent rectangular grid cells are merged.

[0028] In one possible design, generating the target reservoir mesh model based on the mesh coordinate system and the target small-level model includes:

[0029] Multiple unit mesh columns are generated based on the mesh coordinate system and the target small-level model, and each unit mesh column includes multiple hexahedral meshes arranged along the extension direction of the unit mesh column.

[0030] Multiple unit grid columns are merged to generate a target reservoir grid model.

[0031] In one possible design, generating multiple element mesh columns based on the mesh coordinate system and the target small-level model includes:

[0032] Based on the horizontal grid lines, horizontal grid lines, and height grid lines of the grid coordinate system, multiple initial grid cylinders are generated.

[0033] The target small-level model is segmented based on the initial mesh pillars to generate multiple unit mesh pillars.

[0034] In one possible design, after generating multiple element mesh columns based on the mesh coordinate system and the target small-level model, the process further includes:

[0035] Determine whether a reverse fault exists within the unit grid column;

[0036] In response to the absence of reverse faults within the unit grid column, the target sub-layers within the unit grid column are arranged according to the sedimentary sequence of the strata.

[0037] In response to the presence of a reverse fault within the unit grid column, the target sub-layers within the unit grid column are arranged from high to low according to their height coordinate values.

[0038] One possible design also includes:

[0039] Determine whether there are any intersecting target sub-layers within the unit grid column;

[0040] In response to determining that there are intersecting target sub-layers within the unit grid column, the target sub-layers within the unit grid column are rearranged from high to low according to their height coordinate values.

[0041] In one possible design, after generating multiple element mesh columns based on the mesh coordinate system and the target small-level model, the process further includes:

[0042] Add a three-dimensional coordinate index to the hexahedral mesh within the unit mesh column.

[0043] One possible design also includes:

[0044] Determine whether there are duplicate 3D coordinate indices for the hexahedral mesh within the unit mesh column;

[0045] In response to the presence of index repetition, adjust the coordinate index value of any dimension of the hexahedral mesh with the repetitive index.

[0046] Secondly, this application provides a mesh generation device for a reservoir model, comprising:

[0047] The acquisition module is used to acquire the initial reservoir mesh model, which is obtained by dividing the initial reservoir model into tetrahedral meshes;

[0048] An extraction module is used to extract sub-layers from the initial reservoir grid model using a preset layer extraction strategy to generate an initial sub-layer model; the initial sub-layer model is composed of multiple sub-layers, and each sub-layer includes multiple triangular mesh elements;

[0049] The conversion module is used to convert triangular mesh cells in multiple small layers into rectangular mesh cells using a preset mesh processing strategy to generate a target small layer model.

[0050] The model generation module is used to generate a target reservoir grid model based on the grid coordinate system and the target small-level model; the grid type in the target reservoir grid model is a hexahedral grid.

[0051] Thirdly, embodiments of this application provide an electronic device, including: a processor, and a memory communicatively connected to the processor;

[0052] The memory stores computer-executed instructions;

[0053] The processor executes computer execution instructions stored in the memory to implement the method as described in any of the first aspects.

[0054] Fourthly, embodiments of this application provide a computer-readable storage medium storing computer-executable instructions, which, when executed by a processor, are used to implement the method as described in any of the first aspects.

[0055] Fifthly, embodiments of this application provide a computer program product, including a computer program that, when executed by a processor, implements the method described in any of the first aspects.

[0056] This application provides a method, apparatus, device, storage medium, and program product for generating reservoir meshes. The method involves obtaining an initial reservoir mesh model, which is obtained by tetrahedral meshing an initial reservoir model. A preset layer extraction strategy is used to extract sub-layers from the initial reservoir mesh model to generate an initial sub-layer model. The initial sub-layer model consists of multiple sub-layers, each containing multiple triangular mesh elements. A preset mesh processing strategy is used to convert the triangular mesh elements in the multiple sub-layers into rectangular mesh elements to generate a target sub-layer model. A target reservoir mesh model is generated based on the mesh coordinate system and the target sub-layer model. All mesh types in the target reservoir mesh model are hexahedral meshes. This reservoir model mesh generation method refines the initial reservoir mesh model into multiple sub-layers using a preset layer extraction strategy. Then, a preset mesh processing strategy converts all triangular mesh elements in the sub-layers into rectangular mesh elements. Finally, a target reservoir mesh model consisting entirely of hexahedral meshes is obtained based on the mesh coordinate system and the target sub-layer model. By refining the initial reservoir mesh model into multiple smaller layers, even if non-permeable faults and thrust-neck structures exist within the initial reservoir mesh model, mesh processing can be performed effectively. Converting all triangular mesh elements to rectangular mesh elements avoids mesh distortion, improves mesh quality, and ensures the accuracy and reliability of subsequent analyses based on the reservoir mesh model. The reservoir model mesh generation method provided in this application is applicable to both simple, thin-layer reservoirs and complex reservoirs with numerous non-permeable faults and thrust-neck structures, demonstrating broad applicability. Attached Figure Description

[0057] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with this application and, together with the description, serve to explain the principles of this application.

[0058] Figure 1 is a schematic diagram of a hexahedral corner mesh model based on Pillar constraints in the prior art;

[0059] Figure 2 is a schematic diagram of the fault structure of a reservoir model in the prior art;

[0060] Figure 3 shows an application scenario of the reservoir model mesh generation method provided in the embodiments of this application;

[0061] Figure 4 is a flowchart of the mesh generation method for the reservoir model provided in the embodiments of this application;

[0062] Figure 5 is a schematic diagram of a partial initial reservoir grid model provided in an embodiment of this application;

[0063] Figure 6 is a mesh generation effect diagram of a reservoir model with complex normal fault intersection structure provided in the embodiment of this application;

[0064] Figure 7 shows the mesh generation effect of a reservoir model with multi-level X-type and Y-type complex reverse faults provided in the embodiment of this application;

[0065] Figure 8 is a mesh generation effect diagram of a reservoir model with a complex reverse thrust fault intersecting frame provided in the embodiment of this application;

[0066] Figure 9 is a schematic diagram of the stratigraphic deposition model provided in the embodiments of this application;

[0067] Figure 10 is a schematic diagram of any small layer with hybrid mesh cells provided in an embodiment of this application;

[0068] Figure 11 is a schematic diagram of the small-level triangular mesh unit segmentation provided in an embodiment of this application;

[0069] Figure 12 is a schematic diagram of geometric inconsistencies in small-level rectangular mesh cells provided in an embodiment of this application;

[0070] Figure 13 is a schematic diagram of the initial grid column provided in an embodiment of this application;

[0071] Figure 14 is a schematic diagram of the unit grid column provided in the embodiment of this application;

[0072] Figure 15 is a schematic diagram of the small-level arrangement within the unit grid column provided in the embodiment of this application;

[0073] Figure 16 is a schematic diagram of the adjustment of the small-level intersection of the unit grid columns provided in the embodiment of this application;

[0074] Figure 17 is a schematic diagram of the reindexing of the hexahedral mesh within the cell grid column provided in the embodiment of this application;

[0075] Figure 18 is a schematic diagram of the structure of the mesh generation device for the reservoir model provided in the embodiment of this application;

[0076] Figure 19 is a schematic diagram of the structure of the electronic device provided in the embodiment of this application.

[0077] The accompanying drawings illustrate specific embodiments of this application, which will be described in more detail below. These drawings and descriptions are not intended to limit the scope of the concept in any way, but rather to illustrate the concept of this application to those skilled in the art through reference to particular embodiments. Detailed Implementation

[0078] Exemplary embodiments will now be described in detail, examples of which are illustrated in the accompanying drawings. When the following description relates to the drawings, unless otherwise indicated, the same numbers in different drawings denote the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with this application. Rather, they are merely examples of apparatuses and methods consistent with some aspects of this application as detailed in the appended claims.

[0079] It should be noted that in the embodiments of this application, certain software, components, models and other existing solutions in the industry may be mentioned. These should be regarded as exemplary and are only intended to illustrate the feasibility of implementing the technical solution of this application. However, it does not mean that the applicant has used or necessarily used the solution.

[0080] To clearly understand the technical solution of this application, the solutions of the prior art will be described in detail first.

[0081] Three-dimensional reservoir models and three-dimensional geological models play a crucial role in the field of oil exploration and development, serving as the foundation for development scheme design at each stage, remaining oil potential tapping, and reservoir numerical dynamic analysis. Reservoir mesh modeling technology, as the core of reservoir modeling and simulation, has a significant impact on the accuracy of modeling and simulation due to the selection of mesh generation methods and mesh types.

[0082] In existing technologies, reservoir meshes are mainly divided into two categories: structured meshes and unstructured meshes. Among them, structured meshes, especially the Pillar-constrained hexahedral corner mesh model, are widely used due to their high computational efficiency, convenient storage, and natural suitability for geostatistical simulations based on unique grid index mapping and numerical simulator calculations based on finite difference. Figure 1 illustrates an example of a Pillar-constrained hexahedral corner mesh model. However, while the Pillar-grid-based hexahedral corner mesh modeling method can strictly guarantee the consistency between the mesh model and internal cross-sections, it often leads to severe mesh distortion when dealing with complex fault relationships. For example, in reservoir models with multi-level X-type or Y-type intersecting faults (refer to Figure 2) and non-permeable faults, this often affects the accuracy and reliability of subsequent analysis and processing. Specifically, the Pillar-grid-based hexahedral corner mesh modeling technique mainly relies on the grid column morphology and is suitable for thin reservoirs with relatively simple fault relationships. In reality, however, reservoir geological structures are often complex and varied, with numerous non-permeable faults and thrust-nappe structures. These complex geological structures make Pillar's hexahedral corner mesh modeling method difficult to adapt, resulting in poor mesh quality that cannot meet the needs of high-precision reservoir modeling and simulation.

[0083] Therefore, in addressing the technical problems of the existing technologies, to generate a high-precision reservoir mesh model and ensure the accuracy and reliability of subsequent analysis and processing based on the reservoir mesh model, the following steps are taken: First, an initial reservoir mesh model with preliminary mesh generation is obtained. Then, this initial reservoir mesh model is further processed to obtain a more regular mesh model. Specifically, when refining the initial reservoir mesh model, it is divided into multiple sub-layers, and the meshes on each sub-layer are processed sequentially until all meshes in each sub-layer are converted into rectangular meshes. Finally, the multiple sub-layers are merged to obtain a reservoir mesh model where all meshes are hexahedral meshes.

[0084] Figure 3 illustrates an application scenario of the reservoir model mesh generation method provided in an embodiment of this application. As shown in Figure 3, the application scenario provided in this embodiment includes a user terminal 100, a mesh generation device 101, and a database 102. Specifically, when it is necessary to mesh the reservoir model, the user terminal 100 sends a mesh generation command to the mesh generation device 101. After receiving the mesh generation command, the mesh generation device 101 retrieves an initial reservoir mesh model from the database 102. This initial reservoir mesh model is obtained by tetrahedral meshing the initial reservoir model. A preset layer extraction strategy is used to extract small layers from the initial reservoir mesh model to generate an initial small layer model. A preset mesh processing strategy is used to convert triangular mesh elements in multiple small layers into rectangular mesh elements to generate a target small layer model. A target reservoir mesh model is generated based on the mesh coordinate system and the target small layer model. All mesh types in the target reservoir mesh model are hexahedral meshes. Finally, the mesh generation device 101 sends the target reservoir mesh model to the user terminal 100 to complete the mesh generation task of the reservoir model.

[0085] The technical solution of this application and how the technical solution of this application solves the above-mentioned technical problems are described in detail below with specific embodiments. These specific embodiments can be combined with each other, and the same or similar concepts or processes may not be described again in some embodiments. The embodiments of this application will now be described with reference to the accompanying drawings.

[0086] Figure 4 is a flowchart of a reservoir model mesh generation method provided in an embodiment of this application. As shown in Figure 4, the execution subject of this embodiment is a reservoir model mesh generation device. This reservoir model mesh generation device can be implemented by a computer program, or by a medium storing the relevant computer program, such as a USB flash drive and / or optical disc. Alternatively, it can be implemented by a physical device integrating or installing the relevant computer program, such as a chip or electronic device. The electronic device can be a computer or a server, etc. The reservoir model mesh generation method provided in this embodiment includes the following steps:

[0087] S201. Obtain the initial reservoir mesh model. The initial reservoir mesh model is obtained by dividing the initial reservoir model into tetrahedral meshes.

[0088] Specifically, an initial reservoir mesh model is first obtained for subsequent mesh processing. This initial reservoir mesh model is obtained by tetrahedral meshing of the initial reservoir model. Refer to Figure 5, which is a schematic diagram of part of the initial reservoir mesh model.

[0089] The initial reservoir model is a three-dimensional geological model that can describe the shape, size, and distribution of underground reservoirs.

[0090] It should be noted that when constructing the initial reservoir model, a three-dimensional geological model was built using geological modeling software and based on collected geological parameter information. This geological parameter information includes parameters such as geological stratigraphy, lithology, porosity, and permeability. This geological parameter information may be obtained from geological exploration, drilling, logging, and seismic data. Furthermore, this initial reservoir model includes fault planes and bedding planes, and the intersections between fault planes, bedding planes, and bedding planes are strictly geometrically and topologically consistent. Here, topological and geometric consistency means that the geological surfaces such as faults and bedding planes are strictly watertight.

[0091] The initial reservoir mesh model is the model after tetrahedral meshing of the initial reservoir model. Tetrahedral meshing is a method of dividing three-dimensional space into multiple tetrahedral units.

[0092] Optionally, the initial reservoir model can be meshed using the meshing function in meshing software or geological modeling software.

[0093] Optionally, the number of tetrahedral meshes generated can be set based on actual needs.

[0094] S202. A preset layer extraction strategy is used to extract sub-layers from the initial reservoir grid model to generate an initial sub-layer model. The initial sub-layer model consists of multiple sub-layers, and each sub-layer includes multiple triangular mesh elements.

[0095] Specifically, after obtaining the initial reservoir mesh model, sub-layers are extracted from the initial reservoir mesh model based on a preset layer extraction strategy, and each sub-layer includes multiple triangular mesh elements. An initial sub-layer model is then generated based on these multiple sub-layers. This step divides the initial reservoir mesh model and extracts multiple sub-layers to facilitate subsequent mesh processing.

[0096] The preset layer extraction strategy is a pre-defined strategy for extracting sub-layers from the initial reservoir grid model.

[0097] Optionally, the preset layer extraction strategy may be set based on geological layers, lithology, or other geological features; or it may be based on a preset thickness and extracted at intervals along the height direction of the initial reservoir grid model; or it may be based on a preset length and extracted at intervals along the length direction of the initial reservoir grid model. Optionally, different extraction methods may be used for different geological strata in the initial reservoir grid model. In this embodiment, the preset layer extraction strategy is not specifically limited.

[0098] The small layer is a cross-section in the initial reservoir grid model, and the triangular mesh cells of the small layer are obtained by cutting the tetrahedral mesh in the initial reservoir grid model.

[0099] The initial small-level model is a 3D model obtained based on all the small levels after the small levels have been extracted. Optionally, multiple small levels retain their original 3D spatial positions.

[0100] Optionally, small layers can be extracted from the initial reservoir grid model using geological modeling software.

[0101] S203. Using a preset mesh processing strategy, triangular mesh cells in multiple small layers are converted into rectangular mesh cells to generate the target small layer model.

[0102] Specifically, after obtaining the initial small-level model, triangular mesh elements in multiple small levels are converted into rectangular mesh elements based on a preset mesh processing strategy, thereby obtaining the target small-level model. All triangular mesh elements are converted into rectangular mesh elements to facilitate the subsequent generation of hexahedral mesh elements.

[0103] The preset mesh processing strategy is a pre-defined strategy for converting triangular mesh cells into rectangular mesh cells. Optionally, this mesh processing strategy may include steps such as re-meshing, interpolation, or resampling. In this embodiment, the preset mesh processing strategy is not limited; it only needs to be able to convert triangular mesh cells into rectangular mesh cells.

[0104] Alternatively, rectangular grid cells can be generated using a rectangular grid generation algorithm, such as quadtree decomposition or a rectangular grid generator.

[0105] Optionally, after the triangular mesh cells are converted into rectangular mesh cells, the boundaries of the rectangular mesh cells are smoothed to reduce discontinuities and noise.

[0106] The target small-level model is a model composed of multiple small levels after the mesh transformation is completed.

[0107] S204. Generate the target reservoir grid model based on the grid coordinate system and the target small-level model; the grid type in the target reservoir grid model is a hexahedral grid.

[0108] Specifically, after obtaining the target sub-level model, a target reservoir grid model is generated based on the grid coordinate system and the target sub-level model, and the grid type of the target reservoir grid model is all hexahedral mesh. This hexahedral mesh is easy to perform numerical simulation and flow simulation, which is beneficial for subsequent reservoir numerical dynamic analysis, remaining oil potential tapping and other studies.

[0109] The grid coordinate system is a three-dimensional coordinate system used to generate the target reservoir grid model in conjunction with the target sub-level model. Optionally, the target sub-level model can be placed in the grid coordinate system to generate grid cylinders, and then the grid cylinders can be merged to obtain the target sub-level model. Alternatively, the target sub-level model can be segmented using the grid coordinate system, and the segmented target sub-level model can be stretched and merged to obtain grid cylinders, which can then be merged to obtain the target sub-level model. Specifically, in this embodiment, no specific limitation is made on how to generate the target reservoir grid model using the grid coordinate system and the target sub-level model.

[0110] Optionally, the grid coordinate system includes an X-axis, a Y-axis, and a Z-axis. Optionally, the X-axis and Y-axis can be in the latitude and longitude direction, and the Z-axis in the height direction. Optionally, one of the X-axis and Y-axis can be set to the direction of oil and gas source, and the Z-axis in the height direction, thus better reflecting the actual situation of the reservoir.

[0111] Optionally, the grid size or scale of the grid coordinate system can be flexibly adjusted according to the size of the actual reservoir model. For example, the grid size along the X-axis is 10 meters, the grid size along the Y-axis is 1 meter, and the grid size along the Z-axis is 10 meters, etc. In this embodiment, the grid size is not specifically limited.

[0112] The target reservoir mesh model is obtained by processing the target sub-layer model, and all mesh types in the target reservoir mesh model are hexahedral meshes. Specifically, in this embodiment, the size of the hexahedral mesh is not limited and can be flexibly set according to actual needs.

[0113] Alternatively, a hexahedral mesh can be generated using a hexahedral mesh generation algorithm based on the target small-level model. Examples include sweeping or mapping methods.

[0114] Optionally, after generating the hexahedral mesh, the mesh is traversed to determine if there are any mesh intersections or interferences. When hexahedral mesh interference exists, adjustments are made to ensure the continuity and integrity of the hexahedral mesh. Optionally, the hexahedral mesh is traversed and checked using a mesh inspection tool.

[0115] The reservoir mesh generation method provided in this embodiment obtains an initial reservoir mesh model, which is obtained by dividing the initial reservoir model into tetrahedral meshes. A preset layer extraction strategy is used to extract sub-layers from the initial reservoir mesh model to generate an initial sub-layer model. The initial sub-layer model consists of multiple sub-layers, and each sub-layer includes multiple triangular mesh elements. A preset mesh processing strategy is used to convert the triangular mesh elements in the multiple sub-layers into rectangular mesh elements to generate a target sub-layer model. A target reservoir mesh model is generated based on the mesh coordinate system and the target sub-layer model; the mesh type in the target reservoir mesh model is all hexahedral mesh. The reservoir mesh generation method in this embodiment refines the initial reservoir mesh model into multiple sub-layers through a preset layer extraction strategy. Then, a preset mesh processing strategy converts all triangular mesh elements in the sub-layers into rectangular mesh elements. Finally, a target reservoir mesh model consisting entirely of hexahedral meshes is obtained based on the mesh coordinate system and the target sub-layer model. By refining the initial reservoir mesh model into multiple smaller layers, even if non-permeable faults and thrust-neck structures exist within the initial reservoir mesh model, mesh processing can be performed effectively. All triangular mesh elements are converted into rectangular mesh elements, avoiding mesh distortion, improving mesh quality, and ensuring the accuracy and reliability of subsequent analyses based on the reservoir mesh model. The mesh generation method for the reservoir model in this application is applicable to both simple, thin-layer reservoirs and complex reservoirs with numerous non-permeable faults and thrust-neck structures, demonstrating broad applicability.

[0116] Among them, thrust-over structures include landforms such as reverse faults, folds, and faults.

[0117] Specifically, by applying the reservoir model mesh generation method of this embodiment, the mesh generation task of reservoir models with various complex terrains can be completed well. The specific mesh generation effects are shown in Figures 6 to 8. For example, Figure 6 shows the mesh generation effect of a reservoir model with complex normal fault intersection structures. Figure 7 shows the mesh generation effect of a reservoir model with multi-level X-type and Y-type complex reverse faults. Figure 8 shows the mesh generation effect of a reservoir model with a complex reverse-thrust fault intersection framework.

[0118] As an optional implementation, based on any of the above embodiments, a preset layer extraction strategy is used to extract sub-layers from the initial reservoir grid model to generate an initial sub-layer model, including: determining the sedimentary mode corresponding to each stratum using the strata represented in the initial reservoir grid model as the dimension; extracting sub-layers of each stratum along the height direction of the initial reservoir grid model based on the sedimentary mode of each stratum and the corresponding preset layer extraction interval value; and generating the initial sub-layer model based on the extracted sub-layers of each stratum.

[0119] Specifically, when using a pre-defined layer extraction strategy to extract sub-layers from an initial reservoir grid model to generate an initial sub-layer model, the following steps are taken: First, the sedimentary patterns corresponding to each stratum are determined using the strata represented in the initial reservoir grid model as dimensions. Then, based on the sedimentary patterns of each stratum and the corresponding pre-defined layer extraction interval values, sub-layers of each stratum are extracted along the height direction of the initial reservoir grid model. Finally, the initial sub-layer model is generated based on the sub-layers extracted from each stratum.

[0120] In the initial reservoir grid model, the strata are characterized by a dimension that determines the corresponding sedimentary pattern for each stratum. Specifically, a corresponding sedimentary pattern is set for each stratum in the initial reservoir grid model. As shown in Figure 9, the sedimentary patterns of the strata include a scale pattern, a parallel top pattern, and a parallel bottom pattern. This sedimentary pattern guides the horizontal projection shape of the extracted bedding planes from the corresponding strata, making the extracted bedding planes more consistent with the actual strata. Since different strata have different sedimentary patterns, bedding planes are extracted from different strata according to different sedimentary patterns to achieve more accurate extraction.

[0121] The preset layer extraction interval is a pre-set value used to measure the distance between two adjacent small layers during the extraction process, for example, extracting a small layer every 10 meters. Optionally, different strata have different thicknesses, therefore, the preset layer extraction interval value is also different for each stratum. Specifically, in this embodiment, the preset layer extraction interval value is not limited, and it can be flexibly set according to actual conditions.

[0122] Optionally, the extracted small layers can be smoothed to reduce geometric discontinuities caused by the segmentation.

[0123] Optionally, attribute data at a smaller level can be interpolated to ensure data continuity and accuracy.

[0124] Optionally, the extracted and processed strata are stacked in depositional order to form an initial stratum model. Furthermore, it is necessary to ensure the correct topological relationships between adjacent strata to guarantee the integrity of the initial stratum model.

[0125] Optionally, before extracting the sub-layers, the distances from the vertices of each tetrahedral unit to the top and bottom interfaces of the stratum are obtained using a discrete smooth interpolation algorithm. This assists in the subsequent accurate extraction of the sub-layers of each stratum, thereby ensuring the accuracy and reliability of the initial sub-layer model.

[0126] The reservoir model mesh generation method provided in this embodiment employs a preset layer extraction strategy to extract sub-layers from an initial reservoir mesh model to generate an initial sub-layer model. This includes: determining the sedimentary pattern corresponding to each stratum using the strata represented in the initial reservoir mesh model as the dimension; extracting sub-layers along the height direction of the initial reservoir mesh model based on the sedimentary patterns of each stratum and the corresponding preset layer extraction interval values; and generating an initial sub-layer model based on the extracted sub-layers. The method in this embodiment, by extracting sub-layers from each stratum based on its sedimentary pattern and preset layer extraction interval values, can obtain an initial sub-layer model that more closely reflects actual geological conditions.

[0127] As an optional implementation, based on any of the above embodiments, a preset mesh processing strategy is used to convert triangular mesh elements in multiple sub-layers into rectangular mesh elements to generate a target sub-layer model. This includes: converting triangular mesh elements in multiple sub-layers into hybrid mesh elements using a preset mesh transformation strategy. The hybrid mesh elements include rectangular mesh elements far from the fault region and non-rectangular mesh elements near the fault region. A preset mesh regularization strategy is used to perform mesh regularization processing on the non-rectangular mesh elements to generate the target sub-layer model. Multiple sub-layers of the target sub-layer model are all composed of rectangular mesh elements.

[0128] Specifically, when using a preset mesh processing strategy to convert triangular mesh elements in multiple sub-layers into rectangular mesh elements to generate the target sub-layer model, firstly, a preset mesh transformation strategy is used to convert triangular mesh elements in multiple sub-layers into hybrid mesh elements, as shown in Figure 10. Figure 10 is a schematic diagram of any sub-layer with hybrid mesh elements. Then, a preset mesh regularization strategy is used to perform mesh regularization processing on the non-rectangular mesh elements to generate the target sub-layer model composed entirely of rectangular mesh elements.

[0129] Specifically, the preset mesh processing strategies include preset mesh transformation strategies and preset mesh regularization strategies.

[0130] The preset mesh conversion strategy is a pre-defined strategy for converting small-level triangular mesh cells into hybrid mesh cells.

[0131] It should be noted that the hybrid mesh element includes rectangular mesh elements far from the fault region and non-rectangular mesh elements close to the fault region. Specifically, in areas far from the fault region, all triangular mesh elements are converted to rectangular mesh elements using a preset mesh transformation strategy, while in areas close to or far from the fault region, all triangular mesh elements are converted to hybrid mesh elements using the same preset strategy. Furthermore, after mesh regularization, the mesh of the smaller layers in the fault region is arranged in a stepped pattern.

[0132] Optionally, triangular mesh cells can be converted into hybrid mesh cells using rectangular mesh cells in the mesh coordinate system. For example, by comparing rectangular mesh cells with triangular mesh cells, a triangular mesh cell is converted into a rectangular mesh cell when all four corners of the rectangular mesh cell fall within a triangular mesh cell. Optionally, a mesh transformation algorithm can be used to convert triangular mesh cells into rectangular mesh cells, such as quadtree decomposition or a rectangular mesh generator. Specifically, in this embodiment, no specific limitation is made on how to convert triangular mesh cells into hybrid mesh cells.

[0133] Optionally, before converting triangular mesh cells into hybrid mesh cells, triangular mesh cells in multiple sub-layers are analyzed, including their geometry, topological relationships, and coordinate information.

[0134] The preset mesh regularization strategy is a pre-defined strategy for regularizing non-rectangular mesh elements, that is, converting all non-rectangular mesh elements into rectangular mesh elements. Converting all non-rectangular mesh elements into rectangular mesh elements facilitates the subsequent generation of high-quality hexahedral mesh elements, thereby ensuring the accuracy and precision of subsequent numerical simulations based on the reservoir model.

[0135] Optionally, non-rectangular grid cells include, but are not limited to, trapezoids, parallelograms, triangles, pentagons, hexagons, etc.

[0136] Optionally, non-rectangular mesh cells can be converted to rectangular mesh cells using a mesh coordinate system. For example, if the area of ​​a non-rectangular mesh cell exceeds half the area of ​​a rectangular mesh cell, the non-rectangular mesh cell is replaced with a rectangular mesh cell; otherwise, the non-rectangular mesh cell is deleted. Alternatively, all non-rectangular mesh cells can be converted to rectangular mesh cells using re-division algorithms, mesh optimization algorithms, etc., which may involve steps such as re-dividing the mesh cells, adjusting their shapes, and interpolating attribute data. This embodiment does not specifically limit the method for converting non-rectangular mesh cells to rectangular mesh cells.

[0137] Optionally, after converting all triangular mesh cells in the small layer into rectangular mesh cells, the coordinate information of the rectangular mesh cells is traversed to determine whether there are any discontinuities, so as to ensure the continuity and accuracy of the generated target small layer model.

[0138] It should be noted that due to the presence of faults or reverse faults in the initial reservoir model, there are certain height differences at the fault or reverse fault locations within the same stratum. Therefore, after extracting the sub-layers, fault regions will exist within these sub-layers, and discontinuities such as fractures or breaks will be present within these fault regions.

[0139] The target small-level model is obtained by processing the initial small-level model through a preset mesh transformation strategy and a preset mesh regularization strategy. All mesh types in the small levels of the target small-level model are rectangular meshes.

[0140] The reservoir model mesh generation method provided in this embodiment includes the following steps when converting triangular mesh elements in multiple sub-layers into rectangular mesh elements using a preset mesh processing strategy to generate a target sub-layer model: First, a preset mesh transformation strategy is used to convert triangular mesh elements in multiple sub-layers into hybrid mesh elements. The hybrid mesh elements include rectangular mesh elements far from the fault region and non-rectangular mesh elements near the fault region. Then, a preset mesh regularization strategy is used to regularize the non-rectangular mesh elements to generate the target sub-layer model. All sub-layers of the target sub-layer model are composed of rectangular mesh elements. The reservoir model mesh generation method in this embodiment converts all triangular mesh elements in the sub-layers into rectangular mesh elements through preset mesh transformation and mesh regularization strategies. This facilitates the subsequent generation of high-quality hexahedral mesh elements, avoids mesh distortion, ensures mesh quality, and thus guarantees the accuracy and precision of subsequent numerical simulations based on the reservoir model.

[0141] As an optional implementation, based on any of the above embodiments, a preset mesh transformation strategy is used to convert triangular mesh units in multiple small layers into hybrid mesh units, including: dividing the triangular mesh units of the small layers into hybrid mesh units using a unit rectangular mesh based on a mesh coordinate system. The unit rectangular mesh of the mesh coordinate system is formed by the intersection of multiple horizontal grid lines and multiple horizontal grid lines.

[0142] Specifically, when converting triangular mesh cells in multiple small layers into hybrid mesh cells using a preset mesh transformation strategy, firstly, the triangular mesh cells in the small layers are divided using a unit rectangular mesh based on the mesh coordinate system to obtain hybrid mesh cells. This unit rectangular mesh is formed by the intersection of multiple horizontal and vertical mesh lines. The triangular mesh cells are divided into rectangular mesh cells and non-rectangular mesh cells using the unit rectangular mesh. Specifically, the division of triangular mesh cells in the small layers using the unit rectangular mesh in the mesh coordinate system is shown in Figure 11. Figure 11 is a schematic diagram of the triangular mesh cell division.

[0143] Specifically, in this embodiment, a grid coordinate system is established with the X-axis parallel to the oil and gas direction, the Y-axis perpendicular to the X-axis, and the Z-axis representing the height direction of the reservoir model. Horizontal grid lines are parallel to the X-axis, and horizontal grid lines are parallel to the Y-axis.

[0144] Optionally, the area of ​​the unit rectangular grid can be adjusted so that the area of ​​the unit rectangular grid is smaller than that of the triangular grid unit. This divides the triangular grid units into smaller, more regular rectangular grid units and non-rectangular grid units using the unit rectangular grid. For example, if the coordinates of all four corners of a unit rectangular grid are located inside one or more triangular grid units, the unit rectangular grid is retained; otherwise, the original grid shape is retained. Alternatively, the area of ​​the unit rectangular grid can be set to be larger than that of the triangular grid unit, thereby converting multiple triangular grid units into larger rectangular grid units.

[0145] Optionally, an algorithm is used to compare and analyze the triangular mesh cells with the unit rectangular mesh and then cut them. This process may involve subdividing, merging, or adjusting the shape of the triangular mesh cells to ensure that the cut mesh cells match the unit rectangular mesh.

[0146] The reservoir model mesh generation method provided in this embodiment employs a preset mesh transformation strategy to convert triangular mesh units in multiple sub-layers into hybrid mesh units. This includes: dividing the triangular mesh units in the sub-layers using a unit rectangular mesh based on a mesh coordinate system to obtain hybrid mesh units. The unit rectangular mesh in the mesh coordinate system is formed by the intersection of multiple horizontal and vertical mesh lines. The method in this embodiment, by dividing the triangular mesh units based on a mesh coordinate system, can conveniently and efficiently convert triangular mesh units into hybrid mesh units.

[0147] As an optional implementation, based on any of the above embodiments, a preset mesh processing strategy is used to perform mesh regularization processing on non-rectangular mesh cells to generate a target small-scale model, including: determining whether a non-rectangular mesh cell covers the center point of a unit rectangular mesh in the mesh coordinate system. In response to a non-rectangular mesh cell covering the center point of a unit rectangular mesh in the mesh coordinate system, the non-rectangular mesh cell is replaced with a unit rectangular mesh to generate the target small-scale model.

[0148] Specifically, when using a preset mesh processing strategy to regularize non-rectangular mesh cells to generate the target small-level model, the process first determines whether the non-rectangular mesh cell covers the center point of the unit rectangular mesh in the mesh coordinate system. If so, the non-rectangular mesh cell is considered a valid mesh that needs to be retained, and it is replaced with a unit rectangular mesh cell. This allows for the retention of non-rectangular mesh cells with larger areas, reducing errors in subsequent reservoir models.

[0149] It should be noted that the center point of each unit rectangular grid can be calculated based on the coordinates of the four corners of each unit rectangular grid.

[0150] Alternatively, the coverage can be determined by comparing the boundary coordinates of the non-rectangular grid cells with the coordinates of the center point of the unit rectangular grid.

[0151] Optionally, when it is determined that a non-rectangular mesh cell does not cover the center point of a unit rectangular mesh in the mesh coordinate system, the non-rectangular mesh cell can be directly deleted, or it can be further modified into a smaller rectangular mesh. In this embodiment, this is not limited.

[0152] Optionally, in other embodiments, it can also be determined whether to replace it with a unit rectangular grid by calculating the area of ​​the non-rectangular grid cell and the area of ​​the corresponding unit rectangular grid. For example, when the area of ​​the non-rectangular grid cell exceeds 1 / 2 or 1 / 3 of the area of ​​the unit rectangular grid, the non-rectangular grid cell is replaced by the unit rectangular grid of the grid coordinate system.

[0153] Optionally, in other embodiments, rectangular grid cells with the same area as the unit rectangular grid can be obtained by expanding the non-rectangular grid cells outward. Further, when the rectangular grid cell covers the center point of the unit rectangular grid in the grid coordinate system, the non-rectangular grid cell is converted into a rectangular grid cell using interpolation extrapolation based on the coordinate information of the four corners of the unit rectangular grid.

[0154] The reservoir model mesh generation method provided in this embodiment employs a preset mesh processing strategy to regularize non-rectangular mesh cells to generate a target small-scale model. This includes: determining whether a non-rectangular mesh cell covers the center point of a unit rectangular mesh in the mesh coordinate system; and replacing the non-rectangular mesh cell with a unit rectangular mesh in response to the non-rectangular mesh cell covering the center point of the unit rectangular mesh in the mesh coordinate system to generate the target small-scale model. This method, by determining and replacing non-rectangular mesh cells, achieves mesh regularization, thereby generating a target small-scale model with a high-quality mesh.

[0155] As an optional implementation, based on any of the above embodiments, it further includes: deleting the non-rectangular mesh cell in response to the non-rectangular mesh cell not covering the center point of the unit rectangular mesh of the mesh coordinate system.

[0156] Specifically, when it is determined that a non-rectangular mesh cell does not cover the center point of the unit rectangular mesh in the mesh coordinate system, the non-rectangular mesh cell is deleted. This ensures that the converted mesh consists entirely of rectangular mesh cells and reduces the difficulty of mesh conversion.

[0157] It should be noted that non-rectangular mesh cells that do not cover the center point of the unit rectangular mesh are considered to have less information value, and therefore deleting them will not have a significant impact on the accuracy of the target small-scale model. These deleted non-rectangular mesh cells will be removed from the model and will no longer participate in subsequent simulations and analyses.

[0158] Optionally, after deleting non-rectangular mesh cells, geometric and topological verifications are performed on the model to ensure the continuity of the obtained target small-level model.

[0159] The reservoir model mesh generation method provided in this embodiment further includes: deleting non-rectangular mesh elements in response to the non-rectangular mesh elements not covering the center point of the unit rectangular mesh in the mesh coordinate system. This method, by deleting non-rectangular mesh elements that do not cover the center point of the unit rectangular mesh, ensures that the converted meshes are all rectangular meshes, reducing the difficulty of mesh conversion.

[0160] As an optional implementation, based on any of the above embodiments, after replacing the non-rectangular mesh cells with unit rectangular meshes, the method further includes: determining whether the vertex coordinates of adjacent sides of adjacent rectangular mesh cells are consistent. In response to the inconsistency of vertex coordinates of adjacent sides of adjacent rectangular mesh cells, the vertices of adjacent sides of adjacent rectangular mesh cells are merged.

[0161] Specifically, after replacing non-rectangular mesh elements with unit rectangular meshes, geometric discontinuities may exist between the rectangular mesh elements in the fault region. Therefore, to address this issue, after replacing the non-rectangular mesh elements with unit rectangular meshes, it is determined whether the vertex coordinates of adjacent edges of adjacent rectangular mesh elements are consistent. If the vertex coordinates of adjacent edges of adjacent rectangular mesh elements are inconsistent, a geometric inconsistency exists. In this case, the vertices of adjacent edges of adjacent rectangular mesh elements are merged, thereby resolving the geometric inconsistency between adjacent rectangular mesh elements. Specifically, the geometric inconsistency of adjacent rectangular mesh elements is illustrated in Figure 12. Figure 12 is a schematic diagram of geometric inconsistency among small-scale rectangular mesh elements.

[0162] Optionally, an algorithm can be used to compare whether the vertex coordinates of adjacent sides of each pair of adjacent rectangular grid cells are consistent. If the vertex coordinates are completely consistent, it means that the adjacent rectangular grid cells are geometrically continuous; if they are inconsistent, there may be geometric discontinuities or misalignments.

[0163] Optionally, when merging adjacent rectangular grid cells, the coordinates of the merged vertex are obtained by averaging or other suitable calculation methods based on the coordinates of the adjacent vertices. For example, the average of the coordinates of the two vertices can be taken as the coordinates of the merged point. The merged vertex coordinates are then updated in the grid data, along with other grid information related to that vertex, such as the side length and area of ​​the grid cell.

[0164] Optionally, after merging the vertices of adjacent rectangular grid cells, the merged rectangular grid cells are processed by interpolation or smoothing to ensure the continuity and accuracy of the data.

[0165] Optionally, after merging vertices of adjacent rectangular mesh cells, the geometric information of the target small-level model is updated, including the size, shape, and adjacency of the rectangular mesh cells.

[0166] The reservoir model mesh generation method provided in this embodiment, after replacing non-rectangular mesh cells with unit rectangular meshes, further includes: determining whether the vertex coordinates of adjacent edges of adjacent rectangular mesh cells are consistent. In response to inconsistent vertex coordinates of adjacent edges of adjacent rectangular mesh cells, the vertices of adjacent edges of adjacent rectangular mesh cells are merged. The method of this embodiment can effectively detect and resolve geometric inconsistencies between adjacent rectangular mesh cells.

[0167] As an optional implementation, based on any of the above embodiments, a target reservoir grid model is generated based on a grid coordinate system and a target sub-level model, including: generating multiple unit grid columns based on the grid coordinate system and the target sub-level model, wherein each unit grid column includes multiple hexahedral grids arranged along the extension direction of the unit grid column. The multiple unit grid columns are then merged to generate the target reservoir grid model.

[0168] Specifically, when generating the target reservoir mesh model based on the grid coordinate system and the target sub-level model, firstly, multiple element mesh columns are generated based on the grid coordinate system and the target sub-level model. Then, these multiple element mesh columns are merged to obtain a target reservoir mesh model where all mesh types are hexahedral.

[0169] Among them, the unit grid column is a columnar structure that includes multiple hexahedral grids.

[0170] Optionally, there are several ways to generate multiple unit mesh columns based on the mesh coordinate system and the target small-level model. For example, a column can be generated based on the mesh lines of the three dimensions of the mesh coordinate system, and then the rectangular mesh cells of the small levels in the target small-level model can be placed sequentially and spaced within the column to form a unit mesh column. Alternatively, multiple small levels can be divided using the mesh coordinate lines of a certain dimension of the mesh coordinate system to form a unit mesh column. In this embodiment, no specific limitation is made on how to generate multiple unit mesh columns based on the mesh coordinate system and the target small-level model.

[0171] Optionally, adjacent cell grid columns can be merged based on their endpoints, midpoints, or other locations.

[0172] The reservoir model mesh generation method provided in this embodiment generates a target reservoir mesh model based on a mesh coordinate system and a target sub-level model. This includes: generating multiple unit mesh columns based on the mesh coordinate system and the target sub-level model, where each unit mesh column comprises multiple hexahedral meshes arranged along its extension direction. The multiple unit mesh columns are then merged to generate the target reservoir mesh model. This method, by generating multiple unit mesh columns containing hexahedral meshes using a mesh coordinate system and a target sub-level model, and then merging these multiple unit mesh columns, can easily obtain the target reservoir mesh model.

[0173] As an optional implementation, based on any of the above embodiments, multiple unit mesh columns are generated based on a mesh coordinate system and the target small-scale model, including: generating multiple initial mesh columns based on horizontal grid lines, horizontal grid lines, and height grid lines of the mesh coordinate system. The target small-scale model is then segmented based on the initial mesh columns to generate multiple unit mesh columns.

[0174] Specifically, when generating multiple unit mesh columns based on the mesh coordinate system and the target small-scale model, firstly, multiple initial mesh columns are generated in three-dimensional space based on the horizontal grid lines, horizontal grid lines, and height grid lines of the mesh coordinate system. Then, the target small-scale model is segmented based on the initial mesh columns to obtain multiple unit mesh columns. Specifically, Figure 13 is a schematic diagram of the initial mesh column provided in an embodiment of this application; Figure 14 is a schematic diagram of the unit mesh column provided in an embodiment of this application. The schematic diagram of the initial mesh column generation is shown in Figure 13; the schematic diagram of the unit mesh column is shown in Figure 14.

[0175] Specifically, in this embodiment, the target small-level model is placed in a grid coordinate system such that the initial grid cylinder passes through the small level and through the four corners of the rectangular grid cell. Thus, the rectangular grid cells in two adjacent small levels form two faces of a hexahedral grid, while the other four faces are formed by the sides of the initial grid cylinder. In this embodiment, the height grid line passes through the four corners of the rectangular grid cell as an example.

[0176] It should be noted that when generating the initial mesh cylinder, the spacing of the horizontal grid lines, the horizontal grid lines, and the height grid lines needs to be set appropriately so that the initial mesh cylinder can match the rectangular grid cells in the smaller layers, thereby ensuring the generation effect of the hexahedral mesh.

[0177] The reservoir model mesh generation method provided in this embodiment generates multiple element mesh columns based on a mesh coordinate system and a target sub-level model. This includes generating multiple initial mesh columns based on horizontal transverse mesh lines, horizontal longitudinal mesh lines, and height mesh lines within the mesh coordinate system. The target sub-level model is then segmented based on these initial mesh columns to generate more element mesh columns. This method, by using a mesh coordinate system to generate initial mesh columns and then segmenting the target sub-level model to obtain element mesh columns, can simply and efficiently produce element mesh columns with multiple hexahedral meshes.

[0178] As an optional implementation, based on any of the above embodiments, after generating multiple unit grid columns based on the grid coordinate system and the target sub-layer model, the method further includes: determining whether a reverse fault exists within the unit grid column. If no reverse fault exists within the unit grid column, the target sub-layers within the unit grid column are arranged according to the sedimentary order of the strata. If a reverse fault exists within the unit grid column, the target sub-layers within the unit grid column are arranged from highest to lowest according to their height coordinate values.

[0179] Specifically, since reverse faults may exist in the reservoir model, to ensure the orderly arrangement of target sub-layers within the unit grid column, after generating multiple unit grid columns based on the grid coordinate system and the target sub-layer model, it is necessary to determine whether a reverse fault exists within the unit grid column. If not, the target sub-layers within the unit grid column are arranged according to the sedimentary order of the strata; if so, the target sub-layers within the unit grid column are arranged from high to low according to their height coordinate values. Specifically, Figure 15 is a schematic diagram of the sub-layer arrangement within the unit grid column provided in the embodiment of this application, and the arrangement diagram of the sub-layers is shown in Figure 15. As shown in the second small figure in Figure 15, if a reverse fault exists within the unit grid column, and the sub-layers are still arranged according to the original sedimentary order, it will cause a disordered arrangement that cannot reflect the true strata situation.

[0180] Among them, reverse faults are a type of fault in geological structures, mainly formed by horizontal compression and gravity. They refer to faults in which the hanging wall slides upward along the fault plane while the footwall relatively subsides.

[0181] It should be noted that the presence of reverse faults causes vertical displacement and overlap of strata, and arranging them according to the traditional sedimentary sequence cannot accurately reflect the true situation of the strata. However, arranging them according to their height coordinates from high to low can more accurately reflect the impact of reverse faults on the strata, providing a more reliable basis for subsequent simulations and analyses.

[0182] It should be noted that the initial reservoir model was established based on known geological exploration data, including specific information about reverse faults, such as the three-dimensional coordinate range of the reverse faults. Therefore, by comparing the horizontal and vertical coordinate values ​​of the unit grid column with the three-dimensional coordinate range of the reverse fault, it is possible to determine whether a reverse fault exists within the unit grid column.

[0183] The reservoir model mesh generation method provided in this embodiment, after generating multiple unit mesh columns based on the mesh coordinate system and the target sub-layer model, further includes: determining whether a reverse fault exists within the unit mesh column. If no reverse fault exists within the unit mesh column, the target sub-layers within the unit mesh column are arranged according to the sedimentary sequence of the strata. If a reverse fault exists within the unit mesh column, the target sub-layers within the unit mesh column are arranged from high to low according to their height coordinate values. This method, by arranging the target sub-layers according to the sedimentary sequence or height coordinate values ​​based on fault conditions, ensures the generation effect of the unit mesh columns, effectively improving the accuracy and reliability of reservoir simulation, and providing strong support for subsequent oil and gas exploration and development.

[0184] As an optional implementation, based on any of the above embodiments, the method further includes: determining whether there are intersecting target sub-layers within the unit grid column. In response to determining that there are intersecting target sub-layers within the unit grid column, the target sub-layers within the unit grid column are rearranged from high to low according to their height coordinate values.

[0185] Specifically, Figure 16 is a schematic diagram of the adjustment of intersecting small layers of the unit mesh column provided in the embodiment of this application. After generating the unit mesh column, there may still be cases of small layer intersections due to the arrangement of small layers, as shown in the left figure of Figure 16. Therefore, it is necessary to determine whether there are intersecting target small layers within the unit mesh column. When it is determined that there are intersecting target small layers within the unit mesh column, it indicates that there is a problem with the arrangement of small layers. Therefore, the target small layers within the unit mesh column are rearranged from high to low according to their height coordinate values. Specifically, the effect of rearranging the target small layers within the unit mesh column is shown in the right figure of Figure 16.

[0186] It should be noted that the presence of reverse faults may cause vertical overlap between adjacent sub-level grid cells within a single grid column (i.e., sub-level intersection). Optionally, the height coordinates of adjacent sub-levels are sequentially determined from high to low within each grid column to ascertain whether they intersect. For example, if the height coordinate value of an upper sub-level is lower than that of a lower sub-level, it indicates that the sub-levels intersect.

[0187] The reservoir model mesh generation method provided in this embodiment further includes: determining whether there are intersecting target sub-layers within a unit mesh column. In response to determining that there are intersecting target sub-layers within a unit mesh column, the target sub-layers within the unit mesh column are rearranged according to their height coordinate values ​​from high to low. The method in this embodiment checks each unit mesh column to confirm whether there are intersecting target sub-layers. When intersecting target sub-layers are found, the target sub-layers within that unit mesh column are rearranged to avoid errors or intersections in the subsequently generated hexahedral mesh.

[0188] As an optional implementation, based on any of the above embodiments, after generating multiple unit mesh columns based on the mesh coordinate system and the target small-level model, the method further includes: adding a three-dimensional coordinate index to the hexahedral mesh within the unit mesh column.

[0189] Specifically, after generating multiple cell grid columns based on the grid coordinate system and the target sub-layer model, it is necessary to add a three-dimensional coordinate index to the hexahedral mesh within each cell grid column. The three-dimensional coordinate index provides precise spatial location information for each hexahedral mesh, enabling rapid location of a specific hexahedral mesh during subsequent analysis and simulation.

[0190] Optionally, when adding a three-dimensional coordinate index to a hexahedral mesh, the index is determined based on its position in three-dimensional space. This three-dimensional coordinate index may include the vertex coordinates, corner coordinates, or center point coordinates of the hexahedral mesh.

[0191] Optionally, the determined three-dimensional coordinate index information can be stored for quick access and reference in subsequent simulation and analysis processes.

[0192] Optionally, after adding the 3D coordinate index, check whether the added 3D coordinate index is accurate and consistent, specifically including verifying whether the index information matches the actual position of the hexahedral mesh, and whether there are conflicts or redundancies between the indexes.

[0193] The reservoir model mesh generation method provided in this embodiment, after generating multiple element mesh columns based on the mesh coordinate system and the target small-level model, further includes adding a three-dimensional coordinate index to the hexahedral mesh within the element mesh column. This method, by adding a three-dimensional coordinate index to the hexahedral mesh, facilitates rapid location of a specific hexahedral mesh during subsequent analysis and simulation.

[0194] As an optional implementation, based on any of the above embodiments, the method further includes: determining whether there is index duplication in the three-dimensional coordinate indices of the hexahedral mesh within the unit mesh column. In response to the presence of index duplication, adjusting the coordinate index value of any dimension of the hexahedral mesh with the duplicated index.

[0195] Specifically, since reverse faults may exist in the initial reservoir model, these faults may lead to the same 3D coordinate indices for different hexahedral meshes within a unit grid column. Therefore, after adding 3D coordinate indices to the hexahedral meshes, it is necessary to determine whether there are duplicate 3D coordinate indices for the hexahedral meshes within the unit grid column. If duplicate indices exist, the coordinate index value of any dimension of the hexahedral mesh with duplicate indices is adjusted. This ensures that the 3D coordinate index of each hexahedral mesh is unique, thereby ensuring the accuracy of subsequent reservoir model analysis. Specifically, Figure 17 is a schematic diagram of hexahedral mesh re-indexing within a unit grid column provided in the embodiment of this application. Referring to Figure 17, Figure 17 is a schematic diagram of hexahedral mesh re-indexing.

[0196] Specifically, in this embodiment, by comparing the height coordinate index values ​​of the hexahedral mesh within each unit mesh column, if the height coordinate indexes of two different hexahedral meshes are the same, it indicates that there is a case of index duplication.

[0197] Optionally, once duplicate indices are identified, the horizontal x-coordinate index value, horizontal y-coordinate index value, or height coordinate index value of the hexahedral mesh with the duplicate index is adjusted. For example, the coordinate index value can be adjusted by increasing or decreasing a single value.

[0198] Specifically, in this embodiment, adjusting the height coordinate index value of a repeatedly indexed hexahedral mesh is taken as an example. For instance, if the height coordinate index value of the third and fourth hexahedral meshes along the height direction is both 3, then the height coordinate index value of the fourth and subsequent n hexahedral meshes is increased by 3, and the height coordinate index value of the fourth hexahedral mesh becomes 6. This continues until the height coordinate index value of the fifth hexahedral mesh is 7. This ensures that the 3D coordinate index of each hexahedral mesh within a unit mesh column is unique. The specific reindexing result is shown in the right figure of Figure 17, where it is clearly evident that each hexahedral mesh has a unique index number.

[0199] Optionally, the index repeatability is rechecked after the index value is modified to ensure that the adjusted index is unique and does not introduce new duplicates.

[0200] The reservoir model mesh generation method provided in this embodiment further includes: determining whether there are duplicate three-dimensional coordinate indices in the hexahedral meshes within a unit mesh column. In response to the presence of duplicate indices, the coordinate index value of any dimension of the hexahedral mesh with duplicate indices is adjusted. The method of this embodiment ensures that the three-dimensional coordinate index of each hexahedral mesh is unique, thereby ensuring the accuracy of subsequent analyses based on the reservoir model.

[0201] Figure 18 is a schematic diagram of the structure of a reservoir model mesh generation device provided in an embodiment of this application. As shown in Figure 18, the reservoir model mesh generation device provided in this embodiment is located in an electronic device. The reservoir model mesh generation device 30 provided in this embodiment includes: an acquisition module 31, an extraction module 32, a conversion module 33, and a model generation module 34.

[0202] Specifically, module 31 is used to acquire an initial reservoir mesh model, which is obtained by dividing the initial reservoir model into tetrahedral meshes. Module 32 is used to extract sub-layers from the initial reservoir mesh model using a preset layer extraction strategy to generate an initial sub-layer model. The initial sub-layer model consists of multiple sub-layers, and each sub-layer includes multiple triangular mesh elements. Module 33 is used to convert the triangular mesh elements in the multiple sub-layers into rectangular mesh elements using a preset mesh processing strategy to generate a target sub-layer model. Module 34 is used to generate a target reservoir mesh model based on the mesh coordinate system and the target sub-layer model. All mesh types in the target reservoir mesh model are hexahedral meshes.

[0203] Optionally, when extracting sub-layers from the initial reservoir grid model using a preset layer extraction strategy to generate an initial sub-layer model, the extraction module 32 is specifically used for: determining the sedimentary pattern corresponding to each stratum using the strata represented in the initial reservoir grid model as the dimension; extracting sub-layers of each stratum along the height direction of the initial reservoir grid model based on the sedimentary pattern of each stratum and the corresponding preset layer extraction interval value; and generating an initial sub-layer model based on the extracted sub-layers of each stratum.

[0204] Optionally, when the conversion module 33 converts triangular mesh elements in multiple sub-layers into rectangular mesh elements using a preset mesh processing strategy to generate a target sub-layer model, it specifically performs the following: It converts triangular mesh elements in multiple sub-layers into hybrid mesh elements using a preset mesh conversion strategy. The hybrid mesh elements include rectangular mesh elements far from the fault region and non-rectangular mesh elements near the fault region. A preset mesh regularization strategy is then used to perform mesh regularization processing on the non-rectangular mesh elements to generate the target sub-layer model. All sub-layers of the target sub-layer model are composed of rectangular mesh elements.

[0205] Optionally, when the conversion module 33 converts triangular mesh units in multiple small layers into hybrid mesh units using a preset mesh conversion strategy, it specifically performs the following: It divides the triangular mesh units in the small layers using a unit rectangular mesh based on the mesh coordinate system to obtain hybrid mesh units. The unit rectangular mesh in the mesh coordinate system is formed by the intersection of multiple horizontal grid lines and multiple vertical grid lines.

[0206] Optionally, when the conversion module 33 performs mesh regularization processing on non-rectangular mesh cells using a preset mesh processing strategy to generate a target small-scale model, it is specifically used to: determine whether the non-rectangular mesh cells cover the center point of the unit rectangular mesh in the mesh coordinate system. In response to the non-rectangular mesh cells covering the center point of the unit rectangular mesh in the mesh coordinate system, the non-rectangular mesh cells are replaced with unit rectangular meshes to generate the target small-scale model.

[0207] Optionally, the conversion module 33 is also configured to: delete the non-rectangular mesh cell in response to the non-rectangular mesh cell not covering the center point of the unit rectangular mesh in the mesh coordinate system.

[0208] Optionally, after replacing non-rectangular mesh cells with unit rectangular mesh cells, the conversion module 33 is further configured to: determine whether the vertex coordinates of adjacent sides of adjacent rectangular mesh cells are consistent. In response to the inconsistency of vertex coordinates of adjacent sides of adjacent rectangular mesh cells, the vertices of adjacent sides of adjacent rectangular mesh cells are merged.

[0209] Optionally, the model generation module 34, when generating the target reservoir grid model based on the grid coordinate system and the target sub-level model, is specifically used to: generate multiple unit grid columns based on the grid coordinate system and the target sub-level model, wherein each unit grid column includes multiple hexahedral grids arranged along the extension direction of the unit grid column. The multiple unit grid columns are then merged to generate the target reservoir grid model.

[0210] Optionally, the model generation module 34, when generating multiple unit mesh columns based on the mesh coordinate system and the target sub-level model, specifically generates multiple initial mesh columns based on the horizontal grid lines, horizontal grid lines, and height grid lines of the mesh coordinate system. The target sub-level model is then segmented based on these initial mesh columns to generate multiple unit mesh columns.

[0211] Optionally, after generating multiple unit grid columns based on the grid coordinate system and the target sub-layer model, the model generation module 34 is further configured to: determine whether a reverse fault exists within the unit grid column; in response to the absence of a reverse fault within the unit grid column, arrange the target sub-layers within the unit grid column according to the sedimentary order of the strata; in response to the presence of a reverse fault within the unit grid column, arrange the target sub-layers within the unit grid column from high to low according to their height coordinate values.

[0212] Optionally, the model generation module 34 is further configured to: determine whether there are intersecting target sub-layers within the cell grid column. In response to determining that there are intersecting target sub-layers within the cell grid column, the target sub-layers within the cell grid column are rearranged from high to low according to their height coordinate values.

[0213] Optionally, after generating multiple unit mesh columns based on the mesh coordinate system and the target small-level model, the model generation module 34 is also used to add three-dimensional coordinate indexes to the hexahedral mesh within the unit mesh columns.

[0214] Optionally, the model generation module 34 is further configured to: determine whether there are duplicate three-dimensional coordinate indices in the hexahedral mesh within the unit mesh column; and adjust the coordinate index value of any dimension of the hexahedral mesh with duplicate indices in response to the existence of duplicate indices.

[0215] It should be noted that the apparatus provided in this application embodiment can implement all the method steps implemented in the above method embodiment and can achieve the same technical effect. Here, the parts that are the same as those in the method embodiment and the beneficial effects will not be described in detail.

[0216] Figure 19 is a schematic diagram of the structure of an electronic device provided in an embodiment of this application. As shown in Figure 19, the electronic device 40 provided in this embodiment includes: a processor 41 and a memory 42 that is communicatively connected to the processor.

[0217] The memory 42 stores computer-executed instructions; the processor 41 executes the computer-executed instructions stored in the memory 42 to implement the reservoir model mesh generation method provided in any of the above embodiments.

[0218] The program may include program code, which includes computer-executable instructions. Memory 42 may include high-speed RAM, and may also include non-volatile memory, such as at least one disk storage device.

[0219] In this embodiment, the memory 42 and the processor 41 are connected via a bus. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus, etc. Buses can be categorized as address buses, data buses, control buses, etc. For ease of illustration, only a single straight line is used in Figure 19, but this does not imply that there is only one bus or one type of bus.

[0220] This application also provides a computer-readable storage medium, including computer-executable instructions stored therein. When executed by a processor, these instructions are used to implement the mesh generation method for the reservoir model provided in any of the above embodiments. For example, the computer-readable storage medium may be a ROM, random access memory (RAM), CD-ROM, magnetic tape, floppy disk, and optical data storage device. Non-transitory computer-readable storage media may be any available medium or data storage device accessible to the processor, including but not limited to magnetic memory (e.g., floppy disk, hard disk, magnetic tape, magneto-optical disk (MO), etc.), optical memory (e.g., CD, DVD, BD, HVD, etc.), and semiconductor memory (e.g., ROM, EPROM, EEPROM, non-volatile memory (NAND FLASH), solid-state drive (SSD), etc.).

[0221] This application also provides a computer program product, including a computer program that, when executed by a processor, implements the mesh generation method for the reservoir model provided in any of the above embodiments.

[0222] In the several embodiments provided in this application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative; for instance, the division of modules is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple modules may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be indirect coupling or communication connection through some interfaces, devices, or modules, and may be electrical, mechanical, or other forms.

[0223] The modules described as separate components may or may not be physically separate. The components shown as modules may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to implement the solution of this embodiment according to actual needs.

[0224] Furthermore, the functional modules in the various embodiments of this application can be integrated into one processing unit, or each module can exist physically separately, or two or more modules can be integrated into one unit. The unit composed of the above modules can be implemented in hardware or in the form of hardware plus software functional units.

[0225] The integrated modules described above, implemented as software functional modules, can be stored in a computer-readable storage medium. These software functional modules, stored in a storage medium, include several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) or processor to execute some steps of the methods of the various embodiments of this application.

[0226] It should be understood that the aforementioned processor can be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), etc. A general-purpose processor can be a microprocessor or any conventional processor. The steps of the disclosed method can be directly manifested as execution by a hardware processor, or execution by a combination of hardware and software modules within the processor.

[0227] The memory may include high-speed RAM, and may also include non-volatile storage (NVM), such as at least one disk storage device, and may also be a USB flash drive, external hard drive, read-only memory, disk or optical disc, etc.

[0228] The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus, etc. Buses can be categorized as address buses, data buses, control buses, etc. For ease of illustration, the buses shown in the accompanying drawings are not limited to a single bus or a single type of bus.

[0229] The aforementioned storage medium can be implemented from any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic storage, flash memory, magnetic disk, or optical disk. The storage medium can be any available medium accessible to general-purpose or special-purpose computers.

[0230] An exemplary storage medium is coupled to a processor, enabling the processor to read information from and write information to the storage medium. Alternatively, the storage medium can be an integral part of the processor. The processor and storage medium can reside in an Application Specific Integrated Circuit (ASIC). Alternatively, the processor and storage medium can exist as discrete components in an electronic control unit or main control device.

[0231] Those skilled in the art will understand that all or part of the steps of the above-described method embodiments can be implemented by hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium. When executed, the program performs the steps of the above-described method embodiments; and the aforementioned storage medium includes various media capable of storing program code, such as ROM, RAM, magnetic disks, or optical disks.

[0232] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of this application, and are not intended to limit them. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features therein. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of this application.

Claims

1. A method for generating a mesh for a reservoir model, characterized in that, include: An initial reservoir mesh model is obtained, which is obtained by dividing the initial reservoir model into tetrahedral meshes. A preset layer extraction strategy is used to extract sub-layers from the initial reservoir grid model to generate an initial sub-layer model; the initial sub-layer model is composed of multiple sub-layers, and each sub-layer includes multiple triangular mesh elements. A preset mesh processing strategy is used to convert triangular mesh cells in multiple small layers into rectangular mesh cells to generate the target small layer model; A target reservoir grid model is generated based on the grid coordinate system and the target small-level model; all grid types in the target reservoir grid model are hexahedral grids.

2. The method according to claim 1, characterized in that, The step of extracting sub-layers from the initial reservoir grid model using a preset layer extraction strategy to generate an initial sub-layer model includes: The sedimentary patterns corresponding to each stratum are determined using the strata represented in the initial reservoir grid model as the dimension; Based on the sedimentary patterns of each stratum and the corresponding preset layer extraction interval values, small layers of each stratum are extracted along the height direction of the initial reservoir grid model. The initial small-level model is generated based on the extracted small-levels of each stratum.

3. The method according to claim 1 or 2, characterized in that, The step of converting triangular mesh cells in multiple small layers into rectangular mesh cells using a preset mesh processing strategy to generate a target small layer model includes: A preset mesh transformation strategy is used to convert triangular mesh units in multiple small layers into hybrid mesh units; the hybrid mesh units include rectangular mesh units far from the fault region and non-rectangular mesh units close to the fault region; The non-rectangular mesh cells are processed using a preset mesh regularization strategy to generate a target small-level model; the multiple small levels of the target small-level model are all composed of rectangular mesh cells.

4. The method according to claim 3, characterized in that, The step of converting triangular mesh cells in multiple small layers into hybrid mesh cells using a preset mesh transformation strategy includes: The triangular mesh cells of the small layer are divided by a unit rectangular mesh based on the grid coordinate system to obtain a hybrid mesh cell; the unit rectangular mesh of the grid coordinate system is formed by the intersection of multiple horizontal grid lines and multiple horizontal grid lines.

5. The method according to claim 4, characterized in that, The step of using a preset mesh processing strategy to perform mesh regularization processing on the non-rectangular mesh cells to generate the target small-scale model includes: Determine whether the non-rectangular mesh cell covers the center point of the unit rectangular mesh in the mesh coordinate system; In response to the non-rectangular mesh cell covering the center point of the unit rectangular mesh in the mesh coordinate system, the non-rectangular mesh cell is replaced with the unit rectangular mesh to generate the target small-scale model.

6. The method according to claim 5, characterized in that, Also includes: In response to the non-rectangular mesh cell not covering the center point of the unit rectangular mesh of the mesh coordinate system, the non-rectangular mesh cell is deleted.

7. The method according to claim 5 or 6, characterized in that, After replacing the non-rectangular mesh cells with the unit rectangular mesh, the method further includes: Determine whether the vertex coordinates of adjacent sides of adjacent rectangular grid cells are consistent; In response to the inconsistency of vertex coordinates of adjacent sides of adjacent rectangular grid cells, the vertices of adjacent sides of adjacent rectangular grid cells are merged.

8. The method according to any one of claims 1-7, characterized in that, The generation of the target reservoir grid model based on the grid coordinate system and the target small-level model includes: Multiple unit mesh columns are generated based on the mesh coordinate system and the target small-level model, and each unit mesh column includes multiple hexahedral meshes arranged along the extension direction of the unit mesh column. Multiple unit grid columns are merged to generate a target reservoir grid model.

9. The method according to claim 8, characterized in that, The generation of multiple element mesh columns based on the mesh coordinate system and the target small-level model includes: Based on the horizontal grid lines, horizontal grid lines, and height grid lines of the grid coordinate system, multiple initial grid cylinders are generated. The target small-level model is segmented based on the initial mesh pillars to generate multiple unit mesh pillars.

10. The method according to claim 8 or 9, characterized in that, After generating multiple element mesh columns based on the mesh coordinate system and the target small-level model, the process further includes: Determine whether a reverse fault exists within the unit grid column; In response to the absence of reverse faults within the unit grid column, the target sub-layers within the unit grid column are arranged according to the sedimentary sequence of the strata. In response to the presence of a reverse fault within the unit grid column, the target sub-layers within the unit grid column are arranged from high to low according to their height coordinate values.

11. The method according to claim 10, characterized in that, Also includes: Determine whether there are any intersecting target sub-layers within the unit grid column; In response to determining that there are intersecting target sub-layers within the unit grid column, the target sub-layers within the unit grid column are rearranged from high to low according to their height coordinate values.

12. The method according to any one of claims 8-11, characterized in that, After generating multiple element mesh columns based on the mesh coordinate system and the target small-level model, the process further includes: Add a three-dimensional coordinate index to the hexahedral mesh within the unit mesh column.

13. The method according to claim 12, characterized in that, Also includes: Determine whether there are duplicate 3D coordinate indices for the hexahedral mesh within the unit mesh column; In response to the presence of index repetition, adjust the coordinate index value of any dimension of the hexahedral mesh with the repetitive index.

14. A mesh generation device for a reservoir model, characterized in that, include: The acquisition module is used to acquire the initial reservoir mesh model, which is obtained by dividing the initial reservoir model into tetrahedral meshes; An extraction module is used to extract sub-layers from the initial reservoir grid model using a preset layer extraction strategy to generate an initial sub-layer model; the initial sub-layer model is composed of multiple sub-layers, and each sub-layer includes multiple triangular mesh elements; The conversion module is used to convert triangular mesh cells in multiple small layers into rectangular mesh cells using a preset mesh processing strategy to generate a target small layer model. The model generation module is used to generate a target reservoir grid model based on the grid coordinate system and the target small-level model; the grid type in the target reservoir grid model is a hexahedral grid.

15. An electronic device, characterized in that, include: A processor, and a memory communicatively connected to the processor; The memory stores computer-executed instructions; The processor executes computer execution instructions stored in the memory to implement the method as described in any one of claims 1-13.

16. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer-executable instructions, which, when executed by a processor, are used to implement the method as described in any one of claims 1-13.

17. A computer program product, characterized in that, Includes a computer program that, when executed by a processor, implements the method of any one of claims 1-13.

Images