Reservoir structure modeling grid partitioning method and device

Abstract

The application provides a reservoir structure modeling grid division method and device, and relates to the technical field of reservoir geological modeling. The method comprises the following steps: obtaining modeling parameters; the modeling parameters are used for determining geological discontinuum data of a target work area; initial points of the geological discontinuum data are established according to the modeling parameters; extension parameters are received, and regional point distribution of the target work area is determined according to the extension parameters and the initial points; grid division is performed based on the regional point distribution, and a reservoir structure modeling grid division result is obtained. In the embodiment of the application, fractures, faults and wells are considered as initial points of grid division, and geological and development factors can be well characterized. By setting extension parameters for straight wells and fracture geological bodies, grid densification can be naturally realized, and high simulation accuracy can be realized with a small number of grids.

Description

Technical Field

[0001] This invention relates to the field of reservoir geological modeling technology, and in particular to a method and apparatus for reservoir structural modeling mesh generation. Background Technology

[0002] This section is intended to provide background or context for the embodiments of the invention set forth in the claims. The description herein is not an admission that it is prior art simply because it is included in this section.

[0003] Reservoir geological modeling combined with numerical simulation is an effective method to guide the rational development of oil fields. Geological modeling technology generates digital geological bodies by rationally processing data from geophysical exploration, well logging, and geology. The accuracy of its description of geological features determines the accuracy of subsequent production capacity prediction, making it a necessary technology for achieving economical oilfield exploitation.

[0004] Geological models are typically represented as grid-based data volumes, specifically including Cartesian grids, corner grids, triangular prism grids, and perpendicular bisection (PEBI) grids. Cartesian grids are pairwise orthogonal in all three dimensions, ensuring simulation accuracy, but their ability to describe geological volumes is limited, making them suitable only for conceptual modeling. Corner grids are an improvement on Cartesian grids, allowing for some stretching in both the horizontal and vertical directions. While the grid cells remain hexahedral, this comes at the cost of some orthogonality. Triangular prism grids have a triangular basic configuration in the two-dimensional plane and a rectangular configuration in the longitudinal section. While highly flexible, their skewness is often difficult to control, leading to significant errors in differential discretization schemes. Perpendicular bisection (PEBI) grids are polyhedral elements with arbitrary polygonal basic configurations in the two-dimensional plane. The grid is star-shaped and suitable for differential schemes. Furthermore, the centroids of any adjacent grid cells are perpendicular to their corresponding edges, thus ensuring the accuracy of flow scheme calculations and making them suitable for reservoir numerical simulations.

[0005] In vertically bisected mesh modeling, the current mainstream meshing method is to use Delaunay triangulation followed by taking the dual mesh. The specific process can be summarized as follows: 1) Place points in the partitioned region; 2) Perform Delaunay triangulation on the existing point set; 3) Take the dual mesh from the triangular mesh; 4) Smooth the mesh. This process can create K-orthogonal meshes with standard shapes, ensuring the accuracy of flow calculations to a certain extent. However, in numerical simulations, the accuracy depends not only on the mesh shape. On one hand, geological discontinuities such as faults, pinch-outs, and large-scale fractures often have complex morphologies, and their permeability differs significantly from the matrix mesh. Traditional vertically bisected meshing methods often ignore the influence of geological discontinuities. On the other hand, the number of meshes is also a key factor affecting the solution. Too many meshes, i.e., overly dense meshes, while ensuring high solution accuracy, often result in excessively large matrices, placing a heavy burden on the solver. Too few meshes result in overly sparse discrete points, leading to lower solution accuracy. Summary of the Invention

[0006] This invention provides a method and apparatus for mesh generation in reservoir structural modeling. By considering fractures, faults, wells, etc. as the starting points for mesh generation, geological and development factors can be well characterized. By setting extension multiples for geological bodies such as vertical wells and fractures, mesh refinement can be naturally achieved, thereby achieving high simulation accuracy with a smaller number of meshes.

[0007] In a first aspect, embodiments of the present invention provide a reservoir structural modeling mesh generation method, the method comprising:

[0008] Obtain modeling parameters; the modeling parameters are used to determine the geological discontinuity data of the target work area; establish the initial points of the geological discontinuity data according to the modeling parameters; receive extension parameters, and determine the regional layout of the target work area according to the extension parameters and the initial points; perform grid subdivision based on the regional layout to obtain the reservoir structure modeling grid subdivision result.

[0009] Secondly, embodiments of the present invention also provide a reservoir structure modeling mesh generation device, the device comprising:

[0010] The acquisition module is used to acquire modeling parameters; these modeling parameters are used to determine the geological discontinuities data of the target work area.

[0011] The geological module is used to establish the initial points of the geological discontinuity data based on the modeling parameters; the density module is used to receive the extension parameters and determine the regional layout of the target work area based on the extension parameters and the initial points; the meshing module is used to perform meshing based on the regional layout to obtain the reservoir structure modeling meshing results.

