Method and system for restoring palaeogeomorphology in complex fault-block region based on structural backstripping

Abstract

The present application belongs to the technical field of oil and gas exploration and development, and relates to a complex fault block area paleogeomorphology recovery method and system based on structure backstripping, the method steps being: obtaining a structure map of a target layer to be recovered, identifying faults in the structure map, and sequentially numbering the faults according to a set direction to obtain a numbered structure map; dividing the vertical fault throw of each fault according to the calculated vertical fault throw of each fault in the numbered structure map, and dividing the horizontal fault throw of each fault according to the measured horizontal fault throw of each fault in the numbered structure map; setting a reference surface, respectively calculating the distance TR from the top surface structure map after fault throw recovery to the reference surface and the distance BR from the bottom surface structure map after fault throw recovery to the reference surface, and obtaining the compacted paleogeomorphology surface by subtracting the distance TR from the distance BR. The present application quantitatively recovers the paleogeomorphology of the target layer during the sedimentary period, and can effectively predict the sedimentary system and the hydrocarbon generation center.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of oil and gas exploration and development, and relates to a palaeogeomorphology restoration technology, in particular to a palaeogeomorphology restoration method and system based on structural backstripping in a complex fault block area. BACKGROUND

[0002] In the analysis of sedimentary basins, the palaeogeomorphology restoration technology plays an important role. It provides key support for oil and gas exploration and development through the research on the geomorphology environment and its changes in the geological history, and has important guiding significance in terms of the source direction, spatial changes of sedimentary area and sedimentary center, sedimentary facies distribution, hydrocarbon source rock development, reservoir distribution and oil and gas migration and accumulation.

[0003] The existing palaeogeomorphology restoration methods mainly include the following categories: methods based on stratigraphic thickness, such as the impression method and residual thickness method, which infer the palaeogeomorphology features by analyzing the thickness changes of strata; sequence stratigraphy methods, such as the layer flattening method, which restores the geomorphology through the geometric shape and relationship of stratigraphic sequences; sedimentology methods, which reconstruct the palaeogeomorphology through the transformation characteristics of sedimentary facies; basin simulation methods, which simulate and restore the geomorphology based on three-dimensional geological models; and seismic geomorphology methods, which depict the palaeogeomorphology through seismic section technology.

[0004] Although the above methods have certain application value in palaeogeomorphology restoration, they still have significant limitations in complex fault block areas. In complex fault block areas, due to the development of faults, the target layer is divided into multiple different fault blocks in the plane, and there are differences in horizontal and vertical fault throw between different fault blocks. However, the current restoration methods are difficult to effectively solve the fault throw problem caused by fault activity when dealing with complex fault block areas, thereby resulting in low palaeogeomorphology restoration accuracy. This accuracy limitation directly affects the reservoir development and the effectiveness of oil and gas accumulation in complex fault block areas, and thus reduces the practical application value in the field of oil and gas exploration. SUMMARY

[0005] In view of the above problems such as low palaeogeomorphology restoration accuracy in the prior art, the application provides a palaeogeomorphology restoration method and system based on structural backstripping in a complex fault block area, which can improve the palaeogeomorphology restoration accuracy in a complex fault block area and provide a scientific basis for reservoir development prediction and oil and gas accumulation research.

[0006] In order to achieve the above purpose, in a first aspect, the application provides a palaeogeomorphology restoration method based on structural backstripping in a complex fault block area, and the steps are as follows:

[0007] A data acquisition step: acquiring a structure map of a target layer to be restored, wherein the structure map includes a top surface structure map and a bottom surface structure map;

[0008] Fault numbering step: identify faults in the structure map, and sequentially number the faults in a set direction to obtain a numbered structure map;

[0009] Fault distance recovery step: according to the calculated vertical fault distance of each fault in the numbered structure map, recover the vertical fault distance of each fault by segmenting the block, and according to the measured horizontal fault distance of each fault in the numbered structure map, recover the horizontal fault distance of each fault by segmenting the block;

[0010] Paleogeomorphology recovery step: set a reference surface, calculate the distance TR from the top surface structure map after fault distance recovery to the reference surface, and the distance BR from the bottom surface structure map after fault distance recovery to the reference surface, and obtain the paleogeomorphology shape surface after compaction by subtracting the distance TR from the distance BR.

[0011] In some embodiments, it further includes an outlier rejection step: rejecting outliers in the paleogeomorphology shape surface after compaction, the outliers being points that do not obviously conform to the trend.

[0012] In some embodiments, it further includes a compaction correction step: calculating the compression coefficient of the target layer according to the porosity of the stratum at the time of deposition and the porosity of the target layer, and multiplying the calculated compression coefficient by the paleogeomorphology shape surface to obtain the paleogeomorphology after compaction correction.

[0013] In some embodiments, in the fault numbering step, the method of identifying faults in the structure map is:

[0014] Graph construction step: constructing a fault polygon plane graph, and constructing a top surface polygon layer and a bottom surface polygon layer in the plane graph;

[0015] Copy step: copying the fault polygon in the top surface structure map to the top surface polygon layer, and copying the polygon in the bottom surface structure map to the bottom surface polygon layer;

[0016] Identification step: in the top surface fault polygon and the bottom surface fault polygon, fault polygons with similar lengths and positions are considered as the same fault.

[0017] In some embodiments, in the fault numbering step, the method of sequentially numbering the faults in a set direction to obtain a numbered structure map is:

[0018] Layer construction step: constructing a fault numbering layer in the plane graph;

[0019] Numbering step: sequentially numbering the identified faults in a set direction, and placing them in the fault numbering layer;

[0020] Copy step: copying the fault numbering layer after numbering to the top surface structure map and the bottom surface structure map to obtain a numbered structure map.

[0021] In some embodiments, the step of restoring the fault throw comprises:

[0022] a vertical fault throw restoring step of calculating the vertical fault throw of each fault in the structural map, and restoring the vertical fault throw of each fault according to the calculated vertical distance of each fault;

[0023] a horizontal fault throw restoring step of measuring the horizontal fault throw of each fault in the structural map, and restoring the horizontal fault throw of each fault according to the measured horizontal distance of each fault.

[0024] In some embodiments, in the vertical fault throw restoring step, the method of calculating the vertical fault throw of each fault is as follows: for the ith fault, the vertical fault throw of the ith fault is calculated by subtracting the value of the corresponding point on the upper plate of the ith fault from the value of the corresponding point on the lower plate of the ith fault.

