Paleo-tidal range recovery method, device, equipment and medium based on medium tidal plateau thickness

Abstract

The invention relates to a paleo-tidal range recovery method, device, equipment and medium based on the thickness of a medium tidal plateau, in particular to the field of paleo-tidal range recovery in a sedimentary environment, and the paleo-tidal range recovery method comprises the following steps: obtaining the deposition of the medium tidal plateau in different deposition periods, and counting the thickness of the corresponding deposition of the medium tidal plateau; calculating a preliminary paleo-tidal range by using the thickness of the middle tidal plateau deposition; acquiring an ancient landform by using a source-sink system; obtaining the paleo-water depth based on the Milankovitch cycle; calculating and counting an offshore distance by using the obtained ancient landform; and calibrating the preliminary paleo-tidal range by using the paleo-landform, the paleo-water depth and the offshore distance to obtain a paleo-tidal range value, and then drawing to obtain a paleo-tidal range plane graph. According to the method provided by the invention, preliminary restoration of the paleo-tidal range is carried out by applying the relationship between the deposition thickness of the tidal plateau and the tidal range, numerical fitting is carried out on the paleo-tidal range, the paleo-landform, the paleo-water depth and the offshore distance of the research area, the paleo-tidal range is corrected, and finally plane restoration of the paleo-tidal range of the research area is realized.

Description

Technical Field

[0001] The present invention relates to the field of paleo-tidal range restoration in a sedimentary environment, and specifically to a paleo-tidal range restoration method, device, equipment, and medium based on the thickness of a mid-tidal flat, and in particular to a paleo-tidal range quantitative restoration method, device, equipment, and medium based on the thickness of a mid-tidal flat. Background Art

[0002] In different geological periods, barrier-lagoon systems are an indispensable part of the coastal environment. Today, barrier systems account for about 15% of the world's ocean coastline. Ancient barrier coastal sedimentary systems contain important coal and oil and gas resources.

[0003] Coastal barriers usually develop parallel to the coastline, located between the nearshore and open marine environments, dominated by tidal action, and include shore, barrier and back-barrier sediments. The back-barrier is composed of swamps, tidal flats and lagoons. Its formation and preservation are controlled by changes in sediment supply, changes in accommodating space and sea level fluctuations.

[0004] The barrier-type coastal sedimentary system mainly includes shore, barrier island, tidal flat and lagoon. Among them, the barrier-lagoon sedimentary system can be further divided into barrier island, lagoon, tidal channel, tidal delta, etc. Compared with the river and delta sedimentary systems, the hydrodynamics of the barrier coastal sedimentary system is more complex. The distribution of sedimentary bodies is controlled and transformed by multi-directional water flows such as tides and waves. Tidal deltas and tidal sand bars are easily developed under the direct control of tidal action. The size of the tidal range directly reflects the strength of the tidal action. The size of the tidal range has a controlling effect on the scale of sand body development. Therefore, the quantitative characterization of the paleo-tidal range is crucial for the study of the barrier-lagoon sedimentary system, which is conducive to predicting the spatial distribution of sand bodies in oil and gas reservoirs.

[0005] Currently, there are many studies on modern tidal range analysis, which mainly use modern sedimentary case dissection, flume physical simulation and numerical simulation. However, modern sedimentary case dissection has great limitations. Rock outcrops can provide the two-dimensional structure of tidal sediments, but the internal three-dimensional structure information provided is insufficient. Physical simulation is limited by many conditions. On the one hand, the research scale is small, and on the other hand, the running time is limited, so that there are large errors in the experimental results. Numerical simulation is highly operational and suitable for various tidal deposition scenarios. It can provide strong evidence for studying the morphological dynamics and stratigraphic patterns of sediment transport, sedimentary erosion and other spatial and temporal scales. At the same time, common numerical simulations can be divided into sediment simulations based on the Navier-Stokes equation and sediment transport process simulations based on the diffusion equation. However, numerical simulations based on sediment transport patterns are highly practical and effective tools for predicting morphological development.

[0006] Based on numerical simulation, current research believes that river flow and tidal amplitude (tidal range) have a significant impact on the deposition range and thickness of tidal sand bodies, but quantitative research has not been achieved. At the same time, the particle size parameters of modern tidal sediment samples were analyzed, and sensitive indicators reflecting tidal dynamic changes were proposed. The relationship between the thickness of modern tidal flat sedimentary phases and the tidal range was established, realizing the quantitative characterization of the tidal range.

[0007] However, due to the influence of factors such as global sea level, climate change and tectonic movement, there are differences between modern tidal ranges and ancient tidal ranges, which makes it impossible to intuitively and scientifically restore the ancient tidal range in the study area, and it is also impossible to clarify the planar distribution of the ancient tidal range during the sedimentary period. It is also not conducive to the study of sand body development, distribution and reservoir characteristics under tidal action, and the distribution of favorable reservoirs. Summary of the invention

[0008] In view of the problems existing in the prior art, the purpose of the present invention is to provide a method, device, equipment and medium for restoring the paleotidal range based on the thickness of the mid-tidal flat, so as to solve the problem that due to the influence of factors such as global sea level, climate change and tectonic movement, there are differences between modern tidal range and paleotidal range, which makes it impossible to intuitively and scientifically restore the paleotidal range of the study area, and it is also impossible to clarify the planar distribution of the paleotidal range during the sedimentation period, which is also not conducive to the study of sand body development, distribution and reservoir characteristics, and favorable reservoir distribution under the action of tides.

[0009] To achieve this object, the present invention adopts the following technical solutions:

[0010] In a first aspect, the present invention provides a method for restoring the paleotidal range based on the thickness of the mid-tidal flat, the method comprising:

[0011] Obtain mid-tidal flat sediments in different sedimentary periods and calculate the thickness of the corresponding mid-tidal flat sediments;

[0012] The thickness of mid-tidal flat sediments was used to calculate the preliminary paleotidal range;

[0013] Using source-sink systems to obtain paleo-geomorphology;

[0014] Obtaining paleo-water depth based on Milankovitch cycles;

[0015] The distance from the shore is calculated and counted using the acquired paleo-geomorphology;

[0016] The paleotidal range was calibrated using paleo-geomorphology, paleo-water depth and offshore distance to obtain the paleotidal difference, which was then plotted to obtain a paleotidal range map.

[0017] The method provided by the present invention performs preliminary restoration of the paleo-tidal range by using the relationship between the mid-tidal flat sediment thickness and the tidal range, numerically fits the paleo-tidal range with the paleo-geomorphology, paleo-water depth and offshore distance of the study area, corrects the paleo-tidal range, and ultimately achieves plane restoration of the paleo-tidal range in the study area.

[0018] As a preferred technical solution of the present invention, the acquisition of mid-tidal flat deposits in different depositional periods includes: confirming the mid-tidal flat deposits based on drilling core analysis and logging phase analysis;

[0019] Preferably, if dry cracks, cross-laminae, veins, lenses and wavy bedding are identified in the drilling core analysis, the corresponding deposits are mid-tidal flat deposits;

[0020] Preferably, in the logging phase analysis, if the logging curve response is a high-amplitude GR, and the curve shape is sawtooth and peak-shaped, the corresponding deposition is a mid-tidal flat deposition.

[0021] Preferably, the counting of the thickness of the corresponding mid-tidal flat deposits includes: determining the range of the mid-tidal flat deposits by logging phase characteristics, determining the thickness of the mid-tidal flat deposits in combination with drilling core analysis, converting the thickness of the mid-tidal flat deposits according to a scale, and counting the thickness of the mid-tidal flat deposits in different periods.

[0022] As a preferred technical solution of the present invention, the method of calculating the preliminary paleotidal range using the thickness of the mid-tidal flat sediments includes: substituting the thickness of the mid-tidal flat sediments into the tidal range calculation formula to obtain the preliminary paleotidal range.

[0023] As a preferred technical solution of the present invention, the use of the source-sink system to obtain paleo-geomorphology includes: on the basis of the source-sink system division, by analyzing the differences in fault activity and the evolution law of sedimentary stratum thickness in different periods, clarifying the differential sedimentation law and performing differential sedimentation correction on the basis of the sedimentary base level, and then restoring the erosion amount of the erosion zone and the over-stripping zone through the principle of conservation of total mass and the formation dip restoration method, and obtaining the paleo-geomorphology framework of the depositional period, and then establishing the correspondence between different watersheds and sedimentary areas on this basis, thereby restoring the paleo-geomorphology of the prototype basin in the depositional period.

[0024] As a preferred technical solution of the present invention, the method of obtaining paleo-water depth based on Milankovitch cycles includes: calculating single-point paleo-water depth using cobalt and lanthanum elements based on lithofacies-sensitive logging data of rock samples, and according to the Milankovitch astronomical cycle theory, obtaining the sedimentation rate of the depth domain spectrum through evolution harmonic power spectrum analysis and Monte Carlo simulation of de-noised logging data to calculate the evolution average spectral mismatch, and then establishing and fitting the sedimentation rate-paleo-water depth formula using the single-point paleo-water depth and sedimentation rate to calculate the paleo-water depth of the entire sequence of wells.

[0025] As a preferred technical solution of the present invention, the method of calculating and statistically analyzing the offshore distance using the acquired paleo-geomorphology includes: converting the paleo-geomorphology into elevation data, and then calculating and statistically analyzing the offshore distance in combination with the high tide line distribution.

[0026] As a preferred technical solution of the present invention, the method of calibrating the preliminary paleotidal range using paleo-geomorphology, paleo-water depth and offshore distance to obtain the paleotidal range value and then drawing the paleotidal range plane map comprises:

[0027] The correlation between the preliminary paleotidal range and paleotopography, paleowater depth and offshore distance was analyzed respectively, and the paleotidal range was corrected according to the relationship between the mid-tidal flat sediment thickness and paleowater depth and offshore distance to obtain the paleotidal difference. Then, the paleotidal difference was used to draw a plane contour map to obtain a paleotidal range plane map.