[0012] Thirdly, embodiments of the present invention also provide an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the above-described reservoir structure modeling mesh generation method.

[0013] Fourthly, embodiments of the present invention also provide a computer-readable storage medium storing a computer program for performing the above-described reservoir structure modeling mesh generation method.

[0014] The embodiments of this invention bring the following beneficial effects: These embodiments provide a reservoir structure modeling grid partitioning scheme, which includes: acquiring modeling parameters; using the modeling parameters to determine geological discontinuity data for the target work area; establishing initial points for the geological discontinuity data based on the modeling parameters; receiving extension parameters; determining the regional distribution of points for the target work area based on the extension parameters and the initial points; and performing grid partitioning based on the regional distribution to obtain the reservoir structure modeling grid partitioning result. These embodiments of the invention can effectively characterize geological and development factors by considering fractures, faults, wells, etc., as initial points for grid partitioning; and by setting extension parameters for geological bodies such as vertical wells and fractures, grid refinement can be naturally achieved, thereby achieving high simulation accuracy with a smaller number of grids.

[0015] Other features and advantages of the invention will be set forth in the description which follows, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention are realized and obtained in accordance with the structures particularly pointed out in the description, claims and drawings.

[0016] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, preferred embodiments are described below in detail with reference to the accompanying drawings. Attached Figure Description

[0017] To more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the specific embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0018] Figure 1 A flowchart of the reservoir structure modeling mesh generation method provided in this embodiment of the invention;

[0019] Figure 2 This is a well unit layout diagram provided in an embodiment of the present invention;

[0020] Figure 3This is a large-scale fracture element layout diagram for a horizontal well provided in an embodiment of the present invention;

[0021] Figure 4 This is a fault element layout diagram provided in an embodiment of the present invention;

[0022] Figure 5 A diagram illustrating the frontier advancement steps provided for embodiments of the present invention;

[0023] Figure 6 A flowchart of mesh generation provided for an embodiment of the present invention;

[0024] Figure 7 This is a schematic diagram of the electronic device system composition structure provided in an embodiment of the present invention;

[0025] Figure 8 A structural block diagram of a reservoir structure modeling mesh generation device provided in an embodiment of the present invention;

[0026] Figure 9 A structural block diagram of another reservoir structure modeling mesh generation device provided in an embodiment of the present invention;

[0027] Figure 10 A structural block diagram of a geological module provided in an embodiment of the present invention;

[0028] Figure 11 A grid diagram of an embodiment of the present invention is provided;

[0029] Figure 12 A mesh diagram of another embodiment provided for the present invention. Detailed Implementation

[0030] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0031] Based on this, in order to better improve simulation accuracy during the vertical bisecting grid generation process, this invention provides a grid generation algorithm that can consider complex geological bodies and generate a density-controllable grid. The embodiments of this invention provide a reservoir structure modeling grid generation method and apparatus. This method is based on a vertical bisecting grid framework to construct a reservoir grid system that is easy to numerically simulate, belonging to a density-controllable grid that considers geological bodies. Using reservoir geological parameters measured in the field, a grid model of the area can be easily established; by considering fractures, faults, wells, etc., as the starting points for grid generation, geological and development factors can be well characterized; by setting extension factors for geological bodies such as vertical wells and fractures, grid refinement can be naturally achieved, thereby achieving high simulation accuracy with a smaller number of grids.

[0032] To facilitate understanding of this embodiment, a detailed description of a reservoir structural modeling mesh generation method disclosed in this embodiment of the invention will be provided first.

[0033] This invention provides a method for mesh generation in reservoir structural modeling, see [link to relevant documentation]. Figure 1 The flowchart shown is a method for mesh generation in reservoir structural modeling. The method includes the following steps:

[0034] Step S102: Obtain modeling parameters.

[0035] In this embodiment of the invention, modeling parameters are used to determine the geological discontinuities data of the target work area. Geological discontinuities include special or complex morphologies such as faults, pinch-outs, and large-scale fractures. Geological discontinuity data is used to describe the location, morphology, etc., of various types of geological discontinuities.

[0036] Basic modeling parameters may include: simulation area boundary, well location coordinates, rectangular fracture coordinates, contour map, and fault hanging wall and footwall strike lines, which are used as basic input parameters for modeling.

[0037] Step S104: Establish the initial points for geological discontinuity data based on the modeling parameters.

[0038] In this embodiment of the invention, initial points are established for various geological discontinuities. These initial points are used for subsequent mesh creation and mesh generation. The initial points can be used to describe the approximate location and shape of the geological discontinuities.

[0039] Step S106: Receive extension parameters and determine the area layout of the target work area based on the extension parameters and the initial point.

[0040] In this embodiment of the invention, the extension parameter can be a pre-set parameter used to define the mesh density. After determining the initial point of the discontinuity data, points are placed in the target work area excluding the discontinuities to obtain the area point layout of the target work area.

[0041] Step S108: Based on the regional point layout, perform grid subdivision to obtain the reservoir structure modeling grid subdivision result.

[0042] In this embodiment of the invention, after obtaining the initial point and the area distribution, the mesh is divided using a meshing method to obtain the meshing result. Based on the meshing result, the final mesh can be obtained.

[0043] This invention provides a reservoir structure modeling mesh generation scheme, which includes: acquiring modeling parameters; using the modeling parameters to determine geological discontinuity data of the target work area; establishing initial points for the geological discontinuity data based on the modeling parameters; receiving extension parameters; determining the regional distribution of points in the target work area based on the extension parameters and the initial points; and performing mesh generation based on the regional distribution to obtain the reservoir structure modeling mesh generation result. This invention can effectively characterize geological and development factors by considering fractures, faults, wells, etc., as initial points for mesh generation; and by setting extension parameters for geological bodies such as vertical wells and fractures, mesh refinement can be naturally achieved, thus achieving high simulation accuracy with a smaller number of meshes.

[0044] In one embodiment, the following steps may also be performed before obtaining the modeling parameters:

[0045] Obtain the geological parameters of the target work area; generate modeling parameters based on the geological parameters.

[0046] In this embodiment of the invention, the geological parameters of the target work area include at least one or more of the following: geophysical data, well logging data, geological mapping data, seismic data, drilling data, and development plan data. Basic modeling parameters are obtained using data from on-site geophysical exploration, well logging, geological mapping, seismic data, drilling data, and development plan data.

[0047] In one embodiment, the geological discontinuity data includes at least one or more of the following: vertical well unit data, horizontal well fracture unit data, and fault unit data; establishing the initial points for the geological discontinuity data based on the modeling parameters can be performed according to the following steps:

[0048] The initial points for vertical well unit data are established based on the modeling parameters; the initial points for horizontal well fracture unit data are established based on the modeling parameters; and the initial points for fault unit data are established based on the modeling parameters.

[0049] In this embodiment of the invention, the horizontal well fracture unit data may include any one or a combination of horizontal well data and fracture data. Initial points are constructed for the vertical well unit data, horizontal well fracture unit data, and fault unit data according to their respective characteristics to preserve specific geological and development features.

[0050] In one embodiment, establishing the initial points for vertical well unit data based on modeling parameters can be performed using the following steps:

[0051] The initial point for the vertical well unit data is established based on the position and shape parameters of the vertical well unit.

[0052] In an embodiment of the present invention, see Figure 2 The well unit layout diagram shown illustrates that, for vertical wells, the location of the vertical well can be determined based on the modeling parameters. Based on the location of the vertical well, the following can be determined: Figure 2 The black dots represent the initial points of the vertical well unit data, determined by multiple vertices along the well edge in the diagram.

[0053] In one embodiment, establishing the initial points for horizontal well fracture element data based on modeling parameters can be performed according to the following steps:

[0054] The initial points of the horizontal well fracture elements are established based on their location and shape parameters.

[0055] In this embodiment of the invention, the location and shape parameters of horizontal well fracture elements can be determined based on modeling parameters. Figure 3 The horizontal lines in the diagram represent horizontal wells. Multiple points are marked on the horizontal wells. Based on the marked points on the horizontal wells, symmetrical auxiliary points are determined on both sides of the horizontal wells. These points on the horizontal wells and the auxiliary points are used as the initial points of the horizontal wells. The method for determining the initial points of fractures is similar to that for horizontal wells.

[0056] It should be noted that, as Figure 3 In the process, initial points that are close to horizontal wells and fractures can be merged during grid generation because of their proximity. For multiple initial points at the same location or within a certain distance range, only one needs to be used.

[0057] In one embodiment, establishing the initial points of fault element data based on modeling parameters can be performed according to the following steps:

[0058] The initial points of the fault element data are established based on the location and shape parameters of the fault element.

[0059] In this embodiment of the invention, the location and shape parameters of the fault element can be determined based on the modeling parameters. See also Figure 4 Points are symmetrically placed on both sides of the fault line. If the fault is an intersecting fault, points are placed radially at the intersection with a certain radius, with two circles of points. Points are also placed in mirror image with a certain width on both sides of the intersection of each circle with the fault line to obtain the initial points of the fault unit data.

[0060] In one embodiment, determining the area layout of the target work area based on the extension parameters and the initial point can be performed according to the following steps:

[0061] Determine discrete edge group data based on the initial point; determine pre-insertion points based on discrete edge group data and extension parameters; generate guard circle parameters for pre-insertion points based on extension parameters; determine the region layout based on guard circle parameters and pre-insertion points.