[0025] The method of restoring the vertical fault throw of each fault according to the calculated vertical distance of each fault is as follows: for the ith fault, the vertical fault throw of the ith fault is restored by moving all the fault blocks in the direction opposite to the set direction of the ith fault upward along the vertical direction by the calculated vertical distance of the ith fault.

[0026] In some embodiments, in the horizontal fault throw restoring step, the method of measuring the horizontal fault throw of each fault in the structural map is as follows: for the ith fault, the horizontal fault throw of the ith fault is measured directly on the top surface structural map.

[0027] The method of restoring the horizontal fault throw of each fault according to the measured horizontal distance of each fault is as follows: for the ith fault, the horizontal fault throw of the ith fault is restored by moving all the fault blocks in the direction opposite to the set direction of the ith fault along the vertical direction of the ith fault by the measured horizontal distance of the ith fault.

[0028] In the second aspect, the present application provides a system for restoring the paleogeomorphology of a complex fault block area based on structural backstripping, which is used to implement the method for restoring the paleogeomorphology of a complex fault block area based on structural backstripping according to the first aspect of the present application, and comprises:

[0029] a data acquisition module configured to acquire a structural map of a target layer to be restored;

[0030] a fault numbering module configured to identify the faults in the structural map and sequentially number the faults according to a set direction to obtain a numbered structural map;

[0031] a fault throw restoring module configured to calculate the vertical fault throw and the horizontal fault throw of each fault in the numbered structural map, restore the vertical fault throw of each fault according to the vertical distance of each fault, and restore the horizontal fault throw of each fault according to the horizontal distance of each fault.

[0032] The paleogeomorphology recovery module sets a reference surface, calculates a distance TR from the top surface structure map after the fault distance recovery to the reference surface and a distance BR from the bottom surface structure map after the fault distance recovery to the reference surface, and obtains the paleogeomorphology shape surface after compaction by subtracting the distance TR from the distance BR.

[0033] In some embodiments, the method further comprises an abnormal point elimination module and a compaction correction module, the abnormal point elimination module is used for eliminating abnormal points in the paleogeomorphology shape surface after compaction, and the compaction correction module is used for calculating a compression coefficient of the target layer according to the porosity of the stratum during deposition and the porosity of the target layer, multiplying the calculated compression coefficient by the paleogeomorphology shape surface to obtain the paleogeomorphology after compaction correction.

[0034] Compared with the prior art, the method has the advantages and positive effects that:

[0035] (1) The method and system for recovering paleogeomorphology in a complex fault block area based on structure backstripping provided by the application fully considers the vertical fault distance and the horizontal fault distance, simultaneously calculates and measures the horizontal fault distance and the vertical fault distance of each fault, and quantitatively recovers the horizontal fault distance and the vertical fault distance of each fault in a certain order according to the fault blocks to obtain a paleogeomorphology quantitative characterization map in a geological history period, thereby eliminating the influence of the vertical fault distance and the horizontal fault distance on the paleogeomorphology, improving the paleogeomorphology recovery accuracy in the target layer deposition period, and improving the accuracy of the paleogeomorphology recovery.

[0036] (2) The method and system for recovering paleogeomorphology in a complex fault block area based on structure backstripping provided by the application further eliminates abnormal points in the paleogeomorphology after fault distance recovery, removes abnormal points caused by slight differences in the ranges of different surfaces, sampling intervals and the like, and further improves the paleogeomorphology recovery accuracy in the target layer deposition period.

[0037] (3) The method and system for recovering paleogeomorphology in a complex fault block area based on structure backstripping provided by the application further considers the influence of sediment compaction on the paleogeomorphology, performs compaction correction on the paleogeomorphology after fault distance recovery, further improves the paleogeomorphology recovery accuracy in the target layer deposition period, and can more accurately recover the paleogeomorphology in the target layer deposition period. BRIEF DESCRIPTION OF DRAWINGS

[0038] Figure 1 The flowchart of the method for recovering paleogeomorphology in a complex fault block area based on structure backstripping according to the first aspect of the application is shown in the figure;

[0039] Figure 2 The flowchart of the method for identifying faults in a structure map according to the embodiment of the application is shown in the figure;

[0040] Figure 3 The flowchart of the method for sequentially numbering faults in a set direction to obtain a numbered structure map according to the embodiment of the application is shown in the figure;

[0041] Figure 4 This is a flowchart illustrating the displacement recovery process according to an embodiment of the present invention;

[0042] Figure 5 This is a structural block diagram of the paleogeographic restoration system for complex fault-block areas based on tectonic back-exfoliation, as described in the second aspect of the present invention.

[0043] Figure 6 This is a flowchart of the paleogeographic restoration method for complex fault-block areas based on tectonic back-exfoliation, as described in the third aspect of the present invention.

[0044] Figure 7 This is a structural block diagram of the paleogeographic restoration system for complex fault-block areas based on tectonic back-exfoliation, as described in the fourth aspect of the present invention.

[0045] Figure 8 This is a flowchart of the paleogeographic restoration method for complex fault-block areas based on tectonic back-exfoliation, as described in the fifth aspect of the present invention.

[0046] Figure 9 This is a structural block diagram of the paleogeographic restoration system for complex fault-block areas based on tectonic back-exfoliation, as described in the sixth aspect of the present invention.

[0047] Figure 10 This is a flowchart of the paleogeographic restoration method for complex fault-block areas based on tectonic back-exfoliation, as described in the seventh aspect embodiment of the present invention;

[0048] Figure 11 This is a structural block diagram of the paleogeographic restoration system for complex fault-block regions based on tectonic back-exfoliation, as described in the eighth aspect embodiment of the present invention.

[0049] Figure 12 This is a top surface structural diagram of the target layer in an embodiment of the present invention;

[0050] Figure 13 for Figure 12 A three-dimensional display diagram;

[0051] Figure 14 This is a bottom surface structural diagram of the target layer in an embodiment of the present invention;

[0052] Figure 15 for Figure 14 A three-dimensional display diagram;

[0053] Figure 16 This is a top surface structural diagram of the target layer with fault numbers, according to an embodiment of the present invention.

[0054] Figure 17 This is a bottom structural diagram of the target layer with fault numbers, according to an embodiment of the present invention.

[0055] Figure 18 This is a structural diagram of the top surface of the target layer after the vertical displacement of fault F1 has been restored according to an embodiment of the present invention.

