CN113640870B - Method, system and equipment for recovering sedimentary paleo-topography of karst ablation interface - Google Patents
Method, system and equipment for recovering sedimentary paleo-topography of karst ablation interface Download PDFInfo
- Publication number
- CN113640870B CN113640870B CN202110851311.9A CN202110851311A CN113640870B CN 113640870 B CN113640870 B CN 113640870B CN 202110851311 A CN202110851311 A CN 202110851311A CN 113640870 B CN113640870 B CN 113640870B
- Authority
- CN
- China
- Prior art keywords
- thickness
- stratum
- paleo
- deposition
- residual
- 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
- 238000012876 topography Methods 0.000 title claims abstract description 58
- 238000002679 ablation Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000008021 deposition Effects 0.000 claims abstract description 81
- 238000011084 recovery Methods 0.000 claims abstract description 24
- 230000015572 biosynthetic process Effects 0.000 claims description 52
- 238000010606 normalization Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 13
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 7
- 239000004215 Carbon black (E152) Substances 0.000 abstract description 2
- 229930195733 hydrocarbon Natural products 0.000 abstract description 2
- 150000002430 hydrocarbons Chemical class 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 8
- 230000008859 change Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 208000035126 Facies Diseases 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 235000017166 Bambusa arundinacea Nutrition 0.000 description 1
- 235000017491 Bambusa tulda Nutrition 0.000 description 1
- 241001330002 Bambuseae Species 0.000 description 1
- 235000015334 Phyllostachys viridis Nutrition 0.000 description 1
- 238000005299 abrasion Methods 0.000 description 1
- 239000011425 bamboo Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005137 deposition process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 238000009304 pastoral farming Methods 0.000 description 1
- 238000002407 reforming Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
-
- 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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a method, a system and equipment for recovering paleo-topography deposited on a karst ablation interface, wherein the method comprises the steps of obtaining paleo-topography recovery data of a target karst area, wherein the paleo-topography recovery data comprises impression stratum thickness and residual stratum thickness; acquiring the relative original stratum thickness according to the impression stratum thickness and the residual stratum thickness; judging whether the residual stratum is affected by deposition, if so, correcting the thickness of the relative original stratum to obtain a deposition paleo-landform value; and recovering the deposition paleo-landform plane according to the deposition paleo-landform value. The invention can extract the deposition information expressed in the paleo-topography, thereby being convenient for obtaining a more comprehensive paleo-topography recovery result and better guiding the exploration and development work of the karst type hydrocarbon reservoir.
Description
Technical Field
The application relates to the technical field of paleo-topography restoration, in particular to a method, a system and equipment for restoring a deposited paleo-topography of a karst ablation interface.
Background
Karst weathered crust reservoirs in carbonate reservoirs are one of important reservoir types, and a large number of karst weathered crust reservoirs are found in areas such as Tarim basin, erdos basin and Sichuan basin, and large-scale oil and gas fields such as Tahe, jing edge, weifar and Anyue are found. In karst weathering crust reservoirs, ancient landforms are often important reservoir control factors, are closely related to the distribution of reservoirs, and have very important significance for the exploration and development of oil and gas.
At present, an impression method is mainly used in the ancient landform research of the reservoir, and is focused on researching the characteristics of the degraded ancient landform, namely the karst ancient landform characteristics, and the main research focus is on the transformation effect of the karst effect on the reservoir. However, in many carbonate reservoirs, the sedimentary facies also have important control effects on reservoirs, paleo-topography is subject to the results of the comprehensive effects of original sedimentary facies, post-reforming and the like, while the traditional paleo-topography research can ignore sedimentary information, and only using karst paleo-topography is likely to obtain one-sided results, so that deviation can be predicted, which is unfavorable for guiding the exploration and development of karst type reservoirs, therefore, the inventor considers that the paleo-topography recovery mode of the existing karst weathered crust reservoirs still needs further improvement.
Disclosure of Invention
In order to overcome the defects of the prior art, the application provides a method, a system and equipment for recovering the paleo-topography deposited on a karst ablation interface, which can extract deposition information (paleo-topography deposited) expressed in paleo-topography, thereby facilitating obtaining a more comprehensive paleo-topography recovery result and better guiding the exploration and development work of a karst type hydrocarbon reservoir.
In a first aspect, the present application provides a method for restoring a sedimentary paleo-topography of a karst ablation interface, the method comprising:
obtaining paleomorphic restoration data of a target karst region, wherein the paleomorphic restoration data comprises impression stratum thickness and residual stratum thickness;
acquiring the relative original stratum thickness according to the impression stratum thickness and the residual stratum thickness;
judging whether the residual stratum is affected by deposition, if so, correcting the thickness of the relative original stratum to obtain a deposition paleo-landform value;
and recovering the deposition paleo-landform plane according to the deposition paleo-landform value.
Optionally, before the step of obtaining the paleomorphic restoration data of the target karst region, the method further includes:
constructing an upper and lower recovery datum plane of an ablation interface of a target karst region based on stratum division, wherein the upper and lower recovery datum plane comprises an overlying datum plane and an underlying datum plane;
obtaining the thickness of the impression stratum above the ablation interface according to the depth of the overlying reference surface and the depth of the ablation interface;
and obtaining the thickness of the residual stratum below the degraded interface according to the depth of the degraded interface and the depth of the underlying reference surface.
Optionally, obtaining a relative original formation thickness according to the impression formation thickness and the residual formation thickness includes:
performing standardization treatment on the impression stratum thickness and the residual stratum thickness;
based on the normalized impression formation thickness and the residual formation thickness, the relative original formation thickness is calculated according to the following specific calculation formula:
wherein H is 0 Representing the relative original formation thickness; h r Representing the residual formation thickness after the normalization treatment; h i The thickness of the impression layer after normalization treatment is shown.
Optionally, determining whether the residual formation is affected by deposition includes:
acquiring a seismic profile of a target karst region, and analyzing the seismic profile according to preset rule conditions, wherein the preset rule conditions comprise: the form of the underlying reference surface is uneven, the underlying reference surface is degraded, and the carbonate bench deposition environment is the development of a hillside and beach phase and/or the overburden point of the residual stratum;
if the preset rule condition appears in the seismic profile, determining that the residual stratum is affected by deposition;
and if the preset rule condition does not appear in the seismic profile, determining that the residual stratum is not affected by deposition.
Optionally, after determining whether the residual formation is affected by deposition, the method further comprises: if yes, selecting a construction reference surface, wherein the construction reference surface is used for reflecting the relief form of the landform before the deposition of the residual stratum; and subtracting the depth corresponding to the construction reference plane from the depth corresponding to the underlying reference plane to obtain the underlying stratum thickness, and carrying out standardization treatment on the underlying stratum thickness.
Optionally, the sedimentary paleo-topography value is obtained after correcting the relative original stratum thickness, and the related calculation formula is as follows:
P＝H s *H 0
wherein P represents a deposition paleo-landform value, H s Representing the underburden thickness; h 0 Representing the relative original formation thickness.
Optionally, after determining whether the residual formation is affected by deposition, the method further comprises: if not, determining the relative original stratum thickness as a deposited paleo-physical value.
Optionally, the normalization processing adopts a Min-Max normalization method, and a specific calculation formula is as follows:
wherein H is nor Is a normalized value; h * Is the value to be normalized; h min Is the minimum value of the corresponding thickness; h max Is the maximum value of the corresponding thickness.
In a second aspect, the present application provides a deposited paleo-topography retrieval system for a karst ablation interface, the system comprising:
the data acquisition module is used for acquiring paleomorphic restoration data of the target karst area, wherein the paleomorphic restoration data comprises impression stratum thickness and residual stratum thickness;
the relative calculation module is used for obtaining the relative original stratum thickness according to the impression stratum thickness and the residual stratum thickness;
the judging module is used for judging whether the residual stratum is affected by deposition, and if so, correcting the thickness of the relative original stratum to obtain a deposition paleomorphic value;
and the recovery module is used for recovering the deposition paleo-landform plane according to the deposition paleo-landform value.
In a third aspect, the present application provides a computer device, which adopts the following technical scheme:
a computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor implementing the steps of a method of restoring a paleo-deposited topography of the karst ablation interface when the computer program is executed.
The application has the following beneficial technical effects: according to the impression thickness of the target karst region and the thickness of the residual stratum, the relative original stratum thickness is obtained, and the influence of differential ablation caused by the structure on the residual thickness can be recovered; when the residual stratum is affected by the sedimentary inheritance, the relative original stratum thickness is corrected to obtain a sedimentary paleo-topography value, and sedimentary phase change caused by the sedimentary inheritance is eliminated, so that the sedimentary paleo-topography of the target karst area is recovered conveniently, a more comprehensive paleo-topography recovery result is obtained, and the method is helpful for guiding the exploration and development of karst oil and gas reservoirs.
Drawings
FIG. 1 is a method flow diagram of one embodiment of a method for restoring a sedimentary paleo-topography of a karst ablation interface provided by the present invention;
FIG. 2 is a schematic drawing showing the selection of three reference planes at the ablation interface of a target karst region provided by the invention;
FIG. 3 is a functional block diagram of one embodiment of a deposited paleo-topography retrieval system of a karst ablation interface provided by the present invention.
Detailed Description
The present application is described in further detail below in conjunction with figures 1-3.
The embodiment of the application discloses a method for recovering a sedimentary paleo-topography of a karst denudation interface, referring to fig. 1, the method comprises the following steps:
s11: and obtaining paleomorphic restoration data of the target karst region, wherein the paleomorphic restoration data comprises a impression stratum thickness and a residual stratum thickness.
S12: based on the impression formation thickness and the residual formation thickness, a relative raw formation thickness is obtained.
S13: and judging whether the residual stratum is affected by deposition, and if so, correcting the thickness of the relative original stratum to obtain a deposition paleomorphic value.
S14: and recovering the deposition paleo-landform plane according to the deposition paleo-landform value.
According to the method, the device and the system, the relative original stratum thickness is obtained according to the impression thickness and the residual stratum thickness of the target karst area, and the influence of differential ablation caused by construction on the residual thickness can be recovered; when the residual stratum is affected by the sedimentary inheritance, the relative original stratum thickness is corrected to obtain a sedimentary paleo-topography value, and sedimentary phase change caused by the sedimentary inheritance is eliminated, so that the sedimentary paleo-topography of the target karst area is recovered conveniently, a more comprehensive paleo-topography recovery result is obtained, and the method is helpful for guiding the exploration and development of karst oil and gas reservoirs.
Before step S11 of this embodiment obtains paleo-topography recovery data of a target karst region, the method for recovering paleo-topography deposited on a karst ablation interface of this embodiment further includes the steps of:
s101: and constructing an upper and lower recovery datum plane of an ablation interface of the target karst region based on stratum division, wherein the upper and lower recovery datum plane comprises an overlying datum plane and an underlying datum plane.
S102: and obtaining the thickness of the impression stratum above the ablation interface according to the depth of the overlying reference surface and the depth of the ablation interface.
S103: and obtaining the thickness of the residual stratum below the degraded interface according to the depth of the degraded interface and the depth of the underlying reference surface.
In this embodiment, the upper and lower recovery reference surfaces include an overlying reference surface and a underlying reference surface, are both equal-time interfaces, and are easy to identify on a seismic section; further, the requirement of selecting an overlying datum plane is consistent with that of a traditional impression method; the requirements for selecting the underlying reference surface include: the stratum between the underlying datum plane and the ablation interface is a complete continuous deposition gyratory, and belongs to the same secondary gyratory; the distance between the underlying datum plane and the ablation interface is between one to two tertiary convolutions; the underlying datum surface is a deposition transition surface, or weathered ablation interface, see fig. 2.
In this embodiment, the two-dimensional or three-dimensional seismic data may be used to track the depth of the corresponding stratum, so that the depths of the overlying reference plane, the underlying reference plane, and the ablation interface may be obtained respectively; further, referring to FIG. 2, the thickness of the impression layer above the ablation interface is obtained by the difference between the formation depth of the ablation interface and the formation depth of the overlying reference plane; the residual formation thickness below the ablation interface is obtained by the difference between the formation depth of the underlying reference surface and the formation depth of the ablation interface.
In a specific application example, the top of the lamp shadow set in the high-abrasion area in the middle of the Sichuan basin is taken as a research object, and under the influence of the motion curtain of the tung-gulf, the two lamp sections and the four lamp sections of the lamp shadow set of the Sichuan basin are all lifted and exposed, wherein the motion curtain of the tung-gulf after the deposition of the four lamp sections has relatively strong effect, and the large-scale degradation exists at the top of the four lamp sections, so that a karst weathering crust reservoir at the top of the lamp shadow set is formed. The zone reservoir control factors include sedimentary facies and karst effects.
According to the deposition evolution process on the top interface of the lamp shadow set, the deposition period of the grazing bamboo temple set-cang lang shop set of the chile system is taken as a compensation deposition period, the degraded paleo-geomorphic form of the lamp shadow set is basically filled up and supplemented, the structural movement of the region in the period is relatively stable, and the structure movement is easy to track in seismic data. The selected overlying datum plane is the top boundary of the lang paving group (the bottom boundary of the temple group of Longwang).
According to the layer sequence division and deposition characteristics of the lamp shadow sets, consider that the third lamp section and the fourth lamp section are a complete sea-intrusion sea-exit continuous deposition loop (a secondary loop); in addition, the three sections of the lamp act as a sea-going formation of mudstone deposits, the formation thickness of which is relatively stable and thin in the basin, resulting in an easy tracking of the full area of the three bottom borders of the lamp. And selecting the three-section bottom interface of the lamp as the underlying reference surface.
Further, the depth of the top boundary of the cang wave pavement group is subtracted from the depth of the top boundary of the lamp shadow group to obtain the thickness of the impression stratum; and subtracting the depth of the top boundary of the lamp shadow set from the depth of the bottom boundary of the three sections of the lamps to obtain the residual stratum thickness.
In step S12 of this embodiment, the relative original formation thickness is obtained according to the impression formation thickness and the residual formation thickness, which includes the following steps:
s121: the impression formation thickness and the residual formation thickness are normalized.
S122: based on the normalized impression formation thickness and the residual formation thickness, the relative original formation thickness is calculated according to the following specific calculation formula:
wherein H is 0 Representing the relative original formation thickness; h r Representing the residual formation thickness after the normalization treatment; h i The thickness of the impression layer after normalization treatment is shown.
In this embodiment, the normalization process adopts a Min-Max normalization method, and the specific calculation formula is as follows:
wherein H is nor Is a normalized value; h * To be standardA value of the melting; h min Is the minimum value of the corresponding thickness; h max Is the maximum value of the corresponding thickness.
It should be noted that, when the carbonate terraced land is degraded, the high portions of the original formation have a larger degraded amount, and the corresponding stamp thickness of the high portions is relatively smaller, and the value of the residual thickness of the high portions is "increased" by the formula (1); similarly, the value of "reduced" residual thickness is used at the low portion of the construction using equation (1). Thus, the effect of differential ablation due to the formation on the residual thickness can be recovered using equation (1), with the result obtained by equation (1) being relative to the original formation thickness.
In step S13 of this embodiment, it is determined whether the residual stratum is affected by deposition, including the steps of:
s131: acquiring a seismic profile of a target karst region, and analyzing the seismic profile according to preset rule conditions, wherein the preset rule conditions comprise: the form of the underlying reference surface is uneven, the underlying reference surface is degraded, and the carbonate bench deposition environment is the development of the beach phase and/or the overburden point of the residual stratum.
S132: if a preset rule condition appears in the seismic profile, determining that the residual stratum is affected by deposition.
S133: if the preset rule condition does not appear in the seismic profile, determining that the residual stratum is not affected by the deposition.
If the underlying reference surface is uneven in the seismic section of the target karst area, the underlying reference surface is degraded, the carbonate bench deposition environment is the beach concurrence and/or the overburden point appears in the residual stratum, the residual stratum is affected by the deposition inheritance. If the underlying datum plane in the seismic section of the target karst region is flat and the carbonate bench starts to develop an interface, or the sedimentary inheritance does not affect the original stratum thickness, the residual stratum is not affected by the sedimentary inheritance.
After determining whether the residual stratum is affected by deposition in step S13 of the present embodiment, the method for restoring the paleo-topography of the karst ablation interface of the present embodiment further includes: if not, determining the thickness of the relative original stratum as a deposited paleo-topography value; if yes, selecting a structural reference surface which is used for reflecting the relief form of the landform before the deposition of the residual stratum; and subtracting the depth corresponding to the construction reference plane from the depth corresponding to the underlying reference plane to obtain the thickness of the underlying stratum, and carrying out standardization treatment on the thickness of the underlying stratum.
It should be noted that, referring to fig. 2, the thickness from the structural reference plane to the underlying reference plane can reflect the relief form of the relief before the deposition of the residual stratum, the low value of the thickness of the underlying stratum represents the low position of the relief, and the high value of the thickness of the underlying stratum represents the high position of the relief, after the influence of the later-stage structure is removed, the more even the structural reference plane is, the better the form is; under the condition of ensuring that the relief of the ancient landform can be reflected, the closer the distance between the structural reference surface and the underlying reference surface is, the better the closer the distance is. In this embodiment, the normalization of the underlying formation thickness uses a Min-Max normalization method.
Correcting the thickness of the relative original stratum to obtain a sedimentary paleo-topography value, wherein the related calculation formula is as follows:
P＝H s *H 0 (3)
wherein P represents a deposition paleo-landform value, H s Representing the underburden thickness; h 0 Representing the relative original formation thickness.
It should be noted that, the thickness of the relative original stratum obtained in step S12 only recovers the area differential degradation, and the carbonate bench deposition is affected by the topography before deposition, where the carbonate bench deposition is raised in the deposition process, is easier to continue to develop into a hillside, while the area with low-lying topography and deeper water body is developed to fill up the dense carbonate deposition mainly, so that the deposition phase change caused by the topography before deposition is often difficult to be reflected, and therefore, the paleo-topography feature before the residual stratum deposition needs to be found to correct the relative original stratum thickness.
In step S14 of this embodiment, recovering the deposited paleo-feature plane according to the deposited paleo-feature value, including: and drawing a deposition paleo-landform plane according to the deposition paleo-landform value, wherein a relatively large value indicates a deposition paleo-landform high position, and a relatively small value indicates a deposition paleo-landform low position.
According to the embodiment, the upper and lower recovery datum planes of the target karst region ablation interface are determined, so that the paleo-topography recovery data of the target karst region, namely the impression stratum thickness and the residual stratum thickness, are conveniently obtained, the relative original stratum thickness is obtained according to the impression stratum thickness and the residual stratum thickness of the target karst region, and the influence of differential ablation caused by the structure on the residual thickness can be recovered; when the residual stratum is not affected by the sedimentary inheritance, the relative original stratum thickness is the sedimentary paleo-landform value; when the residual stratum is affected by the sedimentary inheritance, the relative original stratum thickness is corrected to obtain a sedimentary paleo-topography value, and sedimentary phase change caused by the sedimentary inheritance is eliminated, so that the sedimentary paleo-topography of the target karst area is recovered conveniently, a more comprehensive paleo-topography recovery result is obtained, and the method is helpful for guiding the exploration and development of karst oil and gas reservoirs.
It should be understood that the sequence number of each step in the foregoing embodiment does not mean that the execution sequence of each process should be determined by the function and the internal logic, and should not limit the implementation process of the embodiment of the present invention.
The embodiment also provides a deposition paleo-topography recovery system of the karst ablation interface, which corresponds to the deposition paleo-topography recovery method of the karst ablation interface in the above embodiment one by one. Referring to fig. 3, the depositional paleo-physical restoration system of the karst ablation interface includes a data acquisition module 301, a relative calculation module 302, a judgment module 303, and a restoration module 304. The functional modules are described in detail as follows:
the data acquisition module 301 is configured to acquire paleomorphic restoration data of a target karst region, where the paleomorphic restoration data includes a stamp stratum thickness and a residual stratum thickness;
a relative calculation module 302, configured to obtain a relative original formation thickness according to the impression formation thickness and the residual formation thickness;
the judging module 303 is configured to judge whether the residual stratum is affected by deposition, and if yes, correct the thickness of the relative original stratum to obtain a deposition paleomorphic value;
the restoration module 304 is configured to restore the deposition paleo-feature plane according to the deposition paleo-feature value.
According to the method, the device and the system, the thickness of the original stratum is calculated according to the thickness of the impression and the thickness of the residual stratum in the paleomorphic restoration data, and the influence of differential ablation caused by the structure on the residual thickness can be restored; when the residual stratum is not affected by the sedimentary inheritance, the thickness of the relative original stratum is taken as a sedimentary paleo-topography value; when the residual stratum is affected by the sedimentary inheritance, the relative original stratum thickness is corrected to obtain a sedimentary paleo-topography value, and sedimentary phase change caused by the sedimentary inheritance is eliminated, so that the sedimentary paleo-topography of the target karst area is recovered conveniently, a more comprehensive paleo-topography recovery result is obtained, and the method is helpful for guiding the exploration and development of karst oil and gas reservoirs.
Specific limitations regarding the deposited paleo-topography restoration system of the karst ablation interface can be found in the above limitations on the deposited paleo-topography restoration method of the karst ablation interface, and are not described in detail herein. The modules in the deposited paleo-topography restoration system of the karst ablation interface can be implemented in whole or in part by software, hardware, and combinations thereof. The above modules may be embedded in hardware or may be independent of a processor in the computer device, or may be stored in software in a memory in the computer device, so that the processor may call and execute operations corresponding to the above modules.
The present embodiment also provides a computer device, which may be a server. The computer device includes a processor, a memory, a network interface, and a database connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The database of the computer equipment is used for storing paleo-geomorphic restoration data, relative original stratum thickness, deposited paleo-geomorphic values, deposited paleo-geomorphic restoration planes and other information. The network interface of the computer device is used for communicating with an external terminal through a network connection. The computer program is executed by a processor to realize a method for recovering the paleo-deposition of a karst denudation interface, and the processor executes the computer program to realize the following steps: s11: and obtaining paleomorphic restoration data of the target karst region, wherein the paleomorphic restoration data comprises a impression stratum thickness and a residual stratum thickness. S12: based on the impression formation thickness and the residual formation thickness, a relative raw formation thickness is obtained. S13: and judging whether the residual stratum is affected by deposition, and if so, correcting the thickness of the relative original stratum to obtain a deposition paleomorphic value. S14: and recovering the deposition paleo-landform plane according to the deposition paleo-landform value.
It will be apparent to those skilled in the art that, for convenience and brevity of description, only the above-described division of the functional units and modules is illustrated, and in practical application, the above-described functional distribution may be performed by different functional units and modules according to needs, i.e. the internal structure of the apparatus is divided into different functional units or modules to perform all or part of the above-described functions. The foregoing are all preferred embodiments of the present application, and are not intended to limit the scope of the present application in any way, therefore: all equivalent changes in structure, shape and principle of this application should be covered in the protection scope of this application.
Claims (6)
1. A method for recovering the paleomorphic deposit of karst ablation interface is characterized in that: the method comprises the following steps:
obtaining paleomorphic restoration data of a target karst region, wherein the paleomorphic restoration data comprises impression stratum thickness and residual stratum thickness;
acquiring the relative original stratum thickness according to the impression stratum thickness and the residual stratum thickness;
judging whether the residual stratum is affected by deposition, if so, correcting the thickness of the relative original stratum to obtain a deposition paleo-landform value;
restoring a deposition paleo-landform plane according to the deposition paleo-landform value;
wherein, the step of judging whether the residual stratum is affected by deposition is as follows:
acquiring a seismic profile of a target karst region, and analyzing the seismic profile according to preset rule conditions, wherein the preset rule conditions comprise: the form of the underlying reference surface is uneven, the underlying reference surface is degraded, and the carbonate bench deposition environment is the development of a hillside and beach phase and/or the overburden point of the residual stratum;
if the preset rule condition appears in the seismic profile, determining that the residual stratum is affected by deposition;
if the preset rule condition does not appear in the seismic profile, determining that the residual stratum is not affected by deposition;
judging whether the residual stratum is affected by deposition, if so, selecting a structural reference surface, wherein the structural reference surface is used for reflecting the relief form of the residual stratum before deposition; subtracting the depth corresponding to the construction reference plane from the depth corresponding to the underlying reference plane to obtain the underlying stratum thickness, and carrying out standardization treatment on the underlying stratum thickness;
correcting the relative original stratum thickness to obtain a sedimentary paleo-topography value, wherein the related calculation formula is as follows:
P=
wherein P represents a sedimentary paleo-topography value,representing the underburden thickness; />Representing the relative original formation thickness;
and judging whether the residual stratum is affected by deposition, and if not, determining the thickness of the relative original stratum as a deposition paleo-physical value.
2. The method for restoring the paleo-topography deposited on the karst ablation interface of claim 1, wherein the method comprises the steps of: before obtaining the paleomorphology restoration data of the target karst region, the method further comprises:
constructing an upper and lower recovery datum plane of an ablation interface of a target karst region based on stratum division, wherein the upper and lower recovery datum plane comprises an overlying datum plane and an underlying datum plane;
obtaining the thickness of the impression stratum above the ablation interface according to the depth of the overlying reference surface and the depth of the ablation interface;
and obtaining the thickness of the residual stratum below the degraded interface according to the depth of the degraded interface and the depth of the underlying reference surface.
3. The method for restoring the paleo-topography deposited on the karst ablation interface of claim 1, wherein the method comprises the steps of: obtaining a relative original formation thickness from the impression formation thickness and the residual formation thickness, comprising:
performing standardization treatment on the impression stratum thickness and the residual stratum thickness;
based on the normalized impression formation thickness and the residual formation thickness, the relative original formation thickness is calculated according to the following specific calculation formula:
wherein,representing the relative original formation thickness; />Representing the residual formation thickness after the normalization treatment; />The thickness of the impression layer after normalization treatment is shown.
4. A method for restoring a paleo-topography deposited on a karst ablation interface according to claim 1 or 3, wherein: the normalization processing adopts a Min-Max normalization method, and a specific calculation formula is as follows:
wherein,is a normalized value; />Is the value to be normalized; />Is the minimum value of the corresponding thickness; />Is the maximum value of the corresponding thickness.
5. The utility model provides a deposit paleo-topography recovery system of karst denudation interface which characterized in that: the system is applied to the method of claim 1, the system comprising:
the data acquisition module is used for acquiring paleomorphic restoration data of the target karst area, wherein the paleomorphic restoration data comprises impression stratum thickness and residual stratum thickness;
the relative calculation module is used for obtaining the relative original stratum thickness according to the impression stratum thickness and the residual stratum thickness;
the judging module is used for judging whether the residual stratum is affected by deposition, and if so, correcting the thickness of the relative original stratum to obtain a deposition paleomorphic value;
and the recovery module is used for recovering the deposition paleo-landform plane according to the deposition paleo-landform value.
6. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor, when executing the computer program, implements the steps of the method for restoring the paleo-physical form of a karst ablation interface according to any one of claims 1 to 4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110851311.9A CN113640870B (en) | 2021-07-27 | 2021-07-27 | Method, system and equipment for recovering sedimentary paleo-topography of karst ablation interface |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110851311.9A CN113640870B (en) | 2021-07-27 | 2021-07-27 | Method, system and equipment for recovering sedimentary paleo-topography of karst ablation interface |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113640870A CN113640870A (en) | 2021-11-12 |
| CN113640870B true CN113640870B (en) | 2024-02-23 |
Family
ID=78418560
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110851311.9A Active CN113640870B (en) | 2021-07-27 | 2021-07-27 | Method, system and equipment for recovering sedimentary paleo-topography of karst ablation interface |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113640870B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105740574A (en) * | 2016-03-03 | 2016-07-06 | 中国石油大学(北京) | Palaeo-geomorphic cross section recovery method |
| CN107589470A (en) * | 2017-08-24 | 2018-01-16 | 西南石油大学 | Palaeokarst Landform restoration methods in small area |
| CN110045431A (en) * | 2019-04-19 | 2019-07-23 | 成都理工大学 | A kind of new marine deposit stratum ancient landform bearing calibration based on gamma data body |
| CN111624651A (en) * | 2019-02-28 | 2020-09-04 | 中国石油天然气股份有限公司 | Reservoir prediction method and device based on ancient landform constraint |
| CN112083483A (en) * | 2020-09-17 | 2020-12-15 | 北京中恒利华石油技术研究所 | Sedimentary ancient landform restoration method based on rock physical model de-compaction correction |
| CN112327355A (en) * | 2020-10-28 | 2021-02-05 | 中海石油(中国)有限公司深圳分公司 | Ancient landform restoration method based on source-sink system |
-
2021
- 2021-07-27 CN CN202110851311.9A patent/CN113640870B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105740574A (en) * | 2016-03-03 | 2016-07-06 | 中国石油大学(北京) | Palaeo-geomorphic cross section recovery method |
| CN107589470A (en) * | 2017-08-24 | 2018-01-16 | 西南石油大学 | Palaeokarst Landform restoration methods in small area |
| CN111624651A (en) * | 2019-02-28 | 2020-09-04 | 中国石油天然气股份有限公司 | Reservoir prediction method and device based on ancient landform constraint |
| CN110045431A (en) * | 2019-04-19 | 2019-07-23 | 成都理工大学 | A kind of new marine deposit stratum ancient landform bearing calibration based on gamma data body |
| CN112083483A (en) * | 2020-09-17 | 2020-12-15 | 北京中恒利华石油技术研究所 | Sedimentary ancient landform restoration method based on rock physical model de-compaction correction |
| CN112327355A (en) * | 2020-10-28 | 2021-02-05 | 中海石油(中国)有限公司深圳分公司 | Ancient landform restoration method based on source-sink system |
Non-Patent Citations (1)
| Title |
|---|
| 朱正平 等.川中高石梯—磨溪地区灯四段古地貌恢复及其对储层的控制作用.《中国石油勘探》.2019,第24卷(第6期),第731-736页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113640870A (en) | 2021-11-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111766630B (en) | Restoration method for ancient landform of basin | |
| Lin et al. | Drainage density, slope angle, and relative basin position in Japanese bare lands from high-resolution DEMs | |
| CN105118090A (en) | Adaptive point-cloud filtering method for complex terrain structure | |
| CN110189409B (en) | PLAXIS-based rapid true three-dimensional geological modeling method and system | |
| CN110208862B (en) | Seismic inversion method based on mixed high-order fractional-order ATpV sparse regularization | |
| CN104656131A (en) | Method for improving seismic horizon automatic tracking precision based on waveform recovery technique | |
| Jordan | Adaptive smoothing of valleys in DEMs using TIN interpolation from ridgeline elevations: An application to morphotectonic aspect analysis | |
| CN101877063A (en) | Sub-pixel characteristic point detection-based image matching method | |
| CN115861247A (en) | High-resolution remote sensing image contour multistage regularization method, system and application | |
| CN113640870B (en) | Method, system and equipment for recovering sedimentary paleo-topography of karst ablation interface | |
| CN111476723A (en) | Method for recovering lost pixels of remote sensing image with failed L andsat-7 scanning line corrector | |
| CN114578431A (en) | Lithologic restoration method for ancient landform | |
| CN112083483A (en) | Sedimentary ancient landform restoration method based on rock physical model de-compaction correction | |
| CN113487733A (en) | Method for quickly establishing three-dimensional model of bottom surface of covering layer | |
| CN111948716A (en) | Method for calculating depth of submarine water channel by using seismic data | |
| CN110441814B (en) | Seismic profile iteration layer flattening method for deep water basin | |
| CN116299669A (en) | Recovery method and system for micro paleo-landform evolution process of well-free delta deposition area | |
| CN112904423B (en) | Engraving method and device for land-phase river channel | |
| CN114076982B (en) | Karst paleo-landform restoration method and device based on waveform characteristic difference | |
| CN113792356B (en) | Internal-discharge landform dynamic near-natural remodeling method based on landform feature preservation | |
| CN112965142B (en) | Method, device, equipment and storage medium for generating karst ancient political view | |
| CN108875615B (en) | Deposition area remote sensing identification method and device, electronic equipment and storage medium | |
| CN112904413B (en) | Method, device and storage medium for drawing karst ancient plot | |
| CN112882096B (en) | Method and device for predicting salt window of salt-containing basin and computer storage medium | |
| CN116677374A (en) | Sand body subdivision comparison method and system based on gradient |
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 |