DEM and GIS technology-based ancient basin deposition center reconstruction method

Abstract

The invention discloses an ancient basin deposition center reconstruction method based on DEM and GIS technologies, and belongs to the technical field of marine science and marine geology and hydrology, and the method comprises the steps: firstly obtaining data of a to-be-researched region, and constructing a digital elevation model based on the data; performing water system analysis based on the digital elevation model to obtain a paleowater system overall distribution map, a basin area and a lithology distribution map; reconstructing the overall palaeowater system distribution diagram to obtain a palaeowater system reconstruction diagram; generating source-sink data based on a source-sink relationship in the paleowater system reconstruction graph; and obtaining a relational expression between the river channel and the fan body area based on the source sink data, and obtaining the scale of the deposition center based on the relational expression between the river channel and the fan body area. The invention provides an efficient and accurate paleo-water system reconstruction method, through combination of DEM and GIS technologies, a reconstruction graph of a paleo-water system can be rapidly generated, and correlation analysis and estimation are carried out. The method has a wide application prospect in the fields of geology and hydrology.

Description

Technical Field

[0001] The present invention belongs to the technical fields of marine science, marine geology and hydrology, and particularly relates to a method for reconstructing the sedimentation center of an ancient basin based on DEM and GIS technologies. Background Art

[0002] The reconstruction of the ancient water system of the sedimentation center of an ancient basin is of great significance for understanding geological history, sedimentation processes and environmental evolution. Traditional methods for reconstructing the ancient water system of the sedimentation center of an ancient basin rely on geological surveys and manual drawing, with low efficiency and limited accuracy. With the development of digital elevation model (DEM) and geographic information system (GIS) technologies, it has become possible to reconstruct and restore the ancient water system of the sedimentation center of an ancient basin using these technologies. However, there is no solution to address the problems existing in the above-mentioned prior art. Summary of the Invention

[0003] To solve the above technical problems, the present invention provides a method and system for reconstructing the sedimentation center of an ancient basin based on DEM and GIS technologies to address the problems existing in the above-mentioned prior art.

[0004] In a first aspect, to achieve the above object, the present invention provides a method for reconstructing the sedimentation center of an ancient basin based on DEM and GIS technologies, including the following steps:

[0005] Obtain data of the area to be studied, and construct a digital elevation model based on the data;

[0006] Perform water system analysis based on the digital elevation model to obtain an overall distribution map of the ancient water system, the drainage area, and a lithology distribution map;

[0007] Reconstruct the overall distribution map of the ancient water system to obtain a reconstructed map of the ancient water system;

[0008] Generate source-sink data based on the source-sink relationship in the reconstructed map of the ancient water system;

[0009] Obtain a relational expression between the area of the river channel and the fan body based on the source-sink data, and obtain the scale of the sedimentation center based on the relational expression between the area of the river channel and the fan body.

[0010] Optionally, the data of the area to be studied includes: artificial seismic data volume, seismic data interpretation results, and geological data, where the geological data includes well logging curves, well logging interpretation results, and formation velocity data.

[0011] Optionally, the process of performing water system analysis based on the digital elevation model to obtain an overall distribution map of the ancient water system, the drainage area, and a lithology distribution map includes:

[0012] Generate an overall distribution map of the ancient water system of the area to be studied based on the digital elevation model;

[0013] Calculate the drainage area of each basin in the overall distribution map of the ancient water system;

[0014] Generate a lithology distribution map of the source area based on the geological data.

[0015] Optionally, the process of generating the overall distribution map of the ancient water system of the area to be studied based on the digital elevation model includes:

[0016] Perform depression filling analysis, flow direction analysis, flow analysis, map algebra analysis, river link analysis, river network grading analysis, raster river network vectorization analysis, basin domain analysis, and catchment area analysis based on the data to obtain the overall distribution map of the ancient water system. Among them, basin domain analysis is performed to obtain a basin domain analysis map.

[0017] Optionally, the process of calculating the drainage area of each basin in the overall distribution map of the ancient water system includes:

[0018] The area of each basin includes: the source area, the channel, and the confluence area. Among them, the source area is the part of the head and tributaries, the channel is the main stream part of the river, and the confluence area is the part where the alluvial fan formed in front of the mountain by the river channel converges, and this part includes sedimentary bodies.

[0019] Optionally, the process of generating a lithology distribution map of the source area based on the geological data includes:

[0020] Generate a seismic attribute map based on the artificial seismic data volume, and obtain the corresponding results based on the correlation analysis between the attribute distribution in the seismic attribute map and the lithology interpretation on the well log in the well log interpretation results;

[0021] Based on the corresponding results, obtain the lithology distribution map through the attribute map.

[0022] Optionally, the process of generating source-sink data based on the source-sink relationship in the ancient water system reconstruction map includes:

[0023] Obtain a gully shape statistical table based on the ancient water system reconstruction map and the basin domain analysis map;

[0024] Obtain source-sink data based on the gully shape statistical table;

[0025] Among them, the source-sink data includes: the source area, the length of the channel part, the river transportation distance, the river curvature, and the fan area of the convergence area.

[0026] In a second aspect, the present invention also provides an ancient basin sediment center reconstruction system based on DEM and GIS technologies, for implementing an ancient basin sediment center reconstruction method based on DEM and GIS technologies. The system includes:

[0027] A data acquisition module, configured to acquire data of the area to be studied and construct a digital elevation model based on the data;

[0028] A water system analysis module, configured to perform water system analysis based on the digital elevation model to obtain an overall distribution map of the ancient water system, a drainage area, and a lithology distribution map;

[0029] An ancient water system reconstruction module, configured to reconstruct the overall distribution map of the ancient water system to obtain a reconstructed map of the ancient water system;

[0030] A source-sink relationship generation module, configured to generate source-sink data based on the source-sink relationship in the reconstructed map of the ancient water system;

[0031] A sedimentation center scale calculation module, configured to generate a relational expression between the river channel and the fan area based on the source-sink data, and calculate the scale of the sedimentation center based on the relational expression.

[0032] In a third aspect, the present invention further provides a computer terminal device, including:

[0033] One or more processors;

[0034] A memory, coupled to the processor, configured to store one or more programs;

[0035] When the one or more programs are executed by the one or more processors, the one or more processors implement the steps of an ancient basin sedimentation center reconstruction method based on DEM and GIS technologies.

[0036] In a fourth aspect, the present invention further provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of an ancient basin sedimentation center reconstruction method based on DEM and GIS technologies are implemented.

[0037] Compared with the prior art, the present invention has the following advantages and technical effects:

[0038] An ancient basin sedimentation center reconstruction method and system based on DEM and GIS technologies provided by the present invention first obtain data of a region to be studied, construct a digital elevation model based on the data; perform water system analysis based on the digital elevation model to obtain an overall distribution map of the ancient water system, a drainage area, and a lithology distribution map; reconstruct the overall distribution map of the ancient water system to obtain a reconstructed map of the ancient water system; generate source-sink data based on the source-sink relationship in the reconstructed map of the ancient water system; obtain a relational expression between the river channel and the fan area based on the source-sink data, and obtain the scale of the sedimentation center based on the relational expression between the river channel and the fan area. The present invention provides an efficient and accurate ancient water system reconstruction method. By combining DEM and GIS technologies, a reconstructed map of the ancient water system can be quickly generated, and relevant analysis and estimation can be performed. This method has broad application prospects in the fields of geology and hydrology. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] The accompanying drawings, which form a part of this application, are used to provide a further understanding of this application. The schematic embodiments and descriptions thereof of this application are used to explain this application and do not constitute an improper limitation of this application. In the drawings:

[0040] Figure 1 It is an example diagram of an embodiment of the present invention;

[0041] Figure 2 It is a schematic diagram of the method of an embodiment of the present invention. Detailed implementation manners

[0042] It should be noted that, without conflict, the embodiments in this application and the features in the embodiments may be combined with each other. The following will detail this application with reference to the accompanying drawings and in combination with the embodiments.

[0043] It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although the logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in a different order from that here.

[0044] Embodiment 1

[0045] As Figure 1-2 shown, in this embodiment, a method for reconstructing the sedimentation center of an ancient basin based on DEM and GIS technologies is provided, including:

[0046] Obtain the data of the area to be studied and construct a digital elevation model based on the data;

[0047] Conduct a river system analysis based on the digital elevation model to obtain the overall distribution map of the ancient river system, the drainage area, and the lithology distribution map;

[0048] Reconstruct the overall distribution map of the ancient river system to obtain the reconstructed map of the ancient river system;

[0049] Generate source-sink data based on the source-sink relationship in the reconstructed map of the ancient river system;

[0050] Obtain the relationship formula between the river channel and the fan area based on the source-sink data, and obtain the scale of the sedimentation center based on the relationship formula between the river channel and the fan area.

[0051] Specifically, it includes the following steps:

[0052] 1. Collect the geological data and relevant literature of the study area.

[0053] 2. Export the DEM using Petrel software.

[0054] 3. Use Arcgis software to conduct a river system analysis and generate a river system map, a drainage area, and a source area lithology distribution map.

[0055] 4. Generate a planar distribution map of the water system based on the DEM and reconstruct the ancient water system.

[0056] 5. Study the source-sink relationship of the ancient water system and generate a statistical chart of the gully shape.

[0057] 6. Analyze the characteristics of the river channels and fans and generate a statistical table.

[0058] 7. Establish a relationship between the areas of the river channels and fans, analyze the results in combination with the literature, and estimate the scale of the deposition center.

[0059] As an implementation manner in this embodiment, the data of the area to be studied includes: artificial seismic data volume, seismic data interpretation results, and geological data, where the geological data includes logging curves, logging interpretation results, and formation velocity data.

[0060] As an implementation manner in this embodiment, the process of obtaining the overall distribution map of the ancient water system, the drainage area, and the lithology distribution map based on the digital elevation model includes:

[0061] Generate the overall distribution map of the ancient water system of the area to be studied based on the digital elevation model;

[0062] Calculate the drainage area of each basin in the overall distribution map of the ancient water system;

[0063] Generate the lithology distribution map of the source area based on the geological data.

[0064] Specifically, it includes the following steps:

[0065] 1. Data collection and preparation

[0066] 1.1 Search for literature

[0067] Input: Research topic

[0068] Output: Relevant literature

[0069] Obtain information such as the geological background and historical water system distribution of the research area by searching for relevant literature.

[0070] 1.2 Export DEM from Petrel

[0071] Input: Artificial seismic data volume, seismic data interpretation results, geological data (including logging curves, logging interpretation results, formation velocity data, etc.)

[0072] Output: Digital elevation model (DEM)

[0073] Export the formation data containing depth (i.e., elevation) information of the study area using Petrel software, import it into ArcGIS software, and through creating TIN (abbreviation for gridded data in the software) - converting TIN to raster - clipping raster (clip the study area you need) - adjusting the symbol system (the clipped DEM may need to adjust the symbol system to make it beautiful, and the symbol system can be adjusted through the layer properties settings) - outputting contour lines - outputting the final required digital elevation model (DEM).

[0074] As an implementation manner in this embodiment, the process of generating the overall distribution map of the ancient water system in the area to be studied based on the digital elevation model includes:

[0075] Conduct depression filling analysis, flow direction analysis, flow volume analysis, map algebra analysis, river linking analysis, river network grading analysis, raster river network vectorization analysis, basin domain analysis, and catchment area analysis based on the data to obtain the overall distribution map of the ancient water system. Among them, basin domain analysis is conducted to obtain the basin domain analysis map.

[0076] Specifically, it includes the following steps:

[0077] 2. Ancient water system analysis

[0078] 2.1 Generating the planar distribution map of the ancient water system using Arcgis

[0079] Input: DEM

[0080] Output: Overall distribution map of the ancient water system

[0081] Using Arcgis software, generate the overall distribution map of the ancient water system in the study area based on DEM. The specific generation process is as follows: depression filling analysis - flow direction analysis - flow volume analysis - map algebra analysis - raster calculator (extracting hydrological data with a certain flow volume) - river linking analysis - river network grading analysis - raster river network vectorization analysis (the river can be smoothed) - basin domain analysis (inputting flow direction) - catchment area analysis (inputting flow direction + river linking) - obtaining the overall distribution map of the ancient water system. The above operation process can be realized in ArcGIS software.

[0082] As an implementation manner in this embodiment, in the process of calculating the basin area of each basin in the overall distribution map of the ancient water system, it includes:

[0083] The areas of each basin include: source area, channel, and confluence area. Among them, the source area is the part of the headwaters and tributaries, the channel is the part of the main stream of the river, and the confluence area is the part where the alluvial fan formed in front of the mountain by the river channel converges, and this part includes sedimentary bodies.

[0084] Specifically, it includes the following steps:

[0085] 2.2 Generating the basin area using Arcgis

[0086] Input: Ancient river system map

[0087] Output: Watershed area

[0088] Based on the ancient river system map, the area of each watershed is calculated using Arcgis. The specific calculation process is as follows: In this article, the specific definitions of the source area, channels, and confluence area are first set in the research context, as shown in the figure. The green shaded area S1 in the figure represents the mountain source, including the tributaries of the source, reflecting the coverage area of the source part; the bold blue solid line L represents the channel, that is, the main stream part of the river, and its length is L; while L’ represents the straight-line distance of the channel main stream on the plane, which can be used to calculate the bending degree of the main stream; finally, the yellow shaded area S2 represents the confluence area, indicating the sediment body of the alluvial fan body formed by the river before the mountain.

[0089] As an implementation method in this embodiment, the process of generating a lithology distribution map of the source area based on the geological data includes:

[0090] Generate a seismic attribute map based on the artificial seismic data volume, and obtain the corresponding results through the correlation analysis of the attribute distribution in the seismic attribute map and the lithology interpretation on the well log in the well log interpretation results;

[0091] Based on the corresponding results, obtain the lithology distribution map through the attribute map.

[0092] Specifically, it includes the following steps:

[0093] 2.3 Generate a lithology distribution map of the source area based on the new seismic attributes (see appendix): Note that the attributes mentioned in the appendix are used here, which are the attributes derived by myself

[0094] Input: Geological data and seismic data

[0095] Output: Lithology distribution map of the source area

[0096] Using Python software, generate a lithology distribution map of the source area according to the geological data. Different values in the seismic attribute map correspond to different lithologies. After obtaining the attribute distribution map, the corresponding results can be obtained by using the correlation analysis between the attributes and the lithology interpretation on the well log; using the corresponding results, the lithology map can be obtained through the attribute map.

[0097] The derivation process of the new seismic attributes includes:

[0098] The RMS attribute, that is, the Root Mean Square, is an important indicator for measuring the signal strength or energy. In seismology, the RMS value is often used to compare the intensity differences of different seismic events or the energy sizes of different seismic data. In addition, RMS can also be used for processing such as noise filtering and signal enhancement of seismic signals.

[0099] For a time - series signal x(t), the RMS can be calculated by the following formula:

[0100] RMS = sqrt((1 / N)*sum(x(t)^2))

[0101] Where x(t) represents the amplitude value of the time - series signal at time t, and N represents the number of sampling points in the time - series signal. The calculation process of this formula is as follows: First, square the amplitude value of each sampling point of the signal, then sum all the squared values and divide by the number of sampling points N, and finally perform a square - root operation on the result to obtain the RMS property. The larger the RMS value, the greater the strength or energy of the signal.

[0102] Envelope usually refers to the envelope curve of the seismic signal. It is the upper limit of the amplitude of the seismic signal and describes the contour of the seismic signal changing with time. There are various methods for extracting the envelope curve, and the most commonly used method is to extract it by performing a Hilbert transform on the seismic signal.

[0103] Assume the input signal is x(t), where t represents time. Then, perform a Hilbert transform on x(t) to obtain a complex signal:

[0104] z(t)=x(t)+j*y(t)

[0105] Where j represents the imaginary unit, y(t) represents the Hilbert transform of x(t), and z(t) is the analytic complex signal, which contains the amplitude and phase information of the original signal x(t).

[0106] Subsequently, the magnitude of the complex signal can be calculated:

[0107] A(t)=sqrt(x(t)^2 + y(t)^2)

[0108] A(t) is the envelope curve of the seismic signal.

[0109] Therefore, the formula for extracting the Envelope property of the seismic signal is:

[0110] Envelope(t)=sqrt(x(t)^2 + y(t)^2)

[0111] Where x(t) represents the original seismic signal, and y(t) represents the Hilbert transform of x(t).

[0112] Taking the Envelope property body as the input signal and extracting the RMS property of this signal can generate a new property EMS (i.e., Env_RMS appearing in the text). The formula for the new seismic property is:

[0113] EMS = sqrt((1 / N) * sum(envelope^2))

[0114] Where Envelope represents the envelope attribute body of the input signal, and N represents the number of samples of the input signal.

[0115] When substituting the Envelope formula into the EMS formula, the following combined formula can be obtained:

[0116] EMS = sqrt((1 / N) * sum((sqrt(x(t)^2 + y(t)^2))^2))

[0117] Further simplification gives:

[0118] EMS = sqrt((1 / N) * sum(x(t)^2 + y(t)^2))

[0119] It can also be written as:

[0120] EMS = sqrt((1 / N) * (sum(x(t)^2) + sum(y(t)^2)))

[0121] Where N represents the number of samples of the input signal, x(t) is the seismic signal, and y(t) represents the Hilbert transform of x(t).

[0122] The formula EMS = sqrt((1 / N) * sum((sqrt(x(t)^2 + y(t)^2))^2)) can be used as a new seismic attribute to measure the intensity or energy of the seismic signal.

[0123] This formula is based on the extraction of the envelope of the seismic signal to calculate the EMS attribute of the signal. First, sum the squared amplitudes of the envelope and divide by the number of sampling points, and then take the square root. For different data and application scenarios, it is necessary to select a suitable envelope extraction method and adjust and optimize the formula to obtain more accurate seismic attribute values. Using the root mean square value of the envelope amplitude is more convenient in practice. The RMS value can be calculated within a set window, and the resolution of the result can be optimized by adjusting the window length.

[0124] Python implementation program for EMS attribute:

[0125] In Python, the NumPy library can be used to implement the calculation of the EMS attribute. The following is an example code:

[0126] import numpy as np

[0127] def calculate_EMS(signal):

[0128] # Perform Hilbert transform on the signal

[0129] hilbert_transform = np.imag(np.hilbert(signal))

[0130] # Calculate EMS attribute

[0131] envelope = np.sqrt(signal**2 + hilbert_transform**2)

[0132] EMS = np.sqrt(np.mean(envelope**2))

[0133] return EMS

[0134] # Example usage

[0135] # Assume there is a seismic signal stored in a NumPy array named "signal"

[0136] signal = np.array([0.1, 0.2, 0.3, 0.4, 0.3, 0.2, 0.1])

[0137] EMS_value = calculate_EMS(signal)

[0138] print("EMS attribute value:", EMS_value)

[0139] # In this example, first, the np.hilbert() function is used to perform the Hilbert transform on the input signal to obtain the Hilbert transform signal. Next, the envelope is calculated according to the EMS formula, and the mean of the sum of squares is calculated. Finally, the square root is taken to obtain the EMS attribute value.

[0140] As an implementation manner in this embodiment, reconstruct the overall distribution map of the ancient water system to obtain a reconstructed map of the ancient water system;

[0141] Specifically, it includes the following steps:

[0142] 3. Reconstruction of the ancient water system

[0143] Input: The planar distribution map of the ancient water system generated in step 2.1

[0144] Output: Reconstructed map of the ancient water system

[0145] According to the planar distribution map of the water system, combined with the geological background and historical documents, reconstruct the distribution map of the ancient water system. This planar map not only includes the distribution of rivers but also the information of lithology distribution. The specific formation of the map needs to be obtained manually with the help of drawing software. This reconstructed map is a more information-rich ancient water system reconstruction map than the information in step 2.1.

[0146] As an implementation method in this embodiment, the process of generating source-sink data based on the source-sink relationship in the ancient water system reconstruction map includes:

[0147] Obtain the gully shape statistical table based on the ancient water system reconstruction map and the basin domain analysis map;

[0148] Obtain the source-sink data based on the gully shape statistical table;

[0149] Among them, the source-sink data includes: source area, length of the channel part, river transportation distance, river curvature, and fan area of the convergence area.

[0150] Specifically, it includes the following steps:

[0151] 4. Research on the source-sink of the ancient water system

[0152] 4.1 Generate the gully shape statistical table of the geological section by Arcgis and expert experience

[0153] Input: DEM map, ancient water system reconstruction map, and basin domain analysis map obtained in 2.1

[0154] Output: Gully shape statistical table

[0155] Using Arcgis software and drawing software, with the help of expert experience, divide the positions of the mountain passes and the shapes, depths, and widths of the mountain pass gullies within the scope of the ancient water system basin according to the DEM map, ancient water system reconstruction map, and basin domain analysis map obtained in 2.1 to obtain the gully shape statistical table.

[0156] 4.2 Research on the source-sink of the ancient water system

[0157] Input: Gully shape statistical table

[0158] Output: Source-sink parameters

[0159] Based on the gully shape statistical table, and Figure 1 , obtain the source-channel-sink parameters, and the parameters include source area, length of the channel part, river transportation distance, river curvature, and fan area of the convergence area.

[0160] As an implementation method in this embodiment, obtain the relationship formula between the river channel and the fan area based on the source-sink data, and obtain the scale of the deposition center based on the relationship formula between the river channel and the fan area.

[0161] Specifically, it includes the following steps:

[0162] 5. Result Analysis and Estimation

[0163] 5.1 Relationship between River Channel and Fan Area

[0164] Input: The above source-channel-confluence parameters

[0165] Output: Relationship

[0166] Based on the statistical data of source-channel-confluence, establish the fitting relationship between the source area (S1), the length of the channel part (L), the river transportation distance (H), the river curvature (W) and the convergent fan area (S2).

[0167] The relationship is a fitting polynomial, specifically as follows:

[0168] S2 = a * S1 + b * L + c * H + d * W

[0169] Among them, a, b, c, d are constants, which can be obtained by fitting the statistical data.

[0170] 5.2 Estimate the Scale of the Deposition Center

[0171] Input: Fitting relationship

[0172] Output: Scale of the deposition center

[0173] Based on the results of paleo-water system reconstruction, lithology distribution results and river channel fan analysis results, determine the location of the deposition center and the total area S of the deposition center formed by multiple convergent fans.

[0174] Among them, the deposition center is formed by multiple convergent fans. The development position and area of the fan have been obtained in step 5.

[0175] S = ΣS2 i (i = 1, 2, 3...)

[0176] Based on this, a method for reconstructing the deposition center of a paleo-basin based on DEM and GIS technologies provided by an embodiment of the present invention first obtains the data of the area to be studied, constructs a digital elevation model based on the data; conducts water system analysis based on the digital elevation model to obtain the overall distribution map of the paleo-water system, the basin area and the lithology distribution map; reconstructs the overall distribution map of the paleo-water system to obtain the reconstructed map of the paleo-water system; generates source-sink data based on the source-sink relationship in the reconstructed map of the paleo-water system; obtains the relationship between the river channel and the fan area based on the source-sink data, and obtains the scale of the deposition center based on the relationship between the river channel and the fan area. The present invention provides an efficient and accurate method for reconstructing the paleo-water system. By combining DEM and GIS technologies, it can quickly generate the reconstructed map of the paleo-water system and conduct relevant analysis and estimation. This method has broad application prospects in the fields of geology and hydrology.

[0177] Example Two

[0178] In this embodiment, a computer terminal device is provided, including:

[0179] One or more processors;

[0180] A memory, coupled to the processor, for storing one or more programs;

[0181] When the one or more programs are executed by the one or more processors, the one or more processors implement the methods in the above embodiments.

[0182] In this embodiment, a computer-readable storage medium is further provided, on which a computer program is stored. When the computer program is executed by a processor, the methods in the above embodiments are implemented.

[0183] In this embodiment, an electronic device is further provided, including a memory and a processor. A computer program is stored in the memory, and the processor is configured to run the computer program to execute the methods in the above embodiments.

[0184] The above programs can run in a processor, or can also be stored in a memory (or referred to as a computer-readable medium). The computer-readable medium includes permanent and non-permanent, removable and non-removable media, and information storage can be implemented by any method or technology. The information can be computer-readable instructions, data structures, program modules, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassette tapes, magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information accessible by a computing device.

[0185] These computer programs can also be loaded onto a computer or other programmable data processing device, so that a series of operation steps are executed on the computer or other programmable device to generate computer-implemented processing. Thus, the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in a process Figure 1 One process or multiple processes and / or blocks Figure 1 One block or multiple blocks, and corresponding different steps can be implemented by different modules.

[0186] In this embodiment, such a device or system is provided. The system is called a paleo-basin sediment center reconstruction system based on DEM and GIS technologies, and includes:

[0187] A data acquisition module, configured to acquire data of the area to be studied and construct a digital elevation model based on the data;

[0188] A water system analysis module, configured to perform water system analysis based on the digital elevation model to obtain an overall distribution map of the paleo-water system, a drainage area, and a lithology distribution map;

[0189] A paleo-water system reconstruction module, configured to reconstruct the overall distribution map of the paleo-water system to obtain a reconstructed paleo-water system map;

[0190] A source-sink relationship generation module, configured to generate source-sink data based on the source-sink relationship in the reconstructed paleo-water system map;

[0191] A sediment center scale calculation module, configured to generate a relational expression of the river channel and fan body areas based on the source-sink data and calculate the scale of the sediment center based on the relational expression.

[0192] This system or device is used to implement the functions of the method in the above embodiment. Each module in this system or device corresponds to each step in the method, and those that have been described in the method will not be repeated here.

[0193] Through the above implementation manner, the problem of reconstructing the paleo-basin sediment center based on DEM and GIS technologies in the related art is solved, so that it can be ensured that the present invention provides an efficient and accurate paleo-water system reconstruction method. By combining DEM and GIS technologies, a reconstructed paleo-water system map can be quickly generated, and relevant analysis and estimation can be performed. This method has broad application prospects in the fields of geology and hydrology.