[0056] Figure 19 is a target layer bottom surface structure map after vertical fault throw recovery of the F1 fault; Figure 18

[0057] Figure 20 is a target layer bottom surface structure map after vertical fault throw recovery of the F1 fault;

[0058] Figure 21 is a target layer bottom surface structure map after vertical fault throw recovery of the F1 fault; Figure 20

[0059] Figure 22 is a target layer top surface structure three-dimensional display map after vertical fault throw recovery of the F1 fault, the F2 fault and the F3 fault;

[0060] Figure 23 is a target layer bottom surface structure three-dimensional display map after vertical fault throw recovery of the F1 fault, the F2 fault and the F3 fault;

[0061] Figure 24 is a target layer top surface structure map after horizontal fault throw recovery of the F1 fault;

[0062] Figure 25 is a target layer bottom surface structure map after horizontal fault throw recovery of the F1 fault;

[0063] Figure 26 is a target layer top surface structure three-dimensional display map after horizontal fault throw recovery of the F1 fault, the F2 fault and the F3 fault;

[0064] Figure 27 is a target layer bottom surface structure three-dimensional display map after horizontal fault throw recovery of the F1 fault, the F2 fault and the F3 fault;

[0065] Figure 28 is a compaction paleogeomorphology map;

[0066] Figure 29 is a porosity-depth statistical relationship map of a certain area;

[0067] Figure 30 is a compaction correction coefficient plane distribution map of a certain area;

[0068] Figure 31 is a compaction-corrected paleogeomorphology map of a certain area;

[0069] Figure 32 is a target layer bottom surface structure map after vertical fault throw recovery of the F1 fault; Figure 31

[0070] In the figure, 1 is a data acquisition module, 2 is a fault number module, 3 is a fault throw recovery module, 4 is a paleogeomorphology recovery module, 5 is an abnormal point elimination module, and 6 is a compaction correction module. DETAILED DESCRIPTION​​​

[0071] The present invention will now be described in detail through exemplary embodiments. However, it should be understood that, without further description, elements, structures, and features in one embodiment may be advantageously incorporated into other embodiments.

[0072] Complex fault-block regions, due to fault development, divide the target layer into multiple distinct fault blocks on the plane, with differences in horizontal and vertical fault displacements between these blocks. However, current restoration methods struggle to effectively address the fault displacement issues caused by fault activity when dealing with complex fault-block regions, resulting in low accuracy in paleogeographic reconstruction. This invention provides a method and system for paleogeographic reconstruction of complex fault-block regions based on tectonic stripping. It fully considers both vertical and horizontal fault displacements, simultaneously calculating and measuring the horizontal and vertical displacements of each fault. By dividing the fault blocks in a specific order, the horizontal and vertical displacements of each fault are quantitatively reconstructed, resulting in a quantitative paleogeographic representation map of geological history. This eliminates the influence of vertical and horizontal fault displacements on paleogeography, resolving the planar position errors caused by horizontal fault displacements and the vertical position errors caused by vertical fault displacements, thus achieving high accuracy in paleogeographic reconstruction. The following detailed description, in conjunction with the accompanying drawings, provides a method and system for paleogeographic reconstruction of complex fault-block regions based on tectonic stripping.

[0073] See Figure 1 According to a first aspect of the present invention, a method for paleogeographic restoration of complex fault-block areas based on tectonic back-exfoliation is provided, the steps of which are as follows:

[0074] S1. Data acquisition steps: Obtain the structural map of the target layer to be recovered, which includes a top structural map and a bottom structural map.

[0075] S2. Fault numbering steps: Identify the faults in the structural map and number them sequentially according to the set direction to obtain a numbered structural map.

[0076] Specifically, see Figure 2 The method for identifying interrupted layers in the construction graph is as follows:

[0077] S211. Graphic construction steps: Construct a fault polygonal planar graphic, and construct a top polygonal layer and a bottom polygonal layer in the planar graphic;

[0078] S212, Copying steps: Copy the fault polygons in the top surface construction diagram to the top surface polygon layer, and copy the polygons in the bottom surface construction diagram to the bottom surface polygon layer.

[0079] S213. Identification steps: In the top and bottom fault polygons, fault polygons with similar lengths and positions are considered as the same fault.

[0080] Specifically, referring to Figure 3 , the method for sequentially numbering the faults in the set direction to obtain the numbered structure map is:

[0081] S221, layer construction step: constructing a fault numbering layer in the planar graph;

[0082] S222, numbering step: sequentially numbering the identified faults in the set direction and placing them in the fault numbering layer;

[0083] S223, copying step: copying the fault numbering layer after numbering to the top and bottom structure maps to obtain the coded structure map.

[0084] S3, fault throw recovery step: recovering the vertical fault throw of each fault according to the calculated vertical fault throw of each fault in the numbered structure map, and recovering the horizontal fault throw of each fault according to the measured horizontal fault throw of each fault in the numbered structure map.

[0085] Specifically, in some embodiments, referring to Figure 4 , the fault throw recovery step includes:

[0086] S31, vertical fault throw recovery step: calculating the vertical fault throw of each fault in the structure map, and recovering the vertical fault throw of each fault according to the calculated vertical distance of each fault.

[0087] Specifically, the method for calculating the vertical fault throw of a fault is: for the ith fault, taking the values of the corresponding points on the upper and lower walls on both sides of the ith fault, subtracting the corresponding point value of the upper wall from the corresponding point value of the lower wall to obtain the vertical fault throw of the ith fault.

[0088] Specifically, the method for recovering the vertical fault throw of a fault according to the calculated vertical distance of the fault is: for the ith fault, fixing the outermost edges of the adjacent two set directions, moving all the fault blocks in the opposite direction of the set direction to the ith fault upward by the calculated vertical fault throw distance of the ith fault to recover the vertical fault throw of the ith fault.

[0089] S32, horizontal fault throw recovery step: measuring the horizontal fault throw of a fault in the structure map, and recovering the horizontal fault throw of each fault according to the measured horizontal fault throw.

[0090] Specifically, the method for measuring the horizontal fault throw of a fault in the structure map is: for the ith fault, directly measuring the horizontal fault throw of the ith fault on the top structure map.

[0091] Specifically, the method for recovering the horizontal fault displacement of the fault according to the measured horizontal fault displacement segment block is as follows: for the i-th fault, fixing the most edge of the two adjacent set directions, moving all the fault blocks in the set direction opposite to the i-th fault along the direction perpendicular to the i-th fault by the measured horizontal fault displacement distance of the i-th fault.

[0092] It should be noted that the order of the above two steps can be interchanged.

[0093] S4, a palaeogeomorphology recovery step: setting a reference surface, respectively calculating the distance TR from the top surface structure map after fault displacement recovery to the reference surface and the distance BR from the bottom surface structure map after fault displacement recovery to the reference surface, and obtaining the compaction palaeogeomorphology surface by subtracting the distance TR from the distance BR, i.e., the recovered palaeogeomorphology.

[0094] The above-mentioned palaeogeomorphology recovery method based on structure backstripping in the complex fault block area of the present application, when recovering the palaeogeomorphology in the complex fault block area, fully considers the vertical fault displacement and the horizontal fault displacement, simultaneously calculates and measures the vertical fault displacement and the horizontal fault displacement of each fault, and quantitatively recovers the vertical fault displacement and the horizontal fault displacement of each fault according to a certain order and by fault block, so as to obtain the quantitative characterization map of the palaeogeomorphology in the geological history period, eliminate the influence of the vertical fault displacement and the horizontal fault displacement on the palaeogeomorphology, and improve the palaeogeomorphology recovery precision in the sedimentary period of the target layer, and the palaeogeomorphology recovery accuracy is high.

[0095] Referring to Figure 5 The second aspect embodiment of the present application provides a palaeogeomorphology recovery system based on structure backstripping in a complex fault block area, which is used for realizing the palaeogeomorphology recovery method based on structure backstripping in a complex fault block area in the first aspect of the present application, and comprises:

[0096] A data acquisition module 1 acquires a structure map of a target layer to be recovered;

[0097] A fault numbering module 2 identifies the faults in the structure map, and sequentially numbers the faults according to a set direction to obtain a numbered structure map;

[0098] A fault displacement recovery module 3 calculates the vertical fault displacement and the horizontal fault displacement of each fault in the numbered structure map, recovers the vertical fault displacement of each fault according to the vertical distance of each fault by fault block, and recovers the horizontal fault displacement of each fault according to the horizontal fault displacement of each fault by segment block;

[0099] A palaeogeomorphology recovery module 4 sets a reference surface, respectively calculates the distance TR from the top surface structure map after fault displacement recovery to the reference surface and the distance BR from the bottom surface structure map after fault displacement recovery to the reference surface, and obtains the compaction palaeogeomorphology surface by subtracting the distance TR from the distance BR.

[0100] The above-mentioned complex fault block area paleogeomorphology recovery system based on structure backstripping fully considers the vertical fault throw and horizontal fault throw of the fault, simultaneously calculates and measures the horizontal fault throw and vertical fault throw of each fault, and quantitatively recovers the horizontal fault throw and vertical fault throw of each fault according to a certain order and in a fault block, so that the quantitative characterization map of the paleogeomorphology in the geological history period is obtained, the influence of the vertical fault throw and horizontal fault throw of the fault on the paleogeomorphology is eliminated, the recovery accuracy of the paleogeomorphology in the sedimentary period of the target layer is improved, and the paleogeomorphology recovery accuracy is high.

[0101] Referring to Figure 6 In a third aspect, the present application provides a complex fault block area paleogeomorphology recovery method based on structure backstripping, which comprises the following steps:

[0102] S1, a data acquisition step: acquiring a structure map of a target layer to be recovered, wherein the structure map comprises a top surface structure map and a bottom surface structure map.

[0103] S2, a fault numbering step: identifying the faults in the structure map, and sequentially numbering the faults according to a set direction to obtain a numbered structure map.

[0104] S3, a fault throw recovery step: recovering the vertical fault throw of each fault in a fault block according to the calculated vertical fault throw of each fault in the numbered structure map, and recovering the horizontal fault throw of each fault in a fault block according to the measured horizontal fault throw of each fault in the numbered structure map.

[0105] S4, a paleogeomorphology recovery step: setting a reference surface, calculating the distance TR from the top surface structure map after fault throw recovery to the reference surface and the distance BR from the bottom surface structure map after fault throw recovery to the reference surface, and obtaining the paleogeomorphology shape surface after compaction by subtracting the distance TR from the distance BR.

[0106] S5, an abnormal point elimination step: eliminating the abnormal points in the paleogeomorphology shape surface after compaction, wherein the abnormal points are points that do not obviously conform to the trend. Since the calculation of the paleogeomorphology shape surface after compaction involves arithmetic operations of multiple surfaces, there may be slight differences in the ranges and sampling intervals of different surfaces, which may lead to the appearance of abnormal points. Therefore, the paleogeomorphology shape surface after recovery is subjected to abnormal point elimination, and the abnormal points that do not obviously conform to the trend are deleted. It should be noted that if there are no points that do not conform to the trend, this step can be omitted.

[0107] It should be noted that in the embodiment, in step S2, the method of identifying the faults in the structure map and sequentially numbering the faults in the set direction to obtain the numbered structure map is the same as the method of restoring the paleogeomorphology in the complex fault block area based on structure backstripping described in the first aspect of the present application, and will not be repeated here. In step S3, the method of calculating the vertical fault throw of the fault and restoring the vertical fault throw of the fault according to the vertical distance calculated by the fault is the same as the method of restoring the paleogeomorphology in the complex fault block area based on structure backstripping described in the first aspect of the present application, and the method of measuring the horizontal fault throw of the fault in the structure map and restoring the horizontal fault throw of the fault according to the measured horizontal fault throw is the same as the method of restoring the paleogeomorphology in the complex fault block area based on structure backstripping described in the first aspect of the present application, and will not be repeated here.

[0108] The above-mentioned method for restoring paleogeomorphology in a complex fault block area based on structure backstripping fully considers the vertical fault throw and the horizontal fault throw of the fault, simultaneously calculates and measures the horizontal fault throw and the vertical fault throw of each fault, and quantitatively restores the horizontal fault throw and the vertical fault throw of each fault according to a certain order to obtain a paleogeomorphology quantitative characterization map in a geological history period, thereby eliminating the influence of the vertical fault throw and the horizontal fault throw of the fault on the paleogeomorphology, improving the accuracy of the paleogeomorphology restoration in the target layer deposition period, and improving the accuracy of the paleogeomorphology restoration. At the same time, when an abnormal point appears, the abnormal point obviously not meeting the trend is deleted, and the accuracy of the paleogeomorphology restoration is further improved. Compared with the method for restoring paleogeomorphology in a complex fault block area based on structure backstripping described in the first aspect of the present application and the system for restoring paleogeomorphology in a complex fault block area based on structure backstripping described in the second aspect of the present application, the accuracy of the paleogeomorphology restoration is higher.

[0109] Referring to Figure 7 The fourth aspect of the present application provides a system for restoring paleogeomorphology in a complex fault block area based on structure backstripping, which is used to realize the method for restoring paleogeomorphology in a complex fault block area based on structure backstripping described in the third aspect of the present application, and comprises:

[0110] The data acquisition module 1 acquires a structure map of a target layer to be restored;

[0111] The fault numbering module 2 identifies the faults in the structure map and sequentially numbers the faults in the set direction to obtain a numbered structure map;

[0112] The fault throw restoration module 3 calculates the vertical fault throw and the horizontal fault throw of each fault in the numbered structure map, restores the vertical fault throw of each fault according to the vertical distance of each fault, and restores the horizontal fault throw of each fault according to the horizontal fault throw of each fault;

[0113] The paleogeomorphology recovery module 4 sets a reference surface, calculates a distance TR from the top surface structure map after the fault distance is recovered to the reference surface, and a distance BR from the bottom surface structure map after the fault distance is recovered to the reference surface, and obtains a paleogeomorphology shape surface after compaction by subtracting the distance TR from the distance BR.

[0114] The abnormal point elimination module 5 is used for eliminating abnormal points in the paleogeomorphology shape surface after compaction.

[0115] The above-mentioned complex fault block area paleogeomorphology recovery system based on structure backstripping fully considers the vertical fault distance and the horizontal fault distance, simultaneously calculates and measures the horizontal fault distance and the vertical fault distance of each fault, and quantitatively recovers the horizontal fault distance and the vertical fault distance of each fault in a certain order according to the fault block, so as to obtain a paleogeomorphology quantitative characterization map in a geological history period, eliminate the influence of the vertical fault distance and the horizontal fault distance on the paleogeomorphology, improve the paleogeomorphology recovery precision in the target layer deposition period, and improve the paleogeomorphology recovery accuracy. Meanwhile, when abnormal points appear, the abnormal points obviously not meeting the trend are deleted, so as to further improve the paleogeomorphology recovery accuracy. Compared with the complex fault block area paleogeomorphology recovery method based on structure backstripping in the first aspect and the complex fault block area paleogeomorphology recovery system based on structure backstripping in the second aspect, the paleogeomorphology recovery accuracy is higher.

[0116] Reference Figure 8 The fifth aspect embodiment of the present application provides a complex fault block area paleogeomorphology recovery method based on structure backstripping, and the steps are as follows:

[0117] S1, a data acquisition step: acquiring a structure map of a target layer to be recovered, wherein the structure map includes a top surface structure map and a bottom surface structure map.

[0118] S2, a fault numbering step: identifying faults in the structure map, and sequentially numbering the faults according to a set direction to obtain a numbered structure map.

[0119] S3, a fault distance recovery step: recovering the vertical fault distance of each fault in the numbered structure map according to the calculated vertical fault distance of each fault, and recovering the horizontal fault distance of each fault according to the measured horizontal fault distance of each fault.

[0120] S4, a paleogeomorphology recovery step: setting a reference surface, calculating a distance TR from the top surface structure map after the fault distance is recovered to the reference surface, and a distance BR from the bottom surface structure map after the fault distance is recovered to the reference surface, and obtaining a paleogeomorphology shape surface after compaction by subtracting the distance TR from the distance BR.

[0121] S5, a compaction correction step: calculating a compression coefficient of the target layer according to the porosity during stratum deposition and the porosity of the target layer, multiplying the calculated compression coefficient by the paleogeomorphology shape surface to obtain a paleogeomorphology after compaction correction.

[0122] Specifically, by correlating the core porosity in the well with the depth, the initial porosity of the study area and the relationship between the porosity and the depth are established. According to the existing research results, the formation porosity decreases regularly with the increase of the depth, which is expressed by the following formula:

[0123]

[0124] In the formula, is the porosity of the target layer, is the porosity when the formation is deposited, is a constant, is the burial depth of the target layer.

[0125] It is generally believed that the stock price is incompressible during the compaction process of the target layer, therefore, the compression coefficient is expressed by the following formula:

[0126]

[0127] In the formula, is the compression coefficient. It should be noted that the compression coefficient is a value that changes with the plane position.

[0128] It should be noted that in this embodiment, in step S2, the method of identifying the faults in the structure map and sequentially numbering the faults according to the set direction to obtain the numbered structure map is the same as the method of restoring the paleogeomorphology of the complex fault block area based on structure backstripping described in the first aspect of the present application, which will not be repeated here. In step S3, the method of calculating the vertical fault throw of the fault and restoring the vertical fault throw of the fault according to the vertical distance calculated by the fault is the same as the method of restoring the paleogeomorphology of the complex fault block area based on structure backstripping described in the first aspect of the present application, the method of measuring the horizontal fault throw of the fault in the structure map and restoring the horizontal fault throw of the fault according to the measured horizontal fault throw is the same as the method of restoring the paleogeomorphology of the complex fault block area based on structure backstripping described in the first aspect of the present application, which will not be repeated here.

[0129] The above-mentioned method for restoring paleogeomorphology in a complex fault block area based on structural backstripping fully considers the vertical fault throw and the horizontal fault throw, simultaneously calculates and measures the horizontal fault throw and the vertical fault throw of each fault, and quantitatively restores the horizontal fault throw and the vertical fault throw of each fault in a certain order according to the fault blocks, so as to obtain a quantitative characterization map of the paleogeomorphology in a geological history period, eliminate the influence of the vertical fault throw and the horizontal fault throw on the paleogeomorphology, improve the restoration accuracy of the paleogeomorphology in the sedimentary period of the target layer, and improve the accuracy of the paleogeomorphology restoration. Meanwhile, the influence of compaction on the paleogeomorphology is considered, the paleogeomorphology after the fault throw restoration is corrected according to compaction, the influence of compaction on the paleogeomorphology is eliminated, and the accuracy of the paleogeomorphology restoration is further improved. Compared with the method for restoring paleogeomorphology in a complex fault block area based on structural backstripping in the first aspect and the system for restoring paleogeomorphology in a complex fault block area based on structural backstripping in the second aspect, the accuracy of the paleogeomorphology restoration is higher.

[0130] Referring to Figure 9 In a sixth aspect, the present application provides a system for restoring paleogeomorphology in a complex fault block area based on structural backstripping, which is used to implement the method for restoring paleogeomorphology in a complex fault block area based on structural backstripping in the fifth aspect, and comprises:

[0131] A data acquisition module 1 is configured to acquire a structure map of a target layer to be restored.

[0132] A fault numbering module 2 is configured to identify faults in the structure map, and sequentially number the faults according to a set direction to obtain a numbered structure map.

[0133] A fault throw recovery module 3 is configured to calculate the vertical fault throw and the horizontal fault throw of each fault in the numbered structure map, recover the vertical fault throw of each fault according to the vertical distance of each fault, and recover the horizontal fault throw of each fault according to the horizontal distance of each fault.

[0134] A paleogeomorphology restoration module 4 is configured to set a reference surface, calculate a distance TR from a top surface structure map after the fault throw recovery to the reference surface and a distance BR from a bottom surface structure map after the fault throw recovery to the reference surface, and obtain a paleogeomorphology shape surface by subtracting the distance TR from the distance BR.

[0135] A compaction correction module 6 is configured to calculate a compression coefficient of the target layer according to the porosity of strata during sedimentation and the porosity of the target layer, multiply the calculated compression coefficient by the paleogeomorphology shape surface to obtain the paleogeomorphology after compaction correction.

[0136] The above-mentioned complex fault block area paleogeomorphology recovery system based on structure backstripping fully considers the vertical fault throw and horizontal fault throw in the recovery of the complex fault block area paleogeomorphology, simultaneously calculates and measures the horizontal fault throw and vertical fault throw of each fault, and quantitatively recovers the horizontal fault throw and vertical fault throw of each fault in a certain order according to the fault block, to obtain a paleogeomorphology quantitative characterization map in a geological history period, eliminate the influence of the vertical fault throw and horizontal fault throw on the paleogeomorphology, improve the paleogeomorphology recovery accuracy of the target layer in the sedimentary period, and improve the paleogeomorphology recovery accuracy. Meanwhile, the influence of compaction on the paleogeomorphology is considered, the paleogeomorphology after the fault throw recovery is corrected, the influence of compaction on the paleogeomorphology is eliminated, and the paleogeomorphology recovery accuracy is further improved. Compared with the complex fault block area paleogeomorphology recovery method based on structure backstripping in the first aspect of the present application and the complex fault block area paleogeomorphology recovery system based on structure backstripping in the second aspect of the present application, the paleogeomorphology recovery accuracy is higher.

[0137] Reference Figure 10 The seventh aspect of the present application provides a complex fault block area paleogeomorphology recovery method based on structure backstripping, which comprises the following steps:

[0138] S1, a data acquisition step: acquiring a structure map of a target layer to be recovered, wherein the structure map comprises a top surface structure map and a bottom surface structure map.

[0139] S2, a fault numbering step: identifying faults in the structure map, and sequentially numbering the faults according to a set direction to obtain a numbered structure map.

[0140] S3, a fault throw recovery step: recovering the vertical fault throw of each fault in a segmented block according to the calculated vertical fault throw of each fault in the numbered structure map, and recovering the horizontal fault throw of each fault in a segmented block according to the measured horizontal fault throw of each fault in the numbered structure map.

[0141] S4, a paleogeomorphology recovery step: setting a reference surface, calculating a distance TR from the top surface structure map after the fault throw recovery to the reference surface and a distance BR from the bottom surface structure map after the fault throw recovery to the reference surface, and obtaining a paleogeomorphology shape surface by subtracting the distance TR from the distance BR.

[0142] S5, an abnormal point elimination step: eliminating abnormal points in the paleogeomorphology shape surface after compaction, wherein the abnormal points are points that do not obviously conform to the trend. Since the calculation of the paleogeomorphology shape surface after compaction involves multiple surface arithmetic operations, there may be slight differences in the ranges and sampling intervals of different surfaces, which may lead to the appearance of abnormal points. Therefore, the paleogeomorphology shape surface obtained after the paleogeomorphology recovery is subjected to abnormal point elimination, and the abnormal points that do not obviously conform to the trend are deleted. It should be noted that if there are no points that do not conform to the trend, this step can be omitted.

[0143] S6, compaction correction step: calculate the compression coefficient of the target layer according to the porosity of the stratum and the porosity of the target layer, and multiply the calculated compression coefficient by the paleogeomorphology surface to obtain the paleogeomorphology after compaction correction.

[0144] It should be noted that in the present embodiment, in step S2, the method of identifying the faults in the structure map and sequentially numbering the faults in the set direction to obtain the numbered structure map is the same as the method of restoring the paleogeomorphology of the complex fault block area based on structure backstripping described in the first aspect of the present application, and will not be repeated here. In step S3, the method of calculating the vertical fault throw of the fault and calculating the vertical fault throw of the fault according to the vertical distance of the fault is the same as the method of restoring the paleogeomorphology of the complex fault block area based on structure backstripping described in the first aspect of the present application, and the method of measuring the horizontal fault throw of the fault in the structure map and restoring the horizontal fault throw of the fault according to the measured horizontal fault throw is the same as the method of restoring the paleogeomorphology of the complex fault block area based on structure backstripping described in the first aspect of the present application, and will not be repeated here.

[0145] The above-mentioned method for restoring paleogeomorphology of complex fault block area based on structure backstripping of the present application, when restoring the paleogeomorphology of the complex fault block area, on the one hand, fully considers the vertical fault throw and the horizontal fault throw of the fault, and calculates and measures the horizontal fault throw and the vertical fault throw of each fault, and quantitatively restores the horizontal fault throw and the vertical fault throw of each fault according to a certain order, to obtain a paleogeomorphology quantitative characterization map in the geological history period, eliminates the influence of the vertical fault throw and the horizontal fault throw of the fault on the paleogeomorphology, improves the accuracy of the paleogeomorphology restoration of the target layer in the sedimentary period, and has high accuracy of paleogeomorphology restoration. On the other hand, when an abnormal point appears, the abnormal point that obviously does not conform to the trend is deleted, to further improve the accuracy of the paleogeomorphology restoration. On the other hand, the influence of compaction on the paleogeomorphology is also considered, and the paleogeomorphology after fault throw restoration is subjected to compaction correction, to eliminate the influence of compaction on the paleogeomorphology, and further improve the accuracy of the paleogeomorphology restoration. Compared with the method for restoring paleogeomorphology of complex fault block area based on structure backstripping of the first aspect, the third aspect and the fifth aspect of the present application and the system for restoring paleogeomorphology of complex fault block area based on structure backstripping of the second aspect, the fourth aspect and the sixth aspect of the present application, the accuracy of the paleogeomorphology restoration is higher.

