CN115508884B - Method, device and system for restoring overlapped basin area structure and application thereof - Google Patents
Method, device and system for restoring overlapped basin area structure and application thereof Download PDFInfo
- Publication number
- CN115508884B CN115508884B CN202110696240.XA CN202110696240A CN115508884B CN 115508884 B CN115508884 B CN 115508884B CN 202110696240 A CN202110696240 A CN 202110696240A CN 115508884 B CN115508884 B CN 115508884B
- Authority
- CN
- China
- Prior art keywords
- horizon
- regional
- difference value
- distance difference
- horizons
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000010586 diagram Methods 0.000 claims abstract description 70
- 238000010276 construction Methods 0.000 claims description 31
- 238000012545 processing Methods 0.000 claims description 21
- 238000011084 recovery Methods 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 11
- 238000005253 cladding Methods 0.000 claims 1
- 238000010030 laminating Methods 0.000 claims 1
- 239000010410 layer Substances 0.000 description 64
- 239000007789 gas Substances 0.000 description 14
- 230000015572 biosynthetic process Effects 0.000 description 13
- 238000005755 formation reaction Methods 0.000 description 13
- 150000003839 salts Chemical class 0.000 description 10
- 238000004590 computer program Methods 0.000 description 7
- 238000004220 aggregation Methods 0.000 description 5
- 230000002776 aggregation Effects 0.000 description 5
- 238000011161 development Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 238000012876 topography Methods 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000019771 cognition Effects 0.000 description 3
- 230000008021 deposition Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 208000010392 Bone Fractures Diseases 0.000 description 1
- 208000007356 Fracture Dislocation Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000001364 causal effect Effects 0.000 description 1
- 239000011247 coating layer Substances 0.000 description 1
- 230000009918 complex formation Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004898 kneading Methods 0.000 description 1
- 238000003475 lamination Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method, a device and a system for restoring regional structures of overlapping basins and application thereof, wherein the method comprises the steps of selecting at least two regional standard horizons in a regional seismic profile structure interpretation chart; respectively leveling the selected regional standard horizons to obtain a leveled regional seismic section structure interpretation map; further, determining at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon according to the at least two flattened regional seismic profile structure interpretation diagrams; and finally, recovering the structure of the overlapped basin area according to the distance difference value between at least one reference horizon and the reference horizon. By the method, the influence of the late structural deformation in the overlying stratum on the early stratum in the underlying stratum can be reduced, the false structure and the false thickness of the stratum caused by the late structural deformation are greatly reduced, and the specific form of the early stratum with different structural periods is restored.
Description
Technical Field
The invention relates to the technical field of petroleum and natural gas exploration and development, in particular to a method, a device and a system for restoring a regional structure of a superimposed basin and application thereof.
Background
In the oil and gas industry, reservoirs are enriched within a particular location of a set of formations in the subsurface today, but the location of the enriched reservoir is constantly moving and adjusting during the geologic history. Since the entire process of oil and gas production, migration and aggregation is well understood in modern oil and gas exploration. Therefore, the basin structure evolution needs to be recovered and analyzed in time, and the structure characteristics of different periods in the geological history and the burial depth condition of the target layer are known. Such research is necessary and urgent, especially for ancient folded basins with long evolution times.
The overlapped basin refers to a basin with a complex structure, wherein the original basin is longitudinally overlapped in sequence due to different deformation of different structural layers under the action of complex dynamics mechanisms of multi-period, multi-direction and multi-scale changes. The basin in each period has its own independent prototype, its structure pattern, geometric pattern and deformation mode are different from those in other periods, and the superposition of different periods shows the gradual reconstruction and superposition evolution process of prototype basin. Therefore, the analysis of the dynamics mechanism of the overlapped basin and the evolution history thereof and the recovery of the stress state of the basin in each period are beneficial to researching the causal mechanism of the geometric and kinematic characteristics of the basin construction in different periods, and have important significance for guiding the oil and gas exploration deployment.
Disclosure of Invention
Based on the knowledge of the prior art, the inventor of the present application generally deploys two-dimensional seismic data in overlapping basins, and usually adopts horizon tracking and key layer leveling of two-dimensional sections to implement early construction morphology. However, complex formation zones such as salt formations, steep formation zones, etc. exist in the basin, and very significant false formations and false formation thicknesses can occur with simple horizon flattening, which severely affects later hydrocarbon exploration studies. And the inventors have found that the above phenomenon is caused by the fact that late complex construction affects the recovery of early construction, and that if the late horizons are simply leveled, large errors occur in the underlying early horizons.
The present invention has been made in view of the above problems, and it is an object of the present invention to provide a method, apparatus, system and application thereof for restoring a superimposed basin area configuration that overcomes or at least partially solves the above problems.
In a first aspect, an embodiment of the present invention provides a method for restoring a overlapped basin area structure, which may include:
Selecting at least two regional standard horizons from the regional seismic profile construction interpretation map;
Respectively carrying out horizon flattening treatment on the selected regional standard horizons to obtain a flattened regional seismic profile structure interpretation map;
Determining at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon according to at least two flattened regional seismic section structure interpretation graphs;
and recovering the overlapped basin area structure according to at least one reference horizon and the distance difference value corresponding to the reference horizon.
Optionally, the determining, according to the at least two flattened regional seismic profile construction interpretation graphs, at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon may include:
vectorizing at least two flattened regional seismic section structure explanatory diagrams to obtain at least two flattened regional seismic section vector diagrams;
Determining at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon in at least two flattened regional seismic section vector diagrams;
and the distance difference value corresponding to the reference horizon is the distance difference value between the reference horizon determined by the later-stage regional standard horizon and the corresponding reference horizon determined by the earlier-stage regional standard horizon in the leveled regional seismic section vector diagram.
Optionally, determining the distance difference value corresponding to the reference horizon may include:
Determining a plurality of coordinate points representing curve forms corresponding to the reference horizon in the regional seismic section vector diagram corresponding to the advanced regional standard horizon;
And taking the vertical distance from each coordinate point to the curve corresponding to the reference horizon in the regional seismic profile vector diagram corresponding to the early regional standard horizon in the vertical direction as the distance difference value corresponding to the reference horizon.
Optionally, the restoring the overlapping basin area structure according to at least one of the reference horizons and the distance difference value corresponding to the reference horizon may include:
Shifting the lower layer of the reference layer up by the distance difference value corresponding to the reference layer in the regional seismic profile vector diagram corresponding to the advanced regional standard layer based on the distance difference value corresponding to the reference layer;
And overlapping the moved early reference layer and the lower layer thereof with the upper layer of the later reference layer.
Optionally, after the lamination process, the method may further include: and filling the overlapped vector images to generate the overlapped basin geological evolution section.
Optionally, selecting at least two regional standard horizons in the regional seismic profile structure interpretation graph may include: selecting at least two regional standard horizons according to a preset standard;
The preset criteria include at least one of: the method comprises the steps of widely distributed horizons, horizons with gentle paleolandform corresponding to geological periods, horizons with strong seismic reflection surface amplitude and horizons with extended regions.
In a second aspect, an embodiment of the present invention provides a device for restoring a superimposed basin area structure, which may include:
the selecting module is used for selecting at least two regional standard horizons from the regional seismic profile structure interpretation graph;
the leveling module is used for respectively leveling the selected regional standard horizons to obtain a leveled regional seismic section structure interpretation graph;
The determining module is used for determining at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon according to at least two flattened regional seismic section structure interpretation graphs;
And the recovery module is used for recovering the overlapped basin area structure according to at least one reference horizon and the distance difference value corresponding to the reference horizon.
In a third aspect, embodiments of the present invention provide a superimposed basin area construction restoration system, which may include: a first subsystem and a second subsystem;
the first subsystem is used for selecting at least two regional standard horizons from the regional seismic profile construction interpretation map; respectively carrying out horizon flattening treatment on the selected regional standard horizons to obtain flattened regional seismic section structure interpretation diagrams;
The second subsystem is used for determining at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon according to at least two flattened regional seismic profile construction interpretation graphs; and recovering the overlapped basin area structure according to at least one reference horizon and the distance difference value corresponding to the reference horizon.
In a fourth aspect, an embodiment of the present invention provides an application of a stacked basin geological evolution section obtained by the stacked basin region structure restoration method according to the first aspect in geological modeling.
In a fourth aspect, embodiments of the present invention provide a geologic modeling system that may include a superimposed basin region structure restoration device as described in the second aspect.
The technical scheme provided by the embodiment of the invention has the beneficial effects that at least:
The embodiment of the invention provides a method, a device and a system for restoring regional structures of a superimposed basin and application thereof, wherein the method comprises the steps of firstly selecting at least two regional standard horizons from a regional seismic profile structure interpretation chart; respectively leveling the selected regional standard horizons to obtain a leveled regional seismic section structure interpretation map; further, determining at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon according to the at least two flattened regional seismic profile structure interpretation diagrams; and finally, recovering the structure of the overlapped basin area according to the distance difference value between at least one reference horizon and the reference horizon. By the method, the influence of the late structural deformation in the overlying stratum on the early stratum in the underlying stratum can be reduced, the false structure and the false thickness of the stratum caused by the late structural deformation are greatly weakened, and the actual early stratum concrete form of different structural periods is restored.
Furthermore, the shape of the reference layer in the underlying stratum refers to the shape of the reference layer of the overlying stratum, so that the deformation of the reference layer is real, a plurality of layers below the underlying stratum reference layer are determined to be real in deformation at the same time, and the stratum recovery is more accurate. And a detailed data basis is provided for cognition and favorable zone evaluation of oil and gas aggregation rules in oil and gas exploration.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims thereof as well as the appended drawings.
The technical scheme of the invention is further described in detail through the drawings and the embodiments.
Drawings
The accompanying drawings are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate the invention and together with the embodiments of the invention, serve to explain the invention. In the drawings:
FIG. 1 is a flow chart of a method for restoring a superimposed basin area configuration provided in an embodiment of the present invention;
FIG. 2 is a flow chart of another method for restoring a superimposed basin area configuration provided in an embodiment of the present invention;
FIG. 3 is an illustrative view of a regional seismic profile configuration provided in an embodiment of the invention;
FIG. 4 is an explanatory view of a regional seismic profile configuration for the region standard horizon (T 3x1) of FIG. 3 after it has been leveled;
FIG. 5 is a vector diagram of the vectorized area seismic section of FIG. 4;
FIG. 6 is a schematic diagram of early and late stage structural evolution profiles versus a determined reference horizon (P 2l1) provided in an embodiment of the present invention;
FIG. 7 is a schematic diagram of determining a plurality of coordinate points characterizing a curve morphology in comparison of late and early structural evolution profiles according to an embodiment of the present invention;
FIG. 8 is a schematic diagram showing the effect of the early structural evolution section provided in the present invention after the overall displacement of all stratum horizons according to the distance differences between the reference horizons;
FIG. 9 is a schematic diagram showing the overall effect of the early stage structure evolution profile after the early stage structure evolution profile is replaced and adjusted according to the embodiment of the present invention;
FIG. 10 is a schematic diagram of a cross-section of geological evolution of a superimposed basin after filling and rendering according to an embodiment of the present invention;
FIG. 11 is a schematic structural diagram of a folded basin region structure restoration device according to an embodiment of the present invention.
Detailed Description
Exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
< Explanation of terms >
Two-dimensional earthquake: the two-dimensional seismic exploration method is to arrange a measuring line on the ground, perform seismic exploration construction along each measuring line, collect seismic wave information reflected by the underground stratum back to the ground, and then obtain a Zhang Dezhen section through electronic computer processing. Two-dimensional seismic exploration techniques are commonly used to determine hydrocarbon reservoir locations, to determine hydrocarbon drilling locations, and to display subsurface geologic formations in two dimensions (length and depth) in a geologically interpreted seismic profile as if a knife were cut down from the surface. Meanwhile, tens of intersecting two-dimensional measuring lines are used together, so that the fluctuation situation of the earth surface before the deposition in a certain geological period in the underground can be compiled. If it is found which places are likely to store oil and gas, it can be determined to be the oil and gas drilling well site.
Regional seismic profile construction interpretation: the interpretation result diagram of the regional seismic section structure is seismic geological data which can find out the change of the structure form of the underground geologic body (including geological phenomena such as stratum horizon, fault and the like) and combine the section plane with the plane for spatial interpretation.
The construction tape: geological structure phenomena such as faults, folds and the like which have common causes and internal relations exist in a certain range.
Salt structure (salt structure): the salt deposit buried deeply under the ground is caused by the load of the overlying giant thick cover layer to cause the density inversion so that the salt flows and floats upwards, and the salt layer is forced to be extruded upwards into the overlying deposit in some sections, so that the structure deformation of rock stratum and surrounding rock, such as salt dome (salt bottom) and rock stratum flow kneading folds, is the salt structure.
High steep structural band: fracture dislocation, which is less strongly likely to occur in the fracture of the structure, is not prominent with respect to the gentle structure, and is called a gentle structure, whereas it is called a steep structure, and its constituent structural band is a steep structural band.
False construction: in the seismic profile construction interpretation, under the conditions of salt dome, local igneous rock mass, reverse-buried fault lower disc, steep dip angle of overlying strata or rapid thickness change of the overlying strata, etc., the structural deformation or structural artifact of the underlying strata is caused on the profile. In the embodiment of the invention, after the standard stratum is leveled in the seismic profile structure interpretation chart, the underlying stratum is partially recessed into a V shape, which is caused by the fact that the overlying stratum has an incorrect structure caused by a steep structural band and is not actually present.
Horizon: fully-defined as a formation horizon refers to a particular location in the formation sequence. There are many kinds of stratum horizons, such as lithology horizons with special lithology characteristics, fossil horizons with special fossils, chronologic horizons with specific times, seismic horizons, electrical logging horizons, etc. Thus, the horizon of a formation may be a boundary of a formation unit, or may be a marker layer belonging to a specific age.
Construction lattice (tectonic framework): refers to a structural framework for controlling the spatial layout of various geologic bodies in a region, and the structural framework often forms a certain structural style. In areas where deformation of the multi-stage structure occurs, the primary structure is generally the primary structure constituting the basic structure lattice thereof. The spatial arrangement format of the different building units or building strips is either on the regional building scale or at some regional building evolution stage.
Constructing a grid section: refers to a profile of a line across a certain zone of construction (e.g., basin, area break, etc.).
Example 1]
The embodiment of the invention provides a method for restoring the structure of a region of a superimposed basin, which can comprise the following steps with reference to fig. 1:
And S11, selecting at least two regional standard horizons in the regional seismic profile structure interpretation graph.
The above-described regional seismic section structural explanatory diagram in the embodiment of the invention is a structural explanatory diagram of a structural grid section. Since the present invention is directed to a method of zone construction restoration for overlapping basins, in this embodiment is a seismic profile measured on one line across the basin or zone of construction, such as a grid seismic profile across a Sichuan basin. The regional standard horizon refers to a stratum horizon which is widely distributed in the whole overlapped basin, clear in wave resistance relation, strong in traceability, less in fault influence and gentle in stratum thickness change in the regional seismic section structure interpretation diagram.
In this embodiment, the interpretation map of the seismic profile of the region may be a bitmap image in order to facilitate image processing of the interpretation map of the seismic profile of the region. Of course, if the steps described below in this embodiment can be performed, the image format may be other formats, which are not particularly limited in this embodiment of the present invention.
And step S12, respectively carrying out horizon leveling treatment on the selected regional standard horizons to obtain a leveled regional seismic section structure interpretation graph.
In the step, the selected regional standard horizon is used as a reference, and horizon flattening processing operation is carried out on the regional seismic profile structure interpretation map. Specifically, the layer leveling process in this embodiment is a data processing process performed by leveling the time of a reflection horizon, which is subjected to seismic interpretation, to a given time position on a seismic section or in a seismic data volume, and maintaining the amplitude value of the reflection horizon unchanged.
In this embodiment, when the horizon flattening process is performed on the above-mentioned region standard horizon, the top time and the ground time of the region standard horizon need to be set to be as large as possible, and at least the earliest stratum horizon and the latest stratum horizon which satisfy the interpretation of the structure of the seismic section of the region can be displayed in the seismic section.
In the embodiment of the invention, as the horizon flattening processing is carried out on the interpretation graph of the regional seismic profile structure, a plurality of false structures are generated in the underlying horizon profile below the regional standard horizon, and the underlying horizon is locally concave into a V shape, which are all false structures caused by the existence of high-steep structural bands and the like in the overlying stratum and are not in practice. The object of the invention is to eliminate the above-mentioned false construction and the corresponding false layer thickness.
And S13, determining at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon according to the at least two flattened regional seismic section structure interpretation diagrams.
Firstly, the section comparison is carried out on the flattened regional seismic section structure interpretation map obtained in the step S12, the reference layers in at least two sections are determined, and then the distance difference between the reference layers is calculated.
And S14, recovering the structure of the overlapped basin area according to at least one reference horizon and the distance difference value corresponding to the reference horizon.
In the step, old stratum (stratum layer) below the reference layer in the later-stage flattened regional seismic section structure interpretation map is deleted, and stratum below the reference layer in the earlier-stage flattened regional seismic section structure interpretation map is replaced by stratum below the reference layer, namely, the earlier-stage reference layer and the lower layer thereof are overlapped with the upper coating layer of the later-stage reference layer, so that the recovery of the overlapped basin region structure is completed. The shape of the reference layer in the underlying stratum refers to the shape of the reference layer of the overlying stratum, so that the deformation of the reference layer is real, a plurality of layers below the underlying stratum reference layer are determined to be deformed real at the same time, and the stratum recovery is more accurate.
According to the method for restoring the regional structure of the overlapped basin provided by the embodiment of the invention, at least two regional standard horizons are selected from an interpretation chart of the regional seismic profile structure; respectively leveling the selected regional standard horizons to obtain a leveled regional seismic section structure interpretation map; further, determining at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon according to the at least two flattened regional seismic profile structure interpretation diagrams; and finally, recovering the structure of the overlapped basin area according to the distance difference value between at least one reference horizon and the reference horizon. By the method, the influence of the late structural deformation in the overlying stratum on the early stratum in the underlying stratum can be reduced, the false structure and the false thickness of the stratum caused by the late structural deformation are greatly weakened, and the actual early stratum concrete form of different structural periods is restored.
Furthermore, a detailed data foundation is provided for cognition and favorable zone evaluation of oil and gas aggregation rules in oil and gas exploration.
In a specific embodiment, in the step S11, at least two regional standard horizons are selected from the regional seismic profile structure interpretation map, specifically, at least two regional standard horizons are selected according to a preset standard; the preset criteria include at least one of: the method comprises the steps of widely distributed horizons, horizons with gentle paleolandform corresponding to geological periods, horizons with strong seismic reflection surface amplitude and horizons with extended regions.
It should be noted that, according to the above criteria, generally, more than three area standard levels can be selected in the interpretation chart of the area seismic section structure, and in the embodiment of the present invention, the restoration method of the area structure claimed in the present invention can be implemented only by selecting two area standard levels. Of course, three or more region standard levels may be selected, when step S12 is executed, leveling processing needs to be performed on the three or more region standard levels, when step S13 is executed, three or more flattened structure explanatory diagrams need to be selected to determine a reference level included in the structure explanatory diagrams, and when a distance difference corresponding to the reference level is determined, an average value of the distance differences may be determined by weighting and averaging.
Of course, the regional standard layers in this embodiment are sequenced in the geological history, and the distance difference corresponding to the reference layer needs to be determined by the section flattened in the early and late stages when the step S13 is performed, which should not be ambiguous to those skilled in the art.
Example 2]
In another specific embodiment, referring to fig. 2, the overlapping basin area configuration restoration method may specifically include the following steps:
and S21, selecting at least two regional standard horizons in the regional seismic profile structure interpretation graph.
The term explanation and detailed description in this step may refer to step S11 in the above embodiment, and will not be repeated here. Referring to fig. 3, an explanatory diagram of a structure of a seismic section of a region provided in this embodiment is shown, and two region standard horizons are selected according to the preset selection rule in this step, in this embodiment, a section of a tri-stack whisker-river group (T 3x1) and a section of a tri-stack feijian-guan group (T 1f 4) are selected, and in this embodiment, the bottom of the stratum horizon in the explanatory diagram is described.
And S22, respectively carrying out horizon leveling treatment on the selected regional standard horizons to obtain a leveled regional seismic section structure interpretation graph.
The term explanation and detailed description in this step may refer to step S12 in the above embodiment, and will not be repeated here. The horizon leveling in this embodiment is developed in layer sequence stratigraphy and geophysical prospecting technology, and assuming that the original thickness of each layer sequence is unchanged, taking a deposition reference plane or a maximum flooding plane as a reference, selecting a reference top surface and a bottom surface of a comparison layer sequence, subtracting the top surface time from the bottom surface time, namely leveling the top surface, and regarding the leveled surface as a lake plane during paleo-deposition to obtain the morphology of the bottom surface, namely the relative paleo-topography before layer sequence stratum deposition.
Specifically, the horizon leveling method of the seismic data is to correct a certain layer to a reference plane at any moment after explaining the layer, and all reflections above and below the layer are corrected with corresponding time. The effect of structural deformation can be removed by performing the horizon flattening process. Horizon flattening can be divided into horizon flattening sections and horizon flattening slices, which are all information such as amplitude, frequency, phase, velocity sections, slices and the like which are picked up and displayed at certain time intervals in a given time window, seismic data and interpreted construction data, and the sections which can be obtained after flattening processing are horizon paleo-geomorphic restoration sections. The seismic section after the paleo-topography recovery of a certain horizon is equivalent to the form of recovering the horizon when the horizon is deposited, namely the paleo-topography recovery is embodied. The contact relation and the structural development history of each structural layer can be studied by utilizing the horizon paleo-topography recovery profile.
Referring to fig. 4, in this step, the selected section (T 3x1) of the triax whisker river set is used as a reference, and the horizon leveling process is performed on the structural evolution profile, so as to obtain a structural explanatory diagram of the leveled regional seismic profile, and it can be seen from fig. 4 that the underlying horizon is partially recessed into a V shape, which is due to the incorrect structure caused by the presence of a steep structural band in the overlying stratum, and is practically absent.
And S23, carrying out vectorization processing on the at least two flattened regional seismic section structure explanatory diagrams to obtain at least two flattened regional seismic section vector diagrams.
Referring to fig. 5, the flattened regional seismic profile structural interpretation map obtained in step S22 is vectorized, and a flattened regional seismic profile vector map shown in the map is obtained. And carrying out vectorization processing on the bitmap image, and finding out the bottom of the stratum layer to represent the trend and evolution process of the stratum, so that the subsequent recovery processing is more convenient.
And step S24, determining at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon in the at least two leveled regional seismic section vector diagrams.
Referring to fig. 5, the early (T 1f 4) and late (T 3x1) flattened regional seismic profile vector diagrams are compared to determine at least one corresponding reference horizon, which in this embodiment is the next two-stack mountain section (P 2l1). From fig. 5, it can be seen that in the early and late leveled vector diagrams, the formation horizon yields below P 2l1 are very different, whereas in the geological history knowledge the yields should be approximately the same. The reason why the occurrence of the underburden is greatly changed is that the section of the overburden is greatly influenced by the structural band, and a great amount of false structural phenomena are generated when the standard horizons in the early (T 1f 4) area are leveled. Whereas the reference horizon (P 2l1) determined in this embodiment is not particularly dramatic in early and late leveled cross-sectional views (which may be understood herein as the vector diagram shown in FIG. 5) due to ease of tracking, the reference horizon (P 2l1) may be used to reduce or eliminate errors in late formation recovery from early formation.
In a more specific embodiment, determining the distance difference corresponding to the reference horizon, as shown with reference to fig. 7, may include the following two steps:
Determining a plurality of coordinate points representing curve forms corresponding to the reference layers in the regional seismic section vector diagram corresponding to the advanced regional standard layers;
And taking the vertical distance from each coordinate point to a curve corresponding to a reference horizon in the regional seismic profile vector diagram corresponding to the early regional standard horizon in the vertical direction as a distance difference value corresponding to the reference horizon.
The distance difference value corresponding to the reference horizon is the distance difference value between the reference horizon determined by the later-stage regional standard horizon and the corresponding reference horizon determined by the earlier-stage regional standard horizon in the leveled regional seismic section vector diagram.
As shown in connection with fig. 7, it is assumed that the distance between each pair of vertical correspondence points (o 1 at the upper late reference horizon coordinate point, its XY coordinates being Xo 1、Yo1, o 2 at the lower early reference horizon coordinate point, its XY coordinates being Xo 2、Yo2) is H, including H1, … … H500, … …, H1000, … …, hn. Let the point of the vertical upward shift of the early reference horizon o2 be P, xp, yp be the X, Y coordinates of the point P, then:
Xp＝Xo1 (1)
Yp＝Yo2+H (2)
Therefore Xp and Yp are coordinate points of vertical upward movement of the lower early reference layer, and n coordinate points are upward movement, so that a new upward movement layer can be formed by connecting, and the new upward movement layer should be overlapped with the upper late reference layer.
By analogy, the coordinates of each point of the earlier reference layer for each horizon below should be (1) and (2), and the coordinates of each horizon up-shift should be connected to form a new horizon, as shown in fig. 8.
And S25, based on the distance difference value corresponding to the reference horizon, moving the lower horizon of the reference horizon upwards by the distance difference value corresponding to the reference horizon in the regional seismic profile vector diagram corresponding to the advanced regional standard horizon.
Step S26, overlapping the shifted early reference layer and the lower layer thereof with the upper layer of the later reference layer.
Referring to fig. 9, old horizons below reference horizon P 2l1 in the late phase are removed, replaced by new horizons, horizons above reference horizon P 2l1 are unchanged, a new late-stage structure evolution profile can be formed, and the false structure phenomenon is significantly reduced compared with fig. 4 and 5.
And step S27, filling the superimposed vector diagram to generate a superimposed basin geological evolution section.
In order to reduce the influence of the phenomenon on the construction recovery in the embodiment of the invention, the early horizons in the construction section are vertically and upwards displaced through the reference horizons and then replaced to the corresponding horizons after the later horizons are leveled. Therefore, the influence of late structural deformation can be reduced by multiple times of replacement, and the authenticity of structural evolution is ensured.
The method obviously reduces the influence of the late structure on the early structure recovery in the superimposed section, largely eliminates the phenomena of false structure and false stratum thickness in the structure evolution, and recovers the basic outline and evolution rule of the regional structure in the geological history as much as possible. Provides a detailed data foundation for understanding the oil gas aggregation rule and evaluating the favorable zone in oil gas exploration.
Based on the same inventive concept, the embodiment of the invention further provides a device for restoring the overlapping basin area structure, and referring to fig. 11, the device may include: the working principle of the selecting module 11, the leveling module 12, the determining module 13 and the recovering module 14 is as follows:
The selection module 11 is configured to select at least two regional standard horizons in the regional seismic profile construction interpretation map. Specifically, the selection module 11 selects at least two regional standard horizons according to a preset standard; the preset criteria include at least one of: the method comprises the steps of widely distributed horizons, horizons with gentle paleolandform corresponding to geological periods, horizons with strong seismic reflection surface amplitude and horizons with extended regions.
The leveling module 12 is configured to perform horizon leveling processing on the selected regional standard horizons, so as to obtain a flattened regional seismic profile structure interpretation map.
The determining module 13 is configured to determine at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon according to the at least two flattened interpretation graphs of the seismic profile structure of the region. Specifically, the determining module 13 performs vectorization processing on at least two flattened regional seismic section structure explanatory diagrams to obtain at least two flattened regional seismic section vector diagrams; determining at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon in the at least two leveled regional seismic section vector diagrams; the distance difference value corresponding to the reference horizon is the distance difference value between the reference horizon determined by the later-stage regional standard horizon and the corresponding reference horizon determined by the earlier-stage regional standard horizon in the leveled regional seismic section vector diagram.
More specifically, the determining module 13 determines the distance difference value corresponding to the reference horizon, which may include: determining a plurality of coordinate points representing curve forms corresponding to the reference layers in the regional seismic section vector diagram corresponding to the advanced regional standard layers; and taking the vertical distance from each coordinate point to a curve corresponding to a reference horizon in the regional seismic profile vector diagram corresponding to the early regional standard horizon in the vertical direction as a distance difference value corresponding to the reference horizon.
The restoration module 14 is configured to restore the overlapping basin area structure according to at least one reference horizon and a distance difference value corresponding to the reference horizon. Specifically, the recovery module 14 ascends the distance difference value corresponding to the reference horizon from the lower horizon of the reference horizon in the area seismic profile vector diagram corresponding to the advanced area standard horizon based on the distance difference value corresponding to the reference horizon; and overlapping the moved early reference layer and the lower layer thereof with the upper layer of the later reference layer.
Of course, after the restoration module 14 performs restoration processing, the superimposed vector image is subjected to filling processing to generate a superimposed basin geological evolution section.
The specific manner in which the respective modules perform the operations in relation to the folded basin area construction restoration device in the above-described embodiments has been described in detail in relation to the embodiments of the method, and will not be described in detail herein.
Based on the same inventive concept, the embodiment of the invention also provides a superimposed basin region construction restoration system, which can comprise: a first subsystem and a second subsystem;
The first subsystem is used for selecting at least two regional standard horizons from the regional seismic profile construction interpretation map; respectively leveling the selected regional standard horizons to obtain a leveled regional seismic section structure interpretation map; the second subsystem is used for determining at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon according to the at least two flattened regional seismic section structure interpretation graphs; and recovering the overlapped basin area structure according to the at least one reference horizon and the distance difference value corresponding to the reference horizon.
Further, the first subsystem may be a Linux system, and the second subsystem may be a Windows system. In the above-described different data processing, different software modules may be used, for example, a seismic workstation based on the Linux system, and a vectoring software module Coreldraw based on the Windows system, which is not particularly limited in this embodiment. It should be noted that, the first subsystem and the second subsystem may be installed on different servers, or may be installed on the same server, and the above scheme in this embodiment may be implemented only by performing a system switching when executing the above steps.
The specific manner in which the operations of the various subsystems are performed in relation to the folded basin area construction restoration system of the above-described embodiments has been described in detail in relation to the method embodiments and will not be described in detail herein.
Based on the same inventive concept, the embodiment of the invention also provides an application of the overlapped basin geological evolution section obtained by the overlapped basin region structure recovery method in geological modeling.
Based on the same inventive concept, the embodiment of the invention also provides a geological modeling system which comprises the overlapping basin region structure recovery device.
In the technical field of petroleum development, oil reservoir geological modeling and numerical simulation are two important means for researching oil reservoir exploration and development, and play an important role in oil reservoir exploration and development. The oil deposit geological modeling is to build a geological model of the oil deposit, and the numerical simulation is to solve an oil deposit mathematical model by using a computer, simulate underground oil-water flow and give oil-water distribution at a certain moment so as to predict the oil deposit dynamics. The superimposed basin address evolution section obtained in the embodiment has great significance for geological modeling, and provides a new direction for cognition of oil and gas aggregation rules in oil and gas exploration.
The method for restoring the overlapped basin area structure in the above embodiment is applied to the geologic modeling and the geologic modeling system, wherein the specific manner of performing the operations by each part is described in detail in the embodiment related to the method, and the implementation can be referred to the implementation of the foregoing method, and the repetition is omitted.
It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, magnetic disk storage, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It will be apparent to those skilled in the art that various modifications and variations can be made to the present invention without departing from the spirit or scope of the invention. Thus, it is intended that the present invention also include such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.
Claims (8)
1. A method for restoring the structure of a region of a superimposed basin, comprising:
Selecting at least two regional standard horizons from the regional seismic profile construction interpretation map;
Respectively carrying out horizon flattening treatment on the selected regional standard horizons to obtain a flattened regional seismic profile structure interpretation map;
vectorizing at least two flattened regional seismic section structure explanatory diagrams to obtain at least two flattened regional seismic section vector diagrams;
Determining at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon in at least two flattened regional seismic section vector diagrams; the distance difference value corresponding to the reference horizon is the distance difference value between the reference horizon determined by the later-stage regional standard horizon and the corresponding reference horizon determined by the earlier-stage regional standard horizon in the leveled regional seismic section vector diagram;
Shifting the lower layer of the reference layer up by the distance difference value corresponding to the reference layer in the regional seismic profile vector diagram corresponding to the advanced regional standard layer based on the distance difference value corresponding to the reference layer;
overlapping the moved early reference horizon and the underlying horizon with the overlying layer of the late reference horizon to restore the overlapping basin area structure.
2. The method of claim 1, wherein determining a distance difference value for the reference horizon comprises:
Determining a plurality of coordinate points representing curve forms corresponding to the reference horizon in the regional seismic section vector diagram corresponding to the advanced regional standard horizon;
And taking the vertical distance from each coordinate point to the curve corresponding to the reference horizon in the regional seismic profile vector diagram corresponding to the early regional standard horizon in the vertical direction as the distance difference value corresponding to the reference horizon.
3. The method of claim 1, further comprising, after the laminating process: and filling the overlapped vector images to generate the overlapped basin geological evolution section.
4. A method according to any one of claims 1 to 3, wherein selecting at least two regional standard horizons in the regional seismic profile interpretation comprises: selecting at least two regional standard horizons according to a preset standard;
The preset criteria include at least one of: the method comprises the steps of widely distributed horizons, horizons with gentle paleolandform corresponding to geological periods, horizons with strong seismic reflection surface amplitude and horizons with extended regions.
5. A superimposed basin area construction restoration device, comprising:
the selecting module is used for selecting at least two regional standard horizons from the regional seismic profile structure interpretation graph;
the leveling module is used for respectively leveling the selected regional standard horizons to obtain a leveled regional seismic section structure interpretation graph;
The determining module is used for carrying out vectorization processing on at least two flattened regional seismic section structure explanatory diagrams to obtain at least two flattened regional seismic section vector diagrams; determining at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon in at least two flattened regional seismic section vector diagrams; the distance difference value corresponding to the reference horizon is the distance difference value between the reference horizon determined by the later-stage regional standard horizon and the corresponding reference horizon determined by the earlier-stage regional standard horizon in the leveled regional seismic section vector diagram;
The recovery module is used for moving the lower horizon of the reference horizon up by the distance difference value corresponding to the reference horizon in the regional seismic profile vector diagram corresponding to the advanced regional standard horizon based on the distance difference value corresponding to the reference horizon; overlapping the moved early reference horizon and the underlying horizon with the overlying layer of the late reference horizon to restore the overlapping basin area structure.
6. A superimposed basin area construction restoration system, comprising: a first subsystem and a second subsystem;
the first subsystem is used for selecting at least two regional standard horizons from the regional seismic profile construction interpretation map; respectively carrying out horizon flattening treatment on the selected regional standard horizons to obtain flattened regional seismic section structure interpretation diagrams;
the second subsystem is used for carrying out vectorization processing on at least two flattened regional seismic section structure explanatory diagrams to obtain at least two flattened regional seismic section vector diagrams; determining at least one corresponding reference horizon and a distance difference value corresponding to the reference horizon in at least two flattened regional seismic section vector diagrams; based on the distance difference value corresponding to the reference horizon, moving the lower horizon of the reference horizon upwards by the distance difference value corresponding to the reference horizon in the regional seismic profile vector diagram corresponding to the advanced regional standard horizon; overlapping the moved early reference horizon and the lower horizon thereof with the upper cladding of the later reference horizon to restore the overlapped basin area structure;
and the distance difference value corresponding to the reference horizon is the distance difference value between the reference horizon determined by the later-stage regional standard horizon and the corresponding reference horizon determined by the earlier-stage regional standard horizon in the leveled regional seismic section vector diagram.
7. Use of a superimposed basin geological evolution profile obtained by a superimposed basin region structure restoration method according to any of the claims 1-4 in geological modeling.
8. A geologic modeling system comprising the superimposed basin area structure restoration device of claim 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110696240.XA CN115508884B (en) | 2021-06-23 | 2021-06-23 | Method, device and system for restoring overlapped basin area structure and application thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110696240.XA CN115508884B (en) | 2021-06-23 | 2021-06-23 | Method, device and system for restoring overlapped basin area structure and application thereof |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115508884A CN115508884A (en) | 2022-12-23 |
| CN115508884B true CN115508884B (en) | 2024-06-07 |
Family
ID=84499605
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110696240.XA Active CN115508884B (en) | 2021-06-23 | 2021-06-23 | Method, device and system for restoring overlapped basin area structure and application thereof |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115508884B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102200588A (en) * | 2011-03-22 | 2011-09-28 | 成都理工大学 | Method for analyzing waveform similarity body curvature of seismic data |
| CN107015290A (en) * | 2017-03-13 | 2017-08-04 | 西北大学 | A kind of method that reworked garden basin primary deposit looks are recovered |
-
2021
- 2021-06-23 CN CN202110696240.XA patent/CN115508884B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102200588A (en) * | 2011-03-22 | 2011-09-28 | 成都理工大学 | Method for analyzing waveform similarity body curvature of seismic data |
| CN107015290A (en) * | 2017-03-13 | 2017-08-04 | 西北大学 | A kind of method that reworked garden basin primary deposit looks are recovered |
Non-Patent Citations (6)
| Title |
|---|
| New Application of Horizon Flattening Technology with Seismic Data for Exploration and Development in Fula Sub-basin, Muglad Basin, Sudan;Aixiang Liu 等;Proceedings of the International Field Exploration and Development Conference 2020;20210615;第3081–3087页 * |
| 一种古构造恢复方法探讨――以乌马营构造为例;张立勤 等;矿物岩石;20050630;第25卷(第02期);第93-98页 * |
| 双界面地震层拉平的古地貌恢复技术及应用 — 以鄂尔多斯盆地天环坳陷为例;刘永涛 等;石油地球物理;20190630;第54卷(第3期);第656-666页 * |
| 叠合盆地差异构造变形与油气聚集;汤良杰 等;地学前缘;20090731;第16卷(第04期);第13-22页 * |
| 层拉平技术在地震解释中的深化应用;郭雯 等;石油地球物理勘探;20201031;第55卷(第05期);第1110-1120页 * |
| 滨里海盆地盐下构造假象识别及真实形态恢复;王波 等;中国石油勘探;20080229(第01期);第63-68页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115508884A (en) | 2022-12-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3329307B1 (en) | Assignment of systems tracts | |
| US10995592B2 (en) | Method and system for analyzing the uncertainty of subsurface model | |
| Van Wees et al. | Geomechanics response and induced seismicity during gas field depletion in the Netherlands | |
| Radwan | Three-dimensional gas property geological modeling and simulation | |
| CN109870719B (en) | Well position layout method, device and system for carbonate rock tight thin reservoir | |
| US11506807B2 (en) | Methods and systems for simulation gridding with partial faults | |
| EP2288944A2 (en) | Virtual petroleum system | |
| WO2018080745A1 (en) | Method and system for enhancing meshes for a subsurface model | |
| WO2009148706A1 (en) | 3d visualization of 2d geophysical data | |
| Tyler et al. | Modeling heterogeneities in fluvial domains: a review of the influence on production profiles | |
| US11409024B2 (en) | Methods and systems for generating simulation grids via zone by zone mapping from design space | |
| US11163093B2 (en) | Constructing stratigraphic images of sediment distribution in a subsurface of the earth | |
| Lobatskaya et al. | GIS-based analysis of fault patterns in urban areas: a case study of Irkutsk city, Russia | |
| CN112130210B (en) | Well position optimization method and device based on ancient and modern geological elements | |
| Alpak et al. | Dynamic impact and flow-based upscaling of the estuarine point-bar stratigraphic architecture | |
| CN115508884B (en) | Method, device and system for restoring overlapped basin area structure and application thereof | |
| Zhang et al. | Turbidite facies identification and analysis using geophysical methods in the salt diapir domain of Lower Congo basin | |
| CA2724117A1 (en) | Litho-facies interpolation | |
| CN113985476B (en) | Stratum ancient structure recovery method, device, equipment and readable storage medium | |
| Shodmonovna | Structural 3D Modeling and Creation of a ThreeDimensional Grid | |
| Weijermans et al. | Integrated Reservoir Characterization and Simulation to Optimize Recovery from a Mature Carboniferous North Sea Gas Field with Water Influx | |
| Potapova et al. | Integrated approach to the correlation polyfacies objects deposits | |
| Rashad et al. | An iterative workflow introduced to model complex heterogeneous reservoirs: an example from NEAG-2 Field, north Western Desert, Egypt | |
| Meirovitz | From Rocks to Models: Use of Outcrop Characterization for Improved Reservoir Performance Prediction | |
| Ma et al. | Constructing 3D Model Framework and Change of Support in Data Mapping |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant |