Geological modeling method and system for simulating central beach based on riverway distribution characteristics

Abstract

The invention discloses a geological modeling method and system for simulating a central beach based on river distribution characteristics, and the method comprises the following steps: transmitting the development characteristic data of a braided river reservoir in a research area to PETREL modeling software, and building a braided river reservoir river microfacies model according to the development characteristic data; obtaining a water flow azimuth angle of each grid in the braided river reservoir river channel microfacies model, and establishing an azimuth angle three-dimensional model body of the braided river reservoir river channel microfacies; taking the azimuth angle three-dimensional model body as a modeling constraint condition, setting a development rule for simulating a central beach microfacies, and establishing a central beach model in the braided river reservoir channel microfacies model according to the development characteristic data; therefore, the central beach position and the central beach distribution direction can be accurately simulated.

Description

Technical Field

[0001] The invention relates to the technical field of oil reservoir description, and in particular to a geological modeling method and system for simulating a central shoal based on river channel distribution characteristics. Background Art

[0002] Braided river is an important type of continental sedimentary system. Due to the strong heterogeneity of its reservoir, the spatial superposition pattern and distribution law of sand bodies are very complex. It is still difficult to establish a three-dimensional geological model with high reliability for braided river reservoirs. Configuration element analysis is an important research technical means for reservoir characterization and a hot topic in current research. Many scholars have conducted in-depth research on the architecture of braided river reservoirs and have done a lot of work on braided river reservoir modeling. Among them, the channel configuration unit and the mid-channel bar configuration unit are the most important configuration units of braided river reservoirs. The mid-channel configuration unit has stronger homogeneity and its reservoir quality is better than that of the channel configuration unit. The quantitative study of its distribution characteristics is of great significance.

[0003] Commonly used phase modeling methods include target-based simulation, multi-point geostatistics, and sequential indicator simulation. Through a detailed comparison of the basic theories and algorithmic processes of several modeling methods, each modeling method has its own advantages and disadvantages, and should be reasonably selected according to different situations and actual needs in combination with the applicable conditions, advantages and disadvantages of each method. Among them, the target-based simulation takes the target body as the object, and on the basis of the prior geometric shape of the geological body, randomly simulates the position of the center point of the geological body, assigns the geometric shape and other attributes of the geological body at this point, and finally establishes a three-dimensional model of the spatial distribution of the geological body. However, in the process of simulating the heart bar, the existing target-based simulation method is easy to simulate the heart bar position in the floodplain, resulting in the simulation results being seriously inconsistent with the actual geological development law. Summary of the invention

[0004] The present invention provides a geological modeling method and system for simulating a mid-bar based on river channel distribution characteristics, which can accurately simulate the mid-bar position and its distribution direction.

[0005] In a first aspect, a geological modeling method for simulating a central shoal based on river channel distribution characteristics is provided, comprising the following steps:

[0006] The development characteristic data of the braided river reservoir in the study area were transferred to the PETREL modeling software, and a channel microfacies model of the braided river reservoir was established based on the development characteristic data;

[0007] Obtaining the water flow azimuth of each grid in the braided river reservoir channel microfacies model, and establishing a three-dimensional azimuth model of the braided river reservoir channel microfacies;

[0008] The azimuth three-dimensional model body is used as a modeling constraint condition, and the development rules for simulating the center-bar microfacies are set, and a center-bar model is established in the braided river reservoir channel microfacies model according to the development characteristic data.

[0009] According to the first aspect, in a first possible implementation of the first aspect, the development characteristic data of the braided river reservoir in the study area includes the scope of the study area, well spacing, river channel geological data and mid-shoal geological data;

[0010] The river geological data include: river direction, river width, river thickness, river amplitude, river wavelength and river ratio;

[0011] The core-bar geological data include: core-bar shape, core-bar width, core-bar thickness, core-bar length and core-bar ratio.

[0012] According to the first aspect, in a second possible implementation of the first aspect, the step of “establishing a braided river reservoir channel microfacies model according to the development characteristic data” specifically includes the following steps:

[0013] Building a three-dimensional network grid, and setting the number of grids and the size of a single grid of the three-dimensional network grid;

[0014] Taking floodplain mud as the channel modeling background, based on the Fluvsim algorithm and according to the channel geological data in the development characteristic data, a braided river reservoir channel microfacies model including floodplain mud microfacies is established on the three-dimensional network framework.