[0028] In a second aspect, the present invention provides a paleotidal range restoration device based on the thickness of a mid-tidal flat, the paleotidal range restoration device comprising:

[0029] The module for obtaining the thickness of the mid-tidal flat sediments is used to obtain the mid-tidal flat sediments in different sedimentary periods and to calculate the thickness of the corresponding mid-tidal flat sediments;

[0030] The module for obtaining the preliminary paleotidal range is used to calculate the preliminary paleotidal range using the thickness of the mid-tidal flat sediments;

[0031] Paleomorphology acquisition module, used to obtain paleomorphology using source-sink system;

[0032] Paleo-water depth acquisition module, used to obtain paleo-water depth based on Milankovitch cycles;

[0033] The offshore distance acquisition module is used to calculate and count the offshore distance using the acquired paleo-geomorphology;

[0034] The paleotidal difference acquisition module is used to calibrate the preliminary paleotidal difference using paleo-geomorphology, paleo-water depth and offshore distance to obtain the paleotidal difference, and then draw a paleotidal difference plan.

[0035] In a third aspect, the present invention provides an electronic device, the electronic device comprising:

[0036] at least one processor; and a memory communicatively coupled to the at least one processor;

[0037] The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the paleotidal range restoration method as described in the first aspect.

[0038] In a fourth aspect, the present invention provides a computer storage medium, wherein the computer storage medium stores computer executable instructions, and when the computer executable instructions are executed by a processor, the paleotidal range recovery method as described in the first aspect is implemented.

[0039] Compared with the prior art solutions, the present invention has the following beneficial effects:

[0040] (1) Based on the statistics of the thickness of the mid-tidal flat, the present invention uses the relationship between modern tidal flat sedimentation and tidal range to perform preliminary quantitative restoration of the tidal range, and then further corrects it according to the fitting results of the paleo-tidal range, paleo-geomorphology, paleo-water depth and offshore distance, and finally realizes the restoration of the paleo-tidal range plane. The restoration of the paleo-tidal range provides a new idea and method for the study of barrier-lagoon sedimentary systems, and has the characteristics of advancement, exploration practicality and wide applicability, and can be widely used in sedimentological research and petroleum exploration.

[0041] (2) Based on the quantitative restoration method of paleotidal range based on the thickness of the mid-tidal flat, the strength and scope of tidal effects in the study area are clarified, laying a foundation for the study of barrier-lagoon sedimentary systems. The strength of paleotidal range in different structural units in the study area and its temporal and spatial evolution can be intuitively presented in the paleotidal range plane diagram, which has a predictive effect on the distribution of subsequent tidal sand bodies and thus guides oil exploration. The use of core data for paleotidal range restoration is a highly operational, reliable, simple and intuitive method that can solve the problems of barrier-tidal flat sedimentary system research and favorable reservoir prediction involved in oil and gas exploration, and has a wide range of uses and broad market prospects. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] Figure 1 It is a flow chart of a method for restoring the paleotidal range based on the thickness of the mid-tidal flat provided by an embodiment of the present invention;

[0043] Figure 2 is a schematic diagram of a paleotidal range recovery device based on the thickness of a mid-tidal flat provided by an embodiment of the present invention;

[0044] Figure 3 is a schematic diagram of an electronic device provided by an embodiment of the present invention;

[0045] Figure 4 is a drilling core identification map of the mid-tidal flat deposits in Example 1 of the present invention;

[0046] Figure 5 It is a 3400-4600m logging curve identification diagram of the mid-tidal flat deposits in Example 1 of the present invention;

[0047] Figure 6 It is a 4600-5000m logging curve identification diagram of the mid-tidal flat deposits in Example 1 of the present invention;

[0048] Figure 7 is the result relationship between the mid-tidal flat sediment thickness and the preliminary paleo-tidal range of a single well obtained in Example 1 of the present invention;

[0049] Figure 8 is the ancient landform restoration map obtained in Example 1 of the present invention;

[0050] Fig. 9 This is the paleo-water depth restoration map obtained in Example 1 of the present invention;

[0051] Fig.10 is a statistical diagram of the offshore distance obtained in Example 1 of the present invention;

[0052] Fig.11 is a schematic diagram of the correlation between paleo-geomorphology and paleo-water depth in Example 1 of the present invention;

[0053] Fig.12 This is a schematic diagram of the correlation between the paleo-tidal range and the paleo-water depth in Example 1 of the present invention;

[0054] Fig.13 is a schematic diagram of the correlation between the paleo-tidal range and the offshore distance in Example 1 of the present invention;

[0055] Fig.14 This is the plane restoration map of the paleotidal range obtained in Example 1 of the present invention.

[0056] The present invention is further described in detail below. However, the following examples are only simplified examples of the present invention and do not represent or limit the scope of protection of the present invention. The scope of protection of the present invention shall be subject to the claims. DETAILED DESCRIPTION

[0057] To better illustrate the present invention and facilitate understanding of the technical solution of the present invention, typical but non-limiting embodiments of the present invention are as follows:

[0058] This embodiment provides a method for restoring the paleo-tidal range based on the thickness of the mid-tidal flat. The flow chart is as follows: Figure 1 As shown, the details are as follows:

[0059] Obtain mid-tidal flat sediments in different sedimentary periods and calculate the thickness of the corresponding mid-tidal flat sediments;

[0060] The thickness of mid-tidal flat sediments was used to calculate the preliminary paleotidal range;

[0061] Using source-sink systems to obtain paleo-geomorphology;

[0062] Obtaining paleo-water depth based on Milankovitch cycles;

[0063] The distance from the shore is calculated and counted using the acquired paleo-geomorphology;

[0064] The paleotidal range was calibrated using paleo-geomorphology, paleo-water depth and offshore distance to obtain the paleotidal difference, which was then plotted to obtain a paleotidal range map.

[0065] Wherein, obtaining mid-tidal flat deposits in different depositional periods includes: confirming the mid-tidal flat deposits based on drilling core analysis and logging phase analysis.

[0066] If dry cracks, interlaced laminae, vein-like, lenticular and wavy bedding are identified in the drilling core analysis, the corresponding deposits are mid-tidal flat deposits.

[0067] Wherein, in the logging phase analysis, if the logging curve response is a high-amplitude GR, and the curve shape is sawtooth and peak-shaped, the corresponding deposition is a mid-tidal flat deposition.

[0068] The statistics of the thickness of the corresponding mid-tidal flat deposits include: determining the range of the mid-tidal flat deposits by logging phase characteristics, determining the thickness of the mid-tidal flat deposits in combination with drilling core analysis, converting the thickness of the mid-tidal flat deposits according to the scale, and counting the thickness of the mid-tidal flat deposits in different periods.

[0069] The method of calculating the preliminary paleotidal range using the thickness of the mid-tidal flat sediments comprises: substituting the thickness of the mid-tidal flat sediments into the tidal range calculation formula to obtain the preliminary paleotidal range.

[0070] Among them, the tidal range calculation formula is to find indicators that sensitively reflect the changes in tidal flat dynamics through the grain size parameter analysis of modern tidal flat columnar sediment samples, to divide the sedimentary phases, and to obtain the relationship between the thickness of the tidal flat sedimentary phase and the tidal range based on the grain size digital characteristics. For example, according to the existing technology (Guo Yanxia, ​​Fan Daidu, Zhao Juan. Grain size characteristics and sedimentary phase division of tidal flat sequences - taking Andong Shoal in Hangzhou Bay as an example [J]. Frontiers of Marine Geology, 2004, 20(5): 9-14. DOI: 10.3969 / j.issn.1009-2722.2004.05.003.) the tidal range calculation formula is obtained.

[0071] Among them, the use of the source-sink system to obtain paleo-geomorphology includes: on the basis of the source-sink system division, by analyzing the differences in fault activity and the evolution law of sedimentary stratum thickness in different periods, clarifying the differential sedimentation law and making differential sedimentation correction based on the sedimentary base level, and then restoring the erosion amount of the erosion zone and the over-stripping zone through the principle of conservation of total mass and the formation dip restoration method, and obtaining the paleo-geomorphology framework of the depositional period, and then establishing the correspondence between different watersheds and sedimentary areas on this basis, and then restoring the paleo-geomorphology of the prototype basin in the depositional period. Specifically, the stratigraphic interface is obtained using seismic data, the planation surface is selected as the sedimentary datum, the sedimentary datum is flattened, the differential settlement in the late stage of deposition is corrected, the erosion zone, overdenudation zone (area where overlap and erosion coexist) and overlap zone are identified according to the seismic reflection characteristics, the erosion amount in the overdenudation zone is restored by the stratigraphic trend method, and the erosion amount in the erosion zone is restored by the stratigraphic trend method and volume inversion method, the differential settlement and the erosion amount in the erosion zone and overdenudation zone are integrated to restore the paleo-geomorphic pattern of the depositional period, and finally, the provenance-sedimentary area correspondence of different structural-sedimentary units is established on the basis of the source-sink system division, and the restoration of the prototype basin in the depositional period is improved.

[0072] Among them, the method of obtaining paleowater depth based on Milankovitch cycle includes: calculating single-point paleowater depth based on lithofacies sensitive logging data of rock samples using cobalt and lanthanum elements, and according to Milankovitch astronomical cycle theory, obtaining the sedimentation rate of the depth domain spectrum through evolution harmonic power spectrum analysis and Monte Carlo simulation of noise-reduced logging data to calculate the evolution average spectral mismatch, and then using the single-point paleowater depth and sedimentation rate to establish and fit the sedimentation rate-paleowater depth formula to calculate the paleowater depth of the entire sequence of wells.

[0073] Among them, the abundance of cobalt and lanthanum elements was obtained by analyzing and testing the rock samples of a single well, and the single-point sedimentation rate and paleo-water depth were calculated using cobalt and lanthanum elements. The sedimentation rate-paleo-water depth formula was established and fitted using the single-point paleo-water depth and sedimentation rate to calculate the paleo-water depth of the entire sequence of wells, and then restore the temporal and spatial distribution of paleo-water depth. Sedimentation rate V S The calculation formula is:

[0074]

[0075] Where V0 is the sedimentation rate of modern lakes under approximate conditions, N Co is the abundance of cobalt in normal lake sediments, S Co is the abundance of cobalt in the sample, T Co is the abundance of cobalt in terrigenous clastic rocks, and t is the ratio of the lanthanum content in the sample to the average abundance of lanthanum in terrigenous clastic rocks.

[0076] The calculation formula for single-point paleo-water depth is:

[0077]

[0078] Where h is the single point paleo-water depth, V S is the deposition rate.

[0079] The method of calculating and counting the offshore distance using the acquired paleo-geomorphology includes: converting the paleo-geomorphology into elevation data, and then calculating and counting the offshore distance in combination with the high tide line distribution.

[0080] The paleotidal difference is obtained by calibrating the preliminary paleotidal difference using paleo-geomorphology, paleo-water depth and offshore distance, and includes:

[0081] The correlation between the preliminary paleotidal range and paleotopography, paleowater depth and offshore distance was analyzed respectively, and the paleotidal range was corrected according to the relationship between the mid-tidal flat sediment thickness and paleowater depth and offshore distance to obtain the paleotidal difference. Then, the paleotidal difference was used to draw a plane contour map to obtain a paleotidal range plane map.

[0082] Specifically, the restored paleo-geomorphic elevation data, paleo-water depth, offshore distance and paleo-tidal range of multiple single wells from the same period were subjected to correlation analysis to determine the paleo-tidal range.

[0083] In the present invention, for blocks with scarce paleo-tidal range data, corrections can be made based on paleo-geomorphology and paleo-water depth.

[0084] In the present invention, in the process of paleo-tidal range restoration, calculation of offshore distance and acquisition of paleo-tidal range plane distribution map and the like, conventional processing software in the art can be selected for auxiliary implementation, such as Petrel software, Global Mapper software, Shuanghu, ArcGIS software and the like.

[0085] Furthermore, the present invention provides a paleotidal range recovery device based on the thickness of the mid-tidal flat, such as Figure 2 As shown, the paleotidal range recovery device comprises:

[0086] The module for obtaining the thickness of the mid-tidal flat sediments is used to obtain the mid-tidal flat sediments in different sedimentary periods and to calculate the thickness of the corresponding mid-tidal flat sediments;

[0087] The module for obtaining the preliminary paleotidal range is used to calculate the preliminary paleotidal range using the thickness of the mid-tidal flat sediments;

[0088] Paleomorphology acquisition module, used to obtain paleomorphology using source-sink system;

[0089] Paleo-water depth acquisition module, used to obtain paleo-water depth based on Milankovitch cycles;

[0090] The offshore distance acquisition module is used to calculate and count the offshore distance using the acquired paleo-geomorphology;

[0091] The paleotidal difference acquisition module is used to calibrate the preliminary paleotidal difference using paleo-geomorphology, paleo-water depth and offshore distance to obtain the paleotidal difference, and then draw a paleotidal difference plan.

[0092] Regarding the device in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0093] Further, the present invention provides an electronic device, which is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device may also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present invention described and / or required herein.

[0094] like Figure 3As shown, the electronic device 10 includes at least one processor 11, and a memory connected to the at least one processor 11, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc., wherein the memory stores a computer program that can be executed by at least one processor, and the processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 to the random access memory (RAM) 13. In the RAM 13, various programs and data required for the operation of the electronic device 10 can also be stored. The processor 11, the ROM 12, and the RAM 13 are connected to each other through a bus 14. The I / O interface 15 is also connected to the bus 14.

[0095] A number of components in the electronic device 10 are connected to the I / O interface 15, including: an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices through a computer network such as the Internet and / or various telecommunication networks.

[0096] The processor 11 may be a variety of general and / or special processing components with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any appropriate processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as a paleotidal range recovery method based on the thickness of the mid-tidal flat.

[0097] In some embodiments, the paleotidal range restoration method based on the thickness of the mid-tidal flat may be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as a storage unit 18. In some embodiments, part or all of the computer program may be loaded and / or installed on the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the paleotidal range restoration method based on the thickness of the mid-tidal flat described above may be performed. Alternatively, in other embodiments, the processor 11 may be configured to execute the paleotidal range restoration method based on the thickness of the mid-tidal flat by any other appropriate means (e.g., by means of firmware).

[0098] Various implementations of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), systems on chips (SOCs), load programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various implementations can include: being implemented in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0099] Computer programs for implementing the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, so that when the computer program is executed by the processor, the functions / operations specified in the flow chart and / or block diagram are implemented. The computer program may be executed entirely on the machine, partially on the machine, partially on the machine and partially on a remote machine as a stand-alone software package, or entirely on a remote machine or server.

[0100] In the context of the present invention, a computer-readable storage medium may be a tangible medium that may contain or store a computer program for use by or in combination with an instruction execution system, device or equipment. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices or equipment, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. A more specific example of a machine-readable storage medium may include an electrical connection based on one or more lines, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0101] To provide interaction with a user, the systems and techniques described herein may be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and a pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices may also be used to provide interaction with the user; for example, the feedback provided to the user may be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form (including acoustic input, voice input, or tactile input).

[0102] The systems and techniques described herein may be implemented in a computing system that includes backend components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes frontend components (e.g., a user computer with a graphical user interface or a web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such backend components, middleware components, or frontend components. The components of the system may be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0103] A computing system may include a client and a server. The client and the server are generally remote from each other and usually interact through a communication network. The client and server relationship is generated by computer programs running on the corresponding computers and having a client-server relationship with each other. The server may be a cloud server, also known as a cloud computing server or cloud host, which is a host product in the cloud computing service system to solve the defects of difficult management and weak business scalability in traditional physical hosts and VPS services.

[0104] The server provided in this embodiment includes: a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the above-mentioned paleo-tidal range recovery method based on the thickness of the mid-tidal flat is implemented.

[0105] Unless otherwise specifically stated, terms such as processing, computing, calculating, determining, displaying, etc. may refer to the actions and / or processes of one or more processing or computing systems, or similar devices, which operate and convert data represented as physical (e.g., electronic) quantities within registers or memories of a processing system into other data similarly represented as physical quantities within memories, registers, or other such information storage, transmission, or display devices of the processing system. Information and signals may be represented using any of a variety of different techniques and methods. For example, data, instructions, commands, information, signals, bits, symbols, and chips mentioned throughout the above description may be represented by voltages, currents, electromagnetic waves, magnetic fields or particles, light fields or particles, or any combination thereof.

[0106] Those skilled in the art will also appreciate that the various illustrative logic blocks, modules, circuits, and algorithm steps described in conjunction with the embodiments of the present invention may all be implemented as electronic hardware, computer software, or a combination thereof. In order to clearly illustrate the interchangeability between hardware and software, the various illustrative components, blocks, modules, circuits, and steps described above are generally described around their functions. Whether such functions are implemented as hardware or software depends on the specific application and the design constraints imposed on the entire system. A skilled person may implement the described functions in an alternative manner for each specific application, but such implementation decisions should not be interpreted as departing from the scope of protection of the present disclosure.

[0107] The steps of the method or algorithm described in conjunction with the embodiments herein may be directly embodied as hardware, a software module executed by a processor, or a combination thereof. The software module may be located in a RAM memory, a flash memory, a ROM memory, an EPROM memory, an EEPROM memory, a register, a hard disk, a mobile disk, a CD-ROM, or any other form of storage medium well known in the art. An exemplary storage medium is connected to the processor so that the processor can read information from the storage medium and can write information to the storage medium. Of course, the storage medium may also be an integral part of the processor. The processor and the storage medium may be located in an ASIC. The ASIC may be located in a user terminal. Of course, the processor and the storage medium may also be present in a user terminal as discrete components.

[0108] For software implementation, the techniques described in the present invention can be implemented with modules (e.g., procedures, functions, etc.) that perform the functions described in the present application. These software codes can be stored in a memory unit and executed by a processor. The memory unit can be implemented within the processor or outside the processor. In the latter case, it is coupled to the processor in a communication manner via various means, which are well known in the art.

[0109] Further, in order to illustrate the recognition effect that can be achieved by the paleotidal range restoration method based on the thickness of the mid-tidal flat provided by the present invention, an actual example is used for illustration, as follows:

[0110] Example 1

[0111] This embodiment provides a method for restoring the paleotidal range based on the thickness of the mid-tidal flat. The specific process is as follows:

[0112] 1) Based on drilling core analysis and logging phase analysis, identify tidal flat deposits in different sedimentary periods

[0113] Based on the drilling core, the lithology of the mid-tidal flat sediments is mainly silt-fine sandstone, and is mostly transported by bed sand load and suspended load. Through the drilling core, typical sedimentary structures such as dry cracks, cross-lamina, veins, lenses, and wavy bedding can be identified ( Figure 4 , a is W1 well, 3473.01-3486.9m, wavy bedding, b is W2 well, 4563.96-4565.07m, wavy bedding, c is W3 well, 3960.28-3961.30m, lenticular bedding, d is W4 well, 4083.96m, feathery bedding, e is W5 well, 3651.7-3651.91m, vein-like bedding), combined with the logging curve response of the middle tidal flat sedimentation, the curve shape is sawtooth, sharp peak type ( Figure 5 and Figure 6 ) and then identify tidal flat deposits in different sedimentary periods.

[0114] 2) Through drilling core analysis and logging phase analysis, the thickness of tidal flat sediments in different periods was counted

[0115] The core analysis is mainly based on the well logging phase analysis, and the thickness of the mid-tidal flat sediments is converted according to the map scale. The thickness of the mid-tidal flat sediments in different periods is statistically analyzed (the results are shown in Figure 2). Figure 7 ).

[0116] 3) Preliminary calculation of the initial tidal range using the modern tidal flat-tidal range relationship

[0117] Substitute the thickness of the mid-tidal flat sediment into the tidal range calculation formula: y = 1.8314x + 157.31 (cm), where x is the thickness of the mid-tidal flat sediment and y is the initial tidal range value, and calculate the initial tidal range of a single well (see the results). Figure 7 ).

[0118] 4) “Three-area and five-step” paleo-geomorphological restoration based on source-sink system

[0119] Paleogeomorphology restoration is based on the division of source-sink system. By analyzing the differences in fault activity and the evolution law of sedimentary stratum thickness in different periods, the differential sedimentation law and sedimentary base level are clarified to make differential sedimentation correction. The erosion amount of erosion zone and over-stripping zone is restored through the principle of total material conservation and stratum dip restoration method, and the paleogeomorphology framework of the sedimentary period is improved. On this basis, the corresponding relationship between different basins and sedimentary areas is established to restore the paleogeomorphology of the prototype basin during the sedimentary period (the results are shown in Figure 8 ).

[0120] Specifically, the stratigraphic interface is obtained using seismic data, the planation surface is selected as the sedimentary datum, the sedimentary datum is flattened, the differential settlement in the late stage of deposition is corrected, the erosion zone, overdenudation zone (area where overlap and erosion coexist) and overlap zone are identified according to the seismic reflection characteristics, the erosion amount in the overdenudation zone is restored by the stratigraphic trend method, and the erosion amount in the erosion zone is restored by the stratigraphic trend method and volume inversion method, the differential settlement and the erosion amount in the erosion zone and overdenudation zone are integrated to restore the paleo-geomorphic pattern of the depositional period, and finally, the provenance-sedimentary area correspondence of different structural-sedimentary units is established on the basis of the source-sink system division, and the restoration of the prototype basin in the depositional period is improved.

[0121] 5) Quantitative restoration of paleo-water depth based on Milankovitch cycles

[0122] The single-point paleo-water depth was calculated using cobalt and lanthanum elements. According to the Milankovitch astronomical cycle theory, the sedimentation rate of the depth domain spectrum was obtained by analyzing the evolution harmonic power spectrum of the de-noised logging data and calculating the evolution average spectrum mismatch through Monte Carlo simulation. The sedimentation rate-paleo-water depth formula was established and fitted using the single-point paleo-water depth and sedimentation rate to calculate the paleo-water depth of the entire sequence of wells, and then restore the temporal and spatial distribution of the paleo-water depth (results are shown in Fig. 9 ).

[0123] 6) Calculate and count the distance from the shore using the restored paleo-geomorphology

[0124] The paleo-geomorphology restored in the Petrel software was exported as elevation data, and this elevation data was imported into the GlobalMapper software. Combined with the high tide line distribution, the offshore distance of a single well was calculated and counted (the results are shown in Fig.10 ).

[0125] 7) Use paleo-geomorphology, paleo-water depth, offshore distance and tidal range data for fitting correction

[0126] The restored tidal range was used to establish correlation analysis with paleo-geomorphology, paleo-water depth and offshore distance, and paleo-tidal range correction was performed to determine the paleo-tidal difference value (the corrected curve is shown in Fig.11 , Fig.12 and Fig.13), specifically, the restored paleogeomorphic elevation data, paleowater depth, offshore distance and paleotidal range of multiple single wells in the same period were subjected to correlation analysis to determine the paleotidal difference. The correlation analysis is referred to Example 1 (the corrected curve is shown in Fig.11 , Fig.12 and Fig.13 ). The correlation between paleo-geomorphic elevation and paleo-water depth is: y = -0.0989x + 342.93, where x is the paleo-geomorphic elevation and y is the paleo-water depth; the correlation between paleo-tidal range and paleo-water depth is: y = 8.799x - 11.044, where x is the paleo-tidal range and y is the paleo-water depth; the paleo-tidal range and offshore distance is: y = 0.2841x, where x is the offshore distance and y is the paleo-tidal range.

[0127] 8) Using the scattered values ​​of the paleo-tidal range, a planar distribution map of the paleo-tidal range is presented

[0128] The paleotidal difference point data were imported into the Petrel software to draw the plane contour map and realize the plane restoration of the paleotidal difference in the study area (the results are shown in Fig.14 ).

[0129] The preferred embodiments of the present invention are described in detail above. However, the present invention is not limited to the specific details in the above embodiments. Within the technical concept of the present invention, a variety of simple modifications can be made to the technical solution of the present invention, and these simple modifications all belong to the protection scope of the present invention.

[0130] It should also be noted that the various specific technical features described in the above specific embodiments can be combined in any suitable manner without contradiction. In order to avoid unnecessary repetition, the present invention will not further describe various possible combinations.

[0131] In addition, various embodiments of the present invention may be arbitrarily combined, and as long as they do not violate the concept of the present invention, they should also be regarded as the contents disclosed by the present invention.

Claims

1. A method for restoring the paleotidal range based on the thickness of the mid-tidal flat, characterized in that: The paleotidal range restoration method comprises: Obtain mid-tidal flat sediments in different sedimentary periods and calculate the thickness of the corresponding mid-tidal flat sediments; The thickness of mid-tidal flat sediments was used to calculate the preliminary paleotidal range; Using source-sink systems to obtain paleo-geomorphology; Obtaining paleo-water depth based on Milankovitch cycles; The distance from the shore is calculated and counted using the acquired paleo-geomorphology; The paleotidal range was calibrated using paleo-geomorphology, paleo-water depth and offshore distance to obtain the paleotidal difference, which was then plotted to obtain a paleotidal range map.

