Geological modeling method and system for simulating the center bank based on river channel distribution characteristics

Abstract

The present invention discloses a geological modeling method and system for simulating a center bar based on river channel distribution characteristics. The method comprises the following steps: transmitting development characteristic data of a braided river reservoir in a study area to PETREL modeling software, and establishing a channel microfacies model of the braided river reservoir according to the development characteristic data; obtaining the water flow azimuth of each grid in the channel microfacies model of the braided river reservoir, and establishing an azimuth three-dimensional model body of the channel microfacies of the braided river reservoir; using the azimuth three-dimensional model body as a modeling constraint condition, setting development rules 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. Therefore, the present invention can accurately simulate the center bar position and center bar distribution direction.

Description

Technical Field

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

[0002] Braided rivers are an important type of terrestrial sedimentary system. Due to the strong heterogeneity of their reservoirs and the complex spatial stacking and distribution patterns of sand bodies, establishing reliable three-dimensional geological models of braided river reservoirs remains challenging. Architectural element analysis is a key research technique for reservoir characterization and a current hotspot. Numerous scholars have conducted in-depth studies on braided river reservoir architecture and conducted extensive research on braided river reservoir modeling. The channel and mid-channel-bar architectural units are the most important architectural units in braided river reservoirs. The mid-channel-bar architectural unit exhibits greater homogeneity and superior reservoir quality compared to the channel architectural unit, making quantitative research on its distribution characteristics of great significance.

[0003] Common facies modeling methods include target-based simulation, multi-point geostatistics, and sequential indicator simulation. A detailed comparison of the basic theories and algorithmic processes of several modeling methods reveals that each has its own advantages and disadvantages, and should be appropriately selected based on the specific circumstances and actual needs, taking into account the applicable conditions, advantages, and disadvantages of each method. Target-based simulation, among others, uses the target body as the object. Based on the prior knowledge of the geologic body's geometry, it randomly simulates the position of the geologic body's center point, assigns the geologic body's geometry and other attributes at that point, and ultimately establishes a three-dimensional model of the geologic body's spatial distribution. However, existing target-based simulation methods tend to simulate the center bar within the floodplain during the process of simulating the center bar, resulting in simulation results that are seriously inconsistent with the actual geological development patterns. Summary of the Invention

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

[0005] In a first aspect, a geological modeling method for simulating a mid-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 braided river reservoir channel microfacies model 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 development rules for simulating the center-bar microfacies are set. 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 center-bar geological data include center-bar shape, center-bar width, center-bar thickness, center-bar length and center-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 based on 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 each grid in 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 a three-dimensional azimuth model of the braided river reservoir channel microfacies” specifically includes the following steps:

[0016] Obtaining the flow dynamic 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 river boundary lines 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 boundary lines on both sides, the water flow azimuth corresponding to each grid between the water flow dynamic axis and the river boundary lines on both sides is calculated using the Kriging interpolation algorithm;

[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 “using the azimuth three-dimensional model as a modeling constraint, setting development rules 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” specifically includes the following steps:

[0021] The braided river reservoir channel microfacies model is used as the core-bar modeling background, and the development rules of the simulated core-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 bar-bar microfacies is: searching for channel microfacies in the braided river reservoir channel microfacies model including the channel and floodplain mud, simulating the bar-bar microfacies inside the searched channel microfacies, and the simulated bar-bar microfacies does not contact the floodplain mud.

[0024] Secondly, a geological modeling system for simulating the center bank based on the distribution characteristics of the river channel 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, which is in communication with the channel facies modeling module and is used to obtain the water flow azimuth of each grid in the braided river reservoir channel microfacies model and establish a three-dimensional azimuth model of the braided river reservoir channel microfacies; and

[0027] The center-bar modeling module is in communication 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 center bar based on the river channel distribution characteristics 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 executes the computer program, the geological modeling method for simulating the heart bar based on the river channel distribution characteristics as described above is implemented.

[0030] Compared with existing technologies, the present invention offers the following advantages: By obtaining the flow azimuth of each grid cell in a braided river reservoir channel microfacies model, a three-dimensional azimuth model of the braided river reservoir channel microfacies is constructed. The azimuth model then constrains the distribution direction of the barbed channel. Simultaneously, the barbed channel microfacies development rules are established to construct a barbed channel model that characterizes the true distribution characteristics of the microfacies within the channel. This method provides a valuable reference for establishing geological models of braided river reservoirs, providing a reliable three-dimensional geological model for accurately simulating the location and distribution of the barbed channel, and for later identifying favorable areas. BRIEF DESCRIPTION OF THE DRAWINGS

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

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

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

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

[0035] Figure 5 This 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 of the river direction established based on the water flow dynamic axis of the present invention;

[0038] Figure 8 The present 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 the present invention is not intended to be limited to those embodiments. On the contrary, it is intended to cover variations, modifications, and equivalents within 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 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 with reference to the accompanying drawings and specific embodiments.

[0041] Note: The following example 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, sequence, etc. Those skilled in the art can apply the concepts of the present invention to construct more embodiments not described 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 mid-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 braided river reservoir channel microfacies model based on 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 as a modeling constraint, setting development rules 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 modeling complex river channels, also known as the "object-based simulation method." This modeling method uses target objects as simulation units and primarily describes the spatial distribution of various discrete geological features, effectively depicting the geometric morphology of river channels.

[0047] Target-based simulation uses a target object as the model. Based on the a priori geometry of the geological body, the center point of the geological body is randomly simulated. Geometric properties and other attributes are assigned at this point, ultimately creating a three-dimensional model of the geological body's spatial distribution. Although during the bar microfacies simulation process, modelers assign a range of bar distribution directions based on the provenance of the study area, and the simulation results agree well with the given range and the overall provenance of the study area, the bar distribution direction differs significantly from the local channel flow direction for individual or several bar and surrounding channels. Generally speaking, the primary direction of a single bar aligns with the flow direction of its corresponding local channel, and the secondary direction of the bar runs nearly parallel to that of the corresponding local channel. Furthermore, bar formations typically develop in the middle of a channel, their most typical location. In rare cases, due to hydrodynamic conditions, sudden changes in riverbed morphology, and sediment properties, they may develop at the channel edge, resulting in rare cases where portions of the bar are adjacent to the floodplain. However, existing target-based simulation algorithms often simulate many intact bar locations within the floodplain, resulting in simulation results that are seriously inconsistent with actual geological development patterns. These two issues pose challenges to the development of accurate, geologically accurate braided river reservoir models.

[0048] To address these issues, the present invention provides a geological modeling method for simulating channel-shoals based on channel distribution characteristics. By obtaining the water flow azimuth for each grid cell in a braided river reservoir channel microfacies model, a three-dimensional azimuth model of the braided river reservoir channel microfacies is constructed. The azimuth model then constrains the distribution direction of the channel-shoal. Simultaneously, the development rules for the channel-shoal microfacies are defined, and a channel-shoal model is established to characterize the true distribution characteristics of the microfacies within the channel. This method provides a valuable reference for establishing geological models of braided river reservoirs, providing a reliable three-dimensional geological model for accurately simulating the location and distribution of channel-shoals and for later identifying 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 center-bar geological data include center-bar shape, center-bar width, center-bar thickness, center-bar length and center-bar ratio.

[0052] The step of “S100, establishing a braided river reservoir channel microfacies model based on 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 each grid in 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 obtained 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 As shown in a; the number of grids is 150 (i direction) * 250 (j direction) * 40 (k direction) = 1,500,000, as shown in 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 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 is triangularly distributed, with the minimum value of river channel width being 300m, the mode value being 400m, and the maximum value of river channel width being 500m; the range of river channel thickness is triangularly distributed, with the minimum value of river channel thickness being 3m, the mode value being 6m, and the maximum value of river channel thickness being 10m; the given value of river channel proportion is 50%, and the river channel morphology is flat on top and convex on the bottom, as shown in Table 1 below.

[0057] Table 1

[0058]

[0059] The simulation of the core bar requires geological parameters, and a core bar ratio of 17% is assumed. The core bar profile is given as a flat-bottomed, convex-topped shape, as shown in Table 2, which is consistent with geological understanding. The core bar orientation is not given a specific value or range; instead, it is constrained by the azimuth 3D model. The core bar width exhibits a triangular distribution, with a minimum of 150 m, a mode of 250 m, and a maximum of 350 m. The core bar width-to-length ratio also exhibits a triangular distribution, with a minimum of 1.2, a mode of 2.2, and a maximum of 3. The core bar thickness also exhibits a triangular distribution, with a minimum of 2 m, a mode of 3 m, and a maximum of 5 m.

[0060] Table 2

[0061]

[0062] Therefore, based on the channel geological data in the development characteristic data, the channel microfacies model was established using the target-based simulation method with the floodplain mud as the background, such as Figure 3 As shown in 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%, as shown in 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 flow dynamic 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 river boundary lines 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 boundary lines on both sides, the water flow azimuth corresponding to each grid between the water flow dynamic axis and the river boundary lines on both sides is calculated using the Kriging interpolation algorithm;

[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 based on 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 in the vertical direction (every 4 grids) to obtain the dynamic axis of the river flow, with a total of 9 slices.

[0069] The present invention further explains that the embodiment of the present invention is merely an example. The degree of refinement of individual grid planes and vertical step sizes depends on the actual exploration and development requirements 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, finer vertical slices can be cut than those in the embodiment, for example, with a step size of 0.1m. Finer grids can more accurately simulate the flow direction within the river channel, but the model calculations place higher demands on the computer itself. Secondly, even if the grid scale is fine enough to depict thin sand layers, the movable and recoverable reserves within the sand bodies are extremely small, and from an economic perspective, they are not within the scope of actual oil and gas field development. However, fine grids do play an important guiding role in accurately depicting the geological characteristics of the reservoir space and clarifying the distribution of sand bodies between underground wells. On the other hand, the vertical grid step size can also be larger, for example, a step size of 5m, to 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 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 As shown in a. Figure 5Taking the area A shown in the figure 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 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°, and 291.4°. Figure 5 As shown in b. Figure 6 As shown in Figure a, the water flow azimuth data variation characteristics of grids at different locations in a river channel with a depth of 0m are shown. The step size of the single grid is 20m*20m in the plane and 0.5m in the vertical direction. From top to bottom, the azimuth angle dataset of the internal grid of the river channel is calculated in 40 vertical plane slices. The azimuth angle variation characteristics of each grid on the river channel section are obtained from multiple plane slices, as shown in Figure 4. Figure 6 Finally, the azimuth angle three-dimensional model of the river microfacies was established by combining 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, using the azimuth three-dimensional model body as a modeling constraint, setting the development rules 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" specifically includes the following steps:

[0073] The braided river reservoir channel microfacies model is used as the core-bar modeling background, and the development rules of the simulated core-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 containing river channels and floodplain mud, simulating 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 braided river reservoir channel microfacies model, such as 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 channel microfacies.

[0080] When executing this statement, if you want to simulate a channel microfacies, first search for the channel microfacies within the channel microfacies model that includes the channel and floodplain mud as the background. Then simulate the channel within the channel; otherwise, do not simulate a channel within any other sedimentary microfacies. This statement also applies to braided river sedimentary systems that contain 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 simulated heart-bar location and heart-bar simulation direction. Therefore, when setting the heart-bar scale, the existing geological research results are directly adopted. This patent embodiment only lists the heart-bar scale range and does not further explain how to obtain this data. However, it should be noted that this data is authentically derived from modern sedimentary and outcrop sedimentary 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 development location characteristics of the heart bar. The target-based simulation method is used to establish the heart bar model (including three microfacies types: floodplain, river channel, and heart bar), as shown in the following example: Figure 8 From the simulation results, it is found that the method proposed in this patent constrains the main direction of the center bank 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 8The bar follows the development rules of the present invention and develops only within the channel. This avoids the problem of truncating the bar from the center at the boundary between the channel and floodplain microfacies when only providing the simulation of the bar direction interval, resulting in a large area of ​​the bar contacting the floodplain microfacies, which is seriously inconsistent with 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, which is in communication with the channel facies modeling module and is used to obtain the water flow azimuth of each grid in the braided river reservoir channel microfacies model and establish a three-dimensional azimuth model of the braided river reservoir channel microfacies; and

[0086] The center-bar modeling module is in communication 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 range of the bar distribution direction according to the source direction of the study area in the process of the bar microfacies simulation, the simulation results were consistent with the given direction range and the overall source direction of the study area. However, for one or several bar and its surrounding river channels, the bar distribution direction was quite different from the flow direction of the local river channel. Generally speaking, the main direction of a single bar is consistent with the flow direction of the corresponding local river channel, such as Figure 7 As shown in the figure, the sub-direction of the bar is nearly parallel to the sub-direction of the corresponding local river channel. At the same time, the bar usually develops in the middle of the river channel, which is the most typical development location of the bar. In a few cases, due to the influence of hydrodynamic conditions, sudden changes in riverbed morphology and sediment characteristics, it may develop at the edge of the river channel, resulting in part of the bar adjacent to the floodplain in very rare cases. However, in the process of simulating the bar, the existing target-based simulation algorithm is prone to simulating many complete bar locations 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 model of braided river reservoirs that conforms to geological knowledge. To this end, the present invention proposes a geological modeling method for simulating 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, it can implement the steps of each of the above method embodiments. 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 that can carry computer program code, recording medium, USB flash drive, mobile hard disk, magnetic disk, optical disk, computer memory, read-only memory (ROM), random access memory (RAM), 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, computer-readable media do 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), 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. The processor is the control center of a computer device and connects various parts of the entire computer device using various interfaces and lines.

[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, at least one application required for a 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] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, servers, or computer program products. Thus, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware. Furthermore, 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 magnetic disk storage and optical storage) containing computer-usable program code.

[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 flowcharts and / or block diagrams, as well as combinations of processes and / or blocks in the flowcharts and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, 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 flowcharts and / or block diagrams. 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 that can direct a computer or other programmable data processing device to work 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 The 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 operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device 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 may make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if such changes and modifications fall within the scope of the claims and their equivalents, the present invention is intended to include such changes and modifications.

Claims

1. A geological modeling method for simulating the center bank 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 braided river reservoir channel microfacies model 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; Taking the azimuth three-dimensional model body as a modeling constraint condition and setting the development rules of the simulated center-bar microfacies, a center-bar model is established in the braided river reservoir channel microfacies model according to the development characteristic data; 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 flow dynamic 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 river boundary lines 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 boundary lines on both sides, the water flow azimuth corresponding to each grid between the water flow dynamic axis and the river boundary lines on both sides is calculated using the Kriging interpolation algorithm; 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.

2. The geological modeling method for simulating a mid-shoal based on river channel distribution characteristics according to claim 1, characterized in that: The development characteristic data of the braided river reservoir 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 center-bar geological data include center-bar shape, center-bar width, center-bar thickness, center-bar length and center-bar ratio.

3. The geological modeling method for simulating a central shoal based on 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 each grid in 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 a mid-shoal based on river channel distribution characteristics according to claim 1, characterized in that: The step of "using the azimuth three-dimensional model as a modeling constraint, setting the development rules for simulating the center-bank microfacies, and establishing a center-bank model in the braided river reservoir channel microfacies model according to the development characteristic data" specifically includes the following steps: The braided river reservoir channel microfacies model is used as the core-bar modeling background, and the development rules of the simulated core-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.

5. The geological modeling method for simulating the center bar 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 containing 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.

6. A geological modeling system for simulating the center bank based on the distribution characteristics of the river channel, characterized by: 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, which is in communication with the channel facies modeling module and is used to obtain the water flow azimuth of each grid in the braided river reservoir channel microfacies model and establish a three-dimensional azimuth model of the braided river reservoir channel microfacies; and a center-bar modeling module, communicatively connected to the water flow azimuth module, for using the azimuth three-dimensional model as a modeling constraint, setting development rules for simulating center-bar microfacies, and establishing a center-bar model in the braided river reservoir channel microfacies model according to the development characteristic data; 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 flow dynamic 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 river boundary lines 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 boundary lines on both sides, the water flow azimuth corresponding to each grid between the water flow dynamic axis and the river boundary lines on both sides is calculated using the Kriging interpolation algorithm; 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.

7. 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 a mid-shoal based on river channel distribution characteristics according to any one of claims 1 to 5 is implemented.

8. An electronic device 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 geological modeling method for simulating the center bar based on the distribution characteristics of the river channel according to any one of claims 1 to 5 is implemented.

Images