[0015] According to the first aspect, in a third possible implementation of the first aspect, the step of “obtaining the water flow azimuth of each grid in the braided river reservoir channel microfacies model and establishing an azimuth three-dimensional model of the braided river reservoir channel microfacies” specifically includes the following steps:

[0016] Obtaining the hydrodynamic axes of different single-layer river channels in the vertical direction and the boundary lines of the river channels on both sides in the braided river reservoir channel microfacies model;

[0017] Obtain the water flow azimuth corresponding to each grid passing through the water flow dynamic axis and the boundary lines of the river channels on both sides in each single-layer river channel;

[0018] According to the water flow azimuth corresponding to each grid passing through the water flow dynamic axis and the river channel boundary lines on both sides, the Kriging interpolation algorithm is used to calculate the water flow azimuth corresponding to each grid between the water flow dynamic axis and the river channel boundary lines on both sides;

[0019] According to the water flow azimuth corresponding to each grid in different single-layer river channels, a three-dimensional azimuth model of the braided river reservoir channel microfacies is established.

[0020] According to the first aspect, in a fourth possible implementation of the first aspect, the step of "taking the azimuth three-dimensional model body as a modeling constraint condition, setting the development rules of the simulated center-bar microfacies, and establishing a center-bar model in the braided river reservoir channel microfacies model according to the development characteristic data" specifically includes the following steps:

[0021] The channel microfacies model of the braided river reservoir is used as the background for the barbed-bar modeling, and the development rules for simulating the barbed-bar microfacies are set;

[0022] The azimuth three-dimensional model body is used to constrain the distribution direction of the bar, and based on the Fluvsim algorithm and the bar geological data in the development characteristic data, a bar model including bar microfacies, river channel microfacies and floodplain mud microfacies is established.

[0023] According to the first aspect, in a fifth possible implementation of the first aspect, the development rule of the center-bar microfacies is: searching for channel microfacies in the braided river reservoir channel microfacies model including the channel and floodplain mud, simulating the center-bar microfacies inside the searched channel microfacies, and the simulated center-bar microfacies does not contact the floodplain mud.

[0024] In a second aspect, a geological modeling system for simulating a central shoal based on river channel distribution characteristics is provided, comprising:

[0025] The channel facies modeling module is used to transmit the development characteristic data of the braided river reservoir in the study area to the PETREL modeling software, and establish the channel microfacies model of the braided river reservoir based on the development characteristic data;

[0026] A water flow azimuth module is connected to the channel phase modeling module for obtaining the water flow azimuth of each grid in the braided river reservoir channel microfacies model and establishing a three-dimensional azimuth model of the braided river reservoir channel microfacies; and

[0027] The center-bar modeling module is communicatively connected with the water flow azimuth module, and is used to use the azimuth three-dimensional model body as a modeling constraint condition, set the development rules of the simulated center-bar microfacies, and establish a center-bar model in the braided river reservoir channel microfacies model according to the development characteristic data.

[0028] In a third aspect, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the geological modeling method for simulating the mid-bar based on the distribution characteristics of the river channel as described above is implemented.

[0029] In a fourth aspect, an electronic device is provided, comprising a storage medium, a processor, and a computer program stored in the storage medium and executable on the processor, wherein when the processor runs the computer program, the method for simulating the central shoal based on the distribution characteristics of the river channel as described above is implemented as a geological modeling method.

[0030] Compared with the prior art, the advantages of the present invention are as follows: by obtaining the water flow azimuth of each grid in the braided river reservoir channel microfacies model, an azimuth three-dimensional model of the braided river reservoir channel microfacies is established, and then the distribution direction of the center bar is constrained by the azimuth model. At the same time, the center bar microfacies development rules are set, and the center bar model is established to characterize the real distribution characteristics of the microfacies inside the river channel. The method described in the present invention provides a reference method for establishing a braided river reservoir geological model, and provides a reliable three-dimensional geological model for accurately simulating the center bar position and its distribution direction, as well as for later finding favorable areas. BRIEF DESCRIPTION OF THE DRAWINGS

[0031] Figure 1 It is a schematic flow chart of a geological modeling method for simulating a central shoal based on river channel distribution characteristics of the present invention;

[0032] Figure 2 A schematic diagram of a three-dimensional network grid established for the present invention;

[0033] Figure 3 The channel microfacies model of the braided river reservoir established for the present invention and its microfacies ratio histogram;

[0034] Figure 4 This is a schematic diagram of the water flow dynamic axis obtained by the present invention based on the river channel microfacies model at every 2.5m vertical slice;

