Method, system and medium for converting sedimentary simulation results into three-dimensional geological models

Abstract

The present invention provides a method for converting sedimentation simulation results into a three-dimensional geological model, including: S1. Obtaining a time domain model (nx, ny, nk); S2. Creating a three-dimensional matrix (nx, ny, nz); S3. ik = nk; S4. Comparing the delta terrain elevation data at the current moment and the previous moment, judging whether (ix, iy) is deposited; S5. If (ix, iy) is deposited, then deposit; if (ix, iy) is not deposited, then erosion; S6. Judging whether ix reaches a maximum value, if yes, jump to S7; if no, assign ix = i x + 1 and jump to S4; S7. Judging whether iy reaches a maximum value, if yes, jump to S8; if no, assign iy = i y + 1 and jump to S4; S8. Judging whether the time step number is traversed, if yes, save the result; if no, assign ik = ik-1 and jump to S4. The present invention can convert two-dimensional sedimentation simulation results into a three-dimensional geological model.

Description

Technical Field

[0001] The present invention relates to the technical field of sedimentation simulation and three-dimensional geological model, and in particular to a method, system and medium for converting sedimentation simulation results into a three-dimensional geological model. Background Art

[0002] Three-dimensional geological modeling is a new technology for geological research that uses computer technology to combine spatial information management, geological interpretation, spatial analysis and prediction, geostatistics, entity content analysis, and graphic visualization tools in a three-dimensional environment. It actually combines geological, logging, geophysical data and various interpretation results or conceptual models to generate a three-dimensional quantitative stochastic model.

[0003] Sedimentation simulation is an important experimental means and technical method for sedimentology theory research, which can be divided into numerical simulation and physical simulation. Physical simulation is an indoor simulation of the physical process of sedimentation. By simulating the sedimentation conditions at that time, the sedimentation process of natural sediments is restored in the laboratory.

[0004] The output results of current sedimentation simulation software are mostly in two-dimensional form and do not have the characteristics of spatial three-dimensional data volumes. Three-dimensional data can better and intuitively display simulation results or prediction results. Summary of the invention

[0005] The present invention provides a method, system and medium for converting sedimentation simulation results into a three-dimensional geological model, and converts the two-dimensional sedimentation simulation results into a three-dimensional geological model to better and intuitively display the simulation results or prediction results.

[0006] A first aspect of an embodiment of this specification discloses a method for converting sedimentation simulation results into a three-dimensional geological model, comprising:

[0007] S1. Obtain a known time domain model (nx, ny, nk), where nx represents the number of sampling points in the x direction; ny represents the number of sampling points in the y direction, and nk represents the number of time sampling points;

[0008] S2. Create a three-dimensional matrix (nx, ny, nz) and assign the value of the three-dimensional matrix to null for saving the depth domain geological model; nz represents the number of sampling points in the depth direction converted to the geological model;

[0009] S3. Assignment: Set a variable ik as the time pointer, and set variables ix and iy as the x and y direction pointers respectively; the initial value of ik is set to nk, and the value range is [1,nk]; the initial values ​​of ix and iy are set to 1, and the value ranges are [1,nx], [1,ny] respectively;

[0010] S4. Compare the delta terrain elevation data at the current moment (ix, iy, ik) with the delta terrain elevation data at the previous moment (ix, iy, ik-1) to determine whether there is deposition at (ix, iy);

[0011] S5. If there is deposition at (ix, iy), then a deposition process is performed: the three-dimensional matrix (ix, iy, ik-1~ik) is assigned with sediment grain size data; if there is no deposition at (ix, iy), then an erosion process is performed: the value at the three-dimensional matrix (ix, iy, ik-1~ik) is modified to null to achieve an erosion effect;

[0012] S6. Determine whether ix has reached the maximum value. If so, jump to S7; if not, assign ix = ix + 1 and jump to S4;

[0013] S7. Determine whether iy reaches the maximum value. If so, jump to S8; if not, assign: iy = iy + 1, and jump to S4;

[0014] S8. Determine whether the traversal of the time step is complete. If so, save the result to obtain the depth domain geological model; if not, assign: ik = ik-1, and jump to S4.

[0015] A second aspect of the embodiments of this specification discloses a system for converting sedimentation simulation results into a three-dimensional geological model, including:

[0016] The time domain model acquisition module is used to obtain a known time domain model (nx, ny, nk), where nx represents the number of sampling points in the x direction; ny represents the number of sampling points in the y direction; and nk represents the number of time sampling points.

[0017] The three-dimensional matrix creation module is used to create a three-dimensional matrix (nx, ny, nz), whose value is null, and is used to save the depth domain geological model; nz represents the number of sampling points in the depth direction converted into a geological or geophysical model;

[0018] The assignment module is used to assign values: create pointers ik, ix, and iy, where the initial value of ik is set to nk and the value range is set to [1, nk]; the initial values ​​of ix and iy are equal to 1, and the value ranges are [1, nx] and [1, ny] respectively;

[0019] A comparison module is used to compare the delta terrain elevation data at the current moment (ix, iy, ik) with the delta terrain elevation data at the previous moment (ix, iy, ik-1) to determine whether there is deposition at (ix, iy);

[0020] The deposition or erosion module is used to perform deposition processing if deposition occurs at (ix, iy): assign sediment grain size data to the three-dimensional matrix (ix, iy, ik-1~ik); if no deposition occurs at (ix, iy), erosion processing is performed: assign null to the three-dimensional matrix (ix, iy, ik-1~ik) to achieve an erosion effect;

[0021] The first judgment module is used to judge whether ix reaches the maximum value. If so, it jumps to the second judgment module; if not, it assigns: ix=ix+1 and jumps to the comparison module;

[0022] The second judgment module is used to judge whether iy reaches the maximum value. If so, it jumps to the third judgment module; if not, it assigns: iy=iy+1 and jumps to the comparison module;

[0023] The third judgment module is used to judge whether the time steps have been traversed. If so, the result is saved to obtain the depth domain geological model; if not, the value is assigned: ik = ik-1, and jump to the comparison module.

[0024] A third aspect of the embodiments of this specification discloses a computer-readable storage medium, which stores computer instructions. When a computer reads the computer instructions, the computer executes the method of converting sedimentation simulation results into a three-dimensional geological model as described above.

[0025] The embodiments of this specification can at least achieve the following beneficial effects:

[0026] First, the sedimentation simulation software is used to simulate the sedimentation, and the sedimentation simulation results are obtained, that is, the time domain model. Then, according to the calculated intermediate delta morphological elevation data, combined with the grain size distribution data, the erosion or deposition is calculated from old to new in time, and then the corresponding data processing is performed to obtain the three-dimensional depth domain lithology model. After converting the lithology model according to the corresponding relationship, the three-dimensional depth domain velocity model can be obtained. The purpose of converting the two-dimensional sedimentation simulation results into a three-dimensional geological model is achieved, and the simulation results or prediction results can be better and intuitively displayed. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0028] Figure 1 The figure is a flow chart of the method for converting sedimentation simulation results into a three-dimensional geological model involved in the present invention.

[0029] Figure 2 It is a schematic diagram showing the simulation results of the Delft3d hydrodynamic simulation software involved in the present invention.

[0030] Figure 3 It is a schematic diagram of the delta lithology simulation involved in the present invention. DETAILED DESCRIPTION

[0031] In the following, only some exemplary embodiments are briefly described. As those skilled in the art will appreciate, the described embodiments may be modified in various ways without departing from the spirit or scope of the embodiments of the present invention. Therefore, the drawings and descriptions are considered to be exemplary and non-restrictive in nature.

[0032] The disclosure below provides many different embodiments or examples to implement different structures of the embodiments of the present invention. In order to simplify the disclosure of the embodiments of the present invention, the parts and settings of specific examples are described below. Of course, they are only examples, and the purpose is not to limit the embodiments of the present invention. In addition, the embodiments of the present invention can repeat reference numbers and / or reference letters in different examples, and this repetition is for the purpose of simplification and clarity, which itself does not indicate the relationship between the various embodiments and / or settings discussed.

[0033] The embodiments of the present invention are described in detail below with reference to the accompanying drawings.

[0034] The basic idea of ​​geological model generation is to gradually build up the surface sedimentation status at different times from old to new in time. The specific idea is to use satellite photos as an example. Assume that there are n satellite photos at different times, and assume that these satellite photos can be used to extract elevation data and surface sedimentary lithology data that meet the accuracy requirements. Then, according to time, the latest photo is numbered 1 and placed at the top layer, and the other photos are placed downward in sequence, and the old photos are numbered n and placed at the bottom; then the deposition or erosion relationship is calculated based on the elevation numbers in these n satellite photos, and finally, combined with the surface sedimentary lithology data, it is possible to build a three-dimensional data body of lithology deposited underground from the initial photo No. 1 to the photo No. n in time; the larger the number of photos n, the higher the accuracy of the established geological model. The detailed steps are as follows:

[0035] like Figure 1 As shown, the first aspect of the embodiment of this specification discloses a method for converting sedimentation simulation results into a three-dimensional geological model, including:

[0036] S1. Obtain a known time domain model (nx, ny, nk), where nx represents the number of sampling points in the x direction; ny represents the number of sampling points in the y direction; and nk represents the number of time sampling points;

[0037] S2. Create a three-dimensional matrix (nx, ny, nz), whose value is null, for saving the depth domain geological model; nz represents the number of sampling points in the depth direction converted into the geological model;

[0038] S3. Assignment: define pointers ik, ix, and iy; the initial value of ik (represents a pointer, traversing from the head to the tail of nk. The "i" in "ik" represents the first letter of index.) is set to nk (the initial value of ik is the initial deposition time, that is, it points to the bottom of the 3D time domain model), and the value range is [1, nk]; the initial values ​​of ix and iy are set to 1, and the value ranges are [1, nx] and [1, ny] respectively, including two nested loops.

[0039] for(i=1;i<=nx;i++)

[0040] for(int j=1;j<=ny;j++)

[0041] {…}

[0042] Its function is to traverse the entire plane (calculation work area); all points that make up the plane are traversed.

[0043] ik is the age pointer, which points to the number of the latest photo from the oldest. ix is ​​the x-direction pointer, which points to the nx end from the nx head; iy is the y-direction pointer, which points to the ny end from the ny head. Step 1: ik points to number 1, then ix and iy are nested in a loop, and when ik equals 1, the entire xy plane is traversed, the deposition or erosion relationship is processed, and the geological model is refreshed. ik = ik-1, and the next cycle is performed. )

[0044] S4. Compare the delta terrain elevation data at the current moment (ix, iy, ik) with the delta terrain elevation data at the previous moment (ix, iy, ik-1) to determine whether there is deposition at (ix, iy);

[0045] S5. If there is deposition at (ix, iy), then deposition is performed: the three-dimensional matrix (ix, iy, ik-1~ik) is assigned with sediment grain size data; if there is no deposition at (ix, iy), then erosion is performed: the value at the three-dimensional matrix (ix, iy, ik-1~ik) is modified to null to achieve the erosion effect;

[0046] S6. Determine whether ix has reached the maximum value. If so, jump to S7; if not, assign ix = ix + 1 and jump to S4;

[0047] S7. Determine whether iy reaches the maximum value. If so, jump to S8; if not, assign: iy = iy + 1, and jump to S4;

[0048] S8. Determine whether the time steps have been traversed. If so, save the result to obtain the depth domain geological model; if not, assign the value: ik = ik-1, and jump to S4.

[0049] A second aspect of the embodiments of this specification discloses a system for converting sedimentation simulation results into a three-dimensional geological model, including:

[0050] The time domain model acquisition module is used to obtain a known time domain model (nx, ny, nk), where nx represents the number of sampling points in the x direction; ny represents the number of sampling points in the y direction; and nk represents the number of time sampling points.

[0051] The three-dimensional matrix creation module is used to create a three-dimensional matrix (nx, ny, nz), whose value is null, and is used to save the depth domain geological model; nz represents the number of sampling points in the depth direction converted into a geological or geophysical model;

[0052] Assignment module, define pointers ik, ix, iy; the initial value of ik is set to nk, and the value range is [1, nk]; the initial values ​​of ix and iy are set to 1, and the value ranges are [1, nx], [1, ny] respectively;

[0053] A comparison module is used to compare the delta terrain elevation data at the current moment (ix, iy, ik) with the delta terrain elevation data at the previous moment (ix, iy, ik) to determine whether there is deposition at (ix, iy);

[0054] The deposition or erosion module is used to perform deposition if deposition occurs at (ix, iy): assign sediment grain size data to the three-dimensional matrix (ix, iy, ik-1~ik); if no deposition occurs at (ix, iy), erosion is performed: modify the value at the three-dimensional matrix (ix, iy, ik-1~ik) to null to achieve an erosion effect;

[0055] The first judgment module is used to judge whether ix reaches the maximum value. If so, it jumps to the second judgment module; if not, it assigns: ix=ix+1 and jumps to the comparison module;

[0056] The second judgment module is used to judge whether iy reaches the maximum value. If so, it jumps to the third judgment module; if not, it assigns: iy=iy+1 and jumps to the comparison module;

[0057] The third judgment module is used to judge whether the time step pointer ik points to the initial moment (traversal is completed). If so, the result is saved to obtain the depth domain geological model (three-dimensional geological model); if not, the value is assigned: ik = ik-1, and jump to the comparison module.

[0058] A third aspect of the embodiments of this specification discloses a computer-readable storage medium, which stores computer instructions. When a computer reads the computer instructions, the computer executes the method of converting sedimentation simulation results into a three-dimensional geological model as described above.

[0059] The technical concept of the present invention is as follows:

[0060] The deposition simulation results can be obtained through existing deposition simulation software. This manual takes Delft3d as an example for explanation, but in practice it can also be achieved through other deposition simulation software.

[0061] The process is as follows:

[0062] The data required for the algorithm to be input is output by the third-party software Delft3d.

[0063] The data output by Delft3d are terrain elevation data and sedimentary lithology data.

[0064] Assume that the simulation time length of Delft3d is 1 to n years, and assume that a time snapshot is output every year, and the output snapshot is named after the simulation era.

[0065] Then, the output of Delft3d is the sediment lithology time slices numbered 1 to n (that is, Figure 2 -abcde), and the slices contain high-precision elevation information.

[0066] Define a three-dimensional empty data volume (expressed in the programming code as applying for an empty three-dimensional data volume [(nx,ny,nk)—nx,ny represents a plane, nk is the number of time slices], which is used to store sedimentation simulation data for subsequent calculations), put the data into the three-dimensional data volume from old to new, put the data deposited first at the bottom, and put the data deposited later at the top, which conforms to the actual deposition time contact relationship; that is, put the data of the first year at the bottom of the three-dimensional data volume first, and then put the data of the second year, the third year... and finally the nth year). The data stored in the three-dimensional data volume is time domain data.

[0067] After the s3-s8 process, the data in the time domain three-dimensional data volume is used as input to obtain a new space domain three-dimensional data volume, which is the actual geological situation underground. (nx, ny, nz) represents the nx horizontal, ny vertical, and nz depth coordinates in space. Figure 3 (as shown in abcd).

[0068] Delft3d is a hydrodynamic simulation software developed by the Delft Institute in the Netherlands. Delft3d can be used to study hydrodynamics and sediment problems in various water bodies such as rivers, lakes, deep-water reservoirs and near-shore oceans.

[0069] Although Delft3d is a comprehensive software with functions such as hydrodynamic calculation, wave calculation, sediment transport and deposition, riverbed and seabed topography prediction, etc., it is not designed for geological modeling or geophysical modeling. Secondly, Delft3d hydrodynamic simulation software uses the finite difference method based on solving the Navier-Stokes equation (hereinafter referred to as the NS equation) to calculate the simulation results, so the calculation results are a collection of two-dimensional slice information at the time sampling interval. The simulation results are shown as follows: Figure 2 As shown, it only shows the two-dimensional value of the model at a certain moment. Figure 2 The figure shows the change in underwater terrain elevation. Similarly, other two-dimensional information such as lithology, water flow velocity, and water level can be obtained. Figure 2 a is the result of simulating delta deposition after 60,000 time steps. Figure 2 be are the results after calculating 120,000, 180,000, and 240,000 time steps respectively, and the calculation results are only in two-dimensional form. Figure 2 This is a display of Delft3d hydrodynamic simulation results, including: a. 60,000 Δt simulation results; b. 120,000 Δt iterative simulation results; c. 180,000 Δt simulation results; d. 240,000 Δt simulation results; e. 300,000 Δt simulation results, Δt represents the time step, or the interval between sampling points during simulation calculation.

[0070] The above functions of Ddelft3d make it possible to make geological models and geophysical models. If you want to realize the geological and geophysical modeling of this software, you need to calculate the two-dimensional data in the time series in accordance with the sedimentological characteristics and convert the time series slice data into a three-dimensional numerical model in the depth domain.

[0071] According to the principle of Delft3d deposition simulation software, the deposition simulation calculation is performed in the time domain, so the final calculation result can actually be described as a data body with a "book" structure. The data body can be described as (nx×ny)×nk, where nt represents the number of time samplings, that is, the total number of calculation steps; nx represents the number of sampling points in the x direction; ny represents the number of sampling points in the y direction, so the data obtained is a set of two-dimensional spatial data slices indexed by the time domain.

[0072] This two-dimensional spatial data slice set records the model data parameters at nk moments within the sedimentation simulation operation time, but does not have the characteristics of a spatial three-dimensional data body. It is necessary to calculate the sedimentation simulation process from the past to the present to obtain a geological model or geophysical model with hydrodynamic characteristics.

[0073] Here we hope to obtain a geological model in the depth domain, so as to study the subsequent delta geological model forward modeling, seismic phases, and attributes. It is necessary to convert the time domain data into depth domain data, that is, to convert (nx×ny)×nk into nx×ny×nz, where nt has no direct correlation with the value of nz. The value of nz is specified by the user. Generally, the value of nz×dz is greater than the maximum vertical thickness of the research model, where nz represents the number of sampling points in the depth direction converted to the geological or geophysical model. Here, it is assumed that the velocity is high where the sediment particles are large in the geological model, and the velocity is low where the sediment particles are small (consistent with the traditional wave impedance, such as high wave impedance of sandstone and low wave impedance of mudstone).

[0074] When the simulation is completed, the erosion or deposition is calculated from the old to the new in time according to the calculated intermediate delta morphological elevation data and the grain size distribution data. Then, the corresponding data processing is performed to obtain the three-dimensional depth domain lithology model. After converting the lithology model according to the corresponding relationship, the three-dimensional depth domain velocity model can be obtained. The logic flow chart is as follows: Figure 1 As shown, according to the above method, Figure 2 After the simulation experiment in the above, the three-dimensional lithology model is obtained. Figure 3 shown.

[0075] The lithology modeling logic takes into account the erosion of the river channel. The lithology is calculated by interpolation from the bottom to the top moment by moment. If the lithology velocity body is needed, it is only necessary to add a corresponding relationship between the lithology and the velocity to obtain the velocity body, which can be used for geophysical forward modeling research.

[0076] Figure 3 The delta shown in (a. delta sedimentation simulation results; b. river channel erosion; c. lithology profile along the provenance direction; d. lithology profile perpendicular to the provenance direction; e. lithology root mean square amplitude; f. lithology average value) is a multi-stage delta simulated and calculated. During the simulation period, the water level gradually rises. The lithology profiles AA' and BB' at the cross section are shown in Figure 1. Figure 3 As shown in c and d, Figure 3 c is the section along the source direction, Figure 3 d is the cross section perpendicular to the provenance direction. The A end is the provenance direction, the A' end is the basin direction, Figure 3 d is the section perpendicular to the provenance direction, on which lithological features such as estuary bar sandstone, channel sand, and delta plain mudstone can be seen, proving that the simulation results obtained in this experimental simulation are consistent with the principles of sedimentology.

[0077] Following the method of geophysical stratigraphic interpretation, the top and bottom interfaces of the delta are traced and interpreted, and then the root mean square amplitude of the sediment lithology is calculated (e.g. Figure 3e) The red area indicates that the lithology is highly differentiated and has undergone relatively large changes, and is very sensitive to thin sand bodies; the blue area indicates that the lithology is single and has not changed. Figure 3 f represents the direct accumulation and statistics of lithology ( Figure 3 e and f are non-seismic attribute extractions, where only the root mean square value and statistical lithology value are directly calculated for the lithology. Red indicates areas with thick sandstone, and green and blue indicate areas with low sandstone content. Figure 3 e and f can be used to compare and study the relationship between lithology and seismic attributes. Since the sedimentary simulation results can be displayed in three dimensions after processing, the delta lithology can be observed in sections and slices like seismic data. Figure 3 It is proved that this method can well convert the Delft3d hydrodynamic simulation results into a geological model after the sedimentation simulation process processing algorithm calculation, which is very realistic and has hydrodynamic characteristics. The geophysical model obtained by using the geophysical method to interpret the geological model is very consistent with the actual geological sedimentary characteristics, has very rich lithological variation characteristics, and the river channel is fresh and visible. This geological and geophysical model (depth domain geological model) generation method can provide new methods and ideas for oil and gas exploration, and has the prospect of practical application in oil and gas exploration.

[0078] The above-described embodiments are used to illustrate the present invention, not to limit the present invention, so changes in the exemplified values ​​or replacement of equivalent elements should still fall within the scope of the present invention.

[0079] From the above detailed description, it can be understood by those skilled in the art that the present invention can indeed achieve the aforementioned objectives and is in compliance with the provisions of the Patent Law.

[0080] Although the preferred embodiments of the present invention have been described, those skilled in the art may make additional changes and modifications to these embodiments once they are aware of the basic creative concepts. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications that fall within the scope of the present invention. The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. It should be pointed out that any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

[0081] It should be noted that the above description of the relevant process is only for example and explanation, and does not limit the scope of application of this specification. For those skilled in the art, various modifications and changes can be made to the process under the guidance of this specification. However, these modifications and changes are still within the scope of this specification.

[0082] The basic concepts have been described above. Obviously, for those of ordinary skill in the art who have read this application, the above invention disclosure is only for example and does not constitute a limitation of this application. Although not explicitly stated here, those of ordinary skill in the art may make various modifications, improvements and amendments to this application. Such modifications, improvements and amendments are suggested in this application, so such modifications, improvements and amendments still belong to the spirit and scope of the exemplary embodiments of this application.

[0083] At the same time, the present application uses specific words to describe the embodiments of the present application. For example, "one embodiment", "an embodiment", and / or "some embodiments" refer to a certain feature, structure or characteristic related to at least one embodiment of the present application. Therefore, it should be emphasized and noted that "one embodiment" or "an embodiment" or "an alternative embodiment" mentioned twice or more in different positions in this specification does not necessarily refer to the same embodiment. In addition, some features, structures or characteristics in one or more embodiments of the present application can be appropriately combined.

[0084] In addition, it will be appreciated by those skilled in the art that various aspects of the present application may be illustrated and described by a number of patentable categories or situations, including any new and useful combination of processes, machines, products or substances, or any new and useful improvements thereto. Therefore, various aspects of the present application may be implemented entirely by hardware, entirely by software (including firmware, resident software, microcode, etc.), or by a combination of hardware and software. The above hardware or software may all be referred to as "units", "modules" or "systems". In addition, various aspects of the present application may take the form of a computer program product embodied in one or more computer-readable media, wherein computer-readable program code is contained therein.

[0085] The computer program code required for the operation of each part of the application can be written in any one or more programming languages, including object-oriented programming languages ​​such as Java, Scala, Smalltalk, Eiffel, JADE, Emerald, C++, C#, VB.NET, Python, etc., such as C programming language, VisualBasic, Fortran2103, Perl, COBOL2102, PHP, ABAP, such as Python, Ruby and Groovy dynamic programming language or other programming languages. The program code can be run completely on the user's computer, or run on the user's computer as an independent software package, or run partly on the user's computer and partly on the remote computer, or run completely on the remote computer or server. In the latter case, the remote computer can be connected to the user's computer through any network form, such as a local area network (LAN) or a wide area network (WAN), or connected to an external computer (e.g., via the Internet), or in a cloud computing environment, or used as a service such as software as a service (SaaS).

[0086] In addition, unless explicitly stated in the claims, the order of the processing elements and sequences described in this application, the use of alphanumeric characters, or the use of other names are not intended to limit the order of the processes and methods of this application. Although the above disclosure discusses some invention embodiments that are currently considered useful through various examples, it should be understood that such details are only for illustrative purposes, and the attached claims are not limited to the disclosed embodiments. On the contrary, the claims are intended to cover all modifications and equivalent combinations that are consistent with the essence and scope of the embodiments of this application. For example, although the implementation of the various components described above can be embodied in a hardware device, it can also be implemented as a pure software solution, for example, installation on an existing server or mobile device.

[0087] Similarly, it should be noted that in order to simplify the description disclosed in this application and thus help understand one or more embodiments of the invention, in the above description of the embodiments of the application, multiple features are sometimes combined into one embodiment, figure or description thereof. However, this method of the application should not be interpreted as reflecting the intention that the claimed object requires more features than those explicitly stated in each claim. On the contrary, the subject of the invention should have fewer features than the above single embodiment.

Claims

1. A method for converting sedimentary simulation results into a three-dimensional geological model, characterized in that: include: S1. Obtain a known time domain model (nx, ny, nk), where nx represents the number of sampling points in the x direction; ny represents the number of sampling points in the y direction, and nk represents the number of time sampling points; S2. Create a three-dimensional matrix (nx, ny, nz) and set the initial value of the 3D matrix to null to save the depth domain geological model; nz represents the number of sampling points in the depth direction converted to the geological model; S3. Set a variable ik as the time pointer, and set variables ix and iy as the x-direction pointer and y-direction pointer respectively; the initial value of ik is nk, and the value range of ik is [1,nk]; the initial values ​​of ix and iy are set to 1, and the value ranges are [1,nx] and [1,ny] respectively; S4. Compare the delta terrain elevation data at the current moment (ix, iy, ik) with the delta terrain elevation data at the previous moment (ix, iy, ik-1) to determine whether there is deposition at (ix, iy); S5. If there is deposition at (ix, iy), then deposition processing is performed: the three-dimensional matrices (ix, iy, ik-1) and (ix, iy, ik) are assigned sediment grain size data; if there is no deposition at (ix, iy), then erosion processing is performed: the values ​​at the three-dimensional matrices (ix, iy, ik-1) and (ix, iy, ik) are modified to null to achieve an erosion effect; S6. Determine whether ix has reached the maximum value. If so, jump to S7; if not, assign ix=ix+1 and jump to S4; S7. Determine whether iy reaches the maximum value. If so, jump to S8; if not, assign iy=iy+1 and jump to S4; S8. Determine whether the time steps have been traversed. If so, save the result to obtain the depth domain geological model; if not, assign: ik=ik-1, and jump to S4.

2. A system for converting sedimentary simulation results into a three-dimensional geological model, characterized in that: include: The time domain model acquisition module is used to obtain the known time domain model (nx, ny, nk), where nx represents the number of sampling points in the x direction; ny represents the number of sampling points in the y direction, and nk represents the number of time sampling points; The three-dimensional matrix creation module is used to create a three-dimensional matrix (nx, ny, nz) and set the initial value of the 3D matrix to null to save the depth domain geological model; nz represents the number of sampling points in the depth direction converted to the geological model; The assignment module is used to set a variable ik as a time pointer, and set variables ix and iy as x-direction pointer and y-direction pointer respectively; the initial value of ik is nk, and the value range of ik is [1,nk]; the initial values ​​of ix and iy are set to 1, and the value ranges are [1,nx] and [1,ny] respectively; The comparison module is used to compare the delta terrain elevation data at the current moment (ix, iy, ik) with the delta terrain elevation data at the previous moment (ix, iy, ik-1) to determine whether there is deposition at (ix, iy); The deposition or erosion module is used to perform deposition if deposition occurs at (ix, iy): the three-dimensional matrices (ix, iy, ik-1) and (ix, iy, ik) are assigned sediment grain size data; If there is no deposition at (ix, iy), then erosion is performed: the values ​​at (ix, iy, ik-1) and (ix, iy, ik) of the three-dimensional matrix are modified to null to achieve the erosion effect; The first judgment module is used to judge whether ix reaches the maximum value. If so, it jumps to the second judgment module; if not, it assigns: ix=ix+1 and jumps to the comparison module; The second judgment module is used to judge whether iy reaches the maximum value. If so, it jumps to the third judgment module; if not, it assigns: iy=iy+1 and jumps to the comparison module; The third judgment module is used to judge whether the time steps have been traversed. If so, the result is saved to obtain the depth domain geological model; if not, the value is assigned: ik=ik-1, and jump to the comparison module.

3. A computer-readable storage medium, characterized in that: The storage medium stores computer instructions. When a computer reads the computer instructions, the computer executes the method for converting sedimentation simulation results into a three-dimensional geological model as claimed in claim 1.

Images