[0062] In this embodiment of the invention, the extension parameter can be a pre-given constant field data, which is used to determine the subsequent grid generation density and to determine approximately how many grids can be drawn in a region. Too many grids will result in a large amount of computation, while too few grids will affect the accuracy of the calculation. The grid generation density can be controlled by the extension parameter.

[0063] Based on the extension parameters, a smooth lattice can be generated in the remaining area of ​​the target work area, excluding geological discontinuities. It should be noted that the extension parameters can be used to determine the point density within the lattice. Then, combined with the initial points, discrete edge group data is determined.

[0064] Discrete edge group data includes lines connecting the outer points of geological bodies such as wells, faults, and fractures. See also Figure 5 See Figure 5 1. First, select two points on one side of the line and draw the perpendicular bisector of this line segment. Determine the length of the perpendicular bisector using the extension parameters. Draw a circle (protective circle) at the endpoint (insertion point) of the perpendicular bisector. See [link / reference]. Figure 5 In step 2, if there are no other points inside the circle, then connecting them forms a triangle. See [reference needed]. Figure 5 If there are other points, 3 in the middle. Figure 5 4. Connect the other points and the endpoints of the line segment, see [reference]. Figure 5 5 in the middle, then continue on both the left and right sides... see Figure 5 16. The vertices of the resulting multiple triangles will be used as the region points.

[0065] In one embodiment, mesh generation based on region-based point distribution can be performed according to the following steps:

[0066] Delaunay triangulation is performed based on regional point distribution; the vertices of each triangle are traversed, and the centroids of the triangles sharing a common vertex are connected in a clockwise direction to obtain the mesh subdivision result for reservoir structural modeling.

[0067] In this embodiment of the invention, after obtaining the reservoir structure modeling mesh subdivision result, interpolation is performed based on the contour data of each layer to complete the structure modeling.

[0068] The execution process of this method is illustrated below with a specific embodiment. A two-dimensional mesh generation method for reservoir structural modeling follows this process:

[0069] 1) Obtain geological parameters:

[0070] Using data from geophysical exploration, well logging, geological mapping, seismic surveys, drilling, and development plans conducted on-site, basic modeling parameters are obtained, including: simulation area boundaries, well location coordinates, rectangular fracture coordinates, contour maps, and fault strike lines on the hanging wall and footwall, which are used as basic input parameters for modeling.

[0071] 2) Basic unit layout:

[0072] To preserve specific geological and development characteristics, the embodiments of this invention establish the following basic unit configurations, which form the skeleton of the entire grid structure, as follows:

[0073] ① Vertical well unit

[0074] The vertical well unit is an equilateral polygon, and the initial point layout is as follows: Figure 2 As shown by the black dot, the line connecting the points on the outer side is the initial advance edge of the well unit, the dashed line is the corresponding Delaunay triangular mesh, and the solid line mesh is the final well mesh.

[0075] ② Horizontal wells and large-scale fracture units

[0076] The final form of horizontal wells and large-scale fracture elements is a combination of rectangular elements. Based on the properties of Delaunay triangles, the grid layout of the horizontal well is designed as follows: Figure 3 As shown in the figure, the middle point is the Delaunay triangle skeleton point. Points are symmetrically distributed on both sides of the skeleton point along the perpendicular line of the skeleton point. The dual mesh is the horizontal well and large-scale fracture element to be obtained.

[0077] ③ Fault unit

[0078] In this invention, a dimensionality reduction method is used to handle faults. Similar to horizontal well / large-scale fracture elements, the fault element processing is also based on the properties of the Delaunay triangle, and the fault point layout format is designed as follows ( Figure 4 ): Points are symmetrically arranged on both sides of the fault line. If the fault is an intersecting fault, points are arranged radially at the intersection with a certain radius. The number of circles of points is two. On both sides of the intersection of each circle with the fault line, points are arranged in a mirror image with a certain width.

[0079] 3) Generate density-constrained field

[0080] Based on actual needs, input appropriate density parameters as edge extension ratios (indirectly reflecting grid density information). This density information applies to the outermost edge of the basic configuration of vertical wells and fractures.

[0081] 4) Use the discrete edge-front advancement method to place points in the remaining region.

[0082] This invention defines a discrete edge-front advancement method, which can generate a smooth Delaunay lattice based on density information and has the property of absolute convergence. Figure 5 The deployment process of the forward propulsion method is demonstrated, and the specific steps of this method are explained below:

[0083] ① Establish discrete edge groups by connecting the outer points of geological bodies such as wells, faults, and fractures to form the initial discrete edge groups;

[0084] ② Select the smallest edge, and according to the pre-set extension parameters for the geological body, extend it a certain distance from the midpoint along the perpendicular bisector to form a pre-insertion point;

[0085] ③ Based on the extension distance in step ②, form a protective circle around the pre-insertion point, within which no other point can be inserted;

[0086] ④ Determine the validity of the pre-insertion point, specifically: a) If there are other points within the protection radius of the point, then form a list of all points within the protection radius and merge the point with the nearest point in the list; b) If the extension line of the point intersects the edge of the discrete edge group of any stage, then the point is invalid; c) If the point is outside the boundary, then the point is invalid.

