A sandstone thickness mapping method and system, an electronic device and a storage medium

Abstract

The present disclosure relates to a sandstone thickness mapping method, system, electronic device and storage medium. The mapping method comprises: obtaining logging data of each well point in a preset area, and determining sandstone and mudstone lithological characteristic data according to the logging data; determining sandstone curve shape and sandstone thickness of each well point according to the sandstone and mudstone lithological characteristic data, and mapping a sedimentary facies pattern of the preset area through the sandstone curve shape and the sandstone thickness; converting the sedimentary facies pattern to a constraint plane for determining the trend change of a sedimentary facies belt in a three-dimensional geological framework model plane range; generating plane sandstone thickness data through the constraint plane based on a geostatistical stochastic simulation algorithm, and mapping a sandstone thickness map of the preset area according to the plane sandstone thickness data. The sandstone thickness is generated by combining the geological statistical principle and the sedimentary facies pattern, and the sandstone thickness map is mapped according to the sandstone thickness, so that high-precision prediction and visual expression of the sandstone thickness are realized.

Description

Technical Field

[0001] This disclosure belongs to the field of geological exploration technology, and in particular relates to a method, system, electronic device and storage medium for drawing sandstone thickness maps. Background Technology

[0002] In the field of oil and gas exploration and development, geological structure and sedimentary characteristics are important research tasks. The spatial distribution and thickness of sandstone bodies directly affect the accumulation and distribution of oil and gas. Sandstone thickness mapping is one of the key technologies for assessing oil and gas exploration potential and distribution, and accurate sandstone thickness mapping has significant guiding significance for oil and gas exploration and development. The technical background of sandstone thickness mapping involves the interdisciplinary application of multiple disciplines such as geology, geophysics, and computer science. Geologically, the formation and evolution of sandstone bodies require analysis using stratigraphy, sedimentology, and paleogeography. By reconstructing the sedimentary environment, the original thickness and distribution range of sandstone bodies can be inferred. Furthermore, stratigraphic contact relationships, lithological variations, and paleogeographic features are also important bases for sandstone thickness mapping.

[0003] Currently, the most common method for compiling sandstone thickness maps is the grid interpolation method (including inverse distance weighted interpolation, Kriging interpolation, spline interpolation, etc.). However, in practical applications, the drawing of sandstone thickness maps also needs to consider the superimposed effects of various geological factors, such as the distribution of sedimentary facies zones, erosion, and tectonic movements. These factors may lead to changes in sandstone thickness, or even the disappearance of sand bodies. Summary of the Invention

[0004] To address the aforementioned issues, this disclosure provides a method, system, electronic device, and storage medium for drawing sandstone thickness maps based on sedimentary facies constraints. By combining geostatistical principles and sedimentary facies models to generate sandstone thickness, and drawing sandstone thickness maps based on the sandstone thickness, high-precision prediction and visualization of sandstone thickness are achieved.

[0005] To address the aforementioned technical problems, the first aspect of this invention proposes a method for drawing sandstone thickness maps based on sedimentary facies constraints, the method comprising:

[0006] Obtain logging data from each well point in a preset area, and determine the lithological characteristics of sandstone and mudstone based on the logging data;

[0007] Based on the lithological characteristics data of the sandstone and mudstone, the sandstone curve morphology and sandstone thickness of each well point are determined, and the sedimentary facies model diagram of the preset area is drawn based on the sandstone curve morphology and sandstone thickness.

[0008] Within the three-dimensional geological framework model plane, the sedimentary facies pattern diagram is transformed into a constrained plane used to determine the trend changes of sedimentary facies zones;

[0009] Based on a geostatistical stochastic simulation algorithm, planar sandstone thickness data is generated by constraining the constrained plane, and a sandstone thickness map of the preset area is drawn based on the planar sandstone thickness data.

[0010] According to a preferred embodiment of the present invention, the step of acquiring logging data from each well point in a preset area and determining sandstone and mudstone lithological characteristic data based on the logging data includes:

[0011] Obtain the natural gamma logging curve data of each sampling point in each well point of the preset area;

[0012] The mud content of the corresponding sampling point is calculated based on the maximum and minimum values ​​of the natural gamma logging curve data, and is used as the lithological characteristic data of the sandstone and mudstone.

[0013] According to a preferred embodiment of the present invention, the method further includes:

[0014] Obtain the sandstone thickness of the target layer at each well point;

[0015] The maximum thickness of the sandstone is used as the bottom surface of the frame model, and the minimum thickness of the sandstone is used as the top surface of the frame model to generate a three-dimensional geological frame model.

[0016] According to a preferred embodiment of the present invention, the step of converting the sedimentary facies pattern diagram to a constrained plane for determining the trend changes of sedimentary facies zones within the three-dimensional geological framework model plane includes:

[0017] Obtain the facies zones in the sedimentary facies pattern diagram and number the facies zones;

[0018] Within the plane of the three-dimensional geological framework model, a constraint plane is generated based on the sedimentary facies pattern diagram and the numbering.

[0019] According to a preferred embodiment of the present invention, the step of generating planar sandstone thickness data by constraining the plane using a geostatistical stochastic simulation algorithm includes:

[0020] Based on the aforementioned geostatistical stochastic simulation algorithm, the constrained plane is used as a constraint condition to generate the sandstone thickness distribution, sandstone thickness randomness, and sandstone thickness heterogeneity, which are then used as the planar sandstone thickness data.

[0021] According to a preferred embodiment of the present invention, the step of drawing a sedimentary facies model diagram of the preset area using the sandstone curve morphology and sandstone thickness includes:

[0022] Based on the sandstone curve morphology and sandstone thickness, the dominant sedimentary facies of the target layer at each well point are determined.

[0023] Based on the dominant sedimentary facies of each well point target layer, the corresponding well point target layers are divided into sedimentary facies, and a sedimentary facies pattern diagram of the preset area is drawn by combining the sedimentary facies of each well point.

[0024] According to a preferred embodiment of the present invention, the step of drawing a sandstone thickness map of the preset area based on the planar sandstone thickness data includes:

[0025] Based on the plane sandstone thickness data, draw a contour map of sandstone thickness in the preset area.

[0026] To address the aforementioned technical problems, a second aspect of the present invention proposes a sandstone thickness map drawing system based on sedimentary facies constraints, the drawing system comprising:

[0027] The data analysis module is used to acquire logging data from various well points in a preset area; and to determine the lithological characteristics of sandstone and mudstone based on the logging data.

[0028] The graphics drawing module is used to determine the sandstone curve shape and sandstone thickness of each well point based on the sandstone and mudstone lithological characteristic data, and to draw the sedimentary facies pattern diagram of the preset area based on the sandstone curve shape and sandstone thickness.

[0029] The constraint plane conversion module is used to convert the sedimentary facies pattern diagram to a constraint plane for determining the trend changes of sedimentary facies zones within the plane range of the three-dimensional geological framework model.

[0030] The sandstone thickness generation module is used to generate planar sandstone thickness data based on a geostatistical stochastic simulation algorithm and constrained by the constrained plane.

[0031] The graphics drawing module is also used to draw a sandstone thickness map of the preset area based on the plane sandstone thickness data.

[0032] To address the aforementioned technical problems, a third aspect of the present invention provides an electronic device, comprising:

[0033] Processor; and

[0034] A memory storing computer-executable instructions, which, when executed, cause the processor to perform the method described in any of the above embodiments.

[0035] To address the aforementioned technical problems, a fourth aspect of the present invention provides a computer storage medium, wherein the computer storage medium stores one or more programs, which, when executed by a processor, implement the method described in any of the above embodiments.

[0036] Compared with existing technologies, this disclosure has the following advantages: This disclosure determines the lithological characteristics of each well point using logging data, and statistically obtains the sandstone curve morphology and sandstone thickness at each well point based on the lithological characteristics data. A sedimentary facies model diagram of the area is then drawn. This sedimentary facies model diagram is converted into a constraint plane within the three-dimensional geological framework model plane, allowing the user to determine the trend changes of sedimentary facies zones. Using this constraint plane, a geostatistical stochastic simulation algorithm is used to predict the sandstone thickness distribution and determine the planar sandstone thickness data. Finally, a sandstone thickness map is drawn using this planar sandstone thickness data. This solution, by combining geostatistical principles and sedimentary facies models, achieves high-precision prediction and visualization of sandstone thickness. Specific beneficial effects include:

[0037] (1) Improve prediction accuracy: By introducing sedimentary facies models, the spatial distribution characteristics of sandstone thickness can be captured more accurately. The constraints of sedimentary facies models make the interpolation process consider not only the distance between data points, but also the influence of the sedimentary environment on the distribution of sandstone thickness, thereby improving the accuracy of prediction.

[0038] (2) Enhanced geological interpretation capabilities: The method of this invention allows geologists to incorporate their expertise and understanding of sedimentary processes into the creation of sandstone thickness maps. This combination of geostatistics and sedimentology provides deeper insights into geological interpretation and contributes to a better understanding of the reservoir characteristics of oil and gas reservoirs.

[0039] (3) Optimize exploration strategies: Accurate sandstone thickness maps provide key geological information for oil and gas exploration, which helps exploration teams to develop more effective exploration strategies, such as determining the best drilling locations and assessing potential oil and gas resources.

[0040] (4) Reduced exploration risks: By reducing prediction uncertainty, this invention helps to reduce risks in the oil and gas exploration and development process. Accurate sandstone thickness prediction can help avoid high-cost drilling activities in areas with low resource potential.

[0041] (5) Improved resource assessment efficiency: The automated device of this invention can quickly generate sandstone thickness maps, improving the efficiency of resource assessment. This rapid response capability is particularly important in the rapidly changing oil and gas market, helping companies to adjust exploration and development plans in a timely manner.

[0042] Other features and advantages of this disclosure will be set forth in the description which follows, and will be apparent in part from the description, or may be learned by practicing the disclosure. The objects and other advantages of this disclosure may be realized and obtained by means of the structures pointed out in the description, claims and drawings. Attached Figure Description

[0043] To more clearly illustrate the technical solutions in the embodiments of this disclosure or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of this disclosure. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0044] Figure 1 A schematic flowchart of a method for drawing a sandstone thickness map according to an embodiment of the present disclosure is shown;

[0045] Figure 2 This illustrates single-well logging data for a preset area according to an embodiment of the present disclosure;

[0046] Figure 3 A sedimentary facies pattern diagram obtained by integrating data from multiple wells according to an embodiment of this disclosure is shown;

[0047] Figure 4 A constrained planar diagram obtained by converting a depositional phase mode diagram according to an embodiment of the present disclosure is shown;

[0048] Figure 5 A sandstone thickness map obtained according to an embodiment of the present disclosure is shown, based on a geostatistical stochastic simulation algorithm and constrained plane constraints.

[0049] Figure 6 A schematic flowchart of a method for drawing a sedimentary facies pattern diagram of a preset region using sandstone curve morphology and sandstone thickness, according to an embodiment of the present disclosure, is shown.

[0050] Figure 7 A block diagram of a sandstone thickness mapping system according to an embodiment of the present disclosure is shown;

[0051] Figure 8 A schematic diagram of an electronic device structure according to an embodiment of the present disclosure is shown. Detailed Implementation

[0052] To make the objectives, technical solutions, and advantages of the embodiments of this disclosure clearer, the technical solutions of the embodiments of this disclosure will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this disclosure, and not all embodiments. Based on the embodiments of this disclosure, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this disclosure.

[0053] The same reference numerals in the accompanying drawings denote the same or similar elements, components, or parts, and therefore, repeated descriptions of the same or similar elements, components, or parts may be omitted below. It should also be understood that although terms such as first, second, third, etc., indicating numbers may be used herein to describe various devices, elements, components, or parts, these devices, elements, components, or parts should not be limited by these terms. That is, these terms are only used to distinguish one from another. For example, a first device may also be referred to as a second device, without departing from the essential technical solution of the invention. Furthermore, the terms "and / or" and "and / or" refer to all combinations including any one or more of the listed items.

[0054] Please see Figure 1 , Figure 1 This is a schematic diagram of a method for drawing sandstone thickness maps provided by the present invention, as shown below. Figure 1 As shown, the method includes:

[0055] S11. Obtain logging data for each well point in the preset area, and determine the lithological characteristics of sandstone and mudstone based on the logging data.

[0056] In this embodiment, well logging data refers to formation information data obtained through well logging technology, which reflects the physical and chemical properties of the formation. Well logging technology is a key technology in oil exploration and development. It uses specialized measuring instruments to measure various physical parameters of the formation along the wellbore, such as resistivity, acoustic velocity, and natural gamma, thereby inferring key information such as rock type, porosity, and permeability of the formation.

[0057] In this embodiment, the sandstone and mudstone include: Mudstone: mainly composed of clay minerals, with possibly very small amounts of sand and silt in some areas. Clay minerals mainly include kaolinite, illite, montmorillonite, etc. Sandstone: mainly composed of quartz sand grains, feldspar, and other clastic rock particles. The content and size of quartz, feldspar, and other mineral particles determine the texture and properties of the sandstone.

[0058] In this embodiment, the lithological characteristics of sandstone and mudstone can be reflected through logging curves during the well logging process. For example: Apparent resistivity: Mudstone has a relatively low apparent resistivity value, typically in the range of 8–15 Ωm. Sandstone has a relatively high apparent resistivity value, which gradually increases with increasing grain diameter and decreasing clay content. Natural gamma: Mudstone has a high natural gamma value, typically in the range of 120–140 API. Sandstone has a relatively low natural gamma value, which gradually decreases with increasing grain diameter. Density: Mudstone has a relatively low density value, typically in the range of 2.35–2.45 g / cm³. Sandstone has a relatively high density value, which gradually increases with increasing grain diameter and decreasing clay content.

[0059] Specifically, the natural gamma logging curve data of each sampling point in each well point of the preset area are obtained; the mud content of the corresponding sampling point is calculated based on the maximum and minimum values ​​of the natural gamma logging curve data, which is used as the lithological characteristic data of sandstone and mudstone.

[0060] The method for calculating clay content using natural gamma rays is based on the relationship between gamma ray intensity and rock radioactivity. A commonly used method and steps for calculating clay content using natural gamma rays are as follows:

[0061] The maximum (GRmax) and minimum (GRmin) values ​​of the natural gamma logging curve for the target formation are determined. This method assumes that the minimum natural gamma value (GRmin) represents pure sandstone, while the minimum natural gamma value (GRmax) represents 100% pure mudstone.

[0062] Calculate the clay content using the following formula:

[0063] Vsh represents the mud content.

[0064] GR represents the mud content of the sampling point data. A cutoff value for mud content (Vsh) is given. For example, if Vsh is set to 50%, then Vsh > 50% is mudstone, and Vsh < 50% is sandstone. This classification is applied to the data of each mud content sampling point. Of course, Vsh can also be set to 40%, which can be determined according to the local mudstone and sandstone conditions. This scheme does not impose any special restrictions on this.

[0065] like Figure 2 As shown, the collection and analysis of single-well and sedimentary facies data: well logging data, core samples, stratigraphic correlation and sedimentary facies data were collected in the study area, and the lithology of sandstone and mudstone was interpreted according to the stratigraphy or sand group.

[0066] S12. Determine the sandstone curve morphology and sandstone thickness at each well point based on the sandstone and mudstone lithological characteristics data, and draw the sedimentary facies model diagram of the preset area based on the sandstone curve morphology and sandstone thickness.

[0067] In this embodiment, the sandstone curve morphology is determined using lithological characteristic data of sandstone and mudstone. For example, if the clay content gradually increases in a certain well section, the natural gamma ray curve may exhibit a bell-shaped shape, meaning the curve value gradually rises from a low value to a high value. Correspondingly, the sandstone curve may exhibit the opposite shape, meaning the amplitude gradually decreases. This usually indicates that the sandstone layer becomes finer upwards or that the clay content gradually increases, resulting in a bell-shaped sandstone curve. If the clay content remains stable or does not change much in a certain well section, the natural gamma ray curve may exhibit a box-shaped shape, meaning the curve value fluctuates within a certain range but remains generally stable. In this case, the sandstone curve may exhibit a relatively stable shape, indicating that the sandstone layer is relatively uniform within that well section, resulting in a box-shaped sandstone curve. Conversely, if the clay content gradually decreases in a certain well section, the natural gamma ray curve may exhibit a funnel-shaped shape. Correspondingly, the sandstone curve may gradually rise, indicating that the sandstone layer becomes coarser downwards or that the clay content gradually decreases, resulting in a funnel-shaped sandstone curve. If the mud content changes rapidly within the well section, the sandstone curve will be finger-shaped or tooth-shaped. A finger-shaped curve may indicate mud interlayers or lenses in the sandstone layer, while a tooth-shaped curve may indicate the alternation of mud and sandstone.

[0068] In this embodiment, sandstone thickness can be determined using natural gamma data, resistivity, or acoustic waves. Sandstone typically has a lower natural gamma value than mudstone. Therefore, on well logging curves, sandstone layers usually appear as regions with low natural gamma values. By measuring the thickness of these regions, the thickness of the sandstone can be preliminarily estimated. Sandstone typically has a higher resistivity than mudstone. On well logging curves, sandstone layers typically appear as regions with high resistivity. Similarly, by measuring the thickness of these regions, information about the sandstone thickness can also be obtained. Acoustic waves typically propagate faster in sandstone than in mudstone. Therefore, on well logging curves, sandstone layers typically appear as regions with high acoustic velocities.

[0069] In this embodiment, taking clay content as an example, after determining the clay content threshold, the sandstone thickness can be calculated by measuring the length of the well section corresponding to the sandstone layer on the logging curve. This usually requires the use of specialized logging interpretation software or manual calculation.

[0070] In this embodiment, after determining the sandstone curve morphology and sandstone thickness, a sedimentary facies model diagram was drawn. A sedimentary facies model diagram is a comprehensive and highly generalized graphical representation of sedimentary processes and their products in a specific sedimentary environment. The following are some common sedimentary facies model diagrams and their brief descriptions: 1. Glacial Facies Model Diagram: This diagram mainly shows the characteristics of sediments under glacial action, including clastic material transported by glaciers and sediments formed after glacial melting. These sediments are typically poorly sorted, have low roundness, and often contain glacial till unique to glaciers. 2. Piedmont-Alluvial Facies Model Diagram: This diagram mainly shows sedimentary deposits formed by temporary flood flows at the valley outlet, i.e., alluvial fans. Alluvial fan sediments typically have a gradually decreasing grain size from the fan apex to the fan edge, and the fan apex often contains larger gravels and boulders. 3. Fluvial Facies Model Diagram: Fluvial facies model diagrams can be divided into braided river, meandering river, and reticulated river sedimentary models, depending on the river type. Braided River Depositional Model: This model illustrates the sedimentary characteristics of braided rivers, which have multiple channels, bifurcations, and convergences. Sediments are typically dominated by sand and gravel, with significant variations in grain size. Meandering River Depositional Model: This model shows the stable riverbed of a meandering river, characterized by channels, oxbow lakes, and natural dikes. Sediment grain size gradually decreases from the riverbed to the floodplain, and wavy bedding and cross-bedding are commonly observed. Reticulated River Depositional Model: This model illustrates the meandering, numerous, and fixed channels of a reticulated river. Sediments are typically dominated by fine-grained clastic rocks, and horizontal and massive bedding are commonly observed. IV. Lacustrine Facies Model: Lacustrine facies models can be categorized based on lake type, including terrigenous clastic lakes, biogenic lakes, and lacustrine-marsh lakes. The sedimentary characteristics of these lakes differ, but they typically include fine-grained clastic rocks such as mudstone and siltstone, and often exhibit horizontal and wavy bedding. V. Deltaic Facies Model Diagram: This diagram primarily illustrates the sedimentary deposition system formed at the confluence of rivers and oceans (or lakes). Deltaic sediments typically exhibit a gradual decrease in grain size from the deltaic plain to the prodelta and often contain sedimentary structures such as mouth bars and distal sand bars. VI. Other Sedimentary Facies Model Diagrams: In addition to the common sedimentary facies model diagrams mentioned above, there are also models for barrier islands, tidal flats, lagoons, estuarine bays, marine facies, gravity flow sediments, carbonate sedimentary facies, and reef facies. These models respectively demonstrate the sedimentary characteristics and sedimentary structures under different sedimentary environments.

[0071] like Figure 3 As shown, based on the sandstone, curve morphology, and thickness data interpreted from single wells, the dominant sedimentary facies of the target intervals are determined, and sedimentary facies are divided for each target interval in each single well. Sedimentary facies pattern diagrams are drawn by combining data from multiple wells. Table 1 shows the statistical data on sandstone thickness from single wells in the preset area.

[0072] hashtag target segment X coordinate Y coordinate Sandstone thickness Sword and Fairy 4 N1-I 16609960.30 4200317.50 16.2 Xian 14 N1-I 16610243.20 4198598.00 31.8 Xian 103 N1-I 16610921.20 4199433.10 29.05 Immortal 105 N1-I 16611018.40 4198691.70 20.5 Xianbei 1 N1-I 16610698.80 4203187.20 23.92 Xiandong 7 N1-I 16613063.80 4199227.20 35.82 Xiandong 8 N1-I 16611761.50 4199436.50 24.7 Xiandong 3 N1-I 16612141.70 4197929.20 34.2 … … … … …

[0073] The sedimentary facies model of the preset area was obtained by plotting the sandstone thickness at each well point, such as... Figure 3 As shown, the sedimentary facies model diagram is a braided river sedimentary model diagram.

[0074] S13. Within the plane of the three-dimensional geological framework model, convert the sedimentary facies pattern diagram to a constrained plane used to determine the trend changes of sedimentary facies zones.

[0075] In this embodiment, the sandstone thickness of the target layer at each well point can be obtained. The maximum thickness of the sandstone is used as the bottom surface of the framework model, and the minimum thickness is used as the top surface of the framework model to generate a three-dimensional geological framework model. Specifically, the sandstone thickness data of the target layer for each well is statistically obtained, and a three-dimensional framework model of the study area is established: the sandstone thickness data of the target layer for each well is statistically obtained, and according to the size of the mapping range, the maximum thickness is used as the bottom surface of the three-dimensional framework model, and the minimum thickness is used as the top surface of the three-dimensional framework model to establish a three-dimensional geological framework model for subsequent geostatistical sandstone thickness modeling.

[0076] In this embodiment, the constraint plane of the deposition pattern diagram is obtained by editing based on the planar range of the three-dimensional frame model. For example... Figure 4 As shown, the number of facies zones in the sedimentary model planar map is counted. For example, if there are 3 facies zones, then facies zone 1 is assigned a value of 1, facies zone 2 is assigned a value of 2, and facies zone 3 is assigned a value of 3. According to this method, the sedimentary facies planar map is converted into a constrained planar map with corresponding values. This map reflects the changes in the trend of sedimentary facies zones.

[0077] Specifically, the facies zones in the sedimentary facies pattern diagram are obtained and numbered; within the plane of the three-dimensional geological framework model, a constraint plane is generated based on the sedimentary facies pattern diagram and the numbering.

[0078] S14. Based on the geostatistical stochastic simulation algorithm, generate planar sandstone thickness data by constraining the plane, and draw a sandstone thickness map of the preset area based on the planar sandstone thickness data.

[0079] In this embodiment, as Figure 5 As shown, geostatistical stochastic simulations constrained by sedimentary facies were used to obtain planar sandstone thickness data: within a three-dimensional geological framework model, geostatistical stochastic simulation algorithms were employed, combined with planar constraints from sedimentary facies models, to predict the distribution of sandstone thickness. In this process, the aforementioned sedimentary facies models served as constraints, guiding the analysis of spatial variability in sandstone thickness. In areas where actual sandstone thickness data was unavailable, the predicted sandstone thickness conformed to the characteristics of sedimentary facies distribution.

[0080] In this embodiment, based on the geostatistical stochastic simulation algorithm, the constrained plane is used as a constraint condition to generate sandstone thickness distribution, sandstone thickness randomness, and sandstone thickness heterogeneity as planar sandstone thickness data.

[0081] In this embodiment, a sandstone thickness contour map of a preset area is drawn based on planar sandstone thickness data. Contour maps, as a graphical representation, are easily understood by non-experts. This facilitates effective communication and exchange between geologists and experts from other fields, policymakers, and the public. The sandstone thickness contour map visually displays the spatial distribution of sandstone thickness through changes in color, lines, or symbols. This visualization method enables geologists and engineers to quickly understand the changing trends and distribution characteristics of sandstone thickness.

[0082] In this embodiment, the present disclosure determines the lithological characteristics of each well point using logging data, and statistically obtains the sandstone curve morphology and sandstone thickness at each well point based on the lithological characteristics data. A sedimentary facies model diagram of the area is then drawn. This sedimentary facies model diagram is converted into a constraint plane within the three-dimensional geological framework model plane, allowing the user to determine the trend changes of sedimentary facies zones. Using this constraint plane, a geostatistical stochastic simulation algorithm is used to predict the sandstone thickness distribution and determine the planar sandstone thickness data. Finally, a sandstone thickness map is drawn using this planar sandstone thickness data. This solution, by combining geostatistical principles and sedimentary facies models, achieves high-precision prediction and visualization of sandstone thickness.

[0083] like Figure 6 As shown, step S12 involves drawing a sedimentary facies model diagram of the preset area using sandstone curve morphology and sandstone thickness, including:

[0084] S21. Based on the sandstone curve morphology and sandstone thickness, determine the dominant sedimentary facies of the target layer at each well point.

[0085] In this embodiment, the grain size variation, lithological assemblage, and sedimentary environment of sediments are determined based on the anomalous amplitude, shape (e.g., box-shaped, bell-shaped, funnel-shaped), and smoothness of the spontaneous potential curve. The conductivity, porosity, and permeability of the rocks are determined based on the trend of the apparent resistivity curve. By comparing the sandstone thickness at different well points, areas with well-developed sandstone and areas with thinner sandstone are identified. The relationship between sandstone thickness variation and sedimentary facies is analyzed using sedimentary facies model diagrams. Finally, by combining the analysis results of curve morphology and sandstone thickness, as well as geological background and sedimentary environment factors, the dominant sedimentary facies at each well point's target layer is comprehensively determined.

[0086] S22. Based on the dominant sedimentary facies at each well point target layer, the corresponding well point target layers are divided into sedimentary facies, and a sedimentary facies model diagram of the preset area is drawn by integrating the sedimentary facies of each well point.

[0087] In this embodiment, the identified sedimentary facies are classified according to their sedimentary environment and characteristics, such as terrestrial facies, marine facies, and transitional marine-terrestrial facies. Within each facies group, specific sedimentary facies are further subdivided based on differences in secondary environment and sediment characteristics, such as fluvial facies, deltaic facies, and lacustrine facies. For complex sedimentary systems, sedimentary subfacies and microfacies can also be further subdivided.

[0088] Please see Figure 7 , Figure 7 This invention provides a sandstone thickness map drawing system, which includes: a data analysis module 12, a graphics drawing module 13, a constraint plane conversion module 14, and a sandstone thickness generation module 15.

[0089] In this embodiment, the data analysis module 12 is used to acquire logging data of each well point in a preset area; and to determine the lithological characteristics of sandstone and mudstone based on the logging data.

[0090] In this embodiment, the graphics drawing module 13 is used to determine the sandstone curve shape and sandstone thickness of each well point based on the sandstone and mudstone lithological characteristic data, and to draw a sedimentary facies pattern diagram of a preset area based on the sandstone curve shape and sandstone thickness.

[0091] In this embodiment, the constraint plane conversion module 14 is used to convert the sedimentary facies pattern diagram to a constraint plane for determining the trend changes of sedimentary facies zones within the plane range of the three-dimensional geological framework model.

[0092] In this embodiment, the sandstone thickness generation module 15 is used to generate planar sandstone thickness data by constraining the plane based on a geostatistical random simulation algorithm.

[0093] In this embodiment, the graphics drawing module 13 is also used to draw a sandstone thickness map of a preset area based on the sandstone thickness data of the plane.

[0094] In this embodiment, the data acquisition module 11 is specifically used to acquire the natural gamma logging curve data of each sampling point in each well point of the preset area.

[0095] In this embodiment, the data analysis module 12 is specifically used to calculate the mud content of the corresponding sampling point based on the maximum and minimum values ​​of the natural gamma logging curve data, as lithological characteristic data of sandstone and mudstone.

[0096] In this embodiment, the drawing system further includes: a frame model generation module, used to obtain the sandstone thickness of the target layer at each well point; using the maximum thickness of the sandstone thickness as the bottom surface of the frame model and the minimum thickness of the sandstone thickness as the top surface of the frame model to generate a three-dimensional geological frame model.

[0097] In this embodiment, the constraint plane conversion module 14 is specifically used to obtain the facies zones in the sedimentary facies pattern diagram and number the facies zones; within the plane range of the three-dimensional geological framework model, a constraint plane is generated based on the sedimentary facies pattern diagram and the numbering.

[0098] In this embodiment, the sandstone thickness generation module 15 is specifically used to generate sandstone thickness distribution, sandstone thickness randomness, and sandstone thickness heterogeneity as planar sandstone thickness data by using a geostatistical random simulation algorithm with a constrained plane as a constraint condition.

[0099] In this embodiment, the graphics drawing module 13 is specifically used to determine the dominant sedimentary facies of the target layer of each well point based on the sandstone curve morphology and sandstone thickness; to divide the corresponding target layer of the well point into sedimentary facies based on the dominant sedimentary facies of the target layer of each well point; and to draw a sedimentary facies pattern diagram of the preset area by integrating the sedimentary facies of each well point.

[0100] In this embodiment, the graphics drawing module 13 is specifically used to draw a sandstone thickness contour map of a preset area based on the sandstone thickness data of the plane.

[0101] like Figure 8 As shown, this embodiment of the invention provides an electronic device, including a processor 1110, a communication interface 1120, a memory 1130, and a communication bus 1140, wherein the processor 1110, the communication interface 1120, and the memory 1130 communicate with each other through the communication bus 1140.

[0102] Memory 1130 is used to store computer programs;

[0103] The processor 1110, when executing the program stored in the memory 1130, implements any of the above drawing methods.

[0104] The electronic device provided in this embodiment of the invention includes a processor 1110 that executes a program stored in a memory 1130 to obtain well logging data for each well point in a preset area, and determines sandstone and mudstone lithological characteristic data based on the well logging data; determines the sandstone curve morphology and sandstone thickness for each well point based on the sandstone and mudstone lithological characteristic data, and draws a sedimentary facies pattern diagram for the preset area based on the sandstone curve morphology and sandstone thickness; converts the sedimentary facies pattern diagram to a constraint plane for determining the trend change of sedimentary facies zones within the three-dimensional geological framework model plane range; generates planar sandstone thickness data based on a geostatistical stochastic simulation algorithm through constraints on the constraint plane, and draws a sandstone thickness map for the preset area based on the planar sandstone thickness data.

[0105] The communication bus 1140 mentioned in the above electronic device can be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus, etc. This communication bus 1140 can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, it is shown in the figure with only one thick line, but this does not indicate that there is only one bus or one type of bus.

[0106] The communication interface 1120 is used for communication between the above-mentioned electronic device and other devices.

[0107] The memory 1130 may include random access memory (RAM) or non-volatile memory, such as at least one disk storage device. Optionally, the memory 1130 may also be at least one storage device located remotely from the aforementioned processor 1110.

[0108] The processor 1110 mentioned above can be a general-purpose processor 1110, including a central processing unit (CPU), a network processor (NP), etc.; it can also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components.

[0109] This invention provides a computer-readable storage medium storing one or more programs that can be executed by one or more processors 1110 to implement the drawing method of any of the above embodiments.

[0110] In the above embodiments, implementation can be achieved, in whole or in part, through software, hardware, firmware, or any combination thereof. When implemented in software, it can be implemented, in whole or in part, as a computer program product. A computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a computer, all or part of the flow or function according to the embodiments of the present invention is generated. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions can be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another. For example, computer instructions can be transmitted from one website, computer, server, or data center to another website, computer, server, or data center via wired (e.g., coaxial cable, fiber optic, digital subscriber line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.) means. The computer-readable storage medium can be any available medium that a computer can access or a data storage device such as a server or data center that integrates one or more available media. The available medium can be a magnetic medium (e.g., floppy disk, hard disk, magnetic tape), an optical medium (e.g., DVD), or a semiconductor medium (e.g., a solid-state drive (SSD)).

[0111] Although the present disclosure has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present disclosure.

Claims

1. A method for drawing sandstone thickness maps based on sedimentary facies constraints, characterized in that, The drawing method includes: Obtain logging data from each well point in a preset area, and determine the lithological characteristics of sandstone and mudstone based on the logging data; Based on the lithological characteristics data of the sandstone and mudstone, the sandstone curve morphology and sandstone thickness of each well point are determined, and the sedimentary facies model diagram of the preset area is drawn based on the sandstone curve morphology and sandstone thickness. Within the three-dimensional geological framework model plane, the sedimentary facies pattern diagram is transformed into a constrained plane used to determine the trend changes of sedimentary facies zones; Based on a geostatistical stochastic simulation algorithm, planar sandstone thickness data is generated by constraining the constrained plane, and a sandstone thickness map of the preset area is drawn based on the planar sandstone thickness data.

2. The method according to claim 1, characterized in that, The step of acquiring logging data from each well point in a preset area and determining sandstone and mudstone lithological characteristics based on the logging data includes: Obtain the natural gamma logging curve data of each sampling point in each well point of the preset area; The mud content of the corresponding sampling point is calculated based on the maximum and minimum values ​​of the natural gamma logging curve data, and is used as the lithological characteristic data of the sandstone and mudstone.

3. The method according to claim 1, characterized in that, The method further includes: Obtain the sandstone thickness of the target layer at each well point; The maximum thickness of the sandstone is used as the bottom surface of the frame model, and the minimum thickness of the sandstone is used as the top surface of the frame model to generate a three-dimensional geological frame model.

4. The method according to claim 3, characterized in that, The process of converting the sedimentary facies pattern diagram to a constrained plane for determining the trend changes of sedimentary facies zones within the three-dimensional geological framework model plane includes: Obtain the facies zones in the sedimentary facies pattern diagram and number the facies zones; Within the plane of the three-dimensional geological framework model, a constraint plane is generated based on the sedimentary facies pattern diagram and the numbering.

5. The method according to claim 1, characterized in that, The method based on geostatistical stochastic simulation, which generates planar sandstone thickness data through the constraint plane, includes: Based on the aforementioned geostatistical stochastic simulation algorithm, the constrained plane is used as a constraint condition to generate the sandstone thickness distribution, sandstone thickness randomness, and sandstone thickness heterogeneity, which are then used as the planar sandstone thickness data.

6. The method according to claim 1, characterized in that, The step of drawing a sedimentary facies model diagram of the preset area using the sandstone curve morphology and sandstone thickness includes: Based on the sandstone curve morphology and sandstone thickness, the dominant sedimentary facies of the target layer at each well point are determined. Based on the dominant sedimentary facies of each well point target layer, the corresponding well point target layers are divided into sedimentary facies, and a sedimentary facies pattern diagram of the preset area is drawn by combining the sedimentary facies of each well point.

7. The method according to any one of claims 1 to 6, characterized in that, The step of drawing a sandstone thickness map of the preset area based on the planar sandstone thickness data includes: Based on the plane sandstone thickness data, draw a contour map of sandstone thickness in the preset area.

8. A sandstone thickness map drawing system based on sedimentary facies constraints, characterized in that, The drawing system includes: The data analysis module is used to acquire logging data from various well points in a preset area; and to determine the lithological characteristics of sandstone and mudstone based on the logging data. The graphics drawing module is used to determine the sandstone curve shape and sandstone thickness of each well point based on the sandstone and mudstone lithological characteristic data, and to draw the sedimentary facies pattern diagram of the preset area based on the sandstone curve shape and sandstone thickness. The constraint plane conversion module is used to convert the sedimentary facies pattern diagram to a constraint plane for determining the trend changes of sedimentary facies zones within the plane range of the three-dimensional geological framework model. The sandstone thickness generation module is used to generate planar sandstone thickness data based on a geostatistical stochastic simulation algorithm and constrained by the constrained plane. The graphics drawing module is also used to draw a sandstone thickness map of the preset area based on the plane sandstone thickness data.

9. An electronic device, characterized in that, include: processor; as well as A memory storing computer-executable instructions, which, when executed, cause the processor to perform the method according to any one of claims 1-7.

10. A computer storage medium, characterized in that, in, The computer storage medium stores one or more programs that, when executed by a processor, implement the method of any one of claims 1-7.

Images