[0194] The above is only a preferred specific implementation manner of the present application, but the protection scope of the present application is not limited thereto. Any changes or substitutions that can be easily thought of by those skilled in the art within the technical scope disclosed in the present application should be covered by the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A method for reconstructing the sedimentary center of an ancient basin based on DEM and GIS technology, characterized in that: The following steps are involved: Obtain data of the area to be studied and build a digital elevation model based on the data; Based on the digital elevation model, water system analysis is performed to obtain the overall distribution map of the ancient water system, the drainage area and the lithology distribution map; Reconstructing the overall distribution map of the ancient water system to obtain a reconstruction map of the ancient water system; Generate source-sink data based on the source-sink relationship in the ancient river system reconstruction map; Based on the source-sink data, a relationship between the channel and the fan area is obtained, and based on the relationship between the channel and the fan area, the scale of the sedimentation center is obtained.

2. The method according to claim 1, characterized in that The data of the area to be studied include: artificial seismic data volume, seismic data interpretation results and geological data, wherein the geological data include well logging curves, well logging interpretation results and formation velocity data.

3. The method according to claim 2, characterized in that The process of analyzing the water system based on the digital elevation model to obtain the overall distribution map of the ancient water system, the drainage area and the lithology distribution map includes: Generate an overall distribution map of the ancient water system in the area to be studied based on the digital elevation model; Calculate the drainage area of ​​each drainage basin in the overall distribution map of the ancient drainage system; A lithology distribution map of the source area is generated based on the geological data.

4. The method according to claim 3, characterized in that The process of generating the overall distribution map of the ancient water system in the area to be studied based on the digital elevation model includes: Based on the data, we conduct depression filling analysis, flow direction analysis, flow analysis, map algebra analysis, river link analysis, river network classification analysis, grid river network vectorization analysis, basin analysis and watershed analysis to obtain the overall distribution map of the ancient water system. Among them, basin analysis is performed to obtain a basin analysis map.

5. The method according to claim 3, characterized in that: The process of calculating the drainage area of ​​each drainage basin in the overall distribution diagram of the ancient drainage system includes: The area of ​​each river basin includes: source area, channel and sink area, among which the source area is the source and tributary part, the channel is the main stream part of the river, and the sink area is the converging part of the alluvial fan formed by the river in front of the mountain, which includes sedimentary bodies.

6. The method according to claim 5, characterized in that The process of generating a lithology distribution map of a source area based on the geological data includes: Generate a seismic attribute map based on the artificial seismic data volume, and obtain corresponding results based on the correlation analysis between the attribute distribution in the seismic attribute map and the lithology interpretation on the well logging in the well logging interpretation result; Based on the corresponding results, the lithology distribution map is obtained through the attribute map.

7. The method according to claim 6, characterized in that Based on the source-sink relationship in the ancient river system reconstruction map, the process of generating source-sink data includes: Obtaining a valley shape statistics table based on the ancient water system reconstruction map and the basin analysis map; Obtain source-sink data based on valley shape statistics; The source-sink data include: source area, length of channel section, river transport distance, river curvature and convergence sector area.

8. A paleo-basin sediment center reconstruction system based on DEM and GIS technology, characterized in that: The system comprises: A data acquisition module is used to acquire data of the area to be studied and to construct a digital elevation model based on the data; A water system analysis module is used to perform water system analysis based on the digital elevation model to obtain an overall distribution map of the ancient water system, a drainage area, and a lithology distribution map; The ancient water system reconstruction module is used to reconstruct the overall distribution map of the ancient water system and obtain the ancient water system reconstruction map; Source-sink relationship generation module, used to generate source-sink data based on the source-sink relationship in the ancient river system reconstruction map; The sediment center scale calculation module is used to generate a relationship between the river channel and the fan area based on the source-sink data, and calculate the scale of the sediment center based on the relationship.

9. A computer terminal device, characterized in that: include: one or more processors; A memory, coupled to the processor, for storing one or more programs; When the one or more programs are executed by the one or more processors, the one or more processors implement the various steps of the ancient basin sedimentation center reconstruction method based on DEM and GIS technology as described in any one of claims 1-7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, each step of the method for reconstructing the sediment center of an ancient basin based on DEM and GIS technology as described in any one of claims 1 to 7 is implemented.

Images