[0035] Figure 5 It is a schematic diagram of obtaining different grid azimuth data based on a single-period river channel in the present invention;

[0036] Figure 6 It is a schematic diagram of the change characteristics of the river section and the azimuth angle data of each grid on the plane established based on the water flow dynamic axis of the present invention;

[0037] Figure 7 It is a schematic diagram of the azimuth three-dimensional model body about the direction of the river channel established by the present invention based on the water flow dynamic axis;

[0038] Figure 8 The invention establishes a heart beach model and a heart beach hollowing model based on direction and rule constraints. DETAILED DESCRIPTION

[0039] Reference will now be made in detail to specific embodiments of the present invention, examples of which are illustrated in the accompanying drawings. Although the present invention will be described in conjunction with specific embodiments, it will be understood that it is not intended to limit the present invention to the embodiments described. On the contrary, it is intended to cover changes, modifications and equivalents included in the spirit and scope of the present invention as defined by the appended claims. It should be noted that the method steps described herein can all be implemented by any functional block or functional arrangement, and any functional block or functional arrangement can be implemented as a physical entity or a logical entity, or a combination of the two.

[0040] In order to enable those skilled in the art to better understand the present invention, the present invention is further described in detail below in conjunction with the accompanying drawings and specific implementation methods.

[0041] Note: The example to be introduced below is only a specific example, and is not intended to limit the embodiments of the present invention to the following specific steps, values, conditions, data, sequences, etc. Those skilled in the art can apply the concept of the present invention to construct more embodiments not mentioned in this specification by reading this specification.

[0042] See also Figure 1 As shown, an embodiment of the present invention provides a geological modeling method for simulating a central shoal based on river channel distribution characteristics, comprising the following steps:

[0043] S100, transmitting the development characteristic data of the braided river reservoir in the study area to the PETREL modeling software, and establishing a channel microfacies model of the braided river reservoir according to the development characteristic data;

[0044] S200, obtaining the water flow azimuth of each grid in the braided river reservoir channel microfacies model, and establishing a three-dimensional azimuth model of the braided river reservoir channel microfacies;

[0045] S300, using the azimuth three-dimensional model body as a modeling constraint condition, and setting a development rule for simulating the center-bar microfacies, and establishing a center-bar model in the braided river reservoir channel microfacies model according to the development characteristic data;

[0046] Specifically, in this embodiment, the Fluvsim algorithm is a simulation method proposed by Deutsch for building complex river channels, also known as the "target-based simulation method". This modeling method uses the target object as the simulation unit, mainly describes the spatial distribution of various discrete geological features, and can well describe the geometric morphology of the river channel.

[0047] Target-based simulation takes the target body as the object. On the basis of the prior geometric shape of the geological body, the position of the center point of the geological body is randomly simulated, and the geometric shape and other attributes of the geological body are assigned at this point, and finally a three-dimensional model of the spatial distribution of the geological body is established. Although in the process of the microfacies simulation of the heart bar, the modeling personnel give the interval of the heart bar distribution direction according to the source direction of the study area, and the simulation results are consistent with the given direction interval and the overall direction of the source of the study area, but for one or several heart bars and their surrounding rivers, the heart bar distribution direction is quite different from the local river flow direction. Generally speaking, the main direction of a single heart bar is consistent with the corresponding local river flow direction, and the heart bar secondary direction is nearly parallel to the corresponding local river secondary direction. At the same time, the heart bar is usually developed in the middle of the river channel, which is the most typical development position of the heart bar; in a few cases, due to the influence of hydrodynamic conditions, sudden changes in riverbed morphology and sediment characteristics, it may develop on the edge of the river channel, resulting in a few cases where part of the heart bar is adjacent to the floodplain. However, the existing target-based simulation algorithms tend to simulate many complete mid-bar locations in the floodplain during the mid-bar simulation process, resulting in serious inconsistency between the simulation results and the actual geological development laws. In terms of the simulation effect of the existing algorithms, the existence of the above two problems poses a challenge to the establishment of an accurate geological modeling of braided river reservoirs that conforms to geological knowledge.

[0048] Therefore, in response to the above problems, the present invention provides a geological modeling method for simulating the center bar based on the distribution characteristics of the river channel. By obtaining the water flow azimuth of each grid in the river channel microfacies model of the braided river reservoir, an azimuth three-dimensional model of the river channel microfacies of the braided river reservoir is established, and then the azimuth model is used to constrain the center bar distribution direction. At the same time, the center bar microfacies development rules are set, and a center bar model is established to characterize the real distribution characteristics of the microfacies inside the river channel. The method described in the present invention provides a reference method for establishing a geological model of a braided river reservoir, and provides a reliable three-dimensional geological model for accurately simulating the location of the center bar and its distribution direction, as well as for later finding favorable areas.

[0049] Preferably, in another embodiment of the present application, the development characteristic data of the braided river reservoir in the study area includes the scope of the study area, well spacing, river channel geological data and mid-shoal geological data;

[0050] The river geological data include: river direction, river width, river thickness, river amplitude, river wavelength and river ratio;

[0051] The core-bar geological data include: core-bar shape, core-bar width, core-bar thickness, core-bar length and core-bar ratio.

[0052] The step of “S100, establishing a braided river reservoir channel microfacies model according to the development characteristic data” specifically includes the following steps:

[0053] Building a three-dimensional network grid, and setting the number of grids and the size of a single grid of the three-dimensional network grid;

[0054] Taking floodplain mud as the channel modeling background, based on the Fluvsim algorithm and according to the channel geological data in the development characteristic data, a braided river reservoir channel microfacies model including floodplain mud microfacies is established on the three-dimensional network framework.

[0055] Specifically, in this embodiment, the single grid step size and the number of grids in the three-dimensional network framework established by the present invention are comprehensively acquired based on the scope of the study area, well spacing and development characteristic data. The geological characteristics include: river direction, width, thickness, amplitude, wavelength and proportion. The study area can be set as a regular square, about 3km from east to west and about 5km from north to south. There are 20 wells in the study area, with an average well spacing of about 870m. In the three-dimensional grid model established in the present invention, the single grid step size is 20m*20m*0.5m, such as Figure 2 a; the number of grids is 150 (i direction) * 250 (j direction) * 40 (k direction) = 1,500,000, such as Figure 2 As shown in b, the model is enlarged 10 times in the vertical direction for demonstration effect.

[0056] The range of river channel direction is triangularly distributed, with the minimum value of river channel direction being -80°, the mode value of river channel direction being 0°, and the maximum value of river channel direction being 80°; the range of river channel amplitude is triangularly distributed, with the minimum value of river channel amplitude being 600m, the mode value of river channel amplitude being 650m, and the maximum value of river channel amplitude being 700m; the range of river channel wavelength value is triangularly distributed, with the minimum value of river channel wavelength being 2000m, the mode value of river channel wavelength being 3000m, and the maximum value of river channel wavelength being 4000m; the range of river channel width value is triangularly distributed, with the minimum river channel width value being 300m, the mode value being 400m, and the maximum river channel width value being 500m; the range of river channel thickness value is triangularly distributed, with the minimum river channel thickness value being 3m, the mode value being 6m, and the maximum river channel thickness value being 10m; the given value of the river channel proportion is 50%, and the river channel morphology is flat on the top and convex on the bottom, as shown in Table 1 below.

[0057] Table 1

[0058]

[0059] The simulation of the heart bar requires geological parameters, and the given heart bar ratio is 17%. The given heart bar profile is flat-bottomed and convex-topped, as shown in Table 2, which is consistent with geological knowledge. The heart bar direction does not give a specific value or value range. The heart bar direction is constrained by the azimuth three-dimensional model body. The heart bar width is triangularly distributed, with a minimum of 150m, a mode of 250m, and a maximum of 350m; the ratio of the heart bar width to length is triangularly distributed, with a minimum of 1.2, a mode of 2.2, and a maximum of 3. The heart bar thickness is triangularly distributed, with a minimum of 2m, a mode of 3m, and a maximum of 5m.

[0060] Table 2

[0061]

[0062] Therefore, according to the river geological data in the development characteristic data, the river microfacies model was established by using the target-based simulation method with the floodplain mud as the background, such as Figure 3 a. The results of the river channel microfacies model are statistically analyzed. The proportion of developed floodplain microfacies and developed river channel microfacies in this model is 50%. Figure 3 As shown in b.

[0063] Preferably, in another embodiment of the present application, the step of "S200, obtaining the water flow azimuth of each grid in the braided river reservoir channel microfacies model, and establishing an azimuth three-dimensional model of the braided river reservoir channel microfacies" specifically includes the following steps:

[0064] Obtaining the hydrodynamic axes of different single-layer river channels in the vertical direction and the boundary lines of the river channels on both sides in the braided river reservoir channel microfacies model;

[0065] Obtain the water flow azimuth corresponding to each grid passing through the water flow dynamic axis and the boundary lines of the river channels on both sides in each single-layer river channel;

[0066] According to the water flow azimuth corresponding to each grid passing through the water flow dynamic axis and the river channel boundary lines on both sides, the Kriging interpolation algorithm is used to calculate the water flow azimuth corresponding to each grid between the water flow dynamic axis and the river channel boundary lines on both sides;

[0067] According to the water flow azimuth corresponding to each grid in different single-layer river channels, a three-dimensional azimuth model of the braided river reservoir channel microfacies is established.

[0068] Specifically, in this embodiment, the water flow dynamic axis of each river channel is obtained according to the river channel microfacies model, such as Figure 4 The model has a single grid step of 0.5m in the vertical direction, with a total of 40 slices from top to bottom. Due to space limitations, this patent uses a step of 2.5m (every 4 grids apart) in the vertical direction to obtain the dynamic axis of the river flow, with a total of 9 slices.

[0069] The present invention further illustrates that the embodiment of the present invention is only an example, and the fineness of a single grid plane and vertical step length depends on the actual exploration and development needs of the oil and gas field, the size of the work area, and the computing power of the computer. According to the method proposed in the present invention, a slice that is finer vertically than the embodiment can be cut out, for example, the step length is 0.1m; a finer grid can more accurately simulate the flow direction of the river channel, but the model operation requires a higher computer itself; secondly, even if the grid scale is fine enough to depict a thin sand layer, the movable and recoverable reserves inside the sand body are extremely small, and from the perspective of economic benefits, it is not within the scope of actual exploitation of the oil and gas field. However, the fine grid can indeed play an important guiding role in the fine depiction of the geological characteristics of the reservoir space and the clarification of the distribution of sand bodies between underground wells. On the other hand, the vertical grid step length scale can also be larger, for example, the step length is 5m, which can meet the exploration scale requirements.

[0070] From top to bottom, they are the water flow dynamic axes when the thickness is 0m (K=1), such as Figure 4 a. When the thickness is 2.5m (K=5), the water flow dynamic axis is as follows: Figure 4 b. When the thickness is 5m (K=10), the water flow dynamic axis is as follows: Figure 4 c. When the thickness is 7.5m (K=15), the water flow dynamic axis is as follows: Figure 4 d. When the thickness is 10m (K=20), the water flow dynamic axis is as follows: Figure 4 e. When the thickness is 12.5m (K=25), the water flow dynamic axis is as follows: Figure 4 f. When the thickness is 15m (K=30), the water flow dynamic axis is as follows: Figure 4 g. When the thickness is 17.5m (K=35), the water flow dynamic axis is as follows: Figure 4 h. When the thickness is 20m (K=40), the water flow dynamic axis is as follows: Figure 4 As shown in i.

[0071] According to the river flow dynamic axis and river boundary line obtained from different single layers, the main direction of the river is determined. The single-period river flow dynamic axis and the river boundary lines on both sides are as follows: Figure 5 a. Figure 5Taking the area A shown in as an example, the azimuth data of grid A is obtained based on the water flow dynamic axis, and its value is measured to be 312.94°; at the same time, the azimuth data of the inner boundary line of the river channel in area A (the right boundary line in the figure) is measured, and the measurement result is 339.3°. Then, when the grid azimuths corresponding to the water flow dynamic axis and the inner boundary line of the river channel in area A have been obtained, the present invention uses the Kriging algorithm to interpolate the remaining 6 grid azimuths between the two lines, and the grid azimuth values ​​from the water flow dynamic axis to the inner boundary of the river channel are: 310.54°, 311.21°, 311.52°, 312.29°, 315.69°, 326.55°. Similarly, the azimuth of the hydrodynamic axis position in area B was measured to be 299.49°, and the grid azimuth data of the inner boundary of the river channel was 289.64°. The Kriging algorithm was used to interpolate the grid azimuths between the two lines, and the 6 grid azimuth data from the hydrodynamic axis to the inner boundary of the river channel were obtained as follows: 294.71°, 293.05°, 291.22°, 290.39°, 290.47°, 291.4°, as shown in Figure 2. Figure 5 As shown in b. Figure 6 As shown in a, it shows the change characteristics of water flow azimuth data of grids at different positions in a river channel with a depth of 0m. The step length of the single grid is 20m*20m in the plane and 0.5m in the vertical direction. From top to bottom, the azimuth angle data set of the internal grids of the river channel in 40 plane slices in the vertical direction is calculated, and the azimuth angle change characteristics of each grid on the river channel section are obtained from multiple plane slices, as shown in Figure 6 Finally, the azimuth angle three-dimensional model of the river microfacies was established by using the water flow azimuth data of 1.5 million grids in the study area, as shown in Figure 7 shown.