[0146] Referring to Figure 11 The eighth aspect of the present application provides a system for restoring paleogeomorphology of complex fault block area based on structure backstripping, which is used to realize the method for restoring paleogeomorphology of complex fault block area based on structure backstripping of the seventh aspect of the present application, and comprises:

[0147] The data acquisition module 1 acquires a structure map of a target layer to be restored;

[0148] The fault numbering module 2 identifies the faults in the structure map, and sequentially numbers the faults in the set direction to obtain a numbered structure map;

[0149] The fault distance recovery module 3 calculates the vertical fault distance and the horizontal fault distance of each fault in the numbered structure map, recovers the vertical fault distance of each fault according to the vertical distance of each fault, and recovers the horizontal fault distance of each fault according to the horizontal fault distance of each fault;

[0150] The paleogeomorphology recovery module 4 sets a reference surface, calculates the distance TR from the top surface structure map after fault distance recovery to the reference surface and the distance BR from the bottom surface structure map after fault distance recovery to the reference surface, and obtains the paleogeomorphology curve by subtracting the distance TR from the distance BR;

[0151] The abnormal point elimination module 5 is used to eliminate abnormal points in the paleogeomorphology curve after compaction.

[0152] The compaction correction module 6 calculates the compression coefficient of the target layer according to the porosity of the stratum during deposition and the porosity of the target layer, and obtains the paleogeomorphology after compaction correction by multiplying the calculated compression coefficient by the paleogeomorphology curve.

[0153] The above-mentioned complex fault block area paleogeomorphology recovery method based on structure backstripping can, in the recovery of the paleogeomorphology of a complex fault block area, on the one hand, fully consider the vertical fault distance and the horizontal fault distance of the fault, simultaneously calculate and measure the horizontal fault distance and the vertical fault distance of each fault, quantitatively recover the horizontal fault distance and the vertical fault distance of each fault in a certain order, and obtain a paleogeomorphology quantitative characterization map of a geological history period, thereby eliminating the influence of the vertical fault distance and the horizontal fault distance of the fault on the paleogeomorphology, improving the accuracy of the paleogeomorphology recovery of the target layer deposition period, and improving the accuracy of the paleogeomorphology recovery. On the other hand, when abnormal points appear, the abnormal points that do not obviously conform to the trend are deleted, and the accuracy of the paleogeomorphology recovery is further improved. On the other hand, the influence of compaction on the paleogeomorphology is also considered, the paleogeomorphology after fault distance recovery is subjected to compaction correction, the influence of compaction on the paleogeomorphology is eliminated, and the accuracy of the paleogeomorphology recovery is further improved. Compared with the complex fault block area paleogeomorphology recovery method based on structure backstripping in the first aspect, the third aspect and the fifth aspect of the present application and the complex fault block area paleogeomorphology recovery system based on structure backstripping in the second aspect, the fourth aspect and the sixth aspect of the present application, the accuracy of the paleogeomorphology recovery is higher.

[0154] The effectiveness of the above-mentioned complex fault block area paleogeomorphology recovery method and system based on structure backstripping is described below in combination with specific embodiments.

[0155] Embodiment: Taking the paleogeomorphology recovery of a target layer in a certain research area as an example. The recovery steps are as follows:

[0156] S1, data acquisition: the top surface structure map of the target layer collected or drawn in a certain research area is shown in Figure 12 and Figure 13 , and the bottom surface structure map of the target layer collected or drawn is shown in Figure 14 and Figure 15.

[0157] S2, Fault Numbering: Identify the faults in the top and bottom structure maps, and sequentially number the faults from west to east to obtain the numbered structure maps. The numbered top structure map is shown in Figure 16 , and the numbered bottom structure map is shown in Figure 17 , Figure 16 , Figure 17 In the above, there are three faults, respectively numbered F1 fault, F2 fault, and F3 fault.

[0158] S3, Fault Displacement Recovery

[0159] S31, Vertical Fault Displacement Recovery: In the top structure map, read the values of the corresponding points on the hanging wall and foot wall of F1 fault. The value of the corresponding point on the foot wall E1 is -3595, and the value of the corresponding point on the hanging wall E2 is -3695. The vertical fault displacement of F1 fault is calculated as E1-E2=-3595-(-3695)=100m. Fix the edges of the top and bottom structure maps of the target layer in the westmost and northmost positions, and move all the fault blocks east or south of F1 fault in the vertical upward direction by a distance of 100m to recover the vertical fault displacement of F1 fault. The top structure map after F1 fault recovery is shown in Figure 18 and Figure 19 , and the bottom structure map after F1 fault recovery is shown in Figure 20 and Figure 21 . Then, sequentially recover the vertical fault displacement of F2 fault and F3 fault in the top and bottom structure maps to obtain the top structure map after vertical fault displacement recovery (see Figure 22 ) and the bottom structure map after vertical fault displacement recovery (see Figure 23 ).

