CN112965142B - Method, device, equipment and storage medium for generating karst ancient political view - Google Patents
Method, device, equipment and storage medium for generating karst ancient political view Download PDFInfo
- Publication number
- CN112965142B CN112965142B CN202110174230.XA CN202110174230A CN112965142B CN 112965142 B CN112965142 B CN 112965142B CN 202110174230 A CN202110174230 A CN 202110174230A CN 112965142 B CN112965142 B CN 112965142B
- Authority
- CN
- China
- Prior art keywords
- karst
- information
- value
- slope
- target layer
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Processing Or Creating Images (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The application provides a method, a device, equipment and a storage medium for generating a karst ancient apparent picture. According to the method, geological data information is obtained firstly, and according to the geological data information, a target layer in a seismic data work area and depth information of each construction unit in the target layer are determined. And calculating the normalized elevation value of the karst basin and the normalized elevation value of the karst slope based on the depth information of each construction unit in the target layer, and finally generating the karst ancient apparent picture of the target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope. According to the technical scheme, the depth information of the karst basin, the karst slope and the karst elevation in the target layer is processed, and the purpose of more accurately generating the karst ancient apparent map is achieved.
Description
Technical Field
The application relates to the technical field of exploration and development, in particular to a method, a device, equipment and a storage medium for generating a karst ancient scenic map.
Background
The karst ancient landform is a three-dimensional relief form with height change formed by erosion and other geological transformation in a certain geological history period, and the karst ancient landform directly controls the development of a karst reservoir stratum, so that the karst ancient landform is essential for research on recovery of the karst ancient landform.
In the prior art, there are many common methods for restoring karst ancient landforms, taking a residual thickness method as an example, specifically, after the landform to be restored is degraded and begins to be deposited on an overlying stratum as an equal-time surface, a certain special layer section in the stratum is selected as a reference surface and leveled, and the obtained residual thickness on the reference surface is distributed into a restored ancient landform shape.
However, in practical applications, in tilted monoclinic formations caused by later tectonic movements, only the tilted side is subjected to weathering and erosion for a long time, and the geologic structure is characterized in that the formations are tilted, and epibiotic karsts are not necessarily developed on the target layer, so that the ancient geomorphology recovered by the method is not accurate.
Disclosure of Invention
The application provides a method, a device, equipment and a storage medium for generating a karst ancient apparent picture. The method is used for solving the problem that the prior art can not accurately recover the karst ancient landform morphology.
In a first aspect, an embodiment of the present application provides a method for generating a karst paleo-apparent map, including:
acquiring geological data information;
according to the geological data information, determining a target layer in a seismic data work area and depth information of each construction unit in the target layer, wherein the construction unit of the target layer comprises a karst basin, a karst plateau and a karst slope located between the karst basin and the karst plateau;
calculating a normalized elevation value of the karst basin and a normalized elevation value of the karst slope based on the depth information of each construction unit in the target layer;
and generating the karst ancient apparent picture of the target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope.
In one possible design of the first aspect, the determining, from the geological data information, a target layer in a seismic data work area and depth information of each structural unit in the target layer includes:
constructing a seismic data work area according to the three-dimensional pre-stack time migration data volume and the seismic interpretation data in the geological data information;
determining a denudation tip vanishing point and an upper overtopping point of a target layer in the seismic data work area;
carrying out plane position calibration according to the ablation tip vanishing point and the upper overtaking point, and determining the target layer;
and calculating the depth information of each construction unit in the target layer according to the drilling and logging data and the layered data in the geological data information.
In this possible design, the calculating depth information of each structural unit in the target layer according to the drilling log data and the layered data in the geological data information includes:
analyzing the drilling and logging data and the layered data in the geological data information to determine the current residual thickness value distribution information of each construction unit in the target layer;
determining depth information of the karst basin, depth information of the karst plateau, original residual thickness information of the karst basin and current residual thickness information of the karst slope according to current residual thickness value distribution information of each construction unit in the target layer;
determining the denuded thickness information of each statistical point in the karst slope according to the original residual thickness information of each statistical point in the karst basin and the current residual thickness information of each statistical point in the karst slope;
and determining the depth information of the karst slope according to the information of the denuded thickness at each statistical point in the karst slope, the minimum depth value of the denuded point and the maximum depth value of the upper overtaking point.
Optionally, the determining the depth information of the karst slope according to the information of the eroded thickness at each statistical point in the karst slope, the minimum depth value of the point where the erosion is quenched and the maximum depth value of the upper overtaking point includes:
determining the maximum value and the minimum value of the denudated thickness of the karst slope according to the denudated thickness information of each statistical point in the karst slope;
determining a maximum difference value of the eroded thickness of the karst slope according to the maximum value of the eroded thickness and the minimum value of the eroded thickness;
determining the maximum depth difference value of the karst slope according to the minimum depth value of the point where the erosion is extinguished and the maximum depth value of the upper overtaking point;
and determining the depth information of the karst slope according to the denudated thickness information, the maximum denudated thickness difference and the maximum depth difference at each statistical point in the karst slope.
In another possible design of the first aspect, the determining the normalized elevation value of the karst basin and the normalized elevation value of the karst slope based on the depth information of each construction unit in the target layer includes:
determining the maximum depth value of the karst plateau and the minimum depth value of the karst basin according to the depth information of each construction unit in the target layer;
determining the current elevation difference of the target layer according to the maximum depth value of the karst elevation and the minimum depth value of the karst basin;
calculating a normalized elevation value of the karst basin according to the depth information of each statistic point in the karst basin and the current elevation difference by taking the current elevation difference as a reference for normalization processing;
and calculating the normalized elevation value of the karst slope according to the depth information of the karst slope and the current elevation difference by taking the current elevation difference as a reference for normalization processing.
In yet another possible design of the first aspect, the generating the karst paleo-apparent map of the target layer according to the normalized elevation values of the karst basin and the normalized elevation values of the karst slope includes:
determining karst paleotopographic change information of the target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope;
identifying ancient water system distribution information in the karst ancient landform according to the karst ancient landform change information, the water system development characteristics and configuration information among the landform characteristics;
and generating the karst ancient political view based on the karst ancient geomorphic change information and the ancient water system distribution information.
In a second aspect, an embodiment of the present application provides a device for generating a karst ancient geomorphology map, including: the device comprises an acquisition module, a determination module, a processing module and a generation module;
the acquisition module is used for acquiring geological data information;
the determining module is used for determining a target layer in a seismic data work area and depth information of each construction unit in the target layer according to the geological data information, wherein the construction unit of the target layer comprises a karst basin, a karst plateau and a karst slope located between the karst basin and the karst plateau;
the processing module is used for calculating the normalized elevation value of the karst basin and the normalized elevation value of the karst slope based on the depth information of each construction unit in the target layer;
and the generation module is used for generating the karst ancient apparent map of the target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope.
In a possible design of the second aspect, the determining module is specifically configured to:
constructing a seismic data work area according to the three-dimensional pre-stack time migration data volume and the seismic interpretation data in the geological data information;
determining a denudation tip vanishing point and an upper overtopping point of a target layer in the seismic data work area;
carrying out plane position calibration according to the ablation tip vanishing point and the upper overtaking point, and determining the target layer;
and calculating the depth information of each construction unit in the target layer according to the drilling and logging data and the layered data in the geological data information.
In this possible design, the determining module is configured to calculate depth information of each structural unit in the target layer according to the drilling log data and the layered data in the geological data information, and specifically includes:
the determining module is specifically configured to:
analyzing the drilling and logging data and the layered data in the geological data information to determine the current residual thickness value distribution information of each construction unit in the target layer;
determining depth information of the karst basin, depth information of the karst plateau, original residual thickness information of the karst basin and current residual thickness information of the karst slope according to current residual thickness value distribution information of each construction unit in the target layer;
determining the denuded thickness information of each statistical point in the karst slope according to the original residual thickness information of each statistical point in the karst basin and the current residual thickness information of each statistical point in the karst slope;
and determining the depth information of the karst slope according to the information of the denuded thickness at each statistical point in the karst slope, the minimum depth value of the denuded point and the maximum depth value of the upper overtaking point.
Optionally, the determining module is configured to determine the depth information of the karst slope according to the eroded thickness information at each statistical point in the karst slope, the minimum depth value of the point where the erosion is quenched and extinguished, and the maximum depth value of the upper overtaking point, and specifically includes:
the determining module is specifically configured to:
determining the maximum value and the minimum value of the denudated thickness of the karst slope according to the denudated thickness information of each statistical point in the karst slope;
determining a maximum difference value of the eroded thickness of the karst slope according to the maximum value of the eroded thickness and the minimum value of the eroded thickness;
determining the maximum depth difference value of the karst slope according to the minimum depth value of the point where the erosion is extinguished and the maximum depth value of the upper overtaking point;
and determining the depth information of the karst slope according to the denudated thickness information, the maximum denudated thickness difference and the maximum depth difference at each statistical point in the karst slope.
In another possible design of the second aspect, the processing module is specifically configured to:
determining the maximum depth value of the karst plateau and the minimum depth value of the karst basin according to the depth information of each construction unit in the target layer;
determining the current elevation difference of the target layer according to the maximum depth value of the karst elevation and the minimum depth value of the karst basin;
calculating a normalized elevation value of the karst basin according to the depth information of each statistic point in the karst basin and the current elevation difference by taking the current elevation difference as a reference for normalization processing;
and calculating the normalized elevation value of the karst slope according to the depth information of the karst slope and the current elevation difference by taking the current elevation difference as a reference for normalization processing.
In yet another possible design of the second aspect, the generating module is specifically configured to:
determining karst paleotopographic change information of the target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope;
identifying ancient water system distribution information in the karst ancient landform according to the karst ancient landform change information, the water system development characteristics and configuration information among the landform characteristics;
and generating the karst ancient political view based on the karst ancient geomorphic change information and the ancient water system distribution information.
In a third aspect, an embodiment of the present application provides a computer device, including: a processor, a memory, a display, and a computer program stored on the memory and executable on the processor;
the processor, when executing the computer program instructions, implements the first aspect and the method of generating a karst ancient appearance map provided by each possible design.
In a fourth aspect, embodiments of the present application provide a computer-readable storage medium having stored therein computer program instructions for implementing the method for generating a karst paleo-apparent map provided by the first aspect and various possible designs when executed by a processor.
In a fifth aspect, embodiments of the present application provide a computer program product comprising a computer program which, when executed by a processor, is configured to implement the method for generating the karst paleo-apparent map provided by the first aspect and each possible design.
The embodiment of the application provides a method, a device, equipment and a storage medium for generating a karst ancient scenic map. According to the method, geological data information is obtained, and according to the geological data information, a target layer in a seismic data work area and depth information of each construction unit in the target layer are determined. And calculating the normalized elevation value of the karst basin and the normalized elevation value of the karst slope based on the depth information of each construction unit in the target layer, and finally generating the karst ancient apparent picture of the target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope. According to the technical scheme, the depth information of the karst basin, the karst slope and the karst elevation in the target layer is processed, so that the karst ancient apparent map is generated more accurately.
Drawings
FIG. 1 is a flow chart of a first embodiment of a method for generating a karst paleo-apparent map provided in an embodiment of the present application;
fig. 2A is a flowchart of a second method for generating a karst paleography provided in an embodiment of the present application;
FIG. 2B is a schematic diagram of a target layer of an ancient apparent karst map provided by an embodiment of the application;
FIG. 2C is a present elevation difference view of a target layer of an ancient apparent view of a karst provided by an embodiment of the present application;
FIG. 3 is a flow chart of a third embodiment of a method for generating a karst paleo-apparent map provided in an embodiment of the present application;
fig. 4 is a flowchart of a fourth embodiment of a method for generating a karst paleography provided in an embodiment of the present application;
FIG. 5 is a general flow chart of an embodiment of a method for generating a karst paleography provided by an embodiment of the present application;
fig. 6 is a schematic structural diagram of a device for generating a karst paleography according to an embodiment of the present application;
fig. 7 is a schematic structural diagram of a computer device according to an embodiment of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are some embodiments of the present application, but not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
Before introducing the embodiments of the present application, a description will be given of a background and an application scenario of the present application.
Karst paleography is a three-dimensional topographic relief form with high and low changes formed by erosion and other geological changes in a certain geological historical period. The karst ancient landform directly controls the development of the karst reservoir, so the karst ancient landform is essential to research on the karst ancient landform.
The common karst ancient landform restoration method comprises the following steps: residual thickness methods, impression methods, sequence stratigraphy methods, sedimentology methods, horizon leveling methods, etc., but each of these methods has some geological condition pertinence in application. For example, the residual thickness method is to deposit an overburden stratum into an isochronal surface after the degradation of the landform to be restored is finished, select a specific interval in the stratum as a reference surface, flatten the reference surface, and distribute the residual thickness on the surface into an ancient landform shape; the method comprises the steps of selecting an isochronous reference surface from an overburden stratum deposited after the weathering and ablation of a stratum to be restored is finished, and restoring the ancient landform form shape semi-quantitatively by utilizing the mirror image relationship between the stratum deposited at the later stage and an ablation surface; the stratigraphic sequence stratigraphic framework method is characterized in that after a regional isochronal plane is selected to be leveled by building a stratigraphic sequence stratigraphic framework of an overlying stratum, layers starting from the bottom surface of a denudation plane are connected by a smooth curve to represent an ancient landform form; the sedimentology method is to understand the ancient terrain characteristics of each region through geological mapping and the denudation degree of each region, analyze sedimentary facies and ancient environment, and qualitatively recover the ancient landform; the method for leveling the horizon is characterized in that the original thicknesses of the sequence are assumed to be unchanged, the reference top and the reference bottom of the comparison sequence are selected in the three-dimensional seismic data volume by taking a reference surface or a maximum flooding surface as a reference, the bottom surface and the top surface are subtracted, namely, the top surface is leveled, the leveled surface is taken as a lake plane during ancient deposition, and the obtained bottom surface form is a relative ancient landform before deposition.
However, under some complicated geological conditions, the above methods have certain limitations, and cannot accurately restore the ancient landform. If the tilted monoclinic stratum caused by later tectonic movement only one tilted side is subjected to weathering and corrosion for a long time, the geologic structure is characterized in that the stratum is tilted, and the epibiotic karst is not necessarily developed on the target layer, so that the ancient landform recovered by the method is not accurate.
In order to solve the technical problems, the technical conception process of the inventor is as follows: the inventor finds that by taking the karst ancient landform with the half moat structure as the substrate as an example, if the karst ancient landform is divided into a karst basin, a karst plateau and a karst slope by a certain relation, and the three parts and the associated information between the three parts are analyzed and processed by combining geological data information, a more accurate karst ancient apparent image map compared with the prior art can be obtained.
The technical solution of the present application will be described in detail below with reference to specific examples. It should be noted that the following specific embodiments may be combined with each other, and the same or similar concepts or processes may not be described in detail in some embodiments.
Fig. 1 is a flowchart of a first embodiment of a method for generating a karst paleo-apparent graph according to an embodiment of the present application. As shown in fig. 1, the method for generating the karst paleography may include the following steps:
and 11, acquiring geological data information.
In this embodiment, the restoration of the karst ancient landform with the half moat structure as the base is described as an example.
In this step, in order to generate the karst ancient landform map, geological data information related to the karst ancient landform needs to be collected, acquired and sorted.
Optionally, the geological data information may include: geological background data, well drilling rock cores and slices, well logging lithology data, well logging data, layering data, a three-dimensional pre-stack time migration data volume, matched seismic interpretation data (horizon and fault interpretation) and the like.
And step 12, determining a target layer in the seismic data work area and depth information of each construction unit in the target layer according to the geological data information.
In the step, a seismic data work area is constructed according to the acquired geological data information, and the seismic data work area is used as an implementation basis for generating a karst ancient and ancient apparent map.
In one possible implementation, the seismic data work area may be an off-the-shelf work area collected during collaboration with the oil field, such as a Petrel work area.
Optionally, the region to be recovered, i.e., the target zone, may be obtained by processing the seismic data work area.
Specifically, the target layer may be composed of each structural unit, and depth information of each structural unit may be determined according to geological data information.
Wherein the construction unit of the target layer comprises a karst basin, a karst plateau and a karst slope between the karst basin and the karst plateau.
And step 13, calculating the normalized elevation value of the karst basin and the normalized elevation value of the karst slope based on the depth information of each construction unit in the target layer.
In this step, the construction units in the target layer are a karst basin, a karst plateau and a karst slope, respectively. According to the depth information of the karst basin, the karst elevation and the karst slope, the normalized elevation value of the karst basin and the normalized elevation value of the karst slope are determined, and the normalized elevation value of the karst basin and the normalized elevation value of the karst slope can be used as main parameters for generating the karst ancient apparent map.
In one possible implementation, the maximum depth value of the karst elevation and the minimum depth value of the karst basin may be determined according to the determined depth information of each construction unit, and then the current elevation difference of the target layer may be determined based on the maximum depth value of the karst elevation and the minimum depth value of the karst basin. Calculating a normalized elevation value of the karst basin according to the depth information of each statistic point in the karst basin and the current elevation difference by taking the current elevation difference as a reference for normalization processing; and calculating the normalized elevation value of the karst slope according to the depth information of the karst slope and the current elevation difference.
And 14, generating an ancient apparent karst map of the target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope.
In the step, when the normalized elevation value of the karst basin and the normalized elevation value of the karst slope are determined, the elevation data corresponding to the two normalized elevation values are arranged and input into a preset model, so that a three-dimensional landform change map of the karst ancient landform can be generated, corresponding ancient water system distribution information is configured on the three-dimensional landform change map, and the karst ancient landform map of the target layer can be obtained through certain drawing elements of modification, coloring and complete response.
In one possible design, the model may be a software Petrel.
According to the method for generating the karst ancient and political map, geological data information is obtained, and according to the geological data information, the depth information of a target layer in a seismic data work area and each construction unit in the target layer is determined. And calculating the normalized elevation value of the karst basin and the normalized elevation value of the karst slope based on the depth information of each construction unit in the target layer, and finally generating the karst ancient apparent picture of the target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope. According to the technical scheme, the depth information of the target layer and each construction unit in the target layer is determined, the depth information is processed, the karst ancient political view is generated more accurately, and the error problem existing in the prior art when the karst ancient geomorphology with the semi-through structure as the substrate is recovered is solved.
On the basis of fig. 1, fig. 2A is a flowchart of a second method for generating a karst paleography provided in the embodiment of the present application. As shown in fig. 2A, the step 12 can be implemented by:
and step 21, constructing a seismic data work area according to the three-dimensional pre-stack time migration data body and the seismic interpretation data in the geological data information.
In this scenario, fig. 2B is a schematic diagram of a target layer of an ancient political view of a karst provided in an embodiment of the present application; fig. 2C is a present elevation difference view of a karst ancient apparent map of a target layer provided by an embodiment of the present application. The embodiments of the present application will be described with reference to fig. 2B and 2C.
In this step, the seismic data work area may be constructed by obtaining a three-dimensional prestack time migration data volume and seismic interpretation data in the geological data information.
In one possible implementation, the seismic data work area may be an off-the-shelf work area collected through a collaborative process with an oil field, such as a Petrel work area. The Petrel work area is constructed by a three-dimensional prestack time migration data body and seismic interpretation data in geological data information.
And step 22, determining a denudation tip vanishing point and an upper overtaking point of the target layer in the seismic data work area.
In this step, two key location points are identified one by one on the survey line in each direction within the seismic data work area. The two location points are respectively a denuded tip vanishing point and an upper overtaking point (i.e. an upper overtaking point on the formation).
Optionally, the ablation tip vanishing point is point P1 in fig. 2B, and the upper overtopping point is point P2 in fig. 2B. Wherein. Wherein, the P1 point and the P2 point are points divided by the ablation boundary in the seismic data work area.
And step 23, calibrating the plane position according to the ablation tip vanishing point and the upper overtaking point, and determining the target layer.
In this step, a region to be restored, i.e., a target layer, is determined based on the ablation tip vanishing point and the upper overtaking point.
In one possible implementation, the ablation tip vanishing point and the upper overtopping point are plane-position-calibrated, i.e. the P1 point and the P2 point are connected by a smooth curve, and the construction units are distinguished on the plane.
Specifically, the part except the point P1 can be directly divided into a karst bulge (karst plateau), namely an exposed ablation surface; outside the point P2 are epigenetic karst weakly developed or undeveloped regions (karst basins); the section between point P1 and point P2 is the area of incomplete erosion (the karst slope).
The construction unit of the target layer is composed of a karst plateau, a karst basin and a karst slope.
And 24, calculating the depth information of each construction unit in the target layer according to the drilling and logging data and the layered data in the geological data information.
In this step, based on the determined target zone and the drilling and logging data and the layered data in the geological data information, the depth information of the karst plateau, the karst basin and the karst slope is determined.
Specifically, this step can be implemented as follows:
and step 1, analyzing drilling and logging data and layered data in the geological data information to determine current residual thickness value distribution information of each construction unit in a target layer.
Wherein the current residual thickness value distribution information of each construction unit includes: current residual thickness information of the karst basin, current residual thickness information of the karst plateau, and current residual thickness information of the karst slope.
And step 2, determining the depth information of the karst basin, the depth information of the karst plateau, the original residual thickness information of the karst basin and the current residual thickness information of the karst slope according to the current residual thickness value distribution information of each construction unit in the target layer.
Optionally, the original residual thickness information of the karst basin is determined according to the current residual thickness value distribution information of the karst basin, as shown in fig. 2B, the original residual thickness information of the karst basin may be obtained in the following manner:
in particular, well bores (i.e., statistical points) in undecorated areas of the well log data and the layered data may be collected, with well bores located in the karst basin, and with well bores of a number that may be presentIs n, HiThe original residual thickness for the ith well, then HaveThe calculation formula for the average of the original residual thickness of the karst basin is as follows (ignoring the wells not located in the karst basin):
Have＝(H1+H2……+Hn)/n
optionally, the current residual thickness information of the karst slope is determined according to the current residual thickness value distribution information of the karst slope, as shown in fig. 2B, the current residual thickness information of the karst slope may be obtained in the following manner:
specifically, well bores in the area suffering from degradation in the well bore log data and the layered data may be collected, the well bores may be located in a karst slope, and the number of well bores may be n, hiThe current residual thickness for the ith well is expressed as follows (ignoring wells not in the karst slope):
hi＝h1，h2，……hn
and 3, determining the denuded thickness information of each statistical point in the karst slope according to the original residual thickness information of each statistical point in the karst basin and the current residual thickness information of each statistical point in the karst slope.
Alternatively, the average of the original residual thickness of the karst basin may be taken as the original residual thickness at each statistical point in the karst slope.
In one possible implementation, the denuded thickness at each statistical point in the karst slope may be obtained by subtracting the current residual thickness at each statistical point in the karst slope from the average of the original residual thicknesses of the karst basin, respectively.
Specifically, the denuded thickness at each statistical point in the karst slope may be BiIs shown, in which:
Bi＝Have–hi，i＝1,2,……n
and 4, determining the depth information of the karst slope according to the denuded thickness information of each statistical point in the karst slope, the minimum depth value of the denuded point and the maximum depth value of the upper overtaking point.
Optionally, a denudation tip vanishing point and an upper overtopping point may be determined on each survey line of the target layer, and then a minimum depth value is determined from the denudation tip vanishing points, and a maximum depth value is determined from the upper overtopping points. Then, the depth information of the karst slope is determined by utilizing the denuded thickness information, the minimum depth value of the denuded point and the maximum depth value of the upper overtaking point, and the method can be realized by the following specific method:
firstly, determining the maximum value and the minimum value of the denuded thickness of the karst slope according to the denuded thickness information of each statistical point in the karst slope.
In particular, the thickness B of the denuded zone at a statistical point of the karst slopeiIn the method, the maximum value of the denudation thickness of the karst slope is obtained and recorded as Bimax(ii) a Minimum value of denuded thickness, noted Bimin。
And secondly, determining the maximum difference of the denuded thickness of the karst slope according to the maximum denuded thickness and the minimum denuded thickness.
In particular, the maximum difference in ablated thickness of the karst slope may be denoted as Δ B, where Δ B ═ Bimax-Bimin。
And thirdly, determining the maximum depth difference value of the karst slope according to the minimum depth value of the point where the erosion is quenched and the maximum depth value of the upper overtaking point.
Specifically, the maximum depth difference of the karst slope may be represented by Δ P, and the minimum depth value of the point where the erosion is peaked off may be represented by P1minThe maximum depth value of the upper overtoint may be represented by P2maxAnd (4) showing. The calculation formula can be:
ΔP＝P2max–P1min
and finally, determining the depth information of the karst slope according to the denudated thickness information, the maximum denudated thickness difference value and the maximum depth difference value of each statistical point in the karst slope.
In one possible design, the depth information of the karst slope may be a depth range value for which the denudation thickness is converted into the karst slope, using Bi-newAnd (4) showing. It calculatesThe formula can be:
Bi-new＝(Bi/ΔB)*ΔP
according to the method for generating the karst ancient scenic map, the seismic data work area is constructed according to the three-dimensional prestack time migration data body and the seismic interpretation data in the geological data information. And then determining a denudation tip vanishing point and an upper overtaking point of the target layer in the seismic data work area, and calibrating the plane position based on the denudation tip vanishing point and the upper overtaking point to determine the target layer. And calculating the depth information of each construction unit in the target layer according to the drilling and logging data and the layered data in the geological data information. According to the technical scheme, the depth information of the target layer and each construction unit in the target layer is determined, and a foundation is provided for generating the paleo-karst apparent map.
On the basis of the above embodiment, fig. 3 is a flowchart of a third embodiment of a method for generating a karst ancient political view provided in the embodiment of the present application. As shown in fig. 3, the step 13 can be implemented by:
and step 31, determining the maximum depth value of the karst plateau and the minimum depth value of the karst basin according to the depth information of each construction unit in the target layer.
Optionally, from the depth information of the karst basins, the thickness value corresponding to the lowest point of the karst basin in the target layer, i.e. the minimum depth value of the karst basin, may be determined, as shown in fig. 2C, where the lowest point of the karst basin is denoted by G1. In one possible implementation, the minimum thickness value S (G1) of the karst basin may be 0, with the horizontal line at G1 as the horizontal reference line.
Optionally, according to the depth information of the karst elevation, a thickness value corresponding to the highest point of the karst elevation in the target layer, that is, the maximum depth value of the karst elevation, may be determined, as shown in fig. 2C, where the highest point of the karst elevation is represented by G2.
And step 32, determining the current elevation difference of the target layer according to the maximum depth value of the karst elevation and the minimum depth value of the karst basin.
Alternatively, the current elevation difference of the target layer may be determined by subtracting the minimum depth value of the karst basin from the maximum depth value of the karst elevation,current elevation difference GnowThe calculation formula of (a) is as follows:
Gnow＝S(G2)-S(G1)
wherein, S (G2) may be the altitude where the G2 point is located, i.e., the maximum depth value of the karst plateau, and S (G1) may be the altitude where the G1 point is located, i.e., the minimum depth value of the karst basin.
And step 33, taking the current elevation difference as a normalization treatment as a reference, and calculating the normalized elevation value of the karst basin according to the depth information of each statistic point in the karst basin and the current elevation difference.
In this step, the current elevation difference of the target layer is used as a reference for normalization, specifically, S (G2) is 1, S (G1) is 0, and the depth information G of each statistical point in the karst basin is madeiVarying between 0 and 1.
In one possible implementation, G isnowPerforming normalization processing to obtain normalized elevation value G of karst basini-newThe calculation formula of (a) is as follows:
Gi-new＝(Gi–Gnow)*100％
and step 34, taking the current elevation difference as a normalization treatment as a reference, and calculating the normalized elevation value of the karst slope according to the depth information of the karst slope and the current elevation difference.
In this step, G isnowPerforming normalization processing based on the depth information B of the karst slopei-newAnd the current elevation difference GnowDetermining a normalized elevation value B of the karst slopei-newThe calculation formula is as follows:
*Bi-new＝(Bi-new-Gnow)*100％
according to the method for generating the karst ancient field image, the maximum depth value of the karst plateau and the minimum depth value of the karst basin are determined according to the depth information of each construction unit in the target layer. And then determining the current elevation difference of the target layer according to the maximum depth value of the karst elevation and the minimum depth value of the karst basin. Finally, on the basis of normalization processing of the current elevation difference, calculating a normalized elevation value of the karst basin according to the depth information of each statistic point in the karst basin and the current elevation difference; and calculating the normalized elevation value of the karst slope according to the depth information of the karst slope and the current elevation difference. According to the technical scheme, the normalized elevation value of the karst slope and the normalized elevation value of the karst basin are determined through the depth information of each construction unit, and data information is provided for the subsequent generation of the karst ancient and ancient apparent map.
On the basis of the foregoing embodiment, fig. 4 is a flowchart of a fourth embodiment of a method for generating a karst paleography provided in the embodiment of the present application. As shown in fig. 4, the step 14 can be implemented by:
and step 41, determining the karst paleotopographic change information of the target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope.
In the step, the data corresponding to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope are input into Petrel software after being sorted, the boundary range of the seismic data work area is defined, and a three-dimensional landform change diagram can be displayed in a 3D window by utilizing triangular contour interpolation.
Wherein, the karst ancient landform change information of the target layer can be the three-dimensional landform change map.
And 42, identifying ancient water system distribution information in the karst ancient landform according to the karst ancient landform change information, the water system development characteristics and the configuration information among the landform characteristics.
In the above steps, only the three-dimensional landform change map is acquired, and the ancient water system distribution information needs to be added to the three-dimensional landform change map to perfect and generate the subsequent karst ancient apparent map.
In one possible implementation, possible ancient water system distribution characteristics can be identified on the recovered three-dimensional ancient landform change map through the concept of 'ancient future theory' according to the current landform and the corresponding water system development characteristics. And then, qualitatively identifying main ancient water system distribution information according to the ancient landform and topography characteristics, the distribution of karst valleys, valley lands and karst depressions on the plane and the mutual configuration relationship.
Further, the identified ancient water system distribution information is added to the three-dimensional landform change map.
And 43, generating a karst ancient political view map based on the karst ancient landform change information and the ancient water system distribution information.
In this step, the three-dimensional landform change map carrying the ancient water system distribution information is further refined.
In a possible implementation, a three-dimensional geomorphic change map carrying ancient water system distribution information is modified and colored, and each ancient geomorphic unit is marked by combining drilling and logging data and relative elevation change, so that corresponding drawing elements are perfected, and a final karst ancient geomorphology map is obtained.
According to the method for generating the karst ancient landform map, the karst ancient landform change information of the target layer is determined according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope. And then identifying ancient water system distribution information in the karst ancient landform according to the karst ancient landform change information, the water system development characteristics and the configuration information among the landform characteristics. And finally, generating the karst ancient landform map based on the karst ancient landform change information and the ancient water system distribution information. In the technical scheme, the three-dimensional landform change map is determined based on the normalized elevation value of the karst basin and the normalized elevation value of the karst slope, ancient water system distribution information is added and perfected on the three-dimensional landform change map, and the karst ancient land apparent map can be generated more accurately.
On the basis of the above embodiments, fig. 5 is a general flow chart of an embodiment of a method for generating a karst ancient appearance map provided in the embodiments of the present application. The technical solutions related to the embodiments of the present application are briefly summarized, and details of the implementation principle are not described herein. As shown in fig. 5, the general flow chart includes the following steps:
step 1, acquiring geological data information;
step 2, determining a target layer and a construction unit corresponding to the target layer according to geological data information;
step 3, calculating a normalized elevation value of a karst basin and a normalized elevation value of a karst slope in the construction unit;
step 4, generating a three-dimensional landform change diagram relative to terrain change by utilizing the normalized elevation value of the karst basin and the normalized elevation value of the karst slope;
step 5, adding ancient water system distribution information to the three-dimensional landform change diagram;
and 6, perfecting the three-dimensional landform change map added with the ancient water system distribution information to generate a karst ancient political map.
According to the method for generating the karst ancient scenic map, the target layer and the construction unit corresponding to the target layer are determined according to the obtained geological data information. And then, calculating the normalized elevation value of the karst basin and the normalized elevation value of the karst slope in the construction unit, and generating a three-dimensional landform change map relative to terrain change by using the normalized elevation value of the karst basin and the normalized elevation value of the karst slope. And finally, adding ancient water system distribution information to the three-dimensional landform change map, and finally modifying the three-dimensional landform change map added with the ancient water system distribution information, so that the karst ancient landform map taking the half-moat structure as the substrate is generated more accurately.
On the basis of the foregoing method for generating a karst paleo-topography map, fig. 6 is a schematic structural diagram of a device for generating a karst paleo-topography map provided in the embodiment of the present application. As shown in fig. 6, the generating means includes: an acquisition module 61, a determination module 62, a processing module 63 and a generation module 64.
An obtaining module 61, configured to obtain geological data information;
the determining module 62 is configured to determine a target layer in the seismic data work area and depth information of each construction unit in the target layer according to the geological data information, where the construction unit of the target layer includes a karst basin, a karst plateau, and a karst slope between the karst basin and the karst plateau;
the processing module 63 is configured to calculate a normalized elevation value of the karst basin and a normalized elevation value of the karst slope based on the depth information of each construction unit in the target layer;
and the generating module 64 is used for generating the karst ancient apparent picture of the target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope.
In one possible design of the embodiment of the present application, the determining module 62 is specifically configured to:
constructing a seismic data work area according to a three-dimensional prestack time migration data volume and seismic interpretation data in geological data information;
determining a denudation tip vanishing point and an upper overtopping point of a target layer in a seismic data work area;
carrying out plane position calibration according to the ablation tip vanishing point and the upper overtaking point to determine a target layer;
and calculating the depth information of each construction unit in the target layer according to the drilling and logging data and the layered data in the geological data information.
In this possible design, the determining module 62 is configured to calculate depth information of each structural unit in the target layer according to the drilling log data and the layered data in the geological data information, specifically:
the determining module 62 is specifically configured to:
analyzing drilling and logging data and layered data in the geological data information to determine current residual thickness value distribution information of each construction unit in a target layer;
determining depth information of a karst basin, depth information of a karst plateau, original residual thickness information of the karst basin and current residual thickness information of a karst slope according to current residual thickness value distribution information of each construction unit in a target layer;
determining the denuded thickness information of each statistical point in the karst slope according to the original residual thickness information of each statistical point in the karst basin and the current residual thickness information of each statistical point in the karst slope;
and determining the depth information of the karst slope according to the denuded thickness information of each statistical point in the karst slope, the minimum depth value of the denuded point vanishing point and the maximum depth value of the upper overtaking point.
Optionally, the determining module 62 is configured to determine the depth information of the karst slope according to the eroded thickness information at each statistical point in the karst slope, the minimum depth value of the point where the erosion is quenched and eroded, and the maximum depth value of the upper overtaking point, where:
the determining module 62 is specifically configured to:
determining the maximum value and the minimum value of the denudated thickness of the karst slope according to the denudated thickness information of each statistical point in the karst slope;
determining the maximum difference value of the denuded thickness of the karst slope according to the maximum value of the denuded thickness and the minimum value of the denuded thickness;
determining the maximum depth difference value of the karst slope according to the minimum depth value of the point where the erosion is extinguished and the maximum depth value of the upper overtaking point;
and determining the depth information of the karst slope according to the denudated thickness information, the maximum denudated thickness difference value and the maximum depth difference value at each statistical point in the karst slope.
In another possible design of the embodiment of the present application, the processing module 63 is specifically configured to:
determining the maximum depth value of the karst plateau and the minimum depth value of the karst basin according to the depth information of each construction unit in the target layer;
determining the current elevation difference of the target layer according to the maximum depth value of the karst elevation and the minimum depth value of the karst basin;
calculating a normalized elevation value of the karst basin according to the depth information of each statistic point in the karst basin and the current elevation difference by taking the current elevation difference as a reference for normalization processing;
and taking the current elevation difference as a reference for normalization processing, and calculating the normalized elevation value of the karst slope according to the depth information of the karst slope and the current elevation difference.
In yet another possible design of the embodiment of the present application, the generating module 64 is specifically configured to:
determining karst ancient landform change information of a target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope;
identifying ancient water system distribution information in the karst ancient landform according to the karst ancient landform change information, the water system development characteristics and the configuration information among the landform characteristics;
and generating the karst ancient landform map based on the karst ancient landform change information and the ancient water system distribution information.
The device for generating the karst ancient topographic map provided by the embodiment of the application can be used for executing the technical scheme for generating the karst ancient topographic map in the embodiment, the implementation principle and the technical effect are similar, and the description is omitted.
It should be noted that the division of the modules of the above apparatus is only a logical division, and the actual implementation may be wholly or partially integrated into one physical entity, or may be physically separated. And these modules can be realized in the form of software called by processing element; or may be implemented entirely in hardware; and part of the modules can be realized in the form of calling software by the processing element, and part of the modules can be realized in the form of hardware. For example, the processing module 63 may be a separate processing element, or may be integrated into a chip of the apparatus, or may be stored in a memory of the apparatus in the form of program code, and a processing element of the apparatus calls and executes the functions of the processing module 63. Other modules are implemented similarly. In addition, all or part of the modules can be integrated together or can be independently realized. The processing element here may be an integrated circuit with signal processing capabilities. In implementation, each step of the above method or each module above may be implemented by an integrated logic circuit of hardware in a processor element or an instruction in the form of software.
Fig. 7 is a schematic structural diagram of a computer device according to an embodiment of the present application. As shown in fig. 7, the computer apparatus may include: a processor 71, a memory 72 and a display 73.
Processor 71 executes computer-executable instructions stored in memory, causing processor 71 to perform the aspects of the embodiments described above. The processor 71 may be a general-purpose processor including a central processing unit CPU, a Network Processor (NP), and the like; but also a digital signal processor DSP, an application specific integrated circuit ASIC, a field programmable gate array FPGA or other programmable logic device, discrete gate or transistor logic, discrete hardware components.
A memory 72 is coupled to the processor 71 via the system bus and communicates with each other, the memory 72 storing computer program instructions.
The display 73 may be a user interface for displaying the karst ancient landscape. The user interface may include graphics, text, icons, video, and any combination thereof.
The system bus may be a Peripheral Component Interconnect (PCI) bus, an Extended Industry Standard Architecture (EISA) bus, or the like. The system bus may be divided into an address bus, a data bus, a control bus, and the like. For ease of illustration, only one thick line is shown, but this does not mean that there is only one bus or one type of bus. The transceiver is used to enable communication between the database access device and other computers (e.g., clients, read-write libraries, and read-only libraries). The memory may comprise Random Access Memory (RAM) and may also include non-volatile memory (non-volatile memory), such as at least one disk memory.
The computer device provided in the embodiment of the present application may be used to execute the technical solution of the karst ancient scenic map generation method in the above embodiment, and the implementation principle and the technical effect are similar, which are not described herein again.
The embodiment of the application further provides a chip for running the instructions, and the chip is used for executing the technical scheme of the karst ancient political view generation method in the embodiment.
The embodiment of the present application further provides a computer-readable storage medium, where computer instructions are stored in the computer-readable storage medium, and when the computer instructions are run on a computer, the computer is caused to execute the technical solution of the karst paleo-map generation method according to the above embodiment.
The embodiments of the present application further provide a computer program product, where the computer program product includes a computer program, where the computer program is stored in a computer-readable storage medium, and at least one processor may read the computer program from the computer-readable storage medium, where the at least one processor, when executing the computer program, may implement the technical solution of the method for generating the karst ancient scenic map in the above embodiments.
Finally, it should be noted that: the above embodiments are only used for illustrating the technical solutions of the present application, and not for limiting the same; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present application.
Claims (9)
1. A method for generating a karst ancient political view is characterized by comprising the following steps:
acquiring geological data information;
according to the geological data information, determining a target layer in a seismic data work area and depth information of each construction unit in the target layer, wherein the construction unit of the target layer comprises a karst basin, a karst plateau and a karst slope located between the karst basin and the karst plateau;
calculating a normalized elevation value of the karst basin and a normalized elevation value of the karst slope based on the depth information of each construction unit in the target layer;
generating an ancient apparent karst map of the karst of the target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope;
generating the karst ancient apparent map of the target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope, wherein the generating comprises the following steps:
determining karst paleotopographic change information of the target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope;
identifying ancient water system distribution information in the karst ancient landform according to the karst ancient landform change information, the water system development characteristics and configuration information among the landform characteristics;
and generating the karst ancient political view based on the karst ancient geomorphic change information and the ancient water system distribution information.
2. The method of claim 1, wherein determining a target interval in a seismic data work area and depth information for each constitutional unit in the target interval based on the geological data information comprises:
constructing a seismic data work area according to the three-dimensional pre-stack time migration data volume and the seismic interpretation data in the geological data information;
determining a denudation tip vanishing point and an upper overtopping point of a target layer in the seismic data work area;
carrying out plane position calibration according to the ablation tip vanishing point and the upper overtaking point, and determining the target layer;
and calculating the depth information of each construction unit in the target layer according to the drilling and logging data and the layered data in the geological data information.
3. The method of claim 2, wherein calculating depth information for each of the construction units in the target formation based on the drilling log data and the stratigraphic data in the geological data information comprises:
analyzing the drilling and logging data and the layered data in the geological data information to determine the current residual thickness value distribution information of each construction unit in the target layer;
determining depth information of the karst basin, depth information of the karst plateau, original residual thickness information of the karst basin and current residual thickness information of the karst slope according to current residual thickness value distribution information of each construction unit in the target layer;
determining the denuded thickness information of each statistical point in the karst slope according to the original residual thickness information of each statistical point in the karst basin and the current residual thickness information of each statistical point in the karst slope;
and determining the depth information of the karst slope according to the information of the denuded thickness at each statistical point in the karst slope, the minimum depth value of the denuded point and the maximum depth value of the upper overtaking point.
4. The method of claim 3, wherein determining the depth information of the karst slope according to the eroded thickness information at each statistical point in the karst slope, the point minimum depth value of the eroded pinches, and the maximum depth value of the upper overtop includes:
determining the maximum value and the minimum value of the denudated thickness of the karst slope according to the denudated thickness information of each statistical point in the karst slope;
determining a maximum difference value of the eroded thickness of the karst slope according to the maximum value of the eroded thickness and the minimum value of the eroded thickness;
determining the maximum depth difference value of the karst slope according to the minimum depth value of the point where the erosion is extinguished and the maximum depth value of the upper overtaking point;
and determining the depth information of the karst slope according to the denudated thickness information, the maximum denudated thickness difference and the maximum depth difference at each statistical point in the karst slope.
5. The method of any one of claims 1-4, wherein determining the normalized elevation value for the karst basin and the normalized elevation value for the karst slope based on the depth information for each of the construction units in the target layer comprises:
determining the maximum depth value of the karst plateau and the minimum depth value of the karst basin according to the depth information of each construction unit in the target layer;
determining the current elevation difference of the target layer according to the maximum depth value of the karst elevation and the minimum depth value of the karst basin;
calculating a normalized elevation value of the karst basin according to the depth information of each statistic point in the karst basin and the current elevation difference by taking the current elevation difference as a reference for normalization processing;
and calculating the normalized elevation value of the karst slope according to the depth information of the karst slope and the current elevation difference by taking the current elevation difference as a reference for normalization processing.
6. An apparatus for generating a karst ancient political view, comprising: the device comprises an acquisition module, a determination module, a processing module and a generation module;
the acquisition module is used for acquiring geological data information;
the determining module is used for determining a target layer in a seismic data work area and depth information of each construction unit in the target layer according to the geological data information, wherein the construction unit of the target layer comprises a karst basin, a karst plateau and a karst slope located between the karst basin and the karst plateau;
the processing module is used for calculating the normalized elevation value of the karst basin and the normalized elevation value of the karst slope based on the depth information of each construction unit in the target layer;
the generation module is used for generating a karst paleo-apparent map of the target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope;
the generation module is specifically configured to:
determining karst paleotopographic change information of the target layer according to the normalized elevation value of the karst basin and the normalized elevation value of the karst slope;
identifying ancient water system distribution information in the karst ancient landform according to the karst ancient landform change information, the water system development characteristics and configuration information among the landform characteristics;
and generating the karst ancient political view based on the karst ancient geomorphic change information and the ancient water system distribution information.
7. A computer device comprising a processor, a memory, a display and a computer program stored on the memory and executable on the processor, the processor when executing the program implementing a method of generating a karst historic site map according to any one of claims 1 to 5.
8. A computer-readable storage medium having stored thereon computer-executable instructions for implementing a method of generating a karst paleo-map as claimed in any one of claims 1 to 5 when executed by a processor.
9. A computer program product comprising a computer program for implementing a method of generating a karst paleo-apparent map as claimed in any one of claims 1 to 5 when executed by a processor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110174230.XA CN112965142B (en) | 2021-02-09 | 2021-02-09 | Method, device, equipment and storage medium for generating karst ancient political view |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110174230.XA CN112965142B (en) | 2021-02-09 | 2021-02-09 | Method, device, equipment and storage medium for generating karst ancient political view |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112965142A CN112965142A (en) | 2021-06-15 |
| CN112965142B true CN112965142B (en) | 2022-04-12 |
Family
ID=76284275
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110174230.XA Active CN112965142B (en) | 2021-02-09 | 2021-02-09 | Method, device, equipment and storage medium for generating karst ancient political view |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112965142B (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105372716A (en) * | 2015-10-28 | 2016-03-02 | 中国石油大学(华东) | Evaluation method for distribution of supergene karst reservoir of carbonatite |
| CN105510993A (en) * | 2015-12-23 | 2016-04-20 | 中国石油大学(北京) | Foreland basin deep buried and compressed type complex gypsum-salt rock identification and distribution prediction method |
| CN105956238A (en) * | 2016-04-22 | 2016-09-21 | 中国石油天然气股份有限公司 | Palaeogeomorphic restoration method and apparatus |
| CN107589470A (en) * | 2017-08-24 | 2018-01-16 | 西南石油大学 | Palaeokarst Landform restoration methods in small area |
| CN108267783A (en) * | 2018-01-02 | 2018-07-10 | 中国石油天然气集团公司 | A kind of method, apparatus and system of determining buried-hill trap |
| CN108734779A (en) * | 2017-04-17 | 2018-11-02 | 中国石油化工股份有限公司 | A kind of palaeogeomorphology |
| CN110824557A (en) * | 2019-10-22 | 2020-02-21 | 中国石油天然气股份有限公司 | Quantitative karst ancient landform recovery method based on well-seismic combination |
| CN111239818A (en) * | 2020-02-12 | 2020-06-05 | 成都理工大学 | Ancient landform analysis method based on three-dimensional dip angle attribute body correction |
| CN111650642A (en) * | 2020-06-08 | 2020-09-11 | 中国石油大学(华东) | Industrialized process for recovering relatively micro ancient landform by using three-dimensional earthquake |
| CN112305602A (en) * | 2019-08-01 | 2021-02-02 | 中国石油天然气股份有限公司 | Carbonate reservoir prediction method based on prestack multi-attribute and ancient landform fusion technology |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105425316A (en) * | 2015-10-28 | 2016-03-23 | 中国石油化工股份有限公司 | Ancient landform restoration method based on impression method |
| CN108196302A (en) * | 2017-11-28 | 2018-06-22 | 中国石油天然气股份有限公司 | Dolomite fracture hole develops reservoir "sweet spot" area's Forecasting Methodology and device |
-
2021
- 2021-02-09 CN CN202110174230.XA patent/CN112965142B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105372716A (en) * | 2015-10-28 | 2016-03-02 | 中国石油大学(华东) | Evaluation method for distribution of supergene karst reservoir of carbonatite |
| CN105510993A (en) * | 2015-12-23 | 2016-04-20 | 中国石油大学(北京) | Foreland basin deep buried and compressed type complex gypsum-salt rock identification and distribution prediction method |
| CN105956238A (en) * | 2016-04-22 | 2016-09-21 | 中国石油天然气股份有限公司 | Palaeogeomorphic restoration method and apparatus |
| CN108734779A (en) * | 2017-04-17 | 2018-11-02 | 中国石油化工股份有限公司 | A kind of palaeogeomorphology |
| CN107589470A (en) * | 2017-08-24 | 2018-01-16 | 西南石油大学 | Palaeokarst Landform restoration methods in small area |
| CN108267783A (en) * | 2018-01-02 | 2018-07-10 | 中国石油天然气集团公司 | A kind of method, apparatus and system of determining buried-hill trap |
| CN112305602A (en) * | 2019-08-01 | 2021-02-02 | 中国石油天然气股份有限公司 | Carbonate reservoir prediction method based on prestack multi-attribute and ancient landform fusion technology |
| CN110824557A (en) * | 2019-10-22 | 2020-02-21 | 中国石油天然气股份有限公司 | Quantitative karst ancient landform recovery method based on well-seismic combination |
| CN111239818A (en) * | 2020-02-12 | 2020-06-05 | 成都理工大学 | Ancient landform analysis method based on three-dimensional dip angle attribute body correction |
| CN111650642A (en) * | 2020-06-08 | 2020-09-11 | 中国石油大学(华东) | Industrialized process for recovering relatively micro ancient landform by using three-dimensional earthquake |
Non-Patent Citations (2)
| Title |
|---|
| "川东南—黔西北桐湾Ⅲ幕岩溶古地貌恢复";冯伟明 等;《地质论评》;20170930;第63卷(第5期);第1270-1280页 * |
| "鄂尔多斯盆地靖边气田的古地貌定量恢复新方法";乔博 等;《天然气勘探与开发》;20181231;第41卷(第4期);第32-37页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112965142A (en) | 2021-06-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Natali et al. | Modeling Terrains and Subsurface Geology. | |
| EP3217192B1 (en) | Method for looking for oil and gas reservoir on the basis of trap-3d software | |
| CA2912626C (en) | 3d trap evaluation method of searching for oil-gas reservoir | |
| CN105593908A (en) | Identifying matching properties between a group of bodies representing a geological structure and a table of properties | |
| CA2724115A1 (en) | 3d visualization of 2d geophysical data | |
| US10387583B2 (en) | Rotations from gradient directions | |
| CN109598068B (en) | Ancient structure constraint modeling method, device and equipment | |
| CN112965142B (en) | Method, device, equipment and storage medium for generating karst ancient political view | |
| Sahoo et al. | Creating three-dimensional channel bodies in LiDAR-integrated outcrop characterization: A new approach for improved stratigraphic analysis | |
| CN111983671A (en) | Shallow water lake basin reservoir prediction method and device based on micro-ancient landform restoration | |
| US11385369B2 (en) | Three-dimensional, stratigraphically-consistent seismic attributes | |
| CN112800518B (en) | Stratum surface model correction method utilizing adjacent stratum cross-correlation constraint | |
| CA2724117A1 (en) | Litho-facies interpolation | |
| CN113240771B (en) | Automatic generation method and device of stratum isopachous map | |
| US20160116618A1 (en) | Method for determining a stacked model describing architectural elements | |
| CN113269879B (en) | Automatic construction method and device for three-dimensional model of geologic body | |
| CN112346117B (en) | Reservoir characteristic prediction method and device based on seismic attribute fusion | |
| CN113269852B (en) | Stratum partition map and stratum isobath mapping method and device | |
| CN113495303B (en) | Method and device for recovering denudation quantity | |
| CN112666607B (en) | Method and device for inverting thickness distribution of yellow soil layer by gravity | |
| CN114076982A (en) | Karst ancient landform restoration method and device based on waveform characteristic difference | |
| CN116047595A (en) | Method for describing micro-amplitude structure | |
| CN117950034A (en) | Balance profile acquisition method and device | |
| CN112904423A (en) | Method and device for carving continental river | |
| CN114994748A (en) | Ancient landform restoration method and device |
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 | ||
| GR01 | Patent grant |