2. The method for restoring the paleotidal range according to claim 1, characterized in that: The obtaining of mid-tidal flat deposits in different depositional periods comprises: confirming the mid-tidal flat deposits based on drilling core analysis and logging phase analysis; Preferably, if dry cracks, cross-laminae, veins, lenses and wavy bedding are identified in the drilling core analysis, the corresponding deposits are mid-tidal flat deposits; Preferably, in the logging phase analysis, if the logging curve response is a high amplitude GR, and the curve shape is sawtooth and peak-shaped, the corresponding deposition is a mid-tidal flat deposition; Preferably, the counting of the thickness of the corresponding mid-tidal flat deposits comprises: determining the range of the mid-tidal flat deposits by logging phase characteristics, determining the thickness of the mid-tidal flat deposits in combination with drilling core analysis, converting the thickness of the mid-tidal flat deposits according to a scale, and counting the thickness of the mid-tidal flat deposits in different periods.

3. The method for restoring the paleotidal range according to claim 1, characterized in that: The method of calculating the preliminary paleotidal range using the thickness of the mid-tidal flat sediments comprises: substituting the thickness of the mid-tidal flat sediments into a tidal range calculation formula to obtain the preliminary paleotidal range.

4. The method for restoring the paleotidal range according to claim 1, characterized in that: The method of obtaining paleo-geomorphology by using the source-sink system includes: based on the division of the source-sink system, analyzing the differences in fault activities and the evolution law of sedimentary stratum thickness in different periods, clarifying the differential sedimentation law and making differential sedimentation correction based on the sedimentary base level, and then restoring the erosion amount of the erosion zone and the over-stripping zone through the principle of conservation of total mass and the formation dip restoration method, and obtaining the paleo-geomorphology framework of the depositional period, and then establishing the correspondence between different watersheds and sedimentary areas on this basis, and then restoring the paleo-geomorphology of the prototype basin in the depositional period.

5. The method for restoring the paleotidal range according to claim 1, characterized in that: The method of obtaining paleo-water depth based on Milankovitch cycles includes: calculating single-point paleo-water depth using cobalt and lanthanum elements based on lithofacies-sensitive logging data of rock samples, and according to the Milankovitch astronomical cycle theory, calculating the evolutionary average spectral mismatch through evolution harmonic power spectrum analysis and Monte Carlo simulation of noise-reduced logging data to obtain the sedimentation rate of the depth domain spectrum, and then establishing and fitting the sedimentation rate-paleo-water depth formula using the single-point paleo-water depth and sedimentation rate to calculate the paleo-water depth of the entire sequence of wells.

6. The method for restoring the paleotidal range according to claim 1, characterized in that: The method of calculating and counting the offshore distance using the acquired paleo-geomorphology includes: converting the paleo-geomorphology into elevation data, and then calculating and counting the offshore distance in combination with the high tide line distribution.

7. The method for restoring the paleotidal range according to claim 1, characterized in that: The method of calibrating the preliminary paleotidal range using paleo-geomorphology, paleo-water depth and offshore distance to obtain the paleotidal range value and then drawing the paleotidal range plane map includes: The correlation between the preliminary paleotidal range and paleotopography, paleowater depth and offshore distance was analyzed respectively, and the paleotidal range was corrected according to the relationship between the mid-tidal flat sediment thickness and paleowater depth and offshore distance to obtain the paleotidal difference. Then, the paleotidal difference was used to draw a plane contour map to obtain a paleotidal range plane map.

8. A paleotidal range recovery device based on the thickness of the mid-tidal flat, characterized in that: The paleotidal range recovery device comprises: The module for obtaining the thickness of the mid-tidal flat sediments is used to obtain the mid-tidal flat sediments in different sedimentary periods and to calculate the thickness of the corresponding mid-tidal flat sediments; The module for obtaining the preliminary paleotidal range is used to calculate the preliminary paleotidal range using the thickness of the mid-tidal flat sediments; Paleomorphology acquisition module, used to obtain paleomorphology using source-sink system; Paleo-water depth acquisition module, used to obtain paleo-water depth based on Milankovitch cycles; The offshore distance acquisition module is used to calculate and count the offshore distance using the acquired paleo-geomorphology; The paleotidal difference acquisition module is used to calibrate the preliminary paleotidal difference using paleo-geomorphology, paleo-water depth and offshore distance to obtain the paleotidal difference, and then draw a paleotidal difference plan.

9. An electronic device, characterized in that: The electronic device comprises: at least one processor; and a memory communicatively coupled to the at least one processor; The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the paleotidal range restoration method described in any one of claims 1 to 7.

10. A computer storage medium, characterized in that: The computer storage medium stores computer executable instructions, and when the computer executable instructions are executed by a processor, the paleotidal range restoration method described in any one of claims 1 to 7 is implemented.

Images