CN116401933A - Path simulation method and system applied to multipoint geostatistical modeling - Google Patents
Path simulation method and system applied to multipoint geostatistical modeling Download PDFInfo
- Publication number
- CN116401933A CN116401933A CN202111623704.0A CN202111623704A CN116401933A CN 116401933 A CN116401933 A CN 116401933A CN 202111623704 A CN202111623704 A CN 202111623704A CN 116401933 A CN116401933 A CN 116401933A
- Authority
- CN
- China
- Prior art keywords
- array
- simulation
- work area
- path
- data template
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 81
- 238000000034 method Methods 0.000 title claims abstract description 80
- 238000004458 analytical method Methods 0.000 claims description 13
- 238000010276 construction Methods 0.000 claims description 10
- 238000003491 array Methods 0.000 abstract description 5
- 230000007812 deficiency Effects 0.000 abstract description 3
- 238000012549 training Methods 0.000 description 33
- 230000008569 process Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 239000004576 sand Substances 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 206010063385 Intellectualisation Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/12—Symbolic schematics
Abstract
The invention provides a path simulation method and a path simulation system applied to multipoint geostatistical modeling, wherein the method is characterized in that related arrays are set on the basis of the grid scale of a required work area, a set data module is utilized to traverse grids of the work area, data of the arrays are decided according to a set mode, a new array is built on the basis of the distribution condition of surrounding values of the grids of central points in a data template, element values of the new array are decided according to the element corresponding condition of the array, the sequence of a simulation path is formed, the simulation of a known data enrichment area is realized, then a known data deficiency area is simulated, finally a simulation idea of a non-known data area is simulated, the problem of excessively high simulation result uncertainty in the prior art is solved, the multipoint geostatistical simulation path has an intelligent search function, the uncertainty caused by the random path is reduced in an auxiliary mode on the basis of a multipoint geostatistical data template, and the precision of multipoint geostatistical modeling is improved.
Description
Technical Field
The invention relates to the technical field of geological research in oil and gas field development, in particular to a path simulation method applied to multipoint geostatistical modeling and application.
Background
Reservoir modeling is a core technology for quantitatively describing reservoir development and construction data, and a method for realizing modeling generally comprises two major categories of deterministic modeling and stochastic modeling. The multipoint geostatistical modeling is a leading edge research direction of the current stochastic modeling method, a certain data template (data template) is adopted, an analog grid (work area) and a training image are scanned, a data event (data event) and a training pattern (training pattern) are obtained, and then a partial area or a full area of the data event is replaced according to corresponding probability or similarity, so that the purpose of building a geological model according to the training image is achieved. A pattern vector distance-based multipoint geostatistical modeling method as disclosed in CN108986217a constructs a random path for simulation implementation by aiming at an acquired training pattern library; traversing points to be simulated in the random path, and scanning simulation reality by using a data template to obtain data events at the points to be simulated; taking the training pattern class corresponding to the training pattern prototype with the total angle value closest to the total angle value of the data event as the training pattern class most similar to the data event; obtaining a training pattern most similar to the data event; and then replace the data event with the most similar training pattern.
Therefore, because the simulation path adopted by the current multi-point geostatistical modeling method is a random path, although the modeling effect of the multi-point geostatistical modeling method is improved to a certain extent compared with that of the traditional two-point geostatistical modeling method, the quantity and quality of information contained in each data template on the random path cannot be judged when the geostatistical modeling method is adopted in practical application, so that the simulation path cannot reflect the reservoir analysis process of a geologist, and the underlying scientific knowledge of reservoir characteristics cannot be embodied. Therefore, the current multipoint geostatistical modeling method still has difficulty in reducing uncertainty caused by random simulation paths despite the introduction of data templates and training images.
The information disclosed in the background section of the invention is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person skilled in the art.
Disclosure of Invention
To solve the above problems, the present invention provides a path simulation method applied to multipoint geostatistical modeling, and in one embodiment, the method includes:
step S1, after a grid of a work area to be simulated is obtained, establishing a first array and a second array with equal dimensions according to the grid scale of the work area;
step S2, traversing all grids of the work area by using a data template corresponding to the requirement, and judging the distribution condition of surrounding values of the grids of the central point in the data template;
step S3, assigning a value to the first array according to a set strategy, and further calculating a new third array based on the distribution condition of the surrounding values of the central point grid in the data template;
s4, comparing the center point of the third array element with the corresponding element in the second array, and deciding the value of the element of the third array according to the existence condition of the value of the corresponding element in the second array;
and S5, outputting the elements of the third array according to the sequence from small to large to generate a corresponding target path array, and further simulating a work area according to the sequence of the target path array to perform multipoint geostatistical simulation.
Preferably, in one embodiment, in said step S2,
when traversing all grids of the work area by using the data template, judging whether values exist around the central point in the data template, and if so, recording the known points of the corresponding data template when the current grid is the central point.
Further, in one embodiment, in the step S3, 1 (0, 1) initial values are assigned to the first array using a random number generating function.
Specifically, in one embodiment, in the step S3, the known points recorded when the grid is used as the center point are subtracted from the initial value assigned by the first array to obtain a new third array.
Further, in the step S4, the center point of the third array element is compared with the corresponding element in the second array, and if the corresponding element in the second array has a value, an effective minimum value is assigned to the corresponding element according to a preset rule.
Based on other aspects of the method described in any one or more of the embodiments above, the present invention also provides a storage medium having stored thereon program code that can implement the method described in any one or more of the embodiments above.
Based on the application aspect of the method in any one or more of the above embodiments, the present invention further provides a path simulation system applied to multipoint geostatistical modeling, where the system includes:
the basic array construction module is configured to establish a first array and a second array with equal dimensions according to the grid scale of the work area after the grid of the work area to be simulated is obtained;
the traversal analysis module is configured to traverse all grids of the work area by using a data template corresponding to the requirement, and judge the distribution condition of the surrounding values of the central point grids in the data template;
the logic array construction module is configured to assign a value to the first array according to a set strategy, and further calculate a new third array based on the distribution condition of values around the central point grid in the data template;
the logic array assignment module is configured to compare the center point of the third array element with the corresponding elements in the second array, and to determine the value of the element of the third array according to the existence condition of the value of the corresponding element in the second array;
and the target path determining module is configured to output the elements of the third array in the order from small to large to generate a corresponding target path array, so as to simulate the work area further in the order of the target path array to perform multipoint geostatistical simulation.
Specifically, in one embodiment, the traversal analysis module is configured to:
when traversing all grids of the work area by using the data template, judging whether values exist around the central point in the data template, and if so, recording the known points of the corresponding data template when the current grid is the central point.
Further, in one embodiment, the logic array construction module is specifically configured to: and subtracting the known points recorded when the grid is taken as the center point from the initial value assigned by the first array to obtain a new third array.
Optionally, in one embodiment, the logic array assignment module is specifically configured to:
and comparing the center point of the third array element with the corresponding elements in the second array, and if the corresponding elements in the second array have values, assigning an effective minimum value to the corresponding elements according to a preset principle.
Compared with the closest prior art, the invention has the following beneficial effects:
the method is based on the grid scale setting related array of the required work area, the pattern grid of the work area is traversed by utilizing the data module, a new array is constructed based on the distribution condition of surrounding values of the central point grid in the data template, element values of the new array are decided according to element corresponding conditions of the array, a sequence of a simulation path is formed, the simulation of a known data enrichment area, then a known data deficiency area and finally a simulation idea without the known data area are realized, the problem of excessively high simulation result uncertainty in the prior art is solved, and the method can be effectively applied to the multi-point geostatistical modeling of two-dimensional and three-dimensional patterns of different types of work areas; the multi-point geostatistical simulation path has an intelligent search function, on the basis of a multi-point geostatistical data template, uncertainty caused by a random path is reduced in an auxiliary mode, accuracy of multi-point geostatistical modeling is improved, and geometric shapes and distribution characteristics of the geologic body displayed by a training image are better reproduced.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
The accompanying drawings are included to provide a further understanding of the invention, and are incorporated in and constitute a part of this specification, illustrate the invention and together with the embodiments of the invention, serve to explain the invention, without limitation to the invention. In the drawings:
FIG. 1 is a flow chart of a path simulation method for multi-point geostatistical modeling according to one embodiment of the present invention;
FIG. 2 is a schematic flow diagram of a path modeled in accordance with an embodiment of the present invention as applied to multi-point geostatistical modeling;
FIGS. 3 (a) and 3 (b) are a river phase training image and a well data schematic diagram, respectively, of a two-dimensional example;
FIGS. 4 (a), 4 (b) and 4 (c) are graphs showing the results obtained by simulation of two-dimensional examples based on sequential paths, random paths and paths obtained by embodiments of the present invention, respectively;
FIG. 5 (a) is an example of a training image and grid plot of a three-dimensional instance of historical multi-point geostatistical modeling;
FIGS. 5 (b), 5 (c) and 5 (d) are three-dimensional and raster pattern comparison results of three-dimensional example simulations based on sequential paths, random paths and paths of the present invention, respectively;
FIG. 6 is a schematic diagram of a path modeling system for multi-point geostatistical modeling provided by an embodiment of the present invention.
Detailed Description
The following will explain the embodiments of the present invention in detail with reference to the drawings and examples, so that the practitioner of the present invention can fully understand how to apply the technical means to solve the technical problems, achieve the implementation process of the technical effects, and implement the present invention according to the implementation process. It should be noted that, as long as no conflict is formed, each embodiment of the present invention and each feature of each embodiment may be combined with each other, and the formed technical solutions are all within the protection scope of the present invention.
Although a flowchart depicts operations as a sequential process, many of the operations can be performed in parallel, concurrently, or at the same time. The order of the operations may be rearranged. The process may be terminated when its operations are completed, but may have additional steps not included in the figures. The processes may correspond to methods, functions, procedures, subroutines, and the like.
The terms "first," "second," and the like may be used herein to describe various elements, but these elements should not be limited by these terms, and these terms are used merely to distinguish one element from another. The term "and/or" as used herein includes any and all combinations of one or more of the associated listed items. When an element is referred to as being "connected" or "coupled" to another element, it can be directly connected or coupled to the other element or intervening elements may be present.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises" and/or "comprising," when used herein, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
Reservoir modeling is a core technology for quantitative reservoir description and can be divided into two methods, namely deterministic modeling and stochastic modeling. The multipoint geostatistical modeling is a leading edge research direction of the current stochastic modeling method, and adopts a certain data template (data template), scans an analog grid (work area) and a training image to obtain a data event (data event) and a training pattern (training pattern), and replaces a partial area or a whole area of the data event according to corresponding probability or similarity, so that the purpose of building a geologic model according to the training image is achieved.
In the multi-point geostatistical phase modeling, the determination of the value of a point to be estimated, the reduction of the uncertainty of modeling, and the success of the multi-point geostatistical modeling method are always key. The multipoint geostatistical modeling method can be divided into a multipoint geostatistical modeling method based on probability and a multipoint geostatistical modeling method based on similarity according to the modeling value basis, such as a multipoint geostatistical modeling method based on pattern vector distance disclosed in CN108986217A, and a random path realized by simulation is constructed by aiming at an acquired training pattern library; traversing points to be simulated in the random path, and scanning simulation reality by using a data template to obtain data events at the points to be simulated; taking the training pattern class corresponding to the training pattern prototype with the total angle value closest to the total angle value of the data event as the training pattern class most similar to the data event; obtaining a training pattern most similar to the data event; and then replace the data event with the most similar training pattern.
Therefore, the modeling effect of the multi-point geostatistical modeling method is improved to a certain extent compared with that of the traditional two-point geostatistical modeling method, but because the simulation path adopted by the current multi-point geostatistical modeling method is a random path, the quantity and quality of information contained in each data template on the random path cannot be judged, so that the simulation path cannot reflect the reservoir analysis process of a developed geologist, and the reservoir analysis process of the developed geologist implies a methodology for scientifically knowing the reservoir characteristics. Therefore, the current multipoint geostatistical modeling method still has difficulty in reducing uncertainty caused by random simulation paths despite the introduction of data templates and training images.
In order to solve the problems, the invention provides a path simulation method and application applied to multipoint geostatistical modeling, and the invention optimizes a simulation path on the basis of a framework of multipoint geostatistical modeling based on similarity, so that the multipoint geostatistical simulation path has an intelligent search function, a simulation path capable of simulating the reservoir analysis thought and process of a geologist is formed, uncertainty brought by a random simulation path to the multipoint geostatistical modeling is reduced, and the precision of the multipoint geostatistical modeling is improved.
The detailed flow of the method of embodiments of the present invention is described in detail below based on the attached drawing figures, where the steps shown in the flowchart of the figures may be performed in a computer system containing, for example, a set of computer executable instructions. Although a logical order of steps is depicted in the flowchart, in some cases the steps shown or described may be performed in a different order than presented.
Example 1
Fig. 1 is a schematic flow chart of a path simulation method for multipoint geostatistical modeling according to an embodiment of the present invention, and as can be seen with reference to fig. 1, the method includes the following steps.
Step S110, after the grids of the work area to be simulated are obtained, establishing a first array and a second array with equal dimensions according to the grid scale of the work area;
step S120, traversing all grids of the work area by using a data template corresponding to the requirement, and judging the distribution condition of surrounding values of the grids of the central point in the data template;
step S130, assigning a value to the first array according to a set strategy, and further calculating a new third array based on the distribution condition of the surrounding values of the central point grid in the data template;
step S140, comparing the center point of the third array element with the corresponding elements in the second array, and deciding the value of the third array element according to the existence of the value of the corresponding element in the second array;
and step S150, outputting the elements of the third array in the order from small to large to generate a corresponding target path array, and further simulating the work area according to the order of the target path array to perform multipoint geostatistical simulation.
The simulation path formed based on the embodiment can be used in a full-type similarity-based multi-point geostatistical modeling algorithm, and when the multi-point geostatistical modeling algorithm generates a random path, the simulation path is applied to access the simulation path, and then the multi-point geostatistical simulation is performed according to the simulation path of the invention, as shown in fig. 2.
The simulation path obtained by the invention realizes multipoint geostatistical simulation, improves the random path of multipoint geostatistical modeling, can push the multipoint geostatistical modeling to traverse all the non-simulated grids, and can also enable the multipoint geostatistical modeling to be carried out according to the set sequence when traversing the non-simulated grids. The principle of the invention is to simulate the thought of a geologist for analyzing a reservoir, so that the multipoint geostatistical random modeling process can simulate the thought and the method of the geologist for analyzing the reservoir, and according to the concentration degree of known data, the simulation sequence of non-simulation points is arranged according to a descending order, so as to realize the intellectualization of a simulation path.
In practical application, in step S110, 1 array (denoted as Sim) with dimensions equal to the grid and 1 array (denoted as SimFa) for storing known data are established according to the grid scale of the work area to be simulated; in the execution process, the work area can be divided into grid blocks in a three-dimensional space to form a three-dimensional grid body formed by stacking the grid blocks. And setting the dimensions of the grid in three directions XYZ, and equally dividing the dimensions of the work area XYZ in three directions to establish a three-dimensional grid body. If there is well-traversed and well-log interpreted phase data in the grid, the value of the grid is a known value. If the heterogeneity of the reservoir in the work area is strong, the size of the grid XYZ in three directions can be properly reduced;
the steps of building Sim and simFa arrays are the same as the steps of building the work area grids, but Sim arrays do not require preset values, simFa arrays require first identifying the grids traversed by the well track, and then setting the values of these grids to known values.
Further, in the step S120, when traversing all grids of the work area with the data template, it is determined whether there is a value around the center point inside the data template, and if there is a value, the known points of the data template corresponding to the current grid being the center point are recorded.
When the method is actually applied, all grids of a work area are traversed by using a set data template, whether values exist around the center point in the data template or not is judged, and the known points of the data template corresponding to the center point are recorded; the data template can be a three-dimensional grid body with m x n taking u as a center point, m and n are generally odd numbers, when traversing the grid of the work area, judging whether the grid in the range of the template has a value by adopting a judgment statement, if so, counting, and then assigning a value to the center point
In an alternative embodiment, in the step S130, 1 (0, 1) initial values are assigned to the first array using a random number generating function.
Specifically, in one embodiment, in the step S130, the known points recorded when the grid is used as the center point are further subtracted from the initial value assigned by the first array to obtain a new third array.
In actual application, a random number generating function is adopted to assign 1 (0, 1) initial value to the Sim array, and the initial value is used for subtracting the known point number recorded in the step 2 by taking the grid as a central point to obtain 1 new array (marked as newSim);
the researchers consider that since Sim is a random value between (0, 1), after counting the grids that are central points, if there is a value, it is at least an integer value of 1 or more, and the two are subtracted, the more the number of known points in the data template is, the smaller the value is, and the earlier the ordering in the order path array is, the more the simulation is performed, so that the region with more known information is simulated first.
Further, in the step S140, the center point of the third array element is compared with the corresponding element in the second array, and if the corresponding element in the second array has a value, an effective minimum value is assigned to the corresponding element according to a preset rule.
In practical application, comparing the element center point of the newSim array with the corresponding element in the SimFa array, and if the corresponding element in the SimFa array has a value, assigning an effective minimum value (for example, -999); in the step, assigning a value to a grid corresponding to a value in the SimFa array in the newSim array, wherein the value is only a minimum value for distinguishing from a small value existing in the step, so that a computer program can accurately judge that the points are known points so as to directly skip in a subsequent path array and not simulate again;
further, each array element is output in the order of newSim array elements from small to large, and 1 path array (denoted as order) is generated as a target path array; based on this, the simulation work area is simulated geostatistically in the order of the resulting target path array (except for 999 elemental points). Based on the logic of the above embodiment, it is possible to order the simulation areas according to the degree of data concentration, to simulate the areas with known data enrichment (more wells), to simulate the areas with known data deficiency (less wells), and to simulate the areas without known data (no wells).
The simulation path obtained by the scheme of the invention improves the random path of the multipoint geostatistical modeling, so that the multipoint geostatistical random modeling process can simulate the thought and method of reservoir analysis by geologist development, firstly simulate the area with rich data (more wells), then simulate the area with insufficient data (less wells), finally simulate the area without data (no wells), and enable the multipoint geostatistical simulation path to have an intelligent search function, and improve the precision of the multipoint geostatistical modeling on the basis of a multipoint geostatistical data template.
Case example:
1. two-dimensional model instance
Establishing a training image TI (grid scale: 100 x 1 x 45), establishing a simulation grid (the same as the training image), wherein 7 wells exist in a simulation work area, the number of set condition points is 2, and performing conditional simulation on the simulation work area, as shown in fig. 3, fig. 3 (a) is a river phase training image of a two-dimensional example, wherein a dark color part is river sand, light gray is mud, and the river sand is a stacked structure of top and bottom protrusions; FIG. 3 (b) is well data, dark colored river sand, light gray mud;
fig. 4 (a) is a model implementation obtained based on sequential path simulation in a two-dimensional example, fig. 4 (b) is a model implementation obtained based on random path simulation, and fig. 4 (c) is a model implementation obtained based on path simulation of the present invention, and it can be seen that the model obtained by path simulation of the present invention can embody potential information of training images more closely to knowledge of geologist.
2. Three-dimensional model instance
A training image TI (mesh size: 155×94×40) is created, a simulation mesh (same as the training image) is created, 32 wells exist in the simulation work area, the number of condition points is set to 3, and the simulation work area is subjected to conditional simulation. A training image and a raster image of a three-dimensional example are shown in fig. 5 (a). The three-dimensional example is shown in fig. 5 (b), 5 (c) and 5 (d), respectively, based on a three-dimensional diagram and a raster diagram of sequential path, random path and path simulation of the present invention.
Based on the simulation results of the two-dimensional and three-dimensional examples, the invention can effectively show that the thought and the method for carrying out reservoir analysis by a geology development geologist are included in the multi-point geostatistical modeling process by optimizing the simulation path, thus reducing the uncertainty in the modeling process, improving the modeling precision, better reproducing the geometric shape and the distribution characteristics of the geologic body displayed by the training image, and being feasible and reliable to carry out phase modeling by applying the method.
For the foregoing method embodiments, for simplicity of explanation, the methodologies are shown as a series of acts, but one of ordinary skill in the art will appreciate that the present invention is not limited by the order of acts, as some steps may, in accordance with the present invention, occur in other orders or concurrently. Further, those skilled in the art will also appreciate that the embodiments described in the specification are all preferred embodiments, and that the acts and modules referred to are not necessarily required for the present invention.
It should be noted that in other embodiments of the present invention, the method may also be combined with any one or some of the above embodiments to obtain a new path simulation method, so as to implement optimization and support for the multi-point geostatistical modeling technique.
It should be noted that, based on the method in any one or more of the foregoing embodiments of the present invention, the present invention further provides a storage medium, where a program code is stored, where the program code can implement the method in any one or more of the foregoing embodiments, and when executed by an operating system, can implement the path simulation method applied to multi-point geostatistical modeling as described above.
Example two
The method is described in detail in the embodiments disclosed in the present application, and the method of the present application may be implemented by using various types of devices or systems, so based on other aspects of the method described in any one or more of the embodiments, the present application also provides a path simulation system applied to multi-point geostatistical modeling, where the system is used to perform the path simulation method applied to multi-point geostatistical modeling described in any one or more of the embodiments. Specific examples are given below for details.
Specifically, fig. 6 shows a schematic structural diagram of a path simulation system applied to multipoint geostatistical modeling according to an embodiment of the present invention, and as shown in fig. 6, the system includes:
the basic array construction module is configured to establish a first array and a second array with equal dimensions according to the grid scale of the work area after the grid of the work area to be simulated is obtained;
the traversal analysis module is configured to traverse all grids of the work area by using a data template corresponding to the requirement, and judge the distribution condition of the surrounding values of the central point grids in the data template;
the logic array construction module is configured to assign a value to the first array according to a set strategy, and further calculate a new third array based on the distribution condition of values around the central point grid in the data template, and the new third array is used as a logic array of a decision simulation sequence;
the logic array assignment module is configured to compare the center point of the third array element with the corresponding elements in the second array, and to determine the value of the element of the third array according to the existence condition of the value of the corresponding element in the second array;
and the target path determining module is configured to output the elements of the third array in the order from small to large to generate a corresponding target path array, so as to simulate the work area further in the order of the target path array to perform multipoint geostatistical simulation.
Further, in one embodiment, the traversal analysis module is configured to:
when traversing all the non-simulated grids in the work area by using the data template, judging whether values exist around the central point in the data template, and if so, recording the known points of the corresponding data template when the current grid is the central point.
Specifically, in one embodiment, the logic array construction module is specifically configured to: and subtracting the known points recorded when the grid is taken as the center point from the initial value assigned by the first array to obtain a new third array.
In an alternative embodiment, the logic array assignment module specifically assigns 1 (0, 1) initial values to the first array using a random number generation function.
Further, in one embodiment, the logic array assignment module is specifically configured to:
and comparing the center point of the third array element with the corresponding elements in the second array, and if the corresponding elements in the second array have values, assigning an effective minimum value to the corresponding elements according to a preset principle.
The embodiment of the invention provides a path simulation system applied to multipoint geostatistical modeling, and each module or unit structure can independently or in combination operate according to actual analysis requirements and calculation requirements so as to realize corresponding technical effects.
It is to be understood that the disclosed embodiments are not limited to the specific structures, process steps, or materials disclosed herein, but are intended to extend to equivalents of these features as would be understood by one of ordinary skill in the relevant arts. It is also to be understood that the terminology used herein is for the purpose of describing particular embodiments only, and is not intended to be limiting.
Reference in the specification to "one embodiment" or "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment of the invention. Thus, the appearances of the phrase "one embodiment" in various places throughout this specification are not necessarily all referring to the same embodiment.
Although the embodiments of the present invention are described above, the embodiments are only used for facilitating understanding of the present invention, and are not intended to limit the present invention. Any person skilled in the art can make any modification and variation in form and detail without departing from the spirit and scope of the present disclosure, but the scope of the present disclosure is still subject to the scope of the appended claims.
Claims (10)
1. A path simulation method applied to multi-point geostatistical modeling, the method being for multi-point geostatistical modeling, the method comprising:
step S1, after a grid of a work area to be simulated is obtained, establishing a first array and a second array with equal dimensions according to the grid scale of the work area;
step S2, traversing all grids of the work area by using a data template corresponding to the requirement, and judging the distribution condition of surrounding values of the grids of the central point in the data template;
step S3, assigning a value to the first array according to a set strategy, and further calculating a new third array based on the distribution condition of the surrounding values of the central point grid in the data template;
s4, comparing the center point of the third array element with the corresponding element in the second array, and deciding the value of the element of the third array according to the existence condition of the value of the corresponding element in the second array;
and S5, outputting the elements of the third array according to the sequence from small to large to generate a corresponding target path array, and further simulating a work area according to the sequence of the target path array to perform multipoint geostatistical simulation.
2. The method according to claim 1, wherein, in said step S2,
when traversing all grids of the work area by using the data template, judging whether values exist around the central point in the data template, and if so, recording the known points of the corresponding data template when the current grid is the central point.
3. The method according to claim 1, wherein in the step S3, 1 (0, 1) initial values are assigned to the first array by using a random number generating function.
4. The method according to claim 1, wherein in the step S3, the known points recorded when the grid is centered are subtracted from the initial value assigned by the first array to obtain a new third array.
5. The method according to claim 1, wherein in the step S4, the center point of the third array element is compared with the corresponding element in the second array, and if the corresponding element in the second array has a value, a valid minimum value is assigned to the corresponding element according to a predetermined rule.
6. A storage medium having stored thereon program code for implementing the method of any of claims 1 to 5.
7. A path simulation system for application in multipoint geostatistical modeling, characterized in that the system is adapted to perform the method as claimed in claims 1-5, the system comprising:
the basic array construction module is configured to establish a first array and a second array with equal dimensions according to the grid scale of the work area after the grid of the work area to be simulated is obtained;
the traversal analysis module is configured to traverse all grids of the work area by using a data template corresponding to the requirement, and judge the distribution condition of the surrounding values of the central point grids in the data template;
the logic array construction module is configured to assign a value to the first array according to a set strategy, and further calculate a new third array based on the distribution condition of values around the central point grid in the data template;
the logic array assignment module is configured to compare the center point of the third array element with the corresponding elements in the second array, and to determine the value of the element of the third array according to the existence condition of the value of the corresponding element in the second array;
and the target path determining module is configured to output the elements of the third array in the order from small to large to generate a corresponding target path array, so as to simulate the work area further in the order of the target path array to perform multipoint geostatistical simulation.
8. The system of claim 7, wherein the traversal analysis module is configured to:
when traversing all grids of the work area by using the data template, judging whether values exist around the central point in the data template, and if so, recording the known points of the corresponding data template when the current grid is the central point.
9. The system of claim 7, wherein the logic array construction module is specifically configured to: and subtracting the known points recorded when the grid is taken as the center point from the initial value assigned by the first array to obtain a new third array.
10. The system of claim 7, wherein the logic array assignment module is specifically configured to:
and comparing the center point of the third array element with the corresponding elements in the second array, and if the corresponding elements in the second array have values, assigning an effective minimum value to the corresponding elements according to a preset principle.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111623704.0A CN116401933A (en) | 2021-12-28 | 2021-12-28 | Path simulation method and system applied to multipoint geostatistical modeling |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111623704.0A CN116401933A (en) | 2021-12-28 | 2021-12-28 | Path simulation method and system applied to multipoint geostatistical modeling |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116401933A true CN116401933A (en) | 2023-07-07 |
Family
ID=87018460
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111623704.0A Pending CN116401933A (en) | 2021-12-28 | 2021-12-28 | Path simulation method and system applied to multipoint geostatistical modeling |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116401933A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117291138A (en) * | 2023-11-22 | 2023-12-26 | 全芯智造技术有限公司 | Method, apparatus and medium for generating layout elements |
-
2021
- 2021-12-28 CN CN202111623704.0A patent/CN116401933A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117291138A (en) * | 2023-11-22 | 2023-12-26 | 全芯智造技术有限公司 | Method, apparatus and medium for generating layout elements |
| CN117291138B (en) * | 2023-11-22 | 2024-02-13 | 全芯智造技术有限公司 | Method, apparatus and medium for generating layout elements |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10467357B2 (en) | Geobody continuity in geological models based on multiple point statistics | |
| Zhang | Filter-based training pattern classification for spatial pattern simulation | |
| US9140821B2 (en) | Ordered multipoint geostatistics simulation using non-symmetric search mask | |
| CA2875496C (en) | System and method for optimizing the number of conditioning data in multiple point statistics simulation | |
| WO2009138290A2 (en) | Multi-point reservoir modelling | |
| WO2010057505A1 (en) | A deterministic version of the multiple point geostatistics simulation / reconstruction method wxth. the simulated / reconstructed values are directly taken from the training images without prior estimation of the conditional | |
| EP2880470A2 (en) | Hybrid method of combining multipoint statistic and object-based methods for creating reservoir property models | |
| CN107204040B (en) | Multipoint geostatistical modeling method and device and computer storage medium | |
| EP2880469A2 (en) | Method for editing a multi-point facies simulation | |
| CN104504754B (en) | A kind of method and device of oil and gas reservoir multi-point statistic modeling | |
| CN114092697B (en) | Building facade semantic segmentation method with attention fused with global and local depth features | |
| CN116401933A (en) | Path simulation method and system applied to multipoint geostatistical modeling | |
| US8422745B2 (en) | Method of modelling an underground heterogeneous medium from multipoint statistics | |
| CN110927793B (en) | Reservoir prediction method and system based on sequential random fuzzy simulation | |
| US20150317419A1 (en) | Local direct sampling method for conditioning an existing reservoir model | |
| EP3320450B1 (en) | Improved geobody continuity in geological models based on multiple point statistics | |
| CN116226623B (en) | Mark layer division method and device based on SegNet segmentation model and computer equipment | |
| CN113408109B (en) | Hybrid parallel method for multipoint geostatistical random simulation process | |
| CN114491898A (en) | Multi-point geostatistical modeling method and system | |
| CN116484698A (en) | Method for predicting soil pollutant distribution of small sample by fusing data enhancement and neural network | |
| Arpat | SIMPAT: stochastic simulation with patterns | |
| CN117973199A (en) | Method and device for modeling generation countermeasure network sedimentary facies based on single training image | |
| Izzo | Exploring the possible role of parametric design in addressing the compositional problem: theory and practice of an evolving approach toward the urban design perspective. | |
| Li et al. | A New Multiple Grid method in SIMPAT: a proposal | |
| Dubrule et al. | [5] 3 Geostatistics: A Way Towards Data Integration and Uncertainty Quantification |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |