Method and system for recovering ancient landform of syngenetic fracture development area

Abstract

The invention discloses a syngenetic fracture development area ancient landform recovery method and system, and belongs to the technical field of oil-gas exploration and development. According to the method, syngenetic fracture and horizon are explained by using three-dimensional seismic data, a three-dimensional geological framework model is established by combining the explained syngenetic fracture and horizon data, gridded seismic bedding plane and syngenetic fracture plane polygon data are further obtained, subtraction is performed on top and bottom bedding plane data to obtain the apparent thickness of the stratum, and the thickness of the stratum is calculated. According to the method, inclination angle correction is carried out on the apparent thickness of the stratum to obtain the true thickness of the stratum, abnormal value areas in seismic true thickness data are filtered and removed through syngenetic fracture plane polygons, finally plane interpolation is carried out on the data to obtain an accurate and reliable ancient landform plane graph, and recovery of the true thickness of the stratum on the periphery of the syngenetic fracture is achieved through the process. The method greatly improves the accuracy and precision of ancient landform recovery of the syngenetic fracture massive development region, and is of great significance to fine exploration and development and well location deployment.

Description

Technical Field

[0001] The present invention belongs to the technical field of oil and gas exploration and development, and particularly relates to a method and system for restoring paleogeomorphology in a syn - genetic fault - developed area. Background Technique

[0002] Paleogeomorphology is comprehensively affected by various geological conditions such as weathering, erosion, tectonic movement, and karstification. Restoring the true paleogeomorphology is of great significance for analyzing the distribution of sand bodies and the characteristics of tectonic evolution. Paleogeomorphology restoration refers to using seismic and logging data to restore the geomorphology of a selected stratum to before the deposition of that stratum, which is an important link in basin analysis, sedimentary facies analysis, and reservoir prediction.

[0003] At present, common paleogeomorphology restoration methods include: residual thickness method (residual thickness method), impression method, filling and leveling method, horizon flattening method, back - stripping method, etc. The basic principle of the residual thickness method is to select a certain stratum interface under the erosion surface as the reference surface, and indirectly reflect the paleogeomorphology by the residual thickness between the reference surface and the top surface of the weathered crust. The impression method is based on the principle of sedimentary compensation, and uses the thickness from the overlying stratum of the erosion surface to its horizontal reference surface to mirror - reflect the general shape of the paleogeomorphology. The areas with large thickness are low - potential areas of the paleogeomorphology, and the areas with small thickness are high - potential areas of the paleogeomorphology. The principles of the other several technical methods are basically similar to the impression method.

[0004] In recent years, a variety of paleogeomorphology methods have been developed or improved to be applicable to different geological conditions. For example, the Chinese invention patent application document with the application publication number CN108267797A discloses a method for restoring the paleogeomorphology in a syn - genetic normal - fault - developed area within a complex fault block. This method calculates the growth index of the syn - genetic normal fault, and within the control range of the down - thrown block of the syn - genetic normal fault, corrects the thickness of the down - thrown block by the ratio of the thickness of the down - thrown block to the growth index, eliminating the paleogeomorphology illusion caused by the thickening of the down - thrown - block strata due to the gradual activity of the syn - genetic normal fault during paleogeomorphology restoration, making the paleogeomorphology more reasonable. However, this method is only applicable to the syn - genetic normal - fault - developed area. For the reverse - fault - developed area, due to the existence of the superposition area of the hanging - wall and foot - wall strata, the thickness cannot be corrected by the growth index. Summary of the Invention

[0005] In view of the above problems, on the one hand, the present invention provides a method for restoring paleogeomorphology in a syn - genetic fault - developed area, and the restoration method includes the following steps:

[0006] Obtain the seismic horizon data and fault interpretation data in the depth domain of the target layer according to the 3D seismic data;

[0007] Set the stratum framework model parameters according to the seismic horizon data and fault interpretation data and establish a 3D geological framework model;

[0008] Obtain the syn - sedimentary fault plane polygon and seismic horizon data by using the three - dimensional geological framework model; the seismic horizon data includes the seismic top interface horizon and the seismic bottom interface horizon.

[0009] Obtain the formation apparent thickness of the target layer by using the seismic horizon data, and calculate the formation dip angle θ of the seismic top interface horizon of the target layer. Use the formation dip angle θ to correct the formation apparent thickness to obtain the true formation thickness.

[0010] Compile the fault polygon corresponding to the paleogeomorphology according to the syn - sedimentary fault plane polygon, and perform planar interpolation on the true formation thickness in combination with the fault polygon to obtain the paleogeomorphology plan.

[0011] Furthermore, obtaining the seismic horizon data and fault interpretation data in the depth domain of the target layer from the 3D seismic data includes the following steps:

[0012] Perform closed - loop interpretation on the structure of the syn - sedimentary fault development area of the target layer by using 3D seismic data.

[0013] Obtain the seismic horizon data and fault interpretation data in the depth domain of the target layer according to the interpretation results.

[0014] Among them, the seismic horizon data and fault interpretation data include the formation thickness of the hanging wall and the footwall of the syn - sedimentary fault.

[0015] Furthermore, setting the formation framework model parameters and establishing a three - dimensional geological framework model according to the seismic horizon data and fault interpretation data includes the following steps:

[0016] By means of 3D geological modeling, determine the geological framework model parameters by using the syn - sedimentary fault interpretation data and establish a fault framework model.

[0017] Determine the fault spatial structure relationship of the target layer according to the fault framework model.

[0018] Establish a three - dimensional geological framework model including the syn - sedimentary fault framework model according to the fault spatial structure relationship, and obtain the fault surface and seismic horizon data according to the syn - sedimentary fault framework model.

[0019] Furthermore, the geological framework model parameters include the seismic horizon grid cell size, the seismic horizon filtering radius at the periphery of the fault surface, the horizon smoothness, and the spatial interpolation algorithm.

[0020] The fault spatial structure relationship includes the contact relationships of the main fault, derivative faults, and secondary faults of the syn - sedimentary fault.

[0021] Furthermore, the specific method of obtaining the syn - sedimentary fault plane polygon and seismic horizon data by using the syn - sedimentary fault framework model is as follows:

[0022] Draw a polygon with the intersection line of the fracture surface and the bedding plane as the center line of the syn - sedimentary fracture plane polygon, and this polygon is the syn - sedimentary fracture plane polygon.

[0023] Use the syn - sedimentary fracture plane polygon to filter the seismic horizon data at the fracture periphery. After confirming that the filtered seismic horizon data meets the requirements of rationality and accuracy, obtain the gridded seismic horizon data.

[0024] Further, the specific method for obtaining the apparent thickness of the target layer formation using the seismic horizon data is as follows:

[0025] Subtract the seismic top - interface horizon from the seismic bottom - interface horizon to obtain the apparent thickness of the target layer formation.

[0026] Further, calculating the dip angle θ of the seismic top - interface horizon of the target layer and using the dip angle θ to correct the apparent thickness of the formation to obtain the true thickness of the formation includes the following steps:

[0027] Grid the seismic top - interface horizon and calculate the dip angle θ of each grid cell.

[0028] Use the dip angle θ to perform dip correction on the apparent thickness of the formation to obtain the true thickness of the formation.

[0029] Further, the following formula is used for performing dip correction on the apparent thickness of the formation using the dip angle θ:

[0030] True thickness of the formation = Apparent thickness of the formation * cosθ.

[0031] Further, specifically compiling the fracture polygon corresponding to the paleogeomorphology according to the syn - sedimentary fracture plane polygon is: Based on the mid - line of the syn - sedimentary fracture plane polygon, re - compile the paleogeomorphology fracture polygon.

[0032] In a second aspect, the present invention proposes a paleogeomorphology restoration system in a syn - sedimentary fracture development area, including:

[0033] A first acquisition unit for obtaining seismic horizon data and fracture interpretation data in the depth domain of the target layer according to 3D seismic data;

[0034] A building unit for setting formation framework model parameters according to the seismic horizon data and fracture interpretation data and building a syn - sedimentary fracture framework model;

[0035] A second acquisition unit for obtaining a syn - sedimentary fracture plane polygon and seismic horizon data using the syn - sedimentary fracture framework model; the seismic horizon data includes the seismic top - interface horizon and the seismic bottom - interface horizon;

[0036] A calculation and correction unit for obtaining the apparent thickness of the target layer formation using the seismic horizon data, calculating the formation dip angle θ of the seismic top interface horizon of the target layer, and correcting the apparent thickness of the formation using the formation dip angle θ to obtain the true thickness of the formation;

[0037] A restoration unit for preparing a fault polygon corresponding to the paleogeomorphology according to the syn-depositional fault plane polygon, and performing planar interpolation on the true thickness of the formation in combination with the fault polygon to obtain a paleogeomorphology plan view.

[0038] In a third aspect, the present invention provides an electronic device, including a processor, a communication interface, a memory, and a communication bus. Among them, the processor, the communication interface, and the memory complete communication with each other through the communication bus;

[0039] The memory stores a computer program;

[0040] The processor is configured to implement the syn-depositional fault development area paleogeomorphology restoration method when executing the program stored in the memory.

[0041] In a fourth aspect, the present invention provides a computer-readable storage medium storing a computer program, and when the computer program is run, it executes the syn-depositional fault development area paleogeomorphology restoration method.

[0042] Advantages of the present invention:

[0043] The present invention interprets syn-depositional faults and horizons using 3D seismic data, establishes a 3D geological framework model based on the jointly interpreted syn-depositional fault and horizon data, further obtains gridded seismic horizon and syn-depositional fault plane polygon data, then subtracts the top and bottom horizon data to obtain the apparent thickness of the formation, corrects the apparent thickness of the formation for dip to obtain the true thickness of the formation, filters out the abnormal value area in the seismic true thickness data using the syn-depositional fault plane polygon, and finally performs planar interpolation on this data to obtain an accurate and reliable paleogeomorphology plan view. The above process realizes the restoration of the true formation thickness around the syn-depositional fault, greatly improves the accuracy and precision of paleogeomorphology restoration in areas with a large number of syn-depositional faults, and is of great significance for fine exploration and development and well location deployment.

[0044] Other features and advantages of the present invention will be described in the subsequent specification, and part of them will become obvious from the specification or be understood by implementing the present invention. The objectives and other advantages of the present invention can be achieved and obtained through the structures pointed out in the specification, claims, and drawings. Description of the Drawings

[0045] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the accompanying drawings required for the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other accompanying drawings can be obtained based on these drawings.

[0046] Figure 1 The flowchart of the paleogeomorphological restoration method in the contemporaneous fault development area proposed by the present invention is shown;

[0047] Figure 2 The pattern diagrams of contemporaneous normal faults and contemporaneous reverse faults proposed in the present invention are shown;

[0048] Figure 3 The schematic diagram of the formation thickness profile after filtering the fault polygons in the present invention is shown;

[0049] Figure 4 The schematic diagram of the calculation result of the apparent thickness of the contemporaneous fault formation in the present invention is shown;

[0050] Figure 5 The schematic diagram of the restored true paleogeomorphological feature profile in the present invention is shown;

[0051] Figure 6 The three-dimensional diagram of the restored true paleogeomorphological feature in the present invention is shown;

[0052] Figure 7 The cross-sectional diagram of the development of contemporaneous faults in the upper member of the lower Ganchaigou Formation (E3 2 ) in the Qaidam Basin in the embodiment of the present invention is shown;

[0053] Figure 8 The framework model and horizon data of the contemporaneous faults in the Qaidam Basin in the embodiment of the present invention are shown;

[0054] Figure 9 The fault polygons and the seismic horizon data after grid division in the embodiment of the present invention are shown;

[0055] Figure 10a The plan view of the apparent formation thickness in the embodiment of the present invention is shown;

[0056] Figure 10b The plan view of the formation dip angle at the top surface of the target layer in the embodiment of the present invention is shown;

[0057] Figure 11 The final paleogeomorphological plan view in the embodiment of the present invention is shown;

[0058] Figure 12 The schematic diagram of the paleogeomorphological restoration system in the contemporaneous fault development area proposed in the present invention is shown;

[0059] Figure 13 Shows a schematic diagram of an electronic device proposed in an embodiment of the present invention. Specific embodiments

[0060] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0061] The present invention proposes a method for restoring paleogeomorphology in the coeval fault development area, as Figure 1 shown, including the following steps:

[0062] S1: Obtain seismic horizon data and fault interpretation data in the depth domain of the target layer according to 3D seismic data; specifically:

[0063] First, use 3D seismic data to perform closed interpretation on the structure of the coeval fault development area of the target layer;

[0064] Obtain high-precision seismic horizons of the top and bottom interfaces of the target layer and coeval fault interpretation data according to the interpretation results.

[0065] Among them, the seismic horizon data and fault interpretation data are horizon data and fault data interpreted from the 3D seismic data volume. For example, it can be data manually traced according to the seismic event characteristics of 4×8 bins. If it is time-domain seismic data, the time-domain horizons can be converted to the depth domain through time-depth conversion. If it is depth-domain seismic horizons and fault data, they are directly used. In the present invention, to simplify the process steps, the conventional time-depth domain conversion steps are omitted, and the subsequent steps are uniformly implemented according to the depth-domain result data.

[0066] S2: Set the formation framework model parameters according to the seismic horizon data and fault interpretation data and establish a 3D geological framework model; the 3D geological framework model includes a coeval fault framework model and a non-coeval fault framework model;

[0067] Specifically, it includes the following steps:

[0068] Through 3D geological modeling means, use the coeval fault interpretation data to determine the geological framework model parameters, and establish a fault framework model according to the geological framework model parameters;

[0069] Determine the fault spatial structure relationship of the target layer according to the fault framework model; the fault spatial structure relationship includes contact relationships such as coeval fault main faults, derivative faults, and secondary faults;

[0070] Establish a 3D geological framework model based on the fracture spatial structure relationship. The 3D geological framework model includes a syndepositional fracture framework model and a non-syndepositional fracture framework model. According to the syndepositional fracture framework model, obtain fracture surfaces with reasonable structure and high spatial consistency.

[0071] The parameters of the geological framework model are used to define the spatial structure and accuracy of the geological framework model, and at least include the following four types: the size of the seismic horizon grid unit, the filtering radius of the seismic horizons at the periphery of the fracture surface, the smoothness of the horizons, and the spatial interpolation algorithm. The spatial interpolation algorithm can be inverse distance weighting or Kriging, etc.

[0072] Among them, the definition of a syndepositional fracture is a fracture with different thicknesses in the hanging wall and footwall. The thickness of the hanging wall is the subtraction thickness of the top and bottom horizons of the hanging wall of the fracture, and the thickness of the footwall is the subtraction thickness of the top and bottom horizons of the footwall of the fracture. To determine whether a fracture is a syndepositional fracture, first, it is necessary to determine the development of the fracture at this position based on seismic data interpretation, and second, based on the top and bottom horizons interpreted from seismic data, judge the thicknesses of the hanging wall and footwall of the fracture. A syndepositional normal fault is one where the hanging wall drops and the footwall rises, and the formation thicknesses of the hanging wall and footwall are inconsistent; a syndepositional reverse fault is one where the hanging wall rises and the footwall drops, and the formation thicknesses of the hanging wall and footwall of the fracture are inconsistent, specifically as Figure 2 shown Figure 2 shows the schematic diagrams of syndepositional normal faults and syndepositional reverse faults.

[0073] S3: Use the syndepositional fracture framework model to obtain syndepositional fracture plane polygons and seismic horizon data. Specifically:

[0074] The syndepositional fracture plane polygon is the position where the intersection of the fracture surface and the horizon in 3D is projected onto the plane, that is, the intersection line of the fracture surface and the horizon. For the convenience of daily mapping, draw a polygon around this intersection line, and this intersection line is the center line of the syndepositional fracture plane polygon;

[0075] Use the syndepositional fracture plane polygon to filter the seismic horizon data at the periphery of the fracture and check the rationality and accuracy of the filtered seismic horizon data to obtain the gridded seismic horizon data.

[0076] As Figure 3 shown, shows the schematic diagram of the formation thickness profile after filtering the seismic horizon data at the periphery of the fracture using the fracture polygon. For a normal fault, the hanging wall drops and the footwall rises; for a reverse fault, the hanging wall of the fracture rises and the footwall drops.