[0072] Preferably, in another embodiment of the present application, the step of "S300, taking the azimuth three-dimensional model body as a modeling constraint condition, and setting the development rules of the simulated center-bar microfacies, and establishing a center-bar model in the braided river reservoir channel microfacies model according to the development characteristic data" specifically includes the following steps:

[0073] The channel microfacies model of the braided river reservoir is used as the background for the barbed-bar modeling, and the development rules for simulating the barbed-bar microfacies are set;

[0074] The azimuth three-dimensional model body is used to constrain the distribution direction of the bar, and based on the Fluvsim algorithm and the bar geological data in the development characteristic data, a bar model including bar microfacies, river channel microfacies and floodplain mud microfacies is established.

[0075] Preferably, in another embodiment of the present application, the development rule of the center-bar microfacies is: searching for channel microfacies in the braided river reservoir channel microfacies model including the channel and floodplain mud, simulating the center-bar microfacies inside the searched channel microfacies, and the simulated center-bar microfacies does not contact the floodplain mud.

[0076] Specifically, in this embodiment, the channel microfacies model of the braided river reservoir is Figure 3 As shown in a, it is set as the background of the heart beach model, and the simulation rules are set during the heart beach simulation process. The rule execution statements are as follows:

[0077] M bar =If(Facies=Channel,Replace the Channel,Replace none other

[0078] facies)

[0079] Among them, M bar It refers to the median-shoal model, If is the conditional judgment statement to be executed, Facies refers to the microfacies type, and Channel refers to the river microfacies.

[0080] During the execution of this statement, if you want to simulate the center-bar microfacies, first find the channel microfacies in the channel microfacies model that contains the channel and floodplain mud as the background, and simulate the center-bar inside the channel, otherwise do not simulate the center-bar inside any other sedimentary microfacies. This statement is also applicable to the case where the braided river sedimentary system contains other sedimentary microfacies.

[0081] At the same time, before simulating the heart-bar model, it is necessary to first obtain the heart-bar geological data of the study area. This patent focuses on describing the heart-bar simulation position and heart-bar simulation direction. Therefore, when setting the heart-bar scale, the existing geological research results in the early stage are directly adopted. The embodiment of this patent only lists the heart-bar scale range, and no further explanation is given on how to obtain these data. However, it should be noted that this data is truly derived from modern sedimentation, outcrop sedimentation measurements, and verified by underground drilling core data and oil and gas well production curves.

[0082] According to the geological data of the heart bar, Table 2 above, the azimuth three-dimensional model body is used to constrain the direction characteristics of the heart bar, and the conditional statement is used to constrain the location characteristics of the heart bar development. The target-based simulation method is used to establish the heart bar model (including three microfacies types: floodplain, river channel, and heart bar), such as Figure 8 From the simulation results, it is found that the method proposed in this patent constrains the main direction of the center bar according to the direction of river flow, especially in areas with high river curvature, which is consistent with the actual geological characteristics, such as Figure 8As shown. The bar follows the development rules of the present invention and develops only inside the river channel, which avoids the problem that when only the simulation of the bar direction interval is provided, the bar is cut off from the center at the boundary between the river channel and the floodplain microfacies, resulting in a large area of ​​the bar contacting the floodplain microfacies, which is seriously inconsistent with the actual geological understanding.

[0083] The embodiment of the present invention further provides a geological modeling system for simulating a mid-shoal based on river channel distribution characteristics, comprising:

[0084] The channel facies modeling module is used to transmit the development characteristic data of the braided river reservoir in the study area to the PETREL modeling software, and establish the channel microfacies model of the braided river reservoir based on the development characteristic data;

[0085] A water flow azimuth module is connected to the channel phase modeling module for obtaining the water flow azimuth of each grid in the braided river reservoir channel microfacies model and establishing a three-dimensional azimuth model of the braided river reservoir channel microfacies; and

[0086] The center-bar modeling module is communicatively connected with the water flow azimuth module, and is used to use the azimuth three-dimensional model body as a modeling constraint condition, set the development rules of the simulated center-bar microfacies, and establish a center-bar model in the braided river reservoir channel microfacies model according to the development characteristic data.

[0087] Because the modelers gave the interval of the distribution direction of the bar according to the source direction of the study area in the process of shoal microfacies simulation, the simulation results are consistent with the given direction interval and the overall direction of the source of the study area. However, for one or several shoals and their surrounding rivers, the distribution direction of the bar is quite different from the flow direction of the local river. Generally speaking, the main direction of a single bar is consistent with the flow direction of the corresponding local river, such as Figure 7 As shown, the sub-direction of the heart bar is nearly parallel to the sub-direction of the corresponding local river channel. At the same time, the heart bar is usually developed in the middle of the river channel, which is the most typical development position of the heart bar; in a few cases, due to the influence of hydrodynamic conditions, sudden changes in riverbed morphology and sediment characteristics, it may develop on the edge of the river channel, resulting in that in very rare cases, part of the heart bar is adjacent to the floodplain. However, in the process of simulating the heart bar, the existing target-based simulation algorithm is prone to simulate many complete heart bar positions in the floodplain, resulting in the simulation results being seriously inconsistent with the actual geological development laws. In terms of the simulation effect of the existing algorithm, the existence of the above two problems poses a challenge to the establishment of an accurate geological modeling of braided river reservoirs that conforms to geological knowledge. To this end, the present invention proposes a geological modeling method for simulating the heart bar based on the distribution characteristics of the river channel to solve the above two problems.

[0088] Specifically, this embodiment corresponds one-to-one to the above method embodiment, and the functions of each module have been described in detail in the corresponding method embodiment, so they will not be repeated here.

[0089] Based on the same inventive concept, an embodiment of the present application further provides a computer-readable storage medium on which a computer program is stored. When the computer program is executed by a processor, all or part of the method steps of the above method are implemented.

[0090] The present invention implements all or part of the process in the above method, and can also be completed by instructing related hardware through a computer program. The computer program can be stored in a computer-readable storage medium. When the computer program is executed by a processor, the steps of each method embodiment can be implemented. Among them, the computer program includes computer program code, and the computer program code can be in source code form, object code form, executable file or some intermediate form. The computer-readable medium may include: any entity or device capable of carrying computer program code, recording medium, USB flash drive, mobile hard disk, disk, optical disk, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electric carrier signal, telecommunication signal and software distribution medium. It should be noted that the content contained in the computer-readable medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, the computer-readable medium does not include electric carrier signals and telecommunication signals.

[0091] Based on the same inventive concept, an embodiment of the present application also provides an electronic device, including a memory and a processor, wherein the memory stores a computer program running on the processor, and when the processor executes the computer program, all or part of the method steps in the above method are implemented.

[0092] The processor may be a central processing unit (CPU), other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc. The processor is the control center of a computer device, and uses various interfaces and lines to connect various parts of the entire computer device.

[0093] The memory can be used to store computer programs and / or modules, and the processor realizes various functions of the computer device by running or executing the computer programs and / or modules stored in the memory, and calling the data stored in the memory. The memory can mainly include a program storage area and a data storage area, wherein the program storage area can store an operating system, an application required for at least one function (such as a sound playback function, an image playback function, etc.); the data storage area can store data created according to the use of the mobile phone (such as audio data, video data, etc.). In addition, the memory can include a high-speed random access memory, and can also include a non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a smart memory card (SmartMed i aCard, SMC), a secure digital (Secure Digital, SD) card, a flash card (Flash Card), at least one disk storage device, a flash memory device, or other volatile solid-state storage devices.

[0094] It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as methods, systems, servers or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage and optical storage, etc.) containing computer-usable program codes.

[0095] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), servers, and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as a combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0096] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0097] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0098] Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include these modifications and variations.

Claims

1. A geological modeling method for simulating the center shoal based on the distribution characteristics of the river channel, characterized in that: The following steps are involved: The development characteristic data of the braided river reservoir in the study area were transferred to the PETREL modeling software, and a channel microfacies model of the braided river reservoir was established based on the development characteristic data; Obtaining the water flow azimuth of each grid in the braided river reservoir channel microfacies model, and establishing a three-dimensional azimuth model of the braided river reservoir channel microfacies; The azimuth three-dimensional model body is used as a modeling constraint condition, and the development rules for simulating the center-bar microfacies are set, and a center-bar model is established in the braided river reservoir channel microfacies model according to the development characteristic data.

