A Basin Numerical Simulation Method and System

Abstract

The present invention discloses a basin numerical simulation method and system. The method includes: obtaining parameters of each geological horizon of a target basin; in response to a corresponding instruction, taking the corresponding geological horizon as a processing object, executing a first rectangular different geological event range increasing strategy, or executing a second rectangular different geological event range increasing strategy, or executing a first polygonal different geological event range increasing strategy, or executing a second polygonal different geological event range increasing strategy, or executing a corresponding different geological event range increasing strategy and making corresponding modifications to the depth values of the corresponding different geological event ranges on the geological horizon to obtain a geological horizon to be numerically simulated; performing block numerical simulation calculations on the geological horizon. The system: when a processor executes a computer program stored in a memory, the above method is implemented. According to the present invention, the problem that the existing basin simulation software has an unreasonable processing method for geological horizons with different geological events can be solved.

Description

Technical Field

[0001] The present invention belongs to the field of basin simulation, and more particularly, relates to a method and system for numerical simulation of basins. Background Art

[0002] Basin simulation refers to quantitatively simulating the formation and evolution of hydrocarbon-bearing basins, the generation, migration, and accumulation of hydrocarbons by a computer in terms of time and space based on physical and chemical geological mechanisms, so as to reveal the essence of the oil and gas laws in the basin.

[0003] The technology of numerical simulation of basins has developed extremely rapidly and made rapid progress in the past decade or so. The wide application of basin simulation technology has enabled basin analysis to develop towards quantification, dynamization, and automated mapping. Nowadays, basin simulation technology is not only a form of expression of geological processes, but also an indispensable means for studying various dynamic parameters. Therefore, it has received extensive attention from basin analysis experts and petroleum geologists.

[0004] As a long-term technical problem, basin simulation faces difficulties and challenges in many aspects such as the accurate solution of simulation algorithms, the consideration of geological factors such as faults and diagenesis, the restoration of intermittent mutation processes of oil and gas migration, and the reconstruction of paleo-hydrodynamic processes. Strengthening the application of three-dimensional geological attribute modeling and structural modeling technologies and implementing interactive simulation at key stages guided by the oil and gas migration and accumulation laws that conform to geological laws will be the main direction for the future development of basin simulation technology. In addition, integrating the latest progress and new technologies of disciplines such as petroleum geology, mathematical geology, and computer science into basin simulation will also greatly promote the progress of this technology.

[0005] In terms of the industrial application of software, there are mainly three basin simulation software products active in the international commercial software market currently: PetroMod of the German Institute of Organic Geochemistry (IES), which consists of three relatively independent systems, namely the cross-sectional two-dimensional oil and gas system analysis software PetroFlow, the planar two-dimensional oil and gas system analysis software Finesse, and the sedimentation analysis software Sedpak; the TemisPack (two-dimensional), Genex (one-dimensional), and Temis 3D (three-dimensional) series software of the French Institute of Petroleum (IFP); BasinMod of the American Platte River Associates (PRA), whose main products are BasinMod 1-D and BasinMod 2-D. The above-mentioned software is comprehensive in content, advanced in technology, and highly commercialized, and has its own unique features in many aspects of solving basin analysis problems.

[0006] The research and development of domestic basin simulation software reached its peak around the 1990s, and more than 10 related software were developed, such as BASIMS (Research Institute of Petroleum Exploration and Development, PetroChina), PRES (CNOOC Research Center), BIAS (former Earth Software Company), GEMDASS, and PASS. In terms of the overall level of commercialized basin simulation software, there is still a gap between domestic and foreign, and continuous efforts are needed to tackle key problems.

[0007] China's basin simulation technology was developed based on tracking Western technologies in the early 1980s, and its development process is generally similar to that in the West. In the first 10 years, one-dimensional models were mainly used, focusing on studying the three histories of the basin, namely geological history, thermal history, and hydrocarbon generation history, and most were in the experimental application stage; in the following 10 years, two-dimensional models were mainly used, focusing on studying the hydrocarbon expulsion history and migration and accumulation history of oil and gas, and entered the actual application stage comprehensively; currently, it is developing towards three-dimensional models and three-phase multi-component migration.

[0008] In basin simulation, different regions of the same geological horizon may have different sedimentary types. Regarding this difficult problem, current mainstream domestic and foreign software do not provide good solutions. Either a geological horizon is divided into multiple regions for separate independent calculations, losing the overall meaning of the basin; or a geological horizon is turned into multiple geological horizons, each with only one geological type, responsible for one region, and the sedimentary thickness outside the region is set to zero, which does not conform to the actual stratigraphic sedimentation process and is also inaccurate in age division. Therefore, how to give a reasonable solution to the problem of different geological events in the same geological horizon is of great significance for basin simulation. Summary of the Invention

[0009] The purpose of the present invention is to solve the problem that the existing basin simulation software has an unreasonable processing method for geological horizons with different geological events.

[0010] To achieve the above purpose, the present invention provides a basin numerical simulation method and system.

[0011] According to the first aspect of the present invention, a basin numerical simulation method is provided, and the method includes the following steps:

[0012] Obtain the parameters of each geological horizon of the target basin, and the parameters of the geological horizon include age, shape, geological event type, different geological event ranges, and different geological event types;

[0013] In response to an instruction to increase a rectangular different geological event range for the target first-type geological horizon, execute a predetermined first rectangular different geological event range increasing strategy on the target first-type geological horizon to obtain the geological horizon to be numerically simulated and calculated, where the first-type geological horizon is the geological horizon with an empty geological event type and different geological event ranges;

[0014] In response to an instruction to increase the rectangular different geological event ranges for the target second type of geological horizon, a predetermined second rectangular different geological event range increasing strategy is executed on the target second type of geological horizon to obtain a geological horizon to be numerically simulated. The second type of geological horizon is a geological horizon with non-empty geological event types and different geological event ranges.

[0015] In response to an instruction to increase the polygonal different geological event ranges for the target first type of geological horizon, a predetermined first polygonal different geological event range increasing strategy is executed on the target first type of geological horizon to obtain a geological horizon to be numerically simulated.

[0016] In response to an instruction to increase the polygonal different geological event ranges for the target second type of geological horizon, a predetermined second polygonal different geological event range increasing strategy is executed on the target second type of geological horizon to obtain a geological horizon to be numerically simulated.

[0017] In response to an instruction to increase the local erosion thickness for the target geological horizon, according to the type of the target geological horizon and the shape of the corresponding different geological event ranges to be increased, a corresponding different geological event range increasing strategy is executed on the target geological horizon, and the depth values of the corresponding different geological event ranges on the target geological horizon are modified accordingly to obtain a geological horizon to be numerically simulated.

[0018] Perform block numerical simulation calculations on the geological horizon to be numerically simulated.

[0019] Preferably, the execution of the predetermined first rectangular different geological event range increasing strategy on the target first type of geological horizon includes:

[0020] Based on the rectangular different geological event ranges to be increased, the overall range of the target first type of geological horizon is cut into corresponding numbers of rectangular cuts.

[0021] Based on the geological event ranges and geological event types corresponding to each obtained rectangular cut, corresponding different geological event ranges and different geological event types are added to the target first type of geological horizon.

[0022] Based on the different geological event ranges added to the target first type of geological horizon, corresponding different geological event ranges and different geological event types are added to the remaining geological horizons of the target basin, and the different geological event type is the geological event type of the corresponding geological horizon.

[0023] Preferably, the execution of the predetermined second rectangular different geological event range increasing strategy on the target second type of geological horizon includes:

[0024] Judge the intersection of the different geological event ranges of the rectangle to be added with each different geological event range of the target second type of geological horizon:

[0025] If the intersection of the different geological event ranges of the rectangle to be added and a different geological event range is empty, the length of the range of the intersection is less than the length of one grid, or the width of the range of the intersection is less than the width of one grid, it is judged that there is no intersection between the two, and skip;

[0026] If the different geological event ranges of the rectangle to be added contain a different geological event range, the length of the area where the latter is not included by the former is less than the length of one grid, or the width of the area where the latter is not included by the former is less than the width of one grid, it is judged that the former contains the latter, and update the different geological event type corresponding to this different geological event range to the different geological event type corresponding to the different geological event ranges of the rectangle to be added;

[0027] If the length of the area where a different geological event range is included by the different geological event ranges of the rectangle to be added is greater than the length of one grid and the width is greater than the width of one grid, it is judged that there is an intersection between the two, and delete the different geological event range and different geological event type carried by the corresponding area on the target second type of geological horizon, and execute the first rectangle different geological event range addition strategy on this corresponding area.

[0028] Preferably, the execution of the predetermined first polygon different geological event range addition strategy on the target first type of geological horizon includes:

[0029] Give the different geological event ranges of the polygon to be added on the overall range of the target first type of geological horizon;

[0030] Perform hierarchical four-way cutting on the overall range with the different geological event ranges of the polygon to be added:

[0031] If the proportion of the area corresponding to a geological event type in an area within the overall range after one cut reaches a predetermined first ratio, stop the hierarchical four-way cutting of this area, mark this area as this geological event type, and add the corresponding different geological event range and different geological event type to this area; add the corresponding different geological event range and different geological event type to the corresponding areas of the remaining geological horizons of the target basin, and this different geological event type is the geological event type of the corresponding geological horizon;

[0032] If the number of grid rows corresponding to a region within the overall range after the first segmentation is less than the predetermined number of grid rows or the number of grid columns is less than the predetermined number of grid columns, stop the hierarchical four-way segmentation of this region, mark this region as the geological event type with a relatively large proportion in the corresponding region of this region, and add corresponding different geological event ranges and different geological event types to this region; add corresponding different geological event ranges and different geological event types to the corresponding regions of the remaining geological horizons of the target basin, and this different geological event type is the geological event type of the corresponding geological horizon.

[0033] Preferably, the implementation of the predetermined second polygon different geological event range increasing strategy for the target second type of geological horizon includes:

[0034] Judge the intersection of the polygon different geological event range to be added and each different geological event range of the target second type of geological horizon:

[0035] If the intersection of the polygon different geological event range to be added and a different geological event range is empty or the area proportion of the intersection within this different geological event range is less than the predetermined second ratio, skip it;

[0036] If the polygon different geological event range to be added contains a different geological event range or the area proportion of the intersection of the two within this different geological event range is greater than the predetermined third ratio, update the different geological event type corresponding to this different geological event range to the different geological event type corresponding to the polygon different geological event range to be added;

[0037] If the area proportion of the intersection of the polygon different geological event range to be added and the different geological event range within a different geological event range is less than the third ratio, delete the different geological event range and different geological event type carried by the corresponding area on the target second type of geological horizon corresponding to this different geological event range, and implement the first polygon different geological event range increasing strategy for this corresponding area.

[0038] Preferably, the corresponding different geological event range increasing strategy implemented for the target geological horizon includes:

[0039] If the target geological horizon is the first type of geological horizon and the different geological event range to be added is a rectangle, implement the first rectangle different geological event range increasing strategy for the target geological horizon;

[0040] If the target geological horizon is the second type of geological horizon and the different geological event range to be added is a rectangle, implement the second rectangle different geological event range increasing strategy for the target geological horizon;

[0041] If the target geological horizon is a first type of geological horizon and the range of different geological events to be added is a polygon, perform the first polygon different geological event range addition strategy on the target geological horizon;

[0042] If the target geological horizon is a second type of geological horizon and the range of different geological events to be added is a polygon, perform the second polygon different geological event range addition strategy on the target geological horizon.

[0043] Preferably, the first ratio, the second ratio, and the third ratio are 90%, 10%, and 90% respectively.

[0044] Preferably, the performing of the block numerical simulation calculation on the geological horizons to be numerically simulated includes:

[0045] Cut each geological horizon into corresponding blocks based on their respective ranges of different geological events;

[0046] Take each block as a separate geological horizon data;

[0047] Divide the geological horizon data of the geological horizons that have corresponding relationships with each other among all geological horizons into a group, and perform numerical simulation calculation on the geological horizon data of this group.

[0048] Preferably, the performing of the numerical simulation calculation on a group of geological horizon data includes a finite element solution step;

[0049] The finite element solution step includes:

[0050] Obtain the target intermediate calculation result corresponding to each grid point of each block;

[0051] Obtain a finite element matrix based on the obtained target intermediate calculation results;

[0052] Solve the finite element matrix to obtain a solution result;

[0053] Write each obtained solution result back to the corresponding grid point of the corresponding block.

[0054] According to the second aspect of the present invention, there is provided a basin numerical simulation system, which includes a processor and a memory. When the processor executes the computer program stored in the memory, the above-mentioned any basin numerical simulation method is implemented.

[0055] The beneficial effects of the present invention are as follows:

[0056] The basin numerical simulation method of the present invention divides geological horizons with different geological events and sets corresponding geological events for each obtained block. During the basin simulation process, the geological horizons are simulated in blocks and combined simulation is carried out in the steps that require combined calculation, and finally the overall basin simulation process is completed. The basin numerical simulation method of the present invention, on the one hand, ensures the overall significance of the basin, and on the other hand, realizes accurate simulation calculation for different geological event regions of the geological horizons, thus effectively solving the problem that the existing basin simulation software has an unreasonable processing method for geological horizons with different geological events.

[0057] The basin numerical simulation system of the present invention and the above-mentioned basin numerical simulation method belong to a general inventive concept, so they have the same beneficial effects as the above-mentioned basin numerical simulation method.

[0058] Other features and advantages of the present invention will be described in detail in the following specific implementation section. BRIEF DESCRIPTION OF THE DRAWINGS

[0059] By describing the exemplary embodiments of the present invention in more detail in conjunction with the drawings, the above and other objects, features, and advantages of the present invention will become more apparent. Among them, in the exemplary embodiments of the present invention, the same reference numerals generally represent the same components.

[0060] Figure 1 Shows the implementation flowchart of the basin numerical simulation method according to Embodiment 1 of the present invention;

[0061] Figure 2 Shows the geological horizon data structure diagram according to Embodiment 1 of the present invention;

[0062] Figure 3 Shows the schematic diagram of the division of different geological event ranges of rectangles in a single region according to Embodiment 1 of the present invention;

[0063] Figure 4 Shows the schematic diagram of the division of different geological event ranges of rectangles in multiple regions according to Embodiment 1 of the present invention;

[0064] Figure 5 Shows the schematic diagram of the division of different geological event ranges of polygons in a single region according to Embodiment 1 of the present invention;

[0065] Figure 6 Shows the schematic diagram of the division of different geological event ranges of polygons in multiple regions according to Embodiment 1 of the present invention;

[0066] Figure 7 Shows the schematic diagram of the newly added local erosion thickness according to Embodiment 1 of the present invention;

[0067] Figure 8Shows a block diagram based on different geological event ranges according to Embodiment 1 of the present invention;

[0068] Figure 9 Shows a schematic diagram of the correspondence between the block data and the row and column numbers of the original geological horizons according to Embodiment 1 of the present invention. Detailed implementation mode

[0069] The preferred embodiments of the present invention will be described in more detail below. Although the preferred embodiments of the present invention are described below, it should be understood that the present invention can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided to make the present invention more thorough and complete, and to fully convey the scope of the present invention to those skilled in the art.

[0070] Embodiment 1: Figure 1 Shows a flowchart of the implementation of the basin numerical simulation method according to an embodiment of the present invention. Refer to Figure 1 The basin numerical simulation method according to an embodiment of the present invention includes the following steps:

[0071] Obtain the parameters of each geological horizon of the target basin, where the parameters of the geological horizon include age, shape, geological event type, different geological event ranges, and different geological event types;

[0072] In response to an instruction to add a rectangular different geological event range to the target first type of geological horizon, execute a predetermined first rectangular different geological event range addition strategy on the target first type of geological horizon to obtain a geological horizon to be numerically simulated, where the first type of geological horizon is a geological horizon with an empty geological event type and different geological event ranges;

[0073] In response to an instruction to add a rectangular different geological event range to the target second type of geological horizon, execute a predetermined second rectangular different geological event range addition strategy on the target second type of geological horizon to obtain a geological horizon to be numerically simulated, where the second type of geological horizon is a geological horizon with a non-empty geological event type and different geological event ranges;

[0074] In response to an instruction to add a polygonal different geological event range to the target first type of geological horizon, execute a predetermined first polygonal different geological event range addition strategy on the target first type of geological horizon to obtain a geological horizon to be numerically simulated;

[0075] In response to an instruction to add a polygonal different geological event range to the target second type of geological horizon, execute a predetermined second polygonal different geological event range addition strategy on the target second type of geological horizon to obtain a geological horizon to be numerically simulated;

[0076] In response to an instruction to increase the local erosion thickness of a target geological horizon, according to the type of the target geological horizon and the shape of the corresponding different geological event ranges to be increased, corresponding different geological event range increasing strategies are executed on the target geological horizon, and the depth values of the corresponding different geological event ranges on the target geological horizon are modified accordingly, so as to obtain a geological horizon to be numerically simulated and calculated;

[0077] Perform block numerical simulation and calculation on the geological horizon to be numerically simulated and calculated.

[0078] Specifically, in the embodiments of the present invention, Figure 2 shows the geological horizon data structure diagram of the embodiments of the present invention. Refer to Figure 2 , the parameters of each geological horizon include age, griddata, geoevent, diffpartextent, and diffpartgeoevent. Among them, age represents the formation time of the geological horizon, with the unit of million years; griddata is represented by a set of xyz-direction gridded coordinates, where xy of each grid coordinate point represents the geographical location, and z represents the depth in the current basin; the values of geoevent include four types: deposition, erosion, virtual layer, and sedimentary hiatus. Diffpartextent is a set of rectangular ranges or polygonal ranges, and each shape represents the range of a different geological event; diffpartgeoevent is a set of geological event types, and each geological event type represents the event type of a different geological event. The value range of this type is the same as that of geoevent, also including four types: deposition, erosion, virtual layer, and sedimentary hiatus.

[0079] Further, in the embodiments of the present invention, the execution of the predetermined first rectangular different geological event range increasing strategy on the target first type of geological horizon includes:

[0080] Cut the overall range of the target first type of geological horizon based on the rectangular different geological event range to be increased, so as to obtain a corresponding number of rectangular cuts;

[0081] Add corresponding different geological event ranges and different geological event types to the target first type of geological horizon based on the geological event ranges and geological event types corresponding to each obtained rectangular cut;

[0082] Add corresponding different geological event ranges and different geological event types to the remaining geological horizons of the target basin based on the different geological event ranges added to the target first type of geological horizon. The different geological event type is the geological event type of the corresponding geological horizon.

[0083] Further, in the embodiments of the present invention, the execution of the predetermined second rectangular different geological event range increasing strategy for the target second type of geological horizon includes:

[0084] Judging the intersection of the rectangular different geological event range to be added and each different geological event range of the target second type of geological horizon:

[0085] If the intersection of the rectangular different geological event range to be added and a different geological event range is empty, the length of the range of the intersection is less than the length of one grid, or the width of the range of the intersection is less than the width of one grid, it is judged that there is no intersection between the two, and it is skipped;

[0086] If the rectangular different geological event range to be added contains a different geological event range, the length of the area where the latter is not included by the former is less than the length of one grid, or the width of the area where the latter is not included by the former is less than the width of one grid, it is judged that the former contains the latter, and the different geological event type corresponding to this different geological event range is updated to the different geological event type corresponding to the rectangular different geological event range to be added;

[0087] If the length of the area where a different geological event range is included by the rectangular different geological event range to be added is greater than the length of one grid and the width is greater than the width of one grid, it is judged that there is an intersection between the two, and the different geological event range and different geological event type carried by the corresponding area on the target second type of geological horizon corresponding to this different geological event range are deleted, and the first rectangular different geological event range increasing strategy is executed for this corresponding area.

[0088] Further, in the embodiments of the present invention, the execution of the predetermined first polygonal different geological event range increasing strategy for the target first type of geological horizon includes:

[0089] Give the polygonal different geological event range to be added on the overall range of the target first type of geological horizon;

[0090] Perform step-by-step four equal divisions on the overall range with the polygonal different geological event range to be added;

[0091] If the proportion of the area corresponding to a geological event type in an area within the overall range after one division reaches a predetermined first ratio, stop the step-by-step four equal division of this area, mark this area as this geological event type, and add the corresponding different geological event range and different geological event type to this area; add the corresponding different geological event range and different geological event type to the corresponding areas of the remaining geological horizons of the target basin, and this different geological event type is the geological event type of the corresponding geological horizon;

[0092] If the number of grid rows corresponding to a region within the overall range after the first segmentation is less than the predetermined number of grid rows or the number of grid columns is less than the predetermined number of grid columns, stop the hierarchical four-way segmentation of this region, mark this region as the geological event type with a larger proportion in the corresponding region of this region, and add corresponding different geological event ranges and different geological event types to this region; add corresponding different geological event ranges and different geological event types to the corresponding regions of the remaining geological horizons of the target basin, and this different geological event type is the geological event type of the corresponding geological horizon.

[0093] Furthermore, in the embodiments of the present invention, the execution of the predetermined second polygon different geological event range increasing strategy for the target second type of geological horizon includes:

[0094] Perform an intersection judgment on the polygon different geological event range to be added and each different geological event range of the target second type of geological horizon:

[0095] If the intersection of the polygon different geological event range to be added and a different geological event range is empty or the area ratio of the intersection within this different geological event range is less than the predetermined second ratio, skip it;

[0096] If the polygon different geological event range to be added contains a different geological event range or the area ratio of the intersection between the two within this different geological event range is greater than the predetermined third ratio, update the different geological event type corresponding to this different geological event range to the different geological event type corresponding to the polygon different geological event range to be added;

[0097] If the area ratio of the intersection of the polygon different geological event range to be added and the different geological event range within a different geological event range is less than the third ratio, delete the different geological event range and the different geological event type carried by the corresponding area on the target second type of geological horizon corresponding to this different geological event range, and perform the first polygon different geological event range increasing strategy on this corresponding area.

[0098] Furthermore, in the embodiments of the present invention, the execution of the corresponding different geological event range increasing strategy for the target geological horizon includes:

[0099] If the target geological horizon is the first type of geological horizon and the different geological event range to be added is a rectangle, perform the first rectangle different geological event range increasing strategy on the target geological horizon;

[0100] If the target geological horizon is the second type of geological horizon and the different geological event range to be added is a rectangle, perform the second rectangle different geological event range increasing strategy on the target geological horizon;

[0101] If the target geological horizon is a first type of geological horizon and the range of different geological events to be added is a polygon, execute the first polygon different geological event range adding strategy on the target geological horizon;

[0102] If the target geological horizon is a second type of geological horizon and the range of different geological events to be added is a polygon, execute the second polygon different geological event range adding strategy on the target geological horizon.

[0103] Furthermore, in the embodiments of the present invention, the first ratio, the second ratio, and the third ratio are 90%, 10%, and 90% respectively.

[0104] Specifically, in the embodiments of the present invention, for a geological horizon, there may not be only a single type of geological event in its entire area. For example, there may be some areas with deposition and some areas with erosion. For such a situation, it is necessary to add different geological event ranges and different geological event types to represent these local areas. Different geological event ranges and different geological event types appear in pairs, representing the range and geological type of a local area respectively. The specific definitions and range division methods of different geological events include rectangles and polygons, and whether there are existing different geological event ranges and other situations. The range division and definition of different geological events are divided into four types, and the processing method for each type is as follows:

[0105] I. Initial geological horizon, adding a rectangular different geological event range:

[0106] Figure 3 Shows a schematic diagram of the rectangular different geological event range division in a single area of the embodiments of the present invention. Refer to Figure 3 For a geological horizon, the geological event type is deposition, and the initial different geological event range and different geological event type are empty, indicating that there are no different geological events yet. Now there is a local range in this geological horizon that needs to be set as erosion. First, based on this local range, perform rectangular cutting on the entire range of this geological horizon:

[0107] For Figure 3For the erosion range in case a, after cutting, 5 rectangular ranges can be obtained. Then, add 5 different geological event ranges to this geological horizon. Each different geological event range is the range of one of the rectangles, and add 5 corresponding different geological event types, namely 1 erosion and 4 depositions, to represent the geological event types corresponding to these 5 rectangular ranges. At the same time, add these 5 different geological event ranges to each of the other geological horizons, and add 5 corresponding geological event types, with the values being the original geological event types of each geological horizon. For example, if the geological event type of one of the other geological horizons is sedimentary hiatus, then add 5 different geological event types with the value of sedimentary hiatus to this geological horizon.

[0108] For Figure 3 For the erosion ranges in cases b, c, d, and e, when some boundaries of the erosion range are close to the boundaries of the geological horizon range, that is, the distance length is less than the length of a geological horizon grid or the distance width is less than the width of a geological horizon grid, the part close to the boundary cannot be cut into independent rectangles. Therefore, 4, 3, 3, and 2 rectangular ranges are obtained respectively. The method of rectangular cutting is similar to that in case a.

[0109] Case b will generate 4 different geological event ranges and different geological event types, namely 1 erosion and 3 depositions. At the same time, add these 4 different geological event ranges to each of the other geological horizons in case b, and add 4 corresponding different geological event types, with the values being the original geological event types of each geological horizon.

[0110] Case c will generate 3 different geological event ranges and different geological event types, namely 1 erosion and 2 depositions. At the same time, add these 3 different geological event ranges to each of the other geological horizons in case c, and add 3 corresponding different geological event types, with the values being the original geological event types of each geological horizon.

[0111] Case d will generate 3 different geological event ranges and different geological event types, namely 1 erosion and 2 depositions. At the same time, add these 3 different geological event ranges to each of the other geological horizons in case d, and add 3 corresponding different geological event types, with the values being the original geological event types of each geological horizon.

[0112] Case e will generate 2 different geological event ranges and different geological event types, namely 1 erosion and 1 deposition. At the same time, add these 2 different geological event ranges to each of the other geological horizons in case e, and the values are the geological event types of each geological horizon itself.

[0113] The purpose of adding the same range of different geological events to other geological horizons is to keep the division of the range of different geological events in each geological horizon consistent, so as to perform block-by-block calculations according to the range in the follow-up. Since this newly added different geological event only applies to the current geological horizon, the values of the newly added different geological event types in other geological horizons are still the original geological event types of each geological horizon.

[0114] II. For the geological horizons with existing ranges of different geological events, add rectangular ranges of different geological events:

[0115] Figure 4 The schematic diagram of the division of the rectangular range of different geological events in multiple regions according to the embodiment of the present invention is shown. Refer to Figure 4 , when the range of different geological events and the type of different geological events in the geological horizon are not empty, the newly added range of different geological events needs to be compared one by one with the existing ranges of different geological events. Figure 4 In , the current geological horizon already has a total of 9 different geological event ranges from a to i. At this time, if a new geological event range is added, the newly added geological event range will be judged for intersection with the different geological event ranges a to i.

[0116] If there is no intersection, or the length of the intersection range is less than the length of one grid, or the width of the intersection range is less than the width of one grid, it is considered that there is no intersection and it will be directly skipped, such as the different geological event range d;

[0117] If the newly added geological event range can entirely contain an existing range of different geological events, or the length of the unincluded part is less than the length of one grid, or the width of the unincluded part is less than the width of one grid, it is considered to be entirely contained, and the type of different geological events corresponding to this existing range of different geological events will be modified to the type of the newly added different geological event. For example, for the different geological event range f, the corresponding type of different geological events needs to be changed to the type of the newly added different geological event;

[0118] If the ranges intersect, such as the different geological event ranges a, b, c, e, g, h, i, first delete the current range of different geological events and the type of different geological events of this existing range, and then, in the way of adding rectangular ranges of different geological events to the initial geological horizon as described above, divide the existing range into new rectangular blocks, and generate the geological event ranges and types corresponding to the new rectangular blocks.

[0119] For example, for different geological event range a, which has a corner intersection with the newly added different geological event ranges, delete the different geological event ranges and different geological event types corresponding to different geological event range a of the current geological horizon, and then according to the aforementioned rectangular range division method, three new different geological event ranges a1, a2, and a3 can be generated, and add the different geological event ranges and different geological event types corresponding to the newly added different geological event ranges a1, a2, and a3. Among them, the values of the different geological event ranges are the ranges of these three regions respectively, the value of the different geological event type corresponding to different geological event range a1 is the type of the newly added geological event, and the values of the different geological event types corresponding to different geological event ranges a2 and a3 are the values of the different geological event types corresponding to different geological event range a. For other geological horizons, similarly delete the different geological event ranges and different geological event types corresponding to different geological event range a, and add the different geological event ranges and different geological event types of the three regions a1, a2, and a3. Among them, the values of the different geological event ranges are the ranges of these three regions respectively, and the values of the different geological event types are the values of the different geological event types corresponding to different geological event range a of each geological horizon.

[0120] For example, for different geological event range b, which has a horizontal bar intersection with the newly added different geological event ranges, delete the different geological event ranges and different geological event types corresponding to different geological event range b of the current geological horizon, and then according to the aforementioned rectangular range division method, two new different geological event ranges b1 and b2 can be generated, and add the different geological event ranges and different geological event types corresponding to the newly added different geological event ranges b1 and b2. Among them, the values of the different geological event ranges are the ranges of these two regions respectively, the value of the different geological event type corresponding to different geological event range b1 is the type of the newly added geological event, and the value of the different geological event type corresponding to different geological event range b2 is the value of the different geological event type corresponding to different geological event range b. For other geological horizons, similarly delete the different geological event ranges and different geological event types corresponding to different geological event range b, and add the different geological event ranges and different geological event types of the two regions b1 and b2. Among them, the values of the different geological event ranges are the ranges of these two regions respectively, and the values of the different geological event types are the values of the different geological event types corresponding to different geological event range b of each geological horizon.

[0121] For example, different geological event ranges h and i have a corner intersection with the newly added different geological event ranges, and the bottom boundary of the intersection part is narrow, less than the length of a geological horizon grid in the corresponding direction. Therefore, this part is regarded as unable to be cut into independent rectangles. According to the method of adding rectangular different geological event ranges to the initial geological horizon, the different geological event range h is divided into different geological event ranges h1 and h2, and the different geological event range i is divided into different geological event ranges i1 and i2.

[0122] The processing methods for other different geological event ranges are the same as those for different geological event ranges h and i. After all the processing is completed, the different geological event ranges and different geological event types corresponding to different geological event ranges a, b, c, e, g, h, i are deleted, and the different geological event ranges and different geological event types corresponding to different geological event ranges a1, a2, a3, b1, b2, c1, c2, c3, e1, e2, g1, g2, h1, h2, i1, i2 are newly added.

[0123] III. Initial geological horizon, adding polygonal different geological event ranges:

[0124] When the newly added different geological event range is a polygon, it is necessary to divide the area where the polygon is located into four equal parts level by level until the area of a certain geological type in the divided area is greater than 90% of the current divided area. This divided area as a whole is recorded as this geological type, forming a pair of different geological event ranges and different geological event types; or if the number of rows and columns of the next-level divided area is less than 2 rows or 2 columns, the four equal division is no longer continued, and the geological type with the largest area proportion in the current divided area is recorded as the overall geological type of the current divided area, forming a pair of different geological event ranges and different geological event types.

[0125] Figure 5 It shows a schematic diagram of the division of polygonal different geological event ranges in a single area of an embodiment of the present invention. Refer to Figure 5, for a geological horizon with a sedimentary geological event type, the initial different geological event ranges and different geological event types are empty, indicating that there are no different geological events in this geological horizon. Now, a local polygon range in this geological horizon needs to be set as erosion, as shown in step1 in the figure. First, divide the rectangular range of the geological horizon into four equal parts, as shown in step2 in the figure, to obtain four sub-regions a, b, c, and d; then judge these four regions in sequence. It is found that the area of the different geological event range in region a exceeds 90% of the area of region a. Then, set the entire region a as erosion, as shown in step3 in the figure. That is, add 1 different geological event range to the current geological horizon, with the value being the range of region a, and at the same time add a different geological event type, with the type being erosion; regions b, c, and d do not meet the 90% condition at this time. Further divide regions b, c, and d into four equal parts respectively, as shown in step4 in the figure, and then judge whether there is a situation where the area of one geological type exceeds 90% in the regions after the division of regions b, c, and d. As shown in step5 in the figure, for region b, there are 2 regions that meet the conditions after division. One region is that the area of the original geological type exceeds 90%, and the other region is that the area of the newly added different geological event exceeds 90%. Then, add 2 different geological event ranges to the current geological horizon, with the values being the areas of these two regions respectively, and at the same time add 2 different geological event types, with the types being sedimentation and erosion respectively; for regions c and d, there is 1 region that meets the conditions after division, that is, the area of erosion exceeds 90%. Then, add 2 different geological event ranges to the current geological horizon, with the values being the areas of these two regions respectively, and at the same time add 2 different geological event types, with the values both being erosion; next, continue to divide the remaining regions that do not meet the area requirements into four equal parts, as shown in step6 in the figure, until all regions meet the area ratio requirements or the number of rows and columns at the next level is less than 2. Add the geological types and ranges of these regions to the different geological event ranges and different geological event types of the geological horizon.

[0126] When adding different geological event ranges and different geological event types to the current geological horizon, the same range of different geological event ranges is also added to each other geological horizon, and the values of the corresponding different geological event types are the original geological event types of each geological horizon, so as to keep the division of different geological event ranges of each geological horizon always consistent.

[0127] IV. For a geological horizon with existing different geological event ranges, adding a polygonal different geological event range:

[0128] When the different geological event ranges and different geological event types of a geological horizon are not empty, the newly added polygonal different geological event range needs to be compared with the existing different geological event ranges one by one.

[0129] Figure 6The figure shows a schematic diagram of the division of different geological event ranges of polygons in multiple regions. Refer to Figure 6 , the current geological horizon already has 6 different geological event ranges from a to f. At this time, if a new polygon with different geological event ranges is added, the newly added polygon with different geological event ranges will be judged for intersection with the different geological event ranges from a to f.

[0130] If there is no intersection, or the proportion of the intersection area is less than 10%, it will be skipped; if the newly added polygon with different geological event ranges can wholly contain an existing different geological event range, or the proportion of the intersection area is greater than 90%, then the type of the different geological event corresponding to the existing different geological event range will be modified to the type of the newly added different geological event; if the ranges intersect and no geological type has an area proportion exceeding 90%, then the different geological event range and the type of the different geological event corresponding to the existing different geological event range will be deleted, and then the existing different geological event range will be equally divided into four levels step by step in the manner of adding a polygon with different geological event ranges to the initial geological horizon until the area of a certain geological type in the divided area is greater than 90% of the current divided area, or the number of rows or columns of the next-level divided area is less than 2 and no further equal division is carried out, and the geological type with the largest area proportion in the current divided area will be recorded as the overall geological type of the current divided area.

[0131] For example, for the different geological event range a, it intersects with a corner of the newly added different geological event range. As shown in step2 in the figure, the area occupied by the original geological type in the different geological event range a exceeds 90%. For this area, the original different geological event range and the type of the different geological event of the different geological event range a can be used, and there is no need to add or modify the different geological event range and the type of the different geological event. This area is processed.

[0132] For example, for the different geological event range e, it intersects with the newly added different geological event range, and the proportion of the area of the newly added different geological event range in the different geological event range e is greater than 90%. Then, the type of the different geological event corresponding to the different geological event range e can be modified to the type of the newly added different geological event, and there is no need to add or modify the different geological event range. This area is processed.

[0133] For the different geological event ranges a and e, they are existing different geological event ranges and no new ranges are divided. Therefore, the different geological event ranges have not changed, and no changes need to be made to other geological horizons. After comparing the different geological event ranges from a to f with the newly added different geological event range one by one, the processing results are shown in step3 in the figure. There are still different geological event ranges b, c, d, and f that have not been processed and need to start equal division into four levels, as shown in step4.

[0134] The method of successively dividing each region into four equal parts refers to the method of successively dividing the different geological event ranges of the polygon by increasing the initial geological horizon as described above. First, temporarily record the different geological event types corresponding to this region in the current geological horizon and each of the other geological horizons, and then delete the different geological event ranges and different geological event types corresponding to this region in all geological horizons. Next, divide this region into four equal parts successively according to the method of successively dividing the different geological event ranges of the polygon by increasing the initial geological horizon, generating new geological event ranges and types.

[0135] For example, for different geological event range b, after dividing it into four equal parts, judge these four regions in turn. In 3 regions, the area proportion of a certain geological type exceeds 90%. As shown in step5, 2 of them are the original geological event types of different geological event range b, and 1 is the type of newly added different geological event. Add 3 different geological event ranges to the current geological horizon, with the values being the areas of these 3 regions, and add 3 different geological event types. Among them, 2 values are the geological event types of themselves, and 1 is the type of newly added different geological event. For each of the other geological horizons, add 3 identical different geological event ranges and add 3 different geological event types, and the values are all the values of the different geological event types corresponding to the different geological event range b of each geological horizon.

[0136] For example, for different geological event range c, after dividing it into four equal parts, judge these four regions in turn. In 1 region, the area proportion of a certain geological type exceeds 90%. As shown in step5, it is the type of newly added different geological event. Add 1 different geological event range to the current geological horizon, with the value being the area of this region, and add 1 different geological event type, with the value being the type of newly added different geological event. For each of the other geological horizons, add 1 identical different geological event range and add 1 different geological event type, and the value is the value of the different geological event type corresponding to the different geological event range c of each geological horizon.

[0137] The processing methods of different geological event ranges d and f are similar to those of different geological event ranges b and c. When the processing of different geological event ranges b, c, d, and f is completed, as shown in step6, there are still some blank grids left unprocessed. Further divide them into four equal parts until all grids are processed.

[0138] Specifically, in the embodiments of the present invention, in addition to directly dividing the local geological event range, the local geological event range can also be divided by inputting the local erosion thickness. Figure 7 Shows a schematic diagram of the newly added local erosion thickness in the embodiments of the present invention. Refer to Figure 7, there are already three geological horizons h1, h2, and h3. Now, a local erosion thickness is added to the geological horizons h1, h2, and h3. That is, within the range corresponding to the local erosion thickness of these three geological horizons, local different geological events occur. At this time, on the one hand, the range of local geological events needs to be increased, and on the other hand, the depth values of the geological horizons h1, h2, and h3 in the local range need to be changed. The method of increasing the range of local geological events is the same as the method of increasing the range of different geological events described above. According to the range type of the local erosion thickness, whether it is a rectangle or a polygon, and then according to whether there is already a range of different geological events, it corresponds to the above four situations for corresponding processing. For changing the depth values of the geological horizons h1, h2, and h3 in the local range, the processing method is as Figure 7 shown. For the geological horizon h2, within the range corresponding to the local erosion, the depth value of each point is the depth value of the geological horizon h3 at the same position plus the thickness value of the local erosion thickness at the same position. That is, adding this local erosion thickness to the geological horizon h3 is the depth of the geological horizon h2 at this position. The local geological type corresponding to the geological horizon h2 within this local erosion thickness range is a virtual layer; for the geological horizon h1, within the range corresponding to the local erosion, the depth value of each point is the depth value of the geological horizon h3 at the same position. That is, all the geological horizons within this range are eroded away to the position of the geological horizon h3. The local geological type corresponding to the geological horizon h1 within this local erosion thickness range is an erosion layer. The newly generated geological horizons h1 and h2 are shown in the figure.

[0139] Furthermore, in the embodiment of the present invention, the block numerical simulation calculation of the geological horizons to be numerically simulated includes:

[0140] Cutting each geological horizon into corresponding blocks based on their respective different geological event ranges;

[0141] Regarding each block as a separate geological horizon data;

[0142] Dividing the geological horizon data of all geological horizons that have a corresponding relationship with each other into a group, and performing numerical simulation calculation on the geological horizon data of this group.

[0143] Furthermore, in the embodiment of the present invention, the numerical simulation calculation of a group of geological horizon data includes a finite element solution step;

[0144] The finite element solution step includes:

[0145] Obtaining the target intermediate calculation result corresponding to each grid point of each block;

[0146] Obtaining a finite element matrix based on the obtained target intermediate calculation result;

[0147] Solve the finite element matrix to obtain a solution result;

[0148] Write each obtained solution result back to the corresponding grid points of the corresponding block.

[0149] Specifically, in the embodiments of the present invention, for the same horizon, the geological types in different block areas may be different. During the basin simulation calculation, the calculation results obtained for different geological types are different, so it is necessary to perform block calculation.

[0150] If there is no range of different geological events, no block division is performed, and this set of geological horizons is directly used as a whole for simulation calculation; if there is a range of different geological events, this set of geological horizons is divided according to the range of different geological events, and each block corresponds to the geological horizon data within a different range of geological events.

[0151] Figure 8 shows the block diagram based on different ranges of geological events in the embodiments of the present invention. Refer to Figure 8 , the shape data of a geological horizon is divided according to different ranges of geological events. After division, each part is used as a separate geological horizon data object, whose age is the same as that of the original geological horizon, the shape is the grid coordinate points contained in this block area, the geological event type is the different geological event types corresponding to this block in the original geological horizon, and the different ranges of geological events and different geological event types are empty.

[0152] After each geological horizon is divided into blocks, the geological horizon data under the same block are grouped together for simulation calculation. For example, Figure 8 in the splitting method of, a total of 6 groups are formed, and each group includes 4 geological horizon blocks.

[0153] Next, the block data is merged for finite element calculation. The process of basin simulation includes a step of finite element solution, which requires all the data of each geological horizon to be put together to form a large matrix for solution. Therefore, it is necessary to obtain the current intermediate calculation results of each geological horizon from each geological horizon block and merge them to form a matrix. The order of the elements in the matrix is first in the order of geological horizons, then in the row numbers of the coordinate grid points in each geological horizon, and finally in the column numbers of the grid coordinate points. To obtain the intermediate calculation result corresponding to a certain coordinate grid point, it is necessary to find the block where it is located and its position in the block. When dividing the geological horizons into blocks, each block additionally records the position of the starting coordinate point of this block in the grid of the original geological horizon. As Figure 9 shown, the position of the starting coordinate point of block b in the grid of the original geological horizon is row m and column n, and the number of rows and columns of block b are p rows and q columns respectively. Then the row and column range of block b in the grid of the original geological horizon is from row m to row m + p and from column n to column n + q.

[0154] When obtaining the intermediate calculation results of a coordinate point, it is necessary to traverse each block to find the block whose range can include the row and column numbers of this coordinate point, and then calculate the position of this coordinate point in this block. For example, if you want to find the coordinate point with row number i and column number j, after traversing each sub-block, it is found that the range of b can include this coordinate point, and its relative position in block b is the (i - m)-th row and (j - n)-th column of b. After obtaining the intermediate results of all grid points of all geological horizons, perform finite element matrix solution and obtain the calculation results.

[0155] Next, update the calculation results to each block. After this finite element calculation step, it is necessary to continue to complete other basin simulation calculation processes in blocks, so it is necessary to write the calculation results of this step back to the corresponding geological horizon blocks again.

[0156] Similar to the method of corresponding row and column numbers in block data merging, traverse each finite element calculation result. Through their positions in the matrix, it is possible to obtain which geological horizon they belong to and their row and column positions in this geological horizon. Through the original row and column numbers of the geological horizon, obtain the corresponding block and the specific row and column positions in the block, write the finite element calculation results back to the coordinate points in these blocks, and then complete the remaining basin simulation calculations.

[0157] Next, obtain the basin simulation results from each block. Similar to the method of corresponding row and column numbers in block data merging, when it is necessary to obtain the calculation results at a certain position, through the original row and column numbers of the geological horizon, obtain the corresponding block and the specific position in the block, and then obtain the corresponding calculation results from the block.

[0158] So far, the implementation process of different geological events in the same horizon in basin simulation is completed.

[0159] The implementation process of the basin numerical simulation method of the embodiment of the present invention is the initial geological horizon, adding rectangular different geological event ranges; geological horizons with existing different geological event ranges, adding rectangular different geological event ranges; initial geological horizon, adding polygonal different geological event ranges; geological horizons with existing different geological event ranges, adding polygonal different geological event ranges; adding local erosion thickness; performing block simulation based on different geological event ranges.

[0160] According to the basin numerical simulation method of the embodiment of the present invention, correct simulation calculation results can be obtained for each different geological event area of the geological horizon; block simulation of polygonal different geological event areas is realized, and the accuracy of processing different geological event area ranges is improved.

[0161] Embodiment 2: Based on the basin numerical simulation method proposed in Embodiment 1, an embodiment of the present invention correspondingly proposes a basin numerical simulation system, which includes a processor and a memory. When the processor executes the computer program stored in the memory, the basin numerical simulation method proposed in Embodiment 1 is implemented.

[0162] The embodiments of the present invention have been described above. The above description is exemplary and not exhaustive, and is also not limited to the disclosed embodiments. Many modifications and variations are obvious to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A basin numerical simulation method, characterized in that, Including: Obtaining parameters of each geological horizon of a target basin, where the parameters of the geological horizon include age, shape, geological event type, different geological event ranges, and different geological event types; In response to an instruction to increase a rectangular different geological event range for a target first-type geological horizon, performing a predetermined first rectangular different geological event range increasing strategy on the target first-type geological horizon to obtain a geological horizon to be numerically simulated, where the first-type geological horizon is a geological horizon with an empty geological event type and different geological event ranges; In response to an instruction to increase a rectangular different geological event range for a target second-type geological horizon, performing a predetermined second rectangular different geological event range increasing strategy on the target second-type geological horizon to obtain a geological horizon to be numerically simulated, where the second-type geological horizon is a geological horizon with a non-empty geological event type and different geological event ranges; In response to an instruction to increase a polygonal different geological event range for a target first-type geological horizon, performing a predetermined first polygonal different geological event range increasing strategy on the target first-type geological horizon to obtain a geological horizon to be numerically simulated; In response to an instruction to increase a polygonal different geological event range for a target second-type geological horizon, performing a predetermined second polygonal different geological event range increasing strategy on the target second-type geological horizon to obtain a geological horizon to be numerically simulated; In response to an instruction to increase the local erosion thickness for a target geological horizon, according to the type of the target geological horizon and the shape of the corresponding different geological event range to be increased, performing the corresponding different geological event range increasing strategy on the target geological horizon, and correspondingly modifying the depth values of the corresponding different geological event ranges on the target geological horizon to obtain a geological horizon to be numerically simulated; Performing block numerical simulation calculations on the geological horizon to be numerically simulated.

2. The basin numerical simulation method according to claim 1, characterized in that, The performing a predetermined first rectangular different geological event range increasing strategy on the target first-type geological horizon includes: Cutting the overall range of the target first-type geological horizon based on the rectangular different geological event range to be increased to obtain a corresponding number of rectangular cuts; Adding the corresponding different geological event ranges and different geological event types to the target first-type geological horizon based on the geological event ranges and geological event types corresponding to each rectangular cut obtained; Adding the corresponding different geological event ranges and different geological event types to the remaining geological horizons of the target basin based on the different geological event ranges added to the target first-type geological horizon, where the different geological event type is the geological event type of the corresponding geological horizon.

3. The basin numerical simulation method according to claim 2, characterized in that, The performing a predetermined second rectangular different geological event range increasing strategy on the target second-type geological horizon includes: Judging the intersection of the rectangular different geological event range to be increased with each different geological event range of the target second-type geological horizon: If the intersection of the different geological event ranges of the rectangle to be added and a different geological event range is empty, the length of the range of the intersection is less than the length of one grid, or the width of the range of the intersection is less than the width of one grid, it is determined that there is no intersection between the two, and it is skipped; If the different geological event range of the rectangle to be added contains a different geological event range, the length of the area of the latter that is not contained by the former is less than the length of one grid, or the width of the area of the latter that is not contained by the former is less than the width of one grid, it is determined that the former contains the latter, and the different geological event type corresponding to the different geological event range is updated to the different geological event type corresponding to the different geological event range of the rectangle to be added; When the length of the area of a different geological event range contained by the different geological event range of the rectangle to be added is greater than the length of one grid and the width is greater than the width of one grid, it is determined that there is an intersection between the two, and the different geological event range and the different geological event type carried by the corresponding area on the target second type of geological horizon of the different geological event range are deleted, and the first rectangle different geological event range addition strategy is executed for the corresponding area.

4. The basin numerical simulation method according to claim 3, characterized in that, The execution of the predetermined first polygon different geological event range addition strategy for the target first type of geological horizon includes: Giving the polygon different geological event range to be added on the overall range of the target first type of geological horizon; Performing hierarchical four-way division on the overall range with the polygon different geological event range to be added: If the proportion of the area corresponding to a geological event type in a region within the overall range after one division reaches a predetermined first ratio, stop performing hierarchical four-way division on this region, mark this region as this geological event type, and add the corresponding different geological event range and different geological event type to this region; add the corresponding different geological event range and different geological event type to the corresponding regions of the remaining geological horizons of the target basin, and this different geological event type is the geological event type of the corresponding geological horizon; If the number of grid rows corresponding to a region within the overall range after one division is less than the predetermined number of grid rows or the number of grid columns is less than the predetermined number of grid columns, stop performing hierarchical four-way division on this region, mark this region as the geological event type with the larger proportion of the corresponding area in this region, and add the corresponding different geological event range and different geological event type to this region; add the corresponding different geological event range and different geological event type to the corresponding regions of the remaining geological horizons of the target basin, and this different geological event type is the geological event type of the corresponding geological horizon.

5. The basin numerical simulation method according to claim 4, characterized in that, The execution of the predetermined second polygon different geological event range addition strategy for the target second type of geological horizon includes: Judging the intersection of the polygon different geological event range to be added and each different geological event range of the target second type of geological horizon: If the intersection of the polygon different geological event range to be added and a different geological event range is empty or the area proportion of the intersection within this different geological event range is less than a predetermined second ratio, it is skipped; If the different geological event range of the polygon to be added includes a different geological event range or the proportion of the area of the intersection of the two within the different geological event range is greater than a predetermined third ratio, update the different geological event type corresponding to the different geological event range to the different geological event type corresponding to the different geological event range of the polygon to be added; If the proportion of the area of the intersection of the different geological event range of the polygon to be added within a different geological event range and the different geological event range is less than the third ratio, delete the different geological event range and the different geological event type of the corresponding area on the target second type of geological horizon corresponding to the different geological event range, and execute the first polygon different geological event range addition strategy on the corresponding area.

6. The basin numerical simulation method according to claim 5, characterized in that, The corresponding different geological event range addition strategy executed on the target geological horizon includes: If the target geological horizon is the first type of geological horizon and the different geological event range to be added is a rectangle, execute the first rectangle different geological event range addition strategy on the target geological horizon; If the target geological horizon is the second type of geological horizon and the different geological event range to be added is a rectangle, execute the second rectangle different geological event range addition strategy on the target geological horizon; If the target geological horizon is the first type of geological horizon and the different geological event range to be added is a polygon, execute the first polygon different geological event range addition strategy on the target geological horizon; If the target geological horizon is the second type of geological horizon and the different geological event range to be added is a polygon, execute the second polygon different geological event range addition strategy on the target geological horizon.

7. The basin numerical simulation method according to claim 6, characterized in that, The first ratio, the second ratio and the third ratio are 90%, 10% and 90% respectively.

8. The basin numerical simulation method according to claim 7, characterized in that, The block numerical simulation calculation of the geological horizon to be numerically simulated includes: Based on their respective different geological event ranges, each geological horizon is divided into corresponding blocks; Each block is used as a separate geological horizon data; The geological horizon data with corresponding relationships among all geological horizons are divided into a group, and numerical simulation calculation is performed on the geological horizon data of this group.

9. The basin numerical simulation method according to claim 8, characterized in that, Performing numerical simulation calculation on a group of geological horizon data includes a finite element solution step; The finite element solution step includes: Obtain the target intermediate calculation result corresponding to each grid point of each block; Obtain the finite element matrix based on the obtained target intermediate calculation result; Solve the finite element matrix to obtain the solution result; Write each obtained solution result back to the corresponding grid point of the corresponding block.

10. A basin numerical simulation system, characterized in that, Including a processor and a memory, when the processor executes the computer program stored in the memory, it implements the basin numerical simulation method according to any one of claims 1-9.

Images