Ancient landform recovery method and system based on balanced profile structure recovery

Abstract

The invention belongs to the technical field of sedimentary basin analysis, and particularly relates to an ancient landform restoration method and system based on equilibrium profile structure restoration. The method comprises the following steps: acquiring drilling data of a target area, and establishing an isochronous stratigraphic framework based on the drilling data of the target area; according to the isochronous stratigraphic framework, determining a target stratum bottom interface and a maximum ocean flooding surface of the target area; determining a first stratum deposition true thickness of the target area according to a target stratum bottom interface and a maximum ocean flooding surface of the target area; the first stratum deposition true thickness is the stratum deposition true thickness of the target area under the in-situ coordinates. According to the balance profile structure recovery principle and method, the stratum thickness of the target layer is recovered to the real deposition thickness and returned to the original position in the deposition process, then the original ancient landform in the deposition process is recovered, the reasonability and accuracy of sedimentary facies research and reservoir sand body research are improved, and the exploration success rate is increased.

Description

Technical Field

[0001] This disclosure belongs to the field of sedimentary basin analysis technology, and specifically relates to a paleogeographic restoration method and system based on balanced profile structural restoration. Background Technology

[0002] Paleogeography is an important factor controlling the sedimentary facies, reservoir sand body development and distribution in oil and gas basins, and to a certain extent controls the reservoir-seal combination and oil reservoir distribution. Therefore, paleogeographic research and paleogeographic reconstruction are important research contents in oil and gas exploration.

[0003] There are three main categories of paleogeographic restoration methods: The first category is based on the thickness of strata, such as the residual thickness method (Dai Jinyou, 2005) and the impression method (Yuan Zhen, 2013); the second category is based on sedimentary sequence stratigraphy, such as Zhao Junxing (2001) using vertical and horizontal lithological changes to qualitatively identify paleogeographic features, and Chen Hongde (2005) using high-resolution sequence stratigraphy to restore paleogeographic features; the third category is based on the flattening of geological and seismic layers, such as the balanced profile restoration method (Zhang Hui, 2018), the seismic layer flattening method (Du Yang, 2019), the geological layer flattening method (Yan Haijun, 2016), seismic attribute analysis and seismic slicing method (Chen Hui, 2016), and the double-interface seismic layer flattening method (Liu Yongtao, 2019), etc. The core of the above paleogeographic restoration methods is to obtain a reference surface through analysis and estimate the stratigraphic thickness between the reference surface and the target layer. This stratigraphic thickness is then corrected through a series of adjustments to obtain the true thickness of the sedimentary layer. The paleogeographic features at the time of deposition are then obtained by mirroring or imprinting the true thickness of the sedimentary layer. Correction methods include: erosion correction (Xian Benzhong, 2017), compaction correction (Guo Qiulin, 1998), residual correction (Deng Xingliang, 202), and paleowater depth correction (Jiang Zhenglong, 2009). Dou Luxing (2018) and Zheng Rongcai (2007) analyzed the influence of tectonics on paleogeographic features and performed fault removal and fold removal analyses on the stratigraphic thickness. However, all these thickness estimates and corrections were performed on the basis of in-situ stratigraphy and did not restore the original position of the sedimentary layer. The restored paleogeographic features deviate in position and are deformed in shape compared to the original paleogeographic features. Summary of the Invention

[0004] To address the above problems, this disclosure provides a paleogeographic restoration method based on balanced profile structural restoration, the method comprising:

[0005] Obtain drilling data for the target area and establish an isochronous stratigraphic framework based on the drilling data for the target area;

[0006] The target stratum bottom interface and maximum flood level of the target area are determined based on the isochronous stratigraphic framework.

[0007] The first true thickness of the sedimentary layer in the target area is determined based on the bottom interface of the target layer and the maximum flooding surface; the first true thickness of the sedimentary layer is the true thickness of the sedimentary layer in the target area under in-situ coordinates.

[0008] The correction amount of the new coordinates relative to the original coordinates is determined, and the paleogeographic map of the target area is determined based on the first and second true thicknesses of the sedimentary layers in the target area; the second true thickness of the sedimentary layers is the true thickness of the sedimentary layers in the target area under the new coordinates.

[0009] Preferably, before establishing an isochronous stratigraphic framework based on drilling data from the target area, the process includes: preprocessing the drilling data;

[0010] The preprocessing includes:

[0011] Perform single-well vertical consistency processing on the target area;

[0012] Alternatively, the target area can be subjected to multi-well plane standardization processing.

[0013] Preferably, establishing an isochronous stratigraphic framework based on drilling data from the target area includes:

[0014] Several target wells are selected consecutively based on the sediment source direction of the target area;

[0015] The target well is divided into three-level sequence stratigraphy, and an isochronous stratigraphic framework is established; the isochronous stratigraphic framework includes the target stratigraphic framework, the isochronous stratigraphic framework, and the sequence stratigraphic framework.

[0016] Preferably, establishing an isochronous stratigraphic framework based on drilling data from the target area further includes:

[0017] Based on the relationship between rock porosity and depth in the target area, the rock compaction coefficient of the target area is determined;

[0018] And / or, determine the core ripple wavelength of the target layer in the target region.

[0019] Preferably, determining the target layer bottom interface and maximum flooding level of the target area based on the isochronous stratigraphic framework includes:

[0020] The target layer and the maximum flooding surface of each target well in the isochronous stratigraphic framework are marked onto the seismic profile reflection phase axis of the target area;

[0021] The type of the seismic profile of the target area is adjusted to a lithological profile, and the target layer bottom interface and the maximum flooding surface of the target area are determined based on the reflection characteristics of the target layer bottom interface and the maximum flooding surface on the reflection phase axis in the seismic profile.

[0022] Preferably, determining the true thickness of the first stratigraphic deposit in the target area based on the bottom interface of the target layer and the maximum flooding surface includes:

[0023] The temporal apparent thickness of the stratigraphy in the target area is determined based on the bottom interface of the target layer and the first maximum flooding surface.

[0024] The time-domain apparent thickness of the target area is converted to depth to obtain the depth-domain apparent thickness of the target area.

[0025] Terrain correction is performed on the apparent depth thickness of the target region to obtain the true depth thickness of the target region.

[0026] The depth domain true thickness of the target area is corrected to determine the first stratum sedimentary true thickness of the target area.

[0027] Preferably, determining the temporal apparent thickness of the stratigraphy of the target area based on the bottom interface of the target layer and the first maximum flooding surface includes:

[0028] The seismic interpretation interface of the target region at the target layer is extended to the erosion zone of the target region at the target layer;

[0029] The temporal apparent thickness of the target area is determined based on the difference between the bottom interface of the target layer and the first maximum flood level.

[0030] Preferably, the time-domain apparent thickness of the target area is converted to depth to obtain the depth-domain apparent thickness of the target area, including:

[0031] Construct the target layer velocity field of the target region;

[0032] The depth domain apparent thickness of the target region is determined based on the target layer velocity field and the time domain apparent thickness of the strata.

[0033] Preferably, terrain correction is performed on the apparent depth thickness of the target region to obtain the true depth thickness of the target region, including:

[0034] Extract the seismic dip angle attribute of the maximum flood surface in the target area;

[0035] Based on the earthquake dip angle attribute, the apparent depth thickness of the target area is dip-corrected to determine the true depth thickness of the target area.

[0036] Preferably, the depth domain true thickness of the target region is corrected, including:

[0037] Determine the residual correction amount, and perform residual correction on the true thickness of the depth domain based on the residual correction amount;

[0038] The compaction thickness is determined based on the rock compaction coefficient of the target area, and the true thickness of the depth domain is compacted and corrected based on the compaction thickness.

[0039] The anomalous body thickness correction amount is determined based on the deposition distance of the target area, and the anomalous body thickness is corrected for the true thickness of the depth domain based on the anomalous body thickness correction amount.

[0040] Preferably, determining the correction amount of the new coordinates relative to the original coordinates includes:

[0041] Construct an equilibrium evolution profile of the target region, and determine the pre-recovery profile line on the equilibrium evolution profile;

[0042] The equilibrium evolution profile of the target region is restored, and the restored profile line is redefined.

[0043] The amount of scatter point expansion is determined by comparing the profile lines before and after restoration.

[0044] The scattered point scaling is then meshed in a plane to obtain the plane scaling mesh data;

[0045] The correction amount is determined based on the deviation angle of the restored profile line and the gridded data of the planar expansion and contraction.

[0046] Preferably, determining the correction amount based on the deviation angle of the restored profile line and the planar stretching mesh data includes:

[0047] Determine the cosine and sine trigonometric function values ​​of the angle between the restored profile line and the due north direction;

[0048] The new coordinates relative to the original coordinates are determined by multiplying the cosine trigonometric function value, the sine trigonometric function value, and the planar scaling grid data, respectively, in terms of the horizontal and vertical coordinate corrections.

[0049] Preferably, determining the paleogeographic map of the target area based on the first and second stratigraphic true thicknesses includes:

[0050] The true thicknesses of the first and second strata in the target area are three-dimensionally meshed to obtain paleogeographic surface data of the target area.

[0051] The paleogeographic surface data of the target area are processed to obtain a paleogeographic map of the target area.

[0052] This disclosure also proposes a paleogeomorphology restoration system based on balanced profile structural reconstruction, the system comprising:

[0053] A module is established to acquire drilling data of the target area and establish an isochronous stratigraphic framework based on the drilling data of the target area.

[0054] The first determining module is used to determine the target layer bottom interface and the maximum flooding surface of the target area based on the isochronous stratigraphic framework.

[0055] The second determining module is used to determine the first true thickness of the stratigraphic sedimentary layer in the target area based on the bottom interface of the target layer and the maximum flooding surface of the target area; the first true thickness of the stratigraphic sedimentary layer is the true thickness of the stratigraphic sedimentary layer in the target area under in-situ coordinates.

[0056] The third determining module is used to determine the correction amount of the new coordinates relative to the in-situ coordinates, and to determine the paleogeographic map of the target area based on the first and second true thicknesses of the sedimentary layers; the second true thickness of the sedimentary layers is the true thickness of the sedimentary layers of the target area under the new coordinates.

[0057] Preferably, the system further comprises:

[0058] The preprocessing module is used to preprocess the drilling data;

[0059] The preprocessing includes:

[0060] Perform single-well vertical consistency processing on the target area;

[0061] Alternatively, the target area can be subjected to multi-well plane standardization processing.

[0062] Preferably, the establishment module is used to establish an isochronous stratigraphic framework based on drilling data of the target area, including:

[0063] The establishment module is used to continuously select several target wells based on the sediment source direction of the target area;

[0064] The target well is divided into three-level sequence stratigraphy, and an isochronous stratigraphic framework is established; the isochronous stratigraphic framework includes the target stratigraphic framework, the isochronous stratigraphic framework, and the sequence stratigraphic framework.

[0065] Preferably, the establishment module is used to establish an isochronous stratigraphic framework based on drilling data of the target area, and further includes:

[0066] The establishment module is used to determine the rock compaction coefficient of the target area based on the relationship between rock porosity and depth in the target area;

[0067] And / or, determine the core ripple wavelength of the target layer in the target region.

[0068] Preferably, the first determining module is used to determine the target layer bottom interface and the maximum flooding level of the target area based on the isochronous stratigraphic framework, including:

[0069] The first determining module is used to calibrate the target layer and the maximum flooding surface of each target well in the isochronous stratigraphic framework onto the seismic profile reflection phase axis of the target area;

[0070] The type of the seismic profile of the target area is adjusted to a lithological profile, and the target layer bottom interface and the maximum flooding surface of the target area are determined based on the reflection characteristics of the target layer bottom interface and the maximum flooding surface on the reflection phase axis in the seismic profile.

[0071] Preferably, the second determining module is used to determine the first true thickness of the sedimentary layer in the target area based on the bottom interface of the target layer and the maximum flooding surface, including:

[0072] The second determining module is used to perform time-depth conversion on the time-domain apparent thickness of the target area to obtain the depth-domain apparent thickness of the target area;

[0073] Terrain correction is performed on the apparent depth thickness of the target region to obtain the true depth thickness of the target region.

[0074] The depth domain true thickness of the target area is corrected to determine the first stratum sedimentary true thickness of the target area.

[0075] Preferably, the third determining module is used to determine the correction amount of the new coordinates relative to the original coordinates, including:

[0076] The third determining module is used to construct the equilibrium evolution profile of the target region and determine the pre-recovery profile line on the equilibrium evolution profile.

[0077] The equilibrium evolution profile of the target region is restored, and the restored profile line is redefined.

[0078] The amount of scatter point expansion is determined by comparing the profile lines before and after restoration.

[0079] The scattered point scaling is then meshed in a plane to obtain the plane scaling mesh data;

[0080] The correction amount is determined based on the deviation angle of the restored profile line and the gridded data of the planar expansion and contraction.

[0081] Preferably, the third determining module is used to determine the paleogeographic map of the target area based on the first and second stratigraphic true thicknesses, including:

[0082] The third determining module is used to perform three-dimensional meshing on the true thickness of the first and second strata in the target area to obtain paleogeographic surface data of the target area.

[0083] The paleogeographic surface data of the target area are processed to obtain a paleogeographic map of the target area.

[0084] This disclosure also proposes a paleogeographic restoration device based on balanced profile structural restoration, comprising:

[0085] Processor and memory;

[0086] The processor invokes the computer program stored in the memory to execute the paleogeographic restoration method based on balanced profile structure restoration described above.

[0087] This disclosure also proposes a computer-readable storage medium.

[0088] The computer-readable storage medium stores a computer program that, when executed by a processor, enables the processor to perform the paleogeographic restoration method based on balanced profile structure restoration described above.

[0089] This disclosure has the following beneficial effects:

[0090] This disclosure establishes a rapid, reasonable, effective, and highly accurate method for restoring pre-depositional paleomorphological features. Based on geological background analysis, well logging petrophysical analysis, and seismic data interpretation, it utilizes the principle and method of balanced profile structural restoration to restore the thickness of the target layer to its true depositional thickness and return it to its original position at the time of deposition. This, in turn, restores the original paleomorphological features at the time of deposition, improves the rationality and accuracy of sedimentary facies studies and reservoir sand body studies, and ultimately increases the exploration success rate.

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

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

[0093] Figure 1 A diagram illustrating the paleogeographic restoration method based on balanced profile structure restoration in an embodiment of this disclosure is shown.

[0094] Figure 2 This document shows a detailed flowchart of the paleogeographic restoration method based on balanced profile structure restoration in an embodiment of this disclosure.

[0095] Figure 3 This diagram illustrates a paleogeographic reconstruction system based on balanced profile reconstruction in an embodiment of this disclosure.

[0096] Figure 4 A diagram of a paleogeographic restoration device based on balanced profile reconstruction is shown in an embodiment of this disclosure. Detailed Implementation

[0097] Example embodiments will now be described more fully with reference to the accompanying drawings. However, example embodiments can be implemented in many forms and should not be construed as limited to the examples set forth herein; rather, these embodiments are provided to make this disclosure more comprehensive and complete, and to fully convey the concept of the example embodiments to those skilled in the art. The described features, structures, or characteristics can be combined in any suitable manner in one or more embodiments. In the following description, numerous specific details are provided to give a full understanding of embodiments of this disclosure. However, those skilled in the art will recognize that the technical solutions of this disclosure can be practiced with one or more of the specific details omitted, or other methods, components, apparatus, steps, etc., can be employed. In other instances, well-known technical solutions are not shown or described in detail to avoid obscuring various aspects of this disclosure.

[0098] Furthermore, the accompanying drawings are merely illustrative of this disclosure and are not necessarily drawn to scale. The same reference numerals in the drawings denote the same or similar parts, and therefore repeated descriptions of them will be omitted. Some block diagrams shown in the drawings are functional entities and do not necessarily correspond to physically or logically independent entities. These functional entities may be implemented in software, in one or more hardware units or integrated circuits, or in different network and / or processor devices and / or microcontroller devices.

[0099] The flowchart shown in the attached diagram is merely an illustrative example and does not necessarily include all steps. For example, some steps may be broken down, while others may be combined or partially combined; therefore, the actual execution order may change depending on the specific circumstances.

[0100] The terms "first," "second," etc., used in the specification, claims, and accompanying drawings of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of this application described herein can be implemented, for example, in orders other than those illustrated or described herein.

[0101] Furthermore, the terms “comprising” and “having”, and any variations thereof, are intended to cover non-exclusive inclusion, such that a process, method, system, product, or device that includes a series of steps or sub-modules is not necessarily limited to those steps or sub-modules that are explicitly listed, but may include other steps or sub-modules that are not explicitly listed or that are inherent to such process, method, product, or device.

[0102] To accurately reconstruct paleomorphological features, study the relationship between paleomorphology and hydrocarbon accumulation, and rationally evaluate the development and distribution patterns of reservoir sand bodies, this paper proposes a paleomorphological restoration method and system based on balanced profile structural restoration, utilizing drilling and seismic data. This method was applied to reconstruct the pre-Cretaceous paleomorphological features of the Kuqa Depression in the Tarim Basin, providing a basis for sedimentary facies studies, reservoir in-depth research, and well location deployment for the Cretaceous Yageliemu Formation and Shushanhe Formation in this area.

[0103] like Figure 1 As shown, this disclosure proposes a paleogeographic restoration method based on balanced profile structural restoration, the method comprising:

[0104] S1 performs preprocessing on the drilling data;

[0105] S2 acquires drilling data for the target area and establishes an isochronous stratigraphic framework based on the drilling data for the target area;

[0106] S3 determines the target layer bottom interface and maximum flooding level of the target area based on the isochronous stratigraphic framework;

[0107] S4 determines the first true thickness of the stratigraphic sedimentary layer in the target area based on the bottom interface of the target layer and the maximum flooding surface; the first true thickness of the stratigraphic sedimentary layer is the true thickness of the stratigraphic sedimentary layer in the target area under in-situ coordinates.

[0108] S5 determines the correction amount of the new coordinates relative to the original coordinates, and determines the paleogeographic map of the target area based on the first and second true thicknesses of the sedimentary layers in the target area; the second true thickness of the sedimentary layers is the true thickness of the sedimentary layers in the target area under the new coordinates.

[0109] Specifically, the preprocessing in S1 includes:

[0110] Perform single-well vertical consistency processing on the target area;

[0111] Alternatively, the target area can be subjected to multi-well plane standardization processing.

[0112] Specifically, the isochronous stratigraphic framework established in S2 based on drilling data of the target area includes:

[0113] Several target wells are selected consecutively based on the sediment source direction of the target area;

[0114] The target well is divided into three-level sequence stratigraphy, and an isochronous stratigraphic framework is established; the isochronous stratigraphic framework includes the target stratigraphic framework, the isochronous stratigraphic framework, and the sequence stratigraphic framework.

[0115] Specifically, the isochronous stratigraphic framework established in S2 based on drilling data from the target area also includes:

[0116] Based on the relationship between rock porosity and depth in the target area, the rock compaction coefficient of the target area is determined;

[0117] And / or, determine the core ripple wavelength of the target layer in the target region.

[0118] Specifically, S3 determines the target layer bottom interface and maximum flooding level of the target area based on the isochronous stratigraphic framework, including:

[0119] The target layer and the maximum flooding surface of each target well in the isochronous stratigraphic framework are marked onto the seismic profile reflection phase axis of the target area;

[0120] The type of the seismic profile of the target area is adjusted to a lithological profile, and the target layer bottom interface and the maximum flooding surface of the target area are determined based on the reflection characteristics of the target layer bottom interface and the maximum flooding surface on the reflection phase axis in the seismic profile.

[0121] In this embodiment, seismic profiles are a common way to display seismic data in petroleum exploration. Seismic data has a set of data acquisition and processing procedures. Here, the processed data is displayed in two dimensions in Geoeast software. Lithological profiles are an interpretation and processing result of seismic profiles, which are obtained by 90-degree phase conversion of seismic profiles.

[0122] Specifically, S4 determines the true thickness of the first stratigraphic deposit in the target area based on the bottom interface of the target layer and the maximum flooding surface, including:

[0123] The temporal apparent thickness of the stratigraphy in the target area is determined based on the bottom interface of the target layer and the first maximum flooding surface.

[0124] The time-domain apparent thickness of the target area is converted to depth to obtain the depth-domain apparent thickness of the target area.

[0125] Terrain correction is performed on the apparent depth thickness of the target region to obtain the true depth thickness of the target region.

[0126] The depth domain true thickness of the target area is corrected to determine the first stratum sedimentary true thickness of the target area.

[0127] Specifically, determining the temporal apparent thickness of the stratigraphy in the target area based on the target layer bottom interface and the first maximum flooding surface includes:

[0128] The seismic interpretation interface of the target region at the target layer is extended to the erosion zone of the target region at the target layer;

[0129] The temporal apparent thickness of the target area is determined based on the difference between the bottom interface of the target layer and the first maximum flood level.

[0130] Specifically, the time-domain apparent thickness of the target region is converted to depth to obtain the depth-domain apparent thickness of the target region, including:

[0131] Construct the target layer velocity field of the target region;

[0132] The depth domain apparent thickness of the target region is determined based on the target layer velocity field and the time domain apparent thickness of the strata.

[0133] Specifically, terrain correction is performed on the apparent depth thickness of the target region to obtain the true depth thickness of the target region, including:

[0134] Extract the seismic dip angle attribute of the maximum flood surface in the target area;

[0135] Based on the earthquake dip angle attribute, the apparent depth thickness of the target area is dip-corrected to determine the true depth thickness of the target area.

[0136] Specifically, the depth domain true thickness of the target region is corrected, including:

[0137] Determine the residual correction amount, and perform residual correction on the true thickness of the depth domain based on the residual correction amount;

[0138] The compaction thickness is determined based on the rock compaction coefficient of the target area, and the true thickness of the depth domain is compacted and corrected based on the compaction thickness.

[0139] The anomalous body thickness correction amount is determined based on the deposition distance of the target area, and the anomalous body thickness is corrected for the true thickness of the depth domain based on the anomalous body thickness correction amount.

[0140] Specifically, determining the correction amount of the new coordinates relative to the original coordinates in S5 includes:

[0141] Construct an equilibrium evolution profile of the target region, and determine the pre-recovery profile line on the equilibrium evolution profile;

[0142] The equilibrium evolution profile of the target region is restored, and the restored profile line is redefined.

[0143] The amount of scatter point expansion is determined by comparing the profile lines before and after restoration.

[0144] The scattered point scaling is then meshed in a plane to obtain the plane scaling mesh data;

[0145] The correction amount is determined based on the deviation angle of the restored profile line and the gridded data of the planar expansion and contraction.

[0146] Specifically, determining the correction amount based on the deviation angle of the restored profile line and the planar stretching mesh data includes:

[0147] Determine the cosine and sine trigonometric function values ​​of the angle between the restored profile line and the due north direction;

[0148] The new coordinates relative to the original coordinates are determined by multiplying the cosine trigonometric function value, the sine trigonometric function value, and the planar scaling grid data, respectively, in terms of the horizontal and vertical coordinate corrections.

[0149] Specifically, the process in S5 of determining the paleogeographic map of the target area based on the first and second stratigraphic true thicknesses includes:

[0150] The true thicknesses of the first and second strata in the target area are three-dimensionally meshed to obtain paleogeographic surface data of the target area.

[0151] The paleogeographic surface data of the target area are processed to obtain a paleogeographic map of the target area.

[0152] The purpose of this disclosure is to establish a rapid, reasonable, effective, and highly accurate method for restoring pre-depositional paleomorphological features. Based on geological background analysis, well logging petrophysical analysis, and seismic data interpretation, this method utilizes the principle and method of balanced profile structural restoration to restore the thickness of the target layer to its true depositional thickness and return it to its original position at the time of deposition. This, in turn, restores the original paleomorphological features at the time of deposition, improves the rationality and accuracy of sedimentary facies studies and reservoir sand body studies, and ultimately increases the exploration success rate.

[0153] like Figure 2 As shown, this disclosure is implemented in the following manner:

[0154] 1) Drilling data analysis, processing, and isochronous formation framework establishment

[0155] The analysis, processing, and establishment of the isochronous formation framework for drilling data includes four parts: preprocessing of logging data, establishment of drilling framework comparison profiles, target processing of logging parameters in the profiles, and statistics of core ripple wavelengths.

[0156] Pre-processing of well logging data includes single-well vertical consistency processing and multi-well planar standardization processing. The single-well curves measured multiple times are processed for consistency by selecting the thick mudstone section above the target layer. The lithological sections with obvious and widely distributed logging curve characteristics are selected as marker layers to establish curve frequency distribution maps. The average value of the marker layer curves is used as the standard value to calibrate the curves of other wells and establish a unified lithological and logging comparison relationship.

[0157] Establish a drilling framework comparison profile, continuously select wells according to the direction of sediment source, and divide the target layer into three-level sequences according to the three-level sequence division method of sequence stratigraphy. Establish a framework of target layer, isochronous layer, and sequence stratigraphy, and determine the first maximum flooding surface (lake flooding surface) above the target layer.

[0158] Wavelength statistics of core ripple marks: the wavelength of core ripple marks in the target layer is statistically analyzed through drilling core observation.

[0159] 2) High-precision seismic horizon interpretation

[0160] The synthetic record calibration clearly identifies the bottom interface of the target layer and the seismic reflection phase axis of the first maximum flooding surface above it. High-pass filtering is applied to the seismic data to improve seismic resolution, and the seismic data is phase-modified by 90 degrees to match the seismic profile phase axis with the isochronous grid profile. Tracking parameters are set to automatically pick the target layer horizon. The process includes two parts: well-seismic joint calibration and automatic tracking and interpretation of the target layer and the maximum flooding surface horizon.

[0161] Well-seismic joint calibration involves calculating the reflection coefficient based on the well logging acoustic and density curves, extracting the time-shifted seismic trace wavelet, estimating the well composite seismic record by convolving the seismic trace and the reflection coefficient, matching and comparing the seismic record with the well-side seismic trace, and calibrating the target layer and the maximum flooding surface of each well in the isochronous stratigraphic framework onto the seismic profile reflection phase axis.

[0162] Stratification tracking interpretation involves pre-processing seismic data with high-pass filtering and 90-degree phase conversion to improve data accuracy and convert it into lithological profiles. This allows the seismic profiles to be compared with stratigraphic framework profiles, clearly identifying the seismic phase axis reflection characteristics of the target stratum's bottom interface and the maximum flooding surface. The stratigraphic tracking amplitude variation range is set for automatic stratigraphic tracking throughout the entire work area, completing the automatic interpretation of the seismic strata at the target stratum's bottom interface and the maximum flooding surface.

[0163] 3) Acquisition and correction of paleogeographic reconstruction data

[0164] The apparent thickness of the stratigraphy in the time domain is calculated by subtracting the bottom interface of the target layer from the first maximum flood level. The apparent thickness in the time domain is then converted to depth to obtain the apparent thickness in the depth domain. The topographically corrected apparent thickness is converted to the true thickness. The true thickness is then corrected to obtain the true thickness of the stratigraphy reconstructed from the in-situ paleogeography. The process includes seven parts: obtaining the apparent thickness of the stratigraphy, correcting the erosion amount, obtaining the true thickness, correcting the residuals, correcting the compaction, correcting the anomalies, and correcting the paleowater depth.

[0165] Apparent stratigraphic thickness is obtained by seismic horizon interpretation and tracing to obtain time-domain planar interpretation data of the bottom interface of the target layer and the maximum flooding surface. The time-domain apparent stratigraphic thickness is obtained by subtracting the maximum flooding surface data from the bottom interface horizon data.

[0166] For erosion correction, since erosion may exist in some areas, the trend surface extension method is used to extend the seismic interpretation interface of the target layer to the erosion area in the direction of the trend, construct the seismic interpretation layer of the target layer in the erosion area, restore the stratum thickness in the erosion area, and return to the previous step to recalculate and obtain the apparent thickness of the strata in the time domain after erosion correction.

[0167] True thickness is obtained by constructing the velocity field of the target layer based on the time-depth relationship established by well-seismic calibration, and calculating the apparent thickness in the depth domain from the layer velocity and the apparent thickness in the time domain. The dip angle attribute of the maximum flooding surface is extracted, and the true thickness of the formation is obtained by dip angle correction of the apparent thickness in the depth domain.

[0168] Residual correction involves statistically analyzing the formation thickness from the maximum flooding surface to the bottom boundary of the target layer and the estimated true formation thickness at the depth domain at the previous well point. The difference between these values ​​is used to calculate the residual correction amount for the estimated well point thickness. The residual correction amount is then gridded to obtain gridded data. Finally, the gridded residual correction amount is summed with the true formation thickness from the previous step to obtain the true formation thickness after residual correction.

[0169] Compaction correction is performed. The rock compaction coefficient is obtained by fitting the relationship between rock porosity and depth using cross plot analysis. The compacted thickness is calculated by combining the true thickness of the formation after residual correction with the compaction coefficient. The true thickness of the formation after compaction correction is obtained by combining the compacted thickness with the true thickness of the formation after residual correction calculated in the previous step.

[0170] Anomaly correction is performed based on sedimentary principles: deeper paleocurrents, lower elevations, and higher clay content result in larger logging gamma values. Waveform indicator phase-controlled gamma ray inversion is conducted using a selected lithology-sensitive gamma curve to obtain the inverted gamma data volume. The root-mean-square (RMS) gamma of the target layer is extracted, and the average gamma is calculated. Correction coefficients are derived from the RMS and average gamma values. These coefficients are used to estimate the thickness correction for the lithological anomaly. The corrected thickness is then calculated by combining this correction with the compacted true formation thickness calculated in the previous step to obtain the corrected true formation thickness.

[0171] Paleowater depth correction: Based on the ripple mark wavelengths obtained from the above statistics, the paleowater depth of each well point is restored using the Diem paleowater depth restoration formula. The paleowater depth is then gridded to obtain gridded paleowater depth data. The true thickness of the formation after correction is calculated using the paleowater depth gridded data and the anomaly body calculated in the previous step to obtain the true thickness data after paleowater depth correction.

[0172] 4) Balanced profile restoration

[0173] Multiple equilibrium evolution profiles were constructed within the study area. Equilibrium restoration techniques were used to remove faults and folds from the profiles, restoring the stretching and contraction of the target layer during its depositional period. Base points were established along the profile lines, and the stretching and contraction at these base points were estimated. The stretching and contraction at each point along the profile line were then interpolated using these base points. The process comprises four parts: profile line selection, equilibrium profile restoration, base point stretching and contraction estimation, and profile line stretching and contraction interpolation.

[0174] For profile selection, choose a profile that runs through the entire work area, either along the direction of maximum paleostress or perpendicular to the main structural trend.

[0175] Balanced profile restoration involves constructing stable zones (where faults and folds are not developed) along the profile line, i.e., areas with stable and flat terrain, and establishing nail points as the zero point for structural balance restoration. The balance restoration technique is used to remove faults and folds from the profile. Based on mechanical properties, the fault is pushed back and the strata on both sides of the fault are connected to remove the fault. After removing the fault, the target layer interface is straightened and restored to a horizontal position to remove folds.

[0176] For estimation of the expansion and contraction of the base points, multiple base points are constructed on the profile line before restoration, and the expansion and contraction of each base point are obtained by comparing with the profile line after restoration.

[0177] The profile line expansion / contraction interpolation involves linearly interpolating the expansion / contraction amounts of each scattered point (pin point, foundation point) along the profile line to obtain an array of scattered expansion / contraction amounts for the profile line.

[0178] 5) Coordinate system conversion for paleogeographic mapping

[0179] The process involves three parts: planar scaling of the profile lines, mapping of the scatter points to obtain the planar scaling, calculation of the angle between the profile lines and true north, and the values ​​of the cosine and sine trigonometric functions of the angle. The planar scaling is then corrected using cosine and sine to obtain the corrected vertical and horizontal coordinates. Finally, the true thickness of the sedimentary layers is corrected for paleowater depth under the new coordinate system, and paleogeographic mapping is performed.

[0180] The scaling factor is meshed in a plane. Given the work area boundary and the grid density, the scaling factor is processed by the minimum energy grid method of the DoubleFox software to obtain the scaling factor plane mesh data.

[0181] The new coordinate system is corrected by measuring the angle between the profile line and the due north direction based on the direction of the profile line on the map, calculating the cosine and sine trigonometric function values ​​of the angle, multiplying the gridded data of the scaling amount with the cosine and sine trigonometric function values ​​to obtain the gridded corrected longitudinal and abscissa coordinates, and adding the original coordinates with the corrected longitudinal and abscissa coordinates to obtain the new corrected coordinates.

[0182] Paleomorphological mapping was achieved by combining the true thickness data of sedimentary strata under the new coordinate system and the in-situ coordinate system into a three-dimensional gridded data. The data was then drawn using Shuanghu contour line drawing software to obtain the paleomorphological surface data and to obtain a three-dimensional paleomorphological map.

[0183] Example

[0184] 1) Drilling data analysis, processing, and isochronous formation framework establishment

[0185] The analysis, processing, and establishment of the isochronous formation framework for drilling data includes four parts: preprocessing of logging data, establishment of drilling framework comparison profiles, target processing of logging parameters in the profiles, and statistics of core ripple wavelengths.

[0186] Pre-processing of well logging data includes single-well longitudinal consistency processing and multi-well planar standardization processing. The single-well curves measured multiple times are processed for consistency by selecting the thick mudstone section above the target layer. The lithological sections with obvious and widely distributed logging curve characteristics are selected as marker layers to establish curve frequency distribution maps. The average value of the marker layer curves is used as the standard value to calibrate the curves of other wells and establish a unified lithological and logging comparison relationship.

[0187] Establish a drilling framework correlation profile, continuously select wells based on sediment source direction, and divide the target layer into three-level sequences according to the three-level sequence stratigraphy principle. Establish an isochronous sequence stratigraphic framework for the target layer and determine the first maximum flooding surface (lacustrine flooding surface) above the target layer.

[0188] Target processing of well logging parameters was performed, and cross-plot analysis was conducted to determine the relationship between rock porosity and depth. The rock compaction coefficient was fitted based on Guo Qiulin's classical exponential relationship formula, with the fitted formula being: φ=φ0e(-Ch), where h is the formation burial depth in meters, φ is the rock porosity at depth h in percentage, φ0 is the initial porosity at the surface, and C is the compaction coefficient. Cross-plot analysis of well logging curves and drilling lithology was conducted to select the lithology-sensitive well logging curve. Through analysis, the gamma curve was selected as the lithology-sensitive curve for this study.

[0189] Wavelength statistics of core ripple marks: The wavelength of core ripple marks in the target layer is statistically analyzed by observing and recording core samples from drilling. It is denoted as λ and the unit is m.

[0190] 2) High-precision seismic horizon interpretation

[0191] The synthetic record calibration clearly identifies the bottom interface of the target layer and the seismic reflection phase axis of the first maximum flooding surface above it. High-pass filtering is applied to the seismic data to improve seismic resolution, and the seismic data is phase-modified by 90 degrees to match the seismic profile phase axis with the isochronous grid profile. Tracking parameters are set to automatically pick the target layer horizon. The process includes two parts: well-seismic joint calibration and automatic tracking and interpretation of the target layer and the maximum flooding surface horizon.

[0192] Well-seismic joint calibration involves calculating the reflection coefficient based on the well logging acoustic and density curves, extracting the time-shifted seismic trace wavelet, estimating the well composite seismic record by convolving the seismic trace and the reflection coefficient, matching and comparing the seismic record with the well-side seismic trace, and calibrating the target layer and the maximum flooding surface of each well in the isochronous stratigraphic framework onto the seismic profile reflection phase axis.

[0193] Stratification tracking interpretation involves pre-processing seismic data with high-pass filtering and 90-degree phase conversion to improve data accuracy and convert it into lithological profiles. This allows the seismic profiles to be compared with stratigraphic framework profiles, clearly identifying the seismic phase axis reflection characteristics of the target stratum's bottom interface and the maximum flooding surface. The stratigraphic tracking amplitude variation range is set for automatic stratigraphic tracking throughout the entire work area, completing the automatic interpretation of the seismic strata at the target stratum's bottom interface and the maximum flooding surface.

[0194] 3) Acquisition and correction of paleogeographic reconstruction data

[0195] The apparent thickness of the stratigraphy in the time domain is calculated by subtracting the bottom interface of the target layer from the first maximum flood level. The apparent thickness in the time domain is then converted to depth to obtain the apparent thickness in the depth domain. The topographically corrected apparent thickness is converted to the true thickness. The true thickness is then corrected to obtain the true thickness of the stratigraphy reconstructed from the in-situ paleogeography. The process includes seven parts: obtaining the apparent thickness of the stratigraphy, correcting the erosion amount, obtaining the true thickness, correcting the residuals, correcting the compaction, correcting the anomalies, and correcting the paleowater depth.

[0196] The apparent thickness of the strata is obtained by seismic horizon interpretation and tracing. The time-domain planar interpretation data of the bottom interface (TM, ms) and the maximum flooding surface (Tf, ms) of the target layer are obtained. The bottom interface horizon data is subtracted from the maximum flooding surface horizon data, HT = TM - Tf, to obtain the time-domain apparent thickness of the strata (HT, ms).

[0197] For erosion correction, since erosion may exist in some areas, the trend surface extension method is used to extend the seismic interpretation interface of the target layer to the erosion area in the direction of the trend, construct the seismic interpretation layer of the target layer in the erosion area, restore the stratum thickness in the erosion area, and return to the previous step to recalculate and obtain the apparent thickness of the strata in the time domain after erosion correction.

[0198] True thickness is obtained by constructing a velocity field based on the time-depth relationship established by well-seismic calibration. The target layer velocity (V, m / s) is multiplied by the apparent thickness in the time domain (HT) to obtain the apparent thickness in the depth domain, HD = HT * V. The maximum flooding surface seismic dip angle attribute (α) is extracted, and the apparent thickness in the depth domain is corrected for dip angle to obtain the true thickness of the formation (HR, m), HR = HD * cos(α).

[0199] Residual correction involves statistically analyzing the formation thickness (Hw, m) from the maximum flooding surface to the bottom boundary of the target layer and the estimated true formation thickness (HR) at the previous well point. The difference is then used to calculate the residual correction amount (HΔ, m) for the well point thickness, where HΔ = Hw - HR. The residual correction amount is then gridded to obtain gridded data. Finally, the residual correction amount of the gridded data is summed with the true formation thickness from the previous step to obtain the true formation thickness after residual correction (HRS, m).

[0200] Compaction correction: The rock compaction coefficient (C) is obtained by fitting the relationship between rock porosity and depth using cross plot analysis. The compacted thickness (HC, m) is calculated by combining the residual corrected true thickness of the formation with the compaction coefficient, where HC = HRS * C. The compacted thickness is then combined with the residual corrected true thickness of the formation calculated in the previous step to obtain the compacted corrected true thickness of the formation (HY, m), where HY = HC + HRS.

[0201] Anomaly correction is performed based on sedimentary principles: deeper paleocurrents, lower elevations, and higher clay content result in larger logging gamma values. Waveform indicator phase control simulation gamma inversion is then performed using the selected lithology-sensitive gamma curves to obtain the inverted gamma data volume. The root mean square gamma of the target layer is extracted, and the average gamma (API) is calculated. The correction factor (A) is obtained from the root mean square gamma (GrRMS) and the average gamma (Grave), where A = GrRMS(i,j) / Grave, and GrRMS(i,j) represents the root mean square gamma value at point (i,j). The thickness correction amount of the lithological anomaly is estimated from the correction factor, where Ha = (HY_ave - HY(i,j)) * A, HY_ave represents the average thickness of the formation after compaction correction, and HY(i,j) represents the thickness of the formation after compaction correction at point (i,j). The thickness correction amount of the anomaly is combined with the true thickness of the formation after compaction correction calculated in the previous step to obtain the true thickness of the formation after anomaly correction (HG,m), where HG = HY(i,j) + Ha.

[0202] Paleodepth correction is performed based on the wave mark wavelengths (λ, m) obtained from the above statistics, using the Diem paleodepth reconstruction formula:

[0203]

[0204] The paleowater depths at each well point are recovered, and the paleowater depth gridded data is obtained by gridding. The true thickness of the formation after correction of the anomaly body calculated in the previous step is then calculated to obtain the true thickness data (Hend, m) after paleowater depth correction, where Hend = HG + h. max .

[0205] Where: h max To restore the ancient water depth, m; λ is the orbital velocity of water particles, m / s; λ is the wavelength of the ripple mark, m; D is the diameter of the sediment particles, m; ρ is the density of the water, g / cm³. 3 ;ρ s Sediment density (g / cm³) 3 .

[0206] 4) Balanced profile restoration

[0207] Multiple equilibrium evolution profiles were constructed within the study area. Equilibrium restoration techniques were used to remove faults and folds from the profiles, restoring the stretching and contraction of the target layer during its depositional period. Base points were established along the profile lines, and the stretching and contraction at these base points were estimated. The stretching and contraction at each point along the profile line were then interpolated using these base points. The process comprises four parts: profile line selection, equilibrium profile restoration, base point stretching and contraction estimation, and profile line stretching and contraction interpolation.

[0208] Section line selection: Select a section line that runs through the entire work area, either along the direction of maximum paleostress or perpendicular to the main structural trend.

[0209] Balanced profile restoration involves constructing stable zones (where faults and folds are not developed) along the profile line, i.e., areas of stable terrain, and establishing nail points as the zero point for structural balance restoration. The balance restoration technique is used to remove faults and folds from the profile. Based on mechanical properties, the fault is pushed back to connect the strata on both sides of the fault to remove the fault. After removing the fault, the target layer interface is straightened and restored to a horizontal position to remove folds.

[0210] For estimation of the expansion and contraction of the base points, multiple base points are constructed on the profile line before restoration, and the expansion and contraction of each base point (L, m) is obtained by comparing it with the profile line after restoration.

[0211] The profile line expansion and contraction amount interpolation is performed by linearly interpolating the expansion and contraction amounts of each scattered point (pin point, foundation point) along the profile line to obtain the profile line scattered point expansion and contraction amount array L(n).

[0212] 5) Coordinate system conversion for paleogeographic mapping

[0213] The process involves three parts: planar scaling of the profile lines, mapping of the scatter points to obtain the planar scaling, calculation of the angle between the profile lines and true north, and the values ​​of the cosine and sine trigonometric functions of the angle. The planar scaling is then corrected using cosine and sine to obtain the corrected vertical and horizontal coordinates. Finally, the true thickness of the sedimentary layers is corrected for paleowater depth under the new coordinate system, and paleogeographic mapping is performed.

[0214] The scaling factor is planarized. Given the work area boundary and the grid density, the scaling factor scattered points are processed by the minimum energy grid method of the Shuanghu software to obtain the scaling factor planarized grid data L(i,j).

[0215] The new coordinate system calibration involves measuring the angle (β) between the profile line and true north on the map, calculating the cosine (cos(β)) and sine (sin(β)) trigonometric function values ​​of the angle, and multiplying the gridded data of the scaling factor with the cosine and sine trigonometric function values ​​to obtain the gridded calibrated vertical and horizontal coordinates.

[0216] XΔ(i,j)=L(i,j)*sin(β), YΔ(i,j)=L(i,j)*cos(β),

[0217] The original coordinates of the calculated point are added to the corrected ordinates to obtain the new corrected coordinate system.

[0218] XN(i,j)=XΔ(i,j)+XO(i,j), YN(i,j)=YΔ(i,j)+YO(i,j).

[0219] Paleomorphological mapping was achieved by combining the true thickness data of sedimentary strata under the new coordinate system and the in-situ coordinate system to form a three-dimensional gridded data (X, Y, Hend). The data was then drawn using Shuanghu contour mapping software to obtain the paleomorphological surface data and to obtain a three-dimensional paleomorphological map.

[0220] like Figure 3 As shown, this disclosure also proposes a paleogeomorphological restoration system based on balanced profile structural restoration, the system comprising:

[0221] Preprocessing module 10 is used to preprocess the drilling data;

[0222] Module 20 is used to acquire drilling data of the target area and establish an isochronous stratigraphic framework based on the drilling data of the target area.

[0223] The first determining module 30 is used to determine the target layer bottom interface and the maximum flooding surface of the target area based on the isochronous stratigraphic framework.

[0224] The second determining module 40 is used to determine the first true thickness of the stratigraphic sedimentary layer in the target area based on the bottom interface of the target layer and the maximum flooding surface of the target area; the first true thickness of the stratigraphic sedimentary layer is the true thickness of the stratigraphic sedimentary layer in the target area under in-situ coordinates.

[0225] The third determining module 50 is used to determine the correction amount of the new coordinates relative to the original coordinates, and to determine the paleogeographic map of the target area based on the first and second true thicknesses of the sedimentary layers in the target area; the second true thickness of the sedimentary layers is the true thickness of the sedimentary layers in the target area under the new coordinates.

[0226] like Figure 4 As shown, corresponding to the paleogeographic restoration method based on balanced profile structure restoration provided above, this disclosure also provides a paleogeographic restoration device based on balanced profile structure restoration. Since the embodiment of this device is similar to the embodiment of the method described above, the description is relatively simple. For relevant details, please refer to the description in the embodiment section above. The device described below is merely illustrative. The device may include: a processor 1, a memory 2, a communication bus (i.e., the aforementioned device bus), and a lookup engine. The processor 1 and the memory 2 communicate with each other through the communication bus and communicate with external systems through a communication interface. The processor 1 can call logical instructions in the memory 2 to execute the paleogeographic restoration method based on balanced profile structure restoration.

[0227] Furthermore, the logical instructions in the aforementioned memory 2 can be implemented as software functional units and sold or used as independent products, and can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this disclosure, in essence, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this disclosure. The aforementioned storage medium includes various media capable of storing program code, such as memory chips, USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0228] On the other hand, this disclosure also provides a processor-readable storage medium storing a computer program 3, which, when executed by a processor 1, is implemented to perform the paleogeographic restoration method based on balanced profile structure restoration provided in the above embodiments.

[0229] The processor-readable storage medium can be any available medium or data storage device that the processor 1 can access, including but not limited to magnetic memory (e.g., floppy disk, hard disk, magnetic tape, magneto-optical disk (MO)), optical memory (e.g., CD, DVD, BD, HVD), and semiconductor memory (e.g., ROM, EPROM, EEPROM, non-volatile memory (NAND FLASH), solid-state drive (SSD)).

[0230] Those skilled in the art should understand that, despite the detailed description of this disclosure with reference to the foregoing embodiments, modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this disclosure.

Claims

1. A paleogeographic restoration method based on balanced profile structural restoration, characterized in that, The method comprises: Obtain drilling data for the target area and establish an isochronous stratigraphic framework based on the drilling data for the target area; The target stratum bottom interface and maximum flood level of the target area are determined based on the isochronous stratigraphic framework. The first true thickness of the sedimentary layer in the target area is determined based on the bottom interface of the target layer and the maximum flooding surface; the first true thickness of the sedimentary layer is the true thickness of the sedimentary layer in the target area under in-situ coordinates. The correction amount of the new coordinates relative to the original coordinates is determined, and the paleogeographic map of the target area is determined based on the first and second true thicknesses of the sedimentary layers in the target area; the second true thickness of the sedimentary layers is the true thickness of the sedimentary layers in the target area under the new coordinates.

2. The paleogeographic restoration method based on balanced profile structure restoration according to claim 1, characterized in that, Before establishing an isochronous stratigraphic framework based on drilling data from the target area, the process includes: preprocessing the drilling data; The preprocessing includes: Perform single-well vertical consistency processing on the target area; Alternatively, the target area can be subjected to multi-well plane standardization processing.

3. The paleogeographic restoration method based on balanced profile structure restoration according to claim 1, characterized in that, Establishing an isochronous stratigraphic framework based on drilling data from the target area includes: Several target wells are selected consecutively based on the sediment source direction of the target area; The target well is divided into three-level sequence stratigraphy, and an isochronous stratigraphic framework is established; the isochronous stratigraphic framework includes the target stratigraphic framework, the isochronous stratigraphic framework, and the sequence stratigraphic framework.

4. The paleogeographic restoration method based on balanced profile structure restoration according to claim 1, characterized in that, Establishing an isochronous stratigraphic framework based on drilling data from the target area also includes: Based on the relationship between rock porosity and depth in the target area, the rock compaction coefficient of the target area is determined; And / or, determine the core ripple wavelength of the target layer in the target region.

5. The paleogeographic restoration method based on balanced profile structure restoration according to claim 1, characterized in that, Determining the target layer bottom interface and maximum flooding level of the target area based on the isochronous stratigraphic framework includes: The target layer and the maximum flooding surface of each target well in the isochronous stratigraphic framework are marked onto the seismic profile reflection phase axis of the target area; The type of the seismic profile of the target area is adjusted to a lithological profile, and the target layer bottom interface and the maximum flooding surface of the target area are determined based on the reflection characteristics of the target layer bottom interface and the maximum flooding surface on the reflection phase axis in the seismic profile.

6. The paleogeographic restoration method based on balanced profile structure restoration according to claim 1, characterized in that, The true thickness of the first stratigraphic layer in the target area is determined based on the bottom interface of the target layer and the maximum flooding surface, including: The temporal apparent thickness of the stratigraphy in the target area is determined based on the bottom interface of the target layer and the first maximum flooding surface. The time-domain apparent thickness of the target area is converted to depth to obtain the depth-domain apparent thickness of the target area. Terrain correction is performed on the apparent depth thickness of the target region to obtain the true depth thickness of the target region. The depth domain true thickness of the target area is corrected to determine the first stratum sedimentary true thickness of the target area.

7. The paleogeographic restoration method based on balanced profile structure restoration according to claim 6, characterized in that, Determining the temporal apparent thickness of the stratigraphy of the target region based on the bottom interface of the target layer and the first maximum flooding surface includes: The seismic interpretation interface of the target region at the target layer is extended to the erosion zone of the target region at the target layer; The temporal apparent thickness of the target area is determined based on the difference between the bottom interface of the target layer and the first maximum flood level.

8. The paleogeographic restoration method based on balanced profile structure restoration according to claim 6, characterized in that, The time-domain apparent thickness of the target region is converted to depth to obtain the depth-domain apparent thickness of the target region, including: Construct the target layer velocity field of the target region; The depth domain apparent thickness of the target region is determined based on the target layer velocity field and the time domain apparent thickness of the strata.

9. The paleogeographic restoration method based on balanced profile structure restoration according to claim 6, characterized in that, To obtain the true depth domain thickness of the target region, terrain correction is performed on the apparent depth domain thickness of the target region, including: Extract the seismic dip angle attribute of the maximum flood surface in the target area; Based on the earthquake dip angle attribute, the apparent depth thickness of the target area is dip-corrected to determine the true depth thickness of the target area.

10. The paleogeographic restoration method based on balanced profile structure restoration according to claim 6, characterized in that, Correcting the true depth of the target region in the depth domain includes: Determine the residual correction amount, and perform residual correction on the true thickness of the depth domain based on the residual correction amount; The compaction thickness is determined based on the rock compaction coefficient of the target area, and the true thickness of the depth domain is compacted and corrected based on the compaction thickness. The anomalous body thickness correction amount is determined based on the deposition distance of the target area, and the anomalous body thickness is corrected for the true thickness of the depth domain based on the anomalous body thickness correction amount.

11. The paleogeographic restoration method based on balanced profile structure restoration according to claim 1, characterized in that, The determination of the correction amount for the new coordinates relative to the original coordinates includes: Construct an equilibrium evolution profile of the target region, and determine the pre-recovery profile line on the equilibrium evolution profile; The equilibrium evolution profile of the target region is restored, and the restored profile line is redefined. The amount of scatter point expansion is determined by comparing the profile lines before and after restoration. The scattered point scaling is then meshed in a plane to obtain the plane scaling mesh data; The correction amount is determined based on the deviation angle of the restored profile line and the gridded data of the planar expansion and contraction.

12. The paleogeographic restoration method based on balanced profile structure restoration according to claim 11, characterized in that, The step of determining the correction amount based on the deviation angle of the restored profile line and the planar stretching mesh data includes: Determine the cosine and sine trigonometric function values ​​of the angle between the restored profile line and the due north direction; The new coordinates relative to the original coordinates are determined by multiplying the cosine trigonometric function value, the sine trigonometric function value, and the planar scaling grid data, respectively, in terms of the horizontal and vertical coordinate corrections.

13. The paleogeographic restoration method based on balanced profile structure restoration according to claim 1, characterized in that, The process of determining the paleogeographic map of the target area based on the true thickness of the first and second stratigraphic layers includes: The true thicknesses of the first and second strata in the target area are three-dimensionally meshed to obtain paleogeographic surface data of the target area. The paleogeographic surface data of the target area are processed to obtain a paleogeographic map of the target area.

14. A paleogeographic restoration system based on balanced profile structural reconstruction, characterized in that, The system includes: A module is established to acquire drilling data of the target area and establish an isochronous stratigraphic framework based on the drilling data of the target area. The first determining module is used to determine the target layer bottom interface and the maximum flooding surface of the target area based on the isochronous stratigraphic framework. The second determining module is used to determine the first true thickness of the stratigraphic sedimentary layer in the target area based on the bottom interface of the target layer and the maximum flooding surface of the target area; the first true thickness of the stratigraphic sedimentary layer is the true thickness of the stratigraphic sedimentary layer in the target area under in-situ coordinates. The third determining module is used to determine the correction amount of the new coordinates relative to the original coordinates, and to determine the paleogeographic map of the target area based on the first and second true thicknesses of the sedimentary layers in the target area; the second true thickness of the sedimentary layers is the true thickness of the sedimentary layers in the target area under the new coordinates.

15. The paleogeographic restoration system based on balanced profile structure restoration according to claim 14, characterized in that, The system further comprises: The preprocessing module is used to preprocess the drilling data; The preprocessing includes: Perform single-well vertical consistency processing on the target area; Alternatively, the target area can be subjected to multi-well plane standardization processing.

16. The paleogeographic restoration system based on balanced profile structure restoration according to claim 14, characterized in that, The establishment module is used to establish an isochronous stratigraphic framework based on drilling data of the target area, including: The establishment module is used to continuously select several target wells based on the sediment source direction of the target area; The target well is divided into three-level sequence stratigraphy, and an isochronous stratigraphic framework is established; the isochronous stratigraphic framework includes the target stratigraphic framework, the isochronous stratigraphic framework, and the sequence stratigraphic framework.

17. The paleogeographic restoration system based on balanced profile structure restoration according to claim 14, characterized in that, The establishment module is used to establish an isochronous stratigraphic framework based on drilling data of the target area, and also includes: The establishment module is used to determine the rock compaction coefficient of the target area based on the relationship between rock porosity and depth in the target area; And / or, determine the core ripple wavelength of the target layer in the target region.

18. The paleogeographic restoration system based on balanced profile structure restoration according to claim 14, characterized in that, The first determining module is used to determine the target layer bottom interface and the maximum flood level of the target area based on the isochronous stratigraphic framework, including: The first determining module is used to calibrate the target layer and the maximum flooding surface of each target well in the isochronous stratigraphic framework onto the seismic profile reflection phase axis of the target area; The type of the seismic profile of the target area is adjusted to a lithological profile, and the target layer bottom interface and the maximum flooding surface of the target area are determined based on the reflection characteristics of the target layer bottom interface and the maximum flooding surface on the reflection phase axis in the seismic profile.

19. The paleogeographic restoration system based on balanced profile structure restoration according to claim 14, characterized in that, The second determining module is used to determine the true thickness of the first stratigraphic sedimentary layer in the target area based on the bottom interface of the target layer and the maximum flooding surface, including: The second determining module is used to perform time-depth conversion on the time-domain apparent thickness of the target area to obtain the depth-domain apparent thickness of the target area; Terrain correction is performed on the apparent depth thickness of the target region to obtain the true depth thickness of the target region. The depth domain true thickness of the target area is corrected to determine the first stratum sedimentary true thickness of the target area.

20. The paleogeographic restoration system based on balanced profile structure restoration according to claim 14, characterized in that, The third determining module is used to determine the correction amount of the new coordinates relative to the original coordinates, including: The third determining module is used to construct the equilibrium evolution profile of the target region and determine the pre-recovery profile line on the equilibrium evolution profile. The equilibrium evolution profile of the target region is restored, and the restored profile line is redefined. The amount of scatter point expansion is determined by comparing the profile lines before and after restoration. The scattered point scaling is then meshed in a plane to obtain the plane scaling mesh data; The correction amount is determined based on the deviation angle of the restored profile line and the gridded data of the planar expansion and contraction.

21. The paleogeographic restoration system based on balanced profile structure restoration according to claim 14, characterized in that, The third determining module is used to determine the paleogeographic map of the target area based on the true thickness of the first and second stratigraphic layers, including: The third determining module is used to perform three-dimensional meshing on the true thickness of the first and second strata in the target area to obtain paleogeographic surface data of the target area. The paleogeographic surface data of the target area are processed to obtain a paleogeographic map of the target area.

22. A paleogeographic restoration device based on balanced profile structural restoration, characterized in that, include: Processor and memory; The processor invokes the computer program stored in the memory to execute the paleogeographic restoration method based on balanced profile structure restoration as described in any one of claims 1 to 13.

23. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, enables the processor to perform the paleogeographic restoration method based on balanced profile structure restoration as described in any one of claims 1 to 13.

Images