[0087] ⑤ If the point is valid, then consider its extended edge as the parent edge, and connect the point to the two endpoints of the parent edge to form two new discrete edges;

[0088] ⑥ Regardless of the validity of a point, once a discrete edge (mother edge) is selected, it cannot be reused and must be removed from the edge group;

[0089] ⑦ Repeat steps ② through ⑥ until the edge group is empty.

[0090] 5) Perform Delaunay triangulation on the entire point set;

[0091] 6) Traverse each triangle vertex and connect the centroids of the triangles sharing a common vertex in a clockwise direction. The resulting grid is the final grid.

[0092] 7) Interpolate the contour data of each layer to complete the structural modeling.

[0093] See Figure 6 The mesh generation flowchart shown below, in one embodiment, can be implemented according to the following steps:

[0094] A. Obtain geological information, fracture diagnosis information, and well location information for the area, and import them into the modeling program:

[0095] A1. Obtain three-dimensional geological structural information of the area, including three-dimensional geological coordinates, fault strike, geological stratification information, and contour maps of each stratum;

[0096] A2. Obtain the crack diagnosis information for this area, i.e., the equivalent rectangular crack coordinates;

[0097] A3. Obtain well information for the area, including target coordinates, well diameter, and well trajectory for vertical and horizontal wells;

[0098] B. After importing the partitioning program, begin the calculations in the following order:

[0099] B1. Based on well data, fault data, fracture data, and boundary data, establish initial points;

[0100] B2. Enter the extension length value for each source item (i.e., well, fracture, etc.);

[0101] B3. Conduct forward deployment and complete regional deployment;

[0102] B4. Perform Delaunay triangulation and obtain the dual mesh.

[0103] C. Export the results and perform visualization processing.

[0104] Assuming a rectangular region containing three intersecting large-scale cracks, mesh this region. After performing the above steps to mesh the individual cells, the result might be as follows: Figure 11 A mesh diagram of one embodiment is shown.

[0105] Assuming a construction area with a complex boundary, containing nine vertical wells (corresponding to WELL1, WELL2, ..., WELL9 in the diagram), and four faults with significant layer misalignment between the fault layers, meshing this area might yield the following results: Figure 12 Another embodiment of the mesh diagram is shown.

[0106] This invention provides a method and apparatus for mesh generation in reservoir structural modeling. When performing mesh generation, this invention does not require manually inserted mesh cells beforehand, and it has the advantages of high speed and high reliability compared to other polygon mesh generation algorithms. Furthermore, it preserves the orthogonality of the mesh while explicitly representing the geological body.

[0107] This invention also provides a reservoir structure modeling mesh generation device, as described in the following embodiments. Since the principle behind this device's problem-solving is similar to that of the reservoir structure modeling mesh generation method, its implementation can be found in the implementation of the reservoir structure modeling mesh generation method; repeated details will not be elaborated further. See [link to documentation]. Figure 8 The diagram shows a structural block diagram of a reservoir structure modeling mesh generation device, which includes:

[0108] The acquisition module 81 is used to acquire modeling parameters; the modeling parameters are used to determine the geological discontinuity data of the target work area; the geology module 82 is used to establish the initial points of the geological discontinuity data according to the modeling parameters; the density module 83 is used to receive extension parameters and determine the regional distribution of the target work area according to the extension parameters and the initial points; the meshing module 84 is used to perform meshing based on the initial points and the regional distribution to obtain the reservoir structure modeling meshing results.

[0109] In one embodiment, see Figure 9 The diagram shows another reservoir structure modeling mesh generation device, which also includes a parameter module 85 for:

[0110] Obtain the geological parameters of the target work area; the geological parameters include at least one or more of the following: geophysical data, well logging data, geological mapping data, seismic data, drilling data, and development plan data; generate modeling parameters based on the geological parameters.

[0111] In one embodiment, see Figure 10 The geological module structure diagram shown indicates that the geological discontinuity data includes at least one or more of the following: vertical well unit data, horizontal well fracture unit data, and fault unit data; the geological module includes:

[0112] The first initial point unit 821 is used to establish the initial point of the vertical well unit data according to the modeling parameters; the second initial point unit 822 is used to establish the initial point of the horizontal well fracture unit data according to the modeling parameters; and the third initial point unit 823 is used to establish the initial point of the fault unit data according to the modeling parameters.

[0113] In one embodiment, the first initial point unit is specifically used to: establish the initial point of the vertical well unit data based on the position parameters and shape parameters of the vertical well unit.

[0114] In one embodiment, the second initial point unit is specifically used to: establish the initial point of the horizontal well fracture unit based on the position and shape parameters of the horizontal well fracture unit.

[0115] In one embodiment, the third initial point unit is specifically used to: establish the initial point of the fault unit data based on the position parameters and shape parameters of the fault unit.

[0116] In one embodiment, the density module is specifically used for: determining discrete edge group data based on the initial point; determining pre-insertion points based on the discrete edge group data and extension parameters; generating guard circle parameters for the pre-insertion points based on the extension parameters; and determining the region layout based on the guard circle parameters and the pre-insertion points.

[0117] In one embodiment, the subdivision module is specifically used for: performing Delaunay triangulation based on regional point distribution; traversing each triangle vertex, and connecting the centroids of the triangles sharing a common vertex in a clockwise direction to obtain the reservoir structure modeling mesh subdivision result.

[0118] Based on the same inventive concept, this invention also provides an embodiment of an electronic device for implementing all or part of the above-described reservoir structure modeling mesh generation method. This electronic device specifically includes the following:

[0119] The device comprises a processor, memory, a communications interface, and a bus; wherein the processor, memory, and communications interface communicate with each other via the bus; the communications interface is used to realize information transmission between related devices; the electronic device can be a desktop computer, tablet computer, or mobile terminal, etc., and this embodiment is not limited to these. In this embodiment, the electronic device can be implemented with reference to the embodiments for implementing the above-mentioned reservoir structure modeling mesh generation method and the embodiments for implementing the above-mentioned reservoir structure modeling mesh generation device, the contents of which are incorporated herein by reference, and repeated details will not be described again.

[0120] Figure 7 This is a schematic diagram of the system composition structure of an electronic device provided in an embodiment of the present invention. Figure 7 As shown, the electronic device 70 may include a processor 701 and a memory 702; the memory 702 is coupled to the processor 701. It is worth noting that... Figure 7 This is an example; other types of structures can also be used to supplement or replace this structure to achieve telecommunications functions or other functions.

[0121] In one embodiment, the functionality of the reservoir structure modeling mesh generation method can be integrated into the processor 701. The processor 701 can be configured to perform the following controls:

[0122] Obtain modeling parameters; modeling parameters are used to determine the geological discontinuity data of the target work area; establish the initial points of the geological discontinuity data based on the modeling parameters; receive extension parameters, and determine the regional layout of the target work area based on the extension parameters and the initial points; perform grid subdivision based on the regional layout to obtain the reservoir structure modeling grid subdivision results.

[0123] As can be seen from the above, the electronic device provided in the embodiments of the present invention can effectively characterize geological and development factors by considering fractures, faults, wells, etc. as the initial points for grid partitioning; by setting extension parameters for geological bodies such as vertical wells and fractures, grid densification can be naturally achieved, thereby achieving high simulation accuracy with a smaller number of grids.

[0124] In another embodiment, the reservoir structure modeling mesh generation device can be configured separately from the processor 701. For example, the reservoir structure modeling mesh generation device can be configured as a chip connected to the processor 701, and the function of the reservoir structure modeling mesh generation method can be realized through the control of the processor.

[0125] like Figure 7 As shown, the electronic device 70 may further include: a communication module 703, an input unit 704, an audio processing unit 705, a display 706, and a power supply 707. It is worth noting that the electronic device 70 does not necessarily need to include these components. Figure 7 All components shown; in addition, the electronic device 70 may also include Figure 7 For components not shown, please refer to existing technologies.

[0126] like Figure 7 As shown, processor 701, sometimes also referred to as controller or operation control, may include a microprocessor or other processor device and / or logic device, which receives input and controls the operation of various components of electronic device 70.

[0127] The memory 702 may be, for example, one or more of a cache, flash memory, hard drive, removable media, volatile memory, non-volatile memory, or other suitable device. It may store the aforementioned failure-related information, and also store a program for executing that information. The processor 701 may execute the program stored in the memory 702 to perform information storage or processing, etc.

[0128] Input unit 704 provides input to processor 701. Input unit 704 may be, for example, a keypad or touch input device. Power supply 707 provides power to electronic device 70. Display 706 displays images and text. Display may be, for example, an LCD display, but is not limited to this.

[0129] The memory 702 can be a solid-state memory, such as a read-only memory (ROM), random access memory (RAM), a SIM card, etc. It can also be a memory that retains information even when power is off, can be selectively erased, and contains more data; examples of this type of memory are sometimes referred to as EPROMs. The memory 702 can also be some other type of device. The memory 702 includes a buffer memory 7021 (sometimes called a buffer). The memory 702 may include an application / function storage unit 7022 for storing application programs and function programs or processes for executing the operation of the electronic device 70 via the processor 701.

[0130] The memory 702 may also include a data storage unit 7023 for storing data, such as contacts, digital data, pictures, sounds, and / or any other data used by the electronic device. The driver storage unit 7024 of the memory 702 may include various drivers for the electronic device's communication functions and / or for performing other functions of the electronic device (such as messaging applications, address book applications, etc.).

[0131] The communication module 703 is a transmitter / receiver that transmits and receives signals via the antenna 708. The communication module (transmitter / receiver) 703 is coupled to the processor 701 to provide input signals and receive output signals, which is the same as in a conventional mobile communication terminal.

[0132] Based on different communication technologies, multiple communication modules 703 can be configured in the same electronic device, such as cellular network modules, Bluetooth modules, and / or wireless LAN modules. The communication module (transmitter / receiver) 703 is also coupled to a speaker 709 and a microphone 710 via an audio processing unit 705 to provide audio output via the speaker 709 and receive audio input from the microphone 710, thereby realizing typical telecommunications functions. The audio processing unit 705 may include any suitable buffer, decoder, amplifier, etc. Additionally, the audio processing unit 705 is also coupled to a processor 701, enabling on-device recording via the microphone 710 and on-device playback of stored audio via the speaker 709.

[0133] In embodiments of the present invention, a computer-readable storage medium is also provided for implementing all steps of the reservoir structure modeling mesh generation method in the above embodiments. The computer-readable storage medium stores a computer program, which, when executed by a processor, implements all steps of the reservoir structure modeling mesh generation method in the above embodiments. For example, when the processor executes the computer program, it implements the following steps:

[0134] Obtain modeling parameters; modeling parameters are used to determine the geological discontinuity data of the target work area; establish the initial points of the geological discontinuity data based on the modeling parameters; receive extension parameters, and determine the regional layout of the target work area based on the extension parameters and the initial points; perform grid subdivision based on the regional layout to obtain the reservoir structure modeling grid subdivision results.

[0135] As can be seen from the above, the computer-readable storage medium provided in the embodiments of the present invention can effectively characterize geological and development factors by considering fractures, faults, wells, etc. as the initial points for grid partitioning; by setting extension parameters for geological bodies such as vertical wells and fractures, grid densification can be naturally achieved, thereby achieving higher simulation accuracy with a smaller number of grids.

[0136] While this invention provides the method operation steps as described in the embodiments or flowcharts, more or fewer operation steps may be included based on conventional or non-inventive labor. The order of steps listed in the embodiments is merely one possible execution order among many and does not represent the only possible execution order. In actual device or client product execution, the methods shown in the embodiments or drawings can be executed sequentially or in parallel (e.g., in a parallel processor or multi-threaded processing environment).

[0137] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0138] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0139] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0140] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0141] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between embodiments can be referred to interchangeably. Each embodiment focuses on describing the differences from other embodiments. In particular, the system embodiments are basically similar to the method embodiments, so the description is relatively simple; relevant parts can be referred to the descriptions in the method embodiments.

[0142] In this document, relational terms such as "first" and "second" are used merely to distinguish one entity or operation from another, without necessarily requiring or implying any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. The terms "upper," "lower," etc., indicating orientation or positional relationships based on the orientation or positional relationships shown in the accompanying drawings, are used only for the convenience of describing the invention and for simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore should not be construed as limiting the invention.

[0143] Unless otherwise explicitly specified and limited, the terms "installation," "connection," and "linking" should be interpreted broadly. For example, they can refer to fixed connections, detachable connections, or integral connections; they can refer to mechanical connections or electrical connections; they can refer to direct connections or indirect connections through an intermediate medium; and they can refer to the internal connection between two components. Those skilled in the art can understand the specific meaning of the above terms in this invention according to the specific circumstances.

[0144] It should be noted that, unless otherwise specified, the embodiments and features described in the present invention can be combined with each other. The present invention is not limited to any single aspect, nor to any single embodiment, nor to any combination and / or substitution of these aspects and / or embodiments. Each aspect and / or embodiment of the present invention can be used alone, or in combination with one or more other aspects and / or other embodiments.

[0145] Finally, it should be noted that the above-described embodiments are merely specific implementations of the present invention, used to illustrate the technical solutions of the present invention, and not to limit it. The scope of protection of the present invention is not limited thereto. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that any person skilled in the art can still modify or easily conceive of changes to the technical solutions described in the foregoing embodiments within the technical scope disclosed in the present invention, or make equivalent substitutions for some of the technical features; and these modifications, changes, or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be covered within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A method for mesh generation in reservoir structural modeling, characterized in that, include: Obtain modeling parameters; The modeling parameters are used to determine the geological discontinuities in the target work area; The initial points of the geological discontinuity data are established based on the modeling parameters; Receive extension parameters, which are used to define the density of subsequent grid generation. Determine the regional layout of the target work area based on the extension parameters and the initial point, which includes: determining discrete edge group data based on the initial point, wherein the discrete edge group data includes the connection lines of the outer points of geological bodies such as wells, faults, and fractures. First, select two points on one side of the line connecting the outer points, and draw the perpendicular bisector of the line segment connecting the two points. The length of the perpendicular bisector is determined by the extension parameter. Draw a protective circle at the endpoint of the perpendicular bisector, which is the pre-insertion point. If there are no other points in the protective circle, a triangle is formed. If there are other points, connect the other points to the endpoints of the line segment. Then, continue to determine the discrete edge group data on the left and right sides of the line connecting the outer points. The vertices of the multiple triangles obtained in the end are used as the region points. If the extension line of the pre-insertion point intersects the edge of the discrete edge group in any stage, or if the pre-insertion point is outside the boundary, then the pre-insertion point is invalid. Based on the area distribution, a Delaunay triangular mesh is generated. By traversing the vertices of each triangle and connecting the centroids of the triangles sharing a common vertex in a clockwise direction, the mesh generation result for reservoir structural modeling is obtained.

2. The method according to claim 1, characterized in that, Before obtaining the modeling parameters, the following is also included: Obtain the geological parameters of the target work area; the geological parameters include at least one or more of the following: geophysical data, well logging data, geological mapping data, seismic data, drilling data, and development plan data; Modeling parameters are generated based on the geological parameters.

3. The method according to claim 1, characterized in that, The geological discontinuity data includes at least one or more of the following: vertical well unit data, horizontal well fracture unit data, and fault unit data; The initial points for establishing the geological discontinuity data based on the modeling parameters include: The initial points for the vertical well unit data are established based on the modeling parameters. The initial points for the horizontal well fracture element data are established based on the modeling parameters. The initial points of the fault unit data are established based on the modeling parameters.

4. The method according to claim 3, characterized in that, The initial points for the vertical well unit data are established based on the modeling parameters, including: The initial point of the vertical well unit data is established based on the position and shape parameters of the vertical well unit.

5. The method according to claim 3, characterized in that, The initial points for establishing the horizontal well fracture element data based on the modeling parameters include: The initial point of the horizontal well fracture element is established based on its position and shape parameters.

6. The method according to claim 3, characterized in that, Establishing the initial points of the fault unit data based on the modeling parameters includes: The initial point of the fault element data is established based on the location and shape parameters of the fault element.

7. A reservoir structure modeling mesh generation device, characterized in that, include: The acquisition module is used to obtain modeling parameters; The modeling parameters are used to determine the geological discontinuities in the target work area; The geological module is used to establish the initial points of the geological discontinuity data based on the modeling parameters; The density module receives extension parameters, which define the density of subsequent mesh generation. Based on these extension parameters and the initial point, it determines the regional point layout of the target work area. This includes: determining discrete edge group data based on the initial point, where the discrete edge group data includes lines connecting the outer edges of geological bodies such as wells, faults, and fractures; first, selecting two points on one side of the outer edge line, drawing a perpendicular bisector connecting these two points, determining the length of the perpendicular bisector based on the extension parameters, and drawing a protective circle at the endpoint of the perpendicular bisector used as the pre-insertion point. If there are no other points within the protective circle, a triangle is formed. If there are other points, they are connected to the endpoints of the line segment. Then, the steps following the determination of the discrete edge group data continue on both sides of the outer edge line, using the vertices of the resulting triangles as the regional point layout. If the extension line of the pre-insertion point intersects the edge of any stage of the discrete edge group, or if the pre-insertion point is outside the boundary, then the pre-insertion point is invalid. The meshing module is used to perform Delaunay triangular meshing based on the points in the region. It traverses each triangle vertex and takes the centroid of the triangles with common vertices in a clockwise direction to obtain the meshing result of reservoir structure modeling.

8. The apparatus according to claim 7, characterized in that, It also includes a parameter module for: Obtain the geological parameters of the target work area; the geological parameters include at least one or more of the following: geophysical data, well logging data, geological mapping data, seismic data, drilling data, and development plan data; Modeling parameters are generated based on the geological parameters.

9. The apparatus according to claim 7, characterized in that, The geological discontinuity data includes at least one or more of the following: vertical well unit data, horizontal well fracture unit data, and fault unit data; the geological module includes: The first initial point unit is used to establish the initial point of the vertical well unit data according to the modeling parameters; The second initial point unit is used to establish the initial point of the horizontal well fracture unit data according to the modeling parameters; The third initial point unit is used to establish the initial point of the fault unit data based on the modeling parameters.

10. The apparatus according to claim 9, characterized in that, The first initial point unit is specifically used for: The initial point of the vertical well unit data is established based on the position and shape parameters of the vertical well unit.

11. The apparatus according to claim 9, characterized in that, The second initial point unit is specifically used for: The initial point of the horizontal well fracture element is established based on its position and shape parameters.

12. The apparatus according to claim 9, characterized in that, The third initial point unit is specifically used for: The initial point of the fault element data is established based on the location and shape parameters of the fault element.

13. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the reservoir structure modeling mesh generation method according to any one of claims 1 to 6.

14. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that executes the reservoir structure modeling mesh generation method according to any one of claims 1 to 6.

Images