[0160] S32, Horizontal Fault Displacement Recovery: In the top structure map, directly measure the horizontal fault displacement of F1 fault as 150m. Fix the edges of the top and bottom structure maps of the target layer in the westmost and northmost positions, and move all the fault blocks east or south of F1 fault in the vertical direction of F1 fault by a distance of 150m to obtain the top structure map after F1 fault recovery (see Figure 24 , and the bottom structure map after F1 fault recovery (see Figure 25 ). Then, sequentially recover the horizontal fault displacement of F2 fault and F3 fault in the top and bottom structure maps to obtain the top structure map after horizontal fault displacement recovery (see Figure 26 ) and the bottom structure map after horizontal fault displacement recovery (see Figure 27 ).

[0161] S4, Palaeogeomorphology Recovery: Set the reference surface RS as a constant 0, calculate the distance TR from the top structure map after fault displacement recovery to the reference surface RS, and the distance BR from the bottom structure map after fault displacement recovery to the reference surface RS. The palaeogeomorphology surface is obtained by subtracting the distance TR from the distance BR (seeFigure 28 ).

[0162] S5, compaction correction: by correlating the core porosity data with depth in the well (see Figure 29 ), the initial porosity and the porosity-depth relationship in the study area are established.

[0163]

[0164] The compression coefficient calculation formula is:

[0165]

[0166] In this embodiment, the initial porosity is set to 49%, and the target layer is buried at a depth of 2,000 m. The compression coefficient contour map calculated by the compression coefficient calculation formula is shown in FIG. 6. Figure 30 .

[0167] The paleogeomorphology after compaction correction is shown in FIG. 7 and FIG. 8. Figure 31 and Figure 32 .

[0168] The above embodiments are used to explain the present application, but not to limit the present application, and any modification and change made to the present application within the spirit and protection scope of the claims of the present application shall fall into the protection scope of the present application.

Claims

1. A method for paleogeomorphology restoration of a complex fault-block area based on structural backstripping, characterized in that, The steps are: Data acquisition step: obtain the structure map of the target layer to be recovered, which includes the top surface structure map and the bottom surface structure map; Fault numbering step: identify the faults in the structure map, and sequentially number the faults according to the set direction to obtain a numbered structure map; Fault throw recovery step: calculate the vertical throw of each fault in the numbered structure map, and divide the blocks to recover the vertical throw of each fault according to the measured horizontal throw of each fault in the numbered structure map; Paleogeomorphology recovery step: set a reference surface, calculate the distance TR from the top surface structure map after fault throw recovery to the reference surface, and the distance BR from the bottom surface structure map after fault throw recovery to the reference surface, and obtain the paleogeomorphology shape surface after compaction by subtracting the distance TR from the distance BR.

2. The method for palaeogeomorphology restoration based on the construction back-stripping of complex fault-block area according to claim 1, characterized in that, It also includes an abnormal point elimination step: eliminating abnormal points in the paleogeomorphology shape surface after compaction, which are points that do not obviously conform to the trend.

3. The method for palaeogeomorphology restoration based on structural backstripping of complex fault-block area according to claim 1 or 2, characterized in that, It also includes a compaction correction step: calculating the compression coefficient of the target layer according to the porosity during stratigraphic deposition and the porosity of the target layer, and multiplying the calculated compression coefficient by the paleogeomorphology shape surface to obtain the paleogeomorphology after compaction correction.

4. The method for palaeogeomorphology restoration based on the construction back-stripping of complex fault-block area according to claim 1, characterized in that, In the fault numbering step, the method for identifying faults in the structure map is: Graph construction step: constructing a fault polygon plane graph, and constructing a top surface polygon layer and a bottom surface polygon layer in the plane graph; Copy step: copying the fault polygon in the top surface structure map to the top surface polygon layer, and copying the polygon in the bottom surface structure map to the bottom surface polygon layer; Identification step: regarding fault polygons with similar lengths and positions as the same fault in the top surface fault polygon and the bottom surface fault polygon.

5. The method for palaeogeomorphology restoration based on the construction back-stripping of complex fault-block area according to claim 4, characterized in that, In the fault numbering step, the method for sequentially numbering the faults according to the set direction to obtain a numbered structure map is: Layer construction step: constructing a fault numbering layer in the plane graph; Numbering step: sequentially numbering the identified faults according to the set direction and placing them in the fault numbering layer; Copy step: copying the numbered fault numbering layer to the top surface structure map and the bottom surface structure map to obtain a numbered structure map.

6. The method for palaeogeomorphology restoration based on the construction back-stripping of complex fault-block area according to claim 1, wherein, The fault throw recovery step includes: Vertical throw recovery step: calculating the vertical throw of each fault in the structure map, and dividing the blocks to recover the vertical throw of each fault according to the calculated vertical distance of each fault; Horizontal throw recovery step: measuring the horizontal throw of a fault in the structure map, and dividing the blocks to recover the horizontal throw of each fault according to the measured horizontal throw.

7. The method for palaeogeomorphology restoration based on the construction back-stripping of complex fault-block area according to claim 6, characterized in that, In the vertical throw recovery step, the method for calculating the vertical throw of a fault is: for the ith fault, taking the values of the corresponding points on the upper and lower sides of the ith fault, subtracting the value of the corresponding point on the upper side from the value of the corresponding point on the lower side to obtain the vertical throw of the ith fault; The method for recovering the vertical throw of a fault by dividing the blocks according to the calculated vertical distance of the fault is: for the ith fault, fixing the most edge of the adjacent two set directions, and moving all the fault blocks in the direction opposite to the set direction to recover the vertical throw of the ith fault.

8. The method for palaeogeomorphology restoration based on the construction back-stripping of complex fault-block area according to claim 6, wherein, In the horizontal fault throw recovery step, the method for measuring the horizontal fault throw of the fault in the structure map is: for the ith fault, directly measuring the horizontal fault throw of the ith fault on the top surface structure map; The method for recovering the horizontal fault throw of the fault according to the segmented block is: for the ith fault, fixing the two most edge blocks with the set direction, moving all the blocks with the opposite direction of the set direction to the ith fault by the measured horizontal fault throw of the ith fault.

9. A system for restoring paleogeomorphology of a complex fault-block region based on tectonic backstripping, for implementing the method for restoring paleogeomorphology of a complex fault-block region based on tectonic backstripping according to any one of claims 1 to 8, characterized in that, Comprise: a data acquisition module for acquiring a structure map of a target layer to be recovered; a fault numbering module for identifying the faults in the structure map and sequentially numbering the faults according to a set direction to obtain a numbered structure map; a fault throw recovery module for calculating the vertical fault throw and the horizontal fault throw of each fault in the numbered structure map, recovering the vertical fault throw of each fault according to the segmented blocks, and recovering the horizontal fault throw of each fault according to the segmented blocks; a paleogeomorphology recovery module for setting a reference surface, calculating the distance TR from the top surface structure map after fault throw recovery to the reference surface and the distance BR from the bottom surface structure map after fault throw recovery to the reference surface, and obtaining the paleogeomorphology surface after compaction by subtracting the distance TR from the distance BR.

10. The system for palaeogeomorphology restoration based on structural backstripping of complex fault-block areas according to claim 9, wherein, Further comprising an abnormal point elimination module and a compaction correction module, the abnormal point elimination module is used for eliminating abnormal points in the paleogeomorphology surface after compaction; the compaction correction module is used for calculating the compression coefficient of the target layer according to the porosity during stratum deposition and the porosity of the target layer, multiplying the calculated compression coefficient by the paleogeomorphology to obtain the paleogeomorphology after compaction correction.

Images