Numerical simulation method and device for fracture-developed zone, electronic equipment and storage medium

Abstract

The embodiment of the invention discloses a numerical simulation method and device for a fracture-developed zone, electronic equipment and a storage medium. The method comprises the steps that target number distribution obeyed by the number of weakening points on a target surface in a geologic body numerical model is determined according to the geologic body numerical model corresponding to a geologic body model inlaid with a fracture development zone, and the weakening points are used for representing the fracture development zone; according to the target number distribution, weakening point probability distribution obeyed by the probability that any point on the target surface is the weakening point is obtained, and according to the weakening point probability distribution and the target number distribution, the total probability that any point on the target surface is the weakening point is obtained; and obtaining a numerical simulation result of the fracture-developed zone according to the total probability. According to the technical scheme provided by the embodiment of the invention, the accuracy of numerical simulation of the fracture-developed zone can be ensured under a relatively small calculation amount.

Description

Technical Field

[0001] The embodiments of the present invention relate to the field of geophysical exploration technology, and in particular to a numerical simulation method, device, electronic equipment and storage medium for a fracture development zone. Background Art

[0002] Numerical simulation of fracture development zones is crucial in many aspects, such as numerical simulation of three-dimensional geological body seismic forward modeling, oil and gas exploration, and understanding of underground fracture response laws.

[0003] It should be noted that the scale of cracks is generally centimeters to decimeters in width and meters in extension length, while the scale of conventional faults is meters or even kilometers in fault distance and extension length.

[0004] Due to the small scale of fractures, the current numerical simulation schemes used for fracture development zones have problems such as low simulation accuracy or large amount of calculation, which need to be improved. Summary of the invention

[0005] The embodiments of the present invention provide a numerical simulation method, device, electronic device and storage medium for a fracture development zone, so as to ensure the accuracy of the numerical simulation of the fracture development zone with a relatively small amount of calculation.

[0006] According to one aspect of the present invention, a numerical simulation method for a fracture development zone is provided, which may include:

[0007] For a geological body numerical model corresponding to a geological body model inlaid with fracture development zones, determining a target quantity distribution obeyed by the number of weak points on a target surface in the geological body numerical model, wherein the weak points are used to characterize the fracture development zones;

[0008] According to the target quantity distribution, the probability distribution of the weak point that any point on the target surface is a weak point is obtained, and according to the weak point probability distribution and the target quantity distribution, the total probability that any point on the target surface is a weak point is obtained;

[0009] According to the total probability, the numerical simulation results of the fracture development zone are obtained.

[0010] According to another aspect of the present invention, a numerical simulation device for a fracture development zone is provided, which may include:

[0011] A target quantity distribution determination module is used to determine the target quantity distribution obeyed by the number of weak points on the target surface in the geological body numerical model corresponding to the geological body model inlaid with the fracture development zone, wherein the weak points are used to characterize the fracture development zone;

[0012] A total probability obtaining module is used to obtain the probability distribution of weak points obeyed by the probability of any point on the target surface being a weak point according to the target quantity distribution, and to obtain the total probability of any point on the target surface being a weak point according to the weak point probability distribution and the target quantity distribution;

[0013] The numerical simulation result obtaining module can be used to obtain the numerical simulation results of the fracture development zone according to the total probability.

[0014] According to another aspect of the present invention, there is provided an electronic device, which may include:

[0015] at least one processor; and

[0016] a memory communicatively connected to at least one processor; wherein,

[0017] The memory stores a computer program that can be executed by at least one processor. The computer program is executed by at least one processor so that the numerical simulation method of the fracture development zone provided by any embodiment of the present invention is implemented when the at least one processor executes the computer program.

[0018] According to another aspect of the present invention, a computer-readable storage medium is provided, on which computer instructions are stored, and the computer instructions are used to enable a processor to implement the numerical simulation method of the fracture development zone provided by any embodiment of the present invention when executed.

[0019] The technical solution of the embodiment of the present invention is aimed at the geological body numerical model corresponding to the geological body model inlaid with the fracture development zone, by determining the target number distribution obeyed by the number of weakening points on the target surface in the geological body numerical model for characterizing the fracture development zone; then, according to the target number distribution, the probability distribution of weakening points obeyed by the probability that any point on the target surface is a weakening point is obtained, and according to the probability distribution of weakening points and the target number distribution, the total probability that any point on the target surface is a weakening point is obtained; finally, the numerical simulation result of the fracture development zone is obtained according to the total probability. The above technical solution, by considering the small-scale fracture development zone as the weakening point presented in the local mean geological body, adopts the random probability statistical algorithm to express the comprehensive weakening effect of the fracture development zone, and converts the fracture end point into the contribution of the weakening point of the geological body numerical model, thereby ensuring the accuracy of the numerical simulation of the fracture development zone with a small amount of calculation.

[0020] It should be understood that the contents described in this section are not intended to identify the key or important features of the embodiments of the present invention, nor are they intended to limit the scope of the present invention. Other features of the present invention will become easily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0022] Figure 1 It is a flow chart of a numerical simulation method of a fracture development zone provided by an embodiment of the present invention;

[0023] Figure 2 is a flow chart of another numerical simulation method for fracture development zones provided by an embodiment of the present invention;

[0024] Figure 3 is a flow chart of another numerical simulation method for a fracture development zone provided by an embodiment of the present invention;

[0025] Figure 4 is a schematic diagram of a cross-section T and an orthogonal plane S in a fracture numerical model exemplified in another method for numerical simulation of a fracture development zone provided in an embodiment of the present invention;

[0026] Figure 5 is a schematic diagram of a fracture development zone and a random weakening point in a corresponding fracture numerical model exemplified in another method for numerical simulation of a fracture development zone provided by an embodiment of the present invention;

[0027] Figure 6 is a schematic diagram of a forward modeling profile effect of a fracture numerical model exemplified in another method for numerical simulation of a fracture development zone provided in an embodiment of the present invention;

[0028] Figure 7 It is a structural block diagram of a numerical simulation device for a fracture development zone provided by an embodiment of the present invention;

[0029] Figure 8 It is a schematic diagram of the structure of an electronic device for implementing the numerical simulation method of the fracture development zone according to an embodiment of the present invention. DETAILED DESCRIPTION

[0030] In order to enable those skilled in the art to better understand the scheme of the present invention, the technical scheme in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present invention.

[0031] It should be noted that the terms "first", "second", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. The situations of "target", "original", etc. are similar and will not be repeated here. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units that are clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0032] Before introducing the embodiment of the present invention, the research and development ideas of the embodiment of the present invention are briefly described to better understand why the embodiment of the present invention can ensure the accuracy of numerical simulation of fracture development zones with a smaller amount of calculation.

[0033] In the case of using linear geometric figures to express cracks, in order to accurately express cracks, it is necessary to set up a crack numerical model that matches the crack scale, that is, it is necessary to set up a crack numerical model with a grid size of meters, but this will bring a huge amount of calculation to the numerical simulation process of the crack development zone, which is difficult to implement. Therefore, in order to ensure feasibility, it is possible to consider implementing the numerical simulation process of the crack development zone in a geological body numerical model with a grid size of kilometers, but this will face the difficulty of how to accurately express cracks.

[0034] On this basis, after many studies, it was found that for geological bodies, the development of cracks leads to a comprehensive weakening effect on the characteristics of geological bodies, seismic waves are discrete samples of the crack development zone, and the crack endpoints are numerical singular points in the wave field. Therefore, the embodiment of the present invention adopts a scale equivalence method, regards the small-scale crack development zone as a weakening point presented in the local mean geological body, uses a random probability statistical algorithm to express the comprehensive weakening effect of the crack development zone, and converts the crack endpoints into the contribution of the weakening point of the geological body numerical model, thereby accurately expressing the crack development zone with a small amount of calculation, thereby providing a numerical model with clear scale targeting for subsequent seismic wave field forward simulation.

[0035] Figure 1It is a flow chart of a numerical simulation method of a fracture development zone provided in an embodiment of the present invention. This embodiment can be applied to the case of numerical simulation of fracture development zones, and in particular, can be applied to the case of regarding fracture development zones as weak points in the numerical model of geological bodies, and realizing the numerical simulation process of fracture development zones by using the numerical model of geological bodies. This method can be executed by the numerical simulation device of fracture development zones provided in an embodiment of the present invention, and the device can be implemented by software and / or hardware, and the device can be integrated in an electronic device, and the electronic device can be various user terminals or servers.

[0036] See also Figure 1 The method of the embodiment of the present invention specifically includes the following steps:

[0037] S110, for a geological body numerical model corresponding to a geological body model inlaid with fracture development zones, determining a target quantity distribution obeyed by the number of weakened points on a target surface in the geological body numerical model, wherein the weakened points can be used to characterize the fracture development zones.

[0038] The geological body numerical model can be understood as the result of digitizing the geological body model. In the embodiment of the present invention, the geological body model is inlaid with fracture development zones, and the fracture development zones can be characterized by weakening points (or random weakening points) in the geological body numerical model. Since the geological body numerical model here can be used to realize the numerical simulation of the fracture development zone, it can also be called a fracture numerical model.

[0039] The target plane can be understood as the plane currently being processed in the numerical model of the geological body. For example, it can be any plane in the target direction or a plane that meets the preset conditions. The target direction can be the horizontal direction (i.e., the horizontal direction) or the longitudinal direction, which can be set according to actual conditions and is not specifically limited here.

[0040] The target number distribution can be understood as the distribution obeyed by the number of weakened points on the target surface, that is, the number of weakened points can be represented by random numbers. In combination with the application scenarios that may be involved in the embodiments of the present invention, the type of target number distribution can be a normal distribution. Determine the target number distribution.

[0041] S120. According to the target quantity distribution, a probability distribution of a weak point obeyed by the probability that any point on the target surface is a weak point is obtained, and according to the weak point probability distribution and the target quantity distribution, a total probability that any point on the target surface is a weak point is obtained.

[0042] The probability distribution of weak points can be understood as the distribution of the probability that any point on the target surface is a weak point. In combination with the application scenarios that may be involved in the embodiments of the present invention, the type of the probability distribution of weak points can be a mean distribution. The probability distribution of weak points is determined according to the target quantity distribution.

[0043] There are many points on the target surface, and the total probability can be understood as the probability that any of these points is a weak point. According to the probability distribution of weak points and the distribution of target quantity, the total probability is obtained.

[0044] S130. According to the total probability, the numerical simulation result of the fracture development zone is obtained.

[0045] Among them, since the weak points on the target surface are used to characterize the fracture development zone, the numerical simulation results of the fracture development zone can be obtained according to the total probability. It can be understood that since the grid size corresponding to the numerical model of the geological body does not match the fracture scale, the numerical simulation results here are only numerical simulation results, not numerical solution results, thus ensuring both the amount of calculation and accuracy.

[0046] The technical solution of the embodiment of the present invention is aimed at the geological body numerical model corresponding to the geological body model inlaid with the fracture development zone, by determining the target number distribution obeyed by the number of weakening points on the target surface in the geological body numerical model for characterizing the fracture development zone; then, according to the target number distribution, the probability distribution of weakening points obeyed by the probability that any point on the target surface is a weakening point is obtained, and according to the probability distribution of weakening points and the target number distribution, the total probability that any point on the target surface is a weakening point is obtained; finally, the numerical simulation result of the fracture development zone is obtained according to the total probability. The above technical solution, by considering the small-scale fracture development zone as the weakening point presented in the local mean geological body, adopts the random probability statistical algorithm to express the comprehensive weakening effect of the fracture development zone, and converts the fracture end point into the contribution of the weakening point of the geological body numerical model, thereby ensuring the accuracy of the numerical simulation of the fracture development zone with a small amount of calculation.

[0047] Figure 2 It is a flow chart of another numerical simulation method of fracture development zone provided in an embodiment of the present invention. This embodiment is optimized based on the above-mentioned technical solutions. In this embodiment, optionally, the target surface is any surface in the target direction, and for the vertical direction perpendicular to the target direction, according to the total probability, the numerical simulation result of the fracture development zone is obtained, including: according to the total probability, the vertical quantity distribution obeyed by the number of weakened points on any vertical surface in the vertical direction in the numerical model of the geological body is obtained; according to the target quantity distribution and the vertical quantity distribution, the numerical simulation result of the fracture development zone is obtained. Among them, the explanation of the terms that are the same or corresponding to the above-mentioned embodiments will not be repeated here.

[0048] See also Figure 2 The method of this embodiment may specifically include the following steps:

[0049] S210. For a geological body numerical model corresponding to a geological body model inlaid with fracture development zones, determine a target quantity distribution obeyed by the number of weakened points on any target surface in a target direction in the geological body numerical model, wherein the weakened points can be used to characterize the fracture development zones.

[0050] The target direction can be understood as the direction currently being processed, and the target plane can be understood as any plane in the target direction. For example, when the target direction is a horizontal direction, the target plane can also be called a cross-sectional plane; when the target direction is a longitudinal direction, the target plane can also be called an orthogonal plane.

[0051] S220. According to the target quantity distribution, obtain the probability distribution of weak points to which the probability of any point on the target surface being a weak point obeys, and according to the weak point probability distribution and the target quantity distribution, obtain the total probability that any point on the target surface is a weak point.

[0052] S230. For a vertical direction perpendicular to the target direction, according to the total probability, obtain a vertical quantity distribution obeyed by the number of weak points on any vertical plane in the geological body numerical model located in the vertical direction.

[0053] The vertical direction can be understood as a direction perpendicular to the target direction. The vertical plane can be understood as any plane in the vertical direction. The vertical quantity distribution can be understood as the distribution obeyed by the number of weak points located on the vertical plane in the numerical model of the geological body. In combination with the application scenarios that may be involved in the embodiments of the present invention, the type of vertical quantity distribution can be a binomial distribution. According to the total probability, the vertical quantity distribution is obtained.

[0054] S240. According to the target quantity distribution and the vertical quantity distribution, the numerical simulation results of the fracture development zone are obtained.

[0055] Among them, according to the target quantity distribution in the target direction and the vertical quantity distribution in the vertical direction, the two cooperate with each other to obtain the numerical simulation results of the fracture development zone.

[0056] The technical solution of the embodiment of the present invention obtains the vertical quantity distribution obeyed by the number of weakened points on any vertical plane in the vertical direction in the numerical model of the geological body through total probability, and then combines the target quantity distribution and the vertical quantity distribution, that is, combines the quantity distribution in two directions, to determine the numerical simulation results, thereby further improving the accuracy of the numerical simulation results.

[0057] An optional technical solution, based on the target quantity distribution and the vertical quantity distribution, obtains the numerical simulation results of the fracture development zone, including:

[0058] The standard deviation of the target quantity distribution and the standard deviation of the vertical quantity distribution are taken as the numerical simulation results of the fracture development zone.

[0059] Among them, as mentioned above, since the grid size corresponding to the numerical model of the geological body does not match the fracture scale, the numerical simulation results to be obtained are only numerical simulation results. Therefore, the standard deviation of the target quantity distribution and the standard deviation of the vertical quantity distribution can be used as numerical simulation results, which will help to determine the directional characteristics of the fracture development zone through numerical simulation results in the future.

[0060] On this basis, optionally, the numerical simulation method of the fracture development zone mentioned above also includes:

[0061] According to the numerical relationship between the standard deviation of the target quantity distribution and the standard deviation of the vertical quantity distribution, the directional characteristics of the fracture development zone are obtained.

[0062] Among them, the size of the standard deviation can reflect the difference in the number of weakening points, and in some geological bodies, there are more transverse cracks, while in some geological bodies, there are more longitudinal cracks. Therefore, according to the numerical relationship between the standard deviations of the two quantity distributions, the directional characteristics of the fracture development zone can be obtained, and the directional characteristics can characterize the directionality of the fracture. Usually, the directional characteristics correspond to the direction with a larger standard deviation. For example, here, the target direction is the horizontal direction. When the standard deviation of the target quantity distribution is greater than the standard deviation of the vertical quantity distribution, the directional characteristics are horizontal characteristics, otherwise they are longitudinal characteristics.

[0063] Figure 3 It is a flow chart of another numerical simulation method of a fracture development zone provided in an embodiment of the present invention. This embodiment is optimized based on the above-mentioned technical solutions. In this embodiment, optionally, determining the target quantity distribution obeyed by the number of weakened points on the target surface in the numerical model of the geological body includes: obtaining the fracture porosity of the numerical model of the geological body; for the target quantity distribution to be determined, obtaining the first mean and standard deviation of the target quantity distribution according to the fracture porosity, and determining the target quantity distribution according to the first mean and standard deviation; wherein the target quantity distribution characterizes the distribution obeyed by the number of weakened points on the target surface in the numerical model of the geological body. The explanations of the terms that are the same or corresponding to the above-mentioned embodiments are not repeated here.

[0064] See also Figure 3 The method of this embodiment may specifically include the following steps:

[0065] S310 , for a geological body numerical model corresponding to the geological body model inlaid with fracture development zones, obtain the fracture porosity of the geological body numerical model.

[0066] Fracture porosity can be used to measure the number of fractures in the numerical model of the geological body, which exists because fracture development zones are embedded in the geological body model. Obtain fracture porosity.

[0067] S320, for the target quantity distribution to be determined, obtaining a first mean and a standard deviation of the target quantity distribution according to the fracture porosity, and determining the target quantity distribution according to the first mean and the standard deviation;

[0068] The target quantity distribution represents the distribution obeyed by the number of weak points on the target surface in the numerical model of the geological body, and the weak points can be used to characterize the fracture development zone.

[0069] Among them, a first mean and a standard deviation corresponding to the target quantity distribution are obtained according to the fracture porosity, and then the target quantity distribution is determined according to the first mean and the standard deviation.

[0070] S330. According to the target quantity distribution, a probability distribution of a weak point obeyed by the probability that any point on the target surface is a weak point is obtained, and according to the weak point probability distribution and the target quantity distribution, a total probability that any point on the target surface is a weak point is obtained.

[0071] S340. According to the total probability, the numerical simulation results of the fracture development zone are obtained.

[0072] The technical solution of the embodiment of the present invention reversely infers the first mean and standard deviation through the fracture porosity, and then determines the target quantity distribution according to the first mean and standard deviation, thereby achieving accurate determination of the target quantity distribution.

[0073] An optional technical solution is to obtain the probability distribution of weak points obeyed by the probability that any point on the target surface is a weak point according to the target quantity distribution, including:

[0074] Obtaining a multiple of a second mean value relative to the first mean value, wherein the second mean value represents the mean value of the number of all weak points in the numerical model of the geological body;

[0075] According to the multiple and the target quantity distribution, the probability distribution of the weak point that any point on the target surface is a weak point is obtained.

[0076] The second mean can be understood as the mean of all the weak points in the numerical model of the geological body. The multiple can be understood as the multiple of the second mean relative to the first mean, that is, the result of the second mean / first mean. The multiple is obtained, and then the probability distribution of the weak points is obtained according to the multiple and the target number distribution.

[0077] In order to better understand the above-mentioned technical solutions as a whole, the following is an exemplary description of them in combination with specific examples. For example, for a fracture model data volume with 1000×1000×1000 points and horizontal fracture development, see Figure 4 , Figure 5 and Figure 6, assuming that the number m of weakened points on any cross section T follows a normal distribution with a first mean of ε and a standard deviation of δ:

[0078] m:N(ε,δ 2 )iid

[0079] Among them, ε and δ can be calculated according to the fracture porosity. Here, the fracture porosity is set to 0.01ε%, and the second mean of the number of all weakening points in the fracture model data body is 100ε, that is, the multiple in the above text is 100, so the fracture development intensity in the fracture numerical model can be controlled by ε.

[0080] Furthermore, the probability p(m) that any point on the cross section T is a weakening point follows a uniform distribution:

[0081]

[0082] Furthermore, the total probability P that any point on the cross section T is a weakening point is:

[0083]

[0084] The 10000 is the result of 100×100 in p(m).

[0085] Therefore, the number n of weak points on any orthogonal surface S follows a binomial distribution:

[0086] n:B(n,100×100,P)

[0087] Among them, the mean E(n) = ε, the standard deviation

[0088] The difference in the number of weak points on the cross-section T and the orthogonal plane S (δVSσ(n)) represents the directional characteristics of the crack, that is, the directional characteristics of the crack can be determined by the numerical relationship between δ and σ(n).

[0089] On this basis, the calculation example is as follows:

[0090] For a fracture model data volume with 1000×1000×1000 points and horizontal fracture development, the number m of weakened points on any cross section T is assumed to obey a normal distribution with a first mean of ε=300 and a standard deviation of δ=300. The number of all weakened points in the fracture model data volume has a second mean of 30000 (i.e., 100*300), and the fracture porosity is 3% (i.e., 0.01*300%).

[0091] The total probability that any point on the cross section T is a weakening point is The mean value of the number of weak points on any orthogonal surface S is E(n) = 300, and the standard deviation is

[0092] Therefore, the difference in the number of weakened points on each cross-section (δ=300) is much greater than the difference in the number of weakened points on each orthogonal plane (σ(n)=17.3), that is, the directional characteristics of the fracture development zone are horizontal characteristics.

[0093] Similarly, for longitudinally developed cracks, the same method can be used to establish a crack numerical model.

[0094] The above examples solve the equivalence problem of small-scale fracture numerical models, thus laying a foundation for providing the accuracy of fracture forward simulation response laws, improving the understanding of the seismic response characteristics of fractured oil and gas reservoirs, and providing effective support for the exploration and development of fractured oil and gas reservoirs.

[0095] Figure 7 This is a structural block diagram of a numerical simulation device for fracture development zones provided in an embodiment of the present invention. The device is used to execute the numerical simulation method for fracture development zones provided in any of the above embodiments. The device and the numerical simulation method for fracture development zones in the above embodiments belong to the same inventive concept. For details not fully described in the embodiments of the numerical simulation device for fracture development zones, reference can be made to the embodiments of the numerical simulation method for fracture development zones described above. Figure 7 The device may specifically include: a target quantity distribution determination module 410, a total probability acquisition module 420 and a numerical simulation result acquisition module 430.

[0096] The target quantity distribution determination module 410 is used to determine the target quantity distribution obeyed by the number of weak points on the target surface in the geological body numerical model corresponding to the geological body model inlaid with the fracture development zone, wherein the weak points are used to characterize the fracture development zone;

[0097] A total probability obtaining module 420 is used to obtain a probability distribution of a weak point obeyed by the probability that any point on the target surface is a weak point according to the target quantity distribution, and to obtain a total probability that any point on the target surface is a weak point according to the weak point probability distribution and the target quantity distribution;

[0098] The numerical simulation result obtaining module 430 is used to obtain the numerical simulation result of the fracture development zone according to the total probability.

[0099] Optionally, the target surface is any surface in the target direction, and for a vertical direction perpendicular to the target direction, the numerical simulation result obtaining module 430 may include:

[0100] A vertical quantity distribution obtaining submodule is used to obtain the vertical quantity distribution obeyed by the number of weak points on any vertical surface in the vertical direction in the geological body numerical model according to the total probability;

[0101] The submodule for obtaining numerical simulation results is used to obtain the numerical simulation results of the fracture development zone according to the target quantity distribution and the vertical quantity distribution.

[0102] On this basis, optional, numerical simulation results are obtained from submodules, including:

[0103] The numerical simulation result obtaining unit is used to take the standard deviation of the target quantity distribution and the standard deviation of the vertical quantity distribution as the numerical simulation results of the fracture development zone.

[0104] On this basis, optionally, the numerical simulation device for the fracture development zone further includes:

[0105] The directional characteristic obtaining module is used to obtain the directional characteristics of the fracture development zone according to the numerical relationship between the standard deviation of the target quantity distribution and the standard deviation of the vertical quantity distribution.

[0106] Optionally, based on any of the above devices, the target quantity distribution is a normal distribution; and / or,

[0107] The probability distribution of weak points is a mean distribution; and / or,

[0108] The vertical quantity distribution is the binomial distribution.

[0109] Optionally, the target quantity distribution determination module 410 may include:

[0110] A fracture porosity acquisition unit, used to acquire fracture porosity of a geological body numerical model;

[0111] a target quantity distribution determination unit, for obtaining a first mean and a standard deviation of the target quantity distribution according to the fracture porosity for the target quantity distribution to be determined, and determining the target quantity distribution according to the first mean and the standard deviation;

[0112] The target quantity distribution represents the distribution obeyed by the number of weak points on the target surface in the numerical model of the geological body.

[0113] On this basis, optionally, the total probability obtaining module 420 may include:

[0114] A multiple acquisition unit, used for acquiring a multiple of a second mean value relative to a first mean value, wherein the second mean value represents a mean value of the number of all weak points in the numerical model of the geological body;

[0115] The weakening point probability distribution obtaining unit is used to obtain the weakening point probability distribution obeyed by the probability that any point on the target surface is a weakening point according to the multiple and the target quantity distribution.

[0116] The numerical simulation device of the fracture development zone provided by the embodiment of the present invention determines the target quantity distribution obeyed by the number of weak points on the target surface in the geological body numerical model for characterizing the fracture development zone through a target quantity distribution determination module for the geological body numerical model corresponding to the geological body model inlaid with the fracture development zone; then, through a total probability acquisition module, according to the target quantity distribution, the probability distribution of weak points obeyed by the probability that any point on the target surface is a weak point is obtained, and according to the weak point probability distribution and the target quantity distribution, the total probability that any point on the target surface is a weak point is obtained; finally, through a numerical simulation result acquisition module, the numerical simulation result of the fracture development zone is obtained according to the total probability. The above device regards the small-scale fracture development zone as a weak point presented in the local mean geological body, adopts a random probability statistical algorithm, expresses the comprehensive weakening effect of the fracture development zone, and converts the fracture end point into the contribution of the weak point of the geological body numerical model, thereby ensuring the accuracy of the numerical simulation of the fracture development zone under a small amount of calculation.

[0117] The numerical simulation device for fracture development zones provided in the embodiment of the present invention can execute the numerical simulation method for fracture development zones provided in any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the execution method.

[0118] It is worth noting that in the embodiment of the numerical simulation device for the fracture development zone mentioned above, the various units and modules included are only divided according to functional logic, but are not limited to the above division, as long as the corresponding functions can be achieved; in addition, the specific names of the functional units are only for the convenience of distinguishing each other, and are not used to limit the scope of protection of the present invention.

[0119] Figure 8 A schematic diagram of the structure of an electronic device 10 that can be used to implement an embodiment of the present invention is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present invention described and / or required herein.

[0120] like Figure 8As shown, the electronic device 10 includes at least one processor 11, and a memory connected to the at least one processor 11, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc., wherein the memory stores a computer program that can be executed by at least one processor, and the processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 to the random access memory (RAM) 13. In the RAM 13, various programs and data required for the operation of the electronic device 10 can also be stored. The processor 11, the ROM 12, and the RAM 13 are connected to each other through a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0121] A number of components in the electronic device 10 are connected to the I / O interface 15, including: an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices through a computer network such as the Internet and / or various telecommunication networks.

[0122] The processor 11 may be a variety of general and / or special processing components with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any appropriate processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as a numerical simulation method for a fracture development zone.

[0123] In some embodiments, the numerical simulation method of the fracture development zone may be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as a storage unit 18. In some embodiments, part or all of the computer program may be loaded and / or installed on the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the numerical simulation method of the fracture development zone described above may be performed. Alternatively, in other embodiments, the processor 11 may be configured to execute the numerical simulation method of the fracture development zone by any other appropriate means (e.g., by means of firmware).

[0124] Various implementations of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), systems on chips (SOCs), load programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various implementations can include: being implemented in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0125] Computer programs for implementing the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, so that when the computer program is executed by the processor, the functions / operations specified in the flow chart and / or block diagram are implemented. The computer program may be executed entirely on the machine, partially on the machine, partially on the machine as a stand-alone software package and partially on a remote machine, or entirely on a remote machine or server.

[0126] In the context of the present invention, a computer-readable storage medium may be a tangible medium that may contain or store a computer program for use by or in combination with an instruction execution system, device or equipment. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices or equipment, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. A more specific example of a machine-readable storage medium may include an electrical connection based on one or more lines, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0127] To provide interaction with a user, the systems and techniques described herein may be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and a pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices may also be used to provide interaction with the user; for example, the feedback provided to the user may be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form (including acoustic input, voice input, or tactile input).

[0128] The systems and techniques described herein may be implemented in a computing system that includes backend components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes frontend components (e.g., a user computer with a graphical user interface or a web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such backend components, middleware components, or frontend components. The components of the system may be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0129] A computing system may include a client and a server. The client and the server are generally remote from each other and usually interact through a communication network. The client and server relationship is generated by computer programs running on the corresponding computers and having a client-server relationship with each other. The server may be a cloud server, also known as a cloud computing server or cloud host, which is a host product in the cloud computing service system to solve the defects of difficult management and weak business scalability in traditional physical hosts and VPS services.

[0130] It should be understood that the various forms of processes shown above can be used to reorder, add or delete steps. For example, the steps described in the present invention can be executed in parallel, sequentially or in different orders, as long as the desired results of the technical solution of the present invention can be achieved, and this document does not limit this.

[0131] The above specific implementations do not constitute a limitation on the protection scope of the present invention. It should be understood by those skilled in the art that various modifications, combinations, sub-combinations and substitutions can be made according to design requirements and other factors. Any modification, equivalent substitution and improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims

1. A numerical simulation method for a fracture development zone, characterized in that: include: For a geological body numerical model corresponding to a geological body model inlaid with fracture development zones, determining a target quantity distribution obeyed by the number of weak points on a target surface in the geological body numerical model, wherein the weak points are used to characterize the fracture development zones; According to the target quantity distribution, a probability distribution of a weak point obeyed by the probability that any point on the target surface is the weak point is obtained, and according to the weak point probability distribution and the target quantity distribution, a total probability that any point on the target surface is the weak point is obtained; According to the total probability, the numerical simulation result of the fracture development zone is obtained.

2. The method according to claim 1, characterized in that: The target surface is any surface in the target direction. For a vertical direction perpendicular to the target direction, the numerical simulation result of the fracture development zone is obtained according to the total probability, including: According to the total probability, a vertical quantity distribution obeyed by the number of weak points on any vertical plane in the vertical direction in the numerical model of the geological body is obtained; According to the target quantity distribution and the vertical quantity distribution, the numerical simulation result of the fracture development zone is obtained.

3. The method according to claim 2, characterized in that The step of obtaining the numerical simulation result of the fracture development zone according to the target quantity distribution and the vertical quantity distribution includes: The standard deviation of the target quantity distribution and the standard deviation of the vertical quantity distribution are taken as the numerical simulation results of the fracture development zone.

4. The method according to claim 3, characterized in that Also includes: According to the numerical relationship between the standard deviation of the target quantity distribution and the standard deviation of the vertical quantity distribution, the directional characteristics of the fracture development zone are obtained.

5. The method according to any one of claims 2 to 4, characterized in that: The target quantity distribution is a normal distribution; and / or, The probability distribution of the weak points is a mean distribution; and / or, The vertical quantity distribution is a binomial distribution.

6. The method according to claim 1, characterized in that The step of determining the target quantity distribution to which the number of weak points on the target surface in the geological body numerical model obeys includes: Obtaining fracture porosity of the numerical model of the geological body; For a target quantity distribution to be determined, obtaining a first mean and a standard deviation of the target quantity distribution according to the fracture porosity, and determining the target quantity distribution according to the first mean and the standard deviation; The target quantity distribution represents the distribution obeyed by the quantity of the weak points on the target surface in the numerical model of the geological body.

7. The method according to claim 6, characterized in that The method of obtaining the probability distribution of weak points obeyed by the probability that any point on the target surface is the weak point according to the target quantity distribution includes: Obtaining a multiple of a second mean value relative to the first mean value, wherein the second mean value represents the mean value of the number of all weak points in the numerical model of the geological body; According to the multiple and the target quantity distribution, a probability distribution of weak points is obtained, to which the probability that any point on the target surface is the weak point obeys.

8. A numerical simulation device for fracture development zone, characterized in that: include: A target quantity distribution determination module is used to determine, for a geological body numerical model corresponding to a geological body model inlaid with a fracture development zone, a target quantity distribution obeyed by the number of weak points on a target surface in the geological body numerical model, wherein the weak points are used to characterize the fracture development zone; A total probability obtaining module, used to obtain, according to the target quantity distribution, a probability distribution of a weak point obeyed by the probability that any point on the target surface is the weak point, and to obtain, according to the weak point probability distribution and the target quantity distribution, a total probability that any point on the target surface is the weak point; The numerical simulation result obtaining module is used to obtain the numerical simulation result of the fracture development zone according to the total probability.

9. An electronic device, characterized in that: include: at least one processor; as well as a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor performs the numerical simulation method for the fracture development zone according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the numerical simulation method of the fracture development zone as described in any one of claims 1 to 7 when executed.

Images