2. The geological modeling method for simulating the central shoal based on the river channel distribution characteristics according to claim 1, characterized in that: The development characteristic data of braided river reservoirs in the study area include the scope of the study area, well spacing, river channel geological data and mid-shoal geological data; The river geological data include: river direction, river width, river thickness, river amplitude, river wavelength and river ratio; The core-bar geological data include: core-bar shape, core-bar width, core-bar thickness, core-bar length and core-bar ratio.

3. The geological modeling method for simulating the central shoal based on the river channel distribution characteristics according to claim 1, characterized in that: The step of "establishing a braided river reservoir channel microfacies model based on the development characteristic data" specifically includes the following steps: Building a three-dimensional network grid, and setting the number of grids and the size of a single grid of the three-dimensional network grid; Taking floodplain mud as the channel modeling background, based on the Fluvsim algorithm and according to the channel geological data in the development characteristic data, a braided river reservoir channel microfacies model including floodplain mud microfacies is established on the three-dimensional network framework.

4. The geological modeling method for simulating the central shoal based on the river channel distribution characteristics according to claim 1, characterized in that: The step of "obtaining the water flow azimuth of each grid in the braided river reservoir channel microfacies model and establishing the azimuth three-dimensional model of the braided river reservoir channel microfacies" specifically includes the following steps: Obtaining the hydrodynamic axes of different single-layer river channels in the vertical direction and the boundary lines of the river channels on both sides in the braided river reservoir channel microfacies model; Obtain the water flow azimuth corresponding to each grid passing through the water flow dynamic axis and the boundary lines of the river channels on both sides in each single-layer river channel; According to the water flow azimuth corresponding to each grid passing through the water flow dynamic axis and the river channel boundary lines on both sides, the Kriging interpolation algorithm is used to calculate the water flow azimuth corresponding to each grid between the water flow dynamic axis and the river channel boundary lines on both sides; According to the water flow azimuth corresponding to each grid in different single-layer river channels, a three-dimensional azimuth model of the braided river reservoir channel microfacies is established.

5. The geological modeling method for simulating the central shoal based on the river channel distribution characteristics according to claim 1, characterized in that: The step of "taking the azimuth three-dimensional model body as a modeling constraint condition, setting the development rules of the simulated center-bar microfacies, and establishing the center-bar model in the braided river reservoir channel microfacies model according to the development characteristic data" specifically includes the following steps: The channel microfacies model of the braided river reservoir is used as the background for the barbed-bar modeling, and the development rules for simulating the barbed-bar microfacies are set; The azimuth three-dimensional model body is used to constrain the distribution direction of the bar, and based on the Fluvsim algorithm and the bar geological data in the development characteristic data, a bar model including bar microfacies, river channel microfacies and floodplain mud microfacies is established.

6. The geological modeling method for simulating the central shoal based on the river channel distribution characteristics according to claim 1, characterized in that: The development rule of the center-bar microfacies is: searching for channel microfacies in the braided river reservoir channel microfacies model including the channel and floodplain mud, simulating the center-bar microfacies inside the searched channel microfacies, and the simulated center-bar microfacies does not contact the floodplain mud.

7. A geological modeling system for simulating the central shoal based on the distribution characteristics of the river channel, characterized in that: include: The channel facies modeling module is used to transmit the development characteristic data of the braided river reservoir in the study area to the PETREL modeling software, and establish the channel microfacies model of the braided river reservoir based on the development characteristic data; A water flow azimuth module is connected to the channel phase modeling module for obtaining the water flow azimuth of each grid in the braided river reservoir channel microfacies model and establishing a three-dimensional azimuth model of the braided river reservoir channel microfacies; and The center-bar modeling module is communicatively connected with the water flow azimuth module, and is used to use the azimuth three-dimensional model body as a modeling constraint condition, set the development rules of the simulated center-bar microfacies, and establish a center-bar model in the braided river reservoir channel microfacies model according to the development characteristic data.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the geological modeling method for simulating the center bar based on the distribution characteristics of the river channel as described in any one of claims 1 to 6 is implemented.

9. An electronic device comprising a storage medium, a processor, and a computer program stored in the storage medium and executable on the processor, characterized in that: When the processor runs the computer program, it implements the geological modeling method for simulating the center bar based on the river channel distribution characteristics as described in any one of claims 1 to 6.

Images