[0077] S4: Use the seismic horizon data to obtain the apparent thickness of the target layer formation; the seismic horizon data includes the seismic top interface horizon and the seismic bottom interface horizon;

[0078] Specifically, the initial top and bottom formation apparent thickness of the target layer is obtained by subtracting the seismic top interface horizon from the seismic bottom interface horizon. The formula for the seismic apparent thickness is: formation apparent thickness = seismic bottom interface horizon - seismic top interface horizon. Among them, the seismic bottom interface horizon is the data obtained after interpolation of the seismic bottom interface horizon data within the geological framework model grid, and the seismic top interface horizon is the data obtained after interpolation of the seismic top interface horizon data within the grid of the geological framework model.

[0079] As Figure 4 shown, the calculation results of the formation apparent thickness in the cases of syn-sedimentary normal faults and syn-sedimentary reverse faults in the present invention are shown. Among them, a thinning section appears between the hanging wall and the footwall of the syn-sedimentary normal fault, and a thickening section appears between the hanging wall and the footwall of the syn-sedimentary reverse fault.

[0080] S5: Calculate the formation dip angle of the seismic top interface horizon of the target layer, and correct the formation apparent thickness using the formation dip angle θ to obtain the true formation thickness;

[0081] After meshing the seismic top interface horizon, calculate the formation dip angle θ of each grid cell;

[0082] Use the formation dip angle θ to perform dip angle correction on the formation apparent thickness to obtain the true formation thickness.

[0083] The formula for the true formation thickness is: true formation thickness = formation apparent thickness * cosθ.

[0084] S6: Compile the fault polygon corresponding to the paleogeomorphology according to the syn-sedimentary fault plane polygon obtained in step S3; specifically, the paleogeomorphology fault polygon can be re-compiled based on the midline of the syn-sedimentary fault plane polygon.

[0085] S7: Combine the fault polygon and perform planar interpolation on the true formation thickness to obtain the paleogeomorphology plan. As Figure 5 and 6 shown, the restored true paleogeomorphology characteristic profile and three-dimensional diagram obtained by using the restoration method of the present invention are respectively shown, and the faults disappear.

[0086] The numerical accuracy at the fault periphery of this paleogeomorphology plan is relatively high, and it can be more accurately used to analyze the reservoir sedimentary facies and the distribution characteristics of sand bodies.

[0087] To more clearly illustrate the paleogeomorphology restoration method in the syn-sedimentary fault development area proposed by the present invention, this embodiment takes the study of the paleogeomorphology restoration of syn-sedimentary faults in the upper member of the Lower Ganchaigou Formation (E3 2 ) in the Qiekerike area of the Qaidam Basin as an example to elaborate on the above method in detail. As Figure 7As shown, multiple groups of syndepositional faults have developed in this area, such as the three groups of syndepositional reverse faults shown in the figure. In this embodiment, relevant paleogeomorphic restoration is carried out relying on the newly processed prestack depth migration seismic data of Qiekelike, which specifically includes the following steps:

[0088] Step 1: Use 3D seismic data to obtain seismic horizon data and fault interpretation data in the depth domain of the target layer. Focus on interpreting the top and bottom horizons of the target layer section and the syndepositional faults in the study area E3 2 The top and bottom horizons of the target layer section and the syndepositional faults are interpreted. The bin distribution of the interpretation is 10×20 (i.e., 150m×300m) according to lines and traces.

[0089] Step 2: Set the parameters of the stratigraphic framework model according to the seismic horizon data and fault interpretation data in the depth domain of the target layer, and establish a 3D geological framework model;

[0090] Specifically, establish a fault system model according to the fault interpretation data. First, determine the spatial structural relationship of the faults, and focus on determining the contact relationships such as the main faults, derivative faults, and secondary faults of the syndepositional faults;

[0091] Establish a 3D geological framework model including the syndepositional fault framework model according to the structural relationship, and obtain a fault surface with reasonable structure and high spatial consistency.

[0092] The syndepositional fault framework model of Qiekelike in the Qaidam Basin and the seismic horizon data established in this embodiment are as Figure 8 shown. The figure shows the original data for establishing the 3D geological framework model, that is, the seismic horizon data. The grid-like data is the 3D seismic interpretation horizon data, and the fault data is shown as a fault surface. Since the fault has a main fault and its derivative faults, and they intersect with each other, the contact relationship between the faults needs to be determined for the fault surface data, so that the derivative fault stops extending to the main fault.

[0093] Step 3: Use the results of the 3D geological framework model to obtain the planar polygon of the syndepositional fault and the seismic horizon data.

[0094] Use the planar polygon of the syndepositional fault to filter the seismic horizon data around the fault, and check the rationality and accuracy of the filtered data to obtain the gridded seismic horizon data, specifically as Figure 9 shown. The seismic horizon data includes the seismic bottom interface horizon and the seismic top interface horizon. This figure shows the seismic horizon data obtained from the 3D geological framework model. The red line data is actually a line composed of polygons. This polygon is the fault polygon, which is the line where the fault intersects the plane. When the polygon becomes the horizon data of this layer, the unreasonable data inside the polygon is hollowed out.

[0095] Step 4: Obtain the apparent thickness of the target formation using the seismic horizon data obtained in Step 3. Specifically, in this embodiment, the apparent thickness of the top and bottom strata of the target formation is obtained by subtracting the seismic top horizon from the seismic bottom horizon, and the obtained apparent formation thickness is as Figure 10a shown. Since the data of the seismic top horizon and the bottom horizon of the seismic interface have been filtered using the fault polygon previously, the abnormal areas of thinning or thickening around the syn-depositional faults in the formation thickness map are filtered out.

[0096] The formula for the apparent thickness of the top and bottom strata of the target formation is: Apparent formation thickness = Seismic bottom horizon - Seismic top horizon.

[0097] Step 5: Calculate the formation dip angle of the seismic top horizon of the target formation, and correct the apparent formation thickness to obtain the true formation thickness. In this embodiment, the formation dip angle plan of the seismic top horizon of the target formation obtained is as Figure 10b shown. The formation dip angle of the generated horizon of the formation top interface in this figure. Since the apparent formation thickness is obtained by subtracting the seismic top horizon from the seismic bottom horizon, if the formation dip angle is 0°, the subtraction is the true thickness. However, if the formation dip angle is 45 degrees, the apparent thickness obtained by subtracting the upper and lower horizons is significantly greater than the thickness when the dip angle is 0°. Therefore, it is necessary to calculate the dip angle for dip angle correction in order to convert the apparent formation thickness into the true thickness.

[0098] After gridifying the seismic top horizon, calculate the formation dip angle θ of each grid unit, and perform dip angle correction on the initial apparent formation thickness to obtain the true formation thickness. There are still problems with abnormal values of thinning or thickening around the faults in this formation thickness, but the true formation thickness values in other areas far from the faults are reliable and true.

[0099] The formula for the true formation thickness is: True formation thickness = Apparent formation thickness * cosθ.

[0100] Step 6: Compile the fault polygon corresponding to the paleogeomorphology based on the syn-depositional fault plane polygon obtained in Step 3. Based on the midline of the syn-depositional fault plane polygon, recompile the paleogeomorphology fault polygon.

[0101] Step 7: Combine the paleogeomorphology fault polygon, perform planar interpolation on the true formation thickness, and obtain the final accurate and reliable paleogeomorphology plan. The final paleogeomorphology plan is as Figure 11 shown, where the black line is the determined final fault polygon. First, correct the apparent formation thickness to obtain the true formation thickness through the above steps, and then jointly interpolate using the fault polygon and the true thickness data to fill in the data that was hollowed out and deleted due to the abnormal value area. The numerical accuracy around the faults in this paleogeomorphology plan is relatively high, and it can be used more accurately to analyze the reservoir sedimentary facies and sand body distribution characteristics.

[0102] Based on the same inventive concept of the present invention, another exemplary embodiment of the present invention provides a paleogeomorphology restoration system in a syn-sedimentary fault development area, as Figure 12 shown, including:

[0103] A first acquisition unit 1201, configured to acquire seismic horizon data and fault interpretation data in the depth domain of the target layer according to 3D seismic data;

[0104] A building unit 1202, configured to set formation framework model parameters and build a 3D geological framework model according to the seismic horizon data and fault interpretation data;

[0105] A second acquisition unit 1203, configured to acquire syn-sedimentary fault plane polygons and seismic horizon data by using the 3D geological framework model; the seismic horizon data includes a seismic top interface horizon and a seismic bottom interface horizon;

[0106] A calculation and correction unit 1204, configured to acquire the apparent thickness of the target layer formation by using the seismic horizon data, calculate the formation dip angle θ of the seismic top interface horizon of the target layer, and correct the apparent thickness of the formation by using the formation dip angle θ to obtain the true thickness of the formation;

[0107] A restoration unit 1205, configured to compile fault polygons corresponding to the paleogeomorphology according to the syn-sedimentary fault plane polygons, and perform planar interpolation on the true thickness of the formation in combination with the fault polygons to obtain a paleogeomorphology plan view.

[0108] Based on the same inventive concept, another exemplary embodiment of the present invention provides an electronic device. As Figure 13 shown, the electronic device includes at least one processor 1301, at least one communication interface 1302, at least one memory 1303, and at least one communication bus 1304; wherein, the processor 1301, the communication interface 1302, and the memory 1303 communicate with each other through the communication bus 1304;

[0109] The memory 1303 stores a computer program;

[0110] The processor 1301, when executing the program stored in the memory 1303, implements the paleogeomorphology restoration method in the syn-sedimentary fault development area.

[0111] Optionally, the communication interface may be the interface of a communication module, such as the interface of a GSM module; the processor may be a processor CPU, or a specific integrated circuit ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement the embodiments of the present invention. The memory may include high-speed RAM memory and may also include non-volatile memory, such as at least one disk memory. Among them, the memory stores a program, and the processor invokes the program stored in the memory to execute some or all of the above method embodiments.

[0112] Based on the same inventive concept, an embodiment of the present application further provides a computer-readable storage medium storing a computer program, which when run, implements some or all of the above method embodiments. Optionally, the storage medium may be a non-transitory computer-readable storage medium. For example, the non-transitory computer-readable storage medium may be ROM, random access memory (RAM), CD-ROM, magnetic tape, floppy disk, and optical data storage devices, etc.

[0113] Although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: they can still modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements for some of the technical features; and these modifications or replacements 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 the present invention.

Claims

1. A method for restoring paleogeomorphology in the synsedimentary fault development area, characterized in that The restoration method includes the following steps: Obtain seismic horizon data and fault interpretation data in the depth domain of the target layer based on 3D seismic data; Set the parameters of the stratigraphic framework model according to the seismic horizon data and fault interpretation data, and establish a 3D geological framework model; Obtain syn-depositional fault plane polygons and seismic horizon data by using the 3D geological framework model; the seismic horizon data includes the seismic top interface horizon and the seismic bottom interface horizon; Obtain the apparent thickness of the strata in the target layer by using the seismic horizon data, calculate the dip angle θ of the seismic top interface horizon of the target layer, and correct the apparent thickness of the strata by using the dip angle θ to obtain the true thickness of the strata; Compile fault polygons corresponding to the paleogeomorphology according to the syn-depositional fault plane polygons, and perform planar interpolation on the true thickness of the strata in combination with the fault polygons to obtain a paleogeomorphology plan view.

2. The paleogeomorphology restoration method for the contemporaneous fault development area according to claim 1, characterized in that The step of obtaining seismic horizon data and fault interpretation data in the depth domain of the target layer based on 3D seismic data includes the following steps: Perform closed-loop interpretation on the structure of the syn-depositional fault development area of the target layer by using 3D seismic data; Obtain seismic horizon data and fault interpretation data in the depth domain of the target layer according to the interpretation results; Among them, the seismic horizon data and fault interpretation data include the thickness of the strata in the hanging wall and the footwall of the syn-depositional fault.

3. The paleogeomorphic restoration method for the contemporaneous fault development area according to claim 2, wherein The step of setting the parameters of the stratigraphic framework model according to the seismic horizon data and fault interpretation data and establishing a 3D geological framework model includes the following steps: By means of 3D geological modeling, determine the parameters of the geological framework model by using syn-depositional fault interpretation data and establish a fault framework model; Determine the fault spatial structure relationship of the target layer according to the fault framework model; Establish a 3D geological framework model including a syn-depositional fault framework model according to the fault spatial structure relationship, and obtain fault surfaces and seismic horizon data according to the syn-depositional fault framework model.

4. The paleogeomorphic restoration method for the contemporaneous fault development area according to claim 3, characterized in that The parameters of the geological framework model include the size of the seismic horizon grid unit, the filtering radius of the seismic horizon around the fault surface, the smoothness of the horizon, and the spatial interpolation algorithm; The fault spatial structure relationship includes the contact relationship between the main fault, the derivative fault, and the secondary fault of the syn-depositional fault.

5. The paleogeomorphology restoration method for the contemporaneous fault-developed area according to claim 3, wherein Specifically, obtaining syn-depositional fault plane polygons and seismic horizon data by using the syn-depositional fault framework model is as follows: Draw a polygon with the intersection line of the fault surface and the horizon as the center line of the syn-depositional fault plane polygon, which is the syn-depositional fault plane polygon; Filter the seismic horizon data around the fault by using the syn-depositional fault plane polygon, and obtain the gridded seismic horizon data after confirming that the filtered seismic horizon data meets the requirements of reasonableness and accuracy.

6. The paleogeomorphic restoration method for the contemporaneous fault-developed area according to claim 1, wherein Specifically, obtaining the apparent thickness of the strata in the target layer by using the seismic horizon data is as follows: Subtract the seismic bottom interface horizon from the seismic top interface horizon to obtain the apparent thickness of the strata in the target layer.

7. The paleogeomorphic restoration method for the contemporaneous fault-developed area according to claim 1 or 6, characterized in that The steps of calculating the dip angle θ of the seismic top interface horizon of the target layer and correcting the apparent thickness of the strata by using the dip angle θ to obtain the true thickness of the strata include the following steps: Grid the seismic top interface horizon and calculate the dip angle θ of each grid unit; Correct the apparent thickness of the strata by using the dip angle θ to obtain the true thickness of the strata.

8. The paleogeomorphology restoration method for the contemporaneous fault development area according to claim 7, characterized in that, The following formula is used to correct the apparent thickness of the strata by using the dip angle θ: True formation thickness = apparent formation thickness * cosθ.

9. The paleogeomorphology restoration method for the contemporaneous fault-developed area according to claim 1, characterized in that, Specifically, the fault polygon corresponding to the paleogeomorphology is compiled according to the synsedimentary fault plane polygon as follows: Based on the midline of the synsedimentary fault plane polygon, the paleogeomorphology fault polygon is re-compiled.

10. A paleogeomorphology restoration system in a synsedimentary fault development area, characterized in that, It includes: A first acquisition unit, configured to acquire seismic horizon data and fault interpretation data in the depth domain of the target layer according to 3D seismic data; A building unit, configured to set formation framework model parameters and build a synsedimentary fault framework model according to the seismic horizon data and fault interpretation data; A second acquisition unit, configured to acquire a synsedimentary fault plane polygon and seismic horizon data by using the synsedimentary fault framework model; the seismic horizon data includes the seismic top interface horizon and the seismic bottom interface horizon; A calculation and correction unit, configured to acquire the apparent formation thickness of the target layer by using the seismic horizon data, calculate the formation dip angle θ of the seismic top interface horizon of the target layer, and correct the apparent formation thickness by using the formation dip angle θ to obtain the true formation thickness; A restoration unit, configured to compile a fault polygon corresponding to the paleogeomorphology according to the synsedimentary fault plane polygon, and perform planar interpolation on the true formation thickness in combination with the fault polygon to obtain a paleogeomorphology plan view.

11. An electronic device, characterized in that, It includes a processor, a communication interface, a memory, and a communication bus. Among them, the processor, the communication interface, and the memory complete communication with each other through the communication bus; The memory stores a computer program; The processor is configured to, when executing the program stored in the memory, implement the method for restoring the paleogeomorphology in the synsedimentary fault development area according to any one of claims 1-9.

12. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is run, it executes the method for restoring the paleogeomorphology in the synsedimentary fault development area according to any one of claims 1-9.

Images